全文快速搜索:   高级搜索

  中国石油大学学报(自然科学版)  2017, Vol. 41 Issue (4): 54-64  DOI:10.3969/j.issn.1673-5005.2017.04.007
0

引用本文 [复制中英文]

陈海峰, 王凤启, 王民, 等. 基于变系数ΔlogR技术的烃源岩TOC精细评价[J]. 中国石油大学学报(自然科学版), 2017, 41(4): 54-64. DOI: 10.3969/j.issn.1673-5005.2017.04.007.
[复制中文]
CHEN Haifeng, WANG Fengqi, WANG Min, et al. TOC evaluation using variable-coefficient ΔlogR model[J]. Journal of China University of Petroleum (Edition of Natural Science), 2017, 41(4): 54-64. DOI: 10.3969/j.issn.1673-5005.2017.04.007.
[复制英文]

基金项目

国家自然科学基金项目(41402110);中央高校基本科研业务费专项(14CX05017A)

作者简介

陈海峰(1974-),男,副教授,博士,研究方向为地球物理测井与油气资源评价。E-mail:chf_bj@126.com

文章历史

收稿日期:2016-04-26
基于变系数ΔlogR技术的烃源岩TOC精细评价
陈海峰1,2 , 王凤启1,2 , 王民3 , 于惠宇4     
1. 东北石油大学地球科学学院, 黑龙江大庆 163318;
2. 东北石油大学非常规油气成藏与开发省部共建重点实验室, 黑龙江大庆 163318;
3. 中国石油大学非常规油气与新能源研究院, 山东青岛 266580;
4. 中国石油勘探开发研究院西北分院, 甘肃兰州 730020
摘要: 海拉尔盆地乌尔逊凹陷主力烃源岩层(南一段)烃源岩非均质性强, 总有机碳(TOC)波动明显, 利用有限的测试样品表征整套烃源岩层的生烃潜力具有较强的不确定性; 烃源岩层的测井响应表现为高自然伽马、高声波时差、中等电阻率, 易于识别, 但是烃源岩测井评价参数的通用性差, 难以建立起具有普遍适用性的TOC预测公式。针对上述问题, 采用变系数ΔlogR技术评价该区烃源岩TOC, 将ΔlogR技术中的经验参数视为待定系数, 从模型参数地质意义出发, 通过针对单井和“分区、分相”预测模型中的待定参数, 完成全区120口井TOC测井评价, 利用测井计算TOC实现烃源岩分级评价, 通过“钻井定频率、层序定展布”的方法, 得到南一段不同丰度级别烃源岩的厚度。结果表明, 变系数ΔlogR技术预测TOC误差平均为16.6%, 比传统方法平均降低19.9%;南一段差烃源岩厚度一般15~20 m, 中等烃源岩厚度一般10~20 m, 优质烃源岩的厚度一般20~45 m; 南一段差烃源岩、中等烃源岩、优质烃源岩的体积分别为0.46×1011、1.08×1011和1.59×1011 m3
关键词: 测井    烃源岩    总有机碳    变系数ΔlogR技术    乌尔逊凹陷    海拉尔盆地    
TOC evaluation using variable-coefficient ΔlogR model
CHEN Haifeng1,2 , WANG Fengqi1,2 , WANG Min3 , YU Huiyu4     
1. Geoscience College, Northeast PetroleumUniversity, Daqing 163318, China;
2. Province Key Laboratory of Unconventional Oil and Gas Reservoir and Development, Northeast Petroleum University, Daqing 163318, China;
3. Research Institute of Unconventional Petroleum and Renewable Energy in China University of Petroleum, Qingdao 266580, China;
4. Research Institute of Petroleum Exploration and Development China Northwest Branch, Lanzhou 730020, China
Abstract: The source rocks in the First Member of Nantun Formation of Wuerxun sag, Hailaer Basin are characterized by strong organic heterogeneity with large variations in total organic carbon(TOC) values, resulting in the organic abundance of source rocks which cannot be represented by limited and discrete test samples. Source-rock sections can be easily recognized by high natural gamma ray, high interval transit time and medium resistivity logging response, but it is difficult to establish a universal TOC prediction formula, because the key parameters in ΔlogR model vary from well to well. Unlike in conventional approaches where ΔlogR is approximated empirically, an new approach that assumes ΔlogR as unknow parameters (k and a) is applied. In the variable-coefficient ΔlogR model, "k" is predicted according to single well, whereas "a" is predicted in the light of zones and facies. This variable-coefficient ΔlogR model is applied to 120 wells. Sources rocks in the First Member of Nantun Formation are divided into three grades according to the relationship between measured TOC and Hydrogen index (IH). TOC limit for poor source rocks, medium source rocks and excellent source rocks are 1.0% and 2.0%, thus logging-derived TOC is used to classify source rocks. Within the sequence stratigraphic framework, the depth of source rocks of various grades is determined by calculating thickness ratio of various-grade source rocks in system tracts of single well and by using the Kriging interpolation between all the wells. The results show that logging calculated TOC by the variable-coefficient ΔlogR model agrees well with the measured values. The TOC predicting error by the variable-coefficient ΔlogR model is about 16.6%, with 19.9% reduction than the traditional ΔlogR model. For the three rock types, poor source rocks, medium source rocks, and excellent source rocks, their depths are 15-20 m, 10-20 m, and 20-45 m; and their volumes are 0.46×1011m3, 0.46×1011m3, and 1.08×1011m3, respectively. These three types account for 14.7%, 34.5%, and 50.8% of total source rocks in volume.
Keywords: well logging    source rocks    total organic carbon(TOC)    variable-coefficient ΔlogR model    Wuerxun sag    Hailaer Basin    

