Two examples of spectral approach to source depth estimation in gravity and magnetics

Fourier transformation of gravity and magnetic anomalies from space to the frequency domain provides aready method for source depth estimation since progressively deeper sources are indicated by the lower frequency components of an anomaly. Two examples are presented here to demonstrate the efficacy of this approach. One example pertains to the interpretation of vertical intensity magnetic anomalies over a layered ultramafic body for estimating the thickness of relatively non-magnetic layer in it and map the dispositions of the layers. The other example is the estimation of crustal thickness in the Indian region from the Bouguer anomaly map. The spectral method is shown, by these two examples, to provide a rapid and elegant tool for the source depth estimation for magnetic and gravity data.     

Continuous wavelet transform,theoretical aspects and application to aeromagnetic data at the Huanghua Depression,Dagang Oilfield,China

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.     

Curie Point Depths of Albania Inferred from Ground Total Field Magnetic Data

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, z_t, 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, z₀, 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, z_b, was calculated from z_b=2z₀−z_t. 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.     

Depth Estimation from Gravity Data Using the Maximum Entropy Method (MEM) and the Multi Taper Method (MTM)

The depth determination from the gravity data in frequency domain is carried out using the classical fast Fourier transform (FFT) method utilizing scaling properties of ensemble of anomalous source. The problem of calculating power spectrum from the FFT is well described in the literature. Here, the application of other high-resolution methods of power spectrum calculation, such as maximum entropy method (MEM) and multi-taper method (MTM) are explored to estimate depth to anomalous sources. At the outset, the FFT, the MEM and the MTM are tested on synthetic gravity data, generated for different types of synthetic models and then all these methods are applied to the field gravity data of the Bengal basin. The MTM with scaling is found to be superior for providing the detailed subsurface information rather than the MEM and the FFT methods in the case of synthetic as well as field examples.     

Mapping of magnetic basement in Central India from aeromagnetic data for scaling geology

The Central Indian region has a complex geology covering the Godavari Graben, the Bastar Craton (including the Chhattisgarh Basin), the Eastern Ghat Mobile Belt, the Mahanadi Graben and some part of the Deccan Trap, the northern Singhbhum Orogen and the eastern Dharwar Craton. The region is well covered by reconnaissance‐scale aeromagnetic data, analysed for the estimation of basement and shallow anomalous magnetic sources depth using scaling spectral method. The shallow magnetic anomalies are found to vary from 1 to 3 km, whereas magnetic basement depth values are found to vary from 2 to 7 km. The shallowest basement depth of 2 km corresponds to the Kanker granites, a part of the Bastar Craton, whereas the deepest basement depth of 7 km is for the Godavari Basin and the southeastern part of the Eastern Ghat Mobile Belt near the Parvatipuram Bobbili fault. The estimated basement depth values correlate well with the values found from earlier geophysical studies. The earlier geophysical studies are limited to few tectonic units, whereas our estimation provides detailed magnetic basement mapping in the region. The magnetic basement and shallow depth values in the region indicate complex tectonic, heterogeneity, and intrusive bodies at different depths, which can be attributed to different thermo‐tectonic processes since Precambrian.     

Sensitivity analysis of the parameters of earthquake recurrence time power law scaling

The stability of the power law scaling of earthquake recurrence time distribution in a given space–time window is investigated, taking into account the magnitude of completeness and the effective starting time of aftershock sequences in earthquake catalogs from Southern California and Japan. A new method is introduced for sampling at different distances from a network of target events. This method allows the recurrence times to be sampled many times on the same area. Two power laws with unknown exponents are assumed to govern short- and long-recurrence-time ranges. This assumption is developed analytically and shown to imply simple correlation between these power laws. In practice, the results show that this correlation structure is not satisfied for short magnitude cutoffs (m _c = 2.5, 3.5, 4.5), and hence the recurrence time distribution departs from the power law scaling. The scaling parameters obtained from the stack of the distributions corresponding to different magnitude thresholds are quite different for different regions of study. It is also found that significantly different scaling parameters adjust the distribution for different magnitude thresholds. In particular, the power law exponents decrease when the magnitude cutoff increases, resulting in a slower decrease of the recurrence time distribution, especially for short time ranges. For example, in the case of Japan, the exponent p₂ of the power law scaling at large recurrence times follows roughly the relation: , where m _c is the magnitude cutoff. In case of Southern California, it is shown that Weibull distribution provides a better alternative fit to the data for moderate and large time scales.     

Depth determination from a non-stationary magnetic profile for scaling geology

The conventional spectral analysis method for interpretation of magnetic data assumes stationary spatial series and a white‐noise source distribution. However, long magnetic profiles may not be stationary in nature and source distributions are not white. Long non‐stationary magnetic profiles can be divided into stationary subprofiles following Wiener filter theory. A least‐squares inverse method is used to calculate the scaling exponents and depth values of magnetic interfaces from the power spectrum. The applicability of this approach is demonstrated on non‐stationary synthetic and field magnetic data collected along the Nagaur–Jhalawar transect, western India. The stationarity of the whole profile and the subprofiles of the synthetic and field data is tested. The variation of the mean and standard deviations of the subprofiles is significantly reduced compared with the whole profile. The depth values found from the synthetic model are in close agreement with the assumed depth values, whereas for the field data these are in close agreement with estimates from seismic, magnetotelluric and gravity data.     

DEPTH ESTIMATION OF MAGNETIC SOURCES BY MEANS OF FOURIER AMPLITUDE SPECTRA*

The Fourier transform of a square-shaped section of a magnetic survey, digitized in a square grid, forms a rectangular matrix of coefficients which can be condensed to a series of average amplitudes dependent only on their frequency and no longer on the direction of the respective partial waves. These average amplitudes together represent a spectrum which–when plotted in a semilogarithmic coordinate system (log amplitude versus frequency)–often shows straight segments which decrease with increasing frequency. By continuing the given field downwards these straight segments become horizontal at a certain depth, the so-called “white depth”. This white depth may be used as a first estimate for the depth of magnetic sources producing the respective part of the field. It is shown that the sources which correspond to such use of the white depth can be expected to be “randomly distributed with some positive autocorrelation”. As an example for such a depth estimation the interpretation of the aeromagnetic survey of NW-Germany by a relief in 8–16 km depth is given. The relief divides the subsurface in an upper nonmagnetic layer and a lower layer with magnetization M= 2 Am^?1.     

Study of effects of focal depth on the characteristics of Rayleigh waves using finite difference method

This paper presents a relationship between the focal depth in terms of Rayleigh-wave wavelength and the dominant frequency of Rayleigh waves generated in a homogeneous half-space. Rayleigh waves were simulated using a (2, 4) staggered grid P-SV wave finite difference algorithm with VGR-stress imaging technique as a free surface boundary condition. VGR is an acronym for vertical grid-size reduction. The simulated seismic responses using P-wave and SV-wave sources at different focal depths revealed Rayleigh-wave generation up to certain focal depth only for the considered frequency bandwidth. A shift of normalized spectral shape of Rayleigh wave towards lower frequency with increasing focal depth was inferred. Largest spectral amplitude was obtained in the wavelength for which the ratio of focal depth to the wavelength of Rayleigh wave was around 0.17 in the case of P-wave source and 0.9 in the case of SV-wave source. An exponential decrease of spectral amplitude of Rayleigh wave with the departure of the ratio of focal depth to Rayleigh wave wavelength from the above mentioned values was obtained.     

Curie-point Depth from Spectral Analysis of Magnetic Data in Erciyes Stratovolcano (Central TURKEY)

Curie-point depth and heat flow values of the Erciyes region are determined to identify the thermal regime of the Central Anatolia by applying the spectral analysis method to the magnetic anomaly data. To compute the spectrum of the data, the magnetic anomaly of the region is transformed into 2-D Fourier domain to attain the average Curie depth. This method is useful in determining the top boundary of magnetic anomaly sources and reveals the Curie depth as 13.7 km in the study area. The obtained results imply a high thermal gradient (42.3°C km^?1) and corresponding heat flow values (88.8 mWm^?2) in the research area. Using the temperature value measured at borehole drilled by the General Directorate of Mineral Research and Exploration of Turkey (MTA), the values for the thermal gradient and heat flow value were computed as 50.7°C km^?1, 106.5 mWm^?2. From the heat flow value, the Curie-point depth was determined as 11.4 km in this region. It is concluded from the obtained values that the region has very high geothermal potential caused by partial melting of the lower crust.     

