左銀輝， 張旺， 李兆影， 李佳蔚， 郝情情， 胡杰

1 “油氣藏地質與開發(fā)工程國家重點實驗室”,成都理工大學， 成都 610059 2 成都理工大學能源學院， 成都 610059 3 中國石油川慶鉆探工程有限公司地質勘探開發(fā)研究院， 成都 610051 4 中國科學院地質與地球物理研究所，巖石圈演化實驗室， 北京 100029 5 中國冶金地質總局礦產(chǎn)資源研究院， 北京 100025

?

查干凹陷中、新生代構造-熱演化史

左銀輝1,2， 張旺1,2， 李兆影3， 李佳蔚4， 郝情情5， 胡杰2

1 “油氣藏地質與開發(fā)工程國家重點實驗室”,成都理工大學， 成都 610059 2 成都理工大學能源學院， 成都 610059 3 中國石油川慶鉆探工程有限公司地質勘探開發(fā)研究院， 成都 610051 4 中國科學院地質與地球物理研究所，巖石圈演化實驗室， 北京 100029 5 中國冶金地質總局礦產(chǎn)資源研究院， 北京 100025

查干凹陷是中國內蒙古銀根-額濟納旗盆地中最具有勘探潛力的凹陷.為了揭示其構造-熱演化歷史，本文利用35個磷灰石裂變徑跡和119個鏡質體反射率數(shù)據(jù)，采用耦合反演的方法恢復了查干凹陷白堊紀以來的熱歷史.結果顯示查干凹陷白堊系具有高的古地溫梯度，并且查干凹陷經(jīng)歷了地溫梯度快速增加階段(K1b-K1s)，地溫梯度高峰階段(K1y)，高地溫延續(xù)階段(K2w)和熱沉降階段(Cz)四個構造-熱演化階段.此外，基于熱史恢復結果，正演獲得查干凹陷高的古地溫梯度有利于烴源巖成熟生烴，早白堊世的高古地溫梯度控制著該地區(qū)的烴源巖有機質的熱演化.本文的研究成果可以為下一步油氣資源評價和勘探?jīng)Q策提供基礎依據(jù).

查干凹陷； 中、新生代； 構造-熱演化； 磷灰石裂變徑跡； 鏡質體反射率

1 引言

我國是石油消耗的大國之一，現(xiàn)在一半以上的原油(>2×108t)依賴于進口.傳統(tǒng)的大型油田都面臨著減產(chǎn)(大慶油田、勝利油田、中原油田老區(qū)等).油氣產(chǎn)量的減少嚴重制約著國民經(jīng)濟的高速發(fā)展，為了解決這一問題，在陸上，加大中生代斷陷盆地的資源調查，尋找油氣產(chǎn)量的接替區(qū)，包括內蒙古銀根-額濟納旗盆地(中生界有效沉積巖分布面積為10.4×104km2)、二連盆地(約10.0×104km2)、海拉爾盆地(約7.1×104km2).勘探已經(jīng)證實在中生代斷陷盆地具有一定量的油氣資源，海拉爾盆地探明儲量已經(jīng)超過10×108t，二連盆地和海拉爾盆地的年產(chǎn)量都已經(jīng)達到上百萬噸.內蒙古三大盆地中銀根-額濟納旗盆地(簡稱為銀-額盆地)勘探程度最低，2009年開始，中國石化集團公司中原油田分公司以中生代沉積地層最厚、認為最具勘探潛力的查干凹陷作為突破口，尋找白堊系的油氣資源，截至2013年底，已發(fā)現(xiàn)石油儲量5000余萬噸，建立了吉祥和如意兩個油田，證實查干凹陷具有較好的勘探前景.但是，關于查干凹陷的構造-熱演化的研究相對薄弱，已經(jīng)嚴重制約查干凹陷生烴史、排烴史、運聚史和成藏史等油氣方面的正確認識，直接影響著油氣勘探?jīng)Q策及勘探成功率.

沉積盆地構造-熱演化史的研究意義在于全面地評價盆地內烴源巖生、排、運、聚等對應的地質時間.進行構造-熱歷史的研究是正確認識油氣成藏及其演化的重要科學依據(jù)之一，也是重要的基礎課題之一.

近年來，隨著查干凹陷油氣勘探的不斷推進，鉆井數(shù)量的增加，為研究查干凹陷中、新生代的構造-熱演化提供保障.查干凹陷是內蒙古銀-額盆地最具勘探前景的凹陷，其構造-熱歷史的定量研究成果將指導銀-額盆地其他構造單元的油氣資源評價，對整個盆地的油氣勘探具有重要意義.本文將利用磷灰石裂變徑跡和鏡質體反射率古溫標耦合反演的方法恢復查干凹陷的構造-熱歷史，并探討構造-熱演化與沉積構造演化、巖漿活動的耦合關系及其對烴源巖熱演化的作用.

2 地質背景

