2. 中国石油天然气集团有限公司天然气成藏与开发重点实验室, 河北廊坊 065007
2. Key Laboratory of Natural Gas Accumulation and Development, CNPC, Langfang 065007, China
吐哈盆地位于新疆东部, 是新疆三大含油气盆地之一。盆地经历了二叠纪-三叠纪断拗充填阶段和侏罗纪以来的前陆盆地演化阶段, 两大构造层分别发育以中二叠统桃东沟群湖相暗色泥岩和中-下侏罗统水西沟群煤系为主的烃源岩[1-2]。以往的勘探研究对象多侧重于中-下侏罗统水西沟群煤系烃源岩。
近年来, 随着勘探技术水平的提升, 吐哈探区勘探领域延伸至前侏罗系地层之中[3]。前侏罗系地层中目前已经发现多套烃源岩或潜在烃源岩地层, 与其它几套烃源岩地层相比, 中二叠统桃东沟群烃源岩埋深相对较浅、厚度大、在全区内稳定分布, 有望成为吐哈探区前侏罗系地层勘探领域的“主战场”。
尽管目前中二叠统桃东沟群泥质烃源岩虽然已经被证实为有效烃源岩[4], 但由于埋藏较深, 地震资料品质差, 钻井揭露少, 导致其勘探程度较低, 对烃源岩地化特征及生烃潜力等方面认识不足, 这将会影响到对盆地内前侏罗系油气资源前景的认识和勘探部署的决策。本文基于有机地球化学相关测试结果和测井资料, 对比现有预测烃源岩的测井解释方法预测结果, 优选出适合评价吐哈盆地桃东沟群烃源岩的测井解释方法, 进而系统地评价吐哈盆地桃东沟群烃源岩。
1 烃源岩的测井评价方法烃源岩测井评价始于20世纪40年代, 经过多年的积累和发展, 已经形成了多种用于反演烃源岩有机质丰度、成熟度和生烃潜量等地化参数的成熟方法[5]。
1.1 有机质丰度有机质丰度是指沉积物中的有机质含量, 常用总有机碳含量来评价有机质丰度[6-8]。利用测井技术计算得到的烃源岩的有机质丰度, 一直以来都是国内外研究的重点, 其计算原理基于烃源岩具有的高伽马(GR)、高声波时差(DT)、高电阻率(RD)和低密度(DEN)的“三高一低”特征。应用较为广泛的预测模型和方法包括ΔlogR模型及其改进型、多元回归参数法和地震反演预测法。
1.1.1 ΔlogR模型(传统ΔlogR法)ΔlogR模型由PASSEY等[9]于1990年提出, 其原理是将可以反映岩石孔隙度的声波时差曲线和电阻率曲线反向叠合, 并以细颗粒的非烃源岩为基线, 将两条曲线的幅度差定义为ΔlogR[8], 即:
$ \Delta \log R=\log \left(\frac{R}{R_{\text {基 }}}\right)+0.02 \times\left(D_{\mathrm{T}}-D_{\mathrm{T}_{\text {基 }}}\right) $ | (1) |
式中: R, R基, DT, DT基分别为电阻率、电阻率曲线基线值、声波时差值、声波时差曲线基线值, 其中R基与DT基为定值。
ΔlogR与总有机碳含量(TOC)呈线性关系, 并且为有机质成熟度的函数, 因此TOC可以表示为:
$ T_{\mathrm{OC}}=\Delta \log R \times 10^{1.5374-0.994 R_{\mathrm{o}}} $ | (2) |
式中: TOC为总有机碳含量, Ro为镜质体反射率[9]。采用该模型在鄂尔多斯盆地北部的杭锦旗地区评价有机质丰度, 取得了良好的应用效果[10]。
1.1.2 ΔlogR模型改进型(优化和广义ΔlogR法)由于ΔlogR模型在陆相及中、深层烃源岩的TOC计算结果误差较大, 所以国内学者对该模型进行了大量的修改, 其中优化ΔlogR模型和广义ΔlogR模型的应用范围较广泛。
优化ΔlogR模型由边雷博等[11]于2018年依据一一映射原理提出, 该模型将不同范围的声波时差及其对应的电阻率与260~460μs/m声波时差及其对应的电阻率一一映射, 从而建立起未知叠合系数与声波时差为260~460μs/m时的最佳叠合系数的映射关系, 因此ΔlogR可表示为:
$ \Delta \log R=\frac{\log \frac{R}{R_{\text {基 }}}}{1.05}+1.16 \frac{A_{\mathrm{C}}-A_{\mathrm{C}_{\mathrm{H}}}}{A_{\mathrm{C}_{\max }}-A_{\mathrm{C}_{\min }}} $ | (3) |
式中: ACmax=460μs/m, ACmin=260μs/m, AC和AC基为定值。
将利用(3)式得到的ΔlogR值代入(2)式, 即可获得有机碳含量的值, 该模型在渤海湾盆地东濮凹陷西南部洼陷沙河街组沙三段中亚段烃源岩的有机碳预测中, 获得了良好的应用效果。
广义ΔlogR模型由胡慧婷等[12]于2016年提出, 为了消除误差, 用GR曲线代替成熟度参数。王祥等[13]于2020年对该模型进行修改, 将GR和DEN一起使用, 因为GR值远大于DEN值, 故直接建立的模型公式主要受GR影响, 因此, 首先对GR取自然对数值, 再建立模型, 即:
$ \begin{gathered} T_{\mathrm{OC}}=\left(a \times \log G_{\mathrm{R}}+b \times D_{\mathrm{EN}}+c\right) \times \\ \Delta \log R+\Delta T_{\mathrm{OC}} \end{gathered} $ | (4) |
式中: a, b, c为拟合系数。利用该模型精确地预测了松辽盆地徐家围子断陷安达地区和渤中凹陷西部的烃源岩TOC值。
1.1.3 多元回归参数法该方法选取与TOC相关性较强的测井参数(如GR、DEN、DT等), 以实测TOC值为因变量, 测井参数为自变量进行一元、二元、三元等线性回归。以往的研究表明, 选择的参数越多, 其相关性就越高, 即:
$ \begin{gathered} T_{\mathrm{OC}}=a \times G_{\mathrm{R}}+b \times D_{\mathrm{T}}+c \times D_{\mathrm{EN}}+ \\ d \times S_{\mathrm{P}}+e \times R_{\mathrm{D}}+f \end{gathered} $ | (5) |
式中: a, b, c, d, e, f均为拟合系数; SP为自然电位值。该模型精确评价了太行山东麓的页岩气储层的地化特征[14]。
1.1.4 地震叠前反演法该方法通过三维地震资料获取了纵横波速度、纵波阻抗和横波阻抗等弹性数据, 经过岩石物理交会分析研究, 发现纵波速度和DEN与TOC的相关性最高, 因此TOC可表示为:
$ T_{\mathrm{OC}}=k+m \times v_{\mathrm{P}}+n \times D_{\mathrm{EN}} $ | (6) |
式中: k, m, n均为拟合系数; vP为纵波波速; DEN为密度。利用该模型计算渤中凹陷烃源岩TOC值取得了较好的效果[15]。
1.2 有机质成熟度常用于确定有机质成熟度的指标包括镜质体反射率(Ro)和岩石最高热解峰温(Tmax), 在陆相烃源岩中, 镜质体反射率的使用效果优于岩石最高热解峰温。利用现有测井参数计算Ro的方法主要为多元回归参数法, 以往的研究表明, Ro与深度存在经验公式:
$ R_{\rm o}=a \times H+b $ | (7) |
式中: H为埋深; a, b为拟合系数[16]。该公式在不同的地区的应用效果差异很大, 后续学者对镜质体反射率、常规测井参数和深度进行相关性分析后发现[17], 镜质反射率与深度、声波时差的相关性较好, 因此, 镜质体反射率可表示为:
$ R_{\mathrm{o}}=a \times H+b \times D_{\mathrm{T}}+c $ | (8) |
式中: H为埋深; DT为声波时差。
1.3 生烃潜量生烃潜量是评价烃源岩生烃能力的一个重要的参数, 主要通过热解实验获取[18]。大量的实验数据表明, 生烃潜量随着有机碳含量的增大而增大, 因此, 可以通过两者的统计关系直接获得, 以往研究发现, 生烃潜量PG与TOC值存在着如下关系:
$ P_{\mathrm{G}}=a \times \mathrm{e}^{b \times T _\mathrm{OC}} $ | (9) |
该模型的精度实际应用效果一般, 薛冰等[14]将实测生烃潜量与常规测井参数、TOC、Ro值进行相关性分析, 发现生烃潜量与Ro、DT和SP的相关性高, 得出:
$ P_{\mathrm{G}}=a \times R_{\mathrm{o}}+b \times D_{\mathrm{T}}+c \times S_{\mathrm{P}}+d $ | (10) |
研究区中二叠统烃源岩勘探程度较低, 钻遇中二叠统的钻井有15口(图 1), 其中常规测井数据及岩心录井资料齐全的钻井共有7口, 且这7口井中二叠统地层的沉积相均为半深湖—滨浅湖相。本次研究基于岩石热解等相关地球化学试验, 获得包括62个TOC数据点以及43个镜质体反射率和生烃潜量的数据点。
本文根据上述7口钻井的实测数据, 采用数理分析软件, 拟合得到上文中所述模型在研究区内的参数, 再将预测得到的参数与实测数据进行对比, 选择适合评价本研究区的测井方法。
2.1 有机碳含量由于研究区内地震资料品质较差, 地震叠前反演预测工作难以进行。基于此, 本文采用ΔlogR模型及其改进型和多元回归参数法反演区内烃源岩的TOC。
应用ΔlogR模型及其改进型, 将RD曲线和DT曲线进行叠加, 非烃源岩段中RD曲线和DT曲线相互平行或者相交, 相交处所对应的值为R基与DT基值, 利用对应的公式计算TOC值, 并与实测值进行对比。前人研究表明, 在使用多元回归参数法时, 选择的参数越多, 其相关性越高[19], 因此, 本次研究选择了声波时差、电阻率、自然电位、密度和伽马等测井参数进行拟合, 模型公式如表 1所示。
基于实测数据, 利用数理分析软件, 拟合出上述各方法的模型公式及相关系数(表 1)。在此基础上, 对比上述方法的应用结果, 发现多元回归参数法的应用效果远优于ΔlogR模型及其改进型(图 2)。因此, 本文采用多元回归参数法评价烃源岩TOC的值。
基于前人研究结果和本文1.2节的内容, 本次研究分析了研究区桃东沟群烃源岩Ro与深度、DT值的相关性, 结果如图 3a和图 3b所示。结果显示, Ro与深度的相关性较好, Ro与DT具有弱正相关关系, 这表明可以利用埋深和DT建立预测Ro模型。因此本次采用多元参数回归的方法, 分别建立了Ro与深度以及Ro与深度和DT值之间的关系式如下:
$ \left\{\begin{array}{l} R_{\mathrm{o}}=0.0001 \times H+0.4326 \\ R^{2}=0.7422 \end{array}\right. $ | (11) |
$ \left\{\begin{array}{l} R_{\mathrm{o}}=0.0001 \times H+0.0025 \times D_{T}+0.2446 \\ R^{2}=0.7766 \end{array}\right. $ | (12) |
因为(12)式的应用效果优于(11)式, 所以本文利用(12)式解释区内烃源岩有机质成熟度。
2.3 生烃潜量前人研究表明, 常用于计算生烃潜量的方法有两种, 一种是建立生烃潜量与TOC之间的直接关系, 另一种为多元回归参数法。因研究区勘探程度较低, 目的层钻遇井少, 地化参数的主要研究对象为露头剖面。本文对研究区露头与钻井实测的生烃潜量与TOC进行拟合(图 4), 发现生烃潜量与TOC之间存在如下关系:
$ \left\{\begin{array}{l} P_{\mathrm{G}}=2.0274 \times \mathrm{e}^{0.1374 \times T_{\mathrm{OC}}} \\ R^{2}=0.7393 \end{array}\right. $ | (13) |
前人对实测生烃潜量与常规测井参数、TOC、Ro值进行相关性分析, 发现生烃潜量与Ro、DT和SP的相关性高, 因此本文采用三元参数回归的方法, 得到生烃潜量与Ro、DT和SP的关系式如下:
$ \left\{\begin{aligned} P_{\mathrm{G}}=&-0.094 D_{\mathrm{T}}-0.1069 S_{P}+6.9486 R_{\mathrm{o}}+\\ & 14.1963 \\ R^{2}=& 0.3780 \end{aligned}\right. $ | (14) |
分别将(13)式与(14)式计算得到的生烃潜量预测值与实测值进行相关性分析(图 5), 其中, 公式(13)中TOC值为采用本文2.1节中的多元回归参数模型反演所得的数值。对比上述两种方法的相关性不难发现, 上述两种方法在研究区内的实际应用效果一般, 但基于热解实验数据所得到的模型优于利用回归参数建立的模型。因此本文采用(13)式预测区内烃源岩生烃潜量。
由吐哈盆地构造分区图可知, 目前钻遇中二叠统地层的钻井大都分布在台南凹陷和了墩隆起, 其中, 连22井和英17井的部分数据为收集得到。本次研究利用测井数据反演台南凹陷和了墩隆起的中二叠统烃源岩地化参数。
为了验证优选所得的烃源岩地化参数评价模型应用时的精确性, 本文利用钻井实测数据与不同方法得到的预测数据计算上述模型的平均误差率(表 2), 结果显示: 在研究区内, 利用多元回归参数法计算TOC的平均误差率为18.81%;利用二元线性回归法计算Ro的平均误差率为13.40%;利用(13)式计算PG的平均误差率为32.15%, 表明优选所得模型在评价烃源岩地化参数时具有较高的准确性。
基于测井方法预测结果, 分别绘制该区两口井中二叠统烃源岩综合地球化学剖面, 结果如图 6和图 7所示, 吐哈盆地中二叠统桃东沟群烃源岩地化参数频率直方图如图 8所示。
由图 6和图 7可以看出: 位于台南凹陷的鲁南1井地层变化较小, 烃源岩层主要为黑色泥岩, 为优质烃源岩; 位于了墩隆起的南湖1井所在地层变化明显, 主要烃源岩层段为多套深灰色泥岩, 为差-较好-好烃源岩。
由图 8可知, 吐哈盆地中二叠统桃东沟群烃源岩的TOC值多数大于1.0%, 镜质体反射率主要分布范围为0.5%~1.3%, PG主要分布范围为2.0~6.0mg/g。但是, 台南凹陷烃源岩TOC值和镜质体反射率要高于了墩隆起, 两者PG值相差不大, 台南凹陷略高。
根据我国陆相生油岩的有机质丰度和成熟度评价标准[21], 就TOC值而言, 吐哈盆地中二叠统桃东沟群烃源岩整体属于优质烃源岩; 就PG值而言, 吐哈盆地中二叠统桃东沟群烃源岩大多属于较优质烃源岩。因此, 吐哈盆地中二叠统桃东沟群烃源岩属于优质-较优质烃源岩, 成熟度为低成熟-成熟阶段, 主要为低成熟阶段。台南凹陷的烃源岩的优质烃源岩比例明显高于了墩隆起; 台南凹陷烃源岩处于低成熟-成熟阶段, 而了墩隆起的烃源岩处于低成熟阶段。
统计研究区内7口资料齐全钻井的TOC、镜质体反射率和生烃潜量的平均值(表 3), 结果显示: 台南凹陷钻井单井烃源岩TOC与PG值均高于了墩隆起的南湖1井, 成熟度差异不大, 均处于低成熟阶段。这表明台南凹陷的烃源岩的有机质丰度优于了墩隆起。
基于吐哈盆地桃东沟群烃源岩地球化学参数的实测数据, 对比目前应用较为广泛的烃源岩测井评价方法, 优选出多元回归参数法应用于总有机碳含量的计算, 将基于声波时差与埋深所建立的模型用于镜质体反射率的计算, 将基于有机碳含量与生烃潜量所建立的模型用于生烃潜量的评价。
将上述模型方法应用于吐哈盆地桃东沟群烃源岩地化参数评价, 结果表明, 研究区内桃东沟群烃源岩的总有机碳含量几乎均大于1.0%, 生烃潜量为2.0~6.0mg/g, 为优质-较优质烃源岩; 镜质体反射率主要分布范围为0.5%~1.3%, 处于低成熟-成熟阶段。
吐哈盆地桃东沟群烃源岩地化参数的纵向与横向对比结果表明: 台南凹陷烃源岩有机质丰度优于了墩隆起。
[1] |
李志军, 黄先律, 杨飞. 利用三维地震资料分析胜北构造带构造演化与油气聚集[J]. 石油物探, 1999, 38(2): 66-72. LI Z J, HUANG X L, YANG F. Analysis of structural evolution and hydrocarbon accumulation in Shengbei structural belt using 3D seismic data[J]. Geophysical Prospecting for Petroleum, 1999, 38(2): 66-72. |
[2] |
李文厚, 周立发, 柳益群, 等. 吐哈盆地沉积格局与沉积环境的演变[J]. 新疆石油地质, 1997, 18(2): 135-141. LI W H, ZHOU L F, LIU Y Q, et al. Sedimentary pattern and evolution of sedimentary environment in Turpan Hami Basin[J]. Xinjiang Petroleum Geology, 1997, 18(2): 135-141. |
[3] |
梁世君. 吐哈探区油气勘探成果及潜力[J]. 新疆石油地质, 2020, 41(6): 631-641. LIANG S J. Oil and gas exploration achievements and potential in Tuha exploration area[J]. Xinjiang Petroleum Geology, 2020, 41(6): 631-641. |
[4] |
高岗, 梁浩, 李华明, 等. 吐哈盆地石炭系—下二叠统烃源岩地球化学特征[J]. 石油勘探与开发, 2009, 36(5): 583-592. GAO G, LIANG H, LI H M, et al. Geochemical characteristics of Carboniferous Lower Permian source rocks in Tuha basin[J]. Petroleum Exploration and Development, 2009, 36(5): 583-592. DOI:10.3321/j.issn:1000-0747.2009.05.006 |
[5] |
袁东山, 王国斌, 汤泽宁, 等. 测井资料评价烃源岩方法及其进展[J]. 石油天然气学报, 2009, 31(4): 192-194. YUAN D S, WANG G B, TANG Z N, et al. Methods and progress of evaluating source rocks with logging data[J]. Journal of Oil and Gas Technology, 2009, 31(4): 192-194. DOI:10.3969/j.issn.1000-9752.2009.04.046 |
[6] |
黄世伟, 张廷山, 王顺玉, 等. 赤水地区上三叠统须家河组源岩特征及天然气成因[J]. 地球科学与环境学报, 2005, 27(3): 19-22. HUANG S W, ZHANG T S, WANG S Y, et al. Source rock characteristics and gas genesis of Upper Triassic Xujiahe Formation in Chishui area[J]. Journal of Earth Sciences and Environment, 2005, 27(3): 19-22. DOI:10.3969/j.issn.1672-6561.2005.03.005 |
[7] |
朱建伟, 赵刚, 刘博, 等. 油页岩测井识别技术及应用[J]. 吉林大学学报(地球科学版), 2012, 42(2): 289-295. ZHU J W, ZHAO G, LIU B, et al. The logging identification technology and application of oil shale[J]. Journal of Jilin University(Earth Science Edition), 2012, 42(2): 289-295. |
[8] |
王攀, 彭苏萍, 杜文凤, 等. 基于测井参数的煤系烃源岩总有机碳含量预测模型[J]. 煤炭学报, 2017, 42(5): 1266-1276. WANG P, PENG S P, DU W F, et al. Prediction model of total organic carbon content in coal measure source rocks based on logging parameters[J]. Journal of China Coal Society, 2017, 42(5): 1266-1276. |
[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(12): 1777-1794. |
[10] |
赵桂萍, 李良. 杭锦旗地区基于测井响应特征的泥质烃源岩有机质丰度评价研究[J]. 石油物探, 2016, 55(6): 879-886. ZHAO G P, LI L. Study on the evaluation of the organic matter abundance of argillaceous source rocks based on logging response characteristics in Hangjinqi area[J]. Geophysical Prospecting for Petroleum, 2016, 55(6): 879-886. DOI:10.3969/j.issn.1000-1441.2016.06.013 |
[11] |
边雷博, 柳广弟, 孙明亮, 等. 优化的ΔlogR技术及其在中—深层烃源岩总有机碳含量预测中的应用[J]. 油气地质与采收率, 2018, 25(4): 40-45. BIAN L B, LIU G D, SUN M L, et al. Optimized ΔlogR technique and its application to prediction of total organic carbon content in middle deep source rocks[J]. Petroleum Geology and Recovery Efficiency, 2018, 25(4): 40-45. |
[12] |
胡慧婷, 苏瑞, 刘超, 等. 广义ΔlogR技术预测陆相深层烃源岩有机碳含量方法及其应用[J]. 天然气地球科学, 2016, 27(1): 149-155. HU H T, SU R, LIU C, et al. The method and application of generalized ΔlogR technology in predicting organic carbon content in continental deep source rocks and its application[J]. Natural Gas Geoscience, 2016, 27(1): 149-155. |
[13] |
王祥, 马劲风, 张新涛, 等. 一种考虑密度因素的广义ΔlogR法预测总有机碳含量——以渤中凹陷西南部陆相深层烃源岩为例[J]. 地球物理学进展, 2020, 35(4): 1471-1480. WANG X, MA J F, ZHANG X T, et al. A generalized ΔlogR method considering density factor to predict total organic carbon content: A case study of continental deep source rocks in southwest Bozhong sag[J]. Progress in Geophysics, 2020, 35(4): 1471-1480. |
[14] |
薛冰, 刘度, 罗洪浩, 等. 太行山东麓页岩储层地球化学参数测井计算方法研究[J]. 煤炭科学技术, 2020, 48(7): 325-333. XUE B, LIU D, LUO H H, et al. Study on logging calculation method of geochemical parameters of shale reservoir in eastern foot of Taihang[J]. Coal Science and Technology, 2020, 48(7): 325-333. |
[15] |
马昭军, 唐建明, 徐天吉. 多波多分量地震勘探技术研究进展[J]. 勘探地球物理进展, 2010, 33(4): 247-253. MA Z J, TANG J M, XU T J. Research progress of multiwave and multicomponent seismic exploration technology[J]. Progress in Exploration Geophysics, 2010, 33(4): 247-253. |
[16] |
卢庆治, 胡圣标, 郭彤楼, 等. 川东北地区异常高压形成的地温场背景[J]. 地球物理学报, 2005, 48(5): 1110-1116. LU Q Z, HU S B, GUO T L, et al. Background of geothermal field formed by abnormal high pressure in northeastern Sichuan[J]. Chinese Journal of Geophysics, 2005, 48(5): 1110-1116. |
[17] |
信石印. 巴布亚盆地主力烃源岩特征及有效性分析[J]. 石油实验地质, 2018, 40(6): 843-848. XIN S Y. Characteristics and effectiveness analysis of main source rocks in Papua Basin[J]. Petroleum Geology and Experiment, 2018, 40(6): 843-848. |
[18] |
李贤庆, 侯读杰, 胡国艺, 等. 鄂尔多斯盆地中部地区下古生界碳酸盐岩生烃潜力探讨[J]. 矿物岩石地球化学通报, 2002, 21(3): 152-157. LI X Q, HOU D J, HU G Y, et al. Discussion on hydrocarbon generation potential of lower Paleozoic Carbonate rocks in the central area of Ordos Basin[J]. Bulletin of Mineralogy, Petrology and Geochemistry, 2002, 21(3): 152-157. |
[19] |
杜江民, 张小莉, 钟高润, 等. 致密油烃源岩有机碳含量测井评价方法优选及应用——以鄂尔多斯盆地延长组长7段烃源岩为例[J]. 地球物理学进展, 2016, 31(6): 2526-2533. DU J M, ZHANG X L, ZHONG G R, et al. Optimization and application of logging evaluation method for organic carbon content of tight oil source rocks: A case study of Chang 7 member of Yanchang Formation in Ordos Basin[J]. Progress in Geophysics, 2016, 31(6): 2526-2533. |
[20] |
刘春艳. 新疆三塘湖盆地石炭系烃源岩地球化学特征研究[D]. 兰州: 兰州大学, 2009 LIU C Y. Geochemical characteristics of Carboniferous source rocks in Santanghu basin, Xinjiang[D]. Lanzhou: Lanzhou University, 2009 |
[21] |
尚慧芸, 李晋超. 陆相生油岩有机质的丰度及类型[J]. 石油学报, 1981, 2(增刊1): 8-16. SHANG H Y, LI J C. The abundance and types of organic matter in terrestrial sourcerocks[J]. Acta Petrolei Sinica, 1981, 2(S1): 8-16. |