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

    低階煤孔隙結(jié)構(gòu)定量表征及瓦斯吸附放散特性

    2021-11-10 03:22:14蔣靜宇程遠(yuǎn)平
    煤炭學(xué)報 2021年10期
    關(guān)鍵詞:孔容低階煤樣

    蔣靜宇,程遠(yuǎn)平,張 碩

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

    在已探明煤炭儲量中,低階煤廣泛分布在我國西北、華北和東北地區(qū)[1],其種類主要為褐煤和次煙煤(包括長焰煤、不黏煤和弱黏煤)[2],富含煤層氣資源[3]。在我國低階煤礦區(qū)如鐵法礦區(qū)[4]、沈北礦區(qū)[5]等發(fā)生多起煤與瓦斯突出事故。由于煤是內(nèi)部結(jié)構(gòu)復(fù)雜的多孔介質(zhì)[6],氣體在煤儲層中賦存和運(yùn)移受吸附和擴(kuò)散過程的控制。而煤巖孔隙結(jié)構(gòu)對上述過程有關(guān)鍵控制作用[7-8],其孔隙數(shù)量及空間分布的非均質(zhì)性影響礦井瓦斯異常富集和瓦斯災(zāi)害[9]。開展孔隙結(jié)構(gòu)表征及瓦斯放散特性研究對低階煤瓦斯抽采和災(zāi)害防治意義重大。

    國際純粹與應(yīng)用化學(xué)聯(lián)合會(IUPAC)孔隙分類標(biāo)準(zhǔn)為:大孔 (>50 nm)、介孔 (2~50 nm)和微孔 (≤2 nm)[10]。ZHAO等[11]研究表明隨著煤階增加,比表面積呈現(xiàn)先減小后增大的趨勢。JIANG等[12]采用壓汞(MIP)和N2/CO2吸附法研究了煤化作用對納米孔隙演化的影響,表明第3次煤化躍遷對微孔發(fā)育具有關(guān)鍵扭轉(zhuǎn)作用,而第5次煤化躍遷顯著提高介孔孔容。LI等[13]采用N2/CO2吸附法和小角度X射線散射(SAXS)法表征煤中納米孔隙,表明閉孔數(shù)量占比超過94%。CO2吸附法可測定小于1.5 nm 的微孔,N2吸附法可測定小于300 nm的介孔,SAXS可無損表征煤中納米孔隙[14-15]。尚缺少將MIP,N2/CO2吸附法和SAXS表征手段結(jié)合綜合表征低階煤的孔隙結(jié)構(gòu)的相關(guān)研究。

    YAO等[7]基于N2吸附試驗結(jié)果,應(yīng)用Frenkel-Halsey-Hill (FHH)模型計算得到分形維數(shù)D1和D2,D1值越大表明孔隙表面越不規(guī)則,瓦斯吸附量越高,而D2值越高,瓦斯吸附能力越差。FU等[3]將N2吸附法和MIP測得的煤樣孔隙參數(shù)與分形理論結(jié)合,指出低階煤孔隙結(jié)構(gòu)有利于煤層氣富集但不利于煤層氣滲流。YI等[16]研究表明酸處理后煤樣微孔更加發(fā)育,且促進(jìn)許多封閉介孔打開,導(dǎo)致D2增加,比表面積亦增加。宋曉夏等[17]基于SAXS法研究構(gòu)造煤的孔隙結(jié)構(gòu),表明煤中微孔比例增加時,最可幾孔徑減小,孔隙表面分形維數(shù)增大。

    平均隨機(jī)鏡質(zhì)體反射率(Ro)在0.5%~0.6%時,發(fā)生第1次煤化躍變(FCJ),煤的分子結(jié)構(gòu)和孔隙會發(fā)生顯著改變[18],繼而影響吸附瓦斯能力[19]。煤化作用和顯微組分對煤的瓦斯吸附能力有重要影響,朗格繆爾壓力(PL)主要受顯微組分的控制,而朗格繆爾體積(VL)則受煤化程度的影響較大[20]。張慧杰等[21]研究表明,煤中介孔和大孔分布特征是影響煤體瓦斯放散特性的主要因素。有關(guān)低階煤的孔隙結(jié)構(gòu)及瓦斯放散特性,有關(guān)學(xué)者做了較好的前期研究,而尚缺少孔隙結(jié)構(gòu)變化與瓦斯吸附/放散特性相關(guān)性的深入研究。系統(tǒng)研究低階煤孔隙結(jié)構(gòu)和瓦斯吸附/放散特性對低階煤突出防治有重要意義。

    筆者以4組低階煤和2組中階煤煤樣為研究對象,采用N2/CO2吸附法、SAXS法和甲烷吸附/解吸試驗,引入分形理論,研究孔隙和瓦斯放散特性,為煤礦低階煤瓦斯災(zāi)害防治提供借鑒。

    1 取樣和試驗方法

    1.1 取樣與參數(shù)測定

    試驗煤樣分別取自遼寧省沈北煤田蒲河煤礦、鐵法礦區(qū)小青煤礦和大興煤礦、沈北紅陽二礦及河南省義馬礦區(qū)新義煤礦,分別標(biāo)記為PH01,XQ02,DX03,DX04,HY05和XY06。根據(jù)國際化學(xué)聯(lián)合會標(biāo)準(zhǔn)(ASTM,2007),采用5E-MAG6600全自動工業(yè)分析儀完成工業(yè)分析。根據(jù)中國標(biāo)準(zhǔn)GB/T 476—2008《煤中碳和氫的測定方法》,采用德國Elementar公司的元素分析儀進(jìn)行元素分析。參照國際標(biāo)準(zhǔn)ISO 7404-5,1984《煤的巖相分析法》,使用顯微鏡光度計(Zeiss,Germany)測定煤樣有機(jī)顯微組分和平均隨機(jī)鏡質(zhì)體反射率(Ro),煤的高壓等溫吸附試驗參考GB/T 19560—2008《煤的高壓等溫吸附試驗方法》。試驗主要結(jié)果見表1。

    通常認(rèn)為低階煤是Ro,max<0.65%的煤[2]。由表1可知,6組煤樣Ro測定結(jié)果為0.43%~1.09%,結(jié)合文獻(xiàn)[2]的分類方法,煤樣PH01,XQ02,DX03,DX04屬于低階煤,HY05和XY06屬于中階煤,作為對比煤樣。低階煤含水率為5.45%~15.12%,尤其是褐煤PH01含水率高達(dá)15.12%。中階煤含水率相對較低,為1.30%~3.03%。低階煤的揮發(fā)分V為28.04%~42.50%,均高于中階煤揮發(fā)分(20.67%~21.78%)。低階煤的含硫量小于中階煤,低階煤S元素質(zhì)量分?jǐn)?shù)為0.47%~0.83%,中階煤S元素質(zhì)量分?jǐn)?shù)為1.67%~2.88%。低階煤樣中鏡質(zhì)組(V)質(zhì)量分?jǐn)?shù)為69.3%~96.08%,總體上高于中階煤樣中鏡質(zhì)組質(zhì)量分?jǐn)?shù)(66.42%~75.93%)。

    1.2 試驗方法及原理

    采用N2/CO2吸附法和小角度X射線散射(SAXS)法綜合表征低階煤微觀孔隙結(jié)構(gòu);開展甲烷吸附/解吸試驗,結(jié)合孔隙試驗結(jié)果,分析低階煤瓦斯賦存和放散特性。

    (1)物理吸附法:采用N2/CO2吸附對煤樣的介孔和微孔進(jìn)行測定。試驗設(shè)備是全自動氣體吸附分析儀(AUTOSORB-IQ2)。N2吸附試驗是在77 K下進(jìn)行,測定孔徑范圍是<300 nm。樣品粒徑為0.20~0.25 mm,試驗前干燥處理。CO2吸附試驗是在溫度273 K、壓力低于3 039.75 Pa進(jìn)行,測定孔徑為0.35~1.50 nm。

    氣體吸附法是根據(jù)物理吸附等溫線的形態(tài)變化來獲得吸附劑的表面性質(zhì)與孔徑分布等參數(shù)[22]。溫度不變,相對壓力P/P0改變時,煤中氣體吸附量會隨之發(fā)生改變[23-24]。N2/CO2吸附試驗通過記錄不同的壓力下的氣體吸附量得到吸附等溫線,采用BJH模型和DFT模型可計算得到煤的孔徑分布、孔容和比表面積等參數(shù)[25]。

    (2)小角度X射線散射(SAXS)試驗:試驗在中國科學(xué)院高能物理研究所北京同步輻射中心小角散射試驗站完成。煤樣粒徑<0.074 mm,測試孔隙半徑為1~30 nm(直徑為2~60 nm)。首先記錄X射線波長、相機(jī)長度以及背底。完成背底測量后進(jìn)行樣品測量,用于對散射強(qiáng)度I(q)扣除背底處理[26]。接著將二維圖像轉(zhuǎn)換為散射強(qiáng)度I與距離r數(shù)據(jù)。散射花樣圖可表達(dá)為散射矢量q的變量[27]。通過分析q-I(q)曲線,可得到煤樣孔徑分布和分形維數(shù)等參數(shù)[28]。由于散射體的復(fù)雜性,實際散射結(jié)果會產(chǎn)生一定偏離[28]。針對偏離的校正方法參考李志宏等[29]提出的方法。對多分散體系來說,當(dāng)分散體之間的距離大于其本身的尺寸,X射線散射強(qiáng)度可近似為單分散體散射強(qiáng)度的累加[17,28]。使用ASTS軟件可計算得到煤樣的孔徑分布,其基于球形多分散體的分布函數(shù)[30]為

    (1)

    式中,P(R)為散射體半徑的體積分布;R為球半徑nm;N(R)為半徑為R的粒子的相對數(shù)量系統(tǒng)(π/qmax

    (3)甲烷吸附/解吸試驗:實驗室篩分粒徑為0.20~0.25 mm、質(zhì)量為50 g的煤樣,干燥脫氣處理后進(jìn)行甲烷等溫(30 ℃)吸附試驗,可得到VL和PL。甲烷解吸試驗步驟為:篩分出粒徑為1~3 mm、質(zhì)量為50 g的煤樣后將樣品放到密封罐,抽真空處理。接著讓樣品吸附甲烷氣體,直至達(dá)到吸附平衡后瞬間釋放壓力,統(tǒng)計120 min內(nèi)不同時間點的甲烷解吸量。將解吸量換算成標(biāo)準(zhǔn)狀態(tài)下解吸量后可獲得煤樣解吸瓦斯量與時間的關(guān)系曲線。

    基于煤樣瓦斯解吸曲線圖,采用單孔擴(kuò)散模型[31]可計算得到煤體瓦斯擴(kuò)散系數(shù)(D),該模型表達(dá)式為

    (2)

    式中,Vt為t時間內(nèi)單位質(zhì)量顆粒煤的累計擴(kuò)散瓦斯總量,cm3/g;V∞為無窮大時刻的單位質(zhì)量顆粒煤的累計擴(kuò)散瓦斯總量,cm3/g;D為擴(kuò)散系數(shù),cm2/s;rp為瓦斯擴(kuò)散粒徑,cm。

    對于初始階段(t<600 s),Vt/V∞<0.5的單孔擴(kuò)散模型[31],式(2)可以簡化為

    (3)

    2 結(jié)果分析與討論

    2.1 N2/CO2吸附試驗結(jié)果

    2.1.1N2吸附試驗結(jié)果和孔隙特征

    依據(jù)77 K溫度下測定的相對壓力及相應(yīng)吸附量繪制得到的N2等溫吸附線分別選用BJH模型和DFT模型可得到煤樣孔徑分布曲線,如圖1所示。

    圖1 N2吸附試驗孔徑分布曲線Fig.1 Pore size distribution curves of N2 adsorption

    BJH模型分析結(jié)果顯示煤樣孔徑呈多峰分布,低階煤孔徑主要分布在2~8 nm,煤樣PH01除外,其最發(fā)育孔徑多為1~5 nm;而2個中階煤樣最發(fā)育孔徑為2~4 nm。DFT模型分析結(jié)果表明試驗煤樣孔徑分布規(guī)律基本上和BJH分析結(jié)果基本保持一致,但低階煤DX03和中階煤XY06最發(fā)育孔徑主要為4~36 nm的介孔。表2列出了基于N2吸附試驗結(jié)果得到的孔隙結(jié)構(gòu)參數(shù)。

    表2 N2吸附孔隙結(jié)構(gòu)參數(shù)(孔徑0.7~300 nm)

    低階煤與中階煤的孔容及比表面積存在差異,低階煤BJH總孔容為(4.46~40.84)×10-3cm3/g大于中階煤總孔容(5.76~9.13)×10-3cm3/g;低階煤的BET比表面積為1.16~38.14 m2/g,大于中階煤BET比表面積(0.91~2.35 m2/g)。低階煤樣PH01和XQ02中介孔(2~50 nm)孔容較大,遠(yuǎn)大于2個中階煤的介孔孔容。低階煤樣的大孔(>50 nm)孔容和中階煤大孔孔容大小規(guī)律不明顯。單位質(zhì)量煤體的孔容在一定程度上能反映煤體的瓦斯吸附能力,為研究低變質(zhì)程度煤的微孔、介孔孔容大小與變質(zhì)程度的關(guān)系,基于本研究N2吸附試驗及相關(guān)文獻(xiàn)數(shù)據(jù)[11-12,32]繪制了低中階煤樣孔容與Ro的關(guān)系曲線,如圖2所示。

    圖2 低中階煤基于N2吸附試驗的孔容與Ro關(guān)系曲線Fig.2 Relationship between pore volume and Ro based on N2 adsorption experiment of low-medium rank coals

    隨著煤化程度的提高,孔容總體呈現(xiàn)先減小后增大的趨勢。且相關(guān)研究顯示,低階煤的介孔孔容和比表面積普遍高于煤化程度較低的中階煤,且隨著煤階的升高,總孔容和比表面積呈先減小后增大的趨勢[11, 33],而介孔是影響煤體瓦斯放散特性的主要因素之一[21],因此低階煤中瓦斯?jié)B流過程較部分中階煤更為通暢,其瓦斯放散能力較高于部分中階煤。

    2.1.2CO2吸附試驗結(jié)果和孔隙特征

    采用CO2吸附法測得煤樣的CO2等溫吸附曲線如圖3所示。

    圖3 煤樣CO2等溫吸附曲線Fig.3 CO2 adsorption isotherm of coal samples

    除HY05煤樣外,煤樣等溫吸附曲線形狀較為近似于朗格繆爾曲線,在低相對壓力區(qū)吸附量快速增加、吸附線明顯向上突起,而在高壓區(qū)近似趨于直線(圖3)。HY05煤樣的單位質(zhì)量CO2吸附體積明顯低于其他煤樣,表明其微孔發(fā)育最差。利用ASiQwin軟件,基于DFT模型可計算煤體微孔孔容、最可幾孔徑及比表面積。CO2吸附法測得的煤樣基于孔容的孔徑分布曲線如圖4所示。

    圖4 CO2吸附法測得煤樣微孔孔徑分布Fig.4 Pore size distribution of measured by CO2 adsorption

    6組煤樣孔徑分布均呈多峰分布,孔徑多集中分布在0.40~0.65 nm和0.85 nm,這與李陽等[34]試驗結(jié)果一致。煤樣在0.90~1.50 nm孔徑分布基本相似,而在0.35~0.70 nm的孔徑分布存在波動性,且此孔徑范圍內(nèi)處于FCJ附近的煤樣其孔容波動性較為明顯 (圖4中虛線圈)。有研究表明,Ro分別在0.5%~0.6%,1.1%~1.2%,1.4%~1.5%,2.2%~2.3%,2.8%~2.9%,3.7%~3.9%[35]及4.2%~4.8%[36]時會發(fā)生煤化躍變。這可以解釋圖4圓圈內(nèi)的數(shù)據(jù)點均屬于FCJ附近的煤樣。當(dāng)Ro<1.4%時,煤樣微孔受芳香族結(jié)構(gòu)控制,官能團(tuán)的解離和聚合作用也會影響微孔演化[37],F(xiàn)CJ對絲炭的孔隙結(jié)構(gòu)影響不大,但對原煤的孔隙結(jié)構(gòu)影響明顯[38]。表3列出了CO2吸附試驗測得的孔隙結(jié)構(gòu)參數(shù)。

    表3 CO2吸附試驗孔隙結(jié)構(gòu)參數(shù)(孔徑0.35~1.50 nm)

    煤樣的最可幾(最發(fā)育)孔徑在0.47~0.57 nm,數(shù)據(jù)較離散,但中階煤XY06最發(fā)育孔徑最小,為0.47 nm。低階煤的微孔孔容和比表面積相對較大,低階煤的微孔孔容為0.055~0.064 cm3/g,總體上高于中階煤的微孔孔容(0.026~0.060 cm3/g)。基于CO2吸附試驗結(jié)果以及國內(nèi)外文獻(xiàn)相關(guān)數(shù)據(jù)[11-12,20, 39],繪制低中階煤樣微孔孔容與Ro的關(guān)系曲線,如圖5所示。

    圖5 低中階煤基于CO2吸附試驗的孔容與Ro關(guān)系曲線Fig.5 Relationship between pore volume and Ro based on CO2 adsorption of low-medium rank coals

    低階煤的微孔孔容和比表面積高于部分煤化程度較低的中階煤,且總體上隨煤化程度的提高呈現(xiàn)先減小后增大的趨勢。微孔孔容可反映煤體的氣體吸附能力,因而低階煤的瓦斯吸附能力較高于部分煤化程度較低的中階煤。

    2.2 小角X射線散射測試結(jié)果

    由于N2/CO2吸附法需要利用氣體分子探針進(jìn)入開放型孔隙,導(dǎo)致煤中超微孔和閉孔難以測定,SAXS可彌補(bǔ)上述不足[40]。SAXS可進(jìn)行特大晶胞物質(zhì)的結(jié)構(gòu)分析以及測定粒度在幾十納米以下超細(xì)粉末粒子大小、形狀及分布[41]。選取3個低階煤PH01,XQ02和DX03及中階煤XY06煤樣的SAXS二維干涉花樣成像和扣除背底后的散射強(qiáng)度曲線如圖6所示。

    圖6 煤樣SAXS散射花樣成像結(jié)果和扣除背底后的散射強(qiáng)度曲線Fig.6 Imaging results of SAXS scattering pattern of coal samples and scattering intensity curves after subtracting background

    煤樣的散射強(qiáng)度隨散射矢量的增加而降低,依據(jù)Porod和Guinier定理,可得到煤樣散射曲線[28]。根據(jù)Porod和Guinier散射曲線可分析煤樣的孔徑分布規(guī)律。使用ASTS軟件可計算得到煤樣2~60 nm的孔徑分布曲線,如圖7所示。

    圖7 SAXS分析的煤樣孔徑分布Fig.7 Pore size distribution of coals analyzed by SAXS

    4個煤樣孔徑分布基本呈現(xiàn)雙峰分布,低階煤PH01,XQ02和 DX03的介孔較中階煤XY06介孔更發(fā)育。從分布峰的位置可以看出,3組低階煤樣的最發(fā)育孔徑在10 nm左右,而中階煤XY06最發(fā)育孔徑在35 nm左右,遠(yuǎn)大于N2吸附法測得低階煤和中階煤的最發(fā)育孔徑(分別為2.92~3.26,1.63~2.99 nm)。表明煤中可能含有一些封閉或氣體分子探針無法測到的孔[37]。氣體吸附法有可能測不到煤中閉孔或部分口小腔大的半開孔,造成煤體孔隙率計算可能偏小,影響煤層氣儲量或煤體瓦斯含量的計算。

    2.3 低階煤孔隙結(jié)構(gòu)的分形特征

    2.3.1基于N2吸附試驗的孔隙分形特征

    基于N2吸附試驗數(shù)據(jù)應(yīng)用FHH模型,可計算孔隙的分形維數(shù)D[42],計算式為

    (4)

    式中,V為平衡壓力P下吸附的氣體分子體積,mL/g;V0為使用BET方程計算的單層體積,mL/g。

    通過線性擬合N2吸附等溫線的ln(ln(P0/P))和ln(V/V0)可得到斜率S,根據(jù)式D=3+S計算得到分形維數(shù)D。煤樣分形維數(shù)計算結(jié)果如圖8所示。圖8中曲線均在ln(ln(P0/P))=-0.5處表現(xiàn)出突變,該處對應(yīng)孔徑約為5 nm,前后2處分為2部分,分別得到分形維數(shù)D1(>5 nm)和D2(<5 nm)。

    圖8 基于N2吸附試驗結(jié)果的煤樣分形維數(shù)Fig.8 Fractal dimension of coal samples based on results of N2adsorption experiment

    表4列出了煤樣分形維數(shù)D1和D2計算結(jié)果。

    表4 基于N2吸附試驗結(jié)果計算的煤樣分形維數(shù)

    試驗煤樣D1和D2分別為2.44~2.92和2.03~2.61。表面越光滑,分形維數(shù)越接近于2,表面越粗糙,分形維數(shù)則趨近于3[43]。分形維數(shù)D1和D2分別代表不同孔徑段的孔隙復(fù)雜程度(D1>5 nm,D2<5 nm)。研究表明D1反映煤基質(zhì)中毛細(xì)凝聚及滲流情況,D1越大,孔隙結(jié)構(gòu)的非均質(zhì)性越強(qiáng),氣/液表面張力越大,甲烷吸附量越低;而D2越大,則表面越不規(guī)則,為瓦斯吸附提供更多空間[7, 44]。除了DX04煤樣,低階煤的分形維數(shù)D2值(2.25~2.61)均大于中階煤的D2值(2.04和2.10),表明低階煤<5 nm的吸附孔隙發(fā)育、表面較為復(fù)雜。而分形維數(shù)D1和變質(zhì)程度沒有明顯規(guī)律,可能煤階不是影響滲流孔隙的主要因素。

    2.3.2基于SAXS試驗的孔隙分形特征

    分形物體在分形區(qū)的小角X射線散射強(qiáng)度符合[28]:

    I(q)=I0q-α

    (5)

    其中,α為衰減指數(shù),在lnI(q)-lnq曲線的線性區(qū)域斜率為d,則α= -d,α值為0~4。當(dāng)3<α<4時,物體表現(xiàn)為表面分形,分形維數(shù)DS=6-α;當(dāng)0<α<3時,物體表現(xiàn)為孔/質(zhì)量分形,分形維數(shù)Dp/m=α。如圖9所示,分形維數(shù)可利用SAXS試驗數(shù)據(jù)作lnI(q)-lnq圖獲得,通過線性擬合得到斜率d,所得α值均處于3~4,然后可計算得到分形維數(shù)DS。

    圖9 基于小角X射線散射試驗結(jié)果的煤樣分形維數(shù)Fig.9 Fractal dimension of coal sample based on result of small-angle X-ray scattering experiment

    表5列出了4組煤樣的衰減指數(shù)α和分形維數(shù)DS。α值為3.29~3.59,DS為2.41~2.72。

    分形維數(shù)描述孔隙表面光滑度和平整度,體現(xiàn)著煤樣孔隙的發(fā)育復(fù)雜程度[28]。DS值越接近3表明其表面越粗糙,孔隙發(fā)育更加復(fù)雜。低階煤PH01,XQ02和DX03的DS值分別為2.44,2.72和2.71,均大于XY06的DS值(2.41) (表5),表明低階煤樣的孔隙結(jié)構(gòu)更為復(fù)雜。這主要是受FCJ的影響,煤化躍變導(dǎo)致煤樣微觀分子結(jié)構(gòu)中脂肪族和羥基的伸縮性降低[35],因而XY06煤樣的DS值最小。

    表5 煤樣基于小角X射線散射試驗孔隙分形維數(shù)

    2.4 低階煤孔隙結(jié)構(gòu)對瓦斯吸附及放散特性的影響

    6個煤樣的瓦斯吸附等溫曲線如圖10所示。在恒定溫度和不同壓力條件下的瓦斯等溫吸附曲線均符合朗格繆爾曲線特征。

    圖10 低階煤樣的瓦斯等溫吸附曲線Fig.10 Gas adsorption isotherms of low-rank coal samples

    6個煤樣瓦斯吸附參數(shù)測定結(jié)果見表6。

    表6 煤樣瓦斯吸附特性和基于單孔擴(kuò)散模型的擴(kuò)散系數(shù)

    表6中,VL值越大表明煤體吸附瓦斯能力越強(qiáng)。除DX04煤樣(VL值為21.04 m3/t)外,低階煤的VL值為31.14~35.14 m3/t,總體上大于中階煤的VL值(15.80~28.21 m3/t)?;谠囼灉y得的VL、微孔孔容和比表面積,繪制了VL分別與微孔孔容和比表面積的關(guān)系曲線,如圖11所示。

    圖11 低中階煤VL分別與微孔孔容和比表面積關(guān)系曲線Fig.11 Relationships between VL and micropore volume and specific surface area of low-medium rank coals respectively

    由圖11可知,試驗煤樣微孔孔容和比表面積均與朗格繆爾體積VL存在正相關(guān)關(guān)系,即微孔孔容和比表面積越大,煤體吸附瓦斯能力越強(qiáng)。微孔在很大程度上決定煤體的瓦斯吸附能力[2]。試驗煤樣低階煤的微孔較中階煤發(fā)育,瓦斯吸附能力總體上高于試驗中階煤煤樣,表明低階煤易于瓦斯的儲存和富集。

    瓦斯解吸試驗分別設(shè)定了初始吸附平衡壓力為1,2,3,4,5 MPa。6個煤樣(干燥基)不同瓦斯壓力條件下120 min內(nèi)的累計瓦斯解吸曲線如圖12所示。

    不同吸附平衡瓦斯壓力條件下,120 min內(nèi)解吸瓦斯量由大到小的順序為XQ02>PH01> DX03>XY06>DX04>HY05。前3組低階煤的累積瓦斯解吸量較大,且初始瓦斯解吸速度快(圖12(a)~(c)),明顯高于中階煤的初始瓦斯解吸速度及累積瓦斯解吸量(圖12(e),(f))。試驗低階煤樣瓦斯放散能力總體上強(qiáng)于中階煤樣,表明低階煤體孔隙連通性較好,能夠在較短的時間內(nèi)解吸大量瓦斯。

    將煤體假設(shè)為半徑相等且表面光滑的均質(zhì)球狀顆粒或孔隙聚合而成的多孔介質(zhì)?;谖狡胶鈮毫? MPa時煤粒累計解吸瓦斯量與時間的關(guān)系曲線(圖12),引入單孔擴(kuò)散模型[31],可計算6個煤樣的常擴(kuò)散系數(shù)(D),計算結(jié)果見表6。除了DX04煤樣,低階煤樣的D為(13.43~32.05)×10-7cm2/s,高于試驗中階煤樣的常擴(kuò)散系數(shù)(12.77×10-7cm2/s和13.07×10-7cm2/s),低階煤的介孔和大孔孔容整體上高于試驗中階煤的介孔和大孔孔容。為研究介孔、大孔孔容與瓦斯放散特性的關(guān)系,繪制試驗煤樣介孔、大孔孔容與煤樣擴(kuò)散系數(shù)D的關(guān)系曲線,如圖13所示。

    圖12 煤樣不同吸附平衡壓力下瓦斯解吸曲線Fig.12 Gas desorption curves of different coal samples under different equilibrium pressures

    試驗煤樣的擴(kuò)散系數(shù)隨著總孔容減小具有波動性減小的趨勢,擴(kuò)散系數(shù)與介孔孔容和大孔孔容存在正相關(guān)性(圖13)。DX04煤樣的擴(kuò)散系數(shù)較小則是由于其大孔孔容較小。試驗結(jié)果表明,低階煤體相對于低階煤體具有較高的瓦斯放散能力,易于低階煤瓦斯突出過程的突出連續(xù)發(fā)展。

    2.5 低階煤分子結(jié)構(gòu)演化對瓦斯吸附放散特性的影響

    隨著煤化程度的提高,煤體微觀分子結(jié)構(gòu)會發(fā)生改變,可導(dǎo)致煤基質(zhì)孔隙結(jié)構(gòu)發(fā)生改變。采用傅里葉變換紅外光譜(FTIR)技術(shù)可表征煤樣微觀分子結(jié)構(gòu),4組代表性煤樣的FTIR譜圖如圖14[45]所示。

    根據(jù)吸收波段的不同可劃分出煤樣分子結(jié)構(gòu)中的羥基(3 600~3 000 cm-1)、脂肪族(3 000~2 700 cm-1)、含氧官能團(tuán)(1 800~1 000 cm-1)和芳香族(900~700 cm-1),如圖14所示。

    圖14 煤樣傅里葉變換紅外光譜譜圖[45]Fig.14 Fourier transform infrared spectrum of coal[45]

    對各波段譜峰進(jìn)行Gaussian分峰擬合后可計算出各官能團(tuán)的譜峰面積。在獲得峰面積基礎(chǔ)上,可計算出芳香度fa、脂氫比Hal/H、單碳環(huán)數(shù)(R/C)u、脂肪族結(jié)構(gòu)參數(shù)(CH2/CH3,CH2與CH3結(jié)構(gòu)在紅外光譜中的面積比)和含氧官能團(tuán)結(jié)構(gòu)參數(shù)′C′,計算結(jié)果見表7。

    其中fa和(R/C)u能在一定程度上反映煤的大分子結(jié)構(gòu)中芳香族結(jié)構(gòu)的信息,CH2/CH3可以反映煤中脂肪族側(cè)鏈的長度。從表7可知,YMXY煤樣的fa最大,表明中階煤樣擁有較高的芳構(gòu)化程度。盡管低階煤TFXQ的fa相對較大,但其單碳環(huán)數(shù)(R/C)u偏小。而Hal/H值可反映出低階煤TFDX和SBPH的脂肪族含量比2個中階煤樣高。低階煤SBPH的CH2/CH3和′C′值均最大,反映其擁有較長的脂肪族側(cè)鏈及含氧官能團(tuán)側(cè)鏈。

    表7 基于FTIR結(jié)果計算得到的煤樣分子結(jié)構(gòu)參數(shù)[45]

    研究表明當(dāng)Ro<1.4%時,煤體微孔主要受到脂肪族官能團(tuán)的控制[36]。在第1次煤化躍變(FCJ)之前,煤中芳香族含量較低,含氧官能團(tuán)和脂肪族官能團(tuán)占比較高,微觀分子結(jié)構(gòu)中存在較多的長側(cè)鏈,分子結(jié)構(gòu)較為松散[2,20],因而其微孔孔容和比表面積較大;而在FCJ之后,煤樣微孔孔容不斷減小(圖5),是由于煤分子中的長側(cè)鏈逐漸分解(CH2/CH3及′C′逐漸減小),芳香族占比逐漸增大,煤體分子結(jié)構(gòu)變得緊實,導(dǎo)致微孔孔容和比表面積均減小,進(jìn)而降低煤體瓦斯吸附能力。隨著煤化程度的增加,煤體介孔孔容總體上呈減小的趨勢,且在FCJ前的孔容減小趨勢較為明顯,擬合直線斜率較大(圖2)。表明低階煤的介孔較部分中階煤更發(fā)育。

    試驗低階煤樣的水分為5.45%~15.12%,明顯高于中階煤的水分(1.03%~3.03%,表1)。這驗證了在FCJ過程中煤分子發(fā)生縮聚,煤基質(zhì)不斷壓實,介孔孔容、煤中水分含量和孔隙率隨之降低,孔隙連通性變差[19],導(dǎo)致FCJ后中階煤體瓦斯放散能力的降低。采用SAXS無損測出的低階煤樣品2~28 nm的介孔孔容均明顯高于中階煤XY06的介孔孔容(圖7)。對比N2吸附試驗結(jié)果,表明煤中可能存在氣體吸附法無法測得的封閉或半封閉孔隙。煤儲層在采動或卸壓等條件下,這些孔隙很可能打開,導(dǎo)致瓦斯流動性增強(qiáng)。

    試驗低階煤樣的微孔、介孔孔容和比表面積總體上大于中階煤,宏觀表現(xiàn)為低階煤較中階煤具有更高的吸附能力和瓦斯放散特性,易發(fā)生瓦斯富集,在采動等影響下,易導(dǎo)致瓦斯的快速解吸,從而有利于突出的連續(xù)性發(fā)展。綜上所述,低階煤具有初始瓦斯解吸速度快的特點,雖有利于抽采瓦斯,但受采動影響條件下,井下低階煤層易于發(fā)生瓦斯異常涌出或突出災(zāi)害事故。

    3 結(jié) 論

    (1)N2吸附試驗結(jié)果顯示,低階煤煤樣介孔孔徑主要分布在2~8 nm,其總孔容和BET比表面積相對較大,分別為(4.46~40.84)×10-3cm3/g 和1.16~38.14 m2/g。CO2吸附試驗表明,低階煤煤樣≤2 nm 微孔孔容和比表面積相對較大,低階煤的微孔孔容為0.055~0.064 cm3/g,總體上高于中階煤的微孔孔容(0.026~0.060 cm3/g)。隨著煤化程度的提高,試驗低階煤樣的微孔和介孔孔容均呈現(xiàn)波動性減小的趨勢。

    (2)SAXS試驗結(jié)果表明,低階煤樣最發(fā)育孔徑主要分布在10 nm。這與N2吸附法測定結(jié)果(最發(fā)育孔徑為3 nm左右)存在較大差異,可能是煤中存在閉孔或氣體吸附法無法測到的孔。

    (3)基于N2吸附的分形維數(shù)D1和D2分別為2.44~2.92和2.03~2.61。除了DX04煤樣,低階煤的D2值(2.25~2.61)均大于中階煤的D2值(2.04和2.10),表明低階煤<5 nm的吸附孔隙較為發(fā)育。而分形維數(shù)D1和變質(zhì)程度沒有明顯規(guī)律,可能煤階不是影響滲流孔隙的主要因素?;赟AXS試驗結(jié)果測得孔隙分形維數(shù)DS為2.41~2.72,低階煤的DS值(2.44~2.72)均大于中階煤XY06的DS值(2.41),表明低階煤樣的孔隙結(jié)構(gòu)更為復(fù)雜。

    (4)甲烷吸附/解吸試驗表明,除DX04煤樣外,低階煤的VL值為31.14~35.14 m3/t,大于中階煤的VL值(15.80~28.21 m3/t)。表明低階煤樣瓦斯吸附/解吸能力總體上強(qiáng)于中階煤樣。低階煤的微孔和介孔較為發(fā)育,微孔孔容與VL存在正相關(guān)關(guān)系,介孔孔容與擴(kuò)散系數(shù)(D)存在正相關(guān)關(guān)系。而VL和D反映煤體的吸附、放散特性,因此低階煤具有較強(qiáng)的瓦斯吸附和放散能力。

    (5)低階煤分子結(jié)構(gòu)松散,含有較多長側(cè)鏈的官能團(tuán),微孔和介孔發(fā)育較好,瓦斯吸附和擴(kuò)散能力較強(qiáng)。低階煤具有初始瓦斯解吸速度快特點,雖有利于煤層氣的開發(fā),但是受采動影響下構(gòu)造賦存區(qū)的低階煤層易于發(fā)生瓦斯異常涌出和瓦斯突出。

    猜你喜歡
    孔容低階煤樣
    一種微球狀流化床催化劑及其制備方法和應(yīng)用
    能源化工(2022年3期)2023-01-15 02:26:43
    山西低階煤分布特征分析和開發(fā)利用前景
    預(yù)氧化對煤比表面積及孔徑分布影響的實驗研究
    煤礦安全(2020年10期)2020-11-02 08:27:50
    勘 誤
    一類具低階項和退化強(qiáng)制的橢圓方程的有界弱解
    突出煤樣孔徑分布研究
    煤(2019年4期)2019-04-28 06:51:38
    Extended Fisher-Kolmogorov方程的一類低階非協(xié)調(diào)混合有限元方法
    談?wù)劽簶又苽渲袘?yīng)注意的問題
    帶旋轉(zhuǎn)孔容腔瞬態(tài)演化與建模方法研究
    桃園井田煤的孔隙結(jié)構(gòu)特征研究
    国产av又大| 日韩中文字幕视频在线看片| 成人三级做爰电影| 国产在线精品亚洲第一网站| 国产精品99久久99久久久不卡| 亚洲专区中文字幕在线| 宅男免费午夜| 久久精品国产亚洲av高清一级| 国产欧美日韩综合在线一区二区| 天堂8中文在线网| 亚洲精品一卡2卡三卡4卡5卡| 亚洲成人免费av在线播放| 国产在线免费精品| 丝袜人妻中文字幕| 国产成人啪精品午夜网站| 国产亚洲精品一区二区www | 黑人巨大精品欧美一区二区蜜桃| 中文字幕人妻丝袜一区二区| 欧美精品一区二区大全| 国产单亲对白刺激| 美国免费a级毛片| 欧美+亚洲+日韩+国产| 18禁国产床啪视频网站| 久久天躁狠狠躁夜夜2o2o| 在线av久久热| 啦啦啦中文免费视频观看日本| 亚洲成人免费电影在线观看| 日韩欧美一区二区三区在线观看 | 色精品久久人妻99蜜桃| 大片电影免费在线观看免费| 咕卡用的链子| 久久精品国产a三级三级三级| 99在线人妻在线中文字幕 | 久久人妻av系列| 久热爱精品视频在线9| 精品一区二区三卡| 亚洲av日韩在线播放| 一本色道久久久久久精品综合| 黄色毛片三级朝国网站| 12—13女人毛片做爰片一| 欧美av亚洲av综合av国产av| 一本一本久久a久久精品综合妖精| 这个男人来自地球电影免费观看| 国产人伦9x9x在线观看| 亚洲欧美日韩另类电影网站| 80岁老熟妇乱子伦牲交| av超薄肉色丝袜交足视频| 国产亚洲欧美精品永久| 亚洲va日本ⅴa欧美va伊人久久| 丁香六月天网| 丝袜美腿诱惑在线| 国产成人免费观看mmmm| 王馨瑶露胸无遮挡在线观看| av线在线观看网站| 国产男女超爽视频在线观看| 三级毛片av免费| 啦啦啦中文免费视频观看日本| 国产成人欧美| 欧美激情极品国产一区二区三区| 最黄视频免费看| 欧美成人免费av一区二区三区 | 国产精品亚洲一级av第二区| 亚洲视频免费观看视频| 日韩大码丰满熟妇| 精品第一国产精品| 国产一区二区三区视频了| 人人妻人人爽人人添夜夜欢视频| 中文欧美无线码| 欧美久久黑人一区二区| 黄色成人免费大全| 国产不卡av网站在线观看| 日韩人妻精品一区2区三区| 国产亚洲午夜精品一区二区久久| 国产精品免费视频内射| 宅男免费午夜| 波多野结衣一区麻豆| 亚洲精品久久午夜乱码| 亚洲专区中文字幕在线| 亚洲人成77777在线视频| 国产精品秋霞免费鲁丝片| 一区二区三区国产精品乱码| 精品亚洲乱码少妇综合久久| 国产有黄有色有爽视频| 大码成人一级视频| 黄色视频在线播放观看不卡| 少妇被粗大的猛进出69影院| 国产黄色免费在线视频| 50天的宝宝边吃奶边哭怎么回事| 每晚都被弄得嗷嗷叫到高潮| 久久精品人人爽人人爽视色| 色婷婷av一区二区三区视频| 免费在线观看视频国产中文字幕亚洲| 少妇精品久久久久久久| 精品免费久久久久久久清纯 | 手机成人av网站| 精品人妻在线不人妻| 色精品久久人妻99蜜桃| 9色porny在线观看| 在线亚洲精品国产二区图片欧美| 成人亚洲精品一区在线观看| 国产单亲对白刺激| 国产精品一区二区精品视频观看| 亚洲精品一二三| 成人手机av| 亚洲 欧美一区二区三区| 天天躁夜夜躁狠狠躁躁| 欧美日韩亚洲综合一区二区三区_| 欧美黄色淫秽网站| 老熟妇乱子伦视频在线观看| 亚洲精品一二三| 如日韩欧美国产精品一区二区三区| 一级毛片女人18水好多| 日韩欧美国产一区二区入口| 免费观看av网站的网址| 午夜久久久在线观看| 午夜两性在线视频| 国产一区二区三区综合在线观看| 亚洲人成伊人成综合网2020| 国产精品熟女久久久久浪| 午夜福利乱码中文字幕| 久久精品亚洲av国产电影网| 日日夜夜操网爽| 亚洲精品在线观看二区| 纵有疾风起免费观看全集完整版| 久久国产精品人妻蜜桃| 999久久久国产精品视频| 精品久久久久久久毛片微露脸| 最新在线观看一区二区三区| 手机成人av网站| 黄片小视频在线播放| 欧美人与性动交α欧美精品济南到| 真人做人爱边吃奶动态| tube8黄色片| 久久久久久久久久久久大奶| 国产精品一区二区免费欧美| 亚洲精品在线观看二区| 天天操日日干夜夜撸| 伦理电影免费视频| 午夜精品国产一区二区电影| 99国产综合亚洲精品| 国产亚洲欧美精品永久| 亚洲精品在线美女| 国产成人精品在线电影| 高清av免费在线| 亚洲专区国产一区二区| 亚洲久久久国产精品| 日本wwww免费看| 欧美 亚洲 国产 日韩一| 精品一品国产午夜福利视频| 免费人妻精品一区二区三区视频| 精品少妇内射三级| 国产精品麻豆人妻色哟哟久久| 欧美精品一区二区大全| 精品国产乱子伦一区二区三区| 12—13女人毛片做爰片一| 丝袜美腿诱惑在线| 国产片内射在线| 飞空精品影院首页| 国产亚洲精品久久久久5区| 男女高潮啪啪啪动态图| 18禁观看日本| 亚洲国产精品一区二区三区在线| 久久人人爽av亚洲精品天堂| 精品亚洲乱码少妇综合久久| 亚洲精华国产精华精| xxxhd国产人妻xxx| 无遮挡黄片免费观看| 男女边摸边吃奶| 啦啦啦视频在线资源免费观看| tube8黄色片| 最新的欧美精品一区二区| 国产av国产精品国产| 亚洲美女黄片视频| 最近最新免费中文字幕在线| 日日爽夜夜爽网站| 韩国精品一区二区三区| 中文字幕人妻丝袜一区二区| 女人精品久久久久毛片| 热re99久久国产66热| 欧美 日韩 精品 国产| 日韩欧美免费精品| 久久午夜综合久久蜜桃| 亚洲少妇的诱惑av| 国产黄色免费在线视频| 精品国产国语对白av| 久久精品国产亚洲av香蕉五月 | 窝窝影院91人妻| 欧美av亚洲av综合av国产av| 91国产中文字幕| 久久人妻av系列| 国产无遮挡羞羞视频在线观看| 啦啦啦在线免费观看视频4| 国产无遮挡羞羞视频在线观看| 日韩三级视频一区二区三区| 后天国语完整版免费观看| 久久ye,这里只有精品| 另类精品久久| 十八禁网站免费在线| 精品久久久久久电影网| 日韩免费av在线播放| 久久中文看片网| 黄片大片在线免费观看| 国产精品亚洲av一区麻豆| 制服人妻中文乱码| 女性被躁到高潮视频| 亚洲欧美色中文字幕在线| 久久久久精品人妻al黑| 一级,二级,三级黄色视频| 两人在一起打扑克的视频| 人人澡人人妻人| 侵犯人妻中文字幕一二三四区| 丝袜喷水一区| 久久国产精品大桥未久av| 欧美成狂野欧美在线观看| 男女高潮啪啪啪动态图| 中文字幕高清在线视频| 在线观看www视频免费| 日韩成人在线观看一区二区三区| 午夜福利免费观看在线| 两性夫妻黄色片| 法律面前人人平等表现在哪些方面| 一边摸一边抽搐一进一出视频| 日韩视频在线欧美| 欧美日韩福利视频一区二区| 欧美日韩精品网址| 老司机影院毛片| 亚洲欧洲精品一区二区精品久久久| 乱人伦中国视频| 一边摸一边抽搐一进一小说 | 亚洲欧美精品综合一区二区三区| 一个人免费在线观看的高清视频| 国产精品一区二区在线不卡| 国产野战对白在线观看| 国产av又大| 女性被躁到高潮视频| 岛国毛片在线播放| svipshipincom国产片| 老汉色∧v一级毛片| 女人久久www免费人成看片| 欧美在线一区亚洲| 丝袜美腿诱惑在线| 日本一区二区免费在线视频| 美女扒开内裤让男人捅视频| 国产亚洲av高清不卡| 欧美黑人精品巨大| 国产黄频视频在线观看| 黑人猛操日本美女一级片| 丝袜在线中文字幕| 久久中文字幕人妻熟女| www.熟女人妻精品国产| 亚洲国产成人一精品久久久| 免费观看av网站的网址| 亚洲精品国产色婷婷电影| 免费在线观看影片大全网站| 精品久久久久久久毛片微露脸| 欧美成狂野欧美在线观看| av天堂久久9| 一区二区三区激情视频| 国产精品香港三级国产av潘金莲| 久久天躁狠狠躁夜夜2o2o| 久久久久精品国产欧美久久久| 欧美日韩av久久| 欧美精品高潮呻吟av久久| 国产亚洲精品久久久久5区| 精品国产一区二区久久| 天堂8中文在线网| 在线观看www视频免费| 亚洲av美国av| 两性夫妻黄色片| 最近最新中文字幕大全电影3 | 亚洲国产欧美一区二区综合| 黑人猛操日本美女一级片| 精品少妇一区二区三区视频日本电影| 日韩欧美一区二区三区在线观看 | 日本av免费视频播放| 久久久久网色| 少妇裸体淫交视频免费看高清 | 免费在线观看黄色视频的| 在线观看66精品国产| 性色av乱码一区二区三区2| 国产精品电影一区二区三区 | 十八禁人妻一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 777米奇影视久久| 好男人电影高清在线观看| 国产黄频视频在线观看| 一级a爱视频在线免费观看| 免费观看人在逋| 国产av国产精品国产| 宅男免费午夜| www.精华液| 丝袜喷水一区| 色婷婷av一区二区三区视频| 少妇粗大呻吟视频| 欧美日韩成人在线一区二区| 国产精品1区2区在线观看. | 久久午夜亚洲精品久久| 麻豆成人av在线观看| 在线看a的网站| 9191精品国产免费久久| 99国产精品免费福利视频| 多毛熟女@视频| 无遮挡黄片免费观看| 人人妻人人澡人人看| 中文亚洲av片在线观看爽 | 国产精品免费大片| 国产精品麻豆人妻色哟哟久久| 后天国语完整版免费观看| 亚洲欧美一区二区三区久久| 精品国产一区二区三区四区第35| 中国美女看黄片| 91精品国产国语对白视频| 精品一区二区三区视频在线观看免费 | 51午夜福利影视在线观看| 欧美久久黑人一区二区| 十八禁高潮呻吟视频| 久久免费观看电影| 99热国产这里只有精品6| 汤姆久久久久久久影院中文字幕| 国产高清国产精品国产三级| 午夜精品国产一区二区电影| a级毛片在线看网站| 日韩一区二区三区影片| 人妻 亚洲 视频| 精品久久久久久电影网| 在线观看免费视频日本深夜| 精品国产一区二区三区久久久樱花| 国产高清国产精品国产三级| 久久人人97超碰香蕉20202| 久久精品国产亚洲av高清一级| 我的亚洲天堂| 91成人精品电影| 男人操女人黄网站| 岛国在线观看网站| 美女高潮到喷水免费观看| 法律面前人人平等表现在哪些方面| 天天躁夜夜躁狠狠躁躁| 麻豆av在线久日| 咕卡用的链子| 亚洲久久久国产精品| 亚洲avbb在线观看| 自线自在国产av| 欧美日韩亚洲综合一区二区三区_| 精品第一国产精品| 涩涩av久久男人的天堂| 久久亚洲真实| 亚洲欧美一区二区三区久久| 国产91精品成人一区二区三区 | 亚洲五月色婷婷综合| 啦啦啦 在线观看视频| 欧美国产精品一级二级三级| 波多野结衣av一区二区av| 中文字幕人妻熟女乱码| 久久99热这里只频精品6学生| 一个人免费看片子| 一个人免费在线观看的高清视频| 午夜成年电影在线免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产熟女午夜一区二区三区| 97人妻天天添夜夜摸| 成在线人永久免费视频| 精品一区二区三区av网在线观看 | 国产精品久久久久久人妻精品电影 | 国产在线一区二区三区精| 久久性视频一级片| 久久精品亚洲熟妇少妇任你| 啦啦啦 在线观看视频| 亚洲综合色网址| 飞空精品影院首页| 久久久久久久久久久久大奶| 一本综合久久免费| 欧美日本中文国产一区发布| 欧美国产精品一级二级三级| av线在线观看网站| 丝袜人妻中文字幕| 女人爽到高潮嗷嗷叫在线视频| 午夜日韩欧美国产| 不卡一级毛片| 中文字幕人妻丝袜制服| av一本久久久久| 欧美精品啪啪一区二区三区| 亚洲色图av天堂| 精品视频人人做人人爽| 香蕉国产在线看| 久久精品国产综合久久久| 1024香蕉在线观看| 精品熟女少妇八av免费久了| 老熟妇乱子伦视频在线观看| 亚洲中文av在线| av一本久久久久| 男女午夜视频在线观看| 国产成人影院久久av| 欧美人与性动交α欧美软件| 99精品欧美一区二区三区四区| av又黄又爽大尺度在线免费看| 亚洲国产欧美日韩在线播放| 国产伦人伦偷精品视频| 午夜免费鲁丝| 国产淫语在线视频| 午夜老司机福利片| 9热在线视频观看99| 91成年电影在线观看| 搡老乐熟女国产| 亚洲国产毛片av蜜桃av| 成年人黄色毛片网站| 无人区码免费观看不卡 | 波多野结衣av一区二区av| 嫁个100分男人电影在线观看| netflix在线观看网站| 日本黄色日本黄色录像| 精品欧美一区二区三区在线| 久久婷婷成人综合色麻豆| 91av网站免费观看| 精品国产乱子伦一区二区三区| 久久精品熟女亚洲av麻豆精品| 国产极品粉嫩免费观看在线| 精品福利观看| 日韩熟女老妇一区二区性免费视频| 亚洲欧洲精品一区二区精品久久久| 99国产极品粉嫩在线观看| 黄片大片在线免费观看| 欧美人与性动交α欧美软件| 美女午夜性视频免费| 亚洲人成伊人成综合网2020| 天堂俺去俺来也www色官网| 日本a在线网址| 欧美日韩国产mv在线观看视频| 日本精品一区二区三区蜜桃| 欧美中文综合在线视频| 女人高潮潮喷娇喘18禁视频| 成年人黄色毛片网站| 九色亚洲精品在线播放| 欧美亚洲 丝袜 人妻 在线| 十分钟在线观看高清视频www| 伊人久久大香线蕉亚洲五| 亚洲国产av影院在线观看| 国产精品麻豆人妻色哟哟久久| 法律面前人人平等表现在哪些方面| 国产一卡二卡三卡精品| 啦啦啦中文免费视频观看日本| 久久久国产精品麻豆| 国产在视频线精品| 午夜福利在线观看吧| 日本黄色视频三级网站网址 | 国产一区二区 视频在线| 亚洲熟妇熟女久久| 国产成+人综合+亚洲专区| 天天操日日干夜夜撸| 日韩中文字幕欧美一区二区| 女警被强在线播放| 少妇精品久久久久久久| svipshipincom国产片| 国产精品二区激情视频| 男男h啪啪无遮挡| 天堂中文最新版在线下载| 亚洲专区国产一区二区| 中文字幕人妻熟女乱码| 欧美精品高潮呻吟av久久| 一级毛片电影观看| 国产欧美亚洲国产| 无限看片的www在线观看| av超薄肉色丝袜交足视频| 欧美日韩黄片免| 一区二区av电影网| 日韩视频在线欧美| 人人澡人人妻人| 露出奶头的视频| 成年人午夜在线观看视频| 国产黄色免费在线视频| 久久人人97超碰香蕉20202| videosex国产| av不卡在线播放| 这个男人来自地球电影免费观看| 在线观看www视频免费| 无人区码免费观看不卡 | 国产日韩一区二区三区精品不卡| 亚洲久久久国产精品| netflix在线观看网站| cao死你这个sao货| 黄色 视频免费看| 操美女的视频在线观看| 人妻 亚洲 视频| 啦啦啦中文免费视频观看日本| 99国产综合亚洲精品| www日本在线高清视频| 亚洲av国产av综合av卡| 欧美+亚洲+日韩+国产| 亚洲综合色网址| 高清视频免费观看一区二区| 久久精品国产综合久久久| 久久久久久久久久久久大奶| 中文字幕制服av| 欧美一级毛片孕妇| 99re6热这里在线精品视频| 久久精品国产综合久久久| 老鸭窝网址在线观看| 欧美 日韩 精品 国产| 亚洲精品乱久久久久久| 久久久欧美国产精品| 日日摸夜夜添夜夜添小说| 中文字幕人妻丝袜制服| 天堂俺去俺来也www色官网| 亚洲人成伊人成综合网2020| 国产高清激情床上av| 午夜免费鲁丝| 亚洲欧美一区二区三区黑人| 超碰97精品在线观看| 久久久久久亚洲精品国产蜜桃av| 宅男免费午夜| 91九色精品人成在线观看| 麻豆av在线久日| 好男人电影高清在线观看| 2018国产大陆天天弄谢| 日韩三级视频一区二区三区| 757午夜福利合集在线观看| 91国产中文字幕| 精品亚洲成国产av| 在线播放国产精品三级| 国产一区二区三区在线臀色熟女 | 亚洲专区字幕在线| 国产精品成人在线| 少妇精品久久久久久久| 国产成人av激情在线播放| 免费观看人在逋| 激情视频va一区二区三区| 十分钟在线观看高清视频www| 亚洲精品一二三| 亚洲欧美一区二区三区久久| 十八禁人妻一区二区| 女人高潮潮喷娇喘18禁视频| 最新美女视频免费是黄的| 成在线人永久免费视频| 黄网站色视频无遮挡免费观看| av片东京热男人的天堂| 精品乱码久久久久久99久播| 一本大道久久a久久精品| 大香蕉久久成人网| 亚洲国产毛片av蜜桃av| 久久99热这里只频精品6学生| 亚洲国产精品一区二区三区在线| 久久婷婷成人综合色麻豆| 自拍欧美九色日韩亚洲蝌蚪91| 老汉色av国产亚洲站长工具| 亚洲午夜精品一区,二区,三区| 欧美日韩福利视频一区二区| 欧美黑人精品巨大| 大片免费播放器 马上看| 亚洲精品自拍成人| 黑人巨大精品欧美一区二区蜜桃| 亚洲欧洲精品一区二区精品久久久| 午夜福利在线免费观看网站| 国产日韩欧美视频二区| 国产精品亚洲av一区麻豆| 男女床上黄色一级片免费看| www日本在线高清视频| 最黄视频免费看| 久热这里只有精品99| 亚洲国产欧美网| 十八禁人妻一区二区| 丝袜喷水一区| 精品国内亚洲2022精品成人 | 亚洲色图 男人天堂 中文字幕| 巨乳人妻的诱惑在线观看| 啦啦啦中文免费视频观看日本| 交换朋友夫妻互换小说| 亚洲av成人不卡在线观看播放网| 亚洲精品国产色婷婷电影| 欧美日韩亚洲综合一区二区三区_| 妹子高潮喷水视频| 欧美另类亚洲清纯唯美| 露出奶头的视频| 日本黄色视频三级网站网址 | 手机成人av网站| 精品国产一区二区久久| 国产精品免费一区二区三区在线 | 亚洲三区欧美一区| 久久精品国产99精品国产亚洲性色 | 又大又爽又粗| 一个人免费看片子| 国精品久久久久久国模美| 中亚洲国语对白在线视频| 日本五十路高清| 淫妇啪啪啪对白视频| 亚洲七黄色美女视频| 亚洲 国产 在线| 搡老乐熟女国产| 日本a在线网址| 好男人电影高清在线观看| 久久午夜综合久久蜜桃| 51午夜福利影视在线观看| 午夜精品国产一区二区电影| 久久久久精品国产欧美久久久| 中文字幕另类日韩欧美亚洲嫩草| 成人av一区二区三区在线看| av一本久久久久| 免费观看av网站的网址| 久久久久久久国产电影| 高清在线国产一区| 一夜夜www| 久久久久久久国产电影| 制服人妻中文乱码| 国产成人精品无人区| 99国产极品粉嫩在线观看| 午夜免费鲁丝| 一级,二级,三级黄色视频| 波多野结衣av一区二区av| 啦啦啦视频在线资源免费观看| 99riav亚洲国产免费| 美女午夜性视频免费| 9191精品国产免费久久| 欧美日韩视频精品一区| 人妻 亚洲 视频| 欧美人与性动交α欧美软件| 国产国语露脸激情在线看| 午夜福利,免费看|