Origin analysis of paleo-overpressure of source rock based on logging response characteristics: Case study of the source rocks in the 7th member of Yanchang Formation of Triassic, Ordos Basin, China

  • Zeyang XU , 1, 2 ,
  • Zhen LIU , 1, 2 ,
  • Jingzhou ZHAO 3, 4 ,
  • Jiacheng DANG 5 ,
  • Jun LI 3, 4 ,
  • Ziyi TANG 3, 4
Expand
  • 1. State Key Laboratory of Petroleum Resources and Prospecting,China University of Petroleum(Beijing),Beijing 102249,China
  • 2. College of Geosciences,China University of Petroleum (Beijing),Beijing 102249,China
  • 3. Provincial Key Lab of Hydrocarbon Accumulation Geology,Xi’an Shiyou University,Xi’an 710065,China
  • 4. College of Earth Sciences & Engineering,Xi’an Shiyou University,Xi’an 710065,China
  • 5. Research Institute of Exploration and Development,PetroChina Xinjiang Oilfield Company,Karamay 834000,China

Received date: 2023-07-01

  Revised date: 2023-09-11

  Online published: 2023-11-24

Supported by

The National Natural Science Foundation of China(41502132)

Abstract

The paleo-pressure during the accumulation period exerts a significant influence on the distribution of oil and gas, yet its origin is considerably more challenging to identify and analyze compared to present pressure. In this paper, a quantitative analysis model was proposed to investigate the existence of paleo-overpressure in the 7th member of Yanchang Formation (Chang 7 Member) of Ordos Basin and its logging response mechanism, with a focus on the influence of organic matter content. The origin of paleo-overpressure and its significance in hydrocarbon accumulation within the Chang 7 Member of the study area were discussed by employing the inversion analysis based on the classical origin identification template and the forward analysis based on the basin modeling. The findings indicate that the abundant organic matter content in Chang 7 Member exerts a significant influence on logging parameters. The quantitative analysis reveals that the increment of acoustic time difference and the reduction of density induced by organic matter are 2.56-169.78 μs/m and 0.01-0.61 g/cm3, respectively. Based on the data of acoustic time difference and density with the removal of the contribution of organic matter, and used the methods of calculated porosity contrast and sonic velocity-density crossplot, it was concluded that there is paleo-overpressure in Chang 7 Member, and the fluid expansion and pressurization caused by hydrocarbon generation is the main origin of overpressure during the accumulation period. Combined with basin modeling, it was pointed out that the overpressure due to hydrocarbon generation matches the thermal evolution of source rock in time and hydrocarbon distribution in space, which further supports the understanding that the overpressure results from the fluid expansion by hydrocarbon generation. It was considered that paleo-pressure is one of the important controlling factors of reservoir formation, and hydrocarbon accumulation is distributed in the local low excess pressure area in the high excess pressure area.

Cite this article

Zeyang XU , Zhen LIU , Jingzhou ZHAO , Jiacheng DANG , Jun LI , Ziyi TANG . Origin analysis of paleo-overpressure of source rock based on logging response characteristics: Case study of the source rocks in the 7th member of Yanchang Formation of Triassic, Ordos Basin, China[J]. Natural Gas Geoscience, 2023 , 34(11) : 1950 -1960 . DOI: 10.11764/j.issn.1672-1926.2023.09.009

0 引言

超压是含油气盆地中较为普遍存在的地质现象1-2,其成因及分布一直是石油地质研究人员关注的热点及难点3-5。最初研究地层压力主要是为了满足工程上的需要,以确保钻井安全并节约成本6-8。随着油气成藏动力学的发展,地下异常压力逐渐成为探讨油气运移、充注及保存的重要研究内容之一9-11。在以往的油气勘探中,人们倾向于关注现今地层压力分布特征。相应地,多种基于地球物理响应的异常压力定量计算模型也相继被提出12-15,并在许多地区得到成功应用。然而,从动态分析的角度看,成藏时期的古压力特征对于油气分布同样存在重要影响16,而且对其识别及成因分析的难度更大。此外,现有的超压成因判识方法大都基于不同类型测井参数在超压层段的异常响应17-18。然而受岩石成分、流体性质以及井壁扩径等多种因素的干扰,测井参数具有很强的多解性19-20。尤其在具有生烃能力的暗色泥岩层段中,高声波时差、低密度的固体有机质将较大程度地影响超压成因判识结果。
延长组7段(长7段)是鄂尔多斯盆地三叠系主要的烃源岩层段之一,其发育一套深湖—半深湖相暗色泥页岩21。盆地西南部城壕、镇泾等地区现今长7段地层压力普遍等于或小于对应深度处的静水压力22,前人几乎很少关注其压力特征,甚至在钻井、录井及试油等工程资料中也鲜有压力方面的报道。然而,从测井资料来看,长7段存在明显的声波时差和密度异常。通常情况下,人们很容易把这种富有机质层段的测井参数异常简单地归结为骨架组成的差异所致23。当然,也有学者认为长7段底部的声波时差异常增大是由于历史时期的不均衡压实造成的24。李超等25较早指出鄂尔多斯盆地长7段泥岩测井声波时差异常与有机质含量存在较大关系,并提出了相应的校正方法。事实上,文献[24]的结论就是在校正基础上得出的。然而,仅对声波时差进行校正似乎难以体现出超压对测井参数的影响程度,可能还需引入密度这类体积属性进行对比分析。笔者目前尚未看到对鄂尔多斯盆地长7段的测井密度响应量进行定量计算的相关研究。另外,以文献[17]的发表为分水岭,如今我国许多学者对大多数沉积盆地超压成因的认识已经或正在发生转变。也有不少学者推测长7段存在过古超压,且其成因为生烃增压26-28。但该认识主要是前人根据简单的推测得出的,大都没有经过严格的论证。
徐泽阳等29在对松辽盆地长垣地区白垩系青山口组一段超压进行分析时,提出了剔除高声波时差、低密度固体有机质影响的校正方法。该方法虽然同时考虑了声波时差和密度2类参数与有机质含量的关系,但其计算模型中的有机质部分仅考虑了固体组分,这可能导致出现校正不彻底的问题。因此,本文在上述方法的基础上,同时考虑了固体有机质以及烃类流体含量对声波时差和密度的影响,并提出了相应的声波时差增量及密度减量定量计算模型。根据计算结果,进一步讨论了研究区长7段是否存在或曾经存在过超压。此外,采用赵靖舟等17总结的多个超压成因判识方法,本文选取计算孔隙度对比法、声波—密度交会图法以及地质综合分析法对鄂尔多斯盆地三叠系长7段超压成因进行讨论。根据讨论结果,选取相应的恢复模型对长7段古压力演化过程进行了模拟,进一步佐证了古超压存在的认识,并探讨了成藏期古压力分布特征对烃类聚集的影响。

1 地质背景

研究区位于鄂尔多斯盆地西南部(图1),研究层段主要为三叠系延长组7段和8段。该地区自延长组沉积以来,先后经历了印支、燕山及喜马拉雅多期构造运动30-31。其中,印支期是烃源岩、储集层形成时期,燕山期是主要的油气成藏期,喜马拉雅期则对应盆地的整体抬升期。长7段属于上三叠统延长组中有机质含量最高、厚度最大的有利烃源岩层段。长8段储集层主要发育三角洲前缘分流河道砂体,从物性上来看部分储集层已达致密,覆压孔隙度及渗透率平均值分别为8%~15%和(0.05~0.5)×10-3 μm2
图1 研究区概况

(a)研究区位置;(b)关键井位图;(c)上三叠统地层柱状图

Fig.1 Overview of the study area

一般认为长8段油藏是由长7段石油“倒灌”形成的,这种情况往往需要较强的运移动力。因此,运聚期的流体动力条件很可能是决定烃类分布及富集程度的重要因素之一。此外,烃类流体需要较强的充注动力才能进入长8段这类致密或近乎致密的储集层中。然而,从现有的压力资料及文献来看,研究层段上下范围内地层压力以常压为主,局部存在微弱负压。这说明研究区曾经很可能存在过超压,但在后期演化过程中逐渐散失。

2 有机质测井响应量分析

在正常压实趋势下,随着深度的增大,声波时差和中子值将逐渐减小,而密度和电阻率值将逐渐增大。但在研究区长7段测井剖面中,声波时差、密度、中子及电阻率皆显示出与正常压实趋势相反的变化规律(图2)。一般情况下,出现测井曲线异常时,人们往往优先考虑到地层中可能存在异常压力。然而如前文所述,研究区长7段现今地层压力为常压或微弱负压,此处的异常似乎与超压无关。观察测井曲线及实测TOC含量在剖面上的变化特征后可以看出,4类测井曲线与TOC含量都存在一定的正相关关系。其中,电阻率曲线与TOC含量的关系最好。这说明上述异常可能是由暗色泥岩中的有机质引起的。而在里189井中,埋深2 200 m处的声波时差和密度曲线存在明显峰值,但TOC含量并没有发生相应的变化[图2(a)]。这或许暗示有机质含量并不一定是测井参数异常的唯一主控因素。参考相关钻井及录井资料,在该部位并未出现卡钻、扩径、井壁垮塌以及钻井液漏失等异常情况。屈童等28通过包裹体分析及盆地模拟,推测长7段存在过古超压,并指出其可能是造成测井参数异常的重要因素。但这一研究过程中存在明显矛盾,其在盆地模拟时认为古超压是生烃增压所致,但在利用测井参数估计古过剩压力时又选取了偏向于欠压实成因的等效深度模型。出现这一矛盾可能是因为前人分析古超压成因时选取的测井曲线类型过于单一。
图2 研究区典型井长7段测井响应特征

(a)里189井;(b)坪114井;(c)午121井

Fig.2 Log response characteristics of the Chang 7 Member in the typical wells of the study area

为讨论有机质含量对声波时差和密度的影响量,本文将烃源岩划分为非有机质骨架、固体有机质、烃类流体和地层水4部分。则单位体积的烃源岩多组分等效模型可表示为:
t = t m a V m a + t o m V o m + t h V h + t w V w
ρ = ρ m a V m a + ρ o m V o m + ρ h V h + ρ w V w
V m a + V o m + V h + V w = 1
V o m = φ T O C 1 - φ
V h + V w = φ
V h = S h × φ
在未考虑有机质含量时,一般将岩石看作由非有机质骨架和孔隙流体2部分构成,且孔隙流体一般被认为主要是地层水。故仅考虑岩石骨架和地层水2部分时,单位体积非烃源岩物理体积模型可表示为:
t ' = t m a V m a ' + t w V w '
ρ ' = ρ m a V m a ' + ρ w V w '
V m a ' + V w ' = 1
V w ' = φ
需要注意的是,孔隙度相同的情况下,考虑和不考虑有机质含量时的非有机质骨架体积和地层水体积的含义是不同的,其相互关系可参考图3。将式(1)式(10)联立,则可得到由有机质含量引起的声波时差增量和密度减量定量计算模型:
t T O C = t - t ' = φ T O C 1 - φ t o m - t m a + S h φ t h - t w
ρ T O C = ρ ' - ρ = φ T O C 1 - φ ρ m a - ρ o m + S h φ ρ w - ρ h
图3 考虑和不考虑有机质含量情况下各组分体积含义对比

