en
×

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。
作者简介:

侯玉培(1984—),女,山东肥城人,助理研究员,硕士,从事油气藏数值模拟方法研究。E-mail:houyupei.slyt@sinopec.com。

中图分类号:TE319

文献标识码:A

文章编号:1009-9603(2021)05-0124-07

DOI:10.13673/j.cnki.cn37-1359/te.2021.05.015

参考文献 1
梁文福.油田开发智能应用系统建设成果及展望[J].大庆石油地质与开发,2019,38(5):283-289.LIANG Wenfu.Constructed achievements and prospects of the in⁃ telligent application system for the oilfield development[J].Petro⁃ leum Geology & Oilfield Development in Daqing,2019,38(5):283-289.
参考文献 2
杜金虎,时付更,杨剑锋,等.中国石油上游业务信息化建设总体蓝图[J].中国石油勘探,2020,25(5):1-8.DU Jinhu,SHI Fugeng,YANG Jianfeng,et al.Overall blueprint of information construction of PetroChina upstream business[J].Chi⁃ na Petroleum Exploration,2020,25(5):1-8.
参考文献 3
SUTTON R P,FARSHAD F.Evaluation of empirically derived PVT properties for Gulf of Mexico crude oils[C].SPE 13172,1990.
参考文献 4
BEAL C.The viscosity of air,water,natural gas,crude oil and its associated gases at oil field temperatures and pressures[C].SPE 946094,1946.
参考文献 5
CARLAW H S,JEGER J C.Conduction of heat in solids[M].Ox⁃ ford:Clarendon Press,1959:29-111.
参考文献 6
秦积舜,李爱芬.油层物理学[M].东营:石油大学出版社,2003:76-80.QIN Jishun,LI Aifen.Petrophysics[M].Dongying:Press of Univer⁃ sity of Petroleum,China,2003:76-80.
参考文献 7
姜瑞忠,刘秀伟,崔永正,等.非稳态窜流多段压裂水平井井底压力分析[J].油气地质与采收率,2019,26(5):86-95.JIANG Ruizhong,LIU Xiuwei,CUI Yongzheng,et al.Bottomhole pressure analysis of multistage fractured horizontal well during unsteady crossflow[J].Petroleum Geology and Recovery Efficien⁃ cy,2019,26(5):86-95.
参考文献 8
BEGGS H D,ROBINSON J R.Estimating the viscosity of crude oil systems[C].SPE 5434,1975.
参考文献 9
纪禄军,赵金洲,杨承孝,等.IPR方程通式中Vogel参数与流动效率R的计算[J].河南石油,2000,(4):13-15,60.JI Lujun,ZHAO Jinzhou,YANG Chengxiao.Calculation of Vogel parameter and flow efficiency R in IPR general formula[J].Henan Petroleum,2000,(4):13-15,60.
参考文献 10
VOGEL J V.Inflow performance relationships for solution gas drive wells[C].SPE 1476,1968.
参考文献 11
张琪.采油工程原理与设计[M].东营:石油大学出版社,2000:1-15.ZHANG Qi.Principle and design of oil production engineering [M].Dongying:Press of University of Petroleum,China,2000:1-15.
参考文献 12
FETKOVICH M J.The isochronal testing of oil wells[C].SPE 4529,1973.
参考文献 13
JONES L G,BLOUNT E M.Use of short term multiple rate flow tests to predict performance of wells having turbulence[C].SPE 6133,1976.
参考文献 14
CINCO LEY H,SAMANIEGO F V,DOMINGUEZ N A.Transient pressure behaviour for a well with a finite-conductivity vertical fracture[C].SPE 6014,1976.
参考文献 15
CINCO L H,RAMEY H J,MILLER F G.Pseudo-skin factors for partially-penetrating directionally-drilled wells[C].SPE 5589,1975.
参考文献 16
吴星晔.储层-井筒-井口一体化动态分析方法研究[D].成都:西南石油大学,2017.WU Xingye.The dynamic analysis method of reservoir-wellboewellhead integration[D].Chengdu:Southwest Petroleum Universi⁃ ty,2017.
参考文献 17
WASP E J,KENNY J P,GANDHI R L.Solid-liquid flow slurry pipe line transportation[M].Houston:Gulf Publishing Company,1979:25-122.
参考文献 18
SVEDEMAN S J.Criteria for sizing multiphase flowlines for ero⁃ sive/corrosive service[C].SPE 26569,1994.
参考文献 19
李学军.气井冲蚀腐蚀临界速度的计算方法[J].石油石化节能,1991,(4):94.LI Xuejun.Calculation method of critical velocity of erosion and corrosion in gas well[J].Energy Conservation in Petroleum and Petrochemical Industry,1991,(4):94.
参考文献 20
TURNER R G.Analysis and prediction of minimum flow rate for the continuous removal of liquids from gas well[C].SPE 2198,1969.
目录contents

    摘要

    油气田开发生产系统包含油藏、井筒、管网三大生产环节。传统数值模拟方法主要模拟油藏中流体的流动,而油藏-井筒-管网一体化耦合模拟同时反映油藏供液能力、井筒举升能力和管网集输能力,对准确预测生产态势、实现生产系统全流程优化具有重要意义。通过研究井筒-管网模型建立方法及流体多相管流相关式适应性,实现了井筒、管网内流动模拟和流动保障模拟;通过节点链接建立一体化模型,耦合关键节点的流动关系,形成了油藏-井筒-管网一体化耦合模拟方法,实现了油气田开发全流程模拟。基于埕岛油田西A区块,根据海上油田实际生产需求和约束条件,优化生产制度、配产配注等,基于一体化模型进行全流程优化,最大程度提高了油气田开发系统产量效益。

    Abstract

    The reservoir,wellbore and pipe network are major links in a development and production system of oil and gas fields. Traditional numerical simulation methods mainly simulate the fluid flow in the reservoir,while the coupling simula- tion based on reservoir-wellbore-pipe network integration reflects fluid supply capacity of the reservoir,lifting capacity of the wellbore,and gathering and transportation capacity of the pipe network at the same time. It is of great significance to ac- curate prediction of the production situation and the whole process optimization of a production system. By studying the con- struction method of a wellbore-pipe network model and the adaptability of multiphase pipe flow correlations,we realized the flow simulation and flow assurance simulation in the wellbore and the pipe network. An integrated model was built through node links,and the flow relationships between key nodes were coupled to develop a coupling simulation method based on reservoir-wellbore-pipe network integration,enabling the whole process simulation of oil and gas field develop- ment. In the West A Block of Chengdao Oilfield,the production system,proration,injection allocation,and other aspects were optimized according to the actual production demand and constraints of offshore oil fields. In addition,the whole pro- cess optimization based on the integration model maximized the production benefit of the oil and gas field development sys- tem.

  • 油气田开发是一个高度连续的生产系统,涵盖了流体从油藏到井底的渗流、从井底到井口的井筒内垂直管流和井口到集输管线的地面集输管流等生产环节,各环节相互依存、相互制约。目前主要通过数值模拟方法模拟地下油藏的供液能力,缺少对井筒举升能力和地面管网集输能力的研究,并且井筒-管网系统中流体流动型态复杂,难以精确模拟流体流动规律。油藏-井筒-管网沿程相互制约的节点多,各节点遵守质量守恒和能量守恒,难以进行一体化模型耦合。油气田开发过程中,流体从油藏到地面的流动需要在协调统一的环境下才能高效运行,如何在局部最优基础上进行全流程优化,获得效率最优的系统开发方案。为此,开展油藏-井筒-管网的一体化耦合模拟和优化方法研究,集成地下油藏、井筒及地面管网等生产系统,实现整个油气田开发系统高效运行及立体优化。一体化的油藏管理作为智能油田建设的核心,随着智能油田的建设需求日益迫切[1-2],因此研究油藏-井筒-管网一体化耦合模拟方法对智能油田建设具有重要意义。

  • 1 模型建立及模拟方法

  • 一体化模型包括3个基本模型,即油藏模型、井筒模型和管网模型。油藏模型建立以及模型拟合方法相对比较成熟,本次研究重点为井筒和地面管网的模型建立及模拟方法。

  • 1.1 井筒模型建立及举升动态模拟

  • 通过建立井身结构模型、热传导模型,优选多相管流相关式,构建反映井筒举升能力的井筒模型,完成举升动态模拟。

  • 1.1.1 井筒模型建立

  • 井筒模型主要包括基于井筒管柱结构、举升方式和举升参数的井身结构模型及考虑流体能量传递的热传导模型。

  • 井身结构模型  实际井身结构复杂,通过优选井筒内关键节点,设置参数进行模拟(表1),建立精细井身结构模型(图1)。

  • 表1 不同类型井的井身结构关键节点参数

  • Table1 Key node parameters of casing programs in different types of wells

  • 图1 电潜泵井井身结构模型示意

  • Fig.1 Casing program model of well equipped with electric submersible pump

  • 热传导模型  假定传递的能量是由流体温度下降时产生的热量获得,考虑流体能量传递的热传导模型依据文献[3-5]建立。

  • 1.1.2 举升动态模拟

  • 储层及井类型、流体类型、举升方式、流体流动方式、完井方式和防砂方式主要影响井筒内流体多相管流流动型态、流体流入动态。流体流入动态产能计算公式优选及拟合影响油藏模型供液能力,油藏模型建立及模拟依据文献[6-15]。受压力、温度和流速影响,井筒中流体流动是多相管流,具有泡流、雾流、段塞流等多种复杂流动型态,其影响着井筒的举升能力。根据实际井筒内流体流动特点,对不同多相管流相关式适应性进行分析[16-19] (表2),结合井类型、类别和流体特征,优选多相管流相关式,建立反映生产实际的多相管流模型,精确描述井筒内流体流动型态(图2)。在此基础上,以井底为节点,将油藏流入动态作为输入边界条件,分析日产液量与垂直管流压力的关系,生成举升动态曲线(图3),为一体化模型提供井筒举升动态模拟的基础数据。

  • 表2 多相管流相关式适应性分析

  • Table2 Adaptive analysis of multiphase pipe flow correlations

  • 图2 井筒内流体流动型态

  • Fig.2 Flow pattern in wellbore

  • 图3 井筒举升动态曲线

  • Fig.3 Dynamic curve during wellbore lifting

  • 1.2 管网模型建立及流动保障模拟

  • 1.2.1 管网模型建立

  • 管网模型包括地面集输管网以及各种地面设备(如分离器、加热器、压缩机、外输泵等)、管线内设备(如各种阀门、管线内分离装置等)。通过定义井、关键设备、管线实际地理位置,确定管线真实拓扑结构、管线长度、高程差,建立管网模型。基于优选水平多相管流相关式如Beggs and Brill相关式,给定井口日产液量、含水率、气油比、外输终端压力和温度,计算地面管线沿程节点流量、压力和温度分布,实现管网节点流动模拟。

  • 1.2.2 流动保障模拟

  • 基于管道内流体流动状态,对比临界流量、压力、温度等指标,预测可能发生的出砂、结蜡、水合物、管道积液等异常问题,实现管网内实时预警、诊断分析和优化。

  • 砂粒临界运移速度  油井出砂是砂岩油层开采过程中的常见问题之一。当井筒内流体速度大于砂粒临界运移速度时,砂粒被提升到井口,造成油井出砂[16]。砂粒临界运移速度与砂粒粒径、形状、密度和流体密度有关。

  • 砂粒临界运移速度模型[17] 为:

  • V=-3N+9N2+gr2ρfρs-ρf(0.014476+0.19841r)0.5ρf(0.011607+0.19881r)
    (1)
  • 临界侵蚀速度[19-20] 可以表示为:

  • Ve=Cρm
    (2)
  • 临界携液速度  气井井筒积液主要是由于地层压力下降,导致气体流速降低,没有足够的能量把井筒中产生的液滴携带出井口时,液滴将在井底形成积液,严重影响正常生产[16]。根据 Tuner 模型[20],气体能拖动液滴的最小速度为:

  • vt=20.4×σ1/4ρ1-ρg1/4ρg1/2
    (3)
  • 1.3 一体化模型建立及耦合模拟

  • 按照节点的流入流出关系将油藏模型、井筒模型、管网模型等链接,建立一体化模型,耦合关键节点(如井底、井口和分离器等)流动关系,实现油藏、井筒、管网等模型上、下游各节点模拟。

  • 根据实际生产需求,建立不同类型的模型,耦合模型之间传递的变量也各不相同(表3)。

  • 油藏-井筒-管网一体化模型中各子系统之间耦合传递是动态耦合。运用油藏模型计算地层压力和饱和度,传递给井筒模型,流入动态曲线(IPR) 与井筒模型的垂直管流曲线交点即井底流压和日产液量,根据流体流速、管径、粗糙度计算管网模型的沿程压力降,压力损失误差限制在管网,直到达到管网模型耦合迭代收敛,进而实现一体化系统平衡。

  • 表3 不同耦合模型间传递的变量

  • Table3 Variables passed between different coupling models

  • 油藏模型为:

  • 2pr2+1r×pr=ϕμCtK×pt
    (4)
  • 初始条件为:

  • t=0,p=pi
    (5)
  • 外边界条件为:

  • r=re,pr=0
    (6)
  • 内边界条件为:

  • r=rw,rpr=qμ2πKh
    (7)
  • 井筒模型为:

  • -dpdz=ρ1HL+ρg1-HLgsinθ+λGv2DA1-ρ1HL+ρg1-HLvvsgp
    (8)
  • 管网模型为:

  • Δpγ=λld×v22g
    (9)
  • 其中:

  • λ=0.31644Re4
    (10)
  • 以井底为耦合节点为例,对于计算产量与实际产量存在误差的井,因井口日产液量、压力等数据和管网数据均为实测值,在一体化模型拟合时不需调整管网模型,通过调整油藏模型参数或井筒模型参数来拟合井底流压和日产液量。可调整的油藏模型参数有地层压力和采液指数;井筒模型的调整,可检查举升曲线是否存在外推情况,检查泵效或校正多相管流相关式等。

  • 2 一体化全流程优化

  • 油藏-井筒-管网一体化全流程优化是基于油藏-井筒-管网一体化模型,从油气田开发系统多节点、全过程、多目标进行全系统优化,最大程度提高系统生产效率,节约运行成本。

  • 一体化全流程优化包括优化目标确定、约束条件设置和控制变量优选 3 个关键环节,整体思路是根据实际生产设备的能力和生产需求,设定节点约束条件、系统约束条件,明确敏感性控制变量,利用最优算法进行系统求解或优化,在此过程中,反复调整敏感性控制变量进行迭代,直到满足整个系统约束条件,达到系统最优目标。

  • 优化目标确定  全流程包括单井、注采源汇端、全局等多层次目标,可以优化产量、效益、热能和集输管线路径等多个目标。

  • 约束条件设置  约束条件包括系统约束、井约束及举升约束等,全流程优化可以同时满足多种约束条件,对于不同的系统和举升方式,依据生产实际选择相应的约束条件(表4)。

  • 表4 约束条件设置

  • Table4 Constraint setting

  • 控制变量优选  控制变量是井筒内可控设备的控制参数,控制生产井日产液量或注入井日注水量、日注气量。控制变量的选取要满足约束条件对控制变量的敏感性。不同举升方式的控制变量如表5所示。对于自喷井,控制设备是油嘴,控制变量是井口节流压差;对于抽油机井,控制设备是泵,控制变量是冲次;对于气举井,控制设备是注气设备或气举阀,控制变量是注气速率或井口节流压差。

  • 表5 控制变量优选

  • Table5 Selection of control variables

  • 油藏-井筒-管网一体化全流程优化结果包含油藏、单井、节点、管道、分离器等全流程各节点日产液量、压力、含水率、气油比等参数,具有全面性和实时性特点。

  • 3 实例应用

  • 基于埕岛油田西 A 区块开展油藏-井筒-管网一体化耦合模拟方法的应用,研究油藏-井筒-管网一体化生产潜力,查找单井生产瓶颈,优化生产制度及科学合理的配产配注方案,实现油藏-井筒-管网一体化全流程最优产油量。

  • 3.1 单井油藏-井筒-管网一体化潜力分析

  • 以埕岛油田西 A 区块电潜泵井 35I 井为例,某时刻35I井泵频率为35 Hz,日产液量为46 m3 /d。首先建立该井油藏-井筒-管网一体化模型,在模型拟合时,给定井口日产液量、含水率、气油比,调整影响采液指数的参数如渗透率、有效厚度、表皮系数等,拟合井底流压和日产液量(图4);通过调整电泵磨损系数和气体分离效率来拟合泵吸入口、排出口压力(图5),调整总导热系数拟合温度剖面(图5),日产液量、压力、温度拟合好后的油藏-井筒-管网一体化模型反映了实际生产条件,可以用来进行一体化潜力分析。

  • 图4 埕岛油田西A区块35I井协调点拟合

  • Fig.4 Coordination point fitting of Well35I in Chengdao West A Block

  • 图5 埕岛油田西A区块35I井压力、温度剖面拟合

  • Fig.5 Pressure/temperature profile fitting of Well35I in Chengdao West A Block

  • 井筒方面的单井潜力主要体现在垂直管流曲线上,矿场实际中电潜泵的泵频率可调范围最大为 60 Hz,利用一体化模型对泵频率进行敏感性分析,井筒方面单井生产潜力为46~95 m3 /d。油藏方面潜力主要体现在流入动态曲线上,日产液量范围受合理生产压差、最小井底流压、最大采液量、采油速度等影响因素限制。采用合理生产压差条件下油藏模型计算结果,油藏方面单井生产潜力为 50 m3 /d。

  • 一体化模型拟合好后,以单井日产油量最优为目标、以泵频率为控制变量进行优化。单井一体化潜力为油藏潜力和井筒-管网潜力的最小值,单井一体化生产潜力为 50 m3 /d,对应泵频率为 38 Hz。因此,制约单井生产的瓶颈为井筒,提高泵频率至 38 Hz,可小幅度提液;当单井日产液量达到 50 m3 /d 后,油藏产能成为制约单井生产的瓶颈,下步可通过采取增产措施、优化注采等提高油藏供液能力,实现单井提质增效。

  • 3.2 注采系统一体化全流程优化

  • 海上油藏在开发初期,平台处理能力有限,应用全流程优化方法,满足处理能力限制的前提下协调整个注采系统,优化生产、注入系统配产配注,达到整个平台产油量最优目标。

  • 以生产系统作为主模型,以注入系统作为关联模型,建立埕岛油田西 A 区块注采井组的一体化模型,生产系统有 2个油藏即油藏 A 和 B,其产量通过井口管汇节点 A 汇入节点 B,最终经分离器节点汇入分离器中(图6a);注入系统的注水量经节点C、节点B和节点A汇入A油藏,通过节点A实现了油藏A 注采系统之间的关联(图6b)。

  • 设定一体化模型的优化目标是满足平台处理能力下日产油量最优;设置生产系统地层压力和注入系统注入压力为约束条件;设置生产井、注水井井口节流压差为控制变量,以调整生产井日产液量和注水井日注水量。

  • 优化计算中,利用序列二次规划最优算法,根据生产系统产量和时间步长,通过产量计算亏空量,传递给注水系统,以累积注水量作为注水系统约束条件,结合生产、注水系统约束条件,调整井口节流压差,寻求平台最优产油量。

  • 对比不考虑平台能力的单系统配产配注方案和一体化优化配产配注方案的日产油量(图7a)、日产液量(图7b)和地层压力(图7c)发现,配产配注方案目标同样是日产油量,单系统配产配注方案日产油量比一体化优化配产配注方案高,但其日产液量远超过平台最大处理能力,地层压力也不合理;而一体化优化配产配注方案同时考虑了地层压力、井筒条件、地面注入压力、地面处理能力,更科学、全面、贴近生产实际,实现了井组的高效开发。

  • 图6 埕岛油田西A区块注采井组的一体化模型

  • Fig.6 Integrated model of injection-production well group in Chengdao West A Block

  • 图7 埕岛油田西A区块单系统和一体化优化配产配注方案对比

  • Fig.7 Comparison between single-system production and injection program and integrated program with optimized proration and injection allocation in Chengdao West A Block

  • 4 结论

  • 在井筒模型、管网模型建立及模拟基础上,通过对多系统、多节点耦合模拟,形成了油藏-井筒管网一体化耦合模拟方法。油藏-井筒-管网一体化耦合模拟消除了单系统模拟时边界条件影响,用更精确、更高频率的井筒及管网实时数据来校正一体化模型,更准确地模拟整个开发生产系统,实现了更贴近生产实际的油气田开发全流程模拟。

  • 基于油藏-井筒-管网一体化耦合模拟,考虑油气田开发系统实际生产目标、约束条件和控制变量,从整个系统多节点、全过程、多目标进行油藏-井筒-管网全流程优化。

  • 基于埕岛油田西 A 区块,应用一体化全流程优化查找单井生产瓶颈,优化生产制度、优化井组配产配注,实现了油藏-井筒-管网全流程最优,提高油气田开发系统智能化管理水平。

  • 符号解释

  • A ——管流横截面积,m2

  • C ——常数,一般取值为100~250;

  • C t——综合压缩系数,MPa-1

  • D ——管流直径,m;

  • g ——重力加速度,cm/s2,取值为980;

  • Gν——混合物质量流量,kg/s;

  • h ——油藏厚度,m;

  • H L——持液率;

  • K ——油藏渗透率,mD;

  • l ——管线长度,m;

  • N ——流体动态黏度,mPa·s;

  • p ——地层压力,MPa;

  • p i ——原始地层压力,MPa;

  • p ——沿程压差,MPa;

  • q ——日产液量,m3 /d;

  • r ——颗粒半径,cm;

  • r e——泄油半径,m;

  • r w——井筒半径,m;

  • Re ——雷诺数;

  • t ——时刻,s;

  • μ——原油黏度,mPa·s;

  • v ——流体流速,m/s;

  • v sg——气体流速,m/s;

  • V ——临界运移速度,cm/s;

  • V e——临界侵蚀速度,cm/s;

  • v t ——气体能拖动液滴的最小速度,m/s;

  • z ——井筒内深度,m;

  • θ——管流倾斜角度,(°);

  • ρf ——流体密度,g/cm3

  • ρs ——砂粒密度,g/cm3

  • ρg——气体密度,g/cm3

  • ρl ——液相密度,g/cm3

  • ρm——混合液密度,g/cm3

  • γ——重度,N/m3

  • σ——气液表面张力,N/m;

  • ϕ——孔隙度;

  • λ——沿程阻力系数。

  • 参考文献

    • [1] 梁文福.油田开发智能应用系统建设成果及展望[J].大庆石油地质与开发,2019,38(5):283-289.LIANG Wenfu.Constructed achievements and prospects of the in⁃ telligent application system for the oilfield development[J].Petro⁃ leum Geology & Oilfield Development in Daqing,2019,38(5):283-289.

    • [2] 杜金虎,时付更,杨剑锋,等.中国石油上游业务信息化建设总体蓝图[J].中国石油勘探,2020,25(5):1-8.DU Jinhu,SHI Fugeng,YANG Jianfeng,et al.Overall blueprint of information construction of PetroChina upstream business[J].Chi⁃ na Petroleum Exploration,2020,25(5):1-8.

    • [3] SUTTON R P,FARSHAD F.Evaluation of empirically derived PVT properties for Gulf of Mexico crude oils[C].SPE 13172,1990.

    • [4] BEAL C.The viscosity of air,water,natural gas,crude oil and its associated gases at oil field temperatures and pressures[C].SPE 946094,1946.

    • [5] CARLAW H S,JEGER J C.Conduction of heat in solids[M].Ox⁃ ford:Clarendon Press,1959:29-111.

    • [6] 秦积舜,李爱芬.油层物理学[M].东营:石油大学出版社,2003:76-80.QIN Jishun,LI Aifen.Petrophysics[M].Dongying:Press of Univer⁃ sity of Petroleum,China,2003:76-80.

    • [7] 姜瑞忠,刘秀伟,崔永正,等.非稳态窜流多段压裂水平井井底压力分析[J].油气地质与采收率,2019,26(5):86-95.JIANG Ruizhong,LIU Xiuwei,CUI Yongzheng,et al.Bottomhole pressure analysis of multistage fractured horizontal well during unsteady crossflow[J].Petroleum Geology and Recovery Efficien⁃ cy,2019,26(5):86-95.

    • [8] BEGGS H D,ROBINSON J R.Estimating the viscosity of crude oil systems[C].SPE 5434,1975.

    • [9] 纪禄军,赵金洲,杨承孝,等.IPR方程通式中Vogel参数与流动效率R的计算[J].河南石油,2000,(4):13-15,60.JI Lujun,ZHAO Jinzhou,YANG Chengxiao.Calculation of Vogel parameter and flow efficiency R in IPR general formula[J].Henan Petroleum,2000,(4):13-15,60.

    • [10] VOGEL J V.Inflow performance relationships for solution gas drive wells[C].SPE 1476,1968.

    • [11] 张琪.采油工程原理与设计[M].东营:石油大学出版社,2000:1-15.ZHANG Qi.Principle and design of oil production engineering [M].Dongying:Press of University of Petroleum,China,2000:1-15.

    • [12] FETKOVICH M J.The isochronal testing of oil wells[C].SPE 4529,1973.

    • [13] JONES L G,BLOUNT E M.Use of short term multiple rate flow tests to predict performance of wells having turbulence[C].SPE 6133,1976.

    • [14] CINCO LEY H,SAMANIEGO F V,DOMINGUEZ N A.Transient pressure behaviour for a well with a finite-conductivity vertical fracture[C].SPE 6014,1976.

    • [15] CINCO L H,RAMEY H J,MILLER F G.Pseudo-skin factors for partially-penetrating directionally-drilled wells[C].SPE 5589,1975.

    • [16] 吴星晔.储层-井筒-井口一体化动态分析方法研究[D].成都:西南石油大学,2017.WU Xingye.The dynamic analysis method of reservoir-wellboewellhead integration[D].Chengdu:Southwest Petroleum Universi⁃ ty,2017.

    • [17] WASP E J,KENNY J P,GANDHI R L.Solid-liquid flow slurry pipe line transportation[M].Houston:Gulf Publishing Company,1979:25-122.

    • [18] SVEDEMAN S J.Criteria for sizing multiphase flowlines for ero⁃ sive/corrosive service[C].SPE 26569,1994.

    • [19] 李学军.气井冲蚀腐蚀临界速度的计算方法[J].石油石化节能,1991,(4):94.LI Xuejun.Calculation method of critical velocity of erosion and corrosion in gas well[J].Energy Conservation in Petroleum and Petrochemical Industry,1991,(4):94.

    • [20] TURNER R G.Analysis and prediction of minimum flow rate for the continuous removal of liquids from gas well[C].SPE 2198,1969.

  • 参考文献

    • [1] 梁文福.油田开发智能应用系统建设成果及展望[J].大庆石油地质与开发,2019,38(5):283-289.LIANG Wenfu.Constructed achievements and prospects of the in⁃ telligent application system for the oilfield development[J].Petro⁃ leum Geology & Oilfield Development in Daqing,2019,38(5):283-289.

    • [2] 杜金虎,时付更,杨剑锋,等.中国石油上游业务信息化建设总体蓝图[J].中国石油勘探,2020,25(5):1-8.DU Jinhu,SHI Fugeng,YANG Jianfeng,et al.Overall blueprint of information construction of PetroChina upstream business[J].Chi⁃ na Petroleum Exploration,2020,25(5):1-8.

    • [3] SUTTON R P,FARSHAD F.Evaluation of empirically derived PVT properties for Gulf of Mexico crude oils[C].SPE 13172,1990.

    • [4] BEAL C.The viscosity of air,water,natural gas,crude oil and its associated gases at oil field temperatures and pressures[C].SPE 946094,1946.

    • [5] CARLAW H S,JEGER J C.Conduction of heat in solids[M].Ox⁃ ford:Clarendon Press,1959:29-111.

    • [6] 秦积舜,李爱芬.油层物理学[M].东营:石油大学出版社,2003:76-80.QIN Jishun,LI Aifen.Petrophysics[M].Dongying:Press of Univer⁃ sity of Petroleum,China,2003:76-80.

    • [7] 姜瑞忠,刘秀伟,崔永正,等.非稳态窜流多段压裂水平井井底压力分析[J].油气地质与采收率,2019,26(5):86-95.JIANG Ruizhong,LIU Xiuwei,CUI Yongzheng,et al.Bottomhole pressure analysis of multistage fractured horizontal well during unsteady crossflow[J].Petroleum Geology and Recovery Efficien⁃ cy,2019,26(5):86-95.

    • [8] BEGGS H D,ROBINSON J R.Estimating the viscosity of crude oil systems[C].SPE 5434,1975.

    • [9] 纪禄军,赵金洲,杨承孝,等.IPR方程通式中Vogel参数与流动效率R的计算[J].河南石油,2000,(4):13-15,60.JI Lujun,ZHAO Jinzhou,YANG Chengxiao.Calculation of Vogel parameter and flow efficiency R in IPR general formula[J].Henan Petroleum,2000,(4):13-15,60.

    • [10] VOGEL J V.Inflow performance relationships for solution gas drive wells[C].SPE 1476,1968.

    • [11] 张琪.采油工程原理与设计[M].东营:石油大学出版社,2000:1-15.ZHANG Qi.Principle and design of oil production engineering [M].Dongying:Press of University of Petroleum,China,2000:1-15.

    • [12] FETKOVICH M J.The isochronal testing of oil wells[C].SPE 4529,1973.

    • [13] JONES L G,BLOUNT E M.Use of short term multiple rate flow tests to predict performance of wells having turbulence[C].SPE 6133,1976.

    • [14] CINCO LEY H,SAMANIEGO F V,DOMINGUEZ N A.Transient pressure behaviour for a well with a finite-conductivity vertical fracture[C].SPE 6014,1976.

    • [15] CINCO L H,RAMEY H J,MILLER F G.Pseudo-skin factors for partially-penetrating directionally-drilled wells[C].SPE 5589,1975.

    • [16] 吴星晔.储层-井筒-井口一体化动态分析方法研究[D].成都:西南石油大学,2017.WU Xingye.The dynamic analysis method of reservoir-wellboewellhead integration[D].Chengdu:Southwest Petroleum Universi⁃ ty,2017.

    • [17] WASP E J,KENNY J P,GANDHI R L.Solid-liquid flow slurry pipe line transportation[M].Houston:Gulf Publishing Company,1979:25-122.

    • [18] SVEDEMAN S J.Criteria for sizing multiphase flowlines for ero⁃ sive/corrosive service[C].SPE 26569,1994.

    • [19] 李学军.气井冲蚀腐蚀临界速度的计算方法[J].石油石化节能,1991,(4):94.LI Xuejun.Calculation method of critical velocity of erosion and corrosion in gas well[J].Energy Conservation in Petroleum and Petrochemical Industry,1991,(4):94.

    • [20] TURNER R G.Analysis and prediction of minimum flow rate for the continuous removal of liquids from gas well[C].SPE 2198,1969.