CN111936888B - 用于倾斜正交晶介质的波场传播子 - Google Patents

Abstract

系统和方法包括：接收烃储层的储层数据；接收与波场传播子的选择相关的指示；利用傅立叶有限变换和有限差分应用波场传播子，以对与倾斜正交晶介质相关联的波场建模，所述倾斜正交晶介质表示包括所述烃储层的地下区域；以及结合所述波场传播子处理储层数据，以生成用于在包括烃储层并含有结构或地层特征的地下区域上方的地震勘探的输出，所述结构或地层特征有助于烃的呈现、偏移或累积。

Description

用于倾斜正交晶介质的波场传播子

技术领域

本公开大体上涉及分析地震数据，且更具体来说，涉及利用波场传播子来处理地震测量。

背景技术

本章节旨在向读者介绍可能与本公开的各个方面相关的本领域的各个方面，这些方面在下文进行描述和/或要求保护。此讨论被认为有助于为读者提供背景知识，以便于更好地理解本公开的各个方面。因此，应理解，将从这个角度阅读这些陈述，而不是作为对现有技术的承认。

地震勘测包括通过将声能向下发送到地面并且记录从地下区域内的地质层返回的反射声能来生成地球的地下区域的图像或地图。在地震勘测期间，能源置于可能包括烃沉积物的地球表面区域上或上方的各个位置处。每当激活源时，源生成地震(例如，声波)信号，所述地震信号向下行进穿过地球、被反射并在其返回时使用安置于地球的地下区域上或上方的一个或多个接收器进行记录。由接收器记录的地震数据随后可以用于创建对应的地下区域的图像或轮廓。

发明内容

下文阐述本文所公开的某些实施例的概述。应理解，呈现这些方面仅仅是为了向读者提供这些某些实施例的简要概述，并且这些方面无意于限制本公开的范围。实际上，本公开可以包括下文可能未阐述的多个方面。

地震处理可以包括例如地震数据偏移、全波形反演处理、逆时偏移，以及利用地下的精确(以及通常复杂)表示的其它技术。生成的地下的表示可以称为波场传播子(或波传播子)，并且各种技术可以用作波传播方法以生成用于地震处理的不同类型的波场传播子。在一些实施例中，这些传播方法可以包括伪谱(PS)波传播方法、倾斜横向各向同性(TTI)波传播方法(诸如，TTI傅立叶有限差分(FFD)波传播方法)，以及倾斜正交晶FFD(TORFFD)波传播方法(其例如可以是傅立叶有限变换和有限差分的链算子)。如将在下文更详细地论述，使用TORFFD波传播方法可以相对于TTI波传播方法允许包括烃储层并含有结构或地层特征的地下区域的更高分辨率图像，所述结构或地层特征有助于经由处理接收到的地震数据呈现、偏移或累积生成的烃，同时至少改进了使用伪谱波传播方法所需的处理时间。

附图说明

在阅读以下详细描述后并且在参考附图后可以更佳地理解本公开的各个方面，在附图中：

图1是根据本文呈现的实施例的可以基于对经由地震勘测系统获取的地震数据的分析来执行的各个过程的流程图；

图2是根据本文呈现的实施例的海洋环境中的海洋勘测系统的示意图；

图3是根据本文呈现的实施例的海洋环境中的第二海洋勘测系统的示意图；

图4是根据本文呈现的实施例的计算系统的框图，所述计算系统可以基于经由图2的海洋勘测系统和/或图3的第二海洋勘测系统获取的数据来执行本文所描述的操作；

图5是根据本文呈现的实施例的用于经由图4的计算系统生成地震图像的方法的流程图；

图6说明根据本文呈现的实施例的用于由图4的计算系统应用的两个倾斜正交晶建模技术的伪声波场快照；

图7说明根据本文呈现的实施例的由图4的计算系统利用TTI波传播方法生成的第一地震图像的示例；

图8说明根据本文呈现的实施例的由图4的计算系统在倾斜正交晶介质上利用伪谱波传播方法生成的第一地震图像的示例；以及

图9说明根据本文呈现的实施例的由图4的计算系统在倾斜正交晶介质上利用FFD波传播方法生成的第一地震图像的示例。

具体实施方式

下文将描述一个或多个特定实施例。为了提供对这些实施例的简要描述，在说明书中并未描述实际实施方案的所有特征。应理解，在任何此种实际实施方案的开发中，如在任何工程或设计项目中，都必须做出许多特定于实施方案的决策，以实现开发人员的特定目标，诸如遵守可能在实施方案之间变化的系统相关和业务相关的约束。此外，应了解，这种开发工作可能复杂且耗时，但是对于受益于本公开的普通技术人员而言，这将是设计、制作和制造的常规任务。

地震数据可以提供关于描述的有价值信息，诸如，地球的地下区域内的烃沉积物的位置和/或变化。不同处理技术，诸如对于特定选择的介质(例如，对于作为波场传播子的地球的地下区域的选定表示)应用不同的波传播方法导致通过处理地震数据生成的不同图像和/或用于生成通过处理地震数据生成的不同图像的不同时间量。因此，可能存在特定选择的介质(或波场传播子)可能受益于特定类型的波传播方法的情况。本申请详述至少一种波传播方法，所述波传播方法可以有益地与倾斜正交晶(TOR)介质结合使用。

通过介绍的方式，可以使用各种地震勘测系统和技术获取地震数据，相对于图2和图3论述所述地震勘测系统和技术中的两个。不管利用的地震数据收集技术，在获取地震数据之后，计算系统可以分析所获取的地震数据并且可以使用地震数据分析的结果(例如，地震图、地质构造图等)来在烃勘探和生产行业中执行各种操作。例如，图1说明详述可以基于对所获取的地震数据的分析而进行的各种过程的方法10的流程图。尽管以特定次序描述方法10，但是应注意，方法10可以以任何合适次序执行。

现在参考图1，在框12处，可以基于分析的地震数据来确定与相应地震勘测相关联的地球的地下区域内的烃沉积物的位置和特性。在一个实施例中，可以分析获取的地震数据，以生成说明地下区域内的各种地质构造的地图或轮廓。基于烃沉积物的识别的位置和特性，在框14处，可以勘探地下区域的某些位置或部分。也就是说，烃勘探组织可以使用烃沉积物的位置来确定地下区域表面处钻入地球的位置。因此，烃勘探组织可以使用烃沉积物的位置和特性以及相关联的覆盖层来确定沿着其钻入地球的路径、如何钻入地球等。

在勘探装备已放入地下区域内之后，在框16处，可以经由自然自流井、人工升举井等产生存储于烃沉积物中的烃。在框18处，可以经由运输车、管线等将产生的烃运输到炼油厂等。在框20处，可以根据各种精炼程序来处理产生的烃以使用烃开发不同产物。

应注意，关于方法10论述的过程可以包括其它合适的过程，这些过程可以基于经由一个或多个地震勘测获取的地震数据中所指示的烃沉积物的位置和特性。因此，应理解，上述过程并不意图描绘可以在确定地下区域内的烃沉积物的位置和特性之后执行的过程的详尽列表。

考虑到前述内容，图2是海洋勘测系统22(例如，用于结合图1的框12使用)的示意图，所述海洋勘测系统可以用于获取关于海洋环境中的地球的地下区域的地震数据(例如，波形)。通常，可以在海洋24中或在位于海底28下方的地球的地下区域26上方的其它水体中进行使用海洋勘测系统22的海洋地震勘测。

海洋勘测系统22可以包括船30、一个或多个地震源32、(地震)拖缆34、一个或多个(地震)接收器36，和/或可以帮助获取表示地球的地下区域26内的地质构造的地震图像的其它装备。船30可以拖曳地震源32(例如，气枪阵列)，所述地震源可以产生指向海底28的能量，诸如声波(例如，地震波形)。船30还可以拖曳具有接收器36(例如，水听器)的拖缆34，所述接收器36可以获取地震波形，所述地震波形表示在从地下区域26内的各种地质构造(例如盐丘、断层、褶皱等)反射之后由地震源32输出的能量。另外，尽管海洋勘测系统22的描述用一个地震源32(图2中表示为气枪阵列)和一个接收器(在图2中表示为一组水听器)进行描述，但是应注意，海洋勘测系统22可以包括多个地震源32和多个接收器36。通过相同方式，尽管海洋勘测系统22的以上描述用一个拖缆34进行描述，但是应注意，海洋勘测系统22可以包括类似于拖缆34的多个拖缆。另外，额外的船30可以包括额外的源32、拖缆34等以执行海洋勘测系统22的操作。

