全波形反演在缝洞型储层速度建模中的应用

速度是地震偏移成像准确与否的关键所在.全波形反演综合利用地震波场运动学和动力学信息,能够得到相比传统速度建模方法更高频的成分.全波形反演的理论比较成熟,但实际应用成功的例子相对较少,特别是对于陆上地震资料.塔里木盆地地震地质条件复杂,为了实现缝洞型储层的准确成像,本文开展了针对目标靶区的全波形反演精细速度建场研究.采用一种时间域分层多尺度全波形反演流程:首先通过层析成像建立初始速度模型;其次利用折射波反演浅层速度模型;最后利用反射波反演中深层速度模型.偏移成像结果表明基于全波形反演的速度建模技术能有效改善火成岩下伏构造的成像精度,显示了全波形反演在常规陆上采集资料的应用潜力.     

基于反射地震数据的时频域包络反演

包络信号含有丰富的低频分量,即使在地震数据缺失低频条件下,包络目标函数也能有效缓解全波形反演的周期跳跃现象.但是,当初始速度模型较为平滑时,观测数据中的反射地震事件在模拟数据中没有与之相对应的波形,导致包络反演初期无法很好地利用反射波信号进行速度建模.本文提出基于反射地震数据的时频域包络反演方法,通过结合反射波全波形反演理论,构建反射波时频域包络目标函数,来提高包络反演的速度建模精度.本文首先利用Gabor变换获取时频域地震数据,并提取振幅信息,即为时频域包络信号.然后,推导反射波时频域包络反演的伴随震源和梯度算子.Marmousi模型数据测试结果表明,基于反射地震数据的时频域包络反演方法可以为全波形反演提供一个较好的初始速度模型.     

基于反射地震数据的时频域包络反演

包络信号含有丰富的低频分量,即使在地震数据缺失低频条件下,包络目标函数也能有效缓解全波形反演的周期跳跃现象.但是,当初始速度模型较为平滑时,观测数据中的反射地震事件在模拟数据中没有与之相对应的波形,导致包络反演初期无法很好地利用反射波信号进行速度建模.本文提出基于反射地震数据的时频域包络反演方法,通过结合反射波全波形反演理论,构建反射波时频域包络目标函数,来提高包络反演的速度建模精度.本文首先利用Gabor变换获取时频域地震数据,并提取振幅信息,即为时频域包络信号.然后,推导反射波时频域包络反演的伴随震源和梯度算子.Marmousi模型数据测试结果表明,基于反射地震数据的时频域包络反演方法可以为全波形反演提供一个较好的初始速度模型.     

利用全波形反演提高逆时偏移成像的精度

逆时偏移对速度模型有较强的依赖性,速度模型精度高低决定着偏移效果的好坏,速度模型越准确,偏移效果越好.因此,如何获取高精度的速度模型是实现有效逆时偏移的关键步骤.全波形反演方法利用叠前地震波场的运动学和动力学信息重建地下速度结构能够得到高精度的速度模型.为了验证全波形反演能否提高逆时偏移效果,本论文以OVERTHRUST速度模型为例,在频率域进行全波形反演,并实现了速度模型的更新.逐步缩小了初始速度模型与真实速度模型的残差,提高了速度模型的准确度.然后使用初始OVERTHRUST速度模型以及经过反演的速度模型分别进行逆时偏移,应用结果表明,使用经过全波形反演后的速度模型进行逆时偏移,其成像效果较好.因此,逆时偏移结合全波形反演可以有效的提高成像精度.     

VSP直达波和反射波波动方程走时联合反演

地面资料全波形反演采用低波数回折波敏感核和高波数偏移等时线恢复地下模型的长波长和短波长分量.当回折波的穿透深度有限时,很容易陷入局部极值.波动方程反射波走时反演采用反射波的透射敏感核更新速度模型.当浅层速度存在较大误差时,仍无法得到满意的反演结果.VSP资料中直达波是一种透射波,穿透深度大,可用于井旁背景速度场建模.本文发展了针对VSP观测方式的波动方程走时联合反演(直达波和反射波)方法.采用全局优化参数反演法实现VSP上下行波和纵横波分离,构建了基于下行直达波和上行反射波走时残差的混合目标函数,推导了相应的伴随状态方程和梯度公式,给出了背景速度和反射系数分步反演流程.理论和实际VSP资料应用表明：直达波和反射波走时联合反演可以得到运动学特征精确的速度模型,为后续全波形反演提供可靠的初始模型;联合反演比单独直达波/反射波走时反演的精度高,比直达波和反射波走时级联反演的耗时少.     

地震反射波反演二阶优化方法及其应用

常规地震数据大都缺乏低频成分与长炮检距信号,经典的全波形反演方法不易获取中、深层弹性参数模型的长波长分量,波动方程反射波形反演作为替代方法近来受到极大关注.然而,现有的反射波形反演方法几乎都采用梯度类的一阶优化算法,收敛性和精度都有待提高.本文在二阶优化理论框架下,推导弹性参数背景与扰动模型的反射波敏感核、泛函梯度以及海森算子,揭示海森矩阵对泛函梯度的去模糊化作用和改善反演的工作机制.推覆体模型合成数据实验表明,相比于常用的共轭梯度法,利用近似海森矩阵的高斯-牛顿法明显提升了反射波形反演的收敛性与宽谱建模能力.在东海实例中,本文方法超越常用的反射走时层析技术,通过改善中、深层偏移速度建模,支撑逆时偏移高分辨率刻画长江坳陷内部复杂的断裂系统,改善了深部基底的成像质量.     

黏声方程Q值反射波反演

地震波在非弹性介质中的衰减效应常用品质因子Q度量.相对准确的Q模型对提高强衰减介质中地震波成像的质量至关重要.本文提出了黏声介质反射波形反演(QRWI)方法来重建地下宏观Q模型.在缺乏大偏移距和低频地震数据时,该方法以黏声波方程为波场传播引擎,利用反射波核函数对模型中深部的敏感性去提取背景Q值.当速度高、低波数成分均已知时,基于波形拟合的QRWI可以获得较高分辨率的反演结果.由于地下介质速度的高波数扰动很难准确估计,本文通过引入峰值频移目标函数,极大地降低了QRWI对速度高波数成分的依赖.理论合成数据实验结果表明,本文方法反演得到的宏观Q模型可以满足衰减补偿逆时偏移成像的要求.     

起伏地形下的高精度反射波走时层析成像方法

全球造山带及中国大陆中西部普遍具有强烈起伏的地形条件.复杂地形条件下的地壳结构成像问题像一面旗帜引领了当前矿产资源勘探和地球动力学研究的一个重要方向.深地震测深记录中反射波的有效探测深度可达全地壳乃至上地幔顶部,而初至波通常仅能探测上地壳浅部.为克服和弥补初至波探测深度的不足,本文基于前人对复杂地形条件下初至波成像的已有研究成果,采用数学变换手段将笛卡尔坐标系的不规则模型映射到曲线坐标系的规则模型,并将快速扫描方法与分区多步技术相结合,发展了反射波走时计算和射线追踪的方法.进而利用反射波走时反演,实现起伏地形下高精度的速度结构成像,从而为起伏地形下利用反射波数据高精度重建全地壳速度结构提供了一种全新方案.数值算例从正演计算精度、反演中初始模型依赖性、反演精度、纵横向分辨率以及抗噪性等方面验证了算法的正确性和可靠性.     

基于Hinge损失函数的垂向全变差约束全波形反演

