Mechanistic flow modelling in pipes - What difference does it make?

Report 47 Downloads 40 Views
Mechanistic flow modelling in pipes Pablo Adames September 11, 2012

Mechanistic flow modelling in pipes

Problem context

gathering systems trunk lines Network

tie ins

empirical correlations

Steady state software

mechanistic models Hydrodynamic models

inclination range

wells unied

90 to 0◦

-10 to +30◦ 90◦

Mechanistic flow modelling in pipes

The problem

What difference does it make to select one type of hydrodynamic model over the other?

Mechanistic flow modelling in pipes

Basic concepts To understand how gas-liquid flow models were developed: Flow patterns Slip and holdup Angle of inclination Model types: Analytical ⇔ mechanistic ⇔ empirical Fluid properties

Mechanistic flow modelling in pipes

Horizontal gas-liquid flow patterns

Mechanistic flow modelling in pipes

Flow pattern observation

In the following slides we will look at early flow pattern observations

Mechanistic flow modelling in pipes

Horizontal dispersed bubble UofC, circa 1970

Mechanistic flow modelling in pipes

Horizontal stratified wavy UofC, circa 1970

Mechanistic flow modelling in pipes

Horizontal slug flow UofC, circa 1970

Mechanistic flow modelling in pipes

Horizontal slug flow II Higher liquid loading, more frequent slugs

Mechanistic flow modelling in pipes

Horizontal annular mist UofC, circa 1970

Mechanistic flow modelling in pipes

Slip and holdup in multiphase pipe flow

Fundamental to understanding multiphase flow Holdup: phase fraction Slip: relative phase velocity

Holdup and slip change in a flowing system Their changes are interrelated There is no equivalent concept in single phase flow

Mechanistic flow modelling in pipes

Holdup as a fraction

Stratified flow Concept applicable to all flow patterns Gas phase flows on top Liquid flows in the bottom Area fraction of liquid, EL : L EL = AGA+A L

Mechanistic flow modelling in pipes

The Hydrodynamic Slip vslip = vG − vL The lighter phase will generally use energy more effectively to travel along faster. . . vslip > 0 But sometimes when descending. . . vslip < 0

Mechanistic flow modelling in pipes

Input and in situ phase fractions Why does phase fraction change in a pipe?

A recipe for phase fraction change: Ingredients 1 2 3 4 5 6

Inmiscible phases (they don’t blend) A pipeline Inertial forces (phase motion) Gravitational forces Dissipative forces (friction) Residence time

Hydrodynamic phase separation Input phase fraction changes downstream Each phase uses energy differently

Mechanistic flow modelling in pipes

Input vs. in situ phase fractions An highway analogy

INPUT SECTION: Input ratio of trucks to cars:

3 3

= 1.0

Input fraction of trucks to trafic:

3 3+3

= 0.5

in situ SECTION: in situ ratio of trucks to cars:

6 3.5

in situ fraction of trucks to trafic:

= 1.71 6 6+3.5

= 0.63

Mechanistic flow modelling in pipes

Input vs. in situ fraction Analysis and conclusion

Top view of highway: Road ≈ pipeline Trucks ≈ liquid Cars ≈ gas

The in situ fraction of the slower moving vehicle/fluid is greater than its input fraction There is hydrodynamic retention in steady state of the heavier phase This is known as the holdup or slip effect

Mechanistic flow modelling in pipes

The first flow pattern maps

Mandhane et al. horizontal

Aziz et al. vertical up

Return Source: Engineering Data Book, Gas Processors Suppliers Association, 2004. 12th Edition FPS, Tulsa, Oklahoma

Mechanistic flow modelling in pipes

Types of flow models

Analytical: built with first principles Empirical: built from observations alone Mechanistic: built from general laws and observations

Mechanistic flow modelling in pipes

Empirical versus mechanistic Why bother?

Mechanistic flow models Empirical flow correlations Developed by correlating dimensionless numbers

Developed from physical laws

Extrapolation is uncertain

Closed with empirical correlations

Interpolation can be very good

Reduced dependence on range of data

Can be simple to solve by hand

Usually computers needed to solve

Mechanistic flow modelling in pipes

Empirical versus mechanistic Why bother?

Mechanistic flow models Empirical flow correlations Developed by correlating dimensionless numbers

Developed from physical laws

Extrapolation is uncertain

Closed with empirical correlations

Interpolation can be very good

Reduced dependence on range of data

Can be simple to solve by hand

Usually computers needed to solve

Mechanistic flow modelling in pipes

Empirical versus mechanistic Why bother?

Mechanistic flow models Empirical flow correlations Developed by correlating dimensionless numbers

Developed from physical laws

Extrapolation is uncertain

Closed with empirical correlations

Interpolation can be very good

Reduced dependence on range of data

Can be simple to solve by hand

Usually computers needed to solve

Mechanistic flow modelling in pipes

History of the multiphase flow models

Source: Shippen, M., Steady-State Multiphase Flow -Past, Present, and Future, with a Perspective on Flow Assurance. Energy & Fuels, 2012. 26: p. 4145-4157.

Mechanistic flow modelling in pipes

