CN105089582B - 基于井下流量控制设备的油藏数值模拟方法及装置 - Google Patents

基于井下流量控制设备的油藏数值模拟方法及装置 Download PDF

Info

Publication number
CN105089582B
CN105089582B CN201510284080.2A CN201510284080A CN105089582B CN 105089582 B CN105089582 B CN 105089582B CN 201510284080 A CN201510284080 A CN 201510284080A CN 105089582 B CN105089582 B CN 105089582B
Authority
CN
China
Prior art keywords
well
flow
reservoir
numerical
models
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
CN201510284080.2A
Other languages
English (en)
Other versions
CN105089582A (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.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas 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 China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201510284080.2A priority Critical patent/CN105089582B/zh
Publication of CN105089582A publication Critical patent/CN105089582A/zh
Application granted granted Critical
Publication of CN105089582B publication Critical patent/CN105089582B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明提供了一种基于井下流量控制设备的油藏数值模拟方法及装置,该油藏数值模拟方法包括:获取生产井、注入井、油藏及流量控制设备数据,确定各井、油藏及流量控制设备的动态参数及静态参数;根据每口井的井身结构和流量控制设备创建多井段模型;根据油藏特征及流体属性建立油藏渗流模型,所述油藏渗流模型包括:油气水三相渗流方程、外边界条件以及约束条件;将所述多井段模型及油藏渗流模型耦合,生成数值模型;根据油藏的静态参数建立地质模型,并根据原始地层压力分布或者压力与深度的关系,以及岩石、流体属性对所述地质模型进行初始化;根据所述数值模型及初始化后的地质模型建立带有稀疏系数矩阵的非线性方程组并求解所述非线性方程组。

Description

基于井下流量控制设备的油藏数值模拟方法及装置
技术领域
[0001] 本发明涉及油藏数值模拟技术,尤其涉及一种基于井下流量控制设备的油藏数值 模拟方法及装置。
背景技术
[0002] 长期以来,油藏数值模拟作为油气田开发过程中一项可信赖的技术,在油藏工程 师做大型决策、评价油藏以及诊断和改善油藏产能时发挥了不可替代的作用。油藏数值模 拟中任何井模型的基本目的都是提供准确的源汇项。近年来,随着中国水驱陆相油田开始 普遍进入高含水和高采出程度的“双特高”开发阶段,特低渗、低产、低丰度、高粘度边际油 田的开发逐步得到重视,井眼轨迹和井身结构复杂以及安装了井下监测、调控设备的多分 支井、水平井和波状井成为现阶段老油田挖潜增产、新油田经济开发的重要技术。但是现有 的商业化油藏数值模拟软件中的井模型是一个常量源汇项,采用固定的井底流压或产量作 为内边界条件,难以模拟井下流量控制设备的影响。
[0003] 针对以上难题,赵国忠等(赵国忠,孙巍,何鑫•基于分层注水数学模型的油藏数值 模拟•东北石油大学学报,2012,36 (6) :82-87.)于2〇12年针对大庆油田分层注水工艺提出 了考虑水嘴影响的油藏数值模拟方法,主要是通过建立考虑启动压差、水嘴直径等嘴损特 征的分层注水模型,计算分层注水井各段嘴前流压,实现嘴损方程与井一网格压力方程的 親合。这种方法的局限性在于:①对于井身结构复杂的水平井和多分支井不适用;②对于井 筒内存在多个水嘴的注水井不适用;③分层注水工艺需要使用封隔器,模型中未考虑其影 响;④仅考虑了水嘴这一种井下流量控制设备,仅适用于注水井,对于油井尚不适用。
发明内容
[0004]本发明提供一种基于井下流量控制设备的油藏数值模拟方法,以对采取不同井型 分层注水、分层采油工艺的油藏进行高精度数值模拟,为科学的使用井下流量控制设备提 供依据,达到提高采收率的目的。
[0005]为了实现上述目的,本发明提供一种基于井下流量控制设备的油藏数值模拟方 法,所述的油藏数值模拟方法包括:
[0006]获取生产井、注入井、油藏及流量控制设备数据,确定各井、油藏及流量控制设备 的动态参数及静态参数;
[0007]根据每口井的井身结构和流量控制设备创建多井段模型;
[0008]根据油藏特征及流体属性建立油藏渗流模型,所述油藏渗流模型包括:油气水三 相渗流方程、外边界条件以及约束条件;
[0009]将所述多井段模型及油藏渗流模型耦合,生成数值模型;
[0010]根据油藏的静态参数建立地质模型,并根据原始地层压力分布或者压力与深度的 关系,以及岩石、流体属性对所述地质模型进行初始化;
[0011]根据所述数值模型及初始化后的地质模型建立带有稀疏系数矩阵的非线性方程 组并求解所述非线性方程组。
[0012]在一实施例中,根据每口井的井身结构和流量控制设备创建多井段模型,包括: [0013]根据井眼轨迹,流量控制设备的安装使用位置以及所需模拟精度确定所述多井段 模型的段数及节点数;
[00M]根据所述流量控制设备的使用方式确定各井段的分支数目、流体在该井段的流入 位置及流出口位置;
[0015]根据所述段数、节点数、各井段的分支数目、流体在该井段的流入位置及流出口位 置确定各井段、节点以及分支的位置,搭建多井段模型的骨架;
[0016]根据单井、油藏及所述流量控制设备的静态参数确定每井段、分支及节点的深度; [0017]根据油管直径和套管直径及粗糙度确定各井段的直径、截面积、体积及粗糙度; [0018]根据模拟精度、井身结构和流量控制设备动静态参数确定最长段长度和最小段长 度,使多井段模型内的各井段均短于最长段长度,均长于最小段长度。
[0019]在一实施例中,在根据每口井的井身结构和流量控制设备创建多井段模型之后, 所述的油藏数值模拟方法包括:根据流入、流出口位置、井底流压、每段、分支及节点的深 度、油管直径和套管直径及粗糙度确定各段的直径、截面积、体积及粗糙度选择井筒压力计 算方法。
[0020]在一实施例中,将所述多井段模型及油藏渗流模型耦合,生成数值模型,包括: [0021]根据所述多井段模型,建立各井段的各分支内的油气水三相流动方程、压力方程;
[0022] 将所述多井段模型作为内边界条件与油藏渗流模型耦合,建立综合数学模型;
[0023] 对所述综合数学模型进行离散化,形成数值模型。
[0024]在一实施例中,根据流量方程计算某一相流体从网格块到井段节点的流动,所述 流量方程如下:
[0025] qPj = TwjMpj (Pj+Hcj-Pn-Hnc)
[0026]其中,qPj是p相在网格块与井段连接j处的流量,Twj是连接传导率,MPj是p相在网格 块与井段连接j处的流度,P」是网格块压力,H。」是网格块到连接处的深度差,Pn是段节点处 的压力,Hn。是段节点到连接处的深度差。
[0027]在一实施例中,所述静态参数包括:每口井的井位坐标、完井方式、射孔深度、井 头、井深、井斜及井轨迹数据,油藏目的层位的构造形态、测井解释数据、地震解释数据、断 层数据,目的层位的孔隙度、渗透率、含油饱和度及原始地层压力参数,地层内岩石及流体 的高压物性资料,流量控制设备的类型、安装位置、设备尺寸及控制方式。
[0028]在一实施例中,所述动态参数包括:单井生产动态数据、吸水剖面测试数据、产油 剖面测试数据、井下流量控制设备的调节参数及井下流量控制设备监测数据。
[0029]为了实现上述目的,一种基于井下流量控制设备的油藏数值模拟装置,所述的油 藏数值模拟装置包括:
[0030]参数获取单元,用于获取生产井、注入井、油藏及流量控制设备数据,确定各井、油 藏及流量控制设备的动态参数及静态参数;
[0031]多井段模型创建单元,用于根据每口井的井身结构和流量控制设备创建多井段模 型;
[0032]油藏渗流模型生成单元,用于根据油藏特征及流体属性建立油藏渗流模型,所述 油藏渗流模型包括:油气水三相渗流方程、外边界条件以及约束条件;
[0033]数值模型生成单元,用于将所述多井段模型及油藏渗流模型耦合,生成数值模型; [0034]地质模型生成单元,用于根据油藏的静态参数建立地质模型,并根据原始地层压 力分布或者压力与深度的关系,以及岩石、流体属性对所述地质模型进行初始化;
[0035]稀疏矩阵生成单元,用于根据所述数值模型及初始化后的地质模型建立带有稀疏 系数矩阵的非线性方程组并求解所述非线性方程组。
[0036]在一实施例中,所述多井段模型创建单元包括:
[0037]第一信息确定模块,用于根据井眼轨迹,流量控制设备的安装使用位置以及所需 模拟精度确定所述多井段模型的段数及节点数;
[0038]第二信息确定模块,用于根据所述流量控制设备的使用方式确定各井段的分支数 目、流体在该井段的流入位置及流出口位置;
[0039]第三信息确定模块,用于根据所述段数、节点数、各井段的分支数目、流体在该井 段的流入位置及流出口位置确定各井段、节点以及分支的位置,搭建多井段模型的骨架; [0040]第四信息确定模块,用于根据单井、油藏及所述流量控制设备的静态参数确定每 井段、分支及节点的深度;
[0041]第五信息确定模块,用于根据油管直径和套管直径及粗糙度确定各井段的直径、 截面积、体积及粗糙度;
[0042]第六信息确定模块,用于根据模拟精度、井身结构和流量控制设备的动静态参数 确定最长段长度和最小段长度,使多井段模型内的各井段均短于最长段长度,均长于最小 段长度。
[0043]在一实施例中,所述油藏数值模拟装置还包括:计算方法选取单元,用于根据流 入、流出口位置、井底流压、每段、分支及节点的深度、油管直径和套管直径及粗糙度确定各 段的直径、截面积、体积及粗糙度选择井筒压力计算方法。
[0044]在一实施例中,所述数值模型生成单元包括:
[0045]方程建立模块,用于根据所述多井段模型,建立各井段的各分支内的油气水三相 流动方程、压力方程;
[0046]数学模型建立模块,用于将所述多井段模型作为内边界条件与油藏渗流模型耦 合,建立综合数学模型;
[0047]数值模型生成模块,用于对所述综合数学模型进行离散化,形成数值模型。
[0048]在一实施例中,所述油藏数值模拟装置还包括:流动信息计算单元,用于根据流量 方程计算某一相流体从网格块到井段节点的流动,所述流量方程如下:
[0049] qPj = TwjMpj (Pj+Hcj-Pn-Hnc)
[0050]其中,qpj是p相在网格块与井段连接j处的流量,Twj是连接传导率,Mpj是p相在网格 块与井段连接j处的流度,Pj是网格块压力,知是网格块到连接处的深度差,pn是段节点处 的压力,Hn。是段节点到连接处的深度差。
[0051]在一实施例中,所述静态参数包括:每口井的井位坐标、完井方式、射孔深度、井 头、井深、井斜及井轨迹数据,油藏目的层位的构造形态、测井解释数据、地震解释数据、断 层数据,目的层位的孔隙度、渗透率、含油饱和度及原始地层压力参数,地层内岩石及流体 的高压物性资料,流量控制设备的类型、安装位置、设备尺寸及控制方式。
[0052]在一实施例中,所述动态参数包括:单井生产动态数据、吸水剖面测试数据、产油 剖面测试数据、井下流量控制设备的调节参数及井下流量控制设备监测数据。
[0053]本发明实施例的有益效果在于,通过本发明,实现了对采取不同井型分层注水、分 层采油工艺的油藏进行高精度数值模拟,为科学的使用井下流量控制设备提供了依据,达 到提高了采收率的目的。针对带有井下流量控制设备的直井、水平井及复杂结构井,本发明 能够模拟水嘴、封隔器、栗等多种井下流量控制设备,可适用于生产井和注入井。本发明考 虑井筒内的管流特征以及流量控制设备的压损特征计算压力分布,提高了油藏数值模拟技 术中井的模拟精度。
附图说明
[0054]为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本 发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以 根据这些附图获得其他的附图。
[0055]图1为本发明一实施例的油藏数值模拟方法流程图;
[0056]图2为本发明实施例的多井段模型创建方法流程图;
[0057]图3为本发明实施例中考虑井下流量控制设备的多段井模型基本骨架示意图 [0058]图4为本发明实施例的数值模型生成方法流程图;
[0059]图5为本发明实施例中生产/注入井油套环空内带有封隔器和流量控制设备时的 流动示意图;
[0060]图6为本发明另一实施例的油藏数值模拟方法流程图;
[0061]图7为本发明实施例的油藏数值模拟装置的结构框图;
[0062]图8为本发明实施例的多井段模型创建单元702的结构框图;
[0063]图9为本发明实施例的数值模型生成单元704的结构框图。
具体实施方式
[0064]下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完 整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于 本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例,都属于本发明保护的范围。
[0065]为了解决现有技术中存在的问题,本发明实施例提供一种基于井下流量控制设备 的油藏数值模拟方法,如图1所示,所述的油藏数值模拟方法包括:
[0066]步骤101:获取生产井、注入井、油藏及流量控制设备数据,确定各井、油藏及流量 控制设备的动态参数及静态参数;
[0067]步骤102:根据每口井的井身结构和流量控制设备创建多井段模型;
[0068]步骤103:根据油藏特征及流体属性建立油藏渗流模型,所述油藏渗流模型包括: 油气水三相渗流方程、外边界条件以及约束条件;
[0069]步骤104:将所述多井段模型及油藏渗流模型耦合,生成数值模型;
[0070]步骤105:根据油藏的静态参数建立地质模型,并根据原始地层压力分布或者压力 与深度的关系,以及岩石、流体属性对所述地质模型进行初始化;
[0071]步骤106:根据所述数值模型及初始化后的地质模型建立带有稀疏系数矩阵的非 线性方程组并求解所述非线性方程组。
[0072]由图1所示的流程可知,本发明首先根据每口井的井身结构和流量控制设备创建 多井段模型,根据油藏特征及流体属性建立油藏渗流模型,将多井段模型及油藏渗流模型 耦合为数值模型,然后根据油藏的静态参数建立地质模型并初始化,最后根据所述数值模 型及初始化后的地质模型建立带有稀疏系数矩阵的非线性方程组并求解所述非线性方程 组。通过本发明,可以实现对采取不同井型分层注水、分层采油工艺的油藏进行高精度数值 模拟,为科学的使用井下流量控制设备提供了依据,达到提高了采收率的目的。
[0073]步骤101中,静态参数包括:每口井的井位坐标、完井方式、射孔深度、井头、井深、 井斜及井轨迹数据,油藏目的层位的构造形态(顶面、底面构造图)、测井解释数据、地震解 释数据、断层数据,目的层位的孔隙度、渗透率、含油饱和度及原始地层压力参数,地层内岩 石和流体的高压物性资料(相渗曲线、毛冠压力曲线、PTV资料),流量控制设备的类型、安装 位置、设备尺寸和控制方式。
[0074]动态参数包括:单井生产动态数据(产油量、产水量、产气量、气油比、含水率、油 压、套压、井底流压、配注量)、吸水剖面测试数据、产油剖面测试数据、井下流量控制设备的 调节参数、井下流量控制设备监测数据。
[0075]如图2所示,步骤102具体实施时,包括如下步骤:
[0076]步骤201:根据井眼轨迹饱括井段长度和井身结构),流量控制设备的安装使用位 置以及所需模拟精度确定所述多井段模型的段数及节点数。
[0077]多井段模型的段数及节点数越多模拟精度越高,同时计算时间越长,数据丰富且 区块内安装有井下流量控制设备的井数不多时,可适当增加段数及节点数,以提高模拟精 度。
[0078]步骤202:根据所述流量控制设备的使用方式确定各井段的分支数目、流体在该井 段的流入位置及流出口位置。
[0079]井段用以表征井筒内的流动,分支用以表征井下流量控制设备和油套环空中的流 动。分支的数目需要根据模拟精度的需要来确定。一个井段存在任意数量的出口,并且可以 利用井段和分支形成环形流动路径,环空中流体可以被封隔器堵塞,也可以流经多个流量 控制设备。
[0080]步骤2〇3:根据所述段数、节点数、各井段的分支数目、流体在该井段的流入位置及 流出口位置确定各井段、节点以及分支的位置,搭建多井段模型的骨架,图3为本发明实施 例中考虑井下流量控制设备的多段井模型基本骨架示意图,其中标号丨为流量控制设备,标 号2为封隔器,标号3为油套环空,标号4为油管。
[0081]多段井模型中,顶部井段相对特殊,该井段的节点对应井底流压测试深度,该井段 的压力即为井底流压。该井段属于初始段不需要计算压力损失,其压力方程被井生产方式 提供的边界条件取代。
[0082]步骤204:根据单井、油藏及所述流量控制设备的静态参数确定每井段、分支及节 点的深度,包括垂直深度和补心海拔。
[0083]步骤2〇5:根据油管直径和套管直径及粗糙度确定各井段的直径、截面积、体积及 粗糙度,这些数据用以计算井段内的压力分布。
[0084]步骤2〇6:根据模拟精度、井身结构和流量控制设备动静态参数确定最长段长度和 最小段长度,使多井段模型内的各井段均短于最长段长度,均长于最小段长度。
[0085] —实施例中,在根据每口井的井身结构和流量控制设备创建多井段模型之后,还 需要根据流入、流出口位置、井底流压、每段、分支及节点的深度、油管直径和套管直径及粗 糙度确定各段的直径、截面积、体积及粗糙度选择井筒压力计算方法,具体实施时:
[0086]首先,需要检测提供的参数所能满足的井筒压力计算方法,备选方法包括:①基于 油气水三相均匀流动模型计算井筒压力分布;②基于井筒管流表(VFP)计算井筒压力分布; ③基于漂移流量滑移模型计算三相非均匀流动时的井筒压力分布;④基于流量控制设备压 损模型计算流经设备后的压力分布。
[0087]基于油气水三相均匀流动模型的计算方法假设各相之间没有滑移,没有速度差, 压力计算主要根据段内平均流体密度。基于VFP表计算方法主要是为模型计算出所需的举 升曲线,可以处理垂直井和水平井,可以考虑气举、井底泵、气体压缩机以及地面油嘴的影 响,具备黑油和组分流体的计算能力。基于漂移流量滑移模型的计算方法则可以考虑各相 之间的滑移。
[0088]其次,根据所需模拟选择井筒压力计算方法,四种压力计算方法中,方法④是必须 具备的,其它三种方法的模拟精度为:②〉③〉①。
[0089] —实施例中,如图4所示,步骤104中,将所述多井段模型及油藏渗流模型耦合,生 成数值模型,包括:
[0090]步骤401:根据所述多井段模型,建立各井段的各分支内的油气水三相流动方程、 压力方程;
[0091]步骤4〇2:将所述多井段模型作为内边界条件与油藏渗流模型耦合,建立综合数学 模型;
[0092]步骤4〇3:对所述综合数学模型进行离散化,形成数值模型。
[0093]图5为本发明实施例中生产/注入井油套环空内带有封隔器和流量控制设备时的 流动示意图,其中51代表设备分支,52代表油管段。可以根据流量方程计算某一相流体从网 格块到井段节点的流动,流量方程如下:
[0094] Qpj = TwjMpj (Pj+Hcj—Pn—Hnc)
[0095] 其中,qpj是p相在网格块与井段连接j处的流量,Twj是连接传导率,MPj是p相在网格 块与井段连接j处的流度,Pj是网格块压力,Hcj是网格块到连接处的深度差,Pn是段节点处 的压力,Hn。是段节点到连接处的深度差。
[0096] —实施例中,图1所示的步骤105中,需要设置非线性分析控制参数:设置计算精 度、迭代求解次数初始化计算数据(初始化地质模型的三维压力及饱和度场数据和时间数 据)。
[0097]步骤106中,根据所述数值模型及初始化后的地质模型建立带有稀疏系数矩阵的 非线性方程组并求解所述非线性方程组时,时间参数更新及初始化赋值如下:计算完某一 时间步后,把这一时间步计算的压力与饱和度数据赋值到下一时间步的各个网格及井段、 分支节点,作为下一时间步的初始值。
[0098]图1、图2及4所示的基于井下流量控制设备的油藏数值模拟方法可以通过图6进行 概括,分为左右两个分支流程,图6中已对两个分支进行了描述,在此不再赘述D [0099]本发明实施例的有益效果在于,通过本发明,实现了对采取不同井型分层注水、分 层采油工^的油藏进行高精度数值模拟,为科学的使用井下流量控制设备提供了依据,达 到提高了采收率的目的。针对带有井下流量控制设备的直井、水平井及复杂结构井,本发明 能够模拟水嘴、封隔器、栗等多种井下流量控制设备,可适用于生产井和注入井。本发明考 虑井筒内的管流特征以及流量控制设备的压损特征计算压力分布,提高了油藏数值模拟技 术中井的模拟精度。
[0100]本发明实施例提供一种基于井下流量控制设备的油藏数值模拟装置,如图7所示, 所述的油藏数值模拟装置包括:参数获取单元7〇1,多井段模型创建单元702,油藏渗流模型 生成单元703,数值模型生成单元704,地质模型生成单元705及稀疏矩阵生成单元706。
[0101]参数获取单元701用于获取生产井、注入井、油藏及流量控制设备数据,确定各井、 油藏及流量控制设备的动态参数及静态参数;
[0102]多井段模型创建单元702用于根据每口井的井身结构和流量控制设备创建多井段 模型;
[0103] 油藏渗流模型生成单元703用于根据油藏特征及流体属性建立油藏渗流模型,所 述油藏渗流模型包括:油气水三相渗流方程、外边界条件以及约束条件;
[0104] 数值模型生成单元7〇4用于将所述多井段模型及油藏渗流模型耦合,生成数值模 型;
[0105] 地质模型生成单元7〇5用于根据油藏的静态参数建立地质模型,并根据原始地层 压力分布或者压力与深度的关系,以及岩石、流体属性对所述地质模型进行初始化;
[0106]稀疏矩阵生成单元706用于根据所述数值模型及初始化后的地质模型建立带有稀 疏系数矩阵的非线性方程组并求解所述非线性方程组。
[01 07]其中,静态参数包括:每口井的井位坐标、完井方式、射孔深度、井头、井深、井斜及 井轨迹数据,油藏目的层位的构造形态(顶面、底面构造图)、测井解释数据、地震解释数据、 断层数据,目的层位的孔隙度、渗透率、含油饱和度及原始地层压力参数,地层内岩石和流 体的高压物性资料(相渗曲线、毛冠压力曲线、PTV资料),流量控制设备的类型、安装位置、 设备尺寸和控制方式。
[0108] 动态参数包括:单井生产动态数据(产油量、产水量、产气量、气油比、含水率、油 压、套压、井底流压、配注量)、吸水剖面测试数据、产油剖面测试数据、井下流量控制设备的 调节参数、井下流量控制设备监测数据。
[0109] 一实施例中,如图8所示,所述多井段模型创建单元702包括:第一信息确定模块 801,第二信息确定模块802,第三信息确定模块8〇3,第四信息确定模块804,第五信息确定 模块805,第六信息确定模块806。
[0110]第一信息确定模块801用于根据井眼轨迹,流量控制设备的安装使用位置以及所 需模拟精度确定所述多井段模型的段数及节点数;
[0111] 第二信息确定模块802用于根据所述流量控制设备的使用方式确定各井段的分支 数目、流体在该井段的流入位置及流出口位置;
[0112] 第三信息确定模块803用于根据所述段数、节点数、各井段的分支数目、流体在该 井段的流入位置及流出口位置确定各井段、节点以及分支的位置,搭建多井段模型的骨架;
[0113] 第四信息确定模块8〇4用于根据单井、油藏及所述流量控制设备的静态参数确定 每井段、分支及节点的深度;
[0114] 第五信息确定模块8〇5用于根据油管直径和套管直径及粗糙度确定各井段的直 径、截面积、体积及粗糙度;
[0115] 第六信息确定模块8〇6用于根据模拟精度、井身结构和流量控制设备的动静态参 数确定最长段长度和最小段长度,使多井段模型内的各井段均短于最长段长度,均长于最 小段长度。
[0116] —实施例中,所述油藏数值模拟装置还包括:计算方法选取单元,用于根据流入、 流出口位置、井底流压、每段、分支及节点的深度、油管直径和套管直径及粗糙度确定各段 的直径、截面积、体积及粗糙度选择井筒压力计算方法。
[0117] 一实施例中,如图9所示,所述数值模型生成单元704包括:方程建立模块901,数学 模型建立模块902及数值模型生成模块903。
[0118]方程建立模块901用于根据所述多井段模型,建立各井段的各分支内的油气水三 相流动方程、压力方程;
[0119]数学模型建立模块902用于将所述多井段模型作为内边界条件与油藏渗流模型串禹 合,建立综合数学模型;
[0120]数值模型生成模块9〇3用于对所述综合数学模型进行离散化,形成数值模型。
[0121] —实施例中,所述油藏数值模拟装置还包括:流动信息计算单元,用于根据流量方 程计算某一相流体从网格块到井段节点的流动,所述流量方程如下:
[0122] qpj = TwjMpj(Pj+Hcj—Pn-Hnc)
[0123]其中,qpj是p相在网格块与井段连接j处的流量,Twj是连接传导率,Mpj是p相在网格 块与井段连接j处的流度,Pj是网格块压力,H。」是网格块到连接处的深度差,Pn是段节点处 的压力,Hn。是段节点到连接处的深度差。
[0124]本发明实施例的有益效果在于,通过本发明,实现了对采取不同井型分层注水、分 层采^工艺的油藏进行高精度数值模拟,为科学的使用井下流量控制设备提供了依据,达 到提高了采收率的目的。针对带有井下流量控制设备的直井、水平井及复杂结构井,本发明 能够模拟水嘴、封隔器、栗等多种井下流量控制设备,可适用于生产井和注入井。本发明考 虑井筒内的管流特征以及流量控制设备的压损特征计算压力分布,提高了油藏数值模拟技 术中井的模拟精度。
[0125]本领域内的g术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序 产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实 施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机 可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产 品的形式。
[0126]本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程 图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流 程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序 指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产 生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实 %在流fe图一|王或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。 [0127]这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特 定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指 令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或 多个方框中指定的功能。
[0128]这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计 算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或 其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一 个方框或多个方框中指定的功能的步骤。
[0129]本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例 的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员, 依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内 容不应理解为对本发明的限制。

