CN111856580B - 一种超远偏移距obs数据初至波能量增强方法及处理终端 - Google Patents
一种超远偏移距obs数据初至波能量增强方法及处理终端 Download PDFInfo
- Publication number
- CN111856580B CN111856580B CN202010735221.9A CN202010735221A CN111856580B CN 111856580 B CN111856580 B CN 111856580B CN 202010735221 A CN202010735221 A CN 202010735221A CN 111856580 B CN111856580 B CN 111856580B
- Authority
- CN
- China
- Prior art keywords
- data
- component
- obs
- arrival
- wave
- 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 24
- 238000012545 processing Methods 0.000 title claims abstract description 18
- 230000002708 enhancing effect Effects 0.000 claims abstract description 11
- 239000013535 sea water Substances 0.000 claims description 7
- 238000007781 pre-processing Methods 0.000 claims description 5
- 238000009792 diffusion process Methods 0.000 claims description 3
- 238000012958 reprocessing Methods 0.000 abstract 1
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000000926 separation method Methods 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
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
- G01V1/305—Travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/41—Arrival times, e.g. of P or S wave or first break
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
本发明涉及一种超远偏移距OBS数据初至波能量增强方法及处理终端,所述方法包括如下步骤:步骤1:根据OBS数据确定OBS在海底的位置和炮点的位置,并计算偏移距,得到包含偏移距的OBS数据;步骤2:从OBS数据中获得P分量和Z分量;步骤3:拾取P分量上初至波的旅行时,基于初至旅行时预设时窗,分别截取一段P分量数据和Z分量数据;步骤4:计算P分量数据的自相关和P分量数据与Z分量数据的互相关;步骤5:根据P分量数据和Z分量数据按对应公式得到初至波能量增强后的数据;步骤6:对比P分量数据、Z分量数据和初至波能量增强后的数据,优化时窗大小进行重新处理。本发明有效提高拾取超远偏移距初至波旅行时的可信度。
Description
技术领域
本发明涉及地震资料处理技术领域,具体是一种超远偏移距OBS数据初至波能量增强方法及处理终端。
背景技术
在海洋深部地壳结构研究中,通过OBS(中文称之为海底地震仪)探测是其中最重要的手段之一,通过在海面处的震源激发、设置在海底的OBS接收的方式来记录海底舌部的反射和折射信号,根据反射和折射信号研究海洋深部地壳结构。为了能够更好记录到来自海洋深部壳幔结构处的反射和折射信号,要求炮点(即震源)和接收点(即OBS)的间距(即偏移距)非常大,可达几百公里。
在获得反射和折射信号进行室内分析时,通常是通过初至波旅行时来反演海洋深部地壳结构,而初至波旅行时的可信度和分辨率与拾取的初至旅行时的数量及分布范围息息相关。因此,拾取分布范围广的初至波旅行时由于其可信度和分辨率较高,反演的海洋深部地壳结构效果也更好。实际获得数据中,一般在60km(公里)以内,拾取的初至波旅行时比较可靠,而在超远偏移距(>60km)地方,拾取初至波旅行时误差较大,不利于反演。目前,通常偏移距>60km可以视为超远偏移距,但其他起点的偏移距(例如>70km)也可以视为超远偏移距,本发明的超远偏移距是指偏移距>60km。
为了获得分布范围更广(也即偏移距更远)的初至波旅行时,也即需要增强初至波能量,以使得初至波能够在较远偏移距处也能够被拾取到,现有技术中,通常考虑在有限的中远偏移距范围内来增加射线数。例如,万奎元等提出利用在沉积层中经过多次反射的二次Pg折射波来增加反演的分辨率;Bharadwaj提出利用地震干涉方法把折射的一阶多次波转换为初至折射波;申请号为201810648783.2的中国发明专利提出基于镜像原理把折射多次波转换为虚拟OBS点上的初至波。在这些现有方案中,均未考虑超远偏移距处的初至波。
发明内容
针对现有技术的不足,本发明的目的之一是提供一种超远偏移距OBS数据初至波能量增强方法,其能够解决超远偏移距处初至波能量增强问题;
本发明的目的之二是提供一种处理终端,其能够解决超远偏移距处初至波能量增强问题。
实现本发明的目的之一的技术方案为:一种超远偏移距OBS数据初至波能量增强方法,包括如下步骤:
步骤1:获得OBS数据,根据OBS数据确定OBS在海底的位置以及炮点的位置,并计算偏移距,得到包括偏移距的OBS数据;
步骤2:从包括偏移距的OBS数据中获得P分量和Z分量;
步骤3:拾取P分量上初至波的旅行时,基于初至波旅行时预设一个时窗,通过时窗分别在P分量上初至波的旅行时和Z分量上初至波的旅行时各截取一段P分量数据和Z分量数据;
步骤4:按公式①、②分别计算P分量数据的自相关Rz(τ)和P分量数据与Z分量数据的互相关Rpz(τ):
式中,x(t)表示通过时窗截取到的Z分量数据,y(t)表示通过时窗截取到的P分类数据,T表示时窗宽度;
根据公式③计算标定算子cal:
步骤5:按公式④计算得到初至波能量增强的数据U:
式中,P和Z分别表示P分量数据和Z分量数据,ρ表示海水密度,c表示海水声速。
进一步地,所述根据OBS数据确定OBS在海底的位置,其实现过程包括如下步骤:
从OBS数据中拾取直达波的旅行时,通过最小二乘法对直达波的旅行时进行处理得到OBS在海底的位置。
进一步地,所述步骤2中,获得P分量和Z分量后,还包括分别对P分量和Z分量进行预处理,预处理包括噪音衰减、球面扩散和反褶积。
进一步地,所述步骤3中,时窗位于P分量和Z分量上初至波的旅行时上初至信号特征对应较好的偏移距处,截取得到对应的P分量数据和Z分量数据。
进一步地,在执行完步骤5之后,还包括以下步骤:
步骤6:对比P分量数据、Z分量数据和初至波能量增强后的数据U,以确定是否优化步骤3中的时窗大小,若是,重新选择时窗大小后,重复步骤3-步骤6,若否,则结束处理,以使得最终的初至波U合格。
实现本发明的目的之二的技术方案为:一种处理终端,其包括,存储器,用于存储程序指令;
处理器,用于运行所述程序指令,以执行所述超远偏移距OBS数据初至波能量增强方法的步骤。
本发明的有益效果为:本发明通过P分量数据和Z分量数据的组合来增强OBS数据中超远偏移距处初至波的能量,并且也间接提高了近偏移距处初至波的信噪比,能够提高拾取初至波旅行时的可信度。
附图说明
图1为本发明较佳实施例的流程图;
图2为拾取初至波的旅行时的示意图;
图3是某实际工区中OBS数据的P分量数据;
图4是同一工区中OBS数据的Z分量数据;
图5为对与图3和图4同一OBS数据按本发明处理得到的初至波;
图6为处理终端的结构示意图。
具体实施方案
下面,结合附图以及具体实施方案,对本发明做进一步描述:
如图1-图5所示,一种超远偏移距OBS数据初至波能量增强方法,包括如下步骤:
步骤1:获得OBS数据,以确定OBS在海底的位置,根据OBS海底的位置以及炮点的位置计算偏移距,得到包括偏移距的OBS数据。本步骤中,具体从OBS数据中拾取直达波的旅行时,通过最小二乘法对直达波的旅行时进行处理得到OBS在海底的位置。
步骤2:从包括偏移距的OBS数据中获得P分量和Z分量,分别对P分量和Z分量进行预处理,预处理包括噪音衰减、球面扩散和反褶积等。
采集到的OBS数据中,一般会记录4个分量,分别是水听器(水检检波器)分量、两个水平分量以及一个垂直分量,垂直分量通过陆检检波器接收得到。水听器分量即是P分量,两个水平分量分别是X分量和Y分量,垂直分量为Z分量,其中,P分量和Z分量记录了较多的初至波。
步骤3:拾取P分量上初至波的旅行时,对于初至波不清晰的地方,根据经验或已有信息进行大致延伸得到。基于初至波的旅行时预设一个时窗,运用该时窗分别在P分量上和Z分量上各截取一段P分量数据和Z分量数据。其中,时窗位于P分量和Z分量上初至信号特征较好的偏移距处,通过在初至信号特征较好的偏移距处设置时窗,从而截取对应的P分量数据和Z分量数据。
如图2所示,图2为拾取初至波的旅行时的示意图,图中左边的矩形框表示时窗,该时窗位于初至信号特征对应较好的偏移距处,右边的矩形框表示初至波不清晰的区域,对于此处不清晰的初至波,可根据经验或已有信息进行大致延伸得到。
步骤4:由于P分量和Z分量分别为水检检波器和陆检检波器接收的地震信号,两种检波器的仪器响应不同、与海底耦合关系不同,因此,在组合之前需要对P分量和Z分量数据进行标定处理,以使得两者分量数据在振幅、频率和相位特征一致。
按公式①、②分别计算P分量数据的自相关Rz(τ)和计算P分量数据与Z分量数据的互相关Rpz(τ):
式中,x(t)表示通过时窗截取到的Z分量数据,y(t)表示通过时窗截取到的P分类数据,T表示时窗宽度。
根据公式③计算标定算子cal:
标定算子cal为维纳滤波算子。
步骤5:按公式④计算得到初至波能量增强的数据U:
式中,P和Z分别表示P分量数据和Z分量数据,ρ表示海水密度,c表示海水声速,也即是声波在海水中的传播速度,ρ和c可通过海洋测绘得到,通常可从海洋测绘的CTD数据中获得。
数据U也即是由P分量数据和Z分量数据合并后的数据。
优选地,在执行完步骤5之后,还可以执行步骤6:
步骤6:对比P分量数据、Z分量数据和合并得到的数据U,以确定是否优化步骤3中的时窗大小,若是,重新选择时窗大小后,重复步骤3-步骤6,若否,则结束处理,以使得最终的初至波能量得到增强的数据U合格。
本发明通过P分量数据和Z分量数据的组合来增强OBS数据中超远偏移距处初至波的能量,并且也间接提高了近偏移距处初至波的信噪比,能够提高拾取初至波旅行时的可信度。
图3是某实际工区中OBS数据的P分量数据,图4是同一工区中OBS数据的Z分量数据。图3和图4可看出,OBS数据的最大偏移距超过150km,虽近偏移距处的初至波非常清晰,但在超远偏移距上(>60km)初至波的拾取比较困难,如位于偏移距100-160km处的方框。图5为对与图3和图4同一OBS数据按本发明处理得到的初至波能量增强的数据U,从图5可明显看出,在超远偏移距处,初至波的能量得到非常明显提升,如位于偏移距100-160km处的方框。
如图6所示,本发明还涉及一种实现以上方法的实体装置的处理终端100,其包括,
存储器101,用于存储程序指令;
处理器102,用于运行所述程序指令,以执行所述超远偏移距OBS数据初至波能量增强方法中的步骤。
本说明书所公开的实施例只是对本发明单方面特征的一个例证,本发明的保护范围不限于此实施例,其他任何功能等效的实施例均落入本发明的保护范围内。对于本领域的技术人员来说,可根据以上描述的技术方案以及构思,做出其它各种相应的改变以及变形,而所有的这些改变以及变形都应该属于本发明权利要求的保护范围之内。
Claims (6)
1.一种超远偏移距OBS数据初至波能量增强方法,其特征在于,包括如下步骤:
步骤1:获得OBS数据,根据OBS数据确定OBS在海底的位置以及炮点的位置,并计算偏移距,得到包括偏移距的OBS数据;
步骤2:从包括偏移距的OBS数据中获得P分量和Z分量;
步骤3:拾取P分量上初至波的旅行时,基于初至波的旅行时预设一个时窗,通过时窗分别在P分量上初至波的旅行时和Z分量上初至波的旅行时各截取一段P分量数据和Z分量数据;
步骤4:按公式①、②分别计算P分量数据的自相关Rz(τ)和P分量数据与Z分量数据的互相关Rpz(τ):
式中,x(t)表示通过时窗截取到的Z分量数据,y(t)表示通过时窗截取到的P分量 数据,T表示时窗宽度;
根据公式③计算标定算子cal:
步骤5:按公式④计算得到初至波能量增强的数据U:
式中,P和Z分别表示P分量数据和Z分量数据,ρ表示海水密度,c表示海水声速。
2.根据权利要求1所述的超远偏移距OBS数据初至波能量增强方法,其特征在于,所述根据OBS数据确定OBS在海底的位置,其实现过程包括如下步骤:
从OBS数据中拾取直达波的旅行时,通过最小二乘法对直达波的旅行时进行处理得到OBS在海底的位置。
3.根据权利要求1所述的超远偏移距OBS数据初至波能量增强方法,其特征在于,所述步骤2中,获得P分量和Z分量后,还包括分别对P分量和Z分量进行预处理,预处理包括噪音衰减、球面扩散和反褶积。
4.根据权利要求1所述的超远偏移距OBS数据初至波能量增强方法,其特征在于,所述步骤3中,时窗位于P分量和Z分量上初至波的旅行时上初至信号特征对应较好的偏移距处,截取得到对应的P分量数据和Z分量数据。
5.根据权利要求1所述的超远偏移距OBS数据初至波能量增强方法,其特征在于,在执行完步骤5之后,还包括以下步骤:
步骤6:对比P分量数据、Z分量数据和初至波能量增强后的U,以确定是否优化步骤3中的时窗大小,若是,重新选择时窗大小后,重复步骤3-步骤6,若否,则结束处理,以使得最终的初至波能量增强后的数据U合格。
6.一种处理终端,其包括,
存储器,用于存储程序指令;
处理器,用于运行所述程序指令,以执行如权利要求1-5任一项所述超远偏移距OBS数据初至波能量增强方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010735221.9A CN111856580B (zh) | 2020-07-27 | 2020-07-27 | 一种超远偏移距obs数据初至波能量增强方法及处理终端 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010735221.9A CN111856580B (zh) | 2020-07-27 | 2020-07-27 | 一种超远偏移距obs数据初至波能量增强方法及处理终端 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111856580A CN111856580A (zh) | 2020-10-30 |
CN111856580B true CN111856580B (zh) | 2022-02-15 |
Family
ID=72947737
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010735221.9A Expired - Fee Related CN111856580B (zh) | 2020-07-27 | 2020-07-27 | 一种超远偏移距obs数据初至波能量增强方法及处理终端 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111856580B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101545981A (zh) * | 2008-03-28 | 2009-09-30 | 中国石油天然气集团公司 | 可控震源地震数据零相位子波最小相位化方法 |
CN103675910A (zh) * | 2013-11-29 | 2014-03-26 | 中国石油天然气集团公司 | 一种水陆检波器地震数据标定因子反演方法 |
US20140200816A1 (en) * | 2013-01-14 | 2014-07-17 | Westerngeco L.L.C. | Seismic data processing |
CN107678063A (zh) * | 2017-09-25 | 2018-02-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于等级相关分析的多分量转换波裂缝预测方法 |
CN109696705A (zh) * | 2018-12-18 | 2019-04-30 | 中国石油天然气集团有限公司 | 深海obs上下行波场成像数据校正方法及装置 |
CN110687606A (zh) * | 2019-10-25 | 2020-01-14 | 长安大学 | 一种海底节点地震仪三分量定向校正方法 |
CN111273349A (zh) * | 2020-02-19 | 2020-06-12 | 广州海洋地质调查局 | 一种用于海底浅部沉积层的横波速度提取方法及处理终端 |
CN111352160A (zh) * | 2020-03-19 | 2020-06-30 | 中国科学院地质与地球物理研究所 | 一种海底地震仪自动重定位装置及方法 |
-
2020
- 2020-07-27 CN CN202010735221.9A patent/CN111856580B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101545981A (zh) * | 2008-03-28 | 2009-09-30 | 中国石油天然气集团公司 | 可控震源地震数据零相位子波最小相位化方法 |
US20140200816A1 (en) * | 2013-01-14 | 2014-07-17 | Westerngeco L.L.C. | Seismic data processing |
CN103675910A (zh) * | 2013-11-29 | 2014-03-26 | 中国石油天然气集团公司 | 一种水陆检波器地震数据标定因子反演方法 |
CN107678063A (zh) * | 2017-09-25 | 2018-02-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于等级相关分析的多分量转换波裂缝预测方法 |
CN109696705A (zh) * | 2018-12-18 | 2019-04-30 | 中国石油天然气集团有限公司 | 深海obs上下行波场成像数据校正方法及装置 |
CN110687606A (zh) * | 2019-10-25 | 2020-01-14 | 长安大学 | 一种海底节点地震仪三分量定向校正方法 |
CN111273349A (zh) * | 2020-02-19 | 2020-06-12 | 广州海洋地质调查局 | 一种用于海底浅部沉积层的横波速度提取方法及处理终端 |
CN111352160A (zh) * | 2020-03-19 | 2020-06-30 | 中国科学院地质与地球物理研究所 | 一种海底地震仪自动重定位装置及方法 |
Non-Patent Citations (1)
Title |
---|
南海北部陆坡崎岖海底区地震成像:OBS旅行时反演;刘斌;《物探与化探》;20161215;第40卷(第06期);第1244-1249页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111856580A (zh) | 2020-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6738715B2 (en) | Method for attenuating noise in seismic data | |
RU2319982C2 (ru) | Определение характеристик наклонных трещин с помощью обменных отраженных волн в сочетании с процессом последовательного исключения влияния вышележащих слоев | |
AU2002310037B2 (en) | Multiple suppression for ocean bottom seismic data | |
AU2009213054B2 (en) | Method for combining pressure and motion seismic signals from streamers where sensors are not at a common depth | |
EP2286277B1 (en) | Method for attenuating low frequency noise in a dual-sensor seismic streamer | |
CA2558000C (en) | Method, media, and signals for processing seismic data to obtain a velocity field | |
US5696734A (en) | Method and system for eliminating ghost reflections from ocean bottom cable seismic survey signals | |
AU2002324926A1 (en) | Method for attenuating noise in seismic data | |
US7069150B2 (en) | Method for optimizing migration fields using time slice analysis | |
AU2002310037A1 (en) | Multiple suppression for ocean bottom seismic data | |
US7110323B2 (en) | Method, system and apparatus for interpreting seismic data using duplex waves | |
AU2004318849B2 (en) | Generalized 3D surface multiple prediction | |
CA2497296C (en) | Removal of noise from seismic data using improved radon transformations | |
US8417458B2 (en) | Removing ground-roll from geophysical data | |
CN111856580B (zh) | 一种超远偏移距obs数据初至波能量增强方法及处理终端 | |
CN113514889B (zh) | 一种提升海洋深反射地震数据中低频信号能量的处理方法 | |
CN112415593B (zh) | 一种检波方法 | |
CN112824939B (zh) | 地震数据的剩余静校正方法、装置、设备及可读存储介质 | |
US11435492B2 (en) | Method for improving seismic acquisitions utilising active ultralight seismic detection systems | |
CN114859409B (zh) | 一种大洋岩石圈间断面信息的获取方法及装置 | |
Nasralla | Distributed Acoustic Sensing (DAS) As a New Tool For Subsurface Imaging | |
Guo et al. | Using Wide-Angle Reflection Wave Technology to Detect Seismic Inversion Data of Complex Geological Structure Zone for a Future Smart World | |
WO2022164700A1 (en) | Refraction seismic system and method for generating a velocity model | |
CN117991382A (zh) | 一种远探测声波测井的伪反射识别和分离方法及*** | |
CN117452500A (zh) | 一种基于反射波的近地表结构反演及静校正方法及*** |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220215 |