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

    高瓦斯煤層沖擊地壓發(fā)生條件與影響因素1)

    2017-07-03 15:00:25尹萬(wàn)蕾潘一山李忠華
    力學(xué)學(xué)報(bào) 2017年3期
    關(guān)鍵詞:摩擦角煤體塑性

    尹萬(wàn)蕾 潘一山 李忠華

    ?(遼寧工程技術(shù)大學(xué)力學(xué)與工程學(xué)院沖擊地壓研究院,遼寧阜新123000)?(遼寧大學(xué)物理學(xué)院,沈陽(yáng)110036)

    高瓦斯煤層沖擊地壓發(fā)生條件與影響因素1)

    尹萬(wàn)蕾?,2)潘一山?,?,3)李忠華?

    ?(遼寧工程技術(shù)大學(xué)力學(xué)與工程學(xué)院沖擊地壓研究院,遼寧阜新123000)?(遼寧大學(xué)物理學(xué)院,沈陽(yáng)110036)

    針對(duì)高瓦斯煤層沖擊地壓?jiǎn)栴},用解析方法得到?jīng)_擊地壓發(fā)生條件,分析了主要影響因素對(duì)滿(mǎn)足沖擊地壓發(fā)生條件的臨界塑性區(qū)半徑和臨界應(yīng)力的影響規(guī)律.結(jié)合五龍礦開(kāi)采實(shí)際情況對(duì)影響高瓦斯煤層沖擊地壓的煤的模量比、煤層瓦斯孔隙壓力、支護(hù)應(yīng)力和內(nèi)摩擦角4個(gè)因素做了對(duì)比分析.研究發(fā)現(xiàn):高瓦斯煤層在巷道掘進(jìn)面附近由于存在開(kāi)挖面空間效應(yīng),掘進(jìn)面前方尚未開(kāi)挖的煤體對(duì)巷道變形起到了限制作用,減少了沖擊地壓的發(fā)生,隨著掘進(jìn)面向前推進(jìn),后方一定距離范圍內(nèi)的巷道支護(hù)應(yīng)力增大.隨著瓦斯解吸滲流的進(jìn)行,巷道壁處孔隙壓力降低,巷道沖擊地壓危險(xiǎn)性明顯提高,此時(shí)提高支護(hù)應(yīng)力,沖擊危險(xiǎn)性有所降低.高瓦斯煤層巷道發(fā)生沖擊地壓的臨界塑性區(qū)半徑和臨界應(yīng)力隨模量比、瓦斯孔隙壓力的增大而快速減小,隨支護(hù)應(yīng)力的增大而增大,臨界塑性區(qū)半徑隨內(nèi)摩擦角的增大而增大,臨界應(yīng)力與內(nèi)摩擦角不是單調(diào)函數(shù)關(guān)系,存在一個(gè)極小值點(diǎn),當(dāng)內(nèi)摩擦角小于此極小值時(shí),臨界應(yīng)力隨內(nèi)摩擦角增大而減??;當(dāng)內(nèi)摩擦角大于此極小值時(shí),臨界應(yīng)力隨內(nèi)摩擦角增大而增大.

    沖擊地壓,高瓦斯煤層,模量比,孔隙壓力,支護(hù)應(yīng)力,內(nèi)摩擦角

    引言

    對(duì)沖擊地壓發(fā)生理論的研究已有近百年歷史,產(chǎn)生了強(qiáng)度理論、剛度理論、能量理論、沖擊傾向性理論、失穩(wěn)理論和組合理論.業(yè)內(nèi)學(xué)者對(duì)沖擊地壓的發(fā)生機(jī)理與類(lèi)型劃分[14]進(jìn)行了大量研究,但這些理論研究很少涉及煤層瓦斯對(duì)沖擊地壓發(fā)生的影響作用.目前,中國(guó)高瓦斯煤層發(fā)生沖擊地壓的礦井已不在少數(shù),如撫順老虎臺(tái)礦、阜新五龍礦和王營(yíng)礦、邯鄲陶二礦、豐城建新礦、平頂山十礦和十二礦等.但對(duì)高瓦斯煤層沖擊地壓的專(zhuān)門(mén)研究卻較少.潘一山等[58]最先系統(tǒng)地分析了瓦斯對(duì)煤體失穩(wěn)破壞的影響,建立了瓦斯煤層沖擊地壓數(shù)學(xué)模型;李鐵等[910]通過(guò)大量微震、瓦斯監(jiān)測(cè)和現(xiàn)場(chǎng)調(diào)查,發(fā)現(xiàn)沖擊地壓震前、震后和同震都伴隨瓦斯異常涌出的現(xiàn)象,對(duì)以往“沖擊地壓只是忽略或沒(méi)有瓦斯作用的煤體突出”觀點(diǎn)提出質(zhì)疑,認(rèn)為煤炭深部開(kāi)采條件下沖擊地壓與瓦斯密切相關(guān),高壓瓦斯氣體極有可能參與了沖擊地壓的孕育,存在一種含氣多孔介質(zhì)和儲(chǔ)氣構(gòu)造在開(kāi)挖卸荷和高壓吸附瓦斯解吸膨脹耦合作用下誘發(fā)的沖擊地壓.王振[11]分析了瓦斯對(duì)煤體力學(xué)性質(zhì)的影響以及瓦斯?jié)B流和裂隙擴(kuò)展的規(guī)律,并研究了沖擊地壓與瓦斯突出在不同階段相互誘發(fā)轉(zhuǎn)化的條件.馬海峰[12]研究了應(yīng)力場(chǎng)對(duì)瓦斯場(chǎng)的影響,并建立了氣固耦合模型.董飛亞[13]研究了高瓦斯煤層沖擊地壓發(fā)生機(jī)理.王磊[14]研究了應(yīng)力場(chǎng)和瓦斯場(chǎng)的耦合效應(yīng).徐晨陽(yáng)[15]總結(jié)了高瓦斯煤層沖擊地壓的特征和機(jī)理,并分析了研究中存在的問(wèn)題.王濤等[16]從多角度、多因素方面分析沖擊地壓的發(fā)生引起礦體震動(dòng)對(duì)瓦斯的吸附影響,并分析了涌出的原因與條件.梁冰等[1718]根據(jù)煤巖變形與瓦斯?jié)B流的作用機(jī)理,提出了煤體與瓦斯突出的固流耦合失穩(wěn)理論,建立了含瓦斯煤體的本構(gòu)關(guān)系,同時(shí)給出了材料參數(shù)的確定方法.俞善炳[1920]建立了煤與瓦斯突出的一維流動(dòng)模型,給出了破碎啟動(dòng)準(zhǔn)則,并討論了突出啟動(dòng)過(guò)程,研究了大型突出相應(yīng)于恒穩(wěn)推進(jìn)的情況,分析了煤與瓦斯突出的重要無(wú)量綱參數(shù)以及判據(jù).趙陽(yáng)升[21]研究了沖擊地壓與突出的統(tǒng)一機(jī)制,煤體沖擊傾向、孔隙瓦斯壓力及圍巖應(yīng)力對(duì)突出的影響,提出了突出的數(shù)學(xué)模型.藍(lán)航等[22]統(tǒng)計(jì)了近年來(lái)我國(guó)沖擊地壓發(fā)生的區(qū)域、條件、特點(diǎn)以及防治手段和效果,分析了地質(zhì)因素與開(kāi)采技術(shù)因素的影響.崔乃鑫等[23]研究了含瓦斯煤層實(shí)施鉆屑法過(guò)程中瓦斯的影響,采用經(jīng)典理論推得了檢測(cè)含瓦斯煤層沖擊地壓的鉆屑量指標(biāo).彭蘇萍等[24]根據(jù)煤層瓦斯與常規(guī)砂巖賦存天然氣機(jī)理對(duì)比結(jié)果,提出了以煤層裂隙為探測(cè)目標(biāo)的煤層瓦斯富集AVO技術(shù)預(yù)測(cè)理論.梁盛開(kāi)等[25]針對(duì)煤與瓦斯突出的預(yù)測(cè)問(wèn)題提出了多因素綜合評(píng)價(jià)的新途徑.宋真龍等[26]分析了不同瓦斯壓力環(huán)境中煤樣單軸壓縮與循環(huán)加載過(guò)程中能量集聚與耗散關(guān)系,得出高瓦斯礦井深部開(kāi)采進(jìn)行沖擊傾向性評(píng)價(jià)需考慮瓦斯因素的結(jié)論.劉??h等[27]采用非線性理論——突變理論,研究煤與瓦斯突出機(jī)理,得到延期突出的滯后現(xiàn)象是由地應(yīng)力、煤體中的瓦斯、煤的物理力學(xué)特性及外力作用(爆破)等因素共同作用的結(jié)果.張志剛等[28]建立了瓦斯氣體滲流的非線性滲流方程,結(jié)合試驗(yàn)驗(yàn)證得到含瓦斯煤體內(nèi)非線性附加阻力來(lái)源于基質(zhì)膨脹、邊界層影響、吸附的動(dòng)態(tài)傳質(zhì)過(guò)程的結(jié)論.郭為等[29]通過(guò)頁(yè)巖基質(zhì)解吸--擴(kuò)散--滲流耦合實(shí)驗(yàn),利用有限差分法求解解吸--擴(kuò)散--滲流耦合數(shù)學(xué)模型,對(duì)基巖氣體流動(dòng)影響因素進(jìn)行了分析,得到了頁(yè)巖基質(zhì)的滲透率、擴(kuò)散系數(shù)、解吸附常數(shù)等影響因素對(duì)基質(zhì)氣體的流量和壓力傳播規(guī)律.劉文超等[30]利用相似變量變換方法和基于空間坐標(biāo)變換的有限差分法研究啟動(dòng)壓力梯度的低滲透多孔介質(zhì)非達(dá)西滲流模型邊界問(wèn)題,對(duì)考慮內(nèi)邊界變壓力情況下啟動(dòng)壓力梯度一維低滲透多孔介質(zhì)非達(dá)西滲流動(dòng)邊界模型進(jìn)行求解,得到研究低滲透多孔介質(zhì)非穩(wěn)態(tài)滲流問(wèn)題需考慮動(dòng)邊界的影響,并分析了模型求解的正確性.研究結(jié)果表明,游離瓦斯產(chǎn)生孔隙壓力,以體積力作用于煤體;吸附瓦斯影響煤體的力學(xué)性質(zhì),宏觀上表現(xiàn)為彈性模量降低,煤體強(qiáng)度減弱.在高瓦斯煤層中,應(yīng)力場(chǎng)與瓦斯場(chǎng)共同作用于煤體,彈性勢(shì)能與瓦斯內(nèi)能大量積聚,瓦斯煤體變形系統(tǒng)處于非穩(wěn)定平衡狀態(tài),遇外部擾動(dòng)失穩(wěn)而發(fā)生沖擊地壓,稱(chēng)為高瓦斯煤層沖擊地壓.

    本文在以上研究的基礎(chǔ)上,以高瓦斯厚煤層中圓形斷面巷道為分析對(duì)象,研究高瓦斯煤層沖擊地壓發(fā)生條件及其影響因素,為進(jìn)一步研發(fā)高瓦斯煤層沖擊地壓預(yù)測(cè)與防治技術(shù)提供理論依據(jù).

    1 基本方程及其解析解

    根據(jù)文獻(xiàn)[31],沖擊地壓是采礦活動(dòng)形成的地下煤巖變形系統(tǒng)在外部擾動(dòng)下失穩(wěn)而發(fā)生的動(dòng)力現(xiàn)象.煤巖變形系統(tǒng)的組成材料為煤和巖石,由于煤巖材料具有應(yīng)變軟化性質(zhì),在外部載荷作用下,煤巖變形系統(tǒng)中部分區(qū)域(塑性變形區(qū))的煤巖材料變成了應(yīng)變軟化的非穩(wěn)定材料,整個(gè)煤巖變形系統(tǒng)處于非穩(wěn)定平衡狀態(tài),在外部擾動(dòng)下系統(tǒng)失穩(wěn)而發(fā)生沖擊地壓.系統(tǒng)失穩(wěn)時(shí)滿(mǎn)足的條件稱(chēng)為沖擊地壓發(fā)生條件,根據(jù)擾動(dòng)響應(yīng)判別準(zhǔn)則得到?jīng)_擊地壓發(fā)生條件.由擾動(dòng)響應(yīng)判別準(zhǔn)則,煤巖變形系統(tǒng)在外載荷P作用下產(chǎn)生的塑性變形區(qū)特征深度為ρ,對(duì)于外載荷的一個(gè)微小擾動(dòng)ΔP,塑性變形區(qū)特征深度產(chǎn)生增量Δρ,如果Δρ→∞,則系統(tǒng)將失穩(wěn)而發(fā)生沖擊地壓,即沖擊地壓發(fā)生條件為

    巷道為典型的地下開(kāi)挖形成的煤巖變形系統(tǒng).巷道斷面形狀根據(jù)煤礦井下地質(zhì)條件和生產(chǎn)的要求而選擇,如矩形、直墻拱形等.雖然一般很少采用圓形,但由于圓形斷面巷道結(jié)構(gòu)簡(jiǎn)單,便于分析,且其圍巖塑性變形區(qū)特征深度與其他斷面巷道偏差不大[32],因此本文以高瓦斯厚煤層中圓形斷面巷道為分析對(duì)象,研究高瓦斯煤層沖擊地壓發(fā)生條件及其影響因素.相應(yīng)地將滿(mǎn)足沖擊地壓發(fā)生條件的塑性區(qū)半徑定義為“臨界塑性區(qū)半徑”,將滿(mǎn)足沖擊地壓發(fā)生條件的煤體應(yīng)力定義為“臨界應(yīng)力”.臨界塑性區(qū)半徑和臨界應(yīng)力是沖擊地壓發(fā)生條件的具體體現(xiàn).

    在原巖應(yīng)力為P、瓦斯孔隙壓力為Pg的煤層中開(kāi)挖半徑為a的圓形斷面巷道.在巷道壁處作用的支護(hù)應(yīng)力為Ps、瓦斯孔隙壓力為Pa.沿巷道軸線方向取單位長(zhǎng)度計(jì)算,在不計(jì)體力的條件下,為軸對(duì)稱(chēng)平面應(yīng)變問(wèn)題,建立極坐標(biāo)系,如圖1所示.

    假設(shè)原巖應(yīng)力較大,在巷道周邊已經(jīng)出現(xiàn)外半徑為ρ的塑性變形區(qū).假設(shè)b為巷道開(kāi)挖的影響區(qū)域,r=b處的煤體應(yīng)力為P,瓦斯孔隙壓力為Pg.

    圖1 圓形斷面巷道分析模型Fig.1 Analysis model of circular section roadway

    1.1 瓦斯壓力分布規(guī)律

    由文獻(xiàn)[6]知,巷道開(kāi)挖后,巷道圍巖中瓦斯流動(dòng)為軸對(duì)稱(chēng)徑向滲流,單位時(shí)間流過(guò)單位長(zhǎng)度半徑為r的柱面的瓦斯流量Q為

    式中,p為r處的瓦斯孔隙壓力,K為滲透系數(shù).由式(2)得

    式中,c1為積分常數(shù).

    由式(3)和邊界條件p(a)=Pa,p(b)=Pg,得

    將式(4a)和式(4b)代入式(3)得瓦斯壓力分布規(guī)律為1.2煤體骨架平衡方程

    對(duì)于巖石類(lèi)材料,由修正的太沙基有效應(yīng)力原理,作用于煤體骨架上的有效應(yīng)力與孔隙壓力p的關(guān)系為

    圓形斷面巷道圍巖外載荷作用下產(chǎn)生的徑向應(yīng)力分量為σr、環(huán)向應(yīng)力分量為σθ、軸向應(yīng)力分量為,徑向有效應(yīng)力分量為、環(huán)向有效應(yīng)力分量為、軸向有效應(yīng)力分量為,則由式 (6)得

    忽略體力情況下的軸對(duì)稱(chēng)平面應(yīng)變問(wèn)題的平衡方程為

    將式(5b),式(6a),式(6b)代入式(7)得

    1.3 幾何方程與位移分布規(guī)律

    令徑向應(yīng)變分量為εr、環(huán)向應(yīng)變分量為εθ、徑向位移為μ,則幾何方程

    假設(shè)巷道圍巖不可壓縮,ua為巷道壁處的徑向位移,′為有效應(yīng)力強(qiáng)度,為應(yīng)變強(qiáng)度,則

    1.4 本構(gòu)方程與損傷變量

    假設(shè)彈性區(qū)無(wú)損傷,損傷變量D=0,則彈性區(qū)本構(gòu)方程為

    式中,ˉE稱(chēng)為等效彈性模量.

    在彈性區(qū)與塑性區(qū)交界r=ρ處,滿(mǎn)足Mohr-Coulomb屈服條件

    假設(shè)塑性區(qū)損傷線性演化,且 D(εc)= 0,(λ為單軸壓縮應(yīng)力應(yīng)變曲線峰后降模量,εu為峰后應(yīng)力降低至0時(shí)對(duì)應(yīng)的應(yīng)變),則損傷演化方程為

    由于損傷,有效應(yīng)力分量分別變?yōu)?σ′θ=則塑性區(qū)本構(gòu)方程變?yōu)?/p>

    將式(15)代入式(16),得

    1.5 彈性區(qū)應(yīng)力分布規(guī)律

    彈性區(qū)ρ≤r≤b,將式(12)代入式(8)得

    代入式(12)和式(10a),得

    1.6 塑性區(qū)應(yīng)力分布規(guī)律

    塑性區(qū)a≤r<ρ,將式(17)代入式(8)得

    代入式(17)和式(10a),得

    1.7 沖擊地壓發(fā)生條件

    由r=ρ徑向應(yīng)力連續(xù)條件,令b→∞,得煤體應(yīng)力P與塑性區(qū)半徑ρ的關(guān)系式

    由于臨界塑性區(qū)半徑和臨界應(yīng)力是沖擊地壓發(fā)生條件的具體體現(xiàn),首先推導(dǎo)臨界塑性區(qū)半徑和臨界應(yīng)力的表達(dá)式.根據(jù)沖擊地壓發(fā)生的擾動(dòng)響應(yīng)判別準(zhǔn)則[13],由式(1)和式(22)得臨界塑性區(qū)半徑

    代入式(22),得臨界應(yīng)力

    2 高瓦斯煤層沖擊地壓發(fā)生的影響因素分析

    高瓦斯煤層巷道開(kāi)挖前,瓦斯孔隙壓力處處相等,等于初始孔隙壓力Pg.巷道開(kāi)挖后在巷道壁處形成自由面,煤層瓦斯逐漸解吸,沿徑向巷道方向流動(dòng).巷道形成的時(shí)刻巷道壁處孔隙壓力等于初始孔隙壓力,即Pa0=Pg;隨著瓦斯不斷解吸,向巷道壁方向不斷流動(dòng),巷道壁周?chē)紫秹毫χ饾u降低,即Pa<Pg;直到最終時(shí)刻Pa=Pamin為止,巷道圍巖處處達(dá)到動(dòng)態(tài)平衡.

    以阜新五龍礦高瓦斯厚煤層掘進(jìn)的運(yùn)輸平巷為例,分析高瓦斯煤層沖擊地壓發(fā)生的影響因素.該巷道為矩形斷面,橫截面積為24m2.根據(jù)文獻(xiàn)[32],矩形斷面巷道圍巖塑性區(qū)深度比相同橫截面積的圓形斷面略大,當(dāng)圓形斷面巷道半徑取a=3m時(shí),兩種不同斷面巷道的塑性區(qū)深度相等.按圓形斷面計(jì)算所得結(jié)論與實(shí)際情況之間的偏差可以忽略.

    實(shí)測(cè)數(shù)據(jù)為:E=2.5GPa,λ/E=1.5,α =0.5,C=2.5MPa,φ=35?;巷道埋深H0=880m,上覆巖層平均容重27kN/m3,則煤體應(yīng)力P=23.76MPa;采用的支護(hù)方式近似為彈性支護(hù),支護(hù)應(yīng)力初始值為 Ps0=0,最終值為 Psmax=1.2MPa;初始瓦斯孔隙壓力Pg=5MPa,巷道壁處孔隙壓力最終降至Pamin=0.5MPa后保持基本穩(wěn)定,此時(shí)支護(hù)應(yīng)力也近似達(dá)到最終值.

    在巷道形成的初始時(shí)刻,Ps=Ps0=0,Pa=Pa0=Pg=5MPa.由式 (23)和式 (24)計(jì)算得到臨界塑性區(qū)半徑ρcr=1.1a=3.3m,臨界煤體應(yīng)力Pcr=1.34σc=12.8MPa.因?yàn)?P=23.76MPa> Pcr,所以該巷道存在沖擊地壓危險(xiǎn).但是,實(shí)際上巷道在掘進(jìn)初期并沒(méi)有發(fā)生沖擊地壓,原因在于在掘進(jìn)面附近存在開(kāi)挖面空間效應(yīng)[33],掘進(jìn)面前方尚未開(kāi)挖的煤體對(duì)巷道變形起到了限制作用,使得掘進(jìn)面附近不具備發(fā)生沖擊地壓條件.

    隨著掘進(jìn)面向前推進(jìn),后方一定距離范圍內(nèi)的巷道支護(hù)應(yīng)力達(dá)到最大值Ps=Psmax=1.2MPa,同時(shí)隨著瓦斯解吸滲流的進(jìn)行,巷道壁處孔隙壓力降低到最小值Pa=Pamin=0.5MPa.由式(23)和式(24)計(jì)算得到臨界塑性區(qū)半徑ρcr=1.36a=4.07m,臨界煤體應(yīng)力 Pcr=1.74σc=16.67MPa.P=23.76MPa>Pcr,所以該巷道仍然存在沖擊地壓危險(xiǎn).原因在于支護(hù)應(yīng)力較小.如果提高支護(hù)應(yīng)力,當(dāng)達(dá)到 Ps=3.83MPa,臨界應(yīng)力 Pcr=23.77MPa>P,則不會(huì)發(fā)生沖擊地壓.

    2.1 臨界塑性區(qū)半徑的影響因素分析

    由式(23)可知,臨界塑性區(qū)半徑的主要影響因素為模量比、瓦斯孔隙壓力、支護(hù)應(yīng)力和內(nèi)摩擦角.采用控制變量法,分別對(duì)各個(gè)因素對(duì)臨界塑性區(qū)半徑的影響程度進(jìn)行分析.在其他因素?cái)?shù)值不變(實(shí)測(cè)值)的條件下,由式(23)可以得到單個(gè)因素對(duì)臨界塑性區(qū)半徑的影響規(guī)律,如圖2~圖5所示.

    圖2 模量比對(duì)臨界塑性區(qū)半徑的影響Fig.2 Influenc of modulus ratio on the critical plastic zone radius

    圖3 瓦斯孔隙壓力對(duì)臨界塑性區(qū)半徑的影響Fig.3 Influenc of gas pore pressure on the critical plastic zone radius

    圖4 支護(hù)應(yīng)力對(duì)臨界塑性區(qū)半徑的影響Fig.4 E ff ect of supporting stress on the critical plastic zone radius

    圖5 內(nèi)摩擦角對(duì)臨界塑性區(qū)半徑的影響Fig.5 Influenc of angle of internal friction on the critical plastic zone radius

    由圖2可見(jiàn),臨界塑性區(qū)半徑隨模量比λ/E的增大而快速減小.當(dāng)λ/E=0時(shí),為理想彈塑性情況,ρcr/a→∞,表明理想彈塑性煤巖應(yīng)力應(yīng)變曲線峰后沒(méi)有出現(xiàn)軟化現(xiàn)象,不可能發(fā)生沖擊地壓;當(dāng)λ/E→∞時(shí),為脆性煤巖情況,ρcr/a→1,表明脆性煤巖應(yīng)力應(yīng)變曲線峰后快速跌落,巷道壁處剛進(jìn)入塑性變形即發(fā)生沖擊地壓;當(dāng)0<λ/E<∞時(shí),為具有塑性軟化的情況,∞>ρcr/a>1,煤巖越脆硬沖擊傾向性越大,越容易發(fā)生沖擊地壓.采用冪函數(shù)進(jìn)行曲線擬合,得

    由圖3可見(jiàn),臨界塑性區(qū)半徑隨瓦斯孔隙壓力的增大而減小,但降低的幅度不大,孔隙壓力越大越容易發(fā)生沖擊地壓.采用線性函數(shù)進(jìn)行曲線擬合,得=-0.0319Pa+1.5444.一般情況下pa<5,則第2項(xiàng)比第1項(xiàng)大兩個(gè)數(shù)量級(jí),因此可以近似認(rèn)為瓦斯孔隙壓力對(duì)臨界塑性區(qū)半徑幾乎沒(méi)有影響,可以忽略,則=常數(shù),能夠滿(mǎn)足現(xiàn)場(chǎng)實(shí)際要求.

    由圖4可見(jiàn),臨界塑性區(qū)半徑隨支護(hù)應(yīng)力的增大而增大,支護(hù)應(yīng)力越大越不容易發(fā)生沖擊地壓,加強(qiáng)支護(hù)可有效降低沖擊地壓的發(fā)生.采用線性函數(shù)進(jìn)行曲線擬合,得 =0.0642Ps+1.2959.一般情況下P<2,則第2項(xiàng)比第1項(xiàng)大一個(gè)數(shù)量級(jí),因

    s此支護(hù)應(yīng)力對(duì)臨界塑性區(qū)半徑的影響不可忽略.提高支護(hù)應(yīng)力可以使臨界塑性區(qū)半徑增大,有效降低沖擊地壓發(fā)生的可能性,所以加強(qiáng)支護(hù)是防治沖擊地壓的一項(xiàng)重要措施.

    由圖5可見(jiàn),臨界塑性區(qū)半徑隨內(nèi)摩擦角的增大而增大.表明內(nèi)摩擦角越大越不容易發(fā)生沖擊地壓.當(dāng) φ → 0時(shí)采用二次函數(shù)進(jìn)行曲線擬合,得0.0001φ+1.5444.當(dāng) φ的單位取為弧度時(shí),=0.164φ2+0.0057φ +1.2981.

    2.2 臨界應(yīng)力的影響因素分析

    由式(24)可知,臨界應(yīng)力的主要影響因素為模量比、瓦斯孔隙壓力、支護(hù)應(yīng)力和內(nèi)摩擦角.采用控制變量法,分別對(duì)各個(gè)因素對(duì)臨界應(yīng)力的影響程度進(jìn)行分析.在其他因素?cái)?shù)值不變(實(shí)測(cè)值)的條件下,由式(24)可以得到單個(gè)因素對(duì)臨界應(yīng)力的影響規(guī)律,如圖6~圖9所示.

    圖6 模量比對(duì)臨界應(yīng)力的影響Fig.6 Influenc of modulus ratio on critical stress

    圖7 瓦斯孔隙壓力對(duì)臨界應(yīng)力的影響Fig.7 Influenc of gas pore pressure on critical stress

    圖8 支護(hù)應(yīng)力對(duì)臨界應(yīng)力的影響Fig.8 E ff ect of supporting stress on critical stress

    圖9 內(nèi)摩擦角對(duì)臨界應(yīng)力的影響Fig.9 Influenc of angle of internal friction on critical stress

    由圖6可見(jiàn),臨界應(yīng)力隨模量比λ/E的增大而快速減小.當(dāng)λ/E=0時(shí),Pcr/σc→∞,表明理想彈塑性煤巖應(yīng)力應(yīng)變曲線沒(méi)有峰后軟化現(xiàn)象,不可能發(fā)生沖擊地壓;當(dāng)λ/E→∞時(shí)表明脆性煤巖應(yīng)力應(yīng)變曲線峰后快速跌落,巷道壁處剛進(jìn)入塑性變形即發(fā)生沖擊地壓;當(dāng)0<λ/E<∞時(shí)表明煤巖越脆硬沖擊傾向性越大,越容易發(fā)生沖擊地壓.采用冪函數(shù)進(jìn)行曲線擬合,得

    由圖7可見(jiàn),臨界應(yīng)力隨瓦斯孔隙壓力的增大而減小,但降低的幅度不大.表明孔隙壓力越大越容易發(fā)生沖擊地壓.采用線性函數(shù)進(jìn)行曲線擬合,得一般情況下 Pa<5,則第2項(xiàng)比第1項(xiàng)大約一個(gè)數(shù)量級(jí),因此可以近似認(rèn)為瓦斯孔隙壓力對(duì)臨界載荷影響較小,但不可忽略.

    由圖8可見(jiàn),臨界應(yīng)力隨支護(hù)應(yīng)力的增大而增大.表明支護(hù)應(yīng)力越大越不容易發(fā)生沖擊地壓,加強(qiáng)支護(hù)可有效降低沖擊地壓的發(fā)生.采用線性函數(shù)進(jìn)行曲線擬合,得=0.2894Ps+1.4431.一般情況下Ps<2,因此支護(hù)應(yīng)力對(duì)臨界載荷的影響不可忽略.

    當(dāng)P=Pcr時(shí),由式(24)得

    加強(qiáng)支護(hù)提高支護(hù)應(yīng)力,使Ps≥Psmax時(shí),不會(huì)發(fā)生沖擊地壓.

    =0得到.按本文數(shù)據(jù),φm≈ 35?.當(dāng) φ < φm時(shí),臨界應(yīng)力隨內(nèi)摩擦角增大而減?。划?dāng)φ>φm時(shí),臨界應(yīng)力隨內(nèi)摩擦角增大而增大.采用二次函數(shù)進(jìn)行曲線擬合,得=0.0418φ2-2.7418φ+57.551.當(dāng) φ的單位取為弧度時(shí)=137.22φ2-157.09φ+57.551.

    3 結(jié)論

    (1)煤的模量比、煤層瓦斯孔隙壓力、支護(hù)應(yīng)力和內(nèi)摩擦角等對(duì)高瓦斯煤層沖擊地壓是否發(fā)生具有重要影響.其中煤的模量比、支護(hù)應(yīng)力和內(nèi)摩擦角影響較大.煤層瓦斯孔隙壓力雖然影響較小,但不可忽略.

    (2)高瓦斯煤層掘進(jìn)巷道,在掘進(jìn)面附近由于存在開(kāi)挖面空間效應(yīng),掘進(jìn)面前方尚未開(kāi)挖的煤體對(duì)巷道變形起到了限制作用,使得掘進(jìn)面附近不具備發(fā)生沖擊地壓條件.隨著掘進(jìn)面向前推進(jìn),后方一定距離范圍內(nèi)的巷道支護(hù)應(yīng)力增大,同時(shí)隨著瓦斯解吸滲流的進(jìn)行,巷道壁處孔隙壓力降低,巷道存在沖擊地壓危險(xiǎn).如果提高支護(hù)應(yīng)力,則沖擊危險(xiǎn)性降低.

    (3)高瓦斯煤層巷道發(fā)生沖擊地壓的臨界塑性區(qū)半徑、臨界應(yīng)力隨模量比λ/E的增大而快速減小,表明煤巖越脆硬沖擊傾向性越大,越容易發(fā)生沖擊地壓.

    (4)臨界塑性區(qū)半徑、臨界應(yīng)力隨瓦斯孔隙壓力的增大而減小,但降低的幅度不大,孔隙壓力越大越容易發(fā)生沖擊地壓.

    (5)臨界塑性區(qū)半徑、臨界應(yīng)力隨支護(hù)應(yīng)力的增大而增大,支護(hù)應(yīng)力越大越不容易發(fā)生沖擊地壓,加強(qiáng)支護(hù)可有效降低沖擊地壓的發(fā)生.

    (6)臨界塑性區(qū)半徑隨內(nèi)摩擦角的增大而增大,內(nèi)摩擦角越大越不容易發(fā)生沖擊地壓.臨界應(yīng)力與內(nèi)摩擦角不是單調(diào)函數(shù)關(guān)系,存在一個(gè)極小值點(diǎn).當(dāng)內(nèi)摩擦角小于此極小值時(shí),臨界應(yīng)力隨內(nèi)摩擦角增大而減小;當(dāng)內(nèi)摩擦角大于此極小值時(shí),臨界應(yīng)力隨內(nèi)摩擦角增大而增大.

    1潘一山,李忠華,章夢(mèng)濤.我國(guó)沖擊地壓分布、類(lèi)型、機(jī)理及防治研究.巖石力學(xué)與工程學(xué)報(bào),2003,22(11):1844-1851(Pan Yishan,Li Zhonghua,Zhang Mengtao.Distribution,type,mechanism and prevention of rockburst in China.Chinese Journal of Rock Mechanics and Engineering,2003,22(11):1844-1851(in Chinese))

    2郭德勇,韓德馨.煤與瓦斯突出黏滑機(jī)理研究.煤炭學(xué)報(bào),2003,28(6):598-602(Guo Deyong,Han Dexin.The stick-slip mechanism ofcoalandgasoutburst.JournalofChinaCoalSociety,2003,28(6):598-602(in Chinese))

    3李新元.“圍巖--煤體”系統(tǒng)失穩(wěn)破壞及沖擊地壓預(yù)測(cè)的探討.中國(guó)礦業(yè)大學(xué)學(xué)報(bào),2000,29(6):633-636(Li Xinyuan.Discussion on unstability and damage of surrounding rock-coal mass system and prediction of rockburst.Journal of China University of Mining&Technology,2000,29(6):633-636(in Chinese))

    4竇林名,何學(xué)秋.煤礦沖擊礦壓的分級(jí)預(yù)測(cè)研究.中國(guó)礦業(yè)大學(xué)學(xué)報(bào).2007,36(6):717-722(Dou Linming,He Xueqiu.Technique of classificatio forecasting rock burst in coal mines.Journal of China University of Mining&Technology.2007,36(6):717-722(in Chinese))

    5潘一山,李忠華,唐鑫.阜新礦區(qū)深部高瓦斯礦井沖擊地壓研究.巖石力學(xué)與工程學(xué)報(bào),2005,24(s1):5202-5205(Pan Yishan,Li Zhonghua,Tang Xin.Research on rockburst for highly gassy deep mine in Fuxin Coal mining area.Journal of Rock Mechanics and Engineering,2005,24(s1):5202-5205(in Chinese))

    6李忠華.高瓦斯煤層沖擊地壓發(fā)生理論研究及應(yīng)用.[博士論文].阜新:遼寧工程技術(shù)大學(xué)力學(xué)與工程學(xué)院,2007(Li Zhonghua.Study on theory and application of high gassy seam rockburst.[PhD Thesis].Fuxin:Liaoning Technical University Liaoning Technical University,2007(in Chinese))

    7潘一山.煤與瓦斯突出、沖擊地壓復(fù)合動(dòng)力災(zāi)害一體化研究.煤炭學(xué)報(bào),2016,41(1):105-112(Pan Yishan.Integrated study on compound dynamic disaster of coal-gas outburst and rockburst.Journal of China Coal Society,2011,41(1):105-112(in Chinese))

    8 Pan YS,Li ZH.Analysis of rock structure stability in coal mines.International Journal for Numerical and Analytical Methods in Geomechanics,2010,29(10):1045-1063

    9李鐵,蔡美峰,王金安等.深部開(kāi)采沖擊地壓與瓦斯的相關(guān)性探討.煤炭學(xué)報(bào),2005,30(5):562-567(Li Tie,Cai Meifeng,Wang Jin’an,et al.Discussion on relativity between rockburst and gas in deep exploitation.Journal of China Coal Society,2005,30(5):562-567(in Chinese))

    10李鐵,梅婷婷,李國(guó)旗等.“三軟”煤層沖擊地壓誘導(dǎo)煤與瓦斯突出力學(xué)機(jī)制研究.巖石力學(xué)與工程學(xué)報(bào),2011,30(6):1283-1288(Li Tie,Mei Tingting,Li Guoqi,et al.Mechanism study of coal and gas outburst induced by rockburst in “three-soft”coal seam.Chinese Journal of Rock Mechanics and Engineering,2011,30(6):1283-1288(in Chinese))

    11王振,尹光志,胡千庭等.高瓦斯煤層沖擊地壓與突出的誘發(fā)轉(zhuǎn)化條件研究.采礦與安全工程學(xué)報(bào),2010,27(4):572-575(Wang Zhen,Yin Guangzhi,Hu Qianting,et al.Inducing and transforming conditions from rockburst to coal gas outburst in a high gassy coal seam.Journal of Mining&Safety Engineering,2010,27(4):572-575(in Chinese))

    12馬海峰.高瓦斯煤層采動(dòng)應(yīng)力對(duì)采場(chǎng)瓦斯壓力的影響研究.[碩士論文].淮南:安徽理工大學(xué)能源安全學(xué)院,2011(Ma Haifeng.The research on the e ff ect of mining stress on gas pressure in stop in high gassy seam.[Master Thesis].Huainan:School of Energy and Safety Anhui University of Science and Technology,2011(in Chinese))

    13董飛亞.深部高瓦斯煤層沖擊地壓預(yù)測(cè)防治技術(shù)研究.[碩士論文].焦作:河南理工大學(xué)能源科學(xué)與工程學(xué)院,2012(Dong Feiya.Research on rockburst prediction and prevention and control technology in deep coal seam with high gas.[Master Thesis].Jiaozuo:Energy Science and Engineering Henan Polytechnic University,2012(in Chinese))

    14王磊.應(yīng)力場(chǎng)和瓦斯場(chǎng)采動(dòng)耦合效應(yīng)研究.[博士論文].淮南:安徽理工大學(xué)能源與安全學(xué)院,2010(Wang Lei.Research of caving coupling e ff ects of stress fiel and gas field[PhD Thesis].Huainan:School of Energy and Safety Anhui University of Science and Technology,2010(in Chinese))

    15徐晨陽(yáng).高瓦斯煤層沖擊地壓特征機(jī)理研究現(xiàn)狀.山西大同大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,30(1):69-72(Xu Chenyang.The mechanism of rock bursts characteristics of gas coal seam.Journal of Shanxi Datong University(JCR Science Edition),2014,30(1):69-72(in Chinese))

    16王濤,王曌華,劉華博等.沖擊地壓后瓦斯異常涌出條件及致災(zāi)原因分析,煤炭學(xué)報(bào),2014,39(2):371-376(Wang Tao,Wang Zhaohua,Liu Huabo,et al.Discussion about the mechanism of gas disaster induced by coal bump.Journal of China Coal Society,2014,39(2):371-376(in Chinese))

    17梁冰,章夢(mèng)濤,潘一山.煤和瓦斯突出的固流藕合失穩(wěn)理論.煤炭學(xué)報(bào),1995,20(5):492-496(Liang Bing,Zhang Mengtao,Pan Yishan.et al.Theory of instability of fl w fixatio coupling for coal and gas outburst.Journal of China Coal Society,1995,20(5):492-496(in Chinese))

    18梁冰,王泳嘉,章夢(mèng)濤.含瓦斯煤的內(nèi)時(shí)本構(gòu)關(guān)系及其參數(shù)的實(shí)驗(yàn)研究.固體力學(xué)學(xué)報(bào),1996,17(3):229-234(Liang Bing,Wang Yongjia,Zhang Mengtao.The study of endochronic constitutive equations of coal e ff ected by gas and the determining of parameters by experiment.Chinese Journal of Solid Mechanics,1996,17(3):229-234(in Chinese))

    19俞善炳.煤與瓦斯突出的一維流動(dòng)模型和啟動(dòng)判據(jù).力學(xué)學(xué)報(bào),1992,24(4):418-431(Yu Shanbing.One-dimensional fl w model for coal-gas out bursts and initiation criterion.Acta Mechanica Sinica,1992,24(4):418-431(in Chinese))

    20俞善炳.恒穩(wěn)推進(jìn)的煤與瓦斯突出.力學(xué)學(xué)報(bào),1988,20(2):97-106(Yu Shanbing.Steady advance of coal and gas bursts.Acta Mechanica Sinica,1988,20(2):97-106(in Chinese))

    21趙陽(yáng)升.瓦斯壓力在突出中作用的數(shù)值模擬研究.巖石力學(xué)與工程學(xué)報(bào),1993,12(4):328-337(Zhao Yangsheng.The numerical simulation on e ff ect of methane pressure to outburst.Chinese Journal of Rock Mechanics and Engineering,1993,12(4):328-337(in Chinese))

    22藍(lán)航,齊慶新,潘俊鋒等.我國(guó)煤礦沖擊地壓特點(diǎn)及防治技術(shù)分析.煤炭科學(xué)技術(shù),2011,39(1):11-36(Lan Hang,Qi Qingxin,Pan Junfeng,et al.Analysis on features as well as prevention and control technology of mine strata pressure bumping in China.Coal Science and Technology,2011,39(1):11-36(in Chinese))

    23崔乃鑫,李忠華,潘一山.考慮瓦斯影響的煤層沖擊地壓鉆屑量指標(biāo)研究.遼寧工程技術(shù)大學(xué)學(xué)報(bào),2006,25(2):192-193(Cui Naixin,Li Zhonghua,Pan Yishan.Study on index of drilling bits for coalbed rockburst influence by gas.Journal of Liaoning Technical University,2006,25(2):192-193(in Chinese))

    24彭蘇萍,高云峰,楊瑞召等.AVO探測(cè)煤層瓦斯富集的理論探討和初步實(shí)踐——以淮南煤田為例.地球物理學(xué)報(bào),2005,48(6):1475-1486(Peng Suping,Gao Yunfeng,Yang Ruizhao,et al.Theory and application of AVO for detection of coalbed methane——A case from the Huainan coalfield Chinese Journal of Geophysics,2005,48(6):1475-1486(in Chinese))

    25梁盛開(kāi),羅楊陽(yáng),曹瓊.煤礦瓦斯突出的多因素綜合預(yù)測(cè).固體力學(xué)學(xué)報(bào),2010,31(S1):177-179(Liang Shengkai,Luo Yangyang,Cao Qiong.Prediction of coal mine’s coal and gas outburst based on multi-factors probability.Chinese Journal of Solid Mechanics,2010,31(S1):177-179(in Chinese))

    26宋真龍,韓佩博,李文璞等.瓦斯對(duì)沖擊性煤樣能量耗散的影響.煤炭學(xué)報(bào),2015,40(4):843-849(Song Zhenlong,Han Peibo,Li Wenpu,et al.Impact of energy dissipation of coal samples with rockburst tendency from gas in its failure process.Journal of China Coal Society,2015,40(4):843-849(in Chinese))

    27劉??h,鮮學(xué)福,姜德義.煤與瓦斯延期突出機(jī)理及其預(yù)測(cè)預(yù)報(bào)的研究.巖石力學(xué)與工程學(xué)報(bào),2001,21(5):647-650(Liu Bao-xian,Xian Xuefu,Jiang Deyi.Research on the mechanism and prediction of delay outburst of coal and gas.Chinese Journal of Rock Mechanics and Engineering,2001,21(5):647-650(in Chinese))

    28張志剛,程波.含瓦斯煤體非線性滲流模型.中國(guó)礦業(yè)大學(xué)學(xué)報(bào),2015,44(3):453-459(Zhang Zhigang,Cheng Bo.Study of a nonlinear seepage model of coal containing gas.Journal of China University of Mining&Technology,2015,44(3):453-459(in Chinese))

    29郭為,胡志明,左羅等.頁(yè)巖基質(zhì)解吸--擴(kuò)散--滲流耦合實(shí)驗(yàn)及數(shù)學(xué)模型.力學(xué)學(xué)報(bào),2015,47(6):916-922(Guo Wei,Hu Zhiming,Zuo Luo,et al.Gas desorption--di ff usion--seepage coupled experiment of shale matrix and mathematic model.Chinese Journal of Theoretical and Applied Mechanics,2015,47(6):916-922(in Chi-nese))

    30劉文超,姚軍,陳掌星等.低滲透多孔介質(zhì)滲流動(dòng)邊界模型的解析與數(shù)值解.力學(xué)學(xué)報(bào),2015,47(6):605-612(Liu Wenchao,Yao Jun,Chen Zhangxin,et al.Research on analytical and numerical solutions of a moving boundary model of seepage fl w in low-permeable porous media.Chinese Journal of Theoretical and Applied Mechanics,2015,47(6):605-612(in Chinese))

    31潘一山.沖擊地壓發(fā)生和破壞過(guò)程研究.[博士論文].北京:清華大學(xué),1999(Pan Yishan.Study on rockburst Initiation and failure propagation.[PhD Thesis].Beijing:Tsinghua University,1999(in Chinese))

    32蔡光華,陸海軍,陳寶銀等.矩形和直墻拱形斷面圍巖巷道破壞的數(shù)值模擬研究.武漢工業(yè)學(xué)院學(xué)報(bào),2011,30(1):74-78(Cai Guanghua,Lu Haijun,Chen Baoyin,et al.Research on numerical simulation of surrounding rock laneway failure under rectangular and straight wall arch.Journal of Wuhan Polytechnic University,2011,30(1):74-78(in Chinese))

    33侯公羽.基于開(kāi)挖面空間效應(yīng)的圍巖--支護(hù)相互作用機(jī)制.巖石力學(xué)與工程學(xué)報(bào),2011,30(1):2871-2877(Hou Gongyu.Interaction mechanism between surrounding rock and support based on spatial e ff ect of excavation face.Chinese Journal of Rock Mechanics and Engineering,2011,30(1):2871-2877(in Chinese))

    CONDITION AND INFLUENCE FACTOR OF ROCK BURST IN HIGH GASSY COAL SEAM1)

    Yin Wanlei?,2)Pan Yishan?,?,3)Li Zhonghua??(Research Institute of Rockburst,School of Mechanics and Engineering,Liaoning Technical University,F(xiàn)uxin 123000,Liaoning,China)?(School of Physics,Liaoning University,Shenyang 110036,China)

    Aiming at the problem of rock burst in high gassy coal seam,the occurrence conditions of rock burst are obtained by analytical analysis,the influenc rule of main factors on the radius and the critical stress of the critical plastic zone are analyzed.In connection with Wulong mining practice,the e ff ects of the coal modulus ratio,gas pore pressure,support stress,and internal friction angle on rock burst of in high gassy coal seam are analyzed in comparison.The results show that,non excavated solid coal plays a limiting role in the deformation of the roadway because of the spatial e ff ect with excavation face nearby roadway heading face about high gassy coal seam.It reduces the danger of rock burst,the support stress of the roadway is increased in a certain distance from the rear along with the excavating face ahead.At the same time,along with the gas desorption seepage,the pore pressure of the tunnel wall is decreased,and the risk of rockburst is increased,at this point to improve the support stress,the impact of risk is reduced.The critical plastic zone radius and the critical stress of the high gassy coal seam tunnel decrease rapidly with the increase of the modulus ratio λ/E and the pore pressure,and they increase with the increase of the support stress,the radius of the critical plastic zone increases with the increase of the internal friction angle,the relationship between the critical stress and the internal friction angle is not monotone,there exists a minimum value,when the internal friction angle is less than the minimum value,the critical stress decreases with the increase of the internal friction angle;when the internal friction angle is larger,the results are opposite.

    rock burst,high gassy coal seam,modulus ratio,pore pressure,support stress,internal friction angle

    TD324

    :A

    10.6052/0459-1879-16-302

    2016–10–28 收稿,2017–03–12 錄用,2017–03–16 網(wǎng)絡(luò)版發(fā)表.

    1)國(guó)家重點(diǎn)研發(fā)計(jì)劃(2016YFC0801401),國(guó)家自然科學(xué)基金(51174107,51374119),國(guó)家青年自然科學(xué)基金(51404131)資助項(xiàng)目.

    2)尹萬(wàn)蕾,博士研究生,主要從事沖擊地壓預(yù)測(cè)防治研究.E-mail:ywl696@126.com

    3)潘一山,教授,主要從事固體力學(xué)和巖石力學(xué)方面的教學(xué)與研究工作.E-mail:panyish cn@sina.com

    尹萬(wàn)蕾,潘一山,李忠華.高瓦斯煤層沖擊地壓發(fā)生條件與影響因素.力學(xué)學(xué)報(bào),2017,49(3):716-725

    Yin Wanlei,Pan Yishan,Li Zhonghua.Condition and influenc factor of rock burst in high gassy coal seam.Chinese Journal of Theoretical and Applied Mechanics,2017,49(3):716-725

    猜你喜歡
    摩擦角煤體塑性
    應(yīng)用摩擦角,巧解動(dòng)力學(xué)問(wèn)題
    基于應(yīng)變梯度的微尺度金屬塑性行為研究
    硬脆材料的塑性域加工
    注熱井周?chē)后w蠕變過(guò)程的滲透率變化規(guī)律模擬研究
    鈹材料塑性域加工可行性研究
    借助摩擦角 快解勻速運(yùn)動(dòng)問(wèn)題
    正交試驗(yàn)下煤體滲透性影響因素評(píng)價(jià)
    煤炭工程(2019年4期)2019-05-05 02:30:16
    以“中央廚房”為突破口探索時(shí)政報(bào)道的融煤體之路——以浙江之聲為例
    石英玻璃的熱輔助高效塑性域干磨削
    摩擦角在平衡問(wèn)題中的應(yīng)用
    男人爽女人下面视频在线观看| 久久国产亚洲av麻豆专区| 两性夫妻黄色片| 成人手机av| 一区在线观看完整版| 在线观看免费午夜福利视频| 日本av免费视频播放| 成年人黄色毛片网站| 亚洲熟女精品中文字幕| 这个男人来自地球电影免费观看| 满18在线观看网站| 狠狠狠狠99中文字幕| 亚洲av欧美aⅴ国产| 精品国产乱码久久久久久小说| 汤姆久久久久久久影院中文字幕| 国产高清国产精品国产三级| 国产日韩欧美视频二区| 成年动漫av网址| 性色av乱码一区二区三区2| 亚洲熟女精品中文字幕| 99国产精品99久久久久| 在线亚洲精品国产二区图片欧美| 久久女婷五月综合色啪小说| bbb黄色大片| 我的亚洲天堂| 18禁观看日本| 老司机影院毛片| 视频区图区小说| 国产99久久九九免费精品| 热re99久久精品国产66热6| 国产欧美日韩一区二区三 | 国产成人av教育| 一级毛片女人18水好多| 国产欧美日韩一区二区三区在线| 一个人免费在线观看的高清视频 | 黄色视频在线播放观看不卡| 18禁国产床啪视频网站| 成人黄色视频免费在线看| 成人18禁高潮啪啪吃奶动态图| 国产精品1区2区在线观看. | 亚洲成人免费电影在线观看| 亚洲国产精品999| 日韩大码丰满熟妇| 日本av手机在线免费观看| 亚洲中文字幕日韩| 亚洲国产日韩一区二区| 热99re8久久精品国产| 一级毛片女人18水好多| 亚洲 国产 在线| 久久久水蜜桃国产精品网| 香蕉国产在线看| 淫妇啪啪啪对白视频 | 亚洲av欧美aⅴ国产| 欧美一级毛片孕妇| 亚洲欧美精品综合一区二区三区| 久久久欧美国产精品| 久久免费观看电影| 制服人妻中文乱码| 蜜桃国产av成人99| 国产精品一区二区在线观看99| 99热国产这里只有精品6| 窝窝影院91人妻| 12—13女人毛片做爰片一| 一二三四社区在线视频社区8| 久久九九热精品免费| 久久综合国产亚洲精品| 亚洲 国产 在线| 91成年电影在线观看| 国产麻豆69| 色播在线永久视频| 亚洲第一av免费看| 国产福利在线免费观看视频| 中文字幕最新亚洲高清| 国产亚洲一区二区精品| 91成年电影在线观看| 日本vs欧美在线观看视频| 久久精品国产综合久久久| 精品亚洲成a人片在线观看| 99精国产麻豆久久婷婷| 国产亚洲午夜精品一区二区久久| 久久久久久免费高清国产稀缺| 亚洲男人天堂网一区| 国产一区二区三区av在线| 精品一区二区三区av网在线观看 | 免费黄频网站在线观看国产| 纵有疾风起免费观看全集完整版| 欧美精品av麻豆av| 国产成人欧美| 90打野战视频偷拍视频| 日本猛色少妇xxxxx猛交久久| 久久久久国产精品人妻一区二区| 18禁国产床啪视频网站| 99热国产这里只有精品6| 99久久综合免费| 好男人电影高清在线观看| 日韩三级视频一区二区三区| 成人国语在线视频| 丝袜脚勾引网站| 久热这里只有精品99| 亚洲激情五月婷婷啪啪| 国产免费视频播放在线视频| 午夜影院在线不卡| 色精品久久人妻99蜜桃| 日日摸夜夜添夜夜添小说| 亚洲专区中文字幕在线| 亚洲欧美一区二区三区黑人| 丰满迷人的少妇在线观看| 蜜桃国产av成人99| 久久人妻熟女aⅴ| 2018国产大陆天天弄谢| 日韩三级视频一区二区三区| 可以免费在线观看a视频的电影网站| 亚洲午夜精品一区,二区,三区| 国产av一区二区精品久久| 免费人妻精品一区二区三区视频| 亚洲av片天天在线观看| 午夜福利乱码中文字幕| 美女主播在线视频| 老司机亚洲免费影院| 日韩欧美一区视频在线观看| 久久香蕉激情| 老司机影院成人| 国产欧美日韩一区二区三 | 狂野欧美激情性bbbbbb| 91精品三级在线观看| 久久久精品国产亚洲av高清涩受| 手机成人av网站| 18禁国产床啪视频网站| 日韩大片免费观看网站| 国产日韩欧美视频二区| 免费高清在线观看日韩| 高清在线国产一区| 在线亚洲精品国产二区图片欧美| 99精国产麻豆久久婷婷| 国产一区二区 视频在线| 成人av一区二区三区在线看 | 老司机在亚洲福利影院| 国产一区二区三区综合在线观看| 狂野欧美激情性xxxx| 日韩 亚洲 欧美在线| 99精品欧美一区二区三区四区| 日韩有码中文字幕| 亚洲午夜精品一区,二区,三区| 亚洲精品一二三| 精品熟女少妇八av免费久了| 国产av又大| 国产成人免费无遮挡视频| 国产日韩欧美亚洲二区| 欧美日韩成人在线一区二区| 日日摸夜夜添夜夜添小说| av在线老鸭窝| 一区在线观看完整版| 女人爽到高潮嗷嗷叫在线视频| 免费少妇av软件| 国产成人影院久久av| 丰满迷人的少妇在线观看| 自线自在国产av| 免费观看人在逋| 法律面前人人平等表现在哪些方面 | 欧美+亚洲+日韩+国产| 中文字幕人妻丝袜制服| xxxhd国产人妻xxx| 精品国产一区二区三区四区第35| 亚洲人成电影观看| 国产成人精品久久二区二区免费| 久久久欧美国产精品| 19禁男女啪啪无遮挡网站| 午夜福利在线观看吧| 国产成人av激情在线播放| 777久久人妻少妇嫩草av网站| 搡老岳熟女国产| 欧美精品人与动牲交sv欧美| 欧美少妇被猛烈插入视频| 午夜两性在线视频| 成在线人永久免费视频| svipshipincom国产片| 老汉色∧v一级毛片| 国产有黄有色有爽视频| 欧美少妇被猛烈插入视频| 三上悠亚av全集在线观看| 欧美激情高清一区二区三区| 91大片在线观看| 成人国产一区最新在线观看| 交换朋友夫妻互换小说| 国产成人精品久久二区二区免费| 亚洲第一青青草原| 爱豆传媒免费全集在线观看| 悠悠久久av| 免费女性裸体啪啪无遮挡网站| 中文字幕最新亚洲高清| av天堂在线播放| 亚洲精华国产精华精| 久久99热这里只频精品6学生| 午夜福利视频在线观看免费| 成年人黄色毛片网站| 国产成+人综合+亚洲专区| 欧美午夜高清在线| 成人国产av品久久久| 久久久久久人人人人人| 精品福利永久在线观看| 看免费av毛片| 国产免费av片在线观看野外av| 国产极品粉嫩免费观看在线| 亚洲熟女毛片儿| 亚洲av电影在线进入| 熟女少妇亚洲综合色aaa.| 999精品在线视频| 亚洲七黄色美女视频| 久久久久久久久免费视频了| 男女无遮挡免费网站观看| 麻豆国产av国片精品| 亚洲va日本ⅴa欧美va伊人久久 | 淫妇啪啪啪对白视频 | 久久久精品国产亚洲av高清涩受| 国产精品99久久99久久久不卡| 亚洲中文av在线| 777米奇影视久久| 999久久久精品免费观看国产| 免费观看人在逋| 淫妇啪啪啪对白视频 | 日本撒尿小便嘘嘘汇集6| 狂野欧美激情性bbbbbb| 久久精品国产亚洲av高清一级| 日本av手机在线免费观看| 欧美97在线视频| 99精品久久久久人妻精品| 色老头精品视频在线观看| 国产欧美亚洲国产| 久久久久精品国产欧美久久久 | 香蕉国产在线看| 老司机深夜福利视频在线观看 | 一区二区三区乱码不卡18| 宅男免费午夜| 老司机午夜福利在线观看视频 | 亚洲第一av免费看| 日韩制服骚丝袜av| 999精品在线视频| 国产成人精品无人区| 黄色a级毛片大全视频| 欧美97在线视频| 国产极品粉嫩免费观看在线| 国产无遮挡羞羞视频在线观看| 亚洲欧美精品自产自拍| 国产免费一区二区三区四区乱码| 欧美日本中文国产一区发布| 伦理电影免费视频| 日韩电影二区| e午夜精品久久久久久久| 黑人巨大精品欧美一区二区mp4| 国产一区二区三区在线臀色熟女 | 俄罗斯特黄特色一大片| 啦啦啦在线免费观看视频4| 一区福利在线观看| av又黄又爽大尺度在线免费看| 另类精品久久| 1024视频免费在线观看| 一本综合久久免费| 久久精品人人爽人人爽视色| 精品少妇久久久久久888优播| 亚洲国产看品久久| 成人国语在线视频| 乱人伦中国视频| 91国产中文字幕| 国产精品熟女久久久久浪| 捣出白浆h1v1| 国产色视频综合| √禁漫天堂资源中文www| 午夜激情久久久久久久| 王馨瑶露胸无遮挡在线观看| 亚洲精品一区蜜桃| 国产成人精品在线电影| 国产成人一区二区三区免费视频网站| 久久午夜综合久久蜜桃| 精品卡一卡二卡四卡免费| 精品人妻一区二区三区麻豆| 视频区欧美日本亚洲| 免费不卡黄色视频| 国产亚洲精品一区二区www | 91老司机精品| 亚洲国产成人一精品久久久| 99国产综合亚洲精品| 中文精品一卡2卡3卡4更新| 亚洲欧美精品综合一区二区三区| 久久中文看片网| 男女床上黄色一级片免费看| 后天国语完整版免费观看| 国产成人a∨麻豆精品| 欧美精品高潮呻吟av久久| 亚洲黑人精品在线| 国产99久久九九免费精品| 两个人免费观看高清视频| 国产97色在线日韩免费| 精品人妻1区二区| 中文字幕人妻熟女乱码| 亚洲精品中文字幕在线视频| 1024视频免费在线观看| 成年人午夜在线观看视频| 国产一区二区三区在线臀色熟女 | 欧美日韩亚洲综合一区二区三区_| 午夜免费鲁丝| 亚洲精品一二三| 日本猛色少妇xxxxx猛交久久| 亚洲成人免费av在线播放| 在线观看免费午夜福利视频| 国产成人欧美| 69精品国产乱码久久久| 国产av国产精品国产| 在线观看免费视频网站a站| 成人国语在线视频| tube8黄色片| 日本一区二区免费在线视频| 中文字幕制服av| 一区二区三区精品91| 9热在线视频观看99| 久久人妻福利社区极品人妻图片| 黑人操中国人逼视频| 在线天堂中文资源库| 久久影院123| 亚洲精品久久成人aⅴ小说| 欧美黄色淫秽网站| 丝瓜视频免费看黄片| 一区在线观看完整版| 韩国高清视频一区二区三区| 久久精品亚洲av国产电影网| 一级黄色大片毛片| 亚洲精华国产精华精| 一边摸一边抽搐一进一出视频| 少妇裸体淫交视频免费看高清 | 亚洲精品粉嫩美女一区| 亚洲精品av麻豆狂野| 亚洲国产精品999| 亚洲欧洲精品一区二区精品久久久| 久久精品熟女亚洲av麻豆精品| 爱豆传媒免费全集在线观看| 新久久久久国产一级毛片| 国产成人一区二区三区免费视频网站| 国产深夜福利视频在线观看| 极品人妻少妇av视频| 老熟妇仑乱视频hdxx| 国产一区二区三区av在线| 午夜两性在线视频| 亚洲美女黄色视频免费看| 国产日韩欧美视频二区| 成人手机av| 国产亚洲一区二区精品| 亚洲色图 男人天堂 中文字幕| 国产视频一区二区在线看| 电影成人av| 午夜激情av网站| 欧美精品人与动牲交sv欧美| 久久国产精品影院| 丝袜人妻中文字幕| 日韩中文字幕欧美一区二区| 久久精品亚洲av国产电影网| 老司机靠b影院| tocl精华| 欧美日韩黄片免| 在线观看免费日韩欧美大片| 欧美亚洲日本最大视频资源| 欧美日韩一级在线毛片| 午夜视频精品福利| 国产麻豆69| 老熟女久久久| 99久久精品国产亚洲精品| 国产在线一区二区三区精| 精品国产一区二区久久| 亚洲av成人不卡在线观看播放网 | avwww免费| 亚洲国产欧美在线一区| 日韩电影二区| 多毛熟女@视频| 午夜视频精品福利| 岛国毛片在线播放| 国产日韩欧美在线精品| 狠狠狠狠99中文字幕| 午夜91福利影院| 国产淫语在线视频| 大型av网站在线播放| 男女高潮啪啪啪动态图| 伊人久久大香线蕉亚洲五| a级毛片在线看网站| 午夜日韩欧美国产| 1024香蕉在线观看| 成人国产av品久久久| 丰满迷人的少妇在线观看| 精品人妻一区二区三区麻豆| 国产精品偷伦视频观看了| 老司机亚洲免费影院| 国产精品.久久久| 精品一区在线观看国产| 日本91视频免费播放| 麻豆国产av国片精品| av免费在线观看网站| 99精国产麻豆久久婷婷| 我的亚洲天堂| 久热这里只有精品99| 免费不卡黄色视频| 黄网站色视频无遮挡免费观看| 母亲3免费完整高清在线观看| 丝袜脚勾引网站| 最近最新免费中文字幕在线| 在线观看www视频免费| 国产熟女午夜一区二区三区| 国产深夜福利视频在线观看| 狂野欧美激情性xxxx| 日韩视频在线欧美| 国产片内射在线| 亚洲国产精品999| 久久99热这里只频精品6学生| 美女高潮到喷水免费观看| 色播在线永久视频| 国产免费福利视频在线观看| 老熟妇乱子伦视频在线观看 | 丰满少妇做爰视频| 男女下面插进去视频免费观看| 老熟妇乱子伦视频在线观看 | 女性被躁到高潮视频| 母亲3免费完整高清在线观看| 亚洲精品成人av观看孕妇| 国产色视频综合| 在线亚洲精品国产二区图片欧美| 色婷婷av一区二区三区视频| 免费少妇av软件| 法律面前人人平等表现在哪些方面 | 国产精品偷伦视频观看了| 亚洲精品中文字幕在线视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲成人免费av在线播放| 欧美一级毛片孕妇| 久久天堂一区二区三区四区| 亚洲色图综合在线观看| 精品国产一区二区三区久久久樱花| 欧美人与性动交α欧美软件| 午夜福利在线观看吧| 啦啦啦啦在线视频资源| 日本一区二区免费在线视频| 欧美精品高潮呻吟av久久| 国产成人啪精品午夜网站| a级片在线免费高清观看视频| 久久久久国产一级毛片高清牌| 美女中出高潮动态图| 国产片内射在线| 捣出白浆h1v1| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产精品一区二区三区在线| 亚洲欧美激情在线| 精品熟女少妇八av免费久了| av在线老鸭窝| 久久ye,这里只有精品| 成年动漫av网址| av又黄又爽大尺度在线免费看| 成人18禁高潮啪啪吃奶动态图| 男人爽女人下面视频在线观看| videosex国产| 黄色怎么调成土黄色| 日本欧美视频一区| 亚洲一区二区三区欧美精品| 精品久久久精品久久久| 99精品欧美一区二区三区四区| 久久精品久久久久久噜噜老黄| 动漫黄色视频在线观看| 1024视频免费在线观看| 在线十欧美十亚洲十日本专区| 欧美日韩亚洲综合一区二区三区_| 亚洲 欧美一区二区三区| 深夜精品福利| 久久久精品94久久精品| 精品亚洲成a人片在线观看| 69精品国产乱码久久久| 亚洲精品一卡2卡三卡4卡5卡 | 麻豆国产av国片精品| 亚洲精品在线美女| 成人手机av| 另类精品久久| 国产一区二区三区av在线| 韩国高清视频一区二区三区| 视频在线观看一区二区三区| 黄色视频在线播放观看不卡| 亚洲人成77777在线视频| 欧美精品av麻豆av| 国产精品久久久av美女十八| 热99re8久久精品国产| 女人被躁到高潮嗷嗷叫费观| tocl精华| 婷婷色av中文字幕| 国产欧美日韩一区二区精品| 欧美日韩黄片免| 波多野结衣一区麻豆| 啦啦啦啦在线视频资源| 免费高清在线观看视频在线观看| 9热在线视频观看99| 大型av网站在线播放| 91精品伊人久久大香线蕉| 久久天堂一区二区三区四区| 高清av免费在线| 免费黄频网站在线观看国产| 免费不卡黄色视频| 男人舔女人的私密视频| 亚洲天堂av无毛| 热re99久久精品国产66热6| 青草久久国产| 黄网站色视频无遮挡免费观看| 18在线观看网站| 亚洲成人免费电影在线观看| 免费高清在线观看视频在线观看| 秋霞在线观看毛片| a在线观看视频网站| 亚洲三区欧美一区| 久久性视频一级片| 啦啦啦啦在线视频资源| 成年人黄色毛片网站| 欧美日韩黄片免| 亚洲天堂av无毛| 飞空精品影院首页| 中文字幕人妻熟女乱码| 亚洲全国av大片| 国产精品秋霞免费鲁丝片| 亚洲全国av大片| 在线 av 中文字幕| 99国产精品一区二区蜜桃av | 国产成人免费无遮挡视频| 99国产精品一区二区蜜桃av | 天堂俺去俺来也www色官网| 菩萨蛮人人尽说江南好唐韦庄| 午夜福利一区二区在线看| 黄色毛片三级朝国网站| 国产福利在线免费观看视频| 午夜视频精品福利| 国产激情久久老熟女| www.av在线官网国产| 国产老妇伦熟女老妇高清| 69精品国产乱码久久久| 国产xxxxx性猛交| 人人妻人人添人人爽欧美一区卜| 精品免费久久久久久久清纯 | 一边摸一边做爽爽视频免费| 国产麻豆69| 久久天堂一区二区三区四区| 天天躁日日躁夜夜躁夜夜| 丝袜人妻中文字幕| 国产亚洲av片在线观看秒播厂| 丰满迷人的少妇在线观看| 久久人妻福利社区极品人妻图片| 在线观看www视频免费| 少妇的丰满在线观看| 男女床上黄色一级片免费看| 热99国产精品久久久久久7| 亚洲美女黄色视频免费看| 午夜影院在线不卡| 久久国产精品人妻蜜桃| 国产一区二区三区综合在线观看| 亚洲成人手机| 一级片免费观看大全| 老熟女久久久| av网站在线播放免费| 久久毛片免费看一区二区三区| 天堂8中文在线网| 肉色欧美久久久久久久蜜桃| 丝袜在线中文字幕| 午夜激情久久久久久久| 成人免费观看视频高清| 精品福利观看| 日韩一卡2卡3卡4卡2021年| 亚洲第一青青草原| av不卡在线播放| 国产一区二区三区av在线| 日韩 欧美 亚洲 中文字幕| 国产成人免费无遮挡视频| 男女边摸边吃奶| 精品亚洲成a人片在线观看| 人人澡人人妻人| 国产成人精品无人区| 成人国产一区最新在线观看| 欧美大码av| 精品国产乱码久久久久久男人| 亚洲精品中文字幕在线视频| 免费久久久久久久精品成人欧美视频| 成人亚洲精品一区在线观看| 黄网站色视频无遮挡免费观看| 成年女人毛片免费观看观看9 | 天天躁夜夜躁狠狠躁躁| 男女高潮啪啪啪动态图| 日韩欧美免费精品| 午夜视频精品福利| 久久久久精品国产欧美久久久 | 国产成人a∨麻豆精品| 日韩大码丰满熟妇| 国产精品久久久av美女十八| 日日摸夜夜添夜夜添小说| 久久精品成人免费网站| 波多野结衣av一区二区av| xxxhd国产人妻xxx| 一级,二级,三级黄色视频| 在线观看人妻少妇| 一区二区三区四区激情视频| 视频在线观看一区二区三区| tube8黄色片| 黄色毛片三级朝国网站| 18在线观看网站| 亚洲第一av免费看| 国产成人精品久久二区二区免费| 亚洲五月婷婷丁香| 日韩电影二区| 在线观看舔阴道视频| 欧美精品一区二区免费开放| 国产精品香港三级国产av潘金莲| 99九九在线精品视频| 国产欧美日韩一区二区精品| av视频免费观看在线观看| 又黄又粗又硬又大视频| 一进一出抽搐动态| 国产伦理片在线播放av一区| 成年人午夜在线观看视频| 99久久精品国产亚洲精品| 国产欧美日韩精品亚洲av|