Fig.3 The volume meaning comparison of each component in the rocks with and without organic matter content

式(1)式(11)中: t t '分别为考虑有机质和不考虑有机质时的声波时差,μs/m; ρ ρ '分别为考虑有机质和不考虑有机质时的岩石密度,g/cm3 t T O C为有机质含量引起的声波时差增量,μs/m; ρ T O C为有机质含量引起的密度减量,g/cm3 t m a t o m t h t w分别为非有机质骨架、固体有机质、烃类流体和地层水的声波时差,分别取地区经验常数223.39 μs/m、550 μs/m、780 μs/m和625 μs/m2532 ρ m a ρ o m ρ h ρ w分别为非有机质骨架、固体有机质、烃类流体和地层水的密度,分别取经验值2.68 g/cm3、1.2 g/cm3、0.8 g/cm3和1.05 g/cm3 V m a V o m V h V w分别为考虑有机质含量情况下单位体积中非有机质骨架、固体有机质、烃类流体和地层水所占体积,无量纲; V m a ' V w '分别为不考虑有机质含量情况下单位体积中非有机质骨架和地层水所占体积,无量纲; φ为泥岩孔隙度,无量纲; φ T O C为固体有机质在泥岩骨架中的体积占比,无量纲; S h为烃类流体在孔隙流体中的体积占比,即含烃饱和度,取区域平均值0.5833
注意到,应用式(11)式(12)计算时,还需要确定 φ V o m。当计算烃源岩孔隙度时,本文推荐使用徐泽阳等29提出的富有机质泥岩孔隙度计算模型,即:
φ A C = 1 / C P t - 1 - φ T O C t m a - φ T O C t o m t f - 1 - φ T O C t m a - φ T O C t o m
φ D E N = 1 - φ T O C ρ m a + φ T O C ρ o m - ρ 1 - φ T O C ρ m a + φ T O C ρ o m - ρ f
t f = t h S h + t w 1 - S h
ρ f = ρ h S h + ρ w 1 - S h
式(13)式(16)中: φ A C为考虑有机质含量时用声波时差计算出的泥岩孔隙度,无量纲; φ D E N为考虑有机质含量时用密度计算出的泥岩孔隙度,无量纲; t f为孔隙流体声波时差,μs/m; ρ f为孔隙流体密度,g/cm3 C P为压实校正系数,取经验值1.5;其余参数意义与式(1)式(11)一致。
为避免混淆,此处需强调式(11)式(12)中的 φ分别是指在考虑有机质的情况下声波时差或密度计算出的孔隙度,二者并非真实的泥岩孔隙度。在正常压实情况下,计算值与实测值可能是相近的;而在超压段中,两者存在差异与否可能与超压的成因有很大关系。因此,在本文的相关计算中,一般是将式(13) φ A C代入式(11) φ中,将式(14) φ D E N代入式(12) φ中,进而完成对 φ的消元处理。
在确定 V o m时,本文首先参考前人在鄂尔多斯盆地的经验公式34-35计算有机质含量:
Δ L g R = Δ L g R / R B L + 0.006   1 t - t B L
w T O C = Δ L g R × 10 2.296 - 0.168   8 × 11
式中: Δ L g R为声波时差曲线与电阻率曲线的间距,无量纲; R为泥岩电阻率,Ω·m; R B L为纯泥岩段对应的电阻率,Ω·m; t为泥岩声波时差,μs/m; t B L为纯泥岩段对应的声波时差,μs/m; w T O C为固体有机质在泥岩骨架中的质量占比,%。
此处需要指出的是,电阻率同时受到了压力与有机质含量的影响,似乎 Δ L g R法难以在剔除压力影响的前提下估计有机质含量。但实际上, Δ L g R法中的计算模型及参数是根据实测有机质含量与电阻率之间的经验关系确定的,因此从统计学角度来看,压力的影响量在此处已经作为系统量被统一剔除。事实上,在诸多存在超压的地区, Δ L g R法普遍取得了令人较为满意的结果36-38。当然,作为一种较为简单的统计学方法,这种方法的精度会受到有机质含量以外的复杂因素影响。LIU等39引入了SVR模型对延长组张家滩页岩TOC进行了估算并取得了较为理想的结果。SVR模型本质上是一种更复杂的统计学方法,需要更多实测资料支撑计算机学习,其预测精度与测试资料的丰富程度呈正相关关系,显然更加适用于TOC的预测。遗憾的是,本文研究所掌握的资料不足以支撑使用这个方法,但笔者推荐在条件允许的情况下采用更加准确的TOC预测方法。
此外,在应用本文所建的模型时,需要使用体积百分比。可根据下式进行转换:
V o m = 1 - φ ρ m a ρ o m k w T O C / 100
式中: w T O C为固体有机质在泥岩骨架中的质量占比,%; ρ m a为非有机质泥岩骨架密度,2.68 g/cm3 ρ o m为固体有机质密度,取1.2 g/cm3 φ为孔隙度,无量纲; k为有机碳转化系数,与有机质类型和成熟度有关32,研究区Ⅰ型干酪根取1.25。
应用式(11)式(14)式(17)式(18)对68个泥岩样品点进行计算,结果显示长7段泥页岩有机质含量引起的声波时差增量在2.56~169.78 μs/m之间,主要分布在2.56~80 μs/m之间;有机质含量引起的密度减量在0.01~0.61 g/cm3之间,主要分布在0.01~0.4 g/cm3之间。通过对比总异常量与有机质引起的增量或减量分布特征可以看出,声波时差总异常量主峰范围明显高于有机质引起的增量范围,而密度总异常量分布与有机质引起的增量范围相似(图4)。这说明除有机质含量之外还有其他因素引起声波时差异常,但密度的异常似乎仅与有机质含量有关。虽然笔者目前还未掌握研究区长7段流体包裹体压力分析的相关资料,但前人应用盆地模拟40、包裹体分析41等手段,已经证实了邻区长7段曾存在过古超压。由于古超压往往会撑开连通孔并导致微裂缝发育,而这些微裂缝及其对测井参数的影响在撤去超压后往往难以彻底消失,因此笔者推测研究区声波时差这类对微裂缝较为敏感的参数出现异常与古超压的存在关系较大。
图4 研究区长7段有机质含量对声波时差及密度影响量分布特征

Fig.4 Distribution characteristics of the influence of organic matter content in Chang 7 Member of the study area on acoustic time difference and density

3 古超压成因分析

在上述定量计算的基础上,本文对声波时差和密度受有机质含量的影响量进行了校正,并进一步应用校正后结果进行了古超压成因的分析。
HERMANRUD等42指出超压段声波时差、密度、电阻率、中子等与孔隙度密切相关的测井数据会出现明显异常。然而,仅根据测井参数是否发生反转来判识古超压成因显然是不够的。笔者较为认可BOWERS43对测井资料的划分方案。其中,声波时差和电阻率这类测井数据被划分为传导属性,密度和中子这类测井数据被划分为体积属性。一般认为,不均衡压实将导致传导属性和体积属性发生同幅度的响应;而流体膨胀将通过增加微裂缝、改变孔喉结构的方式较大程度地影响传导属性,但对体积属性的影响较小。
由于各测井参数物理意义及单位各不相同,不同物理参数之间的异常幅度往往难以对比。为解决这一问题,TINGAY44将声波时差及密度数据同时换算为相应的孔隙度进行对比分析。在不均衡压实成因的超压段中,传导属性与体积属性计算出的孔隙度异常幅度一致;在流体膨胀或传导成因的超压段中,体积属性计算的孔隙度异常幅度相对较低或几乎没有异常。这种归一化的办法避免了主观影响且操作简单,得到了许多学者的认可31745。为对比有机质含量校正前后的超压成因判识方法,除使用式(13)式(14)计算孔隙度外,本文还应用Willy公式计算了不考虑有机质含量情况下得到的孔隙度:
φ A C ' = 1 / C P t - t m a t w - t m a
φ D E N ' = ρ m a - ρ ρ m a - ρ w
式中: φ A C '为不考虑有机质含量时用声波时差计算出的泥岩孔隙度,无量纲; φ D E N '为不考虑有机质含量时用密度计算出的泥岩孔隙度,无量纲;其余参数意义与式(1)式(14)一致。
图5展示了里97井长7段及其相邻泥岩层段校正前后声波时差孔隙度和密度孔隙度的差异。可以看出校正前,2类孔隙度都较大幅度地偏离正常压实趋势,表现出不均衡压实的特征。得益于页岩油勘探事业的发展,城壕地区长7段有许多实测的泥岩物性资料。据统计,长7段岩石样品孔隙度介于0.85%~9.57%之间,其中泥页岩样品平均孔隙度为1.55%46。虽然校正后的密度孔隙度与实测孔隙度均值仍有一定差距,但基本已经落在实测值分布区间内。而校正后的声波孔隙度与实测值存在较大程度的差异。笔者认为,历史时期流体膨胀导致的裂缝是上述现象出现的重要原因。
图5 研究区里95井有机质含量校正前后计算孔隙度对比剖面

Fig.5 Calculated porosity comparison profile before and after organic matter content correction in Well Li 95 in the study area

除此之外,还有一种更为直接的办法,即将体积属性与传导属性做交会分析。一般选择声波速度与密度进行对比。加载过程中声波速度与密度沿加载曲线同步变化。因不均衡压实仅为压实过程的停顿,故声波速度与密度仍然落在加载曲线上;由于卸载成因超压段中声波速度变化幅度远大于密度,故而数据点多落在加载曲线之下[图6(a)]。校正前后,研究区长7段超压数据点呈现出不同的分布特征。校正前,长7段超压点落在反映不均衡压实的加载曲线上;校正后,长7段大部分投影点位于流体膨胀区域[图6(b)]。
图6 研究区典型井声波速度—密度交会图

(a)声波速度—密度交会判识图版47-49;(b)里95井和元429井声波速度—密度交会图

Fig.6 Velocity-density crossplot of the typical wells in the study area

根据3种方法综合分析认为,剔除有机质含量影响后,传导属性仍然存在明显的异常。这些现象证实了测井参数的异常并不是单单受有机质含量控制,其预示长7段可能存在古超压,且揭示出古超压成因为生烃作用导致的流体膨胀。这一认识为后续古压力恢复工作中计算模型的选取提供了重要依据。

4 古压力恢复及成藏意义

