Ultrasonic transducer system and method using broadband system responses

Abstract

A transducer system with a transducer and circuitry for applying a waveform to excite the transducer during an excitation period. The applying circuitry also comprises circuitry for changing a frequency of the waveform during the excitation period.

Claims

1. A transducer system, comprising: an ultrasonic transducer; and processing circuitry configured to apply a waveform to excite the ultrasonic transducer during an excitation period, the applied waveform having a changing ultrasonic frequency during the excitation period; wherein the applied waveform is generated by the processing circuitry as a pulse train including a first plurality of sequential sets of pulses, wherein each of the first plurality of sequential sets of pulses includes two or more pulses, wherein each of the sets of pulses of the first plurality of sequential sets of pulses has a different respective ultrasonic frequency, and wherein all pulses of each set of pulses of the first plurality of sequential sets of pulses have the same respective ultrasonic frequency and the same amplitude.

2. The transducer system of claim 1, wherein the sets of pulses provide a chirp signal in the pulse train.

3. The transducer system of claim 1, wherein the processing circuitry is further configured to receive an echo waveform in response to excitation of the ultrasonic transducer, and to determine an absolute time of flight between the pulse train and the echo waveform.

4. The transducer system of claim 3, wherein the processing circuitry is further configured to measure a plurality of correlation measures between the pulse train and the echo waveform.

5. The transducer system of claim 3, wherein the processing circuitry is further configured to measure a first plurality of correlation measures between a first instance of a pulse train and a respective first echo waveform, and to measure a second plurality of correlation measures between a second instance of a pulse train and a respective second echo waveform, wherein the second plurality of correlation measures is lesser than the first plurality of correlation measures.

6. The transducer system of claim 1: wherein the ultrasonic transducer is a first ultrasonic transducer of the transducer system and produces a first transmitted signal in response to a first pulse train, wherein the pulse train including the first plurality of sequential sets of pulses is the first pulse train; wherein the transducer system further comprises: a second ultrasonic transducer for producing a first received echo waveform in response to the transmitted signal; and wherein the processing circuitry is further configured to determine an absolute time of flight between the first pulse train and the first received echo waveform.

7. The transducer system of claim 6, wherein the processing circuitry is further configured to apply a second pulse train including a second plurality of sequential sets of pulses to excite the second ultrasonic transducer, wherein the second ultrasonic transducer produces a second transmitted signal in response to the second pulse train, and wherein the second pulse train has a changing ultrasonic frequency such that each of the sets of pulses of the second plurality of sequential sets of pulses has a different respective ultrasonic frequency; and wherein the first ultrasonic transducer is configured to produce a second received echo waveform in response to the second transmitted signal.

8. The transducer system of claim 7, wherein the processing circuitry is further configured to determine an absolute time of flight between the second pulse train and the second received echo waveform.

9. The transducer system of claim 8, wherein the processing circuitry is further configured to determine a relative time of flight between the first received echo waveform and the second received echo waveform.

10. The transducer system of claim 6, wherein the processing circuitry is further configured to determine a velocity of flow in response to time of flight between the first pulse train and the first received echo waveform.

11. The transducer system of claim 1, wherein, for each set of pulses of the first plurality of sequential sets of pulses other than a sequentially first set of pulses and a sequentially last set of pulses, the respective ultrasonic frequency of the set of pulses differs from the respective ultrasonic frequency of the immediately preceding set of pulses and the respective ultrasonic frequency of the immediately succeeding set of pulses by a same amount.

12. The transducer system of claim 11, wherein the same amount is determined based on dividing a desired frequency band by one less than a total number of sets of pulses in the first plurality of sequential set of pulses.

13. The transducer system of claim 11, wherein, for each set of pulses of the first plurality of sequential sets of pulses other than the sequentially first set of pulses, the respective ultrasonic frequency of the set of pulses is greater than the respective ultrasonic frequency of the immediately preceding set of pulses.

Description

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING

(1) FIG. 1 illustrates a diagram of a flow meter transducer system.

(2) FIG. 2 illustrates a sequence diagram of the excitation and response waveforms of the flow meter transducer system of FIG. 1.

(3) FIG. 3 illustrates a flowchart of a method of operation of the processor of FIG. 1.

(4) FIG. 4 depicts a sequence of pulses with changing (e.g., increasing) frequencies according to a preferred embodiment.

(5) FIGS. 5A through 5G illustrate a digital sample set of the transmitted pulse train PT and a digital sample set of the received echo waveform EVM, with different time shifted correlation measures taken for the differing time-shifted positions of the sets.

(6) FIG. 6 illustrates a generally sinusoidal curve SC as a depiction of an analog correlation between the UPS and DNS waveforms.

(7) FIG. 7 illustrates a partial plot of three cycles from a correlation of UPS and DNS waveforms, where the illustrated partial plots show a main lobe ML, a left side lobe LSL and a right side lobe RSL.

(8) FIG. 8 illustrates four frequency response plots for a first flowmeter system such as system 10 in FIG. 1, where each plot corresponds to a different temperature.

(9) FIG. 9 illustrates four frequency response plots for a second flowmeter system such as system 10 in FIG. 1, where each plot corresponds to a different temperature.

DETAILED DESCRIPTION OF EMBODIMENTS

(10) FIG. 1 illustrates a diagram of a flow meter transducer system 10. System 10 as described in this initial paragraph is known in the art, but it also is included and improved upon in connection with the preferred embodiments as detailed hereafter. System 10 includes a pipe 12 through which a material, such as water or gas, may flow, and for sake of reference in FIG. 1 the flow direction is shown from right to left. Within pipe 12 are located a first transducer TR.sub.1 and a second transducer TR.sub.2. Proximate transducer TR.sub.1 is a reflector R.sub.1, and proximate transducer TR.sub.2 is a reflector R.sub.2. Each reflector R.sub.x is oriented so as to communicate a signal with its respective proximate transducer TR.sub.x and also to communicate a signal to the other reflector. As shown by a dashed line inside pipe 12 in FIG. 1, therefore, transducer TR.sub.1 communicates a signal to reflector R.sub.1, which reflects that signal to reflector R.sub.2, which reflects that signal to transducer TR.sub.2. Likewise, transducer TR.sub.2 communicates a signal to reflector R.sub.2, which reflects that signal to reflector R.sub.1, which reflects that signal to transducer TR.sub.1. Lastly, a processor 14, such as a digital signal processor, microprocessor, microcontroller, or some other electronic circuitry, receives a clock signal from a clock 16, and processor 14 is coupled to both transducer TR.sub.1 and transducer TR.sub.2 for exciting either transducer TR.sub.x to transmit a signal and to process a correspondingly received signal by the other transducer, as further explored below. Clock 16 is typically a low-power (e.g., power consumption ˜140 μA/MHz) crystal oscillator with speeds in the range of 1-24 MHz.

(11) FIG. 2 illustrates a sequence diagram of the flow meter transducer system 10 of FIG. 1, where the sequence as illustrated and first described is also known in the art, but it also is included and improved upon in connection with the preferred embodiments as detailed hereafter. In general, FIG. 2 represents a first transmission in a first direction from first transducer TR.sub.1 to second transducer TR.sub.2, followed by a second transmission in second direction, opposite the first, from second transducer TR.sub.2 to first transducer TR.sub.1. For sake of reference, the first direction is shown and referred to as downstream (DNS) in the context of pipe 12 in FIG. 1, and the second direction is shown and referred to as upstream (UPS) in the context of pipe 12 in FIG. 1.

(12) Looking at the DNS timeline, at a time t.sub.0, processor 14, either directly or via additional circuitry (not shown), applies an excitation pulse train to transducer TR.sub.1 which in response transmits a corresponding ultrasonic pulse train signal, where both the applied and transmitted pulse trains consist of a number of pulses, such as 10 to 40 such pulses, shown beginning therefore at time t.sub.0 and ending at time t.sub.1. The transmitted pulses are directed to reflector R.sub.1, which reflects them to reflector R.sub.2, toward transducer TR.sub.2. Time passes as these signals pass along a channel that includes the portion of the interior of pipe 12, and any material inside that pipe, that is generally along the dashed line shown in FIG. 1. This time is referred to in the art as time of flight (TOF). Thus, the DNS TOF occurs between times t.sub.0 and t.sub.2, or it may be measured relative to some known offset from either of those times.

(13) At time t.sub.2, second transducer TR.sub.2 begins to respond to the pulses transmitted by first transducer TR.sub.1, as shown by the received signal that begins at time t.sub.2. Processor 14, being coupled to second transducer TR.sub.2, is operable to sample this received signal. Toward this end, processor 14 is typically programmed to begin capturing signals from second transducer TR.sub.2 at a time slightly before the expected response at time t.sub.2. While the received signal is analog, the captured samples are typically in digital form, so included in the sampling is an analog-to-digital (ADC) conversion, either as part of processor 14 or an additional intermediate device (not shown). Moreover, the sampling rate is commonly, but not necessarily, a multiple of (e.g., four or five times) of the pulse frequency f.sub.PLS or of the entire pulse train frequency or of the nominal resonance frequency. For each sample, a sample amplitude S.sub.x is stored, and also preferably stored is a corresponding sample time st.sub.x of when the sample was taken. Note that sample time st.sub.x can be an actual time or based on a count of either time or a number of samples (which thereby can represent time based on the sample period). For sake of reference, therefore, each sample may be represented by the data pair of (S.sub.x, st.sub.x). In any event, ideally, the received signal at time t.sub.2 will generally increase toward a peak amplitude at time t.sub.3 and decay thereafter, as the transmitted pulses create an oscillatory response in, and due to the resistive/inductive/capacitive (RLC) nature of, the receiving transducer TR.sub.2.

(14) At time t.sub.5, which is preferably some time after the amplitude of the received signal at transducer TR.sub.2 has sufficiently declined, processor 14 reverses the communication direction of system 10, such as via a multiplexer (MUX) operation. Thus, at a later time t.sub.5, the above process repeats, but in the reverse (i.e., UPS) direction, that is, from second transducer TR.sub.2 to first transducer TR.sub.1. From the preceding discussion, therefore, one skilled in the art will appreciate that at time t.sub.5, processor 14 applies an excitation pulse train at frequency f.sub.PLS to second transducer TR.sub.2, causing it to commence emission of a corresponding pulse train, where both the excitation and resultant transmission pulse train consist of the same number of pulses (e.g., 10 to 40) as for the DNS pulses and hence through time t.sub.6, directed to reflector R.sub.2, which reflects them to reflector R.sub.1, toward transducer TR.sub.1. Following the TOF of these pulses, at time t.sub.7, first transducer TR.sub.1 begins to respond to the pulses transmitted by first transducer TR.sub.2, with the transducer TR.sub.1 received signal again detected, converted to digital, and sampled by processor 14, increasing toward a peak amplitude at time t.sub.8, and decaying thereafter to a zero or near-zero level by time t.sub.9.

(15) Given the timing of FIG. 2, processor 14 is operable (e.g., via appropriate programming) to determine the UPS TOF, the DNS TOF, and the relative difference of the UPS and DNS TOF, using preferred embodiment methodologies detailed later. From these measures, the flow rate through pipe 12 may be calculated, for example, based on the following Equation 1:

(16) $\begin{array}{cc}v=\frac{L}{2}\times \left(\frac{1}{{\mathrm{TR}}_{12}}-\frac{1}{{\mathrm{TR}}_{21}}\right)=\frac{L}{2}\times \left(\frac{{\mathrm{TR}}_{21}-{\mathrm{TR}}_{12}}{{\mathrm{TR}}_{21}{\mathrm{TR}}_{12}}\right)& \mathrm{Equation}1\end{array}$
where,

(17) L is the length of the channel path between first transducer TR.sub.1 and second transducer TR.sub.2;

(18) TR.sub.12 is the DNS TOF;

(19) TR.sub.21 is the UPS TOF; and

(20) v is the velocity of flow.

(21) Note from Equation 1, therefore, that flow velocity is directly related to the DNS and UPS TOF, and the difference between them. Thus, processor 14 may measure a set of UPS and DNS TOF, or multiple sets, and accordingly determine flow velocity based on those measures. Moreover, the accuracy of the TOF measures directly influences the accuracy of the velocity determination.