Claims (14)

1. 一种基于井下流量控制设备的油藏数值模拟方法,其特征在于,所述的油藏数值模 拟方法包括: 获取生产井、注入井、油藏及流量控制设备数据,确定各井、油藏及流量控制设备的动 态参数及静态参数; 根据每口井的井身结构和流量控制设备创建多井段模型; 根据油藏特征及流体属性建立油藏渗流模型,所述油藏渗流模型包括:油气水三相渗 流方程、外边界条件以及约束条件; / 将所述多井段模型及油藏渗流模型耦合,生成数值模型; 根据油藏的静态参数建立地质模型,并根据原始地层压力分布或者压力与深度的关 系,以及岩石、流体属性对所述地质模型进行初始化; 根据所述数值模型及初始化后的地质模型建立带有稀疏系数矩阵的非线性方程组并 求解所述非线性方程组。
2.根据权利要求1所述的油藏数值模拟方法,其特征在于,根据每口井的井身结构和流 量控制设备创建多井段模型,包括: 根据井眼轨迹,流量控制设备的安装使用位置以及所需模拟精度确定所述多井段模型 的段数及节点数; 根据所述流量控制设备的使用方式确定各井段的分支数目、流体在该井段的流入位置 及流出口位置; 根据所述段数、节点数、各井段的分支数目、流体在该井段的流入位置及流出口位置确 定各井段、节点以及分支的位置,搭建多井段模型的骨架; 根据单井、油藏及所述流量控制设备的静态参数确定每井段、分支及节点的深度; 根据油管直径和套管直径及粗糙度确定各井段的直径、截面积、体积及粗糙度; 根据模拟精度、井身结构和流量控制设备动静态参数确定最长段长度和最小段长度, 使多井段模型内的各井段均短于最长段长度,均长于最小段长度。
3. 根据权利要求1所述的油藏数值模拟方法,其特征在于,在根据每口井的井身结构和 流量控制设备创建多井段模型之后,所述的油藏数值模拟方法包括:根据流入、流出口位 置、井底流压、每段、分支及节点的深度、油管直径和套管直径及粗糙度确定各段的直径、截 面积、体积及粗糙度选择井筒压力计算方法。
4. 根据权利要求1所述的油藏数值模拟方法,其特征在于,将所述多井段模型及油藏渗 流模型耦合,生成数值模型,包括: 根据所述多井段模型,建立各井段的各分支内的油气水三相流动方程、压力方程; 将所述多井段模型作为内边界条件与油藏渗流模型耦合,建立综合数学模型; 对所述综合数学模型进行离散化,形成数值模型。
5.根据权利要求4所述的油藏数值模拟方法,其特征在于,根据流量方程计算某一相流 体从网格块到井段节点的流动,所述流量方程如下: Qpj 二 TwjMpj (Pj+Hcj_Pn_Hnc) 其中,qPj是P相在网格块与井段连接j处的流量,Twj是连接传导率,MPj是p相在网格块与 井段连接j处的流度,Pj是网格块压力,Hcj是网格块到连接处的深度差,?„是段节点处的压 力,H„。是段节点到连接处的深度差。
6.根据权利要求1-5中任一项所述的油藏数值模拟方法,其特征在于,所述静态参数包 括:每口井的井位坐标、完井方式、射孔深度、井头、井深、井斜及井轨迹数据,油藏目的层位 的构造形态、测井解释数据、地震解释数据、断层数据,目的层位的孔隙度、渗透率、含油饱 和度及原始地层压力参数,地层内岩石及流体的高压物性资料,流量控制设备的类型、安装 位置、设备尺寸及控制方式。
7.根据权利要求1-5中任一项所述的油藏数值模拟方法,其特征在于,所述动态参数包 括:单井生产动态数据、吸水剖面测试数据、产油剖面测试数据、井下流量控制设备的调节 参数及井下流量控制设备监测数据。
8.—种基于井下流量控制设备的油藏数值模拟装置,其特征在于,所述的油藏数值模 拟装置包括: 参数获取单元,用于获取生产井、注入井、油藏及流量控制设备数据,确定各井、油藏及 流量控制设备的动态参数及静态参数; 多井段模型创建单元,用于根据每口井的井身结构和流量控制设备创建多井段模型; 油藏渗流模型生成单元,用于根据油藏特征及流体属性建立油藏渗流模型,所述油藏 渗流模型包括:油气水三相渗流方程、外边界条件以及约束条件; 数值模型生成单元,用于将所述多井段模型及油藏渗流模型耦合,生成数值模型; 地质模型生成单元,用于根据油藏的静态参数建立地质模型,并根据原始地层压力分 布或者压力与深度的关系,以及岩石、流体属性对所述地质模型进行初始化; 稀疏矩阵生成单元,用于根据所述数值模型及初始化后的地质模型建立带有稀疏系数 矩阵的非线性方程组并求解所述非线性方程组。
9.根据权利要求8所述的油藏数值模拟装置,其特征在于,所述多井段模型创建单元包 括: 第一信息确定模块,用于根据井眼轨迹,流量控制设备的安装使用位置以及所需模拟 精度确定所述多井段模型的段数及节点数; 第二信息确定模块,用于根据所述流量控制设备的使用方式确定各井段的分支数目、 流体在该井段的流入位置及流出口位置; 第三信息确定模块,用于根据所述段数、节点数、各井段的分支数目、流体在该井段的 流入位置及流出口位置确定各井段、节点以及分支的位置,搭建多井段模型的骨架; 第四信息确定模块,用于根据单井、油藏及所述流量控制设备的静态参数确定每井段、 分支及节点的深度; 第五信息确定模块,用于根据油管直径和套管直径及粗糙度确定各井段的直径、截面 积、体积及粗糙度; 第六信息确定模块,用于根据模拟精度、井身结构和流量控制设备的动静态参数确定 最长段长度和最小段长度,使多井段模型内的各井段均短于最长段长度,均长于最小段长 度。
10.根据权利要求8所述的油藏数值模拟装置,其特征在于,所述油藏数值模拟装置还 包括:计算方法选取单元,用于根据流入、流出口位置、井底流压、每段、分支及节点的深度、 油管直径和套管直径及粗糙度确定各段的直径、截面积、体积及粗糙度选择井筒压力计算 方法。
11. 根据权利要求8所述的油藏数值模拟装置,其特征在于,所述数值模型生成单元包 括: 方程建立模块,用于根据所述多井段模型,建立各井段的各分支内的油气水三相流动 方程、压力方程; 数学模型建立模块,用于将所述多井段模型作为内边界条件与油藏渗流模型耦合,建 立综合数学模型; 数值模型生成模块,用于对所述综合数学模型进行离散化,形成数值模型。
12. 根据权利要求11所述的油藏数值模拟装置,其特征在于,所述油藏数值模拟装置还 包括:流动信息计算单元,用于根据流量方程计算某一相流体从网格块到井段节点的流动, 所述流量方程如下: qpj = TwjMpj (Pj+Hcj—Pn—Hnc) 其中,qpj是P相在网格块与井段连接j处的流量,Twj是连接传导率,Mpj是p相在网格块与 井段连接j处的流度,Pj是网格块压力,是网格块到连接处的深度差,段节点处的压 力,Hn。是段节点到连接处的深度差。 I3 •根据权利要求8-12中任一项所述的油藏数值模拟装置,其特征在于,所述静态参数 包括:每口井的井位坐标、完井方式、射孔深度、井头、井深、井斜及井轨迹数据,油藏目的层 位的构造形态、测井解释数据、地震解释数据、断层数据,目的层位的孔隙度、渗透率、含油 饱和度及原始地层压力参数,地层内岩石及流体的高压物性资料,流量控制设备的类型、安 装位置、设备尺寸及控制方式。
14.根据权利要求8-12中任一项所述的油藏数值模拟装置,其特征在于,所述动态参数 包括:单井生产动态数据、吸水剖面测试数据、产油剖面测试数据、井下流量控制设备的调 节参数及井下流量控制设备监测数据。
CN201510284080.2A 2015-05-28 2015-05-28 基于井下流量控制设备的油藏数值模拟方法及装置 Active CN105089582B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510284080.2A CN105089582B (zh) 2015-05-28 2015-05-28 基于井下流量控制设备的油藏数值模拟方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510284080.2A CN105089582B (zh) 2015-05-28 2015-05-28 基于井下流量控制设备的油藏数值模拟方法及装置

