• <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)定窗
    免费在线观看视频国产中文字幕亚洲| av超薄肉色丝袜交足视频| 制服诱惑二区| av网站在线播放免费| 99国产精品99久久久久| 久久精品国产a三级三级三级| 日韩欧美一区二区三区在线观看 | 十分钟在线观看高清视频www| 50天的宝宝边吃奶边哭怎么回事| 亚洲三区欧美一区| 美女高潮喷水抽搐中文字幕| 69av精品久久久久久| 午夜老司机福利片| 在线观看舔阴道视频| 欧美+亚洲+日韩+国产| 国产一区二区三区综合在线观看| 国产欧美日韩综合在线一区二区| 久久青草综合色| 老司机在亚洲福利影院| 欧美午夜高清在线| 十八禁人妻一区二区| 校园春色视频在线观看| 在线十欧美十亚洲十日本专区| av天堂在线播放| 岛国毛片在线播放| 亚洲人成伊人成综合网2020| 9191精品国产免费久久| 婷婷精品国产亚洲av在线 | 精品久久久精品久久久| 色在线成人网| 男女下面插进去视频免费观看| 国产麻豆69| 在线观看免费午夜福利视频| 男人舔女人的私密视频| 日韩中文字幕欧美一区二区| 久久久久久免费高清国产稀缺| 大型黄色视频在线免费观看| 国产蜜桃级精品一区二区三区 | 中文字幕另类日韩欧美亚洲嫩草| 中文字幕人妻丝袜一区二区| 久久国产精品男人的天堂亚洲| 美女国产高潮福利片在线看| 一本大道久久a久久精品| 中文字幕精品免费在线观看视频| 欧美日韩精品网址| 精品电影一区二区在线| 一边摸一边做爽爽视频免费| 日韩制服丝袜自拍偷拍| 俄罗斯特黄特色一大片| 成熟少妇高潮喷水视频| 女性被躁到高潮视频| 精品久久久久久,| 国产有黄有色有爽视频| 国产亚洲欧美精品永久| 国产亚洲精品久久久久久毛片 | 怎么达到女性高潮| 日韩欧美国产一区二区入口| 精品免费久久久久久久清纯 | 亚洲av美国av| 亚洲人成电影免费在线| 无限看片的www在线观看| 18禁观看日本| 黄频高清免费视频| 欧美一级毛片孕妇| 女同久久另类99精品国产91| 日韩三级视频一区二区三区| 日韩一卡2卡3卡4卡2021年| 亚洲国产欧美网| 人人妻,人人澡人人爽秒播| 午夜福利免费观看在线| 精品久久蜜臀av无| 一二三四在线观看免费中文在| 国产av一区二区精品久久| 一级,二级,三级黄色视频| 精品欧美一区二区三区在线| 欧美日韩av久久| 一区二区日韩欧美中文字幕| 国产精品成人在线| 久久国产亚洲av麻豆专区| 欧美+亚洲+日韩+国产| 午夜福利在线免费观看网站| 老司机午夜福利在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久久人人人人人| 国产一区二区三区在线臀色熟女 | 9191精品国产免费久久| 国产精品香港三级国产av潘金莲| 国产精品久久久久久人妻精品电影| 美女扒开内裤让男人捅视频| 国产亚洲一区二区精品| 国产一区二区三区在线臀色熟女 | 亚洲va日本ⅴa欧美va伊人久久| 免费少妇av软件| 欧美精品人与动牲交sv欧美| 人人澡人人妻人| 丝袜在线中文字幕| 中文字幕另类日韩欧美亚洲嫩草| 成人影院久久| 伊人久久大香线蕉亚洲五| 露出奶头的视频| 亚洲免费av在线视频| 又黄又粗又硬又大视频| 丰满的人妻完整版| 午夜精品久久久久久毛片777| 啦啦啦在线免费观看视频4| 丝瓜视频免费看黄片| 成在线人永久免费视频| 国产精品欧美亚洲77777| 又紧又爽又黄一区二区| 天堂√8在线中文| 女警被强在线播放| 欧美性长视频在线观看| 久久精品国产亚洲av香蕉五月 | 乱人伦中国视频| 精品国产乱码久久久久久男人| 欧美成人免费av一区二区三区 | 免费女性裸体啪啪无遮挡网站| 欧美一级毛片孕妇| 久热爱精品视频在线9| 欧美黄色片欧美黄色片| 久久久精品国产亚洲av高清涩受| 一边摸一边抽搐一进一出视频| 午夜福利一区二区在线看| 亚洲五月婷婷丁香| 天天躁日日躁夜夜躁夜夜| 无遮挡黄片免费观看| 亚洲黑人精品在线| 亚洲精品国产一区二区精华液| 高清在线国产一区| 少妇裸体淫交视频免费看高清 | 成人18禁高潮啪啪吃奶动态图| 在线观看免费视频日本深夜| 亚洲av日韩精品久久久久久密| 首页视频小说图片口味搜索| 国产精品1区2区在线观看. | 精品久久久久久久久久免费视频 | 国产伦人伦偷精品视频| av一本久久久久| 亚洲精品成人av观看孕妇| 精品电影一区二区在线| 精品乱码久久久久久99久播| 男女免费视频国产| 欧美人与性动交α欧美精品济南到| 国产精品亚洲av一区麻豆| 欧美久久黑人一区二区| 久久久久久免费高清国产稀缺| 视频区欧美日本亚洲| 成人国语在线视频| 国产精品98久久久久久宅男小说| 成人亚洲精品一区在线观看| 成人影院久久| 色在线成人网| 亚洲熟女精品中文字幕| 美国免费a级毛片| 亚洲精品久久成人aⅴ小说| 国产色视频综合| 热re99久久国产66热| 女人高潮潮喷娇喘18禁视频| 久久中文看片网| 亚洲一区中文字幕在线| 免费在线观看影片大全网站| 巨乳人妻的诱惑在线观看| 久久久久精品人妻al黑| 亚洲国产欧美一区二区综合| 99久久国产精品久久久| 一区二区三区国产精品乱码| 一边摸一边做爽爽视频免费| 搡老岳熟女国产| 黄色视频不卡| 超色免费av| 免费av中文字幕在线| 在线十欧美十亚洲十日本专区| 精品亚洲成国产av| 91成年电影在线观看| 久9热在线精品视频| 两个人看的免费小视频| 久久国产精品影院| 一区福利在线观看| 亚洲七黄色美女视频| 精品福利永久在线观看| 久久这里只有精品19| 少妇猛男粗大的猛烈进出视频| 国产无遮挡羞羞视频在线观看| 一个人免费在线观看的高清视频| 欧美日韩亚洲国产一区二区在线观看 | 欧美亚洲日本最大视频资源| 国产成人精品久久二区二区91| 黄色a级毛片大全视频| 免费看a级黄色片| 国产精品久久视频播放| 国产亚洲精品一区二区www | 无限看片的www在线观看| 亚洲精品一二三| 视频区图区小说| 亚洲午夜精品一区,二区,三区| 男女床上黄色一级片免费看| 欧美丝袜亚洲另类 | 欧美日韩精品网址| 国产精品自产拍在线观看55亚洲 | 国产精品.久久久| 久久久久视频综合| 精品一区二区三区视频在线观看免费 | 日韩欧美三级三区| 亚洲欧美激情综合另类| 美女视频免费永久观看网站| a级毛片在线看网站| 亚洲 欧美一区二区三区| 十八禁高潮呻吟视频| 最近最新免费中文字幕在线| 美国免费a级毛片| 精品乱码久久久久久99久播| 欧美丝袜亚洲另类 | 黑人巨大精品欧美一区二区mp4| 90打野战视频偷拍视频| 正在播放国产对白刺激| 美女午夜性视频免费| 欧美色视频一区免费| 国产精品久久久久成人av| 高清视频免费观看一区二区| 成年女人毛片免费观看观看9 | 色精品久久人妻99蜜桃| 两性午夜刺激爽爽歪歪视频在线观看 | 99国产综合亚洲精品| 成人影院久久| videos熟女内射| 91国产中文字幕| 欧美乱色亚洲激情| 亚洲av成人av| 久久影院123| 久久国产乱子伦精品免费另类| 青草久久国产| 一级a爱视频在线免费观看| 国产免费现黄频在线看| 欧美精品av麻豆av| 国产又色又爽无遮挡免费看| 欧美黄色淫秽网站| av线在线观看网站| 午夜精品久久久久久毛片777| 亚洲午夜精品一区,二区,三区| 国产一区二区激情短视频| 黑人猛操日本美女一级片| 日本五十路高清| 国产成人影院久久av| 一本大道久久a久久精品| 天堂√8在线中文| 欧美中文综合在线视频| av中文乱码字幕在线| 青草久久国产| 亚洲精品成人av观看孕妇| 精品国产乱码久久久久久男人| 人人妻人人添人人爽欧美一区卜| 亚洲第一av免费看| 精品久久久久久电影网| 国产欧美日韩一区二区三区在线| 国产精品偷伦视频观看了| 超碰成人久久| av有码第一页| 我的亚洲天堂| 午夜福利,免费看| 大码成人一级视频| videos熟女内射| 亚洲一卡2卡3卡4卡5卡精品中文| 精品一区二区三区视频在线观看免费 | 久久久久久久午夜电影 | 久久久久久免费高清国产稀缺| 在线看a的网站| 啦啦啦视频在线资源免费观看| 日韩欧美一区二区三区在线观看 | 成人国产一区最新在线观看| 国产欧美日韩一区二区精品| 欧美日韩精品网址| 日本a在线网址| 丝袜美足系列| 国产精品九九99| 精品第一国产精品| 亚洲精品一二三| 成熟少妇高潮喷水视频| 丁香欧美五月| 国产国语露脸激情在线看| 黄网站色视频无遮挡免费观看| 18禁美女被吸乳视频| 91精品国产国语对白视频| 黄片小视频在线播放| 国产激情久久老熟女| 乱人伦中国视频| 男人操女人黄网站| 国产精品一区二区精品视频观看| 亚洲av第一区精品v没综合| 啦啦啦视频在线资源免费观看| 国产91精品成人一区二区三区| 91麻豆av在线| 国产精品久久久av美女十八| bbb黄色大片| 999久久久精品免费观看国产| 日韩人妻精品一区2区三区| 老司机深夜福利视频在线观看| 视频区欧美日本亚洲| 91老司机精品| 日韩欧美三级三区| 精品第一国产精品| 欧美激情 高清一区二区三区| 免费久久久久久久精品成人欧美视频| 午夜福利欧美成人| 丝袜美腿诱惑在线| 久久精品国产综合久久久| 一二三四社区在线视频社区8| svipshipincom国产片| 人人妻,人人澡人人爽秒播| 色精品久久人妻99蜜桃| 中文字幕av电影在线播放| 欧美+亚洲+日韩+国产| 老汉色av国产亚洲站长工具| 可以免费在线观看a视频的电影网站| 中文亚洲av片在线观看爽 | 无限看片的www在线观看| 国产精品免费一区二区三区在线 | 精品久久久精品久久久| 精品国产一区二区三区四区第35| 视频在线观看一区二区三区| 久久国产精品人妻蜜桃| 大型黄色视频在线免费观看| 亚洲七黄色美女视频| 久久国产精品人妻蜜桃| 日韩欧美国产一区二区入口| 交换朋友夫妻互换小说| 欧美久久黑人一区二区| 老熟女久久久| 一级,二级,三级黄色视频| 一级毛片精品| 国产真人三级小视频在线观看| 久久精品成人免费网站| 人妻丰满熟妇av一区二区三区 | 亚洲一码二码三码区别大吗| 高潮久久久久久久久久久不卡| 久久精品国产a三级三级三级| 中文字幕人妻熟女乱码| 亚洲全国av大片| 国产熟女午夜一区二区三区| 国产精品综合久久久久久久免费 | 亚洲中文av在线| 黄色a级毛片大全视频| 亚洲,欧美精品.| 国产高清videossex| 精品视频人人做人人爽| svipshipincom国产片| 国产主播在线观看一区二区| √禁漫天堂资源中文www| 久久国产精品男人的天堂亚洲| 成人国语在线视频| 日本欧美视频一区| www.精华液| 在线观看免费视频网站a站| 国产亚洲欧美精品永久| 久久精品国产亚洲av香蕉五月 | 99在线人妻在线中文字幕 | 亚洲欧美一区二区三区黑人| 美女国产高潮福利片在线看| a级毛片在线看网站| 欧美日韩精品网址| 18在线观看网站| 在线观看免费视频日本深夜| 欧美精品一区二区免费开放| 久久婷婷成人综合色麻豆| 久久午夜亚洲精品久久| 一进一出抽搐gif免费好疼 | 一本综合久久免费| 一夜夜www| 欧美 亚洲 国产 日韩一| 国产亚洲精品久久久久5区| 日韩免费高清中文字幕av| 叶爱在线成人免费视频播放| 欧美国产精品va在线观看不卡| 亚洲av第一区精品v没综合| 免费在线观看视频国产中文字幕亚洲| 69av精品久久久久久| a级毛片黄视频| 亚洲av欧美aⅴ国产| 水蜜桃什么品种好| 大香蕉久久成人网| 亚洲国产欧美日韩在线播放| 国产成人一区二区三区免费视频网站| 男女午夜视频在线观看| 大码成人一级视频| av视频免费观看在线观看| 丰满人妻熟妇乱又伦精品不卡| 精品国产美女av久久久久小说| 欧美日韩亚洲国产一区二区在线观看 | 一夜夜www| 757午夜福利合集在线观看| 99国产精品一区二区蜜桃av | 欧美人与性动交α欧美精品济南到| 丝袜在线中文字幕| x7x7x7水蜜桃| 搡老岳熟女国产| 欧美另类亚洲清纯唯美| 宅男免费午夜| 黄色怎么调成土黄色| 捣出白浆h1v1| 国产成人精品久久二区二区免费| 成在线人永久免费视频| 亚洲国产精品合色在线| 国产黄色免费在线视频| 亚洲国产欧美日韩在线播放| 老司机午夜十八禁免费视频| 巨乳人妻的诱惑在线观看| 国产精品免费大片| 日韩欧美三级三区| 久久 成人 亚洲| 视频区欧美日本亚洲| 91九色精品人成在线观看| 午夜亚洲福利在线播放| 欧美不卡视频在线免费观看 | 丝瓜视频免费看黄片| 亚洲一区二区三区欧美精品| 国产免费男女视频| 久久久久精品人妻al黑| 两个人看的免费小视频| 他把我摸到了高潮在线观看| 久久久久久亚洲精品国产蜜桃av| 又紧又爽又黄一区二区| 国产av精品麻豆| 怎么达到女性高潮| 五月开心婷婷网| 91九色精品人成在线观看| 精品免费久久久久久久清纯 | 制服诱惑二区| av线在线观看网站| av网站在线播放免费| 他把我摸到了高潮在线观看| 国产99久久九九免费精品| 亚洲中文av在线| 美国免费a级毛片| 久久久久久久精品吃奶| 亚洲精品成人av观看孕妇| 午夜精品在线福利| 老熟女久久久| 久久久国产一区二区| 日韩欧美一区二区三区在线观看 | 午夜福利一区二区在线看| xxx96com| 露出奶头的视频| 欧美性长视频在线观看| 日韩大码丰满熟妇| 欧美日韩亚洲高清精品| 久久久国产成人精品二区 | av超薄肉色丝袜交足视频| 黄网站色视频无遮挡免费观看| 99精品久久久久人妻精品| 乱人伦中国视频| 777米奇影视久久| 在线永久观看黄色视频| 亚洲 欧美一区二区三区| 黄色视频,在线免费观看| www日本在线高清视频| 亚洲五月天丁香| 麻豆乱淫一区二区| 三级毛片av免费| 亚洲国产精品sss在线观看 | 在线国产一区二区在线| 亚洲一码二码三码区别大吗| 在线观看www视频免费| 欧美日韩福利视频一区二区| 在线观看一区二区三区激情| 久久影院123| 国产亚洲欧美在线一区二区| 老司机亚洲免费影院| 电影成人av| 欧美日韩黄片免| 老司机午夜十八禁免费视频| 在线观看免费视频网站a站| 91老司机精品| 国产有黄有色有爽视频| 久久国产精品影院| svipshipincom国产片| 亚洲成av片中文字幕在线观看| 在线国产一区二区在线| 啪啪无遮挡十八禁网站| 欧美 日韩 精品 国产| 欧美丝袜亚洲另类 | av欧美777| 9色porny在线观看| 99精国产麻豆久久婷婷| 一二三四在线观看免费中文在| 国产熟女午夜一区二区三区| 亚洲欧美色中文字幕在线| 搡老乐熟女国产| 五月开心婷婷网| www日本在线高清视频| 国产精华一区二区三区| 国产日韩欧美亚洲二区| 中文字幕最新亚洲高清| 一个人免费在线观看的高清视频| 国产乱人伦免费视频| 两性夫妻黄色片| 乱人伦中国视频| 午夜久久久在线观看| 国产欧美日韩精品亚洲av| 女人爽到高潮嗷嗷叫在线视频| 亚洲黑人精品在线| 久久久久久人人人人人| 大片电影免费在线观看免费| 老汉色av国产亚洲站长工具| 在线十欧美十亚洲十日本专区| 亚洲五月婷婷丁香| 国产亚洲精品第一综合不卡| 正在播放国产对白刺激| 午夜福利在线免费观看网站| 亚洲欧美精品综合一区二区三区| 人人妻人人澡人人看| av在线播放免费不卡| 亚洲欧美日韩高清在线视频| 天堂中文最新版在线下载| av不卡在线播放| 亚洲av熟女| 久久国产精品影院| x7x7x7水蜜桃| 国产真人三级小视频在线观看| 亚洲综合色网址| 午夜福利乱码中文字幕| 美女福利国产在线| 99国产精品99久久久久| 大型av网站在线播放| 丝袜在线中文字幕| 成熟少妇高潮喷水视频| 国产欧美亚洲国产| 激情在线观看视频在线高清 | 免费少妇av软件| 人人妻人人添人人爽欧美一区卜| 一夜夜www| 亚洲精品中文字幕在线视频| 亚洲成国产人片在线观看| 欧美老熟妇乱子伦牲交| 亚洲专区字幕在线| 后天国语完整版免费观看| 一级a爱片免费观看的视频| 人人澡人人妻人| 性色av乱码一区二区三区2| 日韩中文字幕欧美一区二区| 成人国语在线视频| 99re6热这里在线精品视频| 两性午夜刺激爽爽歪歪视频在线观看 | 久久这里只有精品19| 国产蜜桃级精品一区二区三区 | a级毛片在线看网站| 黑人巨大精品欧美一区二区蜜桃| 多毛熟女@视频| 国产精品一区二区在线观看99| 97人妻天天添夜夜摸| 日韩有码中文字幕| 久久精品熟女亚洲av麻豆精品| 看免费av毛片| 久久亚洲真实| 超色免费av| 成年版毛片免费区| 国产成人av教育| 精品国产美女av久久久久小说| 亚洲午夜理论影院| 在线观看免费日韩欧美大片| 丁香欧美五月| 亚洲精品av麻豆狂野| av超薄肉色丝袜交足视频| 乱人伦中国视频| 精品第一国产精品| 亚洲av日韩在线播放| 丝袜在线中文字幕| √禁漫天堂资源中文www| 人人妻人人澡人人爽人人夜夜| 麻豆成人av在线观看| 夜夜夜夜夜久久久久| 高清视频免费观看一区二区| 久久久久国产精品人妻aⅴ院 | 免费看a级黄色片| 黄片播放在线免费| 国产高清视频在线播放一区| 国产主播在线观看一区二区| 色尼玛亚洲综合影院| 欧美激情 高清一区二区三区| 午夜免费鲁丝| av有码第一页| 国产高清国产精品国产三级| 91大片在线观看| 国产精品免费一区二区三区在线 | 看黄色毛片网站| 欧美日韩亚洲综合一区二区三区_| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲av成人av| 欧美日韩精品网址| 两人在一起打扑克的视频| 日本精品一区二区三区蜜桃| a级毛片黄视频| aaaaa片日本免费| 亚洲熟妇中文字幕五十中出 | 国产精品久久久人人做人人爽| 久久人人爽av亚洲精品天堂| 国产欧美日韩精品亚洲av| 日日摸夜夜添夜夜添小说| 美女高潮到喷水免费观看| 9191精品国产免费久久| 狠狠婷婷综合久久久久久88av| 欧美性长视频在线观看| 中亚洲国语对白在线视频| 妹子高潮喷水视频| 免费在线观看亚洲国产| 人妻 亚洲 视频| 亚洲专区字幕在线| 日韩有码中文字幕| a级片在线免费高清观看视频| 免费日韩欧美在线观看| 午夜老司机福利片| 高清视频免费观看一区二区| 成年版毛片免费区| 国产97色在线日韩免费| 19禁男女啪啪无遮挡网站| 亚洲精品成人av观看孕妇| 欧美一级毛片孕妇|