Directional dependence of P- and S-wave propagation and polarization in foliated rocks from the Kola superdeep well: Evidence from laboratory measurements and calculations based on TOF neutron diffraction

We have measured P- and S-wave velocities on two amphibolite and two gneiss samples from the Kola superdeep borehole as a function of pressure (up to 600 MPa) and temperature (up to 600 °C). The velocity measurements include compressional (Vp) and shear wave velocities (Vs1, Vs2) propagating in three orthogonal directions which were in general not parallel to inherent rock symmetry axes or planes. The measurements are accompanied by 3D-velocities calculations based on lattice preferred orientation (LPO) obtained by TOF (Time Of Flight) neutron diffraction analysis which allows the investigation of bulk volumes up to several cubic centimetres due to the high penetration depth of neutrons. The LPO-based numerical velocity calculations give important information on the different contribution of the various rock-forming minerals to bulk elastic anisotropy and on the relations of seismic anisotropy, shear wave splitting, and shear wave polarization to the structural reference frame (foliation and lineation). Comparison with measured velocities obtained for the three propagation directions that were not in accordance with the structural frame of the rocks (foliation and lineation) demonstrate that for shear waves propagating through anisotropic rocks the vibration directions are as important as the propagation directions. The study demonstrates that proper measurement of shear wave splitting by means of two orthogonal polarized sending and receiving shear wave transducers is only possible when their propagation and polarization directions are parallel and normal to foliation and lineation, respectively.     

Existence of transverse waves in anisotropic poroelastic media

Out of the four waves in an anisotropic poroelastic medium, two are termed as quasi-transverse waves. The prefix 'quasi' refers to their polarizations being nearly, but not exactly, perpendicular to direction of propagation. In this composite medium, unlike perfectly elastic medium, the propagation of a longitudinal wave along a phase direction may not be accompanied by transverse waves. The existence of a transverse wave in anisotropic poroelastic media is ensured by the two equations restricting the choice of elastic coefficients of porous aggregate as well as fluid–solid coupling. Necessary and sufficient conditions for the existence of transverse waves along the coordinate axes and in the coordinate planes for general anisotropy are discussed. The discussion is extended to the case of orthotropic materials and existence for few specific phase directions is also explored. The conditions for the transverse waves decided on the basis of their apparent polarizations, that is, particle motion being perpendicular to ray direction, are also discussed. For a particular numerical model, the existence of these apparent transverse waves is solved numerically for phase directions in coordinate planes. For general directions of phase propagation, the existence of these transverse waves is checked graphically for the chosen numerical model.     

A comparative study of the design spectra defined by Eurocode 8, UBC, IBC and Turkish Earthquake Code on R/C sample buildings

Earthquake codes have been revised and updated depending on the improvements in the representation of ground motions, soils and structures. These revisions have been more frequently seen in recent years. One of the key changes in earthquake codes has been performed on the design spectra. In this paper, the design spectra recommended by Turkish Earthquake Code and three other well known codes (Uniform Building Code, Eurocode 8, and International Building Code) are considered for comparison. The main purpose of this study is to investigate the differences caused by the use of different codes in the dynamic analysis and seismic verification of given types of buildings located at code defined different sites. The differences in expressions and some important points for elastic and inelastic spectra defined by the codes are briefly illustrated in tables and figures. Periods, base shears, lateral displacements and interstory drifts for the analyzed buildings located at code defined ground type are comparatively presented.     

Integration of reflection seismic and sediment grain-size data from Lake Khubsugul (Northern Mongolia): a reply to Prokopenko and Kendall

Prokopenko and Kendall (J Paleolimnol doi:, 2008) criticise the work presented in Fedotov et al. (J Paleolimnol 39:335–348, 2008), and instead propose an alternative interpretation for the grain-size evolution recorded in the KDP-01 core, retrieved from the central part of Lake Khubsugul. Their interpretation is based (i) on a seismic-stratigraphic re-interpretation of sparker seismic profile khub012 (which they copied from Fedotov et al. (EOS Trans 87:246–250, 2006)), (ii) on the presupposition that changes in lake level are the dominant control on facies distribution in Lake Khubsugul, and (iii) on the invalidation of our age-depth model. In this reply to their comment, we demonstrate that they interpreted seismic artefacts and geometries caused by changes in profile orientation as true stratigraphic features and that the lake-level reconstruction they derive from this interpretation is therefore incorrect. We also demonstrate that their grain-size predictions, which they consider to be predominantly driven by changes in lake level, are inconsistent with the measured sulphate concentration, which is a demonstrated proxy of lake level in Lake Khubsugul, and with the measured grain-size record. Finally, we point out that even if there would be a problem with the age-depth model, this problem would not affect the part of the sedimentary sequence discussed in Fedotov et al. (J Paleolimnol 39:335–348, 2008).     

基于Matlab的断裂带温泉水地球化学特征及地震活动性研究

