CN113445972A - 一种致密油藏分层注水测调周期确定方法及系统 - Google Patents

一种致密油藏分层注水测调周期确定方法及系统 Download PDF

Info

Publication number
CN113445972A
CN113445972A CN202110932795.XA CN202110932795A CN113445972A CN 113445972 A CN113445972 A CN 113445972A CN 202110932795 A CN202110932795 A CN 202110932795A CN 113445972 A CN113445972 A CN 113445972A
Authority
CN
China
Prior art keywords
injection
layer
water
water injection
pressure
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
Application number
CN202110932795.XA
Other languages
English (en)
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.)
Xian Shiyou University
Original Assignee
Xian Shiyou University
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 Xian Shiyou University filed Critical Xian Shiyou University
Priority to CN202110932795.XA priority Critical patent/CN113445972A/zh
Publication of CN113445972A publication Critical patent/CN113445972A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
    • E21B43/20—Displacing by water
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B21/00—Methods or apparatus for flushing boreholes, e.g. by use of exhaust air from motor
    • E21B21/08—Controlling or monitoring pressure or flow of drilling fluid, e.g. automatic filling of boreholes, automatic control of bottom pressure
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00—Survey of boreholes or wells
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00—Survey of boreholes or wells
    • E21B47/06—Measuring temperature or pressure
    • 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
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Forestry; Mining

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Business, Economics & Management (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Agronomy & Crop Science (AREA)
  • Animal Husbandry (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Geometry (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Mechanical Engineering (AREA)
  • Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)

Abstract

一种致密油藏分层注水测调周期确定方法及系统，包括以下步骤：建立分层注水多层节点分析模型；对分层注水多层节点分析模型用非线性方程组数值求解方法中的逐步搜索法进行求解；通过分层注水多层节点分析模型，得到此工况下的两层注入压力和实际注入量，再根据单井配注误差的定义，设计配注量‑实际配注量/设计配注量，预测出配注误差；根据地层压力随时间的下降趋势，代入分层注水多层节点分析模型，得到两层注入压力和实际注入量随时间的变化趋势，将其与配注量要求对比，预测出配注测调周期。本发明对明确低渗透油藏分层注水特征与规律，以及提高分注工艺效果提供了相应的理论依据，可以进行分层注水分注配注误差、配注合格率、测调周期的预测，有效的指导分层注水工艺实践，提搞注水效益。

Description

一种致密油藏分层注水测调周期确定方法及系统
技术领域
本发明属于石油工程技术领域，特别涉及一种致密油藏分层注水测调周期确定方法及系统。
背景技术
目前国内分层注水技术历经固定式分注、钢丝投捞式分注、电缆测调式分注，进入智能分注的第四个阶段，即在装备和工具已经智能化的发展趋势下，如何针对注水开发实际高效用好此类工具，与之配套的分层注水注采井的压力系统理论研究至为关键，其首要核心问题是分注配注量的确定。这也上最新的利用大数据及人工智能分析进行精细注水优化研究中，多层多向注入量和产量劈分技术涉及到的内容。
分层注水配注量的确定方法有基于静态参数的单因素有效厚度法(H)及地层系数法(KH)、基于生产动态的多因素的劈分系数法，基于驱替表征的方法，和分层注水压力系统分析方法。其中基于压力系统分析的方法又有分注节点分析，完全油藏数值模拟，和管流与渗流耦合的方法。节点分析法基于流出和流入节点的质量、动量守恒原理，与其它方法相比，简洁实用。国内外对节点分析在注水工艺及分层注水工艺中的应用研究还较少。
致密油藏储层注水开发特征鲜明，一是非均质性强，层间和层内吸水差异性大；二是单井和单层注水量小，分注人工或井下自动测试调配难度大。为了进一步提高其分注开发效果，需进行分注配注误差和分注测调周期预测的研究。
发明内容
本发明的目的是提供一种致密油藏分层注水测调周期确定方法及系统，解决了现有技术中存在的问题。
为实现上述目的，本发明采用以下技术方案：
一种致密油藏分层注水测调周期确定方法，包括以下步骤：
在单层笼统注水节点模型的基础上建立分层注水多层节点分析模型；
对分层注水多层节点分析模型进行求解，直至井底节点处流入和流出压力差满足第一层井底节点流入和流出压力差要求最大值为止；
通过分层注水多层节点分析模型，得到此工况下的两层注入压力和实际注入量，再根据单井配注误差的定义，设计配注量-实际配注量/设计配注量，预测出配注误差；
根据地层压力随时间的下降趋势，代入分层注水多层节点分析模型，得到两层注入压力和实际注入量随时间的变化趋势，将其与配注量要求对比，预测出配注测调周期。
进一步的，单层笼统注水井节点分析中，压降系统的划分包括：地面水平管流压降△p1，地下垂直管流压降△p2，嘴流压降△p3和注入水在地下的渗流压降△p4四部分；其中，管流压降用达西公式计算，嘴流压降用嘴流压降模型计算，注入水在地下的渗流压降用平面径向流压降模型计算；分层注水井节点模型的压降系统划分的区别在于井筒管流压降、嘴流压降、渗流压降都按多层各分为多部分，分别为△p2-1、…、△p2-m，△p3-1、…、△p3-m，△p4-1、△p4-m，m为分层注水的层数；
进一步的，各部分压降的具体算式如下：
Figure BDA0003211728030000021
式中，Q为注水量,m3/d；Lp为地面管线长度,m；Dp为地面管线内径,m；Ap为地面管线内截面积,m2；
Figure BDA0003211728030000022
式中，Lw、Hw分别为井深和垂深,m；Dw为油管内径,m；Aw为油管内截面积,m2；λp、λw分别为地面管线和油管内流动摩阻系数，无量纲；下标1、m分别指第一层和第m层，m＝1，2，…，mmax；mmax为分层注水最大层数，无量纲；
Figure BDA0003211728030000023
式中Dn为配水嘴内径，m；An为配水嘴过流面积，对可调味水嘴，如“U”型节流口，过流面积m2；Cd为可调水嘴流量系数，Cd(x)＝-3.1004x4+8.7432x3-8.8037x2+3.64x+0.2635，无量纲；x为水嘴开度，0-1.0，无量纲；
Figure BDA0003211728030000031
式中，h为储层有效厚度,m；Kw为水相渗透率,μm2；rw为井筒半径,m；re为注水波及半径,m；s为表皮因子,无量纲；ρw为水相密度,kg/m3；μw为水相粘度,mPa·s；Bw为水相体积系数,无量纲；
分层注水节点分析模型的建立，以第1层井底流压为求解点，流入动态方程如式、第1层流出动态方程、第2至mmax层流出动态方程，连同多层的注入量间的关系式，一起构成多层分注注水井的节点分析模型，如式(5)所示；
Figure BDA0003211728030000032
式中，pjwf-in为流入节点井底压力，Pa；pjwf-out为流出节点井底压力，Pa；pr为地层压力,MPa。
进一步的，采用非线性方程组数值求解方法中的逐步搜索法求解，先假设配注量劈分系数(上层注入流量与总注入流量之比)rsfj和总注入量Q，其值rsfj∈(0，1),Q∈[0,1600]，再代入上层流出动态方程计算，得到总注入量和Q、上层注入量Q1和注入压力pjwf1，接着代入下层流出动态方程求得下层注入量Q2和注入压力pjwf2，从小至大逐个试算搜索，直到井底节点处流入和流出压力差满足要求ε为止，计算模型如式(6)所示；
Figure BDA0003211728030000033
式中，Qj为总配注量，m3/d；rsf-m为第m层配注量劈分系数，无量纲；ε为第一层井底节点流入和流出压力差要求最大值，MPa。
进一步的，配注误差，表达式如下：
EQ＝[Q(t2)-Q(t1)]/Q(t1) (7)
如果两个时刻仅地层压力发生了变化，则有：
[Q，Q1，…，Qm，pjwf1，…，pjwfm](t)＝fw[(流体物性:ρw,Bw,μw)，(地面管线:Lp,εp,Dp)，(井筒1:Lw1,Hw1,εw,Dw,rw)，(井筒2:Lw2,Hw2,)，…，(井筒m:Lwm,Hwm,)，(配水嘴1:2Rn1,hn1+2Rn1,xn1)，(配水嘴2:2Rn2,hn2+2Rn2,xn2)，…，(配水嘴m:2Rnm,hnm+2Rnm,xn2)，(储层1:h1,re1,pr1(t),Kw1,s1,pff,)，(储层2:h2,re2,pr2(t),Kw2,s2)，…，(储层m:hm,rem,prm(t),Kwm,sm)，(注水工艺:pj,Qmax)]；(8)
相应配注合格率的计算式为：
RQ＝1-EQ (9)。
进一步的，测调周期的定义为，两次测调配注误差大于要求时的时间间隔，表达式如下：
CYQ＝tn-t1,当EQ＝[Q(tn)-Q(t1)]/Q(t1)＜EQ要求时,n＝2,3,4,... (10)
[Q，Q1，…，Qm，pjwf1，…，pjwfm](t)＝fz[(流体物性:ρw,Bw,μw)，(地面管线:Lp,εp,Dp)，(井筒1:Lw1,Hw1,εw,Dw,rw)，(井筒2:Lw2,Hw2,)，…，(井筒m:Lwm,Hwm,)，(配水嘴1:2Rn1,hn1+2Rn1,xn1)，(配水嘴2:2Rn2,hn2+2Rn2,xn2)，…，(配水嘴m:2Rnm,hnm+2Rnm,xn2)，(储层1:h1,re1,pr1(t),Kw1,s1,pff,)，(储层2:h2,re2,pr2(t),Kw2,s2)，…，(储层m:hm,rem,prm(t),Kwm,sm)，(注水工艺:pj,Qmax)])]；(11)。
进一步的，设计配注量为前次测调配注量；实际配注量为后次测调配注量。
进一步的，一种致密油藏分层注水测调周期确定系统，包括：
模型建立模块，用于在单层笼统注水节点模型的基础上建立分层注水多层节点分析模型；
求解模块，用于对分层注水多层节点分析模型进行求解，直至井底节点处流入和流出压力差满足第一层井底节点流入和流出压力差要求最大值为止；
配注误差预测模块，用于通过分层注水多层节点分析模型，得到此工况下的两层注入压力和实际注入量，再根据单井配注误差的定义，设计配注量-实际配注量/设计配注量，预测出配注误差；
配注测调周期预测模块，用于根据地层压力随时间的下降趋势，代入分层注水多层节点分析模型，得到两层注入压力和实际注入量随时间的变化趋势，将其与配注量要求对比，预测出配注测调周期。
与现有技术相比，本发明有以下技术效果：
本发明，考虑分层注水工艺中配水间、注水井井口、注水井底配水器、注水井储层各位置处的压力和流量关系，建立了分层注水注水井节点分析模型，可用于2层到多层的分层注水压力系统的分析，对模型的求解采用非线性方程组数值求解方法中的逐步搜索法，此法相对于传统图解法在分层层数>2层时可快速求解，理论上可进行无数层的求解，所以在工程实际计算中更具优越性。再将节点模型结合分注测调采集数据，建立了分注配注误差和分注测调周期的预测方法，该方法可以进行分层注水分注配注误差、配注合格率、测调周期的预测，这三项参数可应用于有效指导分层注水工艺实践，如：分层注水工艺的选井依据、分层注水工艺的应用界限、分层注水工艺的测调预警等。本发明对明确低渗透油藏分层注水特征与规律提供了相应的理论依据，可明显提高分层注水效益。
附图说明
图1分层注水注水井节点分析中压降系统的划分示意图；
图2分注配注周期计算结果图。
具体实施方式
以下结合附图对本发明进一步说明：
请参阅图1至图2，
(1)分层注水多层节点分析模型及其求解
节点分析模型体现的是流入和流出节点的流体质量和动量守恒，其各段的压降计算采用稳态的管流和稳态的渗流数学模型。单层笼统注水井节点分析中，压降系统的划分包括：地面水平管流压降△p1，地下垂直管流压降△p2，嘴流压降△p3，和注入水在地下的渗流压降△p4四部分。其中，管流压降用达西公式计算，嘴流压降用嘴流压降模型(分固定水嘴和可调水嘴)[10,17]计算，注入水在地下的渗流压降用平面径向流压降模型计算。分层注水井节点模型是在单层笼统注水节点模型的基础上建立的，其压降系统划分的区别在于井筒管流压降、嘴流压降、渗流压降都按多层各分为多部分，即分别为△p2-1、…、△p2-m，△p3-1、…、△p3-m，△p4-1、△p4-m，m为分层注水的层数，如图1所示。
各部分压降的具体算式如下：
Figure BDA0003211728030000061
式中，Q为注水量,m3/d；Lp为由配水间到井口的注水管线长度,m；Dp为注水管线内径,m；Ap为地面管线内截面积,m2；
Figure BDA0003211728030000062
式中，Lw、Hw分别为井深和垂深,m；Dw为油管内径,m；Aw为油管内截面积,m2；λp、
λw分别为地面管线和油管内流动摩阻系数，无量纲；下标1、m分别指第一层和第m层，m＝1，2，…，mmax；mmax为分层注水最大层数，无量纲。
Figure BDA0003211728030000063
式中Dn为配水嘴内径，m；An为配水嘴过流面积，对可调味水嘴，如“U”型节流口，过流面积，m2；Cd为可调水嘴流量系数，Cd(x)＝-3.1004x4+8.7432x3-8.8037x2+3.64x+0.2635，无量纲；x为水嘴开度，0-1.0，无量纲；
Figure BDA0003211728030000064
式中，h为储层有效厚度,m；Kw为水相渗透率,μm2；rw为井筒半径,m；re为注水波及半径,m；s为表皮因子,无量纲；ρw为水相密度,kg/m3；μw为注入水的粘度,mPa·s；Bw为注入水的体积系数,无量纲；
分层注水节点分析模型的建立，以第1层井底流压为求解点，流入(井筒内管流)动态方程如式、第1层流出(储层中渗流)动态方程、第2至mmax层流出(储层中渗流)动态方程，连同多层的注入量间的关系式，一起构成多层分注注水井的节点分析模型，如式(5)所示。
Figure BDA0003211728030000071
式中，pjwf-in为流入节点井底压力，Pa；pjwf-out为流出节点井底压力，Pa；pr为地层压力,MPa；。
对节点模型的求解,传统方法采用图解法，即绘出某节点处的流入和流出动态曲线，两线交点对应的压力和流量即为满足系统平衡时的注入压力和注水量。在双层注水节点分析模型中，井底节点处为一条流入两条流出，故需考虑两条流出动态曲线的叠加，可又无法确定上下层注入量，所以图解法在分层注水时节点模型中实际无法求解。
对此，本专利对节点分析模型的求解采用非线性方程组数值求解方法中的逐步搜索法。原理是先假设配注量劈分系数(上层注入流量与总注入流量之比)rsfj和总注入量Q，其值rsfj∈(0，1),Q∈[0,1600]，再代入上层流出动态方程计算，得到总注入量和Q、上层注入量Q1和注入压力pjwf1，接着代入下层流出动态方程求得下层注入量Q2和注入压力pjwf2，从小至大逐个试算搜索，直到井底节点处流入和流出压力差满足要求ε为止，计算模型如式(6)所示。
Figure BDA0003211728030000072
式中，Qj为总配注量，m3/d；rsf-m为第m层配注量劈分系数，无量纲；ε为第一层井底节点流入和流出压力差要求最大值，MPa。对节点分析模型采用非线性方程组数值求解方法中的逐步搜索法，更简便快捷，且在>2层的多层分注时更具有优越性。
(2)分层注水配注误差和配注合格率的预测
在储层参数和注水工艺参数确定的情况下，通过分注注水井节点分析模型，可计算得到此工况下的两层注入压力和实际注入量。再根据单井配注误差的定义，[设计配注量(前次测调配注量)-实际配注量(后次测调配注量)]/设计配注量(前次测调配注量)，即可预测出配注误差，表达式如下：
EQ＝[Q(t2)-Q(t1)]/Q(t1) (7)
如果两个时刻仅地层压力发生了变化，则有：
[Q，Q1，…，Qm，pjwf1，…，pjwfm](t)＝fw[(流体物性:ρw,Bw,μw)，(地面管线:Lp,εp,Dp)，(井筒1:Lw1,Hw1,εw,Dw,rw)，(井筒2:Lw2,Hw2,)，…，(井筒m:Lwm,Hwm,)，(配水嘴1:2Rn1,hn1+2Rn1,xn1)，(配水嘴2:2Rn2,hn2+2Rn2,xn2)，…，(配水嘴m:2Rnm,hnm+2Rnm,xn2)，(储层1:h1,re1,pr1(t),Kw1,s1,pff,)，(储层2:h2,re2,pr2(t),Kw2,s2)，…，(储层m:hm,rem,prm(t),Kwm,sm)，(注水工艺:pj,Qmax)]. (8)
相应配注合格率的计算式为：
RQ＝1-EQ (9)
(3)分层注水测调周期的预测
根据地层压力随时间的下降趋势，代入分注注水井节点分析模型，可计算得到两层注入压力和实际注入量随时间的变化趋势，将其与配注量要求对比，就可预测出合适的配注测调周期。测调周期的定义为，两次测调配注误差大于要求时的时间间隔，表达式如下：
CYQ＝tn-t1,当EQ＝[Q(tn)-Q(t1)]/Q(t1)＜EQ要求时,n＝2,3,4,... (10)
[Q，Q1，…，Qm，pjwf1，…，pjwfm](t)＝fz[(流体物性:ρw,Bw,μw)，(地面管线:Lp,εp,Dp)，(井筒1:Lw1,Hw1,εw,Dw,rw)，(井筒2:Lw2,Hw2,)，…，(井筒m:Lwm,Hwm,)，(配水嘴1:2Rn1,hn1+2Rn1,xn1)，(配水嘴2:2Rn2,hn2+2Rn2,xn2)，…，(配水嘴m:2Rnm,hnm+2Rnm,xn2)，(储层1:h1,re1,pr1(t),Kw1,s1,pff,)，(储层2:h2,re2,pr2(t),Kw2,s2)，…，(储层m:hm,rem,prm(t),Kwm,sm)，(注水工艺:pj,Qmax)])]. (11)
分注测调周期预测的基础上还可进行测调井分类及其测调预警，方法为：将研究区块注水井的储层参数，代入分注注水井节点分析模型，可求出各井的注入压力和配注量，将其与区块的设计配注量对比，即可得到哪些井无需测调哪些井需要测调并且在何时测调，由此就可对分注井进行分类并进行测调预警。
实例1分层注水配注误差、配注合格率的预测
进行实例计算，输入的已知参数取自长庆油田G83区L42-20井的资料，包括《L42-20井压力测试成果报告》、《L42-20井分注调配测试成果报告》和《G83区长4+5长6油藏基础数据》，具体如表1所示，其中两层水嘴开度全设为1.0，注水井地层压力随时间能增加。
表1输入的已知参数
Figure BDA0003211728030000091
Figure BDA0003211728030000101
分注配注误差计算结果如表2所示：
表2分注配注误差计算结果
Figure BDA0003211728030000102
实例2分层注水测调周期的预测
进行实例计算，输入的已知参数如表3所示，其中两层水嘴开度全设为1.0，假设地层压力与时间的拟合关系为pr＝pr0+0.01t，EQ要求＝50％。
表3输入的已知参数
Figure BDA0003211728030000103
Figure BDA0003211728030000111
分注配注周期计算结果如表4和附图2所示：
表4分注配注周期计算结果
Figure BDA0003211728030000112
可见，计算得到此时该井的测调周期为40～50天。

