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

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

3.
多测线变偏移距VSP地震各向异性反演   总被引:1,自引:0,他引:1       下载免费PDF全文
田鑫  洪启宇  郑需要 《地震学报》2015,37(2):266-277
本文从一般弱各向异性介质参数反演中得到了使用两条相互正交的变偏VSP测线计算弱各向异性(WA)参数的反演公式. 如果仅仅使用qP波, 则可以确定9个独立的WA参数, 这9个WA参数可以完全地描述井中接收点在两个剖面内介质的各向异性性质. 通过对走时曲线进行最小二乘条件下的三次样条光滑, 可以获得慢度矢量的垂直分量和水平分量. 如果介质是横向非均匀介质, 则水平分量仅仅被用作反演时的约束条件. 为了获得偏振矢量, 本文引入质心计算方法, 该方法计算简单、 稳定, 而且不存在解的奇异问题. 在寻找与一般各向异性介质最接近的高对称性正交各向异性介质和TTI各向异性介质时, 使用qP波各向异性坐标变换方法和最小二乘求解方法, 得到了与一般各向异性介质最接近的正交各向异性和TTI各向异性参数及其对称轴方向参数的计算公式. 使用这些方法, 对瓜哇海地区布设的两条相互正交的变偏VSP测线数据进行各向异性反演, 获得了井中10个接收点处介质的WA参数. 数值计算和实际资料反演表明, 本文所使用的反演方法能够准确地得到VSP井中接收点处介质的WA参数, 这是地震勘探中研究地壳介质各向异性性质最直接和最可靠的方法.   相似文献   

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

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

6.
得益于结合最优化算法和完整的波动方程,全波形反演已成为地震勘探最前沿的研究方向,正逐渐发展成为获取地球内部信息的重要工具.但是,在初始速度模型不够精准的条件下,利用纯反射波进行波形反演来更新光滑的背景速度模型是非常困难的.本文展示了一种利用纯反射波通过交替进行最小二乘逆时偏移和波形反演的方法来成功地同时反演反射系数和背景速度的方法.其波形反演中用于更新背景速度的层析梯度是通过基于Born近似的正演模拟来分离反射波波场得到的,这样可以很好地更新光滑的背景速度模型.之后,常规全波形反演把上一步得到的光滑背景速度模型作为初始模型,利用纯反射波获得最终的高保真、高分辨率的速度模型.进一步分析揭示了高分辨率的速度模型是由于在全波形反演中应用了反褶积成像条件.相对于其他类型的地震波,反射波带有地下更深部的信息,因此,反射波波形反演能够进一步提高我们对地球内部的成像精度.  相似文献   

7.
本文研究了巴黎盆地一个井场的多偏移距VSPs剪切波分裂,并证实沉积构造至少必须有正交的各向异性对称,认为引起这种正交对称的最可能的原因是薄互层(或岩性)各向异性和裂隙各向异性的聚形(combination)。过渡带(靠近剪切波点状奇异区(pointsingularities)发生剪切波偏振和时滞的快速变化是造成异常观测波形的原因,观测到剪切波偏振在近点状奇异性处的传播方向上有很大的变化这个结果可能是第一次明确无误地识别出沉积盆地存在着点状奇异性的效应.  相似文献   

8.
本文使用炮并行和区域分解(物理上分割模型,使用基于MPI的分布式存储架构的计算集群,节约单个CPU内核的内存使用量,快速进行正演数值模拟)两种并行算法.该方法的每一步迭代都能确保近似海森矩阵的正定,因此,算法稳健.将时间正向传播的炮波场和反向逆时间传播的残差波场(伴随波场)进行零延迟互相关计算,得到误差泛函的梯度,然后对梯度乘以一个预条件算子,从而加快反演的收敛速度.通过抛物线搜索方法而估计步长,使用L-BFGS算法(限定内存的BFGS算法)求解模型的更新量,进行二维时间空间域弹性波全波形反演.将该反演方法应用到Marmousi2弹性波理论模型,分别反演Marmousi2理论模型的纵波速度、横波速度以及密度等三个参数.我们分别使用截止频率为2 Hz、5 Hz、10 Hz和20 Hz四个阶段的低通巴特沃斯滤波器,采用多尺度的策略,从理论模型数据的低频分量开始反演,将低频分量的反演结果作为高频分量反演时的初始模型,然后依次反演数据的高频分量.理论模型数值试验反演所得到的结果证实:二维时间空间域弹性波全波形反演计算灵活,适用于各种观测系统,能够方便地对地震数据进行加时窗;二维时间空间域弹性波全波形反演所得纵波速度模型的分辨率最高,横波速度模型的分辨率次之,密度模型的分辨率稍微差些.  相似文献   

9.
时间二阶积分波场的全波形反演   总被引:4,自引:4,他引:0       下载免费PDF全文
陈生昌  陈国新 《地球物理学报》2016,59(10):3765-3776
通过对波场的时间二阶积分运算以增强地震数据中的低频成分,提出了一种可有效减小对初始速度模型依赖性的地震数据全波形反演方法—时间二阶积分波场的全波形反演方法.根据散射理论中的散射波场传播方程,推导出时间二阶积分散射波场的传播方程,再利用一阶Born近似对时间二阶积分散射波场传播方程进行线性化.在时间二阶积分散射波场传播方程的基础上,利用散射波场反演地下散射源分布,再利用波场模拟的方法构建地下入射波场,然后根据时间二阶积分散射波场线性传播方程中散射波场与入射波场、速度扰动间的线性关系,应用类似偏移成像的公式得到速度扰动的估计,以此建立时间二阶积分波场的全波形迭代反演方法.最后把时间二阶积分波场的全波形反演结果作为常规全波形反演的初始模型可有效地减小地震波场全波形反演对初始模型的依赖性.应用于Marmousi模型的全频带合成数据和缺失4Hz以下频谱成分的缺低频合成数据验证所提出的全波形反演方法的正确性和有效性,数值试验显示缺失4Hz以下频谱成分数据的反演结果与全频带数据的反演结果没有明显差异.  相似文献   

10.
芦俊  石瑛  杨春颖 《地球物理学报》2018,61(8):3310-3323
针对裂缝各向异性介质,本文提出一种非正交假设下的矢量波场分离方法.本文首先对多分量地震勘探中常见的波型泄漏现象进行了数学描述,提出在纵、横波波场分离的同时应该考虑恢复纵、横波的矢量振幅.为了对裂缝方位角与各向异性系数进行定量预测,本文将矢量波场分离拆分成三个步骤来实施:第一步,用Z、R两分量的仿射坐标系变换分离ZR平面内的P波投影与SV波;第二步,用ZR平面内的P波投影与T分量的仿射坐标系变换分离P波与SH波;第三步,用纯净的SV波与SH波的成像剖面分离快慢横波,并预测裂缝发育参数.模型数据与实际数据的试验结果表明,本文提出的纵、横波波场分离方法能够获得完整的矢量振幅信息,并提供裂缝预测的精度.  相似文献   

11.
多分量地震波波场分解研究   总被引:17,自引:2,他引:15       下载免费PDF全文
提出了一种有效的多分量地震波波场分解方法, 该方法假设在相邻道集波场的振幅、相位变化不大的情况下, 根据纵波和转换波在传播中相对于不同的极化矢量和慢度得到的垂向和水平向分量, 在最小方差约束下从多分量输入数据中获取纵波和转换波的偏振方向和振幅并重构纵波和转换波, 从而分离纵波和转换波. 这种方法精度较高, 不仅能有效分解一次反射波的多分量波场, 多次反射的纵波和转换波也能被恰当地分解. 该方法不需要输入地层速度等参数信息, 减少了人为误差的影响. 通过对合成资料的处理, 验证了方法的可行性和正确性.  相似文献   

12.
张盼  邢贞贞  胡勇 《地球物理学报》2019,62(10):3974-3987
在常规地震采集中,被动源地震波场往往被视为噪声而去除,这就造成了部分有用信息的丢失.在目标区进行主动源和被动源弹性波地震数据的多分量混合采集,并对两种数据进行联合应用,使其在照明和频带上优势互补,能显著提高成像和反演的质量.本文针对两种不同类型的主被动源混采地震数据,分别提出了相应的联合全波形反演方法.首先,针对主动源与瞬态被动源弹性波混采地震数据,为充分利用被动源对深部照明的优势,同时有效压制被动震源点附近的成像异常值,提出了基于动态随机组合的弹性波被动源照明补偿反演策略.然后,针对低频缺失主动源与背景噪声型被动源弹性波混采地震数据,为充分利用被动源波场携带的低频信息,并避免对被动源的定位和子波估计,提出了基于地震干涉与不依赖子波算法的弹性波主被动源串联反演策略.最后,分别将两种方法在Marmousi模型上进行反演测试.结果说明,综合利用主动源和被动源弹性波混采地震数据,不仅能增强深部弹性参数反演效果,还能更好地构建弹性参数模型的宏观结构,并有助于缓解常规弹性波全波形反演的跳周问题.  相似文献   

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

14.
转换波方位各向异性裂缝检测技术研究及应用   总被引:1,自引:0,他引:1  
HTI裂缝各向异性介质中,转换波随方位角的变化比较复杂,目前还没有解析公式可以表达其变化特征,只能通过物理实验或数值模拟来分析其应用的可能性.数值模拟结果表明,转换波在裂缝各向异性介质中传播时,其R分量和T分量的振幅属性都具有明显的方位各向异性特征,R分量振幅方位各向异性拟合椭圆的长轴方向指示裂缝方位,这与纵波方位各向异性特征相似;根据P波AVAZ方位各向异性分析原理,对转换波R分量振幅方位各向异性曲线进行方位椭圆拟合,寻找椭圆的长轴方向,即裂缝主方位,再由P波AVAZ技术中振幅响应与炮检方向和裂缝走向之间的夹角关系式得到裂缝的发育密度,从而构建转换波方位各向异性AVAZ裂缝检测技术.该技术已用于川西新场气田某区块的裂缝储层预测,取得了较好的应用效果.  相似文献   

15.
高精度的成像需要准确的速度场信息,波形反演被认为是目前具有最高分辨率的速度反演方法之一.计算效率是目前全波形反演需要考虑的一个主要问题.为了很好地解决计算效率的问题,本文引入了一种高效的无串扰同时源反演方法,并详细介绍了其原理与计算流程.同时,基于此多震源同时反演方法,本文拓展出实波场反演、虚波场反演及复波场反演的三种反演策略,进而分析复波场的实部和虚部对全波形反演的影响.实验表明,相比于复波场全波形反演,无论是实波场反演还是虚波场反演,反演分辨率有所降低;相比于实波场反演,虚波场反演对初始模型的依赖性较小,目标函数的非线性较弱;最后,通过使用组合反演策略,即初始阶段采用虚波场反演,中后期阶段采用复波场反演,不仅可以降低反演的非线性,而且能够保证高精度建模.  相似文献   

16.
低频成分缺失和地下速度强烈变化会导致严重的周期跳现象,是地震数据全波形反演的难题.通过对地震数据加时间阻尼和时间积分降主频处理,提出了一种可有效去除周期跳现象的多主频波场时间阻尼全波形反演方法.由浅到深的速度不准确会造成波形走时失配和走时失配的累积.浅部速度的准确反演可有效地减小深部波形走时失配与周期跳现象.对地震数据施加时间阻尼得到时间阻尼数据,利用不同阻尼值的时间阻尼地震数据实现由浅到深的全波形反演.低主频波场的周期跳现象相对高主频波场的要弱.对地震波场进行不同阶的时间积分以得到不同主频的波场,把低主频波场的全波形反演结果作为高主频波场全波形反演的初始模型.应用缺失4 Hz以下频谱成分的二维盐丘模型合成数据验证所提出的全波形反演方法的正确性和有效性,数值试验结果显示多主频波场的时间阻尼全波形反演方法对缺失低频成分地震数据和地下速度强烈变化具有很好的适应性.  相似文献   

17.
多分量地震数据中低频缺失是弹性波全波形反演中的一大难题,低频的缺失导致全波形反演无法有效恢复介质的长波长成分进而使反演陷入局部极值。为此,本文提出了反演介质纵横波速度长波长分量的弹性波包络反演方法。该方法利用包络算子具有的解调多分量数据中隐含的低频信息的能力,构造多分量地震数据的包络目标函数进行反演,用以恢复地下介质纵横波速度的长波长成分。一系列数值试验表明,即使在多分量地震数据中缺失低频信息、并且初始模型缺少先验信息的情况下,这种弹性波包络反演方法能够有效降低波形反演的非线性,可以为后续的常规弹性波全波形反演或者深度偏移提供足够精确的初始模型,且该方法对横波速度长波长分量的重建尤为有效。Mamousi-2模型的高精度纵横波速度的反演结果表明,利用该方法反演的纵横波速度作为常规弹性波全波形反演的初始模型,可以显著提高反演结果的精度。此外,本文对弹性波包络反演方法的适用性也进行了初步的研究与讨论。  相似文献   

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

19.
为克服各向异性弹性波动方程正演模拟的局限,本文研究了各向异性介质拟声波方程的交错网格有限差分数值解法.首先,从VTI介质胡克定律和qP-qSV波频散关系两种思路出发,通过声假设近似,给出了两种不同形式的VTI介质一阶拟声波方程,并通过引入波场的伪速度分量,推导了一种新的VTI介质一阶应力-速度方程,并通过旋转坐标系将其推广到TTI介质中;其次,构造了一阶拟声波方程的交错网格高阶有限差分格式,并推导了相应的PML边界条件;最后,对本文方法中固有的qSV人为干扰波的产生机制和压制方法进行了简单讨论.数值结果表明:3种一阶拟声波方程在运动学和动力学上是等价的,相对于各向异性弹性波正演模拟,其节省了内存,提高了计算效率;各向异性因素会影响反射波旅行时和振幅等波场特征,在后续的处理、反演和解释中不可忽略;VTI介质HESS模型的逆时偏移结果也验证了本文方法的合理性.   相似文献   

20.
高斯波包反射走时速度反演方法   总被引:1,自引:1,他引:0       下载免费PDF全文
李辉  殷俊锋  王华忠 《地球物理学报》2017,60(10):3916-3933
扰动高斯波包理论指出,在Gabor域描述模型的扰动成分,且入射波场为短时宽带信号时,扰动波场可在时间域通过高斯波包算子描述.在此基础上通过拟合反射波的走时,提出一种速度反演方法.反射波走时残差利用地震道局部波形的互相关函数表示,以走时残差的二范数作为目标函数,优化目标函数实现对速度场的反演.基于一阶Born近似,利用扰动高斯波包理论推导出目标函数对速度场的梯度是本文理论部分的核心内容.梯度包括两部分:正传的背景波场与反传的扰动高斯波包之间的互相关,反传的背景波场和正传的扰动高斯波包之间的互相关.梯度表达式中背景波场和扰动波场均利用高斯波包算子模拟.计算梯度的具体算法中,如何模拟扰动波场,以及如何计算反射波的走时残差是两个要点,文中对此做了详细的讨论.数值实验进一步阐述了反演的实现策略,实验结果表明高斯波包反射走时速度反演方法和实现策略有效可行,并得到了理想的反演结果.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号