An improved spectral turning-bands algorithm for simulating stationary vector Gaussian random fields

We propose a spectral turning-bands approach for the simulation of second-order stationary vector Gaussian random fields. The approach improves existing spectral methods through coupling with importance sampling techniques. A notable insight is that one can simulate any vector random field whose direct and cross-covariance functions are continuous and absolutely integrable, provided that one knows the analytical expression of their spectral densities, without the need for these spectral densities to have a bounded support. The simulation algorithm is computationally faster than circulant-embedding techniques, lends itself to parallel computing and has a low memory storage requirement. Numerical examples with varied spatial correlation structures are presented to demonstrate the accuracy and versatility of the proposal.     

平稳地震动过程具分段均布随机频率的谱拟合模型

本文首先针对平稳随机地震动谱拟合问题建立了一个统一的谱拟合模型,在此基础上,提出了一新的具分段均布随机频率的拟合模型,并与几类常用的谱拟合模型进行了比较研究,指出拟合模型的频率特性对拟合功率谱与反应谱具有显著的影响。算例表明,当结构的自振频率与拟合模型的分量频率接近时,计算的响应值明显偏大,对于小阻尼情况,这一现象更为明显。     

Response of base-isolated buildings to random excitations described by the Clough-Penzien spectral model

The stationary response of base-isolated buildings subjected to earthquake excitation is studied. The frequency content of earthquake input is described by the Clough-Penzien spectral model. The response parameters of interest are (1) the root-mean-square (RMS) displacement σ_x of the basement relative to the foundation (i.e. shear deformation of the isolation system) and (2) the ratio (σ_a/σ_a0) of the RMS value of the absolute acceleration at the roof of the isolated structure over the corresponding value when the isolation system is locked. The variation of these response parameters with the effective frequency f₀ of the base-isolated structure is investigated. As input, earthquakes with moment magnitudes M = 7-3 and M = 6-0 are considered. The acceleration spectra corresponding to these two earthquake sizes have pronouncedly different frequency content over the frequency range 0–1-1–0 Hz which is of primary importance for base-isolated structures. An important conclusion that comes from these analyses is that confidence in the effectiveness of a base-isolated system should be based primarily on its capacity to absorb/dissipate energy and less on its influence in shifting the fundamental period of the structure out of the range of dominant earthquake energy.     

Unification of periodic fluctuations in high latitude geophysical phenomena

The oscillatory nature of wave-particle interaction is proposed to be the source of non-local fluctuations in the electron energy. The energized electron flux when precipitated in the lower ionosphere gives rise to various geophysical phenomena. The time scale of energy oscillation is estimated and its validity is discussed with reference to observed time scales of various geophysical phenomena.     

Enhancement of radon signals in geophysical studies with the track technique

The existence of a suspected geological fault has been confirmed using Solid State Nuclear Track Detectors (SSNTDs) by measuring radon concentration variations in the upper soil above its inferred position. The results obtained prompted us to increase the natural radon signal in the soil, using an additional radon source; this enhancement technique, has been experimentally checked with SSNTD detectors.On leave from Faculté des Sciences, Laboratoire des D.S.T.N. Université de Dakar, Dakar-Fann, Sénégal     

Wavelet-vaguelette decomposition of spatiotemporal random fields

A wavelet-based orthogonal decomposition of the solution to stochastic differential/pseudodifferential equations of parabolic type is derived in the cases of random initial conditions and random forcing. The family of spatiotemporal models considered can represent anomalous diffusion processes when the spatial operator involved is a fractional or multifractional pseudodifferential operator. The results obtained are applied to the generation of the sample paths of Gaussian spatiotemporal random fields in the family studied.     

基于褶积模型的地球物理反演模型增强

地球物理数据在采集和处理过程中,由于存在噪声、模型误差、以及数据离散化误差等系统误差,导致了异常体边界模糊和模型分辨率降低等一些不可避免的不良系统退化效应的产生.本文提出了一种新的地球物理反演模型增强方法,通过消除反演估计模型中的系统误差,压制模型中的不良系统退化效应,增强反演模型的分辨率.文章从理论上分析了数据中存在的系统误差对模型求解的影响,提出了一个新的系统误差褶积退化模型,并根据该模型提出了一种基于混合范数总变分正则化的盲反褶积模型增强算法.最后,文章通过1D线性反演增强试验和2D大地电磁反演增强试验,验证了所提出的地球物理系统退化模型的正确性,以及盲反褶积增强算法的有效性.试验结果表明,方法可以有效地提高反演参数模型的分辨率.     

基于频谱分解的碳酸盐岩储层识别

