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

    頁巖中甲烷等溫吸附量計算問題及方法改進

    2018-03-01 05:58:14舒向偉朱忠云方子和
    石油實驗地質(zhì) 2018年1期
    關(guān)鍵詞:等溫甲烷頁巖

    方 帆,孫 沖,舒向偉,朱忠云,方子和

    (1.中國石化 江漢油田分公司 勘探開發(fā)研究院,武漢 430223; 2.長江大學(xué) 地球科學(xué)學(xué)院,武漢 430100)

    頁巖氣是賦存于頁巖地層中的一種非常規(guī)天然氣資源,頁巖的含氣性評價是一項重要的基礎(chǔ)性工作,準(zhǔn)確、客觀地分析和評價頁巖地層的含氣性,對于頁巖氣資源潛力分析以及勘探開發(fā)部署都具有重要指導(dǎo)意義[1]。在頁巖吸附性能的研究中,等溫吸附實驗是最常用的手段。對于頁巖氣來說,主要以吸附態(tài)和游離態(tài)2種賦存形式存在:游離態(tài)氣體賦存于納米級孔隙和微裂縫中,吸附態(tài)氣體賦存于有機質(zhì)及黏土礦物孔隙中[2-3]。

    目前頁巖等溫吸附測量主要有2種方法:容量法和重量法。容量法是根據(jù)物質(zhì)守衡原理,測量吸附前后甲烷氣體摩爾數(shù)的變化來計算吸附量;重量法根據(jù)力學(xué)平衡原理,測量吸附前后樣品重量的變化來計算吸附量。近年來隨著非常規(guī)天然氣的蓬勃發(fā)展,實驗研究人員通過對頁巖和煤層的等溫吸附研究發(fā)現(xiàn),與煤層不一樣的是,大量頁巖樣品等溫吸附曲線高壓段常常出現(xiàn)下拐現(xiàn)象,這偏離了傳統(tǒng)的單分子層吸附理論,無論是容量法還是重量法都存在此問題[4-7]。基于容量法,筆者對頁巖等溫吸附實驗及計算方面存在的問題加以分析,并進一步探討改進目前的實驗方法。

    1 目前問題

    1.1 實驗曲線存在下拐

    在化石能源領(lǐng)域,等溫吸附實驗較早應(yīng)用于煤層氣,隨后擴展至頁巖氣領(lǐng)域。目前頁巖等溫吸附實驗方法和設(shè)備都是借鑒套用煤層氣實驗方法和設(shè)備。實驗室測定方法主要是借鑒《GB/T 19560—2008煤的高壓等溫吸附試驗方法》[8],實驗設(shè)備較有代表性的是美國Terra Tek公司的ISO-300型等溫吸附儀。由于頁巖儲層通常比較深,儲層溫度高、壓力高,同時頁巖的有機質(zhì)含量低,吸附能力低,加之頁巖等溫吸附實驗儀器和方法存在差異,各實驗室測得的結(jié)果也不盡相同。無論是容量法還是重量法,頁巖等溫吸附實驗都常會遇到“倒吸”的現(xiàn)象,即中高壓段吸附量不升反降、出現(xiàn)下拐,測得的Langmuir體積和Langmuir壓力甚至出現(xiàn)負值,這與實際不相符,造成實驗測試數(shù)據(jù)失去應(yīng)用價值。

    1.2 過剩吸附與絕對吸附區(qū)分

    頁巖等溫吸附應(yīng)用的基本理論為單分子層吸附理論,由該理論建立Langmuir方程模型。等溫吸附實驗測得的氣體吸附量為過剩吸附量(也稱Gibbs吸附量、相對吸附量等),而理論吸附模型中的氣體吸附量為絕對吸附量。根據(jù)Gibbs定義,過剩吸附量和絕對吸附量間的關(guān)系為[9]:

    Vabs=Vex/(1-ρg/ρa)

    (1)

    式中:Vabs為平衡壓力下甲烷的絕對吸附量,cm3/g;Vex為平衡壓力下甲烷的過剩吸附量,cm3/g;ρa為甲烷吸附相的密度,g/cm3;ρg為平衡壓力下甲烷氣相的密度,g/cm3。

    煤巖等溫吸附曲線較少出現(xiàn)下拐情況,原因一是煤巖吸附量大,通常達十幾至幾十m3/t,頁巖吸附量相對較少,一般不超過5 m3/t[10-12];二是煤巖埋深較淺,實驗壓力不高,通常不到10 MPa,而頁巖埋深約2 000~4 000 m,地層壓力高,實驗壓力也高,通常高達20 MPa以上。從甲烷PVT相圖來看(圖1),煤層氣通常處于氣相區(qū)和低壓超臨界區(qū),頁巖氣則處于高壓超臨界區(qū)(涪陵地區(qū)頁巖地層溫度40~90 ℃),頁巖吸附曲線正是在中高壓段出現(xiàn)下拐。

    對于煤層氣來說,因常處于氣相區(qū)和低壓超臨界區(qū)(<10 MPa),實驗壓力不高,ρg較小(表1),因此Vabs與Vex相差不大,無需考慮這兩者之間的差別。

    對于頁巖氣來說,隨著實驗壓力的升高,甲烷氣相密度ρg也隨壓力升高而不斷趨近于ρa(通常認為ρa在0.340~0.424 g/cm3之間),此時需要將過剩吸附量轉(zhuǎn)換為絕對吸附量。

    1.3 吸附相密度

    由于頁巖等溫吸附需要區(qū)分過剩吸附量和絕對吸附量,這就需要知道吸附相密度ρa,但至今尚無任何方法直接測定超臨界條件下吸附相密度ρa。理論上認為氣體吸附是由于氣體分子在吸附劑表面放熱失去平均移動能而形成吸附相,失去平均移動能的分子將會由于分子之間的作用力使相互之間的距離到達勢能最小值,類似于氣體液化。因其具有準(zhǔn)液體的性質(zhì),考慮到液體的難以壓縮性,可認為ρa在一定溫度和壓力范圍內(nèi)(如50 MPa以下)是一定值[13]。

    表1 甲烷氣體密度與壓力關(guān)系

    當(dāng)然也有研究人員認為吸附相密度隨不同吸附劑和壓力變化,它不是固定值而是動態(tài)變化的,隨著壓力的增大,甲烷吸附相密度ρa不斷增大,但高壓時變化不大并且趨于定值[14];即便如此,我們?nèi)钥蓪⒓淄槲较嗝芏纫暈槎ㄖ担驗槲覀冄芯靠紤]的是高壓段,低壓段吸附相密度考慮與否對吸附曲線的影響甚小。

    根據(jù)甲烷PVT相圖,甲烷ρa值處于臨界密度(10.13 mol/L)到常壓沸點液體密度(26.56 mol/L)之間,即0.163~0.424 g/cm3,略低于常壓沸點液體甲烷的密度。胡濤等人基于L-F方程,利用活性炭等溫吸附實驗,得到10~50 ℃甲烷吸附相密度在21.718~23.592 mol/L[13],經(jīng)換算在0.347~0.377 g/cm3。Dubinin通過范德華方程以及2個經(jīng)驗公式[15],得到了吸附態(tài)甲烷的密度為0.371 g/cm3。Haydel等在較低的壓力(6.895 MPa)和溫度(40 ℃)下使用無機壁材SiO2進行實驗[16],得到吸附相密度為0.374 g/cm3。Menon提出吸附相的密度等于該物質(zhì)液相的密度。Tsai等將液相密度具體化,認為吸附相密度等于該物質(zhì)臨界點處的密度。對于這種觀點,Mavor通過實驗發(fā)現(xiàn),吸附相密度應(yīng)采用0.101 MPa和-161 ℃條件下液態(tài)甲烷的密度0.422 3 g/cm3[17]。Ambrose結(jié)合分子模擬和Langmuir方程,得出吸附相密度0.340 g/cm3[18]。

    總之,大量研究傾向于吸附相密度可能在0.340~0.424 g/cm3。

    1.4 吸附相體積的影響

    現(xiàn)有頁巖等溫吸附實驗中,計算甲烷吸附量時自由體積一直采用恒定值,而甲烷在頁巖顆粒或有機質(zhì)表面吸附必然會占據(jù)部分空間體積,因此自由體積不是一個恒量。如果不考慮其他因素的影響,隨著吸附量的增加,其自由體積是不斷減小的。若自由體積偏大,則游離態(tài)甲烷氣體量偏大,導(dǎo)致計算吸附量偏小,因此吸附量計算時需要考慮吸附相體積的影響。

    1.5 自由體積測量

    自由空間體積是指樣品缸裝入頁巖樣品后頁巖顆粒之間的空隙、頁巖顆粒內(nèi)部空隙、樣品缸剩余的自由空間、連接管和閥門內(nèi)部空間體積的總和,簡稱自由體積。這個自由空間是相對甲烷來說的,即實際要測量的是甲烷的自由活動空間體積,按理應(yīng)該用甲烷氣來測量,但由于甲烷與巖石接觸存在吸附而不能采用,于是實驗中自由空間體積是采用氦氣來標(biāo)定的。氦氣的分子直徑(0.26 nm)比甲烷的分子直徑(0.38 nm)小[19],隨著壓力增大,氦氣分子能進入更小的納米級喉道孔隙,而甲烷分子不能進入,這會導(dǎo)致計算的自由空間體積偏大。筆者以重慶涪陵JYA井龍馬溪組頁巖樣品為例:樣品量220.80 g,氦氣測得實驗樣品自由體積108.815 cm3,實驗結(jié)果Langmuir體積1.65 cm3/g,Langmuir壓力1.56 MPa。在其所有實驗數(shù)據(jù)不變的條件下,僅將自由體積作微小的變動,研究頁巖吸附量的變化情況(表2,圖2)。

    結(jié)果表明,自由體積量不到1%的變化,可以引起Langmuir體積幾倍或十幾倍的變化。測量的自由體積量若偏小,吸附曲線形態(tài)上移;若偏大,曲線形態(tài)則下移。正因為自由體積測量的微小變化對吸附曲線形態(tài)以及吸附量有極為重大的影響,甲烷自由體積的準(zhǔn)確測量尤為重要。

    表2 自由體積變動對Langmuir體積(壓力)的影響

    圖2 自由體積對等溫吸附曲線的影響

    國外也有研究人員研究了微小孔隙結(jié)構(gòu)對自由體積(孔隙度)測量的影響,以及微小孔隙分布對吸附量的影響,認為用氦氣來測量甲烷的自由空間體積存在瑕疵:氦氣測得的自由體積值偏大,導(dǎo)致頁巖吸附量值偏小。自由體積偏差大小受到巖石微孔隙結(jié)構(gòu)和分布的影響,因此認為頁巖吸附曲線下拐根源在于頁巖微孔隙的影響[20]。

    自由體積的測量還存在一個問題,同一樣品在不同壓力下測得的自由體積量是不同的,自由體積該如何選擇也是問題。張慶玲用氦氣標(biāo)定2塊頁巖樣品的自由空間體積[6],在2 MPa和10 MPa平衡壓力下測得的自由體積相差2 cm3,自由體積相對變化量達6%。不同的自由體積計算得到不同的吸附量。

    1.6 計量精度及缸體體積的穩(wěn)定性

    頁巖等溫吸附實驗對設(shè)備的計量要求非常高。溫控要精準(zhǔn),波動范圍要控制在±0.3 ℃以內(nèi);參考缸和樣品缸壓力傳感器量程要高(不低于5 000 psi),精度不低于0.1%,要成對校準(zhǔn)比對,高壓時兩表絕對偏差不能超過5 psi。參考缸和樣品缸體積的穩(wěn)定性要求高,實驗溫度改變時其體積要重新標(biāo)定,精度不能低于0.1 cm3,并每隔一段時間需要重新標(biāo)定。參考缸和樣品缸都是彈性體,考慮到器壁的壓變性,須測定同一溫度不同壓力下缸體體積變化量,高壓下對于容量變化較大的缸體在吸附量計算中還需加以校正。

    2 方法改進

    2.1 計算時考慮過剩吸附與絕對吸附差異

    基于以上問題,頁巖等溫吸附量計算時需將過剩吸附量轉(zhuǎn)化為絕對吸附量。二元Langmuir方程:

    (2)

    式中:P為氣體平衡壓力,MPa;VL為Langmuir體積,即飽和吸附量,cm3/g;PL為Langmuir壓力,即吸附量為飽和吸附量一半時所對應(yīng)的壓力,MPa。

    Langmuir方程變化為:

    (3)

    這里ρg可采用真實氣體狀態(tài)方程求得,也可采用其他方程,如BWR或SRK等計算[21]。這里ρa的獲取有3種途徑:公式法、給定值和擬合值,采用較多的為后兩者。筆者曾采用實驗數(shù)據(jù)來擬合ρa,但擬合值收斂性較差,ρa分布范圍太廣(0.3~5.0 g/cm3),有的實驗數(shù)據(jù)不收斂,無法擬合。對于給定值法,前人選取的定值有0.423g/cm3(對應(yīng)沸點時液相甲烷密度)[22-23]或者0.373 g/cm3(van der Waals密度)[24]。國內(nèi)有的研究人員采用定值0.375 g/cm3或0.421 g/cm3。

    因等溫吸附實驗所得吸附量為過剩吸附量,其曲線下拐屬于正?,F(xiàn)象。將過剩吸附量轉(zhuǎn)化為絕對吸附量后,低壓段形態(tài)變化不大,高壓段其吸附曲線向下彎曲程度有所減緩,吸附曲線上移??梢钥闯?,ρa分別采用0.375 g/cm3或0.421 g/cm3,其絕對吸附曲線形態(tài)相差不大(圖3a)。

    2.2 考慮吸附相體積的影響

    吸附量計算時,以前一直將自由體積作為恒量來計算,因甲烷吸附在巖石顆粒表面占據(jù)了部分體積空間,即自由體積是變化的,隨壓力升高其自由體積減小。各壓力下吸附相體積計算公式為:

    V=16n/ρa

    (4)

    式中:V為吸附相體積,cm3;n為被頁巖吸附的氣體摩爾數(shù),mol。

    由圖3b實例可以看出:考慮了吸附相體積的影響,計算的吸附量又有所增加,絕對吸附曲線繼續(xù)上移,吸附曲線下拐趨勢又進一步得到減緩。

    2.3 自由體積問題的解決方法

    近年來研究人員從不同方面研究了頁巖曲線下拐的問題,有的從頁巖成分、孔隙結(jié)構(gòu)分析,有的從實驗計算誤差角度分析,有的從吸附模型適應(yīng)性分析。一些觀點認為,頁巖中甲烷等溫吸附并不遵循Langmuir理論模型[4,25],也有的研究認為頁巖中甲烷等溫吸附遵循Langmuir理論模型[14,26]。頁巖中甲烷等溫吸附是否遵循Langmuir理論模型,首先實驗數(shù)據(jù)要正確,否則談模型適應(yīng)性無意義。

    實驗中將過剩吸附轉(zhuǎn)換為絕對吸附,又考慮了吸附相體積的影響,絕對吸附量曲線可能還存在下拐現(xiàn)象,即偏離Langmuir理論模型,這時原因很可能出在自由體積的測量上[5,20]。若自由體積測量值偏大,則會導(dǎo)致曲線下拐。容積法多次壓力點“充填—平衡—充填”實驗方法存在計算誤差累計問題,這個誤差是由自由體積多次參與計算累加所引起的[27]。對于隨壓力容積量產(chǎn)生較大變化的缸體,自由體積還要隨壓力修正,因此自由體積作為一個變量,其準(zhǔn)確測量一直是個難題。

    圖3 涪陵頁巖氣田JYA-2頁巖樣品等溫吸附曲線

    通常來說,先有數(shù)據(jù),然后建立理論模型;反過來,如果理論模型正確,也可以反推數(shù)據(jù),然后將反推的數(shù)據(jù)與實測數(shù)據(jù)比對來驗證理論模型。鑒于樣品缸中甲烷自由體積的測量難題,不妨認為頁巖等溫吸附滿足Langmuir理論模型,然后基于二元Langmuir方程,將甲烷初始自由體積設(shè)為未知參數(shù)V0,重新建立新的三元Langmuir方程,這里ρa采用給定值。

    (5)

    式中:f(V0)是一個含有V0的表達式,V0表示甲烷自由體積,cm3。

    利用最小二乘法、迭代技術(shù)或Matlab等軟件,擬合計算得到甲烷初始自由空間體積V0、VL和PL值。計算中還需考慮吸附相體積和缸體彈性變化對自由體積的影響,因而計算復(fù)雜,建議開發(fā)專用數(shù)據(jù)處理模板或軟件。圖4為JYD-2頁巖樣品,基于三元模型得到等溫吸附曲線。

    通過多塊頁巖樣品吸附實驗壓力數(shù)據(jù),利用三元Langmuir模型擬合計算出V0、VL、PL(表3,4)。比較表3中2組自由體積數(shù)據(jù):Langmuir模型計算的自由體積略低于氦氣測量的數(shù)據(jù),模型計算出的是甲烷的自由體積,氦氣測量的是近似甲烷的自由體積,理論上前者應(yīng)該低于后者。前者反映的是初始實驗壓力下甲烷的自由體積,后者是測試壓力下氦氣的自由體積(通常為3~5 MPa),兩者數(shù)據(jù)比對其相對差值很小(平均值0.66%)。前者采用的是Langmuir方程,使用甲烷氣體計算;后者采用的是氣體狀態(tài)方程,使用氦氣測量。采用不同的測試原理,使用不同的測試氣體,2種完全不同的方法有這么好的比對結(jié)果實屬不易。另外利用模型計算的自由體積來計算巖石骨架密度,將其與氦氣測量巖石骨架密度比較:從2種方法計算的巖石骨架密度來看,氦氣法理論上應(yīng)該大些,兩者相對差值在2%以內(nèi),這個結(jié)果完全可以接受。4組多樣品數(shù)據(jù)反映基于Langmuir吸附模型計算結(jié)果與實測結(jié)果非常接近,說明預(yù)設(shè)的頁巖等溫吸附模型正確,即頁巖等溫吸附符合單分子層理論(實驗壓力0~25 MPa)。

    圖4 涪陵頁巖氣田JYD-2頁巖樣品等溫吸附曲線

    表4 模型計算與舊方法實驗結(jié)果

    表3 模型計算與氦氣測量結(jié)果比對

    鑒于以上多塊樣品測試結(jié)果,表明三元Langmuir模型可以解決甲烷自由體積測試難題。表4中實驗樣品PX1至PX4為同一樣品分為4份,分別由ISO-300等溫吸附儀4組缸平行測量。舊方法計算的樣品Langmuir體積為1.60~3.67 cm3/g,最大與最小相差2.07 cm3/g;Langmuir壓力為4.33~6.92 MPa,最大與最小相差2.6 MPa。同樣實驗數(shù)據(jù)由模型計算,樣品Langmuir體積為1.61~2.00 cm3/g,最大與最小相差0.39 cm3/g;Langmuir壓力為2.32~3.81 MPa,最大與最小相差1.5 MPa。從平行樣的數(shù)據(jù)分布區(qū)間來看(圖5),模型計算處理的數(shù)據(jù)更集中,平行性更好,結(jié)果可信度更高。

    以上結(jié)果都是在吸附相密度ρa取定值0.421 g/cm3下獲得的?,F(xiàn)以涪陵龍馬溪組JYB-2頁巖樣品為例,研究吸附相密度改變對實驗結(jié)果的影響(圖6,表5)。圖6為吸附相密度與Langmuir體積的關(guān)系曲線,這里曲線存在一個變點,當(dāng)吸附相密度小于變點時,吸附相密度對Langmuir體積影響較大,密度越小,影響越大;當(dāng)吸附相密度大于變點時,吸附相密度對Langmuir體積影響非常小,即數(shù)學(xué)上尋找曲線由快速下降到突然轉(zhuǎn)緩趨于平衡的轉(zhuǎn)折點。這個變點就是該頁巖樣品吸附相密度的合理選取點。從表5看,若吸附相密度從0.3 g/cm3變化至1.2 g/cm3(為研究需要,將吸附相密度范圍擴大了些),其Langmuir體積逐漸減小,但變化不大。也就是說,實際工作中其實不需要花精力去尋找這個變點,依照前人傾向于吸附相密度可能在0.340~0.424 g/cm3之間的研究成果,無論采用哪個值都可以,Langmuir體積值變化不大,工程應(yīng)用上完全不存在問題。

    圖5 模型計算與舊方法實驗結(jié)果曲線

    圖6 涪陵頁巖氣田JYB-2頁巖樣品吸附相密度對Langmuir體積的影響

    這里也存在一個問題,就是對于一個實驗樣品來說,初始甲烷自由體積應(yīng)該與吸附相密度無關(guān)。只是由于給定的吸附相密度不同,模型擬合計算的自由體積不同。因?qū)嶋H中測量存在不確定度,自由體積的真值很難確定,但鑒于自由體積與吸附相密度的成對使用,使實驗結(jié)果Langmuir體積和Langmuir壓力變化不大,實驗結(jié)果一致性較好(表5)。通過應(yīng)用三元Langmuir模型,模擬計算甲烷初始自由體積,可以很好地校正下拐曲線和異常曲線。

    三元Langmuir模型解決了以前自由體積的測量難題,今后頁巖等溫吸附實驗可以不采用氦氣來測自由體積,這一實驗步驟完全可以取消;但建議保留此法,以便和模型計算結(jié)果有個比較參考。同理也可以用三元Langmuir模型和甲烷吸附實驗來計算頁巖的孔隙度,這個孔隙度是針對甲烷氣的,在工程上來說更合理。

    表5 涪陵頁巖氣田JYB-2頁巖樣品吸附相密度對Langmuir參數(shù)結(jié)果的影響

    3 結(jié)論

    (1)頁巖不同于煤巖,頁巖等溫吸附曲線下拐屬于正?,F(xiàn)象。頁巖實驗壓力高,高壓時需要考慮過剩吸附與絕對吸附的差別,需要將實驗測得的過剩吸附量轉(zhuǎn)化為絕對吸附量。

    (2)自由體積對吸附量的計算有重大影響,自由體積不是恒量,甲烷吸附在巖石顆粒表面占據(jù)了部分體積空間,吸附量計算時需要考慮吸附相體積的影響。對于隨壓力容量變化較大的缸體,考慮到器壁的壓變性,吸附量計算中也需加以校正。

    (3)利用氦氣來測量甲烷自由體積的方法存在缺陷,基于Langmuir理論,可將實驗初始甲烷自由體積設(shè)為未知參數(shù),建立新的三元Langmuir方程,通過數(shù)值迭代技術(shù)或Matlab軟件計算自由體積。

    (4)基于Langmuir理論模型計算的自由體積、巖石骨架密度與氦氣測量的自由體積、巖石骨架密度作比較,比較結(jié)果令人滿意,反映出頁巖中甲烷等溫吸附規(guī)律符合單分子層吸附理論,可以用Langmuir模型來計算吸附量或修正實驗曲線(本文實驗壓力0~25 MPa)。

    [1] 鮑云杰.從煤層含氣量測定技術(shù)發(fā)展看頁巖含氣性評價的發(fā)展方向[J].石油實驗地質(zhì),2014,36(6):762-766.

    BAO Yunjie.Prospect of shale gas evaluation in view of measurement techniques development of coal seam gas content[J].Petro-leum Geology & Experiment,2014,36(6):762-766.

    [2] U.S. Energy Information Administration.Technically recoverable shale oil and shale gas resources:An assessment of 137 shale formation in 41 countries outside the United States[EB/OL].Washington DC:U.S. Department of Energy,2013[2014-09-26].https://www.eia.gov/analysis/studies/worldshalegas/pdf/overview.pdf.

    [3] 張金川,金之鈞,袁明生.頁巖氣成藏機理和分布[J].天然氣工業(yè),2004,24(7):15-18.

    ZHANG Jinchuan,JIN Zhijun,YUAN Mingsheng.Reservoiring mechanism of shale gas and its distribution[J].Natural Gas Industry,2004,24(7):15-18.

    [4] 林臘梅,張金川,韓雙彪,等.泥頁巖儲層等溫吸附測試異常探討[J].油氣地質(zhì)與采收率,2012,19(6):30-32.

    LIN Lamei,ZHANG Jinchuan,HAN Shuangbiao,et al.Study on abnormal curves of isothermal adsorption of shale[J].Petroleum Geology and Recovery Efficiency,2012,19(6):30-32.

    [5] 聶海寬,張金川,馬曉彬,等.頁巖等溫吸附氣含量負吸附現(xiàn)象初探[J].地學(xué)前緣,2013,20(6):282-288.

    NIE Haikuan,ZHANG Jinchuan,MA Xiaobin,et al.A preliminary study of negative adsorption phenomena of shale adsorption gas content by isothermal adsorption[J].Earth Science Frontiers,2013,20(6):282-288.

    [6] 張慶玲.頁巖容量法等溫吸附實驗中異?,F(xiàn)象分析[J].煤田地質(zhì)與勘探,2015,43(5):31-33.

    ZHANG Qingling.The analysis of abnormal phenomena in shale isothermal absorption volumetry test[J].Coal Geology & Exploration,2015,43(5):31-33.

    [7] 俞凌杰,范明,陳紅宇,等.富有機質(zhì)頁巖高溫高壓重量法等溫吸附實驗[J].石油學(xué)報,2015,36(5):557-563.

    YU Lingjie,FAN Ming,CHEN Hongyu,et al.Isothermal adsorption experiment of organic-rich shale under high temperature and pressure using gravimetric method[J].Acta Petrolei Sinica,2015,36(5):557-563.

    [8] 中華人民共和國國家質(zhì)量監(jiān)督檢驗檢疫總局,中國國家標(biāo)準(zhǔn)化管理委員會.GB/T 19560—2008 煤的高壓等溫吸附試驗方法[S].北京:中國標(biāo)準(zhǔn)出版社,2009.

    General Administration of Quality Supervision,Inspection and Quarantine of the People’s Republic of China,Standardization Administration of China.GB/T19560—2008 Experimental method of high-pressure isothermal adsorption to coal[S].Beijing:China Standard Press,2009.

    [9] 周理,周亞平,孫艷,等.超臨界吸附及氣體代油燃料技術(shù)研究進展[J].自然科學(xué)進展,2004,14(6):615-623.

    ZHOU Li,ZHOU Yaping,SUN Yan,et al.Research advances of supercritical adsorption and gas as alternative fuels[J].Progress in Natural Science,2004,14(6):615-623.

    [10] 熊偉,郭為,劉洪林,等.頁巖的儲層特征以及等溫吸附特征[J].天然氣工業(yè),2012,32(1):113-116.

    XIONG Wei,GUO Wei,LIU Honglin,et al.Shale reservoir characteristics and isothermal adsorption properties[J].Natural Gas Industry,2012,32(1):113-116.

    [11] 張志英,楊盛波.頁巖氣吸附解吸規(guī)律研究[J].實驗力學(xué),2012,27(4):492-497.

    ZHANG Zhiying,YANG Shengbo.On the adsorption and desorption trend of shale gas[J].Journal of Experimental Mechanics,2012,27(4):492-497.

    [12] 閆建萍,張同偉,李艷芳,等.頁巖有機質(zhì)特征對甲烷吸附的影響[J].煤炭學(xué)報,2013,38(5):805-811.

    YAN Jianping,ZHANG Tongwei,LI Yanfang,et al.Effect of the organic matter characteristics on methane adsorption in shale[J].Journal of China Coal Society,2013,38(5):805-811.

    [13] 胡濤,馬正飛,姚虎卿.甲烷超臨界高壓吸附等溫線研究[J].天然氣化工,2002,27(2):36-40.

    HU Tao,MA Zhengfei,YAO Huqing.Study on high pressure adsorption isotherms of supercritical methane[J].Natural Gas Chemical Industry,2002,27(2):36-40.

    [14] 周尚文,薛華慶,郭偉,等.基于重量法的頁巖氣超臨界吸附特征實驗研究[J].煤炭學(xué)報,2016,41(11):2806-2812.

    ZHOU Shangwen,XUE Huaqing,GUO Wei,et al.Supercritical isothermal adsorption characteristics of shale gas based on gravimetric method[J].Journal of China Coal Society,2016,41(11):2806-2812.

    [15] DUBININ M M.The potential theory of adsorption of gases and vapors for adsorbents with energetically nonuniform surfaces[J].Chemical Reviews,1960,60(2):235-241.

    [16] HAYDEL J J,Kobayashi R.Adsorption equilibria in methane-propane-silica gel system at high pressures[J].Industrial & EngineeringChemistry Fundamentals,1967,6(4):546-554.

    [17] MENON P G.Adsorption at high pressures[J].Chemical Reviews,1968,68(3):277-294.

    [18] AMBROSE R J,HARTMAN R C,DIZA-CAMPOS M,et al.Shale gas-in-place calculations Part Ⅰ:New pore-scale considerations[J].Society of Petroleum Engineers,2010,17(1):219-229.

    [19] VERMESSE J,VIDAL D,MALBRUNOT P.Gas adsorption on zeolites at high pressure[J].Langmuir,1996,12(17):4190-4196.

    [20] ROSS D J K,BUSTIN R M.Impact of mass balance calculations on adsorption capacities in microporous shale gas reservoirs[J].Fuel,2007,86(17/18):2696-2706.

    [21] 熊健,劉向君,梁利喜.頁巖中超臨界甲烷等溫吸附模型研究[J].石油鉆探技術(shù),2015,43(3):96-102.

    XIONG Jian,LIU Xiangjun,LIANG Lixi.Isothermal adsorption model of supercritical methane in shale[J].Petroleum Drilling Techniques,2015,43(3):96-102.

    [22] GASPARIK M,GHANIZADEH A,BERTIER P,et al.High-pressure methane sorption isotherms of black shales from the Netherlands[J].Energy & Fuels,2012,26(8):4995-5004.

    [23] REXER T F T,BENHAM M J,APLIN A C,et al.Methane adsorption on shale under simulated geological temperature and pressure conditions[J].Energy & Fuels,2013,27(6):3099-3109.

    [24] 張慶玲,曹利戈.煤的等溫吸附測試中數(shù)據(jù)處理問題研究[J].煤炭學(xué)報,2003,28(2):131-135.

    ZHANG Qingling,CAO Lige.Study of data processing in coal sorption isotherm testing[J].Journal of China Coal Society,2003,28(2):131-135.

    [25] TAN Z M,GUBBINS K E.Adsorption in carbon micropores at supercritical temperatures[J].The Journal of Physical Chemistry,1990,94(15):6061-6069.

    [26] 郭為,熊偉,高樹生,等.溫度對頁巖等溫吸附/解吸特征影響[J].石油勘探與開發(fā),2013,40(4):481-485.

    GUO Wei,XIONG Wei,GAO Shusheng,et al.Impact of temperature on the isothermal adsorption/desorption characteristics of shale gas[J].Petroleum Exploration and Development,2013,40(4):481-485.

    [27] 何斌,寧正福,楊峰,等.頁巖等溫吸附實驗及實驗誤差分析[J].煤炭學(xué)報,2015,40(S1):177-184.

    HE Bin,NING Zhengfu,YANG Feng,et al.Shale isothermal adsorption experiment and experimental error analysis[J].Journal of China Coal Society,2015,40(S1):177-184.

    猜你喜歡
    等溫甲烷頁巖
    液氧甲烷發(fā)動機
    EPDM/PP基TPV非等溫結(jié)晶行為的研究
    論煤炭運輸之甲烷爆炸
    水上消防(2020年1期)2020-07-24 09:26:02
    Gas from human waste
    頁巖氣開發(fā)降溫
    能源(2016年1期)2016-12-01 05:10:02
    快速檢測豬鏈球菌的環(huán)介導(dǎo)等溫擴增方法
    納米CaCO3對FEP非等溫結(jié)晶動力學(xué)的影響
    中國塑料(2015年3期)2015-11-27 03:41:54
    酰胺類β成核劑對PP-R非等溫結(jié)晶動力學(xué)影響
    中國塑料(2014年2期)2014-10-17 02:50:59
    基于VB6.0的紅外甲烷檢測儀串行通信的實現(xiàn)
    我國頁巖氣可采資源量初步估計為31萬億m3
    中国国产av一级| 老司机午夜十八禁免费视频| 日本撒尿小便嘘嘘汇集6| 成人18禁高潮啪啪吃奶动态图| 一本大道久久a久久精品| 啦啦啦啦在线视频资源| 免费少妇av软件| 99re6热这里在线精品视频| 麻豆av在线久日| 亚洲免费av在线视频| 淫妇啪啪啪对白视频 | 久久久国产欧美日韩av| 中文字幕人妻丝袜一区二区| 亚洲人成电影观看| 爱豆传媒免费全集在线观看| 亚洲av成人一区二区三| 伊人久久大香线蕉亚洲五| 欧美黑人欧美精品刺激| 久久免费观看电影| 亚洲综合色网址| 在线观看舔阴道视频| 欧美日韩黄片免| 丝袜喷水一区| 国产成人精品在线电影| 大片电影免费在线观看免费| 久久精品亚洲熟妇少妇任你| 国产成人系列免费观看| 免费高清在线观看日韩| 视频区图区小说| 亚洲精品av麻豆狂野| 性色av一级| 国产精品久久久人人做人人爽| 国产成人精品久久二区二区91| 夜夜夜夜夜久久久久| xxxhd国产人妻xxx| 亚洲熟女毛片儿| 国产日韩欧美在线精品| 久久人妻熟女aⅴ| 亚洲免费av在线视频| netflix在线观看网站| 成人国产一区最新在线观看| 日韩,欧美,国产一区二区三区| a在线观看视频网站| 正在播放国产对白刺激| 午夜两性在线视频| 欧美97在线视频| 中文字幕另类日韩欧美亚洲嫩草| e午夜精品久久久久久久| 一区二区三区四区激情视频| 妹子高潮喷水视频| 18禁国产床啪视频网站| 国产精品成人在线| 我的亚洲天堂| 婷婷成人精品国产| 一二三四社区在线视频社区8| 一区二区三区精品91| 久久午夜综合久久蜜桃| 日韩中文字幕欧美一区二区| 亚洲精品久久成人aⅴ小说| 丝袜美腿诱惑在线| 三级毛片av免费| 真人做人爱边吃奶动态| av线在线观看网站| 久久女婷五月综合色啪小说| 国产成人免费观看mmmm| 亚洲自偷自拍图片 自拍| 日本91视频免费播放| 肉色欧美久久久久久久蜜桃| 人妻久久中文字幕网| 欧美性长视频在线观看| 日韩 欧美 亚洲 中文字幕| 国产亚洲精品久久久久5区| 国产区一区二久久| 欧美人与性动交α欧美软件| 欧美在线黄色| 纵有疾风起免费观看全集完整版| 国产精品亚洲av一区麻豆| 亚洲专区字幕在线| 国产淫语在线视频| 人成视频在线观看免费观看| www.精华液| 一二三四在线观看免费中文在| 在线永久观看黄色视频| 欧美亚洲 丝袜 人妻 在线| 久久影院123| 国产精品一二三区在线看| 99re6热这里在线精品视频| 午夜精品久久久久久毛片777| 国精品久久久久久国模美| 午夜精品国产一区二区电影| 丝袜美腿诱惑在线| 免费久久久久久久精品成人欧美视频| 又紧又爽又黄一区二区| 日本精品一区二区三区蜜桃| 18禁黄网站禁片午夜丰满| 午夜日韩欧美国产| 精品亚洲成a人片在线观看| 国产精品av久久久久免费| 免费不卡黄色视频| 亚洲精品国产av成人精品| 亚洲第一av免费看| 久久免费观看电影| 我的亚洲天堂| av天堂在线播放| 精品国产一区二区三区久久久樱花| 国产区一区二久久| 丝袜脚勾引网站| av视频免费观看在线观看| a在线观看视频网站| 亚洲欧美日韩高清在线视频 | 午夜福利一区二区在线看| 久久久久国内视频| 啦啦啦视频在线资源免费观看| 嫩草影视91久久| 香蕉丝袜av| 最新在线观看一区二区三区| 97人妻天天添夜夜摸| 国产91精品成人一区二区三区 | 国产一区二区 视频在线| 蜜桃国产av成人99| 午夜福利在线免费观看网站| av又黄又爽大尺度在线免费看| 九色亚洲精品在线播放| 蜜桃在线观看..| 国产熟女午夜一区二区三区| 丰满少妇做爰视频| 久久精品熟女亚洲av麻豆精品| 啪啪无遮挡十八禁网站| 黄片播放在线免费| 日本五十路高清| 人人妻人人澡人人爽人人夜夜| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜福利视频精品| 丝袜喷水一区| 超色免费av| 精品少妇久久久久久888优播| 亚洲专区中文字幕在线| 性色av乱码一区二区三区2| 国产精品免费大片| 妹子高潮喷水视频| 老司机影院毛片| 91麻豆av在线| 午夜两性在线视频| 久久国产精品大桥未久av| 久久久久国产一级毛片高清牌| 日日摸夜夜添夜夜添小说| av在线app专区| 国产福利在线免费观看视频| 精品一区二区三区av网在线观看 | 亚洲少妇的诱惑av| 欧美在线一区亚洲| 久久人人爽av亚洲精品天堂| av天堂在线播放| 天天操日日干夜夜撸| 国产精品一区二区在线不卡| 日韩制服丝袜自拍偷拍| 亚洲欧美精品自产自拍| 久久综合国产亚洲精品| 色94色欧美一区二区| 亚洲精品一区蜜桃| 香蕉国产在线看| 亚洲综合色网址| 亚洲专区字幕在线| 老司机深夜福利视频在线观看 | 久久99热这里只频精品6学生| 黑人操中国人逼视频| 亚洲av成人不卡在线观看播放网 | 久久久久视频综合| 女人爽到高潮嗷嗷叫在线视频| 丝袜人妻中文字幕| 亚洲伊人久久精品综合| 高清av免费在线| 看免费av毛片| 天天躁夜夜躁狠狠躁躁| 男女床上黄色一级片免费看| 成人黄色视频免费在线看| 脱女人内裤的视频| 欧美+亚洲+日韩+国产| 大香蕉久久网| 69精品国产乱码久久久| 午夜免费鲁丝| 水蜜桃什么品种好| 成人av一区二区三区在线看 | 肉色欧美久久久久久久蜜桃| 超色免费av| 老司机福利观看| 欧美一级毛片孕妇| 99九九在线精品视频| 国产成人欧美| 午夜福利一区二区在线看| 国产伦人伦偷精品视频| 午夜老司机福利片| 青春草亚洲视频在线观看| 亚洲专区国产一区二区| 亚洲精品中文字幕在线视频| 国产又爽黄色视频| 叶爱在线成人免费视频播放| 国产极品粉嫩免费观看在线| 黑人欧美特级aaaaaa片| 国产精品免费视频内射| 久久亚洲精品不卡| 国产精品av久久久久免费| 免费久久久久久久精品成人欧美视频| 女人被躁到高潮嗷嗷叫费观| av超薄肉色丝袜交足视频| av在线播放精品| 精品一区二区三卡| 99香蕉大伊视频| a级毛片黄视频| 成年人免费黄色播放视频| 狂野欧美激情性xxxx| 亚洲三区欧美一区| 国产精品自产拍在线观看55亚洲 | 国产精品秋霞免费鲁丝片| 日韩 亚洲 欧美在线| 国产亚洲欧美在线一区二区| 久久ye,这里只有精品| av天堂在线播放| 又大又爽又粗| 夜夜骑夜夜射夜夜干| 在线精品无人区一区二区三| 亚洲精品乱久久久久久| 国产高清视频在线播放一区 | 亚洲成人免费电影在线观看| 久久人妻熟女aⅴ| 午夜免费鲁丝| 国产主播在线观看一区二区| 少妇人妻久久综合中文| 亚洲精品在线美女| 亚洲九九香蕉| 精品少妇久久久久久888优播| 男女午夜视频在线观看| 一区二区三区精品91| 国产成+人综合+亚洲专区| 亚洲,欧美精品.| 久久久久视频综合| 久久人人97超碰香蕉20202| 国产福利在线免费观看视频| 黄色怎么调成土黄色| 亚洲欧美色中文字幕在线| 国产黄频视频在线观看| 亚洲午夜精品一区,二区,三区| 国产亚洲欧美在线一区二区| 久久午夜综合久久蜜桃| 青春草视频在线免费观看| 日韩有码中文字幕| 国产精品一二三区在线看| 国产免费福利视频在线观看| 91精品国产国语对白视频| 精品一品国产午夜福利视频| 国产精品秋霞免费鲁丝片| 国产日韩欧美视频二区| 最近最新中文字幕大全免费视频| 黑人操中国人逼视频| 国产伦理片在线播放av一区| 美女扒开内裤让男人捅视频| 午夜福利视频在线观看免费| 日本一区二区免费在线视频| 韩国高清视频一区二区三区| 最近最新免费中文字幕在线| 欧美人与性动交α欧美软件| 国产一区二区三区在线臀色熟女 | 视频区欧美日本亚洲| 精品一区在线观看国产| 男人添女人高潮全过程视频| 久久久久久亚洲精品国产蜜桃av| 老熟妇仑乱视频hdxx| 免费女性裸体啪啪无遮挡网站| 精品国产乱子伦一区二区三区 | 国产一区二区在线观看av| 一区二区三区激情视频| 一级,二级,三级黄色视频| 国产在视频线精品| 19禁男女啪啪无遮挡网站| 黑人猛操日本美女一级片| 人人妻,人人澡人人爽秒播| 视频在线观看一区二区三区| 高清欧美精品videossex| www.av在线官网国产| √禁漫天堂资源中文www| 成人国产一区最新在线观看| 99香蕉大伊视频| 下体分泌物呈黄色| 久久毛片免费看一区二区三区| 国产av精品麻豆| 国产一级毛片在线| 69精品国产乱码久久久| 一本色道久久久久久精品综合| 久久久精品94久久精品| 两性午夜刺激爽爽歪歪视频在线观看 | 少妇人妻久久综合中文| 成人影院久久| 亚洲成人手机| 国产一级毛片在线| 日韩三级视频一区二区三区| 成人黄色视频免费在线看| 久久影院123| 亚洲精品成人av观看孕妇| 人人妻人人澡人人爽人人夜夜| 一级毛片电影观看| 国产精品麻豆人妻色哟哟久久| 国产黄频视频在线观看| 曰老女人黄片| 国产三级黄色录像| 成人18禁高潮啪啪吃奶动态图| 波多野结衣av一区二区av| 久久香蕉激情| 777久久人妻少妇嫩草av网站| 日韩熟女老妇一区二区性免费视频| 亚洲avbb在线观看| 动漫黄色视频在线观看| 又大又爽又粗| 亚洲国产精品成人久久小说| 国产在线免费精品| 国产高清国产精品国产三级| 巨乳人妻的诱惑在线观看| 后天国语完整版免费观看| 不卡av一区二区三区| 日韩,欧美,国产一区二区三区| av一本久久久久| 久久狼人影院| 美女扒开内裤让男人捅视频| 亚洲欧美色中文字幕在线| 波多野结衣一区麻豆| 久久久久网色| 日韩欧美免费精品| av一本久久久久| 啦啦啦啦在线视频资源| 欧美激情久久久久久爽电影 | 丝袜美足系列| 国产深夜福利视频在线观看| 国产男女超爽视频在线观看| 免费久久久久久久精品成人欧美视频| 男人爽女人下面视频在线观看| 亚洲va日本ⅴa欧美va伊人久久 | 免费观看人在逋| 熟女少妇亚洲综合色aaa.| 亚洲午夜精品一区,二区,三区| 国产91精品成人一区二区三区 | 可以免费在线观看a视频的电影网站| 999久久久精品免费观看国产| 汤姆久久久久久久影院中文字幕| 麻豆乱淫一区二区| 十八禁人妻一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美精品一区二区大全| 一区二区av电影网| 啦啦啦视频在线资源免费观看| 三上悠亚av全集在线观看| 欧美日韩亚洲综合一区二区三区_| 亚洲综合色网址| 欧美成人午夜精品| 亚洲国产日韩一区二区| 老司机午夜十八禁免费视频| 97人妻天天添夜夜摸| 国产欧美日韩精品亚洲av| 亚洲人成电影观看| 99国产精品99久久久久| 青春草视频在线免费观看| 亚洲美女黄色视频免费看| 丰满少妇做爰视频| 精品乱码久久久久久99久播| 夜夜夜夜夜久久久久| 久久久久网色| 久久中文看片网| 国产成人av教育| 一边摸一边做爽爽视频免费| 久久久国产精品麻豆| 亚洲精品国产色婷婷电影| 久久人人97超碰香蕉20202| 免费观看av网站的网址| 日韩中文字幕视频在线看片| 一区二区三区乱码不卡18| 亚洲精品一区蜜桃| 成年美女黄网站色视频大全免费| 精品第一国产精品| 在线观看人妻少妇| 亚洲美女黄色视频免费看| 亚洲一区中文字幕在线| videosex国产| 9色porny在线观看| 久热这里只有精品99| 可以免费在线观看a视频的电影网站| 中文字幕制服av| 免费高清在线观看视频在线观看| 免费在线观看视频国产中文字幕亚洲 | 精品亚洲成国产av| 成年人黄色毛片网站| 亚洲欧美一区二区三区久久| 亚洲伊人色综图| 麻豆国产av国片精品| 日韩中文字幕欧美一区二区| 美女大奶头黄色视频| 国产成人精品久久二区二区免费| 久久中文看片网| 一级片'在线观看视频| 韩国高清视频一区二区三区| 国产成人精品无人区| 美女福利国产在线| 人妻 亚洲 视频| 少妇的丰满在线观看| 国产成人欧美| 超色免费av| 精品福利观看| 天天添夜夜摸| 亚洲欧美日韩高清在线视频 | 成人国产av品久久久| 男女高潮啪啪啪动态图| 99久久国产精品久久久| 国产97色在线日韩免费| 曰老女人黄片| 久久久国产精品麻豆| 法律面前人人平等表现在哪些方面 | av线在线观看网站| 国产精品 国内视频| 黑人操中国人逼视频| 美女高潮喷水抽搐中文字幕| 久久久久久亚洲精品国产蜜桃av| 999久久久精品免费观看国产| 亚洲欧美一区二区三区黑人| 50天的宝宝边吃奶边哭怎么回事| 久久久国产欧美日韩av| 一区二区三区精品91| 大片免费播放器 马上看| 亚洲精品一区蜜桃| 免费看十八禁软件| 一区二区三区乱码不卡18| 亚洲久久久国产精品| 亚洲国产中文字幕在线视频| 亚洲国产成人一精品久久久| 电影成人av| 99久久综合免费| 亚洲久久久国产精品| 亚洲色图综合在线观看| 亚洲人成77777在线视频| 国产在视频线精品| 成人免费观看视频高清| 久久ye,这里只有精品| 人妻 亚洲 视频| 久久精品久久久久久噜噜老黄| 法律面前人人平等表现在哪些方面 | 午夜免费观看性视频| 国产男女内射视频| 伦理电影免费视频| 在线观看免费视频网站a站| 国产av又大| 狂野欧美激情性bbbbbb| 国产成人精品久久二区二区91| 肉色欧美久久久久久久蜜桃| 久久精品国产亚洲av高清一级| 日日夜夜操网爽| 80岁老熟妇乱子伦牲交| 人妻久久中文字幕网| 日韩欧美国产一区二区入口| tocl精华| 国产精品免费视频内射| 久久99一区二区三区| av电影中文网址| 国产免费现黄频在线看| 91老司机精品| 他把我摸到了高潮在线观看 | 国产精品秋霞免费鲁丝片| 国产一级毛片在线| 亚洲精品久久久久久婷婷小说| 老熟妇乱子伦视频在线观看 | 大码成人一级视频| av在线app专区| 18禁黄网站禁片午夜丰满| 婷婷色av中文字幕| 这个男人来自地球电影免费观看| 在线十欧美十亚洲十日本专区| 久久人妻福利社区极品人妻图片| 亚洲国产欧美网| www.av在线官网国产| 国产亚洲精品一区二区www | 在线看a的网站| 黑人巨大精品欧美一区二区mp4| 如日韩欧美国产精品一区二区三区| 好男人电影高清在线观看| 爱豆传媒免费全集在线观看| 久久亚洲国产成人精品v| 国产av国产精品国产| 久久人妻福利社区极品人妻图片| 又紧又爽又黄一区二区| 女警被强在线播放| 久久狼人影院| 看免费av毛片| 性色av乱码一区二区三区2| 亚洲精品成人av观看孕妇| 亚洲人成电影免费在线| 日本欧美视频一区| 精品一区二区三区av网在线观看 | 亚洲精品日韩在线中文字幕| 亚洲全国av大片| 777米奇影视久久| 亚洲第一青青草原| 国产免费一区二区三区四区乱码| 亚洲精华国产精华精| 久久久久久久大尺度免费视频| 18禁黄网站禁片午夜丰满| 欧美97在线视频| 50天的宝宝边吃奶边哭怎么回事| 99精国产麻豆久久婷婷| 黑人巨大精品欧美一区二区mp4| 精品熟女少妇八av免费久了| 男人添女人高潮全过程视频| 一级毛片女人18水好多| 国产一区二区三区在线臀色熟女 | 视频区图区小说| 欧美在线黄色| 精品国产乱子伦一区二区三区 | 黄网站色视频无遮挡免费观看| 三级毛片av免费| 成年动漫av网址| av天堂久久9| 国产精品一区二区在线观看99| 老司机在亚洲福利影院| 欧美日韩一级在线毛片| 免费黄频网站在线观看国产| 亚洲情色 制服丝袜| 国产成人精品在线电影| 国产精品99久久99久久久不卡| 亚洲精品久久久久久婷婷小说| 国产精品影院久久| 国产高清videossex| 欧美日韩中文字幕国产精品一区二区三区 | av国产精品久久久久影院| 久久中文字幕一级| 久久久精品区二区三区| 国产精品亚洲av一区麻豆| 久久亚洲国产成人精品v| 精品久久久精品久久久| 在线观看免费日韩欧美大片| 97在线人人人人妻| 成人18禁高潮啪啪吃奶动态图| 精品国产一区二区三区久久久樱花| 亚洲专区国产一区二区| 菩萨蛮人人尽说江南好唐韦庄| 两个人免费观看高清视频| 亚洲 欧美一区二区三区| 亚洲欧美精品自产自拍| 国产一级毛片在线| 2018国产大陆天天弄谢| 女人久久www免费人成看片| 9热在线视频观看99| 久久香蕉激情| 久久亚洲精品不卡| 亚洲中文日韩欧美视频| 18在线观看网站| 女人高潮潮喷娇喘18禁视频| 久久影院123| 日韩欧美免费精品| 婷婷色av中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 91精品三级在线观看| www日本在线高清视频| 国产熟女午夜一区二区三区| 国产真人三级小视频在线观看| 精品国内亚洲2022精品成人 | 少妇被粗大的猛进出69影院| 国产免费一区二区三区四区乱码| 18在线观看网站| 免费黄频网站在线观看国产| 久久久久久久精品精品| 黄色视频不卡| 欧美乱码精品一区二区三区| 精品免费久久久久久久清纯 | 欧美日韩视频精品一区| 国产视频一区二区在线看| 丝袜人妻中文字幕| 欧美成人午夜精品| 波多野结衣一区麻豆| 亚洲全国av大片| 国产区一区二久久| 美女午夜性视频免费| svipshipincom国产片| 宅男免费午夜| 啦啦啦中文免费视频观看日本| 亚洲av电影在线进入| 又黄又粗又硬又大视频| 老鸭窝网址在线观看| 国产视频一区二区在线看| 久久香蕉激情| 精品亚洲成国产av| 高清av免费在线| 午夜福利,免费看| 日日爽夜夜爽网站| 久久人妻福利社区极品人妻图片| 9热在线视频观看99| 免费人妻精品一区二区三区视频| 午夜精品国产一区二区电影| 国产区一区二久久| av片东京热男人的天堂| 午夜影院在线不卡| 国产成人免费无遮挡视频| 最新的欧美精品一区二区| 波多野结衣av一区二区av| 啦啦啦中文免费视频观看日本| 久久免费观看电影| 国产亚洲精品第一综合不卡| 免费在线观看日本一区| 男女边摸边吃奶| 亚洲少妇的诱惑av| 久久99热这里只频精品6学生| 色精品久久人妻99蜜桃| 国产亚洲av片在线观看秒播厂| av在线播放精品| 国产黄色免费在线视频| 亚洲第一av免费看| 窝窝影院91人妻| 亚洲精品日韩在线中文字幕| 无遮挡黄片免费观看| 丁香六月欧美| 桃红色精品国产亚洲av| 欧美97在线视频|