Method for determining a mixture density measurement value, method for determining a flow rate measurement value based on a differential pressure measurement value and the mixture density measurement value, and a measuring point for this purpose.

Abstract

The method according to the invention serves to determine a mixture density measurement value ρ. m a gas-laden liquid using a Coriolis mass flow meter with a measuring tube for guiding the gas-laden liquid. The method comprises: determining (210) a mixture density start value based on at least one natural frequency of a bending vibration mode of the measuring tube; determining (220) a mass flow rate; determining (230) a gas volume fraction start value α s based on the initial mixture density value ρ s and a liquid density value ρ l ; Determining (240) a gas volume fraction final value α f based on the initial gas volume fraction value α s , and the mass flow rate; and (250) determining the mixture density value based on the liquid density value and the gas volume fraction final value α f .

Description

[0001] The invention relates to a method for determining a mixture density measurement value, a method for determining a flow rate measurement value based on a differential pressure measurement value and the mixture density measurement value, and a measuring point for this purpose.

[0002] Flow measurements based on a differential pressure reading require a density measurement as an auxiliary variable to determine a volumetric or mass flow rate. A Coriolis mass flow meter, connected in series with the differential pressure meter, is frequently used for density measurement. This is straightforward for homogeneous liquids; however, if the liquid contains free gas bubbles, the determination of a mixture density cannot be directly derived from the natural frequency of a bending mode of the Coriolis flow meter. Instead, an initial mixture density value is first determined based on the natural frequency, and a gas volume fraction is then calculated based on this initial value. The mixture density is subsequently determined based on the gas volume fraction and a value for the density of the liquid.However, with free gas bubbles, decoupling can occur between the gas bubbles and the liquid, rendering the underlying density model, which essentially assumes a homogeneous medium, unsuitable for determining the mixture density. Such a mixture density measurement is therefore not readily suitable as a basis for characterizing the medium. It is thus known to determine a gas volume fraction based on such a mixture density measurement, which depends on a relative deviation of the mixture density from the liquid density. The effect of the decoupling between free gas bubbles and the medium must also be compensated for. A mixture density measurement can then be calculated based on the corrected gas volume fraction. Approaches for determining a mixture density measurement are described, for example, in WO 2023 / 196960 A1 and WO 2021 / 169721 A1.Although satisfactory correction algorithms exist for the effect of decoupling at high flow rates, these leave much to be desired, especially at low flow rates; it is therefore the object of the present invention to remedy this.

[0003] The problem is solved according to the invention by the method for determining a mixture density measurement value according to claim 1, the method for determining a flow measurement value according to claim 10 and the measuring point according to claim 12.

[0004] The inventive method for determining a mixture density measurement value of a gas-laden liquid using a Coriolis mass flow meter with a measuring tube for guiding the gas-laden liquid comprises: Determining a mixture density starting value ρ sbased on at least one natural frequency of a bending vibration mode of the measuring tube; determining a mass flow rate measurement; determining a gas volume fraction starting value α s based on the initial mixture density value ρ s and a liquid density value ρ l Determining a gas volume fraction final value α f based on the initial gas volume fraction value α s , and the mass flow rate measurement; and determining the mixture density measurement based on the liquid density value and the gas volume fraction final value.

[0005] In a further development of the invention, the initial gas volume fraction value α s proportional to the relative deviation of the initial mixture density value ρ s from the liquid density value ρ l , also: αs=ρl−ρsρl

[0006] In a further development of the invention, the gas volume fraction final value α f proportional to the initial gas volume fraction value α sand inversely proportional to a flow-dependent correction term, which is a monotonically increasing function of the mass flow measurement, and in particular assumes a value less than 1.1 when the flow rate is zero.

[0007] In a further development of the invention, the correction term has a sigmoidal function or a logistic function, the argument of which is a linear function of the mass flow rate measurement dm / dt.

[0008] In a further development of the invention, the correction term has a sigmoidal function or a logistic function, the argument of which is a linear function of the reciprocal of the liquid density value ρ. l is.

[0009] In a further development of the invention, the flow-dependent correction term k(dm / dt) has the following form: k(dmdt)=1+(kmax−1)⋅S(dmdtρl), where k max is a monotonically decreasing function of the initial gas volume fraction value, where k max> 1, where k max < 2, and where k(dm / dt) approaches k for maximum mass flow measurements max strives.

