石油地球物理勘探  2020, Vol. 55 Issue (1): 206-216  DOI: 10.13810/j.cnki.issn.1000-7210.2020.01.024
0
文章快速检索     高级检索

引用本文 

王彦国, 黄远生, 张瑾. 基于重力异常增强的改进小子域滤波方法. 石油地球物理勘探, 2020, 55(1): 206-216. DOI: 10.13810/j.cnki.issn.1000-7210.2020.01.024.
WANG Yanguo, HUANG Yuansheng, ZHANG Jin. Gravity interpretation using improved small subdomain filtering of enhanced anomaly. Oil Geophysical Prospecting, 2020, 55(1): 206-216. DOI: 10.13810/j.cnki.issn.1000-7210.2020.01.024.

本项研究受国家重点研发计划项目“铀矿基地深部成矿条件地球物理探测技术研究”(2017YFC0602603)、江西省自然科学基金项目“基于位场多源型tilt-depth法的场源深度快速自动反演研究”(20171BAB213030)、国家自然科学基金项目“基于位场广义梯度张量的欧拉反褶积方法研究”(41504098)和东华理工大学放射性地质与勘探技术国防重点学科实验室开放基金“稳定的反犙滤波方法在深部铀矿勘查中的应用”(RGET1806)联合资助

作者简介

王彦国 副教授, 博士, 1985年生。2008年获吉林大学地球物理学专业学士学位; 2010年和2013年分别获吉林大学固体地球物理学硕士和博士学位。现就职于东华理工大学地球物理与测控技术学院, 主要研究方向为位场数据处理解释及方法研究

王彦国, 江西省南昌市经开区广兰大道418号东华理工大学地球物理与测控技术学院, 330013。Email:wangyg8503@126.com

文章历史

本文于2019年5月10日收到,最终修改稿于同年11月7日收到
基于重力异常增强的改进小子域滤波方法
王彦国 , 黄远生 , 张瑾     
东华理工大学放射性地质与勘探技术国防重点学科实验室, 江西南昌 330013
摘要:小子域滤波是一种常用的重力异常边界拾取手段。传统小子域滤波子域划分重心不稳,易导致异常曲线出现扭曲。为此对小子域滤波的子域剖分模式进行改进,使其更合理地反映不同走向的构造;同时提出了基于迭代差分的稳定增强滤波解决小子域滤波易受高频干扰影响和识别异常界限能力不足的缺陷。模型试验表明,稳定增强异常改进型小子域滤波不仅计算稳定性强,而且检测出的边界与实际场源边界位置一致性较高。在鸭绿江盆地重力数据处理中,增强异常改进小子域滤波相对于常规算法可以获得更丰富的边界细节信息,而且检测出的异常界限与岩性接触带具有较好的一致性,同时计算结果中的负异常区域与研究区低密度岩石的分布也具有一定的对应关系,不仅验证了本文算法的有效性和实用性,还为研究区煤田勘查及地质—地球物理综合解释提供了可靠依据。
关键词小子域滤波    异常界限    随机干扰    增强异常    鸭绿江盆地    
Gravity interpretation using improved small subdomain filtering of enhanced anomaly
WANG Yanguo , HUANG Yuansheng , ZHANG Jin     
Fundamental Science on Radioactive Geology and Exploration Technology Laboratory, East China University of Technology, Nanchang, Jiangxi 330013, China
Abstract: Small subdomain filtering(SSF) is widely used for picking anomaly boundaries of gravity data, but the method could make anomaly curve distorted because of the unreasonable subdomain division.In this paper, the subdomain division pattern of SSF was improved to reflect the structures with diffe-rent strikes.Moreover, a stable anomaly-enhanced filtering method based on iterative differential was proposed to solve the problem that small subdomain filtering is easily affected by high-frequency disturbance and enhance anomaly boundary identification.Compared with the conventional SSF, in model tests, the improved SSF of enhanced anomaly is more stable, and the detected edges are more consistent with boundaries of field sources.In the application to the gravity data of Yalu River Basin, the improved SSF of enhanced anomaly output more detailed results than conventional SSF.The detected boundaries match lithologic contact zones well, and the negative value area can generally reflect the distribution of low-density rocks in research area.It demonstrated the effectiveness and practicability of the method proposed in this paper, and the results facilitate coalfield exploration and geological-geophysical comprehensive interpretation.
Keywords: small subdomain filtering    anomaly boundary    random disturbance    enhanced anomaly    Yalu River Basin    
0 引言