前文对古超压成因的分析主要基于测井参数的响应,属于地质反演分析的范畴。本节在前文结论的基础上,选取基于生烃增压的古压力正演模型对长7段压力演化过程进行了分析。
现在常用的PetroMod软件在恢复古压力时,主要考虑了不均衡压实和生烃增压2类超压成因。软件中定义了深度、有效应力和孔隙度的相互关系,同时考虑了流体膨胀情况下的流体压缩系数。以一维模拟为例,针对不均衡压实成因的超压,软件会根据目的层上下岩层的渗透率是否达到阈值而识别流体是否排出受限,继而通过之前提到的自定义关系进行压力计算。当模拟识别出生烃过程时,软件会根据生烃量、孔隙度和流体压缩系数综合判断增压量,并根据提前定义的破裂压力阈值确定何时泄压。
具体而言,软件中不均衡压实成因超压恢复的核心为Athy模型50
φ = φ 1 + φ 0 - φ 1 e - k σ
式中: φ为孔隙度,无量纲; φ 0为初始孔隙度,无量纲; φ 1为最小孔隙度,无量纲; σ为有效应力,MPa; k为经验系数,无量纲。
而生烃增压的核心模型为:
P = K V φ ρ f C 1 - ρ f ρ k
式中: P为剩余压力,MPa; V为岩石体积,cm3 φ为孔隙度,无量纲; ρ f为流体密度,g/cm3 ρ k为干酪根密度,g/cm3 K为干酪根转化量,g; C为压缩系数,水约为0.5 GPa-1,油约为1 GPa-1
红河105井位于镇泾地区中部,在长7段钻遇约7 m厚暗色泥岩,其平均TOC值为2.59%,平均氢指数为203 mg/g,平均S 1+S 2值为10.22 mg/g,成熟度R O值在0.72%~0.93%之间。为对比不均衡压实作用以及生烃增压作用在古压力演化过程中的作用,本文在对红河105井长7段暗色泥岩演化过程模拟时分别作出不同的设置。第一种是存在生烃过程的设置,将长7段指定为烃源岩层段并将TOC及氢指数输入软件中;第二种是不存在生烃过程的设置,将长7段指定为普通泥岩层段。从恢复结果来看(图7),在侏罗纪长7段孔隙度变化较小,一般稳定在20%~30%之间。这似乎非常符合不均衡压实超压形成条件,但该时期的多次抬升可能不利于这类成因超压的保存。事实上,恢复结果显示在侏罗纪长7段并未出现过超压,其也可能与最大埋深一直较小有关。因为有效应力不变只是出现不均衡压实超压的条件之一,压力异常增大还需要总应力的不断增加。显然,侏罗纪多次抬升使得总应力表现出多次减小的特征。即使在沉降过程中有所增加,总应力也一直无法超过该时期最小孔隙度形成时对应的上覆地层压力。虽然侏罗纪之后埋藏深度开始增大,但孔隙度也发生明显减小,此时不均衡压实的典型特征也随之消失。此外,对比2种不同设置的古压力恢复结果,存在生烃作用的情况下,早白垩世至新近纪地层压力明显高于静水压力,且地层压力与静水压力开始分离的时间刚好对应生烃量开始快速增加的时期;而对于不存在生烃作用的情况,其地层压力始终与静水压力一致。因此,本文认为演化阶段早期长7段经历多次抬升,不均衡压实超压的保存条件不佳;而在演化过程的中后期,超压的主要来源为生烃增压。从烃源岩热演化史来看,白垩纪以来长7段进入生烃阶段。相应的,地层压力开始偏离静水压力。晚白垩世长7段达到生烃高峰的同时,地层剩余压力也达到峰值,此时正对应延长组大多数油藏的形成时期。在生烃强度较大的地区,剩余压力可达5~10 MPa。从形成时间和压力大小来看,生烃增压完全能够满足成藏对运移及充注动力的需求。而现今弱负压—常压的特征,主要是之后长期抬升过程导致的压力散失所致。
图7 研究区红105井长7段地层压力及相关要素演化过程

Fig.7 Evolution process of formation pressure and related factors in Chang 7 Member of Well Hong 105 in the study area

生烃增压可作为成藏动力的另一重要证据是现今油藏与古压力的分布关系。本文根据盆地模拟结果,绘制了镇泾地区长7段古压力系数分布图。可以看出,该地区古压力呈现出“东南高、西北低”的分布特征,而规模较大的油藏基本集中分布在古压力系数大于1.27的地区。同时,烃类的运移方向一般由高剩余压力区指向低剩余压力区,故局部的相对低压区能够形成有利的烃类聚集区。显然,镇泾地区长8段油藏分布基本符合这一特征,其受砂体控制的同时,主要受到了古压力分布的影响(图8)。此外,笔者从式(23)模型推断,软件给出的白垩纪以来尤其是在晚白垩世的剩余压力主要来自于生烃增压。因此,由恢复的古压力与油藏分布相吻合这一特征,也能够印证古超压的成因确实为生烃导致的流体膨胀。
图8 镇泾地区长7段晚白垩世古压力系数与长8段油藏、砂体叠合图

Fig.8 Late Cretaceous paleo-pressure coefficient of Chang 7 Member, Zhenjing area, superposed with reservoir and sand body of Chang 8 Member

5 结论

本文根据鄂尔多斯盆地西南部长7段存在测井参数偏离正常压实趋势的现象,引发了对研究区是否存在古超压的探讨,并提出了一套基于测井参数校正及盆地模拟的古超压识别及成因分析方法。主要取得了以下三方面认识:
(1)本文提出了有机质含量响应定量计算模型,并对泥岩声波时差和密度进行了分析及校正。结果显示由有机质含量引起的声波时差增量为2.56~169.78 μs/m,主要分布在2.56~80 μs/m之间;由有机质含量引起的密度减量在0.01~0.61 g/cm3之间,主要分布在0.01~0.4 g/cm3之间。其中,声波时差总异常量分布峰值明显高于有机质含量引起的声波时差增量分布峰值;密度总异常量分布区间及峰值与有机质含量引起的密度减量相似。这说明有机质含量并非引起测井参数异常的唯一原因,且研究区长7段很可能存在古超压。
(2)剔除有机质含量影响后,结合计算孔隙度对比法以及声波速度—密度交会图法对长7段古超压成因进行了分析。校正前,各方法的超压成因判识结果皆指向不均衡压实成因;校正后,各方法的判识结果都指向流体膨胀成因。
(3)应用盆地模拟软件对长7段古压力进行恢复,认为长7段不具备形成不均衡压实超压的条件。同时,根据烃源岩热演化与超压形成过程在时间上的耦合关系,推断长7段超压主要源自生烃增压。指出成藏期古压力分布与油藏在空间上的耦合关系,认为其既反映古压力对油藏分布的控制作用,又进一步佐证了生烃导致的流体膨胀是古超压形成的主要原因。
1
HUNT J M.Generation and migration of petroleum from abno-rmally pressured fluid compartments[J].AAPG Bulletin,1990, 74(1):1-12.