[0010] In a further development of the invention, the flow-dependent correction term k max a linear function of a power of the initial value of the gas volume fraction α s is, which in particular has the following form: kmax=a⋅αs3+c where e is a negative exponent with a magnitude less than 1, and where a and c are constants.

[0011] In a further development of the invention, the mixture density measurement value is determined as the product of the gas volume fraction final value and the liquid density value.

[0012] In a further development of the invention, the liquid density value is determined in a state of the liquid in which the liquid is free of free gas bubbles.

[0013] A liquid is considered to be free of free gas bubbles in the sense of the invention, in particular, if a fluctuation of the initial mixture density value remains below a threshold value, above which an inhomogeneous medium is to be assumed.

[0014] In a further development of the invention, the liquid, when it is free of free gas bubbles, has a viscosity of no more than 40 cP, in particular no more than 25 cP.

[0015] The inventive method for determining a flow rate measurement value by means of a differential pressure measuring point comprising a differential pressure sensor and a Venturi tube as an induced pressure transmitter for the differential pressure sensor, wherein the Venturi tube is through which a gas-laden liquid flows, wherein the Venturi tube is connected in series with a Coriolis mass flow meter, wherein the method comprises: determining a mixture density measurement value with the Coriolis mass flow meter using a method according to one of the preceding claims; acquiring a differential pressure measurement value by means of the differential pressure sensor, wherein; and determining a flow rate measurement value for the flow through the Venturi tube based on the differential pressure measurement value and the mixture density measurement value.

[0016] The measuring point according to the invention comprises: a Coriolis mass flow meter; a differential pressure measuring device with a Venturi measuring tube as an differential pressure sensor which is connected in series with the Coriolis mass flow meter; and a calculation arrangement for carrying out the method according to the invention for determining a flow measurement value.

[0017] In a further development of the invention, the Coriolis mass flow meter has an inlet opening and an outlet opening for the gas-laden liquid, wherein a flow axis is defined which passes through a center of a cross-section of the inlet opening and a center of a cross-section of the outlet opening, wherein the flow axis deviates from the vertical by no more than 30°, for example no more than 15° and in particular no more than 5°.

[0018] The invention will now be explained in more detail with reference to the embodiment shown in the drawings. It shows: Fig. 1: a schematic representation of an embodiment of a measuring point according to the invention; Fig. 2: a flowchart of an embodiment of the present invention; Fig. 3a: a diagram showing flow-dependent density measurement errors based on a mixture density starting value without further corrections; Fig. 3b: a diagram showing flow-dependent density measurement errors based on a mixture density start value after determining a mixture density measurement value according to Hemp; and Fig. 3c: a diagram showing flow-dependent density measurement errors based on a mixture density starting value with a determination of the mixture density measurement value according to the invention.

[0019] The in Fig. Figure 1 illustrates an embodiment of a measuring point 100 according to the invention, comprising a pipeline 110 in which a Coriolis mass flow meter 120 and a differential pressure flow meter 130 are installed in a vertical orientation in a hydraulic series connection. The Coriolis mass flow meter 120 has two parallel measuring tubes 122, 124 for guiding a medium flowing in the pipeline. Coriolis mass flow meters with only one measuring tube or four measuring tubes are equally suitable for carrying out the invention. The Coriolis mass flow meter 120 further comprises a measuring and operating circuit 142 for exciting and evaluating measuring tube oscillations and for determining at least one initial mixture density value and a mass flow measurement value.The differential pressure flow meter 130 has a Venturi tube 132 as a differential pressure sensor, which is connected in series with the measuring tubes 122, 124 of the Coriolis mass flow meter, and a measuring and operating circuit 143, which is configured to determine at least one flow-dependent differential pressure measurement value. The measuring point 100 further comprises a calculation arrangement 140, which in the simplest case includes the measuring and operating circuits 142, 143 of the two flow meters connected via a communication path, in particular a data bus. Optionally, the calculation arrangement 140 further comprises an additional flow computer 144 which is connected to the two measuring and operating circuits 142, 143 via a communication path, in particular a data bus.If no flow computer is available, it is advantageous for the calculation of the final mixture density value to be performed by the measuring and operating circuit 142 of the Coriolis mass flow meter 120. This final mixture density value is then transmitted to the measuring and operating circuit 143 of the differential pressure flow meter 130, which then calculates a flow rate based on a differential pressure reading and the final mixture density value. If a flow computer 144 is available, process steps that perform calculations based on primary measured variables, such as the initial mixture density value, the mass flow rate reading from the Coriolis mass flow meter, and the differential pressure reading, can be carried out by the flow computer 144.

