TTI介质层析反演成像方法及计算机可读存储介质技术

本发明专利技术公开了一种TTI介质层析反演成像方法及计算机可读存储介质，包括：输入数据；基于TTI介质的各向异性参数(vp0,ε,δ,θ0)其中的各向异性参数ε和δ进行反演，获取精确的参数ε和δ的结果；基于精确的参数ε和δ的结果，对各向异性参数ε、δ和vp0进行反演，获取精确的参数ε、δ和vp0的结果；基于精确的参数ε、δ和vp0的结果，对各向异性参数ε、δ、vp0和θ0进行反演，获取所述各向异性参数(vp0,ε,δ,θ0)的反演结果。其优点在于：通过TTI各向异性介质中不同参数对旅行时的不同敏感程度进行反演，大大提高了TTI介质层析建模的精度，并能够提高偏移成像的精度。

TTI dielectric tomography inversion imaging method and computer readable storage medium

The invention discloses a TTI medium tomography inversion imaging method and a computer readable storage medium, including: input data; inversion of anisotropic parameters E and E based on anisotropic parameters (vp0, e, delta, theta 0) of TTI medium to obtain accurate results of parameters E and delta; and inversion of the anisotropic parameters E and Delta Based on accurate results of parameters E and delta. The anisotropic parameters e, Delta and VP0 are inverted to obtain the exact parameters e, Delta and vp0. Based on the exact parameters e, Delta and vp0, the anisotropic parameters e, delta, VP0 and theta 0 are inverted to obtain the inversion results of the anisotropic parameters (vp0, e, delta, theta 0). Its advantage is that the accuracy of TTI medium tomography modeling and migration imaging can be greatly improved by inversion of different parameters in TTI anisotropic medium to different sensitivity of travel time.

