Scattering pattern analysis and true-amplitude generalized Radon transform migration for acoustic transversely isotropic media with a vertical axis of symmetry

In multi-parameter ray-based anisotropic migration/inversion, it is essential that we have an understanding of the scattering mechanism corresponding to parameter perturbations. Because the complex nonlinearity in the anisotropic inversion problem is intractable, the construction of true-amplitude linearized migration/inversion procedures is needed and important. By using the acoustic medium assumption for transversely isotropic media with a vertical axis of symmetry and representing the anisotropy with P-wave normal moveout velocity, Thomsen parameter δ and anelliptic parameter η, we formalize the linearized inverse scattering problem for three-dimensional pseudo-acoustic equations. Deploying the single-scattering approximation and an elliptically anisotropic background introduces a new linear integral operator that connects the discontinuous perturbation parameters with the multi-shot/multi-offset P-wave scattered data. We further apply the high-frequency asymptotic Green's function and its derivatives to the integral operator, and then the scattering pattern of each perturbation parameter can be explicitly presented. By naturally establishing a connection to generalized Radon transform, the pseudo-inverse of the integral operator can be solved by the generalized Radon transform inversion. In consideration of the structure of this pseudo-inverse operator, the computational implementation is done pointwise by shooting a fan of rays from the target imaging area towards the acquisition system. Results from two-dimensional numerical tests show amplitude-preserving images with high quality.     

Bayesian Markov chain Monte Carlo inversion for anisotropy of PP- and PS-wave in weakly anisotropic and heterogeneous media

A single set of vertically aligned cracks embedded in a purely isotropic background may be considered as a long-wavelength effective transversely isotropy (HTI) medium with a horizontal symmetry axis. The crack-induced HTI anisotropy can be characterized by the weakly anisotropic parameters introduced by Thomsen. The seismic scattering theory can be utilized for the inversion for the anisotropic parameters in weakly anisotropic and heterogeneous HTI media. Based on the seismic scattering theory, we first derived the linearized PP- and PS-wave reflection coefficients in terms of P- and S-wave impedances, density as well as three anisotropic parameters in HTI media. Then, we proposed a novel Bayesian Markov chain Monte Carlo inversion method of PP- and PS-wave for six elastic and anisotropic parameters directly. Tests on synthetic azimuthal seismic data contaminated by random errors demonstrated that this method appears more accurate, anti-noise and stable owing to the usage of the constrained PS-wave compared with the standards inversion scheme taking only the PP-wave into account.     

非均匀正交各向异性特征参数地震散射反演方法研究

在长波长假设条件下,水平层状地层中发育一组垂直排列的裂缝构成了等效正交各向异性介质.各向异性参数与裂缝弱度参数的估算有助于非均匀各向异性介质的各向异性特征描述,而弹性逆散射理论是非均匀介质参数反演的有效途径.基于地震散射理论,我们首先推导了非均匀正交介质中纵波散射系数方程,并通过引入正交各向异性特征参数,提出了一种新颖的正交各向异性方位弹性阻抗参数化方法.为了提高反演的稳定性与横向连续性,我们发展了贝叶斯框架下的正交各向异性方位弹性阻抗反演方法,同时考虑了柯西稀疏约束正则化和平滑模型约束正则化,最终使用非线性的迭代重加权最小二乘策略实现了各向异性特征参数的稳定估算.模型和实际资料处理表明,反演结果与测井解释数据相吻合,证明了该方法能够稳定可靠地从方位叠前地震资料中获取各向异性特征参数,减小参数估算的不确定性,为非均匀正交介质的各向异性预测提供了一种高可靠性的地震反演方法.     

Azimuthal Seismic Amplitude Variation with Offset and Azimuth Inversion in Weakly Anisotropic Media with Orthorhombic Symmetry

Seismic amplitude variation with offset and azimuth (AVOaz) inversion is well known as a popular and pragmatic tool utilized to estimate fracture parameters. A single set of vertical fractures aligned along a preferred horizontal direction embedded in a horizontally layered medium can be considered as an effective long-wavelength orthorhombic medium. Estimation of Thomsen’s weak-anisotropy (WA) parameters and fracture weaknesses plays an important role in characterizing the orthorhombic anisotropy in a weakly anisotropic medium. Our goal is to demonstrate an orthorhombic anisotropic AVOaz inversion approach to describe the orthorhombic anisotropy utilizing the observable wide-azimuth seismic reflection data in a fractured reservoir with the assumption of orthorhombic symmetry. Combining Thomsen’s WA theory and linear-slip model, we first derive a perturbation in stiffness matrix of a weakly anisotropic medium with orthorhombic symmetry under the assumption of small WA parameters and fracture weaknesses. Using the perturbation matrix and scattering function, we then derive an expression for linearized PP-wave reflection coefficient in terms of P- and S-wave moduli, density, Thomsen’s WA parameters, and fracture weaknesses in such an orthorhombic medium, which avoids the complicated nonlinear relationship between the orthorhombic anisotropy and azimuthal seismic reflection data. Incorporating azimuthal seismic data and Bayesian inversion theory, the maximum a posteriori solutions of Thomsen’s WA parameters and fracture weaknesses in a weakly anisotropic medium with orthorhombic symmetry are reasonably estimated with the constraints of Cauchy a priori probability distribution and smooth initial models of model parameters to enhance the inversion resolution and the nonlinear iteratively reweighted least squares strategy. The synthetic examples containing a moderate noise demonstrate the feasibility of the derived orthorhombic anisotropic AVOaz inversion method, and the real data illustrate the inversion stabilities of orthorhombic anisotropy in a fractured reservoir.     

Azimuthally pre-stack seismic inversion for orthorhombic anisotropy driven by rock physics

Based on the long-wavelength approximation, a set of parallel vertical fractures embedded in periodic thin interbeds can be regarded as an equivalent orthorhombic medium. Rock physics is the basis for constructing the relationship between fracture parameters and seismic response. Seismic scattering is an effective way to inverse anisotropic parameters. In this study, we propose a reliable method for predicting the Thomsen’s weak anisotropic parameters and fracture weaknesses in an orthorhombic fractured reservoir using azimuthal pre-stack seismic data. First, considering the influence of fluid substitution in mineral matrix, porosity, fractures and anisotropic rocks, we estimate the orthorhombic anisotropic stiffness coefficients by constructing an equivalent rock physics model for fractured rocks. Further, we predict the logging elastic parameters, Thomsen’s weak parameters, and fracture weaknesses to provide the initial model constraints for the seismic inversion. Then, we derive the P-wave reflection coefficient equation for the inversion of Thomsen’s weak anisotropic parameters and fracture weaknesses. Cauchy-sparse and smoothing-model constraint regularization taken into account in a Bayesian framework, we finally develop a method of amplitude variation with angles of incidence and azimuth (AVAZ) inversion for Thomsen’s weak anisotropic parameters and fracture weaknesses, and the model parameters are estimated by using the nonlinear iteratively reweighted least squares (IRLS) strategy. Both synthetic and real examples show that the method can directly estimate the orthorhombic characteristic parameters from the azimuthally pre-stack seismic data, which provides a reliable seismic inversion method for predicting Thomsen’s weak anisotropic parameters and fracture weaknesses.     

起伏地表条件下各向异性地震波最短路径射线追踪

在地震波正反演研究中,考虑起伏地表和地震各向异性具有非常重要的理论意义和实际应用价值.本文在前人研究的基础上,将最短路径追踪算法引入到起伏地表各向异性介质模型的地震波走时计算中.模型剖分时,整体模型划分成正方形单元,起伏边界附近以不规则网格逼近,进而采用非规则节点布置实现非规则网格处的最短路径计算.追踪计算中采用Sena群速度近似公式,得到各向异性地震波的走时,实现了复杂地表情况下各向异性介质模型中地震波的射线追踪.理论模型计算结果显示,本文方法能够可靠地应用于复杂各向异性介质模型,具有较高的计算精度.     

基于VTI各向异性介质的频率域海洋可控源电磁三维约束反演

近年来,海洋可控源电磁法(MCSEM)被引入油气勘探领域以降低勘探风险.在海洋环境中,受沉积因素所造成的电阻率各向异性的影响,地电模型往往会非常复杂.为更好地反映地下电性结构,本文实现了基于VTI各向异性介质的频率域海洋可控源电磁三维反演.其中,正演采用基于Yee氏交错网格的三维有限差分算法,所形成的离散线性系统通过大规模并行矩阵直接求解器(MUMPS)进行求解.反演采用基于不等式约束的有限内存BFGS(L-BFGS)算法.最后,利用VTI各向异性介质合成数据,分别进行了电阻率各向异性覆盖层和电阻率各向异性高阻层的三维反演,结果表明:(1)基于并行直接法的MCSEM非常适用于海洋电磁所特有的多场源问题;(2)针对各向异性覆盖层模型进行三维各向异性约束反演,提高了解的可靠性;(3)针对电阻率各向异性高阻层,Inline和broadside数据覆盖的反演结果对异常体位置有很好的反映.     

沙漠绿洲地区夏季地表能量收支的数值模拟

本文在MSPAS(Modified Soil-Plant-Atmosphere Scheme)的基础上,引入了一个有效的晴天大气辐射传输方案,建立了一个能在物理上真实地模拟陆气相互作用及其反馈机制的二维模式MLAIM（Modified Land Atmosphere Interaction Model）.本文利用HEIFE实验的观测资料对MLAIM的模拟结果进行了检验,对其中不合理的部分进行了分析,指出了在干旱半干旱区陆面过程参数修正的必要性,对干旱半干旱区土壤水分传输以及大气近地面层湍流输送的参数化方案进行了改进.改进后的模式能够较好地模拟夏季连续晴天条件下沙漠的地表能量收支,因此,本文利用MLAIM研究了绿洲对其周围沙漠地表能量收支的影响,并对地表能量收支各分量之间的相互作用进行了分析.结果表明,绿洲向其下风向沙漠的水汽输送是导致其上下风向沙漠间地表能量收支差异的最重要的因子.     

Stability of source mechanisms inverted from P‐wave amplitude microseismic monitoring data acquired at the surface

We study the stability of source mechanisms inverted from data acquired at surface and near‐surface monitoring arrays. The study is focused on P‐wave data acquired on vertical components, as this is the most common type of acquisition. We apply ray modelling on three models: a fully homogeneous isotropic model, a laterally homogeneous isotropic model and a laterally homogeneous anisotropic model to simulate three commonly used models in inversion. We use geometries of real arrays, one consisting in surface receivers and one consisting in ‘buried’ geophones at the near‐surface. Stability was tested for two of the frequently observed source mechanisms: strike‐slip and dip‐slip and was evaluated by comparing the parameters of correct and inverted mechanisms. We assume these double‐couple source mechanisms and use quantitatively the inversion allowing non‐double‐couple components to measure stability of the inversion. To test the robustness we inverted synthetic amplitudes computed for a laterally homogeneous isotropic model and contaminated with noise using a fully homogeneous model in the inversion. Analogously amplitudes computed in a laterally homogeneous anisotropic model were inverted in all three models. We show that a star‐like surface acquisition array provides very stable inversion up to a very high level of noise in data. Furthermore, we reveal that strike‐slip inversion is more stable than dip‐slip inversion for the receiver geometries considered here. We show that noise and an incorrect velocity model may result in narrow bands of source mechanisms in Hudson's plots.     

Interval anisotropic parameter estimation from walkaway vertical seismic profiling data

This paper presents a new explicit method for the estimation of layered vertical transverse isotropic (VTI) anisotropic parameters from walkaway VSP data. This method is based on Dix‐type normal moveout (NMO) inversion. To estimate interval anisotropic parameters above a receiver array, the method uses time arrivals of surface‐related double‐reflected downgoing waves. A three‐term NMO approximation function is used to estimate NMO velocity and a non‐hyperbolic parameter. Assuming the vertical velocity is known from zero‐offset VSP data, Dix‐type inversion is applied to estimate the layered Thomsen anisotropic parameters ?, δ above the receivers array. Model results show reasonable accuracy for estimates through Dix‐type inversion. Results also show that in many cases we can neglect the influence of the velocity gradient on anisotropy estimates. First breaks are used to estimate anisotropic parameters within the walkaway receiver interval. Analytical uncertainty analysis is performed to NMO parameter estimates. Its conclusions are confirmed by modelling.     

各向异性TI介质qP反射波走时层析成像

地震走时层析成像是反演地层各向异性参数分布的有效方法,但是关于地震各向异性介质走时层析成像的研究并不多,其技术远远没有达到成熟的阶段.在野外数据采集时,地表反射波观测方式相对井间和垂直地震剖面观测方式的成本更低,利用qP反射波走时反演各向异性参数具有更加广泛的实用价值.本文实现的TI介质地震走时层析成像方法结合了TI介质反射波射线追踪算法、走时扰动方程和非线性共轭梯度算法,它可以对任意强度的TI介质模型进行反演,文中尝试利用qP反射波走时重建TI介质模型的参数图像.利用qP反射波对层状介质模型和块状异常体模型进行走时反演,由于qP波相速度对弹性模量参数和Thomsen参数的偏微分不同,所以可以分别反演弹性模量参数和Thomsen参数.数值模拟结果表明:利用qP反射波可以反演出TI介质模型的弹性模量参数与Thomsen参数,不同模型的走时迭代反演达到了较好的收敛效果,与各向同性介质走时反演结果相比较,各向异性介质走时反演结果具有较好的识别能力.     

DInSAR沙漠形变监测应用研究(英文)

