CN1886678A - 向量3分量3维克希霍夫预堆栈迁移 - Google Patents
向量3分量3维克希霍夫预堆栈迁移 Download PDFInfo
- Publication number
- CN1886678A CN1886678A CNA2004800354655A CN200480035465A CN1886678A CN 1886678 A CN1886678 A CN 1886678A CN A2004800354655 A CNA2004800354655 A CN A2004800354655A CN 200480035465 A CN200480035465 A CN 200480035465A CN 1886678 A CN1886678 A CN 1886678A
- Authority
- CN
- China
- Prior art keywords
- data
- components
- receiver location
- receiver
- picture point
- 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
- 230000005012 migration Effects 0.000 title claims abstract description 58
- 238000013508 migration Methods 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 claims abstract description 76
- 230000008569 process Effects 0.000 claims description 21
- 230000015572 biosynthetic process Effects 0.000 claims description 16
- 238000012937 correction Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 11
- 241001269238 Data Species 0.000 claims description 7
- 238000010008 shearing Methods 0.000 claims description 7
- 230000006835 compression Effects 0.000 claims description 5
- 238000007906 compression Methods 0.000 claims description 5
- 238000005553 drilling Methods 0.000 claims description 4
- 238000013507 mapping Methods 0.000 abstract description 3
- 238000003384 imaging method Methods 0.000 description 12
- 230000007704 transition Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 5
- 150000003839 salts Chemical class 0.000 description 5
- 238000012546 transfer Methods 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 3
- 238000002310 reflectometry Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000007418 data mining Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000003760 hair shine Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000003969 polarography Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000035945 sensitivity Effects 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/284—Application of the shear wave component and/or several components of the seismic signal
-
- 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
-
- 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
- G01V1/305—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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/42—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators in one well and receivers elsewhere or vice versa
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Analysis (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
一种用于迁移通过地下接收机和地面地震源获取的三分量、三维(3-C,3-D)地震数据的设备和方法。这个方法利用完全3分量反射波场。它使用动态的向量能量映射方法来成像反射位置和把每个时间样本仅仅映射到它的反射的图像点。所以,这个方法减少不想要的数据涂污和虚假镜像图像。这个方法克服只使用单分量轨迹或在1C或3-C3-D VSP迁移中的预先旋转的三分量轨迹的弱点和产生更好的3-D图像。
Description
技术领域
本发明涉及一种改进地震迁移的精度的、地球物理勘探的方法。具体地,本发明使用偏移或零偏移探查测量来精确地迁移在三维(3-D)表面地震数据中、在垂直地震剖面(VSP)中和在交叉井(cross-well)中的地震探查数据中呈现的反射器。
背景技术
在表面地震探查中,由地震源注入地球的能量从地下地理特性反映出来,以及被多个接收机记录。这个过程通过使用形成线(2D获取)或覆盖面积(3-D获取)的、源和接收机结构被重复许多次。得到的数据被处理来通过使用被称为迁移的过程产生反射器的图像。
传统的反射地震学利用表面源和接收机来检测来自地下阻抗对比度(impedance contrast)的反射。得到的图像由于在源、反射器和接收机之间的长的和复杂的行进路径,常常在空间精度、分辨率和相干性方面受到损害。具体地,由于地震信号通过具有低的、横向变化的速度的高吸收性靠近表面的天气层的双向通道,地下的图像可能是质量很差的。为了克服这个困难,通常被称为垂直地震剖面(VSP)的技术被开发来在钻孔附近的地下成像。在VSP中,使用表面地震源,以及信号在井下接收机或井下接收机阵列处被接收。这是对于接收机(或接收机阵列)的不同的深度重复进行的。在偏移VSP中,多个间隔开的源被顺序地激活,使得能够成像比起用单个源可能的更大的距离范围。
VSP数据获取可以通过在部分或全部完成钻井后把接收机沿缆线输送到井下而完成。VSP方法的优点在于,数据质量比起表面数据获取好得多。VSP获取也可以通过把接收机阵列输送到井下作为井底组件(BHA)的一部分而完成。这被称为VSP同时钻井。
授权给Wyatt等的美国专利4,627,036给出VSP方法的早先的例子,该专利全部在此引用以供参考。现在参照图1,图上显示典型的、用于地面地震获取的VSP结构。在示例的图上,Vibroseis源11被显示为把能量传递给地球。应当指出,如果想要的话,可以利用任何其它适当的地震源,诸如爆炸。在水下环境中,源可以是气枪或水下震动器。
接收机12被显示为位于在钻孔14中想要的深度。对于接收机12的位置,能量在点16处从地下岩层15被反射。从接收机12产生的输出被记录车17记录。在VSP中,接收机12典型地被移动到新的位置用于每次拍摄,在地音探测器位置之间的距离是某个恒定的距离,诸如50英尺。如果想要的话,可以利用间隔开某个想要的距离的接收机阵列或可以使用多个间隔开的源。
由VSP得到的数据具有图2所示的样子。Wyatt讨论被称为VSP-CDP的处理技术的使用,通过这种技术,诸如图2所示的那样的VSP数据可被堆叠来产生远离井的地下的图像。
地面地震数据的迁移的处理过程被使用来得到比起通过CDP或堆叠方法可得到的那些更好的地下的图像。在迁移时,目的是放置地震反射在它们的适当的空间位置:在表面CDP方法中,另一方面,假设反射发源于在源与接收机之间中途的反射点。通常使用的、用于迁移的方法是克希霍夫方法,其中定义用于地下的速度模型。计算从源到不同的点和从不同的点到接收机的行进时间。反射器的实际的图像是通过组合从多个源接收机对到多个成像点的数据而得到的。如果速度模型是合理地精确的,则信号在正确的图像点处建设性地干涉。这个概念原先是对于地面地震数据开发的。Wiggins(1984)把迁移的使用扩展到其中观察表面不限于是平的水平面的情形。克希霍夫迁移对于VSP数据的使用已由Dillon讨论。
Dillon的教导限于2-D迁移。最近,VSP克希霍夫深度迁移由Bicquart被使用于3-D VSP数据。正如Bicquart指出的,克希霍夫和其它广角迁移方法对于速度误差是敏感的。在地面反射地震学中很难精确地得到速度,因此限制在与陡峭的倾斜有关的结构中克希霍夫迁移的有效性。与地面地震获取相反,在VSP中可以从井探查精确地得到合理地精确的速度。通过良好的速度深度信息,克希霍夫深度迁移产生在井附近的更好的3-D深度图像。然而,在偏移2DVSP和3-DVSP中,源和接收机相对于地下成像点是不对称的。这种不对称性在计算加权因子时需要费很大的工作量。
与地震数据处理尤其是迁移技术的改进同时地,继续开发在获取的数据中相当基本的特性。近年来,多分量地震数据形成获取的地震数据的总量的增加的部分。这样做的理由是认识到传统的单分量地震数据主要响应于在地下在垂直方向上压缩波能量。传统的数据最通常是通过在水下环境中的压缩波源和水下听音器检测器或在地面地震获取中的垂直源和垂直检测器获取的。表示地下的岩石学和流体内容的附加信息是可以从有关剪切波在地下的传播的知识得到的。剪切波的到达最通常是通过具有不同于垂直的取向的接收机检测的。多分量记录的附加优点在于,即使对于压缩能量,接收信号的三个分量的知识可以提供在接收机处接收到能量的方向的指示,以及在该方向的能量总量。
Hokstad推导出用于预堆叠多分量克希霍夫迁移的方程。成像方程是根据粘滞弹性波理论推导得出的。由Hokstad推导的多分量成像方程的数学结构允许计算对于局部入射和散射波型(qP-qP,qP-qS1,qS1-qS1,等等)的所有的组合的分开的图像。
Hokstad的教导的限制在于,它们不解决3-D地震成像的真实世界问题。虽然由Hokstad得出的结果是非常好的,例子限于2D数据,以及不提供处理3-D多分量数据的任何实际的建议。在本发明中解决3-D多分量数据的迁移的问题。
发明内容
本发明是地下的地球地岩层的成像的方法。地震源在一个或多个源位置处被激活,生成地震波而成为地球地岩层。在一个或多个接收机位置处得到三维(3-D)地震数据。接收的地震数据包含有关地球的3-D结构的信息。完成3-C数据的3-D克希霍夫迁移。在克希霍夫迁移中,使用从每个源位置到多个图像点和从每个多个图像点到每个接收机位置的行进时间。源和接收机可以是在地面或地下位置。典型地,三分量基本上互相正交。在电缆实施方案中,用于三个分量的接收机可以是万向接头安装的。当在MWD环境下使用时,接收机可被安装在非旋转的套筒上,该套筒可被夹持在钻井壁。接收机可以是地音探测器或加速度表。
3-D迁移过程包括前向建模和成像操作。前向建模是计算地震波行进时间和在提出的速度模型内3-D空间中每个图像网格处波(射线)-方向角。这些行进时间和射线-方向角是对于每个源和接收机位置计算的。行进时间可以是对于压缩波或是对于剪切波的。
成像操作通过使用加权因子相加记录的反射的地震波能量到在它们的反射的位置。加权因子是射线几何关系、波行进时间、源-接收机孔径、小波相位和其它因素的函数。行进时间被使用来定位在记录的轨迹中的幅度(反射的能量)。迁移的输出图像幅度是一个代表地球物理反射率的标量值。
成像操作可以由现场处理器执行。替换地,处理可以在远端位置在与获取不同的时间完成。来自井地点的数据可以通过包括卫星链路的任何适当的装置或通过互联网连接被发送到远端位置。使得处理器能够接入多分量地震数据和执行3-C 3-D迁移处理的指令可以位于机器可读的存储器装置。这些指令使得处理器能够接入数据和处理数据。
附图说明
通过参考附图可以最好地了解本发明。图上相同的标号表示相同的单元,以及其中:
图1(现有技术)是用于获取VSP地震数据的典型的地面几何关系;
图2(现有技术)是实际的VSP地震数据的图;
图3a和3b显示描绘在本发明的一个实施例中实行的操作的框图;
图4a显示以x,y,z分量的、典型的3-C地震数据;
图4b是在3-D显示中图4a的数据;
图5显示图3a和3b的方法的缺点;
图6a和6b显示描绘在本发明的第二个实施例中实行的操作的框图;
图7a显示具有一个源和一个接收机的综合探测模型;
图7b显示对于图7a的模型的综合的3-C轨迹;
图8a-8c是通过使用传统的3-D迁移方法、在图7b中的数据的迁移的图像的3-D显示;
图9a-9c是通过使用V3-D迁移方法、在图7b中的数据的迁移的图像的3-D显示;
图10a-10c显示来自2D迁移、传统的3-D迁移和全向量3-C,3-D迁移场数据的结果的比较;以及
图11显示3-C万向接头安装的3-C接收机的例子。
具体实施方式
对于本发明,使用图1所示的现有技术系统的修改的版本。接收机12包括间隔开的接收机的阵列。典型地,使用5-80个接收机。每个接收机包括三分量(3-C)接收机。在本发明的一个实施例中,三个分量被表示为H1,H2和Z分量,Z分量是垂直分量,以及H1和H2轴与Z轴正交,并且互相正交。接收机可以是万向接头安装的。这使得在不同的钻孔中易于使用接收机。通过这样的安排,水平分量接收机的机械结构通常是与Z分量接收机的机械结构不同的,因为后者具有沿接收机的运动方向作用的重力。可以使用地音探测器或加速度表。在本发明的替换实施例中,三个接收机在灵敏度方面基本上是相同的,并且沿四面体的顶角取向。接收机的取向是通过使用现有技术已知的任何方法被确定的。
处理的基本部分是使用3-C 3-D向量克希霍夫预堆叠迁移。在实施迁移本身之前先讨论它。3-D预堆叠克希霍夫迁移通常被表示为:
其中M(x)是在3-D位置x处迁移的图像点,W(xs,xr,x)代表加权因子,或与探查几何关系、沿射线路径的速度、和地音探测器孔径有关的幅度补偿函数。xs和xr是在3-D中源和接收机位置。W(xs,xr,x)代表是与记录的反射波场Psr(ts(xs,x)+tr(x,xr))无关的,其中ts(xs,x)和tr(x,xr)是从源到图像位置和从图像位置到源的行进时间。记录的反射波场Psr是满足记录时间t=ts(xs,x)+tr(x,xr)对于源-接收机对来说是恒定的条件的散射能量的叠加。迁移处理过程把在时间t的记录的反射能量重新分布到其中可存在反射器的椭球(为了简化起见,在恒定的速度背景)表面。在每个位置x处迁移图像是对于到每个接收机xr的每个源xs在x处所有的椭球表面的加权的能量WPsr的叠加。传统的3-D预堆叠克希霍夫迁移无方向地映射在t的能量,意味着迁移算子均匀地分布相等的幅度到椭球上所有的点。如上所述,传统的克希霍夫迁移没有考虑从其反射能量的方向。这个特性将造成产生虚假的镜面图像反射。使用向量项允许向量3-C,3-D克希霍夫迁移克服这个问题。在迁移中使用的记录的波场Psr是标量,一个分量数据,或被旋转到三个分量数据的固定的方向。不使用总的反射波场。
向量3-C,3-D克希霍夫预堆叠迁移(V3-D迁移)是基于由公式(1)给出的、传统的克希霍夫累积。主要差别是迁移算子。V3-D迁移处理在每个时间采样点t处的反射波场作为3-D波场向量,而不是如在传统的3-D克希霍夫迁移中的标量值。进入的射线向量(与波前垂直)被引入到公式中。它是代表在接收机处从反射的图像点的射线方向的单位向量。当通过使用波场Psr(ts(xs,x)+tr(x,xr))构建在位置x处的图像时,我们只使用从x起源的波场数据。定向信息是通过使用3分量输入数据的动态极化分析而被确定的。假设在t时反射波场向量是
以及在位置xr处单位射线向量是
则迁移公式(1)成为
其中
A(xs,xr,x)是新的迁移算子。加权因子W(xs,xr,x)没有改变。
数据的处理是通过使用在图3a和3b上显示的、使用公式(1)的本发明的第一实施例完成的。3-C地震数据被旋转到预定的方向,成为1-C数据或把任何1-C(大多数是垂直分量)数据输入到迁移处理。
用于每个源和接收机位置的3-D行进时间表是通过使用用于地下的速度模型被生成101。对于参照图3a-3b讨论的一个方法,把地震数据预先旋转到通常相应于接收机到源方向100的方向。这在以后可被称为传统的方法。生成103对于每个接收机位置的行进时间表。在每个3-D空间图像位置x处对于所有的源-接收机对实行用加权因子的幅度的累积。幅度是在每个源-接收机对的总的行进时间t=ts(xs,x)+tr(x,xr)处的轨迹值。加权因子包括相位校正、射线轨迹几何校正、源-接收机孔径校正、和其它与能量损耗有关的因子。在每个网格点处的迁移图像是在轨迹的所有的贡献完成(109,111)后以上的累积113的求和结果。
以上讨论的处理3-C 3-D数据的方法具有下面讨论的几个缺点。首先,3-C反射数据包含来自所有的方向的能量,包括在如图4a-b所示的3-D空间中在每个接收机位置在每个记录采样时间的向量。使用单个分量数据或旋转在特定的方向旋转的数据将消除在迁移中来自所有的其它的方向的能量。这导致不完全的迁移图像。第二,各个分量迁移结果的代数求和/或使用代数求和的各个分量数据到迁移是不正确的。第三,接收的反射能量是定向的,应当被定向地分布到它的反射方向。在所有的方位角相等地分布反射能量将产生虚假的镜像图像。
图5显示用图3a-3b所示的处理方法的虚假的镜像图像的问题。图5上显示3-C,3-D VSP获取几何关系的平面图。源由250表示,以及三个分量接收机由221,223和225表示。由231表示的、用于从源250传播的地震能量的射线路径在由201表示的真实结构处被反射以及行进到接收机221,223和225。通过上述的方法,假设我们使用来自接收机部件221的反射数据,反射能量被相等地发布到每个位置201和203。203是虚假反射器。理由是从250到201和从201到221的行进时间等于从250到203和从203到221的行进时间。能量分布在传统的第一迁移方法中是非定向的。通过使用任何旋转的或非旋转的单分量数据,得到类似的结果。
为了解决以上讨论的问题,本发明的第二实施例迁移以向量形式的3-C数据和动态地旋转3-C数据。这被显示于图6a-6b。正如在本发明的第一实施例中那样,来自单次拍摄的数据被收集,以及对于这个源的相应的行进时间值ts(xs,x)是通过从速度模型进行前向建模(或从存储器恢复)而得到的。对于在3-D体积中每个输出位置得到行进时间。接着,对于在3-D体积中想要的输出图像点x的网格的选择的接收机位置得到行进时间值tr(x,xr)303。对于每个图像网格位置点x,从源到接收机的总的行进时间是。三个分量幅度是对于这些轨迹在时间t得到的307。3-C数据被旋转到图像到接收机射线方向309。当然,必须记住,射线方向通常相应于对于P波的最大幅度的方向,而对于剪切波,射线方向通常与最大幅度的方向正交。幅度和相位校正被应用到跟踪,以及产品被加到图像网格311。图6a上未示出但在克希霍夫迁移中隐含的是,对于图像网格上的多个图像点完成处理。
进行检验,以查看是否还有踪迹。如果是的话,则处理过程回到303。如果对于这个特定的源位置不再有踪迹,则进行检验,查看是否还有源位置301。如果还有源位置,则处理过程回到301。如果不再有源位置,则输出迁移的图像317。然后通过使用基本上相同的方法处理附加模式。
以上讨论的方法可以通过考虑对于压缩波和剪切波在速度域中的各向异性而被实施。行进时间计算然后通过使用各向异性速度被完成。
射线方向和粒子运动方向的更多的讨论例如在Postma的经典文章中给出。本发明的一个实施例执行对于横向各向异性媒体的3-C 3-D成像。当方位角各向异性(由于应变或断裂)被叠加到TI媒体时,引起稍微更复杂的情形。对于这样的媒体,弹性张量具有正交对称性。最通常类型的地球地岩层具有更复杂的弹性张量。虽然在理论上有可能通过这样的媒体进行射线跟踪(例如,参阅Crampin),但弹性张量的公式是有问题的。
本发明的方法在上面是参照在VSP探查中获取的数据讨论的。在理论上,方法也可被使用于在地面上记录的数据。
为了说明V3-D迁移比起传统的1-C,3-D克希霍夫预堆叠迁移的优点,我们检查记录的小波的脉冲响应。为了简化起见,在图7a所示的恒定的背景速度模型中使用单个源-接收机对。记录的数据包含两个3分量小波,一个在1秒以及另一个在2秒。为了清楚地说明迁移的定向效应,记录的小波的样本幅度在1秒仅仅在x方向,以及在2秒仅仅在y方向,如图7b所示。这分别相应于在x方向和y方向的反射器。这些轨迹的传统的、预先旋转的3-C或1-C数据迁移导致一个或两个椭球面,这取决于数据在迁移之前如何被旋转。在任何情形下,每个椭球面的幅度将具有恒定的幅度。这当然是由于在传统的3-D迁移中数据的非定向分布。能量被均匀分布到满足对于该样本的行进时间条件的所有的椭球表面点。图8a,b和c通过使用在最大能量方向上预先旋转的数据显示在x和y方向上3-D截面显示的迁移图像。图像幅度在x和y截面上是相同的。幅度的极性对于接收机位置是对称的,这意味着,能量映射是无方向性的。源和接收机位置被表示为501和503。
通过利用完全3分量数据,来自V3-D迁移的迁移的数据是不同的。对于在1秒的小波,迁移图像的幅度在x方向达到最大值,以及在y方向是零。对于在2秒的小波,幅度在y方向达到它的最大值,以及在x方向是零。见图9a-c。V3-D迁移在与接收机位置相反的方向上从正的图像产生相反的极性图像。这是正确的,因为在接收机的一个方向上具有正的反射率的反射器的反射响应等同于在接收机的相反的方向上具有负的反射率的反射器的响应。传统的3-D迁移将生成沿整个椭圆的相同极性的图像。
来自场的结果的例子显示于图10a-10c。图10a是在盐包附近获取的传统的2D数据的源-接收机面上的2D迁移。2D迁移图像仅仅是对于源-接收机剖面的3-D反射的2D投影。图10b上的图像是通过使用被预先旋转到最大能量(以上讨论的第一方法)从3-C 3-D VSP克希霍夫迁移得到的。对于3-D盐图像的迁移结果没有正确地解决离线盐反射;盐图像对于源-接收机是位置对称的。这是由于在传统的3-D迁移中反射数据的无方向性映射没有把数据映射到仅仅它的原先的位置。它对称地映射反射能量。用向量3-C,3-D VSP迁移来迁移3-C数据,显示垂直于源-接收机线的、盐体的横向变化(图10c),与当前的地理解释相一致。
本发明的方法是相对于在缆线上输送的接收机组件实行的VSP探查讨论的。然而,这并不是对本发明的方法的限制。本发明的方法也可以通过使用在井底组件(BHA)上输送的三分量接收机和地面地震处理被实行。
参考万向接头安装的接收机。万向接头安装的接收机的例子被显示于图9,其中三个传感器701,703,705被万向接头安装的倾斜的外壳707上。即使在外壳是倾斜的,传感器701仍能够保持垂直取向。这样的结构是必须的,以便得到在固定的参考坐标系统上地震场的三个分量。在钻井内机壳的取向可以由诸如磁场计那样的适当的取向传感器确定。
Baker Hughes公司具有多电平接收机(MLR),它可以被配置成从1到13电平。这大大地加速数据获取。井下接收机可以与其它测井业务相组合地运行,通过缆线或者管道输送,减小到井中的行程数目并且节省吊装时间。在高角度井中,井下接收机可以在钻井管道或螺旋管中被输送,以及也与各种各样的开孔测井业务相组合地运行,大大地减小吊装时间。
上述的3C-3D向量迁移方法可以在通用数字计算机上实施。正如本领域技术人员将会看到的,用于计算机的指令放置在机器可读的存储器装置中,诸如ROM、EPROM、EAROM、快闪存储器和光盘。这些可以是计算机的一部分或可以通过适当的通信信道被链接到计算机,以及甚至可以在远端位置。同样地,以上讨论的那种多分量数据可被存储在计算机或可以通过适当的通信信道被链接到计算机。通信信道可以包括互联网,使得用户能够接入来自一个远端位置的数据以及得到来自另一个远端位置的指令来处理数据。在机器可读的存储器装置上的指令使得计算机能够接入多分量数据和按照上述的方法来处理数据。
虽然上述的公开内容是针对本发明的优选实施例的,但本领域技术人员将会看到各种修改方案。所以上述的公开内容打算包括在所附权利要求的范围和精神内的所有的这样的变化。
Claims (46)
1.一种估计地球地岩层的方法,包括:
(a)激活在至少一个源位置处的地震源以及生成进入到地球地岩层的地震波,
(b)在至少一个接收机位置处接收地震数据的三个分量,一部分接收的地震数据是由生成的地震波与不是在通过至少一个源位置和所述至少一个接收机位置的垂直面上的一部分地球地岩层的交互作用得到的,
(c)通过使用从所述至少一个源位置到地球地岩层中的至少一个图像点的行进时间和从所述至少一个图像点到至少一个接收机位置的行进时间,同时处理地震数据的三个分量。
2.权利要求1的方法,其中所述至少一个接收机位置处在(i)地面位置和(ii)井下位置中的一个位置。
3.权利要求1的方法,其中至少一个源位置处在地面位置。
4.权利要求1的方法,其中地震数据的三个分量基本上互相正交。
5.权利要求1的方法,其中部分地球地岩层还包括反射界面。
6.权利要求1的方法,还包括使用用于得到行进时间的速度模型。
7.权利要求1的方法,还包括执行用于得到行进时间的查找表。
8.权利要求1的方法,其中所述生成的地震波包括(i)压缩波,和(ii)剪切波的至少一个波。
9.权利要求1的方法,其中所述一部分接收的地震波数据对应于(i)压缩波,和(ii)剪切波的至少一个波。
10.权利要求1的方法,其中处理地震数据的三个分量的步骤还包括以特定的方向旋转接收的地震数据的三个分量。
11.权利要求1的方法,其中处理地震数据的三个分量的步骤还包括以特定的方向旋转所述地震数据的三个分量的向量和。
12.权利要求10的方法,其中所述特定的方向是由通过所述至少一个源位置和所述至少一个接收机位置的垂直面定义的水平方向。
13.权利要求10的方法,其中所述特定的方向是由将所述至少一个图像点连接到所述至少一个接收机位置的线定义的水平方向。
14.权利要求11的方法,其中所述特定的方向是由通过所述至少一个源位置和所述至少一个接收机位置的垂直面定义的水平方向。
15.权利要求11的方法,其中所述特定的方向是由将所述至少一个图像点连接到所述至少一个接收机位置的线定义的水平方向。
16.权利要求1的方法,其中所述至少一个图像点还包括多个图像点。
17.权利要求16的方法,其中多个图像点是在输出点的网格上。
18.权利要求1的方法,其中所述处理步骤还包括施加作为(i)幅度校正和(ii)相位校正的至少一项的校正。
19.权利要求18的方法,其中所述处理步骤还包括把在校正后的接收的地震数据的贡献加到在所述至少一个图像点处的输出上。
20.权利要求19的方法,其中所述加上贡献步骤是对于每个至少三个分量的每个分开地完成的。
21.权利要求1的方法,其中至少一个源位置还包括多个附加源位置,所述方法还包括对于多个附加源位置重复进行(a)--(c)。
22.权利要求1的方法,其中至少一个接收机位置还包括多个附加接收机位置,方法还包括对于多个附加接收机位置重复进行(a)--(c)。
23.权利要求19的方法,其中至少一个源位置还包括多个附加源位置,所述方法还包括对于多个附加源位置重复进行(a)--(c)。
24.权利要求19的方法,其中至少一个接收机位置还包括多个附加接收机位置,所述方法还包括对于多个附加接收机位置重复进行(a)--(c)。
25.权利要求20的方法,其中至少一个源位置还包括多个附加源位置,所述方法还包括对于多个附加源位置重复进行(a)--(c)。
26.权利要求28的方法,其中至少一个接收机位置还包括多个附加接收机位置,方法还包括对于多个附加接收机位置重复进行(a)--(c)。
27.一种用于估计地球地岩层的系统,包括:
(a)在至少一个源位置处的地震源,其生成进入所述地球地岩层的地震波,
(b)至少一个传感器,接收在至少一个接收机位置处接收地震数据的三个分量,一部分接收的地震数据是由生成的地震波与不是在通过至少一个源位置和所述至少一个接收机位置的垂直面上的一部分地球地岩层的交互作用得到的,
(c)处理器,其通过使用从至少一个源位置到地球地岩层中的至少一个图像点的行进时间和从至少一个图像点到至少一个接收机位置的行进时间,同时处理接收的地震数据的三个分量。
28.权利要求27的系统,其中至少一个传感器包括在井下工具上的多个间隔开的传感器。
29.权利要求27的系统,其中所述三个分量是基本上互相正交。
30.权利要求27的系统,其中至少一个传感器还包括在井下工具上的三个万向接头安装的传感器。
31.权利要求27的系统,其中处理器使用(i)速度模型和(ii)表格的至少一个,用于得到行进时间。
32.权利要求27的系统,其中处理器以特定的方向旋转所述得到的地震数据的三个分量。
33.权利要求27的系统,其中处理器以特定的方向旋转地震数据的三个分量的向量和。
34.权利要求33的系统,其中所述特定的方向是由通过该至少一个源位置和该至少一个接收机位置的垂直面定义的方向。
35.权利要求27的系统,其中处理器对地震数据的至少三个分量施加(i)幅度校正和(ii)相位校正的至少一项校正。
36.权利要求35的系统,其中所述处理器还把所述得到的地震数据的贡献加到在至少一个图像点处的输出上。
37.权利要求36的系统,其中处理器对于每个所述三个分量分开地加上所述贡献。
38.权利要求27的系统,其中得到所述三个分量还包括使用在被使用来钻井的井底组件的非旋转的套筒上的三个传感器。
39.一种用于估计地球地岩层的系统,包括:
(a)计算机,对于从由在至少一个源位置处的源生成的地震波与地球地岩层的交互作用得到的、在至少一个接收机位置处接收的、多分量地震数据执行向量3-D迁移;
(b)机器可读的存储器装置,包括使得计算机能够执行向量3-D迁移的指令。
40.权利要求39的系统,其中多分量地震数据部分是来自于在由至少一个源和至少一个接收机定义的弧矢平面外的地球地岩层的一部分。
41.权利要求39的系统,其中指令还使得能够以特定的方向旋转多分量地震数据的向量和。
42.权利要求41的系统,其中特定的方向是由将所述至少一个图像点连接到所述至少一个接收机位置的线定义的。
43.权利要求41的系统,其中特定的方向是由所述至少一个图像点连接到所述至少一个源位置的线定义的。
44.权利要求39的系统,其中机器可读的存储器装置是从包含(i)ROM,(ii)EPROM,(iii)EAROM,(iv)快闪存储器,和(v)光盘的组中选择的。
45.权利要求39的系统,其中机器可读的存储器装置处在与计算机的位置不同的位置。
46.权利要求39的系统,其中计算机处在与其中存储多分量数据的位置不同的位置。
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US51510703P | 2003-10-28 | 2003-10-28 | |
US60/515,107 | 2003-10-28 | ||
US57652604P | 2004-06-03 | 2004-06-03 | |
US60/576,526 | 2004-06-03 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1886678A true CN1886678A (zh) | 2006-12-27 |
CN100401101C CN100401101C (zh) | 2008-07-09 |
Family
ID=34555954
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2004800354655A Expired - Fee Related CN100401101C (zh) | 2003-10-28 | 2004-10-28 | 向量3分量3维克希霍夫预堆栈迁移 |
Country Status (7)
Country | Link |
---|---|
US (2) | US7542373B2 (zh) |
EP (1) | EP1678527B1 (zh) |
CN (1) | CN100401101C (zh) |
CA (1) | CA2544254C (zh) |
EA (1) | EA009638B1 (zh) |
NO (1) | NO340025B1 (zh) |
WO (1) | WO2005043194A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102395902A (zh) * | 2009-04-16 | 2012-03-28 | 兰德马克图形公司 | 使用快速面向目标照明计算的地震成像系统及方法 |
CN104216007A (zh) * | 2013-06-05 | 2014-12-17 | 中国石油天然气集团公司 | 二维垂直地震剖面和三维地面数据的空间同步成像处理方法 |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7937224B2 (en) * | 2006-05-17 | 2011-05-03 | Westerngeco L.L.C. | Diplet-based seismic processing |
US7952960B2 (en) * | 2006-10-03 | 2011-05-31 | Bp Corporation North America Inc. | Seismic imaging with natural Green's functions derived from VSP data |
GB2459052B (en) * | 2006-12-26 | 2011-05-25 | Baker Hughes Inc | Imaging near-borehole reflectors using shear wave reflections from a multi-component acoustic tool |
US7508736B2 (en) * | 2007-03-09 | 2009-03-24 | Baker Hughes Incorporated | Vector migration of 1st order free-surface related downgoing multiples from VSP data |
EP2126686A4 (en) * | 2007-03-20 | 2017-10-04 | ION Geophysical Corporation | Apparatus and method for processing geophysical information |
US8010293B1 (en) * | 2007-10-29 | 2011-08-30 | Westerngeco L. L. C. | Localized seismic imaging using diplets |
CA2703588C (en) * | 2007-12-12 | 2015-12-01 | Exxonmobil Upstream Research Company | Method and apparatus for evaluating submarine formations |
WO2009077440A2 (en) * | 2007-12-14 | 2009-06-25 | Shell Internationale Research Maatschappij B.V. | Method of processing data obtained from seismic prospecting |
US8395967B2 (en) * | 2008-06-11 | 2013-03-12 | Baker Hughes Incorporated | Vector migration of virtual source VSP data |
US8699298B1 (en) | 2008-06-26 | 2014-04-15 | Westerngeco L.L.C. | 3D multiple prediction and removal using diplets |
US20100067328A1 (en) * | 2008-09-17 | 2010-03-18 | Andrew Curtis | Interferometric directional balancing |
US20110093203A1 (en) * | 2009-10-21 | 2011-04-21 | Chevron U.S.A. Inc. | System and method for clustering arrivals of seismic energy to enhance subsurface imaging |
US9036450B2 (en) * | 2010-04-12 | 2015-05-19 | Schlumberger Technology Corporation | Generating an image of a subterranean structure |
US8773951B2 (en) * | 2011-03-18 | 2014-07-08 | Chevron U.S.A. Inc. | System and method for seismic imaging with reduced computational cost |
US10977396B2 (en) | 2012-01-13 | 2021-04-13 | Schlumberger Technology Corporation | Determining an elastic model for a geologic region |
WO2014197464A1 (en) * | 2013-06-03 | 2014-12-11 | Westerngeco Llc | Processing of seismic data |
CN104502972B (zh) * | 2014-12-11 | 2017-04-05 | 中国石油集团东方地球物理勘探有限责任公司 | 三分量地震波积分偏移方法和装置 |
CN110161562B (zh) * | 2019-04-29 | 2020-09-18 | 中南大学 | 倾斜横向各向同性介质中的裂缝弱度反演方法及系统 |
US11391855B2 (en) * | 2020-03-13 | 2022-07-19 | Saudi Arabian Oil Company | Developing a three-dimensional quality factor model of a subterranean formation based on vertical seismic profiles |
WO2023059688A1 (en) * | 2021-10-06 | 2023-04-13 | Bp Corporation North America Inc. | Method and apparatus for implementing a high-resolution seismic pseudo-reflectivity image |
Family Cites Families (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4627036A (en) * | 1982-10-08 | 1986-12-02 | Phillips Petroleum Company | Vertical seismic profiling |
US4635238A (en) * | 1984-09-12 | 1987-01-06 | Phillips Petroleum Company | Data processing method for correcting P and S wave seismic traces |
US4802147A (en) * | 1985-05-23 | 1989-01-31 | Mobil Oil Corporation | Method for segregating and stacking vertical seismic profile data in common reflection point bins |
US4809239A (en) * | 1987-07-14 | 1989-02-28 | Schlumberger Technology Corporation | Method for evaluating parameters related to the elastic properties of subsurface earth formations |
US4893694A (en) * | 1988-11-14 | 1990-01-16 | Mobil Oil Corporation | VSP-based method and apparatus for tieing seismic data shot using different types of seismic sources |
US4918670A (en) * | 1989-02-15 | 1990-04-17 | Conoco Inc. | Method for interval velocity analysis and determination of reflector position from before stack seismic data |
US5062086A (en) * | 1990-08-27 | 1991-10-29 | Conoco Inc. | Calculation of raypaths and wavepaths from traveltime tables for the tomographic estimation of transmission velocities |
US5142501A (en) * | 1991-06-13 | 1992-08-25 | Chevron Research And Technology Company | Method of geophysical exploration by analyzing shear-wave polarization directions |
US5111437A (en) * | 1991-08-05 | 1992-05-05 | Exxon Production Research Company | Method for positioning seismic sources to flatten CDP fold in VSP surveys |
US6002642A (en) * | 1994-10-19 | 1999-12-14 | Exxon Production Research Company | Seismic migration using offset checkshot data |
US5790473A (en) * | 1995-11-13 | 1998-08-04 | Mobil Oil Corporation | High fidelity vibratory source seismic method for use in vertical seismic profile data gathering with a plurality of vibratory seismic energy sources |
US6081482A (en) * | 1996-05-09 | 2000-06-27 | 3Dgeo Development, Inc. | Semi-recursive imaging under complex velocity structures |
FR2759172B1 (fr) * | 1997-02-05 | 1999-03-05 | Inst Francais Du Petrole | Methode de traitement de donnees sismiques de puits multi-composantes orientees |
US6553315B2 (en) * | 1997-10-15 | 2003-04-22 | Albin K. Kerekes | Seismic imaging using omni-azimuth seismic energy sources and directional sensing |
US6131694A (en) * | 1998-09-02 | 2000-10-17 | Ahlliburton Energy Services, Inc. | Vertical seismic profiling in a drilling tool |
FR2791438B1 (fr) * | 1999-03-23 | 2001-08-24 | Schlumberger Services Petrol | Procede de controle de qualite d'une migration de donnees sismiques utilisant la transformee de radon generalisee |
US6292754B1 (en) * | 1999-11-11 | 2001-09-18 | Bp Corporation North America Inc. | Vector recomposition of seismic 3-D converted-wave data |
GB2362467B (en) * | 2000-05-18 | 2004-03-31 | Schlumberger Ltd | A method of processing seismic data |
US6950786B1 (en) * | 2000-10-10 | 2005-09-27 | Schlumberger Technology Corporation | Method and apparatus for generating a cross plot in attribute space from a plurality of attribute data sets and generating a class data set from the cross plot |
US6611764B2 (en) * | 2001-06-08 | 2003-08-26 | Pgs Americas, Inc. | Method and system for determining P-wave and S-wave velocities from multi-component seismic data by joint velocity inversion processing |
FR2831961B1 (fr) * | 2001-11-07 | 2004-07-23 | Inst Francais Du Petrole | Methode de traitement de donnees sismiques de puits en amplitude preservee absolue |
US6748330B2 (en) * | 2002-04-10 | 2004-06-08 | Schlumberger Technology Corporation | Method and apparatus for anisotropic vector plane wave decomposition for 3D vertical seismic profile data |
US6904368B2 (en) * | 2002-11-12 | 2005-06-07 | Landmark Graphics Corporation | Seismic analysis using post-imaging seismic anisotropy corrections |
US7508736B2 (en) * | 2007-03-09 | 2009-03-24 | Baker Hughes Incorporated | Vector migration of 1st order free-surface related downgoing multiples from VSP data |
-
2004
- 2004-10-25 US US10/972,880 patent/US7542373B2/en not_active Expired - Fee Related
- 2004-10-28 CA CA2544254A patent/CA2544254C/en not_active Expired - Fee Related
- 2004-10-28 WO PCT/US2004/035877 patent/WO2005043194A1/en active Application Filing
- 2004-10-28 CN CNB2004800354655A patent/CN100401101C/zh not_active Expired - Fee Related
- 2004-10-28 EP EP04810092.9A patent/EP1678527B1/en not_active Expired - Fee Related
- 2004-10-28 EA EA200600830A patent/EA009638B1/ru not_active IP Right Cessation
-
2006
- 2006-04-26 NO NO20061848A patent/NO340025B1/no not_active IP Right Cessation
-
2009
- 2009-05-29 US US12/475,264 patent/US7764571B2/en not_active Expired - Fee Related
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102395902A (zh) * | 2009-04-16 | 2012-03-28 | 兰德马克图形公司 | 使用快速面向目标照明计算的地震成像系统及方法 |
CN102395902B (zh) * | 2009-04-16 | 2015-04-01 | 兰德马克图形公司 | 使用快速面向目标照明计算的地震成像系统及方法 |
CN104216007A (zh) * | 2013-06-05 | 2014-12-17 | 中国石油天然气集团公司 | 二维垂直地震剖面和三维地面数据的空间同步成像处理方法 |
CN104216007B (zh) * | 2013-06-05 | 2017-05-03 | 中国石油天然气集团公司 | 二维垂直地震剖面和三维地面数据的空间同步成像处理方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2005043194A1 (en) | 2005-05-12 |
CN100401101C (zh) | 2008-07-09 |
US20050135189A1 (en) | 2005-06-23 |
EP1678527A1 (en) | 2006-07-12 |
NO20061848L (no) | 2006-06-22 |
CA2544254C (en) | 2013-02-12 |
CA2544254A1 (en) | 2005-05-12 |
EP1678527B1 (en) | 2019-11-27 |
EA200600830A1 (ru) | 2006-10-27 |
EA009638B1 (ru) | 2008-02-28 |
US7542373B2 (en) | 2009-06-02 |
US7764571B2 (en) | 2010-07-27 |
NO340025B1 (no) | 2017-02-27 |
US20090323471A1 (en) | 2009-12-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100401101C (zh) | 向量3分量3维克希霍夫预堆栈迁移 | |
US7355923B2 (en) | Seismic analysis using post-imaging seismic anisotropy corrections | |
US8457899B2 (en) | Method of processing data obtained from seismic prospecting | |
CN102282481B (zh) | 基于地震能见度分析的数据采集和叠前偏移 | |
JP2020522699A (ja) | 地下構造の検出 | |
US8045417B2 (en) | Analyzing 2-D surface and/or borehole seismic data to locate subsurface diffractors | |
EA026344B1 (ru) | Система и способ получения и обработки сейсмических данных о полях упругих волн | |
AU2008216238B2 (en) | Performing 3-D scatterer imaging from 2-D seismic data | |
Bellefleur et al. | Downhole seismic imaging of a massive sulfide orebody with mode-converted waves, Halfmile Lake, New Brunswick, Canada | |
MX2014010606A (es) | Extraccion de datos de cizallamiento sv a partir de datos marinos de la onda p. | |
GB2321968A (en) | Seismic data processing method | |
US5142501A (en) | Method of geophysical exploration by analyzing shear-wave polarization directions | |
US20120269035A1 (en) | Evaluating Prospects from P-Wave Seismic Data Using S-Wave Vertical Shear Profile Data | |
CA1311830C (en) | Method for extending the lateral subsurface coverage in vsp surveys | |
MX2014005103A (es) | Extraccion de datos de cizallamiento sv a partir de datos sismicos de onda p. | |
Mari et al. | Seismic well surveying | |
Liu et al. | Depth imaging of multicomponent seismic data through the application of 2D full‐waveform inversion to P‐and SH‐wave data: SEAM II Barrett model study | |
Ganssle | Seismic explorations | |
US20080109169A1 (en) | Method and system for characterizing seismic reflection points | |
Li et al. | Inversion techniques for P-wave anisotropic parameters from VSP surveys | |
O'brien | Aspects of seismic reflection prospecting for oil and gas | |
Kumar | Seismic data processing: velocity analysis & migration |
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: 20080709 Termination date: 20201028 |
|
CF01 | Termination of patent right due to non-payment of annual fee |