乌尔逊凹陷是海拉尔盆地内油气资源较为富集的地区之一[1-2], 中国第一次资源评价预测油气资源量为5.56×108 t, 第二次资源评价预测油气资源量为1.89×108 t, 第三次资源评价预测油气资源量为3.39×108 t, 三轮油气资源评价结果出入较大, 同时, 目前该区的油气勘探程度较低, 近年来鲜有工业油井被发现, 有必要进一步论证油气资源潜力。烃源岩是油气生成的物质基础[3-5], 烃源岩精细评价是客观评价油气资源潜力的重要前提。在乌尔逊凹陷, 以往的油气资源评价均利用岩心测试地球化学指标(主要是总有机碳含量, TOC)的平均值代表整套烃源岩层的生烃能力[6], 然而该区烃源岩的非均质性强、测试样品有限且离散、不同有机质含量烃源岩的生、排烃能力存在差别[7-8], 导致“平均法”预测的油气资源量具有较大的不确定性。换言之, 研究区新一轮的资源评价须建立在丰富的评价数据和烃源岩分级评价的基础上, 即须开展烃源岩精细评价。在烃源岩精细评价方法中, 最高效的是测井评价技术[9-10]。高分辨率的测井资料不仅可提供丰富的TOC数据[11-13], 评价结果能够指导烃源岩分级评价[14], 结合层序地层分布还能够预测不同丰度级别烃源岩的厚度分布和发育体积[15-16]。烃源岩TOC测井评价技术已经比较成熟, ΔlogR技术、多元回归、神经网络等方法都有较多成功应用的实例[14-21]。其中, ΔlogR技术所需资料容易获取且可操作性较强, 在烃源岩TOC评价中应用最为广泛[17-22]。变系数ΔlogR技术[23-24]将ΔlogR技术内对TOC敏感的经验参数当成待定系数, 依据研究区实际资料厘定待定系数的值, 改善了ΔlogR技术对陆相高泥质、低TOC含量烃源岩TOC的评价效果。笔者以乌尔逊凹陷主力烃源岩层(南一段)为研究对象, 以岩心测试TOC和常规测井曲线为基础, 针对该区烃源岩非均质性强、富含泥质、TOC含量偏低、烃源岩测井评价参数的通用性差的特点, 采用变系数ΔlogR技术评价烃源岩TOC。从模型参数的地质意义出发, 通过基于单井和“分区、分相”预测待定系数的方法, 实现烃源岩TOC的准确预测。

1 基本地质特征

海拉尔盆地位于内蒙古自治区呼伦贝尔盟西南部, 是松辽盆地最重要的接替领域。乌尔逊凹陷南邻贝尔凹陷, 北与红旗、新宝力格凹陷局部连通, 呈南北向展布, 为一个西断东超的箕状断陷。盆地基底为前古生界和古生界, 沉积盖层自下而上为下白垩统的塔木兰沟组、铜钵庙组、南屯组、大磨拐河组、伊敏组, 上白垩统的青元岗组和新生界的第三、第四系。富含有机质的暗色泥岩主要发育在南屯组、大磨拐河组和伊敏组, 其中, 大磨拐河组和伊敏组由于绝大部分暗色泥岩尚未生烃而无法成为有效的烃源岩层[12]。南屯组自下而上又可分为南一段和南二段:南一段时期曾发生大规模湖侵, 深湖—半深湖暗色泥岩大面积发育, 是重要的烃源岩层; 南二段时期水体变浅, 辫状河三角洲砂泥互层发育, 油气储集层较发育。乌尔逊凹陷目前钻井169口, 钻穿南一段井120口, 研究区油流井72口。

2 目标层烃源岩特征

本次将研究区烃源岩精细评价的对象确定为南一段, 原因有二:一是乌尔逊凹陷中机质丰度高、类型好、生排烃能力强的优质烃源岩主要发育在南一段[6], 油源对比表明该区的原油主要来自南一段[1], 生烃模拟显示南一段烃源岩的生、排烃量均超过全区源岩生、排烃总量的90%[15], 即南一段为起控藏作用的主力烃源岩层; 二是南一段钻井中有35口井进行了较为系统的烃源岩取样分析, 开展了高分辨率层序地层研究[26-27], 具备开展烃源岩精细评价的地质基础, 但此项工作还未开展。

2.1 有机地球化学特征

南一段地层埋深1 250~3 250 m, 地层岩性主要为暗色泥岩夹薄层砂岩。暗色泥岩的沉积环境包括湖相和辫状河三角洲相, 颜色以灰黑色和灰色为主, 厚度一般为50~100 m, 最厚超过120 m。暗色泥岩实测TOC值为0.2%~9.8%, 平均值为2.41%, 局部TOC平均为7.68%, 氯仿沥青“A”平均为0.498%, 达到优质烃源岩级别[7]。暗色泥岩有机质主要源于低等水生生物, 典型特征是饱和烃气相色谱以 nC19 为主峰, 饱和烃 δ13C 分布在-29.37‰~-33.58‰[7]。暗色泥岩镜质体反射率 Ro 值为1.0%~1.4%, 处于生烃高峰阶段。暗色泥岩TOC波动频繁, 如W23井(图 1), 该井岩心实测TOC分布在0.2%~5.5%, 其中2 750~2 100 m的烃源岩TOC数值可达到其他深度段的2倍以上, 若利用岩心实测TOC的平均值代表整段烃源岩TOC含量, 则夸大了低TOC烃源岩的发育比例, 掩盖了高TOC烃源岩, 显然不利于客观评价油气资源潜力。

图 1 乌尔逊凹陷南一段烃源岩有机非均质性及测井响应(乌23井) Fig.1 Organic heterogeneity and logging response of source rocks in the First Member of Nantun Formation of Wuerxun sag (well Wu23)
2.2 测井响应特征

烃源岩层识别是测井评价烃源岩TOC的基础, 南一段烃源岩层具有高自然伽马、高声波时差、中等阻率测井响应, 与砂岩层的测井响应区别明显, 烃源岩层易于识别(图 1)。从图 1还可以看出, 烃源岩层的TOC数值与声波时差和电阻率叠合后的“幅度差”具有较好的相关性:TOC较低时, 如2 050~2 070 m, “幅度差”也较低; 随着烃源岩层TOC含量增高, 如2 070~2 100 m, “幅度差”也呈增大趋势。显然, 研究区烃源岩层在测井曲线上不仅容易识别, 而且TOC与“幅度差”呈正相关趋势, 符合利用ΔlogR技术评价TOC的前提。

然而, ΔlogR技术在研究区的应用效果并不理想, 即对大多数井烃源岩TOC预测误差较大。以W23井为例, 经过反复验证, ΔlogR技术[14]中唯一的待定参数Ro取值1.2%, 由图 1中TOC计算结果可知, ΔlogR技术预测TOC与岩心实测TOC除了在个别点比较接近以外, 多数测点两者明显分离, 测井预测TOC与实测TOC的平均相对误差达到38.8%(图 1)。分析认为, ΔlogR技术预测TOC误差大的原因在于其模型参数的通用(针对性)差:①根本原因是研究区烃源岩非均质性强、泥质含量高、TOC含量低, 导致单井之间烃源岩的测井响应并不完全一致, 因此对应的TOC预测公式也不完全相同, 进而采用同一预测公式必然导致TOC预测误差; ②直接原因是ΔlogR技术中模型参数的针对性差, 已有研究表明[18-20], 在模型框架不变的情况下, 通过调整模型中解释参数, TOC预测效果可得到明显改善。

综上分析可知, 乌尔逊凹陷烃源岩层的测井响应明显, 测井特征符合利用ΔlogR技术预测烃源岩层预测TOC, 但是较强的地层非均质性使得ΔlogR技术的应用效果不佳, 依据研究区实测资料修正ΔlogR模型中的解释参数是提高TOC预测精度可行的方法。

3 变系数ΔlogR技术在研究区的应用

采用变系数ΔlogR技术评价南一段烃源岩TOC, 重点探讨模型中两个关键待定系数的确定方法, 通过针对性地确定这两个待定系数, 实现了烃源岩TOC的准确预测。

3.1 模型的方法原理

变系数ΔlogR技术预测烃源岩TOC的公式为

$ \Delta {\rm{log}}R=k({\rm{log}}R-B)+\left( 1-k \right)(\Delta t-B). $ (1)
$ w({\rm{TOC}})=a\Delta {\rm{log}}R+b. $ (2)

式中, R为电阻率, Ω·m; Δt为声波时差, μs/m; B为基线值, 基线位置对应贫有机质的细粒岩处, 具体按照文献[17]中方法确定; w(TOC)为总有机碳; b为有机碳含量的背景值, 根据实测TOC统计得到研究区b取值0.2;ka为两个关键的待定系数, 是重点探讨的对象。

3.2 待定系数的确定 3.2.1 待定系数的确定原则

从待定系数的地质意义出发, 建立待定系数ka的厘定方法。为了便于论述, 选取研究区测井曲线质量好且进行TOC系统取样分析的井, 利用网格搜索法[19]分别得到单井待定系数ka与TOC预测误差的关系曲线(图 2图 3), 记TOC误差最小值对应的ka分别为“最优k值”和“最优a值”。

图 2 待定系数k与TOC预测误差的关系 Fig.2 Relationship beteween TOC prediction errors and parameters k
图 3 待定系数a与TOC预测误差的关系 Fig.3 Relationship between TOC prediction errors and parameter a

