约束稀疏脉冲反演技术在阿格巴达组b1砂层组中的应用.docxVIP

约束稀疏脉冲反演技术在阿格巴达组b1砂层组中的应用 目前，印尼的石油勘探和开发主要集中在尼罗亚三角洲和贝宁盆地的一些地区。研究区位于尼日尔三角洲盆地的核心部位, 地处尼日利亚西南部的红树林沼泽区, 面积约280km2, 距港口城市Warri约40km。研究的目的层为海陆过渡相的阿格巴达组, 形成于新生代的始新世-上新世, 为一套进积的近岸-三角洲前缘沉积, 其沉积厚度大, 储层发育, 几乎储存了三角洲所有的油气资源。 目前研究区最新的地震资料是1995年采集的, 2006年对地震资料进行了重新处理。三维地震工区面积为275km2, 地震道采样率为2ms。三维地震基本覆盖全区, 且地震资料反射波组特征比较明显, 层位易追踪, 资料信噪比及分辨率较高, 基本可以满足储层地震反演的要求。研究区处于西非沼泽区和河流发育带, 地理情况复杂, 目前工区仅有钻井6口, 都是在20世纪60～70年代完钻的, 工区尚处于勘探开发早期。笔者运用Jason软件, 采用约束稀疏脉冲反演技术, 对研究区重点储层阿格巴达组B1砂层组进行了波阻抗反演, 分析了其砂体展布特征, 为储层预测指明了方向。 1 基于约束稀疏脉冲反演的地震声波阻抗实时监控 基于地震资料直接转换的稀疏脉冲反演方法比较完整地保留了地震反射的基本特征, 不存在基于模型反演方法多解性的问题, 能较明显地反映岩相、岩性的空间变化, 在岩性相对稳定的条件下, 能较好地反映储集层的物性变化。约束稀疏脉冲反演方法具有较宽的应用领域, 尤其在勘探初期钻井资料缺少的条件下同样适用。结合钻井资料, 波阻抗反演结果可以确定储集层岩性、物性差异界面, 分析出研究区目的层的砂体展布特征, 为储层预测和评价提供依据。 约束稀疏脉冲反演 (CSSI) 技术采用一个快速的趋势约束脉冲反演算法, 用地震解释层位和测井约束控制波阻抗的趋势和幅值范围, 脉冲算法产生宽带结果, 恢复缺失的部分低频和高频成分, 主要目的是利用约束井资料及地震反射系数建立声波阻抗数据体。它采用的数学模型, 即最小误差函数 (目标函数) 为: J=Σ|ri|p+λqΣ (di-si)q+α2Σ (ti-zi)2(1) 式中,J为目标函数;ri为反射系数样点;λ为实际地震与合成地震记录残差权重因子;di为地震道采样点;si为合成地震道采样点;α为趋势权重因子, 一般取其值为1;ti为根据测井资料定义的波阻抗约束趋势采样点;zi为介入测井约束的最小和最大波阻抗之间的采样;p,q为标准优化因子, 通常取p=1,q=2;i为反演的地震道采样序号。 2 地震层序界面和波阻抗反演技术 工区完钻的井较少, 而且测井项目不一, 资料不齐全。其中AX-1井完钻时间最晚, 其测井资料相对较为完整, 有自然电位、自然伽马、声波时差和密度测井等测井系列, 且有较全的测录井井段。 前期以高分辨率层序地层学为指导, 完成了工区三维地震层序界面的追踪与闭合工作, 对尼日尔三角洲盆地重点储层阿格巴达组地层进行了层序划分对比, 识别出B1、C1、D1、E1、E3、G、H总计7个层序界面。在此基础上, 笔者基于Jason软件平台采用约束稀疏脉冲反演技术, 充分结合地质、地震、测井资料, 对重点层位B1砂层组开展了波阻抗反演研究, 通过建立的波阻抗模型结合测井, 沉积等资料, 分析该层位的砂体展布特征。 该反演方法主要包括以下几个步骤:①测井资料预处理;②层位追踪;③初始地质模型的建立;④迭代计算;⑤反演结果输出。其中, 子波估算与合成地震记录, 初始模型建立, 约束稀疏脉冲反演处理是其反演的关键技术。 2.1 叠后地震正演估算法 地震子波是约束波阻抗反演中的关键因素。子波与模型反射系数褶积产生合成地震数据, 合成地震数据与实际地震资料的误差最小是终止迭代次数的约束条件。根据工区的实际资料, 叠后地震子波估算采用相对实用有效的方法即多道地震统计法, 通过多道记录自相关统计的方法提取子波振幅谱信息, 进而求取零相位、最小相位或常相位子波。得到的子波、合成记录与实际记录频带一致, 与实际地震记录波组关系对应关系良好, 有利于进行储层层位的精细标定。 2.2 建立波阻抗模型 建立波阻抗模型的过程实际上就是把地震界面信息与测井波阻抗结合起来的过程, 对地震而言, 是经正确层序解释起控制作用的波阻抗界面;对测井而言, 是为波阻抗界面间的地层赋予合适的波阻抗信息。建立尽可能接近实际地层情况的波阻抗模型, 是减少其最终结果多解性的根本途径。 初始模型的横向分辨率取决于地震层位解释的精细程度, 纵向分辨率受地震采样率的限制, 为了较多地保留测井的高频信息, 反映薄层的变化细节, 通常要对地震数据进行加密采样。研究区共有6口钻井, 全区均匀分布, 前期三维地震工区划分出7个层序界面, 除一口井缺少