• <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)定窗
    一级片免费观看大全| www.熟女人妻精品国产| 18禁美女被吸乳视频| 欧美日本中文国产一区发布| 在线天堂中文资源库| 1024香蕉在线观看| 国产欧美日韩精品亚洲av| 一本一本久久a久久精品综合妖精| 久久毛片免费看一区二区三区| 国产欧美日韩一区二区三| 欧美午夜高清在线| 99re6热这里在线精品视频| 精品人妻熟女毛片av久久网站| 99国产精品免费福利视频| 亚洲人成77777在线视频| 久久99一区二区三区| 国产精品av久久久久免费| 99国产精品一区二区蜜桃av | www.自偷自拍.com| 欧美人与性动交α欧美软件| 欧美日韩av久久| 一级a爱视频在线免费观看| 黑人巨大精品欧美一区二区mp4| 精品国产超薄肉色丝袜足j| 午夜精品国产一区二区电影| 极品少妇高潮喷水抽搐| 精品一区二区三卡| 亚洲第一青青草原| 日韩熟女老妇一区二区性免费视频| 国产有黄有色有爽视频| 男女免费视频国产| 国产av精品麻豆| 久久人妻福利社区极品人妻图片| 1024香蕉在线观看| 在线观看免费视频网站a站| 欧美+亚洲+日韩+国产| 老熟妇乱子伦视频在线观看| www日本在线高清视频| 黄色成人免费大全| 国产成人欧美| 无遮挡黄片免费观看| 两个人免费观看高清视频| 国产不卡av网站在线观看| 国产三级黄色录像| 黄频高清免费视频| 最近最新免费中文字幕在线| 免费在线观看黄色视频的| 18禁裸乳无遮挡动漫免费视频| 女人精品久久久久毛片| 欧美黄色淫秽网站| 人人澡人人妻人| 少妇猛男粗大的猛烈进出视频| 在线十欧美十亚洲十日本专区| 精品福利观看| 日韩中文字幕欧美一区二区| 狠狠精品人妻久久久久久综合| 国产日韩欧美亚洲二区| 欧美在线黄色| 欧美黄色片欧美黄色片| 日韩熟女老妇一区二区性免费视频| 亚洲av欧美aⅴ国产| 别揉我奶头~嗯~啊~动态视频| 一本色道久久久久久精品综合| 国产精品 欧美亚洲| 亚洲情色 制服丝袜| 一区二区三区乱码不卡18| 水蜜桃什么品种好| 动漫黄色视频在线观看| 欧美亚洲日本最大视频资源| 久久人人97超碰香蕉20202| 国产一卡二卡三卡精品| 午夜福利在线免费观看网站| 国产伦人伦偷精品视频| 汤姆久久久久久久影院中文字幕| 久久性视频一级片| 国产熟女午夜一区二区三区| 老司机午夜福利在线观看视频 | 国产免费视频播放在线视频| 老司机深夜福利视频在线观看| 麻豆成人av在线观看| 丝瓜视频免费看黄片| h视频一区二区三区| 一二三四社区在线视频社区8| 精品国产乱码久久久久久男人| 老鸭窝网址在线观看| 免费在线观看黄色视频的| 老司机福利观看| 看免费av毛片| 五月天丁香电影| 亚洲精品成人av观看孕妇| 国产无遮挡羞羞视频在线观看| www.自偷自拍.com| 国产主播在线观看一区二区| 美女高潮到喷水免费观看| 俄罗斯特黄特色一大片| 99国产综合亚洲精品| 国产精品熟女久久久久浪| 日韩中文字幕欧美一区二区| 飞空精品影院首页| 日韩中文字幕欧美一区二区| 男女无遮挡免费网站观看| 丝袜美腿诱惑在线| 女警被强在线播放| 搡老乐熟女国产| 丰满少妇做爰视频| 一区二区三区激情视频| 久久人人97超碰香蕉20202| 国产伦理片在线播放av一区| 中亚洲国语对白在线视频| 成人国产av品久久久| 国产成人一区二区三区免费视频网站| 国产精品免费大片| 男男h啪啪无遮挡| 999久久久精品免费观看国产| av在线播放免费不卡| 天天添夜夜摸| 色精品久久人妻99蜜桃| 老司机午夜十八禁免费视频| 欧美成人免费av一区二区三区 | 国产av又大| 色尼玛亚洲综合影院| 黑人操中国人逼视频| 亚洲欧美日韩另类电影网站| 天天躁狠狠躁夜夜躁狠狠躁| 久久久国产精品麻豆| 久久久久国内视频| 首页视频小说图片口味搜索| 午夜日韩欧美国产| 精品国产亚洲在线| 午夜福利乱码中文字幕| 中文字幕另类日韩欧美亚洲嫩草| 日本一区二区免费在线视频| 久久久国产一区二区| 国产免费视频播放在线视频| 欧美 日韩 精品 国产| 国产精品免费一区二区三区在线 | 另类精品久久| 久久狼人影院| 免费观看av网站的网址| 欧美日韩视频精品一区| 日本一区二区免费在线视频| 中文字幕最新亚洲高清| 女人爽到高潮嗷嗷叫在线视频| 美女国产高潮福利片在线看| 亚洲精品中文字幕在线视频| 亚洲专区中文字幕在线| 国产不卡一卡二| 久久av网站| 黑丝袜美女国产一区| 视频区欧美日本亚洲| 老鸭窝网址在线观看| 三级毛片av免费| 成人国产av品久久久| 亚洲熟女精品中文字幕| 美女扒开内裤让男人捅视频| 香蕉久久夜色| 精品亚洲成a人片在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 成年女人毛片免费观看观看9 | 99久久精品国产亚洲精品| 日日爽夜夜爽网站| 超碰成人久久| 久久久久久久国产电影| 国产在线免费精品| 久久午夜亚洲精品久久| 国产精品 欧美亚洲| 国产成+人综合+亚洲专区| 老汉色∧v一级毛片| 在线观看一区二区三区激情| 日本黄色日本黄色录像| 女警被强在线播放| 久久国产亚洲av麻豆专区| 久久人人97超碰香蕉20202| 精品少妇黑人巨大在线播放| 国产精品 国内视频| 婷婷丁香在线五月| 在线观看舔阴道视频| 国产视频一区二区在线看| 久久久久久久国产电影| 久久精品国产综合久久久| 亚洲三区欧美一区| kizo精华| 久久免费观看电影| 黑人欧美特级aaaaaa片| 亚洲专区国产一区二区| 美女午夜性视频免费| 亚洲av第一区精品v没综合| 老司机午夜福利在线观看视频 | 老熟妇乱子伦视频在线观看| 飞空精品影院首页| 亚洲欧洲日产国产| 啦啦啦中文免费视频观看日本| 在线观看舔阴道视频| 大码成人一级视频| 热re99久久国产66热| 老熟妇仑乱视频hdxx| 亚洲欧美色中文字幕在线| 一级毛片精品| 国产精品熟女久久久久浪| 免费日韩欧美在线观看| avwww免费| 69av精品久久久久久 | 九色亚洲精品在线播放| 亚洲欧美色中文字幕在线| 午夜福利乱码中文字幕| 欧美 日韩 精品 国产| 一二三四社区在线视频社区8| 视频在线观看一区二区三区| 成人18禁在线播放| 久久午夜亚洲精品久久| 99re6热这里在线精品视频| 一夜夜www| 国产成人免费观看mmmm| 欧美国产精品va在线观看不卡| 国产精品久久久av美女十八| 免费女性裸体啪啪无遮挡网站| 少妇的丰满在线观看| 亚洲专区字幕在线| 久久狼人影院| 久久亚洲真实| 少妇的丰满在线观看| 午夜激情av网站| 深夜精品福利| 麻豆国产av国片精品| 三级毛片av免费| 亚洲精品国产一区二区精华液| 一本综合久久免费| 亚洲精品一卡2卡三卡4卡5卡| 日本av手机在线免费观看| 国产区一区二久久| av一本久久久久| 变态另类成人亚洲欧美熟女 | 日韩欧美三级三区| 欧美人与性动交α欧美精品济南到| 中文字幕高清在线视频| 韩国精品一区二区三区| 亚洲精品国产精品久久久不卡| 精品久久久久久电影网| 欧美日韩国产mv在线观看视频| 午夜视频精品福利| 国产一区有黄有色的免费视频| 在线亚洲精品国产二区图片欧美| 少妇精品久久久久久久| 国产一区二区三区综合在线观看| 免费人妻精品一区二区三区视频| 国产精品一区二区在线不卡| 久久久精品区二区三区| 精品人妻熟女毛片av久久网站| 精品国产乱码久久久久久小说| 国产精品99久久99久久久不卡| 一级黄色大片毛片| 高清av免费在线| 久久久久久久久免费视频了| 免费在线观看视频国产中文字幕亚洲| 欧美+亚洲+日韩+国产| 黄频高清免费视频| 中文字幕最新亚洲高清| 老熟妇仑乱视频hdxx| 精品免费久久久久久久清纯 | 成人国语在线视频| av不卡在线播放| 精品一区二区三卡| 亚洲五月婷婷丁香| 国产精品电影一区二区三区 | 大片免费播放器 马上看| 一边摸一边做爽爽视频免费| 在线观看www视频免费| 两人在一起打扑克的视频| 国产片内射在线| 不卡一级毛片| 成人免费观看视频高清| 精品国产一区二区久久| 国产1区2区3区精品| 黄色成人免费大全| av不卡在线播放| 巨乳人妻的诱惑在线观看| 99国产精品一区二区蜜桃av | 亚洲熟女毛片儿| 国产午夜精品久久久久久| 黄片大片在线免费观看| 精品欧美一区二区三区在线| 两性午夜刺激爽爽歪歪视频在线观看 | 一区二区日韩欧美中文字幕| 成年人黄色毛片网站| 精品国内亚洲2022精品成人 | svipshipincom国产片| 国产精品欧美亚洲77777| 51午夜福利影视在线观看| 国产成人免费观看mmmm| 亚洲午夜精品一区,二区,三区| 十八禁网站免费在线| 久久这里只有精品19| 午夜福利影视在线免费观看| 欧美黑人精品巨大| 久久精品亚洲av国产电影网| 日日摸夜夜添夜夜添小说| 免费黄频网站在线观看国产| 国产伦人伦偷精品视频| 亚洲国产av影院在线观看| 男男h啪啪无遮挡| 国产精品影院久久| 国产一区二区激情短视频| 精品熟女少妇八av免费久了| 天天躁日日躁夜夜躁夜夜| 99国产极品粉嫩在线观看| 岛国毛片在线播放| 国产区一区二久久| 我要看黄色一级片免费的| 丁香六月天网| 免费观看人在逋| 黄色丝袜av网址大全| 国产在视频线精品| 丝袜在线中文字幕| 黄色视频,在线免费观看| 国产欧美日韩综合在线一区二区| 国产成人精品久久二区二区免费| 亚洲美女黄片视频| 免费观看av网站的网址| 最黄视频免费看| 80岁老熟妇乱子伦牲交| 老汉色av国产亚洲站长工具| 18在线观看网站| svipshipincom国产片| 手机成人av网站| 亚洲精品美女久久久久99蜜臀| 无人区码免费观看不卡 | 999久久久国产精品视频| 国产男女超爽视频在线观看| 国产精品久久电影中文字幕 | 可以免费在线观看a视频的电影网站| 国产精品亚洲av一区麻豆| 汤姆久久久久久久影院中文字幕| 亚洲精品国产一区二区精华液| 久久久久网色| 日韩欧美一区视频在线观看| 满18在线观看网站| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲专区字幕在线| 午夜福利视频在线观看免费| 男女下面插进去视频免费观看| 两个人看的免费小视频| 午夜精品国产一区二区电影| 少妇裸体淫交视频免费看高清 | 亚洲伊人久久精品综合| 欧美日韩精品网址| 日韩三级视频一区二区三区| 少妇的丰满在线观看| 欧美成人午夜精品| 大香蕉久久网| 老司机午夜福利在线观看视频 | 国产不卡av网站在线观看| av片东京热男人的天堂| 岛国在线观看网站| 一区二区日韩欧美中文字幕| av有码第一页| 国产视频一区二区在线看| 一区二区日韩欧美中文字幕| 视频区图区小说| 欧美久久黑人一区二区| 水蜜桃什么品种好| 在线观看免费视频网站a站| 久久免费观看电影| 亚洲欧美日韩高清在线视频 | 69精品国产乱码久久久| 国产成人啪精品午夜网站| 精品福利观看| 国产在线免费精品| 欧美精品av麻豆av| www.自偷自拍.com| 51午夜福利影视在线观看| 黄色a级毛片大全视频| 日本五十路高清| 日韩大码丰满熟妇| 又黄又粗又硬又大视频| 天天影视国产精品| 色精品久久人妻99蜜桃| 亚洲精品中文字幕一二三四区 | 国产一区二区在线观看av| 脱女人内裤的视频| 大型黄色视频在线免费观看| 日本精品一区二区三区蜜桃| 久久这里只有精品19| av免费在线观看网站| 一二三四社区在线视频社区8| 99久久国产精品久久久| 超碰成人久久| 1024香蕉在线观看| 成年女人毛片免费观看观看9 | 黑人欧美特级aaaaaa片| 亚洲 欧美一区二区三区| www.熟女人妻精品国产| 麻豆乱淫一区二区| 天天添夜夜摸| 丁香六月欧美| 在线观看免费视频网站a站| 欧美人与性动交α欧美精品济南到| 我要看黄色一级片免费的| 日本黄色日本黄色录像| 成人三级做爰电影| 一级黄色大片毛片| 亚洲中文av在线| 9色porny在线观看| 少妇被粗大的猛进出69影院| 成人国语在线视频| 一边摸一边做爽爽视频免费| 亚洲七黄色美女视频| 午夜福利视频精品| 欧美激情 高清一区二区三区| 亚洲免费av在线视频| 美女国产高潮福利片在线看| 人妻一区二区av| 极品人妻少妇av视频| 午夜两性在线视频| 国产男女超爽视频在线观看| 久久久精品94久久精品| 伊人久久大香线蕉亚洲五| 美女高潮到喷水免费观看| 成人18禁在线播放| 黄色 视频免费看| 母亲3免费完整高清在线观看| 成人精品一区二区免费| 久久毛片免费看一区二区三区| 亚洲免费av在线视频| 亚洲专区字幕在线| 九色亚洲精品在线播放| 亚洲av成人一区二区三| 精品一区二区三卡| 99香蕉大伊视频| 欧美日韩精品网址| 亚洲性夜色夜夜综合| 成人18禁高潮啪啪吃奶动态图| 午夜免费成人在线视频| 国产精品成人在线| 国产精品久久久久成人av| 久久精品国产亚洲av香蕉五月 | 久久亚洲真实| 色播在线永久视频| 精品一品国产午夜福利视频| 日本wwww免费看| 精品福利观看| 俄罗斯特黄特色一大片| 国产成人av激情在线播放| 91成人精品电影| 欧美亚洲 丝袜 人妻 在线| 男男h啪啪无遮挡| 久久中文看片网| 制服诱惑二区| 极品教师在线免费播放| 婷婷丁香在线五月| av线在线观看网站| 欧美人与性动交α欧美精品济南到| 久久久欧美国产精品| 午夜福利视频精品| 国产淫语在线视频| 91成人精品电影| 国产一区二区激情短视频| 午夜久久久在线观看| 国产一区有黄有色的免费视频| 人人妻人人澡人人看| 亚洲中文字幕日韩| 在线看a的网站| 啦啦啦在线免费观看视频4| 另类亚洲欧美激情| 在线av久久热| 成人国产一区最新在线观看| 久久久久国内视频| 国产福利在线免费观看视频| e午夜精品久久久久久久| 久久婷婷成人综合色麻豆| 又大又爽又粗| 两性午夜刺激爽爽歪歪视频在线观看 | 狠狠精品人妻久久久久久综合| 9191精品国产免费久久| 国产三级黄色录像| 久久青草综合色| 99精品欧美一区二区三区四区| 成人影院久久| 国产xxxxx性猛交| 操出白浆在线播放| 麻豆乱淫一区二区| cao死你这个sao货| 王馨瑶露胸无遮挡在线观看| 亚洲精品粉嫩美女一区| av国产精品久久久久影院| 十八禁网站网址无遮挡| 国产精品一区二区在线不卡| 国产精品久久久av美女十八| 欧美日韩福利视频一区二区| 夜夜骑夜夜射夜夜干| 黄色毛片三级朝国网站| 欧美黑人精品巨大| 国产亚洲欧美精品永久| 亚洲国产看品久久| www.精华液| 精品视频人人做人人爽| 狂野欧美激情性xxxx| 狠狠精品人妻久久久久久综合| 日本一区二区免费在线视频| 老汉色∧v一级毛片| 欧美日韩成人在线一区二区| 一边摸一边做爽爽视频免费| 日韩人妻精品一区2区三区| 欧美人与性动交α欧美软件| 国产成人精品在线电影| 操美女的视频在线观看| 黄色毛片三级朝国网站| 国产精品九九99| 成人精品一区二区免费| 色94色欧美一区二区| 色综合婷婷激情| 日韩视频一区二区在线观看| 黑人操中国人逼视频| 国产精品自产拍在线观看55亚洲 | 久久久精品94久久精品| 精品一品国产午夜福利视频| 国产男女内射视频| 国产真人三级小视频在线观看| 婷婷丁香在线五月| 国产xxxxx性猛交| 可以免费在线观看a视频的电影网站| 国产成人系列免费观看| av福利片在线| 老汉色∧v一级毛片| 操出白浆在线播放| 大型黄色视频在线免费观看| 免费在线观看黄色视频的| 在线亚洲精品国产二区图片欧美| 最黄视频免费看| 日韩免费av在线播放| 国产麻豆69| 丝袜喷水一区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲av成人不卡在线观看播放网| av网站免费在线观看视频| 99久久99久久久精品蜜桃| 国产精品自产拍在线观看55亚洲 | 黄色怎么调成土黄色| 亚洲成av片中文字幕在线观看| 国产精品电影一区二区三区 | 欧美日韩中文字幕国产精品一区二区三区 | 国产男女超爽视频在线观看| 露出奶头的视频| 欧美乱妇无乱码| 欧美成人午夜精品| 亚洲五月婷婷丁香| 国产精品亚洲一级av第二区| 汤姆久久久久久久影院中文字幕| 欧美激情久久久久久爽电影 | 欧美亚洲 丝袜 人妻 在线| 激情在线观看视频在线高清 | 久久久久国内视频| 国产精品免费视频内射| 母亲3免费完整高清在线观看| 欧美激情久久久久久爽电影 | 一区二区三区精品91| 国产精品熟女久久久久浪| 国产成人精品在线电影| 亚洲avbb在线观看| 日韩免费av在线播放| 国产成人av教育| 变态另类成人亚洲欧美熟女 | 精品国产乱码久久久久久男人| 色播在线永久视频| videosex国产| 国产在视频线精品| 中文字幕人妻丝袜制服| 五月天丁香电影| 一夜夜www| 99国产精品一区二区三区| 亚洲一码二码三码区别大吗| 老司机在亚洲福利影院| 亚洲精品美女久久av网站| 亚洲av美国av| 精品午夜福利视频在线观看一区 | 国产一区有黄有色的免费视频| 午夜激情久久久久久久| 男女无遮挡免费网站观看| 亚洲欧洲日产国产| 欧美黑人欧美精品刺激| 久久国产亚洲av麻豆专区| 老司机影院毛片| 王馨瑶露胸无遮挡在线观看| 亚洲人成77777在线视频| 无限看片的www在线观看| 国产福利在线免费观看视频| 国产老妇伦熟女老妇高清| 精品国产一区二区三区久久久樱花| 又大又爽又粗| 亚洲精品成人av观看孕妇| 日韩熟女老妇一区二区性免费视频| 这个男人来自地球电影免费观看| 亚洲欧美精品综合一区二区三区| a在线观看视频网站| 丝袜在线中文字幕| 国产成人av教育| 天天影视国产精品| 夫妻午夜视频| a级片在线免费高清观看视频| 高清在线国产一区| 日本黄色日本黄色录像| 天天躁日日躁夜夜躁夜夜| 99国产精品一区二区蜜桃av | 一级毛片女人18水好多| 色94色欧美一区二区| 精品久久久精品久久久| 欧美乱码精品一区二区三区| 99在线人妻在线中文字幕 | 亚洲欧洲日产国产| 日韩欧美三级三区| 精品国内亚洲2022精品成人 | 美女午夜性视频免费| 正在播放国产对白刺激| 99在线人妻在线中文字幕 | 精品卡一卡二卡四卡免费| 亚洲精品中文字幕一二三四区 |