[发明专利]深度域储层地震反演方法及装置

说明书

技术领域

本发明涉及油气勘探技术领域，尤其涉及深度域储层地震反演方法及装置。

背景技术

在石油天然气工业中，地震反演作为储层预测的核心技术，能够将测井纵向上的高分辨率优势和地震横向上的密集采样优势紧密结合起来，估算出地层横向上的岩性、物性等特征信息，目前已经被广泛应用于石油勘探开发的各个阶段。但是常规的时间域地震反演方法对于薄储层的识别能力较差，其测井曲线从深度域换算到时间域的过程中丢失了很多宝贵的高层信息，存在很大的局限性。

近年来，对于地震反演储层预测的最为常见的做法是基于褶积模型的深度域地震资料反演。该方法的首要问题是子波的获取，即通过时深转换将时间域子波转换为深度域子波，然后根据地震道褶积和反子波，消除子波的褶积效应，从而获得一个近似的地层反射系数系列，以便用于反演高分辨率的岩层型剖面。

发明人在实现本发明的过程中，发现上述现有技术存在如下不足：该方法虽然可以保留测井资料的高频信息，但是在操作过程中难以建立进行时深转换所需要的精确三维速度场和提取深度域地震子波，且在地震子波确定时，相位谱的确定容易出现误差，不适宜大规模应用。

发明内容

本发明实施例提供一种深度域储层地震反演方法，用以克服测井曲线在时深转换过程中高频信息易丢失的缺点，且避免了在反演过程中求取深度域子波，包括：

根据时间域叠后地震数据，获得深度域叠后地震数据；

将测井曲线和深度域叠后地震数据进行标定，获得深度域测井曲线模型；

对深度域叠后地震数据进行地层层位解释，并将解释结果数据网格化，获得深度域地质框架模型；

分别计算单井的深度域测井曲线模型和深度域叠后地震数据的分形参数，将每一个地震道的分形参数转换为测井曲线的分形参数体；

根据深度域测井曲线模型和深度域地质框架模型，对测井曲线的分形参数体进行分形随机插值，获得深度域储层地震反演结果。

一个实施例中，所述根据时间域叠后地震数据，获得深度域叠后地震数据，包括：

建立空变速度场，根据空变速度场将时间域叠后地震数据转换为深度域叠后地震数据。

一个实施例中，所述深度域测井曲线模型包括：波阻抗测井曲线模型、自然伽马测井曲线模型和孔隙度测井曲线模型其中之一或任意组合。

一个实施例中，所述根据深度域测井曲线模型和深度域地质框架模型，对测井曲线的分形参数体进行分形随机插值，获得深度域储层地震反演结果，根据如下公式进行计算：

其中，S(a)是插值点a的深度域储层地震反演结果，f_line是分形随机插值的线性分量，f_unine是分形随机插值的随机分量，K是标定系数，H是反应深度域叠后地震数据细节和粗糙度的分维参数，σ是深度域叠后地震数据正态分布的方差，G是服从标准正态分布的Gauss随机变量，ΔX是已知点和插值点之间的距离。

一个实施例中，所述分形随机插值的线性分量f_line采用方位距离加权法实现，其中，

在研究区域内，网格节点z(x,y)的估计值为：

其中，D是以网格节点z(x,y)为原点构成的四个象限内总的等分角数，l、i、j、k是已知点的序号，n是已知点的总数，λ_i是在D范围内第i个已知点的权系数，z_i是在D范围内第i个已知点的参数值，r_j、r_l分别是第j个及第l个等分角内离待估计点最近已知点的距离。

一个实施例中，深度域储层地震反演方法还包括：

在获得深度域储层地震反演结果后，对深度域储层地震反演结果进行精度检验。

一个实施例中，所述在获得深度域储层反演结果后，对深度域储层地震反演结果进行精度检验，包括：

分别抽取全部反演井的每一口井作为验证井，对未被抽取的反演井进行深度域储层地震反演，将未被抽取的反演井的反演结果与验证井的测井曲线进行对比，计算出未被抽取的反演井的反演结果和验证井的测井曲线之间的反演误差。

本发明实施例中还提供了一种深度域储层地震反演装置，包括：

时深转换模块，用于根据时间域叠后地震数据，获得深度域叠后地震数据；

井震联合标定模块，用于将测井曲线和深度域叠后地震数据进行标定，获得深度域测井曲线模型；