Estimation of depth and magnetization of magnetic sources from magnetic data: The linearized continuous inverse problem for 21/2D structures

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 2_1/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 2_1/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.     

谱矩方法在磁源体深度反演中的应用研究

谱矩方法可以对数据的表面形貌做较为细致的描述.它以随机过程为理论基础,用各阶谱矩及统计不变量等具体的参数表征表面的几何形态,算术平均顶点曲率是一种基于四阶谱矩的统计不变量.通常,埋深不同的场源所引起的磁异常尺度不同,从曲率的角度来理解即为磁异常曲面的弯曲程度不同.因此,本文应用算术平均顶点曲率提取磁异常的几何信息,并将所提取的信息用于场源深度的反演.理论上推导了基于谱矩的球状磁源体和板状磁源体的反演公式,得到了场源深度与磁异常、曲率之间的关系式.结合理论模型计算验证了方法的有效性,并与欧拉反褶积方法进行对比.与传统的方法相比,该方法快速简单,无需调节参数,且有较好的反演精度.最后,将该方法用于塔里木盆地航磁异常的反演和解释中,反演出的磁源体深度可满足区域磁异常数据分析和解释的要求,为克拉通沉积盆地磁异常源的深度划分提供丰富的信息.     

利用基于子波估计的频谱校正方法提高Q值估计精度(英文)

用Q值刻画的地震衰减在地震信号处理和解释中具有很广泛的应用。利用反射地震资料进行Q值估计需要解决地震子波和反射系数序列耦合的问题。从反射地震资料中去除反射系数序列的影响,这个过程称为频谱校正。本文提出了一种基于子波估计的求取Q值的方法,进而设计了一个反Q滤波器。该方法利用反射地震资料的高阶统计量进行子波估计,并利用所估计子波实现频谱校正。我们利用合成数据实验给出了质心频移法与频谱比法这两种常用的Q值估计方法在不同参数设置下的性能。人工合成数据和实际数据处理表明,利用本文提出的方法进行频谱校正后,可以得到可靠的Q值估计。经过反Q滤波,地震数据的高频部分得到了有效地恢复。     

The effects of anisotropy of marine magnetic anomalies on the Curie point depth estimates from spectral analysis

Azimuthally averaged power spectra are widely used in the Curie point depth (CPD) estimation with the implicit assumption that the magnetization distribution is random and uncorrelated. However, the marine magnetic anomalies are caused by bands of normal and reverse magnetization and show obvious trends. To investigate the effects of the anisotropy of marine magnetic anomalies on the CPD estimates, we develop 3D fractal striped magnetization models to produce lineated marine magnetic anomalies for the first time. We analyze the spectra anisotropy of the lineated magnetic anomalies of the synthetic fractal striped magnetization models and investigate its effects on the CPD estimates. The synthetic models and actual data show that the spectra of the lineated marine magnetic anomalies are directionally anisotropic. The amplitude response is strong and the slope of the logarithmic spectrum is large in a direction perpendicular to the stripes of magnetic anomalies, whereas the amplitude response is weak and the slope of the logarithmic spectrum is small in a direction parallel to the stripes of magnetic anomalies. The depth estimates in the perpendicular direction are close to the actual values, whereas the depths estimates in the parallel direction are significantly lower than the actual values. The actual marine magnetic anomalies of the South China Sea exhibit an anisotropic power spectrum that is consistent with the spectral anisotropy of magnetic anomalies of the synthetic fractal striped magnetization models.     

The search for weak harmonic signals in a spectrum with application to gravity data

When considering the search for discovery or amplitude estimation of a spectral line with a probabilistic approach, great attention must be paid to the meaning of each step. We give the probability law for the amplitude of a spectral peak in the presence of random noise appearing in a periodogram and discuss the effective probability of the existence of the corresponding wave. We find that the estimated amplitude of a spectral peak is biased and should be corrected when the signal-to-noise ratio is small. As a first application to gravity data, it results in a re-estimation of the gravimetric amplitude factors (delta factors) provided by least-squares tidal analysis. We also estimate the probability of observing a spectral line above a given level in the spectrum of a purely random noise. This allows us to compute for given spectrum the number of peaks expected to overcross the classical levels used in statistical analysis (like nσ, where σ is the standard deviation of the temporal noise distribution and n is an integer with typical values equal to 2 or 3). A specific application to real data is investigating the gravity spectrum derived from a 5 year record of the French superconducting gravimeter and we show that the predicted statistics are indeed in agreement with the observations. We also show the statistical consequence of using longer observing periods to obtain the spectral estimations. The problem of detecting translational motion of the Earth's solid inner core (Slichter modes) in a gravity spectrum is analyzed and the probabilities of having a triplet of random peaks thresholding specific levels in a given frequency window are computed. We show that, in the case of a typical gravity spectrum (1 year of hourly data and a frequency window of 0.03 cycle h⁻¹), the probability of having a random set of three peaks exceeding a level of 3 σ, is very high. This emphasizes the need for a very careful analysis of spectral lines before inferring the existence of a true physical signal.     

A numerical comparison of time and frequency‐domain marine electromagnetic methods for hydrocarbon exploration in shallow water

In shallow water the frequency domain controlled source electromagnetic method is subject to airwave saturation that strongly limits the sensitivity to resistive hydrocarbon targets at depth. It has been suggested that time‐domain CSEM may offer an improved sensitivity and resolution of these deep targets in the presence of the airwave. In order to examine and test these claims, this work presents a side‐by‐side investigation of both methods with a main focus on practical considerations, and how these effect the resolution of a hydrocarbon reservoir. Synthetic noisy data for both time‐domain and frequency domain methods are simulated using a realistic frequency dependent noise model and frequency dependent scaling for representative source waveforms. The synthetic data studied here include the frequency domain response from a compact broadband waveform, the time‐domain step‐response from a low‐frequency square wave and the time‐domain impulse response obtained from pseudo‐random binary sequences. These data are used in a systematic resolution study of each method as a function of water‐depth, relative noise and stacking length. The results indicate that the broadband frequency domain data have the best resolution for a given stacking time, whereas the time‐domain data require prohibitively longer stacking times to achieve similar resolution.     

Least-squares Minimization Approaches to Interpret Total Magnetic Anomalies Due to Spheres

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.     

A stochastic framework for downscaling processes of spatial averages based on the property of spectral multiscaling and its statistical diagnosis on spatio-temporal rainfall fields

Spectral multi-scaling postulates a power-law type of scaling of spectral distribution functions of stationary processes of spatial averages, over nested and geometrically similar sub-regions of the spatial parameter space of a given spatio-temporal random field. Presently a new framework is formulated for down-scaling processes of spatial averages, following naturally from the postulate of spectral multi-scaling, and key ingredients required for its implementation are described. Moreover, results from an extensive diagnostic study are presented, seeking statistical evidence supportive of spectral multi-scaling. Such evidence emerges from two sources of data. One is a 13 year long historical record of radar observations of rainfall in southeastern UK (Chenies radar), with high spatial (2 km) and temporal (5 min) resolution. The other is an ensemble of rain rate fields simulated by a spatio-temporal random pulse model fitted to the historical data. The results are consistent between historical and simulated rainfall data, indicating frequency-dependent scaling relationships interpreted as evidence of spectral multi-scaling across a range of spatial scales.     

金属矿山地地区地震勘探随机噪声的波动方程模拟

消减随机噪声是目前陆地地震勘探数据处理的关键问题之一,分析随机噪声的产生机制及特征是对其进行有效压制的先决条件.本文针对中国南方山地金属矿区的勘探环境,根据随机噪声中包含的自然噪声和人文噪声的发声机理分别确定其噪声源函数,以波动方程作为噪声传播模型对山地地区随机噪声进行建模,将随机噪声作为一个综合波场,并且与实际噪声记录进行比较.随机噪声记录作为时空域的二维随机过程,分别对模拟噪声和实际噪声记录的时间域波形(振动图)特征包括频谱、功率谱密度,相空间轨迹图,统计量特征(能量分布,累积分布,均值,方差,峰度,偏度),和空间域波形(波剖面)特征包括波数谱和统计量特征进行比较,对比结果显示在时空域模拟噪声和实际噪声都有基本相同的性质,证明了本文对随机噪声模拟方法的可行性,为进一步研究随机噪声时空域传播特性以及噪声消除奠定理论基础.