[发明专利]一种地球物理勘探中岩性约束下的孔隙度反演方法

权利要求书

1.一种地球物理勘探中岩性约束下的孔隙度反演方法，其特征在于：所述方法利用地震数据和测井数据，构建孔隙度反演的目标函数；

所述测井数据包括孔隙度、伽马曲线和波阻抗曲线；

所述方法包括：

(1)波阻抗反演；

(2)岩性反演；

(3)在反演的岩性数据的约束下，对所述目标函数进行反演得到储层的孔隙度,

所述步骤(3)包括：

(31)采用步骤(2)得到的岩性数据体对地震数据进行过滤，只保留砂岩部分的数据，然后对保留的数据进行高斯配置协模拟，其公式如下：

$y={\Σ}_{i=1}^n{w}_i{x}_i+{\Σ}_{j=1}^m{u}_j{z}_j---\left(4\right)$

公式(4)就是所述孔隙度反演的目标函数，

其中，y为待估算的孔隙度；

w为i点处的权重,u为j点处的权重，由变差函数求得；

x为井点处孔隙度曲线值，由测井的孔隙度曲线求得；

z为步骤(1)得到的三维波阻抗数据体；

(32)采用模拟退火的算法，对所述目标函数进行反演，得到砂岩中的孔隙度，

其中，岩性数据体是将步骤(1)得到的三维波阻抗数据体按照不同的门槛值重新编码成l或0而得到的。

2.根据权利要求1所述的地球物理勘探中岩性约束下的孔隙度反演方法，其特征在于：所述步骤(1)包括：

(11)利用测井曲线求取岩性曲线；

(12)分岩性对井中的波阻抗数据进行地质统计分析：应用步骤(11)得到的岩性曲线，分别对砂岩和泥岩进行波阻抗直方统计分析得到砂岩和泥岩的阻抗分布规律，然后分别对砂岩和泥岩进行变差函数分析，最后采用线性拟合的方法对砂岩孔隙度曲线和波阻抗曲线的相关关系进行拟合，得到孔隙度与阻抗的关系曲线；

(13)在测井资料和地质层位的约束下，应用步骤(12)得到的砂岩和泥岩的阻抗分布规律对地震数据进行高精度测井约束波阻抗反演，获得三维波阻抗数据体。

3.根据权利要求2所述的地球物理勘探中岩性约束下的孔隙度反演方法，其特征在于：所述步骤(12)中的变差函数分析是采用下面的变差函数实现的：

变差函数是指区域化变量Z(x)在x和(x+h)两点处的增量的半方差：

G(x,h)＝1/2*∑[Z(x)–Z(x+h)]² (1)。

4.根据权利要求3所述的地球物理勘探中岩性约束下的孔隙度反演方法，其特征在于：所述步骤(2)是采用序贯指示模拟实现的，所述序贯指示模拟如下：

将步骤(1)得到的三维波阻抗数据体按照不同的门槛值重新编码成l或0，得到岩性数据体：

假设在位置u处的随机函数Z(u)，对门槛值为z的指示变换可写成：

则

Prob{Z(u)＜＝z|(n)}＝E{I(u，z)|(n)} (3)

对于模拟目标区内的每一类岩性，当它出现在某一位置时，指示变量为l，否则为0。