Publications (2)

Publication Number Publication Date
CN105089582A CN105089582A (zh) 2015-11-25
CN105089582B true CN105089582B (zh) 2017-12-05

Family

ID=54570815

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510284080.2A Active CN105089582B (zh) 2015-05-28 2015-05-28 基于井下流量控制设备的油藏数值模拟方法及装置

Country Status (1)

Country Link
CN (1) CN105089582B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106223939A (zh) * 2016-07-19 2016-12-14 中国石油天然气股份有限公司 一种油藏数值确定方法和装置
CN109891047A (zh) * 2016-09-02 2019-06-14 沙特阿拉伯石油公司 控制碳氢化合物产出
CN108242069B (zh) * 2016-12-23 2021-03-23 中国石油化工股份有限公司 基于地震剖面及时深关系的油藏剖面编绘方法
CN107066718B (zh) * 2017-04-05 2020-05-15 西南石油大学 一种四维动态地应力模拟方法
CN106991236B (zh) * 2017-04-05 2020-05-19 西南石油大学 一种基于四维地应力动态变化的重复压裂选井选层方法
CN107133452B (zh) * 2017-04-18 2019-12-03 中国石油大学(北京) 油藏渗流数值模拟方法及装置
CN107387030B (zh) * 2017-07-26 2020-04-10 中国海洋石油总公司 一种利用本井气的接力举升采油设计方法
CN107832515B (zh) * 2017-11-01 2020-08-28 中国石油大学(北京) 油藏与井筒的耦合模拟方法和装置
CN109142676B (zh) * 2018-08-03 2021-05-18 中国石油集团川庆钻探工程有限公司 一种利用测井渗透率进行地层有效裂缝的识别方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102203638A (zh) * 2008-09-19 2011-09-28 雪佛龙美国公司 用于模拟地质力学储层系统的计算机实现的系统和方法
CN104179499A (zh) * 2013-05-27 2014-12-03 中国石油化工股份有限公司 考虑油藏参数时变的数值模拟方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102203638A (zh) * 2008-09-19 2011-09-28 雪佛龙美国公司 用于模拟地质力学储层系统的计算机实现的系统和方法
CN104179499A (zh) * 2013-05-27 2014-12-03 中国石油化工股份有限公司 考虑油藏参数时变的数值模拟方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于有限元法的油藏开发数值模拟;杨军征等;《新疆石油地质》;20110228;第32卷(第1期);54-55 *
非结构网格油藏数值模拟方法研究;谢海兵等;《石油学报》;20010131;第22卷(第1期);63-65 *