2
张启明,董伟良.中国含油气盆地中的超压体系[J].石油学报,2000,21(6):1-11,127.

ZHANG Q M,DONG W L. Overpressure system in oil and gas basins in China[J].Acta Petrolei Sinica,2000,21(6):1-11,127.

3
李军,唐勇,吴涛,等.准噶尔盆地玛湖凹陷砾岩大油区超压成因及其油气成藏效应[J].石油勘探与开发,2020,47(4):679-690.

LI J,TANG Y,WU T,et al. Overpressure origin and its effects on petroleum accumulation in the conglomerate oil province in Mahu Sag, Junggar Basin, NW China[J]. Petroleum Exploration and Development,2020,47(4):679-690.

4
郝芳,董伟良.沉积盆地超压系统演化、流体流动与成藏机理[J].地球科学进展,2001,16(1):79-85.

HAO F,DONG W L. Evolution of,fluid flow and petroleum accumulation in overpressured systems in sedimentary basins[J]. Advances in Earth Science,2001,16(1):79-85.

5
李明诚,李剑.“动力圈闭”——低渗透致密储层中油气充注成藏的主要作用[J].石油学报,2010,31(5):718-722.

LI M C,LI J. “Dynamic trap”: A main action of hydrocarbon charging to form accumulations in low permeability-tight reservoir[J]. Acta Petrolei Sinica,2010,31(5):718-722.

6
SAYERS C M,JOHNSON G M,DENYER G. Predrill pore-pressure prediction using seismic data[J]. Geophysics,2002,67(4):1286-1292.

7
AHMED A,ELKATATNY S,ALI A, et al. Comparative analysis of artificial intelligence techniques for formation pressure prediction while drilling[J].Arabian Journal of Geosciences,2019, 12:1-13.

8
ZHANG J. Pore pressure prediction from well logs:Methods, modifications,and new approaches[J].Earth-Science Reviews,2011,108(1-2):50-63.

9
郝芳,邹华耀,姜建群.油气成藏动力学及其研究进展[J].地学前缘,2000,7(3):11-21.

HAO F,ZOU H Y,JIANG J Q. Dynamics of petroleum accumulation and its advances[J]. Earth Science Frontiers,2000,7(3):11-21.

10
迟元林,萧德铭,殷进垠.松辽盆地三肇地区上生下储“注入式”成藏机制[J].地质学报,2000,74(4):371-377.

CHI Y L,XIAO D M,YIN J Y.The injection pattern of oil and gas migration and accumulation in the Sanzhao area of Songliao Basin[J]. Acta Geologica Sinica,2000,74(4):371-377.

11
杨智,何生,何治亮,等.准噶尔盆地腹部超压层分布与油气成藏[J].石油学报,2008,29(2):199-205,212.

YANG Z,HE S,HE Z L,et al. Distribution of overpressure stratum and its relationship with hydrocarbon accumulation in the central part of Junggar Basin[J].Acta Petrolei Sinica,2008,29(2):199-205,212.

12
EATON B A. The effect of overburden stress on geopressure prediction from well logs[J]. Journal of Petroleum Technology, 1972, 24(8):929-934.

13
MAGARA K. Importance of aquathermal pressuring effect in Gulf Coast[J]. AAPG Bulletin, 1975, 59(10):2037-2045.

14
FILLIPPONE W R. On the prediction of abnormally pressured sedimentary rocks from seismic data[C]. Offshore Te-chnology Conference,OnePetro,1979.

15
BOWERS G L. Pore pressure estimation from velocity data:Accounting for overpressure mechanisms besides undercompaction[J]. SPE Drilling & Completion,1995,10(2):89-95.

16
卢朝进. 深埋型和抬升型致密砂岩孔隙动力演化特征研究[D].北京:中国石油大学(北京),2020:1-136.

LU Z J. Differential Evolution Characteristics between Buried and Uplift of Tight Sandstone on the Porodynamics[D]. Beijing:China University of Petroleum(Beijing),2020:1-136.

17
赵靖舟,李军,徐泽阳.沉积盆地超压成因研究进展[J].石油学报,2017,38(9):973-998.

ZHAO J Z,LI J,XU Z Y. Advances in the origin of overpressures in sedimentary basins[J]. Acta Petrolei Sinica,2017,38(9):973-998.

18
何生,何治亮,杨智,等.准噶尔盆地腹部侏罗系超压特征和测井响应以及成因[J].地球科学(中国地质大学学报),2009,34(3):457-470.

HE S,HE Z L,YANG Z,et al. Characteristics, Well log responses and mechanisms of overpressures within the Jurassic Formation in the central part of Junggar Basin[J]. Earth Science(Journal of China University of Geosciences),2009,34(3):457-470.

19
HAN D, NUR A, MORGAN D. Effects of porosity and clay content on wave velocities in sandstones[J]. Geophysics, 1986, 51(11):2093-2107.

20
樊洪海,张传进.复杂地层地层孔隙压力求取新技术[J].石油钻探技术,2005,33(5):43-46.

FAN H H,ZHANG C J. New methods for calculation of pore pressure in complex geologic environment[J]. Petroleum Drilling Techniques,2005,33(5):43-46.

21
袁选俊,林森虎,刘群,等.湖盆细粒沉积特征与富有机质页岩分布模式——以鄂尔多斯盆地延长组长7油层组为例[J].石油勘探与开发,2015,42(1):34-43.

YUAN X J,LIN S H,LIU Q,et al. Lacustrine fine-grained sedimentary features and organic-rich shale distribution pattern: A case study of Chang 7 Member of Triassic Yanchang Formation in Ordos Basin, NW China[J]. Petroleum Exploration and Development,2015,42(1):34-43.

22
李士祥. 鄂尔多斯盆地中生界异常低压成因及对成藏的影响[D].成都:成都理工大学,2017:1-119.

LI S X. Effects on Petroleum Accumulation and Genesis of the Mesozoic Abnormal Low Pressure in Ordos Basin[D].Chengdu:Chengdu University of Technology,2017:1-119.

23
郑奎,杨晋玉,胡晓雪,等.胡尖山—姬塬地区长7段油页岩定性定量评价及长9油藏勘探评价[J].长江大学学报(自然科学版),2022,19(5):27-36,80.

ZHENG K,YANG J Y,HU X X,et al. Qualitative and quantitative evaluation of Chang 7 oil shale and exploration evaluation of Chang 9 reservior in Hujianshan-Jiyuan area[J]. Journal of Yangtze University(Natural Science Edition),2022,19(5):27-36,80.

24
马立元,邱桂强,李超,等.鄂尔多斯盆地镇泾地区延长组异常压力演化及其成藏意义[J].地质通报,2020,39(4):503-511.

MA L Y,QIU G Q,LI C,et al. The evolution of abnormal pressure of Yanchang Formation in Zhenjing area of Ordos Basin and its reservoir-forming significance[J]. Geological Bulletin of China,2020,39(4):503-511.

25
李超,张立宽,罗晓容,等.泥岩压实研究中有机质导致声波时差异常的定量校正方法[J].中国石油大学学报(自然科学版),2016,40(3):77-87.

LI C,ZHANG L K,LUO X R,et al. A quantitative method for revising abnormally high sonic data in rich-organic rock during compaction study[J]. Journal of China University of Petroleum(Edition of Natural Science),2016,40(3):77-87.

26
王帅.延长油田长7段致密油生烃增压及其对成藏的控制[D].北京:中国石油大学(北京),2018:1-54.

WANG S. Overpressure Evolution Caused by Hydrocarbon Generation of Chang 7 Section in Yanchang Oilfield and its control Effect on Accumulation[D]. Beijing:China University of Petroleum, Beijing,2018:1-54.

27
李超正.鄂尔多斯盆地长7段致密砂岩储层石油充注有效性研究[D].北京:中国石油大学(北京),2020:1-126.

LI C Z. Research on Effectiveness of Oil Charging in the Tight Sandstone Reservoir of the Chang 7 Member,Ordos Basin[D]. Beijing:China University of Petroleum(Beijing),2020:1-126.

28
屈童,高岗,梁晓伟,等.鄂尔多斯盆地长7段致密油成藏机理分析[J].地质学报,2022,96(2):616-629.

QU T,GAO G,LIANG X W,et al. Analysis of tight oil accumulation mechanism of Chang 7 member in the Ordos Basin[J]. Acta Geologica Sinica,2022,96(2):616-629.

29
徐泽阳,赵靖舟,李军.松辽盆地长垣地区白垩系青山口组一段有机质含量对超压分析的影响及校正方法[J].石油与天然气地质,2019,40(4):938-946.

XU Z Y,ZHAO J Z,LI J. The impact of organic matter content on overpressure analysis and its correction method in the first member of Cretaceous Qingshankou Formation,Placanticline area,Songliao Basin[J]. Oil & Gas Geology,2019,40(4):938-946.

30
王峰,田景春,范立勇,等.鄂尔多斯盆地三叠系延长组沉积充填演化及其对印支构造运动的响应[J].天然气地球科学,2010,21(6):882-889.

WANG F,TIAN J C,FAN L Y,et al. Evolution of sedimentary fillings in Triassic Yanchang Formation and its response to Indosinian movement in Ordos Basin[J]. Natural Gas Geoscience,2010,21(6):882-889.

31
张金山.鄂尔多斯盆地构造演化及含油气展望[J].石油与天然气地质,1982,3(4):304-315.

ZHANG J S. Tectonic evolution and oil-bearing prospect of Ordos Basin[J]. Oil & Gas Geology,1982,3(4):304-315.

32
阎东启,柳波,贾梦成,等.高有机质丰度泥岩压实的计算误区及定量校正方法[J].断块油气田,2018,25(1):29-33.

YAN D Q,LIU B,JIA M C,et al. Misunderstandings of high organic matter abundance mudstone compaction and quantitative correction methods[J]. Fault-Block Oil & Gas Field,2018,25(1):29-33.

33
吴勇,侯雨庭,李会庚,等.高精度地层元素测井在页岩油储层评价中的应用——以鄂尔多斯盆地长7段为例[J].非常规油气,2016,3(2):1-10.

WU Y,HOU Y T,LI H G,et al. Application of high-precision formation elements logging to shale oil reservoir evaluation: A case study of seventh member of Yanchang Formation, Ordos Basin[J]. Unconventional Oil & Gas,2016,3(2):1-10.

