CN117369020A - 海洋拖缆地震勘探老化电缆检波器灵敏度校正方法 - Google Patents
海洋拖缆地震勘探老化电缆检波器灵敏度校正方法 Download PDFInfo
- Publication number
- CN117369020A CN117369020A CN202210778544.5A CN202210778544A CN117369020A CN 117369020 A CN117369020 A CN 117369020A CN 202210778544 A CN202210778544 A CN 202210778544A CN 117369020 A CN117369020 A CN 117369020A
- Authority
- CN
- China
- Prior art keywords
- cable
- seismic
- amplitude value
- sensitivity
- average amplitude
- 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.)
- Pending
Links
- 230000035945 sensitivity Effects 0.000 title claims abstract description 56
- 238000012937 correction Methods 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000002431 foraging effect Effects 0.000 title description 3
- 238000005070 sampling Methods 0.000 claims abstract description 72
- 230000002159 abnormal effect Effects 0.000 claims abstract description 39
- 238000012545 processing Methods 0.000 claims abstract description 6
- 230000032683 aging Effects 0.000 claims abstract description 5
- 230000005284 excitation Effects 0.000 claims abstract description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 3
- 230000002547 anomalous effect Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于海洋地震资料处理领域,具体公开了一种海洋拖缆地震勘探老化电缆检波器灵敏度校正方法,包括:电缆中检波器接收震源激发产生的地震波,生成海洋地震记录;选取时窗;确定异常电缆;计算异常电缆采样点的平均振幅值插值,替换异常电缆采样点的平均振幅值;根据标准电缆采样点的平均振幅值与标准电缆中检波器灵敏度的关系,以及步骤S4中得到的平均振幅值,计算得到异常电缆中检波器的校正灵敏度。本发明有针对性地对变化较大的电缆中检波器的灵敏度进行校正,提高了旧电缆的灵敏度,保证旧电缆的正常使用,且不会对符合要求的新电缆的灵敏度进行错误的校正。本发明适用于海洋拖缆地震勘探。
Description
技术领域
本发明属于海洋地震资料处理领域,具体地说是一种海洋拖缆地震勘探老化电缆检波器灵敏度校正方法。
背景技术
地震勘探是通过人工方法激发地震波,研究地震波在地层中传播的规律与特点,查明地下地质构造,为寻找油气田或其他勘探目标服务的一种物探方法,地震勘探的采集过程就是产生地震波与接收地震波的过程,海洋拖缆地震勘探,使用气枪作为产生地震波的震源,拖拽电缆中安装有特定组合的压电式检波器,压电式检波器将地震波的压力信号转换为电信号并被仪器设备记录下来以供后续处理解释。压电检波器将压力转换为电信号,根据内部电子元件特性,转换具有一比例常数即灵敏度,单位一般为伏每巴(v/bar),不同种类拖缆设备,灵敏度值不同。海洋拖缆地震勘探采集设备随着使用时间的增长,检波器内部元件性能变化,导致检波器灵敏度下降,最终出现不同使用时间的电缆灵敏度不同的现象。
发明内容
本发明的目的,是要提供一种海洋拖缆地震勘探老化电缆检波器灵敏度校正方法,以提高旧电缆的灵敏度,保证旧电缆的正常使用,且不会对符合要求的新电缆的灵敏度进行错误的校正。
本发明为实现上述目的,所采用的技术方法如下:
一种海洋拖缆地震勘探老化电缆检波器灵敏度校正方法,包括依次进行的以下步骤:
S1、电缆中检波器接收震源激发产生的地震波,生成海洋地震记录;每根电缆中包含多个地震通道,每个地震通道是由多组检波器线性组合构成的;
S2、将生成的海洋地震记录加载到地震观测系统中,选取时窗;
S3、针对时窗内每根电缆采样点的振幅进行数据处理,确定异常电缆;
S4、去掉异常电缆采样点的平均振幅值,保留正常电缆采样点的平均振幅值,计算异常电缆采样点的平均振幅值插值,作为估算的异常电缆采样点的平均振幅值,替换真实测量的异常电缆采样点的平均振幅值;
S5、根据标准电缆采样点的平均振幅值与标准电缆中检波器灵敏度的关系,以及步骤S4中得到的平均振幅值,计算得到异常电缆中检波器的校正灵敏度。
作为限定:步骤S1中震源采用气枪震源。
作为限定:步骤S2具体包括依次进行的以下步骤:
S21、将海洋地震记录加载到地震观测系统中,得到炮点水深值和偏移距;
S22、将多炮海洋地震记录中相同电缆中相同地震通道采样点接收的数据进行叠加,叠加成一炮海洋地震记录,通道数不变;
S23、计算地震波由震源激发到达海底反射后被电缆中检波器接收的旅行时间,计算公式为:
其中,D1为偏移距,单位为米,D2为震源到海底的深度,单位为米,V1为声波在水中传播速度,单位为米/秒,Tn为地震波由震源激发到达海底反射后被电缆中检波器接收的旅行时间,单位为秒;
选取地震波由震源激发到达海底反射后被电缆中检波器接收的旅行时间以下500毫秒时窗,截取时窗内的海洋地震记录。
作为限定:步骤S21中得到炮点水深值和偏移距后,选取连续多炮的海洋地震记录,根据标准电缆中检波器的灵敏度,转换地震通道采样点的振幅单位,将振幅单位由毫伏转换为微巴,其转换公式为:
T微巴=T毫伏÷S标准÷1000
其中,T微巴为采样点的振幅单位为微巴的振幅值,T毫伏为采样点的振幅单位为毫伏的振幅值,S标准为标准电缆中检波器的灵敏度,单位为伏/巴。
作为限定:步骤S3具体包括依次进行的以下步骤:
S31、时窗内对每根电缆中每个地震通道采样点振幅值进行均方根计算,再求取每根电缆所有地震通道采样点均方根振幅值的平均值,即为每根电缆采样点的平均振幅值;
S32、对比所有相邻电缆采样点平均振幅值的变化幅度,根据变化幅度的大小确定异常电缆,当变化幅度超过设定的阈值时,则此变化幅度对应的电缆为异常电缆。
作为限定:步骤S32中所有相邻电缆检波器生成的海洋地震记录中时窗内采样点平均振幅值的变化幅度计算公式为:
dB=20log10(a/b)
其中,a为相邻两个电缆中前一个电缆采样点的平均振幅值,b为相邻两个电缆中后一个电缆采样点的平均振幅值,dB为相邻电缆采样点平均振幅值的变化幅度。
作为限定:步骤S4具体为采用二维三次样条函数插值方法,计算异常电缆采样点的平均振幅值插值。
作为限定:步骤S5中异常电缆中检波器的校正灵敏度的计算公式为:
S校正=(T标准×S标准)÷T插值
其中,T标准×S标准=真实测量的异常电缆采样点的平均振幅值,S标准为标准电缆中检波器的灵敏度,单位为v/bar,T标准为标准电缆采样点的平均振幅值,单位为微巴, T插值为采用二维三次样条函数插值方法计算出的平均振幅值插值,单位为微巴,S校正为异常电缆中检波器的校正灵敏度,单位为伏每巴。
本发明由于采用了上述方案,与现有技术相比,所取得的有益效果是:
本发明提供的一种海洋拖缆地震勘探老化电缆检波器灵敏度校正方法,有针对性地对变化较大的电缆中检波器的灵敏度进行校正,提高了旧电缆的灵敏度,保证旧电缆的正常使用,且不会对符合要求的新电缆的灵敏度进行错误的校正;将振幅单位由毫伏转换为微巴,因为在灵敏度有失真的情况下,振幅值并非真实,而转换成的电压值则是真实记录值,因此,提高了灵敏度校正的准确性。
本发明适用于海洋拖缆地震勘探。
附图说明
下面结合附图及具体实施例对本发明作更进一步详细说明。
图1为本实施例检波器的示意图;
图2为本实施例每道地震通道的示意图;
图3为本实施例每根电缆的示意图;
图4为本实施例每根电缆采样点的平均振幅值示意图;
图5为本实施例去掉异常电缆采样点的平均振幅值保留正常电缆采样点的平均振幅值的示意图;
图6为本实施例异常电缆采样点的平均振幅值插值替换平均振幅值示意图;
图7为本实施例插值前后每根电缆采样点的平均振幅值对比图;
图8为本实施例异常电缆中检波器为进行灵敏度校正前的地震记录;
图9为本实施例异常电缆中检波器为进行灵敏度校正后的地震记录;
图10为本实施例异常电缆中检波器为进行灵敏度校正前和校正后的地震通道振幅差值。
具体实施方式
下面结合实施例对本发明作进一步说明,但本领域的技术人员应当理解,本发明并不限于以下实施例,任何在本发明具体实施例基础上做出的改进和等效变化,都在本发明权利要求保护的范围之内。
实施例一种海洋拖缆地震勘探老化电缆检波器灵敏度校正方法
一种海洋拖缆地震勘探老化电缆检波器灵敏度校正方法,包括依次进行的以下步骤:
S1、用气枪震源,释放气枪电缆,激发气枪,震源正常激发200炮,所有电缆中检波器接收震源激发产生的地震波,生成海洋地震记录;每根电缆中包含多个地震通道,每个地震通道是由多组检波器线性组合构成的,本实施例中一根电缆中包含12 道地震通道,检波器的示意图如图1所示,每道地震通道的示意图如图2所示,每根电缆的示意图如图3所示;
S2、将生成的海洋地震记录加载到地震观测系统中,选取时窗;具体包括依次进行的以下步骤:
S21、将海洋地震记录加载到地震观测系统中得到炮点水深值和偏移距,选取连续200炮的海洋地震记录,根据标准电缆中检波器的灵敏度,转换地震通道采样点的振幅单位,将振幅单位由毫伏转换为微巴,其转换公式为:
T微巴=T毫伏÷S标准÷1000
其中,T微巴为采样点的振幅单位为微巴的振幅值,T毫伏为采样点的振幅单位为毫伏的振幅值,S标准为标准电缆中检波器的灵敏度,单位为伏/巴;
S22、将200炮海洋地震记录中相同电缆中相同地震通道采样点接收的数据进行叠加,叠加成一炮海洋地震记录,通道数不变;
S23、计算地震波由震源激发到达海底反射后被电缆中检波器接收的旅行时间,计算公式为:
其中,D1为偏移距,单位为米,D2为震源到海底的深度,单位为米,V1为声波在水中传播速度,单位为米/秒,Tn为地震波由震源激发到达海底反射后被电缆中检波器接收的旅行时间,单位为秒;
选取地震波由震源激发到达海底反射后被电缆中检波器接收的旅行时间以下500毫秒时窗,截取时窗内的海洋地震记录;
S3、针对时窗内每根电缆采样点的振幅进行数据处理,确定异常电缆;具体包括依次进行的以下步骤:
S31、时窗内对每根电缆中每个地震通道采样点振幅值进行均方根计算,再求取每根电缆所有地震通道采样点均方根振幅值的平均值,即为每根电缆采样点的平均振幅值,如图4所示;
S32、对比所有相邻电缆采样点平均振幅值的变化幅度,计算公式为:
dB=20log10(a/b)
其中,a为相邻两个电缆中前一个电缆采样点的平均振幅值,b为相邻两个电缆中后一个电缆采样点的平均振幅值,dB为相邻电缆采样点平均振幅值的变化幅度;
根据变化幅度的大小确定异常电缆,当变化幅度超过设定的阈值时,则此变化幅度对应的电缆为异常电缆;
S4、去掉异常电缆采样点的平均振幅值,保留正常电缆采样点的平均振幅值,如图5所示,采用二维三次样条函数插值方法,计算异常电缆采样点的平均振幅值插值,异常电缆采样点的平均振幅值插值为估算的异常电缆采样点的平均振幅值,替换真实测量的异常电缆采样点的平均振幅值,如图6所示,插值前后每根电缆采样点的平均振幅值对比如图7所示;
S5、根据标准电缆采样点的平均振幅值与标准电缆中检波器灵敏度的关系,以及步骤S4中得到的平均振幅值,计算得到异常电缆中检波器的校正灵敏度,计算公式为:
S校正=(T标准×S标准)÷T插值
其中,T标准×S标准=真实测量的异常电缆采样点的平均振幅值,S标准为标准电缆中检波器的灵敏度,单位为v/bar,T标准为标准电缆采样点的平均振幅值,单位为微巴, T插值为采用二维三次样条函数插值方法计算出的平均振幅值插值,单位为微巴,S校正为异常电缆中检波器的校正灵敏度,单位为伏每巴。
本实施例中异常电缆中检波器进行灵敏度校正前的地震记录如图8所示,进行灵敏度校正后的地震记录如图9所示,异常电缆中检波器为进行灵敏度校正前和校正后的地震通道振幅差值如图10所示。
Claims (8)
1.一种海洋拖缆地震勘探老化电缆检波器灵敏度校正方法,其特征在于,包括依次进行的以下步骤:
S1、电缆中检波器接收震源激发产生的地震波,生成海洋地震记录;每根电缆中包含多个地震通道,每个地震通道是由多组检波器线性组合构成的;
S2、将生成的海洋地震记录加载到地震观测系统中,选取时窗;
S3、针对时窗内每根电缆采样点的振幅进行数据处理,确定异常电缆;
S4、去掉异常电缆采样点的平均振幅值,保留正常电缆采样点的平均振幅值,计算异常电缆采样点的平均振幅值插值,作为估算的异常电缆采样点的平均振幅值,替换真实测量的异常电缆采样点的平均振幅值;
S5、根据标准电缆采样点的平均振幅值与标准电缆中检波器灵敏度的关系,以及步骤S4中得到的平均振幅值,计算得到异常电缆中检波器的校正灵敏度。
2.根据权利要求1所述的海洋拖缆地震勘探老化电缆检波器灵敏度校正方法,其特征在于,步骤S1中震源采用气枪震源。
3.根据权利要求2所述的海洋拖缆地震勘探老化电缆检波器灵敏度校正方法,其特征在于,步骤S2具体包括依次进行的以下步骤:
S21、将海洋地震记录加载到地震观测系统中,得到炮点水深值和偏移距;
S22、将多炮海洋地震记录中相同电缆中相同地震通道采样点接收的数据进行叠加,叠加成一炮海洋地震记录,通道数不变;
S23、计算地震波由震源激发到达海底反射后被电缆中检波器接收的旅行时间,计算公式为:
其中,D1为偏移距,单位为米,D2为震源到海底的深度,单位为米,V1为声波在水中传播速度,单位为米/秒,Tn为地震波由震源激发到达海底反射后被电缆中检波器接收的旅行时间,单位为秒;
选取地震波由震源激发到达海底反射后被电缆中检波器接收的旅行时间以下500毫秒时窗,截取时窗内的海洋地震记录。
4.根据权利要求3所述的海洋拖缆地震勘探老化电缆检波器灵敏度校正方法,其特征在于,步骤S21中得到炮点水深值和偏移距后,选取连续多炮的海洋地震记录,根据标准电缆中检波器的灵敏度,转换地震通道采样点的振幅单位,将振幅单位由毫伏转换为微巴,其转换公式为:
T微巴=T毫伏÷S标准÷1000
其中,T微巴为采样点的振幅单位为微巴的振幅值,T毫伏为采样点的振幅单位为毫伏的振幅值,S标准为标准电缆中检波器的灵敏度,单位为伏/巴。
5.根据权利要求4所述的海洋拖缆地震勘探老化电缆检波器灵敏度校正方法,其特征在于,步骤S3具体包括依次进行的以下步骤:
S31、时窗内对每根电缆中每个地震通道采样点振幅值进行均方根计算,再求取每根电缆所有地震通道采样点均方根振幅值的平均值,即为每根电缆采样点的平均振幅值;
S32、对比所有相邻电缆采样点平均振幅值的变化幅度,根据变化幅度的大小确定异常电缆,当变化幅度超过设定的阈值时,则此变化幅度对应的电缆为异常电缆。
6.根据权利要求5所述的海洋拖缆地震勘探老化电缆检波器灵敏度校正方法,其特征在于,步骤S32中所有相邻电缆检波器生成的海洋地震记录中时窗内采样点平均振幅值的变化幅度计算公式为:
dB=20log10(a/b)
其中,a为相邻两个电缆中前一个电缆采样点的平均振幅值,b为相邻两个电缆中后一个电缆采样点的平均振幅值,dB为相邻电缆采样点平均振幅值的变化幅度。
7.根据权利要求6所述的海洋拖缆地震勘探老化电缆检波器灵敏度校正方法,其特征在于,步骤S4具体为采用二维三次样条函数插值方法,计算异常电缆采样点的平均振幅值插值。
8.根据权利要求6所述的海洋拖缆地震勘探老化电缆检波器灵敏度校正方法,其特征在于,步骤S5中异常电缆中检波器的校正灵敏度的计算公式为:
S校正=(T标准×S标准)÷T插值
其中,T标准×S标准=真实测量的异常电缆采样点的平均振幅值,S标准为标准电缆中检波器的灵敏度,单位为v/bar,T标准为标准电缆采样点的平均振幅值,单位为微巴,T插值为采用二维三次样条函数插值方法计算出的平均振幅值插值,单位为微巴,S校正为异常电缆中检波器的校正灵敏度,单位为伏每巴。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210778544.5A CN117369020A (zh) | 2022-06-30 | 2022-06-30 | 海洋拖缆地震勘探老化电缆检波器灵敏度校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210778544.5A CN117369020A (zh) | 2022-06-30 | 2022-06-30 | 海洋拖缆地震勘探老化电缆检波器灵敏度校正方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117369020A true CN117369020A (zh) | 2024-01-09 |
Family
ID=89402872
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210778544.5A Pending CN117369020A (zh) | 2022-06-30 | 2022-06-30 | 海洋拖缆地震勘探老化电缆检波器灵敏度校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117369020A (zh) |
-
2022
- 2022-06-30 CN CN202210778544.5A patent/CN117369020A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2016200908B2 (en) | Method for acquiring and processing marine seismic data to extract and constructively use the up-going and down-going wave-fields emitted by the source(s) | |
EP3324215B1 (en) | Method of first arrival picking of multi-channel seismic survey data | |
US8964502B2 (en) | Zero offset profile from near-field hydrophones | |
EP0101281A2 (en) | Underwater seismic testing | |
CN111487678B (zh) | 一种确定高分辨率小多道地震最小偏移距和系统延迟的分析方法 | |
GB2440276A (en) | Processing seismic data acquired using three or more over/under streamers and/or sources | |
CN101881836A (zh) | 用于根据地震信号来计算地震属性的方法 | |
US8639442B2 (en) | Identifying invalid seismic data | |
KR102003466B1 (ko) | 다중채널 3차원 지층탐사자료의 너울 영향 보정방법, 이를 이용한 해양 3d 탄성파 탐사방법 | |
MX2014011894A (es) | Metodos y sistemas para prospecciones sismicas terrestres. | |
CN112327362B (zh) | 速度域的海底多次波预测与追踪衰减方法 | |
US11442186B2 (en) | Source-receiver position estimation using direct arrival modeling and inversion | |
CN103454681A (zh) | 评价三维地震观测系统成像效果的方法和设备 | |
WO2010027682A2 (en) | Processing seismic data in common group-center gathers | |
KR100412097B1 (ko) | 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치 | |
CN117369020A (zh) | 海洋拖缆地震勘探老化电缆检波器灵敏度校正方法 | |
CN113075732B (zh) | 一种消除高分辨率小多道地震地层异常起伏的方法 | |
US11307326B2 (en) | Calibration of combined acceleration and pressure sensors | |
CN109425892B (zh) | 地震子波的估计方法及系统 | |
Fricke et al. | A standard quantitative calibration procedure for marine seismic sources | |
AU2013201068B2 (en) | Determination of particle displacement or velocity from particle acceleration measurements | |
CN109490962A (zh) | 一种浅地层剖面非零偏移距的消除方法 | |
US5212668A (en) | Method for quantification of the level of weather and sea noise during marine seismic surveys | |
KR101845965B1 (ko) | 탄성파탐사 자료의 너울 보정 방법 | |
US4813026A (en) | Method for logarithmic analysis of seismic reflection signals |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |