A comparative analysis of seismic response of shallow buried underground structure under incident P,SV and Rayleigh waves

In this study, A time-domain seismic response analysis method and a calculation model of the underground structure that can realize the input of seismic P, SV and Rayleigh waves are established, based on the viscoelastic artificial boundary elements and the boundary substructure method for seismic wave input. After verifying the calculation accuracy, a comparative study on seismic response of a shallow-buried, double-deck, double-span subway station structure under incident P, SV and Rayleigh waves is conducted. The research results show that there are certain differences in the cross-sectional internal force distribution characteristics of underground structures under different types of seismic waves. The research results show that there are certain differences in the internal force distribution characteristics of underground structures under different types of seismic waves. At the bottom of the side wall, the top and bottom of the center pillar of the underground structure, the section bending moments of the underground structure under the incidences of SV wave and Rayleigh wave are relatively close, and are significantly larger than the calculation result under the incidence of P wave. At the center of the side wall and the top floor of the structure, the peak value of the cross-sectional internal force under the incident Rayleigh wave is larger than the calculation result under SV wave. In addition, the floor of the underground structure under Rayleigh waves vibrates in both the horizontal and vertical directions, and the magnification effect in the vertical direction is more significant. Considering that the current seismic research of underground structures mainly considers the effect of body waves such as the shear waves, sufficient attention should be paid to the incidence of Rayleigh waves in the future seismic design of shallow underground structures.     

Rayleigh波入射下相邻山体地震响应边界元法模拟

采用间接边界元法（IBEM）,对Rayleigh波入射下两邻近山体的地震响应进行了定量分析。结果表明:与单个山体在Rayleigh波入射下的地震响应相比,两山地形的地震反应规律更为复杂,反应特征受控于入射波频率和山体间距等因素。总体上看:入射Rayleigh波频率较低时,两山间距对双山地形地震反应影响较大。低频波入射,受邻近山体影响,迎波面山体水平位移峰值约为入射波水平位移的6.3倍,放大效应可达单山作用的1.5倍。且山体山脚处竖向位移的频谱振荡更为剧烈,特定频率下山脚的竖向位移反应可达单山放大效应的1.6倍。较高频波入射时,迎波面一侧山体受邻近山体影响较小,且对波表现出明显的屏障效应,背波面山体地震反应强度被削弱。     

频率-波数域递推计算海底多分量地震记录中的反射系数

在频率-波数域中采用解析法,解出多层条件下海底实测的多分量地震数据分解成上行和下行P波和S波的算法,导出海底各层地震反射系数随入射角变化(简称RVA)的递推计算公式,为海底多波多分量AVO弹性参数的反演及流体因子预测提供基础数据.合成数据的计算结果表明,本文给出的算法能较可靠地从海底多波多分量记录中提取RVA信息.     

THE APPLICATION OF PML BOUNDARY CODITIONS IN THE SIMULATION OF SEISMIC WAVE BY THE HIGH-ORDER STAGGERED-GRID FINITE DIFFERENCES AT THE TRANSVERSELY ISOTROPIC MEDIA