Claims (8)

1.一种致密油藏分层注水测调周期确定方法，其特征在于，包括以下步骤：
在单层笼统注水节点模型的基础上建立分层注水多层节点分析模型；
对分层注水多层节点分析模型进行求解，直至井底节点处流入和流出压力差满足第一层井底节点流入和流出压力差要求最大值为止；
通过分层注水多层节点分析模型，得到此工况下的两层注入压力和实际注入量，再根据单井配注误差的定义，设计配注量-实际配注量/设计配注量，预测出配注误差；
根据地层压力随时间的下降趋势，代入分层注水多层节点分析模型，得到两层注入压力和实际注入量随时间的变化趋势，将其与配注量要求对比，预测出配注测调周期。
2.根据权利要求1所述的一种致密油藏分层注水测调周期确定方法，其特征在于，单层笼统注水井节点分析中，压降系统的划分包括：地面水平管流压降△p1，地下垂直管流压降△p2，嘴流压降△p3和注入水在地下的渗流压降△p4四部分；其中，管流压降用达西公式计算，嘴流压降用嘴流压降模型计算，注入水在地下的渗流压降用平面径向流压降模型计算；分层注水井节点模型的压降系统划分的区别在于井筒管流压降、嘴流压降、渗流压降都按多层各分为多部分，分别为△p2-1、…、△p2-m，△p3-1、…、△p3-m，△p4-1、△p4-m，m为分层注水的层数。
3.根据权利要求2所述的一种致密油藏分层注水测调周期确定方法，其特征在于，各部分压降的具体算式如下：
Figure FDA0003211728020000011
式中，Q为注水量,m3/d；Lp为地面管线长度,m；Dp为地面管线内径,m；Ap为地面管线内截面积,m2；
Figure FDA0003211728020000012
式中，Lw、Hw分别为井深和垂深,m；Dw为油管内径,m；Aw为油管内截面积,m2；λp、λw分别为地面管线和油管内流动摩阻系数，无量纲；下标1、m分别指第一层和第m层，m＝1，2，…，mmax；mmax为分层注水最大层数，无量纲；
Figure FDA0003211728020000021
式中Dn为配水嘴内径，m；An为配水嘴过流面积，对可调味水嘴，如“U”型节流口，过流面积m2；Cd为可调水嘴流量系数，Cd(x)＝-3.1004x4+8.7432x3-8.8037x2+3.64x+0.2635，无量纲；x为水嘴开度，0-1.0，无量纲；
Figure FDA0003211728020000022
式中，h为储层有效厚度,m；Kw为水相渗透率,μm2；rw为井筒半径,m；re为注水波及半径,m；s为表皮因子,无量纲；ρw为水相密度,kg/m3；μw为水相粘度,mPa·s；Bw为水相体积系数,无量纲；
分层注水节点分析模型的建立，以第1层井底流压为求解点，流入动态方程如式、第1层流出动态方程、第2至mmax层流出动态方程，连同多层的注入量间的关系式，一起构成多层分注注水井的节点分析模型，如式(5)所示；
Figure FDA0003211728020000023
式中，pjwf-in为流入节点井底压力，Pa；pjwf-out为流出节点井底压力，Pa；pr为地层压力,MPa。
4.根据权利要求1所述的一种致密油藏分层注水测调周期确定方法，其特征在于，采用非线性方程组数值求解方法中的逐步搜索法求解，先假设配注量劈分系数，上层注入流量与总注入流量之比rsfj和总注入量Q，其值rsfj∈(0，1),Q∈[0,1600]，再代入上层流出动态方程计算，得到总注入量和Q、上层注入量Q1和注入压力pjwf1，接着代入下层流出动态方程求得下层注入量Q2和注入压力pjwf2，从小至大逐个搜索试算，直到井底节点处流入和流出压力差满足要求ε为止，计算模型如式(6)所示；
Figure FDA0003211728020000031
式中，Qj为总配注量，m3/d；rsf-m为第m层配注量劈分系数，无量纲；ε为第一层井底节点流入和流出压力差要求最大值，MPa。
5.根据权利要求1所述的一种致密油藏分层注水测调周期确定方法，其特征在于，配注误差，表达式如下：
EQ＝[Q(t2)-Q(t1)]/Q(t1) (7)
如果两个时刻仅地层压力发生了变化，则有：
[Q，Q1，…，Qm，pjwf1，…，pjwfm](t)＝fw[(流体物性:ρw,Bw,μw)，(地面管线:Lp,εp,Dp)，(井筒1:Lw1,Hw1,εw,Dw,rw)，(井筒2:Lw2,Hw2,)，…，(井筒m:Lwm,Hwm,)，(配水嘴1:2Rn1,hn1+2Rn1,xn1)，(配水嘴2:2Rn2,hn2+2Rn2,xn2)，…，(配水嘴m:2Rnm,hnm+2Rnm,xn2)，(储层1:h1,re1,pr1(t),Kw1,s1,pff,)，(储层2:h2,re2,pr2(t),Kw2,s2)，…，(储层m:hm,rem,prm(t),Kwm,sm)，(注水工艺:pj,Qmax)]；(8)
相应配注合格率的计算式为：
RQ＝1-EQ (9)。
6.根据权利要求1所述的一种致密油藏分层注水测调周期确定方法，其特征在于，测调周期的定义为，两次测调配注误差大于要求时的时间间隔，表达式如下：
CYQ＝tn-t1,当EQ＝[Q(tn)-Q(t1)]/Q(t1)＜EQ要求时,n＝2,3,4,... (10)
[Q，Q1，…，Qm，pjwf1，…，pjwfm](t)＝fz[(流体物性:ρw,Bw,μw)，(地面管线:Lp,εp,Dp)，(井筒1:Lw1,Hw1,εw,Dw,rw)，(井筒2:Lw2,Hw2,)，…，(井筒m:Lwm,Hwm,)，(配水嘴1:2Rn1,hn1+2Rn1,xn1)，(配水嘴2:2Rn2,hn2+2Rn2,xn2)，…，(配水嘴m:2Rnm,hnm+2Rnm,xn2)，(储层1:h1,re1,pr1(t),Kw1,s1,pff,)，(储层2:h2,re2,pr2(t),Kw2,s2)，…，(储层m:hm,rem,prm(t),Kwm,sm)，(注水工艺:pj,Qmax)])]；(11)。
7.根据权利要求1所述的一种致密油藏分层注水测调周期确定方法，其特征在于，设计配注量为前次测调配注量；实际配注量为后次测调配注量。
8.一种致密油藏分层注水测调周期确定系统，其特征在于，基于权利要求1至7任意一项所述的一种致密油藏分层注水测调周期确定方法，包括：
模型建立模块，用于在单层笼统注水节点模型的基础上建立分层注水多层节点分析模型；
求解模块，用于对分层注水多层节点分析模型进行求解，直至井底节点处流入和流出压力差满足第一层井底节点流入和流出压力差要求最大值为止；
配注误差预测模块，用于通过分层注水多层节点分析模型，得到此工况下的两层注入压力和实际注入量，再根据单井配注误差的定义，设计配注量-实际配注量/设计配注量，预测出配注误差；
配注测调周期预测模块，用于根据地层压力随时间的下降趋势，代入分层注水多层节点分析模型，得到两层注入压力和实际注入量随时间的变化趋势，将其与配注量要求对比，预测出配注测调周期。
CN202110932795.XA 2021-08-13 2021-08-13 一种致密油藏分层注水测调周期确定方法及系统 Pending CN113445972A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110932795.XA CN113445972A (zh) 2021-08-13 2021-08-13 一种致密油藏分层注水测调周期确定方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110932795.XA CN113445972A (zh) 2021-08-13 2021-08-13 一种致密油藏分层注水测调周期确定方法及系统

