A seismic reflection from isotropic‐fractured fluid‐saturated layer

Average elastic properties of a fluid‐saturated fractured rock are discussed in association with the extremely slow and dispersive Krauklis wave propagation within individual fractures. The presence of the Krauklis wave increases P‐wave velocity dispersion and attenuation with decreasing frequency. Different laws (exponential, power, fractal, and gamma laws) of distribution of the fracture length within the rock show more velocity dispersion and attenuation of the P‐wave for greater fracture density, particularly at low seismic frequencies. The results exhibit a remarkable difference in the P‐wave reflection coefficient for frequency and angular dependency from the fractured layer in comparison with the homogeneous layer. The biggest variation in behaviour of the reflection coefficient versus incident angle is observed at low seismic frequencies. The proposed approach and results of calculations allow an interpretation of abnormal velocity dispersion, high attenuation, and special behaviour of reflection coefficients versus frequency and angle of incidence as the indicators of fractures.     

中国大陆背景噪声强度时空分布

使用中国地震台网888个宽频带地震台站2010年1月至2011年6月记录的连续波形垂直向分量, 进行了中国大陆背景噪声强度时空分布研究, 得到了中国大陆周期为10 s的背景噪声强度的时空分布图像。 背景噪声强度的时间分布具有一定的季节性和周期性变化特征, 研究表明2011年1~3月间的背景噪声强度最高。 背景噪声强度空间分布显示, 中国大陆背景噪声强度呈现明显的分区特性, 其分布与地质结构无明显相关性, 表明背景噪声场的强度高于场地效应。 其中, 东南沿海地区的背景噪声强度最大, 向内陆地区逐渐衰减, 至青藏高原地区背景噪声强度减至最小。 这种方向性的强度分布很可能与菲律宾和太平洋的潮汐等因素有关, 受印度洋的影响不大; 而到新疆地区, 背景噪声的强度又有一定程度增加, 表明这一地区的背景噪声主要不是来源于中国东南海域, 而是来自于欧亚大陆内部。 本研究得到的背景噪声强度可以为多个地区的微震强度提供基础数据资料, 也可以为背景噪声互相关的地震学研究提供参考。 更为重要的是, 还可以为基于背景噪声提取面波衰减研究提供基础数据, 进而实现从背景噪声中提取面波衰减。     

裂缝、裂隙介质衰减频散研究

为研究裂缝、裂隙介质中波致流引起的衰减,将裂缝看作背景孔隙岩石中非常薄且孔隙度非常高的层状介质,并等价成White周期层状模型.分别考虑不同类型的裂隙和孔隙之间的挤喷流影响,结合改进的Biot方程,推导得到裂缝裂隙介质的刚度与频率的关系.当缝隙中饱含流体时,介质的衰减和速度频散受裂缝、孔隙之间和裂隙、孔隙之间流体流动的显著影响.在低频极限下,裂缝裂隙介质的性质由各向异性Gassmann理论和挤喷流模型获得;而在非常高的频率时,由于缝隙中的压力来不及达到平衡,波致流的影响可忽略.分析表明,裂隙密度主要影响波的衰减,而裂隙纵横比主要控制优势衰减频率和速度显著变化的频率范围;由于不同裂隙的衰减机制不同,衰减和速度频散大小有所差异,但基本趋势相同.     

衰减夹层GPR模拟的时频域全局反射误差

全局反射误差分析是深入研究探地雷达(GPR)吸收边界条件吸收效率的有力工具.基于常规完全匹配层(PML)的标准交错网格有限差分算法必须满足严格的CFL条件限制,即在单位时间步长内,不容许电磁波传播的距离超过单元网格尺寸.为了提高主区域所有网格节点的计算效率,并有效地吸收传播后期出现的低频隐失波,提出基于非分裂递归卷积完全匹配层(UCPML)的旋转交错网格(RSG)GPR正演算法.该算法采用不同的网格交错策略,并在边界条件中引入了吸收低频隐失波的自由可变因子,使得该算法允许选取较大的时间步长,提高了计算效率,并且实现了对低频隐失波的高效吸收.本文首先给出了RSG差分格式,推导了基于UCPML的RSG差分更新方程,探讨了数值色散的稳定性条件,然后以绕射现象严重的衰减夹层数值模拟为例,分别从波场快照、单道波记录、时间域反射误差(TDRE)、频率域反射误差(FDRE)四个方面分析了UCPML与常规PML的全局反射误差,说明了基于UCPML和RSG的GPR正演算法能更有效地吸收低频隐失波.     

Nonstationary deconvolution using maximum kurtosis optimization

Deconvolution is an essential step for high-resolution imaging in seismic data processing. The frequency and phase of the seismic wavelet change through time during wave propagation as a consequence of seismic absorption. Therefore, wavelet estimation is the most vital step of deconvolution, which plays the main role in seismic processing and inversion. Gabor deconvolution is an effective method to eliminate attenuation effects. Since Gabor transform does not prepare the information about the phase, minimum-phase assumption is usually supposed to estimate the phase of the wavelet. This manner does not return the optimum response where the source wavelet would be dominantly a mixed phase. We used the kurtosis maximization algorithm to estimate the phase of the wavelet. First, we removed the attenuation effect in the Gabor domain and computed the amplitude spectrum of the source wavelet; then, we rotated the seismic trace with a constant phase to reach the maximum kurtosis. This procedure was repeated in moving windows to obtain the time-varying phase changes. After that, the propagating wavelet was generated to solve the inversion problem of the convolutional model. We showed that the assumption of minimum phase does not reflect a suitable response in the case of mixed-phase wavelets. Application of this algorithm on synthetic and real data shows that subtle reflectivity information could be recovered and vertical seismic resolution is significantly improved.     

建筑抗震实时混合试验时滞补偿方法研究综述

时滞补偿是建筑抗震实时混合试验的重要问题。文章对时滞补偿方法的研究进行总结,将时滞补偿方法分为基于时间的位移预测时滞补偿方法、基于反馈力的位移预测时滞补偿方法、基于加载系统模型求逆的时滞补偿方法、基于控制理论的时滞补偿方法和基于集成理念的综合补偿方法。     

Modelling of viscoacoustic wave propagation in transversely isotropic media using decoupled fractional Laplacians

A new wave equation is derived for modelling viscoacoustic wave propagation in transversely isotropic media under acoustic transverse isotropy approximation. The formulas expressed by fractional Laplacian operators can well model the constant-Q (i.e. frequency-independent quality factor) attenuation, anisotropic attenuation, decoupled amplitude loss and velocity dispersion behaviours. The proposed viscoacoustic anisotropic equation can keep consistent velocity and attenuation anisotropy effects with that of qP-wave in the constant-Q viscoelastic anisotropic theory. For numerical simulations, the staggered-grid pseudo-spectral method is implemented to solve the velocity–stress formulation of wave equation in the time domain. The constant fractional-order Laplacian approximation method is used to cope with spatial variable-order fractional Laplacians for efficient modelling in heterogeneous velocity and Q media. Simulation results for a homogeneous model show the decoupling of velocity dispersion and amplitude loss effects of the constant-Q equation, and illustrate the influence of anisotropic attenuation on seismic wavefields. The modelling example of a layered model illustrates the accuracy of the constant fractional-order Laplacian approximation method. Finally, the Hess vertical transversely isotropic model is used to validate the applicability of the formulation and algorithm for heterogeneous media.     

含孔隙、裂隙地层随钻多极子声波测井理论

