CN106199699B - 利用传输矩阵法清除鬼波的方法 - Google Patents
利用传输矩阵法清除鬼波的方法 Download PDFInfo
- Publication number
- CN106199699B CN106199699B CN201610509512.XA CN201610509512A CN106199699B CN 106199699 B CN106199699 B CN 106199699B CN 201610509512 A CN201610509512 A CN 201610509512A CN 106199699 B CN106199699 B CN 106199699B
- Authority
- CN
- China
- Prior art keywords
- ghost reflection
- wave detector
- ghost
- reflection
- channel
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 239000011159 matrix material Substances 0.000 title claims abstract description 37
- 238000012546 transfer Methods 0.000 title claims abstract description 15
- 230000005540 biological transmission Effects 0.000 claims abstract description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 12
- 230000005284 excitation Effects 0.000 claims description 7
- 239000013535 sea water Substances 0.000 claims description 6
- 230000007935 neutral effect Effects 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 2
- 229910017435 S2 In Inorganic materials 0.000 claims 1
- 238000012545 processing Methods 0.000 abstract description 5
- 238000001028 reflection method Methods 0.000 abstract description 5
- 230000000694 effects Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000002706 plastid Anatomy 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
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
本发明提出了一种利用传输矩阵法清除鬼波的方法,包括:获取每个检波器所接收到的单道地震数据;建立鬼波传输矩阵方程:Gm=d,利用G矩阵的广义逆求解出不含鬼波的模型数据向量m;m为不含鬼波的模型数据向量;D为包含鬼波的观测数据向量;对于每个检波器所接收到的单道地震数据,重复所述步骤S2以清除单道地震数据中的鬼波噪声,直到将所有检波器接收到的地震道数据处理完毕,最后输出不含鬼波噪声的地震数据。本发明有效避免处理中没有震源子波的问题,去鬼波方法简单,得到的结果分辨率高,频带变宽。
Description
技术领域
本发明涉及地震勘探中反射地震资料处理技术领域,特别涉及一种利用传输矩阵法清除鬼波的方法。
背景技术
鬼波是海上地震勘探中一种非常强的干扰波,它紧紧跟在一次波之后传播,也会被拖缆中的检波器接收到。海上拖缆地震采集一般是将震源和检波器沉放到水下一定深度,当震源激发出子波后,地震波在向下传播(简称下行波)的同时也向上传播(简称上行波)。由于海水面是一个很强的波阻抗界面,当上行波到达海面会产生反射,再向下传播;同理,由地下反射回来的地震波,有的直接到达检波器,有的继续向上传播,经海面反射后到达检波器,这种海面反射波称为虚反射(也称鬼波)。
从地震记录上看,鬼波同相轴与一次波同相轴混合在一起,很难把它们区分开,这对后续的解释工作是不利的。另外,由于鬼波的影响,地震资料的频谱上周期性出现了一系列陷波点,限制了有效频带的宽度,损失了一部分有意义的低频信息和高频信息,特别是一部分对于较深地质体成像非常有意义的低频信息。
目前压制鬼波的主要方法分为两种:
(1)从采集着手,例如:海底电缆采集(OBC)、上下双电缆采集、双传感器采集、斜拖缆采集等,配合对应的处理技术能够在一定程度上压制鬼波。通过采集的方法去除鬼波噪声费用昂贵,后期的处理过程也很复杂。
(2)基于子波反褶积处理的鬼波压制方法,在采集过程中记录下有关震源子波。由于偏移距变化和勘探环境变化,而且同一工区用采集到的唯一震源子波去除鬼波噪声,效果也会打折扣。
发明内容
本发明的目的旨在至少解决所述技术缺陷之一。
为此,本发明的目的在于提出一种利用传输矩阵法清除鬼波的方法,有效避免处理中没有震源子波的问题,去鬼波方法简单,得到的结果分辨率高,频带变宽。
为了实现上述目的,本发明的实施例提供一种利用传输矩阵法清除鬼波的方法,包括如下步骤:
步骤S1,获取每个检波器所接收到的单道地震数据;
步骤S2,建立鬼波传输矩阵方程:Gm=d,利用G矩阵的广义逆求解出不含鬼波的模型数据向量m,
其中,G为鬼波的传输矩阵,G矩阵列向量中1到-1间的零元素由鬼波的延迟时间决定,按照震源和检波器的深度得到鬼波的延迟时间,即震源和检波器所在水下的深度除以海水的速度;m为不含鬼波的模型数据向量;d为包含鬼波的观测数据向量;
步骤S3,对于每个检波器所接收到的单道地震数据,重复所述步骤S2以清除单道地震数据中的鬼波噪声,直到将所有检波器接收到的地震道数据处理完毕,最后输出不含鬼波的模型数据向量m。
进一步,在所述步骤S2中,
所述鬼波传输矩阵方程:Gm=d,表示为:
进一步,在所述步骤S2中,所述利用G矩阵的广义逆求解出不含鬼波的模型数据向量m,表示如下:
m=[GTG]-1GTd。
进一步,在所述步骤S1中,每个检波器所接收到的单道地震数据包括:激发点震源所在水下深度和接收点检波器所在水下深度。
进一步,在所述步骤S3中,对于每个检波器所接收到的单道地震数据,重复所述步骤S2以清除单道地震数据中的鬼波噪声,包括如下步骤:
首先,根据激发点震源所在水下深度重复所述步骤S2,去除震源所产生的鬼波噪声;
然后,根据接收点检波器所在水下深度重复所述步骤S2,去除接收点鬼波噪声,直到将所有检波器接收到的地震道数据处理完毕;
最后,通过上述对激发点震源和接收点检波器中鬼波噪声的清除,输出不含鬼波的模型数据向量。
根据本发明实施例的利用传输矩阵法清除鬼波的方法,只需知道海上地震勘探的震源和检波器的深度,即可对地震资料数据包含的鬼波进行处理,处理结果客观真实。本发明通过采用海洋地震数据鬼波的矩阵类似描述方法,通过过广义逆反演方法的数学运算就能去除地震数据中的鬼波,有效避免处理中没有震源子波的问题,去鬼波方法简单,得到的结果分辨率高,频带变宽。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明实施例的利用传输矩阵法清除鬼波的方法的流程图;
图2(a)至(c)分别为不包含鬼波的模型、包含鬼波的模型和清除鬼波后的模型示意图;
图3(a)和(b)分别为原始炮数据和清除鬼波后的数据示意图。
具体实施方式
下面详细描述本发明的实施例,实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
如图1所示,本发明实施例的利用传输矩阵法清除鬼波的方法,包括如下步骤:
步骤S1,获取每个检波器所接收到的单道地震数据。
在本发明的一个实施例中,每个检波器所接收到的单道地震数据包括:激发点震源所在水下深度和接收点检波器所在水下深度。
步骤S2,建立鬼波传输矩阵方程:Gm=d,表示为:
其中,G为鬼波的传输矩阵,G矩阵列向量中1到-1间的零元素由鬼波的延迟时间决定,按照震源和检波器的深度得到鬼波的延迟时间,即震源和检波器所在水下的深度除以海水的速度。其中,海水速度一般为每秒1500米,实际由地震勘探中测定此参数。例如震源水枪位于水下6米,海水速度1500米每秒,6米/1500米/1000毫秒=4毫秒,即1与-1之间有4个0。
m为不含鬼波的模型数据向量,d为包含鬼波的观测数据向量。
然后,利用G矩阵的广义逆求解出不含鬼波的模型数据向量m。
m=[GTG]-1GTd
步骤S3,对于每个检波器所接收到的单道地震数据,重复步骤S2以清除单道地震数据中的鬼波噪声,直到将所有检波器接收到的地震道数据处理完毕,最后输出不含鬼波噪声的地震数据m。
具体地,首先,根据激发点震源所在水下深度重复步骤S2,去除震源所产生的鬼波噪声;
然后,根据接收点检波器所在水下深度重复步骤S2,去除接收点鬼波噪声,直到将所有检波器接收到的地震道数据处理完毕;
最后,通过上述对激发点震源和接收点检波器中鬼波噪声的清除,输出不含鬼波噪声的地震数据。
下面参考图2和图3对本发明实施例的鬼波清除前后的效果进行说明。
为了测试广义反演去鬼波方法效果,本发明以单反射水平界面模型为例,设炮点检波点深6米,带通子波主频范围3-120赫兹。图2(a)为不含鬼波模型,图2(b)为含有深度为6米的炮点和检波点鬼波的模型,图2(c)为图2(b)去除炮点鬼波和检波点鬼波后显示结果。通过图2(a)至(c)可以很明显看到模型中加入鬼波后频带变窄,分辨率下降;而去除鬼波后模型的频带变宽。
图3(a)示出了原始炮数据,其中,炮点深6米,电缆深6米(检波点深6米)。图3(b)示出了图3(a)的原始数据去除鬼波后的记录,其中右上角是相关记录的频谱。通过图3(a)和(b)的频谱可以看出鬼波清除效果非常明显。
根据本发明实施例的利用传输矩阵法清除鬼波的方法,只需知道海上地震勘探的震源和检波器的深度,即可对地震资料数据包含的鬼波进行处理,处理结果客观真实。本发明通过采用海洋地震数据鬼波的矩阵类似描述方法,通过过广义逆反演方法的数学运算就能去除地震数据中的鬼波,有效避免处理中没有震源子波的问题,去鬼波方法简单,得到的结果分辨率高,频带变宽。
本发明实施例的利用传输矩阵法清除鬼波的方法,通过对实际资料的去噪处理,具有以下有益效果:
(1)采用近似的矩阵来表示鬼波,方法简单,压制鬼波效果好;
(2)该方法几乎可以被应用到任何类型的海洋勘探数据体。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求极其等同限定。
Claims (5)
1.一种利用传输矩阵法清除鬼波的方法,其特征在于,包括如下步骤:
步骤S1,获取每个检波器所接收到的单道地震数据;
步骤S2,建立鬼波传输矩阵方程:Gm=d,利用G矩阵的广义逆求解出不含鬼波的模型数据向量m,
其中,G为鬼波的传输矩阵,G矩阵列向量中1到-1间的零元素由鬼波的延迟时间决定,按照震源和检波器的深度得到鬼波的延迟时间,即震源和检波器所在水下的深度除以海水的速度;m为不含鬼波的模型数据向量;d为包含鬼波的观测数据向量;
步骤S3,对于每个检波器所接收到的单道地震数据,重复所述步骤S2以清除单道地震数据中的鬼波噪声,直到将所有检波器接收到的地震道数据处理完毕,最后输出不含鬼波的模型数据向量m。
2.如权利要求1所述用传输矩阵法清除鬼波的方法,其特征在于,在所述步骤S2中,所述鬼波传输矩阵方程:Gm=d,表示为:
3.如权利要求2所述的利用传输矩阵法清除鬼波的方法,其特征在于,在所述步骤S2中,所述利用G矩阵的广义逆求解出不含鬼波的模型数据向量m,表示如下:
m=[GTG]-1GTd。
4.如权利要求1所述的利用传输矩阵法清除鬼波的方法,其特征在于,在所述步骤S1中,每个检波器所接收到的单道地震数据包括:激发点震源所在水下深度和接收点检波器所在水下深度。
5.如权利要求3所述的利用传输矩阵法清除鬼波的方法,其特征在于,在所述步骤S3中,对于每个检波器所接收到的单道地震数据,重复所述步骤S2以清除单道地震数据中的鬼波噪声,包括如下步骤:
首先,根据激发点震源所在水下深度重复所述步骤S2,去除震源所产生的鬼波噪声;
然后,根据接收点检波器所在水下深度重复所述步骤S2,去除接收点鬼波噪声,直到将所有检波器接收到的地震道数据处理完毕;
最后,通过上述对激发点震源和接收点检波器中鬼波噪声的清除,输出不含鬼波的模型数据向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610509512.XA CN106199699B (zh) | 2016-06-30 | 2016-06-30 | 利用传输矩阵法清除鬼波的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610509512.XA CN106199699B (zh) | 2016-06-30 | 2016-06-30 | 利用传输矩阵法清除鬼波的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106199699A CN106199699A (zh) | 2016-12-07 |
CN106199699B true CN106199699B (zh) | 2018-07-20 |
Family
ID=57462780
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610509512.XA Active CN106199699B (zh) | 2016-06-30 | 2016-06-30 | 利用传输矩阵法清除鬼波的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106199699B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107450103B (zh) * | 2017-09-07 | 2020-04-28 | 中海石油(中国)有限公司 | 一种基于边界积分反算子的鬼波压制方法 |
CN110703332B (zh) * | 2018-07-09 | 2021-06-29 | 中国石油化工股份有限公司 | 一种鬼波压制方法 |
CN110967735A (zh) * | 2018-09-28 | 2020-04-07 | 中国石油化工股份有限公司 | 自适应的鬼波压制方法及*** |
CN109725352B (zh) * | 2018-11-21 | 2020-04-10 | 中国石油天然气集团有限公司 | 地震道的反射时间校正方法及*** |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8116166B2 (en) * | 2007-09-10 | 2012-02-14 | Westerngeco L.L.C. | 3D deghosting of multicomponent or over / under streamer recordings using cross-line wavenumber spectra of hydrophone data |
US8693282B2 (en) * | 2010-05-25 | 2014-04-08 | Westerngeco L.L.C. | Deghosting seismic data |
CN102169189B (zh) * | 2011-01-19 | 2015-01-14 | 中国海洋石油总公司 | 深水层间多次波消除方法 |
US9274239B2 (en) * | 2012-01-13 | 2016-03-01 | Westerngeco L.L.C. | Wavefield deghosting |
US9360576B2 (en) * | 2012-04-27 | 2016-06-07 | Pgs Geophysical As | Methods and apparatus for generating deghosted seismic data |
CN103308943B (zh) * | 2013-05-10 | 2015-11-18 | 中国石油天然气股份有限公司 | 一种海洋地震资料处理中层间多次波衰减的方法及装置 |
CN104360394B (zh) * | 2014-12-01 | 2017-02-22 | 中国海洋石油总公司 | 一种海上变深度缆地震数据鬼波压制的方法 |
-
2016
- 2016-06-30 CN CN201610509512.XA patent/CN106199699B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106199699A (zh) | 2016-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2141514B1 (en) | Method for attenuation of multiple reflections in seismic data | |
CN101907728B (zh) | 采集和处理海洋地震数据以提取和使用上行和下行波场 | |
US6654693B2 (en) | Angle dependent surface multiple attenuation for two-component marine bottom sensor data | |
CN106199699B (zh) | 利用传输矩阵法清除鬼波的方法 | |
US8855934B2 (en) | Noise attenuation of seismic data | |
US20120275267A1 (en) | Seismic Data Processing | |
EP3339910A1 (en) | Device and method for model-based deblending | |
MX2010012702A (es) | Determinacion conjunta de interpolacion y eliminacion de fantasmas de datos sismicos. | |
CN109143363B (zh) | 海洋拖缆双检采集鬼波压制方法及*** | |
US10073183B2 (en) | Methods and systems that attenuate noise in seismic data | |
MX2011008294A (es) | Formacion de imagenes con mediciones vectoriales. | |
WO2010082126A2 (en) | Processing seismic data | |
CN107884828A (zh) | 一种在空间‑频率域基于格林理论的鬼波压制方法 | |
US11994640B2 (en) | Attenuation of low-frequency noise in continuously recorded wavefields | |
US9658354B2 (en) | Seismic imaging systems and methods employing correlation-based stacking | |
CN105301656B (zh) | 一种压制虚反射信号的方法及装置 | |
CN113514889B (zh) | 一种提升海洋深反射地震数据中低频信号能量的处理方法 | |
MX2012005137A (es) | Sistema y tecnica para aumentar la separacion de los sensores del movimiento de particulas en su cable sismico. | |
US10338253B2 (en) | Method of suppressing spectral artefacts of wavefield decomposition caused by imperfect extrapolation | |
US11448791B2 (en) | Marine deterministic notch compensation | |
Zou et al. | Enhancing low-frequency water-column acoustic reflections in marine multichannel seismic data for seismic oceanography | |
CN107589456B (zh) | 一种获取地震数据的方法和装置及计算机可读存储介质 | |
CN105425297A (zh) | 一种压制虚反射信号的方法及装置 | |
US11287542B2 (en) | Seismic data processing including surface multiple modeling for ocean bottom cable systems | |
Teixeira et al. | Implementation of processing techniques in single-channel high-resolution seismic data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |