大地电磁资料精细处理和二维反演解释技术研究(一)——阻抗张量分解与构造维性分析

本文对大地电磁观测阻抗实施一种数学变换--"共轭阻抗变换",发现转换后的观测阻抗和区域阻抗之间存在特定的关系,这种特定关系不受电场局部畸变的影响,而且同样不需要关于地下区域结构维性的假设,在区域结构是三维的情况下也是成立的.对转换后的观测阻抗采用已有的Swift旋转方法即可求得区域主轴方位角,然后采用最优化方法求取区域阻抗相位、振幅以及畸变因子.同时根据转换后的观测阻抗重新定义了不受电场局部畸变影响的构造维性参数.采用合成理论数据验证了新算法的正确性,和Swift、Bahr、GB、相位张量、WAL方法进行了对比分析,并将新方法应用于实测资料的解释,发展了一套消除局部畸变和进行构造维性分析的MT精细资料处理技术.进一步的工作是将新算法推广到多点多频,并采用最优化技术分解,发展一种稳定性好、多测点多频点的阻抗张量分解技术.     

A first arrival detection method for low SNR microseismic signal

Most of the microseismic signals have low signal-to-noise ratio (SNR) due to the strong background noise, which makes it difficult to locate the first arrival time. Both accuracy and stability of conventional methods are poor in this situation. To overcome this problem, here we proposed a new method based on the adaptive Morlet wavelet and principal component analysis process in wavelet coefficients matrix. The three components of microseismic signal make it possible to extract the features in wavelet coefficients domain. Then the reconstructed signal from weighted features presents an obvious first arrival. Tests on synthetic signals and real data provide a solid evidence for its feasibility in low SNR microseismic signal.     

广义S变换及其在大地电磁测深数据处理中的应用

广义S变换是一种优于短时窗傅立叶变换和小波变换的时频分析方法,利用广义S变换能够准确定位大地电磁资料中存在的噪声,通过定义时频窗对噪声进行滤除,从而明显提高阻抗视电阻率与相位的估算质量.本文基于S变换和大地电磁测深资料处理的基本原理,研究了基于广义S变换的大地电磁测深资料的处理流程和方法.对理论模拟信号及实测大地电磁场时间序列数据的处理,证实了方法的有效性.     

Wavelet analysis on transient behaviour of tidal amplitude fluctuations observed by meteor radar in the lower thermosphere above Bulgaria

On the basis of bispectral analysis applied to the hourly data set of neutral wind measured by meteor radar in the MLT region above Bulgaria it was demonstrated that nonlinear processes are frequently and regularly acting in the mesopause region. They contribute significantly to the short-term tidal variability and are apparently responsible for the observed complicated behavior of the tidal characteristics. A Morlet wavelet transform is proposed as a technique for studying nonstationary signals. By simulated data it was revealed that the Morlet wavelet transform is especially convenient for analyzing signals with: (1) a wide range of dominant frequencies which are localized in different time intervals; (2) amplitude and frequency modulated spectral components, and (3) singular, wave-like events, observed in the neutral wind of the MLT region and connected mainly with large-scale disturbances propagated from below. By applying a Morlet wavelet transform to the hourly values of the amplitudes of diurnal and semidiurnal tides the basic oscillations with periods of planetary waves (1.5/20 days), as well as their development in time, are obtained. A cross-wavelet analysis is used to clarify the relation between the tidal and mean neutral wind variability. The results of bispectral analysis indicate which planetary waves participated in the nonlinear coupling with the atmospheric tides, while the results of cross-wavelet analysis outline their time intervals if these interactions are local.     

基于双谱幅值和相位重构的地震子波提取

利用地震记录双谱中包含子波的幅值和相位信息,以及其超强的抗噪声干扰能力,采用一种基于双谱幅值和相位重构的地震子波提取方法,首先提取出子波幅值及相位信息,进而通过傅立叶反变换,使子波得以完全恢复.本文针对双谱相位重构递推公式,提出一种新的初值选取方法,使地震子波估计的稳定性得到了提高.仿真实验证实了该方法的可行性.     

Measurement of interstation phase velocity by wavelet transformation

In this paper,we present wavelet transformation method to measure interstation phase velocity. We use Morlet wavelet function as mother wavelet to filter two seismograms at various period of interest,and correlate the wavelet filtered seismograms to form cross-correlogram. If both wavelet filtered signals are in phase at that period,the phase of the cross-correlogram is a minimum. Using 3-spline interpolation to transform cross-correlation matrix to a phase velocity verse period image,it is convenient for u...     

Advanced methods of spectral analysis of quasiperiodic wave-like processes in the ionosphere: Specific features and experimental results

