CN1181358C - 用于局部浪高确定的方法和在海洋地震信号中的声传感器 - Google Patents
用于局部浪高确定的方法和在海洋地震信号中的声传感器 Download PDFInfo
- Publication number
- CN1181358C CN1181358C CNB998146323A CN99814632A CN1181358C CN 1181358 C CN1181358 C CN 1181358C CN B998146323 A CNB998146323 A CN B998146323A CN 99814632 A CN99814632 A CN 99814632A CN 1181358 C CN1181358 C CN 1181358C
- Authority
- CN
- China
- Prior art keywords
- seismic
- time
- sensor
- arrival
- depth
- 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 60
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 35
- 238000001228 spectrum Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 abstract description 16
- 238000012545 processing Methods 0.000 abstract description 8
- 230000002238 attenuated effect Effects 0.000 abstract 1
- 230000008859 change Effects 0.000 description 11
- 239000012530 fluid Substances 0.000 description 7
- 230000014509 gene expression Effects 0.000 description 7
- 239000013535 sea water Substances 0.000 description 7
- 230000000916 dilatatory effect Effects 0.000 description 6
- 238000005070 sampling Methods 0.000 description 6
- 230000001052 transient effect Effects 0.000 description 6
- 238000012937 correction Methods 0.000 description 5
- 230000006872 improvement Effects 0.000 description 5
- 241001269238 Data Species 0.000 description 4
- 230000008901 benefit Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 102100027841 Acyl-CoA wax alcohol acyltransferase 2 Human genes 0.000 description 2
- DOEADYYICZVJDD-UHFFFAOYSA-N [4-[(4-aminophenyl)diazenyl]phenyl]arsonic acid Chemical compound C1=CC(N)=CC=C1N=NC1=CC=C([As](O)(O)=O)C=C1 DOEADYYICZVJDD-UHFFFAOYSA-N 0.000 description 2
- 108010024239 aromatic amino acid aminotransferase Proteins 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 238000011002 quantification Methods 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000013517 stratification Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000001351 cycling effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 239000003350 kerosene Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008450 motivation Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000001915 proofreading effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- GOLXNESZZPUPJE-UHFFFAOYSA-N spiromesifen Chemical compound CC1=CC(C)=CC(C)=C1C(C(O1)=O)=C(OC(=O)CC(C)(C)C)C11CCCC1 GOLXNESZZPUPJE-UHFFFAOYSA-N 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000725 suspension Substances 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/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
-
- 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/56—De-ghosting; Reverberation compensation
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Oceanography (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
一种联系海洋地震数据获得活动使用的分析地震信号的方法,包括步骤:使用一个水中声传感器(18)检测一个接收机幻影信号;和使用所述接收机幻影信号(30)估计声传感器上方的水柱高度。该方法可以进一步包括步骤:从由与声源不同偏移定位的多个水中声传感器(18)接收的地震信号,辨别到达时间的变化;确定在辨别到达时间变化与同假定声传感器深度分布有关的期望到达时间变化之间的时间差;及把时间差转换成在假定声传感器深度分布与所述声传感器的实际深度分布之间的深度差。该方法为局部浪高和声传感器高程的改进确定而提供,并且允许在地震数据中与局部浪高和地震传感器高程的变化有关的“噪声”,在以后数据处理期间衰减。
Description
本发明涉及一种分析地震信号的方法,更具体地说,涉及一种适于联系海洋地震数据获得活动使用的分析地震信号的方法,该方法提供局部浪高和声传感器深度的改进确定,并且允许在以后的数据处理期间减小在与局部浪高和声传感器深度的变化有关的地震数据中的噪声。
收集地震数据以远程检测地下地质状态,特别是联系对油气的勘探和其生产,如油和天然气。为了收集在海洋环境中的地震数据,诸如***之类的声源用来产生经海水传播并且进入地下地质建造的声信号。在海底处和在不同地质层之间声阻抗的变化使声能的一部分反射并且向海面返回。这些反射信号由声传感器接收,并且处理成建立地下地质的图象。
在海洋环境中,这些声传感器(也叫做地震传感器,经常是称作水听器的压力传感器)一般包含在长管形拖缆内,并且拖在地震勘测船后面。拖缆填充有允许拖缆段近似实现无浮力的煤油或其他浮力材料。拖缆经常带有提供基本抗拉强度和禁止拖缆段延伸的一个或多个内部应力件(如钢缆),同时允许拖缆较柔软并且能够缠绕在地震勘测船上的一个适当直径的鼓上。拖缆拖拉的深度(或“高程”)一般由定位在离地震勘测船最近端部处的反射器调节(例如,见美国专利No.5,357,892),并且由沿拖缆长度以正常间隔放置的叫做吊舱的控制器件调节(例如,见我们发表的PCT国际申请No.WO 98/28636)。
在拖缆中水听器的深度一般在地震勘测船上由附着到吊舱上的深度传感器监视。因为与显著靠近的水听器间隔(如每12.5米一组水听器)相比,吊舱沿拖缆隔得较宽(如每300米),所以一个具体声传感器或一组声传感器的深度一般必须通过从在传感器或传感器组侧的吊舱深度值的内插近似。
这种类型的较粗略深度确定***使得地震勘测船船员难以确定何时某些类型的问题出现在拖缆内。例如,拖缆段一般“平衡”,直到他们近似无浮力。由于在地震勘测船上和在海水中的温度变化,平衡问题(过分正或负的浮力)有时出现。然而,如果能监视在每段中水听器每一个的深度,则有可能确定哪些段正在经受平衡问题,并且在他们影响获得的地震数据或引起拖拉问题之前更正这些问题。
在吊舱上的深度传感器一般检测局部周围水压力,并且把该压力读数转换成深度值。然而,在吊舱处测得的水压力包括随拖拉拖缆恒定变化的两种类型的瞬态条件。第一种瞬态条件是局部浪高,在传感器正上方的局部海平面减去平均海平面。局部浪高的变化也称作波浪。第二种瞬态条件是相对于平均海平面测得的实际拖缆高程(或深度)。实际拖缆高程的变化一般归因于诸如在拖缆段中的正或负浮力之类的力、波浪诱导力、洋流、反射器、吊舱等。在吊舱处的水压力都受这些瞬态条件的影响。为了消除波浪影响,一般在一个延伸时间段上(如在10与100秒之间)平均或滤波测得的水压力值。尽管这种平均或滤波对于吊舱产生更准确的“平均”深度值,但它消除了使用测量深度值补偿具有小于平均时段或滤波长度一半的循环周期的瞬态条件,如波浪,的可能性。
两种类型的“噪声”由流束深度和局部浪高的波动引入到数据中。第一种类型的噪声由幻影效应引起。来自声传感器或声源上方海面的声反射将引起在与传感器或源深度有关的频率下接收声信号的抵消(即“幻影效应”)。幻影是在频谱中在频率F=n/Tg下出现的凹痕,其中n是整数(0、1、2、…),并且幻影周期Tg等于接收机(或源)深度H(至海面的距离)除以海水声传播速度的两倍。深度H(并因此幻影频率F)对于入射角需要校正(如在下面将要更详细讨论的那样)。有两种纪影,一种在源侧引入而一种在接收机侧引入。当接收机或源的深度变化时,幻影陷频的变化出现。这些变化归因于由在流束或源上方浪高的变化造成的流束或源的绝对高程的变化。
为了补偿这种幻影效应,地震传感器一般拖拉在其中第一非零凹痕幻影频率在地震频谱外侧(在约5Hz和约80Hz之间)的深度处,其中在地震勘测期间得到关于感兴趣地质地下层的绝大多数信息。一个消卷积过程能用来补偿由幻影效应引起的接收地震信号的频率独立衰减(即把数据“去幻影”)。然而,在常规地震数据处理过程中,这种消卷积过程假定地震传感器放置在海面下方一个恒定距离。传感器位置与该假定位置的任何偏离将使去幻影过程不适当地操作至某种程度;某些频率将被过放大,而某些频率将保持欠放大。由于在延伸时间段上平均或否则滤波深度值以除去波浪对深度值的影响,所以当传感器正在接收感兴趣的地震数据时,由常规地震数据获得设备提供的深度值不能用来提供描述传感器的实际(和变化)深度值的定制或个别化的去幻影地震数据。
第二种类型的噪声归因于流束绝对高程的变化,流束引起从下面地震反射器接收的地震信号的到达时间的无意移动。由于绝大多数数据分析涉及把设想相同地下位置的多个地震轨迹结合在一起,所以这些时间移动将引起反射器地震图象的模糊。
尽管这两种类型的偏移在其常规意义上没有引入“噪声”(即干扰或掩蔽希望信号的多余信号),但容易理解,他们禁止地下层适当地震成象,并且在其更一般的意义上构成噪声。对于某些类型的地震成象,如分析产生油气储层的延时图象,这些效应可以足以掩蔽能期望从储层流体的撤出引起的地震响应的任何变化。代表申请人进行的研究已经得出结论,如果利用常规的地震数据处理方案,则来自只要2米显著浪高(SWH)海的不平海效应就能掩蔽地震响应的任何变化,可能期望该地震响应至少对于某些储层类型与油气生产有关。
在常规海洋地震勘测中,补偿局部海面高度变化进行的唯一偿试涉及补偿由潮汐效应造成的平均海平面的变化。对于波浪影响或流束深度值的短循环时间变化没有进行校正地震数据的偿试。尽管众所周知,如果在海浪大时段期间获得地震数据,则将严重降低地震数据的质量,但通常没有进行偿试,以补偿这些类型的瞬态条件。当地震勘测船船员或其客户的船上观察者决定海上条件太差或者未能满足达成的合同规范时,简单地停止由地震勘测船进行的地震数据获得。简单地迫使客户与在大海浪状态期间获得的地震数据比在小海浪状态期间获得的地震数据噪声多(即具有较低的质量)的事实相处。
地震数据获得承包人具有在“有问题”天气条件下获得地震数据的显著动机,因为他们一般不补偿由坏天气导致的停工时间,而在坏天气期间花费的时间量能容易地在整个活动时段的10%与50%之间的范围内。一些地震数据获得承包人特别积极,差不多在坏天气下继续地震数据获得活动。这当地震勘测船正在获得多客户数据时特别可信。一般“投机”地获得多客户数据,地震承包人支付获得活动的成本并且然后试图通过许可对获得地震数据的访问偿还这些成本和获利。一些承包人显然相信,在以后数据处理期间能除去(或至少掩蔽坏天气的影响或者客户可能不会认识到数据实际上有多少噪声。这种情形在过去一直比较复杂,因为客户在获得地震数据时没有独立确定海浪状态的方法。
因此本发明一个目的在于,提供一种分析地震信号的改进方法。
本方法描述实施例的一个优点在于,允许联系海洋地震数据获得活动确定局部浪高和声传感器高程。
本发明描述实施例的进一步的优点在于,提供一种用来直接从地震数据确定局部浪高的客观方法。
本发明描述实施例的另一个优点在于,在没有常规基于水压力的深度传感器的情况下可以确定各个声传感器或声传感器阵列的高程。
本发明描述实施例的另外一个优点在于,在以后数据处理期间可以衰减由局部浪高变化和/地震传感器高程变化引入到地震数据中的“噪声”。
本发明涉及一种分析在海洋地震勘测期间由响应声源的操作的多个水中地震传感器获得的地震信号的方法。该方法对于至少一些信号的每一个包括步骤:在信号内选择一个时间窗口,该时间窗口得到在信号中表示的清晰事件;在所述窗口中由信号的振幅/频谱确定接收机幻影陷频;及从所述接收机幻影陷频导出在产生信号的传感器上方的水柱高度估计。该方法可以进一步包括步骤:从由与声源不同偏移定位的多个水中声传感器接收的地震信号,辨别到达时间的变化;确定在辨别到达时间变化与同假定声传感器深度分布有关的期望到达时间变化之间的时间差;及把时间差转换成在假定声传感器深度分布与所述声传感器的实际深度分布之间的深度差。该方法为局部浪高和声传感器高程的改进确定而提供,并且允许在地震数据中与局部浪高和地震传感器高程的变化有关的“噪声”,在以后数据处理期间衰减。
图1是进行海洋地震勘测的地震勘测船的示意图;
图2是地震拖缆一部分和来自图1的海面的放大图;
图3是从单个声脉冲接收的地震轨迹的图片;
图4是在海底反射已经对准近似相同到达时间之后来自图3的地震轨迹的一个窗口部分;
图5是曲线图,表示振幅与来自图4的开窗口轨迹的频谱的关系;
图6是曲线图,画出水底反射到达时间与偏离距离的关系及一条传播时间外移曲线;及
图7是曲线图,表示对于图3中所示一组地震轨迹按照本发明确定的地震传感器深度和局部海面高度。
图1表示进行海洋地震勘测的地震勘测船和有关设备的示意图。一个常规海洋地震勘测***10包括一艘地震勘测船12、一个声源14、一根或多根拖缆16(每根包含多个水听器18、一种类型的声传感器)、及记录设备20。当使用地震勘测***10获得地震数据时,声源14(一般一个或多个***)产生一个声脉冲,该声脉冲经海水22传播,并且由海底24和由在其中层的声阻抗变化的各种地质层之间的界面26部分反射。
当使用这种类型的地震勘测***10获得地震数据时,记录的最显著的反射事件一般是直接水底到达28。直接水底到达28也能看作是由声源14产生的直接反射离开海底24且由水听器18直接接收的声脉冲的一个第一部分28。声脉冲的一个第二部分30反射离开海底24,反射离开水听器18上方的海面32,及然后由水听器接收。该第二部分30(称作接收机幻影信号)相对于第一部分28不同相,并且将以接收机幻影陷频导致接收声信号的破坏性干扰,如以上讨论的那样。声脉冲的一个第三部分34在反射离开海底24和由水听器18接收之前在声源14的正上方反射离开海面32。该第三部分34类似地相对于第一部分28不同相,并且以源幻陷频导致接收声信号的破坏性干扰。尽管第一部分28、第二部分30、和第三部分34表示成射线路径,但这仅仅为了便利,并且应该理解声信号一般以球形方式扩散,并且类似的射线路径存在于声源14与水听器18的每一个之间。地震勘测船12一般也装有一个允许监视船下面的水深度的深度发声器36。
图2表示在时间的特定瞬间来自图1的地震拖缆16和海面32的放大部分。海面32由一根平均海平面38基准虚线叠置。平均海平面38一般补偿潮汐效应,并且将理想地代表如果在该区域的海浪完全平静则存在的海平面。海面高度变化40量化在水听器18正上方的海面32高于或低于该平均海平面38的距离。类似地,拖缆16和水听器18由一根希望拖缆深度42基准虚线叠置。该希望拖缆深度42一般表示在平均海平面38下面打算拖拉拖缆16的深度。然而,即使反射器和吊舱运行正常,也要理解,拖缆16波动到一定程度,并且在拖缆内的水听器18不会精确地定位在希望拖缆深度42处。传感器深度变化44量化水听器18高于或低于希望拖缆深度42的距离。在图2中为了说明目的已经夸大了海面高度变化40和传感器深度变化44的程度。也表示在图2中的是一个偏移角46,偏移角46也称作θ,表示在直接水底到达28射线路径与一根竖直射线路径之间的角度。该偏移角46以后用来描述在本发明最佳实施例中进行的计算。
图3是地震轨迹图片48,表示已经由在单根拖缆16上的水听器18从由一个声源14产生的单个声脉冲接收的多条地震轨迹。熟悉本专业的技术人员将理解,这些地震轨迹反映由水听器18的每一个在时间上测量的水压力振幅的变化。在该地震轨迹图片48中第一和最明显的事件是直接水底到达50。在这种情况下,由声源14产生的一个声脉冲到达海底24并且反射到离地震勘测船12最近的水听器18所需要的时间近似为4.2秒(在该例子中水深度约为3100米)。该同一脉冲的直接水底到达部分到达离地震勘测船12最远的水听器18所需要的时间刚超过5秒。注意,海上状态是约2米的显著浪高,并且当获得该特定地震轨迹图片48时,在远偏移处的直接水底到达射线路径的入射角约为35度。
图4表示在已经时差校正直接水底到达50从而他们近似同时对准并且已经选择每条轨迹的相同32毫秒部分(包括直接水底到达50)之后来自图3的地震轨迹的开窗口部分52。没有发生来自图3的开窗口部分的“延伸”。轨迹还没有正常时差校正,试图创建一条等效零偏移轨迹。时间窗口仅已经时间移动,从而直接水底到达50对于轨迹的每一条近似出现在窗口中的相同地方。
应该理解,当正在获得地震数据时,在水听器18正上方的局部海面高度恒定地变化。该问题的时间变化本质使对于整个地震轨迹得到的幻影凹痕模糊。为了解决该问题,选择一个短时间窗口(最好小于约250毫秒),以在声传感器附近在相对短的时间段上提供水柱高度的“快相”。在该最佳实施例中,这是刚好在由水听器18接收到直接水底到达50之前和之后的时间段。然而,这些短时间窗口的使用能使另外的子波提取法(能用来确定幻影陷频)失效,由于(例如)静止子波假设的失效。振幅与频谱的关系最好代之以用来对于选择时间窗口的每一个确定接收机幻影陷频。
图5是曲线图,表示对于来自图4的地震轨迹(具有0至800米偏移的轨迹)的多个开窗口部分振幅与频谱的关系。图5表示源陷频58如何对于轨迹的每一条虚拟地恒定(因为地震轨迹来自单个声源脉冲,所以源陷频只按到达入射角的余弦变化),同时接收机陷频56显著变化,因为在拖缆16中水听器18的每一个上方的水柱高度主要由于波浪影响而不同。
如果源深度靠近接收机深度,则两个幻影凹痕在频率上靠近,并且可能难以区分他们和进行准确的拾波。这最简单地通过把源和接收机放置在不同深度处解决。需要的深度间隔不能简单地定义,因为它取决于主要由浪高确定的凹痕变化。如果把源和接收机在深度上隔开最大波峰对波谷浪高(显著浪高)的至少一半,那么凹痕重叠从来就是不可能的。
也能使用用来区分源和接收机凹痕的其他方法。源陷频58的变化一般小于接收机陷频56的变化,因为***悬浮***“骑在”波浪上,保持源幻影时段更一致。而且,通过工作在公共射击区域,源陷频58对于与相同射击有关的所有接收机将是恒定的(尽管它按入射角的余弦变化)。视在幻影凹痕时段随偏移按入射与竖直的角的余弦减小。该时差曲线需要在转换成局部海面高度之前,或者经使用偏移和水深度计算、或者通过经拾波拟合一根最好吻合的光滑曲线(如一根余弦型曲线或一根二阶或三阶多项式曲线)从凹痕拾波除去。偏移也将影响颠倒海面点的横向位置,由于反射点虚拟地已不在接收机的正上方。如下面讨论的那样,容易校正带有偏移的横向移动。
如以上讨论的那样,计算声传感器上方水柱高度的接收机幻影频率的转换十分简明。如果选择第一非零接收机幻影陷频,则能按如下估计水柱的高度:
WCH=(ATV)/(2*RGNF*COS(THETA))
其中:
WCH=水柱高度
ATV=声传播速度
RGNF=第一非零接收机幻影陷频
THETA=到达入射角
如果例如,一个具有400米偏移的水听器18具有2度的到达入射角(THETA=2度,THETA的余弦是.9994),则在160Hz(RGNF=160)处发现第一非零接收机幻影陷频,并且海水的声传播速度是每秒1460米(ATV=1460),然后通过如下计算可以估计水听器上方的水柱高度:
WCH=(1460)/(2*160*.9994)=4.6米
源和接收机的近似深度、水深度、及偏移距离一般将是已知的,允许十分准确地计算THETA值。
尽管以上公式提供了水听器18上方水柱高度的估计,但对于是比水听器更靠近源的WCH*TAN(THETA)的点,更精确地估计通过水听器的一个水平面上方的水柱高度,因为这是在其处假定声脉冲的第二部分30已经反射离开海面32的横向位置。
以上方法允许由确定的接收机幻影陷频56确定在一个具体声传感器附近的水柱高度的估计。如果这些计算对于在一根拖缆16中的多个不同水听器18进行,并且进行水听器18是水平的假设,则也可以确定局部海面状态。这种方法的改进也允许确定在水听器18的实际位置与该水平对准假设之间的偏离。
如果对于布置在离声源不同偏移处的一组声传感器能辨别到达时间的变化,并且如果能确定在这些辨别时间变化、与如果声传感器以一种具体深度分布定位则期望的时间变化之间的时间差,那么有可能把这些确定的时间差转换成在假定传感器分布与实际深度分布之间的深度差。
图6是综合图表60,表示期望的水底反射到达时间时差曲线62和实际水底反射到达时间拾波64。在到达时间拾波64与传播时间时差曲线62之间的偏离程度在图6中已经大地夸大,以便更清楚地描述该过程。传播时间时差曲线62使对于直接水底反射28的期望到达时间与偏移距离联系起来(在声源14与水听器18之间的距离)。期望到达时间与如果传感器18的每一个精确地布置在希望拖缆深度42处则期望的到达时间相对应。
到达时间拾波64对于布置在该具体偏移处的水听器18表示实际直接水底到达时间,并且这些时间拾波由在地震轨迹图片48中所示的地震轨迹形成。尽管有可能使用其他技术辨别这些反射器到达时间(或反射器到达时间的等效变化),如通过使用确定在相邻轨迹或相邻轨迹的开窗口部分之间的到达时间位移的相关技术,但本申请人已经发现,通过对于直接水底到达反***确地拾波子波开始时间能得到良好的结果。最好,应该使到达时间拾波64最接近0.1毫秒,这相当于近似为15厘米的精度。如果以大于0.1毫秒的抽样间隔记录地震信号,则可能必须内插记录值以便以足够的精度适当地选择事件到达时间。应该除去在接收地震信号中的任何波涛噪声,最好用最小相位低阻滤波器。而且最好尽可能靠近子波开始拾波(第一到达、第一波峰),从而变化幻影到达不会干扰时间拾波。如果使用相关技术确定事件到达时间的变化,则最好严格限制相关的轨迹的部分,以避免源和/或接收机幻影干扰事件到达时间差的准确确定。
传播时间时差曲线62对于一种给定声传感器深度分布画出到达时间随变化偏移的期望变化。在该描述的实施例中,对于该方法的虚拟所有拖拉拖缆有关的应用情况就是这样,假定的声传感器深度分布是水平的,即平行于希望拖缆深度42。关于源/接收机偏移距离、源深度、平均接收机深度、反射点深度、及/或声传输速度(或诸速度)的估计或假设能用来约束传播时间时差曲线62。
在最简单的模型中,假定拖缆16是水平的,假定海底24是水平的,假定海水22具有单一的声传播速度,并因此期望的水底传播时间时差曲线62期望是简单的双曲线(即直接水底到达将表示地球物理学家称作“正常时差”或“双曲线时差”的时差)。然后能选择最好拟合实际水底反射到达时间拾波56的双曲线,并且用来在期望反射器到达时间与实际反射器到达时间之间(或等效于,在期望变化与反射器到达时间的实际变化之间)的时间差。在实际中,一条三阶多项式曲线一般拟合为到达时间拾波56,以产生传播时间时差曲线62。三阶多项式系数允许以上辨别的简单假设的某一程度的变化。
也能进行关于适当传播时间时差曲线62的较复杂假设。例如,来自深度发声器36的深度测量能用来产生海底表面的轮廓。该海底轮廓及声源14和水听器18的设计布局能用来产生一条期望的传播时间时差曲线62。基于该模型计算的传播时间时差曲线62能进一步拟合到海底到达拾波62上。可以类似地模型化影响期望传播时间时差值的物理参数,如缓慢变化的海底结构、水速度随深度的变化、及子波的直接影响,以确定一条适当的传播时间时差曲线62。通过查找跨过从不同声脉冲得到的结果的相关,也有可能辨别这些参数。
使用如下公式能计算水听器18高于或低于希望拖缆深度42的距离:
AASE=EASE+((ARAT-ERAT)*ATV)/COS(THETA)
其中:
AASE=实际声传感器高程
EASE=期望声传感器高程
ARAT=实际反射到达时间
ERAT=期望反射到达时间
ATV=期望传播速度
THETA=到达入射角
如果例如,在2000米偏移处的水听器18具有5米(EASE-5)的期望高程,对于2000米偏移来自传播时间时差曲线54的期望水底反射到达时间是4.4050秒(ERAT 4.4050),从由水听器获得的地震轨迹辨别的实际水底反射到达时间是4.4046秒(ARAT 4.4046),海水的声传播速度是每秒1460米(ATV=1460),及到达入射角是20度(THETA=20度,并且THETA的余弦=.9397),那么可以按如下计算实际声传感器高程:
AASE=-5+((4.4046-4.4050)*1460)/.9397=-5.6米
尽管除在源和接收机之外的不同深度处的任何地震事件能用在该方法中(如在感兴趣地质区域的主要地震反射器),但水底反射一般是在地震记录中的最强反射,并且经常有可能把该事件与其他干扰反射隔离。海底也常常较平(或至少在感兴趣的距离上具有较均匀的坡度),并且这简化选择适当传播时间时差曲线62的过程。然后通过把估计水柱高度校正传感器高于或低于希望拖缆深度42的距离,可以把实际声传感器高程用来“微调”以上讨论的局部海面高度的估计。
图7是综合图表66,表示对于图3中所示的地震轨迹的一些按照本发明确定的局部海面高度68和拖缆深度70。这些轨迹具有300至725米的偏移,底部刻度表示离开最靠近地震勘测船的传感器的相对距离。该曲线表示在当直接海底到达到达水听器18的每一个时的时刻周围一个短时间段内在传感器上方的海平面高度和在平均海平面下面的地震传感器深度的“快相”画面。对于较大偏移的快相值的时间因此是比对于较小偏移的快相值的时间晚一秒的部分。因为对于局部水柱高度和在希望拖缆深度与实际拖缆深度之间的偏差已经调节了计算局部海面高度68的值,通过简单地在较大偏移范围上平均计算的海面高度,可以确定平均海平面38。对于具有偏移的横向位移已经进行校正(上述的WCH*TAN(THETA)校正),尽管在这些小到达角度(<13度)处,这形成微小差别。打算的拖缆深度42是5米,把吊舱设置成控制正或负1米的深度,并且能看到,按照本发明方法确定的实际拖缆深度70落在这些漂移极限内。
尽管联系点接收机拖缆探测设备(即记录由每单个水听器接收的地震信号)已经描述了发明的方法,并且对于这种类型的数据可以得到较好的解,但该方法也可以用于常规组形成数据(即在记录之前把由一组相邻水听器接收的地震信号加在一起)。在至少一些条件下,当使用单传感器数据时,能期望局部海面高度和拖缆深度由发明的方法以大于正或负50米的精度确定。
发明的方法可以用来在海洋地震数据获得活动期间监视局部浪高和地震传感器的高程。这可能是有用的,例如,监视海面状态以保证与关于天气条件的合同职责相符。也可以用来监视海面状态对诸如拖缆、反射器和吊舱之类的海中设备的影响。确定的传感器深度值也能用来控制传感器***件,如反射器和吊舱。
该方法在以后地震数据处理用途中也可能是有用的。可以重新处理在非理想天气条件下获得的地震数据,并且对于每个地震传感器(或地震传感器组)可以确定局部浪高和地震传感器高程。该信息能用来估计在获得地震数据时存在的海面条件和水听器如何准确地保持在其希望深度处。使用该信息能进一步重新处理地震数据,以便更好地把数据去幻影(即调谐消卷积算子,以便描述在地震传感器与海面之间的实际距离,而不是在地震传感器与海面之间的假定距离),并且对于由在实际高程处而不是在希望高程处的地震传感器引起的时间偏移校正反射器到达时间。
因为使用用常规地震数据获得技术和设备获得的地震数据,能实现该方法,所以该方法能用来显著改进在非最优天气条件下获得的地震数据的成象特性。大多数地震勘测获得在最后不能辨别潜在油气矿区的区域上的大量数据。然而,在这些没有表示显著资源的区域中,便利的可能是,使用发明的方法重新处理地震数据,以便“精调”研究下的矿区的地震图象。
尽管以上描述的本发明的实施例利用直接水底到达(并因此确定在直接水底到达到达每个水听器时的声传感器高程和局部浪高),但本发明不限于这种到达的使用。特别是,在某些情况下,产生水盐度界面反射,并且这种反射能用来代替直接水底到达。另外,其他地震反射器和折射器能与发明的方法一起使用。尽管对于这些其他地震反射器的到达时间一般不容易拾波,并且幻影凹痕不比与直接水底到达有关的那些明显,但如果对于一条具体地震轨迹能确定两个或多个延迟局部浪高和声传感器高程,则这打开了外推诸值以改进用于感兴趣的具体地质目标的图象的可能性。如果我们例如对位于5,000米深度处的油气储层感兴趣,并且由直接海底到达确定的对于一条具体轨迹的声传感器高程和局部浪高值分别是-5.6米和1.2米,及当由位于3,000米处的区域地质反射器确定时对于该轨迹的声传感器高程和局部浪高值分别是-5.8米和1.0米,则可能希望分别使用-5.9米和0.9米的外推值时间移动该轨迹和使之去幻影。以这种方式,可能有可能更准确地成象感兴趣的地质目标。
在某些情况下,特别是对于较大偏移,首先记录的地震事件可能是一种折射的而不是反射的到达。折射到达从源到海床,沿海床并且然后到声接收器传播。这些折射到达将呈现相对于水底到达50讨论的相同类型的接收机幻影凹痕,尽管到达角度θ(theta)是常数。对于折射到达的传播时间时差曲线62作为一阶近似也可以是线性的,而不是双曲线。通过对使用水底反射的上述最佳实施例方法进行较小变更,这些折射到达也可以用来确定局部浪高和地震传感器高程。最好与发明的方法一起使用的地震反射器或折射器事件在地震记录中早于其他类型的到达,以限制由“伪幻影”引起的干扰和由另外的传播模式引起的其他类型的声干扰。
尽管本发明的描述实施例利用诸如***之类的常规地震源作为声源,并且利用诸如水听器之类的常规地震传感器作为声传感器,但本发明不限于这样的器件的使用。也有可能使用一个或多个专用非地震声源,如声波发射器或扬声器。常规***主要与发明的方法一起良好地工作,因为他们一般产生一个宽带声信号。如果专用声源与该方法一起使用,则可能最好把他们放在拖缆16的前部、中部和末端,因为该方法在偏移角度46近似为45度或更小时显得工作最好。
类似地,有可能使用打算接收直接海底声反射到达的专用非地震声传感器。常规水听器与发明的方法一起工作得特别好,因为他们是高度准确的并且具有良好的高频响应。然而,对于相同的空间抽样间隔可能不必得到拖缆深度和局部浪高测量,因为水听器放置在拖缆16中。如果专用传感器和处理设备用来确定传感器上方的局部浪高和传感器的高程,则可能最好以比水听器18大的空间抽样间隔隔开专用传感器,并且以显著比在常规地震数据获得中使用的2或3毫秒抽样间隔小的抽样间隔(如每0.1毫秒)抽样声信号的反射部分。地震信号类似地可能不是用来分析地质地下层的“地震数据”。它可以代之以联系进行上述测量缓慢地获得。尽管使用常规水听器18(水压力传感器)描述了本发明的最佳实施例,但该方法以类似方式与听地器(粒子速度传感器)、加速度表、和其他类型的声传感器一起工作。接收机幻影陷频56直接与直接水底到达28和接收机幻影信号30在到达时间上的差有关。也可以使用利用接收机幻影信号30传播时间测量以估计水柱高度的发明方法的另外实施例。
尽管联系常规拖拉海洋拖缆地震数据获得描绘了本发明的描述实施例,但该方法不限于这种类型的环境。它可以类似地用在海底电缆、竖直采用的地震传感器组、和其他类型的地震数据获得***中。
Claims (11)
1.一种分析在海洋地震勘测期间响应声源的操作从多个地震传感器接收的被接收的地震信号的方法,该方法包括步骤:
在一个特定的地震信号内选择一个时间窗口,该窗口得到在该信号中表示的有关的清晰事件,所述特定的地震信号从一个特定的传感器被接收到;
计算由地震传感器产生的地震信号的振幅/频谱,其中所述振幅/频谱包括一个源幻影陷频和一个接收机幻影陷频;
通过根据接收机幻影陷频和源幻影陷频的相对可变性对所述接收机幻影陷频与所述源幻影陷频进行区分,在所述窗口中由所述地震信号的振幅/频谱确定接收机幻影陷频;及
从所述接收机幻影陷频导出在产生所述特定地震信号的所述特定传感器上方的水柱估计高度。
2.根据权利要求1的方法,其中导出所述水柱估计高度的步骤包括:
从地震传感器接收的地震信号辨别到达时间的变化;
确定在被辨别的到达时间变化和与一个假设的地震传感器的深度分布有关的期望的到达时间变化之间的时间差;
将所述的时间差转换成在所述的假设地震传感器的深度分布和所述的地震传感器的实际深度分布之间的深度差。
3.根据权利要求2所述的方法,其中辨别到达时间变化的步骤包括从所述地震传感器接收的所述地震信号辨别与一个普通反射体有关的到达时间的步骤,以及其中确定时间差的步骤包括的步骤为:对于所述假定地震传感器深度分布,通过把到达时间的期望变化与偏移的变化联系起来的一条传播时间时差曲线,确定在所述被辨别的到达时间与期望时间之间的时间差。
4.根据权利要求2所述的方法,其中所估计的水柱高度与所述深度差结合,以估计在海洋地震勘测时所述地震传感器附近相对于平均海平面的局部浪高。
5.根据权利要求2所述的方法,还包括利用水柱的估计高度产生局部海面状态的显示的步骤。
6.一种实现海洋地震勘测的方法,该方法包括步骤:
分析在所述海洋地震勘测期间响应声源的操作由多个水中地震传感器产生的地震信号,该分析步骤包括步骤:
在一个特定的地震信号内选择一个时间窗口,该窗口得到在该信号中表示的有关的清晰事件,所述特定的传感器用于产生所述的地震信号;
计算由所述声传感器产生的地震信号的振幅/频谱,其中所述振幅/频谱包括一个源幻影陷频和一个接收机幻影陷频;
在所述窗口中由所述地震信号的振幅/频谱确定接收机幻影陷频,其中确定接收机幻影陷频包括根据接收机幻影陷频和源幻影陷频的相对可变性对所述接收机幻影陷频和所述源幻影陷频进行区分;及
从所述接收机幻影陷频导出在产生所述特定地震信号的所述特定传感器上方的水柱估计高度。
7.根据权利要求6的方法,其导出所述估计高度包括:
从地震传感器接收的地震信号辨别到达时间的变化;
确定在被辨别的到达时间变化和与一个假设的地震传感器的深度分布有关的期望的到达时间变化之间的时间差;
将所述的时间差转换成在所述的假设地震传感器的深度分布和所述的地震传感器的实际深度分布之间的深度差。
8.根据权利要求7所述的方法,其中:
辨别到达时间变化的步骤包括从所述地震传感器接收的所述地震信号辨别与一个普通反射体有关的到达时间的步骤;
确定时间差的步骤包括的步骤为:对于所述假定地震传感器深度分布,通过把到达时间的期望变化与偏移的变化联系起来的一条传播时间时差曲线,确定在所述辨别到达时间与期望时间之间的时间差。
9.根据权利要求7所述的方法,其中所估计的水柱高度与所述深度差结合,以估计在海洋地震勘测时所述地震传感器附近相对于平均海平面的局部浪高。
10.根据权利要求9所述的方法,其中所述深度差用于控制一个反射器和一个控制所述地震探测器深度的吊舱。
11.根据权利要求10所述的方法,其中所述声源布置在一个其深度不同于所述地震传感器的位置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GB9828066.2 | 1998-12-18 | ||
GBGB9828066.2A GB9828066D0 (en) | 1998-12-18 | 1998-12-18 | Seismic signal analysis method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1330771A CN1330771A (zh) | 2002-01-09 |
CN1181358C true CN1181358C (zh) | 2004-12-22 |
Family
ID=10844603
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB998146323A Expired - Fee Related CN1181358C (zh) | 1998-12-18 | 1999-11-18 | 用于局部浪高确定的方法和在海洋地震信号中的声传感器 |
Country Status (10)
Country | Link |
---|---|
US (1) | US6681887B1 (zh) |
EP (1) | EP1147435B1 (zh) |
CN (1) | CN1181358C (zh) |
AU (1) | AU760590C (zh) |
CA (1) | CA2352816A1 (zh) |
DE (1) | DE69933395D1 (zh) |
EA (1) | EA200100676A1 (zh) |
GB (2) | GB9828066D0 (zh) |
NO (1) | NO332303B1 (zh) |
WO (1) | WO2000037968A1 (zh) |
Families Citing this family (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2379741B (en) * | 2001-09-18 | 2003-11-19 | Westerngeco Ltd | Method for reducing the effect of Sea-surface ghost reflections |
GB2389183B (en) * | 2002-05-28 | 2006-07-26 | Westerngeco Ltd | Processing seismic data |
GB2393513A (en) * | 2002-09-25 | 2004-03-31 | Westerngeco Seismic Holdings | Marine seismic surveying using a source not having a ghost at a non-zero frequency |
GB2400662B (en) | 2003-04-15 | 2006-08-09 | Westerngeco Seismic Holdings | Active steering for marine seismic sources |
GB2405473B (en) * | 2003-08-23 | 2005-10-05 | Westerngeco Ltd | Multiple attenuation method |
US7336561B2 (en) * | 2004-09-07 | 2008-02-26 | Pgs Americas, Inc. | System for attenuation of water bottom multiples in seismic data recorded by pressure sensors and particle motion sensors |
US20060133202A1 (en) * | 2004-12-22 | 2006-06-22 | Tenghamn Stig R L | Motion sensors in a marine seismic streamer |
US7660191B2 (en) * | 2005-07-12 | 2010-02-09 | Westerngeco L.L.C. | Methods and apparatus for acquisition of marine seismic data |
US20070223308A1 (en) * | 2006-03-21 | 2007-09-27 | Frivik Svein A | Methods of range selection for positioning marine seismic equipment |
FR2909457B1 (fr) * | 2006-12-01 | 2009-01-23 | Thales Sa | Procede d'elimination de sources fantomes pour un sonar passif comportant plusieurs antennes lineaires. |
US8593907B2 (en) * | 2007-03-08 | 2013-11-26 | Westerngeco L.L.C. | Technique and system to cancel noise from measurements obtained from a multi-component streamer |
US8077543B2 (en) * | 2007-04-17 | 2011-12-13 | Dirk-Jan Van Manen | Mitigation of noise in marine multicomponent seismic data through the relationship between wavefield components at the free surface |
US8060314B2 (en) * | 2007-04-19 | 2011-11-15 | Westerngeco L. L. C. | Updating information regarding sections of a streamer that are in a body of water |
US9279899B2 (en) * | 2007-07-18 | 2016-03-08 | Westerngeco L.L.C. | System and technique to estimate physical propagation parameters associated with a seismic survey |
US7554881B2 (en) * | 2007-07-31 | 2009-06-30 | Westerngeco L.L.C. | Determining seawater conductivity for the compensation of electromagnetics-based survey data |
US8019549B2 (en) * | 2008-12-10 | 2011-09-13 | Honeywell International Inc. | Event-based power management for seismic sensors |
US7957222B2 (en) * | 2007-10-05 | 2011-06-07 | Honeywell International, Inc. | Acoustic communication and control for seismic sensors |
US8964500B2 (en) * | 2007-10-05 | 2015-02-24 | Honeywell International Inc. | Communication in a seismic sensor array |
US8456948B2 (en) * | 2008-06-28 | 2013-06-04 | Westerngeco L.L.C. | System and technique to obtain streamer depth and shape and applications thereof |
US9207349B2 (en) | 2008-08-27 | 2015-12-08 | Pgs Geophysical As | Methods and apparatus for determining seismic streamer array geometry and seismic sensor response |
US20100054080A1 (en) * | 2008-08-27 | 2010-03-04 | Guillaume Cambois | Determining Seismic Streamer Array Geometry And Seismic Sensor Response Using Dual Sensor Seismic Streamer Arrays |
US20120147700A1 (en) * | 2010-12-14 | 2012-06-14 | Svein Arne Frivik | Determining Streamer Depth and Sea Surface Profile |
US20130023131A1 (en) * | 2011-07-20 | 2013-01-24 | Cggveritas Services Sa | Connecting part between lead-in and marine streamer and method |
US20130028048A1 (en) * | 2011-07-25 | 2013-01-31 | Soellner Walter | Methods and apparatus for seismic imaging which accounts for sea-surface variations |
US9081118B2 (en) * | 2011-11-21 | 2015-07-14 | Cggveritas Services Sa | Device and method for computing depth velocity variations |
US9103943B2 (en) | 2011-11-28 | 2015-08-11 | Fugro-Geoteam As | Acquisition and processing of multi-source broadband marine seismic data |
US9766355B2 (en) * | 2013-05-09 | 2017-09-19 | Robert H. Brune | Use of vector rotational measurements and vector pressure gradient measurements to enhance spatial sampling of dual-sensor water bottom seismic data |
AU2014280832B2 (en) * | 2013-06-10 | 2016-08-11 | Downunder Geosolutions Pty Ltd | Seismic data spectrum restoring and broadening |
US20140379266A1 (en) * | 2013-06-25 | 2014-12-25 | Westerngeco L.L.C. | Processing survey data containing ghost data |
US10598807B2 (en) * | 2014-02-18 | 2020-03-24 | Pgs Geophysical As | Correction of sea surface state |
GB2525072B (en) * | 2014-02-18 | 2020-08-19 | Pgs Geophysical As | Correction of sea surface state |
US10274624B2 (en) * | 2015-09-24 | 2019-04-30 | Magseis Ff Llc | Determining node depth and water column transit velocity |
CN105785435B (zh) * | 2016-03-29 | 2018-08-03 | 中国石油天然气集团公司 | 地震物理模型实验***及其模拟数据采集方法和装置 |
US10324210B2 (en) | 2016-06-30 | 2019-06-18 | Schlumberger Technology Corporation | Method and apparatus for determining rough sea topography during a seismic survey |
US10571589B2 (en) * | 2016-08-17 | 2020-02-25 | Pgs Geophysical As | Constraint of dithering of source actuations |
US20190018124A1 (en) * | 2017-03-22 | 2019-01-17 | AquaFusion, Ltd. | Echo sounding apparatus and echo sounding method |
US11815641B2 (en) | 2020-12-04 | 2023-11-14 | Pgs Geophysical As | Composite far offset impulsive source activations for marine seismic surveying and processing |
CN113552628B (zh) * | 2021-07-20 | 2023-08-15 | 黄河勘测规划设计研究院有限公司 | 地震涌浪高度计算方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4353121A (en) * | 1980-07-24 | 1982-10-05 | Fairfield Industries, Inc. | High resolution, marine seismic stratigraphic system |
EP0057675B1 (en) * | 1980-07-24 | 1987-02-25 | Fairfield Industries, Inc. | High resolution, marine seismic stratigraphic system |
US4375089A (en) * | 1980-10-27 | 1983-02-22 | Western Geophysical Company Of America | Depth transducer system for a seismic streamer cable |
US4520467A (en) * | 1982-03-18 | 1985-05-28 | Shell Oil Company | Marine seismic system |
US5142498A (en) * | 1991-08-28 | 1992-08-25 | Exxon Production Research Company | Controlled phase marine source array |
US5691957A (en) * | 1994-06-30 | 1997-11-25 | Woods Hole Oceanographic Institution | Ocean acoustic tomography |
NO971151D0 (no) * | 1996-05-21 | 1997-03-12 | Western Atlas Int Inc | Fremgangsmåte for måling av vanndybde over en bunnkabel |
-
1998
- 1998-12-18 GB GBGB9828066.2A patent/GB9828066D0/en not_active Ceased
-
1999
- 1999-11-18 CN CNB998146323A patent/CN1181358C/zh not_active Expired - Fee Related
- 1999-11-18 WO PCT/GB1999/003837 patent/WO2000037968A1/en active IP Right Grant
- 1999-11-18 EP EP99973497A patent/EP1147435B1/en not_active Expired - Lifetime
- 1999-11-18 AU AU11704/00A patent/AU760590C/en not_active Ceased
- 1999-11-18 GB GB9927300A patent/GB2344889B/en not_active Expired - Fee Related
- 1999-11-18 US US09/868,132 patent/US6681887B1/en not_active Expired - Fee Related
- 1999-11-18 CA CA002352816A patent/CA2352816A1/en not_active Abandoned
- 1999-11-18 EA EA200100676A patent/EA200100676A1/ru unknown
- 1999-11-18 DE DE69933395T patent/DE69933395D1/de not_active Expired - Lifetime
-
2001
- 2001-06-18 NO NO20013030A patent/NO332303B1/no not_active IP Right Cessation
Also Published As
Publication number | Publication date |
---|---|
EP1147435B1 (en) | 2006-09-27 |
AU1170400A (en) | 2000-07-12 |
DE69933395D1 (de) | 2006-11-09 |
GB2344889A (en) | 2000-06-21 |
EA200100676A1 (ru) | 2001-12-24 |
US6681887B1 (en) | 2004-01-27 |
GB9927300D0 (en) | 2000-01-12 |
AU760590B2 (en) | 2003-05-15 |
GB2344889B (en) | 2001-04-04 |
WO2000037968A1 (en) | 2000-06-29 |
AU760590C (en) | 2004-07-15 |
CA2352816A1 (en) | 2000-06-29 |
NO20013030L (no) | 2001-08-20 |
NO20013030D0 (no) | 2001-06-18 |
CN1330771A (zh) | 2002-01-09 |
NO332303B1 (no) | 2012-08-20 |
EP1147435A1 (en) | 2001-10-24 |
GB9828066D0 (en) | 1999-02-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1181358C (zh) | 用于局部浪高确定的方法和在海洋地震信号中的声传感器 | |
US9696443B2 (en) | Method and device for processing seismic data | |
US7684281B2 (en) | System for combining signals of pressure sensors and particle motion sensors in marine seismic streamers | |
US7957906B2 (en) | Method for attenuating low frequency noise in a dual-sensor seismic streamer | |
US10690792B2 (en) | Amplitude-versus-angle analysis for quantitative interpretation | |
US9477000B2 (en) | System and method for the removal of shallow water multiples using a hybrid multi-channel prediction method | |
MX2014003060A (es) | Sistemas y metodos de filtracion de dominio de frecuencia y discriminacion de dominio de espacio-tiempo de datos sismicos. | |
AU2014274533B2 (en) | Methods and systems for optimizing generation of seismic images | |
US9733375B2 (en) | Method and device for alternating depths marine seismic acquisition | |
RU2282877C2 (ru) | Способ корректировки сейсмических данных при морской сейсмической разведке | |
EA016522B1 (ru) | Способ оценки фактора качества формации | |
NO346705B1 (no) | Dempning av støy ved skuddgjentakelse ved marin seismisk kartlegging av undergrunnen | |
AU2015205965B2 (en) | Methods and systems that combine wavefields associated with generalized source activation times and near-continuously recorded seismic data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
C10 | Entry into substantive examination | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20041222 Termination date: 20121118 |