查干凹陷位于銀-額盆地東部查干德勒蘇坳陷的中部，凹陷西臨西尼凸起，東與楚干凸起和白云凹陷相隔，西南為木巴圖隆起，東南緊靠狼山．構造活動相當復雜，其呈不規(guī)則的“菱形”展布，北東長60 km，北西寬40 km，勘探面積約2000 km2，是一個呈西北斷、東南超的單斷箕狀結構的中生代斷陷盆地，是銀-額盆地中最富油氣勘探潛力的凹陷．根據(jù)基底起伏、斷裂系統(tǒng)解釋成果及構造演化特征，查干凹陷劃分為“兩凹一凸”的構造格局，即西部次凹、東部次凹和毛敦次凸，各個構造單位又包括多個次一級構造單元(圖1a，b)．受區(qū)域構造運動的控制，查干凹陷僅發(fā)育白堊系和新生界地層，包括早白堊世巴音戈壁組(分兩段)、蘇紅圖組(分兩段)、銀根組、晚白堊世烏蘭海組及新生界(圖1c).沉積蓋層厚度西部次凹明顯比東部次凹的大，西部沉積中心在額很洼陷Y1-Y2井區(qū)，厚度達到5500 m；東部次凹沉積中心在罕塔廟洼陷北部地區(qū)，最大厚度為3500 m(圖2).由于燕山和喜山等造山運動，研究區(qū)在銀-額盆地整體構造演化控制的基礎之上，具體可以劃分為經(jīng)歷了3期構造階段：(1)早白堊世巴音戈壁組-銀根組沉積時期為裂陷階段，該時期斷裂活動強烈，伴隨多期火山活動，沉積一套中基性火山巖與碎屑巖的組合；(2)晚白堊世烏蘭蘇海組沉積時期為坳陷期，發(fā)育一套河流相為主的地層；(3)新生代為擠壓抬升期，凹陷發(fā)生局部的擠壓沖斷現(xiàn)象，發(fā)育一組逆沖斷層和反轉構造，局部地區(qū)接受新生界沉積.

圖1 查干凹陷構造分區(qū)圖(a)；構造剖面圖(b)；地層柱狀圖(c)Fig.1 (a) Structural unit division of the Chagan sag；(b) Structural profile map；(c) Stratigraphic column map

圖2 查干凹陷主要地層厚度圖(a—c)及沉積蓋層厚度圖(d) (a)蘇一段；(b)巴二段；(c)巴一段.Fig.2 Isopach map of the main strata (a—c) and sedimentary cover of the Chagan sag (d) (a)The 1st member of Suhongtu Formation；(b)The 2nd member of Bayingebi Formation; (c)The 1st member of Bayingebi Formation.

3 方法及基本參數(shù)

3.1 方法與原理

總的來說，沉積盆地構造-熱歷史的研究方法逐漸從單一古溫標到多溫標耦合反演的方向發(fā)展.

本文利用磷灰石裂變徑跡和鏡質體反射率古溫標進行約束，多溫標耦合反演相結合的方法恢復研究區(qū)的構造-熱歷史：首先利用古溫標數(shù)據(jù)反演得到樣品經(jīng)歷的古地溫梯度(或樣品經(jīng)歷的古地溫演化曲線)、剝蝕量及對應的地質時間；再以古地溫梯度及對應的地質時間(或樣品經(jīng)歷的古地溫演化曲線)為約束，結合盆地的構造演化史擬合出樣品經(jīng)歷的地溫梯度演化曲線，并以反演得到的剝蝕量為基礎，結合單井現(xiàn)今地層厚度恢復地層埋藏史；以地溫梯度演化曲線和埋藏史為基礎，利用盆地模擬軟件，采用正演的方法對單井的溫度史及生烴史進行模擬，模擬得到的溫標與實測的古溫標進行對比，如果模擬值與實測值具有很好的擬合度，則認為假設的地溫梯度演化曲線是可行的，如果擬合度較差，則修改地溫梯度演化曲線，直到模擬值與實測值具有很好的擬合度為止，此時的地溫梯度演化曲線則為該井經(jīng)歷的熱史.

3.2 古溫標及基本參數(shù)

在熱歷史模擬計算中需要的參數(shù)包括古溫標數(shù)據(jù)和基礎地質數(shù)據(jù).

(1) 古溫標數(shù)據(jù)

該凹陷古溫標數(shù)據(jù)包括14口井35個磷灰石裂變徑跡數(shù)據(jù)(表1)和15口井119個鏡質體反射率數(shù)據(jù)(表2，圖3).其中，鏡質體反射率數(shù)據(jù)隨深度變化具有較好的線性關系(圖3)，暗示受同一地溫梯度控制；14口井35個磷灰石裂變徑跡數(shù)據(jù)是在中國科學院高能物理研究所核物理實驗室測試完成的，這些數(shù)據(jù)分布在銀根組、蘇二段、蘇一段、巴二段和巴一段.35個磷灰石裂變徑跡年齡分布在5.5±0.7～229±16 Ma，磷灰石裂變徑跡長度分布在10.6±1.7～13.2±1.8 μm，其中14個磷灰石裂變徑跡年齡比地層年齡大(圖4)，徑跡年齡代表物源區(qū)的構造-熱事件；其他磷灰石裂變徑跡年齡比地層年齡小(圖4)，徑跡年齡代表查干凹陷的構造-熱事件，可以用來恢復查干凹陷的構造-熱歷史.從徑跡年齡與深度的關系可以得出，在3620 m徑跡年齡為0 Ma，即該凹陷磷灰石裂變徑跡完全退火的深度為3620 m.根據(jù)查干凹陷平均地溫梯度33.6 ℃/km(左銀輝等，2013b)，計算得到磷灰石裂變徑跡完全退火開始的溫度為130 ℃，這與前人研究成果(Gleadow et al．，1983)相當.

(2) 基礎地質參數(shù)