频谱分解技术是将地震信号从时间域转换到频率域,分析频率对不同尺度地质体的振幅、相位响应特征的一项技术.频谱分解能够得到高于传统分辨率的解释结果,提高刻画储层分布的能力.本文详述了短时傅里叶变换、连续小波变换和S变换的数学原理及适用性:短时傅里叶变换使用固定时窗,不能根据信号的变化调整分辨率,只适合分析分段平稳或近似平稳的信号;连续小波变换使用移动的、尺度可变的小波作为时窗,具有多分辨率特点,但是实际中选择能反映信号特征的小波函数不易;S变换使用频率的倒数来调节时窗,具有多分辨率特征,对数据处理的适应性较强.将这三种方法分别应用于碳酸盐岩储层发育区,利用靠近地震主频的35 Hz分频剖面,分析了不同时窗大小的短时傅里叶变换效果,不同类型小波的连续小波变换效果,并对比了不同频谱分解算法对储层的描述精度.通过分析得出分频剖面比常规地震剖面更有利于储层识别,且S变换效果最好.     

Modal decomposition method for stationary response of non-classically damped systems

The stationary response of multi-degree-of-freedom non-classically damped linear systems subjected to stationary input excitation is studied. A modal decomposition procedure based on the complex eigenvectors and eigenvalues of the system is used to derive general expressions for the spectral moments of response. These expressions are in terms of cross-modal spectral moments and explicitly account for the correlation between modal responses; thus, they are applicable to structures characterized with significant non-classical damping as well as structures with closely spaced frequencies. Closed form solutions are presented for the important case of response to white-noise input. Various quantities of response of general engineering interest can be obtained in terms of these spectral moments. These include mean zero-crossing rate and mean, variance and distribution of peak response over a specified duration. Examples point out several instances where non-classical damping effects become significant and illustrate the marked improvement of the results of this study over conventional analysis based on classical damping approximations.     

地震复谱分解技术及其在烃类检测中的应用

谱分解技术在地震解释领域已得到广泛应用,但常用的谱分解方法存在两方面的不足.一是时间分辨率低,难以对薄层进行刻画;二是在烃类检测中多解性强,难以区分流体类型.为了改善该问题,本文提出一种基于地震复谱分解技术的烃类检测方法.复谱分解是指用一个包含多个不同频率Ricker子波的复子波库对地震道进行分解,从而得到时变子波频率和相位信息的过程.借助稀疏反演技术复谱分解可以获得高分辨率的时频能量谱和时频相位谱.本文首先通过拟合算例验证了复谱分解方法刻画薄层的能力以及求取子波频率和相位的准确性.然后利用基于Kelvin-Voigt模型的黏弹波动方程数值模拟对衰减引起子波相位改变的原因进行了分析.最后通过实际资料应用展示了本文方法在储层预测中的高时间分辨率优势,验证了利用子波相位信息识别气藏的有效性.     

lp norm inverse spectral decomposition and its multi-sparsity fusion interpretation

Spectral decomposition has been widely used in the detection and identification of underground anomalous features(such as faults,river channels,and karst caves)...     

On estimation of high frequency geophysical signals in Earth rotation by complex demodulation

Polar motion and UT1 contain small geophysical signals (<0.1 milliarcsecond) driven by the daily cycle in the solar heating. These signals are either irregular or quasi-harmonic therefore need to be expressed by the time series and studied using the time domain techniques. This work is focused on the complex demodulation method which seems to be particularly useful for studying high frequency geophysical signals in Earth rotation. After general introduction of the method and its properties, we describe its application for modeling diurnal and subdiurnal components of polar motion and UT1 and of the corresponding excitation functions. We also derive dynamical equations which can be used for direct comparison of the signals demodulated from the Earth rotation and geophysical excitation data.     

平稳随机激励下随机结构非线性动力响应分析

为了在结构强震灾变过程模拟中正确反映结构参数的随机性并准确界定结构失效模式,提出一种平稳随机激励下随机结构非线性动力响应分析的简化方法.首先运用点估计法将复合非线性随机系统转化为一系列平稳随机激励下的确定性非线性结构,然后再应用成熟的等效线性化方法计算出平稳随机激励下非线性结构动力响应的各阶矩.文末进行了算例分析与验证...     

分频混色技术在高精度地震解释中的应用

分频解释技术已经越来越多的应用到地震储层预测和油藏描述中,成为地震解释的一种常规技术手段,分频解释技术中分频RGB颜色融合显示方法因其独有的视觉效果也开始受到关注.本文介绍了分频混色技术在高精度地震解释中的应用以及在应用过程中应注意的几点问题,包括优选高精度频谱分解方法、RGB颜色融合技术应用中频率选择的原则等.在此基础上,进一步探讨了分频数据RGB混色方法能有效解决地震沉积学中地层切片精度不满足薄储层预测的问题.RGB颜色融合显示方法能充分发掘高精度地震数据宽频的优势,得到精度更高的储层预测结果.     

