重力场向下延拓Milne法

重力场向下延拓能够突出局部和浅部的异常信息,分离叠加的异常特征.但是向下延拓通常具有过程不稳定、下延深度小、结果不准确等问题.针对向下延拓所存在的不足,本文利用重力场及其垂向一阶导数,基于辛普森(Simpson)求积公式,推导出重力场向下延拓米尔尼(Milne)公式.将本文向下延拓方法应用于模型数据,向下延拓模型结果及误差曲线表明,相对于向下延拓快速傅里叶变换(FFT)法和积分迭代法,向下延拓Milne法的深度更大,相对误差更小;相对模型值,向下延拓Milne法能够获得稳定且准确的结果.对加拿大乃查科(Nechako)盆地地区实测航空重力数据进行本文方法向下延拓验证,处理结果表明,相对于实测异常,本文方法向下延拓结果能够很好还原实测数据,并且在进一步向下延拓中反映原始异常的趋势,增强局部和细小异常信息.     

Depth Estimates for Slingram Electromagnetic Anomalies from Dipping Sheet-like Bodies by the Normalized Full Gradient Method

The Normalized Full Gradient (NFG) method was proposed in the mid–1960s and was generally used for the downward continuation of the potential field data. The method eliminates the side oscillations which appeared on the continuation curves when passing through anomalous body depth. In this study, the NFG method was applied to Slingram electromagnetic anomalies to obtain the depth of the anomalous body. Some experiments were performed on the theoretical Slingram model anomalies in a free space environment using a perfectly conductive thin tabular conductor with an infinite depth extent. The theoretical Slingram responses were obtained for different depths, dip angles and coil separations, and it was observed from NFG fields of the theoretical anomalies that the NFG sections yield the depth information of top of the conductor at low harmonic numbers. The NFG sections consisted of two main local maxima located at both sides of the central negative Slingram anomalies. It is concluded that these two maxima also locate the maximum anomaly gradient points, which indicates the depth of the anomaly target directly. For both theoretical and field data, the depth of the maximum value on the NFG sections corresponds to the depth of the upper edge of the anomalous conductor. The NFG method was applied to the in-phase component and correct depth estimates were obtained even for the horizontal tabular conductor. Depth values could be estimated with a relatively small error percentage when the conductive model was near-vertical and/or the conductor depth was larger.     

Application of Normalized Full Gradient Method to Self Potential (SP) Data

Recently, Normalized Full Gradient (NFG) method has widespread applications to natural potential fields, especially in gravity and magnetic. In this study, usage of NFG in Self-Potential (SP) data evaluation is tested. Results are compared to other SP interpretation methods. The NFG method is applied to synthetic and field SP data. As a consequence of application of the method to the anomalies of spherical, cylindrical and vertical sheet models, whose theoretical structures are explicit, the structures were found very close to their actual locations. In order to see the capability of the method in detecting the number of sources, NFG method was applied to different spherical models at different depths and locations. The least-squares inverse solution was applied to the same models and NFG method was found more powerful in detecting model structure. Sensitivity of NFG method for application to noisy data is also tested. An anomaly is generated by adding a random noise to two close sphere SP anomalies. The method seems to work for the two close spheres at high S/N ratio. Then, NFG method was applied to two field examples. The first one is the cross section taken from the SP anomaly map of the Ergani-Süleymanköy (Turkey) copper mine. The depth of the mineral deposit at that site was found about 38 m from the ground level. This result is well matched to previous studies. NFG was also applied to SP data from Seferihisar Izmir (Western Turkey) geothermal field and the location of the point source was determined. The field data of this site have already been modeled by the thermoelectric source (coupling) solution method. When these two methods are compared, they seem to support each other. It is concluded that the NFG method works perfectly when the structure model is simple. It is observed that natural potential sources close to earth’s surface are identified by the method more accurately at greater harmonics, while deep sources are identified at lesser harmonics. It produces reasonable results for noisy multi-source models than the other parameter identification methods (inverse solution, power spectrum, etc.).     

位场向下延拓的波数域广义逆算法

位场向下延拓是位场数据处理和反演中的重要运算，但是它的不稳定性影响了它在许多处理和反演方法技术中的应用.本文通过把位场向下延拓视为向上延拓的反问题，得到向下延拓的褶积型线性积分方程，再利用Fourier变换矩阵的正交对称特性，并结合矩阵的奇异值分解和广义逆原理，提出了一种稳定的不需要进行求逆运算的位场向下延拓广义逆方法——波数域广义逆算法，解决了位场大深度向下延拓的不稳定性问题.把这种方法用于三维理论模型数据和实际磁场数据的向下延拓获得了理想的结果.     