图3是陆地勘测系统38(例如，用于结合图1的框12使用)框图，所述陆地勘测系统可以用于获得关于非海洋环境中的地球的地下区域26的信息。陆地勘测系统38可以包括陆基地震源40和陆基接收器44。在一些实施例中，陆地勘测系统38可以包括多个陆基地震源40以及一个或多个陆基接收器44和46。实际上，出于论述目的，陆地勘测系统38包括陆基地震源40和两个陆基接收器44和46。陆基地震源40(例如，地震振动器)可以安置于在所关注地下区域26上方的地球的表面42上。陆基地震源40可以产生指向地球的地下区域26的能量(例如，声波、地震波形)。在到达地下区域26内的各种地质构造(例如盐丘、断层、褶皱)之后，由陆基地震源40输出的能量可以从地质构造中反射出来，并且由一个或多个陆基接收器(例如，44和46)获取或记录。

在一些实施例中，陆基接收器44和46可以分散在地球的表面42上以形成网格状图案。因此，响应于能量经由地震源40指向地下区域26，每个陆基接收器44或46可以接收反射的地震波形。在一些情况下，由地震源40产生的一个地震波形可以从不同地质构造反射回来并且由不同接收器接收。例如，如图3中所示，地震源40可以示出可以指向地下区域26的能量作为地震波形48。第一接收器44可以接收地震波形48从一个地质构造的反射，并且第二接收器46可以接收地震波形48从不同地质构造的反射。因此，第一接收器44可以接收反射的地震波形50，并且第二接收器46可以接收反射的地震波形52。

不论如何获取地震数据，计算系统(例如，用于结合图1的框12使用)可以分析由接收器36、44、46获取的地震波形，以确定关于地下区域26内的地质结构、烃沉积物的位置和特性等的地震信息。图4是此计算系统60的示例的框图，所述计算系统可以执行各种数据分析操作以分析由接收器36、44、46获取的地震数据，以确定结构和/或预测地下区域26内的地质构造的地震特性。

现在参考图4，计算系统60可以包括通信组件62、处理器64、存储器66、存储装置68、输入/输出(I/O)端口70，以及显示器72。在一些实施例中，计算系统60可以省略显示器72、通信组件62和/或输入/输出(I/O)端口70中的一个或多个。通信组件62可以是可以促进接收器36、44、46、一个或多个数据库74、其它计算装置，和/或其它具有通信能力的装置之间的通信的无线或有线通信组件。在一个实施例中，计算系统60可以经由网络组件、数据库74等接收接收器数据76(例如，地震数据、地震图等)。计算系统60的处理器64可以分析或处理接收器数据76，以确定与地球的地下区域26内的地质构造有关的各种特征。

处理器64可以是能够执行计算机可执行代码的任何类型的计算机处理器或微处理器。处理器64还可以包括可以执行下文描述的操作的多个处理器。存储器66和存储装置68可以是可以用作介质以存储处理器可执行代码、数据等的任何合适的制品。这些制品可以表示计算机可读介质(例如，任何适当形式的存储器或存储装置)，所述计算机可读介质可以存储由处理器64用于执行当前公开的技术的处理器可执行代码。通常，处理器64可以执行软件应用程序，所述软件应用程序包括根据本文描述的实施例处理经由地震勘测的接收器获取的地震数据的程序。

存储器66和存储装置68也可以用于存储数据、数据分析、软件应用程序等。存储器66和存储装置68可以表示可以存储由处理器64使用以执行本文中描述的各种技术的处理器可执行代码的非暂时性计算机可读介质(例如，任何合适形式的存储器或存储装置)。应注意，非暂时性仅指示介质是有形的，而不是信号。

