CN102105900A - 在各向异性介质中传播伪声准p波的方法 - Google Patents
在各向异性介质中传播伪声准p波的方法 Download PDFInfo
- Publication number
- CN102105900A CN102105900A CN2009801290649A CN200980129064A CN102105900A CN 102105900 A CN102105900 A CN 102105900A CN 2009801290649 A CN2009801290649 A CN 2009801290649A CN 200980129064 A CN200980129064 A CN 200980129064A CN 102105900 A CN102105900 A CN 102105900A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mtd
- mfrac
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 230000001902 propagating effect Effects 0.000 title claims abstract description 30
- 238000009825 accumulation Methods 0.000 claims abstract description 6
- 238000013508 migration Methods 0.000 claims description 26
- 230000005012 migration Effects 0.000 claims description 26
- 230000006870 function Effects 0.000 claims description 16
- 238000003384 imaging method Methods 0.000 claims description 10
- 230000000644 propagated effect Effects 0.000 claims description 7
- 238000013500 data storage Methods 0.000 claims description 3
- 230000003111 delayed effect Effects 0.000 claims description 2
- 238000005286 illumination Methods 0.000 claims description 2
- 238000010606 normalization Methods 0.000 claims description 2
- 238000004088 simulation Methods 0.000 description 20
- 208000037918 transfusion-transmitted disease Diseases 0.000 description 11
- 238000010586 diagram Methods 0.000 description 6
- 239000006185 dispersion Substances 0.000 description 6
- 230000005284 excitation Effects 0.000 description 6
- 238000009826 distribution Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 2
- 240000004760 Pimpinella anisum Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 229910003460 diamond Inorganic materials 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 238000007598 dipping method Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000012432 intermediate storage Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000004613 tight binding model Methods 0.000 description 1
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
-
- 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/67—Wave propagation modeling
- G01V2210/675—Wave equation; Green's functions
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
一种用于在具有可变倾斜的各向异性介质中保持稳定并且不限于弱各向异性条件的伪声准P波传播的计算机实现方法。该方法包括获取感兴趣的地下区域的地震勘探体积,并且确定所述地震勘探体积的模拟几何体系。该方法进一步包括将所述模拟几何体系用于初始条件来传播至少一个波场通过所述地震勘探体积,防止沿着所述地震勘探体积的对称轴累积能量,以及通过利用有限准S波速度来保证应力-应变关系中的正刚度系数,从而产生稳定波场。该方法包括利用所述稳定波场来生成感兴趣的地下区域的地下图像。
Description
技术领域
本发明一般涉及使用地震信号的地球物理勘探,尤其涉及在可变倾斜各向异性介质中传播伪声准P波并将传播的波场用于地下特性表征的方法。
背景技术
由于与应力有关的矿物和缺陷的择优排序,在许多油气勘探区域(例如,墨西哥湾、北海、和西非海岸)中处处可以观察到各向异性。在这些区域中,往往可以将岩石特性表征成对称轴垂直或倾斜的横向各向同性(“TI”)介质。各向异性介质中的波传播呈现与各向同性介质中的波传播不同的运动学和动力学,因此,需要各向异性模拟和偏移方法来适当地成像储层用于油气勘探。
然而,三维(“3D”)各向异性地震模拟和偏移是计算强度大的任务。与完全弹性的现有技术解决方案相比,基于分散关系的模拟和偏移是计算效率高的替代品。在一种现有技术的方法(Alkhalifah(2000))中,引入了垂直横向各向同性(“VTI”)介质的伪声近似。在那个现有技术方法的近似中,沿着垂直对称轴将剪切波的相速设置成零。如Grechka等人(2004)所述,这种简化不会消除其它方向上的剪切波。根据Alkhalifah的近似,已经为VTI介质中的地震模拟和偏移提出了若干时空域伪声偏微分方程(PDE)(Alkhalifah,2000;Zhou等人,2006;以及Du等人,2008)。PDE的这些系统在运动学上非常接近于涉及矢量场的完全弹性解决方案。
作为VTI介质的一种扩展,正如在与背斜结构和/或逆冲岩席相关联的区域中观察到的那样,可以倾斜TI介质的对称轴(“TTI”)。Zhou等人(2006)通过应用绕对称轴的旋转,将他们的VTI伪声方程推广到用于2D TTI介质的系统。因此,准SV波的相速在与倾斜轴平行或垂直的方向上是零。Lesage等人根据相同相速近似,进一步将Zhou的TTI系统从2D推广到3D。然而,由于倾斜和/或某些岩石特性的迅速横向变化(当垂直速度大于水平速度时),这些现有技术的伪声模拟和偏移方法可能在数值上变得不稳定,并且导致不稳定的波传播。
本领域的普通技术人员应该懂得,各向同性介质中的平面波极化矢量与缓慢矢量平行(P波)或正交(S波)。除了特定传播方向以外,在各向异性介质中没有纯粹的纵波和剪切波。由于这种原因,在各向异性波理论中,快模往往被称为“准P”波,而慢模往往被称为“准S1”和“准S2”。
发明内容
本发明为各向异性介质提供了伪声模拟方法和伪声偏移方法。本发明实施例的一些方面包括一种用于在可变倾斜各向异性介质中保持稳定并且不限于弱各向异性条件的伪声准P波传播的计算机实现方法。该方法还包括获取感兴趣的地下区域的地震勘探体积,并且确定所述地震勘探体积的模拟几何体系。该方法进一步包括:将所述模拟几何体系用于初始条件来传播至少一个波场通过所述地震勘探体积,防止沿着对称轴累积能量,以及通过利用小的有限准S波速度来保证应力-应变关系中的正刚度系数,从而产生稳定波场。该方法还包括利用所述稳定波场来生成感兴趣的地下区域的地下图像。
本发明的另一个实施例包括一种地球物理地震偏移方法,包含下列步骤:建立与地震勘探体积相对应的地震数据集和速度/各向异性模型;以及对于每个共炮/接收器记录,将边界条件设置成包括来自源位置的激励。该实施例进一步包括按照伪声波方程或其等效物前向传播波场:
其中,
Vs0是准SV波的垂直速度,Vp0是准P波的垂直速度,θo是TI介质中对称轴相对于垂直方向的倾斜角,φo是对称轴的方位角,ε、δ是汤姆森各向异性参数,P是标量波场,以及Q是辅助函数。该实施例还包括,对于每个共炮/接收器记录,将边界条件设置成反向传播记录的炮记录,并且按照上面的伪声波方程来反向传播地震数据。该实施例包括应用诸如计算的前向波场和反向波场或它们的等效格林函数之间的互相关(但不限于此)的成像条件来导出地下图像。
本发明的另外实施例还包括如下步骤:通过逆时偏移(RTM)、高斯束偏移、Kirchhoff偏移或其它基于波方程的偏移来传播波场或计算格林函数。
本发明的另外实施例还包括如下步骤:作为选项,除了互相关之外,还应用涉及照明标准化和/或反射角域道集生成和/或相位-振幅补偿的成像条件。
本发明的另外实施例还包括如下步骤:处理共炮/接收器信号,并且在其它相关域中以及以其它模拟和偏移形式来传播波场,所述其它相关域包括(但不限于)共偏置、共方位角和共反射角,所述其它模拟和偏移形式包括(但不限于)延迟炮、平面波和相位编码。
本发明的另外实施例还包括如下步骤:取代汤姆森参数,使用诸如正常时差速度、水平速度的其它等效术语来传播波场或计算格林函数。
本发明的另外实施例包括一种地球物理地震偏移方法,包含如下步骤:建立与地震勘探体积相对应的地震数据集和速度/各向异性模型;以及对于每个共炮/接收器记录,将边界条件设置成包括来自源位置的激励。该实施例还包括按照下列伪声波方程或其等效物来前向传播波场:
该实施例进一步包括:对于每个共炮/接收器记录,将边界条件设置成反向传播记录的炮记录,并且按照上面的伪声波方程反向传播地震数据。该实施例包括应用诸如(但不限于)计算的前向波场和反向波场或它们的等效格林函数之间的互相关的成像条件来导出地下图像。
在地球物理地震偏移中,本发明的不同实施例可以利用其它伪声波方程来前向传播波场。例如,本发明的一个实施例包括按照下面的伪声波方程或其等效物来前向传播波场:
其中,ω是角频率。
用于地球物理地震偏移的本发明的进一步实施例包括按照下面的伪声波方程或其等效物来前向传播波场:
其中,Q和R是辅助函数。
本发明的另一个实施例包括一种地球物理地震偏移方法,包含如下步骤:建立与地震勘探体积相对应的地震数据集和速度/各向异性模型;以及对于每个共炮/接收器记录,将边界条件设置成包括来自源位置的激励。该实施例还包括按照用于倾斜介质的伪声波方程及其等效公式来前向传播波场:
是准P波的垂直速度,是准P波的正常时差速度,η=(ε-δ)/(1+2δ)是Alkhalifah-Tsvankin各向异性参数(按照汤姆森各向异性参数ε和δ表达),P是标量波场,以及U、V、Q和R是辅助函数。该实施例进一步包括:对于每个共炮/接收器记录,将边界条件设置成反向传播记录的炮记录,并且按照上面的伪声波方程来反向传播地震数据。该实施例包括应用诸如计算的前向波场和反向波场或它们的等效格林函数之间的互相关的成像条件来导出地下图像。
用于地球物理地震偏移的本发明的不同实施例可以利用其它伪声波方程为倾斜介质前向传播波场。例如,本发明的一个实施例包括按照用于倾斜介质的伪声波方程及其等效公式来前向传播波场:
其中,
是准P波的垂直速度,是准P波的正常时差速度,a是剪切波与P波速度比的平方,η=(ε-δ)/(1+2δ)是Alkhalifah-Tsvankin各向异性参数(按照汤姆森各向异性参数ε和δ表达),以及Q和R是辅助函数。
用于地球物理地震偏移的本发明的进一步实施例包括:按照用于倾斜介质的伪声波方程及其等效公式或它们的导出公式/等效物来前向传播波场:
其中,F是标量波场。
本发明的另一个实施例包括一种地球物理地震偏移方法,包含如下步骤:建立与地震勘探体积相对应的地震数据集和速度/各向异性模型;以及对于每个共炮/接收器记录,将边界条件设置成包括来自源位置的激励。该实施例还包括按照伪声波方程或其导出公式/等效物来前向传播波场:
该实施例进一步包括:对于每个共炮/接收器记录,将边界条件设置成反向传播记录的炮记录,并且按照上面的伪声波方程反向传播地震数据。该实施例包括应用诸如(但不限于)计算的前向波场和反向波场或它们的等效格林函数之间的互相关的成像条件来导出地下图像。
本发明的另一个实施例包括一种地球物理地震模拟方法,包含如下步骤:建立与地震勘探体积相对应的速度/各向异性模型;以及对于每个炮,设置波场的初始条件。该实施例还包括按照伪声波方程或其等效物来前向传播波场:
其中,w(t)是源函数,以及是源位置的矢量。可以改变源项及其插入形式而不影响支配PDE。
用于地球物理地震模拟的本发明的另一个实施例包括:按照用于倾斜介质的伪声波方程(方程5)及其等效公式来前向传播波场。
用于地球物理地震模拟的本发明的另一个实施例包括:按照用于倾斜介质的伪声波方程(方程6)及其等效公式来前向传播波场。
还应当理解,本发明打算与通常包括电子配置的系统一起使用,所述电子配置包括至少一个处理器、用于存储程序代码或其它数据的至少一个存储设备、可选的视频监视器或其它显示设备(即,液晶显示器)以及至少一个输入设备。所述处理器最好是能够显示图像和处理复杂数学算法的基于微处理器或微控制器的平台。所述存储设备可以包括用于存储在与本发明相关联的特定处理期间生成或使用的事件或其它数据的随机访问存储器(RAM)。所述存储设备还可以包括用于存储用于本发明的控制和处理的程序代码的只读存储器(ROM)。
一个这样的实施例包括一种被配置成进行在可变倾斜各向异性介质中保持稳定并且不限于弱各向异性条件的伪声准P波传播的系统。所述系统包括:数据存储设备,其具有包括感兴趣的地下区域的地震勘探体积的计算机可读数据;以及处理器,其被配置和安排成执行存储在处理器可访问存储器中的机器可执行指令以便执行一种方法。这个特定实施例的方法包括:确定所述地震勘探体积的模拟几何体系;以及将所述模拟几何体系用于初始条件来传播至少一个波场通过所述地震勘探体积,防止沿着所述地震勘探体积内的各向异性区域的对称轴累积能量,以及保证应力-应变关系中的正刚度系数,从而产生稳定波场。该方法进一步包括:利用所述稳定波场来生成感兴趣的地下区域的地下图像。
通过参照形成本说明书的一部分、相同标号在各个图形中表示相应部件的附图对如下描述和所附权利要求书加以研究,本发明的这些和其它目的、特征、和特性,以及操作方法、结构的相关元件的功能、部件的组合、和制造成本将变得更加显而易见。然而,不言而喻,这些附图只用于例示和描述的目的,而无意作为限制本发明的定义。正如用在说明书和权利要求书中的那样,除非上下文另有明确指明,单数形式“一个”、“一种”、和“该”也包括复数指示物。
附图说明
图1是例示依照本发明一个或多个实施例的方法的流程图;
图2是例示依照本发明一个或多个实施例的方法的流程图;
图3是例示依照本发明一个或多个实施例的方法的流程图;
图4例示了按照现有技术的Alkhalifah近似的示范性波传播模拟,其中f=1;
图5例示了按照本发明一个实施例的示范性波传播模拟;
图6例示了按照本发明一个实施例的示范性波传播模拟,其中Vso/Vpo=0.01;
图7例示了按照现有技术的Alkhalifah近似的示范性波传播模拟,其中Vso=0;
图8例示了按照现有技术的Alkhalifah近似的示范性相速分布;
图9例示了按照现有技术的Alkhalifah近似的示范性群速分布;
图10例示了本发明一个实施例的示范性相速分布,其中Vso/Vpo=0.01;
图11例示了本发明一个实施例的示范性群速分布,其中Vso/Vpo=0.01;
图12例示了按照利用一阶5×5PDF系统的本发明一个实施例,具有可变倾斜对称轴的介质中的示范性波传播模拟;
图13例示了用在本发明一个实施例中的几何体系的示意图;以及
图14是执行依照本发明实施例的方法的系统的实施例的示意性例示图。
具体实施方式
本发明的一个实施例例示在图1中,其中,流程图10描述了传播在具有可变倾斜的各向异性介质中保持稳定的准P波的方法。本发明不限于弱各向异性条件。这个特定实施例包括获取感兴趣的地下区域的地震勘探体积12,并且确定地震勘探体积的模拟几何体系14。该实施例进一步包括将模拟几何体系用于初始条件来传播至少一个波场通过所述地震勘探体积,防止沿着所述地震勘探体积的对称轴累积能量,以及利用有限准S波速度来保证应力-应变关系中的正刚度系数,从而产生稳定波场16。然后可以利用所述稳定波场来生成感兴趣的地下区域的地下图像18。
本领域的普通技术人员应该懂得,本发明的不同实施例可以为各向异性介质提供伪声模拟方法或伪声偏移方法。例如,图2例示了用于具有可变倾斜的各向异性介质中的波传播的伪声模拟方法的一个实施例的流程图,其中该方法不局限于弱各向异性条件。该实施例包括获取感兴趣的地下区域的地震勘探体积22,并且确定所述地震勘探体积的模拟几何体系24。该实施例还包括将模拟几何体系用于初始条件来传播至少一个波场通过所述地震勘探体积,其中,沿着所述地震勘探体积的对称轴的人为准剪切波速度大于或等于零,从而防止沿着对称轴累积能量,从而产生稳定波场26。然后可以利用稳定波场来生成感兴趣的地下区域的地下图像28。
图3例示了可以用于伪声偏移的本发明另一个实施例的流程图。该实施例包括获取感兴趣的地下区域的地震勘探体积32,并且确定地震勘探体积的模拟几何体系34。该实施例还包括将模拟几何体系用于初始条件来传播至少一个波场通过地震勘探体积,其中,准剪切波能量不沿着地震勘探体积的对称轴累积,从而产生稳定波场36。然后可以利用稳定波场来生成感兴趣的地下区域的地下图像38。
本发明相对于传统的声各向异性模拟和偏移具有若干优点。本发明在具有可变倾斜的TI介质中提供了稳定的波传播方式,因此可以获得仿真波场传播和反射率图像。现有技术的伪声模拟和偏移方法基于Alkhalifah近似,在所述Alkhalifah近似中,沿着对称轴将剪切波的相速设置成零。尽管现有技术的方法可以用在恒定倾斜的TI介质中,但在倾斜变化可以在对称轴附近局部集中高能量的区域中,零速度剪切波可以使波传播不稳定(即,振幅变成无界)。图4示出了在可变倾斜介质中(例如,在背斜结构的顶部附近)现有技术的方法(f=1)是不稳定的40。相反,图5示出了基于本发明的波传播(f=0.98)在相同介质中保持稳定42。另外,本发明可以提供控制剪切波与P波速度比以优化模拟和偏移的结果的灵活性。例如,可以将剪切波与P波速度比设置成接近实际值,以近似弹性波传播中的运动学。而且,在某些岩石中,相对于对称轴,垂直速度可以大于水平速度。在这样的情况下,基于Alkhalifah近似的波方程将导致负刚度矩阵,从而与数值实现算法无关地产生不稳定波场。本发明可以使用有限剪切波速度来保证应力-应变关系中的正刚度系数,从而生成稳定波传播。
在基于Alkhalifah近似的现有技术方法中,使剪切波相速等于零并不会消除剪切波。而是高能量集中在对称轴附近。唯一例外是剪切波处处都消失的椭圆各向异性(即,ε=δ)。在本发明中,将垂直剪切波速度从等于零中解脱出来,因此,能量较少集中在对称轴附近。由于即使满足椭圆各向异性条件也存在附加交叉导数,所以剪切波将不会消失。
就计算成本而言,本发明的实施例所使用的PDE与现有技术的方法相比涉及要计算附加空间导数项。在具有可变倾斜的区域中,需要与非零Vs0相关联的附加工作负荷以实现地震模拟和偏移所需的稳定性和可靠性。在倾斜几乎恒定或非常平缓的区域中,可以跳过所述附加工作负荷。
对于本领域的普通技术人员来说清楚的是,可以不偏离本发明的范围地以许多方式改变上面的实施例。例如,本领域的普通技术人员可以明显看出,在本发明中,可以将不同初始条件或边界条件或PDE的不同线性组合方便地用在模拟和偏移中。
在本发明的一个实施例中,各向异性模拟方法包括:建立与地震勘探体积相对应的速度和各向异性模型;设置诸如源激励的初始条件;按照方程[1]或其等效物,在具有倾斜对称轴或垂直对称轴的横向各向同性介质中传播波。对于正向模拟,在方程[1]或方程[2]中,需要在方程的右侧引入形式的源函数,其中是源位置,以及w(t)是源小波。
在本发明的上述实施例中,与现有技术方法中将f四舍五入成1的近似相反,方程[3]中的垂直剪切波速度可以非零(因此,f可以不同于1)。因此,在本发明中,剪切波的相速在与对称轴平行和垂直的两个方向上都可以非零。在具有可变倾斜的介质中,准剪切波的有限速度可以避免往往发生在对称轴附近的高能量的局部集中。本发明不需要弱各向异性假设。
利用上面的PDE,可以为各向异性介质导出本发明的其它实施例。如果倾斜角θo=0,那么上面的PDE简化成3D VTI系统,类似地,当φo=0时,简化成3D HTI系统。用于3D VTI介质的二阶3×3系统可以采取方程6中的形式。通过用在方程[1]中给出的f1和f2取代在方程[6]中给出的g1和g2,可将PDE的这种系统推广到用于倾斜TI介质的其等效公式。
作为f=1时的方程[1]或[2]或[4]中的PDE的替代物,方程[5]中的PDE的一阶5×5系统是双曲线的,并且在具有可变倾斜的TI介质中是稳定的。本发明的这个实施例是对称双曲线的(适定的,甚至对于可变系数)。这个系统也可推广到可变倾斜TTI。3D中的PDE的上述完整一阶5×5系统在2D中减小成4×4。
如上所述,本发明的另外实施例还提供了伪声偏移方法。一个实施例包括如下步骤:建立与地震勘探体积相对应的地震数据集和速度/各向异性模型;设置波传播的边界条件;按照方程[1]方程[2]、方程[4]或方程[6]或它们的等效物,在各向异性介质中分开传播来自源激励的波和记录的地震数据;以及应用诸如(但不限于)两个传播波场之间的互相关的成像条件来获取地下图像。可以应用不同的初始和/或边界条件,而不影响本发明的范围。用于传播源小波的示范性边界条件(例如,根据方程[1])如下:
以及用于地震数据的逆时外推的边界条件如下:
其中,w(t’)是源函数,xs是源位置,以及D(x,y,xs,ys;t)是要偏移的炮记录。
下列例子例示了本发明的进一步实施例:
1.为VTI介质中的准P波建立四阶分散关系
vso未设置成零的用于VTI介质的Tsvankin相速关系导致分散关系:
ω4-Bω2+C=0 [12]
其中:
ωqP±对应于准P波;ωqS±对应于准SV波。
2.为VTI介质中的准P波建立四阶PDE
3.为VTI介质建立PDE的二阶3×3系统
设:
其中,F(x,y,z,t)是满足方程[15]的波场。假设初始条件为:
导致:
设:
然后,通过方程[4]让方程[15]等效于PDE的二阶3×3系统。
图6示出了将上面的PDE用于本发明的那个特定实施例的VTI介质中的波前传播。与基于现有技术方法(例示在图7中)的波前相比,外部qP波前(图6中的44和图7中的48)保持几乎相同,但内部qSV波前(图6中的46和图7中的50)具有不同于菱形的形式。图8和图9分别示出了按照Alkhalifah(现有技术)近似的相速和群速。相反,图10和图11分别示出了按照本发明实施例的相速和群速。与Alkhalifah近似相比,qSV波的相速从沿着对称轴等于零中解脱出来。因此,群速的最大值或高能量不像现有技术的方法中那样沿着对称轴集中。通过应用绕倾斜轴的旋转,可对恒定倾斜的TTI介质应用相同观察。
4.为VTI介质建立PDE的一阶5×5系统
Du等人(2008)为vs0=0的VTI介质提出了下列PDE的一阶2×2系统:
其中,g1和g2在方程[6]之后给出,p是标量波场,以及q是辅助函数。新波场P以及新辅助函数U、V、Q和R通过下式定义:
然后方程5是PDE的完整一阶5×5系统。这个系统可以通过对称化而显示成双曲线的。设:
然后:
其中,
在2D中,消除变量V,并删除第三个方程,得出一阶4×4系统。
图12示出了在可变倾斜介质中受这种一阶PDE支配的稳定波前传播。
5.在TTI介质中建立四阶分散关系
通过解脱沿着对称轴Vso=0(或f=1)的Alkhalifah近似,可以从Tsvankin相速关系(2001)中导出下列方程:
其中,θ0是TI介质中对称轴相对于垂直方向的倾斜角,以及φo是对称轴的方位角。认识到:
可以导出下列四阶分散关系:
6.在TTI介质中建立四阶PDE
将上述四阶分散关系的两侧乘以标量波场P,并且将频率-波数算子转换到时空域,TTI/VTI介质的四阶PDE采取方程[8]的形式。
7.为TTI介质建立PDE的二阶2×2系统
上述TTI介质的四阶伪声PDE可以经由方程[3],通过2×2时空域PDE系统来求解:
PDE的2×2系统在水平速度vh和正常时差(NMO)速度vn方面也可以采取等效形式。
在作为特殊情况的2D介质中,通过下列简化偏导数算子,上面的PDE仍然保持有效:
8.为TTI介质建立PDE的二阶3×3系统
作为替代物,也可以使用不同的线性组合,通过方程[4]中的3×3时空域PDE或其等效物来求解TTI介质的四阶伪声PDE。
本发明的实施例可以在诸如现场可编程门阵列(FPGA)、图形处理单元(GPU)、存储单元的协处理器加速架构或通用计算机上实现。本发明提供了编有执行本发明的方法的指令的装置和通用计算机和/或协处理器,以及编码执行本发明的方法的指令的计算机可读介质。在图14中示意性地例示了执行本发明实施例的系统。该系统52包括数据存储设备或存储器54。可以使存储的数据用于诸如可编程通用计算机的处理器56。处理器56可以包括诸如显示器58和图形用户界面60的接口组件。图形用户界面(GUI)可以用于显示数据和处理后的数据产物,以及允许用户在实现该方法的各个方面的选项当中作出选择。数据可以经由总线62直接从数据获取设备或从中间存储设备或处理设施(未示出)传送给系统52。
尽管为了例示的目的,根据当前认为最实用的优选实施例对本发明作了详细描述,但应该明白,这样的细节仅仅为了那个目的,本发明不局限于所公开的实施例,而是相反,打算涵盖在所附权利要求书的精神和范围之内的所有修改和等效安排。例如,尽管本文提到一台计算机,但这台计算机可以包括通用计算机、专用计算机、编程成执行这些方法的ASIC(专用集成电路)、计算机阵列或网络、或其它适用计算设备。作为进一步的例子,还应该明白,本发明设想,可以尽可能地将任何实施例的一个或多个特征与任何其它实施例的一个或多个特征组合。
Claims (15)
1.一种用于在可变倾斜各向异性介质中保持稳定并且不限于弱各向异性条件的伪声准P波传播的计算机实现方法,所述方法包括:
获取感兴趣的地下区域的地震勘探体积;
确定所述地震勘探体积的模拟几何体系;
将所述模拟几何体系用于初始条件来传播至少一个波场通过所述地震勘探体积,防止沿着所述地震勘探体积内的各向异性区域的对称轴累积能量,以及利用有限准S波速度来保证应力-应变关系中的正刚度系数,从而产生稳定波场;以及
利用所述稳定波场来生成感兴趣的地下区域的地下图像。
2.如权利要求1所述的方法,其中,传播至少一个波场通过所述地震勘探体积包括:人为准剪切波速度沿着所述地震勘探体积内的各向异性区域的对称轴大于或等于零。
3.如权利要求1所述的方法,其中,传播至少一个波场通过所述地震勘探体积包括:限制准剪切波沿着所述地震勘探体积的所述对称轴累积。
4.如权利要求1所述的方法,其中,多个波场被传播通过所述地震勘探体积。
5.如权利要求1所述的方法,其中,传播至少一个波场通过所述地震勘探体积包括:利用逆时偏移、基于波方程的偏移、高斯束偏移或Kirchhoff偏移之一。
6.如权利要求1所述的方法,包括:前向和反向传播波场通过所述地震勘探体积,以及将成像条件应用于前向和反向波场或等效格林函数以导出地下图像。
7.如权利要求6所述的方法,其中,应用成像条件的步骤包括:应用前向和反向波场或等效格林函数之间的互相关来导出地下图像。
8.如权利要求7所述的方法,其中,应用成像条件的步骤包括:照明标准化、反射角域道集生成和相位-振幅补偿中的至少一个。
9.如权利要求1所述的方法,其中,成像输出几何体系包括共偏置、共方位角和共反射角域。
10.如权利要求1所述的方法,其中,传播至少一个波场包括:利用延迟炮、平面波和相位编码中的至少一个。
11.如权利要求1所述的方法,其中,传播至少一个波场包括:利用正常时差速度、水平速度和汤姆森参数中的至少一个。
12.一种配置成执行在可变倾斜各向异性介质中保持稳定并且不限于弱各向异性条件的伪声准P波传播的系统,所述系统包含:
数据存储设备,具有包括感兴趣的地下区域的地震勘探体积的计算机可读数据;以及
处理器,被配置和安排成执行存储在处理器可访问的存储器中的机器可执行指令以便执行包含如下步骤的方法:
确定所述地震勘探体积的模拟几何体系;
将所述模拟几何体系用于初始条件来传播至少一个波场通过所述地震勘探体积,防止沿着所述地震勘探体积内的各向异性区域的对称轴累积能量,以及利用有限准S波速度来保证应力-应变关系中的正刚度系数,从而产生稳定波场;以及
利用所述稳定波场来生成感兴趣的地下区域的地下图像。
13.如权利要求12所述的系统,其中,传播至少一个波场通过所述地震勘探体积包括:人为准剪切波速度沿着所述地震勘探体积内的各向异性区域的对称轴大于或等于零。
14.如权利要求12所述的系统,其中,传播至少一个波场通过所述地震勘探体积包括:限制准剪切波沿着所述地震勘探体积的所述对称轴累积。
15.如权利要求1所述的系统,其中,多个波场被传播通过所述地震勘探体积。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US8467208P | 2008-07-30 | 2008-07-30 | |
US61/084,672 | 2008-07-30 | ||
PCT/US2009/050220 WO2010014379A2 (en) | 2008-07-30 | 2009-07-10 | Method for propagating pseudo acoustic quasi-p waves in anisotropic media |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102105900A true CN102105900A (zh) | 2011-06-22 |
CN102105900B CN102105900B (zh) | 2013-03-20 |
Family
ID=41610909
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009801290649A Expired - Fee Related CN102105900B (zh) | 2008-07-30 | 2009-07-10 | 在各向异性介质中传播伪声准p波的方法 |
Country Status (7)
Country | Link |
---|---|
EP (1) | EP2308025A4 (zh) |
CN (1) | CN102105900B (zh) |
AU (1) | AU2009276976B2 (zh) |
BR (1) | BRPI0916637A2 (zh) |
CA (1) | CA2730017C (zh) |
EA (1) | EA201170261A1 (zh) |
WO (1) | WO2010014379A2 (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015149237A1 (en) * | 2014-03-31 | 2015-10-08 | Prad Research And Development Limited | Subsurface formation modeling with integrated stress profiles |
CN106597531A (zh) * | 2015-10-16 | 2017-04-26 | 中国石油化工股份有限公司 | 含垂直裂缝的页岩的波场传播特征的正演模拟方法 |
CN107102353A (zh) * | 2017-05-08 | 2017-08-29 | 厦门大学 | 基于高阶差分方法的弹性波方程逆时偏移成像方法 |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2797434C (en) | 2010-05-12 | 2017-09-19 | Shell Internationale Research Maatschappij B.V. | Seismic p-wave modelling in an inhomogeneous transversely isotropic medium with a tilted symmetry axis |
US8385151B2 (en) * | 2010-06-24 | 2013-02-26 | Chevron U.S.A. Inc. | Reverse time migration with absorbing and random boundaries |
CA2806651A1 (en) | 2010-07-28 | 2012-02-02 | Cggveritas Services Sa | 3-d harmonic-source reverse time migration systems and methods for seismic data analysis |
US9207342B2 (en) * | 2012-03-09 | 2015-12-08 | Chevron U.S.A. Inc. | Correction of shear log for elastic anisotropy |
AU2012387658B2 (en) | 2012-08-17 | 2016-06-23 | Landmark Graphics Corporation | Systems and methods for imaging seismic data |
US20150185347A1 (en) * | 2013-12-30 | 2015-07-02 | Chevron U.S.A. Inc. | System and method of mitigating instabilities in a pseudoacoustic wave propagator |
US10215869B2 (en) | 2015-03-30 | 2019-02-26 | Chevron U.S.A. Inc. | System and method of estimating anisotropy properties of geological formations using a self-adjoint pseudoacoustic wave propagator |
CN111999770B (zh) * | 2020-09-03 | 2024-01-16 | 中国地质科学院 | 一种tti介质转换ps波精确束偏移成像方法及系统 |
CN113341455B (zh) * | 2021-06-24 | 2024-02-09 | 中国石油大学(北京) | 一种粘滞各向异性介质地震波数值模拟方法、装置及设备 |
CN117233838B (zh) * | 2023-09-20 | 2024-04-05 | 长江大学 | 一种二维vti介质中的弹性准纵横波场分离和逆时偏移成像方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR810000683B1 (ko) * | 1976-05-27 | 1981-06-20 | 죠오지프 에스 이미리이 | 지진파 데이타 수집장치 |
GB2278920B (en) | 1993-06-07 | 1996-10-30 | Geco As | Method of determining earth elastic parameters in anistropic media |
GB9906995D0 (en) * | 1998-09-16 | 1999-05-19 | Geco Prakla Uk Ltd | Seismic detection apparatus and related method |
US6819111B2 (en) * | 2002-11-22 | 2004-11-16 | Baker Hughes Incorporated | Method of determining vertical and horizontal resistivity, and relative dip in anisotropic earth formations having an arbitrary electro-magnetic antenna combination and orientation with additional rotation and position measurements |
JP4105139B2 (ja) * | 2004-09-30 | 2008-06-25 | 石油資源開発株式会社 | 3次元画像生成方法及びプログラム |
KR100660562B1 (ko) * | 2005-09-30 | 2006-12-22 | (주)지 텍크 | 지반 표층부 3차원 초고해상도 s파 탐사장치 및 탐사방법 |
US9519072B2 (en) * | 2006-05-11 | 2016-12-13 | Schlumberger Technology Corporation | Method and apparatus for locating gas hydrate |
-
2009
- 2009-07-10 BR BRPI0916637A patent/BRPI0916637A2/pt not_active IP Right Cessation
- 2009-07-10 CN CN2009801290649A patent/CN102105900B/zh not_active Expired - Fee Related
- 2009-07-10 WO PCT/US2009/050220 patent/WO2010014379A2/en active Application Filing
- 2009-07-10 CA CA2730017A patent/CA2730017C/en not_active Expired - Fee Related
- 2009-07-10 AU AU2009276976A patent/AU2009276976B2/en active Active
- 2009-07-10 EA EA201170261A patent/EA201170261A1/ru unknown
- 2009-07-10 EP EP09803352.5A patent/EP2308025A4/en not_active Withdrawn
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015149237A1 (en) * | 2014-03-31 | 2015-10-08 | Prad Research And Development Limited | Subsurface formation modeling with integrated stress profiles |
GB2539592A (en) * | 2014-03-31 | 2016-12-21 | Logined Bv | Subsurface formation modeling with integrated stress profiles |
US10386523B2 (en) | 2014-03-31 | 2019-08-20 | Schlumberger Technology Corporation | Subsurface formation modeling with integrated stress profiles |
GB2539592B (en) * | 2014-03-31 | 2020-12-30 | Logined Bv | Subsurface formation modeling with integrated stress profiles |
CN106597531A (zh) * | 2015-10-16 | 2017-04-26 | 中国石油化工股份有限公司 | 含垂直裂缝的页岩的波场传播特征的正演模拟方法 |
CN106597531B (zh) * | 2015-10-16 | 2019-10-29 | 中国石油化工股份有限公司 | 含垂直裂缝的页岩的波场传播特征的正演模拟方法 |
CN107102353A (zh) * | 2017-05-08 | 2017-08-29 | 厦门大学 | 基于高阶差分方法的弹性波方程逆时偏移成像方法 |
Also Published As
Publication number | Publication date |
---|---|
BRPI0916637A2 (pt) | 2019-03-26 |
AU2009276976A1 (en) | 2010-02-04 |
WO2010014379A2 (en) | 2010-02-04 |
CN102105900B (zh) | 2013-03-20 |
WO2010014379A3 (en) | 2010-03-25 |
EA201170261A1 (ru) | 2011-08-30 |
AU2009276976B2 (en) | 2015-07-16 |
EP2308025A2 (en) | 2011-04-13 |
CA2730017A1 (en) | 2010-02-04 |
EP2308025A4 (en) | 2016-11-16 |
CA2730017C (en) | 2017-03-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102105900A (zh) | 在各向异性介质中传播伪声准p波的方法 | |
US8332156B2 (en) | Method for propagating pseudo acoustic quasi-P waves in anisotropic media | |
US11092707B2 (en) | Determining a component of a wave field | |
US20170299745A1 (en) | Prestack egs migration method for seismic wave multi-component data | |
Yan et al. | Elastic wave-mode separation for VTI media | |
EP2335093B1 (en) | Estimation of soil properties using waveforms of seismic surface waves | |
US9158017B2 (en) | Seismic imaging apparatus utilizing macro-velocity model and method for the same | |
EP3256886B1 (en) | Efficient and stable absorbing boundary condition in finite-difference calculations | |
US20120275267A1 (en) | Seismic Data Processing | |
EP2787372A2 (en) | Device and method for stable least-squares reverse time migration | |
US20100110832A1 (en) | Fourier finite-difference migration for three dimensional tilted transverse isotropic media | |
US20180210101A1 (en) | System and method for seismic inversion | |
US20120263015A1 (en) | Device and method for calculating 3d reverse time migration in tilted orthorhombic media | |
US11199641B2 (en) | Seismic modeling | |
Jia et al. | Superwide-angle one-way wave propagator and its application in imaging steep salt flanks | |
Oren et al. | 3D anisotropic elastic time-reverse imaging of surface-recorded microseismic data | |
EP2184620B1 (en) | Fourier finite-difference migration for three dimensional tilted transverse isotropic media | |
Tang et al. | A fast RTM implementation in TTI media | |
Tang et al. | Fourier finite-difference wave-equation migration in tilted transversely isotropic media with an improved solution for coefficient estimation | |
Yan et al. | 3D elastic wave mode separation for TTI media | |
Fathalian et al. | Viscoacoustic VTI and TTI wave equations and their application for anisotropic reverse time migration: Constant-Q approximation | |
Bakker | A stable one-way wave propagator for VTI media | |
Qin et al. | Pure quasi-P-wave modeling and imaging using an approximated space-domain pseudo-differential operator in the VTI media | |
EP3092520B1 (en) | Determining a component of a wave field in an anisotropic subsurface for migration | |
Zhan | Reverse-time Migration in Tilted Transversely Isotropic Media with Decoupled Equations |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130320 |
|
CF01 | Termination of patent right due to non-payment of annual fee |