This application claims the priority benefit of Taiwan application serial no. 96120961, filed on Jun. 11, 2007. All disclosure of the Taiwan application is incorporated herein by reference.
1. Field of the Invention
The present invention relates to a method for measuring flow velocity of fluid, and more specifically, relates to an acoustic method for measuring flow velocity of fluid.
2. Description of Related Art
In current medical projects, ultrasound is applied to monitor organs and blood vessels in human body, and the Doppler Effect of ultrasound may be applied to measure the blood flow velocity of a portion of a certain blood vessel, which serves as important information to diagnose blood vessel choke and aneurysm.
According to the above description, the conventional method for calculating the blood flow velocity V may be calculated by the following equation. Doppler shift frequency, and the equation of Doppler shift frequency is shown as the equation (1) below:
In equation (1), fd is Doppler shift frequency, Doppler shift frequency, λ is the wavelength of ultrasonic beam 20a, and θ is Doppler angle which is the angle between the traveling direction E1 of the acoustic beam 20a and the direction of the blood flow velocity V. fd can be acquired via measurement, and λ is a given parameter.
The method of measuring the Doppler angle θ includes observing the blood vessel 10 image formed by ultrasound with naked eyes and manually marking a direction marker M in the sample volume S. The direction marker M is parallel to the blood vessel 10 shown in
However, the Doppler angle θ measured by the above method may be affected by manual observation and result in relatively bigger error. Secondly, the blood vessel 10 shown in
Bd=2|fmax−fd| (2)
Wherein, Bd is Doppler bandwidth. In equation (3), W is the diameter of emitting source of the ultrasonic transducer 20, while F is the focal depth of the ultrasound acoustic beam 20a. W and F are given parameters, and the maximum Doppler bandwidth fmax and Doppler shift bandwidth fd can be obtained from
However, the Doppler shift bandwidth fd in
The present invention provides an acoustic method for measuring the flow velocity of fluid, which can calculate the flow velocity of fluid more accurately.
The present invention provides an acoustic method for measuring flow velocity of fluid comprises the following steps. First, a sound wave is transmitted from a transduction surface of a transducer to a fluid, wherein the transduction surface has at least two signal transduction sections Dn, and n is a positive integer representing the number of the signal transduction sections. Next, when the sound wave is scattered by the fluid and becomes a scattered signal, the scattered signal is received from these transduction sections Dn. According to the scattered signal received by these transduction sections Dn, a plurality of maximum Doppler frequencies fn are calculated. Next, the flow velocity of the fluid is calculated based on these maximum Doppler frequencies fn.
The present invention uses at least two maximum Doppler frequencies to measure the flow velocity of fluid. Compared to the prior arts, the present invention can measure the flow velocity of fluid without knowing the Doppler shift frequency. Therefore the present invention can obtain the flow velocity of fluid accurately without being affected by the shift of Doppler shift frequency.
In order to the make the aforementioned and other objects, features and advantages of the present invention comprehensible, a preferred embodiment accompanied with figures is described in detail below.
The acoustic method for measuring flow velocity of fluid of the present embodiment includes the following steps. First, sound wave U is transmitted to fluid 202 from a transduction surface 112 of the transducer 110, fluid 202, for example, flows in a tube 204. Fluid 202 has the characteristic of scattering sound, wherein the fluid 202 may be blood, water bubble or other fluid that can scatter sound. The transduction surface 112 has at least two signal transduction sections D1 and D2, and sound wave U is emitted from the signal transduction sections D1, D2 through vibrating the transduction surface 112. In addition, sound wave U may be a focused acoustic beam, and the tube 204 may pass through the focal point A of sound wave U, as shown in
In the present embodiment, the area of the signal transduction section D2 may be smaller than the area of the signal transduction section D1, and the signal transduction section D2 is located inside the signal transduction section D1. That is, the range of the signal transduction section D1 covers the whole signal transduction section D2, and the signal transduction section D1 and the signal transduction section D2 overlaps with each other, and the width W2 of the signal transduction section D2 is smaller than the width W1 of the signal transduction section D1.
The transducer 110 may be an annular array transducer, and the shapes of the signal transduction section D1 and the signal transduction section D2 thereof, for example, are circular, as shown in
Since the transducer 110 (or 110′) has two signal transduction sections D1, D2 (or D1′, D2′), therefore the sound wave U includes two acoustic beams U1 and U2 respectively generated by the signal transduction sections D1, D2. Wherein the traveling directions of the acoustic beams U1 and U2 are the same. That is, the acoustic beam U1 and the acoustic beam U2 are pointed to the fluid 204 in tube 204 at the same Doppler angle, i.e. the angles between the traveling directions of the acoustic beams U1, U2 and the tube 204 are the same. The acoustic beam U1 and the acoustic beam U2 overlap each other, and in the focal point A, the width B1 of the acoustic beam U1 is smaller than the width B2 of the acoustic beam U2. Therefore a sample volume S1 can be designated according to the acoustic beam U1, and a sample volume S2 can be designated according to the acoustic beam U2, wherein the sample volume S1 and the sample volume S2 overlap each other, and the width W4 of the sample volume S2 is larger than the width W3 of the sample volume S1, as shown in
Again, referring to
Next, according to the scattered signal received by the signal transduction sections D1, D2, the maximum Doppler frequencies f1 and f2 are calculated. More specifically, when the signal processor 130 receives these electric signals, the signal processor 130 calculates the maximum Doppler frequency f2 according to the scattered signal received by the signal transduction section D2, and calculate the maximum Doppler frequency f1 according to the scattered signal received by the signal transduction section D1.
In the present embodiment, the signal processor 130 includes a spectrum analyzer 132 and a calculator 134. When the signal processor 130 receives the electric signals from the signal receiving generator 120, the spectrum analyzer 132 analyzes the scattered signal received by the signal transduction sections D1 and D2 to generate a Doppler spectrum graph, as shown in
Referring to
Of course, the spectrum analyzer 132 may also use methods other than Fourier transformation to generate the Doppler spectrum graph, therefore herein the emphasized the Doppler spectrum graph shown in
In the present embodiment, a plurality of frequency data of the scattered signal can be obtained from the sample volumes S1, S2 first, for example, the frequency data of 64 to 128 times scattered signal is obtained. Next, the data of these frequencies is subject to Fourier transformation to obtain a plurality of Doppler spectrums. Next, these Doppler spectrums are averaged and normalized to reduce noise, and then more accurate spectrums SP1 and SP2 are further obtained.
After the spectrum analyzer 132 generates the spectrum graph, two maximum Doppler frequencies f1 and f2 are obtained according to the spectrum, wherein the maximum Doppler frequency f1 is obtained according to the scattered signal received by the signal transduction section D1, and the maximum Doppler frequency f2 is obtained according to the scattered signal received by the signal transduction section D2. More specifically, the spectrum analyzer 132 may first set a predetermined threshold T2 (as shown in
Next, the flow velocity of fluid 202 is calculated based on these maximum Doppler frequencies f1 and f2. More specifically, after the spectrum analyzer 132 obtains the maximum Doppler frequencies f1 and f2, the data of the maximum Doppler frequencies f1 and f2 are transmitted to the calculator 134. Next, the calculator 134 may calculate the flow velocity of fluid 202 according to the values of the maximum Doppler frequencies f1 and f2.
In the present embodiment, the flow velocity of fluid 202 is calculated by using the following equation (4):
Wherein V is the flow velocity of fluid 202, λ is the wavelength of the sound wave U, and F is the focal length of the sound wave U. W1 is the width of the signal transduction section D1, and W2 is the width of the signal transduction section D2. λ, F, W1 and W2 are known parameters, therefore if only the maximum Doppler frequencies f1 and f2 are measured, then V in equation (4) can be calculated. That is, the flow velocity of fluid 202 can be calculated.
The equation (4) may be derived as follows:
Bd1=2|f1−fd1| (5a)
Bd2=2|f2−fd2| (5b)
Wherein Bd1 is Doppler bandwidth of spectrum SP1, and Bd2 is Doppler bandwidth of spectrum SP2. The θ is Doppler angle, i.e. the angle between flow velocity of fluid 202 and the traveling direction of sound wave U. fd1 is the frequency to which the maximum intensity value of spectrum SP1 corresponds, while fd2 is the frequency to which the maximum intensity value of spectrum SP2 corresponds. That is, fd1 is an axial component flow velocity (i.e. V cos θ) of fluid 202 measured in the sample volume S1, while fd2 is an axial component flow velocity of fluid 202 measured in the sample volume S2.
Since the axial component flow velocities of fluid 202 measured in sample volume S1 and S2 are measured under the condition of using the same Doppler angle, that is, fd1 and fd2 are same or approximate. Therefore, after subtracting equation (5a) with (5b), the equation (7) may be obtained as shown below:
f1−f2=½(Bd1−Bd2) (7)
Next, subtracting equations (6a) and (6b) may yield equation (8), as shown below:
The Equation (8) may also be expressed as the following equation (9):
Next, after substituting the equation (7) into equation (9), the following equation (10) may be deduced:
The lateral component flow velocity V sin θ of fluid 202 may be calculated using the equation (10). Next, after substituting the equation (10) into equation (6a), the following equation (11) may be deduced through a simple mathematical calculation:
It can be known from the equation of Doppler shift frequency that the axial component flow velocity V cos θ of fluid 202 satisfies the following equation (12), as shown below:
By substituting equation (5a) and equation (11) into the equation (12), the following equation (13) may be deduced as shown below:
Equation (10) and equation (13) may be combined to obtain equation (4) using simple trigonometric function calculation. Therefore, the flow velocity V of fluid 202 can be obtained using the equation (4). In addition, equation (10) and equation (13) may be combined to deduce equation (14) using simple trigonometric function calculation, as shown below:
It can be known from the equation (14) that the Doppler angle θ can also be calculated using the maximum Doppler frequencies f1 and f2, and then the flow direction of fluid 202 in tube 204 can be obtain or the extending direction of tube 204 can be deduced. For example, the equation (14) of the present embodiment may be applied to measure blood flow direction or deduce the extending direction of blood vessel.
In addition, it can be known from above description, although the flow velocity V of fluid 202 can be obtained directly through calculation according to the maximum Doppler frequencies f1 and f2, however in the present embodiment, the Doppler angle θ may also be calculated first according to equation (14). Then, the flow velocity V of fluid 202 may be calculated according to the calculated Doppler angle θ. More specifically, the Doppler shift frequencies fd1 or fd2 may be calculated by using the sample volume S1 or S2 first. Next, the Doppler shift frequency equation, i.e. equation (1) may be used by the calculator 134 to calculate the axial component flow velocity V cos θ of fluid 202 through the Doppler shift frequencies fd1 or fd2. Next, after the maximum Doppler frequencies f1 and f2 are measured, the calculator 134 may calculate the Doppler angle θ according to equation (14). Next, using the known axial component flow velocity V cos θ of fluid 202 and the Doppler angle θ, the flow velocity V of fluid 202 may be calculated.
It should be noted that the acoustic method for measuring flow velocity of fluid revealed in above
In addition, it should be noted that the acoustic method for measuring flow velocity of fluid of the present embodiment may be applied in medical field, for example, in measuring blood flow velocity, or may be applied in industry. In addition, depending upon the different application fields, sound wave U may be ultrasound or regular sound that can be sensed by human ears.
To sum up, the present invention uses a transducer to measure at least two maximum Doppler frequencies, and further calculates the flow velocity of fluid and Doppler angle. Compared to the prior arts, the present invention can measure flow velocity of fluid and Doppler angle without knowing the Doppler shift frequency. Therefore, the present invention can measure the flow velocity of fluid accurately without being influenced by the shift of the Doppler shift frequency, and obtain the Doppler angle. Moreover, when the maximum Doppler frequencies are known, the equation (4) may be used to calculate the flow velocity of fluid, and the equation (14) may be used to calculate the Doppler angle. Thus, according to an embodiment of the present invention, the flow velocity of the fluid may be precisely calculated by a simple mathematical calculation, and also the Doppler angle may be obtained more accurately.
It will be apparent to those skilled in the art that various modifications and variations can be made to the structure of the present invention without departing from the scope or spirit of the invention. In view of the foregoing, it is intended that the present invention cover modifications and variations of this invention provided they fall within the scope of the following claims and their equivalents.
Number | Date | Country | Kind |
---|---|---|---|
96120961 A | Jun 2007 | TW | national |
Number | Name | Date | Kind |
---|---|---|---|
5928153 | Chiang et al. | Jul 1999 | A |
6142944 | Li et al. | Nov 2000 | A |