一种TI介质纯qP波正演方法及其在逆时偏移中的应用

基于Tsvankin提出的精确频散关系,利用近似展开的方法,推导出解耦合的TTI介质纯qP波近似方程,并将方程中的偏微分算子分解成一个laplace算子和一个标量算子,用于代表qP波的精确传播方向,构建时间域二阶纯qP波方程.此推导过程无需设置横波速度为零,能够更加精确地描述qP波的运动学特征.这个方程相比于求解波数域二阶解耦qP波方程,计算效率高,存储需求小;相比于基于Alkhalifah频散关系推导的时间域二阶纯qP波方程,假象干扰压制好,数值误差小,更具一般性.但此方法求解波矢量时采用波场梯度一阶渐近近似,会造成垂直于对称轴方向的波场振幅不准确.为了较正振幅,将椭圆分解方法应用于此方程中,构建纯qP波椭圆分解方程,使得振幅更加均衡,并与Xu等提出的方程比较分析,应用本文构建的纯qP波椭圆分解方程得到的波场振幅值更加准确.本文首先选取了均匀TI介质模型进行了qP波正演模拟,并抽取波场单道波形进行振幅分析,验证了本文构建的纯qP波方程和纯qP波椭圆分解方程的正确性及有效性;然后选取BP TTI模型进行了qP波正演模拟,将其qP波正演结果和均匀TI介质模型振幅分析结果相结合,突出了本文构建的纯qP波椭圆分解方程的优势及适应性;最后选取逆冲模型和BPTTI模型,应用本文构建的纯qP波椭圆分解方程对其进行逆时偏移成像,验证了本文构建的纯qP波椭圆分解方程在逆时偏移中的可行性和适用性.     

基于新的惩罚因子算法的波场重构反演

全波形反演是一种高精度的反演方法,其目标函数是一个强非线性函数,易受局部极值影响,而且反演过程计算量较大.波场重构反演是近几年提出的一种改进的全波形反演理论.该反演方法通过将波动方程作为惩罚项引入到目标函数中,通过拓宽解的寻找空间减弱了局部极小值的影响,而且反演过程不需要计算伴随波场,提高了计算效率.但该反演方法一直缺少准确的惩罚因子算法,直接影响到该方法的准确度.本文将波场重构反演拓展到时间域并利用梯度法进行波场重构.频率域的惩罚因子用来加强波动方程的约束,而时间域惩罚因子表现为调节模拟波场和实际波场的权重因子.为此,我们根据约束优化理论,在波动方程准确以及重构波场与反演参数解耦的假设下,提出以波动方程为目标函数的新的惩罚因子算法.根据波形反演在应用时普遍存在的噪音干扰、子波错误和低频信息缺失的情况下,应用部分Sigsbee2A模型合成数据对本文提出的算法进行实验.数值实验结果表明：基于新的惩罚因子算法,在其他信息不准确的情况下,波场重构反演可以给出高精度的反演结果.     

基于射线追踪的成像域走时层析反演及其应用

层析成像是现今地震勘探开发处理中的常用手段,针对传统走时层析反演中角度覆盖问题,利用正则化约束加以解决,利用角道集的拾取拟合来判断实验模型的准确性,同时,不再利用原始炮记录拾取剩余时差,利用深度残差与走时残差的线性关系提高走时残差的拾取精度。以三层洼陷模型和实际资料进行试算,结果表明:本方法可以较好地反演更新偏移速度场,得到更精确的偏移成像结果,对于小尺度地质构造可以精细刻画。      

实际单井微测井资料Q值估计方法的改进

近地表对地震波有强烈的吸收衰减作用,严重影响地震数据的分辨率,利用微测井资料对近地表吸收结构进行估算,可以为高分辨率地震勘探提供技术支撑.采用频谱比等主流方法对实际单井微测井数据进行Q估计时,观测到的微测井数据需要满足微测井Q估计的假设条件和Q衰减理论模型,但是由于采集因素的影响,实际微测井数据会存在“限幅”、初至时间异常、激发子波不一致等问题,造成单井微测井Q估计不稳定、不可靠,甚至会估计出负的Q值.本文在梳理实际单井微测井数据存在问题的基础上,对这几个问题采用针对性策略进行了处理,缓解了实际单井微测井数据的缺陷对Q估计带来的影响,实际数据的Q估计结果验证了本文方法的可靠性和有效性.     

基于双变网格算法的地震波正演模拟

为了适应对局部复杂模型的精细模拟,本文实现了可变网格算法,对速度场进行局部加密,从空间上有效地提高模拟精度同时又降低计算机内存需求.但是在数值模拟中,由于稳定性条件的限制,当空间网格变化时,时间稳定性仍然必须满足最短波长的原则,从而增加了时间汁算量.为了配合空间可变网格技术,本文对时间层计算也进行了局部变化,提出了双变...     

弹性矢量波成像域联合层析反演

在地震弹性矢量波场框架下,推导了多波联合层析速度反演方程以及走时残差与角道集剩余曲率的转换关系式,提出了一种利用成像域角道集更新P波、S波速度的走时层析反演方法.其实现过程可以概括为:将弹性波多分量数据作为输入,基于高斯束实现矢量波场成像并提取角道集,利用层析反演方程求解慢度更新量,最终获得多波联合反演结果.模型试算和实际资料处理验证了该方法的反演效果,能够为弹性矢量波联合深度偏移提供高质量的叠前速度场.     

近地表速度建模研究现状及发展趋势

速度是反映地下构造和岩石物性的一个重要参数。近地表速度的精度直接影响勘探区域的地震资料静校正、速度分析以及最终成像的效果。目前常用的近地表建模方法和技术多基于高频近似的射线理论,不能满足当前近地表精细建模的需求。通过对微测井、折射波法、面波法、走时层析以及全波形反演等近地表建模技术的全面调研,总结了它们的适应性、优缺点及研究应用现状,指出联合走时层析与波形反演的技术方法在时间域分步骤、多尺度的反演策略是近地表高精度建模的有效手段和发展方向,能有效提高建模精度,适应高精度成像要求。该方法在近地表的矿产普查、工程物探、油气勘探等领域存在广泛的应用前景。     

地震资料多次波压制方法研究进展及展望

在反射波地震勘探中,发育的多次波会严重降低地震资料处理的精度,进而误导解释人员给出错误的解释结果.因此,多次波压制方法研究一直以来都是地震资料处理中不可忽视的重要环节.目前,将多次波压制方法分为基于特征差异的滤波方法和基于波动理论的预测相减法是一种较为流行的分类方法.前者计算效率较高,但是多次波压制精度有限.后者理论先进,不需要或者极少需要地层先验信息,多次波压制精度较高,但是往往伴随高计算成本.通常,采用单一方法往往难以满足复杂实际资料处理的需求,汇二者所长的综合法则可以针对实际地震资料的特点,通过不同方法组合解决实际应用问题.据上述分类方法,首先详细论述了每一类方法中经典方法的原理、研究进展及发展方向等;然后介绍了目前实际应用中的综合法多次波压制实例;最后,在回顾传统多次波压制方法的研究进程基础上重点介绍了近年来出现的几种新方法,并展望多次波压制的发展前景.     

表层多次波成像方法技术研究

在处理表层多次波比较发育的地震数据资料时,和常规将多次波作为噪音进行压制和去除的思路不同,本文提出利用表层多次波来进行地震构造成像的方法技术流程.首先基于互相关技术将表层多次波转化为准一次波,在利用准一次波进行偏移成像之前,提出应用局部倾斜叠加变换的方法,压制由于其特殊的构建过程而常常伴随准一次波的背景干涉噪音,提高准一次波数据的信噪比;然后应用扩展的分步傅里叶偏移方法对去噪后的准一次波数据进行成像,此偏移成像方法在实现分步傅里叶方法时选取多个背景参考速度,不仅继承常规分步傅里叶偏移方法稳定高效,没有有限差分计算存在的频散影响的特点,而且可以增强常规分步傅里叶偏移方法处理复杂陡倾构造的能力.通过数值模拟数据测试,说明本文提出的这一套利用表层多次波进行复杂地下构造成像的方法技术的可行性和有效性.     

实现复频移完全匹配层吸收边界条件的递推卷积方法

完全匹配层吸收边界（PML）已经被证明是非常有效的边界吸收技术,对体波和面波的吸收都具有非常好的效果,已经被广泛应用于弹性波的数值模拟中。但是在某些情况下传统的PML技术还是存在一定的问题,比如对掠射情况下的体波和窄区域自由表面条件下的面波的吸收等等。在坐标变换中采用复频移拉伸函数的复频移PML可以有效地改善PML边界条件的吸收性能。基于弹性波一阶速度-应力方程,推导了复频移PML的递推卷积实现方法,并采用交错网格高阶有限差分法对其进行了数值模拟,与传统的PML进行了对比。结果表明：传统的PML对掠射情况下的体波和窄区域自由表面条件下的面波吸收不足,会产生虚假反射,影响真实波场;而基于递推卷积的复频移PML算法能够有效地改善困难情况下的吸收效果,并且在实现过程中不用分裂变量,应用更加方便简单。计算卷积时采用递推的形式,推导过程直观易懂,易于编程,而且不会增加计算量,存储量也没有太大的变化。