Key information

Title: Method and tool for planning and dimensioning subsea pipelines for produced fluids
Application number: 20190663
Case type: National
Status: 15.03.2021 Patent meddelt (B1)
Filed: 24.05.2019
Effective date: 24.05.2019
Publicly available: 25.11.2020
Next annual fee due: 31.05.2027
Applicant: LEDAFLOW TECHNOLOGIES DA (NO)
Owner: LEDAFLOW TECHNOLOGIES DA (NO)
Inventor: Jørn Kjølaas (NO)
Agent: ACAPO ONSAGERS AS (NO)
Granted: 15.03.2021
Patent number: 345508
Expiry date: 24.05.2039

Abstract and drawing

Disclaimer: This text has been machine-scanned and may contain errors – please refer to "Publications" for legally binding content.

The invention relates to a computer implemented method and tool for determining pressure-drop in multiphase pipeline flow where the effective surface roughness of liquid film coated sections of the inner pipeline wall is assumed to be equal to the maximum hydraulic roughness. The maximum hydraulic roughness is further assumed to be proportional to a maximum stable droplet size, where K is a correlation coefficient. The invention further relates to applying the computer implemented method for designing a pipelinebased fluid transport system for transport of multiphase fluids.

Publications

Latest published version

Documents

Expand

Download files (34)

Date

Date	Doc. No.	Process	Case number	In/Out	Journal description	To/from
22.03.2021	15-01	Saksbehandling	20190663(345508)	UT	PT Registreringsbrev Nasjonal Patent (15) (PT20190663)
22.01.2021	14-01	Saksbehandling	20190663	UT	PT Meddelelse om patent	ONSAGERS AS
25.11.2020	13-05	Saksbehandling	20190663	INN	P26649NO00 Besvarelse 25 nov 2020	ONSAGERS AS
25.11.2020	13-04	Saksbehandling	20190663	INN	P26649NO00 Norske krav	ONSAGERS AS
25.11.2020	13-03	Saksbehandling	20190663	INN	P26649NO00 Endret fremstilling november 20 - renskrevet kopi	ONSAGERS AS
25.11.2020	13-02	Saksbehandling	20190663	INN	P26649NO00 Endret fremstilling november 20 - med markerte endringer	ONSAGERS AS
25.11.2020	13-01	Saksbehandling	20190663	INN	Korrespondanse (Hovedbrev inn)	ONSAGERS AS
13.10.2020	12-01	Saksbehandling	20190663	INN	Sammendrag
25.08.2020	11-01	Saksbehandling	20190663	UT	Realitet patent	ONSAGERS AS
20.05.2020	10-04	Saksbehandling	20190663	INN	Beskrivelse merket	ONSAGERS AS
20.05.2020	10-03	Saksbehandling	20190663	INN	Hovedbrev	ONSAGERS AS
20.05.2020	10-02	Saksbehandling	20190663	INN	Beskrivelse ren	ONSAGERS AS
20.05.2020	10-01	Saksbehandling	20190663	INN	Korrespondanse (Hovedbrev inn)	ONSAGERS AS
01.04.2020	08-02	Saksbehandling	20190663	UT	PT20190663_Priority
01.04.2020	08-01	Saksbehandling	20190663	UT	Brev prioritetsdokument	ONSAGERS AS
30.03.2020	09-01	Saksbehandling	20190663	INN	Korrespondanse (Hovedbrev inn)	ONSAGERS AS
06.01.2020	07-04	Saksbehandling	20190663	UT	mothold D1
06.01.2020	07-03	Saksbehandling	20190663	UT	mothold D2
06.01.2020	07-02	Saksbehandling	20190663	UT	mothold D6
06.01.2020	07-01	Saksbehandling	20190663	UT	Generelt Brev Patent	ONSAGERS AS
12.12.2019	06-03	Saksbehandling	20190663	INN	02-PT_report 112826	ONSAGERS AS
12.12.2019	06-02	Saksbehandling	20190663	INN	01-Realitet patent	ONSAGERS AS
12.12.2019	06-01	Saksbehandling	20190663	INN	Korrespondanse (Hovedbrev inn)	ONSAGERS AS
12.12.2019	05-02	Saksbehandling	20190663	UT	PT_report 11:28:26
12.12.2019	05-01	Saksbehandling	20190663	UT	Realitet patent	ONSAGERS AS
17.07.2019	04-01	Saksbehandling	20190663	UT	Infobrev til oppfinner	ONSAGERS AS
17.07.2019	03-04	Saksbehandling	20190663	INN	Erklæring Søkers rett til oppfinnelse	ONSAGERS AS
17.07.2019	03-03	Saksbehandling	20190663	INN	Hovedbrev	ONSAGERS AS
17.07.2019	03-02	Saksbehandling	20190663	INN	Fullmakt	ONSAGERS AS
17.07.2019	03-01	Saksbehandling	20190663	INN	Korrespondanse (Hovedbrev inn)	ONSAGERS AS
03.06.2019	02-01	Saksbehandling	20190663	UT	Formalia 1	ONSAGERS AS
24.05.2019	01-03	Saksbehandling	20190663	INN	P26649NO00 Figures	ONSAGERS AS
24.05.2019	01-02	Saksbehandling	20190663	INN	Beskrivelse	ONSAGERS AS
24.05.2019	01-01	Saksbehandling	20190663	INN	Søknadsskjema Patent	ONSAGERS AS

Description, claims and drawings

Disclaimer: This text has been machine-scanned and may contain errors – please refer to "Publications" for legally binding content.

Description

Method and tool for planning and dimensioning subsea pipelines for produced fluids

Field of invention

The invention relates to a computer implemented method and tool for determining pressure-drop in multiphase pipeline flow, and further to a computer implemented method for designing a pipeline-based fluid transport system for transport of multiphase fluids. The invention is suited for planning and dimensioning subsea pipelines transporting produced hydrocarbon fluids from oil and gas reservoirs.

Background

The capital expenses associated with developing and installing production and transportation facilities required to extract oil and gas from reserves in reservoirs beneath a seabed often requires reserves containing relatively huge volumes of extractable reserves to make the project profitable at a competitive oil and gas price. Additional and more marginal fields/discoveries often cannot justify these capital costs.

However, nearby a major field/discovery there are often discovered more economic marginal volumes. Such additional reserves may have their capital costs for developing the reserves significantly reduced by connecting their production wells to the production and/or transportation facilities developed for the major reserve. As a result, it is often laid out a web of subsea pipelines (may be denoted as tie-backs in the literature) on the seabed connecting production wells of several additional fields/discoveries in a region to a common riser of major fields/discoveries for bringing the produced hydrocarbon fluids to an offshore production facility on the surface, or alternatively to a hub or template for passing the produced hydrocarbons to a transportation pipeline to an onshore production facility, etc. Typically, the system contains one or more well(s) from one or several templates to the reservoir(s) and pipelines from the templates to the topside/onshore production facilities (e.g. platforms, separators, installation), such as illustrated schematically in figure 6 (the figure is a facsimile of figure 14 of US 2010/0059221).

The number of pipelines and pipeline lengths typically involved in development of offshore oil and gas fields makes the design, more precisely optimization of the pipeline design, an essential part of the development of oil fields. There may be considerable savings in terms of material savings, reduced need for energy to transport the produced fluids, increased transportation capacity of produced fluids, and capital costs associated with the instalment of the pipeline by finding the optimized pipeline path (where to lay the pipeline in the subsea terrain) and pipeline dimensions (diameter, pump capacity, valves and so on). One key parameter in this regard is the prediction of the pressure drop of the produced fluid flow over various segments of the pipeline(s) being designed.

The produced fluids in oil and gas extraction behave typically as multiphase fluid flows which may take many different regimes such as slug flow, bubbly flow, stratified flow, annular flow, and/or churn flow, and involving complex irregular interactions between the flowing phases inducing pressure drops, deposits, liquid accumulation, unstable flows, etc. This complexity and irregularity of the multiphase produced fluid requires knowledge of the expected fluid behaviour along the entire length of the pipeline(s) being designed.

Such prior knowledge can only be attained by computational fluid dynamics.

Ideally, these computational fluid dynamical calculations should be performed with three-dimensional a priori mathematical descriptions of the laws of nature governing the multiphase flow to obtain a satisfactorily robust and reliable prediction of the fluid behaviour along the pipeline. However, the pipeline lengths (may be tens of kilometres long) and number of pipelines typically involved in the development of oil and gas exploration make such “full blooded” computational fluid dynamical models far too demanding in computing power to be practically applicable for decision processes in oil and gas industry.

Prior art

One approach for considerably reducing the need for computing power to solve the numerical model in computational fluid dynamics is approximating the involved flow physics such that the mathematical description may be given in less than three space dimensions. For example, assuming the multiphase flow being fully developed (having stable flow volumes of the continuous fluid phases) in a pipeline segment enables forming a one-dimensional mathematical description of the laws of nature governing the multiphase fluid flow in the pipeline segment.

An example of such approach is known from document EP 2926 278 which teaches a method for determining flow parameters of a multiphase flow in a pipeline segment by feeding in the pipeline inclination (relative to the earth gravitational field), pipeline inner diameter and the flow volume rates of the continuous phases of the multiphase flow (typically oil, water and natural gas), and then run a computer implemented one-dimensional mathematical/numerical model based on Eulerformulated transport equations for conservation of mass, momentum, turbulence fields, energy, etc. for all fluid phases and fields of the multiphase flow. The mathematical/numerical model is closed by applying standard wall function relations to estimate the fluxes of mass, momentum and energy etc. as boundary conditions at the fluid phase interfaces. Another example of one-dimensional approach for predicting multiphase fluid flow behaviour is the commercially available computer-implemented simulation tool sold under the trademark OLGA® by Schlumberger.

However, for gas-dominated flows, it is observed that the pressure-drop may be significantly under-predicted by one-dimensional computational fluid dynamic calculations. An example of such observation is presented in figure 1. The figure shows predicted pressure drop for three-phase gas-dominated pipeline flows compared with experimental measurements presented in Kjølaas [1] and in Kjølaas et al. [2]. The prediction of the pressure drop was made with a commercially available one-dimensional computational fluid dynamic program sold under the trademark “LedaFlow 2.4”. Figure 1 shows measured pressure-drops for a threephase flow having water cuts (marked by “WC” in the figure) of 0, 25, 50, and 100 %, respectively, as function of the gas superficial velocity (Marked as “USG” on the figure) in an 8-inch pipeline inclined at 5 ° (relative to the horizontal) and fed with gas, oil and water at a pressure of 60 bar. As seen on the figure, especially for high water cuts and high superficial gas velocities, the pressure drop may be underpredicted by as much as 20 – 30 %, which is a highly unfortunate inaccuracy for dimensioning pipelines for wet gas transport.

Laurinat et al. [3] describes a rigorous model of the phenomenon of the wall film, where the liquid distribution is modelled along the pipe perimeter by writing massand momentum conservation equations for the circumferential transport of the liquid along the wall. They concluded that droplet mass transfer through atomization and deposition was not sufficient to counteract gravitational drainage of the film. They thus argued that asymmetries in the film surface roughness cause secondary flow recirculation and give rise to circumferential interfacial stresses, which become significant for thin films near the top of the pipe. A similar approach is known from Bonizzi & Andreussi [4], who were able to obtain good agreement with measured film distributions by accounting for wave spreading in addition to droplet deposition and gravity.

Document WO 2018/191736 discloses a method for determining the pressure drop in low-liquid load multiphase pipeline flows which includes applying the viscosity of the liquid to estimate the effective roughness of a liquid film on the interior wall of the pipeline and then applying the effective roughness of the liquid film to determine the pressure-drop. The method is described as suited for annular flows in substantially vertical pipes at high gas rates.

Wallis [5] and Belt and Portela [6], both proposed that the interfacial friction coefficient for thin films can be represented by a relation determining the interfacial friction coefficient as a function of the ratio of liquid film thickness and pipeline diameter. Asali & Hanratty [7] proposed a somewhat similar model, but one that scales very differently with the gas density and pipe diameter. Van Rossum et al. [8] proposed that the onset of entrainment for low film Reynolds numbers (outside the disturbance wave regime) was well correlated by a critical Weber number. Biberg et al. [9] has proposed a model where they assumed that the hydraulic roughness kS equals a constant multiplied by the film thickness δ.

The above propositions of Wallis [5], Belt et al. [6], Asali & Hanratty [7], van Rossum et al. [8], and Biberg et al. [9] have been tested experimentally by the present inventor and found to have limited success with respect to the modelling of thin films. Specifically, the idea that the film thickness is given by some Reynolds number does not work for high liquid viscosities, because the film thickness then becomes unreasonably large.

There is a need for an improved approach to determine the pressure-drop in multiphase pipeline flow, and in particular for gas-dominated multiphase pipeline flows.

Objective of the invention

The main objective of the invention is to provide a computer implemented method and tool for determining the pressure-drop in non-vertical pipeline-based transport systems for multiphase fluids, and in particular transport systems for gas-dominated multiphase fluids.

Another objective of the invention is to provide a computer implemented method and tool for determining pipeline dimensions and penstock/pipeline paths for a nonvertical pipeline-based transport system for multiphase fluids, and in particular extracted wet gas fluids from a production well to a production facility.

Description of the invention

Experiments made by the present inventor and theoretical considerations indicate that the problem with under-prediction of the pressure-drop in gas-dominated pipeline flows is due to liquid at the bottom of the pipeline being “torn off” into liquid droplets and then entrained in and carried along with the relatively rapidly flowing gas phase. These liquid droplets deposit on the pipeline wall somewhat further downstream and is transformed into a relatively slow-moving film of liquid on the pipeline wall. The interaction between the liquid film and the flowing gas makes the liquid film surface considerably rougher than the surface roughness of the (dry) pipeline wall. Without being bound by theory, it is assumed that it is this relatively large surface roughness at liquid coated pipeline wall segments as compared to dry pipeline wall which causes the under-prediction of the pressuredrop. Figure 2 is a drawing schematically illustrating a cross-sectional cut view of a non-vertical pipeline as seen along its centre axis. A “pond” (5) of liquid at the lower part of the pipeline (1) circumference is releasing droplets (4) which are entrained in the gas flow (3) and deposited in the pipeline wall (1) and forming a liquid film (2).

The present invention utilises an approach for improving the pressure-drop prediction in the multiphase fluid flow which:

a) does not try to resolve the full details of the liquid film on the pipeline wall, but instead assumes that the wall film thickness is primarily driven by droplet deposition and entrainment,

b) assumes the action of gravity and secondary flows can be ignored, i.e. the present invention assumes it is unnecessary to have a detailed representation of the wall film, it is sufficient to be able to predict the average surface roughness of the liquid film and determining the frictional shear stress therefrom, and

c) assumes the average effective surface roughness of the liquid film is equal to the hydraulic roughness, ks, which can be determined by the following postulates:

1) The hydraulic roughness, ks of a thin liquid film scales linearly with the size of the ripples, dripple, on the film:

2) When the ripples on the film become sufficiently large, they are destroyed/ atomized by the turbulent fluctuations in the gas which limits the ripple size to the value

3) where is the maximum stable droplet size as determined by a droplet model, and

4) the gas flow velocity is usually sufficiently large to effectively maintain the size of the ripples of the liquid film at

The above assumptions may be summarised as: The effective surface roughness, keff, of liquid film coated sections of the inner pipeline wall is assumed to be equal to the maximum hydraulic roughness, which is assumed to be proportional to a maximum stable droplet size, '.

where K is a correlation coefficient.

Thus, in a first aspect, the invention relates to a computer implemented method for determining pressure-drop in a multiphase flow in a non-vertical pipeline having an inner wall facing the multiphase flow, where multiphase flow at least comprises a continuous gas phase forming a gas/inner wall interface and a continuous liquid phase, and

wherein the method comprises:

applying a computational fluid dynamic model describing the multiphase pipeline flow, and

solving the computational fluid dynamic model to determine the pressuredrop in the multiphase flow,

characterized in that the computational fluid dynamic model is adapted, for the gas/inner wall interface, to:

applying an effective surface roughness, keff ; determined by the relation: keff = , where K is a correlation coefficient and is a maximum stable droplet size as determined by a droplet model, and

applying the effective surface roughness, keff, to determine a shear-stress, τg, and then applying the shear-stress, τg, as a boundary condition determining the gas flow resistance at the gas/inner wall interface.

The method of the invention according to the first aspect applies a relatively coarse, but still often acceptably accurate, approximation of the effective surface roughness, keff, by assuming that the entire surface area of the inner wall of the pipeline in contact with the continuous gas phase is covered by a liquid film.

However, for some applications, such as e.g. high gas load flows, the liquid droplet generation and subsequent deposition of liquid droplets may not be sufficient to form a liquid film on the entire inner wall being in contact with the continuous gas phase. A more accurate approximation will thus be obtained by applying a computational fluid dynamical model which is adapted to determine which parts of the gas/inner wall interface which are covered by a liquid film, and which parts of the gas/inner wall interface which are not covered by a liquid film.

Thus, in a second aspect, the invention relates to a computer implemented method for determining pressure-drop in a multiphase flow in a non-vertical pipeline having an inner wall facing the multiphase flow, where multiphase flow at least comprises a continuous gas phase forming a gas/inner wall interface and a continuous liquid phase, and

wherein the method comprises:

applying a computational fluid dynamic model describing the multiphase pipeline flow, and

solving the computational fluid dynamic model to determine the pressuredrop in the multiphase flow,

characterized in that the computational fluid dynamic model is adapted, for the gas/inner wall interface, to:

determining a first fraction of the gas/inner wall interface considered being covered by a liquid film and a second fraction of the gas/inner wall interface considered not being covered by a liquid film,

for the first fraction; applying an effective surface roughness, keffi, determined by the relation: , where K is a correlation coefficient and a maximum stable droplet size as determined by a droplet model, and applying the effective surface roughness, keffi, to determine a shearstress, τg l, and

for the second fraction; applying an effective surface roughness, keff 2, defined to be equal to a surface roughness of a dry inner wall of the pipeline, and applying the effective surface roughness, keff,2, to determine a shear-stress, τg 2, and applying the shear-stresses, τg1 and τg2 as boundary conditions determining the gas flow resistance at the gas/inner wall interface.

As used herein, the term “multiphase flow” is a simultaneous flow of one or more continuous liquid phases and at least one continuous gaseous phase. The multiphase flow may exist as a homogeneous mixture or as independently flowing fluid phases which may attain complex interactions causing flow phenomena such as slug flow, mists, emulsions, or bubble flow. The pressure-drop in multiphase flows may be significantly higher than pressure-drops found in single-phase flows.

The term “multiphase flow in a non-vertical pipeline” as used herein encompasses any gas and liquid containing multiphase flow in a pipeline where a liquid phase may be expected to flow along the lower part (relative to the earth gravitational field) of the inner wall of the pipeline. I.e. multiphase flows in pipelines where the pipeline does not ascend/descend too steep to enable forming a substantially continuous liquid phase at the lower part of the pipeline circumference. In practice, the term “non-vertical” encompasses pipeline flows where the pipeline inclination angle (towards the horizon) is in the range of from - 75 ° to 75 °, where a positive inclination angle corresponds to the pipeline ascending in the flow direction, a negative inclination angle corresponds to the pipeline descending in the flow direction, and inclination 0 ° corresponds to the pipeline being horizontal relative to the earth gravitation field. The pipeline encompasses any form of tube/pipe being applied to transport or move a multiphase fluid flow, such as e.g. tie-backs, oil and gas production bore tubes, transport pipes in pharmaceutical and/or chemical industrial process facilities etc.

The term “computational fluid dynamic model” as used herein encompasses any computer implemented numerical model able to solve the laws of nature governing the multiphase fluid flow inside a pipeline, and which enables determining which part of the surface of the inner wall of the pipeline which is in contact with the at least one continuous gas phase of the multiphase flow. The invention according to the first and second aspect may apply any known computational fluid dynamic model having the above-given properties which applies a three-, two- or onedimensional mathematical description of the laws of nature governing the fluid flow.

The at least one continuous gas phase may, and often will, contain a multitude of entrained liquid droplets torn of the at least one continuous liquid phase (typically being present at the bottom of the pipeline). Some of these entrained liquid droplets will, depending on the gas flow velocity and other flow characteristics, deposit on the gas/inner wall interface. For surface areas of the inner wall exposed to relatively heavy liquid droplet deposition, the deposited liquid droplets will form a liquid film covering this part of the inner wall. However, for surface areas of the inner wall not exposed to relatively heavy liquid droplet deposition, there will be no effective liquid film formation, and, in some cases, the inner pipeline wall may even be dry. Thus, the term “gas/inner wall interface” as used herein encompasses both surface area of the inner pipeline wall facing a continuous gas phase being covered by a liquid film and surface area of the inner pipeline wall facing a continuous gas phase not being covered by a liquid film.

The invention according to the first and second aspect applies to all multiphase flows in pipelines involving a gas flow as a separate phase (and not solely as small dispersed bubbles inside the liquid). The invention according to the first and second aspect is especially suited for predicting the pressure-drop in pipeline-based transport systems where the gas volume fraction of the multiphase flow, measured as the ratio of the flow volume rate of the continuous gas phase over the sum of the flow volume rate of all present continuous liquid phases, is higher than 50 %, preferably higher than 75 %, and most preferably higher than 90 %. Such multiphase flows have typically a relatively large continuous gas phase occupying a significant volume (the gas zone) of the pipeline such that the flow resistance at the gas – inner wall interface may significantly influence the pressure-drop of the multiphase flow. Especially if the gas – inner wall interface is covered by a liquid film.

An example embodiment of a computational fluid dynamic model suited for being used in the first and the second aspect of the invention is a "hybrid model" approach, which considers the multiphase fluid to consist of two or more continuous fluid phases being separated by large-scale interfaces and further consisting of one or more dispersed fluid phases such as e.g. bubbles/droplets of one fluid suspended in a continuous phase of another fluid. The hybrid model approach considers each of the continuous and each of the dispersed fluid phases in the multiphase flow to be a separate fluid field. The velocities of the continuous fluid fields are obtained by solving the respective momentum equations, while the velocities of the dispersed fields are calculated from empirical "slip relations", which are expressions that describe the velocity difference between the dispersed fields and the associated continuous field. The simplest version of a slip relation is to simply assume that the velocity difference is zero, meaning that the dispersed field has the same velocity as the associated continuous field. The purpose of using slip relations instead of momentum equations for the dispersed fields is to reduce the computational cost of solving the equations, thus reducing the execution time for simulations.

An example embodiment of a hybrid computational fluid dynamic model suited for use in the first and the second aspect of the invention may advantageously describe the multiphase pipeline flow as a set of fluid fields representing:

- a set of continuous fluid phases being separated by large-scale interfaces, at least comprising a continuous liquid phase and a continuous gas phase being separated by a large-scale gas/liquid interface, and

- a set of dispersed phases, at least comprising a dispersed phase of liquid droplets entrained in the continuous gas phase,

and wherein

- each continuous fluid field is described by Eulerian formulated volume and ensemble averaged turbulent transport equations which are solved for each continuous fluid field, and

- each of the dispersed phases is described by "slip relations", which are expressions that describe the velocity difference between the dispersed fields and the associated continuous field.

The hybrid computational fluid dynamic model enables determining the production rate of liquid droplets being torn-off from a continuous liquid phase and being entrained in the continuous gas phase. This information may, in one example embodiment, be applied to determine which parts of the gas/wall interface being covered by a liquid film and which parts of the gas/inner wall interface not being covered by a liquid film. Hybrid computational fluid dynamic models are inherently “made ready” for determining a liquid droplet concentration profile,

the continuous gas phase as a function of vertical distance above the large-scale interface separating the continuous liquid phase and the continuous gas phase, and is thus especially suited for use in the second aspect of the invention.

If the pipeline section is non-vertical, which in practice means the pipeline has an inclination angle in the range of from -75° to 75°, where 90° is vertical relative to the earth gravitational field, there will usually be a continuous liquid phase at the lower part of the pipeline along the length of the section which interacts with the flowing continuous gas phase such that liquid droplets are being formed, entrained by and carried along with the relatively fast flowing continuous gas phase. The resulting liquid droplets in the continuous gas phase will gradually hit and be adsorbed by the inner wall and/or hit and fall back and be absorbed by the continuous liquid phase at the bottom of the pipeline and thus lowering the liquid droplet concentration with increasing height (inside the pipeline) above the liquid phase. The liquid droplets hitting the inner wall at the gas/wall interface will wet the inner wall and may with sufficiently high droplet concentration form a film of liquid at the gas/wall interface covering the entire inner wall. The invention according to the first aspect of the invention is suited for such cases. On the other hand, if the liquid droplet concentration in the continuous gas phase is too low to form a liquid film covering the entire inner wall at the gas/wall interface, there will typically be a region of the inner wall (at the gas/wall interface) not having a liquid film and a region being covered by a liquid film. The invention according to second aspect is suited for such cases but may of course also be applied for cases where the entire inner wall at the gas/wall interface is covered by a liquid film. In the physical world, there will usually be a gradual transition from completely wetted (i.e. covered by a liquid film) inner wall to a dry wall. However, in practice, the transition from wet to dry inner wall may be well approximated by introducing a critical droplet concentration criterium determining which part of the inner wall at the gas/inner wall interface is covered by a liquid film and which part is to be considered being dry.

Thus, in one example embodiment, the invention according to the second aspect further comprises applying a hybrid computational fluid dynamic model adapted to determine a liquid droplet concentration profile, as a function of vertical distance above a large-scale interface separating a continuous liquid phase and a continuous gas phase and which further assumes that the inner wall of the pipeline is covered by a liquid film when g

cut-off threshold, and that the inner wall of the pipeline is dry when

Analyses of experimental flow data indicate that a reasonable value for is around 10'5 for hydrocarbon liquids, while a significantly lower value of around 10'9 should be applied for aqueous liquids. The reason for this variation is believed to be caused by differences in the attraction forces between the wall and the liquid. Specifically, aqueous liquids are observed to very easily stick to the steel walls of the pipe, and this adhesive behaviour facilitates build-up of a wall film at lower droplet concentrations than hydrocarbon liquid, which does not exhibit this behaviour.

An example embodiment for determination of the critical droplet concentration

suitable for use in the invention according to second aspect is given by the following expression:

where σgl is the gas-liquid surface tension.

The term “a surface roughness of a dry inner wall of the pipeline” as used herein, is the surface roughness of the inner wall of the pipeline in dry condition, i.e. without a liquid film or deposited liquid droplets. The surface roughness of a (dry) pipeline wall is a material property of material of the pipeline, which usually is a steel.

The term “correlation coefficient K” as used herein, is a coefficient introduced to make the assumption that the average effective surface roughness of liquid-film covered gas/inner wall interface is proportional to a maximum stable droplet size

(of the same liquid as constituting the liquid film). The correlation

coefficient K can be determined experimentally by applying experimentally measured pressure-drops for various multiphase flows in a pipeline segment and then applying a computational fluid dynamical to model the same multiphase flows in the pipeline segment to tune which effective wall roughness at the gas/inner wall interface which gives a match between the measured and calculated pressure-drop.

The term “a predetermined maximum droplet size, ” as used herein is a maximum stable droplet size as determined from a droplet model which typically balances the destructive forces of turbulence on the droplet and the stabilising force of the liquid of the droplet. There are known several liquid droplet models in the literature, all of which may be applied in the present invention according to the first and second aspect.

An example embodiment of a droplet model suitable for use in the method according to the first aspect of the invention, is a model developed by Arai et al.

[11]:

where σ* = σgl (1 k ■ Ca ), where σgl is the surface tension between the liquid and the gas, k is a correlation coefficient, and Ca is a capillary number defined to be:

where μd is the dynamic viscosity of the liquid droplet, pg is volumetric density of the continuous gas phase, and ε is the energy dissipation rate per unit mass given by the relation:

where Sg is the gas perimeter, τg is the gas-wall shear stress, ug is the gas velocity, Si is the interface perimeter, τί is the interfacial shear stress, Δu is the gas-liquid slip velocity, and Dhg is the hydraulic diameter for the continuous gas phase. The σ*, or more precisely, the content of the parenthesis, (1 k ■ Ca), may be considered being a viscosity correction term.

This model contains the maximum droplet size on both sides of the equality sign, requiring an iterative method to obtain a value. To simplify the calculation, one may approximate the droplet size in equation (4) using the droplet size model derived by Hinze [10]: dHinze = 0.725 so that equation (4) becomes:

The Hinze droplet size model (dHinze) is equivalent to the Arai model (equation (3)) for the limiting case of zero viscosity. By applying the Hinze model in equation (4),

can be calculated explicitly, without iterations.

Arai et al. [11] pointed out that the maximum droplet size did not increase indefinitely with the viscosity in their stir tank experiments, like their model suggests. Instead it seemed to level off at a certain point. Arai et al. reasoned that the cause of this was that when the viscosity became large enough, re-establishing the spherical form after deformation would be delayed by the viscous forces, allowing the turbulence to further deform already deformed droplets, facilitating the break-up process. Arai et al. did however not attempt to model this phenomenon.

Figure 3 is a diagram showing the ratio between the surface roughness and dHinze plotted against the capillary number Ca. The black dotted line is to the left of the breaking point marked cut-off according to the droplet model of Arai et al. when the correlation coefficient k = 33, i.e. which corresponds to eqn. (7) below. The experimental data in figure 3 indicates that above a certain capillary number, the effective roughness stops increasing, which is qualitatively consistent with the observations by Arai et al. This effect is taken into account in the present example embodiment for droplet model by introducing a cut-off in the viscosity correction term equal to a value of 6.6. The figure contains some experimental data points (at the upper right of the figure) which exceeds this cut-off, but these are ignored in order to keep the droplet model simple. With these corrections, σ* may be given with a viscosity correction term as follows:

The viscosity correction by Arai et al. including the present introduced cut-off criterium and experimentally determined correlation coefficient, k , may be introduced into eqn. (1) to give:

The viscosity correction may in other words be expressed by using the Hinze droplet model to determine the maximum droplet size and expressing the correlation coefficient, K, in the assumption summarised in eqn. (1) as a relation containing the capillary number, Ca. The above expression, eqn. (7), for the effective surface roughness has been shown to function well in correcting computational fluid dynamic predictions of pressure-drop in multiphase pipeline flows for two-phase flows.

For three phase flows, the relatively rapid flow rate of the gas phase usually effectively mixes the two liquid phases and forms an emulsion-like state where one liquid is suspended and dispersed as droplets in the other. This is expected to affect the viscosity and the surface tension of the liquid (emulsion). It is also observed that when such suspension forms, the viscosity effect on the effective surface roughness does not level-off. The reason for this observation may be that the dynamic behaviour of liquids with a dispersion-enhanced viscosity differs from that of pure liquids. To account for this observation, the present inventor proposes to treat the effect of the pure viscosity and the dispersion effect separately, only imposing a limit on the effect of the pure liquid viscosity. This leads to the following expression for the effective wall roughness ks:

where is the effective dynamic viscosity of the liquid mixture and is the dynamic viscosity of the major continuous liquid phase, an s the capillary number Ca calculated by using

The dynamic viscosity of a liquid mixture, may be determined by using an emulsion model. Different fluid systems display different behaviours in this respect, and the emulsion model must be selected accordingly. One such model was developed by Brinkman [12] which sOtates that:

where φ is the intrinsic volumetric fraction of liquid droplets suspended in the major continuous liquid phases. In the special case where one liquid (with volume fraction ai) is completely dispersed into another liquid (with volume fraction 012), φ is given by:

Computational fluid dynamic models describing multiphase flows involving liquids and gases typically determines the droplet concentration of liquid suspensions and may thus be adapted to incorporate the emulsion model of Brinkman to determine the effective dynamic viscosity of the mixed liquid film.

The gas-liquid surface tension of the liquid film formed by the suspension is for simplicity assumed to be a simple arithmetic average of the surface tension values for both liquid phases using the respective average volumetric phase fractions on and 012 as weighting parameters:

The average phase fractions ai and 012 are parameters typically being determined by computational fluid dynamic models able to handle multiphase fluid flows.

An example embodiment suited for use in the method according to the first invention for converting the average effective surface roughness to the shear stress at the gas-wall interface in the gas zone is applying a friction factor model that accounts for the wall roughness, such as e.g. the Håland model [13] to determine a gas friction coefficient fg from the effective surface roughness, by the relation:

Here Reg is the gas Reynolds number and Dhg is the hydraulic diameter for the continuous gas phase. The gas Reynolds number is defined as:

Here, pg is the gas density, ug is the gas velocity, and pg is the gas viscosity. The hydraulic diameter for the continuous gas phase Dhg may be calculated as:

where αg is the volume fraction of the gas zone, D is the pipe diameter, and Sg is the length of the gas-wall perimeter.

And then apply the gas friction coefficient fg to calculate the gas shear stress τg according to the relation:

For very low liquid loads flows, e.g. flows where USL < 0.01 m/s, the liquid phase at the bottom of the pipeline may become very thin and narrow. Here “USL” is superficial velocity of the liquid phase, i.e. the volumetric liquid flow divided by the cross-sectional area of the pipeline segment. Since this liquid phase is the source of the droplet field, the droplets need to travel laterally to deposit on the side walls, and if the lateral distance is very large, the droplets may not be able to reach the side walls. This effect is in the example embodiment of the invention taken into account by applying a correction which says that the width of the droplet field cannot be greater than a constant multiplied with the width of the bottom film. The best results were obtained by selecting a value of 7.0 for this constant. If the width of the droplet field is smaller than the pipe diameter, we assume that the part of the wall outside this range is dry. This effect is however typically unimportant, unless the liquid rate is very small (USL ≤ 0.01 m/s).

In a third aspect, the invention relates to a computer implemented method for designing a pipeline-based fluid transport system for transport of multiphase fluids, where the method comprises:

applying the method according to the first or the second aspect of the invention to numerically predict the pressure-drop when loading an intended flow rate volume of the multiphase fluid through a set of pipeline segments having different inner diameters and inclination angles towards the earth gravitational field, and