运用位场数据识别地质体边界是对地球物理资料进行地质解释的一项重要内容。近年来,在利用位场资料推算断裂构造的水平位置及深度方面已出现了许多方法。其中基于位场梯度的边界识别方法算法[1]简单、易于实现,是一类常用的方法(如水平总梯度[2-3]、垂向导数[4-5]、Theta map[6]、Tilt angle[7]、归一化标准差[8]等)。但是这类方法对区域边界信息不敏感,且易受随机干扰的影响。Euler反褶积法[9-11]也可以用于边界识别,该方法是通过计算位场一阶导数和点位坐标实现的,导数计算结果或计算点位稍有偏差,便会在反演中产生较大的偏离,从而极大地影响边界反演精度。利用重力归一化总梯度及其相位同样可以确定场源深度和边界水平位置[12], 但此法受谐波数、向下延拓的影响,场源深度和水平位置的估计受到了一定限制,因此精度不高。

小子域滤波法是杨高印[13]提出的一项位场数据处理技术,是对滑动平均法的一种改进,该方法能够增强区域异常的界限特征。夏玲燕等[14]将小子域滤波和Tilt梯度相结合,许海红等[15]将小子域滤波和水平总梯度相结合,均获得了较良好的边界识别效果。马涛等[16]针对小子域滤波存在偏心问题对梯度保持滤波法进行改进,使滤波结果更加稳定可靠;肖锋等[17]采用“田”字型子域划分模式,使计算更便捷;马国庆等[18]针对小子域滤波存在异常曲线扭曲的缺陷,提出了优化小子域滤波,更好地展示了区域异常间的界限。段晓旭[19]采用加密点距的方式解决小子域存在异常曲线扭曲的问题,并提出点位移动法增强异常梯度带。蔡钟等[20]针对小子域滤波存在边界异常失真的现象,提出了五/六边形子域滤波法,在一定程度上改善了小子域滤波的效果。许海红等[21]对上述多种小子域滤波方法进行了对比分析,给出了不同方法的优缺点及适用条件。但是上述改进方法仅能对具有明显梯级带特征的断裂构造进行识别,而无法检测具有弱梯级带或非梯级带特征的场源边界。为此,张凤旭等[22-23]提出了三方向小子域滤波进行断裂构造识别,提高了边界识别能力;Jiang等[24]对小子域滤波的剖分模式进行了改进,并实现了基于梯度张量的小子域滤波,同样获得了丰富的边界信息;Tai等[25]在三方向小子域滤波基础上,提出了加权小子域滤波法,使断裂构造的划分更快捷、更精确。但这三种改进方法均采用了差分或导数计算,因此计算更易受到干扰的影响。

针对上述方法的优缺点,本文提出了新的子域剖分模式,使其可以更准确地检测出不同方向的异常界限; 同时为了消除干扰对小子域滤波的影响和提高异常分辨率,给出了一种新的滤波技术——基于迭代差分的稳定增强滤波。

1 改进小子域滤波 1.1 传统小子域滤波

传统小子域滤波[13]是将滑动窗口按中心点的不同侧面划分出8个子域,并将8个子域均方差最小子域的平均值作为滤波输出,图 1为传统小子域滤波子域剖分模式的示意图。该方法可以较好地保留异常间的界限,但还存在着以下不足:①虽然对位场异常进行了低通滤波处理,然而不同子域对噪声的敏感程度不同,从而导致滤波结果中出现明显的“虚假”异常[26];②8个子域的重心均偏离窗口中心点,易造成异常形态不规则扭曲,从而产生虚假信息[15, 18];③不能较好地检测出梯级带不明显或非梯级带特征的异常界限,且检测结果存在较大偏差[22]

图 1 传统小子域滤波5×5滑动窗口及分解成的8个子域[13]
1.2 改进型小子域滤波

针对传统小子域滤波存在的问题,本文对子域剖分模式进行改进,改进型小子域滤波与传统方法的区别在于剖分子域的个数不同。以5×5窗口为例,传统剖分是按中心点的不同侧面划分出8个子域(图 1),而本文方法的剖分是将5×5窗口剖分出21个子域,每个子域包含9个数据点(图 2)。从图 2中可以看出,改进后的子域划分方式可以检测任意方向上的异常变化,能更好地体现出不同走向上的构造,获得更为准确、丰富的异常界限。该方法的计算流程与传统方法一致,即将窗口内21个子域均方差最小的子域数据平均值作为窗口中心点的滤波结果,以紧缩梯级带作为场源边界的识别依据。

图 2 改进型小子域滤波5×5滑动窗口及分解成的21个子域

对于N×N窗口来说,子域剖分个数为$\frac{(N+1)(N+9)}{4}$,每个子域的数据点为$\frac{(N+1)^{2}}{4}$。也就是说,随着窗口增大,子域剖分个数也随之增加,而传统小子域滤波及其他改进方法的子域剖分个数大多是固定的,与窗口大小无关。因此本文提出的子域剖分模式则更能有效地检测异常的细节变化。

2 异常增强滤波技术

传统滤噪技术一般是低通滤波。低通滤波在消除噪声时,对低频有效信号同样也有一定的消减,从而导致异常分辨率降低。因此,在原异常中表现为非梯度带的异常界限,小子域滤波仍不能有效识别。鉴于此,本文提出了一种稳定增强滤波技术,该方法在滤除噪声的同时,可以提高对异常的分辨能力。

2.1 基本原理

设位场异常u的波谱为U,低通滤波算子为φ,0≤φ≤1,滤波后的波谱可以表示为

$ U_{\varphi}=U \varphi $ (1)

虽然上述低通滤波可以滤除随机干扰,但对低频有效信息也存在一定的消减作用。因此,采用迭代补偿模式完成对低频成分的弥补,迭代过程如下:

(1) 将UUφ的差值利用滤波算子φ进行修正,修正结果为

$ U_{\varphi}^{(1)}=\left(U-U_{\varphi}\right) \varphi+U_{\varphi}=U \varphi[1+(1-\varphi)] $ (2)

(2) 将UUφ(1)的差值再次利用φ进行修正,修正结果为

$ \begin{aligned} U_{\varphi}^{(2)} &=\left(U-U_{\varphi}^{(1)}\right) \varphi+U_{\varphi}^{(1)} \\ &=U \varphi\left[1+(1-\varphi)+(1-\varphi)^{2}\right] \end{aligned} $ (3)

(3) 按照如上迭代补偿模式,第n次修正结果为

$ U_{\varphi}^{(n)}=\left(U-U_{\varphi}^{(n-1)}\right) \varphi+U_{\varphi}^{(n-1)} $ (4)

根据式(2)~式(4),得到迭代通式

$ \begin{aligned} U_{\varphi}^{(n)}=& U_{\varphi}\left[1+(1-\varphi)+(1-\varphi)^{2}+\cdots+\right.\\ &\left.(1-\varphi)^{n}\right]=U\left[1-(1-\varphi)^{n+1}\right] \end{aligned} $ (5)

上述迭代模式事实上与侯重初[27]提出的补偿圆滑滤波在原理上是一致的。由式(5)可以看出,采用迭代法得到的波谱Uφ(n)等价于在波谱U中加入了一个与迭代次数n相关的低通滤波1-(1-φ)n+1,因此迭代法仍是起到压制干扰的作用,而无法提高异常的分辨率。

为了提高异常分辨率,对迭代结果进行如下处理:(1)设初始迭代次数为n的滤波结果uφ(n)=F-1[Uφ(n)],迭代次数间隔为Δn的滤波结果uφ(nn)=F-1[Uφ(nn)],则一阶迭代差分的结果为

$ \begin{aligned} \mathrm{d} u_{\varphi}^{1 , n} &=u_{\varphi}^{(n+\Delta n)}-u_{\varphi}^{(n)}=\mathrm{F}^{-1}\left[U_{\varphi}^{(n+\Delta n)}-U_{\varphi}^{(n)}\right] \\ &=\mathrm{F}^{-1}\left\{U(1-\varphi)^{n+1}\left[1-(1-\varphi)^{\Delta n}\right]\right\} \end{aligned} $ (6)

式中F-1表示反傅里叶变换。

(2) 令迭代次数间隔为2Δn的滤波结果为uφ(n+2Δn)=F-1[Uφ(n+2Δn)],则二阶迭代差分的结果为

$ \begin{aligned} \mathrm{d} u_{\varphi}^{2, n} &=2 u_{\varphi}^{(n+\Delta n)}-\left[u_{\varphi}^{(n)}+u_{\varphi}^{(n+2 \Delta n)}\right] \\ &=\mathrm{F}^{-1}\left\{U(1-\varphi)^{n+1}\left[1-(1-\varphi)^{\Delta n}\right]^{2}\right\} \end{aligned} $ (7)

(3) 仿照上述流程,第m阶迭代差分结果为

$ \mathrm{d} u_{\varphi}^{m , n}=\mathrm{F}^{-1}\left\{U(1-\varphi)^{n+1}\left[1-(1-\varphi)^{\Delta n}\right]^{m}\right\} $ (8)

由式(8)可知,迭代差分m次的滤波因子Φ=(1-φ)n+1[1-(1-φ)Δn]m事实上是一个低通滤波ΦL=[1-(1-φ)Δn]m与高通滤波ΦH=(1-φ)n+1的组合。

2.2 滤波特性

令滤波算子φ=exp(-wh),m=1,Δn=10,其中w表示波数,h=10Δx,Δx=0.1km。图 3a给出了上文提及的滤波算子的滤波曲线。可以看出,低通滤波φ(向上延拓算子)虽然对高频成分具有较好的压制作用,但对中、低频的有用信号一定程度上也有消减作用;相对于低通滤波φ,增强滤波器中的低通滤波ΦL可以更好地保留低频成分;稳定增强滤波Φ在中、低频则趋近于高通滤波算子ΦH,而在高频处则趋于低通滤波ΦL,即滤波算子Φ为一带通滤波器,对中频成分起到了放大作用,而对高频成分起到了压制作用。也就是说,稳定增强滤波后的位场异常在理论上不仅具有较强的稳定性,而且还可以提高原异常的分辨能力,因此,将这种迭代差分滤波称为稳定增强滤波。

图 3 稳定增强滤波算子滤波响应曲线 (a)不同滤波算子的滤波响应曲线;(b)不同n时的滤波响应曲线;(c)不同Δn时的滤波响应曲线;(d)不同m时的滤波响应曲线

图 3b~图 3d是不同参数的增强滤波响应曲线。可以看出,随着初始迭代次数n的增加(图 3b),滤波算子的主峰向高频移动,但幅值却随之降低;随着迭代次数间隔Δn的增加(图 3c),增强滤波算子频带变宽,主峰同时向高频移动;随着迭代差分阶数m的增加(图 3d),滤波算子频带变窄且主峰向低频移动。因此选择较少的初始迭代次数n、较小的迭代次数间隔Δn和较大的迭代差分次数m,会使滤波结果更稳定。

理论上,稳定增强滤波器中的低通滤波ΦL、迭代初始次数n、迭代间隔Δn及迭代差分阶数m的选择具有随意性,但为确保计算结果的稳定性及异常分辨能力,建议φ选用向上延拓不少于5倍点距的延拓算子,选取n=1、Δn < 50及m < 5。后文中的模型及实例部分低通滤波算子均是选取向上延拓10倍的延拓算子,m=n=1,Δn=10。

3 理论模型分析

为了验证本文方法的有效性和优越性,建立了由4个埋深不同、尺度不同、密度不同的重力异常体(表 1)组成的复杂模型进行试验。同时,在模型的正演重力异常中添加了随机干扰,以检验方法的稳定性。图 4是该模型产生的重力异常和分别添加1%、3%噪声的重力异常及利用本文方法对图 4c进行稳定增强滤波处理的结果。可以看出,理论重力异常(图 4a)可以通过梯度带大致识别异常体A、B、C的边界位置,但梯度带较宽,边界位置不明确;模型体D的异常受叠加场影响,圈闭极值与模型体中心不对应,在异常图主要以等值线同向扭曲为主。在重力异常含1%随机噪声时(图 4b),异常体D所对应的异常更加模糊。当重力异常含3%噪声时(图 4c),难以通过异常梯度带识别所有模型边界。稳定增强滤波后的重力异常(图 4d)不仅等值线相对于图 4c更加圆滑,且异常分辨率也得到明显提高,尤其对规模小、埋深大的异常体D。

表 1 理论模型异常体的参数

图 4 叠加模型产生的理论重力异常(a)、含1%噪声的重力异常(b)、含3%噪声的重力异常(c)及其稳定增强滤波结果(d)

图 5是无噪数据在不同窗口下的传统小子域滤波与本文改进方法滤波处理结果。可以看出,当计算窗口较小时(如3×3),传统小子域(图 5a上)与改进小子域滤波(图 5a下)仅对埋深较浅的异常体C的边界异常进行了有效紧缩;随着窗口尺度的增加,无论是传统小子域滤波还是改进小子域滤波,地质体边界异常被紧缩得更加明显,如9×9窗口时,除异常体D缺失部分边界,其他三个异常体的边界异常均得到了有效增强。由图 5可以明显发现,传统小子域滤波结果在模型体角点处存在明显的等值线扭曲现象,而改进型小子域滤波的处理结果并不存在此现象,能更好地识别异常体的边界。

