地球物理学报  2014, Vol. 57 Issue (5): 1554-1567   PDF    
利用自然伽马测井计算准噶尔盆地沉积层生热率及其热流贡献
饶松1,2, 朱传庆3, 廖宗宝4, 姜光政1,2, 胡圣标1, 汪集旸1    
1. 中国科学院地质与地球物理研究所 岩石圈演化国家重点实验室, 北京 100029;
2. 中国科学院大学, 北京 100049;
3. 中国石油大学(北京), 北京 102249;
4. 中海油研究总院, 北京 100027
摘要:沉积层放射性生热的热流贡献(沉积层热流)是沉积盆地大地热流的重要组成部分,能够有效促进中国西部“冷”盆深层-超深层烃源岩的增温和热演化.本文利用不同的自然伽马(GR)-生热率(A)经验关系式分别计算了准噶尔盆地不同构造单元16口钻孔共6120个沉积层生热率,通过与实测生热率的统计对比,确定了适用于研究区的GR-A经验关系,建立了准噶尔盆地地层生热率柱,据此计算了研究区沉积层热流贡献,并以盆参2井为例定量分析了沉积层热流的增温效应.结果表明,准噶尔盆地沉积层平均生热率为1.179±0.339 μW·m-3,总体上随着时代变老,沉积层生热率呈现出递减趋势.准噶尔盆地沉积层热流平均为7.9±4.9 mW·m-2,约占地壳热流的29.2%和大地热流的19.6%,区域上与盆地沉积层厚度大体一致,表现为中央坳陷最高,北天山山前冲断带变化较大,陆梁隆起和西部隆起次之,东部隆起和乌伦古坳陷最低.沉积层热流能够有效增高深层—超深层烃源层受热温度,促进有机质热演化,如在考虑和忽略沉积层生热的两种情况下计算的盆参2井下侏罗统三工河组烃源岩底部(5300 m)温度差异最大为7.3 ℃,这显然对于地温梯度小、主体油气藏埋深大的准噶尔盆地油气资源评价和勘探目标优选具有重要意义.
关键词生热率     沉积层热流     自然伽马测井     准噶尔盆地     深层—超深层    
Heat production rate and heat flow contribution of the sedimentary formation in Junggar Basin, northwest China
RAO Song1,2, ZHU Chuan-Qing3, LIAO Zong-Bao4, JIANG Guang-Zheng1,2, HU Sheng-Biao1, WANG Ji-Yang1    
1. State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China;
2. University of Chinese Academy of Sciences, Beijing 100049, China;
3. China University of Petroleum (Beijing), Beijing 102249, China;
4. CNOOC Research Institute, Beijing 100027, China
Abstract: As an important part of the terrestrial heat flow, radioactive heat flow production of sedimentary formation can warm the deep and ultra-deep source rocks and promote their maturity evolution effectively in the‘cold’basins in western China which are characterized by low geothermal gradient and low terrestrial heat flow. 6,120 heat production rates of sedimentary formation in 16 boreholes among six tectonic units in Junggar Basin were acquired using different empirical relationships between natural gamma ray (GR) and heat production rate (A). By comparison with the measured heat production rates, the applicable GR-A empirical relationship to the study area was determined. Then, the heat production rate columns of the sedimentary formation in Junggar Basin were established and the heat flow contribution of the sediments was calculated. Finally, as an example, its warming effect of Well Pc2 was calculated quantitatively. The results showed that the average heat production rate of the sediments in Junggar Basin was 1.179±0.339 μW·m-3 and it decreased gradually as the formation age became older. The Quaternary formation was with the maximum heat production rate while the Carboniferous formation minimum. The average heat flow contribution of the sediments was 7.9±4.9 mW·m-2, which was approximately 29.2% of average crustal heat flow and 19.6% of average terrestrial heat flow. The heat flow contribution of sediments was consistent with the thickness of sediments and showed the following characteristics. (1) Relatively high heat flow contribution of sediments was confined to the Central Depression. (2) Southern Depression varied dramatically on heat flow contribution of the sediments. (3) The Luliang Uplift and Western Uplift were characterized by lower heat flow contribution of the sediments. (4) The lowest heat flow of sediments occurred in the Eastern Uplift and Wulungu Depression. The heat flow contribution of sediments can warm effectively the deep and ultra-deep source rocks in Junggar Basin and promote maturity evolution of the organic matter. For example, the temperature in the bottom of J1s source rock (5300 m) in Well Pc2 was 7.3 ℃ higher in the condition of considering the heat production of sediments. Obviously, it was very significant to the hydrocarbon resource evaluation and exploration targets optimization in Junggar Basin which was characterized by low geothermal gradient and deep reservoirs.
Key words: Heat production rate     Heat flow of sedimentary formation     GR logging     Junggar Basin     Deep and ultra-deep formation    

1 引言

地壳岩石中含有60多种不稳定核素,这些自然放射性核素衰变释放出的热能是地球内部的主导热源.但并不是所有的放射性元素都对生热有贡献,放射性元素生热必须具备三个条件:①具有足够的丰度;②生热量大;③半衰期与地球年龄相当.在地球发展的整个历史时期中,能为地球提供大量热能的放射性同位素只有铀、钍和钾这三种元素.放射性元 素的时空分布对地球内部温度场的影响很大,地表 热流中来自放射性元素生热的热流贡献可以达到 30%~40% 甚至50%(邦特巴思,1988邱楠生等,2004).