In the realm of the numerical simulation, finite difference method and finite element method are more intuitive and effective than other simulation methods. In the process of simulating seismic wave propagation, the finite differences method is widely used because of its high computational efficiency and the advantage of the algorithm is more efficient. With the demand of precision, more and more researchers have proposed more effective methods of finite differences, such as the high-order staggered-grid finite differences method, which can restore the actual process of wave propagation on the premise of ensuring accuracy and improving the efficiency of operation. In the past numerical simulation of seismic wave field, different models of isotropic medium are mostly used, but it is difficult to reflect the true layer situation. With the research demand of natural seismology and seismic exploration, the research on anisotropic media is more and more extensive. Transversely isotropic(TI)media can well simulate the seismic wave propagation in the formation medium, such as gas-bearing sandstone, mudstone, shale et al., the character of TI media is reflected by introducing the Thomsen parameters to reflect its weak anisotropy of vertical direction by using Thomson parameter. Therefore, studying the process of seismic wave propagation in TI media can restore the true information of the formation to the greatest extent, and provide a more reliable simulation basis for the numerical simulation of seismic wave propagation. In the geodynamic simulation and the numerical simulation of the seismic wave field, under the limited influence of the calculation area, if no boundary conditions are added, a strong artificial boundary reflection will be generated, which greatly reduces the validity of the simulation. In order to minimize the influence of model boundaries on the reflection of seismic waves, it is often necessary to introduce absorbing boundary conditions. At present, there are three types of absorption boundary conditions: one-way wave absorption boundary, attenuation absorption boundary, and perfectly matched layer(PML)absorption boundary. In terms of numerical simulation of seismic waves, the boundary absorption effect of PML is stronger than the first two, which is currently the most commonly used method, and it also represents the cutting-edge development direction of absorption boundary technology. The perfectly matched layer absorbing boundary is effectively applied to eliminating the reflective waves from model boundaries, but for transversely isotropic medium, the effect of the absorbing is not very well. For this reason, the elastic dynamic wave equations in transversely isotropic media are derived, and we describe a second-order accurate time, tenth-order accurate space, formulation of the Madariaga-Virieux staggered-grid finite difference methods with the perfectly matched layer(PML)are given. In addition, we have established vertical transversely isotropic(VTI)media and arbitrary inclined tilted transversely isotropic(TTI)media models, using a uniform half-space velocity model and a two-layer velocity model, respectively. By combining the actual geoscience background, we set the corresponding parameters and simulation conditions in order to make our model more research-oriented. When setting model parameters, different PML thickness, incident angle, source frequency and velocity layer models were transformed to verify the inhibition of boundary reflection effect by PML absorption boundary layer. The implementations of this simulation show that the formula is correct and for the transversely isotropic(TI)media of any angular symmetry axis, when the thickness of the PML layer reaches a certain value, the seismic wave reflection effect generated by the artificial boundary can be well suppressed, and the absorption effect of PML is not subject to changes in incident angle and wave frequency. Therefore, the results of our study indicate that our research method can be used to simulate the propagation process of seismic waves in the transversely isotropic(TI)media without being affected by the reflected waves at the model boundary to restore the actual formation information and more valuable geological research.     

Variations in the dynamic response of structures founded on piles induced by obliquely incident SV waves

Although the seismic actions generally consist of a combination of waves, which propagates with an angle of incidence not necessarily vertical, the common practice when analyzing the dynamic behavior of pile groups is based on the assumption of vertically incident wave fields. The aim of this paper is to analyze how the angle of incidence of SV waves affects the dynamic response of pile foundations and piled structures. A three-dimensional boundary element-finite element coupling formulation is used to compute impedances and kinematic interaction factors corresponding to several configurations of vertical pile groups embedded in an isotropic homogeneous linear viscoelastic half-space. These results, which are provided in ready-to-use dimensionless graphs, are used to determine the effective dynamic properties of an equivalent single-degree-of-freedom oscillator that reproduces, within the range where the peak response occurs, the response of slender and nonslender superstructures through a procedure based on a substructuring model. Results are expressed in terms of effective flexible-base period and damping as well as maximum shear force at the base of the structure. The relevance and main trends observed in the influence of the wavefront angle of incidence on the dynamic behavior of the superstructure are inferred from the presented results. It is found that effective damping is significantly affected by the variations of the wave angle of incidence. Furthermore, it comes out that the vertical incidence is not always the worst-case scenario.     

Earthquake site effect modeling in the Granada basin using a 3-D indirect boundary element method