I/O端口70可以是可以耦合到诸如输入装置(例如，键盘、鼠标)、传感器、输入/输出(I/O)模块等的其它外围组件的接口。I/O端口70可以使计算系统60能够经由I/O端口70与海洋勘测系统22、陆地勘测系统38等中的其它装置通信。

显示器72可以描绘与正由处理器64处理的软件或可执行代码相关联的可视化。在一个实施例中，显示器72可以是能够从计算系统60的用户接收输入的触摸显示器。显示器72还可以用于查看和分析所获取的地震数据的分析结果，以确定地下区域26内的地质构造、地下区域26内的烃沉积物的位置和特性、与地下区域26内的一个或多个井相关联的地震特性的预测等。显示器72可以是任何合适类型的显示器，诸如，液晶显示器(LCD)、等离子体显示器，或有机发光二极管(OLED)显示器等。除了经由显示器72描绘本文所描述的可视化之外，还应注意，计算系统60还可以经由例如纸(例如，经由打印)等的其它有形元件来描绘可视化。

考虑到前述内容，本文中描述的本技术也可以使用超级计算机来执行，所述超级计算机采用多个计算系统60、云计算系统等在多个计算系统60上分布要执行的过程。在这种情况下，作为超级计算机的一部分运行的每个计算系统60可能不包括作为计算系统60的一部分列出的每个组件。例如，每个计算系统60可能不包括显示器72，因为多个显示器72可能不可用于被设计为连续地处理地震数据的超级计算机。

在执行各种类型的地震数据处理之后，计算系统60可以将分析结果存储在一个或多个数据库74中。数据库74可以经由通信组件62通信地耦合到可以向计算系统60传输数据以及从计算系统60接收数据的网络。另外，数据库74可以存储关于地下区域26的信息，诸如关于地下区域26的先前地震图、地质样本数据、地震图像等。

尽管已经关于计算系统60讨论了上述组件，但是应注意，类似的组件可以组成计算系统60。此外，计算系统60也可以是海洋勘测系统22或陆地勘测系统38的一部分，且因此可以监视和控制源32或40、接收器36、44、46等的某些操作。此外，应注意，所列出的组件被提供为示例组件，并且本文描述的实施例不限于参考图4描述的组件。

在一些实施例中，计算系统60可以基于经由上述接收器接收的地震数据来生成地下区域26的二维表示或三维表示。另外，可以组合与多个源/接收器组合相关联的地震数据，以创建可以延伸一定距离的地下区域26的近乎连续的轮廓。在二维(2-D)地震勘测中，接收器位置可以沿着单条线放置，而在三维(3-D)勘测中，接收器位置可以以网格图案分布在整个表面上。因此，当地球层存在于记录位置的正下方时，2-D地震勘测可以提供地球层的横截面图片(垂直切片)。另一方面，3-D地震勘测可以创建可以对应于地下区域26的3-D图片的数据“立方体”或体积。

另外，4-D(或延时)地震勘测可以包括在3-D勘测期间多次获取的地震数据。使用在不同时间获取的不同地震图像，计算系统60可以比较两个图像以识别地下区域26中的变化。

在任何情况下，地震勘测可以由非常大量的个别地震记录或地震道构成。因此，计算系统60可以被采用以分析所获取的地震数据，以获得表示地下区域26的图像并确定烃沉积物的位置和特性。为此，可以使用多种地震数据处理算法来从所获取的地震数据中去除噪声，偏移预处理的地震数据，识别多个地震图像之间的偏移，对准多个地震图像等。

在计算系统60分析所获取的地震数据之后，地震数据分析的结果(例如，地震图、地震图像、地质构造图等)可以用于在烃勘探和生产行业中执行各种操作。例如，如上所述，所获取的地震数据可以用于执行图1的方法10，图1详述可以基于对所获取的地震数据的分析而进行的各种过程。

