CN113029892B - 基于区域统计规律的油水相对渗透率曲线合理性评价方法 - Google Patents
基于区域统计规律的油水相对渗透率曲线合理性评价方法 Download PDFInfo
- Publication number
- CN113029892B CN113029892B CN202010186607.9A CN202010186607A CN113029892B CN 113029892 B CN113029892 B CN 113029892B CN 202010186607 A CN202010186607 A CN 202010186607A CN 113029892 B CN113029892 B CN 113029892B
- Authority
- CN
- China
- Prior art keywords
- oil
- relative permeability
- water
- saturation
- evaluation
- 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.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 288
- 230000035699 permeability Effects 0.000 title claims abstract description 267
- 238000000034 method Methods 0.000 title claims abstract description 105
- 238000011156 evaluation Methods 0.000 claims abstract description 159
- 230000008569 process Effects 0.000 claims abstract description 18
- 239000012071 phase Substances 0.000 claims description 143
- 239000011148 porous material Substances 0.000 claims description 17
- 239000011435 rock Substances 0.000 claims description 16
- 238000005259 measurement Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 239000008346 aqueous phase Substances 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 238000013278 delphi method Methods 0.000 claims description 3
- 238000012847 principal component analysis method Methods 0.000 claims description 3
- 238000011160 research Methods 0.000 claims description 2
- 238000011161 development Methods 0.000 description 10
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000011158 quantitative evaluation Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012854 evaluation process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
- G01N13/04—Investigating osmotic effects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Dispersion Chemistry (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
基于区域统计规律的油水相对渗透率曲线合理性评价方法,包括:确定表征油水相对渗透曲线六个评价参数;利用层次分析法和熵权法组合确定六个评价参数权重;结合区域已有油水相对渗透曲线分别建立六个油水相对渗透评价参数与孔隙度和渗透率函数关系,选取拟合最好的函数关系对油水相渗曲线评价参数预测;计算区域统计规律与实测油水相对渗透曲线六个评价参数误差;结合待评价对象油水相对渗透曲线六个评价参数误差和权重计算综合误差指数;确定待评价对象油水相对渗透曲线评价等级;该方法在开展油水相对渗透曲线合理性评价时综合考虑表征油水相对渗透曲线评价参数、区域统计特征及每个评价参数权重,评价结果受主观因素影响较小,评价结果可靠。
Description
技术领域
本发明属于油气田勘探开发技术领域,尤其涉及一种基于区域统计规律的油水相对渗透率曲线合理性评价方法。
背景技术
油水相对渗透率曲线能够反映油田开发过程中油和水地下渗流的复杂过程,是油田开发生产过程中开展油田开发动态分析、计算油井产能以及开发指标预测中重要的基础资料,确定合理的油田油水相渗曲线对油田高效开发十分重要。海上油田开发受到开发成本的限制,实测油水相对渗透率曲线数量有限,甚至有些油田和区块没有实测的油水相对渗透率曲线,油田在开发生产过程中,通常采用类比法确定油田油水相对渗透率曲线的合理性,该方法属于定性评价,类比过程中不同技术人员应用指标以及标准不一致,影响油田技术人员的开发决策。在这种背景下,开展油水相对渗透率曲线合理性评价具有重要而现实的意义:一方面能够对已有的油水相对渗透率曲线质量进行分析,另一方面能够为无实测油水相对渗透率曲线的油田选取合适的油水相对渗透率曲线提供依据。因此急需建立一种针对油水相对渗透率曲线进行定量评价的方法,指导油田确定合理的油水相对渗透率曲线。为此,申请人设计基于区域统计规律的油水相对渗透率曲线合理性评价方法,解决以上技术问题。
发明内容
本发明目的在于提供一种基于区域统计规律的油水相对渗透率曲线合理性评价方法,以解决针对油水相对渗透率曲线进行定量评价的方法,指导油田确定合理的油水相对渗透率曲线的技术问题。
为实现上述目的,本发明的基于区域统计规律的油水相对渗透率曲线合理性评价方法的具体技术方案如下:
一种基于区域统计规律的油水相对渗透率曲线合理性评价方法,其包括以下步骤:
第一步,确定表征油水相对渗透曲线的六个评价参数,该六个评价参数分别为束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数;
第二步,利用层次分析法和熵权法组合确定束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数六个评价参数的权重;
第三步,开展油水相对渗透率曲线区域统计规律研究,结合区域已有油水相对渗透曲线分别建立六个油水相对渗透评价参数与孔隙度和渗透率的函数关系,选取拟合最好的函数关系开展油水相渗曲线评价参数预测,根据岩石物性参数确定待评价对象油水相对渗透率曲线的理想评价参数;
第四步,计算区域统计规律与实测油水相对渗透曲线六个评价参数误差;
第五步,结合待评价对象油水相对渗透曲线六个评价参数误差和权重计算综合误差指数;
第六步,确定待评价对象油水相对渗透曲线的评价等级。
前述的基于区域统计规律的油水相对渗透率曲线合理性评价方法,其中,
所述第一步是根据霍纳普油水相对渗透率曲线公式确定能够表征油水相对渗透率曲线的评价参数,包括:
油水相对渗透率曲线反映的是油相和水相相对渗透率随含水饱和度变化的关系曲线,霍纳普采用指数形式来表征油水相对渗透率曲线:
式中,Kro——油相相对渗透率,f;Kromax——束缚水饱和度下油相相对渗透率,f;Sw——油相相对渗透率,f;Sor——残余油饱和度,f;Swi——束缚水饱和度,f;m——油相指数,f;Krw——水相相对渗透率,f;Krwmax——残余油饱和度下水相相对渗透率,f;n——水相指数,f;
由霍纳普公式可以看出,油水相对渗透率曲线形态主要受Kromax、Sor、Swi、m、Krwmax以及n等六个参数控制,因此采用这六个参数作为评价参数。
前述的基于区域统计规律的油水相对渗透率曲线合理性评价方法,其中,
所述第二步,结合层次分析法和熵权法组合确定束缚水饱和度下油水相对渗透率曲线六个评价参数的权重,包括:在对具有多个评价参数的评价对象进行综合评价时,首先需要确定每个评价参数的权重;确定权重的方法为主观赋权法、客观赋权法和组合集成赋权法;该主观赋权法根据决策者的经验对每个评价参数的权重赋值,其包括层次分析法、德尔菲法以及因素成对比较法;该客观赋权法根据数理统计方法对各个评价参数值客观数据进行分析确定权重,包括主成分分析法、均方差权重法以及熵权法;合集成赋权法将将主、客观赋权法所得的权重值用集成的方法形成最终的评价参数权重值,其选用层次分析法与熵权法相结合的方法确定评价参数权重,计算综合权重值时选用加权系数进行控制,每个评价参数的权重定义如下:
w1=ε×α1+(1-ε)×β1 (3)
w2=ε×α2+(1-ε)×β2 (4)
w3=ε×α3+(1-ε)×β3 (5)
w4=ε×α4+(1-ε)×β4 (6)
w5=ε×α5+(1-ε)×β5 (7)
w6=ε×α6+(1-ε)×β6 (8)其中,ε为加权系数,加权系数为0.5;w1、w2、w3、w4、w5和w6分别为束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的综合权重;α1、α2、α3、α4、α5和α6分别为利用层次分析法确定的束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的权重;β1、β2、β3、β4、β5和β6分别为利用熵权法确定的束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的权重。
前述的基于区域统计规律的油水相对渗透率曲线合理性评价方法,其中,
所述第三步结合区域已有油水相对渗透率曲线建立六个油水相对渗透评价参数与孔隙度和渗透率的函数关系,包括:油水相对渗透率曲线受到岩石孔隙特征的影响,孔隙特征的度量指标包括孔隙度和渗透率等参数,在建立区域统计规律时,不同区域油水相对渗透率曲线评价参数与孔隙特征的度量参数关系不同;分别建立每个评价参数与孔隙度和渗透率的关系,该关系包括指数关系、线性关系、对数关系和幂指数关系,选取拟合最好的函数关系开展油水相对渗透率曲线评价参数预测,根据岩石物性参数确定待评价对象油水相对渗透率曲线的理想评价参数。
前述的基于区域统计规律的油水相对渗透率曲线合理性评价方法,其中,
所述第四步,计算区域统计规律与实测油水相对渗透曲线六个评价参数的误差是,计算出待评价对象油水相对渗透率曲线的束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数等六个评价参数,计算待评价对象油水相对渗透率曲线的六个评价参数的误差,计算公式如下:
error1、error2、error3、error4、error5和error6分别为待评价油水相对渗透率曲线束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的计算误差;
Kromax″、Sor″、Swi″、m″、Krwmax″和n″分别为根据步骤三区域统计规律确定的待评价油水相对渗透率曲线束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的计算值;
Kromax′、Sor′、Swi′、m′、Krwmax′和n′分别为待评价油水相对渗透率曲线束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的实际值。
前述的基于区域统计规律的油水相对渗透率曲线合理性评价方法,其中,
所述第五步,结合待评价对象油水相对渗透曲线六个评价参数误差和权重计算综合误差指数是,根据步骤二确定的六个评价参数的权重和步骤三确定的六个评价参数的误差,计算综合误差指数,具体包括:利用公式计算待评价相渗曲线的综合误差指数:
式中A为待评价对象油水相对渗透率曲线的综合误差指数;
前述的基于区域统计规律的油水相对渗透率曲线合理性评价方法,其中,
所述第六步,根据综合误差指数确定评价结果,当A≤1时认为待评价对象油水相对渗透率曲线的评价结果为好,当1<A≤1.5时认为待评价对象油水相对渗透率曲线的评价结果为较好,当1.5<A≤2时认为待评价对象油水相对渗透率曲线的评价结果为一般,当A>2时认为待评价对象油水相对渗透率曲线的评价结果为较差。
本发明的基于区域统计规律的油水相对渗透率曲线合理性评价方法具有以下优点:本方法在开展油水相对渗透曲线合理性评价时综合考虑了表征油水相对渗透曲线的评价参数、区域的统计特征以及每个评价参数的权重,评价结果受主观因素的影响较小,评价结果可靠。根据区域统计规律结合待评价对象的岩石物性参数确定油水相对渗透率曲线的理想形态,对比理想形态与实际实验结果的差异,确定油水相对渗透曲线的合理性。
附图说明
图1为本发明基于区域统计规律的油水相对渗透率曲线合理性评价方法的流程图。
图2为本发明六个评价参数的评价权重图。
图3为本发明海上油田某区块束缚水饱和度下油相相对渗透率与孔隙度关系曲线图。
图4为本发明海上油田某区块束缚水饱和度下油相相对渗透率与空气渗透率关系曲线图。
图5为本发明海上油田某区块残余油饱和度下水相相对渗透率与孔隙度关系曲线图。
图6为本发明海上油田某区块残余油饱和度下水相相对渗透率与空气渗透率关系曲线图。
图7为本发明海上油田某区块残余油饱和度与孔隙度关系曲线图。
图8为本发明海上油田某区块残余油饱和度与空气渗透率关系曲线图。
图9为本发明海上油田某区块束缚水饱和度与孔隙度关系曲线图。
图10为本发明海上油田某区块束缚水饱和度与空气渗透率关系曲线图。
图11为本发明海上油田某区块油相指数与孔隙度关系曲线图。
图12为本发明海上油田某区块油相指数与空气渗透率关系曲线图。
图13为本发明海上油田某区块水相指数与孔隙度关系曲线图。
图14为本发明海上油田某区块水相指数与空气渗透率关系曲线图。
具体实施方式
为更好地了解本发明的目的、结构及功能,下面结合附图,对本发明基于区域统计规律的油水相对渗透率曲线合理性评价方法做进一步详细描述。
本发明的基于区域统计规律的油水相对渗透率曲线合理性评价方法的具体实施流程,如图1所示:
第一步,确定表征油水相对渗透曲线的六个评价参数,该六个评价参数分别为束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数;
根据霍纳普油水相对渗透率曲线公式确定能够表征油水相对渗透率曲线的评价参数,包括:
油水相对渗透率曲线反映的是油相和水相相对渗透率随含水饱和度变化的关系曲线,霍纳普采用指数形式来表征油水相对渗透率曲线:
式中,Kro——油相相对渗透率,f;Kromax——束缚水饱和度下油相相对渗透率,f;Sw——油相相对渗透率,f;Sor——残余油饱和度,f;Swi——束缚水饱和度,f;m——油相指数,f;Krw——水相相对渗透率,f;Krwmax——残余油饱和度下水相相对渗透率,f;n——水相指数,f;
由霍纳普公式可以看出,油水相对渗透率曲线形态主要受Kromax、Sor、Swi、m、Krwmax以及n等六个参数控制,因此采用这六个参数作为评价参数。
第二步,利用层次分析法和熵权法组合确定束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数六个评价参数的权重;
结合层次分析法和熵权法组合确定束缚水饱和度下油水相对渗透率曲线六个评价参数的权重,包括:在对具有多个评价参数的评价对象进行综合评价时,首先需要确定每个评价参数的权重;确定权重的方法为主观赋权法、客观赋权法和组合集成赋权法;该主观赋权法根据决策者的经验对每个评价参数的权重赋值,其包括层次分析法、德尔菲法以及因素成对比较法;该客观赋权法根据数理统计方法对各个评价参数值客观数据进行分析确定权重,包括主成分分析法、均方差权重法以及熵权法;合集成赋权法将将主、客观赋权法所得的权重值用集成的方法形成最终的评价参数权重值,其选用层次分析法与熵权法相结合的方法确定评价参数权重,计算综合权重值时选用加权系数进行控制,每个评价参数的权重定义如下:
w1=ε×α1+(1-ε)×β1 (3)
w2=ε×α2+(1-ε)×β2 (4)
w3=ε×α3+(1-ε)×β3(5)
w4=ε×α4+(1-ε)×β4 (6)
w5=ε×α5+(1-ε)×β5 (7)
w6=ε×α6+(1-ε)×β6 (8)
其中,ε为加权系数,加权系数为0.5;w1、w2、w3、w4、w5和w6分别为束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的综合权重;α1、α2、α3、α4、α5和α6分别为利用层次分析法确定的束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的权重;β1、β2、β3、β4、β5和β6分别为利用熵权法确定的束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的权重;
本实施例选取具有一定开发生产经验和实验经验的专家对六个评价参数之间的重要性进行打分,构造判断(成对比较)矩阵,根据判断(成对比较)矩阵确定主观权重,六个评价参数主观权重如图2所示。
本实施例选取待分析的海上油田区块,筛选出具有代表性的油水相对渗透率曲线50条,统计每条相渗曲线的六个评价参数和岩石孔隙特征参数,岩石孔隙特征参数包括孔隙度和渗透率数据,绘制成表,采用熵权法计算客观权重,结合主观权重和客观权重计算得到综合权重,如图2所示。
第三步,结合区域已有油水相对渗透曲线分别建立六个油水相对渗透评价参数与孔隙度和渗透率的函数关系,选取拟合最好的函数关系开展油水相渗曲线评价参数预测;
结合区域已有油水相对渗透率曲线建立六个油水相对渗透评价参数与孔隙度和渗透率的函数之间的关系,包括:油水相对渗透率曲线受到岩石孔隙特征的影响,孔隙特征的度量指标包括孔隙度和渗透率等参数,在建立区域统计规律时,不同区域油水相对渗透率曲线评价参数与孔隙特征的度量参数关系不同;分别建立每个评价参数与孔隙度和渗透率的关系,该关系包括指数关系、线性关系、对数关系和幂指数关系,选取拟合最好的函数关系开展油水相对渗透率曲线评价参数预测;
本实施例利用统计的油水相对渗透率曲线的六个评价参数和岩石孔隙特征参数数据,选取48个样本分别建立六个评价参数与孔隙度和渗透率之间的函数关系,建立的函数关系包括指数关系、线性关系、对数关系和乘幂关系四种,海上油田某区块六个评价参数与岩石孔隙特征参数之间的回归关系如图3至图14所示。
以R2最大为原则,确定六个评价参数与岩石孔隙特征参数选取拟合最好的函数关系。选取结果如下:束缚水饱和度下油相相对渗透率选取空气渗透率建立乘幂关系、残余油饱和度下水相相对渗透率空气渗透率建立乘幂关系、残余油饱和度选取空气渗透率建立对数关系、束缚水饱和度选取空气渗透率建立乘幂关系、油相指数选取空气渗透率建立乘幂关系、水相指数选取空气渗透率建立乘幂关系。具体表达式为:
Kromax″=0.0646×Kair 0.2899 (16)
Krwmax″=0.1138×Kair 0.1855 (17)
Sor=0.0207Ln(Kair)+0.1801 (18)
Swi″=0.8736×Kair -0.174 (19)
m″=13.676×Kair -0.255 (20)
n″=0.9124×Kair 0.2409 (21)
式中,Kair——空气渗透率,10-3μm2。
第四步,计算区域统计规律与实测油水相对渗透曲线六个评价参数误差;计算区域统计规律与实测油水相对渗透曲线六个评价参数的误差是,计算出待评价对象油水相对渗透率曲线的束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数等六个评价参数,计算待评价对象油水相对渗透率曲线的六个评价参数的误差,计算公式如下:
error1、error2、error3、error4、error5和error6分别为待评价油水相对渗透率曲线束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的计算误差;
Kromax″、Sor″、Swi″、m″、Krwmax″和n″分别为根据步骤三区域统计规律确定的待评价油水相对渗透率曲线束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的计算值;
Kromax′、Sor′、Swi′、m′、Krwmax′和n′分别为待评价油水相对渗透率曲线束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的实际值;
计算得到样本49的六个评价参数的误差分别为103%、86.0%、94.0%、54.0%、74.0%、69.0%,而样本50的六个评价参数的误差分别为233.0%、323.0%、286.0%、201.0%、326.0%、342.0%。
第五步,结合待评价对象油水相对渗透曲线六个评价参数误差和权重计算综合误差指数;
结合待评价对象油水相对渗透曲线六个评价参数误差和权重计算综合误差指数是,根据步骤二确定的六个评价参数的权重和步骤三确定的六个评价参数的误差,计算综合误差指数,具体包括:利用公式计算待评价相渗曲线的综合误差指数:
式中A为待评价对象油水相对渗透率曲线的综合误差指数;
利用确定的六个评价参数与岩石孔隙特征参数关系对待评价2个样本评价参数进行预测,计算预测结果与实验测得结果之间的误差,结合每个评价参数的权重得到待评价对象油水相对渗透率曲线的综合误差指数。计算得到样本49的综合误差指数为0.8,而样本50的综合误差指数为2.8。
第六步,综合误差指数确定评价结果,当A≤1时认为待评价对象油水相对渗透率曲线的评价结果为好,当1<A≤1.5时认为待评价对象油水相对渗透率曲线的评价结果为较好,当1.5<A≤2时认为待评价对象油水相对渗透率曲线的评价结果为一般,当A>2时认为待评价对象油水相对渗透率曲线的评价结果为较差。
在本实施例中,综合评价认为样本49的评价结果为好,该曲线被应用于数值模拟中,样本50的评价结果较差,在进行该区域数值模拟研究时,未考虑使用该曲线。
本发明核心发明点是:1、根据霍纳普公式确定了定量评价油水相对渗透率曲线形态的评价参数。采用主观和客观方法确定每个评价参数的权重。2、基于区域的已有样本点数据建立了孔隙度、渗透率与评价参数之间的关系。选择拟合程度最好的曲线对评价对象的合理曲线形态进行预测。对比合理曲线与评价对象的差异,计算综合误差指数,从而实现定量评价的目的。
根据区域统计规律结合待评价对象的岩石物性参数确定油水相对渗透率曲线的理想形态,对比理想形态与实际实验结果的差异,确定油水相对渗透曲线的合理性;目前实际应用时主要首先选取类比油田的相渗曲线,对比实验相渗曲线与类比油田相渗曲线形态差异,差异较小时就认为是合理的,属于定性评价,本发明目的是将这个评价过程定量化;该确定油水相对渗透曲线的合理性是根据霍纳普油水相对渗透率曲线公式确定能够表征油水相对渗透曲线的评价参数,结合主观和客观方法确定每个评价参数的权重,利用区域统计规律计算出每个评价参数的相对误差,根据权重和相对误差计算得出的综合误差指数,确定待评价对象油水相对渗透曲线的合理性评价结果。
可以理解,本发明是通过一些实施例进行描述的,本领域技术人员知悉的,在不脱离本发明的精神和范围的情况下,可以对这些特征和实施例进行各种改变或等效替换。另外,在本发明的教导下,可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此,本发明不受此处所公开的具体实施例的限制,所有落入本申请的权利要求范围内的实施例都属于本发明所保护的范围内。
Claims (6)
1.一种基于区域统计规律的油水相对渗透率曲线合理性评价方法,其特征在于,包括以下步骤:
第一步,确定表征油水相对渗透率曲线的六个评价参数,该六个评价参数分别为束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数;
第二步,利用层次分析法和熵权法组合确定束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数六个评价参数的权重;
第三步,开展油水相对渗透率曲线区域统计规律研究,结合区域已有油水相对渗透率曲线分别建立六个油水相对渗透率曲线评价参数与孔隙度和空气渗透率的函数关系,选取拟合最好的函数关系开展油水相渗透率曲线评价参数预测,根据岩石孔隙特征的度量指标确定待评价对象油水相对渗透率曲线的理想评价参数;
第四步,计算区域统计规律与实测油水相对渗透率曲线六个评价参数误差;
第五步,结合待评价对象油水相对渗透率曲线六个评价参数误差和权重计算综合误差指数;
第六步,确定待评价对象油水相对渗透率曲线的评价结果;
所述第四步,计算区域统计规律与实测油水相对渗透率曲线六个评价参数的误差是,计算出待评价对象油水相对渗透率曲线的束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数六个评价参数,计算待评价对象油水相对渗透率曲线的六个评价参数的误差,计算公式如下:
error1、error2、error3、error4、error5和error6分别为待评价对象油水相对渗透率曲线束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的计算误差;
Kromax”、Sor”、Swi”、m”、Krwmax”和n”分别为根据步骤三区域统计规律确定的待评价对象油水相对渗透率曲线束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的计算值;
Kromax'、Sor'、Swi'、m'、Krwmax'和n'分别为待评价对象油水相对渗透率曲线束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的实际值。
2.根据权利要求1所述的基于区域统计规律的油水相对渗透率曲线合理性评价方法,其特征在于,
所述第一步是根据霍纳普油水相对渗透率曲线公式确定能够表征油水相对渗透率曲线的评价参数,包括:
油水相对渗透率曲线反映的是油相和水相相对渗透率随含水饱和度变化的关系曲线,霍纳普采用指数形式来表征油水相对渗透率曲线:
式中,Kro——油相相对渗透率;Kromax——束缚水饱和度下油相相对渗透率;Sw——含水饱和度;Sor——残余油饱和度;Swi——束缚水饱和度;m——油相指数;Krw——水相相对渗透率;Krwmax——残余油饱和度下水相相对渗透率;n——水相指数;
由霍纳普油水相对渗透率曲线公式可以看出,油水相对渗透率曲线形态主要受Kromax、Sor、Swi、m、Krwmax以及n六个参数控制,因此采用这六个参数作为评价参数。
3.根据权利要求1所述的基于区域统计规律的油水相对渗透率曲线合理性评价方法,其特征在于,
所述第二步,结合层次分析法和熵权法组合确定束缚水饱和度下油水相对渗透率曲线六个评价参数的权重,包括:在对具有多个评价参数的待评价对象进行综合评价时,首先需要确定每个评价参数的权重;确定权重的方法为主观赋权法、客观赋权法和组合集成赋权法;该主观赋权法根据决策者的经验对每个评价参数的权重赋值,其包括层次分析法、德尔菲法以及因素成对比较法;该客观赋权法根据数理统计方法对各个评价参数值客观数据进行分析确定权重,包括主成分分析法、均方差权重法以及熵权法;组合集成赋权法将主、客观赋权法所得的权重值用集成的方法形成最终的评价参数权重值,其选用层次分析法与熵权法相结合的方法确定评价参数权重,计算综合权重值时选用加权系数进行控制,每个评价参数的权重定义如下:
w1=ε×α1+(1-ε)×β1 (3)
w2=ε×α2+(1-ε)×β2 (4)
w3=ε×α3+(1-ε)×β3 (5)
w4=ε×α4+(1-ε)×β4 (6)
w5=ε×α5+(1-ε)×β5 (7)
w6=ε×α6+(1-ε)×β6 (8)
其中,ε为加权系数,加权系数为0.5;w1、w2、w3、w4、w5和w6分别为束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的综合权重;α1、α2、α3、α4、α5和α6分别为利用层次分析法确定的束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的权重;β1、β2、β3、β4、β5和β6分别为利用熵权法确定的束缚水饱和度下油相相对渗透率、残余油饱和度下水相相对渗透率、残余油饱和度、束缚水饱和度、油相指数和水相指数的权重。
4.根据权利要求1所述的基于区域统计规律的油水相对渗透率曲线合理性评价方法,其特征在于,
所述第三步结合区域已有油水相对渗透率曲线建立六个油水相对渗透率曲线评价参数与孔隙度和空气渗透率的函数关系,包括:油水相对渗透率曲线受到岩石孔隙特征的影响,岩石孔隙特征的度量指标包括孔隙度和空气渗透率参数,在建立区域统计规律时,不同区域油水相对渗透率曲线评价参数与孔隙特征的度量参数关系不同;分别建立每个评价参数与孔隙度和空气渗透率的关系,该关系包括指数关系、线性关系、对数关系和幂指数关系,选取拟合最好的函数关系开展油水相对渗透率曲线评价参数预测,根据岩石孔隙特征的度量指标确定待评价对象油水相对渗透率曲线的理想评价参数。
6.根据权利要求1所述的基于区域统计规律的油水相对渗透率曲线合理性评价方法,其特征在于,
所述第六步,根据综合误差指数确定评价结果,当A≤1时认为待评价对象油水相对渗透率曲线的评价结果为好,当1<A≤1.5时认为待评价对象油水相对渗透率曲线的评价结果为较好,当1.5<A≤2时认为待评价对象油水相对渗透率曲线的评价结果为一般,当A>2时认为待评价对象油水相对渗透率曲线的评价结果为较差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010186607.9A CN113029892B (zh) | 2020-03-17 | 2020-03-17 | 基于区域统计规律的油水相对渗透率曲线合理性评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010186607.9A CN113029892B (zh) | 2020-03-17 | 2020-03-17 | 基于区域统计规律的油水相对渗透率曲线合理性评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113029892A CN113029892A (zh) | 2021-06-25 |
CN113029892B true CN113029892B (zh) | 2022-12-13 |
Family
ID=76458570
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010186607.9A Active CN113029892B (zh) | 2020-03-17 | 2020-03-17 | 基于区域统计规律的油水相对渗透率曲线合理性评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113029892B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116072232B (zh) * | 2021-12-29 | 2024-03-19 | 中国石油天然气集团有限公司 | 一种相对渗透率曲线确定方法、装置、设备和存储介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102468997A (zh) * | 2011-07-01 | 2012-05-23 | 中国人民解放军国防科学技术大学 | 一种多维网络安全指标体系稳定性评估方法 |
CN102915406A (zh) * | 2012-10-26 | 2013-02-06 | 中国石油大学(华东) | 径向流条件下油水相对渗透率曲线的计算方法 |
CN105740617A (zh) * | 2016-01-28 | 2016-07-06 | 中国电子科技集团公司第十研究所 | 一体化天馈系统品质度的衡量方法 |
CN108593514A (zh) * | 2018-03-26 | 2018-09-28 | 中国石油化工股份有限公司 | 基于储层物性的油水相对渗透率表征处理方法 |
CN109356564A (zh) * | 2018-11-08 | 2019-02-19 | 中国海洋石油集团有限公司 | 一种基于三元动态参数的水驱油相渗时变规律表征方法 |
CN109653725A (zh) * | 2018-09-13 | 2019-04-19 | 山东鼎维石油科技有限公司 | 基于沉积微相和岩石相的混积储层水淹程度测井解释方法 |
CN109800521A (zh) * | 2019-01-28 | 2019-05-24 | 中国石油大学(华东) | 一种基于机器学习的油水相对渗透率曲线计算方法 |
CN110288258A (zh) * | 2019-07-02 | 2019-09-27 | 中国石油化工股份有限公司 | 一种高含水油藏剩余油挖潜方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1744144A1 (de) * | 2005-07-12 | 2007-01-17 | Martechnic GmbH | Verfahren und Einrichtung zur Messung der Belastung einer Atmosphäre mit Aerosolen und Stäuben |
KR101209857B1 (ko) * | 2009-02-20 | 2012-12-10 | 삼성코닝정밀소재 주식회사 | 유리 표면 이물 검사 장치 및 방법 |
CN102358069B (zh) * | 2011-08-09 | 2013-09-25 | 成立 | 遮光玻璃板 |
CN105547948A (zh) * | 2016-01-27 | 2016-05-04 | 重庆川仪分析仪器有限公司 | 粉尘浓度在线监测用自动校准装置 |
CN105486620B (zh) * | 2016-02-02 | 2018-08-31 | 广东美的制冷设备有限公司 | 一种空气净化器、激光粉尘浓度传感器及其检测方法 |
CN106093350B (zh) * | 2016-08-10 | 2018-05-04 | 中国石油天然气股份有限公司 | 确定非均质碳酸盐岩储层饱和度指数的方法 |
-
2020
- 2020-03-17 CN CN202010186607.9A patent/CN113029892B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102468997A (zh) * | 2011-07-01 | 2012-05-23 | 中国人民解放军国防科学技术大学 | 一种多维网络安全指标体系稳定性评估方法 |
CN102915406A (zh) * | 2012-10-26 | 2013-02-06 | 中国石油大学(华东) | 径向流条件下油水相对渗透率曲线的计算方法 |
CN105740617A (zh) * | 2016-01-28 | 2016-07-06 | 中国电子科技集团公司第十研究所 | 一体化天馈系统品质度的衡量方法 |
CN108593514A (zh) * | 2018-03-26 | 2018-09-28 | 中国石油化工股份有限公司 | 基于储层物性的油水相对渗透率表征处理方法 |
CN109653725A (zh) * | 2018-09-13 | 2019-04-19 | 山东鼎维石油科技有限公司 | 基于沉积微相和岩石相的混积储层水淹程度测井解释方法 |
CN109356564A (zh) * | 2018-11-08 | 2019-02-19 | 中国海洋石油集团有限公司 | 一种基于三元动态参数的水驱油相渗时变规律表征方法 |
CN109800521A (zh) * | 2019-01-28 | 2019-05-24 | 中国石油大学(华东) | 一种基于机器学习的油水相对渗透率曲线计算方法 |
CN110288258A (zh) * | 2019-07-02 | 2019-09-27 | 中国石油化工股份有限公司 | 一种高含水油藏剩余油挖潜方法 |
Non-Patent Citations (2)
Title |
---|
油藏综合相渗曲线拟合方法;戴胜群等;《油气藏评价与开发》;20110626(第03期);全文 * |
相对渗透率曲线标准化方法评价;缪飞飞等;《断块油气田》;20131219(第06期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113029892A (zh) | 2021-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109838230B (zh) | 油藏水淹层的定量评价方法 | |
CN107038516B (zh) | 一种中渗复杂断块油藏水驱开发效果定量评价方法 | |
CN111425193B (zh) | 一种基于聚类分析测井岩石物理相划分的储层可压性评价方法 | |
CN113190967B (zh) | 一种基于多目标筛选的细支卷烟材料组配方法 | |
CN105630743A (zh) | 一种光谱波数的选择方法 | |
CN110685651B (zh) | 一种多层合采气井产量劈分方法及系统 | |
RU2709047C1 (ru) | Способ адаптации гидродинамической модели продуктивного пласта нефтегазоконденсатного месторождения с учетом неопределенности геологического строения | |
CN112070356A (zh) | 一种基于rf-lssvm模型预测混凝土抗碳化性能的方法 | |
CN113029892B (zh) | 基于区域统计规律的油水相对渗透率曲线合理性评价方法 | |
CN114048973A (zh) | 一种基于组合赋权的河流水质等级评价方法 | |
CN112253087A (zh) | 一种基于多源测井资料的生物扰动储集层物性计算方法 | |
CN111339719A (zh) | 一种缝洞型油藏注水数值模拟自动历史拟合的方法 | |
CN111476448B (zh) | 一种多层合采直井动态产量劈分方法 | |
CN114580261A (zh) | 可采储量预测模型建立方法、可采储量预测方法和装置 | |
CN107133879B (zh) | 一种筛选相似油田的方法 | |
CN110486008B (zh) | 一种径向复合油藏的参数解释方法及系统 | |
CN112983377A (zh) | 低渗透油藏直井重复压裂增产改造潜力评价及选井方法 | |
CN109655394B (zh) | 一种孔喉结构参数约束下的核磁共振t2谱渗透率计算方法 | |
CN116950654A (zh) | 一种低渗致密砂岩气藏开发效果评价方法 | |
CN111488666A (zh) | 一种气藏水平井产能预测模型建立、预测方法及装置 | |
CN111441766B (zh) | 一种分层试气无阻流量预测方法 | |
CN111155980A (zh) | 一种水流优势通道识别方法及装置 | |
CN115705452A (zh) | 整装砂岩油藏开发中后期的新型采收率预测方法 | |
CN112459775A (zh) | 一种剩余油分散程度表征方法 | |
CN108572129B (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 |