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

    拉曼光譜在油氣地質(zhì)應(yīng)用中的研究進(jìn)展

    2022-09-05 03:26:44柳少波魯雪松范俊佳馬行陟桂麗黎
    光譜學(xué)與光譜分析 2022年9期
    關(guān)鍵詞:特征研究

    劉 強(qiáng), 柳少波, 2, 魯雪松, 2*, 范俊佳, 2, 田 華, 2, 馬行陟, 2, 桂麗黎, 2

    1. 中國石油勘探開發(fā)研究院, 北京 100083 2. 中國石油天然氣集團(tuán)公司盆地構(gòu)造與油氣成藏重點(diǎn)實(shí)驗(yàn)室, 北京 100083

    引 言

    1928年, 印度科學(xué)家Raman發(fā)現(xiàn)拉曼效應(yīng)[1], 但直至20世紀(jì)70年代左右, 拉曼光譜才逐漸應(yīng)用到地質(zhì)研究中。 1969年, Griffith[2]介紹拉曼光譜在礦物研究中的應(yīng)用[2]。 1970年, Tuinstra與Koening[3]發(fā)現(xiàn)石墨拉曼光譜中存在無序峰(Disordered band, D峰)、 石墨峰(Graphite Band, G峰)。 1975年, Rosasco等[4]應(yīng)用拉曼光譜分析流體包裹體。 隨著儀器、 技術(shù)的進(jìn)步, 拉曼光譜大量應(yīng)用到地質(zhì)研究中。 國內(nèi)外學(xué)者利用拉曼光譜結(jié)合冷熱臺、 合成包裹體、 水熱金剛石壓腔(hydrothermal diamond anvil cell, HDAC)、 熔融毛細(xì)硅管(fused silica capillary capsules, FSCCs)、 高壓可視反應(yīng)腔(high-pressure optical cell, HPOC)、 掃描電鏡(scanning electron microscope, SEM)等儀器或技術(shù), 在不同溫壓條件下, 對礦物成分與結(jié)構(gòu)[5-6]、 有機(jī)質(zhì)成熟度[7-21]、 氣體組分[22-23]、 烴類[24-27]、 氣體密度[23, 28-40]、 溶液中陰陽離子類型及濃度[41-54]、 氣體溶解度與擴(kuò)散系數(shù)[55-66]、 天然氣水合物[67-68]、 同位素相對豐度[69-73]等開展大量研究工作, 并建立了相應(yīng)的定量/半定量分析方法與圖版。

    前人已發(fā)表較多關(guān)于拉曼光譜在礦物學(xué)研究[74]、 流體包裹體分析[1, 75-78]、 壓力標(biāo)定[79-81]、 有機(jī)質(zhì)熱成熟度測試[9, 82]、 天然氣水合物研究[68]、 氣體同位素分析[83]等方面的綜述性文章, 并對樣品與儀器差異性、 熒光、 拉曼散射截面參數(shù)、 標(biāo)準(zhǔn)物質(zhì)光譜等問題都進(jìn)行較為詳細(xì)討論[1, 78]。 但前人綜述中論述對象相對單一, 部分研究成果并不適用于實(shí)際油氣地質(zhì)研究, 同時(shí)大量研究成果并未得到廣泛應(yīng)用。 本文主要從油氣地質(zhì)研究內(nèi)容及需求出發(fā), 結(jié)合拉曼光譜定性、 定量分析原理, 將激光拉曼光譜在油氣地質(zhì)研究中的應(yīng)用分為礦物與流體包裹體分析、 有機(jī)質(zhì)成熟度計(jì)算與流體包裹體壓力恢復(fù)三方面, 并結(jié)合相應(yīng)的應(yīng)用實(shí)例, 指出實(shí)際測試分析過程中樣品、 測試條件、 數(shù)據(jù)處理可能存在的一些問題, 以及進(jìn)一步研究的方向, 以促進(jìn)拉曼光譜進(jìn)一步應(yīng)用到油氣地質(zhì)研究中。

    1 拉曼光譜分析原理

    拉曼光譜基于物質(zhì)對光的非彈性散射效應(yīng): 當(dāng)固定頻率的光照射物質(zhì), 發(fā)生彈性碰撞, 頻率不變?yōu)槿鹄€(Rayleigh); 頻率發(fā)生改變的為非彈性散射(拉曼效應(yīng): Raman), 其中頻率增大(ν1+νm)為反斯托克(anti-Stokes)散射, 頻率減小(νl-νm)為斯托克(Stokes)散射[1]。 拉曼光譜研究常用斯托克散射, 但也有學(xué)者利用反斯托克散射(coherent anti-Stocks Raman scatting, CARS)進(jìn)行研究[84]。

    拉曼光譜定性分析原理:νm表示分子基本的旋轉(zhuǎn)、 振動或晶格頻率, 僅取決于分子本身的固有振動和轉(zhuǎn)動能級結(jié)構(gòu), 與入射光的頻率無關(guān)。 只有分子振動過程中極化率發(fā)生改變的物質(zhì)才具拉曼活性, 單原子離子(Na+、 Cl-等)不具拉曼活性[75]。 由于多原子分子的不同振動也存在差異, 具有不同的拉曼特征峰, 所以, 拉曼光譜可以定性識別具有拉曼活性的物質(zhì)(表1和表2)。

    表1 流體包裹體中常見離子和礦物拉曼特征峰峰位[1]

    物質(zhì)的振動頻率與物質(zhì)的質(zhì)量有關(guān), 可以利用拉曼光譜識別不同的同位素組成。 雙原子分子的振動頻率計(jì)算公式[式(1)][85]

    (1)

    式(1)中,ν為頻率;K為鍵的力常數(shù);m1和m2分別為相對原子質(zhì)量。

    拉曼光譜定量分析原理: 拉曼散射的面積、 強(qiáng)度與分子的濃度、 入射光強(qiáng)度等因素有關(guān)。 根據(jù)Placzek的極化率理論, 拉曼特征峰面積的計(jì)算公式[75]

    (2)

    表2 流體包裹體中常見氣體拉曼特征峰峰位、 相對拉曼散射截面系數(shù)(RRSCS)(Σ和σ)[77, 88]

    拉曼峰強(qiáng)度的計(jì)算公式[式(3)][85]

    (3)

    式(3)中,IR為拉曼散射強(qiáng)度;I0為入射強(qiáng)度;N為特定狀態(tài)下的散射分子數(shù);ν為激光頻率;α為分子的極化率;δ為振動幅度。

    由式(2)和式(3)可以看出, 拉曼散射的面積和強(qiáng)度都與N(分子數(shù))有關(guān), 可以利用拉曼峰強(qiáng)度比或峰面積比來求取特定物質(zhì)的相對含量[75]。

    2 拉曼光譜在礦物與流體包裹體分析中的應(yīng)用

    礦物與流體包裹體分析是拉曼光譜在油氣地質(zhì)研究中的基礎(chǔ)應(yīng)用, 可以確定成巖-成藏流體類型及組成、 成巖-成藏作用過程, 主要是依據(jù)拉曼定量與定性分析原理進(jìn)行研究。 拉曼光譜結(jié)合冷熱臺、 熔融毛細(xì)硅管(FSCCs)、 高壓可視反應(yīng)腔(HPOC)、 人工合成包裹體、 水熱金剛石壓腔(HDAC)、 掃描電鏡(SEM)、 電子探針(electron probe microanalysis, EPMA)等儀器、 技術(shù), 可進(jìn)行實(shí)時(shí)、 原位、 特定溫壓下拉曼光譜測試。 主要應(yīng)用有巖礦鑒定、 流體包裹體組分分析、 水巖相互作用過程、 地質(zhì)樣品同位素研究(表3)。

    表3 拉曼光譜在油氣地質(zhì)研究中主要應(yīng)用

    2.1 巖礦鑒定

    對于通過常規(guī)光學(xué)顯微鏡無法區(qū)分、 粒徑較小的礦物, 如包裹體中的子礦物, 利用拉曼特征峰峰位和強(qiáng)度, 結(jié)合標(biāo)準(zhǔn)礦物圖庫, 可以快速準(zhǔn)確識別不同的礦物類型。 類質(zhì)同象系列礦物、 同質(zhì)多相礦物由于元素組成、 分子結(jié)構(gòu)不同, 可以利用拉曼光譜特征進(jìn)行區(qū)分, 如正長石與鈉長石、 α-石英與β-石英[6, 86]。 利用拉曼光譜進(jìn)行面掃、 三維重構(gòu), 可以識別不同礦物的空間分布與非均質(zhì)性。 例如: 何佳樂等[86]利用激光拉曼大面積掃描快速成像的方法對兩片變質(zhì)巖薄片進(jìn)行礦物分析并定名; Wille等[87]將激光拉曼與陰極發(fā)光、 電子背散衍射(electron backscatter diffraction, EBSD)、 電子探針(EPMA)、 掃描電鏡(SEM)結(jié)合對天然錫石的環(huán)帶中的微量元素及晶體取向進(jìn)行微觀表征。

    2.2 流體包裹體組分分析

    油氣儲層中常見的流體包裹體有純液相包裹體、 氣液兩相包裹體、 富氣體包裹體、 含子礦物包裹體、 含有機(jī)質(zhì)包裹體與烴類包裹體, 常見組成物質(zhì)有輕烴、 發(fā)熒光烴類、 非烴氣體、 水、 鹽類、 子礦物和瀝青[98]。 通過拉曼特征峰峰位與拉曼參數(shù)可以無損的進(jìn)行流體包裹體成分定性識別、 定量/半定量表征, 這些是顯微測溫法、 爆裂法、 激光剝蝕電感耦合等離子質(zhì)譜(laser ablation-inductively coupled plasma-mass spectrometry, LA-ICP-MS)等不具有的優(yōu)勢[1]。

    (1)對于氣液兩相包裹體和單相氣體包裹體中的氣相, 可以通過不同氣體的特征峰面積或者強(qiáng)度比值來判斷氣體組成和相對含量, 從而區(qū)分不同類型及期次流體包裹體[1]。 Biehl等[99]利用拉曼光譜, 計(jì)算鞍狀白云石、 硬石膏、 螢石中流體包裹體的氣體組成, 并以此區(qū)分原生與次生流體包裹體, 并結(jié)合巖相學(xué)、 穩(wěn)定同位素、 稀土元素?cái)?shù)據(jù), 研究熱液白云化作用及硫酸鹽熱還原反應(yīng)(thermochemical sulfate reduction, TSR)對碳酸鹽巖次生孔隙形成過程中的影響。

    (2)富氣體流體包裹體中水的相對含量, 同樣可以通過水峰與氣體特征峰面積或強(qiáng)度比值進(jìn)行計(jì)算[61]。

    圖1 純水拉曼光譜(a)、 NaCl濃度(mol·kg-1 H2O)與特征峰強(qiáng)度比(I3 425/I3 260)關(guān)系(b)(據(jù)Caumon等修改[50])

    (4)烴類包裹體: 張鼐等[24-25]、 陳勇等[26-27]對鏈烷烴、 芳香烴、 環(huán)烷烴、 不飽和烴拉曼光譜特征開展大量研究工作, 結(jié)果表明: 飽和烴拉曼特征峰與含碳數(shù)無關(guān), 只與烴類分子結(jié)構(gòu)和基團(tuán)特征有關(guān), 拉曼光譜只能區(qū)分C1-C4的烷烴, 對于烴類組成復(fù)雜的包裹體只能進(jìn)行定性的判別, 無法準(zhǔn)確判定烴類類型。 含不飽和烴及其衍生物(原油、 瀝青)的油氣包裹體在激光激發(fā)下產(chǎn)生熒光, 會掩蓋拉曼信號, 可以通過紫外拉曼[101]、 紅外拉曼[102]和相干反斯托克散射(CARS)[84]進(jìn)行測試, 在樣品處理過程中可加入熒光抑制劑、 多次累積減少熒光背景的方法進(jìn)行測試[1]。

    (5)含子礦物/瀝青流體包裹體: 子礦物與瀝青, 可以通過特征峰拉曼位移進(jìn)行識別[1]。

    2.3 水-巖相互作用過程

    結(jié)合冷熱臺、 水熱金剛石壓腔(DAC)、 熔融毛細(xì)硅管(FSCCs)和高壓可視反應(yīng)腔(HPOC)可以對特定溫壓下, 對水-巖相互作用過程進(jìn)行原位、 實(shí)時(shí)監(jiān)測, 從而研究化學(xué)反應(yīng)過程、 物質(zhì)的相態(tài)方程、 礦物的形成與轉(zhuǎn)化條件, 可用于油氣地質(zhì)中成巖-成藏演化過程研究。 劉詩琦等[92]利用激光拉曼和水熱金剛石壓腔(HDAC)對碳酸鹽巖-流體反應(yīng)進(jìn)行模擬, 通過ICH3COOH/IH2O和ICO2/IH2O變化情況, 判斷不同溫壓下的碳酸鹽巖溶蝕-沉淀規(guī)律。 Du等[91]利用激光拉曼在南海對天然氣水合物的形成進(jìn)行原位研究, 從而可以判斷天然氣水合物的形成過程和形成時(shí)的溫度、 壓力, 對天然氣水合物的勘探與開發(fā)具有重要意義。

    2.4 物質(zhì)同位素研究

    3 拉曼光譜在有機(jī)質(zhì)成熟度計(jì)算中的應(yīng)用

    沉積有機(jī)質(zhì)在埋藏過程中, 溫度作用會使化學(xué)組成與結(jié)構(gòu)發(fā)生不可逆變化, 雜原子與氫被消耗, 分子結(jié)構(gòu)從無序向有序變化, 這種分子結(jié)構(gòu)變化可以通過有機(jī)質(zhì)拉曼光譜特征進(jìn)行表征[82]。 除埋藏過程中地溫加熱外, 巖漿侵入及構(gòu)造活動所產(chǎn)生的異常熱事件同樣可以使有機(jī)質(zhì)結(jié)構(gòu)發(fā)生相應(yīng)的變化[82]。 因此, 有機(jī)質(zhì)拉曼光譜可以用于沉積盆地中的有機(jī)質(zhì)成熟度測定和有機(jī)質(zhì)地質(zhì)歷史經(jīng)歷最高溫度(Tmax)計(jì)算[8, 12](表3), 從而可以用于恢復(fù)熱史、 油氣成藏期次約束[103]。 拉曼光譜計(jì)算有機(jī)質(zhì)Tmax主要應(yīng)用于構(gòu)造運(yùn)動、 變質(zhì)作用研究, 本文不做過多介紹。

    有機(jī)質(zhì)的拉曼光譜由一級區(qū)域(1 000~1 800 cm-1)和二級區(qū)域(2 400~3 500 cm-1)組成[圖2(a)]。 一級區(qū)域包含兩個(gè)主峰; 無序峰(D峰, 約1 340~1 360 cm-1)和石墨峰(G峰, 約1 580 cm-1)。 除D峰、 G峰外, 在一級區(qū)域擬合過程中, 不同學(xué)者提出一些特征峰, 最常見的為D2峰(約1 610 cm-1)、 D3峰(約1 500 cm-1)、 D4峰(約1 200 cm-1)、 D5(約1 260 cm-1)和D6(約1 440 cm-1); 二級區(qū)域(2 400~3 500 cm-1)中有機(jī)質(zhì)拉曼特征峰是一級區(qū)域特征峰的和頻、 倍頻, 一般在無煙煤階段(VRo>2%)才能夠區(qū)分[82]。

    國內(nèi)外學(xué)者對不同成熟度煤[7, 10-11, 13, 20]、 瀝青[11-12]、 干酪根[7, 15, 17, 104]、 沉積巖中分散有機(jī)質(zhì)[16, 93]和古生物[14, 19]開展了大量拉曼光譜研究, 主要可以得出以下幾點(diǎn)認(rèn)識:

    (1)不同類型、 成熟度有機(jī)質(zhì)樣品, 均能利用拉曼參數(shù)來計(jì)算成熟度。 常用拉曼參數(shù)有: G峰半高寬(full width at half maximum of graphite band, G-FWHM)、 G-D峰間距(Raman band separation, RBS)、 峰強(qiáng)比(R1:ID/IG)、 峰面積比(R2:AD/AG)、 D峰半高寬(D-FWHM)、 半高寬比(D-FWHM/G-FWHM)和D峰或G峰峰位, 多個(gè)拉曼參數(shù)組合RaMM[82]。 不同類型樣品在選取合適參數(shù)時(shí)可能存在差異, Zhou等[12]對瀝青(BRo: 0.43%~4.81%)、 肖賢明等[21]對有機(jī)質(zhì)中的礦物-有機(jī)質(zhì)復(fù)合體(mineral-organic aggregation, MOA) (BRo: 1.71%~4.57%)進(jìn)行拉曼光譜測試, 發(fā)現(xiàn)與成熟度擬合關(guān)系最好的是D峰半高寬(D-FWHM), 其次是D-G峰間距(RBS); Henry等[93]對泥巖中有機(jī)質(zhì)與煤(VRo: 0.5%~6%)進(jìn)行測試, 發(fā)現(xiàn)5個(gè)拉曼參數(shù)中, G峰半高寬(G-FWHM)與VRo等擬合關(guān)系最好,R2為0.96, D-G峰間距(RBS)最差,R2為0.62。 實(shí)際油氣地質(zhì)研究中可以利用RBS計(jì)算成熟度, 同時(shí)用G/D峰半高寬(G/D-FWHM)進(jìn)行約束。

    (2)不同拉曼參數(shù)用于表征有機(jī)質(zhì)成熟度有明顯的分階段現(xiàn)象[圖2(b)], 不同成熟階段需要選取不同的拉曼參數(shù)進(jìn)行表征[11-13, 20, 104]。 Wilkins等[13]和Zhang等[20]分別測試鏡質(zhì)體反射率0.38%~2.53%和0.38%~3.52%的煤, 發(fā)現(xiàn)在鏡質(zhì)體反射率1.2%和1.3%處, 拉曼參數(shù)出現(xiàn)明顯的拐點(diǎn)。 劉德漢等[11]和Zhou等[12]對高過成熟度煤和瀝青拉曼光譜分析中發(fā)現(xiàn), 當(dāng)Ro值<3.5%時(shí), RBS與瀝青成熟度有較好的擬合結(jié)果, 當(dāng)Ro>3.5%時(shí), 特征峰強(qiáng)度比(R1)與瀝青成熟度有較好的擬合結(jié)果。 不同成熟度階段, 拉曼參數(shù)擬合結(jié)果的差異性可能與有機(jī)質(zhì)結(jié)構(gòu)發(fā)生明顯的變化有關(guān)[20], Ro在1.3%~3.52%范圍內(nèi), 烷基發(fā)生斷裂, 多環(huán)芳烴結(jié)構(gòu)不斷增加; 隨著成熟度繼續(xù)增加, 3~5環(huán)芳香烴和無定形碳結(jié)構(gòu)減少[20]。

    圖2 瀝青拉曼光譜(a)、 RBS與不同類型有機(jī)質(zhì)成熟度之間的關(guān)系(b)(數(shù)據(jù)據(jù)文獻(xiàn)[11-12, 19, 21, 104])

    (3)不同學(xué)者數(shù)據(jù)處理方法、 擬合方法不同, 如去基線方式、 擬合公式選取、 擬合峰個(gè)數(shù)、 位置, 導(dǎo)致拉曼參數(shù)與有機(jī)質(zhì)成熟的擬合效果不同, 同時(shí)不同實(shí)驗(yàn)室的測試結(jié)果不具可對比性[16, 93]。 Henry等提出應(yīng)直接讀取有機(jī)質(zhì)拉曼光譜參數(shù), 可以進(jìn)行更快、 更實(shí)際的應(yīng)用, 同時(shí)減少人為誤差[82]。

    4 拉曼光譜在流體包裹體壓力恢復(fù)中的應(yīng)用

    拉曼光譜反映共價(jià)鍵的振動能量, 從而提供關(guān)于鍵強(qiáng)度和分子結(jié)構(gòu)信息, 以及影響共價(jià)鍵的物理(例如: 密度、 分子質(zhì)量)和化學(xué)(例如: 組成)參數(shù)[28], 這是利用拉曼光譜特征進(jìn)行壓力恢復(fù)的基礎(chǔ)。 拉曼光譜進(jìn)行流體包裹體壓力恢復(fù)主要有兩種方法: 氣體特征峰位移法和氣體溶解度法[80](表3)。 前人對利用礦物[5-6]、 鹽水溶液[79, 94-96]、 液態(tài)烴類[97]拉曼位移變化量計(jì)算壓力開展大量研究工作, 但由于測試誤差較大[6, 79], 并不適用于油氣地質(zhì)研究。

    4.1 氣體拉曼特征峰位移計(jì)算壓力

    國內(nèi)外學(xué)者利用冷熱臺、 高壓可視反應(yīng)腔(HPOC)對單一氣體體系(CH4[23, 28, 30-32, 35, 39-40, 105], CO2[29-30, 33-34, 36-39, 106], N2[28, 38-39, 106]和H2[23]等)、 混合氣體體系: 二元體系(CH4-CO2[30, 38, 39], CH4-N2[28, 39], CO2-N2[38-39], CH4-H2[23], CH4-H2O[59], CO2-H2O[33])、 多元體系(CO2-CH4-N2[39], CH4-H2O-NaCl[61], CO2-H2O-NaCl[37], CH4-N2-H2O與CH4-N2-CO2-CO[40])已經(jīng)開展了大量的研究工作, 建立了拉曼參數(shù)-壓力-密度三者之間的關(guān)系式及圖版[圖3(b)], 主要取得以下幾點(diǎn)認(rèn)識:

    (1)CH4的特征峰峰位隨壓力、 密度增加從高頻位置向低頻位置移動[39]; CO2費(fèi)米雙峰(12CO2ν1: 1 285 cm-1左右,12CO22ν2: 1 388 cm-1左右)峰間距隨壓力、 密度增加而增加[39]。

    (2)混合氣體中不同氣體相互影響, 但不同氣體拉曼光譜特征受影響程度存在差異性[28, 30, 38-39]。 相同壓力條件下, CH4-N2體系中N2含量增加, CH4拉曼位移變化量逐漸減小[28, 39]; CH4-CO2體系中CO2含量增加, CO2費(fèi)米雙峰峰間距增加, 而CH4拉曼位移變化不規(guī)則, CO2含量<50%時(shí), CO2含量對CH4位移偏移量影響較小, 可以忽略。 純氣體、 富含氣體包裹體中除了含有不同組分氣體, 同樣含有少量水。 前人[32-33, 37]研究表明, 測試溫度<60 ℃, CO2-H2O體系、 CH4-CO2體系中含有的少量的水對于利用CH4和CO2計(jì)算壓力無影響, 隨溫度增加, 氣體中氣態(tài)水增加導(dǎo)致計(jì)算壓力偏差變大。

    (3)溫度會影響氣體拉曼位移變化。 Zhang等[35]和Le等[39]對不同溫度下CH4拉曼光譜研究結(jié)果顯示: 相同壓力和氣體組成條件下, 隨溫度升高, CH4的ν1拉曼特征峰強(qiáng)度變小、 面積變小、 拉曼位移變大。 Wang等[35]不同溫度下CO2拉曼光譜研究結(jié)果顯示: 相同密度CO2氣體, 隨溫度升高, CO2的費(fèi)米雙峰距間距變小。

    圖3 螢石中單相包裹體拉曼光譜(a)、 CH4特征峰拉曼位移與壓力的關(guān)系(b)(數(shù)據(jù)據(jù)文獻(xiàn)[30-32, 39, 105, 107-108])

    (4)測試數(shù)據(jù)需要進(jìn)行校正。 Lamadrid等[36]系統(tǒng)性采用不同類型儀器、 激光器、 光柵對低密度CO2進(jìn)行拉曼測試, 提出CO2密度計(jì)算之間的差異是由于使用了不同類型的拉曼儀器, 具有不同的配置(激光波長和光柵數(shù)), 以及實(shí)驗(yàn)校準(zhǔn)過程導(dǎo)致。

    4.2 氣體溶解度計(jì)算壓力

    由于室溫下流體包裹體中氣體溶解度相對較小, 同時(shí)可能受氣泡中氣體的影響, 因此, 應(yīng)當(dāng)測試流體包裹體在均一狀態(tài)下的氣體溶解度, 通過氣體-NaCl-H2O體系相態(tài)方程計(jì)算壓力。

    氣體溶解度的測定, 主要利用人工合成流體包裹體[55-56, 61-63]、 高壓可視反應(yīng)腔(HPOC)合成不同的氣體-H2O-NaCl體系, 在不同溫壓下進(jìn)行激光拉曼測試, 通過相態(tài)方程[109-111], 計(jì)算出對應(yīng)測試條件下的氣體溶解度, 從而建立拉曼參數(shù)與氣體溶解度的關(guān)系, 利用這種關(guān)系進(jìn)一步可以計(jì)算氣體在水中的擴(kuò)散系數(shù), 研究天然氣氣藏的破壞與保存情況。 前人已經(jīng)利用這種方法研究了CH4[55-56, 59, 61-63](圖4)、 CO2[58, 64-65, 112]和H2S[66]的溶解度, 主要取得以下幾點(diǎn)認(rèn)識:

    (1)常用于計(jì)算氣體溶解度的拉曼參數(shù)為特征峰強(qiáng)度比、 峰面積; 由于CH4, CO2和H2O具有不同分子結(jié)構(gòu)振動產(chǎn)生的特征峰, 因此不同學(xué)者在選擇不同特征峰強(qiáng)度比、 峰面積比建立的擬合公式存在一定的差異, 如Caumon等[64]在研究CO2-H2O(NaCl)體系時(shí), 建立了ACO2/AH2Ob,ACO2/AH2Os,ACO2-1 380/AH2Ob,ACO2-1 380/AH2Os,ICO2/IH2Ob和ICO2/IH2Os與溶解度之間的關(guān)系式, 不同峰強(qiáng)度、 受鹽度影響不同,ACO2-1 380/AH2Ob在不同鹽溶液中計(jì)算CO2需要進(jìn)行校正, 而ACO2-1 380/AH2Os則可以忽略鹽度的影響。

    (2)不同拉曼參數(shù)受測試方式、 溫壓、 鹽度和數(shù)據(jù)處理等影響不同。 水的相對散射截面參數(shù)隨溫度升高、 鹽度降低而減小, 相同氣體與水的特征峰峰面積比下, 氣體溶解度計(jì)算偏小[62, 113]。 壓力對水峰峰形影響較小, 一般不考慮影響[62, 94]。

    (3)一般應(yīng)當(dāng)選擇氣體典型特征峰(CH4: 2 910 cm-1; CO2: 1 380 cm-1)與水的O—H伸縮振動峰(2 800~3 800cm-1)面積比值, 并帶入相應(yīng)溫度、 壓力下的擬合公式, 計(jì)算氣體溶解度。 在確定氣體組成、 氣體溶解度、 鹽度數(shù)據(jù)、 均一溫度后, 可以利用氣體-H2O-NaCl體系相態(tài)方程準(zhǔn)確計(jì)算流體包裹體均一時(shí)的壓力。

    圖4 不同溫度、 鹽度條件下CH4濃度(mol·kg-1 H2O)與特征峰面積比(PAR)系(據(jù)文獻(xiàn)[55-56, 59, 62]修改)

    5 拉曼光譜在油氣地質(zhì)研究中存在的問題

    5.1 樣品問題

    地質(zhì)樣品經(jīng)歷復(fù)雜的地質(zhì)環(huán)境, 由復(fù)雜的體系構(gòu)成, 與實(shí)驗(yàn)室模擬的單一、 多元體系往往存在一定的差異性, 是制約拉曼光譜在地質(zhì)中應(yīng)用的主要因素。 如: 部分礦物和烴類包裹體在激光拉曼光譜測試過程中產(chǎn)生熒光干擾; 礦物雙折射率[45], 不同類型陰陽離子[51]對水峰峰形的影響不同, 都會導(dǎo)致流體包裹體鹽度計(jì)算和氣體在水中溶解度計(jì)算存在一定誤差; 不同類型有機(jī)質(zhì)、 有機(jī)質(zhì)中不同的顯微組分具有不同的熱物性, 相同地質(zhì)條件下, 熱演化程度不同, 在實(shí)際應(yīng)用過程中需要進(jìn)行換算。

    5.2 測試流程

    拉曼光譜測試過程中主要包括受拉曼光譜校準(zhǔn)、 儀器參數(shù)、 溫壓條件等客觀因素影響。 利用拉曼光譜進(jìn)行定量化的分析, 必須進(jìn)行較準(zhǔn)才能獲取準(zhǔn)確的測試數(shù)據(jù), 如標(biāo)樣、 硅片、 氖燈等校準(zhǔn)。 對鹽水體系、 氣體組分、 氣體密度和有機(jī)質(zhì)成熟度研究發(fā)現(xiàn), 不同的激光波長、 光柵數(shù)測試條件下, 物質(zhì)的拉曼光譜特征不同[18, 36, 50]; 不同溫度和壓力條件下, 物質(zhì)的相對拉曼散射截面參數(shù)(RRSCS)會發(fā)生變化, 如隨著溫度升高, 水的O—H伸縮振動帶的拉曼散射截面(RRSCS)線性減小[110], 這會影響包裹體的鹽度及氣體在水中溶解度計(jì)算。

    5.3 數(shù)據(jù)處理方式

    拉曼光譜數(shù)據(jù)處理主要包括去基線和分峰擬合, 拉曼參數(shù)通常受人為因素影響, 如去基線方式、 擬合峰的個(gè)數(shù)、 位置和強(qiáng)度等[16, 18]。 不同學(xué)者在定量計(jì)算中引入的不同的拉曼參數(shù), 如有機(jī)質(zhì)成熟計(jì)算中G峰半高寬(G-FWHM)、 G-D峰間距(RBS)、 峰高比(R1: Dh/Gh)、 峰面積比(R2、 D/G)、 D峰半高寬(D-FWHM)、 半高寬比(D-FWHM/G-FWHM)、 D峰或G峰峰位。 不同的數(shù)據(jù)處理方式和拉曼參數(shù)選取在一定程度影響實(shí)驗(yàn)數(shù)據(jù)的可信度與可對比度。

    6 未來研究方向

    結(jié)合實(shí)際油氣地質(zhì)研究需要、 目前拉曼光譜研究進(jìn)展與實(shí)際應(yīng)用中存在的問題, 認(rèn)為拉曼光譜研究未來發(fā)展主要圍繞標(biāo)準(zhǔn)化、 復(fù)雜體系與復(fù)雜條件、 新技術(shù)及多種技術(shù)方法聯(lián)用三個(gè)方面展開。

    (1)拉曼光譜測試標(biāo)準(zhǔn)化是前人大量的研究成果的系統(tǒng)性梳理, 從而使這些成果能夠快速應(yīng)用于油氣地質(zhì)研究, 主要包括: 建立標(biāo)準(zhǔn)礦物、 流體的拉曼特征峰數(shù)據(jù)庫; 挑選和制備校準(zhǔn)標(biāo)樣, 用于不同實(shí)驗(yàn)室、 實(shí)驗(yàn)儀器的校準(zhǔn); 規(guī)范拉曼光譜測試樣品前處理方式、 測試條件、 數(shù)據(jù)處理方法, 建立標(biāo)準(zhǔn)的拉曼參數(shù)圖版。

    (2)復(fù)雜體系、 復(fù)雜條件下的拉曼光譜研究是為進(jìn)一步貼近復(fù)雜的地質(zhì)環(huán)境、 地質(zhì)流體體系, 如多元混合氣體組分含量、 含混合氣體組分包裹體壓力恢復(fù)、 混合氣體在復(fù)雜鹽水體系中溶解度; 不同溫壓下, 氣體在鹽水中的擴(kuò)散系數(shù)。 因此, 需要利用水熱金剛石壓腔(HDAC)、 冷熱臺、 熔融毛細(xì)硅管(FSCCs)和高壓可視反應(yīng)腔(HPOC)開展低溫-高溫、 低壓-高壓環(huán)境與復(fù)雜體系(鹽-水體系、 氣-水體系、 氣-鹽-水體系和氣-氣體系等)原位激光拉曼測試。

    (3)紫外激光拉曼、 紅外激光拉曼、 相干反斯托克拉曼散射(CARS)等可以用于研究發(fā)熒光的烴類包裹體、 礦物, 擴(kuò)展拉曼光譜的應(yīng)用范圍。 與掃描電鏡(SEM)、 電子探針(EPMA)等技術(shù)、 方法聯(lián)用, 可以開展更加微觀的研究工作, 同時(shí)與物質(zhì)的化學(xué)組成數(shù)據(jù)相結(jié)合, 也能更好解釋拉曼光譜特征。

    猜你喜歡
    特征研究
    抓住特征巧觀察
    FMS與YBT相關(guān)性的實(shí)證研究
    2020年國內(nèi)翻譯研究述評
    遼代千人邑研究述論
    新型冠狀病毒及其流行病學(xué)特征認(rèn)識
    視錯覺在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    如何表達(dá)“特征”
    不忠誠的四個(gè)特征
    EMA伺服控制系統(tǒng)研究
    抓住特征巧觀察
    18禁在线无遮挡免费观看视频| 一二三四中文在线观看免费高清| 国产女主播在线喷水免费视频网站| 日本黄大片高清| 亚洲欧美精品自产自拍| 日本黄大片高清| 男人舔奶头视频| 日本av免费视频播放| 美女主播在线视频| 在线观看免费日韩欧美大片 | 黑人猛操日本美女一级片| 亚洲一区二区三区欧美精品| 男女啪啪激烈高潮av片| 欧美精品人与动牲交sv欧美| 日韩中字成人| 久久久久久久久久人人人人人人| 观看美女的网站| 黑人高潮一二区| 日韩av不卡免费在线播放| 2022亚洲国产成人精品| 高清午夜精品一区二区三区| 欧美精品人与动牲交sv欧美| av福利片在线| 在线观看三级黄色| 亚洲国产精品成人久久小说| a级一级毛片免费在线观看| 亚洲第一区二区三区不卡| 日韩制服骚丝袜av| 一本大道久久a久久精品| 涩涩av久久男人的天堂| 成人毛片60女人毛片免费| 又大又黄又爽视频免费| av国产久精品久网站免费入址| 性色av一级| 成人综合一区亚洲| 三级国产精品片| 超碰97精品在线观看| 日韩在线高清观看一区二区三区| 久久热精品热| 只有这里有精品99| 视频区图区小说| 观看av在线不卡| 日韩人妻高清精品专区| 三级国产精品欧美在线观看| freevideosex欧美| 亚洲美女黄色视频免费看| 精品久久久久久电影网| 亚洲色图综合在线观看| 伊人久久国产一区二区| 午夜福利在线观看免费完整高清在| 99热网站在线观看| 能在线免费看毛片的网站| 亚洲图色成人| 天堂中文最新版在线下载| 一级毛片 在线播放| 午夜福利视频精品| av在线app专区| 夫妻性生交免费视频一级片| 国产一区亚洲一区在线观看| 国产精品熟女久久久久浪| 在线观看一区二区三区激情| 这个男人来自地球电影免费观看 | 美女脱内裤让男人舔精品视频| 天堂中文最新版在线下载| 国产精品成人在线| 我的女老师完整版在线观看| 欧美另类一区| 国产片特级美女逼逼视频| av女优亚洲男人天堂| 夜夜爽夜夜爽视频| 亚洲av国产av综合av卡| 久久99热这里只频精品6学生| 天堂俺去俺来也www色官网| a级片在线免费高清观看视频| 日韩av不卡免费在线播放| 亚洲精品国产色婷婷电影| 少妇人妻精品综合一区二区| 色5月婷婷丁香| 国产黄频视频在线观看| 久久韩国三级中文字幕| 极品少妇高潮喷水抽搐| 国产免费福利视频在线观看| 精品一区在线观看国产| 一级黄片播放器| 久久av网站| 日本与韩国留学比较| 偷拍熟女少妇极品色| 啦啦啦啦在线视频资源| av福利片在线观看| 亚州av有码| 亚洲av免费高清在线观看| 亚洲av成人精品一区久久| 久久久久国产网址| 亚洲中文av在线| 99九九线精品视频在线观看视频| 丰满少妇做爰视频| 啦啦啦啦在线视频资源| 丁香六月天网| kizo精华| 伦精品一区二区三区| 女性生殖器流出的白浆| 国产精品久久久久久精品古装| 性色av一级| 欧美 日韩 精品 国产| 观看免费一级毛片| 在线观看免费视频网站a站| 女人久久www免费人成看片| 中文天堂在线官网| 精品国产一区二区久久| 七月丁香在线播放| 日韩 亚洲 欧美在线| 人人妻人人看人人澡| 欧美亚洲 丝袜 人妻 在线| 色视频在线一区二区三区| 如何舔出高潮| 日本黄大片高清| 一级毛片电影观看| 三级国产精品片| 人妻制服诱惑在线中文字幕| 久久久久久久久大av| 中文字幕精品免费在线观看视频 | 男人爽女人下面视频在线观看| 国产日韩欧美在线精品| 亚洲国产精品国产精品| 亚洲av男天堂| 国产亚洲午夜精品一区二区久久| 国产日韩一区二区三区精品不卡 | 五月伊人婷婷丁香| 国产男女内射视频| 春色校园在线视频观看| 国产又色又爽无遮挡免| 国产精品秋霞免费鲁丝片| 最近手机中文字幕大全| 中国美白少妇内射xxxbb| 国产 一区精品| 美女福利国产在线| 80岁老熟妇乱子伦牲交| 极品人妻少妇av视频| 国内精品宾馆在线| 国产在线视频一区二区| 777米奇影视久久| 亚洲精品第二区| 黄色配什么色好看| 一边亲一边摸免费视频| 亚洲av不卡在线观看| 国产精品久久久久久久电影| 青春草亚洲视频在线观看| 极品教师在线视频| 亚洲欧美精品自产自拍| 高清不卡的av网站| 国产伦精品一区二区三区四那| 国产真实伦视频高清在线观看| 久久久久国产网址| 久久久久久伊人网av| 亚洲欧美清纯卡通| 亚洲,欧美,日韩| 哪个播放器可以免费观看大片| 免费观看性生交大片5| 精品亚洲乱码少妇综合久久| 亚洲精品乱码久久久久久按摩| 观看美女的网站| 老熟女久久久| 亚洲av中文av极速乱| 黄色毛片三级朝国网站 | 黑人猛操日本美女一级片| 91精品国产九色| 插逼视频在线观看| 在线观看美女被高潮喷水网站| 亚洲不卡免费看| 亚洲成人av在线免费| av福利片在线| 高清黄色对白视频在线免费看 | 午夜免费鲁丝| 国产精品伦人一区二区| 岛国毛片在线播放| 久久久亚洲精品成人影院| 99久久中文字幕三级久久日本| 亚洲综合精品二区| 新久久久久国产一级毛片| 色94色欧美一区二区| 七月丁香在线播放| 熟女人妻精品中文字幕| 熟女av电影| av播播在线观看一区| 国产高清有码在线观看视频| 免费黄网站久久成人精品| 成年美女黄网站色视频大全免费 | 偷拍熟女少妇极品色| 亚洲精品中文字幕在线视频 | 97超视频在线观看视频| 国产精品女同一区二区软件| 麻豆乱淫一区二区| 91精品一卡2卡3卡4卡| 国产一区二区三区av在线| 一级二级三级毛片免费看| 国产日韩欧美视频二区| 久热这里只有精品99| 亚洲激情五月婷婷啪啪| 视频中文字幕在线观看| 成年av动漫网址| 中文字幕人妻丝袜制服| 日日撸夜夜添| 国产免费一区二区三区四区乱码| 纯流量卡能插随身wifi吗| 尾随美女入室| 亚洲激情五月婷婷啪啪| 国产成人精品婷婷| 王馨瑶露胸无遮挡在线观看| 中文字幕制服av| 日日摸夜夜添夜夜添av毛片| 午夜av观看不卡| 婷婷色综合www| 日韩成人伦理影院| 啦啦啦在线观看免费高清www| 性色av一级| 亚洲第一区二区三区不卡| 视频中文字幕在线观看| 国产乱人偷精品视频| 国产成人精品久久久久久| 人人澡人人妻人| videos熟女内射| 波野结衣二区三区在线| 亚洲久久久国产精品| 久久毛片免费看一区二区三区| 国产在线免费精品| 亚洲av欧美aⅴ国产| 国产成人免费观看mmmm| 成人特级av手机在线观看| 性色avwww在线观看| 九九爱精品视频在线观看| 少妇高潮的动态图| 亚洲天堂av无毛| 丰满乱子伦码专区| 女的被弄到高潮叫床怎么办| 日韩av在线免费看完整版不卡| 国模一区二区三区四区视频| 亚洲欧美日韩东京热| 在线观看三级黄色| 久久久久久人妻| 三级国产精品欧美在线观看| 插阴视频在线观看视频| 成人亚洲精品一区在线观看| 人人澡人人妻人| 精品一区二区三卡| 久久精品国产自在天天线| av女优亚洲男人天堂| 国产黄片美女视频| 人人澡人人妻人| 久久狼人影院| 自拍偷自拍亚洲精品老妇| av在线播放精品| 国产老妇伦熟女老妇高清| 日本wwww免费看| 又大又黄又爽视频免费| 免费观看性生交大片5| 免费看av在线观看网站| 日日摸夜夜添夜夜添av毛片| 街头女战士在线观看网站| 天堂8中文在线网| a级毛色黄片| 一个人看视频在线观看www免费| 亚洲天堂av无毛| 国产在视频线精品| 国产日韩欧美在线精品| 久久这里有精品视频免费| a 毛片基地| 亚洲性久久影院| 寂寞人妻少妇视频99o| 一级毛片 在线播放| 毛片一级片免费看久久久久| av线在线观看网站| 欧美激情极品国产一区二区三区 | 国产精品熟女久久久久浪| 内射极品少妇av片p| 一级片'在线观看视频| 久久精品国产亚洲av天美| 91aial.com中文字幕在线观看| 下体分泌物呈黄色| 国产熟女欧美一区二区| 国产精品久久久久久av不卡| 少妇高潮的动态图| 国产av国产精品国产| 人妻 亚洲 视频| 99久久中文字幕三级久久日本| 纯流量卡能插随身wifi吗| 欧美亚洲 丝袜 人妻 在线| 国产片特级美女逼逼视频| 全区人妻精品视频| 精品久久久精品久久久| 亚洲经典国产精华液单| 黄色怎么调成土黄色| 国产午夜精品一二区理论片| 高清av免费在线| 伊人久久国产一区二区| 一区二区av电影网| 国产精品久久久久久精品古装| 寂寞人妻少妇视频99o| 午夜91福利影院| 大香蕉久久网| 校园人妻丝袜中文字幕| 国产探花极品一区二区| 久久人人爽av亚洲精品天堂| 国产又色又爽无遮挡免| 精品一区在线观看国产| 午夜免费男女啪啪视频观看| 夫妻午夜视频| av视频免费观看在线观看| 我的女老师完整版在线观看| 熟妇人妻不卡中文字幕| 精品久久久久久久久亚洲| 亚洲精品成人av观看孕妇| 纵有疾风起免费观看全集完整版| 麻豆成人午夜福利视频| 边亲边吃奶的免费视频| 国语对白做爰xxxⅹ性视频网站| 黄色配什么色好看| 黄色毛片三级朝国网站 | 性色avwww在线观看| 伦理电影大哥的女人| 嫩草影院新地址| 国产在线男女| 亚洲精品国产色婷婷电影| 在线观看免费视频网站a站| 国产精品人妻久久久影院| 99久久人妻综合| 国产 一区精品| 天堂俺去俺来也www色官网| 亚洲精品自拍成人| 日韩 亚洲 欧美在线| 午夜福利网站1000一区二区三区| 日韩欧美精品免费久久| 国产白丝娇喘喷水9色精品| 在线观看美女被高潮喷水网站| 日韩在线高清观看一区二区三区| 国产日韩一区二区三区精品不卡 | 色94色欧美一区二区| 日日爽夜夜爽网站| www.av在线官网国产| 久久人人爽av亚洲精品天堂| 内射极品少妇av片p| 久久人人爽av亚洲精品天堂| 免费看光身美女| 色视频在线一区二区三区| 国产毛片在线视频| 亚洲精品久久午夜乱码| 精品酒店卫生间| 日本av手机在线免费观看| 国产伦在线观看视频一区| 成人漫画全彩无遮挡| 蜜桃在线观看..| 国产亚洲午夜精品一区二区久久| 国产精品国产av在线观看| 国产男女内射视频| 人人妻人人看人人澡| a级毛片免费高清观看在线播放| 大又大粗又爽又黄少妇毛片口| 又大又黄又爽视频免费| 精品亚洲乱码少妇综合久久| 色吧在线观看| 夜夜爽夜夜爽视频| 免费不卡的大黄色大毛片视频在线观看| 精品亚洲成a人片在线观看| 亚洲av.av天堂| 国产深夜福利视频在线观看| 一级毛片 在线播放| 丁香六月天网| 又粗又硬又长又爽又黄的视频| 99热这里只有是精品50| 欧美精品国产亚洲| 成人午夜精彩视频在线观看| av一本久久久久| 国产成人午夜福利电影在线观看| 偷拍熟女少妇极品色| 麻豆成人av视频| 一级二级三级毛片免费看| 日韩成人伦理影院| 国产高清有码在线观看视频| 国产精品一区二区在线观看99| 久久青草综合色| 亚洲久久久国产精品| 亚洲精品久久午夜乱码| 色94色欧美一区二区| 18禁在线无遮挡免费观看视频| 国产男女超爽视频在线观看| 国模一区二区三区四区视频| 制服丝袜香蕉在线| 中文乱码字字幕精品一区二区三区| 我的老师免费观看完整版| 热99国产精品久久久久久7| 国产免费福利视频在线观看| 欧美xxxx性猛交bbbb| 在线观看av片永久免费下载| 99热全是精品| 少妇人妻一区二区三区视频| 在线亚洲精品国产二区图片欧美 | 精品国产国语对白av| 99精国产麻豆久久婷婷| 欧美 亚洲 国产 日韩一| 我的老师免费观看完整版| 精品人妻熟女av久视频| 久久99蜜桃精品久久| 热re99久久精品国产66热6| 午夜日本视频在线| 在线观看www视频免费| 日本色播在线视频| 亚洲美女搞黄在线观看| 国产成人免费观看mmmm| 免费人成在线观看视频色| 免费久久久久久久精品成人欧美视频 | 久久99一区二区三区| 日韩亚洲欧美综合| 日韩av不卡免费在线播放| 51国产日韩欧美| 婷婷色麻豆天堂久久| 少妇熟女欧美另类| 3wmmmm亚洲av在线观看| videos熟女内射| 18禁裸乳无遮挡动漫免费视频| 免费少妇av软件| 国国产精品蜜臀av免费| 在线精品无人区一区二区三| 久久国产亚洲av麻豆专区| 亚洲成色77777| 婷婷色综合www| 91久久精品电影网| 秋霞在线观看毛片| 日韩伦理黄色片| av在线app专区| 欧美精品国产亚洲| 日日摸夜夜添夜夜添av毛片| 黄色毛片三级朝国网站 | 最近的中文字幕免费完整| 这个男人来自地球电影免费观看 | 三级经典国产精品| 亚洲精品国产av成人精品| 精品国产露脸久久av麻豆| 在线 av 中文字幕| 国产黄频视频在线观看| 欧美成人午夜免费资源| 欧美精品国产亚洲| 日日摸夜夜添夜夜添av毛片| 91精品伊人久久大香线蕉| 中文字幕人妻丝袜制服| 亚洲精品第二区| 一区在线观看完整版| 老熟女久久久| 日韩av免费高清视频| a级片在线免费高清观看视频| 综合色丁香网| 人体艺术视频欧美日本| 国产黄色视频一区二区在线观看| 欧美精品一区二区免费开放| 伊人久久国产一区二区| 国产精品三级大全| 不卡视频在线观看欧美| 插阴视频在线观看视频| 男人舔奶头视频| 欧美成人精品欧美一级黄| 久久久久久伊人网av| 成年女人在线观看亚洲视频| 天美传媒精品一区二区| 九九在线视频观看精品| a 毛片基地| 深夜a级毛片| 免费黄网站久久成人精品| 亚州av有码| 免费看不卡的av| 日日爽夜夜爽网站| 成年美女黄网站色视频大全免费 | a级片在线免费高清观看视频| 日日爽夜夜爽网站| 男女免费视频国产| 免费在线观看成人毛片| 欧美日韩在线观看h| 国产精品秋霞免费鲁丝片| 中文字幕免费在线视频6| 国产成人91sexporn| 精品一品国产午夜福利视频| 少妇裸体淫交视频免费看高清| 国产精品.久久久| 亚洲精品国产av成人精品| 国产免费福利视频在线观看| 国产精品国产三级国产专区5o| 又大又黄又爽视频免费| 男女无遮挡免费网站观看| 久久综合国产亚洲精品| 熟女电影av网| 人妻系列 视频| 国产欧美日韩综合在线一区二区 | 日韩不卡一区二区三区视频在线| 久久狼人影院| 91在线精品国自产拍蜜月| 国产高清不卡午夜福利| 免费黄色在线免费观看| 搡女人真爽免费视频火全软件| 午夜影院在线不卡| 国产综合精华液| 99视频精品全部免费 在线| 国产 精品1| 日本av免费视频播放| 中文在线观看免费www的网站| 久久久国产一区二区| 伦精品一区二区三区| 日韩欧美精品免费久久| 汤姆久久久久久久影院中文字幕| 熟妇人妻不卡中文字幕| 夫妻午夜视频| 久久久久久久久久久免费av| 国产美女午夜福利| 亚洲欧美成人综合另类久久久| 免费少妇av软件| 色视频在线一区二区三区| 夜夜看夜夜爽夜夜摸| 青春草国产在线视频| 亚洲伊人久久精品综合| 国产精品国产av在线观看| 男男h啪啪无遮挡| 久久久久久久久久久久大奶| 如何舔出高潮| 2018国产大陆天天弄谢| 免费看光身美女| 一本一本综合久久| 精华霜和精华液先用哪个| 少妇精品久久久久久久| 成人漫画全彩无遮挡| 99视频精品全部免费 在线| 韩国av在线不卡| 菩萨蛮人人尽说江南好唐韦庄| 99久久精品一区二区三区| 国产成人免费无遮挡视频| 中文天堂在线官网| 国产爽快片一区二区三区| 欧美精品国产亚洲| 黄色配什么色好看| 黄色视频在线播放观看不卡| 久久精品国产自在天天线| 久久精品熟女亚洲av麻豆精品| 久久久久久久精品精品| 国产成人精品久久久久久| 少妇熟女欧美另类| 成人国产av品久久久| 全区人妻精品视频| 热re99久久精品国产66热6| 亚洲精品乱久久久久久| 一个人看视频在线观看www免费| 日本欧美国产在线视频| 婷婷色综合大香蕉| 秋霞伦理黄片| 亚洲婷婷狠狠爱综合网| 国产精品国产三级专区第一集| av黄色大香蕉| 中国国产av一级| 欧美日韩综合久久久久久| 亚洲精品第二区| 七月丁香在线播放| 欧美激情国产日韩精品一区| 在线 av 中文字幕| 日韩精品有码人妻一区| 午夜免费男女啪啪视频观看| 精品久久久精品久久久| 欧美日韩在线观看h| 精品一区二区三卡| 在线观看免费高清a一片| 久久久亚洲精品成人影院| 国产乱人偷精品视频| 欧美性感艳星| 寂寞人妻少妇视频99o| 中文字幕亚洲精品专区| 国产亚洲午夜精品一区二区久久| 老司机影院毛片| 日本黄色片子视频| 亚洲精品色激情综合| 久久久久国产网址| 欧美日韩av久久| 极品人妻少妇av视频| 国产伦在线观看视频一区| 女性被躁到高潮视频| 少妇丰满av| 国产精品人妻久久久久久| 99久久精品国产国产毛片| 日韩中字成人| 另类亚洲欧美激情| 一级毛片我不卡| 久久99蜜桃精品久久| 男人添女人高潮全过程视频| 大片免费播放器 马上看| 在线观看免费日韩欧美大片 | 91久久精品电影网| 久久精品国产鲁丝片午夜精品| 一级毛片久久久久久久久女| 日韩欧美 国产精品| 自线自在国产av| 亚洲欧洲精品一区二区精品久久久 | 国产成人一区二区在线| 国产老妇伦熟女老妇高清| 九九久久精品国产亚洲av麻豆| 尾随美女入室| 男人爽女人下面视频在线观看| 99久久精品热视频| 欧美人与善性xxx| 国产片特级美女逼逼视频| 国产一区二区在线观看av| 高清视频免费观看一区二区| 日韩亚洲欧美综合| 九九久久精品国产亚洲av麻豆| 高清视频免费观看一区二区| a 毛片基地| 女人精品久久久久毛片| 99re6热这里在线精品视频| 国产片特级美女逼逼视频| 亚洲伊人久久精品综合| 精品人妻熟女毛片av久久网站| 久久精品夜色国产| av又黄又爽大尺度在线免费看| 岛国毛片在线播放|