Backgrounds of adaptive Fourier transform and analytical wavelet transform have been briefly described in comparison with traditional Fourier transform using a time window. As an example, all three transforms are used to analyze quasiperiodic wave-like processes in the ionosphere, which accompanied the passage of the solar terminator and rocket launch from the Plesetsk site. The advantages of adaptive Fourier transform and analytical wavelet transform as compared to traditional Fourier transform, which make it possible to reliably detect wave-like disturbances against a background of noise at a signal-to-noise ratio not less than 0.1, have been demonstrated.     

Improvement of MT data processing using stationary and coherence tests

The great majority of data processing methods for magnetotelluric measurements are based on an average Fourier spectral analysis to estimate either the transfer function or the coherence function. This assumes that the time-series data are stationary. The aim of this paper is to present a preprocessing method able to extract the stationary and coherent part of the original signals. The practitioner can then apply the magnetotelluric method of his choice to these new data. This preliminary data sorting is done in four steps: (i) slow drifts are eliminated using a high-pass filter whose cut-off frequency is determined by an iterative procedure. Based on run tests, this procedure is also able to remove segments with non-independent samples in the time series; (ii) non-stationary segments are eliminated after band-pass filtering; (iii) non-coherent segments are eliminated before spectral analysis; (iv) the impedance tensor value is then retained, at a given frequency, only if the signals are coherent. This preprocessing method was tested on the simplest, but still used, magnetotelluric method which uses only two field components, and it was found that the average resistivity standard deviation decreased significantly from 14.6 Ωm without sorting to 8.6 Ωm after sorting.     

Wavelet power spectrum analysis of heterogeneities from sonic velocity logs

The continuous wavelet transform (CWT) is used to evaluate local variations in the power-law exponents of sonic log data. The resulting wavelet spectrum can be compared with the corresponding global estimates obtained by conventional Fourier transform methods. In Fourier analysis, the fundamental tool used to characterize a fluctuating velocity distribution is the power spectrum. It represents the energy contained in each wavenumber and thus provides information regarding the importance of each scale of heterogeneity. However, important spatial information regarding the location of events becomes implicit in the phase angle of the Fourier transform. In this paper, it is shown how the square of the amplitude of the wavelet transform is related to the Fourier spectrum and how spatial information can be expressed in an explicit manner. Using the conservation of energy, it is shown that the average wavelet power spectrum over the total depth range is equal to the global power spectrum. A Gaussian wavelet is chosen to realize the wavelet transform. Two synthetic sonic logs with exponential and von Karman correlation functions are used to demonstrate the potential of the suggested analysis. Furthermore, the wavelet transform is applied to the KTB (Continental Deep Drilling Program) sonic log data. The wide range of applications of the CWT shows that this transform is a natural tool for characterizing the structural properties of underground heterogeneities. It offers the possibility of separating the multiscale components of heterogeneities.     

小波变换与地震信号特征分析

地震波的瞬时信号,如瞬时振幅、瞬时频率及瞬时相位等,是研究地球介质的重要参数。根据小波定义,对Morlet 小波进行了修正,并对修正后的小波形态进行了深入讨论。理论分析表明,小波变换效果受到整形参数、小波长度、中心频率、频带宽度及小波个数等参数的制约,特别是整形参数与小波中心频率及频带之间关系对小波变换起到决定性作用。在地震波信号的实际处理中,可选取恰当的整形参数,同时采用合适的小波中心频率以避免小波变换对信号产生的遗漏和冗余。文中给出了实际地震记录处理的示例。     

利用小波变换提取面波群速度

选取了几种常见的小波母函数,分别提取了同一理论下的面波数据的群速度,并与理论群速度进行对比,结果表明Morlet小波提取面波群速度的效果最好.此外,将Morlet小波与常用的多重滤波提取群速度的结果进行了比较,结果表明: ① 多重滤波法非常依赖高斯滤波系数α的取值,α的取值应随面波周期的增大而减小;② 在α取值得当的前提下,在20—35 s周期范围内多重滤波法提取面波群速度的相对误差比Morlet小波小,在周期大于35 s时,两者相对误差相近; ③ 合适的α值的选取需在不同周期段耗费大量时间进行大量试验,这说明多重滤波法不具备自适应性;而采用小波变换分析短周期信号时,时间窗变窄,频率窗变长,当分析长周期信号时,时间窗变长,频率窗变窄,具有对信号的自适应性,这是小波变换相比多重滤波法的最大优点.      

Hilbert-Huang变换在提取地震信号动力特性中的应用

H ilbert-Huang变换(HHT)是一种处理非线性、非平稳信号的新方法。它通过经验模态分解将信号分解为有限个固有模态函数,并对每个固有模态函数进行H ilbert变换得到H ilbert谱。本文将这种方法应用于地震信号动力特性的提取,有效地获得了信号能量的时频分布,量化提取了中心频率、瞬时相位、瞬时能量、H ilbert能量、最大振幅对应的时频分布等动力特性,并与Fourier变换、小波变换等进行了比较,显示了HHT的优势以及对于进一步实现结构分析和控制的重要意义。     

