[发明专利]一种地球物理勘探中的地震反演方法

说明书

本发明公开了一种地球物理勘探中地震反演方法，获取地震工区内的测井数据，对每个地震层相应点的力和力矩数据进行采集和编辑；将每层地震数据进行层层解释，并将解释结果数据网格化，获得基于深度的地质框架模型；根据确定的地质框架模型，确定测井曲线，并对测井曲线的分形参数体进行分形随机插值，获得初始测井曲线；根据地震资料确定实际测井曲线，并根据比较结果，对初始测井曲线进行修正，根据修正后的曲线确定地震反演模型。本发明通过获取地震工区内的测井数据，对每个地震层相应点的力和力矩数据进行采集和编辑，引入力、力矩各个参量，结合具体的土质环境建立勘探模型，并通过对最终值的修正，确定最终的反演模型。

技术领域

本发明涉及地震反演技术领域，具体领域为一种地球物理勘探中地震反演方法。

背景技术

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

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

测井约束地震反演技术在岩性识别和储层预测方面发挥了重要作用，通过将空间上采样密集的地震资料和垂向上高分辨率的测井资料联合起来，将界面反射型的地震资料反演出反映地层的波阻抗数据体，提高储层的识别精度，并能将储层的厚度变化情况更精细的刻画出来，尤其是对薄储层的识别方面，通过测井约束提高分辨率，能够识别常规地震不能识别的薄层。

目前的反演方法主要是基于振幅的反演方法，即认为地震振幅强弱的变化反映了储层的空间变化，并且大多都是基于模型的反演，模型的建立过程主要是在地震解释层位的约束下，测井属性的空间内差和外推，整个建模的过程，没有考虑地质体在横向上的变化情况，因此，虽然基于模型的反演方法能够提高对薄层的识别能力，但对储层的空间尖灭点识别会存在误差。

因此，现有技术中的地震反演方法依据模型以及波形识别进行反演，并未考虑其实际的物理因素及受力因素。

发明内容

本发明的目的在于提供一种地球物理勘探中地震反演方法，以解决模型中未考虑应力因素的技术问题。

为实现上述目的，本发明提供如下技术方案：包括：

步骤a，获取地震工区内的测井数据，对每个地震层相应点的力和力矩数据进行采集和编辑；

步骤b，将每层地震数据进行层层解释，并将解释结果数据网格化，获得基于深度的地质框架模型；

步骤c，根据确定的地质框架模型，确定测井曲线，并对测井曲线的分形参数体进行分形随机插值，获得初始测井曲线；

步骤d，根据地震资料确定实际测井曲线，并根据比较结果，对初始测井曲线进行修正，根据修正后的曲线确定地震反演模型；

在上述步骤a中，每个地震层的若干采样点按照下述模型建立：其中一采样点建立模型为：

其中，为该地震层采样点的检测力和力矩零偏值，为采样点处自身的重力和重力矩，和为采样点处所受到的外力和外力矩，其中，力矩为采样点相对于地震层空洞的中心线的力矩；

设该地震层空洞的中心坐标系原点到采样点质心的矢量为则

式(1)在基于地震层空洞的中心的测量坐标系中的详细表达式为：