待定系数k能够改变声波时差和电阻率在ΔlogR中的相对比重, 由于声波时差主要对有机质中的干酪根响应, 电阻率主要对烃类响应, 因此k的地质意义是调整模型识别干酪根和烃的能力[20]。此外, 待定系数k值还具有削弱干扰因素, 突出TOC测井响应的作用[23]。由于烃源岩生、排烃过程本身具有较强的非均质性, 加上干扰TOC测井响应的因素异常复杂(包括岩性、物性、流体矿化度等), 导致单井之间的“最优k值”不尽一致, 甚至出现较大差异。如图 2所示, Wd3井(图 2(a))和Wd1井(图 2(c))代表三角洲前缘相, 两者的岩性组合和TOC值均相当, 将Wd1井“最优k”值应用于Wd3井TOC预测, 则Wd3井TOC预测误差从18.6%增加至35.1%, 误差增加16.5%;同理, S33井(图 2(b))和S35井(图 2(d))代表湖相沉积, 两者岩性组合和TOC值也相当, 将S33井“最优k”值应用于S35井TOC预测, 则S35井TOC预测误差从16.3%增加至21.4%, 误差增加5.1%。此外, Wd1井(氯仿沥青A平均值0.028%)有机质中烃的含量明显比Wd3井(氯仿沥青A平均0.011%)、S33井(氯仿沥青A平均0.021%)和S35井(氯仿沥青A平均0.015%)高, 这很可能是该井待定系数k值偏高的原因。

总之, 受烃源岩非均质性的影响, 单井之间“最优k值”不尽一致, 更重要的是, TOC预测误差对k值比较敏感, “借用”其他井的k值预测TOC很容易导致较大的误差, 因此将待定系数k的厘定原则确定为依据单井计算。

受有机质来源、保存、演化等影响, 相同ΔlogR值代表的TOC值可能不同, 如TOC数值相同的未成熟烃源岩与成熟烃源岩, 后者通常因含烃而表现出更大的电阻率和ΔlogR值, 因此待定系数a的地质意义是反映ΔlogR代表TOC大小的能力[24]。同一层位的、邻井之间的烃源岩的沉积环境、有机质来源、保存、演化程度等通常相近, 进而相同ΔlogR值代表的TOC值相当。如图 3所示, 将均代表三角洲前缘的Wd1井(图 3(c))的“最优a值”应用于Wd3井(图 3(a))TOC预测, S33井TOC预测误差从18.9%增加至19.5%, 误差仅增加0.6%;将均代表湖相的S35井(图 3(d))的“最优a值”应用于S33(图 3(b))井TOC预测, S33井TOC预测误差从16.8%增加至20.1%, 误差仅增加3.3%。Wd1井和S35井烃源岩沉积环境不同, 两者的“最优a值”差别较大, 互用彼此的“最优a值”则会导致较大的TOC误差。可见, 尽管不同单井之间“最优a值”也存在波动, 但是利用同一沉积环境中邻井的k值预测TOC并不会引入较大的TOC误差, 因此将待定系数a的厘定原则确定为依据周边井确定。

3.2.2 待定系数k的预测

前人[24]通过理论分析和实例论证, 建立了基于单井“logRt”交会图预测待定系数k的方法。该方法统称所有干扰TOC测井响应的因素为“干扰因素”, 指出待定系数k的实质是“干扰因素”导致声波时差曲线(t)和对数电阻率曲线(logR)的变化比例, 即“logRt”交会图的斜率。待定系数k的求取步骤为:①提取出目标层中与烃源岩互层发育的粉砂岩、粉砂质泥岩层; ②针对提取出的粉砂岩、粉砂质泥岩层, 以对数电阻率(logR)为横坐标, 以声波时差(Δt)为纵坐标, 构建交会图; ③利用最小二乘法拟计算交会图的斜率, 计算得到k值, 记为“预测k值”。

图 4为研究区35口TOC系统取样井中“预测k值”与“最优k值”的对比结果, 可见“logRt”交会图法得到的“预测k值”与基于单井实测TOC优化计算的“最优k值”非常吻合, 两者相关度的平方(R2)达到0.92。同时, 从典型井k与TOC预测误差关系曲线看(图 2(a)~图 2(d)), “预测k值”与“最优k值”对应的TOC误差非常接近, 表明“预测k值”计算TOC效果可靠, 切实可行。根据35口井统计, “预测k值”导致TOC相对误差平均增加4.2%, 可信度较高。

图 4 单井预测k值与最优k值的对比关系 Fig.4 Comparition between predicted k and ideal k for 35 single wells
3.2.3 待定系数a的预测

图 3分析表明待定系数a在同一区带、同一沉积环境的邻间之间变化不大, 说明“分区、分相”预测a值是可行的。为此, 依据南一段的展布形态和烃源岩发育环境的差异, 将南一段分为南部洼陷区和北部洼陷区, 将烃源岩沉积环境分为湖相和辫状河三角洲相。表 1统计了不同区带、相带中具有实测TOC井的待定系数a的优化计算结果。

