首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 625 毫秒
1.
The current local wavenumber methods for the interpretation of magnetic anomalies compute the locations of geological bodies by solving complex matrices. Presently, such methods require to know the structural index, which is a parameter that represents the source type. The structural index is hard to know in real data; consequently, the precision of current methods is low. We present the fast local wavenumber (FLW) method, and define the squared sum of the horizontal and vertical local wavenumbers as the cumulative local wavenumber. The FLW method is the linear combination of the umulative local wavenumberand other wavenumbers, and is used to compute the locations and structural index of the source without a priori information and matrix solution. We apply the FLW method to synthetic magnetic anomalies, and the results suggest that the FLW method is insensitive to background and oblique magnetization. Next, we apply the FLW method to real magnetic data to obtain the location and structural index of the source.  相似文献   

2.
This paper presents a new inversion method for the interpretation of 2D magnetic anomaly data, which uses the combination of the analytic signal and its total gradient to estimate the depth and the nature (structural index) of an isolated magnetic source. However, our proposed method is sensitive to noise. In order to lower the effect of noise, we apply upward continuation technique to smooth the anomaly. Tests on synthetic noise-free and noise corrupted magnetic data show that the new method can successfully estimate the depth and the nature of the causative source. The practical application of the technique is applied to measured magnetic anomaly data from Jurh area, northeast China, and the inversion results are in agreement with the inversion results from Euler deconvolution of the analytic signal.  相似文献   

3.
We presented using the correlation coefficient of the analytic signal of real data and the analytic signal of synthetic data generated by the assumed source to estimate the structural index and the depth of the source. First, we assumed that the causative sources are located at different locations in the underground and the structural index of the assumed source is changed from 0 to 3, and then we separately compute the correlation coefficients of the analytic signal of the measured data and the analytic signal of the anomaly generated by each assumed source, the correlation coefficient can get the maximum value when the location and structural index of the assumed source are consistent with the real source. We tested the correlation coefficient method on synthetic noise-free and noise-corrupted magnetic anomalies, and the inversion results indicate that the new method can successfully finish the inversion of magnetic data. We also applied it to measured magnetic data, and we obtain the structural index and the location of the source.  相似文献   

4.
Recent improvements in the local wavenumber approach have made it possible to estimate both the depth and model type of buried bodies from magnetic data. However, these improvements require calculation of third‐order derivatives of the magnetic field, which greatly enhances noise. As a result, the improvements are restricted to data of high quality. We present an alternative method to estimate both the depth and model type using the first‐order local wavenumber approach without the need for third‐order derivatives of the field. Our method is based on normalization of the first‐order local wavenumber anomalies and provides a generalized equation to estimate the depth of some 2D magnetic sources regardless of the source structure. Information about the nature of the sources is obtained after the source location has been estimated. The method was tested using synthetic magnetic anomaly data with random noise and using three field examples.  相似文献   

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

6.
Magnetic anomalies are often disturbed by the magnetization direction, so we can’t directly use the original magnetic anomaly to estimate the exact location and geometry of the source. The 2D analytic signal is insensitive to magnetization direction. In this paper, we present an automatic method based on the analytic signal horizontal and vertical derivatives to interpret the magnetic anomaly. We derive a linear equation using the analytic signal properties and we obtain the 2D magnetic body location parameters without giving a priori information. Then we compute the source structural index (expressing the geometry) by the estimated location parameters. The proposed method is demonstrated on synthetic magnetic anomalies with noise. For different models, the proposed technique can both successfully estimate the location parameters and the structural index of the sources and is insensitive to noise. Lastly, we apply it to real magnetic anomalies from China and obtain the distribution of unexploited iron ore. The inversion results are consistent with the parameters of known ore bodies.  相似文献   

