CN112987100B - 预测受层序界面控制的弱振幅岩溶层中储层的方法和装置 - Google Patents
预测受层序界面控制的弱振幅岩溶层中储层的方法和装置 Download PDFInfo
- Publication number
- CN112987100B CN112987100B CN201911298585.9A CN201911298585A CN112987100B CN 112987100 B CN112987100 B CN 112987100B CN 201911298585 A CN201911298585 A CN 201911298585A CN 112987100 B CN112987100 B CN 112987100B
- Authority
- CN
- China
- Prior art keywords
- seismic data
- stratum
- amplitude
- weak amplitude
- weak
- 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 56
- 239000011435 rock Substances 0.000 claims description 18
- 238000003860 storage Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000000513 principal component analysis Methods 0.000 claims description 6
- 238000012847 principal component analysis method Methods 0.000 abstract 1
- 230000015572 biosynthetic process Effects 0.000 description 28
- 238000005755 formation reaction Methods 0.000 description 28
- 238000010586 diagram Methods 0.000 description 14
- 238000004590 computer program Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 6
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000004062 sedimentation Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000000547 structure data Methods 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/32—Transforming one recording into another or one representation into another
- G01V1/325—Transforming one representation into another
-
- 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/40—Transforming data representation
- G01V2210/48—Other transforms
-
- 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
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是实施例提供的储层预测模块的结构示意图;
图8是实施例提供的电子设备的示意图;
其中:11-原始数据获取模块,111-地层特性获取单元,112-地震数据获取单元,12-拉平处理模块,13-主成分获取模块,14-弱振幅地震数据获取模块,15-储层预测模块,151-划分单元,152-振幅方差计算单元,153-反拉平单元,154-储层预测单元。
具体实施方式
下面结合附图和实施例对本申请作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释相关发明,而非对该发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与有关发明相关的部分。
图1是实施例提供的预测受层序界面控制的弱振幅岩溶层中储层的方法的流程图。如图1所示,本说明书实施例提供的前述方法包括步骤S101-S105。
S101:获取包括弱振幅岩溶层和相邻围岩层的地震数据的时间域地震数据体。
为了预测受到层序界面控制的弱振幅岩溶层的储层,需要提取地震数据总体中与弱振幅岩溶层有关的地震数据体,以及对前述弱振幅岩溶层的地震数据造成影响的、与相邻围岩层有关的地震数据体,并将二者作为一整体进行处理分析。请注意,前述技术术语“地震数据总体”指的是经过检波器接收所有的地震数据组成的地震数据体,其代表了待勘探区域的整体地层特性。
实际应用中,获取包括振幅岩溶层和相邻围岩层的地震数据的地震域地震数据体,需要识别出代表弱振幅岩溶层的地震数据,以及代表相邻围岩层的地震数据;随后再将两种数据分别提取形成时间域地震数据体。
一实施例中,相邻围岩层的的地震数据可以是上覆地层是的地震数据。
图2是一个实施例提供的获取弱振幅岩溶层和相邻围岩层的时间域地震数据体的过程。在一个实施例中,获取前述的时间域地震数据体的步骤包括S1011和S1012。
S1011:根据测井资料,确定层序地层和层序地层界面。
测井资料是对待勘探区域进行实地钻井获得的数据。测井资料可以是对钻进获得的岩芯进行特征识别获得的数据,也可以是根据其他已知的可表示地层属性的数据。以根据岩心进行特征识别获得层序地层和层序地层界面为例,因为相邻地层之间会发生较为明显的岩性突变,因此,可以通过岩相分析确定层序地层和对应的层序地层界面。
实际应用中,可以比对地层沉积走向关键剖面上的多个井的测井资料,识别区域性稳定界面特征,并将区域性稳定界面特征作为确定层序地层和层序地层界面的条件之一。具体的,使用单个钻井的测井资料可以对全井段进行高精度、连续定量的分析,获得小区块的层序地层和层序地层界面特征;使用多个钻井的测井资料可以进行比对,以分析地层侧向加积过程和形态。
S1012:根据层序地层和层序地层界面从地震数据总体中获取弱振幅岩溶层和相邻围岩层的时间域地震数据体。
步骤S1012是参照根据测井资料已经确定的地层层序界面和地层层特征,对地震数据总体进行分析,确定其中对应弱振幅岩溶层和围岩层的地震数据,并以前述地震数据组成时间域地震数据体的过程。
本领域技术人员已知,因为地震振幅数据可以反映地层层序特征,测井资料可以反映地层层序特征,所以二者之间具有对应的特征关联。因此,可以通过前述测井资料获得的层序地层和层序地层界面对地震数据总体进行相应的标定,获得不同地层的地震数据体。
具体的,根据层序地层和层序地层界面,从地震数据总体中获取弱振幅岩溶层和相邻围岩层的时间域地震数据体的过程包括步骤S201-S203。
S201:根据层序地层确定地震波传播的实深关系。
在已经确定各个层序地层的岩性的情况下,通过查找岩性的地震波传播特性、各个层序地层的厚度、地震波在不同层序地层界面的传播特性可以确定地震波在地层中传播的实深关系,即确定大体对应各个层序地层的地波数据。
S202:根据前述实深关系,地震数据总体中的波峰数据或波谷数据确定层序界面数据。
因为地震波在相邻地层交界界面处会有较强的反射或折射,因此对应地层界面处的地震波数据的振幅较大或者较小,所以根据地震数据总体中的波峰数据或者波谷数据可以确定层序界面数据;根据层序界面数据可以大体确定相邻层序地层的分界区域。
S203:根据层序地层和层序界面数据从地震数据总体中获取弱振幅岩溶层和围岩层对应的时间域地震数据体。
在确定层序界面数据后,即可以前述的层序界面数据为基础,并参照各个层序地层的岩性确定弱振幅岩溶层和围岩层分别对应的地震数据,继而获得前述地层组成的时间域地震数据体。
实际应用中,可以通过人工分析确定地层的重要层序,然后对于重要层序之间的中间年代地层采用自动追踪,获得对应的地层界面数据,继而确定不同地层对应的地震数据体。
当然,在其他实施例中,如果地震数据分析人员的经验较为丰富,能够较为准确的识别出不同地层对应的地震振幅特征,也可以由地震分析人员直接从地震数据总体中获取弱振幅岩溶层和围岩层对应的时间域地震数据体。
S102:以时间域地震数据体的层序界面数据为约束,对时间域地震数据体进行拉平处理,得到地层域地震数据体。
步骤S102中,将时间域地震数据体转换为地层域地震数据体,实际是以沉积地层的形成原理为依据,对地层进行特性还原,使得弱振幅岩溶层和围岩层的层序界面转换成平面的过程。
S103:根据地层域地震数据体获得主成分地震数据。
步骤S103是对地层域地震数据体进行处理,以确定融合在其中的代表地层特性、地层层序界面的主成分地震数据的过程。
因为地震反射同相轴的产状、相位等信息与地层和沉积过程有强相关性,而同一地层不同区域的沉积过程大概率相同,所以可以通过一小块区域既可以确定地震数据中的主成分地震数据。基于前述分析,本说明书实施例可以采用核主成分分析的方法提取地层域地震数据体中的主成分地震数据;对应的前述的主成分地震数据为核主成分地数据。
确定核主成分地震数据的方法包括:以一个地震道为中心,选取若干个相邻道形成一个子体,将子体作为一个核主成分单元,提取其中的主成分作为地震数据体的主成分核主成分地震数据。
S104:根据地层域地震数据体与主成分地震数据获得地层域弱振幅地震数据体。
步骤S104是将主成分地震数据从地层域数据体中去除,以获得剩余的地层域弱振幅地震数据体的过程。
在主成分地震数据为核主成分地震数据的情况下,步骤S104可以将地层域地震数据体按照核主成分地震数据的大小进行划分,用按照核主成分地震数据分别减去核主成分地震数据而获得弱振幅地震数据,并以前述弱振幅地震数据形成地层域弱振幅地震数据体。
S105:根据地层域弱振幅地震数据体预测弱振幅岩溶层的储层。
因为地层域弱振幅地震数据体是去除代表地层特性界面特性的数据,而仅是地层本身组成成分特性的数据,所以对其进行分析可以直接确定不同区域的差别特性,继而通过差别特性预测弱振幅岩溶层的储层。
本申请实施例中,实现步骤S105的方法可以有两种,分别为下述的方法一和方法二。
方法一
图3是一实施例提供的根据地层域弱振幅地震数据体预测弱振幅岩溶层中储层方法的流程图。如图3所示,前述预测方法包括步骤S301-S304。
S301:对地层域弱振幅地震数据体进行划分,得到地层域弱振幅地震数据子体。
步骤S301是采用一定大小的窗口,对地层域弱振幅地震数据体进行划分,形成地层域弱振幅地震数据子体的过程。
S302:计算地层域弱振幅地震数据子体内弱振幅数据的振幅方差。
S303:将地层域弱振幅地震数据子体的振幅方差恢复为时间域弱振幅地震数据子体的振幅方差。
在前述步骤S102中,将时间域地震数据进行了拉平处理,以求得地层域地震数据体;而前述的地层域弱振幅地震数据子体的振幅方差仍对应地层域地震数据体,并不能用于表征储层的实际空间分布特性,因此需要时间域地震数据体和地层域地震数据体的关系,将前述的振幅方差恢复到时间域弱振幅地震数据子体。
S304:根据时间域弱振幅数据子体的振幅方差预测弱振幅岩溶层的储层。
方法二
图4是另一实施例提供的根据地层域弱振幅地震数据体预测弱振幅岩溶层中储层方法的流程图。如图4所示,前述的方法包括步骤S401-S404。
S401:将地层域弱振幅地震数据体恢复为时间域弱振幅地震数据子体。
在前述步骤S102中,将时间域地震数据进行了拉平处理,以求得地层域地震数据体,并最终得到了地层域弱振幅地震数据体;而地层域弱振幅地震数据体并不能用于表征储层的实际空间分布特性,因此需要将地层域弱振幅地震数据体恢复为时间域弱振幅地震数据体。
S402:对时间域弱振幅地震数据子体进行划分,得到时间域弱振幅地震数据子体。
步骤S402是采用一定大小的窗口,对时间域弱振幅地震数据体进行划分,形成时间域弱振幅地震数据子体的过程。
S403:计算时间域弱振幅地震数据子体内弱振幅数据的振幅方差。
S404:根据时间域弱振幅地震数据子体的振幅方差预测弱振幅岩溶层的储层。
综合前文步骤S101-S105可知,本说明书实施例提供的预测受层序界面控制的弱振幅岩溶层中储层的方法,采用拉平处理技术将时间域地震数据体转换为地层域地震数据体,并在地层域地震数据体基础上获取和层序界面强反射相关的主成分地震数据;通过主成分地震数据和地层域地震数据体的比对确定了地层域弱振幅地震数据体,地层域弱振幅地震数据反应了弱振幅岩溶层不同区域的差异,可以用于识别其中的储层。
基于同一发明构思,本说明书实施例还提供一种预测受层序界面控制的弱振幅岩溶层中储层的装置,可以用于实现上述实施例所描述的方法,如下面的实施例。由于预测受层序界面控制的弱振幅岩溶层中储层的装置解决问题的原理与前述方法相似,因此处理地震数据的装置的具体方法可以参见前述方法内容,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的系统较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图5是实施例提供的预测受层序界面控制的弱振幅岩溶层中储层的装置的结构示意图。如图5所示,前述装置包括原始数据获取模块11、拉平处理模块12、主成分获取模块13、弱振幅地震数据获取模块14和储层预测模块15。
其中:原始数据获取模块11用于获取包括弱振幅岩溶层和相邻围岩层的地震数据的时间域地震数据体;拉平处理模块12用于以时间域地震数据体中的层序界面数据为约束,对时间域地震数据体进行拉平处理,得到地层域地震数据体;主成分获取模块13用于根据地层域地震数据体获得主成分地震数据;弱振幅地震数据获取模块14用于根据地层域地震数据体和主成分地震数据获得地层域弱振幅地震数据体;储层预测模块15用于根据地层域弱振幅地震数据体预测弱振幅岩溶层中的储层。
图6是实施例中的原始数据获取模块的结构示意图。请参见图6,在一个具体实施例中,原始数据获取模块11包括地层特性获取单元111和地震数据获取单元112。地层特性获取单元111用于根据测井资料,确定层序地层和层序地层界面;地震数据体获取单元用于根据层序地层和层序地层界面从地震数据总体中获取弱振幅岩溶层和相邻围岩层的时间域地震数据体。
在一个具体实施例中,地震数据体获取单元包括时深关系确定子单元、层序界面数据确定子单元和地震数据体获取子单元。时深关系确定子单元用于根据层序地层确定地震波传播的时深关系;层序界面数据确定子单元用于根据时深关系、地震数据中体中的波峰数据或波谷数据确定层序界面数据;地震数据体获取子单元,用于根据层序地层和层序界面数据从地震数据总体中提取弱振幅岩溶层和围岩层对应的时间域地震数据体。
在一个实施例中,主成分地震数据包括核主成分地震数据;主成分获取模块13获取主成分地震数据包括:对地层域数据做核主成分分析,得到核主成分地震数据。
图7是实施例提供的储层预测模块的结构示意图。在一个实施例中,储层预测模块15包括划分单元151、振幅方差计算单元152、反拉平单元153和储层预测单元154。划分单元151用于对地层域弱振幅地震数据体进行划分,得到地层域弱振幅地震数据子体;振幅方差计算单元152用于计算地层域弱振幅地震数据子体内弱振幅数据的振幅方差;反拉平单元153用于将地层域弱振幅地震数据子体的振幅方差恢复为时间域弱振幅地震数据子体的振幅方差;储层预测单元154用于根据时间域弱振幅地震数据子体的振幅方差预测弱振幅岩溶层的储层。
在另外一个实施例中,储层预测模块15包括反拉平单元153、划分单元151、振幅方差计算单元152和储层预测单元154。划分单元151用于对时间域弱振幅地震数据子体进行划分,得到时间域弱振幅地震数据子体;振幅方差计算单元152用于计算时间域弱振幅地震数据子体内弱振幅数据的振幅方差;储层预测单元154用于根据时间域弱振幅地震数据子体的振幅方差预测弱振幅岩溶层的储层。
除了提供前述的方法和装置外,本实施还提供一种实现前述方法的电子设备,以及存储实现前述方法程序的存储介质。
图8是实施例提供的电子设备的示意图。如图8所示,电子设备包括处理器21和存储器22,存储器22和处理器21二者电性连接。
实际应用中,存储器22可以是固态存储器,例如,只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器,其即使在断电时也保存信息,可被选择性地擦除且设有更多数据,该存储器的示例有时被称为EPROM等。存储器还可以是计算机设备领域已知的其他存储器。
在一应用中,处理器21可以加载存储在存储器22或者其他与电子设备连接的设备存储的程序,实现前文提及的预测受层序界面控制的弱振幅岩溶层中储层的方法。
请参照图8,本实施例提供的电子设备中,除了前述的处理器21和存储器22外,还包括包括输入部件23和输出部件24。
输入部件23用于获取待融合的近地表结构模型,以及各种加权参数参数。
输出部件24用于输出融合后近地表结构数据。
此外,电子设备还应当包括电源25;可能的情况下,还可以包括实现和其他电子设备联系的通信模块26。
本申请的实施例还提供一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时实现上述实施例中的预测受层序界面控制的弱振幅岩溶层中储层的方法的全部步骤,并且在执行上述方法时可以达到前述的效果。
本领域内的技术人员应明白,本发明的实施例可提供为方法、装置、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(装置)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器21以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器21执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器22中,使得存储在该计算机可读存储器22中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上描述仅为本说明书的较佳实施例以及对所运用技术原理的说明。本领域技术人员应当理解,本说明书中所涉及的发明范围,并不限于上述技术特征的特定组合而成的技术方案,同时也应涵盖在不脱离上述发明构思的情况下,由上述技术特征或其等同特征进行任意组合而形成的其它技术方案。例如上述特征与本说明书中公开的但不限于具有类似功能的技术特征进行互相替换而形成的技术方案。
Claims (6)
1.一种预测受层序界面控制的弱振幅岩溶层中储层的方法,其特征在于,包括:
根据测井资料,确定层序地层和层序地层界面;
根据所述层序地层确定地震波传播的时深关系;
根据所述时深关系、地震数据总体中的波峰数据或波谷数据确定层序界面数据;
根据所述层序地层和所述层序界面数据从所述地震数据总体中获取弱振幅岩溶层和围岩层对应的时间域地震数据体;
以所述时间域地震数据体中的层序界面数据为约束,对所述时间域地震数据体进行拉平处理,得到地层域地震数据体;
根据所述地层域地震数据体获得主成分地震数据;
根据所述地层域地震数据体和所述主成分地震数据获得地层域弱振幅地震数据体;
根据所述地层域弱振幅地震数据体预测所述弱振幅岩溶层中的储层;
其中根据所述地层域弱振幅地震数据体预测所述弱振幅岩溶层中的储层,包括:
对所述地层域弱振幅地震数据体进行划分,得到地层域弱振幅地震数据子体;
计算所述地层域弱振幅地震数据子体内弱振幅数据的振幅方差;
将所述地层域弱振幅地震数据子体的振幅方差恢复为时间域弱振幅地震数据子体的振幅方差;
根据所述时间域弱振幅地震数据子体的振幅方差预测所述弱振幅岩溶层的储层;
或包括:将所述地层域弱振幅地震数据体恢复为时间域弱振幅地震数据子体;
对所述时间域弱振幅地震数据子体进行划分,得到时间域弱振幅地震数据子体;
计算所述时间域弱振幅地震数据子体内弱振幅数据的振幅方差;
根据所述时间域弱振幅地震数据子体的振幅方差预测所述弱振幅岩溶层的储层。
2.根据权利要求1所述的方法,其特征在于,所述主成分地震数据包括核主成分地震数据;根据所述地层域地震数据体获得主成分地震数据,包括:
对所述地层域地震数据体做核主成分分析,得到所述核主成分地震数据。
3.一种预测受层序界面控制的弱振幅岩溶层中储层的装置,其特征在于,包括:
原始数据获取模块,用于获取包括弱振幅岩溶层和相邻围岩层的地震数据的时间域地震数据体;
拉平处理模块,用于以所述时间域地震数据体中的层序界面数据为约束,对所述时间域地震数据体进行拉平处理,得到地层域地震数据体;
主成分获取模块,用于根据所述地层域地震数据体获得主成分地震数据;
弱振幅地震数据获取模块,用于根据所述地层域地震数据体和所述主成分地震数据获得地层域弱振幅地震数据体;
储层预测模块,用于根据所述地层域弱振幅地震数据体预测所述弱振幅岩溶层中的储层;
所述原始数据获取模块包括:
地层特性获取单元,用于根据测井资料,确定层序地层和层序地层界面;
地震数据体获取单元,用于根据所述层序地层和所述层序地层界面从地震数据总体中获取所述弱振幅岩溶层和相邻所述围岩层的所述时间域地震数据体;
所述地震数据体获取单元包括:
时深关系确定子单元,用于根据所述层序地层确定地震波传播的时深关系;
层序界面数据确定子单元,用于根据所述时深关系、所述地震数据总体中的波峰数据或波谷数据确定层序界面数据;
地震数据体获取子单元,用于根据所述层序地层和所述层序界面数据从所述地震数据总体中获取所述弱振幅岩溶层和所述围岩层对应的所述时间域地震数据体;
所述储层预测模块包括:
划分单元,用于对所述地层域弱振幅地震数据体进行划分,得到地层域弱振幅地震数据子体;
振幅方差计算单元,用于计算所述地层域弱振幅地震数据子体内弱振幅数据的振幅方差;
反拉平单元,用于将所述地层域弱振幅地震数据子体的振幅方差恢复为时间域弱振幅地震数据子体的振幅方差;
储层预测单元,用于根据所述时间域弱振幅地震数据子体的振幅方差预测所述弱振幅岩溶层的储层;
或包括:反拉平单元,用于将所述地层域弱振幅地震数据恢复为时间域弱振幅地震数据子体;
划分单元,用于对所述时间域弱振幅地震数据子体进行划分,得到时间域弱振幅地震数据子体;
振幅方差计算单元,计算所述时间域弱振幅地震数据子体内弱振幅数据的振幅方差;
储层预测单元,用于根据所述时间域弱振幅地震数据子体的振幅方差预测所述弱振幅岩溶层的储层。
4.根据权利要求3所述的装置,其特征在于,所述主成分地震数据包括核主成分地震数据;
所述主成分获取模块获取所述主成分地震数据包括:对所述地层域地震数据体做核主成分分析,得到所述核主成分地震数据。
5.一种存储介质,其特征在于:所述存储介质存储有多条指令;所述指令适用于由处理器加载并执行如权利要求1-2任一项所述的预测受层序界面控制的弱振幅岩溶层中储层的方法。
6.一种电子设备,其特征在于,包括:存储器和处理器;
所述存储器存储有多条指令;所述指令适用于由所述处理器加载并执行如权利要求1-2任一项所述的预测受层序界面控制的弱振幅岩溶层中储层的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911298585.9A CN112987100B (zh) | 2019-12-17 | 2019-12-17 | 预测受层序界面控制的弱振幅岩溶层中储层的方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911298585.9A CN112987100B (zh) | 2019-12-17 | 2019-12-17 | 预测受层序界面控制的弱振幅岩溶层中储层的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112987100A CN112987100A (zh) | 2021-06-18 |
CN112987100B true CN112987100B (zh) | 2023-09-26 |
Family
ID=76341901
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911298585.9A Active CN112987100B (zh) | 2019-12-17 | 2019-12-17 | 预测受层序界面控制的弱振幅岩溶层中储层的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112987100B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6226596B1 (en) * | 1999-10-27 | 2001-05-01 | Marathon Oil Company | Method for analyzing and classifying three dimensional seismic information |
CN102465699A (zh) * | 2010-11-08 | 2012-05-23 | 中国石油化工股份有限公司 | 碳酸盐岩储层预测方法 |
CN103527184A (zh) * | 2013-10-28 | 2014-01-22 | 北京大学 | 一种白云岩储层的预测方法和系统 |
CN105425292A (zh) * | 2015-11-20 | 2016-03-23 | 中国石油天然气集团公司 | 一种油气预测方法及装置 |
CN105676285A (zh) * | 2016-01-23 | 2016-06-15 | 北京诺克斯达石油科技有限公司 | 一种基于pca的地震长短旋回分析法 |
CN109425899A (zh) * | 2017-09-01 | 2019-03-05 | 中国石油天然气股份有限公司 | 一种碳酸盐岩断层破碎带分布的预测方法及装置 |
CN110095812A (zh) * | 2019-05-09 | 2019-08-06 | 中国海洋石油集团有限公司 | 一种针对浅层气及底辟微裂隙影响的地震振幅补偿方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
BRPI0910752A2 (pt) * | 2008-05-22 | 2015-09-29 | Exxonmobil Upstream Res Co | métodos para transformar um volume de dados sísmicos adquirido em um exame sísmico, para produzir hidrocarbonetos de uma região de subsuperfície, para fundir superfícies e para explorar hidrocarbonetos, e, produto de programa de computador |
NO2917770T3 (zh) * | 2012-11-08 | 2018-06-02 |
-
2019
- 2019-12-17 CN CN201911298585.9A patent/CN112987100B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6226596B1 (en) * | 1999-10-27 | 2001-05-01 | Marathon Oil Company | Method for analyzing and classifying three dimensional seismic information |
CN102465699A (zh) * | 2010-11-08 | 2012-05-23 | 中国石油化工股份有限公司 | 碳酸盐岩储层预测方法 |
CN103527184A (zh) * | 2013-10-28 | 2014-01-22 | 北京大学 | 一种白云岩储层的预测方法和系统 |
CN105425292A (zh) * | 2015-11-20 | 2016-03-23 | 中国石油天然气集团公司 | 一种油气预测方法及装置 |
CN105676285A (zh) * | 2016-01-23 | 2016-06-15 | 北京诺克斯达石油科技有限公司 | 一种基于pca的地震长短旋回分析法 |
CN109425899A (zh) * | 2017-09-01 | 2019-03-05 | 中国石油天然气股份有限公司 | 一种碳酸盐岩断层破碎带分布的预测方法及装置 |
CN110095812A (zh) * | 2019-05-09 | 2019-08-06 | 中国海洋石油集团有限公司 | 一种针对浅层气及底辟微裂隙影响的地震振幅补偿方法 |
Non-Patent Citations (2)
Title |
---|
Seismic characterization of a carbonate reservoir in Tarim Basin;Liu Y et al.;Geophysics;第85卷(第5期);B177-B188 * |
基于层控蚁群层位自动追踪的Wheeler转换方法及其应用;殷文 等;中国石油大学学报(自然科学版);第41卷(第01期);51-59 * |
Also Published As
Publication number | Publication date |
---|---|
CN112987100A (zh) | 2021-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA3076280C (en) | Method and apparatus for determining thin interlayers | |
US10859725B2 (en) | Resource production forecasting | |
Rovelli et al. | Can seismic waves be trapped inside an inactive fault zone? The case study of Nocera Umbra, central Italy | |
CN102736107B (zh) | 一种能量约束非均质储层厚度识别系统 | |
NO337423B1 (no) | Fremgangsmåte for oppdatering av en geologisk modell med seismiske data | |
CN104755960A (zh) | 基于盆地建模对用于处理地震数据的速度模型的改进 | |
CN104502966A (zh) | 一种薄储层的预测方法及系统 | |
US20140336940A1 (en) | Estimation of q-factor in time domain | |
Cultrera et al. | Azimuth‐dependent amplification of weak and strong ground motions within a fault zone (Nocera Umbra, central Italy) | |
Calderoni et al. | Investigating spectral estimates of stress drop for small to moderate earthquakes with heterogeneous slip distribution: Examples from the 2016–2017 Amatrice earthquake sequence | |
Klin et al. | ER3D: a structural and geophysical 3-D model of central Emilia-Romagna (northern Italy) for numerical simulation of earthquake ground motion | |
CN107884830A (zh) | 沉积地层旋回划分方法及装置 | |
CN112987100B (zh) | 预测受层序界面控制的弱振幅岩溶层中储层的方法和装置 | |
CN112505754A (zh) | 基于高精度层序格架模型的井震协同划分沉积微相的方法 | |
CN112130209B (zh) | 岩溶储层预测方法及装置 | |
CN113806674B (zh) | 古河道纵向尺度的量化方法、装置、电子设备及存储介质 | |
Chin et al. | Earthquakes and seismic hazard in southern New Caledonia, southwest Pacific | |
CN111075440B (zh) | 不整合油气藏规模预测方法以及装置 | |
Pendrel et al. | Estimating porosity from 3D seismic inversion and 3D geostatistics | |
CN114859407A (zh) | 火山岩储层声学特征参数的确定方法及装置 | |
Lancieri et al. | Simulated shaking maps for the 1980 Irpinia earthquake, Ms 6.9: Insights on the observed damage distribution | |
CN111815769A (zh) | 逆冲推覆构造带构造的建模方法、计算设备及存储介质 | |
Mantovani et al. | Seismic noise-based strategies for emphasizing recent tectonic activity and local site effects: The Ferrara Arc, Northern Italy, case study | |
Weimann et al. | Unconventional Variables that Impact on the Fracture Height Growth on Vaca Muerta Formation and its Implications for Fracture Simulation and Well Productivity | |
Salem | Subsidence mechanisms and tectonic control on seismic stratigraphy and reservoir characterisation of the Kotla Graben and adjacent highs, West Central Sirt Basin, Libya |
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 |