In this work, we develop the indirect boundary element method (IBEM) to simulate the seismic site response in a realistic, large-scale 3-D sedimentary basin. Most previous applications of boundary element method have used full-space Green’s functions for wave propagation between element points. We use half-space Green’s functions, which include the seismic wavefield interactions at the free surface and require only the boundary elements of the basin interface. In this way, the size of the matrix equation for solution in the IBEM can be reduced to approximately a quarter of that using full-space Green’s functions. The site response modeling of the Granada basin in southern Spain using the IBEM shows that the basin-induced scattering waves were identified as propagating back and forth inside the basin. The scattered waves also generate surface waves that are weakly propagated outside of the basin. The wave propagation inside and outside of the basin shows different patterns. We observe that the scattered wave is locally amplified, and its propagation direction deviates from that of the incident waves propagation direction. Therefore, the computed seismic response in the basin could provide us with good estimates of the seismic motion.     

基于混合波场地震动输入技术的近海场地地震反应分析方法

基于人工边界子结构模型,提出一种利用混合波场实现近海场地中地震P波和SV波垂直输入的方法.该方法中用于波动输入的混合波场由计算模型两侧截断边界的自由波场和底面边界的入射波场构成,避免了不规则近海场地的自由波场求解.同时采用基于声流体单元的流固耦合算法模拟场地-海水动力相互作用,利用流体介质人工边界和黏弹性人工边界单元模...     

地震波斜入射对高面板坝地震反应的影响

为研究地震波斜入射对高面板坝地震反应的影响,根据地震波动入射理论,采用FORTRAN进行波动荷载的编程计算,并与大型通用有限元软件ADINA相结合,实现基于黏弹性人工边界的地震波斜入射,研究P波和SV波分别以不同角度入射对高面板堆石坝地震反应的影响。结果表明,地震波斜入射时大坝地震动反应与垂直入射时明显不同,常规垂直入射的结果偏于不安全,因此在高面板坝地震反应分析和抗震设计中应考虑地震波斜入射的影响。     

台风激发的第二类地脉动特征及激发模式分析

0.003~1 Hz频段的地脉动主要来源于海浪运动与固体地球的耦合作用,台风引起的强烈海浪运动往往可使地脉动能量显著增强.由于涉及大气-海洋-固体地球三个圈层之间的复杂动量传递与耦合过程,迄今为止,关于台风激发地脉动的具体源区位置及激发机制尚存在争议.本文选取日本、中国东南沿海及台湾地区的地震台站波形连续记录,研究了2008年台风"森拉克"和"黑格比"激发地脉动的时频特征,开展相应数值模拟,并与观测数据进行了对比分析研究.结果表明台风激发第二类地脉动存在两种主要模式:(1)近岸源区激发,即台风引起波浪入射至海岸反射并与后续来波相互作用形成驻波作用于海底而激发;(2)台风中心附近源区激发,即台风中心移动过程中不同时期激发的同频率波浪相向传播、相互作用产生驻波作用于海底而激发,源区位置主要集中于台风中心左后方.此外,结合波浪再分析数据、台风风场特征,我们进一步对第二类地脉动激发过程中的影响因素进行了分析,发现:第一种模式激发的地脉动与近岸源区波浪场强度、观测点至源区距离及台风中心至海岸线距离等因素相关;而第二种模式激发的地脉动则主要受台风中心附近波浪场的频率成分与传播方向影响.     

地震波模拟中多轴复频移近似完全匹配层的吸收效果及稳定性

在进行地震波模拟计算的过程中用有限的计算区域模拟地下无限空间,需要进行边界截断.为了在边界处不产生虚假反射影响模拟结果,需要引入吸收边界条件.本文采用的近似完全匹配层是一种新型非分裂完全匹配层,计算效率较高.同时相比于其他非分裂完全匹配层,其还具有不改变方程的形式、易于实现等优势.但是当入射波角度较大,边界吸收效果变弱,且残留在边界中的能量使近似完全匹配层变得极其不稳定.多轴复频移近似完全匹配层的提出就是为了改善对大角度入射波的吸收并且提高边界的稳定性.通过实验模拟和矩阵特征值灵敏度来研究多轴复频移近似完全匹配层的吸收效果及稳定性.结果表明该方法不仅能够吸收掠入波,而且对常规入射波的吸收也得到提升,同时拥有更好的稳定性.     

