CN115327626A - 一种用于三维ati介质的弹性波场矢量分解方法及系统 - Google Patents

一种用于三维ati介质的弹性波场矢量分解方法及系统 Download PDF

Info

Publication number
CN115327626A
CN115327626A CN202210960164.3A CN202210960164A CN115327626A CN 115327626 A CN115327626 A CN 115327626A CN 202210960164 A CN202210960164 A CN 202210960164A CN 115327626 A CN115327626 A CN 115327626A
Authority
CN
China
Prior art keywords
dimensional
ati
wave field
medium
elastic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210960164.3A
Other languages
English (en)
Other versions
CN115327626B (zh
Inventor
左佳卉
赵杨
张乐乐
陈科
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN202210960164.3A priority Critical patent/CN115327626B/zh
Publication of CN115327626A publication Critical patent/CN115327626A/zh
Application granted granted Critical
Publication of CN115327626B publication Critical patent/CN115327626B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)

Abstract

本发明涉及一种用于三维ATI介质的弹性波场矢量分解方法及系统，包括：对观测坐标系进行旋转，使得观测坐标系的垂向坐标轴平行于各向异性介质的倾斜对称轴方向，得到新观测坐标系；基于设置的三维ATI模型的弹性参数和Thomsen参数，计算三维VTI介质条件下的刚度矩阵，建立三维VTI介质弹性波动方程；在新观测坐标系下，采用三维VTI介质弹性波动方程描述三维ATI介质，确定空间域的三维ATI分解算子；在每一次时间迭代过程中，对弹性波场进行正反傅里叶变换，并基于设置的三维ATI模型的弹性参数和Thomsen参数，得到空间域的辅助波场；根据空间域的三维ATI分解算子和辅助波场，得到P波场和S波场，本发明可广泛用于地球物理数据处理和弹性参数反演领域中。

Description

一种用于三维ATI介质的弹性波场矢量分解方法及系统
技术领域
本发明涉及地球物理数据处理和弹性参数反演领域，特别是关于一种用于三维ATI(Arbitrary Transverse Isotropy，任意横向各向同性)介质的弹性波场矢量分解方法及系统。
背景技术
实际地下介质普遍存在各向异性特征，通常由定向排列的垂直裂缝或具有周期性的薄互层引起。当地震波穿过这些介质时，由于波长大于裂缝宽度或者薄层厚度，常出现地震各向异性特征。地层各向异性对称轴沿垂直方向时，可采用VTI(VerticalTransverseIsotropy，垂直横向各向同性)模型近似。若经过强烈的构造运动，对称轴方向发生改变，不再沿着垂直方向，当倾斜对称轴在水平面的投影具有一定方位时，各向异性介质可以采用任意倾向的TI介质进行近似，此时对称轴的空间取向可以为任意方向，简写为ATI。面对这种越来越复杂的各向异性油气储层，现有的地震勘探方法已经无法满足当前高精度油气勘探的需求。
目前，工业界常用的地震偏移方法主要包括射线类和波动方程类两类，基于弹性波动方程的逆时偏移方法由于不受倾角和偏移孔径的限制，可以较好地保存反映各向异性地层结构的地震波振幅和相位等信息，所以逐渐应用到实际勘探生产中。该方法可以在地下重建P波场和S波场，更符合真实的地震波传播过程，应用成像条件可以得到不同波模式的成像结果。但是在弹性各向异性地层中，耦合传播的P波场和S波场将产生串扰噪声，降低成像剖面的质量。特别是针对复杂的三维ATI介质来说，更需要研究一种高效精确的弹性波场解耦方法，以提高弹性逆时偏移的成像精度。
传统弹性波场分解基于亥姆霍兹定理，通过旋度和散度算子实现P波场和S波场的分离，但是该算子认为波的偏振方向和传播方向一样，所以只能应用在各向同性介质中。面对各向异性介质，通过求解克里斯托弗方程可以计算实际波的偏振方向，进而实现P波场和S波场的分离。许多学者也基于这个思路，研究了弹性波场的矢量分解方法，通过在二维非均匀VTI介质中，利用非稳态滤波实现P波场和S波场的分离。另外一种在均匀VTI介质中的弹性波场分解方法，使分解前后的波场具有相同的振幅、相位和物理单位。随后低秩近似的引入大大简化了数值计算的过程，同时，有学者将该方法推广到二维非均匀VTI介质中。
然而，当面对更复杂的各向异性介质时，例如三维ATI介质，这些适用于VTI模型的弹性波场分解方法由于无法考虑对称轴的倾角和方位信息，不能对弹性波场进行有效分离。另外，随着实际资料采集数量的增加，这些分解方法的计算成本依然很高，始终无法满足实际三维生产的要求。
发明内容
针对上述问题，本发明的目的是提供一种用于三维ATI介质的弹性波场矢量分解方法及系统，能够考虑对称轴的倾角和方位信息，且满足实际三维生产要求。
为实现上述目的，本发明采取以下技术方案：第一方面，提供一种用于三维ATI介质的弹性波场矢量分解方法，包括：
对观测坐标系进行旋转，使得观测坐标系的垂向坐标轴平行于各向异性介质的倾斜对称轴方向，得到新观测坐标系；
基于设置的三维ATI模型的弹性参数和Thomsen参数，计算三维VTI介质条件下的刚度矩阵，建立三维VTI介质弹性波动方程，进行波场延拓；
在新观测坐标系下，采用三维VTI介质弹性波动方程描述三维ATI介质，确定空间域的三维ATI分解算子；
在每一次时间迭代过程中，对弹性波场进行正反傅里叶变换，并基于设置的三维ATI模型的弹性参数和Thomsen参数，得到空间域的辅助波场；
根据空间域的三维ATI分解算子和辅助波场，进行弹性波场分解，得到P波场和S波场。
进一步地，所述三维ATI模型的弹性参数包括纵波速度vp、横波速度vs和地层密度ρ，所述Thomsen参数包括纵波各向异性ε、纵横波各向异性过度参数δ和横波各向异性γ。
进一步地，所述刚度矩阵和Thomson参数之间的关系为：
X33＝ρvp 2
C44＝ρvs 2
C11＝(1+2ε)ρvp 2
C66＝(1+2γ)ρvs 2
Figure BDA0003792694940000021
其中，C11、C13、C33、C44、C66为刚度矩阵中的独立参数；ρ为密度；vp、vs为沿着介质各向异性对称轴方向的P波速度和S波速度。
进一步地，所述空间域的三维ATI分解算子
Figure BDA0003792694940000022
为：
Figure BDA0003792694940000023
其中，
Figure BDA0003792694940000031
分别为对x轴、y轴和z轴的空间一阶偏导；r为系数；θ为倾斜对称轴在新观测坐标系中x-z平面内沿z轴顺时针方向的夹角；
Figure BDA0003792694940000032
为倾斜对称轴在新观测坐标系中x-y平面内的投影沿x轴顺时针方向的夹角。
进一步地，所述在每一次时间迭代过程中，对弹性波场进行正反傅里叶变换，并基于设置的三维ATI模型的弹性参数和Thomsen参数，得到空间域的辅助波场，包括：
根据得到的三维ATI分解算子，构建满足三维ATI介质的泊松方程；
在每一次时间迭代过程中，对弹性波场进行正反傅里叶变换，并基于设置的三维ATI模型的弹性参数和Thomsen参数以及构建的满足三维ATI介质的泊松方程，计算空间域的辅助波场。
进一步地，所述空间域的辅助波场w为：
Figure BDA0003792694940000033
其中，R为系数，u1、u2、u3、u4、u5、u6和u7分别为弹性波场u进行傅里叶变换后得到的波数域形式U构建的参数U1、U2、U3、U4、U5、U6和U7分别进行反傅里叶变换后得到的对应的空间域形式。
进一步地，所述P波场up和S波场us为：
Figure BDA0003792694940000034
Figure BDA0003792694940000035
第二方面，提供一种用于三维ATI介质的弹性波场矢量分解系统，包括：
坐标系旋转模块，用于对观测坐标系进行旋转，使得观测坐标系的垂向坐标轴平行于各向异性介质的倾斜对称轴方向，得到新观测坐标系；
刚度矩阵计算模块，用于基于设置的三维ATI模型的弹性参数和Thomsen参数，计算三维VTI介质条件下的刚度矩阵，建立三维VTI介质弹性波动方程，进行波场延拓；
三维ATI介质描述模块，用于在新观测坐标系下，采用三维VTI介质弹性波动方程描述三维ATI介质，确定空间域的三维ATI分解算子；
辅助波场确定模块，用于在每一次时间迭代过程中，对弹性波场进行正反傅里叶变换，并基于设置的三维ATI模型的弹性参数和Thomsen参数，得到空间域的辅助波场；
弹性波场分解模块，用于根据空间域的三维ATI分解算子和辅助波场，进行弹性波场分解，得到P波场和S波场。
第三方面，提供一种处理设备，包括计算机程序指令，其中，所述计算机程序指令被处理设备执行时用于实现上述用于三维ATI介质的弹性波场矢量分解方法对应的步骤。
第四方面，提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序指令，其中，所述计算机程序指令被处理器执行时用于实现上述用于三维ATI介质的弹性波场矢量分解方法对应的步骤。
本发明由于采取以上技术方案，其具有以下优点：
1、在复杂三维各向异性介质中，本发明同时考虑各向异性以及倾角和方位角的影响，基于特征分析推导得到适用于三维ATI介质的矢量波场分解，通过数值实例的结果，说明本发明能够有效得到分解干净的P波场和S波场。
2、本发明构建泊松方程后，能够避免直接通过求解泊松方程计算三维ATI介质的辅助波场，而采用一种新的方式，大大提高了计算效率以及在三维ATI介质中弹性波场矢量分解的应用。
3、本发明通过给定不同的倾角和方位角参数，带入推导的三维ATI分解算子的公式，适用于大多数的TI介质，也可以退化到三维各向同性的情况，根据介质的弹性参数对本构方程的影响，针对不同弹性介质特征进行三维弹性波分解计算。
4、本发明可以应用在其它的地球物理数据处理工作中，例如分离经过倾斜薄互层间多次反射的P波和S波，以及分离经过裂隙发育页岩气储层的混合弹性波等等，应用领域及其广泛。
综上所述，本发明可以广泛应用于地球物理数据处理和弹性参数反演领域中。
附图说明
通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。在整个附图中，用相同的附图标记表示相同的部件。在附图中：
图1是本发明一实施例提供的方法流程示意图。
具体实施方式
下面将参照附图更详细地描述本发明的示例性实施方式。虽然附图中显示了本发明的示例性实施方式，然而应当理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了能够更透彻地理解本发明，并且能够将本发明的范围完整地传达给本领域的技术人员。
应理解的是，文中使用的术语仅出于描述特定示例实施方式的目的，而无意于进行限制。除非上下文另外明确地指出，否则如文中使用的单数形式“一”、“一个”以及“所述”也可以表示包括复数形式。术语“包括”、“包含”、“含有”以及“具有”是包含性的，并且因此指明所陈述的特征、步骤、操作、元件和/或部件的存在，但并不排除存在或者添加一个或多个其它特征、步骤、操作、元件、部件、和/或它们的组合。文中描述的方法步骤、过程、以及操作不解释为必须要求它们以所描述或说明的特定顺序执行，除非明确指出执行顺序。还应当理解，可以使用另外或者替代的步骤。
尽管可以在文中使用术语第一、第二、第三等来描述多个元件、部件、区域、层和/或部段，但是，这些元件、部件、区域、层和/或部段不应被这些术语所限制。这些术语可以仅用来将一个元件、部件、区域、层或部段与另一区域、层或部段区分开。除非上下文明确地指出，否则诸如“第一”、“第二”之类的术语以及其它数字术语在文中使用时并不暗示顺序或者次序。因此，以下讨论的第一元件、部件、区域、层或部段在不脱离示例实施方式的教导的情况下可以被称作第二元件、部件、区域、层或部段。
需要说明的是：本发明实施例中提到的P波为一种纵波，粒子振动方向和波前进方平行，在所有地震波中，前进速度最快，也最早抵达。P波能在固体、液体或气体中传递。S波为一种横波，前进速度仅次于P波，粒子振动方向垂直于波的前进方向。
本发明实施例提供的用于三维ATI介质的弹性波场矢量分解方法及系统，考虑了复杂各向异性参数以及对称轴倾角和方位角，能够高效率地获得弹性波场分解得到的P波场和S波场。
实施例1
如图1所示，本实施例提供一种用于三维ATI介质的弹性波场矢量分解方法，包括以下步骤：
1)根据观测坐标系下三维ATI介质的倾角和方位角信息，对观测坐标系进行旋转，使得观测坐标系的垂向坐标轴平行于各向异性介质的倾斜对称轴方向，得到新观测坐标系。
具体地，此时，原观测坐标系下的三维ATI介质可看为新观测坐标系下的三维VTI介质。
2)基于设置的三维ATI模型的弹性参数和Thomsen参数，计算三维VTI介质条件下的刚度矩阵Cijkl，建立三维VTI介质弹性波动方程，进行波场延拓，具体为：
2.1)基于PREM地球物理模型(初步地球参考模型)，设置三维ATI模型(任意倾向横向各向同性模型)的弹性参数和Thomsen参数(汤姆森参数)。
具体地，三维ATI模型的弹性参数包括纵波速度vp、横波速度vs和地层密度ρ，Thomsen参数包括纵波各向异性ε、纵横波各向异性过度参数δ和横波各向异性γ。
具体地，本发明可以参照地球物理勘探中熟知的PREM地球物理模型来设置参数，也不仅限于该模型，例如可以参照智利和阿根廷中部区域的远震p波层析成像模型或iasp91地球物理模型。
2.2)基于设置的三维ATI模型的弹性参数和Thomsen参数，计算三维VTI介质条件下的刚度矩阵Cijkl。
具体地，根据弹性张量的Voigt对称性，四阶张量Cijkl中的下角标ij(kl)＝11、22、33、12、23、31可分别对应为i(j)＝1、2、3、4、5、6，于是原本81个独立常数可降为36个，记为刚度矩阵Cijkl。
刚度矩阵和Thomson参数之间的关系为：
C33＝ρvp 2 (1)
C44＝ρvs 2 (2)
C11＝(1+2ε)ρvp 2 (3)
C66＝(1+2γ)ρvs 2 (4)
Figure BDA0003792694940000061
其中，C11、C13、C33、C44、C66为刚度矩阵中的独立参数；ρ为密度；vp、vs为沿着介质各向异性对称轴方向的P波速度和S波速度。
2.3)根据计算的刚度矩阵Cijkl，建立三维VTI介质弹性波动方程，进行波场延拓，模拟波场传播的过程。
具体地，将计算的刚度矩阵Cijkl带入三维VTI介质弹性波动方程中建立针对该模型的弹性波动方程，输入震源的三分量位移数据ux、uy、uz，从零时刻开始，模拟波场传播的过程。
建立的弹性波动方程为：
Figure BDA0003792694940000062
Figure BDA0003792694940000063
Figure BDA0003792694940000064
其中，ux、uy、uz为三分量位移；C11、C13、C33、C44、C66为刚度矩阵中的独立参数；
Figure BDA0003792694940000065
分别为二阶和混合空间偏导。
3)在新观测坐标系下，采用三维VTI介质弹性波动方程描述三维ATI介质，确定空间域的三维ATI分解算子。
具体地，在新观测坐标系下，采用三维VTI介质弹性波动方程描述三维ATI介质，针对三维ATI介质，定义θ为倾斜对称轴在新观测坐标系中x-z平面内沿z轴顺时针方向的夹角，
Figure BDA0003792694940000071
为倾斜对称轴在新观测坐标系中x-y平面内的投影沿x轴顺时针方向的夹角，则空间域的三维ATI分解算子
Figure BDA0003792694940000072
为：
Figure BDA0003792694940000073
其中，
Figure BDA0003792694940000074
分别为对x轴、y轴和z轴的空间一阶偏导；r为系数，且：
Figure BDA0003792694940000075
4)在每一次时间迭代过程中，对弹性波场进行正反傅里叶变换，并基于设置的三维ATI模型的弹性参数和Thomsen参数，得到空间域的辅助波场，具体为：
4.1)基于亥姆霍兹定理，根据得到的三维ATI分解算子
Figure BDA0003792694940000076
构建满足三维ATI介质的泊松方程：
Figure BDA0003792694940000077
其中，u为弹性波场，w为辅助波场。
具体地，根据亥姆霍兹定理，任何弹性波场均可分解为一个无源场和一个无旋场，即
Figure BDA0003792694940000078
若已知三维ATI分解算子
Figure BDA0003792694940000079
则需要通过计算一个辅助波场w，进行弹性波场分解得到分解的P波场和S波场。
4.2)在每一次时间迭代过程中，对弹性波场u进行傅里叶变换，得到弹性波场u的波数域形式U：
U＝-(r1kx 2+r2ky 2+r3kz 2+r4kxky+r5kxkz+r6kykz)W (12)
其中，W为辅助波场w的频率域表示形式，r1、r2、r3、r4、r5、r6均为系数，且：
Figure BDA00037926949400000710
Figure BDA00037926949400000711
r3＝sin2θ+r2cos2θ (15)
r4＝2 sinθcosθ(1-cos2θ-r2sin2θ) (16)
Figure BDA00037926949400000712
Figure BDA00037926949400000713
4.3)根据弹性波场u的波数域形式U，构建参数U1、U2、U3、U4、U5、U6和U7：
Figure BDA00037926949400000714
Figure BDA0003792694940000081
Figure BDA0003792694940000082
Figure BDA0003792694940000083
Figure BDA0003792694940000084
Figure BDA0003792694940000085
Figure BDA0003792694940000086
其中，kx、ky和kz分别为x、y和z方向的波数。
4.4)对参数U1、U2、U3、U4、U5、U6和U7进行反傅里叶变换，得到对应的空间域形式u1、u2、u3、u4、u5、u6和u7。
4.5)基于设置的三维ATI模型的弹性参数和Thomsen参数，根据构建的三维ATI介质的泊松方程以及得到的参数U1、U2、U3、U4、U5、U6和U7及其对应的空间域形式u1、u2、u3、u4、u5、u6和u7，计算空间域的辅助波场w。
具体地，为避免直接求解该泊松方程，提高计算效率，采用快速算法确定辅助波场w，在空间域的辅助波场w的计算公式为：
Figure BDA0003792694940000087
其中，u1、u2、u3、u4、u5、u6、u7分别为参数U1、U2、U3、U4、U5、U6、U7的时间空间域表示形式；R为弹性系数。且：
Figure BDA0003792694940000088
5)采用矢量分解法，根据得到的空间域的三维ATI分解算子
Figure BDA0003792694940000089
和空间域的辅助波场w，进行弹性波场分解，得到P波场up和S波场us：
Figure BDA00037926949400000810
Figure BDA00037926949400000811
其中，·表示点乘积；×表示叉乘积。
上述步骤3)的具体推导过程为：
将三维VTI介质弹性波动方程变换到频率波数域写为矩阵形式：
Figure BDA00037926949400000812
Figure BDA0003792694940000091
其中，
Figure BDA0003792694940000092
为三分量位移ux、uy、uz傅里叶变换后的形式；ω为角频率；kx、ky、kz为三个对称轴方向的波数。
在频散关系
Figure BDA0003792694940000093
(v为相速度)的基础上，将等式(30)两边同除波数的模
Figure BDA0003792694940000094
后，简写为：
Figure BDA0003792694940000095
其中，A为公式(30)中左侧的3×3矩阵，
Figure BDA0003792694940000096
为列向量。
对公式(31)进行特征分析，矩阵A的特征值为ρω2，可以表示不同波的相速度
Figure BDA0003792694940000097
矩阵A的特征向量为
Figure BDA0003792694940000098
表示波的偏振方向。
因此在空间域的偏振方向公式为：
Figure BDA0003792694940000099
Figure BDA00037926949400000910
Figure BDA00037926949400000911
其中，a1、a2、a3分别为P波偏振方向、SH波偏振方向和SV波偏振方向。
三维VTI分解算子公式为：
Figure BDA00037926949400000912
针对三维ATI介质，将观测坐标系进行旋转，使得旋转后的新观测坐标系的z轴沿倾斜对称轴方向，在三维空间下的坐标旋转公式为：
Figure BDA00037926949400000913
构建新观测坐标系下的空间偏导数
Figure BDA0003792694940000101
因此得到三维任意倾向TI介质的分解算子
Figure BDA0003792694940000102
计算公式为上述公式(9)。
上述步骤4)的具体推导过程为：
根据泊松方程求解空间域的辅助波场w，对上述公式(11)展开并进行傅里叶变换至频率波数域得到上述公式(12)。
令系数
Figure BDA0003792694940000103
因此上述公式(12)简化为：
-f(ε,δ)U＝W (37)
基于弱各向异性介质的假设，将系数f(ε,δ)在ε＝0和δ＝0处进行Taylor(泰勒)展开为：
Figure BDA0003792694940000104
其中，R、E均为弹性系数，且弹性系数R为上述公式(27)，弹性系数E为：
Figure BDA0003792694940000105
将系数f(ε,δ)带入上述公式(37)整理得到：
Figure BDA0003792694940000106
代表了频率域中与不同波数进行计算后的弹性波场，最后变换至空间域，得到辅助波场w的计算公式(11)。
实施例2
本实施例提供一种用于三维ATI介质的弹性波场矢量分解系统，包括：
坐标系旋转模块，用于根据观测坐标系下三维ATI介质的倾角和方位角信息，对观测坐标系进行旋转，使得观测坐标系的垂向坐标轴平行于各向异性介质的倾斜对称轴方向，得到新观测坐标系。
刚度矩阵计算模块，用于基于设置的三维ATI模型的弹性参数和Thomsen参数，计算三维VTI介质条件下的刚度矩阵，建立三维VTI介质弹性波动方程，进行波场延拓。
三维ATI介质描述模块，用于在新观测坐标系下，采用三维VTI介质弹性波动方程描述三维ATI介质，确定空间域的三维ATI分解算子。
辅助波场确定模块，用于在每一次时间迭代过程中，对弹性波场进行正反傅里叶变换，并基于设置的三维ATI模型的弹性参数和Thomsen参数，得到空间域的辅助波场。
弹性波场分解模块，用于根据空间域的三维ATI分解算子和辅助波场，进行弹性波场分解，得到P波场和S波场。
实施例3
本实施例提供一种与本实施例1所提供的用于三维ATI介质的弹性波场矢量分解方法对应的处理设备，处理设备可以是用于客户端的处理设备，例如手机、笔记本电脑、平板电脑、台式机电脑等，以执行实施例1的方法。
所述处理设备包括处理器、存储器、通信接口和总线，处理器、存储器和通信接口通过总线连接，以完成相互间的通信。存储器中存储有可在处理设备上运行的计算机程序，处理设备运行计算机程序时执行本实施例1所提供的用于三维ATI介质的弹性波场矢量分解方法。
在一些实现中，存储器可以是高速随机存取存储器(RAM：Random AccessMemory)，也可能还包括非不稳定的存储器(non-volatile memory)，例如至少一个磁盘存储器。
在另一些实现中，处理器可以为中央处理器(CPU)、数字信号处理器(DSP)等各种类型通用处理器，在此不做限定。
实施例4
本实施例提供一种与本实施例1所提供的用于三维ATI介质的弹性波场矢量分解方法对应的计算机程序产品，计算机程序产品可以包括计算机可读存储介质，其上载有用于执行本实施例1所述的用于三维ATI介质的弹性波场矢量分解方法的计算机可读程序指令。
计算机可读存储介质可以是保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是但不限于电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意组合。
上述各实施例仅用于说明本发明，其中各部件的结构、连接方式和制作工艺等都是可以有所变化的，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。