[0020] An embodiment of the method according to the invention will now be described using the following examples: Fig. 2 explained.

[0021] The method 200 according to the invention comprises determining 210 a mixture density start value based on at least one natural frequency of a bending vibration mode of the at least one measuring tube. For this purpose, a bending vibration mode of the at least one measuring tube is excited, and the associated natural frequency is determined, for example, by amplitude maximization or phase control. Details of this are known to those skilled in the art. The method 200 further comprises determining 220 a mass flow rate measurement value dm / dt, which is carried out in a known manner based on a phase or time difference between the vibration signals of an inlet-side and an outlet-side vibration sensor.

[0022] To determine a correct mixture density measurement value, a gas volume fraction starting value α is first determined. s based on the initial mixture density value ρ s and a liquid density value ρ l, where a deviation of the initial mixture density value from the liquid density value ρ l as being caused exclusively by a gas volume fraction of the mixture with a density of zero. For pressure ranges up to a few MPa, this assumption leads to satisfactory results. The initial gas volume fraction α s In this case, this corresponds to a relative density deviation of the initial mixture density value ρ. s from the liquid density value ρ l , also: αs=ρl−ρsρl

[0023] The method according to the invention 200 further comprises determining 240 a gas volume fraction final value α f based on the initial gas volume fraction value α s and the mass flow rate measurement value dm / dt. The gas volume fraction final value α f The process is initially proportional to the starting value α of the gas volume fraction. s applied, and then divided by a correction factor k, which takes the following form: k(dmdt,αs)=1+(kmax(αs)−1)⋅S(dmdt), where k max (α s ) a linear function of a power of the initial value of the gas volume fraction α s is, which in particular has the following form: kmax=a⋅αse+c where k max > 1, where k max < 2.

[0024] The flow rate-dependent factor S(dm / dt) has, in particular, a sigmoidal or logistic function whose argument is a linear function of the mass flow rate measurement. Furthermore, the argument of the sigmoidal or logistic function is preferably a linear function of the reciprocal of the liquid density value ρ. l This means that the resulting argument is approximately proportional to the volumetric flow rate or flow velocity. In principle, any type of sigmoidal function is suitable, for example: S(x)=11+e−x+C, S(x)=tanh(x)+C, or S(x)=x1+x2+C, where C is a constant and where x is a linear function of a quotient of mass flow rate measurement dm / dt and liquid density value ρ l is, so x=k⋅dm / dtρl+z, where k and z are constants.

[0025] One embodiment of the sigmoidal function has the following form: S(dm / dtρl)=11+e−x, where 3 < k < 6 and -6 < z < -3.

[0026] In a final step, the mixture density measurement value ρ is determined. m by multiplying the final value of the gas volume fraction α f with the liquid density value ρ m .

[0027] Optionally, the method 200 can be further developed for flow measurement based on a flow-dependent differential pressure measurement value, in particular on a Venturi tube as a differential pressure transmitter, which is arranged in series with the Coriolis mass flow meter, with which the mixture density measurement value ρm is determined. In this case, the acquisition of a differential pressure measurement Δp (260) and the calculation of a flow rate measurement (270) as a function of the differential pressure measurement Δp and the mixture density measurement ρ follow. m .

[0028] Fig. 3a, Fig. 3b and Fig. 3c show relative density measurement errors (ρ mix - ρ t ) / ρ t as a function of the mass flow rate for different gas loadings, where the gas loadings cover the same value ranges in all three diagrams. Here, ρ denotes t the actual or true mixture density and ρ mix the mixture density measurement value determined according to the respective method.

[0029] In Fig. 3a is ρ mix = ρ s The relative density measurement error of the initial mixture density value is thus shown. It is obvious that the initial mixture density value ρ sis unsuitable as a mixture density measurement value over the entire flow range, since the mixture density ρ s is systematically underestimated with increasing mass flow.

[0030] In Fig. 3b is ρ mix = ρ S Figure 1 shows the relative density measurement error of a mixture density measurement where the correction term k has no flow dependence, i.e., when the sigmoid function S is permanently set to 1. In this case, the correction term has the maximum value k = max, independent of the flow rate. This corresponds to Hemp's approach. However, this approach is also unsuitable for use as a mixture density measurement over the entire flow rate range.

