Analytical study of ground motion caused by seismic wave propagation across faulted rock masses

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.     

基于深度域地震成像的中沙海槽盆地东北部结构构造研究

文章基于叠前深度偏移地震成像分析, 结合前人重磁反演等地球物理资料成果, 对中沙海槽盆地东北部结构构造进行探索研究。研究表明: 中沙海槽盆地东北部发育新生代地层, 厚度在1500~2500m之间, 地层层速度为1500~ 3500m·s^-1, 不存在中生界沉积层, 盆地基底为海山或前寒武系基岩及岩浆岩隆起。中沙海槽盆地新生代陆源海相烃源岩丰富, 盆地新生界地层厚度大, 盆地凹陷、断裂发育, 具有一定的油气勘探潜力。     

Numerical manifold method for dynamic nonlinear analysis of saturated porous media

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.     

Synthetic SH seismograms in a laterally varying medium by the discrete wavenumber method

Summary. We present a new method to calculate the SH wavefield produced by a seismic source in a half-space with an irregular buried interface. The diffracting interface is represented by a distribution of body forces. The Green's functions needed to solve the boundary conditions are evaluated using the discrete wavenumber method. Our approach relies on the introduction of a periodicity in the source-medium configuration and on the discretization of the interface at regular spacing. The technique developed is applicable to boundaries of arbitrary shapes and is valid at all frequencies. Some examples of calculation in simple configurations are presented showing the capabilities of the method.     

中国大陆科学钻探先导孔零偏VSP资料解释

用六级三分量检波器在中国大陆科学钻探先导孔中实施了零偏VSP测量.数据处理结果表明,中国大陆科学钻探孔区超高压变质岩石的地震波速度主要介于4500～7000m/s之间,显著高于一般的沉积岩地区,而且随深度变化不明显.声波测井速度系统地稍低于VSP层速度,可能是由于井壁处岩石的完整性受到破坏而造成的.地震波速度与岩石密度和岩性存在明确的对应关系,榴辉岩的密度和地震波速度均显著高于片麻岩类岩石;由榴辉岩退变生成的斜长角闪岩类岩石,其密度和地震波速度均呈现出较大的变化,主要与其退变质程度有关;超基性岩中的裂隙系统导致其密度和地震波速度大幅度下降.由于榴辉岩与其他岩石类型之间存在较大的波阻抗差异,因此用零偏VSP资料标定该区地震波的地质层位是有效的.关于地震波反射的原因,通过综合研究地震波(包括反射纵波、上行转换横波、井筒波)的特征、岩石速度和密度分布以及井径变化,认为主要是岩性分界面、韧性剪切带和断裂(带),但还有一些因素尚待进一步研究.     

Restricted Three-Body Problem: An Approximation of its General Solution – Part One – the Manifold of Symmetric Periodic Solutions

This paper gives the results of a programme attempting to exploit ‘la seule bréche’ (Poincaré, 1892, p. 82) of non-integrable systems, namely to develop an approximate general solution for the three out of its four component-solutions of the planar restricted three-body problem. This is accomplished by computing a large number of families of ‘solutions précieuses’ (periodic solutions) covering densely the space of initial conditions of this problem. More specifically, we calculated numerically and only for μ = 0.4, all families of symmetric periodic solutions (1st component of the general solution) existing in the domain D:(x ₀ ∊ [−2,2],C ∊ [−2,5]) of the (x ₀, C) space and consisting of symmetric solutions re-entering after 1 up to 50 revolutions (see graph in Fig. 4). Then we tested the parts of the domain D that is void of such families and established that they belong to the category of escape motions (2nd component of the general solution). The approximation of the 3rd component (asymmetric solutions) we shall present in a future publication. The 4th component of the general solution of the problem, namely the one consisting of the bounded non-periodic solutions, is considered as approximated by those of the 1st or the 2nd component on account of the `Last Geometric Theorem of Poincaré' (Birkhoff, 1913). The results obtained provoked interest to repeat the same work inside the larger closed domain D:(x ₀ ∊ [−6,2], C ∊ [−5,5]) and the results are presented in Fig. 15. A test run of the programme developed led to reproduction of the results presented by Hénon (1965) with better accuracy and many additional families not included in the sited paper. Pointer directions construed from the main body of results led to the definition of useful concepts of the basic family of order n, n = 1, 2,… and the completeness criterion of the solution inside a compact sub-domain of the (x ₀, C) space. The same results inspired the ‘partition theorem’, which conjectures the possibility of partitioning an initial conditions domain D into a finite set of sub-domains D i that fulfill the completeness criterion and allow complete approximation of the general solution of this problem by computing a relatively small number of family curves. The numerical results of this project include a large number of families that were computed in detail covering their natural termination, the morphology, and stability of their member solutions. Zooming into sub-domains of D permitted clear presentation of the families of symmetric solutions contained in them. Such zooming was made for various values of the parameter N, which defines the re-entrance revolutions number, which was selected to be from 50 to 500. The areas generating escape solutions have being investigated. In Appendix A we present families of symmetric solutions terminating at asymptotic solutions, and in Appendix B the morphology of large period symmetric solutions though examples of orbits that re-enter after from 8 to 500 revolutions. The paper concludes that approximations of the general solution of the planar restricted problem is possible and presents such approximations, only for some sub-domains that fulfill the completeness criterion, on the basis of sufficiently large number of families.     

全球Ms≥7￣3／4级地震的定向迁移

对全球尺度的６条大地震带内１９００～１９９０年中１８４次Ｍｓ≥７３／４级地震进行了沿地震带方向定向迁移的分析，获得了全球统一的地震定向迁移规律，总体是由西向东，迁移速度由７００ｋｍ／ａ变为１５０ｋｍ／ａ，此现象可以作多种暂态地球动力作用过程的推论，如以大西洋脊间歇式张裂引起上地幔软流物质自西向东运动，呈现纵波式的振荡传播；也可解释为非洲板块、阿拉伯板块和印度板块自西南向东北对欧亚地震带依次的推压引起向东的应变波的传播；太平洋脊两侧洋底板块向西北和东北两侧的斜向推压，可能是造成两侧地震带地震向北迁移的触发源