(22) As would be expected, if flow is in one direction, then the TOF in that direction should be less than the TOF in the opposite direction. Moreover, if there is no flow in pipe 12 (or it is empty), then the UPS and DNS TOF should be equal. Such expectations, however, represent ideal devices and conditions. The present inventors have recognized, however, that various factors will influence these considerations and may cause errors in flow measurement calculations. For example, each of transducers TR.sub.1 and TR.sub.2 has a nominal resonance frequency, so that nominally each also should be excited at a frequency at or closest to that frequency. However, based on manufacturing variability as well as potential changes over time, the actual resonance frequency of a transducer may depart from its nominal value, thereby creating an impedance mismatch. As a result, exciting one transducer with a signal based on the nominally expected resonant frequency will cause a less than optimal response in the received signal of the other transducer. Aging also may affect each transducer and also the channel in pipe 12, which is necessarily part of the overall bi-directionally analyzed system. For example, sedimentation may form or change, thereby changing the overall system resonance frequency. Once more, therefore, if a non-resonance frequency excitation signal or pulse is transmitted through the system, then the response at the receiving transducer, including its sensitivity to any impedance mismatch between the transmitting and receiving transducer, will be less than optimal and, hence, more prone to measure/detection errors. Other factors affecting resonance can include variability of devices for the same manufacturer or for devices from different manufacturers. Still other factors are the medium through which the waves pass and the temperature to which the system (and particularly the transducers) are exposed. In view of these considerations and others as may be ascertained by one skilled in the art, the preferred embodiments implement additional improvements in system 10, as further explored below.

(23) FIG. 3 illustrates a flowchart of a method 30 of operation of processor 14 from system 10, as may be implemented with appropriate instructions stored in the processor and/or computing-device readable media, including hardware, firmware, software, or any combination thereof. By way of introduction, method 30 includes an iterative methodology to excite one transducer at a time, where the excitation signal is a pulse waveform having a series of pulses and the pulse frequency is changed (e.g., increased, but alternatively could be decreased). As a first example, the change in pulse frequency is shown and described as a discrete change in frequency at different points in time, aligned with the period of pulses; in an alternative preferred embodiment example provided thereafter, the change in frequency can be continuous. In either event, a pulse train is generated that sweeps across a range of frequencies, and thereafter, an echo waveform is sampled and correlation methodologies are used to determine the DNS TOF, the UPS TOF, and the difference between those TOFs, after which the velocity and volume of flow determination is made from Equation 1. For clarity, such concepts are first introduced in connection with FIG. 3, with additional details and Figures provided thereafter.

(24) Method 30 commences with a step 32, which establishes a direction index d, which as better understood from the remaining discussion will facilitate frequency pulse transmissions for a first direction (i.e., d=1) followed later by frequency pulse transmissions for a second direction. (i.e., d=2). Thus, in step 32, the direction index d is initialized to a value of one. Next, method 30 continues to a step 34.

(25) Step 34 sets an initial value of a frequency increment value f.sub.INC, where as demonstrated later f.sub.INC increases the frequency at which each pulse in a total number of pulse sets NPS is transmitted. Moreover, in general, the entirety of the NPS pulse sets sweeps across broad band range of frequencies, also shown in step 34 as BW SWEEP. Hence, in step 34, f.sub.INC=BW SWEEP/(NPS−1). For sake of a numeric example, therefore, consider an instance where the total number of pulse sets NPS to be transmitted by a transducer TR.sub.x is 11, with each pulse set having PPS pulses per set (e.g., PPS=4). Consider further that the bandwidth to be swept by the entirety of the NPS=11 pulse sets is 80 kHz. In such a case, then f.sub.INC=BW SWEEP/(NPS−1)=80 kHz/(11−1)=8 kHz. Next, method 30 continues to step 36.

(26) Step 36 initializes a pulse frequency value f.sub.PLS, which as shown later indicates a frequency at which each pulse in a pulse set is to be transmitted. Since step 36 is to initialize f.sub.PLS, where in a preferred embodiment the first (i.e., initialized) value is set at one end of the entire bandwidth through which the total NPS*PPS pulses will sweep. To achieve this setup, step 36 sets f.sub.PLS equal to the nominal, or estimated, or approximated frequency resonance of system 10 minus ½ BW SWEEP. As a numeric example, therefore, assume that system 10 is specified (e.g., by manufacturer or empirical testing) to have a nominal, estimated, or approximated resonance frequency of 160 kHz, and recall the example above where BW SWEEP=80 kHz. Thus, in step 36, f.sub.PLS=160 kHz−½(80 kHz)=120 kHz. Next, method 30 continues to step 38.

(27) Step 38 initializes a pulse set counter PSC to a value of one. As shown below, the pulse set counter PSC provides an index to cause the successive output of each pulse set, up to the total of NPS pulse sets to be transmitted by a transducer TR.sub.x. Next, method 30 continues to step 40.

(28) In step 40, the transducer TR.sub.d, meaning the transducer of system 10 with the index d as its subscript, transmits PPS pulses at the frequency f.sub.PLS. Thus, for the first instance of step 40, that is, where PSC=1, and because step 32 initialized d=1, then transducer TR.sub.1 transmits PPS pulses at the pulse frequency f.sub.PLS established in step 36. By way of illustration, therefore, FIG. 4 depicts a sequence of pulses, each with a respective frequency off f.sub.1, f.sub.2, . . . , f.sub.11, where NPS*PPS (e.g., 11*4=44) pulses are represented given the above exemplary numbers of an implementation of method 30 to transmit a total of 44 pulses. For the first instance of step 40 in method 30, therefore, then transducer TR.sub.1 transmits PPS pulses at frequency f.sub.PLS=f.sub.1. Thus, for the numeric examples given earlier, for this first instance, then f.sub.1=120 kHz. Next, method 30 continues to step 42.

(29) In step 42, a condition is evaluated to determine whether the pulse set count PSC equals the total number pulse sets NPS to be transmitted in a given pulse sequence. If step 42 does not find that PSC=NPS, then method 30 continues to step 44, whereas if PSC=NPS, then method 30 continues to step 46.