在一些实施例中，可以结合例如图5中说明的方法78的地震处理方案，生成地震数据分析的结果。如所说明，在地震数据的任何解释、与期望的勘探目标相一致的任何进一步图像增强、从处理后的地震数据中的属性生成、地震数据根据需要的重新解释，以及钻探前景或其它地震勘测应用的确定和/或生成之前，方法78包括地震处理序列，所述地震处理序列包括步骤80中的地震数据收集、步骤82中的地震数据编辑、步骤84中的初始处理，以及步骤86中的信号处理、调节和成像(例如，其可以包括成像截面或体积的产生，以及例如结合所说明的偏移过程使用的在本文中描述的波传播方法)。由于处理方法78，可以识别烃在地下区域26内的位置。用于地震建模的技术通常且更具体地是波传播方法，所述波传播方法可以有益地结合倾斜正交晶(TOR)介质使用以改进地震数据的处理结果，以便允许更好地检测地下区域26内的烃。

在一个实施例中，生成用于TOR介质的双向标量波场传播子。TOR介质的波场传播子可以在对应介质中采用伪声P波分散关系。用于TOR介质的波场传播子组合快速傅立叶变换(FFT)和有限差分(FD)运算，从而产生倾斜正交晶FFD(TORFFD)波传播方法，所述TORFFD波传播方法作为标量算子还可以克服TOR介质中存在的qP波和qSV波的耦合。

当用于TOR介质时，类似于用于TOR介质的伪谱(PS)波传播方法，TORFFD波传播方法(以及合成的波场传播子)可以例如利用3D空间中的每波场采样的两个点。然而，通过利用TORFFD波传播方法代替PS波传播方法会提高效率，例如，与用于PS波传播方法的30个FFT相比，利用TORFFD波传播方法包括在每个时间步长中计算2个正向FFT、3个逆FFT和1个FD运算。此外，相对于PS波传播方法，使用TORFFD波传播方法允许更大的时间步长，例如，库朗数高达0.555以及处理时间加速(例如，高达PS波传播方法的大约两倍、大约三倍、大约四倍、大约五倍或大约六倍的速度)，同时实现了相对于PS波传播方法的总精度提高。

下文描述本TORFFD波传播方法和相关联的波传播子的理论基础。用于地震建模和逆时偏移(例如，用于成像地震数据的技术)的声波方程式表示为：

如在方程式1中所阐述，p(x,t)是地震压力波场并且v(x)是传播速度。方程式1的声压解可以表示为：

在恒速介质的情况下，借助于FFT，方程式2可以作为解运算。在具有可变速度的模型中，方程式2可以通过用v(x)代替v来提供近似值。此外，在TOR的情况下，方程式2右侧上的项v(x)|k|可以用声学近似值f(v,k)代替，从而得出：

如在方程式3中所阐述，θ是相对于垂直测量到的倾角，是TI对称轴在水平面中的投影与原始X坐标之间的方位角，ψ是围绕TI轴线的旋转角，ε₁、ε₂、δ₁、δ₂和δ₃是正交晶介质的扩展Thomsen参数，v_z是沿着TI对称轴的P波相速度，方向 而“快”速度沿着方向 并且/>和/>表示在与新轴线对准的旋转坐标系中评估的波数。

当使用方程式2在非均匀介质中传播波时，FFT可能不再直接用于从波数域回到空间域的傅立叶逆变换。作为解决方案，可以利用傅立叶有限差分(FFD)作为FFT和FD算子的链算子，以便在每个时间步长处，FFD方法首先使用全局参考参数传播当前波场，并且然后根据局部模型参数的变化通过有限差分应用波场校正。用于TOR介质的波场传播子的本技术可以通过将方程式2右手侧上的空间波数混合域项变换成两个项的乘积来利用类似的方法：

如方程式4中所阐述，_v0表示参考参数。方程式4右手侧上的第一项仅仅是k的函数。因此，我们可以在计算过程中使用逆FFT返回到空间域。可以通过将波数域中的相移变换成空间有限差分来如下近似求解剩余项G(v,k)：

如方程式5中所阐述，i,j＝1,2,3；i≠j；i,j≠n且a,b_n,c_pn,c_mn和d_n从围绕k＝0的泰勒展开式确定，这些如表1中说明为空间中的局部模型参数的变量并且用作4阶FD系数(其中i＝1,2,3)：

表1

表2提供表1中出现的变量的其它表达式，其中w_0i具有与w_i类似的表达式，并且所有空间变量v由对应的常数模型参数v₀替代。

表2

