一种一体化油藏数值模拟方法
技术领域
本发明涉及油气田开发领域油藏数值模拟应用方法改进与扩展,具体涉及一种一体化油藏数值模拟方法。
背景技术
当前油藏数值模拟技术主要以油藏压力、开采速度与油水产量等开发指标为研究对象,缺乏采油工程数据约束,油藏数值模拟结果与实际存在偏差,例如模拟出的油藏局部压力、单井生产指数等不准确;在油藏方案制定或生产调整数值模拟研究中,采油工程参数与地面系统数据没有得到应用,导致油藏方案计划产量与抽油机、管网系统生产能力不一致。
发明内容
本发明的目的在于提供一种一体化油藏数值模拟方法。
本发明的技术方案是:一种一体化油藏数值模拟方法,其特征在于,包括如下步骤:
(1)建立油藏数值模型,通过数值矩阵迭代求解,对井底流压、单井静压与含水率等参数进行逐一拟合计算,输出所有井静压与生产指数:
产量计算公式:Qp,j=TwjMp,j(Pj-Pw-Hwj),
其中,Pj代表某射孔连接静压,即为单井静压,Pw代表井底流压,Hwj代表某射孔连接静水压头,单位均为:MPa;Mp,j代表流度,单位:μm2/(Mpa·s);Twj代表某连接井传导率,单位:m3;Qp,j代表某射孔连接产量,单位:m3/d;
其中,J代表生产指数,单位m3/(MPa·d),
其中,Dx和Dy代表X和Y方向网格尺寸,Kx和Ky代表X和Y方向渗透率;S代表表皮系数;Kh代表地层系数;θ为井筒与网格连接夹角;c为单位转换系数英制单位为0.001127,米制为0.008527,实验室单位制为3.6;K为相对渗透率,无量纲;kpj是某一流体的相对渗透率;
Bp,j为体积系数,无量纲;μ为粘度,单位cP,rw为井直径;
(2)建立井筒模型,以步骤(1)的单井静压Pj与生产指数J为输入数据,
选择流体管流Beggs&Brown多相流公式进行压降与持液率计算,其中压降公式为:
其中,P为管段内流动的平均绝对压力,单位:Pa;
G为流体质量流量,计算公式G=ρmωmA,单位:kg/s;ρm代表气液混合物的密度;X为管段长度,单位:m;
Hl为截面持液率,无量纲;λ为两相混输水力摩阻系数,无量纲;
ρl为液相密度,单位:kg/m3;ρg为气相密度,单位:kg/m3;
ωm为气液混合物速度,单位:m/s;ωsg为气相折算速度,单位:m/s;
D为管内径,单位:m;θ为管段倾角,单位:rad;
g为重力加速度,单位:m/s2;
根据抽油机冲程与冲次等参数,以调整泵漏失参数为主,对井口压力与产量进行回归与模型校正;校正后,通过人工举升方程与井筒多相流方程代表此生产制度下的流体运动特征;
其中,回归与模型校正步骤如下:
1)井筒模型边界条件输入:步骤(1)计算的单井静压、采液指数与井口压力,其中,采液指数是指单位生产压差下油井的日产液量,
2)校正前井筒模型,
3)对比产液量与观测值,
4)若对比产液量与观测值在合理误差范围以内,则完成井筒后管网模型;若两者不在合理误差范围以内,则调整地层参数和抽油机滑脱系数后重复回归与模型校正步骤2)、3)、4)直至控制在合理误差范围以内;
(3)建立管网模型,根据井口压力、井口产量、地面系统温度与压力等参数对管网模型进行校正,使之能代表实际生产状态;同时,耦合井筒与管网模型,再次对井筒模型进一步校验;
其中,校验步骤如下:
1)管网模型边界条件输入:各井产液量、外输压力、加热炉出口温度与步骤(2)中计算的压降,
2)校正前管网模型,
3)对比井口压力与观测值,
4)若对比井口压力与观测值在合理误差范围以内,则校正完成;若两者不在合理误差范围以内,则调整流动相关式和摩擦系数系数后重复校验步骤2)、3)、4)直至控制在合理误差范围以内;
(4)对油藏与井筒管网模型进行一体化耦合并不断进行迭代求解,分析剩余油分布特征并提出调整生产制度建议。
进一步的技术方案,步骤(4)中代入地面、井筒对油藏的约束。
本发明的有益效果:
1、本发明通过地下油藏到地面管网,在数据上与流程上实现连续,不再是各专业模型独立运行,而是联立求解,互为约束;由于油藏能力决定了对人工举升的技术要求,人工举升则限制了单井产量;单井产量决定了对管网系统的设计要求,而管网系统也影响了油藏开发方案执行;所以本发明提供的方案能够更真实地反映出流体运移特征,从而提高对油藏与生产的认识。
2、本发明中通过将井筒模型、油藏模型、一体化耦合模型三种模型相结合,通过这种耦合方式,打破了油藏工程、采油工程与地面工程三者之间的专业壁垒,本发明充分利用多学科数据,考虑三大系统之间参数的约束性,以油藏上单井生产指数与井筒上生产指数为对象,通过迭代式相互校验,确定出一个更合理的生产指数;同时考虑地面管网压力与井口压力条件,利用泵漏失量与气体滑脱等参数校正井筒模型,确定单井产量与生产制度匹配性;再将油藏数值模型与井筒管网模型相耦合计算,充分考虑了油藏流体供给能力、抽油机系统与管网系统约束,从而量化了从地下油藏到地面管网的一系列生产指标,给出最优抽油机冲程、冲次等,提高了生产工艺制度调整的科学性;强化了对管网系统外输能力的认识,为老区改造中管网更换、加热炉与增压泵设计提供了更有效的论证。
附图说明
图1为生产井流入与流出动态交互图,
图2为井筒模型校正流程图,
图3井筒模型校正结果图,
图4管网模型校正流程图,
图5管网模型校正结果图,
图6为常规油藏数值模拟结果图,
图7为一体化耦合模拟结果。
具体实施方式
下面通过非限制性实施例,进一步阐述本发明,理解本发明。
本发明为一种一体化油藏数值模拟方法,包括如下步骤:
(1)建立油藏数值模型,通过数值矩阵迭代求解,对井底流压、单井静压与含水率等参数进行逐一拟合计算,输出所有单井静压与生产指数:
井名 |
油藏单井静压(MPa) |
液相生产指数(sm<sup>3</sup>/d/bar) |
NH201 |
8.8277 |
0.038777 |
NH2-1 |
9.039 |
1.88 |
NH2-6 |
12.027 |
0.094043 |
NH2-7 |
8.197 |
0.245831 |
NH2-8 |
8.754 |
3.857953 |
NH2-2 |
7.211 |
0.1 |
NH2-14 |
8.376 |
80.75948 |
NH2-13 |
7.909 |
0.722834 |
NH2-12 |
8.848 |
3.35 |
NH6 |
9.961 |
0.353728 |
NH2-17 |
83 |
2 |
NH2-9 |
8.117 |
0.07 |
NH2-20 |
7.346 |
1.12 |
NH2-P1 |
8.37 |
2.83 |
NH2-29 |
10.91 |
3.18 |
NH2-19 |
10.686 |
0.8 |
NH2-21 |
10.39 |
4.71 |
NH2-23 |
10.047 |
0.36 |
NH2-10 |
9.4021 |
7.61 |
NH2-30 |
10.69 |
4.28 |
NH2-41 |
12.706 |
7.22 |
NH2-42 |
12.85 |
6.3 |
NH2-43 |
14.041 |
1.003846 |
NH2-45 |
14.331 |
0.53 |
(2)建立井筒模型,以步骤(1)计算得到的单井静压Pj与生产指数J为输入数据,选择流体管流Beggs&Brown多相流公式进行压降计算,其中压降公式为:
其中,P为管段内流动的平均绝对压力,单位:Pa;
G为流体质量流量,计算公式G=ρmωmA,单位:kg/s;ρm代表气液混合物的密度;X为管段长度,单位:m;
Hl为截面持液率,无量纲;λ为两相混输水力摩阻系数,无量纲;
ρl为液相密度,单位:kg/m3;ρg为气相密度,单位:kg/m3;
ωm为气液混合物速度,单位:m/s;ωsg为气相折算速度,单位:m/s;
D为管内径,单位:m;θ为管段倾角,单位:rad;
g为重力加速度,单位:m/s2;
根据抽油机冲程与冲次等参数,以调整泵漏失参数为主,如图1,对井口压力与产量进行回归与模型校正;
如图2,回归与模型校正步骤如下:
1)井筒模型边界条件输入:步骤(1)计算的单井静压、采液指数与井口压力,其中,采液指数是指单位生产压差下油井的日产液量,
2)校正前井筒模型,
3)对比产液量与观测值,
4)若对比产液量与观测值在合理误差范围以内,则完成井筒后管网模型;若两者不在合理误差范围以内,则调整地层参数和抽油机滑脱系数后重复模型校正步骤2)、3)、4)直至控制在合理误差范围以内;
井筒模型校正结果如图3,当井筒内泵漏失系数越大,单井产量就越低,在同一深度处,井底流压就越高;
(3)如图4,建立管网模型,根据井口压力、井口产量、地面系统温度与压力等参数对管网模型进行校正,使之能代表实际生产状态;同时,耦合井筒与管网模型,再次对井筒模型进一步校验;
其中,校验步骤如下:
1)管网模型边界条件输入:各井产液量、外输压力、加热炉出口温度与步骤(2)中计算的压降,
2)校正前管网模型,
3)对比井口压力与观测值,
4)若对比井口压力与观测值在合理误差范围以内,则校正完成;若两者不在合理误差范围以内,则调整流动相关式和摩擦系数系数后重复校验步骤2)、3)、4)直至控制在合理误差范围以内;
主要拟合结果如图5:不同的流动方程计算出的压力结果截然不同,其中Bregg&Brill Revised更能准确地描述管网压力分布特征;
(4)对油藏与井筒管网模型进行一体化耦合并不断进行迭代求解,代入地面、井筒对油藏的约束,分析剩余油分布特征并提出调整生产制度建议。
实施例1:以某采油厂某新投入开发区块为例,根据以上步骤,进行了油藏与井筒管网一体化模拟研究,分别建立了油藏模型、井筒模型与管网模型,并在2018年3月20日生产时间点进行了一体化模拟运算与生产分析。
利用一体化耦合模型,精确预测产量变化,提出调整措施和方向。
井筒模型、油藏模型、一体化耦合模型运算结果特点如下表:
将一体化模拟结果与单油藏模型数值模拟结果进行了对比,发现两者具有明显不同的产量变化趋势,常规油藏数值模拟过程更乐观的估计了稳产时间。
图6为常规油藏数值模拟结果,即模拟在不考虑井筒与管网条件下,井日产油量,图6中,地层压力出现下降,但本年度产量仍可维持稳定,该方法建议维持当前液量;
图7为一体化耦合模拟结果,即模拟在考虑井筒与管网约束条件下,井日产油量,图7中,地层压力下降,本年度下半年产量将出现下滑,因此采用该拟合方法建议降低产液量生产。
根据分析结果,矿场根据一体化耦合模型优化建议,积极加快水井投注工作,并主动对8口油井降低冲程、冲次生产,保持合理采油速度生产,实现区块可持续稳产。
实施例2:模拟显示某地区北部管网沿程压力相对较高、集输能量充足,南部沿程压力普遍偏低,需要优选合适站点,部署增压泵,提升集输能力。根据模拟建议在管网南部新增部署了增压泵2台,使集输生产平稳运行。
分析结果指导建立该井区地面管网模型,能够精确预测沿程压力和温度,实现了由信息化管网到智能化管网的升级。
综上所述,本发明油藏与井筒、管网一体化耦合模拟三大系统之间参数的约束性,以油藏上单井生产指数与井筒上生产指数为对象,通过迭代式相互校验,力求确定出一个更合理的生产指数;同时考虑地面管网压力与井口压力条件,利用泵漏失量与气体滑脱等参数校正井筒模型,确定单井产量与生产制度匹配性;再将油藏数值模型与井筒管网模型相耦合计算,这样充分考虑了油藏流体供给能力、抽油机系统与管网系统约束,从而量化了从地下油藏到地面管网的一系列生产指标,本发明模拟方法能够更真实地反映出流体运动特征,从而提高对油藏与生产的认识。