(30) In step 44, two operations occur to prepare the inner looping portion of method 30 for its next iteration, that is, to facilitate a next successive pulse set transmission at an increased (or decreased) frequency f.sub.PLS, relative to the pulse frequency from the immediately preceding iteration of the method inner loop. More specifically, step 44 increments the pulse set counter PSC. Additionally, step 44 increases the current frequency f.sub.PLS by the incremental value, f.sub.INC, determined earlier in step 34. Using the numeric examples from above, therefore, recall that f.sub.INC=8 kHz, and note for the first iteration of the first inner loop of method 30, that f.sub.PLS=120 kHz. Thus, for the first iteration where step 44 is reached, then the pulse set count PSC is incremented from one to two, and the pulse frequency f.sub.PLS in increased from 120 kHz to 128 kHz (i.e., f.sub.PLS=f.sub.PLS+f.sub.INC=120 kHz+8 kHz). Following step 44, method 30 returns to a next iteration of step 40, with the immediately-preceding adjustments from step 44. Thus, in step 40, transducer TR.sub.d (e.g., TR.sub.1) transmits a next successive pulse set in a sequence of NPS pulse sets, at the now-increased frequency of f.sub.PLS. Thus, for the second instance of step 40, the second pulse set shown in the time sequence of FIG. 4 is transmitted by transducer TR.sub.1, at a frequency of f.sub.2, where with the numeric examples given, f.sub.2=128 kHz.

(31) From the above, one skilled in the art will appreciate that the method 30 initialization steps and inner loop of steps 40, 42, and 44 perform an iterative methodology so that a first transducer TR.sub.d transmits a total of NPS pulse sets, with each set having PPS pulses at a respective different frequency f.sub.PLS. Moreover, in the example provided, the respective different frequency f.sub.PLS increases linearly among the total NPS pulse sets. The following Table 1 illustrates the entirety of the pulses, again for the example of NOMINAL RESONANCE=160 kHz, and BW SWEEP=80 kHz (and hence, f.sub.INC=8 kHz).

(32) TABLE-US-00001 TABLE 1 ƒ.sub.d ƒ.sub.PLS ƒ.sub.1 120 kHz ƒ.sub.2 128 kHz ƒ.sub.3 136 kHz ƒ.sub.4 144 kHz ƒ.sub.5 152 kHz ƒ.sub.6 160 kHz ƒ.sub.7 168 kHz ƒ.sub.8 176 kHz ƒ.sub.9 184 kHz ƒ.sub.10 192 kHz .sub.ƒ11 200 kHZ
Table 1, therefore, summarizes a general pattern of the linear increase of each pulse set relative to the pulse set(s) either/or preceding or following it, so that the entirety of the NPS pulse sets sweep across an entire broad band desired bandwidth, as may be established with the value, BW SWEEP, centered approximately about an estimated NOMINAL RESONANCE of system 10. Note that an increasing frequency signal as created by the entire pulse train is sometimes referred to as a multi-tone signal. As detailed later, such an approach increases the chances of exciting the system at or near its actual resonance frequency, even though initial that resonance frequency is unknown, and it also results in an improved signal detection in the resultant echo waveform, particularly where correlation is used for one or both of determining the TOF for either the UPS or DNS waveforms, or for determining the delta TOF as between those two.

(33) Having illustrated one approach to providing a multi-tone broad band signal sweep to excite a transducer, note that in an alternative preferred embodiment the manner of creating varying frequency pulses may be adjusted by one skilled in the art. For example, either of the variables PSC of PPS may be adjusted. Indeed, PPS may be set to one, so that the excitation frequency is continuously changed with respect to time, which is sometimes referred to as a chirp signal. As another example, rather than a linear difference between the frequency of successive pulse sets, an exponential difference may be implemented. In general, therefore, the excitation signal may be represented as in the following Equations 2.1 and 2.2:

(34) $\begin{array}{cc}x=\sin \left(2\pi \left(f_1+\frac{k}{2}t\right)t\right)\mathrm{where},k=\frac{f_2-f_1}{T}& \mathrm{Equation}2.1\end{array}$

(35) T is the time it takes to sweep from f.sub.1 to f.sub.2.
excitation pulse=sign(x)  Equation 2.2
Equation 2.2 is intended to create a square wave, where the pulse goes either positive or negative, based on the sign (i.e., positive or negative) of the value x in Equation 2.1. Thus, as the sinusoid of x from Equation 2.1 rises and falls, at the same frequency so shall the square wave excitation pulse resulting from Equation 2.2. Still further, note that alternative preferred embodiments excite the transducer with transducer excitation waveforms other than a pulse train, that is, other types of signals that are periodic and have a known amplitude and frequency, where such signal is further modified per the preferred embodiments to have a frequency that is changing (i.e., either increasing or decreasing) during the excitation period (e.g., between t.sub.0 and t.sub.1 in FIG. 2). For example, an alternative excitation waveform could be a continuous sinusoidal signal of a known amplitude and (changing) frequency coming from a resonant circuit. In any event, therefore, the term waveform in the context of transducer excitation preferably includes such other variations, where for a periodic signal the shape of the period is the same during each period, with it understood that such shape is either expanding or contracting in the time dimension with the change of frequency. Lastly, note also that the sequencing of frequency, whether discrete or continuous, may be controlled and implemented in numerous manners. For example, in one approach, the sequencing may be represented by predetermined data that is stored in a table in, or accessible by, processor 14, so as to apply a corresponding frequency-changing excitation signal to transducers TR.sub.1 and TR.sub.2. In another example, the sequencing may be calculated by processor 14.

(36) Returning to method 30 and its step 42, note the condition of the step is satisfied once all NPS*PPS=44 pulses in a sequence for a transducer TR.sub.d have been transmitted. For example, once the last pulse at frequency f.sub.11 of the 44 pulses of Table 1 have been transmitted, then step 42 is again reached, and in this instances PSC=NPS=11, so method 30 continues to steps 46 and 48. In step 46, the transducer TR.sub.opp, meaning the transducer opposite of transducer TR.sub.d in system 10, receives an echo waveform in response to the pulses transmitted by transducer TR.sub.d. Thus, when transducer TR.sub.1 transmits pulses, then step 46 samples the echo waveform received by transducer TR.sub.2, and in opposite fashion, when transducer TR.sub.2 transmits pulses, then step 46 samples the echo waveform received by transducer TR.sub.1. Thereafter, step 48 determines the absolute TOF.sub.d-opp as between the transmitting transducer TR.sub.d and the receiving transducer TR.sub.opp, where a discussion for a preferred embodiment of such a determination is provided later and, by way of introduction, is preferably achieved by correlating the transmitted waveform with a digital sampling of the received echo waveform. Next, method 30 continues to step 50.

