[发明专利]一种湖相碳酸盐岩孔隙度的计算方法

权利要求书

1.一种湖相碳酸盐岩孔隙度的计算方法，其特征在于，包括以下步骤：

步骤S1：选定湖相碳酸盐岩储层，在地层框架约束下构建湖相碳酸盐岩的波阻抗模型，获取湖相碳酸盐岩的地震反演模型；具体包括以下步骤：

步骤S11：选定湖相碳酸盐岩储层，获取湖相碳酸盐岩的全叠加地震数据、测井数据、地层速度模型和层位数据，其中测井数据包括波阻抗参数；

步骤S12：根据湖相碳酸盐岩的地下层位的层序关系并利用层位数据构建湖相碳酸盐岩的地层层位框架；

步骤S13：根据测井数据开展地层速度与波阻抗参数的交会统计分析，获得湖相碳酸盐岩的地层速度与波阻抗参数的岩石物理拟合关系，将湖相碳酸盐岩的地层速度模型转化为初始波阻抗模型；

步骤S14：对湖相碳酸盐岩的井信息进行地质统计学分析，并据此建立湖相碳酸盐岩的变差函数；

步骤S15：在湖相碳酸盐岩的地层层位框架下，以所述初始波阻抗模型为协克里金约束量，结合湖相碳酸盐岩的变差函数，利用协克里金算法更新波阻抗模型获得湖相碳酸盐岩的更新波阻抗模型；

步骤S16：对湖相碳酸盐岩的更新波阻抗模型数据进行0-10Hz低通滤波处理，获得湖相碳酸盐岩的低频地震反演模型，其矩阵形式记为M；

步骤S2：通过多道空间约束地震反演方法和湖相碳酸盐岩的地震反演模型提高湖相碳酸盐岩地震反演波阻抗的空间表征能力，获取湖相碳酸盐岩的三维波阻抗参数；

步骤S3：利用BP神经网络算法实现湖相碳酸盐岩三维波阻抗参数与孔隙度参数的转化，获取湖相碳酸盐岩的孔隙度参数。

2.如权利要求1所述的湖相碳酸盐岩孔隙度的计算方法，其特征在于，所述步骤S2包括以下步骤：

步骤S21：根据全叠加地震数据和测井数据提取地震子波，获得子波矩阵W，并构建地震反演的褶积算子P，其中褶积算子的公式如下：

P＝kron(I,W) (式1)

式中，kron(I,W)表示I和W的Kroneker积，I表示单位矩阵；

步骤22：通过多道空间约束的地震反演目标方程建立三维地震数据与三维波阻抗参数的映射关系，其中多道空间约束的地震反演目标方程的表达式为：

式中，S表示全叠加地震数据的矩阵；

m为待反演的地下地层界面的反射系数向量；

λ₁为低频地震反演模型的约束系数；

M为低频地震反演模型的矩阵形式；

K＝kron(I,C)，其中t为终止时间，t₀为起始时间，τ为时间采样点；

λ₂为多道空间约束系数；

||m||_2,1为多道空间约束项；

步骤23：对多道空间约束的地震反演目标方程进行最小化求解，获得m的反演结果，并计算湖相碳酸盐岩的波阻抗数据I_P，其中波阻抗数据I_P的表达式为：

式中，t为终止时间，t₀为起始时间，τ为时间采样点。

3.如权利要求2所述的湖相碳酸盐岩孔隙度的计算方法，其特征在于，所述步骤S3包括以下步骤：

步骤S31：对所述三维波阻抗参数进行条件化处理，包括按照设计的采样率对空间进行插值和利用扩散滤波处理方法除去随机噪音；

步骤S32：对湖相碳酸盐岩的测井孔隙度数据进行条件化处理，包括0-90Hz带通滤波和空间采样粗化处理，使其空间采样率与对三维波阻抗参数进行空间插值所采用的采样率保持一致；

步骤S33：利用BP神经网络的算法进行波阻抗参数与孔隙度参数的非线性映射计算，将湖相碳酸盐岩的波阻抗参数转化为湖相碳酸盐岩的孔隙度参数。