CN112965109A - 一种强反射屏蔽下的储层的预测方法及装置 - Google Patents
一种强反射屏蔽下的储层的预测方法及装置 Download PDFInfo
- Publication number
- CN112965109A CN112965109A CN202110164198.7A CN202110164198A CN112965109A CN 112965109 A CN112965109 A CN 112965109A CN 202110164198 A CN202110164198 A CN 202110164198A CN 112965109 A CN112965109 A CN 112965109A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- curve
- strong reflection
- frequency
- seismic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 79
- 230000004044 response Effects 0.000 claims abstract description 21
- 238000001228 spectrum Methods 0.000 claims description 57
- 238000001914 filtration Methods 0.000 claims description 22
- 238000010586 diagram Methods 0.000 claims description 21
- 230000015654 memory Effects 0.000 claims description 15
- 239000004576 sand Substances 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 230000035945 sensitivity Effects 0.000 claims description 7
- 238000010835 comparative analysis Methods 0.000 claims description 3
- 239000003208 petroleum Substances 0.000 abstract description 2
- 229910052799 carbon Inorganic materials 0.000 description 12
- 230000008859 change Effects 0.000 description 8
- 239000003610 charcoal Substances 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 125000005626 carbonium group Chemical group 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 230000005856 abnormality Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003415 peat Substances 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
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
本发明提供了一种强反射屏蔽下的储层的预测方法及装置,属于石油勘探开发技术领域。该方法包括:确定储层的地震响应特征完全被其下方强反射界面的强反射所覆盖时,储层与其下方强反射界面之间的距离λ;根据距离λ和储层的平均速度v,获得频率fd;利用频率fd进行AC曲线重构,获得重构后的AC曲线;利用重构后的AC曲线和DEN曲线得到波阻抗曲线;利用波阻抗曲线进行反演得到地震反演数据体,根据所述地震反演数据体进行储层预测。本发明能够准确识别地震响应特征完全淹没在其下方强反射界面的强反射中的储层。
Description
技术领域
本发明涉及一种强反射屏蔽下的储层的预测方法及装置,属于石油勘探开发技术领域。
背景技术
现有技术中,常采用地震振幅来进行储层预测,原理是目的层储层与其上覆、下覆泥岩地层的波阻抗差异,会造成目的层储层在地震剖面振幅上形成异常,进行地震解释时,这些振幅异常被认为是油气聚集的标志。但是,当目的层储层下方临近一个强的波阻抗界面时,将会在油藏附近形成一个极强的地震反射波异常,目的层储层的反射波会受到这个极强的反射波的干扰被屏蔽,使得目的层储层形成的振幅异常变得难以识别,也更无法识别储层段内的单砂体。因此采用常规的地震振幅方法已无法有效识别这种强反射屏蔽背景下的储层的分布范围。
目前也有针对这种强反射屏蔽背景下的储层进行预测的方法,有些是在强反射模型指导下剔除强反射屏蔽,例如申请公布号为CN106199715A的发明专利申请文件中,公开了一种基于模型的强反射分离方法和装置,该方法通过建立强反射模型的正演数据,从原始地震数据上减除正演的强反射振幅,得到储层段的地震数据体;有些是利用多子波分解与重构方法在频率域把地震数据分解为多个不同频率的数据体,优选出能够反映薄储层的高频数据体,例如申请公布号为CN107132573A的发明专利申请文件中,公开了一种应用子波分解重构技术识别强阻抗屏蔽下岩性油藏的方法,在频率域对子波进行筛选,优选参数重构出新的地震道,得到子波分解重构之后的地震数据体,在此基础上完成储层预测。
上述方案中都是通过优选子波来得到分离强反射的地震数据体,不同的是前一方案是通过寻找最佳子波得到强反射的地震响应,然后从原始地震数据中直接减去强反射的地震响应得到储层段的地震数据体;后一方案是通过寻找与储层对应频率的子波,直接得到反映储层的地震数据体。这两种方案均存在以下问题:得到的地震数据体只是对地震数据体进行数学公式的运用,并没有充分利用井上的信息,导致最终得到的储层预测精度低。
发明内容
本发明的目的在于提供一种强反射屏蔽下的储层的预测方法及装置,用以解决目前强反射屏蔽背景下的储层预测精度低的问题。
为实现上述目的,本发明提供了一种强反射屏蔽下的储层的预测方法,该方法包括以下步骤:
(1)获取目的层段的地震数据、测井数据和已钻井资料;
(2)根据所述地震数据、测井数据和已钻井资料,分析储层与其下方强反射界面之间的距离不同时,其下方强反射界面的地震反射波对储层的地震反射波的影响程度,确定储层的地震响应特征完全被其下方强反射界面的强反射所覆盖时,储层与其下方强反射界面之间的距离λ;
(3)根据所述距离λ和储层的平均速度v,获得频率fd;
(4)利用所述频率fd进行AC曲线重构,获得重构后的AC曲线,重构过程为:
对AC曲线的频率谱进行低通滤波处理,得到AC曲线频率谱的低频部分,其中,低通滤波的截止频率为频率fd;
选取一种能够区分岩性的测井曲线作为敏感曲线,对所述敏感曲线的频率谱进行高通滤波处理,得到敏感曲线频率谱的高频部分,其中,高通滤波的截止频率为频率fd;
根据AC曲线频率谱的低频部分和敏感曲线频率谱的高频部分,得到重构后的AC曲线频率谱,进而得到重构后的AC曲线;
(5)利用所述重构后的AC曲线和DEN曲线得到波阻抗曲线;
(6)利用所述波阻抗曲线进行反演得到地震反演数据体,根据所述地震反演数据体进行储层预测。
本发明还提供了一种强反射屏蔽下的储层的预测装置,该装置包括处理器和存储器,所述处理器执行由所述存储器存储的计算机程序,以实现上述的强反射屏蔽下的储层的预测方法。
本发明的有益效果是:本发明针对地震响应特征完全淹没在其下方强反射界面强反射中的储层,首先,确定了储层的地震响应特征完全被其下方强反射界面的强反射所覆盖时,储层与其下方强反射界面之间的距离λ以及相应的频率fd;然后,将频率fd作为低通滤波和高通滤波的截止频率进行AC曲线重构,由于重构后的AC曲线既包含原AC曲线的低频部分又包含敏感曲线的高频部分,因此重构后的AC曲线既能反映波阻特征又能反映储层;最后,利用重构后的AC曲线和DEN曲线得到波阻抗曲线,并基于该波阻抗曲线进行反演得到地震反演数据体,该地震反演数据体就是将储层从其下方强反射界面强反射中分离出来后的地震反演数据体,从而根据该地震反演数据体就能准确识别储层,也就是说,本发明在纵向上充分利用测井曲线构建出不改变波阻特征且能突出反映储层内单一砂体的AC曲线,进而基于构建的AC曲线得到波阻抗曲线,然后横向上充分利用波阻抗曲线的地震响应与地震波形的对应关系进行外推,从而预测储层的平面分布,储层预测精度高。
进一步地,在上述强反射屏蔽下的储层的预测方法及装置中,获得所述距离λ的过程为:
利用测井数据生成合成地震记录,按照使合成地震记录的波形与井旁地震道的波形达到最佳相似的原则,生成一个井旁道子波;
根据已钻井资料分别构建储层的波阻抗模型、储层下方强反射界面的波阻抗模型、并通过改变储层与其下方强反射界面之间的距离,分别构建不同距离下储层与其下方强反射界面叠加的波阻抗模型;
将上述每个波阻抗模型分别与所述井旁道子波进行褶积运算,得到各波阻抗模型对应的波形反射图,通过对各波形反射图进行对比分析,得到储层的地震响应特征完全被其下方强反射界面的强反射所覆盖时,储层与其下方强反射界面之间的距离λ。
进一步地,在上述强反射屏蔽下的储层的预测方法及装置中,所述频率fd为:fd=v/4λ。
进一步地,在上述强反射屏蔽下的储层的预测方法及装置中,得到重构后的AC曲线频率谱的过程为:首先将敏感曲线频率谱的高频部分转换到AC曲线频率谱的高频部分的值域范围内,得到重构后的AC曲线频率谱的高频部分,然后将AC曲线频率谱的低频部分与重构后的AC曲线频率谱的高频部分进行叠加,得到重构后的AC曲线频率谱。
进一步地,在上述强反射屏蔽下的储层的预测方法及装置中,通过对测井曲线进行直方图分析,选取一种能将砂岩和泥岩区分开来的测井曲线作为所述敏感曲线。
进一步地,在上述强反射屏蔽下的储层的预测方法及装置中,所述敏感曲线为GR曲线。
进一步地,在上述强反射屏蔽下的储层的预测方法及装置中,采用最小平方法生成所述井旁道子波。
进一步地,在上述强反射屏蔽下的储层的预测方法及装置中,利用所述波阻抗曲线进行波形指示反演得到地震反演数据体,在地震反演数据体上追踪解释单砂体油藏对应的同相轴,形成一个解释层位,沿解释层位提取地震反演数据体的振幅信息,形成振幅分布平面图;根据已钻井信息,对振幅分布平面图进行综合地质解释,得到单砂体油藏的平面展布范围。
附图说明
图1是本发明方法实施例中的强反射屏蔽下的储层的预测方法流程图;
图2-1是本发明方法实施例中井旁道子波的时间与振幅的关系图;
图2-2是本发明方法实施例中井旁道子波的振幅与频率的关系图;
图2-3是本发明方法实施例中井旁道子波的周期与频率的关系图;
图3-1是本发明方法实施例中储层的波阻抗模型、石炭系的波阻抗模型以及两者叠加的波阻抗模型示意图;
图3-2是本发明方法实施例中各不同距离(3米、8米、12米、16米和22米)下的储层、石炭系以及两者叠加的波形反射图;
图4-1是本发明方法实施例中GR曲线的直方图分析结果图;
图4-2是本发明方法实施例中DEN曲线的直方图分析结果图;
图4-3是本发明方法实施例中IMP曲线的直方图分析结果图;
图4-4是本发明方法实施例中AC曲线的直方图分析结果图;
图5-1是本发明方法实施例中过A井的地震剖面图;
图5-2是本发明方法实施例中过A井的反演剖面图;
图6-1是本发明方法实施例中的地震最大波谷振幅属性图;
图6-2是本发明方法实施例中的反演最大波谷振幅属性图;
图7是本发明装置实施例中的强反射屏蔽下的储层的预测装置结构示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。
方法实施例:
本实施例的强反射屏蔽下的储层的预测方法(以下简称储层预测方法),如图1所示,该方法包括以下步骤:
步骤1、获取目的层段的地震数据、测井数据和已钻井资料;
步骤2、根据获取的地震数据、测井数据和已钻井资料,分析储层与其下方强反射界面之间的距离不同时,其下方强反射界面的地震反射波对储层的地震反射波的影响程度,确定储层的地震响应特征完全被其下方强反射界面的强反射所覆盖时,储层与其下方强反射界面之间的距离λ;具体包括以下步骤:
(1)利用测井数据计算得到反射系数序列,对反射系数序列使用理论子波进行褶积运算生成合成地震记录,将合成地震记录的波形与井旁地震道的波形进行对比,按照使合成地震记录的波形与井旁地震道的波形达到最佳相似的原则,生成一个井旁道子波。
其中,井旁道子波采用最小平方法生成,具体为:式中,t0、ts分别为目的层储层的起止时间,Ri(t)为实际井旁地震道的波形,Si(t)为反射系数序列,bi(t)为子波,*为褶积符号,通过令E取最小值导出求井旁道子波的方程,再求出井旁道子波。
(2)根据已钻井资料分别构建储层的波阻抗模型、储层下方强反射界面的波阻抗模型、并通过改变储层与其下方强反射界面之间的距离,分别构建不同距离下储层与其下方强反射界面叠加的波阻抗模型;
(3)将上述每个波阻抗模型分别与步骤(1)中获得的井旁道子波进行褶积运算,得到各波阻抗模型对应的波形反射图,通过对各波形反射图进行对比分析,得到储层的地震响应特征完全被其下方强反射界面的强反射所覆盖时,储层与其下方强反射界面之间的距离λ。
本实施例中,强反射界面是指强波阻抗界面,即界面上下存在着较大的波阻抗差异。由于波阻抗通过密度和速度计算得到,一般当界面上下速度差在2000左右时,即认为该界面是强反射界面。例如石炭系基岩(凝灰岩)的速度、密度值很高,沉积岩(砂岩、泥岩)的速度、密度值相对比较低。
下面以储层下方的强反射界面为石炭系且储层与石炭系之间为泥岩为例,对步骤(2)和步骤(3)进行详细说明:
首先对已钻井资料进行统计分析,得到已知油藏目的层储层与其下方石炭系之间的距离,以及目的层储层上下的砂岩、泥岩、石炭系信息,具体包括:砂岩(包括油层和水层)的平均速度、平均密度、平均厚度;上下泥岩的平均速度、平均密度、平均厚度;石炭系的平均速度、平均密度;
其次,根据上述统计得到的参数,分别构建储层的波阻抗模型(即砂岩的波阻抗模型)、石炭系的波阻抗模型、并通过改变储层与其下方石炭系之间的距离(即改变储层与石炭系之间泥岩的厚度),分别构建不同距离下(即不同泥岩厚度下)储层与其下方石炭系叠加的波阻抗模型(简称叠加的波阻抗模型),例如设定5个不同的泥岩厚度,每个泥岩厚度对应一个叠加的波阻抗模型,从而构建得到5个叠加的波阻抗模型,这5个模型对应的泥岩厚度各不相同;
接着,将储层的波阻抗模型与步骤(1)中得到的井旁道子波进行褶积运算,得到储层的波阻抗模型对应的波形反射图(简称储层的波形反射图);将石炭系的波阻抗模型与步骤(1)中得到的井旁道子波进行褶积运算,得到石炭系的波阻抗模型对应的波形反射图(简称石炭系的波形反射图);将各不同泥岩厚度下叠加的波阻抗模型分别与步骤(1)中得到的井旁道子波进行褶积运算,得到各不同泥岩厚度下叠加的波阻抗模型对应的波形反射图(简称叠加的波形反射图);
由于储层和石炭系的波阻抗模型参数不会发生变化,因此储层和石炭系的波形反射图不会发生变化,而叠加的波阻抗模型会随着储层与石炭系之间泥岩厚度的变化而相应改变,因此叠加的波形反射图也会随着泥岩厚度的变化而响应改变,据此,结合储层和石炭系的波形反射图,对每种泥岩厚度对应的叠加的波形反射图进行分析,就能看出随着储层与石炭系间泥岩厚度的变化,储层下方石炭系的强反射波形对储层的反射波形的干涉程度,进而能得到储层的地震响应特征完全被其下方石炭系的强反射所覆盖时,储层与其下方石炭系之间的距离λ。
步骤3、根据距离λ和储层的平均速度v,获得频率fd;其中,fd=v/4λ。
频率fd为储层的地震响应特征完全被其下方强反射界面的强反射所覆盖时,储层与其下方强反射界面之间的距离λ对应的频率。
本实施例中,根据地震分辨率的极限为四分之一波长原则,令fd=v/4λ。
步骤4、利用频率fd进行AC曲线重构,获得重构后的AC曲线,使重构后的AC曲线既能反映波阻特征又能反映储层;重构过程具体如下:
(1)对AC曲线的频率谱进行低通滤波处理,得到AC曲线频率谱的低频部分,其中,低通滤波的截止频率为频率fd;
具体地,对AC曲线进行傅里叶变换(FT),得到AC曲线的频率谱F(f)AC,利用低通滤波器公式:进行过滤运算得到AC曲线频率谱的低频部分F(f)AClow=C1(f)×F(f)AC,其中,C1(f)为低通滤波函数,f为频率,fd为步骤3中获得的频率。
(2)通过对测井曲线(例如声波时差曲线(即AC曲线)、密度曲线(即DEN曲线)、自然电位曲线(即SP曲线)、自然伽马曲线(即GR曲线)、IMP曲线(即波阻抗曲线)等)进行直方图分析或者交汇图分析,选取一种能够区分储层的测井曲线作为敏感曲线,对敏感曲线的频率谱进行高通滤波处理,得到敏感曲线频率谱的高频部分,其中,高通滤波的截止频率为频率fd;
具体地,对敏感曲线进行傅里叶变换(FT),得到敏感曲线的频率谱F(f)MG,利用高通滤波器公式:进行过滤运算得到敏感曲线频率谱的高频部分F(f)MGhigh=C2(f)×F(f)MG,其中,C2(f)为高通滤波函数;f为频率;fd为步骤3中获得的频率。
(3)根据AC曲线频率谱的低频部分F(f)AClow和敏感曲线频率谱的高频部分F(f)MGhigh,得到重构后的AC曲线频率谱F(f)ACls;
具体为:通过拟合得到AC曲线频率谱的高频部分与敏感曲线频率谱的高频部分之间的拟合关系,然后利用该拟合关系将敏感曲线频率谱的高频部分转换到AC曲线频率谱的高频部分的值域范围内,得到重构后的AC曲线频率谱的高频部分,然后将AC曲线频率谱的低频部分与重构后的AC曲线频率谱的高频部分进行叠加,得到重构后的AC曲线频率谱。
(4)将重构后的AC曲线频率谱F(f)ACls转换到时间域,得到重构后的AC曲线ACgj。
具体地,对重构后的AC曲线频率谱F(f)ACls开展反傅里叶变换(LFT),得到重构后的AC曲线ACgj。
步骤5、利用重构后的AC曲线ACgj和DEN曲线得到波阻抗曲线AIgj;其中,波阻抗曲线AIgj=ACgj×DEN。
步骤6、利用波阻抗曲线AIgj进行反演得到地震反演数据体,根据地震反演数据体进行储层预测。
具体地,综合使用波阻抗曲线AIgj、层位和地震数据开展波形指示反演,得到将储层从石炭系强反射中分离出来后的地震反演数据体,在该地震反演数据体上追踪解释单砂体油藏对应的同相轴,形成一个解释层位,沿解释层位提取地震反演数据体的振幅信息,在平面上对不同的振幅值赋予不同的颜色,形成振幅分布平面图;根据已钻井信息,对振幅分布平面图进行综合地质解释,得到单砂体油藏的平面展布范围。
其中,波形指示反演方法为现有技术,包括以下步骤:在贝叶斯框架下,通过似然函数分布和先验分布得到后验概率分布,并将波阻抗曲线AIgj作为目标函数,不断扰动模型参数,使后验概率分布函数最大时的解随机有效实现,取多次有效实现的均值作为期望值输出:
式中:T(x0)为未知点的值,T(xi)为波形优选的已知样本点的值,λi为第i个已知样本点对未知样本点的权重,n为优选样本本点的个数。
综上所述,本实施例的储层预测方法,对于地震响应特征完全淹没在其下方强反射界面强反射中的储层,首先,确定了储层的地震响应特征完全被其下方强反射界面的强反射所覆盖时,储层与其下方强反射界面之间的距离λ以及相应的频率fd;然后,将频率fd作为低通滤波和高通滤波的截止频率进行AC曲线重构,使重构后的AC曲线既能反映波阻特征又能反映储层;接着,利用重构后的AC曲线和DEN曲线得到波阻抗曲线,并基于该波阻抗曲线进行反演得到地震反演数据体,该地震反演数据体就是将储层从其下方强反射界面强反射中分离出来后的地震反演数据体,在地震反演数据体上追踪反映储层的同相轴,沿地震解释层位提取反演振幅值,选取振幅较强的区域作为识别出的油藏区域,预测出油藏的平面分布范围,也就是说,本实施例的储层预测方法,在纵向上充分利用测井曲线构建出不改变波阻特征且能突出反映储层内单一砂体的AC曲线,进而基于构建的AC曲线得到波阻抗曲线,然后横向上运用波形指示高分辨反演,充分利用波阻抗曲线的地震响应与地震波形的对应关系进行外推,从而预测储层的平面分布,储层预测精度高。
下面基于某工区的地震资料,验证本实施的储层预测方法的有效性。
工区三维地震资料的工区面积30km2,该工区内目的层储层下方的强反射界面为石炭系,该工区已钻有A1一口探井,B1、B2二口评价井。A1井在石炭系强反射界面上方4米处钻遇了一层2米油层和在石炭系上方7米处钻遇另一层2.2米砂体,在后期油藏滚动扩边时依据目的层对应的振幅属性在A1井高部位部署了B1、B2井,然而这二口井却失利了。证实受石炭系强反射的影响,仅根据目的层对应的振幅属性无法准确预测油层的平面分布。
下面利用本实施例的储层预测方法对被石炭系强反射淹没的储层分布范围进行预测,具体过程如下:
(1)确定储层的地震响应特征完全被其下方石炭系的强反射所覆盖时,储层与其下方石炭系之间的距离λ;
其中,生成的井旁道子波如图2-1至2-3所示,图2-1为井旁道子波的时间与振幅的关系图,图2-2为井旁道子波的振幅与频率的关系图,图2-3为井旁道子波的周期与频率的关系图。
已知油藏目的层储层与石炭系的距离,以及砂岩(包括油层和水层)、泥岩、石炭系信息的统计结果如下:
水层:平均速度3050m/s,平均密度2.18g/cm3;油层:平均速度2700m/s,平均密度2.1g/cm3;储层(包括水层与油层)平均厚度2米;泥岩:平均速度3310m/s,平均密度2.27g/cm3;水层与油层间平均厚度2米;储层与石炭系强反射界面间泥岩平均厚度4米;石炭系:平均速度5500m/s,平均密度2.5g/cm3。
根据上述参数,分别构建储层的波阻抗模型(即砂岩的波阻抗模型)、石炭系的波阻抗模型、并通过改变储层与其下方石炭系之间的距离(即改变储层与石炭系之间泥岩的厚度),分别构建泥岩厚度为3米、8米、12米、16米、22米的5组叠加的波阻抗模型,见图3-1;然后,将各波阻抗模型分别与生成的井旁道子波进行褶积运算,得到储层的波形反射图(见图3-2的第一列)、石炭系的波形反射图(见图3-2的第二列)、以及各不同泥岩厚度下叠加的波形反射图(见图3-2的第三列)。
从图3-2中可以看出:在石炭系不整合上的单个好储层,只有储层底界与石炭系顶部之间的距离达到22米时,上部好的储层的地震响应特征才不会受到石炭系强反射的影响;而储层底界与石炭系顶部的距离小于12米时,储层的地震响应特征完全被石炭系的强反射所覆盖,由此确定距离λ=12米。
本实施例的工区中A井的储层距石炭系只有4米,小于距离λ=12米,因此储层产生的地震反射波形会淹没于石炭系强反射中,仅根据地震振幅属性无法刻画储层展布。
(2)根据距离λ=12米和储层的平均速度v=2700m/s,获得频率fd=v/4λ=56Hz。
(3)利用频率fd进行AC曲线重构,获得重构后的AC曲线,具体如下:
如图4-1至图4-4所示,对测井曲线GR、IMP、DEN、AC进行直方图分析,图中,颜色较浅的线条代表砂岩,颜色较深的线条代表泥岩,结合图4-1至图4-4可以看出,GR曲线能够将砂岩和泥岩区分开来,因此选取GR曲线作为敏感曲线;
将fd=56Hz作为低通滤波的截止频率,对AC曲线频率谱F(f)AC进行低通滤波,得到AC曲线频率谱的低频部分F(f)AClow;将fd=56Hz作为高通滤波的截止频率,对GR曲线频率谱F(f)GR进行高通滤波,得到GR曲线频率谱的高频部分F(f)GRhigh;
根据AC曲线频率谱的低频部分F(f)AClow和GR曲线频率谱的高频部分F(f)GRhigh得到重构后的AC曲线频率谱F(f)ACls;对重构后的AC曲线频率谱F(f)ACls开展反傅里叶变换(LFT),得到重构后的AC曲线ACgj。
(5)根据公式得到波阻抗曲线AIgj=ACgj×DEN。
(6)利用波阻抗曲线AIgj进行反演得到地震反演数据体,根据地震反演数据体进行储层预测。
过A井的地震剖面图如图5-1所示,利用本实施例的储层预测方法获得的过A井的反演剖面图如图5-2所示,结合图5-1和图5-2可以看出,本实施例的储层预测方法能够将将石炭系顶上方储层段内各储层从石炭系强反射中分离出来,并且从图中可以看出,储层信息淹没在石炭系顶强反射对应的强波谷,被分为二个同相轴,对应着二套砂体,其中油层对应于下方同相轴上。
地震最大波谷振幅属性图如图6-1所示,利用本实施例的储层预测方法获得的反演最大波谷振幅属性图如图6-2所示,从图6-2中能清晰地看出油层单砂体的平面展布范围,A井分布于油藏平面分布范围内,而评价井B1、B2在油藏平面分布范围外,图6-2预测出在图6-1上无法预测的圈闭群,并优选一个圈闭布署了C1、C2二口探井,C1、C2井也钻遇了这套油藏,进一步证实了本实施例的储层预测方法的准确性。
综上所述,本实施例基于某工区的地震资料,首先,通过改变储层与石炭系强反射界面之间的距离,建立不同的波阻抗模型,通过正演模拟分析,对于超覆在石炭系不整合上的储层,确定了当储层底界与石炭系顶部之间的距离小于12米时,储层的地震响应特征完全被石炭系的强反射所覆盖;接着,通过已钻井资料统计得到该工区中A井的储层距石炭系只有4米,进一步确定了储层产生的地震反射波形会淹没于石炭系强反射中;然后,对于该与石炭系强反射界面距离小于12米的储层,通过测井曲线直方图分析出对岩性敏感的GR曲线,根据AC曲线与GR曲线得到重构后的AC曲线,使得重构后的AC曲线在不改变波组对应关系的基础上,也能够分离出石炭系强反射淹没的薄储层信息,并对重构后的AC曲线,纵向上利用测井曲线波阻特征与井点处地震响应的波形建立对应性,横向上运用波形指示高分辨反演,利用地震波形的对应性进行外推,分辨出石炭系上方的薄油层,并在平面上显示A井油藏的平面分布,得到B1、B2在油藏平面分布范围外,不仅如此,还在预测的油藏分布图上预测出了常规振幅属性图上无法预测的圈闭群,对此预测的圈闭群,优选一个圈闭布署了C1、C2二口探井,C1、C2井也钻遇了这套油藏,进一步证实了本实施例的储层预测方法的准确性,说明本实施例中的储层预测方法能够预测被石炭系强反射淹没的储层的分布范围,并能单独准确预测砂组段内油层的信息,并能取得良好的应用效果,具有一定的可靠性。
装置实施例:
本实施例的强反射屏蔽下的储层的预测装置,如图7所示,该装置包括处理器、存储器,存储器中存储有可在处理器上运行的计算机程序,所述处理器在执行所述计算机程序时实现上述方法实施例中的方法。
也就是说,以上方法实施例中的方法应理解为可由计算机程序指令实现强反射屏蔽下的储层的预测方法的流程。可提供这些计算机程序指令到处理器,使得通过处理器执行这些指令产生用于实现上述方法流程所指定的功能。
本实施例所指的处理器是指微处理器MCU或可编程逻辑器件FPGA等的处理装置。
本实施例所指的存储器包括用于存储信息的物理装置,通常是将信息数字化后再以利用电、磁或者光学等方式的媒体加以存储。例如:利用电能方式存储信息的各式存储器,RAM、ROM等;利用磁能方式存储信息的的各式存储器,硬盘、软盘、磁带、磁芯存储器、磁泡存储器、U盘;利用光学方式存储信息的各式存储器,CD或DVD。当然,还有其他方式的存储器,例如量子存储器、石墨烯存储器等等。
通过上述存储器、处理器以及计算机程序构成的装置,在计算机中由处理器执行相应的程序指令来实现,处理器可以搭载各种操作系统,如windows操作系统、linux系统、android、iOS系统等。
Claims (9)
1.一种强反射屏蔽下的储层的预测方法,其特征在于,该方法包括以下步骤:
(1)获取目的层段的地震数据、测井数据和已钻井资料;
(2)根据所述地震数据、测井数据和已钻井资料,分析储层与其下方强反射界面之间的距离不同时,其下方强反射界面的地震反射波对储层的地震反射波的影响程度,确定储层的地震响应特征完全被其下方强反射界面的强反射所覆盖时,储层与其下方强反射界面之间的距离λ;
(3)根据所述距离λ和储层的平均速度v,获得频率fd;
(4)利用所述频率fd进行AC曲线重构,获得重构后的AC曲线,重构过程为:
对AC曲线的频率谱进行低通滤波处理,得到AC曲线频率谱的低频部分,其中,低通滤波的截止频率为频率fd;
选取一种能够区分岩性的测井曲线作为敏感曲线,对所述敏感曲线的频率谱进行高通滤波处理,得到敏感曲线频率谱的高频部分,其中,高通滤波的截止频率为频率fd;
根据AC曲线频率谱的低频部分和敏感曲线频率谱的高频部分,得到重构后的AC曲线频率谱,进而得到重构后的AC曲线;
(5)利用所述重构后的AC曲线和DEN曲线得到波阻抗曲线;
(6)利用所述波阻抗曲线进行反演得到地震反演数据体,根据所述地震反演数据体进行储层预测。
2.根据权利要求1所述的强反射屏蔽下的储层的预测方法,其特征在于,获得所述距离λ的过程为:
利用测井数据生成合成地震记录,按照使合成地震记录的波形与井旁地震道的波形达到最佳相似的原则,生成一个井旁道子波;
根据已钻井资料分别构建储层的波阻抗模型、储层下方强反射界面的波阻抗模型、并通过改变储层与其下方强反射界面之间的距离,分别构建不同距离下储层与其下方强反射界面叠加的波阻抗模型;
将上述每个波阻抗模型分别与所述井旁道子波进行褶积运算,得到各波阻抗模型对应的波形反射图,通过对各波形反射图进行对比分析,得到储层的地震响应特征完全被其下方强反射界面的强反射所覆盖时,储层与其下方强反射界面之间的距离λ。
3.根据权利要求1或2所述的强反射屏蔽下的储层的预测方法,其特征在于,所述频率fd为:fd=v/4λ。
4.根据权利要求1所述的强反射屏蔽下的储层的预测方法,其特征在于,得到重构后的AC曲线频率谱的过程为:首先将敏感曲线频率谱的高频部分转换到AC曲线频率谱的高频部分的值域范围内,得到重构后的AC曲线频率谱的高频部分,然后将AC曲线频率谱的低频部分与重构后的AC曲线频率谱的高频部分进行叠加,得到重构后的AC曲线频率谱。
5.根据权利要求1所述的强反射屏蔽下的储层的预测方法,其特征在于,通过对测井曲线进行直方图分析,选取一种能将砂岩和泥岩区分开来的测井曲线作为所述敏感曲线。
6.根据权利要求1所述的强反射屏蔽下的储层的预测方法,其特征在于,所述敏感曲线为GR曲线。
7.根据权利要求2所述的强反射屏蔽下的储层的预测方法,其特征在于,采用最小平方法生成所述井旁道子波。
8.根据权利要求1所述的强反射屏蔽下的储层的预测方法,其特征在于,利用所述波阻抗曲线进行波形指示反演得到地震反演数据体,在地震反演数据体上追踪解释单砂体油藏对应的同相轴,形成一个解释层位,沿解释层位提取地震反演数据体的振幅信息,形成振幅分布平面图;根据已钻井信息,对振幅分布平面图进行综合地质解释,得到单砂体油藏的平面展布范围。
9.一种强反射屏蔽下的储层的预测装置,其特征在于,该装置包括处理器和存储器,所述处理器执行由所述存储器存储的计算机程序,以实现如权利要求1-8任一项所述的强反射屏蔽下的储层的预测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110164198.7A CN112965109B (zh) | 2021-02-05 | 2021-02-05 | 一种强反射屏蔽下的储层的预测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110164198.7A CN112965109B (zh) | 2021-02-05 | 2021-02-05 | 一种强反射屏蔽下的储层的预测方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112965109A true CN112965109A (zh) | 2021-06-15 |
CN112965109B CN112965109B (zh) | 2023-05-12 |
Family
ID=76274784
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110164198.7A Active CN112965109B (zh) | 2021-02-05 | 2021-02-05 | 一种强反射屏蔽下的储层的预测方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112965109B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002073240A1 (en) * | 2001-03-13 | 2002-09-19 | Conoco Phillips Company | Method and process for prediction of subsurface fluid and rock pressures in the earth |
CN105223608A (zh) * | 2015-08-07 | 2016-01-06 | 中国石油大学(华东) | 一种含煤强屏蔽缝洞型储层的地震预测与描述方法 |
CN105353411A (zh) * | 2015-10-27 | 2016-02-24 | 中国石油大学(华东) | 超深层煤上储层砂体刻画与描述方法 |
CN105445800A (zh) * | 2015-11-13 | 2016-03-30 | 中国石油化工股份有限公司 | 一种厚层砂体顶部分异岩性油藏的识别方法 |
WO2016084359A1 (en) * | 2014-11-26 | 2016-06-02 | Ricoh Company, Ltd. | Imaging device, object detector and mobile device control system |
CN106707338A (zh) * | 2016-11-18 | 2017-05-24 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种强屏蔽下储层高精度预测方法 |
CN106970422A (zh) * | 2017-03-02 | 2017-07-21 | 中国石油化工股份有限公司 | 一种在三类avo亮点特征储层区识别非亮点油藏的方法 |
CN111522062A (zh) * | 2019-04-04 | 2020-08-11 | 中国海洋石油集团有限公司 | 基于火山岩屏蔽量化分析的下伏地层振幅补偿法 |
-
2021
- 2021-02-05 CN CN202110164198.7A patent/CN112965109B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002073240A1 (en) * | 2001-03-13 | 2002-09-19 | Conoco Phillips Company | Method and process for prediction of subsurface fluid and rock pressures in the earth |
WO2016084359A1 (en) * | 2014-11-26 | 2016-06-02 | Ricoh Company, Ltd. | Imaging device, object detector and mobile device control system |
CN105223608A (zh) * | 2015-08-07 | 2016-01-06 | 中国石油大学(华东) | 一种含煤强屏蔽缝洞型储层的地震预测与描述方法 |
CN105353411A (zh) * | 2015-10-27 | 2016-02-24 | 中国石油大学(华东) | 超深层煤上储层砂体刻画与描述方法 |
CN105445800A (zh) * | 2015-11-13 | 2016-03-30 | 中国石油化工股份有限公司 | 一种厚层砂体顶部分异岩性油藏的识别方法 |
CN106707338A (zh) * | 2016-11-18 | 2017-05-24 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种强屏蔽下储层高精度预测方法 |
CN106970422A (zh) * | 2017-03-02 | 2017-07-21 | 中国石油化工股份有限公司 | 一种在三类avo亮点特征储层区识别非亮点油藏的方法 |
CN111522062A (zh) * | 2019-04-04 | 2020-08-11 | 中国海洋石油集团有限公司 | 基于火山岩屏蔽量化分析的下伏地层振幅补偿法 |
Non-Patent Citations (2)
Title |
---|
刘淑华等: "声波曲线重构反演技术在储层预测中的应用", 《天然气工业》 * |
王鹏飞等: "利用正演模拟优选春光区块地震同相轴预测储层", 《石油地质与工程》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112965109B (zh) | 2023-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108802812B (zh) | 一种井震融合的地层岩性反演方法 | |
EP0896677B1 (en) | Method for inverting reflection trace data from 3-d and 4-d seismic surveys | |
EP1963886B1 (fr) | Methode d'evaluation quantitative des pressions de fluides et de detection des surpressions d'un milieu souterrain | |
EP3963374B1 (en) | System and method for lateral statistical estimation of rock and fluid properties in a subsurface formation | |
EP2075597A2 (en) | Spectral conditioning for surface seismic data | |
CN116520419B (zh) | 一种热流体裂缝通道识别方法 | |
CN107728230B (zh) | 一种基于地球物理三联技术的岩性油藏预测方法 | |
CN108508489B (zh) | 一种基于波形微变化匹配的地震反演方法 | |
CN104330826A (zh) | 一种去除复杂地表条件下多种噪音的方法 | |
CN113219531A (zh) | 致密砂岩气水分布的识别方法及装置 | |
CN111077578B (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 | |
CN112965109B (zh) | 一种强反射屏蔽下的储层的预测方法及装置 | |
CN115327627A (zh) | 多信息融合的致密砂岩气分流河道表征方法及装置 | |
CN105158803A (zh) | 一种选取优势储层的方法 | |
Khattri et al. | Seismic discriminants of stratigraphy derived from Monte Carlo simulation of sedimentary formations | |
CN111025393B (zh) | 针对含薄煤层地层的储层预测方法、装置、设备及介质 | |
CN113514884A (zh) | 一种致密砂岩储层预测方法 | |
CN113848593A (zh) | 一种定量预测含煤地层中岩浆岩侵蚀区的方法 | |
CN107831537B (zh) | 高砂地比浊积岩储层描述方法 | |
CN105989240A (zh) | 砂层预测方法及装置 | |
CN116774285B (zh) | 基于特征曲线重构的薄互层预测方法、装置、设备及介质 | |
Rotimi et al. | Characterizing Geometrical Anisotropy of Petrophysical properties in the middle Shahejie formation, Liaohe Depression, China | |
Ohakwere-Eze et al. | Petrophysical evaluation and lithology delineation using cross-plots analysis from some onshore wells in the Nigerian-delta, west Africa | |
Melani et al. | Case History The use of variational mode decomposition in assisting sedimentary cyclicity analysis: A case study from an Albian carbonate reservoir, Campos Basin, southeast Brazil |
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 |