CN110056346B - 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法 - Google Patents

一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法 Download PDF

Info

Publication number
CN110056346B
CN110056346B CN201910310522.4A CN201910310522A CN110056346B CN 110056346 B CN110056346 B CN 110056346B CN 201910310522 A CN201910310522 A CN 201910310522A CN 110056346 B CN110056346 B CN 110056346B
Authority
CN
China
Prior art keywords
oil
oil reservoir
well
original
dimensional
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
Application number
CN201910310522.4A
Other languages
English (en)
Other versions
CN110056346A (zh
Inventor
江艳平
窦松江
李涛
张东星
彭飞飞
李伟男
邢卫东
乐涛涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201910310522.4A priority Critical patent/CN110056346B/zh
Publication of CN110056346A publication Critical patent/CN110056346A/zh
Application granted granted Critical
Publication of CN110056346B publication Critical patent/CN110056346B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects

Abstract

本发明公开了一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法,包括基础数据收集整理、确定油藏油水界面海拔深度及油柱高度、确定参与油藏原始含水饱和度分析和模拟的井位、构建原始含水饱和度随油柱高度的趋势变化函数、建立油藏三维含水饱和度变化趋势体和以井数据为硬数据、三维含水饱和度趋势体为趋势约束,建立油藏三维原始含水饱和度模型。本发明目的在于通过构建饱和度趋势变化函数,合理表征油藏内原始含水饱和度随油柱高度趋势变化的地质规律,从而提高原始含水饱和度井间预测精度,实现高精度原始含水饱和度模型的建立,满足油藏开发阶段油藏精细研究和数值模拟研究的需要。

Description

一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法
技术领域
本发明涉及一种油藏三维原始含水饱和度,具体为一种基于趋势变化函 数的油藏三维原始含水饱和度模拟方法,属于石油开发技术领域。
背景技术
在油田开发阶段油藏描述研究工作中,需要通过油藏数值模拟来分析当 前并预测未来不同时间段的含水饱和度(Sw)变化及剩余油的分布情况,因 此需要为油藏数值模拟研究提供高精度的油藏三维原始Sw模型。特别是在油 田开发中后期,剩余油分布高度复杂,对油藏研究越来越精细,因此对油藏 原始Sw模型的精度要求也越来越高。
目前针对油藏三维原始Sw建模方法还很少,主要是建立常数均质Sw模 型,即设定油水界面上下分别为一常数值,由于油藏内储层为非均质性变化, 其原始Sw受储层孔隙结构影响较大,且由于油水因为密度差异在垂向上存在 重力分异作用,因此油藏内原始Sw也呈现非均质性变化,显然常数模型不能 体现油藏内Sw的非均质性分布,其精度远远不能满足油藏精细开发的需求; 其次,少数是利用经验公式法来模拟原始Sw,由于不同类型油藏地质成因不 同,储层及流体性质差别也较大,用统一的经验公式不足以表征不同类型油藏Sw分布的个性特征,因此其精度也不高。
发明内容
本发明的目的就在于为了解决油田开发阶段,特别是油田开发中后期, 油藏精细研究和数值模拟研究而提供一种基于趋势变化函数的油藏三维原始 含水饱和度模拟方法。
本发明通过以下技术方案来实现上述目的:一种基于趋势变化函数的油 藏三维原始含水饱和度模拟方法,包括以下步骤:
1)基础数据收集整理和分析;
2)确定油藏油水界面深度及油柱高度;
3)确定参与油藏原始含水饱和度分析和模拟的井位;
4)构建原始含水饱和度随油柱高度变化的趋势变化函数;
5)建立油藏三维含水饱和度变化趋势体;
6)以井数据为硬数据、三维含水饱和度趋势体为趋势约束,建立油藏三 维原始含水饱和度模型。。
优选的,所述步骤2)中,确定油柱高度及油藏油水界面深度;包括以下 几个方面:
①根据油藏剖面图中油水界面所在的海拔深度,确定为油水界面深度 (OWC);
②根据油藏含油面积图中,含油范围最高部位和最低部位深度差,确定 为油柱高度(H)。
优选的,所述步骤3)中,确定参与油藏原始含水饱和度分析和模拟的井 位,包括以下几个方面:
①油藏内探井、评价井数较多,分布位置具原始Sw区域代表性,则选取 油藏开发前的完钻井;
②油藏内探井、评价井数少,且分布位置不具代表性,需首先分析油藏 开发后水体的推进情况,进而分析油藏勘探评价阶段的完钻井和投入开发后 的新钻井的含水饱和度(Sw)的变化情况,从而筛选出能代表油藏原始含水 饱和度的井位(简称“原始Sw井”)。
优选的,所述步骤4)中,在油藏油水界面上下,特别是油水界面之上的含 油范围内,原始Sw随油柱高度增加,呈现规律性降低的趋势,采用构建Sw 随油柱高度趋势变化函数,来表征油藏内原始Sw在垂向上的规律性变化,主 要包括以下几个方面:
①确定油柱0m处Sw值:根据测井解释油气水层的含水饱和度标准,确 定油水界面处即油柱高度0m处的Sw值;
②确定油藏最高部位Sw值:根据相渗资料、高部位井的Sw值以及区块 油藏经验认识,确定油藏最高部位及油柱最大高度处的Sw值;
③单井Sw值统计:统计原始Sw井的Sw特征值和油柱高度(H),油水界 面之上油柱高度为正值,油水界面之下油柱高度为负值。对于钻遇多个油砂 体的井的Sw特征值,可采用以下厚度加权平均算法进行计算,该井油柱高度 (H)取多个油砂体顶界距油水界面距离的算术平均值。
Figure BDA0002031304950000031
式中:Sw—某口井平均原始含水饱和度值;hn—某个油砂体厚度,单位m; Swn—某个油砂体平均原始含水饱和度值。
④建立Sw趋势变化函数:以油柱高度(H)为横坐标、原始含水饱和度(Sw) 为纵坐标,将以上数据进行拟合分析,从而建立Sw随H变化的趋势函数(SwT)。
优选的,所述步骤5)中:将二维趋势变化函数转化为受油水界面控制的 三维Sw变化趋势体,来表征油藏内部Sw的规律性变化。主要包括以下几个 方面:
①建立三维油柱高度体:根据步骤2)中确定的油藏油水界面深度(OWC), 联合油藏三维海拔深度体(Z),利用以下公式计算建立油藏三维油柱高度体。
Hv=Zv-OWC
式中:Hv—某油藏三维油柱高度体,m;OWC—某油藏油水界面海拔深度,m; Zv—某油藏三维海拔深度体,m;
②建立三维Sw趋势体:利用以上建立的油藏三维油柱高度体(Hv),应 用步骤4)中得到的原始Sw趋势变化函数SwT,计算建立R油田KI-3油藏三 维Sw趋势体(SwTv)。
优选的,所述步骤6)中,以井数据为硬数据、三维Sw变化趋势体为趋势 约束,在相控的基础上,建立油藏三维原始含水饱和度模型,所得到的Sw模 型既忠实于井点Sw数据,又大大提高了井间Sw预测精度。主要包括以下几 个方面:
①Sw井数据网格化:将步骤3)中确定的原始Sw井的测井解释Sw数据 粗化到三维地质模型井网格中,使原始SW井的测井解释Sw值作为硬数据参 与Sw模拟,确保模拟得到的原始Sw模型的井点数据与实际数据一致。
②相控:即沉积微相或岩相约束,将油藏三维地质建模过程中建立的沉 积微相或岩相模型以分区设定的方式,应用到Sw模拟过程中。其主要方法为: 首先通过原始Sw井的Sw数据统计分析,确定各微相或岩相的原始Sw数值分 布范围,然后在模拟过程中选用分相带区设定的控制方式,依次设定各微相 或岩相的原始Sw数值分布范围。通过相控的方法,确保油藏内不同微相、特 别是岩相中Sw的地质规律性差异。
③三维Sw变化趋势体约束:将步骤5)中构建的三维Sw变化趋势体以三 维趋势(Trend)约束的方式,应用到Sw模拟过程中,确保模拟得到的原始 Sw模型的Sw数值分布整体符合Sw变化趋势体表征的变化规律,从而合理预 测井间的原始Sw值,联合原始Sw井数据的硬约束和相控表征的不同地质体 差异,通过选择合适的算法,实现高精度原始Sw饱和度模型的建立。
本发明的有益效果是:该基于饱和度趋势变化函数的油藏含水饱和度模 拟方法设计合理,对油藏原始Sw模拟提供了一种切实有效的方法,弥补了以 往常数模型和经验公式模型的弊端,能合理表征油藏内油水垂向分布的地质 规律,大大提高了原始含水饱和度井间预测精度,井间吻合度达到90%以上, 实现了高精度原始含水饱和度模型的建立,满足了油田开发阶段油藏精细研 究和数值模拟研究的需求,对油藏含油气预测和新井部署也具有一定的指导 意义。同时,该方法适应于各种油藏类型,摆脱了油藏流体资料条件限制, 方法操作较为简单,具有较强的适用性和推广价值,适用于油藏开发阶段、 油藏三维地质建模过程中的油藏三维原始含水饱和度模拟,为油藏精细研究 和数值模拟研究的提供高精度地质模型,为油藏井位部署、剩余油挖潜、井 网调整提供支撑。
附图说明
图1为本发明流程示意图;
图2为本发明R油田主力层油藏剖面示意图;
图3为本发明R油田KI-3砂层组含油面积示意图;
图4为本发明R油田KI油组油层电性、含油性下限示意图;
图5为本发明R油田R1井KI-3砂层组相对渗透率曲线示意图;
图6为本发明R油田KI-3油藏H与原始Sw趋势函数拟合分析示意图;
图7为本发明R油田KI-3油藏三维海拔深度体示意图;
图8为本发明R油田KI-3油藏油柱高度体示意图;
图9为本发明R油田KI-3油藏三维Sw趋势体示意图;
图10为本发明R油田KI-3油藏岩相模型示意图;
图11为本发明R油田KI-3油藏原始含水饱和度模型示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行 清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而 不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做 出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一
R油田是一典型的断背斜构造油藏,地层产状与构造走向一致,主要开发 层系为KI油组,主要含油层位为KI和KII油组。KI-3砂层组是该油田主力 层之一,为具有独立油水系统的油藏,属中高孔-中高渗储层。该油田于2015 年4月26日投入开发,天然能量开发,目前综合含水38.1%,采出程度10.27%, 整体处于中含水、中低采出程度阶段。
请参阅图1~11,一种基于趋势变化函数的油藏三维原始含水饱和度模拟 方法,包括以下步骤:
1)基础数据收集整理
包括钻井、地质、地震、测井、化验、试油、试井、生产动态等基础及综 合解释研究成果数据。钻井数据重点为钻井井位坐标、补心海拔、完钻井时 间等信息;地质资料重点为含油面积图、油藏剖面图、地层对比图等综合研 究成果图件;测井数据重点为小层解释数据表、综合解释成果曲线等;生产 动态数据重点为各井历史日产油、日产水、含水等情况。
2)确定油藏油水界面深度及油柱高度
①根据油藏剖面图中油水界面所在的海拔深度,确定为油水界面深度 (OWC)。从图2为R油田油藏剖面图,从图中可直接读取KI-3油藏油水界 面深度为-685m。同时,图3为R油田K1-1砂体顶界含油面积图,由于该构 造图为油砂体顶界构造,因此图中油藏边界对应的深度即为油藏油水界面深 度,因此也可确定该油藏油水界面深度为-685m。
②根据油藏含油面积图中,油藏最高部位和最低部位深度差,确定为油柱 高度(H)。图3为R油田K1-1砂体顶界含油面积图,含油面积最高部位深 度为-570m,最低部位深度为-685m,从而确定R油田K1-1油藏油柱高度为115m。
3)确定参与油藏原始含水饱和度分析和模拟的井位
①油藏内探井、评价井数较多,分布位置具原始含水饱和度区域代表性, 则选取油藏开发前的完钻井;
②油藏内探井、评价井数少,且分布位置不具代表性,需首先分析油藏开 发后水体的推进情况,进而分析油藏勘探评价阶段的完钻井和投入开发后的 新钻井的含水饱和度(Sw)的变化情况,从而筛选出能代表油藏原始含水饱 和度的井位(简称“原始Sw井”)。
R油田在投入开发前仅有R1、R2和R1-1共3口井,其中R1井为探井,R2、 R1-1井为评价井,且R1-1井在油水边界之外,显然含油范围内这2口井对整 个油藏原始SW分布不具有区域代表性,需增加油藏开发后的部分井。KI-3砂 层组作为该油田主力层之一,在2015年4月投产时已动用,由于边水能量充 足,采用天然能量开发。投产后边水推进较快,低部位井含水上升快,产量 递减大,特别是构造低部位的R1-4、R1-8、R1-12、R1-2、R1-18、R1-19、R1-20 井,投产4个月,平均含水由初期的12.8%上升至25.6%,因此低部位在投产 后含水推进迅速,低部位投产后的新完钻井已不能代表KI-3油藏原始Sw。因 此在低部位要取该油藏投入开发前,即2015年4月26日之前的完钻井作为 参与原始饱和度分析和模拟的井位。而构造中部位井R1-9、R1-13、R1、R1-14、 R1-15、R1-10井在投产三年时间内含水变化不大,高部位井R1-11、R1-23井 截止到目前一直不含水,因此中部位需取投产三年时间内即2018年5月份之 前的井,高部位的完钻井均可取。从表1可以看到,R1-21、R1-22井为投产 3-4个月时间在低部位新钻的加密调整井,由于受低部位边水推进影响,这两 口井的Sw均不能油藏的原始Sw,因此不能参与油藏原始含水饱和度分析和模 拟;而R1-23井是投产1年半内在高部位新钻的加密调整井,该部位离边水 远,不受边水影响,可参与油藏原始含水饱和度分析和模拟;位于油水边界 之外的投产前的完钻井需全部参与模拟。因此该油藏除R1-21、R1-22井外, 其余井均需参与油藏原始含水饱和度分析和模拟,共确定19口原始Sw井。
4)构建原始含水饱和度随油柱高度的趋势变化函数
在油藏油水界面上下特别是油水界面之上的含油范围内,原始Sw随油柱高 度增加,呈现规律性降低的趋势,采用构建Sw随油柱高度趋势变化函数,来 表征油藏内原始Sw在垂向上的规律性变化,主要包括以下几个方面:
①确定油柱0m处Sw值:根据测井解释油气水层的含水饱和度标准,确 定油水界面处即油柱高度0m处的Sw值。从图4可以看出,KI-3油藏油层划 分标准为电阻率≥20.0Ω.m、声波时差≥215μs/m、孔隙度≥8%、含水饱和 度≤50%,因此确定油柱高度0m处的Sw=0.5;
②确定油藏最高部位Sw值:根据相渗资料、高部位井的Sw值以及区块 油藏经验认识,确定油藏最高部位及油柱最大高度处的Sw值。从图5可以看 出,该油藏束缚水饱和度为0.28,因此确定该区Sw≥0.28,而高部位井的Sw 值逐渐趋于0.28,且从该区块经验角度分析,该区油藏圈闭幅度高,油藏顶 部几乎无可动水,因此确定油藏最高部位,即最大油柱高度115m处Sw=0.28。
③单井Sw值统计:统计原始Sw井的Sw特征值和油柱高度(H),油水界 面之上油柱高度为正值,油水界面之下油柱高度为负值。对于钻遇多个油砂 体的井的Sw特征值,采用厚度加权平均算法即公式1进行计算,该井油柱高 度(H)取多个油砂体顶界距油水界面距离的算术平均值。R油田KI-3砂层组 油层发育,大部分井钻遇2个及以上油砂体,将每口井钻遇的油砂体的H进 行算术平均、Sw进行厚度加权平均值,作为各井的H和Sw特征值;R1-5井 无有效储层,所以无Sw特征值,从而得到以下R油田单井H和原始Sw统计 表。
R油田单井油柱高度H和原始Sw统计表
Figure BDA0002031304950000091
Figure BDA0002031304950000101
④建立Sw趋势变化函数:以油柱高度(H)为横坐标、原始含水饱和度 (Sw)为纵坐标,将以上数据进行拟合分析,从而建立Sw随H变化的趋势 函数(SwT)。将以上三步得到的数据点通过拟合分析(图6),从而确立R油 田KI-3油藏原始Sw随H变化的趋势函数(SwT)如下:
SwT=-2E-12H6+7E-10H5-7E-08H4+2E-06H3+0.0001H2-0.0094H+0.5,该公式 相关系数达到0.95以上,说明该函数拟合精度高。
5)建立油藏三维含水饱和度变化趋势体
①建立三维油柱高度体:根据步骤2)中确定的油藏油水界面深度(OWC) -685m,联合油藏三维海拔深度体(Z)(见图7,图中白色线为油水界面深度 线),利用公式计算建立R油田KI-3油藏三维油柱高度体(Hv)(见图8),图 8中0值线即油柱高度为0,即为油水界面;
②建立三维Sw趋势体:利用以上建立的油藏三维油柱高度体(Hv),应 用步骤4)中得到的原始Sw趋势变化函数SwT,计算建立R油田KI-3油藏三 维Sw趋势体(SwTv)(见图9)。从而将二维趋势变化函数转化为受油水界面 控制的三维Sw变化趋势体,来表征KI-3油藏内部原始Sw的规律性变化。
6)以井数据为硬数据、三维含水饱和度趋势体为趋势约束,建立油藏三 维原始含水饱和度模型。
①Sw井数据网格化:将步骤3)中确定的19口原始Sw井的测井解释Sw 数据粗化到三维地质模型井网格中,使原始Sw井的测井解释Sw值作为硬数 据参与Sw模拟,确保模拟得到的原始Sw模型的井点数据与实际数据一致;
②相控:该油藏三维地质模型建立的相模型为岩相模型(见图10),通过 原始Sw井的Sw数据统计分析,确定砂岩相原始Sw数值分布为0.28~0.8, 泥岩相原始Sw数值为常数值1。按照该数值区间,在原始Sw模拟过程中选 用分岩相区设定的控制方式,依次设定对各相进行数值范围设定,确保油藏 内岩相中Sw的地质规律性差异;
③三维Sw变化趋势体约束:将步骤5)中构建的三维Sw变化趋势体以 三维趋势(Trend)约束的方式,应用到Sw模拟过程中,确保模拟得到的原 始Sw模型的整体Sw数值分布符合Sw变化趋势体表征的变化规律,从而合 理预测井间的原始Sw值,在模拟联合原始Sw井数据的硬约束和相控表征的 不同地质体差异,采用常用的序贯高斯模拟的随机算法,建立了高精度R油 田KI-3砂层组原始Sw模型(见图11)。
工作原理:在使用该基于饱和度趋势变化函数的油藏含水饱和度模拟方 法时,应用该方法完成了R油田KI-3油藏三维原始含水饱和度模型的建立, 利用该模型进行数值模拟研究,明确了目前剩余油分布情况,指导开发方案 调整、优化井位部署、措施优选以及剩余油挖潜等工作。通过数值模拟研究 和生产动态分析目前油藏水体推进较为缓慢,构造低部位含水高平均48.5%, 中部位井含水低约15%左右,高部位井不含水,同时考虑到局部岩性变化的遮 挡作用,因此在中高部位剩余油富集区提出加密调整井位2口,分别为R1-24和R1-25井。其中R1-24井位于构造高部位,水体未推进;R1-25井位于构造 中部位,向低部位储层变差,干层具遮挡作用,且邻井井距大,采出程度低, 水体未波及,目前两口井均已实施。从模型误差分析上看,两口井在本方法 所的建立原始Sw饱和度模型中预测的顶部砂体Sw分别为0.288和0.335,实 际钻遇顶部砂体Sw为0.286和0.352,误差分别为0.7%和4.8%,误差均在 5%以内;从单井生产效果上看,两口井初期单井日产均在500BBL以上且不含 水,与数值模拟预测结果一致,取得了较好的效果。可见该方法合理有效, 所建立的油藏三维原始含水饱和度模型精度高,能合理预测井间原始含水饱 和度分布,对油藏精细研究和开发具有直接的指导意义。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节, 而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实 现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且 是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨 在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。 不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实 施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起 见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也 可以经适当组合,形成本领域技术人员可以理解的其他实施方式。

