共查询到20条相似文献，搜索用时 15 毫秒
1.
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. 相似文献
2.
A new best estimate methodology is proposed and oriented towards the determination of parameters related to a magnetic field
anomaly produced by a simple geometric-shaped model or body such as a thin dike and horizontal cylinder. This approach is
mainly based on solving a system of algebraic linear equations for estimating the three model parameters, e.g., the depth
to the top (center) of the body (z), the index parameter or the effective magnetization angle (θ) and the amplitude coefficient or the effective magnetization intensity (k). The utility and validity of this method is demonstrated by analyzing two synthetic magnetic anomalies, using simulated
data generated from a known model with different random errors components and a known statistical distribution. This approach
was also examined and applied to two real field magnetic anomalies from the United States and Brazil. 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 such an approach is found to be in high accordance with that obtained from drilling information. The
advantages of such a proposed method over other existing interpretative techniques are clarified, where it can be generalized
to be automatically applicable for interpreting other geological structures described by mathematical formulations. 相似文献
3.
A new method is proposed to interpret magnetic anomalies due to a thin dike, a sphere, and a fault like structure, where depth, horizontal location, effective magnetization intensity and effective magnetization inclination of a buried structure are simultaneously obtained. The proposed method is based on Fair function minimization and also on stochastic optimization modeling. This new technique was firstly tested on a theoretical synthetic data randomly generated by a chosen statistical distribution from a known model with different random noises components. This mathematical simulation shows a very close agreement between the assumed and the estimated parameters. The applicability and validity of this method are thereafter applied to magnetic anomaly data taken from United States, Australia, India, and Brazil. The agreement between the results obtained by the new method and those obtained by other interpretative methods is good and comparable. Moreover, the depth obtained by such a method is found to be in high accordance with that obtained from drilling information. 相似文献
4.
A new approach is proposed in order to interpret field self-potential (SP) anomalies related to simple geometric-shaped models such as sphere, horizontal cylinder, and vertical cylinder. This approach
is mainly based on solving a set of algebraic linear equations, and directed towards the best estimate of the three model
parameters, e.g., electric dipole moment, depth, and polarization angle. Its utility and validity are demonstrated through
studying and analyzing synthetic self-potential anomalies obtained by using simulated data generated from a known model and
a statistical distribution with different random errors components. Being theoretically tested and proven, this approach has
been consequently applied on two real field self-potential anomalies taken from Colorado and Turkey. A comparable and acceptable
agreement is obtained between the results derived by the new proposed method and those deduced by other interpretation methods.
Moreover, the depth obtained by such an approach is found to be very close to that obtained by drilling information. 相似文献
5.
Giovanna Berrino Gennaro Corrado Umberto Riccardi 《Journal of Volcanology and Geothermal Research》2008
Bathymetric and gravity surveys were carried out from 1988 to 1994, in the Gulf of Naples (Southern Italy) to offshore extend the already existing Bouguer anomaly map. In order to improve the knowledge of the structural setting beneath the active Neapolitan volcanoes (Vesuvio, Campi Flegrei and Ischia), 862 stations were surveyed within the isobath of 400 m; at the same time, and about 2000 on-land gravity values were also collected. A new Bouguer anomaly map spanning the whole volcanic region was drawn from the final data set. Gravity anomalies were referred to the new absolute gravity station in Naples and computed according to 1980 Geodetic Reference System. Finally, a density value of 2200 kg/m3 was used in the computation of the Bouguer and terrain effects. We carried out the inversion of the gravity anomalies adopting a 2.5-D modelling along selected profiles crossing the investigated area. The interpretative models were constrained to data obtained from deep wells and other geophysical investigations. 相似文献
6.
重力异常对地壳横向密度变化敏感,而无约束重力反演得到的密度模型其垂向分辨能力往往不理想.为了改善反演结果的垂向分辨率,本文参考已有先验分层模型,基于贝叶斯原理,提出了一种重震联合反演的新策略,可实现多种参考模型和复杂加权参数条件下的最大后验概率估计.理论模型测试结果表明,对于深度加权、多参考模型约束等多种问题,本文提出的新方法都可以稳健地获得最优化的模型参数.本文同时以中国地震科学台阵在龙门山地区及周边的一维接收函数分层模型和地震层析成像结果为参考,通过此方法对该区的重力异常进行反演,获得了该区的高精度三维密度结构,其水平分辨率优于10 km,垂直分辨率优于5 km.结合四条通过汶川和芦山地震震中的剖面进行分析后发现,反演得到的密度结构模型在过强震震源区位置横向变形显著,其揭示的分层地壳结构和变形模式与地表已知断裂构造具有相关性.本文提出的重震联合反演新策略,可为研究潜在强震风险源区的地壳结构和物性特征提供有效的科技方法支撑. 相似文献
7.
A geophysical interpretative method is proposed to depth, amplitude coefficient and geometrical shape factor determination
of a buried structure from an observed gravity anomaly related to a cylinder or a sphere-like structure.The method is based
on nonlinearly constrained mathematical modelling and also on stochastic optimization approaches. The proposed interpretative
method first has been tested on theoretical synthetic models with different random errors at a certain depth, where a very
close agreement has been observed between assumed and evaluated parameters. Subsequent field data have been considered for
which the interpreted results by other methods are available for comparison. The agreement between the obtained results by
the proposed technique and by other geophysical methods is good. A statistical analysis has been also carried out to demonstrate
the accuracy and the precision of the suggested interpretative method. 相似文献
8.
Ayd?n Büyüksaraç 《Physics of the Earth and Planetary Interiors》2007,160(1):86-95
The residual aeromagnetic and gravity anomalies of inner East Anatolia, surveyed by the Mineral Research and Exploration (MTA) of Turkey, display complexities. Some faults, which are known and new lineaments, are drawn from maxspot map derived from the location of the horizontal gradient of gravity anomalies. Tectonic lineaments of inner East Anatolia exhibit similarities to the direction of East Anatolian Fault Zone. Anticlockwise rotation, approximately −30°, defined from disorientations of aeromagnetic anomalies. The lineaments obtained from maxspots map produced from the gravity anomalies and disoriented aeromagnetic anomalies are in-line with the mobilistic system revealed by the palaeomagnetic data. These Alpine age continental rotations caused westward wrenching of the global lithosphere and led to significant tectonic reactivation and deformations. GPS measurements, current tectonic knowledge and the results of the evaluation of potential field data were combined in a base map to demonstrate similarities. 相似文献
9.
RADU BOTEZATU 《Geophysical Prospecting》1970,18(Z1):800-815
One of the most important aspects in interpretation of gravity and magnetic anomalies, which generally represent the sum of many effects from the subsoil, is that these effects cannot be adequately separated into individual effective anomalies by means of known procedures. Moreover, the concept itself of anomaly separation must be expanded in a certain sense, because the addition of these effects in the mapped anomalies can arise from many structures situated in the subsoil in different relative positions to one another. A combined procedure is proposed, based on Fourier analysis utilising both periodic complex functions and nonperiodic ones, by use of which the number of individual effective anomalies cumulated in the mapped anomalies is obtained. In order to illustrate the use of this procedure and the advantages of obtaining more valuable informations on geological causes of anomalies, one example from Romania is given. 相似文献
10.
A new and simple method based on a nonlinearly mathematical optimization concept has been proposed in this research to interpret magnetic anomalies due to vertical faults and thin dikes. This proposed interpretative method consists of three main steps. The first step is to formulate nonlinearly constrained optimization problems to describe the geophysical problems related to the studied structures. The second step is to suggest an interior penalty function in order to convert these nonlinearly constrained optimization problems into nonlinearly unconstrained optimization ones. The third step is to solve the converted nonlinearly unconstrained optimization problems by using the famous Hooke and Jeevess algorithm in order to estimate the geophysical parameters of the studied structures such as: depth, amplitude coefficient, and index parameter. The Hooke and Jeevess algorithm is purposely chosen for being robust and also its application to magnetic data converges rapidly towards the optimal estimation of parameters. This method was first tested on theoretical models with different random noise, where a very close agreement was obtained between the assumed and evaluated parameters. The validity of this new method was also tested on practical field examples taken from Australia, India, United States, and Brazil, where available magnetic data existed and was previously analyzed by different interpretative methods. The agreement between the results obtained by our developed method and those obtained by the other geophysical methods is good. The advantages of this newly proposed method, compared with the other published interpretative methods, also have been discussed and demonstrated. 相似文献
11.
Structural interpretations of negative gravity anomalies at Campi Flegrei and Campania are reviewed. The interpretations are not uniquely determined, as different models correspond to equally good anomaly fits.These models are examined using a new method based on the interpretation of the apparent density map computed from the gravity anomaly map.As regards the Campi Flegrei anomaly, a funnel-shaped caldera 3D model is readily obtained. This is more appropriate than previous interpretations, as not only anomaly but also apparent densities are fitted.As about the regional structures, the used approach suggests that the more pertinent model is a large magmatic intrusion, 9 to 14 km deep.A final scheme is proposed linking this intrusion to the funnel-shaped caldera as the main regional reservoir. 相似文献
12.
Regional–residual separation is essential in gravity and magnetic data interpretation and a variety of techniques have been proposed. Graphical determination of the regional allows geological information to be taken into account. Upward continuation can be used to obtain the regional field either empirically or using some hypothesis about the geology. In some cases, a matched filter can be designed and used to separate deep and shallow sources. Simple low pass filtering has also been used but without much success. Here we propose to use a non-linear filter approach to remove gravity and magnetic anomalies smaller than a given width. This technique attempts to mimic the graphical separation method. The results from synthetic models are presented as well as the results from a case study in eastern Canada and compared to regional gravity and magnetic anomalies obtained by other techniques. Contrary to the regional fields obtained by upward continuation, non-linear filtering does not have any physical meaning. However, its main advantage is that it gives a regional component of the gravity or magnetic field similar to the one obtained from a graphical separation. 相似文献
13.
WAVELET TRANSFORM ANALYSES OF FAULTS DETECTION ON ISOSTATIC GRAVITY ANOMALIES: A CASE STUDY FROM THE QAIDAM BASIN AND ITS ADJACENT AREAS 
下载免费PDF全文
下载免费PDF全文 Isostasy is used to describe a condition to which the Earth's crust and mantle tend, in the absence of disturbing forces. Eliminating the gravity effect of crust, isostatic gravity anomalies contain abundant geological structure information, which can be extracted by edge detection methods of gravity or magnetic anomalies. In order to accurately obtain the edge information, a great variety of methods, such as analytical signal amplitude, tilt angle, theta map θ, etc., have been proposed by domestic and international scholars, and many significant advances have been made in recent days. However, each method has its advantages and disadvantages. Wavelet transform is an effective method developed in recent years. It has the enhanced noise resistance and a feature of multi-scale decomposition, and can be used to identify more detailed information of edge. Here with the aim of demonstrating its effectiveness in faults detection, we established a theoretical geological model, which consists of five geological bodies. The geological bodies with different density present a fault zone and the areas on its sides, as well as two geological bodies with different geological properties. We calculated the gravity anomalies caused by this model, in addition, we added 5% Gaussian noise to the gravity anomalies for a comparative analysis to analyze the effects of wavelet transform on edge detection. Finally, we applied wavelet transform method to the decomposition of isostatic gravity anomalies, obtained 1st to 5th order wavelet transform details of the gravity anomalies, and compared with the well-studied faults in the Qaidam Basin and its adjacent areas. The results obtained by wavelet transform matched well with the known faults, the anomalies of different order denote the location of different fault zone(e.g. the Tanan Fault is nearly invisible on the original and the first order isostasy gravity anomalies map, but is well expressed on the second order isostasy gravity anomalies map; The apparent details of the 4th and 5th indicate that faults in front of the Saishiteng-Xitieshan Shan are deep faults and they are likely to distribute continuously in the deep underground). Besides, we calculated the estimated depth of isostasy gravity anomalies of different order through power spectrum analysis as well, finding that different faults extend to different depth. For example, the Danghe-Nanshan Fault and the Southern Fault in the middle Qilian Shan are 10km in depth approximately, but the faults in front of the Saishiteng-Xitieshan Shan are more than 70km in depth. In addition, we made two comparative studies, the first one is comparing the results mentioned above with the result through wavelet transform of Bouguer gravity anomalies. The second one is comparing with the results through other edge detection methods of isostatic gravity anomalies. In spite of the inconformity between anomalies and the faults to some extent, which is likely caused by the change of lithology or faults distribution in the deep underground, we finally found that:more subtle details induced from faults can be detected from isostatic gravity anomalies by using wavelet transform because of its feature of multi-scale decomposition. The wavelet transform method is proved to be more accurate and reliable(at least in the Qaidam Basin and its adjacent areas)comparing with other methods. 相似文献
14.
V字型密度界面是一类常见的密度界面,如海沟、半地堑以及俯冲带之下的莫霍面,利用重力数据刻画此类密度界面形态对于区域构造研究、油气勘探以及物理海洋学等都具有重要意义.本文首先建立了Lp-范数形式的模型约束函数,并利用正则化原理将其与重力数据误差函数和已知深度约束函数结合形成V字型密度界面反演的目标函数,推导了目标函数的梯度表达式,并以非线性共轭梯度法为核心给出了反演流程.二维简单模型试算结果表明p=5时该方法能准确地刻画V字型密度界面起伏特征,且亦能准确地应用于二维复杂密度界面和三维界面的反演.最后将反演方法应用于挑战者深渊及邻区的实际资料处理之中,利用研究区海底地形数据和沉积层厚度数据对自由空间重力异常逐层剥离而得到莫霍面引起的重力异常,用本文方法对此重力异常进行反演,结果呈现了板块俯冲作用引起的V字型莫霍面起伏特征. 相似文献
15.
A new method is introduced to determine the top and bottom depth of a vertical line using gravity anomalies. For this, gravity at a distance x from the origin and horizontal derivative at that point are utilized. A numerical value is obtained dividing the gravity at point x by horizontal derivative. Then a new equation is obtained dividing the theoretical gravity equation by the derivative equation. In that equation, assigning various values to the depth and length of vertical line, several new numerical values are obtained. Among these values, a curve is obtained for the one that is closest to the first value from attending the depth and length values. The intersection point of these curves obtained by repeating this procedure several times for different points x yield the real depth and length values of the line. The method is tested on two synthetics and field examples. Successful results are obtained in both applications. 相似文献
16.
Transformation between gravimetric and GPS/levelling-derived geoids using additional gravity information 总被引：1，自引：0，他引：1
The transformation from the gravimetric to the GPS/levelling-derived geoid using additional gravity information for the covariance function of geoid height differences has been investigated in a test area in south-western Canada. A “corrector surface” model, which accounts for datum inconsistencies, long-wavelength geoid errors, vertical network distortions and GPS errors, has been constructed using least-squares collocation. The local covariance function of geoid height differences is usually obtained from residual values between the GPS/levelling and gravimetric geoid heights after the elimination of all known systematic distortions. If additional gravity data (in the form of gravity anomalies) are available, the covariance function of geoid height differences can be determined by the following steps: (1) transforming the GPS/levelling-derived geoid heights into gravity anomalies; (2) forming differences between the computed in step 1 and given gravity anomalies; (3) determining the parameters of the local covariance function of the gravity anomaly differences; (4) constructing an analytical covariance model for the geoid height differences from the covariance function of the gravity anomaly differences using the parameters derived in step 3. The advantage of the proposed method stems from the great number of gravity data used to derive the empirical covariance function. A comparison with the least-squares adjustment shows that the standard deviation of the residuals of the predicted geoid height differences with respect to the control point values decreases by 2.4 cm. 相似文献
17.
Mean free-air gravity anomalies are often needed in geodesy for gravity field modelling. Two possible ways of compiling the mean free-air gravity anomalies are discussed. One way is via simple Bouguer gravity anomalies and the second and more time consuming way is via refined Bouguer gravity anomalies. In flat areas the differences between using any of the two ways should not be significant. In the mountains however, every effect introducing a high dependency, such as e.g. terrain effect, can negatively affect the interpolation process. In fact, a numerical experiment conducted in one part of Rocky Mountains revealed large and systematic differences. The effect of these differences on the geoid model is more then two meters in the test area. Our investigation shows that this bias is caused by the location of gravity measurement points, chosen mostly on hill-tops. At such points, the terrain correction to gravity is systematically larger than the mean value of the correction. Therefore, it is not possible to prevent the mean free-air gravity anomalies obtained from simple Bouguer gravity anomalies from having a systematic bias. One can see this bias as a result of the aliasing effect because the simple Bouguer gravity anomalies in the mountains contain a higher frequency signal (terrain effect) that is, according to sampling theorem, impossible to reconstruct by sparse measured gravity data, see e.g. (Goos et al., 2003). Therefore, the more rigorous way of computing the mean free-air gravity anomalies is via refined Bouguer gravity anomalies. 相似文献
18.
基于全球EGM2008自由空气重力异常模型,本文计算了蒙古及周边地区的布格重力异常和AiryHeiskanen均衡重力异常.在此基础上,本文采用Crust 1.0地壳模型为参考,通过重力正演方法,对蒙古及周边地区不同深度地壳密度结构模型的重力异常进行了计算,并对得到的正演布格重力异常与实际重力异常进行了对比和分析.研究结果表明:蒙古西部杭爱山地区与阿尔泰山地区的构造变形差异性明显,现今均衡重力异常中杭爱山周边没有明显的均衡异常高值区,而阿尔泰山地区西南方向存在均衡重力异常高值分布,分析与新构造运动密切相关;Crust 1.0模型给出的壳幔横向密度不均匀体分布对于计算Moho面起伏引起的重力异常作用明显;Crust 1.0给出的地壳内界面变形可以反映深大活动断裂的深部构造变形.研究结果对于认识蒙古东西部构造特征差异,以及现今西部活动断裂的地球物理场特征具有参考意义,也可以为进一步应用Crust 1.0模型为参考开展三维密度结构反演提供一定帮助. 相似文献
19.
A new technique for delineation of subsurface features using gravity eigenimage is proposed to separate the gravitational anomalies from its background. The singular value decomposition and multifractal method has been combined and tested on several synthetic gravity anomalies and also applied to field gravity dataset of Vindhyan basin in central India which is a promising area for petroleum exploration. The eigenimage of gravity data helps to understand the relationship between the geological structure and source of anomaly. Fault structure is the major structure and could be the principal contributor that differentiates the regional and local geo-anomalies in the studied region. 相似文献
20.
Gravity anomaly reference fields, required e.g. in remove-compute-restore (RCR) geoid computation, are obtained from global
geopotential models (GGM) through harmonic synthesis. Usually, the gravity anomalies are computed as point values or area
mean values in spherical approximation, or point values in ellipsoidal approximation. The present study proposes a method
for computation of area mean gravity anomalies in ellipsoidal approximation (‘ellipsoidal area means’) by applying a simple
ellipsoidal correction to area means in spherical approximation. Ellipsoidal area means offer better consistency with GGM
quasigeoid heights. The method is numerically validated with ellipsoidal area mean gravity derived from very fine grids of
gravity point values in ellipsoidal approximation. Signal strengths of (i) the ellipsoidal effect (i.e., difference ellipsoidal
vs. spherical approximation), (ii) the area mean effect (i.e., difference area mean vs. point gravity) and (iii) the ellipsoidal
area mean effect (i.e., differences between ellipsoidal area means and point gravity in spherical approximation) are investigated
in test areas in New Zealand and the Himalaya mountains. The impact of both the area mean and the ellipsoidal effect on quasigeoid
heights is in the order of several centimetres. The proposed new gravity data type not only allows more accurate RCR-based
geoid computation, but may also be of some value for the GGM validation using terrestrial gravity anomalies that are available
as area mean values. 相似文献