浑善达克沙地是我国主要的沙尘暴源区之一,近代浑善达克沙地的生态环境脆弱,沙化日益严重,已成为困扰北京沙尘的主要源头之一,严重影响了京津、乃至整个华北地区的生态安全以及当地生态、社会和经济的可持续发展,这一问题引起广泛的关注。DInSAR技术是重要的监测沙漠高程变化的技术,本文利用DInSAR技术,在对比分析6景ERS-2和EnviSat雷达图像基础上,对2004年10月26日和2005年10月11日两幅EnviSatASAR影像进行了干涉处理,经过配准、平地效应改正、滤波、相位解缠和地理编码,得到高程形变模型。根据研究结果,大部分地区高程是降低的,少数地区高程是升高的,这与2005年春天沙尘暴高强度发生的情况相一致,推广利用DInSAR技术实现沙漠地表形变监测具有重要作用。     

电磁勘探中各向异性研究现状和展望

电各向异性广泛存在于自然界中,已成为电磁勘探资料解释中不可忽视的因素.特别是在一些沉积岩地区,由于层理发育导致地下介质电阻率随电流方向发生变化,表现出很强的各向异性导电性.此时,利用各向同性模型进行电磁数据解释将引起很大的误差.本文回顾了电磁勘探中电各向异性研究的历史,简要介绍电各向异性成因、数值模拟中的数学描述及电各向异性介质中电磁场正反演模拟方法等研究现状,总结电各向异性在大地电磁、可控源电磁、航空电磁和感应测井等方面的最新研究进展及未来的挑战.随着多维电磁数据采集及正反演模拟技术的进步和计算能力的提高,给电各向异性信息提取和应用带来了新的契机.电各向异性研究将在矿产资源勘查、油气存储与运移特征分析、地下水和地热、环境工程地质调查、大地构造及地质灾害预测等方面发挥积极作用.     

天祝-永登地区大地电磁资料的一维各向异性反演

对分布在甘肃天祝、永登一带的10个大地电磁测点资料进行了一维各向异性介质反演研究，并与二维对称各向异性介质反演结果进行了对比。结果表明在大多数情况下一维各向异性反演有可能揭示出地壳介质电导率各向异性结构的主要特征；同时利用某些平均阻抗和特征阻抗的一维各向同性介质反演亦可获得较为可信的电性分层界面。反演所获得之介质电导率可能为各主方向电导率主值的某种平均。     

Retrieval of source parameters of an event of the 2000 West Bohemia earthquake swarm assuming an anisotropic crust

We propose an inversion scheme for retrieval of characteristics of seismic point sources, which in contrast to common practice, takes into account anisotropy. If anisotropy is neglected during inversion, the moment tensors retrieved from seismic waves generated by sources situated in anisotropic media may be biased. Instead of the moment tensor, the geometry of the source is retrieved directly in our inversion; if necessary, the moment tensor can be then determined from the source geometry aposteriori. The source geometry is defined by the orientation of the slip vector and the fault normal as well as the strength of the event given by the size of the slip and the area of the fault. This approach allows direct interpretation of the source geometry in terms of shear and tensile faulting. It also makes possible to identify volumetric source changes that occur during rupturing. We apply the described algorithm to one event of the 2000 West Bohemia earthquake swarm episode. For inversion we use information of the direct P waves. The structure is approximated by three different models determined from travel-time observations. The models are inhomogeneous isotropic, inhomogeneous anisotropic, and homogeneous anisotropic. For these models we obtain seismic moments M_T = 3.2 − 3.8 × 10¹⁴ Nm and left-lateral near-vertical oblique normal faulting on a N-S trending rupture surface. The orientation of the rupture surface is consistent with fault-plane solutions of earlier studies and with the spatial distribution of other events during this swarm. The studied event seems to be accompanied by a small amount of crack opening. The amount of crack opening is slightly reduced when the inhomogeneous anisotropic model is assumed, but it persists. These results and additional independent observations seem to indicate that tensile faulting occurs as a result of high fluid pressure.     

Research Note: Insights into the data dependency on anisotropy: an inversion prospective

While velocity contrasts are responsible for most of the events recorded in our data, the long wavelength behavior of the velocity model is responsible for the geometrical shape of these events. For isotropic acoustic materials, the wave dependency on the long (wave propagation) and short (scattering) wavelength velocity components is stationary with the propagation angle. On the other hand, in representing a transversely isotropic with a vertical symmetry axis medium with the normal moveout velocity, the anellepticity parameter η, the vertical scaling parameter δ, and the sensitivity of waves vary with the polar angle for both the long and short wavelength features of the anisotropic dimensionless medium parameters (δ and η). For horizontal reflectors at reasonable depths, the long wavelength features of the η model is reasonably constrained by the long offsets, whereas the short wavelength features produce very week reflections at even reasonable offsets. Thus, for surface acquired seismic data, we could mainly invert for smooth η responsible for the geometrical shape of reflections. On the other hand, while the δ long wavelength components mildly affects the recorded data, its short wavelength variations can produce reflections at even zero offset, with a behavior pattern synonymous to density. The lack of the long wavelength δ information will mildly effect focusing but will cause misplacement of events in depth. With low enough frequencies (very low), we may be able to recover the long wavelength δ using full waveform inversion. However, unlike velocity, the frequencies needed for that should be ultra‐low to produce long‐wavelength scattering‐based model information as δ perturbations do not exert scattering at large offsets. For a combination given by the horizontal velocity, η, and ε, the diving wave influence of η is absorbed by the horizontal velocity, severely limiting the η influence on the data and full waveform inversion. As a result, with a good smooth η estimation, for example, from tomography, we can focus the full waveform inversion to invert for only the horizontal velocity and maybe ε as a parameter to fit the amplitude. This is possibly the most practical parametrization for inversion of surface seismic data in transversely isotropic with vertical symmetry axis media.     

复杂山区初至波层析反演静校正

提高静校正精度是取得复杂山区良好地震成像的一个重要条件.而建立在水平折射面假设基础之上折射波静校正方法,无论是假设前提还是实际应用效果,都不适应于地表剧烈起伏,速度纵、横向变化大的复杂区.为此本文提出使用初至波层析反演静校正方法,即利用地震记录中初至旅行时反演出表层速度模型,计算出炮点和检波点的静校正量.通过正演模拟数据和实际资料的验证,很好的解决了复杂地表引起的静校正问题.     

Simultaneous elastic parameter inversion in 2-D/3-D TTI medium combined later arrival times

Traditional traveltime inversion for anisotropic medium is, in general, based on a “weak” assumption in the anisotropic property, which simplifies both the forward part (ray tracing is performed once only) and the inversion part (a linear inversion solver is possible). But for some real applications, a general (both “weak” and “strong”) anisotropic medium should be considered. In such cases, one has to develop a ray tracing algorithm to handle with the general (including “strong”) anisotropic medium and also to design a non-linear inversion solver for later tomography. Meanwhile, it is constructive to investigate how much the tomographic resolution can be improved by introducing the later arrivals. For this motivation, we incorporated our newly developed ray tracing algorithm (multistage irregular shortest-path method) for general anisotropic media with a non-linear inversion solver (a damped minimum norm, constrained least squares problem with a conjugate gradient approach) to formulate a non-linear inversion solver for anisotropic medium. This anisotropic traveltime inversion procedure is able to combine the later (reflected) arrival times. Both 2-D/3-D synthetic inversion experiments and comparison tests show that (1) the proposed anisotropic traveltime inversion scheme is able to recover the high contrast anomalies and (2) it is possible to improve the tomographic resolution by introducing the later (reflected) arrivals, but not as expected in the isotropic medium, because the different velocity (qP, qSV and qSH) sensitivities (or derivatives) respective to the different elastic parameters are not the same but are also dependent on the inclination angle.     

基于交叉梯度结构约束的大地电磁主轴各向异性并行三维反演

各向异性介质对大地电磁观测数据的影响往往不可忽略,因此需要提高大地电磁各向异性三维反演的可靠性和有效性.为了满足大地电磁各向异性三维反演的需求,本文研究了一种基于交叉梯度结构约束的大地电磁主轴各向异性并行三维反演算法.根据大地电磁平面波理论假设,正演方程采用背景场与二次场分离的计算方式,二次场利用交错网格有限差分法求解.由于各向异性反演的多解性,本文将各向异性介质简化为主轴各向异性,并在此基础上进一步采用有限内存拟牛顿LBFGS法实现三维各向异性反演.为了提高各向异性反演的分辨率,反演目标函数中引入交叉梯度项,利用先验的结构信息,对三个方向的电阻率参数进行结构约束,最终的反演进一步利用MPI(Message Passing Interface,消息传递接口)技术实现分频并行计算,测试结果显示并行接近线性加速比.     

从高频垂测反射回波散射谱反演电离层不均匀结构方法初探

在背景各向异性不均匀的随机起伏电离层中高频电波散射理论的基础上,进行了不同漂移速度方向下高频垂测反射回波的相位和幅度功率谱的数值计算．结合理论分析,提出了根据非常波和寻常波的相位和幅度功率谱的谱线振荡形态以及变化趋势的分析来估计电离层不均匀体的谱结构以及漂移速度的大小及方向等参量的新方法．