(37) In step 50, a condition is evaluated to determine whether pulse sequences have been transmitted in both the UPS and DNS directions. Specifically, the direction index d is evaluated to determine if it has reached a value of two. If d is less than two, method 30 continues to step 52, where d is incremented and method 30 returns back to step 36, thereby creating an outer loop in method 30. Specifically, step 36 re-establishes f.sub.PLS down to the initial frequency used for pulse f.sub.1 in FIG. 4, and step 38 re-establishes the pulse set count to a value of 1. Next, the inner loop of method 30 will repeat, but here for transmission in the opposite direction as compared to the previous NPS*PPS pulses that were communicated from one transducer to the other. Thus, where the first set of NPS*PPS pulses are transmitted in a first direction from transducer TR.sub.1 to transducer TR.sub.2, then a second iteration (i.e., for d=2) causes a second set of NPS*PPS pulses to be transmitted in a second direction from transducer TR.sub.2 to transducer TR.sub.1, where again each pulse set, after the first pulse set, has a frequency that linearly increases (by f.sub.INC) from the immediately preceding pulse set. Once the NPS*PPS pulses are thereby transmitted in this second direction, steps 46 and 48 again occur, this time sampling the echo waveform received by transducer TR.sub.1 and determining the absolute TOF as between transducer TR.sub.2 and transducer TR.sub.1 Next, method 30 continues to step 54.

(38) Step 54 determines ΔTOF, that is, the difference in TOF determinations made by the two preceding iterations of step 48, that is, the difference between the UPS and DNS TOF determinations. As detailed below, a preferred embodiment for the step 54 determination correlates the samplings of the received DNS and UPS echo waveforms with each other, where the time shift as between the two represents the delta TOF value. As also later explained, because a multitoned (e.g., chirp) signal is used in creating each of the echo waveforms, a more robust and accurate correlation measure is achieved, thereby improving the accuracy and reliability of the delta TOF. Further in this regard, following step 54, step 56 determines flow velocity v, as preferably achieved by substituting the values of UPS and DNS TOF from step 48, and delta TOF from step 54, into Equation 1. Moreover, as each of those substituted values is improved by the processes described herein, then ultimately so is there an improvement in the step 56 velocity v determination.

(39) FIGS. 5A through 5G illustrate a graphical representation of a correlation process used in a preferred embodiment to accomplish step 48 (and step 54) described above. Recall that step 48 determines absolute TOF, that is, TOF as between a pulse train PT transmitted by one transducer TR.sub.d and its corresponding echo waveform received at the opposing transducer TR.sub.opp. By way of introduction to these Figures, therefore, the correlation is a time-shifting step-by-step comparison of two signals to one another, where for each step the comparison is at a different time-shifted position as compared to the preceding/following step. Correlation, if reaching the proper and desired result, determines a point in time when the two compared signals are closest/best aligned in time. As demonstrated by the following discussion, such a time may be used as either the completion, or reference, for determining the absolute TOF period. Also relevant to this discussion as a component of TOF is that it is known that some amount of time will be known to elapse as between the time the pulse train PT starts (e.g., t.sub.0 in FIG. 2) and the time the echo waveform starts to appear (e.g., t.sub.2 in FIG. 2); for sake of discussion, this time is referred to as a preliminary offset time pt.sub.ofs. As also detailed later, note that correlation may be expressed mathematically in various forms, such as an integration (or dot product) over time of one sample set relative to the other, where a maximum value in the mathematical result provides an indication of best alignment.

(40) Looking first to FIG. 5A, it is intended to illustrate a digital sample set having 400 samples of the transmitted pulse train PT and a digital sample set having 400 samples of the received echo waveform EVM. For references, the first and last sample in each set are numbered 1 and 400 respectively, in these Figures. Note that having the same number of samples is used to simplify the present example, whereas in a preferred embodiment the length (i.e., number of samples) is greater for the echo waveform than for the pulse train. The vertical direction in FIG. 5A is intended to depict a time alignment, as further evident below, whereby a sample in one sample set is correlated with a sample in the other sample set that is in the same vertical time alignment. These same representations occur in the remaining FIGS. 5B through 5G, and in the latter Figures each represents a case where the two sample sets are further time shifted relative to one another, representing different time position comparisons of the sample set as part of the correlation process. Looking in more detail to FIG. 5A, following the preliminary offset time offset pt.sub.ofs, the PT sample set and the EVM sample set are aligned in a starting correlation time ct.sub.0. At this time alignment, a correlation measure is made from each sample set with the sample set aligned in the same vertical position; in FIG. 5A, however, none of the 400 samples in the PT sample set are sample aligned with any of the 400 samples in the EVM sample set. Hence, the mathematic implementation of comparison will yield a zero output (ignoring any noise), as each sample set is compared to a zero entity in the same vertical alignment.

(41) Looking next to FIG. 5B, it represents the sample set EVM shifted one time slot to the left, relative to FIG. 5A, thereby occurring at a sequentially-next correlation time ct.sub.1. At this new time alignment, again a correlation measure is made from each sample set with the sample set aligned in the same vertical position; as was the case in FIG. 5A, however, again in FIG. 5B none of the 400 samples in the PT sample set is vertically time aligned with any of the 400 samples in the EVM sample set. Thus, the mathematic comparison will yield a zero output, as again each sample set is compared to a zero entity in the same vertical alignment. Indeed, note that FIG. 5C again illustrates a same result, as FIG. 5C depicts the sample set EVM shifted another single time slot to the left, relative to FIG. 5B, thereby occurring at a sequentially-next correlation time ct.sub.2. As with all other instances in the correlation, at a given time alignment, again a correlation measure is made as between the sample sets and data in the same vertically-aligned position, but once more in FIG. 5C none of the 400 samples in the PT sample set is vertically time aligned with any of the 400 samples in the EVM sample set, thereby again yielding a zero output. One skilled in the art will appreciate that this process will repeat for the next two correlation times (ct.sub.3 and ct.sub.4, not shown), as for those time slots there still will be no time slot overlap as between the two sample sets.

