[发明专利]一种自动提取地质体埋藏深度的方法

权利要求书

1.一种自动提取地质体埋藏深度的方法，其基于地球物理异常能谱自动提取地质体埋藏深度的方法，其特征在于，其包括以下步骤：

第一步、原始数据获取：首先利用重力仪或磁力仪野外采集地球物理数据；

第二步、对上述获得地球物理数据进行预处理，得到径向平均对数能谱；

第三步、通过径向平均对数能谱不同波段进行线性拟合计算，和拟合结果的自动优化选择；

第四步、依据第三步线性拟合的斜率信息，实现地质体埋藏深度的自动提取，斜率的负值即为地质体埋藏深度。

2.如权利要求1所述的自动提取地质体埋藏深度的方法，其中第三步包括以下步骤：

1)输入径向平均对数能谱及数据区参考地质体埋深；

2)根据参考地质体埋深计算拟合波段参考；在实际拟合时，可对该参考区适当扩大；

3)在拟合波段内有N个数据采样点，自第一个数据采样点开始，连续选取M个数据进行函数拟合并给出拟合误差，其中M大于等于函数拟合所需最小数据量，对线性函数拟合，M≥3；N、M为正整数；

4)按步骤3)，依次对拟合波段内以第二个数据点及以后的数据点作为拟合波段数据起点，遍历拟合波段内所有可能用于拟合的连续数据，并给出拟合结果及拟合误差函数值；

5)选取步骤3)、4)拟合函数误差最小值对应结果作为最优的地质体埋深计算结果，并给出拟合使用的相应波段信息。

3.如权利要求2所述的自动提取地质体埋藏深度的方法，其特征在于，拟合波段对应于研究区参考地质体埋深对应的主要频谱波段。

4.如权利要求2所述的自动提取地质体埋藏深度的方法，其特征在于，平均对数能谱拟合计算根据地质体特征具有线性或分形拟合形式。

5.如权利要求2所述的自动提取地质体埋藏深度的方法，其特征在于，拟合优化计算的误差函数e具有以下形式：e＝σ(f(k)-d)/(k₂-k₁)，其中f(k)为拟合函数，d为实际径向平均对数能谱，k1，k2是参与拟合运算所取的波数范围，σ(f(k)-d)是在[k1，k2]波数段内拟合数据与径向平均对数能谱数据之差的标准方差。

6.如权利要求2-4所述的自动提取地质体埋藏深度的方法，其特征在于，拟合计算根据拟合形式可自动选取拟合波段内的所有可能的连续采样点进行计算。

7.如权利要求6所述的自动提取地质体埋藏深度的方法，其特征在于，所有可能的连续采样点任意位于拟合波段内满足拟合计算最小采样点数量以上的连续数据。

8.如权利要求2或5所述的自动提取地质体埋藏深度的方法，其特征在于，自动提取的地质体埋藏深度对应拟合优化计算的误差函数值最小。