[发明专利]一种重磁干扰区的方差系数干扰校正方法

说明书

技术领域

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

背景技术

重磁勘探方法是两种重要的地球物理勘探方法，在油气及金属、非金属矿产勘探中应用广泛。地表人文干扰对地面重力磁力勘探造成资料的畸变影响，需要对这种干扰畸变进行校正才能获得地下地质体的不受干扰的真实资料数据。现有重力、磁力勘探规范中(《地面重力勘探技术规程(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=0,|dz|≤1.25σ(|dz|σ-1.25)2×0.05,1.25σ<|dz|<3σ0.15+(|dz|σ-3.0)1.5×0.106,3σ≤|dz|<6.0σ0.70,|dz|≥6.0σ]]>