THERMOMETRIC FLOW MEASUREMENT

Abstract: According to the present invention, two temperature meters are mounted externally on a pipe within which fluid flows. The temperature meters deliver electrical signals to a computer or corresponding device. The mutual spacing between the temperature meters is known. The temperature signals are sampled a number of times, m, and a transfer function is calculated, this function having n terms, where m>n, so that the least square method can be applied. Preferably, it is assumed that the elements of the transfer function have small second difference significance, that the first and the last elements have low values, and that all elements are positive.

Description: Thermo etric Flow Measurement The invention relates to a method according to the preamble of Claim 1, and to apparatus for carrying out the method. There are many situations in which it is necessary to be able to measure the flow of fluid in a pipe, in the absence of permanent installations. For example, there may be a need to measure the flow of a medium in dis¬ trict heating networks at any one of a number of acces¬ sible points, in order to be able to localize a suspect¬ ed leakage. There may also be a need to measure the flow of medium in a central heating pipe in a building structure, because of a suspected short circuit or a stoppage in the system. There are many examples of different situa- tions in which there is a need to measure the flow of medium in a pipe or conduit, without performing work on the actual pipe or conduit itself. DE-A-3 423 966 teaches a method of measuring the flow of a medium in a conduit with the aid of two sequentially mounted, quick-response thermometers. In this case, there is generated a temperature-change front which travels through the conduit, for insta ;e by causing an upstream heating element to generate a heat pulse, or by changing the setting of a mixer valve. When the distance between the thermometers is known, the flow rate of the medium can be calculated, by measuring the time lapse between the occurrences of the front at respective temperatures. GB-A-2 159 631 teaches a measuring device which is based on similar principles, although in this case one thermo¬ meter is replaced with a quick-response heat source, the activation of which constitutes a starting point and wherein a signal from a downstream thermometer is de- tected in order to determine the time of flight or tran¬ sit time. Sϋ-A-710 004 teaches an arrangement in which two thermo¬ couples are mutually spaced in the path of a liquid medium flowing in through conduit. The thermocouples are connected in series and the frequency spectrum of the resultant composite signal is analyzed. This is said to show a sharply defined minimum for a given frequency, from which the rate of flow of the medium can be estab- lished. US-A-4 604 904 deals with measuring the correlationship between different parameters that occur spontaneously along a pipe (capacitance, ultrasonic, optical deteσ- tion, thermodetection, electrical charging, conductivi¬ ty), and proposes that a plurality of σapaσitive measur¬ ing sensors are placed along the pipe at distances from one another which will enable the terminal length of the measuring apparatus along the pipe to be compensated for. One industrially used, standard flow meter is based on generating a homogenous magnetic field through an insu¬ lating pipe, in the transverse direction of the flow. According to Faraday's law, a voltage will then occur perpendicularly to both the field and to the flow direc¬ tion. Such meters provide a result which indicates the volumetric flow, independent of the nature of the flow. Such flow meters, however, are both complicated and expensive. It is an object of the present invention to provide a method and apparatus for measuring the flow of a medium, which, with respect to the quality of the measurements, is able to compete with magnetic measurements, but which is simpler, less expensive and easier to apply or to install, but which will nevertheless provide a satisfac- tory result. It has been found, surprisingly, that this object can be achieved by effecting correlation measurement with temperature sensors which are mounted in thermal contact with the outer surface of a standard pipe, wherein the only work required is that of ensuring that the pipe surface is sufficiently clean to provide a good contact. In certain cases, it may be appropriate to form a simple recess or shallow in the pipe surface. The signals obtained from the sensors are sampled and the sampling results stored. The correlation between the signal- series for different time shifts is determined and the highest correlation factor chosen. The time shift corresponding to this highest correlation factor then constitutes, together with the distance, an expression for the velocity of the flowing medium. When the cross- sectional area of the conduit or pipe is known, the rate of flow can be determined. When making a comparison between the measurements with the aid of a magnetic flow meter, it was found that extremely good agreement is obtained with highly satis¬ factory linearity. When measuring correlation on spontaneously occurring temperature fluctuations, by detecting the fluctuations from outside the pipe or conduit, it has been found necessary to work with relatively long sampling series, in order to obtain an overdefined system of measurements which can be processed in accordance with the least squares method. It is not possible to effect a simple transit time measurement, as when generating a gradient artificially and measuring with the aid of sensors mounted in the flow of, and consequently it is necessary to work with statistic correlation methods. With the exception of extremely simple cases in which rapid variations occur, it is also necessary to work with a calculated transfer function or transverse correlation whose maximum value will give the desired transit time measurement, while using physically reasonable secondary conditions. The invention wilα. now be described in more detail with reference to an exemplifying embodiment thereof and also with reference to the accompanying drawings. Figure 1 illustrates a system for measuring the flow of medium in a pipe, in accordance with the invention. Figure 2 illustrates diagramiαatically an amplifying and filtering chain for obtaining a temperature signal. Figure 3 is a block schematic of a two-channel system constructed in accordance with Figure 2. Figures 4 A-D illustrate examples of the correlation result measured and calculated with different secondary conditions. Figures 5 A-E illustrate examples of comparisons between measurements carried simultaneously with magnetic flux intervals and in accordance with the invention, in which the fluid was a liquid, namely water. As illustrated in Figure 1, the inventive measuring process is effected with the aid of two measuring heads 1, 2 which are mutually spaced along a pipe. The measuring result is obtained from a central processing unit 3 which includes an indicator window 4. The mea¬ suring heads are appropriately given different colours, for instance red and blue, so that the meter will show when the current passes from "red to blue" or vice versa. Figure 2 illustrates a measuring circuit which functions to measure the variations in temperature that occur externally along the pipe. A thermistor T is located adjacent the measuring head. A resistor R is connected, together with the thermistor, between earth and a con¬ stant voltage. The voltage of the junction is the magnitude whose variations are to be studied. Although the measuring head may incorporate more or less all of the electronics, it is preferred at present for the measuring head to incorporate only the thermistor, said measuring head being connected to the remainder of the circuit by means of a screened, twisted pair of conduc- tors, so that the i pedence-adapting amplifier 11 and the resistor R 1 will be spaced from the thermistor. High frequency disturbances are eliminated by a lowpass filter R2C1. The variable resistor R3 is provided in order to bring the voltage variations generated by the temperature variations to a suitable scale. The time constant R3C2 may be given relatively long duration, e.g. a duration of 100 s. The amplifiers 12, 13 and 14 form a pushpull coupled stage together with the amplifi¬ ers 15, 16 and 17, whereas R4C3 eliminates high frequen- cy signals and an output signal can be collected from the output of the amplifier 18. The various amplifiers may be standard amplifiers (e.g. type 741). The output signals are then sampled with a sampling frequency of 1- 50 Hz. For practical reasons, it is suitable to construct the system in a slightly more complicated way, since it is desired to automate the setting of the resistor R3, so that even the most signi icant bit is preferably filled in an AD-converter. In the block schematic shown in Figure 3, the stages A-E in Figure 2 are corresponded by duplicated boxes, so that signals from the thermistors can be taken sequentially to the computer MPϋ. The Analog Mux 1 of a first multiplex circuit delivers analog signals to the analog-digital-converter. The Analog Mux 2 of a second multiplex circuit adjusts the resistor R under the control of the computer, so that the converter will neither obtain a signal which is so large as to result in data overflow, nor a signal which is so small as to reduce the number of significant bits. By controlling the remaining units (Mux 1, Mux 2, Mux 3), the signal unit MPU is able to detect the volt¬ ages across the thermistors, for instance each second, and to take up, e.g., 1024 values each with 8 bits. It is then suitable to perform the relevant calculations for one such time series and to take in 10 new values, while eliminating the 10 oldest values, and to perform a new calculation. # CALCULATION OF THE TRANSFER FUNCTION Conventions: W column-matrices: upper cases in bold. w. lower case letters with indexes: matrices ele¬ ments. H remaining matrices: in bold and underlined. Signs: M number of points in X and Y respectively. m - M-l. Tl, T2 time series with temperature data. X, Y Tl and T2 respectively differentiated twice (*£, y^ 0<i<m). N number of points in A. n = N-l. A auto-covariance-function for X (a., -n≤i≤n). 1 C cross-covariance-function between X and Y (c, 0<u<). 1 H transfer function (h. , o≤i≤n) . 1 #.1 PROCESSING OF MEASURED D ^ The time series Tl and T2 respectively are obtained by measuring the temperatures at locations l and 2 at discrete points in time. The auto-covariance and cross- covariance-functions can only be calculated from static time series, i.e. time series in which the waiting value is independent of time. Since temperature data measured over a relatively short period of time (some minutes) may well include trends, and thus not be static, it is necessary to differentiate data prior to analysis. We differentiated twice with a "soft" differentiation method which did not only take into account the nearest lying points: Tl . = — (T +2T . -2T . -T ) i 8 l 1+2 i+1 i-1 1-2 ' All analyses have been made on the time series 0 1 2 „, y, , y„, ... y ), where 0 l J 2 m . , y. were calculated in accordance with 1 1 X. = T2. 1 1 #.2 THE TRANSFER FUNCTION Our basic assumption is that Y is linearly dependent on X in accordance with the transfer function H. The following relationship applies when disturbances can be ignored X *i= ∑ h k x d _ k (H-def) ic=0 In the transfer function H, the truth that h. equals 0 when i<0 will apply for each physical system. The opposite would imply that the system is responsive to future events. We also assume that the effect of a separate event decays and looses significance after a given period of time, i.e. h.=0 when i>n. i #.3 AUTO- AND CROSS-COVARIANCE In the case of a time series, the term variance for a stocastic variable can be broadened to include the expression auto-covariance-function. In addition to the variance (=a ), the expression also includes certain information concerning possible periodicity, among other things. As with the case when calculating mean values and vari¬ ances, we are unable to produce the precise auto-.covari- ance-function; it must be estimated. The formula is 7 r+J a,- = X X ■ kι -i (A-def) r where r and j can be randomly selected. In order to be able to calculate an auto-covariance, it is necessary for the time series to be static. A static time series implies that the expectance value of a. is dependent on j in the above definition. Note that a . = a.. -1 1 In the case of two time series, there is defined the cross-covariance-function which provides information relating to the coupling between two time series. The cross-covariance, C, is estimated with the aid of the formula Ci - -|∑ ϊikXjt-i (C-def) #.4 LEAST SQUARES-ADAPTATION OF H When X and Y are known and H is sought, equation (H-def) provides an overdetermined equation system n n-1 l n x n+_.1, x π n+l n m m-1 m-n m or X H = Y (H-over) Because of disturbances in X and Y, there is no transfer function which will precisely solve equation (H-over), although the best approximation is found through the least squares method. The present demonstration is performed on the basis hereof. The problem is one of finding the matrix H which mini¬ mizes RES, the sum of the squared residuals. From the equation (H-def): in matrix form: RES = (X H - Y) T (X H - Y) = T = (X H) (X H) - (X H) Y - Y ( X H) + Y Y = m m m = H X X H - 2H Y + const. (RES) where the rule for the transponat of a matrix T T T ( B) ■=!& has been used to obtain the last line. This sum has its minimum when a RES β 0 for all i/g^ 0<i≤n at, a RES 3Λ, = (tfV.1) i + (x*z m i - 2 ∑ T Y) where (...) 1. denotes the i:th element in the vector. These N equations together give the well known formula for the least square method: Z*Σ B = Z*Y (H-from-XY) A closer examination of the matrices X TX and XTX show that the matrices elements are comprised of elements from the auto-covariance-functions and the cross-covari¬ ance-functions. m (X T Y ) i = ∑ y k χ k -i = ( m+l- ) c i -*Σ r Y = (m+l-n) C m (X T X) itj = ∑ k .i x k . j = (m+l- ) .a i-J k=n -> X X = (m+l-n) where a„ a „ a „ a 0 -1 -2 -n a a a a 1 0 -1 1-n A * = a a a ..... a Λ 2 1 0 2-n • • ■ • * • • • • • • • • a a , a Λ a„ n n-1 n-2 0 The equation can thus be written as H = C (H-from-AC) #.5 DIRECT DERIVATION FROM THE DEFINITIONS The equation (H-from-AC) can thus be derived directly from the definitions (H-def), (A-def) and (C-def). These relationships can be written in a simpler manner, with the aid of the convolution operator, *: Y = H * X A * X * X C = Y * X The cross-covariance-function can then be written as C = Y * X = = (H * X) * X = = H * (X * X) = = H * where we have utilized that the convolution operator is associative, i.e. (i*g)*h=f*(g*h). This expression can be rewritten in surama form. n or in matrix form, C * - *** h H which is identical to equation (H-from-AC) from the preceding section. #.6 CALCULATING H WITH THE AID OF THE EQUATION SYSTEM The equation system (H-from-AC) may, of course, be at¬ tacked directly with standard methods in order to re¬ solve H. When this was carried out in one experiment (pipe D, v=0.88 m/s) , hi, 0<i<n according to Figure 4A was obtained. The H obtained in this manner gave no useful information whatsoever. Disturbances in measure¬ ment data dominated totally. The following values are those which can be expected of the transfer function for this system: 1. h.≥O for 0<i<n. 1 2. h. small for i=0 and i->n. 1 3. A quiet sequence of the h-values. 4. The h-values grow monotonously towards a maximum and then decrease monotonously towards zero. In order to obtain a physically reasonable transfer function, we solve the equation system (H-from-AC) with the secondary condition that the sum of h:s second- differences shall be small. This condition provides a "quiet" sequence of the h-values. The difference can be written as D. , .=h. , -h. l+l ,l l+l l The second order difference will then be D i+l,i ~ °i,i-l (h ι + χ- v - (h - - - - h. , - 2*h. - h. , 1+1 1 l-l The form factor is the sum of the second order differ¬ ences in the square 12+1 r> FORM = ∑ L ) = i=-l n+1 where 6 -4 1 0 0 0 -4 6 -4 1 0 0 1 -4 6 -4 1 0 Q = 0 1 -4 6 -4 0 0 0 1 -4 6 0 0 0 0 1 -4 6 The equation (RES) is modified by adding the above form factor, multiplied with a weighting factor, β. The new expression is + const. As before, the problem is one of finding the H which minimizes RES. Analogously with the above derivation (RES) - (H-from-XY) there is obtained the equation system (X T X + β Q)H = X T Y which for a static time series can be replaced with ( + β <J)H = C (H-from-ACQ) Figures 4B-C illustrate analyses of the same data as Figure 4A, although with the above equation. The Fig¬ ures show analyses with different values of β: 0.001, 1.0 and 100. A clear improveraeat is already obtained when β==0.001, and when β=100 the result is clearly acceptable, even though the configuration of the curve is still not exactly what one would expect from a physi¬ cal aspect. 4.7 THE EQUATION IN THE PROGRAM HCOV Originally, however, another equation is ierived for determining H. With a starting point fr t equation (H-from-AC), derived directly from the definitions ac¬ cording to section #.5, there was obtained RES = + α Q)H - 2 H T £ £ + const. The H which minimizes RES is then given by the equation (& T & + a Q)H = Λ (H-HCOV) It can be shown that this expression is closely related to equation (H-from-ACQ) . When multiplying both sides T -1 with (A ) from the left there is obtained T —1 T T —1 T ( ) (h + α Q)H = (A ) TC -> (A + β (A T ) *"1 Q)H = C which shows that equation (H-HCOV) used in our analyses T -1 only differs from equation (H-from-ACQ) in that (A ) Ω is used instead of Q. #.8 DETERMINING HEAT PULSE VELOCITY The purpose of the analyses is to obtain a measurement of the time required for a heat pulse to travel through a given distance in the pipe. In order to achieve this, we adapt a second degree polynomial to the five points in H which lie nearest the top, Figure 4D. λ , i.e. max the λ for which h is maximum (λ is a continuous vari¬ able), can be easily obtained from the curve h(λ) . In order to determine the velocity of the heat pulse, we also require: DTIME the time between two mutually se¬ quential measurements. DIST the distance between the measuring points on the pipe. The calculated velocity is v _ DIST λ^ r DTIME The calculations can be carried out on a computer, for instance a conventional Personal Computer. Since the actual algorithms for numerical calculation are well known in themselves, they can be taken from standard literature, such as Press et al, "Numerical Recipes" (Cambridge University Press, 1988). Comparison tests have been carried out with the aid of the aforedescribed method. The comparisons were carried out with the aid of a flow meter manufactured by Endress & Hauser. A series of flows were arranged in pipes of five different internal diameters, ranging from between 12.5 mm and 35.9 mm. The rates of f -*v were varied from 0.3 m/s to 1.06 m/s, with the Re α ..Id's number varying from between 1.100 to 15.000. The line¬ arity was found to be good and the accuracy was as good as one could expect in view of the fact that the inter- nal diameters of the pipes were not known to desired degrees of accuracy. The graphs shown in Figures 5A-5C have been drawn on the basis of pipes whose respective diameters were 35.9 mm, 27.2 mm, 21.6 mm, 16.0 mm and 12.5 mm (estimated values on the basis of measured external diameters and measured wall thicknesses at the ends of respective pipes). The thermometrically measured values are plotted along the y-axis of the graphs, and it will be seen from the graphs that the linearity is very good and that the measurements are clearly reliable.

