弹性波场和稳恒电流场的Frechet导数计算方法

在地球物理参数反演中，观测的物理场量，例如波场中的位移或电流场中的电位，对于介质参数的编数的计算是非常重要的。本文采用自伴随算子和格林函数的方法给出了弹性波场和稳恒电流场问题的Ｆｒｅｃｈｅｔ导数解析式。该式表明，Ｆｒｅｃｈｅｔ导数可由格林函数Ｇ及其梯度ΔＧ求得。把任意场源正演问题的解和Ｆｒｅｃｈｅｔ导数的计算统一起来，并提出用有限元方法计算任意背景场下的格林函数，从而得到Ｆｒｅｃｈｅｔ导数的数值解     

黄土塬煤层采空区弹性波多波波场特征

黄土塬已知煤层采空区的弹性波多波波场研究表明,反射波、散射波、天然源面波和人工源面波等各种弹性波波场显示出表征性特点。借助于这种波场特征,对于煤层采空区的地震勘探具有一定的指向意义。     

裂隙各向异性介质2.5维弹性波场数值模拟

研究地震波在各向异性介质中的传播规律是勘探地震学领域的一个重要内容。本文首先讨论了２．５维波动方程，其次，在２．５维空间，用有限元数值模拟方法模拟了裂隙各向异性介质波场，结果用波场切片做了显示。最后对模拟结果做了基本分析，这里讨论的模拟方法也适用于其它类型的得各向异性介。     

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

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

基于解耦传播的波场分解方法在VTI介质弹性波逆时偏移中的应用

由于弹性波逆时偏移更符合实际情况,而且转换波的成像结果有更高的分辨率,因此弹性波逆时偏移越来越受到重视。然而,弹性波逆时偏移需要多波多分量数据,为了减少成像结果中的串扰假象,在逆时偏移过程中进行P波和S波分解就变得非常必要。结合基于向量的激发振幅成像条件,我们把基于解耦传播的波场分解方法应用到弹性波逆时偏移中,并对比了其在各向同性介质和横向各向同性（VTI）介质中的应用效果。结果说明,该方法可以在各向同性介质中完全分解P波和S波,并保留波场的向量信息。虽然在VTI介质中有较小的分解残余,但是该分解残余不会在逆时偏移结果中产生明显的串扰;因此,这种波场分解方法可以应用于各向同性介质和VTI介质弹性波逆时偏移。该方法是在时间空间域实现的,可以在波场传播过程中直接对P波和S波进行分离,应用方便,计算效率高。与利用Helmholtz分解的弹性波逆时偏移相比,该方法避免了在PS波成像结果中的极性反转问题。在复杂Hess VTI模型的逆时偏移结果中,高速岩体和断层的成像清晰,甚至是两个低速薄夹层也能较好成像;这说明该方法对复杂介质具有较好的适应性。PS波成像结果中的各向异性体成像清晰,说明各向异性介质弹性波逆时偏移可以对传统逆时偏移不能很好成像的构造进行成像。     

煤矿巷道内地震勘探波场特征

为了研究巷道内地震波的波场特征,识别出波场中的各种波,利用交错网格高阶有限差分正演方法模拟了煤矿巷道内点震源激发及平面波震源激发产生的弹性波场。通过波场分析得出煤矿巷道内的共炮点道集中除了有常规地面地震中的各种波外,还有在巷道内上下界面的多次反射波和多次折射转换波,是巷道内地震勘探的主要干扰波。利用平面波震源在巷道内激发产生的波场中,巷道下方的空洞在平面波震源激发的地震记录上,其多次反射绕射波互相叠合在一起,表现为“糖葫芦”状的明显异常特征。这些都为巷道内地震勘探资料处理和解释提供了重要依据。     

高阶交错网格有限差分弹性波场模拟的精度分析

交错网格波场数值模拟是目前地震正演中广泛使用的方法,为对比分析不同阶数的差分格式下产生的计算效率和精度差异,重新推导了弹性波方程的4种时间4阶、空间2N阶的差分公式及系数,并计算了他们的稳定性条件。利用这4种差分格式进行弹性波场数值模拟,对比分析了波场快照、合成地震记录及CPU时间。结果表明：时间4阶、空间6+6阶精度的交错网格有限差分方法在进行地震波场数值模拟时具有较高的计算精度和计算效率。     

叠前弹性波逆时深度偏移及波场分离技术探讨

这里探讨了Sun和McMechan提出的叠前弹性波标量逆时深度偏移方法,即在地表附近,对地表接收到的弹性波波场分量进行波场逆时延拓,然后分别对在波场延拓过程中通过波场分离得到的纯纵、横波反射波场进行偏移。在模型中,每个网格点的纵、横波成像时间,为震源到该网格点的纵波初至时间。利用纵波速度模型和声波传播方程的有限差分解,对反射纵波进行逆时延拓和成像,利用横波速度模型和声波传播方程的有限差分解,对转换横波进行逆时延拓和成像,并在偏移前进行极性校正。然后,在Sun和McMechan偏移成像方法的基础上,提出了直接偏移成像方法,即对地表接收到的弹性波波场分量在整个模型范围内进行波场逆时延拓。在波场延拓过程中,对符合成像条件的网格点进行波场分离、偏移成像和转换波极性校正。最后,通过数值模拟对二种偏移成像方法的结果进行了比较。     

各向异性介质储层数值模拟及波场特征研究

各向异性介质储层中波场特征的研究,有利于对多波资料的准确解释。在弹性波方程基础上,结合数值模拟技术与AVO技术,对多种介质模型和第二类含气砂岩储层模型进行数值模拟与波场特征描述。以表征介质各向异性程度的Thomsen参数为出发点,研究各参数对地震波场动力学和运动学特征的影响,并对第二类含气砂岩介质中模拟得到的正演记录,以及各向异性参数对波场的影响程度进行分析,进一步加深了对该类储层各向异性特征的认识,这为各向异性介质中的参数反演与地震资料解释提供了有力的参考依据。     

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

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

地下管线弹性波探测波场的建模与分析

在地下管线弹性波探测中,由位于地下圆孔边界上有限时长的力源激发形成瞬态场,应用波函数展开法对此波场进行建模。以波函数和待定系数组成的级数序列描述探测波场中的势函数,根据界面处应力和位移连续的边界条件求解待定系数,从而得到波场中弹性波的解析解。据此进行数值研究,分析探测波场中弹性波的传播特征,说明了散射波的类型以及能量对比关系,并总结了反射信号幅值与地下管线物理参数之间由声阻抗匹配程度支配的反射信号强度变化规律,而管线半径对反射信号强度的影响要大于物理参数的作用。最后以实验结果验证上述结论。     

井间地震弹性波传播特征数值模拟

采用弹性波方程数值模拟和Zoeppritz方程解相结合的方式,研究了井间地震固-固界面上以及低速和高速夹层中弹性波场特征及传播规律.给出了井间地震弹性波的主要类型,其中包含一种特殊面波,该面波由激发井和接收井内的地层界面处产生的二次源形成,并说明了反射波以广角反射占优势;利用正演模拟记录对广角反射现象进行了分析,得出了在临界角处广角反射波的振幅、相位突变,极性反转,这与反射系数变为复数相对应;提出利用临界反射波特征和临界反射点位置来确定地层速度.震源在高速夹层中激发时,其中传播的弹性波场简单;而在低速夹层中,则弹性波场非常复杂.低速夹层中的导波存在明显的频散现象,其振幅和频率与低速层的速度和厚度密切相关.提出利用导波可以确定低速层的横向变化和连通性.     

涟邵矿区的槽波特征

介绍了涟邵矿区槽波测量概况,讨论了该矿区槽波特征,指出了槽波地震探测技术能解决的问题。     

黏弹介质槽波超前探测数值模拟与波场分析

目前常用槽波超前探测掘进工作面前方隐伏断层,而实际煤层具有黏弹性,对槽波有吸收衰减作用。为研究黏弹煤层介质中槽波超前探测的波场特征和传播规律,基于Kelvin-Voigt一阶速度-应力方程,建立了三维含断层煤系地质模型,进行了三维有限差分数值模拟。结果表明：黏弹介质中反射槽波能量衰减大,其传播与衰减特征比完全弹性介质更符合实际煤层反射槽波传播特征。煤层Q值越大槽波超前探测效果越好,煤层Q值过小时,槽波超前探测效果差,可利用x分量的反射横波与z分量的P-S波超前探测。断层落差小于煤厚时,y和z分量的槽波超前探测效果都较好;落差大于煤厚时,可利用y分量进行槽波超前探测。断层夹角大于60°时,可利用y和z分量进行槽波超前探测;夹角小于60°时,可利用x和y分量进行槽波超前探测。研究结果可为煤矿槽波超前探测提供理论支持。     

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

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

弹性波CT技术及应用

弹性波ＣＴ技术是一种新的物探方法。本文简单地介绍了该技术的方法原理及影响图像重建的因素，并以实例说明了弹性波ＣＴ技术在煤矿井下的应用效果     

应用弹性波法评价岩体质量

近年来弹性波法在水电、交通、铁道等工程勘察中普遍开展,从应用单一的波速参数发展到现在的系列测试技术,尤以弹性波动力学参数的取得,使测试结果更全面地反映了工程岩体的特性。加之本技术具有快速、简便、经济、便于大面积测试等优点,因而受到愈来愈多的工程技术人员的重视和应用。     

基于波场延拓的多次波压制技术及其改进

水层多次波是海上地震勘探过程中的一个重要干扰因素。为压制水层多次波,本文在倾斜模型条件下采用基本的正向延拓方法去除多次波,取得了较好的压制效果。在此基础上实现了弯曲水底情况下正、反向波场延拓方法压制多次波;经对比发现,反向延拓方法克服了复杂水下介质情况下对格林函数求取的干扰,消除多次波效果更好。