• <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高清观看| 国产精品久久久久久精品电影| 无遮挡黄片免费观看| 久久精品人妻少妇| 三级国产精品欧美在线观看| 亚洲欧美日韩东京热| 欧美日本视频| 成人欧美大片| 老司机深夜福利视频在线观看| 亚洲人与动物交配视频| 男插女下体视频免费在线播放| 国产黄色小视频在线观看| 亚洲精品国产成人久久av| 亚洲一级一片aⅴ在线观看| 精品久久久久久,| 午夜福利在线在线| 身体一侧抽搐| 成年版毛片免费区| 国产午夜精品论理片| 国产在线男女| 日韩一区二区视频免费看| 国产精品一区二区免费欧美| 午夜久久久久精精品| 亚洲av电影不卡..在线观看| 欧美中文日本在线观看视频| 国产成人av教育| aaaaa片日本免费| 非洲黑人性xxxx精品又粗又长| 一本一本综合久久| 久久久国产成人免费| 精品久久久久久久末码| 免费搜索国产男女视频| 欧美性猛交╳xxx乱大交人| 欧美精品国产亚洲| 国产精品一及| 一a级毛片在线观看| 大型黄色视频在线免费观看| 成人毛片a级毛片在线播放| 成人毛片a级毛片在线播放| 亚洲 国产 在线| 午夜福利高清视频| 中文字幕av成人在线电影| 麻豆久久精品国产亚洲av| 国产黄片美女视频| 欧美精品啪啪一区二区三区| 成人欧美大片| 国产亚洲精品av在线| 免费人成在线观看视频色| 亚洲国产精品久久男人天堂| 免费在线观看成人毛片| 干丝袜人妻中文字幕| 中出人妻视频一区二区| 成人鲁丝片一二三区免费| 美女 人体艺术 gogo| 成人性生交大片免费视频hd| 人妻少妇偷人精品九色| 一进一出抽搐动态| 精品人妻熟女av久视频| 亚洲人成伊人成综合网2020| 国产中年淑女户外野战色| 成年人黄色毛片网站| 一本久久中文字幕| 日韩高清综合在线| 免费人成视频x8x8入口观看| 国国产精品蜜臀av免费| 国产成人福利小说| 噜噜噜噜噜久久久久久91| 女人十人毛片免费观看3o分钟| 99国产精品一区二区蜜桃av| 桃红色精品国产亚洲av| 午夜老司机福利剧场| a在线观看视频网站| 日本一二三区视频观看| 国产人妻一区二区三区在| 久久久久性生活片| 亚洲一区高清亚洲精品| 如何舔出高潮| 久久精品国产鲁丝片午夜精品 | 美女被艹到高潮喷水动态| 啦啦啦观看免费观看视频高清| 中文字幕久久专区| 不卡视频在线观看欧美| 99久久九九国产精品国产免费| 免费av不卡在线播放| 日本-黄色视频高清免费观看| 欧美一区二区精品小视频在线| 成人国产一区最新在线观看| 两个人的视频大全免费| av在线观看视频网站免费| 国产 一区精品| 久久国产乱子免费精品| 熟女电影av网| 午夜福利在线观看免费完整高清在 | 欧美成人免费av一区二区三区| 国产亚洲欧美98| 亚洲欧美精品综合久久99| 内地一区二区视频在线| 久久久久久久久久久丰满 | 欧美最新免费一区二区三区| 日韩在线高清观看一区二区三区 | 两个人的视频大全免费| av在线观看视频网站免费| 3wmmmm亚洲av在线观看| 久久中文看片网| 久久久久久国产a免费观看| 少妇人妻一区二区三区视频| 黄片wwwwww| 久久99热这里只有精品18| 午夜a级毛片| 色吧在线观看| 国产成年人精品一区二区| 久久久久久久久大av| 高清日韩中文字幕在线| 亚洲一级一片aⅴ在线观看| АⅤ资源中文在线天堂| 国产黄a三级三级三级人| 国产一区二区三区在线臀色熟女| 成人特级av手机在线观看| 亚洲熟妇熟女久久| 欧美性猛交╳xxx乱大交人| 国产黄色小视频在线观看| 女人十人毛片免费观看3o分钟| 亚洲欧美日韩卡通动漫| 久久久久久久精品吃奶| 久久精品综合一区二区三区| 欧美激情在线99| 久久久色成人| 无人区码免费观看不卡| 一级黄片播放器| 免费看日本二区| 国产av一区在线观看免费| 国产在视频线在精品| 日韩欧美国产在线观看| 成人美女网站在线观看视频| 成人综合一区亚洲| 亚洲最大成人中文| 成年版毛片免费区| 成年女人看的毛片在线观看| 国产乱人视频| ponron亚洲| 久久欧美精品欧美久久欧美| av国产免费在线观看| 97超视频在线观看视频| 五月玫瑰六月丁香| 国产一区二区亚洲精品在线观看| 欧美性猛交╳xxx乱大交人| 91在线精品国自产拍蜜月| 俄罗斯特黄特色一大片| 国产高清有码在线观看视频| 在现免费观看毛片| 久久久久性生活片| 99久久成人亚洲精品观看| 亚洲中文日韩欧美视频| 国产精品久久久久久久电影| 亚洲真实伦在线观看| 久久久久久伊人网av| 国产一区二区激情短视频| 嫩草影院入口| 好男人在线观看高清免费视频| av女优亚洲男人天堂| 亚洲人成伊人成综合网2020| 蜜桃久久精品国产亚洲av| 久久久久久久亚洲中文字幕| 欧美激情国产日韩精品一区| 精品一区二区三区人妻视频| 午夜影院日韩av| 亚洲熟妇熟女久久| 中国美女看黄片| 午夜免费激情av| 99热这里只有是精品在线观看| 九色国产91popny在线| 看片在线看免费视频| 日韩精品中文字幕看吧| 国产三级中文精品| 草草在线视频免费看| 男女边吃奶边做爰视频| 久久久久久久午夜电影| 日本 欧美在线| 亚洲一区高清亚洲精品| 黄色欧美视频在线观看| 亚洲经典国产精华液单| 色在线成人网| 亚州av有码| 亚洲五月天丁香| 乱系列少妇在线播放| 老司机福利观看| 成人一区二区视频在线观看| 久久久久久久久久成人| 中文字幕高清在线视频| 欧美高清性xxxxhd video| 一本精品99久久精品77| 国产老妇女一区| 伊人久久精品亚洲午夜| 国产高清三级在线| 黄色欧美视频在线观看| 国产精品国产高清国产av| 又粗又爽又猛毛片免费看| 国产主播在线观看一区二区| 观看美女的网站| 啦啦啦啦在线视频资源| 中文字幕高清在线视频| 日本黄大片高清| 男人舔女人下体高潮全视频| 久久久久久久午夜电影| 最近在线观看免费完整版| 亚洲图色成人| 特级一级黄色大片| 人人妻,人人澡人人爽秒播| 精品午夜福利在线看| 国产精品久久久久久av不卡| 在线免费观看不下载黄p国产 | 国产三级中文精品| 天堂动漫精品| 欧美xxxx黑人xx丫x性爽| 精品午夜福利视频在线观看一区| 看片在线看免费视频| 在线a可以看的网站| 午夜a级毛片| 嫩草影院新地址| 国产日本99.免费观看| 直男gayav资源| 亚洲欧美日韩东京热| 日本与韩国留学比较| 免费高清视频大片| 久久亚洲真实| 国产中年淑女户外野战色| 国产又黄又爽又无遮挡在线| av福利片在线观看| 3wmmmm亚洲av在线观看| 天美传媒精品一区二区| 欧美丝袜亚洲另类 | 五月伊人婷婷丁香| 少妇的逼好多水| 午夜爱爱视频在线播放| 男女那种视频在线观看| 在线a可以看的网站| 一a级毛片在线观看| 欧美性感艳星| 51国产日韩欧美| 大又大粗又爽又黄少妇毛片口| 在线观看66精品国产| 97人妻精品一区二区三区麻豆| 欧美一区二区亚洲| 在线免费观看不下载黄p国产 | 美女黄网站色视频| 两个人视频免费观看高清| 三级毛片av免费| 99热精品在线国产| 白带黄色成豆腐渣| 无遮挡黄片免费观看| 国产精品女同一区二区软件 | 麻豆成人av在线观看| 草草在线视频免费看| 一区二区三区四区激情视频 | 色噜噜av男人的天堂激情| 我要搜黄色片| 狂野欧美白嫩少妇大欣赏| 天堂影院成人在线观看| 天天一区二区日本电影三级| 人人妻,人人澡人人爽秒播| 五月玫瑰六月丁香| 高清毛片免费观看视频网站| 日本三级黄在线观看| 欧美日韩国产亚洲二区| 久久久久久九九精品二区国产| 国产精品免费一区二区三区在线| 91久久精品国产一区二区成人| 欧美成人免费av一区二区三区| 亚洲av五月六月丁香网| 久久国产乱子免费精品| 成人特级av手机在线观看| 精品欧美国产一区二区三| 中文在线观看免费www的网站| 精品福利观看| 男女下面进入的视频免费午夜| 我要搜黄色片| 免费看av在线观看网站| 国产精品亚洲美女久久久| 日韩欧美国产一区二区入口| 亚洲国产欧美人成| 日本撒尿小便嘘嘘汇集6| 日日夜夜操网爽| 国产一区二区三区视频了| 天堂√8在线中文| 免费观看的影片在线观看| 国内久久婷婷六月综合欲色啪| 日韩国内少妇激情av| 国产精品福利在线免费观看| 欧美性感艳星| 成人二区视频| 特级一级黄色大片| 亚洲精品一区av在线观看| 俄罗斯特黄特色一大片| 桃色一区二区三区在线观看| 中文资源天堂在线| 久久人人精品亚洲av| 天美传媒精品一区二区| 男人狂女人下面高潮的视频| 天天一区二区日本电影三级| 亚洲乱码一区二区免费版| 91久久精品国产一区二区三区| 韩国av在线不卡| 精品国内亚洲2022精品成人| 国产高潮美女av| 精品久久久久久久久av| 日韩亚洲欧美综合| 久久久久久久亚洲中文字幕| 两个人的视频大全免费| 亚洲熟妇熟女久久| 国产真实乱freesex| av黄色大香蕉| 天堂av国产一区二区熟女人妻| 国产高潮美女av| 成人特级av手机在线观看| 美女 人体艺术 gogo| 欧美高清成人免费视频www| 成人毛片a级毛片在线播放| 国产成人福利小说| 久99久视频精品免费| 色综合婷婷激情| 我的老师免费观看完整版| 欧美一区二区国产精品久久精品| 亚洲人与动物交配视频| 赤兔流量卡办理| 午夜福利在线在线| 欧美日韩乱码在线| 中出人妻视频一区二区| 成人国产综合亚洲| 亚洲精品成人久久久久久| 欧美激情在线99| 亚洲国产精品久久男人天堂| 久久精品91蜜桃| 成人特级av手机在线观看| 国产精华一区二区三区| 亚洲不卡免费看| 欧美区成人在线视频| 免费高清视频大片| 久9热在线精品视频| 欧美潮喷喷水| 神马国产精品三级电影在线观看| 免费看美女性在线毛片视频| 桃色一区二区三区在线观看| 日日夜夜操网爽| 中文字幕av成人在线电影| 狂野欧美激情性xxxx在线观看| 国产69精品久久久久777片| 18+在线观看网站| 无遮挡黄片免费观看| 欧美在线一区亚洲| 男女那种视频在线观看| 此物有八面人人有两片| 午夜亚洲福利在线播放| 欧美高清性xxxxhd video| 亚洲第一区二区三区不卡| 国产一区二区亚洲精品在线观看| 韩国av在线不卡| 国内精品一区二区在线观看| 亚州av有码| x7x7x7水蜜桃| 国语自产精品视频在线第100页| 黄色女人牲交| 1024手机看黄色片| 成人综合一区亚洲| 日本与韩国留学比较| 色精品久久人妻99蜜桃| 久久国产精品人妻蜜桃| 此物有八面人人有两片| 69av精品久久久久久| 日韩大尺度精品在线看网址| 人妻夜夜爽99麻豆av| 日韩欧美免费精品| 校园春色视频在线观看| 动漫黄色视频在线观看| a级毛片免费高清观看在线播放| 精品日产1卡2卡| 久久香蕉精品热| 亚洲无线在线观看| 少妇人妻精品综合一区二区 | 色噜噜av男人的天堂激情| a在线观看视频网站| 午夜免费男女啪啪视频观看 | 九九久久精品国产亚洲av麻豆| 男人和女人高潮做爰伦理| av在线亚洲专区| www日本黄色视频网| 成年人黄色毛片网站| 久久久久免费精品人妻一区二区| www日本黄色视频网| 午夜爱爱视频在线播放| 久久精品国产亚洲av涩爱 | 直男gayav资源| 日韩中文字幕欧美一区二区| 91久久精品国产一区二区三区| 91麻豆av在线| 国产极品精品免费视频能看的| 中文字幕熟女人妻在线| 久久国产精品人妻蜜桃| 99国产极品粉嫩在线观看| 精品久久久久久久末码| 成人一区二区视频在线观看| www日本黄色视频网| 日本欧美国产在线视频| 免费av观看视频| 国产高清激情床上av| 日韩 亚洲 欧美在线| 内射极品少妇av片p| 亚洲av.av天堂| 欧美色欧美亚洲另类二区| 一个人看的www免费观看视频| 国产成年人精品一区二区| 在线看三级毛片| 三级男女做爰猛烈吃奶摸视频| 成年女人永久免费观看视频| 欧美中文日本在线观看视频| 日本黄大片高清| 一个人免费在线观看电影| 亚洲av成人精品一区久久| 亚洲美女黄片视频| 成年人黄色毛片网站| 亚洲,欧美,日韩| 97超视频在线观看视频| 成人二区视频| 国产精品一区二区性色av| 女人十人毛片免费观看3o分钟| 国产精品电影一区二区三区| 国产免费av片在线观看野外av| 99精品久久久久人妻精品| 日韩精品中文字幕看吧| 亚洲国产欧洲综合997久久,| 成年女人毛片免费观看观看9| 亚洲色图av天堂| 日本 av在线| 午夜激情福利司机影院| 欧美绝顶高潮抽搐喷水| 亚洲真实伦在线观看| 97热精品久久久久久| 国产精品不卡视频一区二区| 日韩欧美一区二区三区在线观看| 两人在一起打扑克的视频| 免费av不卡在线播放| 日本免费a在线| 欧美最黄视频在线播放免费| 一级a爱片免费观看的视频| 成人av在线播放网站| 在线a可以看的网站| 舔av片在线| 午夜精品一区二区三区免费看| 日本色播在线视频| 99在线视频只有这里精品首页| 伊人久久精品亚洲午夜| 变态另类丝袜制服| 精品99又大又爽又粗少妇毛片 | 欧美在线一区亚洲| 精品久久久久久成人av| 男人舔奶头视频| 99久久精品国产国产毛片| 看免费成人av毛片| 欧美性猛交黑人性爽| 中文字幕人妻熟人妻熟丝袜美| 桃色一区二区三区在线观看| 精品国内亚洲2022精品成人| 色吧在线观看| 人妻少妇偷人精品九色| 国产精品久久久久久久电影| 一个人免费在线观看电影| a级一级毛片免费在线观看| 国产精品人妻久久久影院| 如何舔出高潮| av在线观看视频网站免费| 在线免费十八禁| 国产高清不卡午夜福利| 欧美成人性av电影在线观看| 国产欧美日韩精品一区二区| 高清日韩中文字幕在线| 成人亚洲精品av一区二区| 国产三级中文精品| 真人一进一出gif抽搐免费| 变态另类成人亚洲欧美熟女| 岛国在线免费视频观看| 日日摸夜夜添夜夜添av毛片 | 国产精品伦人一区二区| 欧美+亚洲+日韩+国产| 日韩亚洲欧美综合| 小说图片视频综合网站| 欧美成人性av电影在线观看| 日日摸夜夜添夜夜添av毛片 | 国产成人a区在线观看| 热99在线观看视频| 色av中文字幕| 午夜福利18| 乱系列少妇在线播放| 国产午夜福利久久久久久| 成人无遮挡网站| 欧美激情国产日韩精品一区| 久久久久久国产a免费观看| 亚洲人成网站高清观看| 国产蜜桃级精品一区二区三区| 色av中文字幕| 美女免费视频网站| 麻豆av噜噜一区二区三区| 中文字幕av成人在线电影| 99在线视频只有这里精品首页| 中文亚洲av片在线观看爽| 男女视频在线观看网站免费| 欧美绝顶高潮抽搐喷水| 国国产精品蜜臀av免费| 黄色欧美视频在线观看| 亚洲av免费在线观看| 久久精品国产亚洲av天美| 一a级毛片在线观看| 一级av片app| ponron亚洲| 麻豆国产av国片精品| av.在线天堂| 一区二区三区免费毛片| 看十八女毛片水多多多| 久久6这里有精品| 特级一级黄色大片| 欧美日韩中文字幕国产精品一区二区三区| 欧美一级a爱片免费观看看| 亚洲综合色惰| 亚洲美女搞黄在线观看 | 亚洲成人中文字幕在线播放| 午夜激情福利司机影院| 日韩人妻高清精品专区| 免费搜索国产男女视频| 亚洲va日本ⅴa欧美va伊人久久| 有码 亚洲区| 日日摸夜夜添夜夜添小说| 少妇的逼好多水| 日本一本二区三区精品| 日韩欧美三级三区| 国产大屁股一区二区在线视频| 久久99热这里只有精品18| 亚洲第一电影网av| 在线观看一区二区三区| 亚洲精品456在线播放app | 亚洲成人久久爱视频| 欧美日韩黄片免| 日本一二三区视频观看| 在线a可以看的网站| videossex国产| 婷婷精品国产亚洲av| 亚洲五月天丁香| 三级男女做爰猛烈吃奶摸视频| 欧美成人a在线观看| av福利片在线观看| 精品久久久久久久末码| 精品福利观看| 一个人观看的视频www高清免费观看| 黄色配什么色好看| 国产精品久久久久久av不卡| 国产精品98久久久久久宅男小说| 亚洲黑人精品在线| 亚洲av二区三区四区| 亚洲黑人精品在线| 黄色日韩在线| 一个人观看的视频www高清免费观看| 亚洲av五月六月丁香网| 一个人观看的视频www高清免费观看| 免费无遮挡裸体视频| 99久久无色码亚洲精品果冻| 免费无遮挡裸体视频| 亚洲五月天丁香| 麻豆成人av在线观看| 亚洲精品影视一区二区三区av| 中文亚洲av片在线观看爽| 国产亚洲精品久久久com| 国产一区二区三区视频了| 亚洲四区av| 99久久久亚洲精品蜜臀av| 国产精品一区www在线观看 | 美女cb高潮喷水在线观看| 永久网站在线| 可以在线观看的亚洲视频| 男女啪啪激烈高潮av片| 中文字幕精品亚洲无线码一区| 嫩草影视91久久| 免费在线观看成人毛片| 午夜福利18| 老司机福利观看| 草草在线视频免费看| 在线播放国产精品三级| 免费av毛片视频| 日本一本二区三区精品| 变态另类丝袜制服| 少妇的逼水好多| 亚洲狠狠婷婷综合久久图片| 国产午夜精品论理片| av视频在线观看入口| 三级男女做爰猛烈吃奶摸视频| 禁无遮挡网站| 大型黄色视频在线免费观看| 99热这里只有是精品在线观看| 日本与韩国留学比较| 亚洲av不卡在线观看| 国产蜜桃级精品一区二区三区| 日韩欧美国产在线观看| 国产精品免费一区二区三区在线| 免费在线观看日本一区| 亚洲精品国产成人久久av| 伦精品一区二区三区| 日韩av在线大香蕉| 小说图片视频综合网站| 日韩亚洲欧美综合| 亚洲精品亚洲一区二区| 女人被狂操c到高潮| 久久亚洲精品不卡| 又黄又爽又刺激的免费视频.| 国产一区二区在线av高清观看|