Claims: Claims 1. A method for measuring the rate of flow of a fluid through a pipe, comprising measuring temperature changes at at least two different locations at mutually known distances along the pipe, c h a r a c t e r i z e d by mounting at two said positions, without contact with the fluid, temperature sensors which produce electrical temperature signals which are sampled a large number of times m, wherein one samp *ling series, x,1,....xm is taken from the one sensor and y1„ ...y J m is taken from the other sensor, whereafter a transfer function H=h . h ,...h is calculated, where n is smaller than m 0 1 n and the elements of which fulfil the equation wherein H is determined by the least squares method; determining a maximum for the transfer function H, from which there is obtained the number of sampling intervals to which transport of the fluid over the known distance corresponds, from which the mean flow rate is deter¬ mined. 2. A method according to Claim 1, c h a r a c ¬ t e r i z e d in that the overdefined equation system is resolved while using the secondary condition that second order differences between successive elements of the transfer function are small. 3. A method according to Claim 1 or 2, c h a r a c ¬ t e r i z e d in that the overdefined equation system is solved under the secondary condition that the first and the last elements have low values and that all are positive. 4. An arrangement for applying the method according to Claim 1, comprising two temperature sensors (T) , c h a r a c t e r i z e d in that the arrangement also includes a computer, an analog/digital converter, an amplifier having variable amplification (D), whose output signal is coupled to the analog/digital con¬ verter, a first multiplex circuit coupled to the comput¬ er for controlling the variable amplification, such that its output signal will be compatible with the analog/ digital converter, and a second multiplex circuit cou- pled to the computer so as to conduct alternately to the analog/digital converter signals from the two tempera¬ ture sensors to the analog/digital converter.

Inventors: STENLUND LARS EMANUEL, STENLUND, LARS, EMANUEL
Applicants: STENLUND LARS EMANUEL, STENLUND, LARS, EMANUEL
Application year: 1992
Main classification: G01F1
IPC codes: G01F1, G01P5