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

    超臨界態(tài)甲烷密度研究

    2021-03-31 00:03:26張新賓宋黨育李云波齊永安范順利
    煤田地質(zhì)與勘探 2021年1期
    關(guān)鍵詞:游離態(tài)狀態(tài)方程超臨界

    張新賓,宋黨育,李云波,齊永安,范順利

    超臨界態(tài)甲烷密度研究

    張新賓1,2,宋黨育1,2,李云波1,2,齊永安1,2,范順利1,2

    (1. 河南理工大學(xué) 資源環(huán)境學(xué)院,河南 焦作 454003;2. 中原經(jīng)濟(jì)區(qū)煤層(頁巖)氣河南省協(xié)同創(chuàng)新中心,河南 焦作 454003)

    超臨界甲烷密度不僅是非常規(guī)天然氣資源量計(jì)算與吸附能力測定的重要參數(shù),同時(shí)也是衡量超臨界甲烷擴(kuò)散效率與溶解有機(jī)物能力的主要指標(biāo)。通過對比分析各種氣體狀態(tài)方程的適用性,認(rèn)為基于亥姆霍茲能量基本狀態(tài)方程可以準(zhǔn)確計(jì)算0~30 MPa、270~360 K條件下甲烷的密度。利用Microsoft Office Excel編寫了甲烷密度的計(jì)算程序,與NIST(美國國家標(biāo)準(zhǔn)與技術(shù)研究院)商業(yè)軟件計(jì)算結(jié)果相比,誤差小于0.05%。分析了0~100 MPa、270~360 K范圍內(nèi)甲烷密度的變化規(guī)律。結(jié)果表明,甲烷密度隨壓力增大而增大,在低于30 MPa時(shí)增速較大且對溫度的敏感性較強(qiáng),高于30 MPa時(shí)增速逐漸變緩,且敏感性減弱。在煤層原位條件下隨著埋深的增大,甲烷密度隨溫度升高而減小,隨壓力增大而增大;在溫度與壓力共同作用下,甲烷密度呈先增速不變、近似線性增加,后增速逐漸減小、凸曲線形增加的變化規(guī)律。游離態(tài)甲烷密度受溫度的影響比吸附態(tài)甲烷小,是深部煤層氣資源增量的主要貢獻(xiàn)者。研究結(jié)果為深部煤層氣賦存及其潛力預(yù)測提供了基礎(chǔ)參數(shù)。

    超臨界甲烷;甲烷密度;氣體狀態(tài)方程;亥姆霍茲能量;煤層氣資源量

    超臨界流體是指溫度和壓力都超過臨界點(diǎn)時(shí)的流體。經(jīng)典超臨界流體理論認(rèn)為,超臨界流體是結(jié)構(gòu)均勻的單相流體,其蒸發(fā)潛熱和表面張力消失,不再發(fā)生相變,也不存在氣液共存界限[1]。嚴(yán)格意義上講超臨界流體既不是液體也不是氣體,在有的溫壓區(qū)間內(nèi)具有氣體的高擴(kuò)散性,有的溫壓條件下具有液體的溶解性,同時(shí)在臨界點(diǎn)附近溫度或壓力的微小變化可能造成密度的急劇改變。因此,密度是判斷超臨界流體性質(zhì)的主要依據(jù)[2]。

    甲烷的臨界溫度是–82.6℃,臨界壓力是4.6 MPa,當(dāng)溫度和壓力都超過臨界值時(shí)即為超臨界態(tài)[3]。在正常地溫與壓力梯度下,當(dāng)埋藏深度超過500 m時(shí)煤層中的甲烷已達(dá)到超臨界狀態(tài)。但是由于儲層溫度遠(yuǎn)遠(yuǎn)超過甲烷的臨界溫度,當(dāng)儲層壓力超過甲烷的超臨界壓力后,甲烷的密度是否會發(fā)生明顯變化,超臨界態(tài)甲烷對煤層的吸附能力會產(chǎn)生哪些影響?前人對這些問題并沒有非常明確的認(rèn)識。嚴(yán)格來說超臨界態(tài)甲烷的密度及其隨壓力、溫度的變化規(guī)律與理想氣體明顯不同[4],在深部非常規(guī)天然氣資源評價(jià)中,游離態(tài)甲烷資源量的計(jì)算需要準(zhǔn)確的甲烷密度參數(shù)。

    1978年,國際純粹與應(yīng)用化學(xué)聯(lián)合會(IUPAC)提出了基于狀態(tài)方程推導(dǎo)超臨界甲烷熱物性的方法,同時(shí)列出了部分溫度與壓力下甲烷的熱力學(xué)數(shù)據(jù)[5]。理想氣體狀態(tài)方程主要是在假設(shè)體積不變的情況下推導(dǎo)出密度公式,然而理想氣體狀態(tài)方程無法描述液相區(qū)和密度較高的氣相區(qū)的性質(zhì)[6]。為得到真實(shí)氣體的狀態(tài)方程,可對理想氣體狀態(tài)方程進(jìn)行修正,數(shù)學(xué)上最簡單的方式是定義壓縮因子[7-8]。在甲烷高壓等溫吸附實(shí)驗(yàn)中目前采用固定的溫度、壓力點(diǎn)下甲烷的壓縮系數(shù)校正自由體積中的甲烷量[9],與理想氣體狀態(tài)方程和利用壓縮因子得到的真實(shí)氣體狀態(tài)方程相比,R-K方程在計(jì)算甲烷密度時(shí)的應(yīng)用范圍和計(jì)算精度方面都有了一定的提高,但其在高壓下和預(yù)計(jì)物質(zhì)的液相密度時(shí)誤差較大[10]。如果可以快速計(jì)算不同溫度、壓力下游離態(tài)甲烷的密度,可以更加準(zhǔn)確地測定煤的吸附能力。目前可以從一些文獻(xiàn)中查閱部分壓力和溫度下甲烷密度,也可以利用一些商業(yè)軟件計(jì)算不同溫度和壓力下甲烷的密度和其他熱力學(xué)參數(shù),如ASPEN和REFPROP軟件,分別為美國AspenTech公司和NIST(美國國家標(biāo)準(zhǔn)與技術(shù)研究院)開發(fā)的商業(yè)軟件。也有學(xué)者通過實(shí)驗(yàn)測定了不同溫度、壓力下超臨界甲烷密度[11-15]。從行業(yè)應(yīng)用角度而言,需要根據(jù)超臨界甲烷密度隨溫度與壓力的變化規(guī)律,快速計(jì)算不同溫度、壓力下超臨界甲烷密度。

    論文基于亥姆霍茲能量守恒定律,利用亥姆霍茲能量形式的基本方程,通過將無量綱亥姆霍茲能量分解為理想氣體行為和殘余流體行為兩部分,采用優(yōu)化經(jīng)驗(yàn)熱力學(xué)方程結(jié)構(gòu)的方法,建立了270~360 K、0~30 MPa的溫度與壓力區(qū)間內(nèi)甲烷密度的計(jì)算方法。利用Excel軟件,采用迭代法、單變量求解和宏操作,編制了不同溫度、壓力下超臨界甲烷密度的直接計(jì)算程序,在此基礎(chǔ)上分析超臨界甲烷密度隨溫度與壓力的變化規(guī)律。

    1 特殊區(qū)域氣體狀態(tài)方程

    在化工和能源技術(shù)領(lǐng)域,工藝和設(shè)備的設(shè)計(jì)需要準(zhǔn)確的流體熱物理參數(shù)。然而,由于統(tǒng)計(jì)熱力學(xué)在理論上還沒有達(dá)到足夠的精確與簡捷,在大多數(shù)情況下人們必須使用經(jīng)驗(yàn)公式,這些公式是基于覆蓋整個(gè)流體區(qū)域的相應(yīng)精確實(shí)驗(yàn)數(shù)據(jù)得出的。1989年,U. Setzmann等[16]提出了優(yōu)化經(jīng)驗(yàn)熱力學(xué)方程結(jié)構(gòu)的新思路,在考慮物理依賴性的綜合函數(shù)表達(dá)式,即所謂“數(shù)據(jù)庫”的基礎(chǔ)上,同時(shí)優(yōu)化了方程的結(jié)構(gòu)和長度。

    特殊區(qū)域氣體狀態(tài)方程是基于亥姆霍茲能量形式的基本方程而建立的新狀態(tài)方程[17]。在熱力學(xué)中,亥姆霍茲能是一個(gè)熱力學(xué)勢,用在恒定的溫度和體積下從封閉熱力系統(tǒng)能得到最大“有用”功。對于這樣一個(gè)系統(tǒng),亥姆霍茲能差值的負(fù)數(shù)等于溫度和體積保持不變下可逆等溫過程的最大功輸出,即=–。新狀態(tài)方程采用將其轉(zhuǎn)化為無量綱形式的亥姆霍茲能量=/(),即為無量綱形式的亥姆霍茲能量??蓪⑵浞纸鉃槔硐霘怏w行為o和殘余流體行為r兩部分:

    理想氣體部分的無量綱亥姆霍茲能量與經(jīng)驗(yàn)方程關(guān)聯(lián):

    式(1)中的重點(diǎn)是確定亥姆霍茲能量殘余部分的函數(shù)形式。首先是建立方程的通用表達(dá)式,該表達(dá)式起著“數(shù)據(jù)庫”的作用。為了使最終方程盡可能簡單,“數(shù)據(jù)庫”以和表示為一個(gè)簡單的多項(xiàng)式:

    式中:d的范圍是1~9,寬度為1;t的范圍是–0.5~10,寬度為0.25。

    這意味著整個(gè)“數(shù)據(jù)庫”包含387個(gè)項(xiàng)。然后再根據(jù)式(3),采用一種特殊優(yōu)化程序確定有限項(xiàng)的最佳組合,將線性優(yōu)化過程與所有數(shù)據(jù)(線性和非線性)的非線性擬合迭代組合[18],最終得到亥姆霍茲能量殘余部分的函數(shù)形式[17]:

    式中:n為方程的系數(shù);d、t為指數(shù)。

    這樣就建成了270~360 K、0~30 MPa條件下的有效狀態(tài)方程:

    wj,k為方差為Q的零均值白噪聲序列;ζk為變速因子;Φk為氣隙磁量;φu,k∈Φk,為第u個(gè)導(dǎo)磁爪與磁極間的氣隙磁通量。

    與方程推導(dǎo)過程中采用的實(shí)驗(yàn)數(shù)據(jù)集相比,利用特殊區(qū)域氣體狀態(tài)方程計(jì)算得到的甲烷密度的不確定度為:0<≤12 MPa時(shí),Δ/≤±0.02%;12<≤30 MPa時(shí),Δ/≤±0.05%。在范圍和精度方面都滿足了超臨界甲烷密度的計(jì)算要求。

    2 超臨界甲烷密度的計(jì)算

    考慮到超臨界甲烷的溫度、壓力范圍和對密度計(jì)算精確度的要求,通過對比文獻(xiàn)中不同的狀態(tài)方程,認(rèn)為采用亥姆霍茲能量方程與熱力學(xué)參數(shù)的關(guān)系式可以更加準(zhǔn)確地計(jì)算超臨界甲烷密度。

    2.1 計(jì)算方法與過程

    根據(jù)亥姆霍茲能量方程與熱力學(xué)參數(shù)的關(guān)系式,壓力與無量綱形式的亥姆霍茲能量具有關(guān)系:

    對式(4)求偏導(dǎo)得:

    聯(lián)立式(6)、式(7)得:

    式中:c取值162.66 kg/m3;c取值190.564 K;m為比氣體常數(shù),J/(kg·K)。其中m=/,為氣體的摩爾質(zhì)量,g/mol;n、d、t取值見表1。

    將對應(yīng)參數(shù)代入式(8)后即可計(jì)算得到特定溫度和壓力下的甲烷密度。

    表1 甲烷狀態(tài)方程中的參數(shù)ni、di、ti取值

    2.2 計(jì)算實(shí)現(xiàn)過程

    式(8)可轉(zhuǎn)化為一元九次方程?;诎⒇悹柖ɡ?即各項(xiàng)系數(shù)經(jīng)過有限次四則運(yùn)算和乘方與開方運(yùn)算無法求解5次及以上的一元高次方程),對于一元九次方程并沒有通用的代數(shù)解法和求根公式,可以采用迭代法求近似解的方式求解。迭代法也稱輾轉(zhuǎn)法,是一種不斷用變量舊值遞推新值的過程。迭代算法是用計(jì)算機(jī)解決問題的一種基本方法,利用計(jì)算機(jī)運(yùn)算速度快、適合做重復(fù)性操作的特點(diǎn),對一組指令或一定步驟重復(fù)執(zhí)行,在每次執(zhí)行這組指令或步驟時(shí),都從變量的原值推出一個(gè)新值[19]。

    為簡化計(jì)算過程,同時(shí)考慮到Microsoft Office Excel(簡稱Excel)強(qiáng)大的計(jì)算能力和便捷性,本文采用其中的單變量求解功能求解方程。在Excel運(yùn)算中,單變量求解是根據(jù)目標(biāo)值不斷調(diào)整變量,直至計(jì)算結(jié)果與目標(biāo)值無限接近時(shí),才能確定變量值[20]。

    將式(8)轉(zhuǎn)化為:

    甲烷密度計(jì)算方法如圖1所示,計(jì)算步驟為:①首先,在Excel中輸入溫度、壓力常數(shù)與相關(guān)計(jì)算參數(shù)以便方程引用;②選擇任一單元格,輸入任意數(shù)值作為變量密度的初始值;③選擇任一單元格輸入式(9),公式必須引用密度的初始值;④設(shè)置迭代計(jì)算參數(shù),如迭代次數(shù)和最大誤差;⑤選擇單變量求解功能,在目標(biāo)單元格中引用輸入公式的單元格,并將目標(biāo)值設(shè)置為0,在可變單元格編輯框中引用輸入變量密度的單元格,即可實(shí)現(xiàn)自動迭代計(jì)算。當(dāng)?shù)螖?shù)達(dá)到設(shè)置次數(shù)或者目標(biāo)值誤差小于設(shè)置誤差時(shí)計(jì)算結(jié)束。當(dāng)計(jì)算結(jié)束時(shí),此時(shí)密度變量單元格的最新數(shù)值便是甲烷密度值。

    為簡化操作過程,將一系列操作過程錄制為宏并設(shè)置控件。若要獲得指定溫度和壓力下的甲烷密度,只需在輸入框中輸入溫度和壓力,計(jì)算機(jī)便會在自動運(yùn)行一系列程序后將甲烷密度的計(jì)算結(jié)果直接顯示在輸出框中。

    圖1 超臨界甲烷密度計(jì)算方法演示界面

    2.3 結(jié)果誤差分析

    為檢驗(yàn)本方法計(jì)算結(jié)果的準(zhǔn)確性,計(jì)算得到0~30 MPa,270、300、330、360 K四個(gè)溫度下的甲烷密度,以NIST商業(yè)軟件的計(jì)算結(jié)果為參考值進(jìn)行誤差分析,結(jié)果如圖2所示。相對誤差=(m–n)/n,其中,n為NIST軟件計(jì)算的甲烷密度,m為本文計(jì)算得出的甲烷密度。

    由圖2可知,各個(gè)溫度下隨著壓力的增大誤差都呈現(xiàn)上下波動的趨勢,誤差絕對值總體小于0.05%。在壓力為18 MPa之前,各個(gè)溫度下的誤差都較小,均小于0.012 5%;在18 MPa之后,各個(gè)溫度下的誤差都有所增大,其中,270 K溫度下誤差最大,但最大誤差都不超過0.05%。

    3 甲烷密度的變化規(guī)律

    為了分析不同溫度、壓力下甲烷的密度變化特征,采用本論文方法計(jì)算得出0~30 MPa,270、300、330和360 K 4個(gè)溫度下的甲烷密度(圖3中紅色虛線左側(cè)),同時(shí)利用NIST得到30~100 MPa下相應(yīng)溫度的甲烷密度(圖3中紅色虛線右側(cè)),結(jié)果如圖3所示。

    圖2 不同壓力、溫度下本文方法計(jì)算的密度值與NIST方法計(jì)算的密度值之間的誤差

    由甲烷密度隨壓力的變化曲線(圖3a)可以看出,隨著壓力的增大,各個(gè)溫度下的甲烷密度都呈現(xiàn)單調(diào)增加的趨勢,且隨壓力升高密度增加幅度越來越低。在臨界壓力(4.6 MPa)前后,密度沒有明顯的突變;在0~100 MPa、270~360 K范圍內(nèi)溫度對密度的影響遠(yuǎn)遠(yuǎn)小于壓力,隨著壓力的增大,不同溫度間的甲烷密度差異呈現(xiàn)先增大后減小的趨勢,在約30 MPa時(shí)密度差異最大而后開始減小(圖3a)。由圖3b可以看出,在低壓階段(約30 MPa之前)溫度對甲烷密度的影響較為顯著,隨壓力升高,溫度越低密度增加的速度越快,高壓階段(約30 MPa之后)不同溫度下密度隨壓力增加的幅度趨于相近。甲烷密度增量的變化規(guī)律以30 MPa為界大致分為兩個(gè)階段(圖3b):在30 MPa之前,密度增量隨壓力增大呈先增加后減小的拋物線形,且溫度越低密度增加幅度越大;在30 MPa之后,隨壓力升高甲烷密度增量繼續(xù)減小,溫度越低密度增量減小幅度越大。

    圖3 不同溫度、壓力下的甲烷密度和密度增量

    參考華北盆地的平均地溫梯度與正常儲層壓力梯度[21-24],得到埋深在100~3 000 m范圍內(nèi)的溫度和壓力分別為289.7~363.0 K、1~30 MPa,同時(shí)計(jì)算不同埋深對應(yīng)的溫度和壓力下的甲烷密度,如圖4所示。

    由圖4可知,隨著煤層埋藏深度逐漸增加,煤層壓力與溫度逐漸增大,甲烷密度在1 500 m以淺呈近線性增加,深度超過1 500 m以后增加幅度逐漸變緩。密度的變化同時(shí)受控于溫度與壓力,壓力增大使密度升高,溫度升高使密度降低,溫度與壓力的雙重作用導(dǎo)致密度隨埋深增加出現(xiàn)先呈近線性增加而后增幅變緩的變化規(guī)律。在埋藏深度增加的過程中,對密度的變化起主導(dǎo)作用的是壓力,但是隨埋深增加,壓力的影響逐漸變?nèi)酢?/p>

    由甲烷密度隨埋藏深度的變化規(guī)律可知,隨煤層埋深的增加,游離態(tài)甲烷的資源量逐漸增加。例如2 000 m埋深處甲烷的密度是500 m埋深的3.6倍,即假設(shè)煤層孔隙率不變,游離態(tài)資源量增加2.6倍。煤層對甲烷的吸附能力同樣受溫度、壓力雙重因素控制,埋藏深度超過1 000 m以后溫度對吸附量的影響超過壓力,從而出現(xiàn)吸附量降低的轉(zhuǎn)換帶[25]。所以煤層埋深超過1 000 m以后對資源量增量的貢獻(xiàn)主要來源于密度增加的游離態(tài)甲烷,這種效應(yīng)在以游離態(tài)為主的頁巖氣中尤為明顯。在埋藏深度超過2 000 m以后,隨深度繼續(xù)增加,游離態(tài)甲烷密度增加幅度變緩,對資源量增加的貢獻(xiàn)變?nèi)酢?/p>

    圖4 不同埋深下甲烷密度隨溫度、壓力的變化關(guān)系

    4 結(jié)論

    a. 基于亥姆霍茲能量基本方程的經(jīng)驗(yàn)狀態(tài)方程可以描述超臨界甲烷的熱力學(xué)行為,利用Microsoft Office Excel編寫的超臨界甲烷密度計(jì)算程序計(jì)算得出的270~360 K、0~30 MPa溫度與壓力范圍內(nèi)的甲烷密度,其準(zhǔn)確度符合要求。此方法程序準(zhǔn)確便捷,可以與煤的高壓甲烷等溫吸附實(shí)驗(yàn)數(shù)據(jù)計(jì)算程序相結(jié)合,更加準(zhǔn)確地計(jì)算儲層條件下游離態(tài)甲烷資源量和煤層中超臨界甲烷的吸附能力。

    b. 在270~360 K、0~100 MPa條件下,甲烷的密度隨著壓力的增大呈現(xiàn)單調(diào)增加的趨勢,密度隨壓力增加的速率在30 MPa時(shí)出現(xiàn)拐點(diǎn),高于30 MPa后不僅甲烷密度增加緩慢,同時(shí)密度對溫度的敏感性出現(xiàn)反轉(zhuǎn),在30 MPa前溫度越低密度隨壓力的增量越大,在高于30 MPa后溫度越低密度隨壓力的增量越小。

    c. 在平均地溫梯度與正常儲層壓力梯度下,隨煤層埋藏深度的增加,甲烷密度在約1 500 m以淺呈現(xiàn)近線性增加,而后隨埋深增加密度增加幅度逐漸變緩。與吸附態(tài)相比,游離態(tài)甲烷受溫度的影響較小,游離態(tài)甲烷資源量是深部煤層氣資源增量的主要貢獻(xiàn)者。

    [1] BANUTI D T. Crossing the widom-line-supercritical pseudo- boiling[J]. The Journal of Supercritical Fluids,2015,98:12–16.

    [2] 王艷,徐進(jìn)良,李文,等. 超臨界流體密度波動機(jī)理的分子動力學(xué)模擬[J]. 科學(xué)通報(bào),2020,65(17):1694–1704. WANG Yan,XU Jinliang,LI Wen,et al.Molecular dynamics study of mechanism of density fluctuation in supercritical fluid(in Chinese)[J]. Chinese Science Bulletin,2020,65(17):1694–1704.

    [3] 周理,李明,周亞平. 超臨界甲烷在高表面活性炭上的吸附測量及其理論分析[J]. 中國科學(xué):B輯,2000,30(1):49–56. ZHOU Li,LI Ming,ZHOU Yaping. Adsorption measurement and theoretical analysis of supercritical methane on high surface activated carbon[J]. Science in China(Series B),2000,30(1):49–56.

    [4] 梁洪彬,戚志林,向祖平,等. 超臨界頁巖氣密度特征研究[J]. 油氣藏評價(jià)與開發(fā),2017,7(5):74–79. LIANG Hongbin,QI Zhilin,XIANG Zuping,et al.Research on the characteristics of the density of supercritical shale gas[J]. Reservoir Evaluation and Development,2017,7(5):74–79.

    [5] ANGUS S,ARMSTRONG B,De REUCK K M. International thermodynamic tables of fluid state methane[M]. Oxford:International Union of Pure and Applied Chemistry,1978.

    [6] 童景山. 化工熱力學(xué)[M]. 北京:清華大學(xué)出版社,1995. TONG Jingshan. Chemical engineering thermodyna-miss[M]. Beijing:Tsinghua University Press,1995.

    [7] 楊兆彪,秦勇,高弟,等. 超臨界條件下煤層甲烷視吸附量、真實(shí)吸附量的差異及其地質(zhì)意義[J]. 天然氣工業(yè),2011,31(4):13–16. YANG Zhaobiao,QIN Yong,GAO Di,et al.Differences between apparent and true adsorption quantity of coalbed methane under supercritical conditions and their geological significance[J]. Natural Gas Industry,2011,31(4):13–16.

    [8] 孫紅明,傅雪海,丁永明,等. 煤層氣壓縮因子解析計(jì)算模型的建立[J]. 新疆石油地質(zhì),2014,35(2):215–220. SUN Hongming,F(xiàn)U Xuehai,DING Yongming,et al.Analytical calculation model for coalbed methane compressibility factor[J]. Xinjiang Petroleum Geology,2014,35(2):215–220.

    [9] 中國國家標(biāo)準(zhǔn)化管理委員會. 煤的高壓等溫吸附實(shí)驗(yàn)方法:GB/T 19560—2004[S]. 北京:中國標(biāo)準(zhǔn)出版社,2008. China National Standardization Management Committee. Experimental method of high-pressure isothermal adsorption to coal:GB/T 19560—2004[S]. Beijing:Standards Press of China,2008.

    [10] 時(shí)鈞,汪家鼎,余國琮,等. 化學(xué)工程手冊(第二版)上冊[M]. 北京:化學(xué)工業(yè)出版社,1996. SHI Jun,WANG Jiading,YU Guocong,et al. Chemical engineering handbook(2rdEdition)[M]. Beijing:Chemical Industry Press,1996.

    [11] ACHTERMANN H J,BOSE T K,ROGENER H,et al. Precise determination of the compressibility factor of methane,nitrogen,and their mixtures from refractive index measurements[J]. International Journal of Thermophysics,1986,7(3):709–720.

    [12] KLEINRAHM R,DUSCHEK W,WAGNER W,et al. Measurement and correlation of the(pressure,density,temperature) relation of methane in the temperature range from 273.15 K to 323.15 K at pressures up to 8 MPa[J]. The Journal of Chemical Thermodynamics,1988,20(5):621–631.

    [13] PIEPERBECK N,KLEINRAHM R,WAGNER W,et al. Results of(pressure,density,temperature) measurements on methane and on nitrogen in the temperature range from 273.15 K to 323.15 K at pressures up to 12 MPa using a new apparatus for accurate gas-density measurements[J]. The Journal of Chemical Thermodynamics,1991,23(2):175–194.

    [14] 熊杰,武雅珍,韓佳昊,等. 基于模糊控制方法的加氣站甲烷擴(kuò)散模型選擇[J]. 能源與環(huán)保,2019,41(7):57–62.XIONG Jie,WU Yazhen,HAN Jiahao,et al. Selection of methane diffusion model in gas station based on fuzzy control method[J]. China Energy and Environmental Protection,2019,41(7):57–62.

    [15] ACHTERMANN H J,HONG Jianguo,WAGNER W,et al. Refractive index and density isotherms for methane from 273 to 373 K and at pressures up to 34 MPa[J]. Journal of Chemical & Engineering Data,1992,37(4):414–418.

    [16] SETZMANN U,WAGNER W. A new method for optimizing the structure of thermodynamic correlation equations[J]. International Journal of Thermophysics,1989,10(6):1103–1126.

    [17] WAGNER W,SPAN R. Special equations of state for methane,argon,and nitrogen for the temperature range from 270 to 350 K at pressures up to 30 MPa[J]. International Journal of Thermophysics,1993,14(4):699–725.

    [18] SETZMANN U,WAGNER W. A new equation of state and tables of thermodynamic properties for methane covering the range from the melting line to 625 K at pressures up to 100 MPa[J]. Journal of Physical & Chemical Reference Data,1991,20(6):1061–1155.

    [19] 柳輝. 解非線性方程的牛頓迭代法及其應(yīng)用[J]. 重慶工學(xué)院學(xué)報(bào)(自然科學(xué)版),2007,21(8):95–98. LIU Hui.Newton iteration method for nonlinear equation solutions and its application[J]. Journal of Chongqing Institute of Technology(Natural Science Edition),2007,21(8):95–98.

    [20] 杜茂康,周玉敏,等. Excel與數(shù)據(jù)處理[M]. 北京:電子工業(yè)出版社,2009.DU Maokang,ZHOU Yumin,CAO Huiying,et al. Excel and data processing[M]. Beijing:Publishing House of Electronics Industry,2009.

    [21] 胡圣標(biāo),何麗娟,汪集旸. 中國大陸地區(qū)大地?zé)崃鲾?shù)據(jù)匯編(第三版)[J]. 地球物理學(xué)報(bào),2001,44(5):611–626. HU Shengbiao,HE Lijuan,WANG Jiyang. Compilation of heat flow data in the China continental area(3rdEdition)[J]. Chinese Journal of Geophysics,2001,44(5):611–626.

    [22] 王鈞,黃尚瑤,黃歌山,等. 中國地溫分布的基本特征[M]. 北京:地震出版社,1990. WANG Jun,HUANG Shangyao,HUANG Geshan,et al. The basic characteristics of China’s geothermal distribution[M]. Beijing:Seismological Press,1990.

    [23] 張鵬,王良書,劉紹文,等. 南華北盆地群地溫場研究[J]. 地球物理學(xué)進(jìn)展,2007,22(2):604–608. ZHANG Peng,WANG Liangshu,LIU Shaowen,et al.Geothermal field in the south Huabei basins[J]. Progress in Geophysics,2007,22(2):604–608.

    [24] 鐘玲文. 中國煤儲層壓力特征[J]. 天然氣工業(yè),2003,23(5):132–134. ZHONG Lingwen.Pressure characteristics of coal reservoirs in China[J]. Natural Gas Industry,2003,23(5):132–134.

    [25] 秦勇,申建,王寶文,等. 深部煤層氣成藏效應(yīng)及其耦合關(guān)系[J]. 石油學(xué)報(bào),2012,33(1):48–54. QIN Yong,SHEN Jian,WANG Baowen,et al. Accumulation effects and coupling relationship of deep coalbed methane[J]. Acta Petrolei Sinica,2012,33(1):48–54.

    Study on density of the supercritical methane

    ZHANG Xinbin1,2, SONG Dangyu1,2, LI Yunbo1,2, QI Yong’an1,2, FAN Shunli1,2

    (1. School of Resources and Environment, Henan Polytechnic University, Jiaozuo 454003, China; 2. Collaborative Innovation Center of Coalbed Methane and Shale Gas for Central Plains Economic Region, Jiaozuo 454003, China)

    The density of supercritical methane is not only an important parameter for calculation of unconventional natural gas resources and adsorption capacity of methane in coal, but also a main indicator for measurement of the diffusion efficiency of supercritical methane and the ability to dissolve organic matter. Through analysis and comparison of the applicability of various gas state equations, it is realized that the Helmholtz energy state equation can accurately calculate the density of methane under the conditions of 0-30 MPa and 270-360 K. A calculation program for methane density was written using Microsoft Office Excel, and the error was less than 0.05% compared with NIST commercial software. Finally, the changing of methane density was analyzed. In the range of 0-100 MPa and 270-360 K, the density of methane increases with the increase of pressure. The increase rate is greater before 30 MPa and the sensitivity to temperature is strong. Afterwards, the growth rate gradually slows down, while the sensitivity decreases. Under the in-situ condition of the coal seam, with the increase of buried depth, the density increases with temperature and decreases with pressure. Under the combined action of temperature and pressure, the density of methane first increases approximately linearly, and the growth rate gradually decreases.The results provide basic parameters for the occurrence and potential prediction of deep CBM.

    supercritical methane; methane density; gas state equation; Helmholtz energy; CBM resources

    移動閱讀

    語音講解

    P618.11

    A

    1001-1986(2021)01-0137-06

    2020-09-11;

    2020-12-14

    國家自然科學(xué)基金項(xiàng)目(U2003104);河南省自然科學(xué)基金項(xiàng)目(202300410186);河南理工大學(xué)青年骨干教師資助計(jì)劃項(xiàng)目(2018XQG-04)

    張新賓,1994年生,男,河南周口人,碩士,從事煤地質(zhì)及非常規(guī)天然氣地質(zhì)學(xué)方面的研究工作. E-mail:758320669@qq.com

    宋黨育,1971年生,男,河北獻(xiàn)縣人,博士,教授,從事煤層氣(瓦斯)地質(zhì)與煤地球化學(xué)方面的研究工作. E-mail:dangyusong@hpu.edu.cn

    張新賓,宋黨育,李云波,等.超臨界態(tài)甲烷密度研究[J]. 煤田地質(zhì)與勘探,2021,49(1):137–142. doi:10.3969/j.issn. 1001-1986.2021.01.014

    ZHANG Xinbin,SONG Dangyu,LI Yunbo,et al. Study on density of the supercritical methane[J]. Coal Geology & Exploration,2021,49(1):137–142. doi: 10.3969/j.issn.1001-1986.2021.01.014

    (責(zé)任編輯 范章群)

    猜你喜歡
    游離態(tài)狀態(tài)方程超臨界
    金屬及其化合物知識答疑
    超臨界CO2在頁巖氣開發(fā)中的應(yīng)用研究進(jìn)展
    云南化工(2021年5期)2021-12-21 07:41:20
    游離植物甾醇研發(fā)創(chuàng)新管理應(yīng)用
    LKP狀態(tài)方程在天然氣熱物性參數(shù)計(jì)算的應(yīng)用
    煤氣與熱力(2021年6期)2021-07-28 07:21:30
    農(nóng)村初中數(shù)學(xué)建構(gòu)式生態(tài)課堂中“游離態(tài)”學(xué)生成因及對策
    基于隨機(jī)與區(qū)間分析的狀態(tài)方程不確定性比較
    600MW超臨界機(jī)組熱經(jīng)濟(jì)性定量分析
    用狀態(tài)方程模擬氨基酸水溶液的熱力學(xué)性質(zhì)
    1200MW等級超超臨界機(jī)組可行性研究
    混合星物質(zhì)的狀態(tài)方程和奇異夸克物質(zhì)的穩(wěn)定窗
    久久精品久久久久久噜噜老黄| 欧美性感艳星| 蜜臀久久99精品久久宅男| 人人妻人人添人人爽欧美一区卜 | 高清av免费在线| 日韩不卡一区二区三区视频在线| 午夜福利在线在线| 在线天堂最新版资源| 黄色配什么色好看| 欧美+日韩+精品| 91精品一卡2卡3卡4卡| 日韩国内少妇激情av| 女人久久www免费人成看片| 久久热精品热| 亚洲怡红院男人天堂| 亚洲精品色激情综合| 精品一区在线观看国产| 天美传媒精品一区二区| 成人无遮挡网站| 亚洲国产精品成人久久小说| 韩国av在线不卡| 香蕉精品网在线| 免费黄频网站在线观看国产| 五月开心婷婷网| av国产免费在线观看| 国产综合精华液| 亚洲精品成人av观看孕妇| av专区在线播放| 大话2 男鬼变身卡| 国产亚洲午夜精品一区二区久久| 国产精品久久久久久精品电影小说 | 岛国毛片在线播放| 亚洲av成人精品一二三区| 久久99热这里只频精品6学生| 久久精品久久久久久久性| 欧美日本视频| 国产一区二区在线观看日韩| av卡一久久| 特大巨黑吊av在线直播| 国产探花极品一区二区| 亚洲欧美日韩无卡精品| 久久综合国产亚洲精品| 中文字幕亚洲精品专区| 99视频精品全部免费 在线| 成人国产麻豆网| 卡戴珊不雅视频在线播放| 极品教师在线视频| 黑丝袜美女国产一区| 少妇人妻 视频| 国产免费视频播放在线视频| av黄色大香蕉| 成人亚洲欧美一区二区av| 搡女人真爽免费视频火全软件| 久久 成人 亚洲| 纯流量卡能插随身wifi吗| 久久鲁丝午夜福利片| 亚洲美女黄色视频免费看| 亚洲国产精品一区三区| 亚洲成人一二三区av| 伦精品一区二区三区| 街头女战士在线观看网站| 舔av片在线| 青春草国产在线视频| 欧美日本视频| 精品久久久久久久久亚洲| 国产高清有码在线观看视频| 高清在线视频一区二区三区| 两个人的视频大全免费| 欧美另类一区| 国产精品蜜桃在线观看| 欧美区成人在线视频| 高清日韩中文字幕在线| 人妻少妇偷人精品九色| 亚洲,欧美,日韩| 亚洲熟女精品中文字幕| 午夜老司机福利剧场| 在线观看人妻少妇| 制服丝袜香蕉在线| 欧美少妇被猛烈插入视频| 韩国av在线不卡| 午夜免费男女啪啪视频观看| 日本av免费视频播放| 久久99热这里只频精品6学生| 久久久久久久大尺度免费视频| 久久久久国产网址| 亚洲一级一片aⅴ在线观看| 久久国内精品自在自线图片| 久久国产精品男人的天堂亚洲 | 人妻系列 视频| 亚洲欧美成人精品一区二区| 99热国产这里只有精品6| 六月丁香七月| 亚洲人成网站在线播| 成人亚洲欧美一区二区av| 一本色道久久久久久精品综合| 国产精品熟女久久久久浪| 大又大粗又爽又黄少妇毛片口| 久久久久久九九精品二区国产| 国产亚洲最大av| 国产视频内射| av福利片在线观看| 久久久精品94久久精品| 午夜精品国产一区二区电影| 国产白丝娇喘喷水9色精品| 精品一区二区免费观看| 亚洲欧美中文字幕日韩二区| 不卡视频在线观看欧美| 18禁裸乳无遮挡动漫免费视频| 五月天丁香电影| 久久毛片免费看一区二区三区| 在线观看免费日韩欧美大片 | 国产精品爽爽va在线观看网站| 国产精品国产三级国产专区5o| 伦精品一区二区三区| 久久精品国产自在天天线| av国产久精品久网站免费入址| 又黄又爽又刺激的免费视频.| av在线蜜桃| 久久久精品94久久精品| 亚洲精品乱码久久久v下载方式| 色网站视频免费| 亚洲精品色激情综合| 精品人妻熟女av久视频| 伊人久久精品亚洲午夜| 国产精品人妻久久久影院| 五月天丁香电影| 国产精品人妻久久久久久| 亚洲欧美精品专区久久| 亚洲国产最新在线播放| 日韩在线高清观看一区二区三区| 国产精品av视频在线免费观看| 插阴视频在线观看视频| 免费看不卡的av| 亚洲欧美日韩无卡精品| 国产精品嫩草影院av在线观看| 午夜免费观看性视频| 欧美日韩在线观看h| 国产免费视频播放在线视频| 国产欧美日韩精品一区二区| 在线观看一区二区三区| 丝袜脚勾引网站| 久久av网站| 亚洲av成人精品一二三区| 亚洲av免费高清在线观看| 人人妻人人看人人澡| 美女主播在线视频| 午夜老司机福利剧场| 日韩一区二区视频免费看| 亚洲欧美精品专区久久| av国产久精品久网站免费入址| 大又大粗又爽又黄少妇毛片口| 蜜臀久久99精品久久宅男| 国产免费视频播放在线视频| 欧美成人午夜免费资源| 久久久久国产网址| 婷婷色av中文字幕| 老女人水多毛片| 国产精品人妻久久久久久| 男的添女的下面高潮视频| 国产一区二区在线观看日韩| 国产精品女同一区二区软件| 简卡轻食公司| 免费大片黄手机在线观看| 七月丁香在线播放| 午夜激情福利司机影院| 欧美日韩一区二区视频在线观看视频在线| 在线观看国产h片| 久久久久久久大尺度免费视频| 国产精品久久久久久久电影| 综合色丁香网| 美女中出高潮动态图| 在线观看免费日韩欧美大片 | 一本—道久久a久久精品蜜桃钙片| 最近的中文字幕免费完整| 久久鲁丝午夜福利片| 日韩三级伦理在线观看| 欧美成人精品欧美一级黄| 日本-黄色视频高清免费观看| 亚洲精品国产成人久久av| 亚洲精品国产av蜜桃| 日本-黄色视频高清免费观看| 成人18禁高潮啪啪吃奶动态图 | 最后的刺客免费高清国语| 人妻 亚洲 视频| 成人黄色视频免费在线看| 国产免费视频播放在线视频| 免费看不卡的av| 草草在线视频免费看| 欧美精品亚洲一区二区| 九九在线视频观看精品| 观看美女的网站| 免费黄频网站在线观看国产| 一区二区三区精品91| 美女内射精品一级片tv| 亚洲久久久国产精品| 亚州av有码| 韩国高清视频一区二区三区| 久久人妻熟女aⅴ| 久久综合国产亚洲精品| av专区在线播放| 99热国产这里只有精品6| 中文乱码字字幕精品一区二区三区| 亚洲,欧美,日韩| 亚州av有码| 一级av片app| 精品熟女少妇av免费看| 天天躁日日操中文字幕| 色5月婷婷丁香| 日韩大片免费观看网站| 午夜精品国产一区二区电影| 国产精品国产三级专区第一集| 丰满少妇做爰视频| 国产一区二区三区综合在线观看 | 欧美人与善性xxx| 成人国产av品久久久| 国产精品久久久久久久电影| av在线播放精品| 青春草视频在线免费观看| 久久国产乱子免费精品| 少妇熟女欧美另类| 午夜福利网站1000一区二区三区| 六月丁香七月| 内射极品少妇av片p| xxx大片免费视频| 特大巨黑吊av在线直播| 自拍偷自拍亚洲精品老妇| 午夜福利影视在线免费观看| 国产成人精品福利久久| 青春草亚洲视频在线观看| 精品99又大又爽又粗少妇毛片| 午夜福利网站1000一区二区三区| 91精品伊人久久大香线蕉| 小蜜桃在线观看免费完整版高清| 国产精品人妻久久久影院| 国产精品久久久久久久久免| 国产黄色免费在线视频| 麻豆成人午夜福利视频| 熟女电影av网| 少妇猛男粗大的猛烈进出视频| av视频免费观看在线观看| 亚洲av男天堂| 99热这里只有精品一区| 街头女战士在线观看网站| 男人和女人高潮做爰伦理| 亚洲欧美日韩东京热| 亚洲欧洲日产国产| 夜夜爽夜夜爽视频| 国内揄拍国产精品人妻在线| 久久精品久久久久久久性| 亚洲av日韩在线播放| 国产免费视频播放在线视频| 十八禁网站网址无遮挡 | 国产 精品1| 国产高清不卡午夜福利| 综合色丁香网| 久久亚洲国产成人精品v| 五月玫瑰六月丁香| videos熟女内射| 久久国产乱子免费精品| 人妻少妇偷人精品九色| 中文资源天堂在线| 亚洲熟女精品中文字幕| 亚洲,一卡二卡三卡| 欧美3d第一页| 国产精品不卡视频一区二区| 永久免费av网站大全| 亚洲精品色激情综合| 国产淫片久久久久久久久| 少妇的逼水好多| 我的老师免费观看完整版| 在线播放无遮挡| 成人毛片a级毛片在线播放| 亚洲精品国产成人久久av| 建设人人有责人人尽责人人享有的 | 国产中年淑女户外野战色| 国产av一区二区精品久久 | 丰满迷人的少妇在线观看| 嘟嘟电影网在线观看| 熟女av电影| 婷婷色综合大香蕉| 午夜福利视频精品| 中文在线观看免费www的网站| 欧美三级亚洲精品| 亚洲自偷自拍三级| 欧美xxxx性猛交bbbb| 色5月婷婷丁香| 精品国产三级普通话版| 国产中年淑女户外野战色| 寂寞人妻少妇视频99o| av国产久精品久网站免费入址| 日韩国内少妇激情av| 日韩强制内射视频| 赤兔流量卡办理| 亚洲自偷自拍三级| 99热这里只有是精品50| 七月丁香在线播放| 欧美成人午夜免费资源| 视频区图区小说| 超碰97精品在线观看| 久久久亚洲精品成人影院| 亚洲精品国产av成人精品| 性高湖久久久久久久久免费观看| 亚洲,一卡二卡三卡| 91狼人影院| 欧美少妇被猛烈插入视频| 国产成人a区在线观看| 亚洲国产精品专区欧美| 国产欧美日韩一区二区三区在线 | 精品人妻偷拍中文字幕| 精品人妻熟女av久视频| 在线观看av片永久免费下载| 久久久久网色| 中文天堂在线官网| 精品人妻熟女av久视频| 国产乱人视频| 国产成人aa在线观看| 永久网站在线| 成年人午夜在线观看视频| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美清纯卡通| 亚洲电影在线观看av| 99久久综合免费| 久久精品国产亚洲av涩爱| 国产精品爽爽va在线观看网站| 精品视频人人做人人爽| 精品午夜福利在线看| 内地一区二区视频在线| 成年美女黄网站色视频大全免费 | 国产av码专区亚洲av| 成人特级av手机在线观看| 91精品伊人久久大香线蕉| 久久久欧美国产精品| 日韩制服骚丝袜av| 亚洲美女视频黄频| 狂野欧美白嫩少妇大欣赏| 亚洲中文av在线| 一级片'在线观看视频| 欧美 日韩 精品 国产| 身体一侧抽搐| 一区二区三区四区激情视频| 青春草视频在线免费观看| 久久精品人妻少妇| 国产综合精华液| 高清在线视频一区二区三区| 欧美最新免费一区二区三区| 亚洲精品乱久久久久久| 成年人午夜在线观看视频| 91狼人影院| 草草在线视频免费看| 成年人午夜在线观看视频| 国产精品99久久久久久久久| 日本vs欧美在线观看视频 | a 毛片基地| 国产成人精品福利久久| 国产欧美另类精品又又久久亚洲欧美| 超碰97精品在线观看| 自拍偷自拍亚洲精品老妇| 三级国产精品片| 亚洲精品视频女| 最近的中文字幕免费完整| 精品人妻偷拍中文字幕| 亚洲欧美日韩另类电影网站 | 日本欧美国产在线视频| 高清黄色对白视频在线免费看 | 精品一品国产午夜福利视频| 黑人猛操日本美女一级片| 亚洲国产色片| 国产女主播在线喷水免费视频网站| 国内精品宾馆在线| 交换朋友夫妻互换小说| 少妇 在线观看| 亚洲精品视频女| 99热全是精品| 亚洲欧美日韩无卡精品| 亚洲精品国产av蜜桃| 最新中文字幕久久久久| 久久精品国产鲁丝片午夜精品| 亚洲美女视频黄频| 国产色爽女视频免费观看| 日韩三级伦理在线观看| 在线观看免费视频网站a站| 亚洲欧洲日产国产| 欧美老熟妇乱子伦牲交| 久久久久久久精品精品| h日本视频在线播放| av卡一久久| 国产男女超爽视频在线观看| 高清欧美精品videossex| 免费观看在线日韩| 水蜜桃什么品种好| 亚洲欧美精品专区久久| 午夜精品国产一区二区电影| 亚洲欧美成人综合另类久久久| 国产深夜福利视频在线观看| 我的女老师完整版在线观看| 联通29元200g的流量卡| 伦理电影免费视频| 亚洲在久久综合| 亚洲精品久久午夜乱码| 亚洲经典国产精华液单| 国产精品久久久久成人av| 亚洲第一区二区三区不卡| 一本一本综合久久| 亚洲av国产av综合av卡| 91精品国产九色| 最新中文字幕久久久久| 国产 一区 欧美 日韩| 大码成人一级视频| 久久久精品免费免费高清| 国产精品一区二区在线不卡| 久久久久国产网址| 91精品国产国语对白视频| 黄色日韩在线| 2021少妇久久久久久久久久久| 1000部很黄的大片| 人妻系列 视频| 久久国产乱子免费精品| 久久久久久久久久久丰满| 最近最新中文字幕免费大全7| 亚洲第一av免费看| 久久久精品94久久精品| 在线精品无人区一区二区三 | 欧美xxⅹ黑人| 国产美女午夜福利| 国产一区二区三区av在线| 亚洲精品国产av蜜桃| 亚洲综合色惰| 亚洲丝袜综合中文字幕| 欧美一区二区亚洲| 亚洲av二区三区四区| 夜夜骑夜夜射夜夜干| 国产高潮美女av| 51国产日韩欧美| 国产精品成人在线| 国产爽快片一区二区三区| 高清日韩中文字幕在线| 99久久精品一区二区三区| 在线观看国产h片| 成人高潮视频无遮挡免费网站| 卡戴珊不雅视频在线播放| 国产av码专区亚洲av| 在线免费十八禁| 国产亚洲5aaaaa淫片| 极品教师在线视频| 色吧在线观看| 亚洲人成网站在线观看播放| 天美传媒精品一区二区| 国产高清不卡午夜福利| 国产精品一区www在线观看| 免费高清在线观看视频在线观看| 亚洲国产成人一精品久久久| 最近最新中文字幕免费大全7| 久久婷婷青草| 人人妻人人添人人爽欧美一区卜 | 日韩成人av中文字幕在线观看| 三级经典国产精品| 久久久久网色| 国产精品无大码| 亚洲精品一区蜜桃| 一边亲一边摸免费视频| 日日摸夜夜添夜夜添av毛片| av卡一久久| 男女免费视频国产| 性色av一级| 草草在线视频免费看| 国产免费又黄又爽又色| 美女xxoo啪啪120秒动态图| 成人免费观看视频高清| 五月天丁香电影| 一区二区三区精品91| 熟妇人妻不卡中文字幕| 大码成人一级视频| 国产成人一区二区在线| 欧美日韩国产mv在线观看视频 | 美女视频免费永久观看网站| av免费在线看不卡| 久久毛片免费看一区二区三区| av国产精品久久久久影院| 国产熟女欧美一区二区| 少妇熟女欧美另类| 国产伦精品一区二区三区视频9| 国产永久视频网站| 高清不卡的av网站| 全区人妻精品视频| 在线免费观看不下载黄p国产| 少妇人妻一区二区三区视频| 精品亚洲成国产av| 中国国产av一级| 黄色日韩在线| 国产黄色视频一区二区在线观看| 久久精品国产自在天天线| 国产一区有黄有色的免费视频| 一级毛片电影观看| 成人免费观看视频高清| 人妻少妇偷人精品九色| 3wmmmm亚洲av在线观看| 精品一区二区三卡| 久久青草综合色| 永久网站在线| 亚洲精品中文字幕在线视频 | 成人特级av手机在线观看| 九草在线视频观看| 国产精品熟女久久久久浪| 亚洲精品,欧美精品| 国产乱人偷精品视频| 男女边吃奶边做爰视频| 国产在线免费精品| 国产老妇伦熟女老妇高清| 一区二区三区精品91| 高清在线视频一区二区三区| 国产亚洲最大av| 亚洲经典国产精华液单| 狂野欧美激情性bbbbbb| 亚洲av成人精品一二三区| 大又大粗又爽又黄少妇毛片口| 少妇裸体淫交视频免费看高清| 网址你懂的国产日韩在线| 亚洲国产色片| 久久精品国产自在天天线| 欧美97在线视频| 黄色怎么调成土黄色| 成年美女黄网站色视频大全免费 | 精品人妻熟女av久视频| 99热国产这里只有精品6| 美女中出高潮动态图| 啦啦啦中文免费视频观看日本| 亚洲第一区二区三区不卡| 性色av一级| 97超视频在线观看视频| 91狼人影院| 亚洲欧美清纯卡通| 国产一区二区三区av在线| 欧美日韩亚洲高清精品| 精品国产三级普通话版| 最近最新中文字幕大全电影3| 啦啦啦在线观看免费高清www| 亚洲av中文字字幕乱码综合| 欧美另类一区| 国产淫语在线视频| 久久综合国产亚洲精品| 九草在线视频观看| 精品久久久精品久久久| 国产免费一级a男人的天堂| 亚洲欧美一区二区三区国产| 只有这里有精品99| 美女高潮的动态| 91aial.com中文字幕在线观看| 在线亚洲精品国产二区图片欧美 | 边亲边吃奶的免费视频| 国产精品伦人一区二区| 自拍欧美九色日韩亚洲蝌蚪91 | 视频区图区小说| 亚洲第一av免费看| 久久97久久精品| 22中文网久久字幕| 精品亚洲乱码少妇综合久久| 舔av片在线| 免费观看在线日韩| 欧美 日韩 精品 国产| 久久久国产一区二区| 少妇高潮的动态图| 亚洲精品第二区| 久久久成人免费电影| 在线播放无遮挡| 日本一二三区视频观看| 五月伊人婷婷丁香| 搡老乐熟女国产| 国产有黄有色有爽视频| 婷婷色av中文字幕| 亚洲内射少妇av| 波野结衣二区三区在线| 建设人人有责人人尽责人人享有的 | 亚洲,一卡二卡三卡| 国产永久视频网站| 极品少妇高潮喷水抽搐| 综合色丁香网| 久久久久久九九精品二区国产| 老熟女久久久| 久久综合国产亚洲精品| 亚洲av日韩在线播放| 免费黄网站久久成人精品| 久久ye,这里只有精品| 精品熟女少妇av免费看| 欧美极品一区二区三区四区| 激情五月婷婷亚洲| 欧美日韩一区二区视频在线观看视频在线| 91精品国产九色| 免费av中文字幕在线| 熟女电影av网| 久久久久久久大尺度免费视频| 欧美少妇被猛烈插入视频| 人体艺术视频欧美日本| 一级av片app| 夫妻午夜视频| 一级毛片黄色毛片免费观看视频| 亚洲av综合色区一区| 色哟哟·www| 久久久久久久久久人人人人人人| 三级国产精品片| 午夜激情久久久久久久| 人人妻人人澡人人爽人人夜夜| 国产视频内射| 日韩av免费高清视频| 大陆偷拍与自拍| 晚上一个人看的免费电影| 亚洲精品国产成人久久av| 日本-黄色视频高清免费观看| 欧美日韩亚洲高清精品| 国产视频首页在线观看| 日本黄大片高清| 久久精品国产亚洲av天美| 美女xxoo啪啪120秒动态图| 国产高潮美女av| 国产大屁股一区二区在线视频| 国产高清有码在线观看视频|