补偿向下延拓方法研究及应用

位场向下延拓是重、磁处理和解释的常用方法,但其不稳定性限制了其在资料处理及反演中的应用.本文基于补偿圆滑滤波思想以及空间域向下延拓迭代法,通过逐次补偿的办法实现位场的稳定向下延拓.同时,在频率域空间给出了该下延方法的频率域响应因子,并讨论了其低通滤波特性,理论模型和实际位场资料试验表明该方法向下延拓稳定性具有较高的延拓精度.将其应用于重力密度界面反演中,改进反演的稳定性,实际莫霍界面反演表明下延因子具备实用性.     

位场向下延拓组合滤波器的设计和应用

一 重、磁异常向下延拓能将水平方向迭加的异常分开,并能突出异常细节,从而丰富异常的解释内容,提高异常解释的可靠性.本文从对引起向下延拓不稳定的各种因素的讨论入手,结合对一些向下延拓方法在频率域的分析,提出了一种与已知资料频谱特征相“匹配”的组合滤波器。经理论模型和实例考核,说明该滤波器既能保证向下延拓有较高的精度,又有较强的抗干扰能力。     

Atmospheric effects on satellite gravity gradiometry data

Atmospheric masses play an important role in precise downward continuation and validation of satellite gravity gradiometry data. In this paper we present two alternative ways to formulate the atmospheric potential. Two density models for the atmosphere are proposed and used to formulate the external and internal atmospheric potentials in spherical harmonics. Based on the derived harmonic coefficients, the direct atmospheric effects on the satellite gravity gradiometry data are investigated and presented in the orbital frame over Fennoscandia. The formulas of the indirect atmospheric effects on gravity anomaly and geoid (downward continued quantities) are also derived using the proposed density models. The numerical results show that the atmospheric effect can only be significant for precise validation or inversion of the GOCE gradiometric data at the mE level.     

位场向下延拓的改进迭代维纳滤波法

根据维纳滤波理论导出的位场向下延拓滤波器为最佳下延滤波器,但因其实现需要已知待求位场和噪声的功率谱而在实际应用中受到限制.针对该问题,本文首先提出一种基于位场径向平均功率谱的位场噪声水平估计方法,进而利用偏差准则求取正则化参数,实现位场正则化向下延拓;然后将位场正则化下延结果的功率谱作为待求位场功率谱的估计初值,采用带修正项的迭代维纳滤波方法来更新对待求位场功率谱的估计,最后提出本文的位场向下延拓改进迭代维纳滤波方法.基于理论重力模型数据及航磁实测数据进行了向下延拓对比试验,结果表明,改进迭代法具有较好的收敛性,且下延精度优于Tikhonov正则化法和递增型维纳滤波法.     

位场向下延拓的最小曲率方法

针对位场向下延拓的不适定性,我们将位场向下延拓视为向上延拓的反问题,提出以位场最小曲率作为约束条件来求解稳定的下延位场.我们将剖面位场向上延拓表达式用傅里叶矩阵的形式表示,以矩阵乘法形式给出延拓的表达式,同时向待反演的下延位场引入最小曲率约束,得到向下延拓的最小曲率解,并利用正交变换给出了更为简洁的频率域解.随后,利用Kronecker积将上述全部结果拓展至三维位场,给出了三维位场向下延拓的最小曲率解.此外,我们将位场数据的填充、扩充问题与向下延拓问题统筹考虑,提出一种新的向下延拓迭代格式,该算法面向实际资料处理需求、无须预扩充或填补数据.下延迭代时,对原始数据直接向下延拓,而空白部分利用上一次下延位场估计的上延值替代其空白值并对其向下延拓,直至获得最小曲率约束下稳定的向下延拓结果.同时,我们也讨论了利用改进L曲线和广义交叉验证(GCV)计算正则参数最优估计的问题.对理论模型和实际航空重力资料进行了向下延拓检验,处理结果表明位场向下延拓的最小曲率方法解能满足实际位场资料对向下延拓处理的需求,具有较高的下延精度.     

海域航空重力快速构建区域大地水准面

