基于双相介质弹性波波场分离的逆时偏移方法

随着地震勘探的发展,以单相介质模型为基础的声波逆时偏移和弹性波逆时偏移很难比较精确地描述地下介质和满足复杂的工程要求,尤其是在处理地下含流体多孔介质时,所以需要引入双相介质模型。首先基于Biot理论双相介质一阶速度-应力弹性波方程,利用散度场和旋度场的波场分离理论将固相和流相中的纵横波场解耦分离,并由波场快照图得出固相与流相中的纵波波场与横波波场均完成分离,而快纵波与慢纵波相互耦合均存在于纵波波场中。采用高阶交错网格有限差分方法,结合PML吸收边界条件,推导了具有波场分离的双相介质弹性波正传与反传的偏移延拓算子,结合互相关成像条件,可以得到双相介质弹性波固相与流相中的各分量偏移成像剖面。最后对模型使用该双相介质弹性波波场分离的逆时偏移方法进行偏移成像试算并取得成功,效果较常规双相介质弹性波逆时偏移成像效果有了明显提升。     

基于修正伪Helmholtz分解算子的横向各向同性介质纵横波波场分离

纵横波波场分离是弹性波逆时偏移成像中的关键步骤,可以有效消除纵横波串扰和成像畸变,提高成像精度。在各向异性介质中,常用的波场分离法有非平稳滤波器或者低秩近似的方法,但由于使用多次傅里叶变换,导致波场分离计算成本较高。借鉴基于Helmhotz算子构建各向同性介质纵横波解耦方程的思路,提出消除振幅畸变的修正伪Helmholtz分解算子,通过待定系数方法求解横向各向同性介质中P波和S波2种模式的表达式,并将其变换为一阶解耦的伪弹性波方程,实现时空域纵横波波场解耦。通过简单模型的波场分离测试,得到分离后的P波与S波波场,验证了该方法的波场分离的有效性。此外,将解耦得到的矢量P波与矢量S波应用到弹性逆时偏移之中,利用矢量点积互相关成像条件得到清晰的弹性波逆时偏移成像结果,说明本方法在复杂介质中有着较好的适用性,同时可以有效地应用到VTI介质弹性逆时偏移之中。     

基于源波场重建与波场分解的逆时偏移

逆时偏移以其高精度、适于强变速地质体和不受反射层倾角限制的特点,已经成为处理盐丘、高陡倾角反射层和推覆构造等复杂地质体的标准成像手段。然而,其成像条件需要具有同一时间上的源波场与检波点波场,这在大规模数据上会导致存储的不足与I /O 时间过长。这种成像条件也会在波场传播路径上产生假象。鉴于此,提出了一套可以有效解决这些问题的逆时偏移方案。对于存储与I / O 问题,用源波场重建技术来使源波场与检波点波场在时间上同步,即在源波场正演过程中,仅需存储成像区域边界的波场和最后两层时间的波场作为源波场重建的边界条件和初始条件,从而大幅度节约存储与I /O。对于假象问题,利用上下左右行波分解成像条件与Laplace 滤波技术结合,以达到去除假象的目的。经测试,该方案取得较好的成像效果。其意义在于改进了逆时偏移去假象的方案,使其可以有效地对复杂地质情况进行高质量的成像,并通过源波场重建技术来减少实现逆时偏移所需要的存储与I /O 时间。     

基于矢量波场分离的TI介质逆时偏移

基于两步法分离各向异性矢量波波场,第一步,将各向异性波动方程进行相似变换(使其具有明确直观的运动学和动力学特征),使其投影在各向同性方向上;第二步,各向异性介质传播算子波场分离基础上的波形校正,将投影于各向同性方向上的波场校正回各向异性方向上,最终达到各向异性波场分离的结果。笔者首先对VTI介质进行试算,以证明其可行性;此外,将地表数据进行分解,即所采集到的多分量地震数据在进行偏移成像之前先进行数据分解,得到纯净的qP波数据、qSV波数据以及qSH波数据,通过地表数据的分解将成像过程变成一个可控过程,通过对地表数据输入的控制得到我们所希望的成像结果分量,可以减少无谓的计算;对Hess模型正演炮记录进行了数据分解,并对其进行了逆时偏移成像以验证理论的可行性。     

基于行波分离的VTI介质逆时偏移

