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沉积层较厚的集中地区, 表现出了实际记录较衰减关系值偏大的现象, 认为同样是由于厚沉积层对地震动加速度反应谱长周期的放大作用导致的。 研究了震中距差别不大的情况下, 场地类型与沉积层厚度对反应谱特征周期的影响, 对比基岩台站与软弱地基土层台站的强震动记录反应谱, 发现软弱土层台站的土层对地震动有一定的放大作用, 导致中长周期地震动被放大, 对比位于沉积层较薄的隆起区台站与位于沉积层较厚的凹陷区台站强震动记录反应谱, 发现厚的沉积层不仅对反应谱长周期有放大的作用, 同时也会使得反应谱特征周期值变大。     

基于数学建模的潜在地震破裂面源检测方法研究

大多数地震破裂面源检测方法都是通过简化地震震源,将地震震源表示成线源或者点源,无法有效描述地震带地震破裂面源产状和大小,不适用地震震级较大的情况下地震危险性检测。因此提出基于数学建模的潜在地震破裂面源检测方法,在地震震级较大时仍能检测出地震危险性概率。选取适宜的地震基岩水平峰值加速度衰减关系,分析地震震级、破裂长度、破裂宽度相互关系,确定地震引起的潜在地震破裂面源大小,计算给定地震动小于在场点处产生地震动的概率,将该概率同地震动加速度衰减关系结合,得到地震动年超越概率,分析地震危险性。经过实验检测发现,所提方法检测出的年超越概率与峰值加速度、最大震级有关,该概率能精准表示地震带地震破裂面源产状和大小,说明该方法检测地震危险性是合理的。     

Attenuation Characteristics of Earthquake Ground Motion for Large Volume Airgun

In order to further deepen the understanding of seismic wave propagation characteristics induced by the large volume airgun source, experimental data from multiple fixed excitation points in Fujian Province were used to obtain the equivalent single excitation high signal-to-noise ratio velocity and displacement records through linear stacking and simulation techniques. Then the peak displacements of different epicentral distances were used to calculate the equivalent magnitude of the airgun source excitation at different fixed excitation points so as to establish the attenuation relationship between equivalent magnitude,epicenter distance and velocity peak. Our results show that:① Within 270 km of epicentral distance,for the large-volume airgun's single shot,the peak velocity range is about 700-4 nm/s,and the peak displacement range is about 200. 0-0. 2 nm;② The equivalent magnitude of the P-wave from the airgun source with a total capacity of 8,000 in 3 is 0. 181-0. 760,and the equivalent magnitude of the S-wave is 0. 294-0. 832. By contrast,the equivalent magnitude of the P-wave from the airgun source with a total capacity of 12,000 in 3 is 0. 533-0. 896,and the equivalent magnitude of the S-wave is 0. 611-0. 946. The S-wave energy is greater than the P-wave energy, and the excitation efficiency varies greatly with different excitation environment;③ The peak velocity increases with the equivalent magnitude,and decreases with the epicentral distance. The vertical component of the P-wave peak velocity is the largest among those three components,while the S-wave has the smallest vertical component and similar horizontal components. Hence,our research can provide an important basis for the quantitative judgment of the seismic wave propagation distance using the airgun and the design of the observation system in deep exploration or monitoring with airgun.     

地球内核顶部300km速度和衰减各向异性的区域变化

衰减结构是地球内核的重要性质,它可以与地球内核的速度结构结合,对内核的形成和演化机制提供更全面的信息.本文系统收集了1991年到2014年全球、区域和临时地震台网的PKPDF和PKPBC数据,研究了澳大利亚、非洲和太平洋中部下方内核顶部300km的速度和衰减各向异性结构.速度结果表明,澳大利亚下方内核的速度没有明显的各向异性,但是非洲和太平洋中部下方的内核具有明显的各向异性,且非洲的速度各向异性强于太平洋中部.同时,相对于AK135模型,澳大利亚的平均速度快0.5%,而非洲和太平洋中部的平均速度与参考模型没有明显差异.对于内核的衰减结构,我们得到以下结果:1)在东西方向,内核顶部200km左右的区域,澳大利亚的衰减最强(Q值在400左右),非洲和太平洋中部的Q值分别在600和500左右.2)澳大利亚下方的内核衰减没有明显的各向异性,非洲和太平洋中部下方的内核衰减存在明显的各向异性.此外,内核在非洲地区的衰减各向异性强于太平洋中部的各向异性.3)最后,内核中三个区域的速度和衰减具有良好的相关性,即高/低速对应于高/低衰减.考虑到以上结果以及三个区域的位置,我们认为内核顶部的速度和衰减结构都存在区域变化,而不是简单的半球变化.这种区域变化很可能是由于核幔边界热结构的不均一性和内核耦合,使得内核顶部的不同区域在形成过程中受不同的变形影响,从而形成铁晶体不同的生长和排列,引发了不同的各向异性特征.     

2018年2月12日河北永清4.3级地震强震动记录分析

2018年2月12日河北省廊坊市永清县发生4.3级地震,震中位于廊固凹陷内的河西务断裂附近,国家强震动台网中心获取74组强震动记录,本文对触发台站进行场地分类,结果显示大部分为Ⅰ类场地(51.35%),其次为Ⅱ类场地(35.14%)和Ⅲ类场地(13.51%)。根据5个典型台站加速度时程记录及反应谱分析本次地震强震动特征,并对各台站峰值加速度反应谱与方位角和震中距的关系进行分析。采用克里格插值方法绘制地震仪器烈度分布图,由于台站分布不均,导致极震区附近烈度影响场计算缺值现象,本文通过拟合本次地震强震动记录得到该地区地震动衰减关系,在空白区域建立空间随机假设台站进行补点插值,解决了计算缺值问题,为缺少台站记录的震区提供准确快速制作烈度分布图的思路,为震害调查和地震应急救援提供重要依据。     

蓖齿眼子菜的光合速率及影响因素

本文研究了蓖齿眼子菜(Potamogeon pectinatus L)在滇池的光合速率、体内叶绿素季节变化及光和水温对其光合速率影响,并根据实验,计算出了蓖齿眼子菜生长的光补偿点,研究结果表明,蓖齿眼子菜1～3月份的净光合产氧量最大,为1.82～1.83mg/(g·h),6月份最低,为0.63 mg/(g·h);植物体内Chl.a和Chl.b平均含量分别为0.25～1.15mg/(g·h)和0.28～0.85mg/(g·h),蓖齿眼子菜可以通过调节体内Chl.a/Chl.b值来适应不同季节生长;蓖齿眼子菜有较宽的水温适应范围,10℃时的净光合产氧量可达0.88mg/(g·h),最高值出现在水温25℃,为1.16mg/(g·h),而30℃以上的水温已表现出对其生长的不利影响;蓖齿跟子菜生长的光补偿点为358～1256Lx,并随温度上升而逐渐增加.     

Predictive optimal control for seismic analysis of non-linear and hysteretic structures

In this paper, an effective active predictive control algorithm is developed for the vibration control of non-linear hysteretic structural systems subjected to earthquake excitation. The non-linear characteristics of the structural behaviour and the effects of time delay in both the measurements and control action are included throughout the entire analysis (design and validation). This is very important since, in current design practice, structures are assumed to behave non-linearly, and time delays induced by sensors and actuator devices are not avoidable. The proposed algorithm focuses on the instantaneous optimal control approach for the development of a control methodology where the non-linearities are brought into the analysis through a non-linear state vector and a non-linear open-loop term. An autoregressive (AR) model is used to predict the earthquake excitation to be considered in the prediction of the structural response. A performance index that is quadratic in the control force and in the predicted non-linear states, with two additional energy related terms, and that is subjected to a non-linear constraint equation, is minimized at every time step. The effectiveness of the proposed closed-open loop non-linear instantaneous optimal prediction control (CONIOPC) strategy is presented by the results of numerical simulations. Since non-linearity and time-delay effects are incorporated in the mathematical model throughout the derivation of the control methodology, good performance and stability of the controlled structural system are guaranteed. Copyright © 1999 John Wiley & Sons, Ltd.