[发明专利]一种基于解析信号振幅约束的欧拉反褶积方法

说明书

本发明公开了一种基于解析信号振幅约束的重力欧拉反褶积方法，主要包括以下步骤：步骤一：解析信号振幅圈定参与计算的异常值范围；步骤二：根据异常值范围赋予合适的构造系数N值；步骤三：矩形窗口扫描构建欧拉反褶积反演矩阵求解；步骤四：统计方法剔除离散解值。本发明通过创建一种模拟数据模型，得到反演结果与模型是吻合的，并以此对比规则形体的边缘检测法与欧拉反褶积结果，显示解析信号振幅与欧拉反褶结果具有很好的吻合性，由此提供了一种圈定范围‑分布赋N值‑反演计算‑统计剔除离散值的“四步法”欧拉反褶积方法，模拟数据反演表明该方法能够有效的反演出浅深、相邻等异常体的空间分布，能够适用于复杂区域的重力数据反演计算。

技术领域

本发明涉及重力异常源的深度位置演算分析技术领域，具体涉及一种基于解析信号振幅约束的重力欧拉反褶积方法。

背景技术

重力资料在金属矿藏、油气田勘探、盆地研究、造山带构造及深部壳幔研究中都发挥着重要作用，重力数据具有很好的稳定性，横向分辨能力强，而由于不同深度层的密度不均匀性均对重力异常值产生影响，且不同深度不同规模的异常体在地表产生的异常特征有可能是相同的，故重力勘探的纵向分辨能力是较弱的，这时提高重力资料的纵向分辨能力就显得尤其重要。

因此，能够反演出重力异常源的深度位置，能够在很大程度上上很好的改善重力资料的纵向分辨能力，扩宽重力资料的勘探能力和应用广度。

发明创造内容

为解决上述问题，本发明提出一种基于解析信号振幅约束的重力欧拉反褶积方法，通过反演出重力异常源的深度位置，能够很好的改善重力资料的纵向分辨能力。

本发明所采用的技术方案如下：一种基于解析信号振幅约束的重力欧拉反褶积方法，主要包括以下步骤：

步骤一：解析信号振幅圈定参与计算的异常值范围；

步骤二：根据异常值范围赋予合适的构造系数N值；

步骤三：矩形窗口扫描构建欧拉反褶积反演矩阵求解；

步骤四：统计方法剔除离散解值。

本方案进一步的，上述步骤一中通过不断改变解析信号振幅阈值，圈定不同范围，反映不同埋深异常体计算范围。

本方案进一步的，上述步骤二中N的赋值是根据解析信号振幅圈定范围的形态，自动识别异常体形体特征进行赋值；N的赋值采用表格1所列数值。

表格1

本方案进一步的，上述步骤三中的具体方法为：

采用正方形窗口移动方法进行反演矩阵的构造，设定计算点是(x_i,y_i)，其反演矩阵为公式(1)所示：

其中：

其中Δf_(i-m,j-m),x,表示(i-m,j-m)点的重力异常对x的导数，相应的下标为y,z时，表示对y,z的导数，B为背景场，N为构造系数,2m+1表示正方形窗口边的数据量。

本方案进一步的，上述步骤四中剔除离散值的方法为通过计算均值和均方差进行剔除，具体计算方法为：

1)先按公式(2)计算每一圈定区域均值：

2)再按照公式(3)计算均方差：

3)满足公式(4)的计算值予以保留，不满足的则予以剔除:

与现有技术相比较，本发明的有益效果是：