基于Walkaway VSP的Thomsen各向异性参数估计与应用(英文)

Thomsen各向异性参数的求取对于正确的时深转换和深度域偏移成像处理至关重要。相比其它各向异性参数估算方法,从VSP资料中更容易获取准确的各向异性参数用于地面地震偏移成像。本文分析研究了利用Walkaway VSP资料估算VTI介质Thomsen各向异性参数的方法,该方法基于VTI介质近偏移距动校正公式利用Walkaway VSP近偏移距初至信息求取各向异性参数δ;基于各向异性介质纵波速度Thomsen近似公式采用射线追踪时差扫描方法求取各向异性参数ε。数值模型正演表明利用该方法估算的各向异性参数误差较小。利用塔里木盆地8个方位的Walkaway VSP实际资料求取了该区深度域Thomsen各向异性参数ε和δ值,同时结合地面三维地震资料建立了较为准确的各向异性深度一速度模型用于叠前深度偏移成像,进一步提高了碳酸盐岩储层的成像精度,减小了目标地质体的深度误差。     

基于各向异性梯度的裂缝密度反演

裂缝密度是评价裂缝性储层的重要参数. 本文提出一种利用各向异性梯度计算裂缝密度的新方法,即运用裂缝介质等效理论将高陡倾角裂缝介质等效为横向各向同性介质,得到裂缝介质AVO响应特征,从而得到裂缝介质的各向异性梯度,之后由该梯度与裂缝密度的关系得到裂缝介质的裂缝密度. 实际处理中可以根据地震数据直接求取裂缝介质的各向异性梯度,最后求取裂缝介质的裂缝密度. 模型试算表明, 该方法能得到裂缝介质准确的裂缝密度. 反演结果与初始模型基本一致验证了该方法的正确性; 另外地震记录添加一定的信噪比后也能反演得到准确的裂缝密度,证明了该方法的稳定性.      

基于Born敏感核函数的VTI介质多参数全波形反演

本文基于VTI介质拟声波方程,利用散射积分原理,在Born近似下导出了速度与各向异性参数的敏感核函数,同时结合作者前期研究提出的矩阵分解算法实现了一种新的VTI介质多参数全波形反演方法.矩阵分解算法通过对核函数-向量乘进行具有明确物理含义的向量-标量乘分解累加运算实现目标函数一阶方向或二阶方向的直接求取,从而避免了庞大核函数矩阵与Hessian矩阵的存储,该方法同时可以大大降低常规全波形反演在计算二阶方向时的庞大计算量.为了克服不同参数对波场影响程度的不同,本文利用作者前期在VTI介质射线走时层析成像研究中提出的分步反演策略实现了多参数联合全波形反演.理论模型实验表明,本文提出的基于Born敏感核函数的各向异性矩阵分解全波形反演方法可以获得较好的多参数反演结果.     

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

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

地壳分层各向异性介质接收函数及其粒子群反演

地壳不同深度介质的地震各向异性是研究地壳不同深度范围变形方式的重要依据.鉴于地壳介质的复杂性,如何从远震体波接收函数中提取不同深度的各向异性参数仍是一个有待深入研究的课题.在已有研究的基础上,本文利用广义反射-透射系数矩阵方法计算的合成地震图,研究了复杂地壳分层各向异性介质的接收函数随反方位角(back azimuth)变化及不同层位各向异性参数对接收函数波场的影响,为各向异性介质接收函数的解释提供了新的理论依据.通过引入粒子群优化理论,发展了分层各向异性介质接收函数全局反演算法.数值及观测数据的验证结果表明,在各向同性速度模型确定的前提下,我们的方法能够可靠地提取地壳分层各向异性参数;在反演中引入曲波变换去噪技术,对于正确解析不同层位的各向异性参数具有重要价值.     

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

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

多组垂直裂缝诱导各向异性的HTI近似及其特征分析

在裂缝诱导各向异性理论研究中通常使用等效HTI介质来近似多组裂缝所引起的综合效应.由于构造运动的复杂性,多组裂缝普遍存在于地壳与油气储层中.为了研究多组裂缝的地震属性特征,分析常用的等效HTI模型对于多组裂缝近似精度及附加裂缝对介质属性特征的影响,本文利用线性滑移模型进行了多组垂直裂缝的单斜各向异性等效介质理论计算,并利用空间搜索方法求取与其最为接近的HTI介质各向异性弹性参数.重点研究了在两种各向异性介质中纵波速度、快慢横波速度和极化特征及其差异,量化分析附加裂缝对于地震属性如速度、极化方向和走时等的影响,研究对附加裂缝敏感的地震属性.此研究结果和方法为进一步研究多组裂缝的反演及识别方法提供基础,同时对于将高阶对称性各向异性介质中已存在的计算方法应用于低阶对称性时的适用程度、精度分析及相关方法研究具有重要作用.     