34
杜江民,张小莉,钟高润,等.致密油烃源岩有机碳含量测井评价方法优选及应用——以鄂尔多斯盆地延长组长7段烃源岩为例[J].地球物理学进展,2016,31(06):2526-2533.

DU J M,ZHANG X L,ZHONG G R,et al. Analysis on the optimization and application of well logs indentification methods for organic carbon content in source rocks of the tight oil: Illustrated by the example of the source rocks of Chang 7 member of Yanchang Formation in Ordos Basin[J].Progress in Geophysics,2016,31(6):2526-2533.

35
柏杨. 致密油烃源岩TOC测井定量预测模型优选及适用规律研究[D].石家庄:河北地质大学,2021:1-73.

BAI Y. Optimization and Application of TOC Logging Quantitative Prediction Model for Tight Oil Source Rocks[D]. Shijiazhuang:Hebei GEO University,2021:1-73.

36
蔡冬梅,赵弟江,彭靖松,等.优质烃源岩识别及其多属性反演技术定量评价——以渤海海域辽东南洼陷为例[J].石油地球物理勘探,2018,53(2):330-338,222-223.

CAI D M,ZHAO D J,PENG J S,et al. High-quality hydrocarbon source rock identification and quantitative evaluation with multi-attributein-version: A case study of Liaodong south sub-sag, Bohai Bay Basin[J]. Oil Geophysical Prospecting,2018,53(2):330-338,222-223.

37
李思,代红霞,胡燕婷,等.基于广义ΔlgR技术的致密气源岩有机碳测井评价[J].长江大学学报(自科版),2018,15(23):36-41,6.

LI S,DAI H X,HU Y T,et al. Evaluation of organic carbon logging of tight gas source rocks basedon generalized-ΔlgR technology[J]. Journal of Yangtze University(Natural Science Edition),2018,15(23):36-41,6.

38
樊妍,王绪龙,向才富,等.准噶尔盆地东部二叠系平地泉组烃源岩富集规律与主控因素[J].现代地质,2022,36(4):1105-1117.

FAN Y,WANG X L,XIANG C F,et al. Enrichment patterns and main controlling factors of source rocks in the Permian Pingdiquan Formation,Eastern Junggar Basin[J]. Geoscience,2022,36(4):1105-1117.

39
LIU X, LEI Y, LUO X, et al. TOC determination of Zhangjiatan shale of Yanchang Formation, Ordos Basin, China, using support vector regression and well logs[J]. Earth Science Informatics,2021,14:1033-1045.

40
黄彦杰. 鄂尔多斯盆地槐树庄区长7致密油成藏机理与富集规律[D].西安:西安石油大学,2020:1-102.

HUANG Y J. Accumulation Mechanism and Enrichment Law of Chang 7 Tight Oil in Huaishuzhuang Area,Ordos Basin[D].Xi’an:Xi’an Shiyou University,2020:1-102.

41
韩载华. 鄂尔多斯盆地七里村油田长7烃源岩与成藏动力研究[D].西安:西安石油大学,2020:1-87.

HAN Z H. Study on Source Rock and Reservoir Forming Dynamic of Chang 7 Member in Qiicun Oilfield,Ordos Basin[D]. Xi’an:Xi’an Shiyou University,2020:1-87.

42
HERMANRUD C,WENSAAS L,TEIGE G M G,et al.Shale Porosities from Well Logs on Haltenbanken(Offshore Mid-Norway) Show No Influence of Overpressuring[J].AAPG Memoir,1998,70:87-104.

43
BOWERS G L. Detecting high overpressure[J]. The Leading Edge, 2002, 21(2):174-177.

44
TINGAY M R P, HILLIS R R, SWARBRICK R E, et al. Origin of overpressure and pore-pressure prediction in the Ba-ram province,Brunei[J].AAPG Bulletin,2009,93(1):51-74.

45
侯志强,张书平,李军,等.西湖凹陷中部西斜坡地区超压成因机制[J].石油学报,2019,40(9):1059-1068,1124.

HOU Z Q,ZHANG S P,LI J,et al. Genetic mechanism of overpressures in the west slope of central Xihu sag[J]. Acta Petrolei Sinica,2019,40(9):1059-1068,1124.

46
梁晓伟,关梓轩,牛小兵,等.鄂尔多斯盆地延长组7段页岩油储层储集性特征[J].天然气地球科学,2020,31(10):1489-1500.

LIANG X W,GUAN Z X,NIU X B,et al. Reservoir characteristics of shale oil in Chang 7 Member of Yanchang Formation,Ordos Basin[J]. Natural Gas Geoscience,2020,31(10):1489-1500.

47
DASGUPTA S, CHATTERJEE R, MOHANTY S P. Magnitude, mechanisms, and prediction of abnormal pore pressure using well data in the Krishna-Godavari Basin, east coast of India[J]. AAPG Bulletin, 2016, 100(12):1833-1855.

48
O’CONNOR S, SWARBRICK R, LAHANN R. Geologically‐driven pore fluid pressure models and their implications for petroleum exploration. Introduction to thematic set[J]. Geofluids, 2011, 11(4):343-348.

49
TINGAY M R P, MORLEY C K, LAIRD A, et al. Evidence for overpressure generation by kerogen-to-gas maturation in the northern Malay Basin[J]. AAPG Bulletin, 2013, 97(4):639-672.

50
ATHY L F. Density, porosity, and compaction of sedimentary rocks[J]. AAPG Bulletin,1930, 14(1):1-24.

Outlines

/