(42) Looking next to FIG. 5D, it represents the sample set EVM shifted three time slots to the left, relative to FIG. 5C, thereby occurring at a correlation time ct.sub.5. At this new time alignment, again a correlation measure is made from each sample set with the sample set aligned in the same vertical position; as was the case in the preceding FIGS. 5A-5C. In FIG. 5D, however, note that the first sample in the EVM sample slot (i.e., sample 1) aligns with the last sample in the PT sample set (i.e., sample 400). Thus, the mathematic comparison from the alignment of these single samples will yield an output greater than zero, but inasmuch as the remaining 399 samples of each set are not aligned to a sample in the opposing set, the total correlation result will still be a relatively low number.

(43) The above single time slot shift and correlation measure continues, as represented in the remaining FIGS. 5E, 5F, and 5G. By way of an arbitrary example, therefore, FIG. 5E represents the sample set EVM sample set shifted fifteen time slots to the left, relative to FIG. 5D, thereby occurring at a correlation time ct.sub.16, thus, in FIG. 5E, more data samples in the respective sets are aligned into respective time slots, which will result in an increased correlation measure, as compared to FIG. 5D. FIG. 5F represents the EVM sample set shifted numerous additional time slots to the left, relative to FIG. 5E, thereby occurring at a correlation time ct.sub.404. FIG. 5F, therefore, is intended to illustrate the point in time when the two sample sets are completely aligned. Again for this time slot alignment, a correlation measure is taken, and the resultant value will ultimately prove to be the highest value of all eventually-completed correlation measures. This highest value, therefore, will confirm that the best match occurred at time ct.sub.404, so the 404 time slots thereby indicated may be added to the preliminary offset time pt.sub.ofs, thereby providing a total absolute TOF as between the pulse train PT and the resulting echo waveform EVM. Further in this regard, lastly FIG. 5G indicates the ongoing correlation process, as the EVM sample set continues to shift to the left for additional time slots, and at each a correlation measure is obtained so as to fully evaluate the potential time alignments of the two sample sets relative to one another. As the shifting continues in this manner, however, the respective correlation measures will be less than that of the fully aligned samples in FIG. 5F.

(44) Having represented and described correlation in general, note that it may be implemented in various manners as ascertainable by one skilled in the art. In one preferred embodiment the following Equation 3 is used, as may be readily programmed into the operation of processor 14:
corr(k)=Σ.sub.i=1.sup.l.sup.2r.sub.1(i+k)r.sub.2(i)  Equation 2
where, r.sub.1 is the received (either UPS or DNS) waveform; r.sub.2 is the reference (either DNS, UPS or TX) waveform; l.sub.1 is the length (number of samples) in r.sub.1; l.sub.2 is the length (number of samples) in r.sub.2; and k=(0 . . . (l.sub.1−l.sub.2), where l.sub.1>l.sub.2).

(45) Having described various preferred embodiment aspects, a key benefit is now observed by returning to FIG. 5F. In FIG. 5F, note that the samples in the pulse train PT sample set necessarily reflect the changing frequency of the chirp signal illustrated in FIG. 4. Somewhere along the changing frequencies of the chirp signal, one frequency will match or most closely approximate the actual resonant frequency of system 10. Thus, even though system 10 has an unknown or changing resonant frequency, the sweeping of various frequencies as shown in FIG. 4 will eventually cause a period in the pulse train to match, or get the closest to, the unknown system resonant frequency. Thus, for the time slot(s) where that frequency occurs, in that time slot, or in one of those time slots if the same frequency occurs in multiple time slots, the EVM sample set will have its largest magnitude, that is, the peak response to the resonant frequency; hence, this may be considered the resonance-occurring time slot. Note further, however, that the sweeping chirp frequency not only ensures a match or near match to the system resonant frequency, but it also improves the determinability of the resonance-occurring time slot. Specifically, with respect to time slots before or after that resonant time slot, as the exciting frequency in each of those slots moves away from the resonant frequency, then the EVM response will immediately start to reduce as compared to the EVM response in the resonance-occurring time slot. Indeed because the preferred embodiment uses the multi-toned (e.g., chirp) pulse train, that reduction will occur not only in the EVM signal, but likewise it also will appear in the correlation measures, that is, the correlation for the resonance-occurring time slot peak will be readily higher than the next highest peaks around it. Thus, correlating to the resonance-occurring time slot should be readily resolvable per the preferred embodiments, so that only a single peak is determined in the correlation, thereby increasing the accuracy of the correlation. A more accurate correlation in turn provides a corresponding accuracy in the determination of the proper resonance-occurring time slot. As introduced earlier, therefore, the timing of the resonance-occurring time slot is then readily added to the preliminary offset time pt.sub.ofs, thereby providing the absolute TOF in the direction being analyzed.

(46) Also in a preferred embodiment, the method illustrated in FIGS. 5A through 5G, and the evaluation of Equation 3, are used to determine the relative TOF, that is, delta TOF is determinable by correlating the UPS and DNS echo waveforms. Thus, in this case, r.sub.1 is one of the UPS or DNS echo waveforms, while r.sub.2 is the other of those echo waveforms. In other words, as those waveforms will be similar as part of the same overall system response, then correlating them with one another will determine the difference in time as between the occurrence of one and the other—hence, delta TOF.