含流体裂缝介质中地震波场数值模拟

油气勘探开发实践证明,裂缝常常是油气藏存储的空间或运移的通道,因此,裂缝各向异性介质中地震波场的研究越来越倍受关注,国内外很多岩石物理学者、地球物理专家等对裂缝信息的描述提出了很多理论认识与方法技术.本文根据Eshelby-Cheng各向异性裂缝介质模型理论,求取各向异性裂缝介质的弹性参数,并建立Eshelby-Cheng各向异性裂缝介质的波动方程,利用时间错格伪谱法对含流体裂缝介质进行数值模拟,模拟结果表明,采用时间错格伪谱法能有效解决各向异性介质的波场传播,利用时间错格有限差分算子替代普通的差分算子来求解时间导数,利用快速傅氏变换求解空间导数,大大提高了正演模拟的计算精度与计算效率.并且与各向同性介质相比,地震波在含流体裂缝各向异性介质中的传播要复杂得多,各向同性介质层中的波是纯的,其横波不会发生分裂,而在各向异性介质层中,横波将发生分裂.     

永登地震过程中各向异性介质参数演化特征

运用合成地震图方法研究了1995年7月22日永登Ms5.8地震过程中各向异性介质参数动态演化特征.研究结果表明,各向异性介质参数的动态演化特征可以描述地震的孕震过程;各向异性介质参数的微观变化反映了孕震区域应力场的宏观变化;通过研究各向异性介质参数的变化有可能监测地震孕震过程.     

水平单轴应力与横波反射系数

地球物理勘探中,地下应力的估测是非常重要的.本文利用应力诱导的弱各向异性的各向异性参数,借助已有的具有水平对称轴横向各向同性介质顾及高阶小量的横波反射系数近似公式,针对下介质为各向同性受单轴应力作用情况随应力的增加求取应力精度降低的问题,从横波反射系数近似公式的推导出发,联合SH波与SV波反射系数给出新的改进的单轴应力与横波反射系数的直观表达式;计算并对比了受单轴应力作用的各向同性介质的横波近似和严格的反射系数,并利用理论得到的横波反射系数对应力进行了估值.结果表明,利用改进后横波反射系数与应力的直观表达式,当应力增加时,应力估值精度有了明显的改善,这为利用横波反射数据直观估测应力提供理论依据.     

Research on anisotropic parameters by syn-thetic seismogram

Introduction The anisotropy of the Earth crust medium is a common phenomenon(Crampin,1984).More and more observation materials and study results have proved that on most of the Earth′s upper crust is spread fluid-filled cracks and microcracks,which are aligned according to the contempo-rary stress-field.Such distributions of aligned cracks show effectively anisotropy to seismic wave and the phenomenon is called extensive-dilatancy anisotropy(EDA)(Crampin et al,1984).At the same time,it is …     

Seismogram Analysis of Indonesian Earthquakes at DAV Observation Station

The S-wave velocity across the earth structure under Indonesia for Indonesia earthquakes has been investigated through seismogram analysis, simultaneously in the time domain and three Cartesian components. The data were recorded at DAV observational station, the Philippines. The main data set is the seismogram comparison between the measured and synthetic seismogram, instead of travel time data, as commonly used in other seismological research. The synthetic seismogram is calculated using the GEMINI method, which is equivalent to Mode Summation. The above seismogram comparison shows that the global earth mantle of PREMAN gives a deviating synthetic seismogram and has earlier arrival times than those of the measurement. The gradient of β_h in the upper mantle layers is altered into a positive, rather than negative slope as stated in the PREMAN model, and negative corrections are imposed to the zero order of the polynomials coefficients in all earth mantle layers. The excellent fitting, as well as travel time or waveform, is obtained from the surface waves of Love and Rayleigh, surface wave to the S and SS mantle waves as well as the core reflected waves. This result expresses that part of the earth mantle, due to a collision between India and Asia tectonic released zones, has a negative anomaly in S-wave velocity and vertical anisotropy in all of the earth mantle layers.     

Seismogram Analysis and Fitting of Three Earthquakes,West Papua,Recorded at Stations ENH and KMI,China

The S wave velocity structure of the earth below Eastern Southeast Asia has been investigated by analyzing the seismogram from surface wave to multiple depth waves in the time domain and three Cartesian components simultaneously. The wave passes across the front area of subduction zone between the Philippine plate and the Asian plate. The main data are waveform comparisons, instead of the arrival times. The synthetic seismogram is calculated using the GEMINI method. The synthetic seismogram constructed by PREMAN global earth model deviates greatly from the measured one. To solve this problem, corrections are needed for the β speed structure. Corrections cover the gradient change of βh, which turns from negative to positive in upper mantle layers as in the PREMAN, change of earth crust depth and change of zero order coefficients of β velocity function in all earth mantle layers. So, the fitting is obtained, as well as the arrival time or the waveform of Love and Rayleigh surface waves, the S wave and the repetitive depth waves ScS2 and ScS3. This result reveals that the Southeast Asia, being stretched due to tectonic release, has a mantle in some parts with negative anomaly of S wave velocity and vertical anisotropy in all earth mantle layers.     