单位体积岩石中所含放射性元素在单位时间内衰变所释放的能量称为岩石放射性生热率,它是描述岩石热物理性质的基本参数之一,也是研究一个地区深部热状态、岩石圈热-流变学结构等深部物理特征以及构造-热演化等动力学过程必不可少的有效参数(胡圣标等,1994Hu et al., 1994Wang et al., 2000汪集旸等,2001Shi et al., 2002Liu et al., 2004He et al., 2008).对于含油气盆地,沉积层生热率数据的准确获取,对盆地现今地温场计算、热史恢复、烃源岩热演化模拟和油气资源评价等都具有重要意义(汪集旸和汪辑安,1986王良书等, 19911996何丽娟,1999He et al., 20012002).

目前,确定地壳岩石放射性生热率的方法主要包括三种:一是通过地球化学方法(ICP-MS和XRF)直接测量岩石样品中放射性产热元素U、Th、K元素的含量,结合岩石密度测量结果,采用一些学者提出的经验公式来计算岩石样品放射性生热率(Birch,1954Rybach,1976Wollenberg and Smith, 1987),这种方法获得的数据质量高,是目前最常用的确定岩石生热率的方法,但需要在研究区采集相当数量的有代表性的岩样,且测试耗时较长,成本较高;二是根据地震波速估算岩石生热率(王良书等,1996Rybach,1986),这种方法获取的数据质量不高,主要用于地壳深部(中、下地壳)岩石生热率的估算;三是采用地球物理测井方法直接深入钻孔测量井壁岩石中天然放射性强度(自然伽马测井)或U、Th、K元素的含量(自然伽马能谱测井)换算成岩石生热率,最常见的是通过自然伽马(GR)与生热率(A)的相关关系,将自然伽马测井值直接换算成岩石生热率(Rybach,1986Bucker et al., 1996骆淼等,2008Bucker et al., 2001Norden and Forster, 2006).在GR-A相关关系可靠的基础上,这种方法获得的数据质量依然很高,且方便快捷,弥补了地球 化学方法耗时长、成本高和需要系统岩石样品的弊端.

作者收集了准噶尔盆地不同构造单元16口代表性钻孔自然伽马测井曲线,采用不同学者提出的GR-A经验关系式,分别计算了6120个沉积层生热率,通过与实测生热率的对比,确定了适用于研究区的GR-A经验关系,据此建立了准噶尔盆地地层生热率柱.此外,结合地震转换波探测揭示的沉积层结构,计算了研究区沉积层生热的热流贡献,以盆参2井为例定量分析了沉积层热流的增温效应.本次工作是关于准噶尔盆地沉积层生热率及其热流贡献的首次系统性研究,是研究区现今地温场和岩石圈热结构最新研究进展的一部分,该研究成果将为盆地深层—超深层烃源岩热史研究、热演化模拟和油气资源评价提供可靠的地热学依据.

2 地质背景

准噶尔盆地位于新疆北部,处于东经81°~92°、北纬43°~48°之间,是我国西部重要的含油气盆地. 盆地平面形状呈南宽北窄的近三角形,东西长700 km,南北宽370 km,面积13.4×104 km2,沉积层最大厚度可达15 km(陈业全和王伟锋,2004).从现今大地构造位置上看,准噶尔盆地位于哈萨克斯坦板块、西伯利亚板块及塔里木板块的三角地带,是哈萨克斯坦板块的东延部分,为三面被古生代缝合线包围的晚石炭世到第四纪发展起来的典型的多旋回叠合盆地(陈发景等,2005张朝军等,2006),具有由早古生界褶皱基底和前寒武系结晶基底组成的“双层基底”结构特点(邵学钟等,2008赵俊猛等,2008),经历了多期构造-热演化(潘长春等,1997Wang et al., 2001Qiu et al., 2005).以整体构造演化特征为背景,根据盆地各区构造演化和含油气方面的差异性,准噶尔盆地可以划分为6个一级构造单元共44个二级构造单元(杨海波等,2004),如图 1所示.

图 1 准噶尔盆地构造单元划分(据杨海波等,2004)与计算沉积层生热率钻井、地震转换波测线分布(据邵学钟等,2008) I—乌伦古坳陷(Wulungu Depression);Ⅱ—陆梁隆起(Luliang Uplift);Ⅲ—中央坳陷(Central Depression); Ⅳ—西部隆起(Western Uplift);Ⅴ—北天山山前坳陷(Southern Depression);Ⅵ—东部隆起(Eastern Uplift). Fig. 1 Distribution of structural units(after Yang et al., 2004),the calculated boreholes of heat production rate and seismic converted wave survey line(after Shao et al., 2008)in Junggar Basin

关于现今地温场和岩石圈热结构最新研究表明(饶松等,2013),准噶尔盆地现今地温梯度平均21.3± 3.7 ℃·km-1,大地热流平均42.5±7.4 mW·m-2,表现为低地温梯度、低大地热流的“冷”盆特征.准噶尔盆地大地热流与地温梯度分布规律基本一致,主要受控于基底的构造形态,东部隆起最高,陆梁隆起次之,乌伦古凹陷、中央凹陷和西部隆起较低,北天山山前坳陷最低.准噶尔盆地地壳热流介于18.8~ 26.0 mW·m-2,地幔热流介于16.5~23.7 mW·m-2,壳幔热流比值变化于0.79~1.58,属于典型的“冷壳冷幔”型热结构.

3 原理与方法 3.1 自然伽马测井原理

自然伽马测井是最常见的放射性测井方法之一,记录的是地层中所含放射性元素衰变时放出的伽马射线强度,它既可以在套管井中进行测量,也可以在空井和油基泥浆井中进行测量,由于其广泛的应用范围使其成为最常见的地球物理测井方法之一(张胜业和潘玉铃,2004).目前,自然伽马测井不仅在油气勘探中用于划分岩性、确定泥质含量、地层对 比,而且在古环境、古气候恢复方面日益受到关注(宋子齐等,2009李建红和周伦先,2008陈中红等,2004).

自然伽马测井的原理是:γ射线探测器探测到地层的γ射线,并将γ射线变换成电脉冲信号(每一道γ射线变换成一个电脉冲信号),此脉冲信号送入井下的放大器进行放大,然后送入地面的放大器再次进行放大(其原因是脉冲信号通过电缆之后会有所衰减).由于脉冲信号中混合一些干扰信号,需经过鉴别器进行鉴别,排除干扰,将一些有畸变的信号送入整形器进行整形.归一后的波形送入计数率计电路,此电路把脉冲信号变换成与单位时间内脉冲数成正比的电位差,记录仪连续记录电位差,最后得到自然伽马测井曲线.自然伽马测井最常 用的单位是API(美国石油学会—American Petroleum Institute推荐使用).

3.2 利用自然伽马测井计算生热率的理论基础和方法

自然伽马测井记录的是沉积地层中来自于U、Th和K总的伽马强度,其对生热率A/(μW·m-3)的贡献可以用下式描述(Rybach,1976):

式中,ρ是岩石密度(单位为kg·m-3);CU、CTh、CK分别为放射性元素U、Th、K的含量,单位分别为ppm(×106)、ppm(×106)、%.

由于自然伽马测井仪对U、Th和K的敏感性与方程(1)中相应的生热系数具有相似的比例,且地壳岩石中Th/U、K/U比例是恒定的,因此在特定的钻孔深度间隔内,自然伽马测井计数与生热率之间必定存在着简单的线性关系,这便构成了利用自然伽马测井直接计算岩石生热率的理论基础.

1986年,Rybach通过实测生热率与自然伽马计数的对比,最早提出GR-A的线性关系:

1996年,Bucker and Rybach(1996)基于大量钻孔岩心实测资料,通过最小二乘法拟合,对GR-A的线性关系做了进一步修正,这也是目前使用最为广泛的经验关系式(Bucker et al., 2001Norder and Forster, 2006):

骆淼等(2008)通过对中国大陆科学钻探(CCSD)主孔100~5000 m自然放射性测井(自然伽马和自然伽马能谱)与实测生热率数据的相关分析,建立了新的GR-A线性关系:

3.3 研究区GR-A经验关系

截止到目前,准噶尔盆地沉积层生热率研究极少,仅邱楠生(2002)实测了盆地8口钻孔共45个岩心样品的生热率.由于取芯深度和测井深度之间可能存在的差异,以及自然伽马测井中统计涨幅现象的存在,因此很难通过仅有的45个实测生热率值和自然伽马测井计数来直接建立或者判识研究区GR-A的经验关系.

收集了准噶尔盆地16口代表性钻孔的自然伽马测井曲线,其测量间隔均为0.125 m.图 2为分别采用以上不同经验关系式计算结果的对比,三者变化趋势完全一样,只是计算获得的生热率值存在差异,其中利用Bucker and Rybach(1996)经验关系式计算得到的生热率最大,Rybach(1986)次之,骆淼等(2008)最小;对于生热率值较大的岩性段(如莫深1井6783~6800 m泥岩段、石120井2452~2476 m泥岩段),三者差异较大,而对于生热率值较小的岩性段(如莫深1井6890~6900 m凝灰岩段、石120井2486~2498 m砂砾岩段),三者差异较小.

图 2 准噶尔盆地莫深1井6750~6900 m(a)和石120井2450~2600 m(b)自然伽马与计算生热率曲线 Fig. 2 The curves of GR logging and calculated heat production rate in 6750~6900 m interval of Well Ms1(a) and 2450~2600 m interval of Well S120(b)in Junggar Basin

在保证数据有效分析的前提下,为了提高数据分析的效率,作者对自然伽马测井曲线进行1/80抽 析,即每10 m取一个自然伽马测井值,如对于莫深1井,分别提取500、510、520 m,…,7480、7490、7500 m 深度共701个自然伽马测井值.然后,分别采用以上不同经验关系式,各计算了16口钻孔共6120个沉积层生热率.准噶尔盆地沉积层生热率计算值和实测值的统计对比表明,利用Bucker and Rybach(1996)经验关系式计算得到的6120个生热率与45个实测生热率,具有非常相近的平均值和标准差(前者为1.179±0.339 μW·m-3,后者为1.137±0.331 μW·m-3邱楠生,2002),如表 1所示,以及最相似的频率分布特征(图 3).因此,Bucker and Rybach(1996)经验关系式最适用于准噶尔盆地沉积层生热率的计算.

图 3 准噶尔盆地利用自然伽马测井计算的沉积层生热率与实测生热率统计直方图 (a)为利用Rybach(1986)经验关系式计算的生热率;(b)为利用Bucker and Rybach(1996)经验关系式计算的生热率; (c)为利用骆淼等(2008)经验关系式计算的生热率;(d)为邱楠生(2002)实测沉积层生热率. Fig. 3 Histograms of measured and calculated heat production rate by different GR—A empirical relationships (a),(b) and (c)is calculated heat production rate by GR—A empirical relationship of Rybach(1986)Bucker and Rybach(1996) and Luo et al.(2008),respectively.(d)is measured heat production rate by Qiu(2002).
4 计算结果