全波形反演可以为叠前深度偏移成像提供更高精度的速度模型,但该方法具有较强的非线性,对初始速度模型的依赖性较强,尤其是在实际应用中,地质条件复杂多变,速度变化不连续,增加了反演非线性程度,常常使反演陷入局部极小值,影响反演的精度.全变差约束在图像去噪领域应用广泛,属于非光滑约束,在去噪过程中能有效的保留图像的不连续界面和边缘信息.本文提出基于Hinge损失函数的垂向全变差约束全波形反演方法,在全变差约束的基础上,利用Hinge损失函数控制模型的更新方向,并使用原-对偶混合梯度算法进行求解,给出这一优化问题的迭代格式,有效提高了对地下不连续界面的重构精度,同时也降低反演对初始速度模型的依赖程度.数值算例证明：与常规全波形反演方法相比,基于全变差约束的全波形反演方法可以有效的重构速度模型中的不连续界面,尤其对高速体边缘的重构效果更明显,但该方法对初始速度模型的依赖性仍然较强;基于Hinge损失函数的垂向全变差约束全波形反演方法降低了对初始速度模型的依赖程度,可以从一个较差的初始模型通过循环迭代的方式最终得到同样精确的速度模型,较好的重构了高速体边缘和不连续界面.     

跨孔雷达全波形反演成像方法的研究

跨孔雷达全波形反演是一种使用全波形信息反演两钻孔之间地下信息的层析成像技术.常规的层析成像反演大部分采用射线追踪方法,其中基于初至时的射线追踪方法可以反演出速度剖面（介电常数剖面）,基于最大振幅的层析成像可以反演出衰减剖面（电导率剖面）.常规射线追踪方法有许多不足,究其原因是该方法仅使用了小部分的信号信息.为了进一步提高成像分辨率,本文全面推导了全波形跨孔雷达层析成像反演方法,该方法利用雷达波全幅度相位信息能够反演出地下高分辨率的介电常数和电导率图像.本文通过基于局域网的分布式并行算法,有效地解决了巨量数据正演计算问题.文中首先建立了基于单轴各向异性介质完全匹配层的时间域有限差分二维正演算法,进而通过应用包括时间维度在内的全波场信息与残场逆向传播的全波场信息乘积来计算梯度方向,通过求取以步长为自变量的目标函数的极值确定步长公式,并提出以第一次介电常数反演作为同步反演的初始模型,能够有效提高收敛速度.本文对多组模型进行成像实验,取得了较好的反演效果.     

地震全波形反演及其探测壳-幔结构的研究进展

地球介质的物理属性(如速度)隐含了地球内部结构和演化信息,因而获取高精度的地下介质模型是探测内部结构的重要目标之一.全波形反演方法利用地震观测记录的振幅、相位等全波信息来重建速度和密度等介质参数模型,被认为是建立高精度地下介质模型的最有潜力的研究方法之一,其不但是石油勘探地震成像的主要研究方向,而且在结构地震学获得了较...     

一种基于拟相位信息目标泛函的反射波波形反演方法（英文）

反射波波形反演(RFWI)通过交替更新模型的中低波数和高波数成分,为传统的全波形反演(FWI)提供较为准确的初始速度场。但是由于低频信息的缺失和振幅信息的复杂性,传统的反射波波形反演存在极强的非线性。将振幅与相位信息的分离能够有效提高RFWI的线性程度,然而传统的相位计算方法具有严重的相位跳变(wrapping)问题。基于此,本文首先给出一种拟相位信息的计算方法,利用地震数据的包络信息,以避免传统相位信息的跳变问题。然后利用拟相位信息构建反射波波形反演的目标泛函,并给出相应的伴随震源和梯度表达式。理论分析证明了拟相位信息的计算方法能够保证相位信息的稳定和准确性,并且基于拟相位信息的目标泛函具有更广的全局收敛域。将本文方法应用到部分Sigsbee2A模型并与传统的目标泛函对比结果表明:与传统的全波形反演方法相结合,新的RFWI能够在初始模型不准确情况下得到较为准确的反演结果。除此之外,本文方法对含噪音数据和低频数据缺失情况具有一定的适用性。     

基于振幅加权的拖缆与海底地震资料联合波形反演方法

南黄海新近系沉积层与其下方中-古生代海相地层之间速度差异大,存在强波阻抗界面,导致接收的浅部和中深部地层的反射波振幅差异极大.一般全波形反演方法未考虑不同界面反射波振幅的差异,对南黄海浅部沉积层速度的反演效果好,而对中深部地层速度的反演效果差.为此,本文提出一种基于振幅加权的拖缆与海底地震资料联合波形反演方法.该方法对中深部弱反射波数据施加较大的权重,对浅部强反射波数据施加较小的权重,从而均衡浅、中和深部地层反射波的振幅,使得自浅至深所有地层速度均能得到较好反演效果.同时,该方法利用海底地震资料和拖缆地震记录积分的低频成分反演低波数背景速度场,利用拖缆资料的高频成分提高反演模型的分辨率,充分发挥了拖缆和海底地震资料的优势互补作用,获得了比仅使用单一资料更好的反演结果.根据地震剖面和钻井资料,建立了南黄海中部隆起的速度模型.利用该理论模型的合成地震记录对本文方法进行了测试和分析,验证了本文方法的正确性和有效性.     

数据域特征波反射走时反演

地震波的运动学信息(走时、斜率等)通常用于宏观速度建模.针对走时反演方法,一个基本问题是走时拾取或反射时差的估计.对于成像域反演方法,可以通过成像道集的剩余深度差近似计算反射波时差.在数据域中,反射地震观测数据是有限频带信号,如果不能准确地确定子波的起跳时间,难以精确地确定反射波的到达时间.另一方面,如果缺乏关于模型的先验信息,则很难精确测量自地下同一个反射界面的观测数据同相轴和模拟数据同相轴之间的时差.针对走时定义及时差测量问题,首先从叠前地震数据的稀疏表达出发,利用特征波场分解方法,提取反射子波并估计局部平面波的入射和出射射线参数.进一步,为了实现自动和稳定的走时拾取,用震相的包络极值对应的时间定义反射波的到达时,实现了立体数据中间的自动生成.理论上讲,利用包络极值定义的走时大于真实的反射波走时,除非观测信号具有无限带宽(即delta脉冲).然而,走时反演的目的是估计中-大尺度的背景速度结构,因此走时误差导致的速度误差仍然在可以接受的误差范围内.利用局部化传播算子及特征波聚焦成像条件将特征波数据直接投影到地下虚拟反射点,提出了一种新的反射时差估计方法.既避免了周期跳跃现象以及串层等可能性,又消除了振幅因素对时差测量的影响.最后,在上述工作基础之上,提出了一种基于特征波场分解的新型全自动反射走时反演方法(CWRTI).通过对泛函梯度的线性化近似,并用全变差正则化方法提取梯度的低波数部分,实现了背景速度迭代反演.在理论上,无需长偏移距观测数据或低频信息、对初始模型依赖性低且计算效率高,可以为后续的全波形反演提供可靠的初始速度模型.理论和实际资料的测试结果证明了本文方法的有效性.     

全波形反演初始模型建立策略研究综述

全波形反演充分利用叠前地震波场中的运动学和动力学信息估计地球内部介质的弹性参数,是一种高精度的速度建模方法和油藏描述的手段.但是对于实际资料,尤其是陆地资料,全波形反演很难收敛到正确的结果上,往往无法实现其建立高精度速度模型的目标,需要较准确的初始模型来降低误差泛函的非线性性.研究与调研认为,建立全波形反演初始模型的方法主要有三种:走时层析、偏移速度分析和Laplace域反演.本文首先分析了全波形反演对初始模型的严重依赖性,其次对全波形反演初始模型建立策略进行简要的介绍与分析,最后给出了建立全波形反演初始模型的认识及建议.     

反射地震数据的逐层波形反演

