• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    順層鉆孔預抽煤層瓦斯精準布孔模式及工程實踐

    2022-04-02 01:41:32林海飛季鵬飛孔祥國李樹剛竇桂東
    煤炭學報 2022年3期
    關鍵詞:布孔順層煤體

    林海飛,季鵬飛,孔祥國,李樹剛,竇桂東,李 可

    (1.西安科技大學 安全科學與工程學院,陜西 西安 710054;2.西部礦井開采及災害防治教育部重點實驗室,陜西 西安 710054;3.煤炭行業(yè)西部礦井瓦斯智能抽采工程研究中心,陜西 西安 710054;4.陜西彬長小莊礦業(yè)有限公司,陜西 咸陽 710054)

    我國大部分礦區(qū)煤層瓦斯含量較高,瓦斯災害防控壓力大,制約了礦井安全生產。鉆孔預抽煤層瓦斯技術是降低煤層瓦斯含量,預防瓦斯事故發(fā)生,實現煤炭安全開采的根本措施。其中,順層鉆孔預抽煤層瓦斯技術在礦井瓦斯治理方面應用效果顯著。針對煤層瓦斯賦存的區(qū)域性特點,關鍵是精準確定順層鉆孔布置參數,實現煤層瓦斯定時定點預抽達標,對防治礦井瓦斯災害具有重要意義。

    順層鉆孔預抽煤層瓦斯時的影響因素眾多,目前多集中于研究煤層瓦斯壓力、抽采負壓等單因素對其影響,而實際上各因素間存在交互作用。鑒于此,筆者針對彬長礦區(qū)煤層瓦斯賦存特點,采用響應面法設計順層鉆孔參數優(yōu)化方案,建立含瓦斯煤體流固耦合理論模型,利用COMSOL軟件建立三維鉆孔預抽模型,研究不同煤層瓦斯賦存條件下布孔參數單因素及交互因素對瓦斯抽采效果的影響規(guī)律,并對布孔參數進行優(yōu)化設計,以期實現煤層瓦斯精準預抽,并形成一套靈活、簡單、多選擇并且適用于彬長礦區(qū)煤層瓦斯賦存特點的順層預抽鉆孔參數選擇方法。

    1 含瓦斯煤體流固耦合抽采模型

    1.1 基本假設

    根據煤層瓦斯賦存特點,提出如下假設:① 煤體均質且各向同性;② 煤體骨架變形、孔隙壓力及吸附解吸效應引起的變形遠小于煤體自身尺度;③ 視煤體為理想線彈性體,應力應變滿足廣義胡克定律;④ 瓦斯擴散服從Fick定律,滲流服從Darcy定律;⑤ 預抽過程忽略煤體中水分、灰分和賦存溫度變化的影響。

    1.2 含瓦斯煤體滲透率動態(tài)演化方程

    滲透率是研究煤體瓦斯運移規(guī)律的關鍵參數之一。瓦斯抽采過程中,煤層原始應力狀態(tài)被破壞,引起煤體變形,導致孔隙率發(fā)生改變。而滲透率與孔隙率相關,任意時刻煤體滲透率可通過Kozeny-Carman方程與孔隙率建立聯系。另外,外部應力、孔隙壓力和吸附解吸效應均對煤體孔隙產生影響,考慮以上3個因素建立了含瓦斯煤體滲透率動態(tài)演化模型:

    (1)

    式中,為任意時刻煤體滲透率,m;為煤體初始滲透率,m;為任意時刻煤體孔隙率,%;為煤體初始孔隙率,%;為有效應力系數;0為煤體初始體應變;為初始瓦斯壓力,MPa;為煤基質體積模量,MPa;為煤體初始吸附應變;為任意時刻煤體的體應變;為任意時刻瓦斯平衡壓力,MPa;為任意時刻煤體吸附應變。

    1.3 煤體瓦斯擴散滲流控制方程

    原始狀態(tài)煤層瓦斯處于吸附-解吸動態(tài)平衡狀態(tài),吸附態(tài)瓦斯在任意平衡壓力時的含量,可由廣義Langmuir方程表述為

    (2)

    式中,為平衡壓力時的吸附態(tài)瓦斯含量,kg/m;為標況下的瓦斯密度,kg/m;為煤體密度,kg/m;,為Langmuir體積/壓力常數,m/kg,MPa。

    鉆孔預抽煤層瓦斯時原始平衡狀態(tài)被打破,吸附態(tài)瓦斯解吸到煤體微孔隙并在濃度梯度作用下以Fick擴散為主的方式運移到裂隙中;而煤體裂隙中游離態(tài)瓦斯在壓力梯度作用下以Darcy滲流的方式在裂隙中運移。綜合考慮瓦斯擴散與滲流過程,構建了煤體瓦斯運移方程:

    (3)

    式中,為任意時刻吸附態(tài)瓦斯含量,kg/m;為瓦斯擴散系數,m/s;為解吸擴散系數,反映吸附態(tài)瓦斯解吸并向裂隙擴散的難易程度,s;為瓦斯摩爾質量,g/mol;為摩爾氣體常數,J/(mol·K);為煤體溫度,K;為瓦斯動力黏性系數,Pa·s。

    1.4 含瓦斯煤體變形控制方程

    外力作用下煤體產生應力應變,瓦斯吸附解吸效應引起煤體膨脹/收縮變形,而孔隙壓力作用于孔隙壁同樣使煤體變形。根據彈性力學理論,建立考慮有效應力的平衡方程、幾何方程和應力應變方程,得到含瓦斯煤體應力-應變關系式為

    (4)

    式中,為剪切模量,MPa;,,,為,方向位移分量;,為張量標記符號,分別取值1,2,3;為泊松比;為方向瓦斯壓力,MPa;為煤體的體積模量,MPa;為方向吸附解吸應變;為方向凈體力分量。

    通過耦合瓦斯擴散方程、滲流方程和含瓦斯煤體變形方程(式(1)~(4)),再輔以切合實際的初始/邊界條件,可得到考慮煤體滲透率動態(tài)演化的瓦斯抽采流固耦合理論模型。

    2 模擬方案及模擬結果

    2.1 模擬方案

    ..幾何模型與定解條件

    陜西彬長礦區(qū)位于黃隴煤田中部,主采4號煤,煤厚0.15~43.87 m,平均厚度10.64~19.39 m,埋深477.80~995.42 m,平均埋深705.00 m。本次研究共收集大佛寺等8個礦井瓦斯參數。其中,埋深瓦斯含量為1.19~8.32 m/t,透氣性系數為0.04~5.34 m/(MPa·d);瓦斯Langmuir體積、壓力常數依次為19.84~23.06 m/t,0.54~0.95 MPa。

    根據彬長礦區(qū)煤層瓦斯實際賦存特征、開采條件和地應力狀態(tài),建立196.0 m×50.0 m×14.0 m的瓦斯預抽三維模型,如圖1所示。其中,煤層走向長50.0 m,傾向長196.0 m,煤厚14.0 m;煤層走向方向設置50.0 m×3.5 m的巷道煤壁暴露面,巷道氣壓設為標準大氣壓。模擬區(qū)域上覆巖層平均容重為24.5 kN/m,依據垂直地應力等于上覆巖層平均容重與煤層埋深之積,計算得煤層垂直地應力為17.28 MPa,即模型頂部施加均布載荷17.28 MPa,模擬上覆巖層應力;模型底部為固定邊界,四周為滑動邊界,巷道煤壁暴露面為自由邊界。煤層布置平行雙順層鉆孔(模擬最小抽采鉆孔單元瓦斯預抽效果),孔深180.0 m,封孔長度15.5 m,封孔區(qū)域設為內壁不流通邊界條件。模型采用自由四面體進行非結構化網格劃分,并在鉆孔、巷道煤壁處做加密處理,共劃分226 124~297 387個網格,網格質量0.517 3~0.535 1。瓦斯預抽模型涉及的物性參數由現場實測、實驗室測試所得,具體參數見表1。

    圖1 煤層瓦斯預抽三維模型Fig.1 3D model of gas pre-drainage in coal seam

    表1 煤層物理特性參數

    煤層瓦斯含量、透氣性系數、抽采負壓、鉆孔直徑和布孔間距5個因素的取值按模擬方案依次設置。

    ..響應面法模擬方案

    響應面分析法具有試驗次數少、預測精度高等優(yōu)點,在優(yōu)化設計中不僅可以得到響應目標與設計變量的關系,還可得到設計變量最優(yōu)組合,應用Central Composite響應面分析法,開展方案設計。研究表明,煤層瓦斯含量、透氣性系數、抽采負壓、鉆孔直徑及布孔間距5個因素對鉆孔預抽煤層瓦斯效果影響較顯著,依據前述彬長礦區(qū)瓦斯賦存特征和現場抽采經驗等,開展5因素5水平響應面方案設計,以(0,±1,±)編碼(表2),共50組優(yōu)化方案(表3)。

    表2 設計因素編碼與水平

    ..響應目標

    應用響應面分析法,首先應明確響應目標,目前確定布孔參數的指標主要為有效抽采半徑,該指標主要用于確定布孔間距,對于抽采負壓、鉆孔直徑等抽采工藝參數,目前沒有明確的選擇方法,其參數選擇多依賴于現場抽采經驗。相同預抽時間下,反映了瓦斯抽采效果的好壞,值越大說明抽采效果越好,但該指標無法判別鉆孔間有無抽采盲區(qū)或抽采提前達標,忽略了多鉆孔聯合抽采時的抽采疊加效應,對鉆孔間煤層瓦斯抽采情況判斷不夠精準。

    /不僅考慮了抽采負壓、鉆孔直徑、布孔間距等布孔參數對煤層瓦斯抽采效果的影響,同時考慮了多鉆孔聯合抽采的疊加效應,可判斷鉆孔間有無抽采盲區(qū)或抽采提前達標。在抽采疊加效應作用下,最優(yōu)布孔間距應略大于2,為保證鉆孔間煤層瓦斯預抽達標(預抽結束時鉆孔間最大瓦斯壓力≤預抽達標時煤層殘余瓦斯壓力),將/作為響應目標。/<1,表明達標時長小于預期時長,可提前預抽達標;/=1,表明達標時長等于預期時長,規(guī)定時間內恰好預抽達標;/>1,表明達標時長大于預期時長,若要預抽達標需相應延長預抽時間(圖2)。此外,/還能反映布孔參數組合的可靠性、合理性,/≤1且越接近1,表明布孔方案越合理。筆者對/進行顯著性、交互作用分析,并優(yōu)化確定適用于不同煤層瓦斯賦存特征的布孔參數,以期提高瓦斯抽采效率,實現煤層瓦斯精準預抽目標,詳細流程如圖3所示。

    表3 試驗方案與響應結果

    確定/時,需判定煤層瓦斯抽采達標情況,假設中煤體為各向同性均質體,忽略煤層的層理/節(jié)理對抽采的影響,因此/均在兩鉆孔水平方向進行判定,但受相鄰鉆孔相互擾動產生抽采疊加效應的影響,最大瓦斯壓力不一定出現在兩鉆孔中心線中點這一固定位置,因此沿鉆孔水平方向布置測線,監(jiān)測兩鉆孔水平間距內的最大瓦斯壓力。選用《煤礦安全規(guī)程》瓦斯抽采率指標作為判定依據,其判定條件為:煤層瓦斯預抽率>30%,即煤層殘余瓦斯含量≤70%原始瓦斯含量。而煤層殘余瓦斯壓力與殘余瓦斯含量間的換算關系,可由Langmuir方程表述如下

    圖2 布孔參數優(yōu)化示意Fig.2 Schematic diagram of hole parameters optimization

    圖3 煤層瓦斯精準預抽流程Fig.3 Accurate pre-drainage process of coal seam gas

    (5)

    式中,為殘余瓦斯含量,m/t;為游離態(tài)瓦斯含量,m/t;為吸附態(tài)瓦斯含量,m/t。

    已知煤層瓦斯含量時,可由式(5)反算煤層殘余瓦斯壓力(達標瓦斯壓力),進而確定/。

    根據表2給出的布孔參數模擬方案,當煤層瓦斯含量為0.243,3.000,5.000,7.000,9.757 m/t時,若預抽達標,瓦斯壓力需降至0.010 2,0.138 4,0.248 5,0.377 1,0.594 4 MPa以下,并以此作為響應目標的判定標準。而煤層透氣性系數可通過式(6)與滲透率建立聯系,以便于開展數值模擬。

    =2

    (6)

    式中,為透氣性系數,m/(MPa·d);為單位換算系數,包括氣體壓力的單位由MPa換算成Pa、瓦斯流速的單位由m/s換算成m/d兩部分;為標準狀況下大氣壓力,0.101 3 MPa。

    經計算,當煤層透氣性系數分別為:0.122,1.500,2.500,3.500和4.878 m/(MPa·d)時,對應的滲透率為:0.304×10,3.750×10,6.250×10,8.750×10,12.196×10m。

    2.2 模擬結果

    ..煤層瓦斯預抽效果分析

    以各影響因素中編碼值為0的參數為例,分析煤層瓦斯預抽效果。不同預抽時間下煤層瓦斯壓力分布如圖4所示。由圖4可知,預抽100 d時垂直鉆孔方向的抽采達標區(qū)域以鉆孔為中心呈類圓狀分布,相鄰鉆孔間的抽采疊加效應不明顯;隨著預抽時間延長,抽采達標區(qū)域不再孤立存在,相鄰鉆孔間的抽采達標區(qū)域逐漸連通并形成復合抽采達標區(qū);預抽400 d時抽采疊加效應越明顯,復合抽采達標區(qū)域面積進一步擴大。

    圖4 順層鉆孔預抽煤層瓦斯效果Fig.4 Effect of borehole pre-drainage of coal seam gas

    ..響應模型及其適用性分析

    預抽模擬時長設定為400 d,周期相對較長,受篇幅限制,僅以100,175,250,325和400 d為例,詳細分析單因素/交互因素對煤層瓦斯預抽效果的影響規(guī)律,并開展順層鉆孔預抽煤層瓦斯精準布孔參數的優(yōu)化設計,試驗結果見表3。

    對5組試驗結果進行多元回歸擬合,通過對比各模型顯著性檢測和相關性檢驗的數據,預抽100,175,250,325和400 d時/響應模型均確定為五元二階多項式,響應模型及方差分析見表4,5。方差分析顯示:各模型無顯著影響概率()均小于0.000 1,模型回歸效果極顯著。

    表4 不同預抽時間下Pmax/Pb的響應模型

    表5 Pmax/Pb響應模型的方差分析

    續(xù)表

    表6 響應模型誤差分析

    為進一步驗證模型可靠性,以/的試驗值和模型預測值為橫/縱坐標繪制散點圖,如圖5所示。

    由圖5可知,各模型散點均位于=附近,表明模型擬合度較好,可利用該模型進行因素間交互作用分析和順層鉆孔布置參數的優(yōu)化設計。

    3 煤層瓦斯精準預抽影響因素分析

    3.1 單因素對瓦斯精準預抽的影響

    影響煤層瓦斯精準預抽的因素分為地質因素和工程因素,兩者共同決定煤層瓦斯的預抽效果。5個因素中,為地質因素,~為工程因素。對比表5中5因素的值(檢驗結果)可知,不同預抽時間下,/對各因素敏感性均為:>>>>,對/影響最為顯著,若要消除抽采空白帶同時不過度抽采,合理確定順層鉆孔的布孔間距最重要,該結論恰好表明/很好地考慮到布孔間距的確定問題。

    圖5 響應目標Pmax/Pb試驗值與預測值Fig.5 Experimental and predicted values of response target Pmax/Pb

    為直觀分析單因素對/的影響規(guī)律,固定各因素中編碼值為0的參數,并以各自變量水平值為橫軸,響應值為縱軸,擬合/隨各因素變化的曲線,如圖6所示。

    圖6 單因素對Pmax/Pb的影響Fig.6 Influence of single factor on Pmax/Pb

    由圖6可知,或固定不變時,隨著預抽時間延長,單位時間內/的下降量逐漸減小。為5 m/t時,相鄰預抽時間/下降量依次為:0.134 2,0.083 7,0.060 5,0.047 4;為2.5 m/(MPa·d)時,相鄰預抽時間/下降量依次為:0.127 9,0.088 8,0.060 9,0.046 2。是因為煤層瓦斯含量一定時,預抽初期瓦斯壓力與抽采負壓形成的壓差較大,瓦斯?jié)B流動力充足,減小速度快,相應的/下降量越大,而的減小降低了其與抽采負壓形成的壓差,預抽后期瓦斯?jié)B流動力匱乏,減小速度放緩,相應的/下降量減小。煤層瓦斯含量的高低決定了滲流動力的大小,而透氣性系數的高低決定了瓦斯?jié)B流阻力的大小,其作用原理與煤層瓦斯含量類似,在此不再贅述。預抽時間固定不變時,隨著,增大,/逐漸下降,兩者與/呈負相關,預抽250,325,400 d時抽采達標(/≤1)的煤層瓦斯含量范圍依次為:>6.477 4,>4.881 3,>4.093 6 m/t;煤層透氣性系數范圍依次為:>3.418 9,>2.394 1,>2.050 9 m/(MPa·d),即煤層瓦斯含量越高或透氣性系數越大時抽采達標時間越短。

    不同預抽時間下,隨著,增大,/變化不明顯。固定不變時,隨著預抽時間延長,單位時間內/的下降量逐漸減?。划旑A抽時間固定不變時,/隨增加而增大,但/增長速率有所放緩。是因為布孔間距設置過大時,相鄰鉆孔間抽采疊加效應被削弱,瓦斯?jié)B流動力匱乏,單位時間內僅有少量瓦斯被抽出,減小速度慢,/下降量也越??;隨著預抽時間延長,鉆孔間抽采疊加效應逐漸增強,下降速度有所提升,/下降量逐漸增加。

    3.2 交互因素對瓦斯精準預抽的影響

    ..地質因素交互

    /不僅受單因素影響,同時受因素間交互作用影響。,對/的影響非常顯著且均屬地質因素,如圖7所示。由圖7可知,不同預抽時間下與的響應曲面較陡、響應等值線扭曲程度大且分布較密集,兩因素水平范圍內的/分別位于1.10~1.50,1.00~1.40,0.90~1.40,0.80~1.30,0.80~1.20,說明地質因素對/影響非常顯著,對/數值大小影響較明顯并具有一定的時間效應。另外,不同預抽時間下,/的響應等值線均聚集在右下角,表明即使,處于交互作用下,,對/的影響程度仍然顯著,/隨,增加而減速下降,在與的交互作用下,煤層瓦斯含量越高且煤層透氣性系數越大時,抽采效果越好,抽采達標時間越短。

    圖7 地質因素交互作用對Pmax/Pb的影響規(guī)律Fig.7 Influence of interaction between geological factors on Pmax/Pb

    ..工程因素交互

    圖8表明了工程因素交互對/的影響。由圖8(a)可知,不同預抽時間下,與的響應等值線呈圓形分布,兩因素水平范圍內存在極值點使得/達到最小,相鄰響應等值線在數值上僅差0.005,說明抽采負壓與鉆孔直徑的交互作用不顯著,兩者間的交互作用對煤層瓦斯預抽效果影響較小,并與單因素分析結果相互佐證。

    由圖8(b),(c)可知,不同預抽時間下/均沿方向變化較明顯,說明抽采負壓、鉆孔直徑對/的影響程度小于布孔間距,在與抽采負壓、鉆孔直徑交互作用下,布孔間距對/的影響仍然顯著,這與單因素分析結果一致。工程因素之間的交互作用同樣有一定時間效應,隨著預抽時間延長,/不斷降低。

    ..地質因素與工程因素交互

    單因素分析發(fā)現,,對/的影響程度最低,導致,與其他因素之間的交互作用并不明顯。因此,地質因素與工程因素間的交互作用著重分析與,與,如圖9所示。

    圖8 工程因素交互作用對Pmax/Pb的影響Fig.8 Influence of interaction between engineering factors on Pmax/Pb

    圖9 地質因素與工程因素交互對Pmax/Pb的影響Fig.9 Influence of interaction between geological factors and engineering factors on Pmax/Pb

    當布孔間距設置過大時,相鄰鉆孔間的抽采疊加效應被削弱,因此,相同或條件下,隨著布孔間距增大,/逐漸增大。布孔間距對/的影響同樣有時間效應,隨著預抽時間延長,相同布孔間距條件下,不同煤層瓦斯含量的煤層所能達到的/差距逐漸增大;相同煤層瓦斯含量條件下,擴大布孔間距增大/的現象有所改善,這說明布孔間距設置過大時,應延長預抽時間以實現預抽達標。對比圖9發(fā)現,在布孔間距與地質因素的交互作用下,/的響應等值線均沿布孔間距方向變化最為明顯且分布相對密集,說明布孔間距對/的影響程度高于煤層瓦斯含量、煤層透氣性系數,布孔間距對/影響更顯著,與單因素分析結果相互佐證。

    4 順層鉆孔精準布孔方法及應用

    4.1 鉆孔預抽煤層瓦斯精準布孔方法

    分析地質/工程因素各單因素及交互因素對煤層瓦斯預抽效果的影響發(fā)現,地質因素和預抽時間對/影響均比較大。因此,順層鉆孔布置參數的優(yōu)化問題應根據煤層瓦斯賦存的區(qū)域性特點及工程因素隨時間的變化規(guī)律,開展“分時分區(qū)”式布孔,而不是傳統(tǒng)“以點代面”式布孔。

    為實現煤層瓦斯的精準抽采,提高順層預抽鉆孔的瓦斯抽采效率,提出“分時分區(qū)”式順層鉆孔預抽煤層瓦斯精準布孔方法(圖10)。

    圖10 順層鉆孔預抽煤層瓦斯精準布孔方法Fig.10 Accurate borehole placement methodfor pre-draining coal seam gas by drilling along the seam

    所謂“分時分區(qū)”式順層鉆孔預抽煤層瓦斯精準布孔方法具體包括“分時”和“分區(qū)”兩步,根據采掘進度合理安排預抽時間,實現精準的第1步“分時”,保證預抽按時達標;掌握工作面瓦斯賦存規(guī)律的前提下,對工作面瓦斯地質單元進行等級劃分,形成若干個抽采目標區(qū),并依據各目標區(qū)實際瓦斯賦存特點,優(yōu)選合理抽采工藝參數,做到精準的第2步“分區(qū)”,保證各目標區(qū)預抽達標;兩步結合形成適用于不同煤層瓦斯賦存特征的順層鉆孔預抽煤層瓦斯精準布孔方法。

    4.2 順層鉆孔布孔參數優(yōu)化結果

    利用Design-Expert軟件中的優(yōu)化模塊,在5因素水平范圍內開展布孔參數優(yōu)化(表2),依次確定不同預抽時間(100,175,250,325,400 d)下,適用于彬長礦區(qū)煤層瓦斯賦存條件的最優(yōu)化布孔參數(抽采負壓、鉆孔直徑、布孔間距)和與之對應的/。隨后利用Suffer軟件對不同預抽時間下適用于彬長礦區(qū)煤層瓦斯賦存條件的最優(yōu)化布孔參數(抽采負壓、鉆孔直徑、布孔間距)和與之對應的/依次進行等值線化處理,最終繪制了“分區(qū)分時”的最優(yōu)布孔參數等值線分布,如圖11所示。

    圖11依次代表了預抽100,175,250,325,400 d時適用于彬長礦區(qū)煤層瓦斯賦存條件的最優(yōu)化布孔參數。以圖11(a)為例,進行詳細介紹。圖11(a)橫坐標表示了煤層瓦斯含量,縱坐標表示了煤層透氣性系數,由橫縱坐標圈定的范圍代表了彬長礦區(qū)實際煤層瓦斯賦存情況。另外,圖中黑、藍、紅3種不同顏色的等值線,依次代表了適用于彬長礦區(qū)不同煤層瓦斯賦存條件下的最優(yōu)化抽采負壓、鉆孔直徑、布孔間距。

    已知區(qū)域煤層地質因素和預抽時間的前提下,便可由圖11快捷、準確地確定出適用于不同煤層瓦斯賦存特點的最優(yōu)化布孔參數。主要步驟包括:① 明確預抽時間;② 確定分時分區(qū)式布孔參數選擇圖;③ 在圖11中標注試驗工作面具體位置;④ 根據已標注出的具體位置,確定試驗工作面的最優(yōu)化布孔參數。以煤層瓦斯含量5.0 m/t,煤層透氣性系數3.0 m/(MPa·d),預抽時間100 d為例,通過以上3條信息可準確確定圖11(a)中的一點“W”,根據“W”附近3種不同顏色的等值線可確定適用于該工況條件下的最優(yōu)化布孔參數依次為:布孔間距3.92 m、抽采負壓25.3 kPa、鉆孔直徑85 mm。

    現場實際應用時,首先應明確試驗工作面的計劃預抽時間和煤層瓦斯實際賦存狀況,依據煤層瓦斯實際賦存情況將工作面劃分為若干個抽采目標區(qū)域;根據預抽時間確定圖11中與之對應的布孔參數優(yōu)化圖;根據各分區(qū)煤層瓦斯含量和透氣性系數在布孔參數優(yōu)化圖上找到與各分區(qū)瓦斯賦存特征相對應的“定點”或“定區(qū)域”,進而快捷、準確地確定適用于該分區(qū)的最優(yōu)化布孔參數。

    由圖11可知,不同預抽時間下,隨著地質因素的改變,布孔間距優(yōu)選值的變化規(guī)律性最明顯;當煤層瓦斯含量和透氣性系數均比較小時,預抽100,175,250 d時/>1,說明相鄰鉆孔間存在抽采空白帶,是由于試驗方案中布孔間距的水平下限設定過大,此時若要提高布孔參數選取優(yōu)度,避免出現抽采空白帶,布孔間距應略小于3.50 m。此外,隨著預抽時間延長,布孔參數選取優(yōu)度逐步提高,預抽325,400 d時/均等于1,因此還可適當延長預抽時間來消除抽采空白帶,以保證煤層瓦斯預抽達標。

    圖11 分時分區(qū)式布孔參數選擇Fig.11 Time-sharing partition hole layout parameter selection

    4.3 工程應用

    ..預抽鉆孔設計

    陜西彬長某礦屬于高瓦斯礦井,煤層傾角2°~3°,試驗工作面走向長2 108 m,傾向長196 m,回采高度14 m,煤層原始瓦斯含量為5.31~6.22 m/t,煤層透氣性系數為1.54~2.05 m/(MPa·d),工作面預留360 d開展采前預抽。由圖11可知,不同預抽時間下,試驗工作面布孔參數選取優(yōu)度較好,其中布孔間距對瓦斯預抽效果影響最大且優(yōu)選值變化規(guī)律性最明顯。為此,預抽100,175,250,325和400 d時,在試驗區(qū)域最優(yōu)化布孔參數取值范圍內沿布孔間距變化最顯著方向選取~等4個點,擬合建立優(yōu)選布孔參數與預抽時間函數關系,進而確定預抽360 d時的布孔參數優(yōu)化方案,如圖12所示(圖中,,,,,,,,為擬合常數;為預抽時間,d)。

    由圖12知,預抽時間與優(yōu)選抽采負壓、鉆孔直徑均呈二次函數關系、與優(yōu)化布孔間距符合冪函數關系;預抽360 d,工作面優(yōu)選布孔參數范圍為:抽采負壓23.9~22.8 kPa、鉆孔直徑102.4~114.1 mm、布孔間距5.33~6.03 m,可以發(fā)現抽采負壓、鉆孔直徑的優(yōu)選值差距不大,對/影響較小,因此,準確設定順層預抽鉆孔的布孔間距前提下,兩者可根據施工條件做出適當調整。

    圖12 試驗工作面最優(yōu)化布孔參數Fig.12 Optimal hole layout parameters of test working face

    由于試驗區(qū)域大部分煤層瓦斯含量為6.00 m/t,透氣性系數為1.95 m/(MPa·d),最終順層預抽鉆孔布置參數選擇為:布孔間距6.0 m(與相鄰已預抽工作面布孔間距一致)、抽采負壓24.0 kPa、鉆孔直徑113 mm。另外,試驗工作面采高較大,水平鉆孔無法解決整個工作面的瓦斯預抽問題,但前述布孔參數優(yōu)化方案已給出合理的布孔參數(已優(yōu)化出最小抽采鉆孔單元),因此保證最大孔間距不超6.0 m前提下,可適當調整鉆孔仰角以滿足整個工作面的瓦斯預抽。最終試驗區(qū)域兩巷分別施工高中低3層含仰角鉆孔。其中,回風巷鉆孔水平間距1.50 m,垂直間距0.40 m,仰角依次為3°,5°,7°;運輸巷鉆孔水平間距1.50 m,仰角依次為2°,4°,6°,兩巷鉆孔的孔深180.0 m,封孔長度15.5 m。

    ..預抽效果分析

    圖13 試驗工作面瓦斯預抽效果Fig.13 Effectiveness of pre-drainage at test working face

    試驗區(qū)域兩巷煤層瓦斯預抽效果如圖13所示,由圖13可知,預抽期間回風巷/運輸巷抽采瓦斯體積分數維持在2.0%~7.0%,2.0%~9.0%;兩巷瓦斯抽采混合流量與瓦斯體積分數呈相反變化趨勢。當預抽時間達到190 d時,運輸巷瓦斯抽采混合流量出現徒增,這是由于試驗工作面的部分區(qū)域開展了水力割縫、注液態(tài)CO/注氮驅替煤體瓦斯試驗,一定程度上提高了抽采混合流量,相應降低了抽采瓦斯體積分數。為明確優(yōu)化后的順層鉆孔預抽煤層瓦斯效果,預抽結束后,在僅開展本煤層瓦斯預抽區(qū)域測定了煤層殘余瓦斯含量并計算了煤層瓦斯預抽率。經實測,僅開展煤層鉆孔瓦斯預抽360 d時預抽率為33.50%左右,說明預抽已達標,一定程度上保證了試驗工作面定時定點預抽達標,確保了回采工作的順利進行。

    5 結 論

    (1)建立了考慮煤體滲透率動態(tài)演化、瓦斯擴散滲流及煤體變形等影響的含瓦斯煤體流固耦合抽采模型;多鉆孔聯合抽采時,應考慮相鄰鉆孔間的抽采疊加效應,隨著預抽時間延長,垂直鉆孔方向的抽采達標區(qū)域逐漸連通并形成復合抽采達標區(qū)域。

    (2)鉆孔預抽煤層瓦斯效果與地質因素、工程因素和預抽時間有關;不同預抽時間下,鉆孔間最大瓦斯壓力與達標壓力比(/)對各因素敏感性依次為布孔間距、煤層瓦斯含量、透氣性系數、鉆孔直徑、抽采負壓,其中布孔間距對/影響程度最大。

    (3)不同預抽時間下,地質因素中煤層瓦斯含量和透氣性系數的交互作用對/的影響非常顯著,對/數值大小影響較大;工程因素中抽采負壓與鉆孔直徑的交互作用不明顯,對/數值大小影響較??;布孔間距與地質因素的交互作用下,布孔間距對/的影響仍然是顯著的;各因素之間的交互作用對/影響均存在一定時間效應。

    (4)綜合考慮地質因素的空間分布規(guī)律及工程因素隨時間變化規(guī)律,提出了煤層“分時分區(qū)”式順層鉆孔預抽煤層瓦斯精準布孔方法,可劃分抽采目標區(qū)并精準確定各目標區(qū)布孔參數;試驗工作面利用優(yōu)化布孔參數布置順層預抽鉆孔,在規(guī)定時間內預抽達標,實現了煤層瓦斯精準預抽。

    猜你喜歡
    布孔順層煤體
    彎曲載荷下不同布孔形式篩管極限彎矩分析*
    石油機械(2023年10期)2023-10-17 03:09:52
    V形起爆網路分析及其應用
    采礦技術(2021年2期)2021-04-17 07:07:16
    消除鉆孔瓦斯抽采盲區(qū)空白帶的布孔方式研究
    基于三維數值模擬的含軟弱夾層順層巖質邊坡開挖穩(wěn)定性研究
    礦產勘查(2020年2期)2020-12-28 00:24:34
    注熱井周圍煤體蠕變過程的滲透率變化規(guī)律模擬研究
    預應力錨索在公路順層巖質邊坡中的應用
    工程與建設(2019年4期)2019-10-10 01:45:46
    紅砂巖順層邊坡監(jiān)測及變形破壞探究
    正交試驗下煤體滲透性影響因素評價
    煤炭工程(2019年4期)2019-05-05 02:30:16
    以“中央廚房”為突破口探索時政報道的融煤體之路——以浙江之聲為例
    傳媒評論(2017年7期)2017-11-08 03:11:15
    試論露天深孔臺階爆破布孔工藝
    亚洲欧洲日产国产| 亚洲国产欧美网| 成人亚洲欧美一区二区av| 久久久久久久国产电影| 永久网站在线| 国产成人精品无人区| 一区福利在线观看| 男人舔女人的私密视频| 在线免费观看不下载黄p国产| 国产欧美亚洲国产| 少妇被粗大的猛进出69影院| 国产色婷婷99| 少妇 在线观看| 午夜91福利影院| 午夜福利影视在线免费观看| 国产精品久久久久久av不卡| 亚洲成人手机| 卡戴珊不雅视频在线播放| 制服人妻中文乱码| av天堂久久9| 天堂8中文在线网| 王馨瑶露胸无遮挡在线观看| av电影中文网址| 久久午夜福利片| 国产精品一区二区在线不卡| 搡老乐熟女国产| 国产人伦9x9x在线观看 | 国产成人av激情在线播放| 久久精品久久精品一区二区三区| 久久久国产精品麻豆| 中文字幕亚洲精品专区| www.精华液| 最近最新中文字幕免费大全7| 国产欧美日韩综合在线一区二区| 欧美精品国产亚洲| 日韩中文字幕视频在线看片| 欧美少妇被猛烈插入视频| 精品第一国产精品| 国产激情久久老熟女| 国产97色在线日韩免费| 国产精品一国产av| 麻豆av在线久日| 老鸭窝网址在线观看| 国产av一区二区精品久久| 国产精品麻豆人妻色哟哟久久| 精品第一国产精品| 日韩一卡2卡3卡4卡2021年| a级毛片在线看网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久精品久久精品一区二区三区| 国产男女内射视频| 一区二区av电影网| 国产精品成人在线| 久久人妻熟女aⅴ| 精品国产一区二区三区四区第35| 国产片特级美女逼逼视频| 亚洲av电影在线观看一区二区三区| 国产97色在线日韩免费| 波多野结衣一区麻豆| 男女下面插进去视频免费观看| 亚洲在久久综合| 秋霞伦理黄片| 人妻一区二区av| 国产乱人偷精品视频| 国产精品久久久久久久久免| 妹子高潮喷水视频| 蜜桃国产av成人99| 国产成人欧美| 侵犯人妻中文字幕一二三四区| 色播在线永久视频| 卡戴珊不雅视频在线播放| 少妇猛男粗大的猛烈进出视频| 欧美日韩视频高清一区二区三区二| 亚洲精品国产色婷婷电影| 欧美少妇被猛烈插入视频| 精品一区二区免费观看| 这个男人来自地球电影免费观看 | 女人被躁到高潮嗷嗷叫费观| 国产 精品1| 丝袜脚勾引网站| 下体分泌物呈黄色| 亚洲欧美清纯卡通| 嫩草影院入口| 毛片一级片免费看久久久久| 亚洲成人手机| 人妻人人澡人人爽人人| 视频区图区小说| 日本免费在线观看一区| 亚洲精品美女久久av网站| 国产男女超爽视频在线观看| 久久人人97超碰香蕉20202| 亚洲国产欧美日韩在线播放| 久久 成人 亚洲| 日韩大片免费观看网站| 国产成人免费观看mmmm| 男女免费视频国产| 日韩电影二区| 18禁国产床啪视频网站| 国产精品久久久久久久久免| 丰满饥渴人妻一区二区三| 精品国产国语对白av| 菩萨蛮人人尽说江南好唐韦庄| 青草久久国产| 中文字幕另类日韩欧美亚洲嫩草| 免费黄频网站在线观看国产| 极品少妇高潮喷水抽搐| 侵犯人妻中文字幕一二三四区| 人人澡人人妻人| 热re99久久国产66热| 七月丁香在线播放| 亚洲精品日韩在线中文字幕| 一区二区三区精品91| 美女高潮到喷水免费观看| 熟女av电影| 久久久国产欧美日韩av| 免费日韩欧美在线观看| av国产久精品久网站免费入址| 久久久久人妻精品一区果冻| 黄片无遮挡物在线观看| 老司机影院毛片| 99国产综合亚洲精品| 男女无遮挡免费网站观看| 国产亚洲av片在线观看秒播厂| 国产乱来视频区| 黄片播放在线免费| 少妇的丰满在线观看| 精品国产露脸久久av麻豆| 成年女人在线观看亚洲视频| 成年女人在线观看亚洲视频| 性色av一级| av电影中文网址| 亚洲欧美一区二区三区久久| 在线观看免费高清a一片| 亚洲av综合色区一区| 性色av一级| 丝袜在线中文字幕| 亚洲,欧美精品.| 赤兔流量卡办理| 久久久久久人妻| 十八禁高潮呻吟视频| 国产在线免费精品| 涩涩av久久男人的天堂| 中文字幕精品免费在线观看视频| 免费不卡的大黄色大毛片视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 成人毛片a级毛片在线播放| 亚洲男人天堂网一区| 伦精品一区二区三区| 免费黄网站久久成人精品| 亚洲欧美日韩另类电影网站| 中文乱码字字幕精品一区二区三区| 午夜日韩欧美国产| 高清不卡的av网站| 各种免费的搞黄视频| 免费播放大片免费观看视频在线观看| 丝袜人妻中文字幕| 日韩不卡一区二区三区视频在线| 国产福利在线免费观看视频| 亚洲精品aⅴ在线观看| 宅男免费午夜| 国产亚洲精品第一综合不卡| 老女人水多毛片| a 毛片基地| 精品少妇内射三级| 啦啦啦中文免费视频观看日本| 18禁裸乳无遮挡动漫免费视频| 久久精品亚洲av国产电影网| 久久精品亚洲av国产电影网| 国产熟女欧美一区二区| 日本黄色日本黄色录像| 亚洲成av片中文字幕在线观看 | 国产伦理片在线播放av一区| 波多野结衣av一区二区av| 国产一级毛片在线| 涩涩av久久男人的天堂| 另类精品久久| 男人爽女人下面视频在线观看| 亚洲精品成人av观看孕妇| 99久久综合免费| 精品国产国语对白av| 黄色毛片三级朝国网站| 亚洲精品日本国产第一区| 国产精品蜜桃在线观看| 日本91视频免费播放| 亚洲第一青青草原| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品自拍成人| 最新中文字幕久久久久| 有码 亚洲区| 久久韩国三级中文字幕| 亚洲四区av| 欧美 日韩 精品 国产| 国产又爽黄色视频| 久久精品国产亚洲av涩爱| 亚洲精品美女久久久久99蜜臀 | 2018国产大陆天天弄谢| 欧美少妇被猛烈插入视频| 久久午夜福利片| 国产精品久久久久久精品古装| 狠狠精品人妻久久久久久综合| 国产精品久久久av美女十八| 日韩在线高清观看一区二区三区| 久久久国产精品麻豆| 最黄视频免费看| 制服丝袜香蕉在线| 汤姆久久久久久久影院中文字幕| 亚洲精品乱久久久久久| 一区二区三区精品91| 久久这里只有精品19| 国产精品 国内视频| 女人久久www免费人成看片| 国产熟女午夜一区二区三区| 90打野战视频偷拍视频| 成人亚洲精品一区在线观看| 伦理电影大哥的女人| 亚洲国产精品999| 国产成人精品一,二区| 欧美+日韩+精品| 欧美日韩一区二区视频在线观看视频在线| 亚洲av中文av极速乱| 欧美精品av麻豆av| 丰满少妇做爰视频| 另类精品久久| 18+在线观看网站| 色94色欧美一区二区| 国产福利在线免费观看视频| 少妇人妻 视频| 国产国语露脸激情在线看| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品在线美女| 精品99又大又爽又粗少妇毛片| 韩国高清视频一区二区三区| 自线自在国产av| 亚洲三级黄色毛片| 午夜福利,免费看| 老司机影院毛片| 精品人妻偷拍中文字幕| 日本色播在线视频| 最近最新中文字幕免费大全7| 男女国产视频网站| 91在线精品国自产拍蜜月| 9191精品国产免费久久| 丝袜脚勾引网站| 欧美成人精品欧美一级黄| 女人精品久久久久毛片| 免费av中文字幕在线| 亚洲精品国产av成人精品| 丰满饥渴人妻一区二区三| 在线观看免费视频网站a站| 亚洲经典国产精华液单| 亚洲精品自拍成人| 欧美中文综合在线视频| 激情视频va一区二区三区| 97在线人人人人妻| 午夜免费观看性视频| 伦精品一区二区三区| 国产 精品1| 国产一区有黄有色的免费视频| 亚洲中文av在线| 成年美女黄网站色视频大全免费| 免费观看a级毛片全部| av网站在线播放免费| 啦啦啦在线免费观看视频4| 久久国内精品自在自线图片| 亚洲精品一区蜜桃| 秋霞伦理黄片| 久久女婷五月综合色啪小说| 街头女战士在线观看网站| 男女边摸边吃奶| 日本猛色少妇xxxxx猛交久久| 天美传媒精品一区二区| 超碰成人久久| 另类精品久久| 日韩视频在线欧美| 欧美激情极品国产一区二区三区| 成人毛片60女人毛片免费| 久久精品久久久久久噜噜老黄| 国产精品亚洲av一区麻豆 | 男女啪啪激烈高潮av片| 夫妻午夜视频| 国产成人精品在线电影| 欧美日本中文国产一区发布| 青草久久国产| 日韩成人av中文字幕在线观看| 中国国产av一级| 日韩欧美精品免费久久| 老汉色∧v一级毛片| 欧美97在线视频| 看十八女毛片水多多多| 成人黄色视频免费在线看| 亚洲内射少妇av| 国产免费一区二区三区四区乱码| 一级a爱视频在线免费观看| 亚洲,欧美精品.| 中文欧美无线码| 国产av国产精品国产| 久久久精品区二区三区| 超碰成人久久| 麻豆精品久久久久久蜜桃| 免费少妇av软件| 免费高清在线观看日韩| 超色免费av| 欧美另类一区| 亚洲五月色婷婷综合| 少妇被粗大的猛进出69影院| 美女大奶头黄色视频| 成人国产av品久久久| 99国产综合亚洲精品| 国产精品一国产av| 搡老乐熟女国产| 日韩一区二区三区影片| 秋霞伦理黄片| 男女免费视频国产| 在线观看三级黄色| 下体分泌物呈黄色| 精品一品国产午夜福利视频| 永久免费av网站大全| 一级片'在线观看视频| 色婷婷久久久亚洲欧美| 欧美精品高潮呻吟av久久| 精品人妻偷拍中文字幕| 国产欧美日韩综合在线一区二区| 欧美日韩视频高清一区二区三区二| videossex国产| 久久久久久久亚洲中文字幕| 秋霞伦理黄片| 亚洲人成电影观看| av在线观看视频网站免费| 国产亚洲欧美精品永久| 秋霞在线观看毛片| 女的被弄到高潮叫床怎么办| 国产精品不卡视频一区二区| 亚洲国产日韩一区二区| 一区二区三区精品91| 成人国产av品久久久| 国产精品一二三区在线看| 日韩在线高清观看一区二区三区| 五月开心婷婷网| 免费人妻精品一区二区三区视频| 搡女人真爽免费视频火全软件| 亚洲精品久久午夜乱码| 在线 av 中文字幕| 欧美av亚洲av综合av国产av | 欧美激情极品国产一区二区三区| 国产有黄有色有爽视频| 大片电影免费在线观看免费| 亚洲色图综合在线观看| 国产又爽黄色视频| 69精品国产乱码久久久| 精品久久久精品久久久| 激情视频va一区二区三区| 国产亚洲精品第一综合不卡| 国产毛片在线视频| 精品酒店卫生间| 老鸭窝网址在线观看| 99热网站在线观看| 久久久精品94久久精品| 日韩免费高清中文字幕av| 国产精品久久久久久精品古装| 一级爰片在线观看| 在线天堂中文资源库| 母亲3免费完整高清在线观看 | 有码 亚洲区| 日本欧美国产在线视频| 晚上一个人看的免费电影| 国产伦理片在线播放av一区| 成人毛片60女人毛片免费| 色婷婷av一区二区三区视频| 久久99热这里只频精品6学生| 欧美精品亚洲一区二区| 久久精品国产鲁丝片午夜精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品国产三级专区第一集| 爱豆传媒免费全集在线观看| av片东京热男人的天堂| 久久久久久免费高清国产稀缺| 国产成人精品久久久久久| 成人黄色视频免费在线看| 97在线视频观看| 不卡av一区二区三区| 九九爱精品视频在线观看| 亚洲,欧美,日韩| 视频在线观看一区二区三区| 国产男人的电影天堂91| 一级毛片我不卡| 久久久久久人人人人人| 黄片小视频在线播放| 电影成人av| 欧美少妇被猛烈插入视频| 高清欧美精品videossex| 青春草视频在线免费观看| 国产高清不卡午夜福利| 2018国产大陆天天弄谢| 欧美日韩国产mv在线观看视频| 亚洲欧美色中文字幕在线| 亚洲中文av在线| 看非洲黑人一级黄片| 黑人猛操日本美女一级片| 国产片特级美女逼逼视频| 美女午夜性视频免费| 国产av精品麻豆| 中国三级夫妇交换| 久久精品国产综合久久久| 久久国产精品男人的天堂亚洲| 精品国产一区二区三区久久久樱花| 免费久久久久久久精品成人欧美视频| 热re99久久国产66热| videos熟女内射| 免费高清在线观看日韩| 搡女人真爽免费视频火全软件| 日日爽夜夜爽网站| 欧美日韩亚洲高清精品| 国产精品嫩草影院av在线观看| 亚洲精品国产av蜜桃| av.在线天堂| 久久久国产一区二区| 亚洲四区av| 国产伦理片在线播放av一区| 亚洲一区二区三区欧美精品| xxxhd国产人妻xxx| 妹子高潮喷水视频| 免费大片黄手机在线观看| 日韩中字成人| 国产av国产精品国产| 伊人久久大香线蕉亚洲五| 五月开心婷婷网| 欧美另类一区| 亚洲欧美日韩另类电影网站| 欧美+日韩+精品| 波多野结衣一区麻豆| 免费观看性生交大片5| 天天躁夜夜躁狠狠躁躁| 久久精品人人爽人人爽视色| 欧美另类一区| 女人高潮潮喷娇喘18禁视频| 欧美日韩视频精品一区| 亚洲欧美日韩另类电影网站| 永久免费av网站大全| 在现免费观看毛片| 日韩一区二区视频免费看| 日韩不卡一区二区三区视频在线| 亚洲精品av麻豆狂野| 国产爽快片一区二区三区| 国产精品一区二区在线观看99| 日韩一卡2卡3卡4卡2021年| 色婷婷久久久亚洲欧美| 在线观看免费视频网站a站| 成人国产麻豆网| 99九九在线精品视频| 久久久久久久亚洲中文字幕| 人成视频在线观看免费观看| 国产精品免费大片| 熟妇人妻不卡中文字幕| 国产成人精品久久二区二区91 | 国产一区二区三区综合在线观看| 最黄视频免费看| 高清欧美精品videossex| 亚洲国产欧美在线一区| 国产一区二区 视频在线| 亚洲少妇的诱惑av| 一区二区三区精品91| 一区二区三区乱码不卡18| 少妇被粗大的猛进出69影院| 一本—道久久a久久精品蜜桃钙片| 一级黄片播放器| 91午夜精品亚洲一区二区三区| www日本在线高清视频| av.在线天堂| 午夜福利视频精品| 在线精品无人区一区二区三| 国产视频首页在线观看| 精品一品国产午夜福利视频| 日韩中文字幕视频在线看片| 日日啪夜夜爽| 日韩制服骚丝袜av| 国产精品久久久久久久久免| 亚洲天堂av无毛| 久久久久久免费高清国产稀缺| 1024视频免费在线观看| 国产精品嫩草影院av在线观看| 日韩视频在线欧美| 91久久精品国产一区二区三区| 午夜福利,免费看| 好男人视频免费观看在线| 天天操日日干夜夜撸| 亚洲成人一二三区av| 亚洲av男天堂| 秋霞伦理黄片| 啦啦啦视频在线资源免费观看| 制服丝袜香蕉在线| 超碰97精品在线观看| 女人精品久久久久毛片| 亚洲精华国产精华液的使用体验| 大片免费播放器 马上看| 男女下面插进去视频免费观看| 啦啦啦啦在线视频资源| 色婷婷久久久亚洲欧美| 亚洲美女黄色视频免费看| 新久久久久国产一级毛片| 欧美少妇被猛烈插入视频| 黄色一级大片看看| 亚洲一区中文字幕在线| 欧美日韩精品成人综合77777| 亚洲欧美色中文字幕在线| 男女高潮啪啪啪动态图| 国产在线免费精品| 欧美日韩av久久| 欧美日韩精品成人综合77777| 午夜影院在线不卡| 色视频在线一区二区三区| 欧美激情高清一区二区三区 | 女人高潮潮喷娇喘18禁视频| 91久久精品国产一区二区三区| 女性生殖器流出的白浆| 丝袜脚勾引网站| 久久精品国产自在天天线| 欧美av亚洲av综合av国产av | 妹子高潮喷水视频| 亚洲成国产人片在线观看| 久久av网站| 日韩免费高清中文字幕av| 只有这里有精品99| 老女人水多毛片| 两个人看的免费小视频| 99国产综合亚洲精品| 麻豆精品久久久久久蜜桃| 国产成人免费观看mmmm| 亚洲,欧美精品.| 美国免费a级毛片| av在线观看视频网站免费| 男女午夜视频在线观看| 黑人欧美特级aaaaaa片| 久久人妻熟女aⅴ| 婷婷色av中文字幕| 中文精品一卡2卡3卡4更新| 捣出白浆h1v1| 亚洲欧美中文字幕日韩二区| 免费播放大片免费观看视频在线观看| 午夜福利乱码中文字幕| 毛片一级片免费看久久久久| 在线观看三级黄色| 日韩,欧美,国产一区二区三区| 在线 av 中文字幕| 免费在线观看黄色视频的| 久久精品aⅴ一区二区三区四区 | 久久精品国产亚洲av涩爱| 老司机亚洲免费影院| 国产一区二区三区av在线| 免费少妇av软件| 亚洲,一卡二卡三卡| 欧美精品人与动牲交sv欧美| 韩国精品一区二区三区| 搡老乐熟女国产| 中文字幕另类日韩欧美亚洲嫩草| 超碰97精品在线观看| 中文精品一卡2卡3卡4更新| 九草在线视频观看| 捣出白浆h1v1| 下体分泌物呈黄色| 久久这里只有精品19| 18禁国产床啪视频网站| 超碰成人久久| 免费在线观看黄色视频的| 欧美变态另类bdsm刘玥| 久久午夜福利片| 在线观看免费日韩欧美大片| 在线观看一区二区三区激情| 国产av国产精品国产| 久久久亚洲精品成人影院| 极品少妇高潮喷水抽搐| 欧美日韩视频高清一区二区三区二| 黑人猛操日本美女一级片| 成人亚洲精品一区在线观看| 国产乱来视频区| 国产一区亚洲一区在线观看| 岛国毛片在线播放| 亚洲五月色婷婷综合| 亚洲成国产人片在线观看| 久久鲁丝午夜福利片| 亚洲欧美一区二区三区国产| 亚洲av.av天堂| 亚洲精品美女久久久久99蜜臀 | 夜夜骑夜夜射夜夜干| 亚洲美女黄色视频免费看| 久久狼人影院| 亚洲色图综合在线观看| 日韩一本色道免费dvd| 国产野战对白在线观看| 人人澡人人妻人| 久久人妻熟女aⅴ| 天天躁夜夜躁狠狠躁躁| 久久午夜综合久久蜜桃| 在线看a的网站| 欧美日韩亚洲高清精品| 多毛熟女@视频| 国产熟女欧美一区二区| 欧美日韩亚洲高清精品| 国产亚洲最大av| 永久网站在线| 黑人巨大精品欧美一区二区蜜桃| 黑人欧美特级aaaaaa片| 交换朋友夫妻互换小说| av在线app专区| 91在线精品国自产拍蜜月| 亚洲av国产av综合av卡| 巨乳人妻的诱惑在线观看| 一二三四中文在线观看免费高清| 亚洲人成电影观看| 美女中出高潮动态图| 嫩草影院入口| 国产精品一国产av| 考比视频在线观看| 欧美 亚洲 国产 日韩一| 亚洲av成人精品一二三区| 伦精品一区二区三区|