Mechanistic flow pattern maps Angle of inclination dependency

Remember the Mandhane and Aziz et al flow pattern maps? First Flow Maps

What happens between horizontal and vertical? Only the mechanistic flow pattern maps answer that question

Let’s see an example using the Xiao et al Mechanistic model Between -10 and +10 degrees of inclination with the horizontal

Mechanistic flow modelling in pipes

The angle sensitivity of flow patterns

Mechanistic flow modelling in pipes

The operating line and the flow pattern map(s)?

Mechanistic flow modelling in pipes

Concluding on angle dependency

Would you use a single empirical flow pattern map for the whole pipeline again?

Mechanistic flow modelling in pipes

Flow model challenge What is the nature of the modelling challenge? Too many variables for a simple cause effect analysis

There are good models but the search for understanding is still on. . . Three-phase flow patterns Heavy oils Sand transport Non Newtonian rheologies Transient real time simulation Complex fluids

Mechanistic flow modelling in pipes

Flow model challenge What is the nature of the modelling challenge? Too many variables for a simple cause effect analysis

There are good models but the search for understanding is still on. . . Three-phase flow patterns Heavy oils Sand transport Non Newtonian rheologies Transient real time simulation Complex fluids

Mechanistic flow modelling in pipes

Flow model challenge What is the nature of the modelling challenge? Too many variables for a simple cause effect analysis

There are good models but the search for understanding is still on. . . Three-phase flow patterns Heavy oils Sand transport Non Newtonian rheologies Transient real time simulation Complex fluids

Mechanistic flow modelling in pipes

Flow model challenge What is the nature of the modelling challenge? Too many variables for a simple cause effect analysis

There are good models but the search for understanding is still on. . . Three-phase flow patterns Heavy oils Sand transport Non Newtonian rheologies Transient real time simulation Complex fluids

Mechanistic flow modelling in pipes

Flow model challenge What is the nature of the modelling challenge? Too many variables for a simple cause effect analysis

There are good models but the search for understanding is still on. . . Three-phase flow patterns Heavy oils Sand transport Non Newtonian rheologies Transient real time simulation Complex fluids

Mechanistic flow modelling in pipes

Flow model challenge What is the nature of the modelling challenge? Too many variables for a simple cause effect analysis

There are good models but the search for understanding is still on. . . Three-phase flow patterns Heavy oils Sand transport Non Newtonian rheologies Transient real time simulation Complex fluids

Mechanistic flow modelling in pipes

Effect of increasing water cut Inclination +0.25◦ , VSL = 0.05m/s, VsG = 1.3m/s, WC=1%

Mechanistic flow modelling in pipes

Effect of increasing water cut Inclination +0.25◦ , VSL = 0.05m/s, VsG = 1.3m/s, WC=5%

Mechanistic flow modelling in pipes

Effect of increasing water cut Inclination +0.25◦ , VSL = 0.05m/s, VsG = 1.3m/s, WC=10%

Mechanistic flow modelling in pipes

Effect of increasing inclination to +1◦ Inclination +1.0◦ , VSL = 0.1m/s, VsG = 1.4m/s, WC=10%

Mechanistic flow modelling in pipes

Effect of increasing the gas superficial velocity Inclination +1.0◦ , VSL = 0.1m/s, VsG = 1.7m/s, WC=10%

Mechanistic flow modelling in pipes

Effect of increasing the gas superficial velocity Inclination +1.0◦ , VSL = 0.1m/s, VsG = 2.5m/s, WC=10%

Mechanistic flow modelling in pipes

Main idea of unit cell model Fixed frame of reference

Mechanistic flow modelling in pipes

Main idea of unit cell model Moving frame of reference

Mechanistic flow modelling in pipes

A case Study

Now let us consider a real life case. . . Comparing published field measurements versus different model predictions.

Frigg to St. Fergus UK

Return http://www.uk.total.com/pdf/Library/PUBLICATIONS/Library-StFergusBrochure.pdf

Frigg to St. Fergus Slug catcher

http://www.uk.total.com/pdf/Library/PUBLICATIONS/Library-StFergusBrochure.pdf

Mechanistic flow modelling in pipes

Frigg to St. Fergus Description

d i = 774 mm (30.5 inch) Frigg to MCP01: 188.4km, 15 point elevation profile MCP01 to St. Fergus: 175km, 15 point elevation profile Detailed gas composition Specified variables: Pdowstream , Tupstream Calculated variable: Pupstream

Mechanistic flow modelling in pipes

Frigg to St. Fergus Description

d i = 774 mm (30.5 inch) Frigg to MCP01: 188.4km, 15 point elevation profile MCP01 to St. Fergus: 175km, 15 point elevation profile Detailed gas composition Specified variables: Pdowstream , Tupstream Calculated variable: Pupstream

Mechanistic flow modelling in pipes

Frigg to St. Fergus Description

d i = 774 mm (30.5 inch) Frigg to MCP01: 188.4km, 15 point elevation profile MCP01 to St. Fergus: 175km, 15 point elevation profile Detailed gas composition Specified variables: Pdowstream , Tupstream Calculated variable: Pupstream

