首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 109 毫秒
1.
We have developed three different least-squares approaches to determine successively: the depth, magnetic angle, and amplitude coefficient of a buried sphere from a total magnetic anomaly. By defining the anomaly value at the origin and the nearest zero-anomaly distance from the origin on the profile, the problem of depth determination is transformed into the problem of finding a solution of a nonlinear equation of the form f(z)=0. Knowing the depth and applying the least-squares method, the magnetic angle and amplitude coefficient are determined using two simple linear equations. In this way, the depth, magnetic angle, and amplitude coefficient are determined individually from all observed total magnetic data. The method is applied to synthetic examples with and without random errors and tested on a field example from Senegal, West Africa. In all cases, the depth solutions are in good agreement with the actual ones.  相似文献   

2.
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.  相似文献   

3.
A simple and fast determination of the limiting depth to the sources may represent a significant help to the data interpretation. To this end we explore the possibility of determining those source parameters shared by all the classes of models fitting the data. One approach is to determine the maximum depth-to-source compatible with the measured data, by using for example the well-known Bott–Smith rules. These rules involve only the knowledge of the field and its horizontal gradient maxima, and are independent from the density contrast.Thanks to the direct relationship between structural index and depth to sources we work out a simple and fast strategy to obtain the maximum depth by using the semi-automated methods, such as Euler deconvolution or depth-from-extreme-points method (DEXP).The proposed method consists in estimating the maximum depth as the one obtained for the highest allowable value of the structural index (Nmax). Nmax may be easily determined, since it depends only on the dimensionality of the problem (2D/3D) and on the nature of the analyzed field (e.g., gravity field or magnetic field). We tested our approach on synthetic models against the results obtained by the classical Bott–Smith formulas and the results are in fact very similar, confirming the validity of this method. However, while Bott–Smith formulas are restricted to the gravity field only, our method is applicable also to the magnetic field and to any derivative of the gravity and magnetic field. Our method yields a useful criterion to assess the source model based on the (∂f/∂x)max/fmax ratio.The usefulness of the method in real cases is demonstrated for a salt wall in the Mississippi basin, where the estimation of the maximum depth agrees with the seismic information.  相似文献   

4.
We have developed a new numerical method to determine the shape (shape factor), depth, polarization angle, and electric dipole moment of a buried structure from residual self-potential (SP) anomalies. The method is based on defining the anomaly value at the origin and four characteristic points and their corresponding distances on the anomaly profile. The problem of shape determination from residual SP anomaly has been transformed into the problem of finding a solution to a nonlinear equation of the form q = f (q). Knowing the shape, the depth, polarization angle and the electric dipole moment are determined individually using three linear equations. Formulas have been derived for spheres and cylinders. By using all possible combinations of the four characteristic points and their corresponding distances, a procedure is developed for automated determination of the best-fit-model parameters of the buried structure from SP anomalies. The method was applied to synthetic data with 5% random errors and tested on a field example from Colorado. In both cases, the model parameters obtained by the present method, particularly the shape and depth of the buried structures are found in good agreement with the actual ones. The present method has the capability of avoiding highly noisy data points and enforcing the incorporation of points of the least random errors to enhance the interpretation results.  相似文献   

5.
In this paper, we describe a non‐linear constrained inversion technique for 2D interpretation of high resolution magnetic field data along flight lines using a simple dike model. We first estimate the strike direction of a quasi 2D structure based on the eigenvector corresponding to the minimum eigenvalue of the pseudogravity gradient tensor derived from gridded, low‐pass filtered magnetic field anomalies, assuming that the magnetization direction is known. Then the measured magnetic field can be transformed into the strike coordinate system and all magnetic dike parameters – horizontal position, depth to the top, dip angle, width and susceptibility contrast – can be estimated by non‐linear least squares inversion of the high resolution magnetic field data along the flight lines. We use the Levenberg‐Marquardt algorithm together with the trust‐region‐reflective method enabling users to define inequality constraints on model parameters such that the estimated parameters are always in a trust region. Assuming that the maximum of the calculated gzz (vertical gradient of the pseudogravity field) is approximately located above the causative body, data points enclosed by a window, along the profile, centred at the maximum of gzz are used in the inversion scheme for estimating the dike parameters. The size of the window is increased until it exceeds a predefined limit. Then the solution corresponding to the minimum data fit error is chosen as the most reliable one. Using synthetic data we study the effect of random noise and interfering sources on the estimated models and we apply our method to a new aeromagnetic data set from the Särna area, west central Sweden including constraints from laboratory measurements on rock samples from the area.  相似文献   