[0031] In Fig. 3c is ρ mix = ρ m . Accordingly, the relative density measurement error of the mixture density measurement value determined according to the invention is ρ mshown. Accordingly, the latter is suitable to serve as a mixture density measurement value over the entire flow range. QUOTES INCLUDED IN THE DESCRIPTION

[0000] This list of documents cited by the applicant was automatically generated and is included solely for the reader's convenience. The list is not part of the German patent or utility model application. The DPMA accepts no liability for any errors or omissions. Cited patent literature

[0000] WO 2023 / 196960 A1

[0002] WO 2021 / 169721 A1

[0002]

Claims

[1] Method for determining a mixture density measurement ρ m a gas-laden liquid by means of a Coriolis mass flow meter with a measuring tube for guiding the gas-laden liquid, comprising: Determine (210) a mixture density start value based on at least one natural frequency of a bending vibration mode of the measuring tube; Determining (220) a mass flow rate measurement; Determining (230) a gas volume fraction starting value α s based on the initial mixture density value ρ s and a liquid density value ρ l Determining (240) a gas volume fraction final value α f based on the initial gas volume fraction value α s , and the mass flow rate measurement; and (250) Determining the mixture density measurement value based on the liquid density value and the gas volume fraction final value α f . [2] Method according to claim 1, wherein the initial gas volume fraction value α s is proportional to the relative deviation of the initial mixture density value ρ s from the liquid density value ρ l , also: αs=ρl−ρsρl [3] Method according to claim 1 or 2, wherein the final gas volume fraction α f proportional to the initial gas volume fraction value α s and inversely proportional to a flow-dependent correction term, which is a monotonically increasing function of the mass flow measurement, and in particular assumes a value less than 1.1 when the flow rate is zero. [4] Method according to claim 3, wherein the correction term has a sigmoid function or a logistic function whose argument is a linear function of the mass flow measurement. [5] Method according to claim 3 or 4, wherein the correction term has a sigmoid function or a logistic function whose argument is a linear function of the reciprocal of the liquid density value. [6] Method according to any one of claims 3 to 5, wherein the flow-dependent correction term k(dm / dt) has the following form: k(dmdt)=1+(kmax−1)⋅S(dmdtρl), where k max a monotonically decreasing function of the initial value of the gas volume fraction α s is, where k max > 1, where k max < 2, and where k(dm / dt) for maximum mass flow measurements versus k max strives. [7] Method according to claim 6, wherein the flow-dependent correction term k max a linear function of a power of the initial value of the gas volume fraction, which in particular has the following form: kmax=a⋅αs3+c where e is a negative exponent with a magnitude less than 1, and where a and c are constants. [8] Method according to one of the preceding claims, wherein the mixture density measurement value is determined as the product of the gas volume fraction final value and the liquid density value. [9] Method according to any of the preceding claims, wherein the liquid density value is determined in a state of the liquid in which the liquid is free of free gas bubbles. [10] Method according to any of the preceding claims, wherein the liquid, when free of free gas bubbles, has a viscosity of not more than 40 cP, in particular not more than 25 cP. [11] Method for determining a flow measurement value using a differential pressure measuring point comprising a differential pressure sensor and a Venturi tube as an differential pressure transmitter for the differential pressure sensor, wherein the Venturi tube is through which a gas-laden liquid flows, wherein the Venturi tube is connected in series with a Coriolis mass flow meter, the method comprising: Determining a mixture density measurement value with the Coriolis mass flow meter using a method according to one of the preceding claims; Acquiring (260) a differential pressure measurement Δp using the differential pressure sensor; and Determine (270) a flow rate measurement for the flow through the Venturi tube based on the differential pressure measurement Δp and the mixture density measurement. [12] Measuring point, comprising a Coriolis mass flow meter; a differential pressure gauge with a Venturi measuring tube, as a differential pressure sensor which is connected in series with the Coriolis mass flow meter; and a calculation arrangement for carrying out the method according to claim 10. [13] Measuring point according to claim 11, wherein the Coriolis mass flow meter has an inlet opening and an outlet opening for the gas-laden liquid, wherein a flow axis is defined which passes through a center of a cross-section of the inlet opening and a center of a cross-section of the outlet opening, wherein the flow axis deviates from the vertical by no more than 30°, for example no more than 15° and in particular no more than 5°.

Images