Abstract:
A method of measuring the flow rate of a fluid through a conduit by means of a sonic flow meter employing the steps of measuring the speed of sound transmission through the fluid from an upstream point in the conduit to an opposed downstream point, measuring the speed of sound transmission through the fluid from the downstream point to the upstream point, subtracting the two measurements to obtain a liquid flow rate indication, adding the two measurements to obtain the fluid sound velocity, measuring the temperature of the fluid in the conduit, determining a measurement of the fluid viscosity from the fluid speed of sound transmission and the temperature, and correcting the determined fluid velocity utilizing the detected fluid viscosity.

Description:
BACKGROUND OF THE INVENTION 
     1. Field of the Invention 
     This invention relates to the field of sonic flow meters of the type particularly described as ultrasonic, single-path, transit time flow meters such as described in U.S. Pat. No. 3,720,105. The invention is particularly directed to the field of ultrasonic flow meters having computational apparatus for correcting sonically determined flow rates for viscosity, to provide corrected flow rates. 
     2. Description of the Prior Art 
     A convenient and commonly used means of measuring liquid flow is by the use of sonic flow meters. These are devices which transmit sound pulses through the fluid flowing in a conduit. In the most commonly practiced means of employing ultrasonics for measuring fluid flow rate, a sonic path is provided between an upstream and downstream location spaced on opposite sides of a conduit. In the preferred method of practicing the invention the rate of travel of the sound through the liquid is measured in both directions, that is, with the fluid flow component and against the fluid flow component. By subtracting such measurements, which effectively cancels out the rate of sound travel through the fluid itself, the rate of the fluid flow can be accurately determined. 
     The application of ultrasonic flow meters to measurement of such products as water which change little in viscosity, has been successfully employed. Difficulty, however, is experienced in the use of ultrasonic flow meters to measure volumetric flow rates of fluids which vary considerably in viscosity. This is particularly a problem in attempts to utilize ultrasonic flow meters in the petroleum industry wherein the viscosity of the hydrocarbon fluids varies over a great range. 
     An object of the present invention is to provide means of utilizing ultrasonic flow meters for measuring the flow rates of fluids having variable viscosities. 
     Another object of this invention is to provide means employing sonic measurements to determine the viscosity of a fluid flowing in a conduit. 
     Another object of this invention is to provide means employing empirically derived relationships for calculating corrected flow rates of fluids wherein the measurements required for determining the flow rate are supplied as upstream and downstream sonic measurements along with the fluid temperature measurement. 
     These general objects as well as other and more specific objects of the invention will be fulfilled in the following description and claims, taken in conjunction with the attached drawings. 
     SUMMARY OF THE INVENTION 
     A method is provided for measuring the flow rate V of a fluid through a conduit having an inside dimensional constant D and a fluid temperature T employing an ultrasonic flow meter. The method is accomplished by measuring the speed of sound transmission through the liquid from a downstream point in the conduit to an opposed upstream point separated by a sonic path length L to provide a value f A . Next, the speed of sound transmission is measured through the fluid in the opposite direction, that is, from the upstream point to the downstream point to provide a value f B . By adding f A  to f B , a value C is obtained which is indicative of the fluid sound velocity. A corrected fluid sound velocity C&#39; is obtained by the formula 
     
         [L(f.sub.A +f.sub.B)]/2N 
    
     where N is a multiplier factor. Subtracting f A  from f B  provides a value V from which the fluid velocity V&#39; can be calculated. The kinematic viscosity V k  of the fluid is calculated utilizing values C&#39; and T and an empirically derived algorithm. The Reynolds number R e  of the fluid is calculated by the formula DV/V k  and a constant K 2  is derived from the Reynolds number employing an empirically determined relationship. A constant K 1  is derived for the formula 2ND/L 2 . The uncorrected fluid flow velocity V&#39; is derived from V/K 1 , the corrected flow rate V&#34; by the formula V&#39;/K 2 . 
    
    
     DESCRIPTION OF THE DRAWINGS 
     FIG. 1 is a block diagram of the components employed in the invention and showing diagrammatically the measurement of sound along a sonic path in a conduit. 
     FIG. 2 is a logarithmic graph showing the typical relationship between the constant K 2  and the Reynolds number for a variety of standard pipe sizes. 
    
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
     In the drawing and the specification the following nomenclature is employed: 
     VCO=Voltage Controlled Oscillator 
     f A  =Upstream VCO frequency 
     f B  =Downstream VCO frequency 
     T=Liquid Temperature 
     C=Velocity of sound in the liquid as indicated by a sum of f A  and f B . 
     C&#39;=Velocity of sound in the liquid corrected for the parameters by which f A  and f B  are measured. 
     V=Velocity of the fluid as indicated by f B  -f A . 
     V&#39;=Velocity of the fluid corrected for the parameters by which f A  and f B  are measured. 
     V&#34;=Velocity of the fluid corrected for variations in fluid velocity and viscosity. 
     D=Dimensional constant derived from the conduit cross-sectional shape and size and the angle of acoustic transmission direction with respect to the conduit axis. For a circular conduit and a 45° angle of acoustic transmission, D is equal to the inside diameter of the conduit. 
     L=Length of the sonic path 
     R e  =Reynolds number 
     K 1  =Constant used to correct the sonic meter measurement V into units of flow based on physical conditions. 
     K 2  =Constant used to correct the sonic meter measurement V&#39; for fluid velocity and viscosity. 
     Referring to the drawings, the basic apparatus for practicing the invention is illustrated in block diagram and schematic form in FIG. 1. The object of the invention is to determine the flow rate of the fluid 10 flowing through conduit 12. The invention is particularly directed towards such measurement when the fluid 10 is of the type having variable viscosity, such as hydrocarbon fluids. Positioned in conduit 12 is a first, upstream, sonic energy transmitting and receiving station 14. A second sonic transmitting and receiving station 16 is placed downstream of station 14 and on the opposite side of conduit 10 from station 14. Between stations 14 and 16 a sonic path of length L is established. 
     The apparatus including sonic stations 14 and 16 is well known and is of the type commonly employed in sonic flow meters, and for more information as to the structure and operation of such devices, reference may be made to U.S. Pat. No. 3,720,105. 
     Affixed to conduit 12 is a temperature measuring device 18 by which the temperature T of the fluid 10 is detected. 
     While various means may be employed for transmitting and receiving sound pulses between stations 14 and 16, a commonly employed arrangement is the use of voltage controlled oscillators (VCO) 20 and 22. A measure of the speed of sound transmission along the path L between the second sonic station 16 and the first sonic station 14 is given by voltage controlled oscillator 20 which provides a signal output indicated by f A , that is, the upstream voltage controlled oscillator frequency. In like manner, a measure of the speed of sound transmission between upstream station 14 and downstream station 16 is provided by voltage controlled oscillator 22 which provides an output f B  representative of the downstream voltage controlled oscillator frequency. The signal f A  is proportional to the speed of sound transmission through liquid 10 from downstream point 16 in the conduit to the opposed upstream point 14. The signal at f B  is proportional to the speed of sound transmission through the fluid 10 from the upstream point 14 to the downstream point 16. The signals f A  and f B  are added in the circuit 24 providing, at the output, a signal C which is indicative of the fluid sound velocity, that is, the speed of travel of sound in the fluid irrespective of the direction or speed of movement of the fluid. The signal C is treated in circuit 25 to correct for the physical parameter of measurement by the formulae LC/2N, to provide the corrected sound velocity signal C&#39;. By subtracting f A  from f B  in subtraction circuit 26 a signal V is provided which is dependent on the fluid flow velocity. In effect, the subtraction circuit 26 cancels the effect of the speed of sound transmission of the fluid itself leaving only the component attributable to the movement of fluid 10 through conduit 12 and points 14 and 16. 
     In the computational circuitry of FIG. 1, process constants are inserted. D is a constant representative of size and dimensions of conduit 12 and the angle of the sonic flow path. L is the sonic path length between points 14 and 16. These values may be inserted into the circuitry, either in analog or digital form, by the circuitry components 28 and 30. An additional process constant is a multiplier factor N which is typically selected as a function of sound velocity C and transducer spacing L to yield the largest possible values for f A  and f B  within the limitations of VCO&#39;s 20 and 22. This corrective multiplier is placed into the circuitry utilizing component 32. 
     An important aspect of the invention is the method by which the kinematic viscosity of the liquid 10 may be determined and the use of the determined kinematic viscosity in correcting the flow velocity to provide a viscosity corrected fluid flow rate measurement. 
     Circuitry 34, which may be in the form of a programmed computer chip using an empirically derived algorithm, is employed for computing the kinematic viscosity of the liquid. The input into circuitry 34 is the sound velocity C&#39; of the fluid and the temperature T. It has been unexpectedly discovered that fluid kinematic viscosity is relatively linear on a logarithmic scale as a function of the speed of sound transmission of the fluid and of the temperature. It has been determined that when the fluid 10 is a hydrocarbon type liquid, such as crude oil, gasoline, propane, diesel fuel, etc. an algorithm employed in circuitry 34 may take the following form: ##EQU1## where K 3 , K 4 , K 5  and K 6  are empirically determined constants, the values of which depend on the nature of the fluids to be metered and units employed in V k , T, and C&#39;. 
     A commonly employed method of correcting a flow meter of the sonic type employs the Reynolds number of the liquid. This may be obtained in circuitry 36 using the formula DV&#39;/V k . That is, the Reynolds number is the ratio of the measured flow rate divided by the kinematic viscosity of the fluid. Since the flow of fluid through a conduit is not at the same rate at all distances from the inside wall of the conduit due to the drag imposed by the conduit walls, a constant is required to reflect the average velocity along the line between the sonic transducers at points 14 and 16 versus the average velocity of the fluid in the cross-section of the conduit. This is the constant K 2  and is obtained from the Reynolds number R e . This information in the calibration of a meter reading is normally taken from a chart supplied with each sonic flow meter model. 
     FIG. 2 shows a typical empirically derived chart for a sonic flow meter wherein the constant K 2  is the ordinate and the Reynolds number R e  is the abscissa. This chart has been plotted for a sonic flow meter where conduit 12 is a standard pipe size with calibrations made for pipes of 3&#34;, 6&#34;, 12&#34;, 24&#34;, 48&#34;, and 96&#34; nominal pipe size. For each sonic meter the computational circuitry of element 38 provides the correct output K 2  for the Reynolds number input. 
     To correct the measured flow rate signal V for the physical parameters employed in the measurement, a constant K 1  is required. This constant is obtained in circuitry 40 by the formula 2ND/L 2 . 
     Circuit 41 corrects the detected flow velocity V by the formula V/K 1  to obtain V&#39;. The value V&#39; is fed to circuit 36 for use in calculating the Reynolds number R e  as previously described, and in addition, is fed to circuit 42 where the fluid velocity V&#39; is corrected by the formula V&#39;/K 2  to obtain the final fluid velocity V&#34; which is fully corrected for physical measurement parameters and fluid viscosity. 
     All of the circuitry required can be accomplished using state of the art technology and microcomputer processes which are not a part of the present invention since they are well within the skill of the art. 
     The invention has been described and illustrated based on the assumption that conduit 12 has a circular internal cross-section. While this is the typical environment in which sonic flow meters are employed, it can be seen that the same principles of the invention are applicable to conduits 12 having square, rectangular, elliptical, or any other shape cross-section. 
     While the invention has been described with a great degree of particularity it is manifest that many changes may be made in the details of construction and the arrangement of components without departing from the spirit and scope of this disclosure. It is understood that the invention is not limited to the embodiments set forth herein for purposes of exemplification, but is to be limited only by the scope of the attached claim or claims, including the full range of equivalency to which each element thereof is entitled.