One-Dimensional qP-Wave Velocity Model of the Upper Crust for the West Bohemia/Vogtland Earthquake Swarm Region

The western part of the Bohemian Massif (West Bohemia/Vogtland region) is characteristic in the relatively frequent recurrence of intraplate earthquake swarms and in other manifestations of past-to-recent geodynamic activity. In this study we derived 1D anisotropic qP-wave model of the upper crust in the seismogenic West Bohemia/Vogtland region by means of joint inversion of two independent data sets - travel times from controlled shots and arrival times from local earthquakes extracted from the WEBNET seismograms. We derived also simple 1-D P-wave and S-wave isotropic models. Reasons for deriving these models were: (a) only simplified crustal velocity models, homogeneous half-space or 1D isotropic layered models of this region, have been derived up to now and (b) a significant effective anisotropy of the upper crust in the region which was indicated recently by S-wave splitting. Both our anisotropic qP-wave and isotropic P-and S-wave velocity models are constrained by four layers with the constant velocity gradient. Weak anisotropy for P-waves is assumed. The isotropic model is represented by 9 parameters and the anisotropic one is represented by 24 parameters. A new robust and effective optimization algorithm - isometric algorithm - was used for the joint inversion. A two-step inversion algorithm was used. During the first step the isotropic P- and S-wave velocity model was derived. In the second step, it was used as a background model and the parameters of anisotropy were sought. Our 1D models are adequate for the upper crust in the West Bohemia/Vogtland swarm region up to a depth of 15 km. The qP-wave velocity model shows 5% anisotropy, the minimum velocity in the horizontal direction corresponds to an azimuth of 170°. The isotropic model indicates the V_P/V_S ratio variation with depth. The difference between the hypocentre locations based on the derived isotropic and anisotropic models was found to be several hundreds of meters.     

怀来、 延庆及其周边地区地壳 介质各向异性研究

提出了利用人工爆破P波走时反演地壳介质方位各向异性参数的方法. 在假定介质是弱各向异性介质的情况下, 使用扰动理论得到了线性化的反演公式, 其中待反演的弱各向异性参数是P波走时的线性函数. 如果在反演公式中参考走时取相同震中距接收点的P波平均走时, 那么所获得的弱各向异性参数与参考介质速度的选取无关. 反演得到的弱各向异性参数可以看作是不同震中距和不同深度范围内介质的等效弱各向异性参数. 等效弱各向异性参数在一定程度上反映了不同深度范围内水平方向相速度随方位的变化. 这种变化可能是不同时期构造应力作用的结果. 2007年中国地震局在首都圈怀来地区实施了一次大吨位人工爆破实验, 以爆破点为中心, 布设了高密度的地震观测台网和台阵. 台站相对于爆破点具有360°的全方位覆盖, 所得到的地震记录数据为研究怀来、 延庆地区地壳介质P波方位各向异性提供了必要条件. 我们通过走时反演获得了与水平方位相关的弱各向异性参数, 并对弱各向异性参数进行坐标变换, 得到了能够直观描述岩石弱各向异性的具有水平对称轴的横向各向同性介质, 给出了对应的3个独立弱各向异性参数及其对称轴方位, 讨论了介质各向异性与构造应力场的关系. 结果表明该地区地壳介质存在明显的方位各向异性, 其最大值约为4.6%.      

Seismic velocity estimation from post-critical wide-angle reflections in layered structures

The travel time inversion of wide-angle seismic data is a technique commonly used in the deep seismic sounding. We propose an application of this technique to a smaller scale of a sedimentary layer, where the characteristics of seismic observations changes significantly. Field observations confirmed by synthetic analysis recognize the dominant amplitudes of wide-angle post-critical reflections. A case study is presented in this paper, of a joint interpretation of conventional reflection seismic with reflection imaging, combined with the wide-angle travel time inversion of additional full-spread observations. A joint interpretation results in a precise recognition of the seismic velocity distribution, that is further used for the seismic depth conversion with the uncertainty analysis of the depth of the reflecting horizons. Despite the salt layer in the studied structure this method is able to precisely recognize the seismic velocities of the sub-salt structures.     

