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

    滇東北筇竹寺組極低級(jí)變質(zhì)作用及古地溫研究

    2025-03-04 00:00:00劉建清趙瞻何利陸俊澤冉敬胡志中
    沉積學(xué)報(bào) 2025年1期

    摘 要 【目的】極低級(jí)變質(zhì)作用處于石油—濕氣—干氣以及泥煤—煙煤—無(wú)煙煤的演化階段,與古地溫有著極其密切的關(guān)系,因而,極低級(jí)變質(zhì)作用及古地溫的研究對(duì)于油氣、煤等能源礦產(chǎn)的勘探具有重要的現(xiàn)實(shí)意義。與已取得勘探突破的上揚(yáng)子地區(qū)龍馬溪組相比,下寒武統(tǒng)牛蹄塘組(筇竹寺組)頁(yè)巖有機(jī)質(zhì)豐度和熱演化程度更高,具有更大的沉積厚度和更廣泛的分布面積,頁(yè)巖氣資源潛力巨大,是中國(guó)頁(yè)巖氣勘探開(kāi)發(fā)的另一個(gè)重點(diǎn)層系。目前,筇竹寺組這套高成熟度烴源巖尚未開(kāi)展系統(tǒng)的極低級(jí)變質(zhì)作用及古地溫研究,有必要揭示古地溫等技術(shù)參數(shù)及其對(duì)油氣成藏的影響?!痉椒ā吭谄拭鎸?shí)測(cè)基礎(chǔ)上,通過(guò)伊/蒙混層礦物結(jié)構(gòu)有序度、伊利石結(jié)晶度和烴源巖鏡質(zhì)體反射率等技術(shù)被首次開(kāi)展了滇東北地區(qū)筇竹寺組極低級(jí)變質(zhì)作用及古地溫研究。【結(jié)果】(1)筇竹寺組粉砂質(zhì)泥巖中伊/蒙混層礦物中伊利石含量達(dá)95%。自生伊利石及伊/蒙混層礦物中伊利石結(jié)晶度介于0.24~0.35;由瀝青反射率計(jì)算的烴源巖鏡質(zhì)體反射率(Ro)介于2.22%~2.65%,低于揚(yáng)子?xùn)|南緣(Ro介于3.18%~3.43%)及海相干酪根天然氣生成成熟度上限(Ro=3%);(2)滇東北地區(qū)筇竹寺經(jīng)歷最高古地溫確定在200 ℃~227.57 ℃,低于揚(yáng)子?xùn)|南緣(211 ℃~246 ℃)及烴類(lèi)死亡最高溫度(230 ℃~300 ℃)。早白堊世—新近紀(jì)最大埋深期古地溫梯度介于2.53~2.91 ℃/100m,與揚(yáng)子?xùn)|南緣接近;(3)綜合考慮各變質(zhì)作用劃分標(biāo)準(zhǔn)及方案,滇東北地區(qū)筇竹寺組分別屬于極低級(jí)變質(zhì)A帶—B帶、近變質(zhì)帶、低級(jí)近變帶—高級(jí)近變帶,系葡萄石—綠纖石相?!窘Y(jié)論】滇東北地區(qū)筇竹寺期處于由沉積凹陷向隆起過(guò)渡的有利指向區(qū),該區(qū)筇竹寺組經(jīng)歷了較低的最高古地溫,且有機(jī)質(zhì)成熟度相對(duì)于揚(yáng)子其他地區(qū)成熟度更低,具有油氣形成條件和物質(zhì)基礎(chǔ)。目前尚未發(fā)現(xiàn)古油藏出露,后期構(gòu)造保存條件較好,滇東北地區(qū)筇竹寺組的油氣遠(yuǎn)景應(yīng)予以關(guān)注。

    關(guān)鍵詞 滇東北;筇竹寺組;伊利石;鏡質(zhì)體反射率;極低級(jí)變質(zhì)作用;古地溫

    第一作者簡(jiǎn)介 劉建清,男,1969年出生,教授級(jí)高級(jí)工程師,巖石學(xué),E-mail: ljianqing0813@163.com

    通信作者 趙瞻,男,教授級(jí)高級(jí)工程師,E-mail: zhaozhan1982@163.com

    中圖分類(lèi)號(hào) P618.13 文獻(xiàn)標(biāo)志碼 A

    0 引言

    極低級(jí)變質(zhì)作用是變質(zhì)地質(zhì)學(xué)研究的前沿領(lǐng)域之一,指很低溫的一種變質(zhì)作用,即從成巖作用向低級(jí)變質(zhì)的轉(zhuǎn)變過(guò)程,與石油、天然氣、煤等能源礦產(chǎn)密切相關(guān)。在成巖—極低級(jí)變質(zhì)作用階段,石油演化階段從石油—濕氣—干氣變化,煤則經(jīng)歷由泥煤—煙煤—無(wú)煙煤的變化。這些過(guò)程有著較多評(píng)價(jià)指標(biāo),例如可依據(jù)伊利石結(jié)晶度KI(°Δ2θ)0.42和0.25將極低級(jí)變質(zhì)作用劃分為成巖帶、近變質(zhì)帶和淺變質(zhì)帶。因此,通過(guò)對(duì)沉積盆地成巖—極低級(jí)變質(zhì)作用研究,對(duì)于明確成巖—變質(zhì)界限,闡明盆地—構(gòu)造—熱演化史,指導(dǎo)油氣勘探,具有重要的理論和現(xiàn)實(shí)意義[1?5]。

    國(guó)際上極低級(jí)變質(zhì)作用的主要標(biāo)志有:(1)標(biāo)志礦物及組合;(2)伊利石和伊/蒙轉(zhuǎn)換,如伊利石結(jié)晶度、多型,伊/蒙混層可膨脹層的百分比,伊利石有序度、晶疇大小、b0值等;(3)有機(jī)質(zhì)參數(shù),如鏡質(zhì)體反射率、筆石反射率、牙形刺顏色蝕變指數(shù)、孢粉顏色等[6?15]。在油氣的生成和轉(zhuǎn)變過(guò)程中,含油氣盆地古地溫及古地溫梯度是熱史、生烴史研究的重要基礎(chǔ),不僅可以重塑油藏歷史演化過(guò)程,而且可以預(yù)測(cè)其發(fā)展演化趨勢(shì)。近年來(lái)的研究,可依據(jù)上述指標(biāo)開(kāi)展烴源巖埋藏過(guò)程古地溫定量計(jì)算,精準(zhǔn)標(biāo)定其經(jīng)歷的最高古地溫,為熱史、生烴史研究奠定基礎(chǔ)。

    中上揚(yáng)子地區(qū)下寒武統(tǒng)筇竹寺組(牛蹄塘組)普遍發(fā)育一套黑色烴源巖,厚度大、有機(jī)碳含量高,是四川盆地眾多下古生界油氣藏的烴源巖,受到石油地質(zhì)界的廣泛關(guān)注。以往對(duì)于該套烴源巖成熟度的工作,主要根據(jù)鏡質(zhì)體反射率(Ro)進(jìn)行成熟、過(guò)成熟的定性評(píng)價(jià),未曾開(kāi)展較為系統(tǒng)的低級(jí)變質(zhì)作用及其相關(guān)熱歷史的研究。滇東北昭通地區(qū)位于上揚(yáng)子地區(qū),筇竹寺組烴源巖鏡質(zhì)體反射率低于揚(yáng)子?xùn)|南緣[16],成熟度較低,更具油氣勘探的現(xiàn)實(shí)意義。該區(qū)沉積盆地的研究缺乏古地溫資料,制約了盆地生烴史的研究。在滇東北昭通地區(qū)鋅廠溝筇竹寺組剖面測(cè)量基礎(chǔ)上,開(kāi)展了筇竹寺組黑色頁(yè)巖黏土礦物組成、伊利石結(jié)晶度、鏡質(zhì)體反射率研究及與變質(zhì)程度和古地溫定量計(jì)算,旨在揭示該區(qū)筇竹寺組烴源巖經(jīng)歷的最大古地溫,為盆地生烴史的研究提供基礎(chǔ)資料,深化該區(qū)筇竹寺組的研究,服務(wù)于油氣勘探。

    1 地質(zhì)背景

    1.1 地質(zhì)構(gòu)造背景

    根據(jù)云南省地質(zhì)志的劃分,研究區(qū)在一級(jí)構(gòu)造單位為揚(yáng)子準(zhǔn)地臺(tái),二級(jí)構(gòu)造單元為滇東臺(tái)褶帶[17]。經(jīng)過(guò)晉寧運(yùn)動(dòng)之后的加里東、海西、印支、燕山、喜馬拉雅各構(gòu)造期,該區(qū)沉積了震旦系、古生界、中生界、新生界的巨厚蓋層(圖1)。區(qū)內(nèi)出露的震旦系—白堊系均為整合或平行不整合接觸,下白堊統(tǒng)與新近系之間為角度不整合接觸,表明該區(qū)沉積蓋層于晚燕山—早喜山期褶皺定型。其中上古生界與下古生界的廣西運(yùn)動(dòng)造成的平行不整合是區(qū)內(nèi)最大的沉積間斷[18]。

    早寒武世筇竹寺期是上揚(yáng)子地區(qū)早古生代最大的一次海侵,普遍發(fā)育了厚度較大的黑色頁(yè)巖。滇東北地區(qū)位于上揚(yáng)子康滇古陸東緣,筇竹寺期接受來(lái)自該古陸的剝蝕和沉積。同時(shí),受新元古代裂谷盆地的持續(xù)影響,筇竹寺組沉積基底凹凸不平,自該區(qū)至鎮(zhèn)雄一帶筇竹寺組沉積厚度300~500 m,各地黑色頁(yè)巖累計(jì)厚度30~60 m,沉降中心位于鎮(zhèn)雄一帶。自晚寒武世開(kāi)始,在華南加里東幕式造山作用背景下,上揚(yáng)子地區(qū)呈現(xiàn)隆、凹相間的格局,發(fā)育一系列的古隆起。滇東北昭通地區(qū)位于其中的黔中隆起,區(qū)內(nèi)奧陶系—志留系缺失或發(fā)育不完整,而往北漸入黔北凹陷[19]。

    1.2 剖面特征

    為開(kāi)展下寒武統(tǒng)筇竹寺組的極低級(jí)變質(zhì)作用及相關(guān)最高古地溫研究,實(shí)測(cè)了云南省昭通市昭陽(yáng)區(qū)鋅廠溝筇竹寺組剖面。剖面位于金沙江東岸一支溝(鋅廠溝)(剖面起點(diǎn)坐標(biāo):103°13′46.21″ E,27°34′13.29″ N;終點(diǎn)坐標(biāo):103°14′04.03″ E,27°34′10.60″ N),見(jiàn)頂、底,底與下寒武統(tǒng)麥地坪組厚層粉晶白云巖為平行不整合接觸,頂與中寒武統(tǒng)滄浪鋪組紫紅色泥巖呈整合接觸。剖面共分29層,0層麥地坪組,29層為滄浪鋪組,除7~8層局部風(fēng)化和掩蓋較強(qiáng)外,其余露頭連續(xù)出露,剖面質(zhì)量較好。剖面縱向上兩分明顯(圖2),底部(1~5層)為灰黑色—黑色中厚層狀—薄層狀炭質(zhì)泥巖、炭質(zhì)粉砂巖、含炭質(zhì)泥巖,為烴源巖層(圖3),沉積相分析為深水陸棚環(huán)境,厚19.57 m;上部(6~28層)為灰—深灰色泥質(zhì)粉砂巖、深灰色泥巖、鈣泥質(zhì)粉砂巖夾粉砂質(zhì)泥巖,風(fēng)化面略帶淺黃灰色、淺黃綠色,普遍發(fā)育水平層理,局部(27層)見(jiàn)浪成對(duì)稱(chēng)波痕(圖3)。其中17層、27層深灰色泥巖為較差烴源巖層,厚度分別為6.46 m和13.17 m。沉積相主要為淺水陸棚,局部(17層、27層)為深水陸棚環(huán)境,反映上部有兩次短暫的水體加深過(guò)程[19]。

    2 樣品采集與分析方法

    為開(kāi)展滇東北地區(qū)筇竹寺組的極低級(jí)變質(zhì)作用及相關(guān)熱歷史研究,在鋅廠溝筇竹寺組剖面采集泥頁(yè)巖進(jìn)行黏土礦物X射線衍射分析、伊利石結(jié)晶度分析,對(duì)筇竹寺組底部黑色泥質(zhì)粉砂巖及中上部深灰色頁(yè)巖采集樣品進(jìn)行瀝青反射率分析。樣品在四川省科源工程技術(shù)測(cè)試中心(頁(yè)巖氣評(píng)價(jià)與開(kāi)采四川省重點(diǎn)實(shí)驗(yàn)室)完成。

    黏土礦物的提取、制備及分析流程按《SYT5163-208:沉積巖中黏土礦物和常見(jiàn)非黏土礦物X射線衍射分析方法》進(jìn)行開(kāi)展。黏土礦物X射線衍射分析在DMX-Ⅲ型X射線衍射儀上完成,測(cè)試條件為Cu靶,電壓35 kV,電流15 mA,步寬0.02° 2θ,掃描范圍2°~36° 2θ。黏土礦物含量利用衍射峰強(qiáng)度計(jì)算方法獲得[20]。

    伊利石結(jié)晶度測(cè)試儀器為X’pert Powder,測(cè)試條件為Cu靶,電壓40 kV,電流40 mA;發(fā)射狹縫與散射狹縫均為1°,接收狹縫0.3 mm,掃描速度:2θ為2°/min;采樣步寬:2θ為0.02°;掃描范圍5°~45°(自然定向片,飽和乙二醇片)。

    干酪根分選、制備流程依據(jù)國(guó)家標(biāo)準(zhǔn)GB/T19144—2010《沉積巖中干酪根分離方法》開(kāi)展。瀝青反射率儀器為Axio Scope.A1顯微鏡、Jamp;M分光光度計(jì),放大倍數(shù)500 倍,反射率量程0.1%~10%。測(cè)試條件:環(huán)境溫度23±3 ℃,相對(duì)濕度小于70%RH,光線波長(zhǎng)546 nm±5 nm。標(biāo)準(zhǔn)物質(zhì):釓鎵石榴石,浸油反射率1.72%,氧化鋯,浸油反射率3.16%(在23 ℃±1 ℃,546 nm波長(zhǎng)的綠波下,折射率為1.518 0±0.000 4)。測(cè)試執(zhí)行標(biāo)準(zhǔn)為SY/5124—2012《沉積巖中鏡質(zhì)體反射率測(cè)定方法》[5]。

    3 分析結(jié)果

    3.1 巖石學(xué)特征

    鋅廠溝剖面筇竹寺組泥質(zhì)粉砂巖、粉砂質(zhì)泥巖主要由石英、鉀長(zhǎng)石、斜長(zhǎng)石、方解石、白云石、黏土礦物組成(表1)。全巖X射線衍射顯示,石英含量33%~45%,平均40%,鉀長(zhǎng)石含量0~14%,平均6.6%,斜長(zhǎng)石含量0~30%,平均14.4%,方解石含量0~18%,平均6.3%,白云石含量0~15%,平均6.60%,黏土礦物含量12%~39%,平均26.30%。顯微鏡下,石英呈棱角狀或次棱角狀,粒徑0.05 mm左右,顆粒間充填泥質(zhì)、碳質(zhì)等成分;長(zhǎng)石顆粒,棱角狀,粒徑0.05 mm左右,局部可見(jiàn)聚片雙晶;方解石和白云石呈泥晶—微晶集合體,為沉積期沉淀的碳酸鹽礦物;“泥質(zhì)”呈基底式膠結(jié)石英碎屑,呈隱晶質(zhì)狀或重結(jié)晶形成雛晶狀,為沉積期沉淀并經(jīng)成巖—極低級(jí)轉(zhuǎn)變黏土礦物和自生黏土礦物,其中夾雜少量碳質(zhì)成分[19](圖4)。鏡下筇竹寺組泥質(zhì)粉砂巖、粉砂質(zhì)泥巖見(jiàn)清晰的層理構(gòu)造(S0),無(wú)置換性變質(zhì)面理(S1),長(zhǎng)石、方解石顆粒亦無(wú)明顯的壓扁、拉長(zhǎng)等變質(zhì)作用的特征,反映構(gòu)造擠壓作用并不強(qiáng)烈,巖石尚處于成巖—極低級(jí)變質(zhì)作用,未達(dá)到淺變質(zhì)作用階段。

    鋅廠溝剖面筇竹寺組泥頁(yè)巖黏土礦物分析結(jié)果如表2 所示。黏土礦物主要為伊利石,含量58%~94%,次為綠泥石(含量1%~35%)及伊/蒙混層(2%~6%),另含少量蒙脫石(0~7%)及高嶺石(0~16%)。伊/蒙混層中伊利石占95%,蒙脫石占5%。伊利石結(jié)晶度KI(°Δ2θ)0.24~0.35。

    筇竹寺組物源來(lái)自康滇古陸古—中元古界東川群、會(huì)理群淺變質(zhì)綠片巖相火山—沉積巖系[18],最低級(jí)的變質(zhì)泥質(zhì)礦物為絹云母及白云母,而鋅廠溝剖面黏土礦物組成為經(jīng)歷成巖作用及極低級(jí)變質(zhì)作用的伊利石、伊/蒙混層、高嶺石及綠泥石等形成溫度、壓力較低的礦物,最高級(jí)的變質(zhì)礦物為伊利石,與源區(qū)黏土礦物明顯不同。同時(shí),從物源區(qū)綠片巖相的變質(zhì)程度可知,其所夾碳質(zhì)Ro可超過(guò)4%。在化學(xué)風(fēng)化過(guò)程中鉀長(zhǎng)石、絹云母、白云母可徹底轉(zhuǎn)化為高嶺石,蒙脫石在風(fēng)化帶相當(dāng)穩(wěn)定,可帶至沉積盆地,黑云母可轉(zhuǎn)化為蒙脫石。而基性巖風(fēng)化過(guò)程中鐵、鎂、鈣質(zhì)大多呈溶液狀態(tài)被帶走[21]。滇東北地區(qū)筇竹寺烴源巖瀝青反射率換算的Ro 明顯低于4%,表明源區(qū)泥質(zhì)礦物經(jīng)歷了較為徹底的化學(xué)風(fēng)化作用生成如高嶺石、蒙脫石等新的礦物,搬運(yùn)至康滇古陸東緣沉積、埋藏,經(jīng)歷成巖和極低級(jí)變質(zhì)。換言之,源區(qū)黑云母、白云母、絹云母等成熟度較高古老黏土礦物在到達(dá)沉積盆地時(shí)已轉(zhuǎn)化為高嶺石、蒙脫石等成熟度較低的黏土礦物,源區(qū)高成熟度黏土礦物對(duì)沉積盆地黏土礦物成熟度的影響可以忽略。

    3.2 有機(jī)質(zhì)成熟度特征

    鏡質(zhì)體反射率是確定烴源巖有機(jī)質(zhì)成熟度的有效指標(biāo)。由于鏡質(zhì)組分源于高等植物,該方法明顯受到生物演化及母質(zhì)來(lái)源類(lèi)型的限制。碳酸鹽烴源巖或下古生界烴源巖中缺乏陸相高等植物來(lái)源的鏡質(zhì)體,一般采用海相鏡質(zhì)體或原生瀝青來(lái)代替[22]。鋅廠溝筇竹寺組剖面瀝青反射率Rran 介于2.81%~3.45%(表3)。劉德漢等[23]針對(duì)我國(guó)海相高演化碳酸鹽巖烴源巖提出了瀝青反射率(Rran)與等效鏡質(zhì)體反射率(Ro)的換算公式:

    Ro=0.668×Rran+0.346 (1)

    式中:Ro為鏡質(zhì)體反射率,Rran為瀝青體反射率。這一方法也廣泛應(yīng)用于寒武系—志留系等缺乏陸相高等植物來(lái)源地層烴源巖鏡質(zhì)體反射率的轉(zhuǎn)換計(jì)算,并取得較好的應(yīng)用效果[16]。計(jì)算表明,鋅廠溝剖面筇竹寺組烴源巖鏡質(zhì)體反射率介于2.22%~2.65%(表3)。

    4 溫度計(jì)算結(jié)果

    4.1 伊/蒙混層溫度計(jì)算結(jié)果

    伊(I)/蒙(S)間層礦物中伊利石含量與溫度的相關(guān)性,相關(guān)學(xué)者研究結(jié)果如表4所示。鋅廠溝剖面筇竹寺組泥頁(yè)巖伊/蒙混層礦物伊利石層占95%,綜合表4各學(xué)者的研究結(jié)果,該區(qū)伊/蒙混層礦物形成溫度應(yīng)在200 ℃以上。

    王強(qiáng)[28]曾應(yīng)用Chi et al.[29]根據(jù)Harvey et al.[27]的原始數(shù)據(jù)建立了伊/蒙混層中伊利石含量與形成溫度的線性表達(dá)式開(kāi)展黏土礦物研究,取得較好成效,該表達(dá)式如下:

    Y=0.342+2.737×10-3T (2)

    式中:Y 為伊/蒙混層中伊利石的百分含量(%),T 表示伊/蒙混層礦物的形成溫度(℃)[30]。計(jì)算結(jié)果如表2中T1所示,為222.14 ℃。

    4.2 伊利石結(jié)晶度溫度計(jì)算結(jié)果

    Weaver[31]首先觀察到,伊利石的X射線衍射峰的銳度隨著成巖及變質(zhì)程度的增高而增高,在此基礎(chǔ)上提出了伊利石結(jié)晶度的概念。當(dāng)前國(guó)際上通用Kübler[10?11]提出的測(cè)量伊利石第一級(jí)底面反射峰001峰的半高寬(又稱(chēng)Kübler指數(shù))來(lái)反映伊利石的結(jié)晶度,用KI(°Δ2θ)表示,KI 值越小,峰形越尖銳,伊利石的結(jié)晶度越好。

    胡大千等[5]、畢先梅等[32]歸納總結(jié)了伊利石結(jié)晶度與形成溫度之間的關(guān)系,提出了成巖—極低級(jí)變質(zhì)作用的劃分方案:晚成巖帶(高級(jí)成巖帶),伊利石結(jié)晶度介于1.00~0.42,溫度上限200 ℃;低級(jí)近變質(zhì)帶(極低級(jí)變質(zhì)帶A),伊利石結(jié)晶度介于0.42~0.30;高級(jí)近變質(zhì)帶(極低級(jí)變質(zhì)帶B),伊利石結(jié)晶度變化介于0.30~0.25,溫度上限300 ℃~350 ℃。鋅廠溝剖面筇竹寺組泥質(zhì)粉砂巖伊利石結(jié)晶度介于0.24~0.35(表2),據(jù)此分析,該區(qū)泥質(zhì)粉砂巖處于極低級(jí)變質(zhì)帶A帶—B帶,伊利石形成溫度介于200 ℃~300 ℃。

    Chi et al.[33]研究了Wairakei地?zé)嵯到y(tǒng)沉積物中各種組分與結(jié)構(gòu)的黏土礦物,包括碎屑伊利石、綠泥石、伊/蒙混層礦物等,他們對(duì)該區(qū)自生伊利石的結(jié)晶度值和測(cè)井溫度進(jìn)行了系統(tǒng)研究,表明兩者之間有明顯的線性關(guān)系,即隨著母巖溫度或深度的增加,自生的伊利石結(jié)晶度值KI(°Δ2θ)降低。在實(shí)測(cè)數(shù)據(jù)基礎(chǔ)上,提出了以下線性表達(dá)式[30]:

    T=249-89.3×KI (R2=0.92) (3)

    式中:T 表示伊利石形成溫度(℃),KI 表示伊利石結(jié)晶度(°Δ2θ),R2表示相關(guān)系數(shù)。這個(gè)公式適用范圍為90 ℃~230 ℃。對(duì)鋅廠溝剖面筇竹寺組泥頁(yè)巖伊利石結(jié)晶度利用式(3)計(jì)算表明,伊利石形成溫度T2介于217.75 ℃~227.57 ℃(表2)。

    4.3 鏡質(zhì)體反射率溫度計(jì)算結(jié)果

    Frey[6]研究認(rèn)為:成巖帶,Ro小于2.0%,溫度上限200℃;近變質(zhì)帶,Ro為2.0%~4.0%,溫度上限300 ℃~350 ℃;淺變質(zhì)帶,Ro 大于4.0%,溫度大于300 ℃~350 ℃。早寒武世,由于缺乏高等植物,鏡質(zhì)體反射率由瀝青反射率換算[21]。鋅廠溝剖面筇竹寺組黑色頁(yè)巖瀝青反射率換算成鏡質(zhì)體反射率介于2.22%~2.65%(表3)。據(jù)此,認(rèn)為該區(qū)筇竹寺組烴源巖處于近變質(zhì)帶,粗略估計(jì)經(jīng)歷的溫度上限為200 ℃ ~300 ℃。鏡質(zhì)體反射率3%被認(rèn)為是海相干酪根天然氣生成成熟度上限[34]。在揚(yáng)子?xùn)|南緣,筇竹寺組(牛蹄塘組)烴源巖鏡質(zhì)體反射率介于3.18%~3.43%[16],高于3%的極限值,而在鋅廠溝剖面鏡質(zhì)體反射率低于3%。因此,開(kāi)展滇東北地區(qū)筇竹寺組頁(yè)巖氣勘探具有現(xiàn)實(shí)意義。顯微鏡下顯示,鋅廠溝剖面干酪根組分呈黑色(圖5),但并未完全呈深黑色,表明可能其仍具一定的生烴潛力。一般而言,干酪根顯微組分由腐泥組、殼質(zhì)組、鏡質(zhì)組、惰質(zhì)組組成。腐泥組原始有機(jī)質(zhì)為藻類(lèi)及其他低等水生生物及細(xì)菌,為腐泥化產(chǎn)物;殼質(zhì)組原始有機(jī)質(zhì)為陸生植物孢粉、花粉、角質(zhì)層、樹(shù)脂、蠟和木栓層等;鏡質(zhì)組原始有機(jī)質(zhì)來(lái)自高等植物的結(jié)構(gòu)和無(wú)結(jié)構(gòu)木質(zhì)纖維;惰質(zhì)組原始有機(jī)質(zhì)來(lái)源于森林火災(zāi)后絲炭化的木質(zhì)纖維或再沉積有機(jī)質(zhì)。寒武紀(jì)時(shí)期無(wú)高等植物,干酪根顯微組分無(wú)鏡質(zhì)組,由于時(shí)代老,也不是惰質(zhì)組和殼質(zhì)組。該區(qū)干酪根顯微組分形態(tài)為不均勻棉絮狀結(jié)構(gòu),可能為較典型的腐泥無(wú)定形體。

    Barker et al. [13]統(tǒng)計(jì)研究了鏡質(zhì)體反射率與古地溫之間的相互關(guān)系,提出了鏡質(zhì)體反射率地質(zhì)溫度計(jì):

    In(Ro%)=0.007 8Tmax-1.2 (4)

    Barker et al.[35]確立了鏡質(zhì)體反射率與古地溫的關(guān)系,給出的地質(zhì)溫度計(jì):

    T=((In(Ro%)+1.26)/0.008 11 (5)

    Barker et al.[14]修正后的鏡質(zhì)體反射率地質(zhì)溫度計(jì):

    T=(In(Ro%)+1.68)/0.012 4 (6)

    Mullis et al.[36]建立的鏡質(zhì)體反射率地質(zhì)溫度計(jì):

    T=(Ro%+8.623 8)/0.054 3[5] (7)

    上述式(4)至式(7)古溫度計(jì)算結(jié)果,分別以T3 至T6 列于表3??梢钥闯?,T3 介于256.27 ℃~278.82 ℃ ,T4 介于254.18 ℃ ~275.90 ℃ ,T5 介于199.91 ℃~214.10 ℃,T6介于199.76 ℃~207.63 ℃。

    5 討論

    5.1 最高古地溫及地溫梯度

    鋅廠溝剖面筇竹寺組伊/蒙混層礦物中伊利石含量為95%,根據(jù)Nadeau et al. [24]、Hofman et al. [25]、Jenning et al.[26]、Harvey et al.[27]關(guān)于伊/蒙混層礦物中伊利石含量與溫度關(guān)系的相關(guān)性描述(表4),認(rèn)為該區(qū)泥質(zhì)粉砂巖中伊/蒙混層礦物形成溫度在200 ℃以上。

    濰北凹陷研究表明應(yīng)用黏土礦物相關(guān)參數(shù)可較有效計(jì)算恢復(fù)古地溫[28]。式(2)、式(3)是這些方法在滇東北地區(qū)的又一實(shí)踐。式(2)根據(jù)伊/蒙混層礦物中伊利石含量計(jì)算其形成溫度T1為222.14 ℃。式(3)根據(jù)伊利石結(jié)晶度計(jì)算伊利石形成溫度T2介于217.75 ℃~227.57 ℃。

    應(yīng)用鏡質(zhì)體反射率開(kāi)展烴源巖最大古地溫計(jì)算及溫度的取舍,已在內(nèi)蒙古錫林郭勒有較成功的實(shí)踐[5]。T3至T6為依據(jù)這些方法和烴源巖鏡質(zhì)體反射率計(jì)算的古地溫,其中T3和T4溫度較高,且較為接近,溫度介于254.18 ℃~278.82 ℃。T5和T6溫度較低,介于199.76℃~214.10℃,與式(2)、式(3)及依據(jù)伊/蒙混層礦物結(jié)構(gòu)有序度估算的溫度(200℃以上)接近。注意到式(4)、式(6)均是Barker et al.[13?14]提出,式(6)是式(4)修正后的公式,且式(6)與式(2)、式(3)、式(7)及伊/蒙混層礦物結(jié)構(gòu)有序度估算的溫度接近和吻合,因此,式(6)及(2)、式(3)、式(7)計(jì)算結(jié)果及伊/蒙混層礦物結(jié)構(gòu)有序度估算的溫度可能更可信,式(4)、式(5)計(jì)算溫度可能偏高。

    將式(2)、式(3)、式(6)、式(7)計(jì)算結(jié)果及伊/蒙混層礦物結(jié)構(gòu)有序度估算的溫度綜合考慮,分析鋅廠溝剖面筇竹寺組泥頁(yè)巖經(jīng)歷的最大古地溫介于200 ℃~227.57 ℃,并未經(jīng)歷太高的古地溫。該區(qū)寒武系—下白堊統(tǒng)均為整合或平行不整合接觸,奧陶系—志留系受黔中隆起影響部分缺失,地層序列為連續(xù)的沉積或水平抬升;下白堊統(tǒng)與新近系之間為角度不整合接觸,蓋層褶皺之后,地表再次下降接受新近系沉積;之后,地表抬升至現(xiàn)今形態(tài),因此,早白堊世—新近紀(jì)為最大埋深期。根據(jù)該區(qū)奧陶—志留系部分缺失及其余地層厚度①,按地表平均溫度20 ℃[16]估算,該區(qū)筇竹寺組早白堊世—新近紀(jì)最大埋深期古地溫梯度介于2.53~2.91 ℃/100 m。

    5.2 極低級(jí)變質(zhì)作用

    對(duì)極低級(jí)變質(zhì)作用的理解及它與成巖作用和低級(jí)變質(zhì)作用間界線的認(rèn)識(shí),目前沒(méi)有取得共識(shí),而且,石油地質(zhì)學(xué)者、煤巖地質(zhì)學(xué)者、沉積和巖石學(xué)者的理解,也相差甚遠(yuǎn),所用術(shù)語(yǔ)也極不統(tǒng)一,如晚成巖作用(帶)、埋藏變質(zhì)作用(帶),及上述低級(jí)近變質(zhì)作用(帶)、高級(jí)變質(zhì)作用(帶)、近變質(zhì)作用(帶)、淺變質(zhì)作用(帶)等。但成巖階段—極低級(jí)變質(zhì)—淺變質(zhì)作用是個(gè)連續(xù)的過(guò)程,在巖石類(lèi)型、黏土礦物組成、巖石結(jié)構(gòu)構(gòu)造特征、鏡質(zhì)體反射率等方面呈連續(xù)的變化,因此,可根據(jù)連續(xù)變化特點(diǎn),選擇代表性學(xué)者研究成果開(kāi)展極低級(jí)變質(zhì)作用及油氣地質(zhì)意義研究。

    一般而言,隨著變質(zhì)級(jí)的增加,在變質(zhì)泥巖中出現(xiàn)二八面體礦物系列,即蒙脫石—伊/蒙混層—伊利石—白云母,而在變質(zhì)鎂鐵質(zhì)巖石中出現(xiàn)三八面體系列礦物,即蒙脫石—綠泥石/蒙脫石混層—綠泥石;高嶺石呈高嶺石—地開(kāi)石—葉臘石的變化;其他標(biāo)志也同步變化,微構(gòu)造由S0層理發(fā)展為S1面理;巖石類(lèi)型由頁(yè)巖、泥巖變化為筆狀構(gòu)造泥巖、瓦板巖、千枚巖;鏡質(zhì)體反射率由0.5%經(jīng)2%增至4%。

    Frey et al.[9]根據(jù)伊利石結(jié)晶度、伊/蒙混層中伊利石的含量、蒙脫石—綠泥石系列的變化、高嶺石—葉臘石系列的變化、巖石類(lèi)型、微構(gòu)造特征將成巖階段至淺變質(zhì)階段劃分為早期成巖帶、晚期成巖帶、低級(jí)近變帶、高級(jí)近變帶、淺變質(zhì)帶。滇東北地區(qū)筇竹寺組伊利石結(jié)晶度介于0.24~0.35,伊/蒙混層伊利石含量95%,黏土礦物組合高嶺石0~16%、綠泥石1%~35%、伊利石58%~94%,仍有少量早成巖帶—低級(jí)近變帶的黏土礦物存在,換算的鏡質(zhì)體反射率2.22%~2.65%,巖石為泥質(zhì)粉砂巖,未出現(xiàn)置換性面理S1,根據(jù)Frey et al.[9]的研究,將其劃分在低級(jí)近變帶—高級(jí)近變帶。索書(shū)田等[37]將成巖階段—淺變質(zhì)階段劃分為沸石相(lt;200 ℃)、葡萄石—綠纖石相(200 ℃~370 ℃)、綠片巖相(gt;370 ℃),本次計(jì)算筇竹寺組經(jīng)歷的最大古地溫為200 ℃~227.57 ℃,屬葡萄石—綠纖石相。而按照畢先梅等[32]、Frey[6]伊利石結(jié)晶度及鏡質(zhì)體反射率的劃分,則屬極低級(jí)變質(zhì)A帶—B帶、近變質(zhì)帶,巖石未達(dá)到淺變質(zhì)階段。

    5.3 油氣地質(zhì)意義

    滇東北昭通地區(qū)筇竹寺組烴源巖鏡質(zhì)體反射率(Ro)介于2.22%~2.65%,處于近變質(zhì)帶及干氣階段,明顯低于揚(yáng)子?xùn)|南緣(Ro介于3.18%~3.43%)[16]及海相干酪根天然氣生成成熟度上限(Ro=3%);顯示其經(jīng)歷的最大古地溫(200℃~227.57℃)也低于揚(yáng)子?xùn)|南緣筇竹寺組經(jīng)歷的最高古地溫(211 ℃~246 ℃)[16]及烴類(lèi)生烴死亡溫度(230 ℃~300 ℃)[31],表明該區(qū)烴源巖可能尚具微弱的生烴能力,這和揚(yáng)子其他地區(qū)截然不同的。筇竹寺組(牛蹄塘組)烴源巖的研究,目前更多地集中在揚(yáng)子?xùn)|南緣,對(duì)于滇東北這套較低成熟度的烴源巖,目前尚未予以應(yīng)有的關(guān)注。該區(qū)晚侏羅世—古近紀(jì)最大埋深期古地溫梯度介于2.53~2.91 ℃/100m,與揚(yáng)子?xùn)|南緣這一時(shí)期古地溫梯度(2.61 ℃/100m)接近[16],滇東北地區(qū)現(xiàn)今地溫梯度(2.10 ℃/100m)[38],其古地溫梯度的變化顯示滇東北地區(qū)晚侏羅世—古近紀(jì)古地溫梯度高于現(xiàn)今地溫梯度。該區(qū)烴源巖成熟度及經(jīng)歷的古地溫較低,或與該區(qū)處于黔中隆起,奧陶紀(jì)—志留紀(jì)部分地層缺失、蓋層總厚度減少以及筇竹寺組最大埋深減少、經(jīng)歷的最大古地溫和烴源巖熱演化降低有關(guān)。

    該區(qū)新近系茨營(yíng)組與下白堊統(tǒng)之間為角度不整合接觸,震旦系—下白堊統(tǒng)蓋層褶皺及斷裂系統(tǒng)在晚燕山期—早喜山期形成,是區(qū)內(nèi)最強(qiáng)烈的一次構(gòu)造運(yùn)動(dòng)。褶皺作用導(dǎo)致圈閉的形成,斷裂的疏導(dǎo)和破壞可導(dǎo)致油氣的運(yùn)移、調(diào)整、散失和破壞,并可能造成油氣藏的破壞,在地表形成一些古油藏。值得關(guān)注的是,滇東北地區(qū)褶皺及斷裂系統(tǒng)目前未有筇竹寺組古油藏的發(fā)現(xiàn)和報(bào)道,明顯不同于揚(yáng)子?xùn)|南緣斷裂已破壞油氣藏,有麻江、金沙巖孔眾多下古生界古油藏分布的特點(diǎn)。這是否意味著該區(qū)構(gòu)造作用較弱,筇竹寺組有較好的構(gòu)造保存條件?滇東北地區(qū)筇竹寺組是康滇古陸東緣寒武紀(jì)初期海侵沉積,處于由凹陷向隆起的過(guò)渡地帶,是油氣生成和運(yùn)移的有利指向區(qū),具有油氣生成基礎(chǔ)和物質(zhì)條件。同時(shí),按不同學(xué)者劃分方案分別屬于極低級(jí)變質(zhì)A帶—B帶、近變質(zhì)帶、低級(jí)近變帶—高級(jí)近變帶、葡萄石—綠纖石相。滇東北筇竹寺組烴源巖具較低成熟度,經(jīng)歷的最大古地溫較低,未達(dá)到淺變質(zhì)階段,對(duì)于油氣藏和圈閉可能尚具較弱的補(bǔ)給能力,但其較好的構(gòu)造保存條件及有利的油氣運(yùn)移指向區(qū),使得這一地區(qū)筇竹寺組的勘探前景值得關(guān)注。

    此次極低級(jí)變質(zhì)作用及古地溫和地溫梯度的研究,為該區(qū)筇竹寺組油氣成藏史的研究奠定了一定的基礎(chǔ)。未來(lái)可開(kāi)展滇東北其他烴源巖地層古地溫研究,獲取其經(jīng)歷的最大古地溫及各時(shí)代古地溫梯度,并結(jié)合本文的研究和各時(shí)代地層厚度、古水深、剝蝕量,建立該區(qū)筇竹寺組沉降史、熱史及生烴史。綜合構(gòu)造演化歷史,分析油氣生成、調(diào)整、運(yùn)移、破壞等成藏史,預(yù)測(cè)油氣發(fā)生、發(fā)展規(guī)律。

    6 結(jié)論

    (1) 滇東北地區(qū)筇竹寺組泥頁(yè)巖黏土礦物主要為伊利石,少量綠泥石、伊/蒙混層礦物、蒙脫石及高嶺石。伊/蒙混層礦物中伊利石含量達(dá)95%。自生伊利石及伊/蒙混層礦物中伊利石結(jié)晶度介于0.24~0.35;由瀝青反射率計(jì)算的烴源巖鏡質(zhì)體反射率(Ro)介于2.22%~2.65%,低于揚(yáng)子?xùn)|南緣(Ro介于3.18%~3.43%)及海相干酪根天然氣生成成熟度上限(Ro=3%)。該區(qū)部分烴源巖尚具微弱的生烴能力。

    (2) 滇東北地區(qū)筇竹寺最高古地溫確定在200 ℃~227.57 ℃,低于揚(yáng)子?xùn)|南緣(211 ℃~246 ℃)及烴類(lèi)死亡最高溫度(230 ℃~300 ℃),部分烴源巖尚具微弱生烴能力。該區(qū)晚早白堊世至新近紀(jì)最大埋深期古地溫梯度介于2.53~2.91 ℃/100 m。筇竹寺組經(jīng)歷的最高古地溫低于揚(yáng)子?xùn)|南緣,古地溫梯度與揚(yáng)子?xùn)|南緣接近。

    (3) 滇東北地區(qū)筇竹寺組可分別劃分為極低級(jí)變質(zhì)A帶—B帶、近變質(zhì)帶、低級(jí)近變帶—高級(jí)近變帶、葡萄石—綠纖石相,巖石均未達(dá)到淺變質(zhì)階段。

    (4) 研究區(qū)筇竹寺組經(jīng)歷的最高古地溫較低與位于黔中隆起北緣奧陶系—志留系部分沉積缺失、蓋層總厚度及經(jīng)歷的最大埋深降低有關(guān)。筇竹寺組沉積處于由凹陷向隆起過(guò)渡的有利指向區(qū),經(jīng)歷了較低的最高古地溫和極低級(jí)變質(zhì)作用,具較低的成熟度,具有油氣形成的條件和物質(zhì)基礎(chǔ),無(wú)破壞性古油藏顯示,構(gòu)造保存條件較好,滇東北筇竹寺組的油氣遠(yuǎn)景應(yīng)予以關(guān)注。

    參考文獻(xiàn)(References)

    [1] 沈其韓,耿元生,宋會(huì)俠. 加強(qiáng)極低級(jí)變質(zhì)作用研究[J]. 巖石礦物學(xué)雜志,2018,37(2):342-348.[Shen Qihan, Geng Yuansheng,Song Huixia. Strengthening the study of very low-grade metamorphism[J]. Acta Petrologica et Mineralogica, 2018, 37(2):342-348.]

    [2] 湯艷,張?jiān)迄i,齊先茂,等. 很低級(jí)變質(zhì)作用研究及其在沉積盆地中的應(yīng)用[J]. 巖石礦物學(xué)雜志,2015,34(3):353-364.[TangYan, Zhang Yunpeng, Qi Xianmao, et al. A study of the very lowgrademetamorphism and its application to the sedimentary basin[J]. Acta Petrologica et Mineralogica, 2015, 34(3): 353-364.]

    [3] Jiang W T, Peacor D R, Essene E J. Clay minerals in the MacadamsSandstone, California: Implications for substitution of H3O+and H2O and metastability of illite[J]. Clays and Clay Minerals,1994, 42(1): 35-45.

    [4] Essene E J, Peacor D R. Clay mineral thermometry: A critical perspective[J]. Clays and Clay Minerals, 1995, 43(5): 540-553.

    [5] 胡大千,韓春元,馬瑞,等. 內(nèi)蒙古錫林郭勒地區(qū)上古生界極低級(jí)變質(zhì)作用:伊利石和鏡質(zhì)體反射率的證據(jù)[J]. 巖石學(xué)報(bào),2012,28(9):3042-3050.[Hu Daqian, Han Chunyuan, Ma Rui, etal. The very low grade metamorphism in the Upper Paleozoic inXinlingol area of Inner Mongolia, NE China: Evidence from studiesof illite and vitrinite reflectance[J]. Acta Perologica Sinica,2012, 28(9): 3042-3050.]

    [6] Frey M. Very low-grade metamorphism of clastic sedimentaryrocks[M]//Frey M. Low temperature metamorphism. Glasgow:Blackie and Son Ltd. , 1987: 9-58.

    [7] Pollastro R M. The illite/smectite geothermometer-concepts,methodology, and application to basin history and hydrocarbongeneration[M]//Nuccio V F, Barker C E. Applications of thermalmaturity studies to energy exploration. Rocky Mountain Section:SEPM, 1990: 1-18.

    [8] Huang W L, Longo J M, Pevear D R. An experimentally derivedkinetic model for smectite-to illite conversion and its use as a geothermometer[J]. Clays and Clay Minerals, 1993, 41(2): 162-177.

    [9] Frey M, Robinson D. Low-grade metamorphism[M]. Oxford:Blackwell Science, 1999: 10-26.

    [10] Kübler B. La cristallinité de l’llite et les zones tout à faitsupérieures du métamorphism[M]//de la Baconnière. Etages tectoniques,colloque de neuchatel 1966. Neuchatel, Switzerland:Université Neuchátel, à la Baconnière, 1967: 105-121.

    [11] Kübler B. Anchimétamorphisme et schistosité[J]. Bull CentreRech Pau-SNPA, 1967, 1: 259-278.

    [12] Kübler B. Evaluation quantitative du métamorphisme par la cristallinitéde l’illite[J]. Bull Centre Rech Pau-SNPA, 1968, 2:385-397.

    [13] Barker C E, Pawlewicz M J. The correlation of vitrinite reflectancewith maximum temperature in humic organic matter[M]//Buntebarth G, Stegena L. Paleogeothermics: Evaluation of geothermalconditions in the geological past. Berlin: Springer,1986: 79-93.

    [14] Barker C E, Pawlewicz M J. Calculation of vitrinite reflectancefrom thermal histories and peak temperatures: A comparison ofmethods[M]//Mukhopadhyay P K, Dow W G. Vitrinite reflectanceas a maturity parameter: Applications and limitations.Washington: American Chemical Society, 1994: 216-229.

    [15] Stasiuk L D, Lockhart G D, Nassichuk W W, et al. Thermal maturityevaluation of dispersed organic matter inclusions fromkimberlite pipes, Lac de Gras, Northwest Territories, Canada[J].International Journal of Coal Geology, 1999, 40(1): 1-25.

    [16] 楊平,謝淵,汪正江,等. 黔北震旦系燈影組流體活動(dòng)與油氣成藏期次[J]. 石油勘探與開(kāi)發(fā),2014,41(3):313-322,335.[Yang Ping, Xie Yuan, Wang Zhengjiang, et al. Fluid activityand hydrocarbon accumulation period of Sinian Dengying Formationin northern Guizhou, South China[J]. Petroleum Explorationand Development, 2014, 41(3): 313-322, 335.]

    [17] 云南省地質(zhì)礦產(chǎn)局. 云南省區(qū)域地質(zhì)志[M]. 北京:地質(zhì)出版社,1982:18-23.[Yunnan Provincial Bureau of Geology andMineral Resources. Regional geology of Yunnan province[M].Beijing: Geological Publishing House, 1982: 18-23.]

    [18] 劉建清,何利,陳風(fēng)霖,等. 滇東北鹽津地區(qū)中三疊統(tǒng)關(guān)嶺組底部綠豆巖年代學(xué)及地球化學(xué)研究[J]. 巖石學(xué)報(bào),2021,37(7):2245-2255.[Liu Jianqing, He Li, Chen Fenglin, et al. Studieson the chronology and geochemistry of the green pisolites atthe bottom of the Middle Triassic Guanling Formation in Yanjinarea, northeastern Yunnan province[J]. Acta Petrologica Sinica,2021, 37(7): 2245-2255.]

    [19] 劉建清,何利,何平,等. 康滇古陸東緣筇竹寺組地球化學(xué)特征及意義:以云南省昭通市昭陽(yáng)區(qū)鋅廠溝剖面為例[J]. 沉積學(xué)報(bào),2021,39(5):1305-1319.[Liu Jianqing, He Li, He Ping, etal. Geochemical characteristics and significance of the QiongzhusiFormation on the eastern margin of the ancient Kangding-Yunnan Land: Taking the Xinchanggou section of Zhaoyang district,Zhaotong city, Yunnan province as an example[J]. ActaSedimentologica Sinica, 2021, 39(5): 1305-1319.]

    [20] 趙明,陳小明,季峻峰,等. 山東昌濰古近系原型盆地粘土礦物的成巖演化與古地溫[J]. 巖石學(xué)報(bào),2006,22(8):2195-2204.[Zhao Ming, Chen Xiaoming, Ji Junfeng, et al. Diageneticand paleogeothermal evolution of the day minerals in the PaleogeneChangwei prototype basin of Shandong province, China[J]. Acta Petrologica Sinica, 2006, 22(8): 2195-2204.]

    [21] 馮增昭. 沉積巖石學(xué)[M]. 北京:石油工業(yè)出版社,1982:6.[Feng Zengzhao. Sedimentary petrology[M]. Beijing: PetroleumIndustry Press, 1982: 6.]

    [22] 古偉,陳義才,王長(zhǎng)城,等. 川東涪陵地區(qū)長(zhǎng)興組烴源巖地球化學(xué)特征及生烴能力評(píng)價(jià)[J]. 天然氣勘探與開(kāi)發(fā),2013,36(2):23-53.[Gu Wei, Chen Yicai, Wang Changcheng, et al. Geochemicalcharacteristics of source rock in Changxing Formationand its hydrocarbon–generating potential, Fuling area, eastern Sichuan[J]. Natural Gas Exploration amp; Development, 2013, 36(2):23-53.]

    [23] 劉德漢,史繼揚(yáng). 高演化碳酸鹽烴源巖非常規(guī)評(píng)價(jià)方法探討[J]. 石油勘探與開(kāi)發(fā),1994,21(3):113-115.[Liu Dehan, ShiJiyang. Discussion on unconventional evaluation methods ofhighly evolved carbonate source rocks[J]. Petroleum Explorationand Development, 1994, 21(3): 113-115.]

    [24] Nadeau P H, Reynolds Jr C R. Burial and contact metamorphismin the Mancos shale[J]. Clays and Clay Minerals, 1981,29(4): 249-259.

    [25] Hoffman J, Hower J. Clay mineral assemblages as low grademetamorphic geothermometers: Application to the thrust faulteddisturbed belt of Montana, U. S. A. [M]//Schoole P A, SchlugerP R. Aspects of diagenesis. Tulsa: SEPM Society for SedimentaryGeology, 1979: 55-79.

    [26] Jennings S, Thompson G R. Diagenesis of Plio-Pleistocene sedimentsof the Colorado River Delta, southern California[J]. Journalof Sedimentary Research, 1986, 56(1): 89-98.

    [27] Harvey C C, Browne P R L. Mixed-layer clay geothermometryin the Wairakei geothermal field, New Zealand[J]. Clays andClay Minerals, 1991, 39(6): 614-621.

    [28] 王強(qiáng). 粘土礦物地溫計(jì)研究與應(yīng)用[D]. 青島:中國(guó)石油大學(xué)(華東),2004:9.[Wang Qiang. Clay minerals research and application:As geothermometers[D]. Qingdao: China Universityof Petroleum (East China), 2004: 9.]

    [29] Chi Ma, Browne P R L. Alteration mineralogy of sediments inthe Huka Falls Formation of the Te Mihi area, Walrakel[C]//Proceedingsof the 13th Zealand geothermal workshop. Wairakei,1991: 185-191.

    [30] 王丹萍,楊風(fēng)麗,江裕標(biāo),等. 粘土礦物地溫計(jì)在沉積盆地?zé)嵫莼芯恐械膽?yīng)用[J]. 四川地質(zhì)學(xué)報(bào),2010,30(3): 351-355.[Wang Danping, Yang Fengli, Jiang Yubiao, et al. The applicationof clay mineral geothermometer to the study of sedimentarybasin thermal evolution[J]. Sichuan Geological Journal, 2010, 30(3): 351-355.]

    [31] Weaver C E. Possible uses of clay minerals in search for oil[J].Bulletin of the American Association of Petroleum Geologists,1960, 44(9): 1505-1508.

    [32] 畢先梅,莫宣學(xué). 成巖—極低級(jí)變質(zhì)—低級(jí)變質(zhì)作用及有關(guān)礦產(chǎn)[J]. 地學(xué)前緣,2004,11(1):287-294.[Bi Xianmei, MoXuanxue. Transition from diagenesis to low-grade metamorphismand related minerals and energy resources[J]. Earth ScienceFrontiers, 2004, 11(1): 287-294.]

    [33] Chi Ma, Browne P R L, Harvey C C. Crystallinity of subsurfaceclay minerals in the Te Mihi sector of Wairakei geothermal system,New Zealand[C]//Proceedings of the 14th New Zealandgeothermal workshop. Wairakei, 1992: 267-272.

    [34] 陳建平,趙文智,王招明,等. 海相干酪根天然氣生成成熟度上限與生氣潛力極限探討:以塔里木盆地研究為例[J]. 科學(xué)通報(bào),2007,52(增刊1):95-100.[Chen Jianping, Zhao Wenzhi,Wang Zhaoming, et al. Upper limit of maturity and hydrocarbongeneration potential of marine kerogen natural gas: Take theTarim Basin as an example[J]. Chinese Science Bulletin, 2007,52(Suppl. 1): 95-100.]

    [35] Barker C E, Goldstein R H. Fluid-inclusion technique for determiningmaximum temperature in calcite and its comparison tothe vitrinite reflectance geothermometer[J]. Geology, 1990, 18(10): 1003-1006.

    [36] Mullis J, Wolf M, Ferreio M, et al. Temperature determinationthrough fluid inclusion microthermometry and vitrinite reflectancevalues in the Diagenetic-and Anchi–Zones. Conf. Abstr.[C]//EUG XI. Strassbourg, 2001: 30.

    [37] 索書(shū)田,游振東,周漢文. 極低級(jí)變質(zhì)作用和極低級(jí)變質(zhì)帶綜述[J]. 地質(zhì)科技情報(bào),1995,14(1):1-8. [Suo Shutian,You Zhendong, Zhou Hanwen. Very-low-grade metamorphismand metamorphic belt: A review[J]. Geological Science andTechnology Information, 1995, 14(1): 1-8.]

    [38] 袁玉松,馬永生,胡圣標(biāo),等. 中國(guó)南方現(xiàn)今地?zé)崽卣鱗J]. 地球物理學(xué)報(bào),2006,49(4):1118-1126.[Yuan Yusong, Ma Yongsheng,Hu Shengbiao, et al. Present-day geothermal characteristicsin South China[J]. Chinese Journal of Geophysics, 2006, 49(4): 1118-1126.]

    基金項(xiàng)目:中國(guó)地質(zhì)調(diào)查局油氣調(diào)查專(zhuān)項(xiàng)(DD20190080)

    一级a做视频免费观看| 久久 成人 亚洲| av在线播放精品| 曰老女人黄片| 免费黄频网站在线观看国产| 亚洲三级黄色毛片| 插逼视频在线观看| 久久精品国产自在天天线| 亚洲精品国产成人久久av| 日本欧美国产在线视频| 一个人免费看片子| 欧美日韩一区二区视频在线观看视频在线| 久久精品久久久久久噜噜老黄| 成人18禁高潮啪啪吃奶动态图 | 免费观看性生交大片5| 好男人视频免费观看在线| 一区二区三区免费毛片| 国产高清有码在线观看视频| 国产视频首页在线观看| 亚洲精品视频女| 久久久久久久国产电影| 精品99又大又爽又粗少妇毛片| 免费看av在线观看网站| 在线亚洲精品国产二区图片欧美 | 久久97久久精品| 久久久久久久精品精品| 亚洲精品乱久久久久久| 亚洲美女黄色视频免费看| 久久99精品国语久久久| 日日啪夜夜撸| 蜜臀久久99精品久久宅男| 纵有疾风起免费观看全集完整版| 中文天堂在线官网| 亚洲激情五月婷婷啪啪| 亚洲精品久久久久久婷婷小说| 69精品国产乱码久久久| 日本wwww免费看| 十分钟在线观看高清视频www | 99热全是精品| 国产91av在线免费观看| 亚洲第一av免费看| 欧美成人精品欧美一级黄| 99热这里只有是精品在线观看| 插阴视频在线观看视频| 亚洲va在线va天堂va国产| 国产精品久久久久成人av| 精品一区二区三区视频在线| 亚洲国产色片| 国产亚洲一区二区精品| 香蕉精品网在线| 亚洲国产精品国产精品| 国产一区有黄有色的免费视频| 观看免费一级毛片| 欧美变态另类bdsm刘玥| freevideosex欧美| 少妇裸体淫交视频免费看高清| 秋霞在线观看毛片| 免费看不卡的av| 天天躁夜夜躁狠狠久久av| 国产成人精品婷婷| 午夜91福利影院| 亚洲情色 制服丝袜| a级片在线免费高清观看视频| 天堂俺去俺来也www色官网| 久久精品夜色国产| 在线观看www视频免费| 久久久久久人妻| 国产色爽女视频免费观看| 日韩视频在线欧美| 久久毛片免费看一区二区三区| 不卡视频在线观看欧美| 99视频精品全部免费 在线| 一级,二级,三级黄色视频| 深夜a级毛片| 国产成人免费观看mmmm| 搡老乐熟女国产| 亚洲情色 制服丝袜| 蜜桃久久精品国产亚洲av| 国产深夜福利视频在线观看| 久久精品熟女亚洲av麻豆精品| av有码第一页| 亚洲欧洲国产日韩| 亚洲精品国产av成人精品| 观看免费一级毛片| 婷婷色麻豆天堂久久| 久久久国产一区二区| 五月天丁香电影| av在线老鸭窝| 欧美日韩综合久久久久久| 亚洲av成人精品一二三区| 国产高清国产精品国产三级| 欧美 日韩 精品 国产| 亚洲在久久综合| 这个男人来自地球电影免费观看 | 欧美精品人与动牲交sv欧美| 免费人妻精品一区二区三区视频| 视频中文字幕在线观看| 麻豆成人午夜福利视频| 中文字幕久久专区| 日韩人妻高清精品专区| 黑人猛操日本美女一级片| 51国产日韩欧美| 人人妻人人看人人澡| 交换朋友夫妻互换小说| √禁漫天堂资源中文www| 国产精品熟女久久久久浪| tube8黄色片| 大话2 男鬼变身卡| 免费久久久久久久精品成人欧美视频 | 新久久久久国产一级毛片| 一个人看视频在线观看www免费| 亚洲精品aⅴ在线观看| 亚洲av男天堂| 下体分泌物呈黄色| 亚洲三级黄色毛片| www.av在线官网国产| a级一级毛片免费在线观看| 美女福利国产在线| 男女边摸边吃奶| 各种免费的搞黄视频| 黄色怎么调成土黄色| 99久久精品国产国产毛片| 最近的中文字幕免费完整| 亚洲欧美一区二区三区国产| 久久久久久人妻| 欧美变态另类bdsm刘玥| 国产成人精品无人区| 欧美精品一区二区免费开放| 草草在线视频免费看| 国产视频首页在线观看| 下体分泌物呈黄色| 国产精品一二三区在线看| 亚洲成人一二三区av| 久久久a久久爽久久v久久| 尾随美女入室| 亚洲第一区二区三区不卡| av免费在线看不卡| 嫩草影院入口| 在线免费观看不下载黄p国产| 日本91视频免费播放| 在线精品无人区一区二区三| 男女免费视频国产| 日韩精品有码人妻一区| 日本猛色少妇xxxxx猛交久久| 91久久精品国产一区二区成人| 免费黄色在线免费观看| 亚洲图色成人| 欧美变态另类bdsm刘玥| 亚洲国产精品一区二区三区在线| av福利片在线| a级毛片免费高清观看在线播放| 亚洲,一卡二卡三卡| 中文字幕精品免费在线观看视频 | 午夜福利在线观看免费完整高清在| 欧美日韩综合久久久久久| 亚洲精品日韩在线中文字幕| 国产精品国产三级国产av玫瑰| 国国产精品蜜臀av免费| 国产综合精华液| 2021少妇久久久久久久久久久| 欧美xxxx性猛交bbbb| 亚洲图色成人| 色5月婷婷丁香| 欧美日韩一区二区视频在线观看视频在线| 国产欧美日韩综合在线一区二区 | 人妻少妇偷人精品九色| 十八禁高潮呻吟视频 | 日韩强制内射视频| 高清欧美精品videossex| 色视频www国产| 赤兔流量卡办理| 国产黄片美女视频| 日韩伦理黄色片| 免费观看性生交大片5| 极品教师在线视频| 精品人妻熟女毛片av久久网站| 久久99蜜桃精品久久| 9色porny在线观看| 国产日韩欧美视频二区| 少妇的逼水好多| 91午夜精品亚洲一区二区三区| 色吧在线观看| 国产精品一区二区在线不卡| 欧美精品亚洲一区二区| 夜夜爽夜夜爽视频| 在线看a的网站| 午夜福利在线观看免费完整高清在| 18+在线观看网站| 日本黄大片高清| 国产精品国产三级专区第一集| 久久久久久久久大av| 天美传媒精品一区二区| 最近的中文字幕免费完整| 99久久精品一区二区三区| 22中文网久久字幕| 国产视频首页在线观看| 国产欧美日韩一区二区三区在线 | 丝瓜视频免费看黄片| 中文资源天堂在线| 久久这里有精品视频免费| 中文字幕人妻丝袜制服| 国产精品一区二区在线观看99| 日韩一本色道免费dvd| 亚洲av成人精品一二三区| 亚洲av电影在线观看一区二区三区| 性色avwww在线观看| 天堂俺去俺来也www色官网| 国产日韩欧美在线精品| 亚洲国产日韩一区二区| 交换朋友夫妻互换小说| 一区二区av电影网| 18禁在线播放成人免费| 色婷婷久久久亚洲欧美| 内地一区二区视频在线| 亚洲人成网站在线观看播放| 国产熟女欧美一区二区| 国产真实伦视频高清在线观看| 免费不卡的大黄色大毛片视频在线观看| 91在线精品国自产拍蜜月| 成年人午夜在线观看视频| av黄色大香蕉| 国产精品人妻久久久久久| 午夜日本视频在线| 桃花免费在线播放| 国产精品女同一区二区软件| 精品一区在线观看国产| 国产视频首页在线观看| 热99国产精品久久久久久7| 日韩视频在线欧美| 国产有黄有色有爽视频| 久久精品久久久久久噜噜老黄| 国产综合精华液| 久久久久久人妻| 在线播放无遮挡| 亚洲激情五月婷婷啪啪| 国产亚洲午夜精品一区二区久久| 一级a做视频免费观看| 国产高清国产精品国产三级| 大片免费播放器 马上看| 热99国产精品久久久久久7| 亚洲精品中文字幕在线视频 | 久久 成人 亚洲| a级毛色黄片| 欧美变态另类bdsm刘玥| 亚洲美女黄色视频免费看| 精品国产国语对白av| 国产高清三级在线| 午夜福利在线观看免费完整高清在| 我的女老师完整版在线观看| 午夜老司机福利剧场| 亚洲图色成人| 免费av不卡在线播放| 午夜久久久在线观看| 少妇被粗大的猛进出69影院 | 成人亚洲精品一区在线观看| 一区二区av电影网| 18禁在线无遮挡免费观看视频| 另类精品久久| 国产在线视频一区二区| 国产精品无大码| 午夜免费鲁丝| 国产日韩欧美视频二区| 成人综合一区亚洲| 你懂的网址亚洲精品在线观看| 18禁裸乳无遮挡动漫免费视频| 啦啦啦视频在线资源免费观看| 国产欧美日韩一区二区三区在线 | 在线 av 中文字幕| 最新中文字幕久久久久| 中文精品一卡2卡3卡4更新| 国产精品久久久久久av不卡| 国产在视频线精品| 色哟哟·www| 国产亚洲午夜精品一区二区久久| 国产精品三级大全| 久久99热6这里只有精品| 国产午夜精品一二区理论片| 日本色播在线视频| 青春草亚洲视频在线观看| 超碰97精品在线观看| 精品99又大又爽又粗少妇毛片| 午夜免费男女啪啪视频观看| 夜夜看夜夜爽夜夜摸| 各种免费的搞黄视频| 日本爱情动作片www.在线观看| 亚洲精品乱码久久久v下载方式| 免费大片黄手机在线观看| 久久国产乱子免费精品| 久久久久国产精品人妻一区二区| 亚洲欧美成人精品一区二区| 亚洲一区二区三区欧美精品| 女人精品久久久久毛片| 极品人妻少妇av视频| 一二三四中文在线观看免费高清| 极品教师在线视频| 中国美白少妇内射xxxbb| 国产精品成人在线| 少妇的逼好多水| 国产视频首页在线观看| 国内精品宾馆在线| 亚洲不卡免费看| 美女脱内裤让男人舔精品视频| 99热国产这里只有精品6| 日本爱情动作片www.在线观看| 韩国高清视频一区二区三区| 自线自在国产av| 麻豆成人av视频| 最近中文字幕2019免费版| 精品久久久久久久久av| 亚洲精品国产av成人精品| 成年人午夜在线观看视频| 精品一区二区三区视频在线| 最新中文字幕久久久久| 亚洲国产av新网站| 国产永久视频网站| 亚洲怡红院男人天堂| 亚洲国产欧美在线一区| 成人无遮挡网站| 欧美成人精品欧美一级黄| 狂野欧美白嫩少妇大欣赏| 在线观看免费视频网站a站| 国内少妇人妻偷人精品xxx网站| 成人国产av品久久久| 9色porny在线观看| 日日啪夜夜撸| 国产成人精品婷婷| 免费在线观看成人毛片| 赤兔流量卡办理| 3wmmmm亚洲av在线观看| 91精品伊人久久大香线蕉| 亚洲人成网站在线播| 高清视频免费观看一区二区| 黄色日韩在线| 777米奇影视久久| 18禁在线无遮挡免费观看视频| 欧美变态另类bdsm刘玥| 在线观看av片永久免费下载| 亚洲情色 制服丝袜| 香蕉精品网在线| 欧美人与善性xxx| 国产午夜精品久久久久久一区二区三区| av天堂久久9| 午夜影院在线不卡| 人妻夜夜爽99麻豆av| 啦啦啦中文免费视频观看日本| 成年av动漫网址| 夜夜看夜夜爽夜夜摸| 九九久久精品国产亚洲av麻豆| 日韩亚洲欧美综合| 欧美97在线视频| 亚洲av在线观看美女高潮| 成人二区视频| av专区在线播放| 大话2 男鬼变身卡| 黑人猛操日本美女一级片| 下体分泌物呈黄色| 十八禁网站网址无遮挡 | 日韩中文字幕视频在线看片| 免费黄色在线免费观看| 亚洲欧美日韩东京热| 亚洲欧洲精品一区二区精品久久久 | 国产91av在线免费观看| 欧美老熟妇乱子伦牲交| 夜夜看夜夜爽夜夜摸| 国产亚洲最大av| 一级黄片播放器| 午夜影院在线不卡| 成人黄色视频免费在线看| 日日摸夜夜添夜夜爱| 亚洲久久久国产精品| 亚洲伊人久久精品综合| 国产亚洲av片在线观看秒播厂| 精品国产一区二区久久| 天堂俺去俺来也www色官网| 国产精品国产三级专区第一集| 久久久久网色| 久久热精品热| 国产精品99久久99久久久不卡 | 成人综合一区亚洲| 成人午夜精彩视频在线观看| 久久精品国产亚洲av天美| 97超碰精品成人国产| 国产男女超爽视频在线观看| 国产国拍精品亚洲av在线观看| 99久久精品国产国产毛片| 国产毛片在线视频| 亚洲天堂av无毛| 亚洲,一卡二卡三卡| 日本91视频免费播放| 欧美三级亚洲精品| 国产精品久久久久久精品电影小说| 九九在线视频观看精品| 久久久久网色| 欧美三级亚洲精品| av不卡在线播放| 国产成人精品无人区| 一本大道久久a久久精品| 欧美三级亚洲精品| 亚洲va在线va天堂va国产| 亚洲第一区二区三区不卡| 久久青草综合色| 九草在线视频观看| 欧美日韩国产mv在线观看视频| www.av在线官网国产| 精品少妇内射三级| 国产在线免费精品| 亚洲综合精品二区| 人体艺术视频欧美日本| av女优亚洲男人天堂| 性色av一级| 美女大奶头黄色视频| 亚洲精品成人av观看孕妇| 汤姆久久久久久久影院中文字幕| 久久久久国产精品人妻一区二区| 午夜av观看不卡| 国产精品久久久久久久久免| 成人黄色视频免费在线看| 久久久久网色| 久久这里有精品视频免费| 久久久精品免费免费高清| 成人免费观看视频高清| 成人黄色视频免费在线看| 51国产日韩欧美| 久久女婷五月综合色啪小说| 亚洲精品久久午夜乱码| 狠狠精品人妻久久久久久综合| 久久热精品热| 美女cb高潮喷水在线观看| 日韩精品免费视频一区二区三区 | 欧美xxⅹ黑人| 亚洲精品aⅴ在线观看| 日本av手机在线免费观看| 伊人久久国产一区二区| 91久久精品电影网| 国产免费视频播放在线视频| 日韩在线高清观看一区二区三区| 丰满人妻一区二区三区视频av| 日本黄色片子视频| 亚洲av二区三区四区| 久久人人爽人人爽人人片va| 男人狂女人下面高潮的视频| 精品国产一区二区三区久久久樱花| 国产欧美另类精品又又久久亚洲欧美| 青春草国产在线视频| 亚洲人成网站在线观看播放| 观看av在线不卡| 亚洲欧美中文字幕日韩二区| 在线精品无人区一区二区三| 亚洲自偷自拍三级| 免费看日本二区| 色网站视频免费| 看非洲黑人一级黄片| 成人无遮挡网站| 欧美变态另类bdsm刘玥| 在线观看三级黄色| 欧美丝袜亚洲另类| 亚洲av在线观看美女高潮| 亚洲精品自拍成人| 亚洲欧美一区二区三区国产| 人人妻人人澡人人爽人人夜夜| 人人妻人人澡人人看| 亚洲精品成人av观看孕妇| 亚洲一区二区三区欧美精品| 亚洲精品日韩av片在线观看| 久久精品久久久久久久性| 美女主播在线视频| 亚洲人成网站在线播| 最近最新中文字幕免费大全7| 久久久久久久国产电影| 日韩中文字幕视频在线看片| 80岁老熟妇乱子伦牲交| √禁漫天堂资源中文www| 欧美激情国产日韩精品一区| 久久ye,这里只有精品| 欧美激情极品国产一区二区三区 | 免费观看性生交大片5| 久久毛片免费看一区二区三区| 久久精品久久精品一区二区三区| 高清午夜精品一区二区三区| 黄色视频在线播放观看不卡| 国产真实伦视频高清在线观看| av线在线观看网站| 男女边吃奶边做爰视频| 成年人免费黄色播放视频 | 少妇被粗大的猛进出69影院 | 日本午夜av视频| 国产深夜福利视频在线观看| 色哟哟·www| 免费看光身美女| 精品国产国语对白av| 边亲边吃奶的免费视频| 极品教师在线视频| 亚州av有码| 国产成人午夜福利电影在线观看| 亚洲欧美日韩卡通动漫| 精品熟女少妇av免费看| 日本黄色片子视频| 七月丁香在线播放| 亚洲精品一二三| 热re99久久国产66热| 欧美精品一区二区大全| 成人二区视频| 黑人高潮一二区| 国产午夜精品一二区理论片| 国产成人a∨麻豆精品| 亚洲精品aⅴ在线观看| 婷婷色av中文字幕| 丰满人妻一区二区三区视频av| 18+在线观看网站| 热re99久久国产66热| 黄色日韩在线| 精品卡一卡二卡四卡免费| 中国美白少妇内射xxxbb| 高清视频免费观看一区二区| xxx大片免费视频| 精品亚洲乱码少妇综合久久| 日日啪夜夜撸| 久久久久精品性色| 纵有疾风起免费观看全集完整版| 亚洲四区av| 亚洲国产精品一区二区三区在线| 五月天丁香电影| 久久人妻熟女aⅴ| 能在线免费看毛片的网站| 黑丝袜美女国产一区| 久久青草综合色| 久久久久国产网址| 亚洲色图综合在线观看| 日日摸夜夜添夜夜添av毛片| 26uuu在线亚洲综合色| 国产黄色视频一区二区在线观看| 最近中文字幕2019免费版| 成人漫画全彩无遮挡| 两个人免费观看高清视频 | 久久国产亚洲av麻豆专区| 国产探花极品一区二区| 国产亚洲91精品色在线| 国产日韩欧美亚洲二区| 午夜日本视频在线| 国产老妇伦熟女老妇高清| 亚洲经典国产精华液单| 久久精品久久精品一区二区三区| 免费人成在线观看视频色| 日韩伦理黄色片| 国产日韩欧美视频二区| 亚洲第一区二区三区不卡| 精品人妻熟女av久视频| 精华霜和精华液先用哪个| 亚洲国产欧美日韩在线播放 | 亚洲婷婷狠狠爱综合网| 乱码一卡2卡4卡精品| av福利片在线| 久热久热在线精品观看| 国产欧美日韩综合在线一区二区 | 午夜精品国产一区二区电影| 视频区图区小说| 国产精品久久久久久精品古装| h日本视频在线播放| 中文精品一卡2卡3卡4更新| 国产日韩欧美视频二区| 欧美精品一区二区大全| 亚洲va在线va天堂va国产| 嫩草影院新地址| 国产欧美日韩一区二区三区在线 | 国产精品99久久久久久久久| 精品久久国产蜜桃| 久久久精品免费免费高清| 亚州av有码| 青春草视频在线免费观看| 七月丁香在线播放| 又爽又黄a免费视频| 亚洲av.av天堂| 久久久久久人妻| 国产真实伦视频高清在线观看| 欧美3d第一页| 国产一区二区在线观看av| 人人妻人人看人人澡| 夜夜骑夜夜射夜夜干| 观看av在线不卡| 老女人水多毛片| av国产精品久久久久影院| 在线观看av片永久免费下载| 超碰97精品在线观看| 国产精品人妻久久久影院| 午夜免费观看性视频| 久久精品久久精品一区二区三区| 国内少妇人妻偷人精品xxx网站| 亚洲va在线va天堂va国产| 一区在线观看完整版| 久热这里只有精品99| 青春草亚洲视频在线观看| 寂寞人妻少妇视频99o| 熟妇人妻不卡中文字幕| 黑人高潮一二区| 亚洲中文av在线| 性色avwww在线观看| 午夜日本视频在线| 国产精品麻豆人妻色哟哟久久| 精品久久久久久久久av| 精品国产一区二区三区久久久樱花| 男人狂女人下面高潮的视频| 美女国产视频在线观看| 亚洲性久久影院| 男人爽女人下面视频在线观看| 男的添女的下面高潮视频| 新久久久久国产一级毛片| 久久99热这里只频精品6学生| 久久久久久久久久人人人人人人| 国产精品伦人一区二区| 国产亚洲午夜精品一区二区久久| 亚洲国产日韩一区二区| 午夜激情福利司机影院| 五月玫瑰六月丁香| 日日摸夜夜添夜夜添av毛片| 国产精品久久久久久精品电影小说|