惠民商河油田沙三上亚段薄储层预测方法研究

乔军芳 王志海

摘 要：针对商河油田沙三上亚段的薄互层砂岩储层难度大的问题，本文采用测井约束反演结合地震属性的方法，有效的解决了难题。利用地震属性方法刻画砂体的边界，利用地质统计学反演预测砂体纵向展布特征，然后利用刻画的砂体边界对反演结果进行相控，完成砂体的预测，该方法对薄储层的预测具有较好的效果。

关键词：地震属性;属性融合;地質统计学反演;薄储层预测

商河油田地理位置位于山东省商河县境内，构造上位于惠民凹陷的临南洼陷的北部缓坡带。沙三上亚段为主力勘探开发层系，工区东部的C44、C42块已探明开发，C142断阶带中部为空白区，图1为该区沙三上亚段油气分布图。近期C142-41井沙三上亚段薄层试油效果较好，展现该断阶内沙三上亚段仍具备一定勘探潜力。沙三上亚段为浊积砂岩体沉积，受沉积控制，浊积砂体平面上连片分布、横向不连通，纵向上相互叠置，多数井钻遇的单砂体只有几米厚，制约本区沙三段岩性油藏的进一步勘探，及平面上油水关系矛盾。本文通过采用地震属性分析和测井约束反演两种方法相结合，描述该区沙三上亚段储层，解决了油水关系矛盾，取得很好的成果。

1 储层预测流程与方法

储层预测主要通过描述两个方面来表征储层的分布特征，一是砂体纵向上的分布特征（包括砂体厚度和纵向叠置形态）;二是砂体的横向连续性（描述砂体可能的边界）。由于砂体较薄，反演结果对砂体边界的刻画效果较差，而地震属性往往效果较好。本次首先运用地震属性分析方法预测砂体在平面上的展布特征及刻画砂体的边界，再进行测井约束反演刻画砂体纵向的分布特征，计算砂体厚度，最后利用地震属性分析的结果对反演的结果进行相控，得到最后的砂体厚度预测结果。

2 储层预测

2.1 地震属性分析

地震属性是储层预测的常规手段之一，是由叠前或叠后地震数据，经过数学变换而导出的有关地震波的几何学、运动学、动力学或统计学特征。依据地震波的运动学和动力学等特征将地震属性分为：振幅、频率、相位、能量、波形、衰减、相关和比率等八大类，每一大类包含几至二十几类不等。地震属性是表征地震波形态、运动学特征、动力学特征和统计特征的物理量，有着明确的物理意义。地震属性分析主要应用：特殊地质异常体（包括断裂、缝洞、河道等）的识别、岩性的识别、地层层序识别（包括沉积相带、地层不整合等）及含油气性的识别。本文采用多属性融合数，预测砂体展布特征，利用地震相和瞬时相位技术预测砂体边界。图3为S31砂组1+2+3小层融合结果，在C142-41附近存在一厚度中心，C142-7井附近则砂体不发育，在两口井之间发育一相对独立厚度中心。

2.2 地震反演

地震波阻抗反演是目前运用最广的储层预测技术。常用反演方法有稀疏脉冲反演、随机反演、模型反演。稀疏脉冲反演结果与地震资料的振幅、频率、相位等特征具有较好的对应关系，稳定性好，对测井资料数量无明确要求，适用于大多数情况，不足之处是分辨率低，无法分辨薄储层。随机反演对储层分辨能力较高，明显高于稀疏脉冲等反演方法，尤其是同时分辨薄储层和厚储层方面具有明显的优点，具有精细储层描述的能力，反演结果与测井数据吻合程度很高，并与实际地震反射资料及构造因素有较好的相关性，但要求井数据较多，且分布较均匀，井少时统计规律不强，容易产生假象。

2.2.1 反演处理

本区钻井较多，井控程度高，单砂体较薄，适用随机反演，所以本次采用了分辨率较高地质统计学随机反演。本次随机反演采用Jason的StatMC 模块进行，主要流程包括3个步骤：数据准备、随机反演函数分析、随机反演。数据准备包括前期完成的井震标定、提取子波、模型建立、稀疏脉冲反演的结果，还需要计算岩性曲线、井曲线重新采样（与模型采样一致）。函数分析主要包括概率密度函数（PDF）和变差函数的分析，在分析完两类函数之后，完成了随机反演的10个实现，并对5个实现统计计算，最后得到平均波阻抗数据体。最终反演结果的过井剖面反演结果较好展示了该区沉积特征，中间厚两侧变薄，反演结果与井吻合较好。随机反演结果分辨率相对较高，细节方面刻画的更清楚。

2.2.2 储层反演可靠性分析

本次反演过程中加强了对反演过程的质量控制，波阻抗反演剖面很好的反映了岩性变化特征。①反演结果较地震有较高的分辨率，在纵向单砂组可以区分，横向上砂体明显不连续;②反演数据体与伽玛曲线、声波曲线对应较好，砂岩阻抗值明显高于围岩阻抗值，差值明显，展布形态与该地区的地质沉积特征吻合，反演结果的精度得到提高，反演结果较真实的反映了砂体的发育和展布情况;③反演完成后，将标定结果投影到反演剖面上，井波阻抗与反演剖面吻合很好，钻遇砂体反应明显，达到埋藏深度和厚度误差较小的目的。

3 储层预测结果

3.1 砂体分布预测及厚度计算

反演厚度预测图（反演）

反演厚度预测图（相控）

本次依据反演的结果来绘制砂体厚度图，然后横向的连续性则依据之前的属性分析结果进行相控，完成砂体厚度图的绘制。图4为反演得到砂体厚度图，图5为经过井控和相控得到砂体厚度图。从结果看，C142-7与C142-41井之间存在一砂体发育区，且与东部C142-41砂体不连通，解决了平面上油水关系矛盾的问题。

3.2 目标评价

根据预测结果，在C142-7井和C142-41井之间的砂体1上倾部位部署了一口评价井，该井钻遇沙三上砂体9m/4层，油层4.5m/2层，初产8t/d，效果较好，图6为目标井位反演剖面图。

综上所述，本区沙三上亚段主要发育薄互层的浊积砂岩体，仅仅通过测井约束反演难以有效的刻画砂体的展布特中，在反演的基础上，利用地震属性刻画的砂体边界进行相控处理，非常有效，也得到了钻井的验证。

参考文献：

[1]杨国权.地震属性融合技术在河道砂体识别中的应用[J].物探化探计算技术，2016（6）：833.

[2]张林科.基于地震属性的地震相分析思路[J].工程地球物理学报，2010（6）：694.

[3]段新意.相控地质统计学反演方法研究及其应用[J].中国石油学会2019年物探技术研讨会.

作者简介：

乔军芳（1971- ），女，汉族，山东东营人，工程师，主要从事地震综合研究工作。