(47) As a further refinement to the preceding, it is noted that correlation of the UPS and DNS echo waveforms, particularly due to the use of the multi-tone excitation signals, is well-suited to identify the particular time slot during which the DNS and UPS waveforms are most closely aligned, as the amplitude of the correlation signal during that time slot should be considerably larger than that of the amplitude in the adjacent time slots. However, within that single time slot, recall that multiple samples have been taken, as the sampling rate is a multiple of the approximate resonance frequency. In this regard, therefore, FIG. 6 illustrates a generally sinusoidal curve SC (shown in dotted line) as a depiction of an analog correlation between the UPS and DNS waveforms. However, in the preferred embodiment, digital sampling and correlation of those samples are implemented, so circles along curve SC represent the samples used in the correlation, where assume for sake of example that the sampling rate is approximately four times the frequency of the waveform, so each illustrated waveform period includes four samples (i.e., four circles). Further, curve SC has a peak amplitude at a time t.sub.p during a period, T.sub.max, so that the peak correlation therefore occurs during T.sub.max. Note, however, that peak time, t.sub.p, is at a distance, δ, from the nearest sample time st.sub.x. Thus, the preferred embodiment determined delta TOF not solely from T.sub.max, but also from the offset of the distance, δ from the peak time, t.sub.p. To achieve this step, the preferred embodiment performs a three point correlation and an interpolation of points during T.sub.max, such as a cosine interpolation, so as to approximate curve SC and thereby determine the peak time, t.sub.p, thereby further establishing the distance, δ, to the nearest sample time. For example, such a correlation and an interpolation may be evaluated from three samples (and their respective sample times) that occur during T.sub.max, as shown in the following Equations 4, 5, 6, 7, and 8:
corr(k)=Σ.sub.i=1.sup.N-kr.sub.1(i+k)r.sub.2(i)  Equation 4
where, r.sub.1 is the received (either UPS or DNS) waveform; r.sub.2 is the received (either DNS or UPS) waveform; N is the length of the UPS and DNS waveforms; and k is the offset between the UPS and DNS waveforms.
Equation 4 provides the correlation in the following Equation 5, so as to provide the three samples, Z.sub.1, Z.sub.0, and Z.sub.1, as further described below:

(48) $\begin{array}{cc}{Z}_n=\mathrm{corr}\left(\hat{k}-n\right)\mathrm{where},n=\left(-1,0,1\right)& \mathrm{Equation}5\\ \omega ={\cos }^{-1}\left(\frac{\left(Z_{-1}+Z_1\right)}{2*Z_0}\right)& \mathrm{Equation}6\\ \phi ={\tan }^{-1}\left(\frac{\left(Z_{-1}-Z_1\right)}{2*Z_0*\sin \left(\omega \right)}\right)& \mathrm{Equation}7\\ \delta =\frac{-\phi }{\omega }& \mathrm{Equation}8\end{array}$
where, Z.sub.0 is the sample closest to t.sub.p; Z.sub.−1 is the sample immediately before Z.sub.0; and Z.sub.1 is the sample after Z.sub.0.

(49) The above-described correlation based methodology accounts for sample slips within a cycle (e.g., T.sub.max) by ensuring Z.sub.0 is always greater than Z.sub.−1 and Z.sub.+1. If Z.sub.0 is not greater than Z.sub.1 and Z.sub.+1 then “n” in Equation 5 can be moved on either direction (−2, −1, 0) or (0, 1, 2) until the condition is satisfied. It is to be noted that delta TOF can be greater than a cycle at high flows, in which case a correlation-based technique cannot correct for a cycle slip. This situation is dependent on the resolvability between the correlation peak and the next highest peaks around it. Hence, the ratio between the correlation peak and the neighboring peak can be used as a measure to evaluate the performance of different excitation pulse techniques, where in the preferred embodiment the use of a multi-tone signal has been found to be very robust and immune to channel affecting issues (e.g., temperature, change in medium) and noise, so as to mitigate the possibility of sample slips and permit the use of the correlation-based technique.

(50) Further in regard to the preceding, FIG. 7 illustrates a partial plot of three cycles from a correlation of UPS and DNS waveforms as described above, where the illustrated partial plots show a main lobe ML, a left side lobe LSL and a right side lobe RSL. The main lobe ML has a peak at line PL.sub.ML, which is readily discernable as distinguishably greater than the side lobe peak line PL.sub.SL that exists for the side lobes LSL and RSL, whichever is higher. Thus, the preferred embodiment manner of frequency sweeping across a certain bandwidth creates a response in the TOF correlation that renders the peak correlation magnitude detectable with less chance of cycle slip, as compared to inadvertently mistaking a side lobe peak as the maximum correlation peak as could occur if a side lobe is close to the amplitude of the center cycle, in which case then the system could confuse the main lobe to be either of the side lobes and thereby causing a cycle slip. Indeed, the ratio between the correlation peak at the main lobe ML and the neighboring peaks to either side (i.e., side lobes LSL and RSL, whichever is greater) can be used as a measure to evaluate the performance of different excitation pulse techniques. Moreover, the ratio also may indicate a benefit in adjusting transmit power so as to improve the ratio, as the ratio provides a measure whether the transmit scheme has provided the needed discernibility in the received signal. The ratio also provides a measure of confidence in the computed TOF. Further in this regard, the following Table 2 depicts in each of its numbered rows an example of excitation of a transducer, where rows 1 and 2 illustrate a prior art single frequency transducer excitation, while rows 3-11 illustrate a preferred embodiment chirp excitation, at the listed center frequency and bandwidth sweep about the center, in respective frequencies. The last column, for each row, lists the ratio between the magnitude of the main lobe ML and the largest side lobe magnitude, which is commonly referred to as the peak to side lobe ratio. Note in this regard, therefore, that row 5 of the Table 2 illustrates a relatively large variation, particularly as compared to the single tone examples in rows 1 and 2 of the Table. Thus, the preferred embodiments demonstrate ample improvement in peak correlation detection, thereby reducing cycle slip and improving the related TOF determination.

(51) TABLE-US-00002 TABLE 2 Peak to Side Lobe Ratio (dB) Raw signal Input (approx. Row Tone Freq (KHz) Pulses 800 codes) 1 Single 160 40 0.069 2 Single 200 40 0.097 3 Chirp 160, BW 20 KHz 40 0.153 4 Chirp 160, BW 40 KHz 40 0.497 5 Chirp 160, BW 80 KHz 40 1.830 6 Chirp 180, BW 20 KHz 40 0.160 7 Chirp 180, BW 40 KHz 40 0.488 8 Chirp 180, BW 80 KHz 40 1.547 9 Chirp 180, BW 20 KHz 40 0.171 10 Chirp 180, BW 40 KHz 40 0.330 11 Chirp 180, BW 80 KHz 40 0.895
Indeed, each of the chip tone excitations yields an improved result over the single frequency excitation approach, with various degrees of improvement as shown in Rows 3-11 of Table 2. Lastly, while the preceding has been described with respect to delta TOF, note that the same preferred embodiment aspects also may be applied to absolute TOF.