6.
Possible configurations of the magnetic field in the outer magnetosphere during geomagnetic polarity reversals are investigated by considering the idealized problem of a magnetic multipole of order m and degree n located at the centre of a spherical cavity surrounded by a boundless perfect diamagnetic medium. In this illustrative idealization, the fixed spherical (magnetopause) boundary layer behaves as a perfectly conducting surface that shields the external diamagnetic medium from the compressed multipole magnetic field, which is therefore confined within the spherical cavity. For a general magnetic multipole of degree n, the non-radial components of magnetic induction just inside the magnetopause are increased by the factor 1 + [(n + 1)/n] relative to their corresponding values in the absence of the perfectly conducting spherical magnetopause. An exact equation is derived for the magnetic field lines of an individual zonal (m = 0), or axisymmetric, magnetic multipole of arbitrary degree n located at the centre of the magnetospheric cavity. For such a zonal magnetic multipole, there are always two neutral points and n – 1 neutral rings on the spherical magnetopause surface. The two neutral points are located at the poles of the spherical magnetopause. If n is even, one of the neutral rings is coincident with the equator; otherwise, the neutral rings are located symmetrically with respect to the equator. The actual existence of idealized higher-degree (n > 1) axisymmetric magnetospheres would necessarily imply multiple (n + 1) magnetospheric cusps and multiple (n) ring currents. Exact equations are also derived for the magnetic field lines of an individual non-axisymmetric magnetic multipole, confined by a perfectly conducting spherical magnetopause, in two special cases; namely, a symmetric sectorial multipole (m = n) and an antisymmetric sectorial multipole (m = n – 1). For both these non-axisymmetric magnetic multipoles, there exists on the spherical magnetopause surface a set of neutral points linked by a network of magnetic field lines. Novel magnetospheric processes are likely to arise from the existence of magnetic neutral lines that extend from the magnetopause to the surface of the Earth. Finally, magnetic field lines that are confined to, or perpendicular to, either special meridional planes or the equatorial plane, when the multipole is in free space, continue to be confined to, or perpendicular to, these same planes when the perfectly conducting magnetopause is present.Also Honorary Research Associate, Rutherford Appleton Laboratory, Chilton, Didcot, Oxon OX11 0QX, UK and Visiting Reader in Physics. University of Sussex, Falmer, Brighton BN1 9QH, UK  相似文献   

7.
Summary A method called the zero-line method for interpreting the total field magnetic anomalies of spherical bodies is described. The contour of zero-anomaly is controlled by a linear equation of the typeAx 2 +By 2 +Cxy +Dx +Ey +F = 0, where the coefficientsA toF are related to the depth of the sphere and the parameters relating to its position and magnetisation. From the coordinates (x, y) of the various points lying on the zero-contour, five normal equations for the above equation are set and solved for the various coefficients. The various parameters relating to the position and magnetisation of the sphere are then determined from these coefficients. It is expected that similar methods can be developed for interpreting magnetic anomalies of other three dimensional bodies also.  相似文献   

8.
We have developed a least-squares method to determine simultaneously the depth and the width of a buried thick dipping dike from residualized magnetic data using filters of successive window lengths. The method involves using a relationship between the depth and the half-width of the source and a combination of windowed observations. The relationship represents a family of curves (window curves). For a fixed window length, the depth is determined for each half-width value by solving one nonlinear equation of the form f (z) = 0 using the least-squares method. The computed depths are plotted against the width values representing a continuous curve. The solution for the depth and the width of the buried dike is read at the common intersection of the window curves. The method involves using a dike model convolved with the same moving average filter as applied to the observed data. As a result, this method can be applied to residuals as well as to measured magnetic data. Procedures are also formulated to estimate the amplitude coefficient and the index parameter. The method is applied to theoretical data with and without random errors. The validity of the method is tested on airborne magnetic data from Canada and on a vertical component magnetic anomaly from Turkey. In all cases examined, the model parameters obtained are in good agreement with the actual ones and with those given in the published literature.  相似文献   

9.
Nonparametric inverse methods provide a general framework for solving potential‐field problems. The use of weighted norms leads to a general regularization problem of Tikhonov form. We present an alternative procedure to estimate the source susceptibility distribution from potential field measurements exploiting inversion methods by means of a flexible depth‐weighting function in the Tikhonov formulation. Our approach improves the formulation proposed by Li and Oldenburg (1996, 1998) , differing significantly in the definition of the depth‐weighting function. In our formalism the depth weighting function is associated not to the field decay of a single block (which can be representative of just a part of the source) but to the field decay of the whole source, thus implying that the data inversion is independent on the cell shape. So, in our procedure, the depth‐weighting function is not given with a fixed exponent but with the structural index N of the source as the exponent. Differently than previous methods, our choice gives a substantial objectivity to the form of the depth‐weighting function and to the consequent solutions. The allowed values for the exponent of the depth‐weighting function depend on the range of N for sources: 0 ≤N≤ 3 (magnetic case). The analysis regarding the cases of simple sources such as dipoles, dipole lines, dykes or contacts, validate our hypothesis. The study of a complex synthetic case also proves that the depth‐weighting decay cannot be necessarily assumed as equal to 3. Moreover it should not be kept constant for multi‐source models but should instead depend on the structural indices of the different sources. In this way we are able to successfully invert the magnetic data of the Vulture area, Southern Italy. An original aspect of the proposed inversion scheme is that it brings an explicit link between two widely used types of interpretation methods, namely those assuming homogeneous fields, such as Euler deconvolution or depth from extreme points transformation and the inversion under the Tikhonov‐form including a depth‐weighting function. The availability of further constraints, from drillings or known geology, will definitely improve the quality of the solution.  相似文献   

10.
The magnetic anomaly caused by a buried dike is separated into its even and odd components, which have a simple symmetry with respect to the origin. These values are integrated up to the half-maximum abscissa for the even component, and the maximum abscissa for the odd component. The integration nomograms are generated using various values to the half-width and depth in the theoretical anomaly equations. These nomograms are used to determine the half-width and depth to the top of the dike for the field anomaly. The method also includes the determination of the index parameter (Q) and the amplitude coefficient (P). An example using theoretical data shows the effectiveness of the present method.  相似文献   

11.
The ‘depth from extreme points’ method is an important tool to estimate the depth of sources of gravity and magnetic data. In order to interpret gravity gradient tensor data conveniently, formulas for the tensor data form regarding depth from the extreme points method were calculated in this paper. Then, all of the gradient tensor components were directly used to interpret the causative source. Beyond the gzz component, also the gxx and gyy components can be used to obtain depth information. In addition, the total horizontal derivative of the depth from extreme points of the gradient tensor can be used to describe the edge information of geologic sources. In this paper, we investigated the consistency of the homogeneity degree calculated by using the different components, which leads to the calculated depth being confirmed. Therefore, a more integrated interpretation can be obtained by using the gradient tensor components. Different synthetic models were used with and without noise to test the new approach, showing stability, accuracy and speed. The proposed method proved to be a useful tool for gradient tensor data interpretation. Finally, the proposed method was applied to full tensor gradient data acquired over the Vinton Salt Dome, Louisiana, USA, and the results are in agreement with those obtained in previous research studies.  相似文献   

12.
Ground total magnetic field data of Albania were used to produce estimates of the Curie point isotherm. The strategy followed was to estimate the depth to the bottom of the deepest magnetic sources. Firstly, the average depth to the top of the deepest crustal block, zt, was computed by linear fitting to the second lowest-frequency segment of the azimuthally averaged power spectrum of the total magnetic field data. Then, the depth to the centroid of the deepest crustal block, z0, was computed by linear fitting to the lowest-frequency segment of the azimuthally averaged power spectrum of a distribution of magnetic dipoles. Finally, the depth to the bottom, the inferred Curie point depth, zb, was calculated from zb=2z0−zt. Curie depth estimates for Albania vary from about 17 to 25.5 km (below sea level). These results are consistent with the depths inferred by extrapolating geothermal gradient and heat-flow values, suggesting that the Curie point depth analysis is useful to estimate the regional thermal structure. It also suggests that the approach was valid and that ground total magnetic field data can be used for this purpose.  相似文献   

13.
Interpretation of magnetic anomalies of dikes using correlation factors   总被引：1，自引：0，他引：1  
The magnetic anomaly due to a buried dike consists of the sum of two easily separated elementary functions. These functions, which have simple symmetry, are called even and odd functions. The correlation factors (r 0,1 for the even andr 0,2 for the odd function) between least-squares residual anomalies from even and odd functions are computed. Correlation values are used to determine the depth to the top and the half-width of the dike. The method also includes the determination of the index parameter and the amplitude coefficient. The validity of the method is tested against a theoretical and a field example where the parameters of the latter were determined by other investigators in comparing the results.  相似文献   

14.
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.  相似文献   

15.
We have developed a least-squares minimization approach to depth determination of a buried ore deposit from numerical horizontal gradients obtained from self-potential (SP) data using filters of successive window lengths (graticule spacings). The problem of depth determination from SP gradients has been transformed into the problem of finding a solution to a nonlinear equation of the form f(z)=0. Formulas have been derived for vertical and horizontal cylinders and spheres. Procedures are also formulated to estimate the electrical dipole moment and the polarization angle. The method is applied to synthetic data with and without random noise. Finally, the validity of the method is tested on two field examples. In both cases, the depth obtained is found to be in a very good agreement with that obtained from drilling information.  相似文献   

