重力异常曲化平的改进插值—迭代法

doi:10.11720/wtyht.2022.1403

重力异常曲化平的改进插值—迭代法

杨婧^,, 郭良辉^,

中国地质大学(北京) 地球物理与信息技术学院,北京 100083

A modified interpolation-iteration method for gravity anomaly continuation from undulating surface to plane

YANG Jing^,, GUO Liang-Hui^,

School of Geophysics and Information Technology, China University of Geosciences,Beijing 100083,China

通讯作者: 郭良辉(1980-),男,教授,博导,主要从事地球物理数据处理反演新方法及应用研究工作。Email:guo_lianghui@163.com

责任编辑: 王萌

收稿日期: 2021-07-22 修回日期: 2021-10-18

基金资助:

国家重点研发计划项目“月球内部圈层结构与演化过程的研究”(2021YFA0715101)
国家自然科学基金面上项目“全波形与重力联合反演方法及在南岭—武夷—东南沿海成矿带岩石圈三维精细结构成像的应用”(41974101)

中央高校基本科研业务费专项资金

地质过程与矿产资源国家重点实验室科技部专项经费

Received: 2021-07-22 Revised: 2021-10-18

作者简介 About authors

杨婧(1999-),女,硕士研究生,主要从事重震数据处理和反演研究工作。Email: 2010200015@cugb.edu.cn

摘要

重力异常曲化平即将起伏观测面上的重力异常延拓到平面上,从而为频率域处理和反演提供平坦观测面的重力异常数据。本文在常规插值—迭代法基础上,给出重力异常曲化平的改进插值—迭代法,即在异常迭代修正过程中引入起伏观测面修正因子,加快曲化平迭代收敛,促进曲化平效果提升。理论模型试验表明本文方法适用于观测面起伏较大、延拓跨度大的复杂条件下的稳定、有效曲化平,效果优于常规插值—迭代法。川滇地区实际重力异常数据试验表明本文曲化平方法有效增强了异常信号和细节特征,为后续处理和解释提供可靠数据。

关键词： 曲化平; 重力异常; 插值—迭代法; 川滇地区

Abstract

Gravity anomaly continuation from undulating surface to plane can provide gravity data on a flat horizontal plane for frequency-domain data processing and inversion. Based on the theory of conventional interpolation-iteration methods, this study proposed a modified interpolation-iteration method by introducing a correction factor of the undulating observation surface in the iteration and correction process. The improved method accelerated the iterative convergence speed and promoted the continuation effects. The theoretical model-based tests show that this method can be used to achieve stable and effective large-span gravity anomaly continuation from greatly undulating surface to plane. The continuation results presented by this method are better than those obtained using conventional interpolation-iteration methods. The application of Bouguer gravity anomaly data of the Sichuan-Yunnan region demonstrates that the modified interpolation-iteration method effectively enhanced anomalous signals and details and can provide reliable data for subsequent processing and interpretation.

Keywords： continuation from undulating surface to plane; gravity anomaly; interpolation-iteration method; Sichuan-Yunnan region

