基于横波分裂的青藏高原多圈层各向异性研究进展

有关青藏高原横波分裂的各向异性研究已经开展了近30年,在理论方法和实际应用方面取得了重要进展,并获取了大量的横波分裂测量结果,为认识青藏高原壳幔各向异性变形特征和动力学机制提供了重要的依据.本文首先介绍了地震各向异性的来源与应用,随后回顾了横波分裂分析方法的发展,简述了各种横波分裂方法的原理,最后通过总结近30年来青藏高原上地壳、整个地壳和上地幔横波分裂各向异性研究成果,系统分析了青藏高原壳幔各向异性变形特征.基于各横波分裂结果的对比分析来看,XKS波分裂测量结果最为稳定,近震直达S波分裂测量结果次之,而Pms波分裂测量结果相对离散,往往相同区域内不同的研究结果差异较大,主要原因可能是相比XKS波和近震直达S波,Pms波的信噪比较低,次要原因可能是各研究在方法和处理分析等方面的差异.     

各向异性介质控制束偏移及海上数据应用实例

随着勘探区域逐渐从陆地过渡到海洋,勘探目标逐渐趋于复杂化,高精度成像方法已经成为海洋油气勘探的瓶颈技术。高斯束偏移是一种灵活且高效的深度域偏移方法,对实际资料成像具有较好的适应性。该文发展了一种适应于海洋观测系统的高精度高斯束偏移方法,首先将海上接收的共偏移距地震记录进行加窗局部倾斜叠加,通过数学变换将共炮域公式推广到共偏移距域,再从炮点和检波点分别进行射线追踪,最后采用数据驱动的方式进行成像。一方面,考虑到地下介质的各向异性,引入了各向异性射线追踪方程;另一方面,根据有效信号和干扰信号在τ-p域中的相干性差异,在高斯束偏移过程中对地震信号进行控制,降低偏移剖面中的随机噪声,提高同相轴的连续性,最终实现了一种VTI介质共偏移距域数据驱动控制束偏移理论方法。在实现算法的基础上,通过各向异性洼陷模型、修改的SEG/Hess VTI模型及海上实际资料成像试处理,结果表明:各向异性参数对共偏移距道集中的大偏移距信息成像质量改善明显;当地层各向异性不能忽略时,新方法能够更加准确地恢复地下的复杂构造;新方法能够在一定程度上提高低信噪比数据的偏移成像效果。     

Generalized solution to the anisotropic Mandel's problem

Existing solutions to Mandel's problem focus on isotropic, transversely isotropic, and orthotropic materials, the last two of which have one of the material symmetry axes coincide with the vertical loading direction. The classical plane strain condition holds for all these cases. In this work, analytical solution to Mandel's problem with the most general matrix anisotropy is presented. This newly derived analytical solution for fully anisotropic materials has all the three nonzero shear strains. Warping occurs in the cross sections, and a generalized plane strain condition is fulfilled. This solution can be applied to transversely isotropic and orthotropic materials whose material symmetry axes are not aligned with the vertical loading direction. It is the first analytical poroelastic solution considering mechanical general anisotropy of elasticity. The solution captures the effects of material anisotropy and the deviation of the material symmetry axes from the vertical loading direction on the responses of pore pressure, stress, strain, and displacement. It can be used to match, calibrate, and simulate experimental results to estimate anisotropic poromechanical parameters. This generalized solution is capable of reproducing the existing solutions as special cases. As an application, the solution is used to study the responses of transversely isotropic and orthotropic materials whose symmetry axes are not aligned with the vertical loading direction. Examples on anisotropic shale rocks show that the effects of material anisotropy are significant. Mandel-Cryer's effects are highly impacted by the degree of material anisotropy and the deviation of the material symmetry axes from the vertical loading direction.     

Thermoplasticity and strain localization in transversely isotropic materials based on anisotropic critical state plasticity

Geomaterials such as soils and rocks are inherently anisotropic and sensitive to temperature changes caused by various internal and external processes. They are also susceptible to strain localization in the form of shear bands when subjected to critical loads. We present a thermoplastic framework for modeling coupled thermomechanical response and for predicting the inception of a shear band in a transversely isotropic material using the general framework of critical state plasticity and the specific framework of an anisotropic modified Cam–Clay model. The formulation incorporates anisotropy in both elastic and plastic responses under the assumption of infinitesimal deformation. The model is first calibrated using experimental data from triaxial tests to demonstrate its capability in capturing anisotropy in the mechanical response. Subsequently, stress‐point simulations of strain localization are carried out under two different conditions, namely, isothermal localization and adiabatic localization. The adiabatic formulation investigates the effect of temperature on localization via thermomechanical coupling. Numerical simulations are presented to demonstrate the important role of anisotropy, hardening, and thermal softening on strain localization inception and orientation. Copyright © 2016 John Wiley & Sons, Ltd.     

中国大陆地壳与上地幔地震各向异性研究

地壳与上地幔各向异性研究在地球动力学领域有广泛的应用和重要的科学意义,其研究有助于许多地质和地球物理基本问题的解释.中国大陆的各向异性研究,在地壳与上地幔两个尺度上都有许多重要进展,并在基本理论和岩石实验研究中取得进步.文章认为,各向异性研究揭示的科学问题将会是广泛的和深远的.     

Vp/Vs Anisotropy and Implications for Crustal Composition Identification and Earthquake Prediction

The ratio of P- to S-wave velocities (Vp/Vs) is regarded as one of the most diagnostic properties of natural rocks. It has been used as a discriminant of composition for the continental crust and provides valuable constraints on its formation and evolution processes. Furthermore, the spatial and temporal changes in Vp/Vs before and after earthquakes are probably the most promising avenue to understanding the source mechanics and possibly predicting earthquakes. Here we calibrate the variations in Vp/Vs in dry, anisotropic crustal rocks and provide a set of basic information for the interpretation of future seismic data from the Wenchuan earthquake Fault zone Scientific Drilling (WFSD) project and other surveys. Vp/Vs is a constant (φ0) for an isotropic rock. However, most of crustal rocks are anisotropic due to lattice-preferred orientations of anisotropic minerals (e.g., mica, amphibole, plagioclase and pyroxene) and cracks as well as thin compositional layering. The Vp/Vs ratio of an anisotropic rock measured along a selected pair of propagation-vibration directions is an apparent value (φij) that is significantly different from the value for its isotropic counterpart (φ0). The usefulness of apparent Vp/Vs ratios as a diagnostic of crustal composition depends largely on rock seismic anisotropy. A 5% of P- and S-wave velocity anisotropy is sufficient to make it impossible to determine the crustal composition using the conventional criteria (Vp/Vs<1.756 for felsic rocks, 1.756l.944 fluid-tidied porous/fractured or partially molten rocks) if the information about the wave propagation-polarization directions with respect to the tectonic framework is unknown. However, the variations in Vp/Vs measured from borehole seismic experiments can be readily interpreted according to the orientations of the ray path and the polarization of the shear waves with respect to the present-day principal stress directions (I.e., the orientation of cracks) and the frozen fabric (I.e., foliation and lineation).     

晶质岩石的地震波性质及其地质、地球物理意义

人类有关地球内部物质成份、结构和物理状态的认识绝大多数来自于地震波的资料,而地震波资料的正确解释又离不开岩石地震波性质的高温高压实验研究。地壳和地幔岩石的地震波性质和各向异性是当今世界地学前沿研究的一项重要内容,该领域横跨了地震学、岩石物理学和构造地质学,在过去一段时间里非常活跃并取得许多重要的成果。本文仅将作者近年来研究岩石地震波性质（例如,岩石波速滞后性,地震波速随围压的变化规律,纵、横波速之间的关系,泊松比等）的部分进展做一简扼的综述,并谨以此纪念已故中国科学院院士张文佑先生诞辰100周年。     

Arbitrary Lagrangian Eulerian based finite element analysis of cone penetration in soft clay

This paper presents a numerical model for the analysis of cone penetration in soft clay based on the finite element method. The constitutive behaviour of the soil is modelled by modifying an elastic, perfectly-plastic soil model obeying Von-Mises yield criterion to take into account the strain-softening, rate dependent behaviour of soft clay. Since this is a problem involving large soil deformations, the analysis is carried out using an Arbitrary Lagrangian Eulerian method where the quality of the mesh is preserved during penetration. The variation of cone resistance is examined with various parameters such as rigidity index of the soil, in situ stress anisotropy and roughness at the cone–soil interface, which influence the penetration resistance of the cone. A theoretical correlation has been developed incorporating these parameters and the results have been compared with previous correlations based on the cavity expansion theory, finite element method and strain path method. With the increase in strain-softening, relative brittleness of the soil increases and the penetration resistance is significantly reduced. With the rising strain-rate dependency, penetration resistance increases but this increase is independent of the degree of brittleness of the soil.     

Anisotropic Mean Traveltime Curves: A Method to Estimate Anisotropic Parameters from 2D Transmission Tomographic Data

We present the mathematical deduction and properties of the mean traveltime curves for homogeneous elliptical anisotropic media. These curves generalize their isotropic counterparts which have been introduced in the past as a simple data quality analysis technique at the pre-inversion stage for 2D transmission experiments, allowing the inference of prior velocity models to gain stability at the tomographic inversion. Also, the anisotropy parameters (maximum velocity, anisotropic direction and ratio) are shown to affect the shape of these curves. The degree of asymmetry of the anisotropic mean traveltime curves (displacement of the mean time and standard deviation minima from the middle of the gathering line) is related to the direction of anisotropy which can then be visually estimated. Least squares’ fitting of the anisotropic theoretical models to their experimental counterparts is an effective method to estimate at the pre-inversion stage a macroscopic elliptical anisotropic velocity model, valid at the scale of the experiment, and able to match the experimental mean traveltime distribution. Sensitivity analysis has shown that the mean curve is less prone to errors than the standard deviation curve. Parameter identification from the standard deviation curve becomes unstable for noise levels higher than 5%; data errors produce smearing of the value of the estimated anisotropy ratio and wrong directions of anisotropy biased towards zero degrees. Also, identification from the mean traveltime curve becomes stable when the maximum velocity is well constrained. Finally, this methodology is illustrated with the application to the Grimsel data set. Performing MTC analysis is always recommended since it does not need high numerical requirements, and as shown in the sensibility analysis section, errors in data can be misinterpreted as geological anisotropies. J.L. Fernández Martínez is a visiting professor at UC Berkeley, Department of Civil and Environmental Eng., CA 94720-1710.     

A destructuration theory and its application to SANICLAY model

Many natural clays have an undisturbed shear strength in excess of the remoulded strength. Destructuration modeling provides a means to account for such sensitivity in a constitutive model. This paper extends the SANICLAY model to include destructuration. Two distinct types of destructuration are considered: isotropic and frictional. The former is a concept already presented in relation to other models and in essence constitutes a mechanism of isotropic softening of the yield surface with destructuration. The latter refers to the reduction of the critical stress ratio reflecting the effect of destructuration on the friction angle, and is believed to be a novel proposition. Both the types depend on a measure of destructuration rate expressed in terms of combined plastic volumetric and deviatoric strain rates. The SANICLAY model itself is generalized from its previous form by additional dependence of the yield surface on the third isotropic stress invariant. Such a generalization allows to obtain as particular cases simplified model versions of lower complexity including one with a single surface and associative flow rule, by simply setting accordingly parameters of the generalized version. A detailed calibration procedure of the relatively few model constants is presented, and the performance of three versions of the model, in descending order of complexity, is validated by comparison of simulations to various data for oedometric consolidation followed by triaxial undrained compression and extension tests on two structured clays. Copyright © 2009 John Wiley & Sons, Ltd.