Modeling the seismic response of basin-edge structure by a sequential waveform method

A sequential waveform method is developed to simulate the seismic response of basin-edge structure excited by a plane incident P-wave. The full procedure involves: (a) a previous parameterization of the investigated model using the seismic wave velocities and depths of the sedimentary stratifications; (b) an input motion determined from the records at stations installed on hard rock; (c) forward computation of the P-SV elastic wave field by means of a two-dimensional finite difference (FD) method; (d) the optimization of the model vector using simulated annealing technique and comparing the simulated seismic response of the tested structure with the observed wave field; (e) the correction of the initial model by trial-and-error by testing the differences between synthetics and observed data, and (f) the final solution obtained by iteration using the conjugate gradient algorithm. The search of an optimal basin-edge model has been parallel processed by varying the shapes and velocities of strata on the basis of the fitting of relative timing, amplitude and phase between the output and the observed data. The input motion and sensitivity have been checked and the validity of the method has been demonstrated by numeric analysis. Using the teleseismic records generated by 7 earthquakes recorded at 26 broadband seismic stations, we have studied the seismic velocity structure of the southern edge of the Jiyang depression located in the Bohai Bay basin, northern China. Two cross sections show an agreement between the velocity results and the geological sections available in the region. In addition, we obtain evidence of three hidden faults under the sections and features that suggest major extensions at the Paleogene.     

Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.     

斯通利波和渗透率的关系分析

理论和实际资料证明,在渗透性地层,低频斯通利波出现的频散和能量衰减与地层渗透率有关,因此可利用斯通利波的传播速度和能量衰减等特性来估算地层渗透率。斯通利波的能量和渗透率负相关。因为斯通利波能量代表了斯通利波的有关信息,因此可以认定斯通利波和渗透率负相关。     

近百年中国东部夏季降水的时空变率

利用中国东部25°N以北28个站1880-1999年夏季季降水序列,用旋转复经验正交函数（RCEOF）方法,研究了中国东部地区百年干湿的时空演变规律。结果表明,夏季降水空间变率大值区依次为:长江中下游地区、淮河流域、江南、华北、西南及东北。除西南外的5个关键区大体上反映了从6月到8月夏季雨带自南向北椎进所滞留的地区。旋转空间位相分布揭示了长江中下游地区、江南、东北的旱涝异常主要表现为驻波振动特征;而淮河流域、华北、西南地区显示出降水异常信号具有部分的行波特征。尤其第4空间模显示出旱涝异常信号从东北南部可沿着黄淮下游传到长江下游地区。对于近百年中国东部地区夏季于湿变化,长江中下游地区、淮河流域、华北及东北四个地区都存在20-25年时间尺度的周期振荡;长江中下游地区及华北地区都存在准60年时间尺度的振荡周期;东北地区主要表现出36年时间尺度的振荡周期;淮河流域存在明显的70-80年时间尺度的振荡周期;华北地区存在的11年时间尺度的振荡周期恰好与太阳黑子活动的11年周期相一致。在年代际时间尺度（包括次年代际时间尺度）上,长江中下游、淮河流域及华北地区的夏季降水的变化与太阳活动有显著的正相关。     

Simulation of artificial random field considering the temporal-spatial variation

Introduction For the seismic design of special structures such as nuclear power station, marine platform, long-span bridge and dam, generally the time-history response analysis of the structure under seismic excitation is imperative, which was coded in most seismic design codes. The earthquake records suitable for the seismic situation and site condition are necessary to be used as the seismic input in the dynamic analysis of structures. As a result of the limited observational condition of st…