CN101487899B - 三维双传感器拖缆数据中的波场分离方法 - Google Patents
三维双传感器拖缆数据中的波场分离方法 Download PDFInfo
- Publication number
- CN101487899B CN101487899B CN200910002704.1A CN200910002704A CN101487899B CN 101487899 B CN101487899 B CN 101487899B CN 200910002704 A CN200910002704 A CN 200910002704A CN 101487899 B CN101487899 B CN 101487899B
- Authority
- CN
- China
- Prior art keywords
- record
- wave
- uplink
- vibration velocity
- vertical particle
- 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 91
- 230000009977 dual effect Effects 0.000 title claims abstract description 16
- 238000000926 separation method Methods 0.000 title abstract description 16
- 238000001228 spectrum Methods 0.000 claims abstract description 83
- 239000002245 particle Substances 0.000 claims abstract description 74
- 230000009021 linear effect Effects 0.000 claims abstract description 9
- 230000008569 process Effects 0.000 claims description 38
- 238000005070 sampling Methods 0.000 claims description 25
- 238000006243 chemical reaction Methods 0.000 claims description 20
- 238000010606 normalization Methods 0.000 claims description 10
- 230000008901 benefit Effects 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims description 3
- 239000013598 vector Substances 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 19
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 238000009825 accumulation Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000002360 explosive Substances 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000010891 electric arc Methods 0.000 description 1
- 238000010892 electric spark Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000009329 sexual behaviour Effects 0.000 description 1
- 230000002123 temporal effect 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- 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
- 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
来自双传感器拖缆数据的压力记录和垂直质点振速记录被变换到主测线波数域。在每个主测线波数上一系列缩放滤波器被应用到经过变换的垂直质点振速记录,其中缩放滤波器系列中的每个针对不同的交叉拖缆波数范围、并在其中地震事件近似线性的主测线道的块中被计算。压力频谱和经过缩放的垂直质点振速频谱被组合来分离上行和下行波场分量。被分离的上行和下行波场分量被逆变换回到时空域。
Description
技术领域
本发明一般地涉及地球物理勘探的领域。更具体地,本发明涉及海洋地震数据处理。
背景技术
在油气工业中,地球物理勘探一般被用于帮助对地下岩层的搜索和评估。地球物理勘探技术产生(关于)地球地下结构的知识,这对寻找和提取有价值的矿产资源,特别是烃类沉积如石油和天然气很有用。地球物理勘探的一种广为人知的技术是地震勘测。在基于陆地的地震勘测中,地震信号在地表上或附近产生,然后向下传播到地表之下。在海洋地震勘测中,地震信号也可以经过在地表下覆盖的水层向下传播。地震能量源被用来产生地震信号,该信号在传播到地球内部之后,至少部分被地下的地震反射体反射。这种地震反射体典型地是地下岩层之间的分界层,这些岩层具有不同弹性属性,特别是声波速率和岩石密度,其导致分界层上的不同声阻抗。反射的地震能量被地表上或附近、在覆盖的水体中、或凿洞中的已知深度的地震传感器(也被称为地震接收器)检测到并被记录。
用于在陆地地震勘测中产生地震信号的合适地震源可以包括爆炸物或振动器。海洋地震勘测典型地使用由船只拖拽、并被周期性地激活来产生声波场的水下地震源。产生波场的地震源可以是几种类型,包括小爆炸弹药量、电火花或电弧、海洋振动器、以及典型地,气枪。地震源枪可以是水枪、蒸汽枪、且最典型地,气枪。典型地,海洋地震源不只包括单个源部件,而是空间分布的阵列的源部件。这种布置对于气枪由其如此,气枪是海洋地震源当前最常用的形式。
地震传感器的合适类型包括质点振速传感器,特别是在陆地勘测中,以及水压传感器(典型地,水压梯度传感器),特别是在海洋勘测中。有时,质点加速度传感器代替或在质点振速传感器之外被使用。质点振速传感器和水压传感器分别在地震检波器和水听器领域普遍为人所知。地震传感器可以单独部署,更一般地部署为传感器阵列。此外,压力传感器和质点振速传感器可以在海洋勘测中一起部署,成对或阵列成对地布置。
在典型的海洋地震勘测中,地震勘测船在水面上行进,典型地以5节的速度,并装载了地震获取装置,例如航海控制器、地震源控制器、地震传感器控制器,以及记录装置。地震源控制装置使由地震勘测船在水体中拖行的地震源在选定的时间启动。地震拖缆,也被称为地震缆线,是由拖行地震源的地震勘测船或另一地震勘测船在水体中拖行的类似缆线的伸长结构。常规地,地震拖缆包括压力传感器例如水听器,但建议地震拖缆在水听器之外包括水质点振速传感器例如地震检波器或质点加速度传感器例如加速计。压力传感器和质点运动传感器被接近地部署,沿着地震缆线成对或阵列成对地布置。
执行勘测得到的结果地震数据被处理,以产生与勘测区域的地质结构和地下岩层属性相关的信息。被处理的地震数据被处理用于显示和分析这些地下岩层中的潜在烃类含量。地震数据处理的目标是从地震数据中提取尽可能多关于地下岩层的信息,以充分描绘地质地表下(结构)。为了识别地球的地表下可能找到石油积聚的位置,大量金额的资金花费在搜集、处理和解析地震数据上。从记录的地震数据构造定义感兴趣的地下地层的反射层,提供了地球深度或时间上的图像。
地球地表下的结构的图像被产生,以使解析器(interpreter)选择最有可能具有石油积聚的位置。为了验证石油的存在,必须钻井。钻井以确定石油储藏是否存在,是非常昂贵和耗时的任务。为此,存在持续的需求来改善对地震数据的处理和显示,从而产生地球地表下结构的图像,其可以提高解析器的能力,无论该解析是由电脑或人完成,以估算在地球地表下特定位置存在石油积聚的可能性。
双传感器拖缆反射地震数据由压力场和垂直质点振速场记录构成。对地震数据的处理链的中心要素是它被分离为仅包括压力波场的上行和下行分量的记录。该分离可以在将数据变换到频率-波数(f-kx-ky)域,分别取压力记录的频率-波数频谱和垂直质点振速记录的频率-波数频谱的经过缩放的形式之间的差值及之和两者,并将得到的频谱除以2后执行。(注意,通过简单使用前面的缩放滤波器的逆,可以可替换地通过分别取垂直质点振速记录的频率-波数频谱和压力记录的频率-波数频谱的经过缩放的形式之间的差值及之和,并将得到的频谱除以2,得到垂直质点振速波场的上行和下行分量)。从频率-波数域逆变换回时空域生成期望的上行和下行波场分量。在该过程中,只有垂直质点振速记录(或可替换地,只有压力记录)被缩放改变。对于非易散性(non-evanescent)能量,其在频率-波数域中被实滤波器缩放,该滤波器对于给定频率随着增加的波数而系统地增加。但是,在与拖缆交叉方向上的空间假频(aliasing)在海洋地震勘测中很普遍。在与拖缆交叉方向假频的情况下,能量被限制(wrap)到较低的交叉拖缆波数ky。随后,如果没有考虑这些波数中的这些卷绕效应,则缩放滤波器从太低的,错误的波数计算。于是,垂直质点振速记录(或可替换地,压力记录)中的假频能量(aliased energy)被一直太低(或分别地,太高)的滤波器系数缩放。
初始记录中上行和下行波场分量的叠加在对应的频率-波数频谱中引起接收器虚反射(ghost)凹槽(notch)的特定模式。每当记录的能量在例如压力场的频谱中的特定频率-波数组合处被抵消,对应的记录的能量在垂直质点振速的频谱中最大。但是,这种对应引起在压力记录(或可替换地,垂直质点振速记录)的频率波数频谱中的凹槽处上行和下行波场分量的不正确分离。在这些位置,对于垂直质点振速记录(或可替换地,压力记录),假频能量被不正确地缩放,从而得到的分离的波场分量被不正确地计算。
从而,需要一种方法来分离3维双传感器拖缆地震数据中的上行和下行波场分量,其可以适当地处理与拖缆交叉方向上的假频能量。
发明内容
本发明是用于分离三维双传感器拖缆数据中的上行和下行波场分量的方法,该数据可能具有与拖缆交叉方向上的假频能量。来自拖缆数据的压力记录和垂直质点振速记录被变换到波数域。在每个主测线波数一系列缩放滤波器被应用到经过变换的垂直质点振速记录,其中一系列缩放滤波器中的每个针对不同的交叉拖缆波数范围、在其中所有地震事件近似线性的主测线道的块中被计算。压力频谱和经过缩放的垂直质点振速频谱被组合,以分离上行和下行波场分量。被分离的上行和下行波场分量被逆变换回时空域。
附图说明
通过参考下列详细描述和附图,本发明及其优势可以被更容易地理解,在附图中:
图1是说明用于分离三维双传感器拖缆数据中的上行和下行波场分量的本发明的第一实施例的处理步骤的流程图;
图2是用于分离三维双传感器拖缆数据中的上行和下行波场分量的本发明的第二实施例的初始处理步骤的流程图;
图3是说明本发明的一个实施例的用于处理来自图2的窗口的中间处理步骤的流程图;并且
图4是说明本发明的一个实施例的处理来自图3的谱比的中间处理步骤的流程图。
图5是说明如图1-4所示用于分离三维双传感器拖缆数据中的上行和下行波场分量的本发明的第二实施例的最后处理步骤的流程图;
图6示出了本发明的方法中使用的示例性第三板的框图;
图7示出了从图6的示例性第三板中的条创建的掩模的框图;
图8示出了针对不同的最大去假频(dealias)阶次L,与本发明的方法的结果相比的参考道的幅度频谱的图;
图9示出了针对不同的最大去假频阶次L,与本发明的方法的更多结果相比的参考道的幅度频谱;
图10示出了由本发明的方法估算的信号以及它与去假频阶次L=0的参考道的差值的图;并且
图11示出了由本发明的方法估算的信号以及它与去假频阶次L=5的参考道的差值的图。
尽管本发明结合其优选实施例被描述,可以理解本发明不限于此。相反,本发明旨在覆盖所附权利要求书定义的本发明的范围内的所有替换、修改和等价物。
具体实施方式
在可以为波场分离正确缩放垂直质点振速记录之前,其频谱需要被合适地去假频。实现该目的的传统方法是在与拖缆交叉方向上的道内插,以降低拖缆间隔并由此增加交叉拖缆Nyquist波数。本发明的方法提供了用于适当处理波场分离期间的假频能量的替换方案,而不需要明显的道内插。
本发明的方法使用快速傅里叶变换(FFT)的周期性属性。但是,不需要为波场分离计算内插道。替代地,在每个主测线波数kx处若干个缩放滤波器被应用到垂直质点振速记录的初始频谱。这些缩放滤波器中的每个是针对不同的交叉拖缆波数范围计算的,并且仅作用于频谱中该部分的能量,该部分是,或在假频能量的情况下应该是,在该波数范围内。本发明的方法中的波场分离在主测线维度上定义的只包括近似线性的地震事件的块中工作。仅为了描述的简单和清楚,在下列讨论中针对单个块来描述过程。
图1-5是说明用于波场分离的本发明的实施例的流程图。图1和2-5分别示出了本发明的方法的两个实施例。图2示出了本发明的方法的第二实施例的初始步骤,图3和图4示出了图2所示方法的进一步的中间步骤,且图5示出了图2-4中示出的方法的最后步骤。
图6-11说明了参考图1-5讨论的流程图中描述的一些步骤。图6和图7分别说明了谱比、第三板及其掩模。图8和图9说明了针对不同的最大去假频阶次,得自本发明的方法的结果的幅度频谱的比较。图10和图11分别说明了针对去假频阶次L=0和5由本发明的方法估算的信号。
图1是说明用于分离在三维双传感器拖缆数据中的上行和下行波场分量的本发明的方法的第一实施例的处理步骤。拖缆数据在与拖缆交叉方向上可能具有假频能量,其可以被本发明的方法处理。
在步骤11中,来自双传感器拖缆数据的压力记录和垂直质点振速记录从时空(t-x-y)域变换到主测线波数(t-kx-y)域,其中t表示时间,x和y分别为主测线和交叉拖缆(联络测线)维度,且kx是主测线波数。
在步骤12中,在每个主测线波数kx上一系列缩放滤波器被应用到来自步骤11的经过变换的垂直质点振速记录。于是,这些缩放滤波器被应用到恒定kx的切片。这些缩放滤波器系列中的每个针对不同的交叉拖缆波数范围、并在其中所有地震事件近似线性的主测线道的块中被计算。这些缩放滤波器中的每个仅作用于频率-波数频谱中适当部分的能量,所述部分在非假频能量(unaliased energy)的对应波数范围中,或应该在假频能量的对应波数范围内。
在步骤13中,压力频谱和来自步骤12的经过缩放的垂直质点振速频谱被组合,以分离上行和下行波场分量。
在步骤14中,来自步骤13的分离的上行和下行波场分量被逆变换回时空(t-x-y)域。
图2的流程图示出用于分离三维双传感器拖缆数据中的上行和下行波场分量的本发明的方法的第二实施例的初始处理步骤。拖缆数据可能具有与拖缆交叉方向上的假频能量。图2在如上参考图1讨论的第一实施例的讨论上扩展。
在步骤20中,压力和垂直质点振速记录在时空(t-x-y)域中被分为块。每个块包括交叉拖缆(y)维度中的所有道以及时间(t)维度中的所有道采样点。但是,块在主测线(x)维度上被限制,使得块中的所有地震事件近似线性。从而,块在主测线维度上可重叠。当块重叠时,它们必须在波场分离后被适当地组合。
在步骤21中,来自步骤20的块被选择。
在步骤22中,在步骤21中选择的块中的压力和垂直质点振速记录从t-x-y域变换到t-kx-y域。
在步骤23中,都来自步骤22的具有恒定kx的两个切片,一个是经过变换的垂直质点振速记录的,一个是对应的压力记录的,被选择。
在步骤24中,在步骤23中选择的具有恒定kx的两个切片被分为若干个窗口,所述窗口包括交叉拖缆维度中所有的道,但不一定包括时间维度中的所有道采样点。这些窗口可在时间维度上重叠。当窗口重叠时,它们必须在波场分离后被合适地组合。
在步骤25中,都来自步骤24的两个对应的窗口,一个来自压力记录,一个来自垂直质点振速记录,被选择。
在步骤26中,最大去假频阶次,由L指定的整数,被选择。初始交叉拖缆波数范围对应于L=0。
在步骤27中,在步骤25中选择的窗口被发送到图3的步骤31用于进一步的处理。
在步骤28中,确定在步骤23选择的具有恒定kx的切片中,是否还有任何更多的窗口要在步骤25中被选择。如果还有更多的窗口要被选择,则过程回到步骤25。如果没有更多的窗口要被选择,则过程进入下一步骤29。
在步骤29中,确定是否还有更多的具有恒定kx的切片要在步骤23中被选择。如果还有更多的切片要被选择,则过程回到步骤23。如果没有更多的切片要被选择,则过程进入下一步骤30。
在步骤30中,确定是否还有更多的块要在步骤21中被选择。如果还有更多的块要被选择,则过程回到步骤21。如果没有更多的块要被选择,则过程对于图2中的流程图结束。现在过程进入图3中的流程图。
图3的流程图示出本发明的一个实施例的用于处理来自图2的窗口的处理步骤。每个这样的窗口经历以下的流程。
在步骤31中,两个选择的窗口从图2的步骤27得到,一个来自压力记录,一个来自垂直质点振速记录。
在步骤32中,为当前kx以及期望的交叉拖缆波数范围计算要在波场分离中应用的缩放滤波器,所述期望的交叉拖缆波数范围是数据窗口给出的Nyquist波数范围的L倍。
在步骤33中,来自步骤31的垂直质点振速记录的窗口在交叉拖缆维度中,以其道数的L倍被零道填充。这里,L是在图2的步骤26中选择的最大去假频阶次。
在步骤34中,来自步骤33的垂直质点振速记录经过填充的窗口在时间维度中,以其道中采样数量的L倍被零采样填充,生成第一扩展记录。这里,L是在图2的步骤26中选择的最大去假频阶次。
在步骤35中,来自步骤34的第一扩展记录在每两个实况道(livetrace)之间被周期性地置零L个道,生成第二扩展记录。
在步骤36中,分别来自步骤34和35的垂直质点振速记录的第一和第二扩展窗口从主测线波数(t-kx-y)域变换到频率-波数(f-kx-ky)域。分别地,来自步骤34的第一扩展记录的频谱在下面被称为“第一频谱”,而来自步骤35的第二扩展记录的频谱被称为“第二频谱”。频率和波数范围现在具有在图2的步骤25中选择的窗口的初始频谱中存在的采样的(L+1)倍数量的采样。分别地,在初始频谱的频率波数域中的采样数量对于频率维度被称为nf,而对于交叉拖缆波数维度被称为nk。
在步骤37中,来自步骤36的第一频谱和第二频谱两者的频率范围都被缩小,仅考虑所有波数最内部的nf个抽样。
在步骤38中,少量白噪声被加到来自步骤37的被缩小的第二频谱。增加白噪声是为了防止在下一步骤39中第二频谱作为除数使用时太小的问题。
在步骤39中,取来自步骤37的第一频谱和来自步骤38的第二频谱的比率。第一和第二频谱的结果比率在下面被称为“第三板”。
在步骤40中,来自步骤32的缩放滤波器和来自步骤39的第三板被发送到图4的步骤41用于进一步的处理。过程对于图3的流程图结束,且过程现在进入图4的流程图。
图4是说明用于处理来自图3的谱比的本发明的方法的中间处理步骤的流程图。
在步骤41中,缩放滤波器和第三板从图3的步骤40得到。
在步骤43中,每两个对应的条组合为一个掩模,将具有正L索引的条移到初始正波数范围,并将具有负L索引的条移到初始负波数范围。如果L是奇数,交换两条的位置。在步骤42和43中,(L+1)个掩模被创建,每个具有nf乖nk个采样。
在步骤45中,每两个对应的条被组合为一个滤波器,将具有正L索引的条移到初始正波数范围,并将具有负L索引的条移到初始负波数范围。如果L是奇数,交换两条的位置。在步骤44和45中,(L+1)个滤波器被创建,对应于在步骤42和43中创建的(L+1)个掩模,且作用于掩模的不同波数范围上。
在步骤46中,由在相同频率波数采样处的所有掩模的采样构成的所有矢量被规格化。该规格化是确保在分离的波场中不引入额外能量所必需的。
在步骤47中,对应的滤波器和掩模乘在一起,以创建各缩放滤波器系列。
在步骤48中,在步骤47中创建的各缩放滤波器系列被发送到图5的步骤51用于进一步的处理。
图5是说明如图1-4所示的用于分离三维双传感器拖缆数据中的上行和下行波场分量的本发明的方法的第二实施例的最后处理步骤的流程图。
在步骤51中,各缩放滤波器系列从图4的步骤48得到。
在步骤52中,来自步骤51的各缩放滤波器系列被应用到来自图2的步骤25的垂直质点振速记录的窗口的初始频谱。
在步骤53中,对来自步骤52的垂直质点振速记录的窗口的所有得到的经过滤波的频谱求总和。该总和生成该窗口中的垂直质点振速记录的经过缩放的频谱。
在步骤54中,来自图2的步骤25的窗口中的压力记录的频谱和来自步骤53的窗口中的垂直质点振速记录的经过缩放的频谱被组合。这产生窗口中上行和下行波场的频谱。
在步骤55中,来自步骤54的上行和下行波场的所有频谱从频率-波数(f-kx-ky)域逆变换到主测线波数(t-kx-y)域。
在步骤56中,针对恒定kx的切片中的来自图2的步骤24的所有窗口,来自步骤55的逆变换的频谱被组合。该组合产生切片中的上行和下行波场。
在步骤57中,来自步骤56的,包括从来自图2的步骤23的切片计算的上行和下行波场的所有切片在块中组合。该组合产生t-kx-y域中的块中的上行和下行波场的记录。
在步骤58中,来自步骤57的上行和下行波场分量从主测线波数(t-kx-y)域逆变换回时空(t-x-y)域。该变换产生块中的时空域中的上行和下行波场分量。
在步骤59中,从步骤58得到的块中的上行和下行波场分量被组合,产生最后的时空域中的上行和下行波场分量。
通过也为压力记录对每个窗口计算掩模系列、并在规格化步骤之前将它与用于垂直质点振速记录的系列相组合,本发明的方法可以进一步稳定。这样,两种记录的共模信号被强化,并且随机噪声的影响被降低。本方法目的在于从其在较低频率的比率预测假频和非假频能量的比率。本发明的方法在时间窗口中的应用是被期望的,因为本方法假设有限数量的局部线性事件。
上面描述的本发明的方法针对其中垂直质点振速记录被缩放而压力记录没有被缩放的实施例。可替换地,本发明的方法包括可替换的实施例,其中压力记录被缩放而垂直质点振速记录没有。上述讨论可以直接被修改,以包括该可替换的实施例。仅需要交换对压力记录和垂直质点振速记录的引用。
本发明在上面被作为一种方法讨论,仅用于说明性的目的,但也可以被实现为一种系统。本发明的系统优选地以计算机工具,特别是数字计算机,和其它常规数据处理装置一起来实现。这样的数据处理装置,在该领域广为人知,将包括计算机处理装置的任意合适组合或网络,包括但不限于,硬件(处理器、临时和永久存储设备、以及任何其它合适的计算机处理装置)、软件(操作系统、应用程序、数学程序库、以及任何其它合适的软件)、连接(电、光、无线或另外的)、以及外围设备(输入输出设备例如键盘、点击设备和扫描仪;显示设备例如监视器和打印机;存储媒介例如磁盘和硬盘,以及任何其它合适的装置)。
另外,注意到在上面的图5的步骤54中的分离上行和下行波场分量,可以在上面的图5的步骤52中将缩放滤波器系列应用到窗口中的垂直质点振速频谱之后的任意点被执行。步骤顺序的任意得到的改变在本发明的方法之内。
创建掩模的过程以及提出的方法的效果在后面的图中示出。使用的数据是使用距离为25m的23条拖缆来向前建模的(forward modeled)。该模型由在与拖缆交叉方向上具有60°倾斜而在拖缆方向上没有倾斜的单个反射器构成。后面的图示出了kx=0的切片的结果。没有窗口化在时间维度被应用,因为数据仅包括单个反射事件和对应的接收器虚反射。
图6是说明如本发明的方法中使用的示例第三板的图示。图6示出了在对每个波数减少到nf个抽样之后的示例第三板。用来构造掩模的条61被垂直线62分开。这些条从L=-5到L=5编号和标注。
图7是说明从图6中的示例第三板中的条创建的掩模的图示。图7示出了被重排列为掩模71直到L=4的图6的条。示出了在规格化之后的掩模,该规格化是图4中的步骤46。每个掩模具有和垂直质点振速记录的初始频谱相同的采样数量。
图8-11示出了对于道数量11,建模的无虚反射(ghost-free)压力场(参考)和估算的压力场之间的比较。图8示出了参考道的幅度频谱81与本发明的方法对于不同的最大去假频阶次L的结果82比较的图。图9示出了参考道的幅度频谱91与本发明的方法对于不同的最大去假频阶次L的更多结果92比较的图。选择的去假频阶次越高,直到其估算道的频谱等于参考道的频谱的频率越高。
图10示出了对于去假频阶次L=0,由本发明的方法估算的信号101以及它与参考道的差102的图。图11示出了对于去假频阶次L=5,本发明的方法估算的信号111以及它与参考道的差112的图。
应该理解,上面仅是对本发明的特定实施例的详细描述,并且在不偏离本发明的范围的前提下,可以根据本公开对公开的实施例进行多种改变、修改和替换。因此上面的描述不是要限制本发明的范围。而是,本发明的范围只被所附权利要求及其等价物确定。
Claims (24)
1.一种用于分离三维双传感器拖缆数据中的上行和下行波场分量的方法,包括:
将来自拖缆数据的压力记录和垂直质点振速记录变换到主测线波数域;以及
在每个主测线波数将一系列缩放滤波器应用到经过变换的垂直质点振速记录,其中针对不同的交叉拖缆波数范围并在主测线道的块中计算缩放滤波器系列中的每个,在所述主测线道的块中所有地震事件近似线性;
组合经过变换的压力记录和经过缩放的垂直质点振速记录,以分离上行和下行波场分量;以及
将分离的上行和下行波场分量逆变换回时空域。
2.如权利要求1所述的方法,其中,将来自拖缆数据的压力记录和垂直质点振速记录变换到主测线波数域包括:
分别从安装在拖缆中的压力记录器和垂直质点振速记录器取回压力记录和垂直质点振速记录;
在拖缆数据中选择多个块,每个块包括交叉拖缆维度中的所有道、时间维度中的所有道采样点、以及主测线维度中的仅足够的道,使得每个块中的所有地震事件是近似线性的;
将每个块中的压力记录和垂直质点振速记录变换到主测线波数域;
在每个经过变换的块中选择具有恒定主测线波数的多个切片;
将每个切片分为多个窗口,每个窗口包括交叉拖缆维度中的所有道;以及
对每个窗口执行以下步骤:
为窗口选择最大去假频阶次L;
为主测线波数和交叉拖缆波数的选择范围计算缩放滤波器;
用零道和零采样填充窗口,产生第一扩展记录;
用零采样周期性地置零第一扩展记录,以便产生第二扩展记录;
将第一和第二扩展记录变换到频率波数域,分别产生第一和第二频谱;
缩小第一和第二频谱;以及
将第三板计算为第一频谱与第二频谱的比率。
3.如权利要求2所述的方法,其中,交叉拖缆波数的选择范围是窗口的Nyquist波数范围的L倍。
4.如权利要求2所述的方法,其中,用零道和零采样填充窗口包括:
在交叉拖缆维度中用零道以窗口中道数量的L倍填充窗口;以及
在时间维度中用零采样以道中采样数量的L倍填充窗口。
5.如权利要求2所述的方法,其中,用零采样周期性地置零第一扩展记录包括:
周期性地在每两个实况道之间置零L个道,以产生第二扩展记录。
6.如权利要求2所述的方法,其中,缩小第一和第二频谱包括:
将第一和第二频谱的频率范围缩小到最内部的采样,其频率范围大小等于经过变换的垂直质点振速记录中的频率数量。
7.如权利要求2所述的方法,其中,取第一频谱和第二频谱的比率包括:
在除之前将少量白噪声加到第二频谱。
8.如权利要求2所述的方法,其中,应用一系列的缩放滤波器包括:
将第三板分为2(L+1)条,从-L到L编号;
将每两个对应编号的条组合为掩模;
将缩放滤波器分为2(L+1)条,从-L到L编号;
将每两个对应编号的条组合为滤波器;
规格化掩模;
将对应的滤波器和掩模乘在一起,生成一系列缩放滤波器;以及
将缩放滤波器系列应用到经过变换的垂直质点振速记录。
9.如权利要求8所述的方法,其中,每个条在交叉拖缆波数维度中包含的采样数量是经过变换的垂直质点振速记录中的一半。
10.如权利要求8所述的方法,其中,将每两个对应编号的条组合为滤波器包括:
将具有正L的条移到正波数范围;
将具有负L的条移到负波数范围;以及
如果L是奇数,交换条的位置。
11.如权利要求8所述的方法,其中,规格化掩模包括:
规格化相同频率波数采样处在所有掩模上的采样的所有矢量。
12.如权利要求8所述的方法,其中组合经过变换的压力记录和经过缩放的垂直质点振速记录以分离上行和下行波场分量包括:
对所有经过滤波的经过转换的垂直质点振速记录求总和,生成窗口中的经过缩放的垂直质点振速频谱;
组合窗口中的压力频谱和经过缩放的垂直质点振速频谱,产生窗口中的上行和下行波场分量的频谱;
将上行和下行波场分量的频谱从频率波数域逆变换到主测线波数域;
将所有窗口中的上行和下行波场分量的经过逆变换的频谱组合为恒定kx的切片,生成切片中的上行和下行波场分量,其中kx是主测线波数;
将所有切片中的上行和下行波场分量组合为块,生成块中的上行和下行波场分量;
将块中的上行和下行波场分量逆变换到时空域;以及:
组合所有块中的上行和下行波场分量,生成分离的上行和下行波场分量。
13.一种用于分离三维双传感器拖缆数据中的上行和下行波场分量的系统,包括:
用于将来自拖缆数据的压力记录和垂直质点振速记录变换到主测线波数域的装置;和
用于在每个主测线波数将一系列缩放滤波器应用到经过变换的垂直质点振速记录的装置,其中针对不同的交叉拖缆波数范围并在主测线道的块中计算缩放滤波器系列中的每个,在所述主测线道的块中所有地震事件近似线性;
用于组合经过变换的压力记录和经过缩放的垂直质点振速记录,以分离上行和下行波场分量的装置;以及
用于将分离的上行和下行波场分量逆变换回时空域的装置。
14.如权利要求13所述的系统,其中,用于将来自拖缆数据的压力记录和垂直质点振速记录变换到主测线波数域的装置包括:
用于分别从安装在拖缆中的压力记录器和垂直质点振速记录器取回压力记录和垂直质点振速记录的装置;
用于在拖缆数据中选择多个块的装置,每个块包括交叉拖缆维度中的所有道、时间维度中的所有道采样点、以及主测线维度中的仅足够的道,使得每个块中的所有地震事件是近似线性的;
用于将每个块中的压力记录和垂直质点振速记录变换到主测线波数域的装置;
用于在每个经过变换的块中选择具有恒定主测线波数的多个切片的装置;
用于将每个切片分为多个窗口的装置,其中每个窗口包括交叉拖缆维度中的所有道;以及
用于对每个窗口执行下述的装置:
用于为该窗口选择最大去假频阶次L的装置;
用于为主测线波数和交叉拖缆波数的选择范围计算缩放滤波器的装置;
用于用零道和零采样填充窗口,产生第一扩展记录的装置;
用于用零采样周期性地置零第一扩展记录,以便产生第二扩展记录的装置;
用于将第一和第二扩展记录变换到频率波数域,分别产生第一和第二频谱的装置;
用于缩小第一和第二频谱的装置;以及
用于将第三板计算为第一频谱与第二频谱的比率的装置。
15.如权利要求14所述的系统,其中,交叉拖缆波数的选择范围是窗口的Nyquist波数范围的L倍。
16.如权利要求14所述的系统,其中,用于用零道和零采样填充窗口的装置包括:
用于在交叉拖缆维度中用零道以窗口中道数量的L倍填充窗口的装置;以及
用于在时间维度中用零采样以道中采样数量的L倍填充窗口的装置。
17.如权利要求14所述的系统,其中,用于用零采样周期性地置零第一扩展记录的装置包括:
用于周期性地在每两个实况道之间置零L个道以产生第二扩展记录的装置。
18.如权利要求14所述的系统,其中,用于缩小第一和第二频谱的装置包括:
用于将第一和第二频谱的频率范围缩小到最内部的采样,其频率范围大小等于经过变换的垂直质点振速记录中的频率数量的装置。
19.如权利要求14所述的系统,其中,用于取第一频谱和第二频谱的比率的装置包括:
用于在除之前将少量白噪声加到第二频谱的装置。
20.如权利要求14所述的系统,其中,应用一系列的缩放滤波器的装置包括:
用于将第三板分为编号为-L到L的2(L+1)条的装置;
用于将每两个对应编号的条组合为掩模的装置;
用于将缩放滤波器分为编号为-L到L的2(L+1)条的装置;
用于将每两个对应编号的条组合为滤波器的装置;
用于规格化掩模的装置;
用于将对应的滤波器和掩模乘在一起,生成一系列缩放滤波器的装置;以及
用于将缩放滤波器系列应用到经过变换的垂直质点振速记录的装置。
21.如权利要求20所述的系统,其中,每个条在交叉拖缆波数维度中包含的采样数量是经过变换的垂直质点振速记录中的一半。
22.如权利要求20所述的系统,其中,用于将每两个对应编号的条组合为滤波器的装置包括:
用于将具有正L的条移到正波数范围的装置;
用于将具有负L的条移到负波数范围的装置;以及
用于如果L是奇数,交换条的位置的装置。
23.如权利要求20所述的系统,其中,用于规格化掩模的装置包括:
用于规格化相同频率波数采样处在所有掩模上的采样的所有矢量的装置。
24.如权利要求20所述的系统,其中,用于组合经过变换的压力记录和经过缩放的垂直质点振速记录以分离上行和下行波场分量的装置包括:
用于对所有经过滤波的经过变换的垂直质点振速记录求总和,生成窗口中的经过缩放的垂直质点振速频谱的装置;
用于组合窗口中的压力频谱和经过缩放的垂直质点振速频谱,产生窗口中的上行和下行波场分量的频谱的装置;
用于将上行和下行波场分量的频谱从频率波数域逆变换到主测线波数域的装置;
用于将所有窗口的上行和下行波场分量的经过逆变换的频谱组合为恒定kx的切片,生成切片中的上行和下行波场分量的装置,其中kx是主测线波数;
用于将所有切片中的上行和下行波场分量组合为块,生成块中的上行和下行波场分量的装置;
用于将块中的上行和下行波场分量逆变换到时空域的装置;以及:
用于组合所有块中的上行和下行波场分量,生成分离的上行和下行波场分量的装置。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/009,440 US7646672B2 (en) | 2008-01-18 | 2008-01-18 | Method for wavefield separation in 3D dual sensor towed streamer data with aliased energy in cross-streamer direction |
US12/009440 | 2008-01-18 | ||
US12/009,440 | 2008-01-18 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101487899A CN101487899A (zh) | 2009-07-22 |
CN101487899B true CN101487899B (zh) | 2013-08-21 |
Family
ID=40671432
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910002704.1A Expired - Fee Related CN101487899B (zh) | 2008-01-18 | 2009-01-19 | 三维双传感器拖缆数据中的波场分离方法 |
Country Status (10)
Country | Link |
---|---|
US (1) | US7646672B2 (zh) |
EP (1) | EP2081054A3 (zh) |
CN (1) | CN101487899B (zh) |
BR (1) | BRPI0903062B1 (zh) |
CA (1) | CA2648682C (zh) |
EA (1) | EA014279B1 (zh) |
EG (1) | EG24976A (zh) |
MX (1) | MX2009000695A (zh) |
MY (1) | MY147535A (zh) |
SG (1) | SG154388A1 (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8582397B2 (en) * | 2009-01-06 | 2013-11-12 | Therataxis, Llc | Creating, directing and steering regions of intensity of wave propagation in inhomogeneous media |
US8239135B2 (en) * | 2009-05-07 | 2012-08-07 | Pgs Geophysical As | Method for calculation of seismic attributes from seismic signals |
US8208342B2 (en) * | 2009-09-14 | 2012-06-26 | Pgs Geophysical As | Method for combining signals of pressure and particle motion sensors in marine seismic streamers |
US8902699B2 (en) * | 2010-03-30 | 2014-12-02 | Pgs Geophysical As | Method for separating up and down propagating pressure and vertical velocity fields from pressure and three-axial motion sensors in towed streamers |
CN101893720B (zh) * | 2010-07-02 | 2012-09-05 | 中国科学院地质与地球物理研究所 | 一种地震波的矢量波场分离与合成的方法和系统 |
US8456950B2 (en) | 2010-07-30 | 2013-06-04 | Pgs Geophysical As | Method for wave decomposition using multi-component motion sensors |
WO2013019609A1 (en) | 2011-07-29 | 2013-02-07 | Hexagon Metrology, Inc. | Coordinate measuring system data reduction |
US9405027B2 (en) * | 2012-01-12 | 2016-08-02 | Westerngeco L.L.C. | Attentuating noise acquired in an energy measurement |
US9423518B2 (en) | 2012-02-09 | 2016-08-23 | Pgs Geophysical As | Method for processing dual-sensor streamer data with anti-alias protection |
US10459097B2 (en) * | 2012-11-19 | 2019-10-29 | Pgs Geophysical As | Methods and systems for extrapolating wavefields |
US10338250B2 (en) | 2013-03-14 | 2019-07-02 | Pgs Geophysical As | Method of removing incoherent noise |
US9322944B2 (en) | 2013-03-15 | 2016-04-26 | Pgs Geophysical As | Wavefield regularization by 3-D wavefield decomposition for geophysical data |
US20150276955A1 (en) * | 2013-11-06 | 2015-10-01 | Robert H. Brune | Method and System for Extending Spatial Wavenumber Spectrum Of Seismic Wavefields On Land Or Water Bottom Using Rotational Motion |
US9791580B2 (en) | 2014-04-17 | 2017-10-17 | Pgs Geophysical As | Methods and systems to separate wavefields using pressure wavefield data |
CN110554383B (zh) * | 2019-09-04 | 2021-04-06 | 中国科学院电子学研究所 | 用于微波频段的mimo环形阵列方位向成像方法及装置 |
US11573346B2 (en) * | 2021-04-15 | 2023-02-07 | Saudi Arabian Oil Company | Determining a seismic quality factor for subsurface formations for marine vertical seismic profiles |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5995904A (en) * | 1996-06-13 | 1999-11-30 | Exxon Production Research Company | Method for frequency domain seismic data processing on a massively parallel computer |
CN1404582A (zh) * | 2000-01-21 | 2003-03-19 | 施鲁博格控股有限公司 | 用于地震波场分离的系统和方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4486865A (en) | 1980-09-02 | 1984-12-04 | Mobil Oil Corporation | Pressure and velocity detectors for seismic exploration |
US4594693A (en) | 1983-11-04 | 1986-06-10 | Mobil Oil Corporation | Seismic trace interpolation using f-k filtering |
US5621699A (en) | 1995-07-07 | 1997-04-15 | Pgs Ocean Bottom Seismic, Inc. | Apparatus and method of calibrating vertical particle velocity detector and pressure detector in a sea-floor cable with in-situ passive monitoring |
FR2743897B1 (fr) | 1996-01-23 | 1998-04-10 | Geophysique Cie Gle | Procede de prospection sismique marine au moyen d'un couple de capteurs hydrophone et geophone |
US5677892A (en) | 1996-08-14 | 1997-10-14 | Western Atlas International, Inc. | Unaliased spatial trace interpolation in the f-k domain |
US5617372A (en) | 1996-08-14 | 1997-04-01 | Western Atlas International, Inc. | Unaliased spatial trace interpolation in the f-k domain |
GB9906456D0 (en) | 1999-03-22 | 1999-05-12 | Geco Prakla Uk Ltd | Method and system for reducing effects of sea surface ghost contamination in seismic data |
AU777935B2 (en) | 1999-06-21 | 2004-11-04 | Westerngeco Seismic Holdings Limited | 3-D seismic trace extrapolation and interpolation |
US7835225B2 (en) * | 2006-10-11 | 2010-11-16 | Pgs Geophysical As | Method for attenuating particle motion sensor noise in dual sensor towed marine seismic streamers |
-
2008
- 2008-01-18 US US12/009,440 patent/US7646672B2/en active Active
- 2008-12-19 SG SG200809502-8A patent/SG154388A1/en unknown
-
2009
- 2009-01-05 MY MYPI20090028A patent/MY147535A/en unknown
- 2009-01-05 EG EG2009010017A patent/EG24976A/xx active
- 2009-01-09 CA CA2648682A patent/CA2648682C/en not_active Expired - Fee Related
- 2009-01-16 EA EA200900030A patent/EA014279B1/ru not_active IP Right Cessation
- 2009-01-16 EP EP09150793A patent/EP2081054A3/en not_active Withdrawn
- 2009-01-16 MX MX2009000695A patent/MX2009000695A/es active IP Right Grant
- 2009-01-19 BR BRPI0903062-0 patent/BRPI0903062B1/pt not_active IP Right Cessation
- 2009-01-19 CN CN200910002704.1A patent/CN101487899B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5995904A (en) * | 1996-06-13 | 1999-11-30 | Exxon Production Research Company | Method for frequency domain seismic data processing on a massively parallel computer |
CN1404582A (zh) * | 2000-01-21 | 2003-03-19 | 施鲁博格控股有限公司 | 用于地震波场分离的系统和方法 |
Also Published As
Publication number | Publication date |
---|---|
BRPI0903062B1 (pt) | 2019-12-03 |
MY147535A (en) | 2012-12-31 |
CA2648682A1 (en) | 2009-07-18 |
US7646672B2 (en) | 2010-01-12 |
EA200900030A1 (ru) | 2009-08-28 |
SG154388A1 (en) | 2009-08-28 |
CN101487899A (zh) | 2009-07-22 |
BRPI0903062A2 (pt) | 2010-11-23 |
AU2008261111A1 (en) | 2009-08-06 |
EG24976A (en) | 2011-03-28 |
EP2081054A2 (en) | 2009-07-22 |
MX2009000695A (es) | 2009-08-12 |
EA014279B1 (ru) | 2010-10-29 |
EP2081054A3 (en) | 2011-03-02 |
CA2648682C (en) | 2014-09-30 |
US20090185444A1 (en) | 2009-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101487899B (zh) | 三维双传感器拖缆数据中的波场分离方法 | |
CN110462445B (zh) | 地球物理深度学习 | |
AU2018201492B2 (en) | Method and system for separating seismic sources in marine simultaneous shooting acquisition | |
EP2420865B1 (en) | Method for separating seismic sources in marine seismic surveys | |
US9151856B2 (en) | Separating interfering signals in seismic data | |
AU2012345565B2 (en) | Separation of simultaneous source data | |
US9229123B2 (en) | Method for handling rough sea and irregular recording conditions in multi-sensor towed streamer data | |
CN101539634A (zh) | 通过防混淆、防泄漏傅立叶变换内插地震数据的方法 | |
EP2360495B1 (en) | DIP-based corrections for data reconstruction in three-dimensional surface-related multiple prediction | |
US8315124B2 (en) | System and method for suppression of seismic multiple reflection signals | |
US20150066374A1 (en) | Seismic data processing with frequency diverse de-aliasing filtering | |
US20220342103A1 (en) | Noise Attenuation Methods Applied During Simultaneous Source Deblending and Separation | |
US11255992B2 (en) | Deblending method using patterned acquisition seismic data | |
WO2023076779A1 (en) | Separation of seismic sources by joint interpolation and deblending |
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: 20130821 Termination date: 20210119 |
|
CF01 | Termination of patent right due to non-payment of annual fee |