图 5 不同大小窗口下的传统小子域滤波(上)和改进小子滤波(下)结果对比 (a)滤波窗口3×3;(b)滤波窗口5×5;(c)滤波窗口7×7;(d)滤波窗口9×9

图 6是不同噪声水平的重力异常在不同窗口下的传统小子域滤波与改进小子域滤波处理结果对比。可以看出,由于随机噪声对不同方向子域的影响程度不同,导致了两种方法的滤波结果中均存在着多条不连续分布的异常紧缩带,且无法确定哪些紧缩带是有效的;而且,噪声水平越高,小子域滤波处理结果中的虚假信息就越多,处理结果的可靠性越差。

图 6 含1%及3%噪声时不同窗口下的传统小子域滤波(上)与本文滤波(下)结果对比 (a)含1%噪声、窗口5×5;(b)含1%噪声、窗口9×9;(c)含3%噪声、窗口5×5;(d)含3%噪声、窗口9×9

图 7是含3%噪声重力异常经过稳定增强滤波后的传统小子域滤波和改进型小子域滤波处理结果对比。相对于图 6来说,稳定增强滤波后的小子域滤波更加稳定,不但不包含明显的虚假紧缩带,更能有效地检测出地质体边界,尤其异常体D的边界也有明显的展示。

图 7 含3%噪声重力异常经稳定增强滤波后的传统小子域滤波和改进小子域滤波结果对比 (a)传统小子域滤波,窗口5×5;(b)传统小子域滤波,窗口9×9;(c)改进小子域滤波,窗口5×5;(d)改进小子域滤波,窗口9×9
4 实际重力资料应用

为了验证稳定增强异常改进型小子域滤波对实际资料的处理效果,选取吉林省南部鸭绿江盆地实测重力数据(点距和线距均为1km)进行试验。鸭绿江盆地位于中朝板块东北缘,二级大地构造单元隶属于辽东台隆区,盆地主体为太子河—浑江坳陷,盆地内除志留系、泥盆系地层缺失外,其他时代的岩石均有出露(图 8)。研究区存在一系列逆冲推覆构造,使得地层关系十分复杂、连续性极差及地层倒置现象非常普遍[28],这给地质体边界的确定甚至地质解释都带来了较大难度。研究区内各个时代的岩石出露较全,奥陶系与石炭系间存在明显的密度差(表 2),这为方法试验的可靠性提供了较好的物性基础。

图 8 鸭绿江盆地地质图及工区位置

表 2 鸭绿江盆地地层/岩石密度统计表 单位:g·cm-3

从盆地内布格重力异常图(图 9a)可以看出,反映地质体边界位置或断裂构造的梯级带异常梯度较平缓,因而地质体边界不易直接根据重力异常进行精确解释。较小型地质体的边界受区域异常的影响较大,在异常图中主要表现为异常等值线突然变宽或变窄以及同形扭曲等非梯级带特征,地质体的边界位置确定难度较大。从图 9b可以看出,经稳定增强处理的重力异常图的分辨率获得了明显提高,不仅异常梯级带连续性增强,且异常呈现更明显的NE走向,与地质图(图 8)中的岩石分布呈现出了很强的一致性。图 10是布格重力异常的传统小子域和改进型小子域滤波结果,可以看出,对于原始布格重力异常,传统小子域滤波(图 10上)对大型异常梯度带进行了较好的紧缩,然而紧缩带的连续性较差,且存在明显的无规律性弯曲;而改进型小子域滤波(图 10下)同样也反映了大型断裂的位置,且边界异常信息更丰富,紧缩带的连续性也更强。对于稳定增强异常,无论是经传统小子域滤波(图 11上)还是改进型小子域滤波(图 11下)处理的重力数据对边界的识别能力显然都获得了大幅度的提升,边界异常信息更丰富,而改进型小子域滤波结果显示出更丰富的构造边界信息,且紧缩带的走向和连续性也更好。需要指出的是,5×5窗口的小子域滤波(图 11a)数据中的一些紧缩带在9×9窗口的小子域滤波中力图(图 11b)上没有显示。如在大安镇的东南侧,5×5窗口的小子域滤波数据中存在两条NE走向的异常紧缩带,但9×9窗口的小子域滤波数据中并未被有效识别出来,且检测出的紧缩带连续性也较差。这说明在实际资料处理中,需要结合不同窗口的小子域滤波数据进行地质解释。另外,新方法获得的异常紧缩带与不同时代的岩性界限(图 8)具有良好的对应性;增强异常的改进小子域滤波数据的负异常分布也较好地反映了研究区中—新生代地层和燕山期花岗岩等低密度岩石的分布。

