CN110059411A - A kind of two-fluid slug flow calculation methodologies of coupling pressure equation - Google Patents
A kind of two-fluid slug flow calculation methodologies of coupling pressure equation Download PDFInfo
- Publication number
- CN110059411A CN110059411A CN201910316222.7A CN201910316222A CN110059411A CN 110059411 A CN110059411 A CN 110059411A CN 201910316222 A CN201910316222 A CN 201910316222A CN 110059411 A CN110059411 A CN 110059411A
- Authority
- CN
- China
- Prior art keywords
- phase
- fluid
- equation
- flow
- slug
- 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
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 103
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000008878 coupling Effects 0.000 title claims abstract description 17
- 238000010168 coupling process Methods 0.000 title claims abstract description 17
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 17
- 238000004364 calculation method Methods 0.000 title claims abstract description 15
- 239000007788 liquid Substances 0.000 claims abstract description 48
- 230000005514 two-phase flow Effects 0.000 claims abstract description 23
- 239000012071 phase Substances 0.000 claims description 93
- 230000004323 axial length Effects 0.000 claims description 12
- 239000007791 liquid phase Substances 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 6
- 238000004088 simulation Methods 0.000 abstract description 7
- 239000003208 petroleum Substances 0.000 abstract description 4
- 238000004519 manufacturing process Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- 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/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Software Systems (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention relates to a kind of two-fluid slug flow calculation methodologies of coupling pressure equation, it is characterised in that comprising steps of establishing the pressure equilibrium equation of compressible fluid according to the Conservation Relationship of pipeline deformation at any time and fluid flow through conduit velocity variations；Establish the two-fluid model of gas liquid two-phase flow slug flowing, wherein the two-fluid model includes continuity equation, the equation of momentum and the energy equation of gas liquid two-phase flow slug flowing；The equation group being made of pressure equilibrium equation, continuity equation, the equation of momentum and energy equation is solved, the state parameter that the slug flow of gas liquid two-phase flow changes over time is obtained.The present invention can carry out the accurate acquisition of gas liquid two-phase flow compressible fluid slug flow fluid plug frequency and liquid plug amount, be mainly used in petroleum industry production in solves the problems, such as slug flow brought by flow instability institute must progress slug flow field simulation calculating.
Description
Technical field
The present invention relates to a kind of two-fluid slug flow calculation methodologies of coupling pressure equation, are related to gas-liquid two-phase slug flow
Dynamic calculating field.
Background technique
Slug flow is gas-liquid two-phase flow pattern common in petroleum industry gas-liquid gathering line, because fluid flowing or landform etc. because
Occurs a kind of Liquid Flow kenel of gas-liquid gap flowing in the pipeline that element is induced.The slug flow that liquid plug and air mass are alternately present
Flowing in the duct can cause pressure in pipeline significant fluctuation occur, larger to the equipment impact of pipe outlet, so that
There are great unstability for fluid flowing in pipeline.Accurate simulation predicts slug flowing, realizes to liquid slug frequency and liquid plug
The accurate trapping of amount is of great significance to flow instability brought by slug flowing is solved in petroleum industry produces.
Meanwhile the accurately dynamic rule of simulation slug flow, and the important technology for improving process flow, designing slug catcher process equipment
It supports.
Slug flowing, because its gas-liquid two-phase interval flows in the duct, show to be alternately present in the duct it is single-phase with it is more
Phase, so that its transition simulation process is complicated.Currently, traditional slug flow model is complex, need to the liquid in slug unit
Plug establishes conservation equation with liquid film area respectively, is solved.In addition, traditional two-fluid model is merely by the equation of gas state
The problems such as update density relies on merely pressure Poisson's equation renewal speed, system suitability existing for Lai Jinhang pressure correction.
Summary of the invention
In view of the above-mentioned problems, the object of the present invention is to provide a kind of two-fluid slug stream calculation sides of coupling pressure equation
Method realizes accurate simulation and the prediction of gas liquid two-phase flow slug flow by establishing the two-fluid model of coupling pressure equation.
To achieve the above object, the present invention takes following technical scheme: a kind of two-fluid slug flow of coupling pressure equation
Calculation method, comprising the following steps:
According to the Conservation Relationship of pipeline deformation at any time and fluid flow through conduit velocity variations, compressible fluid is established
Pressure equilibrium equation；
Establish the two-fluid model of gas liquid two-phase flow slug flowing, wherein the two-fluid model includes gas-liquid two-phase
Flow slug flow dynamic continuity equation, the equation of momentum and energy equation；
The equation group being made of pressure equilibrium equation, continuity equation, the equation of momentum and energy equation is solved, gas-liquid is obtained
The state parameter that the slug flow of two-phase flow changes over time.
In above-described embodiment, the pressure equilibrium equation of compressible fluid:
Wherein, i is fluids within pipes phase number, αiFor the phase fraction of i fluid phase, ρiFor the density of i fluid phase, uiFor i stream
The speed of body phase, ciFor i fluid phase velocity of wave, p is pipeline pressure, and t is the time, and x is pipeline axial length.
In the various embodiments described above, the continuity equation of gas liquid two-phase flow slug flowing:
Wherein, i is fluids within pipes phase number, αiFor the phase fraction of i fluid phase, ρiFor the density of i fluid phase, uiFor i stream
The speed of body phase, t are the time, and x is pipeline axial length.
In the various embodiments described above, the equation of momentum of gas liquid two-phase flow slug flowing:
Wherein, i is fluids within pipes phase number, αiFor the phase fraction of i fluid phase, ρiFor the density of i fluid phase, uiFor i stream
The speed of body phase, p are pipeline pressure, pglThe pressure between gas-liquid two-phase interface, ГiFor the frictional resistance of i fluid phase and tube wall, ГglFor
Active force of the gas-liquid two-phase in pipeline axial direction, inclination angle of the θ between pipeline and horizontal line, g is acceleration of gravity, and t is the time, and x is
Pipeline axial length.
In the various embodiments described above, the energy equation of gas liquid two-phase flow slug flowing:
Wherein, i is fluids within pipes phase number, αiFor the phase fraction of i fluid phase, ρiFor the density of i fluid phase, uiFor i stream
The speed of body phase, p are pipeline pressure, EiFor i fluid phase specific internal energy, HiFor the ratio total enthalpy of i fluid phase, θ is between pipeline and horizontal line
Inclination angle, g is acceleration of gravity, and t is the time, and x is pipeline axial length.
In the various embodiments described above, the state parameter that the slug flow of gas liquid two-phase flow changes over time includes that slug flow causes
Pressure, liquid holdup, gas phase flow velocity or liquid phase flow rate.
The invention adopts the above technical scheme, which has the following advantages:
1, the present invention is by establishing coupling pressure equation with two-fluid model, it can be achieved that the standard of gas liquid two-phase flow slug flow
It really simulates and predicts, improve traditional two-fluid model and update density or simple dependence pressure Poisson by state equation merely
Equation renewal speed, Lai Jinhang pressure correction calculate existing system suitability problem；
2, the present invention is by establishing the two-fluid model of coupling pressure equation, the dynamic simulation meter of slug flow in Lai Jinhang pipeline
The free switching, it can be achieved that simulation of single-phase and two-phase flow is calculated, while precision height can be obtained and be easy to convergent analog result；
To sum up, the present invention can carry out the standard of gas liquid two-phase flow compressible fluid slug flow fluid plug frequency and liquid plug amount
Really acquisition is mainly used in and solves the problems, such as that the institute of flow instability brought by slug flowing must carry out in petroleum industry production
Slug flow field simulation calculate, to improve process flow, be aided with slug catching apparatus type selecting and design etc..
Detailed description of the invention
Fig. 1 is that the embodiment of the present invention is calculated the formation of pipeline hydraulic slug at any time and develops schematic diagram；
Fig. 2 is that t=40s pipeline phase fraction and pipeline pressure schematic diagram data is calculated in the embodiment of the present invention；
Fig. 3 is that t=40s pipeline gas/liquid phase velocity schematic diagram data is calculated in the embodiment of the present invention.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiments of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, ordinary skill people
Member's every other embodiment obtained without making creative work, shall fall within the protection scope of the present invention.
The two-fluid slug flow calculation methodologies of coupling pressure equation provided by the invention, include the following steps:
1, according to the Conservation Relationship of pipeline deformation at any time and fluid flow through conduit velocity variations, compressible fluid is established
Pressure equilibrium equation:
Wherein, i is fluids within pipes phase number (i=1,2), αiFor the phase fraction of i fluid phase, ρiFor the close of i fluid phase
Degree, uiFor the speed of i fluid phase, ciFor i fluid phase velocity of wave, p is pipeline pressure, and t is the time, and x is pipeline axial length.
2, on the basis of above-mentioned compressible fluid pressure equilibrium equation, couple gas-liquid two-phase fluid continuity equation,
The equation of momentum and energy equation establish the two-fluid model of gas liquid two-phase flow slug flowing, wherein two-fluid model includes connecting
Continuous property equation, the equation of momentum and energy equation, specifically:
Continuity equation:
Wherein, i is fluids within pipes phase number (i=1,2)；αiFor the phase fraction of i fluid phase, ρiFor the close of i fluid phase
Degree, uiFor the speed of i fluid phase, t is the time, and x is pipeline axial length；
The equation of momentum:
Wherein, i is fluids within pipes phase number (i=1,2), αiFor the phase fraction of i fluid phase, ρiFor the close of i fluid phase
Degree, uiFor the speed of i fluid phase, p is pipeline pressure, pglThe pressure between gas-liquid two-phase interface, ГiFor i fluid phase and tube wall
Frictional resistance, ГglFor gas-liquid two-phase pipeline axial direction active force, if the force direction is contrary with fluid flow rate before it
Symbol is "+", its preceding symbol is "-" if the force direction is identical as fluid flow rate direction, and θ is between pipeline and horizontal line
Inclination angle, g are acceleration of gravity, and t is the time, and x is pipeline axial length；
Energy equation:
Wherein, i is fluids within pipes phase number (i=1,2), αiFor the phase fraction of i fluid phase, ρiFor the close of i fluid phase
Degree, uiFor the speed of i fluid phase, p is pipeline pressure, EiFor i fluid phase specific internal energy, HiFor the ratio total enthalpy of i fluid phase, θ is pipeline
Inclination angle between horizontal line, g are acceleration of gravity, and t is the time, and x is pipeline axial length；
3, it can choose single order centered difference discrete equation and divide uniform grid, solve by pressure equilibrium equation, continuity
Changing over time for the slug flow of gas liquid two-phase flow finally can be obtained in the equation group of equation, the equation of momentum and energy equation
State parameter, including liquid holdup (Fig. 1), pressure (Fig. 2) and gas-liquid flow speed data (Fig. 3) etc..
That verifies the two-fluid slug flow calculation methodologies of coupling pressure equation of the invention below by specific embodiment can
With property and accuracy.
As shown in Figure 1, the increase of simulated time at any time, the slug energy of flow in pipeline is enough to be accurately captured.It can be in figure
It is apparent from the formation and development process of slug flowing, it was demonstrated that calculation method of the invention can capture in flow process
Slug flowing.After simulating 40s, full packages shows slug flow and stablizes the state occurred；
As shown in Fig. 2, pressure change embodies being remarkably decreased for manifold pressure caused by slug flow resistance；
As shown in figure 3, the variation of gas phase and liquid phase flow rate shows the feature of typical waterpower slug flow.
From the result of embodiment can be seen that the present invention can accurately capture in multiphase pipe flow the generation of slug flow and
Development process, the variation including parameters such as pressure caused by slug flow, liquid holdup, gas phase flow velocity, liquid phase flow rates, availability is high,
Accuracy is good.The parameter of the present embodiment is shown in Table 1, example boundary condition selection outlet level pressure 105Pa。
Table 1
| Variable | Numerical value | Variable | Numerical value |
| Caliber | 0.078m | Length | 37m |
| Gas superficial velocity | 0.2m/s | Liquid phase specific speed | 0.488 |
| Density of gas phase | 1.205kg/m3 | Density of liquid phase | 1000kg/m3 |
| Gas phase viscosity | 1.796×10-5Pa·s | Liquid phase viscosity | 1.39×10-3Pa·s |
| Initial time entrance gas phase content | 0.048 | Initial time entrance liquid phase content | 0.952 |
According to the disclosure and teachings of the above specification, those skilled in the art in the invention can also be to above-mentioned embodiment party
Formula carries out change and modification appropriate.Therefore, the invention is not limited to the specific embodiments disclosed and described above, to this
Some modifications and changes of invention should also be as falling into the scope of the claims of the present invention.In addition, although this specification
In use some specific terms, these terms are merely for convenience of description, does not limit the present invention in any way.
Claims (6)
1. a kind of two-fluid slug flow calculation methodologies of coupling pressure equation, it is characterised in that comprising steps of
According to the Conservation Relationship of pipeline deformation at any time and fluid flow through conduit velocity variations, the pressure of compressible fluid is established
Equilibrium equation；
Establish the two-fluid model of gas liquid two-phase flow slug flowing, wherein the two-fluid model includes gas liquid two-phase flow
Continuity equation, the equation of momentum and the energy equation of slug flowing；
The equation group being made of pressure equilibrium equation, continuity equation, the equation of momentum and energy equation is solved, gas-liquid two-phase is obtained
The state parameter that the slug flow of flowing changes over time.
2. the two-fluid slug flow calculation methodologies of coupling pressure equation according to claim 1, which is characterized in that compressible
The pressure equilibrium equation of fluid:
Wherein, i is fluids within pipes phase number, αiFor the phase fraction of i fluid phase, ρiFor the density of i fluid phase, uiFor i fluid phase
Speed, ciFor i fluid phase velocity of wave, p is pipeline pressure, and t is the time, and x is pipeline axial length.
3. the two-fluid slug flow calculation methodologies of coupling pressure equation according to claim 1, which is characterized in that gas-liquid two
Mutually flow the dynamic continuity equation of slug flow:
Wherein, i is fluids within pipes phase number, αiFor the phase fraction of i fluid phase, ρiFor the density of i fluid phase, uiFor i fluid phase
Speed, t is the time, and x is pipeline axial length.
4. according to claim 1 to the two-fluid slug flow calculation methodologies of 3 described in any item coupling pressure equations, feature exists
In the equation of momentum of gas liquid two-phase flow slug flowing:
Wherein, i is fluids within pipes phase number, αiFor the phase fraction of i fluid phase, ρiFor the density of i fluid phase, uiFor i fluid phase
Speed, p is pipeline pressure, pglThe pressure between gas-liquid two-phase interface, ГiFor the frictional resistance of i fluid phase and tube wall, ГglFor gas-liquid
Active force of the two-phase in pipeline axial direction, inclination angle of the θ between pipeline and horizontal line, g are acceleration of gravity, and t is the time, and x is pipeline
Axial length.
5. according to claim 1 to the two-fluid slug flow calculation methodologies of 3 described in any item coupling pressure equations, feature exists
In the energy equation of gas liquid two-phase flow slug flowing:
Wherein, i is fluids within pipes phase number, αiFor the phase fraction of i fluid phase, ρiFor the density of i fluid phase, uiFor i fluid phase
Speed, p is pipeline pressure, EiFor i fluid phase specific internal energy, HiFor the ratio total enthalpy of i fluid phase, θ inclining between pipeline and horizontal line
Angle, g are acceleration of gravity, and t is the time, and x is pipeline axial length.
6. according to claim 1 to the two-fluid slug flow calculation methodologies of 3 described in any item coupling pressure equations, feature exists
In the state parameter that the slug flow of gas liquid two-phase flow changes over time includes pressure caused by slug flow, liquid holdup, gas phase stream
Speed or liquid phase flow rate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910316222.7A CN110059411A (en) | 2019-04-19 | 2019-04-19 | A kind of two-fluid slug flow calculation methodologies of coupling pressure equation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910316222.7A CN110059411A (en) | 2019-04-19 | 2019-04-19 | A kind of two-fluid slug flow calculation methodologies of coupling pressure equation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110059411A true CN110059411A (en) | 2019-07-26 |
Family
ID=67319721
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910316222.7A Pending CN110059411A (en) | 2019-04-19 | 2019-04-19 | A kind of two-fluid slug flow calculation methodologies of coupling pressure equation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110059411A (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130317791A1 (en) * | 2012-04-26 | 2013-11-28 | Conocophillips Company | Hydrodynamic slug flow model |
| CN109635405A (en) * | 2018-12-05 | 2019-04-16 | 中国石油大学(北京) | A kind of Multiphase Flow transient state calculation method and system based on space conservation |
-
2019
- 2019-04-19 CN CN201910316222.7A patent/CN110059411A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130317791A1 (en) * | 2012-04-26 | 2013-11-28 | Conocophillips Company | Hydrodynamic slug flow model |
| CN109635405A (en) * | 2018-12-05 | 2019-04-16 | 中国石油大学(北京) | A kind of Multiphase Flow transient state calculation method and system based on space conservation |
Non-Patent Citations (1)
| Title |
|---|
| 赵志勇等: "气液两相段塞流中双流体模型的分析", 《管道技术与设备》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105302997B (en) | A kind of column separation-based on three-dimensional CFD makes the analogy method of water hammer up | |
| Newell Jr et al. | Analysis of combined free and forced convection for fully developed laminar flow in horizontal tubes | |
| CN106979808B (en) | Ultrasonic and target type flowmeter combined wet natural gas flow measuring method | |
| CN207937329U (en) | A kind of simulated sea bottom mixed transporting tube is influenced experimental provision by gas-liquid-solid multiphase flow flow velocity erosion | |
| Dang et al. | Experimental study on void fraction, pressure drop and flow regime analysis in a large ID piping system | |
| Taştan et al. | Effects of intake geometry on the occurrence of a free-surface vortex | |
| CN107122516A (en) | A kind of fire extinguishing system loses determination method from segment pipe along stroke pressure | |
| KR101988557B1 (en) | Apparatus for calculating pressure lost in pipeline and method thereof | |
| CN106650120A (en) | Measuring and calculating method for leakage rate of liquid hydrocarbon pipeline | |
| CN108197070A (en) | Natural gas not exclusively blocks pipeline method for numerical simulation | |
| CN114444343A (en) | Virtual height numerical simulation method for wet natural gas venturi tube flowmeter | |
| Vatani et al. | Experimental examination of gas-liquid two-phase flow patterns in an inclined rectangular channel with 90 bend for various vertical lengths | |
| Xu et al. | Experimental and simulation study of the effect of gravity on the solid-liquid two-phase flow and erosion of ball valve | |
| CN110348112A (en) | A kind of adaptive constrained procedure calculating natural gas vent pipeline maximum flow rate | |
| CN110059411A (en) | A kind of two-fluid slug flow calculation methodologies of coupling pressure equation | |
| CN114861568B (en) | Similar modeling method for spray evaporation two-phase flow process | |
| Wang et al. | The Effect of Spray-Sprinkler Design Parameters on Energy Consumption Using Computational Fluid Dynamics | |
| Irikura et al. | Numerical simulation of slugging of stagnant liquid at a V-shaped elbow in a pipeline | |
| CN107687876A (en) | A kind of measurement apparatus and measuring method for measuring fluid flow | |
| Guo et al. | Influence of offset terminal elbow shape on the discharge capacity of a high-rise building drainage system | |
| Lakehal et al. | Transition of gas-liquid stratified flow in oil transport pipes | |
| Ballil et al. | Investigation of turbulent airflow in a circular pipe using different orifice plate geometries | |
| CN217133883U (en) | Bernoulli demonstration experiment device | |
| Li et al. | An improved smoothed particle hydrodynamics method for modeling multiphase flows | |
| CN117993330B (en) | Curved surface river non-orthogonal water blocking normal form calculation method |
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: 20190726 |