CN116341423B - 油水两相流滑移速度模型的计算方法 - Google Patents

油水两相流滑移速度模型的计算方法 Download PDF

Info

Publication number
CN116341423B
CN116341423B CN202310618582.9A CN202310618582A CN116341423B CN 116341423 B CN116341423 B CN 116341423B CN 202310618582 A CN202310618582 A CN 202310618582A CN 116341423 B CN116341423 B CN 116341423B
Authority
CN
China
Prior art keywords
oil
water
speed
flow
phase
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
CN202310618582.9A
Other languages
English (en)
Other versions
CN116341423A (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.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum 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 Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN202310618582.9A priority Critical patent/CN116341423B/zh
Publication of CN116341423A publication Critical patent/CN116341423A/zh
Application granted granted Critical
Publication of CN116341423B publication Critical patent/CN116341423B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Lubricants (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明公开了油水两相流滑移速度模型的计算方法，包括以下步骤：S1.建立油水滑移速度的基本公式；S2.结合粘度的影响，修改油水滑移速度的基本公式；S3.根据油水流型分别建立油水滑移速度模型，并确立模型参数。本发明公开的油水两相流滑移速度模型的计算方法，通过建立油水两相流滑移速度模型有效地解决了现有技术中由于没有考虑油水粘度和流态对油水滑移速度的影响，容易产生与真实值相差甚远的问题。

Description

油水两相流滑移速度模型的计算方法
技术领域
本发明涉及井筒多相流领域，具体为油水两相流滑移速度模型的计算方法。
背景技术
近年来，页岩油气作为一种非常规油气资源，因其储量丰富和能源清洁的特点而成为全球油气勘探开发的热点。在油井生产过程中井筒流体通常是以油气水三相的形式流动，但在许多时候为了简便通常会研究油水两相流动。油水滑移现象在井筒多相流中十分常见，可以通过油水滑移速度计算持水率的大小，持水率结合持气率则可以确定井筒内混合密度、混合粘度、压降等参数，同时在测井工作中持水率有着对井筒状况分析的作用，对油气田的开发有着重要的价值。
目前国内外学者对油水滑移速度的研究多是基于滑移机理并结合实验现象，建立滑移速度模型，这些模型中的系数是通过实验数据拟合确定的。但在实际的研究中可以发现，这些模型存在一定的局限性，具体表现为这类模型在特定的实验数据下能够取得较好的准确性，但在一般的实验条件下则表现出误差超过60%的特点，无法准确的描述滑移速度的大小。在油水滑移速度模型的研究中，许多模型只考虑了浮力、表面张力等参数的影响，没有考虑油水粘度和流态对油水滑移速度的影响，容易产生与真实值相差甚远的情况。
发明内容
本发明提供油水两相流滑移速度模型的计算方法，包括以下步骤：S1.建立油水滑移速度的基本公式；S2.结合粘度的影响，修改油水滑移速度的基本公式；S3.根据油水流型分别建立油水滑移速度模型，并确立模型参数。
进一步地，所述S1步骤中建立的油水滑移速度的基本公式为：，式中，Vsow为油相真实速度，m/s；Vm为油水混合速度，m/s；V∞为液滴终端上升速度，m/s；H0为持油率，无因次；C0为系数，无因次；k为指数，无因次。
进一步地，所述S2步骤中修改后的油水滑移速度的基本公式为：，其中，μγ为相对粘度，mPa·s；Vsow为油相真实速度，m/s；Vm为油水混合速度，m/s；V∞为液滴终端上升速度，m/s；H0为持油率，无因次；C0为系数，无因次；k为指数，无因次；f（μγ）为相对粘度μγ的计算公式。
进一步地，所述μγ为相对粘度的计算公式为：，式中，μ0为油相粘度，mPa·s；μw为水相粘度，mPa·s。
进一步地，所述S1步骤中液滴终端上升速度V∞依据油水比不同以及计算位置不同，可划分为：油滴流、水滴流以及油水段塞流。
进一步地，所述油滴流的计算公式为：，式中，ρw为水相密度，kg/m3；ρ0为油相密度，kg/m3；g为重力加速度，m/s2；σ为油水表面张力，N/m；b为系数，无因次。
进一步地，所述水滴流的计算公式为：，式中，ρw为水相密度，kg/m3；ρ0为油相密度，kg/m3；g为重力加速度，m/s2；σ为油水表面张力，N/m；b为系数，无因次。
进一步地，所述油水段塞流的计算公式为：，式中，D为管径，m；ρw为水相密度，kg/m3；ρ0为油相密度，kg/m3；g为重力加速度，m/s2；b为系数，无因次。
进一步地，所述S3步骤中通过分别对油滴流、水滴流以及油水段塞流进行拟合，确定修改后的滑移速度公式中的参数。
进一步地，所述油水段塞流的油水滑移速度的基本公式为：，/>，式中，D为管径，m；μγ为相对粘度，mPa·s；Vsow为油相真实速度，m/s；ρw为水相密度，kg/m3；ρ0为油相密度，kg/m3；g为重力加速度，m/s2；Vm为油水混合速度，m/s；V∞为液滴终端上升速度，m/s；H0为持油率，无因次。
进一步地，所述水滴流的油水滑移速度的基本公式为：，/>，式中，μγ为相对粘度，mPa·s；Vsow为油相真实速度，m/s；ρw为水相密度，kg/m3；ρ0为油相密度，kg/m3；g为重力加速度，m/s2；Vm为油水混合速度，m/s；V∞为液滴终端上升速度，m/s；H0为持油率，无因次；σ为油水表面张力，N/m。
进一步地，所述油滴流的油水滑移速度的基本公式为：，/>，式中，μγ为相对粘度，mPa·s；Vsow为油相真实速度，m/s；ρw为水相密度，kg/m3；ρ0为油相密度，kg/m3；g为重力加速度，m/s2；Vm为油水混合速度，m/s；V∞为液滴终端上升速度，m/s；H0为持油率，无因次；σ为油水表面张力，N/m。
进一步地，计算结果的收敛值范围为10-10以内。
本发明提供油水两相流滑移速度模型的计算方法，有效地解决了现有技术中由于没有考虑油水粘度和流态对油水滑移速度的影响，容易产生与真实值相差甚远的问题。
附图说明
图1为本发明提供油水两相流滑移速度模型的计算方法的模型建立流程图；
图2为本发明提供油水两相流滑移速度模型的计算方法的Davarzani 和 Miller(1983)实验数据与预测数据比较示意图；
图3为本发明提供油水两相流滑移速度模型的计算方法的Flores(1997)实验数据与预测数据比较示意图；
图4为本发明提供油水两相流滑移速度模型的计算方法的Lucas G P（2001）实验数据与预测数据比较示意图；
图5为本发明提供油水两相流滑移速度模型的计算方法的王振亚(2010)实验数据与预测数据比较示意图。
具体实施方式
以下结合附图对本发明的实施方法进行详细说明，所描述的仅为部分实施例，并非全部实施例，为了清楚的目的，在附图及说明中省略了与本发明无关的表示及描述。
如图1所示，本发明提供油水两相流滑移速度模型的计算方法，包括以下步骤：包括以下步骤：S1.建立油水滑移速度的基本公式；S2.结合粘度的影响，修改油水滑移速度的基本公式；S3.根据油水流型分别建立油水滑移速度模型，并确立模型参数。
其中，S1步骤中建立的油水滑移速度的基本公式为：，式中，Vsow为油相真实速度，m/s；Vm为油水混合速度，m/s；V∞为液滴终端上升速度，m/s；H0为持油率，无因次；C0为系数，无因次；k为指数，无因次。
S2步骤中修改后的油水滑移速度的基本公式为：，其中，μγ为相对粘度，mPa·s。相对粘度μγ的计算公式为：/>，式中，μ0为油相粘度，mPa·s；μw为水相粘度，mPa·s。
其中，液滴终端上升速度V∞依据油水比不同以及计算位置不同，可划分为：油滴流、水滴流以及油水段塞流。油滴流的液滴终端上升速度计算公式为：；水滴流的液滴终端上升速度计算公式为：/>；油水段塞流的液滴终端上升速度计算公式为：/>，式中，D为管径，m；ρw为水相密度，kg/m3；ρ0为油相密度，kg/m3；g为重力加速度，m/s2；σ为油水表面张力，N/m；b为系数，无因次。
为确定滑移速度模型中参数f(μr)，C0，b，k的大小，采用包维尔优化算法，代入油相表观速度、水相表观速度、持水率和滑移速度等实验数据对滑移速度模型进行拟合，进行1000次迭代，并将收敛值控制在10-10范围内，分别对水滴流、油水段塞流和油滴流三种滑移速度公式进行拟合。
建立油滴流的滑移速度方程：
，
；
建立水滴流的滑移速度方程：
，
；
建立油水段塞流的滑移速度方程：
，
，式中，μγ为相对粘度，mPa·s；Vsow为油相真实速度，m/s；ρw为水相密度，kg/m3；ρ0为油相密度，kg/m3；g为重力加速度，m/s2；Vm为油水混合速度，m/s；V∞为液滴终端上升速度，m/s；H0为持油率，无因次；σ为油水表面张力，N/m；D为管径，m。
实施例：
选择Davarzani和Miller(1983)、Flores(1997)、Lucas G P（2001）和王振亚(2010)四组油水两相流实验数据，根据油滴流、水滴流和油水段塞流三种流态对实验数据分成三类。确定实验条件中的表面张力、油粘度、水粘度、油密度、水密度和管径等实验条件，并选择数据中的油相表观速度、水相表观速度、持水率和混合速度代入到所建立的三类油水滑移速度方程中计算，得到公式预测的滑移速度大小。将公式预测的滑移速度与实验测得的滑移速度进行比较，分析滑移速度模型的准确性。
实验数据如图2、图3、图4、图5所示，随着持水率的增大，变化趋势基本保持一致，同时误差分别为5.2%，6.3%，4.8%和5.4%，预测结果与实验结果误差较小，模型准确度较高。
以上所述仅是本发明的优选实施方式，应当理解本发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。

Claims (6)

1.油水两相流滑移速度模型的计算方法，其特征在于，包括以下步骤：
S1.建立油水滑移速度的基本公式；
S2.结合粘度的影响，修改滑移速度公式；
S3.根据油水流型分别建立油水滑移速度模型，并确立模型参数；
所述S1步骤中建立的油水滑移速度的基本公式为：，式中，Vsow为油相真实速度，m/s；Vm为油水混合速度，m/s；V∞为液滴终端上升速度，m/s；H0为持油率，无因次；C0为系数，无因次；k为指数，无因次；
所述S2步骤中修改后的滑移速度公式为：，其中，μγ为相对粘度，mPa·s；
所述μγ的计算公式为：f（μγ）=，式中，μ0为油相粘度，mPa·s；μw为水相粘度，mPa·s；
所述S1步骤中液滴终端上升速度V∞依据油水比不同以及计算位置不同，划分为：油滴流、水滴流以及油水段塞流；
所述油滴流的计算公式为：，式中，ρw为水相密度，kg/m3；ρ0为油相密度，kg/m3；g为重力加速度，m/s2；σ为油水表面张力，N/m；b为系数，无因次；
所述水滴流的计算公式为：；
所述油水段塞流的计算公式为：，式中，D为管径，m。
2.根据权利要求1所述的油水两相流滑移速度模型的计算方法，其特征在于，所述S3步骤中通过分别对油滴流、水滴流以及油水段塞流进行拟合，确定修改后的滑移速度公式中的参数。
3.根据权利要求2所述的油水两相流滑移速度模型的计算方法，其特征在于，所述油水段塞流的滑移速度的方程为：
，/>。
4.根据权利要求2所述的油水两相流滑移速度模型的计算方法，其特征在于，所述水滴流的滑移速度的方程为：
，/>。
5.根据权利要求2所述的油水两相流滑移速度模型的计算方法，其特征在于，所述油滴流的滑移速度的方程为：
，/>。
6.根据权利要求3-5任意一条所述的油水两相流滑移速度模型的计算方法，其特征在于，计算结果的收敛值范围为10-10以内。
CN202310618582.9A 2023-05-30 2023-05-30 油水两相流滑移速度模型的计算方法 Active CN116341423B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310618582.9A CN116341423B (zh) 2023-05-30 2023-05-30 油水两相流滑移速度模型的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310618582.9A CN116341423B (zh) 2023-05-30 2023-05-30 油水两相流滑移速度模型的计算方法

Publications (2)

Publication Number Publication Date
CN116341423A CN116341423A (zh) 2023-06-27
CN116341423B true CN116341423B (zh) 2023-08-01

Family

ID=86880765

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310618582.9A Active CN116341423B (zh) 2023-05-30 2023-05-30 油水两相流滑移速度模型的计算方法

Country Status (1)

Country Link
CN (1) CN116341423B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117454063B (zh) * 2023-12-26 2024-03-12 西南石油大学 一种井筒油-气-水多相流流态判别及持水率计算方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114991747A (zh) * 2022-05-26 2022-09-02 长江大学 一种基于环形阵列探针测量的页岩油产量解释方法
CN115906695A (zh) * 2022-11-18 2023-04-04 西南石油大学 一种多信息融合的产出剖面测井解释优化方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8196395B2 (en) * 2009-06-29 2012-06-12 Lightsail Energy, Inc. Compressed air energy storage system utilizing two-phase flow to facilitate heat exchange
US10041019B2 (en) * 2014-10-10 2018-08-07 Continental Automotive Systems, Inc. Drilling fluid system
CN108509703B (zh) * 2018-03-22 2022-01-28 中国石油大学(华东) 一种气藏状态参数随钻数值反演分析方法
CN111948341B (zh) * 2020-07-27 2021-09-21 中国石油大学(北京) 固液两相条件下气泡上升速度的测试装置及测试方法
CN113361161B (zh) * 2021-06-02 2022-07-29 中国石油大学(华东) 基于数字岩心的页岩油水两相流动模拟方法、装置及存储介质
CN114609010B (zh) * 2022-03-02 2024-02-02 中国石油大学(华东) 一种页岩储层油水相对渗透率的测定方法及装置
CN115270488A (zh) * 2022-08-07 2022-11-01 西南石油大学 一种预测页岩油储层自吸采收率的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114991747A (zh) * 2022-05-26 2022-09-02 长江大学 一种基于环形阵列探针测量的页岩油产量解释方法
CN115906695A (zh) * 2022-11-18 2023-04-04 西南石油大学 一种多信息融合的产出剖面测井解释优化方法

Also Published As

Publication number Publication date
CN116341423A (zh) 2023-06-27

Similar Documents

Publication Publication Date Title
CN116341423B (zh) 油水两相流滑移速度模型的计算方法
Hasan et al. A study of multiphase flow behavior in vertical wells
Zhou et al. A new model for predicting gas-well liquid loading
Schmidt et al. Severe slugging in offshore pipeline riser-pipe systems
CN104895560A (zh) 一种深水测试井筒压力、温度场模拟及水合物预测方法
CN106600442B (zh) 一种深水气井测试临界携液条件优化方法
Belt et al. Comparison of commercial multiphase flow simulators with experimental and field databases
CN112966418B (zh) 一种海洋天然气水合物开采立管动力特性分析方法
Wang et al. Prediction of the critical gas velocity of liquid unloading in a horizontal gas well
EP3289319A1 (en) Multiphase flow meters and related methods
Spedding et al. Three phase oil-water-gas horizontal co-current flow: I. Experimental and regime map
Hewitt Gas-liquid flow
CN115345090A (zh) 一种欠饱和煤层气储层压降漏斗动态传播的计算方法
Duan et al. Experimental and computational investigations on severe slugging in a catenary riser
CN113468826B (zh) 基于真实液膜分布的页岩气水平井临界携液流量预测方法
GB2573653A (en) Improved flow measurement
CN111648759A (zh) 一种页岩气水平井积液位置判断方法
CN115936258A (zh) 一种井筒结垢动态沉积堵塞预测模型构建方法
Zhou et al. Nonlinear adaptive observer for managed pressure drilling system
US20080262808A1 (en) Method for dimensioning industrial installations where a two-phase gas-liquid mixture flows in an intermittent regime
CN115290509A (zh) 浸润调控驱油过程中毛细管压力曲线及动态表征模型的确定方法
Cornish The vertical multiphase flow of oil and gas at high rates
Syikilili et al. Effect of pipeline inclination on multiphase flow
Capovilla et al. Experimental investigation of vertical high-velocity two-phase flows in large-diameter pipes
CN105160056A (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
GR01 Patent grant
GR01 Patent grant