【技术实现步骤摘要】
TTI介质层析反演成像方法及计算机可读存储介质
本专利技术涉及地震勘探领域，更具体地，涉及一种TTI介质层析反演成像方法及计算机可读存储介质。
技术介绍
传统的地震偏移成像技术大都是基于地球介质为各向同性的假设，但是大量的观察和实验表明地球内部介质的各向异性是普遍存在的，如果运用基于各向同性的假设的偏移成像方法来处理各向异性探区的地震数据，将会严重影响速度分析的精度，继而降低地震偏移成像的准确性。各向异性可分为VTI垂直横向各向异性、TTI倾斜横向各向异性和HTI方位各向异性。其中，VTI介质目前研究的比较多，而对TTI和HTI则研究的较少。而对应到地震勘探中的层析反演速度建模领域中，由于TTI介质各参数量级不统一的问题，算法一直都在致力于研究不同反演策略对TTI多参数反演的影响，因而造成一定程度上TTI介质层析速度建模研究的滞后。因此，有必要基于不同TTI各向异性参数对旅行时的敏感性，开发一种TTI介质层析反演成像方法。公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
本专利技术提出了一种TTI介质层析反演成像方法及计算机可读存储介质，其通过研究不同TTI各向异性参数对旅行时的敏感性，采取参数化手段和策略进行层析反演对速度模型进行迭代更新，提供一种精确的TTI各向异性的速度模型，为偏移成像打下良好基础。根据本专利技术的一方面，提出了一种TTI介质层析反演成像方法，包括：输入数据；基于TTI介质的各向异性参数(vp0,ε,δ,θ0)其中的各向异性参数ε和δ进行反演，获取精确的参数ε和δ的结果；基于所述精确的参数ε和δ的结果，对所述各向异性参数ε、δ和vp0进行反演，获取精确的参数ε、δ和vp0的结果；基于所述精确的参数ε、δ和vp0的结果，对所述各向异性参数ε、δ、vp0和θ0进行反演，获取所述各向异性参数(vp0,ε,δ,θ0)的反演结果；其中，vp0为沿对称轴方向的qp波相速度；ε和δ为Thomsen参数；θ0为波前法线方向与对称轴的夹角。优选地，对所述各向异性参数ε和δ反演时，所述各向异性参数vp0和θ0设为0，进行多次迭代获取所述各向异性参数ε和δ的反演结果。优选地，对所述各向异性参数ε、δ和vp0反演时，所述各向异性参数θ0设为0，进行多次迭代获取所述各向异性参数ε、δ和vp0的反演结果。优选地，对所述各向异性参数ε、δ、vp0和θ0反演，进行多次迭代获取所述各向异性参数ε、δ、vp0和θ0的反演结果。根据本专利技术的另一方面，提出了一种计算机可读存储介质，其上存储计算机程序，其中，所述程序被处理器执行时实现以下步骤：基于TTI介质的各向异性参数(vp0,ε,δ,θ0)其中的各向异性参数ε和δ进行反演，获取精确的参数ε和δ的结果；基于所述精确的参数ε和δ的结果，对所述各向异性参数ε、δ和vp0进行反演，获取精确的参数ε、δ和vp0的结果；基于所述精确的参数ε、δ和vp0的结果，对所述各向异性参数ε、δ、vp0和θ0进行反演，获取所述各向异性参数(vp0,ε,δ,θ0)的反演结果；其中，vp0为沿对称轴方向的qp波相速度；ε和δ为Thomsen参数；θ0为波前法线方向与对称轴的夹角。优选地，对所述各向异性参数ε和δ反演时，所述各向异性参数vp0和θ0设为0，进行多次迭代获取所述各向异性参数ε和δ的反演结果。优选地，对所述各向异性参数ε、δ和vp0反演时，所述各向异性参数θ0设为0，进行多次迭代获取所述各向异性参数ε、δ和vp0的反演结果。优选地，对所述各向异性参数ε、δ、vp0和θ0反演，进行多次迭代获取所述各向异性参数ε、δ、vp0和θ0的反演结果。根据本专利技术的一种TTI介质层析反演成像方法及计算机可读存储介质，其优点在于：通过TTI各向异性介质中不同参数对旅行时的不同敏感程度进行反演，大大提高了TTI介质层析建模的精度，并能够提高偏移成像的精度。本专利技术的方法和计算机可读存储介质具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本专利技术的特定原理。附图说明通过结合附图对本专利技术示例性实施例进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施例中，相同的附图标记通常代表相同部件。图1示出了TTI介质各向异性参数对旅行时的敏感度图。图2示出了根据本专利技术的一种TTI介质层析反演成像方法的步骤的流程图。图3示出了实际模型中TTI介质各向异性参数中参数vp0的成像示意图。图4a和图4b分别示出了根据本专利技术的一个示例性实施例的ε和δ双参数反演结果示意图。图5a、图5b和图5c分别示出了根据本专利技术的一个示例性实施例的vp0、ε和δ三参数反演结果示意图。图6a、图6b、图6c和图6d分别示出了根据本专利技术的一个示例性实施例的vp0、ε、δ和θ0四参数反演结果示意图。具体实施方式下面将参照附图更详细地描述本专利技术。虽然附图中显示了本专利技术的优选实施例，然而应该理解，可以以各种形式实现本专利技术而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本专利技术更加透彻和完整，并且能够将本专利技术的范围完整地传达给本领域的技术人员。由于TTI各向异性参数(vp0,ε,δ,θ0)对旅行时的影响不同，反演更新量的尺度也会不同，因此对这些参数的敏感度分析是必要的。旅行时与TTI介质各向异性参数的关系为：其中，vp0为沿对称轴方向的qp波相速度；ε和δ为Thomsen参数；θ0为波前法线方向与对称轴的夹角。将式(1)中的右边求导项移至左边，变形得到：图1示出了TTI介质各向异性参数对旅行时的敏感度图。通过式(2)可以绘制TTI介质各向异性参数对旅行时的敏感度的曲线，如图1所示，可以看出，TTI介质各向异性参数对旅行时的敏感度是不同的，其中，ε最大，δ次之，vp0和θ0最不敏感，因此本专利技术决定先对最敏感的ε和δ两个参数进行反演策略设置。根据本专利技术一种TTI介质层析反演成像方法，可以包括：输入数据；基于TTI介质的各向异性参数(vp0,ε,δ,θ0)其中的各向异性参数ε和δ进行反演，获取精确的参数ε和δ的结果；基于精确的参数ε和δ的结果，对各向异性参数ε、δ和vp0进行反演，获取精确的参数ε、δ和vp0的结果；基于精确的参数ε、δ和vp0的结果，对各向异性参数ε、δ、vp0和θ0进行反演，获取所述各向异性参数(vp0,ε,δ,θ0)的反演结果。其中，vp0为沿对称轴方向的qp波相速度；ε和δ为Thomsen参数；θ0为波前法线方向与对称轴的夹角。根据TTI介质各向异性参数对旅行时的敏感程度的不同，提出针对性的反演策略，能够大大提高TTI介质层析建模的精度。地球物理界中，层析反演是一种已知介质各向异性参数在可能不精确的前提下，射线穿越介质的计算旅行时与观测旅行时(即实际旅行时)存在一定残差，这部分残差是由介质各向异性参数不精确引起的，因此可以根据介质各向异性参数沿射线路径的积分即是旅行时残差这一理论对介质参数进行反演。因此根据以上已推导出的各向异性参数与旅行时残差的本文档来自技高网...

