多次透射公式实用形式稳定性分析

本文以二维出平面波动方程为例,利用数值方法研究了多次透射公式的实用形式的稳定问题,结果表明：多次透射公式的实用形式与波动方程的集中质量显式有限元结合,不但存在零频漂移,而且同样存在高频振荡失稳,而与中心差分方法结合仅存在零频漂移失稳。     

多次透射公式的一种稳定实现措施

基于对双曲偏微分议程数值解稳定性的GKS（Gustafsson,Kreiss and Sundstrom）准则的物理解释，本文提出了在波动模拟中多次透射公式飘移失稳的一个简单措施，并用波源问题和散射问题的数值试验检验了建议的稳定措施的有效性。     

海上多次波的联合衰减法

地震勘探尤其是海上地震勘探中存在着各类多次波,由于多次波的存在严重影响了速度分析、叠加、偏移成像等地震资料处理,海上多次波主要有全程多次波和层间多次两大类,为了压制海上不同类型的多次波,本文首先分析波场外推、预测反褶积和拉东变换衰减不同多次波的理论基础,然后联合采用这些方法,分别衰减全程多次波和层间多次波;即：首先对炮记录或者接收点记录进行波场外推,建立海底多次波模型,预测并减去全程多次波,然后利用预测反褶积衰减掉周期性明显的多次波,最后将数据转换到τ-p域,用拉冬变换根据同一时间多次波和有效波在速度等方面的差异,进一步分离层间多次波和剩余的全程多次波,并在该域中切除分离出的多次波,从而实现联合多次波衰减处理.通过对悉尼海区和里海等实际地震资料的处理证明,文提出的联合多次波衰减方法在海洋地震资料的处理中有着广泛的应用,联合衰减多次波处理流程具有快速、简洁、易于实现的特点,经处理后的地震资料有效信号损失小、保福性好的特点,有利于速度分析、叠加、叠前偏移等的后续处理工作.     

多次透射公式的一种高频失稳机制

本文利用一维双曲型偏微分方程组初边值问题数值稳定性的GKS定理的物理解释和推广，分析了多次透射公式在数值实现中可能出现的一种高频失稳机制，即在多维离散网格，沿某一空间方向的外行简谐波，由于与其他空间方向节点运动的耦合效应可以使能量传播方向反向。文中通过数值试验对这一失稳机制作了初步验证。     

线性拉东域预测反褶积在海洋多次波去除中的应用

海洋多次波是海洋地震资料处理中最难去除的噪声,目前时空域预测反褶积、SRME、高精度拉冬切除等方法均能有效衰减部分多次波,但对于强反射界面(如硬海底,海底以下高速层等)产生的多次波,单独采用以上方法仍会产生能量较强的残余多次波,针对该问题,根据多次波与一次波周期性差异在线性拉东域比时空域更为明显且有利于压制长周期多次波的特点,可以将线性拉东域预测反褶积技术作为SRME、高精度拉东切除技术的补充,取长补短,联合应用消除深水海底多次波.本文对线性拉东域预测反褶积的基本原理进行了阐述,讨论了如何选取参数,并对其应用范围进行了分析.将线性拉东域预测反褶积应用于南海某凹陷海洋资料的多次波衰减处理中,有效的消除了残余海底多次波,验证了该方法的实用性.     

混合域高分辨率抛物Radon变换及在衰减多次波中的应用

高分辨率Radon变换存在计算效率和分辨率不能兼得的困境.时间域算法可以获得很高的分辨率,但计算效率非常低;频率域算法具有良好计算效率,但分辨率不理想.为此发展了混合域高分辨率抛物Radon变换,即对频率域抛物Radon变换引入时变的稀疏权.本文给出了一种新的混合域高分辨率抛物Radon变换实现方法,并将该算法应用于叠前数据衰减多次波.文中给出了Radon变换和衰减多次波的流程.理论和实际数据算例表明本文方法既有较高的分辨率又有很高的计算效率.     

基于混合有限元格式的完美匹配层与多次透射公式人工边界比较研究

介绍了完美匹配层(PML)人工边界可以吸收不同频率和任意角度入射波的原理以及PML人工边界的构造方法. 在此基础上,将PML人工边界应用于地震波动数值模拟的速度应力混合有限元格式中,探讨了PML应用的可行性,并通过数值试验研究了PML人工边界的反射率,比较了PML人工边界与多次透射公式(MTF)人工边界应用于体波和面波模拟中数值反射的差异,对两种边界的透射效果进行了分析. 结果表明, 尽管数值离散后PML人工边界不再保持完美匹配特性,但PML人工边界在近场波动数值模拟中可获得比MTF人工边界更为理想的吸收效果,在角点透射、大角度掠射情形下尤为明显;PML人工边界在混合有限元格式的数值算法中,未见失稳等不良反应,比MTF人工边界有更好的稳定性;在合理选择参数的情况下,PML人工边界的运算量可接受.      

表面多次波衰减的研究现状与进展

表面多次波是一种相干干扰,在海洋地球物理勘探中表面多次波尤为明显.本文给出了表面多次波衰减方法的研究现状的系统阐述,并分析了表面多次波衰减研究的发展趋势.表面多次波衰减主要有滤波法和波动理论方法.对于滤波法给出了广泛应用的抛物Radon变换法的实例,指出了抛物Radon 变换法存在滤波时能量泄漏的问题.同时本文分析了基于波动理论的处理优势和存在的问题,该方法可以处理任意介质的情形,但是对于数据的采集要求比较严格,另外在进行表面多次波衰减处理以前,要进行一系列的预处理,这些预处理的误差直接决定着多次波衰减的效果.最后,本文给出了表面多次波衰减的发展趋势——将滤波法和波动理论法结合来实现多次波的衰减.两类方法如何有效地结合是该类方法一个重要的研究方向.     

在成像空间中衰减多次波方法研究

在偏移后成像空间中的共成像点道集中可以对多次波进行衰减,对于给定的偏移速度模型,一次波与多次波在叠前偏移后的共成像点道集中具有不同的动校时差,这样我们就可以使用类似于偏移前衰减多次波的方法将一次波和多次波进行分离.本文在成像空间中应用抛物Radon变换分离多次波和有效波,由于每个共成像点道集都包含了复杂三维波场传播效应,所以本文方法具有处理三维数据和复杂地下构造的能力.相比于SRME以及传统Radon变换衰减多次波方法,本文方法能够在保持较小的计算量的同时,保证了衰减多次波的准确性.通过对模型数据试算和对实际数据的处理验证了本文方法在叠前时间偏移后衰减多次波的能力,并取得了很好的成像效果.     

基于频率漂移法的地质透射雷达波的衰减层析成像

本文介绍一种电磁波衰减层析成像的方法。这种成象方法利用电磁波振幅谱重心频率的下移作为投影函数来重建本征衰减系数（目标函数）的空间分布。作者已经将帝种方法应用于处理地质雷达的预测资料。对于大多数天然介质而言。当电磁波信号的频率不高于过滤频率（介质电导率与电容率之比）时，电磁波的衰减会随着频率的增高而增强。因而，接收到的电磁波信号的振幅谱重心频率将发生向低频方向的飘移。这种下移正比于本征衰减系数沿射线     

Multi—transmitting formula for attenuating waves

The MTF is extended to case of attenuating incident wave by introducing an attenuation coefficient.The reflection coefficients of this modified MTF and MTF are evaluated and compared when an attenuating wave impinges on the boundary,and the results demonstrate that MTF can be used to absorb slightly attenuating waves and the modi-fied MTF is more capable of absorbing heavily attenuating waves than MTF.The accuracy of modified MTF is also tested by numerical examples of fluid saturated porous media.     