applying the determined pressure-drops to optimise the pipeline diameter and the pipeline trajectory/path in the terrain/seabed with aim of minimising material use in the pipeline-based fluid transport system and energy consumption for operating the pipeline-based fluid transport system.

The method according to the first aspect of the invention has the advantage of being significantly more accurate than known prior art one-dimensional computer implemented methods and also significantly less demanding in computing power for numerically predicting the pressure-loss in a pipeline segment. This enables a significantly more extensive use of the method to optimizing and trouble-shooting projected pipeline-based fluid transport systems than what is possible with known prior art models. For example, the under-prediction in pressure-drops shown in figure 1, would if applied to build a pipeline based on those predictions, result in a pipeline-based transport system lacking pumping power of around 1 MW for a 20 km long pipeline. There may be truly vast savings in costs, material consumption (pipeline steel, number of pressure boosting pumps, valves etc.), and energy consumption during operation of the pipeline-based fluid transport system by being able to predict relatively accurately the fluid flow behaviour in the pipeline system to identify a good balance need for pipeline dimensions (diameter etc.) and pipeline length (compromise between steepness and detour lengths in rugged terrain, etc.).

List of figures

Figure 1 is a diagram showing measured pressure drops in a pipeline segment for various flows as compared to predicted pressure-drop in a similar pipeline segment made by a prior art computational fluid dynamic model available under the trademark LedaFlow 2.4. The experimental measurements are presented in Kjølaas et al. [2].

Figure 2 is a drawing schematically illustrating the assumed droplet induced formation of a liquid film on the gas/inner wall interface.

Figure 3 is a diagram showing the ratio between surface roughness at the gas/inner wall interface and the droplet size determined by Hinze droplet model [10] plotted against capillary number. The dotted line marked “Equation (13)” corresponds to eqn. (7) herein.

Figures 4a) and 4b) display diagrams showing measured and predicted pressure drops for different two-phase flows at different gas flow velocities. The markers represent measured values, the dashed lines represent the predictions made by a prior art computational fluid dynamic model, while the solid lines represent the predictions made by a similar computational fluid dynamic model adapted according to the present invention.

Figure 5 displays diagrams showing measured and predicted normalized pressure drops for different three-phase flows at different gas flow velocities and water cuts. The markers represent the measured values, the dashed lines represent the predictions made by a prior art computational fluid dynamic model, while the solid lines represent the predictions made by a similar computational fluid dynamic model adapted according to the present invention.

Figure 6 is a facsimile of figure 14 of US 2010/0059221 showing an example of a production system.

Verification of the invention

The invention is demonstrated by first providing an example embodiment of a simple mathematical model for a two-phase flow showing how the mathematical model may be adapted to incorporate the correction of the shear forces according to the invention at the gas/wall interface. And further by presenting a comparison of predictions of the pressure drop made by a commercial computational fluid dynamic model marketed under the trademark “LedaFlow 2.4” and by a new version of the same commercial computational fluid dynamic model having incorporated the correction of the shear forces according to the invention and which is just being released under the trademark “LedaFlow 2.5”. The predictions are compared to experimental measurements of the pressure drop.

Example embodiment of a mathematical model having incorporated the invention The example embodiment of the computational fluid dynamical model is a so-called "steady-state point model", which predicts the pressure drop and phase fractions given a set of input parameters, which are:

1. Pipe inner diameter D, wall surface roughness ks and inclination θ.

2. The volumetric flow rates for each continuous fluid phase.

3. The thermodynamic properties of each continuous fluid phase (comprises at least densities, viscosities and surface tensions of the fluid phases).

The example embodiment of computational fluid dynamic model is called a "point model" because it only provides predictions locally in a system. By "locally", we mean a region small enough such that we can assume that all the fluid properties are constant in that region. The example embodiment further assumes the following:

4. The flow consists of one liquid phase and one gas phase (two-phase flow), both of which are considered incompressible.

5. The flow is stratified with liquid droplets, which means that the liquid flows as a continuous liquid phase at the bottom of the pipe, and possibly also as droplets entrained in the gas above.

6. The interface between the continuous liquid phase at the bottom of the pipe and the gas is assumed to be flat.

7. The flow is steady and fully developed, meaning that the flow does not change in time or in space. Mathematically, this is equivalent to assuming that all temporal and spatial derivatives in the mass/momentum equations are zero.

8. The liquid droplets travel at the same velocity as the gas.

The mass balance equations for the example may be given as:

where Qg and Qi are the volumetric gas and liquid flow rates, A is the cross section area of the pipe, αg, α1 and αdroplet are the volumetric fractions of the gas, the liquid film at the bottom of the pipe, respectively, and the liquid droplets, ug is the gas velocity, and ui is the velocity of the continuous liquid phase at the bottom of the pipe.

The momentum equations for the gas zone (continuous gas phase and liquid droplets) and the continuous liquid phase for this model are:

Here, dp/dx is the pressure gradient, pg is the gas density, pi is the liquid density, Sg is the gas-wall perimeter, Si is the film-wall perimeter, Si is the width of the gas liquid interface, τg is the gas-wall shear stress, xi is the liquid film shear stress, and τi is the shear stress between the gas and liquid at the interface. From geometrical considerations, it can be shown that the parameters α1 and S1 are related as follows for a flat continuous liquid layer in a circular pipe:

Furthermore, the relationship between the interface width Si and the film perimeter Si is:

Finally, we have that the sum of the wall perimeters Sg and Si must equal the total pipe perimeter.

The shear stresses can be expressed as follows:

Here, fg, fi and fi are the friction factors for the gas, the liquid film and the interface, and these friction factors require closure laws.

The example model incorporates the inventive determination of the gas flow resistance by applying the Håland model [13] as closure law for determination of the shear stress at the gas/inner wall interface, i.e. to determine the gas friction coefficient,fg, from the effective surface roughness, as defined in relations (1) or preferably by relations (11) to (13) given above.

In most prior art models, the keff parameter takes on the value for a dry pipeline wall, which is part of the user input to the model. However, in the example model, the value of keff is determined as given in equations (1) and (3) to (7) given above to take into account for droplets depositing on the wall, creating a liquid film at the gas/inner wall interface which has a roughness that is significantly different from the effective surface roughness of a dry pipe wall.

For the continuous liquid phase, the example model applies the same closure law, which in this case may be given as:

Here Rei is the Reynolds number of the continuous liquid, which is usually expressed as:

where:

However, for the liquid-wall shear stress, kef.i, takes on the value provided as user input for the surface roughness of the pipeline wall.

There are many proposed closure laws for the interfacial friction factor in the literature which may be applied to close the model over the gas-liquid interface. The example model applies a correlation proposed by Andritsos & Hanratty [14]:

Here, hi is the centre line depth of the continuous liquid phase.

To calculate the droplet fraction adroplet, another closure law is needed. The example model applies a model suggested by Ishii & Mishima [15] to calculate the total droplet fraction adroplet:

where

and

The example model is now complete as an example embodiment of the invention according to the first aspect.

However, the example model may further comprise a model for determination of the droplet concentration and thus be an example embodiment of the invention according to the second aspect. The example embodiment, for the sake of simplicity, assumes that the droplet concentration drops linearly with the vertical distance above the gas-liquid interface, and that it reaches a value of zero at the top of the pipe and that the pipe is shaped as a rectangular duct, the local concentration of droplets may be given by:

where y is the vertical distance from the interface, and hg is the distance from the gas-liquid interface to the top of the pipe.

Eqn. (32) may then be applied together with eqn. (2) to determine whether the effective surface roughness, keff ; at the gas/inner wall interface is to be determined by relation (1) or preferably by relations (12) to (14) given above, or to be defined to be equal to a surface roughness of a dry inner wall of the pipeline.

Comparison of measured and predicted pressure drops

The invention according to the invention is further verified by a comparison of predicted pressure-drops for various two-phase flows in a pipeline segment by a prior art commercially available computational fluid dynamic model as compared to the predicted pressure-drops by a lust released upgraded version of a similar computational fluid dynamic model having incorporated the inventive approach for determining the shear-stresses at the gas/inner wall interface. The predictions are compared to measured pressure-drops for different.

The experimental measurements of the pressure-drop were conducted at the SINTEF Multiphase Flow Laboratory applying their “Large Scale Loop Facility” adapted with a 94-meter-long 8"pipe with an inclination angle 2.5°. The nominal pressure was 60 bara, yielding a gas density of 67 kg/m3. For the experiments, nitrogen was used as the gas phase and Exxsol D60 as oil phase. For the aqueous phase, we used regular tap water with NaOH for corrosion protection, with and without glycerol. The purpose of adding glycerol was to increase the viscosity of the aqueous phase, emulating MEG injection. In the experiments with glycerol, the volumetric concentration was in the range 70-74%. The glycerol experiments were conducted at temperatures 23°C and 45°C yielding viscosities of about 42 and 14 cP, respectively, while the experiments without glycerol were conducted at 30°C. When changing the temperature, we also adjusted the pressure such that the gas density was kept the same in all the experiments. Additional details regarding the experimental set-up and execution are presented in Kjølaas et al. [1, 2],

