文|大大怪

编辑|大大怪

前言

寻求地球物理响应获得的物理量与地下储层参数的对应关系通常是非常复杂的，只能通过间接途径获得，随着研究介质的复杂性变高，许多问题的求解变得繁琐而不易进行。

研究区在扇三角洲沉积背景下发育大量薄互层砂体，不仅出现大量条带状河道砂体，而且由于研究区距离物源较远，在物源匮乏区或匮乏期发育滨浅湖沉积，出现一定数量的滨浅湖砂坝透镜状砂体，导致研究区砂体叠置样式复杂，进一步导致研究区油藏复杂程度增加，储层预测难度也进一步增加。

单纯的从实际地震资料入手，在目标油区很难总结出有效的振幅、波形规律以预测相对高渗条带，因而考虑采用地震正演方法对薄砂体及砂体叠置开展地震响应研究，为更精准的预测砂体提供依据。

基于机制模型的地震正演方法是研究薄砂层反射波特征的重要手段之一，地震正演模拟主要包括物理模拟和数值模拟两大类:

物理模拟是对与实际地质体成比例的模型进行探测，模拟野外采集，通过处理获得其地震响应; 数值模拟是基于弹性介质波场传播理论在计算机环境下模拟实际地质模型的地球物理响应方法。

虽然在数值模拟地震正演方面的研究较多，但多数研究以地层、断层或单层砂体为对象，对其位置、厚度、形态变化等的地震响应开展正演分析以及基于沉积和储层发育特征的砂体叠置样式的地震响应研究较少。

因此笔者依据机制模型的地震正演结果，总结和梳理地震响应反射特征，筛选出对砂体厚度和分布比较敏感的地震属性，对砂体叠置样式等机制模型的地震响应反射特征进行研究。

研究区地质概况

苏北盆地位于江苏省中北部，内部划分为 4 个呈近东西向展布的二级构造单元，由南向北分别为东台坳陷、建湖隆起、盐阜坳陷和滨海隆起。

高邮凹陷位于苏北盆地南部东台坳陷中部，南为通扬隆起，北接柘垛低凸起，东南靠吴堡低凸起与溱潼凹陷相连，西接菱塘桥低凸起与金湖凹陷相隔，东西长度约为 100 km，南北宽度约为 25 30 km，面积为 2 670km2 ，呈北东向长条形分布，新生界地层沉积厚度为7 000 m，是苏北盆地沉降最深的一个凹陷。高邮凹陷构造单元可划分为南断阶、深凹带和北斜坡 3 个部分。

高邮凹陷深凹带是江苏油田最主要的含油气单元之一，具有新生代沉积厚度大，生储盖组合优异，油气富集等特点。戴南组是高邮凹陷深凹带主要含油气层段之一，沉积于吴堡运动造成的不整合面之上，与下伏阜四段不整合( 假整合) 接触，从下而上分别为戴一段( E2 d1 ) 和戴二段( E2 d2 ) ，主要发育扇三角洲、滨浅湖沉积。

机制模型构建原则

高邮凹陷深凹带戴南组沉积时期，同生断裂发育，构造活动剧烈，因此主要发育扇三角洲沉积，但在马家嘴油田西部等地区由于断裂活动较弱、波浪作用较强，形成滨浅湖沉积，并且在马家嘴油田出现滨浅湖与扇三角洲在纵向上互层现象，形成砂体叠置样式的多样性。因此将复杂多样的砂体叠置样式进行分类总结，从而指导机制模型的构建。

机制模型构建合理与否直接决定在理想状态下砂体叠置分布模式的地球物理响应特征是否具有代表性，更进一步决定实际油田实际资料由地震相向沉积相、岩相的解释准确与否。

首先从测井参数特征总结入手，通过各口井测井数据的统计分析得到泥岩背景及储层的声学参数，从而确定机制模型中各个单元的声学参数值设定; 其次通过测井参数统计，得到储层条带纵向上的厚度范围及条数，为机制模型如何构建及构建的单元厚度及条数奠定基础。