(52) In another aspect of a preferred embodiment, with correlation implemented in determining absolute and differential TOF, a search and track approach is implemented to improve efficiency by reducing the number of correlation operations performed over time. Specifically, returning to the example of FIGS. 5A through 5G, in that approach, wherein each sample set included 440 samples, then a first total of 880 correlation measures are taken, so as to compare the two sets relative to one another at each relative possible time slot combination. This approach is referred to as a search correlation. Once the TOF is so determined (either absolute or differential), however, efficiency is further achieved by reducing the number of time slot positions that are evaluated, for subsequent instances when TOF is later determined for subsequent measures. For example, for such subsequent correlation measure, a second total of only ±N (e.g., N=30) time slots around the earlier determined peak time may be evaluated, with there being therefore a strong likelihood that the maximum peak will fall within those second total of time slots. This approach is referred to as a track correlation. Thus, the initial search correlation is later followed by the track correlation. Note further that various criteria may be implemented which, if satisfied, returns the correlation methodology back from track correlation (i.e., limited sample slots) to search correlation (i.e., all sample slots in a sample set). For example, if a peak magnitude found using track correlation is sufficiently different (e.g., by some percentage) from the peak magnitude found prior using a search correlation, then the method can revert back to a search correlation. As another example, time may be a criterion, whereby a search correlation is first performed, and thereafter for a period of time all subsequent correlations are track correlations until the time elapses, whereupon again the correlation reverts to a search correlation and so forth for like subsequent time periods.

(53) Given the preceding, note that a preferred embodiment method 30 applies an excitation signal to a transducer system in a manner that sweeps, such as via a chirp signal, across a broad band range of frequencies, thereby also facilitating a sufficient resonance response in a transducer system that is well detected using correlation for absolute and delta TOF. Note, therefore, that the preferred embodiment is effective, for example, in a transducer system that has either an unknown or changing system response. For example, FIG. 8 illustrates four frequency response plots for a single flowmeter system such as system 10 in FIG. 1, where each plot corresponds to frequency response at a different respective temperature. Frequency is shown along the horizontal axis and detected amplitude response along the vertical axis. For the given system, therefore, it may be seen that for all four plots, a first peak, suggesting a first frequency of resonance, occurs as the frequency sweeps in the range of 160 to 170 kHz. For the same system, a second frequency of resonance, occurs as the frequency sweeps in the range of 210 to 220 kHz. The preferred embodiment chirp signal, therefore, will allow this system to be excited at either (or both) such frequencies, thereby providing the opportunity of proper TOF measurement at such frequencies. Further in this regard, FIG. 9 illustrates a comparable set of four plots also for a single flowmeter system such as system 10 in FIG. 1, but where the transducers have been replaced with another set of transducers, from a same manufacturer. By comparing FIG. 9 to FIG. 8, however, note that the location of peaks (i.e., the frequency response) can change considerably, even among different devices from a same manufacturer and the variation across temperature is also different. Thus, if under the prior art a single excitation frequency were used directed at a system with a response as shown in FIG. 8, that same frequency would be less ideal if used for a system with a response as shown in FIG. 9. In contrast, with the preferred embodiment sweep across a bandwidth of frequencies, the preferred embodiment is equally able to properly and efficiently excite the systems of either FIG. 8 or FIG. 9. Moreover, method 30 may be repeated at various intervals, based on a change of condition, where such a condition can be time, environmental (e.g., temperature), or still others. As a result, in response to changes in the field once a system 10 is implemented, the preferred embodiment sweep will still encompass a frequency (or frequencies) that will maximally excite the system and produce a measurable and accurate response, so as to adopt to changes in the field, without undue manual testing or the rigidity of a single excitation frequency.

(54) From the above, the preferred embodiments are shown to provide a transducer system and method using broad band excitation pulse frequencies for such transducers. The preferred embodiment also may use correlation for TOF measures, further enhanced by the use of the broad band excitation. The preferred embodiments, as a result, provide numerous advantages over the prior art. For example, the preferred embodiments have been empirically shown to provide accurate TOF measures even in noisy environments (i.e., relatively low signal-to-nose ratio (SNR)), with a strong improvement in reducing cycle slips as compared to a single frequency excitation system. Likewise, the preferred embodiments have been empirically shown to provide accurate TOF measures in attenuating media (e.g., methane), also with a strong improvement in reducing cycle slips as compared to a single frequency excitation system. As another example, the preferred embodiments are useful for transducer systems with unknown or changing resonance structure as long as the frequency variation is within the chirp frequency range, as the chirp tone will cause an efficient response regardless of the unknown or unknowability of actual resonant frequency, including if that frequency changes over time, temperature, or other condition. As another example, as long as the system has a wide bandwidth, the preferred embodiment eliminates the need for calibrating the excitation pulse sequence or searching for an excitation frequency at the manufacturer's site, hence significantly reducing calibration overhead. As another result, therefore, real-time in the field calibration (either automatically or manually) is eliminated. As yet another example, even if multiple flowmeters do not have similar frequency responses, the preferred embodiment works similarly on all the flowmeters, hence benefiting the manufacturer by significantly reducing the calibration time. As yet another example, cycle slip issues are also reduced by the preferred embodiment, in response to a bigger correlation peak difference and using the difference in absolute TOF to shift either the upstream or downstream data before computing delta TOF, while the number of correlation computations may be reduced using the search and track aspects of the preferred embodiment, thereby improving the processing demands. Further, the preferred embodiments have been shown in a favorable implementation with respect to flow meters, but note numerous aspects may apply to other systems. For example, preferred embodiment aspects may be applied to transducer applications other than flow meters. As another example, the preferred embodiment is shown by example in a system with two transducers, but various aspects also may be applied to a single transducer, wherein it is excited with a frequency and transmits a pulse train, after which it then responds to the reflection of that pulse train. In view of the above, therefore, while various alternatives have been provided according to the disclosed embodiments, still others are contemplated and yet others can ascertained by one skilled in the art. Given the preceding, therefore, one skilled in the art should further appreciate that while some embodiments have been described in detail, various substitutions, modifications or alterations can be made to the descriptions set forth above without departing from the inventive scope, as is defined by the following claims.