现阶段,各向同性介质假设已经不能满足复杂地质构造条件下的成像需求,忽略介质各向异性影响会使速度估计不准确,成像精度降低.为此,笔者从VTI介质一阶速度一应力方程出发,利用Poynting矢量特性对地震波进行方向行波分离,得出了基于上下行波分离的互相关成像条件,并将其应用在复杂VTI介质中.通过层状模型和SEG/EAGE岩丘模型测试验证了该算法的有效性以及在成像上的优越性;另外,当偏移速度不准确时,本文方法仍能够有效压制偏移噪声,成像质量更佳.     

弹性波叠前逆时偏移

利用弹性波一阶速度—应力方程,采用交错网格高阶有限差分法进行震源波场的模拟和接收波场的逆时外推,并结合了能流密度矢量的互相关成像条件,实现了弹性波叠前逆时偏移,且能很好的消除叠前逆时偏移中的假象干扰.笔者通过一个二层模型和一个金属矿模型进行了叠前逆时偏移成像,结果证明所述方法可以有效地消除假象干扰并能对复杂陡倾角界面进行清晰地成像.     

VTI介质中井间地震三分量数据的波场特征

用波动方程的有限差分数值解法,模拟了井间地震观测到的VTI介质中传播的波场,合成井间地震共炮点的三分量记录,对比识别了野外实际井间地震三分量观测记录上的快纵波、慢纵波、快横波和慢横波,分析了各类波的传播速度和偏振特性,为井间地震数据的波场分离、不同类型波场的成像以及纵横波和快慢波的联合解释提供依据。      

基于互相关成像条件的VTI介质保幅高斯束逆时偏移

保幅高斯束逆时偏移是一种兼具高斯束偏移灵活、高效和逆时偏移高精度特性的深度域成像算法,不仅能够消除偏移过程中由于入射角变化和几何扩散引起的振幅误差,而且具有面向目标成像的能力.其以Kirchhoff积分为基础,采用高斯束的加权积分来构建格林函数,进而实现正反向波场的延拓,最后通过相应的成像条件进行成像.这里从基于相速度...     

逆时偏移双向照明波场分离成像条件

通过对波场分离互相关成像条件与炮-检点双向照明互相关成像条件进行研究,将两种方法在噪声消除和补偿成像振幅畸变方面的优势进行结合,提高成像分辨率,改进现有成像条件。其主要步骤是将震源和检波点波场分离成它们的单程波传播分量,采用双向照明补偿互相关成像条件对需要的波场进行成像。通过对模型试算,验证了笔者提出的成像条件能有效的压制噪音,补偿成像振幅。     

旋转交错网格 TTI 介质波场模拟与波场分解

为了得到 TTI 介质中传播的纯纵波和纯横波波场,采用旋转交错网格高阶有限差分法对 TTI 介质中传播的地震波场进行了数值模拟,并利用 qP 波、qSV 波、SH 波的极化方向相互垂直的特性,来实现 qP 波、qSV 波、SH 波的分解。通过对均匀介质模型和复杂层状介质模型分解前后波场快照和地震记录的对比,验证了本文方法的正确性。     

VTI介质中的地层弹性参数反演

在VTI介质中，用非双曲线可以很好地近似P－P波中长排列情况下的旅行时；而对P－SV波来说，应用非双曲线却不能很好地近似P－SV波的施行时。对P－P波、P－SV波的旅行时进行三项泰勒展开时，它们的系数却包含了地震波的各向异性信息及其垂直速度信息，这为地震波的弹性参数反演提供了基础。实验分析表明，在中等排列长度下可以精确地得到了P－P波，P－SV波旅行时的三项泰勒展式系数，从而可以精确地进行地层弹性参数反演。     

基于EMD的云南东川蒋家沟泥石流年输沙量多时间尺度分析

基于经验模态分解(EMD)方法,对云南东川蒋家沟泥石流1965年至2004年40年的年输沙量时间序列进行多时间尺度分析,该时间序列分解成3个本征模函数(IMF)分量和趋势项(Res)分量。分析表明：蒋家沟泥石流的年输沙量存在多尺度的周期性波动,分别是2~4 a、6~10 a和17~21 a的准周期。通过对各IMF分量和Res趋势项分析发现,自1985年之后,波动的振幅开始变大,表明蒋家沟泥石流的年输沙量变化很大。降雨对蒋家沟泥石流的年输沙量及周期性波动有一定程度的影响。人类活动改变环境的同时,也为蒋家沟泥石流提供固体物质来源,这使得蒋家沟泥石流的年输沙量周期性波动变得更加复杂。