在确定储层、泥岩的地球物理及几何参数后进行机制模型的构建，机制模型构建遵循由等效模型向叠置模式过渡的原则。通过等效模型的建立及地震正演响应分析，总结储层、泥岩厚度、泥岩隔层、韵律变化等的地震响应差异及规律，因此需要建立以下几类模型:

( 1) 砂体楔形模型。首先需明确实际地震频率下砂体的调谐厚度，剔除因调谐厚度造成的地震振幅的影响。然后由于不同类型沉积其韵律特征不一致，落实不同韵律的砂体其正演响应特征是否一致，因此需要建立不同韵律的砂体模型与均值模型对比，分析其地震响应差异，此模型的正演响应分析结果决定进行叠置模式机制模型正演时是否考虑韵律。

( 2) 实时动态正演模型。通过改变砂体的相对位置，分析砂体连续性/不连续性的二维机制模型的响应特征。细分为砂体横向叠置正演模式，分析横向叠置时砂体边界的识别精度; 横纵向多期砂体叠置正演模型，分析砂体叠置不同宽度与不同高程差叠置下的地震响应特征及砂体识别精度。

( 3) 砂体叠置样式等效机制模型。基于砂体的主控沉积因素差异造成的砂体叠置样式的不同，建立等效的机制模型，分析砂泥厚度变化的地震响应特征，筛选表征砂泥变化的敏感地震属性。

( 4) 在多种正演模型分析基础上，基于实际井的砂体分布特征建立机制模型，分析模型正演地震响应与实际地震的对应关系。

地球物理参数评价

基于测井解释结果，分别获取砂泥岩的测井参数，统计储层与泥岩的波阻抗参数，为模型正演提供正演设置参数，力求与地下实际情况的地球物理参数接近，通过大量数据统计，砂岩与泥岩的岩石物理数据存在明显区别，因而可以开展模型正演研究。

并统计研究区单层砂体厚度解释结果，分布区间主要为 1 4 m，大于 8 m 的较少，因此建立机制模型时单层砂体小于 10 m。

地震正演等效机制模型

根据储层预测与砂体边界刻画等需求，分别建立多种类型的机制模型进行正演，包括砂体韵律性、砂体调谐厚度、砂体叠置、隔夹层厚度等多种模型，从而总结对应的地震响应特征规律。

砂体韵律模型：基于油田不同沉积类型储层的韵律性不同，分别建立均值模型和韵律模型进行机制模型正演，分析不同模型的地震响应差异。

从沉积上分析，当单砂体的沉积韵律不同时单砂体内物性也会有所差异，砂体内部韵律性是否产生明显的地震响应以及是否对砂体预测产生明显影响等问题对储层预测精度具有重要意义。因而针对单砂体的韵律性，分别建立均值模型和韵律模型进行地震正演，分析不同韵律模型的地震响应差异。

通过不同韵律模型的正演对比可知，在调谐厚度之内，韵律模型与均值模型正演的地震波形、振幅变化规律基本一致。砂体较薄时，韵律模型正演结果与均质模型正演结果差异小，当砂体逐渐变厚，虽然正演合成记录波形纵向上的位置会随着砂体变厚而逐渐变大，但是从砂体厚度－地震能量相关图上分析，二者的能量变化规律趋势完全一致。

因而可以看出韵律变化对薄层砂体地震响应的影响较小，并且在后期机制模型构建中，可以忽略韵律对地震响应的影响。因此建立机制正演模型时采用均值模型来表征基于沉积模式的地球物理响应特征。

砂体调谐厚度模型：在进行机制模型地震正演前，需要结合实际地震资料的子波频率、不同储层单元的声学参数对地震可识别的最小厚度( λ/4) 进行判别。建立泥岩背景下厚度从 1 80 m 的楔形体模型，由于实际地震资料的主频约为 20 Hz，为了便于对比分析，设置子波从 12 35 Hz，分别进行正演，并提取 ＲMS 振幅属性，以判别 20 Hz 时地震数据的识别尺度。

利用不同频率时正演的地震信号分析，提取 ＲMS 振幅，并通过判断不同频率时 ＲMS 振幅的拐点，得到不同频率时地震信号识别砂体的调谐厚度范围，例如当 20Hz 时，地震可识别砂体的厚度为 40 m。

通过研究区各储层单元的几何参数统计，砂组的地层厚度范围平均为 40 m，单砂体的厚度小于 10 m，均小于调谐厚度，尤其是含油砂体，远低于调谐厚度，因此在实际资料的基础上直接预测砂体相对困难。

所以考虑以多层相对集中发育的砂体组成的砂组为目标开展储层预测，但砂组中砂体存在多种多样的叠置方式，可能对地震响应产生较大影响，因而开展砂体叠置方式机制模型正演研究。

砂体叠置模型：由于研究区戴南组发育的砂体叠置方式多样，砂体厚度、隔夹层厚度以及联通关系复杂，因而采用动态实时正演的方法开展研究。动态实时正演是基于抽象及理论化沉积模型的正演，根据沉积模式，形成符合沉积规律的沉积模式体( 二维、三维) ，利用基于地质单元体的动态正演方法可以实时的根据沉积模式对地质体进行修改式“拖动”。

对应的地震响应也发生实时变化，可以更加直观的观察和分析地震响应随地质模式的变化而变化的规律特征。研究区主要储层类型为水下分流河道砂体和滨浅湖滩坝砂体，因而实时动态正演过程中采用砂体透镜状分布的地质模型，主要模拟河道和滩坝砂体叠置的地震响应特征。

水下分流河道地质剖面上为顶平底凸、短轴状、河道砂中间厚两端减薄; 滩坝砂体为顶凸底平状，中间厚两端减薄。因此以研究区沉积砂体发育特征为基础，建立不同参数不同叠合模式的河道砂体的动态正演模型。动态正演过程中可以灵活改变模拟参数及叠置状态，便于分析不同叠置模式砂体的响应特征，确定识别模式，分析敏感属性，指导储层预测研究。

结语

( 1) 通过对研究区测井参数等的统计，确定储层、泥岩的地球物理及几何参数，其中单砂体厚度小于 10 m，砂岩声波时差为 257 μs·m－1 、波阻抗为 9844 m·s －1 ·g·cm－3 ，泥岩声波时差为 290 μs·m－1 、波阻抗为 8 828 m·s －1 ·g·cm－3 。

( 2) 韵律模型正演分析显示韵律变化对薄层砂体地震响应的影响较小，可以忽略韵律对地震响应的影响。砂体调谐厚度模型正演分析显示当研究区主频为 20 Hz 时，地震可识别砂体的厚度为 40 m，单砂体的厚度小于 10 m，远低于调谐厚度，在实际资料的基础上直接预测砂体是非常困难。

( 3) 地震属性对于横向叠置范围小于 50 m 砂体，以及大于 1 /4 砂体厚度高程的叠置特征能够进行分辨。并且泥岩隔层在泥砂岩薄互层的岩性油藏开发中具有重要意义，因而通过砂岩隔层组合等效模型进行正演，可以建立隔层与砂体厚度变化的 4类模型。

( 4) 峰谷为多套砂泥复合的地震响应，过井剖面的正演与单井结果一致，建议薄互层砂泥岩油藏研究尺度为以砂体组合模式开展储层预测研究。

参考文献：

[1] 赵晨光．薄互层地震反射波的特征分析［J］．石油地球物理勘探，1986，21( 1) : 32-46．

[2] 潘兴祥．复杂地质体的正演模拟［J］．油气地球物理，2003，1( 3) : 53-55．

[3] 高丽坤，林春明，姚玉来，等．苏北盆地高邮凹陷古近系戴南组沉积相及沉积演化［J］． 沉 积 学 报，2010，28( 4) : 706-716．

[4] 周健，林春明，张霞，等．江苏高邮凹陷古近系戴南组一段物源体系和沉积相［J］．古地理学报，2011，13( 2) :161-174．

[5] 刘金华，吴立峰，吴艳梅，等．高邮凹陷黄珏—马家嘴地区古近系戴南组透镜状砂体成因分析［J］．天然气地球科学，2015，26( 2) ，261-268．