Teleseismic receiver function using stacking and smoothing of multi seismic-records at a single station

Receiver functions (RFs) obtained using teleseismic wave records at a seismic station and synthetic seismograms indicate that RF with a single teleseismic wave record is related to the selection of record section and to the calculating parameters of the RF. The scatter noise contained in the seismogram also affects the quality of RF. A new method for calculating receiver function, stacking and smoothing multi-seismic records in a single station, is presented in this paper. The RF results using some records and some synthetic seismograms with different noises indicate prominent mantle discontinuity and thus prove that the method is effective and satisfied.     

Seismogram Analysis of Earthquake C060394F,South Java: S-wave Velocity Structure

The S wave velocity structure between the hypocenter of C060394F earthquake,South Java and a series of observatory stations located in Australia and South-East Asia have been investigated through seismogram analysis in the time domain and the three Cartesian components.The synthetic seismogram is constructed from the PREMAN global earth model.Seismogram comparison between the measured and synthetic seismograms shows large discrepancies.A correction to the S wave velocity structure is needed to solve these d...     

S-wave splitting intensity analysis and inversion

Analysing S-wave splitting has become a routine step in processing multicomponent data. Typically, this analysis leads to determining the principal directions of a transversely isotropic medium with a horizontal symmetry axis, which is assumed to be responsible for azimuthal anisotropy, and to the time delays between the fast and slow S-waves. These parameters are commonly estimated layer-by-layer from the top. Errors in layer stripping occurring in shallow layers might propagate to deeper layers. We propose a method for S-wave splitting analysis and compensation that consists of inverting interval values of splitting intensity to obtain a model of anisotropic parameters that vary with time and/or depth. Splitting intensity is a robust attribute with respect to structural variations and is commutative, which means that it can be summed along a ray (or throughout a sensitivity kernel volume) and can be linearly related to anisotropic perturbations at depth. Therefore, it is possible to estimate anisotropic properties within a geological formation (e.g. the reservoir) by analysing the differences of splitting intensity measured at the top and at the bottom of the layer. This allows us to avoid layer stripping, in particular, for shallow layers where anisotropic parameters are difficult to estimate due to poor coverage, and it makes S-wave splitting analysis simpler to apply. We demonstrate this method on synthetic and real data. Because the splitting intensity attribute shows usefulness in S-wave splitting analysis in transversely isotropic media, we extend the splitting intensity theory to lower symmetry classes. It enables the characterization of tilted transversely isotropic and tilted orthorhombic media, opening new opportunities for anisotropic model building.     

垂向裂隙各向异性煤层地震波响应

煤层中存在的裂隙会导致介质表现为各向异性,本文以HTI型煤层为例,结合各向异性介质弹性矩阵和各向异性裂隙理论,推导出不同充填物的垂直裂隙中各向异性参数表达式,将其应用于地震波响应分析;通过改进的交错网格差分法和各向异性Christoffel方程波场分解法,得到地震波合成记录和分解后的P波和SV波记录;将Thomsen群速度与相速度公式,经过坐标轴旋转变换,得到HTI型煤层中不同各向异性参数的地震波速度响应表达式;建立不同类型煤层地质模型,分析了裂隙密度、裂隙充填物以及煤层厚度等参数变化时的地震波响应特征.研究结果为分析垂向裂隙各向异性薄煤层地震波传播规律提供工具,为选用相应地震数据进行地震波各向异性参数反演提供依据.     

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

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

分时窗提取地震子波及在合成地震记录中的应用

提出了利用地震和测井资料精确提取井旁地震子波的分时窗提取地震子波方法，将此方法用于合成地震记录的制作，提高了合成地震记录与地震剖面的吻合度和分辨率，文中详细介绍了该方法的具体实现步骤，并给出了模型处理分析和实例分析。     

Elastic wave propagation in inhomogeneous anisotropic media

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｍｅｄｉａｉｎｅａｒｔｈａｒｅｑｕｉｔｅｃｏｍｐｌｅｘ．Ｔｈｅｒｅｅｘｉｓｔｓｅｖｅｒａｌｕｎｃｏｎｔｉｎｕｏｕｓｐｌａｉｎｓ．Ｎｏｒｍａｌｙ,ｉｔｉｓｕｓｅｄｔｏａｐｐｒｏｘｉｍａｔｅｒｅａｌｍｅｄｉｕｍｗｉｔｈｌａｙ...