7.
本文提出归一化总水平导数法,通过对总水平导数进行空间归一化计算实现了异常体水平位置和深度的估计,此外还推导出基于归一化总水平导数的欧拉反褶积法来估算地下地质体的空间位置,两种方法反演结果的相互验证可有效地提高反演结果的可信度.理论模型试验证明空间归一化总水平导数法和归一化总水平导数欧拉反褶积法均能有效地完成异常体的水平位置和深度的估计,所获得的位置参数与理论值相一致.在利用归一化总水平导数法进行磁异常解释时,对数据进行化磁极计算可得到更加准确的结果.将其应用于实际航磁数据的解释,获得了岩脉的大致分布特征.  相似文献   

8.
In this paper, we present a case study on the use of the normalized source strength (NSS) for interpretation of magnetic and gravity gradient tensors data. This application arises in exploration of nickel, copper and platinum group element (Ni‐Cu‐PGE) deposits in the McFaulds Lake area, Northern Ontario, Canada. In this study, we have used the normalized source strength function derived from recent high resolution aeromagnetic and gravity gradiometry data for locating geological bodies. In our algorithm, we use maxima of the normalized source strength for estimating the horizontal location of the causative body. Then we estimate depth to the source and structural index at that point using the ratio between the normalized source strength and its vertical derivative calculated at two levels; the measurement level and a height h above the measurement level. To discriminate more reliable solutions from spurious ones, we reject solutions with unreasonable estimated structural indices. This method uses an upward continuation filter which reduces the effect of high frequency noise. In the magnetic case, the advantage is that, in general, the normalized magnetic source strength is relatively insensitive to magnetization direction, thus it provides more reliable information than standard techniques when geologic bodies carry remanent magnetization. For dipping gravity sources, the calculated normalized source strength yields a reliable estimate of the source location by peaking right above the top surface. Application of the method on aeromagnetic and gravity gradient tensor data sets from McFaulds Lake area indicates that most of the gravity and magnetic sources are located just beneath a 20 m thick (on average) overburden and delineated magnetic and gravity sources which can be probably approximated by geological contacts and thin dikes, come up to the overburden.  相似文献   

9.
Balanced edge detection filters can recognize the edges of the shallow and deep bodies simultaneously, and are commonly used in the edge detection of potential field data. In this paper, we present using the balanced edge detection filters to estimate source locations, and derive two linear equations based on the balanced edge detection filters that can estimate the locations of the source without any priori information about the nature (structural index) of the source. The proposed methods are demonstrated on synthetic gravity anomalies, and the inversion results show that the proposed methods can successfully estimate location parameters of the sources. I also apply the proposed methods to real magnetic data, and the inversion results estimated by the proposed methods are consistent with the results estimated by the other similar method.  相似文献   

10.
We presented a new method for interpreting 2D magnetic data, called direct analytic signal (DAS) method, which directly used the analytic signal of magnetic anomaly to compute the depth and the structural index of the source. The DAS method needs only the computation of the first order derivatives of magnetic anomaly, so that the inversion results are more stable than the results obtained by the other existing analytic signal methods. The DAS method is tested on synthetic magnetic data with and without noise, and the DAS method can successfully obtain the depth and the structural index of the source. We also applied the DAS method to interpret a real magnetic data over a shallow geological source whose source parameters are known from closely drilling information, and the inversion results are in accord with the true values.  相似文献   

11.
针对天然大地电磁场信号在人文活动密集地区易受噪声干扰的问题,本文提出利用两个同步测点天然电磁场时间序列之间的单位脉冲响应,合成本地点受干扰时段的数据,从而去除大地电磁噪声.首先,选择高信噪比时段的数据,采用最小二乘法,估算本地点与参考点之间的单位脉冲响应,再根据卷积定律,结合参考磁场合成本地点的磁场和电场.最后用合成数据替换含噪声时段数据,实现时间域去噪.实测高信噪比数据和含噪数据的处理结果表明,该方法可以高精度合成本地点磁场与电场信号,有效去除本地点电场和磁场噪声,包括相关噪声,提高大地电磁数据质量.  相似文献   

