地球内部各圈层介质的地球各向异性与地球动力学

地球内部各圈层的结构和构造在本质上是不均匀的各向异性的,各向异性的研究不论是对深化地球本体的认识,还是在资源、能源、灾害和环境以及全球变化的研究中均具有重要作用。本文基于这一科学问题的本质性和重要性着重讨论了以下５个问题：１）各向异性介质中地震波场效应与判据;２）复杂介质地震各向异性－Ｓ波分裂的成因分析;３）当今Ｓ波分裂监测技术的进展与探索;４）地震各向异性在不同科学问题研究中的效能;５）地震各向     

地震各向异性研究进展

简要介绍了各向异性介质中地震波传播理论的最新研究进展，总结了地震各向异性在地球动力学，地震灾害监测以及资源探查等方面的应用研究新进展。目前，主要通过数值模拟和物理模拟来研究地球内部典型各向异性介质中波的传播规律；利用横波分裂，Pn波的方位变化，面波成像方法研究地球内部的结构和物质交换，应力场的变化以及地震的孕育和发生过程。在资源探查中，地震各向异性理论主要用于解决测线不闭合，中长排列动校，时深转换以及裂缝分布等问题。     

运用地球动力合成作用原理宏观预测全球火山、地震

控制地球运动或形成全球性受力状态的天然力为地球动力,包括地心引力(重力)、地球等速自转力、地球公转变速惯性力和星际引力(增量)4种力。(1)重力即地心引力。地球自转前的地心引力为原始重力P_0,地球自转后的地心引力     

东亚地震各向异性——2006年西太平洋地球物理会议专题介绍

2006年7月24—27日在北京召开了2006年西太平洋地球物理会议,本文简介了这次会议的概况,介绍了东亚地震各向异性专题中的最新研究进展。     

东亚大陆地震的地球自转动力观测资料分析

东亚大陆大地震的活动带走向、活动方式、震源主压应力方向、总迁移方向,沿纬度的分布和发震频度随时间的变化,均与地球自转速率变化有成因联系。本文从地球自转加速、减速、匀速的变化趋势进行这方面的观测资料分析,证明地球速率变化是东亚大陆地震的基本动力来源     

上地幔各向异性研究综述及在川滇地区的研究思路

简要介绍了地球介质各向异性研究的发展历程，重点介绍了国内外有关上地幔各向异性的研究工作及成果；阐述了上地幔各向异性的起源及其地球动力学意义；较详细地分析了上地幔各向异性研究的各种方法及其优缺点；论述了川滇地区各向异性研究的意义和目的，并提出了整个研究工作的内容和步骤。     

"空间技术揭示的古吉拉特地震的地球系统过程"国际研讨会

“空间技术揭示的古吉拉特地震的地球系统过程”国际研讨会于2003年1月27－29日在印度坎普尔举行。来自印度、美国、中国、德国和俄罗斯约80位专家参加了该研讨会。会计收到约70篇研究报告。文章介绍了会议概况、学术报告和专题辩论的情况。     

Complex frequency-shifted multi-axial perfectly matched layer for frequency-domain seismic wavefield simulation in anisotropic media

Seismic anisotropy has an important influence on seismic data processing and interpretation. Although the frequency-domain seismic wavefield simulation has a problem of solving the large scale linear sparse matrix due to the computational limitations, it has some advantages over the time-domain seismic wavefield simulation including efficient inversion using only a limited number of frequency components and easy implementation of multiple sources. To accurately simulate seismic wave propagation in the frequency domain, we also need to choose the absorbing boundary conditions to absorb artificial reflections from edges of the model as we do in the time domain. Compared with the classical boundary conditions including the perfectly matched layer and complex frequency-shifted perfectly matched layer, the complex frequency-shifted multi-axial perfectly matched layer has been proven to effectively suppress the unwanted reflections at grazing incidence and solve the instability problem in the time-domain seismic numerical modelling in anisotropic elastic media. In this paper, we propose to extend the complex frequency-shifted multi-axial perfectly matched layer absorbing boundary condition to the frequency-domain seismic wavefield simulation in anisotropic elastic media. To test the validity of our proposed algorithm, we compare the results (snapshots and seismograms) of the frequency-domain seismic wavefield simulation with those of the time-domain modelling. The model studies indicate that the complex frequency-shifted multi-axial perfectly matched layer absorbing boundary condition is stable in the frequency-domain seismic wavefield simulation in anisotropic media, and provides better absorbing performance than the complex frequency-shifted perfectly matched layer boundary condition.     

Construction of regional and local seismic anisotropic structures from wide-angle seismic data: crustal deformation in the southeast of China

We present a method to obtain spatial distributions of seismic anisotropy associated with regional stress and local faulting in the crust from wide-angle seismic data. The method contains three steps. The first step consists of obtaining radial- and transverse-component seismic sections using a pre-stack depth migration algorithm from the S-wave velocity model determined by conventional interpretation of picked intra-crustal seismic events. In the second step, we compute time delays between split shear-waves and polarizations of fast split shear-waves by minimizing the transverse-component seismic energy. The time delay and polarization in each layer are derived using a layer-stripping method. The final step is to estimate the average splitting parameters along the whole profile. Thus, the average time delay and polarization can be regarded as caused by the effects owing to regional structure and stress fields, whereas the residual values of the splitting parameters are considered to be related to local structures and local faulting. Our method allows us to construct multi-layer anisotropic images, which may later be interpreted in terms of intra-layer coupling/decoupling or deformation. We present results from a set of three-component seismic data acquired by a controlled source experiment in the southeast region of China. The results demonstrate that the average polarizations and time delays are consistent with the direction and strength of the stress field, and their lateral variations related to local anisotropy match the spatial distribution of surface faulting crossing the acquisition seismic profile.     

Looking for layered anisotropic structures in the mantle beneath the northern Apennines

Competing geodynamic scenarios proposed for northern Apennines (Italy) make very different predictions for the orientation of strain in the upper mantle. Constraints on the pattern are offered by observations of seismic anisotropy. Previous study of the anisotropy beneath the northern Apennines used birefringence of core-refracted shear waves (SKS phases), and demonstrated the presence of two domains: Tuscan and Adria. In the transition between the two domains, across the Apennines orogen, anisotropy measurements reflect a complex deep structure. To define better the upper-mantle structure beneath this area we analyze seismological data recorded by a set of seismic stations that operated for 3 years, between 2003 and 2006, located in the outer part of the Apennines belt, in the Adria terrane, collected by the RETREAT Project. Directionally distributed sets of SKS records were inverted for layered anisotropic structures with a well-tested method, adding new results to previous hypotheses for this area. New data analysis argues for two-layer anisotropy for sites located on the Apennines wedge and also one site in the Tuscan terrane. Beneath the wedge an upper layer with nearly north-south fast polarization pervades the lithospheric mantle, while at depth a nearly NW–SE Apennines-parallel direction is present in the lower layer. Beneath Tuscany a shallower NW–SE direction and a deeper E–W one suggest the deeper strain from active slab retreat, with a mantle-wedge circulation (i.e. an east–west corner flow), overlain by an Apennines-parallel fast polarization that could be a remnant of lower-crust deformation.     

Any P-wave kinematic algorithm for vertical transversely isotropic media can be ADAPTed to moderately anisotropic media of arbitrary symmetry type

Most of P-wave anisotropic kinematic algorithms (modeling, processing, and inversion) have been developed for the case of Transverse Isotropy (TI). Does it mean that when dealing with more complex symmetry types (Arbitrarily tilted TI, orthorhombic, monoclinic or even triclinic), all these algorithms are irrelevant? In fact, not at all. It has recently been demonstrated that in 2D geometry any qP-wave TI kinematic algorithm can be simply generalized to the case of monoclinic symmetry using the so-called Azimuthally Dependent Anisotropy Parameter Transformation (ADAPT), assuming moderate anisotropy. The extension of the technique to the case of arbitrary anisotropy type (triclinic) is achieved in this paper. The method is successfully checked for seismic modeling in a full 2D model with layers of contrasted anisotropy types and with arbitrary vertical and horizontal velocity variations (non-constant gradient). Typically, the approximate travel times using ADAPT differ from the exact travel times by a few milliseconds for total travel times of the order of a few seconds. Applications to seismic processing are also described. The simplicity of the procedure and the generality of the applicability of the ADAPT recipe are striking and very convenient for practical applications. They certainly deserve further analysis.     

鄂尔多斯块体北缘与西缘地区地壳各向异性特征

本研究使用内蒙古自治区数字测震台网2010年1月至2017年10月区域小地震的波形记录资料,采用SAM方法,进行了地壳剪切波分裂的分析,得到鄂尔多斯块体北缘与西缘地区地壳介质地震各向异性的初步研究结果.根据15个台站161个有效地震记录的分析,鄂尔多斯块体北缘与西缘地区的快剪切波平均偏振方向为NE44.4°±38.4°,慢剪切波平均时间延迟为1.7±1.6 ms·km^-1.研究区域的快剪切波偏振显示出两个优势方向,一个是NE方向,另一个是近NS方向.区内的逆冲凸起与走滑正倾断层构造对剪切波分裂产生了直接的影响,造成了剪切波分裂参数的复杂分布,反映了剪切波分裂参数受到区域应力和构造共同作用的影响.鄂尔多斯块体北缘的快波偏振特征有NE和近NS两个优势偏振方向,其东区与西区的快剪切波偏振表现出明显不同的特征.东区的第一快剪切波优势偏振方向为NE,第二快剪切波优势偏振方向为近NS;西区的第一快剪切波优势偏振方向为近EW,第二快剪切波优势偏振方向为近NS.鄂尔多斯块体北缘的区域背景主压应力方向可能总体上为近NS方向,但空间分布有差异,东区NE方向的优势偏振与西区近EW方向的优势偏振更可能反映了断裂与构造的影响.鄂尔多斯块体西缘的快剪切波偏振特征显示出非常清楚的NE向的优势偏振方向,近NS向的优势偏振方向则不太明显,反映出该地区复杂构造对各向异性分布的影响.慢波时间延迟呈现出西低东高的特点,时间延迟的高值出现在鄂尔多斯块体北缘的东部,时间延迟的这种西低东高的各向异性强度变化,可能反映了区域构造活动西强东弱的特性.

     

Advances in the deep tectonics and seismic anisotropy of the Lijiang-Xiaojinhe fault zone in the Sichuan-Yunnan Block,Southwestern China

The Sichuan-Yunnan Block (SYB) is located at the SE margin of the Qinghai-Tibetan Plateau (TP). Under the influence of the southeastward movement of material originated from the TP, intense crustal deformation, frequent seismic activity, and complex geological structures are observed in the SYB. The Lijiang-Xiaojinhe fault (LXF) goes through the central part of the SYB, dividing it into two blocks from north to south, and forming an intersecting fault system with the surrounding faults. This paper firstly introduces the morphology and the nature of the LXF, the distribution of the regional surface displacements and the focal mechanisms, and then analyzes the medium deformation and the effects of faults. Moreover, according to the regional tectonics and geophysical patterns, the paper discusses the characteristics of the north-south blocks of the SYB and the abrupt change of deep structure along the LXF zone. Since seismic anisotropy is an essential property for detecting crustal stress, deep structures and dynamical mechanisms, this paper is dedicated to the advances in seismic anisotropy at different depths and different scales in the study area. There are noteworthy differences in the anisotropic features between the north part and the south part of the SYB, possibly associated with a clear boundary adjacent to the LXF. Such phenomenon suggests some close correlation between anisotropic zoning boundary and the LXF, although this boundary is not consistent with the LXF in strike. The results from the deformation of the crust and the upper mantle elucidate the distribution patterns of the crust-mantle coupling in the north part and the crust-mantle decoupling in the south part, even though this conclusion needs to be further verified by more studies. Presently, the scientific understanding of the deep tectonics and the media deformation around the “generalized” LXF i.e. the LXF with the Jinpingshan fault on its eastern side, is still insufficient, and related equivocal topics deserve more in-depth studies.     

湖北地区上地幔各向异性及其动力学意义