储层岩石中普遍存在孔隙与裂隙,对钻井中的测井声波产生重要影响.基于孔、裂隙介质弹性波理论,导出了随钻声波测井的井孔声场表达式.据此考察了地层裂隙密度与含气饱和度的变化时井孔内随钻多极子模式波（斯通利波、弯曲波和螺旋波）的速度、衰减与灵敏度以及地层纵、横波的响应特征.裂隙密度与含气饱和度对模式波的速度频散与衰减都有影响,且两参数的值越大,影响越大.具体来说,速度对裂隙密度更敏感,而衰减对含气饱和度更敏感.具有"艾里相"特征的随钻偶极和四极子波在地层含气时产生强烈衰减,可以作为判断地层含气的一个明显指示.理论模拟与实际测井数据分析结果符合较好.     

Some remarks on Q-compensated sparse deconvolution without knowing the quality factor Q

The subsurface media are not perfectly elastic, thus anelastic absorption, attenuation and dispersion (aka Q filtering) effects occur during wave propagation, diminishing seismic resolution. Compensating for anelastic effects is imperative for resolution enhancement. Q values are required for most of conventional Q-compensation methods, and the source wavelet is additionally required for some of them. Based on the previous work of non-stationary sparse reflectivity inversion, we evaluate a series of methods for Q-compensation with/without knowing Q and with/without knowing wavelet. We demonstrate that if Q-compensation takes the wavelet into account, it generates better results for the severely attenuated components, benefiting from the sparsity promotion. We then evaluate a two-phase Q-compensation method in the frequency domain to eliminate Q requirement. In phase 1, the observed seismogram is disintegrated into the least number of Q-filtered wavelets chosen from a dictionary by optimizing a basis pursuit denoising problem, where the dictionary is composed of the known wavelet with different propagation times, each filtered with a range of possible values. The elements of the dictionary are weighted by the infinity norm of the corresponding column and further preconditioned to provide wavelets of different values and different propagation times equal probability to entry into the solution space. In phase 2, we derive analytic solutions for estimates of reflectivity and Q and solve an over-determined equation to obtain the final reflectivity series and Q values, where both the amplitude and phase information are utilized to estimate the Q values. The evaluated inversion-based Q estimation method handles the wave-interference effects better than conventional spectral-ratio-based methods. For Q-compensation, we investigate why sparsity promoting does matter. Numerical and field data experiments indicate the feasibility of the evaluated method of Q-compensation without knowing Q but with wavelet given.     

2012年唐山4.8级地震强震动记录特征及场地相关特性分析

2012年5月28日河北省唐山市古冶区与滦县交界发生4.8级地震, 国家强震动台网中心在河北、 天津和北京的94个强震动台站记录到了本次地震的加速度。 本文给出了获取记录的强震动台站分布及强震动记录结果, 统计了强震动记录数量随震中距的变化, 给出了3个较小震中距台站记录到的加速度时程; 绘制了空间地震动峰值加速度等值线图及周期0.2 s、 2.0 s加速度反应谱值的等值线, 发现峰值加速度等值线与长周期加速度反应谱等值线极值分布具有明显地域差异, 分析认为是由于厚沉积层对长周期地震动具有放大作用造成的。 通过强震动记录与适用于本区的三个衰减关系对比, 分析了此次地震的峰值加速度衰减特征, 同时研究了周期0.2 s、 2.0 s加速度反应谱值的衰减特征, 周期2.0 s反应谱值随震中距的衰减与衰减关系能较好地对应, 然而在震中距100～130 km沉积层较厚的集中地区, 表现出了实际记录较衰减关系值偏大的现象, 认为同样是由于厚沉积层对地震动加速度反应谱长周期的放大作用导致的。 研究了震中距差别不大的情况下, 场地类型与沉积层厚度对反应谱特征周期的影响, 对比基岩台站与软弱地基土层台站的强震动记录反应谱, 发现软弱土层台站的土层对地震动有一定的放大作用, 导致中长周期地震动被放大, 对比位于沉积层较薄的隆起区台站与位于沉积层较厚的凹陷区台站强震动记录反应谱, 发现厚的沉积层不仅对反应谱长周期有放大的作用, 同时也会使得反应谱特征周期值变大。