Claims (10)

1.一种用于三维ATI介质的弹性波场矢量分解方法，其特征在于，包括：
对观测坐标系进行旋转，使得观测坐标系的垂向坐标轴平行于各向异性介质的倾斜对称轴方向，得到新观测坐标系；
基于设置的三维ATI模型的弹性参数和Thomsen参数，计算三维VTI介质条件下的刚度矩阵，建立三维VTI介质弹性波动方程，进行波场延拓；
在新观测坐标系下，采用三维VTI介质弹性波动方程描述三维ATI介质，确定空间域的三维ATI分解算子；
在每一次时间迭代过程中，对弹性波场进行正反傅里叶变换，并基于设置的三维ATI模型的弹性参数和Thomsen参数，得到空间域的辅助波场；
根据空间域的三维ATI分解算子和辅助波场，进行弹性波场分解，得到P波场和S波场。
2.如权利要求1所述的一种用于三维ATI介质的弹性波场矢量分解方法，其特征在于，所述三维ATI模型的弹性参数包括纵波速度vp、横波速度vs和地层密度ρ，所述Thomsen参数包括纵波各向异性ε、纵横波各向异性过度参数δ和横波各向异性γ。
3.如权利要求2所述的一种用于三维ATI介质的弹性波场矢量分解方法，其特征在于，所述刚度矩阵和Thomson参数之间的关系为：
C33＝ρvp 2
C44＝ρvs 2
C11＝(1+2ε)ρvp 2
C66＝(1+2γ)ρvs 2
Figure FDA0003792694930000011
其中，C11、C13、C33、C44、C66为刚度矩阵中的独立参数；ρ为密度；vp、vs为沿着介质各向异性对称轴方向的P波速度和S波速度。
4.如权利要求2所述的一种用于三维ATI介质的弹性波场矢量分解方法，其特征在于，所述空间域的三维ATI分解算子
Figure FDA0003792694930000012
为：
Figure FDA0003792694930000013
其中，
Figure FDA0003792694930000014
分别为对x轴、y轴和z轴的空间一阶偏导；r为系数；θ为倾斜对称轴在新观测坐标系中x-z平面内沿z轴顺时针方向的夹角；
Figure FDA0003792694930000015
为倾斜对称轴在新观测坐标系中x-y平面内的投影沿x轴顺时针方向的夹角。
5.如权利要求4所述的一种用于三维ATI介质的弹性波场矢量分解方法，其特征在于，所述在每一次时间迭代过程中，对弹性波场进行正反傅里叶变换，并基于设置的三维ATI模型的弹性参数和Thomsen参数，得到空间域的辅助波场，包括：
根据得到的三维ATI分解算子，构建满足三维ATI介质的泊松方程；
在每一次时间迭代过程中，对弹性波场进行正反傅里叶变换，并基于设置的三维ATI模型的弹性参数和Thomsen参数以及构建的满足三维ATI介质的泊松方程，计算空间域的辅助波场。
6.如权利要求5所述的一种用于三维ATI介质的弹性波场矢量分解方法，其特征在于，所述空间域的辅助波场w为：
Figure FDA0003792694930000021
其中，R为系数，u1、u2、u3、u4、u5、u6和u7分别为弹性波场u进行傅里叶变换后得到的波数域形式U构建的参数U1、U2、U3、U4、U5、U6和U7分别进行反傅里叶变换后得到的对应的空间域形式。
7.如权利要求6所述的一种用于三维ATI介质的弹性波场矢量分解方法，其特征在于，所述P波场up和S波场us为：
Figure FDA0003792694930000022
Figure FDA0003792694930000023
8.一种用于三维ATI介质的弹性波场矢量分解系统，其特征在于，包括：
坐标系旋转模块，用于对观测坐标系进行旋转，使得观测坐标系的垂向坐标轴平行于各向异性介质的倾斜对称轴方向，得到新观测坐标系；
刚度矩阵计算模块，用于基于设置的三维ATI模型的弹性参数和Thomsen参数，计算三维VTI介质条件下的刚度矩阵，建立三维VTI介质弹性波动方程，进行波场延拓；
三维ATI介质描述模块，用于在新观测坐标系下，采用三维VTI介质弹性波动方程描述三维ATI介质，确定空间域的三维ATI分解算子；
辅助波场确定模块，用于在每一次时间迭代过程中，对弹性波场进行正反傅里叶变换，并基于设置的三维ATI模型的弹性参数和Thomsen参数，得到空间域的辅助波场；
弹性波场分解模块，用于根据空间域的三维ATI分解算子和辅助波场，进行弹性波场分解，得到P波场和S波场。
9.一种处理设备，其特征在于，包括计算机程序指令，其中，所述计算机程序指令被处理设备执行时用于实现权利要求1-7中任一项所述的用于三维ATI介质的弹性波场矢量分解方法对应的步骤。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序指令，其中，所述计算机程序指令被处理器执行时用于实现权利要求1-7中任一项所述的用于三维ATI介质的弹性波场矢量分解方法对应的步骤。
CN202210960164.3A 2022-08-11 2022-08-11 一种用于三维ati介质的弹性波场矢量分解方法及系统 Active CN115327626B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210960164.3A CN115327626B (zh) 2022-08-11 2022-08-11 一种用于三维ati介质的弹性波场矢量分解方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210960164.3A CN115327626B (zh) 2022-08-11 2022-08-11 一种用于三维ati介质的弹性波场矢量分解方法及系统

Publications (2)

Publication Number Publication Date
CN115327626A true CN115327626A (zh) 2022-11-11
CN115327626B CN115327626B (zh) 2024-06-18

Family

ID=83922251

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210960164.3A Active CN115327626B (zh) 2022-08-11 2022-08-11 一种用于三维ati介质的弹性波场矢量分解方法及系统

Country Status (1)

Country Link
CN (1) CN115327626B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117910279A (zh) * 2024-03-18 2024-04-19 西北工业大学 基于空间坐标旋转变换的三维声场快速处理方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NO20031618D0 (no) * 2002-04-10 2003-04-09 Schlumberger Technology Bv Fremgangsmåte og apparat for anisotropisk vektorplan bölgedekomposisjon for3D vertikale seismiske profildata
CN111158047A (zh) * 2020-03-04 2020-05-15 中国石油大学(北京) 一种三维弹性波场矢量分解法、装置及计算机存储介质

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NO20031618D0 (no) * 2002-04-10 2003-04-09 Schlumberger Technology Bv Fremgangsmåte og apparat for anisotropisk vektorplan bölgedekomposisjon for3D vertikale seismiske profildata
CN111158047A (zh) * 2020-03-04 2020-05-15 中国石油大学(北京) 一种三维弹性波场矢量分解法、装置及计算机存储介质

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
姚振岸;孙成禹;邓小凡;伍敦仕;: "任意空间取向TI介质弹性波速度特征分析", 石油地球物理勘探, no. 04, 15 August 2017 (2017-08-15) *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117910279A (zh) * 2024-03-18 2024-04-19 西北工业大学 基于空间坐标旋转变换的三维声场快速处理方法
CN117910279B (zh) * 2024-03-18 2024-05-28 西北工业大学 基于空间坐标旋转变换的三维声场快速处理方法

Also Published As

Publication number Publication date
CN115327626B (zh) 2024-06-18

Similar Documents

Publication Publication Date Title
Xu et al. Accurate simulations of pure quasi-P-waves in complex anisotropic media
Min et al. Improved frequency-domain elastic wave modeling using weighted-averaging difference operators
Song et al. Modeling of pseudoacoustic P-waves in orthorhombic media with a low-rank approximation
Jo et al. An optimal 9-point, finite-difference, frequency-space, 2-D scalar wave extrapolator
Yan et al. Elastic wave-mode separation for VTI media
Zhang et al. A stable TTI reverse time migration and its implementation
EP3028071B1 (en) Method and device for the generation and application of anisotropic elastic parameters in horizontal transverse isotropic (hti) media
US9341729B2 (en) Amplitude contrast seismic attribute
CN110058303B (zh) 声波各向异性逆时偏移混合方法
EP3772665A1 (en) Methods and devices using effective elastic parameter values for anisotropic media
Thiel et al. Comparison of acoustic and elastic full‐waveform inversion of 2D towed‐streamer data in the presence of salt
CN111158047B (zh) 一种三维弹性波场矢量分解法、装置及计算机存储介质
Carcione et al. Numerical test of the Schoenberg-Muir theory
Geng et al. Gabor‐frame‐based Gaussian packet migration
CN115327626B (zh) 一种用于三维ati介质的弹性波场矢量分解方法及系统
Wang et al. Elastic wave‐mode separation in 2D transversely isotropic media using optical flow
Zhong et al. Elastic reverse time migration method in vertical transversely isotropic media including surface topography
Raknes et al. Challenges and solutions for performing 3D time-domain elastic full-waveform inversion
CA2693247C (en) Seismic attributes for structural analysis
CN116755148A (zh) 正交各向异性介质多方位反射波走时反演方法
US20190302289A1 (en) Wavefield Propagator for Tilted Orthorhombic Media
Wang et al. Dipole sonic response in deviated boreholes penetrating an anisotropic formation
Baev Imaging of layered media in inverse scattering problems for an acoustic wave equation
CN112698400A (zh) 反演方法、反演装置、计算机设备和计算机可读存储介质
Shin et al. Frequency‐domain waveform modelling and inversion for coupled media using a symmetric impedance matrix

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant