首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 171 毫秒
1.
—We have developed a least-squares minimization approach to depth determination from magnetic data. By defining the anomaly value T(0) at the origin and the anomaly value T(N) at any other distance (N) on the profile, the problem of depth determination from magnetic data has been transformed into finding a solution to a nonlinear equation of the form f(z)=0. Formulas have been derived for a sphere, horizontal cylinder, dike, and for a geologic contact. Procedures are also formulated to estimate the effective magnetization intensity and the effective magnetization inclination. A scheme for analyzing the magnetic data has been formulated for determining the model parameters of the causative sources. The method is applied to synthetic data with and without random errors. Finally, the method is applied to two field examples from Canada and Arizona. In all cases examined, the estimated depths are found to be in goodagreement with actual values.  相似文献   

2.
This paper presents a neural network approach to determine 2D inverse modeling of a buried structure from gravity anomaly profile. The results of the applied neural network method are compared with the results of two other methods, least-squares minimization and the simple method. Sphere, horizontal cylinder and vertical cylinder and their gravity effects are considered as the synthetic models and the synthetic data, respectively. The synthetic data are also corrupted with noise to evaluate the capability of the methods. Then the Dehloran bitumen map in Iran is chosen as a real data application. Anomaly value of the cross-section, which is taken from the gravity anomaly map of Dehloran bitumen, is very close to those obtained from these methods.  相似文献   

3.
We have developed a least‐squares minimization approach to depth determination using numerical second horizontal derivative anomalies obtained from magnetic data with filters of successive window lengths (graticule spacings). The problem of depth determination from second‐derivative magnetic anomalies has been transformed into finding a solution to a non‐linear equation of the form, f(z) = 0. Formulae have been derived for a sphere, a horizontal cylinder, a dike and a geological contact. Procedures are also formulated to estimate the magnetic angle and the amplitude coefficient. We have also developed a simple method to define simultaneously the shape (shape factor) and the depth of a buried structure from magnetic data. The method is based on computing the variance of depths determined from all second‐derivative anomaly profiles using the above method. The variance is considered a criterion for determining the correct shape and depth of the buried structure. When the correct shape factor is used, the variance of depths is less than the variances computed using incorrect shape factors. The method is applied to synthetic data with and without random errors, complicated regionals, and interference from neighbouring magnetic rocks. Finally, the method is tested on a field example from India. In all the cases examined, the depth and the shape parameters are found to be in good agreement with the actual parameters.  相似文献   

4.
埋藏在6000m以下的地震反射异常体是礁体还是火成岩体?或者是海底滑塌体?在高质量重磁电资采集的基础上,通过对研究区域内的地层岩石物性细致研究,利用CEMP资料界面约束二维反演、重磁力异常剥离技术和LCT重磁力地震联合反演解释技术细致研究了异常体的物性参数,从密度、磁化率、电阻率、速度等多种特征推测判断了深层异常体地质属性,为钻井部署提供了间接依据。  相似文献   

5.
6.
The Hankel transform theorem can be applied to the inversion of gravity data for the buried sphere, the horizontal cylinder, and the vertical rod. This new approach leads to exact solutions of the transforms for the assumed bodies. A comparison with the classical procedure by Fourier transform reveals that for the sphere and the vertical rod, the Hankel transform is preferable.  相似文献   

7.
The inductive response of a conducting horizontal cylinder embedded in a uniform earth is studied using numerical results obtained for an analytical solution for the problem of a conducting cylinder buried in a homogeneous earth for the case of a uniform inducing field. A check of the validity of the numerical results is made by a comparison with analogue model measurements for a number of cases. Numerical results for a range of cylinder radii (a = 1–10 km), depths of burial (d= 0–4 km), conductivity contrasts (σ2= 10?2-10 Sm?1), and source frequencies (f= 10?1-10?4 Hz) of interest in the interpretation of magnetotelluric field measurements are presented. The results indicate that for a uniform inducing field the conductivity and depth of burial of a horizontal cylindrical inhomogeneity are best determined through a measurement of the amplitudes Hy, Hz and Ex and the phases φy and Ψx.  相似文献   

8.
A method is presented for determining bounds of the properties of axial symmetric bodies from a finite number of gravity and magnetic observations based on Parker's theory of ideal bodies. Bounds on the density contrast and the intensity of magnetization are calculated as a function of depth to the top of the anomalous source, restricting the range of smallest possible solutions to fit the data. The model studied is approximated by an array of vertical annuli cylinders, each of uniform density and magnetization. Linear programming algorithms based on the ideal body theory were used to calculate the distribution of these parameters within the body. Simultaneous inversion of gravity and magnetic data is performed assuming a constant ratio between the density contrast and the intensity of magnetization and that a common body is responsible for both observed fields. The parameter k(|J|/δp) provides information about the rock type of the structure. Interpretation of gravity and aeromagnetic data from Darnley Bay, NWT, Canada, indicated the presence of a shallow ultrabasic intrusion.  相似文献   

9.
The basic theoretical foundations of the spectral analysis of gravity and magnetic data based on the approximation approach have been considered and corresponding expressions implementing various statements of the F-approximation method have been derived [Strakhov and Kerimov, 1999; 2000; 2001]. The author considered the algorithms and computer implementation of the F-approximation of anomalous gravity and magnetic fields presented in the results of testing with the use of a number of model examples and data on the accuracy of the model field approximation earlier on in [Kerimov, 2003]. In this paper, I consider methodological features of constructing the F-approximation from the detailed gravity and magnetic surveys, anomalous gravity and magnetic fields reconstructed with the use of the F -approximation, and examples of the field transforms obtained with the use of this approach.  相似文献   

10.
In this paper, the coefficients of Parkinson vectora, b and the in-phase partA r, Br and quadrature-phase partA i, Bi of transter functions are calculated by using two kinds of data processing methods, respectively. The results of both methods are close to each other. All these six parameters above are low in value, showing that the electric structure in Heze area is rather uniform. It seems that the anomalous changes appeared in both results, and the values of parameters decrease before and after Heze earthquake.  相似文献   

11.
An interpretative method based on a nonlinearly mathematical optimization concept has been developed in this paper, in order to interpret self-potential anomalies (SP) due to horizontal cylinder, vertical cylinder, sphere and sheet-like structures. This interpretative method comprises three main steps. The first step is to formulate mathematically a nonlinearly constrained minimization problem (NCMP) to describe the geophysical problem related to the studied structure. The second one is to suggest an interior penalty function in order to convert the nonlinearly constrained minimization problem (NCMP) into a nonlinearly unconstrained minimization one (NUMP). The third step is to solve the converted nonlinearly unconstrained minimization problem (NUMP) by the well-known Hooke and Jeeves direct search algorithm in order to estimate the geophysical parameters of the studied structure, i.e., depth, polarization angle, electric dipole moment (magnitude of polarization) and geometric shape factor. The Hooke and Jeeves direct search algorithm is purposely chosen for being robust and its application to SP data allows a rapid convergence towards the optimal estimate of parameters. This interpretative method was first tested on theoretical synthetic models with different random noise, where a very close agreement was obtained between assumed and evaluated parameters.The validity of the proposed interpretative method is also tested on practical field examples taken from Turkey, India and Germany, where available SP data existed and was previously analyzed by different interpretative methods. The agreement between the results obtained by the developed method and those obtained by other published methods is good.Acknowledgment Authors would like to thank Dr. I. Othman Director General of the Atomic Energy Commission of Syria for his interest and continuous encouragement to achieve this work. Special thanks to the reviewers for their constructive suggestions aimed at enhancing the quality of this paper.  相似文献   

12.
The paper presents data of electrical conductivity monitoring for the upper crust obtained with the use of a controllable electromagnetic source by the dipole, vertical electromagnetic, and transient electromagnetic sounding methods, as well as variations in the geomagnetic field and the Lake Baikal water level. With the help of spectral, spectral time, and multivariate analyses, main spatial and temporal patterns of variations in geophysical parameters are determined. The present study revealed anomalous effects in the behavior of the canonical coherence (parameter K) and aggregated signal (parameter A) and other anomalous variations in geophysical parameters related to earthquakes in the Baikal rift zone.  相似文献   

13.
During the last 15 years, more attention has been paid to derive analytic formulae for the gravitational potential and field of polyhedral mass bodies with complicated polynomial density contrasts, because such formulae can be more suitable to approximate the true mass density variations of the earth (e.g., sedimentary basins and bedrock topography) than methods that use finer volume discretization and constant density contrasts. In this study, we derive analytic formulae for gravity anomalies of arbitrary polyhedral bodies with complicated polynomial density contrasts in 3D space. The anomalous mass density is allowed to vary in both horizontal and vertical directions in a polynomial form of \(\lambda =ax^m+by^n+cz^t\), where mnt are nonnegative integers and abc are coefficients of mass density. First, the singular volume integrals of the gravity anomalies are transformed to regular or weakly singular surface integrals over each polygon of the polyhedral body. Then, in terms of the derived singularity-free analytic formulae of these surface integrals, singularity-free analytic formulae for gravity anomalies of arbitrary polyhedral bodies with horizontal and vertical polynomial density contrasts are obtained. For an arbitrary polyhedron, we successfully derived analytic formulae of the gravity potential and the gravity field in the case of \(m\le 1\), \(n\le 1\), \(t\le 1\), and an analytic formula of the gravity potential in the case of \(m=n=t=2\). For a rectangular prism, we derive an analytic formula of the gravity potential for \(m\le 3\), \(n\le 3\) and \(t\le 3\) and closed forms of the gravity field are presented for \(m\le 1\), \(n\le 1\) and \(t\le 4\). Besides generalizing previously published closed-form solutions for cases of constant and linear mass density contrasts to higher polynomial order, to our best knowledge, this is the first time that closed-form solutions are presented for the gravitational potential of a general polyhedral body with quadratic density contrast in all spatial directions and for the vertical gravitational field of a prismatic body with quartic density contrast along the vertical direction. To verify our new analytic formulae, a prismatic model with depth-dependent polynomial density contrast and a polyhedral body in the form of a triangular prism with constant contrast are tested. Excellent agreements between results of published analytic formulae and our results are achieved. Our new analytic formulae are useful tools to compute gravity anomalies of complicated mass density contrasts in the earth, when the observation sites are close to the surface or within mass bodies.  相似文献   

14.
In the present study, the deterministic chaotic behaviour of interplanetary magnetic field (IMF) under various geomagnetic conditions of low and high solar active periods was analyzed, using the time series of IMF |B| and Bz, by employing chaotic quantifiers like, Lyapunov exponent, Tsallis entropy, correlation dimension, and non-linear prediction error. We have investigated whether the chaotic behaviour of interplanetary magnetic field would modify, when it produces major geomagnetic storms, and how it depends on the phase of solar activity. The yearly average values of Lyapunov exponent for the time series of IMF |B| and Bz, show solar flux dependence, whereas those values of entropy, correlation dimension and non-linear prediction error had no significant solar flux dependence. The yearly average values of entropy for quiet periods are higher compared to those values for major storm periods belonging to low/high solar active conditions, for both the time series |B| and Bz.  相似文献   

15.
The gravity gradient tensor (GGT) is deduced from products of second-order derivatives of the gravitational potential. A new method based on the invariants of the GGT has been proposed in this research to interpret gravity data due to sphere, infinite horizontal cylinder and semi-infinite vertical cylinder. The method estimates the depth of these simple causative sources from the multiplication of the maximum of the gravity vertical component by the maximum value of the invariants I 1 to I 2 ratio. To show the reliability and correctness of the estimated depths on 3-D models, the method has been tested using theoretical data with and without random noise. In addition, I have applied the method to a field-data example in Texas, USA and the depth obtained by the present method is compared with those published in the literature.  相似文献   

16.
The correlation discovered by Hide and Malin between the variable parts of the Earth's gravitational field and magnetic field (suitably displaced in longitude) was tentatively and qualitatively explained by them in terms of the influence on both fields of irregularities (or “surface bumps”) at the core-mantle interface. In this paper, a quantitative analysis of this phenomenon is developed, through study of an idealised problem in which conducting fluid occupying the region z < η(x) flows over the surface z = η(x) in the presence of a magnetic field (B0,0,0), the whole system rotating with angular velocity (0,0,Ω). It is assumed that |η′(x)| « 1 so that perturbation methods are applicable. Determination of the magnetic potential in the “mantle” region z < η(x) requires solution of the full hydromagnetic problem in the fluid. It is shown that three wave modes are excited, two of which (for values of the parameters of the problem of geophysical interest) have a boundary layer character. Phase interactions between these modes lead to a shift and a distortion of the magnetic pattern relative to the gravitational pattern. The correlation between the gravitational potential and the magnetic potential (shifted by a distance x0) is determined on the plane z = d (d a? |η|) as a function of x0/d and the curves obtained are qualitatively similar to that based on the observed data; the maximum correlation obtained varies between 0.67 and 1, depending on values of the parameters of the problem, and is about 0.72 for reasonable estimates of these parameters in the geophysical context.  相似文献   

17.
Gravity and magnetic anomalies may have the same source but it is always difficult to achieve correlated solutions if interpretations are carried out separately. Therefore it is useful to invert both anomalous gravity and magnetic profiles at the same time, so that the solution mav be more constrained. Existing inversion techniques do not lend themselves equally to this kind of extension, since the parameters—such as density and susceptibility contrasts—should not be related. The algorithms more easily adaptable to simultaneous inversion are those which look for the shape and the position of anomalous bodies with constant unknown density or susceptibility contrasts. In this study, we use the generalized inverse method in the 2 1/2 D case. Parameters are the coordinates of the vertices of the polygonal cross-sections of the anomalous bodies and both density and susceptibility contrasts of each body. The two types of profile to be processed must be superposable but may have different lengths, spacings, and coordinate origins. Both synthetic examples and field data from geologically known areas have been processed, and it appears that the simultaneous inversion technique may provide an important additional degree of control in the interactive interpretation process.  相似文献   

18.
—?The aeromagnetic data of Macedonia and Thrace were used to produce Curie point estimates. The data were high pass filtered to remove components arising from topography and magnetic core fields which were not adequately modeled by a DGRF. The depth to the centroid, z 0, of the deepest distribution of the magnetic dipoles was obtained by computing a least-squares fit to the lowest-frequency segment of the azimuthally averaged log power spectrum. The average depth to the top of the deepest crustal block was computed as the depth to the top, z t , of the second lowest-frequency segment of the spectrum. The depth to the bottom of the deepest magnetic dipoles, the inferred Curie point depth, was then calculated from z b =2z 0???z t . The Curie depth estimates for Macedonia and Thrace range between 11.2 and 17.3?km. These results are consistent with the depths inferred by extrapolating known geothermal gradient and heat-flow values.  相似文献   

19.
—The aeromagnetic data of the island of Crete were inverted to produce Curie point estimates. The data were high-pass filtered to remove components arising from topography and magnetic core fields which were not adequately modeled by IGRF. The depth to the centroid, z 0?, of the deepest distribution of the magnetic dipoles, was obtained by computing a least-squares fit to the lowest-fre quency segment of the azimuthally averaged log power spectrum. The depth to the top of the deepest crustal block was computed as the depth, z t ?, to the centroid of the second deepest distribution, using the second lowest-frequency segment of the spectrum. The depth to the bottom of the deepest dipoles, the inferred Curie point depth, is then z b = 2z 0?z t ?. The Curie depth estimates range between 24 and 28 km. This is in accordance with the depths inferred by extrapolating heat-flow values measured in boreholes.  相似文献   

20.
有限长圆柱体磁异常场全空间正演方法   总被引:2,自引:0,他引:2       下载免费PDF全文
在经典位场理论中,许多简单形体位场异常难以通过积分得到全空间的解析式.圆柱体是一类很重要的理论模型体,常用于模拟圆柱状地质体或非地质体(如管线),但目前还不能用解析公式正演有限长圆柱体在三维空间里的磁异常,而多是采用近似简化为有限长磁偶极子或线模型代替.对于有限长圆柱体,特别是半径相对于上顶埋深较大时,这种近似的误差不可忽略.本文利用共轭复数变量替换法,推导出有限长圆柱体在全空间的引力位一阶、二阶导数,利用Poisson关系得到磁异常正演公式,进而利用有限长圆柱体磁异常正演公式求解管状体的磁异常,得到不同磁化方向、不同大小的管线产生的磁场的特征,并将其推广到截面为椭圆的情况.最后通过模拟计算定量给出了将圆柱体近似为线模型的条件.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号