用优化聚束滤波方法消除低信噪比地震资料中的多次波

针对常规方法如f k,Radon变换不能很好解决的消除低信噪比地震勘探资料的多次波问题,本文提出了优化聚束滤波方法,采用具有静态权的自适应聚束滤波器,并调整聚束滤波器设计的约束准则,去除低信噪比复杂实际资料中的多次波．实际资料处理结果表明,优化聚束滤波方法改善了最小方差无偏聚束滤波方法不能较好地处理低覆盖低信噪资料处理的局限性,在实际处理中计算量增加不大．     

离散介质中地震波传播的数值模拟

伪谱法应用付里叶变换和有限差分求解波动方程，本利用伪谱法模拟了三维离散介质中地震波的传播。结果表明：(1)该方法的计算精度高；(2)即使在均匀、各向同性的离散介质中，地震波的传播也表现出方位各向异性的特点。     

五大连池火山构造地震波特征的初步研究

本文选择五大连池火山区微震活动的模拟、数字地震记录,研究地震波特征,对火 山构造地震进行地震波型分类,分析地震波动力学、运动学特性,以期对我省五大连池火山 区构造地震波有一个初步认识．     

低频大气声重波的非线性Schr?dinger方程

本文采用多重时空尺度展开的方法,对于不可压缩等温大气,导出了二维情况下低频声重波的非线性Schr?dinger方程,讨论了声重波的非线性性质.在垂直方向上,声重波具有可传播性和幅度随高度增长的特性;在水平方向上,声重波能表现为孤子包络调制的波列.     

各向异性弹性波有限元叠前逆时偏移

利用有限元法和最小走时射线追踪的界面点法,实现了各向异性弹性波的叠前逆时偏移.理论模型资料的偏移结果清晰准确,证实了该系统的可靠性.通过对各向异性模型资料的各向同性偏移处理,发现常规偏移剖面存在较大误差,地质体的垂向深度和横向位置与实际模型有偏离.     

对波达波夫和Pride震电波方程组的对比分析

用Biot介质参数说明了波达波夫震电波方程组中弹性动力学 参数的含义，解释了第一类和第二类震电效应的意义，在忽略第一类震电效应条件下将该方 程组与Pride方程组进行比较，说明了二者在描述第二类震电效应方面的异同点. 同时指出 ：波达波夫方程组忽略了流体与固体的耦合质量；方程中的黏性耗散项丢掉了一个孔隙度因 子，依据该方程组计算出的弹性波和转换电场的幅度将偏大；边界条件之一存在错误，会影 响对波在界面上的反射透射规律的描述.     

各向异性弹性波有限元叠前逆时偏移

利用有限元法和最小走时射线追踪的界面点法，实现了各向异性弹性波的叠前逆时偏移.理论模型资料的偏移结果清晰准确，证实了该系统的可靠性.通过对各向异性模型资料的各向同性偏移处理，发现常规偏移剖面存在较大误差，地质体的垂向深度和横向位置与实际模型有偏离.     

流体饱和多孔介质波动方程的有限元解法

讨论了流体饱和多孔介质中波传播问题的有限元解法,首先在Ｂｉｏｔ理论的基础 上,概述了数学问题的提法,然后提出了一种新型简便的人工边界上的无反射边界条件,同时 给出了有人工边界时流体饱和多孔介质波动方程的有限元计算公式．数值试验的结果表明, 本文提出的无反射边界条件和数值方法均很有效．     

流体饱和多孔介质波动方程的有限元解法

讨论了流体饱和多孔介质中波传播问题的有限元解法,首先在Ｂｉｏｔ理论的基础 上,概述了数学问题的提法,然后提出了一种新型简便的人工边界上的无反射边界条件,同时 给出了有人工边界时流体饱和多孔介质波动方程的有限元计算公式．数值试验的结果表明, 本文提出的无反射边界条件和数值方法均很有效．