成像域立体层析速度反演方法技术

本发明专利技术提供一种成像域立体层析速度反演方法，包括：步骤1，对输入的初始速度场做叠前深度偏移，得到深度偏移剖面，并抽取角度域共成像点道集；步骤2，在角度域共成像点道集中，拾取与之对应位置处的剩余曲率与深度残差，并转换为走时残差；步骤3，计算立体层析数据空间；步骤4，在当前速度模型中进行射线追踪，得到模型中走时信息和反演核函数；步骤5，利用得到的立体层析核函数以及真实的数据空间与模型中数据空间的残差，建立反演方程组，并计算得到模型的更新量，完成此次迭代。该成像域立体层析速度反演方法可以获得更为准确的立体层析数据空间，提高了反演的精度与稳定性，在今后立体层析反演的实践中具有更好的利用价值。

Stereo Tomography Velocity Inversion Method in Imaging Domain

The invention provides a stereo tomographic velocity inversion method in imaging domain, which includes: step 1, pre-stack depth migration of the input initial velocity field, obtaining depth migration profile, and extracting common imaging point gathers in angle domain; step 2, in common imaging point gathers in angle domain, picking up residual curvature and depth residuals at corresponding positions, and converting them into travel time residuals; step 3, calculating. Stereo tomographic data space; Step 4, ray tracing in current velocity model to get travel time information and inversion kernels; Step 5, using the obtained stereo tomographic kernels and the residuals of real data space and data space in the model, establish inversion equations, and calculate the updates of the model to complete the iteration. The inversion method of stereo tomographic velocity in imaging domain can obtain more accurate stereo tomographic data space, improve the accuracy and stability of inversion, and have better application value in future practice of stereo tomographic inversion.

