CN109425900A - 一种地震储层预测方法 - Google Patents
一种地震储层预测方法 Download PDFInfo
- Publication number
- CN109425900A CN109425900A CN201710791183.7A CN201710791183A CN109425900A CN 109425900 A CN109425900 A CN 109425900A CN 201710791183 A CN201710791183 A CN 201710791183A CN 109425900 A CN109425900 A CN 109425900A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- seismic
- coefficient
- evaluating
- porosity
- 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
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. analysis, for interpretation, for correction
- 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/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
-
- 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
- G01V2210/6244—Porosity
Abstract
本发明涉及一种地震储层预测方法。该地震储层预测方法:1)利用地震资料,对目标层地震构造进行解释,得到待预测区目标层地震层位解释数据;根据已完钻井的储层厚度、储层孔隙度数据,定义厚度加权平均孔隙度,计算储层评价系数;2)提取地震属性,并与储层评价系数进行回归分析,筛选出具有相关性的地震属性作为优选地震属性,建立优选地震属性与储层评价系数的相关关系;3)根据相关关系将优选地震属性网格转化为储层评价系数网格,再经校正处理,即得反映待预测区储层发育分布图。该方法建立储层评价系数与地震属性之间的关系,通过地震属性优化,实现了对储层的预测;该方法提高了储层预测的精度,满足了储层评价的需要。
Description
技术领域
本发明属于地震数据的处理领域,具体涉及一种地震储层预测方法。
背景技术
地震储层预测基础理论的应用可以追溯到上个世纪初,1927年地震反射理论应用到油气勘探中,经过几十年的发展,自20世纪80年代以来,储层预测技术逐步发展成熟,在储层的岩性、空间展布形态、物性和含油性预测中发挥着重要作用。目前,地震储层预测已成为油气勘探开发的主导技术,并在实际应用中不断发展完善。
地震储层预测主要是通过分析地震波的速度、振幅、相位、频率、波形等参数的变化来预测储集岩层的分布范围、储层特征等。岩性厚度、储层物性和充填在其中的流体性质的空间变化,造成了地震反射波速度、密度、波阻抗、振幅、相位、频率、波形等的相应变化,这些变化即是目前储层地震预测的主要依据。
目前国内外的地震储层预测方法主要有三类:(Ⅰ)基于地震反演的储层预测方法;(Ⅱ)基于地震属性分析的储层预测方法;(Ⅲ)AVO分析技术。储层预测的内容包括:①预测储集体的几何特征,如分布范围、储层厚度、储集空间大小、储集体的空间展布特征等;②预测储层物性,如孔隙度等;③预测储层含油气性。其中评价储层好坏的关键指标是储层厚度和储层孔隙度,储层厚度越大、储层孔隙度越高则储层评价越好,即储集性能越好。
申请公布号为CN104950331A的专利公开了一种砂泥岩储层的孔隙度与泥岩含量的地震预测方法,该方法是采用叠前反演,纵、横波阻抗联合反演孔隙度和泥质含量,来提高地震预测孔隙度和泥质含量的计算精度。
常规的储层预测是将储层厚度预测和孔隙度预测分别进行的,在预测过程中发现,储层密度、速度、波阻抗等地震属性的变化不仅与储层厚度变化有关,同时与孔隙度变化关系密切;也即在预测储层厚度时,其依据的地震属性变化同时也包含着储层孔隙度变化的因素,因而这些地震属性的变化并不能准确反映储层厚度的变化,预测储层孔隙度的情况与此相同。由于储层厚度和储层孔隙度均导致了地震属性的变化,只预测其中一种,预测结果往往有偏差,并且预测成果在应用中也不能直观的表达储层的储集性能。例如在目标评价和井位部署时,常常遇到这样的问题:单从预测储层厚度图上看,某区域储层厚度最大,可是从预测孔隙度图上看,储层物性好(孔隙度值较大)的又是另外一个区域;只有当某个区域其储层厚度和储层孔隙度均比较发育时,则该区域评价为储集性能优越。
由于储层预测因素不单一,往往造成储层预测的精度欠佳,不能有效指导油气藏的勘探开发部署。
发明内容
本发明的目的在于提供一种地震储层预测方法,从而解决现有预测方法采用单一因素,储层厚度、储层物性两种参数之间相互影响,不能准确评价储集性能的问题。
为实现上述目的,本发明的地震储层预测方法的技术方案是:
一种地震储层预测方法,包括以下步骤:
1)利用地震资料,对目标层地震构造进行解释,得到待预测区目标层地震层位解释数据;
根据已完钻井的储层厚度、储层孔隙度数据,根据式(1)计算厚度加权平均孔隙度,根据式(2)计算储层评价系数;
式(1)~式(3)中,为目标层厚度加权平均孔隙度,%;hi为目标层内各单层储层厚度,m;Φi为目标层内各单层储层孔隙度,P为储层评价系数,H为储层总厚度;
2)从步骤1)所得目标层地震层位解释数据提取地震属性,并与储层评价系数进行回归分析,筛选出具有相关性的地震属性作为优选地震属性,建立优选地震属性与储层评价系数的相关关系;
3)根据相关关系将优选地震属性网格转化为储层评价系数网格,再经校正处理,即得反映待预测区储层发育分布图。
厚度加权平均孔隙度是将孔隙度按厚度加权算数平均,是应用加权平均算法得到的。加权平均算法是利用若干个按照厚度排列起来的同一变量(孔隙度)的观测值,并以厚度为权数,计算出观测值的加权算术平均数。
本发明的地震储层预测方法,提出用储层评价系数来进行储层预测,其既表征厚度又体现孔隙度,将以往储层厚度、储层物性在储层评价时无法对应的关系,转化为储层评价系数与地震属性一一对应的单因素关系,提高了储层预测的精度,满足了储层评价的需要。该方法通过对待预测工区目标层常规地震层位解释数据进行层属性提取,建立储层评价系数和地震属性之间的关系,从而预测储层发育区,为油气藏的勘探开发提供了快捷、有效的技术手段,其预测吻合率高、研究周期短、成本低。
步骤2)中,优选地震属性与储层评价系数的回归系数不低于60。选取相关性好的地震属性,可建立起良好的相关关系,提高储层预测的适用性和可靠性。
步骤3)中,所述校正处理包括以下步骤:
a)利用井点处残差,使用残差网格对储层评价系数网格进行校正,得到校正后的储层评价系数网格;
b)对校正后的储层评价系数网格进行趋势面分析处理,得到反映待预测区储层发育趋势数据;
c)将待预测区储层发育趋势数据通过平滑滤波处理成图,即得所述待预测区储层发育分布图。
通过上述校正处理,可有效消除随机性因素,得到直观、有效的储层发育成果图。
根据储层发育成果图,标出储层评价系数等值线,确定油气藏有利区,完成储层预测。
附图说明
图1为本发明的地震储层预测方法的流程图;
图2为本发明实施例的合成记录标定图;
图3为本发明实施例的补偿中子-纵波阻抗-测井解释孔隙度的交汇图;
图4为本发明实施例的地震属性提取图;
图5为本发明实施例的储层厚度与提取的地震属性的相关分析图;
图6为本发明实施例的储层孔隙度与提取的地震属性的相关分析图;
图7为本发明实施例的储层评价系数与提取的地震属性的相关分析图;
图8为本发明实施例的储层评价系数与弧长地震属性的交汇图;
图9为本发明实施例中最终预测的有利储层发育区。
具体实施方式
下面结合某工区内储层预测的实施例和附图对本发明的实施方式作进一步说明。
实施例
本实施例的地震储层预测方法,流程示意图如图1所示,包括以下步骤:
1)收集待预测工区的分析化验资料以及已完钻井的钻井、测井资料,建立该区地震解释数据库;
2)分析井周围地震道数据的主频,确定子波的频率并选取零相位子波;通过调节子波的主频和相位,使合成记录与地震道达到最佳匹配,如图2所示,建立合理的井震关系,得到合成记录地震层位标定结果,从而形成地质构造解释和储层预测的基础;
3)储层地球物理特征分析
3.1储层测井相应特征
利用待测工区测井声波时差和岩石密度计算对应的纵波阻抗,通过如图3所示的补充中子-纵波阻抗-测井解释孔隙度的交汇图分析可知,利用纵波阻抗能够很好的反映储层物性变化,纵波阻抗值越低,储层孔隙度越大,物性越好;非储层具有高波阻抗,低自然伽马和高电阻特征,而储层具有低波阻抗的特征;
3.2储层地震反射特征
地震振幅和频率对储层和非储层可以很好的区分,待预测工区目的层储层的振幅响应为强振幅、变强度、连续波形,非储层石膏、泥岩为强波谷振幅;
波阻抗波形的负半周代表了储层,其波形宽窄及振幅强弱与储层厚度具有较好的对比关系;储层表现为低波阻抗,波谷最小振幅对应于储层的中心位置,且储层厚度越大波形越宽、振幅越大;相反,储层厚度变薄则振幅变弱、波形变窄;
波阻抗越低、波形越宽、振幅越强则储层孔隙度越大、储层厚度越大,储层越发育;反之,波阻抗越高、波形越窄、振幅越弱则储层孔隙度越小、储层厚度越小,储层越不发育;
4)在待预测工区内,利用建立的地震解释库,根据步骤2)所得合成记录地震层位标定结果,对目标层进行层位追踪和构造精细解释,得到地震层位解释数据,解释密度达到25m×25m;
5)统计待预测工区内测井资料可靠的已完钻井目的层各个小层的厚度和测井解释孔隙度,按照式(1)计算目标层厚度加权平均孔隙度,按照式(2)计算目标层储层评价系数:
式(1)~式(3)中,为目标层厚度加权平均孔隙度,%;hi为目标层单层储层厚度,m;Φi为单层储层孔隙度,P为储层评价系数,H为储层总厚度;
6)从目标层地震层位解释数据提取地震属性,并与储层评价系数进行交汇分析,筛选出具有相关性的地震属性作为优选地震属性,建立优选地震属性与储层评价系数的相关关系,具体包括:
6.1以目标层地震解释结果为时窗中心,依据目标层储层厚度及地震资料固有的分辨率确定时窗大小,沿目标层提取整个工区的包括最小振幅、负振幅积分、平均谷值、弧长在内多种地震属性(如图4所示);
6.2以储层厚度、储层孔隙度、储层评价系数与步骤6.1提取的多种地震属性进行交汇分析,结果分别如图5、图6、图7所示;
对比交汇分析结果,储层评价系数与提取的地震属性的回归系数(88.1)大于单独用储层厚度或储层孔隙度和提取的地震属性的回归系数(分别为66.0、65.3);
6.3选取与储层评价系数具有良好相关性的地震属性为优选地震属性,建立相关关系,由图7、图8的结果可知,储层评价系数与弧长的回归系数达到67.4,两者之间具有较好的正相关关系;
7)根据相关关系将优选地震属性网格转化为储层评价系数网格,再经校正处理,即得反映待预测区储层发育分布图,具体包括:
7.1根据相关关系将弧长网格转化为储层评价系数网格,计算井点处残差,使用残差网格校正得到校正后的储层评价系数网格;
7.2利用多项式函数拟合法或构造图多次平滑法,对校正后的储层评价系数网格进行趋势面分析处理,得到反映待预测工区目标层储层发育趋势数据;
7.3将待预测工区目标层储层发育趋势数据通过平滑滤波消除其随机因素,得到平滑后的待预测区数据,利用该数据成图,即得反映待预测工区目标层储层发育成果图;
8)利用储层发育成果图,生成待预测工区目标层储层评价系数等值线,圈定发育范围,如图9所示的;将储层评价系数不小于1.6的范围作为油气藏有利区,从而实现待预测工区目标层储层预测。
本发明的地震储层预测方法,通过提取目的层地震属性,建立储层评价系数与地震属性之间的关系,通过地震属性优化,实现了对储层的预测,为油气藏的有效勘探开发部署提供了科学依据。如图9所示,据此研究成果部署的K806井,在储层评价系数预测图中位于储层评价系数大于1.6带内,获得了较好的实钻效果,目的层钻油层达到22m,从而证实了本发明的适用性和可靠性。
Claims (4)
1.一种地震储层预测方法,其特征在于,包括以下步骤:
1)利用地震资料,对目标层地震构造进行解释,得到待预测区目标层地震层位解释数据;
根据已完钻井的储层厚度、储层孔隙度数据,根据式(1)计算厚度加权平均孔隙度,根据式(2)计算储层评价系数;
式(1)~式(3)中,为目标层厚度加权平均孔隙度,%;hi为目标层内各单层储层厚度,m;Φi为目标层内各单层储层孔隙度,P为储层评价系数,H为储层总厚度;
2)从步骤1)所得目标层地震层位解释数据提取地震属性,并与储层评价系数进行回归分析,筛选出具有相关性的地震属性作为优选地震属性,建立优选地震属性与储层评价系数的相关关系;
3)根据相关关系将优选地震属性网格转化为储层评价系数网格,再经校正处理,即得反映待预测区储层发育分布图。
2.如权利要求1所述的地震储层预测方法,其特征在于,步骤2)中,优选地震属性与储层评价系数的回归系数不低于60。
3.如权利要求1所述的地震储层预测方法,其特征在于,步骤3)中,所述校正处理包括以下步骤:
a)利用井点处残差,使用残差网格对储层评价系数网格进行校正,得到校正后的储层评价系数网格;
b)对校正后的储层评价系数网格进行趋势面分析处理,得到反映待预测区储层发育趋势数据;
c)将待预测区储层发育趋势数据通过平滑滤波处理成图,即得所述待预测区储层发育分布图。
4.如权利要求1~3任一项所述的地震储层预测方法,其特征在于,根据储层发育成果图,标出储层评价系数等值线,确定油气藏有利区,完成储层预测。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710791183.7A CN109425900B (zh) | 2017-09-05 | 2017-09-05 | 一种地震储层预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710791183.7A CN109425900B (zh) | 2017-09-05 | 2017-09-05 | 一种地震储层预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109425900A true CN109425900A (zh) | 2019-03-05 |
CN109425900B CN109425900B (zh) | 2021-02-05 |
Family
ID=65514130
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710791183.7A Active CN109425900B (zh) | 2017-09-05 | 2017-09-05 | 一种地震储层预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109425900B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109884708A (zh) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | 基于构造导向滤波的古生界潜山储层地震预测方法 |
CN112198551A (zh) * | 2020-09-12 | 2021-01-08 | 北京恺标技术发展有限公司 | 一种通用的定量化评价储层裂缝强度的新方法 |
CN112946751A (zh) * | 2021-03-17 | 2021-06-11 | 成都捷科思石油天然气技术发展有限公司 | 一种预测不同储层平面分布的方法及系统 |
CN112946754A (zh) * | 2019-12-10 | 2021-06-11 | 中国石油天然气集团有限公司 | 储层孔隙度预测方法及装置 |
CN113393335A (zh) * | 2020-03-13 | 2021-09-14 | 中国石油天然气集团有限公司 | 基于多地震属性优化的储层油气预测方法及装置 |
CN113625340A (zh) * | 2021-08-02 | 2021-11-09 | 中海石油(中国)有限公司 | 一种基于属性敏感度选择储层敏感地震属性组合的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08320381A (ja) * | 1995-05-25 | 1996-12-03 | Hitachi Ltd | 地震波伝播の表示方法 |
CN102147479A (zh) * | 2011-01-11 | 2011-08-10 | 中国海洋石油总公司 | 一种储层空间物性参数的建模方法 |
CN104360415A (zh) * | 2014-10-31 | 2015-02-18 | 中国石油化工股份有限公司 | 一种致密砂岩储层裂缝识别的方法 |
CN104678434A (zh) * | 2013-12-02 | 2015-06-03 | 中国石油化工股份有限公司 | 一种预测储层裂缝发育参数的方法 |
CN106526669A (zh) * | 2016-09-19 | 2017-03-22 | 中国石油化工股份有限公司 | 一种页岩油气藏的地震储层预测方法 |
-
2017
- 2017-09-05 CN CN201710791183.7A patent/CN109425900B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08320381A (ja) * | 1995-05-25 | 1996-12-03 | Hitachi Ltd | 地震波伝播の表示方法 |
CN102147479A (zh) * | 2011-01-11 | 2011-08-10 | 中国海洋石油总公司 | 一种储层空间物性参数的建模方法 |
CN104678434A (zh) * | 2013-12-02 | 2015-06-03 | 中国石油化工股份有限公司 | 一种预测储层裂缝发育参数的方法 |
CN104360415A (zh) * | 2014-10-31 | 2015-02-18 | 中国石油化工股份有限公司 | 一种致密砂岩储层裂缝识别的方法 |
CN106526669A (zh) * | 2016-09-19 | 2017-03-22 | 中国石油化工股份有限公司 | 一种页岩油气藏的地震储层预测方法 |
Non-Patent Citations (2)
Title |
---|
郑和荣等: "《中国中西部四大盆地碎屑岩油气成藏体系》", 31 August 2016, 地质出版社 * |
黄德利等: "肇州油田扶余油层优质储层概率分布", 《大庆石油地质与开发》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109884708A (zh) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | 基于构造导向滤波的古生界潜山储层地震预测方法 |
CN112946754A (zh) * | 2019-12-10 | 2021-06-11 | 中国石油天然气集团有限公司 | 储层孔隙度预测方法及装置 |
CN112946754B (zh) * | 2019-12-10 | 2024-03-01 | 中国石油天然气集团有限公司 | 储层孔隙度预测方法及装置 |
CN113393335A (zh) * | 2020-03-13 | 2021-09-14 | 中国石油天然气集团有限公司 | 基于多地震属性优化的储层油气预测方法及装置 |
CN112198551A (zh) * | 2020-09-12 | 2021-01-08 | 北京恺标技术发展有限公司 | 一种通用的定量化评价储层裂缝强度的新方法 |
CN112946751A (zh) * | 2021-03-17 | 2021-06-11 | 成都捷科思石油天然气技术发展有限公司 | 一种预测不同储层平面分布的方法及系统 |
CN112946751B (zh) * | 2021-03-17 | 2023-11-10 | 成都捷科思石油天然气技术发展有限公司 | 一种预测不同储层平面分布的方法及系统 |
CN113625340A (zh) * | 2021-08-02 | 2021-11-09 | 中海石油(中国)有限公司 | 一种基于属性敏感度选择储层敏感地震属性组合的方法 |
CN113625340B (zh) * | 2021-08-02 | 2023-09-19 | 中海石油(中国)有限公司 | 一种基于属性敏感度选择储层敏感地震属性组合的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109425900B (zh) | 2021-02-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104502969B (zh) | 河道砂岩性油藏的识别方法 | |
CN109061765B (zh) | 非均质薄砂岩互层油藏的圈闭评价方法 | |
CN109425900A (zh) | 一种地震储层预测方法 | |
WO2019062655A1 (zh) | 薄夹层的确定方法和装置 | |
CN102109611B (zh) | 一种利用低频地震属性预测油藏优质储层的方法 | |
CN106951660A (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
CN108802812A (zh) | 一种井震融合的地层岩性反演方法 | |
CN104657523A (zh) | 砂砾岩综合地质建模方法 | |
CN105510993A (zh) | 前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法 | |
CN103454686A (zh) | 基于地层切片的小尺度沉积相进行储层预测的方法及系统 | |
CN105445800A (zh) | 一种厚层砂体顶部分异岩性油藏的识别方法 | |
CN104502966A (zh) | 一种薄储层的预测方法及系统 | |
CN108415075A (zh) | 一种少井条件下的储层预测方法 | |
CN105527653A (zh) | 一种基于地质信息的虚拟井构建方法 | |
CN111352172B (zh) | 一种用井震联合法获取铀异常在砂体中空间分布位置的方法 | |
CN107728230B (zh) | 一种基于地球物理三联技术的岩性油藏预测方法 | |
CN107132573A (zh) | 一种应用子波分解重构技术识别强阻抗屏蔽下岩性油藏的方法 | |
CN105607120A (zh) | 一种基于时移测井的地震相约束构建初始模型的方法 | |
CN107966730A (zh) | 识别断陷湖盆陡坡带近岸水下扇沉积相带的方法 | |
Schneider et al. | Interpretation of fractured zones using seismic attributes—Case study from Teapot Dome, Wyoming, USA | |
CN106353806B (zh) | 一种检测生物礁储层连通性的方法 | |
CN105116449A (zh) | 一种弱反射储层的识别方法 | |
CN112505754B (zh) | 基于高精度层序格架模型的井震协同划分沉积微相的方法 | |
Changzi et al. | Seismic prediction of sweet spots in the Da'anzhai shale play, Yuanba area, the Sichuan Basin | |
CN115857047B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |