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

    多孔介質(zhì)中甲烷水合物聚散過程的交流阻抗譜響應(yīng)特征

    2016-08-17 01:42:57金學(xué)彬邢蘭昌劉昌嶺鄭金吾
    天然氣工業(yè) 2016年3期
    關(guān)鍵詞:水合物甲烷飽和度

    金學(xué)彬 陳 強(qiáng) 邢蘭昌 劉昌嶺 鄭金吾

    1.中國石油大學(xué)(華東)信息與控制工程學(xué)院 2.國土資源部天然氣水合物重點(diǎn)實(shí)驗(yàn)室·青島海洋地質(zhì)研究所3.海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室海洋礦產(chǎn)資源評(píng)價(jià)與探測(cè)技術(shù)功能實(shí)驗(yàn)室

    多孔介質(zhì)中甲烷水合物聚散過程的交流阻抗譜響應(yīng)特征

    金學(xué)彬1,2陳 強(qiáng)2,3邢蘭昌1劉昌嶺2,3鄭金吾1

    1.中國石油大學(xué)(華東)信息與控制工程學(xué)院2.國土資源部天然氣水合物重點(diǎn)實(shí)驗(yàn)室·青島海洋地質(zhì)研究所3.海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室海洋礦產(chǎn)資源評(píng)價(jià)與探測(cè)技術(shù)功能實(shí)驗(yàn)室

    金學(xué)彬等.多孔介質(zhì)中甲烷水合物聚散過程的交流阻抗譜響應(yīng)特征.天然氣工業(yè),2016,36(3):120-127.

    電學(xué)特性實(shí)驗(yàn)是探索天然氣水合物(以下簡(jiǎn)稱水合物)生成分解過程及其動(dòng)力學(xué)特征的有效技術(shù)手段。較之于過去單一的僅對(duì)電阻率參數(shù)進(jìn)行測(cè)試的實(shí)驗(yàn)方法,交流阻抗譜測(cè)試方法可以獲得一定頻率范圍內(nèi)被測(cè)體系的阻抗(包括電阻、電抗)特性變化規(guī)律,從而有助于更全面深入地認(rèn)識(shí)水合物的電學(xué)性質(zhì)特征。為此,設(shè)計(jì)了甲烷水合物交流阻抗譜測(cè)試系統(tǒng),以甲烷水合物—多孔介質(zhì)—濃度為3.5%的鹽水為研究體系,測(cè)試了水合物在多孔介質(zhì)內(nèi)合成與分解過程中研究體系的阻抗譜,并分析了阻抗譜特征參數(shù)隨水合物飽和度變化的規(guī)律。結(jié)果表明:①采用激發(fā)源的不同頻率測(cè)試的結(jié)果差異明顯,在低測(cè)試頻率(0.1~10.6 Hz)下,測(cè)試體系的水合物飽和度與阻抗譜參數(shù)之間未表現(xiàn)出明顯的規(guī)律性關(guān)系,而在中高測(cè)試頻率(101.8 Hz~1 MHz)下 ,阻抗譜參數(shù)隨著測(cè)試體系水合物飽和度的增加而增加;②當(dāng)測(cè)試頻率為101.8 Hz、1 114 Hz和10 746 Hz時(shí),在一定的誤差允許范圍內(nèi),測(cè)試體系的阻抗譜參數(shù)與水合物飽和度之間符合單調(diào)遞增的線性關(guān)系;③阻抗幅值的頻散特性與水合物飽和度之間存在著明顯的對(duì)應(yīng)關(guān)系。

    天然氣水合物 電學(xué)特性 合成與分解過程 交流阻抗譜 測(cè)試方法 多孔介質(zhì) 水合物飽和度 頻散特性

    天然氣水合物(以下簡(jiǎn)稱水合物)在海底沉積物中形成時(shí),其儲(chǔ)層的理化特性會(huì)出現(xiàn)高電阻率、小聲波時(shí)差和低自然電位幅度等異?,F(xiàn)象。其中,電阻率異?,F(xiàn)象不僅可以用來判斷含水合物的儲(chǔ)層位置,也可以結(jié)合Archie公式來估算水合物的飽和度[1-7]。水合物電學(xué)特性實(shí)驗(yàn)研究不僅是探索水合物生成分解過程及其動(dòng)力學(xué)特征的有效技術(shù)手段,也能夠?yàn)橐巴獾乃衔锟碧胶唾Y源評(píng)價(jià)提供基礎(chǔ)和校正數(shù)據(jù)。Buffett等[8]利用電阻率的變化速度來研究水合物形成的動(dòng)力學(xué)過程,認(rèn)為電阻法測(cè)量可以量化水合物的成核規(guī)模。Spangenberg等[9-11]通過四電極法測(cè)量了由溶解氣在均勻玻璃珠表面形成甲烷水合物的電阻率,指出所形成的水合物樣品的電阻率增大指數(shù)和含水飽和度在雙對(duì)數(shù)坐標(biāo)下不一定是線性關(guān)系,即存在非阿爾奇現(xiàn)象。陳強(qiáng)等[12]應(yīng)用自行設(shè)計(jì)的實(shí)驗(yàn)裝置進(jìn)行了含水合物沉積物體系電學(xué)特性的研究,主要對(duì)多孔介質(zhì)中CO2水合物的生成與分解過程進(jìn)行分析討論,結(jié)果表明,水合物生成過程中多孔介質(zhì)體系阻抗增大,分解過程中阻抗減小,且阻抗變化與系統(tǒng)的溫度、壓力變化相互對(duì)應(yīng)。陳玉鳳等[13]對(duì)南海沉積物中水合物飽和度與電阻率的關(guān)系進(jìn)行了實(shí)驗(yàn)研究,結(jié)果表明對(duì)于含水合物的沉積物,其雙對(duì)數(shù)坐標(biāo)系的電阻率增大指數(shù)和含水飽和度并不是阿爾奇公式所描述的線性關(guān)系,其飽和度指數(shù)n不是定值,而隨含水飽和度Sw的增加而增加,并認(rèn)為體系電阻率與水合物在孔隙中的分布狀態(tài)有關(guān)。

    截至目前,不少研究者針對(duì)多孔介質(zhì)中水合物的電學(xué)特性開展了一定的工作,但是大多數(shù)的研究局限于對(duì)電阻率參數(shù)進(jìn)行測(cè)試,即采用固定頻率/直流電源的信號(hào)作為激勵(lì)源對(duì)水合物體系進(jìn)行測(cè)試。然而,現(xiàn)有的研究結(jié)果表明,影響含水合物沉積物電學(xué)性質(zhì)的因素眾多、機(jī)理復(fù)雜,不同頻率/類型的激勵(lì)源所獲得的電阻率數(shù)據(jù)及其反演的水合物飽和度尚存在一定的不確定性[14]。近年來,交流阻抗技術(shù)逐漸被應(yīng)用到地球科學(xué)研究領(lǐng)域。交流阻抗測(cè)試方法可以獲得一定頻率范圍內(nèi)被測(cè)體系的阻抗(包括電阻、電抗)特性變化規(guī)律,從而有助于更全面深入地認(rèn)識(shí)被測(cè)物質(zhì)電學(xué)性質(zhì)的特征。Du Frane等[15]在高壓反應(yīng)釜中測(cè)量了20 Hz~2 MHz頻率范圍內(nèi)甲烷水合物的阻抗譜,在對(duì)阻抗譜進(jìn)行分析的基礎(chǔ)上給出了甲烷水合物的電導(dǎo)率。綜上所述,電學(xué)參數(shù)是水合物的一個(gè)極其重要的性質(zhì)參數(shù)。但是由于受水合物特殊的實(shí)驗(yàn)條件以及復(fù)雜的生成機(jī)理影響,含水合物沉積物電學(xué)性質(zhì)的變化規(guī)律較為復(fù)雜,目前被廣泛采用的電阻率參數(shù)還不足以充分刻畫其電學(xué)性質(zhì);此外,采用交流阻抗法對(duì)含水合物沉積物全阻抗特性進(jìn)行研究仍處于探索階段。

    針對(duì)以上問題,筆者提出用交流阻抗法對(duì)含水合物沉積物的全阻抗特性進(jìn)行實(shí)驗(yàn)研究,測(cè)試水合物在多孔介質(zhì)內(nèi)合成與分解過程中體系的阻抗譜,探索阻抗譜特征參數(shù)隨水合物飽和度變化的規(guī)律,為進(jìn)一步開發(fā)天然氣水合物飽和度測(cè)量方法和探索水合物生成分解過程的動(dòng)力學(xué)機(jī)理提供參考。

    1 實(shí)驗(yàn)技術(shù)與方法

    1.1技術(shù)方法原理

    筆者提出對(duì)含水合物多孔介質(zhì)體系進(jìn)行全阻抗電學(xué)特性(電阻特性和電容特性)研究,采用電化學(xué)測(cè)試技術(shù)中的交流阻抗法來實(shí)現(xiàn)對(duì)一定頻率范圍內(nèi)的體系阻抗譜的測(cè)試。交流阻抗法指通過控制電化學(xué)系統(tǒng)的電流(或系統(tǒng)的電勢(shì))在小幅度的條件下隨時(shí)間按正弦規(guī)律變化,同時(shí)測(cè)量相應(yīng)的系統(tǒng)電勢(shì)(或電流)隨時(shí)間的變化,或者直接測(cè)量系統(tǒng)的交流阻抗(或?qū)Ъ{),進(jìn)而分析電化學(xué)系統(tǒng)的反應(yīng)機(jī)理、計(jì)算系統(tǒng)的相關(guān)參數(shù)[16]。

    首先以電化學(xué)工作站與新型高壓反應(yīng)釜為基礎(chǔ)搭建了阻抗譜測(cè)試實(shí)驗(yàn)系統(tǒng),然后以甲烷水合物—多孔介質(zhì)—濃度為3.5%的鹽水為研究體系,獲取了甲烷水合物在多孔介質(zhì)內(nèi)合成及分解過程中飽和度和阻抗譜參數(shù)等數(shù)據(jù)。在分析處理實(shí)驗(yàn)數(shù)據(jù)的基礎(chǔ)上,評(píng)價(jià)不同激勵(lì)源頻率的響應(yīng)特征,并對(duì)某些測(cè)試頻率下的阻抗譜參數(shù)與水合物飽和度數(shù)據(jù)進(jìn)行曲線擬合,探討兩者之間的規(guī)律性關(guān)系。

    1.2實(shí)驗(yàn)裝置

    圖1 天然氣水合物阻抗譜測(cè)量系統(tǒng)示意圖

    實(shí)驗(yàn)系統(tǒng)(圖1)主要由以下幾部分組成:高壓反應(yīng)釜、低溫恒溫箱、溫度壓力測(cè)量傳感器、數(shù)據(jù)采集器、電化學(xué)工作站以及用于控制電化學(xué)工作站并進(jìn)行數(shù)據(jù)存儲(chǔ)、顯示和處理的計(jì)算機(jī)。

    反應(yīng)釜的最大容積為942.5 mL,耐壓20 MPa,其外壁為雙層結(jié)構(gòu),夾層允許循環(huán)水流動(dòng)以對(duì)反應(yīng)釜進(jìn)行制冷,內(nèi)壁襯以具有絕緣和耐腐蝕作用的聚四氟乙烯。反應(yīng)釜的上端蓋可以拆卸,在上端蓋上鉆孔以便于連接壓力、溫度傳感器以及電極。低溫恒溫槽采購自杭州雪中炭恒溫技術(shù)有限公司,最低溫度可達(dá)-30 ℃,控溫精度為±0.05 ℃。高壓反應(yīng)釜的氣體壓力可通過壓力傳感器來記錄,壓力傳感器的量程為0~40 MPa,精度為0.1級(jí)。從反應(yīng)釜上端插入釜內(nèi)的Pt100熱電阻記錄反應(yīng)體系的溫度,其量程為-50~260 ℃,精度為0.1級(jí)。整個(gè)系統(tǒng)的溫度、壓力由數(shù)據(jù)采集系統(tǒng)Agilent 34972A采集顯示。電化學(xué)工作站采用德國Zahner公司的IM6,其可實(shí)現(xiàn)的最大頻率范圍為:10 mHz~8 MHz。電化學(xué)工作站的控制、阻抗譜參數(shù)的采集與處理均通過配以IM6專用軟件(Thales Software Package)的計(jì)算機(jī)完成。

    1.3實(shí)驗(yàn)步驟

    實(shí)驗(yàn)材料包括甲烷氣體、去離子水、氯化鈉和天然海砂,具體規(guī)格如表1所示。

    表1 實(shí)驗(yàn)材料及其規(guī)格表

    實(shí)驗(yàn)分為合成過程與分解過程兩個(gè)階段,具體實(shí)驗(yàn)步驟如下:

    1)用去離子水清洗高壓反應(yīng)釜并烘干,將洗凈烘干后的天然海砂裝填于高壓反應(yīng)釜,依據(jù)溫度傳感器的長(zhǎng)度確定海砂填充量;使用去離子水配置氯化鈉溶液,濃度為3.5%。

    2)將氯化鈉溶液緩慢加到海砂里,并確保此時(shí)的海砂處于水飽和狀態(tài)(海砂孔隙度為40%);將工作電極、輔助電極和Pt100溫度傳感器插入釜內(nèi)海砂中,蓋上上端蓋并進(jìn)行固定。

    3)通過電極引線將工作電極和輔助電極與IM6主機(jī)連接,在軟件中對(duì)阻抗譜測(cè)試環(huán)境進(jìn)行設(shè)置:掃描模式選擇從高頻往低頻掃描,頻率范圍設(shè)為0.1 Hz~1 MHz,測(cè)量頻率范圍內(nèi)每一個(gè)數(shù)量級(jí)的取點(diǎn)數(shù)以66 Hz為界限,分別設(shè)置為4個(gè)點(diǎn)(下限)和8個(gè)點(diǎn)(大于 66 Hz);將溫度傳感器和壓力傳感器通過安捷倫數(shù)據(jù)采集器與計(jì)算機(jī)相連接,并對(duì)數(shù)據(jù)采集通道進(jìn)行配置。

    4)對(duì)系統(tǒng)抽真空,注入甲烷氣體至設(shè)定壓力8.5 MPa,室溫下放置24 h,以使甲烷在海砂孔隙水中充分?jǐn)U散和溶解;設(shè)置恒溫箱溫度為20 ℃,對(duì)被測(cè)反應(yīng)系統(tǒng)進(jìn)行恒溫測(cè)試,同時(shí)采集記錄溫度、壓力和阻抗譜數(shù)據(jù),溫度、壓力數(shù)據(jù)采集頻率為0.1 Hz,阻抗譜每隔15 min采集一次。

    5)設(shè)置恒溫箱溫度為2 ℃,開始水合物生成過程;當(dāng)反應(yīng)體系的溫度、壓力基本不變時(shí),合成過程結(jié)束。采集記錄整個(gè)過程反應(yīng)體系內(nèi)的溫度、壓力和阻抗譜數(shù)據(jù)。

    6)水合物分解過程采取分步升溫的方式進(jìn)行,即逐漸升高溫度控制器設(shè)置溫度,每設(shè)置一個(gè)溫度值后等待足夠長(zhǎng)時(shí)間使系統(tǒng)達(dá)到新的相平衡(如40 h,使溫度壓力持續(xù)保持平穩(wěn)),然后采集1 h的阻抗譜數(shù)據(jù),同時(shí)對(duì)體系溫度、壓力數(shù)據(jù)進(jìn)行采集記錄。

    2 實(shí)驗(yàn)結(jié)果與討論

    2.1體系的溫度、壓力及飽和度變化趨勢(shì)

    假設(shè)甲烷水合物的化學(xué)式為CH4?nH2O,n為水合指數(shù),對(duì)n取值為6.0進(jìn)行相關(guān)計(jì)算[17],則水合物飽和度的計(jì)算式為[18]:

    式中Sh表示多孔介質(zhì)中水合物飽和度;p1表示系統(tǒng)初始?jí)毫?,Pa;Vg表示反應(yīng)釜?dú)庀囿w積,m3;p2表示水合物生成(分解)過程中的系統(tǒng)壓力,Pa;Z1表示系統(tǒng)初始狀態(tài)下氣體的壓縮因子;R表示摩爾氣體常數(shù),8.314 J/(mol·℃);T1表示系統(tǒng)初始溫度,℃;Z2表示水合物生成(分解)過程中氣體的壓縮因子;T2表示水合物生成(分解)過程中的系統(tǒng)溫度,℃;Mw表示水的摩爾質(zhì)量,18 g/mol;mw表示反應(yīng)體系初始加入的鹽水總質(zhì)量。

    圖2是水合物合成實(shí)驗(yàn)過程中,反應(yīng)體系內(nèi)溫度、壓力及飽和度的變化趨勢(shì)圖。如圖2-a所示,合成過程前后歷經(jīng)大約300 h,在水合物生成實(shí)驗(yàn)開始階段,反應(yīng)釜內(nèi)的溫度受恒溫箱的控制而迅速降低,與此同時(shí)反應(yīng)釜內(nèi)壓力也隨之降低。當(dāng)溫度降到約4 ℃時(shí),壓力開始急劇下降,說明反應(yīng)體系內(nèi)甲烷水合物開始生成。之后,伴隨著水合物的大量生成,反應(yīng)體系內(nèi)水合物飽和度持續(xù)增加,當(dāng)反應(yīng)進(jìn)行到約第99 h時(shí)飽和度達(dá)到22.13%。

    圖2 甲烷水合物合成過程溫度、壓力及飽和度變化曲線圖

    為了更加清晰地展示在水合物合成過程中反應(yīng)體系內(nèi)發(fā)生的一系列特征反應(yīng)現(xiàn)象,將合成過程中的前50 h數(shù)據(jù)進(jìn)行細(xì)節(jié)放大,即如圖2-b所示的合成過程初期反應(yīng)體系內(nèi)溫度、壓力及飽和度變化曲線。溫度、壓力共同標(biāo)示著合成過程的不同階段,從開始生成階段到大量生成階段,再到緩慢進(jìn)行階段,最后趨于穩(wěn)定。大量生成階段對(duì)應(yīng)于一個(gè)明顯的溫度特征現(xiàn)象,即當(dāng)水合物大量生成時(shí)會(huì)在短時(shí)間內(nèi)散發(fā)出大量的熱,從而引起體系內(nèi)溫度明顯上升。由于恒溫箱的溫度控制作用,此溫度上升現(xiàn)象持續(xù)時(shí)間不長(zhǎng),體系內(nèi)的溫度此后逐漸下降至設(shè)定溫度。伴隨著整個(gè)合成反應(yīng)過程的進(jìn)行,體系內(nèi)水合物飽和度持續(xù)增大。

    水合物分解實(shí)驗(yàn)過程采用了分步升溫的方式,即通過控制體系的溫度來控制水合物的飽和度,繼而測(cè)試不同飽和度狀態(tài)下體系的阻抗譜,為后續(xù)分析水合物飽和度與阻抗譜特征參數(shù)之間的定量關(guān)系提供基礎(chǔ)。圖3是水合物分解過程中,反應(yīng)體系內(nèi)溫度、壓力隨著時(shí)間的變化趨勢(shì)圖。如圖3所示,溫度在恒溫箱的調(diào)節(jié)下,每上升一個(gè)梯度,將會(huì)打破體系內(nèi)原有的相平衡,從而引起一定量的水合物發(fā)生分解,緊接著體系內(nèi)壓力升高,于是再一次達(dá)到新的平衡。表2列出了分解過程測(cè)試的不同狀態(tài)點(diǎn)對(duì)應(yīng)的溫度、壓力及飽和度數(shù)據(jù)。

    圖3 甲烷水合物分解過程體系內(nèi)溫度、壓力變化曲線圖

    表2 分解過程測(cè)試的狀態(tài)點(diǎn)相關(guān)參數(shù)表

    2.2飽和度與阻抗譜特征參數(shù)的關(guān)系

    圖4是水合物合成過程于不同測(cè)試頻率下反應(yīng)體系阻抗譜參數(shù)變化圖。圖4中,0.1 Hz、1.1 Hz與10.6 Hz的結(jié)果是以主坐標(biāo)軸(左縱軸)展示的,其余頻率的結(jié)果是以次坐標(biāo)軸展示的。如圖4所示,在低測(cè)試頻率(0.1 ~10.6 Hz)下,隨著反應(yīng)體系狀態(tài)的變化,測(cè)試體系的阻抗譜參數(shù)變化很明顯,即測(cè)試體系阻抗譜參數(shù)能夠靈敏地反映出體系內(nèi)飽和度的變化;但是,在實(shí)驗(yàn)開始的一段時(shí)間內(nèi),出現(xiàn)了阻抗譜參數(shù)值波動(dòng)且呈下降趨勢(shì)的現(xiàn)象,推測(cè)是由以下兩個(gè)方面引起的:①水合物合成過程的初期,反應(yīng)體系內(nèi)各相態(tài)還較為混亂,因此會(huì)導(dǎo)致測(cè)試體系阻抗值產(chǎn)生瞬時(shí)波動(dòng);②在低測(cè)試頻率下,體系的容抗效應(yīng)會(huì)突顯出來。本實(shí)驗(yàn)采用的測(cè)量體系的容抗效應(yīng)可以類比成平行板電容器模型來討論,據(jù)報(bào)道[19],當(dāng)實(shí)驗(yàn)中電極總面積與電極距離固定不變時(shí),被測(cè)樣品的綜合表觀介電常數(shù)成為影響電容量變化的決定因素,而介質(zhì)的綜合表觀介電常數(shù)主要由介質(zhì)溫度、各相飽和度、密度、礦化度、鹽度以及多孔介質(zhì)特性等共同決定。因此容抗會(huì)存在一定的不確定性。當(dāng)測(cè)試頻率增大到1 114 Hz時(shí),可以發(fā)現(xiàn),隨著體系內(nèi)合成反應(yīng)的進(jìn)行,檢測(cè)到的測(cè)試體系阻抗值呈單調(diào)增大的趨勢(shì)。當(dāng)測(cè)試頻率增大到邊界值附近時(shí),在較長(zhǎng)的一段時(shí)間內(nèi),阻抗譜參數(shù)值波動(dòng)明顯,可能是因?yàn)橐环矫鎸?shí)驗(yàn)初期,反應(yīng)體系還未達(dá)到穩(wěn)定狀態(tài);另一方面在邊界值附近頻率點(diǎn),實(shí)驗(yàn)設(shè)備可能出現(xiàn)了不可預(yù)測(cè)的干擾。

    圖4 水合物合成過程中不同測(cè)試頻率下的阻抗譜變化曲線圖

    圖5是水合物分解過程于不同測(cè)試頻率下反應(yīng)體系阻抗譜參數(shù)與水合物飽和度對(duì)應(yīng)關(guān)系圖。圖5中,0.1 Hz、1.1 Hz與10.6 Hz的結(jié)果是以主坐標(biāo)軸(左縱軸)展示的,其余頻率的結(jié)果是以次坐標(biāo)軸展示的。

    圖5 水合物分解過程中不同測(cè)試頻率下水合物飽和度與阻抗譜的關(guān)系曲線圖

    如圖5所示,在低測(cè)試頻率(0.1~10.6 Hz)下,測(cè)試體系水合物飽和度與阻抗譜參數(shù)之間未表現(xiàn)出明顯的規(guī)律性關(guān)系;在中高頻段測(cè)試頻率(101.8 Hz~1 MHz)下,阻抗譜參數(shù)隨著測(cè)試體系水合物飽和度的增加而單調(diào)增加,這與已報(bào)道過的含水合物儲(chǔ)層具有高電阻率特點(diǎn)等相關(guān)結(jié)論相吻合[20]。

    將測(cè)試頻率的中高頻段(101.8 Hz~1 MHz)內(nèi)具代表性的頻率點(diǎn)下得到的測(cè)試體系阻抗譜參數(shù)與水合物飽和度之間的關(guān)系進(jìn)行線性擬合,結(jié)果如圖6所示。從圖6可以看出,在一定的誤差允許范圍內(nèi),測(cè)試體系的阻抗譜參數(shù)與水合物飽和度之間符合單調(diào)增的線性關(guān)系。從各個(gè)結(jié)果對(duì)比分析來看,首先,不同測(cè)試頻率點(diǎn)處阻抗譜參數(shù)與水合物飽和度關(guān)系曲線的單調(diào)性和線性度是有差異的,尤其低頻段(0.1~10.6 Hz)與中高頻段(101.8 Hz~1 MHz)差異明顯,這正體現(xiàn)了測(cè)試體系的頻散特性;其次,同頻段內(nèi)的不同頻率點(diǎn)之間的結(jié)果又遵循相似的規(guī)律。綜上所述,通過選擇最為合適的測(cè)試頻率(1 114 Hz)對(duì)反應(yīng)體系的阻抗譜參數(shù)進(jìn)行測(cè)試,應(yīng)用測(cè)試結(jié)果可以推算出體系的水合物飽和度值,即可用于開發(fā)水合物飽和度測(cè)量和計(jì)算方法。

    圖6 水合物分解過程中水合物飽和度與特征頻率點(diǎn)處阻抗的關(guān)系曲線圖

    2.3甲烷水合物合成與分解過程中水合物飽和度與414 180 Hz下阻抗的關(guān)系

    電磁波電阻率測(cè)井常用的低頻工作頻率為400 kHz[21],為了考察此測(cè)試頻率下反應(yīng)體系阻抗與水合物飽和度之間的關(guān)系,對(duì)測(cè)試頻率為414 180 Hz時(shí),甲烷水合物合成與分解過程展開研究。圖7是水合物合成過程中反應(yīng)體系內(nèi)水合物飽和度以及測(cè)試頻率為414 180 Hz下樣品的阻抗值隨時(shí)間變化的趨勢(shì)圖。如圖7所示,隨著體系內(nèi)合成反應(yīng)的進(jìn)行,檢測(cè)到的測(cè)試體系阻抗值呈單調(diào)增大的趨勢(shì),并且可以很好地與水合物飽和度變化曲線對(duì)應(yīng)起來。

    圖7 水合物合成過程中水合物飽和度與阻抗(414 180 Hz)隨時(shí)間變化曲線圖

    將水合物分解過程中測(cè)試頻率為414 180 Hz下反應(yīng)體系阻抗與水合物飽和度之間的關(guān)系進(jìn)行線性擬合(圖8),如圖8所示,在一定誤差允許范圍內(nèi),測(cè)試體系的阻抗與水合物飽和度之間也符合單調(diào)增的線性關(guān)系。該結(jié)果顯示了此頻率段在水合物飽和度測(cè)量上的潛力,這對(duì)后續(xù)開展水合物野外勘探研究具有一定的參考價(jià)值和指導(dǎo)意義。

    圖8 水合物分解過程中水合物飽和度與阻抗(414 180 Hz)的關(guān)系曲線圖

    2.4不同水合物飽和度下測(cè)試體系的頻散特性

    在交變電場(chǎng)的作用下,測(cè)試體系的阻抗隨著測(cè)量頻率的變化而變化,體現(xiàn)了測(cè)試體系的頻散特性,其受飽和度、測(cè)量頻率、溫度、壓力等諸多因素的影響[22],為此對(duì)不同水合物飽和度下測(cè)試體系的頻散特性展開初步研究。圖9給出了測(cè)試體系在不同水合物飽和度狀態(tài)下阻抗幅值以及相角隨著測(cè)試頻率的變化曲線。如圖9-a所示,隨著水合物飽和度的增大,測(cè)試體系孔隙流體的導(dǎo)電能力降低,致使阻抗幅值增大;與此同時(shí),測(cè)試體系阻抗幅值的頻散程度(低頻和高頻下體系阻抗的模之差)也相應(yīng)增加。如圖9-b所示,隨著水合物飽和度的增大,相角的值是相應(yīng)減小的,同時(shí)相角的頻散程度(低頻和高頻下體系相角的絕對(duì)值之差)未呈現(xiàn)出特定的規(guī)律性??梢?,多孔介質(zhì)體系阻抗幅值的頻散特性與水合物飽和度之間存在明顯的對(duì)應(yīng)關(guān)系,受水合物飽和度的影響較大。這可能是因?yàn)榧淄樗衔锏拇嬖冢淖兞硕嗫捉橘|(zhì)體系整體的骨架結(jié)構(gòu),使得孔隙結(jié)構(gòu)變得更為復(fù)雜,從而影響了導(dǎo)電流體在孔隙中的微觀分布狀態(tài),進(jìn)而促使離子交換特性也發(fā)生改變,最終影響了測(cè)試體系的頻散特性[22]。

    圖9 不同水合物飽和度下多孔介質(zhì)體系的頻散特性曲線圖

    3 結(jié)論與展望

    1)在自主研發(fā)的水合物電學(xué)實(shí)驗(yàn)裝置上對(duì)甲烷水合物的生成和分解過程進(jìn)行了實(shí)驗(yàn)測(cè)試研究,獲得了甲烷水合物合成與分解過程中反應(yīng)體系的阻抗譜參數(shù),采用基于最小二乘法的線性擬合方式建立了水合物飽和度與阻抗譜參數(shù)兩者之間的關(guān)系,探討了不同水合物飽和度下測(cè)試體系的頻散特性。

    2)研究發(fā)現(xiàn),使用激發(fā)源的不同頻率測(cè)試的結(jié)果差異明顯,如當(dāng)測(cè)試頻率為10.6 Hz時(shí),測(cè)試體系飽和度與阻抗譜參數(shù)之間未表現(xiàn)出明顯的規(guī)律性關(guān)系,當(dāng)測(cè)試頻率為1 114 Hz時(shí),在一定的誤差允許范圍內(nèi),測(cè)試體系的阻抗譜參數(shù)與水合物飽和度之間符合單調(diào)增的線性關(guān)系;阻抗幅值的頻散特性與水合物飽和度之間存在明顯的對(duì)應(yīng)關(guān)系。

    3)交流阻抗譜測(cè)試技術(shù)是研究水合物生成和分解過程電學(xué)特性的一種有效方法,但尚存一些問題需要解決:研究和尋求更系統(tǒng)的數(shù)據(jù)解釋方法,深入地探討含水合物沉積物體系的頻散特性,揭示飽和度與頻散特性之間的規(guī)律性關(guān)系;從技術(shù)上實(shí)現(xiàn)水合物模擬實(shí)驗(yàn)環(huán)境與野外水合物賦存環(huán)境更為相近,從而使得實(shí)驗(yàn)研究結(jié)果更具有實(shí)際應(yīng)用價(jià)值。

    [1]Collett TS. Well log evaluation of natural gas hydrates[R]. Denver, Colorado: U. S. Department of Energy, 1992.

    [2]Collett TS, Ladd J. Detection of gas hydrate with downhole logs and assessment of gas hydrate Concentrations (saturations) and gas volumes on the Blake Ridge with electrical resistivity log data[C]// Proceedings of the Ocean Drilling Program:Scientifi c Results, 2000, 164: 179-191.

    [3]李小森, 馮景春, 李剛, 王屹. 電阻率在天然 氣水合物三維生成及開采過程中的變化特性模擬實(shí)驗(yàn)[J]. 天然氣工業(yè), 2013, 33(7): 18-23. Li Xiaosen, Feng Jingchun, Li Gang, Wang Yi. An experimental study of resistivity variation in the 3D simulation of methane hydrate generation and production[J]. Natural Gas Industry, 2013, 33(7): 18-23.

    [4]粟科華, 孫長(zhǎng)宇, 李楠, 鐘小禹, 陳光進(jìn). 天然氣水合物三維成藏物模實(shí)驗(yàn)系統(tǒng)的構(gòu)建與檢驗(yàn)[J]. 天然氣工業(yè), 2013, 33(12): 173-178. Su Kehua, Sun Changyu, Li Nan, Zhong Xiaoyu, Chen Guangjin. A large-scale three-dimensional device for simulating natural gas hydrates accumulation and distribution process in sediments[J]. Natural Gas Industry, 2013, 33(12): 173-178.

    [5]王英梅, 吳青柏, 蒲毅彬, 展靜. 溫度梯度對(duì)粗砂中甲烷水合物形成和分解過程的影響及電阻率響應(yīng)[J]. 天然氣地球科學(xué), 2012, 23(1): 19-25. Wang Yingmei, Wu Qingbai, Pu Yibin, Zhan Jing. Effect of temperature gradient on process of methane hydrate formationdissociation and its resistivity changes in coarse sand[J]. Natural Gas Geoscience, 2012, 23(1): 19-25.

    [6]趙洪偉, 刁少波, 業(yè)渝光, 張劍, 劉昌嶺. 多孔介質(zhì)中水合物阻抗探測(cè)技術(shù)[J]. 海洋地質(zhì)與第四紀(jì)地質(zhì), 2005, 25(1): 137-142. Zhao Hongwei, Diao Shaobo, Ye Yuguang, Zhang Jian, Liu Changling. Technique of detecting impedance of hydrate in porous medium[J]. Marine Geology & Quaternary Geology, 2005, 25(1): 137-142.

    [7]Riedel M, Collett TS, Hyndman RD. Gas hydrate concentration estimates from chlorinity, electrical resistivity and seismic velocity [J]. Geological Survey of Canada, 2005, 4934: 28-32.

    [8]Buffett BA, Zatsepina OY. Formation of gas hydrate from dissolved gas in natural porous media[J]. Marine Geology, 2000, 164(1/2): 69-77.

    [9]Spangenberg E, Kulenkampff J. Influe nce of methane hydrate content on electrical sediment properties[J]. Geophysical Research Letters, 2006, 33(24): 1-5.

    [10]Spangenberg E, Kulenkampff J, Naumann R, Erzinger J. Pore space hydrate formation in a glass bead sample from methane dissolved in water[J]. Geophysical Research Letters, 2005, 32(24): 1-4.

    [11]Spangenberg E. Modeling of the influence of gas hydrate content on the electrical properties of porous sediments[J]. Journal of Geophysical Research: Solid Earth, 2001, 106(B4): 6535-6548.

    [12]陳強(qiáng), 業(yè)渝光, 孟慶國, 劉昌嶺. 多孔介質(zhì)中CO2水合物飽和度與阻抗關(guān)系模擬實(shí)驗(yàn)研究[J]. 天然氣地球科學(xué), 2009, 20(2): 249-253. Chen Qiang, Ye Yuguang, Meng Qingguo, Liu Changling. Simulation experiment of the relationship between CO2hydrate saturation and resistance in porous media[J]. Natural Gas Geoscience, 2009, 20(2): 249-253.

    [13]陳玉鳳, 李棟梁, 梁德青, 周雪冰, 吳能友. 南海沉積物天然氣水合物飽和度與電阻率的關(guān)系[J]. 石油學(xué)報(bào), 2013, 34(3): 507-512. Chen Yufeng, Li Dongliang, Liang Deqing, Zhou Xuebing, Wu Nengyou. Relationship between gas hydrate saturation and resistivity in sediments of the South China Sea[J]. Acta Petrolei Sinica, 2013, 34(3): 507-512.

    [14]邢蘭昌, 陳強(qiáng), 劉昌嶺. 基于電化學(xué)阻抗譜測(cè)試方法研究四氫呋喃水合物的生成和分解過程[J]. 巖礦測(cè)試, 2015, 34(6): 704-711. Xing Lanchang, Chen Qiang, Liu Changling. Investigation of formation and dissociation processes of tetrahydrofuran hydrate based on electrochemical impedance spectroscopy[J]. Rock and Mineral Analysis, 2015, 34(6): 704-711.

    [15]Du Frane WL, Stern LA, Weitemeyer KA, Constable S, Pinkston JC, Roberts JJ. Electrical properties of polycrystalline methane hydrate[J]. Geophysical Research Letters, 2011, 38(9): 1-5.

    [16]賈錚, 戴長(zhǎng)松, 陳玲. 電化學(xué)測(cè)量方法[M]. 北京:化學(xué)工業(yè)出版社, 2006. Jia Zheng, Dai Changsong, Chen Ling. The electrochemical measuring methods[M]. Beijing: Chemical Industry Press, 2006.

    [17]劉昌嶺, 業(yè)渝光, 孟慶國. 顯微激光拉曼光譜測(cè)定甲烷水合物的水合指數(shù)[J]. 光譜學(xué)與光譜分析, 2010, 30(4): 963-966. Liu Changling, Ye Yuguang, Meng Qingguo. Determination of hydration number of methane hydrates using micro-laser Raman spectroscopy[J]. Spectroscopy and Spectral Analysis, 2010, 30(4): 963-966.

    [18]胡高偉, 業(yè)渝光, 張劍, 刁少波. 沉積物中天然氣水合物微觀分布模式及其聲學(xué)響應(yīng)特征[J]. 天然氣工業(yè), 2010, 30(3): 120-124. Hu Gaowei, Ye Yuguang, Zhang Jian, Diao Shaobo. Micro-models of gas hydrate and their impact on the acoustic properties of the host sediments[J]. Natural Gas Industry, 2010, 30(3): 120-124.

    [19]杜燕, 何世輝, 黃沖, 馮自平. 多孔介質(zhì)中水合物生成與分解二維實(shí)驗(yàn)研究[J]. 化工學(xué)報(bào), 2008, 59(3): 673-680. Du Yan, He Shihui, Huang Chong, Feng Ziping. Experimental studies of natural gas hydrate formation and dissociation in porous media with 2D experimental system[J]. Journal of Chemical Industry and Engineering, 2008, 59(3): 673-680.

    [20]王飛, 潘保芝. 天然氣水合物儲(chǔ)層的地球物理特性及測(cè)井評(píng)價(jià)[J]. 油氣地球物理, 2010, 8(3): 19-24. Wang Fei, Pan Baozhi. Gas hydrate reservoir geophysical characteristics and logging evaluation[J]. Petroleum Geophysics, 2010, 8(3): 19-24.

    [21]Cook AE, Anderson BI, Rasmus J, Sun K, Li Qiming, Collett TS, et al. Electr ical anisotropy of gas hydrate-bearing sand reservoirs in the Gulf of Mexico[J]. Marine and Petroleum Geology, 2012, 34(1): 72-84.

    [22]肖占山, 徐世浙, 羅延鐘, 王東, 朱世和. 泥質(zhì)砂巖復(fù)電阻率的頻散特性實(shí)驗(yàn)[J]. 高校地質(zhì)學(xué)報(bào), 2006, 12(1): 123-130. Xiao Zhanshan, Xu Shizhe, Luo Yanzhong, Wang Dong, Zhu Shihe. Complex resistivity dispersion properties experiment of shaly sands[J]. Geological Journal of China Universities, 2006, 12(1): 123-130.

    (修改回稿日期2015-11-09編 輯羅冬梅)

    Response characteristics of AC impedance spectroscopy during the formation and dissociation of methane hydrate in porous media

    Jin Xuebin1,2, Chen Qiang2,3, Xing Lanchang1, Liu Changling2,3, Zheng Jinwu1
    (1. College of Information and Control Engineering, China University of Petroleum, Qingdao, Shandong 266580, China; 2. Key Laboratory of Natural Gas Hydrate, Ministry of Land and Resources, Qingdao Institute of Marine Geology, Qingdao, Shandong 266071, China; 3. Laboratory for Marine Mineral Resources & Qingdao National Laboratory for Marine Science and Technology, Qingdao, Shandong 266061, China)
    NATUR. GAS IND. VOLUME 36, ISSUE 3, pp.120-127, 3/25/2016. (ISSN 1000-0976; In Chinese)

    Electrical property test is an effective technique to identify the formation and dissociation process of natural gas hydrate (“the hydrate” for short) and to investigate its dynamic characteristics. The AC impedance spectroscopy testing method is better than the simple ones which only test resistivity parameters, for it can be used to test the variation rules of impedance (e.g. resistivity and reactance) of the tested system in a certain range of frequency. And thus, the electrical characteristics of hydrate can be understood further. In this paper, a system used to test impedance spectroscopy of methane hydrate was designed. Methane hydrate-porous media-brine of 3.5% was prepared as the object in this study. Then, its impedance spectrum was tested during the formation and dissociation of hydrate in porous media and the variation rules of characteristic parameters of impedance spectroscopy with hydrate saturation were analyzed. It is shown that the test results are significantly different when the frequency of excitation source is different. When the test frequency is low (0.1-10.6 Hz), there is no obvious relationship between the hydrate saturation and the impedance spectroscopy parameters of tested system. At middle-high test frequency (101.8 Hz-1 MHz), however, the impedance spectroscopy parameters rise with the increase of hydrate saturation. Besides, when the test frequency is 101.8, 1 114 and 1 0746 Hz, there is a monotonic increasing linear relationship between the hydrate saturation and the impedance spectroscopy parameters of the tested system within a certain range of allowable errors. And finally, there is an obvious correspondence between hydrate saturation and the frequency dispersion of impedance amplitude.

    Natural gas hydrate; Electrical property; Formation and dissociation process; AC impedance spectroscopy; Testing method; Porous media; Saturation of hydrate; Dispersion characteristics

    10.3787/j.issn.1000-0976.2016.03.017

    國土資源部公益性行業(yè)科研專項(xiàng)基金“天然氣水合物實(shí)驗(yàn)測(cè)試技術(shù)研究”(編號(hào):201111026)、國土資源部天然氣水合物重點(diǎn)實(shí)驗(yàn)室開放基金重點(diǎn)項(xiàng)目“含天然氣水合物多孔介質(zhì)的電性參數(shù)頻散特性及飽和度計(jì)算方法研究”(編號(hào):SHW[2014]-ZD-01)。

    金學(xué)彬,1990年生,碩士研究生;從事天然氣水合物電學(xué)測(cè)試研究工作。地址:(266580)山東省青島市黃島區(qū)長(zhǎng)江西路66號(hào)。ORCID:0000-0002-9588-6953。E-mail:jinxbupc@163.com

    猜你喜歡
    水合物甲烷飽和度
    糖臬之吻
    氣井用水合物自生熱解堵劑解堵效果數(shù)值模擬
    液氧甲烷發(fā)動(dòng)機(jī)
    論煤炭運(yùn)輸之甲烷爆炸
    水上消防(2020年1期)2020-07-24 09:26:02
    熱水吞吐開采水合物藏?cái)?shù)值模擬研究
    Gas from human waste
    天然氣水合物保壓轉(zhuǎn)移的壓力特性
    我國海域天然氣水合物試采成功
    制作一個(gè)泥土飽和度測(cè)試儀
    巧用有機(jī)物的不飽和度
    久久久久网色| 91午夜精品亚洲一区二区三区| av福利片在线观看| 日韩国内少妇激情av| 性高湖久久久久久久久免费观看| 一级二级三级毛片免费看| 国产精品一及| 国产探花极品一区二区| 日日摸夜夜添夜夜添av毛片| 亚洲国产成人一精品久久久| 老司机影院成人| 国产黄频视频在线观看| 国产精品蜜桃在线观看| 少妇 在线观看| 中文欧美无线码| 80岁老熟妇乱子伦牲交| 国产成人精品一,二区| 性色av一级| 午夜精品国产一区二区电影| 久久国产乱子免费精品| 成人国产av品久久久| 特大巨黑吊av在线直播| 王馨瑶露胸无遮挡在线观看| 亚洲人与动物交配视频| 狂野欧美激情性bbbbbb| 高清毛片免费看| 五月天丁香电影| av女优亚洲男人天堂| 国产精品一二三区在线看| 国产精品人妻久久久久久| 人妻少妇偷人精品九色| 久久99热这里只频精品6学生| 成人无遮挡网站| 最后的刺客免费高清国语| av不卡在线播放| 亚洲精品国产av成人精品| 国产精品一区www在线观看| 干丝袜人妻中文字幕| 久久精品国产亚洲av天美| 在线观看一区二区三区激情| 亚洲欧美成人精品一区二区| 国产一区二区三区综合在线观看 | 纵有疾风起免费观看全集完整版| 在线观看av片永久免费下载| 欧美亚洲 丝袜 人妻 在线| 久久久久久九九精品二区国产| 午夜福利在线在线| 欧美日韩精品成人综合77777| 亚洲一区二区三区欧美精品| 亚洲精品456在线播放app| 五月天丁香电影| 美女内射精品一级片tv| 国产精品熟女久久久久浪| 少妇 在线观看| 亚洲国产毛片av蜜桃av| 亚洲精品一二三| 亚洲最大成人中文| 成年美女黄网站色视频大全免费 | av国产精品久久久久影院| 97热精品久久久久久| 日韩亚洲欧美综合| 久久人人爽人人片av| 免费黄网站久久成人精品| 中文字幕亚洲精品专区| 久久人人爽人人片av| av一本久久久久| 久久 成人 亚洲| 欧美成人a在线观看| 久久久国产一区二区| 自拍欧美九色日韩亚洲蝌蚪91 | 色视频在线一区二区三区| 亚洲精品乱码久久久久久按摩| 亚洲图色成人| 国产黄色视频一区二区在线观看| 亚洲色图综合在线观看| 日韩人妻高清精品专区| 国内精品宾馆在线| 老师上课跳d突然被开到最大视频| 成人18禁高潮啪啪吃奶动态图 | 欧美3d第一页| 中文资源天堂在线| 麻豆精品久久久久久蜜桃| 少妇人妻一区二区三区视频| 国产高潮美女av| 韩国av在线不卡| 免费观看的影片在线观看| 精品人妻熟女av久视频| 中文字幕人妻熟人妻熟丝袜美| 91狼人影院| 日韩人妻高清精品专区| 中国三级夫妇交换| 久久人人爽人人爽人人片va| 精品一区在线观看国产| 成人影院久久| 国语对白做爰xxxⅹ性视频网站| 国产亚洲一区二区精品| 丝瓜视频免费看黄片| 日韩在线高清观看一区二区三区| 秋霞伦理黄片| av黄色大香蕉| 亚洲中文av在线| 一区二区三区精品91| 一级毛片 在线播放| 亚洲精品国产色婷婷电影| 精品午夜福利在线看| 国产精品女同一区二区软件| 噜噜噜噜噜久久久久久91| 亚洲精品久久午夜乱码| 自拍欧美九色日韩亚洲蝌蚪91 | 蜜桃在线观看..| 国产大屁股一区二区在线视频| 国产中年淑女户外野战色| 一级毛片电影观看| 亚洲精品自拍成人| 丝袜喷水一区| 久久久久网色| 亚洲成人中文字幕在线播放| .国产精品久久| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品456在线播放app| 青春草国产在线视频| 校园人妻丝袜中文字幕| 精品久久久久久电影网| 亚洲自偷自拍三级| 亚洲高清免费不卡视频| 又粗又硬又长又爽又黄的视频| 精品国产露脸久久av麻豆| 一个人看的www免费观看视频| 日韩中字成人| 成人国产av品久久久| 干丝袜人妻中文字幕| 亚洲一级一片aⅴ在线观看| 一区二区三区精品91| 国产真实伦视频高清在线观看| 欧美精品一区二区大全| 欧美精品人与动牲交sv欧美| 欧美日韩在线观看h| 一本—道久久a久久精品蜜桃钙片| xxx大片免费视频| 国产v大片淫在线免费观看| 久久久久久久精品精品| freevideosex欧美| 91久久精品电影网| 亚洲色图av天堂| 久久毛片免费看一区二区三区| 只有这里有精品99| 肉色欧美久久久久久久蜜桃| 最近手机中文字幕大全| 热re99久久精品国产66热6| 国产在线视频一区二区| 亚洲精品日韩av片在线观看| 青青草视频在线视频观看| 欧美3d第一页| 国产伦理片在线播放av一区| 国产精品人妻久久久久久| 国产精品偷伦视频观看了| 亚洲国产精品一区三区| 久久精品国产亚洲网站| 国产欧美日韩精品一区二区| 亚洲三级黄色毛片| 中文资源天堂在线| 国产精品人妻久久久久久| 亚洲不卡免费看| 国产高清三级在线| 国产伦理片在线播放av一区| 国产精品一区二区三区四区免费观看| 黑人猛操日本美女一级片| av在线播放精品| 韩国av在线不卡| 国产大屁股一区二区在线视频| 国产探花极品一区二区| 街头女战士在线观看网站| 精品国产露脸久久av麻豆| 免费av不卡在线播放| 日日撸夜夜添| 免费看av在线观看网站| 十八禁网站网址无遮挡 | 最新中文字幕久久久久| 中文资源天堂在线| 国产高潮美女av| 国产成人精品婷婷| 天天躁日日操中文字幕| 日韩欧美精品免费久久| 国产欧美日韩一区二区三区在线 | 久久久亚洲精品成人影院| 久久久久久伊人网av| 国产伦精品一区二区三区视频9| 久久久久久伊人网av| 日本av手机在线免费观看| 欧美性感艳星| 干丝袜人妻中文字幕| 另类亚洲欧美激情| 人妻制服诱惑在线中文字幕| 我的老师免费观看完整版| 久久精品人妻少妇| av线在线观看网站| 一级二级三级毛片免费看| 内地一区二区视频在线| 欧美日本视频| 欧美激情极品国产一区二区三区 | 久久精品国产亚洲av涩爱| 黄色一级大片看看| 免费观看在线日韩| 建设人人有责人人尽责人人享有的 | 国产熟女欧美一区二区| videossex国产| 久久 成人 亚洲| 男女下面进入的视频免费午夜| 国产欧美另类精品又又久久亚洲欧美| 天美传媒精品一区二区| 久久精品久久精品一区二区三区| 国产中年淑女户外野战色| 丝袜喷水一区| 中文字幕制服av| 亚洲av综合色区一区| 欧美高清性xxxxhd video| 在线 av 中文字幕| 毛片一级片免费看久久久久| 欧美国产精品一级二级三级 | 精品酒店卫生间| 国产 一区精品| 爱豆传媒免费全集在线观看| 国产亚洲5aaaaa淫片| 在线观看三级黄色| 国产精品久久久久成人av| 蜜桃亚洲精品一区二区三区| 亚洲欧美日韩东京热| 亚洲精品久久久久久婷婷小说| 精品一区在线观看国产| 亚洲国产色片| 久久久久久久国产电影| 美女中出高潮动态图| 制服丝袜香蕉在线| 亚洲av福利一区| 嫩草影院入口| 成人二区视频| 亚洲精品aⅴ在线观看| 久久久久视频综合| 狂野欧美激情性bbbbbb| 一级毛片黄色毛片免费观看视频| 国产高清有码在线观看视频| 另类亚洲欧美激情| 一本色道久久久久久精品综合| 麻豆精品久久久久久蜜桃| 久久国产精品大桥未久av | 人妻少妇偷人精品九色| 色网站视频免费| 久久久精品94久久精品| 久久精品国产a三级三级三级| 亚洲精品视频女| 国模一区二区三区四区视频| 亚洲在久久综合| 国产亚洲91精品色在线| 97在线人人人人妻| 亚洲精品,欧美精品| 国产在视频线精品| 色综合色国产| 深爱激情五月婷婷| 国产大屁股一区二区在线视频| 免费人妻精品一区二区三区视频| 婷婷色综合www| 内射极品少妇av片p| 国产男人的电影天堂91| 成年av动漫网址| 精品人妻偷拍中文字幕| av国产久精品久网站免费入址| 中国美白少妇内射xxxbb| 18禁裸乳无遮挡动漫免费视频| 妹子高潮喷水视频| 午夜日本视频在线| 在线 av 中文字幕| 全区人妻精品视频| 六月丁香七月| 人妻一区二区av| av又黄又爽大尺度在线免费看| 一级av片app| 边亲边吃奶的免费视频| 亚洲va在线va天堂va国产| 高清视频免费观看一区二区| 多毛熟女@视频| 又爽又黄a免费视频| 只有这里有精品99| 欧美精品亚洲一区二区| 欧美xxⅹ黑人| 日韩免费高清中文字幕av| 久久精品国产亚洲av涩爱| 一级a做视频免费观看| 欧美日韩亚洲高清精品| 你懂的网址亚洲精品在线观看| 久久久久久久精品精品| 国产精品国产av在线观看| 我的女老师完整版在线观看| 老熟女久久久| 美女主播在线视频| 日本黄色片子视频| 亚洲av中文av极速乱| 国产亚洲精品久久久com| 色综合色国产| 人人妻人人添人人爽欧美一区卜 | 国产成人免费观看mmmm| 蜜桃亚洲精品一区二区三区| 国产永久视频网站| 我要看黄色一级片免费的| 免费大片黄手机在线观看| 国产视频首页在线观看| 精品视频人人做人人爽| 日韩中文字幕视频在线看片 | 五月玫瑰六月丁香| 老师上课跳d突然被开到最大视频| 精华霜和精华液先用哪个| 国产人妻一区二区三区在| 老司机影院毛片| 日韩成人av中文字幕在线观看| 人妻 亚洲 视频| 亚洲怡红院男人天堂| a级毛色黄片| 欧美精品一区二区大全| 精品一区二区三区视频在线| 美女脱内裤让男人舔精品视频| 国产高清有码在线观看视频| 国产探花极品一区二区| 综合色丁香网| 99久久精品一区二区三区| 一级毛片 在线播放| 国产91av在线免费观看| 国产一区二区在线观看日韩| 晚上一个人看的免费电影| 国产亚洲5aaaaa淫片| 国产精品一区www在线观看| 国产精品久久久久久久电影| 能在线免费看毛片的网站| 国产精品一区二区三区四区免费观看| 天美传媒精品一区二区| 亚洲久久久国产精品| 黑人猛操日本美女一级片| 亚洲第一av免费看| 在线观看三级黄色| 国产色婷婷99| 国产国拍精品亚洲av在线观看| 久久精品国产a三级三级三级| 免费观看a级毛片全部| 国产黄色视频一区二区在线观看| 欧美精品人与动牲交sv欧美| 一本一本综合久久| 日本免费在线观看一区| 在线观看免费高清a一片| kizo精华| 一级二级三级毛片免费看| 成人午夜精彩视频在线观看| 日韩大片免费观看网站| 婷婷色av中文字幕| 国产精品一区二区在线观看99| 丰满人妻一区二区三区视频av| 国产成人一区二区在线| 黄色日韩在线| 日本午夜av视频| 美女xxoo啪啪120秒动态图| 在线亚洲精品国产二区图片欧美 | 国产黄片视频在线免费观看| 国产色爽女视频免费观看| 国产精品国产三级专区第一集| 啦啦啦视频在线资源免费观看| tube8黄色片| 男人舔奶头视频| 91久久精品国产一区二区成人| 精品人妻一区二区三区麻豆| 亚洲精品亚洲一区二区| 卡戴珊不雅视频在线播放| 国产淫语在线视频| 精华霜和精华液先用哪个| 2021少妇久久久久久久久久久| 在线精品无人区一区二区三 | 老女人水多毛片| freevideosex欧美| 六月丁香七月| 国产精品国产av在线观看| 亚洲美女搞黄在线观看| 99热全是精品| 美女中出高潮动态图| 精品久久久久久电影网| 国产人妻一区二区三区在| 国产毛片在线视频| 久久99精品国语久久久| 丰满迷人的少妇在线观看| 国内精品宾馆在线| 国产高清三级在线| 亚洲精品久久久久久婷婷小说| 久久99热这里只频精品6学生| 欧美极品一区二区三区四区| 深夜a级毛片| 亚洲国产欧美人成| 2021少妇久久久久久久久久久| 亚洲精品一二三| 精品久久久久久电影网| 国产中年淑女户外野战色| 欧美日韩综合久久久久久| 国产亚洲午夜精品一区二区久久| 久久精品国产鲁丝片午夜精品| 国产黄色免费在线视频| 午夜福利影视在线免费观看| 中国美白少妇内射xxxbb| 日韩中文字幕视频在线看片 | 少妇高潮的动态图| 久久久精品94久久精品| 直男gayav资源| 日日啪夜夜撸| 亚洲四区av| 三级经典国产精品| 一级二级三级毛片免费看| 日韩视频在线欧美| 视频区图区小说| 一级黄片播放器| 插逼视频在线观看| 久久久久人妻精品一区果冻| 国产一区有黄有色的免费视频| 3wmmmm亚洲av在线观看| 国产成人aa在线观看| 精品国产露脸久久av麻豆| 久久青草综合色| 久久精品夜色国产| 免费黄频网站在线观看国产| 这个男人来自地球电影免费观看 | 三级国产精品片| 熟女电影av网| 五月玫瑰六月丁香| 欧美日韩亚洲高清精品| 国产极品天堂在线| 插阴视频在线观看视频| 人人妻人人澡人人爽人人夜夜| 晚上一个人看的免费电影| 国产精品久久久久久精品电影小说 | 国产女主播在线喷水免费视频网站| 欧美成人a在线观看| 黄色欧美视频在线观看| 在线观看av片永久免费下载| 国产av国产精品国产| 大又大粗又爽又黄少妇毛片口| 精品国产三级普通话版| 建设人人有责人人尽责人人享有的 | 精品国产一区二区三区久久久樱花 | 成人漫画全彩无遮挡| 男女无遮挡免费网站观看| 亚洲av.av天堂| 国产无遮挡羞羞视频在线观看| 两个人的视频大全免费| 大片免费播放器 马上看| 亚洲,欧美,日韩| 亚洲内射少妇av| 老女人水多毛片| 成人国产麻豆网| 韩国高清视频一区二区三区| 亚洲av欧美aⅴ国产| 免费观看的影片在线观看| 国产精品.久久久| 99久久精品一区二区三区| 亚洲欧美成人精品一区二区| 高清不卡的av网站| 国产熟女欧美一区二区| 乱系列少妇在线播放| av在线观看视频网站免费| 极品教师在线视频| 国产午夜精品久久久久久一区二区三区| 国语对白做爰xxxⅹ性视频网站| 国产色婷婷99| 国产欧美日韩精品一区二区| 3wmmmm亚洲av在线观看| 国产美女午夜福利| .国产精品久久| 久久人人爽av亚洲精品天堂 | 大片免费播放器 马上看| 高清在线视频一区二区三区| 欧美三级亚洲精品| 国产av国产精品国产| 中文字幕久久专区| 最新中文字幕久久久久| 99热全是精品| 男的添女的下面高潮视频| 亚洲欧美日韩卡通动漫| 国产伦在线观看视频一区| 国产精品免费大片| 亚洲av成人精品一区久久| 99久久精品热视频| 嘟嘟电影网在线观看| 亚洲激情五月婷婷啪啪| 最近手机中文字幕大全| 成人免费观看视频高清| 午夜免费鲁丝| 亚洲无线观看免费| 国产精品爽爽va在线观看网站| 成年女人在线观看亚洲视频| av免费在线看不卡| 亚洲图色成人| 久久精品国产自在天天线| 婷婷色av中文字幕| 日产精品乱码卡一卡2卡三| 看非洲黑人一级黄片| 精品亚洲乱码少妇综合久久| 亚洲国产欧美人成| 亚洲最大成人中文| 99久久综合免费| 1000部很黄的大片| 亚洲精品中文字幕在线视频 | 我的女老师完整版在线观看| 午夜视频国产福利| 国产成人免费无遮挡视频| 永久网站在线| 麻豆精品久久久久久蜜桃| 亚洲美女视频黄频| 成年女人在线观看亚洲视频| 久久精品久久精品一区二区三区| 老师上课跳d突然被开到最大视频| 久久人妻熟女aⅴ| 秋霞在线观看毛片| 草草在线视频免费看| 国产色爽女视频免费观看| 狂野欧美激情性xxxx在线观看| 伦理电影大哥的女人| 在线 av 中文字幕| 国产精品福利在线免费观看| 国产大屁股一区二区在线视频| 91久久精品国产一区二区三区| 色综合色国产| a 毛片基地| 国产一区二区三区av在线| 新久久久久国产一级毛片| 丰满乱子伦码专区| 伦理电影大哥的女人| 亚洲欧美中文字幕日韩二区| 成人黄色视频免费在线看| 少妇精品久久久久久久| 免费大片18禁| 欧美日韩一区二区视频在线观看视频在线| 18禁裸乳无遮挡免费网站照片| 久久精品人妻少妇| 啦啦啦啦在线视频资源| 高清日韩中文字幕在线| 80岁老熟妇乱子伦牲交| 午夜福利视频精品| 日韩av不卡免费在线播放| 在线亚洲精品国产二区图片欧美 | 看非洲黑人一级黄片| 精品一区二区三卡| 亚洲欧美日韩卡通动漫| 交换朋友夫妻互换小说| 中文字幕免费在线视频6| 亚洲av福利一区| 国产在线视频一区二区| 成年人午夜在线观看视频| 婷婷色综合大香蕉| 亚洲第一av免费看| 91在线精品国自产拍蜜月| 国产精品一区二区性色av| 亚洲精品456在线播放app| 久久久久视频综合| 99re6热这里在线精品视频| 毛片女人毛片| 汤姆久久久久久久影院中文字幕| 又爽又黄a免费视频| 舔av片在线| 一级毛片黄色毛片免费观看视频| 美女高潮的动态| 爱豆传媒免费全集在线观看| 一个人看视频在线观看www免费| 一本久久精品| 国产真实伦视频高清在线观看| 老女人水多毛片| 青青草视频在线视频观看| 永久网站在线| 中文字幕亚洲精品专区| 欧美国产精品一级二级三级 | 国产精品秋霞免费鲁丝片| 国产午夜精品一二区理论片| 99久久综合免费| 2022亚洲国产成人精品| 在线观看人妻少妇| 日韩成人av中文字幕在线观看| 大话2 男鬼变身卡| 午夜精品国产一区二区电影| 黑人猛操日本美女一级片| 国产男女超爽视频在线观看| 人妻系列 视频| 精品久久久久久久久亚洲| 国产大屁股一区二区在线视频| 日韩av免费高清视频| 日日摸夜夜添夜夜添av毛片| 午夜免费鲁丝| 免费看av在线观看网站| 亚洲无线观看免费| 国产黄片视频在线免费观看| 欧美激情极品国产一区二区三区 | 青青草视频在线视频观看| 免费看光身美女| 午夜福利视频精品| 你懂的网址亚洲精品在线观看| 亚洲国产精品成人久久小说| 这个男人来自地球电影免费观看 | 国产精品一区二区三区四区免费观看| 18禁在线无遮挡免费观看视频| 免费av不卡在线播放| 女性被躁到高潮视频| 亚洲精品自拍成人| 青春草国产在线视频| 日本午夜av视频| 久久精品国产a三级三级三级| 中国三级夫妇交换| 中文字幕久久专区| 成人影院久久| 十八禁网站网址无遮挡 | 日本一二三区视频观看| 欧美三级亚洲精品| 精品人妻熟女av久视频| 2022亚洲国产成人精品| 久久久久久人妻| 激情 狠狠 欧美| 国产精品一区二区性色av|