The commercially available computational fluid dynamic model is sold under the trademark “LedaFlow 2.4” by LedaFlow Technologies DA. The LedaFlow 2.4 model describes thus the multiphase fluid as consisting of continuous fluid phases being separated by large-scale interfaces and of dispersed fluid phases suspended in one or more of the continuous fluid phases. The LedaFlow 2.4 applies the typically prior art approach for determining the shear-stresses at the gas/inner wall interface by assuming the effective surface roughness being equal to the surface roughness of a dry pipeline wall.

The upgraded version of the computational fluid dynamic model is going to be marketed under the trademark “LedaFlow 2.5” and is similar to the LedaFlow 2.4 model except that the assumptions presented above for the effect on the shearstresses by the formation of a liquid film on the gas/inner wall are incorporated into the model. The LedaFlow 2.5 model applies the approach according to the second aspect of the invention utilising the critical droplet concentration criteria as given in eqn. (2) above to determine how much of the gas-wall perimeter is dry (with an effective surface roughness equal to keff,2), and how much of the gas-wall perimeter is covered by a liquid film (with an effective surface roughness of equal to keff,1). The LedaFlow 2.5 model is adapted to apply eqn. (7) to determine keff,1 when the multiphase flow is a two-phase flow, and to apply eqn. (8) to determine keff,1 when the multiphase flow is a three-phase flow comprising two immiscible liquids and a gas phase.

Two-phase flow

Figure 4a) shows two diagrams presenting measured and predicted pressure drops plotted against the superficial gas velocity USG. In the top graph, the superficial liquid velocity USL is 0.01 m/s, and in the bottom graph USL equals 0.1 m/s.

Figure 4b) shows two diagrams presenting measured and predicted pressure drops plotted against the superficial liquid velocity USL. In the top graph, the superficial gas velocity USG is 9 m/s, and in the bottom graph USG equals 12 m/s.

The predictions are made by the prior art computational fluid dynamic model LedaFlow 2.4 and by the computational fluid dynamic model (LeadFlow 2.5) adapted according to the present invention. The experimental measurements are presented in Kjølaas et al. [1, 2]. The markers represent measured values, the dashed lines represent the predictions made by LedaFlow 2.4, while the solid lines represent the predictions made by LedaFlow 2.5.

As seen from the figures, the prior art model (LedaFlow 2.4) systematically underpredicts the pressure drop at these conditions, while the model adapted according to the present invention (LedaFlow 2.5) is in much better agreement with the measurements.

Three-phase flow

Figure 5 is a diagram showing the normalized pressure drop plotted against the water cut for three different gas rates (USG=6, 9 and 12 m/s), and two different liquid rates (USL=0.01 and 0.1 m/s). Each graph contains data for one or more fluid systems (with/without glycerol), where the main distinguishing factor is the water viscosity. The markers in the graphs represent the measured values, while the lines are predictions. The dashed lines are predictions obtained with the prior art LedaFlow 2.4 model, while the solid lines are predictions by the model (LedaFlow 2.5) adapted according to the invention. Both models applied the Brinkman emulsion model [12].

Figure 5 shows that also for the three-phase flow, the model adapted according to the invention is significantly more accurate in prediction the pressure-drop than the prior art model which does not compensate for the liquid film effect on the gas flow resistance.

It is noted that the model adapted according to the invention obtains a maximum pressure drop at a water cut of 50 % while the experiments indicate a maximum pressure-drop at a water cut of 80 %. The predictions in the high water cut range are thus low compared to experiments for this fluid system. It is assumed that the main reason for this discrepancy is that an oil-continuous emulsion forms on the wall. However, such behaviour is difficult model and is outside the objective of the hydrodynamic modelling approach of the present invention.

Nevertheless, the verification presented above clearly shows that the above presented assumption of describing the effective surface roughness by a liquid droplet model at the gas/inner wall interface when a liquid film may be assumed formed, is a simple and computational effective correction of the shear-stresses at the gas/inner wall interface obtaining significantly improved pressure-drop predictions as compared to prior art models assuming the shear-stresses are due to the surface roughness of the (dry) pipeline wall.

References

1. Jan Kjølaas, “Large scale experiments on high-rate three-phase low liquid loading flows”, 18th International Conference on Multiphase Production Technology, Cannes 2019.

2. Jan Kjølaas et al., “Pressure drop measurements in low liquid loading threephase flows”, 17th International Conference on Multiphase Production Technology, Cannes 2017.

3. J. Laurinat, T. Hanratty and W. Jepson, “Film thickness distribution for gas– liquid annular flow in a horizontal pipe,” PCH Physicochem. Hydrodyn., vol.

6, p. 179–195, 1985.

4. M. Bonizzi and P. Andreussi, “Prediction of the liquid film distribution in stratifieddispersed gas–liquid flow,” Chemical Engineering Science, vol. 142, p.165–179, 2016.

5. G. Wallis, One-dimensional Two-phase Flow, New York: McGraw-Hill, 1969.

6. R. J. Belt and L. M. Portela, “Prediction of the interfacial shear-stress in vertical annular flow”, International Journal of Multiphase Flow, vol. 35, no. 7, p. 689– 697, 2009.

7. J. Asali and T. Hanratty, “Interfacial Drag and Film Height for Vertical Annular Flow,” AIChE Journal, vol. 31, no. 6, pp. 895-902, 1985.

8. J. J. v. Rossum, “Experimental Investigation of Horizontal Liquid Films,” Chem. Eng. Sci., 11,35 ()., vol. 11, no. 1, pp. 35-52, 1959.

9. D. Biberg, C. Lawrence, G. Staff and H. Holm, “Pressure drop in low liquid loading flows - the effect of a thin liquid film on the pipe wall,” in BHRG, Cannes, France, 2017.

10. J. O. Hinze, “Fundamentals of the hydrodynamic mechanism of splitting in dispersion process,” AIChE, vol. 289, no. 1, 1955.

11. K. Arai, M. Konno, Y. Matunaga and S. Saito, “Effect of Dispersed-phase viscosity on the maximum stable drop size for break up in turbulent flow,” Journal of chemical engineering of Japan, vol. 10, no. 4, 1977.

12. H. Brinkman, “The viscosity of concentrated suspensions and solutions,” J. Chem.Phys., vol. 20, no. 4, pp. 571-584, 1952.

13. Haaland, SE. "Simple and Explicit Formulas for the Friction Factor in Turbulent Flow". Journal of Fluids Engineering, vol 105 no.1, pp. 89–90, 1983

14. N. Andritsos and T. J. Hanratty, "Interfacial instabilities for horizontal gasliquid flows in pipelines," Int. J. Multiphase Flow, vol. 13, pp. 583-603, 1987.

15. M. Ishii and K. Mishima, "Two-fluid model and hydrodynamic constitutive relations," Nuclear Engineering and Design, vol. 82, pp. 107-126, 1984.

Claims

CLAIMS

1. A computer implemented method for determining pressure-drop in a multiphase flow in a non-vertical pipeline having an inner wall facing the multiphase flow, where the multiphase flow at least comprises a continuous gas phase forming a gas/inner wall interface and a continuous liquid phase, and wherein the method comprises:

applying a computational fluid dynamic model describing the multiphase pipeline flow, and

solving the computational fluid dynamic model to determine the pressuredrop in the multiphase flow,

characterized in that the computational fluid dynamic model is adapted, for the gas/inner wall interface, to:

applying an effective surface roughness, keff ; determined by the relation: keff = where K is a correlation coefficient and is a maximum stable droplet size as determined by a droplet model, and

2. A method according to claim 1, wherein the computational fluid dynamic model describes the multiphase pipeline flow as a set of fluid fields representing:

- a set of continuous fluid phases being separated by a large-scale interfaces, at least comprising a continuous liquid phase and a continuous gas phase being separated by a large-scale gas/liquid interface, and

- a set of dispersed phases, at least comprising a dispersed phase of liquid droplets entrained in the continuous gas phase,

and wherein

- each fluid field is described by Eulerian formulated volume and ensemble averaged turbulent transport equations which are solved independently for each fluid field by applying the volume of fluid method for each fluid field representing a continuous phase and the dispersed flow method for each fluid field representing a dispersed fluid phase.

3. A computer implemented method for determining pressure-drop in a multiphase flow in a non-vertical pipeline having an inner wall facing the multiphase flow, where multiphase flow at least comprises a continuous gas phase forming a gas/inner wall interface and a continuous liquid phase, and

wherein the method comprises:

applying a computational fluid dynamic model describing the multiphase pipeline flow, and

solving the computational fluid dynamic model to determine the pressuredrop in the multiphase flow,

characterized in that the computational fluid dynamic model is adapted, for the gas/inner wall interface, to:

for the first fraction; applying an effective surface roughness, kefi, determined by the relation: , where K is a correlation coefficient and a maximum stable droplet size as determined by a droplet model, and applying the effective surface roughness, kef.i, to determine a shearstress, τg ,1, and

for the second fraction; applying an effective surface roughness, keff,2, defined to be equal to a surface roughness of a dry inner wall of the pipeline, and applying the effective surface roughness, keff, 2, to determine a shear-stress, τg 2, and applying the shear-stresses, τg1 and τg2 as boundary conditions determining the gas flow resistance at the gas/inner wall interface.

4. A method according to claim 3, wherein the computational fluid dynamic model describes the multiphase flow as a set of fluid fields representing:

- a set of dispersed phases, at least comprising a dispersed phase of liquid droplets entrained in the continuous gas phase,

and wherein

- each continuous fluid field is described by Eulerian formulated volume and ensemble averaged turbulent transport equations which are solved for each continuous fluid field, and

- each of the dispersed phases is described by "slip relations", which are expressions that describe the velocity difference between the dispersed fields and the associated continuous field.

5. A method according to claim 4, wherein the determination of the first fraction of the gas/inner wall interface considered being covered by a liquid film and the second fraction of the gas/inner wall interface considered not being covered by a liquid film is obtained by:

determining a liquid droplet concentration profile, in the continuous gas phase as a function of vertical distance above the large-scale interface separating the continuous liquid phase and the continuous gas phase, determining a critical droplet concentration by the relation: = — min[max(250 · σgl, 4.0), 9.0], where σgl is the gas-liquid surface tension, and

applying the determined liquid droplet concentration profile, at the gas/inner wall interface and define the part of the gas/inner wall interface where the is satisfied to be the first fraction, and similarly define

the part of the gas/inner wall interface where the is satisfied to be the second fraction.

6. A method according to any preceding claim, wherein is given by the relation:

where σ* = σgl (1 k · Ca ), where σgl is the surface tension between the liquid and the gas, k is a correlation coefficient, and Ca is a capillary number defined to be:

where μd is the dynamic viscosity of the liquid droplet, pg is volumetric density of the continuous gas phase, and ε is the energy dissipation rate per unit mass given by the relation:

where Sg is the gas perimeter, τg is the gas/inner wall shear stress, ug is the gas velocity, Si is the interface perimeter, τί is the interfacial shear stress, Δu is the gasliquid slip velocity, and Dhg is the hydraulic diameter for the continuous gas phase.

7. A method according to any of claims 1 - 5, wherein the multiphase flow is a two-dimensional flow and the effective surface roughness keff for liquid film covered gas/inner wall interface is given by the correlation:

where dHinze is dHinze = σgl is the gas-liquid surface tension, pg

is the gas density, and ε is the energy dissipation rate per unit mass given by the relation:

where Sg is the gas perimeter, τg is the gas/inner wall shear

stress, ug is the gas velocity, Si is the interface perimeter, τί is the interfacial shear stress, Δu is the gas-liquid slip velocity, and Dhg is the hydraulic diameter for the continuous gas phase, and Ca is a capillary number defined to be: Ca =

where is the dynamic viscosity of the liquid droplet.

8. A method according to any of claims 1 - 5, wherein the multiphase flow is a three-dimensional flow comprising a first and a second liquid, and

the correlation coefficient, K , is:

where Ca is the capillary number:

where σeff is an effective surface tension of the mixed first and second liquid:

where α1is the fraction of pipeline inner wall covered by the first

liquid, a2 is the fraction of pipeline inner wall covered by the second liquid, σ1 is the surface tension of the first liquid, σ2 is the surface tension of the second liquid, and

where is the effective dynamic viscosity of the mixed liquid:

where φ is the droplet concentration of the first liquid

being in suspension in second liquid phase, or alternatively, the droplet concentration of second liquid being in suspension in the first liquid phase.

9. A method according to any preceding claim, wherein the non-vertical pipeline has an inclination angle in the range of from - 75° to 75°, where a positive inclination angle corresponds to the pipeline ascending in the flow direction, a negative inclination angle corresponds to the pipeline descending in the flow direction, and an inclination angle of 0° corresponds to the pipeline being horizontal relative to the earth gravitational field.

10. A computer implemented method for designing a pipeline-based fluid transport system for transport of multiphase fluids, where the method comprises: applying the method according to any of claims 1 - 9 to numerically predict the pressure-drop when loading an intended flow rate volume of the multiphase fluid through a set of pipeline segments having different inner diameters and inclination angles towards the earth gravitational field, and

11. A computer program, comprising processing instructions which causes a computer to perform the method according to any of claims 1 – 9 when the instructions are executed by a processing device in the computer.

12. A computer, comprising a processing device and a computer memory, the computer memory is storing a computer program as set forth in claim 11.

KRAV

1. Datamaskinimplementert fremgangsmåte for å bestemme trykkfallet i en flerfasestrømning i en ikke-vertikal rørledning som har en innervegg vendt mot flerfasestrømningen, hvor flerfasestrømningen omfatter minst en kontinuerlig gassfase som danner en gass/innervegg-grenseflate og en kontinuerlig væskefase, og hvor fremgangsmåten omfatter å:

anvende en numerisk strømningsberegningsmodell som beskriver flerfasestrømningen i rørledningen, og

løse den numeriske strømningsberegningsmodellen for å bestemme trykkfallet i flerfasestrømningen,

karakterisert ved at den numeriske strømningsberegningsmodellen er tilpasset til, for gass/innervegg-grenseflaten, å:

anvende en effektiv overflateruhet, keff ; bestemt ved relasjonen:

hvor K er en korrelasjonskoeffisient og er største

stabile dråpestørrelse bestemt av en dråpemodell, og

anvende den effektive overflateruheten, keff, til å bestemme en skjærspenning, τg, og deretter anvende skjærspenningen, τg, som en grensebetingelse som bestemmer gasstrømningsmotstanden ved gass/innervegg-grenseflaten.

2. Fremgangsmåte ifølge krav 1, hvor den numeriske strømningsberegningsmodellen beskriver flerfase-strømningen i rørledningen som et sett av fluidfelt som representerer:

et sett kontinuerlige faser separert av storskal agrenseflater (ENG: large-scale interfaces), som omfatter minst en kontinuerlig væskefase og en kontinuerlig gassfase separert av en storskala gass/væske-grenseflate, og

et sett dispergerte faser, som omfatter minst en dispergert fase av væskedråper revet med av den kontinuerlige gassfasen,

og hvor

hvert fluidfelt blir beskrevet av Eulerformulerte volum- og ensemblemidlede turbulente transportligninger som løses uavhengig for hvert fluidfelt ved å anvende fluidvolummetoden (ENG: volume of fluid method) for hvert fluidfelt som representerer en kontinuerlig fase og den dispergerte strømningsmetoden for hvert fluidfelt som representerer en dispergert fluidfase.

3. Datamaskinimplementert fremgangsmåte for å bestemme trykkfallet i en flerfasestrømning i en ikke-vertikal rørledning som har en innervegg vendt mot flerfasestrømningen, hvor flerfasestrømningen omfatter minst en kontinuerlig gassfase som danner en gass/innervegg-grenseflate og en kontinuerlig væskefase, og hvor fremgangsmåten omfatter å:

anvende en numerisk strømningsberegningsmodell som beskriver flerfasestrømningen i rørledningen, og

løse den numeriske strømningsberegningsmodellen for å bestemme trykkfallet i flerfasestrømningen,

karakterisert ved at den numeriske strømningsberegningsmodellen er tilpasset til, for gass/innervegg-grenseflaten, å:

bestemme en første fraksjon av gass/innervegg-grenseflaten antatt å være dekt av en væskefilm og en andre fraksjon av gass/innervegg-grenseflaten antatt å ikke være dekt av en væskefilm,

for den første fraksjonen, anvende en effektiv overflateruhet, keff, 1, bestemt ved relasjonen: hvor K er en korrelasjonskoeffisient og er største stabile dråpestørrelse bestemt av en dråpemodell, og anvende den effektive overflateruheten, keff.i, til å bestemme en skjærspenning, τg ,1, og

for den andre fraksjonen, anvende en effektiv overflateruhet, keff2, definert å være lik overflateruheten til en tørr innervegg i rørledningen, og anvende den effektive overflateruheten, keff2, til å bestemme en skjærspenning, τg 2, og anvende skjærspenningene, τg ,1 og τg 2 som grensebetingelser som bestemmer gasstrømningsmotstanden ved gass/innervegg-grenseflaten.

4. Fremgangsmåte ifølge krav 3, hvor den numeriske strømningsberegningsmodellen beskriver flerfase-strømningen som et sett av fluidfelt som representerer:

et sett kontinuerlige fluidfaser separert av storskalagrenseflater (ENG: largescale interfaces), som omfatter minst en kontinuerlig væskefase og en kontinuerlig gassfase separert av en storskala gass/væske-grenseflate, og

et sett dispergerte faser, som omfatter minst en dispergert fase av væskedråper revet med av den kontinuerlige gassfasen,

og hvor

hvert fluidfelt blir beskrevet av Eulerformulerte volum- og ensemblemidlede turbulente transportligninger som løses for hvert kontinuerlige fluidfelt, og hver av de dispergerte fasene beskrives av «slipprelasjoner», som er uttrykk som beskriver hastighetsforskjellen mellom de dispergerte feltene og tilhørende kontinuerlige felt.

5. Fremgangsmåte ifølge krav 4, hvor bestemmelsen av den første fraksjonen av gass/vegg-grenseflaten antatt å være dekt av en væskefilm og den andre fraksjonen av gass/innervegg-grenseflaten antatt å ikke være dekt av en væskefilm oppnås ved å:

bestemme en konsentrasjonsprofil for væskedråper, i den kontinuerlige gassfasen som en funksjon av vertikalavstanden over storskalagrenseflaten som separerer den kontinuerlige væskefasen og den kontinuerlige gassfasen, bestemme en kritisk dråpekonsentrasjon, , ved relasjonen:

hvor σgl er gass/væske-

overflatespenningen, og

anvende den bestemte konsentrasjonsprofilen for væskedråper,

ved gass/innervegg-grenseflaten og definere den delen av gass/innervegggrenseflaten hvor er tilfredsstilt til å være den første fraksjonen, og tilsvarende definere den delen av gass/innervegg-grenseflaten hvor

er tilfredsstilt til å være den andre fraksjonen.

6. Fremgangsmåte ifølge ethvert foregående krav, hvor bestemmes av relasjonen:

hvor σ* = σgl (1 k ■ Ca), hvor σgl er overflatespenningen mellom væsken og gassen, k er en korrelasjonskoeffisient, og Ca er et kapillærtall definert å være:

hvor μα er skjærviskositet til væskedråpen, pg er volumetrisk tetthet til den kontinuerlige gassfasen, og ε er energidissipasjonshastigheten per masseenhet gitt ved relasjonen:

hvor Sg er gassperimeteren, τg er gass/innervegg skjærspenningen, ug er gasshastigheten, Si er grenseflateperimeteren, τί er grenseflateskjærspenningen, Δu er gass-væske-slipphastigheten, og Dhg er hydraulisk diameter for den kontinuerlige gassfasen.

7. Fremgangsmåte ifølge ethvert av krav 1 til 5, hvor flerfasestrømningen er en todimensjonal strømning og den effektive overflateruheten keff for væskefilmdekt gass/innervegg-grenseflaten er gitt ved korrelasjonen:

hvor dHinze er dHinze = σgl er gass-væske overflatespenningen,

pg er gasstettheten, og ε er energi dissipasjonshastigheten per masseenhet gitt ved relasjonen:

er gassperimeteren, τg er gass/innerveggskjærspenningen, ug er gasshastigheten, Si er grenseflateperimeteren, τί er grenseflateskjærspenningen, Δu er gass-væske-slipphastigheten, og Dhg er hydraulisk diameter for den kontinuerlige gassfasen, og Ca er kapillærtallet definert

til å være: hvor pd er skjærviskositet til væskedråpen.

8. Fremgangsmåte ifølge ethvert av krav 1 til 5, hvor flerfasestrømningen er en tredimensjonal strømning omfattende en første og en andre væske, og korrelasjonskoeffisienten, K, er:

hvor Ca er kapillærtallet:

hvor er en effektiv overflatespenning til den blandede første og andre væske:

hvor α1 er fraksjonen av den indre rørveggen dekt av den første væsken, α2 er fraksjonen av den indre rørveggen dekt av den andre væsken, σ1 er overflatespenningen til den første væsken, σ2 er overflatespenningen til den andre væsken, og hvor er den effektive skjærviskositeten til den blandede væsken:

hvor φ er dråpekonsentrasjonen til den første væsken i suspensjon i den andre væskefasen, eller alternativt, dråpekonsentrasjonen til den andre væsken i suspensjon i den første væskefasen.

9. Fremgangsmåte ifølge ethvert foregående krav, hvor den ikke-vertikale rørledningen har en skråstillingsvinkel i området fra - 75° til 75°, hvor en positiv skråstillingsvinkel samsvarer med en rørledning som stiger i strømningsretningen, en negativ skråstillingsvinkel samsvarer med en rørledning som synker i strømningsretningen og hvor en skråstillingsvinkel på 0° samsvarer med rørledning som er horisontal relativ til jordas tyngdekraftfelt.

10. Datamaskinimplementert fremgangsmåte for å utforme et rørledningsbasert transportsystem for å transport av flerfasestrømninger, hvor fremgangsmåten omfatter:

anvende en fremgangsmåte ifølge ethvert av krav 1 til 9 for å numerisk predikere trykktapet når man laster en tilsiktet volumhastighet av flerfasestrømningen gjennom et sett av rørledningssegmenter med forskjellige indre diametere og skråstillingsvinkler i forhold til jordens tyngdekraftfelt, og anvende de bestemte trykkfallene til å optimalisere rørledningsdiameter og rørledningstrase i terrenget/sjøbunnen med det formål å minimalisere materialforbruk i det rørledningsbaserte transportsystemet og energiforbruk under drift av det rørledningsbaserte transportsystemet.

11. Datamaskinprogram, omfattende prosesseringsinstruksjoner som får en datamaskin til å utføre fremgangsmåten ifølge et av kravene 1 til 9 når instruksjonene utføres av en prosesseringsenhet i datamaskinen.

12. Datamaskin, omfattende en prosesseringsinnretning og datamaskinminne hvor datamaskinminnet har lagret et datamaskinprogram ifølge krav 11.

Drawings

Priority

No priority

Patent classes

IPC classes

E21B 49/08

G06F 17/50

CPC classes

E21B 49/08

G06F 17/5009

G06F 2217/34

G06F 2217/16

Citations

D. Biberg, C. Lawrence, G. Staff og H. Holm, Pressure drop in low liquid loading flows – the effect of a thin liquid film on the pipe wall, BHR Group, 18th International Conference on Multiphase Production Technology, 7-9 June 2017, Cannes, France, OnePetro, Document ID: BHR-2017-471, https://www.onepetro.org/conference-paper/BHR-2017-471 ()

J. Kjølaas, S.T. Johansen, Y. Ladam, R. Belt, T.J. Danielson og M. Stinessen, Modeling of the droplet field in near-horizontal low liquid loading flows, BHR Group, 15th International Conference on Multiphase Production Technology, 15-17 June 2011, Cannes, France, OnePetro, Document ID: BHR-2011-I1, https://www.onepetro.org/conference-paper/BHR-2011-I1 ()

WO 2014082916 A1 (A1)

WO 2018191736 A1 (A1)

WO 2015110599 A1 (A1)

J. Kjølaas, T.E. Unander, M. Wolden, A. Shmueli og H. Holm, Pressure drop measurements in low liquid loading three-phase flows, BHR Group, 18th International Conference on Multiphase Production Technology, 7-9 June 2017, Cannes, France, OnePetro, Document ID: BHR-2017-455, https://www.onepetro.org/conference-paper/BHR-2017-455 ()

Applicants, owners and inventors

Owners

LEDAFLOW TECHNOLOGIES DA

c/o Kongsberg Oil and Gas Technolog, Drengsrudbekken 12

1383 ASKER

NO (ASKER Municipality, Akershus county)

Org. number: 912158152

Inventors

Jørn Kjølaas

Hårstadringen 17

7092 TILLER

NO (TRONDHEIM Municipality, Trøndelag county)

Agent

ACAPO ONSAGERS AS

Edvard Griegs vei 1

5059 BERGEN

NO (BERGEN Municipality, Vestland county)

Org. number: 932214636

Reference: P26649NO00

Status history

Payment history

Oppositions and appeals

No oppositions or appeals

Administrative limitation

No administrative limitation

Liens and licences

No liens or licenses registered

Gazette entries

No Gazette entries found

Only events from 2024.03.04 and onwards are shown.