【技术保护点】
1.一种TTI介质层析反演成像方法，包括：输入数据；基于TTI介质的各向异性参数(vp0,ε,δ,θ0)其中的各向异性参数ε和δ进行反演，获取精确的参数ε和δ的结果；基于所述精确的参数ε和δ的结果，对所述各向异性参数ε、δ和vp0进行反演，获取精确的参数ε、δ和vp0的结果；基于所述精确的参数ε、δ和vp0的结果，对所述各向异性参数ε、δ、vp0和θ0进行反演，获取所述各向异性参数(vp0,ε,δ,θ0)的反演结果；其中，vp0为沿对称轴方向的qp波相速度；ε和δ为Thomsen参数；θ0为波前法线方向与对称轴的夹角。

【技术特征摘要】
1.一种TTI介质层析反演成像方法，包括：输入数据；基于TTI介质的各向异性参数(vp0,ε,δ,θ0)其中的各向异性参数ε和δ进行反演，获取精确的参数ε和δ的结果；基于所述精确的参数ε和δ的结果，对所述各向异性参数ε、δ和vp0进行反演，获取精确的参数ε、δ和vp0的结果；基于所述精确的参数ε、δ和vp0的结果，对所述各向异性参数ε、δ、vp0和θ0进行反演，获取所述各向异性参数(vp0,ε,δ,θ0)的反演结果；其中，vp0为沿对称轴方向的qp波相速度；ε和δ为Thomsen参数；θ0为波前法线方向与对称轴的夹角。2.根据权利要求1所述的TTI介质层析反演成像方法，其中，对所述各向异性参数ε和δ反演时，所述各向异性参数vp0和θ0设为0，进行多次迭代获取所述各向异性参数ε和δ的反演结果。3.根据权利要求1所述的TTI介质层析反演成像方法，其中，对所述各向异性参数ε、δ和vp0反演时，所述各向异性参数θ0设为0，进行多次迭代获取所述各向异性参数ε、δ和vp0的反演结果。4.根据权利要求1所述的TTI介质层析反演成像方法，其中，对所述各向异性参数ε、δ、vp0和θ0反演，进行多次迭代获取所述各向异性参数ε、δ、vp0和θ0的反演结果。5.一...

【专利技术属性】
技术研发人员：张林，张兵，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11