[发明专利]一种基于正则化反演的沉积盆地基底断裂重力识别方法

说明书

本发明提供了一种基于正则化反演的沉积盆地基底断裂重力识别方法，具体包括：将实测重力异常进行网格化处理得到规则网重力数据，分别沿x和y方向从规则网重力数据中提取2.5维重力异常剖面；利用正则化密度界面反演方法对x和y方向的所有2.5维重力异常剖面进行反演；根据反演结果计算水平总梯度thgr；利用thgr的极大值位置识别盆地基底断裂的位置。本发明针对沉积盆地基底断裂特点，利用对基底断裂更为敏感的正则化密度界面反演进行断裂识别，在正则化密度界面反演方法中建立Lp范数形式的模型约束函数，使密度界面反演方法更能突出基底断裂的位置，能够快速地识别出沉积盆地基底小规模断裂。

技术领域

本发明属于盆地基底断裂位置识别技术领域，具体涉及一种基于正则化反演的沉积盆地基底断裂重力识别方法。

背景技术

断裂构造特征对盆地的形成和演化、盆地的分布范围、沉积中心等具有重要的控制作用，对盆地油气的运移和富集也具有重要影响，断裂解释也一直是沉积盆地构造和油气勘探重力研究中的一个重要内容。现有利用重力异常识别断裂的技术包括总水平导数、解析信号振幅、Tilt导数、Theta图等基本方法，以及在基本方法上发展起来的组合方法，如解析信号振幅倾斜角法、解析信号振幅总水平导数法、倾斜角总水平导数等。

为提高重力异常识别盆地基底断裂的能力，通常采用归一化、边界增强滤波等技术对基本位场边缘识别方法进行增强处理，利用不同类型的导数构建的比值类边缘识别技术可较好的均衡深部和浅部地质体的异常，也是提高深部断裂识别能力的主要措施之一。此外，相关系数利用总水平导数与垂直导数的相关系数来进行边界识别，也是识别深部断裂的常用方法。

尽管以上边缘识别技术能有效的均衡深部场源的重力异常，但用于识别沉积盆地基底断裂时，对于盆地基底小规模断裂的识别效果欠佳，其主要原因在于此类断裂引起的重力异常非常微弱，很难通过直接对重力异常处理而提高对此类断裂的识别能力。尽管可采用提高导数阶次等措施来增强对弱异常的识别能力，但往往会使得计算结果不稳定，采用平滑措施又会使得一些细节模糊，很难在提高识别精度与保证计算稳定之间找到最佳的折衷方案。可见，仅从计算的角度，对于原本非常微弱的重力异常无法处理得到清晰的盆地基底小规模断裂位置识别结果，需要转换思路，寻找一种从重力异常中导出的，包含清晰的断裂位置信息的“异常”，并对其处理得到盆地基底断裂位置。

发明内容

本发明的目的在于解决上述现有技术中存在的难题，提供一种基于正则化反演的沉积盆地基底断裂重力识别方法，利用反演技术从重力异常中导出了一种更贴近断裂位置信息的数据，并利用边缘识别提取技术对该数据进行处理，可以识别沉积盆地基底中现有方法无法识别的小规模断裂。

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

本发明提供了一种基于正则化反演的沉积盆地基底断裂重力识别方法，所述方法包括以下步骤：

步骤1：将实测重力异常进行网格化处理得到规则网重力数据，分别沿x和y方向从规则网重力数据中提取2.5维重力异常剖面；

步骤2：利用正则化密度界面反演方法对x和y方向的所有2.5维重力异常剖面进行反演；

步骤3：根据步骤2的反演结果计算水平总梯度thgr；

步骤4：利用thgr的极大值位置识别盆地基底断裂的位置。

本发明的进一步改进在于：

步骤2：利用正则化密度界面反演方法对x和y方向的所有2.5维重力异常剖面进行反演，具体包括：

步骤2.1，构建正则化密度界面反演方法的目标函数φ(m)；

步骤2.2，采用非线性共轭梯度法求解目标函数φ(m)的极小化问题。

本发明的进一步改进在于：

步骤2.1，构建正则化密度界面反演方法的目标函数，具体为：