图 9 鸭绿江盆地布格重力异常(a)及稳定增强重力异常(b)

图 10 鸭绿江盆地布格重力异常的传统小子域(上)和改进型小子域(下)滤波结果对比 (a)滤波窗口5×5;(b)滤波窗口9×9

图 11 鸭绿江盆地稳定增强重力异常的传统小子域(上)和改进型小子域(下)滤波结果对比 (a)滤波窗口5×5;(b)滤波窗口9×9

图 12a是结合区域地质及前人工作成果绘制的构造分区图。图中给出了大型断裂的位置、低密度岩石的分布区(燕山期花岗岩及石炭—白垩纪地层)以及浑江煤田工作区等信息,以此验证本文方法的可靠性。图 12b是基于图 11下解释的3条大型断裂和17条中小型断裂以及在浑江坳陷内圈定的负异常区域,其中F1为二道江—江源断裂,是龙岗隆起与浑江坳陷的界限断裂;F2为石人—大青沟断裂,是浑江坳陷与老岭隆起的分界断裂;F3为鸭绿江断裂,在研究区东部则是老岭隆起与四道沟坳陷的界限。这三条大型断裂在重力异常中均以大型梯级带为异常特征,其余17条中、小型断裂在重力异常图上表现为弱异常梯级带或等值线突然变宽或变窄以及同形扭曲等异常特征。需要指出的是,在浑江坳陷内,地表出露的低密度岩石分布区仅占30%左右的面积,而本文方法识别出的负异常分布范围可达50%以上,也就是说浑江坳陷内有大面积的低密度体隐伏于高密度岩石之下,其中勘探程度较高的浑江煤田工作区中的一些飞来峰下的煤层(主要分布在石炭、二叠、侏罗纪地层中)已被钻探或采掘所证实[28-29],这也再一次证实了本文方法的实用性,同时该方法也为研究区后续煤田勘查有利区(尤其勘探程度较低的六道沟—三道湖—石人镇一带)和地质—地球物理综合解释的有效手段之一。

图 12 鸭绿江盆地构造分区图(a)和断裂构造分布及浑江坳陷负异常分布区(b)
5 结论

针对传统小子域滤波存在异常曲线扭曲的缺陷,本文提出了一种包含更多子域的新子域剖分模式,可以更准确地检测不同走向构造的界限。针对小子域滤波存在易受干扰影响和难以识别非梯级带特征的断裂构造的缺陷,提出了基于迭代差分的稳定增强滤波技术。该技术在滤除噪声干扰的同时,还可以提高异常的分辨率。模型试验和实例应用表明,改进型小子域滤波比传统算法效果更佳,异常的界限可以得到更加准确、精细的刻画;稳定增强滤波后的重力异常不仅稳定性强,且异常分辨率明显提高;增强异常改进型小子域滤波则能够更加准确、精细地刻画密度异常体的界限,同时还可以通过重力正、负异常的分布进一步提高方法的解释能力。

参考文献
[1]
张恒磊, 刘天佑, 杨宇山. 各向异性标准化方差计算重磁源边界[J]. 地球物理学报, 2011, 54(7): 1921-1927.
ZHANG Henglei, LIU Tianyou, YANG Yushan. Calculation of gravity and magnetic source boundary based on anisotropy normalized variance[J]. Chinese Journal of Geophysics, 2011, 54(7): 1921-1927. DOI:10.3969/j.issn.0001-5733.2011.07.026
[2]
Cordell L, Woodward L A, and James H L, et al.Gravimetric expression of graben faulting in the Espanola Basin[C].New Mexico Geological Society 30th Annual Fall Field Conference Guidebook, Santa Fe Country, 1979, 59-64.
[3]
王万银. 位场总水平导数极值空间变化规律研究[J]. 地球物理学报, 2010, 53(9): 2257-2270.
WANG Wanyin. Spatial variation law of the extreme value positions of total horizontal derivative for potential field data[J]. Chinese Journal of Geophysics, 2010, 53(9): 2257-2270. DOI:10.3969/j.issn.0001-5733.2010.09.027
[4]
Bhattachayya B K. Two-dimensional harmonic analysis as a tool for magnetic interpretation[J]. Geophy-sics, 1965, 30(5): 829-857.
[5]
Wang W Y, Zhang G C, Liang J S. Spatial variation law of vertical derivative zero points potential field data[J]. Applied Geophysics, 2010, 7(3): 197-209. DOI:10.1007/s11770-010-0255-z
[6]
Wijns C, Perez C, Kowalczyk P. Theta map:edge detection in magnetic data[J]. Geophysics, 2005, 70(4): L39-L43. DOI:10.1190/1.1988184
[7]
Miller H G, Singh V. Potential tilt-A new concept for location of potential field sources[J]. Journal of Applied Geophysics, 1994, 32(2-3): 213-217. DOI:10.1016/0926-9851(94)90022-1
[8]
Cooper G R J, Cowan D R. Edge enhancement of potential-field data using normalized statistics[J]. Geophysics, 2008, 73(3): H1-H4. DOI:10.1190/1.2837309
[9]
Thompson D T. EULDPH:A new technique for ma-king computer assisted depth estimates from magnetic data[J]. Geophysics, 1982, 47(1): 31-37. DOI:10.1190/1.1441278
[10]
Reid A B, Allsop J M, Granser H, et al. Magnetic interpretation in three dimensions using Euler deconvolution[J]. Geophysics, 1990, 55(1): 80-91. DOI:10.1190/1.1442774
[11]
Hansen R O, Laura S. Multiple-source Euler deconvolution[J]. Geophysics, 2002, 67(2): 525-535. DOI:10.1190/1.1468613
[12]
张凤琴, 张凤旭, 刘财, 等. 利用重力归一化总梯度及相位法研究断裂构造[J]. 吉林大学学报(地球科学版), 2005, 35(1): 123-127.
ZHANG Fengqin, ZHANG Fengxu, LIU Cai, et al. Study of faulted structure using the normalized full gradient of gravity anomaly and its phase[J]. Journal of Jilin University(Earth Science Edition), 2005, 35(1): 123-127.
[13]
杨高印. 位场数据处理的一项新技术-小子域滤波法[J]. 石油地球物理勘探, 1995, 30(2): 240-244.
YANG Gaoyin. A new technique for potential-field data processing:small subdomain filtering[J]. Oil Geophysical Prospecting, 1995, 30(2): 240-244.
[14]
夏玲燕, 吴汉宁, 柏冠军, 等. 柴达木盆地航磁资料微弱信息增强技术研究及在线性构造识别中的应用[J]. 地球物理学进展, 2008, 23(4): 1058-1062.
XIA Lingyan, WU Hanning, BO Guanjun, et al. Research on enhancing weak signal technology and re-cognition of linear structures using aerial-magnetic data in the Qaidam Basin[J]. Progress in Geophysics, 2008, 23(4): 1058-1062.
[15]
许海红, 袁炳强, 杨高印, 等. 小子域滤波联合总水平导数在重力资料处理中的应用——以内蒙古银额盆地EQ区块断裂构造识别为例[J]. 地质通报, 2014, 33(11): 1853-1860.
XU Haihong, YUAN Bingqiang, YANG Gaoyin, et al. The application of the small sub-domain filtering combined with the total horizontal derivative to gravity data processing:a case study of fault structure identification in EQ block of Yin-E Basin, Inner Mongolia[J]. Geological Bulletin of China, 2014, 33(11): 1853-1860. DOI:10.3969/j.issn.1671-2552.2014.11.024
[16]
马涛, 王铁成, 王雨. 一种改进的网格数据保持梯度滤波方法[J]. 石油地球物理勘探, 2007, 42(2): 198-201.
MA Tao, WANG Tiecheng, WANG Yu. An improved gradient preserved grid data filtering approach[J]. Oil Geophysical Prospecting, 2007, 42(2): 198-201. DOI:10.3321/j.issn:1000-7210.2007.02.015
[17]
肖锋, 吴燕冈, 孟令顺. 位场数据处理中小子域滤波的改进[J]. 石油地球物理勘探, 2010, 45(1): 136-139.
XIAO Feng, WU Yangang, MENG Lingshun. An improvement on small subdomain filtering in potential field data processing[J]. Oil Geophysical Prospecting, 2010, 45(1): 136-139.
[18]
马国庆, 杜晓娟, 李丽丽. 优化小子域滤波方法研究及其应用[J]. 石油地球物理勘探, 2013, 48(4): 658-662.
MA Guoqing, DU Xiaojuan, LI Lili. An optimized small sub-domain filtering method[J]. Oil Geophysical Prospecting, 2013, 48(4): 658662.
[19]
段晓旭.重磁异常小子域滤波算法的分析研究与改进[D].北京: 中国地质大学(北京), 2014.
DUAN Xiaoxu.Research and Improvement of Small Subdomain Filtering of Gravity and Magnetic Anomaly[D].China University of Geosciences(Beijing), Beijing, 2014. http://cdmd.cnki.com.cn/Article/CDMD-11415-1014233762.htm
[20]
蔡钟, 倪小东. 重力异常处理中小子域滤波法的改进[J]. 河北工程大学学报(自然科学版), 2015, 32(4): 47-51.
CAI Zhong, NI Xiaodong. An improvement on small sub-domain filtering method in processing gravity anomalies[J]. Journal of Hebei University of Enginee-ring(Natural Science Edition), 2015, 32(4): 47-51. DOI:10.3969/j.issn.1673-9469.2015.04.011
[21]
许海红, 李玉宏, 袁炳强, 等. 位场数据小子域滤波法处理效果对比与优选[J]. 地质通报, 2018, 37(1): 153-164.
XU Haihong, LI Yuhong, YUAN Bingqiang, et al. The processing effect comparison and optimization of small subdomain filtering method in potential field data[J]. Geological Bulletin of China, 2018, 37(1): 153-164.
[22]
张凤旭, 张凤琴, 刘财, 等. 断裂构造精细解释技术-三方向小子域滤波[J]. 地球物理学报, 2007, 50(5): 1543-1550.
ZHANG Fengxu, ZHANG Fengqin, LIU Cai, et al. A technique for elaborate explanation of faulted structures:three-directional small subdomain filtering[J]. Chinese Journal of Geophysics, 2007, 50(5): 1543-1550. DOI:10.3321/j.issn:0001-5733.2007.05.031
[23]
张凤旭, 张兴洲, 张凤琴, 等. 中国东北地区重力场研究-利用改进的三方向小子域滤波划分主构造线及大地构造单元[J]. 地球物理学报, 2010, 53(6): 1475-1485.
ZHANG Fengxu, ZHANG Xingzhou, ZHANG Fengqin, et al. Study on gravity field in Northeastern China area:Classification of main structure lines and tectonic units using the improved three-directional small subdomain filtering[J]. Chinese Journal of Geophysics, 2010, 53(6): 1475-1485.
[24]
Jiang F Y, Gao L K. Edge enhancement of gravity anomalies and gravity gradient tensors using an improved small sub-domain filtering method[J]. Applied Geophysics, 2012, 9(2): 119-130. DOI:10.1007/s11770-012-0321-9
[25]
Tai Z H, Zhang F X, Zhang F Q, et al. Weighted small subdomain filtering technology[J]. Journal of Applied Geophysics, 2017, 144: 50-56. DOI:10.1016/j.jappgeo.2017.06.012
[26]
王彦国, 张凤旭, 刘财, 等. 位场垂向梯度最佳自比值的边界检测技术[J]. 地球物理学报, 2013, 56(7): 2463-2472.
WANG Yanguo, ZHANG Fengxu, LIU Cai, et al. Edge detection in potential fields using optimal auto-ratio of vertical gradient[J]. Chinese Journal of Geophysics, 2013, 56(7): 2463-2472.
[27]
侯重初. 补偿圆滑滤波法[J]. 石油物探, 1981, 20(2): 22-29.
HOU Zhongchu. Filtering of smooth compensation[J]. Geophysical Prospecting for Petroleum, 1981, 20(2): 22-29.
[28]
殷长建, 周晓东, 陈跃军. 吉南浑江坳陷南缘重力滑覆构造特征[J]. 吉林地质, 1999, 18(3): 28-33.
YIN Changjian, ZHOU Xiaodong, CHEN Yuejun. The feature of gravitation gliding capping structure of the southern edge of depression in the Hunjiang area, Southern Jilin Province[J]. Jilin Geology, 1999, 18(3): 28-33.
[29]
何保.浑江煤田推覆构造特征、演化及找煤远景区预测[D].辽宁沈阳, 东北大学, 2009.
HE Bao.The Feature and Evolution of Nappe Structure and Coal Prospect Area Predication of Hunjiang Coalfield[D].Northeastern University, Shenyang, Liaoning, 2009. http://cdmd.cnki.com.cn/article/cdmd-10145-2010055514.htm