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

    構(gòu)造煤與原生結(jié)構(gòu)煤中甲烷擴(kuò)散性能差異性分析

    2022-04-12 07:08:02孟召平
    煤田地質(zhì)與勘探 2022年3期
    關(guān)鍵詞:擴(kuò)散系數(shù)煤體微孔

    孟召平,張 昆,沈 振

    (1.中國(guó)礦業(yè)大學(xué)(北京) 地球科學(xué)與測(cè)繪工程學(xué)院,北京 100083;2.煤與煤層氣共采國(guó)家重點(diǎn)實(shí)驗(yàn)室,山西 晉城 048012)

    我國(guó)自20 世紀(jì)80 年代以來(lái),采用美國(guó)以煤層氣井水力壓裂和排水降壓采氣為特征的地面煤層氣開發(fā)理論與技術(shù),對(duì)我國(guó)煤層氣資源進(jìn)行了勘探開發(fā)研究[1]。截至2020 年底,在我國(guó)施工各類煤層氣井達(dá)18 000 余口,煤層氣地面開發(fā)產(chǎn)量63.7 億m3、利用量59.4 億m3,我國(guó)煤層氣產(chǎn)業(yè)發(fā)展已進(jìn)入規(guī)?;a(chǎn)階段,初步形成了適宜于原生結(jié)構(gòu)煤的煤層氣勘探開發(fā)技術(shù)體系,但是,在構(gòu)造煤區(qū)煤層氣開發(fā)地質(zhì)理論與技術(shù)研究相對(duì)薄弱,在構(gòu)造煤區(qū)煤層氣開發(fā)仍未取得實(shí)質(zhì)性突破,特別是對(duì)構(gòu)造煤中氣體吸附-解吸和擴(kuò)散-滲流機(jī)理的認(rèn)識(shí)仍然不足,導(dǎo)致煤層氣井排采控制不合理,地面煤層氣井普遍產(chǎn)量低、不穩(wěn)定,現(xiàn)有技術(shù)難以支撐產(chǎn)業(yè)快速發(fā)展的問(wèn)題亟待解決。因此,開展構(gòu)造煤與原生結(jié)構(gòu)煤中甲烷擴(kuò)散性能差異性分析,對(duì)于合理有效開發(fā)我國(guó)構(gòu)造煤區(qū)煤層氣資源具有理論和實(shí)際應(yīng)用意義。

    我國(guó)大多數(shù)沉積盆地在成煤過(guò)程中或成煤期后煤層受到多期次構(gòu)造運(yùn)動(dòng)的作用,造成原生結(jié)構(gòu)煤破壞嚴(yán)重,形成了不同結(jié)構(gòu)類型構(gòu)造煤[1-3]。不同類型的構(gòu)造煤在孔隙結(jié)構(gòu)、構(gòu)造特征、煤物理力學(xué)性質(zhì)、吸附/解吸、擴(kuò)散性能和滲透性等方面明顯地不同于原生結(jié)構(gòu)煤,且存在較大的差異性[4-8]。

    對(duì)于不同煤體結(jié)構(gòu)煤中甲烷吸附/解吸與擴(kuò)散性能,國(guó)內(nèi)外學(xué)者進(jìn)行了大量的實(shí)驗(yàn)研究,揭示了煤層遭受構(gòu)造變形破壞后其煤中氣體吸附/解吸能力受到煤的孔隙結(jié)構(gòu)、溫度、壓力和物理化學(xué)性質(zhì)等多種因素影響規(guī)律[4-10]。如孟召平等[4]通過(guò)對(duì)不同煤體結(jié)構(gòu)煤樣進(jìn)行液氮吸附和甲烷等溫吸附實(shí)驗(yàn),認(rèn)為隨著煤體結(jié)構(gòu)破壞程度的增加,煤中微孔比表面積的升高是造成構(gòu)造煤吸附性增強(qiáng)的主要原因。陳亮等[5]通過(guò)實(shí)驗(yàn)認(rèn)為,微孔對(duì)于比表面積的貢獻(xiàn)率可高達(dá)71%,控制著煤體吸附氣含量。簡(jiǎn)闊[6]、張小東[7]等從溫度、分子結(jié)構(gòu)和煤巖煤質(zhì)的變化角度分析了構(gòu)造煤和原生結(jié)構(gòu)煤吸附甲烷能力的差異性及控制機(jī)理。姜波[8]、屈爭(zhēng)輝[9]等揭示了構(gòu)造煤結(jié)構(gòu)及其煤儲(chǔ)層物性特征以及構(gòu)造煤結(jié)構(gòu)對(duì)瓦斯特性的控制機(jī)理。在此基礎(chǔ)上,國(guó)內(nèi)外學(xué)者進(jìn)一步探討了煤中氣體擴(kuò)散性能及其控制機(jī)理。孟召平等[10]分析了低煤階煤中甲烷吸附、擴(kuò)散能力主要受溫度和壓力條件影響機(jī)制。M.Pillalamarry 等[11]認(rèn)為:構(gòu)造煤初始瓦斯放散速率明顯高于原生結(jié)構(gòu)煤,氣體分子在煤孔壁上的解吸是瞬間完成的,煤中氣體在孔隙中的運(yùn)移需要克服一定的阻力,研究發(fā)現(xiàn),菲克定律能夠很好地描述甲烷在煤等多孔介質(zhì)中的運(yùn)移規(guī)律。單孔擴(kuò)散模型和雙孔擴(kuò)散模型的提出,可以通過(guò)擴(kuò)散系數(shù)定量表征甲烷在孔隙中的運(yùn)移特性[12-14];盡管此2 種模型在邊界條件設(shè)定及使用對(duì)象上仍存在一定偏差,但是可以用來(lái)分析煤中氣體擴(kuò)散能力的大小。構(gòu)造煤中氣體的初始擴(kuò)散系數(shù)要大于原生結(jié)構(gòu)煤,其主要是由于構(gòu)造煤孔隙結(jié)構(gòu)發(fā)育,大孔、中孔和過(guò)渡孔的容積遠(yuǎn)大于原生結(jié)構(gòu)所致[15]。由此可以看出,原生結(jié)構(gòu)煤與構(gòu)造煤中甲烷的吸附/解吸和擴(kuò)散能力將發(fā)生明顯變化,一定程度上影響煤層氣的富集與產(chǎn)出規(guī)律。

    筆者通過(guò)采取沁水盆地東部趙莊井田3 號(hào)煤層原生結(jié)構(gòu)煤與構(gòu)造煤樣品進(jìn)行甲烷等溫吸附實(shí)驗(yàn)、低溫液氮和二氧化碳吸附實(shí)驗(yàn),揭示構(gòu)造煤與原生結(jié)構(gòu)煤的吸附-解吸性能及孔隙結(jié)構(gòu)特征;應(yīng)用單孔和雙孔非穩(wěn)態(tài)擴(kuò)散模型,揭示原生結(jié)構(gòu)煤與構(gòu)造煤中甲烷擴(kuò)散差異性及控制機(jī)理,為我國(guó)構(gòu)造煤區(qū)煤層氣勘探開發(fā)提供理論依據(jù)。

    1 實(shí)驗(yàn)方法與計(jì)算模型

    1.1 樣品選擇及預(yù)處理

    原生結(jié)構(gòu)煤和構(gòu)造煤樣品來(lái)自沁水盆地趙莊井田同一口煤層氣開發(fā)井,構(gòu)造煤包括碎裂結(jié)構(gòu)煤、碎粒結(jié)構(gòu)煤和糜棱結(jié)構(gòu)煤[1,4]。目標(biāo)煤層為二疊系下統(tǒng)山西組3 號(hào)煤層,埋深在844.90~850.75 m。樣品經(jīng)過(guò)繩索取心到達(dá)地面后,立即封入解吸罐中,按照USBM 直接法測(cè)量含氣量。煤巖類型以半亮煤、光亮煤為主。煤樣類型與基礎(chǔ)數(shù)據(jù)見表1。

    表1 樣品工業(yè)分析、顯微組分及煤體結(jié)構(gòu)劃分Table 1 Proximate analysis, macerals and coal structure classification of samples

    不同煤體結(jié)構(gòu)之間煤巖煤質(zhì)及變質(zhì)程度的差異性較小,因此,在研究不同煤體結(jié)構(gòu)的吸附/解吸與擴(kuò)散特性時(shí)可排除這些影響因素。不同煤體結(jié)構(gòu)含氣量變化范圍為10.41~15.45 cm3/g,且含氣量隨著煤破壞程度的增大而呈現(xiàn)升高的趨勢(shì),說(shuō)明構(gòu)造煤中含氣量要高于原生結(jié)構(gòu)煤。

    1.2 等溫吸附/解吸實(shí)驗(yàn)

    為了研究原生結(jié)構(gòu)煤和構(gòu)造煤的吸附/解吸與擴(kuò)散性能,本次實(shí)驗(yàn)研究采用Terra-Tek 公司生產(chǎn)的ISO-300 等溫吸附/解吸實(shí)驗(yàn)儀器,按照高壓容量法進(jìn)行甲烷等溫吸附實(shí)驗(yàn)。按照GB/T 19560-2008《煤的高壓等溫吸附試驗(yàn)方法》將樣品粉碎研磨至60~80 目(0.2~0.3 mm),各樣品質(zhì)量100~120 g,整個(gè)過(guò)程中保持在恒溫油浴中(25、35、45℃),實(shí)驗(yàn)壓力在1~10 MPa。由于篩分過(guò)的煤顆粒的尺寸比煤基質(zhì)的孔徑大多個(gè)數(shù)量級(jí),對(duì)煤中孔隙影響較小[12],因此,忽略樣品制備過(guò)程對(duì)煤基質(zhì)孔隙結(jié)構(gòu)的影響。

    1.3 擴(kuò)散計(jì)算模型

    擴(kuò)散過(guò)程主要發(fā)生在煤基質(zhì)孔隙系統(tǒng),當(dāng)儲(chǔ)層壓力低于臨界解吸壓力時(shí),煤孔隙表面吸附態(tài)CH4分子隨即脫附成游離態(tài),并在濃度差的驅(qū)動(dòng)下由孔隙向裂隙運(yùn)移。由于煤孔隙結(jié)構(gòu)的復(fù)雜性和CH4氣體分子的不規(guī)則熱運(yùn)動(dòng),導(dǎo)致擴(kuò)散是一種效率極低的物質(zhì)運(yùn)移方式。目前,為了對(duì)煤層甲烷擴(kuò)散效應(yīng)進(jìn)行定量描述,通過(guò)采用等溫吸附實(shí)驗(yàn)記錄的甲烷吸附/解吸過(guò)程中氣體壓力隨時(shí)間的變化數(shù)據(jù),結(jié)合單孔或者雙孔數(shù)學(xué)擴(kuò)散模型表征氣體在多孔介質(zhì)中擴(kuò)散性能的大小。單孔模型和雙孔模型都是基于菲克第二定律的非穩(wěn)態(tài)擴(kuò)散模型,考慮了擴(kuò)散過(guò)程中氣體濃度隨時(shí)間的變化關(guān)系:

    式中:D為擴(kuò)散系數(shù),m2/s;r為均質(zhì)球粒半徑,m;c為CH4氣體質(zhì)量濃度,g/cm3;t為擴(kuò)散時(shí)間,s;D/r2為有效擴(kuò)散系數(shù),s-1。

    理想條件下,式(1)的解析解可表示為:

    式中:Vt為在時(shí)間t內(nèi)擴(kuò)散氣體的總量;V∞為無(wú)限時(shí)間內(nèi)的總解吸氣體量;ra為擴(kuò)散路徑長(zhǎng)度。

    式(2)是經(jīng)典的“單孔模型”,一些研究人員已將其應(yīng)用于擬合吸附數(shù)據(jù)以評(píng)估擴(kuò)散率[12],且可以相對(duì)準(zhǔn)確地評(píng)估擴(kuò)散過(guò)程。

    煤作為多孔介質(zhì)固體,其孔隙結(jié)構(gòu)非均質(zhì)性較強(qiáng),在單孔擴(kuò)散模型中的一個(gè)必要條件是不考慮孔隙結(jié)構(gòu)的非均質(zhì)分布,假定煤中甲烷擴(kuò)散邊界條件為均質(zhì)球形和均勻的孔隙分布特征。但是,研究表明煤中孔隙結(jié)構(gòu)分布復(fù)雜性較強(qiáng),通常呈現(xiàn)多峰分布,此時(shí),單孔擴(kuò)散模型就不能較為準(zhǔn)確地模擬甲烷的擴(kuò)散過(guò)程。鑒于單孔擴(kuò)散模型的局限性,許多學(xué)者在單孔擴(kuò)散模型的基礎(chǔ)上提出了可以有效模擬大孔和微孔分布的雙孔擴(kuò)散模型(圖1)。

    圖1 雙孔擴(kuò)散模型概念圖[16]Fig.1 Conceptual diagram of the bidisperse model[16]

    雙孔擴(kuò)散模型模擬儲(chǔ)層擴(kuò)散過(guò)程,將基質(zhì)孔隙分為一個(gè)快速擴(kuò)散的大孔階段和一個(gè)擴(kuò)散極慢的微孔擴(kuò)散階段[17]。對(duì)于大孔擴(kuò)散階段,其數(shù)學(xué)表達(dá)式為:

    微孔擴(kuò)散階段表達(dá)式:

    式中:Va為在時(shí)間t內(nèi)大孔中氣體吸附/解吸的總量,cm3/g;Vi為在時(shí)間t內(nèi)微孔中氣體吸附/解吸的總量,cm3/g;Da為大孔擴(kuò)散系數(shù),m2/s;Di為微孔擴(kuò)散系數(shù),m2/s;Va,∞和Vi,∞分別為在平衡壓力點(diǎn)大孔和微孔中的最大氣體吸附量,cm3/g;ri為微孔半徑,m。

    根據(jù)大孔和微孔在整體擴(kuò)散過(guò)程中對(duì)于甲烷運(yùn)移效率的貢獻(xiàn)程度,Pan Zhejun 等(2010)[17]提出比例系數(shù)β來(lái)表征大孔和微孔在整個(gè)擴(kuò)散階段的貢獻(xiàn)率,其數(shù)學(xué)表達(dá)式如下:

    2 實(shí)驗(yàn)結(jié)果和分析

    2.1 構(gòu)造煤與原生結(jié)構(gòu)煤等溫吸附/解吸實(shí)驗(yàn)結(jié)果

    25℃條件下等溫吸附實(shí)驗(yàn)結(jié)果顯示,甲烷在不同煤體結(jié)構(gòu)煤中的吸附/解吸規(guī)律符合Langmuir 模型(圖2)。

    為了排除不同煤體結(jié)構(gòu)中水分和灰分產(chǎn)率的差異性對(duì)煤吸附性影響,此次研究中選取干燥無(wú)灰基數(shù)據(jù)進(jìn)行分析研究。4 種煤體結(jié)構(gòu)Langmuir 體積VL較高,Langmuir 壓力pL相對(duì)較低,說(shuō)明4 種煤體結(jié)構(gòu)煤有較強(qiáng)的吸附能力。從原生結(jié)構(gòu)煤到糜棱結(jié)構(gòu)煤,隨著煤體破碎程度的增加,煤的吸附性增強(qiáng),順序?yàn)椋涸Y(jié)構(gòu)<碎裂結(jié)構(gòu)<碎粒結(jié)構(gòu)<糜棱結(jié)構(gòu),且Langmuir 體積VL值為25.97~35.59 cm3/g,平均31.93 cm3/g;Langmuir壓力pL值變化介于2.00~3.21 MPa,平均2.59 MPa(表2),且隨著煤體破壞程度增加而減小,說(shuō)明構(gòu)造煤達(dá)到最大吸附量所需能量較低,吸附性增強(qiáng)。煤中甲烷解吸是吸附的逆過(guò)程,原理上吸附曲線和解吸曲線應(yīng)該重合。但是圖2 中顯示,4 種煤體結(jié)構(gòu)對(duì)于甲烷的解吸均存在明顯的滯后現(xiàn)象,且隨著煤破碎程度增加滯后環(huán)有減小的趨勢(shì)。這是由于大孔中毛細(xì)凝聚作用的發(fā)生使得孔隙中凝聚的甲烷凹液面的飽和蒸氣壓力低于水平液面,因此在相同的吸附量,解吸壓力相對(duì)滯后,也說(shuō)明煤層氣開發(fā)過(guò)程中產(chǎn)能主要集中于開發(fā)的中后期。

    圖2 25℃時(shí)不同煤體結(jié)構(gòu)煤中甲烷等溫吸附/解吸曲線Fig.2 CH4 adsorption/desorption curves of different coal structures

    2.2 構(gòu)造煤與原生結(jié)構(gòu)煤的擴(kuò)散特征

    2.2.1 原生結(jié)構(gòu)煤

    在煤層氣開發(fā)過(guò)程中,隨著儲(chǔ)層壓力降低,氣體首先開始解吸隨后進(jìn)入擴(kuò)散階段。因此,為了真實(shí)模擬煤儲(chǔ)層壓力下降過(guò)程中甲烷分子在煤孔隙中的擴(kuò)散效應(yīng),本文選取等溫吸附實(shí)驗(yàn)過(guò)程中甲烷解吸動(dòng)力學(xué)數(shù)據(jù)隨時(shí)間的變化關(guān)系分析不同煤體結(jié)構(gòu)擴(kuò)散特征。圖3 為初始解吸壓力下,實(shí)驗(yàn)樣品缸中壓力變化,可以間接反映不同煤體結(jié)構(gòu)煤中甲烷吸附/解吸和擴(kuò)散特性。圖3 中壓力的變化隨著煤破碎程度的增強(qiáng)而升高,表明構(gòu)造煤中解吸出來(lái)的甲烷要高于原生結(jié)構(gòu)煤,構(gòu)造煤的吸附量升高,同時(shí)構(gòu)造煤甲烷的擴(kuò)散率也要高于原生結(jié)構(gòu)煤。初始解吸壓力下不同煤體結(jié)構(gòu)甲烷擴(kuò)散速率和單、雙孔模型擬合數(shù)據(jù)如圖4 所示,圖中顯示煤中甲烷擴(kuò)散率隨時(shí)間的增加逐漸增大,并在后期逐漸達(dá)到擴(kuò)散平衡狀態(tài)。整個(gè)擴(kuò)散過(guò)程可以劃分為3 個(gè)階段:快速擴(kuò)散階段(S1),衰減階段(S2)和緩慢擴(kuò)散階段(S3)。原生結(jié)構(gòu)煤中甲烷擴(kuò)散效率受控于前期快速擴(kuò)散階段(S1)和衰減階段(S2),這兩個(gè)階段甲烷擴(kuò)散率約占總的擴(kuò)散率的80%。此外,單孔擴(kuò)散模型對(duì)原生結(jié)構(gòu)煤的衰減階段(S2)和緩慢擴(kuò)散階段(S3)擬合度較低,而雙孔擴(kuò)散模型對(duì)整個(gè)擴(kuò)散過(guò)程擬合度均較高。

    圖3 初始解吸壓力條件下甲烷解吸-擴(kuò)散速率變化Fig.3 Variations of the desorption-diffusion rate of methane under the initial desorption pressure

    圖4 煤體結(jié)構(gòu)單、雙孔擴(kuò)散模型數(shù)據(jù)擬合Fig.4 Data fitting diagram of unipore and bidisperse models

    表2 不同煤體結(jié)構(gòu)煤Langmuir 參數(shù)Table 2 Langmuir parameters of different coal structures

    2.2.2 構(gòu)造煤

    構(gòu)造煤中不同煤體結(jié)構(gòu)的擴(kuò)散過(guò)程與原生結(jié)構(gòu)煤相似,包含3 個(gè)階段(圖4)。但是,構(gòu)造煤的擴(kuò)散率在快速擴(kuò)散階段(S1)和衰減階段(S2)都高于原生結(jié)構(gòu)煤,可以達(dá)到90%左右,煤在破壞后擴(kuò)散效率明顯提高。構(gòu)造煤中雙孔擴(kuò)散模型與實(shí)測(cè)數(shù)據(jù)擬合度也相對(duì)較高,可以準(zhǔn)確描述不同破壞程度煤中甲烷擴(kuò)散的動(dòng)態(tài)變化。

    2.3 煤中甲烷擴(kuò)散性能差異性的受控因素

    2.3.1 孔隙結(jié)構(gòu)和氣體壓力

    根據(jù)氣體分子平均自由程與多孔介質(zhì)孔隙大小關(guān)系,通常采用克努森系數(shù)(Kn)表征氣體分子在孔隙中的運(yùn)移效率。依據(jù)克努森系數(shù)的大小,氣體在多孔介質(zhì)中的擴(kuò)散模式可分為[15]:菲克型擴(kuò)散(Kn<0.1)、過(guò)渡型擴(kuò)散(0.1≤Kn≤10) 和克努森型擴(kuò)散(Kn>10)??伺拖禂?shù)定義為:

    式中:d為分子運(yùn)移特征長(zhǎng)度,一般為孔隙直徑,m。λ為分子平均自由程,計(jì)算表達(dá)式為:

    式中:KB為玻爾茲曼常數(shù),1.38×10-23J/K;T為氣體溫度,K;d0為分子有效直徑,nm;p為氣體壓力,MPa。

    圖5 是單孔模型計(jì)算得到的不同煤體結(jié)構(gòu)煤在不同解吸平衡壓力下的有效擴(kuò)散系數(shù)。不同煤體結(jié)構(gòu)煤的有效擴(kuò)散系數(shù)均隨壓力的降低而減小,表明隨著壓力降低甲烷在運(yùn)移過(guò)程中所受阻力逐漸增大。這是由于在不同壓力階段多孔介質(zhì)中甲烷分子平均自由程與孔徑的分布關(guān)系不同造成的。

    圖5 單孔擴(kuò)散模型不同煤體結(jié)構(gòu)煤的有效擴(kuò)散系數(shù)變化規(guī)律Fig.5 Variations of effective diffusion coefficients of different coal structures

    根據(jù)分子平均自由程與孔徑的比值,通常采用Knudsen 系數(shù)表征3 種不同擴(kuò)散阻力的擴(kuò)散模型。隨著煤儲(chǔ)層壓力的降低,甲烷的分子平均自由程增大,擴(kuò)散模型逐漸由菲克型擴(kuò)散、過(guò)渡型擴(kuò)散轉(zhuǎn)變?yōu)榭伺蛿U(kuò)散。初始解吸平衡壓力在6 MPa 左右時(shí),4 種煤體結(jié)構(gòu)煤的擴(kuò)散系數(shù)均最高,其中糜棱結(jié)構(gòu)煤的有效擴(kuò)散系數(shù)最大,為8.76×10-5s-1,其次是碎裂結(jié)構(gòu)煤、碎粒結(jié)構(gòu)煤和原生結(jié)構(gòu)煤,分別為8.73×10-5、7.84×10-5和5.5×10-5s-1。當(dāng)壓力約大于6 MPa 時(shí),分子自由程較小,氣體分子與孔壁碰撞的次數(shù)降低,菲克型擴(kuò)散控制了煤孔隙的整個(gè)擴(kuò)散過(guò)程。隨著壓力的進(jìn)一步降低,分子自由程升高,主要擴(kuò)散模式逐漸向過(guò)渡型擴(kuò)散轉(zhuǎn)變。在1.5 MPa 壓力以下,不同煤體結(jié)構(gòu)的有效擴(kuò)散系數(shù)隨壓力變化存在一定程度波動(dòng)性。研究發(fā)現(xiàn)甲烷在1~6 MPa 的平均自由程在1~5 nm 變化[18-19],對(duì)應(yīng)于克努森系數(shù)計(jì)算的孔徑則是10~50 nm,說(shuō)明在此壓力區(qū)間中孔是過(guò)渡型擴(kuò)散的主要場(chǎng)所。

    隨著氣體壓力和濃度的降低,氣體分子與孔壁分子碰撞逐漸加劇,導(dǎo)致有效擴(kuò)散系數(shù)在1.5~6.0 MPa之間逐漸降低。壓力降低至1.5 MPa 后,氣體分子平均自由程>5 nm,解吸擴(kuò)散主要發(fā)生在較小的孔隙中(<10 nm),擴(kuò)散類型主要為克努森型擴(kuò)散,氣體分子與孔壁的碰撞更加頻繁,氣體運(yùn)移阻力達(dá)到最大。通??伺蛿U(kuò)散階段的擴(kuò)散系數(shù)較低,但圖5 顯示,4個(gè)樣品的有效擴(kuò)散系數(shù)在低壓下均顯示了相似的忽高忽低的波動(dòng)性,例如原生結(jié)構(gòu)煤增高至3.86×10-5s-1,碎裂結(jié)構(gòu)煤6.52×10-5s-1,碎粒結(jié)構(gòu)煤7.52×10-5s-1,糜棱結(jié)構(gòu)煤7.74×10-5s-1。這可能與煤中4 nm 左右的細(xì)瓶頸孔的發(fā)育有關(guān)[19],在平衡壓力下,細(xì)瓶頸孔隙中充滿的氣體處于過(guò)飽和狀態(tài),一旦孔口壓力降低,細(xì)瓶頸孔內(nèi)大量氣體瞬間解吸沖出孔口并造成連鎖效應(yīng),導(dǎo)致孔隙周邊氣體濃度迅速增高,流速加快,擴(kuò)散系數(shù)突然增大。

    構(gòu)造煤的甲烷擴(kuò)散系數(shù)在1.5~6.0 MPa 間高于原生結(jié)構(gòu)煤(圖5),平均甲烷擴(kuò)散系數(shù)是原生結(jié)構(gòu)煤的1.53 倍。在初始解吸壓力約6 MPa 時(shí),擴(kuò)散系數(shù)最高的是碎裂結(jié)構(gòu)煤。隨著壓力的降低,碎裂結(jié)構(gòu)煤的擴(kuò)散系數(shù)急劇下降,在1.5~5.0 MPa 范圍內(nèi),各煤體結(jié)構(gòu)間的擴(kuò)散系數(shù)大小順序?yàn)椋好永饨Y(jié)構(gòu)煤>碎粒結(jié)構(gòu)煤>碎裂結(jié)構(gòu)煤>原生結(jié)構(gòu)煤。初始最大擴(kuò)散系數(shù)表明,碎裂結(jié)構(gòu)煤中直徑大于50 nm 孔隙相對(duì)較多,有利于菲克型擴(kuò)散的發(fā)生。但初始?jí)毫?,碎裂結(jié)構(gòu)煤的擴(kuò)散速率降低,說(shuō)明碎裂結(jié)構(gòu)煤中孔的分布低于碎粒結(jié)構(gòu)煤和糜棱結(jié)構(gòu)煤,當(dāng)氣體壓力低于1.5 MPa 時(shí),隨著破壞強(qiáng)度的增加,構(gòu)造煤中形成更多的超微孔和微孔,相應(yīng)的比表面積明顯增大(圖6),氣體吸附性增強(qiáng),導(dǎo)致甲烷難以從孔隙表面解吸。在低壓階段,此時(shí)各煤樣的氣體擴(kuò)散系數(shù)波動(dòng)異常。這也說(shuō)明低壓下,構(gòu)造煤具有突然增高的解吸率,易發(fā)生瓦斯突出。構(gòu)造煤克努森階段較高的擴(kuò)散系數(shù)說(shuō)明較大的孔隙或封閉的孔隙易被破壞成較小的孔隙(微孔或超微孔),并具有良好的連通性,而且構(gòu)造煤孔隙中氣體運(yùn)移路徑的長(zhǎng)度比原生結(jié)構(gòu)煤短,因此擴(kuò)散系數(shù)相對(duì)較高。

    圖6 不同煤體結(jié)構(gòu)煤的比表面積隨孔徑變化規(guī)律Fig.6 Variations of specific surface area with pore size of different coal structures

    圖7 為雙孔擴(kuò)散模型大孔和微孔中擴(kuò)散系數(shù)隨壓力變化,整體上變化規(guī)律和單孔模型相似,有效擴(kuò)散系數(shù)均隨壓力降低而變小,且微孔中擴(kuò)散系數(shù)明顯低于大孔。

    圖7 雙孔擴(kuò)散模型不同煤體結(jié)構(gòu)煤的有效擴(kuò)散系數(shù)變化規(guī)律Fig.7 Variations of effective diffusion coefficients of coal structure by the bidisperse model

    2.3.2 溫 度

    擴(kuò)散系數(shù)是表征擴(kuò)散能力的量化參數(shù),選取相同壓力下不同溫度解吸數(shù)據(jù),根據(jù)單孔理論擴(kuò)散模型計(jì)算出不同煤體結(jié)構(gòu)在解吸壓力約為6、3 和1 MPa 下的有效擴(kuò)散系數(shù)(圖8)。

    圖8 顯示不同煤體結(jié)構(gòu)在6 MPa 和3 MPa 條件下,有效擴(kuò)散系數(shù)與溫度表現(xiàn)出明顯的正相關(guān)性,溫度越高有效擴(kuò)散系數(shù)越大,45℃條件下明顯高于其他2 個(gè)溫度,說(shuō)明高溫環(huán)境下溫度對(duì)甲烷擴(kuò)散的影響要高于壓力。但是到了低壓區(qū)(1 MPa),有效擴(kuò)散系數(shù)與溫度關(guān)系出現(xiàn)波動(dòng)性,無(wú)單調(diào)變化規(guī)律。在較低壓力下甲烷的擴(kuò)散主要發(fā)生在微孔,其不規(guī)則變化規(guī)律影響了甲烷的擴(kuò)散速率,壓力和孔隙結(jié)構(gòu)對(duì)擴(kuò)散影響超過(guò)了溫度。

    圖8 不同溫度下不同煤體結(jié)構(gòu)煤有效擴(kuò)散系數(shù)Fig.8 Effective diffusion coefficients of different coal structures at different temperatures

    3 結(jié) 論

    a.隨著煤破壞程度增強(qiáng),煤的吸能力增加,Langmuir體積VL值逐漸增大,表現(xiàn)為:原生結(jié)構(gòu)煤<碎裂結(jié)構(gòu)煤<碎粒結(jié)構(gòu)煤<糜棱結(jié)構(gòu)煤,且構(gòu)造煤的解吸效率要高于原生結(jié)構(gòu)煤。

    b.采用了單孔、雙孔擴(kuò)散模型分析不同煤體結(jié)構(gòu)煤中甲烷擴(kuò)散性能及其差異性,發(fā)現(xiàn)雙孔擴(kuò)散模型對(duì)擴(kuò)散數(shù)據(jù)擬合度較高,更適合于模擬構(gòu)造煤中氣體擴(kuò)散效應(yīng);揭示了煤中氣體擴(kuò)散呈3 個(gè)階段變化規(guī)律,即:快速擴(kuò)散階段(S1),衰減階段(S2) 和緩慢擴(kuò)散階段(S3)。

    c.不同煤體結(jié)構(gòu)煤中甲烷的擴(kuò)散系數(shù)均隨壓力的降低而減小,在整個(gè)擴(kuò)散階段構(gòu)造煤的擴(kuò)散系數(shù)高于原生結(jié)構(gòu)煤。在高于5.0 MPa、1.5~5.0 MPa 和低于1.5 MPa 的壓力區(qū)間,氣體擴(kuò)散主要發(fā)生在大孔、中孔和微孔中,分別對(duì)應(yīng)菲克型擴(kuò)散、過(guò)渡型擴(kuò)散和克努森型擴(kuò)散。

    d.構(gòu)造煤與原生結(jié)構(gòu)煤相比,孔徑減小,比表面積增大,孔隙結(jié)構(gòu)復(fù)雜度增強(qiáng);有效擴(kuò)散系數(shù)與溫度表現(xiàn)出明顯的正相關(guān)性,溫度越高有效擴(kuò)散系數(shù)越大,但是到了低壓階段孔隙結(jié)構(gòu)對(duì)煤中甲烷擴(kuò)散性影響高于溫度。

    猜你喜歡
    擴(kuò)散系數(shù)煤體微孔
    注熱井周圍煤體蠕變過(guò)程的滲透率變化規(guī)律模擬研究
    正交試驗(yàn)下煤體滲透性影響因素評(píng)價(jià)
    煤炭工程(2019年4期)2019-05-05 02:30:16
    以“中央廚房”為突破口探索時(shí)政報(bào)道的融煤體之路——以浙江之聲為例
    強(qiáng)疏水性PDMS/PVDF微孔膜的制備及其性能研究
    基于Sauer-Freise 方法的Co- Mn 體系fcc 相互擴(kuò)散系數(shù)的研究
    上海金屬(2015年5期)2015-11-29 01:13:59
    FCC Ni-Cu 及Ni-Mn 合金互擴(kuò)散系數(shù)測(cè)定
    上海金屬(2015年6期)2015-11-29 01:09:09
    膜蒸餾用PDMS/PVDF/PTFE三元共混微孔膜制備
    微孔發(fā)泡塑料中成核劑的研究
    非時(shí)齊擴(kuò)散模型中擴(kuò)散系數(shù)的局部估計(jì)
    深部區(qū)域塊段煤體瓦斯基本參數(shù)測(cè)定技術(shù)研究
    河南科技(2014年3期)2014-02-27 14:06:05
    99久久99久久久精品蜜桃| 国产欧美日韩一区二区精品| 天堂√8在线中文| 久久精品人妻少妇| 国产精品98久久久久久宅男小说| 成人无遮挡网站| 99在线人妻在线中文字幕| 亚洲专区中文字幕在线| 国产高清三级在线| 狂野欧美白嫩少妇大欣赏| 国产三级黄色录像| av国产免费在线观看| 免费在线观看影片大全网站| 国产午夜精品论理片| www.熟女人妻精品国产| 久久久久九九精品影院| 90打野战视频偷拍视频| 久久亚洲精品不卡| 日本黄色视频三级网站网址| 国产免费男女视频| 精品久久久久久久久久免费视频| 97超级碰碰碰精品色视频在线观看| 欧美激情在线99| 免费高清视频大片| 亚洲精品一区av在线观看| 欧美区成人在线视频| 国产亚洲精品综合一区在线观看| 丰满人妻熟妇乱又伦精品不卡| 日韩亚洲欧美综合| 久久久久久国产a免费观看| 一区福利在线观看| 国产精品综合久久久久久久免费| 高潮久久久久久久久久久不卡| 亚洲精华国产精华精| 亚洲狠狠婷婷综合久久图片| 亚洲激情在线av| 国产乱人视频| 成人高潮视频无遮挡免费网站| 99热这里只有是精品50| 首页视频小说图片口味搜索| 亚洲七黄色美女视频| 欧美黑人欧美精品刺激| 2021天堂中文幕一二区在线观| av天堂中文字幕网| 国产v大片淫在线免费观看| 18美女黄网站色大片免费观看| 国产精品日韩av在线免费观看| 国产主播在线观看一区二区| 久久欧美精品欧美久久欧美| 熟女电影av网| 波多野结衣巨乳人妻| 亚洲av电影不卡..在线观看| 午夜精品在线福利| 亚洲欧美日韩卡通动漫| 黄色女人牲交| 村上凉子中文字幕在线| 亚洲中文字幕日韩| 欧美一级a爱片免费观看看| 国产真实乱freesex| 亚洲国产精品成人综合色| 一个人免费在线观看的高清视频| 天天添夜夜摸| 男女床上黄色一级片免费看| 欧美又色又爽又黄视频| 亚洲专区国产一区二区| 成年女人毛片免费观看观看9| 真人一进一出gif抽搐免费| 嫩草影院入口| 超碰av人人做人人爽久久 | 亚洲人成网站在线播放欧美日韩| 久久久久免费精品人妻一区二区| 亚洲久久久久久中文字幕| 黄片大片在线免费观看| 欧美成人性av电影在线观看| 久久久久久大精品| 草草在线视频免费看| 亚洲欧美日韩无卡精品| 99国产极品粉嫩在线观看| 999久久久精品免费观看国产| 欧美激情久久久久久爽电影| 久久精品国产综合久久久| 精品国内亚洲2022精品成人| 亚洲avbb在线观看| 精品无人区乱码1区二区| 国模一区二区三区四区视频| а√天堂www在线а√下载| 国产精品爽爽va在线观看网站| 岛国视频午夜一区免费看| 国产真人三级小视频在线观看| 精品人妻1区二区| 亚洲 欧美 日韩 在线 免费| 少妇高潮的动态图| 好看av亚洲va欧美ⅴa在| 午夜视频国产福利| 久久人妻av系列| 婷婷丁香在线五月| 国产中年淑女户外野战色| 午夜精品在线福利| 性色avwww在线观看| 久久久久国内视频| 日本a在线网址| 少妇的逼好多水| 宅男免费午夜| 久久99热这里只有精品18| 美女免费视频网站| 嫩草影视91久久| 亚洲成人久久性| 动漫黄色视频在线观看| 国产精品1区2区在线观看.| 内射极品少妇av片p| 欧美日韩国产亚洲二区| 一区二区三区国产精品乱码| 国产欧美日韩精品一区二区| 中文字幕久久专区| 欧美一级a爱片免费观看看| 制服人妻中文乱码| 在线播放国产精品三级| 亚洲第一欧美日韩一区二区三区| 麻豆成人午夜福利视频| 国产欧美日韩一区二区精品| 俺也久久电影网| 一边摸一边抽搐一进一小说| 天天躁日日操中文字幕| 一个人免费在线观看电影| 偷拍熟女少妇极品色| 九九热线精品视视频播放| 免费高清视频大片| 国产国拍精品亚洲av在线观看 | 精品久久久久久,| 日韩精品青青久久久久久| 日韩免费av在线播放| 熟女电影av网| 国产精品亚洲一级av第二区| 婷婷亚洲欧美| 国产精品国产高清国产av| 欧美精品啪啪一区二区三区| 免费看光身美女| 国产免费男女视频| 91九色精品人成在线观看| 99久久99久久久精品蜜桃| 国产激情偷乱视频一区二区| 美女被艹到高潮喷水动态| 国产探花在线观看一区二区| 最新美女视频免费是黄的| 尤物成人国产欧美一区二区三区| 欧美中文综合在线视频| 精品日产1卡2卡| 亚洲av电影不卡..在线观看| 国产精品 欧美亚洲| 亚洲欧美精品综合久久99| 黄片大片在线免费观看| www.www免费av| 性色avwww在线观看| 亚洲无线在线观看| 久久久久久国产a免费观看| 五月玫瑰六月丁香| 国产欧美日韩精品亚洲av| 琪琪午夜伦伦电影理论片6080| 有码 亚洲区| 波多野结衣高清作品| 哪里可以看免费的av片| 国产一区二区在线av高清观看| 国产精品爽爽va在线观看网站| 国产精品日韩av在线免费观看| 国产一区二区亚洲精品在线观看| 亚洲精品国产精品久久久不卡| 久久久久久国产a免费观看| 综合色av麻豆| 欧美黑人欧美精品刺激| 五月玫瑰六月丁香| 亚洲精品色激情综合| 美女高潮喷水抽搐中文字幕| 亚洲av美国av| 国产国拍精品亚洲av在线观看 | 色综合站精品国产| 19禁男女啪啪无遮挡网站| 久久久久久国产a免费观看| 最近视频中文字幕2019在线8| 九九久久精品国产亚洲av麻豆| 欧美午夜高清在线| 亚洲精品色激情综合| 国产三级黄色录像| 最近最新中文字幕大全免费视频| 超碰av人人做人人爽久久 | 伊人久久精品亚洲午夜| 亚洲va日本ⅴa欧美va伊人久久| 国内精品一区二区在线观看| 国产精品99久久99久久久不卡| 99视频精品全部免费 在线| 一级黄色大片毛片| 又黄又爽又免费观看的视频| 国产精品乱码一区二三区的特点| 国产亚洲精品一区二区www| 国产成人av激情在线播放| 亚洲天堂国产精品一区在线| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲五月婷婷丁香| 日韩欧美免费精品| 最近视频中文字幕2019在线8| 嫩草影院精品99| 日韩 欧美 亚洲 中文字幕| 国产不卡一卡二| 一级a爱片免费观看的视频| 男人舔女人下体高潮全视频| 啦啦啦免费观看视频1| 亚洲精品成人久久久久久| 在线视频色国产色| av欧美777| 久久国产乱子伦精品免费另类| 女人高潮潮喷娇喘18禁视频| 久久精品人妻少妇| 精品人妻偷拍中文字幕| 午夜福利高清视频| 日韩人妻高清精品专区| 国产高清视频在线播放一区| 悠悠久久av| 动漫黄色视频在线观看| 国产免费一级a男人的天堂| 757午夜福利合集在线观看| 日韩国内少妇激情av| 欧美绝顶高潮抽搐喷水| 欧美不卡视频在线免费观看| 欧美激情久久久久久爽电影| 成年版毛片免费区| h日本视频在线播放| 亚洲电影在线观看av| 一个人看的www免费观看视频| 久久国产精品影院| 国产精品久久久久久久电影 | 日本 av在线| 久久亚洲精品不卡| 十八禁网站免费在线| 乱人视频在线观看| 欧美黑人巨大hd| 欧美色视频一区免费| 好男人在线观看高清免费视频| 国内精品久久久久久久电影| 成人性生交大片免费视频hd| 婷婷丁香在线五月| 亚洲成av人片在线播放无| 99久久久亚洲精品蜜臀av| 热99在线观看视频| 日日干狠狠操夜夜爽| x7x7x7水蜜桃| 免费观看精品视频网站| 亚洲专区国产一区二区| 亚洲欧美日韩无卡精品| 久久久国产精品麻豆| 成人午夜高清在线视频| 别揉我奶头~嗯~啊~动态视频| 国产一级毛片七仙女欲春2| 免费无遮挡裸体视频| 91久久精品电影网| www日本黄色视频网| 亚洲18禁久久av| 极品教师在线免费播放| 99久久无色码亚洲精品果冻| 毛片女人毛片| 美女 人体艺术 gogo| 亚洲人成网站在线播| 99久久精品一区二区三区| 亚洲欧美日韩东京热| 1000部很黄的大片| 午夜福利18| 婷婷六月久久综合丁香| 淫秽高清视频在线观看| 亚洲国产精品成人综合色| 国产亚洲av嫩草精品影院| 熟女人妻精品中文字幕| 亚洲av电影不卡..在线观看| av黄色大香蕉| 国产午夜精品论理片| 日本黄色视频三级网站网址| 国产精品影院久久| 亚洲精品成人久久久久久| 一区福利在线观看| 19禁男女啪啪无遮挡网站| 啪啪无遮挡十八禁网站| 51国产日韩欧美| 免费人成在线观看视频色| 国产免费一级a男人的天堂| 3wmmmm亚洲av在线观看| 亚洲美女视频黄频| 黄色视频,在线免费观看| av视频在线观看入口| 国产亚洲欧美98| 精品一区二区三区av网在线观看| 国产免费男女视频| 国产一区二区三区视频了| 午夜福利欧美成人| 露出奶头的视频| 午夜福利在线观看吧| 99riav亚洲国产免费| 热99re8久久精品国产| 熟女少妇亚洲综合色aaa.| 中文字幕高清在线视频| 国产精品久久久久久久电影 | 国产高清有码在线观看视频| 国产淫片久久久久久久久 | 熟妇人妻久久中文字幕3abv| 国产精品99久久99久久久不卡| ponron亚洲| 男女做爰动态图高潮gif福利片| 久久精品综合一区二区三区| 嫩草影院入口| 非洲黑人性xxxx精品又粗又长| 国产成人av激情在线播放| 欧美日韩瑟瑟在线播放| 国产三级中文精品| 欧美午夜高清在线| 欧美日韩乱码在线| 久99久视频精品免费| 一卡2卡三卡四卡精品乱码亚洲| 免费观看精品视频网站| 久久国产精品影院| 国产午夜精品久久久久久一区二区三区 | ponron亚洲| 一进一出好大好爽视频| 嫁个100分男人电影在线观看| 国产成人a区在线观看| 少妇高潮的动态图| 无遮挡黄片免费观看| 亚洲精品粉嫩美女一区| 母亲3免费完整高清在线观看| 婷婷丁香在线五月| 免费一级毛片在线播放高清视频| 无限看片的www在线观看| 亚洲 国产 在线| 日韩欧美一区二区三区在线观看| 一个人免费在线观看电影| 久久久国产成人免费| 制服人妻中文乱码| 欧美精品啪啪一区二区三区| 欧美+亚洲+日韩+国产| 波多野结衣高清作品| 中亚洲国语对白在线视频| 18禁裸乳无遮挡免费网站照片| 欧美一级毛片孕妇| 亚洲国产欧美网| 国产99白浆流出| 午夜a级毛片| 听说在线观看完整版免费高清| 日韩亚洲欧美综合| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲真实伦在线观看| 色吧在线观看| 嫩草影院入口| 久久精品国产99精品国产亚洲性色| 成人欧美大片| 在线a可以看的网站| av国产免费在线观看| 国产精品精品国产色婷婷| 亚洲欧美日韩高清专用| 日韩欧美国产在线观看| 激情在线观看视频在线高清| 99国产极品粉嫩在线观看| 亚洲熟妇熟女久久| 免费搜索国产男女视频| 免费看日本二区| 中文字幕久久专区| 在线天堂最新版资源| 色噜噜av男人的天堂激情| 欧美日韩乱码在线| 51国产日韩欧美| 欧美最黄视频在线播放免费| 国产亚洲精品一区二区www| 亚洲精品亚洲一区二区| 亚洲av免费在线观看| 午夜老司机福利剧场| 精品一区二区三区av网在线观看| 久久久久久国产a免费观看| 亚洲精品色激情综合| 在线观看日韩欧美| 9191精品国产免费久久| 久久精品国产自在天天线| 亚洲精品影视一区二区三区av| 国产高清三级在线| 神马国产精品三级电影在线观看| 欧美日韩乱码在线| 精品久久久久久久久久久久久| 国产午夜福利久久久久久| 琪琪午夜伦伦电影理论片6080| 99久国产av精品| 欧美区成人在线视频| 长腿黑丝高跟| 色视频www国产| 激情在线观看视频在线高清| 亚洲18禁久久av| 成年免费大片在线观看| 老司机在亚洲福利影院| 88av欧美| 国产日本99.免费观看| 在线十欧美十亚洲十日本专区| 久久99热这里只有精品18| 亚洲人成电影免费在线| 日韩国内少妇激情av| 免费电影在线观看免费观看| 亚洲欧美日韩高清专用| 亚洲午夜理论影院| 精华霜和精华液先用哪个| 十八禁人妻一区二区| 丁香六月欧美| 搡老岳熟女国产| 国产亚洲av嫩草精品影院| 午夜日韩欧美国产| 国产不卡一卡二| 真人一进一出gif抽搐免费| 欧美日本视频| www.色视频.com| 欧美日韩一级在线毛片| 国产伦在线观看视频一区| www.色视频.com| 51国产日韩欧美| h日本视频在线播放| 19禁男女啪啪无遮挡网站| 老司机福利观看| 国产99白浆流出| 久久久久久久精品吃奶| 黑人欧美特级aaaaaa片| 无遮挡黄片免费观看| 美女cb高潮喷水在线观看| 少妇的逼水好多| 国产主播在线观看一区二区| 国产美女午夜福利| 国产毛片a区久久久久| 国产国拍精品亚洲av在线观看 | 99久久成人亚洲精品观看| 97碰自拍视频| 亚洲精华国产精华精| 小蜜桃在线观看免费完整版高清| 国产一区二区在线观看日韩 | 黄色丝袜av网址大全| 嫩草影视91久久| 中亚洲国语对白在线视频| 少妇的逼水好多| 色综合亚洲欧美另类图片| 黄色片一级片一级黄色片| 亚洲最大成人中文| 亚洲18禁久久av| 色综合欧美亚洲国产小说| 精品国产亚洲在线| 亚洲av一区综合| 两人在一起打扑克的视频| 国产精品亚洲一级av第二区| 嫩草影院精品99| 最新美女视频免费是黄的| 欧美日本亚洲视频在线播放| 久久天躁狠狠躁夜夜2o2o| 少妇的逼水好多| 搡女人真爽免费视频火全软件 | 成人特级av手机在线观看| 一卡2卡三卡四卡精品乱码亚洲| 999久久久精品免费观看国产| 日本a在线网址| 精品国内亚洲2022精品成人| 亚洲精品成人久久久久久| 日本黄色片子视频| 亚洲精品一卡2卡三卡4卡5卡| 90打野战视频偷拍视频| 精品久久久久久久久久久久久| 99国产极品粉嫩在线观看| 禁无遮挡网站| 成人高潮视频无遮挡免费网站| 九色成人免费人妻av| 一卡2卡三卡四卡精品乱码亚洲| ponron亚洲| 天天躁日日操中文字幕| 色综合婷婷激情| 日韩欧美三级三区| 亚洲欧美精品综合久久99| 精品电影一区二区在线| 亚洲国产精品成人综合色| 人人妻人人看人人澡| 久久久久久久久大av| 成年女人毛片免费观看观看9| 人妻夜夜爽99麻豆av| 欧美成人a在线观看| 18+在线观看网站| 亚洲av美国av| 精品99又大又爽又粗少妇毛片 | bbb黄色大片| 日本成人三级电影网站| 观看美女的网站| 精品国内亚洲2022精品成人| 日本与韩国留学比较| 久9热在线精品视频| 国语自产精品视频在线第100页| 欧美性感艳星| 怎么达到女性高潮| 九九久久精品国产亚洲av麻豆| 桃色一区二区三区在线观看| 一级黄色大片毛片| 国内毛片毛片毛片毛片毛片| 亚洲在线观看片| 国产伦在线观看视频一区| 欧美黑人欧美精品刺激| 亚洲精品乱码久久久v下载方式 | 午夜精品久久久久久毛片777| 免费无遮挡裸体视频| 黄色视频,在线免费观看| 99热精品在线国产| 黄色女人牲交| 国产单亲对白刺激| 免费在线观看影片大全网站| 国模一区二区三区四区视频| 国产亚洲精品一区二区www| 有码 亚洲区| 亚洲在线自拍视频| 久久婷婷人人爽人人干人人爱| 国产亚洲欧美98| 好男人电影高清在线观看| 久久午夜亚洲精品久久| 国产精品一区二区三区四区免费观看 | 国产老妇女一区| 在线视频色国产色| 中文字幕精品亚洲无线码一区| 国产精品综合久久久久久久免费| 国产精品久久视频播放| 欧美+日韩+精品| 舔av片在线| 国产精品电影一区二区三区| 国产91精品成人一区二区三区| 欧美极品一区二区三区四区| 99国产极品粉嫩在线观看| 欧美日韩一级在线毛片| 国产精品久久久久久久电影 | 国产亚洲精品av在线| 亚洲精品粉嫩美女一区| 亚洲在线自拍视频| 一本精品99久久精品77| 国产高清激情床上av| xxxwww97欧美| 好男人电影高清在线观看| 精品乱码久久久久久99久播| svipshipincom国产片| 村上凉子中文字幕在线| 在线免费观看不下载黄p国产 | 亚洲最大成人手机在线| 在线观看66精品国产| 亚洲七黄色美女视频| 久久精品国产亚洲av涩爱 | 人人妻人人看人人澡| 亚洲精品粉嫩美女一区| 国产色婷婷99| 99热只有精品国产| 国产高清videossex| 亚洲美女视频黄频| 欧美中文综合在线视频| 丁香欧美五月| 美女高潮的动态| 婷婷精品国产亚洲av在线| 一级黄片播放器| 欧美又色又爽又黄视频| 国产中年淑女户外野战色| 又粗又爽又猛毛片免费看| 国产美女午夜福利| 伊人久久大香线蕉亚洲五| 一进一出抽搐动态| 人妻夜夜爽99麻豆av| 欧美成狂野欧美在线观看| 亚洲国产高清在线一区二区三| 国产精品久久久久久人妻精品电影| 国产精品99久久久久久久久| 99国产综合亚洲精品| 久久久久久久久中文| 九九热线精品视视频播放| 一区二区三区国产精品乱码| 国产三级黄色录像| 成人欧美大片| 国产精品一及| 99久久成人亚洲精品观看| 亚洲激情在线av| 观看免费一级毛片| 免费电影在线观看免费观看| 级片在线观看| 国产精品久久久人人做人人爽| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美一区二区三区在线观看| 亚洲精品456在线播放app | 噜噜噜噜噜久久久久久91| 看片在线看免费视频| 国产一区二区在线av高清观看| 成人一区二区视频在线观看| 一二三四社区在线视频社区8| 国产精品一区二区免费欧美| 日本精品一区二区三区蜜桃| 中文字幕久久专区| 亚洲av成人av| 夜夜躁狠狠躁天天躁| 亚洲18禁久久av| 亚洲av成人精品一区久久| 久久精品影院6| 真人做人爱边吃奶动态| 国产精品爽爽va在线观看网站| 国产激情偷乱视频一区二区| 国产精品,欧美在线| 观看免费一级毛片| 中文字幕人妻丝袜一区二区| 国产色婷婷99| 亚洲无线在线观看| 午夜日韩欧美国产| 我要搜黄色片| 欧美日韩精品网址| 天天躁日日操中文字幕| 国产亚洲精品久久久久久毛片| 给我免费播放毛片高清在线观看| 亚洲成人免费电影在线观看| 亚洲五月婷婷丁香| 欧美三级亚洲精品| 久久久国产成人免费| 男人舔奶头视频| 亚洲成人久久爱视频| 国产精品综合久久久久久久免费| 亚洲国产色片| 国产精品久久久久久亚洲av鲁大|