基于动力有限差分法的典型海域自由场二维地震反应分析

以我国近海海域工程场地为研究对象,充分考虑上覆海水的自重影响,构建典型的饱和海底自由场计算模型,运用动力有限差分法开展二维地震反应分析,探讨以不同幅值的SV波、P波作为基底输入条件下上覆海水厚度对海底地震动峰值和反应谱的影响,总结上覆有水和无水场地的地震动结果差异,并分析差异产生的原因。结果表明:当基底输入SV波时,上覆有水场地海床表面峰值加速度小于上覆无水场地地表峰值加速度,海水层厚度对峰值加速度的影响可以忽略;当基底输入P波时,上覆有水场地海床表面峰值加速度大于上覆无水场地地表峰值加速度,且随着海水层厚度的增大,海床表面峰值加速度逐渐减小。      

RESEARCH ON TIME-DISTANCE INTERPRETATION OF MULTIPLES AND GHOSTS BASED ON MARINE SEISMIC DATA

Currently, the study on the active fault in the land areas is relatively mature, while there is still lack of detection and research on active faults in the sea areas. Marine exploration, which is different from land areas, has a prominent problem due to the existence of strong reflecting interfaces such as water surface and seafloor in the sea, thus the recording is often accompanied by interference of multiples on seafloor reflections. In addition, because of the characteristics of marine seismic exploration, the source exciting in the water and the geophone receiving in the water, ghost wave usually can be recorded simultaneously during the reflected wave propagation. This phenomenon makes it difficult to distinguish the effective waves and the noise, and has always plagued the data and seriously affects the quality of records. In the offshore and other regions of complex structures, such as inclined interfaces, it is difficult to eliminate the interference of multiples accurately by traditional multiples suppression methods, which are based on the horizontal interface assumption. This paper combines the sea area seismic data and its acquisition method, uses simplified model to simulate the multiples based on the time-distance analysis of multiples and their ghost wave in inclined interface. The time-distance characteristics of the multiples and their ghost waves from different interfaces(including the inclined interface)are obtained, and they are consistent with the actual records. The multiples time-distance simulation can help to distinguish the causes of reflected waves, summarize the multiple-wave time-distance characteristics from different interfaces(including inclined interfaces), and analyze the relationship between the characteristics of multiple waves and primary waves. In particular, this simulation has a significant effect on characterizing the internal multiples that are difficult to identify due to inconspicuous periodicity and the multiples of the inclined interface which present the phenomenon that the vertex of the time-distance curve is shifted. On this basis, relying on the time-distance analysis of ghost wave, we analyze the travel time difference characteristics between reflected waves and their accompanying ghost waves. The differences of the travel time characteristics of different orders ghost wave and reflected wave are summarized and the symmetry of the travel time difference between inclined interface and horizontal interface of ghost waves and reflected waves is analyzed. We simulate the distraction of the ghost wave event with the event of the reflected wave and analyze the influence of the ghost wave on the sea area seismic records. These results can improve the practical interpretation of seismic data. At last, the time-distance information is used to synthesize sea area seismic records, which can help us carry out the effective data processing and understand the characteristics of the time-distance and velocity of multiples in different interfaces and the layer artifact caused by multiples. This study combines the time-distance simulation of multiples and their ghost wave with conventional seismic data processing to analyze the pre-stack and post-stack features of multiple waves and their ghost waves in the seismic records of the sea area. The results of this study are conducive to the effective identification of multiples in seismic records in the sea, provide a theoretical basis for multi-wave suppression and prediction, and may facilitate the future study of sea-area seismic activity detection.     

Modeling the dynamic process of tsunami earthquake by liquid-solid coupling model

Tsunami induced by earthquake is an interaction problem between liquid and solid.Shallow-water wave equation is often used to modeling the tsunami,and the boundary or initial condition of the problem is determined by the displacement or velocity field from the earthquake under sea floor,usually no interaction between them is consid-ered in pure liquid model.In this study,the potential flow theory and the finite element method with the interaction between liquid and solid are employed to model the dynamic processes of the earthquake and tsunami.For model-ing the earthquake,firstly the initial stress field to generate the earthquake is set up,and then the occurrence of the earthquake is simulated by suddenly reducing the elastic material parameters inside the earthquake fault.It is dif-ferent from seismic dislocation theory in which the relative slip on the fault is specified in advance.The modeling results reveal that P,SP and the surface wave can be found at the sea surface besides the tsunami wave.The surface wave arrives at the distance of 600 km from the epicenter earlier than the tsunami 48 minutes,and its maximum amplitude is 0.55 m,which is 2 times as large as that of the sea floor.Tsunami warning information can be taken from the surface wave on the sea surface,which is much earlier than that obtained from the seismograph stations on land.The tsunami speed on the open sea with 3 km depth is 175.8 m/s,which is a little greater than that pre-dicted by long wave theory,(gh)1/2=171.5 m,and its wavelength and amplitude in average are 32 km and 2 m,respectively.After the tsunami propagates to the continental shelf,its speed and wavelength is reduced,but its amplitude become greater,especially,it can elevate up to 10 m and run 55 m forward in vertical and horizontal directions at sea shore,respectively.The maximum vertical accelerations at the epicenter on the sea surface and on the earthquake fault are 5.9 m/s2 and 16.5 m/s2,respectively,the later is 2.8 times the former,and therefore,sea water is a good shock     

Influences of type of wave and angle of incidence on seismic bending moments in pile foundations

When analysing the seismic response of pile groups, a vertically‐incident wavefield is usually employed even though it does not necessarily correspond to the worst case scenario. This work aims to study the influences of both the type of seismic body wave and its angle of incidence on the dynamic response of pile foundations. To this end, the formulation of SV, SH and P obliquely‐incident waves is presented and implemented in a frequency‐domain boundary element‐finite element code for the dynamic analysis of pile foundations and piled structures. Results are presented in terms of bending moments at cap level of single piles and 3 × 3 pile groups, both in frequency and in time domains. It is found that, in general, the vertical incidence is not the most unfavourable situation. In particular, obliquely‐incident SV waves with angles of incidence smaller than the critical one, a situation in which the mechanism of propagation of the waves in the soil changes and surface waves appear, yield bending moments much larger than those obtained for vertically‐incident wavefields. It is also shown that the influence of pile‐to‐pile interaction on the kinematic bending moments becomes significant for non‐vertical incidence, especially for P and SV waves. Copyright © 2013 John Wiley & Sons, Ltd.     

地震体波斜入射下的断层台阶地震反应分析

应用显式有限元法对倾斜断层台阶在地震体波斜入射下的地震反应进行了数值模拟。模拟结果表明，P波入射时的竖向运动分量和SV波入射时的水平向运动分量在台阶上角点处存在极大值、在台阶下角点处存在极小值，同时在上下台面还出现由角点激发、向台阶外传播的转换体波与面波。这一结果说明，在断层台阶场地建设重大工程时应离开上角点一定距离，并考虑上下水平表面地面运动的差异性。     

P波入射下断层参数对地表地震动的影响

基于应用透射人工边界条件的显式有限元法计算断层破碎带宽度及力学参数变化、地震动入射角变化时二维断层场地模型P波入射下地表地震动场的分布。结果表明:(1)低速度破碎带的存在导致整个场地都有P波转换为SV波的分量,且在断层破碎带的区域出现断层陷波;(2)低速度破碎带的存在使输入场地恒定的能量向破碎带集聚放大,随着破碎带宽度增大或其介质波速降低集聚放大效应增大;(3)场地放大效应是频率相关的,宽度较宽或介质波速较低的断层破碎带对输入地震动中较低的频率成份放大显著;(4)竖向断层破碎带能阻隔斜入射地震P波,随着入射角增加隔震效应更显著。     