扩展到具有TOR介质的情况，为了从p(x,t–Δt)和p(x,t)获得p(x,t+Δt)，可以进行以下步骤。可以进行经由FFT的p(x,t)到(k,t)的变换。接下来，可以/>(k,t)乘以以获得/>(k,t)，然后可以通过使用逆FFT将/>(k,t)变换成q(k,t)。可以利用先前在表1中指出的系数将有限差分应用于q(k,t)，以生成u(x,t)。可以实现经由FFT的u(x,t)到/>(x,t)的变换并且可以应用波数域中的校正项u₀，使得并且/>随后可以经由使用逆FFT执行分别/>(k,t)到u₀(k,t)以及/>(k,t)到u₁(k,t)的变换。随后，可以例如根据u(x,t)的局部速度变化组合u₀(x,t)和u₁(x,t)以产生p(x,t+Δt)←u(x,t+Δt)+2p(x,t)-p(x,t-Δt)，由此q(x,t)、/>(k,t)、u(x,t)、/>(k,t)、u₀(x,t)、/>(k,t)、u₁(x,t)和/>(k,t)是临时函数。另外，v₁表示具有最小速度的位置处的模型参数，G(v₁,k)是方程式5的精确响应，H(v₁,k)是FD的近似值。此外，/>在密度相对较低的位置处用作波数域中的波场的校正项，而/>用作第二(例如，另一)校正项。

前述过程可以由计算系统60作为TORFFD波传播方法来进行。在一些实施例中，前述步骤中的至少一些可以以替代顺序执行或完全省略。此外，尽管步骤可以被理解为由计算系统60的处理器64执行，但是应理解，任何合适的装置或系统，或者包括处理器64的合适装置或系统的组合可以执行上述步骤，诸如在计算系统60外部但通信地耦合到计算系统60的处理单元或计算装置或系统的电路系统，并且步骤的执行可以涉及处理器64与存储在有形机器可读介质上的软件结合或利用所述软件进行操作，以在生成TORFFD波传播方法(以及相关联的波场传播子)时执行上述步骤。

另外，应理解，由处理器64执行的上述TORFFD TORFFD波传播方法表示对常规系统的特定改进，以及具有改进功能的改进计算系统60。具体地，当与通过例如计算机系统执行常规计算机算法来执行地震数据分析和/或在不使用TORFFD波传播方法和/或其相关联的TORFFD波场传播子的情况下执行地震数据分析相比较时，通过计算系统60利用TORFFD波传播方法和/或使TORFFD波场传播子作为计算系统的一部分例如相对于替选的波传播模型以及其相关联的波形传播子提高准确性和效率。

作为表示一些前述优点的视觉辅助，图6说明TOR模型中的波传播子响应的示例。在图6中说明通过PS波传播方法的恒定TOR模型中的伪声波场快照88。同样，另外说明恒定TOR模型中的FFD方法(即，TORFFD波传播方法)伪声波场快照90。如可以观察到，与伪声波场快照88中存在的分散相比，伪声波场快照90的分散减小。另外，伪声波场快照90的生成过程比伪声波场快照88的生成过程快至多六倍。另外，PS波传播方法和PS波传播方法两者的使用对于伪声波传播具有优势，因为不存在qP波和qSv波的耦合。

作为TORFFD波传播方法的优点的第二组示例，图7说明利用TTIFDD波传播方法生成的逆时偏移(RTM)图像92，图8说明利用PS波传播方法生成的RTM图像94，并且图5说明使用上述TORFFD波传播方法生成的RTM图像96。生成RTM图像92和RTM图像96的运行时间为大约20分钟，而生成RTM图像94的运行时间为大约两个小时。另外，可以理解，人们可以观察到RTM图像92中的一些主要事件；然而，RTM图像92中的断层平面似乎没有明确定义。相反，RTM图像96中的断层边界比RTM图像92中的断层边界更清晰，这使得呈现强方位各向异性的能力得到改进。另外，相对于RTM图像92，RTM图像96中的下方或相邻断层事件更清晰地成像。RTM图像94同样似乎相对于RTM图像92提供改进的成像；然而，应注意，生成RTM图像94的运行时间是生成RTM图像96的运行时间的大致6倍。因此，利用TORFFD波传播方法似乎可以增加与TTI FFD波传播方法相关联的处理时间，同时(至少)相对于所利用的PS波传播方法保持或逼近生成的RTM图像中的图像质量。

已经通过示例的方式示出上述特定实施例，并且应理解，这些实施例可能易于进行各种修改和替选形式。应进一步理解，权利要求书并不旨在限于所公开的特定形式，而是覆盖落入本公开的精神和范围内的所有修改、等效物和替选方案。

本文提出并要求保护的技术被引用并应用于具有实际性质的物质对象和具体示例，这些物质对象和具体示例可证明地改进本技术领域，因此不是抽象的、无形的或纯粹理论上的。此外，如果附在本说明书末尾的任何权利要求包含表示为“用于[执行][功能]…的装置”或“用于[执行][功能]...的步骤”的一个或多个要素，则预期这些元素应根据35U.S.C.112(f)进行解释。然而，对于包含以任何其它方式表示的元素的任何权利要求，预期不应根据35 U.S.C.112(f)解释此类元素。

Claims (19)

1.一种用于地震数据处理的系统，包括：

处理器，其中，所述处理器被配置成：

接收烃储层的储层数据；

接收与波场传播子的选择相关的指示；

利用傅立叶有限变换和有限差分应用所述波场传播子，以至少部分地通过执行与所述波场传播子相关联的压力属性到变换后压力属性的初始变换，对与倾斜正交晶介质相关联的波场建模，所述倾斜正交晶介质表示包括所述烃储层的地下区域；以及

结合所述波场传播子处理所述储层数据，以生成用于在包括所述烃储层并且含有结构或地层特征的所述地下区域上方的地震勘探的输出，所述结构或地层特征有助于烃的呈现、偏移或累积。

2.根据权利要求1所述的系统，其中，所述处理器被配置成将快速傅立叶变换应用于所述压力属性，以执行与所述波场传播子相关联的所述压力属性的所述初始变换。

3.根据权利要求1所述的系统，其中，所述处理器被配置成至少部分地通过产生所述变换后压力属性与传播速度值的积以生成第一值，对与所述倾斜正交晶介质相关联的所述波场建模，所述倾斜正交晶介质表示包括所述烃储层的所述地下区域。

4.根据权利要求3所述的系统，其中，所述处理器被配置成至少部分地通过执行所述第一值到变换后第一值的变换，对与所述倾斜正交晶介质相关联的所述波场建模，所述倾斜正交晶介质表示包括所述烃储层的所述地下区域。

5.根据权利要求4所述的系统，其中，所述处理器被配置成将逆快速傅立叶变换应用于所述第一值，以执行所述第一值的所述变换。

6.根据权利要求4所述的系统，其中，所述处理器被配置成至少部分地通过利用预定系数将有限差分运算应用于所述变换后第一值以生成第二值，对与所述倾斜正交晶介质相关联的所述波场建模，所述倾斜正交晶介质表示包括所述烃储层的所述地下区域。

7.根据权利要求6所述的系统，其中，所述处理器被配置成至少部分地通过执行所述第二值到变换后第二值的变换，对与所述倾斜正交晶介质相关联的所述波场建模，所述倾斜正交晶介质表示包括所述烃储层的所述地下区域。

8.根据权利要求7所述的系统，其中，所述处理器被配置成将快速傅立叶变换应用于所述第二值，以执行所述第二值的所述变换。

9.根据权利要求7所述的系统，其中，所述处理器被配置成至少部分地通过将波数域中的校正项应用于所述变换后第二值以生成修改后第二值，对与所述倾斜正交晶介质相关联的所述波场建模，所述倾斜正交晶介质表示包括所述烃储层的所述地下区域。

10.根据权利要求9所述的系统，其中，所述处理器被配置成至少部分地通过执行所述修改后第二值到第三值的变换，对与所述倾斜正交晶介质相关联的所述波场建模，所述倾斜正交晶介质表示包括所述烃储层的所述地下区域。

11.根据权利要求10所述的系统，其中，所述处理器被配置成至少部分地通过将所述波数域中的校正项应用于所述第二值以生成第二修改后第二值，对与所述倾斜正交晶介质相关联的所述波场建模，所述倾斜正交晶介质表示包括所述烃储层的所述地下区域。