Publications (1)

Publication Number Publication Date
CN113445972A true CN113445972A (zh) 2021-09-28

Family

ID=77818621

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110932795.XA Pending CN113445972A (zh) 2021-08-13 2021-08-13 一种致密油藏分层注水测调周期确定方法及系统

Country Status (1)

Country Link
CN (1) CN113445972A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114718548A (zh) * 2022-03-14 2022-07-08 中国海洋石油集团有限公司 一种缆控分注井注水智能优化控制方法
CN115017808A (zh) * 2022-06-06 2022-09-06 西安石油大学 一种基于改进蝴蝶算法优化hkelm的管道冲蚀预测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN205297502U (zh) * 2015-12-09 2016-06-08 中国石油天然气股份有限公司 一种数字式分注累积流量自动调节系统
US20170234121A1 (en) * 2016-02-12 2017-08-17 Saudi Arabian Oil Company Systems and methods for transient-pressure testing of water injection wells to determine reservoir damages
CN111927410A (zh) * 2020-08-06 2020-11-13 中国石油天然气股份有限公司 一种井-藏协同分注测试调节设计方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN205297502U (zh) * 2015-12-09 2016-06-08 中国石油天然气股份有限公司 一种数字式分注累积流量自动调节系统
US20170234121A1 (en) * 2016-02-12 2017-08-17 Saudi Arabian Oil Company Systems and methods for transient-pressure testing of water injection wells to determine reservoir damages
CN111927410A (zh) * 2020-08-06 2020-11-13 中国石油天然气股份有限公司 一种井-藏协同分注测试调节设计方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
J.W.戴莱 等: "《流体动力学》", 31 October 1981, 人民教育出版社 *
宋洋等: "注水井节点系统分析数学模型及其应用", 《重庆科技学院学报(自然科学版)》 *
秦国伟等: "节点分析法在分层注采系统中的应用研究", 《西南石油大学学报(自然科学版)》 *
胡改星 等: "超低渗油藏节点分析模型研究及分层注水合格率影响因素初探", 《2018油气田勘探与开发国际会议（IFEDC 2018）论文集》 *
陈涛: "低渗透油藏分层注水节流阀研究", 《中国优秀硕士学位论文全文数据库 工程科技 1辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114718548A (zh) * 2022-03-14 2022-07-08 中国海洋石油集团有限公司 一种缆控分注井注水智能优化控制方法
CN115017808A (zh) * 2022-06-06 2022-09-06 西安石油大学 一种基于改进蝴蝶算法优化hkelm的管道冲蚀预测方法