Mechanistic flow modelling in pipes

Frigg to St. Fergus Description

d i = 774 mm (30.5 inch) Frigg to MCP01: 188.4km, 15 point elevation profile MCP01 to St. Fergus: 175km, 15 point elevation profile Detailed gas composition Specified variables: Pdowstream , Tupstream Calculated variable: Pupstream

Mechanistic flow modelling in pipes

Frigg to St. Fergus Description

d i = 774 mm (30.5 inch) Frigg to MCP01: 188.4km, 15 point elevation profile MCP01 to St. Fergus: 175km, 15 point elevation profile Detailed gas composition Specified variables: Pdowstream , Tupstream Calculated variable: Pupstream

Frigg to ST. Fergus Measured values

Group 1

Group 2

Group 3

Eq-gas MMsm3d 15.679 20.001 21.308 22.614 27.438 28.846 31.760 33.569 33.400 38.900 40.900 43.800 32.400 33.400 36.100 38.400 38.900 40.900 43.800

Measured values Pup Pdown bara bara 25.0 114.0 88.9 42.0 108.0 66.0 55.0 105.0 50.0 43.0 131.0 88.0 60.5 145.0 84.5 82.0 132.0 50.0 93.0 143.0 50.0 100.4 149.1 48.7 39.7 149.1 109.4 59.4 148.4 88.9 118.4 148.4 30.0 131.0 147.9 16.9 58.1 109.2 51.1 60.6 109.3 48.7 69.8 117.5 47.7 74.5 122.9 48.5 77.1 125.1 48.0 82.9 131.4 48.5 91.2 140.3 49.1

ΔP bar

Group 1: Frigg to St. Fergus; Group 2: Frigg to MCP01; Group 3: MCP01 to St. Fergus

Tup C 47.0 47.0 47.0 47.0 47.0 47.0 47.0 28.0 28.0 29.3 32.6 27.9 5.6 2.0 5.1 5.1 23.3 33.3 45.6

Field view

Results for pressure drop Relative error

Group 1

Group 2

Group 3

EatOli

B&B rev

B&B

-11.86 -1.59 -2.56 -5.59 -5.95 -1.55 -1.72 -1.40 -4.95 -5.17 -17.98 -11.60 -3.82 -1.82 -4.26 -3.44 -3.22 -4.51 -3.97 -5.10

0.93 -10.41 -12.80 -9.77 -18.85 -14.48 -14.94 -14.84 -22.83 -22.67 -28.28 -21.37 -16.56 -14.53 -16.59 -15.93 -15.42 -16.45 -15.70 -15.87

-19.85 -19.15 -19.28 -19.53 -24.80 -8.50 -19.35 -33.96 -28.36 -26.88 -27.27 -20.69 -18.45 -16.01 -17.74 -17.01 -16.45 -17.53 -16.90 -20.41

OliMec ei, % -10.26 2.46 1.62 -2.33 -3.30 1.50 1.07 1.19 -2.43 -2.48 -15.53 -9.24 -0.21 1.81 -0.97 -0.35 -0.11 -0.78 -1.33 -2.09

Xiao XiaoMod OLGAS2P 1.73 6.74 3.80 2.32 -0.66 2.36 1.50 0.30 -2.93 0.21 -16.88 1.30 -1.59 0.49 -2.40 -1.83 -1.66 -2.34 -2.76 -0.65

5.72 9.84 6.71 5.11 1.49 4.19 3.23 1.79 -1.17 1.56 -15.95 2.06 0.30 2.14 -0.97 -0.49 -0.37 -1.14 -1.77 1.17

-9.46 -0.16 -1.11 -4.42 -5.45 -0.82 -1.40 -1.20 -4.44 -5.51 -18.23 -11.76 -3.65 -1.66 -4.26 -3.71 -3.61 -4.27 -4.73 -4.73

Mechanistic flow modelling in pipes

Results for pressure drop Summary

ei |ei|

EatOli % -5.10 5.10

B&B rev % -15.87 15.96

B&B % -20.41 20.41

OliMec % -2.09 3.10

Xiao % -0.65 2.83

XiaoMod OLGAS2P % % 1.17 -4.73 3.47 4.73

Mechanistic flow modelling in pipes

Results for holdup Summary (only two measurements)

Measured Holdup 630 418

EatOli

Measured Holdup 630 418

EatOli

1624.0 1504.0

157.78 288.98

B&B rev 21271.0 18945.0

B&B rev 3276.35 4994.85

B&B 4953.0 4323.0

B&B

OliMec m3 308.0 294.0

OliMec ei, % 686.19 -51.11 1086.35 -26.23

Xiao 321.0 317.0

Xiao -49.05 -23.11

XiaoMod OLGAS2P 432.0 424.0

540.0 520.0

XiaoMod OLGAS2P -31.43 3.47

-14.29 29.34

Mechanistic flow modelling in pipes

Conclusions

1

The angle of inclination is very important for gas-liquid flow

2

Flow pattern maps provide insight into pipeline gas-liquid simulations

3

Mechanistic flow models are safer general purpose options

4

Mechanistic flow models are better holdup predictors

Thank you