专利名称：具有吸收边界和随机边界的逆时迁移的制作方法
技术领域：
本发明涉及波场在与地下区域相关联的地球模型中的逆时迁移(Reverse timemigration)，其中，该地球模型可以包括被限定为具有吸收特征和随机扰动特征两者的至少ー个边界区域。
背景技术：
逆时迁移被广泛接受为ー种用于在地质容积(特别是具有复杂结构的那些)中勘探与生产的优选成像技木。不管其计算成本，逆时迁移例如因最优化数值解算机和使用计算机集群而变得日益可承担得起。近年来，诸如GPU和FPGA的协处理器的可用性已经提供了对于显著改进逆时迁移效率的额外展望。然而，在常规方法中，在执行逆时迁移的同时使用保存的源波场和/或按每ー个时间步长的边界值在完全利用由协处理器提供的计算功率的优点方面仍是ー种限制因素。

发明内容
本发明的ー个方面涉及ー种用于生成与地下区域有关的图像的计算机实现方法。所述方法可以包括以下步骤源波场从初始时间状态至最终时间状态在时间中从初始波场状态起向前传播通过与地下区域相关联的地球模型。对于传播源波场来说，该地球模型可以包括被限定为具有吸收特征和随机扰动特征两者的至少ー个边界区域，以使源波场响应于源波场撞击在所述至少一个边界区域上而被阻尼和随机散射。所述方法可以包括通过使源波场从最终时间状态至初始时间状态在时间中向后传播通过该地球模型来确定向后传播的源波场，以重建初始波场状态。对于确定向后传播的源波场来说，所述至少一个边界区域可以被重新限定为具有助推特征和相同随机扰动特征两者，以使该源波场在所述至少一个边界区域处被放大。所述方法可以包括使接收器波场从最终时间状态起传播通过该地球模型。对于传播接收器波场来说，所述至少一个边界区域可以被重新限定为具有吸收特征，以使接收器波场响应于该接收器波场撞击在所述至少一个边界区域上而被阻尼。所述方法可以包括至少部分地基于向后传播的源波场和接收器波场来生成与该地下区域有关的图像。本发明的另一方面涉及ー种用于生成与地下区域有关的图像的系统。该系统可以包括被配置成执行计算机程序模块的一个或多个处理器。该计算机程序模块可以包括向前传播模块，该向前传播模块被配置成使源波场从初始时间状态至最终时间状态，在时间中从初始波场状态起向前传播通过与地下区域相关联的地球模型。该地球模型可以包括被限定为具有吸收特征和随机扰动特征两者的至少ー个边界区域，以使源波场响应于源波场撞击在所述至少一个边界区域上而被阻尼和随机散射。该计算机程序模块可以包括向后传播模块，该向后传播模块被配置成，通过使源波场从最终时间状态至初始时间状态在时间中向后传播通过该地球模型来确定向后传播的源波场，以重建该初始波场状态。对于确定向后传播的源波场来说，所述至少一个边界区域可以被重新限定为具有助推特征和相同随机扰动特征两者，以使该源波场在所述至少一个边界区域处被放大。该计算机程序模块可以包括接收器波场传播模块，该接收器波场传播模块被配置成使接收器波场从最终时间状态起传播通过该地球模型。对于传播接收器波场来说，所述至少一个边界区域可以被重新限定为具有吸收特征，以使接收器波场响应于接收器波场撞击在所述至少一个边界区域上而被阻尼。该计算机程序模块可以包括成像模块，该成像模块被配置成至少部分地基于向后传播的源波场和接收器波场来生成与该地下区域有关的图像。本发明的又一方面涉及ー种计算机可读存储介质，该计算机可读存储介质上具体实施有指令。该指令可以由处理器执行以进行用于生成与地下区域有关的图像的方法。所述方法可以包括使源波场从初始时间状态至最终时间状态在时间中从初始波场状态起向前传播通过与地下区域相关联的地球模型。对于传播源波场来说，该地球模型可以包括被限定为具有吸收特征和随机扰动特征两者的至少ー个边界区域，以使源波场响应于源波场撞击在所述至少一个边界区域上而被阻尼和随机散射。所述方法可以包括通过使源波场从最终时间状态至初始时间状态在时间中向后传播通过该地球模型来确定向后传播的源波场，以重建该初始波场状态。对于确定向后传播的源波场来说，所述至少一个边界区域可以被重新限定为具有助推特征和相同随机扰动特征两者，以使该源波场在所述至少ー个边 界区域处被放大。所述方法可以包括使接收器波场从最终时间状态起传播通过该地球模型。对于传播接收器波场来说，所述至少一个边界区域可以被重新限定为具有吸收特征，以使接收器波场响应于接收器波场撞击在所述至少一个边界区域上而被阻尼。所述方法可以包括至少部分地基于向后传播的源波场和接收器波场来生成与该地下区域有关的图像。当參照附图考虑下面的描述和所附权利要求书时，本发明的这些和其它目的、特征和特性、以及结构的相关元件与部件的组合的操作方法与功能以及制造的经济性将变得更清楚，其全部形成了本说明书的一部分，其中，相同标号指定各个图中的对应部分。然而，应当明白，附图仅仅是出于例示和描述的目的，而非旨在作为对本发明的限制的限定。如在本说明书和权利要求书中使用的，単数形式“一(a)”、“一(an)”，以及“该/所述(the)”包括多个指示物，除非上下文另外清楚地規定。


图I例示了根据本发明的一个或多个实施例的、被配置成生成与地下区域有关的图像的系统。图2例示了在其中限定了吸收边界的地球模型中向前传播期间的示例性波场状态。图3例示了在其中限定随机边界和吸收边界的地球模型中向前传播期间的示例性波场状态。图4例示了在其中未限定助推边界的地球模型中的示例性重建波场。图5例示了在其中限定助推边界的地球模型中的示例性重建波场。图6例示了根据本发明的一个或多个实施例的、用于生成与地下区域有关的图像的方法。
具体实施方式
图I例示了根据本发明的一个或多个实施例的、被配置成生成与地下区域有关的图像的系统100。示例性实施例提供了用于大量缩减和/或去除存储器存取的新方式。这通过利用地球模型中的吸收、助推和/或随机边界条件保留并且稍后在计算域中恢复源波场来实现。与常规方法相比，这些边界条件可以被用于缩减人工反射波场中的伪像，同时按用于成像的振幅恢复有用信号。此外，可以使用时间中的检查点来刷新波场状态并且去除累积的数值误差。这样，本发明的实施例可以实现改进的信噪比和改进的图像质量。在一个实施例中，系统100包括电子存储部102、用户接ロ 104、ー个或多个信息资源106、ー个或多个处理器108和/或其它组件。在一个实施例中，电子存储部102包括电子地存储信息的电子存储介质。电子存储部102的电子存储介质可以包括与系统100集成地(即，基本上不可去除地)设置的系统存储部，和/或例如经由端ロ(例如,USB端ロ、火线(firewire)端ロ等)或驱动器(例如，盘驱动器等)可去除地连接至系统100的可去除存储部。电子存储部102可以包括ー个或多个光学可读存储介质(例如，光盘等)、磁可读存储介质(例如，磁带、磁硬盘驱动器、软盘驱动器等)、基于电荷的存储介质(例如、EEPROM、RAM等)、固态存储介质(例如，闪存驱动器等)， 和/或其它电子可读存储介质。电子存储部102可以存储软件算法、由处理器108确定的信息、经由用户接ロ 104接收的信息、从信息资源106接收的信息、和/或使得系统100能够如在此所述地起作用的其它信息。电子存储部102可以是系统100内的分离组件、或者电子存储部102可以与系统100的ー个或多个其它组件(例如，处理器108)集成地设置。用户接ロ 104被配置成在系统100与用户之间设置用户可以通过其向系统100提供信息和从系统100接收信息的接ロ。这使得统称为“信息”的数据、结果和/或指令和任何其它可传送项能够在用户与系统100之间传送。如在此使用的，术语“用户”可以指单一个体或者可以协作地工作的个体组。适于包含在用户接ロ 104中的接口装置的示例包括下列中的ー个或多个小键盘、按钮、开关、键盘、旋钮(knob)、杠杆、显示屏、触摸屏、扬声器、麦克风、指示灯、声音报警器和/或打印机。在一个实施例中，用户接ロ 104实际上包括多个分离接ロ。应当明白，其它通信技术(硬布线或无线)也可以被本发明设想为用户接ロ 104。例如，本发明设想用户接ロ 104可以与由电子存储部102提供的可去除存储接ロ集成。在这个实施例中，可以将信息从可去除存储部(例如，智能卡、闪存驱动器、可去除盘等)加载到系统100中，使得用户能够定制系统100的实现。适于与系统100 —起使用作为用户接ロ104的其它示例性输入装置和技术包括但不限于，RS-232端ロ、RF链路、IR链路、调制解调器(电话、线缆或其它)。简单地说，用干与系统100传送信息的任何技术被本发明设想为用户接ロ 104。信息资源106包括与下列有关的ー个或多个信息源关心的地下地质容积、生成关心的地下地质容积的地球模型的处理、关心的地下地质容积的实际地球模型，和/或与关心的地下地质容积有关的其它信息。通过非限制例的方式，信息资源106之一可以包括在关心的地下容积处或附近获取的地震数据、从其导出的信息、和/或有关该获取的信息。这种地震数据可以包括源波场和接收器波场。该地震数据可以包括地震数据的各个迹线(trace )。地震数据的迹线包括记录在ー个检测器处的、表示从源起传播通过关心的地下容积的地震能量的数据。根据地震数据导出的信息例如可以包括速度模型、与用于对地震能量传播通过关心的地下容积建模的射束相关联的射束特性、与用于对地震能量传播通过关心的地下容积建模的射束相关联的Green函数、和/或其它信息。与地震数据的获取有关的信息例如可以包括与以下各项有关的数据地震能量源的位置和/或取向、地震能量的一个或多个检测器的位置和/或取向、能量被源生成并且被引导至关心的地下容积的时间，和/或其它信息。关心的地下容积的地球模型可以表示该关心的容积的、诸如空间结构和成分特性的各种特征。地球模型可以描述具有不同成分的区域之间的接触面和/或显示不同地震能量传播速度的区域之间的接触面。处理器108被配置成在系统100中提供信息处理能力。这样，处理器108可以包括数字处理器、模拟处理器、被设计成处理信息的数字电路、被设计成处理信息的模拟电路、状态机、和/或用于电子地处理信息的其它机构中的ー个或多个。尽管处理器108在图I中被示出为单ー实体，但这仅仅是出于例示性目的。在某些实现中，处理器108可以包括多个处理単元。这些处理单元可以物理地位于同一装置或计算平台内，或者处理器108可以表示协作操作的多个装置的处理功能。如图I中所示，处理器108可以被配置成执行ー个或多个计算机程序模块。该ー 个或多个计算机程序模块可以包括向前传播模块110、向后传播模块112、接收器波场传播模块114、成像模块116和/或其它模块中的ー个或多个。处理器108可以被配置成通过软件；硬件；固件；软件、硬件和/或固件的某ー组合；和/或用于配置处理器108上的处理能力的其它机制来执行模块110、112、114和/或116。应当清楚，尽管模块110、112、114和116在图I中被例示为共同位于单ー处理单元内，但在其中处理器108包括多个处理单元的实现中，模块110、112、114和/或116中的一个或多个可以相对于其它模块远程定位。由于不同的模块110、112、114和/或116中的任ー个都可以提供比所描述的更多或更少的功能，因而下面描述的对由不同的模块110、112、114和/或116提供的功能的描述出于例示性目的，而不是旨在进行限制。例如，模块110、112、114和/或116中的一个或多个可以被消除，并且其ー些或全部功能可以由模块110、112、114和/或116中的其它模块提供。作为另ー示例，处理器108可以被配置成执行可以执行下面归因于模块110、112、114和/或116中的ー个的ー些或全部功能的一个或多个附加模块。向前传播模块110可以被配置成，使源波场在时间上从初始波场状态起向前传播通过与地下区域相关联的地球模型。这包括从在呈现初始波场状态时的初始时间状态至最终时间状态传播源波场。与地球模型相关联的地下区域是有限区域。即，该地下区域可以包括用于三维实现的有限容积、用于ニ维实现的有限区、或者用于其它实现的某其它有限区域。这样，该地球模型被边界围绕。一般而言，当地球模型的边界仅表示为陡峭边缘吋，该边界将充任用于在地球模型中传播的波场的反射器。来自地球模型的边界的反射被视为伪像，即，它们未准确地描绘地下区域中的实际波场行为。为了缩减在地球模型中的传播波场的人为反射的存在，该地球模型可以包括相对于该地球模型的其余部分具有不同条件或特征的影响波场传播特性的一个或多个边界区域。边界区域可以包括沿该地球模型的边界的单ー线，或者沿该地球模型的边界的区。该地球模型可以包括被限定为具有吸收特征和/或随机扰动特征的至少ー个边界区域。这样，源波场可以响应于该源波场撞击在所述至少一个边界区域上而在所述至少一个边界区域处被阻尼和/或随机散射。
根据示例性实施例，吸收特征可以包括衰减功能。一般而言，衰减功能包括将使波场的振幅在指定空间范围(例如，边界区域的跨度)上减小的任何功能。地球模型的边界区域处的这种衰减功能有效地使波场在该波场撞击该边界区域时在该边界区域处被部分或完全吸收。边界区域处的吸收特征允许缩减或消除被该边界区域反射的波场。最終，缩减或消除这种反射将缩减地下区域的图像中的伪像。根据ー些实施例，在向前传播源波场期间，边界区域的吸收特征可以被配置为缩减而非消除边界反射的振幅，因为这些事件以后可以被用于重建初始波场状态。然而，在向后传播接收器波场期间，可以应用边界区域处的更强衰减以便消除(或尽可能多地缩减)接收器波场的边界反射。根据各个实施例，随机扰动特征可以包括归因于边界区域的各个区的随机速度，和/或使波场在撞击边界区域时在该边界区域处被散射的其它特性。这样，即使存在衰减边界区域，在地球模型中幸存的反射也将失真。这种失真的反射将无助于相干的、可识别的图像事件。在示例性实施例中，向前传播模块110还可以被配置成在从初始时间状态向最终时间状态传播源波场的同时保存或者以其它方式保留ー个或多个波场状态。所保留的波场 状态可以由电子存储部102保存并且可从其检索。在向前传播波场期间保存或者以其它方式保留波场状态可允许在向后传播该波场期间缩减和/或消除数值误差，如结合向后传播模块112进ー步讨论的。向后传播模块112可以被配置成确定向后传播的源波场。在示例性实施例中，向后传播的源波场可以通过使源波场从最终时间状态至初始时间状态在时间中向后传播通过该地球模型来确定，以重建该初始波场状态。对于向后传播来说，就如结合向前传播模块110描述的向前传播来说，地球模型的所述至少一个边界区域可以被限定为具有助推特征(而非衰减特征)和相同随机扰动特征。这样，源波场在向后传播期间在所述至少一个边界区域处被放大。在示例性实施例中，该助推特征等于在向前传播期间使用的吸收特征的逆。在一些实施例中，向后传播模块112还可以被配置成在确定向后传播的源波场的同时，利用在向前传播期间保留的波场状态。这种保留的波场状态可以由向后传播模块112从电子存储部102检索。所保留的波场状态可以用于缩减向后传播的源波场中的数值误差(例如，舍入误差)。接收器波场传播模块114可以被配置成使接收器波场从最终时间状态起传播通过地球模型。对于传播接收器波场来说，该地球模型的所述至少一个边界区域可以被限定为具有吸收特征，以使接收器波场响应于该接收器波场撞击在所述至少一个边界区域上而在所述至少一个边界区域处被阻尼。在一些实施例中，接收器波场传播模块114可以被配置成与通过向后传播模块112确定向后传播的源波场同时传播接收器波场。成像模块116可以被配置成至少部分地基于向后传播的源波场和接收器波场来生成与该地下区域有关的图像。根据示例性实施例，成像模块116可以被配置成将由向后传播模块112确定的向后传播的源波场与接收器波场交叉关联以生成图像。基于向后传播的源波场和接收器波场的其它图像生成技术也可被设想并且处于本发明的范围内。图2和3对比在向前波场传播的各个阶段，归因于地球模型中的边界的吸收特征、非吸收特征、随机扰动特征、和非随机扰动特征的示例性效果。图2例如例示了在在其中限定吸收边界的地球模型200中向前传播(在这个示例中，向前传播通常对应于向下方向)期间的示例性波场状态。更具体地说，在图2中，地球模型200的左侧边界202被限定为具有吸收特征，而右侧边界204没有。这样，相对于来自右侧边界204的反射208，来自左侧边界202的反射206具有减小的振幅。图3例示了在其中限定了随机边界和吸收边界的地球模型300中向前传播(在这个示例中，向前传播通常对应于向下方向)期间的示例性波场状态。更具体地说，地球模型300的左侧边界302、下侧边界304、和右侧边界306被限定为具有随机扰动特征。左侧边界302还被限定为具有吸收特征。地球模型300中的反射(如来自左侧边界302的反射308和来自右侧边界306的反射310)因随机扰动特征而随机地失真。然而，相对于因归因于左侧边界302的吸收特征而造成的反射310，反射308的振幅被显著削弱。与常规方法相比，与衰减组合的这种失真可以提供改进的图像质量。图4和5对比在向后传播期间归因于地球模型中的边界的助推特征和非助推特征的示例性效果。更具体地说，图4例示了在其中未限定助推边界的地球模型400中的示例性重建波场。在波场在地球模型400中向前传播(在这个示例中，向后传播通常对应于向上方向)期间，吸收特征仅归因于左侧边界402，然而，非助推特征在向后传播期间归因于左侧 边界402。结果，在地球模型400中可以看到所重建波场的多个部分中的伪像404和振幅损失 406。图5例示了在其中限定助推边界的地球模型500中的示例性重建波场。在波场在地球模型500中向前传播期间，吸收特征归因于左侧边界502，而助推特征在向后传播(在这个示例中，向后传播通常对应于向上方向)期间归因于左侧边界502。相对于图4中例示的重建波场，图5的重建波场缺乏明显的伪像和振幅损失。图6例示了根据本发明的一个或多个实施例的、用于生成与地下区域有关的图像的方法600。下面呈现的方法600的操作g在例示。在一些实施例中，方法600可以利用未描述的ー个或多个附加操作和/或不利用所讨论操作中的ー个或多个来完成。另外，其中图6例示和下面描述的方法600的操作的次序不g在进行限制。在一些实施例中，方法600可以在ー个或多个处理装置(例如，数字处理器、模拟处理器、被设计成处理信息的数字电路、被设计成处理信息的模拟电路、状态机、和/或用于电子地处理信息的其它机构)中实现。该ー个或多个处理装置可以包括响应于电子地存储在电子存储介质上的指令来执行方法600的ー些或全部操作的ー个或多个装置。该ー个或多个处理装置可以包括通过要具体设计用于执行方法600的ー个或多个操作的硬件、固件、和/或软件配置的ー个或多个装置。在操作602，源波场从初始时间状态至最终时间状态，在时间中从初始波场状态起向前传播通过与地下区域相关联的地球模型。该地球模型可以包括被限定为具有吸收特征和随机扰动特征两者的至少ー个边界区域，以使源波场响应于该源波场撞击在所述至少一个边界区域上而被阻尼井随机散射。在示例性实施例中，操作602可以通过向前传播模块110来执行。在操作604，向后传播的源波场，通过使源波场从最终时间状态至初始时间状态在时间中向后传播通过该地球模型来确定，以重建该初始波场状态。所述至少一个边界区域可以被重新限定为具有助推特征和如在操作602中的相同随机扰动特征两者，以使该源波场在所述至少一个边界区域处被放大。在一些实施例中，可以执行向后传播模块112，以执行操作604。在操作606，接收器波场从最终时间状态起传播通过地球模型。所述至少一个边界区域可以被重新限定为具有吸收特征，以使接收器波场响应于该接收器波场撞击在所述至少ー个边界区域上而在所述至少一个边界区域处被阻尼。根据ー些实施例，可以执行接收器波场传播模块114以执行操作606。在操作608，至少部分地基于向后传播的源波场和接收器波场来生成与该地下区域有关的图像。在某些实现中，交叉关联向后传播的源波场和接收器波场以产生图像。根据示例性实施例，可以执行成像模块116以执行操作608。尽管基于当前被认为是最实用且优选的实施例，出于例示的目的对本发明进行了详细描述，但应当明白，这种细节仅用于该目的，并且本发明不限于所公开的实施例，而且正相反，其g在覆盖处于所附权利要求书的精神和范围内的修改例和等同布置。例如，应当 明白，本发明在尽可能的情况下设想可以将任何实施例的一个或多个特征与任何其它实施例的ー个或多个特征相组合。
权利要求
1.一种用于生成与地下区域有关的图像的计算机实现的方法，该方法包括 使源波场从初始时间状态至最终时间状态在时间中从初始波场状态起向前传播通过与地下区域相关联的地球模型，其中，该地球模型包括被限定为具有吸收特征和随机扰动特征两者的至少一个边界区域，以使源波场响应于该源波场撞击在所述至少一个边界区域上而被阻尼和随机散射； 通过使源波场从最终时间状态至初始时间状态在时间中向后传播通过该地球模型来确定向后传播的源波场，以重建该初始波场状态，其中，所述至少一个边界区域被重新限定为具有助推特征和所述随机扰动特征两者，以使该源波场在所述至少一个边界区域处被放大； 使接收器波场从最终时间状态起传播通过该地球模型，其中，所述至少一个边界区域被重新限定为具有吸收特征，以使接收器波场响应于该接收器波场撞击在所述至少一个边界区域上而被阻尼；以及 至少部分地基于向后传播的源波场和接收器波场来生成与该地下区域有关的图像。
2.根据权利要求I所述的方法，其中，使接收器波场传播的步骤与确定向后传播的源 场的步骤是同时执行的。
3.根据权利要求I所述的方法，所述方法还包括在使源波场从初始时间状态向最终时间状态传播的同时，保留一个或多个波场状态。
4.根据权利要求3所述的方法，所述方法还包括在确定向后传播的源波场的同时利用所保留的一个或多个波场状态来缩减数值误差。
5.根据权利要求I所述的方法，其中，该吸收特征包括衰减功能。
6.根据权利要求I所述的方法，其中，助推特征等于吸收特征的逆。
7.根据权利要求I所述的方法，其中，随机扰动特征包括随机速度。
8.根据权利要求I所述的方法，其中，生成图像的步骤包括交叉关联向后传播的源波场和接收器波场。
9.一种用于生成与地下区域有关的图像的系统，该系统包括 一个或多个处理器，被配置成执行计算机程序模块，该计算机程序模块包括 向前传播模块，被配置成使源波场从初始时间状态至最终时间状态在时间中从初始波场状态起向前传播通过与地下区域相关联的地球模型，其中，该地球模型包括被限定为具有吸收特征和随机扰动特征两者的至少一个边界区域，以使源波场响应于该源波场撞击在所述至少一个边界区域上而被阻尼和随机散射； 向后传播模块，被配置成通过使源波场从最终时间状态至初始时间状态在时间中向后传播通过该地球模型来确定向后传播的源波场，以重建该初始波场状态，其中，所述至少一个边界区域被重新限定为具有助推特征和所述随机扰动特征两者，以使该源波场在所述至少一个边界区域处被放大； 接收器波场传播模块，被配置成使接收器波场从最终时间状态起传播通过该地球模型，其中，所述至少一个边界区域被重新限定为具有吸收特征，以使接收器波场响应于该接收器波场撞击在所述至少一个边界区域上而被阻尼；以及 成像模块，被配置成至少部分地基于向后传播的源波场和接收器波场来生成与该地下区域有关的图像。
10.根据权利要求9所述的系统，其中，该接收器波场传播模块还被配置成，与向后传播模块确定向后传播的源波场同时地使接收器波场传播。
11.根据权利要求9所述的系统，其中，该向前传播模块还被配置成，在使源波场从初始时间状态向最终时间状态传播的同时，保留一个或多个波场状态。
12.根据权利要求11所述的系统，其中，该向后传播模块还被配置成，在确定向后传播的源波场的同时利用所保留的一个或多个波场状态来缩减数值误差。
13.根据权利要求9所述的系统，其中，该吸收特征包括衰减功能。
14.根据权利要求9所述的系统，其中，助推特征等于吸收特征的逆。
15.根据权利要求9所述的系统，其中，随机扰动特征包括随机速度。
16.根据权利要求9所述的系统，其中，该成像模块还被配置成在生成图像的同时，交叉关联向后传播的源波场和接收器波场。
17.一种其上具体实施有指令的计算机可读存储介质，该指令可由处理器执行以进行用于生成与地下区域有关的图像的方法，所述方法包括以下步骤 使源波场从初始时间状态至最终时间状态在时间中从初始波场状态起向前传播通过与地下区域相关联的地球模型，其中，该地球模型包括被限定为具有吸收特征和随机扰动特征两者的至少一个边界区域，以使源波场响应于该源波场撞击在所述至少一个边界区域上而被阻尼和随机散射； 通过使源波场从最终时间状态至初始时间状态在时间中向后传播通过该地球模型来确定向后传播的源波场，以重建该初始波场状态，其中，所述至少一个边界区域被重新限定为具有助推特征和所述随机扰动特征两者，以使该源波场在所述至少一个边界区域处被放大； 使接收器波场从最终时间状态起传播通过该地球模型，其中，所述至少一个边界区域被重新限定为具有吸收特征，以使接收器波场响应于该接收器波场撞击在所述至少一个边界区域上而被阻尼；以及 至少部分地基于向后传播的源波场和接收器波场来生成与该地下区域有关的图像。
全文摘要
与地下区域有关的图像可以至少部分地基于向后传播的源波场和接收器波场来生成。源波场可以从初始时间状态至最终时间状态在时间中从初始波场状态起向前传播通过与地下区域相关联的地球模型。向后传播的源波场可以从最终时间状态至初始时间状态在时间中向后传播源波场来确定，以重建该初始波场状态。接收器波场可以从最终时间状态起传播通过地球模型。该地球模型可以包括可以被限定为具有吸收特征、助推特征、随机扰动特征、和/或其它特征中的一个或多个特征的至少一个边界区域。这样，波场可以在所述至少一个边界区域处被阻尼、放大、随机散射、和/或以其它方式改变。这些波场可以被用于构建具有改进信噪比的地下区域图像。
文档编号G01V1/28GK102792186SQ201180013718
公开日2012年11月21日 申请日期2011年4月19日 优先权日2010年6月24日
发明者刘伟 申请人:雪佛龙美国公司