CN116931088B - 海洋高分辨率地震数据涌浪相关静校正量快速估算方法 - Google Patents
海洋高分辨率地震数据涌浪相关静校正量快速估算方法 Download PDFInfo
- Publication number
- CN116931088B CN116931088B CN202311188443.3A CN202311188443A CN116931088B CN 116931088 B CN116931088 B CN 116931088B CN 202311188443 A CN202311188443 A CN 202311188443A CN 116931088 B CN116931088 B CN 116931088B
- Authority
- CN
- China
- Prior art keywords
- submarine
- seismic
- recorded
- record
- prestack
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 42
- 230000003068 static effect Effects 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000012876 topography Methods 0.000 claims abstract description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 8
- 239000013535 sea water Substances 0.000 claims description 24
- 238000004364 calculation method Methods 0.000 claims description 23
- 102100030794 Conserved oligomeric Golgi complex subunit 1 Human genes 0.000 claims description 9
- 101000920124 Homo sapiens Conserved oligomeric Golgi complex subunit 1 Proteins 0.000 claims description 9
- 238000001514 detection method Methods 0.000 claims description 6
- 238000010892 electric spark Methods 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 5
- 102100030797 Conserved oligomeric Golgi complex subunit 2 Human genes 0.000 claims description 4
- 101000920113 Homo sapiens Conserved oligomeric Golgi complex subunit 2 Proteins 0.000 claims description 4
- 101001138030 Homo sapiens Protein Largen Proteins 0.000 claims description 3
- 102100020860 Protein Largen Human genes 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 abstract description 16
- 238000012545 processing Methods 0.000 abstract description 8
- 238000010276 construction Methods 0.000 abstract description 4
- 238000005457 optimization Methods 0.000 abstract description 3
- 238000001303 quality assessment method Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 102100029265 Conserved oligomeric Golgi complex subunit 3 Human genes 0.000 description 1
- 101000770432 Homo sapiens Conserved oligomeric Golgi complex subunit 3 Proteins 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- NMJORVOYSJLJGU-UHFFFAOYSA-N methane clathrate Chemical compound C.C.C.C.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O NMJORVOYSJLJGU-UHFFFAOYSA-N 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 229910052761 rare earth metal Inorganic materials 0.000 description 1
- 150000002910 rare earth metals Chemical class 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000000007 visual effect Effects 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
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及海洋地球物理勘探技术领域，提供了一种海洋高分辨率地震数据涌浪相关静校正量快速估算方法，包括：抽取共偏移距道集记录；得到各叠前炮集记录的海底反射波同相轴实际时距曲线；获取海底水深地形数据；根据所述海底水深地形数据和观测系统参数，计算各叠前炮集记录的海底反射波同相轴理论时距曲线；对所述海底反射波同相轴理论时距曲线与海底反射波同相轴实际时距曲线求差，获取叠前炮集记录各地震道的涌浪相关静校正量。本发明根据海底反射地震记录和观测系统参数对涌浪相关的静校正量进行快速估算，能够有效改善现场处理地震成像剖面的分辨率和信噪比，为外业地震资料现场质量评定和施工方案优化提供高效且可靠的数据支持。
Description
技术领域
本发明涉及海洋地球物理勘探技术领域，具体涉及海洋高分辨率地震数据涌浪相关静校正量快速估算方法。
背景技术
海洋高分辨率多道地震探测是指以电火花震源和小道距地震电缆为基础的高分辨率浅层地震探测和成像技术。电火花震源可以提供高宽频地震子波，进而提高数据纵向分辨率，而小道距地震电缆可以提供更高的数据横向分辨率，因此，海洋高分辨率多道地震探测具备常规地震勘探无法比拟的纵、横向分辨率优势，这也使其被广泛应用于海上工程建设（如跨海通道、滨海核电厂址、海上风电等）和海底新型资源勘查（如天然气水合物、深海稀土等）等方面，实现对海底浅层目标体更为精细、可靠的成像。受作业环境和设备限制，海洋高分辨率多道地震探测通常没有使用地震电缆深度控制水鸟，因此，无法像常规地震勘探那样在外业采集地震资料时，控制并获取地震电缆各检波点的沉放深度，用于后续数据处理。海洋小道距（道间距≤6.25m）地震电缆通常设计为海水中零浮力拖缆，但是，由于受到海水非均质性以及海洋环境涌浪作用的影响，在拖曳作业过程中地震电缆在海水中的沉放深度是在不断变化的，其变化剧烈程度跟涌浪强度成正相关。因此，对于海洋高分辨率地震探测而言，涌浪作用对反射地震成像质量的影响较为明显。在现场数据处理过程中，忽略涌浪校正或涌浪校正精度不高会导致成像剖面模糊、地层成像连续性变差，甚至影响对外业采集地震资料的现场质量控制与评价。中国专利CN111983686A公开了一种基于SEGY的浅地层剖面原始数据的可视化显示方法，采用涌浪滤波器进行涌浪静校正。然而通过涌浪滤波器进行涌浪静校正，无法适应各种采集环境获得精准的静校正量，从而影响成像剖面质量。
发明内容
针对现有技术存在的问题，本发明提供海洋高分辨率地震数据涌浪相关静校正量快速估算方法，该方法能够弥补海洋高分辨率多道地震探测难以装配地震电缆深度控制水鸟的不足，获得的涌浪相关静校正量可应用于地震数据现场处理过程，以补偿涌浪作用导致地震电缆沉放深度不断变化对地震数据的影响，实现对外业数据及时有效的质量控制。
海洋高分辨率地震数据涌浪相关静校正量快速估算方法，所述方法包括：
步骤S1：抽取共偏移距道集记录；
步骤S2：在各共偏移距道集记录上拾取海底反射波同相轴，获得各共偏移距道集记录的海底反射波到时实际拾取值，对所述各共偏移距道集记录的海底反射波到时实际拾取值进行格式转换，获得各叠前炮集记录的海底反射波到时实际拾取值，得到各叠前炮集记录的海底反射波同相轴实际时距曲线；
步骤S3：以在最小共偏移距道集记录上拾取的海底反射波同相轴为基础，结合观测系统参数获取海底水深地形数据；所述最小共偏移距道集记录是指将叠前炮集记录中所有偏移距等于最小偏移距的地震道数据抽取出来，形成的共偏移距道集记录；
步骤S4：根据所述海底水深地形数据和观测系统参数，计算各叠前炮集记录的海底反射波同相轴理论时距曲线；
步骤S5：对所述各叠前炮集记录的海底反射波同相轴理论时距曲线与各叠前炮集记录的海底反射波同相轴实际时距曲线求差，获取叠前炮集记录各地震道的涌浪相关静校正量。
在一种可能的实现方式中，所述步骤S1包括：
将M个叠前炮集记录抽取为N个共偏移距道集记录，并且按偏移距由小到大将各共偏移距道集记录依次定义为COG1、COG2、……、COGN，N为小道距地震电缆检波点总个数，各共偏移距道集记录的地震道道数等于M。
在一种可能的实现方式中，所述步骤S2包括：
分别对各共偏移距道集记录的海底反射波同相轴进行追踪拾取，获得各共偏移距道集记录的海底反射波到时实际拾取值，得到各共偏移距道集记录的海底反射波同相轴实际时距曲线；所述各共偏移距道集记录的海底反射波到时实际拾取值/>记为：其中，/>，/>表示第N个共偏移距道集记录的第M个地震道的海底反射波到时实际拾取值；
将各共偏移距道集记录的海底反射波到时实际拾取值，按照叠前炮集记录对应的数据格式抽取为各叠前炮集记录的海底反射波到时实际拾取值/>，得到各叠前炮集记录的海底反射波同相轴实际时距曲线，所述各叠前炮集记录的海底反射波到时实际拾取值/>记为：/>其中，，/>表示第M个叠前炮集记录的第N个地震道的海底反射波到时实际拾取值。
在一种可能的实现方式中，所述步骤S3包括：
根据最小共偏移距道集记录的海底反射波到时实际拾取值、海水声速、以及最小偏移距参数，获取海底水深地形数据，所述海底水深地形数据的计算公式如下：其中，/>表示第m个叠前炮集记录对应的电火花震源与首个地震电缆检波点中间位置的海底水深值，/>表示最小共偏移距道集记录的第m个地震道的海底反射波到时实际拾取值，/>表示海水声速，/>表示最小偏移距的设定值，/>，M为叠前炮集记录总个数。
在一种可能的实现方式中，在所述根据最小共偏移距道集记录的海底反射波到时实际拾取值、海水声速、以及最小偏移距参数，获取海底水深地形数据之前，还包括：对最小共偏移距道集记录的海底反射波到时实际拾取值/>进行平滑处理，剔除其中的高频成分。
在一种可能的实现方式中，所述步骤S4包括：
根据所述海底水深地形数据以及海水声速，构建地震作业测线对应的海水层速度-深度模型；
根据观测系统参数，将理想的炮检位置关系置于海水层速度-深度模型中，利用地震数值模拟技术计算各叠前炮集记录的海底反射波到时理论计算值，得到各叠前炮集记录的海底反射波同相轴理论时距曲线；所述各叠前炮集记录的海底反射波到时理论计算值/>为：/>其中，，/>表示在无涌浪作用的理想拖曳状态下，第M个叠前炮集记录的第N个地震道的海底反射波到时理论计算值。
在一种可能的实现方式中，所述叠前炮集记录各地震道的涌浪相关静校正量的计算公式为：其中，/>表示第m个叠前炮集记录的第n个地震道对应的涌浪相关静校正量，/>表示第m个叠前炮集记录的第n个地震道的海底反射波到时理论计算值，/>表示第m个叠前炮集记录的第n个地震道的海底反射波到时实际拾取值，/>，M为叠前炮集记录总个数，/>，N为小道距地震电缆检波点总个数。
基于以上发明内容，相对于现有技术，本发明其基于海底反射地震记录和观测系统参数对涌浪相关的静校正量进行快速估算，能够有效改善现场处理地震成像剖面的分辨率和信噪比，为外业地震资料现场质量评定和施工方案优化提供高效且可靠的数据支持。
附图说明
为了更清楚说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。
图1为本发明实施例提供的海洋高分辨率地震数据涌浪相关静校正量快速估算方法的流程示意图；
图2为本发明实施例提供的最小共偏移距道集记录剖面浅部显示图；
图3为本发明实施例提供的最小共偏移距道集记录的海底反射波同相轴实际时距曲线；
图4为本发明实施例提供的各叠前炮集记录的海底反射波同相轴实际时距曲线；
图5为本发明实施例提供的海底水深地形数据图；
图6为本发明实施例提供的第50个叠前炮集记录的海底反射波同相轴理论时距曲线；
图7为本发明实施例提供的各叠前炮集记录的海底反射波同相轴理论时距曲线；
图8为本发明实施例提供的叠前炮集记录各地震道的涌浪相关静校正量的示意图；
图9A-9B为本发明实施例提供的采用本发明方法获取的涌浪相关静校正量校正前后的地震成像剖面对比图。
具体实施方式
为了更好的理解本发明的技术方案，下面结合附图对本发明实施例进行详细描述。
应当明确，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
参见图1，为本发明实施例提供的海洋高分辨率地震数据涌浪相关静校正量快速估算方法的流程示意图。如图1所示，所述方法具体包括：
步骤S1：抽取共偏移距道集记录。
所述抽取共偏移距道集记录是指按照偏移距大小，将各叠前炮集记录中偏移距相同的地震道数据抽取出来，形成共偏移距道集记录。
具体的，所述步骤S1包括：将M个叠前炮集记录抽取为N个共偏移距道集记录，并且按偏移距由小到大将各共偏移距道集记录依次定义为COG1、COG2、……、COGN，N为小道距地震电缆检波点总个数，各共偏移距道集记录的地震道道数等于M。
步骤S2：在各共偏移距道集记录上拾取海底反射波同相轴，获得各共偏移距道集记录的海底反射波到时实际拾取值，对所述各共偏移距道集记录的海底反射波到时实际拾取值进行格式转换，获得各叠前炮集记录的海底反射波到时实际拾取值，得到各叠前炮集记录的海底反射波同相轴实际时距曲线。
具体的，所述步骤S2包括：分别对各共偏移距道集记录的海底反射波同相轴进行追踪拾取，获得各共偏移距道集记录的海底反射波到时实际拾取值，得到各共偏移距道集记录的海底反射波同相轴实际时距曲线；所述各共偏移距道集记录的海底反射波到时实际拾取值/>记为：/>其中，，/>表示第N个共偏移距道集记录的第M个地震道的海底反射波到时实际拾取值；
将各共偏移距道集记录的海底反射波到时实际拾取值，按照叠前炮集记录对应的数据格式抽取为各叠前炮集记录的海底反射波到时实际拾取值/>，得到各叠前炮集记录的海底反射波同相轴实际时距曲线，所述各叠前炮集记录的海底反射波到时实际拾取值/>记为：/>其中，，/>表示第M个叠前炮集记录的第N个地震道的海底反射波到时实际拾取值。
步骤S3：以在最小共偏移距道集记录上拾取的海底反射波同相轴为基础，结合观测系统参数获取海底水深地形数据；所述最小共偏移距道集记录是指将叠前炮集记录中所有偏移距等于最小偏移距的地震道数据抽取出来，形成的共偏移距道集记录。
考虑到海洋高分辨率多道地震探测的最小偏移距参数往往较小，通常设定为0，或者1到2倍道间距，首个地震电缆检波点与电火花震源的相对位置浮动变化较小，因此，海底水深地形数据的获取是以最小共偏移距道集记录COG1的海底反射波到时实际拾取值为基础计算的。
具体的，所述步骤S3包括：根据最小共偏移距道集记录的海底反射波到时实际拾取值、海水声速、以及最小偏移距参数，获取海底水深地形数据，所述海底水深地形数据的计算公式如下：/>其中，/>表示第m个叠前炮集记录对应的电火花震源与首个地震电缆检波点中间位置的海底水深值，/>表示最小共偏移距道集记录的第m个地震道的海底反射波到时实际拾取值，/>表示海水声速，表示最小偏移距的设定值，/>，M为叠前炮集记录总个数。
在一个实施例中，海水声速的取值为海水声速估计值或测量值。
在一个实施例中，为确保最终海底水深地形数据的平滑性，在根据最小共偏移距道集记录的海底反射波到时实际拾取值、海水声速、以及最小偏移距参数，获取海底水深地形数据之前，还要对最小共偏移距道集记录的海底反射波到时实际拾取值/>进行平滑处理，剔除其中的高频成分。
步骤S4：根据所述海底水深地形数据和观测系统参数，计算各叠前炮集记录的海底反射波同相轴理论时距曲线。
具体的，所述步骤S4包括：根据所述海底水深地形数据以及海水声速，构建地震作业测线对应的海水层速度-深度模型。
根据观测系统参数，将理想的炮检位置关系置于海水层速度-深度模型中，利用地震数值模拟技术计算各叠前炮集记录的海底反射波到时理论计算值，得到各叠前炮集记录的海底反射波同相轴理论时距曲线；所述各叠前炮集记录的海底反射波到时理论计算值/>为：/>其中，，/>表示在无涌浪作用的理想拖曳状态下，第M个叠前炮集记录的第N个地震道的海底反射波到时理论计算值。
在一个实施例中，所述观测系统参数包括观测系统的炮点位置、最小偏移距和道间距；所述理想的炮检位置关系是指小道距地震电缆整体水平，且沉放深度恒定。
在一个实施例中，所述地震数值模拟技术包括如射线追踪算法。
步骤S5：对所述各叠前炮集记录的海底反射波同相轴理论时距曲线与各叠前炮集记录的海底反射波同相轴实际时距曲线求差，获取叠前炮集记录各地震道的涌浪相关静校正量。
所述叠前炮集记录各地震道的涌浪相关静校正量的计算公式为：其中，/>表示第m个叠前炮集记录的第n个地震道对应的涌浪相关静校正量，/>表示第m个叠前炮集记录的第n个地震道的海底反射波到时理论计算值，/>表示第m个叠前炮集记录的第n个地震道的海底反射波到时实际拾取值，，M为叠前炮集记录总个数，/>，N为小道距地震电缆检波点总个数。
下面给出一个具体实施例，说明本发明实施例中海洋高分辨率地震数据涌浪相关静校正量快速估算方法的具体应用。在本具体实施例中，共选取了3413个叠前炮集记录，即M=3413，小道距地震电缆的检波点总个数为24，即N=24。电火花震源和首个地震电缆检波点的距离为3.125m，即最小偏移距为3.125m，按照偏移距由小到大将叠前炮集记录抽取为共偏移距道集记录COG1、COG2、COG3、……、COG24。将叠前炮集记录中所有偏移距等于最小偏移距的地震道数据抽取出来，即获得最小共偏移距道集记录COG1。图2为最小共偏移距道集记录COG1剖面浅部显示图，从图2中清晰可见海底反射波同相轴。
分别对24个共偏移距道集记录上的海底反射波同相轴进行追踪拾取，获得各共偏移距道集记录的海底反射波到时实际拾取值。图3为最小共偏移距道集记录COG1的海底反射波同相轴实际时距曲线，显示了最小共偏移距道集记录COG1的海底反射波到时实际拾取值/>。将各共偏移距道集记录的海底反射波到时实际拾取值/>，按照叠前炮集记录对应的数据格式抽取为各叠前炮集记录的海底反射波到时实际拾取值/>。图4为各叠前炮集记录的海底反射波同相轴实际时距曲线。
在本具体实施例中，电火花震源和首个地震电缆检波点的距离为3.125m，即最小偏移距为3.125m，叠前炮集记录总个数M为3413。将上述数据代入所述海底水深地形数据的计算公式，可以求得海底水深地形数据。图5为获取的海底水深地形数据图。
基于图5的海底水深地形数据图和海水声速的估计值，构建地震作业测线对应的海水层速度-深度模型。本具体实施例所选观测系统设计参数具体为：炮间距为6.25m；道间距为3.125m；最小偏移距为3.125m。按照观测系统设计参数将模拟的震源和检波点置于海水层速度-深度模型的海面位置，保持沉放深度恒定。根据地震数值模拟技术，获得小道距地震电缆在理想拖曳情况下，各叠前炮集记录的海底反射波到时理论计算值/>，。
如图6所示，显示了第50个叠前炮集记录的海底反射波同相轴理论时距曲线，图6中正方形实心黑点为第50个叠前炮集记录的海底反射波到时理论计算值，。图7为小道距地震电缆在理想拖曳状态下各叠前炮集记录的海底反射波同相轴理论时距曲线。
将理想情况下的各叠前炮集记录的海底反射波到时理论计算值与各叠前炮集记录的海底反射波到时实际拾取值/>求差，获得叠前炮集记录各地震道的涌浪相关静校正量。图8为叠前炮集记录各地震道的涌浪相关静校正量的示意图。
图9A-图9B为采用本发明方法获取的涌浪相关静校正量校正前后的地震成像剖面对比图，其中，图9A为未采用本发明方法获取的涌浪相关静校正量校正，进行处理得到的地震成像结果，图9B为采用本发明方法获取的涌浪相关静校正量校正后，进行处理得到的地震成像结果。可以看出，本发明方法能够有效改善海底地层成像效果，信噪比、分辨率有所提升，同相轴连续性显著增强。
综上所述，本发明针对海洋高分辨率多道地震探测技术应用过程中，由于缺少对小道距地震电缆沉放深度的调节控制，海洋环境涌浪作用会导致地震成像模糊、地层连续性变差的问题，基于共偏移距道集记录的海底反射波信息和观测系统参数，通过构建海水层速度-深度模型，利用地震数值模拟技术获取小道距地震电缆在理想拖曳状态下接收到的海底反射波同相轴理论时距曲线，结合叠前炮集记录的海底反射波同相轴实际时距曲线，从而实现对涌浪相关静校正量的快速估算。本方法不需要对作业过程中各检波点的实际位置进行求取，计算效率较高，并且能够满足现场地震数据处理对静校正量的精度要求，可以有效提高地震成像剖面的信噪比和分辨率，改善海底地层成像质量，使其能够客观可靠地反映采集数据品质和海底目标体情况，为外业地震资料现场质量评定和施工方案优化提供高效且可靠的数据支持。
以上所述，仅为本发明的具体实施方式，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。本发明的保护范围应以所述权利要求的保护范围为准。
Claims (5)
1.一种海洋高分辨率地震数据涌浪相关静校正量快速估算方法，其特征在于，所述方法包括：
步骤S1：抽取共偏移距道集记录；
步骤S2：在各共偏移距道集记录上拾取海底反射波同相轴，获得各共偏移距道集记录的海底反射波到时实际拾取值，对所述各共偏移距道集记录的海底反射波到时实际拾取值进行格式转换，获得各叠前炮集记录的海底反射波到时实际拾取值，得到各叠前炮集记录的海底反射波同相轴实际时距曲线；
步骤S3：以在最小共偏移距道集记录上拾取的海底反射波同相轴为基础，结合观测系统参数获取海底水深地形数据；所述最小共偏移距道集记录是指将叠前炮集记录中所有偏移距等于最小偏移距的地震道数据抽取出来，形成的共偏移距道集记录；
步骤S4：根据所述海底水深地形数据和观测系统参数，计算各叠前炮集记录的海底反射波同相轴理论时距曲线；
步骤S5：对所述各叠前炮集记录的海底反射波同相轴理论时距曲线与各叠前炮集记录的海底反射波同相轴实际时距曲线求差，获取叠前炮集记录各地震道的涌浪相关静校正量；
所述步骤S3包括：
根据最小共偏移距道集记录的海底反射波到时实际拾取值、海水声速、以及最小偏移距参数，获取海底水深地形数据，所述海底水深地形数据的计算公式如下：其中，/>表示第m个叠前炮集记录对应的电火花震源与首个地震电缆检波点中间位置的海底水深值，/>表示最小共偏移距道集记录的第m个地震道的海底反射波到时实际拾取值，/>表示海水声速，/>表示最小偏移距的设定值，/>，M为叠前炮集记录总个数；
所述步骤S4包括：
根据所述海底水深地形数据以及海水声速，构建地震作业测线对应的海水层速度-深度模型；
根据观测系统参数，将理想的炮检位置关系置于海水层速度-深度模型中，利用地震数值模拟技术计算各叠前炮集记录的海底反射波到时理论计算值，得到各叠前炮集记录的海底反射波同相轴理论时距曲线；所述各叠前炮集记录的海底反射波到时理论计算值/>为：/>其中，，/>表示在无涌浪作用的理想拖曳状态下，第M个叠前炮集记录的第N个地震道的海底反射波到时理论计算值，N为小道距地震电缆检波点总个数。
2.根据权利要求1所述的一种海洋高分辨率地震数据涌浪相关静校正量快速估算方法，其特征在于，所述步骤S1包括：
将M个叠前炮集记录抽取为N个共偏移距道集记录，并且按偏移距由小到大将各共偏移距道集记录依次定义为COG1、COG2、……、COGN，N为小道距地震电缆检波点总个数，各共偏移距道集记录的地震道道数等于M。
3.根据权利要求1所述的一种海洋高分辨率地震数据涌浪相关静校正量快速估算方法，其特征在于，所述步骤S2包括：
分别对各共偏移距道集记录的海底反射波同相轴进行追踪拾取，获得各共偏移距道集记录的海底反射波到时实际拾取值，得到各共偏移距道集记录的海底反射波同相轴实际时距曲线；所述各共偏移距道集记录的海底反射波到时实际拾取值/>记为：其中，/>，/>表示第N个共偏移距道集记录的第M个地震道的海底反射波到时实际拾取值；
将各共偏移距道集记录的海底反射波到时实际拾取值，按照叠前炮集记录对应的数据格式抽取为各叠前炮集记录的海底反射波到时实际拾取值/>，得到各叠前炮集记录的海底反射波同相轴实际时距曲线，所述各叠前炮集记录的海底反射波到时实际拾取值/>记为：/>其中，/>，/>表示第M个叠前炮集记录的第N个地震道的海底反射波到时实际拾取值。
4.根据权利要求3所述的一种海洋高分辨率地震数据涌浪相关静校正量快速估算方法，其特征在于，在所述根据最小共偏移距道集记录的海底反射波到时实际拾取值、海水声速、以及最小偏移距参数，获取海底水深地形数据之前，还包括：对最小共偏移距道集记录的海底反射波到时实际拾取值/>进行平滑处理，剔除其中的高频成分。
5.根据权利要求1所述的一种海洋高分辨率地震数据涌浪相关静校正量快速估算方法，其特征在于，所述叠前炮集记录各地震道的涌浪相关静校正量的计算公式为：其中，/>表示第m个叠前炮集记录的第n个地震道对应的涌浪相关静校正量，/>表示第m个叠前炮集记录的第n个地震道的海底反射波到时理论计算值，/>表示第m个叠前炮集记录的第n个地震道的海底反射波到时实际拾取值，，M为叠前炮集记录总个数，/>，N为小道距地震电缆检波点总个数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311188443.3A CN116931088B (zh) | 2023-09-15 | 2023-09-15 | 海洋高分辨率地震数据涌浪相关静校正量快速估算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311188443.3A CN116931088B (zh) | 2023-09-15 | 2023-09-15 | 海洋高分辨率地震数据涌浪相关静校正量快速估算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116931088A CN116931088A (zh) | 2023-10-24 |
| CN116931088B true CN116931088B (zh) | 2023-12-19 |
Family
ID=88377427
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311188443.3A Active CN116931088B (zh) | 2023-09-15 | 2023-09-15 | 海洋高分辨率地震数据涌浪相关静校正量快速估算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116931088B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109188527A (zh) * | 2018-10-31 | 2019-01-11 | 中国石油化工股份有限公司 | 滩浅海地区快速建立三维近海底速度模型的方法 |
| CN109782355A (zh) * | 2018-12-18 | 2019-05-21 | 中国石油天然气集团有限公司 | Obs检波点漂移的检测方法及装置 |
| KR102003466B1 (ko) * | 2019-02-08 | 2019-07-24 | 한국지질자원연구원 | 다중채널 3차원 지층탐사자료의 너울 영향 보정방법, 이를 이용한 해양 3d 탄성파 탐사방법 |
| CN111487678A (zh) * | 2020-04-20 | 2020-08-04 | 青岛海洋地质研究所 | 一种确定高分辨率小多道地震最小偏移距和系统延迟的分析方法 |
| CN112883564A (zh) * | 2021-02-01 | 2021-06-01 | 中国海洋大学 | 一种基于随机森林的水体温度预测方法及预测系统 |
| CN114460649A (zh) * | 2022-04-14 | 2022-05-10 | 自然资源部第一海洋研究所 | 一种深海近底拖曳式多道地震接收阵列形态重建方法 |
| CN114488280A (zh) * | 2022-01-17 | 2022-05-13 | 自然资源部第二海洋研究所 | 一种校正高分辨率地震不平电缆的方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9103943B2 (en) * | 2011-11-28 | 2015-08-11 | Fugro-Geoteam As | Acquisition and processing of multi-source broadband marine seismic data |
-
2023
- 2023-09-15 CN CN202311188443.3A patent/CN116931088B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109188527A (zh) * | 2018-10-31 | 2019-01-11 | 中国石油化工股份有限公司 | 滩浅海地区快速建立三维近海底速度模型的方法 |
| CN109782355A (zh) * | 2018-12-18 | 2019-05-21 | 中国石油天然气集团有限公司 | Obs检波点漂移的检测方法及装置 |
| KR102003466B1 (ko) * | 2019-02-08 | 2019-07-24 | 한국지질자원연구원 | 다중채널 3차원 지층탐사자료의 너울 영향 보정방법, 이를 이용한 해양 3d 탄성파 탐사방법 |
| CN111487678A (zh) * | 2020-04-20 | 2020-08-04 | 青岛海洋地质研究所 | 一种确定高分辨率小多道地震最小偏移距和系统延迟的分析方法 |
| CN112883564A (zh) * | 2021-02-01 | 2021-06-01 | 中国海洋大学 | 一种基于随机森林的水体温度预测方法及预测系统 |
| CN114488280A (zh) * | 2022-01-17 | 2022-05-13 | 自然资源部第二海洋研究所 | 一种校正高分辨率地震不平电缆的方法 |
| CN114460649A (zh) * | 2022-04-14 | 2022-05-10 | 自然资源部第一海洋研究所 | 一种深海近底拖曳式多道地震接收阵列形态重建方法 |
Non-Patent Citations (1)
| Title |
|---|
| 海洋短排列多道反射地震数据观测系统重定义与沉放深度校正;张圣 等;海洋学报;第42卷(第3期);第145页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116931088A (zh) | 2023-10-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3078991B1 (en) | Method for swell effect and mis-tie correction in high-resolution seismic data using multi-beam echo sounder data | |
| US8705317B2 (en) | Method for imaging of targeted reflectors | |
| CN105629220B (zh) | 一种基于单水听器的深海水声被动测距方法 | |
| Cobra et al. | Geometric distortions in side-scan sonar images: a procedure for their estimation and correction | |
| CN113625337B (zh) | 一种极浅水高精度地震资料快速成像方法 | |
| CN114460649B (zh) | 一种深海近底拖曳式多道地震接收阵列形态重建方法 | |
| CN104268848A (zh) | 一种海洋内波波速监测的方法 | |
| Duarte et al. | Advanced processing for UHR3D shallow marine seismic surveys | |
| Köhler et al. | Potentials and pitfalls of permafrost active layer monitoring using the HVSR method: A case study in Svalbard | |
| CN109212598B (zh) | 基于直达波反演的三维空间二次定位方法 | |
| CN116931088B (zh) | 海洋高分辨率地震数据涌浪相关静校正量快速估算方法 | |
| US20130046472A1 (en) | Method of determining the relative position of two detectors at the bottom of the sea | |
| CN113075732B (zh) | 一种消除高分辨率小多道地震地层异常起伏的方法 | |
| CN114838808A (zh) | 一种基于浅海声信道响应的打桩冲击噪声源级测量系统及方法 | |
| Kim et al. | A simple method of correction for profile-length water-column height variations in high-resolution, shallow-water seismic data | |
| CN117908108B (zh) | 一种海洋地震实时监测系统 | |
| CN115980856B (zh) | 适用复杂海底界面的深拖式多道地震炮检精确定位方法 | |
| CN114722590B (zh) | 基于地球物理模型的随机采集观测系统设计优选方法 | |
| CN113568042B (zh) | 一种提高时移地震数据可重复性的拖缆采集方法和系统 | |
| CN112925020B (zh) | 一种参量阵型浅剖智能能量补偿方法 | |
| CN116466390B (zh) | 一种大型水库诱发地震实时监测定位方法 | |
| CN116660979B (zh) | 一种基于Kaiser时窗积分的OBN资料全波形反演方法 | |
| CN111856581B (zh) | 一种obs时钟漂移校正方法及处理终端 | |
| Cazenave et al. | Quantitative bedform analysis using decimetre resolution swath bathymetry | |
| Barnes et al. | Sea surface source-side static optimization for 4D seismic |
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 |