[发明专利]一种检测储层边界的曲率属性方法

说明书

技术领域

本发明属于油田勘探领域，具体涉及一种检测储层边界的曲率属性方法，针对实际油气田勘探中存在的储层边界问题，利用高精度叠后地震资料主测线(联络测线)的振幅导数结果计算二次曲率属性。

背景技术

曲率属性作为地震属性中的一种，已经为研究碳酸盐岩地区缝洞型储层密集发区带及岩溶储层发育的最有效方法。传统的曲率计算方法多用于描述复杂储层的连续性以及延伸性。但在很多有效区带内圈定的曲率属性带有成片的模糊特征，不能精确描述缝洞型储层的边界特征。也就是说，一般的曲率属性都是描述断裂或者裂缝储层的局部突变点，展现的是储层在整个三维空间的内部连续性，而很少被用于来描述储层的边界特征。

发明内容

本发明的目的在于解决上述现有技术中存在的难题，提供一种检测储层边界的曲率属性方法，在高精度叠后地震资料纵横测线的振幅导数的基础上，开展振幅导数曲率的计算方法，可以发掘叠后地震资料中的曲率属性的边界特征，挖掘出比传统的曲率算法更为准确的储层边界特征，因此可以获得更好储层预测效果。

本发明是通过以下技术方案实现的：

一种检测储层边界的曲率属性方法，所述方法包括以下步骤：

(1)对叠后地震资料进行常规处理，得到叠后地震数据体的主测线与联络测线排列图；

(2)在所述叠后地震数据体的主测线与联络测线排列图中，以一个地震数据点为中心圆点，获取该中心圆点的沿层振幅；

(3)计算该中心圆点的导数值；

(4)计算该中心圆点的振幅曲率；

(5)依次对叠后地震数据体中的每一个地震数据点重复步骤(2)至步骤(4)，得到所有地震数据点的振幅曲率，所有地震数据点的振幅曲率即构成用于储层边界检测的曲率属性体。

所述步骤(2)是这样实现的：

在叠后地震数据体的主测线与联络测线排列图中，假定该中心圆点的主测线的序号为I，联络测线的序号为J，则该中心圆点的振幅值为A0；

中心圆点的右边相邻主测线上的点，即主测线的序号为I，联络测线的序号为J+1的点的振幅值为A1，主测线的序号为I，联络测线的序号为J+2的振幅值为A2；

中心点下方沿着联络测线方向的点，即主测线的序号为I+1，联络测线的序号为J的点的振幅值为A3，主测线的序号为I+2，联络测线的序号为J的点的振幅值为A4；(这5个振幅值都是从叠后地震数据中直接读出来的)

从叠后地震数据体中读取线间距和道间距；所述线间距是两条相邻主测线之间的间距；道间距是两条相邻联络测线之间的间距；假设线间距和道间距均为M米。

所述步骤(3)是这样实现的：

用a表示中心圆点的振幅沿主测线方向变化的二阶导数，b表示中心圆点的振幅沿联络测线方向变化的二阶导数，c表示中心圆点的振幅沿主测线、联络测线变化的偏导数，d表示中心圆点的振幅沿主测线方向变化的一阶导数，e表示中心圆点的振幅沿联络测线变化的一阶导数，f和g为求取二阶导数过程中使用的中间结果，计算公式如下：

d=A1-A0M]]>

e=A3-A0M]]>

f=A2-A1M]]>