[发明专利]应用于实际区全波形反演的地震资料预处理方法及系统

说明书

本发明涉及地球物理勘探技术领域，特别涉及一种应用于实际区全波形反演的地震资料预处理方法及系统。该方法包括：获取地震炮集数据和初始模型数据，并进行坐标旋转；将地震炮集数据进行波场分离，拆分直达波数据和有效波数据，并基于直达波数据提取震源子波；根据稳定性条件和频散条件，对初始模型数据进行网格规则化；采用快速线性插值算法对炮集有效波数据进行重采样和归一化；建立炮集观测系统索引数据；开展实际区叠前全波形反演。本发明通过针对实际区全波形反演的地震资料预处理，很好的解决了现有地震资料预处理方法对实际区全波形反演缺乏针对性的问题，提高了全波形反演在复杂实际区地震资料中的适用性。

技术领域

本发明涉及地球物理勘探技术领域，特别涉及一种应用于实际区全波形反演的地震资料预处理方法及系统。

背景技术

全波形反演是基于全波波动方程正演模拟技术，将正演波场与观测波场进行匹配，将数据残差最小化作为目标函数，从地震数据中反演地球物理参数的方法。利用地震反射波形的全波形反演方法是目前理论上精度最高的地震反演方法，能够充分利用地震波的走时、振幅、频率和相位等全波形信息，是建立高精度速度模型的有效手段。

现如今，全波形反演从二维到三维，从单参数到多参数，从各向同性介质到衰减介质、各向异性等复杂介质，逐步从理论模型走向了实际应用。然而，在实际工区中应用全波形反演还面临诸多的实际问题，集中表现在：观测系统不规则、测网方向不统一、资料数据量庞大、测网间隔不满足差分计算条件等，严重制约了全波形反演方法的实际应用推广。

中国专利申请CN106842295A公开了一种测井信息约束的波形反演方法，包括根据初始速度模型，对测井数据空间插值，构建准确的先验速度模型；利用交错网格有限差分方法合成地震记录，计算剩余波场；建立测井数据约束正则化反演方程，求取速度更新量，更新速度模型，通过地震正演得到的正传波场和反向传播波场并在测井数据的约束下计算梯度方向；对梯度方向进行修正，在时间域中运用共轭梯度的方法实现全波形的计算；通过迭代公式更新速度参数判断是否符合迭代收敛条件，输出结果。

中国专利申请CN109471167A公开了一种针对多震源缺失数据的波场重构反演方法，该方法包括采用声波近似正演算子，对多震源数据进行正演模拟，得到多震源地震数据；对多震源地震数据进行随机去除，以模拟实际生产中地震观测数据存在的缺道；采用波场重构反演方法，对随机去除后的多震源数据进行反演，得到多震源缺失数据的重构波场。

中国专利申请CN105389196A公开了一种利用虚拟化技术进行地震资料处理的方法，该方法包括：分析地震资料处理软件运行资源需求特点，根据资源需求特点将地震资料处理应用分成不同的资源需求类型，将节点分成不同的服务器类型，建立服务器模型；根据地震资料处理资源需求特点优化主机资源池，建立适于地震资料处理应用的虚拟存储和虚拟网络服务；根据地震资料处理服务器模型分别建立虚拟机模板，并优化存储和网络驱动；根据应用需求调度主机资源，部署相应的虚拟机服务器，组成虚拟集群实现地震资料处理软件的虚拟化应用。该方法通过虚拟化技术提高了地震资料处理中心设备整体资源利用率。

发明人发现，现有的地震资料预处理方法大都针对常规地震资料处理，无法针对性解决实际资料全波形反演的数据预处理需求，因此亟待开发一种应用于实际区全波形反演的地震资料预处理方法及系统。

发明内容

针对全波形反演在实际应用上的不足，本发明提供一种应用于实际区全波形反演的地震资料预处理方法及系统，本发明方法通过对地震资料的预处理，提高了全波形反演在复杂实际区地震资料中的适用性。

为实现上述目的，本发明采用以下技术方案：

本发明提供一种应用于实际区全波形反演的地震资料预处理方法，其包括以下步骤：

获取地震炮集数据和初始模型数据，并进行坐标旋转；

将地震炮集数据进行波场分离，拆分直达波数据和有效波数据，并基于直达波数据提取震源子波；