Claims (3)

1.一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法,其特征在于:包括以下步骤:
1)基础数据收集整理;
2)确定油藏油水界面海拔深度及油柱高度;
3)确定参与油藏原始含水饱和度分析和模拟的井位;
4)构建原始含水饱和度随油柱高度的趋势变化函数;
所述步骤4)中,油水界面之上的含油范围内,原始Sw随油柱高度增加,呈现规律性降低的趋势,采用构建Sw随油柱高度趋势变化函数,来表征油藏内原始Sw在垂向上的规律性变化,包括以下几个方面:
①确定油柱0m处Sw值:根据测井解释油气水层的含水饱和度标准,确定油水界面处即油柱高度0m处的Sw值;
②确定油藏最高部位Sw值:根据相渗资料、高部位井的Sw值以及区块油藏经验认识,确定油藏最高部位及油柱最大高度处的Sw值;
③单井Sw值统计:统计原始Sw井的Sw特征值和油柱高度H,油水界面之上油柱高度为正值,油水界面之下油柱高度为负值,对于钻遇多个油砂体的井的Sw特征值,采用厚度加权平均算法进行计算,该井油柱高度H取多个油砂体顶界距油水界面距离的算术平均值;
Figure FDA0003538913050000011
式中:Sw—某口井平均原始含水饱和度值;hn—某个油砂体厚度,单位m;Swn—某个油砂体平均原始含水饱和度值;
④建立Sw趋势变化函数:以油柱高度H为横坐标、原始含水饱和度Sw为纵坐标,将以上数据进行拟合分析,从而建立Sw随H变化的趋势函数SwT;
5)以趋势变化函数为约束,建立油藏三维含水饱和度变化趋势体;
所述步骤5)中:将二维趋势变化函数转化为受油水界面控制的三维Sw变化趋势体,来表征油藏内部Sw的规律性变化,包括以下几个方面:
①建立三维油柱高度体:根据步骤2)中确定的油藏油水界面深度OWC,联合油藏三维海拔深度体Z,利用公式2计算建立油藏三维油柱高度体,
Hv=Zv-OWC
式中:Hv—某油藏三维油柱高度体,m;OWC—某油藏油水界面海拔深度,m;Zv—某油藏三维海拔深度体,m;
②建立三维Sw趋势体:利用以上建立的油藏三维油柱高度体Hv,应用步骤4)中得到的原始Sw趋势变化函数SwT,计算建立R油田KI-3油藏三维Sw趋势体SwTv;
6)以井数据为硬数据、三维含水饱和度趋势体为趋势约束,建立油藏三维原始含水饱和度模型;
所述步骤6)中,以井数据为硬数据、三维Sw变化趋势体为趋势约束,在相控的基础上,建立油藏三维原始含水饱和度模型,所得到的Sw模型既忠实于井点Sw数据,又大大提高了井间Sw预测精度,包括以下几个方面:
①Sw井数据网格化:将步骤3)中确定的原始Sw井的测井解释Sw数据粗化到三维地质模型井网格中,使原始Sw井的测井解释Sw值作为硬数据参与Sw模拟,确保模拟得到的原始Sw模型的井点数据与实际数据一致;
②相控:即沉积微相或岩相约束,将油藏三维地质建模过程中建立的沉积微相或岩相模型以分区设定的方式,应用到Sw模拟过程中,方法为:首先通过原始Sw井的Sw数据统计分析,确定各微相或岩相的原始Sw数值分布范围,然后在模拟过程中选用分相带区设定的控制方式,依次设定各微相或岩相的原始Sw数值分布范围,通过相控的方法,确保油藏内不同微相、特别是岩相中Sw的地质规律性差异;
③三维Sw变化趋势体约束:将步骤5)中构建的三维Sw变化趋势体以三维趋势Trend约束的方式,应用到Sw模拟过程中,确保模拟得到的原始Sw模型的Sw数值分布整体符合Sw变化趋势体表征的变化规律,从而合理预测井间的原始Sw值,联合原始Sw井数据的硬约束和相控表征的不同地质体差异,通过选择合适的算法,实现高精度原始Sw饱和度模型的建立。
2.根据权利要求1所述的一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法,其特征在于:所述步骤2)中,确定油藏油水界面海拔深度及油柱高度,包括以下几个方面:
①根据油藏剖面图中油水界面所在的海拔深度,确定为油水界面海拔深度;
②根据油藏含油面积图中,含油范围最高部位和最低部位深度差,确定为油柱高度。
3.根据权利要求1所述的一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法,其特征在于:所述步骤3)中,确定参与油藏原始含水饱和度分析和模拟的井位,包括以下几个方面:
①油藏内探井、评价井数较多,分布位置具原始含水饱和度区域代表性,则选取油藏开发前的完钻井;
②油藏内探井、评价井数少,且分布位置不具代表性,需首先分析油藏开发后水体的推进情况,进而分析油藏勘探评价阶段的完钻井和投入开发后的新钻井的含水饱和度的变化情况,从而筛选出能代表油藏原始含水饱和度的井位。
CN201910310522.4A 2019-04-17 2019-04-17 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法 Active CN110056346B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910310522.4A CN110056346B (zh) 2019-04-17 2019-04-17 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910310522.4A CN110056346B (zh) 2019-04-17 2019-04-17 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法

