Analytical expressions for Fresnel volumes and interface fresnel zones of seismic body waves. Part 2: Transmitted and converted waves. Head waves

Summary Fresnel volumes and interface Fresnel zones of transmitted and head waves are studied. The relation derived for transmitted waves may also be used for converted reflected waves. Considerable attention is devoted to the penetration of Fresnel volumes across structural interfaces, particularly for head waves.     

Analytical expression for fresnel volumes and interface fresnel zones of seismic body waves. Part 1: Direct and unconverted reflected waves

Summary Fresnel volumes, plane-sectional Fresnel zones and interface Fresnel zones of direct and unconverted reflected waves are studied. Exact analytical expressions for various parameters of Fresnel volumes and Fresnel zones are derived and discussed. Among others, these expressions are related to semi-axes of Fresnel zones, to overshooting and penetration distances, and to the off-ray shifts of Fresnel zones.     

菲涅耳体地震层析成像分辨率研究

层析成像分辨率的研究,不仅可以帮助分析层析方法的反演能力,评价层析反演的效果,还可以帮助指导层析参数设置,优化观测系统设计等.本文对比研究了前人提出的两种菲涅耳体层析成像分辨率的计算方法,并针对其存在的问题进行了优化.文中通过对二维理论模型的定量计算,总结了菲涅耳体地震层析成像分辨率的一些规律,并将其与射线层析的分辨率进行了对比.文中最后给出分辨率研究对观测系统与层析成像的指导作用,尤其是对层析中平滑窗设置、变网格层析中模型剖分原则给出了具体的建议.     

Numerical modeling of global seismic phases and its application in seismic phase identification

Earthquake data include informative seismic phases that require identification for imaging the Earth's structural interior. In order to identify the phases, we created a numerical method to calculate the traveltimes and raypaths by a shooting technique based upon the IASP91 Earth model, and it can calculate the traveltimes and raypaths for not only the seismic phases in the traditional traveltime tables such as IASP91, AK135, but also some phases such as pPcP, pPKIKP, and PPPPP. It is not necessary for this method to mesh the Earth model, and the results from the numerical modeling and its application show that the absolute differences between the calculated and theoretical traveltimes from the ISAP91 tables are less than 0.1 s. Thus, it is simple in manipulation and fast in computation, and can provide a reliable theoretical prediction for the identification of a seismic phase within the acquired earthquake data.     

Comparison of techniques that use the single scattering model to compute the quality factorQ from coda waves

Using simulated coda waves, the resolution of the single-scattering model to extract codaQ (Q _c ) and its power law frequency dependence was tested. The back-scattering model ofAki andChouet (1975) and the single isotropic-scattering model ofSato (1977) were examined. The results indicate that: (1) The inputQ _c models are reasonably well approximated by the two methods; (2) almost equalQ _c values are recovered when the techniques sample the same coda windows; (3) lowQ _c models are well estimated in the frequency domain from the early and late part of the coda; and (4) models with highQ _c values are more accurately extracted from late code measurements.     

The effect of horizontal resolution on gravity waves simulated by the GFDL “SKYHI” general circulation model

To examine the effects of horizontal resolution on internal gravity waves simulated by the 40-level GFDL SKYHI general circulation model, a comparison is made between the 3° and 1° resolution models during late December. The stratospheric and mesospheric zonal flows in the winter and summer extratropical regions of the 1° model are much weaker and more realistic than the corresponding zonal flows of the 3° model. The weaker flows are consistent with the stronger Eliassen-Palm flux divergence (EPFD).The increase in the magnitude of the EPFD in the winter and summer extratropical mesospheres is due mostly to the increase in the gravity wave vertical momentum flux convergence (VMFC). In the summer extratropical mesosphere, the increase in the resolvable horizontal wavenumbers accounts for most of the increase in the gravity wave VMFC. In the winter extratropical mesosphere, the increase of VMFC associated with large-scale eastward moving components also accounts for part of the increase in the gravity wave VMFC.The gravity waves in the summer and winter mesosphere of the 1° model are associated with a broader frequency-spectral distribution, resulting in a more sporadic time-distribution of their VMFC. This broadening is due not only to the increase in resolvable horizontal wavenumbers but also occurs in the large-scale components owing to wave-wave interactions. It was found that the phase velocity and frequency of resolvable small-scale gravity waves are severely underestimated by finite difference approximations.     

Two-dimensional scattering of SH waves due to a discontinuity in bedrock

