# Patent

## Title
Method and system for multi-path ultrasonic flow rate measurement

## Abstract

In one embodiment, a multi-path ultrasonic flow meter for determining the flow rate of a fluid in a conduit is disclosed comprising at least two transducer pairs attached to the conduit at two chord locations, one greater than and one less than a mid-radius chord, wherein the composite ratio the two path velocities to the flow rate is substantially constant over the range of Reynolds numbers. In another embodiment, a method of determining the flow rate of a fluid in a conduit is disclosed comprising the steps of determining a composite velocity by determining a weighted average of a plurality of path velocities, determining a chord velocity ratio based on the path velocities, determining a profile correction factor based on the composite velocity and the chord velocity ratio, and determining the flow rate based on the composite velocity and the profile correction factor.

## Background
BACKGROUND OF THE INVENTION
This invention relates generally to ultrasonic flow rate measurement, and more particularly to multi-path ultrasonic flow meter.
Ultrasonic flow meters are used to determine the flow rate (i.e., mean pipe flow rate (Vm)) of a variety of fluids (e.g., liquids, gases, etc.) and speed of sound in the fluid flowing in pipes of different sizes and shapes. Knowledge of the flow rate of the fluid can enable other physical properties or qualities of the fluid to be determined. For example, in some custody-transfer applications, the flow rate can be used to determine the volume (Q) of a fluid (e.g., oil or gas) being transferred from a seller to a buyer through a pipe to determine the costs for the transaction, where the volume is equal to the flow rate multiplied by the cross sectional area (A) of the pipe. In other applications, the speed of sound can be used to determine the mean molecular weight of a fluid flowing in a pipe to optimize and control chemical or combustion processes.
In one type of ultrasonic flow meter employing transit time flow metering, one or more pairs of ultrasonic transducers can be attached to a pipe (or spool piece attached to a pipeline), where each pair can contain transducers located upstream and downstream from each other forming an ultrasonic path between them. Each transducer, when energized, transmits an ultrasonic signal (e.g., a sound wave) along an ultrasonic path through the flowing fluid that is received by and detected by the other transducer. The path velocity (i.e., path or chord velocity (Vp)) of the fluid averaged along an ultrasonic path can be determined as a function of the differential between (i) the transit time of an ultrasonic signal traveling along the ultrasonic path from the downstream transducer upstream to the upstream transducer against the fluid flow direction, and (2) the transit time of an ultrasonic signal traveling along the ultrasonic path from the upstream transducer downstream to the downstream transducer with the fluid flow direction.
The path velocity (Vp) of the fluid averaged along one or more ultrasonic paths between the pairs of transducers does not necessarily represent the flow rate (Vm) of the fluid across the entire cross section of the pipe since path velocity has a radial distribution wherein the flow velocity across the pipe is not necessarily uniform (i.e., all of the fluid in the pipe does not flow at the same path velocity). This radial distribution of the flow velocity is described by its flow profile (or fluid velocity profile), which is a function of several factors, including the Reynolds number (Re). The Reynolds number (Re) is a dimensionless number that characterizes the behavior of fluids flowing through a pipe by providing a measure of the ratio of the inertial force (i.e., flow rate (Vm)) to the viscous force (i.e., kinematic viscosity (v)) of the fluid flowing in a closed circular pipe with a known diameter (D), and is equal to the product of the flow rate and the pipe diameter divided by the kinematic viscosity of the fluid (Re=(Vm*D)/v). Knowledge of the Reynolds number allows for the determination of the profile of fluid flow in the pipe: (1) for a relatively low Reynolds number (e.g., Re<2300), laminar profile (parabolic shape where the maximum path velocity is at the center of the pipe and the flow velocity gradually decreases to zero at the pipe wall); (2) for a relatively high Reynolds number (e.g., Re>4000), turbulent profile (flattened shape where the flow velocity is virtually constant except for an area near the pipe wall); and (3) for an intermediate Reynolds number (e.g., 2300<Re<4000), transitional profile (combination of laminar and turbulent).
In addition to the Reynolds number, which is generally the principal factor affecting a flow profile and is dependent on the kinematic viscosity of the fluid, the flow profile is also a function of several other factors, including fluid characteristics, pipe roughness, pipe configuration (e.g., straight, bends, elbows, tees, valves, headers, etc.), and disturbances caused by the transducer ports. Accordingly, in order to determine the flow profile for a particular ultrasonic flow meter installation, a number of installation-specific pieces of information would be required, including knowledge of several of the physical properties of the fluid and the pipe, some of which may change over time (e.g., relative pipe wall roughness (ε/D)).
Rather than determining a flow profile, prior art solutions using multiple paths (Pi) and multiple path velocities (Vpi), select the locations of the transducer pairs (i.e., the ultrasonic paths) and the weights (wi) to be assigned to each path velocity (Vpi) in determining a mean flow rate (Vm) based on one or more known mathematical techniques (e.g., Gaussian Quadrature, Chebycheff, Lobatto, etc.) as described in, for example, U.S. Pat. Nos. 3,564,950 and 3,940,985.
While these prior art solutions minimize the need to determine some installation-specific pieces of information required to determine the flow profile, the solutions only provide moderate to high accuracy for a flow profile that requires a high order polynomial required to describe it (i.e., the higher the complexity of the flow profile, the higher the order of the polynomial) unless there are sufficient ultrasonic paths and transducers provided in the ultrasonic flow meter to improve the accuracy. For example, if a 15th order polynomial is required to describe the flow profile, a total of eight ((15+1)/2) ultrasonic paths and transducer pairs can be required for improved accuracy, which would also increase the cost of the flow meter, increase the difficulty of installation of large number of transducers given the space constraints on a pipe, and increase the disturbances to the flow causes by the large number of transducer ports. Furthermore, since several ultrasonic flow meter applications require highly accurate flow rate (Vm) measurements (e.g., custody-transfer applications for liquid hydrocarbons often require an accuracy of 0.15%) that cannot always be achieved by simply increasing the number of ultrasonic paths and transducers in these prior art solutions, complicated calibration schemes using calibration factors that are based on several of the physical properties of the fluid and/or the pipe (e.g., the Reynolds number) are required anyway to achieve the required accuracy, which increases the costs and complexity of the flow meters.
It would be advantageous to provide an ultrasonic flow meter that can provide highly accurate flow rate measurements without the use of complicated flow meter design and calibration schemes.
BRIEF DESCRIPTION OF THE INVENTION
In one embodiment, a multi-path ultrasonic flow meter for determining the flow rate of a fluid in a conduit is disclosed comprising at least two transducer pairs attached to the conduit at two chord locations, one greater than and one less than a mid-radius chord, wherein the composite ratio the two path velocities to the flow rate is substantially constant over the range of Reynolds numbers.
In another embodiment, a method of determining the flow rate of a fluid in a conduit is disclosed comprising the steps of determining a composite velocity by determining a weighted average of a plurality of path velocities, determining a chord velocity ratio based on the path velocities, determining a profile correction factor based on the composite velocity and the chord velocity ratio, and determining the flow rate based on the composite velocity and the profile correction factor.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a cross-section of an ultrasonic flow meter cut along the diameter plane of the conduit in one embodiment of the invention.
FIG. 2 shows the ultrasonic path of FIG. 1 viewed along section A-A of FIG. 1.
FIG. 3 is a graph showing the variance of the ratio (Vpi/Vm) over a range of Reynolds numbers (Re) for various chord locations.
FIG. 4 is a graph showing the variance of the composite ratio ((Vp1+Vp2)/Vm) over a range of Reynolds numbers (Re) for various pairs of chord locations.
FIG. 5 is a graph showing the variance of the composite ratio ((Vp1+Vp2)/Vm) over a range of Reynolds numbers for an exemplary pair of chord locations for various relative pipe wall roughnesses.
FIG. 6 is a perspective view of a multi-path ultrasonic flow meter in one embodiment of the invention.
FIG. 7 is an end view of the multi-path ultrasonic flow meter in FIG. 6.



DETAILED DESCRIPTION OF THE INVENTION
FIG. 1 illustrates an ultrasonic flow meter 1 employing transit time flow metering to determine the flow rate of a fluid. Although FIG. 1 shows a single pair of transducers 40, 50 attached to a conduit 20 forming a single ultrasonic path 30, it will be understood by a person of ordinary skill in the art that two or more pairs of transducers 40, 50, each forming separate ultrasonic paths, can be used to form a multi-path ultrasonic flow meter. Based on the fluid flow direction 22, one transducer 40 can be installed upstream of the downstream transducer 50 on the conduit 20, which can be a section of the pipe or a separate spool piece hot-tapped to the pipeline. The ultrasonic path 30 can be on the center (i.e., diameter (D)) of the conduit 20 or on chordal paths (i.e., path not on the center of the conduit 20). Each transducer 40, 50, when energized, transmits an ultrasonic signal along the ultrasonic path 30 through the flowing fluid that is received by and detected by the other transducer 50, 40.
The path velocity (Vp) of the fluid averaged along an ultrasonic path 30 can be determined as a function of the differential between the upstream transit time (tup) (i.e., the time of an ultrasonic signal traveling along the ultrasonic path 30 from the downstream transducer 50 upstream to the upstream transducer 40 against the fluid flow direction 22) and the downstream transit time (tdn) (i.e., the time of an ultrasonic signal traveling along the ultrasonic path 30 from the upstream transducer 40 downstream to the downstream transducer 50 with the fluid flow direction 22) of the ultrasonic signals. In the presence of fluid flow, the downstream transit time (tdn) is faster (or shorter) that the upstream transit time (tup). Since the transit time differential (Δt) is proportional to the path velocity (Vp) of the fluid, the path velocity (Vp) of the fluid averaged along the ultrasonic path 30 can be determined as a function of the transit time differential (Δt) as well as other known parameters, such as the conduit 20 diameter (D), ultrasonic path length (P), and the angle (θ) formed between the ultrasonic path 30 and the fluid flow direction 22 as shown in FIG. 1:




  
    
      
        
          V
          p
        
        =
        
          
            P
            2
          
          ⁢
          
            (
            
              
                
                  t
                  up
                
                -
                
                  t
                  dn
                
              
              
                
                  (
                  
                    
                      t
                      dn
                    
                    ×
                    
                      t
                      up
                    
                  
                  )
                
                ⁢
                cos
                ⁢
                
                  
                
                ⁢
                θ
              
            
            )
          
        
      
    
    
      
        (
        1
        )
      
    
  
  
    
      
        
          cos
          ⁢
          
            
          
          ⁢
          θ
        
        =
        
          L
          P
        
      
    
    
      
        (
        2
        )
      
    
  
  
    
      
        
          V
          p
        
        =
        
          
            
              P
              2
            
            
              2
              ⁢
              L
            
          
          ⁢
          
            (
            
              
                
                  t
                  up
                
                -
                
                  t
                  dn
                
              
              
                
                  t
                  dn
                
                ×
                
                  t
                  up
                
              
            
            )
          
        
      
    
    
      
        (
        3
        )
      
    
  




FIGS. 1 and 2 illustrate the ultrasonic flow meter 1 of FIG. 1 along the ultrasonic path 30, demonstrating the geometric relationship between the ultrasonic path 30 established by placing transducers 40, 50 at the xi chord location on the conduit 20. The path velocity (Vpi) averaged along a particular ultrasonic path 30 corresponding to xi chord location can be expressed by:




  
    
      
        
          V
          pi
        
        =
        
          
            ∫
            
              y
              ⁡
              
                (
                min
                )
              
            
            
              y
              ⁡
              
                (
                max
                )
              
            
          
          ⁢
          
            V
            ⁢
            
              
            
            ⁢
            
              point
              ⁢
              
                
              
              (
              
                
                  x
                  i
                
                ,
                y
              
              )
            
            ⁢
            
              d
              l
            
          
        
      
    
    
      
        (
        4
        )
      
    
  




where:

Vpoint is the velocity at a particular point along the ultrasonic path 30;

y(min)=−√{square root over (R2−xi2)};  (5)

y(max)=√{square root over (R2−xi2)}; and  (6)

R=radius of the conduit (20)
Furthermore, it is known that the-path velocity (Vpi) averaged along a particular ultrasonic path 30 corresponding to xi chord location for a well-developed turbulent fluid flow (e.g., fluid flow through sufficiently long straight lengths of pipe with minimal disturbances before and after the ultrasonic flow meter 1) can be expressed by:




  
    
      
        
          V
          pi
        
        =
        
          
            ∫
            
              -
              
                
                  
                    R
                    2
                  
                  -
                  
                    x
                    i
                    2
                  
                
              
            
            
              
                
                  R
                  2
                
                -
                
                  x
                  i
                  2
                
              
            
          
          ⁢
          
            
              
                V
                m
              
              [
              
                
              
              ⁢
              
                1
                -
                
                  
                    0.715
                    +
                    
                      1.0751
                      *
                      
                        log
                        (
                        
                          
                            
                              
                                x
                                i
                                2
                              
                              +
                              
                                y
                                2
                              
                            
                          
                          R
                        
                        )
                      
                    
                  
                  
                    log
                    ⁡
                    
                      (
                      
                        
                          
                            0.2703
                            ⁢
                            
                              
                            
                            ⁢
                            ɛ
                          
                          D
                        
                        +
                        
                          5.74
                          
                            R
                            e
                            0.9
                          
                        
                      
                      )
                    
                  
                
              
              ]
            
            ⁢
            
              d
              y
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




where:

Vm is the flow rate (i.e., the mean flow rate);
ε/D is the relative pipe wall roughness;
Re is the Reynolds number; and
R=radius of the conduit 20.
This equation can be rewritten as a ratio of the-path velocity (Vpi) averaged along a particular ultrasonic path 30 to the flow rate (Vm) for xi chord location as follows:




  
    
      
        
          
            V
            pi
          
          
            V
            m
          
        
        =
        
          
            ∫
            
              -
              
                
                  
                    R
                    2
                  
                  -
                  
                    x
                    i
                    2
                  
                
              
            
            
              
                
                  R
                  2
                
                -
                
                  x
                  i
                  2
                
              
            
          
          ⁢
          
            
              [
              
                
              
              ⁢
              
                1
                -
                
                  
                    0.715
                    +
                    
                      1.0751
                      *
                      
                        log
                        (
                        
                          
                            
                              
                                x
                                i
                                2
                              
                              +
                              
                                y
                                2
                              
                            
                          
                          R
                        
                        )
                      
                    
                  
                  
                    log
                    ⁡
                    
                      (
                      
                        
                          
                            0.2703
                            ⁢
                            
                              
                            
                            ⁢
                            ɛ
                          
                          D
                        
                        +
                        
                          5.74
                          
                            R
                            e
                            0.9
                          
                        
                      
                      )
                    
                  
                
              
              ]
            
            ⁢
            
              d
              y
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




This equation can be used to determine the variance of the ratio of the path velocity (Vpi) to the flow rate (Vm) for various xi chord locations (i.e., location of a single pair of transducers 40, 50) over a range of Reynolds numbers (Re). For example, assuming a nominal radius of 1 (R=1) and a relative pipe wall roughness of 10−5(ε/D=10−5), FIG. 3 shows the variance of the ratio (Vpi/Vm) over a range of Reynolds numbers (Re) from 2,000 to 2,000,000 for xi chord locations of 0.00, 0.10, 0.20, 0.30, 0.34, 0.40, 0.50, 0.55, 0.60, 0.70, and 0.80. The results shown in FIG. 3 demonstrate the at a chord location of xi=0.55, the ratio (Vpi/Vm) is substantially constant over the range of Reynolds numbers (Re), including all sub-ranges of Reynolds numbers (Re) (e.g., 2,000 to 500,000). Accordingly, if a single pair of transducers 40, 50 were installed at chord location of xi=0.55, knowledge of the Reynolds number (Re) or kinematic viscosity (v) of the fluid would not necessarily be required to determine the flow rate (Vm) based on the path velocity (Vpi) averaged along the ultrasonic path 30 at that chord location. As used herein, the mid-radius chord (xmr) is the chord location xi where the ratio (Vpi/Vm) is substantially constant over the range of Reynolds numbers (Re). The use of a single pair of transducers 40, 50 and a single ultrasonic path 30 at the mid-radius chord xmr is not preferred as it does not offer redundancy in the event of failure of one or more transducers 40, 50.
The results shown in FIG. 3 can also be used to identify multiple chord locations xi that can be used in conjunction with each other where knowledge of the Reynolds number (Re) or kinematic viscosity (v) of the fluid would not necessarily be required to determine the flow rate (Vm) based on the path velocity (Vpi) averaged along the ultrasonic paths 30 at those chord locations. For example, the results shown in FIG. 3 demonstrate that for chord locations xi greater than the mid-radius chord (e.g., xi>xmr), the ratio of the path velocity (Vpi) to the flow rate (Vm) increases as the Reynolds number (Re) increases, while for chord locations less than the mid-radius chord (xi<xmr), the ratio of the path velocity (Vpi) to the flow rate (Vm) decreases as the Reynolds number (Re) increases.
Based on these relationships, the inventors have discovered that one or more chord locations xi greater than the mid-radius chord xmr can be used along with one or more chord locations xi below the mid-radius chord xmr to provide a composite ratio of the path velocity to the pipe flow rate for these multiple (n) paths (e.g., (Vp1/Vm+Vp2/Vm+ . . . Vpn/Vm) or (Vp1+Vp2+Vpn)/Vm) that is substantially constant over the range of Reynolds numbers (Re). Effectively, the effect of the Reynolds number (Re) on the path velocities (Vpi) for the chord locations greater than the mid-radius chord xmr are cancelled by the effect of the Reynolds number (Re) on the-path velocities (Vpi) for the chord locations less than the mid-radius chord xmr. Accordingly, if pairs of transducers 40, 50 are installed at properly selected multiple chord locations xi, knowledge of the Reynolds number (Re) or kinematic viscosity (v) of the fluid would not necessarily be required to determine the flow rate (Vm) based on the path velocities (Vpi) averaged along the ultrasonic paths 30 at the chord locations. FIG. 4 shows the variance of the composite ratio ((Vp1+Vp2)/Vm) over a range of Reynolds numbers (Re) from 2,000 to 2,000,000 for three exemplary pairs of chord locations (x1=0.75 and x2=0.18; x1=0.70 and x2=0.34; and x1=0.65 and x2=0.45). The results shown in FIG. 4 demonstrate that for each of these pairs of chord locations (x1 and x2), the composite ratio ((Vp1+Vp2)/Vm) is substantially constant over the range of Reynolds numbers (Re).
A person of ordinary skill in the art would understand that several other combinations of chord locations, some greater than and some less than the mid-radius chord xmr can be selected beyond these examples to produce a composite ratio ((Vp1+Vp2+Vpn)/Vm) that is substantially constant over the range of Reynolds numbers (Re). For example, while each of the exemplary pairs of chord locations shown in FIG. 4 include a single chord location less than the mid-radius chord xmr(e.g., x1=0.70) and a single chord location greater than the mid-radius chord (x2=0.34), other combinations could include multiple chord locations greater than and less than the mid-radius chord xmr. Similarly, other combinations could include a different number of chord locations (e.g., two) greater than the mid-radius chord xmr than the number of chord locations (e.g., one) less than the mid-radius chord xmr as long as the composite ratio of the path velocity to the pipe flow rate for these multiple (n) paths is substantially constant over the range of Reynolds numbers (Re).
The inventors have also discovered that the composite ratios (e.g., (Vp1+Vp2)/Vm) are also substantially constant over the range of Reynolds numbers (Re) even when the relative pipe wall roughness (ε/D) changes over time. Based on equation (8), FIG. 5 shows the variance of the composite ratio ((Vp1+Vp2)/Vm) over a range of Reynolds numbers (Re) from 2,000 to 2,000,000 for chord locations (x1=0.70 and x2=0.34) for three exemplary relative pipe wall roughness (ε/D=10−3, 10−5, and 10−8). FIG. 5 demonstrates that even if the relative pipe wall roughness (ε/D) for a particular pipe were to significantly change (i.e., from 10−3 to 10−8) after calibration of an ultrasonic flow meter 1 using the exemplary chord locations, that change would not necessarily require additional calibration since the change would have a negligible impact on the accuracy of the ultrasonic flow meter 1.
In one embodiment of a multi-path ultrasonic flow meter 10 shown in FIGS. 6 and 7, a total of four ultrasonic paths (n=4) (p1, p2, p3, and p4) 31, 32, 33, 34 are formed by four pairs of upstream transducers 41, 42, 43, 44 and downstream transducers 51, 52, 53, 54 at the following exemplary chord locations: xi=0.70, x2=0.34 x3=0.34, x4=−0.70) on the conduit 20, which is in the form of a spool piece. These chord locations xi are selected on both sides of the conduit (one side 0.00<xi<1.00 and the other side 0.00<xi<−1.00) each side having a mid-radius chord (e.g., xmr12=0.55 and xmr34=−0.55), such that each side has at least one chord location xi greater than the mid-radius chord and at least one chord location less than the mid-radius chord. As discussed previously, each set of two pairs of chord locations (i.e., x1 and x2; x3 and x4) can provide a composite ratio of the path velocity to the pipe flow rate for each pair of ultrasonic paths ((Vp1+Vp2)/Vm) that is substantially constant over the range of Reynolds numbers (Re). Furthermore, the entire set of four pairs of chord locations (x1, x2, x3, and x4) can also provide a composite ratio of the path velocity to the pipe flow rate for all ultrasonic paths ((Vp1+Vp2+Vp3+Vp4)/Vm) that is substantially constant over the range of Reynolds numbers (Re).
In order to provide a well-developed turbulent fluid flow for the ultrasonic flow meter 10, the conduit 20 should be installed such that there are sufficiently long straight lengths of pipe with minimal disturbances before and after the ultrasonic flow meter 10. Each of the transducers are connected to the flow meter electronics 60, which can include a microprocessor for determining the path velocities and flow rate, an electronic controller to transmit and receive signals from the transducers, a display for displaying the flow rate, and an I/O port for external communications.
In order to determine the flow rate using the multi-path ultrasonic flow meter 10, measurements can be performed along each of the multiple (n) ultrasonic paths to determine the path velocities (Vp1, Vp2, Vp3, Vp4, Vpn). Weights (wi) can be assigned each-path velocities (Vpi) to determine a composite velocity (Vc) by calculating a weighted average:




  
    
      
        Vc
        =
        
          
            ∑
            
              i
              =
              1
            
            n
          
          ⁢
          
            
              w
              i
            
            ⁢
            
              V
              pi
            
          
        
      
    
    
      
        (
        9
        )
      
    
  




In the embodiment using the four exemplary chord locations (i.e., x1=0.70, x2=0.34 x3=−0.34, x4=−0.70), equal weights can be used (w1=w2=w3=w4=(1/n)=1⁄4=0.25) in determining the composite velocity (Vc):




  
    
      
        Vc
        =
        
          
            ∑
            
              i
              =
              1
            
            4
          
          ⁢
          
            
              (
              0.25
              )
            
            ⁢
            
              V
              pi
            
          
        
      
    
    
      
        (
        10
        )
      
    
  




As evidenced by the equations above, by using the exemplary chord locations, this composite velocity (Vc) can be determined without any knowledge of the Reynolds number (Re) or several of the other factors affecting a flow profile (e.g., other fluid characteristics, pipe roughness, pipe configuration (e.g., straight, bends, elbows, tees, valves, headers, etc.), and disturbances caused by the ultrasonic flow meter transducers). Despite not taking into account the flow profile, this composite velocity (Vc) has proven to be highly accurate, with errors between the measured composite velocity (Vc) and the known and actual reference flow rate (Vr) used for calibration, generally less than one percent for a wide range of Reynolds numbers (Re) as confirmed in laboratory testing. In order to obtain a more accurate flow rate (Vm) beyond the accuracy provided by use of only the composite velocity (Vc), however, a simply and empirically determined profile correction factor (K) can be used to account for the effect of the unknown flow profile.
In one embodiment for calibrating the multi-path ultrasonic flow meter 10, an average measured composite velocity (Vc(avg)) and average actual reference flow rate (Vr(avg)) is determined for a reference flow rate (e.g., Vr=1,000 gallons/min.) over a specified period of time (e.g., 1 minute). The profile correction factor (K) for this reference flow rate (V) can be determined by:




  
    
      
        K
        =
        
          
            V
            
              r
              ⁡
              
                (
                avg
                )
              
            
          
          
            V
            
              c
              ⁡
              
                (
                avg
                )
              
            
          
        
      
    
    
      
        (
        11
        )
      
    
  




This profile correction factor (K) provides the actual error percentage of the composite velocity (Vc) at the reference flow rate (Vr). In addition, in one embodiment for calibrating the multi-path ultrasonic flow meter 10, a chord velocity ratio (Rcv) can also be determined. In the embodiment using the four exemplary chord locations (x1=0.70, x2=0.34 x3=−0.34, x4=−0.70), the chord velocity ratio (Rcv) can be determined by:




  
    
      
        
          R
          cv
        
        =
        
          
            
              V
              
                p
                ⁢
                
                  
                
                ⁢
                1
              
            
            +
            
              V
              
                p
                ⁢
                
                  
                
                ⁢
                4
              
            
          
          
            
              V
              
                p
                ⁢
                
                  
                
                ⁢
                2
              
            
            +
            
              V
              
                p
                ⁢
                
                  
                
                ⁢
                3
              
            
          
        
      
    
    
      
        (
        12
        )
      
    
  




This process can be repeated during calibration for several different reference flow rates (Vr), each time associating a profile correction factor (K) with a composite velocity (Vc) and chord velocity ratio (Rcr). Once a sufficient number of reference flow rates (Vr) have been used during calibration, each providing a profile calibration factor (K) associated with a composite velocity (Vc) and a chord velocity ratio (Rcr), the profile correction factor (K) can be represented as a function of the composite velocity (Vc) and the chord velocity ratio (Rcv) (e.g., K can be looked up in a table when Vc and Rcv are known).

K=f(Vc,Rcv)  (13)

In order to minimize the required number of reference flow rates (Vr) in calibration, interpolation can be used to provide the profile correction factors (K) for combinations of composite velocities (Vc) and chord velocity ratios (Rcv) not encountered during calibration. When this interpolation is complete, this profile correction factor (K) table or function can be stored in the memory of the flow meter electronics 60. Once calibrated, when the multi-path ultrasonic flow meter 10 measures a composite velocity (Vc) and a chord velocity ratio (Rcv), the profile correction factor (K) can be determined based on the stored stable or function and applied to the composite velocity (Vc) to provide a more accurate flow rate (Vm) to be displayed:

Vm=K(Vc)  (14)

As evidenced by the equations above, this simplified calibration can be performed and the flow rate (Vm) can be determined without any knowledge of the Reynolds number (Re) or several of the other factors affecting a flow profile (e.g., other fluid characteristics, pipe roughness, pipe configuration (e.g., straight, bends, elbows, tees, valves, headers, etc.), and disturbances caused by the ultrasonic flow meter transducers). Despite not taking into account the flow profile, this flow rate (Vm) has proven to be highly accurate, with errors between the calibration corrected flow rate (Vm) and the known and actual reference flow rate (Vr) used for calibration, generally less than 0.15 percent over the range of Reynolds numbers (Re) as confirmed in laboratory testing. Furthermore, as discussed above and shown in FIG. 5, even if the relative pipe wall roughness (ε/D) for a particular pipe were to significantly change (i.e., from 10−3 to 10−8) after calibration of an ultrasonic flow meter 10 using the exemplary chord locations, that change would not necessarily require additional calibration since the change would have a negligible impact on the accuracy of the ultrasonic flow meter 10.
This written description uses examples to disclose the invention, including the best mode, and also to enable any person skilled in the art to make and use the invention. The patentable scope of the invention is defined by the claims, and may include other examples that occur to those skilled in the art. Such other examples are intended to be within the scope of the claims if they have structural elements that do not differ from the literal language of the claims, or if they include equivalent structural elements with insubstantial differences from the literal language of the claims.

## Claims

1. A multi-path ultrasonic flow meter for determining the flow rate of a fluid in a conduit comprising:
said conduit having a mid-radius chord location where the ratio of the path velocity of said fluid averaged along said mid-radius chord to said flow rate is substantially constant over a range of Reynolds numbers;
a first transducer pair attached to said conduit forming a first ultrasonic path at a first chord location at a first distance from the center of said conduit, wherein said first distance is greater than the distance of said mid-radius chord to said center;
a second transducer pair attached to said conduit forming a second ultrasonic path at a second chord location at a second distance from said center of said conduit, wherein said second distance is less than said distance of said mid-radius chord to said center;
flow electronics for determining a first path velocity of said fluid averaged along said first ultrasonic path and a second path velocity of said fluid averaged along said second ultrasonic path, wherein the composite ratio of said first path velocity and said second path velocity to said flow rate is substantially constant over said range of Reynolds numbers.

2. The multi-path ultrasonic flow meter of claim 1, further comprising:
a third transducer pair attached to said conduit forming a third ultrasonic path at a third chord location at a third distance from said center of said conduit, wherein said third distance is greater than said distance of said mid-radius chord to said center;
said flow electronics for determining a third path velocity of said fluid averaged along said third ultrasonic path, wherein the composite ratio of said first path velocity, said second path velocity, and said third path velocity to said flow rate is substantially constant over said range of Reynolds numbers.

3. The multi-path ultrasonic flow meter of claim 1, further comprising:
a third transducer pair attached to said conduit forming a third ultrasonic path at a third chord location at a third distance from said center of said conduit, wherein said third distance is less than said distance of said mid-radius chord to said center;
said flow electronics for determining a third path velocity of said fluid averaged along said third ultrasonic path, wherein the composite ratio of said first path velocity, said second path velocity, and said third path velocity to said flow rate is substantially constant over said range of Reynolds numbers.

4. The multi-path ultrasonic flow meter of claim 1, wherein said range of Reynolds numbers is from 2,000 to 2,000,000.

5. The multi-path ultrasonic flow meter of claim 1, wherein said range of Reynolds numbers is from 2,000 to 500,000.

6. A method of determining the chord locations for a multi-path ultrasonic flow meter for determining the flow rate of a fluid in a conduit comprising the steps of:
determining a mid-radius chord location in said conduit where the ratio of the path velocity of said fluid averaged along said mid-radius chord to said flow rate is substantially constant over a range of Reynolds numbers;
determining a first chord location of a first transducer pair attached to said conduit forming a first ultrasonic path at a first distance from the center of said conduit, wherein said first distance is greater than the distance of said mid-radius chord to said center;
determining a second chord location of a second transducer pair attached to said conduit forming a second ultrasonic path at a second distance from said center of said conduit, wherein said second distance is less than said distance of said mid-radius chord to said center;
determining a first path velocity of said fluid averaged along said first ultrasonic path and a second path velocity of said fluid averaged along said second ultrasonic path, wherein said first chord location and said second chord location are selected such that the composite ratio of said first path velocity and said second path velocity to said flow rate is substantially constant over said range of Reynolds numbers.

7. The method of claim 6, further comprising the steps of:
determining a third chord location of a third transducer pair attached to said conduit forming a third ultrasonic path at a third distance from said center of said conduit, wherein said third distance is greater than said distance of said mid-radius chord to said center;
determining a third path velocity of said fluid averaged along said third ultrasonic path, wherein said third chord location is selected such that the composite ratio of said first path velocity, said second path velocity, and said third path velocity to said flow rate is substantially constant over said range of Reynolds numbers.

8. The method of claim 6, further comprising the steps of:
determining a third chord location of a third transducer pair attached to said conduit forming a third ultrasonic path at a third distance from said center of said conduit, wherein said third distance is less than said distance of said mid-radius chord to said center;
determining a third path velocity of said fluid averaged along said third ultrasonic path, wherein said third chord location is selected such that the composite ratio of said first path velocity, said second path velocity, and said third path velocity to said flow rate is substantially constant over said range of Reynolds numbers.

9. The method of claim 6, wherein said range of Reynolds numbers is from 2,000 to 2,000,000.

10. The method of claim 6, wherein said range of Reynolds numbers is from 2,000 to 500,000.

11. A method of determining the flow rate of a fluid in a conduit comprising the steps of:
determining the path velocities of said fluid averaged along a plurality of ultrasonic paths formed by a plurality of transducer pairs at a plurality of chord locations;
determining a composite velocity by determining a weighted average of said path velocities;
determining a chord velocity ratio based on said path velocities;
determining a profile correction factor based on said composite velocity and said chord velocity ratio; and
determining said flow rate based on said composite velocity and said profile correction factor.

12. The method of claim 11, wherein said step of determining said composite velocity uses equal weights for said path velocities.

13. The method of claim 11, wherein said step of determining a profile correction factor based on said composite velocity and said chord velocity ratio comprises the step of retrieving said profile correction factor in a table.

14. The method of claim 11, wherein said step of determining said flow rate based on said composite velocity and said profile correction factor, includes multiplying said composite velocity by said profile correction factor.