长乐—南澳断裂带两侧地壳结构差异的地震-重力联合反演

长乐一南澳断裂带出露于福建沿海地区,由于海陆过渡带的特殊性,地球物理探测受到许多限制,难以获得由陆到海的清晰而准确的深部构造形态.2014年福建省地震局采集了横跨长乐一南澳断裂带的广角反射/折射剖面(HX-6),由于观测系统的缺陷和原始资料信噪比等问题,单纯使用地震数据反演长乐一南澳断裂带的深部地壳结构有很强的不确定性,无法解答断裂带两侧地壳结构存在何种差异,影响了对断裂带构造属性和区域构造演化的正确认识.基于岩石波速和密度有良好的对应关系,地震-重力联合反演可以有效降低多解性.本文采用地震走时拟合和重力异常拟合同步进行的方法,利用最新采集的高质量P波地震走时数据与高精度实测重力数据,反演得到了连城一厦门一金门外海剖面的二维地壳波速-密度结构模型.联合反演结果显示:长乐一南澳断裂带两侧地壳厚度差约3 km,壳内分层结构和上地幔顶部波速-密度无显著变化,推断长乐一南澳断裂带是华南地块正常陆壳和台湾海峡减薄陆壳的分界,本研究结果为进一步研究该区深部构造环境和长乐一南澳断裂带的地球动力学意义提供了新的地球物理学证据.     

Reconstruction of the layer anisotropic elastic parameters and high-resolution fracture characterization from P-wave data: a case study using seismic inversion and Bayesian rock physics parameter estimation

Wide-azimuth seismic data can be used to derive anisotropic parameters on the subsurface by observing variation in subsurface seismic response along different azimuths. Layer-based high-resolution estimates of components of the subsurface anisotropic elastic tensor can be reconstructed by using wide-azimuth P-wave data by combining the kinematic information derived from anisotropic velocity analysis with dynamic information obtained from amplitude versus angle and azimuth analysis of wide-azimuth seismic data. Interval P-impedance, S-impedance and anisotropic parameters associated with anisotropic fracture media are being reconstructed using linearized analysis assuming horizontal transverse anisotropy symmetry. In this paper it is shown how additional assumptions, such as the rock model, can be used to reduce the degrees of freedom in the estimation problem and recover all five anisotropic parameters. Because the use of a rock model is needed, the derived elastic parameters are consistent with the rock model and are used to infer fractured rock properties using stochastic rock physics inversion. The inversion is based on stochastic rock physics modelling and maximum a posteriori estimate of both porosity and crack density parameters associated with the observed elastic parameters derived from both velocity and amplitude versus angle and azimuth analysis. While the focus of this study is on the use of P-wave reflection data, we also show how additional information such as shear wave splitting and/or anisotropic well log data can reduce the assumptions needed to derive elastic parameter and rock properties.     

利用深反射地震数据构建的多阶面波频散曲线反演近地表横波速度结构——以跨班公湖—怒江缝合带深反射地震资料为例

深反射地震剖面法为了获取深部结构特征常常采取大的偏移距采集数据.目前公开发表的相关资料中,鲜有利用深反射地震炮集数据获取近地表的结构特征.为此,本文通过正演测试了相关数据处理流程,即利用有限差分正演了起伏地表模型的大偏移距地震单炮弹性波场特征,通过共检波点域面波信号F-K频谱叠加构建新方法,从深反射地震数据集中提取了高品质的多阶面波频散曲线,再利用多阶面波联合反演获得了近地表的结构特征.在前述正演流程基础上,利用跨越班公湖—怒江缝合带的SinoProbe深反射地震剖面中的实际炮集数据,求取了基阶和一阶瑞利波频散曲线,联合反演后得到近地表横波速度结构.该结果与初至波走时反演获取的纵波速度结构具有较好的一致性,且在近地表的浅层分辨率较纵波速度结构特征更高,而更与已有地质认识相吻合.本文提供的相关数据处理流程表明利用深反射地震炮集数据,也能够获取近地表浅层的横波速度结构.     

华北克拉通北缘(怀来-苏尼特右旗)地壳结构

2009年,中国地质科学院地质研究所与美国俄克拉荷马大学合作实施了一条长453 km的深地震反射、宽角反射与折射、三分量反射地震联合探测剖面. 剖面南起怀来盆地,向北依次穿过燕山造山带西缘、内蒙地轴、白乃庙弧带、温都尔庙杂岩带,到达索伦缝合带. 其中,宽角反射与折射剖面采用8个0.5~1.5 t炸药震源激发,使用300套Texan单分量数字检波器接收,获得了高质量的地震资料. 通过资料分析和处理,识别出沉积层及结晶基底的折射波（Pg）、来自上地壳底界面的反射波（Pcp）,中地壳底界面的反射波（Plp）,莫霍界面的反射波（Pmp）及上地幔顶部的折射波（Pn）等5个震相. 分别采用Hole有限差分层析成像和Rayinvr算法对华北克拉通北缘及中亚造山带南部进行了上地壳P波速度结构成像和全地壳二维射线追踪反演成像. 结果显示：（1）中亚造山带地壳厚度~40 km,变化平缓,低于全球平均造山带地壳平均厚度,可能为造山后区域伸展的结果. 阴山-燕山带附近莫霍明显加深,推测其为燕山期造山过程形成的山根,但该山根很可能在后期被改造. （2）测线中部地壳上部速度较高,对应地表大面积花岗岩出露,而下地壳速度较低,速度梯度低,呈通道状,推测其可能曾为古亚洲洋向南俯冲消亡的主动陆缘,并在碰撞后演变为伸展环境下岩浆侵入的通道. （3）华北克拉通北缘与中亚造山带显示出不同速度变化特征,前者变化相对缓而后者则变化剧烈,二者的分界出现在赤峰-白云鄂博断裂附近.     

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.     

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.     

Crustal structure of northeastern Tibetan plateau and Ordos block: Waveform interpretation of the Maqen-Jingbian seismic refraction profile

The Maqen-Jingbian wide-angle seismic reflection and refraction experiment was carried out in 1998, which aims at determining detailed structure in the crust and top of the upper mantle and understanding structural relation between the northeastern Tibetan plateau and the Ordos block. The 1-D crustal models inferred by waveform inversion show strong variations in crustal structure, which can be classified into four different types: ① an Ordos platform with the Proterozoic crust and two high-velocity layers in the northeast section, ② a transitional crust between the northeastern Tibetan plateau and the Ordos block across the Haiyuan earthquake zone, ③ the Qilian orogenic zone in the central part, and 4 the Qinling orogenic zone in the southwestern section. The Moho depth increases from ~42 km to ~62 km from the NE part to the SW part of the profile. The crystalline crust consists of the upper crust and lower crust in northeastern Tibetan plateau. There is an obviously low P-wave velocity layer dipping northeastward, which is 12–13 km thick, at the bottom of the upper crust in Qinling orogenic zone and Haiyuan earthquake zone. The lower crust is characterized by alternating high and low P-wave velocity layers. Beneath Ordos block, i.e., the NE part of the profile, the crust shows quite a smooth increase in P-wave velocity down to the Moho at a depth of about 42 km.     

Linearized elastic parameter sections1

Contrasts in elastic parameters can be estimated directly from seismic data using offset-dependent information in the PP reflection coefficient. A linear approximation to the PP reflection coefficient including three coefficients is fitted to the data, and relative contrasts in various elastic parameters are obtained from linear combinations of the estimated coefficients. Linearized elastic parameter sections for the contrasts in P-wave impedance, P-wave velocity, density, plane-wave modulus, and the change in bulk modulus and shear modulus normalized with the plane-wave modulus are estimated. If the average P- to S-wave velocity ratio is known, linearized parameter sections including the contrast in the average P- to S-wave velocity ratio and a fluid factor section can be computed. Applied to synthetic data, visual comparison of the estimated and true elastic parameter sections agree qualitatively, and the results are confirmed by an analysis of the standard deviation of the estimated parameters. The parameter sections obtained by inversion of a shallow seismic anomaly in the Barents Sea are promising, but the reliability is uncertain because neither well data nor regional trends are available.     

波形反演遗传算法及其在地震测深数据解释中的应用

为了更好地利用地震测深波形数据,提出了地震体波波形反演的遗传算法. 正演使用能精确快速计算互层结构响应的广义反、透射系数理论地震图算法;反演采用遗传算法,实现了地震体波波形反演的遗传算法. 数值试验表明,该算法具有分辨壳内低速层、高低速薄互层结构和一定的抗噪能力. 青藏高原东北缘泽库、夏河、临洮3炮地震测深P波波形反演,得到了上地壳底部低速层和中、下地壳,以及上地幔顶部薄互层的细结构图象.      

Genetic algorithm in seismic waveform inversion and its application in deep seismic sounding data interpretation

A genetic algorithm of body waveform inversion is presented for better understanding of crustal and upper mantle structures with deep seismic sounding （DSS） waveform data. General reflection and transmission synthetic seismogram algorithm, which is capable of calculating the response of thin alternating high and low velocity layers, is applied as a solution for forward modeling, and the genetic algorithm is used to find the optimal solution of the inverse problem. Numerical tests suggest that the method has the capability of resolving low-velocity layers, thin alternating high and low velocity layers, and noise suppression. Waveform inversion using P-wave records from Zeku, Xiahe and Lintao shots in the seismic wide-angle reflection/refraction survey along northeastern Qinghai-Xizang （Tibeteau） Plateau has revealed fine structures of the bottom of the upper crust and alternating layers in the middle/lower crust and topmost upper mantle.     

TI介质各向异性速度多参数分析

本文采用横向各向同性（TI）介质弹性参数的Anderson表征方式,利用小偏移距同类反射波或转换波信息重建纵横波速,以中长排列同类反射波或转换反射波信息重建各向异性因子图像,分步进行深度域TI介质中P、SV波多参数各向异性速度分析,以理论模型验证了方法的可行性.最后给出宽角反射PP与SS波折合剖面的速度与各向异性因子解释结果.     

Three-dimensional velocity model of the crust of the Bohemian Massif and its effects on seismic tomography of the upper mantle

We have compiled a representative three-dimensional P-velocity model of the crust of the Bohemian Massif (BM) to provide a basis for removing effects of the crustal structure in teleseismic tomography of the upper mantle. The model is primarily based on recently published 2D velocity models from findings of wide-angle refraction and near-vertical reflection seismic profiles of CELEBRATION 2000, ALP 2002, and SUDETES 2003 experiments. The best fitting 3D model of the BM crust (NearNeighbour model) is complemented by velocities according to the reference Earth model at sites where data are sparse, which precludes creating artificial heterogeneities that are products of interpolation method. To test the model, we have performed tomographic inversions of the P-wave travel times measured during the BOHEMA II experiment and compared the results obtained with and without crustal corrections. The tests showed that the presented crustal model decreases magnitudes of velocity perturbations leaking from the crust to the mantle in the western part of the BM. The tomographic images also indicated a highvelocity anomaly in the lower crust or just beneath the crust in the Brunovistullian unit. Such anomaly is not described by our model of the crust since no seismic profile intersects this part of the unit. The tests also indicated that crustal corrections are of the great importance especially for interpretations of the uppermost mantle down to depths of about 100 km.     

Elastic reverse-time migration based on amplitude-preserving P- and S-wave separation

Imaging the PP- and PS-wave for the elastic vector wave reverse-time migration requires separating the P- and S-waves during the wave field extrapolation. The amplitude and phase of the P- and S-waves are distorted when divergence and curl operators are used to separate the P- and S-waves. We present a P- and S-wave amplitude-preserving separation algorithm for the elastic wavefield extrapolation. First, we add the P-wave pressure and P-wave vibration velocity equation to the conventional elastic wave equation to decompose the P- and S-wave vectors. Then, we synthesize the scalar P- and S-wave from the vector Pand S-wave to obtain the scalar P- and S-wave. The amplitude-preserved separated P- and S-waves are imaged based on the vector wave reverse-time migration (RTM). This method ensures that the amplitude and phase of the separated P- and S-wave remain unchanged compared with the divergence and curl operators. In addition, after decomposition, the P-wave pressure and vibration velocity can be used to suppress the interlayer reflection noise and to correct the S-wave polarity. This improves the image quality of P- and S-wave in multicomponent seismic data and the true-amplitude elastic reverse time migration used in prestack inversion.     