The direct boundary integral equation technique is applied to determine the impact on surface amplification caused by an inhomogeneity in a bedrock half-space. The particular soil-rock configuration studied is one in which a soil layer rests on a rock half-space which includes a rock inclusion. The particular rock inclusion considered for this study is a semi-infinite rock layer with its upper boundary bordering the soil layer. Materials are considered viscoelastic except for the section of the rock half-space below the level of the rock inclusion which is considered elastic. A parametric study is performed to determine controlling factors for surface displacement due to a vertically incident shear wave. The study includes varying the stiffness and the thickness of the inclusion for a range of frequencies. Solutions from a one-dimensional analysis are compared with the results of a two-dimensional analysis in order to establish limits inside of which a two-dimensional analysis is required. Copyright © 1999 John Wiley & Sons, Ltd.     

Two-dimensional scattering of in-plane body waves due to a discontinuity in bedrock

The direct boundary integral equation technique is used to study in-plane surface amplification of in-plane seismic body waves for the case of an inhomogeneity in a bedrock half-space. In the studied soil configuration, a soil layer rests on a rock half-space which includes a rock inclusion. The rock inclusion considered is a semi-infinite horizontal rock layer in which its upper boundary borders the soil layer. Materials in the soil–rock configuration are considered viscoelastic except for the section of the rock half-space below the level of the rock inclusion which is considered elastic. A parametric study is performed to determine controlling factors for surface displacement due to in-plane body waves. The study investigates varying the stiffness and the thickness of the rock inclusion for a range of frequencies and wave incidence angles. Anti-plane waves for this type of soil-rock configuration have been addressed in a previous article by Heymsfield (Earthquake Engng. Struct. Dyn. 28 : 841–855 (1999)). Copyright © 1999 John Wiley & Sons, Ltd.     

与规范反应谱相对应的Clough-Penzien模型参数研究

根据我国现行抗震规范（GB50011—2001）的反应谱曲线，对Clough—Penzien模型的参数取值进行了具体研究。采用时间包络函数考虑地震动的非平稳性，根据加速度峰值等效原则确定了谱强度因子S0的表达式，表明谱强因子不仅与地面加速度特性、场地类别有关，而且与结构的动力特性（阻尼比、自振周期）有关。最后对谱强度因子计算做了简化处理，为随机抗震计算分析提供了参考依据。     

Modified image algorithm to simulate seismic channel waves in 3D tunnel model with rugged free surfaces

The existence of rugged free‐surface three‐dimensional tunnel conditions in the coal seams, caused either by geological or mining processes, will inevitably influence wave propagation characteristics when the seismic waves go through the coal mines. Thus, a modified image algorithm has been developed to account for seismic channel waves propagating through this complicated topography with irregular free surfaces. Moreover, the seismic channel waves commonly exhibit damped and dispersive signatures, which is not only because of their own unique sandwich geometry of rock–coal–rock but also because of the viscoelastic behavior of coal. Considering the complexity of programming in three‐dimensional tunnel models with rugged free surfaces, an optimized vacuum grid search algorithm, enabling to model highly irregular topography and to compute efficiently, is also proposed when using high‐order staggered finite‐difference scheme to simulate seismic channel wave propagations in viscoelastic media. The numerical simulations are implemented to investigate the accuracy and stability of the method and the impact of coal's viscoelastic behavior on seismic channel wave propagation characteristics. The results indicate that the automatic vacuum grid search algorithm can be easily merged into high‐order staggered finite‐difference scheme, which can efficiently be applied to calculate three‐dimensional tunnel models with rugged free surfaces in the viscoelastic media. The simulation also suggests that the occurrence of a three‐dimensional tunnel with free surfaces has a remarkable influence on the seismic channel wave propagation characteristics and elastic energy distribution.     

地震层析成像方法综述

本文从地震层析成像的模型参数化、正演、反演、解评估四个方面回顾了地震层析成像方法近期发展历程以及各时段地震层析成像方法的优、缺点。     

Resolution assessment of the three-dimensional acquisition geometry of the Anhui experiment in "Geoscience Yangtze Plan"