我国在海域开展了大规模的航空重力勘探,这些资料对构建高精度大地水准面具有重要价值.基于此,本文提出一种利用海域航空重力测量数据快速构建大地水准面的方法.该方法基于移去-恢复法思想,利用位场最小曲率方法对航空重力数据进行高精度向下延拓并获取相应的扰动位,实现航空重力测量快速构建海域大地水准面.与斯托克斯积分计算相比,采用了处理效率更高的频率域位场转换,解决了向下延拓及垂向积分时航空重力异常数据空白及扩边问题,具有较高的位场转换精度.本文应用EGM2008模拟航空重力数据进行模型验证,计算结果与其给出的水准面的精度相当;同时,也选取GRAV-D计划的航空重力数据进行实际验证,计算结果与xGEOID18B水准面模型精度基本一致.模型验证和实际应用验证了本方法的实用性.     

Structure of the Hellenic Subduction Zone from Gravity Gradient Functions and Seismology

The Normalized Full Gradient (NFG) method has widespread applications in the analysis of potential fields, especially the gravity and magnetic fields. This method is used to identify the lateral and horizontal density variations in the crust and lithosphere. In this study, the NFG method was applied to the gravity data of the Cretan Arc and its surroundings. Because of the tectonic features of the eastern Mediterranean, the Cretan Arc and the neighboring areas are seismically very active. Especially the subduction zone and the complicated crustal features have been defined applying many different geophysical methods. In this study, first the NFG method is tested with synthetic prisms (two cubes). After that, the NFG method was applied to the Bouguer gravity data of the Cretan Arc and its subduction zone (Hellenic subduction zone) and Hellenic subduction zone was defined with the foci depth data (USGS) along the south–north direction. Thus, geometry of the focal depth distribution has been created to determine probable media depths and their localizations. According to the NFG results, vertical structural transitions were observed at a depth ranging between 10 and 180?km. Also, these results were compared with the foci depth model and the other results of the related publications. Finally, some considerations in vertical solution with the NFG method have been presented and locations of the different structures at horizontally have been defined with application of the NFG method.     

基于向上延拓的航空重力向下解析延拓解

位场向下与向上延拓之间存在固有的内在联系,向上延拓解算具有稳定可靠的优良特性,本文据此提出了借助向上延拓信息实现航空重力向下延拓稳定解算的两种方法,分别建立了点对点向下解析延拓模型和最小二乘向下解析延拓模型.其核心思想是,依据泰勒级数展开模型,将位场向下延拓解算过程转换为向上延拓计算和垂向偏导数解算两个步骤,通过第一步的处理有效抑制数据观测噪声对解算结果的干扰,通过第二步的处理成功实现向下延拓反问题的稳定解算,较好地解决了向下延拓解算固有的不适定性问题.分析研究了两种解析延拓模型的计算精度及适用条件,利用超高阶位模型EGM2008建立的模拟标准场数据对两种模型解算结果的合理性和有效性进行了数值验证,证明本文新方法实用易行,具有较高的应用价值.     

积分迭代法的正则性分析及其最优步长的选择

位场积分迭代法是一种解决位场大数据量、大深度向下延拓的有效方法.本文基于Kirsch正则化子理论,推导了积分迭代法对应的正则化滤子函数,并证明积分迭代法为一种求解位场向下延拓不适定反问题的正则化方法.针对积分迭代法迭代步长固定、迭代次数较多影响收敛速度的问题,提出该迭代法最优迭代步长的选择原理.理论模型和实测数据对比分...     

Topographic and atmospheric effects on goce gradiometric data in a local north-oriented frame: A case study in Fennoscandia and Iran

Satellite gradiometry is an observation technique providing data that allow for evaluation of Stokes’ (geopotential) coefficients. This technique is capable of determining higher degrees/orders of the geopotential coefficients than can be achieved by traditional dynamic satellite geodesy. The satellite gradiometry data include topographic and atmospheric effects. By removing those effects, the satellite data becomes smoother and harmonic outside sea level and therefore more suitable for downward continuation to the Earth’s surface. For example, in this way one may determine a set of spherical harmonics of the gravity field that is harmonic in the exterior to sea level. This article deals with the above effects on the satellite gravity gradients in the local north-oriented frame. The conventional expressions of the gradients in this frame have a rather complicated form, depending on the first-and second-order derivatives of the associated Legendre functions, which contain singular factors when approaching the poles. On the contrary, we express the harmonic series of atmospheric and topographic effects as non-singular expressions. The theory is applied to the regions of Fennoscandia and Iran, where maps of such effects and their statistics are presented and discussed.     

AUTOMATIC DECONVOLUTION OF GRAVIMETRIC ANOMALIES*