利用Bucker and Rybach(1996)经验关系式计算的6120个沉积层生热率统计结果如图 3b表 1所示.准噶尔盆地沉积层生热率值相对比较集中,大致变化于0.2~2.2 μW·m-3之间,主体介于0.8~ 1.6 μW·m-3,占计算值的81%.从统计结果来看,准 噶尔盆地沉积层平均生热率为1.179±0.339 μW·m-3. 这一计算结果与邱楠生(2002)45个实测生热率平均值1.137±0.331 μW·m-3非常相近,与倪守城等(1999)对新疆北部地区岩石生热率研究结果也基本一致.总体上,准噶尔盆地沉积层平均生热率与塔里木盆地1.196±0.886 μW·m-31)接近,低于西部柴达木盆地(邱楠生,2002),与中国东部地区特别是东南浙、闽、粤等地相比则更低(赵平等,1995).

表 1 利用不同GR-A经验关系计算的生热率与实测生热率的对比 Table 1 Comparisons of measured and calculated heat production rate by different GR-A empirical relationships
图 4 准噶尔盆地代表性钻井沉积层生热率随深度的变化 Fig. 4 The change of calculated heat production rate of sediments with depth of representative boreholes in Junggar Basin

1)《中国陆域海相沉积盆地现今地温场与岩石圈热结构》,北京:中国科学院地质与地球物理研究所研究报告,2011年11月 图 4所示为准噶尔盆地北天山山前坳陷西湖1、霍10井,中央坳陷莫深1、玛13、玛东1井,西部隆起沙门1、中佳1、金龙6、风城1井,陆梁隆起陆9、石120、陆东1、滴中1井,乌伦古坳陷伦参1井,和东部隆起石树1、台8井共16口钻孔的生热率随深度变化,其统计结果见表 2.由图 4表 2可以看出,准噶尔盆地沉积层生热率总体上具有随着深度增加(时代变老)逐渐降低的趋势,如莫深1井上组合沉积层生热率相对较大,下第三系生热率平均为1.477 μW·m-3;中组合沉积层生热率次之,其中白垩系生热率平均为1.534 μW·m-3,侏罗系生热率平均为1.157 μW·m-3;下组合沉积层生热率最低,其中三叠系生热率平均为1.076 μW·m-3,二叠系生热率平均为0.998 μW·m-3,石炭系生热率 平均为0.872 μW·m-3.再如,金龙6井白垩系平均生热率为1.030 μW·m-3,侏罗系为0.823 μW·m-3,三叠系为0.921 μW·m-3、二叠系和石炭系生热率 很低,分别仅为0.529 μW·m-3和0.440 μW·m-3.

此外,从沉积层生热率的离散程度来看,准噶尔盆地下组合要远远大于中—上组合,中—上组合生热率随深度基本呈线性展布,生热率值突变并不多见,而下组合中沉积层生热率“尖峰状”和“箱状”变化更加普遍,如风城1井4100~4600 m井段、台8 井2650~2880 m井段频繁地出现多个高生热率“尖峰”,陆东1井3900~4400 m井段、4810~5050 m 井段出现了高生热率“箱”.这可能反映了盆地石炭纪—二叠纪构造—热活跃期沉积环境和岩石类型的频繁变化.

表 2 准噶尔盆地代表性钻井沉积层生热率统计 Table 2 The calculated heat production rate of sedimentary formation of representative boreholes in Junggar

对6120个计算生热率值按照层位进行统计平均,建立了准噶尔盆地沉积层生热率柱(图 5).总体上看,准噶尔盆地沉积层生热率随着时代变老,生热率呈现 出递减 趋势:第四系平均生热率最高,为1.309 μW·m-3;新近系平均生热率次之,为1.250 μW·m-3; 古近 系和白垩系平均生热率相当,均为1.238 μW·m-3;侏罗系、三叠系平均生热率相对较低,分别为1.141 μW·m-3和1.071 μW·m-3;二叠系和石炭系生热率最低,平均分别仅为0.997 μW·m-3和0.881 μW·m-3. 准噶尔盆地沉积层生热率随着时代变老(埋深增加)逐渐降低的现象,可能体现了地壳放射性元素U、Th和K40向地壳浅部富集的普遍规律.

图 5 准噶尔盆地不同沉积层生热率平均值柱状图 Fig. 5 The average heat production rate columns of the sedimentary formation in Junggar Basin
5 讨论 5.1 准噶尔盆地沉积层热流贡献

地表热流(q)是地球内热在地表的直接反映,它由两部分组成,即地壳岩石中U、Th、K等放射性元素蜕变产生的热量(地壳热流,qc)和来自地球深部的热量(地幔热流,qm).地壳热流中,来自沉积层放射性生热的热流贡献(沉积层热流,qs)取决于两个因素:沉积层厚度和沉积层放射性元素丰度(沉积层生热率).沉积层热流可以用以下公式计算:

沉积层热流在地壳热流和地表热流中的比例可以分别用f1和f2表示:

式中,qi为沉积层中第i层放射性生热的热流贡献,Ai为沉积层中第i层的生热率,Zi为沉积层中第i层的厚度,n为沉积层的层数.

在准噶尔盆地沉积层生热率柱建立的基础上,结合邵学钟等(2008)采用天然地震转换波测深方法揭示的沉积层结构,计算了4条转换波测线共250 个测点(测点位置如图 1所示)沉积层热流值,通过Kringing插值方法,编制了准噶尔盆地沉积层热 流等值线图(图 6).结果表明,准噶尔盆地沉积层 热流平均值为7.9 mW·m-2,约占平均地壳热 流24.6 mW·m-2的29.2%和平均大地热流42.5 mW·m-2的19.6%.由此看来,准噶尔盆地沉积层热流是地壳热流和地表热流的重要组成部分,一定程度上弥补了盆地深部“冷”幔结构的较低热流.区域上,准噶尔盆地沉积层热流分布特征与盆地沉积层厚度大体一致,表现为中央坳陷最大,平均大于13 mW·m-2,阜康凹陷和沙湾凹陷甚至超过15 mW·m-2;北天山山前冲断带变化较大,霍玛吐背斜带普遍大于13 mW·m-2,阜康断裂带介于 7~13 mW·m-2,齐古断褶带变化于5~13 mW·m-2,四棵树凹陷普遍小于7 mW·m-2;陆梁隆起和西部 隆起次之,介于3~7 mW·m-2,平均约为5 mW·m-2; 东部隆起和乌伦古坳陷沉积层热流很低,普遍小于5 mW·m-2.