According to the actual observation conditions of the Yangtze River valley from Anqing city to Maanshan city, we designed the 3D acquisition geometry, and applied the multi-scale checkerboard semblance analysis to assess the preliminary resolution of the designed observation. The checkerboard semblance tests use the refraction and reflection travel-time simultaneous inversion algorithm to quantitatively provide both resolution assessment of velocity structure and Moho topography. The multi-scale checkerboard semblance recovery results show that while the checkerboard semblance threshold value is 0.5, the preliminary resolution of the designed acquisition geometry is better than 10 km in the upper crust (the depth is less than 10 km), around 15 km in the mid-crust (the depth is 10?25 km), and better than 20 km in the lower crust (the depth is 25?33 km). The preliminary tomographic resolution for the Moho topography is about 20 km in the ray-path coverage area beneath the acquisition geometry. While the checkerboard semblance threshold value is 0.75, the preliminary resolution is 20 km in the upper crust, around 20?25 km in the mid-crust and 25 km in the lower crust. And the preliminary tomographic resolution for the Moho topography is better than 30 km in the ray-path coverage area beneath the acquisition geometry. These non-linear checkerboard tests reveal that the designed acquisition geometry is suitable to image the crustal velocity structure of the Yangtze River valley in the Anhui province.     

On steady and travelling waves over bottom topography in the model of homogeneous flow in a beta-plane channel

Abstract

A spectral low-order model is proposed in order to investigate some effects of bottom corrugation on the dynamics of forced and free Rossby waves. The analysis of the interaction between the waves and the topographic modes in the linear version of the model shows that the natural frequencies lie between the corresponding Rossby wave frequencies for a flat bottom and those applying in the “topographic limit” when the beta-effect is zero. There is a possibility of standing or eastward-travelling free waves when the integrated topograhic effect exceeds the planetary beta-effect.

The nonlinear interactions between forced waves in the presence of topography and the beta-effect give rise to a steady dynamical mode correlated to the topographic mode. The periodic solution that includes this steady wave is stable when the forcing field moves to the West with relatively large phase speed. The energy of this solution may be transferred to the steady zonal shear flow if the spatial scale of this zonal mode exceeds the scale of the directly forced large-scale dynamical mode.     

Simulation of the stratospheric gravity waves generated by the Typhoon Matsa in 2005

The generation of stratospheric gravity waves(GWs) due to typhoon is simulated by using a meso-scale model(WRF) with a typhoon case,the Matsa in 2005.An 8-day model run that covers the major stages of the Matsa’s development reproduces the key features of the typhoon.For example,good agreements in the typhoon’s track,the intensity,and the spiral clouds,as well as mean state of stratosphere,are seen between the simulation results and the observation.Simulation results clearly show that with typhoon propagates northwestward,pronounced stratospheric GWs are generated continuously in the vicinity of Matsa.The GWs exhibit the typical curve-like wave fronts away from the Typhoon Matsa,and propagate preferentially in the upstream of the background winds.These characteristics reflect that the stratospheric GWs are closely associated with the typhoon,and thus the GWs are referred to as Tropical Cyclone related Gravity Waves(TC-GWs).The results also show that these waves should have a rather large horizontal scale so that the outmost wave fronts can be seen at the distance of ～1000 km to the typhoon center in the horizontal plane of 20 km.This is consistent with the phenomenon of stratospheric TC-GWs with ～1000 km horizontal scale disclosed by the previous observational analysis results.     

Revision of seismic design codes corresponding to building damagesin the “5.12” Wenchuan earthquake

A large number of buildings were seriously damaged or collapsed in the 5.12 Wenchuan earthquake. Based on field surveys and studies of damage to different types of buildings, seismic design codes have been updated. This paper briefly summarizes some of the major revisions that have been incorporated into the Standard for classification of seismic protection of building constructions GB50223-2008 and Code for Seismic Design of Buildings GB50011-2001. The definition of seismic fortification class for buildings has been revisited, and as a result, the seismic classifications for schools, hospitals and other buildings that hold large populations such as evacuation shelters and information centers have been upgraded in the GB50223-2008 Code. The main aspects of the revised GB50011-2001 code include: (a) modification of the seismic intensity specified for the Provinces of Sichuan, Shanxi and Gansu; (b) basic conceptual design for retaining walls and building foundations in mountainous areas; (c) regularity of building configuration; (d) integration of masonry structures and precast RC floors; (e) requirements for calculating and detailing stair shafts; and (f) limiting the use of single-bay RC frame structures. Some significant examples of damage in the epicenter areas are provided as a reference in the discussion on the consequences of collapse, the importance of duplicate structural systems, and the integration of RC and masonry structures.     

井间高分辨率纵横波层析成像研究井间油藏

采用视速度-偏振井间波场分离法分离得到纵波和横波波场,使用跨孔地震走时层析成像级联方法计算井间纵波和横波速度结构以及泊松比空间分布.综合纵横波速度、泊松比和测井参数,分析两井间岩性和储层特性.数值模型试验表明该方法垂直分辨率达到主频波长的1/4,并能准确分辨倾斜断层和垂直断层.处理相同纵波资料结果表明本文方法分辨率远远高于以往的井间走时层析成像方法,比井间声波层析成像分辨率更高.使用该方法处理垦71区井间资料,得到同台同源纵波和横波速度结构和泊松比分布.泊松比的低/高变化与测井参数指示的砂/泥岩层基本吻合.根据泊松比分布,参考岩石地震物理学测试数据,可以区分井间介质中的泥岩层、砂岩层,以及砂岩层中的饱水和饱油区,确定储层的连通性,圈定井间的剩余油藏.     

Numerical model studies on the generation and dissipation of the diurnal coastal-trapped waves over the Sakhalin shelf in the Sea of Okhotsk

To clarify the generation and dissipation mechanisms of diurnal coastal-trapped waves (CTWs) over the Sakhalin shelf, a series of numerical experiments were conducted using a three-dimensional tidal model of the Okhotsk Sea with density stratification. The tidal model used has good reproduction owing to the careful fitting to the recent observations. The numerical experiments suggested that diurnal CTWs are primarily (~60%) generated by the conversion of tidal energy at the northern corner of the Sakhalin shelf, and further amplified by vorticity generation due to the water column oscillation from Sakhalin Bay and the influence of Kashevarov Bank. From the observations, it was found that diurnal CTWs are effectively dissipated by the strong spin-down due to bottom friction. The conventional turbulent closure model cannot reproduce the observed damping of diurnal CTWs, which raises a caution in modeling the tidal fields in high-latitude regions where diurnal CTWs exist. To resolve this underestimation of the damping, the vertical eddy viscosity was parameterized using its dependence on the observed major axis length of the diurnal tidal current ellipses, which improves the model reproduction on the damping of diurnal CTWs. The model also suggests that the spin-down effects due to friction associated with the sea-ice cover play an important role in the tidal current reduction in the region where diurnal CTWs exist, as the observations suggested.     

地震的"序列归属"问题与ETAS模型——以唐山序列为例

历史上发生过强震地区的余震活动可能持续较长时间,而余震序列在何时可被看作正常的"背景地震活动",即"序列归属"问题在地球动力学和地震物理中有重要意义.时-空"传染型余震序列"(ETAS)模型可分离"背景"地震和"丛集"地震,并用概率形式表示作为相应事件的可能性,为考察此问题提供了可能.本文以1976年唐山M_S7.8地震序列为例,对唐山地区1970年以来的M_L4.0以上地震进行了时-空ETAS模型拟合,并以2010年以来发生的3次M_S4.0以上地震为例讨论了它们的"序列归属"问题.研究结果显示,3次M_S4.0以上地震的背景地震概率分别为0.72、0.88和0.76,表明它们作为1976年唐山M_S7.8的余震的可能性较低,更可能为背景地震.     

On leakage of energy from turbulence to internal waves in the oceanic mixed layer

In this paper, we address the question of energy leakage from turbulence to internal waves (IWs) in the oceanic mixed layer (OML). If this leakage is substantial, then not only does this have profound implications as far as the dynamics of the OML is concerned, but it also means that the equation for the turbulence kinetic energy (TKE) used in OML models must include an appropriate sink term, and traditional models must be modified accordingly. Through comparison with the experimental data on grid-generated turbulence in a stably stratified fluid, we show that a conventional two-equation turbulence model without any IW sink term can explain these observations quite well, provided that the fluctuating motions that persist long after the decay of grid-generated turbulence are interpreted as being due to IW motions generated by the initial passage of the grid through the stably stratified fluid and not during turbulence decay. We conclude that there is no need to postulate an IW sink term in the TKE equation, and conventional models suffice to model mixing in the OML.     

测井约束反演技术在煤田三维地震岩性勘探中的应用

利用地震技术解决岩性问题在煤田勘探中一直是近年来人们研究的一个热门课题,与其他地震反演技术相比,测井约束反演技术不仅可提高原始地震资料的垂向分辨率,而且将地震波转换为层的信息,因此正在岩性勘探中得到广泛应用．本文介绍了测井约束反演技术的原理及影响因素,并以开滦集团荆各庄矿某采区岩性勘探为例,说明了测井约束反演技术在煤田三维地震岩性勘探中的具体方法和应用效果。