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

    煤與瓦斯突出相似模擬試驗(yàn)的能量原理

    2023-12-11 10:14:16程遠(yuǎn)平楊斯杰
    煤炭學(xué)報(bào) 2023年11期
    關(guān)鍵詞:型煤煤體游離

    程遠(yuǎn)平 , 雷 楊 , 楊斯杰

    (1.煤礦瓦斯治理國家工程研究中心, 江蘇 徐州 221116;2.貴州大學(xué) 礦業(yè)學(xué)院, 貴州 貴陽 550025;3.煤礦瓦斯與火災(zāi)防治教育部重點(diǎn)實(shí)驗(yàn)室, 江蘇 徐州 221116;4.中國礦業(yè)大學(xué) 安全工程學(xué)院, 江蘇 徐州 221116)

    煤與瓦斯突出是地下煤礦開采中常見的動(dòng)力災(zāi)害,其特點(diǎn)是煤巖體系在極短時(shí)間內(nèi)連續(xù)破裂失穩(wěn),并伴隨高壓瓦斯猛烈的噴向采掘空間[1-3]。隨著我國淺部煤炭資源逐漸枯竭,深部資源開采將逐漸成為新常態(tài)[4],而如何應(yīng)對深部煤層更為頻繁和復(fù)雜的突出災(zāi)害是煤礦安全生產(chǎn)面臨的一個(gè)重大挑戰(zhàn)。由于技術(shù)和空間限制,實(shí)際突出過程難以直接觀測和記錄,因此開展煤與瓦斯突出的模擬試驗(yàn)是研究突出機(jī)理和獲取相關(guān)數(shù)據(jù)的重要手段。

    最早的實(shí)驗(yàn)室突出模擬可以追溯到20 世紀(jì)50年代,由蘇聯(lián)學(xué)者霍多特[5]使用CO2在實(shí)驗(yàn)室完成,其基本原理是將突出過程視為由煤中高壓瓦斯驅(qū)動(dòng)的動(dòng)力過程,該試驗(yàn)方案為后來的研究者設(shè)計(jì)模擬試驗(yàn)提供了重要參考?;谠囼?yàn)觀測數(shù)據(jù),一些關(guān)鍵的突出特征逐漸揭示,有效推動(dòng)了突出防治技術(shù)的進(jìn)步。例如,KRAVCHENKO 等[6]使用示波器記錄了25 次突出實(shí)驗(yàn)中煤樣的破壞波速,強(qiáng)調(diào)了充分抽采煤層瓦斯和避免煤層突然暴露的重要性。YARTSEV 等[7]通過多次突出實(shí)驗(yàn)發(fā)現(xiàn),突出激發(fā)是由于煤體暴露面瓦斯壓力突然降低引起的,進(jìn)一步突顯了降低瓦斯壓力的重要性。HARGRAVES 等[8]利用突出煤層的煤樣與不同氣體在不同裝置內(nèi)進(jìn)行突出實(shí)驗(yàn),認(rèn)為突出與工作面前方的應(yīng)力集中和瓦斯壓力梯度密切相關(guān),強(qiáng)調(diào)了應(yīng)力條件的控制。KUROIWA 等[9]使用壓制成型的煤樣進(jìn)行不同氣體壓力下的突出實(shí)驗(yàn),指出瓦斯壓力越大,煤體粉化程度和瓦斯涌出量越大,將瓦斯能量釋放與煤體破碎聯(lián)系起來。BODZIONY 等[10-13]進(jìn)行了較為系統(tǒng)的突出模擬試驗(yàn),探討了煤粉粒度分布、水分及孔隙率對型煤強(qiáng)度的影響,不同實(shí)驗(yàn)氣體(CO2、N2、He)對突出閾值的影響,以及滲透率對突出傳播速度的影響等。

    總體而言,物理相似模擬是研究突出機(jī)理和揭示突出特征的重要方法,這些早期的實(shí)驗(yàn)研究不僅為突出防治提供了指導(dǎo),也為現(xiàn)今的儀器設(shè)計(jì)和試驗(yàn)?zāi)M方案制定提供了重要參考。

    1 煤與瓦斯突出模擬試驗(yàn)的基本原理

    近年來,我國的學(xué)者們研發(fā)了許多的突出模擬試驗(yàn)裝置。袁亮等[2,14-16]研發(fā)了一系列突出模擬試驗(yàn)系統(tǒng),成功模擬了巷道掘進(jìn)揭煤誘導(dǎo)突出的物理過程。李術(shù)才等[17]研發(fā)的大尺度真三軸突出定量物理模擬試驗(yàn)系統(tǒng)考慮了實(shí)際煤層環(huán)境,實(shí)現(xiàn)了三維氣固耦合條件下的突出模擬。文光才等[18]分析了深井環(huán)境對突出模擬系統(tǒng)的要求,研制了深井煤巖瓦斯動(dòng)力災(zāi)害模擬試驗(yàn)系統(tǒng)。許江等[19-20]利用多場耦合煤礦動(dòng)力災(zāi)害大型物理模擬試驗(yàn)系統(tǒng),探究了不同巷道布置方式下煤粉–瓦斯兩相流的動(dòng)力過程;基于該裝置,彭守建等[21]進(jìn)一步揭示了突出煤體的運(yùn)移和沖擊特性。唐巨鵬等[22]研發(fā)了真三軸突出實(shí)驗(yàn)系統(tǒng),分析了突出孕育階段煤體破裂規(guī)律,探討了深部和淺部突出發(fā)生機(jī)理的差異。金侃等[23-24]以研究巷道內(nèi)氣–固兩相流為重點(diǎn)開展了突出煤粉–瓦斯兩相流實(shí)驗(yàn),指出突出煤–瓦斯流動(dòng)的破壞性和輸運(yùn)特性會(huì)因?yàn)槲酵咚沟目焖俳馕艿斤@著影響。張超林等[25]基于模塊化思路研制了多功能煤與瓦斯突出模擬試驗(yàn)系統(tǒng),并以河南龍山煤礦為背景進(jìn)行突出模擬試驗(yàn),結(jié)果表明該試驗(yàn)系統(tǒng)穩(wěn)定可靠,能準(zhǔn)確測定試件滲透率,并真實(shí)模擬再現(xiàn)突出的全過程。尹光志等[26-27]以煤體含水率對突出強(qiáng)度的影響為關(guān)注點(diǎn)開展了模擬試驗(yàn)研究,結(jié)果表明突出強(qiáng)度隨著含水率的增加而下降,2 者近似為二次曲線關(guān)系。張春華等[28]通過搭建煤與瓦斯突出模擬試驗(yàn)臺(tái)研究了石門揭煤引發(fā)的突出現(xiàn)象,分析了地質(zhì)構(gòu)造對突出過程和災(zāi)害防治的影響,結(jié)合模擬試驗(yàn)和數(shù)值模擬得到了含“構(gòu)造包體”煤層突出的原因和規(guī)律。

    王剛等[29]研發(fā)的石門揭煤突出模擬試驗(yàn)系統(tǒng)能夠全面考慮地應(yīng)力、瓦斯壓力和煤體結(jié)構(gòu)對突出現(xiàn)象的影響,同時(shí)能反映地應(yīng)力和瓦斯壓力的變化。郭品坤[30]設(shè)計(jì)了一種真三軸煤與瓦斯突出模擬試驗(yàn)系統(tǒng),并在此基礎(chǔ)上進(jìn)行了相似模擬試驗(yàn),研究了瓦斯壓力、氣體吸附能力、煤體斷裂韌度和滲透率等因素對突出過程的作用,提出突出發(fā)展階段可劃分為加速期、穩(wěn)定期和衰減期。LU Y 等[31]開發(fā)了一種多功能突出模擬試驗(yàn)系統(tǒng),能夠模擬煤巖體復(fù)雜應(yīng)力環(huán)境,并用于研究氣固耦合和氣液固多相耦合條件下的強(qiáng)度變形、滲透性和失穩(wěn)破壞特性。CAO J 等[32]設(shè)計(jì)了一種大型三維模擬試驗(yàn)裝置,能夠模擬突出的整個(gè)過程(包括孕育、激發(fā)、發(fā)展和終止階段),并對突出特征和機(jī)理進(jìn)行了全面研究。聶百勝等[33]開發(fā)了一種中等尺度的突出模擬裝置,該裝置能夠模擬接近采場工作面的開采環(huán)境,實(shí)現(xiàn)在不同地質(zhì)條件和瓦斯參數(shù)下的突出模擬;基于該裝置,還研究了突出過程中的能量演化特征,指出瓦斯膨脹能遠(yuǎn)大于煤體彈性能,是突出的主要能量來源[34]。

    總的來說,煤與瓦斯突出模擬試驗(yàn)裝置種類眾多,設(shè)計(jì)者綜合考慮了多種物理因素,能夠?qū)崿F(xiàn)不同條件下的突出模擬試驗(yàn)。對于許多突出模擬試驗(yàn),其基本方法如圖1 所示,可以簡化為:通過機(jī)械手段主動(dòng)打開突出口,使達(dá)到瓦斯吸附平衡的封閉煤體突然暴露,引發(fā)煤與瓦斯突出。因此,壓力閾值被廣泛應(yīng)用于突出條件判斷,即將突出發(fā)生的條件等同于煤中的瓦斯壓力是否超過特定壓力閾值。隨著科學(xué)技術(shù)的發(fā)展,氣體抽真空/充氣系統(tǒng)、數(shù)據(jù)收集系統(tǒng)(包括收集氣體壓力信號(hào)、微震信號(hào)、各種頻段的聲信號(hào)以及應(yīng)力變化等)、應(yīng)力加載系統(tǒng)、溫度控制系統(tǒng)、模擬巷道采掘系統(tǒng)等被逐漸應(yīng)用于突出的相似模擬試驗(yàn)中,它們旨在盡可能還原實(shí)際突出的物理環(huán)境并捕捉與突出過程相關(guān)的重要信息。

    圖1 實(shí)驗(yàn)室模擬突出的基本原理Fig.1 Fundamental principles for simulating an outburst in the laboratory

    2 突出模擬試驗(yàn)的相似性分析

    2.1 相似準(zhǔn)則

    對于物理相似模擬試驗(yàn),相似性是至關(guān)重要的,即需要確保研究的物理過程在模型和原型之間具有較好的一致性。2 個(gè)物理現(xiàn)象相似性包括幾何相似、運(yùn)動(dòng)相似和動(dòng)力相似等多個(gè)方面,以往絕大部分實(shí)驗(yàn)研究均以下述5 個(gè)方面的相似性為重點(diǎn)關(guān)注對象[7, 10-13, 32, 35-39]。

    (1)幾何相似。

    幾何相似是進(jìn)行模擬試驗(yàn)的基礎(chǔ)之一,它要求模型與原型中各種物理過程所處的物理空間幾何形狀相似。這種比例關(guān)系可表示為

    式中,Kl為幾何相似常數(shù);lp1、lp2為原型線性長度;lm1、lm2為模型線性長度。

    值得注意的是,幾何相似常數(shù)Kl通常是相似準(zhǔn)則中首先考慮和確定的,因?yàn)樗瞧渌鞣N物理縮尺關(guān)系轉(zhuǎn)換的基礎(chǔ)。例如,基于Kl,可確定面積相似常數(shù)KA=Kl2、體積相似常數(shù)KV=Kl3、時(shí)間相似常數(shù)。

    (2)運(yùn)動(dòng)相似。

    運(yùn)動(dòng)相似要求模型和原型中的物理質(zhì)點(diǎn)在相應(yīng)點(diǎn)上的速度和加速度方向相同且成比例,包含速度相似和加速度相似。

    速度相似:

    加速度相似:

    式中,Kv為速度相似常數(shù);Ka為加速度相似常數(shù);vp1、vp2為原型中質(zhì)點(diǎn)的速度,m/s;vm1、vm2為模型中質(zhì)點(diǎn)的速度,m/s;ap1、ap2為原型中質(zhì)點(diǎn)的加速度,m/s2;am1、am2為模型中質(zhì)點(diǎn)的加速度,m/s2。

    根據(jù)運(yùn)動(dòng)學(xué)規(guī)律可進(jìn)一步得到運(yùn)動(dòng)相似比(Kv和Ka)與幾何相似比(Kl)的關(guān)系為

    (3)動(dòng)力相似。

    動(dòng)力相似要求模型與原型中對應(yīng)質(zhì)點(diǎn)所受同種力(例如重力、黏滯力、彈性力等)的方向相同且大小成比例:

    式中,KF為動(dòng)力相似常數(shù);Fp1、Fp2為原型中質(zhì)點(diǎn)所受的力,N;Fm1、Fm2為模型中質(zhì)點(diǎn)所受的力,N。

    (4)非定常相似準(zhǔn)則。

    在突出過程中,煤體內(nèi)部的瓦斯流動(dòng)以及巷道中的煤粉–瓦斯流均屬于非定常流,因此突出模擬試驗(yàn)的Strouhal 數(shù)需要與實(shí)際突出相等[23-24]:

    式 中,Stp為 原 型 中 的Strouhal 數(shù);Stm為 模 型 中 的Strouhal 數(shù)。

    (5)相似材料。

    鑒于構(gòu)造煤與突出事故之間的密切關(guān)聯(lián)[1,40],構(gòu)造煤通常是突出模擬試驗(yàn)中的首選材料。然而,與結(jié)構(gòu)完整、抗擾動(dòng)能力強(qiáng)的原生煤不同,構(gòu)造煤的強(qiáng)度較低,黏聚力也較弱,因此很難直接從中獲取巖心試樣[4,41-42]。為了解決這個(gè)問題,通常使用粉碎的構(gòu)造煤顆粒(有時(shí)加入膠結(jié)劑)制備重構(gòu)巖心試樣,以便廣泛用于突出模擬試驗(yàn)。在這種情況下,實(shí)驗(yàn)材料與原始煤層的密度相似常數(shù)可表示為

    式中,Kρ為密度相似常數(shù);ρp1、ρp2為原型中質(zhì)點(diǎn)的密度,g/cm3;ρm1、ρm2為模型中質(zhì)點(diǎn)的密度,g/cm3。

    此時(shí),可獲得動(dòng)力相似比(KF)與運(yùn)動(dòng)相似比(Ka)的關(guān)系為

    式中,mp、mm分別為原型、模型中的質(zhì)點(diǎn)質(zhì)量,g。

    為了能夠更好地揭示煤與瓦斯突出過程的本質(zhì),以往學(xué)者們開展了大量研究以增加突出模擬試驗(yàn)與實(shí)際突出之間的相似性。例如,袁亮等[15-16]考慮了力學(xué)模型和能量模型建立了一種新的物理模擬準(zhǔn)則,并研發(fā)了高吸附含瓦斯煤相似材料;王漢鵬等[43]研制出一種新型的含瓦斯煤體相似材料(以煤粉為骨料,腐植酸鈉水溶液為膠結(jié)劑,混合壓制成型后干燥),并采用該相似材料制作的型煤成功模擬了突出過程,證明了該相似材料能較好地模擬含瓦斯煤體;張淑同等[44]通過統(tǒng)計(jì)分析煤與瓦斯突出現(xiàn)場資料及現(xiàn)有研究成果,得出了突出煤層和突出孔洞的特征,推導(dǎo)出了基于突出力學(xué)作用機(jī)理的相似準(zhǔn)則,研究了突出相似模擬系統(tǒng)和相似材料參數(shù),建立了固–氣耦合的數(shù)值模型,并進(jìn)一步研究了突出的發(fā)展過程和運(yùn)動(dòng)規(guī)律。

    2.2 幾何相似與相似材料

    在突出的相似模擬試驗(yàn)中,幾何相似與相似材料是被廣泛關(guān)注的2 個(gè)方向,它們是運(yùn)動(dòng)相似與動(dòng)力相似的基礎(chǔ)。總體而言,相似準(zhǔn)則中各相似常數(shù)的取值范圍尚未有統(tǒng)一的約束,不同模擬試驗(yàn)之間可能存在較大差異。

    在進(jìn)行突出模擬試驗(yàn)時(shí),通常使用基于構(gòu)造煤的重制型煤。在以往的試驗(yàn)中,型煤的密度相似常數(shù)Kρ為1.04~1.65,其控制因素主要是型煤的重制過程,尤其是成型壓力。使用構(gòu)造煤的重制型煤而非原始煤層的巖心取樣主要有2 個(gè)原因:

    (1)與結(jié)構(gòu)完整、抗擾動(dòng)能力強(qiáng)的原生煤相比,構(gòu)造煤強(qiáng)度較低、黏聚力較弱,難以直接獲得巖心試樣。為解決此問題,通常將原始構(gòu)造煤樣粉碎(有時(shí)加入膠結(jié)劑),并在一定應(yīng)力條件下壓制成重構(gòu)試樣。這是一種普遍做法,可有效獲取構(gòu)造煤的物理力學(xué)性質(zhì)參數(shù)。

    (2)構(gòu)造煤通常具有較高的瓦斯含量和瓦斯解吸速率,這使得構(gòu)造煤層的突出危險(xiǎn)性普遍大于原生煤層[4,45-46]。許多觀測結(jié)果也證實(shí)了構(gòu)造煤與突出之間存在密切關(guān)系[3,30,40,47]。

    因此,型煤的壓制通常是突出模擬試驗(yàn)的必備步驟。圖2 展示了一些實(shí)驗(yàn)中孔隙率與成型壓力之間的關(guān)系[16,38-39,48-52],可以看到,制樣過程中的成型壓力普遍小于80 MPa,這使得型煤孔隙率大多在10%~40%?;趫D2 中的實(shí)驗(yàn)數(shù)據(jù)可估計(jì)出,要使型煤孔隙率降低到1%,成型壓力可能至少需要256 MPa。顯然,這在目前的實(shí)驗(yàn)條件下是難以達(dá)到的,因此高孔隙率的型煤普遍被用于突出模擬試驗(yàn)。然而,需要指出的是,型煤壓制應(yīng)力的保壓時(shí)間、壓制次數(shù),煤體的幾何性質(zhì)(例如體積、形狀)以及力學(xué)參數(shù)(例如抗壓強(qiáng)度、彈性模量)等都可能影響成型應(yīng)力與成型孔隙率之間的關(guān)系[53-57]。事實(shí)上,由于以往的突出模擬試驗(yàn)中較少關(guān)注于型煤孔隙率,因此相當(dāng)缺乏孔隙率的分布數(shù)據(jù),型煤孔隙率與成型壓力之間的關(guān)系也仍未清晰。圖3展示了部分突出模擬試驗(yàn)中實(shí)驗(yàn)型煤的孔隙率分布[10-13,23,38-39,48,50-52,57-70],其中未指明孔隙率的試驗(yàn)研究參考圖2 中蔡成功[48]得出的成型壓力與孔隙率的關(guān)系式進(jìn)行估算,最終可以得到實(shí)驗(yàn)型煤的孔隙率中位數(shù)約為21.2%(這種估算是保守的,實(shí)際中可能會(huì)更高),這顯著大于實(shí)際煤層條件(通常1%~11%)[49-53]。

    圖2 成型壓力與孔隙率之間的關(guān)系[16,38-39,48-52]Fig.2 Relationship between molding pressure and briquette porosity[16,38-39,48-52]

    圖3 突出模擬試驗(yàn)中實(shí)驗(yàn)型煤的孔隙率范圍[10-13,23,38-39,48,50-52,57-70]Fig.3 Porosity range of experimental briquettes in outburst simulation tests[10-13,23,38-39,48,50-52,57-70]

    表1 匯總了1958—2020 年部分突出模擬裝置的關(guān)鍵參數(shù),這些實(shí)驗(yàn)腔體的體積為2.83×10-4~6.75 m3,橫跨4 個(gè)數(shù)量級(jí)。在現(xiàn)場巷道“寬×高”取4 m×3.5 m的情況下,以往大多數(shù)儀器的幾何相似常數(shù)Kl=10~100,而最近一些大型儀器的Kl=2~10。從實(shí)驗(yàn)結(jié)果上來看,不同尺度實(shí)驗(yàn)儀器的模擬結(jié)果總體上是互相支持的,其中大部分結(jié)論與實(shí)際發(fā)生的大尺度突出能建立較好的聯(lián)系。因此,尺度效應(yīng)對實(shí)驗(yàn)室模擬的影響并不明顯。此外,實(shí)驗(yàn)型煤的成型壓力、保壓時(shí)間、添加劑、粒度分布等沒有統(tǒng)一規(guī)定。以黏結(jié)劑為例,多數(shù)實(shí)驗(yàn)中采用質(zhì)量分?jǐn)?shù)為3%~6%的水,這有助于保持煤原有特性;當(dāng)對試驗(yàn)煤體的物理性質(zhì)有特殊要求時(shí),有機(jī)黏結(jié)劑(如腐植酸和聚乙烯醇,特點(diǎn)為強(qiáng)度低、灰分低)、無機(jī)黏結(jié)劑(如水泥和石膏,特點(diǎn)為強(qiáng)度高、灰分高)或復(fù)合黏結(jié)劑(2 種及2 種以上黏結(jié)劑混合)也常被用于型煤的壓制過程中[54-55]??傮w而言,突出模擬試驗(yàn)中型煤的抗拉強(qiáng)度通常(σt=0~1.2 MPa)顯著小于原生煤(σt=1~20 MPa),而相當(dāng)于構(gòu)造煤(σt=0.2~8.0 MPa)中“軟煤”,因此具有低強(qiáng)度的易突屬性[4,56]。

    表1 1958—2020 年部分突出模擬裝置的物理參數(shù)Table 1 Physical parameters of selected outburst simulation devices from 1958 to 2020

    2.3 初始力學(xué)條件的相似性

    在突出模擬試驗(yàn)中,由于重構(gòu)煤樣的密度(或孔隙率)往往偏大,這會(huì)導(dǎo)致其孔隙率(通常?=10%~40%)超過實(shí)際煤層的通常范圍(?=1%~11%),如圖4所示。同時(shí),實(shí)驗(yàn)型煤的滲透率(普遍超過10-10m2,1 mD=0.987×10-15m2)可能比現(xiàn)場測定的煤層滲透率(通常滲透率k=10-13~10-16m2)高出3 個(gè)量級(jí)以上[4,70-71]。因此,在開展突出模擬試驗(yàn)時(shí),型煤內(nèi)部幾乎沒有瓦斯壓力梯度,這與實(shí)際中煤層工作面前方的瓦斯壓力條件有很大差異,如圖5 所示。在這種情況下,實(shí)驗(yàn)室模擬突出通常依靠機(jī)械方式主動(dòng)打開突出口來誘發(fā)突出,即突然釋放密封空間內(nèi)的高壓瓦斯,利用瓦斯膨脹能驅(qū)動(dòng)突出的發(fā)展。

    圖4 實(shí)際煤層與實(shí)驗(yàn)型煤之間的孔隙差異Fig.4 Pore differences between actual coal seams and experimental briquettes

    實(shí)際采掘過程中遭遇突出通常是因?yàn)橛|碰到了工作面前方的高瓦斯壓力梯度區(qū)域。突出發(fā)生時(shí),煤體暴露表面的高瓦斯壓力梯度會(huì)導(dǎo)致煤體的破碎,破壞類型主要為拉伸破壞,此時(shí)煤體是否發(fā)生破碎取決于暴露表面附近的瓦斯壓力差(P0-P1)是否超過煤的抗拉強(qiáng)度(σt),其力學(xué)條件可表示[49,70,72]為

    式中,P1為環(huán)境壓力,MPa。

    雖然突出模擬試驗(yàn)與實(shí)際突出在能量來源上一致(均由瓦斯充當(dāng)驅(qū)動(dòng)力),但實(shí)驗(yàn)型煤的高滲透率性質(zhì)在一定程度上限制了利用物理模擬試驗(yàn)還原實(shí)際突出激發(fā)過程的能力[18]。具體來說,在突出模擬試驗(yàn)中,當(dāng)打開突出口使得煤體暴露后,可觀測到2 種不同結(jié)果:① 如果煤體暴露表面的瓦斯壓力差超過煤體的抗拉強(qiáng)度(P0-P1≥σt),則突出發(fā)生;② 如果煤體暴露表面的瓦斯壓力差低于煤體的抗拉強(qiáng)度(P0-P1<σt),則瓦斯壓力迅速降為環(huán)境壓力,不發(fā)生突出或僅發(fā)生小型的壓出(由實(shí)驗(yàn)應(yīng)力條件決定)。

    因此,在時(shí)間維度上,煤與瓦斯突出相似模擬的激發(fā)過程通常是瞬時(shí)完成的,即表現(xiàn)為煤中的瓦斯能和應(yīng)力能的同時(shí)釋放,而難以模擬真實(shí)的突出激發(fā)過程或孕育過程——這是由于實(shí)驗(yàn)設(shè)計(jì)(依靠機(jī)械方法主動(dòng)打開突出口)的不可避免結(jié)果[15,73-75]。相比之下,實(shí)際突出發(fā)生之前通常伴隨著采礦活動(dòng)以及煤層內(nèi)瓦斯壓力梯度的變化,因此其激發(fā)過程可能持續(xù)數(shù)小時(shí)、數(shù)天甚至更長時(shí)間,并伴隨著潛在的“預(yù)警信號(hào)”(或稱為“前兆信號(hào)”),如聲信號(hào)、微震信號(hào)和瓦斯異常等。這些與突出激發(fā)過程密切相關(guān)的信號(hào)在實(shí)驗(yàn)室中幾乎從未被捕獲。為了模擬真實(shí)的突出激發(fā)過程,近年來部分煤與瓦斯突出模擬試驗(yàn)已經(jīng)開始考慮模擬巷道開挖對壓力梯度的影響[14-15],但仍受限于實(shí)驗(yàn)型煤過高的孔隙率。尤其是,煤樣孔隙率的增加可能導(dǎo)致游離瓦斯量增加,從而進(jìn)一步影響突出潛能的聚集和消散過程。然而,這個(gè)現(xiàn)象以及背后潛在的問題目前尚未引起充分關(guān)注。接下來將對這個(gè)問題進(jìn)行深入研究。

    3 煤與瓦斯突出模擬試驗(yàn)的能量分析

    3.1 突出過程中的能量轉(zhuǎn)化特征

    突出發(fā)生后,煤層系統(tǒng)會(huì)經(jīng)歷從一個(gè)平衡態(tài)向另一個(gè)平衡態(tài)的轉(zhuǎn)變。以往學(xué)者們對煤與瓦斯突出的能量原理進(jìn)行了深入探討,研究重點(diǎn)包括揭示煤層系統(tǒng)能量變化特性以及突出各個(gè)階段的能量平衡[24,76-79]。普遍認(rèn)為,瓦斯膨脹能和煤體彈性能是煤與瓦斯突出的主要能量來源,這些能量導(dǎo)致了煤體的破碎和拋出;在發(fā)展階段中,瓦斯和破碎的煤/巖體會(huì)混合形成在巷道中快速運(yùn)動(dòng)的氣–固兩相流,該過程還存在著瓦斯殘余動(dòng)能、煤體摩擦、熱量、聲發(fā)射以及振動(dòng)等其他各種能量損耗。因此,總體而言,突出過程中的能量轉(zhuǎn)換關(guān)系[50,80-81]可概括為

    式中,A1為煤體彈性能,J;As為瓦斯膨脹能,J;W1為煤體搬運(yùn)功,J;W2為煤體破碎功,J;W3為其他能量消耗的總和,J。

    煤層釋放的能量與能量釋放區(qū)的體積密切相關(guān)[52,82];在實(shí)際中,該區(qū)域會(huì)隨著突出的發(fā)展而向煤層內(nèi)部移動(dòng),引起內(nèi)部煤體的破碎或損傷[30,52]。因此,煤層釋放的總能量At可進(jìn)一步表示為

    式中,E1為單位體積煤體具有的彈性能,J/m3;Es為單位體積煤體具有的瓦斯膨脹能,J/m3;Vs,1和Vs,2分別為煤體彈性能與瓦斯膨脹能的能量釋放區(qū)體積,m3。

    在大多數(shù)突出模擬試驗(yàn)中,Vs,1和Vs,2是可視為相等,且等價(jià)于試驗(yàn)腔體中煤體的總體積。在應(yīng)力條件下,單位體積煤體的彈性能[52,82]可表示為

    式中,E為彈性模量,MPa;μ為泊松比。

    對于瓦斯的熱力學(xué)過程,考慮到高壓瓦斯在煤孔隙中經(jīng)歷減壓并膨脹,因此等容和等壓變換不適用——主要的爭議在于突出過程究竟是更接近于絕熱過程還是等溫過程[24,38,50,77,83-84]。

    若突出過程中瓦斯與周邊煤巖體之間有充足的熱量交換使得瓦斯溫度基本保持恒定,即可認(rèn)為煤與瓦斯突出過程近似為等溫過程[52,85]:

    式中,WG,i為瓦斯等溫膨脹過程中釋放的能量,J;V0為突出總瓦斯體積,L。

    另一方面,若系統(tǒng)內(nèi)相關(guān)材料的導(dǎo)熱系數(shù)較低且熱力過程進(jìn)行得非常迅速,則瓦斯將幾乎不與周圍環(huán)境進(jìn)行熱交換,此時(shí)可假設(shè)突出過程是絕熱的。該過程中的瓦斯膨脹能[46,69-70,83]可表示為

    式中,WG,a為瓦斯絕熱膨脹過程中釋放的能量,J;γn為絕熱系數(shù)。

    無論突出過程是等溫還是絕熱,瓦斯膨脹能主要與瓦斯壓力以及瓦斯體積2 個(gè)因素相關(guān)。鑒于突出總瓦斯體積(V0)由孔隙中的游離瓦斯體積與煤中吸附瓦斯的解吸量組成,存在

    式中,Vf為初始煤孔隙中游離瓦斯的體積,m3;Vd為吸附瓦斯解吸后參與突出的體積,m3;Af為游離瓦膨脹能,J;Ad為解吸瓦斯膨脹能,J。

    在給定的煤層瓦斯壓力(P0)下,單位體積煤體中的游離瓦斯量(Vf)可根據(jù)氣體狀態(tài)方程獲得,即

    式中,T0為突出前的瓦斯溫度,K;?為孔隙率;T1為突出后的瓦斯溫度,K; ξ為氣體壓縮系數(shù),將煤中的瓦斯近似為理想氣體時(shí)可取 ξ =1。

    解吸瓦斯的參與量將在3.4 節(jié)中進(jìn)一步討論。根據(jù)式(17)可知,對于給定的環(huán)境條件,煤體孔隙率與游離瓦斯量呈正比關(guān)系,進(jìn)而控制著初始瓦斯膨脹能。需要指出的是,鑒于突出過程究竟是更接近絕熱(式(14))還是更接近等溫(式(13))仍存爭議,下文的計(jì)算中將基于這2 種計(jì)算方法得出的能量區(qū)間進(jìn)行分析討論。

    3.2 突出模擬試驗(yàn)的儀器與方案

    為了深入研究突出模擬試驗(yàn)中的能量釋放特征,參考了1988—2018 年突出模擬試驗(yàn)裝置的典型結(jié)構(gòu)特征與試驗(yàn)設(shè)計(jì)(圖6),使用自主研發(fā)的真三軸煤與瓦斯突出模擬試驗(yàn)系統(tǒng)進(jìn)行了一系列模擬試驗(yàn),旨在通過突出觀測數(shù)據(jù)來估算突出能量。

    圖6 1988—2018 年典型突出裝置的示意Fig.6 Schematic diagram of some typical outburst simulation instruments from 1988 to 2018

    圖7 展示了突出模擬系統(tǒng)的組成,主要包括:① 突出模擬試驗(yàn)腔體,尺寸為250 mm×250 mm×250 mm(長×寬×高);② 應(yīng)力加載系統(tǒng),垂直加載應(yīng)力≤80 MPa,左右及前后向水平應(yīng)力≤27 MPa;③ 注氣/真空抽氣系統(tǒng),使用ZJP-30 羅茨真空泵抽真空,并利用高靈敏度壓力表準(zhǔn)確調(diào)節(jié)氣體壓力;④ 恒溫控制系統(tǒng),采用電阻加熱帶包覆腔體的恒溫系統(tǒng),外層加保溫材料;⑤ 數(shù)據(jù)采集系統(tǒng),包括瓦斯壓力傳感器、高速攝影機(jī)、溫度傳感器和數(shù)據(jù)采集儀。

    圖7 真三軸煤與瓦斯突出模擬系統(tǒng)的構(gòu)造示意Fig.7 Test system diagram of triaxial coal and gas outburst simulation systems

    實(shí)驗(yàn)煤樣的基本參數(shù)見表2,煤樣原始孔隙率為4.1%,型煤成型孔隙率為20.08%,該范圍接近以往試驗(yàn)的中位數(shù)(圖3),密度相似常數(shù)Kρ約為1.25,其抗拉強(qiáng)度與較軟的構(gòu)造煤相當(dāng)。此外,為了探究氣體吸附/解吸差異對突出能量的影響,突出模擬試驗(yàn)選擇CO2(強(qiáng)吸附性氣體)、N2(弱吸附性氣體)、He(非吸附性氣體)作為實(shí)驗(yàn)氣體;基于實(shí)驗(yàn)室安全的考慮,未采用易燃易爆的CH4作為實(shí)驗(yàn)氣體,具體方案見表3。以CO2和N2等氣體替代CH4作為突出氣體的做法在以往的研究中是較為普遍的[10-13,34-36,57-58,67,86]。

    表2 實(shí)驗(yàn)煤樣基本物性參數(shù)Table 2 Basic physical parameters of the experimental coal sampe

    表3 煤與瓦斯突出模擬試驗(yàn)方案Table 3 Experimental scheme of coal and gas outburst tests

    開展模擬試驗(yàn)的關(guān)鍵步驟如下:

    ① 制備樣品。將原始煤樣粉碎至粒徑0.25~0.5 mm,加入6%水分后攪拌均勻。② 煤樣成型。將煤樣置于實(shí)驗(yàn)腔體中,用應(yīng)力加載系統(tǒng)壓制成型,壓力垂直應(yīng)力為60 MPa,保壓時(shí)間60 min。③ 腔體準(zhǔn)備。在確認(rèn)實(shí)驗(yàn)腔體密封后,抽真空24 h,然后注入實(shí)驗(yàn)氣體,在指定氣體壓力下保持72 h 以上,使實(shí)驗(yàn)型煤接近氣體飽和狀態(tài)。④ 突出前準(zhǔn)備。對實(shí)驗(yàn)型煤施加σx=σy=σz=5 MPa 的環(huán)境應(yīng)力,通過壓力傳感器檢查煤體內(nèi)部氣體壓力是否達(dá)到規(guī)定的突出壓力,并進(jìn)行調(diào)整。⑤ 啟動(dòng)突出。打開密封堵頭,誘發(fā)突出。⑥ 數(shù)據(jù)收集。利用高速攝影機(jī)記錄突出過程,獲取突出孔洞特征、突出強(qiáng)度、煤體堆積特性和粒度分布等突出特征。

    3.3 煤與瓦斯突出模擬試驗(yàn)結(jié)果

    在實(shí)踐中,瓦斯壓力被廣泛作為突出條件的判定參數(shù),即認(rèn)為突出的發(fā)生與瓦斯壓力超過某個(gè)閾值密切相關(guān)[4,87]。模擬試驗(yàn)的結(jié)果表明(表4),突出的瓦斯壓力閾值呈現(xiàn)隨實(shí)驗(yàn)氣體吸附能力增加而減小的趨勢,CO2為0.35~0.40 MPa、N2為0.40~0.45 MPa、He 為0.45~0.50 MPa,該現(xiàn)象可能與氣體吸附導(dǎo)致的煤強(qiáng)度降低密切相關(guān)[88-89]。

    表4 煤與瓦斯突出模擬試驗(yàn)結(jié)果Table 4 Results of the coal and gas outburst simulation experiments

    相對突出強(qiáng)度(Ro)通常被用來定量評(píng)價(jià)突出的破壞性[24,51,69,79],其定義為

    式中,Mo為突出煤體的質(zhì)量,即突出強(qiáng)度,kg;ML為裝入突出試驗(yàn)腔體的煤體總質(zhì)量,kg。

    結(jié)果表明,相對突出強(qiáng)度Ro隨著突出壓力的增加呈上升趨勢;在相同突出壓力下,CO2試驗(yàn)中的Ro分別是N2和He 的1.96~2.02 倍和1.80~1.82 倍。特別是,結(jié)果表明只有在吸附性氣體CO2和N2的突出模擬試驗(yàn)中可以觀測到煤體層裂現(xiàn)象(圖8(a)、(b)),而在非吸附性氣體He 的突出試驗(yàn)中則幾乎沒有層裂現(xiàn)象(圖8(c))。上述結(jié)果強(qiáng)調(diào)了吸附氣體對突出強(qiáng)度的顯著貢獻(xiàn)以及對層裂發(fā)展的決定性影響。

    為了獲得煤體在突出過程中的破碎特征,每次試驗(yàn)后均對現(xiàn)場突出煤體進(jìn)行了仔細(xì)收集和篩分。圖9展示了瓦斯壓力為0.5~0.7 MPa 時(shí)突出煤體的粒度分布,大多數(shù)煤體(質(zhì)量占比48.57%~78.84%)均在粒徑0.25~0.50 mm,僅有13.65%~17.83%的煤體被破碎到0.25 mm 以下。盡管小顆粒煤體的質(zhì)量占比較小,但以往研究普遍指出具有這些解吸速度極快的煤粉對于促進(jìn)突出發(fā)展有重要作用[24,79,90]。突出煤體的粒度分布數(shù)據(jù)將有助于在下面估算突出過程中的瓦斯解吸量。

    圖9 突出煤體的粒度分布Fig.9 Particle-size distribution of outburst coal

    為了獲得突出煤體的堆積特征,噴出的煤體在試驗(yàn)后均以1 m 為間隔被收集并稱重,結(jié)果如圖10 所示??梢园l(fā)現(xiàn),堆積的突出煤體表現(xiàn)為具有起伏的波浪形分布,這意味著突出過程中煤層的破裂傳播可能具有周期性特征,即暴露表面的層裂在高壓力梯度下規(guī)律性的形成并拋出,這種近似周期性的動(dòng)態(tài)破裂最終使得突出煤體出現(xiàn)局部質(zhì)量聚集,因此質(zhì)量隨突出距離的演化曲線呈現(xiàn)為波浪形。這種力學(xué)破壞的周期性本質(zhì)上是突出能量釋放的周期性,這將在第4 節(jié)中進(jìn)一步討論。

    圖10 突出煤體的質(zhì)量分布特征Fig.10 Mass distribution characteristics of outburst coal

    此外,突出煤體堆積特征(圖10)還展現(xiàn)了突出發(fā)展階段中氣體解吸對突出煤體搬運(yùn)的顯著貢獻(xiàn)。具體而言,與非吸附性氣體He 的突出試驗(yàn)相比,N2和CO2的試驗(yàn)中可以觀測到普遍較遠(yuǎn)的突出距離,并且大量的突出煤體聚集在距離突出口較遠(yuǎn)的位置,證明煤中吸附氣體的參與對煤粉–瓦斯流的運(yùn)動(dòng)起著促進(jìn)作用,導(dǎo)致了更強(qiáng)的搬運(yùn)能力。

    3.4 突出模擬試驗(yàn)中的能量估算

    通過突出模擬試驗(yàn)可以發(fā)現(xiàn),完全由游離氣體參與的He 試驗(yàn)也能發(fā)生突出——高孔隙率型煤為游離氣體提供了良好的空間條件,導(dǎo)致游離氣體膨脹能足以誘發(fā)突出,而不再依賴于吸附氣體的貢獻(xiàn)。下面將對其中的能量釋放特征做進(jìn)一步估算。

    煤孔隙中游離氣體的體積在已知壓力(P0)與孔隙率(?)的前提下可以通過式(17)進(jìn)行估算,而游離氣體膨脹能可以通過式(13)和式(14)進(jìn)行估算。然而,對于煤中吸附氣體的能量貢獻(xiàn),則需要在一定假設(shè)下進(jìn)行估算。煤中吸附氣體參與突出是一個(gè)較為復(fù)雜的過程,涉及到吸附氣體的脫附、在孔隙內(nèi)的擴(kuò)散以及在裂隙系統(tǒng)中的滲流等[88,91-94]。作為一個(gè)總的計(jì)算思路,假設(shè)在突出持續(xù)時(shí)間內(nèi),突出煤體與能量釋放區(qū)內(nèi)的煤體均對解吸氣體參與量有貢獻(xiàn)。在這種假設(shè)下,可能會(huì)高估解吸氣體實(shí)際參與突出的量,但是仍然可以作為與游離氣體參與量的對比。因此,基于煤體的解吸性能測試,參考突出壓力(表4)與突出煤體粒徑分布(圖9),通過計(jì)算試驗(yàn)腔體中煤體與突出煤體在突出持續(xù)時(shí)間內(nèi)的氣體解吸量來估計(jì)突出過程中的解吸氣體參與量Vd,可表示為

    式中,t0為突出持續(xù)時(shí)間,s;γi為給定粒度范圍內(nèi)煤體的質(zhì)量分?jǐn)?shù),%;vi為給定粒度范圍內(nèi)煤體在給定時(shí)間內(nèi)的平均解吸速率,m3/(t·s)。

    為了獲得式(19)中的平均解吸速率vi,開展了平衡壓力為0.5、0.6、0.7 MPa,煤樣粒徑為<0.074、0.074~0.25、0.25~0.50、0.5~1、>1 mm 的N2和CO2解吸實(shí)驗(yàn)。鑒于本文所有突出模擬試驗(yàn)的持續(xù)時(shí)間均<2.5 s,每個(gè)粒徑區(qū)間取前2.5 s 的平均解吸速率(vi),且均由3 次重復(fù)的解吸實(shí)驗(yàn)平均得到。

    計(jì)算參數(shù)列于表5,其中在計(jì)算質(zhì)量分?jǐn)?shù)γi(式(19))時(shí),考慮到突出試驗(yàn)腔體內(nèi)未突出的型煤具有較大的孔隙率,將這些煤體的解吸速率等價(jià)為粒徑 > 1 mm 的部分(實(shí)際上可能達(dá)不到該解吸速率),這種處理會(huì)高估解吸氣體的參與量。另一方面,在突出模擬試驗(yàn)中,施加的圍壓會(huì)擠壓煤體的大孔(減小孔隙空間),所以相對于無應(yīng)力條件下吸附/解吸實(shí)驗(yàn)獲得的結(jié)果,突出模擬試驗(yàn)中的煤體會(huì)具有較低瓦斯的解吸能力。因此,上述一系列估算是較為粗略的,其計(jì)算結(jié)果會(huì)在第4 節(jié)中進(jìn)一步與以往試驗(yàn)研究的計(jì)算結(jié)果進(jìn)行結(jié)合分析。

    表5 突出過程中氣體解吸量的計(jì)算參數(shù)Table 5 Parameters related to initial desorption amount of outburst

    突出強(qiáng)度取決于多種物理因素,包括氣體壓力、煤體性質(zhì)和環(huán)境應(yīng)力等[76,95-97]。利用表5 和式(19),可以估算出突出模擬試驗(yàn)中的氣體能量;煤體彈性能則可根據(jù)式(12)、施加的應(yīng)力條件(σx=σy=σz=5 MPa)以及表2 中的煤體力學(xué)參數(shù)進(jìn)行計(jì)算。結(jié)果表明(表6),對于同種氣體的試驗(yàn),相對突出強(qiáng)度隨著總突出能量(At)的增長近似線性的增大,如圖11 所示。在瓦斯膨脹能(As)中,解吸瓦斯的能量貢獻(xiàn)(Ad)占25.66%~66.02%;具體來說,吸附性氣體N2和CO2試驗(yàn)中釋放的氣體膨脹能比非吸附性氣體He的試驗(yàn)中的增大1.34~2.94 倍,表明了煤中吸附瓦斯解吸參與突出后對總瓦斯能量的顯著貢獻(xiàn)。另一方面,瓦斯膨脹能在總突出能量中(At)的占比為87.50%~95.31%(表6),證明了突出過程是氣體驅(qū)動(dòng)的動(dòng)力過程。值得一提的是,本文使用氣體壓力范圍(0.5~0.7 MPa)以及圍壓范圍(5 MPa)是符合相似準(zhǔn)則的,大量研究也表明了在更高的壓力或者圍壓下仍具有相似的能量分布特征[24,98-100]。

    表6 不同實(shí)驗(yàn)條件下的氣體能量及比較Table 6 Gas energy under different conditions and their comparison

    圖11 突出能量與相對突出強(qiáng)度之間的關(guān)系Fig.11 Correlation between the relative intensity of outburst and outburst energy

    4 討 論

    一些需要關(guān)注的現(xiàn)象是,突出模擬試驗(yàn)的突出壓力閾值通常遠(yuǎn)低于實(shí)際突出的壓力條件,普遍為0.1~0.7 MPa(本文為0.4~0.5 MPa)[45,101-104],如圖12所示。特別是,本文的研究結(jié)果還表明,在較低的應(yīng)力條件下(5 MPa),使用非吸附性氣體也能在較低的壓力下(≥0.4 MPa)誘發(fā)突出[45-46,87]。結(jié)合瓦斯膨脹能表達(dá)式(式(13)、(14))與游離瓦斯計(jì)算式(式(17))可知,實(shí)驗(yàn)型煤較高的孔隙率(10%~40%)不可避免地導(dǎo)致了較高的初始游離瓦斯膨脹能,進(jìn)而導(dǎo)致突出模擬試驗(yàn)可以不依賴于吸附瓦斯提供的突出能量。具體而言,當(dāng)孔隙率發(fā)生變化時(shí),游離瓦斯膨脹能也會(huì)相應(yīng)地發(fā)生變化:

    圖12 突出壓力閾值的范圍Fig.12 Range of outburst pressure thresholds

    式中,Ef,0為初始參考狀態(tài)下單位體積煤體中的游離瓦斯膨脹能,J/m3;?0為初始參考狀態(tài)下煤體孔隙率,%。

    通過比較原位條件(設(shè)?0=1%~11%)與實(shí)驗(yàn)型煤(設(shè)?=10%~40%)條件下的游離瓦斯膨脹能,可以發(fā)現(xiàn)游離瓦斯膨脹能會(huì)隨著實(shí)驗(yàn)型煤孔隙率的增大而顯著增大,甚至可能達(dá)到40 倍。一個(gè)直接的影響是,與實(shí)際突出相比,突出模擬試驗(yàn)中的游離瓦斯膨脹能占比普遍偏高。

    如圖13 所示,基于實(shí)際突出事故結(jié)果的能量分析表明,游離瓦斯膨脹能(Ef)在總瓦斯膨脹能中的占比(Rd)通常不超過20%,而該比例在突出模擬試驗(yàn)中則通常為30%~70%(與試驗(yàn)氣體相關(guān));對于非吸附性氣體He 誘發(fā)的突出,Rd則可視為100%。換言之,實(shí)際突出過程中的主要能量來源是煤中的吸附瓦斯,而突出模擬試驗(yàn)中則普遍由游離瓦斯主導(dǎo),甚至可以完全不依賴于吸附瓦斯解吸后釋放的能量。值得一提的是,以往模擬試驗(yàn)對Rd的估算結(jié)果(30%~70%)與本文的(25.66%~66.02%)幾乎一致,表明了3.4 節(jié)中近似計(jì)算方法引起的誤差在可被接受的范圍內(nèi)。

    若假設(shè)孔隙率變化不影響解吸瓦斯的能量貢獻(xiàn),那么結(jié)合式(13)、(14)和式(17)可得孔隙率變化時(shí)對游離瓦斯膨脹能(Ef)的影響,如圖14 所示。其中,能量閾值是基于突出壓力閾值(取0.4~0.5 MPa)進(jìn)行估算的,等效于成功進(jìn)行突出試驗(yàn)所需的總能量。結(jié)果表明,型煤高孔隙率引起的高初始游離瓦斯膨脹能是突出模擬試驗(yàn)可以不依賴于解吸瓦斯膨脹能的關(guān)鍵——在原始煤層的低孔隙率(通常小于11%)下,難以僅依賴孔隙中的游離瓦斯膨脹能誘發(fā)突出。

    圖14 煤體孔隙率對游離瓦斯膨脹能的影響Fig.14 Proportion of free gas energy in total gas energy

    實(shí)際上,大量觀測均表明突出的發(fā)生與煤體破碎導(dǎo)致的局部瓦斯聚集以及伴隨的瓦斯壓力變化密切相關(guān)[76,81,105],該突出激發(fā)過程的關(guān)鍵在于煤層破裂后為游離瓦斯膨脹能的增加創(chuàng)造了有利條件。在采掘過程中觀測到的突出前兆信號(hào)(如聲信號(hào)、微震和瓦斯流量異常等)也通常被認(rèn)為與工作面前方煤層中應(yīng)力集中引起的局部煤體破裂密切相關(guān)[69,79,100],表明突出的激發(fā)具有一定的延遲性。這些研究結(jié)果都強(qiáng)調(diào)了煤層局部破裂或變形引起的孔隙率增大對突出激發(fā)過程的促進(jìn)作用——孔隙空間增大導(dǎo)致了游離瓦斯膨脹能增大,進(jìn)而為突出激發(fā)提供能量來源。

    因此,從煤層局部破裂失穩(wěn)到突出的發(fā)生是一個(gè)能量持續(xù)積累的過程,而非瞬間完成的。從能量演化的角度來看,當(dāng)局部區(qū)域的應(yīng)力能(E1)釋放引起該區(qū)域的破裂或變形后,孔隙率(?)的增大對應(yīng)著游離瓦斯膨脹能(Ef)的增長——當(dāng)物理?xiàng)l件合適時(shí)(例如煤層滲透性較低導(dǎo)致瓦斯持續(xù)聚集時(shí)),局部能量變化表現(xiàn)為 ΔE1+ΔEf>0,此時(shí)能量演化過程表現(xiàn)為“儲(chǔ)能”,即總的突出潛能(Ef+E1)會(huì)不斷趨近于導(dǎo)致突出發(fā)生的能量閾值條件。對應(yīng)的,在某些環(huán)境下,局部應(yīng)力能(E1)的釋放也可能導(dǎo)致 ΔE1+ΔEf>0,即表現(xiàn)為“耗能”。在這種情況下,該區(qū)域中儲(chǔ)存的突出潛能(Ef+E1)總體上在減小,伴隨著突出危險(xiǎn)性的減小。

    在實(shí)際中,突出激發(fā)過程可能是儲(chǔ)能與耗能過程的不斷疊加。本質(zhì)上,突出的“延遲性”是煤體孔隙中的游離瓦斯在持續(xù)儲(chǔ)能的過程。早在20 世紀(jì)50 年代,蘇聯(lián)的Maкeeвкa科學(xué)研究所就發(fā)現(xiàn)很多突出事件是在震動(dòng)性放炮后幾分鐘或者幾小時(shí)后才發(fā)生的[6-7,106-107]。具體來說,1946—1953 年間的289 次突出中有32%延遲幾分鐘,5%延遲30 min 以上。在現(xiàn)場觀測中[1,5,23,52,89],廣泛報(bào)道了突出災(zāi)害前存在的“前兆信號(hào)”,包括但不限于小尺度的動(dòng)力現(xiàn)象、異響以及瓦斯涌出異常。上述結(jié)果均表明,突出的激發(fā)過程是一個(gè)能量逐步積累的過程。

    圖15 展示了模擬試驗(yàn)與實(shí)際突出中能量演化特征的對比。在實(shí)際煤層情況下(圖15(a)),由于煤層原始低孔隙率對游離瓦斯膨脹的限制,總突出潛能通常較低。因此,在未受外力擾動(dòng)前,煤層儲(chǔ)能一般低于突出的能量閾值。然而,當(dāng)煤層局部區(qū)域因外力作用而破裂時(shí),游離瓦斯量會(huì)隨著煤層孔隙率的增大而顯著增加。根據(jù)圖14 可知,單位體積煤體中的游離瓦斯膨脹能可能增加數(shù)倍甚至幾十倍(與孔隙率變化倍數(shù)相關(guān)),最終發(fā)生突出。

    圖15 突出過程中的能量演化原理Fig.15 Schematic diagram of the energy evolution in the outburst process

    然而,對于模擬試驗(yàn),由于重制型煤的高孔隙率對應(yīng)著很高的滲透率,通常無法完全模擬實(shí)際采掘過程(圖6)。因此,模擬試驗(yàn)通常依賴于通過主動(dòng)打開突出口的方式來誘發(fā)突出(圖2)。在這種情況下,模擬試驗(yàn)并不存在應(yīng)力能局部釋放和瓦斯能量逐漸累積的過程(圖15(b)),這在一定程度上也限制了實(shí)驗(yàn)室對前兆信號(hào)的研究。

    在某種程度上,突出模擬試驗(yàn)中對原始煤樣的取樣、粉碎、篩分和壓制等一系列導(dǎo)致裂隙空間增大的過程與實(shí)際突出之前的激發(fā)過程在結(jié)果上是等效的,即煤體的高孔隙率不可避免地導(dǎo)致了突出模擬試驗(yàn)在初始時(shí)就具有較高的游離瓦斯膨脹能,進(jìn)而可以不依賴于應(yīng)力條件以及吸附瓦斯的貢獻(xiàn)。

    5 結(jié) 論

    (1)突出模擬試驗(yàn)通常使用重制的型煤,重制型煤相對于實(shí)際煤層(通常?=1%~11%)而言具有普遍較高的孔隙率(?=10%~40%,中位數(shù)約為21.2%),這一方面導(dǎo)致了初始較高的游離瓦斯膨脹能(增大幾倍至幾十倍),另一方面導(dǎo)致了高滲透率(高出至少3 個(gè)量級(jí))。

    (2)突出模擬試驗(yàn)的結(jié)果表明,氣體膨脹能可占到總突出能量的87.50%~95.31%,其中初始游離氣體的貢獻(xiàn)占比為1/3~2/3。結(jié)合以往實(shí)驗(yàn)研究的結(jié)果,證明了突出模擬試驗(yàn)的本質(zhì),是將突出過程視為由煤中高壓瓦斯驅(qū)動(dòng)的動(dòng)力過程,即突出試驗(yàn)的成功開展可以不依賴于應(yīng)力條件,甚至可以利用非吸附性氣體(如氦氣)誘發(fā)。

    (3)由于重制型煤具有較高的滲透率,突出模擬試驗(yàn)通常無法還原實(shí)際突出的孕育過程以及激發(fā)過程,表現(xiàn)為瓦斯膨脹能與應(yīng)力能的同時(shí)釋放;相反,在實(shí)際煤層中,突出的發(fā)生往往依賴于局部應(yīng)力能釋放造成的局部煤體破裂,表現(xiàn)為突出激發(fā)的“延遲性”,即應(yīng)力能先于瓦斯膨脹能釋放,前者為后者的持續(xù)增長提供條件??傮w而言,突出模擬試驗(yàn)中型煤的高孔隙率嚴(yán)重限制了突出激發(fā)過程的相似性。

    (4)在未來的研究中,開發(fā)可以降低實(shí)驗(yàn)型煤孔隙率的壓制工藝有望讓模擬試驗(yàn)與實(shí)際突出具有更為接近的能量釋放特征,從而使突出模擬試驗(yàn)在相似性上有突破性的進(jìn)展;更重要的是,低滲透率型煤可以用于實(shí)驗(yàn)探究突出的激發(fā)過程,例如在實(shí)驗(yàn)室中捕捉和研究突出前兆信號(hào),進(jìn)而完善突出災(zāi)害的防治理論與技術(shù)。

    猜你喜歡
    型煤煤體游離
    型煤強(qiáng)度重復(fù)測試波動(dòng)情況統(tǒng)計(jì)分析
    中國煤炭(2024年2期)2024-03-11 13:24:26
    莫須有、蜿蜒、夜游離
    注熱井周圍煤體蠕變過程的滲透率變化規(guī)律模擬研究
    配煤及黏結(jié)劑配方對型煤強(qiáng)度的影響*
    正交試驗(yàn)下煤體滲透性影響因素評(píng)價(jià)
    煤炭工程(2019年4期)2019-05-05 02:30:16
    以“中央廚房”為突破口探索時(shí)政報(bào)道的融煤體之路——以浙江之聲為例
    長焰煤型煤復(fù)合粘結(jié)劑的研究
    超薄游離股前外側(cè)皮瓣修復(fù)足背軟組織缺損
    游離血紅蛋白室內(nèi)質(zhì)控物的制備及應(yīng)用
    游離于翻譯的精確與模糊之間——兼評(píng)第八屆CASIO杯翻譯競賽獲獎(jiǎng)譯文
    国产黄a三级三级三级人| 99久久九九国产精品国产免费| 日本免费一区二区三区高清不卡| 亚洲国产欧洲综合997久久,| 夜夜爽天天搞| 国产精品伦人一区二区| 国产伦在线观看视频一区| 大型黄色视频在线免费观看| 国产一区二区三区在线臀色熟女| 国产精品一区www在线观看| 国产综合懂色| 夫妻性生交免费视频一级片| 久久精品国产清高在天天线| 一个人观看的视频www高清免费观看| 日本黄色片子视频| 亚洲在久久综合| av在线观看视频网站免费| 国产成人福利小说| 日韩亚洲欧美综合| 啦啦啦韩国在线观看视频| 插逼视频在线观看| 婷婷色综合大香蕉| 亚洲最大成人av| 国产精品久久视频播放| 亚洲欧美日韩卡通动漫| 色综合色国产| 国产精品一区二区三区四区免费观看| 国产高清不卡午夜福利| 少妇熟女aⅴ在线视频| 国产成人freesex在线| 国产白丝娇喘喷水9色精品| 亚洲色图av天堂| 国产精品一区二区三区四区免费观看| 亚洲在久久综合| 大香蕉久久网| 免费无遮挡裸体视频| 亚洲人成网站在线播放欧美日韩| 亚洲中文字幕一区二区三区有码在线看| 一进一出抽搐gif免费好疼| 亚洲真实伦在线观看| 日韩精品有码人妻一区| 日韩欧美三级三区| 丰满乱子伦码专区| 狠狠狠狠99中文字幕| 欧美成人一区二区免费高清观看| 99在线视频只有这里精品首页| 欧美丝袜亚洲另类| 晚上一个人看的免费电影| 亚洲性久久影院| av天堂在线播放| 91久久精品国产一区二区三区| 国产精品一区www在线观看| 欧美最黄视频在线播放免费| 午夜福利在线在线| а√天堂www在线а√下载| 国产单亲对白刺激| 黄色一级大片看看| 深夜精品福利| 99久久精品国产国产毛片| 十八禁国产超污无遮挡网站| 国产爱豆传媒在线观看| 91久久精品国产一区二区三区| 热99re8久久精品国产| 赤兔流量卡办理| 日韩av不卡免费在线播放| 九九热线精品视视频播放| 99热6这里只有精品| kizo精华| 亚洲乱码一区二区免费版| 欧美性感艳星| 久久草成人影院| 久久久久国产网址| 久久中文看片网| 日本一二三区视频观看| 赤兔流量卡办理| 只有这里有精品99| 国产精品久久久久久av不卡| 黄色欧美视频在线观看| 精品人妻偷拍中文字幕| a级毛片a级免费在线| 一进一出抽搐动态| 成人高潮视频无遮挡免费网站| 少妇人妻一区二区三区视频| 亚洲精品久久国产高清桃花| 国产成人精品一,二区 | 一级毛片久久久久久久久女| 亚洲国产精品久久男人天堂| 寂寞人妻少妇视频99o| 天堂√8在线中文| av女优亚洲男人天堂| 成人欧美大片| 国产三级在线视频| 亚洲五月天丁香| 精品人妻一区二区三区麻豆| 亚洲国产高清在线一区二区三| 日韩三级伦理在线观看| 一级黄片播放器| 尾随美女入室| 婷婷色综合大香蕉| 国产老妇伦熟女老妇高清| 可以在线观看毛片的网站| 边亲边吃奶的免费视频| 搡老妇女老女人老熟妇| 国产 一区 欧美 日韩| 国产一级毛片七仙女欲春2| 亚洲av成人av| 久久精品国产亚洲网站| 波野结衣二区三区在线| 一级毛片电影观看 | 欧美另类亚洲清纯唯美| 色哟哟哟哟哟哟| 精品久久国产蜜桃| 国产91av在线免费观看| 99久久久亚洲精品蜜臀av| 丰满乱子伦码专区| 亚洲国产精品久久男人天堂| ponron亚洲| or卡值多少钱| 午夜精品国产一区二区电影 | 91在线精品国自产拍蜜月| 一卡2卡三卡四卡精品乱码亚洲| av天堂中文字幕网| 国产精品日韩av在线免费观看| 非洲黑人性xxxx精品又粗又长| 欧美人与善性xxx| 日韩精品有码人妻一区| 人妻少妇偷人精品九色| 伦理电影大哥的女人| 日日撸夜夜添| 精品人妻熟女av久视频| 国产精品国产高清国产av| 亚洲综合色惰| 可以在线观看的亚洲视频| 国产伦一二天堂av在线观看| 成人国产麻豆网| 99久久精品国产国产毛片| 中文字幕av在线有码专区| 美女大奶头视频| 精品久久久久久久末码| 国产午夜精品一二区理论片| 少妇被粗大猛烈的视频| av女优亚洲男人天堂| 色5月婷婷丁香| avwww免费| 亚洲欧美清纯卡通| 一个人看视频在线观看www免费| 热99在线观看视频| 中文字幕精品亚洲无线码一区| 搞女人的毛片| 久久中文看片网| 春色校园在线视频观看| 国产精品野战在线观看| 国产高清激情床上av| 日韩亚洲欧美综合| 不卡视频在线观看欧美| 久久人人爽人人爽人人片va| 亚洲av电影不卡..在线观看| 中文亚洲av片在线观看爽| 国产一区二区在线av高清观看| 日本黄色片子视频| 国产精品一区二区性色av| 美女xxoo啪啪120秒动态图| 亚洲美女搞黄在线观看| 国产亚洲5aaaaa淫片| 好男人视频免费观看在线| 啦啦啦观看免费观看视频高清| 国产极品天堂在线| 久久九九热精品免费| kizo精华| 高清日韩中文字幕在线| 蜜桃亚洲精品一区二区三区| 国内少妇人妻偷人精品xxx网站| 日本爱情动作片www.在线观看| 中文在线观看免费www的网站| 日韩人妻高清精品专区| 午夜老司机福利剧场| av专区在线播放| av国产免费在线观看| 18禁裸乳无遮挡免费网站照片| 国产伦精品一区二区三区四那| 哪个播放器可以免费观看大片| 一边摸一边抽搐一进一小说| av在线亚洲专区| www.色视频.com| 欧美xxxx性猛交bbbb| 岛国毛片在线播放| 久久久久久久午夜电影| 高清在线视频一区二区三区 | 日韩中字成人| 欧美色视频一区免费| 特级一级黄色大片| 亚洲欧洲日产国产| 日本在线视频免费播放| av天堂在线播放| 欧美成人一区二区免费高清观看| 久久久久久久久久成人| 国产精品三级大全| 天堂影院成人在线观看| 精华霜和精华液先用哪个| 男人和女人高潮做爰伦理| 欧美激情国产日韩精品一区| 亚洲av成人av| 亚洲久久久久久中文字幕| 边亲边吃奶的免费视频| 日韩高清综合在线| 日本熟妇午夜| 九色成人免费人妻av| 欧美成人免费av一区二区三区| 亚洲欧美精品自产自拍| 99国产精品一区二区蜜桃av| 中文资源天堂在线| АⅤ资源中文在线天堂| 亚洲美女搞黄在线观看| 亚洲国产高清在线一区二区三| 亚洲18禁久久av| 亚洲av男天堂| 寂寞人妻少妇视频99o| 久久精品夜夜夜夜夜久久蜜豆| 欧美一区二区亚洲| 久久99热6这里只有精品| 亚洲三级黄色毛片| 亚洲中文字幕一区二区三区有码在线看| 国产精品久久久久久精品电影| 欧美成人一区二区免费高清观看| 别揉我奶头 嗯啊视频| 成人高潮视频无遮挡免费网站| 国产一区二区在线观看日韩| 色播亚洲综合网| 久久精品国产鲁丝片午夜精品| av福利片在线观看| 久久久久久伊人网av| 一本久久精品| 成人美女网站在线观看视频| 三级国产精品欧美在线观看| 亚洲国产欧洲综合997久久,| 国产美女午夜福利| 在线国产一区二区在线| 中文欧美无线码| 国产一区二区在线av高清观看| 亚洲美女搞黄在线观看| 国产一区二区三区av在线 | 亚洲成人av在线免费| 听说在线观看完整版免费高清| 久久久久久久久久久丰满| 最好的美女福利视频网| 欧美色视频一区免费| 网址你懂的国产日韩在线| 国产精品日韩av在线免费观看| 99久久成人亚洲精品观看| 日本五十路高清| 日本一二三区视频观看| www.av在线官网国产| 男女边吃奶边做爰视频| 久久久久久久久久久丰满| 国产精品福利在线免费观看| 亚洲欧美中文字幕日韩二区| 国产极品精品免费视频能看的| 女人十人毛片免费观看3o分钟| 看片在线看免费视频| 久久久久久伊人网av| 成人二区视频| 日日啪夜夜撸| 一本一本综合久久| 高清午夜精品一区二区三区 | 国产极品精品免费视频能看的| 久久6这里有精品| 赤兔流量卡办理| 在线免费观看的www视频| 国产单亲对白刺激| 国产精品福利在线免费观看| 99热全是精品| 一区二区三区高清视频在线| 国产亚洲精品久久久com| 岛国毛片在线播放| 两个人视频免费观看高清| 久久久久久久久中文| 国产精品人妻久久久影院| 晚上一个人看的免费电影| 久久久午夜欧美精品| 免费观看的影片在线观看| 成年版毛片免费区| 丝袜美腿在线中文| 岛国毛片在线播放| 日日干狠狠操夜夜爽| 欧美高清成人免费视频www| av又黄又爽大尺度在线免费看 | 一级毛片aaaaaa免费看小| 99热网站在线观看| 别揉我奶头 嗯啊视频| 国产一区二区三区av在线 | 久久精品久久久久久噜噜老黄 | 久久久久久久久久成人| 久久综合国产亚洲精品| 尾随美女入室| 嫩草影院精品99| 日韩视频在线欧美| 中文精品一卡2卡3卡4更新| 精品99又大又爽又粗少妇毛片| 亚洲真实伦在线观看| 国产人妻一区二区三区在| 熟女电影av网| 午夜老司机福利剧场| 国产毛片a区久久久久| 最近视频中文字幕2019在线8| 欧美最黄视频在线播放免费| 国产av不卡久久| 久久精品国产鲁丝片午夜精品| 美女高潮的动态| 熟女电影av网| 性色avwww在线观看| 欧美+日韩+精品| 久久久久久大精品| 国产在线男女| 精品久久国产蜜桃| 精品久久久久久久久久免费视频| 久久中文看片网| 国产精品永久免费网站| 26uuu在线亚洲综合色| 桃色一区二区三区在线观看| 久久99蜜桃精品久久| 亚洲人与动物交配视频| 亚洲精品乱码久久久久久按摩| 亚洲av免费高清在线观看| 欧美区成人在线视频| 欧美zozozo另类| av专区在线播放| 91aial.com中文字幕在线观看| 久久久久久久午夜电影| 五月伊人婷婷丁香| 国产黄片视频在线免费观看| 美女 人体艺术 gogo| 91狼人影院| 亚洲欧洲国产日韩| 日本欧美国产在线视频| 91精品一卡2卡3卡4卡| 久久这里只有精品中国| 麻豆成人午夜福利视频| 内射极品少妇av片p| 丰满人妻一区二区三区视频av| 韩国av在线不卡| 在线a可以看的网站| 亚洲精品国产成人久久av| 亚洲图色成人| 亚洲激情五月婷婷啪啪| 日日摸夜夜添夜夜爱| 大又大粗又爽又黄少妇毛片口| av国产免费在线观看| 国产真实乱freesex| 国产精品.久久久| 国产成人精品婷婷| 女的被弄到高潮叫床怎么办| 高清午夜精品一区二区三区 | 国产精品乱码一区二三区的特点| 岛国在线免费视频观看| 三级经典国产精品| 美女内射精品一级片tv| 精品不卡国产一区二区三区| 午夜福利在线观看吧| .国产精品久久| 夜夜爽天天搞| 美女脱内裤让男人舔精品视频 | 看免费成人av毛片| 99久久无色码亚洲精品果冻| 亚洲成人av在线免费| 97热精品久久久久久| 婷婷精品国产亚洲av| 好男人视频免费观看在线| 国产精品一区二区在线观看99 | 久久精品国产清高在天天线| 99热这里只有是精品50| 亚洲精品粉嫩美女一区| 久久久欧美国产精品| 中文字幕熟女人妻在线| 国产伦一二天堂av在线观看| 日本爱情动作片www.在线观看| 国产探花极品一区二区| 91久久精品国产一区二区成人| 大香蕉久久网| 在线观看午夜福利视频| 日本黄色视频三级网站网址| 日本一二三区视频观看| 欧美精品一区二区大全| 麻豆乱淫一区二区| 97人妻精品一区二区三区麻豆| 久久久国产成人精品二区| 中文欧美无线码| 中出人妻视频一区二区| 国产精品av视频在线免费观看| 99久久人妻综合| 日日撸夜夜添| 联通29元200g的流量卡| 久久久久久久午夜电影| 麻豆成人av视频| 亚洲av第一区精品v没综合| 国内精品一区二区在线观看| 欧美成人一区二区免费高清观看| 欧美三级亚洲精品| 久久久国产成人免费| 欧美精品一区二区大全| 少妇熟女欧美另类| 五月伊人婷婷丁香| 级片在线观看| 99在线人妻在线中文字幕| 可以在线观看毛片的网站| 欧美一区二区国产精品久久精品| 国产亚洲精品av在线| 午夜亚洲福利在线播放| 老熟妇乱子伦视频在线观看| 亚洲精品国产av成人精品| 精品一区二区三区人妻视频| 秋霞在线观看毛片| 欧美日韩精品成人综合77777| 亚洲av电影不卡..在线观看| 国产精品免费一区二区三区在线| 亚洲三级黄色毛片| 黄色视频,在线免费观看| 伦精品一区二区三区| 亚洲七黄色美女视频| 国语自产精品视频在线第100页| 亚洲精品久久国产高清桃花| 99久久久亚洲精品蜜臀av| 热99re8久久精品国产| 亚洲欧美精品专区久久| 国产久久久一区二区三区| 色综合亚洲欧美另类图片| 亚洲人成网站高清观看| 国产成人福利小说| 伦精品一区二区三区| 欧美在线一区亚洲| 小说图片视频综合网站| 亚洲国产色片| 只有这里有精品99| 亚洲精品乱码久久久v下载方式| 国产蜜桃级精品一区二区三区| 亚洲人成网站高清观看| 狂野欧美激情性xxxx在线观看| 欧美人与善性xxx| 色播亚洲综合网| 在线天堂最新版资源| 99在线视频只有这里精品首页| 神马国产精品三级电影在线观看| 国产精品伦人一区二区| 亚洲av.av天堂| 欧美最黄视频在线播放免费| 欧美一区二区亚洲| 变态另类成人亚洲欧美熟女| 99在线视频只有这里精品首页| 欧美又色又爽又黄视频| 精品99又大又爽又粗少妇毛片| 我的女老师完整版在线观看| 日韩一区二区视频免费看| 深爱激情五月婷婷| 伦理电影大哥的女人| 美女cb高潮喷水在线观看| 久久久久久九九精品二区国产| 国产精品久久久久久av不卡| 插逼视频在线观看| 在线免费观看不下载黄p国产| 欧美日韩在线观看h| 国产91av在线免费观看| 欧美精品国产亚洲| 亚洲,欧美,日韩| 六月丁香七月| 午夜福利高清视频| 蜜桃亚洲精品一区二区三区| 国产精品免费一区二区三区在线| 久久久欧美国产精品| 亚洲av熟女| 国产精品久久久久久久电影| 少妇熟女欧美另类| 日产精品乱码卡一卡2卡三| 简卡轻食公司| 国产精品乱码一区二三区的特点| 欧美三级亚洲精品| 国产高清有码在线观看视频| 国产精品电影一区二区三区| 美女国产视频在线观看| 欧美成人免费av一区二区三区| 亚洲成人精品中文字幕电影| 精品久久久久久成人av| 99在线人妻在线中文字幕| 亚洲性久久影院| 欧美成人一区二区免费高清观看| 在线观看av片永久免费下载| 午夜精品国产一区二区电影 | 亚洲成人久久爱视频| 91aial.com中文字幕在线观看| 国产av一区在线观看免费| 国产精品三级大全| 色综合色国产| 亚洲欧美中文字幕日韩二区| 精品久久久久久久久亚洲| 一本精品99久久精品77| 91久久精品国产一区二区三区| 国产亚洲av片在线观看秒播厂 | 亚洲av中文av极速乱| 晚上一个人看的免费电影| 97超视频在线观看视频| kizo精华| 午夜老司机福利剧场| 国产三级中文精品| 欧美高清性xxxxhd video| 99在线人妻在线中文字幕| 不卡视频在线观看欧美| 99久国产av精品国产电影| 村上凉子中文字幕在线| 乱码一卡2卡4卡精品| 中文字幕熟女人妻在线| 亚洲成人av在线免费| 少妇人妻精品综合一区二区 | 免费看av在线观看网站| 99热网站在线观看| 亚洲欧美成人精品一区二区| 最后的刺客免费高清国语| 神马国产精品三级电影在线观看| 久久久久久久久久久丰满| 久久久久久久午夜电影| 波野结衣二区三区在线| 99热6这里只有精品| 99久久精品热视频| 亚洲精华国产精华液的使用体验 | 国产精品综合久久久久久久免费| 国产v大片淫在线免费观看| 久久婷婷人人爽人人干人人爱| 精品国内亚洲2022精品成人| 亚洲av熟女| 我的老师免费观看完整版| 极品教师在线视频| 一边亲一边摸免费视频| 欧美不卡视频在线免费观看| 亚洲国产精品成人久久小说 | 日韩大尺度精品在线看网址| av在线播放精品| 日本黄大片高清| 91精品国产九色| 长腿黑丝高跟| 国内精品一区二区在线观看| 成人性生交大片免费视频hd| 亚洲欧洲国产日韩| 在线国产一区二区在线| 99riav亚洲国产免费| av天堂在线播放| 一本精品99久久精品77| 青青草视频在线视频观看| av.在线天堂| 免费看日本二区| 久久精品影院6| 大香蕉久久网| 最近的中文字幕免费完整| 九九久久精品国产亚洲av麻豆| 看片在线看免费视频| 精品久久久久久久末码| av在线观看视频网站免费| 十八禁国产超污无遮挡网站| 亚洲av第一区精品v没综合| 两个人的视频大全免费| 欧美高清性xxxxhd video| 亚洲精品国产成人久久av| 欧美日韩综合久久久久久| 成年免费大片在线观看| 亚洲一级一片aⅴ在线观看| 国产成人一区二区在线| 1024手机看黄色片| 男女下面进入的视频免费午夜| 可以在线观看毛片的网站| 六月丁香七月| 日本与韩国留学比较| 亚洲av成人精品一区久久| 69人妻影院| 亚洲成人久久性| 欧洲精品卡2卡3卡4卡5卡区| 中国国产av一级| 永久网站在线| 欧美另类亚洲清纯唯美| 日韩强制内射视频| 天堂中文最新版在线下载 | 91aial.com中文字幕在线观看| 精品久久久噜噜| 一区二区三区四区激情视频 | 99热这里只有精品一区| or卡值多少钱| 亚洲第一区二区三区不卡| 亚洲av电影不卡..在线观看| 亚洲精品国产成人久久av| 欧美日韩一区二区视频在线观看视频在线 | 国产日本99.免费观看| 97超视频在线观看视频| 国内精品一区二区在线观看| 美女 人体艺术 gogo| 久久久久久久久久黄片| 国产毛片a区久久久久| 久久人人爽人人片av| 美女黄网站色视频| 伊人久久精品亚洲午夜| 狂野欧美白嫩少妇大欣赏| av免费在线看不卡| 成人高潮视频无遮挡免费网站| 51国产日韩欧美| 日韩欧美三级三区| 大香蕉久久网| 亚洲成a人片在线一区二区| 亚州av有码| 看免费成人av毛片| 韩国av在线不卡| 亚洲av成人av| 在线a可以看的网站| 国产精品伦人一区二区| 色综合站精品国产| 亚洲第一区二区三区不卡| 免费黄网站久久成人精品| 美女xxoo啪啪120秒动态图| 一区二区三区高清视频在线| 中文在线观看免费www的网站| 一区二区三区四区激情视频 | 国产乱人偷精品视频|