二阶同步挤压S变换及其在地震谱分解中的应用

同步挤压S变换是一种处理非平稳信号的新时频变换方法,通过"挤压"信号的S变换结果,可以得到高分辨率的时频谱.但是,当信号的相位随时间呈非线性变化时,同步挤压S变换计算出的瞬时频率会出现误差,并造成时频谱分辨率降低.为了改善其对于该类信号的时频分析效果,我们利用时间和频率的二阶偏导数对瞬时频率计算式进行修正,提出二阶同步挤压S变换.合成信号处理结果表明,二阶同步挤压S变换的分辨率不但明显高于常用时频变换,在信号瞬时频率随时间呈二次或正弦变化的情况下,其时频挤压效果也好于同步挤压S变换.我们将二阶同步挤压S变换应用到天然气地震勘探资料的谱分解当中,结果表明,二阶同步挤压S变换可以很好地检测到与天然气相关的谱异常.因此,二阶同步挤压变换对于地震解释是一种很有潜力的方法.     

石墨-石英模型电导率的数值模拟研究及地球物理应用

地壳高导异常及各向异性成因是地学研究的热点之一,断裂带区域常呈现出高电导率特征,对于石墨富集的断裂带,石墨形态及分布对导电性的影响已经引起关注.本文利用有限元数值模拟方法,构建了石墨-石英三维模型,探讨了石墨含量、排列方式、孔隙度及温度对模型电导率的影响.结果表明,石墨连通前,温度的增加能使模型电导率显著增大,而在石墨连通后,石墨含量的增加、片状石墨沿导电方向的定向排列以及石墨直径与厚度比值的增大,均能导致模型电导率发生显著增加,孔隙度的增大则导致模型电导率减小.同时,这些参数的变化也会使模型电导率突变时的石墨含量阈值发生改变,当石墨排列方向改变时,这一阈值变化显著,对于石墨富集的一些断裂带区域,石墨定向排列可能是引起高导的一种成因.

     

S变换谱分解技术在深反射地震弱信号提取中的应用

在深反射地震资料处理中,当来自深部的有效弱信号和噪声干扰频带差异较小且难以区分时,传统滤波方法的应用会受到限制.谱分解方法是一种使用离散傅里叶变换,基于信号的频率-振幅谱等信息生成高分辨率地震图像的方法,通常用来识别介质物性横向分布特征,处理复杂介质内频谱变化和局部相位的不稳定性等问题,包括定位复杂断层和小尺度断裂等.S变换作为一种新的时频分析方法,具有自动调节分辨率的能力,近些年来被广泛应用到勘探地震、大地电磁等数据处理中,逐渐成为地球物理方法中噪声压制的有效方法之一.与常规石油反射地震资料相比,深反射主动源地震为了探测深部结构信息,常采用大药量激发方式、长排列观测系统等,导致深部有效信号基本湮灭在噪声干扰之中.针对深反射数据特点,本文结合谱分解和S变换技术,首先设计了简单的脉冲函数实验数据,证实S变换方法的有效性,同时说明谱分解方法的效果受所用时频分析方法影响较大,而其中决定分辨能力的变换窗函数的选取尤为重要.在此基础上,分别应用到深反射地震资料的单道和叠加剖面实际数据上,对比分析了传统变换谱分解和S变换谱分解的应用效果,单道资料对比结果表明:相比传统谱分解,S变换谱分解方法具有自动调节分辨率的能力,能够精确的标定深反射地震资料中弱信号不同时刻的频率分量; 叠加剖面资料应用结果表明:由S变换谱分解得到的剖面结果与其他谱分解方法结果整体上具有较高的一致性,同时清晰地刻画出原叠加剖面上被噪声湮灭的低频细节特征,提高了剖面的分辨率及同相轴连续性;对比结果明显看出,Gabor变换谱分解方法得到的结果同相轴较为破碎,分析原因认为这是由Gabor变换的时频分解方法的定长窗函数所致,窗口大小不会随着信号频率的变化来调节长度,只能在处理的过程中根据一定的记录长度范围选取窗函数参数,而S变换谱分解方法在窗函数的选取时,通过时变信号的局部频率特征自动调节窗口长度,能够更好的刻画各个频段的细节特征,在深反射剖面成像应用中效果尤为明显.本文结果表明S变换谱分解技术在深地震叠加剖面上的应用有效地提高了来自深部弱反射信号的信噪比和分辨率,并刻画出了叠加剖面上所不具有的低频细节特征,在实际深反射地震资料处理中能有效保护低频弱信号获得更好的成像效果.本文为深地震反射资料中弱信号的保护处理找到一种有效的方法.