12.
We have developed a method for imaging magnetic data collected for mineral exploration to yield the following structural information: depth, model type (structural index) and susceptibility. The active nature of mineral exploration data requires we derive the structural information from a robust quantity: we propose that the first‐ or second‐order analytic‐signal amplitude is suitably stable. The procedure is to normalize the analytic‐signal amplitude by the peak value and then use non‐linear inversion to estimate the depth and the structural index for each anomaly. In our field example, different results are obtained depending on whether we inverted for the first‐ or second‐order analytic‐signal amplitude. This is probably because the two‐dimensional contact, thin sheet or horizontal cylinder models we have assumed are not appropriate. In cases such as these, when our model assumptions are not correct, the results should not be interpreted quantitatively, but they might be useful for giving a qualitative indication of how the structure might vary. With a priori information, it is possible to assume a model type (i.e. set the structural index) and generate estimates of the depth and susceptibility. These data can then be gridded and imaged. If a contact is assumed, the susceptibility contrast is estimated; for the dike model, the susceptibility‐thickness is estimated; for the horizontal cylinder, the susceptibility‐area is estimated. To emphasize that the results are dependent on our assumed model, we advocate prefixing any derived quantity by the term ‘apparent’.  相似文献   

13.
The interpretation of the gravity and magnetic fields from inclined dikes has been studied with artifical data contaminated by various noise components: base level, linear trend, and random noise. A Gaussian window was applied to the data prior to transformation to reduce the influence of noise as demonstrated by an analysis of the horizontal cylinder. The case of the dike is more complicated due to the fact that its spectrum has a number of zeroes at wavenumbers which are inversely related to the width of the dike. Around these wavenumbers, especially the random noise distorts the spectrum making interpretation ambiguous.  相似文献   

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

15.
Full Tensor Gravity Gradiometry (FTG) data are routinely used in exploration programmes to evaluate and explore geological complexities hosting hydrocarbon and mineral resources. FTG data are typically used to map a host structure and locate target responses of interest using a myriad of imaging techniques. Identified anomalies of interest are then examined using 2D and 3D forward and inverse modelling methods for depth estimation. However, such methods tend to be time consuming and reliant on an independent constraint for clarification. This paper presents a semi‐automatic method to interpret FTG data using an adaptive tilt angle approach. The present method uses only the three vertical tensor components of the FTG data (Tzx, Tzy and Tzz) with a scale value that is related to the nature of the source (point anomaly or linear anomaly). With this adaptation, it is possible to estimate the location and depth of simple buried gravity sources such as point masses, line masses and vertical and horizontal thin sheets, provided that these sources exist in isolation and that the FTG data have been sufficiently filtered to minimize the influence of noise. Computation times are fast, producing plausible results of single solution depth estimates t hat relate directly to anomalies. For thick sheets, the method can resolve the thickness of these layers assuming the depth to the top is known from drilling or other independent geophysical data. We demonstrate the practical utility of the method using examples of FTG data acquired over the Vinton Salt Dome, Louisiana, USA and basalt flows in the Faeroe‐Shetland Basin, UK. A major benefit of the method is the ability to quickly construct depth maps. Such results are used to produce best estimate initial depth to source maps that can act as initial models for any detailed quantitative modelling exercises using 2D/3D forward/inverse modelling techniques.  相似文献   

16.
边界识别是重磁数据解释中的常用方法之一,依据其结果可划分出地质体的水平范围。边界识别结果受地质体埋深及导数计算误差的影响所识别边界与真实边界之间存在一定的差距,且边界识别法无法直观地给出地质体的深度信息。为了获得异常体的水平位置和深度信息,本文提出空间归一化边界识别方法,其对不同深度的边界识别函数进行归一化计算,空间归一化边界识别法的最大值对应于异常体的水平位置和深度。常规边界识别结果的误差随理深的减小而减小,而空间归一化边界识别法是通过最大值来判断地质体的位置,最大值是在地质体处获得,因此归一化边界识别方法所获得的结果是准确的。通过理论模型试验证明归一化边界识别方法能有效地完成异常体的水平位置和深度的计算,所获得的水平位置和深度信息与理论值相一致,为下一步的勘探计划提供了更加可靠的依据。将其应用于实际航磁数据的解释,获得了断裂的具体分布形式。  相似文献   

17.
Procedures are formulated using the correlation factors between successive least-squares residual magnetic anomaly profiles due to long horizontal cylinders for interpreting the three principal anomalies (vertical, horizontal, and total). It is demonstrated that correlation values can be used to determine the depth to the center of the buried structure and the index parameter. Procedures are also formulated to estimate the amplitude coefficient. Two worked examples using theoretical data show the effectiveness of the present method.  相似文献   

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

19.
The way potential fields convey source information depends on the scale at which the field is analysed. In this sense a multiscale analysis is a useful method to study potential fields particularly when the main field contributions are caused by sources with different depths and extents. Our multiscale approach is built with a stable transformation, such as depth from extreme points. Its stability results from mixing, in a single operator, the wavenumber low‐pass behaviour of the upward continuation transformation of the field with the enhancement high‐pass properties of n‐order derivative transformations. So, the complex reciprocal interference of several field components may be efficiently faced at several scales of the analysis and the depth to the sources may be estimated together with the homogeneity degrees of the field. In order to estimate the source boundaries we use another multiscale method, the multiscale derivative analysis, which utilizes a generalized concept of horizontal derivative and produces a set of boundary maps at different scales. We show through synthetic examples and application to the gravity field of Southern Italy that this multiscale behaviour makes this technique quite different from other source boundary estimators. The main result obtained by integrating multiscale derivative analysis with depth from extreme points is the retrieval of rather effective information of the field sources (horizontal boundaries, depth, structural index). This interpretative approach has been used along a specific transect for the analysis of the Bouguer anomaly field of Southern Apennines. It was set at such scales, so to emphasize either regional or local features along the transect. Two different classes of sources were individuated. The first one includes a broad, deep source with lateral size of 45∼50 km, at a depth of 13 km and having a 0.5 structural index. The second class includes several narrower sources located at shallowest depths, ranging from 3–6 km, with lateral size not larger than 5 km and structural indexes ranging from 1–1.5. Within a large‐scale geological framework, these results could help to outline the mean structural features at crustal depths.  相似文献   

20.
In regions where active source seismic exploration is constrained by limitations of energy penetration and recovery, cost and logistical concerns, or regulatory restrictions, analysis of natural source seismic data may provide an alternative. In this study, we investigate the feasibility of using locally‐generated seismic noise in the 2–6 Hz band to obtain a subsurface model via interferometric analysis. We apply this technique to three‐component data recorded during the La Barge Passive Seismic Experiment, a local deployment in south‐western Wyoming that recorded continuous seismic data between November 2008 and June 2009. We find traffic noise from a nearby state road to be the dominant source of surface waves recorded on the array and observe surface wave arrivals associated with this source up to distances of 5 kms. The orientation of the road with respect to the deployment ensures a large number of stationary points, leading to clear observations on both in‐line and cross‐line virtual source‐receiver pairs. This results in a large number of usable interferograms, which in turn enables the application of standard active source processing methods like signal processing, common offset stacking and traveltime inversion. We investigate the dependency of the interferograms on the amount of data, on a range of processing parameters and on the choice of the interferometry algorithm. The obtained interferograms exhibit a high signal‐to‐noise ratio on all three components. Rotation of the horizontal components to the radial/transverse direction facilitates the separation of Rayleigh and Love waves. Though the narrow frequency spectrum of the surface waves prevents the inversion for depth‐dependent shear‐wave velocities, we are able to map the arrival times of the surface waves to laterally varying group and phase velocities for both Rayleigh and Love waves. Our results correlate well with the known geological structure. We outline a scheme for obtaining localized surface wave velocities from local noise sources and show how the processing of passive data benefits from a combination with well‐established exploration seismology methods. We highlight the differences with interferometry applied to crustal scale data and conclude with recommendations for similar deployments.  相似文献   

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

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