CN1271420C - 海洋地震数据处理方法和地震测量方法 - Google Patents
海洋地震数据处理方法和地震测量方法 Download PDFInfo
- Publication number
- CN1271420C CN1271420C CNB028052773A CN02805277A CN1271420C CN 1271420 C CN1271420 C CN 1271420C CN B028052773 A CNB028052773 A CN B028052773A CN 02805277 A CN02805277 A CN 02805277A CN 1271420 C CN1271420 C CN 1271420C
- Authority
- CN
- China
- Prior art keywords
- wave mode
- wave
- seismic energy
- seismic
- layer
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 71
- 238000006243 chemical reaction Methods 0.000 claims abstract description 72
- 230000001902 propagating effect Effects 0.000 claims abstract description 33
- 239000002245 particle Substances 0.000 claims description 71
- 230000005540 biological transmission Effects 0.000 claims description 22
- 238000006073 displacement reaction Methods 0.000 claims description 17
- 230000000644 propagated effect Effects 0.000 claims description 15
- 238000012545 processing Methods 0.000 claims description 14
- 230000011514 reflex Effects 0.000 claims description 8
- 238000011282 treatment Methods 0.000 claims description 7
- 239000002344 surface layer Substances 0.000 abstract 3
- 239000010410 layer Substances 0.000 abstract 2
- 239000011435 rock Substances 0.000 description 57
- 230000003068 static effect Effects 0.000 description 51
- 230000006870 function Effects 0.000 description 16
- 230000008859 change Effects 0.000 description 14
- 238000005516 engineering process Methods 0.000 description 11
- 238000005259 measurement Methods 0.000 description 8
- 230000001133 acceleration Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000001615 p wave Methods 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 239000002360 explosive Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000002310 reflectometry Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- 101100400452 Caenorhabditis elegans map-2 gene Proteins 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 101100001672 Emericella variicolor andG gene Proteins 0.000 description 1
- 244000089486 Phragmites australis subsp australis Species 0.000 description 1
- 235000014676 Phragmites communis Nutrition 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
- 230000009466 transformation 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
- G01V1/286—Mode conversion
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
确定海底一层,特别是表面或近表面层(5)的地震性质的方法,包含使以第一波型传播的地震能量指向该层的边界面,从而在边界面引起地震能量的部分波型转换。例如,当地震能量向上传播穿过表面或近表面层(5)和基岩(6)之间的界面时,由于表面或近表面层(5)和基岩(6)之间地震性质的差异,便可能发生部分波型转换。在本发明中,在接收器处接收地震能量的两种波型,即初始波型和在界面处由波型转换产生的波型。由接收器所采集的地震数据确定两种波型在界面和接收器之间的走时之差。
Description
技术领域
本发明涉及海洋地震数据处理方法,特别是涉及对这样的海洋地震数据的处理方法,即这些数据包括两种地震能量传播波型(mode),一种波型是在地层边界的部分波型转换的结果,本发明的方法提供对静态时移的估计,该静态时移是由在地球表面或近地球表面的与其下面地层有不同地震特性的一个地层产生的。本发明还涉及包括处理以上述方式采集的地震数据的地震测量方法。
背景技术
图1是海洋、海底地震测量布置的示意图。在这一测量布置中,地震能量由震源1发射,震源1是从拖船2悬挂在海面下。地震能量是向下发射,穿入地球内部,并被作为反射体3的地质特征反射。被反射的地震能量向上穿过地球内部,进入海洋并由放在海底上的地震接收器4检测到。通过例如确定从震源1到接收器4的地震能量走时,得到关于地球内部的信息。在图1中只画出了一个反射体3,但实际上在地球内部的若干个地质特征将作为地震能量的部分反射体。再有,在图1中只画出了一个源1和一个接收器4,但实际上海洋地震测量布置可以有一个震源阵列和一个接收器阵列。
地球的地质构造不是均一的。在处理海洋地震数据时的一个问题是经常有一个在表面或近表面的一层5,它的特性与其下伏地质构造6(下文中称作“基岩”)的特性显著不同。例如,如果在地表或近地表的一层不如基岩那样固结,便能发生这种情况。特别是,在表面或近表面层5中的地震能量速度可能显著低于基岩6中的速度,这样,这种地表或近地表层通常被称作“低速层”(或简写为LVL)。与如果地表或近地表层和基岩有完全相同的地震特性因而会记录到的走时相比,这一速度差将产生地震能量走时的偏移,这些走时偏移通常被称作“静态时移”,或只称作“静移(statics)”。
图1中所示低速层5是作为表面层,但它不需要延伸到表面,可以再有一层覆盖该低速层。
由表面或近表面低速层5产生的静态时移依赖于该层的厚度和地震能量穿过该层的传输速度。通常低速层5的厚度和穿过该层的传播速度都会发生横向变化,所以在一个位置的地震接收器处观测到的静态时移大多不同于在另一个位置的接收器处观测到的静态时移。作为一级近似,在一个接收器记录的整个数据集将相对在另一接收器记录的数据提前或延迟一个静态时间偏移。
当处理地震数据时,特别希望考虑静态时移。除非从地震数据中去掉这些静态时移,否则就会存在不确定性,诸如来自较深层的地震事件的到时的变化是由于这些较深地层的测度或横向位置的变化造成的还是简单地由于低速层5中的传播效应造成的。
发明内容
本发明提供一种方法,用于处理包括相应的地震能量第一和第二波型的地震数据,该方法包含处理地震数据以得到以第一波型传播的地震能量和以第二波型传播的地震能量穿过该地层的走时差,所述第一波型地震能量和第二波型地震能量沿着实质上相同的路径穿过该地层。
在一个优选实施例中,第二波型是由第一波型在海底的一个层界面处的部分波型转换产生的。
从图1将看到,地震能量两次穿过低速层5:一次是它从震源1向反射体3传播,再次是它以反射体3向接收器4传播。两次穿过低速层5将造成静态时移;由向下传播的地震能量穿过低速层5造成的时移称作震源一侧的静态时移,由向上传输的地震能量穿过低速层5造成的时称称作接收器一侧的静态时移。本发明提供估计接收器一侧静态时移的方法。
本发明的方法可用于处理先前存在的地震数据。或者,它可以被纳入地震测量方法用于在数据获取时或在其后处理地震数据。
本发明的一个实施例假定基岩6和低速层5的物理性质反差足够大,以致在地震能量向上传播穿过基岩6和低速层5之间的边界时,能发生能量传播的P波型到能量传播的S波型的显著转换以及反过来从S波型到P波型的显著转换。这种波型转换的结果是,对于低速层和基岩边界上方的地震接收器所记录的每个P波型事件,必须有一个相应的S波型事件。该P波型事件和相应的S波型事件将发生在不同时刻,因为在低速层中P波型和S波型和传播速度彼此不相等。然而,由于波型转换是在低速层的下边界发生的,所以P波型和相应的S波型事件之间的时间延迟必须是由于两种波型在低速层中不同速度产生的。P波型事件和相应的S波型事件之间的时间延迟将不会受到基岩性质的显著影响。这样,通过在地震数据中识别P波型事件及相应的S波型事件并确定该P波型事件和相应的S波型事件之间的时间延迟,便能容易地确定穿过低速层的P波形和S波型的走时差。
当向下传播的波在低速层5和基岩6之间的界面受到反射时也可以发生波型转换。本发明可以应用于这一情况,因为当在反射时发生波型转换时,被反射的信号将含有P分量和S分量,这两个分量穿过低速层5时将有不同的走时。
当向下传播的波在低速层5和基岩6之间的界面受到临界折射时也可以发生波型转换。沿界面传播的波将在低速层中激发P波型和S波型两种波,本发明能应用于这些P和S波型。
对于指示地震数据一个方面的任何参数,如压强或质点运动分量(术语“质点运动”包括质点位移、质量速度、质点加速度以及质点位移的更高阶导数),P波型和S波型之一的事件通常比其他波型的事件表现得更强。所以,根据本发明的优选实施例,指示地震数据两个不同方面的两个参数用于定位P波型事件和相应的S波型事件。原则上,如果在低速层的下边界发生足够强的波型转换,则通过对所收到的地震数据的检查,能够确定相应的P波型和S波型事件对的位置。然而,在许多情况下,波型转换后的事件振幅低,于是优选互相关技术或反褶积技术。
与S波静态时移相比,P波静态时移通常较小,这是因为P波速度大于S波速度,因而走时变化较小。再有,在低速层5中发现的欠固结的海底沉积物趋向于它们所造成的P波速度横向变化不如它们造成的S波速度横向变化那样显著。在许多情况下,P波静态时移是如此之小,以致能安全地假定为可以忽略。
本发明使得有可能估计每个接收器处S波静态时移和P波静态时移之差。如果P波静态时移已知、能被估计出来或能被假定为零,则本发明使得有可能得到S波静态时移。
本发明的第二方面提供一种地震测量方法,包含如下步骤:使以第一波型传播的地震能量定向到一个海底层的边界面,从而在边界面发生地震能量向第二波型的部分波型转换;在一个或多个接收器荻取包括第一和第二波型地震能量的地震数据;以及根据上文定义的方法处理这些地震数据。
本发明的第三方面提供一种设备用于处理包括第一和第二波型地震能量的地震数据,该设备包含用于处理地震数据的装置以得到以第一波型传播的地震能量和以第二波型传播的地震能量穿过该层的走时差,所述第一波型地震能量和第二波型地震能量沿着实质上相同的路径穿过该层。
在一个优选实施例中,该设备包含一个可编程数据处理器。
本发明的第四方面提供一种存储介质,含有上文定义的设备的数据处理器中使用的程序。
附图说明
现在将参考附图以举例方式详细描述本发明的优选实施例,其中:
图1是海底地震测量的示意图;
图2示意性显示上行P波在基岩和低速层之间边界的波型转换;
图3显示图2的情况中地震数据道(trace)的示意性举例;
图4显示图3中各数据道的互相关图;
图5(a)图5(b)显示对表1模型产生的合成地震数据;
图6(a)和6(b)显示图5(a)和(b)数据的反褶积;
图7显示PP-PS走时差的理论和实验结果;
图8是具有同炮检距震源和接收器几何分布的地震测量布局示意图;
图9(a)和9(b)显示使用图8的地震测量布局得到的质点速度垂直分量和径向分量;
图10(a)和10(b)显示经过初步处理后的图9(a)和9(b)的质点速度垂直分量和径向分量;
图11示意性显示上行S波在基岩和低速层之间边界的波型转换;
图12(a)示意性显示下行P波在低速层和基岩之间边界反射时的波型转换;
图12(b)示意性显示下行S波在低速层和基岩之间边界反射时的波型转换;
图13(a)示意性显示沿低速层和基岩之间边界传播的P波的波型转换;
图13(b)示意性显示沿低速层和基岩之间边界传播的S波的波型转换;以及
图14是根据本发明的设备的示意方框图。
具体实施方式
现在将参考一个实施例描述本发明的原理,该实施例中由部分波型转换产生一对相应的事件。图2显示上行P波穿过基岩6和低速层5之间界面时的波型转换。上行P波7部分地作为P波7′穿过,但也受到部分转换,成为S波7″。假定上行P波是由地球内部两层之间的深界面反射或由通过较低层的折射产生的。在基岩6和低速层5之间的界面处也将产生折射,透射P波7′的折射角将不同于转换S波7″的折射角。
透射P波7′和转换S波7″二者都将在接收器阵列产生信号。如图2中指出的那样。透射P波7′和转换S波7″是沿不同方向传播,因而将到达不同的接收器位置,而不是在同一接收器。然而,来自相邻转换点的相似的转换S波将到达图2中所示透射P波7′的同一接收器。对于在一特定接收器接收的透射P波和转换S波,将假定转换S波的转换点(在低速层的底面)和P波穿过低速层底面的那一点之间的差异不会显著影响由本发明方法得到的静态时移估计值。
透射P波7′将主要记录在接收器记录的地震能量垂直分量上,而转换S波7″将主要记录在接收器4记录的地震能量径向分量上。(径向方向是投射在海底的震源-接收器方向,这一方向将定义为X方向。)这是因为对于地球内部的典型结构和目标反射体的典型深度,在实际地震测量布局中从地球内部深处向上传播的波趋于与垂直方向成30°或更小的角。横波(S波)主要记录在径向分量上,这是因为横波质点运动垂直于传播方向,与此相比,P波(纵波)的质点运动是沿传播方向。
图3是在图2所示部分波型转换情况下指示接收器4处记录地震数据一些方面的两个参数的示意图。这两个参数是在接收器处测量的质点速度的径向分量(X分量)。和在接收器处测量的质点速度的垂直分量(Z分量)。由于S波7″是作为上行P波7的转换波产生的,如果基岩和低速层之间的界面局部地近似为平面,则所记录的S波将含有与P到达信号相似的子波。S波到达将不与P波到达同时发生,而是在时间延迟dt(与P波到达比较)之后到达,时间延迟dt是由低速层5中P波速度和S波速度之差造成的。这一时间延迟dt正是P波静态时移和S波静态时移之差。这样,如果P波静态时移已知或能准确地估计出来,则只要把时间延迟dt加到P波静态时移上便能得到S波静态时移。
原则上,时间延迟dt可以通过比较P波到时和S波到时的任何方法得到。然而,在本发明的优选实施例中,时间延迟dt是通过对地震接收器记录的地震能量垂直分量和水平分量进行反褶积或互相关确定的。
图4显示通过对图3所示质点速度的水平和垂直分量进行互相关得到的互相关图。可以看到,在零时移处有一个波峰,这是因为在所测量的质点速度的X分量中有一个非零振幅(尽管它主要发生在所测量的质点速度的Z分量中)。在所测量的质点速度垂直分量中出现的任何S波能量也将在互相关图中的零时移处造成波峰。
在互相关图中时间dt处的波峰来自所测量的质点速度垂直分量中的P波到达波峰和所测量的质点速度X分量中的S波到达波峰。为了准确地区分这一波峰和零时移处的波峰,波峰不重叠是人们所希望的(尽管不是至关重要的)。如果在记录数据中的子波有足够高的频率而且低速层5的厚度大于2个子波波长,则波峰将不重叠。如果不能肯定低速层5的厚度足以厚到能防止波峰重叠,通过对所测量的质点速度的X分量和Z分量进行反褶积,时移的确定准确度能够提高。或者,所记录的子波能通过下列技术完全分解成它的P波分量和S波分量;例如共同待决英国专利申请0001355.7号或0003406.6号,Robertsson等在“在陆地表面地震记录中使用密集布设的3C单个传感器组进行的波场分离和近地表性质估计”(勘探地球物理协会(SEG)第70届年会,Calgary(2000))或Curtis等在“在陆地地震勘探中的波场分离和近地表性质估计”(62届EAGE大会,Glasgow,详细摘要(2000)中公开的技术。这种分解技术具有从所测量的质点速度的X分量和Z分量中去掉图3所示“泄漏”波峰的作用。如果再对分解后的记录波场P分量和S分量进行互相关,则与图4所示结果相比,在零时延处的波峰应不复存在而在P-S延时dt处的波峰应得到增强。
如果地震波主要是沿径向传播(例如不是反射回负的径向方向),则可以预期对垂直分量和径向分量的互相关或反褶积会在P-S延时dt处给出正波峰,因为地震能量向下穿过低速层5和基岩6之间界面时其速度一般是增大。所记录地震能量的垂直分量和横向分量之间的互相关或反褶积应检测到不是在径向-垂直平面内的波传播,例如由近地表散射引起的波传播。
在一个优选实施例中,对在一个接收器处对许多不同的爆炸得到的互相关图进行平均,以增大所得到的互相关图的信噪比,称作叠加互相关图。这会增大垂直-径向互相关图或反褶积上时移波峰振幅与垂直-横向互相关图或反褶积上该振幅的比值,如果在低速层和基岩之间的界面局部近似水平而且低速层近似于方位各向同性的话。
本发明不限于使用反褶积或互相关技术来得到穿过低速层的走时差,而是任何适当的技术都能使用。原则上,使两道波形相关的任何算法都能使用。另一种适当技术的一个实例是双相干时延估计方法,如L.Ikelle在“地球物理学”第62卷1947页(1997)描述的那样,使用质点速度的径向分量和垂直分量(或者如下文讨论的那样,使用压强和质点速度径向分量)。
图2显示上行P波发生的波型转换。如图11中所示,如果上行波是S波,会发生类似的效应。当上行S波8穿过基岩6和低速层5之间的界面时,它将作为S波8′部分地透射并部分转换为P波8″。然后转换P波8″和透射S波8′都将穿过低速层5传播到接收器4。因此,这一情况与图2所示情况的唯一差别是现在S波到达的振幅将大于P波到达的振辐。P波到达和相应的S波到达之间的时间延迟将仍然是图2和3的举例中的dt。因此,在地震数据的垂直和水平分量之间的互相关或反褶积应再次表现出在P-S时移处的波峰(即在P波穿过低速层5的走时和S波穿过低速层5的走时差值处)。
上行横波(shear wave)通常将向着界面法线折射。然而,考虑到S波型在从下面进入低速层时其速度通常减小,所以S-P转换波一般将向离开法线的方向折射。当选择应用反褶积或互相关技术的炮检距范围时,对这一点应予以考虑。
本发明的一个应用是从使用海底放置的接收器记录的海洋地震数据中提取P波静态时移和S波静态时移之差。在这种情况下,P波静态时移通常很小,以致它能被忽略,因此,利用本发明的方法从数据中提取出的P-S静态时移近似等于S波静态时移。
为说明本发明的方法,使用B.L.H.Kennet在“分层介质中的地震波传播”(剑桥大学出版社,英国剑桥,1983)中公开的那类反射率模型,对一个1D(一维)模型产生了合成地震数据。所用模型的详细内容示于表1,从中可以看到该模型包含一个100m厚的低速层放在一层水和一层基岩之间。水层假定为有无限深度,而基岩层也假定为有无限深度。在地球深度的一个爆炸型地震能量点源(代表来自更深层的反射)用作地震能量源,它位于低速层和基岩层之间界面以下200m。一个直线形接收器阵列放置在海底。
表1
层 | P速度[m/s] | S速度[m/s] | 密度[Kg/m3] | 厚度[m] |
水 | 1500 | 0 | 1000 | 无限 |
低速层 | 1600 | 500 | 1400 | 100 |
基岩 | 2500 | 1500 | 1700 | 无限 |
源深度:低速层/基岩界面下200m
源类型:爆炸型点源,只产生P波,100Hz瑞克子波
接收器位置:在水/低速层界面,炮检距0-100m,间距5m
由这一模型得到的合成数据示于图5(a)和5(b)。图5(a)显示质点速度的垂直分量,图5(b)显示质点速度的径向分量。
在该模型中假定震源产生P波。由接收器记录的信号含有透射P分量和在基岩与低速层之间界面处发生的波型转换所产生的S分量。如图5(a)中所示,主要是在记录地震数据的垂直分量上可以见到透射P波到达。如图5(b)中所示,在数据的径向分量上可以见到透射P波能量和P-S转换波的到达,特别是在增大炮检距的情况下更是如此。
对于典型的海底地震反射测量,大部分透射P能量将被记录在接收器的垂直分量上。通常,在地震能量从下面进入低速层时其速度急剧降低,所以透射波的射线路径将被强烈地折射向垂直线。与此相反,由基层和低速层之间界面处的转换产生的S波将主要出现在接收器记录的地震能量径向分量上。
假定基岩6和低速层5之间的界面是水平的,从而在界面处P波型向S波型转换过程中没有Sh波形产生。记录数据的切向分量将为零,转换后的S波将是纯Sv波。于是,依照C.A.Langston在“地球物理研究杂志”第84卷4749-47462页(1979)中所述,低速层的理论位移响应能写作:
DV(t)=S(t)*EV(t)
DR(t)=S(t)*ER(t) (1)
在式(1)中,S(t)是入射P波的震源时间函数,DV(t)和DR(t)是数据的垂直分量和径向分量,EV(t)和ER(t)分别是垂直和径向脉冲响应。通过从径向分量反褶积记录数据的垂直分量,便能从数据中准确地提取出透射P波和转换S波之间的走时差。这一处理的基本思想在于所记录数据的垂直分量能被看作震源子波的好的估计,因为对于数据垂直分量的低速层脉冲响应,对于到来的平面P波而言,接近于δ函数(上述Langston文章)。换言之:
P波混响也将在垂直分量上,但它们的振辐衰减很快。所以,这一反褶积结果将包含低速层对于横波能量的传递函数(在地震地震学中称作接收器函数):
在时间域中的结果能被解释为很像一个反射率序列。
因为反褶积基本上是频率域中的除法,它能变得不稳定。所以,优选地,应防止分母取太小的值。这是通过把分母的频率孔洞填充成其最大值的一个分数c来实现的,这是由S.K.Dey-Sarkar和R.A.Wiggins在“地球物理研究杂志”第8卷3633-3641页(1976)中提议的一项技术。能发生的另一个问题是噪声能在反褶积结果中导致没有物理意义的并且在数据中不是清楚存在的高频。以一个高斯函数在频率域中乘该结果能避免这一问题,该高斯函数以零为中心并使用高斯宽度参数a来控制其宽度(从而控制频率成分)。于是下式给出频率域中的稳定的反褶积(上述Langston文章):
这里
这里DR(ω)和DV(ω)分别是数据的径向分量和垂直分量的付立叶变换,*代表复共轭。ER(ω)是反褶积后的频率域径向地球响应,能被直接变换回时间域以给出所希望的接收器函数。φSS(ω)和G(ω)分别是所描述的稳定化的分母和高斯“滤波器”。
由于非垂直入射,在数据的径向分量上将有一些P波能量(图5(b))。这在反褶积结果中也将存在,在那里它应正好在零时刻。因为对反褶积的稳定化,接收器函数脉冲的宽度变为有限,所以它的一部分会绕着接收器函数迹线的端部卷起(图6(a))。在一个优选实施例中,这通过在径向分量的前面以迹线长度n补零并在垂直分量的后面补同样数量的零来防止。于是,反褶积后的迹线(长度2n)将在样品n+1处有其零时,而且将没有反卷的到达(图6(b))。
这样,理想上该接收器函数能被解释为低速层对从下面入射的P波的横波响应。所以,该接收器函数的振幅包含关于低速层介质性质(即速度和密度)以低速层/基岩速度反差的信息。然而,通过填充频率孔洞来稳定反褶积改变了接收器函数的振幅。依照C.J.Ammon在“美国地震学会会报”第81卷2504-2510(1991)中的建议,能恢复接收器函数的真正振辐。稳定反褶积的影响能通过对垂直分量的反褶积从它本身估计出来(使用相同的参数a和c值),已知这一反褶积迹线应在零时刻有一单个波峰,其振幅为1。所以,通过以这一自反褶积的最大值除接收器函数,真振辐能得到恢复。
图7显示由表1模型参数计算的透射P波和转换S波理论走时差与从图5(a)和5(b)的合成接收器迹线中取出的测量时差的比较。在图7中,理论走时差以虚线表示,而从图5(a)和5(b)的合成数据中提取的测量时差以实线表示。将会注意到实线有0.2ms的台阶,这是由于在产生合成地震数据时使用的采样时间间隔为0.2ms。可以看到,对图5(a)和5(b)的合成数据进行反褶积得到的时间延迟与由模型参数计算出的理论走时差符合得很好。
在表1的模型中,数据是由位于低速层5的下边界之下地球内部的震源产生的。然而,上文联系表1模型描述的技术能以完全相同的方式应用于由表面震源发射的并且已在地球内部由较深层折射或反射因而上行穿过本方法所用地球内部模型部分的地震能量所得到的数据。
如上所述,使用质点速度的垂直分量和径向分量作为指示地震数据的两个参数的好处在于P波型事件和S波型事件被自然地分开,以P波型事件主要出现在垂直分量上,而S波(横波)型事件主要出现在径向分量上。然而,本发明不限于使用质点速度的垂直和径向分量作为指示地震数据的两个参数。
在对图5(a)和5(b)的合成数据进行反褶积时,反褶积是对所有炮检距进行的,使用合成道的全时间长度。这是可能的,因为计算合成数据用的地球模型极其简单而且获取数据的几何布局也很简单(一种透射震源-接收器,以震源位于模型的基岩区)。然而,这种简单的直接了当的解决方案通常不能用于实践。在实际测量中,地震数据是使用更复杂的获取数据几何布局,而且在全时间长度记录道上对全部炮检距进行数据反褶积通常需要相当大的处理能力。所以,优选的作法是对选定的时间范围和/或炮检距范围进行反褶积或互相关以减少所需要的处理。还优选的作法是选择时间范围和/或炮检距范围以减小反褶积或互相关结果中的噪声。
如在上文中对反褶积方法理论背景的解释中提到的那样,当压缩波从下方入射到低速层时,反褶积结果基本上是低速层的横波脉冲响应。所以,反褶积和互相关能被看作是确定代表低速层的滤波器的途径,在这种情况中滤波器的相位延迟部分有特别的意义。以这种方式去做,显然在选择的数据窗内不应存在下行能量,如水层混响,以及从震源到接收器的直达波。再有,在选定的数据中不应存在在基岩中转换成横波能量的波(除非如下文中参考图11讨论的那样,把在低速层界面S到P的转换用于寻找P静态时移和S静态时移之差)。以这里考虑的方式,这些类型的地震能量未被低速层滤掉,所以能被这一方法看作噪声。优选地,选择用于反褶积/互相关过程的数据,以尽可能地去掉这一能量。
优选地,用于为反褶积/互相关过程选择数据的时间窗长度要短,以使所需要的处理最少。尽管在理论上若干个相继的上行反射P波事件能被低速层以完全相同的方式滤掉(部分地转换为横波),但在谱估计领域长期以来建立的认识是,取较长的数据窗并不会减小结果的方差,尽管两个分量之间的滤波器关系相同。此外,较长的时间窗潜在地包含某些更多的不希望的到达,导致互相关或反褶积结果中附加的噪声。选择数据径向分量所用时间窗应覆盖其时间延迟对应于最大预计时间延迟的横波到达。
优选地,在选择反褶积/互相关所用的那部分数据时,还对炮检距选择数据。对于平面分层介质,零炮检距数据将不包含任何在低速层边界的P波到S波转换所产生的横波能量,因为对于垂直入射的压缩波将不会发生任何转换。另一方面,具有长炮检距的数据(即炮检距显著地大于反射体深度的数据)将不满足地表一致性静校正假设,静态校正将变为动态校正。所以,优选的作法是从反褶积/互相关数据中排除零炮检距和长炮检距的数据。
对于其低速层和基岩的结构先前已有了解的测量场地,有可能计算出预期发生的最大P-S转换所在炮检距范围。如果是这样,则为反褶积/互相关选择的炮检距范围应该包括预期发生最大P-S转换的炮检距范围。
在另一实施例中,(定标的)压强记录用于代替质点速度的垂直分量。X.Li等在“Lomond数据分析:检波器耦合器与转换波成像”(Edinburgh各向异性计划(应用地震各向异性)研究报告,第7卷,转换波II:案例,185-212页(1999/2000)(中已经观测到,由于使用海底获取电缆得到的质点速度两个分量之间存在的不希望的互耦合,从质点速度径向分量向质点速度垂直分量的少量(<5%)横波泄漏能够发生。这一不希望的互耦合能伴有6ms量级的小相移,所以将作为静态事件出现在反褶积和互相关结果中,并带有6ms量级的走时差“。在压强和质点速度的径向分量之间不存在这种耦合效应,所以由质点速度的径向分量反褶积或互相关压强来计算的接收器函数不会受到这一相移。然而,由于压强记录的标量性质,在反褶积或互相关结果中零时或投影部分将更加显著。
在上述实施例中,已通过对指示所测量的地震数据的参数进行反褶积或互相关来确定时移dt。在波型转换事件有低振幅的情况下,这类技术是有利的。如果低速层5和基岩6之间的性质反差比本发明的方法能够工作所需要的反差还要更加强烈,则P-S转换事件能有如此高的振幅,以致不必要进行反褶积或互相关(即明显地计算接收器函数)来估计走时差。代替的是,可以在地震数据中直接观测到P-S转换震相。
图8显示一个地震测量布局,其中由多个震源-接收器对获取数据,这些震源-接收器对都有相同的炮检距,即对每个震源-接收器对,炮检距的大小和符号都相同。这称作(有符号的)共炮检距剖面采集几何布局。震源发射P波型地震能量,它被反射体3反射回到接收器4。当反射P波向上穿过基岩6和低速层5之间的界面时,波型转换发生,于是在接收器接收的地震能量包含P波型和S波型。震源-接收器对的炮检距选择为对应于预期发生最大P-S转换的范围(该炮检距不是零炮检距,也不是太大的炮检距)。
图9(a)和9(b)分别显示在图8的测量布局中由接收器4测量的质点速度的垂直分量和径向分量。可以看到,反射P事件主要出现在质点速度的垂直分量中,而转换后的S事件主要出现在质点速度的径向分量中。发生了一些轻微的泄漏,所以在质点速度的径向分量中有对应于反射P事件的弱事件,在垂直分量中有对应于转换后的S事件的弱事件。
在质点速度的径向分量中的P-S转换事件有足够大的振幅直接可以看到。这样,在垂直分量中已记录到强的P事件之后短的采样时间处(在0-100ms之后)仔细检查质点速度的径向分量,便可能直接观测到相应的P-S转换事件。为提取P波和相应的转换S波的走时差,可在不同接收器产生的各道之间追踪P事件。它的到时可用作零时基准确定在径向分量中预期出现相应S事件的时间窗起点。时间窗的长度对应于最大的预期横波静态时移。一旦在该时间窗内确定了相应的S事件位置,便能直接从数据中读出P事件和S事件之间的走时差。或者,S事件也能在各道之间追踪到。
在图9(a)和9(b)中所示各道中P事件不是“平直”的,因为它不是在每道中同时发生。图10(a)和10(b)显示图9(a)和9(b)的数据在经过初步处理使P事件“平直”之后的数据。
上面描述的本发明实施例涉及向上传播波场在穿过基岩6和低速层5之间界面时受到部分波型转换时发生的波型转换。然而,已向下传播的波在低速层5和基岩6之间的界面受到反射时也可能发生波型转换。图12(a)和12(b)示意性显示这一情况。
在图12(a)中,P波7向下传播穿过低速层5,并在低速层5和基岩6之间的界面被反射。在反射时发生部分方式转换,这样,除了上行的反射P波7′外,还存在上行波型转换后的S波7″。上文描述的用于确定P波静态时移和S波静态时移之差的技术能同样地应用于图12(a)所示测量布局。该方法能以应用于上行透射(和部分转换)波场时完全相同的方式应用,只是在考虑下行震源波场的反射时得使用在较短炮检距记录的数据。这是因为在给定的炮检距在透射布局中使用的模型中由深部的反射体反射的波场与以同样炮检距到达的在低速层5和基岩6之间界面处反射的波相比,其传播更靠近垂直线。
图12(b)总体上与图12(a)对应,但显示在下行S波场8反射时发生的部分波型转换,给出了上行S波场8′和上行波型转换后的P波场8″。
图13(a)和13(b)显示本发明的又一些实施例。在这些实施例中由向下传播穿过低速层5并在低速层5和基岩6之间的界面处受到临界折射的波场诱发波型转换。
图13(a)显示下行传播穿过低速层5的P波9。角θi表示入射角,即波传播方向和界面法线(或者,在非平面界面的情况下界面的瞬时法线)之间的夹角。当P波9入射在低速层5和基岩间的界面上时,人们通常会预期P波9会部分地透射到基岩6中,如虚线9′所示,以及部分地在界面被反射。透射波9′与界面法线之间的夹角θr由斯涅尔(Snell)定律与入射角θi关联。然而,在基岩6中地震能量的速度通常大于低速层5中地震能量的速度,于是,如果入射角超过某个临界角,斯涅尔定律将预测出sinθr的值大于1。当发生这种情况时,折射波9″将沿着低速层5和基岩6之间界面传播。这种现象称作“临界折射”;类似于光学中的完全内反射现象。临界折射波能在低速层中激发向上传播的波,称作“首波”,这些波被接收器4记录。每当震源-接收器炮检距足够大以致在相关界面上的入射角θi大于出现临界折射的临界角时,便能在地震数据中发现这种首波。在低速层5和基岩6的情况中,临界炮检距小,因此对于大多数接收器位置和大多数爆点,都记录到首波。
作为一个例子,考虑一个低速层5由P波传播速度为α1和S波传播速度为β1的物质构成,该低速层5上覆在P波传播速度为α2和S波传播速度为β2的基岩6上,这里满足下列关系:
α1<α2
β1<β2;以及
β2<α1。
图13(a)显示入射角θi满足sinθi=α1/α2关系的情况。这个入射角称作第一临界角,将称作θcαα。
当下行P波9的入射角满足θi=θcαα时sinθr取值1,意味着θr=90°。所以将不存在传播到基岩6中的透射波9′,代之以出现临界折射P波9″,它沿低速层5和基岩6之间的界面传播。由于这一折射的P波沿界面传播,它将在低速层5内激发P波10,而且这些P波10将产生S波11,这是低速层5和基岩6之间界面处波型转换的结果。图13(a)中示出两对P波10和波型转换后的S波11,尽管沿折射波9″路径的全部点都将发生P和S波10、11的激发。激发出的P波10将以与界面法线夹角θcαα传播,而波型转换后的S波11将以与法线夹角θcβα传播,这里sin(θcβα)=β1/α2。对于等于或大于θcαα的所有入射角都发生临界折射P波。图13(b)总体上对应于图13(a),但显示的是下行传播穿过低速层5的S波12的临界折射。如果入射角θi等于或大于第二临界角θcβα,则下行S波12将部分转换为临界折射P波12″,它在低速层5和基岩6之间的界面传播。临界折射P波又激发P波13和波型转换S波14,它们上行传播穿过低速层5。
在图13(a)和13(b)所示情况中,临界折射P波9″和12″都发射P波能量和S波能量上行穿过低速层5。在接收器4检测P波和S波能量再次产生P波静态时移和S波静态时移之差。S波11是由低速层5和基岩6之间界面处的波型转换形成的,因此,在界面上一点处激发的P波10与相应的波型转换S波11二者穿过低速层的走时差将等于P波型静态时移和S波型静态时移之差。这样,在接收器4的P波10的到时(P波至)和S波11的到时(S波至)的时间差等于P波静态时移和S波静态时移之差。这样,通过检测P波至和相应的S波至之间的时间延迟,便可能确定P和S静态时移之差。该时间延迟能由上述任何方法确定,例如通过对所测量的质点速度的垂直分量和径向分量进行反褶积或互相关,或通过对测量的压强和质点速度径向分量进行反褶积或互相关。
应该指出,存在第三个临界角θcββ,这里sinθcββ=β1/β2。当下行S波以等于或大于第三临界角θcββ的角度入射到低速层5和基层6的界面上时将产生临界折射S波,它沿着低速层5和基岩6之间的界面传播。然而,这个临界折射S波将不在低速层中产生向上传播的P波。
在图13(a)和13(b)的例子中,已经假定P波速度和S波速度在向下穿过低速层和基岩界面时都增大。有时发现不是所有的速度向下穿过界面时都增大,这导致比图13(a)和13(b)所示情况更复杂的情况。然而,每当在低速层5和基岩6之间的界面处发生临界折射,导致激出进入低速层5的P波和波型转换的S波时,应用本发明的方法就是可能的。在接收器4采集的地震能量能由上述任何方法分析以提取发射的P波能量和S波能量的走时差,从而得到P波静态时移和S波静态时移之差。
在本发明的上述举例中,P波型和S波形通过低速层的走时差是从质点速度的两个分量得到的,或是从质点速度的一个分量和压强得到的。在本发明的另一实施例中使用所测量的质点位移而不是质点速度。这样,在另一个实施例中,P波型和S波型穿过低速层的走时差是从质点位移的两个分量得到的,例如通过对质点位移的垂直分量和径向分量进行互相关或反褶积。在又一个实施例中,P波型和S波型穿过低速层的走时差是从质点位移的一个分量和压强得到的,例如通过对质点位移的一个分量和压强进行互相关或反褶积。这些实施例总体上对应于上述针对使用质点速度得到走时差的实施例,这里不再进一步描述。
在本发明的又一个实施例中使用质点加速度而不是质点速度。这样,在这些实施例中,P波型和S波型穿过低速层的走时差是从质点加速度的两个分量得到的,或是从质点加速度的一个分量和压强得到的。例如,这能通过对质点加速度的垂直分量和径向分量进行互相关或反褶积或者对质点加速度的一个分量和压强进行互相关或反褶积来完成。这些实施例总体上对应于上述针对使用质点速度得到走时差的实施例,这里不再进一步描述。原则上能够使用质点位移的更高阶导数,从而能从质点位移的更高阶导数的两个分量或是从质点位移的更高阶导数的一个分量和压强得到P波型和S波型穿过低速层的走时差。
在上述实施例中,接收器4被放置在海底上。然而,原则上本发明不限于此,接收器能放置在发生部分波型转换之后地震能量路径中的任何地方。这样,在上述实施例中,原则上接收器能放在低速层和基岩之间界面以上的任何地方。例如,在4-D时间推移地震测量中,接收器能被埋在海底内。然而,应该指出,如果接收器不是位于海底上,则需要附加处理以确定这两种波型的静态时移之差。
在上文对本发明实施例的描述中,把X方向定义在径向方向,即震源-接收器方向在海底上的投影。一个地震接收器将记录质点速度或质点位移沿两个正交的水平方向的分量,这两个方向可以称作接收器的X轴和Y轴。应该指出,接收器的布设可能没有使其X轴与径向方向对齐,因此可能需要对接收器测量的“原始”水平分量进行附加的投影或转动以从接收器的输出计算径向分量和横向分量。原则上,接收器测量到的原始X分量和Z分量能用于以上述任何方法估计P波型和S波型穿过低速层的走时差,而且如果接收器的X轴与径向方向之间的夹角不是太大的话,这将给出可以接受的结果。
本发明的处理方法可以应用于处理前已存在的地震数据。本发明还能被纳入到这样一个地震测量方法中,在那里所采集的地震数据包括由低速层和基岩之间界面处的部分波型转换所产生的波型转换事件。通过把地震能量指向该界面从而当地震能量穿过该界面或被界面反射时便能得到这类地震数据。
如上文中指出的那样,本发明能应用于在基岩和低速层之间界面处部分波型转换为P波的上传S波,以及在界面处部分波型转换为S皮的上传P波。这两种情况都能用于找出S波和P波穿过低速层的走时差。原则上,P波和S波穿过低速层的走时差依赖于入射波的慢度,即使对于海底和低速层与基岩之间界面均平坦的理想情况也是如此。然而,在实践中,对于用于典型地震测量的震源-接收器炮检距范围和对于低速层的典型物质,发现走时差对入射波慢度的依赖性很弱,而且走时差在效果上为常数,不管入射波的慢度多小。还发现典型的变化在整个接收器范围只有几毫秒。
所以,穿过低速层的走时差不依赖于入射波慢度这一假设是安全的。这一假设等效于地表一致性静校正假设,后者一般认为对大多数静态问题有效。做出这一假设意味着在接收器4记录的波场不需要滤波或分离,即使它是从地球内部的不同反射体以不同的慢度到达的波的混合物,而且还是P波和S波的混合物。再有,不必要为了应用本发明的方法而从数据中提取单个事件。这意味着在反褶积或互相关处理中能使用相对较长的时间窗(秒的量级)。再有,各道能在共接收器域叠加,以改善数据的信噪比。
上文中曾提到,接收器可能不布设成以其X轴与径向方向对齐,在这种情况下,由接收器测量到的“原始”水平分量可能需要被投影或转动以从接收器输出计算径向分量和横向分量。应该指出,如果地震力源安排在一个爆点线内而且当爆点线平行于接收器线直接放置在其上方时才被起爆,那么在这种情况下这一处理过程是不必要的。在这种情况下,原始接收器数据的X分量能被看作是径向分量,而原始接收器数据的Y分量能被看作是横向分量。然而,应当指出,当炮检距从负变为正时在地震数据径向分量上的事件的数据符号发生变化(一维地球假设)。当对来自包括正、负两种炮检距的共接收器道集的反褶积或互相关结果进行叠加时必须考虑这一点,因为忽略这一点能导致信号衰减甚至完全丢失信号。在叠加数据之前必须对一个炮检距采集的数据反转其符号。
上文中,参考互相关、反褶积和双相干性作为关联两个地震数据道的方法举例,描述了本发明。然而,本发明不限于这些方法,任何用于关联两个数据道的适当技术或算法都可以使用。
本发明还可以用于P波和S波的“混响事件”。在这些事件中包含例如穿过近表表层5和基岩6之间的界面或在界面上反射时的部分波型转换。然而,在一个混响事件中,P波和S波并不是从发生波型转换的点直接到达接收器。在混响事件中,P波和S波在近表面边界面受到一次或多次反射,从而在入射到接收器上之前多次穿过近表面层5。
地震数据还可能包含来自部分波型转换的事件,但这波型转换不是发生在近表面层边界而是发生在近表面层内。例如,这可能产生于近表面层内的分层效应。
在上述实施例中,P波或S波事件之一是由部分波型转换产生的。然而,如前文说明的那样,本发明不限于此,而是可以用于任何一对对应的P和S事件,只要这两个事件的振辐和/或波形的差异主要来源于近表面层5。例如,在使用发射P波和S波两种波的震源的地震测量中,所采集的数据将包含P波在近表面层5下边界上的特定点反射产生的事件,还将包含S波在近表面层5下边界上的那个点反射产生的相应事件。本发明可以应用于这些事件。(应该指出,由于水体不支持S波传播,所以发射S和P波型的震源需要放置在海底上、地球表面上或在井孔中)。
图14是根据本发明的可编程设备15的示意性框图。该设备包含一个可编程数据处理器16并带有程序存储器17,例如只读存储器ROM形式的存储器,用于存储程序以控制数据处理器16进行上述任何处理方法。该设备进一步包含非易失读/写存储器18用于存储例如在没有电源时必须保持的任何数据。由随机存取存储器(RAM)19为数据处理器提供“工作”存储器或“便笺式”存储器。所提供的输入接口20用于例如接收命令和数据。所提供的输出接口用于例如显示与方法的进展和结果有关的信息。用于处理的地震数据可通过输入接口20提供,或者可以从机器可读数据存储器22中提取。
用于操作该系统和实现前述方法的程序存储在程序存储器17中,它可以具体实现为半导体存储器,例如有名的ROM型存储器。然而,该程序可以存储在其他任何适当的存储介质中,如磁性数据载体17a(如“软盘”)或CD-ROM 17b。
Claims (33)
1.一种处理包括第一和第二波型地震能量的地震数据的处理方法,该方法包含处理地震数据以得到以第一波型传播的地震能量和以第二波型传播的地震能量穿过一个地层的走时差的步骤,所述第一波型地震能量和第二波型地震能量沿着实质上相同的路径穿过该地层。
2.如权利要求1所述的方法,其中第二波型是由第一波型在海底一层的边界面的部分波型转换产生的。
3.如权利要求2所述的方法,其中边界面是该层的下边界面。
4.如权利要求2或3所述的方法,其中该层是表面层或近表面层。
5.如权利要求1所述的方法,其中处理地震数据的步骤包含:确定指示地震数据第一和第二方面的第一和第二参数;以及处理该第一和第二参数以得到以第一波型传播的地震能量和以第二波型传播的地震能量穿过该层的走时差。
6.如权利要求5所述的方法,其中该处理步骤包含使第一和第二参数彼此互相关,从而得到以第一波型传播的地震能量和以第二波型传播的地震能量穿过该层的走时差。
7.如权利要求5所述的方法,其中该处理步骤包含使第一和第二参数彼此反褶积,从而得到以第一波型传播的地震能量和以第二波型传播的地震能量穿过该层的走时差。
8.如权利要求5所述的方法,其中该处理步骤包含使用双相干时间延迟方法处理第一和第二参数,从而得到以第一波型传播的地震能量和以第二波型传播的地震能量穿过该层的走时差。
9.如权利要求5所述的方法,其中该处理步骤包含由第一参数确定第一波型穿过该层的走时和由第二参数确定第二波型穿过该层的走时。
10.如权利要求5、6、7、8和9中任一项所述的方法,其中该第一参数是质点速度的第一分量,该第二参数是质点速度的第二分量,质点速度的第一和第二分量彼此不平行。
11.如权利要求10所述的方法,其中质点速度的第一分量是垂直分量,质点速度的第二分量是水平分量。
12.如权利要求5、6、7和8中任一项所述的方法,其中该第一参数是质点速度的第一分量,第二参数是压强。
13.如权利要求5、6、7、8和9中任一项所述的方法,其中该第一参数是质点位移的第一分量,该第二参数是质点位移的第二分量,质点位移的第一和第二分量彼此不平行。
14.如权利要求13所述的方法,其中质点位移的第一分量实质上是垂直分量,质点位移的第二分量实质上是水平分量。
15.如权利要求2所述的方法,其中该第一波型是P波波型,由部分波型转换产生的第二波型是S波波型。
16.如权利要求2、3、5、6、7、8、9中任一项所述的方法,其中该第一波型是S波波型,由部分波型转换产生的第二波型是P波波型。
17.如权利要求5所述的方法,其中确定第一和第二参数的步骤包含把所接收的地震数据分解为P分量和S分量。
18.如权利要求17所述的方法,其中该处理步骤包含使地震能量的P分量和S分量彼此互相关,从而得到以P波型传播的地震能量和以S波型传播的地震能量穿过该层的走时差。
19.一种地震测量方法,包含如下步骤:
使以第一波型传播的地震能量指向海底一层的边界面,使得在该边界面发生地震能量转变为第二波型的部分波型转换;
在一个或多个接收器采集包括地震能量第一和第二波型的地震数据;以及
处理所述地震数据以得到以第一波型传播的地震能量和以第二波型传播的地震能量穿过一个地层的走时差,所述第一波型地震能量和第二波型地震能量沿着实质上相同的路径穿过该地层。
20.如权利要求19所述的地震测量方法,其中采集地震数据的步骤包含在放置于该层上方的一个或多个接收器处采集地震数据。
21.如权利要求19或20所述的方法,其中采集地震数据的步骤包含在放置于海底上的一个或多个接收器处采集地震数据。
22.如权利要求19所述的地震测量方法,还包含使地震能量向下指向该层,从而当向下穿过该层的第一波型地震能量在边界面反射时发生部分波型转换。
23.如权利要求19所述的地震测量方法,还包含使地震能量向上指向该层,从而当向上传播的第一波型地震能量在透射穿过边界面时发生部分波型转换。
24.如权利要求19所述的地震测量方法,还包含使地震能量以这样一个角度指向该层边界面,该角度使第一波型地震能量沿边界面传播。
25.如权利要求24所述的地震测量方法,其中该第一波型是P波波型,入射地震能量传播方向和边界面法线方向之间的夹角θi满足sinθi≥α1/α2,其中α1是P波在海底该层中的传播速度,α2是在下伏层中P波传播速度。
26.如权利要求24所述的地震测量方法,其中该第一波型是P波波型,入射地震能量传播方向和边界面法线方向之间的夹角θi满足sinθi≥β1/α2,其中β1是S波在海底该层中的传播速度,α2是在下伏层中P波传播速度。
27.一种地震数据处理设备,该地震数据包括第一和第二波型地震能量,该第二波型是由第一波型在海底一层的边界面部分波转换产生的,该设备包含处理地震数据装置,以得到以第一波型传播的地震能量和以第二波型传播的地震能量穿过该层的走时差,所述第一波型地震能量和第二波型地震能量沿着实质上相同的路径穿过该层。
28.如权利要求27所述的设备,还包含:确定指定地震数据第一和第二方面的第一和第二参数的装置;以及处理该第一和第二参数以得到以第一波型传播的地震能量和以第二波型传播的地震能量穿过该层的走时差的装置。
29.如权利要求28所述的设备,其中处理装置适合于使第一和第二参数彼此互相关,从而得到以第一波型传播的地震能量和以第二波型传播的地震能量穿过该层的走时差。
30.如权利要求28所述的设备,其中处理装置适合于使第一和第二参数彼此反褶积,从而得到以第一波型传播的地震能量和以第二波型传播的地震能量穿过该层的走时差。
31.如权利要求28所述的设备,其中处理装置适合于使用双相干时间延迟方法处理第一和第二参数,从而得到以第一波型传播的地震能量和以第二波型传播的地震能量穿过该层的走时差。
32.如权利要求28所述的设备,其中处理装置适合于由第一参数确定第一波型穿过该层的走时和由第二参数确定第二波型穿过该层的走时。
33.如权利要求27至32任何一个所述的设备,包含一个可编程数据处理器。
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GB0101936A GB0101936D0 (en) | 2001-01-25 | 2001-01-25 | A method of procesing marine seismic data and a method of seismic surveying |
GB0101936.3 | 2001-01-25 | ||
GB0110573.3 | 2001-04-30 | ||
GB0110573A GB0110573D0 (en) | 2001-01-25 | 2001-04-30 | A method of processing marine seismic data and a method of seismic surveying |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1493005A CN1493005A (zh) | 2004-04-28 |
CN1271420C true CN1271420C (zh) | 2006-08-23 |
Family
ID=26245616
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB028052773A Expired - Fee Related CN1271420C (zh) | 2001-01-25 | 2002-01-24 | 海洋地震数据处理方法和地震测量方法 |
Country Status (8)
Country | Link |
---|---|
US (1) | US7164619B2 (zh) |
EP (1) | EP1354221B1 (zh) |
CN (1) | CN1271420C (zh) |
AU (1) | AU2002226539B2 (zh) |
CA (1) | CA2435643A1 (zh) |
EA (1) | EA005475B1 (zh) |
NO (1) | NO334087B1 (zh) |
WO (1) | WO2002059647A1 (zh) |
Families Citing this family (43)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2384068B (en) * | 2002-01-11 | 2005-04-13 | Westerngeco Ltd | A method of and apparatus for processing seismic data |
GB2384557B (en) | 2002-01-25 | 2005-06-29 | Westerngeco Ltd | A method of and apparatus for processing seismic data |
GB2394050B (en) * | 2002-10-07 | 2005-11-23 | Westerngeco Seismic Holdings | Processing seismic data |
WO2004031805A1 (en) * | 2002-10-07 | 2004-04-15 | Westerngeco Seismic Holdings Limited | Processing seismic data |
GB2397378B (en) * | 2003-01-15 | 2005-03-02 | Westerngeco Ltd | Method for retrieving local near-surface material information |
GB2398124B (en) * | 2003-02-08 | 2006-10-25 | Abb Offshore Systems Ltd | Estimating the time of arrival of a seismic wave |
WO2004086094A1 (en) | 2003-03-26 | 2004-10-07 | Westergeco Seismic Holdings Limited | Processing seismic data representative of the acceleration wavefield |
GB0306920D0 (en) * | 2003-03-26 | 2003-04-30 | Westerngeco Seismic Holdings | Processing seismic data |
GB0319201D0 (en) * | 2003-08-15 | 2003-09-17 | Westerngeco Ltd | Multi-component seismic data processing |
US7796468B2 (en) * | 2004-02-26 | 2010-09-14 | Saudi Arabian Oil Company | Prediction of shallow drilling hazards using seismic refraction data |
CN100383557C (zh) * | 2004-06-25 | 2008-04-23 | 大庆油田有限责任公司 | 一种提高地震分辨率的方法 |
US20060106546A1 (en) * | 2004-11-17 | 2006-05-18 | Time Domain Corporation | System and method for evaluating materials using ultra wideband signals |
US20080232193A1 (en) * | 2007-03-20 | 2008-09-25 | Geocyber Solutions, Inc. | Methods for Noise Removal and/or Attenuation from Seismic Data by Wavelet Selection |
US8976624B2 (en) * | 2006-05-07 | 2015-03-10 | Geocyber Solutions, Inc. | System and method for processing seismic data for interpretation |
US20070258323A1 (en) * | 2006-05-07 | 2007-11-08 | Ping An | Method for Seismic Trace Decomposition and Reconstruction Using Multiple Wavelets |
US7773456B2 (en) * | 2006-10-02 | 2010-08-10 | Bp Corporation North America Inc. | System and method for seismic data acquisition |
EP2116870A1 (en) * | 2008-05-05 | 2009-11-11 | Shell Internationale Research Maatschappij B.V. | Method for monitoring a multi-layered system |
CA2740883C (en) * | 2008-10-16 | 2015-02-10 | Schlumberger Canada Limited | Vsp imaging of acoustic interfaces |
CA2747426C (en) | 2009-01-09 | 2017-05-23 | Exxonmobil Upstream Research Company | Hydrocarbon detection with passive seismic data |
US20100312480A1 (en) * | 2009-04-24 | 2010-12-09 | Hansteen Fredrik | Method for monitoring fluid flow in a multi-layered system |
US9025413B2 (en) * | 2009-12-07 | 2015-05-05 | Pgs Geophysical As | Method for full-bandwidth source deghosting of marine seismic streamer data |
US20130114375A1 (en) * | 2010-07-26 | 2013-05-09 | Mark A. Meier | Seismic Acquisition Method for Mode Separation |
US8243548B2 (en) * | 2010-08-27 | 2012-08-14 | Board Of Regents Of The University Of Texas System | Extracting SV shear data from P-wave seismic data |
US8325559B2 (en) | 2010-08-27 | 2012-12-04 | Board Of Regents Of The University Of Texas System | Extracting SV shear data from P-wave marine data |
CN102478665B (zh) * | 2010-11-23 | 2013-07-31 | 中国石油天然气集团公司 | 一种确定地震波入射角和振幅的方法 |
CN102590856A (zh) * | 2011-01-11 | 2012-07-18 | 中国科学院地质与地球物理研究所 | 基于小波频谱分析的位场异常分离方法 |
US20130003499A1 (en) * | 2011-06-28 | 2013-01-03 | King Abdulaziz City For Science And Technology | Interferometric method of enhancing passive seismic events |
US9116256B2 (en) * | 2011-07-18 | 2015-08-25 | Cggveritas Services (U.S.) Inc | Method and device for wave fields separation in seismic data |
AP2014007960A0 (en) * | 2012-03-06 | 2014-09-30 | Univ Texas | Extracting SV shear data from P-wave marine data |
US9476999B2 (en) | 2012-06-04 | 2016-10-25 | Shell Oil Company | Seismic imaging |
US10125599B2 (en) * | 2012-08-02 | 2018-11-13 | Micross Advanced Interconnect Technology Llc | Location of sensors in well formations |
US9429670B2 (en) * | 2012-09-28 | 2016-08-30 | Cgg Services Sa | Plural-depth buried seismic sensors acquisition system and method |
US9846247B2 (en) * | 2013-03-15 | 2017-12-19 | Pgs Geophysical As | Methods and systems for interpolation of multi-component seismic data collected in a marine environment |
WO2015013543A1 (en) * | 2013-07-24 | 2015-01-29 | Digitalglobe, Inc. | Bathymetric techniques using satellite imagery |
WO2017024536A1 (zh) * | 2015-08-11 | 2017-02-16 | 深圳朝伟达科技有限公司 | 一种自动切除地震波波至的方法 |
CN105425287B (zh) * | 2015-11-09 | 2017-10-27 | 中国地质大学(北京) | 地震波的叠前分离方法 |
US10474767B2 (en) | 2016-01-26 | 2019-11-12 | Saudi Arabian Oil Company | Gravity modeling a rifted continental margin |
US10310133B2 (en) * | 2016-10-04 | 2019-06-04 | Fairfield Geotechnologies | Calibration of geophone and hydrophone pairs |
US11327188B2 (en) * | 2018-08-22 | 2022-05-10 | Saudi Arabian Oil Company | Robust arrival picking of seismic vibratory waves |
CN109270574B (zh) * | 2018-10-18 | 2019-09-03 | 国家海洋局第二海洋研究所 | 一种基于多种震源采集数据的联合反褶积方法 |
CN110244355B (zh) * | 2019-07-25 | 2021-06-08 | 西南交通大学 | 一种基于震源断层模型的脉冲地震动模拟方法 |
CN111324968B (zh) * | 2020-03-06 | 2023-03-28 | 西南大学 | 一种倾斜地层隧道工程微震监测传感器布设方法 |
CN115980856B (zh) * | 2023-03-21 | 2023-06-30 | 自然资源部第一海洋研究所 | 适用复杂海底界面的深拖式多道地震炮检精确定位方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4597066A (en) * | 1983-04-20 | 1986-06-24 | Chevron Research Company | Method of seismic processing and displaying simultaneously collected conventional and converted P- or S-wave data |
US4683557A (en) * | 1984-10-05 | 1987-07-28 | Mobil Oil Corporation | Acoustic logging method for identifying subsurface formation boundaries |
US5200928A (en) * | 1991-11-07 | 1993-04-06 | Chevron Research And Technology Company | Method for using mode converted P- to S- wave data to delineate an anomalous geologic structure |
US5917160A (en) | 1994-08-31 | 1999-06-29 | Exxon Production Research Company | Single well system for mapping sources of acoustic energy |
US5610875A (en) | 1995-09-29 | 1997-03-11 | Western Atlas International, Inc. | 3-D converted shear wave rotation with layer stripping |
US6128580A (en) * | 1998-04-17 | 2000-10-03 | Bp Amoco Corporation | Converted-wave processing in many-layered anisotropic media |
US6694261B1 (en) * | 1999-06-07 | 2004-02-17 | Conoco Inc. | Method for identification of shallow water flow hazards using marine seismic data |
-
2002
- 2002-01-24 EA EA200300832A patent/EA005475B1/ru not_active IP Right Cessation
- 2002-01-24 CA CA002435643A patent/CA2435643A1/en not_active Abandoned
- 2002-01-24 WO PCT/GB2002/000280 patent/WO2002059647A1/en not_active Application Discontinuation
- 2002-01-24 US US10/466,980 patent/US7164619B2/en not_active Expired - Lifetime
- 2002-01-24 AU AU2002226539A patent/AU2002226539B2/en not_active Expired
- 2002-01-24 CN CNB028052773A patent/CN1271420C/zh not_active Expired - Fee Related
- 2002-01-24 EP EP02716147A patent/EP1354221B1/en not_active Expired - Lifetime
-
2003
- 2003-07-22 NO NO20033311A patent/NO334087B1/no not_active IP Right Cessation
Also Published As
Publication number | Publication date |
---|---|
EP1354221A1 (en) | 2003-10-22 |
US7164619B2 (en) | 2007-01-16 |
EP1354221B1 (en) | 2013-02-27 |
NO20033311D0 (no) | 2003-07-22 |
WO2002059647A1 (en) | 2002-08-01 |
CN1493005A (zh) | 2004-04-28 |
US20040076077A1 (en) | 2004-04-22 |
AU2002226539B2 (en) | 2006-03-09 |
EA005475B1 (ru) | 2005-02-24 |
NO334087B1 (no) | 2013-12-09 |
CA2435643A1 (en) | 2002-08-01 |
EA200300832A1 (ru) | 2003-12-25 |
NO20033311L (no) | 2003-09-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1271420C (zh) | 海洋地震数据处理方法和地震测量方法 | |
CN107526101B (zh) | 一种获取地震反射波的采集和处理方法 | |
US9128205B2 (en) | Process for creating image gathers | |
CN112840236B (zh) | 通过反演分离不同类型的多个震源 | |
CN1073236C (zh) | 处理具有多次反射噪声的地震数据的方法 | |
US10317554B2 (en) | Noise attenuation via thresholding in a transform domain | |
AU2002226539A1 (en) | A method of processing marine seismic data and a method of seismic surveying | |
EP0443234B1 (en) | Method for multifold vertical seismic profile data acquisition and processing | |
US20120035852A1 (en) | Processing seismic data | |
EP2548052A2 (en) | System and method of 3d salt flank vsp imaging with transmitted waves | |
MX2014010606A (es) | Extraccion de datos de cizallamiento sv a partir de datos marinos de la onda p. | |
Schultz | Time-lapse multicomponent geophone and DAS VSP processing and analysis | |
Mari et al. | 3D seismic imaging of a near-surface heterogeneous aquifer: a case study | |
CN1656388A (zh) | 处理地震数据的方法和装置 | |
US20120269035A1 (en) | Evaluating Prospects from P-Wave Seismic Data Using S-Wave Vertical Shear Profile Data | |
Gómez de la Peña | The origin and tectono-sedimentary structure of the Alboran Basin | |
Khalil et al. | Full-waveform processing and interpretation of kilohertz cross-well seismic data | |
Drijkoningen et al. | Nongeometrically converted shear waves in marine streamer data | |
Lee et al. | Fracture imaging using Image Point transform and midpoint imaging of RVSP data | |
Lau et al. | Integrating streamer and ocean-bottom seismic data for sub-basalt imaging on the Atlantic Margin | |
Digranes et al. | Use of mode-converted waves in marine seismic data to investigate the lithology of the sub-bottom sediments in Isfjorden, Svalbard | |
Yu et al. | Imaging strong lateral heterogeneities across the contiguous US using body-to-surface wave scattering | |
Aarseth | Crustal structure across the western Barents Sea from deep seismic and potential field data. Implications for Caledonian trends and subsequent basin formation | |
Mendes | Seismic tomography5 | |
Çoruh et al. | Composite refraction-reflection stack sections: imaging shallow subsurface features |
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: 20060823 Termination date: 20170124 |
|
CF01 | Termination of patent right due to non-payment of annual fee |