[发明专利]一种基于重力曲率的地质体分布判断方法及处理终端

说明书

本发明涉及一种基于重力曲率的地质体分布判断方法及处理终端，所述方法包括如下步骤：步骤1：获得原始测量的重力异常数据，并将原始测量的重力异常数据进行网格化处理，得到网格化后的重力异常数据；步骤2：构建滑动窗口，对滑动窗口内的所有测点的重力异常数据建立方程组，并转换成矩阵方程；步骤3：依次移动滑动窗口，每次移动一个网格单元，建立对应的矩阵方程，遍历所有网格上的重力异常数据，利用所有矩阵方程，建立整体的矩阵方程，并进行求解，得到各个重力曲率属性的分布关系图；根据各个重力曲率属性的分布关系图，判断地下地质体包括水平位置在内的分布。本发明有效克服噪声的影响，计算出的重力曲率属性能够更准确。

技术领域

本发明涉及油气地质勘探开发地质技术领域，具体是一种基于重力曲率的地质体分布判断方法及处理终端。

背景技术

重力曲率属性包括最大正曲率、最小负曲率、平均曲率、高斯曲率、最大曲率和最小曲率，理论上，通过这些重力曲率属性可以分析地质线性构造特征，从而对地质体分布(例如水平位置分布)进行判断。但现有计算重力曲率属性的过程中，由于受到噪声的影响，导致计算的重力曲率属性不准确，进而影响对地质体分布的判断结果，导致误差较大。目前，还没有较好的噪声影响下基于重力曲率判断地质体分布方法。

相关的重力曲率属性参考文献如下：

Roberts A.Curvature attributes and their application to 3Dinterpreted horizons[J].First Break,2001,19(2):85-100.

Hansen R.O.,DeRidder E.Linear feature analysis for aeromagnetic data[J].Geophysics,71(6):L61-67.

Phillips J.D.,Hansen R.O.,Blakely R.J.The use of curvature inpotential-field interpretation[J].Exploration Geophysics,2007,38:111-119.

Cevallos C.Automatic generation of 3D geophysical models usingcurvatures derived from airborne gravity gradient data[J].Geophysics,2014,79(5):G49-G58.

Li X.Curvature of a geometric surface and curvature of gravity andmagnetic anomalies.Geophysics,80(1):G15-G26.

郭良辉，高锐，孟小红，等.重力异常线性构造分析的正负混合曲率法[J].地球物理学进展，2015，30(3)：1110-1117。

发明内容

针对现有技术的不足，本发明的目的之一提供一种基于重力曲率的地质体分布判断方法，其能够解决克服由于噪声计算曲率不准确而不能很好判断地质体分布的问题；

本发明的目的之二提供一种处理终端，其能够解决克服由于噪声计算曲率不准确而不能很好判断地质体分布的问题。

实现本发明的目的之一的技术方案为：一种基于重力曲率的地质体分布判断方法，包括如下步骤：

步骤1：获得原始测量的重力异常数据，并将原始测量的重力异常数据进行网格化处理，得到网格化后的重力异常数据，