Existing techniques of deconvolution of gravity anomalies are principally based on upward and downward continuation of measured fields. It can be shown that a unique set of linear filters, depending only on geometrical parameters, relates density distribution at a given depth to gravity measured on the surface. A method to compute the filter coefficients is developed. Very accurate reconstitution of theoretical models of intricate shape, prove the validity of the linear relationship. One of these sets of linear filters is applied to a field case of underground quarries.     

Conditionality of inverse solutions to discretised integral equations in geoid modelling from local gravity data

The eigenvalue decomposition technique is used for analysis of conditionality of two alternative solutions for a determination of the geoid from local gravity data. The first solution is based on the standard two-step approach utilising the inverse of the Abel-Poisson integral equation (downward continuation) and consequently the Stokes/Hotine integration (gravity inversion). The second solution is based on a single integral that combines the downward continuation and the gravity inversion in one integral equation. Extreme eigenvalues and corresponding condition numbers of matrix operators are investigated to compare the stability of inverse problems of the above-mentioned computational models. To preserve a dominantly diagonal structure of the matrices for inverse solutions, the horizontal positions of the parameterised solution on the geoid and of data points are identical. The numerical experiments using real data reveal that the direct gravity inversion is numerically more stable than the downward continuation procedure in the two-step approach.     

A refined method of recovering potential coefficients from surface gravity data

Summary A new method for computing the potential coefficients of the Earth's external gravity field is presented. The gravimetric boundary-value problem with a free boundary is reduced to the problem with a fixed known telluroid. The main idea of the derivation consists in a continuation of the quantities from the physical surface to the telluroid by means of Taylor's series expansion in such a way that the terms whose magnitudes are comparable with the accuracy of today's gravity measurements are retained. Thus not only linear, but also non-linear terms are taken into account. Explicitly, the terms up to the order of the third power of the Earth's flattening are retained. The non-linear boundary-value problem on the telluroid is solved by an iteration procedure with successive approximations. In each iteration step the solution of the non-linear problem is estimated by the solutions of two linear problems utilizing the fact that the non-linear boundary condition may be split into two parts; the linear spherical approximation of the gravity anomaly whose magnitude is significantly greater than the others and the non-linear ellipsoidal corrections. Finally, in order to solve the problem in terms of spherical harmonics, the transform method composed of the fast Fourier transform and Gauss Legendre quadrature is theoretically outlined. Immediate data processing of gravity data measured on the physical Earth's surface without any continuation of gravity measurements to a reference level surface belongs to the main advantage of the presented method. This implies that no preliminary data handling is needed and that the error data propagation is, consequently, maximally suppressed.     

The iteration method for downward continuation of a potential field from a horizontal plane

This paper introduces an iteration method for the downward continuation of potential‐field data from a horizontal plane, and compares it with the conventional frequency‐domain method (Fourier transform) using 2D and 3D model tests. The paper evaluates the two methods in terms of the results, i.e. downward‐continuation distance and stability. The iteration method proves to be more stable and able to downward continue the potential for a greater distance than the Fourier transform method.     

A study of the upward continuation of gravity data from a plane surface

Summary A summation method of upward continuation of gravity data has been considered under the assumption that observations are available at regular intervals. The upward continued value has been obtained as the sum of products of individual gravity values and corresponding theoretical coefficients. Besides the usual parameter involving horizontal and vertical distances, the theoretical coefficients have been generalized to be dependent also on i) the order of a low order polynomial assumed to represent the gravity variation around a grid point and ii) the weights assigned to the gravity values at the nearest four grid points used for least-squares determination of the polynomial. Since the observations in practical cases are available over a finite area only, the effect of truncation of the area of the integration has also been discussed separately. The method has been programmed and tested on a three-dimensional model, whose true gravity effects were computed at various levels over a finite area. Upward continued values have been computed under various assumptions about the gravity field in the outside region. Comparisons of these results with the true values indicate that the truncation effect becomes increasingly important for larger values of the ratio of elevation to grid separation and/or when the gravity field is not a local one. It has also been found that the greater is the above ratio, the less important is the effect of weights on the theoretical coefficients and practically vanishes(<10 ^–4 ) when the ratio is greater than5.0.     

基于多尺度边缘约束的重力场信号的向下延拓

利用多尺度边缘检测理论,研究改善重力场信号向下延拓稳定性的新方法,简要介绍了小波多尺度边缘检测的基本原理,以EGM96重力场模型为例分析了地球重力场的多尺边缘。提出了基于多尺度边缘约束的向下延拓方法,并给出了模拟算例。计算结果表明,利用多尺度边缘约束改善向下延拓的稳定性是可行的,在观测值信噪比较差的情况下该方法的效果尤为明显。