基于压缩感知的信号重建方法及在气枪震源信号处理中的应用

在气枪源探测过程中,由于各种干扰因素的影响,导致部分有效信号缺失或受随机干扰严重,为了重构出连续完整的数据,依据气枪源信号在傅里叶变换域中具有稀疏性的特点,构建了一种基于压缩感知(Compressive Sensing,简称CS)的缺失信号重建方法。首先进行数值模拟,并将该方法与传统的插值方法处理效果进行对比,对重建效果进行均方根误差及信噪比(SNR)分析,结果显示:压缩感知方法重建前后信号的波形吻合度高、振幅一致性更强、同相轴清晰且连续性相同,对噪声压制较好。综上表明:该方法重建效果优于传统的三次样条插值方法。将该方法应用于实际资料,结果显示受干扰的有效信号能够得到很好的恢复重建。     

Variations in the amplitude and phase of VLF radiowaves in the ionosphere during the August 1, 2008, solar eclipse

Time variations in the amplitude and phase of signals of the Russian telecommunication station (the frequency is 25 kHz) on the Arkhangelsk—Kharkov path with a length of about 1600 km on the day of the August 1, 2008 solar eclipse (SE) and on the adjacent days are analyzed. Two types of effects are detected. An increase of the signal amplitude by approximately 32% in comparison with the background days and the 2.1 μs time shift of the signal during 2—2.5 h is referred to the first type. Changes in the spectral composition of the quasiperiodic disturbances in the ionosphere presented the second type of the effects. For spectral analysis of the quasiperiodic variations in the amplitude and phase of the radio signal, the window Fourier transform, adaptive Fourier transform, and wavelet transformation were applied simultaneously. In the period of SE and after it, oscillations with periods of 10—15 min (according to the amplitude data) and also about 10 and 18 min (according to the phase data) were intensified. Based on radio signal characteristics, the parameters of ionospheric disturbances are estimated.     

时频域相位滤波在远震接收函数噪声压制中的应用

宽频带地震观测数据中有效信号和干扰噪声经常发生混频效应,常规的频率域滤波方法很难将二者分离.地震波信号属于时变非平稳信号,时频分析方法能够同时得到地震波信号随着时间和频率变化的振幅和相位特征,S变换是其中较为高效的时频分析工具之一.本文以S变换为例,提出了基于相位叠加的时频域相位滤波方法.与传统叠加方法相比,相位叠加方法对强振幅不敏感,对波形一致性相当敏感,更加利于有效弱信号信息的检测.时频域相位滤波方法滤除与有效信号不相干的背景噪声,保留了相位一致的有效信号成分,显著提高了信噪比.运用理论合成的远震接收函数数据和实际的宽频带地震观测数据检验结果显示该方法较传统的带通滤波方法相比,即使在信噪较低且混频严重条件下,时频域相位滤波方法的滤波效果依然很明显,有助于识别能量较弱的有效信号.     

复小波包域的微震信号分析提取方法

提出了一种在复小波包域分析提取微震信号的新方法。此方法先采用复小波包分解同时获得信号的幅值信息和相位信息,再采用一种称为“比值加权”的复小波包重构法,重构出信噪比大大增强的微震信号。实际应用表明这种新方法在提取微震信号,提高信噪比方面效果显著。     

线性SDOF体系共振动力响应的Hilbert谱分析

利用HHT方法分析了线性SDOF体系在共振状态下的强迫动力响应,系统地研究了线性SDOF体系在简谐波、线性调频波与正弦调频波输入下处于共振状态时结构动力响应的H ilbert谱与本征模态函数的特征,并且将H ilbert谱、Morlet小波谱及Fourier幅值谱所蕴含的物理意义及分辨率特性进行了比较。研究结果表明,与Fourier幅值谱相比,H ilbert谱在刻画非平稳时程的某些时变特征方面更加直观,其物理意义更加明确,而与Morlet小波谱相比H ilbert谱则具有较高的分辨率。     

山西阳高～河北容城剖面大地电磁资料的初步解释──阻抗张量分解技术及其应用

介绍了阻抗张量分解技术及其在阳高～容城剖面的大地电磁资料解释中的应用。在某些测点局部电流畸变的影响较严重，用张量分解技术对资料分析后表明，沿剖面的区域异常总体上具二维性质，分解技术处理后的视电阻率和相位资料与地质构造有较好的对应性。事实表明，阻抗张量分解技术是对大地电磁资料进行分析解释的有力工具