1、引 用 格 式：，（）：毕然，朱焕来，杜先利 基于自然伽马测井的徐家围子断陷沉积层生热率研究 甘肃科学学报，（）：基于自然伽马测井的徐家围子断陷沉积层生热率研究毕然，朱焕来，杜先利，（东北石油大学地球科学学院，黑龙江 大庆 ；东北石油大学地热能研究中心，黑龙江 大庆 ）摘要以徐家围子断陷为研究对象，研究生热率特点及其对地温场的影响。利用 口重点钻井内 个自然伽马值，依据自然伽马与生热率（）之间的关系，计算其生热率的数值，建立岩性生热率的对比图，估算研究区放射性物质生热对大地热流的贡献以及对深部地层的增温效应。结果显示，徐家围子的平均生热率为（），岩性是生热率的主要影响因素，流纹岩的生热率最高，
2、平均值为，凝灰岩、安山岩、火山角砾岩的生热率平均值也较高，分别为、，而砾岩的生热率最低，为，砂岩的生热率平均值为，泥岩的生热率平均值为，煤的生热率平均值为。根据地层生热率，估算徐家围子沉积层生热贡献的热流为，约占地表总热流 的，沉积层生热对地温梯度的贡献约为，放射性生热对属于“热盆”的松辽盆地的地温场具有不容忽视的影响。关键词徐家围子断陷；生热率；自然伽马测井；沉积层热流；大地热流中图分类号：；文献标志码：文章编号：（）地球内部主要的热量来源于岩石中的放射性元素衰变释放出的大量热能。自然界中存在 多种极为不稳定的放射性元素，这些元素满足一定丰度、较高的产热量以及半衰期与地球年龄相当这个条件才能
3、明显影响地球内部的热量。而满足这种条件的放射性元素只有铀（）、钍（）、钾（），它们在地壳中所释放的热量对地表热流的贡献率高达 。岩石的放射性生热率（）是最主要的岩石热物性之一，是单位时间内放射性元素衰变而释放的大量热能。由于岩石的生热率等热物性参数对研究地球内部热情况、岩石圈的热结构、地表热流以及油气勘探开发等有重要的指导意义，所以近年来愈发重视对岩石的热物性参数的研究 。目前，有种方法可获得岩石中放射性物质的生热率：（）结合岩石的测量密度，利用地球化学的方法直接测量岩石中铀（）、钍（）、钾（）的含量，根据前人提出的经验公式来计算岩石的放射性物质的生热率 ，虽然这种方法确定的生热率数据质量较好
4、，但需要收集大量的岩石样品，所需要测试的时间久且花费的成本较高；（）利用地震波测量计算岩石中放射性物质的生热率，由于该方法得到的生热率数据质量较差，所以适用于中下层地壳放射性物质生热率的测量计算 ；（）通过地球物理测井的方法测得测井中岩石天然放射性强度（自然伽马测井）或者铀（）、钍（）、钾（）的含量（自然伽马能谱测井），再根据前人提出的自然伽马（）与岩石放射性物质生热率（）之间的经验公式计算得到岩石放射性物质的生热率 ，该方法相较于地球化学方法，得到的生热率数据质量较高，同时简第 卷第期 年月 甘 肃 科 学 学 报 收稿日期：；修回日期：作者简介：毕然（），女，江苏南京人，硕士研究生，研究方
5、向为地热资源评价。：通信作者：朱焕来，：明快速，节省成本。前人研究表明，徐家围子现今地温梯度在 ，平均为 ；大地热流在 之间，平均值为 ；安达凹陷、升平凸起、徐西凹陷、徐东斜坡、徐南凹陷的大地热流值分别为 、。总体来说，地温场的特征主要受基底起伏控制，凹陷区低，隆起区高。前人对徐家围子生热率研究较少，李志安 直接测量岩石样品中放射性产热元素、元素的含量，根据生热率和放射性元素之间的关系计算松辽盆地各地层的生热率，但是对于深部岩层无法直接取样计算。为计算徐家围子的岩石生热率，利用该研究区 口 井 的 自 然 伽 马 值，采 用 、等 及骆淼等 提出的种计算方法。通过与实测数据的对比确定适用于研究
6、区的 的经验关系，此外还计算沉积层放射性物质生热对大地热流的贡献，定量分析沉积层放射性物质对深部地层的增温效应。研究徐家围子沉积层生热率及其热流贡献，是研究该地区现今地温场特征的一部分，可为研究区深层超深层烃源岩热演化、大地热流、能源勘探等提供可靠的热物性参数。地质背景徐家围子断陷是松辽盆地主要的生气断陷之一，该断陷形成于晚侏罗纪早白垩世早期。徐家围子的内部构造主要可划分为：安达次坳、杏山次坳、薄荷台次坳、徐西凹陷和徐东斜坡带等，如图所示。徐家围子断陷是西断东超型箕型断陷，在沙河子组时期沉积后升西断裂带发生强烈活动，之后形成断弯褶皱，褶皱顶部遭受到强烈侵蚀，在营城组沉积之后褶皱再次发育，徐家围
7、子断陷内各个构造单元具有构造运动和火山活动双重成因机制。徐家围子断陷地层由下到上依次发育有石炭二叠系变质岩基底、白垩系下统火石岭组、沙河子组、营城组、登娄库组、泉头组、青山口组、姚家组和嫩江组，如图所示，发育冲积扇相、河流相、湖漫相、湖泊相、扇三角洲相、湖底扇相等类型的沉积相。方法与原理自然伽马测井计算生热率的理论基础是：射线探测器可以探测到地层岩石中、等放射性元素衰变发出的射线，探测器发生电离，将射线转换为电脉冲信号，并将此信号送到井下的前置放大器放大，然后再将其送到地面的放大器再次放大，最后送入地面上的控制面板中的计数率电路，该电路会将脉冲信号转换成单位时间内与脉冲数成正比的电位差，由记录
8、仪器记录连续的电位差，最终得到自然伽马测井曲线 。图徐家围子断陷构造 年，提出岩石生热率与自然伽马值的线性关系为.（）.）。（）年，等 通过大量的实测资料，对生热率和自然伽马值的线性关系做进一步的修改：.（）.）。（）年，骆 淼 等 利 用 中 国 大 陆 科 学 钻 探（）主孔自然放射性测井数据资料，计算出新的岩石生热率和自然伽马值线性关系为.（）.）。（）由于自然伽马测井仪对地层中放射性元素、的灵敏性和公式（）对应的元素含量具有相似的系数，并且地层中放射性元素 、的比值是恒定的，则在一定的钻孔深度间隔，放射性元素 的生 热 率 与 自 然伽马测 井 值 必 然 有 线 性关系。甘 肃 科
9、学 学 报 年第期注：层序级别和层序界面的划分依据 年大庆油田研究院研究成果。图徐家围子断陷地层综合柱状图 收集整理研究区 口钻井数据，选取其中 口重点井，运用公式（）（）种方法计算研究区的岩石生热率。以肇深 井为例，种方法计算得到的生热率如图所示。种计算公式的结果表明生热率的变化趋势是一致的，仅为计算的数据结果不同，式（）计算得到的岩石生热率值最大，式（）计算得到的岩石生热率值次之，式（）计算得出的岩石生热率最小。岩石生热率与岩石岩性密切相关，岩第 卷毕然等：基于自然伽马测井的徐家围子断陷沉积层生热率研究注：为根据式（）计算的生热率；为根据式（）计算的生热率；为根据式（）计算的生热率。图徐家
10、围子断陷肇深 井 曲线图 性不同，岩石生热率有所差异。计算结果为提高数据计算的准确性，每 取一个自然伽马值，如对肇深 井进行取样，分别取 、深度的伽马值，口井总共取 个自然伽马值。不同公式计算出的生热率值与实测值对比如表所列。表中公式（）计 算 得 到 的 生 热 率 的 平 均 值 为（）；公式（）计算得到的生热率的平均值为（）；公式（）计算得到的生热率的平均值为（）。徐家围子断陷生热率统计直方图如图所示。由图可以看出种公式计算得到的生热率数值均表不同公式计算出的生热率值与实测值对比 生热率种不同公式公式（）公式（）公式（）实测值 平均值（）标准值（）图徐家围子断陷生热率统计直方图 在之间，
11、占主体的 左右。而公式（）计算的得到生热率与实测值最为接近，所以公式（）计算得出的结果是可靠的。将公式（）计算得到的 个生热率数值按层位进行统计，得到不同层位的岩石生热率的平均值如图所示。由图可知，青一段的生热率值最高，为；而登一段的生热率数值最低，为。将公式（）计算得到的 个生热率数值按岩性进行统计，得到了不同岩性的生热率的平均值及甘 肃 科 学 学 报 年第期范围如表所列，生热率分布直方图如图所示。主要岩性的生热率分布范围显示，流纹岩的生热率主要分布在；火山角砾岩主要分布在 ；砾 岩 主 要 分 布 在 ；砂岩主要分布在；泥岩主要分布在；煤主要分布在；安山岩主要分布在。图徐家围子断陷不同地
12、层平均生热率 表不同岩性生热率统计 岩性取样个数生热率范围（）生热率平均值（）砾岩 砂岩 泥岩 煤层 火山角砾岩 流纹岩 安山岩 凝灰岩 图徐家围子断陷岩性生热率分布直方图 不同岩性生热率的差别较大，流纹岩的生热率最高，平均值为，凝灰岩、安山岩、火山角砾岩的生热率平均值也较高，分别为、，而砾岩的生热率最低，为 ，砂岩的生热率平均值为 ，泥岩的生热率平均值为 ，煤的生热率平均值为。讨论影响因素（）深度从姚二三段至青一段的生热率数值具有上升趋势，而从泉二段至登四段生热率数值下降，可见生热率和深度的变化并没有明显的线性关系。结合生热率随深度变化图，并参考黄少鹏 的观点，可知生热率和深度并不存在简单的
13、一元线性关系。（）岩性岩石中的自然伽马值与黏土矿物的含量有关，岩石中黏土矿物的含量越低岩石的放射性就越低，所以砾岩的放射性要低于砂岩和泥岩，因而徐家围子断陷沉积层中砾岩的生热率要略低于砂岩，泥岩的生热率高于砂岩。而煤是由生物遗体堆积经物理化学作用转化而成，放射性物质含量较低，生热率的值也较低，故煤层的生热率要明显低于泥岩。放射性物质在火山岩中含量较多，和 在酸性岩石的含量是其他岩性岩石的百倍，同时在酸性盐和中性岩中的含量也不低于其他岩石。徐家围子断陷中含有的火山岩有火山熔岩类、火山碎屑熔岩类、正常火山碎屑岩类、火山沉积碎屑岩类类岩石。火山熔岩类从基性到酸性分别为玄武岩、安山岩、英安岩、流纹岩、
14、粗面岩；火山碎屑熔岩分别为火山角砾岩、火山角砾凝灰岩、凝灰熔岩；正常火山碎屑岩分别为火山集块岩、火山角砾岩、凝灰岩、熔结凝灰岩；火山沉积碎屑岩分别为火山角砾岩、沉凝灰岩、凝灰质砂砾（或者砂、粉砂、泥）岩。所以徐家围子断陷火山岩中流纹岩的生热率最高，安山岩和火山角砾岩的生热率也较高，且远远高于砾岩、砂岩和泥岩的生热率。生热率与热流热流（）可 分 为 地 幔 热 流（）和 地 壳 热 流（），地壳热流（）主要来自沉积层放射性物质（、）衰变所产生的热量；地幔热流（）则来自软流圈传递给岩石圈的固体的热量。沉积层放射性物质的热流贡献主要由沉积层的厚度和沉积层放射性元素的丰度决定。第 卷毕然等：基于自然伽
15、马测井的徐家围子断陷沉积层生热率研究沉积层的热流（）计算公式为，（）其中：是沉积层第层的生热率的数值；是沉积层第层的厚度；是沉积层的层数。根据钻井资料显示，安达凹陷、升平凸起、徐西凹陷、徐东斜坡、徐南凹陷的沉积层平均厚度分别是 、，由式（）计 算 得 到 的 沉 积 层 热 流 贡 献 率 分 别 为、，大 地 热 流 值 分 别 为 、。沉积层放射性物质生热对这个地区的大 地 热流 的 贡 献 率 分 别 为、。沉积层放射性物质生热对地表热流的贡献率与地表热流呈负相关，表明沉积层放射性物质产生的热量对大地热流的影响并不是主要的。控制现今大地热流的主要因素是基地起伏，隆起区的大地热流值较高，但
16、在凹陷地区生热率较高，沉积厚度大，则沉积层放射性物质生热对大地热流的热流贡献将占较高的比重。大地热流值与基底起伏呈正相关，与沉积层放射性物质的生热对大地热流的贡献呈负相关，如图所示。由图可知，大地热流变化曲线图与基底起伏变化趋势相同，隆起区的热流较高，凹陷区的热流较低，安达凹陷区的热流值明显低于升平凸起区的热流值；在高大地热流的隆起区，沉积层放射性生热对大地热流的贡献小，而在低大地热流值的凹陷图徐家围子断陷 剖面热流结构示意图 区，沉积层放射性物质对大地热流的贡献高。因此在研究凹陷地区的地温场特征时，沉积层放射性物质生热对地温场影响更加显著。不同地区沉积层放射性物质生热率对大地热流的贡献结果表明，徐家围子断陷沉积层放射性物质生热对大地热流的贡献平均值是，约占大地热流 的，这种增温意义直观表现在地温梯度上，使得徐家围子断陷现今地温梯度 中 由沉积层放射性物质生热引起。沉积层放射性物质生热的增温效应会明显增加深层（深度大于）和超深层（深度超过）温度，分别能提高 和。根据油气有机成因理论，烃源岩的热演化与时间呈线性关系，与温度呈指数关系，温度每升高，反应速率将提高倍。所以沉积层放射性物质生热