Publications (2)

Publication Number Publication Date
CN110056346A CN110056346A (zh) 2019-07-26
CN110056346B true CN110056346B (zh) 2022-05-17

Family

ID=67319303

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910310522.4A Active CN110056346B (zh) 2019-04-17 2019-04-17 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法

Country Status (1)

Country Link
CN (1) CN110056346B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110929383A (zh) * 2019-10-28 2020-03-27 中国石油大港油田勘探开发研究院 一种源下油气柱高度的计算方法
CN111274694B (zh) * 2020-01-17 2022-10-04 中国石油天然气股份有限公司 一种基于含油边界和油水过渡带约束的饱和度建模方法
CN111561312B (zh) * 2020-04-21 2023-09-26 中国石油天然气股份有限公司 基于常规测井数据的饱和度计算方法及装置
CN111550238B (zh) * 2020-05-07 2022-10-04 中国石油天然气股份有限公司 一种等效j函数约束的饱和度建模方法
CN113803054B (zh) * 2020-06-11 2024-04-09 中国石油化工股份有限公司 油水界面深度确定方法及防止生产井水窜的预警方法
CN111898091B (zh) * 2020-06-22 2023-09-26 中国石油天然气股份有限公司 油水界面的确定方法及装置

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4644283A (en) * 1984-03-19 1987-02-17 Shell Oil Company In-situ method for determining pore size distribution, capillary pressure and permeability
GB0008358D0 (en) * 2000-04-06 2000-05-24 Bastian Juergen Non-breathing power transformer without nitrogen blanket
CN101487390A (zh) * 2009-02-23 2009-07-22 大庆油田有限责任公司 一种确定油层原始含油饱和度的阿尔奇模式方法
CN106204302A (zh) * 2016-07-05 2016-12-07 中国石油天然气股份有限公司 一种计算原始含水饱和度的方法与应用
CN107701180A (zh) * 2017-09-19 2018-02-16 中国石油天然气集团公司 一种基于密闭取心的原始油藏含水饱和度计算方法
CN108241785A (zh) * 2018-01-11 2018-07-03 中国海洋石油集团有限公司 一种非均质性储层饱和度场精细表征方法
CN109611087A (zh) * 2018-12-11 2019-04-12 中国石油大学(北京) 一种火山岩油藏储层参数智能预测方法及系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EA015208B1 (ru) * 2006-09-08 2011-06-30 Бп Эксплорейшн Оперейтинг Компани Лимитед Способ повышения нефтеотдачи продуктивного пласта
US10495774B2 (en) * 2014-01-13 2019-12-03 Schlumberger Technology Corporation Method for estimating irreducible water saturation from mercury injection capillary pressure
CN105445441B (zh) * 2015-12-03 2017-06-13 中国石油天然气股份有限公司 储层含气饱和度确定方法及装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4644283A (en) * 1984-03-19 1987-02-17 Shell Oil Company In-situ method for determining pore size distribution, capillary pressure and permeability
GB0008358D0 (en) * 2000-04-06 2000-05-24 Bastian Juergen Non-breathing power transformer without nitrogen blanket
CN101487390A (zh) * 2009-02-23 2009-07-22 大庆油田有限责任公司 一种确定油层原始含油饱和度的阿尔奇模式方法
CN106204302A (zh) * 2016-07-05 2016-12-07 中国石油天然气股份有限公司 一种计算原始含水饱和度的方法与应用
CN107701180A (zh) * 2017-09-19 2018-02-16 中国石油天然气集团公司 一种基于密闭取心的原始油藏含水饱和度计算方法
CN108241785A (zh) * 2018-01-11 2018-07-03 中国海洋石油集团有限公司 一种非均质性储层饱和度场精细表征方法
CN109611087A (zh) * 2018-12-11 2019-04-12 中国石油大学(北京) 一种火山岩油藏储层参数智能预测方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种储层三维地质建模技术方法-以哈萨克斯坦K油田为例;杨永亮;《承德石油高等专科学校学报》;20160630;第18卷(第13期);第4-7页 *
杏北开发区东部过渡带中块剩余油挖潜潜力及加密方式研究;徐文博;《中国优秀硕士学位论文数据库工程科技1辑》;20190115;全文 *

Also Published As

Publication number Publication date
CN110056346A (zh) 2019-07-26

Similar Documents

Publication Publication Date Title
CN110056346B (zh) 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法
CN109061765B (zh) 非均质薄砂岩互层油藏的圈闭评价方法
CN106875471B (zh) 煤系含或隔水层三维可视化建模方法
US20160376885A1 (en) Method and Apparatus for Performance Prediction of Multi-Layered Oil Reservoirs
CN104809277A (zh) 一种超低渗致密储层水平井地质建模方法
CN105822302A (zh) 一种基于井地电位法的油水分布识别方法
CN109653725A (zh) 基于沉积微相和岩石相的混积储层水淹程度测井解释方法
CN109541685B (zh) 一种河道砂体识别方法
CN110244021B (zh) 一种基于各向异性插值的地层分层方法
CN109783944B (zh) 三维流动单元模型建立方法、装置和设备
CN109143399B (zh) 一种识别碳酸盐岩层序界面的方法
Harris The role of geology in reservoir simulation studies
CN106570262A (zh) 一种储层构型结构的描述方法
CN111706317B (zh) 一种确定加密调整区低渗储层剩余油分布状况的方法
CN111706318B (zh) 一种确定低渗储层剩余油分布状况的方法
Liu et al. The Control Theory and Application for Well Pattern Optimization of Heterogeneous Sandstone Reservoirs
CN113742962B (zh) 一种基于水平井的页岩储层属性三维建模方法
CN110309597A (zh) 基于构型界面的阻流带模型确定方法、装置及存储介质
CN115220101A (zh) 一种深层碎屑岩砂体沉积结构的建模方法
CN111027780B (zh) 油势场模拟方法及装置
Peijs-van Hilten et al. Heterogeneity modeling and geopseudo upscaling applied to waterflood performance prediction of an incised valley reservoir: Countess YY Pool, southern Alberta, Canada
Pandey et al. Geocellular Model for Tertiary Reservoirs in Manhera Tibba Gas Field, Jaisalmer Basin, Rajasthan, India
Li et al. Multiscale modeling of meandering fluvial reservoir architecture based on multiple-point geostatistics: a case study of the minghuazhen formation, yangerzhuang oilfield, bohai bay basin, China
CN115422789B (zh) 一种考虑全过程优化断块油藏水驱采收率预测方法及系统
CN114429016A (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
TA01 Transfer of patent application right

Effective date of registration: 20201130

Address after: 100000 Beijing city Dongcheng District No. 16 Andrew

Applicant after: PetroChina Co.,Ltd.

Address before: 300280 Yili 79, Haibin Street Research Institute District, Dagang District, Tianjin

Applicant before: DAGANG OIL FIELD OF CNPC

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant