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

    黔北震旦系燈影組流體活動(dòng)與油氣成藏期次

    2014-07-25 12:38:28楊平謝淵汪正江李奇艷劉家洪張娣楊智印峰
    石油勘探與開發(fā) 2014年3期
    關(guān)鍵詞:黔北燈影白云石

    楊平,謝淵,汪正江,李奇艷,劉家洪,張娣,楊智,印峰

    (1.國(guó)土資源部成都地質(zhì)礦產(chǎn)研究所;2.國(guó)土資源部沉積盆地與油氣資源重點(diǎn)實(shí)驗(yàn)室;3.中國(guó)石油勘探開發(fā)研究院;4.中國(guó)石化勘探南方分公司)

    0 引言

    黔北地區(qū)及其相鄰的四川盆地東南緣震旦系燈影組目前是南方油氣“下組合”天然氣勘探的重要目的層[1],該地區(qū)燈影組時(shí)代老,且受多期構(gòu)造運(yùn)動(dòng)改造影響,具有成巖和成藏演化史復(fù)雜等特點(diǎn)。多源、多期成藏是中國(guó)南方疊合盆地碳酸鹽巖層系油氣成藏最重要的特征之一[2-3],準(zhǔn)確厘定黔北燈影組油氣成藏期次,深入研究?jī)?chǔ)集層發(fā)育機(jī)理,建立流體活動(dòng)、孔隙演化與油氣成藏之間的耦合關(guān)系,是認(rèn)識(shí)其油氣富集規(guī)律的基礎(chǔ)和尋找優(yōu)質(zhì)儲(chǔ)集層、油氣圈閉的重要手段。

    前人對(duì)黔北燈影組古油藏基本石油地質(zhì)條件進(jìn)行了剖析[4-5],而關(guān)于油氣成藏機(jī)理的研究尚顯不足,存在主要問題為:①流體活動(dòng)期次不清,需建立流體活動(dòng)與儲(chǔ)集層世代膠結(jié)物、溶蝕作用關(guān)系。②需進(jìn)一步研究?jī)?chǔ)集層孔隙演化與各期溶蝕機(jī)制及貢獻(xiàn)。③埋藏史、熱史恢復(fù)缺乏依據(jù)。④目前黔北燈影組儲(chǔ)集層流體包裹體測(cè)溫?cái)?shù)據(jù)少,已知的溫度范圍不能很好地解釋油氣成藏期[6]。

    本文采用流體包裹體、陰極發(fā)光、瀝青反射率、氧同位素分析及薄片觀察等技術(shù)方法,結(jié)合黔北地區(qū)沉積、構(gòu)造演化特點(diǎn),探討黔北燈影組流體幕式活動(dòng)期次及儲(chǔ)集層溶蝕、白云石化機(jī)制與孔隙演化的關(guān)系,在此基礎(chǔ)上分析油氣成藏期次,嘗試恢復(fù)油氣成藏與破壞過程。

    1 石油地質(zhì)背景

    黔北地區(qū)海相碳酸鹽巖層系經(jīng)歷了多期構(gòu)造運(yùn)動(dòng)。不同階段受不同構(gòu)造應(yīng)力場(chǎng)控制,形成了一系列斷裂、褶皺、不整合、古隆起、古斜坡及古巖溶。影響構(gòu)造-沉積格局的關(guān)鍵構(gòu)造事件為都勻運(yùn)動(dòng)、廣西運(yùn)動(dòng)、印支運(yùn)動(dòng)及燕山運(yùn)動(dòng)。都勻運(yùn)動(dòng)發(fā)生于奧陶紀(jì)末,是黔中隆起從水下發(fā)育演變?yōu)殛懮习l(fā)育的轉(zhuǎn)折期[7]。黔中隆起控制黔北地區(qū)及鄰區(qū)震旦系—下古生界油氣的成藏,尤其是加里東—海西期油氣成藏與黔中隆起形成的時(shí)間、范圍等密切相關(guān),形成了麻江奧陶系—志留系古油藏、金沙巖孔[5]和仁懷大灣震旦系燈影組(Zdn)古油藏[8]。金沙巖孔古油藏位于黔中隆起北緣巖孔背斜(見圖1),燈影組上部為一套厚102 m的碳酸鹽臺(tái)地內(nèi)灘相藻屑、砂屑及鮞粒白云巖,靠近頂部發(fā)育厚9.79~20.05 m、有大量瀝青顯示的溶孔型優(yōu)質(zhì)儲(chǔ)集層。瀝青多充填于晶間溶孔、鑄???、粒內(nèi)孔、溶洞及裂縫體系中,儲(chǔ)集層殘余孔隙度為 2.03%~10.85%[9]。金沙巖孔燈影組古油藏是加里東期油氣成藏的重要證據(jù),該古油藏與四川威遠(yuǎn)氣田、慈利南山坪古油藏[10]的成藏條件相似。有機(jī)碳同位素組成及V/Ni值顯示瀝青主要來自下寒武統(tǒng)牛蹄塘組(—C1n)黑色泥頁(yè)巖[9]。

    圖1 黔北燈影組儲(chǔ)集層分布與地層綜合柱狀圖

    2 樣品采集與測(cè)試

    流體包裹體與陰極發(fā)光樣品取自金沙巖孔燈影組上部,采樣剖面位于巖孔綠竹和巖孔箐口,主要巖性為含瀝青鮞粒白云巖、砂屑白云巖及細(xì)晶白云巖,局部含藻及藻屑;并于上述剖面及仁懷大灣等地采集用于常規(guī)薄片分析、固體瀝青及氧同位素分析的樣品,采樣位置及樣品編號(hào)詳見文獻(xiàn)[8-9]。使用 MDS600冷熱臺(tái)儀器(配Carl-Zeiss Axioplan2熒光顯微鏡)進(jìn)行流體包裹體的顯微測(cè)溫和測(cè)鹽,測(cè)溫誤差為±0.1 ℃。使用 CL8200MK5型陰極發(fā)光儀進(jìn)行陰極發(fā)光分析,工作束電壓為10 kV,束電流約為500 A。瀝青反射率檢測(cè)儀器為MPV-SP顯微鏡光度計(jì)。使用MAT253穩(wěn)定同位素比質(zhì)譜儀進(jìn)行氧同位素測(cè)定,精度為±0.2‰。

    3 流體包裹體特征

    金沙巖孔燈影組白云巖中各期次白云石、方解石及自生石英中可見數(shù)量不等的流體包裹體,主要類型為氣液兩相鹽水包裹體,從14個(gè)薄片(總共15個(gè)制片樣品)中獲得了包裹體均一溫度數(shù)據(jù),在 7個(gè)薄片中獲得了14組鹽度及冰點(diǎn)資料(見表1)。這些包裹體主要賦存于較大的自形白云石及自生石英中,少量存于粗晶方解石中,大小為4.1~26.8 μm,氣液比主要為10%,部分為5%。在101個(gè)測(cè)溫?cái)?shù)據(jù)中,石英中有70個(gè),白云石中有24個(gè),方解石中的包裹體較小,很難發(fā)現(xiàn)適用于測(cè)溫的較大包裹體。

    表1 金沙巖孔燈影組儲(chǔ)集層流體包裹體測(cè)溫及測(cè)鹽數(shù)據(jù)

    3.1 白云石內(nèi)的包裹體

    鏡下觀察初步劃分兩類共4期白云石,第1類為微晶和鑲嵌狀半自形微晶白云石,為第1期;第2類為孔洞充填的白云石,根據(jù)成巖序列、白云石晶粒大小等可劃分為3期,即分別為第2、3、4期。第1期白云石中很難發(fā)現(xiàn)可供測(cè)溫的包裹體,孔洞中充填的2、3、4期白云石中流體包裹體發(fā)育,以氣液兩相鹽水包裹體為主,氣液比為5%~10%,根據(jù)流體包裹體的產(chǎn)狀、均一溫度、鹽度及密度可將包裹體劃分為 3期(見圖2)。

    圖2 金沙巖孔燈影組各期膠結(jié)物中流體包裹體均一溫度、鹽度及密度分布

    第Ⅰ期包裹體在沿孔洞邊緣或裂縫分布的第 2期白云石中發(fā)育,白云石顆粒較細(xì),表面較臟,包裹體個(gè)體較小,多呈長(zhǎng)方形、圓狀及次圓狀。5個(gè)可供測(cè)溫的包裹體大小為 4.5~7.9 μm,均一溫度為 101.7~111.4 ℃,平均106.3 ℃,平均密度為0.95 g/cm3。

    第Ⅱ期包裹體賦存于沿第 2期白云石邊緣分布的第3期白云石中,個(gè)別沿裂縫分布,白云石結(jié)晶粗大,包裹體較發(fā)育,多呈圓狀及次圓狀,串珠狀分布,大小為4.1~12.8 μm(共10個(gè)測(cè)定數(shù)據(jù)),均一溫度范圍為126.1~148.9 ℃,平均138.3 ℃,平均密度為0.94 g/cm3,個(gè)別包裹體的鹽度最大可達(dá) 16.7%,為爆發(fā)式流體活動(dòng)的特征。

    第Ⅲ期包裹體主要發(fā)育在沿第2期及第3期白云石邊緣分布的粗晶白云石中,也常見于大型溶孔充填的粗—巨晶白云石中,即賦存于第 4期白云石中,包裹體多呈長(zhǎng)條、次圓及不規(guī)則狀,成群分布,包裹體大小為 6.8~22.1 μm(共 9個(gè)測(cè)定數(shù)據(jù)),均一溫度為178.4~225.1 ℃,平均 206.8 ℃,平均密度 0.86 g/cm3。

    3.2 石英內(nèi)的包裹體

    石英內(nèi)的流體包裹體以氣液兩相為主,氣液比為10%,亦可劃分為3期(見圖2)。

    石英內(nèi)的第Ⅰ期包裹體與白云石中的第Ⅰ期包裹體同期形成,多呈長(zhǎng)方形、圓狀及次圓狀,沿裂隙或成群分布(見圖3a),包裹體大小為6.9~12.3 μm(共5個(gè)測(cè)定數(shù)據(jù)),鹽水包裹體均一溫度為87.1~110.8 ℃,平均密度為0.96 g/cm3。

    圖3 燈影組各期膠結(jié)物中流體包裹體特征

    石英中的第Ⅱ期包裹體與白云石中的第Ⅱ期包裹體同期形成,多呈圓狀及次圓狀串珠、成群或沿裂隙分布(見圖3b),大小為4.7~23.7 μm(共36個(gè)測(cè)定數(shù)據(jù)),均一溫度范圍為130.5~163.0 ℃,平均145.8 ℃,平均密度為0.93 g/cm3。

    石英中的第Ⅲ期包裹體與白云石中的第Ⅲ期包裹體同期形成,多呈圓狀及次圓狀,成群或串珠狀分布(見圖3c),包裹體大小7.1~26.8 μm(共29個(gè)測(cè)定數(shù)據(jù)),其均一溫度略低于同期白云石,為 166.9~218.5 ℃,平均199.1 ℃,平均密度為0.86 g/cm3。

    3.3 方解石內(nèi)的包裹體

    方解石內(nèi)的包裹體數(shù)據(jù)均來自薄片YQ-2-1(見圖3d),樣品中裂隙被亮晶方解石充填成脈。亮晶方解石中流體包裹體較發(fā)育,以氣液兩相鹽水包裹體為主,氣液比為 10%,多成群分布,均一溫度為 95.3~116.4 ℃,平均106.1 ℃(共7個(gè)測(cè)定數(shù)據(jù)),冰點(diǎn)為?3.8~?3.0 ℃(2個(gè)測(cè)定數(shù)據(jù)),鹽度(NaCl質(zhì)量分?jǐn)?shù))為4.9%~6.1%。方解石中低溫與低鹽度包裹體既可能形成于石油初始成藏階段,也可能形成于晚期油氣藏抬升破壞階段。多期流體活動(dòng)或白云石化作用會(huì)改造早期形成的方解石,現(xiàn)今裂隙中充填的方解石應(yīng)是晚期氣藏抬升破壞過程中儲(chǔ)集層流體釋放與地表淡水進(jìn)入后形成的,其形成時(shí)間晚于白云石及石英中的包裹體,相當(dāng)于第Ⅳ期。

    4 地?zé)崽荻扰c熱史

    本次研究采用等效鏡質(zhì)體反射率作為主要古溫標(biāo),另外通過分析燈影組儲(chǔ)集層包裹體均一溫度的分布特征,結(jié)合其在關(guān)鍵構(gòu)造運(yùn)動(dòng)時(shí)刻的埋藏深度和與之對(duì)應(yīng)的包裹體均一溫度來恢復(fù)黔北地區(qū)熱史。埋藏史曲線中地層厚度采用遵義幅1∶200 000地質(zhì)調(diào)查報(bào)告中的震旦系—中生界各組厚度,金之鈞等[11]的研究表明黔中隆起晚侏羅世—早白堊世(J3—K1)構(gòu)造抬升時(shí)刻為距今97 Ma,因此距今97 Ma以來的構(gòu)造抬升史可參考四川盆地東南部丁山 1井埋藏史曲線[12]。由于金沙巖孔上侏羅統(tǒng)—下白堊統(tǒng)已遭剝蝕,其原始地層厚度采用鄰區(qū)桐梓幅及丁山 1井相關(guān)資料,主要依據(jù)如下:①川東南—黔中隆起晚侏羅世—早白堊世構(gòu)造抬升時(shí)刻為距今97 Ma,表明黔北地區(qū)在晚侏羅世—早白堊世仍維持湖相沉積。②遵義幅、桐梓幅及丁山1井所在的綦江幅中下侏羅統(tǒng)厚度變化較小,且具有相似的巖性及沉積環(huán)境。

    瀝青反射率(Rb)測(cè)試結(jié)果表明,金沙巖孔固體瀝青的平均瀝青反射率(Rba)為2.95%~3.86%(6個(gè)薄片),仁懷大灣瀝青平均瀝青反射率為2.95%~3.56%(2個(gè)薄片)。根據(jù)公式Ro= 0.668Rb+ 0.346[13]換算得到金沙巖孔固體瀝青等效鏡質(zhì)體反射率(Ro)值為2.32%~2.92%,平均為2.64%;仁懷大灣固體瀝青Ro值為2.32%~2.72%,平為2.52%。上述資料表明黔北燈影組瀝青熱演化程度較高,儲(chǔ)集層曾經(jīng)歷了高溫?zé)嵫莼饔谩G备鞯貐^(qū)燈影組儲(chǔ)集層瀝青等效鏡質(zhì)體反射率普遍低于上覆牛蹄塘組(見表 2),這可能是由于燈影組儲(chǔ)集層油氣成藏過程中形成的超壓抑制了有機(jī)質(zhì)演化[14]。這種超壓現(xiàn)象也普遍存在于四川威遠(yuǎn)燈影組氣藏的演化過程中[15],因此本研究采用金沙巖孔牛蹄塘組烴源巖平均等效鏡質(zhì)體反射率(Roa=3.77%)作為恢復(fù)古地溫的主要指標(biāo),該值與周邊各剖面牛蹄塘組等效鏡質(zhì)體反射率非常接近,數(shù)據(jù)可信度較高。

    表2 黔北地區(qū)燈影組與牛蹄塘組瀝青反射率

    應(yīng)用Ro值恢復(fù)最高古地溫(或古地?zé)崽荻龋┑姆椒ɑ蚬捷^多,本次研究選用了 3種常用方法:① Karweil-Teichmuliler圖解法,即利用溫度、Ro、有效恒溫時(shí)間 3者間的關(guān)系恢復(fù)最高古地溫。從金沙巖孔地層沉積演化史來看,其有效恒溫時(shí)間約為50 Ma,在距今約97 Ma時(shí)達(dá)到最高溫度,之后開始迅速抬升剝蝕。利用該方法得到晚侏羅世—早白堊世最高溫度為222 ℃。②Hood等[18]認(rèn)為烴源巖成熟度或Ro值主要取決于所經(jīng)歷的最高古地溫和在不低于最高古地溫15 ℃范圍內(nèi)的受熱時(shí)間,提出了用有機(jī)質(zhì)變質(zhì)標(biāo)尺和有效受熱時(shí)間來衡量最高古地溫的方法。根據(jù)Hood圖版及有效受熱時(shí)間(60 Ma)得到最大古地溫為270 ℃。③根據(jù) Shibaoka等[19]提出的古地溫計(jì)算公式:Tc=(lgRo+0.87?0.149lgt0)/0.004 5,計(jì)算得到古地溫為231 ℃,較利用Hood圖版恢復(fù)的最高古地溫低,與利用Karwei1-Teichmuliler圖解法估算的溫度較接近。

    由表 1可見,金沙巖孔孔洞充填白云石中流體包裹體最高均一溫度為225.1 ℃,因此認(rèn)為采用Shibaoka公式或Karweil-Teichmuliler圖解法恢復(fù)古地溫相對(duì)適中。設(shè)地表平均溫度為20 ℃,估算J3—K1(最大埋藏深度7 660 m)對(duì)應(yīng)的地?zé)崽荻葹?.75 ℃/100 m。

    將Ⅰ、Ⅱ、Ⅲ期包裹體均一溫度分別投影在埋藏史-熱史圖中,發(fā)現(xiàn)包裹體均一溫度可以和加里東中晚期(O1—S1)、印支期(P2—T2)及燕山早期(J1—J2)燈影組埋深形成較好的對(duì)應(yīng)關(guān)系,反映 3期埋藏增溫引起的生烴及流體活動(dòng)。3期流體包裹體均一溫度并非連續(xù)分布,而是在特定的溫度范圍出現(xiàn)一定的“斷點(diǎn)”,如第Ⅱ期與第Ⅲ期之間缺乏163.0~166.9 ℃的數(shù)據(jù),這是構(gòu)造抬升造成生烴停滯和流體活動(dòng)急劇減少導(dǎo)致的,因此5 000 m深度對(duì)應(yīng)的溫度約為165 ℃,對(duì)應(yīng)的地?zé)崽荻燃s為2.90 ℃/100 m,加里東中晚期古地?zé)崽荻却蠹s為3.18 ℃/100 m,表明古地?zé)崽荻嚷猿尸F(xiàn)降低的趨勢(shì)(見圖4),現(xiàn)今地?zé)崽荻冉抵?.71 ℃/100 m(丁山1井的實(shí)測(cè)值)。

    圖4 黔北金沙震旦系古油氣藏?zé)崾芳把莼^程

    5 膠結(jié)作用

    根據(jù)陰極發(fā)光原理,結(jié)合幕式流體活動(dòng)特點(diǎn)與鏡下礦物結(jié)構(gòu),按地質(zhì)時(shí)間的先后關(guān)系,將黔北燈影組儲(chǔ)集層膠結(jié)物劃分為6個(gè)世代(見圖5)。第1世代膠結(jié)物為顆粒之間的粉晶白云石,后期被桐灣期表生巖溶部分改造,但仍見連片不規(guī)則殘留,陰極發(fā)光較昏暗,與海水Mn2+含量低有關(guān),表明成巖流體與海水關(guān)系密切。第 2世代膠結(jié)物為桐灣期表生巖溶作用后沉淀的方解石,淺埋過程中發(fā)生重結(jié)晶及白云石化,但保留殘余孔洞及縫隙,為后期有機(jī)酸溶蝕及油氣充注創(chuàng)造了條件。陰極發(fā)光為亮橙色,具環(huán)帶狀結(jié)構(gòu),表明高 Mn2+、Fe2+含量及淡水成因特征。第 3世代膠結(jié)物為圍繞或包圍固體瀝青的粉—細(xì)晶白云石。牛蹄塘組烴源巖在低熟階段形成大量有機(jī)酸及液態(tài)烴進(jìn)入燈影組,有機(jī)酸主要對(duì)第1世代及第2世代膠結(jié)物進(jìn)行溶蝕并重新沉淀結(jié)晶,形成以白云石、自生石英為主的膠結(jié)物,所捕獲的液態(tài)烴及鹽水包裹體均一溫度為87.1~111.4 ℃。受儲(chǔ)集層中桐灣期殘留淡水稀釋作用影響,白云石及方解石內(nèi) Mn2+、Fe2+的含量及陰極發(fā)光強(qiáng)度介于第1世代及第2世代之間,包裹體平均鹽度僅為5.47%。第4世代膠結(jié)物在成因、結(jié)構(gòu)及成分上與第 3世代類似,主要受生烴高峰階段的大量流體活動(dòng)影響而形成,儲(chǔ)集層溫度較第3世代高30~40 ℃。該世代膠結(jié)物主要為白云石及石英,通過對(duì)此前膠結(jié)物的部分交代、溶解與結(jié)晶形成。在少數(shù)較粗大粒間孔中,晶體明亮粗大,大小一般為0.1~2.0 mm,且白云石在偏光顯微鏡下具波狀消光,陰極發(fā)光較強(qiáng),多呈亮橙黃色,可能與深埋藏環(huán)境有機(jī)酸選擇性溶蝕有關(guān)。第 5世代膠結(jié)物形成于深埋環(huán)境下。隨著溫度增加,儲(chǔ)集層中液態(tài)烴裂解,重組分形成的高溫固態(tài)焦瀝青充填部分孔隙,輕組分裂解演化為濕氣。另外在牛蹄塘組烴源巖生氣階段,含氣烴流體進(jìn)入儲(chǔ)集層殘余孔隙。因此,受儲(chǔ)集層液態(tài)烴裂解和烴源巖生氣階段流體活動(dòng)的影響,形成了具有高溫成因的白云石及石英。白云石結(jié)晶粗大,大小一般為0.1~4.0 mm,流體包裹體具有最大的均一溫度與鹽度。第 6世代膠結(jié)物是晚期方解石及硅質(zhì),形成于氣藏破壞階段。由于喜馬拉雅早期構(gòu)造運(yùn)動(dòng)的抬升作用,儲(chǔ)集層溫度與壓力逐漸下降,儲(chǔ)集層流體中的SiO2發(fā)生沉淀[20],同時(shí)淡水進(jìn)入儲(chǔ)集層形成方解石。野外露頭及鏡下觀察發(fā)現(xiàn)部分張裂縫中方解石結(jié)晶粗大,陰極發(fā)光較強(qiáng),呈亮橙色,條帶狀硅質(zhì)普遍沿裂縫充填,且均無瀝青充填,反映了氣藏破壞階段的流體活動(dòng)特征。

    圖5 黔北燈影組儲(chǔ)集層世代膠結(jié)物與多期溶孔

    6 溶蝕作用

    第 1期為震旦紀(jì)晚期桐灣期表生巖溶作用。大氣淡水對(duì)尚未完全固結(jié)成巖的沉積物(藻、砂屑顆粒、藻紋層等)中的不穩(wěn)定礦物如文石、高鎂方解石等進(jìn)行組構(gòu)選擇性溶蝕,形成鑄??住⒘?nèi)溶孔及溶洞,隨后部分被亮晶方解石充填,部分殘余孔隙被保留。

    第2~3期溶蝕作用為埋藏環(huán)境下幕式流體活動(dòng)帶來的有機(jī)酸非組構(gòu)選擇性溶蝕作用,形成了粒內(nèi)、粒間、晶間溶孔及溶縫和溶洞等,并伴隨瀝青的充填。

    第 4期溶蝕作用形成超大孔隙、粒內(nèi)粒間溶孔、晶間溶孔、溶縫和溶洞。其孔隙空間干凈、無瀝青,部分被粗晶白云石及石英充填,表明第 4期溶蝕作用發(fā)生于天然氣充注階段。埋藏環(huán)境下有機(jī)酸溶蝕及自生石英沉淀的化學(xué)方程式為[21]:2KAlSi3O8(鉀長(zhǎng)石)+2H++H2O→AlSiO5(OH)4(高嶺石)+4SiO2+2K+。下寒武統(tǒng)黑色泥巖中普遍含有一定量的鉀長(zhǎng)石、斜長(zhǎng)石及黏土礦物[2],在酸性條件下很容易發(fā)生長(zhǎng)石及黏土礦物溶蝕作用,并形成SiO2沉淀,因此當(dāng)含SiO2的酸性流體進(jìn)入儲(chǔ)集層時(shí)即發(fā)生石英的沉淀,隨后酸性流體造成了早期的碳酸鹽膠結(jié)物重溶與結(jié)晶。

    有機(jī)酸溶蝕(2~4期溶蝕)主要受控于幕式流體活動(dòng),各期溶蝕與充填作用分別發(fā)生于淺埋藏深度段(2 100~2 900 m)、中等埋藏深度段(3 200~5 000 m)及深埋藏深度段(5 200~7 600 m)。各期溶蝕作用之后發(fā)生相應(yīng)的充填作用,并表現(xiàn)為多期次溶蝕-充填復(fù)合疊加,埋藏環(huán)境下各期有機(jī)酸溶蝕和油氣充注主要圍繞早期巖溶殘留的孔洞進(jìn)行,均表現(xiàn)為有機(jī)酸溶蝕—油氣充注—充填3個(gè)過程,這3個(gè)過程幾乎同時(shí)進(jìn)行且可以互相佐證。

    7 白云石特征

    第 1期白云石包括菱型微晶和鑲嵌狀半自形微晶白云石,相當(dāng)于第 1世代膠結(jié)物,陰極發(fā)光較弱且昏暗,常保留原始沉積物的構(gòu)造特點(diǎn),屬于準(zhǔn)同生海底成巖環(huán)境產(chǎn)物。第2期白云石呈粉晶—細(xì)晶,自形—半自形,屬淺—中埋藏環(huán)境產(chǎn)物,相當(dāng)于第3世代膠結(jié)物,包裹體均一溫度為101.7~111.4 ℃,在加里東中晚期烴類流體首次充注儲(chǔ)集層時(shí)形成。第3期白云石呈細(xì)—中晶,多呈自形或半自形,有潔凈環(huán)邊,有時(shí)還具有環(huán)帶構(gòu)造,常破壞原巖組構(gòu),屬中—較深埋藏環(huán)境產(chǎn)物,相當(dāng)于第4世代膠結(jié)物,包裹體均一溫度為126.1~148.9 ℃,在印支期烴類流體第2次充注儲(chǔ)集層時(shí)形成。第4期白云石呈中—粗晶,多為自形晶,與第2~3期白云石不同,晶體充填的孔隙普遍較大,孔壁干凈且無瀝青充填,屬深埋藏環(huán)境產(chǎn)物,相當(dāng)于第5世代膠結(jié)物,包裹體均一溫度為178.4~225.1 ℃,形成于早—中侏羅世深埋藏高溫演化階段。

    依據(jù)Vasconcelos等[22]提出的白云石-水氧同位素溫度分餾方程(1 000 lnαd-w=2.73×106T2+0.26)可以估算白云石形成時(shí)的古溫度(見表 3)。該方程式中最重要的參數(shù)是由白云石氧同位素及白云石形成時(shí)古流體氧同位素計(jì)算的分餾系數(shù)。Veizer等[23]依據(jù)古生代腕足化石及前寒武紀(jì)化石中氧同位素資料,提出古生代海水中氧同位素組成顯著低于現(xiàn)代海水,并隨著地質(zhì)年代變新逐漸變輕。Wallmann等[24]依據(jù)上述資料推算了寒武紀(jì)—現(xiàn)今海水氧同位素組成的變化情況,認(rèn)為震旦紀(jì)晚期—早寒武世(距今約550~540 Ma)5 000 m與3 000 m深度海水的δ18O值(SMOW)約為?7.6‰和?6.2‰。研究區(qū)燈影組的沉積環(huán)境為碳酸鹽臺(tái)地或淺灘,處于弱蒸發(fā)環(huán)境,蒸發(fā)性海水較正常海水具有較高的δ18O值[25],因此本次研究假設(shè)黔北燈影組沉積時(shí)古海水的δ18O值約為?5‰。

    表3 黔北地區(qū)燈影組白云巖(白云石及方解石)氧同位素組成及白云石形成溫度估算

    充填白云石包含了2~4世代膠結(jié)物,而3~4世代白云石膠結(jié)物為埋藏作用產(chǎn)物,與之進(jìn)行氧同位素交換的不再是海水而是地層水,因此本文采用四川威遠(yuǎn)氣田地層水氧同位素主要分布范圍(3.2‰~5.9‰)[27]的初始值進(jìn)行計(jì)算。盡管以3.2‰作為燈影組埋藏階段儲(chǔ)集層地層水的氧同位素組成仍然具有一定的主觀性,但白云石化流體的氧同位素組成每變化1‰,只能導(dǎo)致溫度計(jì)算結(jié)果約 5~8 ℃的變化[28]。對(duì)于最大古埋深接近8 000 m、最大古溫度超過220 ℃的黔北燈影組(見圖4),這一溫度偏差仍然在可以接受的范圍內(nèi),換言之,白云石化流體氧同位素組成取值的局部偏差對(duì)溫度計(jì)算結(jié)果并沒有明顯的影響。

    目前,關(guān)于低溫和高溫白云石化環(huán)境的溫度界限仍然存在爭(zhēng)議。如果以80 ℃[29]作為低溫和高溫白云石化環(huán)境的界限,黔北地區(qū)燈影組中以第 1世代白云石為主的致密粉晶白云巖的形成溫度為 19.7~40.5 ℃(見表 3),因此第 1世代白云石的形成溫度應(yīng)該在20 ℃左右,屬于近地表環(huán)境低溫白云石化流體成因,應(yīng)為原生或準(zhǔn)同生白云石。

    充填白云石的形成溫度與其流體包裹體均一溫度具有非常相似的分布范圍。例如仁懷大灣充填白云石的形成溫度為87.6~109.2 ℃(見表3),與白云石第1期包裹體溫度101.7~111.4 ℃相近;金沙巖孔箐口充填白云石的形成溫度為127.9~156.7 ℃,與白云石第2期包裹體溫度126.1~148.9 ℃基本相同。

    含瀝青溶孔白云巖由于既具有低溫成因的第 1世代白云石,同時(shí)溶孔中大量充填了后期中—深埋藏白云石,因此δ18O值無論取?5‰還是 3.2‰,所估算的溫度均為第 1世代白云石與充填白云石的中間值,也進(jìn)一步說明第 1世代白云石與埋藏階段充填白云石具有不同的形成溫度。

    8 油氣成藏演化過程

    黔北燈影組儲(chǔ)集層充填物流體包裹體均一溫度、6世代膠結(jié)物的特征及次序、4期溶蝕作用及4期白云石揭示了4個(gè)含烴流體活動(dòng)期,表明存在3個(gè)油氣成藏期和1個(gè)氣藏破壞期。

    第 1期為淺埋藏油氣初次充注階段(距今 470~428 Ma)。牛蹄塘組烴源巖進(jìn)入生烴門限,含烴流體在黔中隆起的控制下,沿?zé)粲敖M與牛蹄塘組間的不整合面由隆起的北斜坡向構(gòu)造高部位運(yùn)聚成藏。富含液態(tài)烴、有機(jī)酸及SiO2的流體進(jìn)入桐灣期遭受淡水溶蝕及白云石(方解石)充填后殘余的孔洞,有機(jī)酸溶蝕擴(kuò)大了殘余孔洞,溶蝕對(duì)象包括第1、2世代白云石。該過程中形成了第 3世代白云石與石英,賦存的流體包裹體的均一溫度和鹽度較低,但密度較高。

    第2期為中等埋深油氣二次充注階段(距今252~228 Ma)。牛蹄塘組烴源巖進(jìn)入二次生油階段,處于輕質(zhì)油—濕氣階段,儲(chǔ)集層中第 1期充注的液態(tài)烴還未大量裂解,形成輕質(zhì)油氣藏。在第1次生油過程中,干酪根消耗了大部分產(chǎn)生有機(jī)酸的官能團(tuán),因此二次生油過程中有機(jī)酸的含量急劇減少。富含液態(tài)烴及SiO2的流體主要溶蝕孔隙中的第2、3世代白云石,并形成了第 4世代白云石與石英。賦存的流體包裹體均一溫度較高,鹽度變化大,密度適中。

    第 3期為深埋天然氣充注階段(距今 177~145 Ma)。牛蹄塘組烴源巖隨著埋藏深度及溫度的增加相繼進(jìn)入濕氣及干氣演化階段,同時(shí)儲(chǔ)集層中液態(tài)烴已開始大量裂解成氣,形成氣藏及固體瀝青。富含氣態(tài)烴及SiO2的流體主要溶蝕充填于孔隙中的第2、4世代白云石,形成了超大孔及第5世代粗晶白云石與石英,賦存的流體包裹體均一溫度和鹽度最高,但密度低。

    除上述油氣成藏階段的 3幕流體活動(dòng)外,喜馬拉雅早期也存在 1個(gè)流體活動(dòng)期,期間儲(chǔ)集層烴類、流體逸散,與淡水一起沿構(gòu)造裂縫發(fā)生滲透。在野外地質(zhì)露頭及鏡下均可見沿裂縫先后充填的硅質(zhì)條帶與方解石,根據(jù)方解石形成溫度可以估算天然氣藏破壞的深度。依據(jù)周根陶等[30]提出的氧同位素分餾方程(1 000 lnαc-w=20.6×103/T?34.71)及劉子琦等[26]獲得的貴州中西部大氣降水δ18O值,估算淡水方解石形成的溫度為41.2~86.3 ℃,其形成深度為782~2 446 m。根據(jù)方解石中賦存流體包裹體的均一溫度(95.3~116.4 ℃)及從丁山 1井獲取的井下地?zé)崽荻龋?.71℃/100 m),估算出天然氣藏遭到破壞的深度約為2 779~3 557 m。限于采樣位置及方解石中包裹體數(shù)據(jù)個(gè)數(shù),天然氣藏的實(shí)際破壞溫度與深度范圍可能比估算值寬。研究表明天然氣藏埋深約3 600 m時(shí)(距今42 Ma)開始破壞,儲(chǔ)集層流體釋放,約2 800 m時(shí)淡水進(jìn)入儲(chǔ)集層,天然氣藏已破壞殆盡。因此研究認(rèn)為3 600 m可能是黔北地區(qū)燈影組天然氣藏保存的臨界深度,2 800 m為氣藏保存的最淺深度。

    勘探實(shí)踐表明,丁山1井(燈影組頂埋深3 490 m)燈影組下部(4 578.0~4 603.0 m)水體總礦化度為253.00~333.37 g/L,均為 CaCl2型,表明該井段處于水文地質(zhì)封閉環(huán)境。寒武系陡坡寺組(相當(dāng)于黔北高臺(tái)組)—清虛洞組(2 792~2 819 m)水體礦化度較低,為 7.35~9.61 g/L,為 Na2SO4型,該井段的保存條件曾遭受過破壞。林1井(燈影組頂界埋深2 580 m)燈影組白云巖2 799.18~2 866.55 m層段水體礦化度為29.9 g/L,為NaHCO3型,表明該處燈影組儲(chǔ)集層受到了大氣降水下滲的影響,地層水處于交替停滯帶,保存條件較差。綜上所述,本文提出3 600 m的氣藏保存理想埋深與實(shí)際勘探結(jié)論基本吻合。埋深大于3 600 m時(shí),保存條件相對(duì)較好,水型一般為CaCl2型,埋深小于3 600 m時(shí),保存條件遭受不同程度的破壞。

    9 結(jié)論

    黔北燈影組儲(chǔ)集層中發(fā)育 4期流體包裹體,均一溫度分別為 87.1~111.4 ℃、126.1~163.0 ℃、166.9~225.1 ℃和 95.3~116.4 ℃,反映了兩期油氣充注、1期天然氣充注和晚期天然氣藏破壞過程。利用實(shí)測(cè)瀝青反射率、包裹體均一溫度估算出黔北地區(qū)早古生代—中生代古地?zé)崽荻燃s為3.18~2.75 ℃/100 m。

    燈影組儲(chǔ)集層中發(fā)育6個(gè)世代的膠結(jié)物和4期白云石,經(jīng)歷了4期溶蝕作用。不同埋深環(huán)境下的3期流體活動(dòng)帶來的有機(jī)酸對(duì)早期古巖溶的殘余孔隙及充填物具有復(fù)合疊加溶蝕作用。第 1世代白云石為近地表低溫準(zhǔn)同生成因,第3~5世代膠結(jié)物、第2~4期溶蝕作用、第2~4期白云石及第1~3期包裹體分別對(duì)應(yīng)于3幕流體活動(dòng)與油氣成藏。

    燈影組經(jīng)歷了3個(gè)油氣成藏期,即加里東中晚期(距今470~428 Ma)、印支期(距今252~228 Ma)和燕山早期(距今 177~145 Ma)。燈影組氣藏遭受破壞的時(shí)間約為距今42 Ma,深度約為2 779~3 557 m,3 600 m埋深可能是黔北地區(qū)燈影組天然氣藏保存的理想深度。

    致謝:成都地質(zhì)礦產(chǎn)研究所王劍、余謙及丘東洲研究員為本研究提出了很好的意見,中國(guó)石化無錫石油地質(zhì)研究所等單位在樣品測(cè)試方面給予了支持,在此一并表示感謝。

    符號(hào)注釋:

    Rb——瀝青反射率,%;Rba——平均瀝青反射率,%;Ro——瀝青等效鏡質(zhì)體反射率,%;Roa——平均瀝青等效鏡質(zhì)體反射率,%;Tc——烴源巖經(jīng)歷的最高古地溫,℃;t0——烴源巖熱演化時(shí)間,Ma;αd-w——白云石-水氧同位素分餾系數(shù);T——白云石形成時(shí)的古溫度,℃;αc-w——方解石-水氧同位素分餾系數(shù)。

    [1]魏國(guó)齊, 沈平, 楊威, 等.四川盆地震旦系大氣田形成條件與勘探遠(yuǎn)景區(qū)[J].石油勘探與開發(fā), 2013, 40(2): 129-138.Wei Guoqi, Shen Ping, Yang Wei, et al.Formation conditions and exploration prospects of Sinian large gas fields, Sichuan Basin[J].Petroleum Exploration and Development, 2013, 40(2): 129-138.

    [2]金之鈞.中國(guó)海相碳酸鹽巖層系油氣形成與富集規(guī)律[J].中國(guó)科學(xué): D輯: 地球科學(xué), 2011, 41(7): 910-926.Jin Zhijun.Formation and accumulation of oil and gas in marine carbonate sequences in Chinese sedimentary basins[J].Science China:Series D: Earth Sciences, 2012, 55(3): 368-385.

    [3]趙文智, 沈安江, 胡素云, 等.中國(guó)碳酸鹽巖儲(chǔ)集層大型化發(fā)育的地質(zhì)條件與分布特征[J].石油勘探與開發(fā), 2012, 39(1): 1-12.Zhao Wenzhi, Shen Anjiang, Hu Suyun, et al.Geological conditions and distributional features of large-scale carbonate reservoirs onshore China[J].Petroleum Exploration and Development, 2012,39(1): 1-12.

    [4]王津義, 涂偉, 曾華盛, 等.黔西北地區(qū)天然氣成藏地質(zhì)特征[J].石油實(shí)驗(yàn)地質(zhì), 2008, 30(5): 445-449.Wang Jinyi, Tu Wei, Zeng Huasheng, et al.Geologic characteristics of natural gas accumulations in the northwest of Guizhou[J].Petroleum Geology & Experiment, 2008, 30(5): 445-449.

    [5]趙澤桓, 張桂權(quán), 薛秀麗.黔中隆起下組合古油藏和殘余油氣藏[J].天然氣工業(yè), 2008, 28(8): 39-42.Zhao Zehuan, Zhang Guiquan, Xue Xiuli.Fossil oil pools and residual oil and gas pools in the lower assemblage of Qianzhong uplift[J].Natural Gas Industry, 2008, 28(8): 39-42.

    [6]陶樹, 湯達(dá)禎, 李鳳, 等.黔中隆起北緣金沙巖孔古油藏特征及成藏期次厘定[J].中國(guó)礦業(yè)大學(xué)學(xué)報(bào), 2009, 38(4): 576-581.Tao Shu, Tang Dazhen, Li Feng, et al.Paleo-oil reservoir characteristics and hydrocarbon accumulation period redefinition of Jinsha Yankong in north margin of the central Guizhou uplift[J].Journal of China University of Mining & Technology, 2009, 38(4):576-581.

    [7]戎嘉余, 陳旭, 王懌, 等.奧陶-志留紀(jì)之交黔中古陸的變遷: 證據(jù)與啟示[J].中國(guó)科學(xué): D輯: 地球科學(xué), 2011, 41(10): 1407-1415.Rong Jiayu, Chen Xu, Wang Yi, et al.Northward expansion of Central Guizhou Oldland through the Ordovician and Silurian transition: Evidence and implications[J].Science China: Series D:Earth Sciences, 2011, 41(10): 1407-1415.

    [8]楊平, 汪正江, 賀永忠, 等.貴州仁懷縣震旦系燈影組古油藏成藏條件及油氣地質(zhì)意義[J].地質(zhì)通報(bào), 2012, 31(11): 1894-1901.Yang Ping, Wang Zhengjiang, He Yongzhong, et al.The reservoir forming condition and petroleum geological significance of Dengying Formation paleo-reservoirs in Renhuai, Guizhou[J].Geological Bulletin of China, 2012, 31(11): 1894-1901.

    [9]楊平, 謝淵, 汪正江, 等.金沙巖孔燈影組古油藏瀝青有機(jī)地球化學(xué)特征及油源分析[J].地球化學(xué), 2012, 41(5): 452-465.Yang Ping, Xie Yuan, Wang Zhengjiang, et al.Geochemical characteristics and oil source correlation of Dengying formation paleo-reservoir in Jinsha[J].Geochimica, 2012, 41(5): 452-465.

    [10]趙宗舉, 馮加良, 陳學(xué)時(shí), 等.湖南慈利燈影組古油藏的發(fā)現(xiàn)及意義[J].石油與天然氣地質(zhì), 2001, 22(2): 114-119.Zhao Zongju, Feng Jialiang, Chen Xueshi, et al.Discovery of Dengying formation fossil pool in Cili, Hunan and its significance[J].Oil & Gas Geology, 2001, 22(2): 114-119.

    [11]金之鈞, 袁玉松, 劉全有, 等.J3-K1構(gòu)造事件對(duì)南方海相源蓋成藏要素的控制作用[J].中國(guó)科學(xué): D輯: 地球科學(xué), 2012, 42(12):1791-1801.Jin Zhijun, Yuan Yusong, Liu Quanyou, et al.Controls of Late Jurassic-Early Cretaceous tectonic event on source rocks and seals in marine sequences, South China[J].Science China: Series D: Earth Sciences, 2013, 56(2): 228-239.

    [12]袁海鋒.四川盆地震旦系—下古生界油氣成藏機(jī)理[D].成都: 成都理工大學(xué), 2008: 70-82.Yuan Haifeng.The mechanism of hydrocarbon accumulation of Sinian-Lower Palaeozoic, Sichuan Basin[D].Chengdu: Chengdu University of Technology, 2008: 70-82.

    [13]劉德漢, 史繼揚(yáng).高演化碳酸鹽烴源巖非常規(guī)評(píng)價(jià)方法探討[J].石油勘探與開發(fā), 1994, 21(3): 113-115.Liu Dehan, Shi Jiyang.High evolution of carbonate source rocks unconventional evaluation method[J].Petroleum Exploration and Development, 1994, 21(3): 113-115.

    [14]郝芳, 鄒華耀, 方勇, 等.超壓環(huán)境下有機(jī)質(zhì)熱演化和生烴作用機(jī)理[J].石油學(xué)報(bào), 2006, 27(5): 9-17.Hao Fang, Zou Huayao, Fang Yong, et al.Kinetics of organic matter maturation and hydrocarbon generation in overpressure environment[J].Acta Petrolei Sinica, 2006, 27(5): 9-17.

    [15]袁海鋒, 徐國(guó)盛, 王國(guó)芝, 等.川中地區(qū)震旦系油氣成藏過程的相態(tài)演化與勘探前景[J].成都理工大學(xué)學(xué)報(bào): 自然科學(xué)版, 2009,36(6): 662-668.Yuan Haifeng, Xu Guosheng, Wang Guozhi, et al.Phase evolution during hydrocarbon accumulation and exploration prospect for Sinian reservoirs in Central Sichuan Basin, China[J].Journal of Chengdu University of Technology: Science & Technology Edition,2009, 36(6): 662-668.

    [16]楊劍, 易發(fā)成, 錢壯志.黔北下寒武統(tǒng)黑色巖系古地溫及其指示意義[J].礦物學(xué)報(bào), 2009, 29(1): 87-94.Yang Jian, Yi Facheng, Qian Zhuangzhi.Study of paleotemperature of the lower Cambrian black shale series and their implications in northern Guizhou province[J].Acta Mineralogica Sinica, 2009, 29(1):87-94.

    [17]壇俊穎, 王文龍, 王延斌, 等.中上揚(yáng)子下寒武統(tǒng)牛蹄塘組海相烴源巖評(píng)價(jià)[J].海洋地質(zhì)前沿, 2011, 27(3): 23-27.Tan Junying, Wang Wenlong, Wang Yanbin, et al.Marine source rock assessment for lower Cambrian Niutitang Formation in middle and upper Yangtze region[J].Marine Geology Frontiers, 2011, 27(3):23-27.

    [18]Hood A, Gutjahr C C M, Heacock R L.Organic metamorphism and the generation of petroleum[J].AAPG Bulletin, 1975, 59(6):986-996.

    [19]Shibaoka M, Saxby J D, Taylor G H.Hydrocarbon generation in Gippsland Basin, Australia: Comparison with Cooper Basin,Australia[J].AAPG Bulletin, 1978, 62(7): 1151-1158.

    [20]周利敏, 張德會(huì), 席斌斌.巖石中的滲透率、流體流動(dòng)及熱液成礦作用[J].地學(xué)前緣, 2008, 15(3): 299-310.Zhou Limin, Zhang Dehui, Xi Binbin.Rock permeability, fluid flow and hydrothermal ore-forming processes[J].Earth Science Frontiers,2008, 15(3): 299-310.

    [21]Worden R H, Barclay S A.Internally-sourced quartz cement due to externally-derived CO2in sub-arkosic sandstones, North Sea[J].Journal of Geochemical Exploration, 2000, 69/70: 645-649.

    [22]Vasconcelos C, McKenzie J A, Warthmann R, et al.Calibration of theδ18O paleothermometer for dolomite precipitated in microbial cultures and natural environments[J].Geology, 2005, 33(4): 317-320.

    [23]Veizer J, Ala D, Azmy K, et al.87Sr/86Sr,δ13C andδ18O evolution of Phanerozoic seawater[J].Chemical Geology, 1999, 161(1/2/3): 59-88.

    [24]Wallmann K.The geological water cycle and the evolution of marineδ18O values[J].Geochimica et Cosmochimica Acta, 2001, 65(15):2469-2485.

    [25]Sp?tl C, Burns S J.Formation of18O-depleted dolomite within a marine evaporitic sequence, Triassic Reichenhall Formation,Austria[J].Sedimentology, 1991, 38(6): 1041-1057.

    [26]劉子琦, 李紅春, 徐曉梅, 等.貴州中西部洞穴水系與碳酸鈣的穩(wěn)定同位素意義[J].地質(zhì)論評(píng), 2007, 53(2): 233-241.Liu Ziqi, Li Hongchun, Xu Xiaomei, et al.Stable isotopes of water and carbonate samples from caves in central western Guizhou:Implications of paleoclimate and paleoenvironment[J].Geological Review, 2007, 53(2): 233-241.

    [27]胡作維, 黃思靜, 李志明, 等.川東北地區(qū)三疊系飛仙關(guān)組白云化流體溫度[J].中國(guó)科學(xué): 地球科學(xué), 2012, 42(12): 1817-1829.Hu Zuowei, Huang Sijing, Li Zhiming, et al.Temperatures of dolomitizing fluids in the Feixianguan Formation from the Northeastern Sichuan Basin[J].Science China: Earth Sciences, 2012,42(12): 1817-1829.

    [28]Allan J R, Wiggins W D.Dolomite reservoirs: Geochemical techniques for evaluating origin and distribution[M].Tulsa: AAPG,1993: 19-32.

    [29]Machel H G.Concepts and models of dolomitization: A critical reappraisal[C]//Braithwaite C J R, Rizzi G, Darke G.The geometry and petrogenesis of dolomite hydrocarbon reservoirs.London:Geological Society of London? 2004: 7-63.

    [30]周根陶, 鄭永飛.碳酸鈣-水體系氧同位素分餾系數(shù)的低溫實(shí)驗(yàn)研究[J].地學(xué)前緣, 2000, 7(2): 321-338.Zhou Gentao, Zheng Yongfei.Experimental studies of oxygen isotope fractionations between CaCO3and H2O at low temperatures[J].Earth Science Frontiers, 2000, 7(2): 321-338.

    猜你喜歡
    黔北燈影白云石
    水泥回轉(zhuǎn)窯用白云石質(zhì)耐火材料的研究進(jìn)展
    工業(yè)爐(2023年6期)2023-03-29 05:29:18
    遵化市魏家井白云石礦冶金用白云巖礦地質(zhì)特征
    白云石有序度與流變特征的研究進(jìn)展
    半盞——第八話:燈影下
    黔北山鄉(xiāng)水稻飄香
    元稹與燈影牛肉
    川中震旦系燈影組儲(chǔ)集層形成及演化研究
    青陽(yáng)白云石礦采場(chǎng)回采順序優(yōu)化
    黔北鄉(xiāng)間——尋覓時(shí)光的味道
    黔北煤田官田壩向斜煤層氣賦存及潛力分析
    大型av网站在线播放| 亚洲第一青青草原| 免费高清视频大片| 51午夜福利影视在线观看| 成年人黄色毛片网站| 成熟少妇高潮喷水视频| www.精华液| 国产欧美日韩综合在线一区二区| 又大又爽又粗| 90打野战视频偷拍视频| 人人妻人人澡人人看| 一级毛片女人18水好多| 12—13女人毛片做爰片一| 亚洲国产看品久久| 啦啦啦 在线观看视频| 一二三四在线观看免费中文在| 久久久精品欧美日韩精品| 老汉色av国产亚洲站长工具| 男女做爰动态图高潮gif福利片 | 国产精品久久久av美女十八| 免费看美女性在线毛片视频| 亚洲欧美精品综合久久99| 国产亚洲精品av在线| 亚洲精品国产区一区二| 亚洲精品国产精品久久久不卡| 国产亚洲av嫩草精品影院| 国产精品永久免费网站| 亚洲三区欧美一区| 亚洲男人天堂网一区| 91国产中文字幕| 国产人伦9x9x在线观看| 精品国产亚洲在线| 19禁男女啪啪无遮挡网站| 后天国语完整版免费观看| 国产片内射在线| 国产精品1区2区在线观看.| 男女下面进入的视频免费午夜 | 欧美成人一区二区免费高清观看 | 十八禁网站免费在线| 久久久久九九精品影院| 一级毛片女人18水好多| 久久欧美精品欧美久久欧美| 免费看美女性在线毛片视频| 91麻豆av在线| 女性被躁到高潮视频| 国产男靠女视频免费网站| 露出奶头的视频| 一级作爱视频免费观看| 91在线观看av| 97人妻精品一区二区三区麻豆 | 久久精品91无色码中文字幕| 免费久久久久久久精品成人欧美视频| 又紧又爽又黄一区二区| 免费av毛片视频| 人人澡人人妻人| 亚洲专区字幕在线| 美女免费视频网站| 日本 欧美在线| 欧美成人免费av一区二区三区| 国产精品乱码一区二三区的特点 | 欧美大码av| 亚洲欧美激情在线| 在线天堂中文资源库| 成年人黄色毛片网站| 我的亚洲天堂| 亚洲国产欧美日韩在线播放| 嫩草影院精品99| 日本撒尿小便嘘嘘汇集6| 波多野结衣巨乳人妻| 亚洲 国产 在线| 亚洲精品美女久久av网站| 国产精华一区二区三区| 91精品国产国语对白视频| 国产高清视频在线播放一区| 国产精品自产拍在线观看55亚洲| 性欧美人与动物交配| 亚洲一区二区三区不卡视频| 日本 av在线| 老鸭窝网址在线观看| 欧美国产日韩亚洲一区| 亚洲欧美日韩另类电影网站| 国产成年人精品一区二区| 欧美成狂野欧美在线观看| 精品久久久久久久久久免费视频| 身体一侧抽搐| 琪琪午夜伦伦电影理论片6080| 欧美日韩一级在线毛片| 精品高清国产在线一区| 成人av一区二区三区在线看| 在线观看www视频免费| 很黄的视频免费| 在线观看午夜福利视频| 亚洲狠狠婷婷综合久久图片| 久久精品影院6| 国内毛片毛片毛片毛片毛片| 一级片免费观看大全| 中文字幕人妻熟女乱码| 婷婷丁香在线五月| 欧美成人午夜精品| 国产亚洲欧美在线一区二区| 免费不卡黄色视频| 久热爱精品视频在线9| 欧美老熟妇乱子伦牲交| 人人妻人人澡人人看| 国产亚洲精品av在线| 国产精品久久久av美女十八| 大香蕉久久成人网| 黄片播放在线免费| 国产亚洲精品第一综合不卡| 男女下面插进去视频免费观看| 亚洲国产精品久久男人天堂| 亚洲成人久久性| 精品久久久久久久毛片微露脸| 亚洲一区高清亚洲精品| 9热在线视频观看99| 亚洲成人国产一区在线观看| 老司机午夜十八禁免费视频| 在线观看66精品国产| 亚洲久久久国产精品| 日本一区二区免费在线视频| 老司机深夜福利视频在线观看| 国产午夜福利久久久久久| 婷婷精品国产亚洲av在线| 亚洲精品一区av在线观看| 亚洲国产看品久久| 亚洲男人天堂网一区| 18禁美女被吸乳视频| 后天国语完整版免费观看| 日本三级黄在线观看| 欧美不卡视频在线免费观看 | 国产亚洲精品久久久久久毛片| 99国产精品免费福利视频| 国产成人av教育| 搡老熟女国产l中国老女人| 丁香欧美五月| 麻豆av在线久日| 精品一区二区三区av网在线观看| 国产av在哪里看| 69av精品久久久久久| 国产人伦9x9x在线观看| 午夜福利高清视频| 亚洲国产精品成人综合色| 国产高清videossex| 亚洲精品国产区一区二| 俄罗斯特黄特色一大片| 九色国产91popny在线| 午夜视频精品福利| 性色av乱码一区二区三区2| 亚洲第一电影网av| av天堂在线播放| 久久久久久久精品吃奶| 91精品国产国语对白视频| 国产成人系列免费观看| 中国美女看黄片| 欧美中文综合在线视频| 99精品久久久久人妻精品| 男女下面进入的视频免费午夜 | 成人永久免费在线观看视频| 一本久久中文字幕| 女警被强在线播放| 中国美女看黄片| 日本撒尿小便嘘嘘汇集6| 一级毛片女人18水好多| 老司机在亚洲福利影院| 久久狼人影院| 一级,二级,三级黄色视频| 50天的宝宝边吃奶边哭怎么回事| 亚洲成人精品中文字幕电影| 中国美女看黄片| 国产精品一区二区在线不卡| 亚洲色图综合在线观看| 最近最新中文字幕大全免费视频| 一个人观看的视频www高清免费观看 | 亚洲欧美日韩高清在线视频| 欧美日韩精品网址| 亚洲精品粉嫩美女一区| 久久伊人香网站| 国产真人三级小视频在线观看| 久久精品成人免费网站| 91精品国产国语对白视频| 国产又色又爽无遮挡免费看| 在线观看日韩欧美| 激情在线观看视频在线高清| 一个人免费在线观看的高清视频| xxx96com| 91成人精品电影| 日本 欧美在线| 久久 成人 亚洲| 99re在线观看精品视频| x7x7x7水蜜桃| 999久久久精品免费观看国产| 亚洲性夜色夜夜综合| 亚洲国产精品久久男人天堂| 久久人人精品亚洲av| 1024视频免费在线观看| 99久久99久久久精品蜜桃| av视频免费观看在线观看| 麻豆国产av国片精品| 黑人操中国人逼视频| 黄色女人牲交| e午夜精品久久久久久久| 九色亚洲精品在线播放| 在线天堂中文资源库| 成人亚洲精品av一区二区| 午夜两性在线视频| 亚洲一码二码三码区别大吗| 亚洲av片天天在线观看| 十八禁网站免费在线| 亚洲第一欧美日韩一区二区三区| 18禁观看日本| 国产精华一区二区三区| 午夜福利影视在线免费观看| 在线天堂中文资源库| 美女扒开内裤让男人捅视频| 国产男靠女视频免费网站| 国产成人影院久久av| 亚洲国产日韩欧美精品在线观看 | 男人的好看免费观看在线视频 | 怎么达到女性高潮| 精品国产一区二区三区四区第35| 日韩一卡2卡3卡4卡2021年| 亚洲专区国产一区二区| 两个人视频免费观看高清| 日韩大尺度精品在线看网址 | 两个人看的免费小视频| 精品国产一区二区三区四区第35| 欧美色欧美亚洲另类二区 | 精品一区二区三区四区五区乱码| 侵犯人妻中文字幕一二三四区| 色综合站精品国产| 欧美中文综合在线视频| 婷婷精品国产亚洲av在线| 久久久久九九精品影院| 中文字幕精品免费在线观看视频| 亚洲无线在线观看| 国产99久久九九免费精品| 亚洲片人在线观看| 在线视频色国产色| www.www免费av| 美国免费a级毛片| 欧美黑人欧美精品刺激| 高清在线国产一区| 99re在线观看精品视频| 最新在线观看一区二区三区| 身体一侧抽搐| 99riav亚洲国产免费| 美国免费a级毛片| 十八禁网站免费在线| 亚洲avbb在线观看| 97人妻精品一区二区三区麻豆 | 中文字幕色久视频| 人人妻,人人澡人人爽秒播| 桃红色精品国产亚洲av| 亚洲专区国产一区二区| 免费观看人在逋| 夜夜躁狠狠躁天天躁| 午夜福利欧美成人| 国产精品亚洲av一区麻豆| 老司机靠b影院| 母亲3免费完整高清在线观看| 热re99久久国产66热| 在线免费观看的www视频| 色综合欧美亚洲国产小说| 久久久久久人人人人人| 在线观看免费日韩欧美大片| 亚洲欧美精品综合一区二区三区| 99国产精品一区二区蜜桃av| 国产成+人综合+亚洲专区| 19禁男女啪啪无遮挡网站| 亚洲精品国产色婷婷电影| 激情视频va一区二区三区| 人人妻人人爽人人添夜夜欢视频| 欧美日韩亚洲综合一区二区三区_| 男人的好看免费观看在线视频 | 亚洲人成77777在线视频| 精品第一国产精品| 国产免费av片在线观看野外av| 国产欧美日韩综合在线一区二区| 欧美日本视频| 男人舔女人的私密视频| 好看av亚洲va欧美ⅴa在| 国产精品久久久av美女十八| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品自产拍在线观看55亚洲| 国产欧美日韩综合在线一区二区| 无限看片的www在线观看| 男人舔女人的私密视频| 性色av乱码一区二区三区2| 9色porny在线观看| 国产三级在线视频| а√天堂www在线а√下载| 午夜久久久久精精品| 日本五十路高清| 夜夜爽天天搞| 国产1区2区3区精品| 日韩成人在线观看一区二区三区| 视频在线观看一区二区三区| 亚洲中文日韩欧美视频| 在线视频色国产色| 一区在线观看完整版| 又黄又粗又硬又大视频| 国产精品久久久人人做人人爽| 亚洲色图 男人天堂 中文字幕| 美国免费a级毛片| 欧美黄色片欧美黄色片| 三级毛片av免费| 9热在线视频观看99| 久久久水蜜桃国产精品网| 亚洲成a人片在线一区二区| 又黄又爽又免费观看的视频| 男女做爰动态图高潮gif福利片 | 亚洲中文日韩欧美视频| 妹子高潮喷水视频| 男女床上黄色一级片免费看| 国产国语露脸激情在线看| 欧美日韩瑟瑟在线播放| 精品人妻在线不人妻| 可以在线观看的亚洲视频| 国产亚洲精品一区二区www| 人成视频在线观看免费观看| 欧美中文日本在线观看视频| 美女国产高潮福利片在线看| 美女高潮喷水抽搐中文字幕| tocl精华| 国产亚洲精品一区二区www| 精品一品国产午夜福利视频| 日本 欧美在线| 久久国产亚洲av麻豆专区| 此物有八面人人有两片| 欧美激情高清一区二区三区| 国产精品 欧美亚洲| 97碰自拍视频| a在线观看视频网站| 欧美性长视频在线观看| 18禁国产床啪视频网站| 女警被强在线播放| 中文亚洲av片在线观看爽| 两个人视频免费观看高清| 最近最新中文字幕大全电影3 | 在线观看www视频免费| 91av网站免费观看| 精品一区二区三区av网在线观看| 亚洲专区中文字幕在线| 一夜夜www| 久久久久精品国产欧美久久久| 国语自产精品视频在线第100页| 亚洲av成人av| 女人被狂操c到高潮| 亚洲成人免费电影在线观看| 在线免费观看的www视频| av超薄肉色丝袜交足视频| 国产精品电影一区二区三区| 亚洲国产精品sss在线观看| 91国产中文字幕| 老司机深夜福利视频在线观看| 欧美黑人欧美精品刺激| 免费在线观看黄色视频的| 亚洲欧美激情在线| 国产av又大| 久久精品91蜜桃| 日韩大码丰满熟妇| 欧美日韩黄片免| 给我免费播放毛片高清在线观看| 妹子高潮喷水视频| 精品熟女少妇八av免费久了| 精品电影一区二区在线| 日韩一卡2卡3卡4卡2021年| 久久久久亚洲av毛片大全| 精品无人区乱码1区二区| 欧美+亚洲+日韩+国产| 免费在线观看黄色视频的| 午夜福利高清视频| 极品人妻少妇av视频| 人人妻,人人澡人人爽秒播| 很黄的视频免费| 国产av在哪里看| 国产亚洲欧美98| 伊人久久大香线蕉亚洲五| 热99re8久久精品国产| 亚洲 国产 在线| 国产成人一区二区三区免费视频网站| 国产男靠女视频免费网站| 女人精品久久久久毛片| 欧美av亚洲av综合av国产av| 一卡2卡三卡四卡精品乱码亚洲| 巨乳人妻的诱惑在线观看| 精品午夜福利视频在线观看一区| 男女下面进入的视频免费午夜 | 精品第一国产精品| 亚洲精品美女久久av网站| 窝窝影院91人妻| 国产三级黄色录像| 在线观看舔阴道视频| 欧美激情高清一区二区三区| 亚洲av五月六月丁香网| 99riav亚洲国产免费| 在线观看www视频免费| 欧美成狂野欧美在线观看| 黄色 视频免费看| 欧美日韩一级在线毛片| 免费看a级黄色片| 亚洲成国产人片在线观看| 一级a爱片免费观看的视频| 大码成人一级视频| 在线永久观看黄色视频| 国内精品久久久久精免费| 在线天堂中文资源库| 男男h啪啪无遮挡| 国产成人精品在线电影| 天天一区二区日本电影三级 | 久久精品91蜜桃| 欧美色视频一区免费| 精品电影一区二区在线| 亚洲一区高清亚洲精品| 亚洲片人在线观看| 黄片播放在线免费| 亚洲无线在线观看| 天天一区二区日本电影三级 | 成人av一区二区三区在线看| 日韩一卡2卡3卡4卡2021年| 久久国产乱子伦精品免费另类| 老司机福利观看| av视频免费观看在线观看| 久久热在线av| 亚洲无线在线观看| 一级毛片精品| 99精品欧美一区二区三区四区| 满18在线观看网站| av超薄肉色丝袜交足视频| 美女扒开内裤让男人捅视频| 久久久久亚洲av毛片大全| 午夜免费激情av| 国产精品免费视频内射| 国产成人精品久久二区二区91| 国产精品亚洲美女久久久| 亚洲狠狠婷婷综合久久图片| 91精品国产国语对白视频| 亚洲精品美女久久av网站| 九色国产91popny在线| 久久国产乱子伦精品免费另类| 欧美日韩福利视频一区二区| 免费在线观看影片大全网站| 99riav亚洲国产免费| 久久人妻av系列| 中文字幕人妻丝袜一区二区| 国产99白浆流出| 这个男人来自地球电影免费观看| 不卡一级毛片| 极品人妻少妇av视频| 啦啦啦免费观看视频1| 国产片内射在线| 亚洲av成人av| 成人永久免费在线观看视频| 一级,二级,三级黄色视频| 999久久久精品免费观看国产| 99久久国产精品久久久| 午夜福利欧美成人| 久久午夜综合久久蜜桃| 国产精品日韩av在线免费观看 | 露出奶头的视频| 少妇熟女aⅴ在线视频| 亚洲精华国产精华精| 亚洲精品在线美女| 久久这里只有精品19| 亚洲国产精品sss在线观看| 欧美国产日韩亚洲一区| 亚洲精品粉嫩美女一区| 久久精品亚洲熟妇少妇任你| 高清毛片免费观看视频网站| 黑人巨大精品欧美一区二区蜜桃| 午夜福利成人在线免费观看| 日本精品一区二区三区蜜桃| 精品久久久精品久久久| 亚洲欧美一区二区三区黑人| 国产精品亚洲美女久久久| 久久久水蜜桃国产精品网| 丝袜美足系列| 精品国产美女av久久久久小说| 亚洲成国产人片在线观看| 韩国精品一区二区三区| 99在线视频只有这里精品首页| 免费无遮挡裸体视频| 制服丝袜大香蕉在线| 国产精品电影一区二区三区| 18禁国产床啪视频网站| 日韩欧美在线二视频| 免费少妇av软件| 国产成人一区二区三区免费视频网站| 丝袜在线中文字幕| 国产色视频综合| 一边摸一边做爽爽视频免费| 中文字幕人妻熟女乱码| 在线视频色国产色| 亚洲情色 制服丝袜| 国产亚洲欧美在线一区二区| 欧美不卡视频在线免费观看 | videosex国产| 老司机在亚洲福利影院| 在线观看免费午夜福利视频| 美女午夜性视频免费| 91老司机精品| 少妇裸体淫交视频免费看高清 | 国产精品一区二区在线不卡| 久久久水蜜桃国产精品网| 亚洲中文日韩欧美视频| a级毛片在线看网站| 少妇粗大呻吟视频| 午夜福利成人在线免费观看| 天堂√8在线中文| 欧美 亚洲 国产 日韩一| 一级片免费观看大全| ponron亚洲| 国产亚洲欧美精品永久| 变态另类丝袜制服| 亚洲一区中文字幕在线| 在线观看日韩欧美| 日韩精品免费视频一区二区三区| 成人欧美大片| 叶爱在线成人免费视频播放| 黄网站色视频无遮挡免费观看| aaaaa片日本免费| 丰满的人妻完整版| 看黄色毛片网站| 国产精品久久久久久人妻精品电影| 国产一区二区三区视频了| 亚洲成国产人片在线观看| 免费在线观看完整版高清| 美国免费a级毛片| 午夜免费激情av| 亚洲美女黄片视频| 久久影院123| 亚洲狠狠婷婷综合久久图片| 国产视频一区二区在线看| 性少妇av在线| 亚洲,欧美精品.| 久久精品亚洲精品国产色婷小说| 在线观看免费日韩欧美大片| 亚洲国产毛片av蜜桃av| 日韩大尺度精品在线看网址 | 国产99白浆流出| 成年女人毛片免费观看观看9| 88av欧美| 少妇熟女aⅴ在线视频| 麻豆国产av国片精品| 中文字幕人成人乱码亚洲影| 人妻丰满熟妇av一区二区三区| 午夜亚洲福利在线播放| 老汉色∧v一级毛片| 人妻久久中文字幕网| 亚洲熟妇熟女久久| 精品福利观看| 午夜老司机福利片| 一进一出好大好爽视频| 好看av亚洲va欧美ⅴa在| 午夜福利视频1000在线观看 | 9色porny在线观看| 亚洲精品中文字幕一二三四区| 12—13女人毛片做爰片一| 免费久久久久久久精品成人欧美视频| 色综合站精品国产| 国产视频一区二区在线看| 十分钟在线观看高清视频www| 男女床上黄色一级片免费看| 久久久久国产精品人妻aⅴ院| 在线观看免费午夜福利视频| 国产免费av片在线观看野外av| 在线观看舔阴道视频| 成人手机av| 波多野结衣av一区二区av| 50天的宝宝边吃奶边哭怎么回事| 大型av网站在线播放| 嫩草影视91久久| 黄色丝袜av网址大全| 国产精品精品国产色婷婷| 亚洲 欧美 日韩 在线 免费| 自拍欧美九色日韩亚洲蝌蚪91| 国产av一区二区精品久久| 久久精品国产综合久久久| 国产色视频综合| av在线播放免费不卡| 午夜日韩欧美国产| 黄片大片在线免费观看| 久久国产精品男人的天堂亚洲| 国产成人免费无遮挡视频| 中亚洲国语对白在线视频| 午夜影院日韩av| 久久青草综合色| 大型av网站在线播放| 精品国产超薄肉色丝袜足j| 欧美日韩亚洲国产一区二区在线观看| 午夜福利免费观看在线| 日韩中文字幕欧美一区二区| 老司机福利观看| 美女高潮到喷水免费观看| 欧美不卡视频在线免费观看 | a级毛片在线看网站| 国产亚洲精品第一综合不卡| 免费在线观看日本一区| 午夜福利在线观看吧| 国产三级在线视频| 国产激情欧美一区二区| 麻豆av在线久日| 成人三级黄色视频| 久久国产亚洲av麻豆专区| 高清毛片免费观看视频网站| 国产一区在线观看成人免费| 午夜日韩欧美国产| 如日韩欧美国产精品一区二区三区| 在线观看免费视频网站a站| 欧美乱码精品一区二区三区| 69精品国产乱码久久久| 国产精品久久久久久精品电影 | 国产不卡一卡二| 韩国av一区二区三区四区| 69精品国产乱码久久久|