12.根据权利要求11所述的系统，其中，所述处理器被配置成至少部分地通过执行所述第二修改后第二值到第四值的变换，对与所述倾斜正交晶介质相关联的所述波场建模，所述倾斜正交晶介质表示包括所述烃储层的所述地下区域。

13.根据权利要求12所述的系统，其中，所述处理器被配置成至少部分地通过执行所述第三值和所述第四值分别到变换后第三值和变换后第四值的变换，对与所述倾斜正交晶介质相关联的所述波场建模，所述倾斜正交晶介质表示包括所述烃储层的所述地下区域。

14.根据权利要求13所述的系统，其中，所述处理器被配置成将快速傅立叶变换应用于所述第三值和所述第四值，以执行所述第三值和所述第四值的所述变换。

15.根据权利要求13所述的系统，其中，所述处理器被配置成至少部分地通过将所述变换后第三值与所述变换后第四值组合，对与所述倾斜正交晶介质相关联的所述波场建模，所述倾斜正交晶介质表示包括所述烃储层的所述地下区域。

16.一种由计算系统的处理器实施的用于地震数据处理的方法，包括：

接收烃储层的储层数据；

利用傅立叶有限变换和有限差分应用波场传播子，以通过执行与所述波场传播子相关联的压力属性到变换后压力属性的初始变换，对与倾斜正交晶介质相关联的波场建模，所述倾斜正交晶介质表示包括烃储层的地下区域；以及

结合所述波场传播子处理储层数据，以生成用于在包括烃储层并且含有结构或地层特征的所述地下区域上方的地震勘探的输出，所述结构或所述地层特征有助于烃的呈现、偏移或累积。

17.根据权利要求16所述的方法，包括至少部分地通过以下操作对与所述倾斜正交晶介质相关联的所述波场建模，所述倾斜正交晶介质表示包括所述烃储层的所述地下区域：

产生所述变换后压力属性与传播速度值的积以生成第一值；

执行所述第一值到变换后第一值的变换；

利用预定系数将有限差分运算应用于所述变换后第一值以生成第二值；

执行所述第二值到变换后第二值的变换；

将波数域中的校正项应用于所述变换后第二值以生成修改后第二值；

执行所述修改后第二值到第三值的变换；

将所述波数域中的校正项应用于所述第二值以生成第二修改后第二值；

执行所述第二修改后第二值到第四值的变换；

执行所述第三值和所述第四值分别到变换后第三值和变换后第四值的变换；以及

将所述变换后第三值与所述变换后第四值组合。

18.一种或多种有形的非暂时性机器可读介质，包括被配置成使处理器处理地震数据的指令，其中，处理所述地震数据包括：

利用傅立叶有限变换和有限差分应用波场传播子，以通过执行与所述波场传播子相关联的压力属性到变换后压力属性的初始变换，对与倾斜正交晶介质相关联的波场建模，所述倾斜正交晶介质表示包括烃储层的地下区域；以及

结合所述波场传播子处理储层数据，以生成用于在包括烃储层并且含有结构或地层特征的所述地下区域上方的地震勘探的输出，所述结构或地层特征有助于烃的呈现、偏移或累积。

19.根据权利要求18所述的一种或多种有形的非暂时性机器可读介质，其中，处理所述地震数据还包括至少部分地通过以下操作对与所述倾斜正交晶介质相关联的所述波场建模，所述倾斜正交晶介质表示包括所述烃储层的所述地下区域：

产生所述变换后压力属性与传播速度值的积以生成第一值；

执行所述第一值到变换后第一值的变换；

利用预定系数将有限差分运算应用于所述变换后第一值以生成第二值；

执行所述第二值到变换后第二值的变换；

将波数域中的校正项应用于所述变换后第二值以生成修改后第二值；

执行所述修改后第二值到第三值的变换；

将所述波数域中的校正项应用于所述第二值以生成第二修改后第二值；

执行所述第二修改后第二值到第四值的变换；

执行所述第三值和所述第四值分别到变换后第三值和变换后第四值的变换；以及

将所述变换后第三值与所述变换后第四值组合。