表 1 待定系数a“分区、分相”方法预测结果 Table 1 Predicting results of parameter a for different facies and areas

分析可知:①同一区带同一沉积环境之间a值分布稳定, 平均绝对偏差不超过0.32, 说明分区、分相预测待定系数a是可行的; ②同一区带不同相带之间a值相差较大, 说明分相带预测a值是必要的; ③不同区带同一相带之间a值也存在一定的差异, 说明分区带确定a值也是有实际意义的。实际应用过程中, 对应缺少实测TOC而无法计算“最优a值”的井, 待定系数a取对应区带、相带a的平均值。

3.3 模型应用效果检验

统计变系数ΔlogR技术和传统ΔlogR技术在35口取样井、450个取样点的应用效果, 表 2中列出了部分数据点的计算结果。其中, 变系数ΔlogR技术中的待定系数k利用“logRt”交会图法计算, 待定系数a按照“分区、分相”方法确定; ΔlogR技术中的唯一的待定参数Ro依据实测值确定, 缺少化验数据的井通过Ro-深度的关系的趋势线估算。

表 2 变系数ΔlogR技术与传统ΔlogR技术TOC预测结果 Table 2 TOC calculating results of variable-coefficient ΔlogR technology and traditional ΔlogR technology

分析表明, 变系数ΔlogR技术预测TOC与实测TOC值的平均相对误差为16.6%, 对TOC低值(小于2%)和高值(大于2%)部分预测结果都比较可靠; 而ΔlogR技术预测TOC与实测TOC值的平均相对误差多数超过20%, 平均达到36.5%。变系数ΔlogR技术比传统ΔlogR技术平均降低19.9%, 而且井剖面上看(图 5), 变系数ΔlogR技术预测TOC与岩心实测TOC基本吻合, 测井计算TOC曲线准确地反映了实测TOC的垂向变化。

图 5 W23井南一段烃源岩TOC测井评价结果 Fig.5 TOC predicting results of source rocks in the First Member of Nantun Formation of well Wu23
4 研究区烃源岩精细评价结果

在验证变系数ΔlogR技术可行性的基础上, 将其应用于钻穿南一段且全区覆盖的120口井烃源岩TOC评价。在建立烃源岩分级评价标准的基础上, 依据测井计算TOC实现烃源岩分级评价和体系域中各级别烃源岩发育频率的确定, 采用“钻井定频率、层序控展布”的方法, 落实不同丰度级别烃源岩的厚度分布和体积。

4.1 烃源岩分级评价标准

不同丰度级别烃源岩的生排烃能力不同, 对油气成藏的贡献也不同, 因此开展烃源岩分级评价对客观评价油气资源潜力是必要的[15]。TOC和氢指数(IH)都反映烃源岩的生烃能力, 两者的交会图在烃源岩分级评价中已被证明具有较强的可行性[11]。依据研究区南一段烃源岩实测TOC和氢指数(IH)的关系(图 6), 按照TOC将乌尔逊凹陷烃源岩分为三级:TOC含量介于0.4%~1.0%时, IH随着TOC的增大缓慢增加, 该部分烃源岩的生烃能力较差, 评价为差烃源岩; TOC含量介于1.0%~2.0%时, IH随着TOC的增大迅速增加, 烃源岩生烃能力较强, 评价为中等烃源岩; TOC含量大于2.0%时, IH随着TOC的增大缓慢或趋于稳定, 该部分烃源岩的生烃能力很强, 评价为优质烃源岩。杨涛涛等[12]根据物质平衡原理得到乌尔逊凹陷优质烃源岩TOC值下限为2.0%, 本次确定的优质烃源岩TOC下限值与其一致。

图 6 乌尔逊凹陷南一段烃源岩TOC与IH的关系 Fig.6 Relationship between TOC and IH of sources rocks in the First Member of Nantun Formation
4.2 不同级别烃源岩的厚度分布

研究区南一段由一套完整的三级层序构成, 发育位体系域(LST)、水进体系域(TST)和高位体系域(HST)[27]。在烃源岩分级评价的基础上, 根据测井TOC曲线计算得到单井各体系域内不同丰度界别烃源岩的厚度比例, 即“钻井定频率”。以W23井为例(图 5), 其TOC计算参数见表 1, 该井优质烃源岩在低位体系域中的厚度比例约为5%, 水进体系域中优质烃源岩比例约为75%, 高位体系域中优质烃源岩的厚度比例约为35%, 即优质烃源岩主要发育湖侵体系域和高位域; 中等级别烃源岩在各体系域均有发育, 低位域内中等源岩厚度比例约为50%, 水进体系域内中等源岩厚度比例约为20%, 高位域内中等烃源岩厚度比例约为30%;差烃源岩主要分布在低位域, 厚度比例约为50%, 水进体系域内差源岩厚度比例约为2%, 高位域内差烃源岩厚度比例约为5%。同理, 可得到其他119口井中不同体系域内不同级别烃源岩的发育频率。