New approach to analyzing soil-building systems

A new method of analyzing seismic response of soil-building systems is introduced. The method is based on the discrete-time formulation of wave propagation in layered media for vertically propagating plane shear waves. Buildings are modeled as an extension of the layered soil media by assuming that each story in the building is another layer. The seismic response is expressed in terms of wave travel times between the layers, and the wave reflection and transmission coefficients at layer interfaces. The calculation of the response is reduced to a pair of simple finite-difference equations for each layer, which are solved recursively starting from the bedrock. Compared with commonly used vibration formulation, the wave propagation formulation provides several advantages, including the ability to incorporate soil layers, simplicity of the calculations, improved accuracy in modeling the mass and damping, and better tools for system identification and damage detection.     

频率-空间域CPML和特征分析法组合边界研究

在地震波场数值模拟过程中,边界反射是影响其模拟结果的一个重要因素。实际地下介质具有各向异性特征,传统的完全匹配层边界（PML）对于小入射角地震波具有良好效果,但该方法并不能有效地吸收低频波和大角度入射波。针对VTI介质边界反射的问题,本文提出在频率-空间域有限差分法数值模拟中采用卷积完全匹配层（CPML）和特征分析法的组合边界条件,并对该组合边界条件进行数值模拟实验和边界反射吸收效果分析,验证所提方法是一种可靠的人工吸收边界条件,能够有效地压制波场模拟过程中产生的边界反射。      

Sea surface shape derivation above the seismic streamer

The rough sea surface causes perturbations in the seismic data that can be significant for time‐lapse studies. The perturbations arise because the reflection response of the non‐flat sea perturbs the seismic wavelet. In order to remove these perturbations from the received seismic data, special deconvolution methods can be used, but these methods require, as input, the time varying wave elevation above each hydrophone in the streamer. In addition, the vertical displacement of the streamer itself must also be known at the position of each hydrophone and at all times. This information is not available in conventional seismic acquisition. However, it can be obtained from the hydrophone measurements provided that the hydrophones are recorded individually (not grouped), that the recording bandwidth is extended down to 0.05 Hz and that data are recorded without gaps between the shot records. The sea surface elevation, and also the wave‐induced vertical displacement of the streamer, can be determined from the time‐varying pressure that the sea waves cause in the hydrophone measurements. When this was done experimentally, using a single sensor seismic streamer without a conventional low cut filter, the wave induced pressure variations were easily detected. The inversion of these experimental data gives results for the sea surface elevation that are consistent with the weather and sea state at the time of acquisition. A high tension approximation allows a simplified solution of the equations that does not demand a knowledge of the streamer tension. However, best results at the tail end of the streamer are obtained using the general equation.     

地震波斜入射下层状场地中地下综合管廊地震响应分析

利用黏弹性人工边界和等效地震荷载时域波动输入方法,结合土层和半空间的精确动力刚度矩阵,实现了地震波斜入射下层状场地地下综合管廊地震反应分析,建立了不同场地条件下地下综合管廊分析模型。计算结果表明:地震波倾斜入射情况下,综合管廊结构地震响应与垂直入射时具有显著差异,一般SV波以30°临界角附近入射时结构地震反应最为剧烈;地下综合管廊动应力集中主要分布在管廊角部、中柱上下端;成层土波速结构变化对地下综合管廊地震反应亦具有显著影响。总体上看:当穿越软夹层时管廊结构地震反应更为剧烈,且覆盖层越厚,管廊结构内力幅值越大。因此地下综合管廊结构抗震设计宜考虑地震波倾斜入射及场地土层性质的影响。