本文针对层状介质并结合梯度法波形反演，提出逐层波形反演的方法. 首先给出介质扰动响应的概念，并在此基础上分析了梯度法波形反演方法. 波形反演实质上是将实测地震记录和预测地震记录的波形残差信息转化为实际地质模型和预测地质模型的模型残差信息. 波形反演的优点是利用大量振幅相位信息得到高分辨率的反演结果， 其缺点是运行耗时大；当初始模型和实际模型相差较大时，迭代算法容易陷入局部极小点，这是因为目标函数和初始模型同实际模型间的差异是非线性的关系. 逐层波形反演方法是使自上而下每一层的目标函数最小，这样总的目标函数也是最小的. 利用二分法速度扫描确定每一层速度不仅提高了运算速度也避免了迭代算法陷入局部极小点的问题. 结合介质扰动响应和目标函数值变化可以更为准确迅速地确定每一层速度和该层界面位置.     

基于可视性分析与能量补偿的金属矿弹性波全波形反演

采用弹性波全波形反演方法精确重建深部金属矿多参数模型,建模过程采用基于地震照明的反演策略.首先给出基于照明理论的观测系统可视性定义,利用可视性分析构建新的目标函数,对反演目标可视性较高的炮检对接收到的地震记录在波场匹配时占有更高的权重,确保了参与反演计算中的地震数据的有效性;其次将给定观测系统对地下介质的弹性波场照明强度作为优化因子,根据地震波在波阻抗界面处的能量分配特点,自适应补偿波场能量分布和优化速度梯度,以提高弹性波全波形反演过程的稳定性和反演结果的精度.理论模型和金属矿模型反演试验结果表明,基于可视性分析和能量补偿的反演策略可以使弹性波全波形反演更快地收敛到目标函数的全局极小值,获得适用于金属矿高分辨率地震偏移成像的多参数模型.     

一种时间域全波形反演中的地震子波估计方法

全波形反演具有高精度成像能力,然而由理论走向实际应用还存在很多问题.全波形反演通过匹配波形来更新模型,数据的波形与地震子波有直接的关系.本文介绍了一种针对时间域波形反演的子波估计方法,并将其应用到全波形反演中.由于在频率域子波反演可表示为一个线性优化问题,因此本文先给定一个试探子波,在时间域通过有限差分法正演得到地震波场,并在频率域与观测数据比较分析,反演得到预测子波.此外,本文还简要地介绍了如何使用伴随状态法计算全波形反演的梯度.数值实验证明,本文方法反演得到地震子波与真实子波具有很好的吻合度,并在全波形反演中取得了不错的效果.相比不依赖于子波的方法,提前进行子波估计可提高反演效率.此外,本文方法适用于全波形反演中常使用的多尺度反演策略.     

基于多尺度迭代求解三维频率-空间域声波方程的全波形速度反演研究-理论模型

使用广义最小残量方法迭代求解三维频率-空间域声波方程,反演时使用多尺度、多重网格的策略,探讨了如何快速实现高分辨率的三维频率-空间域迭代法声波全波形速度反演.通过对理论模型进行三维频率-空间域迭代法声波全波形反演数值试验,证实该方法的计算速度快、计算效率高,反演所得速度的分辨率高.从而为基于多尺度迭代求解三维频率-空间域声波方程的全波形速度反演成像打下方法基础.     

扩展成像条件下的最小二乘逆时偏移

逆时偏移(RTM)是复杂介质条件下地震成像的重要手段.因受观测系统限制、上覆地层影响以及波场带宽有限等因素的影响,现行的常规RTM所采用的互相关成像条件通常对地下构造进行模糊成像.最小二乘逆时偏移(LSRTM)通过最小化线性Born近似正演数据和采集数据之间的波形差异,采用梯度类反演算法优化反射系数模型,获得的成像结果具有更高的分辨率和更可靠的振幅保真度.然而,基于波形拟合的LSRTM对背景速度模型的依赖性很强.误差太大的速度模型容易产生周波跳跃现象,导致LSRTM难以获得全局最优解.为了克服这一问题,本文基于扩展模型的思想,在线性Born近似下,推导得到RTM扩展成像条件.并基于最小二乘反演理论,提出扩展成像条件下的LSRTM方法.理论模型试算表明,本文方法不仅可以提供分辨率更高、振幅属性更为可靠的成像结果,而且能够在一定程度上消除速度误差对反演成像的影响.