专利名称：使用快速面向目标照明计算的地震成像系统及方法
技术领域：
地震学用于需要地质信息的勘探、考古研究和工程项目。勘探地震学所提供的数据在与其他可利用的地球物理学、钻井和地质数据联合使用时，能够提供关于岩石类型的结构和分布及其含量的信息。这种信息对于搜寻水、地热储层以及诸如碳氢化合物和矿石的矿床具有极大的帮助。大多数石油公司依靠勘探地震学来选择钻探石油探井的地点。
背景技术：
传统的地震学使用人工产生的地震波来勘查地下结构。地震波从源处向下传入地球并从地下结构之间的分界处反射。地面检波器(receiver)探测并记录反射的地震波以供后续分析。尽管直接检验记录的信号可以察觉到一些大尺度结构，但必须对记录的信号进行处理，以消除畸变并且展现地下图像中细微的细节。现有的各种处理方法不能充分消除畸变，并且需要过长的计算时间。本文公开了改进的系统和方法。
结合下列附图考虑接下来的详细描述，可以更好地理解各个公开的实施例，其中图1显示了例证性的海洋地震勘测环境的侧视图；图2显示了例证性的海洋地震勘测环境的顶视图；图3显示了指定通道接收到的交替激发(flip-flop)炮点所产生的例证性的中心点模式(midpoint pattern)；图4显示了例证性的地震勘测记录系统；图5显示了例证性的踪迹集；图6显示了例证性的三维数据体；图7显示了例证性的炮点几何位置；图8显示了例证性的地震成像方法的流程图；以及图9显示了例证性的成像系统。尽管本发明适用于多种修改及替代形式，但在附图中以示例方式显示了其特定实施例，并在本文中对这些特定实施例进行了详细描述。然而，应理解，附图及详细描述不是要将所公开的实施例限制于所显示的特殊形式，相反，意图是覆盖落入到所附权利要求的实质和范围内的所有修改、等同或替代形式。
具体实施例方式本文公开了使用快速面向目标照明计算技术的多种地震成像系统和方法。计算近似照明值的数据体或“矩阵”，并用其来估计真实反射率值的图像矩阵。从格林函数(Green’ s function)获得照明值，格林函数不是以逐个炮点为基础计算或重新计算，而是以多炮点组进行计算，并与滚动和(rolling-sum)相结合，从而大大降低计算量(computational overhead)。因此，与那些依赖于常规三维波动方程照明的系统相比，本发明所公开的系统和方法能够更加快速地和/或以更好的质量提供目标区域照明值。在至少一些实施例中，本发明所公开的地震成像系统包括至少一个存储装置、 存储器和至少一个处理器。存储装置存储来自于指定勘测区域的地震勘测的炮点道集 (gather)。存储器存储地震成像软件，该地震成像软件在被处理器运行时配置系统，以从存储装置获得多炮点道集并将那些具有重叠的检波器位置的炮点道集分组以形成多炮点道集。对于多炮点道集中的每个炮点和检波器位置，系统同时确定勘测区域中的潜在散射点处的格林函数波场。(如本文进一步详述的，同时获得这些格林函数能够避免大量的重新计算)。系统将格林函数组合，以发现每个潜在散射点处的照明值，并且至少部分基于照明值来计算相应的地震反射率值。系统根据这些反射率值来产生勘测区域的照明，其能够接受进一步的处理或者进行显示以供分析者使用。在一些实施方式中，每个多炮点道集对应于海洋地震勘测中的一条航海线。通过频域内的小波束分解(beamlet decomposition)和直接局部平面波分解，使用单程波动方程传播来确定格林函数，从而可以得到其他的效果。每次对勘测区域的一个深度切片执行每个格林函数的同时确定，并在格林函数传播到下一深度切片之前执行照明值确定。如上文所述，如果在滚动和的基础上计算所需要的求和项(summation term)的其中之一，则可以显著节省照明值确定的计算量。通过将照明值的计算限制在感兴趣的目标区域，而不是对整个勘测体积执行计算，可以实现进一步的节省。下面描述所公开的成像系统和方法的上述和其他变形。通常在陆地上和海上都执行勘探地震。在海上，地震勘测船如图1所示在船的后面配置拖缆。当船向前(沿箭头102的方向)移动时每个拖缆110跟在船100的后面， 每个拖缆包括多个均勻间隔的检波器114。每个拖缆110可以进一步包括可编程换向器 (diverter) 118和可编程深度控制器，其将拖缆从船只航线向外拖到工作偏移距(见图2) 并向下拖到所需要的操作深度(图1)。拖缆110可以长达数千米，通常被构造很多段，每一段的长度为25至100米，包括多达35个检波器或者更加均勻间隔的检波器的组。每个拖缆110包括用于使检波器114 和船100上的地震设备互相连接的电缆布线或光缆布线。数据在检波器114附近数字化， 并通过电缆以每秒钟7 (或更高)百万比特数据的速率传输到船100。如图1所示，地震勘测船100还可以拖带一个或多个源112。源112可以是脉冲源或振动源。海洋地震学中所使用的检波器114通常被称为水下地震检波器(hydrophone)， 通常用压电式换能器构造。可以使用各种合适类型的水下地震检波器，例如盘式水下地震检波器和圆柱形水下地震检波器。源112和检波器114通常配置在海面104以下。船上的处理设备控制源和检波器的操作，并且记录所采集的数据。地震勘测提供用于在海面104以下成像以展现地下结构(诸如位于海底以下108 的结构106)的数据。分析者可以绘制出地下层的形貌并且研究所记录的地震数据的特征， 对数据的研究用来确定石油和/或天然气储层的位置。为了对地下结构106成像，源112发射地震波116，地震波在由于地下结构106(和其他地下结构)而存在声阻抗变化的地方发生反射。通过检波器114的模式来探测反射波。 通过记录(作为时间的函数)已从源112传播到地下结构106再到检波器114的到达地震波116，在进行适当的数据处理之后可以获得地下结构106的图像。图2显示了拖曳一组拖缆110和两个源112的地震勘测船100的俯视图(未按比例)。当船100向前移动时，以所谓的交替激发模式(flip-flop pattern)交替地触发源。 使用可编程换向器来提供拖缆之间大致均勻的间隔。拖缆上的指定位置处的检波器与共同的场文件踪�：�(common field file trace number)或共同的通道(common channel) 202 相关联。图3显示了两个炮点的例证性的源和检波器的俯视图。对于第一炮点而言，在位置302触发一个源，检波器阵列的所示部分处于位置304(用虚线轮廓线显示)。对于第二炮点而言，在位置306触发源，检波器阵列的所示部分处于位置308(用实线轮廓线显示)。 假设目前反映的地下结构是水平的，则到达十二个检波器中的每一个检波器的地震波从源和检波器位置之间的中心点下方的位置反射。因此，第一炮点从十二个中心点308的下方产生反射(用具有垂直交叉阴影线的虚线轮廓线显示)，而第二炮点从十二个中心点310的下方产生反射(用具有水平交叉阴影线的实线轮廓线显示)。作为一个示例，向量312显示地震能量从炮点302到中心点314的传播，相等长度的向量316显示传播到检波器位置的反射地震能量。对于第二炮点306而言，向量318和向量320显示相似的传播路径。应注意，中心点314是被多个炮点“击中”的中心点的其中之一，从而当来自于炮点的信息被处理和组合时，可以获得来自于这些区域的更多的信号能量。地震勘测(用于陆地和海洋两者)通常被设计为提供均勻分布的中心点网格，从而每个中心点均有相当高的平均命中次数。图4显示了例证性的地震勘测记录系统，其具有连接到总线402以将数字信号传送到勘测船100上的数据记录电路406的检波器114。位置信息及其他参数传感器404也连接到数据记录电路406，以使数据记录电路能够存储对解释所记录的数据有益的附加信息。作为例证，这种附加信息可以包括阵列方向信息和速度信息。通用数字数据处理系统408被显示为连接到数据记录电路406，进一步被显示为通过总线402连接至定位装置410和地震源112。处理系统408配置记录电路406、定位装置410和地震源112的操作。记录电路406将高速数据流从检波器114采集到诸如光盘或磁盘存储阵列的非易失性存储介质上。定位装置410(包括可编程换向器和深度控制器) 控制检波器114和源112的位置。图4的地震记录系统可包括未在此处具体显示的另外的元件。例如，每个拖揽110 可以具有用于连接到数据记录电路的独立总线402。处理系统408包括具有图形显示和键盘或其他接受用户输入的方法的用户接口，处理系统408可以进一步包括网络接口，该网络接口用于将存储的地震勘测数据传送到具有用于处理地震勘测数据的强大计算资源的中央计算设施。图5显示了检波器114探测和采样的例证性的地震信号(也被称为“踪迹”)。信号指示随时间变化的地震波能量的某种度量(例如位移、速度、加速度、压强)，以可编程采样率在高分辨率(例如对位)下对其进行数字化。可以以不同的方式对这种信号分组，当这样分组时，其被称为“道集”。例如，“共中心点道集”是在定义区域内具有中心点的一组踪迹。“炮点道集”是对地震源的单次引爆(firing)记录的一组踪迹。“多炮点道集”是一组炮点道集，在海洋地震勘测中通常包括沿航海线记录的全部踪迹。
当保持成图5的格式时，所记录的地震勘测数据用处不大。尽管可以用展现大尺度地下结构的曲线图并排标绘制各种记录波形，但是这种结构是变形的，甚至无法看见细微的结构。因此，对数据进行处理，以产生数据体，即数据值的三维阵列，例如如图6所示。 数据体表示整个勘测区域内的某种地震属性。三维阵列包括大小均勻的单元，每个单元具有表示该单元的地震属性的数据值。可以表示各种地震属性，在一些实施例中，每个单元具有多个数据值，以表示多个地震属性。合适的地震属性的示例包括反射率、声阻抗、声速和密度。体数据格式本身更容易进行计算分析和视觉渲染，因此，可将数据体称为勘测区域的 “三维图像”。图7显示了各个参数是如何与例证性的炮点的几何位置相关联的。地震能量从地震源沿射线702传播到目标分界面704，并沿射线706向检波器反射。在反射点(在别处用 (X，y，ζ)坐标表示，缩写为向量无)，表面具有法向量云。在具有水平面内的χ轴和y轴以及朝上的ζ轴的预定坐标系统中，法向量巧与每个轴形成角度。法向量巧与χ轴之间的角度为 0)(，法向量乃与7轴之间的角度为θ γ，法向量污与ζ轴之间的角度为θη。入射射线702和反射射线706相对于法向量成相等(但相反)的“张”角θρ入射射线702相对于ζ轴形成入射角θ s，反射射线相对于ζ轴形成接收(又叫做散射)角0g。根据下面的关系式，可将入射角和散射角对(θ s · θ g)转换为法向角和张角对(θ n · 0r)民=(乾+《)/2总=(这-之)/2。(1)根据炮点的几何位置的这种理解，我们现在开始讨论照明分析及其在地下区域的成像中的作用。宽泛地说，照明分析包括确定有多少地震能量能够被指定体积在检波器的方向上反射。同样地，照明是反射点上的每次击中的炮点几何位置的函数。照明分析是用于研究图像质量上的采集孔径(acquisition aperture)和重叠结构(overlaying structure)的强大工具。大多数用于预测照明度分布的现有技术是基于射线跟踪模拟。例如参见Bear，G，LU, C.，Lu，R.和Willen，D. ,2000, "The construction of subsurface illumination and amplitude maps via ray tracing", The Leading Edge，19(7)，726—728 ；Muerdter, D.禾口 Ratcliff，D. ，2001，"Understanding subsalt illumination through ray-trace modeling,Part 1 :Simple 2D salt models", The Leading Edge,20(6), 578-594；Muerdter, D. ，Kelly, Μ.禾口 Ratcliff，D. ，2001， "Understanding subsalt illumination through ray-trace modeling, Part 2 :Dipping salt bodies, salt peaks, and nonreciprocity of subsalt amplitude response", The Leading Edge，20(7)，688—687 ；Muerdter, D.禾口 Ratcliff，D. ，2001，"Understanding subsalt illumination through ray-trace modeling,Part3 :Salt ridges and furrows, and impact of acquisition orientation，，，The Leading Edge, 20 (8) , 803-816 ; Schneider, W. A.禾口 Winbow, G. A. , 1999, "Efficient and accurate modeling of 3-D seismic illumination，，，SEG Expanded abstracts 18,633-636。尽管射线跟踪成本并不高，但是由于所涉及的高频近似以及射线跟踪所产生的奇异性问题，所得到的照明图(illumination map)可能在复杂区域中具有大的误差。参见例如 Hoffmann,2001,“Illumination,resolution and image quality of PP-and PS-waves for survey planning”，The Leading Edge，20 (9)，1008-1014。为了对付这些问题，使用小波束分解和直接局部平面波分解来得到频域波场的局部角度域信息。更具体而言，使用mi，R. S.禾口 Chen, L. 2002, "Mapping directional illumination and acquisition-aperture efficacy by beamlet propagator，，，SEG Expanded Abstracts 21,1352 ；以及Xie, X. B.禾口 ffu,R. S. ,2002,"Extracting angle related image from migrated wavefield",Expanded abstracts, SEG 72ndAnnual Meeting，1360-1363所公开的技术来确定作为位置、频率以及入射或散射角的函数的源和检波器波场。这些分解使得能够使用基于波动方程的方法 (特别包括单程波动方程传播子)来执行照明分析。关于这种方法的更多细节可在Xie， X. B, Jin, S禾口 Wu, R. S. , 2003, "Three-dimensional illumination analysis using wave equation based propagator，，，SEG Expanded Abstracts 22，989 中获得。这种方法的现有实施方式已被证明计算量太大并且是缓慢的，原因在于其逐个炮点来计算照明的方法。本发明公开的用于计算照明的方法可以使用下列三个特征中的一个或多个特征来加速照明计算计算目标区域而不是整个数据体的照明；多炮点多炮点地来执行计算， 而不是逐个炮点执行计算；在照明计算中使用滚动和，而不是每次都从头开始重新求和。最终目标是生成地震反射率的数据体，之后作为用来定位储层和矿床的工作的一部分，分析者可以对其进行研究。如Luo，M和Xie，X. B.，2005，“Amplitude recovery from back propagated waves to true scattered waves",Technical Report No. 12,Modeling and Imaging Pro ject, University of California, Santa Cruz, 25-34所述,“真实反身寸率，，
图像可写成
权利要求
1.一种地震成像方法，包括 获得用于勘测区域的多个炮点道集；将具有重叠的检波器位置的炮点道集分组，以形成至少一个多炮点道集； 对于指定的多炮点道集中的每个炮点位置和每个检波器位置，同时确定所述勘测区域中的潜在散射点处的格林函数；组合所述格林函数，以发现每个所述潜在散射点处的照明值； 至少部分基于所述照明值来发现每个所述潜在散射点处的反射率值；以及向用户显示所述勘测区域的图像，其中所述图像是至少部分基于所述反射率值的。
2.根据权利要求1所述的方法，进一步包括对其他多炮点道集重复所述确定操作、组合操作和发现操作；以及累积从多个多炮点道集确定的反射率值，以形成所述勘测区域的所述图像。
3.根据权利要求2所述的方法，其中在海洋地震勘测中，每个多炮点道集表示一个航海线。
4.根据权利要求1所述的方法，其中通过频域内的小波束分解和直接局部平面波分解，用单程波动方程来确定所述格林函数。
5.根据权利要求4所述的方法，其中每次对所述勘测区域的一个深度切片执行每个格林函数的所述同时确定，用于每个深度切片的确定与所述组合相互交错。
6.根据权利要求4所述的方法，其中在执行所述组合之前，在整个目标区域上执行每个格林函数的所述同时确定。
7.根据权利要求1所述的方法，其中指定散射点f处的所述反射率值能够被表示为频率ω、炮点索引s、入射角93和散射角θ g处的反射率分量的总和，其中每个所述反射率分量能够被表示为卷积图像值A^ 為氏, )与照明值D(^ ，《，之比，所述卷积图像值能够被表示为向前传播的源波场与向后传播的接收波场之间的频域乘积。
8.根据权利要求7所述的方法，其中照明值伐元馋《，&)能够被表示为源贡献与检波器贡献之间2(无 ，^ )的乘积的总和
9.根据权利要求8所述的方法，其中所述组合包括通过在滚动和的基础上确定检波器贡献水无馋冬 )来利用不同炮点的重叠的检波器位置。
10.根据权利要求1所述的方法，其中仅对作为所述勘测区域的子集的目标区域中的潜在散射点执行所述组合。
11.一种地震成像系统，包括至少一个存储装置，所述存储装置存储来自于指定勘测区域的地震勘测的炮点道集；存储器，所述存储器存储地震成像软件；以及至少一个处理器，所述处理器连接到所述存储器，以执行所述地震成像软件，其中所述软件将所述至少一个处理器配置为从所述存储装置检索多个炮点道集；将具有重叠的检波器位置的炮点道集分组，以形成至少一个多炮点道集；对于每个多炮点道集中的每个炮点位置和每个检波器位置，同时确定所述勘测区域中的潜在散射点处的格林函数；组合所述格林函数，以发现每个所述潜在散射点处的照明值；至少部分基于所述照明值来发现每个所述潜在散射点处的反射率值；以及生成用于向用户显示的所述勘测区域的图像，其中所述图像是至少部分基于所述反射率值的。
12.根据权利要求11所述的系统，其中所述软件进一步配置所述至少一个处理器，以便累积从多个多炮点道集确定的反射率值，以形成所述勘测区域的所述图像。
13.根据权利要求12所述的系统，其中在海洋地震勘测中，每个多炮点道集表示一个航海线。
14.根据权利要求11所述的系统，其中通过频域内的小波束分解和直接局部平面波分解，用单程波动方程来确定所述格林函数。
15.根据权利要求14所述的系统，其中每次对所述勘测区域的一个深度切片执行每个格林函数的所述同时确定，用于每个深度切片的确定与所述组合相互交错。
16.根据权利要求14所述的系统，其中在执行所述组合之前，在整个目标区域上执行每个格林函数的所述同时确定。
17.根据权利要求11所述的系统，其中指定散射点处的所述反射率值能够被表示为频率ω、炮点索引s、入射角93和散射角θ 8处的反射率分量的总和，其中每个所述反射率分量能够被表示为卷积图像值元心,5, , )与照明值之比，所述卷积图像值氏, )能够被表示为向前传播的源波场与向后传播的接收波场之间的频域乘积。
18.根据权利要求17所述的系统，其中照明值伐瓦坏込，^)能够被表示为源贡献与检波器贡献之间成兄 ,5,的乘积的总和
19.根据权利要求18所述的系统，其中作为组合的一部分，所述格林函数发现照明值， 所述软件配置所述至少一个处理器以通过在滚动和的基础上确定检波器贡献d(元 ，&<9g) 来利用不同炮点的重叠的检波器位置。
20.根据权利要求11所述的系统，其中所述软件配置所述至少一个处理器，以便仅对作为所述勘测区域的子集的目标区域中的潜在散射点组合所述格林函数以发现照明值。
全文摘要
在各种公开的地震成像系统和方法中，使用快速面向目标照明计算技术来获得近似照明值的数据体或“矩阵”。通过这些照明值可以发现“真实反射率”值的图像矩阵。从格林函数获得照明值，格林函数不是以逐个炮点为基础计算或重新计算的，而是以多炮点组进行计算，并与滚动和相结合，从而大大降低计算量。因此，与那些依赖于常规三维波动方程照明的系统相比，所公开的系统和方法能够更加快速地和/或以更好的质量提供目标区域照明。
文档编号G01V1/00GK102395902SQ200980158750
公开日2012年3月28日 申请日期2009年4月16日 优先权日2009年4月16日
发明者罗明秋, 金胜汶 申请人:兰德马克图形公司
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