基于Matlab开发出地下流体分析软件,对比分析了甘肃西秦岭北缘的武山、街子和清水温泉水化学组分特征;初步分析了温泉水水质类型、补给来源、水-岩平衡状态以及循环深度,并讨论了温泉水深循环对地震活动性的影响.研究表明温泉水均为大气降水成因,其水化特征受围岩的控制作用;清水温泉水的成熟度没有武山和街子温泉水高.地下水循环深度较浅的武山温泉和街子温泉水化类型属于Na-HCO3·SO4·Cl,所处断裂的地震活动频繁;地下水循环深度较深的清水温泉水化类型属于Na-SO4·Cl,所处断裂的地震活动性弱.     

二步法射线追踪

本文提出的二步法射线追踪,能够在复杂地质模型下有效、准确地确定射线在炮点的出射角,使得射线恰好到达检波点.射线到达地面的位置是射线在炮点的出射角的函数,称为出射函数.二步法射线追踪通过搜索出射函数的极值点来确定其单调变化的区间;然后通过比较检波点位置与出射函数极值的大小,确定与检波点相对应的出射角究竟落在出射函数的哪一个单调变化的区间里;最后,利用出射函数在局部区域里单调变化的性质,迅速、准确地确定与检波点相对应的出射角.出射函数的极值点,对于共炮集中每一个检波点和检波点中每一个不同方向的入射射线都是适用的,因此,二步法射线追踪只需要进行一次一维全局搜索（极值点搜索）.另外,本文采用B-样条的线性组合来表示速度函数,通过提高B-样条的阶数,可以改善速度函数的光滑程度,从而既满足射线追踪对于速度函数光滑性的要求,又加强射线在炮点的出射角对于射线的控制作用.     

Surface wave tomography for azimuthal anisotropy in a strongly reduced parameter space

Large scale seismic anisotropy in the Earth's mantle is likely dynamically supported by the mantle's deformation; therefore, tomographic imaging of 3-D anisotropic mantle seismic velocity structure is an important tool to understand the dynamics of the mantle. While many previous studies have focused on special cases of symmetry of the elastic properties, it would be desirable for evaluation of dynamic models to allow more general axis orientation. In this study, we derive 3-D finite-frequency surface wave sensitivity kernels based on the Born approximation using a general expression for a hexagonal medium with an arbitrarily oriented symmetry axis. This results in kernels for two isotropic elastic coefficients, three coefficients that define the strength of anisotropy, and two angles that define the symmetry axis. The particular parametrization is chosen to allow for a physically meaningful method for reducing the number of parameters considered in an inversion, while allowing for straightforward integration with existing approaches for modelling body wave splitting intensity measurements. Example kernels calculated with this method reveal physical interpretations of how surface waveforms are affected by 3-D velocity perturbations, while also demonstrating the non-linearity of the problem as a function of symmetry axis orientation. The expressions are numerically validated using the spectral element method. While challenges remain in determining the best inversion scheme to appropriately handle the non-linearity, the approach derived here has great promise in allowing large scale models with resolution of both the strength and orientation of anisotropy.     

用波速层析成像方法提高储层横向预测精度

以数值模拟方法为基础,考察了反射地震层析成像技术对油气储层波速横向检测的效果,然后用渤海油田的地震资料进行了实用试验．对数值模拟结果和实际数据的试验结果进行了分析与评价,认为油气勘探中的反射地震资料具有地震射线近似垂直分布的特点,适合用反射地震层析成像技术进行目标层横向非均匀性检测．     

Approximating constant‐Q seismic propagation in the time domain

In this study, we investigate the accuracy of approximating constant‐Q wave propagation by series of Zener or standard linear solid (SLS) mechanisms. Modelling in viscoacoustic and viscoelastic media is implemented in the time domain using the finite‐difference (FD) method. The accuracy of numerical solutions is evaluated by comparison with the analytical solution in homogeneous media. We found that the FD solutions using three SLS relaxation mechanisms as well as a single SLS mechanism, with properly chosen relaxation times, are quite accurate for both weak and strong attenuation. Although the RMS errors of FD simulations using a single relaxation mechanism increase with increasing offset, especially for strong attenuation (Q = 20), the results are still acceptable for practical applications. The synthetic data of the Marmousi‐II model further illustrate that the single SLS mechanism, to model constant Q, is efficient and sufficiently accurate. Moreover, it benefits from less computational costs in computer time and memory.     

影响b值计算误差的Monte Carlo实验研究

本文假设地震的发生与时间的关系遵循稳态泊松模型,假设地震震级M服从参数为1/β的指数分布,其中β=b/loge,在此基础上生成地震目录.用最小二乘法估算b值,分别从样本数量、震级间隔、震级误差三方面分析b值误差的大小.统计分析结果表明,b值的标准偏差随着样本数量的增多而减小;随着震级间隔的增大而减小;震级误差使b值的标准误差变大.本文可以为利用真实地震目录求b值时选择适当的方法提供参考.