Similar Documents

Publication Publication Date Title
CN110334431A (zh) 一种低渗透致密气藏单井控制储量计算及剩余气分析方法
CN113445972A (zh) 一种致密油藏分层注水测调周期确定方法及系统
CN103670369B (zh) 一种注采井间连通状况的判别方法及装置
CN104060985B (zh) 一种层状油藏调剖堵水堵剂进入深度测试方法及系统
CN105464644B (zh) 一种井口压力监测方法及控制系统
CN106869914B (zh) 一种油层中渗流与井筒内流动相耦合的产能预测方法
CN105205318B (zh) 确定多层多段水平裂缝采油井的总产量的方法和装置
CN110397430A (zh) 一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法
CN109614736B (zh) 一种海上稠油油田蒸汽吞吐开发定向井与直井产能倍数的确定方法
CN112855129A (zh) 海洋天然气水合物系统多气合采模拟装置及方法
CN106285621A (zh) 基于多因素决策海上油田老区新井产能的方法
CN112096355A (zh) 一种应用在分层注水中利用水嘴开度实现智能分注的方法
CN107130952B (zh) 一种三级复合石油增产方法
CN109667564A (zh) 一种海上稠油油田蒸汽吞吐开发定向井产能的确定方法
CN111119815A (zh) 一种均衡驱替的确定生产井配产比例的方法
CN112101710B (zh) 一种基于水驱前缘控制的量化注采平衡调整方法
CN105569625B (zh) 井下蒸汽流量干度分配调控装置
CN107355200B (zh) 一种纳微米颗粒分散体系改善水驱选井方法
CN111927410A (zh) 一种井-藏协同分注测试调节设计方法
CN111520124B (zh) 一种预测水平井流入剖面的方法
CN109522583B (zh) 一种用于多层油藏开发层系组合的非均质性界限确定方法
CN107783940A (zh) 一种分层采油前油层层间干扰的表征方法
CN115653570A (zh) 一种井筒-滑套-储层耦合流动流入动态预测方法和系统
CN111764879B (zh) 一种二元复合驱产剂窜流的评价方法
Wang et al. Flow simulation of a horizontal well with two types of completions in the frame of a wellbore–annulus–reservoir model

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20210928