【技术实现步骤摘要】
成像域立体层析速度反演方法
本专利技术涉及油气勘探地震资料处理领域，特别是涉及到一种成像域立体层析速度反演方法。
技术介绍
地震速度的反演是地震数据处理重要环节。立体层析反演是针对传统反射层析在低信噪比数据中同相轴拾取困难这一问题提出的。该方法重新定义了层析反演的数据空间和模型空间，即利用地震波走时，炮、检点位置以及炮、检点处射线出射方向同时约束速度模型来增强反演的稳定性和精度。在立体层析速度反演中，数据空间拾取的准确性是影响反演效果的最重要因素之一。通常立体层析速度反演所需的数据空间是基于数据域在共炮点与共检波点道集上对局部同相轴做倾斜叠加得到的。这是一种稳定的但是效率较低的数据空间拾取方法，对于二维数据具有可接受的工作量，但对于三维数据则很难接受。并且，利用叠前道集进行倾斜叠加，在交互拾取时对于一次波是有效的，但对于复杂构造中存在的绕射波无法分辨。在立体层析中，获取稳定准确的数据空间是保证反演准确性的前提。因此，针对数据空间拾取不准确与工作量大的问题。为此我们专利技术了一种新的成像域立体层析速度反演方法，解决了以上技术问题。
技术实现思路
本专利技术的目的是提供一种能有效的在成像域来反演速度场中低波数成分的成像域立体层析速度反演方法。本专利技术的目的可通过如下技术措施来实现：成像域立体层析速度反演方法，该成像域立体层析速度反演方法包括：步骤1，对输入的初始速度场做叠前深度偏移，得到深度偏移剖面，并抽取角度域共成像点道集；步骤2，在角度域共成像点道集中，拾取与之对应位置处的剩余曲率与深度残差，并转换为走时残差；步骤3，计算立体层析数据空间；步骤4，在当前速度模型中进行射线追踪，得到模型中走时信息和反演核函数；步骤5，利用得到的立体层析核函数以及真实的数据空间与模型中数据空间的残差，建立反演方程组，并计算得到模型的更新量，完成此次迭代。本专利技术的目的还可通过如下技术措施来实现：在步骤2中，在叠前深度偏移剖面中拾取感兴趣的初始反射点位置(x0,z0)以及构造倾角ξ，并把初始反射点位置对应到角度域共成像点道集中；然后自动拾取剩余深度Δz以及剩余曲率tanφ。在步骤2中，获取炮、检点处的射线参数pr、ps：其中θ为入射角，θs、θr分别为射线在炮点、检波点处的出射角度，v为反射点处速度；由走时残差与深度残差的转化示意图，求取走时残差为：Δt＝2sΔzcosξcosθ其中s表示慢度。在步骤3中，由炮、检点处的射线参数pr、ps，得到炮点和检波点处射线出射时的水平慢度矢量p'sx、p'rx为：p'sx＝ps·sinθ′sp'rx＝pr·sinθ′r其中θ′s、θ′r分别为射线分别在炮点、检波点处的出射角度。在步骤3中，通过下面的修正公式对炮、检点处射线的出射信息进行修正，得到正确的其中ψ＝cosξ·cosθ。在步骤3中，由走时残差与深度残差的转化示意图，得到角道集中走时残差与深度残差的转换关系式为：Δt＝2sΔzcosξcosθ通过修正公式以及转换关系式，得到炮、检点坐标S、R，炮、检点处正确的地表射线出射信息以及正确的双程旅行时t。在步骤4中，利用常速度梯度法进行射线追踪，并用哈密顿Hamiltonian系统对射线路径进行描述：其中x表示射线传播位置信息，n表示射线的方向信息，t表示射线的传播时间，v表示介质中波的传播速度，利用龙格－库塔Runge-Kutta算法求解上式，得到模型中走时信息和反演核函数。在步骤5中，选取最小平方准则建立目标函数：其中Δd＝dtrue-g(m)，CD表示数据空间的协方差矩阵，d表示数据空间，m表示模型空间，g(m)表示数据空间与模型空间关系的非线性算子，由核函数表达，dtrue表示真实数据空间，dmod表示正演模拟得到的数据空间。在步骤5中，根据目标函数，采用局部优化求取极值的方法，得到与其对应的反演方程：其中G为核函数，利用最小二乘正交分解法来求解上述矩阵方程，得到的Δm即为该次迭代模型的更新量。该成像域立体层析速度反演方法还包括，在步骤5之后，步骤6，判断是否进行下一次迭代。在步骤6中，当判断需要进行下一次迭代时，流程返回到步骤1，重复此过程。在步骤6中，当判断不需要进行下一次迭代时，跳出此循环，得到最终的层析速度场。本专利技术中的成像域立体层析速度反演方法，针对传统反射层析在低信噪比数据中同相轴拾取困难这一问题，首先采用在偏移过程中包含地下射线角度信息的DSR(双平方根方程)叠前深度偏移方法提取出共成像点道集，然后从中获取剩余曲率与局部半偏移距射线参数，进而获得走时残差以及相对应炮点和接收点处的射线参数，再通过特定的修正公式，计算得到真实的立体层析数据空间。在初始的成像域道集中进行剩余曲率与构造倾角的拾取，并运用特定的修正公式修正得到更为准确的立体层析数据空间，从而建立立体层析反演方程，在成像域实现了立体层析速度反演，达到较高的速度反演精度。本专利技术利用此方法获得了更为准确的立体层析数据空间，提高了反演的精度与稳定性。这种方法将在今后立体层析反演的实践中具有更好的利用价值。附图说明图1为本专利技术的一具体实施例中立体层析的数据空间与模型空间的示意图；图2为本专利技术的一具体实施例中走时残差与深度残差的转化示意图；图3为本专利技术的一具体实施例中真实速度模型的示意图；图4为本专利技术的一具体实施例中初始速度模型的示意图；图5为本专利技术的一具体实施例中速度场中第400CDP点处速度为真实速度0.85倍时的角度域共成像点道集形态的示意图；图6为本专利技术的一具体实施例中速度场中第400CDP点处速度为真实速度1.0时的角度域共成像点道集形态的示意图；图7为本专利技术的一具体实施例中速度场中第400CDP点处速度为真实速度1.075倍时的角度域共成像点道集形态的示意图；图8为本专利技术的一具体实施例中初始偏移剖面中初始反射点拾取和构造倾角ξ的示意图；图9为本专利技术的一具体实施例中10次迭代后得到的速度场的示意图；图10为本专利技术的一具体实施例中22次迭代后得到的速度场的示意图；图11为本专利技术的一具体实施例中初始速度场偏移剖面的示意图；图12为本专利技术的一具体实施例中立体层析迭代22次后速度场偏移剖面的示意图；图13为本专利技术的一具体实施例中初始角道集的示意图；图14为本专利技术的一具体实施例中立体层析速度场角道集的示意图；图15为本专利技术的一具体实施例中真实速度场角道集的示意图；图16为本专利技术的一具体实施例中距离在3.5km处层析速度、初始速度与真实速度对比图；图17为本专利技术的成像域立体层析速度反演方法的一具体实施例的流程图。具体实施方式为使本专利技术的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合附图所示，作详细说明如下。如图17所示，图17为本专利技术的成像域立体层析速度反演方法的流程图。步骤101，对输入的初始速度场做叠前深度偏移，得到深度偏移剖面，并抽取角度域共成像点道集(ADCIG)。步骤102，在偏移剖面中拾取初始反射点位置以及相应位置处的构造倾角，得到炮、检点处的射线参数。在角度域共成像点道集(ADCIG)中，拾取与之对应位置处的剩余曲率与深度残差，并转换为走时残差。首先在叠前深度偏移剖面中拾取感兴趣的初始反射点位置(x0,z0)以及构造倾角ξ，并把初始反射点位置对应到角度域共成像点道集中。然后自动拾取剩余深本文档来自技高网...

【技术保护点】
1.成像域立体层析速度反演方法，其特征在于，该成像域立体层析速度反演方法包括：步骤1，对输入的初始速度场做叠前深度偏移，得到深度偏移剖面，并抽取角度域共成像点道集；步骤2，在角度域共成像点道集中，拾取与之对应位置处的剩余曲率与深度残差，并转换为走时残差；步骤3，计算立体层析数据空间；步骤4，在当前速度模型中进行射线追踪，得到模型中走时信息和反演核函数；步骤5，利用得到的立体层析核函数以及真实的数据空间与模型中数据空间的残差，建立反演方程组，并计算得到模型的更新量，完成此次迭代。

【技术特征摘要】
1.成像域立体层析速度反演方法，其特征在于，该成像域立体层析速度反演方法包括：步骤1，对输入的初始速度场做叠前深度偏移，得到深度偏移剖面，并抽取角度域共成像点道集；步骤2，在角度域共成像点道集中，拾取与之对应位置处的剩余曲率与深度残差，并转换为走时残差；步骤3，计算立体层析数据空间；步骤4，在当前速度模型中进行射线追踪，得到模型中走时信息和反演核函数；步骤5，利用得到的立体层析核函数以及真实的数据空间与模型中数据空间的残差，建立反演方程组，并计算得到模型的更新量，完成此次迭代。2.根据权利要求1所述的成像域立体层析速度反演方法，其特征在于，在步骤2中，在叠前深度偏移剖面中拾取感兴趣的初始反射点位置(x0,z0)以及构造倾角ξ，并把初始反射点位置对应到角度域共成像点道集中；然后自动拾取剩余深度Δz以及剩余曲率tanφ。3.根据权利要求2所述的成像域立体层析速度反演方法，其特征在于，在步骤2中，获取炮、检点处的射线参数pr、ps：其中θ为入射角，θs、θr分别为射线在炮点、检波点处的出射角度，v为反射点处速度；由走时残差与深度残差的转化示意图，求取走时残差为：Δt＝2sΔzcosξcosθ其中s表示慢度。4.根据权利要求3所述的成像域立体层析速度反演方法，其特征在于，在步骤3中，由炮、检点处的射线参数pr、ps，得到炮点和检波点处射线出射时的水平慢度矢量p′sx、p′rx为：p′sx＝ps·sinθ′sp′rx＝pr·sinθ′r其中θ′s、θ′r分别为射线分别在炮点、检波点处的出射角度。5.根据权利要求4所述的成像域立体层析速度反演方法，其特征在于，在步骤3中，通过下面的修正公式对炮、检点处射线的出射信息进行修正，得到正确的其中ψ＝cosξ·cosθ。6.根据权利要求5所述的成像域立体层析速度反演方法，其特征在于，在步骤3中，...

【专利技术属性】
技术研发人员：田坤，于海铖，王常波，刘立彬，张学涛，步长城，王蓬，徐瑞，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司胜利油田分公司物探研究院，
类型：发明
国别省市：山东,37