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

    應(yīng)力恢復(fù)對采動裂隙巖體滲透性演化的影響

    2014-06-07 05:55:12王文學(xué)隋旺華董青紅
    煤炭學(xué)報(bào) 2014年6期
    關(guān)鍵詞:導(dǎo)水滲透性覆巖

    王文學(xué),隋旺華,董青紅

    (1.中國礦業(yè)大學(xué)資源與地球科學(xué)學(xué)院,江蘇徐州 221008;2.中國礦業(yè)大學(xué)深部巖土力學(xué)與地下工程國家重點(diǎn)實(shí)驗(yàn)室,江蘇徐州 221008)

    應(yīng)力恢復(fù)對采動裂隙巖體滲透性演化的影響

    王文學(xué)1,2,隋旺華1,2,董青紅1,2

    (1.中國礦業(yè)大學(xué)資源與地球科學(xué)學(xué)院,江蘇徐州 221008;2.中國礦業(yè)大學(xué)深部巖土力學(xué)與地下工程國家重點(diǎn)實(shí)驗(yàn)室,江蘇徐州 221008)

    研究了采動巖體應(yīng)力恢復(fù)對斷裂帶滲透性的影響并將其應(yīng)用于分層開采和殘煤復(fù)采導(dǎo)水?dāng)嗔褞У挠?jì)算。通過試驗(yàn)對側(cè)限作用下裂隙巖體在軸向分級加載條件下滲流變化規(guī)律進(jìn)行了研究,以模擬應(yīng)力恢復(fù)作用對裂隙巖體滲透性的影響;依據(jù)實(shí)測的采動裂隙巖體體積膨脹系數(shù)的對數(shù)分布規(guī)律,對采動裂隙巖體滲透性分布規(guī)律進(jìn)行了描述,計(jì)算了應(yīng)力恢復(fù)作用下,裂隙巖體內(nèi)裂隙開度及滲透性演化規(guī)律。試驗(yàn)結(jié)果表明:軸向加載條件下裂隙模塊開度及滲透性均顯著下降,加載初期應(yīng)力較小但其對裂隙模塊的開度及滲透性影響更顯著,后期應(yīng)力增大但其對裂隙試塊的滲透性影響程度降低,在20 MPa壓力作用下裂隙很難完全閉合;理論計(jì)算結(jié)果表明:應(yīng)力恢復(fù)作用下采空區(qū)上覆裂隙巖體內(nèi)下部大開度裂隙滲透性有降低但仍為較好的導(dǎo)水通道,上部裂隙開度較小層段滲透性也有明顯降低。根據(jù)水利水電工程地質(zhì)勘查規(guī)范將滲透系數(shù)小于0.01 m/d的層段視為不透水層,采空區(qū)上覆導(dǎo)水?dāng)嗔褞ё畲蟾叨仍趹?yīng)力恢復(fù)后較初始實(shí)測最大值減小了約26%。依據(jù)研究結(jié)果對傳統(tǒng)的分層開采覆巖導(dǎo)水?dāng)嗔褞Ц叨扔?jì)算公式進(jìn)行了修正,采用修正公式計(jì)算下分層煤或殘煤開采防水(砂)煤柱的留設(shè)方案更為合理。

    煤層復(fù)采;裂隙;應(yīng)力恢復(fù);滲透系數(shù);導(dǎo)水?dāng)嗔褞?/p>

    地下煤層的采動打破圍巖原巖應(yīng)力的平衡狀態(tài),引起圍巖的應(yīng)力重分布、變形、破裂運(yùn)動,這部分巖體稱為采動巖體[1]。采動巖體的應(yīng)力釋放、垮落、破裂等運(yùn)動導(dǎo)致其內(nèi)部裂隙發(fā)育,在采空區(qū)上覆巖體內(nèi)形成垮落帶、斷裂帶及整體彎曲下沉帶。隨著采動時(shí)間、距離的延長,采空區(qū)上覆裂隙巖體內(nèi)應(yīng)力將逐漸恢復(fù),并隨時(shí)間延續(xù),巖體內(nèi)的自重應(yīng)力將逐漸趨向于開采前的原巖應(yīng)力的垂直分量。本文中“應(yīng)力恢復(fù)”是指采動巖體中的重分布應(yīng)力在開采結(jié)束后隨著時(shí)間逐步恢復(fù)到上覆地層的自重應(yīng)力的過程。應(yīng)力恢復(fù)將對裂隙開度有壓縮減小的作用,從而影響到采動巖體的滲透性的變化。

    現(xiàn)階段對采空區(qū)應(yīng)力恢復(fù)規(guī)律的研究是采動巖體力學(xué)研究的薄弱環(huán)節(jié),針對應(yīng)力恢復(fù)對裂隙巖體滲透性演化規(guī)律的研究鮮見報(bào)道。開采覆巖裂隙巖體的后期滲透性演化規(guī)律對厚煤層的下分層開采、殘煤的復(fù)采等覆巖斷裂帶的再次發(fā)育演化規(guī)律有著重要的影響,此外采后采空區(qū)應(yīng)力恢復(fù)過程中裂隙的壓縮減小對突水潰砂的發(fā)生也有著抑制作用。本文首先對在應(yīng)力作用下裂隙巖體滲透性的演化規(guī)律進(jìn)行了室內(nèi)試驗(yàn)研究,在此基礎(chǔ)上,結(jié)合山東某礦采后15 a覆巖裂隙探測資料,對采空區(qū)應(yīng)力恢復(fù)對裂隙巖體滲透性影響的演化進(jìn)行理論分析,該研究結(jié)果將為水體下煤層分層開采及殘煤的復(fù)采的評價(jià)提供參考。

    1 應(yīng)力恢復(fù)對裂隙巖體滲透性作用的試驗(yàn)研究

    采動覆巖破壞裂隙發(fā)育多以大傾角為主,采用鉆孔電視探測統(tǒng)計(jì)分析傾角>45°的裂隙占裂隙發(fā)育總量的80%以上[8],且裂隙開度及數(shù)量分布隨距煤層距離增大呈減小趨勢。裂隙巖體在應(yīng)力恢復(fù)至原巖應(yīng)力的過程中,在泊松效應(yīng)的條件下會對覆巖裂隙產(chǎn)生壓縮閉合作用。通過在水泥砂漿內(nèi)預(yù)制裂隙,制備裂隙模塊,在軸向加載條件下測試滲透性的變化,可以反映裂隙巖體在應(yīng)力恢復(fù)過程的滲透性變化規(guī)律。

    1.1 試驗(yàn)?zāi)P?/p>

    滲透裝置為自行設(shè)計(jì)加工,滲透裝置外壁厚20 mm,材質(zhì)選用普通45號鋼進(jìn)行全淬火處理,以提高硬度防止側(cè)向變形。滲透試驗(yàn)?zāi)P腿鐖D1所示,模型內(nèi)直徑為100 mm、高240 mm,其主要有上部加壓活塞、上下濾水板、底蓋及滲透缸體構(gòu)成,底蓋留取進(jìn)水口并通過裂隙模塊后由加壓活塞出水口排出,軸向荷載通過伺服加載系統(tǒng)完成。最大滲透壓力為0.5 MPa,由穩(wěn)壓水源提供。裂隙模塊采用水泥砂漿按不同配比制備,模塊直徑為100 mm、高200 mm,其內(nèi)部裂隙長度80 mm、高200 mm,裂隙兩側(cè)留取10 mm完整試塊防止向側(cè)壁滲流,裂隙試塊形狀如圖2所示。裂隙制備過程中在水泥砂漿中預(yù)設(shè)不同厚度(代表裂隙開度)、寬80 mm、高200 mm的鋼板,待試塊有一定強(qiáng)度后將鋼板拔出,并將磨具去掉,放在水中養(yǎng)護(hù)21 d,再將試塊放入滲透模型中進(jìn)行試驗(yàn)。試驗(yàn)過程中對進(jìn)水口壓力、出水口壓力、滲流量、軸向載荷以及試塊的應(yīng)力-應(yīng)變等進(jìn)行監(jiān)測。

    圖1 滲透試驗(yàn)裝置示意Fig.1 Schematic graph of the set up for permeability test

    圖2 裂隙試塊示意Fig.2 Schematic of the block with a fracture

    本文選擇裂隙傾角為90°,開度為0.80 mm。完整試塊平均單軸抗壓強(qiáng)度為9.53 MPa的裂隙試塊,在軸向加載過程中的滲透性變化規(guī)律進(jìn)行分析,完整試塊的單軸抗壓試驗(yàn)的應(yīng)力位移曲線如圖3所示。對裂隙試塊在側(cè)限作用下分級施加軸向荷載,直至軸向壓力為20 MPa,每次施加軸向荷載后測定裂隙巖塊的滲透系數(shù)。

    1.2 試驗(yàn)結(jié)果及分析

    圖3 水泥砂漿完整試塊的單軸抗壓強(qiáng)度試驗(yàn)應(yīng)力-位移曲線Fig.3 The curves of stress via displacementof the cemented sand in UCS testing

    試驗(yàn)中對多級荷載作用下不同滲透壓差的滲流量進(jìn)行監(jiān)測,圖4給出了在軸向應(yīng)力小于8 MPa時(shí),不同荷載條件下的滲流量與滲透壓差的關(guān)系,可見二者呈較好的線性關(guān)系。當(dāng)軸向應(yīng)力大于8 MPa時(shí),由于在軸向荷載作用下,裂隙開度變得很小,低滲透壓差下滲流量較小,主要對高滲透壓差條件下滲流量進(jìn)行了觀測。圖4可以看出同一級荷載作用下滲流量隨滲透壓差的增加呈線性增加,且隨著軸向加載強(qiáng)度的增加同一滲透壓差作用下滲流量在逐漸減小。圖5為滲透壓差分別為100~400 kPa時(shí),滲流量與軸向應(yīng)力間的關(guān)系,可以看出裂隙試塊隨著軸向應(yīng)力的增加滲流量逐漸減小,采用y=A e-Bx+C的形式擬合,R2均大于0.96,說明裂隙試塊滲流量與軸向應(yīng)力呈很好的負(fù)指數(shù)關(guān)系。

    圖4 不同等級荷載作用下滲流量與滲透壓差的關(guān)系Fig.4 The relationship between the seepage discharge and seepage pressure under different axial loads

    圖5 滲流量與軸向壓力的關(guān)系Fig.5 The relationship between the seepage and the axial loads

    穩(wěn)定流裂隙滲流符合立方定理,如式(1)所示,根據(jù)式(1)可以推導(dǎo)計(jì)算裂隙的水力開度的變化規(guī)律,如圖6(a)所示。由圖6(a)可以看出裂隙開度隨著軸向荷載的增加逐漸減小,初始加載階段裂隙開度的減小速率較大。對其進(jìn)行擬合發(fā)現(xiàn),裂隙開度隨軸向荷載增加符合負(fù)指數(shù)變化的規(guī)律。

    圖6 軸向荷載與裂隙水力開度、裂隙滲透系數(shù)的關(guān)系Fig.6 The relationship between the fracture aperture,thepermeability coefficientand axial load

    Witherspoon P A等[8]指出裂隙滲透系數(shù)與裂隙開度的關(guān)系滿足

    根據(jù)式(2)及不同荷載對應(yīng)的裂隙水力開度可對其滲透系數(shù)進(jìn)行計(jì)算,應(yīng)力-滲流耦合關(guān)系如圖5所示,可以看出側(cè)限作用下裂隙試塊的滲透系數(shù)變化與軸向應(yīng)力呈負(fù)指數(shù)函數(shù)關(guān)系。裂隙滲透系數(shù)在軸向加載前期迅速減小,后期逐漸變緩。該結(jié)果表明軸向荷載對裂隙試塊的滲透性有著重要的影響,裂隙試塊在應(yīng)力較小的情況下,滲透性隨應(yīng)力的變化敏感;而當(dāng)應(yīng)力較大時(shí),裂隙已經(jīng)閉合到一定程度,再同樣增大軸向應(yīng)力對裂隙試塊開度變化及滲透性的影響逐漸減弱,由此可以推斷,采動裂隙巖體應(yīng)力的恢復(fù)對裂隙巖體內(nèi)開度及滲透性有著重要的影響,應(yīng)力恢復(fù)初始階段對裂隙巖體的影響更為顯著,將導(dǎo)致其滲透性顯著降低,應(yīng)力恢復(fù)后期對裂隙巖體滲透性的進(jìn)一步改變較弱。

    2 開采覆巖裂隙簡化模型及應(yīng)用

    2.1 研究區(qū)概述

    山東省某礦區(qū)3號煤層厚度為7.9~9.2 m,煤層傾角3.5°~8.0°,含煤地層被第四系直接覆蓋,第四系厚度在178.4~196.2 m,平均186.18 m?;鶐r上被2~3層中粗砂覆蓋,單位涌水量0.005 473~0.045 890 L/(s·m),屬弱含水層。煤層采高為2.2 m,采用全垮落長壁開采,頂板垮落充分。煤層開采覆巖裂隙充分發(fā)育時(shí)實(shí)測垮落帶、斷裂帶最大發(fā)育高度分別為7.3和20.8 m。開采區(qū)域裂隙巖體經(jīng)過15 a的下沉變形后,再次采用鉆孔電視、沖洗液漏失量、RQD等方法對該區(qū)域采空區(qū)上覆垮落帶和斷裂帶巖體探測,結(jié)果表明由于上覆地層的沉降、裂隙內(nèi)充填物的充填、破碎巖體遇水膨脹等因素的影響,垮落帶巖體已經(jīng)壓密壓實(shí),裂隙巖體內(nèi)裂隙開度明顯減小甚至閉合,觀測研究表明最大導(dǎo)水?dāng)嗔褞Ц叨葴p小近40%,采動巖體內(nèi)自重應(yīng)力已經(jīng)基本恢復(fù)[9]。

    2.2 斷裂帶裂隙及滲透性分布規(guī)律

    現(xiàn)場實(shí)際觀測表明采空區(qū)垂直方向垮落帶及斷裂帶巖體體積膨脹系數(shù)服從對數(shù)分布[10]。Shao等[11]指出垮落帶破碎巖體體積膨脹系數(shù)垂直方向變化不大,可看為定值,斷裂帶巖體體積膨脹系數(shù)垂直方向服從對數(shù)遞減,斷裂帶與彎曲下沉帶接觸位置體積膨脹系數(shù)為1,但未考慮斷裂帶形成過程中垮落帶壓縮變形量?;诖?本文對斷裂帶形成過程垮落帶的壓縮變形量進(jìn)行考慮,指出采空區(qū)上部裂隙巖體體積膨脹系數(shù)分布規(guī)律如式(3)所示。

    地表產(chǎn)生沉降前,開采煤層空間完全由垮落帶和斷裂帶空隙和裂隙空間組成,如下:

    式中,bz為垂直方向位置z處體積膨脹系數(shù);b為垮落帶初始體積膨脹系數(shù);b′為垮落帶壓縮Δh后體積膨脹系數(shù);hc為垮落帶高度;hf為斷裂帶高度;ck為垂直方向體積膨脹系數(shù)衰減系數(shù);=1為垮落帶和彎曲帶接觸位置體積膨脹系數(shù)。

    均勻孔隙介質(zhì)的孔隙率可由式(7)求得。

    根據(jù)Blake-Kozeny方程式(8),可對孔隙介質(zhì)滲透率k進(jìn)行計(jì)算:

    式中,φ為垮落帶巖塊直徑;k0為垮落帶初始滲透率; n0為初始孔隙率;k為滲透率。

    將式(3)和(7)代入式(8)得到不同高度裂隙巖體滲透系數(shù)與體積膨脹系數(shù)的關(guān)系,如下:

    孔隙介質(zhì)滲透率可以按式(10)可轉(zhuǎn)換為滲透系數(shù),即

    式中,K為滲透系數(shù);γ為流體容重。

    垮落帶高度與初始體積膨脹系數(shù)關(guān)系如式(11)所示:

    孔隙與裂隙介質(zhì)間可以通過對應(yīng)寬度內(nèi)滲流量相等進(jìn)行轉(zhuǎn)換[12],即Qp=Qf,轉(zhuǎn)換過程中將裂隙間完整巖塊視為不透水,孔隙介質(zhì)視為均勻滲透,如圖7所示。如圖7所示假定在垂直方向孔隙介質(zhì)滲透系數(shù)為Kp,在其寬度為s范圍內(nèi)的滲流量Qp可通過式(12)表示:

    孔隙介質(zhì)寬度為s范圍的垂直滲流量,可以通過單裂隙開度為d的滲流量Qf表示,如果單裂隙與水平方向夾角為α則單裂隙垂直方向滲流量可以通過式(13)表示:

    式中,Kf為裂隙介質(zhì)滲流系數(shù);I為水力梯度;s為孔隙介質(zhì)單位寬度即轉(zhuǎn)化后裂隙間距;d為裂隙的開度。

    式(12)與(13)對應(yīng)相等可以得到

    圖7 孔隙介質(zhì)與裂隙介質(zhì)轉(zhuǎn)換關(guān)系Fig.7 The transformational relationship between theporous and fractured medium

    聯(lián)立式(2)和(14)可以求得裂隙介質(zhì)的裂隙開度:

    根據(jù)式(11)計(jì)算垮落帶初始體積膨脹系數(shù)為1.3,將垮落帶巖塊直徑取值為0.5 m,按照式(8)計(jì)算其初始滲透率為3.42×10-5m2,對應(yīng)滲透系數(shù)為342 m/s;γ為1×104N/m3,μ為1×10-3N·s/m,將上述參數(shù)代入式(9)和(10)計(jì)算分析,導(dǎo)水?dāng)嗔褞?nèi)滲透系數(shù)分布規(guī)律如圖8(a)所示。

    圖8 應(yīng)力恢復(fù)前、后導(dǎo)水?dāng)嗔褞?nèi)裂隙開度及滲透性分布規(guī)律Fig.8 The aperture and permeability distribution of the fractures in the water flow fractured zone before and after the stress re-establishment

    假定裂隙傾角為90°,間距為0.4 m,根據(jù)式(15)計(jì)算斷裂帶內(nèi)裂隙開度分布規(guī)律如圖8(a)所示。由圖8(a)可以看出,導(dǎo)水?dāng)嗔褞?nèi)裂隙的開度及滲透性隨著距煤層距離的增加而減小,其滲透性減小更明顯。

    2.3 采后應(yīng)力恢復(fù)對裂隙演化的作用

    巖體受力情況下,在裂隙部位可能會產(chǎn)生大的變形,裂隙甚至閉合[13]。Bandis等[14]通過大量的室內(nèi)試驗(yàn)研究提出法向應(yīng)力對裂隙閉合效應(yīng)的雙曲線計(jì)算方法如式(16)所示:

    式中,σn為法向有效應(yīng)力;dnc為應(yīng)力作用下裂隙閉合度量;dmax為最大閉合度;kni為法向應(yīng)力作用下的裂隙剛度。

    巖體裂隙剛度對其受力變形起著關(guān)鍵作用,裂隙面接觸面積、裂隙開度、裂隙粗糙度、裂隙巖塊的強(qiáng)度、裂隙間充填物的厚度、類型、物理力學(xué)性質(zhì)等均對裂隙的剛度有著重要的影響[15-17]。Gale[18]指出天然裂隙的法向應(yīng)力剛度要大于誘發(fā)裂隙剛度;不同的裂隙剛度測試手段及不同裂隙類型其剛度差別可達(dá)幾個(gè)數(shù)量級,Philippe S等[15]首次采用Schmidt錘對裂隙巖體剛度進(jìn)行現(xiàn)場測試,發(fā)現(xiàn)張開裂隙、縫合線等裂隙剛度的量級為10-1GPa/m,封閉填充裂隙剛度量級為101~102GPa/m;Fransson?.等[19]通過注水試驗(yàn)和注漿數(shù)據(jù)估算淺部隧道開采誘發(fā)裂隙剛度為2~5 GPa/m;Griffith W A等[20]對斷層帶巖體裂隙剛度建立模型分析,將裂隙間距定為0.22 m時(shí)計(jì)算裂隙剛度為28 GPa/m。

    研究區(qū)覆巖以粉砂巖、砂質(zhì)泥巖為主且風(fēng)化程度高,強(qiáng)度較低,單軸抗壓強(qiáng)度在20~40 MPa之間屬于中硬偏軟巖。采動覆巖裂隙為擾動裂隙,裂隙內(nèi)未被充填,其裂隙剛度低,本文將覆巖裂隙剛度取值為上述文獻(xiàn)中裂隙剛度實(shí)測資料最小值10-1GPa/m。

    簡化斷裂帶裂隙分布形式,假設(shè)裂隙按間距0.4 m均勻分布,傾角為90°,單裂隙受力模型如圖8(a)所示。在此,只考慮裂隙在泊松效應(yīng)的作用下產(chǎn)生的水平法向應(yīng)力的作用。

    圖9為斷裂帶裂隙應(yīng)力恢復(fù)受力分析示意,其中σ1為垂直恢復(fù)應(yīng)力即上覆地層自重應(yīng)力,σ1=γH,σn為裂隙面法向應(yīng)力,τ為裂隙面剪切應(yīng)力,σ3為水平應(yīng)力。本文中的裂隙模型σn=σ3=λσ1,λ為泊松比,采動裂隙巖體破碎程度高,巖體質(zhì)量差,根據(jù)工程巖體分級標(biāo)準(zhǔn)GB 50218—94[21]該類巖體泊松比較大,本文計(jì)算時(shí)取值為0.3。

    圖9 裂隙受力分析示意Fig.9 Schematic of the stress analysis on the fracture

    研究區(qū)煤層開采結(jié)束已有15 a的時(shí)間,上覆地層沉降移動變形已基本結(jié)束,采空區(qū)內(nèi)垂直應(yīng)力基本恢復(fù)到自重應(yīng)力的垂直分量。采空區(qū)上覆巖層按平均厚度25 m、松散層厚度平均186.2 m計(jì)算,容重分別取值為22與20 kN/m3。導(dǎo)水?dāng)嗔褞б陨系貙幼灾貞?yīng)力完全作用在采動裂隙巖體上,在斷裂帶內(nèi)部垂直方向下部裂隙承受上部裂隙巖體的垂直荷載,斷裂帶內(nèi)裂隙自上而下逐漸增加,導(dǎo)水?dāng)嗔褞细驳貙幼灾貞?yīng)力為3.82 MPa。其內(nèi)部應(yīng)力增加可采用式(17)計(jì)算:

    式中,H0為斷裂帶上覆地層厚度;HZ為距斷裂帶頂端的距離。

    應(yīng)力恢復(fù)后裂隙巖體內(nèi)裂隙開度為

    根據(jù)式(18)和(19)可以計(jì)算出應(yīng)力恢復(fù)后斷裂帶裂隙開度的分布規(guī)律,通過式(2)和(14)可對應(yīng)力恢復(fù)作用下的裂隙巖體滲透性進(jìn)行反演推算,其結(jié)果如圖8(b)所示。

    由圖8(b)可以看出,應(yīng)力恢復(fù)作用對裂隙巖體的滲透性有明顯的減小作用,裂隙開度減小。導(dǎo)水?dāng)嗔褞虏苛严堕_度雖有減小但仍較大,滲透性較好;導(dǎo)水?dāng)嗔褞喜苛严堕_度減小,滲透性明顯減低,根據(jù)水利水電工程地質(zhì)勘查規(guī)范GB 50487—2008[22]將滲透系數(shù)小于0.01 m/d斷裂帶巖層段視為不透水層,應(yīng)力恢復(fù)后導(dǎo)水?dāng)嗔褞ё畲蟾叨葹?4.3 m,較斷裂帶初始最大高度減小26.0%。該分析結(jié)果小于現(xiàn)場探測結(jié)果40%[9],這說明其它因素如裂隙巖體的遇水膨脹、裂隙內(nèi)碎屑物的充填等對裂隙巖體的后期演化也有著重要影響。

    2.4 應(yīng) 用

    室內(nèi)試驗(yàn)與理論分析研究均表明,裂隙巖體的應(yīng)力恢復(fù)會造成斷裂帶裂隙開度、滲透性的減小。應(yīng)力恢復(fù)很難使裂隙完全壓縮閉合,但對于開度較小、剛度較低的裂隙在應(yīng)力恢復(fù)作用下其滲透性降低更明顯,可能會降低至小于0.01 m/d,屬于弱或不透水,具有良好的隔水性。該層段在下分層煤層開采覆巖破壞計(jì)算時(shí)可作為隔水巖層考慮,斷裂帶發(fā)育高度的再次計(jì)算可考慮從滲透性大于0.01 m/d巖層段疊加計(jì)算分析。

    《建筑物、水體、鐵路及主要井巷煤柱留設(shè)與壓煤開采規(guī)程》[23]中水體下煤層開采導(dǎo)水?dāng)嗔褞в?jì)算公式為

    式中,hf為導(dǎo)水?dāng)嗔褞Ц叨?A,B,C為常數(shù)參數(shù);h為采厚。

    該公式在對分層開采過程中覆巖裂隙破壞高度的計(jì)算中未考慮頂分層采后應(yīng)力恢復(fù)對裂隙閉合作用的影響。實(shí)際很多煤層下分層開采在上分層采后較長時(shí)間才進(jìn)行,該階段上分層采后裂隙巖體應(yīng)力已基本恢復(fù)。

    以頂分層開采應(yīng)力恢復(fù)以后二分層開采覆巖導(dǎo)水裂隙發(fā)育高度計(jì)算為例,對該公式進(jìn)行修正。假設(shè)頂分層采后導(dǎo)水?dāng)嗔褞Э蓾B透高度在應(yīng)力恢復(fù)作用下降低率為η,該值可以根據(jù)上文應(yīng)力恢復(fù)作用下裂隙閉合演化規(guī)律確定,本文該值為26%。因此應(yīng)力恢復(fù)后頂分層導(dǎo)水?dāng)嗔褞Ц叨热缡?19)所示,為下面計(jì)算方便將修正參數(shù)C省略:

    其中,h′f為應(yīng)力恢復(fù)后開采覆巖破壞導(dǎo)水?dāng)嗔褞Ц叨?。該覆巖裂隙破壞高度對應(yīng)著相應(yīng)的煤層開采厚度h′,即

    由式(20)對應(yīng)力恢復(fù)后的導(dǎo)水?dāng)嗔褞Ц叨葘ζ涞刃簩娱_采厚度h′反演計(jì)算:

    在二分層開采覆巖裂隙破壞高度計(jì)算過程中可將h′與二分層開采厚度h疊加代入式(22)計(jì)算,即

    該計(jì)算結(jié)果可更合理的確定下分層煤開采覆巖破壞高度,為防水煤柱的留設(shè)提供更科學(xué)合理的方法。

    此外,采空區(qū)裂隙巖體的應(yīng)力恢復(fù)會降低突水潰砂事故發(fā)生的概率及突水過程的涌水量,開采初期最大導(dǎo)水?dāng)嗔褞Оl(fā)育過程中未誘發(fā)突水潰砂事故,在后期其發(fā)生的可能性也不大。

    3 結(jié) 論

    (1)對側(cè)限作用下裂隙試塊軸向分級加載滲透性變化進(jìn)行了室內(nèi)試驗(yàn)研究,結(jié)果表明,初始加載階段裂隙的開度及其滲透性減小十分顯著,后期應(yīng)力增加但對裂隙的開度及滲透性影響不明顯。側(cè)限條件下平行裂隙面垂向加載對裂隙開度及滲透性的減小作用證實(shí)了垂直應(yīng)力恢復(fù)在泊松效應(yīng)下產(chǎn)生水平應(yīng)力對巖體裂隙的閉合效應(yīng)。

    (2)根據(jù)采空區(qū)上覆裂隙巖體體積膨脹系數(shù)服從對數(shù)分布的規(guī)律,采用孔隙介質(zhì)模型對覆巖裂隙的滲透性進(jìn)行定量性描述,結(jié)果表明其滲透性隨距煤層底板距離的增加急劇降低;利用孔隙介質(zhì)與裂隙介質(zhì)間轉(zhuǎn)換關(guān)系,簡化采空區(qū)上覆裂隙巖體裂隙分布規(guī)律,結(jié)果表明斷裂帶內(nèi)裂隙開度隨距煤層底板距離增加快速減小。

    (3)采空區(qū)上覆裂隙巖體應(yīng)力恢復(fù)后,裂隙巖體的裂隙開度及滲透性均較應(yīng)力恢復(fù)前大幅度減低,下部裂隙開度仍較大且滲透性強(qiáng),上部開度較小裂隙在應(yīng)力恢復(fù)作用下其滲透性可降低至0.01 m/d以下,屬于弱或不透水層,按此分析實(shí)例中的最大導(dǎo)水?dāng)嗔褞Ц叨仍趹?yīng)力恢復(fù)后降低26.0%。裂隙巖體內(nèi)該滲透性降低至小于0.01m/d層段在下分層開采或殘煤復(fù)采過程中可作為隔水層段考慮,并據(jù)此對分層開采覆巖裂隙破壞計(jì)算公式做了修正,采用修正公式對下分層煤層開采覆巖破壞裂隙高度的計(jì)算可更合理的確定防水(砂)煤柱的留設(shè)。

    [1] 錢鳴高,石平五,許家林.礦山壓力與巖層控制[M].徐州:中國礦業(yè)大學(xué)出版社,2010.

    [2] Alan A Compoli,Timothy M Barton,Fred C Van Dyke.Gob and gate road reaction to longwallmining in bump-prone strata[R].Reportof Investigation(United states department of the interior,Bureau ofmines),1993:9445.

    [3] King H J,Whittaker B N.A review of current knowledge on roadway behaviour,especially the problems on which further information is required[A].Proceedings of the Symposium on Strata Control in Roadways[C].London:Inst.Min.Met.1971:73-87.

    [4] Choi D S,McCain D L.Design of longwall systems[J].Trans Soc.Min.Eng.AIME,1980,258:1761-1764.

    [5] Mark C.Analysis of longwall pillar stability(alps):an update.proceedings of theworkshop on coal pillarmechanics and design[R].BuMines:Information Circular,1992:78-93.

    [6] 王作宇,劉鴻泉.采空區(qū)應(yīng)力、覆巖移動規(guī)律與頂?shù)装鍘r體應(yīng)力效應(yīng)的一致性[J].煤礦開采,1993(1):40-46.

    Wang Zuoyu,Liu Hongquan.Consistent of goaf stress,overlying stra-tamovement rule and roof and floor rock mass force effect[J].Coal Mining,1993(1):40-46.

    [7] 徐連滿,潘一山,李忠華,等.深部開采覆巖應(yīng)力變化規(guī)律模擬實(shí)驗(yàn)研究[J].中國地質(zhì)災(zāi)害與防治學(xué)報(bào),2011,22(3):61-66.

    Xu Lianman,Pan Yishan,Li Zhonghua,et al.Similarity simulation experiment of deformation of deep mining-induced overburden stress [J].The Chinese Journal of Geological Hazard and Control,2011, 22(3):61-66.

    [8] Witherspoon P A,Wang Y S Y,Gale JE.Validity of cubic law for fluid flow in a deformed fracture[J].Water Resources Research 1980,16:1016-1024.

    [9] 王文學(xué),隋旺華,董青紅,等.松散層下覆巖裂隙采后閉合效應(yīng)及重復(fù)開采覆巖破壞預(yù)測[J].煤炭學(xué)報(bào),2013,38(10):1723-1729.

    Wang Wenxue,Sui Wanghua,Dong Qinghong,et al.Closure effect ofmining-induced fractures under sand aquifers and prediction of overburden failure due to re-mining[J].Journal of China Coal Society,2013,38(10):1723-1729.

    [10] 郭廣禮,繆協(xié)興,張振南.老采空區(qū)破裂巖體變形性質(zhì)研究[J].科學(xué)技術(shù)與工程,2002,2(5):44-47.

    Guo Guangli,Miao Xiexing,Zhang Zhennan.Research on ruptured rock mass deformation characteristics of Longwall goafs[J].Sci.Technol.Eng.,2002,2(5):44-47.

    [11] Shao H,Jiang S G,Wang L Y,et al.Bulking factor of the strata overlying the gob and a three-dimensional numerical simulation of the air leakage flow field[J].Mining Science and Technology, 2011,21:261-266.

    [12] Liu J,Elsworth D,Brady BH.Linking sress-dependenteffective porosity and hydraulic conductivity fields to RMR[J].International Journal of Rock Mechanics and Mining Sciences,1991,36:581-596.

    [13] Liu H H,Rutquist J,James G B.On the relationship between stress and elastic strain for porous and fractured rock[J].International Journal of Rock Mechanical&Mining Science,2009;46:289-296.

    [14] Bandis S C,Lumsden A C,Barton N R.Fundamentals of rock joint deformation[J].Int.J.Rock Mech.Min.Sci.&Geomech, 1983:249-268.

    [15] Philippe S,Aurélien B,Rodolphe C,et al.In-situ characterization of the effective elasticity of a fault zone,and its relationship to fracture spacing[J].Journal of Structural Geology,2011,33:1541-1553.

    [16] Worthington M H,Lubbe R.The scaling of fracture compliance [J].Geological Society,London,Special Publications,2007,270: 73-82.

    [17] Alireza B,Lanru J.Stress effects on permeability in a fractures rock masswith correlated fracture length and aperture[J].International Journal of Rock Mechanics and Mining Sciences,2008,45: 1320-1334.

    [18] Gale JE.The effects of fracture type(induced versus natural)on the stress-fracture closure-permeability relationship[A].Googman RE,Huze RE,editors.Proceedings of the 23rd Rock Mechanics Symposium[C].Berkeley:AIME,1982:290-298.

    [19] Fransson?,Tsang C F,Rutqvist J,et al.Estimation of deformation and stiffness of fractures close to tunnels using data from single-hole hydraulic testing and grouting[J].In International Journal of Rock Mechanics and Mining Sciences,2010,47(6):887-893.

    [20] Griffith W A,Pablo F S,David D Pollard.Influence of outcrop scale fractures on the effective stiffness of fault damage zone rocks [J].Pure appl.Geophys,2009,166:1595-1627.

    [21] GB 50218—94,工程巖體分級標(biāo)準(zhǔn)[S].

    [22] GB 50487—2008,水利水電工程地質(zhì)勘察規(guī)范[S].

    [23] 國家煤炭工業(yè)局.建筑物、水體、鐵路及主要井巷煤柱留設(shè)與壓煤開采規(guī)程[S].北京:煤炭工業(yè)出版社,2000.

    Coal Industry Bureau of the Peop le’s Republic of China.Regulations for setting coal pillar and mining under buildings,water bodies,railways and main roadways[S].Beijing:China Coal Industry Publishing House,2000.

    Influence of cover stress re-establishment on the permeability evolution of them ining-induced fractures

    WANGWen-xue1,2,SUIWang-hua1,2,DONG Qing-hong1,2
    (1.School ofResourcesand Geosciences,China University ofMining and Technology,Xuzhou 221008,China;2.State Key Laboratory forGeomechanics and Deep Underground Engineering,China University ofMining and Technology,Xuzhou 221008,China)

    This paper presented an investigation of the influence of cover stress re-establishmenton the permeability of mining-induced fractured zone and its application to the calculation of the heightofwater flow fractured zone for slices mining and re-mining of abandoned coal seams.An experiment on the hydraulic conductivity of the fractures under axial loading with lateral restraintwas conducted to simulate influence of cover stress re-establishmenton the permeability of fractured zone.Based on the in-situ measured logarithmic distribution of the bulking factor with the distance from the gob in themining-induced fractured zone,the distribution of permeability in thewater flow fractured zone has been proposed,and the evolutions of permeability and aperture of themining-induced fractures under cover stress re-establishmentwere analyzed by the fractured rockmassmechanicsmethod.The experimental results show that the apertureand the hydraulic conductivity of the fracture decrease more significantly in the earlier stage of the axial loading than that in the later stage,even though the stress is smaller;the later stage of the axial loading plays a less role on the influence of the fracture’s permeability,and the aperture is hardly compressed to be closed under an axial loading of20 MPa.The theoretical calculation results show that the permeability of the fractures with larger aperture in the lower part of the water flow fractured zone has been decreased but still has a good permeable pathway for water flow;and the permeability of the upper part of the water flow fractured zone has been decreased obviously.According to the code for engineering geological investigation ofwater resources and hydropower,the hydraulic conductivity of 0.01 m/d is determined as the standard for impermeable layer,the maximum height of the water flow fractured zone after the cover stress re-establishment is decreased by 26%compared to the original one.The formula for calculating height of water flow fractured zone ismodified based on the research results,which would provide amore suitable approach to determine the size of safety pillars formining the lower slices of coal seam and re-mining the abandoned seams.

    re-mining of coal seam;fracture;cover stress re-establishment;hydraulic conductivity;water flow fractured zone

    TD313

    A

    0253-9993(2014)06-1031-08

    王文學(xué),隋旺華,董青紅.應(yīng)力恢復(fù)對采動裂隙巖體滲透性演化的影響[J].煤炭學(xué)報(bào),2014,39(6):1031-1038.

    10.13225/j.cnki.jccs.2014.0213

    Wang Wenxue,SuiWanghua,Dong Qinghong.Influence of cover stress re-establishment on the permeability evolution of the mining-induced fractures[J].Journal of China Coal Society,2014,39(6):1031-1038.doi:10.13225/j.cnki.jccs.2014.0213

    2014-02-24 責(zé)任編輯:常 琛

    國家自然科學(xué)基金委員會與神華集團(tuán)有限公司聯(lián)合資助項(xiàng)目(51174286);國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973)資助項(xiàng)目(2013CB227903);江蘇高校優(yōu)勢學(xué)科建設(shè)工程資助項(xiàng)目(地質(zhì)資源與地質(zhì)工程)

    王文學(xué)(1985—),男,江蘇徐州人,博士研究生。E-mail:wang603698305@163.com。通訊作者:隋旺華(1964—),男,山東臨沂人,教授,博士生導(dǎo)師。E-mail:suiwanghua@cum t.edu.cn

    猜你喜歡
    導(dǎo)水滲透性覆巖
    不同固化劑摻量對濕陷性黃土強(qiáng)度和滲透性的影響
    視唱練耳課程與作曲技術(shù)理論的交叉滲透性探究
    樂府新聲(2021年4期)2022-01-15 05:50:08
    礦區(qū)開采過程中覆巖移動規(guī)律數(shù)值模擬研究*
    煤礦高強(qiáng)度長壁開采覆巖破壞充分采動及其判據(jù)
    準(zhǔn)東大井礦區(qū)巨厚煤層開采覆巖裂隙分布特征
    阿奇山1號巖體現(xiàn)場水文地質(zhì)試驗(yàn)及滲透性評價(jià)
    紡織基小口徑人工血管材料的體外水滲透性研究
    充填開采覆巖變形破壞規(guī)律研究
    中國煤炭(2016年1期)2016-05-17 06:11:36
    改性滌綸針織物導(dǎo)水透濕性能研究
    導(dǎo)水陷落柱預(yù)注漿治理
    国产成人一区二区在线| 久久精品亚洲av国产电影网| 麻豆av在线久日| 国产有黄有色有爽视频| 大话2 男鬼变身卡| 两性夫妻黄色片| 99国产精品一区二区蜜桃av | 国产成人91sexporn| 日本av手机在线免费观看| 国产片特级美女逼逼视频| 精品亚洲成国产av| 性色av乱码一区二区三区2| 一区在线观看完整版| 国产片内射在线| 亚洲国产精品成人久久小说| 国产又爽黄色视频| 性色av乱码一区二区三区2| 国产一区二区在线观看av| 精品一区二区三区av网在线观看 | 免费av中文字幕在线| 在现免费观看毛片| 国产日韩欧美亚洲二区| 国产精品一区二区在线观看99| 一区在线观看完整版| 黄片播放在线免费| 晚上一个人看的免费电影| 午夜福利,免费看| 汤姆久久久久久久影院中文字幕| 久久久国产一区二区| 99国产精品一区二区三区| 七月丁香在线播放| 91九色精品人成在线观看| 亚洲 国产 在线| 19禁男女啪啪无遮挡网站| 日本av免费视频播放| 丰满迷人的少妇在线观看| 秋霞在线观看毛片| 十八禁网站网址无遮挡| 女警被强在线播放| 午夜福利影视在线免费观看| 午夜福利视频在线观看免费| 欧美精品一区二区大全| 国产精品一国产av| 亚洲免费av在线视频| 一区二区av电影网| 1024视频免费在线观看| 婷婷色综合大香蕉| 久久性视频一级片| 精品国产国语对白av| 国产成人一区二区三区免费视频网站 | 一级片'在线观看视频| 国产伦人伦偷精品视频| 欧美av亚洲av综合av国产av| 大片电影免费在线观看免费| 女性生殖器流出的白浆| 午夜激情av网站| 亚洲欧美精品自产自拍| 国产爽快片一区二区三区| 男女国产视频网站| 久久天躁狠狠躁夜夜2o2o | 狠狠婷婷综合久久久久久88av| 韩国精品一区二区三区| 日日摸夜夜添夜夜爱| av国产精品久久久久影院| 国产精品三级大全| 亚洲欧洲精品一区二区精品久久久| 青春草亚洲视频在线观看| 国产片特级美女逼逼视频| 青春草视频在线免费观看| 成年美女黄网站色视频大全免费| 国产激情久久老熟女| 曰老女人黄片| 久久天堂一区二区三区四区| 婷婷色综合大香蕉| 午夜福利影视在线免费观看| 国产av一区二区精品久久| 天天躁夜夜躁狠狠久久av| 免费av中文字幕在线| 亚洲精品第二区| 亚洲视频免费观看视频| 欧美人与性动交α欧美精品济南到| 日日夜夜操网爽| 亚洲精品国产av蜜桃| 国产在线免费精品| 亚洲免费av在线视频| 2021少妇久久久久久久久久久| 天天影视国产精品| 肉色欧美久久久久久久蜜桃| 欧美另类一区| 国产亚洲一区二区精品| videosex国产| 91麻豆av在线| 亚洲三区欧美一区| 亚洲国产日韩一区二区| 国产97色在线日韩免费| 亚洲欧美日韩另类电影网站| 女人久久www免费人成看片| 女警被强在线播放| 中文乱码字字幕精品一区二区三区| 日韩视频在线欧美| 久久精品亚洲av国产电影网| 亚洲,一卡二卡三卡| 亚洲精品国产av蜜桃| 波野结衣二区三区在线| 高清av免费在线| 一二三四在线观看免费中文在| 午夜福利乱码中文字幕| 亚洲熟女精品中文字幕| 老司机深夜福利视频在线观看 | 高清黄色对白视频在线免费看| 欧美日韩亚洲国产一区二区在线观看 | 国产精品麻豆人妻色哟哟久久| 久久人妻熟女aⅴ| 91麻豆av在线| 国产成人欧美| 国产成人精品在线电影| 亚洲欧洲日产国产| 18禁观看日本| 国产女主播在线喷水免费视频网站| 国产精品免费大片| 欧美成人精品欧美一级黄| 黄色视频不卡| 亚洲综合色网址| 亚洲国产精品成人久久小说| 日韩,欧美,国产一区二区三区| 女警被强在线播放| 水蜜桃什么品种好| 日韩熟女老妇一区二区性免费视频| 极品人妻少妇av视频| 夜夜骑夜夜射夜夜干| 国产xxxxx性猛交| 久久天躁狠狠躁夜夜2o2o | 精品国产一区二区三区久久久樱花| 精品欧美一区二区三区在线| 精品国产超薄肉色丝袜足j| 亚洲国产日韩一区二区| 欧美日韩亚洲国产一区二区在线观看 | 精品久久久久久电影网| 在线观看免费日韩欧美大片| av片东京热男人的天堂| 午夜久久久在线观看| 2018国产大陆天天弄谢| 人人妻人人澡人人看| 男女国产视频网站| 9热在线视频观看99| 极品少妇高潮喷水抽搐| 成人午夜精彩视频在线观看| 日韩熟女老妇一区二区性免费视频| 中文字幕人妻熟女乱码| 50天的宝宝边吃奶边哭怎么回事| 色婷婷久久久亚洲欧美| 蜜桃国产av成人99| a级片在线免费高清观看视频| 99精国产麻豆久久婷婷| 久9热在线精品视频| 色婷婷av一区二区三区视频| 久久国产精品男人的天堂亚洲| 97在线人人人人妻| 精品亚洲乱码少妇综合久久| 久久久国产欧美日韩av| 日韩视频在线欧美| 欧美激情 高清一区二区三区| 日本黄色日本黄色录像| 国产三级黄色录像| 国产成人精品久久二区二区91| 免费高清在线观看视频在线观看| 亚洲欧洲国产日韩| 丰满人妻熟妇乱又伦精品不卡| av网站在线播放免费| 亚洲欧美清纯卡通| 精品亚洲乱码少妇综合久久| 欧美黑人欧美精品刺激| 久久精品亚洲av国产电影网| 一本色道久久久久久精品综合| 黄色毛片三级朝国网站| 欧美精品亚洲一区二区| 一区二区日韩欧美中文字幕| 国产熟女欧美一区二区| 精品久久久久久电影网| 啦啦啦视频在线资源免费观看| 大片免费播放器 马上看| 免费日韩欧美在线观看| 中文乱码字字幕精品一区二区三区| 热re99久久精品国产66热6| 亚洲三区欧美一区| 精品亚洲成a人片在线观看| 国产精品久久久av美女十八| 国产精品久久久人人做人人爽| 国产亚洲精品久久久久5区| 天天操日日干夜夜撸| 中文欧美无线码| 一级毛片 在线播放| 熟女少妇亚洲综合色aaa.| 国产免费视频播放在线视频| 自线自在国产av| 大香蕉久久网| 99香蕉大伊视频| 亚洲色图综合在线观看| 免费在线观看影片大全网站 | 久久久久精品人妻al黑| 中文字幕另类日韩欧美亚洲嫩草| 色网站视频免费| 视频在线观看一区二区三区| 日本色播在线视频| 日韩大片免费观看网站| 九色亚洲精品在线播放| 国产淫语在线视频| 亚洲av电影在线进入| 99精国产麻豆久久婷婷| 亚洲av日韩精品久久久久久密 | 老汉色av国产亚洲站长工具| 久久国产精品男人的天堂亚洲| av在线app专区| 精品福利观看| 老司机在亚洲福利影院| 亚洲午夜精品一区,二区,三区| 精品亚洲成国产av| 宅男免费午夜| 波野结衣二区三区在线| 叶爱在线成人免费视频播放| 精品福利永久在线观看| 欧美在线一区亚洲| 国产高清不卡午夜福利| 国产男女超爽视频在线观看| 99国产综合亚洲精品| 亚洲欧美精品自产自拍| 一区在线观看完整版| 男女国产视频网站| 国产精品一二三区在线看| 日韩视频在线欧美| 久久人妻熟女aⅴ| av福利片在线| 欧美黄色淫秽网站| 啦啦啦中文免费视频观看日本| 又黄又粗又硬又大视频| 国产日韩欧美在线精品| 日本a在线网址| 亚洲精品乱久久久久久| 亚洲欧美一区二区三区久久| 成人国产一区最新在线观看 | 最新的欧美精品一区二区| 亚洲国产精品国产精品| 精品人妻熟女毛片av久久网站| 丝袜人妻中文字幕| 亚洲国产中文字幕在线视频| 国产精品熟女久久久久浪| 欧美黑人精品巨大| 久久免费观看电影| 丰满迷人的少妇在线观看| 丰满饥渴人妻一区二区三| 美国免费a级毛片| 一个人免费看片子| 亚洲av欧美aⅴ国产| 国产精品一区二区免费欧美 | 夜夜骑夜夜射夜夜干| 两个人免费观看高清视频| 又大又爽又粗| 天堂8中文在线网| 热re99久久精品国产66热6| videos熟女内射| 午夜福利,免费看| 永久免费av网站大全| 成人影院久久| 在线观看免费高清a一片| 午夜影院在线不卡| 亚洲美女黄色视频免费看| 少妇猛男粗大的猛烈进出视频| 韩国高清视频一区二区三区| 亚洲成av片中文字幕在线观看| 男女边摸边吃奶| 国产亚洲午夜精品一区二区久久| 亚洲国产精品成人久久小说| 精品一品国产午夜福利视频| 国产成人精品在线电影| 首页视频小说图片口味搜索 | 一区二区日韩欧美中文字幕| 九草在线视频观看| 国产高清不卡午夜福利| 黑丝袜美女国产一区| 欧美日韩一级在线毛片| 国产精品三级大全| 国产高清videossex| 午夜91福利影院| 亚洲精品第二区| 精品亚洲成国产av| 深夜精品福利| 永久免费av网站大全| 只有这里有精品99| 菩萨蛮人人尽说江南好唐韦庄| 欧美日韩亚洲高清精品| 一本大道久久a久久精品| 激情视频va一区二区三区| videos熟女内射| 国产91精品成人一区二区三区 | 大香蕉久久网| 日本av免费视频播放| 大码成人一级视频| 欧美精品亚洲一区二区| 精品亚洲成国产av| 18在线观看网站| 久久青草综合色| 国产亚洲av高清不卡| 五月开心婷婷网| 亚洲国产av新网站| 美女福利国产在线| 精品一区二区三卡| 国产色视频综合| 色播在线永久视频| 无遮挡黄片免费观看| 纵有疾风起免费观看全集完整版| 七月丁香在线播放| 亚洲国产欧美日韩在线播放| 麻豆av在线久日| 在线 av 中文字幕| 国产熟女欧美一区二区| 99国产精品99久久久久| 一二三四社区在线视频社区8| 国产视频首页在线观看| 午夜福利视频在线观看免费| 亚洲欧美成人综合另类久久久| 亚洲成人免费电影在线观看 | 黄色 视频免费看| 亚洲精品日韩在线中文字幕| 一边摸一边抽搐一进一出视频| 香蕉丝袜av| 在线观看免费高清a一片| 精品少妇黑人巨大在线播放| 伦理电影免费视频| 亚洲精品一区蜜桃| 一级,二级,三级黄色视频| 欧美大码av| 悠悠久久av| 人人妻人人爽人人添夜夜欢视频| 国产精品 欧美亚洲| 美女中出高潮动态图| 最新在线观看一区二区三区 | 最新的欧美精品一区二区| 亚洲欧美中文字幕日韩二区| 欧美xxⅹ黑人| 日韩免费高清中文字幕av| 女性生殖器流出的白浆| 亚洲 欧美一区二区三区| 久久av网站| www.熟女人妻精品国产| 久久狼人影院| 亚洲免费av在线视频| av天堂在线播放| 国产成人一区二区三区免费视频网站 | 国产av一区二区精品久久| 欧美日韩亚洲国产一区二区在线观看 | 久久国产精品大桥未久av| 少妇猛男粗大的猛烈进出视频| 国产一区二区三区综合在线观看| 久久久久久久久免费视频了| 亚洲成人国产一区在线观看 | 九草在线视频观看| 亚洲 欧美一区二区三区| 国产真人三级小视频在线观看| 看免费成人av毛片| 精品国产一区二区三区四区第35| 18在线观看网站| 国产高清videossex| 亚洲中文字幕日韩| 黄色视频不卡| 日韩大片免费观看网站| 悠悠久久av| 久久久久久久久久久久大奶| 国产高清videossex| 纯流量卡能插随身wifi吗| 国产在线免费精品| 成人三级做爰电影| 久久国产精品男人的天堂亚洲| 可以免费在线观看a视频的电影网站| 日本91视频免费播放| 性少妇av在线| 女人高潮潮喷娇喘18禁视频| 亚洲伊人久久精品综合| 黑人巨大精品欧美一区二区蜜桃| 久久精品久久久久久噜噜老黄| 久久久久久人人人人人| 麻豆国产av国片精品| 国产一级毛片在线| 亚洲精品一区蜜桃| 又大又黄又爽视频免费| 伊人久久大香线蕉亚洲五| 亚洲欧美成人综合另类久久久| www.自偷自拍.com| 老司机午夜十八禁免费视频| 国产日韩欧美亚洲二区| 成人国语在线视频| 午夜影院在线不卡| 欧美日韩亚洲国产一区二区在线观看 | 国产日韩欧美在线精品| 成人三级做爰电影| 在线观看免费视频网站a站| 男女边吃奶边做爰视频| 久久毛片免费看一区二区三区| 成人影院久久| 亚洲av美国av| 亚洲 欧美一区二区三区| 大香蕉久久网| 久久99热这里只频精品6学生| 看免费成人av毛片| 美女福利国产在线| 欧美日韩综合久久久久久| 一本综合久久免费| 亚洲av日韩在线播放| a级毛片黄视频| 国产成人精品无人区| 777久久人妻少妇嫩草av网站| 中国美女看黄片| 亚洲中文字幕日韩| 免费观看a级毛片全部| 晚上一个人看的免费电影| av福利片在线| 老熟女久久久| av片东京热男人的天堂| 在线观看一区二区三区激情| 亚洲精品乱久久久久久| 老熟女久久久| 美女主播在线视频| 亚洲熟女毛片儿| 大话2 男鬼变身卡| 91精品伊人久久大香线蕉| 老汉色av国产亚洲站长工具| 成人18禁高潮啪啪吃奶动态图| av欧美777| 桃花免费在线播放| 热re99久久精品国产66热6| a级毛片在线看网站| 一边摸一边做爽爽视频免费| 亚洲七黄色美女视频| 一区二区三区乱码不卡18| 亚洲,欧美精品.| 多毛熟女@视频| 免费少妇av软件| 午夜av观看不卡| a级毛片黄视频| 中文字幕另类日韩欧美亚洲嫩草| 啦啦啦视频在线资源免费观看| 一区福利在线观看| 亚洲国产av新网站| 亚洲av综合色区一区| 波多野结衣一区麻豆| www.自偷自拍.com| 狂野欧美激情性xxxx| 久久久精品免费免费高清| 欧美久久黑人一区二区| 国产精品.久久久| 国产在视频线精品| 母亲3免费完整高清在线观看| 麻豆乱淫一区二区| 交换朋友夫妻互换小说| 亚洲欧美清纯卡通| 黑丝袜美女国产一区| 国产av国产精品国产| 赤兔流量卡办理| 一级毛片电影观看| 久久精品国产亚洲av涩爱| 中文字幕av电影在线播放| 岛国毛片在线播放| 国产99久久九九免费精品| 男女床上黄色一级片免费看| 中文精品一卡2卡3卡4更新| 十八禁高潮呻吟视频| 一区在线观看完整版| 亚洲国产成人一精品久久久| 亚洲国产成人一精品久久久| 十八禁网站网址无遮挡| 熟女少妇亚洲综合色aaa.| 久久久国产欧美日韩av| 两人在一起打扑克的视频| 国产又爽黄色视频| 久久久久国产一级毛片高清牌| 男女午夜视频在线观看| 亚洲国产av新网站| 中文字幕亚洲精品专区| 免费高清在线观看视频在线观看| 一级a爱视频在线免费观看| 日韩av在线免费看完整版不卡| 国产亚洲一区二区精品| 日本vs欧美在线观看视频| 日韩视频在线欧美| 国产成人av激情在线播放| 视频区欧美日本亚洲| 欧美精品啪啪一区二区三区 | 女人久久www免费人成看片| 50天的宝宝边吃奶边哭怎么回事| 久久毛片免费看一区二区三区| 成人三级做爰电影| 免费在线观看完整版高清| 久久国产亚洲av麻豆专区| 亚洲欧美一区二区三区黑人| 韩国高清视频一区二区三区| 国产成人91sexporn| 精品一品国产午夜福利视频| av在线播放精品| 免费观看av网站的网址| 亚洲欧美精品综合一区二区三区| 亚洲人成电影免费在线| 日韩 欧美 亚洲 中文字幕| 国产人伦9x9x在线观看| 午夜福利免费观看在线| 高清欧美精品videossex| 亚洲激情五月婷婷啪啪| 嫩草影视91久久| 美女扒开内裤让男人捅视频| 啦啦啦 在线观看视频| 夜夜骑夜夜射夜夜干| 这个男人来自地球电影免费观看| 嫩草影视91久久| 亚洲国产av影院在线观看| 十八禁高潮呻吟视频| 婷婷丁香在线五月| 欧美人与善性xxx| 亚洲欧美清纯卡通| 中文欧美无线码| 波多野结衣av一区二区av| 成人18禁高潮啪啪吃奶动态图| 久久女婷五月综合色啪小说| 777久久人妻少妇嫩草av网站| 在线观看免费午夜福利视频| 精品一品国产午夜福利视频| 黄色毛片三级朝国网站| 亚洲av欧美aⅴ国产| 在线亚洲精品国产二区图片欧美| 免费久久久久久久精品成人欧美视频| 国产在线免费精品| 精品福利观看| 精品国产一区二区久久| 亚洲精品久久成人aⅴ小说| 黄色怎么调成土黄色| 你懂的网址亚洲精品在线观看| 亚洲中文字幕日韩| 亚洲成国产人片在线观看| 人妻一区二区av| 夜夜骑夜夜射夜夜干| 两性夫妻黄色片| 国产日韩欧美在线精品| 一区二区三区四区激情视频| 香蕉丝袜av| 老司机在亚洲福利影院| 汤姆久久久久久久影院中文字幕| 麻豆乱淫一区二区| 青草久久国产| 韩国高清视频一区二区三区| 婷婷色综合大香蕉| 乱人伦中国视频| 精品熟女少妇八av免费久了| 亚洲av综合色区一区| 精品熟女少妇八av免费久了| 极品少妇高潮喷水抽搐| 国产精品三级大全| 99re6热这里在线精品视频| 中文字幕高清在线视频| 亚洲一码二码三码区别大吗| 国产亚洲一区二区精品| 男人操女人黄网站| 国产亚洲精品第一综合不卡| 国产精品亚洲av一区麻豆| 久久亚洲精品不卡| 脱女人内裤的视频| 美女高潮到喷水免费观看| 国产1区2区3区精品| 国产亚洲欧美在线一区二区| 日本欧美国产在线视频| 这个男人来自地球电影免费观看| 免费在线观看影片大全网站 | 日韩免费高清中文字幕av| 国产成人精品无人区| 少妇 在线观看| 久久九九热精品免费| 久久人人97超碰香蕉20202| 亚洲色图综合在线观看| 国产在线免费精品| 精品国产国语对白av| 纵有疾风起免费观看全集完整版| 成年av动漫网址| www.熟女人妻精品国产| 香蕉丝袜av| 老司机靠b影院| 黑人猛操日本美女一级片| 亚洲人成77777在线视频| 久久久精品区二区三区| 久久久国产一区二区| 视频区欧美日本亚洲| a 毛片基地| 麻豆国产av国片精品| 国产野战对白在线观看| 啦啦啦在线观看免费高清www| 蜜桃在线观看..| 91精品国产国语对白视频| 亚洲欧美一区二区三区黑人| 亚洲欧美日韩另类电影网站| 波多野结衣一区麻豆| 国产在线一区二区三区精| 久久久久国产一级毛片高清牌| 日韩,欧美,国产一区二区三区| 亚洲精品一卡2卡三卡4卡5卡 | www.精华液| 99精国产麻豆久久婷婷| 丝袜在线中文字幕| av一本久久久久| 天堂中文最新版在线下载| 成人国产一区最新在线观看 | 国产熟女欧美一区二区| 少妇粗大呻吟视频| 欧美少妇被猛烈插入视频| 少妇裸体淫交视频免费看高清 | 亚洲精品久久成人aⅴ小说| 丰满少妇做爰视频| 亚洲色图 男人天堂 中文字幕| 成人午夜精彩视频在线观看| 国产精品久久久久成人av| 精品熟女少妇八av免费久了|