PDF (1920KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

杨婧, 郭良辉. 重力异常曲化平的改进插值—迭代法[J]. 物探与化探, 2022, 46(1): 123-129 doi:10.11720/wtyht.2022.1403

YANG Jing, GUO Liang-Hui. A modified interpolation-iteration method for gravity anomaly continuation from undulating surface to plane[J]. Geophysical and Geochemical Exploration, 2022, 46(1): 123-129 doi:10.11720/wtyht.2022.1403

0 引言

高精度重力测量可以为区域地质构造研究、资源与能源勘查等提供重要的基础数据。野外起伏观测面采集的重力数据,经过各项校正得到的重力异常仍然是原起伏观测面上的异常,反映地下密度不均匀体在起伏观测面上的重力效应^[1]。频率域算法是地球物理领域的一种常用快速算法,已广泛应用于重力数据处理与反演解释,比如频率域延拓^[1,2]、频率域优化滤波^[3]、频率域界面反演^[4]、频率域三维成像^[5]等。然而,这些频率域算法通常要求重力异常所对应的观测面为平面。在观测面起伏较大情况下,若对起伏面重力数据直接应用常规频率域算法做处理与反演,则将带来不可忽略的处理偏差。因此,在起伏观测面条件下,重力异常数据从起伏观测面化到平面上,即曲化平,是开展频率域精细处理解释的必要步骤。

当前重力数据的曲化平方法较多,Pilkington等^[6]将其分为基于源和基于场的两大类。基于源的方法主要有等效源法^{[7,8,9,10,11,12]},它先对起伏观测面异常反演获得地下等效场源,然后再正演计算等效源在任意水平面引起的异常,该类方法通常计算量较大,适用于中、小规模数据的曲化平处理。基于场的方法有有限调和级数法^[13]、泰勒级数法^[14,15]、有限元法^[16,17]、积分法^[18,19]、逐次逼近法^[20]等,该类方法无需反演场源,通常需要向上延拓或向下延拓到某中间层面进而获取最终平面的异常。向上延拓是稳定的,但向下延拓通常是不稳定的,这是由于它具有放大作用,对高频噪声干扰比较敏感,延拓跨度越大,高频干扰越严重。因此,基于向下延拓的曲化平方法一般迭代收敛速度慢、精度有限,仅适用于延拓跨度小的曲化平。

徐世浙等^[21]提出了位场(重、磁场)曲化平的插值—迭代法,属于基于场的方法类,它计算简单快速,可以实现延拓跨度达10~20倍点距的大规模数据的稳定曲化平。该方法基本思路是,将起伏观测面(A)的顶部水平面(B_m)与底部水平面(B)之间的空间剖分为若干等间距的中间层水平面(B_i),见图1所示;接着,把起伏观测面(A)的位场值放到底部水平面(B)并作为该水平面的位场初始值;然后,通过频率域向上延拓计算出各中间水平面(B_i)的位场值,进而从各中间水平面(B_i)集合的位场值插值出起伏观测面(A)的位场近似值;之后,计算起伏观测面(A)的位场实测值与近似值的残差,并用于对底部水平面(B)的位场值进行修正;如此反复迭代,直至起伏观测面(A)的位场实测值与近似值的残差达到误差容限。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 插值—迭代法曲化平示意

Fig.1 The diagram of interpolation-iteration method for gravity anomaly continuation from undulating surface to plane

然而,在实际应用中,本文发现上述插值—迭代方法^[21]在观测面起伏较大、延拓跨度较大的复杂条件下的曲化平效果还是有限的(见后面理论模型数据试验的图4所示)。为此,本文对上述插值—迭代方法做进一步的改进,在异常迭代修正过程中引入起伏观测面修正因子,加快曲化平迭代收敛,促进曲化平效果提升。本文首先介绍改进的插值—迭代法的方法原理和算法流程,然后利用理论模型和川滇地区实际数据试验验证方法的有效性。

1 方法原理

本文曲化平方法是在徐世浙等^[21]的插值—迭代方法基础上改进得到, 如图1所示,假定起伏观测面A和底部水平面B,曲化平的目标是由起伏观测面A的位场数据p_A通过曲化平处理获得底部水平面B的位场数据p_B。根据起伏观测面的高程最大值给定顶部水平面B_m,进而将底部水平面和顶部水平面之间的空间按照等间距剖分为m层,各中间层水平面为B_i(i=1,2,…,m)。

本文算法流程如下图2所示。首先,将起伏观测面A的实测位场数据p_A直接赋值到底部水平面B上,作为底部水平面B的位场初值 $p_{B}^{(k)}$ (k=0),有

(1)

p_{B}^{(1)} = p_{A},

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 本文改进的插值—迭代法算法流程

Fig.2 The workflow of the modified interpolation-iteration method

接着,利用频率域向上延拓方法^[2]将底部水平面B的位场初值 $p_{B}^{(k)}$ (k=0)分别延拓到各中间层水平面B_i上,得到m个中间层水平面的位场值 $p_{B}^{(k)}$ (k=0),再利用线性插值方法将起伏观测面A上各测点的场值计算出来,得到起伏观测面A的位场近似值 $p_{A}^{(k)}$ 。

然后,计算起伏观测面A的位场实测值p_A与近似值 $p_{A}^{(k)}$ 的残差,引入修正因子S,用于对底部水平面B的位场值进行修正:

(2)

p_{B}^{(k + 1)} = p_{B}^{(k)} + S (p_{A} - p_{A}^{(k)}) 。

在徐世浙等^[21]的插值—迭代方法中,修正因子S取常数,且通常取为1,即观测面上各测点的修正因子相同,导致在相同迭代次数下,观测面起伏变化大的地方曲化平迭代收敛速度慢、效果差些。因此,本文对修正因子S做了改进,采用与起伏观测面相关的修正因子,从而加快曲化平迭代收敛速度,提升曲化平效果。本文的修正因子公式如下:

(3)

S = {(\frac{T - T_{\min}}{T_{\max} - T_{\min}})}^{n},

其中,T为起伏观测面高程,T_min和T_max分别是起伏观测面高程的最小值和最大值。可知,观测面起伏越大,修正因子S将越大,曲化平收敛速度将越快,反之,观测面起伏越小,修正因子S将越小,曲化平收敛速度将越慢。一般,n取非负数,若n取值过大,则修正因子S将越大;反之,n取值越小,修正因子S将越小;当n=0时,修正因子S=1,等效于徐世浙等^[21]的算法。

之后,对式(2)得到的底部水平面B的位场修正值,重复上述两个步骤,依次迭代修正,直到起伏观测面A的位场实测值p_A与近似值 $p_{A}^{(k)}$ 的残差达到指定的误差容限,或者迭代次数达到指定的次数容限,则停止迭代,最后一次迭代得到的底部水平面B的位场修正值 $p_{B}^{(k + 1)}$ 就是最终的曲化平结果。

2 理论模型数据试验

理论模型为一个直立长方体,长、宽、高分别为 2 000、4 000、1 000 m,中心埋深1 500 m,剩余密度为1.0 g/cm³。假定测网E向和N向范围均为0~10 km、点距均为0.1 km,则理论模型在海拔0 km的水平面上产生的理论重力异常如图3a所示,异常等值线为近椭圆状,最小值为0.23 mGal,最大值为12.65 mGal,平均值为2.41 mGal。假定起伏观测面高程如图3b所示,高程最小值为1.46 m,最大值为2 023.55 m,观测面高差是测网点距的20多倍。理论模型在起伏观测面上产生的理论重力异常如图3c所示,异常等值线为不规则形状,且与起伏观测面形状有一定相关,异常最小值为0.23 mGal,最大值为7.52 mGal,平均值为1.88 mGal。因此,受观测面起伏影响,理论重力异常变为复杂,难以直接处理解释,值得优先进行曲化平处理。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 理论模型起伏观测面与重力异常

a—海拔0 m平面的理论重力异常;b—起伏观测面高程;c—起伏观测面的理论重力异常

Fig.3 Undulating observation surface of theoretical model and its gravity anomaly

a—theoretical gravity anomaly at 0 m altitude;b—elevation of undulating observation surface;c—theoretical gravity anomaly at undulating observation surface

图4a和图4b显示了常规插值—迭代方法^[21]迭代20次的曲化平结果及其与海拔0 m平面理论重力异常(图3a)的残差,修正因子中n取0,即S=1。常规插值—迭代法曲化平增强了有效信号,最小值为0.23 mGal,最大值为13.78 mGal,但结果与海拔0 m平面理论重力异常有明显的差别。异常残差特征与起伏观测面比较相似,在观测面平缓区残差较小,起伏区残差较大,残差最小值为-2.15 mGal,最大值为3.73 mGal,标准差为0.39 mGal。因此,在观测面起伏较大、延拓跨度较大的情况下,常规插值—迭代方法的曲化平效果是有限的。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 不同曲化平方法的结果及与海拔0 m平面理论重力异常对比

a、b—常规插值-迭代法曲化平结果及与海拔0 m平面理论重力异常残差;c、d—本文改进的插值-迭代法曲化平结果及与海拔0 m平面理论重力异常残差

Fig.4 Comparison between the results of different interpolation-iteration methods and the gravity forward values of plane surface

a、b—the result of the routine interpolation-iteration method for continuation from undulating surface to plane and its deviation from the values of plane surface;c、d—the result of the modified interpolation-iteration method and its deviation from the values of plane surface

图4c和图4d显示了本文改进的插值—迭代方法迭代20次的曲化平结果及其与海拔0 m平面理论重力异常(图3a)的残差,其中修正因子采用式(3)且经过测试,选择效果较好的n=1.5。

图5显示了剖面A不同曲化平方法结果与起伏观测面、海拔0 m平面的理论重力异常、起伏观测面高程的对比。观察发现,水平面的理论重力异常出现一个异常峰值,对应地下密度异常体,起伏观测面下的重力异常出现两个异常峰值,其位置与起伏观测面相对应,这说明起伏观测面扭曲了地下密度体的异常信息,观测面起伏越高,扭曲作用越大。n的取值与实际研究区观测面起伏程度相关,n=0时,S为常数1,等效于徐世浙等^[21]的算法,其均方根误差为0.1551 mGal;n=1时,其均方根误差0.044 5 mGal;n=1.5时,其均方根误差最小,为0.028 7 mGal;n=2时,其均方根误差0.079 4 mGal,故本文选n=1.5。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 剖面A不同修正因子S曲化平的对比

Fig.5 Different S comparison along profile A

改进的插值—迭代法有效增强了信号,曲化平后其最小值为0.23 mGal,最大值为12.33 mGal,结果与海拔0 m平面理论重力异常差别较小。异常残差最小值为-0.86 mGal,最大值为1.03 mGal,标准差为0.17 mGal,我们分析这一微小差异是由高程差距引起收敛速度不一致产生的。因此,与常规插值—迭代法相比,本文改进的插值—迭代法曲化平结果更接近海拔0 m平面理论重力异常值,效果更优,适用于大规模数据、复杂起伏观测面、大延拓跨度的曲化平。

3 实际数据试验

川滇地区构造上位于印支块体、青藏高原块体和扬子块体的交汇处,区域构造动力作用复杂,在缅甸弧东向俯冲、印度板块藏东构造结NE向楔体挤出和俯冲、青藏高原隆升和SE向挤出等共同作用下,该地区整体向E—SE方向挤出,形成了大陆内部典型的剪切、拉张和推覆构造,表现出地壳变形速率高、断层运动剧烈、强震原地复发周期短等区域特征。因此,川滇地区是研究青藏高原隆升过程与大陆强震孕育机理的天然实验场,也是我国区域防震减灾的重要地区。高精度重力数据是川滇地区深部结构与构造研究的基础数据。由于该地区地形上呈现青藏高原、云贵高原、四川盆地等复杂、多样化的起伏形态,因此对该地区起伏地形面的重力异常数据作常规频率域处理和解释应用之前,值得优先作曲化平处理。

本文从ICGEM官网(http://icgem.gfz-potsdam.de/calcgrid)下载了川滇地区地形高程网格数据和布格重力异常网格数据,东经99°~108°,北纬24°~32°,网格大小均为0.1°×0.1°。其中,布格重力异常数据是由地球重力场模型 (earth gravitational model 2008)的自由空气重力异常数据经过标准的地形校正和中间层校正得到的。本文利用窗口大小为3的均值滤波对研究区地形高程数据做了平滑处理,结果见图6a所示,研究区地形起伏较大,高程最大值为4 882 m,最小值为269.2 m,高差约为 4 612.8 m。由于原始布格重力异常数据存在明显的高频噪声干扰,这主要是由于自由空气重力异常数据和地形数据的分辨率与精度不一致造成的^[22],因此,本文对原始布格重力异常数据作了截止波长为100 km的低通滤波,得到去噪后的川滇地区布格重力异常,如图6b所示,异常最大值为-64.51 mGal,最小值为-495.5 mGal,青藏高原地区异常较低,四川盆地等海拔低地区异常较高,云贵高原介于两者之间。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 川滇地区曲化平前后重力异常对比

a—川滇地区地形高程;b—去噪后的布格重力异常;c—本文改进插值—迭代法曲化平后布格重力异常;d—曲化平前后的异常残差

Fig.6 The Bouguer gravity comparison before and after continuation from undulating surface to plane of Sichuan-Yunnan region

a—the elevation of Sichuan-Yunnan region;b—the denoised Bouguer gravity anomaly;c—the Bouguer gravity anomaly after continuation from undulating surface to plane by using the modified interpolation-iteration method;d—the anomaly deviations before and after continuation from undulating surface to plane

本文应用改进的插值—迭代法对去噪后的布格重力异常作曲化平处理,修正因子采用式(3),n=1.5,经过50次迭代,将其化到海拔0 m的位置,结果如图6c所示。曲化平后的布格重力异常最大值为-64.47 mGal,最小值为-508.9 mGal,山区异常信号有效增强,尤其是在川西地区增强较为明显、异常细节更为突出。图6d显示了曲化平前后的布格重力异常残差,最大值为13.36 mGal,最小值为-8.45 mGal,残差变化主要集中在青藏高原和云贵高原等山区,与地形有较强对应关系。

为进一步探讨地形对重力异常值的影响及曲化平的效果,本文选取一条斜向穿过研究区的剖面B(位置如图6a所示),该剖面地势复杂,地形起伏变化较大。该剖面曲化平前后的重力异常对比(图7)显示,地形起伏较大地区的曲化平效果显著,地形起伏较为平缓地区的曲化平效果较弱,这一结果与前文的分析相符。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 剖面B的对比

Fig.7 The gravity anomaly comparison and the elevation along profile B

本文以常规频率域垂直导数换算为例作说明曲化平的意义。导数换算是重力异常分析解释的一个常用步骤,用于突出浅部场源、分析构造边界等。图8a和图8b显示了曲化平前后的布格重力异常的垂直导数,曲化平前的垂直导数数值范围为-0.002 9~0.002 3 mGal/m,曲化平后的数值范围转为-0.004 4~0.003 7 mGal/m。可见,曲化平后的异常幅值得到明显增强,而且刻画出更多、更清晰的异常变化细节,从而反映出更多的浅部场源细节和构造边界信息。

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8 川滇地区曲化平前后垂向导数对比

a—布格重力异常的垂向导数;b—曲化平后布格重力异常的垂向导数

Fig.8 The vertical derivative of Bouguer gravity comparison before and after continuation from undulating surface to plane of Sichuan-Yunnan region

a—the vertical derivative of Bouguer gravity anomaly;b—the vertical derivative of Bouguer gravity

图8中标注了2008年汶川8级地震、2013年芦山7级地震、2014年鲁甸6.5级地震、2021年漾濞6.4级地震的震中位置,易见经过曲化平这些地震震中周缘地区的重力异常变化特征更为鲜明,这有助于后续的深部结构与构造解释研究。

4 结论与讨论

本文在重力异常曲化平的常规插值—迭代法基础上给出了一种改进的插值—迭代法,即在异常迭代修正过程中引入起伏观测面修正因子,加快曲化平迭代收敛,促进曲化平效果提升,实现适用于观测面起伏较大、延拓跨度较大的复杂条件下重力异常曲化平。理论模型和川滇地区实际数据试验验证了本文方法的有效性,效果优于常规插值—迭代法。本文改进的插值—迭代方法不仅适用于重力异常数据,也适用于磁异常数据以及梯度、三分量等数据。本文曲化平的结果可进一步通过常规平化平、平化曲处理实现任意平面或曲面的延拓。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

曾华霖

. 重力场与重力勘探[M]. 北京: 地质出版社, 2005.

[本文引用: 2]

Zeng H

. Gravity field and gravity exploration[M]. Beijing: Geological Publishing House, 2005.

[本文引用: 2]

[2]

Blakely R

. Potential theory in gravity and magnetic applications [M]. Cambridge: Cambridge University Press, 1995:316-320.

[本文引用: 2]

[3]

Guo

, Meng

, Chen

, et al.

Preferential filtering for gravity anomaly separation

[J]. Computers & Geosciences, 2013,51:247-254.

DOI:10.1016/j.cageo.2012.09.012 URL [本文引用: 1]

[4]

Oldenburg D

The inversion and interpretation of gravity anomalies

[J]. Geophysics, 1974,39(4):526-536.

DOI:10.1190/1.1440444 URL [本文引用: 1]

[5]

Cui

, Guo

A wavenumber-domain iterative approach for rapid 3-D imaging of gravity anomalies and gradients

[J]. IEEE Access, 2019,7(1):34179-34188.

DOI:10.1109/Access.6287639 URL [本文引用: 1]

[6]

Pilkington

, Boulanger

Potential field continuation between arbitrary surfaces-Comparing methods

[J]. Geophysics, 2017,82(3):J9-J25.

DOI:10.1190/geo2016-0210.1 URL [本文引用: 1]

[7]

Dampney C N

The equivalent source technique

[J]. Geophysics, 1969,34(1):39-53.

DOI:10.1190/1.1439996 URL [本文引用: 1]

[8]

David A

Equivalent sources used as an analytic base for processing total magnetic field profiles

[J]. Geophysics, 1973,38(2):339-348.

DOI:10.1190/1.1440344 URL [本文引用: 1]

[9]

Bhattacharyya B

, Chan K

Reduction of magnetic and gravity data on an arbitrary surface acquired in a region of high topographic relief

[J]. Geophysics, 1977,42(7):1411-1430.

DOI:10.1190/1.1440802 URL [本文引用: 1]

[10]

Hansen R

, Miyazaki

Continuation of potential fields between arbitrary surfaces

[J]. Geophysics, 1984,49(6):787-795.

DOI:10.1190/1.1441707 URL [本文引用: 1]

[11]

Pilkington

, Urquhart

W E S

Reduction of potential field data to a horizontal plane

[J]. Geophysics, 1990,55(5):549-555.

DOI:10.1190/1.1442866 URL [本文引用: 1]

[12]

Xia

, Sprowl D

, Adkins-Heljeson

Correction of topographic distortions in potential-field data: A fast and accurate approach

[J]. Geophysics, 1993,58(4):515-523.

DOI:10.1190/1.1443434 URL [本文引用: 1]

[13]

Henderson R

, Cordell

Reduction of unevenly spaced potential field data to a horizontal plane by means of finite harmonic series

[J]. Geophysics, 1971,36(5):1046-1046.

[本文引用: 1]

[14]

Guspi

Frequency-domain reduction of potential field measurements to a horizontal plane

[J]. Geoexploration, 1987,24(2):87-98.

DOI:10.1016/0016-7142(87)90083-4 URL [本文引用: 1]

[15]

Pilkington

, Thurston J

Draping corrections for aeromagnetic data: Line versus grid-based approaches

[J]. Exploration Geophysics, 2001,32(2):95-101.

DOI:10.1071/EG01095 URL [本文引用: 1]

[16]

程振炎

重磁场的有限元法曲化平

[J]. 物探与化探, 1981,5(3):153-158.

[本文引用: 1]

Cheng Z

The finite element method of gravity and magnetic field

[J]. Geophysical and Geochemical Exploration, 1981,5(3):153-158.

[本文引用: 1]

[17]

Xu S

The boundary element method in geophysics

[M]. Tulsa: Society of Exploration Geophysicists, 2001:63-67.

[本文引用: 1]

[18]

Xu S

, Yang C

, Dai S

, et al.

A new method for continuation of 3D potential fields to a horizontal plane

[J]. Geophysics, 2003,68(6):1917-1921.

DOI:10.1190/1.1635045 URL [本文引用: 1]

[19]

徐世浙

位场延拓的积分—迭代法

[J]. 地球物理学报, 2006,49(4):1176-1182.

[本文引用: 1]

Xu S

The interpolation-iteration method for continuation of potential fields

[J]. Chinese Journal of Geophysics, 2006,49(4):1176-1182.

[本文引用: 1]

[20]

刘东甲, 洪天求, 廖旭涛, 等.

位场曲化平积分方程的迭代解

[J]. 地球物理学报, 2012,55(10):3467-3476.

[本文引用: 1]

Liu D

, Hong T

, Liao X

, et al.

Iterative solution of integral equation for potential field continuation from irregular surface to a horizontal plane

[J]. Chinese Journal of Geophysics, 2012,55(10):3467-3476.

[本文引用: 1]

[21]

徐世浙, 余海龙.

位场曲化平的插值—迭代法

[J]. 地球物理学报, 50(6):193-197.

[本文引用: 7]

Xu S

, Yu H

The interpolation-iteration method for potential fields continuation from undulating surface to plane

[J]. Chinese Journal of Geophysics, 2007,50(6):193-197.

[本文引用: 7]

[22]

郭良辉, 孟小红, 石磊, 等.

优化滤波方法及其在中国大陆布格重力异常数据处理中的应用

[J]. 地球物理学报, 2012,55(12):4078-4088.

[本文引用: 1]

Guo L

, Meng X

, Shi

, et al.

Preferential filtering method and its application to Bouguer gravity anomaly of Chinese continent

[J]. Chinese Journal of Geophysics, 2012,55(12):4078-4088.

[本文引用: 1]

2005

... 高精度重力测量可以为区域地质构造研究、资源与能源勘查等提供重要的基础数据.野外起伏观测面采集的重力数据,经过各项校正得到的重力异常仍然是原起伏观测面上的异常,反映地下密度不均匀体在起伏观测面上的重力效应^[1].频率域算法是地球物理领域的一种常用快速算法,已广泛应用于重力数据处理与反演解释,比如频率域延拓^[1,2]、频率域优化滤波^[3]、频率域界面反演^[4]、频率域三维成像^[5]等.然而,这些频率域算法通常要求重力异常所对应的观测面为平面.在观测面起伏较大情况下,若对起伏面重力数据直接应用常规频率域算法做处理与反演,则将带来不可忽略的处理偏差.因此,在起伏观测面条件下,重力异常数据从起伏观测面化到平面上,即曲化平,是开展频率域精细处理解释的必要步骤. ...

... [1,2]、频率域优化滤波^[3]、频率域界面反演^[4]、频率域三维成像^[5]等.然而,这些频率域算法通常要求重力异常所对应的观测面为平面.在观测面起伏较大情况下,若对起伏面重力数据直接应用常规频率域算法做处理与反演,则将带来不可忽略的处理偏差.因此,在起伏观测面条件下,重力异常数据从起伏观测面化到平面上,即曲化平,是开展频率域精细处理解释的必要步骤. ...

2005

1995

... 接着,利用频率域向上延拓方法^[2]将底部水平面B的位场初值

p_{B}^{(k)}

(k=0)分别延拓到各中间层水平面B_i上,得到m个中间层水平面的位场值

p_{B}^{(k)}

(k=0),再利用线性插值方法将起伏观测面A上各测点的场值计算出来,得到起伏观测面A的位场近似值

p_{A}^{(k)}

. ...

Preferential filtering for gravity anomaly separation

2013

The inversion and interpretation of gravity anomalies

1974

A wavenumber-domain iterative approach for rapid 3-D imaging of gravity anomalies and gradients

2019

Potential field continuation between arbitrary surfaces-Comparing methods

2017

... 当前重力数据的曲化平方法较多,Pilkington等^[6]将其分为基于源和基于场的两大类.基于源的方法主要有等效源法^{[7,8,9,10,11,12]},它先对起伏观测面异常反演获得地下等效场源,然后再正演计算等效源在任意水平面引起的异常,该类方法通常计算量较大,适用于中、小规模数据的曲化平处理.基于场的方法有有限调和级数法^[13]、泰勒级数法^[14,15]、有限元法^[16,17]、积分法^[18,19]、逐次逼近法^[20]等,该类方法无需反演场源,通常需要向上延拓或向下延拓到某中间层面进而获取最终平面的异常.向上延拓是稳定的,但向下延拓通常是不稳定的,这是由于它具有放大作用,对高频噪声干扰比较敏感,延拓跨度越大,高频干扰越严重.因此,基于向下延拓的曲化平方法一般迭代收敛速度慢、精度有限,仅适用于延拓跨度小的曲化平. ...

The equivalent source technique

1969

Equivalent sources used as an analytic base for processing total magnetic field profiles

1973

Reduction of magnetic and gravity data on an arbitrary surface acquired in a region of high topographic relief

1977

Continuation of potential fields between arbitrary surfaces

1984

Reduction of potential field data to a horizontal plane

1990

Correction of topographic distortions in potential-field data: A fast and accurate approach

1993

Reduction of unevenly spaced potential field data to a horizontal plane by means of finite harmonic series

1971

Frequency-domain reduction of potential field measurements to a horizontal plane

1987

Draping corrections for aeromagnetic data: Line versus grid-based approaches

2001

重磁场的有限元法曲化平

1981

重磁场的有限元法曲化平

1981

The boundary element method in geophysics

2001

A new method for continuation of 3D potential fields to a horizontal plane

2003

位场延拓的积分—迭代法

2006

位场延拓的积分—迭代法

2006

位场曲化平积分方程的迭代解

2012

位场曲化平积分方程的迭代解

2012

位场曲化平的插值—迭代法

2007

... 徐世浙等^[21]提出了位场(重、磁场)曲化平的插值—迭代法,属于基于场的方法类,它计算简单快速,可以实现延拓跨度达10~20倍点距的大规模数据的稳定曲化平.该方法基本思路是,将起伏观测面(A)的顶部水平面(B_m)与底部水平面(B)之间的空间剖分为若干等间距的中间层水平面(B_i),见图1所示;接着,把起伏观测面(A)的位场值放到底部水平面(B)并作为该水平面的位场初始值;然后,通过频率域向上延拓计算出各中间水平面(B_i)的位场值,进而从各中间水平面(B_i)集合的位场值插值出起伏观测面(A)的位场近似值;之后,计算起伏观测面(A)的位场实测值与近似值的残差,并用于对底部水平面(B)的位场值进行修正;如此反复迭代,直至起伏观测面(A)的位场实测值与近似值的残差达到误差容限. ...

... 然而,在实际应用中,本文发现上述插值—迭代方法^[21]在观测面起伏较大、延拓跨度较大的复杂条件下的曲化平效果还是有限的(见后面理论模型数据试验的图4所示).为此,本文对上述插值—迭代方法做进一步的改进,在异常迭代修正过程中引入起伏观测面修正因子,加快曲化平迭代收敛,促进曲化平效果提升.本文首先介绍改进的插值—迭代法的方法原理和算法流程,然后利用理论模型和川滇地区实际数据试验验证方法的有效性. ...

... 本文曲化平方法是在徐世浙等^[21]的插值—迭代方法基础上改进得到, 如图1所示,假定起伏观测面A和底部水平面B,曲化平的目标是由起伏观测面A的位场数据p_A通过曲化平处理获得底部水平面B的位场数据p_B.根据起伏观测面的高程最大值给定顶部水平面B_m,进而将底部水平面和顶部水平面之间的空间按照等间距剖分为m层,各中间层水平面为B_i(i=1,2,…,m). ...

... 在徐世浙等^[21]的插值—迭代方法中,修正因子S取常数,且通常取为1,即观测面上各测点的修正因子相同,导致在相同迭代次数下,观测面起伏变化大的地方曲化平迭代收敛速度慢、效果差些.因此,本文对修正因子S做了改进,采用与起伏观测面相关的修正因子,从而加快曲化平迭代收敛速度,提升曲化平效果.本文的修正因子公式如下: ...

... 其中,T为起伏观测面高程,T_min和T_max分别是起伏观测面高程的最小值和最大值.可知,观测面起伏越大,修正因子S将越大,曲化平收敛速度将越快,反之,观测面起伏越小,修正因子S将越小,曲化平收敛速度将越慢.一般,n取非负数,若n取值过大,则修正因子S将越大;反之,n取值越小,修正因子S将越小;当n=0时,修正因子S=1,等效于徐世浙等^[21]的算法. ...

... 图4a和图4b显示了常规插值—迭代方法^[21]迭代20次的曲化平结果及其与海拔0 m平面理论重力异常(图3a)的残差,修正因子中n取0,即S=1.常规插值—迭代法曲化平增强了有效信号,最小值为0.23 mGal,最大值为13.78 mGal,但结果与海拔0 m平面理论重力异常有明显的差别.异常残差特征与起伏观测面比较相似,在观测面平缓区残差较小,起伏区残差较大,残差最小值为-2.15 mGal,最大值为3.73 mGal,标准差为0.39 mGal.因此,在观测面起伏较大、延拓跨度较大的情况下,常规插值—迭代方法的曲化平效果是有限的. ...

... 图5显示了剖面A不同曲化平方法结果与起伏观测面、海拔0 m平面的理论重力异常、起伏观测面高程的对比.观察发现,水平面的理论重力异常出现一个异常峰值,对应地下密度异常体,起伏观测面下的重力异常出现两个异常峰值,其位置与起伏观测面相对应,这说明起伏观测面扭曲了地下密度体的异常信息,观测面起伏越高,扭曲作用越大.n的取值与实际研究区观测面起伏程度相关,n=0时,S为常数1,等效于徐世浙等^[21]的算法,其均方根误差为0.1551 mGal;n=1时,其均方根误差0.044 5 mGal;n=1.5时,其均方根误差最小,为0.028 7 mGal;n=2时,其均方根误差0.079 4 mGal,故本文选n=1.5. ...

位场曲化平的插值—迭代法

2007

优化滤波方法及其在中国大陆布格重力异常数据处理中的应用

2012

... 本文从ICGEM官网(http://icgem.gfz-potsdam.de/calcgrid)下载了川滇地区地形高程网格数据和布格重力异常网格数据,东经99°~108°,北纬24°~32°,网格大小均为0.1°×0.1°.其中,布格重力异常数据是由地球重力场模型 (earth gravitational model 2008)的自由空气重力异常数据经过标准的地形校正和中间层校正得到的.本文利用窗口大小为3的均值滤波对研究区地形高程数据做了平滑处理,结果见图6a所示,研究区地形起伏较大,高程最大值为4 882 m,最小值为269.2 m,高差约为 4 612.8 m.由于原始布格重力异常数据存在明显的高频噪声干扰,这主要是由于自由空气重力异常数据和地形数据的分辨率与精度不一致造成的^[22],因此,本文对原始布格重力异常数据作了截止波长为100 km的低通滤波,得到去噪后的川滇地区布格重力异常,如图6b所示,异常最大值为-64.51 mGal,最小值为-495.5 mGal,青藏高原地区异常较低,四川盆地等海拔低地区异常较高,云贵高原介于两者之间. ...

优化滤波方法及其在中国大陆布格重力异常数据处理中的应用

2012

〈

〉