The present invention relates to a method for ascertaining flow of a fluid through a circularly cylindrical measuring tube having a straight, measuring tube, longitudinal axis and an inner diameter DI.
Ultrasonic, flow measuring devices are applied widely in process and automation technology. They permit easy determination of volume flow and/or mass flow in a pipeline.
Known ultrasonic, flow measuring devices frequently work according to the travel-time difference principle. In the travel-time difference principle, the different travel times of ultrasonic waves, especially ultrasonic pulses, i.e. so-called bursts, are evaluated as a function of the direction the waves travel in the flowing liquid. To this end, ultrasonic pulses are sent at a certain angle to the tube axis both with, as well as also counter to, the flow. From the travel-time difference, the flow velocity, and therewith, in the case of known diameter of the pipeline section, the volume flow, can be determined, for example, according to the formula, Q=K*((t1−t2)/(t1*t2)), wherein K is a function of the length of the signal path, the ratio between radial and axial sensor separations, the velocity distribution, respectively the flow profile in the measuring tube, and the cross sectional area, and t1, respectively t2, are the travel times of the signals upstream-, respectively downstream.
In the case of the Doppler principle, ultrasonic waves of a certain frequency are coupled into the liquid and the ultrasonic waves reflected by the liquid are evaluated. From the frequency shift between the coupled and reflected waves, the flow velocity of the liquid can likewise be determined. Reflections in the liquid occur, when small air bubbles or impurities are present in it, so that this principle is applied mainly in the case of contaminated liquids.
The ultrasonic waves are produced, respectively received, with the assistance of so-called ultrasonic transducers. To this end, ultrasonic transducers are placed securely in the tube wall of the relevant pipeline section. There are also clamp on, ultrasonic, flow measuring systems. In such case, the ultrasonic transducers are pressed externally on the wall of the measuring tube. A great advantage of clamp on, ultrasonic, flow measuring systems is that they do not contact the measured medium and can be placed on an already existing pipeline.
A further ultrasonic, flow measuring device working according to the travel-time difference principle is disclosed in U.S. Pat. No. 5,052,230. In such case, the travel time is ascertained by means of short ultrasonic pulses, so-called bursts.
The ultrasonic transducers are normally composed of an electromechanical transducer element, e.g. a piezoelectric element, and a coupling layer. The ultrasonic waves are produced in the electromechanical transducer element and led via the coupling layer to the pipe wall and from there into the liquid in the case of clamp-on-systems, and, in the case of inline systems, via the coupling layer into the measured medium. In such case, the coupling layer is sometimes called a membrane, or diaphragm.
Between the piezoelectric element and the coupling layer, another coupling layer can be arranged, a so called adapting, or matching, layer. The adapting, or matching, layer performs, in such case, the function of transmitting the ultrasonic signal and simultaneously reducing reflection at interfaces between two materials caused by different acoustic impedances.
Both in the case of clamp-on-systems, as well as also in the case of inline systems, the ultrasonic transducers are arranged on the measuring tube in a shared plane, either on oppositely lying sides of the measuring tube, in which case the acoustic signal, projected onto a tube cross section, passes once along a secant through the measuring tube, or on the same side of the measuring tube, in which case the acoustic signal is reflected on the oppositely lying side of the measuring tube, whereby the acoustic signal traverses the measuring tube twice along the secant projected on the cross section through the measuring tube. U.S. Pat. No. 4,103,551 and U.S. Pat. No. 4,610,167 show ultrasonic, flow measuring devices with reflections on reflection surfaces provided therefor in the measuring tube. Also known are multipath systems, which have a number of ultrasonic transducer pairs, which, in each case, form a signal path, along which the acoustic signals pass through the measuring tube. The respective signal paths and the associated ultrasonic transducers lie, in such case, in mutually parallel planes parallel to the measuring tube axis. U.S. Pat. No. 4,024,760 and U.S. Pat. No. 7,706,986 show such multipath systems by way of example. An advantage of multipath systems is that they can measure the profile of the flow of the measured medium in the measuring tube at a plurality of locations and thereby provide highly accurate, measured values for the flow. This is achieved based on, among other things, the fact that the individual travel times along the different signal paths are weighted differently. Disadvantageous in the case of multipath systems is, however, their manufacturing costs, since several ultrasonic transducers and, in given cases, a complex evaluating electronics need to be used.
There are different approaches for weighting the signal paths. The paper “Comparison of integration methods for multipath acoustic discharge measurements” by T. Tresch, T. Staubli and P. Gruber in the handout for 6th International Conference on Innovation in Hydraulic Efficiency Measurements, 30 Jul.-1 Aug. 2006 in Portland, Oreg., USA, compares established methods for weighting the travel times along different signal paths for calculating the flow
DE 10 2005 059 062 B4 and DE 10 2006 030 964 A1 disclose methods for correcting a first flow value of a gaseous fluid flowing through a measuring tube, wherein steam is a component. The concentration of the steam is determined or established by means of temperature and/or velocity of sound and then the concentrations of one or more components of the gaseous fluid are ascertained and the flow value corrected.
U.S. Pat. No. 5,835,884 A discloses determining the average flow velocity of a fluid. In such case, volume flow rate is measured in the laminar range (RN=2000) and in the turbulent range (RN=4000) and the average flow velocity for Reynolds numbers between 2000 and 4000 ascertained between the two values by a logarithmic interpolation method. An application of this method to hydrocarbon containing gas mixtures is not disclosed.
JP 56 140 214 A, U.S. Pat. No. 4,300,400, U.S. Pat. No. 5,546,813, EP 1 113 247 A1 and U.S. Pat. No. 4,331,025 A disclose methods for calculating flow velocity based on a function of Reynolds number Re and radius r. None of these documents is concerned, however, with the problem of measuring gas mixtures and the particular issues arising in such case.
The aforementioned documents are concerned exclusively with measuring flow of a fluid, however, not specially with a gas mixture, in the case of which not only the flow measurement—but, instead, also the composition is of interest and in the case of which individual, ascertained values of measured variables can be taken into consideration for determining the flow measurement and the composition for a determining of further physical variables and properties.
The present method begins, thus, with the object of providing a corresponding method, which overcomes the described problems.
An object of the invention is to provide a method for flow measurement by means of ultrasound, designed especially also for gas mixtures and delivering highly accurate measurement results.
The object is achieved by the subject matter of the independent claim 1. Further developments and embodiments of the invention are provided by the features of the respectively dependent claims.
According to the invention, a travel-time difference measurement of acoustic signals along a signal path is performed in a circularly cylindrical measuring tube having a straight, measuring tube, longitudinal axis and an inner diameter DI. This is accomplished preferably with a suitable ultrasonic, flow measuring device. The travel-time difference measurement of acoustic signals along a signal path between two ultrasonic transducers in the upstream- and downstream directions is known to those skilled in the art.
Serving both as transmitter as well as also receiver are usually ultrasonic transducers, especially electromechanical transducers, e.g. piezoelectric elements, which are suitable to send as well as also to receive the acoustic signal, especially an ultrasonic pulse or one or more ultrasonic waves. If ultrasonic transducers are applied as transmitters and receivers, the acoustic signal can pass along the first signal path back and forth, thus in two directions. Transmitter and receiver are, thus, exchangeable.
Referred to as the signal path, also called an acoustic path, is the path of the acoustic signal, thus e.g. the ultrasonic wave or the ultrasonic pulse, between the transmitter, which transmits the acoustic signal, and the receiver, which receives the acoustic signal. In an embodiment of the invention, the acoustic signal is, such as usual in the case of an inline system, radiated perpendicularly to the membrane. The receiver is then so emplaced in or on the measuring tube that the signal, in turn, strikes perpendicularly on its membrane.
If the signal path is composed of a plurality of straight subsections, thus, if, for example, the acoustic signal is reflected on one or more reflection surfaces, which are interfaces formed e.g. between fluid and measuring tube or a reflector arranged on or in the measuring tube, all straight subsections have the same separation from the measuring tube axis, especially the signal path and therewith all of its subsections j extend in a plane parallel to the measuring tube axis, which separation dj is especially unequal to about a fourth of the inner diameter DI (dj≠DI/4). In a further development of the invention, the signal path lies in a plane, in which the measuring tube axis lies. Projected on a cross section of the measuring tube, the inner diameter DI results, since the separation of all subsections of the signal path from the measuring tube, longitudinal axis is zero. The result of the travel-time difference measurement is an average flow velocity vL.
A further method step in an embodiment of the method of the invention is the ascertaining of the kinematic viscosity vkin of the fluid. The kinematic viscosity νkin is related to the dynamic viscosity μdyn in the following way: νkin=μdyn/ρ. Thus, if the dynamic viscosity μdyn is ascertained and the density ρ is known or itself ascertained, then the kinematic viscosity μkin is at hand.
There are many variants, by which the kinematic viscosity νkin of the gas mixture can be ascertained. Examples include using a table, a mathematical formula or linear interpolation between known values. The kinematic viscosity νkin of the fluid can, in such case, depend on different variables and can be correspondingly ascertained.
If the chemical composition of the gas mixture is known in terms of the individual material quantity fractions xi of its components i in the case of a multicomponent system, for example, via input provided by the user or by, in given cases also separately, ascertaining such, the kinematic viscosity vkin of the fluid is ascertained, for example, via the supplemental input of the temperature T of the fluid. In this regard, a temperature sensor can be provided.
Taking into consideration the material fractional amount xi of the individual components i of the gas mixture, it is to be understood that at least the velocity of sound and the temperature are predetermined or measured, since these variables enable calculation of the material fractional amounts. Thus, the material quantity fraction can be included directly in the calculation of the dynamic viscosity and therewith be taken into consideration.
Since the material quantity fraction can be calculated via the velocity of sound and the temperature of the medium, thus, by taking into consideration the temperature and the velocity of sound in the calculation, the kinematic viscosity of the material quantity fraction can be indirectly taken into consideration.
However, in the case of gas mixtures, especially in the case of biogases, the composition of the gas mixture can vary. In such case, the variable kinematic viscosity can be determined by a so-called, real time measurement. This means that, supplementally to flow, at least one changeable variable is measured repeatedly at a time interval. This is preferably the velocity of sound in the gas mixture, from which, then, in the case of constant temperature and constant pressure, an inference of a change in the material quantity fractions of the gas mixture and/or directly of the kinematic viscosity can be made. A preferred repetition interval lies between 5-500 msec (milliseconds), especially preferably, however, between 10-250 msec.
It is advantageous, when the kinematic viscosity νkin of the gas mixture is ascertained by measuring the temperature of the gas mixture and the velocity of sound c in the gas mixture, as well as from certain variables required for determining material-specific properties.
The dynamic viscosity of the gas mixture results advantageously by specifying and/or measuring
Exactly in the case of gas mixtures with time variable composition, for example, biogas, the temperature of the gas mixture often changes overall. These changes must be taken into consideration in determining the kinematic viscosity. In such case, a one-time measurement can be insufficient for these variables and a measurement repetition at a time interval can be advantageous. The measuring interval for temperature measurement lies, in such case, preferably at a maximum of 5 min, especially between 5 sec and 2 min.
Additionally, also an optional measuring of the pressure and the relative humidity can be repeated at the aforementioned time intervals.
In order to enable an exact measuring with small error, it is advantageous, when the hydrocarbon has a material quantity fraction xi of at least 0.1 with reference to the total mass of the gas mixture.
Alternatively, the kinematic viscosity νkin of the fluid is ascertained as a function of the velocity of sound c in the fluid, the temperature T of the fluid, the absolute pressure p of the fluid and the chemical composition of the fluid. Velocity of sound c in the fluid and temperature T of the fluid can, in such case, be ascertained in known manner by the ultrasonic, flow measuring device, or they can be separately ascertained. In the same way, also the density ρ of the fluid is ascertainable.
These are only some examples without any claim of completeness. It is not intended that other methods of ascertaining the kinematic viscosity νkin of the fluid should therewith be excluded.
Thus, other method steps can precede the method step of ascertaining the modified Reynolds number Remod. Examples include ascertaining the chemical composition of the fluid and/or ascertaining the material quantity fractions xi of the individual components i of the fluid, wherein these can also be predetermined by the user, and/or ascertaining the velocity of sound c in the fluid and/or ascertaining the temperature T of the fluid and/or ascertaining the absolute pressure p in the fluid, wherein then the kinematic viscosity νkin of the fluid is ascertained in suitable manner as a function of one or more of these parameters. In the case of gaseous fluids, ascertaining the absolute pressure p plays a greater role than in the case of liquid fluids, since most of these can be considered, for practical purposes, as incompressible.
If, according to a further development of the invention, the fluid is a gas, especially a biogas, with the components methane, water and carbon dioxide, which biogas also can have other components, such as e.g. nitrogen, oxygen, hydrogen, hydrogen sulfide and/or ammonia, then DE 10 2006 030 964 A1 teaches assuming the relative humidity of the fluid to be 100% or supplementally to provide a humidity measuring unit, in order to ascertain the concentration of water as a function of temperature T and the relative humidity RH and to take such into consideration in determining the concentrations of methane and carbon dioxide. This should likewise be included here.
In the next method step, a modified Reynolds number Remod is ascertained according to the formula Remod=(vL*DI)/νkin, wherein then a second flow velocity vA averaged over the cross sectional area of the measuring tube is ascertained by means of a known function vA=f(Remod) as a function of the modified Reynolds number Remod and, according to a further development of the invention, output by the device. In such case, the function vA=f(Remod) in the sense the present invention does not express a formula in the mathematical sense, but, instead, a proportionality between vA and f(Remod).
Taking into consideration the first average flow velocity vL, the formula for calculating the second average flow velocity vA becomes: VA=f(Remod)*VL
In a variant of the invention, the volume flow QV=vA*(π/4)*DI2 and/or the mass flow QM=QV*ρ with the density ρ of the fluid are/is calculated and then output by the device.
For ascertaining the function vA=f(Remod), there exist, analogously to the ascertaining of the kinematic viscosity νkin, likewise many options. One of these is to investigate the ratio vL/vA as a function of Reynolds number Re, respectively modified Reynolds number Remod, e.g. in a suitable calibration plant, experimentally in greater detail and to keep such in the form of a function f. In the case of constant Reynolds number, vL is proportional to vA: vA=f(Remod)*vL. The relationship vA/vL, versus Remod is generally true for all fluids. Therefore, it is not absolutely necessary to use in the calibration plant the same fluid as in the field.
Applied for performing the method is an ultrasonic, flow measuring device having a circularly cylindrical measuring tube having a straight, measuring tube, longitudinal axis and an inner diameter DI, two ultrasonic transducers for travel-time difference measurement of an acoustic signal along a signal path in the measuring tube and a suitable transmitter unit for evaluating the travel-time difference measurement and for performing the method of the invention, especially a so called inline, ultrasonic, flow measuring device having a measuring- or signal path, which is arranged centrally.
The invention is amenable to numerous forms of embodiment. One thereof will now be explained in greater detail based on the appended drawing, the figures of which show as follows:
The steam fraction is estimated or measured with a humidity measuring unit.
Then, via the measured velocity of sound c and the measured temperature T and, in given cases, the measured pressure p, the dynamic, or also the kinematic, viscosity of the biogas can be ascertained via corresponding algorithms. The formula Remod=(vL*DI)/νkin yields the modified Reynolds number.
From a known relationship vA/vL versus Re, then the flow velocity vA output by the flow measuring device can be corrected as a function of the Reynolds number.
The Reynolds number is obtained via the formula, Re=(vA*DI)/νkin, wherein vA is the flow velocity of the fluid through the measuring tube averaged over the total measuring tube cross section. vA is, thus, the surface integral. vL is, in contrast, the average flow velocity measured along the signal path and, correspondingly, the line integral along the signal path.
In the following based on an example of an algorithm, ascertaining of the dynamic viscosity will now be presented.
Based on this algorithm, one can recognize that the dynamic viscosity is calculable based on the temperature of the biogas and on the material fractional amount of methane in the biogas. In such case, the material quantity fraction is expressed as a molar fraction, respectively volume fraction, in % and temperature in ° C. lies in a range between 0-80° C., wherein the viscosity can be calculated with an accuracy of 0.5%—preferably 0.2%—at 1 bar, to the extent that no foreign gas influence is present.
The density in kg/m3 can be calculated via the following formula
with the predetermined or ascertained pressure being expressed in mbar and the measured temperature in degrees Kelvin.
In such case, K is calculated as follows:
wherein X is scaled between 0-1.
The kinematic viscosity can then be ascertained from the relationship: ν=η/ρ.
The Reynolds number Re exhibits the following dependence:
Re=V.D/ν or Re=ρ.V.D/η
Number | Date | Country | Kind |
---|---|---|---|
10 2011 080 365.3 | Aug 2011 | DE | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/EP2012/064370 | 7/23/2012 | WO | 00 | 1/29/2014 |