CN116699684A - 面向储层预测的地震数据拼接方法、装置、设备及介质 - Google Patents
面向储层预测的地震数据拼接方法、装置、设备及介质 Download PDFInfo
- Publication number
- CN116699684A CN116699684A CN202310686496.1A CN202310686496A CN116699684A CN 116699684 A CN116699684 A CN 116699684A CN 202310686496 A CN202310686496 A CN 202310686496A CN 116699684 A CN116699684 A CN 116699684A
- Authority
- CN
- China
- Prior art keywords
- angle
- gathers
- range
- intersection
- work areas
- 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
- 238000000034 method Methods 0.000 title claims abstract description 65
- 238000012545 processing Methods 0.000 claims abstract description 24
- 238000004590 computer program Methods 0.000 claims description 16
- 238000003860 storage Methods 0.000 claims description 13
- 238000007781 pre-processing Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 230000001629 suppression Effects 0.000 claims description 5
- 238000004891 communication Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 230000009286 beneficial effect Effects 0.000 description 4
- 238000005520 cutting process Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000012958 reprocessing Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000001953 sensory effect Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000000638 solvent extraction Methods 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
-
- 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/282—Application of seismic models, synthetic seismograms
-
- 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/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
-
- 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)
Abstract
本申请实施例公开了一种地震数据拼接方法、装置、电子设备及介质。该方法包括：根据至少两个工区的共反射点道集以及时间域层速度体确定至少两个工区的角道集；根据待反演时窗内至少两个工区中目标层位的角道集的特征信息确定至少两个工区的角道集角度范围，确定各角度范围的交集范围；根据目标工区中目标层位的角道集的特征信息，对交集范围进行划分得到交集子范围；根据交集子范围对至少两个工区的角道集进行处理得到拼接数据体。本申请实施例的技术方案，可以各工区内的角道集按照交集子范围进行叠加得到分角度的叠加体，再对分角度的叠加体进行拼接，减少需要进行处理的数据量，进而缩短地震数据拼接的周期，以便于快速的进行叠前储层预测。
Description
技术领域
本申请涉及石油勘探开发领域，尤其涉及一种地震数据拼接方法、装置、电子设备及介质。
背景技术
随着油田勘探的不断深入，许多主力研究区域已完成三维全覆盖，但是三维地震调查都是分区分块进行的，这样导致资料由于采集条件及处理技术的不同，在成像、能量、相位、频率上均存在明显的不一致性，不利于对研究区开展整体评价，为进行整体评价，需要对地震数据进行拼接处理。
目前常用的拼接方法主要有叠后拼接和叠前重新处理拼接两种方法。叠后拼接成果没有叠前信息，多用于进行构造解释，不能进行叠前储层预测；叠前重新处理拼接涉及叠前大面积、多区块叠前连片重新处理，成本很高、周期很长，不利于快速进行叠前储层预测。
发明内容
本申请提供了一种地震数据拼接方法、装置、电子设备及介质，通过精确的定位各工区拼接所要用到的角道集，减少需要进行处理的数据量，以实现缩短地震数据拼接周期，便于快速进行叠前储层预测的效果。
根据本申请的一方面，提供了一种地震数据拼接方法，所述方法包括：
根据至少两个工区的共反射点道集以及时间域层速度体，确定至少两个工区的角道集；
根据待反演时窗内至少两个工区中目标层位的角道集的特征信息，确定至少两个工区的角道集角度范围，并确定各角度范围的交集范围；
根据待反演时窗内目标工区中目标层位的角道集的特征信息，对所述交集范围进行划分，得到交集子范围；
根据交集子范围对至少两个工区的角道集进行处理，得到拼接数据体。
根据本申请的另一方面，提供了一种地震数据拼接装置，所述装置包括：
角道集确定模块，用于根据至少两个工区的共反射点道集以及时间域层速度体，确定至少两个工区的角道集；
交集范围确定模块，用于根据待反演时窗内至少两个工区中目标层位的角道集的特征信息，确定至少两个工区的角道集角度范围，并确定各角度范围的交集范围；
交集子范围确定模块，用于根据待反演时窗内目标工区中目标层位的角道集的特征信息，对所述交集范围进行划分，得到交集子范围；
拼接数据体确定模块，用于根据交集子范围对至少两个工区的角道集进行处理，得到拼接数据体。
根据本申请的另一方面，提供了一种电子设备，该设备包括：
至少一个处理器；以及
与至少一个处理器通信连接的存储器；其中，
存储器存储有可被至少一个处理器执行的计算机程序，计算机程序被至少一个处理器执行，以使至少一个处理器能够执行本申请任一实施例的地震数据拼接方法。
根据本申请的另一方面，提供了一种计算机可读存储介质，计算机可读存储介质存储有计算机指令，计算机指令用于使处理器执行时实现本申请任一实施例的地震数据拼接方法。
本申请实施例提供了一种地震数据拼接方法，根据至少两个工区的共反射点道集以及时间域层速度体，确定至少两个工区的角道集，可以让人更加直观的了解各工区在各个角度的地震数据的情况；根据待反演时窗内至少两个工区中目标层位的角道集的特征信息，确定至少两个工区的角道集角度范围，并确定各角度范围的交集范围，可以对不属于各角度范围的交集范围之内的角道集的数据进行排除，减小需要进行处理的数据量；根据待反演时窗内目标工区中目标层位的角道集的特征信息，对交集范围进行划分，得到交集子范围，可以进一步缩小需要进行处理的角道集的数据的范围，减小工作量；根据交集子范围对至少两个工区的角道集进行处理，得到拼接数据体。本申请实施例的技术方案，可以精确的定位到各工区拼接所要用到的道集数据，减少需要进行处理的数据量，进而缩短地震数据拼接的周期，以便于快速的进行叠前储层预测。
应当理解，本部分所描述的内容并非旨在标识本申请的实施例的关键或重要特征，也不用于限制本申请的范围。本申请的其它特征将通过以下的说明书而变得容易理解。
附图说明
为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是根据本申请实施例一提供的一种地震数据拼接方法的流程图；
图2是根据本申请实施例二提供的一种地震数据拼接方法的流程图；
图3是根据本申请实施例二提供的一种匹配算子示意图；
图4是根据本申请实施例三提供的一种地震数据拼接方法的流程图；
图5是根据本申请实施例四提供的一种地震数据拼接方法的具体实现流程图；
图6是根据本申请实施例四提供的一种地震数据拼接装置的结构示意图；
图7是实现本申请实施例五提供的一种地震数据拼接方法的电子设备的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分的实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
需要说明的是，本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”、“第三”、“第四”、“实际”、“预设”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本申请的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
实施例一
图1为本申请实施例一提供的一种地震数据拼接方法的流程图，本申请实施例可适用于对叠前储层进行预测的情况。该方法可以由地震数据拼接装置来执行，该地震数据拼接装置可以采用硬件和/或软件的形式实现，该地震数据拼接装置可配置于电子设备中。如图1所示，该方法包括：
S110、根据至少两个工区的共反射点道集以及时间域层速度体，确定至少两个工区的角道集。
其中，共反射点道集指的是同一个地下反射点反射回来的地震记录的集合。时间域层速度体是地球物理学中用于描述地下结构的术语，通常用于勘探和采集地震数据，它是一个三维模型，表示地下不同深度内的地质属性，如波速、密度、衰减、反射率。
在本申请实施例中，为了便于观测和分析地震反射振幅随入射角的变化，可以根据至少两个工区的共反射点道集以及时间域层速度体，将固定炮检距的道集记录转换成固定入射角的道集记录，即确定至少两个工区的角道集。由此，可以让人更加直观的了解各工区在各个角度的地震数据的情况。
可选地，在根据至少两个工区的共反射点道集以及时间域层速度体，确定至少两个工区的角道集之前，该方法还包括：
对至少两个工区的共反射点道集进行预处理；其中，预处理包括畸变切除、剩余多次波压制以及道集拉平中的至少一项。
在本申请实施例中，由于采集条件的差异，得到的共反射点道集中存在有各种问题，因此，在进行角道集的确定之前，需要对采集到的共反射点道集进行预处理。其中，预处理包括以下至少一项：1、畸变切除，即对共反射点道集中发生畸变的部分进行切除处理；2、剩余多次波压制，多次波的存在严重影响到地震剖面的分辨率和信噪比，有效的进行剩余多次波压制可提高剖面的成像品质；3、道集拉平，AVO分析需要平的道集，但很多时候使用解析曲线无法达到这个目标，因此，需要先对共反射点道集进行拉平处理。对共反射点道集进行预处理，可以确保数据的准确性和一致性，从而提高后续处理流程的质量。
S120、根据待反演时窗内至少两个工区中目标层位的角道集的特征信息，确定至少两个工区的角道集角度范围，并确定各角度范围的交集范围。
其中，待反演时窗即为需要进行叠前储层反演的时段，目标层位即为需要进行叠前储层反演的位置。
在本申请实施例中，在确定了待反演时窗和目标层位之后，可以根据待反演时窗内至少两个工区中目标层位的角道集的特征信息，确定至少两个工区的角道集角度范围，并选择至少两个工区中最大的最小角度与至少两个工区中最小的最大角度之间的角度范围作为角度保留范围，即确定各角度范围的交集范围。由此，可以对不属于各角度范围的交集范围之内的角道集进行排除，减小需要进行处理的数据量。
S130、根据待反演时窗内目标工区中目标层位的角道集的特征信息，对交集范围进行划分，得到交集子范围。
在本申请实施例中，在确定了各角度范围的交集范围之后，可以根据待反演时窗内目标工区中目标层位的角道集的特征信息，对交集范围进行划分，得到交集子范围。示例性地，可以根据待反演时窗内目标工区A中目标层位的角道集的波形或者频率，将目标工区A与其他工区的角道集的角度交集范围划分为交集子范围。可以理解的是，在进行交集范围的划分时，可以对交集范围进行平均划分，也可以不进行平均划分，可以进一步缩小需要进行处理的数据的范围，减小工作量。
S140、根据交集子范围对至少两个工区的角道集进行处理，得到拼接数据体。
在本申请实施例中，对交集范围进行划分，得到交集子范围之后，即可根据交集子范围对至少两个工区的属于交集子范围内的角道集进行拼接，得到拼接数据体，该拼接数据体可用于叠前储层预测。
本申请实施例提供了一种地震数据拼接方法，根据至少两个工区的共反射点道集以及时间域层速度体，确定至少两个工区的角道集，可以让人更加直观的了解各工区在各个角度的地震数据的情况；根据待反演时窗内至少两个工区中目标层位的角道集的特征信息，确定至少两个工区的角道集角度范围，并确定各角度范围的交集范围，可以对不属于各角度范围的交集范围之内的角道集进行排除，减小需要进行处理的数据量；根据待反演时窗内目标工区中目标层位的角道集的特征信息，对交集范围进行划分，得到交集子范围，可以进一步缩小需要进行处理的数据的范围，减小工作量；根据交集子范围对至少两个工区的角道集进行处理，得到拼接数据体。本申请实施例的技术方案，可以各工区内的角道集按照交集子范围进行叠加，减少需要进行处理的数据量，进而缩短地震数据拼接的周期，以便于快速的进行叠前储层预测。
实施例二
图2为本申请实施例二提供的一种地震数据拼接方法的流程图，本申请实施例以上述实施例为基础进行优化，未在本申请实施例中详尽描述的方案见上述实施例。如图2所示，本申请实施例的方法具体包括如下步骤：
S210、根据至少两个工区的共反射点道集以及时间域层速度体，确定至少两个工区的角道集。
S220、根据待反演时窗内至少两个工区中目标层位的角道集的特征信息，确定至少两个工区的角道集角度范围，并确定各角度范围的交集范围。
S230、将交集范围划分为至少三个交集子范围；其中，各交集子范围中至少两个角度对应角道集的特征信息满足特征一致性条件。
在本申请实施例中，在对交集范围进行划分时，需要将交集范围划分为至少三个交集子范围。可以理解的是，划分出来的交集子范围越多，在拼接时用于锚定拼接位置的点也就越多，得到的拼接结果也就越准确。其中，各交集子范围中至少两个角度对应角道集的特征信息需要满足特征一致性条件。
具体地，特征信息包括频率，确定至少两个角度对应角道集的特征信息是否满足特征一致性条件的过程包括：
若至少两个角度对应角道集的频率的差值小于预设差值，则确定至少两个角度对应角道集的特征信息满足特征一致性条件；
否则，则确定至少两个角度对应角道集的特征信息不满足特征一致性条件。
示例性地，一个交集子范围中包括五个角度，若其中有两个或两个以上角度对应的角道集的频率的差值小于预设差值，则确定该交集子范围中至少两个角度对应角道集的特征信息满足特征一致性条件，否则，则确定该交集子范围中至少两个角度对应角道集的特征信息不满足特征一致性条件。
S240、针对各工区的角道集，将属于同一交集子范围中的角度对应的角道集进行叠加，得到叠加数据体。
在本申请实施例中，可以针对各工区的角道集，将属于同一交集子范围中的角度对应的角道集进行叠加，得到叠加数据体。示例性地，在通过步骤S230的划分方案对N个待拼接工区进行同一交集子范围中的角度对应的角道集进行叠加时，每个待拼接工区均得到至少3个叠加数据体，最终得到的总叠加数据体为至少3*N个。
S250、将同一交集子范围对应的各工区的叠加数据体进行拼接，得到拼接数据体。
在本申请实施例中，将同一交集子范围对应的各工区的叠加数据体进行拼接，即可得到用于叠前储层预测的拼接数据体。
可选地，将同一交集子范围对应的各工区的叠加数据体进行拼接之前，所述方法还包括：
根据至少两个工区中的目标工区的角道集以及其他工区的角道集确定匹配算子；
基于匹配算子对其他工区的角道集进行处理，以使其他工区的角道集与目标工区的角道集相匹配。
图3示出了一种匹配算子的示意图，在本申请实施例中，可以先将至少两个工区中的目标工区的角道集以及其他工区的角道集通过傅里叶变换，得到目标工区和其他工区在频域内的角道集，并根据目标工区和其他工区在频域内的角道集计算匹配算子。再将目标工区和其他工区在频域内的角道集进行反傅里叶变换，得到原来的目标工区和其他工区的角道集。最后，基于匹配算子对其他工区的角道集进行处理，以使其他工区的角道集与目标工区的角道集相匹配。
本申请实施例提供了一种地震数据拼接方法，根据至少两个工区的共反射点道集以及时间域层速度体，确定至少两个工区的角道集，可以让人更加直观的了解各工区在各个角度的地震数据的情况；根据待反演时窗内至少两个工区中目标层位的角道集的特征信息，确定至少两个工区的角道集角度范围，并确定各角度范围的交集范围，可以对不属于各角度范围的交集范围之内的角道集的数据进行排除，减小需要进行处理的数据量；将所述交集范围划分为至少三个交集子范围，可以提高得到的拼接结果的准确性；其中，各交集子范围中至少两个角度对应角道集的特征信息满足特征一致性条件；针对各工区的角道集，将属于同一交集子范围中的角度对应的角道集进行叠加，得到叠加数据体；将同一交集子范围对应的各工区的叠加数据体进行拼接，得到拼接数据体。本申请实施例的技术方案，可以精确的定位到各工区拼接所要用到的道集数据，减少需要进行处理的数据量，进而缩短地震数据拼接的周期，以便于快速的进行叠前储层预测，同时，提高了拼接结果的准确性。
实施例三
图4为本申请实施例三提供的一种地震数据拼接方法的流程图，本申请实施例以上述实施例为基础进行优化，未在本申请实施例中详尽描述的方案见上述实施例。如图4所示，本申请实施例的方法具体包括如下步骤：
S310、根据至少两个工区的共反射点道集以及时间域层速度体，确定至少两个工区的角道集。
S320、根据待反演时窗内至少两个工区中目标层位的角道集的特征信息，确定至少两个工区的角道集角度范围，并确定各角度范围的交集范围。
S330、若各角度对应角道集的特征信息满足特征一致性条件，且各角度连续，则将连续的角度对应的范围作为一个交集子范围。
在本申请实施例中，若该角度对应角道集的频率的差值小于预设差值，则确定该角度对应角道集的特征信息满足特征一致性条件；否则，则确定该角度对应角道集的特征信息不满足特征一致性条件。若存在连续的满足特征一致性条件的角度，则将连续的角度对应的范围作为一个交集子范围。示例性地，在一个交集范围中，如果角度2-10对应角道集的特征信息都满足特征一致性条件，则将角度2-10作为该交集范围的一个交集子范围。
S340、将其他角道集的特征信息对应的连续角度构成的范围作为一个交集子范围。
示例性地，在一个交集范围中，在确定了角度2-10作为该交集范围的一个交集子范围之后，还需要将其他角道集的特征信息对应的连续角度构成的范围作为一个交集子范围，如角度15-20对应角道集的特征信息也都满足特征一致性条件，则将角度15-20也作为该交集范围的一个交集子范围。
S350、针对各工区的角道集，将属于同一交集子范围中的角度对应的角道集进行叠加，得到叠加数据体。
S360、将同一交集子范围对应的各工区的叠加数据体进行拼接，得到拼接数据体。
本申请实施例提供了一种地震数据拼接方法，根据至少两个工区的共反射点道集以及时间域层速度体，确定至少两个工区的角道集，可以让人更加直观的了解各工区在各个角度的地震数据的情况；根据待反演时窗内至少两个工区中目标层位的角道集的特征信息，确定至少两个工区的角道集角度范围，并确定各角度范围的交集范围，可以对不属于各角度范围的交集范围之内的角道集的数据进行排除，减小需要进行处理的数据量；若各角度对应角道集的特征信息满足特征一致性条件，且各角度连续，则将连续的角度对应的范围作为一个交集子范围；将其他角道集的特征信息对应的连续角度构成的范围作为一个交集子范围；针对各工区的角道集，将属于同一交集子范围中的角度对应的角道集进行叠加，得到叠加数据体；将同一交集子范围对应的各工区的叠加数据体进行拼接，得到拼接数据体。本申请实施例的技术方案，可以精确的定位到各工区拼接所要用到的道集数据，减少需要进行处理的数据量，进而缩短地震数据拼接的周期，以便于快速的进行叠前储层预测。
实施例四
图5为本申请实施例四提供的一种地震数据拼接方法的具体实现流程图，在上述实施例的基础上，提供了一种优选实施例，具体如下：
S401、准备N个工区的共反射点道集、时间域层速度体、以及叠前反演目的层位。
S402、根据需要分别对N个工区的共反射点道集进行畸变切除、剩余多次波压制、道集拉平等道集优化处理。
S403、确定N个工区的目标工区A，其他工区的振幅、频率以及相位特征均往目标工区A匹配。
S404、分别利用N个工区的时间域层速度体和共反射点道集进行转换，获得N个工区的角道集。
S405、根据反演时窗需要，分析N个工区中目的层位的角道集的波形特征、频率特征以及AVO特征，确定N个工区的角道集角度保留的最小值和最大值。
S406、综合分析各工区的角度保留范围，各工区角度保留范围一致，选择N个工区中最大的最小角度和最小的最大角度作为角度保留范围。以目标工区A的待反演时窗的道集波形、频率以及目的层AVO特征确定角度划分方案，至少3个角度。
S407、以目标工区A的角度划分方案(至少3个角度)对N个工区进行分角度叠加。每个工区均得到至少3个叠加数据体，最终得到的总叠加数据体为至少3*N个。
S408、以目标工区A的叠加数据体(至少3个叠加数据体)作为目标，以各工区重叠部分数据为基础，进行叠加数据体匹配处理，包括振幅、相位以及频率的匹配。
S409、各工区匹配完成以后，进行拼接处理，最终得到至少3个拼接数据体。
S410、基于拼接数据体，可进行叠前储层预测，用于研究区域整体储层评价。
本申请实施例具有上述实施例相同的有益效果。
实施例五
图6为本申请实施例四提供的一种地震数据拼接装置的结构示意图，该装置可执行本申请任意实施例所提供的地震数据拼接方法，具备执行方法相应的功能模块和有益效果。如图6所示，所述装置包括：
角道集确定模块510，用于根据至少两个工区的共反射点道集以及时间域层速度体，确定至少两个工区的角道集；
交集范围确定模块520，用于根据待反演时窗内至少两个工区中目标层位的角道集的特征信息，确定至少两个工区的角道集角度范围，并确定各角度范围的交集范围；
交集子范围确定模块530，用于根据待反演时窗内目标工区中目标层位的角道集的特征信息，对所述交集范围进行划分，得到交集子范围；
拼接数据体确定模块540，用于根据交集子范围对至少两个工区的角道集进行处理，得到拼接数据体。
在本申请实施例中，交集子范围确定模块530，包括：
第一交集子范围确定单元，用于将所述交集范围划分为至少三个交集子范围；其中，各交集子范围中至少两个角度对应角道集的特征信息满足特征一致性条件。
在本申请实施例中，交集子范围确定模块530，包括：
第二交集子范围确定单元，第二交集子范围确定单元，包括：
第一交集子范围确定子单元，用于若各角度对应角道集的特征信息满足特征一致性条件，且各角度连续，则将连续的角度对应的范围作为一个交集子范围；
第二交集子范围确定子单元，用于将其他角道集的特征信息对应的连续角度构成的范围作为一个交集子范围。
可选地，所述特征信息包括频率；第一交集子范围确定单元，包括：
特征一致性确定子单元，具体用于：
若至少两个角度对应角道集的频率的差值小于预设差值，则确定至少两个角度对应角道集的特征信息满足特征一致性条件；
否则，则确定至少两个角度对应角道集的特征信息不满足特征一致性条件。
在本申请实施例中，拼接数据体确定模块540，包括：
叠加数据体确定单元，用于针对各工区的角道集，将属于同一交集子范围中的角度对应的角道集进行叠加，得到叠加数据体；
拼接数据体确定单元，用于将同一交集子范围对应的各工区的叠加数据体进行拼接，得到拼接数据体。
在本申请实施例中，将同一交集子范围对应的各工区的叠加数据体进行拼接之前，拼接数据体确定模块540，还包括：
匹配算子确定单元，用于根据至少两个工区中的目标工区的角道集以及其他工区的角道集确定匹配算子；
角道集匹配单元，用于基于匹配算子对其他工区的角道集进行处理，以使其他工区的角道集与目标工区的角道集相匹配。
可选地，在本申请实施例中，根据至少两个工区的共反射点道集以及时间域层速度体，确定至少两个工区的角道集之前，所述装置，还包括：
预处理模块，用于对至少两个工区的共反射点道集进行预处理；其中，预处理包括畸变切除、剩余多次波压制以及道集拉平中的至少一项。
本申请实施例所提供的一种地震数据拼接装置可执行本申请任意实施例所提供的一种地震数据拼接方法，具备执行方法相应的功能模块和有益效果。
实施例五
图7示出了可以用来实施本申请的实施例的电子设备10的结构示意图。电子设备旨在表示各种形式的数字计算机，诸如，膝上型计算机、台式计算机、工作台、个人数字助理、服务器、刀片式服务器、大型计算机、和其它适合的计算机。电子设备还可以表示各种形式的移动装置，诸如，个人数字处理、蜂窝电话、智能电话、可穿戴设备(如头盔、眼镜、手表等)和其它类似的计算装置。本文所示的部件、它们的连接和关系、以及它们的功能仅仅作为示例，并且不意在限制本文中描述的和/或者要求的本申请的实现。
如图7所示，电子设备10包括至少一个处理器11，以及与至少一个处理器11通信连接的存储器，如只读存储器(ROM)12、随机访问存储器(RAM)13等，其中，存储器存储有可被至少一个处理器执行的计算机程序，处理器11可以根据存储在只读存储器(ROM)12中的计算机程序或者从存储单元18加载到随机访问存储器(RAM)13中的计算机程序，来执行各种适当的动作和处理。在RAM 13中，还可存储电子设备10操作所需的各种程序和数据。处理器11、ROM 12以及RAM 13通过总线14彼此相连。输入/输出(I/O)接口15也连接至总线14。
电子设备10中的多个部件连接至I/O接口15，包括：输入单元16，例如键盘、鼠标等；输出单元17，例如各种类型的显示器、扬声器等；存储单元18，例如磁盘、光盘等；以及通信单元19，例如网卡、调制解调器、无线通信收发机等。通信单元19允许电子设备10通过诸如因特网的计算机网络和/或各种电信网络与其他设备交换信息/数据。
处理器11可以是各种具有处理和计算能力的通用和/或专用处理组件。处理器11的一些示例包括但不限于中央处理单元(CPU)、图形处理单元(GPU)、各种专用的人工智能(AI)计算芯片、各种运行机器学习模型算法的处理器、数字信号处理器(DSP)、以及任何适当的处理器、控制器、微控制器等。处理器11执行上文所描述的各个方法和处理，例如地震数据拼接方法。
在一些实施例中，地震数据拼接方法可被实现为计算机程序，其被有形地包含于计算机可读存储介质，例如存储单元18。在一些实施例中，计算机程序的部分或者全部可以经由ROM 12和/或通信单元19而被载入和/或安装到电子设备10上。当计算机程序加载到RAM 13并由处理器11执行时，可以执行上文描述的地震数据拼接方法的一个或多个步骤。备选地，在其他实施例中，处理器11可以通过其他任何适当的方式(例如，借助于固件)而被配置为执行地震数据拼接方法。
本文中以上描述的系统和技术的各种实施方式可以在数字电子电路系统、集成电路系统、场可编程门阵列(FPGA)、专用集成电路(ASIC)、专用标准产品(ASSP)、芯片上系统的系统(SOC)、负载可编程逻辑设备(CPLD)、计算机硬件、固件、软件、和/或它们的组合中实现。这些各种实施方式可以包括：实施在一个或者多个计算机程序中，该一个或者多个计算机程序可在包括至少一个可编程处理器的可编程系统上执行和/或解释，该可编程处理器可以是专用或者通用可编程处理器，可以从存储系统、至少一个输入装置、和至少一个输出装置接收数据和指令，并且将数据和指令传输至该存储系统、该至少一个输入装置、和该至少一个输出装置。
用于实施本申请的方法的计算机程序可以采用一个或多个编程语言的任何组合来编写。这些计算机程序可以提供给通用计算机、专用计算机或其他可编程地震数据拼接装置的处理器，使得计算机程序当由处理器执行时使流程图和/或框图中所规定的功能/操作被实施。计算机程序可以完全在机器上执行、部分地在机器上执行，作为独立软件包部分地在机器上执行且部分地在远程机器上执行或完全在远程机器或服务器上执行。
在本申请的上下文中，计算机可读存储介质可以是有形的介质，其可以包含或存储以供指令执行系统、装置或设备使用或与指令执行系统、装置或设备结合地使用的计算机程序。计算机可读存储介质可以包括但不限于电子的、磁性的、光学的、电磁的、红外的、或半导体系统、装置或设备，或者上述内容的任何合适组合。备选地，计算机可读存储介质可以是机器可读信号介质。机器可读存储介质的更具体示例会包括基于一个或多个线的电气连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦除可编程只读存储器(EPROM或快闪存储器)、光纤、便捷式紧凑盘只读存储器(CD-ROM)、光学储存设备、磁储存设备、或上述内容的任何合适组合。
为了提供与用户的交互，可以在电子设备上实施此处描述的系统和技术，该电子设备具有：用于向用户显示信息的显示装置(例如，CRT(阴极射线管)或者LCD(液晶显示器)监视器)；以及键盘和指向装置(例如，鼠标或者轨迹球)，用户可以通过该键盘和该指向装置来将输入提供给电子设备。其它种类的装置还可以用于提供与用户的交互；例如，提供给用户的反馈可以是任何形式的传感反馈(例如，视觉反馈、听觉反馈、或者触觉反馈)；并且可以用任何形式(包括声输入、语音输入或者、触觉输入)来接收来自用户的输入。
可以将此处描述的系统和技术实施在包括后台部件的计算系统(例如，作为数据服务器)、或者包括中间件部件的计算系统(例如，应用服务器)、或者包括前端部件的计算系统(例如，具有图形用户界面或者网络浏览器的用户计算机，用户可以通过该图形用户界面或者该网络浏览器来与此处描述的系统和技术的实施方式交互)、或者包括这种后台部件、中间件部件、或者前端部件的任何组合的计算系统中。可以通过任何形式或者介质的数字数据通信(例如，通信网络)来将系统的部件相互连接。通信网络的示例包括：局域网(LAN)、广域网(WAN)、区块链网络和互联网。
计算系统可以包括客户端和服务器。客户端和服务器一般远离彼此并且通常通过通信网络进行交互。通过在相应的计算机上运行并且彼此具有客户端-服务器关系的计算机程序来产生客户端和服务器的关系。服务器可以是云服务器，又称为云计算服务器或云主机，是云计算服务体系中的一项主机产品，以解决了传统物理主机与VPS服务中，存在的管理难度大，业务扩展性弱的缺陷。
应该理解，可以使用上面所示的各种形式的流程，重新排序、增加或删除步骤。例如，本申请中记载的各步骤可以并行地执行也可以顺序地执行也可以不同的次序执行，只要能够实现本申请的技术方案所期望的信息，本文在此不进行限制。
上述具体实施方式，并不构成对本申请保护范围的限制。本领域技术人员应该明白的是，根据设计要求和其他因素，可以进行各种修改、组合、子组合和替代。任何在本申请的精神和原则之内所作的修改、等同替换和改进等，均应包含在本申请保护范围之内。
Claims (10)
1.一种地震数据拼接方法，其特征在于，所述方法包括：
根据至少两个工区的共反射点道集以及时间域层速度体，确定至少两个工区的角道集；
根据待反演时窗内至少两个工区中目标层位的角道集的特征信息，确定至少两个工区的角道集角度范围，并确定各角度范围的交集范围；
根据待反演时窗内目标工区中目标层位的角道集的特征信息，对所述交集范围进行划分，得到交集子范围；
根据交集子范围对至少两个工区的角道集进行处理，得到拼接数据体。
2.根据权利要求1所述的方法，其特征在于，根据待反演时窗内目标工区中目标层位的角道集的特征信息，对所述交集范围进行划分，得到交集子范围，包括：
将所述交集范围划分为至少三个交集子范围；其中，各交集子范围中至少两个角度对应角道集的特征信息满足特征一致性条件。
3.根据权利要求1所述的方法，其特征在于，根据待反演时窗内目标工区中目标层位的角道集的特征信息，对所述交集范围进行划分，得到交集子范围，包括：
若各角度对应角道集的特征信息满足特征一致性条件，且各角度连续，则将连续的角度对应的范围作为一个交集子范围；
将其他角道集的特征信息对应的连续角度构成的范围作为一个交集子范围。
4.根据权利要求1所述的方法，其特征在于，所述特征信息包括频率；
确定至少两个角度对应角道集的特征信息是否满足特征一致性条件的过程包括：
若至少两个角度对应角道集的频率的差值小于预设差值，则确定至少两个角度对应角道集的特征信息满足特征一致性条件；
否则，则确定至少两个角度对应角道集的特征信息不满足特征一致性条件。
5.根据权利要求1所述的方法，其特征在于，根据交集子范围对至少两个工区的角道集进行处理，得到拼接数据体，包括：
针对各工区的角道集，将属于同一交集子范围中的角度对应的角道集进行叠加，得到叠加数据体；
将同一交集子范围对应的各工区的叠加数据体进行拼接，得到拼接数据体。
6.根据权利要求5所述的方法，其特征在于，将同一交集子范围对应的各工区的叠加数据体进行拼接之前，所述方法还包括：
根据至少两个工区中的目标工区的角道集以及其他工区的角道集确定匹配算子；
基于匹配算子对其他工区的角道集进行处理，以使其他工区的角道集与目标工区的角道集相匹配。
7.根据权利要求1所述的方法，其特征在于，根据至少两个工区的共反射点道集以及时间域层速度体，确定至少两个工区的角道集之前，所述方法还包括：
对至少两个工区的共反射点道集进行预处理；其中，预处理包括畸变切除、剩余多次波压制以及道集拉平中的至少一项。
8.一种地震数据拼接装置，其特征在于，所述装置包括：
角道集确定模块，用于根据至少两个工区的共反射点道集以及时间域层速度体，确定至少两个工区的角道集；
交集范围确定模块，用于根据待反演时窗内至少两个工区中目标层位的角道集的特征信息，确定至少两个工区的角道集角度范围，并确定各角度范围的交集范围；
交集子范围确定模块，用于根据待反演时窗内目标工区中目标层位的角道集的特征信息，对所述交集范围进行划分，得到交集子范围；
拼接数据体确定模块，用于根据交集子范围对至少两个工区的角道集进行处理，得到拼接数据体。
9.一种电子设备，其特征在于，所述设备包括：
至少一个处理器；以及
与所述至少一个处理器通信连接的存储器；其中，
所述存储器存储有可被所述至少一个处理器执行的计算机程序，所述计算机程序被所述至少一个处理器执行，以使所述至少一个处理器能够执行权利要求1-7中任一项所述的地震数据拼接方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使处理器执行时实现权利要求1-7中任一项所述的地震数据拼接方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310686496.1A CN116699684A (zh) | 2023-06-09 | 2023-06-09 | 面向储层预测的地震数据拼接方法、装置、设备及介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310686496.1A CN116699684A (zh) | 2023-06-09 | 2023-06-09 | 面向储层预测的地震数据拼接方法、装置、设备及介质 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116699684A true CN116699684A (zh) | 2023-09-05 |
Family
ID=87836942
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310686496.1A Pending CN116699684A (zh) | 2023-06-09 | 2023-06-09 | 面向储层预测的地震数据拼接方法、装置、设备及介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116699684A (zh) |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010047858A1 (en) * | 2008-10-20 | 2010-04-29 | Exxonmobil Upstream Research Company | Method for predicting time-lapse seismic timeshifts by computer simulation |
| CN105510965A (zh) * | 2015-12-29 | 2016-04-20 | 中国石油天然气股份有限公司 | 流体识别方法和装置 |
| WO2016083892A2 (en) * | 2014-11-25 | 2016-06-02 | Cgg Services Sa | Estimating a time variant signal representing a seismic source |
| CN107193040A (zh) * | 2017-06-27 | 2017-09-22 | 中国石油天然气股份有限公司 | 深度域合成地震记录的确定方法和装置 |
| CN108693557A (zh) * | 2018-04-27 | 2018-10-23 | 中国石油天然气集团有限公司 | 一种地震数据处理方法及装置 |
| CN112305594A (zh) * | 2019-07-23 | 2021-02-02 | 中国石油天然气股份有限公司 | 非均质储层的油气分布确定方法及系统 |
| CN112379438A (zh) * | 2020-11-02 | 2021-02-19 | 中国石油天然气集团有限公司 | 优化入射角部分叠加地震数据ava特征的方法及装置 |
| CN112817042A (zh) * | 2019-11-15 | 2021-05-18 | 中国石油天然气股份有限公司 | 沉积储层中油气识别方法及装置 |
| CN113933901A (zh) * | 2020-06-29 | 2022-01-14 | 中国石油天然气股份有限公司 | 深度域地震数据拼接方法及装置 |
| US20220043175A1 (en) * | 2020-08-07 | 2022-02-10 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Method and Device for Imaging Diffracted Waves Based on Azimuth-dip Angle Gathers, and Storage Medium |
| CN114924314A (zh) * | 2022-06-24 | 2022-08-19 | 中国石油大学(北京) | 地震剪切横波反演模型构建方法、孔隙度预测方法及装置 |
| CN115267912A (zh) * | 2022-08-02 | 2022-11-01 | 西南石油大学 | 基于压缩感知的avo特征保持去噪方法及系统 |
-
2023
- 2023-06-09 CN CN202310686496.1A patent/CN116699684A/zh active Pending
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010047858A1 (en) * | 2008-10-20 | 2010-04-29 | Exxonmobil Upstream Research Company | Method for predicting time-lapse seismic timeshifts by computer simulation |
| WO2016083892A2 (en) * | 2014-11-25 | 2016-06-02 | Cgg Services Sa | Estimating a time variant signal representing a seismic source |
| CN105510965A (zh) * | 2015-12-29 | 2016-04-20 | 中国石油天然气股份有限公司 | 流体识别方法和装置 |
| CN107193040A (zh) * | 2017-06-27 | 2017-09-22 | 中国石油天然气股份有限公司 | 深度域合成地震记录的确定方法和装置 |
| CN108693557A (zh) * | 2018-04-27 | 2018-10-23 | 中国石油天然气集团有限公司 | 一种地震数据处理方法及装置 |
| CN112305594A (zh) * | 2019-07-23 | 2021-02-02 | 中国石油天然气股份有限公司 | 非均质储层的油气分布确定方法及系统 |
| CN112817042A (zh) * | 2019-11-15 | 2021-05-18 | 中国石油天然气股份有限公司 | 沉积储层中油气识别方法及装置 |
| CN113933901A (zh) * | 2020-06-29 | 2022-01-14 | 中国石油天然气股份有限公司 | 深度域地震数据拼接方法及装置 |
| US20220043175A1 (en) * | 2020-08-07 | 2022-02-10 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Method and Device for Imaging Diffracted Waves Based on Azimuth-dip Angle Gathers, and Storage Medium |
| CN112379438A (zh) * | 2020-11-02 | 2021-02-19 | 中国石油天然气集团有限公司 | 优化入射角部分叠加地震数据ava特征的方法及装置 |
| CN114924314A (zh) * | 2022-06-24 | 2022-08-19 | 中国石油大学(北京) | 地震剪切横波反演模型构建方法、孔隙度预测方法及装置 |
| CN115267912A (zh) * | 2022-08-02 | 2022-11-01 | 西南石油大学 | 基于压缩感知的avo特征保持去噪方法及系统 |
Non-Patent Citations (4)
| Title |
|---|
| VASCONCELOS IVAN ET AL: "Nonlinear extended images via image-domain interferometry", GEOPHYSICS, vol. 75, no. 6, 1 November 2010 (2010-11-01), pages 105, XP001561998, DOI: 10.1190/1.3494083 * |
| 刘仕友 ET AL: "基于密度的含噪声角道集波形聚类地震相分析", 石油物探, vol. 58, no. 5, 25 September 2019 (2019-09-25), pages 773 - 782 * |
| 吴天麒: "逆时偏移角道集优化处理方法", 中国优秀硕士学位论文全文数据库, no. 2021, 15 March 2021 (2021-03-15), pages 011 - 189 * |
| 徐伟 等: "全方位局部角度域镜像加权叠加技术与应用", 2020油气田勘探与开发国际会议论文集, 14 September 2020 (2020-09-14), pages 410 - 413 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang | Frequencies of the Ricker wavelet | |
| US11467303B2 (en) | Identifying geologic features in a subterranean formation using a post-stack seismic diffraction imaging condition | |
| CN109115987B (zh) | 一种基于岩石物理模型的流体因子的评价方法及装置 | |
| Di et al. | Improved estimates of seismic curvature and flexure based on 3D surface rotation in the presence of structure dip | |
| CN110879412A (zh) | 地下横波速度反演方法、装置、计算设备及存储介质 | |
| CN115270076A (zh) | 一种品质因子确定方法、装置、电子设备及存储介质 | |
| CN103628866A (zh) | 获取地层的地层压力系数的方法 | |
| CN116699684A (zh) | 面向储层预测的地震数据拼接方法、装置、设备及介质 | |
| CN108387934B (zh) | 一种裂缝储层预测方法、装置、电子设备及存储介质 | |
| CN112684498A (zh) | 一种基于宽方位地震数据的储层裂缝预测方法及系统 | |
| CN112083490B (zh) | 地震数据噪声衰减方法及装置 | |
| CN112649848B (zh) | 利用波动方程求解地震波阻抗的方法和装置 | |
| CN116736374B (zh) | 多方位地震数据的融合方法、装置、设备和介质 | |
| CN117991355A (zh) | 砂岩厚度及孔隙度预测方法、装置、电子设备及存储介质 | |
| CN118131318A (zh) | 一种无井区低频模型的建立方法、装置、设备及介质 | |
| CN110632661B (zh) | 一种油气储层预测方法、装置、设备及存储介质 | |
| Yang et al. | Fractal Slope‐Based Seismic Wave Detection Method | |
| RU2783367C1 (ru) | Способ определения положения разломов на основе сейсмических данных | |
| CN114002741B (zh) | 叠前深度偏移方法及装置、计算机可读存储介质 | |
| CN116025420A (zh) | 地层状态变化信息的确定方法、装置、设备、介质及系统 | |
| CN117908109A (zh) | 一种盐岩地层速度的确定方法、装置、设备及介质 | |
| CN115561815A (zh) | 基于扩展弹性阻抗的地震储层预测方法及装置 | |
| CN110888163B (zh) | 微断裂发育位置的确定方法和装置 | |
| CN115267908A (zh) | 余弦相似度地震属性计算方法、装置、电子设备及介质 | |
| CN116201539A (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 |