The invention relates to special acoustic measurements in a general level, in which low frequency sound, propagating in the fundamental mode of the flow tube is used. In a more particular, the invention relates to a flow measurement arrangement according to the preamble of an independent claim thereof. The invention relates also to an acoustic flow meter as indicated in the preamble of an independent claim thereof. The invention relates also to an acoustic flow determination method in a flow channel according to the preamble of an independent method claim. The invention relates also to a software program product according to the preamble of an independent claim thereof. The invention relates to a flow measurement system according to the preamble of an independent claim thereof.
It is known that long wavelength sound propagating in the so called “piston mode” in a flow channel of rigid walls is a good choice for acoustic flow metering, since the speed of sound is the sum of the speed of sound at rest plus or minus the mean flow velocity in the measuring volume chosen. [B. Robertson. “Effect of arbitrary temperature and flow profiles on the speed of sound in a pipe”, J. Acoust. Soc. Am., Vol. 62, No. 4, pp. 813-818, October 1977 and B. Robertson, “Flow and temperature profile independence of flow measurements using long acoustic waves”, Transactions of ASME, Vol. 106, pp. 18-20, March 1984].
One of the obstacles of such measurement is how to distinguish the downstream and upstream propagating sound waves from each other. Directional filtering, as described in patents [FI 94909, U.S. Pat. No. 5,770,806], is one possibility. Achieving good results requires both amplitude and phase matching of four microphones over wide temperature range, which is very difficult to achieve in practice.
It is known as such to use differential mode sound for determining the sound velocity from the resonance frequency of a resonator, from the disclosure in the Internet in http://www.acoustics.org/press/150th/Garrett.html.
In addition to solving problems relating to the above-mentioned known techniques, it is an objective of the invention to improve the measurement accuracy and/or to speed up the flow determination. This is achieved by the flow meter arrangement arranged for a flow determination according to the embodiments of the invention.
An acoustic flow meter arrangement, comprising a measuring section according to the invention is characterized in an independent claim thereof. A measuring section according to the invention is characterized in an independent claim thereof. A sound reflector/attenuator according to the invention is characterized in an independent claim thereof. An acoustic flow metering method according to the invention is characterized in an independent method claim. A software program product according to the invention is characterized in an independent claim thereof. An acoustic flow meter system according to the invention is characterized in an independent claim thereof. A measuring section (1) provided with a sound source is characterized in an independent claim thereof. An acoustic flow meter arrangement for differential sound measurement is characterized in an independent claim thereof.
Other embodiments of the invention are illustratively embodied and indicated in the dependent claims. The embodiments of the invention are combinable in suitable part.
In a flow arrangement according to an embodiment of the invention, the flow measurement is arranged to occur in a flow channel comprising a measuring section. According to an embodiment of the invention, one end of the measuring section is designed to act as a reflector, returning a quantifiable echo of the wave entering from the measuring section side. The downstream and upstream sounds are decomposed from the measured sound signals of only two sound sensors based on the echo model. Thereby it becomes possible to determine the difference of the transit times between the sound sensors directly without determining first the downstream and upstream transit times separately. This is unlike any other acoustic transit time measurement using the same measurement volume in both directions. The invention relates also to an acoustic flow meter. The arrangement embodied in the invention can be embodied as a discrete but functional ensemble of the parts of the arrangement. However, in suitable part, such parts can be assembled to form a device, a flow meter. Data acquisition means, processor, a database, memory and a transmission line for transferring measurement related data acquired by said data acquisition means can be at least partly assembled into the same cover for the device. The device can be a part of a system, where there are several flow meter arrangements in combination. According to an embodiment the invention, the system comprises communication means arranged to receive, transmit and/or store flow measurement related data. In an embodiment of the invention, the acoustic part of the arrangement comprises means arranged to generate an impulse to be sent into the measuring section of the arrangement by a first transducer. In an embodiment of the invention, the acoustic part of the arrangement comprises a second transducer means to obtain said impulse for a response to be formed and processed by a processor.
According to a very simple embodiment of the invention, for the reflector can be embodied as a sudden expansion of the flow tube cross-section. Thereby it is recommendable that the expanded tube section is lined with sound absorbing material in order to prevent any sound or echoes from further sections of the flow tube from entering the measurement section.
In a flow arrangement according to an embodiment of the invention, the measuring section can be provided with an inner wall dividing the measuring section in two parallel channels preferably of the same cross-sectional area. However, in an embodiment of the invention the cross-sectional areas are not limited only to the same values. In an embodiment of the invention the dual channel section comprises a piece of a special acoustic transmission line, capable of conveying both common mode sound and differential mode sound. Common mode sound is a combination of identical sound waves in both parallel channels. Differential mode sound is a combination of a certain sound wave in one channel and its contra sound in the other.
According to an embodiment of the invention, this differential mode sound combined with directional decomposition provides an ideal choice for flow measuring, because it is almost perfectly reflected from the cross section of the sudden end of the inner wall. Moreover, disturbing sound from outside of the measuring section gives a contribution to the common mode sound only, not in the differential mode sound. Accordingly a differential sound source, generating certain sound into one of the channels and its contra sound into the other, should be chosen. Likewise, and more importantly, according to an embodiment of the invention, differential sound sensors, insensitive to common mode sound are chosen.
According to an optional embodiment of the invention, a dual acoustic transmission line, can be used with a differential mode sound for determining the sound velocity from the resonance frequency of a resonator first, as known as such from http://www.acoustics.org/press/150th/Garrett.html, but in the embodiment of the invention, intricate phase measurements are made in a combination of a multi-frequency sound emission in the dual line.
The embodiments of the invention are now explained in detail, by referring to the following Figs. illustrating examples of the embodiments of the invention having thereby no intention to limit the invention by any means. In the Figs.,
a illustrates a flow channel geometry to be used for an embodiment of the invention,
b illustrates a flow channel geometry to be used for an embodiment of the invention,
a illustrates an example of a cross section of the measuring section according to an embodiment of the invention,
b illustrates another example of a cross section of the measuring section according to an embodiment of the invention,
c illustrates a further example of a cross section of the measuring section according to an embodiment of the invention,
If otherwise not indicated for a particular Fig. or Figs., the same reference numerals are used in the different Fig(s). to indicate the same kind of parts, although the parts are not necessarily to be exactly the same. The various embodiments are shown as examples form the various embodiments of the invention, which shown as such are combinable in suitable part.
Let us first begin from the well known expressions for the downstream transit time TD and upstream transit time TU of sound propagating over flight path of length b in a flowing medium of sound speed at rest c and mean flow velocity v:
c and v as a function of TD and TU are thus
By defining the transit time sum as TSUM=TU+TD and the transit time difference as TDIF=TU−TD,
equation pair (2) can rewritten as
From the measurement point of view there is a significant difference between expressions (2) and (3). Extracting the best estimate for TDIF directly from the measurement data, as in the preferred embodiments of the invention, instead of extracting separate estimates for TD and TU is expected to lead to a better estimate for v, too. Besides TDIF one has, of course, to be able to determine TSUM also. The transit time difference manifests itself in the phase difference at well-defined multiple frequencies of the two cross-power spectrums related to the downstream and upstream propagating sound waves. For this purpose, instead of emitting sound as bursts, in an embodiment of the invention, stationary multi-frequency sound as periodically repeated pseudorandom sequences are chosen so allowing precise phase measurement at a large number of discrete frequencies simultaneously. This choice also allows signal averaging, meaning that sensor signals from several consecutive periods are sampled and accumulated together in synchronism with the updating steps of the sound emission, before the signals are analyzed. Thereby the signal to noise ratio can be substantially improved.
A complex multi-frequency phase factor multiplied with appropriately chosen statistical weight vector represents a generalized cross power spectrum vector. Its Fourier transform appears as a generalized correlation function, the centre of which is marks the transit time in question. In what follows transit time determination is most often formulated as a multi-frequency phase measurement. The alternative correlation function approach is always implicated even when not explicitly stated.
Referring to
The sound propagation back and forth over a flight path b in the frequency domain is represented by a propagation vector P, the components of which can be expressed as
P
k(1)=i {−γb+ikTSUM(b)πN−1}, (4)
where k is the index of the frequency channel, γ is a single absorption parameter or a group of parameters for allowing some frequency dependence and N is half the number channels in time domain vectors. The argument (1) of Pk(1) should be understood as (1·b) implying that relative path length 1 corresponds to the actual pathlength b.
The sound reflection at the reflector is represented by a reflection vector R, the frequency components of which are expressed as
R
k=ρexp(iβ), (5)
where ρ is a single echo amplitude parameter or a group of parameters for allowing some frequency dependence, and β is a single echo phase parameter or a group of parameters for allowing some frequency dependence. In particular, β should contain a term, proportional to k and/or to k2, for taking into account that the precise location of the reflecting plane may be frequency dependent.
Both FM1 and FM2 can now be expressed as the superposition of the primary sound wave vectors FM1D and FM2D and the echoes FM1U and FM2U, as returned by the reflector
FM1=FM1D+FM1U=[1+P(1+ζ)R]FM1D
FM2=FM2D+FM2U=[1+P(ζ)R]FM2D (6)
The distance of the reflector from the second sound sensor is assumed to be ζb. The sub indexD refers to the primary sound wave and U to the echo, implying that the primary sound wave propagates in downstream direction. However, length-wise flow direction, upstream or downstream, is quite irrelevant, and thus does not limit the embodiments of the invention. The decomposition of both sensor signals into primary wave and echo vectors is now given by
FM1D=[1+P(1+ζ)R]−1FM1FM1U=FM1−FM1D
FM2D=[1+P(ζ)R]FM2FM2U=FM2−FM2D (7)
So far the expressions (7) contain quite many unknown parameters, most importantly TSUM, others related to sound absorption and reflection. Phase factor vector
ph=arg(FM1U
versus frequency channel index k, represents a straight line, the slope of which is proportional to TDIF, if correct values for the parameters have been chosen. The unknown parameters related to absorption and reflection, are determined from the best fit of ph to a straight line model at zero flow velocity. As flow velocity is varied, only TSUM and to a lesser extent the absorption parameters are supposed to change, while the rest of the parameters keep their values determined at zero flow.
For optimum data fitting one has to choose a good weighting scheme for depressing the importance of lower precision channels relative to the higher precision channels. Expression (9) shows one good choice for a weight vector w to be implemented as w2 multiplier vector in the least square error sum for finding the best fit.
An alternative choice for w is obtained from the sound source spectrum and the propagation and echo models without the knowledge of measured sensor signals. As it can be deduced from equation (3), improving the relative accuracy of TSUM beyond the relative accuracy of TDIF makes no big difference in the accuracy of v. Accordingly, TSUM does not have to be among the parameters to be obtained by the best fit of ph, but can be obtained from any other procedure leading to reasonable accuracy. One possibility is to implement directional filtering by forming a downstream vector and an upstream vector by adjusting the amplitudes of FM1 and FM2 for taking into account sensitivity imbalance and the sound absorption from one sensor to the other, by time-shifting these vectors relative to each other, and by subtracting them so that the result is a pure downstream or a pure upstream vector referred to the same point of definition. The time shift between these two directionally filtered vectors corresponds to the flight time from the point of definition to the reflector and back, which is proportional to TSUM. Here again the time shifts needed for filtering are TU and TD, which are not known until both TSUM and TDIF are known, leading to an iterative process, according to an embodiment of the invention. Each new value for TSUM can be introduced for obtaining a new value for TDIF from fitting eq. (8), leading to new values for TU and TD, and so on.
According to another embodiment of the invention, another possibility is to implement another reflector at the opposite end of the measuring section as indicated by the dashed lines 2a and 2b in
ph
DIF
=arg(FM2
ph
SUM
=arg(FM1FM2
associated with TDIF and TSUM over the flight paths from the sound source to both sound sensors. In this case the reflected signals at both ends are retarded relative to the unreflected signals by the sum of upstream and downstream flight times from the sensor to the reflector and back. Therefore, no explicit echo model beyond the assumption that the same model applies at both ends is needed. TDIF and TSUM are obtained from the slopes of the phase vectors versus frequency channel k, by implementing a proper weight vector for emphasizing the higher precision frequencies against the lower precision ones. TSUM can also be obtained from the set of resonance frequencies corresponding to an odd number of wavelengths from one reflector to the other one and back, possibly corrected for the exact phase change at the reflection.
So far no explicit description of the sound reflectors has been given. The basic version of an embodiment of the invention is exemplary embodied in
As an alternative to the reflector/attenuators shown in
According to an embodiment of the invention, in which a well-defined reflector is desired at one end of the measuring section only, the inner wall should not stop at the other end abruptly, but in a gradual manner.
In an alternative embodiment of the invention, a differential sound source or sensor are replaced by a pair of matched ordinary sound sources or sensors, coupled to operate in opposite phase.
Apart from the basic dual channel acoustic transmission line, with two sub-channels, as previously embodied in the
As referring to
According to an embodiment of the invention the sound cancellation does not have to be the sum of two opposite phase waves, but can be accomplished from three, four, or more waves of the same amplitude, but their phases differing by 360/3, 360/4 or 360/(n>4) degrees respectively.
In an embodiment the cylinder formed into the measuring section in
said low frequency sound is emitted (602) in a second phase (+) into the second sub-channel (CH2) of the flow channel,
returning (603) an echo back into the measuring section, from a sudden at least at one end (2a, 2b) discontinuity location of said wall (5), said discontinuity location arranged to operate as a reflector,
detecting (604) said emitted sound from said sound source (3) at a distance (b), by at least one sound sensor (4a, 4b) between the sound source and a reflector end (2a, 2b) and
determining (605) the sound velocity at rest (c) and the flow velocity (v) from the known distance of the sound sensors and the obtained values for sum (TSUM) and difference (TDIF) of the upstream and downstream transit times of sound between the sensor locations.
However, the embodiments with several channels having more than two sub-channels may not bring benefits over the basic dual transmission line, in a normal flow measurement determination. But they can be useful for instance in cases in which redundant measurements are needed at a different frequencies, or the measurement data is needed very rapidly in a short duration measurement for a momentarily flow rate values in a statistically more sufficient grounds for averaging for instance. Such conditions could be at least in theory met in conditions in which chemical reaction or other reason that causes sound absorption in the channel at certain frequencies are expected to occur. Utilisation of such embodiments of the invention, where the measuring section is implemented with the help of many parallel sub-channels, can be particularly useful in large tubes or pipes, or arrays thereof, for providing using high frequency sounds and/or achieving a desired accuracy within a short measuring section.
The digital processor, auxiliary electronics and appropriate temperature and pressure sensors for generating, measuring and/or processing the data into flow readings, to be used for the flow measurement, are not shown in the Figs. for clarity reasons, although are serving a multitude of various needs of various embodiments of the invention.
Number | Date | Country | Kind |
---|---|---|---|
FI2005/000288 | Jun 2005 | FI | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/FI2006/000212 | 6/16/2006 | WO | 00 | 8/13/2008 |