CN116341423A - 油水两相流滑移速度模型的计算方法 - Google Patents
油水两相流滑移速度模型的计算方法 Download PDFInfo
- Publication number
- CN116341423A CN116341423A CN202310618582.9A CN202310618582A CN116341423A CN 116341423 A CN116341423 A CN 116341423A CN 202310618582 A CN202310618582 A CN 202310618582A CN 116341423 A CN116341423 A CN 116341423A
- Authority
- CN
- China
- Prior art keywords
- oil
- water
- phase
- speed
- dimensionless
- 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.)
- Granted
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 125
- 238000004364 calculation method Methods 0.000 title claims abstract description 31
- 230000005514 two-phase flow Effects 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 claims description 39
- 239000012071 phase Substances 0.000 claims description 32
- 230000000630 rising effect Effects 0.000 claims description 21
- 239000008346 aqueous phase Substances 0.000 claims description 18
- 230000001133 acceleration Effects 0.000 claims description 14
- 230000005484 gravity Effects 0.000 claims description 14
- 239000007788 liquid Substances 0.000 claims description 13
- 230000014759 maintenance of location Effects 0.000 claims description 13
- 238000010586 diagram Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000003079 shale oil Substances 0.000 description 1
- 238000009491 slugging Methods 0.000 description 1
Images
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(μγ)为相对粘度μγ的计算公式。
进一步地,所述S1步骤中液滴终端上升速度V∞依据油水比不同以及计算位置不同,可划分为:油滴流、水滴流以及油水段塞流。
进一步地,所述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为指数,无因次。
其中,液滴终端上升速度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 (13)
1.油水两相流滑移速度模型的计算方法,其特征在于,包括以下步骤:
S1.建立油水滑移速度的基本公式;
S2.结合粘度的影响,修改油水滑移速度的基本公式;
S3.根据油水流型分别建立油水滑移速度模型,并确立模型参数。
5.根据权利要求1所述的油水两相流滑移速度模型的计算方法,其特征在于,所述S1步骤中液滴终端上升速度V∞依据油水比不同以及计算位置不同,可划分为:油滴流、水滴流以及油水段塞流。
9.根据权利要求1所述的油水两相流滑移速度模型的计算方法,其特征在于,所述S3步骤中通过分别对油滴流、水滴流以及油水段塞流进行拟合,确定修改后的滑移速度公式中的参数。
13.根据权利要求10-12任意一条所述的油水两相流滑移速度模型的计算方法,其特征在于,计算结果的收敛值范围为10-10以内。
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 true CN116341423A (zh) | 2023-06-27 |
CN116341423B 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) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117454063A (zh) * | 2023-12-26 | 2024-01-26 | 西南石油大学 | 一种井筒油-气-水多相流流态判别及持水率计算方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102365458A (zh) * | 2009-06-29 | 2012-02-29 | 光帆能源公司 | 使用两相流以便于热交换的压缩空气能量存储系统 |
US20160102268A1 (en) * | 2014-10-10 | 2016-04-14 | Continental Automotive Systems, Inc. | Drilling fluid system |
CN108509703A (zh) * | 2018-03-22 | 2018-09-07 | 中国石油大学(华东) | 一种气藏状态参数随钻数值反演分析方法 |
CN111948341A (zh) * | 2020-07-27 | 2020-11-17 | 中国石油大学(北京) | 固液两相条件下气泡上升速度的测试装置及测试方法 |
CN113361161A (zh) * | 2021-06-02 | 2021-09-07 | 中国石油大学(华东) | 基于数字岩心的页岩油水两相流动模拟方法、装置及存储介质 |
CN114609010A (zh) * | 2022-03-02 | 2022-06-10 | 中国石油大学(华东) | 一种页岩储层油水相对渗透率的测定方法及装置 |
CN114991747A (zh) * | 2022-05-26 | 2022-09-02 | 长江大学 | 一种基于环形阵列探针测量的页岩油产量解释方法 |
CN115270488A (zh) * | 2022-08-07 | 2022-11-01 | 西南石油大学 | 一种预测页岩油储层自吸采收率的方法 |
CN115906695A (zh) * | 2022-11-18 | 2023-04-04 | 西南石油大学 | 一种多信息融合的产出剖面测井解释优化方法 |
-
2023
- 2023-05-30 CN CN202310618582.9A patent/CN116341423B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102365458A (zh) * | 2009-06-29 | 2012-02-29 | 光帆能源公司 | 使用两相流以便于热交换的压缩空气能量存储系统 |
US20160102268A1 (en) * | 2014-10-10 | 2016-04-14 | Continental Automotive Systems, Inc. | Drilling fluid system |
CN108509703A (zh) * | 2018-03-22 | 2018-09-07 | 中国石油大学(华东) | 一种气藏状态参数随钻数值反演分析方法 |
CN111948341A (zh) * | 2020-07-27 | 2020-11-17 | 中国石油大学(北京) | 固液两相条件下气泡上升速度的测试装置及测试方法 |
CN113361161A (zh) * | 2021-06-02 | 2021-09-07 | 中国石油大学(华东) | 基于数字岩心的页岩油水两相流动模拟方法、装置及存储介质 |
CN114609010A (zh) * | 2022-03-02 | 2022-06-10 | 中国石油大学(华东) | 一种页岩储层油水相对渗透率的测定方法及装置 |
CN114991747A (zh) * | 2022-05-26 | 2022-09-02 | 长江大学 | 一种基于环形阵列探针测量的页岩油产量解释方法 |
CN115270488A (zh) * | 2022-08-07 | 2022-11-01 | 西南石油大学 | 一种预测页岩油储层自吸采收率的方法 |
CN115906695A (zh) * | 2022-11-18 | 2023-04-04 | 西南石油大学 | 一种多信息融合的产出剖面测井解释优化方法 |
Non-Patent Citations (2)
Title |
---|
周勇等: "基于实验基础上的油水两相流滑脱速度求取方法", 《测井技术》, no. 6, pages 516 - 518 * |
宋红伟等: "水平井及大斜度井油水两相流动模型研究", 《测井技术》, vol. 38, no. 4, pages 384 - 390 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117454063A (zh) * | 2023-12-26 | 2024-01-26 | 西南石油大学 | 一种井筒油-气-水多相流流态判别及持水率计算方法 |
CN117454063B (zh) * | 2023-12-26 | 2024-03-12 | 西南石油大学 | 一种井筒油-气-水多相流流态判别及持水率计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN116341423B (zh) | 2023-08-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hasan et al. | A study of multiphase flow behavior in vertical wells | |
CN116341423B (zh) | 油水两相流滑移速度模型的计算方法 | |
Schmidt et al. | Severe slugging in offshore pipeline riser-pipe systems | |
Zhou et al. | A new model for predicting gas-well liquid loading | |
CN104895560A (zh) | 一种深水测试井筒压力、温度场模拟及水合物预测方法 | |
Xiaoming et al. | Flow pattern and pressure fluctuation of severe slugging in pipeline-riser system | |
CN112966418B (zh) | 一种海洋天然气水合物开采立管动力特性分析方法 | |
Wang et al. | Prediction of the critical gas velocity of liquid unloading in a horizontal gas well | |
Spedding et al. | Three phase oil-water-gas horizontal co-current flow: I. Experimental and regime map | |
CN114818535B (zh) | 一种基于流型转换界限的水平气井持液率计算方法 | |
NO20201135A1 (en) | Improved flow measurement | |
CN115345090A (zh) | 一种欠饱和煤层气储层压降漏斗动态传播的计算方法 | |
Duan et al. | Experimental and computational investigations on severe slugging in a catenary riser | |
CN111648759A (zh) | 一种页岩气水平井积液位置判断方法 | |
Schmidt et al. | Experimental study of two-phase normal slug flow in a pipeline-riser pipe system | |
Zhou et al. | Nonlinear adaptive observer for managed pressure drilling system | |
CN112796749B (zh) | 一种基于毛管数的油藏水驱数值模拟方法 | |
Cornish | The vertical multiphase flow of oil and gas at high rates | |
Capovilla et al. | Experimental investigation of vertical high-velocity two-phase flows in large-diameter pipes | |
Syikilili et al. | Effect of pipeline inclination on multiphase flow | |
Miyoshi et al. | Slug-catcher design for dynamic slugging in an offshore production facility | |
CN105160056A (zh) | 高温高压油气直井两相流射孔完井参数与产能优化方法 | |
Nagoo et al. | Will Gas Lifting in the Heel and Lateral Sections of Horizontal Wells Improve Lift Performance? The Multiphase Flow View | |
Skalle et al. | Experimental study of gas rise velocity and its effect on bottomhole pressure in a vertical well | |
CN117454063B (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 |