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

    掘進(jìn)擾動(dòng)下斷層帶煤體瓦斯異常分布特征研究

    2022-11-09 02:29:16馬桂霞周愛桃劉賽菲
    煤礦安全 2022年10期
    關(guān)鍵詞:壓力梯度均質(zhì)煤體

    周 勇,馬桂霞,周愛桃,劉賽菲

    (1.國(guó)家能源集團(tuán)烏海能源有限責(zé)任公司,內(nèi)蒙古烏海 016099;2.德礦事業(yè)(北京)科技有限公司,北京 100086;3.中國(guó)礦業(yè)大學(xué)(北京)應(yīng)急管理與安全工程學(xué)院,北京 100083)

    煤與瓦斯突出危害性大,突發(fā)性強(qiáng),是煤礦開采中最為嚴(yán)重的煤巖動(dòng)力災(zāi)害之一[1]。而我國(guó)大部分是富含瓦斯的石炭、二疊紀(jì)煤層,高瓦斯和煤與瓦斯突出煤層多,這也導(dǎo)致我國(guó)成為世界上煤與瓦斯突出災(zāi)害最嚴(yán)重的國(guó)家[2]。據(jù)相關(guān)統(tǒng)計(jì)資料表明,煤與瓦斯突出發(fā)生的地點(diǎn)在掘進(jìn)工作面的概率高達(dá)70%~80%[3]。從前人針對(duì)煤與瓦斯突出機(jī)理的研究現(xiàn)狀來看,普遍認(rèn)為[4-8]:受掘進(jìn)應(yīng)力的影響,前方煤體物理力學(xué)性質(zhì)及滲透率發(fā)生改變,進(jìn)而導(dǎo)致瓦斯場(chǎng)發(fā)生改變。煤巖體破壞影響煤基質(zhì)和裂隙內(nèi)部瓦斯運(yùn)移規(guī)律,同樣其也是煤與瓦斯突出的能量來源。當(dāng)煤體內(nèi)部應(yīng)力達(dá)到一定的強(qiáng)度且具有較大的瓦斯壓力梯度,就將具備較高的煤與瓦斯突出性。

    同時(shí)在巷道掘進(jìn)過程中,斷層影響區(qū)發(fā)生煤與瓦斯突出的概率又遠(yuǎn)超其他區(qū)域,且造成的危害更大[9]。因此,對(duì)于斷層帶的存在對(duì)煤與瓦斯突出的影響機(jī)制方面,前人們也進(jìn)行了諸多研究,得到了許多能夠指導(dǎo)現(xiàn)場(chǎng)安全生產(chǎn)的成果。程遠(yuǎn)平等[10]、涂慶毅[11]、舒龍勇[12]、高魁[13]等以斷層附近構(gòu)造煤為研究對(duì)象,從地質(zhì)結(jié)合構(gòu)造演化及應(yīng)力場(chǎng)角度研究構(gòu)造煤綜合因素對(duì)煤與瓦斯突出的控制機(jī)理;姜耀東等[14]、任政等[15]、郭曉強(qiáng)等[16]、王愛文等[17]、郭德勇等[18]以采掘擾動(dòng)誘發(fā)的斷層活化為研究對(duì)象,分析斷層活化過程,提出突出過程中的黏滑失穩(wěn)現(xiàn)象,進(jìn)一步挖掘斷層活化誘發(fā)煤巖動(dòng)力災(zāi)害的黏滑失穩(wěn)機(jī)理;Yao Yanbin 等[19]、魏國(guó)營(yíng)等[20]、Cheng Zhai 等[21]、高魁等[22]從瓦斯地質(zhì)及瓦斯賦存規(guī)律入手,結(jié)合區(qū)域構(gòu)造演化理論,研究了不同型式構(gòu)造活動(dòng)對(duì)瓦斯地質(zhì)特征、瓦斯賦存及煤與瓦斯突出的控制機(jī)理。

    斷層是煤礦開采中最常見的地質(zhì)構(gòu)造,在地質(zhì)學(xué)中也有較多的研究,通常具有特定的結(jié)構(gòu),包括斷層核心區(qū)和損傷區(qū)域[23-24]。斷層核心區(qū)域通常為低滲透充填的斷層泥,通常是阻礙瓦斯運(yùn)移的屏障。而斷層損傷區(qū)內(nèi)部遍布大量的裂縫和斷裂帶,通常是瓦斯運(yùn)移的通道,斷層內(nèi)部結(jié)構(gòu)會(huì)表現(xiàn)出明顯的非均質(zhì)性[25-26]。

    從研究現(xiàn)狀來看,在掘進(jìn)工作面前方遇斷層是發(fā)生煤與瓦斯突出概率最高的地質(zhì)條件,國(guó)內(nèi)外學(xué)者也對(duì)其進(jìn)行了大量的研究。但大多研究或是集中在以構(gòu)造煤為研究對(duì)象,或?qū)⑵浜?jiǎn)化成接觸面,都忽略了斷層復(fù)雜的內(nèi)部結(jié)構(gòu),因此無法準(zhǔn)確系統(tǒng)地描述斷層帶內(nèi)部應(yīng)力場(chǎng)及瓦斯運(yùn)移變化對(duì)誘發(fā)掘進(jìn)面煤與瓦斯突出的影響機(jī)制。由于斷層帶內(nèi)部這種顯著的宏觀各向異性,會(huì)控制內(nèi)部瓦斯的運(yùn)移方向,使瓦斯運(yùn)移主要集中在斷層核兩側(cè)的損傷區(qū)域內(nèi),傾向于平行于斷層面運(yùn)移[27]。為此,基于斷層帶滲透性宏觀上的各向異性,將斷層看成橫向各向同性,揭示煤巷掘進(jìn)前方非均質(zhì)斷層帶瓦斯運(yùn)移規(guī)律及其對(duì)煤與瓦斯突出事故的誘發(fā)機(jī)制。

    1 基于斷層帶結(jié)構(gòu)異性的滲透模型

    1.1 斷層帶結(jié)構(gòu)異性概念模型

    斷層斷裂帶結(jié)構(gòu)往往伴隨長(zhǎng)時(shí)間的構(gòu)造發(fā)育期,是逐漸而非某次劇烈變化形成的,這也是斷裂帶結(jié)構(gòu)在幾何空間上復(fù)雜性的原因。斷層斷裂帶是具有一定寬度、長(zhǎng)度、高度的復(fù)雜內(nèi)部結(jié)構(gòu)的三維地質(zhì)體。對(duì)于斷裂帶的內(nèi)部結(jié)構(gòu)簡(jiǎn)化的概念模型,是在近20 年里發(fā)展起來,是典型的二元結(jié)構(gòu),即損傷區(qū)域和內(nèi)嵌的斷層內(nèi)核[28]。它涉及斷層核的應(yīng)變,且在斷層核周圍是布滿裂縫和斷裂帶的損傷區(qū)域。斷層核是指由斷層滑動(dòng)面、構(gòu)造透鏡體和斷層巖充填部分組成的剪切帶,是強(qiáng)烈區(qū)域應(yīng)變和剪切的結(jié)果[29]。根據(jù)斷層核的數(shù)量不同,又可將斷層分為以龐奇博爾(Punchbowl)斷層為代表的單斷層核斷層破碎帶和以卡沃內(nèi)拉斯(Carboneras)斷層為代表的多斷層核的復(fù)合斷層破碎帶。損傷區(qū)域是指斷層核外側(cè)的低序次構(gòu)造(包括導(dǎo)致斷層斷裂帶的滲透性和力學(xué)特性的非均質(zhì)性的次級(jí)斷層褶皺發(fā)育帶、裂隙和巖脈等)切割圍巖形成的異常區(qū)域[30]。典型斷裂帶構(gòu)造示意圖如圖1。

    圖1 典型斷裂帶構(gòu)造示意圖Fig.1 Schematic diagrams of typical fault zone structure

    斷層結(jié)構(gòu)的非均質(zhì)性是導(dǎo)致斷層內(nèi)部力學(xué)特征及滲透特性非均質(zhì)性的主要原因。重點(diǎn)研究基于斷層結(jié)構(gòu)非均質(zhì)性的各向異性滲透模型及瓦斯運(yùn)移規(guī)律,并進(jìn)一步研究其對(duì)煤與瓦斯突出的影響作用機(jī)制。利用前人在建立數(shù)學(xué)微分方程所時(shí)依據(jù)的合理假設(shè)和定律,在模型建立之前做出以下必要假設(shè):①煤體物理結(jié)構(gòu)是各向同性的,為雙孔隙-彈性連續(xù)體;②孔隙內(nèi)瓦斯為理想氣體,且在煤體中的流動(dòng)遵循達(dá)西定律;③條件是等溫的;④將斷層簡(jiǎn)化為理想的橫向各向同性彈塑性體模型。

    1.2 橫向各向同性斷層滲透率模型

    地層的各向異性通常分為4 個(gè)大類,分別是:三斜各向異性、單斜各向異性、正交各向異性和橫向各向同性[31]。斷層內(nèi)部具有明顯的宏觀各向異性,具有良好的分層結(jié)構(gòu),呈現(xiàn)橫向各向同性特征。橫向各向同性具有以下3 個(gè)特征:①關(guān)于垂直于橫向各向同性平面的軸軸對(duì)稱;②橫向平面是具有對(duì)稱性的無限平面,沿著對(duì)稱軸方向上的力學(xué)及滲透性質(zhì)不同;③在包含對(duì)稱軸的面上,各個(gè)方向的力學(xué)及滲透性質(zhì)相同。其中,具有垂直對(duì)稱軸的橫向各向同性地層稱為VTI 地層,傾斜對(duì)稱軸的橫向各同性地層稱為TTI 地層[32],橫向各向同性地層如圖2。其中VTI 地層是當(dāng)對(duì)稱軸與z 軸夾角θ=0 時(shí)TTI 地層的1 種特殊情況。

    圖2 橫向各向同性地層Fig.2 Transversely isotropic strata

    假設(shè)垂直于斷層方向?yàn)? 個(gè)滲透率的主值方向,同時(shí)由于在每一個(gè)包含垂直于橫向各向同性平面的軸的平面內(nèi),各方向的力學(xué)及滲透性質(zhì)相同,則規(guī)定其滲透率主值為Kz。因此,斷層帶滲透率橫向各向同性的三維問題就可以轉(zhuǎn)化為二維平面問題進(jìn)行研究分析。

    在二維平面空間中,滲透率張量Kij表示為:

    當(dāng)非均質(zhì)地層為VTI 地層時(shí),可將式(1)進(jìn)行對(duì)角化處理,其中Kxy=Kyx,可得到對(duì)角線張量的表達(dá)式:

    式中:Kx、Ky為滲透率的主值,對(duì)應(yīng)主滲流方向。

    透率主值方向與滲透率分量方向的坐標(biāo)轉(zhuǎn)化原理圖如圖3。

    圖3 滲透率主值與滲透率分量的坐標(biāo)轉(zhuǎn)化原理圖Fig.3 Coordinate transformation principle diagram of permeability principal value and permeability component

    根據(jù)轉(zhuǎn)化原理圖,可以求得滲透率主值與滲透率分量的關(guān)系[33]:

    式中:θ 為轉(zhuǎn)化方位角,在已知滲透率主值與對(duì)應(yīng)方位角的情況下可以求得任意坐標(biāo)系下的滲透率分量。

    1.3 含瓦斯煤流固耦合模型

    1.3.1 基于雙重孔隙彈性理論的煤巖變形方程

    考慮煤中吸附瓦斯及游離瓦斯對(duì)煤體變形及有效應(yīng)力的作用[34],雙重孔隙介質(zhì)的應(yīng)力平衡方程為:

    式中:σij,j為有效應(yīng)力張量,MPa;δij為Kroneker數(shù);αf、αm為考慮裂隙和孔隙的有效應(yīng)力系數(shù);pf、pm為裂隙和煤基質(zhì)孔隙中的氣體壓力,MPa;σs,i為氣體吸附膨脹應(yīng)力,MPa;Fi為體積力分量,MPa。

    式中:K 為煤的體積模量,K=E/3(1-2μ),MPa;Km為煤基質(zhì)的體積模量,Km=Em/3(1-2μ),MPa;Ks為煤骨架的體積模量,Ks=Km/[1-3?m(1-μ)/(2-4μ)],MPa;E 為煤體彈性模量,MPa;Em為煤基質(zhì)彈性模量,MPa;?m為煤基質(zhì)孔隙率;μ 為煤的泊松比。

    對(duì)于各向同性連續(xù)介質(zhì)變形的幾何方程為[35]:

    式中:ui,j、uj,i分別為軸向和徑向位移張量;εij為應(yīng)變張量。

    基于均勻各向同性介質(zhì)在三維應(yīng)力作用下的變形遵循廣義胡克定律,含瓦斯煤的變形本構(gòu)方程為[36]:

    式中:G 為煤的剪切模量,MPa,G=E/2(1+υ);εV為煤的體積應(yīng)變;δij為總應(yīng)力,MPa;εm為煤基質(zhì)的變形。

    選取DP 準(zhǔn)則匹配Mohr-Coulomb 準(zhǔn)則作為煤巖體屈服破壞準(zhǔn)則,破壞系數(shù)F 為。

    式中:I1為應(yīng)力張量的第1 不變量;J2為應(yīng)力偏量的第2 不變量;C 為黏聚力;φ 為內(nèi)摩擦角。

    1.3.2 煤基質(zhì)瓦斯擴(kuò)散控制方程

    瓦斯在煤層中的擴(kuò)散主要為Fick 型擴(kuò)散,擴(kuò)散的原因是煤基質(zhì)孔隙中與煤體裂隙中氣體的濃度差。根據(jù)質(zhì)量平衡定律,單位時(shí)間內(nèi)微孔單元的質(zhì)量變化量等于煤微孔單元各方向上基質(zhì)與裂隙間擴(kuò)散輸送的氣體量。因此,煤基質(zhì)氣體質(zhì)量守恒方程可以表示為:

    式中:Q 為煤微孔單元各方向上的基質(zhì)與裂隙間擴(kuò)散輸送的氣體量;t 為時(shí)間,s;mm為煤基質(zhì)孔隙中單位體積含氣量,包括孔隙吸附氣體量和游離氣體量[37]。

    式中:VL為朗繆爾體積常數(shù),m3/kg;Mg為甲烷的摩爾質(zhì)量,kg/mol;PL為L(zhǎng)angmuir 壓力常數(shù),孔隙吸附量是最大吸附量的一半時(shí)的吸附平衡壓力,MPa;Vg為標(biāo)準(zhǔn)條件下甲烷的摩爾體積,m3/mol;ρc為煤的視密度,kg/m3;?m為煤基質(zhì)孔隙率;R 為理想氣體含量,8.314 3 J/(mol·K);T 為煤體溫度,K。

    而煤基質(zhì)瓦斯向裂隙擴(kuò)散的源項(xiàng)通量可由下式表示[38]:

    式中:τ 為吸附時(shí)間,h。

    聯(lián)立式(10)~式(12)可得煤基體瓦斯壓力隨時(shí)間的變化控制方程為[39]:

    1.3.3 瓦斯?jié)B流控制方程

    對(duì)于瓦斯在煤中的運(yùn)移,一般認(rèn)為煤裂隙是瓦斯的滲流通道,裂隙中的瓦斯運(yùn)移可以認(rèn)為遵循達(dá)西定律。根據(jù)質(zhì)量守恒定律,對(duì)于單位體積的煤,單位時(shí)間內(nèi)流出裂縫系統(tǒng)的氣體質(zhì)量與在煤基質(zhì)系統(tǒng)中擴(kuò)散的氣體質(zhì)量之和等于裂縫系統(tǒng)中自由氣體的質(zhì)量隨時(shí)間的變化量。因此,煤層中瓦斯運(yùn)動(dòng)的連續(xù)性方程可以表示為:

    式中:v 為裂縫中氣體的滲流速度,m/s;u 為氣體的動(dòng)力黏度,Pa·s;ρg為氣體密度,kg/m3;?f為煤的孔隙度;kf為煤滲透率,mD(1 mD=40 m2/(MPa2·d)。

    結(jié)合式(12)、式(15)、式(16)代入到式(14)中可以得出裂縫中瓦斯流動(dòng)的控制方程可以表示為:

    1.3.4 掘進(jìn)應(yīng)力影響下多應(yīng)力路徑滲透率模型

    考慮了氣體對(duì)煤體介質(zhì)的力學(xué)響應(yīng)和多孔介質(zhì)力學(xué)性質(zhì)的變化對(duì)氣體運(yùn)移特性的影響。由于煤是雙重孔隙介質(zhì)模型,當(dāng)考慮瓦斯吸附引起的膨脹應(yīng)力時(shí),含瓦斯煤的有效應(yīng)力模型就可表示為:

    式中:σeij、σij分別為有效應(yīng)力張量和總應(yīng)力張量,MPa;σs為吸附膨脹應(yīng)力,MPa。

    σs可以表示為[40]:

    式中:a 為參考?jí)毫ο碌臉O限吸附量,m3/t;b 為煤的吸附平衡常數(shù),MPa-1;ρc為煤的密度,t/m3;Vm為氣體摩爾體積,m3/mol,取22.4×10-3m3/mol。

    當(dāng)受巷道掘進(jìn)引起的集中應(yīng)力影響時(shí),煤體受到屈服應(yīng)力作用,由彈性階段進(jìn)入塑性階段。由于煤體內(nèi)部發(fā)生塑性破壞,裂隙開始擴(kuò)展。在此階段,煤體滲透性呈現(xiàn)突變趨勢(shì),進(jìn)入殘余強(qiáng)度階段后逐漸穩(wěn)定。煤在上述過程中的滲透率演化非常復(fù)雜,很難通過嚴(yán)格的數(shù)學(xué)推導(dǎo)得到滲透率模型。因此,根據(jù)試驗(yàn)結(jié)果和前人的研究結(jié)果,假定裂隙數(shù)目與偏應(yīng)力滿足簡(jiǎn)單的線性關(guān)系,得出在掘進(jìn)應(yīng)力影響下多應(yīng)力路徑的孔隙率及滲透率演化模型[36]:

    式中:kf0、?f0為初始滲透率和孔隙度,mD;p0為煤初始瓦斯壓力;σ1、σ3為最大最小主應(yīng)力;σ0為初始地應(yīng)力;γ′為煤體新生裂隙數(shù)目對(duì)偏應(yīng)力的敏感性;f 為有效煤基質(zhì)變形因子,0~1;εL為煤體Langmuir 應(yīng)變常數(shù);△σ 為外部應(yīng)力變化量,MPa;σ0為初始地應(yīng)力;p0為煤初始瓦斯壓力。

    2 數(shù)值模擬計(jì)算及結(jié)果分析

    2.1 幾何模型

    建立了考慮有效應(yīng)力和煤吸附膨脹變形作用的煤層滲透率方程、非均質(zhì)斷層區(qū)域的各向異性滲透率方程、含瓦斯煤變形控制方程和雙重孔隙的煤體瓦斯運(yùn)移方程的氣固耦合模型。由于所涉及方程較多又具有很強(qiáng)的非線性,因此整個(gè)解算過程中,煤巖應(yīng)力應(yīng)變及屈服過程等問題采用COMSOL 中固體力學(xué)模塊求解,瓦斯擴(kuò)散和滲流過程等問題采用偏微分方程PDE 模塊求解,煤與瓦斯間的相互作用通過模塊間參數(shù)傳遞進(jìn)行實(shí)現(xiàn),然后根據(jù)數(shù)值計(jì)算結(jié)果分析采掘工作面前方煤體受非均質(zhì)斷層影響表現(xiàn)出的力學(xué)及瓦斯運(yùn)移響應(yīng)。模型以含斷層煤層為背景,并通過簡(jiǎn)化煤層結(jié)構(gòu)建立的二維掘進(jìn)面前方遇斷層的幾何模型如圖4。

    圖4 數(shù)值計(jì)算幾何模型及邊界條件Fig.4 Numerical calculation geometric model and boundary conditions

    模型全區(qū)域在垂直方向高30 m,水平方向上寬45 m。煤層區(qū)域位于2 個(gè)巖層之間,煤厚6 m。斷層區(qū)域?qū)?5 m,斷層落差為3 m,斷層傾角為53°,斷層類型為正斷層。圖中自由度數(shù)目為16 556;網(wǎng)格單元數(shù)4 087;邊界單元數(shù)367;頂點(diǎn)元素?cái)?shù)18;最小單元質(zhì)量0.596 7;單元面積比0.036 06;平均單元質(zhì)量為0.875 7。COMSOL 中用質(zhì)量為1 來描述理想單元,在二維幾何中,一般要求最小網(wǎng)格質(zhì)量為0.3 以上,該模型網(wǎng)格劃分滿足要求。

    2.2 邊界條件

    1)應(yīng)力場(chǎng)邊界條件。煤層和頂?shù)装鍍?nèi)部邊界采用自由邊界條件;模型底部邊界采用固定約束;左右兩側(cè)邊界采用輥支承約束,固定橫向位移;上部邊界設(shè)定為應(yīng)力邊界,應(yīng)力為上邊界承受的上覆巖層的自重力。

    2)滲流場(chǎng)邊界條件。工作面前方暴露煤壁采用Dirichlet 邊界條件:瓦斯壓力為0.1 MPa(大氣壓);煤層內(nèi)部初始瓦斯壓力為1.8 MPa;其他邊界均設(shè)定為0 通量的Neumann 邊界條件:np=q(n 為邊界的法向量;p 為外部壓力;q 為流通量)。模型所涉及的其他物理力學(xué)參數(shù)見表1。

    表1 某礦5#煤基本物理力學(xué)參數(shù)Table 1 Basic physical and mechanical parameters of 5#coal in a mine

    2.3 模擬方案

    工況方案:在采掘工作面位于斷層附近處,煤層初始瓦斯壓力為1.8 MPa 時(shí),分析煤體的應(yīng)力、應(yīng)變狀態(tài)以及非均質(zhì)斷層對(duì)煤基質(zhì)瓦斯壓力、裂隙瓦斯壓力以及滲透率的影響變化規(guī)律,并將其與均質(zhì)斷層進(jìn)行對(duì)比分析。最后得出非均質(zhì)斷層對(duì)煤體的應(yīng)力場(chǎng)及瓦斯運(yùn)移響應(yīng)的影響規(guī)律。

    3 煤礦區(qū)煤層掘進(jìn)的多物理場(chǎng)耦合效應(yīng)分析

    3.1 掘進(jìn)工作面前方應(yīng)力及位移物理場(chǎng)分布

    埋深恒定(300 m),煤體溫度恒定(303.15 K),且煤層掘進(jìn)5 m 時(shí),當(dāng)煤層內(nèi)部瓦斯壓力為1.8 MPa 時(shí)上覆巖層及暴露面前方煤體的等效應(yīng)力、垂直應(yīng)力、垂直位移及塑性應(yīng)變變化云圖如圖5。

    圖5 掘進(jìn)暴露面前方各力學(xué)物理參數(shù)分布Fig.5 Distribution of mechanical and physical parameters in front of tunneling exposure surface

    由圖5 可以看出,煤層開挖后,工作面前方煤體及頂?shù)装鍛?yīng)力平衡被打破。受采動(dòng)影響,采掘暴露面前方煤體依次劃分為應(yīng)力卸載帶、應(yīng)力集中帶和原巖應(yīng)力帶。由圖5(a)、圖5(b)可以看出沿采掘工作面前方,采空區(qū)上覆巖層及底板巖層出現(xiàn)應(yīng)力卸壓帶,上覆巖層壓力轉(zhuǎn)移到由采掘工作面前方煤體承擔(dān),出現(xiàn)應(yīng)力集中帶,同時(shí)煤體的破壞呈球殼狀,且距離采掘暴露面大約3~4 m 處出現(xiàn)應(yīng)力峰值,垂直應(yīng)力峰值處于應(yīng)力集中帶內(nèi);由圖5(c)可以看出,隨著采掘工作面從開切眼的不斷推進(jìn),煤層上方覆巖出現(xiàn)不同程度的沉降值,從圖中明顯可以看出煤層頂板沉降值在采空區(qū)上方最大,隨煤體走向逐漸減小。在采掘過程中煤層底板也會(huì)出現(xiàn)一定的下沉量,但與煤層頂板相比明顯較?。挥蓤D5(d)可以看出,受采掘擾動(dòng)影響,暴露面前方煤體的應(yīng)力狀態(tài)會(huì)發(fā)生變化,煤體會(huì)產(chǎn)生屈服破壞,出現(xiàn)塑性變形。從圖中可以看出煤體塑性變形區(qū)域的范圍大概延續(xù)到采掘暴露面前方5 m 處。且在煤層與斷層的交界面處由于兩者力學(xué)特性的強(qiáng)弱差異,出現(xiàn)塑性區(qū)階躍的現(xiàn)象。煤體應(yīng)力、破壞程度的改變以及斷層的影響都會(huì)對(duì)煤體滲透特性產(chǎn)生較大的變化,從而進(jìn)一步控制瓦斯壓力及瓦斯運(yùn)移的變化。

    3.2 掘進(jìn)工作面前方(含斷層帶)瓦斯分布

    為了分析受采掘擾動(dòng)以及非均質(zhì)斷層的影響,煤基質(zhì)、裂隙瓦斯壓力以及煤層滲透率的影響變化規(guī)律,以暴露面中心點(diǎn)起,沿?cái)鄬訋ё呦颍帘┞睹媲胺?0 m 煤層中心點(diǎn)止作為觀測(cè)線,如圖4 中所示黃線位置。煤層開挖5 m 后煤壁前方滲透率,煤基質(zhì)、裂隙瓦斯壓力,塑性應(yīng)變和體積應(yīng)變隨著與暴露面距離變化的演化規(guī)律如圖6。

    圖6 滲透率、裂隙瓦斯壓力、塑性應(yīng)變和體積應(yīng)變從暴露面煤壁到右邊界的變化規(guī)律Fig.6 Variation of permeability, gas pressure in fracture,plastic strain and volumetric strain from exposed coal wall to right boundary

    由圖6 可以看出,受采動(dòng)影響,越靠近暴露面煤壁,煤體塑性應(yīng)變?cè)酱螅w積應(yīng)變?cè)酱?,?dǎo)致煤體裂隙孔隙的突增,使煤體滲透率呈現(xiàn)幾十甚至幾百倍的增長(zhǎng)。沿煤層走向方向,在應(yīng)力集中帶內(nèi),會(huì)受到應(yīng)力集中對(duì)煤體的壓縮變形作用,煤體逐漸被壓實(shí),體積應(yīng)變降低到最小值。隨著煤層走向逐漸到原巖應(yīng)力帶,應(yīng)力逐漸減少,體積應(yīng)變出現(xiàn)回升現(xiàn)象,最終回升到1 個(gè)穩(wěn)定值。煤層滲透率的變化與采掘應(yīng)力及煤體塑性變形之間具有密切的關(guān)系,在鄰近采掘暴露面煤壁附近,未受到水平方向的約束,致使?jié)B透率達(dá)到峰值。在應(yīng)力集中區(qū)內(nèi),由于應(yīng)力集中裂隙發(fā)生閉合,導(dǎo)致滲透率降低,在應(yīng)力峰值處降至最低。同時(shí)在應(yīng)力峰值附近出現(xiàn)屈服破壞,沿煤層走向應(yīng)力逐漸降低,逐漸恢復(fù)至原巖應(yīng)力值,滲透率出現(xiàn)小幅上升。在未受采動(dòng)影響的原巖應(yīng)力區(qū)域,滲透率恢復(fù)至初始值。斷層帶內(nèi)滲透率出現(xiàn)先升高然后激降,再激增,最后緩降的現(xiàn)象。這是由于受到斷層核附近的異常高應(yīng)力集中及其內(nèi)部分布著大量低滲透性斷層泥的影響,致使?jié)B透率出現(xiàn)從斷層帶兩側(cè)損傷區(qū)域到斷層核過渡的驟降現(xiàn)象,使得斷層核成為阻礙瓦斯向采掘空間逸散的封閉空間。斷層核內(nèi)部滲透率為6.438 3×10-18m2明顯低于煤層滲透率,約為煤層初始滲透率的1/7 倍。而作為瓦斯運(yùn)移通道的斷層損傷區(qū)內(nèi)部滲透率卻要高于煤層滲透率。煤基質(zhì)和裂隙瓦斯壓力在暴露面附近均為大氣壓,且都沿著煤層走向增加,最終至原巖應(yīng)力區(qū)恢復(fù)至初始瓦斯壓力。鄰近掘進(jìn)暴露面附近煤基質(zhì)和裂隙都具有較大的瓦斯壓力梯度,且沿煤層走向至卸壓區(qū),兩者瓦斯壓力梯度都會(huì)降低,同時(shí)煤層滲透率較大,安全性較高。應(yīng)力集中帶內(nèi)瓦斯壓力梯度增加,煤基質(zhì)和裂隙瓦斯壓力梯度都大于卸壓區(qū),且煤體受集中應(yīng)力及塑性變形影響,滲透率較低,蘊(yùn)育煤與瓦斯突出危險(xiǎn)性較高。雖然煤基質(zhì)和裂隙具有一致的瓦斯壓力分布,但裂隙系統(tǒng)明顯具有較強(qiáng)的滲流能力,煤基質(zhì)中的瓦斯解吸運(yùn)移速度較慢,可以一直向裂隙系統(tǒng)中提供瓦斯來源,使煤體能夠保持較高的瓦斯壓力梯度。在斷層帶內(nèi),瓦斯壓力梯度呈現(xiàn)先增加后降低的趨勢(shì),且斷層核相較于斷層帶兩側(cè)損傷區(qū)域瓦斯壓力梯度明顯較大。此時(shí),若將斷層帶區(qū)域看成均質(zhì)性斷層,就會(huì)得出發(fā)生煤與瓦斯突出的危險(xiǎn)性概率降低的結(jié)論,但這與現(xiàn)場(chǎng)實(shí)際情況明顯不符。

    3.3 橫向各向同性斷層對(duì)突出的影響作用機(jī)制

    為研究非均質(zhì)斷層與均質(zhì)斷層對(duì)煤與瓦斯突出不同的影響程度,分別采用非均質(zhì)橫向各向同性滲透率模型和均質(zhì)滲透率模型對(duì)斷層滲透特性進(jìn)行描述,將兩者影響下的煤體瓦斯壓力進(jìn)行對(duì)比分析,進(jìn)而得出不同滲透率模型對(duì)突出的影響作用機(jī)制。不同斷層性質(zhì)影響下裂隙瓦斯壓力分布云圖如圖7,不同斷層性質(zhì)影響下煤基質(zhì)瓦斯壓力分布云圖如圖8,不同斷層性質(zhì)影響下裂隙瓦斯壓力及煤基質(zhì)瓦斯壓力變化規(guī)律圖如圖9。

    圖7 不同斷層性質(zhì)影響下裂隙瓦斯壓力分布云圖Fig.7 Diagrams of fracture gas pressure distribution under the influence of different fault properties

    圖8 不同斷層性質(zhì)影響下煤基質(zhì)瓦斯壓力分布云圖Fig.8 Distribution diagrams of coal matrix gas pressure under the influence of different fault properties

    圖9 不同斷層性質(zhì)影響下瓦斯壓力變化規(guī)律圖Fig.9 Gas pressure variation under the influence of different fault properties

    煤層中瓦斯壓力梯度的大小對(duì)煤與瓦斯突出具有至關(guān)重要的作用,瓦斯壓力梯度越大,煤與瓦斯突出事故就越易發(fā)生。圖7(a)、圖8(a)中斷層帶內(nèi)瓦斯壓力等值線沿?cái)鄬觾A向分布,而圖7(b)、圖8(b)中斷層帶內(nèi)瓦斯壓力等值線沿垂直于斷層走向分布。由圖9 可以看出,從掘進(jìn)暴露面至斷層帶左側(cè)損傷區(qū)域范圍內(nèi),考慮均質(zhì)斷層的煤層內(nèi)裂隙及煤基質(zhì)瓦斯壓力大于考慮非均質(zhì)斷層的煤層內(nèi)的裂隙及煤基質(zhì)瓦斯壓力,但進(jìn)入斷層區(qū)域后兩者的瓦斯壓力大小值將迅速實(shí)現(xiàn)反轉(zhuǎn),這是由于斷層核附近高應(yīng)力集中區(qū)的存在,使得斷層核成為1 個(gè)阻礙瓦斯逃逸的封閉空間,在斷層核區(qū)域,由于斷層核的封閉作用,使斷層核內(nèi)煤巖體中儲(chǔ)存的瓦斯、后方損傷區(qū)域以及煤體中的瓦斯向采掘工作面中的擴(kuò)散受到嚴(yán)重阻礙,積聚較高的瓦斯內(nèi)能,從而產(chǎn)生壓力梯度突增的現(xiàn)象,使得在斷層帶內(nèi)非均質(zhì)斷層的瓦斯壓力梯度大于均質(zhì)斷層的瓦斯壓力梯度。在斷層帶區(qū)域范圍內(nèi),非均質(zhì)斷層帶內(nèi)的裂隙瓦斯壓力從1.09 MPa 增加到1.78 MPa,煤基質(zhì)瓦斯壓力從1.11 MPa 增加到1.78 MPa,斷層區(qū)域內(nèi)整體瓦斯壓力梯度為0.043、0.042 MPa/m。且瓦斯壓力梯度在斷層帶內(nèi)呈現(xiàn)由斷層損傷區(qū)先降低至斷層核區(qū)域,而后突增,在緩步降低的趨勢(shì)。均質(zhì)斷層帶內(nèi)的裂隙瓦斯壓力從1.18 MPa 增加到1.72 MPa,煤基質(zhì)瓦斯壓力從1.2 MPa 增加到1.73 MPa,斷層區(qū)域整體瓦斯壓力梯度為0.033 75、0.033 MPa/m,小于非均質(zhì)斷層帶內(nèi)的瓦斯壓力梯度,且瓦斯壓力梯度呈現(xiàn)持續(xù)降低的趨勢(shì)。

    煤層瓦斯壓力梯度的大小取決于瓦斯的擴(kuò)散速度和煤層的滲透性,所以在含斷層的煤層中,斷層影響區(qū)內(nèi)往往會(huì)出現(xiàn)較大的瓦斯壓力梯度。同時(shí)斷層區(qū)域由于斷層核的封閉性能,導(dǎo)致透氣性又較低,因此在巷道掘進(jìn)過程中,煤與瓦斯突出事故經(jīng)常發(fā)生在斷層的影響區(qū)域內(nèi)。為了減少煤與瓦斯突出事故的發(fā)生,應(yīng)盡可能對(duì)斷層內(nèi)部滲透率及力學(xué)性質(zhì)的各向異性進(jìn)行深入研究,更好地掌握斷層對(duì)煤體采動(dòng)破壞及瓦斯運(yùn)移影響變化規(guī)律,并采取一系列有效的措施增加煤層的滲透性,降低煤層的瓦斯壓力梯度。

    4 結(jié) 論

    1)以煤巷掘進(jìn)為背景,結(jié)合受掘進(jìn)應(yīng)力影響下的滲透率、孔隙率動(dòng)態(tài)方程、含瓦斯煤流固動(dòng)態(tài)耦合模型以及基于三角函數(shù)建立的非均質(zhì)斷層帶橫向各向同性滲透率模型,植入COMSOL 數(shù)值模擬軟件進(jìn)行結(jié)算,得出煤巷暴露面前方煤體滲透率和煤基質(zhì)、裂隙瓦斯壓力的演化規(guī)律。

    2)根據(jù)模型結(jié)算結(jié)果,煤層滲透率的變化主要影響因素是采掘應(yīng)力及煤體塑性變形。在應(yīng)力卸壓帶、應(yīng)力集中帶以及斷層帶影響區(qū)域滲透率波動(dòng)較大,斷層核處滲透率降至煤層初始滲透率的1/7,達(dá)到最小值為6.438 3×10-18m2;煤基質(zhì)和裂隙瓦斯壓力從暴露面大氣壓沿煤層走向呈現(xiàn)先增加后穩(wěn)定的趨勢(shì)。瓦斯壓力梯度在暴露面達(dá)到峰值,應(yīng)力集中帶瓦斯壓力梯度高于卸壓帶,且滲透率較低,發(fā)生煤與瓦斯突出危險(xiǎn)強(qiáng)。

    3)采用非均質(zhì)橫向各向同性滲透率模型對(duì)斷層滲透特性進(jìn)行描述時(shí),斷層帶內(nèi)裂隙及煤基質(zhì)瓦斯壓力梯度分別為0.043、0.042 MPa/m,明顯高于用均質(zhì)滲透率模型對(duì)斷層滲透特性描述時(shí)的瓦斯壓力梯度,更容易導(dǎo)致煤與瓦斯突出等動(dòng)力災(zāi)害的發(fā)生。斷層影響區(qū)瓦斯壓力梯度大,透氣性低,進(jìn)一步驗(yàn)證了在巷道掘進(jìn)過程中,斷層影響區(qū)是最易發(fā)生煤與瓦斯突出事故的區(qū)域這一結(jié)論。

    猜你喜歡
    壓力梯度均質(zhì)煤體
    注熱井周圍煤體蠕變過程的滲透率變化規(guī)律模擬研究
    正交試驗(yàn)下煤體滲透性影響因素評(píng)價(jià)
    煤炭工程(2019年4期)2019-05-05 02:30:16
    以“中央廚房”為突破口探索時(shí)政報(bào)道的融煤體之路——以浙江之聲為例
    壓力梯度在油田開發(fā)中的應(yīng)用探討
    Orlicz對(duì)偶混合均質(zhì)積分
    疊加原理不能求解含啟動(dòng)壓力梯度滲流方程
    非均質(zhì)巖心調(diào)堵結(jié)合技術(shù)室內(nèi)實(shí)驗(yàn)
    致密砂巖啟動(dòng)壓力梯度數(shù)值的影響因素
    斷塊油氣田(2014年5期)2014-03-11 15:33:45
    深部區(qū)域塊段煤體瓦斯基本參數(shù)測(cè)定技術(shù)研究
    河南科技(2014年3期)2014-02-27 14:06:05
    汽油機(jī)均質(zhì)充氣壓縮點(diǎn)火燃燒過程的混合氣形成
    日韩精品免费视频一区二区三区| 一本大道久久a久久精品| 亚洲精品日本国产第一区| 天天躁夜夜躁狠狠躁躁| 波野结衣二区三区在线| 精品一区二区免费观看| 国精品久久久久久国模美| 国产日韩欧美视频二区| a级片在线免费高清观看视频| 女人久久www免费人成看片| 超色免费av| 2018国产大陆天天弄谢| 欧美 亚洲 国产 日韩一| 天天躁夜夜躁狠狠躁躁| 国产精品一区二区在线观看99| 视频在线观看一区二区三区| svipshipincom国产片| 高清欧美精品videossex| 男人舔女人的私密视频| 国产成人精品在线电影| 精品一区二区三区四区五区乱码 | 日本色播在线视频| 一级毛片电影观看| 在线免费观看不下载黄p国产| 欧美 日韩 精品 国产| 久久亚洲国产成人精品v| 国产精品久久久久久人妻精品电影 | 欧美精品一区二区大全| 尾随美女入室| 国产精品一国产av| xxxhd国产人妻xxx| 男人爽女人下面视频在线观看| 成人国语在线视频| 美女国产高潮福利片在线看| 91老司机精品| 香蕉国产在线看| 亚洲精品美女久久久久99蜜臀 | 国产精品一二三区在线看| 国产色婷婷99| 欧美xxⅹ黑人| 少妇的丰满在线观看| 亚洲欧美一区二区三区久久| 高清在线视频一区二区三区| 黄色视频不卡| 操美女的视频在线观看| 99久久99久久久精品蜜桃| 免费黄网站久久成人精品| 成年人免费黄色播放视频| 青草久久国产| 黑人猛操日本美女一级片| av卡一久久| 妹子高潮喷水视频| a级片在线免费高清观看视频| 日韩,欧美,国产一区二区三区| 交换朋友夫妻互换小说| 夫妻午夜视频| 国产精品免费视频内射| 男女床上黄色一级片免费看| 国产精品一国产av| 精品少妇内射三级| av在线老鸭窝| 国产男女超爽视频在线观看| 最近最新中文字幕免费大全7| 欧美国产精品一级二级三级| 曰老女人黄片| 热99久久久久精品小说推荐| 天天操日日干夜夜撸| 亚洲国产av新网站| 亚洲免费av在线视频| 永久免费av网站大全| 波野结衣二区三区在线| 国产av精品麻豆| 亚洲欧美一区二区三区黑人| 国产成人精品在线电影| 18禁国产床啪视频网站| 免费看不卡的av| 欧美精品一区二区大全| 又粗又硬又长又爽又黄的视频| 天美传媒精品一区二区| 女的被弄到高潮叫床怎么办| 亚洲一区二区三区欧美精品| 亚洲成色77777| 另类亚洲欧美激情| 大陆偷拍与自拍| 少妇人妻 视频| 天美传媒精品一区二区| av国产久精品久网站免费入址| 中文精品一卡2卡3卡4更新| 1024视频免费在线观看| videosex国产| 亚洲成人av在线免费| 欧美激情高清一区二区三区 | 你懂的网址亚洲精品在线观看| 国产一卡二卡三卡精品 | 亚洲成av片中文字幕在线观看| 国产乱来视频区| 一边摸一边做爽爽视频免费| 成年av动漫网址| 人人妻人人澡人人爽人人夜夜| 欧美亚洲 丝袜 人妻 在线| 国产人伦9x9x在线观看| 亚洲专区中文字幕在线 | 国产1区2区3区精品| 亚洲国产成人一精品久久久| 久久精品人人爽人人爽视色| 又粗又硬又长又爽又黄的视频| 最近手机中文字幕大全| 精品亚洲成国产av| 天天操日日干夜夜撸| 国产精品三级大全| 黄色 视频免费看| 五月开心婷婷网| 日韩制服骚丝袜av| 久久久久网色| 国产精品二区激情视频| 性少妇av在线| www日本在线高清视频| 免费黄频网站在线观看国产| 欧美精品人与动牲交sv欧美| 国产精品女同一区二区软件| 亚洲国产欧美日韩在线播放| 国产av精品麻豆| 亚洲成国产人片在线观看| www日本在线高清视频| 日韩精品免费视频一区二区三区| 天天添夜夜摸| 国产麻豆69| netflix在线观看网站| av在线app专区| 亚洲国产精品一区二区三区在线| 国产欧美日韩综合在线一区二区| 人妻人人澡人人爽人人| 少妇被粗大的猛进出69影院| 美女主播在线视频| 欧美日韩国产mv在线观看视频| 热re99久久精品国产66热6| 大香蕉久久网| 最新在线观看一区二区三区 | 丝袜美足系列| 日韩一区二区视频免费看| 国产精品.久久久| 亚洲精品久久久久久婷婷小说| 一本—道久久a久久精品蜜桃钙片| 啦啦啦在线观看免费高清www| 色播在线永久视频| 久久久精品免费免费高清| 国产福利在线免费观看视频| 久久99热这里只频精品6学生| 久久精品久久久久久久性| 久久亚洲国产成人精品v| av免费观看日本| 亚洲婷婷狠狠爱综合网| 一区二区日韩欧美中文字幕| 中文字幕人妻熟女乱码| 日本欧美国产在线视频| 精品人妻熟女毛片av久久网站| 如日韩欧美国产精品一区二区三区| 国产精品女同一区二区软件| 日本91视频免费播放| 黄频高清免费视频| 综合色丁香网| 无遮挡黄片免费观看| 黄频高清免费视频| 多毛熟女@视频| 日韩中文字幕欧美一区二区 | 午夜影院在线不卡| 丝袜喷水一区| 自线自在国产av| 嫩草影视91久久| 一级爰片在线观看| 国产不卡av网站在线观看| 国产成人精品福利久久| 久久免费观看电影| 色婷婷av一区二区三区视频| xxxhd国产人妻xxx| 老司机在亚洲福利影院| 国产又爽黄色视频| 国产成人精品在线电影| 国产精品一国产av| 一级毛片电影观看| 欧美日韩精品网址| 亚洲国产毛片av蜜桃av| 国产成人精品在线电影| 日韩精品有码人妻一区| 久久久久久久精品精品| 午夜激情av网站| 亚洲精品第二区| 亚洲精品国产区一区二| 国产xxxxx性猛交| 免费看不卡的av| videosex国产| av福利片在线| 亚洲情色 制服丝袜| 一级黄片播放器| 秋霞伦理黄片| 看十八女毛片水多多多| 美女主播在线视频| 久久热在线av| 国产激情久久老熟女| 丝袜美足系列| 亚洲伊人色综图| 人成视频在线观看免费观看| 国精品久久久久久国模美| 女人高潮潮喷娇喘18禁视频| 十八禁人妻一区二区| 少妇猛男粗大的猛烈进出视频| 老汉色∧v一级毛片| 色网站视频免费| 日本午夜av视频| 久久ye,这里只有精品| 中国国产av一级| 亚洲精品日本国产第一区| 亚洲精品在线美女| 日韩精品免费视频一区二区三区| 赤兔流量卡办理| 伊人亚洲综合成人网| 日日撸夜夜添| 男男h啪啪无遮挡| 狠狠精品人妻久久久久久综合| 久久人人97超碰香蕉20202| 日本欧美视频一区| 美女脱内裤让男人舔精品视频| 18在线观看网站| 亚洲婷婷狠狠爱综合网| 狂野欧美激情性bbbbbb| 国产黄色视频一区二区在线观看| 女性被躁到高潮视频| 九色亚洲精品在线播放| 欧美成人精品欧美一级黄| 久久狼人影院| 国产在线免费精品| 国产一区二区三区av在线| av天堂久久9| 国产亚洲av高清不卡| 久久久欧美国产精品| 亚洲激情五月婷婷啪啪| 国产男女内射视频| 欧美人与性动交α欧美软件| 久久久国产欧美日韩av| 亚洲国产日韩一区二区| 男人操女人黄网站| 黄片播放在线免费| 国产无遮挡羞羞视频在线观看| 赤兔流量卡办理| 日韩不卡一区二区三区视频在线| 色婷婷久久久亚洲欧美| 亚洲国产日韩一区二区| 又粗又硬又长又爽又黄的视频| 秋霞在线观看毛片| 亚洲av日韩精品久久久久久密 | 亚洲一码二码三码区别大吗| 十分钟在线观看高清视频www| 亚洲精品乱久久久久久| 国产精品三级大全| 大片免费播放器 马上看| 久热爱精品视频在线9| 久久精品国产a三级三级三级| 欧美黑人精品巨大| 在线天堂最新版资源| 视频区图区小说| 啦啦啦视频在线资源免费观看| 99久久综合免费| 极品少妇高潮喷水抽搐| 天天添夜夜摸| 七月丁香在线播放| 咕卡用的链子| 国产淫语在线视频| 亚洲精品日韩在线中文字幕| 成人国语在线视频| 国产精品99久久99久久久不卡 | 欧美人与性动交α欧美精品济南到| 青春草亚洲视频在线观看| 黄网站色视频无遮挡免费观看| 国产精品国产av在线观看| 精品免费久久久久久久清纯 | 啦啦啦中文免费视频观看日本| 19禁男女啪啪无遮挡网站| 最近中文字幕2019免费版| 观看av在线不卡| 99九九在线精品视频| av在线老鸭窝| 人成视频在线观看免费观看| 国产日韩一区二区三区精品不卡| 欧美精品亚洲一区二区| 久久精品aⅴ一区二区三区四区| 国产免费福利视频在线观看| 日本猛色少妇xxxxx猛交久久| 男人爽女人下面视频在线观看| 免费黄频网站在线观看国产| 母亲3免费完整高清在线观看| av又黄又爽大尺度在线免费看| 我要看黄色一级片免费的| 国产午夜精品一二区理论片| 亚洲欧美精品综合一区二区三区| 亚洲三区欧美一区| 91精品三级在线观看| 亚洲av成人精品一二三区| 一本大道久久a久久精品| 午夜福利一区二区在线看| 如日韩欧美国产精品一区二区三区| 我的亚洲天堂| 日日啪夜夜爽| 欧美日韩亚洲综合一区二区三区_| 中国三级夫妇交换| 亚洲国产最新在线播放| 91精品国产国语对白视频| 丝袜在线中文字幕| 国产在线视频一区二区| 一级,二级,三级黄色视频| 欧美国产精品一级二级三级| 中文乱码字字幕精品一区二区三区| 熟妇人妻不卡中文字幕| 精品人妻熟女毛片av久久网站| 一级,二级,三级黄色视频| 国产麻豆69| 久久久久久久精品精品| 欧美亚洲日本最大视频资源| 日韩熟女老妇一区二区性免费视频| 精品免费久久久久久久清纯 | 丁香六月欧美| 亚洲精品一区蜜桃| 国产熟女欧美一区二区| 熟女少妇亚洲综合色aaa.| 麻豆精品久久久久久蜜桃| 精品酒店卫生间| xxx大片免费视频| 麻豆乱淫一区二区| videosex国产| 女人精品久久久久毛片| 国产片特级美女逼逼视频| 热99久久久久精品小说推荐| 波多野结衣一区麻豆| 99九九在线精品视频| www.自偷自拍.com| 国产免费现黄频在线看| 久久精品亚洲av国产电影网| 精品一区在线观看国产| 国产片特级美女逼逼视频| 精品人妻在线不人妻| 成人国语在线视频| 亚洲图色成人| 亚洲天堂av无毛| 日韩一区二区三区影片| 亚洲精品,欧美精品| 黄频高清免费视频| 人人澡人人妻人| 久久久欧美国产精品| 777米奇影视久久| 日日啪夜夜爽| 成人国产av品久久久| 天天操日日干夜夜撸| 欧美人与性动交α欧美软件| 精品免费久久久久久久清纯 | netflix在线观看网站| 久久鲁丝午夜福利片| 欧美久久黑人一区二区| a级片在线免费高清观看视频| 又大又黄又爽视频免费| 97精品久久久久久久久久精品| 成人亚洲欧美一区二区av| 日日摸夜夜添夜夜爱| 日本黄色日本黄色录像| av国产精品久久久久影院| 人人妻人人爽人人添夜夜欢视频| 午夜激情久久久久久久| 国产又爽黄色视频| 啦啦啦在线观看免费高清www| 美女视频免费永久观看网站| 香蕉国产在线看| 韩国高清视频一区二区三区| 国产精品无大码| 精品国产超薄肉色丝袜足j| 国产精品国产三级专区第一集| 久久精品亚洲熟妇少妇任你| 午夜福利网站1000一区二区三区| 丁香六月欧美| 久久女婷五月综合色啪小说| 好男人视频免费观看在线| 午夜福利免费观看在线| 一级毛片我不卡| 久久久久久久久久久免费av| 91精品伊人久久大香线蕉| 性高湖久久久久久久久免费观看| 人人澡人人妻人| 亚洲第一青青草原| 日本av手机在线免费观看| 亚洲国产精品999| 在线观看国产h片| 黄片小视频在线播放| 成年女人毛片免费观看观看9 | 成年人午夜在线观看视频| 国产毛片在线视频| 亚洲一级一片aⅴ在线观看| 一本久久精品| 成年动漫av网址| 99香蕉大伊视频| 飞空精品影院首页| 亚洲av成人精品一二三区| 自线自在国产av| 视频区图区小说| 伊人久久大香线蕉亚洲五| bbb黄色大片| 日韩制服丝袜自拍偷拍| 亚洲av男天堂| 精品人妻一区二区三区麻豆| 亚洲熟女精品中文字幕| 国产在视频线精品| 天天影视国产精品| 91国产中文字幕| 亚洲综合精品二区| 久久精品熟女亚洲av麻豆精品| 色精品久久人妻99蜜桃| 久久久久精品性色| 亚洲av国产av综合av卡| 国产亚洲一区二区精品| 国产精品无大码| 欧美 日韩 精品 国产| 最黄视频免费看| 999久久久国产精品视频| 国产欧美日韩综合在线一区二区| 亚洲精品在线美女| 伊人亚洲综合成人网| 好男人视频免费观看在线| 国产又爽黄色视频| 日韩av在线免费看完整版不卡| 性少妇av在线| 午夜激情av网站| 久久亚洲国产成人精品v| 99久久精品国产亚洲精品| 丝袜在线中文字幕| 又大又黄又爽视频免费| 九九爱精品视频在线观看| 免费不卡黄色视频| 日本欧美视频一区| 欧美另类一区| 国产在线视频一区二区| 在线观看免费高清a一片| 亚洲人成77777在线视频| 美国免费a级毛片| 制服丝袜香蕉在线| 少妇人妻精品综合一区二区| 香蕉国产在线看| 97人妻天天添夜夜摸| 色网站视频免费| 久久久久人妻精品一区果冻| 九九爱精品视频在线观看| 久久人妻熟女aⅴ| 我的亚洲天堂| 无遮挡黄片免费观看| 国产黄频视频在线观看| 亚洲欧美清纯卡通| 精品人妻熟女毛片av久久网站| 日韩成人av中文字幕在线观看| 97精品久久久久久久久久精品| 亚洲欧美精品自产自拍| 久久久久精品性色| 大香蕉久久成人网| 男人舔女人的私密视频| 成年女人毛片免费观看观看9 | 美女国产高潮福利片在线看| 欧美亚洲 丝袜 人妻 在线| 如何舔出高潮| 最近最新中文字幕免费大全7| 在线亚洲精品国产二区图片欧美| 观看av在线不卡| 极品少妇高潮喷水抽搐| 亚洲av福利一区| 80岁老熟妇乱子伦牲交| 久久久久久久精品精品| 波野结衣二区三区在线| 久久精品亚洲熟妇少妇任你| 婷婷色麻豆天堂久久| 亚洲av综合色区一区| 成人国产麻豆网| 亚洲国产毛片av蜜桃av| 国产激情久久老熟女| 黄片小视频在线播放| 啦啦啦在线免费观看视频4| 精品国产一区二区三区久久久樱花| 青青草视频在线视频观看| 亚洲精品久久成人aⅴ小说| 久久久久精品久久久久真实原创| 女人爽到高潮嗷嗷叫在线视频| 男女床上黄色一级片免费看| 成人国语在线视频| 国产亚洲一区二区精品| 亚洲av国产av综合av卡| 久久久久网色| 在线天堂中文资源库| 国产亚洲欧美精品永久| 熟妇人妻不卡中文字幕| 交换朋友夫妻互换小说| 亚洲欧美清纯卡通| 国产又色又爽无遮挡免| 久久久久精品国产欧美久久久 | 18在线观看网站| 人人澡人人妻人| 日韩av在线免费看完整版不卡| 精品一区二区免费观看| 999久久久国产精品视频| 亚洲精品久久久久久婷婷小说| 少妇的丰满在线观看| 波多野结衣av一区二区av| 少妇的丰满在线观看| 巨乳人妻的诱惑在线观看| 各种免费的搞黄视频| 丝袜喷水一区| 亚洲色图 男人天堂 中文字幕| 亚洲精品日本国产第一区| 亚洲av中文av极速乱| 飞空精品影院首页| 国产成人欧美在线观看 | 热99久久久久精品小说推荐| 在线观看免费高清a一片| 亚洲国产精品一区二区三区在线| 老司机亚洲免费影院| 18在线观看网站| 天天躁夜夜躁狠狠躁躁| 国产激情久久老熟女| 亚洲在久久综合| 国产在视频线精品| 国产极品天堂在线| 日日啪夜夜爽| 久久久久精品国产欧美久久久 | 成年动漫av网址| 日本午夜av视频| av电影中文网址| 伦理电影大哥的女人| 欧美精品av麻豆av| 精品亚洲乱码少妇综合久久| 久久99精品国语久久久| 人妻人人澡人人爽人人| 久热这里只有精品99| 日本色播在线视频| 午夜福利影视在线免费观看| 日韩欧美一区视频在线观看| 婷婷色综合大香蕉| 亚洲欧美中文字幕日韩二区| 亚洲,欧美,日韩| 成年女人毛片免费观看观看9 | 免费高清在线观看日韩| 精品一区二区三区av网在线观看 | 中文字幕色久视频| 久久av网站| 久久精品国产综合久久久| 国产黄色免费在线视频| 久久精品国产a三级三级三级| 精品久久久久久电影网| 久久人人97超碰香蕉20202| 熟女少妇亚洲综合色aaa.| 成年人免费黄色播放视频| 王馨瑶露胸无遮挡在线观看| 十八禁人妻一区二区| 极品少妇高潮喷水抽搐| 国产日韩一区二区三区精品不卡| 久久久久精品人妻al黑| 一边摸一边抽搐一进一出视频| 国产极品天堂在线| 国产日韩欧美亚洲二区| 亚洲情色 制服丝袜| 中文字幕人妻丝袜制服| 亚洲av综合色区一区| 中文精品一卡2卡3卡4更新| 日本黄色日本黄色录像| 老司机亚洲免费影院| 99香蕉大伊视频| 精品一区二区三卡| 亚洲,一卡二卡三卡| 一级片免费观看大全| www.自偷自拍.com| www.熟女人妻精品国产| 国产av精品麻豆| 丰满饥渴人妻一区二区三| 精品久久蜜臀av无| 成年人免费黄色播放视频| 欧美日韩精品网址| 菩萨蛮人人尽说江南好唐韦庄| 精品一区二区三卡| 91精品国产国语对白视频| 老鸭窝网址在线观看| av线在线观看网站| 成人漫画全彩无遮挡| 狠狠婷婷综合久久久久久88av| 久热爱精品视频在线9| 色婷婷av一区二区三区视频| a级片在线免费高清观看视频| 啦啦啦中文免费视频观看日本| 免费高清在线观看日韩| 性色av一级| 国产精品人妻久久久影院| 国产一区二区三区综合在线观看| 我的亚洲天堂| 久久久久久久久久久免费av| 中文字幕精品免费在线观看视频| 男女床上黄色一级片免费看| 9色porny在线观看| 亚洲精品美女久久久久99蜜臀 | 亚洲国产欧美一区二区综合| 不卡av一区二区三区| 久久人妻熟女aⅴ| 超碰成人久久| 男的添女的下面高潮视频| 国产片特级美女逼逼视频| 日韩大码丰满熟妇| 欧美亚洲日本最大视频资源| 国产麻豆69| 欧美97在线视频| 国产人伦9x9x在线观看| 亚洲久久久国产精品| 国产成人a∨麻豆精品| 免费观看av网站的网址| 国产精品久久久久久精品电影小说| 久久人妻熟女aⅴ| 欧美日韩精品网址| 两性夫妻黄色片|