共查询到20条相似文献，搜索用时 203 毫秒
1.
Automatic Interpretation of Magnetic Data Using Euler Deconvolution with Nonlinear Background 总被引：1，自引：0，他引：1
Pawan Dewangan T. Ramprasad M. V. Ramana M. Desa B. Shailaja 《Pure and Applied Geophysics》2007,164(11):2359-2372
The voluminous gravity and magnetic data sets demand automatic interpretation techniques like Naudy, Euler and Werner deconvolution.
Of these techniques, the Euler deconvolution has become a popular choice because the method assumes no particular geological
model. However, the conventional approach to solving Euler equation requires tentative values of the structural index preventing
it from being fully automatic and assumes a constant background that can be easily violated if the singular points are close
to each other. We propose a possible solution to these problems by simultaneously estimating the source location, depth and
structural index assuming nonlinear background. The Euler equation is solved in a nonlinear fashion using the optimization
technique like conjugate gradient. This technique is applied to a published synthetic data set where the magnetic anomalies
were modeled for a complex assemblage of simple magnetic bodies. The results for close by singular points are superior to
those obtained by assuming linear background. We also applied the technique to a magnetic data set collected along the western
continental margin of India. The results are in agreement with the regional magnetic interpretation and the bathymetric expressions. 相似文献
2.
We have developed a least-squares minimization approach to determine the depth and the amplitude coefficient of a buried structure
from residual gravity anomaly profile. This approach is basically based on application of Werner deconvolution method to gravity
formulas due to spheres and cylinders, and solving a set of algebraic linear equations to estimate the two-model parameters.
The validity of this new method is demonstrated through studying and analyzing two synthetic gravity anomalies, using simulated
data generated from a known model with different random error components and a known statistical distribution. After being
theoretically proven, this approach was applied on two real field gravity anomalies from Cuba and Sweden. The agreement between
the results obtained by the proposed method and those obtained by other interpretation methods is good and comparable. Moreover,
the depth obtained by the proposed approach is found to be in very good agreement with that obtained from drilling information. 相似文献
3.
常规欧拉反褶积法中构造指数的选取以及分散解存在较多的问题,提出了基于联立垂向一阶导数与解析信号的欧拉齐次方程的RDAS-Euler反演方法。该方法可以更为精确的估计场源的范围及埋深,且不需考虑构造指数N的影响,避免了因构造指数不当而引起的反演误差。通过对单一地质体及组合地质体模型的实验证明本文方法能有效地完成目标体的反演工作,反演结果与理论值之间的误差小于10%,且相对于常规欧拉反褶积法更加稳定准确,能够更好的得到地质体边界及深度信息。将RDAS-Euler法应用于黑龙江省虎林盆地实测布格重力异常数据,获得了丰富的断裂信息,说明RDASEuler法增强了对断裂平面位置的识别能力。 相似文献
4.
G.R.J. Cooper 《Geophysical Prospecting》2014,62(5):1169-1179
This paper introduces the conversion of Euler's equation from a Cartesian coordinate system to a radial coordinate system, and then demonstrates that for sources of the type 1/rN (where r is the distance to the source, and N is the structural index) it can be solved at each point in space without the need for inversion, for a known structural index. It is shown that although the distance to the source that is obtained from Euler's equation depends on the structural index used, the direction to the source does not. For some models, such as the gravity and magnetic response of a contact, calculation of the analytic signal amplitude of the data is necessary prior to the application of the method. Effective noise attenuation strategies, such as the use of moving windows of data points, are also discussed. The method is applied to gravity and magnetic data from South Africa, and yields plausible results. 相似文献
5.
Multiridge Euler deconvolution 总被引：1，自引：0，他引：1
Potential field interpretation can be carried out using multiscale methods. This class of methods analyses a multiscale data set, which is built by upward continuation of the original data to a number of altitudes conveniently chosen. Euler deconvolution can be cast into this multiscale environment by analysing data along ridges of potential fields, e.g., at those points along lines across scales where the field or its horizontal or vertical derivative respectively is zero. Previous work has shown that Euler equations are notably simplified along any of these ridges. Since a given anomaly may generate one or more ridges we describe in this paper how Euler deconvolution may be used to jointly invert data along all of them, so performing a multiridge Euler deconvolution. The method enjoys the stable and high‐resolution properties of multiscale methods, due to the composite upward continuation/vertical differentiation filter used. Such a physically‐based field transformation can have a positive effect on reducing both high‐wavenumber noise and interference or regional field effects. Multiridge Euler deconvolution can also be applied to the modulus of an analytic signal, gravity/magnetic gradient tensor components or Hilbert transform components. The advantages of using multiridge Euler deconvolution compared to single ridge Euler deconvolution include improved solution clustering, increased number of solutions, improvement of accuracy of the results obtainable from some types of ridges and greater ease in the selection of ridges to invert. The multiscale approach is particularly well suited to deal with non‐ideal sources. In these cases, our strategy is to find the optimal combination of upward continuation altitude range and data differentiation order, such that the field could be sensed as approximately homogeneous and then characterized by a structural index close to an integer value. This allows us to estimate depths related to the top or the centre of the structure. 相似文献
6.
A discussion is given of the requirements, the advantages and the methods to be considered in attempting the quantitative interpretation of gravity or magnetic fields from computed maps of the vertically derived field. The transform which is used here is the first vertical derivative (or vertical gradient) with or without downward continuation, but the computed maps are in fact obtained by a controllable Fourier method in which two kinds of operations can be simultaneously performed in complete independance: on one hand the separation of any part of the data by frequency cut-off, and on the other hand the transformation by vertical derivation or continuation of the part which is retained. Taking as raw data either actual surveys or artificially constructed maps, it is first shown how separated and transformed maps of this type can be efficiently obtained under quite flexible conditions, using a special computer program. It is further seen that for correctly controlled filterings the accuracy of the computed maps actually permits to take them as the basis for quantitative interpretation. To effect this, any one of the conventional methods which make use of equivalent model computations may in the first place be adapted to the interpretation of gradients, with the benefits, however, of an enhanced lateral separation of the anomalies and of a large attenuation of the regional effects. Particularly, the delineation of horizontal contours for even fairly complex models can often be made directly in a sufficiently safe way on the anomalies as they show on the vertical gradient maps. This greatly accelerates the process of determining equivalent model bodies. More special methods of interpretation can also be designed by taking into account first the fact that the vertical derivation of the field amounts to an operation of separation on the field's sources themselves, and in addition the availability of the frequency form of the information as a result of using a Fourier method of transform computation. Trial utilisations of various such interpretation processes of either the conventional or the less conventional type, are presented especially in connection with an interpretation study on transformed maps of parts of an offshore aeromagnetic survey (English Channel). 相似文献
7.
In this paper, I introduce a new approach based on truncated singular value decomposition (TSVD) analysis for improving implementation of grid-based Euler deconvolution with constraints of quasi 2D magnetic sources. I will show that by using TSVD analysis of the gradient matrix of magnetic field anomaly (reduced to pole) for data points located within a square window centered at the maximum of the analytic signal amplitude, we are able to estimate the strike direction and dip angle of 2D structures from the acquired eigenvectors. It is also shown that implementation of the standard grid-based Euler deconvolution can be considerably improved by solving the Euler's homogeneity equation for source location and structural index, simultaneously, using the TSVD method. The dimensionality of the magnetic anomalies can be indicated from the ratio between the smallest and intermediate eigenvalues acquired from the TSVD analysis of the gradient matrix. For 2D magnetic sources, the uncertainty of the estimated source location and structural index is significantly reduced by truncating the smallest eigenvalue.Application of the method is demonstrated on an aeromagnetic data set from the Åsele area in Sweden. The geology of this area is dominated by several dike swarms. For these dolerite dikes, the introduced method has provided useful information of strike directions and dip angles in addition to the estimated source location and structural index. 相似文献
8.
We present a new integrated approach to the interpretation of magnetic basement that is based on recognition of characteristic patterns in distributions and alignments of magnetic source depth solutions above and below the surface of magnetic basement. This approach integrates a quantitative analysis of depth solutions, obtained by 2D Werner deconvolution of the magnetic data, with a qualitative evaluation of the Bouguer gravity anomalies. The crystalline/metamorphic basement and sedimentary cover have different origins, tectonic histories, lithologies and magnetic properties. These differences result in different geometries of magnetic sources associated with faults, fracture zones, igneous intrusions, erosional truncations, subcrop edges and other structural discontinuities. Properly tuned, 2D Werner deconvolution is able to resolve the intra‐sedimentary and intra‐basement magnetic source geometries into distinctly different distributions and alignments of calculated depth solutions. An empirical set of criteria, basement indicators, was developed for identification and correlation of the basement surface. The ambiguity of basement correlation with limited or non‐existent well control, which is common for onshore frontier and offshore explorations, can be reduced by incorporating the Bouguer gravity data into the process of correlation. 相似文献
9.
Structure and State of Stress of the Chilean Subduction Zone from Terrestrial and Satellite-Derived Gravity and Gravity Gradient Data 总被引：1，自引：1，他引：0
B. D. Gutknecht H.-J. Götze T. Jahr G. Jentzsch R. Mahatsente St. Zeumann 《Surveys in Geophysics》2014,35(6):1417-1440
It is well known that the quality of gravity modelling of the Earth’s lithosphere is heavily dependent on the limited number of available terrestrial gravity data. More recently, however, interest has grown within the geoscientific community to utilise the homogeneously measured satellite gravity and gravity gradient data for lithospheric scale modelling. Here, we present an interdisciplinary approach to determine the state of stress and rate of deformation in the Central Andean subduction system. We employed gravity data from terrestrial, satellite-based and combined sources using multiple methods to constrain stress, strain and gravitational potential energy (GPE). Well-constrained 3D density models, which were partly optimised using the combined regional gravity model IMOSAGA01C (Hosse et al. in Surv Geophys, 2014, this issue), were used as bases for the computation of stress anomalies on the top of the subducting oceanic Nazca plate and GPE relative to the base of the lithosphere. The geometries and physical parameters of the 3D density models were used for the computation of stresses and uplift rates in the dynamic modelling. The stress distributions, as derived from the static and dynamic modelling, reveal distinct positive anomalies of up to 80 MPa along the coastal Jurassic batholith belt. The anomalies correlate well with major seismicity in the shallow parts of the subduction system. Moreover, the pattern of stress distributions in the Andean convergent zone varies both along the north–south and west–east directions, suggesting that the continental fore-arc is highly segmented. Estimates of GPE show that the high Central Andes might be in a state of horizontal deviatoric tension. Models of gravity gradients from the Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite mission were used to compute Bouguer-like gradient anomalies at 8 km above sea level. The analysis suggests that data from GOCE add significant value to the interpretation of lithospheric structures, given that the appropriate topographic correction is applied. 相似文献
10.
以等效源及位场物性反演为基础,本文提出一种新的求取重磁场源深度的方法.该方法将一层等效源以一定的间隔从浅部向深部移动,并将等效源作为初始模型进行反演,当反演拟合差最小时,停止反演,此时的等效源底深即为所求场源的中心深度.由于仅需要反演一层等效源,比传统的物性反演计算时间大大减少,并且不需要进行深度加权约束.理论模型数据处理结果表明该方法能够获得较准确的场源深度:以长宽比为7.5的薄板模型为例,深度计算误差约为1个点距(25 m);以长宽比为0.5~1.5的厚板模型为例,深度计算误差小于1个点距(25m).将该方法应用于实测航磁梯度数据,计算的磁源中心深度在200~250m之间,钻井资料显示该异常由埋藏深度在200~300m的闪长岩引起,计算结果与钻井资料较吻合. 相似文献
11.
位场全张量梯度数据以其信息量大、含有更高频的信号成分,能更好地描述小的异常特征等优点在地球物理领域中得到广泛应用.边界检测是位场解释中不可缺少的任务,需要新的边界探测器来处理位场梯度张量数据.为了充分利用位场梯度张量数据的多信息成分,本文定义了方向总水平导数和加强方向总水平导数,并利用其定义新的边界检测器.为了能同时显示不同振幅大小异常的边界,本文对其进行了归一化处理.通过模型试验,证明了归一化方法能更加清晰准确地显示浅部和深部的地质体边界信息.最后将该边界检测方法用于加拿大圣乔治湾实际测得全张量重力梯度数据和中国朱日和地区的磁异常数据中,并得到了较好的边界检测结果. 相似文献
12.
Aspects of the interpretation of measured data on the gravity gradient tensor (GGT) are examined. The problem is posed in relation to the great progress achieved in recent years in the development of instrumentation and the method of GGT measurements on mobile carriers. In our opinion, the new methods of measurement and the new data obtained with their help require the development of new methods of interpretation of potential fields. The paper addresses two methods taking advantage of simultaneous measurements of all components of the GGT and anomalies of the gravitaty field V z. It is shown that the joint analysis of all GGT components can provide independent constraints on the noise level in various components. The method of tensor deconvolution proposed in the paper is a tensor analogue of the Euler method. The method is based on the calculation of invariants and is, therefore, stable with respect to the orientation uncertainties of the measuring system. The method provides means for estimating the structural index and, therefore, is particularly effective in the treatment of fields that contain isometric and/or elongated anomalies. The calculation of invariants and the tensor ratio can also be used for the development of procedures enabling automatic estimation of the axis strike azimuths of elongated anomaly-forming bodies. 相似文献
13.
地质体构造边界位置的确定是位场数据解释中的一项重要工作,现有很多基于位场梯度张量数据的边界检测滤波器,但存在识别边界位置模糊且无法均衡深浅地质体异常的缺点.本文定义了位场数据的三维构造张量,并提出基于位场构造张量的边界滤波器.为了同时显示不同振幅异常的边界位置,对新定义的滤波器进行归一化处理.在高阶均衡滤波器的计算中需要计算位场的垂向高阶导数,本文引入一种计算的稳定算法,基于拉普拉斯方程利用位场水平导数求解垂向导数,可减小垂向导数计算中产生的误差.将定义的滤波器应用到合成的重磁数据中证明了新方法相比传统的滤波器能更加清晰、准确地圈定边界位置,而且针对同时含有正负异常的地质情况,可避免产生额外的错误边界.最后将新的滤波器应用到实测的重磁数据的解释中,结果显示基于构造张量的滤波器可更准确清晰地划分出断裂的边界位置,发现更多的构造细节. 相似文献
14.
Non‐linear stochastic inversion of gravity data via quantum‐behaved particle swarm optimisation: application to Eurasia–Arabia collision zone (Zagros,Iran) 
下载免费PDF全文
下载免费PDF全文 Potential field data such as geoid and gravity anomalies are globally available and offer valuable information about the Earth's lithosphere especially in areas where seismic data coverage is sparse. For instance, non‐linear inversion of Bouguer anomalies could be used to estimate the crustal structures including variations of the crustal density and of the depth of the crust–mantle boundary, that is, Moho. However, due to non‐linearity of this inverse problem, classical inversion methods would fail whenever there is no reliable initial model. Swarm intelligence algorithms, such as particle swarm optimisation, are a promising alternative to classical inversion methods because the quality of their solutions does not depend on the initial model; they do not use the derivatives of the objective function, hence allowing the use of L1 norm; and finally, they are global search methods, meaning, the problem could be non‐convex. In this paper, quantum‐behaved particle swarm, a probabilistic swarm intelligence‐like algorithm, is used to solve the non‐linear gravity inverse problem. The method is first successfully tested on a realistic synthetic crustal model with a linear vertical density gradient and lateral density and depth variations at the base of crust in the presence of white Gaussian noise. Then, it is applied to the EIGEN 6c4, a combined global gravity model, to estimate the depth to the base of the crust and the mean density contrast between the crust and the upper‐mantle lithosphere in the Eurasia–Arabia continental collision zone along a 400 km profile crossing the Zagros Mountains (Iran). The results agree well with previously published works including both seismic and potential field studies. 相似文献
15.
Hasan Yildiz 《Studia Geophysica et Geodaetica》2012,56(1):171-184
Gravity field and steady-state Ocean Circulation Explorer (GOCE) is the first satellite mission that observes gravity gradients
from the space, to be primarily used for the determination of high precision global gravity field models. However, the GOCE
gradients, having a dense data distribution, may potentially provide better predictions of the regional gravity field than
those obtained using a spherical harmonic Earth Geopotential Model (EGM). This is investigated in Auvergne test area using
Least Squares Collocation (LSC) with GOCE vertical gravity gradient anomalies (Tzz), removing the long wavelength part from EGM2008 and the short wavelength part by residual terrain modelling (RTM). The results
show that terrain effects on the vertical gravity gradient are significant at satellite altitude, reaching a level of 0.11
E?tv?s unit (E.U.) in the mountainous areas. Removing the RTM effects from GOCE Tzz leads to significant improvements on the LSC predictions of surface gravity anomalies and quasigeoid heights. Comparison
with ground truth data shows that using LSC surface free air gravity anomalies and quasi-geoid heights are recovered from
GOCE Tzz with standard deviations of 11 mGal and 18 cm, which is better than those obtained by using GOCE EGMs, demonstrating that
information beyond the maximal degree of the GOCE EGMs is present. Investigation of using covariance functions created separately
from GOCE Tzz and terrestrial free air gravity anomalies, suggests that both covariance functions give almost identical predictions. However,
using covariance function obtained from GOCE Tzz has the effect that the predicted formal average error estimates are considerably larger than the standard deviations of
predicted minus observed gravity anomalies. Therefore, GOCE Tzz should be used with caution to determine the covariance functions in areas where surface gravity anomalies are not available,
if error estimates are needed. 相似文献
16.
Dr. Raoul Vajk 《Pure and Applied Geophysics》1951,19(3-4):129-143
Summary Methods for the elimination of the regional effect from gravity data are discussed in this paper. A historical review of the average gradient method used in torsion balance data is given. The connection between the averaging, or grid method, used on the gravity meter data and the average gradient method is pointed out and a criticism of these methods is given. Due to the ambiguity of the interpretation of gravity data no unique determination of the regional effect is possible and no mechanical method for its removal can be found. The determination of a regional gravity effect is an interpretation problem and for this reason it is unavoidably the subject of geological consideration and of the personal judgement of the interpreter. The method of smoothing the isogams and gravity curves is also discussed. A few suggestions are made for breaking the gravity data into two components (regional and residual) by the analysis of the gravity data. Finally the application of the second derivative methods for the elimination of the regional effect is discussed. 相似文献
17.
Gregory N. Tsokas 《Pure and Applied Geophysics》1995,144(1):117-134
The anomalies which are shown on the gravity and aeromagnetic data of the Chalkidhiki peninsula of Northern Greece reveal
the presence of three concealed dense and magnetized bodies. The ophiolitic belt which continues northwards of the geophysical
anomalies suggests that they are caused by structures which belong to it.
Multiple source Werner deconvolution shows that the lower surfaces of buried bodies in Toroneos Gulf and the one onshore are
concave upwards.
The northern concealed body has a width of 16 kilometers, while the southern one is 8 kilometres in width. A three-dimensional
model, well constrained, gives an estimate of the dimensions of the onshore feature.
The edges of the bodies, delineated by the terracing method, are formed by faults along the trend of the Hellenides orogenic
belt. The edges of the body in the outer part of Toroneos Gulf may be associated with seismic swarms that occurred in the
1980s. 相似文献
18.
本文提出了一种全张量磁梯度数据的斜导数边界识别方法,该方法利用全张量磁梯度数据定义了北向斜导数、东向斜导数和垂向斜导数.模型试验表明,北向和东向斜导数无法有效识别磁性异常体边界,异常形态复杂虚假异常较多,而垂向斜导数在高纬度斜磁化或者垂直磁化条件下可以清晰而准确得识别地质体的边界,具有一定的深度识别能力.与现有磁张量梯度的边界识别方法相比,识别效果好,可以有效均衡不同深度地质体的响应.将上述方法应用于大兴安岭地区实测航磁三分量数据转换得到的全张量数据,获得了研究区浅部磁性体的分布,并与三分量数据进行联合解释,取得良好的地质效果. 相似文献
19.
边界识别技术是位场数据解释中一项基本的工作,现有的边界识别方法多存在边界识别结果发散和不能均衡深浅地质体异常的缺点.目前一些均衡边界识别方法会因正负异常同时存在而引起额外的错误边界或者存在人为主观因素去除错误边界信息的缺点.本文充分利用重力位场张量梯度的多信息成分,提出利用水平方向解析信号及其垂向导数与传统的均衡边界识别方法做结合的方式定义新的探测边界的方法.通过理论模型试验证明新方法同传统方法相比,能够更加清晰、准确的圈定出深浅地质体的边界.最后将新的边界识别方法应用到实测重力异常数据解释中,取得了良好的边界识别结果并能够发现更多的构造细节. 相似文献
20.
全张量测量技术是在空中或海上用加载了多个加速度计的移动平台技术测量位场的五个独立分量.各张量分量包含不同方向的地下地质体信息,水平张量分量T_(xx)、T_(yy)、T_(xy)、T_(xz)、T_(yz)通常用于识别和映射与地质构造或地层变化有关的测量区域中的目标,垂直张量分量Tzz用于估计深度.然而,这些分量传统上是彼此分开解释,经常遇到错失关键信息的风险.本文所用全张量欧拉反褶积是在单独z方向的欧拉反演基础上发展而来的,它融合了重力异常垂直分量以及其三个方向导数、水平分量以及其三个方向导数.全张量数据信息得以有效应用的同时,欧拉反褶积结果也比常规欧拉反褶积结果更加收敛.最后,结合美国墨西哥湾地区实测航空FTG数据,用重力梯度张量数据进行联合欧拉三维反演研究,有效的识别岩盖的边界信息,划分岩盖范围,为进一步研究盖层底下深部复杂地质情况提供可靠的解释结果. 相似文献