《Petroleum》：A comparative study of estimating the oil generation window for evaluation and tracking the source rock potential by indirect methods
编辑推荐：
研究人员采用间接方法估算了生油深度与热成熟度。初始计算利用测井导出方程（log-derived equations）完成，目标测井曲线与地震属性（seismic attributes）相耦联。间接确定的参数获得了地球化学数据（geochemical data）
研究人员采用间接方法估算了生油深度与热成熟度。初始计算利用测井导出方程（log-derived equations）完成，目标测井曲线与地震属性（seismic attributes）相耦联。间接确定的参数获得了地球化学数据（geochemical data）的支持，并在整个地震剖面上对烃源岩层段进行了追踪。
**基于间接方法的生油窗估算与烃源岩潜力追踪对比研究解读**
**1. 研究背景与科学问题**
波斯湾盆地侏罗系—白垩系地层序列是全球最富集的油气省之一，其中阿尔布阶（Albian）Kazhdumi组是该区域最重要的烃源岩层段。然而，尽管盆地东部和西部已有大量研究，波斯湾盆地中部区域，尤其是Kazhdumi组，勘探程度仍极低。烃源岩评价所必需的地球化学数据严重匮乏，岩石样品的直接地球化学测试在数据连续性和可获取性方面均存在显著局限。这种数据稀缺性严重制约了该区石油系统的精细评价。传统的烃源岩评价高度依赖岩样地球化学分析，成本高昂且耗时，在海上探区尤为受限。因此，发展不依赖于直接岩石地球化学测试的替代性评价方法，实现生油窗位置与热成熟度（thermal maturity）的间接估算，对于在区域尺度上有效追踪和评价烃源岩潜力具有重要的科学意义与应用价值。本研究旨在利用测井资料与地震数据，通过间接方法在波斯湾盆地五个油气田中评价Kazhdumi组的热成熟度并确定生油窗深度，以应对地球化学数据缺失的问题。
**2. 研究内容与总体结论**
研究人员聚焦于波斯湾盆地五个油气田（A至E），首先利用声波和电阻率测井计算生油门限参数（ΔRRS），以估算生油窗起始深度；随后基于地温梯度和测井资料计算镜质体反射率（Vitrinite Reflectance, VR%）并生成连续测井曲线；继而通过地震属性分析，将井中计算的VR%外推至井间区域，首次实现了直接由地震属性预测烃源岩热成熟度。研究结果显示，研究区内生油窗起始深度介于约1350米至超过2000米之间，其横向变化与研究区主要地质构造特征（尤其是卡塔尔背斜和盐丘构造）高度一致。测井估算的热成熟度结果得到了超过200个Rock-Eval热解数据（其中80个来自Kazhdumi组）、镜质体反射率数据以及埋藏史模拟的综合验证，证实了该间接方法在缺乏地球化学数据情况下的适用性和可靠性。这一方法显著降低了烃源岩评价的时间与经济成本，并可将成熟度预测扩展到未钻井区域，为区域勘探决策提供了高效技术手段。
**3. 关键研究方法**
本研究采用了多方法耦合的技术路线。首先，基于声波测井（DT）和电阻率测井（Rt）累积曲线差异计算ΔRRS，据此确定生油窗起始深度，并依托地温梯度（27°C至34°C）通过经验公式计算VR%连续曲线。其次，采用基于模型的波阻抗（Acoustic Impedance, AI）反演方法，并结合多属性分析筛选与VR%相关性最高的地震属性（如瞬时频率、余弦瞬时相位和反演波阻抗等）。再次，利用概率神经网络（Probabilistic Neural Network, PNN）方法进行监督学习，建立地震属性与VR%之间的非线性映射模型，进而将VR%预测外推至整个地震剖面。此外，研究还运用了一维埋藏史与热成熟度模拟（分别应用Baur 2019和Sweeney & Burnham 1990模型），选取B、C、E三个油气田进行独立的埋藏史评价。最后，在E油田基于70张岩屑薄片开展了沉积环境分析，并结合自然伽马（GR）与声波测井进行岩性与总有机碳（TOC）解释。用于验证的地球化学数据集中包含超过200个Rock-Eval热解样品，其中Kazhdumi组样品80个，另有A、B、E、F四口井的镜质体反射率（Ro%）测试数据。
**4. 主要研究结果**
**4.1 生油窗起始深度估算（The onset of oil generation）**
基于测井资料的ΔRRS计算获得各油气田的生油窗起始深度。A油田的起始深度约1350米，对应Kazhdumi组，处于早期成熟阶段，该深度地层温度约65°C，Tmax为439°C。向东至B、C等油气田，生油窗起始深度增加至超过2000米，对应Surmeh组。生油窗深度的横向变化与区域地质构造背景密切相关，例如A、B油田附近的盐丘构造可导致局部热异常，从而降低生油窗起始深度。这一变化模式与现有地球化学指标（Tmax和Ro%）的分布趋势相吻合。
**4.2 基于测井与地震资料的VR%预测（VR prediction by well logs and seismic data）**
研究人员通过测井资料计算了A至E油田Kazhdumi组的热成熟度（VR%），总体呈现从欠成熟到成熟的连续变化。B、C、E油田以欠成熟特征为主，而A、D油田则处于早期成熟至成熟阶段。基于井中VR%测井曲线和地震多属性分析，研究人员建立了PNN预测模型，将VR%由已知井外推至二维地震剖面，生成了连续的热成熟度体积（VR volume）。模型验证结果显示，预测的VR%值与井中实际计算值具有较高的一致性，证明地震属性能够有效反映烃源岩的热成熟度空间分布。生成的VR剖面清晰揭示了成熟度在横向与垂向上的不均匀性，并识别出局部高热成熟度区域。
**4.3 地球化学测试结果（Geochemical tests results）**
Rock-Eval热解与镜质体反射率（Ro%）测试在B、E等油田的多套层系中完成。结果显示，Kazhdumi组在E油田具有较高的TOC（最高约4.25 wt%）和氢指数（HI，最高约674 mgHC/g TOC），显示出良好的生烃潜力。Ro%数据则整体呈现随深度增加而增大的趋势，例如B油田Surmeh组在2302米处Ro%均值为0.624%，至Dehram组4007米处升至1.252%，与测井估算的VR%变化趋势一致。
**4.4 埋藏史与热成熟度模拟（Burial history and thermal maturity modeling）**
针对B、C、E三个油气田的一维埋藏史模拟表明，生油窗起始深度的模拟结果与基于测井的间接估算值大体一致。受构造位置（距卡塔尔背斜的距离）影响，B、C油田（位于背斜上方）的深度相对较浅，而E油田的深度加大。这一结果进一步验证了间接方法的可靠性，并揭示了构造格局对生油窗空间分布的控制作用。
**4.5 沉积环境分析（Depositional environment）**
E油田的薄片及测井综合分析揭示，Kazhdumi组沉积于一个低能量、水体循环受限的较深水环境，代表该地区均斜缓坡（homoclinal ramp）的最深部位。研究识别出中缓坡、外缓坡和盆地相三种沉积亚环境。其中盆地相沉积富含远洋和浮游生物化石（如Hedbergella和Oligostoginid），且具有高TOC含量，显示E油田Kazhdumi组具有最高的烃源岩潜力，这与测井和地震属性预测所指示的高成熟度区域相吻合，进一步验证了间接评价方法的准确性。
**5. 讨论与结论总结**
讨论部分强调，本研究首次实现了利用地震属性直接预测镜质体反射率（VR%），这一成果无需进一步转换或标定即可使用，极大简化了烃源岩评价流程。测井计算获得的生油窗深度和VR%曲线与地球化学测试结果（Ro%相关系数达0.9）、Rock-Eval参数及埋藏史模拟结果均高度吻合，三重验证确保了间接方法的高可靠性。研究同时指出，生油窗深度和热成熟度的区域变化明显受卡塔尔背斜生长和盐丘侵入等构造活动的控制，构造抬升导致的差异埋深和局部热流异常是不可忽视的因素。该方法成本低、效率高，可扩展至未钻井区域，为在勘探早期或数据稀疏区优选烃源岩有利区带、减少盲目钻探提供了有力工具。
该研究的结论可综合归纳为：通过测井资料估算生油窗深度和热成熟度，并结合地震属性进行空间外推，实现了在缺乏岩样地球化学数据条件下对烃源岩热成熟度的快速、连续评价；区域上，热成熟度和生油窗深度受卡塔尔背斜及盐丘等地质构造控制，生成并标定了生油窗深度图件；测井估算的结果与地球化学数据及埋藏史模型校准良好，证明该方法可作为一种优秀的替代评价方案；首次从地震属性直接建立的VR%预测模型，可在井间区域追踪烃源岩成熟度，极大降低了烃源岩评价的时间成本与勘探风险。最后，研究指出将测井曲线耦联至地震属性，有助于通过在地震剖面上识别高热潜力区来缩小沉积学研究的范围。
