CN113466936B - 一种断层阴影区crp道集获取方法、装置、设备及存储介质 - Google Patents
一种断层阴影区crp道集获取方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN113466936B CN113466936B CN202110745589.8A CN202110745589A CN113466936B CN 113466936 B CN113466936 B CN 113466936B CN 202110745589 A CN202110745589 A CN 202110745589A CN 113466936 B CN113466936 B CN 113466936B
- Authority
- CN
- China
- Prior art keywords
- signal
- data
- fault
- result
- shadow region
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000005012 migration Effects 0.000 claims abstract description 57
- 238000013508 migration Methods 0.000 claims abstract description 57
- 238000000605 extraction Methods 0.000 claims abstract description 45
- 238000007781 pre-processing Methods 0.000 claims abstract description 13
- 238000012545 processing Methods 0.000 claims description 29
- 230000009466 transformation Effects 0.000 claims description 14
- 238000012216 screening Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 abstract description 5
- 230000003313 weakening effect Effects 0.000 abstract description 5
- 238000003384 imaging method Methods 0.000 description 17
- 238000004458 analytical method Methods 0.000 description 10
- 230000001629 suppression Effects 0.000 description 9
- 230000003287 optical effect Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000005286 illumination Methods 0.000 description 4
- 230000001066 destructive effect Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002194 synthesizing 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/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/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- 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
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)
- Image Processing (AREA)
Abstract
本发明公开了一种断层阴影区CRP道集获取方法、装置、设备及存储介质。该方法包括：对原始数据进行前置处理，得到第一数据；对所述原始数据进行优势信号提取，得到第一目标优势信号；根据所述第一目标优势信号进行叠前偏移速度建模，得到叠前偏移速度模型；将第一数据输入所述叠前偏移速度模型得到CRP道集，通过本发明的技术方案，能够通过对断层阴影区开展优势信号提取，获得信噪比真实的地震信号，并以此为标准对断层阴影区外的范围开展自适应信号约束提取。对信噪比横向一致性较为稳定的地震数据开展精细叠前偏移速度模型建立，获得一致性得到改善的数据体，最终实现断层阴影区边界的弱化，获得对整体构造的准确认识。
Description
技术领域
本发明实施例涉及勘探技术领域，尤其涉及一种断层阴影区CRP道集获取方法、装置、设备及存储介质。
背景技术
随着精细高效勘探的要求，构造的精细落实及构造圈闭的准确描述是对物探技术的基本要求和迫切需要。地震激发信号在地下介质中传播时，遇到断层会产生地震响应的变化。因此，地震剖面上常常会见到断层上盘地层为连续反射，而断层下盘存在一个“三角区”，其同相轴成像时间出现“上拉”或者“下拉”的现象或者同相轴出现错断，这种不确定的“上拉”或“下拉”的现象被认为是成像畸变。地质上对此类断层下盘从上至下的一个三角形区域用“断层阴影区”来定义。断层阴影区地震成像表现为信噪比低、同相轴产状畸变以及能量明显减弱的现象，被认为成像不可靠。阴影区的存在极大影响了对构造形态的精确认识、圈闭幅度的准确判断以及对构造高点的确定。
长期以来，针对断层阴影区的构造不能够精确落实的问题，常常被认为是断层下方的速度求取不够准确。因此，开展针对复杂断层高精度速度建模技术(如全波形反演)和高精度成像技术(如逆时偏移)被认为是解决断层阴影区成像的主要处理技术解决方案。对于利用成像技术而依然无法完全消除的断层阴影问题，则是通过开展基于地质模式的构造解释及成图方法的调整，来获得对断层阴影区构造的认识。以上两类解决方法，其实质都缺少地震数据信号的强大支持，是一种地质认识指导下的技术解决方案，存在较大的主观性和多解性。
发明内容
本发明实施例提供一种断层阴影区CRP道集获取方法、装置、设备及存储介质，利用地震采集获得的原始数据，对断层阴影区内及其以外的反射信号进行不同处理阶段的多维度信噪比分析，掌握地质体接收地震信号的照明差异。通过对断层阴影区开展优势信号提取，获得信噪比真实的地震信号，并以此为标准对断层阴影区外的范围开展自适应信号约束提取。对信噪比横向一致性较为稳定的地震数据开展精细叠前偏移速度模型建立，获得一致性得到改善的数据体，最终实现断层阴影区边界的弱化，获得对整体构造的准确认识。
第一方面，本发明实施例提供了一种断层阴影区CRP道集获取方法，包括：
对原始数据进行前置处理，得到第一数据；
对所述原始数据进行优势信号提取，得到第一目标优势信号；
根据所述第一目标优势信号进行叠前偏移速度建模，得到叠前偏移速度模型；
将第一数据输入所述叠前偏移速度模型得到CRP道集。
第二方面，本发明实施例还提供了一种断层阴影区CRP道集获取装置，该装置包括：
处理模块，用于对原始数据进行前置处理，得到第一数据；
第一提取模块，用于对所述原始数据进行优势信号提取，得到第一目标优势信号；
建模模块，用于根据所述第一目标优势信号进行叠前偏移速度建模，得到叠前偏移速度模型；
获取模块，用于将第一数据输入所述叠前偏移速度模型得到CRP道集。
第三方面，本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如本发明实施例中任一所述的断层阴影区CRP道集获取方法。
第四方面，本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如本发明实施例中任一所述的断层阴影区CRP道集获取方法。
本发明实施例通过对原始数据进行前置处理，得到第一数据；对所述原始数据进行优势信号提取，得到第一目标优势信号；根据所述第一目标优势信号进行叠前偏移速度建模，得到叠前偏移速度模型；将第一数据输入所述叠前偏移速度模型得到CRP道集，利用地震采集获得的原始数据，对断层阴影区内及其以外的反射信号进行不同处理阶段的多维度信噪比分析，掌握地质体接收地震信号的照明差异。通过对断层阴影区开展优势信号提取，获得信噪比真实的地震信号，并以此为标准对断层阴影区外的范围开展自适应信号约束提取。对信噪比横向一致性较为稳定的地震数据开展精细叠前偏移速度模型建立，获得一致性得到改善的数据体，最终实现断层阴影区边界的弱化，获得对整体构造的准确认识。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1是本发明实施例一中的一种断层阴影区CRP道集获取方法的流程图；
图1a是本发明实施例一中改善断层阴影区优势信号提取处理方法的流程图；
图1b是本发明实施例一中的成果对比图；
图2是本发明实施例二中的一种断层阴影区CRP道集获取装置的结构示意图；
图3是本发明实施例三中的一种计算机设备的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。同时，在本发明的描述中，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
实施例一
图1为本发明实施例一提供的一种断层阴影区CRP道集获取方法的流程图，本实施例可适用于CRP道集获取的情况，该方法可以由本发明实施例中的CRP道集获取装置来执行，该装置可采用软件和/或硬件的方式实现，如图1所示，该方法具体包括如下步骤：
S110，对原始数据进行前置处理，得到第一数据。
其中，所述前置处理包括但不限于：干扰波压制、鬼波压制、多次波压制、能量及频率补偿和数据规则化等。
示例性的，对原始数据进行前置处理，得到第一数据，例如可以是，对原始数据依次进行干扰波压制、鬼波压制、多次波压制、能量及频率补偿和数据规则化后，得到第一数据。
S120，对所述原始数据进行优势信号提取，得到第一目标优势信号。
示例性的，对所述原始数据进行优势信号提取，得到第一目标优势信号的方式可以为：从所述原始数据中筛选断层阴影区的叠加剖面对应的第一断层数据；对第一断层数据进行分频扫描，得到第一扫描结果；根据所述第一扫描结果确定第一断层阴影区空间范围；对所述第一断层阴影区空间范围内的原始数据进行小波变换，分解成至少一个第一小波域数据；获取每个第一小波域数据对应的第一沿层信噪比属性；根据所述第一沿层信噪比属性对所述第一断层阴影区空间范围进行优化，得到第一目标阴影区空间范围；根据所述第一沿层信噪比属性确定第一振幅比例系数；根据所述第一振幅比例系数从所述第一目标阴影区空间范围中提取第一目标优势信号。
在一个示例的例子中，对包括断层阴影区的叠加剖面进行分频扫描，对分频扫描结果进行沿层信噪比统计、沿层振幅统计、垂向振幅统计，综合分频的沿层信噪比及振幅信息、垂向振幅信息，初步划定断层阴影区空间范围及优势信号范围。其中低频部分横向信噪比及振幅一致性好，为优势频率；高频端信号在断层阴影区内振幅弱、信噪比低，与阴影区外有明显区别，对阴影区成像效果起到破坏作用。对不同处理阶段的原始单炮开展针对优势频带的处理，结合单炮和叠加剖面联合并迭代开展信噪比分析和振幅统计，落实优势频带的范围。基于叠后地震数据开展小波变换，分解成若干小波域数据，并统计相对应的沿层信噪比属性，在断层阴影区外，振幅比例系数A设定为1,在断层阴影影响区域内，给定与信噪比反相关的振幅比例系数A，对不同信号进行比例后再重构获得断层阴影区优势信号提取后的数据。分阶段开展信噪比统计、振幅统计与优势信号提取，针对每步处理得到的阶段成果数据进行优势信号提取，得到第一目标优势信号。
S130，根据所述第一目标优势信号进行叠前偏移速度建模，得到叠前偏移速度模型。
示例性的，根据所述第一目标优势信号进行叠前偏移速度建模，得到叠前偏移速度模型，可以获得成像精度更高的偏移速度模型。
S140，将第一数据输入所述叠前偏移速度模型得到CRP道集。
可选的，在将第一数据输入所述叠前偏移速度模型得到CRP道集之后，还包括：
对所述CRP道集进行叠加，得到第一叠加结果；
对所述CRP道集进行优势信号提取，得到第二目标优势信号；
根据所述第二目标优势信号进行叠加，得到第二叠加结果；
根据所述第一叠加结果和所述第二叠加结果生成断层阴影区构造解释。
示例性的，对所述CRP道集进行优势信号提取，得到第二目标优势信号的方式可以为：从所述CRP道集中筛选断层阴影区的叠加剖面对应的第二断层数据；对第二断层数据进行分频扫描，得到第二扫描结果；根据所述第二扫描结果确定第二断层阴影区空间范围；对所述第二断层阴影区空间范围内的原始数据进行小波变换，分解成至少一个第二小波域数据；获取每个第二小波域数据对应的第二沿层信噪比属性；根据所述第二沿层信噪比属性对所述第二断层阴影区空间范围进行优化，得到第二目标阴影区空间范围；根据所述第二沿层信噪比属性确定第二振幅比例系数；根据所述第二振幅比例系数从所述第二目标阴影区空间范围中提取第二目标优势信号。
可选的，对所述原始数据进行优势信号提取，得到第一目标优势信号，包括：
从所述原始数据中筛选断层阴影区的叠加剖面对应的第一断层数据；
对第一断层数据进行分频扫描，得到第一扫描结果；
根据所述第一扫描结果确定第一断层阴影区空间范围；
对所述第一断层阴影区空间范围内的原始数据进行小波变换，分解成至少一个第一小波域数据；
获取每个第一小波域数据对应的第一沿层信噪比属性；
根据所述第一沿层信噪比属性对所述第一断层阴影区空间范围进行优化，得到第一目标阴影区空间范围；
根据所述第一沿层信噪比属性确定第一振幅比例系数；
根据所述第一振幅比例系数从所述第一目标阴影区空间范围中提取第一目标优势信号。
示例性的，从所述原始数据中筛选断层阴影区的叠加剖面对应的第一断层数据；对第一断层数据进行分频扫描，得到第一扫描结果；根据所述第一扫描结果确定第一断层阴影区空间范围；对所述第一断层阴影区空间范围内的原始数据进行小波变换，分解成至少一个第一小波域数据；获取每个第一小波域数据对应的第一沿层信噪比属性；根据所述第一沿层信噪比属性对所述第一断层阴影区空间范围进行优化，得到第一目标阴影区空间范围；根据所述第一沿层信噪比属性确定第一振幅比例系数；根据所述第一振幅比例系数从所述第一目标阴影区空间范围中提取第一目标优势信号。例如可以是，对包括断层阴影区的叠加剖面进行分频扫描，对分频扫描结果进行沿层信噪比统计、沿层振幅统计、垂向振幅统计，综合分频的沿层信噪比及振幅信息、垂向振幅信息，初步划定断层阴影区空间范围及优势信号范围。其中低频(一般18Hz以内)部分横向信噪比及振幅一致性好，为优势频率；高频端信号在断层阴影区内振幅弱、信噪比低，与阴影区外有明显区别，对阴影区成像效果起到破坏作用。对不同处理阶段的原始单炮开展针对优势频带的处理，结合单炮和叠加剖面联合并迭代开展信噪比分析和振幅统计，落实优势频带的范围。基于叠后地震数据开展小波变换，分解成若干小波域数据，并统计相对应的沿层信噪比属性，在断层阴影区外，振幅比例系数A设定为1,在断层阴影影响区域内，给定与信噪比反相关的振幅比例系数A，对不同信号进行比例后再重构获得断层阴影区优势信号提取后的数据。分阶段开展信噪比统计、振幅统计与优势信号提取，针对每步处理得到的阶段成果数据进行优势信号提取，得到第一目标优势信号。
可选的，根据所述第一扫描结果确定第一断层阴影区空间范围，包括：
对所述第一扫描结果进行沿层信噪比统计，得到第一统计结果；
对所述第一扫描结果进行沿层振幅统计，得到第二统计结果；
对所述第一扫描结果进行垂向振幅统计，得到第三统计结果；
根据所述第一统计结果、第二统计结果和第三统计结果确定第一断层阴影区空间范围。
可选的，对所述CRP道集进行优势信号提取，得到第二目标优势信号，包括：
从所述CRP道集中筛选断层阴影区的叠加剖面对应的第二断层数据；
对第二断层数据进行分频扫描，得到第二扫描结果；
根据所述第二扫描结果确定第二断层阴影区空间范围；
对所述第二断层阴影区空间范围内的原始数据进行小波变换，分解成至少一个第二小波域数据；
获取每个第二小波域数据对应的第二沿层信噪比属性；
根据所述第二沿层信噪比属性对所述第二断层阴影区空间范围进行优化，得到第二目标阴影区空间范围；
根据所述第二沿层信噪比属性确定第二振幅比例系数；
根据所述第二振幅比例系数从所述第二目标阴影区空间范围中提取第二目标优势信号。
可选的，根据所述第二扫描结果确定第二断层阴影区空间范围，包括：
对所述第二扫描结果进行沿层信噪比统计，得到第四统计结果；
对所述第二扫描结果进行沿层振幅统计，得到第五统计结果；
对所述第二扫描结果进行垂向振幅统计，得到第六统计结果；
根据所述第四统计结果、第五统计结果和第六统计结果确定第二断层阴影区空间范围。
在一个示例的例子中，本发明通过对断层发育区开展优势信号提取，首先获得断层阴影区信噪比的提高。参考断层阴影区的优势信号提取标准完成阴影区外优势信号的自适应提取。优势信号以地震信号的频率和能量作为主要参数，通过偏移前CMP道集的优势信号提取及偏移后CRP道集优势信号提取，得到最终改善断层阴影区成像的地震数据成果。如图1a所示，本发明实施例提供的改善断层阴影区优势信号提取处理方法，包括以下步骤：
1、针对断层发育区开展分阶段多维度信噪比统计分析，主要处理步骤包括但不限于：干扰波压制、鬼波压制、多次波压制、能量及频率补偿和数据规则化等。对断层阴影区域内外的资料信噪比进行差异性分析，总结断层阴影区内外地震信号的响应差异；
2、在每个处理步骤开展信噪比分析及振幅统计，对断层阴影区存在的范围进行划定。利用层位解释成果，对叠加剖面开展沿层信噪比分析和沿层振幅提取，通过对多个处理阶段的叠加剖面开展分析，对断层阴影区的影响范围进行划定。选择断层阴影区内外不同位置所对应的原始单炮进行标志层位的信噪比分析和振幅统计。
3、对包括断层阴影区的叠加剖面进行分频扫描，对分频扫描结果进行沿层信噪比统计、沿层振幅统计、垂向振幅统计，综合分频的沿层信噪比及振幅信息、垂向振幅信息，初步划定断层阴影区空间范围及优势信号范围。其中低频(一般18Hz以内)部分横向信噪比及振幅一致性好，为优势频率；高频端信号在断层阴影区内振幅弱、信噪比低，与阴影区外有明显区别，对阴影区成像效果起到破坏作用。对不同处理阶段的原始单炮开展针对优势频带的处理，结合单炮和叠加剖面联合并迭代开展信噪比分析和振幅统计，落实优势频带的范围。基于叠后地震数据开展小波变换，分解成若干小波域数据，并统计相对应的沿层信噪比属性，在断层阴影区外，振幅比例系数A设定为1,在断层阴影影响区域内，给定与信噪比反相关的振幅比例系数A，对不同信号进行比例后再重构获得断层阴影区优势信号提取后的数据。分阶段开展信噪比统计、振幅统计与优势信号提取，针对每步处理得到的阶段成果数据进行优势信号提取，以此用于验证优势信号提取技术对于改善断层阴影区成像的适用效果。信噪比及能量的横向一致性分析结果可以作为优势信号提取的主要质控手段；
4、利用消除断层阴影区的CMP道集数据开展叠前偏移速度建模，可以获得成像精度更高的偏移速度模型，并利用此速度模型开展偏移成像。
5、对最终速度模型偏移成像获得的CRP道集进行叠加，对叠加剖面再次开展信噪比统计和优势信号提取，最终获得能够有效改善断层阴影区成像的地震数据。
6、对CRP道集开展优势信号提取，可以提高断层阴影区的道集信噪比，缩小CRP道集间的横向信噪比差异。优势信号提取标准可以参考井旁道集的AVO特征。优势信号提取的第一要素为精细落实构造，第二要素为振幅保真。
与现有技术相比，本发明实施例在处理过程中分阶段使用了优势信号提取技术，对断层阴影区内及其以外进行自适应优势信号提取，解决了常规处理存在的横向信噪比差异、断层阴影区界限特征明显的问题，具有以下优点：1、本发明分阶段开展信噪比分析与统计，对断层阴影内外信号认识更具有真实性和客观性；2、本发明在偏移前CMP道集(CommonMidpoint，共中心点道集)开展优势信号提取，用此数据进行叠前偏移速度建模，提高了速度模型精度；3、本发明应用于偏移后CRP道集上，获得信噪比特征一致的地震数据成果，而叠加数据成果有益于开展精细构造解释，道集数据可以辅助开展叠前储层预测研究；4、本发明还适用于低信噪比地区的地震资料处理及精细构造解释，可广泛应用于勘探开发的地震资料处理解释工作中。
应用本发明实施例提供的优势信号提取技术，获得信噪比一致性较好的CMP道集、CRP道集及偏移数据。如图1b所示，从地震数据源头实现对断层阴影区有效信号的识别和应用，提高了断层阴影区资料的信噪比和速度谱质量。通过采取优势信号提取缩小了断层阴影区与周边地层成像的品质差异，淡化了模糊区的范围影响，对于开展整体构造解释提供了较好的数据基础。
上述实施例仅用于说明本发明，其中各步骤等都是可以有所变化的，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。
本实施例的技术方案，通过对原始数据进行前置处理，得到第一数据；对所述原始数据进行优势信号提取，得到第一目标优势信号；根据所述第一目标优势信号进行叠前偏移速度建模，得到叠前偏移速度模型；将第一数据输入所述叠前偏移速度模型得到CRP道集，利用地震采集获得的原始数据，对断层阴影区内及其以外的反射信号进行不同处理阶段的多维度信噪比分析，掌握地质体接收地震信号的照明差异。通过对断层阴影区开展优势信号提取，获得信噪比真实的地震信号，并以此为标准对断层阴影区外的范围开展自适应信号约束提取。对信噪比横向一致性较为稳定的地震数据开展精细叠前偏移速度模型建立，获得一致性得到改善的数据体，最终实现断层阴影区边界的弱化，获得对整体构造的准确认识。
实施例二
图2为本发明实施例二提供的一种断层阴影区CRP道集获取装置的结构示意图。本实施例可适用于断层阴影区CRP道集获取的情况，该装置可采用软件和/或硬件的方式实现，该装置可集成在任何提供断层阴影区CRP道集获取功能的设备中，如图2所示，所述断层阴影区CRP道集获取装置具体包括：处理模块210、第一提取模块220、建模模块230和获取模块240。
其中，处理模块，用于对原始数据进行前置处理，得到第一数据；
第一提取模块，用于对所述原始数据进行优势信号提取，得到第一目标优势信号；
建模模块，用于根据所述第一目标优势信号进行叠前偏移速度建模，得到叠前偏移速度模型；
获取模块，用于将第一数据输入所述叠前偏移速度模型得到CRP道集。
可选的，还包括：
第一叠加模块，用于在将第一数据输入所述叠前偏移速度模型得到CRP道集之后，对所述CRP道集进行叠加，得到第一叠加结果；
第二提取模块，用于对所述CRP道集进行优势信号提取，得到第二目标优势信号；
第二叠加模块，用于根据所述第二目标优势信号进行叠加，得到第二叠加结果；
生成模块，用于根据所述第一叠加结果和所述第二叠加结果生成断层阴影区构造解释。
上述产品可执行本发明任意实施例所提供的方法，具备执行方法相应的功能模块和有益效果。
本实施例的技术方案，通过对原始数据进行前置处理，得到第一数据；对所述原始数据进行优势信号提取，得到第一目标优势信号；根据所述第一目标优势信号进行叠前偏移速度建模，得到叠前偏移速度模型；将第一数据输入所述叠前偏移速度模型得到CRP道集，利用地震采集获得的原始数据，对断层阴影区内及其以外的反射信号进行不同处理阶段的多维度信噪比分析，掌握地质体接收地震信号的照明差异。通过对断层阴影区开展优势信号提取，获得信噪比真实的地震信号，并以此为标准对断层阴影区外的范围开展自适应信号约束提取。对信噪比横向一致性较为稳定的地震数据开展精细叠前偏移速度模型建立，获得一致性得到改善的数据体，最终实现断层阴影区边界的弱化，获得对整体构造的准确认识。
实施例三
图3为本发明实施例三中的一种计算机设备的结构示意图。图3示出了适于用来实现本发明实施方式的示例性计算机设备12的框图。图3显示的计算机设备12仅仅是一个示例，不应对本发明实施例的功能和使用范围带来任何限制。
如图3所示，计算机设备12以通用计算设备的形式表现。计算机设备12的组件可以包括但不限于：一个或者多个处理器或者处理单元16，系统存储器28，连接不同系统组件(包括系统存储器28和处理单元16)的总线18。
总线18表示几类总线结构中的一种或多种，包括存储器总线或者存储器控制器，外围总线，图形加速端口，处理器或者使用多种总线结构中的任意总线结构的局域总线。举例来说，这些体系结构包括但不限于工业标准体系结构(ISA)总线，微通道体系结构(MAC)总线，增强型ISA总线、视频电子标准协会(VESA)局域总线以及外围组件互连(PCI)总线。
计算机设备12典型地包括多种计算机系统可读介质。这些介质可以是任何能够被计算机设备12访问的可用介质，包括易失性和非易失性介质，可移动的和不可移动的介质。
系统存储器28可以包括易失性存储器形式的计算机系统可读介质，例如随机存取存储器(RAM)30和/或高速缓存存储器32。计算机设备12可以进一步包括其它可移动/不可移动的、易失性/非易失性计算机系统存储介质。仅作为举例，存储系统34可以用于读写不可移动的、非易失性磁介质(图3未显示，通常称为“硬盘驱动器”)。尽管图3中未示出，可以提供用于对可移动非易失性磁盘(例如“软盘”)读写的磁盘驱动器，以及对可移动非易失性光盘(例如CD-ROM,DVD-ROM或者其它光介质)读写的光盘驱动器。在这些情况下，每个驱动器可以通过一个或者多个数据介质接口与总线18相连。存储器28可以包括至少一个程序产品，该程序产品具有一组(例如至少一个)程序模块，这些程序模块被配置以执行本发明各实施例的功能。
具有一组(至少一个)程序模块42的程序/实用工具40，可以存储在例如存储器28中，这样的程序模块42包括——但不限于——操作系统、一个或者多个应用程序、其它程序模块以及程序数据，这些示例中的每一个或某种组合中可能包括网络环境的实现。程序模块42通常执行本发明所描述的实施例中的功能和/或方法。
计算机设备12也可以与一个或多个外部设备14(例如键盘、指向设备、显示器24等)通信，还可与一个或者多个使得用户能与该计算机设备12交互的设备通信，和/或与使得该计算机设备12能与一个或多个其它计算设备进行通信的任何设备(例如网卡，调制解调器等等)通信。这种通信可以通过输入/输出(I/O)接口22进行。另外，本实施例中的计算机设备12，显示器24不是作为独立个体存在，而是嵌入镜面中，在显示器24的显示面不予显示时，显示器24的显示面与镜面从视觉上融为一体。并且，计算机设备12还可以通过网络适配器20与一个或者多个网络(例如局域网(LAN)，广域网(WAN)和/或公共网络，例如因特网)通信。如图所示，网络适配器20通过总线18与计算机设备12的其它模块通信。应当明白，尽管图中未示出，可以结合计算机设备12使用其它硬件和/或软件模块，包括但不限于：微代码、设备驱动器、冗余处理单元、外部磁盘驱动阵列、RAID系统、磁带驱动器以及数据备份存储系统等。
处理单元16通过运行存储在系统存储器28中的程序，从而执行各种功能应用以及数据处理，例如实现本发明实施例所提供的断层阴影区CRP道集获取方法：
对原始数据进行前置处理，得到第一数据；
对所述原始数据进行优势信号提取，得到第一目标优势信号；
根据所述第一目标优势信号进行叠前偏移速度建模，得到叠前偏移速度模型；
将第一数据输入所述叠前偏移速度模型得到CRP道集。
实施例四
本发明实施例四提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如本申请所有发明实施例提供的断层阴影区CRP道集获取方法：
对原始数据进行前置处理，得到第一数据；
对所述原始数据进行优势信号提取，得到第一目标优势信号；
根据所述第一目标优势信号进行叠前偏移速度建模，得到叠前偏移速度模型；
将第一数据输入所述叠前偏移速度模型得到CRP道集。
可以采用一个或多个计算机可读的介质的任意组合。计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。计算机可读存储介质例如可以是但不限于电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括但不限于无线、电线、光缆、RF等等，或者上述的任意合适的组合。
可以以一种或多种程序设计语言或其组合来编写用于执行本发明操作的计算机程序代码，所述程序设计语言包括面向对象的程序设计语言，诸如Java、Smalltalk、C++，还包括常规的过程式程序设计语言，诸如“C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。
Claims (8)
1.一种断层阴影区CRP道集获取方法，其特征在于，包括：
对原始数据进行前置处理，得到第一数据；
对所述原始数据进行优势信号提取，得到第一目标优势信号；
根据所述第一目标优势信号进行叠前偏移速度建模，得到叠前偏移速度模型；
将第一数据输入所述叠前偏移速度模型得到CRP道集；
对所述原始数据进行优势信号提取，得到第一目标优势信号，包括：
从所述原始数据中筛选断层阴影区的叠加剖面对应的第一断层数据；
对第一断层数据进行分频扫描，得到第一扫描结果；
根据所述第一扫描结果确定第一断层阴影区空间范围；
对所述第一断层阴影区空间范围内的原始数据进行小波变换，分解成至少一个第一小波域数据；
获取每个第一小波域数据对应的第一沿层信噪比属性；
根据所述第一沿层信噪比属性对所述第一断层阴影区空间范围进行优化，得到第一目标阴影区空间范围；
根据所述第一沿层信噪比属性确定第一振幅比例系数；
根据所述第一振幅比例系数从所述第一目标阴影区空间范围中提取第一目标优势信号；
根据所述第一扫描结果确定第一断层阴影区空间范围，包括：
对所述第一扫描结果进行沿层信噪比统计，得到第一统计结果；
对所述第一扫描结果进行沿层振幅统计，得到第二统计结果；
对所述第一扫描结果进行垂向振幅统计，得到第三统计结果；
根据所述第一统计结果、第二统计结果和第三统计结果确定第一断层阴影区空间范围。
2.根据权利要求1所述的方法，其特征在于，在将第一数据输入所述叠前偏移速度模型得到CRP道集之后，还包括：
对所述CRP道集进行叠加，得到第一叠加结果；
对所述CRP道集进行优势信号提取，得到第二目标优势信号；
根据所述第二目标优势信号进行叠加，得到第二叠加结果；
根据所述第一叠加结果和所述第二叠加结果生成断层阴影区构造解释。
3.根据权利要求2所述的方法，其特征在于，对所述CRP道集进行优势信号提取，得到第二目标优势信号，包括：
从所述CRP道集中筛选断层阴影区的叠加剖面对应的第二断层数据；
对第二断层数据进行分频扫描，得到第二扫描结果；
根据所述第二扫描结果确定第二断层阴影区空间范围；
对所述第二断层阴影区空间范围内的原始数据进行小波变换，分解成至少一个第二小波域数据；
获取每个第二小波域数据对应的第二沿层信噪比属性；
根据所述第二沿层信噪比属性对所述第二断层阴影区空间范围进行优化，得到第二目标阴影区空间范围；
根据所述第二沿层信噪比属性确定第二振幅比例系数；
根据所述第二振幅比例系数从所述第二目标阴影区空间范围中提取第二目标优势信号。
4.根据权利要求3所述的方法，其特征在于，根据所述第二扫描结果确定第二断层阴影区空间范围，包括：
对所述第二扫描结果进行沿层信噪比统计，得到第四统计结果；
对所述第二扫描结果进行沿层振幅统计，得到第五统计结果；
对所述第二扫描结果进行垂向振幅统计，得到第六统计结果；
根据所述第四统计结果、第五统计结果和第六统计结果确定第二断层阴影区空间范围。
5.一种断层阴影区CRP道集获取装置，其特征在于，包括：
处理模块，用于对原始数据进行前置处理，得到第一数据；
第一提取模块，用于对所述原始数据进行优势信号提取，得到第一目标优势信号；
建模模块，用于根据所述第一目标优势信号进行叠前偏移速度建模，得到叠前偏移速度模型；
获取模块，用于将第一数据输入所述叠前偏移速度模型得到CRP道集；
所述第一提取模块具体用于：
从所述原始数据中筛选断层阴影区的叠加剖面对应的第一断层数据；
对第一断层数据进行分频扫描，得到第一扫描结果；
根据所述第一扫描结果确定第一断层阴影区空间范围；
对所述第一断层阴影区空间范围内的原始数据进行小波变换，分解成至少一个第一小波域数据；
获取每个第一小波域数据对应的第一沿层信噪比属性；
根据所述第一沿层信噪比属性对所述第一断层阴影区空间范围进行优化，得到第一目标阴影区空间范围；
根据所述第一沿层信噪比属性确定第一振幅比例系数；
根据所述第一振幅比例系数从所述第一目标阴影区空间范围中提取第一目标优势信号；
所述第一提取模块还用于：
对所述第一扫描结果进行沿层信噪比统计，得到第一统计结果；
对所述第一扫描结果进行沿层振幅统计，得到第二统计结果；
对所述第一扫描结果进行垂向振幅统计，得到第三统计结果；
根据所述第一统计结果、第二统计结果和第三统计结果确定第一断层阴影区空间范围。
6.根据权利要求5所述的装置，其特征在于，还包括：
第一叠加模块，用于在将第一数据输入所述叠前偏移速度模型得到CRP道集之后，对所述CRP道集进行叠加，得到第一叠加结果；
第二提取模块，用于对所述CRP道集进行优势信号提取，得到第二目标优势信号；
第二叠加模块，用于根据所述第二目标优势信号进行叠加，得到第二叠加结果；
生成模块，用于根据所述第一叠加结果和所述第二叠加结果生成断层阴影区构造解释。
7.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1-4中任一所述的方法。
8.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现如权利要求1-4中任一所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110745589.8A CN113466936B (zh) | 2021-07-01 | 2021-07-01 | 一种断层阴影区crp道集获取方法、装置、设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110745589.8A CN113466936B (zh) | 2021-07-01 | 2021-07-01 | 一种断层阴影区crp道集获取方法、装置、设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113466936A CN113466936A (zh) | 2021-10-01 |
| CN113466936B true CN113466936B (zh) | 2023-07-18 |
Family
ID=77877277
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110745589.8A Active CN113466936B (zh) | 2021-07-01 | 2021-07-01 | 一种断层阴影区crp道集获取方法、装置、设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113466936B (zh) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5808964A (en) * | 1995-06-06 | 1998-09-15 | Institut Francais Du Petrole | Processing method for obtaining zero-offset seismic data by depth domain stacking |
| CN105866832A (zh) * | 2015-01-20 | 2016-08-17 | 中国石油天然气股份有限公司 | 深层致密砂岩气藏的低级序断层的识别方法和装置 |
| CN106249292A (zh) * | 2015-06-09 | 2016-12-21 | 中国石油化工股份有限公司 | 一种共反射点道集优化处理方法及系统 |
| CN107526101A (zh) * | 2017-07-07 | 2017-12-29 | 中国地质调查局油气资源调查中心 | 一种获取地震反射波的采集和处理方法 |
| CA2974134A1 (en) * | 2016-09-09 | 2018-03-09 | China Petroleum & Chemical Corporation | Method and system for seismic data processing |
| WO2018203265A1 (en) * | 2017-05-04 | 2018-11-08 | Chevron U.S.A. Inc. | System and method for analyzing geologic features using seismic data |
| CN109709605A (zh) * | 2019-02-28 | 2019-05-03 | 中国石油化工股份有限公司 | 火成岩复杂区碳酸盐岩断溶体地震数据成像方法及设备 |
| CN111399056A (zh) * | 2020-04-29 | 2020-07-10 | 西南石油大学 | 一种基于划分方位角滤波预测裂缝强度的方法 |
| CN111582114A (zh) * | 2020-04-29 | 2020-08-25 | 中国海洋石油集团有限公司 | 一种地震断层识别方法、装置、设备和存储介质 |
| CN112305612A (zh) * | 2019-07-23 | 2021-02-02 | 中国海洋石油集团有限公司 | 高分辨率复谱分解时频空间域振幅随偏移距变化校正方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2006104757A (ru) * | 2003-08-19 | 2006-07-27 | Тетрасейс Инк. (Ca) | Способ интерпретации сейсмических данных, полученных с использованием дуплексных волн |
-
2021
- 2021-07-01 CN CN202110745589.8A patent/CN113466936B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5808964A (en) * | 1995-06-06 | 1998-09-15 | Institut Francais Du Petrole | Processing method for obtaining zero-offset seismic data by depth domain stacking |
| CN105866832A (zh) * | 2015-01-20 | 2016-08-17 | 中国石油天然气股份有限公司 | 深层致密砂岩气藏的低级序断层的识别方法和装置 |
| CN106249292A (zh) * | 2015-06-09 | 2016-12-21 | 中国石油化工股份有限公司 | 一种共反射点道集优化处理方法及系统 |
| CA2974134A1 (en) * | 2016-09-09 | 2018-03-09 | China Petroleum & Chemical Corporation | Method and system for seismic data processing |
| WO2018203265A1 (en) * | 2017-05-04 | 2018-11-08 | Chevron U.S.A. Inc. | System and method for analyzing geologic features using seismic data |
| CN107526101A (zh) * | 2017-07-07 | 2017-12-29 | 中国地质调查局油气资源调查中心 | 一种获取地震反射波的采集和处理方法 |
| CN109709605A (zh) * | 2019-02-28 | 2019-05-03 | 中国石油化工股份有限公司 | 火成岩复杂区碳酸盐岩断溶体地震数据成像方法及设备 |
| CN112305612A (zh) * | 2019-07-23 | 2021-02-02 | 中国海洋石油集团有限公司 | 高分辨率复谱分解时频空间域振幅随偏移距变化校正方法 |
| CN111399056A (zh) * | 2020-04-29 | 2020-07-10 | 西南石油大学 | 一种基于划分方位角滤波预测裂缝强度的方法 |
| CN111582114A (zh) * | 2020-04-29 | 2020-08-25 | 中国海洋石油集团有限公司 | 一种地震断层识别方法、装置、设备和存储介质 |
Non-Patent Citations (7)
| Title |
|---|
| Kirchhoff modeling, inversion for reflectivity, and subsurface illumination;Bertrand Duquet;《Geophysics》;第65卷(第04期);第1195-1209页 * |
| 地震构造解释中断层阴影区假断层现象分析――以松辽盆地北部敖古拉断层为例;李操;《新疆石油地质》;第41卷(第2期);第223-247页 * |
| 基于CRP道集优化处理的叠前AVA同步反演技术的应用—以KL9构造区为例;邓吉锋 等;《石油物探》;第58卷(第03期);第461-470页 * |
| 基于卷积神经网络的断层预测方法;刘江 等;《CT理论与应用研究》;第29卷(第05期);第522-533页 * |
| 基于叠前优势道集的断裂识别技术研究;胡天乐 等;《中国化工贸易》(第32期);第322页 * |
| 基于叠前道集小波域Q值求取方法研究;魏文 等;《石油物探》;第50卷(第04期);第355-360页 * |
| 威廉·艾思克罗夫特.《地震反射波勘探》.石油大学出版社,2017,第140-143页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113466936A (zh) | 2021-10-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Li et al. | A method for low-frequency noise suppression based on mathematical morphology in microseismic monitoring | |
| US20160320509A1 (en) | Suppressing near-surface scattered surface waves | |
| AU2015383134B2 (en) | Multistage full wavefield inversion process that generates a multiple free data set | |
| US11467303B2 (en) | Identifying geologic features in a subterranean formation using a post-stack seismic diffraction imaging condition | |
| US10534101B2 (en) | Seismic adaptive focusing | |
| US9753166B2 (en) | P-wave and S-wave separation of seismic data in the presence of statics and irregular geometry | |
| US11561312B2 (en) | Mapping near-surface heterogeneities in a subterranean formation | |
| NO20240301A1 (en) | Random noise attenuation for seismic data | |
| CN113466936B (zh) | 一种断层阴影区crp道集获取方法、装置、设备及存储介质 | |
| CN111025393B (zh) | 针对含薄煤层地层的储层预测方法、装置、设备及介质 | |
| Verma et al. | Data conditioning of legacy seismic using migration-driven 5D interpolation | |
| CN112200069A (zh) | 时频域谱减法和经验模态分解联合的隧道滤波方法及系统 | |
| Wang et al. | Poststack seismic inversion using a patch-based Gaussian mixture model | |
| Zand et al. | Least-squares reverse time migration with shifted total variation regularization | |
| CN112147700A (zh) | 速度异常区的低频模型构建方法及系统 | |
| US11215725B2 (en) | Seismic processing workflow for orthogonal wide azimuth 3D surveys | |
| US11768303B2 (en) | Automatic data enhancement for full waveform inversion in the midpoint-offset domain | |
| Ha et al. | The value of constrained conjugate-gradient least-squares migration in seismic inversion: Application to a fractured-basement play, Texas Panhandle | |
| US11415719B2 (en) | Method of application of polarization filtering on single component seismic data for interface wave noise attenuation | |
| US20230266492A1 (en) | Enhancement of single sensor seismic data | |
| US11880010B2 (en) | Providing seismic images of the subsurface using enhancement of pre-stack seismic data | |
| US20230266493A1 (en) | Seismic denoising by weighted multiplanar reconstruction | |
| Wang et al. | Lithological Reservoir Exploration Technology | |
| CN117784254A (zh) | 一种基于人工智能的多尺度岩性预测方法、电子设备及介质 | |
| Hu et al. | Coherent noise attenuation for passive seismic data based on iterative two-dimensional model shrinkage |
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 |