[发明专利]一种利用位场高阶水平梯度模识别地质体边界的方法

权利要求书

1.一种利用位场高阶水平梯度模识别地质体边界的方法，特点是通过以下技术方案实现：

1)在工区采集重力或磁力资料；

2)对野外采集的重磁资料进行各种常规校正及改正；

3)采用插值方法将重磁资料网格化，形成规则网重磁资料f₀；

4)构建圆滑去噪线性方程组，方程组为：

Σj=-LLΣi=-KKc1+c2i2+c3j2+c4i4+c5i2j2+c6j4=Σj=-LLΣi=-KK1]]>

Σj=-LLΣi=-KKc1i2+c2i4+c3i2j2+c4i6+c5i4j2+c6i2j4=Σj=-LLΣi=-KKi2]]>

Σj=-LLΣi=-KKc1j2+c2i2j2+c3j4+c4i4j2+c5i2j4+c6j6=Σj=-LLΣi=-KKj2]]>

Σj=-LLΣi=-KKc1i4+c2i6+c3i4j2+c4i8+c5i6j2+c6i4j4=Σj=-LLΣi=-KKi4]]>

Σj=-LLΣi=-KKc1i2j2+c2i4j2+c3i2j4+c4i6j2+c5i4j4+c6i2j6=Σj=-LLΣi=-KKi2j2]]>

Σj=-LLΣi=-KKc1j4+c2i2j4+c3j6+c4i4j4+c5i2j6+c6j8=Σj=-LLΣi=-KKj4]]>

式中，c₁～c₆为待求的方程解，K＝(x方向上的圆滑去噪算子点数-1)/2，L＝(y方向上的圆滑去噪算子点数-1)/2；

5)求解线性方程组，得到圆滑去噪滤波算子c₁；

6)对f₀进行圆滑去噪，即用圆滑去噪滤波算子c₁与f₀褶积运算，得到去噪之后的重磁数据f₁：

f₁＝c₁*f₀

7)构建水平导数线性方程组，方程组为：

Σj=-NNΣi=-MMa1i2+a2ij+a3i4+a4i3j+a5i2j2+a6ij3=Σj=-NNΣi=-MMi]]>

Σj=-NNΣi=-MMa1ij+a2j2+a3i3j+a4i2j2+a5ij3+a6j4=Σj=-NNΣi=-MMj]]>

Σj=-NNΣi=-MMa1i4+a2i3j+a3i6+a4i5j+a5i4j2+a6i3j3=Σj=-NNΣi=-MMi3]]>

Σj=-NNΣi=-MMa1i3j+a2i2j2+a3i5j+a4i4j2+a5i3j3+a6i2j4=Σj=-NNΣi=-MMi2j]]>

Σj=-NNΣi=-MMa1i2j2+a2ij3+a3i4j2+a4i3j3+a5i2j4+a6ij5=Σj=-NNΣi=-MMij2]]>

Σj=-NNΣi=-MMa1ij3+a2j4+a3i3j3+a4i2j4+a5ij5+a6j6=Σj=-NNΣi=-MMj3]]>

式中，a₁～a₆为待求的方程解，M＝(x方向上的水平导数算子点数-1)/2，N＝(y方向上的水平导数算子点数-1)/2；

8)求解线性方程组，得到a₃、a₆，x及y方向的高阶水平导数算子分别为：fx(3)=a3Δx3,]]>fy(3)=a6Δy3,]]>

式中：Δx为x方向的间距，Δy为y方向的间距；

9)分别对f₁进行x方向及y方向高阶水平导数计算，即进行褶积运算：

f2=fx(3)*f1]]>

f3=fy(3)*f1]]>

10)求取高阶水平梯度模：

11)绘制高阶水平导数模等值线图或映像图，根据高阶水平导数模的极大值识别地质体边界或断裂的水平位置。