CN112444864A - 薄层高分辨率反演方法及存储介质 - Google Patents
薄层高分辨率反演方法及存储介质 Download PDFInfo
- Publication number
- CN112444864A CN112444864A CN201910801647.7A CN201910801647A CN112444864A CN 112444864 A CN112444864 A CN 112444864A CN 201910801647 A CN201910801647 A CN 201910801647A CN 112444864 A CN112444864 A CN 112444864A
- Authority
- CN
- China
- Prior art keywords
- target layer
- seismic
- frequency
- data volume
- data
- 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 72
- 238000003860 storage Methods 0.000 title claims abstract description 11
- 238000001228 spectrum Methods 0.000 claims abstract description 86
- 238000012545 processing Methods 0.000 claims abstract description 18
- 238000004422 calculation algorithm Methods 0.000 claims description 19
- 230000035699 permeability Effects 0.000 claims description 12
- 238000003672 processing method Methods 0.000 claims description 9
- 238000001453 impedance spectrum Methods 0.000 claims description 6
- 238000013213 extrapolation Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 238000010183 spectrum analysis Methods 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 239000011435 rock Substances 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000000605 extraction Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 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/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
- 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
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为图1中步骤S110的流程示意图。
图3为图1中步骤S120的流程示意图。
图4为图1中步骤S130的流程示意图。
图5为图1中步骤S140的流程示意图。
图6a为工业区内的两个连井W1和井W2井中测得的原始地震数据包括的谱图像。
图6b为工业区内的两个连井W1井和W2井的地震数据频谱。
图7为工业区内的两个连井中的W1井进行曲线重构前后波阻抗与孔隙度和渗透率对比效果图。
图8为工业区内的两个连井W1井和W2井的低频阻抗数据体剖面图。
图9为工业区内的两个连井W1井和W2井的地震频带阻抗数据体的剖面图。
图10为工业区内的两个连井W1井和W2井的高频阻抗数据体的剖面图。
图11为工业区内的两个连井W1井和W2井的高分辨率波阻抗数据体的剖面图。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用技术手段来解决技术问题，并达到相应技术效果的实现过程能充分理解并据以实施。本申请实施例以及实施例中的各个特征，在不相冲突前提下可以相互结合，所形成的技术方案均在本发明的保护范围之内。
另外，附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机中执行。并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
实施例一
请参阅图1，本发明提供一种薄层高分辨率反演方法，所述薄层高分辨率反演方法可以应用于电子设备，所述薄层高分辨率反演方法应用于电子设备时执行步骤S110至步骤S150。
步骤S110：对原始地震数据采用提高分辨率处理方法进行处理，并从处理后的原始地震数据中提取目的层的地震数据频谱和该目的层在零相位处的地震数据体。
需要说明的是，原始地震数据中通常包括频率、相位以及振幅谱，通过采用所述提高分辨率处理方法以获得分辨率更高的地震数据频谱。
其中，所述原始地震数据可以是叠前时间偏移的地震数据，也可以是叠前深度偏移的地震数据，所述提高分辨率处理方法可以是反褶积算法，也可以是混合相位子波有色反褶积算法，也可以是基于相位分解的混合相位子波反褶积算法，还可以是改进的混合相位子波反褶积方法，在此不做具体限定，根据实际需求进行设置即可。
可选的，在本实施例中，所述步骤S110中对原始地震数据采用提高分辨率处理方法进行处理的方式可以是，对原始地震数据采用混合相位子波反褶积方法进行处理。
通过采用上述步骤使获得的处理后的原始地震数据与原始地震数据相比频谱变宽，且分辨率更高。
请结合参阅图2，在本实施例中，上述步骤S110具体可以包括步骤S112-S116。
步骤S112：根据多个层井点确定一控制点，并基于该控制点从原始地震数据中采用新复赛谱方法提取该控制点所对应的目的层的子波。
上述步骤可以是基于用户的操作从多个层井点中选择一个层井点作为控制点，并基于所述原始地震数据获得该控制点对应的原始地震谱，基于所述控制点采用所述新复赛谱技术从所述原始地震谱中该控制点对应的提取目的层子波。
步骤S114：根据各所述层井点的原始地震谱、输入子波振幅谱和期望输出子波振幅谱以及测井数据中的测井曲线，得到目标期望输出子波的振幅谱和相位谱。
具体的，上述步骤可以是，对比各所述层井点的原始地震谱、输入子波振幅谱和期望输出子波振幅谱，并采用测井曲线的高频成分监控各层井点的原始地震谱、输入子波振幅谱和期望输出子波振幅谱以保障获得的目标期望输出子波振幅谱和相位谱在提高分辨率后，其分辨率合理，即获得的目标期望输出子波振幅和相位谱在分辨率和保真度之间达到了最佳平衡状态。
步骤S116：根据所述目的层的子波以及所述目标期望输出子波的振幅谱和相位谱，得到目的层的地震数据频谱和该目的层在零相位处的地震数据体。
具体的，上述步骤可以是，通过采用期望输出子波的振幅谱和相位谱对目的层子波进行谱整形处理得到目的层在地震数据频谱和该目的层在零相位处的地震数据体，进而使获得的目的层的地震数据频谱和该目的层在零相位处的地震数据体具有高分辨率。
步骤S120：根据所述目的层的层位数据构建构造框架模型，在所述构造框架模型内对根据测井数据得到的波阻抗曲线进行插值，从插值后的波阻抗曲线中取低于地震频带的阻抗数据，得到所述目的层的低频阻抗数据体。
其中，测井数据通常包括但不限于密度曲线、孔隙度测井曲线、渗透率测井曲线、伽玛测井曲线以及声波测井曲线等，根据测井数据得到阻抗曲线的方式可以是，对根据测井数据中的一个或多个曲线得到阻抗曲线。
具体的，在本实施例中，所述根据测井数据得到阻抗曲线的步骤具体可以是：
根据所述孔隙度测井曲线和渗透率测井曲线对所述声波测井曲线进行重构得到阻抗曲线。
因此，在本实施例中，通过利用孔隙度、渗透率等更能代表储层物性的测井曲线对声波测井曲线进行重构可以有效放大储层与围岩差异以得到阻抗曲线。此外，通过重构后获得的阻抗曲线与岩性曲线GR进行交汇分析，可以定义叠后阻抗区分不同孔隙度或者渗透率优质储层门槛值，为后期识别奠定基础。进一步地，本申请通过采用孔隙度测井曲线和渗透率测井曲线对所述声波测井曲线进行重构得到阻抗曲线，解决了现有技术中存在的仅采用测井数据中的一种曲线(例如声波测井曲线)进行波阻抗反演时，存在的薄层或者薄互层往往储层与围岩声波时差较为接近直接利用声波测井曲线来反演效果差的问题，
在所述构造框架模型内对所述波阻抗曲线进行插值的方式可以是，在所述构造框架模型内对所述波阻抗曲线采用反距离加权插值法进行插值，也可以是在所述构造框架模型内对所述波阻抗曲线采用克里金插值算法进行插值，也可以是在所述构造框架模型内对所述波阻抗曲线采用改进谢别德法进行插值，还可以是在所述构造框架模型内对所述波阻抗曲线采用移动平均法进行插值，在此不做具体限定，根据实际需求进行设置即可。
请结合图参阅图3，在本实施例中，在所述构造框架模型内对所述波阻抗曲线进行插值的步骤包括步骤S122-S124。
步骤S122：根据克里金插值算法以及阻抗空间属性在空间位置上的变化确定所述波阻抗曲线中对待插值点产生影响的影响范围。
步骤S124：为所述影响范围内的采样点获得待插值点的阻抗值，根据所述待插值点的阻抗值在所述构造框架模型内对所述波阻抗曲线进行插值，得到所述目的层的低频阻抗数据体。
通过上述步骤S122-S124，采用克里金插值算法，以实现考虑阻抗空间属性在空间位置上的变化，确定对一个待插值点有影响的距离范围，然后用该距离范围内的采样点来估计待插值点的阻抗值，在插值过程中，为了充分考虑储层发育特征及展布，可以考虑各向异性插值，即样本点沿着主变程插值半径长，垂直主变程插值半径最小，使得插值方向沿着储层主要展布方向，最终得到井插值后的低频阻抗数据体。
步骤S130：根据所述目的层的地震数据频谱与所述测井数据得到所述目的层的地震数据频谱与所述测井数据之间的匹配算子，通过利用所述匹配算子对所述目的层在零相位处的地震数据体进行反演，得到所述目的层的地震频带阻抗数据体。
其中，上述步骤具体可以包括将测井数据进行谱分析得到测井数据谱，将所述测井数据谱与所述目的层的地震数据频谱进行匹配计算得到所述目的层的地震数据频谱与所述测井数据谱之间的匹配算子，其中，通过采用所述匹配算子对所述测井数据谱进行处理后即可使处理后的测井数据谱与目的层的地震数据频谱匹配；或者通过采用所述匹配算子对所述目的层的地震数据频谱进行处理，以使处理后的地震数据频谱与所述测井数据谱匹配。根据所述匹配算子对所述零相位处的地震数据体进行反演的方式可以是进行有色反演，也可以是进行随机反演，在此不做具体限定，根据实际需求进行设置即可。
请结合参阅图4，在本实施例中，步骤S130具体包括步骤S132-步骤S136。
步骤S132：通过对所述测井数据中的波阻抗进行谱分析得到波阻抗谱。
步骤S134：基于所述目的层的地震数据频谱和所述波阻抗谱得到所述目的层的地震数据频谱与所述测井数据之间的匹配算子。
步骤S136：利用所述匹配算子对所述目的层在零相位处的地震数据体进行有色反演，得到所述目的层的地震频带阻抗数据体。
需要说明的是，有色反演作为一种无井约束反演，在提高分辨率方面并没有井的参与控制，因此其分辨率仍受限于原始地震数据的分辨率，通过采用上述步骤S132-S136，以实现通过井的资料设计滤波因子从而达到反演的目的，整个过程中无子波提取过程，无初始模型约束，因而能够较为客观的反映实际地质情况。
步骤S140：根据所述测井数据中的中高频测井信息和所述原始地震数据中目的层的中频地震信息，得到所述目的层的高频阻抗数据体。
其中，上述步骤S140具体可以是，对所述原始地震数据中目的层的中频地震信息及所述测井数据中的中高频测井信息进行贝叶斯随机反演得到高频阻抗数据体，也可以是，对所述原始地震数据中目的层的中频地震信息及所述测井数据中的中高频测井信息进行有色反演得到目的层的高频阻抗数据体，根据实际需求进行设置即可。
请结合图5，在本实施例中，步骤S140具体可以包括步骤S142-S144。
步骤S142：根据所述原始地震数据中的目的层的中频地震信息指导所述测井数据中的中高频测井信息外推。
步骤S144：采用所述地震数据波形信息替代贝叶斯随机反演算法中变差函数，以得到替代后的贝叶斯随机反演算法。
步骤S146：通过对外推后的中高频测井信息采用所述替代后的贝叶斯随机反演算法进行处理，得到目的层的高频阻抗数据体。
具体的，通过上述步骤S142-S146，以实现根据原始地震数据目的层的中频地震信息指导测井数据中的高频测井信息外推，同时用波形信息代替变差函数描述空间结构，从而确保相同反射结构或者相带内的井合理插值以实现对高频测井信息对应点曲线进行插值得到一个高频插值模型，对该高频插值模型通过随机高频扰动与地震合成记录进行约束，并修正全频带的反演，从而获得高频阻抗数据体。在本实施例中，随机反演过程中由于增加了地震信息指导解决了井间插值问题，同时通过相控克服了高频随机性，进而使得纵横向分辨率均得到了提高。
步骤S150：根据所述目的层的低频阻抗数据体、地震频带阻抗数据体以及高频阻抗数据体得到目的层的高分辨率波阻抗数据体，对高分辨率波阻抗数据体进行反演，以根据反演结果识别所述目的层是否为薄层或薄互层。
具体的，上述步骤S150可以是通过对所述目的层的低频阻抗数据体、地震频带阻抗数据体以及高频阻抗数据体进行求和目的层的高分辨率波阻抗数据体，按照频率从低到高的顺序对所述高分辨率波阻抗数据体进行反演得到反演结果，从而能够基于反演结果识别该目的层为薄层或薄互层。
需要说明的是，由于高分辨率反演遵循不同级次旋回地层分频厚度成像原则，因此按照频段从低到高分步实现高分辨率反演以得到反演结果，从而可以达到根据反演结果得到分频成像的厚度地层，从而达到识别目的层为薄层或薄互层的目的。可以理解，上述进行反演的方式可以是随机反演，也可以是有色反演，在此不作具体限定，根据实际需求进行设置即可。
本申请通过采用上述的提薄层高分辨率反演方法，以实现通过对原始地震数据预处理并提取提高分辨率后的地震数据中目的层的地震数据频谱和该目的层在零相位处的地震数据体，并在此基础上，对测井数据进行曲线重构，从而更能够反映储层物性特征，通过克里金插值并且考虑储层展布方向采用各向异性参数控制插值方向获得井插值后的目的层的低频阻抗模型；基于有色反演技术获得目的层的地震频带阻抗模型，并且反演过程无子波提取无井参与，能够客观真实反映地质体特征；通过贝叶斯随机反演获得目的层的高频波阻抗模型，并且将原始地震数据与测井数据结合利用地震波形代替变差函数从而确保井插值合理性，通过随机高频扰动与地震合成记录的约束，修正了全频带的反演结果，最后将测井低频、地震中频、随机高频反演结果求和得到具有高分辨率的目的层的高分辨率波阻抗数据体，从而能够通过按照频率由低到高的顺序对所述目的层的高分辨率波阻抗数据体进行反演从而能够基于反演结果达到识别目的层为薄层或薄互层的目的，进而便于推广和利用。
实施例二
在本实施例中，以一设定工业区内包括两个连井W1井和W2井为例进行说明。其中，W1井和W2井中测得的原始地震数据包括的谱图像如图6a，对原始地震数据采用提高分辨率处理方法进行处理得到处理后的原始图像如图6b，可见，对原始地震数据处理后得到的地震数据频谱的分辨率明显提高，有利于后期薄层和薄互层识别研究。对处理后的原始地震数据采用混合相位子波反褶积技术进行处理后得到的目的层的地震数据频谱和目的层在零相位的地震数据体。进一步地，对测井数据进行处理得到低频阻抗数据体过程中，根据W1井的测井数据中的孔隙度测井曲线和渗透率测井曲线对声波测井曲线进行重构得到阻抗曲线的重构前后波阻抗与孔隙度和渗透率对比效果图如图7所示，其中，井波阻抗曲线重构前为(IMP-DT)、后重构后为(IMP-DT_Perm)，可见，波阻抗重构后对储层物性反映更好且进一步增大了与围岩的差异，为后期反演提供测井数据基础。进一步地，对W1井和W2井采用克里金插值算法进行井间插值得到的低频阻抗数据体的剖面图如图8所示，可见，克里金插值算法能够保证波阻抗横向展布与储层发育特征基本一致，并且引入各向异性参数后，可以保证测井插值沿着储层展布方向。基于目的层的地震数据频谱和所述测井数据对目的层在零相位处的地震数据体进行有色反演得到W1井和W2井的目的层的地震频带阻抗数据体的剖面如图9所示，可见，通过有色反演方法获得了地震频带相对波阻抗，且有色反演方法无子波提取、无井约束，反演结果客观反演地质体情况。根据原始地震数据中目的层的中频地震信息及测井数据中的中高频测井信息得到W1井和W井的高频阻抗数据体的剖面如图10所示，可见在贝叶斯随机反演中加入地震信息约束后，井间插值合理且高频随机性降低。根据所述低频阻抗数据体、地震频带阻抗数据体以及高频阻抗数据体得到W1井和W2井的目的层的高分辨率波阻抗数据体如图11所示，可见，所述高分辨率波阻抗数据体的垂向分辨率明显提高，横向连续性得到了加强。对上述获得的高分辨率波阻抗数据体进行反演，反演过程遵循不同级次旋回地层分频厚度成像原则以得到反演结果，从而可以根据反演结果达到识别目的层为薄层或薄互层的目的。因此，上述的高分辨率发地震反演方法为薄层和薄互层的识别提供了技术方法和流程，并且在实际生产中也得到了实践的检验，因而具有推广应用价值。
实施例三
本发明还提供一种电子设备，所述电子设备包括存储器和处理器，所述存储器与所述处理器相互之间直接或间接的电性连接，以实现数据的传输或交互，其中，所述电子设备可以是服务器，也可以是终端设备，只要具有数据存储和处理能力即可，在此不做具体限定。
所述存储器中存储有以软件或固件(Firmware)的形式存储于所述存储器中的软件功能模块或存储介质，所述处理器通过运行存储在存储器内的软件功能模块或存储介质，从而执行各种功能应用以及数据处理，即实现本发明实施例中的薄层高分辨率反演方法。
综上，本发明提供的一种薄层高分辨率反演方法及存储介质，方法首先通过对原始数据进行提高分辨率处理，并提取目的层的地震数据频谱和零相位的地震数据体；然后通过根据测井数据重构阻抗曲线以放大储层与围岩波阻抗差异，使得波阻抗携带更多的储层信息，避免了传统波阻抗反演无法区分薄层的缺陷，并在重构阻抗曲线基础上，根据目的层的层位数据构造框架模型，并在框架模型上引入各向异性参数控制样点插值方向并进行差值，使得插值主要沿着储层发育方向得到目的层的低频阻抗数据体；通过所述目的层的地震数据频谱与所述测井数据得到所述目的层的地震数据频谱与所述测井数据之间的匹配算子，并根据匹配算子对所述零相位处的地震数据体进行有色反演获得地目的层的地震频带阻抗数据体，反演过程无子波提取、无井约束，因此能够较为客观的反映实际地质情况；通过根据所述测井数据中的中高频测井信息和所述原始地震数据中的目的层的中频地震信息进行贝叶斯随机反演获得目的层的高频阻抗数据体，反演过程中加入中频地震信息指导中高频测井信息外推，同时利用地震波形替代变差函数，保证井间插值合理性，并通过原始地震数据中目的层的中频地震信息及测井数据中的中高频测井信息随机高频扰动与地震合成记录的约束，修正全频带的反演得到高频阻抗数据体；最后将低频阻抗数据体、地震频带阻抗数据体以及高频阻抗数据体相加得到目的层的高分辨率波阻抗数据体，对高分辨率波阻抗数据体进行反演以得到反演结果，从而能够根据反演结果识别所述目的层为薄层或薄互层，由于根据反演得到的反演结果中，由于目的层遵循不同级次旋回地层分频厚度成像原则进而可以根据反演结果识别目的层为薄层和薄互层过程中，因此，所述薄层高分辨率反演方法为目的层的识别提供过了可靠的依据，进而有效保证基于所述反演结果进行识别得到的识别结果的可靠性。
虽然本发明所揭露的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。
Claims (10)
1.一种薄层高分辨率反演方法，其特征在于，所述方法包括以下步骤：
对原始地震数据采用提高分辨率处理方法进行处理，并从处理后的原始地震数据中提取目的层的地震数据频谱和该目的层在零相位处的地震数据体；
根据所述目的层的层位数据构建构造框架模型，在所述构造框架模型内对根据测井数据得到的波阻抗曲线进行插值，从插值后的波阻抗曲线中取低于地震频带的阻抗数据，得到所述目的层的低频阻抗数据体；
根据所述目的层的地震数据频谱与所述测井数据得到所述目的层的地震数据频谱与所述测井数据之间的匹配算子，通过利用所述匹配算子对所述目的层在零相位处的地震数据体进行反演，得到所述目的层的地震频带阻抗数据体；
根据所述测井数据中的中高频测井信息和所述原始地震数据中目的层的中频地震信息，得到所述目的层的高频阻抗数据体；
根据所述目的层的低频阻抗数据体、地震频带阻抗数据体以及高频阻抗数据体得到目的层的高分辨率波阻抗数据体，对高分辨率波阻抗数据体进行反演，以根据反演结果识别所述目的层是否为薄层或薄互层。
2.根据权利要求1所述的薄层高分辨率反演方法，其特征在于，所述对原始地震数据采用提高分辨率处理方法进行处理，具体包括以下步骤：
对原始地震数据采用混合相位子波反褶积方法进行处理。
3.根据权利要求2所述的薄层高分辨率反演方法，其特征在于，对原始地震数据采用混合相位子波反褶积技术进行处理，得到目的层的地震数据频谱和该目的层在零相位处的地震数据体，具体包括以下步骤：
根据多个层井点确定一控制点，并基于该控制点从原始地震数据中采用新复赛谱方法提取该控制点所对应的目的层的子波；
根据各所述层井点的原始地震谱、输入子波振幅谱和期望输出子波振幅谱以及测井数据中的测井曲线，得到目标期望输出子波的振幅谱和相位谱；
根据所述目的层的子波以及所述目标期望输出子波的振幅谱和相位谱，得到目的层的地震数据频谱和该目的层在零相位处的地震数据体。
4.根据权利要求1所述的薄层高分辨率反演方法，其特征在于，所述测井数据包括孔隙度测井曲线、渗透率测井曲线以及声波测井曲线，所述根据测井数据得到波阻抗曲线，具体包括以下步骤：
根据所述孔隙度测井曲线和渗透率测井曲线对所述声波测井曲线进行重构得到阻抗曲线。
5.根据权利要求1所述的薄层高分辨率反演方法，其特征在于，在所述构造框架模型内对根据测井数据得到的波阻抗曲线进行插值，从插值后的波阻抗曲线中取低于地震频带的阻抗数据，得到所述目的层的低频阻抗数据体，具体包括以下步骤：
根据克里金插值算法以及阻抗空间属性在空间位置上的变化确定所述波阻抗曲线中对待插值点产生影响的影响范围；
为所述影响范围内的采样点获得待插值点的阻抗值，根据所述待插值点的阻抗值在所述构造框架模型内对所述波阻抗曲线进行插值，得到所述目的层的低频阻抗数据体。
6.根据权利要求1所述的薄层高分辨率反演方法，其特征在于，根据所述目的层的地震数据频谱与所述测井数据得到所述目的层的地震数据频谱与所述测井数据之间的匹配算子，通过利用所述匹配算子对所述目的层在零相位处的地震数据体进行反演，得到所述目的层的地震频带阻抗数据体，具体包括以下步骤：
通过对所述测井数据中的波阻抗进行谱分析得到波阻抗谱；
基于所述目的层的地震数据频谱和所述波阻抗谱得到所述目的层的地震数据频谱与所述测井数据之间的匹配算子；
利用所述匹配算子对所述目的层在零相位处的地震数据体进行有色反演，得到所述目的层的地震频带阻抗数据体。
7.根据权利要求1所述的薄层高分辨率反演方法，其特征在于，根据所述测井数据中的中高频测井信息和所述原始地震数据中的目的层的中频地震信息，得到所述目的层的高频阻抗数据体，具体包括以下步骤：
根据所述原始地震数据中的目的层的中频地震信息指导所述测井数据中的中高频测井信息外推；
采用所述地震数据波形信息替代贝叶斯随机反演算法中的变差函数，以得到替代后的贝叶斯随机反演算法；
通过对外推后的中高频测井信息采用所述替代后的贝叶斯随机反演算法进行处理，得到目的层的高频阻抗数据体。
8.根据权利要求1所述的薄层高分辨率反演方法，其特征在于，根据所述目的层的低频阻抗数据体、地震频带阻抗数据体以及高频阻抗数据体得到目的层的高分辨率波阻抗数据体，具体包括以下步骤：
通过对所述低频阻抗数据体、地震频带阻抗数据体以及高频阻抗数据体进行求和，得到目的层的高分辨率波阻抗数据体。
9.根据权利要求1所述的薄层高分辨率反演方法，其特征在于，对高分辨率波阻抗数据体进行反演，具体包括以下步骤：
按照频率从低到高的顺序对所述高分辨率波阻抗数据体进行反演。
10.一种存储介质，其特征在于，所述存储介质存储有计算机程序，所述计算机程序可被一个或多个处理器执行，以实现如权利要求1至9中任意一项所述的薄层高分辨率反演方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910801647.7A CN112444864A (zh) | 2019-08-28 | 2019-08-28 | 薄层高分辨率反演方法及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910801647.7A CN112444864A (zh) | 2019-08-28 | 2019-08-28 | 薄层高分辨率反演方法及存储介质 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112444864A true CN112444864A (zh) | 2021-03-05 |
Family
ID=74740656
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910801647.7A Pending CN112444864A (zh) | 2019-08-28 | 2019-08-28 | 薄层高分辨率反演方法及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112444864A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114861428A (zh) * | 2022-04-26 | 2022-08-05 | 淮北矿业（集团）有限责任公司 | 一种深部软岩识别方法、系统及电子设备 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101201405A (zh) * | 2006-12-12 | 2008-06-18 | 中国石油天然气集团公司 | 一种提高地震数据处理分辨率的方法 |
| CN101206264A (zh) * | 2007-11-08 | 2008-06-25 | 符力耘 | 高分辨率非线性地震波阻抗反演方法 |
| CN101408624A (zh) * | 2007-10-08 | 2009-04-15 | 陶庆学 | 三维地震最佳时窗河道砂体储层预测评价技术 |
| WO2016130945A1 (en) * | 2015-02-13 | 2016-08-18 | Schlumberger Technology Corporation | Diagenetic and depositional rock analysis |
| CN108919344A (zh) * | 2018-03-30 | 2018-11-30 | 北京诺克斯达石油科技有限公司 | 适用于层状介质的分频构形反演方法 |
| CN109061764A (zh) * | 2018-09-07 | 2018-12-21 | 中国石油化工股份有限公司 | 一种分频融合波阻抗反演方法 |
| WO2019062655A1 (zh) * | 2017-09-27 | 2019-04-04 | 中国石油天然气股份有限公司 | 薄夹层的确定方法和装置 |
-
2019
- 2019-08-28 CN CN201910801647.7A patent/CN112444864A/zh active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101201405A (zh) * | 2006-12-12 | 2008-06-18 | 中国石油天然气集团公司 | 一种提高地震数据处理分辨率的方法 |
| CN101408624A (zh) * | 2007-10-08 | 2009-04-15 | 陶庆学 | 三维地震最佳时窗河道砂体储层预测评价技术 |
| CN101206264A (zh) * | 2007-11-08 | 2008-06-25 | 符力耘 | 高分辨率非线性地震波阻抗反演方法 |
| WO2016130945A1 (en) * | 2015-02-13 | 2016-08-18 | Schlumberger Technology Corporation | Diagenetic and depositional rock analysis |
| WO2019062655A1 (zh) * | 2017-09-27 | 2019-04-04 | 中国石油天然气股份有限公司 | 薄夹层的确定方法和装置 |
| CN108919344A (zh) * | 2018-03-30 | 2018-11-30 | 北京诺克斯达石油科技有限公司 | 适用于层状介质的分频构形反演方法 |
| CN109061764A (zh) * | 2018-09-07 | 2018-12-21 | 中国石油化工股份有限公司 | 一种分频融合波阻抗反演方法 |
Non-Patent Citations (8)
| Title |
|---|
| YONGQIANG MA ET AL.: "Study of Seismic Inversion Technology Forecasting Sandstone Waterborne of Coal Roof", 《ADVANCED MATERIALS RESEARCH》 * |
| 张华等: "测井数据与信噪比约束的混合相位子波反褶积", 《石油地球物理勘探》 * |
| 张猛等: "混合相位地震子波提取及应用", 《油气地球物理》 * |
| 杜伟维等: "地震波形指示反演及特征参数模拟在薄储层预测中的应用", 《工程地球物理学报》 * |
| 杨瑞召等: "利用谱蓝化和有色反演分辨薄煤层", 《天然气地球科学》 * |
| 段南: "叠前地震波形指示反演在薄互储层预测中的应用", 《地球物理学进展》 * |
| 胡玮等: "波形指示反演在超深层致密砂岩薄储层中的应用", 《地球物理学进展》 * |
| 陆基孟等, 北京：地质出版社 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114861428A (zh) * | 2022-04-26 | 2022-08-05 | 淮北矿业（集团）有限责任公司 | 一种深部软岩识别方法、系统及电子设备 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1963886B1 (fr) | Methode d'evaluation quantitative des pressions de fluides et de detection des surpressions d'un milieu souterrain | |
| US11435493B2 (en) | Enhanced waveform analysis for target modes of borehole waves | |
| CN105116449B (zh) | 一种弱反射储层的识别方法 | |
| CN108957532B (zh) | 储层预测方法及装置 | |
| Faraji et al. | Integration of core data, well logs and seismic attributes for identification of the low reservoir quality units with unswept gas in the carbonate rocks of the world’s largest gas field | |
| CN110687597B (zh) | 一种基于联合字典的波阻抗反演方法 | |
| US20190265376A1 (en) | Method of extracting intrinsic attentuation from seismic data | |
| CN112444864A (zh) | 薄层高分辨率反演方法及存储介质 | |
| CN117607995A (zh) | 利用自然伽马与波阻抗融合的拟阻抗地震岩性反演方法 | |
| Melani et al. | The use of variational mode decomposition in assisting sedimentary cyclicity analysis: A case study from an Albian carbonate reservoir, Campos Basin, southeast Brazil | |
| CN111077578B (zh) | 岩层分布预测方法和装置 | |
| CN105158803A (zh) | 一种选取优势储层的方法 | |
| CN105676285B (zh) | 一种基于pca的地震长短旋回分析法 | |
| CN105989240B (zh) | 砂层预测方法及装置 | |
| Ma et al. | Coupling spatial and frequency uncertainty analyses in reservoir modeling: example of Judy Creek Reef complex in Swan Hills, Alberta, Canada | |
| CN107831537B (zh) | 高砂地比浊积岩储层描述方法 | |
| CN113589366B (zh) | 基于全波形反演的宽频融合建模方法 | |
| CN115079256B (zh) | 深水浊积水道气藏高分辨率反演方法、装置、介质及设备 | |
| CN114428361B (zh) | 基于低频模型优化的表层岩溶缝洞体预测方法及电子设备 | |
| Zhao et al. | Land broadband seismic exploration based on adaptive vibroseis | |
| Rotimi et al. | Characterizing Geometrical Anisotropy of Petrophysical properties in the middle Shahejie formation, Liaohe Depression, China | |
| Wang et al. | Application of multiple attributes fusion technology in the Su-14 Well Block | |
| CN115728839A (zh) | 一种不整合面之上超覆砂体的精细描述方法 | |
| Pang et al. | Application of seismic attenuation analysis in prediction of carbonate reservoirs | |
| Yang et al. | Geophysical prediction technology based on organic carbon content in source rocks of the Huizhou sag, the south China sea |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210305 |
|
| RJ01 | Rejection of invention patent application after publication |