Also Published As

Publication number Publication date
CN105089582A (zh) 2015-11-25

Similar Documents

Publication Publication Date Title
CN105089582B (zh) 基于井下流量控制设备的油藏数值模拟方法及装置
CN105840187B (zh) 致密性油藏水平井分段压裂产能计算方法
CN104533370B (zh) 压裂水平井油藏、裂缝、井筒全耦合模拟方法
CN104879103B (zh) 一种分层注水效果分析方法
CN103590824B (zh) 经过多段压裂改造后的致密气藏水平井的产能计算方法
Yuan et al. Residual oil distribution characteristic of fractured-cavity carbonate reservoir after water flooding and enhanced oil recovery by N2 flooding of fractured-cavity carbonate reservoir
Hustedt et al. Induced fracturing in reservoir simulations: Application of a new coupled simulator to a waterflooding field example
CN105089595B (zh) 水平压裂裂缝导流作用下的油藏数值模拟方法及装置
Amini et al. Using distributed volumetric sources to predict production from multiple-fractured horizontal wells under non-Darcy-flow conditions
CN104727811A (zh) 鱼骨状水平井分段耦合的产能预测方法
CN103510950A (zh) 一种复杂结构井地层流动与管内流动耦合流动实验系统
Yu et al. Non-Darcy flow numerical simulation of XPJ low permeability reservoir
CN110206522A (zh) 一种页岩气藏压裂水平井压裂液返排模拟方法
Jiang et al. A generic physics-based numerical platform with hybrid fracture modelling techniques for simulating unconventional gas reservoirs
Luo et al. Numerical simulation of the impact of polymer rheology on polymer injectivity using a multilevel local grid refinement method
CN107832540A (zh) 一种致密油藏技术极限井距确定方法
CN103161436B (zh) 一种稠油热采水平井试井解释方法
Pan et al. How to sustain a CO2-thermosiphon in a partially saturated geothermal reservoir: Lessons learned from field experiment and numerical modeling
CN104847314B (zh) 高温高压油气直井单相流射孔完井参数优化方法
Liu et al. Manufacturing method of large-scale fractured porous media for experimental reservoir simulation
CN106469333B (zh) 一种低渗稠油油藏水平井热采压力分布预测方法
CN107169227B (zh) 一种分段压裂水平井的粗网格模拟方法及系统
Lei et al. A dynamic discrete fracture model for fluid flow in fractured low-permeability reservoirs
CN104989385A (zh) 基于表皮系数计算的高温高压油气直井射孔参数优化方法
Zhou et al. Study on Percolation Law in Low Permeability Reservoir of Deformable Media

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
GR01 Patent grant