16.
Abstract

Two open curves with fixed endpoints on a boundary surface can be topologically linked. However, the Gauss linkage integral applies only to closed curves and cannot measure their linkage. Here we employ the concept of relative helicity in order to define a linkage for open curves. For a magnetic field consisting of closed field lines, the magnetic helicity integral can be expressed as the sum of Gauss linkage integrals over pairs of lines. Relative helicity extends the helicity integral to volumes where field lines may cross the boundary surface. By analogy, linkages can be defined for open lines by requiring that their sum equal the relative helicity.

With this definition, the linkage of two lines which extend between two parallel planes simply equals the number of turns the lines take about each other. We obtain this result by first defining a gauge-invariant, one-dimensional helicity density, i.e. the relative helicity of an infinitesimally thin plane slab. This quantity has a physical interpretation in terms of the rate at which field lines lines wind about each other in the direction normal to the plane. A different method is employed for lines with both endpoints on one plane; this method expresses linkages in terms of a certain Gauss linkage integral plus a correction term. In general, the linkage number of two curves can be put in the form L=r + n, |r|≦1J2, where r depends only on the positions of the endpoints, and n is an integer which reflects the order of braiding of the curves.

Given fixed endpoints, the linkage numbers of a magnetic field are ideal magneto-hydrodynamic invariants. These numbers may be useful in the analysis of magnetic structures not bounded by magnetic surfaces, for example solar coronal fields rooted in the photosphere. Unfortunately, the set of linkage numbers for a field does not uniquely determine the field line topology. We briefly discuss the problem of providing a complete and economical classification of field topologies, using concepts from the theory of braid equivalence classes.  相似文献   

17.
Abstract

The structure and collapse of linear three-dimensional magnetic neutral points is studied by varying the four parameters (p, q,j|,j⊥ ) that define, in general, the linear field of a neutral point. The effect of these parameters on both the skeleton structure (i.e. the fan and spine) and the actual field line structure of the null is considered. It is found that one current component (j⊥ ) causes the skeleton structure of the null to fold up from its potential state, whereas the other current component (j |;) causes the field lines to bend. The two other parameters (p,q) determine the potential structure of the null and cause the null to transform from a three-dimensional null to a two-dimensional null and from a positive (type B) null to a negative (type A) null.

To investigate the collapse of three-dimensional nulls, solutions to the linear, low-β ideal magnetohydrodynamic equations are found. It is found that three-dimensional null points can collapse if the field line foot-points are free and energy can propagate into the system.  相似文献   

18.
We use the continuous wavelet transform based on complex Morlet wavelets, which has been developed to estimate the source distribution of potential fields. For magnetic anomalies of adjacent sources, they always superimpose upon each other in space and wavenumber, making the identification of magnetic sources problematic. Therefore, a scale normalization factor, a?n, is introduced on the wavelet coefficients to improve resolution in the scalogram. By theoretical modelling, we set up an approximate linear relationship between the pseudo‐wavenumber and source depth. The influences of background field, random noise and magnetization inclination on the continuous wavelet transform of magnetic anomalies are also discussed and compared with the short‐time Fourier transform results. Synthetic examples indicate that the regional trend has little effect on our method, while the influence of random noise is mainly imposed on shallower sources with higher wavenumbers. The source horizontal position will be affected by the change of magnetization direction, whereas the source depth remains unchanged. After discussing the performance of our method by showing the results of various synthetic tests, we use this method on the aeromagnetic data of the Huanghua depression in central China to define the distribution of volcanic rocks. The spectrum slices in different scales are used to determine horizontal positions of volcanic rocks and their source depths are estimated from the modulus maxima of complex coefficients, which is in good accordance with drilling results.  相似文献   

19.
—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.  相似文献   

20.
The estimation of the depth to the top and bottom of a magnetic source from magnetic data defines a nonlinear inverse problem, while the evaluation of the distribution of magnetization determines a linear inverse problem. In this paper, these interpretation problems are resolved in the continuous case of 21/2D magnetized bodies with lateral magnetization variations. A formulation of the magnetic problem accounting for different directions of remanent and total magnetization vectors and including a more general definition of apparent susceptibility is presented. Differences between 2D and 21/2D formulations are stressed, as regards the anomaly amplitude, shape and zero-level.In order to utilize well-known continuous linear inverse methods, Fréchet derivatives of the data functionals with respect to the depth of the source top and bottom, are analytically described. Thus, using the spectral expansion inverse method (Parker, 1977) and linearizing the problem at several steps of an iterative process, the source depth is obtained within a few iterations, although the starting model is distant from the final solution. The interpretation of an anomaly in the Italian region shows the usefulness of the method.  相似文献   

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

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