CN108241785A - 一种非均质性储层饱和度场精细表征方法 - Google Patents
一种非均质性储层饱和度场精细表征方法 Download PDFInfo
- Publication number
- CN108241785A CN108241785A CN201810025475.4A CN201810025475A CN108241785A CN 108241785 A CN108241785 A CN 108241785A CN 201810025475 A CN201810025475 A CN 201810025475A CN 108241785 A CN108241785 A CN 108241785A
- Authority
- CN
- China
- Prior art keywords
- formula
- saturation
- reservoir
- saturation degree
- coefficient
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明涉及一种非均质性储层饱和度场精细表征方法,包括以下步骤:1)建立饱和度公式;2)在油田孔隙度模型和渗透率模型的基础上,结合油柱高度分布状况,代入步骤1)所得到的饱和度公式,利用地质模型建模软件即实现饱和度场模型的构建。本发明以非均质性储层的实际岩心样品毛管压力资料为基础,综合考虑地质油藏特征,深入分析含油饱和度与毛管压力、储层质量系数的关系,首先建立饱和度经验公式,其次构建饱和度场地质模型,从而实现了饱和度场的合理、定量、精细表征。
Description
技术领域
本发明涉及一种非均质性储层饱和度场精细表征方法,属于油气开发技术领域。
背景技术
储层精细描述是油气开发地质研究中的重要内容,是编制油气藏开发方案的地质基础。饱和度作为储层描述中的关键参数,对饱和度场的合理、定量、精细表征是油气藏高效开发的迫切要求。
目前构建饱和度场地质模型主要有两种方法,一是基于测井解释饱和度的直接赋值法,二是运用压汞资料的J函数法,这两类方法都有一定的适用条件,两者在非均质性储层中应用效果不理想。直接赋值法的结果易受井数的多少,以及测井资料好坏的影响,特别是在非均质储层,以此方法建立的饱和度场模型不确定性较大。J函数方法也具有一定的局限性,对于非均质性较强、油柱高度大的油气藏,往往在构造高部位会出现异常值,不适用于此类油田。因此,针对非均质性储层,需要一种更合理、定量、精细的饱和度表征方法,以获得该类油田的饱和度场。从而深化油气藏认识,更好把握开发规律,降低开发风险,提高开发效率,改善开发效果。
发明内容
针对上述问题,本发明的目的是提供一种非均质性储层饱和度场精细表征方法,该方法能够实现对饱和度场的合理、定量以及精细表征。
为实现上述目的,本发明采用以下技术方案:一种非均质性储层饱和度场精细表征方法,包括以下步骤:
1)建立饱和度公式,具体过程如下:
①建立样品饱和度与毛管压力之间的Lambda函数关系,其表达式为:
Sw=APc B+C (1)
式中,Sw表示含水饱和度;Pc表示油藏条件下的毛管压力;A、B、C为饱和度系数;
毛管压力与油柱高度之间具有下述函数关系:
Pc=H*(ρw-ρo)*g (2)
式中,H表示油柱高度;ρw、ρo分别表示水和油的密度;g为重力加速度。
将式(2)代入式(1),得到:
Sw=A[H*(ρw-ρo)*g]B+C (3)
②将储层质量系数Irq作为约束条件,对饱和度公式中的系数A、B、C进行拟合,其中Irq的计算公式为:
式中,Irq表示储层质量系数;K表示岩石渗透率;表示岩石孔隙度;
储层质量系数Irq与饱和度系数A、B、C具有良好的相关性,其相关关系表示为:
式中,a1、a2、a3、b1、b2、b3、c1、c2、c3表示储层质量系数与饱和度系数A、B、C相关关系式的系数。
③由式(3)和式(5),最终确定饱和度的公式如下:
Sw=A[H*(ρw-ρo)*g]B+C
A=a1Irq 2+a2Irq+a3
B=b1Irq 2+b2Irq+b3
C=c1Irq 2+c2Irq+c3
2)在油田孔隙度模型和渗透率模型的基础上,结合油柱高度分布状况,代入步骤1)所得到的饱和度公式,利用地质模型建模软件即实现饱和度场模型的构建。
所述步骤2)中所利用的地质模型建模软件为Petrel。
本发明由于采取以上技术方案,其具有以下优点:1、本发明以非均质性储层的实际岩心样品毛管压力资料为基础,综合考虑地质油藏特征,深入分析含油饱和度与毛管压力、储层质量系数的关系,首先建立饱和度经验公式,其次构建饱和度场地质模型,从而实现了饱和度场的合理、定量、精细表征。2、应用本发明方法在某裂谷盆地砂岩油藏进行应用验证,目标油田为构造油藏,扇三角洲沉积,沉积物源近,储层岩性分选差异大,从细砂岩至粗砂岩均有分布,孔隙结构复杂,储层非均质性强,油柱高度大。应用本发明方法,所建立的饱和度场精度显著提高,通过与测井解释资料进行对比分析,符合率高达85%。因此,该方法可以对非均质性强、油柱高度大的油田的饱和度场进行精细构建,对这类油藏饱和度场的三维展布及对后期开发生产具有重要的指导意义。
附图说明
图1是储层质量系数与饱和度系数之间的相关性拟合曲线;其中,图(a)、图(b)、图(c)分别是储层质量系数与饱和度系数A、B、C之间的相关性拟合曲线;
图2是由油田孔隙度模型、渗透率模型以及饱和度公式建立饱和度场模型示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明提出了一种非均质性储层饱和度场精细表征方法,包括以下步骤:
1)建立饱和度公式,具体过程如下:
①根据毛管压力曲线的特点,建立样品饱和度与毛管压力之间的Lambda函数关系,其表达式为:
Sw=APc B+C (1)
式中,Sw表示含水饱和度;Pc表示油藏条件下的毛管压力;A、B、C为饱和度系数。
毛管压力与油柱高度之间具有下述函数关系:
Pc=H*(ρw-ρo)*g (2)
式中,H表示油柱高度;ρw、ρo分别表示水和油的密度;g为重力加速度。
将式(2)代入式(1),得到:
Sw=A[H*(ρw-ρo)*g]B+C (3)
②将储层质量系数Irq作为约束条件,对饱和度公式中的系数A、B、C进行拟合,其中Irq可通过下式计算:
式中,Irq表示储层质量系数;K表示岩石渗透率;表示岩石孔隙度;
储层质量系数Irq能综合反映储层的岩性、物性特征,与饱和度系数A、B、C具有良好的相关性(如图1所示),因此,其相关关系可表示为:
式中,a1、a2、a3、b1、b2、b3、c1、c2、c3表示储层质量系数与饱和度系数A、B、C相关关系式的系数。
③由式(3)和式(5),最终确定饱和度的公式如下:
Sw=A[H*(ρw-ρo)*g]B+C
A=a1Irq 2+a2Irq+a3
B=b1Irq 2+b2Irq+b3
C=c1Irq 2+c2Irq+c3
2)在油田孔隙度模型和渗透率模型的基础上,结合油柱高度分布状况,代入步骤1)所得到的饱和度公式,利用地质模型建模软件(例如Petrel)即可实现饱和度场模型的构建(如图2所示)。
上述各实施例仅用于说明本发明,其中方法的实施步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (2)
1.一种非均质性储层饱和度场精细表征方法,包括以下步骤:
1)建立饱和度公式,具体过程如下:
①建立样品饱和度与毛管压力之间的Lambda函数关系,其表达式为:
Sw=APc B+C (1)
式中,Sw表示含水饱和度;Pc表示油藏条件下的毛管压力;A、B、C为饱和度系数;
毛管压力与油柱高度之间具有下述函数关系:
Pc=H*(ρw-ρo)*g (2)
式中,H表示油柱高度;ρw、ρo分别表示水和油的密度;g为重力加速度。
将式(2)代入式(1),得到:
Sw=A[H*(ρw-ρo)*g]B+C (3)
②将储层质量系数Irq作为约束条件,对饱和度公式中的系数A、B、C进行拟合,其中Irq的计算公式为:
式中,Irq表示储层质量系数;K表示岩石渗透率;表示岩石孔隙度;
储层质量系数Irq与饱和度系数A、B、C具有良好的相关性,其相关关系表示为:
式中,a1、a2、a3、b1、b2、b3、c1、c2、c3表示储层质量系数与饱和度系数A、B、C相关关系式的系数。
③由式(3)和式(5),最终确定饱和度的公式如下:
Sw=A[H*(ρw-ρo)*g]B+C
A=a1Irq 2+a2Irq+a3
B=b1Irq 2+b2Irq+b3
C=c1Irq 2+c2Irq+c3
2)在油田孔隙度模型和渗透率模型的基础上,结合油柱高度分布状况,代入步骤1)所得到的饱和度公式,利用地质模型建模软件即实现饱和度场模型的构建。
2.如权利要求1所示的一种非均质性储层饱和度场精细表征方法,其特征在于:所述步骤2)中所利用的地质模型建模软件为Petrel。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810025475.4A CN108241785A (zh) | 2018-01-11 | 2018-01-11 | 一种非均质性储层饱和度场精细表征方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810025475.4A CN108241785A (zh) | 2018-01-11 | 2018-01-11 | 一种非均质性储层饱和度场精细表征方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108241785A true CN108241785A (zh) | 2018-07-03 |
Family
ID=62699564
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810025475.4A Pending CN108241785A (zh) | 2018-01-11 | 2018-01-11 | 一种非均质性储层饱和度场精细表征方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108241785A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110056346A (zh) * | 2019-04-17 | 2019-07-26 | 中国石油大港油田勘探开发研究院 | 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法 |
CN111561312A (zh) * | 2020-04-21 | 2020-08-21 | 中国石油天然气股份有限公司 | 基于常规测井数据的饱和度计算方法及装置 |
CN113705122A (zh) * | 2020-08-26 | 2021-11-26 | 中国石油大学(北京) | 润湿反转损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103278436A (zh) * | 2013-02-01 | 2013-09-04 | 西安石油大学 | 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法 |
CN107044277A (zh) * | 2017-06-06 | 2017-08-15 | 西南石油大学 | 低渗透非均质油藏水平井重复压裂增产潜力评价方法 |
-
2018
- 2018-01-11 CN CN201810025475.4A patent/CN108241785A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103278436A (zh) * | 2013-02-01 | 2013-09-04 | 西安石油大学 | 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法 |
CN107044277A (zh) * | 2017-06-06 | 2017-08-15 | 西南石油大学 | 低渗透非均质油藏水平井重复压裂增产潜力评价方法 |
Non-Patent Citations (2)
Title |
---|
C.R史密斯 等: "《实用油藏工程》", 30 June 1995, 北京:石油工业出版社 * |
徐伟 等: "乌干达Albert盆地J油田复杂碎屑岩储层含水饱和度定量表征", 《中国海上油气》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110056346A (zh) * | 2019-04-17 | 2019-07-26 | 中国石油大港油田勘探开发研究院 | 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法 |
CN110056346B (zh) * | 2019-04-17 | 2022-05-17 | 中国石油天然气股份有限公司 | 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法 |
CN111561312A (zh) * | 2020-04-21 | 2020-08-21 | 中国石油天然气股份有限公司 | 基于常规测井数据的饱和度计算方法及装置 |
CN111561312B (zh) * | 2020-04-21 | 2023-09-26 | 中国石油天然气股份有限公司 | 基于常规测井数据的饱和度计算方法及装置 |
CN113705122A (zh) * | 2020-08-26 | 2021-11-26 | 中国石油大学(北京) | 润湿反转损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统 |
CN113705122B (zh) * | 2020-08-26 | 2022-09-30 | 中国石油大学(北京) | 润湿反转损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110334431A (zh) | 一种低渗透致密气藏单井控制储量计算及剩余气分析方法 | |
CN104965979B (zh) | 一种致密砂岩有效储层识别方法 | |
CN1912655A (zh) | 一种沉积相控制预测砂岩油藏储层属性的方法 | |
CN110056346B (zh) | 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法 | |
CN105404735B (zh) | 裂缝与基质对超低渗透油藏单井产量贡献率的定量评价方法 | |
CN105386751B (zh) | 一种基于油藏渗流模型的水平井测井产能预测方法 | |
CN111257926B (zh) | 一种利用老旧地震资料进行古河谷铀储层预测的方法 | |
CN108241785A (zh) | 一种非均质性储层饱和度场精细表征方法 | |
CN101929973A (zh) | 裂缝储层含油气饱和度定量计算方法 | |
CN111815773B (zh) | 适用于机器学习算法的三维复杂地质模型标签制作方法 | |
CN107448197B (zh) | 一种海相砂岩油藏窜流通道空间发育定量识别的方法 | |
CN106150477A (zh) | 一种确定单井控制储量的方法 | |
CN103510947A (zh) | 滩坝砂岩微相识别模式的建立及其在微相识别中的应用 | |
CN107766978A (zh) | 不规则井网的智能优化方法 | |
CN114427432B (zh) | 一种气藏剩余气开发潜力确定方法 | |
CN104156617B (zh) | 用于多层砂岩气藏气层品质分类表征的六阶段建模方法 | |
CN114549774A (zh) | 一种基于钻孔数据的三维地层建模方法 | |
CN104459790A (zh) | 含油气性盆地有效储层的分析方法和装置 | |
CN110863818A (zh) | 一种剩余油/气分布的描述方法及装置 | |
CN106501146A (zh) | 一种致密油藏物性上限确定方法及装置 | |
CN108932397A (zh) | 一种双重介质油藏数值模拟调参方法 | |
Qinglong | Quantitative characterization and development countermeasures of multi-layer heterogeneous interlayer channeling flow | |
Leemput et al. | Full-field reservoir modeling of Central Oman gas-condensate fields | |
CN111550238B (zh) | 一种等效j函数约束的饱和度建模方法 | |
CN115700321A (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: 20180703 |
|
RJ01 | Rejection of invention patent application after publication |