一种重磁干扰区的方差系数干扰校正方法

摘要

地球物理勘探中重磁干扰区的方差系数干扰校正方法根据重力数据或磁力三列数据计算阀门值，阀门值等于6倍的数据采集误差，计算第一个测点周围3倍采集点距以内的观测数据值的平均值和方差，计算测点观测值与平均值之差，比较与阀门值大小，绝对值大于进行校正，对其余所有观测数据点操作一遍，获得一次迭代后的数据，得到校正后的数据值。本发明校正后重磁异常变得光滑连续、异常变化符合地下场源的异常特征，反映了地下地质体的分布特征。

说明书

技术领域

本发明属于地球物理勘探方法，是一种重磁干扰区的方差系数干扰校正 方法。

背景技术

重磁勘探方法是两种重要的地球物理勘探方法，在油气及金属、非金属 矿产勘探中应用广泛。地表人文干扰对地面重力磁力勘探造成资料的畸变影 响，需要对这种干扰畸变进行校正才能获得地下地质体的不受干扰的真实资 料数据。现有重力、磁力勘探规范中(《地面重力勘探技术规程 (SY/T5819-2002)》；《地面磁法勘探技术规程(SY/T5771-2004)》)，针对干 扰畸变影响提出的校正方法是“去除奇异点”，即删除该单点畸变观测值。而 对于多点、连片的干扰区则没有提出对应的校正方法，实际重磁资料处理中 也没有好的方法，往往无处下手，使得干扰畸变得不到校正，使得干扰区域 乃至更大区域的重磁资料解释失真或错误。

发明内容

本发明目的在于提供一种可提高干扰区重磁资料可利用性和解释的可靠 性，消除重力、磁力勘探中对于多点、连片干扰区畸变影响的重磁干扰区的 方差系数干扰校正方法。

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

1)野外采集重力数据或磁力三列数据，并得到数据采集点距D、野外数 据采集观测误差ε；

步骤1)所述的三列数据是x坐标、y坐标、重力值或磁力值z。

步骤1)所述的点距D的单位为m；观测误差ε单位为重力毫伽或磁力 nT。

2)计算阀门值f，阀门值f等于6倍的数据采集误差；

3)计算第一个测点周围3倍采集点距以内的观测数据值z(0)的平均值 A和方差σ，计算测点观测值z(0)与平均值A之差dz，

dz＝z(0)-A；

4)比较dz的绝对值与阀门值f大小，如果dz的绝对值小于f，则不进 行校正；如果dz的绝对值大于f，则要进行校正，校正后的新数据值z(1)，

z(1)＝z(0)-(1-k)*dz；

式中：A为该统计窗口内数据的平均值，dz为计算测点观测值z(0)与 平均值之差，k为校正系数；

步骤4)所述的校正系数k的取值根据dz的绝对值与方差σ的倍数关系 而定，为0-70％。

步骤4)所述的校正系数k的取值根据以下公式确定：

$k=\left(\begin{array}{c}0,\left|\mathrm{dz}\right|\le 1.25\sigma \\ {(\frac{\left|\mathrm{dz}\right|}{\sigma }-1.25)}^2\times \mathrm{0.05,1.25}\sigma <\left|\mathrm{dz}\right|<3\sigma \\ 0.15+{(\frac{\left|\mathrm{dz}\right|}{\sigma }-3.0)}^{1.5}\times \mathrm{0.106,3}\sigma \le \left|\mathrm{dz}\right|<6.0\sigma \\ 0.70,\left|\mathrm{dz}\right|\ge 6.0\sigma \end{array}\right)$

5)重复步骤3)和步骤4)对其余所有观测数据点操作一遍，获得一次 迭代后的数据x，y，z(1)；

6)重复步骤3)至步骤5)5次，得到校正后的数据值。

本发明使得干扰区域的重磁异常畸变得到校正(见图1)，校正后重磁异 常变得光滑连续、异常变化符合地下场源的异常特征(连续有序、无突跳)， 反映了地下地质体的分布特征。

附图说明

图1为本发明干扰校正前磁异常图；

图2为本发明干扰校正后磁异常图。

具体实施方式

以下结合附图详细说明本发明。

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

1)野外采集重力数据或磁力三列数据，并得到数据采集点距、野外数据 采集观测误差；

所述的三列数据是x坐标、y坐标、重力值或磁力值z。

如已知数据采集点距20m、采集观测误差2.0nT。实施如下：

2)计算阀门值f，阀门值f＝12.0nT。

3)计算第一个测点(4068620，761340，50.0)周围60m以内的观测数 据值z(0)的平均值A和方差σ，计算测点观测值z(0)与平均值A之差dz， z(0)＝50.0nT，A＝65.0nT，σ＝4.0nT，dz＝-15.0nT。

4)比较dz的绝对值与阀门值f大小关系，dz的绝对值15＞f(f＝12)， 则校正后的第一个测点的新数据值z(1)＝50-(1.0-0.215)×(-15)＝61.8。 校正系数k的取值根据dz的绝对值与方差σ的倍数关系而定(3.75σ)，校正 系数k取0.215，k的取值采用以下公式计算：

$k=\left(\begin{array}{c}0,\left|\mathrm{dz}\right|\le 1.25\sigma \\ {(\frac{\left|\mathrm{dz}\right|}{\sigma }-1.25)}^2\times \mathrm{0.05,1.25}\sigma <\left|\mathrm{dz}\right|<3\sigma \\ 0.15+{(\frac{\left|\mathrm{dz}\right|}{\sigma }-3.0)}^{1.5}\times \mathrm{0.106,3}\sigma \le \left|\mathrm{dz}\right|<6.0\sigma \\ 0.70,\left|\mathrm{dz}\right|\ge 6.0\sigma \end{array}\right)$

5)对其余所有观测数据点操作一遍，获得一次迭代后的数据(x，y，z (1))：统计计算点周围3倍采集点距以内的观测数据值z(0)的平均值A和 方差σ，计算测点观测值z(0)与平均值之差dz；比较dz的绝对值与阀门 值f大小关系，如果dz的绝对值小于f，则不进行校正；如果dz的绝对值大 于f，则要进行校正，校正后的新数据值z(1)＝z(0)-(1-k)×dz。

6)依照3)、4)、5)步的方法进行重复迭代操作，迭代至5次。比如第 一个测点5次迭代后的值为(4068620，761340，61.8)

7)完成5次迭代，5次迭代后的(x，y，z(5))即为校正后的数据值， 校正后异常与校正前异常的对比图见图1。