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

摘要

一种利用位场高阶水平梯度模识别地质体边界的方法，原理是首先构建圆滑去噪线性方程组，求取圆滑去噪滤波算子，将该算子与重磁数据褶积，获得去噪之后的重磁数据；构建水平导数线性方程组，求取高阶水平导数模滤波算子，将该算子与去噪之后的重磁数据褶积，获得重磁数据的高阶水平导数，然后求取高阶水平导数模；绘制高阶水平导数模等值线图或映像图，根据高阶水平导数模的极大值识别地质体边界或断裂的水平位置。

主权项

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)绘制高阶水平导数模等值线图或映像图，根据高阶水平导数模的极大值识别地质体边界或断裂的水平位置。