模擬計算中的參數(shù)包括巖性參數(shù)、現(xiàn)今地表溫度數(shù)據(jù)、現(xiàn)今地溫梯度、大地熱流、巖石熱物理參數(shù)、地層分層、地層年齡及主要地質時期的剝蝕量等數(shù)據(jù).查干凹陷地表溫度數(shù)據(jù)、地溫梯度及大地熱流分布情況見左銀輝等(2013b)的研究成果.巖性參數(shù)主要包括巖石的孔隙度、滲透率、各巖層的砂泥巖含量、砂泥巖的壓實曲線等，這些參數(shù)采用查干凹陷實測值.此外，還包括壓實系數(shù)和初始孔隙度等數(shù)據(jù)，則依據(jù)各凹陷的實際數(shù)據(jù)利用Sclater 和Christie(1985)的方法進行回歸得到.古地表溫度取查干凹陷年平均溫度(9 ℃)，并設在地質歷史時期不變.磷灰石裂變徑跡退火模型采用Laslett等(1987)的模型，鏡質體反射率模擬熱歷史采用Easy%Ro模型(Sweeney and Burnham，1990).

圖3 查干凹陷鏡質體反射率與深度的關系圖 (鏡質體反射率數(shù)據(jù)見表2)Fig.3 Vitrinite reflectance (Ro) data versus depth in the Chagan sag (Ro data shown in the Table 2)

序號井名深度(m)地層nρs(105cm2)(Ns)ρi(105cm2)(Ni)ρd(105cm2)(Nd)P(χ2)(%)Age(Ma)(±1σ)L(μm)(N)備注1Y11342.0K1y284.281(936)9.107(1991)13.106(7312)99127±912.2±1.7(104)>2Y11497.0K1y70.301(10)10.504(349)13.155(7312)94.17.8±2.6-<3Y22718.2K1s1114.465(219)9.826(482)13.301(7312)0.1107±1813.2±1.8(19)<4Y31817.7K1s2282.380(533)8.003(1792)13.008(7312)1.580±612.1±1.8(54)<5Y32978.7K1s1340.560(124)12.254(2714)13.252(7312)64.713±1-<6Y33343.8K1b2320.414(23)17.195(955)13.179(7312)85.46.6±1.4-<7Y41057.2K1s2284.153(999)13.707(3297)13.350(7312)42.884±612.1±1.6(104)<8Y41529.5K1s1282.895(988)11.787(4023)13.301(7312)8.768±511.7±1.6(104)<9Y41651.6K1b2275.151(397)14.027(1081)13.301(7312)5.2101±910.6±1.7(37)<10Y41703.8K1b2284.595(411)12.308(1101)13.350(7312)0.298±1110.8±1.7(76)<11Y41832.1K1b1282.740(415)9.110(1380)13.350(7312)2.581±811.2±1.7(31)<12Y53418.0K1b2320.791(63)16.666(1327)12.984(7312)1.113±210.8±1.7(14)<13Y61591.0K1s1225.703(874)13.657(2093)12.886(7312)33.4111±811.7±1.6(101)>14Y61761.5K1s1283.899(813)11.353(2367)13.350(7312)0.995±811.9±1.8(100)<15Y62179.0K1b291.592(77)8.602(416)12.642(7312)049±712.8±2.1(7)<16M31257.6K1s2286.754(1517)17.301(3886)12.764(7312)0107±812.0±1.8(101)<17M31797.8K1s1284.429(974)9.577(2106)12.812(7312)4.6121±911.4±2.0(101)>18M31832.1K1s1283.257(860)9.98(2635)12.861(7312)0.0788±711.5±2.0(121)<19M61147.9K1s2285.752(2044)9.610(3415)13.350(7312)0158±1412.5±1.5(102)>20M61082.6K1s2287.388(1048)9.546(1354)13.301(7312)6.1211±1512.7±1.4(111)>21M8718.3K1y297.299(1727)16.715(3955)13.057(7312)7.8118±812.7±1.5(109)>22M8811.7K1y286.339(1035)14.399(2351)13.106(7312)7.9119±912.3±1.8(105)>23M8916.8K1s2286.655(859)11.072(1429)13.155(7312)40.4162±1212.5±1.4(98)>24M91300.0K1s2124.603(471)11.864(1214)13.131(7312)22.9105±911.6±1.7(73)<25M91402.0K1s2239.131(1420)17.555(2730)13.033(7312)0136±1111.1±2.0(105)>26M9765.2K1y104.880(213)8.959(391)13.277(7312)0.2133±2211.8±1.9(23)>27M112077.4K1s1284.189(616)10.248(1507)13.350(7312)0.5114±1011.9±1.6(80)>28M113296.9K1b2300.322(98)15.251(4646)12.593(7312)43.05.5±0.711.5±1.8(26)<29M113540.7K1b1320.373(86)13.255(3058)12.691(7312)69.77.4±0.911.5±0.9(4)<30L12443.0K1s1172.225(168)10.133(765)13.350(7312)82.161±611.1±2.0(17)<31L13357.2K1b2170.962(188)18.239(3566)13.350(7312)14.715±110.7±1.3(9)<32L11977.5K1s1283.450(1040)10.534(3176)13.350(7312)090±811.1±1.8(113)<33J6752.0K1s2124.449(234)8.519(448)13.106(7312)60141±1412.4±1.4(27)>34J61150.0P288.751(1692)10.209(1974)13.057(7312)97.1229±1611.5±1.3(103)>35CD1679.2K1s2305.184(1255)10.830(2622)13.008(7312)0122±1212.3±1.8(103)>