Joint inversion of travel times and waveforms by means of annealing algorithms. Application to a seismic refraction experiment

An algorithm of annealing is applied to a joint inversion of travel times and waveforms belonging to a synthetically generated seismic refraction experiment. The medium (crust and upper mantle) is modelled by a set of plane stratified layers and a halfspace. The obtained structure (elastic parameters and depth of layers) shows that, in spite of contaminating seismic noise and poor knowledge of the seismic source, annealing methods are a good tool in these kinds of inversion problems. We think that many characteristics of the annealing process described here could be used with real data and more sophisticated media for the crust and upper mantle of the earth than the present example.     

Sensitivity study using a genetic algorithm: inversion of amplitude variations with slowness

A sensitivity study of elastic parameters in amplitude-variation-with-slowness (AVS) for small- and large-offset seismic data is presented. In order to handle the non-linearity associated with waveform or amplitude beyond the critical slowness, an inversion algorithm based on Bayes' theory is used. A genetic algorithm was used to obtain the a posteriori probability density (PPD) function. The sensitivity analysis is performed on synthetic data containing P-wave as well as converted S-wave reflections. Four different two-layer models, which represent the typical range of AVS responses associated with the gas-sands normally encountered in exploration, were used to examine how well the elastic parameters can be inverted for different parametrizations by comparing the PPD functions. The sensitivity study results suggest that including wide-angle data in the inversion can greatly enhance the quality of inversion. The converted S-wave reflections can provide valuable extra information that can be used to extract elastic parameters. The results with noisy data demonstrate that the contrast of density and three velocity ratios can be estimated robustly with wide-angle reflection data.     

富含煤层气储层的叠前地震反演预测方法研究——以山西晋东南测区为例

叠前地震反演技术能够将地震中振幅、偏移距和入射角等多种有效信息与油气的敏感参数结合起来,在有效识别油气藏"甜点"方面虽发挥了重要的作用,但对于煤层气的有效预测仍需研究.针对煤层气与石油天然气的赋存地质条件不同,本文以山西沁水盆地煤层气地震资料为例,借鉴叠前地震反演技术的多参数预测思想,通过模型正演和实测资料处理,得到纵波、横波速度和密度等不同的数据体,利用数据体和弹性参数关系式,采用AkiRichard近似公式与叠前宽角度反演联合的方法,寻找出富含煤层气储层的敏感参数.从而实现富含煤层气储层的有效预测.经与测井曲线结果对比,吻合率较高,说明了该方法具有可行性与实用性.     

各向异性介质qP波传播描述I：伪纯模式波动方程

地球介质相对于地震波波长尺度的定向非均匀性会导致波速的各向异性,进而影响地震波场的运动学与动力学特征．各向异性弹性波动方程是描述该类介质波场传播的基本工具,在正演模拟、偏移成像与参数反演中起着关键作用．为了面向实际应用构建灵活、简便的各向异性波场传播算子,人们一直在寻求简化的各向异性波动方程．本文借鉴各向异性弹性波波型分离思想,通过对平面波形式的弹性波方程(即Christoffel方程)实施一种代表向波矢量方向投影的相似变换,推导出了一种适应任意各向异性介质、运动学上与原始弹性波方程完全等价,在动力学上突出qP波的新方程,即qP波伪纯模式波动方程．文中以横向各向同性(TI)介质为例,给出了相应的qP波伪纯模式波动方程及其声学与各向同性近似,并在此基础上开展了正演模拟和逆时偏移试验,展示了这种描述各向异性波场传播的新方程的特点与优势．