SKS波分裂测量是研究大陆地幔的形变特征、探索大陆动力学和演化过程的重要工具.本文选用湖北数字地震台网(HB台网)和中国数字化宽频带地震台网(CB台网)提供的三分量宽频带数字化地震资料,使用旋转相关法和切向能量最小法,计算得到了湖北地区23个地震台站下方上地幔各向异性参数.结果表明,快波偏振方向较为复杂:大别造山带内,呈平行于大别造山带构造走向的NWW向,且与绝对板块运动方向基本一致,反映印支—燕山期扬子板块与华北板块南北向挤压碰撞造山的过程,而地幔流动则是造成这种南北挤压的力源.扬子板块内NWW向虽与绝对板块运动趋势一致,但比较离散,受地幔对流影响不明显,而NE向各向异性主要与燕山后期之后应力场的转变有关.分析认为,该地区的上地幔各向异性为地幔流动和固结在岩石圈中的"化石各向异性"共同作用的结果.根据上地幔各向异性方向与主要构造带的走向平行,可以作为地壳变形与上地幔连贯变形的证据.整个研究区域得到的分裂延迟时间在0.5~1.35 s之间,厚度在58~155 km之间.     

各向异性介质中的弹性阻抗及其反演

地震反演已成为油藏描述中的重要组成部分.绝大多数的常规地震反演是叠后地震数据体上进行,很少考虑各向异性存在的情况.随着勘探开发的发展,地震各向异性和叠前地震波阻抗反演引起了人们极大关注.本文在各向同性介质中弹性阻抗研究基础上,推导出了各向异性介质中的弹性阻抗方程,提出了地震各向异性介质中用弹性阻抗进行储层参数描述的技术路线和框架,并对反演过程中存在的问题进行了有益探讨.     

基于walkaway VSP下行P波的TTI各向异性参数反演

本文基于弱各向异性（WA）介质的正反演公式和qP波的坐标变换,推导了利用qP波反演任意倾斜对称轴的横向各向同性（TTI）介质的各向异性参数和对称轴方向的公式.理论和数值实验表明,利用2个相互正交的变井源距垂直地震剖面（walkaway VSP）可以完全确定钻井中TTI介质qP波的3个WA参数和对称轴的2个方向参数.我们完成了几个由不同数量剖面组成的walkaway VSP模拟实验,使用TTI模型和一般各向异性模型对模拟数据进行了反演,证明了反演公式的正确性和可靠性.使用这些公式,对来自Java Sea的由3条剖面组成的walkaway VSP观测数据进行了各向异性反演,获得了钻井中接收点处介质的WA参数.

     

Anisotropic permeability in fractured reservoirs from frequency‐dependent seismic Amplitude Versus Angle and Azimuth data

Attempts have previously been made to predict anisotropic permeability in fractured reservoirs from seismic Amplitude Versus Angle and Azimuth data on the basis of a consistent permeability‐stiffness model and the anisotropic Gassmann relations of Brown and Korringa. However, these attempts were not very successful, mainly because the effective stiffness tensor of a fractured porous medium under saturated (drained) conditions is much less sensitive to the aperture of the fractures than the corresponding permeability tensor. We here show that one can obtain information about the fracture aperture as well as the fracture density and orientation (which determines the effective permeability) from frequency‐dependent seismic Amplitude Versus Angle and Azimuth data. Our workflow is based on a unified stiffness‐permeability model, which takes into account seismic attenuation by wave‐induced fluid flow. Synthetic seismic Amplitude Versus Angle and Azimuth data are generated by using a combination of a dynamic effective medium theory with Rüger's approximations for PP reflection coefficients in Horizontally Transversely Isotropic media. A Monte Carlo method is used to perform a Bayesian inversion of these synthetic seismic Amplitude Versus Angle and Azimuth data with respect to the parameters of the fractures. An effective permeability model is then used to construct the corresponding probability density functions for the different components of the effective permeability constants. The results suggest that an improved characterization of fractured reservoirs can indeed be obtained from frequency‐dependent seismic Amplitude Versus Angle and Azimuth data, provided that a dynamic effective medium model is used in the inversion process and a priori information about the fracture length is available.