注：n為測量的磷灰石顆粒數(shù)；ρi為外部探測器中的誘發(fā)徑跡密度；ρs為自發(fā)徑跡密度；ρd為自發(fā)徑跡密度；Ni、Ns和Nd為測量的徑跡數(shù)；P(χ2)為χ2概率；Age±1σ為樣品值裂變徑跡年齡；L為平均徑跡長度；N為測量的徑跡數(shù)；“<”為磷灰石裂變徑跡年齡小于地層年齡；“>”為磷灰石裂變徑跡年齡大于地層年齡；樣品測試在中國科學院高能物理研究所核物理實驗室完成.

表2 查干凹陷鏡質體反射率Table 2 Vitrinite reflectance data of the Chagan sag

巖石熱物理參數(shù)主要包括巖石熱導率、巖石生熱率等參數(shù)(左銀輝等，2013b)，地層分層采用鉆孔實際測量值，各地層底界年齡數(shù)據(jù)具體為：新生界(Cz)65 Ma，上白堊統(tǒng)烏蘭蘇海組(K2w)95 Ma，下白堊統(tǒng)銀根組(K1y)100 Ma，下白堊統(tǒng)蘇紅圖組二段(K1s2)105 Ma，下白堊統(tǒng)蘇紅圖組一段(K1s1)110 Ma，下白堊統(tǒng)巴音戈壁組二段(K1b2)128 Ma，下白堊統(tǒng)巴音戈壁組一段(K1b1)135 Ma(圖1c).

查干凹陷存在3個區(qū)域不整合，具體為蘇二段與銀根組、銀根組與烏蘭蘇海組和烏蘭蘇海組與新生界之間，利用泥巖聲波時差測井計算得到典型井蘇二段與銀根組之間的剝蝕量，利用鏡質體反射率反演獲得典型井銀根組與烏蘭蘇海組之間的剝蝕量，以及利用磷灰石裂變徑跡反演得到烏蘭蘇海組與新生界之間的剝蝕量(表3).

表3 查干凹陷典型井主要地質時期的剝蝕量Table 3 Erosion amounts in the main geological periods of the typical wells in the Chagan sag

圖4 查干凹陷磷灰石裂變徑跡與深度的關系(F為頻率)Fig.4 Apatite fission-track age and length versus depth in the Chagan sag

圖5 Y4井(1529.5 m，K1s1)溫度-時間演化圖 使用Monte Carlo模型模擬了10000條溫度-時間路徑；可以接受的溫度-時間路徑有1668條(左圖酒綠色線條)；好的溫度-時間路徑為189條(左圖粉紅色線條)；粗的黑線條為最佳溫度-時間路徑. (a)裂變徑跡年齡分布圖； (b)裂變徑跡長度分布圖，圖中黑線條為模擬的裂變徑跡.GOF為擬合度.Fig.5 Thermal history of sample Y4 (1529.5 m， K1s1) in well Y4 Ten thousand thermal paths were tried using the Monte Carlo inverse modeling method；1668 paths are acceptable (wine green lines) and 189 paths are good (pink lines)．The thick black line is the best temperature path．(a) Apatite fission-track ages．(b) Apatite fission-track length distribution．The thick black line is modeled Apatite fission-track length．GOF=goodness of fit．

圖6 樣品Y3(1817.7m，K1s2)溫度-時間演化圖 使用Monte Carlo模型模擬了10000條溫度-時間路徑；可以接受的溫度-時間路徑有3389條(左圖酒綠色線條)；好的溫度-時間路徑為2181條(左圖粉紅色線條)；粗的黑線條為最佳溫度-時間路徑.(a)裂變徑跡年齡分布圖；(b)裂變徑跡長度分布圖，圖中 黑線條為模擬的裂變徑跡.GOF為擬合度.Fig.6 Thermal history of sample Y3 (1817.7 m， K1s2) in well Y3 Ten thousand thermal paths were tried using the Monte Carlo inverse modeling method；3389 paths are acceptable (wine green lines) and 2181 paths are good (pink lines)．The thick black line is the best temperature path．(a) Apatite fission-track ages．(b) Apatite fission-track length distribution．The thick black line is modeled apatite fission-track length．GOF=goodness of fit．

圖7 Y4井埋藏史及熱史模擬結果Fig.7 The burial and thermal histories of well Y4

圖8 Y3井埋藏史及熱史模擬結果Fig.8 The burial and thermal histories of well Y3

圖9 查干凹陷典型井熱演化歷史Fig.9 Thermal gradient evolution history of the studied wells in the Chagan sag

4 中、新生代構造-熱演化

4.1 反演-磷灰石裂變徑跡

這次對2個磷灰石裂變徑跡樣品進行了熱史反演(圖5，圖6).在模擬中，以關鍵地質時期的溫度為約束條件，利用Monte Carlo模型進行了10000次熱史反演，獲得一定數(shù)量的最佳溫度-時間路徑，這些路徑代表樣品經(jīng)歷的最有可能的溫度-時間演化路徑.2口井的熱史反演結果均顯示查干凹陷經(jīng)歷了銀根組末期和烏蘭蘇海組沉積中期的構造抬升，并且在銀根組沉積時期溫度達到最高.

4.2 正演-鏡質體反射率

以磷灰石裂變徑跡反演的結果為約束條件，利用鏡質體反射率對Y4井和Y3井進行熱史正演計算.模擬的鏡質體反射率數(shù)據(jù)與實測值具有很好的線性關系(圖7，圖8)，由此可見，模擬結果是可信的.模擬結果顯示Y4井先后在蘇紅圖組沉積末期、銀根組沉積末期及烏蘭蘇海組沉積中期經(jīng)歷了3次較明顯的抬升剝蝕，其中銀根組沉積末期的抬升剝蝕最大，達到1120 m.在烏蘭蘇海組沉積之前為裂陷發(fā)育階段，表現(xiàn)為快速沉降，烏蘭蘇海組沉積時期為坳陷發(fā)育階段，沉積速率則由早到晚逐漸減小.從地溫演化來看，在銀根組沉積末期古地溫達到最大，超過130 ℃.從熱史模擬結果看，Y4井在巴音戈壁組沉積時期地溫梯度在45～49 ℃/km之間，自蘇紅圖組沉積時期開始，地溫梯度逐漸增高，到銀根組沉積末期，地溫梯度達到最大，為56 ℃/km，自烏蘭蘇海組沉積至今，查干凹陷表現(xiàn)為熱沉降階段，地溫梯度呈逐漸下降的趨勢，現(xiàn)今僅為33 ℃/km，但是由于烏蘭蘇海組沉積早期快速沉積和新沉積物小的巖石熱導率造成地溫梯度稍有升高(圖7).Y3井也在蘇紅圖組沉積末期、銀根組沉積末期及烏蘭蘇海組沉積中期經(jīng)歷了3次較明顯的抬升剝蝕，其中銀根組沉積末期的抬升幅度最大.從熱史模擬結果看，Y3井的地溫梯度演化與Y4井相似，都在銀根組地溫梯度達到最大，為55 ℃/km(圖8).

此外，利用以上的方法還模擬了11口井的熱史，模擬結果顯示查干凹陷經(jīng)歷了以下4個熱演化階段(圖9)：(1)巴一段—蘇紅組沉積時期(K1b-K1s)：地溫梯度快速增加階段，地溫梯度由巴音戈壁組沉積開始的42～47 ℃/km逐漸增加至蘇紅圖組沉積末期的46～52 ℃/km；(2)銀根組沉積時期(K1y)：地溫梯度高峰階段，此時地溫梯度達到50～58 ℃/km，具有裂陷構造區(qū)的熱流狀態(tài)；(3)烏蘭蘇海組沉積時期(K2w)：高地溫延續(xù)階段，地溫梯度為39～48 ℃/km，這是由于烏蘭蘇海組較厚的新沉積物具有低的巖石熱導率，使得在烏蘭蘇海組的地溫梯度略有升高，在中晚期構造抬升，地溫梯度又開始降低；(4)新生代(Cz)：熱沉降階段，此階段主要受喜山構造運動的影響，查干凹陷主要處于抬升剝蝕期，新生代沉積較薄，地殼處于均衡調整期，地溫梯度逐漸降低，現(xiàn)今為31～34 ℃/km，具有偏高的古地溫梯度.

5 討論

5.1 構造-熱歷史與沉積構造演化、巖漿活動的耦合關系

通過以上研究表明查干凹陷在白堊紀具有較高的古地溫梯度，這種地熱狀態(tài)與其經(jīng)歷的構造運動及巖漿活動密切相關.

查干凹陷是在古生代褶皺基底上發(fā)展起來的中、新生代斷陷盆地，在中晚侏羅世(160～140 Ma)，蒙古鄂霍茨克洋閉合，西伯利亞板塊與華北板塊開始發(fā)生陸陸碰撞，查干凹陷及鄰區(qū)在兩板塊碰撞早期處于擠壓狀態(tài)，開始抬升，巖石圈加厚，因此，查干凹陷缺失三疊系和侏羅系.早白堊世開始，較重的大洋板片對大陸巖石圈的牽引力使地面巖石圈的應力由強烈的擠壓應力轉變?yōu)槔瓘垜?，中晚侏羅世增厚的巖石圈因重力失穩(wěn)發(fā)生伸展垮塌，查干凹陷開始裂陷，巖石圈減薄，地幔物質開始上涌，地溫梯度增加.此階段，盆地開始接受早白堊世巴音戈壁組的沉積.早白堊世蘇紅圖組沉積時期，區(qū)域拉張應力場進一步強化，伸展斷裂作用加強，并且俯沖的蒙古-鄂霍茨克大洋板片發(fā)生了斷離(底侵作用)，在這共同作用下引發(fā)了多期火山活動，即查干凹陷在蘇紅圖組沉積時期共發(fā)生10期火山活動，蘇紅圖組火山巖厚度接近800 m，此階段，巖石圈的大幅度拉伸減薄，地幔物質上涌，帶來大量地幔熱量，導致早白堊世具有高的地熱背景，并在早白堊世晚期地溫梯度達到最大，這與渤海灣盆地早白堊世高的地熱背景及巖石圈減薄具有一致性(Zuo et al．，2013；Qiu et al．，2014).這種底侵作用使地殼物質發(fā)生了部分熔融，從而導致在火山作用的早期還有大規(guī)模的酸性巖漿作用，巖石圈的伸展跨塌造成了平行于蒙古-鄂霍茨克縫合帶的半地塹盆地群(銀額盆地、二連盆地、東戈壁盆地、海拉爾盆地)初步形成；而且，在145～135 Ma，古太平洋板塊開始沿歐亞東緣向北西俯沖(周新華等，2001)，這種板塊邊界作用引起的板內效應也可能是引起查干凹陷及鄰區(qū)巖石圈拉張的重要因素.由此可見，查干凹陷早白堊世具有板內被動裂谷的特征，形成于巖石圈伸展減薄條件下的被動裂谷環(huán)境，也與查干凹陷在早白堊世具有高的地熱背景相一致.晚白堊世開始，盆地進入坳陷期，巖石圈開始均衡調整，地溫梯度開始降低，但是晚白堊世烏蘭蘇海組較厚的新沉積物具有低的巖石熱導率，使得在烏蘭蘇海組的地溫梯度略有升高.

圖10 查干凹陷典型井巴二段烴源巖底成熟度演化歷史Fig.10 Maturation history of the bottom of the Bayingebi 2 Formation in the Chagan sag

5.2 構造-熱演化與烴源巖熱演化的關系

白堊紀高的地熱狀態(tài)有利于油氣的生成，在查干凹陷的平均古生烴門限在1703 m(表4)；這明顯比冷的地熱狀態(tài)對應的門限深度偏淺，例如渤海海域(新生代盆地)受現(xiàn)今地溫場控制(地溫梯度為31.8 ℃/km)，其生烴門限在2500 m(Zuo et al．，2011).同時，研究中發(fā)現(xiàn)早白堊世銀根組沉積時期的地溫達到最大，此時的古地溫場控制著查干凹陷的烴源巖熱演化歷史，也是查干凹陷生烴結束的時間(圖10).早白堊世銀根組的古地溫場將進一步控制查干凹陷的主生烴時期、主排烴時期及成藏關鍵時期等，再結合儲、蓋、圈閉形成時間等的研究成果，可以指出有利油氣聚集區(qū)帶．因此，查干凹陷的構造-熱歷史的研究成果對該地區(qū)的油氣勘探?jīng)Q策具有重要指導意義．

表4 查干凹陷蘇一段底古生烴門限Table 4 Paleo-hydrocarbon generation thresholds for the bottom of the Suhongtu 1 Formation in the Chagan sag

6 結論

本文利用鏡質體反射率和磷灰石裂變徑跡兩種古溫標耦合反演的方法恢復了查干凹陷的構造-熱歷史，揭示了查干凹陷在白堊紀具有高的地熱狀態(tài)，并且構造-熱演化劃分為以下四個階段：地溫梯度快速增加階段(K1b-K1s)，地溫梯度由巴音戈壁組沉積開始的42～47 ℃/km逐漸增加至蘇紅圖組沉積末期的46～52 ℃/km；地溫梯度高峰階段(K1y)，此時地溫梯度達到50～58 ℃/km，具有裂陷構造區(qū)的熱流狀態(tài)；高地溫延續(xù)階段(K2w)，地溫梯度為39～48 ℃/km；熱沉降階段(Cz)，此階段主要受喜山構造運動的影響，查干凹陷處于抬升剝蝕期，新生代沉積較薄，地殼處于均衡調整期，地溫梯度逐漸降低，現(xiàn)今為31～34 ℃/km，地溫梯度較區(qū)域背景略微偏高.在熱史恢復的基礎上，通過正演獲得了主力烴源巖成熟度的演化過程，發(fā)現(xiàn)查干凹陷的烴源巖有機質熱演化主要受早白堊銀根組沉積時期的古地溫梯度控制，此時烴源巖熱演化程度達到最高，隨后生烴停止.

致謝 感謝中國地質大學(北京)袁萬明教授在磷灰石裂變徑跡測試及數(shù)據(jù)分析工作中的指導.

Crowhurst P V, Green P F, Kamp P J J. 2002. Appraisal of (U-Th)/He apatite thermochronology as a thermal history tool for hydrocarbon exploration: An example from the Taranaki Basin, New Zealand.AAPGBulletin, 86(10): 1801-1819.

Farley K A. 2000. Helium diffusion from apatite: General behavior as illustrated by Durango fluorapatite.JournalofGeophysicalResearch, 105(B2): 2909-2914.

Fernandes P, Rodrigues B, Borges M, et al. 2013. Organic maturation of the Algarve Basin (southern Portugal) and its bearing on thermal history and hydrocarbon exploration.MarineandPetroleumGeology, 46: 210-233.

Gleadow A J W, Duddy I R, Lovering J F. 1983. Fission track analysis: A new tool for the evaluation of thermal historyies and hydrocarbon potential. Aust. Petrol. Explor. Assoc. J., 23: 93-102.

Laslett G M, Green P F, Duddy I R, et al. 1987. Thermal annealing of fission tracks in apatite. 2. A quantitative analysis.ChemicalGeology:IsotopeGeoscienceSection, 65(1): 1-13.

He L J, Wang J Y. 2004. Tectono-thermal modelling of sedimentary basins with episodic extension and inversion, a case history of the Jiyang Basin, North China.BasinResearch, 16(4): 587-599.

Hu S, Fu M, Yang S, Barry P K, et al. 2007. Palaeogeothermal response and record of Late Mesozoic lithospheric thinning in the eastern North China Craton. ∥Zhai M G, Windley B F, Kusky T M, eds. Mesozoic Sub-Continental Lithospheric Thinning Under Eastern Asia. London (Spec).GeologicalSociety, 280: 267-280.

Ketcham R A, Carter A, Donelick R A, et al. 2007. Improved modeling of fission-track annealing in apatite.AmericanMineralogist, 92(5-6): 799-810. McCormack N, Clayton G, Fernandes P. 2007. The thermal history of the Upper Palaeozoic rocks of southern Portugal.MarineandPetroleumGeology, 24(3): 145-150.

McKenzie D. 1978. Some remarks on the development of sedimentary basins.EarthandPlanetaryScienceLetters, 40(1): 25-32.

Menzies M, Xu Y G, Zhang H F, et al. 2007. Integration of geology, geophysics and geochemistry: a key to understanding the North China Craton.Lithos, 96(1-2): 1-21.

Qiu N S, Qin J Z, McInnes BI A, et al. 2008. Tectonothermal evolution of the Northeastern Sichuan Basin: Constraints from apatite and zircon (U-Th)/He ages and vitrinite reflectance data.GeologicalJournalofChinaUniversities(in Chinese), 14(2): 223-230.

Qiu N S, Peter R, Mei Q H, et al. 2009. Application of the (U-Th)/He thermochronometry to the tectono-thermal evolution of sedimentary basin—A case history of Well KQ1 in the Tarim Basin.ChineseJournalofGeophysics(in Chinese), 52(7): 1825-1835, doi: 10.3969/j.issn.0001-5733.2009.07.017.

Qiu N S, Wang J Y, Mei Q H, et al. 2010. Constraints of (U-Th)/He ages on early Paleozoic tectonothermal evolution of the Tarim Basin, China.ScienceChina:EarthSciences, 53(7): 964-976.

Qiu N S, Chang J, Zuo Y H, et al. 2012a. Thermal evolution and maturation of lower Paleozoic source rocks in the Tarim Basin, northwest China.AAPGBulletin, 96(5): 789-821.

Qiu N S, Chang J, Li J W, et al. 2012b. New evidence on the Neogene uplift of South Tianshan constraints from the (U-Th)/He and AFT ages of borehole samples of the Tarim basin and implications for hydrocarbon generation.InternationalJournalofEarthSciences, 101(6): 1625-1643.

Qiu N S, Zuo Y H, Chang J, et al. 2014. Geothermal evidence of Mesozoic-Cenozoic lithosphere thinning in the Jiyang sub-basin, Bohai Bay Basin, eastern North China Craton.GondwanaResearch, 26(3-4): 1079-1092.

Reiners P W, Farley K A, Hickes H J. 2002. He diffusion and (U-Th)/He thermochronometry of zircon: Initial results from Fish Canyon Tuff and Gold Butte.Tectonophysics, 349(1-4): 297-308.

Ren Z L, Liu C Y, Zhang X H, et al. 2000. Recovery and comparative research of thermal history on Jiuqu.ChineseJournalofGeophysics(in Chinese), 43(5): 635-645.

Sahu H S, Raab M J, Kohn B P, et al. 2013. Thermal history of the Krishna-Godavari basin, India: Constraints from apatite fission track thermochronology and organic maturity data.JournalofAsianEarthSciences, 73: 1-20.

Sclater J G, Christie P A F. 1985. Continental stretching: an explanation of the post-Mid-Cretaceous subsidence of the central North Sea Basin.JournalofGeophysicalResearch, 85: 3711-3739.

Shuster D L, Farley K A. 2009. The influence of artificial radiation damage and thermal annealing on helium diffusion kinetics in apatite.GeochimicaetCosmochimicaActa, 73(1): 183-196.

Sleep N H, Snell N S. 1976. Thermal contraction and flexure of mid-continent and Atlantic marginal basins.GeophysicalJournalInternational, 45(1): 125-154. Spiegel C, Kohn B, Belton D, et al. 2009. Apatite (U-Th-Sm)/He thermochronology of rapidly cooled samples: the effect of He implantation.EarthandPlanetaryScienceLetters, 285(1-2): 105-114.

Sweeney J J, Burnham A K. 1990. Evaluation of a simple model of vitrinite reflectance based on chemical kinetics.AAPGBulletin, 74(10): 1559-1570. Wolf R A, Farley K A, Kass D M. 1998. Modeling of the temperature sensitivity of the apatite (U-Th)/He thermochronometer.ChemicalGeology, 148(1-2): 105-114. Zhang J, Song H B, Li J B. 2006. Thermal modeling of the tectonic evolution of the southwest subbasin in the South China Sea.ChineseJournalofGeophysics(in Chinese), 48(6): 1357-1365.

Zhou X H, Yang G H, Yang J H, et al. 2001. Sr-Nd-Pb isotope mapping of the Late Mesozoic volcanic rocks across northern margin of North China Craton and implication to geodynamic processes.Geochimica(in Chinese), 30(1): 10-23.

Zuo Y H, Qiu N S, Pang X Q, et al. 2010. Hydrocarbon kitchen evolution of E2s3source rock of the Bohai Offshore area, North China.ChineseJournalofGeophysics(in Chinese), 53(10): 2415-2426, doi: 10.3969/j.issn.0001-5733.2010.10.015.

Zuo Y H, Qiu N S, Zhang Y, et al. 2011. Geothermal regime and hydrocarbon kitchen evolution of the offshore Bohai Bay basin, North China.AAPGBulletin, 95(5): 749-769.

Zuo Y H, Ma W M, Deng Y X, et al. 2013a. Mesozoic and Cenozoic thermal history and source rock thermal evolution history in the Chagan sag, Inner Mongolia.EarthScience-JournalofChinaUniversityofGeosciences(in Chinese), 38(3): 553-560.

Zuo Y H, Qiu N S, Deng Y X, et al. 2013b. Terrestrial heat flow in the Qagan, Inner Mongolia.ChineseJournalofGeophysics(in Chinese), 56(9): 3038-3050, doi: 10.6038/cjg20130917.

Zuo Y H, Qiu N S, Pang X Q, et al. 2013. Geothermal evidence of the Mesozoic and Cenozoic lithospheric thinning in the Liaohe depression.JournalofEarthScience, 24(4): 529-540.

Zuo Y H, Qiu N S, Li J W, et al. 2015. Meso-Cenozoic tectono-thermal evolution history in Bohai Bay Basin, North China.JournalofEarthScience, 26(3):352-360, doi: 10.1007/s12583-014-0500-0.

附中文參考文獻

邱楠生, 秦建中, McInnes B I A等. 2008. 川東北地區(qū)構造-熱演化探討—來自(U-Th)/He年齡和Ro的約束. 高校地質學報, 14(2): 223-230.

邱楠生, Peter R, 梅慶華等. 2009. (U-Th)/He年齡在沉積盆地構

造-熱演化研究中的應用—以塔里木盆地KQ1井為例. 地球物理學報, 52(7): 1825-1835, doi: 10.3969/j.issn.0001-5733.2009.07.017.

任戰(zhàn)利, 劉池陽, 張小會等. 2000. 酒泉盆地群熱演化史恢復及其對比研究. 地球物理學報, 43(5): 635-645.

張健, 宋海斌, 李家彪. 2006. 南海西南海盆構造演化的熱模擬研究. 地球物理學報, 48(6): 1357-1365.

周新華, 楊國輝, 楊進輝等. 2001. 華北克拉通北緣晚中生代火山巖Sr-Nd-Pb同位素填圖及其構造意義. 地球化學, 30(1): 10-23.

左銀輝, 邱楠生, 龐雄奇等. 2011. 渤海海域沙三段烴源灶演化特征研究. 地球物理學報, 53(10): 2415-2426, doi: 10.3969/j.issn.0001-5733.2010.10.015.

左銀輝, 馬維民, 鄧已尋等. 2013a. 查干凹陷中、新生代熱史及烴源巖熱演化. 地球科學-中國地質大學學報, 38(3): 553-560.

左銀輝, 邱楠生, 鄧已尋等. 2013b. 查干凹陷大地熱流. 地球物理學報, 56(9): 3038-3050, doi: 10.6038/cjg20130917.

(本文編輯 胡素芳)

Mesozoic and Cenozoic tectono-thermal evolution history in the Chagan Sag，Inner Mongolia

ZUO Yin-Hui1，2， ZHANG Wang1，2， LI Zhao-Ying3， LI Jia-Wei4， HAO Qing-Qing5， HU Jie2

1StateKeyLaboratoryofOilandGasGeologyandExploitation，ChengduUniversityofTechnology，Chengdu610059，China2CollegeofEnergyResources，ChengduUniversityofTechnology，Chengdu610059，China3GeologicalExploration&DevelopmentResearchInstitute，CNPC，ChuanqingDrillingEngineeringCompanyLimited，Chengdu610051，China4DivisionofLithosphereEvolution，InstituteofGeologyandGeophysics，ChineseAcademyofSciences，Beijing100029，China5InstituteofMineralResourcesResearch，ChinaMetallurgicalGeologyBureau，Beijing100025，China

The Chagan sag has the greatest oil and gas exploration potential among the sags in the Yingen-Ejinaqi Basin, Inner Mongolia. To reveal the tectono-thermal evolution history of the Chagan sag, this article utilized a combination of forward and inversion methods on the basis of 7 apatite fission track and 119 vitrinite reflectance data to reconstruct the Mesozoic and Cenozoic tectono-thermal evolution history in the Chagan sag. The results show that the Chagan sag had a high geothermal gradient during the Cretaceous, and it experienced the following 4 stages of tectono-thermal evolutions: (1) a rapid geothermal gradient increase stage from the Bayingebi Formation depositional period to the Suhongtu Formation depositional period, during which the geothermal gradient increased to 46～52 ℃/km at the end of the Suhongtu Formation depositional period; (2) a geothermal gradient peak stage during the Yingen Formation depositional period, with maximum geothermal gradient ranged from 50 to 58 ℃/km; (3) a high geothermal gradient continuation stage during the Wulansuhai Formation depositional period, with maximum geothermal gradient ranged from 39 to 48 ℃/km; and (4) a thermal subsidence stage during the Cenozoic, during which the Chagan sag is in the uplift and erosion stage due to the Himalayan movement and the geothermal gradient gradually decreased to 31～34 ℃/km at the present day. Moreover, the tectono-thermal evolution was matched with the tectonic evolution and volcanic activities in the Chagan sag. During the Early Cretaceous, the intraplate rift was developed, lithosphere was thinned and multi-phase intense volcanoes were erupted in the Chagan sag, so that a large amount of energy was released to the surface from the deep crust, resulting in a high geothermal gradient during this period. In addition, the high geothermal gradient during the Cretaceous was favorable for hydrocarbon generation, corresponding to a shallow paleo-generation threshold, and the Early Cretaceous geothermal fields controlled the hydrocarbon generation of the Chagan sag. This work may provide new insights for the understanding of the oil and gas exploration potential of the Chagan sag.

Chagan sag; Mesozoic and Cenozoic; Tectono-thermal evolution; Apatite fission track; Vitrinite reflectance

10.6038/cjg20150714.Zuo Y H, Zhang W, Li Z Y，et al. 2015. Mesozoic and Cenozoic tectono-thermal evolution history in the Chagan Sag，Inner Mongolia.ChineseJ.Geophys. (in Chinese),58(7):2366-2379,doi:10.6038/cjg20150714.

國家自然科學基金項目(41402219，41125010，91114202)資助．

左銀輝，男，1980年生，副教授，博士，主要從事沉積盆地現(xiàn)今地溫場，沉積盆地構造-熱演化，油氣資源評價和地球動力學方面的研究．E-mail:zuoyinhui@tom.com

10.6038/cjg20150714

P314

2014-09-17，2015-05-08收修定稿

左銀輝， 張旺， 李兆影等. 2015. 查干凹陷中、新生代構造-熱演化史.地球物理學報，58(7):2366-2379,