CN106812522A - 基于三维地质模型的储层非均质性研究方法 - Google Patents
基于三维地质模型的储层非均质性研究方法 Download PDFInfo
- Publication number
- CN106812522A CN106812522A CN201510868068.6A CN201510868068A CN106812522A CN 106812522 A CN106812522 A CN 106812522A CN 201510868068 A CN201510868068 A CN 201510868068A CN 106812522 A CN106812522 A CN 106812522A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- permeability
- geological model
- dimensional geological
- heterogeneity
- 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
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种基于三维地质模型的储层非均质性研究方法,该基于三维地质模型的储层非均质性研究方法包括:步骤1,根据三维地质模型的渗透率模拟结果,统计渗透率的各项特征参数;步骤2,计算单个储层的非均质性表征参数;步骤3,计算储层垂向渗透率变异系数,并根据计算结果绘制该储层垂向渗透率变异系数的平面等值图。该基于三维地质模型的储层非均质性研究方法具有创新性、实用性,利于推广,为基于储层三维地质模型的非均质性表征方法提供了切实可行的方法。
Description
技术领域
本发明涉及储层非均质性研究技术领域,特别是涉及到一种基于三维地质模型的储层非均质性研究方法。
背景技术
由于沉积、成岩以及后期构造运动等因素的变化,导致储层的岩性、物性、含油性和电性等方面在三维空间上表现出差异性,这种差异性称为储层非均质性。储层非均质性对油气的成藏以及开采开发、剩余油形成和分布等都具有重要影响。储层非均质性具体可分为层间非均质性、平面非均质性、层内非均质性和孔隙非均质性,其中层间、层内和平面非均质性属于宏观非均质性,而孔隙非均质性属于微观非均质性。在储层宏观非均质性研究方面,基于岩心、物性分析资料、测井资料等,人们通过对渗透率的非均质性分析总结出来一套系统的储层非均质性的定量表征参数,包括渗透率突进系数(Tk)、渗透率均质指数(Kp)、渗透率级差(Jk)等渗透率变异系数(Vk)等等,通过这些参数,很好地表征了储层非均质性特征,但是以上参数计算时采用的渗透率参数多为岩石物性分析参数,而一般取心数据量少、代表性差,对于一个油藏,只能得到一个渗透率变异系数(Vk)值,难以准确表征整个油藏的储层非均质性特征。
随着油藏精细地质建模技术的不断发展和完善,对储层的研究手段、工具有了革命性的进步,研究成果也向立体化、数字化的方向发展。特别是三维地质模型提供了基于网格的表征储层参数(包括孔隙度、渗透率、含油饱和度)的海量数据,成为储层非均质性的最佳表现形式。但是在非均质性表征方面,还是基于传统方法计算的传统参数,不能很好地将模型用于储层非均质性研究中去。
为了充分利用精细三维地质模型提供的海量数据,更加准确地表征储层非均质性,需要提出一套基于精细三维地质模型的储层非均质性表征方法。为此我们发明了一种新的基于三维地质模型的储层非均质性研究方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种可以方便的分析参数的平面变化规律,提供更直观、更丰富的表征的基于三维地质模型的储层非均质性研究方法。
本发明的目的可通过如下技术措施来实现:基于三维地质模型的储层非均质性研究方法,该基于三维地质模型的储层非均质性研究方法包括:步骤1,根据三维地质模型的渗透率模拟结果,统计渗透率的各项特征参数;步骤2,计算单个储层的非均质性表征参数;步骤3,计算储层垂向渗透率变异系数,并根据计算结果绘制该储层垂向渗透率变异系数的平面等值图。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,统计的渗透率的各项特征参数包括最大渗透率Kmax、最小渗透率Kmin、平均渗透率Kave这些参数。
在步骤2中,需要计算的单个储层的非均质性表征参数包括渗透率突进系数Tk、渗透率均质指数Kp、渗透率级差Jk,计算公式为:
式中,Kmax为最大渗透率、Kmin为最小渗透率,Kave为平均渗透率。
在步骤3中,对单个储层的每个平面网格,计算该网格处的垂向渗透率变异系数Vk的公式为:
式中:n为单个储层对应的纵向网格数;kj为每个网格的渗透率值;kave为n个对应网格渗透率的平均值,对单个储层的渗透率变异系数Vk通过平面等值图表征。
该基于三维地质模型的储层非均质性研究方法还包括,在步骤3之后,针对某一类储层,通过步骤1—3计算储层非均质性表征参数;并对储层层间非均质性进行分析。
在对储层层间非均质性进行分析时,首先针对每个单层,计算该层的平均渗透率;然后针对多个储层应用传统公式计算层间非均质参数。
本发明中的基于三维地质模型的储层非均质性研究方法,利用模型的海量渗透率参数,计算表征储层非均质性的各项参数,特别是对于渗透率变异系数(Vk),对于每一个平面网格都能得到一个数值,可以方便的分析该参数的平面变化规律,提供更直观、更丰富的表征方法。该方法也可以针对某一类储层(如某类沉积相)开展非均质性分析,达到精细描述储层非均质性的目的。技术思路清楚、应用简单,突破了以往仅仅根据岩心分析化验资料分析储层非均质性的限制,将三维地质模型的海量数据用于储层非均质性表征,得到了垂向渗透率变异系数(Vk)平面分布等值图,提高了储层非均质性表征精度;本次技术成果与以往相比具有更好的结果,具有创新性、实用性,利于推广,为基于储层三维地质模型的非均质性表征方法提供了切实可行的方法。
附图说明
图1为本发明的基于三维地质模型的储层非均质性研究方法的一具体实施例的流程图;
图2为本发明的一具体实施例中胜海201块Ng53储层渗透率模型的示意图;
图3为本发明的一具体实施例中储层垂向渗透率变异系数(Vk)平面分布等值图;
图4为本发明的一具体实施例中胜海201块Ng53砂体主河道相分布图;
图5为本发明的一具体实施例中胜海201块Ng53砂体主河道相渗透率变异系数分布等值图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图1所示,图1为本发明的基于三维地质模型的储层非均质性研究方法的流程图。
在步骤101,根据三维地质模型的渗透率模拟结果,统计渗透率的各项特征参数,包括最大渗透率Kmax、最小渗透率Kmin、平均渗透率Kave等参数;
在步骤102,按传统的公式计算单个储层的非均质性表征参数:包括渗透率突进系数(Tk)、渗透率均质指数(Kp)、渗透率级差(Jk)等;
在步骤103,对单个储层的每个平面网格,首先采用如下公式计算该网格处的垂向渗透率变异系数(Vk)
式中:n为单个储层对应的纵向网格数;kj为每个网格的渗透率值;kave为n个对应网格渗透率的平均值。
这样对单个储层的渗透率变异系数(Vk)可以通过平面等值图表征,改进了以往单个储层只有一个变异系数的不足。
在步骤104,针对某一类储层,也可以开展以上分析,首先通过储层类型特征,将同类储层挑选出来;通过步骤101—步骤103计算储层非均质性表征参数。
在步骤105,对储层层间非均质性的分析,首先针对每个单层,计算该层的平均渗透率;然后针对多个储层应用传统公式计算层间非均质参数,该计算过程应用了模型提供的海量渗透率数据,比单独使用岩心分析样品数据更具有代表性。
在应用本发明一具体实施例中,胜利油田胜海201区块是胜利海上一个重要区块,该区块已经实现属性模型的精细模拟,基于三维地质模型应用本发明的基于三维地质模型的储层非均质性研究方法,在胜海201块开展储层非均质性研究。图2根据三维地质模型的渗透率属性,计算最大渗透率、最小渗透率及平均渗透率,并按传统的公式计算单个储层的非均质性表征参数:包括渗透率突进系数(Tk)、渗透率均质指数(Kp)、渗透率级差(Jk)等;
图3根据模型渗透率属性按公式(1)计算储层垂向渗透率变异系数(Vk),并根据计算结果绘制该参数平面等值图。根据模型沉积微相特性,选择其中的主河道沉积相,如图4所示,计算该沉积相内储层垂向变异系数,并绘制垂向变异系数等值图,如图5所示。
本发明在胜海201块开展储层非均质性研究中获得应用,取得了较好效果。该区块设计新井33口,覆盖地质储量726×104t,如果全部完钻,每年可新增产能13.3×104t,经济效益显著。
Claims (6)
1.基于三维地质模型的储层非均质性研究方法,其特征在于,该基于三维地质模型的储层非均质性研究方法包括:
步骤1,根据三维地质模型的渗透率模拟结果,统计渗透率的各项特征参数;
步骤2,计算单个储层的非均质性表征参数;
步骤3,计算储层垂向渗透率变异系数,并根据计算结果绘制该储层垂向渗透率变异系数的平面等值图。
2.根据权利要求1所述的基于三维地质模型的储层非均质性研究方法,其特征在于,在步骤1中,统计的渗透率的各项特征参数包括最大渗透率Kmax、最小渗透率Kmin、平均渗透率Kave这些参数。
3.根据权利要求2所述的基于三维地质模型的储层非均质性研究方法,其特征在于,在步骤2中,需要计算的单个储层的非均质性表征参数包括渗透率突进系数Tk、渗透率均质指数Kp、渗透率级差Jk,计算公式为:
式中,Kmax为最大渗透率、Kmin为最小渗透率,Kave为平均渗透率。
4.根据权利要求1所述的基于三维地质模型的储层非均质性研究方法,其特征在于,在步骤3中,对单个储层的每个平面网格,计算该网格处的垂向渗透率变异系数Vk的公式为:
式中:n为单个储层对应的纵向网格数;kj为每个网格的渗透率值;kave为n个对应网格渗透率的平均值,对单个储层的渗透率变异系数Vk通过平面等值图表征。
5.根据权利要求1所述的基于三维地质模型的储层非均质性研究方法,其特征在于,该基于三维地质模型的储层非均质性研究方法还包括,在步骤3之后,针对某一类储层,通过步骤1—3计算储层非均质性表征参数;并对储层层间非均质性进行分析。
6.根据权利要求5所述的基于三维地质模型的储层非均质性研究方法,其特征在于,在对储层层间非均质性进行分析时,首先针对每个单层,计算该层的平均渗透率;然后针对多个储层应用公式计算层间非均质参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510868068.6A CN106812522A (zh) | 2015-12-01 | 2015-12-01 | 基于三维地质模型的储层非均质性研究方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510868068.6A CN106812522A (zh) | 2015-12-01 | 2015-12-01 | 基于三维地质模型的储层非均质性研究方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106812522A true CN106812522A (zh) | 2017-06-09 |
Family
ID=59107557
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510868068.6A Pending CN106812522A (zh) | 2015-12-01 | 2015-12-01 | 基于三维地质模型的储层非均质性研究方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106812522A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107590550A (zh) * | 2017-07-26 | 2018-01-16 | 长江大学 | 有关超低渗透储层油田产能评价及预测的方法 |
CN108843312A (zh) * | 2018-06-04 | 2018-11-20 | 西南石油大学 | 油田储层层内非均质性综合评价方法 |
CN110242291A (zh) * | 2019-07-02 | 2019-09-17 | 中国石油化工股份有限公司 | 一种超强非均质油气储层非均质性表征方法 |
CN113107463A (zh) * | 2021-04-22 | 2021-07-13 | 中国石油天然气集团有限公司 | 一种储层非均质性评价方法 |
CN113266333A (zh) * | 2021-06-29 | 2021-08-17 | 西北大学 | 一种通过挤入饱和co2盐水改善油砂储层渗透率的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5482116A (en) * | 1993-12-10 | 1996-01-09 | Mobil Oil Corporation | Wellbore guided hydraulic fracturing |
CN104318052A (zh) * | 2014-09-16 | 2015-01-28 | 中国石油大学(北京) | 一种低渗透油藏水驱波及系数的评价方法 |
CN104747185A (zh) * | 2015-03-19 | 2015-07-01 | 成都北方石油勘探开发技术有限公司 | 非均质油藏储层综合分类评价方法 |
-
2015
- 2015-12-01 CN CN201510868068.6A patent/CN106812522A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5482116A (en) * | 1993-12-10 | 1996-01-09 | Mobil Oil Corporation | Wellbore guided hydraulic fracturing |
CN104318052A (zh) * | 2014-09-16 | 2015-01-28 | 中国石油大学(北京) | 一种低渗透油藏水驱波及系数的评价方法 |
CN104747185A (zh) * | 2015-03-19 | 2015-07-01 | 成都北方石油勘探开发技术有限公司 | 非均质油藏储层综合分类评价方法 |
Non-Patent Citations (1)
Title |
---|
张荣军等: "《低渗透油藏有利区预测及工程设计以西峰油田延安组油藏为例》", 31 August 2013 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107590550A (zh) * | 2017-07-26 | 2018-01-16 | 长江大学 | 有关超低渗透储层油田产能评价及预测的方法 |
CN108843312A (zh) * | 2018-06-04 | 2018-11-20 | 西南石油大学 | 油田储层层内非均质性综合评价方法 |
CN108843312B (zh) * | 2018-06-04 | 2019-04-30 | 西南石油大学 | 油田储层层内非均质性综合评价方法 |
CN110242291A (zh) * | 2019-07-02 | 2019-09-17 | 中国石油化工股份有限公司 | 一种超强非均质油气储层非均质性表征方法 |
CN113107463A (zh) * | 2021-04-22 | 2021-07-13 | 中国石油天然气集团有限公司 | 一种储层非均质性评价方法 |
CN113266333A (zh) * | 2021-06-29 | 2021-08-17 | 西北大学 | 一种通过挤入饱和co2盐水改善油砂储层渗透率的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104564041B (zh) | 基于开发渗透率下限的低渗透碎屑岩储层有效性评价方法 | |
CN106812522A (zh) | 基于三维地质模型的储层非均质性研究方法 | |
CN105626058B (zh) | 一种确定储层岩溶发育程度的方法及装置 | |
CN104806232B (zh) | 一种确定碳酸盐岩储层孔隙度下限的方法 | |
CN104278991B (zh) | 盐湖相烃源岩有机碳和生烃潜量的多元测井计算方法 | |
CN103983536B (zh) | 一种利用测井曲线获得页岩气含气量的方法 | |
CN103352693B (zh) | 一种基于含油气含量的产能预测系统及其方法 | |
CN104453874B (zh) | 一种基于核磁共振的砂砾岩储层含油饱和度的计算方法 | |
CN103670384B (zh) | 一种储层分类方法及系统 | |
CN102565858B (zh) | 一种多孔介质含水饱和度的计算方法 | |
CN104819923A (zh) | 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法 | |
CN101929973A (zh) | 裂缝储层含油气饱和度定量计算方法 | |
CN107966732B (zh) | 基于空间结构导向的地震属性变化率求取方法 | |
CN103485770B (zh) | 一种基于人工神经网络求取含油饱和度的方法及系统 | |
CN109138975B (zh) | 一种基于时移测井数据的求解相渗特征曲线的新方法 | |
CN106897531A (zh) | 一种低渗透石灰岩储层渗透率的定量评价方法 | |
CN108982320B (zh) | 一种利用粒度参数进行复杂孔隙结构储层渗透率计算方法 | |
CN105064986A (zh) | 利用常规测录井资料构建储层四性关系谱的方法 | |
CN111368245B (zh) | 一种考虑多影响因素的定量计算岩石Biot系数方法 | |
CN106599377A (zh) | 基于测井资料和煤体结构参量定量划分煤体结构的方法 | |
CN109710891A (zh) | 一种基于流动单元分类计算砂岩型铀矿砂体渗透率的方法 | |
CN105487136A (zh) | 基于经验模态分解和能量熵判别的碳酸盐岩储集体测井识别方法 | |
CN112363242A (zh) | 基于测录井融合的储层流体识别方法与装置 | |
CN105114067A (zh) | 一种岩性测井相的方法 | |
CN103412335B (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170609 |
|
RJ01 | Rejection of invention patent application after publication |