在“钻井定频率”的基础上, 结合地震分层和井间插值方法实现无井区烃源岩分布预测, 通过分层序累加的方法获取不同丰度级别烃源岩的厚度分布, 即“层序控展布”。操作方法如下:①以体系域为基本单元, 在同一体系域内, 在已知井间通过克里金插值的方法, 获取无井区不同级别烃源岩的厚度比例, 由体系域厚度和烃源岩的发育频率的乘积得到该体系域中某一级别烃源岩的厚度; ②按此方法得到其他体系域中相应级别烃源岩的厚度; ③对所有体系域同一级别烃源岩厚度进行累加, 得到该级别烃源岩厚度分布。

图 7为南一段不同级别烃源岩的厚度分布, 其中, 差烃源岩具有分布面积大、厚度薄的特征, 分布面积与南一段面积相当, 厚度一般为15~20 m, 凹陷四周相对发育, 向湖盆中心厚度略微增大(图 7(a)); 中等级别烃源岩具有发育面积大、厚度中等的特征, 累计厚度大于10 m的中等烃源岩的分布面积约占南一段面积的80%, 厚度一般为10~25 m, 南一段北部和南部发育程度大致相当(图 7(b)); 优质烃源岩具有发育面积中等、厚度大的特征, 累计厚度超过25 m的优质烃源岩的分布面约占南一段面积的55%, 厚度一般为25~45 m, 局部最厚超过65 m, 南一段北部和南部各发育一个厚度高值区, 北部发育程度略高(图 7(c))。

图 7 乌尔逊凹陷南一段不同级别烃源岩厚度分布 Fig.7 Depth contour map of source rocks of various grades in the First Member of Nantun Formation of Wuerxun Sag
4.3 不同级别烃源岩的体积

对南一段进行1 km×1 km网格划分, 以网格为单元, 对图 7中不同丰度级别烃源岩的厚度和面积进行积分, 得到差烃源岩的体积为0.46×1011 m3, 中等烃源岩的体积为1.08×1011 m3, 优质烃源岩的体积为1.59×1011 m3(图 8)。进一步计算得到南一段烃源岩总体积为3.13×1011 m3, 差级别烃源岩的体积比例为14.7%, 中等级别烃源岩的体积比为34.5%, 优质烃源岩的体积比例为50.8%。

图 8 乌尔逊凹陷南一段不同级别烃源岩体积比例 Fig.8 Volume ratio of source rocks of different grades in the First Member of Nantun Formation
5 结论

(1) 研究区主力烃源岩层(南一段)烃源岩的非均质性强, ΔlogR技术采用经验公式(参数)预测该区TOC平均相对误差达到36.5%, 评价结果无法有效地指导勘探实践。

(2) 变系数ΔlogR技术将传统ΔlogR技术中的经验参数视为待定系数(ka), 通过利用“ΔlogRt”交会图预测待定系数k和“分区、分相”法预测待定系数a, 变系数ΔlogR技术预测TOC与实测TOC的平均相对误差为16.6%, 比传统方法平均降低19.9%, 实现了南一段烃源岩TOC准确评价。

(3) 依据实测TOC与氢指数(IH)的关系将南一段烃源岩分为3级, 应用过程中采用测井计算TOC指导烃源岩分级。以120口单井TOC测井评价结果和精细的地震分层为基础, 通过“钻井定频率、层序定展布”的方法, 得到南一段差烃源岩厚度一般为15~20 m, 中等烃源岩厚度一般为15~20 m, 优质烃源岩的厚度一般为20~45 m, 差烃源岩、中等烃源岩、优质烃源岩的体积分别为0.46×1011、1.08×1011和1.59×1011 m3, 体积比例分别为14.7%、34.5%和50.8%。

参考文献
[1]
吴海波, 李军辉, 刘赫. 乌尔逊-贝尔凹陷岩性-地层油藏形成条件及分布规律[J]. 中南大学学报(自然科学版), 2015, 46(6): 2178-2187.
WU Haibo, LI Junhui, LIU He. Formation condition and distribution rule of lithologic and stratigraphic reservoirs in Wurxun-Beir Sag[J]. Journal of Central South University (Science and Technology), 2015, 46(6): 2178-2187. DOI:10.11817/j.issn.1672-7207.2015.06.028
[2]
平贵东, 吕延防, 范立民, 等. 海拉尔盆地乌尔逊-贝尔凹陷油气富集规律及主控因素分析[J]. 中南大学学报(自然科学版), 2013, 44(10): 4167-4178.
PING Guidong, LV Yanfang, FAN Limin, et al. Rules and main controlling factors of hydrocarbon enrichment of Urxun-Beier Depression, Hailar Basin[J]. Journal of Central South University (Science and Technology), 2013, 44(10): 4167-4178.
[3]
吴红烛, 黄志龙, 柳波, 等. 马朗凹陷哈尔加乌组烃源岩分布和生烃演化与石炭系成藏关系[J]. 中国石油大学学报(自然科学版), 2014, 38(1): 22-31.
WU Hongzhu, HUANG Zhilong, LIU Bo, et al. Relationship between distribution, hydrocarbon generation and evolution of source rock in Haerjiawu formation and Carboniferous hydrocarbon accumulation in Malang sag, Santanghu Basin[J]. Journal of China University of Petroleum (Edition of Natural Science), 2014, 38(1): 22-31.
[4]
孙耀庭, 徐守余, 张世奇, 等. 沾化凹陷多元供烃成藏特征及成藏模式[J]. 中国石油大学学报(自然科学版), 2015, 39(6): 42-49.
SUN Yaoting, XU Shouyu, ZHANG Shiqi, et al. Reservoir characteristics and reservoir forming model of multi element service supply in Zhanhua sag[J]. Journal of China University of Petroleum (Edition of Natural Science), 2015, 39(6): 42-49.
[5]
程亮, 刘德华, 夏志刚, 等. 江汉盆地潜江凹陷新沟嘴组致密油成藏模式与分布规律[J]. 中国石油大学学报(自然科学版), 2015, 39(6): 34-41.
CHENG Liang, LIU Dehua, XIA Zhigang, et al. Accumulation model and distribution regularity of Xingouzui Formation tight oil in Qianjiang sag, Jianghan Basin[J]. Journal of China University of Petroleum (Edition of Natural Science), 2015, 39(6): 34-41.
[6]
刘新颖, 邓宏文, 邸永香, 等. 海拉尔盆地乌尔逊凹陷南屯组优质烃源岩发育特征[J]. 石油实验地质, 2009, 31(1): 68-74.
LIU Xinying, DENG Hongwen, DI Yongxiang, et al. Highquality source rocks of Nantun formation in Wuerxun Depression, the Hailaer Basin[J]. Petroleum Geology & Experiment, 2009, 31(1): 68-74. DOI:10.11781/sysydz200901068
[7]
张林晔, 孔祥星, 张春荣, 等. 济阳坳陷下第三系优质烃源岩的发育及其意义[J]. 地球化学, 2003, 32(1): 35-42.
ZHANG Linye, KONG Xiangxing, ZHANG Chunrong, et al. High quality oil prone source rocks in Jiyang depression[J]. Geochemica, 2003, 32(1): 35-42.
[8]
JARVIE D.Evaluation of hydrocarbon generation and storage in the Barnett shale, Ft.Worth Basin, Texas[R].Texas:Humble Geochemical Services Division, 2004.
[9]
PASSEY Q R, CREANEY S, KULLA J B. A practical model for organic richness from porosity and resistivity logs[J]. AAPG Bulletin, 1990, 74(5): 1777-1794.
[10]
KADKHODAIE-ILKHCHI A, REZAEE M R, RAHIMPOUR-BONAB H. A committee neural network for prediction of normalized oil content from well log data:an example from South Pars Gas Field, Persian Gulf[J]. Journal of Petroleum Science and Engineering, 2009, 65(1/2): 23-32.
[11]
金涛, 高日胜, 高彩霞, 等. 从烃源岩测井评价结果看冀中坳陷饶阳凹陷勘探前景[J]. 天然气地球科学, 2010, 21(3): 406-412.
JIN Tao, GAO Risheng, GAO Caixia, et al. Suggestion of exploration prospect in Raoyang Sag of Jizhong Depression from the logging evaluation results of source rock[J]. Natural Gas Geoscience, 2010, 21(3): 406-412.
[12]
杨涛涛, 范国章, 吕福亮, 等. 烃源岩测井响应特征及识别评价方法[J]. 天然气地球科学, 2013, 24(2): 414-422.
YANG Taotao, FAN Guozhang, LV Fuliang, et al. The logging features and identification methods of source rocks[J]. Natural Gas Geoscience, 2013, 24(2): 414-422.
[13]
郭泽清, 孙平, 刘卫红. 利用ΔlogR技术计算柴达木盆地三湖地区第四系有机碳[J]. 地球物理学进展, 2012, 27(2): 626-633.
GUO Zeqing, SUN Ping, LIU Weihong. The carbon calculation by ΔlogR technology in Sanhu area of Qridos basin[J]. Progress in Geophysics, 2012, 27(2): 626-633.
[14]
杨少春, 王娜, 李明瑞, 等. 鄂尔多斯盆地崇信地区三叠系延长组烃源岩测井评价[J]. 天然气地球科学, 2013, 24(3): 470-476.
YANG Shaochun, WANG Na, LI Mingrui, et al. The logging evaluation of source rocks of Triassic Yanchang Formation in Chongxin area, Ordos Basin[J]. Natural Gas Geoscience, 2013, 24(3): 470-476.
[15]
卢双舫, 马延伶, 曹瑞成, 等. 优质烃源岩评价标准及其应用:以海拉尔盆地乌尔逊凹陷为例[J]. 地球科学-中国地质大学学报, 2012, 37(3): 535-544.
LU Shuangfang, MA Yanling, CAO Ruicheng, et al. Evaluation criteria of high-quality source rocks and its applications:taking the Wuerxun Sag in Hailaer Basin as an example[J]. Earth Science—Journal of China University of Geocsiences, 2012, 37(3): 535-544.
[16]
王宗礼, 李君, 罗强, 等. 利用测井资料精细评价冀中坳陷廊固凹陷烃源岩[J]. 天然气地球科学, 2012, 23(3): 430-437.
WANG Zongli, LI Jun, LUO Qiang, et al. Fine evaluation of source rock with logging data in Langgu sag of Jizhong depression[J]. Natural Gas Geoscience, 2012, 23(3): 430-437.
[17]
张志伟, 张龙海. 测井评价烃源岩的方法及其应用效果[J]. 石油勘探与开发, 2000, 27(3): 84-87.
ZHANG Zhiwei, ZHANG Longhai. A method of source rock evaluation by well logging and its application result[J]. Petroleum Exploration and Development, 2000, 27(3): 84-87.
[18]
雷涛, 周文, 杨艺, 等. 塔里木盆地塔中地区海相烃源岩测井评价方法[J]. 岩性油气藏, 2010, 22(4): 89-94.
LEI Tao, ZHOU Wen, YANG Yi, et al. Logging evaluation method for marine source rocks in middle Tarim Basin[J]. Lithologic Reservoirs, 2010, 22(4): 89-94.
[19]
朱光有, 金强, 张林晔. 用测井信息获取烃源岩的地球化学参数研究[J]. 测井技术, 2003, 27(2): 104-109.
ZHU Guangyou, JIN Qiang, ZHANG Linye. Using log information to analyze the geochemical characteristics of source rocks in Jiyang Depression[J]. Well Logging Technology, 2003, 27(2): 104-109.
[20]
刘超. 测井资料评价烃源岩方法改进及应用[D]. 大庆: 东北石油大学, 2011.
LIU Chao.Evaluating source rock using logging data:improvement and application [D].Daqing: Northeast Petroleum University, 2011.
[21]
胡慧婷, 卢双舫, 刘超, 等. 测井资料计算源岩有机碳含量模型对比及分析[J]. 沉积学报, 2011, 29(6): 1199-1205.
HU Huiting, LU Shuangfang, LIU Chao, et al. Models for calculating organic carbon content from logging information:comparison and analysis[J]. Acta Sedimentologica Sinica, 2011, 29(6): 1199-1205.
[22]
刘超, 卢双舫, 黄文彪, 等. ΔlogR技术改进及其在烃源岩评价中的应用[J]. 大庆石油地质与开发, 2011, 30(3): 27-31.
LIU Chao, LU Shuangfang, HUANG Wenbiao, et al. ΔlogR technique:improvement and its application in source rocks evaluation[J]. Petroleum Geology & Oilfield Development in Daqing, 2011, 30(3): 27-31.
[23]
刘超, 卢双舫, 薛海涛. 变系数ΔlogR方法及其在泥页岩有机质评价中的应用[J]. 地球物理学进展, 2014, 29(1): 312-317.
LIU Chao, LU Shuangfang, XUE Haitao. Variable-coefficient ΔlogR model and its application in shale organic evaluation[J]. Progress in Geophysics, 2014, 29(1): 312-317. DOI:10.6038/pg20140144
[24]
刘超, 印长海, 卢双舫. 变系数ΔLgR烃源岩测井评价技术关键参数厘定方法及应用[J]. 天然气地球科学, 2015, 26(10): 1925-1931.
LIU Chao, YIN Changhai, LU Shuangfang. Predicting key parameters for variable-coefficient ΔLgR logging technique and its application in source rocks evaluation[J]. Natural Gas Geoscience, 2015, 26(10): 1925-1931. DOI:10.11764/j.issn.1672-1926.2015.10.1925
[25]
高军义, 刘志宏, 吴相梅, 等. 海拉尔盆地乌尔逊凹陷构造变形对沉降中心迁移的控制[J]. 吉林大学学报(地球科学版), 2014, 44(1): 15-24.
GAO Junyi, LIU Zhihong, WU Xiangmei, et al. Structural deformation control over the subsidence center, migration of Wuerxun sag in Hailar Basin[J]. Journal of Jilin University(Earth Science Edition), 2014, 44(1): 15-24.
[26]
柳成志, 李军辉. 海拉尔盆地乌尔逊凹陷南屯组层序地层及沉积体系分析[J]. 地层学杂志, 2013, 37(3): 303-312.
LIU Chengzhi, LI Junhui. Sequence stratigraphy and sedimentary system of analysis the nantun formation in Wuerxun depression in Hailaer Basin[J]. Journal of Stratigraphy, 2013, 37(3): 303-312.
[27]
吴海波, 李军辉, 刘赫. 海拉尔盆地乌尔逊-贝尔凹陷层序构成样式及油气成藏模式[J]. 岩性油气藏, 2015, 27(5): 155-160.
WU Haibo, LI Junhui, LIU He. Sequence architecture pattern and hydrocarbon accumulation model of Lower Cretaceous in Wuerxun-Beier Depression, Hailaer Basin[J]. Lithologic Reservoirs, 2015, 27(5): 155-160.