CN112379420A - 高精度弯曲测线叠前时间域成像方法及装置 - Google Patents
高精度弯曲测线叠前时间域成像方法及装置 Download PDFInfo
- Publication number
- CN112379420A CN112379420A CN202011189523.7A CN202011189523A CN112379420A CN 112379420 A CN112379420 A CN 112379420A CN 202011189523 A CN202011189523 A CN 202011189523A CN 112379420 A CN112379420 A CN 112379420A
- Authority
- CN
- China
- Prior art keywords
- seismic
- coordinate
- offset
- dimensional array
- point
- 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.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 79
- 238000005452 bending Methods 0.000 title claims abstract description 47
- 238000005259 measurement Methods 0.000 title abstract description 13
- 238000000034 method Methods 0.000 claims abstract description 36
- 230000005012 migration Effects 0.000 claims abstract description 14
- 238000013508 migration Methods 0.000 claims abstract description 14
- 238000001514 detection method Methods 0.000 claims description 20
- 230000015654 memory Effects 0.000 claims description 19
- 239000000523 sample Substances 0.000 claims description 15
- 238000003860 storage Methods 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000005096 rolling process Methods 0.000 claims description 6
- 239000013317 conjugated microporous polymer Substances 0.000 description 19
- 238000012545 processing Methods 0.000 description 7
- IERHLVCPSMICTF-XVFCMESISA-N CMP group Chemical group P(=O)(O)(O)OC[C@@H]1[C@H]([C@H]([C@@H](O1)N1C(=O)N=C(N)C=C1)O)O IERHLVCPSMICTF-XVFCMESISA-N 0.000 description 5
- 210000003643 myeloid progenitor cell Anatomy 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000003491 array Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 210000003746 feather Anatomy 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 239000000126 substance Substances 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/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- 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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
- G01V2210/74—Visualisation of seismic data
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)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种高精度弯曲测线叠前时间域成像方法及装置，该方法包括：获取弯曲测线地震数据；根据预设炮检距参数组将炮检距分为n组，其中，所述炮检距参数包括：最小炮检距、最大炮检距和炮检距增量；根据所述最小CMP号和所述最大CMP号计算CMP总数；根据炮检距分组个数n以及所述CMP总数建立三维数组；根据所述弯曲测线地震数据中每一个地震道的地震数据生成每一个地震道各自对应的三维数组；根据每一个地震道各自对应的三维数组计算每一个地震道各自对应的旅行时，并根据所述旅行时进行弯曲测线叠前时间偏移。本发明实现了针对弯曲测线高精度的叠前时间域成像。
Description
技术领域
本发明涉及反射波地震数据处理过程中的叠前时间成像技术，具体而言，涉及一种高精度弯曲测线叠前时间域成像方法及装置。
背景技术
地震采集方式通常可以划分为两类：一类是二维工区，地震观测是沿着地面上一条直线进行的，所有激发点和接收点全部被布设在该直线上，炮点激发的地震波向下传播，经地下岩性界面反射后传回地面被检波点接收；另一类是三维工区，地震观测是在地面上一个区域内进行的，所有激发点和接收点等间距地被布设在该区域内，炮点激发的地震波向下传播，经地下岩性界面反射后传回地面被三维平面内所有检波器接收。和二维观测相比，三维勘探能够对侧面波正确成像、能够观测地震属性随方位角变化的规律，所以目前地震勘探大部分工区都是采用三维观测。但是受到自然条件的限制，并不是所有地震勘探都能按照二维或者三位方式完成。比如在复杂山地区，受到地表高程剧烈变化的影响，地震勘探只能在山间洼地开展，地震测线随着地形走向而弯曲；再比如在海面上进行二维观测时，受到洋流的影响，拖缆会沿着洋流方向有一定横向漂移。以上两种观测方式既不是二维也不是三维，本行业把它们归类为弯曲测线地震勘探。
针对于标准的二维和三维地震勘探，地震数据处理分别发展了各自完整地震资料处理流程。克希霍夫叠前时间偏移技术是目前生产中应用最为广泛的时间域成像方法。但是对弯曲地震勘探来说，以前还没有人针对于该类数据专门进行时间域成像方法研究，只是把这类数据近似成二维数据处理，当测线弯曲程度较小时，这样的近似可以被生产上接受，但是测线弯曲程度较大时，还没有正确的成像方法。
发明内容
本发明为了解决上述背景技术中的至少一个技术问题，提出了一种高精度弯曲测线叠前时间域成像方法及装置。
为了实现上述目的，根据本发明的一个方面，提供了一种高精度弯曲测线叠前时间域成像方法，该方法包括：
获取弯曲测线地震数据，所述弯曲测线地震数据按照共中心点道集顺序存放，所述弯曲测线地震数据的卷头中包含：最小CMP号以及最大CMP号，所述弯曲测线地震数据中每一个地震道的道头字包含：炮点X坐标、炮点Y坐标、检波点X坐标、检波点Y坐标、中心点X坐标、中心点Y坐标以及CMP号；
根据预设炮检距参数组将炮检距分为n组，其中，所述炮检距参数包括：最小炮检距、最大炮检距和炮检距增量；
根据所述最小CMP号和所述最大CMP号计算CMP总数；
根据炮检距分组个数n以及所述CMP总数建立三维数组，其中：所述三维数组的维度一的长度等于所述炮检距分组个数n，用于存放炮检距分组信息；所述三维数组的维度二的长度等于6，用于存放CMP号、地震道中点X坐标、地震道中点Y坐标、中心点X坐标、中心点Y坐标以及地震道数；所述三维数组的维度三的长度等于所述CMP总数，用于存放CMP点信息；
根据所述弯曲测线地震数据中每一个地震道的地震数据生成每一个地震道各自对应的三维数组；
根据每一个地震道各自对应的三维数组计算每一个地震道各自对应的旅行时，并根据所述旅行时进行弯曲测线叠前时间偏移。
可选的，所述根据所述弯曲测线地震数据中每一个地震道的地震数据生成每一个地震道各自对应的三维数组，包括：
根据地震道的炮点X坐标、炮点Y坐标、检波点X坐标以及检波点Y坐标计算出该地震道的炮检距；
根据该地震道的炮检距确定该地震道的炮检距分组信息，并将该地震道的炮检距分组信息存放在该地震道对应的三维数组的维度一中。
可选的，所述根据每一个地震道各自对应的三维数组计算每一个地震道各自对应的旅行时，包括：
根据三维数组中的地震道中点X坐标以及三维数组中的地震道中点Y坐标分别计算出炮点到成像点的距离以及检波点到成像点的距离；
根据炮点到成像点的距离和检波点到成像点的距离计算出旅行时。
可选的，所述根据每一个地震道各自对应的三维数组计算每一个地震道各自对应的旅行时，包括：
各个地震道在对同一个成像点计算旅行时时，选区各自炮检距的成像点坐标进行计算。
为了实现上述目的，根据本发明的另一方面，提供了一种高精度弯曲测线叠前时间域成像装置，该装置包括：
弯曲测线地震数据获取单元，用于获取弯曲测线地震数据，所述弯曲测线地震数据按照共中心点道集顺序存放，所述弯曲测线地震数据的卷头中包含：最小CMP号以及最大CMP号，所述弯曲测线地震数据中每一个地震道的道头字包含：炮点X坐标、炮点Y坐标、检波点X坐标、检波点Y坐标、中心点X坐标、中心点Y坐标以及CMP号；
炮检距分组单元，用于根据预设炮检距参数组将炮检距分为n组，其中，所述炮检距参数包括：最小炮检距、最大炮检距和炮检距增量；
CMP总数计算单元，用于根据所述最小CMP号和所述最大CMP号计算CMP总数；
三维数组设置单元，用于根据炮检距分组个数n以及所述CMP总数建立三维数组，其中：所述三维数组的维度一的长度等于所述炮检距分组个数n，用于存放炮检距分组信息；所述三维数组的维度二的长度等于6，用于存放CMP号、地震道中点X坐标、地震道中点Y坐标、中心点X坐标、中心点Y坐标以及地震道数；所述三维数组的维度三的长度等于所述CMP总数，用于存放CMP点信息；
三维数组数据生成单元，用于根据所述弯曲测线地震数据中每一个地震道的地震数据生成每一个地震道各自对应的三维数组；
旅行时计算单元，用于根据每一个地震道各自对应的三维数组计算每一个地震道各自对应的旅行时，并根据所述旅行时进行弯曲测线叠前时间偏移。
为了实现上述目的，根据本发明的另一方面，还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述高精度弯曲测线叠前时间域成像方法中的步骤。
为了实现上述目的，根据本发明的另一方面，还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序在计算机处理器中执行时实现上述高精度弯曲测线叠前时间域成像方法中的步骤。
本发明的有益效果为：本发明根据弯曲测线地震数据中每一个地震道的地震数据生成每一个地震道各自对应的三维数组，进而根据每一个地震道各自对应的三维数组计算每一个地震道各自对应的旅行时，最后根据所述旅行时进行弯曲测线叠前时间偏移，实现了针对弯曲测线高精度的叠前时间域成像。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1是本发明实施例高精度弯曲测线叠前时间域成像方法的流程图；
图2是本发明实施例确定炮检距分组信息的流程图；
图3是本发明实施例计算旅行时的流程图；
图4是水平地表旅行时计算图；
图5是本条测线叠前时间偏移成像结果某一位置处的CIP道集示意图；
图6是本条测线叠前时间偏移成像结果的叠加剖面示意图；
图7是本发明实施例计算机设备示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
Kirchhoff积分法偏移用公式为：
公式(1)可以进行深度偏移，也可以完成时间偏移，当完成时间偏移式时，旅行时计算可以如图4所示，以下公式(2)是经典的双平方根旅行时计算公式。
公式(2)中变量可以参看图4，具体含义如下：
τ0：成像时间；
hs：炮点到成像点的距离os；
hr：检波点到成像点的距离or；
vrms：成像点位置处的均方根速度。
对于标准的二维和三维工区来说，每个成像点都有准确的空间位置，公式(2)表达的旅行时和地震波实际走时完全吻合。
但是对弯曲地震勘探来说，例如海洋环境，由于受到洋流的影响，拖缆(检波点)方向和炮线方向存在一定的夹角，导致炮、检点并不在一条直线上，相应地其中心点也是分布在平面上的一个区域内。弯曲测线中心点并没有分布在一条直线上，而是弥散在一个不规则的二维平面上，如果继续沿用公式(2)计算旅行时就会带来很大的误差。
图1是本发明实施例高精度弯曲测线叠前时间域成像方法的流程图，如图1所示，本发明实施例的高精度弯曲测线叠前时间域成像方法包括步骤S101至步骤S106。
步骤S101，获取弯曲测线地震数据，所述弯曲测线地震数据按照共中心点道集顺序存放，所述弯曲测线地震数据的卷头中包含：最小CMP号(minCMP)以及最大CMP号(maxCMP)，所述弯曲测线地震数据中每一个地震道的道头字包含：炮点X坐标(sx)、炮点Y坐标(sy)、检波点X坐标(rx)、检波点Y坐标(ry)、中心点X坐标(CMPx)、中心点Y坐标(CMPy)以及CMP号。
步骤S102，根据预设炮检距参数组将炮检距分为n组，其中，所述炮检距参数包括：最小炮检距(minOffset)、最大炮检距(maxOffset)和(offsetInc)。
在本发明一个实施例中，本步骤具体分组方法是：
炮检距分组个数(n)＝int((maxOffset-minOffset)/offsetInc)；
分组后第一组炮检距范围：minOffset～minOffset+1*offsetInc；
分组后第二组炮检距范围：minOffset～minOffset+2*offsetInc；
…
分组后第i组炮检距范围：minOffset～minOffset+i*offsetInc；
…
分组后第n组炮检距范围：minOffset～minOffset+n*offsetInc。
步骤S103，根据所述最小CMP号(minCMP)和所述最大CMP号(maxCMP)计算CMP总数。
在本发明实施例中，CMP总数可以由以下公式计算得出：
CMPs＝maxCMP–minCMP+1
其中，CMPs为CMP总数。
步骤S104，根据炮检距分组个数n以及所述CMP总数建立三维数组ARRAY(n,6,CMPs)，其中：所述三维数组ARRAY()的维度一的长度等于所述炮检距分组个数n，用于存放炮检距分组信息；所述三维数组的维度二的长度等于6，用于存放CMP号、地震道中点X坐标、地震道中点Y坐标、中心点X坐标、中心点Y坐标以及地震道数；所述三维数组的维度三的长度等于所述CMP总数，用于存放CMP点信息。
步骤S105，根据所述弯曲测线地震数据中每一个地震道的地震数据生成每一个地震道各自对应的三维数组。
在本发明实施例中，本步骤针对每个地震道先对三维数组ARRAY(n,6,CMPs)进行初赋0值操作，进而对每个地震道的地震数据进行以下操作：
读取本地震道的CMP号：CMPk；
读取本地震道的炮点坐标：sx和sy；
读取本地震道的检波点坐标：rx和ry；
读取本地震道CMP的中心点坐标：CMPx和CMPy；
计算本地震道的中点坐标：
计算本道炮检距：
根据本道炮检距offsetk和步骤S102中的炮检距区间确定本地震道的炮检距分组信息(分组编号k)。
进而本步骤对三维数组ARRAY(n,6,CMPs)进行以下操作，生成地震道对应的三维数组数据：
ARRAY(k,1,CMPk)＝CMPk；
ARRAY(k,2,CMPk)＝ARRAY(k,2,CMPk)+MIDx；
ARRAY(k,3,CMPk)＝ARRAY(k,3,CMPk)+MIDy；
ARRAY(k,4,CMPk)＝CMPx；
ARRAY(k,5,CMPk)＝CMPy；
ARRAY(k,6,CMPk)＝ARRAY(k,6,CMPk)+1；
每个地震道对应的三维数组数据包含以上6个数据，即ARRAY(k,2,CMPk)、ARRAY(k,2,CMPk)、ARRAY(k,3,CMPk)、ARRAY(k,4,CMPk)、ARRAY(k,5,CMPk)和ARRAY(k,6,CMPk)。
在本发明一个实施例中，本发明在对地震数据中所有地震道完成步骤S105的操作后，还对数组ARRAY()存放的地震道中点坐标计算平均值。具体做法用计算机语言表达为：按照第一列和第三例执行两次循环，完成以下操作：
在本发明一个实施例中，由于地震数据分布的不规则性，各个共中心点位置上每一个炮检距区间原始地震数据道个数是不相同的，甚至在某些炮检距上会出现地震道缺失现象，从而导致数组ARRAY(*,2,*)和ARRAY(*,3,*)的部分元素值为0。为解决这一问题，本发明还对数组ARRAY()沿着维度二进行插值，具体插值方法是：一次固定数组ARRAY()第一列元素的一个值(对应一个炮检距)，沿着第三维度(CMP)对ARRAY(*,2,*)和ARRAY(*,3,*)进行多项式拟合。
步骤S106，根据每一个地震道各自对应的三维数组计算每一个地震道各自对应的旅行时，并根据所述旅行时进行弯曲测线叠前时间偏移。
在本发明实施例中，本发明计算旅行时的公式为：
式中：
hs：炮点到成像点的距离；
hr：检波点到成像点的距离；
sx：炮点X坐标，从地震数据道头字中得到；
sy：炮点Y坐标，从地震数据道头字中得到；
rx：检波点X坐标，从地震数据道头字中得到；
ry：检波点Y坐标，从地震数据道头字中得到；
ARRAY(k,2,*)：本地震道所在炮检距某个成像点的X坐标，从本地震道对应的三维数组中得到；
ARRAY(k,3,*)：本地震道所在炮检距某个成像点的Y坐标，从本地震道对应的三维数组中得到。
在本发明实施例中，本发明弯曲测线叠前时间偏移其它步骤和普通地震工区相同。
图2是本发明实施例确定炮检距分组信息的流程图，如图2所示，上述步骤步骤S105的根据所述弯曲测线地震数据中每一个地震道的地震数据生成每一个地震道各自对应的三维数组具体包括步骤S201和步骤S202。
步骤S201，根据地震道的炮点X坐标、炮点Y坐标、检波点X坐标以及检波点Y坐标计算出该地震道的炮检距。
步骤S202，根据该地震道的炮检距确定该地震道的炮检距分组信息，并将该地震道的炮检距分组信息存放在该地震道对应的三维数组的维度一中。
图3是本发明实施例计算旅行时的流程图，如图3所示，上述步骤S106的根据每一个地震道各自对应的三维数组计算每一个地震道各自对应的旅行时，包括步骤S301和步骤S302。
步骤S301，根据三维数组中的地震道中点X坐标以及三维数组中的地震道中点Y坐标分别计算出炮点到成像点的距离以及检波点到成像点的距离。
步骤S302，根据炮点到成像点的距离和检波点到成像点的距离计算出旅行时。
在本发明一个实施例中，上述步骤S106在各个地震道在对同一个成像点计算旅行时时，选区各自炮检距的成像点坐标进行计算。
在本发明一个实施例中，本发明的，特点是弯曲测线每个炮检距要单独拟合成像线的坐标。
本发明实例说明如下：
2018年本发明技术已经在我国大型地震数据处理解释一体化系统(GeoEast)成功开发为工业化程序，并在国外某二维工区进行了试验性应用。例如国内某弯曲测线炮、检点位置及羽角分布范围在12.60-34.70之间，如此大角度的羽角存在导致该条拖缆采集的地震数据严重偏离了炮线方向，形成实际上的弯曲地震测线。
在本发明技术使用之前，这种地震数据一般都近似成二维地震测线并完成后续地震成像处理。图5是本条测线叠前时间偏移成像结果某一位置处的CIP道集，左图对应直线方法结果，右图对应本专利技术结果。由于近似的成像点位置和地震道真正的反射点位置存在较大的偏离，直线方法成像结果道集中只有近炮检距位置是正确的，中远炮检距比实际成像时间小，地震同相轴并没有完全拉平，存在一定的上翘现象。相反本专利方法由于拟合了正确的成像点位置，所有同相轴都正确成像。图6是本条测线叠前时间偏移成像结果的叠加剖面，上图对应直线方法结果，下图对应本专利技术结果。由于图5中本专利技术成像的CIP道集同相轴完全拉平，沿着横向叠加起来其结果的信噪比明显好于直线近似技术。
需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
基于同一发明构思，本发明实施例还提供了一种高精度弯曲测线叠前时间域成像装置，可以用于实现上述实施例所描述的高精度弯曲测线叠前时间域成像方法，如下面的实施例所述。由于高精度弯曲测线叠前时间域成像装置解决问题的原理与高精度弯曲测线叠前时间域成像方法相似，因此高精度弯曲测线叠前时间域成像装置的实施例可以参见高精度弯曲测线叠前时间域成像方法的实施例，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
在本发明实施例中，本发明的高精度弯曲测线叠前时间域成像装置具体包括：
弯曲测线地震数据获取单元，用于获取弯曲测线地震数据，所述弯曲测线地震数据按照共中心点道集顺序存放，所述弯曲测线地震数据的卷头中包含：最小CMP号以及最大CMP号，所述弯曲测线地震数据中每一个地震道的道头字包含：炮点X坐标、炮点Y坐标、检波点X坐标、检波点Y坐标、中心点X坐标、中心点Y坐标以及CMP号；
炮检距分组单元，用于根据预设炮检距参数组将炮检距分为n组，其中，所述炮检距参数包括：最小炮检距、最大炮检距和炮检距增量；
CMP总数计算单元，用于根据所述最小CMP号和所述最大CMP号计算CMP总数；
三维数组设置单元，用于根据炮检距分组个数n以及所述CMP总数建立三维数组，其中：所述三维数组的维度一的长度等于所述炮检距分组个数n，用于存放炮检距分组信息；所述三维数组的维度二的长度等于6，用于存放CMP号、地震道中点X坐标、地震道中点Y坐标、中心点X坐标、中心点Y坐标以及地震道数；所述三维数组的维度三的长度等于所述CMP总数，用于存放CMP点信息；
三维数组数据生成单元，用于根据所述弯曲测线地震数据中每一个地震道的地震数据生成每一个地震道各自对应的三维数组；
旅行时计算单元，用于根据每一个地震道各自对应的三维数组计算每一个地震道各自对应的旅行时，并根据所述旅行时进行弯曲测线叠前时间偏移。
在本发明一个实施例中，所述三维数组数据生成单元包括：
炮检距计算模块，用于根据地震道的炮点X坐标、炮点Y坐标、检波点X坐标以及检波点Y坐标计算出该地震道的炮检距；
炮检距分组信息确定模块，用于根据该地震道的炮检距确定该地震道的炮检距分组信息，并将该地震道的炮检距分组信息存放在该地震道对应的三维数组的维度一中。
在本发明一个实施例中，所述旅行时计算单元包括：
距离计算模块，用于根据三维数组中的地震道中点X坐标以及三维数组中的地震道中点Y坐标分别计算出炮点到成像点的距离以及检波点到成像点的距离；
旅行时计算模块，用于根据炮点到成像点的距离和检波点到成像点的距离计算出旅行时。
在本发明一个实施例中，所述旅行时计算单元，具体用于在各个地震道在对同一个成像点计算旅行时时，选区各自炮检距的成像点坐标进行计算。
为了实现上述目的，根据本申请的另一方面，还提供了一种计算机设备。如图7所示，该计算机设备包括存储器、处理器、通信接口以及通信总线，在存储器上存储有可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述实施例方法中的步骤。
处理器可以为中央处理器(Central Processing Unit，CPU)。处理器还可以为其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现场可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等芯片，或者上述各类芯片的组合。
存储器作为一种非暂态计算机可读存储介质，可用于存储非暂态软件程序、非暂态计算机可执行程序以及单元，如本发明上述方法实施例中对应的程序单元。处理器通过运行存储在存储器中的非暂态软件程序、指令以及模块，从而执行处理器的各种功能应用以及作品数据处理，即实现上述方法实施例中的方法。
存储器可以包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需要的应用程序；存储数据区可存储处理器所创建的数据等。此外，存储器可以包括高速随机存取存储器，还可以包括非暂态存储器，例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。在一些实施例中，存储器可选包括相对于处理器远程设置的存储器，这些远程存储器可以通过网络连接至处理器。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
所述一个或者多个单元存储在所述存储器中，当被所述处理器执行时，执行上述实施例中的方法。
上述计算机设备具体细节可以对应参阅上述实施例中对应的相关描述和效果进行理解，此处不再赘述。
为了实现上述目的，根据本申请的另一方面，还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序在计算机处理器中执行时实现上述高精度弯曲测线叠前时间域成像方法中的步骤。本领域技术人员可以理解，实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)、随机存储记忆体(RandomAccessMemory，RAM)、快闪存储器(Flash Memory)、硬盘(Hard Disk Drive，缩写：HDD)或固态硬盘(Solid-State Drive，SSD)等；所述存储介质还可以包括上述种类的存储器的组合。
显然，本领域的技术人员应该明白，上述的本发明的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种高精度弯曲测线叠前时间域成像方法，其特征在于，包括：
获取弯曲测线地震数据，所述弯曲测线地震数据按照共中心点道集顺序存放，所述弯曲测线地震数据的卷头中包含：最小CMP号以及最大CMP号，所述弯曲测线地震数据中每一个地震道的道头字包含：炮点X坐标、炮点Y坐标、检波点X坐标、检波点Y坐标、中心点X坐标、中心点Y坐标以及CMP号；
根据预设炮检距参数组将炮检距分为n组，其中，所述炮检距参数包括：最小炮检距、最大炮检距和炮检距增量；
根据所述最小CMP号和所述最大CMP号计算CMP总数；
根据炮检距分组个数n以及所述CMP总数建立三维数组，其中：所述三维数组的维度一的长度等于所述炮检距分组个数n，用于存放炮检距分组信息；所述三维数组的维度二的长度等于6，用于存放CMP号、地震道中点X坐标、地震道中点Y坐标、中心点X坐标、中心点Y坐标以及地震道数；所述三维数组的维度三的长度等于所述CMP总数，用于存放CMP点信息；
根据所述弯曲测线地震数据中每一个地震道的地震数据生成每一个地震道各自对应的三维数组；
根据每一个地震道各自对应的三维数组计算每一个地震道各自对应的旅行时，并根据所述旅行时进行弯曲测线叠前时间偏移。
2.根据权利要求1所述的高精度弯曲测线叠前时间域成像方法，其特征在于，所述根据所述弯曲测线地震数据中每一个地震道的地震数据生成每一个地震道各自对应的三维数组，包括：
根据地震道的炮点X坐标、炮点Y坐标、检波点X坐标以及检波点Y坐标计算出该地震道的炮检距；
根据该地震道的炮检距确定该地震道的炮检距分组信息，并将该地震道的炮检距分组信息存放在该地震道对应的三维数组的维度一中。
3.根据权利要求1所述的高精度弯曲测线叠前时间域成像方法，其特征在于，所述根据每一个地震道各自对应的三维数组计算每一个地震道各自对应的旅行时，包括：
根据三维数组中的地震道中点X坐标以及三维数组中的地震道中点Y坐标分别计算出炮点到成像点的距离以及检波点到成像点的距离；
根据炮点到成像点的距离和检波点到成像点的距离计算出旅行时。
4.根据权利要求1所述的高精度弯曲测线叠前时间域成像方法，其特征在于，所述根据每一个地震道各自对应的三维数组计算每一个地震道各自对应的旅行时，包括：
各个地震道在对同一个成像点计算旅行时时，选区各自炮检距的成像点坐标进行计算。
5.一种高精度弯曲测线叠前时间域成像装置，其特征在于，包括：
弯曲测线地震数据获取单元，用于获取弯曲测线地震数据，所述弯曲测线地震数据按照共中心点道集顺序存放，所述弯曲测线地震数据的卷头中包含：最小CMP号以及最大CMP号，所述弯曲测线地震数据中每一个地震道的道头字包含：炮点X坐标、炮点Y坐标、检波点X坐标、检波点Y坐标、中心点X坐标、中心点Y坐标以及CMP号；
炮检距分组单元，用于根据预设炮检距参数组将炮检距分为n组，其中，所述炮检距参数包括：最小炮检距、最大炮检距和炮检距增量；
CMP总数计算单元，用于根据所述最小CMP号和所述最大CMP号计算CMP总数；
三维数组设置单元，用于根据炮检距分组个数n以及所述CMP总数建立三维数组，其中：所述三维数组的维度一的长度等于所述炮检距分组个数n，用于存放炮检距分组信息；所述三维数组的维度二的长度等于6，用于存放CMP号、地震道中点X坐标、地震道中点Y坐标、中心点X坐标、中心点Y坐标以及地震道数；所述三维数组的维度三的长度等于所述CMP总数，用于存放CMP点信息；
三维数组数据生成单元，用于根据所述弯曲测线地震数据中每一个地震道的地震数据生成每一个地震道各自对应的三维数组；
旅行时计算单元，用于根据每一个地震道各自对应的三维数组计算每一个地震道各自对应的旅行时，并根据所述旅行时进行弯曲测线叠前时间偏移。
6.根据权利要求5所述的高精度弯曲测线叠前时间域成像装置，其特征在于，所述三维数组数据生成单元包括：
炮检距计算模块，用于根据地震道的炮点X坐标、炮点Y坐标、检波点X坐标以及检波点Y坐标计算出该地震道的炮检距；
炮检距分组信息确定模块，用于根据该地震道的炮检距确定该地震道的炮检距分组信息，并将该地震道的炮检距分组信息存放在该地震道对应的三维数组的维度一中。
7.根据权利要求5所述的高精度弯曲测线叠前时间域成像装置，其特征在于，所述旅行时计算单元包括：
距离计算模块，用于根据三维数组中的地震道中点X坐标以及三维数组中的地震道中点Y坐标分别计算出炮点到成像点的距离以及检波点到成像点的距离；
旅行时计算模块，用于根据炮点到成像点的距离和检波点到成像点的距离计算出旅行时。
8.根据权利要求5所述的高精度弯曲测线叠前时间域成像装置，其特征在于，所述旅行时计算单元，具体用于在各个地震道在对同一个成像点计算旅行时时，选区各自炮检距的成像点坐标进行计算。
9.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4任一项所述的方法。
10.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序在计算机处理器中执行时实现如权利要求1至4任意一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011189523.7A CN112379420B (zh) | 2020-10-30 | 2020-10-30 | 高精度弯曲测线叠前时间域成像方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011189523.7A CN112379420B (zh) | 2020-10-30 | 2020-10-30 | 高精度弯曲测线叠前时间域成像方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112379420A true CN112379420A (zh) | 2021-02-19 |
| CN112379420B CN112379420B (zh) | 2024-05-28 |
Family
ID=74576087
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011189523.7A Active CN112379420B (zh) | 2020-10-30 | 2020-10-30 | 高精度弯曲测线叠前时间域成像方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112379420B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101900833A (zh) * | 2010-06-02 | 2010-12-01 | 西安石油大学 | 一种地震散射p-p波成像速度分析方法 |
| CN103698809A (zh) * | 2012-09-28 | 2014-04-02 | 中国石油天然气集团公司 | 一种无加速比瓶颈的克希霍夫叠前时间偏移并行方法 |
| WO2017133361A1 (zh) * | 2016-02-02 | 2017-08-10 | 中国矿业大学(北京) | 一种煤炭陷落柱识别方法和装置 |
| CN108333625A (zh) * | 2018-02-11 | 2018-07-27 | 中国石油化工股份有限公司 | 时频域叠前共成像点道集属性三维分析方法 |
| CN109581494A (zh) * | 2018-10-23 | 2019-04-05 | 中国石油天然气集团有限公司 | 叠前偏移方法及装置 |
-
2020
- 2020-10-30 CN CN202011189523.7A patent/CN112379420B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101900833A (zh) * | 2010-06-02 | 2010-12-01 | 西安石油大学 | 一种地震散射p-p波成像速度分析方法 |
| CN103698809A (zh) * | 2012-09-28 | 2014-04-02 | 中国石油天然气集团公司 | 一种无加速比瓶颈的克希霍夫叠前时间偏移并行方法 |
| WO2017133361A1 (zh) * | 2016-02-02 | 2017-08-10 | 中国矿业大学(北京) | 一种煤炭陷落柱识别方法和装置 |
| CN108333625A (zh) * | 2018-02-11 | 2018-07-27 | 中国石油化工股份有限公司 | 时频域叠前共成像点道集属性三维分析方法 |
| CN109581494A (zh) * | 2018-10-23 | 2019-04-05 | 中国石油天然气集团有限公司 | 叠前偏移方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 王华忠;刘少勇;孔祥宁;蔡杰雄;方伍宝;: "大规模三维地震数据Kirchhoff叠前深度偏移及其并行实现", 石油地球物理勘探, no. 03, 15 June 2012 (2012-06-15) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112379420B (zh) | 2024-05-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3055715B1 (en) | System and method for regularizing seismic data | |
| CN105277978B (zh) | 一种确定近地表速度模型的方法及装置 | |
| WO2017106127A1 (en) | Automated near surface analysis by surface-consistent refraction methods | |
| Inglis et al. | A pipeline for structured light bathymetric mapping | |
| CA2876632C (en) | System and method for migration velocity modeling | |
| WO2014164970A2 (en) | Methods and systems for locating seismic events | |
| WO2014138511A1 (en) | System and method for estimating uncertainty for geophysical gridding routines lacking inherent uncertainty estimation | |
| WO2021127382A1 (en) | Full waveform inversion in the midpoint-offset domain | |
| CN112099082B (zh) | 一种共面元共方位角道集的地震回折波走时反演方法 | |
| NO20170282A1 (en) | Method and apparatus for processing seismic data | |
| US5587942A (en) | 3D wave equation migration of a 2D grid of seismic data | |
| AU739128B2 (en) | A method of seismic processing, and in particular a 3D seismic prospection method implementing seismic data migration | |
| CN109581494B (zh) | 叠前偏移方法及装置 | |
| CN109188527B (zh) | 滩浅海地区快速建立三维近海底速度模型的方法 | |
| CN112379420A (zh) | 高精度弯曲测线叠前时间域成像方法及装置 | |
| RU2739336C2 (ru) | Определение положения сейсмического оборудования в буксируемой морской группе | |
| CN113552635B (zh) | 固源瞬变电磁数据融合隧道三维超前预报方法及系统 | |
| CN104570091B (zh) | 一种获得初至波射线的方法和装置 | |
| CN112465983A (zh) | 一种高程地表模型点线状障碍物修正方法和装置 | |
| CN111624655B (zh) | 初至波剩余静校正量的确定方法及装置 | |
| US10067254B2 (en) | Removal of an estimated acquisition effect from a marine survey measurement | |
| Reshef et al. | Interval Velocity Analysis in Complex Areas–Azimuth Considerations | |
| Zhou et al. | Automation of MBES noise reduction: An approach based on seafloor bathymetry features derived from manual editing procedures | |
| EP3004940B1 (en) | Method for enhancing the resolution of a seismic image | |
| CN115184987A (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 |