图 6 准噶尔盆地沉积层热流等值线图 Fig. 6 The contour map of heat flow contribution of the sedimentary formation in Junggar Basin

图 7为准噶尔盆地地震转换波测线Ⅱ和Ⅲ沉积层热流及其在地壳热流和地表热流的比例.测线Ⅱ近东西走向,经过盆地西部隆起车排子凸起、红车断裂带,中央坳陷沙湾凹陷、莫南凸起、阜康凹陷,东部隆起北三台凸起、沙奇凸起、石树沟凹陷、黄草湖凸起等共9个二级构造单元;测线Ⅲ近南北走向,依次穿过北天山山前冲断带齐古断褶带、霍玛吐背斜带,中央坳陷莫南凸起、沙湾凹陷、莫索湾凸起、盆1井西凹陷、达巴松凸起,陆梁隆起夏盐凸起、三个泉凸起、英西凹陷、石英滩凸起,乌伦古坳陷索索泉凹陷、红岩断阶带等共13个二级构造单元.由图 7a可以明显看出,测线Ⅱ从西部隆起→中央坳陷→东部隆起,沉积层热流从低→高→低,沉积层热流在地壳热流和地表热流的比例变化规律与沉积层热流一致,表现为坳陷高、隆起低,其中中央坳陷沉积层热流约15 mW·m-2,占地壳热流比例超过50%,在地表热流中的贡献超过40%.图 5b展示的测线Ⅲ沉积层热流贡献变化规律与测线Ⅱ一致,从北天山山前坳陷→中央坳陷→陆梁隆起→乌伦古坳陷,沉积层热流贡献从低→高→低→低变化.

图 7 准噶尔盆地地震转换波测线Ⅱ和Ⅲ沉积层热流贡献 (a1)、(b1)分别为测线Ⅱ和Ⅲ沉积层热流;(a2)、(b2)分别为测线Ⅱ和Ⅲ沉积层热流在地壳热流中的比值; (a3)、(b3)分别为测线Ⅱ和Ⅲ沉积层热流在地表热流中的比值. Fig. 7 The heat flow contribution of the sedimentary formation of seismic converted wave survey line Ⅱ and Ⅲ in Junggar Basin (a1),(b1)is sedimentary formation heat flow of survey line Ⅱ and Ⅲ,respectively.(a2)、(b2)is the proportion of sedimentary formation heat flow in crustal heat flow of survey line Ⅱ and Ⅲ,respectively.(a3)、(b3)is the proportion of sedimentary formation heat flow in terrestrial heat flow of survey line Ⅱ and Ⅲ,respectively.
5.2 准噶尔盆地沉积层放射性生热的增温意义

以上对不同沉积层放射性生热率的统计及盆地沉积层热流贡献计算结果表明,准噶尔盆地沉积层 放射性生热对大地热流的贡献平均为7.9 mW·m-2,约占地表热流值42.5 mW·m-2的19.6%.这种增温效应直观地体现在地热梯度上,使得现今的平均地热梯度21.3 ℃·km-1中约有3.6 ℃·km-1是由沉积层放射性生热引起的(准噶尔盆地平均岩石热导率为2.209 W·(m·K)-1(饶松等,2013).

准噶尔盆地新生代以来独特的大地构造活动性塑造了盆地地壳厚度大、地温梯度小、主体油气藏埋深大的油气地质背景,沉积层放射性生热的增温效应显然会增加深层(埋深超过4500 m)和超深层(埋深大于6000 m)烃源层的受热温度,促进烃源岩有机质热演化,这对于准噶尔盆地以低地温梯度、低大地热流为特征的“冷盆”的油气资源评价、勘探目标优选无疑意义重大.

作者以准噶尔盆地中央坳陷盆参2井(钻井位置如图 1所示)为例,定量计算沉积层放射性生热的热效应.盆参2井完钻井深5300 m,完钻层位为下侏罗统三工河组(J1s,未穿),自上而下钻遇第四系西域组(Q1x),新近系独山子组(N2d)、塔西河组(N1t)、沙湾组(N1s),古近系安集海河组(E3a)、紫泥泉子组(E1-2z)、上白垩统东沟组(K2d)、下白垩统吐谷鲁群(K1tg),上侏罗统齐古组(J3q),中侏罗统头屯河组(J2t)、西山窑组(J2x),下侏罗统三工河组(J1s).其中齐古组与头屯河组不整合接触,剥蚀厚度约200 m;紫泥泉子组与东沟组不整合接触,剥蚀厚度约260 m(李平平等,2006周路等,2007).钻探结果和油气源对比研究均表明,下侏罗统三工河组和八道湾组(J1b,三工河组以下,盆参2井未钻遇)为区域性有效烃源岩,现今仍处于低成熟—成熟演化阶段(陆鹏庆,2010丁安娜等,1996).

在考虑和忽略沉积层放射性生热的两种情况下,采用一维非稳态、有内源热的热传导方程,分别计算盆参2井三工河组底部温度史和现今温度随深度的变化:

计算过程中,平均地表温度取15 ℃,沉积层热导率和地表热流据饶松等(2013)最新研究成果,盆地热流演化史参考Qiu et al.(2005),整个计算过程在PRAI公司开发的BasinMod 1-D软件中完成.

图 8所示为盆参2井沉积层放射性生热的增温效应.从三工河组底部温度史和现今温度随深度的变化的计算结果可以看出:(1)在考虑和忽略沉积层放射性生热的两种情况下,计算的三工河组底部温度差异从古至今逐渐增大,现今二者差异最大;(2)在考虑和忽略沉积层放射性生热的两种情况下,计算的现今温度差异随深度增加逐渐增大,在井底 5300 m深度二者现今温度分别为115.3 ℃和108.0 ℃,相差达7.3 ℃.按照以上规律,越往深层,沉积层放射性生热的增温效应将越明显,约在8000 m深度,沉积层放射性生热的增温效应将超过10 ℃.按照油气有机成因理论,烃源岩热演化与时间成线性关系,而与温度呈指数关系,即温度每增加10 ℃,反应速率将增加1倍.因此,沉积层放射性生热的增温效应能够有效促进准噶尔盆地深层—超深层烃源岩热演化,很显然,这对于地温梯度小、主体油气藏埋深大的准噶尔盆地的油气资源评价和勘探目标优选具有重要意义.

图 8 准噶尔盆地盆参2井沉积层放射性生热的增温效应 (a)J1s底部温度史;(b)现今温度随深度的变化. Fig. 8 The warming effect of sedimentary formation heat flow and influence to thermal evolution of the J1s source rocks of Well Pc2 in Junggar Basin (a)The temperature history in the bottom of J1s;(b)The change of present geotemperature with depth.
6 结论与认识

沉积层放射性生热的热流贡献(沉积层热流)是 地表热流的重要组成部分,能够有效促进中国西部以低地温梯度、低大地热流为特征的“冷盆”的深层—超深层烃源岩增温和热演化.本文利用自然伽马测井曲线计算了准噶尔盆地沉积层生热率,并据此计算了研究区沉积层放射性生热的热流贡献,得到如下结论和认识:

(1)通过生热率计算值与实测值的统计对比,确定了Bucker and Rybach(1996)GR-A经验关系式同样适用于准噶尔盆地沉积层生热率的计算;

(2)准噶尔盆地沉积层平均生热率为1.179±0.339 μW·m-3,总体上随着时代变老,沉积层生热率呈现出递减趋势,第四系最高,石炭系最低;

(3)准噶尔盆地沉积层热流平均为7.9±4.9 mW·m-2,约占地壳热流的29.2%和地表热流的19.6%,区域上与盆地沉积层厚度大体一致,表现为中央坳陷最高,北天山山前冲断带变化较大,陆梁隆起和西部隆起次之,东部隆起和乌伦古坳陷沉积层热流最低;

(4)沉积层热流能够有效提高深层—超深层烃源层的受热温度,促进有机质热演化.沉积层热流的增温效应,对于地温梯度小、主体油气藏埋深大的准噶尔盆地油气资源评价和勘探目标优选具有重要意义.

致 谢 钻孔自然伽马测井等资料收集中得到中国石油勘探开发研究院、新疆油田勘探开发研究院和新疆油田勘探处的大力帮助,两位匿名审稿专家对本文提出了细致、宝贵的修改意见,在此一并致谢!

参考文献
[1] Birch F. 1954. Heat from radioactivity. Nuclear Geology. New York: Wiley, 148-174.
[2] Bücker C, Rybach L. 1996. A simple method to determine heat production from gamma-ray logs.   Marine and Petroleum Geology, 13(4): 373-375.
[3] Bucker C, Jarrard R D, Wonik T. 2001. Downhole temperature, radiogenic heat production, and heat flow from the CRP-3 drillhole, Victoria Land Basin, Antarctica.   Terra Antartica, 8(3): 151-160.
[4] Buntebarth G. 1988. Geothermics: An Introduction (in Chinese). Yi Z X, Xiong L P, Trans. Beijing: Seismological Press, 15-21.
[5] Chen F J, Wang X W, Wang X W. 2005. Prototype and tectonic evolution of the Junggar Basin, northwestern China. Earth Science Frontiers (in Chinese), 12(3): 77-89.
[6] Chen Y Q, Wang W F. 2004. Structural evolution and pool-forming in Junggar Basin. Journal of the University of Petroleum, China (Edition of Natural Science) (in Chinese), 28(3): 4-8.
[7] Chen Z H, Zha M, Jin Q. 2004. Application of natural gamma ray logging and natural gamma spectrometry logging to recovering paleoenvironment of sedimentary basin. Chinese J. Geophys.   (in Chinese), 47(6): 1145-1150.
[8] Ding A N, Hui R Y, Meng Q X, et al. 1996. Characteristics of source rocks and their hydrocarbon formation in Jurassic system, Junggar Basin, Xinjiang. Petroleum Exploration and Development (in Chinese), 23(3): 11-18.
[9] He L J. 1999. Mutiple tectono-thermal modeling of Liaohe Basin in the Cenozoic. Chinese J. Geophys.   (in Chinese), 42(1): 62-68.
[10] He L J, Wang K L, Xiong L P, et al. 2001. Heat flow and thermal history of the South China Sea.   Physics of the Earth and Planetary Interiors, 126(3-4): 211-220.
[11] He L J, Xiong L P, Wang J Y. 2002. Heat flow and thermal modeling of the Yinggehai Basin, South China Sea.   Tectonophysics, 351(3): 245-253.
[12] He L J, Hu S B, Huang S P, et al. 2008. Heat flow study at the Chinese Continental Scientific Drilling site: Borehole temperature, thermal conductivity, and radiogenic heat production. Journal of Geophysical Research, 113(B2), B02404, doi: 10.1029/2007JB004958.
[13] Hu S B, Wang J Y, Wang Y H. 1994. Deep temperature and lithospheric thickness long the eastern segment of the Heishui—Quanzhou geotraverse. Chinese J. Geophys. (in Chinese), 37(3): 330-337.
[14] Hu S B, Wang J Y. 1994. Crustal heat production and mantle heat flow in southeast China.   Science in China (Series B), 37(10): 1252-1263.
[15] Li J H, Zhou L X. 2008. Response features of gamma ray log and its application in glutenite of Dongying depression.   Journal of Oil and Gas Technology (in Chinese), 30(1): 88-91.
[16] Li P P, Zou H Y, Hao F, et al. 2006. Restoration of eroded strata thickness in Cretaceous/Jurassic unconformity in hinterland of Junggar Basin. Acta Petrolei Sinica (in Chinese), 27(6): 34-38.
[17] Liu S W, Wang L S, Li C, et al. 2004. Thermal-rheological structure of lithosphere beneath the northern flank of Tarim Basin, western China: Implications for geodynamics.   Science in China Series D: Earth Sciences, 47(7): 659-672.
[18] Lu P Q. 2010. Geochemical characteristics of the source rocks in the central depression of Junggar Basin. Journal of Oil and Gas Technology (in Chinese), 32(6): 216-218.
[19] Luo M, Pan H P, Zhao Y G, et al. 2008. Natural radioactivity logs and interpretation from the CCSD main hole.   Earth Science-Journal of China University of Geosciences (in Chinese), 33(5): 661-671.
[20] Ni S C, Man F S, Wang J R, et al. 1999. Characteristics of heat production distribution in northern Xinjiang. Journal of University of Science and Technology of China (in Chinese), 29(4): 408-414.
[21] Norden B, Forster A. 2006. Thermal conductivity and radiogenic heat production of sedimentary and magmatic rocks in the Northeast German Basin.   AAPG Bulletin, 90(6): 939-962.
[22] Pan C C, Zhou Z Y, Fan S F, et al. 1997. Thermal history of Junggar Basin. Geochimica (in Chinese), 26(6): 1-7.
[23] Qiu N S. 2002. Characters of thermal conductivity and radiogenic heat production rate in basins of Northwest China.   Chinese Journal of Geology (in Chinese), 37(2): 196-206.
[24] Qiu N S, Hu S B, He L J. 2004. Principles and Applications on Thermal Regime of Sedimentary Basins (in Chinese). Beijing: Petroleum Industry Press, 22-24.
[25] Qiu N S, Zha M, Wang X L, et al. 2005. Tectono-thermal evolution of the Junggar Basin, NW China: constraints from Ro and apatite fission track modelling.   Petroleum Geoscience, 11(4): 361-372.
[26] Rao S, Hu S B, Zhu C Q, et al. 2013. The characteristics of heat flow and lithospheric thermal structure in Junggar Basin, northwest China. Chinese J. Geophys. (in Chinese), 56(8): 2760-2770.
[27] Rybach L. 1976. Radioactive heat production in rocks and its relation to other petrophysical parameters.   Pure and Applied Geophysics, 114(2): 309-317.
[28] Rybach L. 1986. Amount and significance of radioactive heat sources in sediments.   Collection Colloques et Séminares, 44: 311-322.
[29] Shao X Z, Zhang J R, Fan H J, et al. 2008. The basement structure in Junggar Basin: Deep-sounding by converted waves of earthquakes. Xinjiang Petroleum Geology (in Chinese), 29(4): 439-444.
[30] Shi X B, Zhou D, Qiu X L, et al. 2002. Thermal and rheological structures of the Xisha Trough, South China Sea.   Tectonophysics, 351(4): 285-300.
[31] Song Z Q, Li W F, Tang C J, et al. 2009. Dividing sedimentary facies and reservoir distributions by using natural potential and natural gamma ray logging curves. Progress in Geophys.   (in Chinese), 24(2): 650-656.
[32] Wang J Y, Wang J A. 1986. Mantle heat flow of Liaohe rifted basin in north China. Chinese J. Geophys.  (in Chinese), 29(5): 450-459.
[33] Wang J, Hu S B, Cheng B H, et al. 2001. Predication of the deep temperature in the target area of the China Continental Scientific Drhling. Chinese J. Geophys. (in Chinese), 44(6): 774-782.
[34] Wang L S, Li C, Shi Y S. 1991. Influence of radioactive heat generation within lithosphere on geothermal history of extensional basin.   Petroleum Geology & Experiment (in Chinese), 13(3): 257-264.
[35] Wang L S, Li C, Yang C. 1996. The lithospheric thermal structure beneath Tarim Basin, western China. Chinese J. Geophys.  (in Chinese), 39(6): 794-803.
[36] Wang L S, Li C, Liu F T, et al. 2000. Thermal-rheological structure of the lithosphere beneath two types of basins in eastern and western China.   Science in China Series D: Earth Sciences, 43(1): 200-207.
[37] Wang S J, He L J, Wang J Y. 2001. Thermal regime and petroleum systems in Junggar Basin, northwest China.   Physics of the Earth and Planetary Interiors, 126(3-4): 237-248.
[38] Wollenberg H A, Smith A R. 1987. Radiogenic heat production of crustal rocks: An assessment based on geochemical data.   Geophysical Research Letters, 14(3): 295-298.
[39] Yang H B, Chen L, Kong Y H. 2004. A novel classification of structural units in Junggar Basin.   Xinjiang Petroleum Geology (in Chinese), 25(6): 686-688.
[40] Zhang C J, He D F, Wu X Z, et al. 2006. Formation and evolution of multicycle superimposed basins in Junggar Basin. China Petroleum Exploration (in Chinese), 11(1): 47-58.
[41] Zhang S Y, Pan Y L. 2004. Principles on Applied Geophysics (in Chinese). Wuhan: China University of Geosciences Press, 420-430.
[42] Zhao J M, Huang Y, Ma Z J, et al. 2008. Discussion on the basement structure and property of northern Junggar basin. Chinese J. Geophys. (in Chinese), 51(6): 1767-1775.
[43] Zhao P, Wang J Y, Wang J A, et al. 1995. Characteristics of heat production distribution in SE China. Acta Petrologica Sinica (in Chinese), 11(3): 292-305.
[44] Zhou L, Zheng J Y, Lei D W, et al. 2007. Recovery of eroded thickness of the Jurassic of Chemo palaeouplift in Junggar Basin. Journal of Palaeogeography (in Chinese), 9(3): 243-252.
[45] 邦特巴思 G. 1988. 地热学导论. 易志新, 熊亮萍译. 北京: 地震出版社, 15-21.
[46] 陈发景, 汪新文, 汪新伟. 2005. 准噶尔盆地的原型和构造演化.   地学前缘, 12(3): 77-89.
[47] 陈业全, 王伟锋. 2004. 准噶尔盆地构造演化与油气成藏特征.   石油大学学报(自然科学版), 28(3): 4-8.
[48] 陈中红, 查 明, 金 强. 2004. 自然伽玛及自然伽玛能谱测井在沉积盆地古环境反演中的应用. 地球物理学报, 47(6): 1145-1150.
[49] 丁安娜, 惠荣耀, 孟仟祥等. 1996. 准噶尔盆地侏罗系烃源岩及油气形成特征.   石油勘探与开发, 23(3): 11-18.
[50] 何丽娟. 1999. 辽河盆地新生代多期构造热演化模拟. 地球物理学报, 42(1): 62-68.
[51] 胡圣标, 汪集旸, 汪屹华. 1994. 黑水—泉州地学断面东段深部温度与岩石层厚度. 地球物理学报, 37(3): 330-337.
[52] 李建红, 周伦先. 2008. 东营凹陷砂砾岩自然伽马测井响应特征研究及应用.   石油天然气学报, 30(1): 88-91.
[53] 李平平, 邹华耀, 郝 芳等. 2006. 准噶尔盆地腹部白垩系/侏罗系不整合地层剥蚀厚度的恢复方法.   石油学报, 27(6): 34-38.
[54] 陆鹏庆. 2010. 准噶尔盆地腹部地区烃源岩地球化学特征.   石油天然气学报, 32(6): 216-218.
[55] 骆淼, 潘和平, 赵永刚等. 2008. 中国大陆科学钻探主孔自然放射性测井及其解释.   地球科学: 中国地质大学学报, 33(5): 661-671.
[56] 倪守斌, 满发胜, 王兆荣等. 1999. 新疆北部地区岩石生热率分布特征.   中国科学技术大学学报, 29(4): 408-414.
[57] 潘长春, 周中毅, 范善发等. 1997. 准噶尔盆地热历史. 地球化学, 26(6): 1-7.
[58] 邱楠生. 2002. 中国西北部盆地岩石热导率和生热率特征.   地质科学, 37(2): 196-206.
[59] 邱楠生, 胡圣标, 何丽娟. 2004. 沉积盆地热体制研究的理论与应用. 北京: 石油工业出版社, 22-24.
[60] 饶松, 胡圣标, 朱传庆等. 2013. 准噶尔盆地大地热流特征与岩石圈热结构. 地球物理学报, 56(8): 2760-2770.
[61] 邵学钟, 张家茹, 范会吉等. 2008. 准噶尔盆地基底结构的地震转换波探测. 新疆石油地质, 29(4): 439-444.
[62] 宋子齐, 李伟峰, 唐长久等. 2009. 利用自然电位与自然伽马测井曲线划分沉积相带及储层分布.   地球物理学进展, 24(2): 650-656.
[63] 汪集旸, 汪缉安. 1986. 辽河裂谷盆地地幔热流. 地球物理学报, 29(5): 450-459.
[64] 汪集旸, 胡圣标, 程本合等. 2001. 中国大陆科学钻探靶区深部温度预测. 地球物理学报, 44(6): 774-782.
[65] 王良书, 李成, 施央申. 1991. 岩石圈生热效应对拉张盆地地热史的影响.   石油实验地质, 13(3): 257-264.
[66] 王良书, 李成, 杨春. 1996. 塔里木盆地岩石层热结构特征. 地球物理学报, 39(6): 794-803.
[67] 杨海波, 陈磊, 孔玉华. 2004. 准噶尔盆地构造单元划分新方案.   新疆石油地质, 25(6): 686-688.
[68] 张朝军, 何登发, 吴晓智等. 2006. 准噶尔多旋回叠合盆地的形成与演化.   中国石油勘探, 11(1): 47-58.
[69] 张胜业, 潘玉玲. 2004. 应用地球物理学原理. 武汉: 中国地质大学出版社, 420-430.
[70] 赵俊猛, 黄英, 马宗晋等. 2008. 准噶尔盆地北部基底结构与属性问题探讨. 地球物理学报, 51(6): 1767-1775.
[71] 赵平, 汪集旸, 汪缉安等. 1995. 中国东南地区岩石生热率分布特征.   岩石学报, 11(3): 292-305.
[72] 周路, 郑金云, 雷德文等. 2007. 准噶尔盆地车莫古隆起侏罗系剥蚀厚度恢复.   古地理学报, 9(3): 243-252.