Method and system of acoustic wave measurement of axial velocity distribution and flow rate

Abstract

A method is provided to measure a distribution of axial velocities and a flowrate in a pipe or a vessel. The method comprises selecting a single cross-section at a stable-flow segment in a pipe or a vessel, installing a plurality of acoustic wave sensors along a peripheral wall of the pipe or the vessel to form a plurality of effective sound wave paths; measuring sound wave travelling time on each sound wave path; substituting the measured sound wave travelling time data into the model formulas based on a sound path refraction principle for reconstruction calculation to obtain a distribution of axial velocity in the measured cross-section of the pipe or the vessel, u(x,y); and integrating the distribution of the axial velocity u(x,y) along the cross-section to obtain a flow rate. A system is also provided to measure an axial velocity distribution and a flow rate in a pipe.

Claims

1. A method of acoustic wave measurement of an axial velocity distribution in a pipe or a vessel, comprising: selecting a single cross-section located at a stable-flow segment in the pipe or the vessel; installing a plurality of acoustic wave sensors along a peripheral wall of the pipe or the vessel at the cross-section, wherein a plurality of effective sound wave paths are formed between the acoustic wave sensors; measuring a sound wave travelling time along each sound wave path by the plurality of acoustic wave sensors; and substituting the sound wave traveling time of each sound wave path into a reconstruction equation to obtain the axial velocity distribution u(x,y) on the cross-section of the pipe or the vessel, wherein the axial velocity distribution u(x,y) is an axial velocity at a point (x, y) in a coordinate system, the reconstruction equation is a function that correlates the axial velocity distribution u(x, y) with each distance between two acoustic wave sensors at two ends of a sound wave path, and a sound traveling time along the sound wave path.

2. The method of claim 1, wherein a formula below is used as the reconstruction equation: $\int_{l_{i}} u (x, y) \sqrt{1 + {(y^{'})}^{2}} d x = L_{i} \sqrt{c^{2} - {(\frac{L_{i}}{Δ t_{i}})}^{2}}, i = 1, 2, .Math., N$ wherein, I.sub.i denotes an ith sound wave path, L.sub.i is a distance between the two acoustic wave sensors at two ends of an ith sound wave path, Δt.sub.i is the sound traveling time along the ith sound wave path, N is a number of effective sound wave paths, c is a static sound speed at measuring physical conditions of a medium in the pipe or the vessel and wherein measuring a sound wave travelling time along each sound wave path includes measuring sound wave travelling times along each sound wave path in opposite directions, respectively, and an average of the sound travelling times along the each sound wave path in opposite directions is a sound wave travelling time in an effective sound wave path.

3. The method of claim 2, wherein each acoustic wave sensor emits sound waves in turn, and when one acoustic wave sensor emits a sound wave, the rest of the acoustic wave sensors record the sound wave to measure the sound wave traveling time and wherein all of the acoustic wave sensors emit sound waves with the same frequency.

4. The method of claim 2, wherein every two to three acoustic wave sensors emit sound waves simultaneously with different frequencies for each, and wherein the sound waves are identifiable by filtering.

5. The method of one of claim 2, wherein the acoustic wave sensor has an integrated function of both emitting a sound wave and receiving acoustic signals, or the acoustic wave sensor is a combination of a sound wave emitter and a sound wave receiver.

6. The method of claim 2, wherein the axial velocity distribution u(x,y) is reconstructed by fitting u(x, y) with a polynomial of Taylor series expansion.

7. The method of claim 2, wherein, when the axial velocity distribution is characterized as a free jet flow, the axial velocity distribution u(x,y) along the cross-section is reconstructed with proximate fitting with a Gaussian formula below: $u (x, y) = U e^{- \frac{{(x - x_{0})}^{2} + {(y - y_{0})}^{2}}{2 σ^{2}}} + u_{0}$ where (x.sub.0, y.sub.0) is coordinates of a point of a maximum velocity, U is the maximum velocity of the jet flow, u.sub.0 is a velocity at a far edge of the jet flow, and σ is a expansion width.

8. The method of claim 2, wherein the measured cross-section in the pipe or the vessel is perpendicular or approximately perpendicular to an axis of the pipe or an axis of the vessel.

9. The method of claim 2, further comprising integrating the axial velocity distribution u(x,y) along the cross-section to obtain a flow rate in the pipe or the vessel.

10. The method of claim 2, wherein the plurality of acoustic wave sensors include four, five, six, seven or eight acoustic sensors.

11. The method of claim 2, wherein the plurality of acoustic wave sensors are installed evenly along the peripheral wall of the pipe or the vessel.

12. The method of claim 2, wherein the plurality of acoustic wave sensors are installed unevenly along the peripheral wall of the pipe or the vessel according to characteristics of an axial flow field.

13. The method of claim 2, wherein the cross-section forms an angle of substantially 90 degrees or an angle greater or less than 90 degrees with an axis of the pipe or the vessel.

14. A system to measure an axial velocity distribution and a flow rate in a pipe or in a vessel, comprising a plurality of acoustic wave sensors disposed on a peripheral wall around a single cross-section of the pipe or the vessel; a digital to analog conversion card; an analog to digital conversion card; a measuring computer including: one or more processors, a memory, and a plurality of instructions stored in the memory and executable by the one or more processors; wherein the digital to analog conversion card is connected to the measuring computer and the plurality of acoustic wave sensors, and configured to transfer digital signals coded by the one or more processors to analog acoustic signals, and the analog acoustic signals are emitted by the plurality of acoustic wave sensors; wherein the analog to digital conversion card is connected to the plurality of acoustic wave sensors and the measuring computer and configured to transfer the analog acoustic signals collected by the plurality of acoustic wave sensors into the digital signals, and input the digital signals into the measuring computer; and wherein the one or more processors are configured to control the plurality of acoustic wave sensors to emit sound waves and instruct the acoustic wave sensors to measure a sound wave travelling time on each sound wave path between one acoustic wave sensor to all the other acoustic wave sensors, and determine the axial velocity distribution by substituting the measured sound wave traveling times into a reconstruction formula below to reconstruct the axial velocity distribution u(x,y) along the cross-section in the pipe or the vessel to obtain the axial velocity distribution via reconstruction, $\int_{l_{i}} u (x, y) \sqrt{1 + {(y^{'})}^{2}} d x = L_{i} \sqrt{c^{2} - {(\frac{L_{i}}{Δ t_{i}})}^{2}}, i = 1, 2, .Math., N$ where, I.sub.i is an ith sound wave path, L.sub.i is a distance between two acoustic wave sensors at two ends of an ith sound wave path, Δt.sub.i is an average of two sound wave traveling times along the ith sound wave path, N is a number of effective sound wave paths, c is a static sound speed at a medium in the pipe or the vessel at measuring physical conditions.

15. The system of claim 14, wherein the one or more processors are further configured to integrates the axial velocity distribution u(x,y) along the cross-section to obtain the flow rate in the pipe or in the vessel.

16. The system of claim 14, wherein the axial velocity distribution u(x,y) is reconstructed by fitting with a polynomial of Taylor series expansion for reconstruction with acoustic measurement.

17. The system of claim 14, wherein, when the axial velocity distribution is characterized as a free jet flow, the axial velocity distribution along the cross-section is reconstructed by proximate fitting with a Gaussian formula: $u (x, y) = U e^{- \frac{{(x - x_{0})}^{2} + {(y - y_{0})}^{2}}{2 σ^{2}}} + u_{0}$ where (x.sub.0, y.sub.0) is coordinates of a point of a maximum velocity, U is the maximum velocity of the jet flow, u.sub.0 is a velocity at a far edge of the jet flow, and σ is a expansion width.

18. The system of claim 14, wherein the plurality of acoustic wave sensors include four, five, six, seven or eight acoustic sensors.

19. The system of claim 14, where only one sound wave path is counted as an effective sound wave path between one pair of the acoustic wave sensors, and wherein measuring a sound wave travelling time along each sound wave path including measuring sound wave travelling times along each sound wave path in opposite directions, respectively and an average of sound travelling times along the each sound wave path in opposite directions is a sound wave travelling time in an effective sound wave path.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a schematic diagram of an axial velocity distribution along one dimension x direction, and a sound wave path, illustrating an effect of a fluid velocity on the sound wave path.

(2) FIG. 2 is a flowchart of a method to determine an axial velocity distribution and a flow rate in a pipe or a vessel according to one embodiment of the present disclosure.

(3) FIG. 3 is a perspective partial view of a pipe, illustrating positioning of acoustic sensors at a cross section of the pipe according to one embodiment of the present disclosure.

(4) FIG. 4 is an axial velocity distribution of a preset simulation flow field in a round cross-section pipe, in which the flow field center resides at the geometric center of the pipe.

(5) FIG. 5 shows the reconstruction relative errors of an axial velocity distribution in a flow field having a center at a geometric round center, using six of acoustic wave sensors to simulate measurement in a round cross section pipe with 3 order fitting precision by Taylor series expansion.

(6) FIG. 6 shows the reconstruction relative errors of an axial velocity distribution in a flow field having a center at the geometric round center, using six of acoustic wave sensors to simulate measurement in a round cross section pipe with 4 order fitting precision by Taylor series expansion.

(7) FIG. 7 is an axial velocity distribution of a preset simulation flow field in a round cross-section pipe, which flow center deviated from the geometric center of the pipe.

(8) FIG. 8 shows the reconstruction relative errors of an axial velocity distribution in a flow field having a center deviated from a geometric round center, using 6 of acoustic wave sensors to simulate measurement in a round-cross section pipe with 3 order fitting precision by Taylor series expansion.

(9) FIG. 9 shows the reconstruction relative errors of the axial velocity distribution in a flow field having a center deviated from a geometric round center, using 6 of acoustic wave sensors to simulate measurement in a round cross section pipe with 4 order fitting precision by Taylor series expansion.

(10) FIG. 10 shows the reconstruction relative errors of the axial velocity distribution in a flow field having a center at a geometric round center, using 8 acoustic wave sensors to simulate measurement in a round cross section pipe with 3 order fitting precision by Taylor series expansion.

(11) FIG. 11 shows the reconstruction relative errors of the axial velocity distribution in a flow field having a center deviated from a geometric round center, using 8 groups of acoustic wave sensors to simulate measurement in a round cross section pipe with 4 order fitting precision by Taylor series expansion.

(12) FIG. 12 is a system to measure an axial velocity distribution and a flow rate in a pipe or a vessel according to one embodiment of the present disclosure.

(13) FIG. 13 is a table listing sound wave travelling time on each sound wave path between six acoustic sensors obtained via integration of an axial velocity distribution of a preset simulation field along the sound wave paths, in which a center of the fluid flow field is at a geometrical center of a round cross section of a pipe and six acoustic sensors are used.

(14) FIG. 14 is a table listing sound wave travelling time on each sound wave path between six acoustic sensors obtained via integration of an axial velocity distribution of a preset simulation field along the sound wave paths, in which a center of the fluid flow field is deviated from a geometrical center of a round cross section of a pipe and six acoustic sensors are used.

DETAILED DESCRIPTION

(15) Various aspects and examples of a method and system of an acoustic measurement of an axial velocity distribution and a flow rate in a pipe or in a vessel are described below and illustrated in the associated drawings. Unless otherwise specified, a method and system of an acoustic measurement of an axial velocity distribution and a flow rate in a pipe or in a vessel in accordance with the present teachings, and/or its various components, may contain at least one of the structures, components, functionalities, and/or variations described, illustrated, and/or incorporated herein. Furthermore, unless specifically excluded, the process steps, structures, components, functionalities, and/or variations described, illustrated, and/or incorporated herein in connection with the present teachings may be included in other similar devices and methods, including being interchangeable between disclosed embodiments. The following description of various examples is merely illustrative in nature and is in no way intended to limit the disclosure, its application, or uses. Additionally, the advantages provided by the examples and embodiments described below are illustrative in nature and not all examples and embodiments provide the same advantages or the same degree of advantages.

(16) FIG. 1 is a schematic diagram of an axial velocity distribution along one dimension x direction, and a sound wave path, illustrating an effect of a fluid velocity on the sound wave path. In FIG. 1, the sound wave path is schematically illustrated as a straight line SP0 in the absence of a fluid flow (a gas flow or a liquid flow) on the sound wave path. When a flow with a velocity Ux is present at one dimension X which is in the path of the sound wave, the sound wave path is bent or refracted to become a curve and the bent sound wave path is illustrated as SP1. The extent of the refraction/bend correlates with the direction and magnitude of the flow intersecting the sound wave path. As the length of the sound wave path changes, the sound wave traveling time between two points on the sound wave path changes. In other words, the sound wave travelling time in a sound wave path varies as the flow or velocity of fluid on the sound propagation path varies. Thus, when the sound wave travelling time on the different sound wave paths in a pipe or a vessel is known, the velocity and flowrate of the fluid in the pipe or the vessel can be determined according to the correlation between the sound wave path and the velocity or flow.

(17) FIG. 2 is a method 100 to determine an axial velocity distribution and a flow rate in a pipe or a vessel according to one embodiment of the present disclosure. The method 100 is further illustrated by an example as shown in FIG. 3 and FIG. 12, in which an axial velocity distribution and a flowrate in a pipe is determined. FIG. 3 is a perspective partial view of a pipe 30, illustrating positioning of acoustic wave sensors 22 at a cross section 24 of the pipe 30 according to one embodiment of the present disclosure. A cross section of a pipe or a vessel illustrating the example positioning of the acoustic wave sensors is also illustrated in FIG. 12, in which a system to measure an axial velocity distribution and a flow rate in a pipe or a vessel is shown. Referring FIG. 2 and with further reference to FIG. 3 and FIG. 12, the method 100 selecting a single cross-section 24 at a flow-stable segment SE or a segment of interest in a pipe or a vessel at 102 and installing a plurality of acoustic wave sensors on a peripheral wall 26 of the pipe 30 at the cross-section 24 at 104. The plurality of sound wave paths are formed between the plurality of acoustic sensors 11. The number of the acoustic wave sensors 22 may be determined according to the complicity of the flow field. In some embodiments, four, five, six, seven or eight acoustic wave sensors are installed or disposed on the peripheral wall 26 at the selected cross section 24 of the pipe 30. In one embodiment, acoustic wave sensor 22 may be a unitary sensor integrating with both functions of acoustic signal emission and reception. In another embodiment, the acoustic wave sensor 22 may be a combination of an acoustic signal emitter and an acoustic signal receiver.

(18) At 106, the method 100 may include measuring a sound wave traveling time on all sound wave paths. Referring to FIG. 12, a plurality of effective sound wave paths are formed between the plurality of acoustic wave sensors 22. In some embodiments, the mutually propagated sound wave paths between a pair of acoustic wave sensor are only counted as one effective sound wave path. Correspondingly, the method 100 includes measuring two sound wave traveling times on each sound wave path in opposite directions of the two acoustic wave sensors 22. In one example, the sensor 22 at the point A measuring the sound wave traveling time from the point D to point A and the sensor 22 at the point D measuring the sound wave traveling time from the point A to the point D. The average sound traveling time of the two sound wave traveling times between two points (e.g., A to D and D to A) wherein the two acoustic sensors are disposed is used as a sound wave traveling time of an effective sound wave path.

(19) At 108, the method 100 may include substituting the sound wave traveling time on each sound wave path into a reconstruction formula to reconstruct an axial velocity distribution to obtain the axial velocity distribution u(x,y). In some embodiments, the following formula (1) is used as a reconstruction formula of an axial velocity distribution:

(20) $\begin{matrix} \int_{l_{i}} u (x, y) \sqrt{1 + {(y^{'})}^{2}} d x = L_{i} \sqrt{c^{2} - {(\frac{L_{i}}{Δ t_{i}})}^{2}}, i = 1, 2, .Math., N & (1) \end{matrix}$

(21) where, I.sub.i is an ith sound wave path, L.sub.i is a distance between the two acoustic wave sensors at the ith effective sound wave path, Δt.sub.i is the sound wave traveling time along the ith sound wave path, Nis the number of effective sound wave paths, c is the sound traveling speed in a static medium in the pipe or the vessel at the measuring physical conditions.

(22) In some embodiments, the method 100 may further includes integrating the axial velocity u(x,y) obtained from the reconstruction along the cross-section 24 to determine the flow rate in the pipe 30, for example. In some embodiments, every acoustic wave sensor 22 emits sound wave in turn, and the frequency of the acoustic signal from each acoustic wave sensor 22 is the same. When one acoustic wave sensor 22 emits a sound wave, the rest of other acoustic wave sensors record this sound wave. Typically, the acoustic measurement by the acoustic wave sensor 22 takes approximately 2 seconds to complete. Since the flow field is relatively steady, the effect of such sound emitting method on the measurement precision is very small.

(23) In some embodiments, every two to three acoustic wave sensors 22 emit sound waves simultaneously, and each acoustic wave sensor 22 emits sound wave at different frequencies. For example, referring to FIG. 12, the acoustic wave sensors 22 at the points A, B and C emit the sound waves simultaneously first and then the acoustic wave sensors 22 at the points D, E and F emit the sound waves simultaneously. It should be appreciated that any of the acoustic wave sensors may be formed as a group to emit sound wave simultaneously. Compared with the method that every acoustic wave sensor 22 emits sound wave in turn, this approach takes shorter time to finish one complete measurement on the axial velocity distribution in the pipe or the vessel. However, using different frequencies simultaneously in one measurement may have difficulty for effective identification of acoustic signals from each other.

(24) In practical application, whether to control the acoustic wave sensor to emit the sound wave in turn or control every 2-3 acoustic wave sensors to emit sound wave simultaneously can be determined according to the time required to complete the measurement.

(25) It should be noted that the suitable measurement acoustic frequency can be selected based on the size of the measured object and the properties of the flow medium. For gas or multiphase medium, the attenuation rate of sound wave is proportional to the square of the frequency, hence the sound wave would attenuate quickly with increase of the sound wave frequency. However, the larger the frequency is, the more favorable the condition to achieve high precision of obtaining the sound wave travelling time, and thus high precision of measurement results of the axial velocity distribution. A high measurement sound frequency should be used as long as the sufficient intensity of acoustic signals at the ends of the sound paths can be achieved so as to ensure the measurement accuracy.

(26) In order to diminish the influence of the radial component of the flow velocity within the measured cross-section in the pipe, the average of the two sound wave travelling times along each effective sound wave path in opposite directions may be taken as the sound wave traveling time in the reconstruction.

(27) In some embodiments, 4, 5, 6, 7, or 8 acoustic wave sensors may be evenly installed around the measured cross-section to form six (6), ten (10), fifteen (15), twenty one (21), and twenty eight (28) sound wave effective paths, respectively. In the example illustrated in FIG. 12, the six (6) acoustic sensors 22 are disposed on the peripheral wall at the cross section, which forms 15 sound wave paths, i.e., the total of fifteen (15) sound paths between points A and B, C and B, D and C, E and D, E and F, F and A, E and A, A and C, E and C, F and B, D and B, F and D, D and A, E and B and F and C, respectively.

(28) In some embodiments, the axial velocity distribution u(x,y) is fitted with Taylor series expansion for reconstruction. The 3.sup.rd order Taylor series expansion below may be used for reconstruction, for example,
u(x,y)=C.sub.1+C.sub.2x+C.sub.3y+C.sub.4x.sup.2+C.sub.5xy+C.sub.6y.sup.2+C.sub.7x.sup.3+C.sub.8x.sup.2y+C.sub.9xy.sup.2+C.sub.10y.sup.3+o(x,y) (2)
where o(x,y) is the infinitely small quantity term, C.sub.1, C.sub.2, . . . ,C.sub.10 are the polynomial coefficients to be determined. The precision of the fitting polynomial, and correspondingly the number of the terms, is determined by the complexity of the axial flow field. Correspondingly, the number of effective sound wave paths formed by the acoustic wave sensors should be no less than the number of the fitting polynomial coefficients.

(29) In some embodiments, if the axial velocity distribution u(x,y) accords with characteristics of a free jet flow, a Gaussian polynomial below may be used for reconstruction:

(30) $\begin{matrix} u (x, y) = U e^{- \frac{{(x - x_{0})}^{2} + {(y - y_{0})}^{2}}{2 σ^{2}}} + u_{0}, & (3) \end{matrix}$
where (x.sub.0, y.sub.0) is the coordinates of the point of maximum velocity, U is the maximum velocity magnitude of the jet flow with the Gaussian distribution, u.sub.0 is the base velocity at the far edge of the jet away from the center of the jet flow, and a is the expansion width (corresponding to the standard deviation of the Gaussian function). The number of effective sound wave paths determined by the corresponding number of acoustic wave sensors should be no less than five (5) such that ten (10) propagation sound wave paths are formed between each pair of the acoustic wave sensors.

(31) In some embodiments, the cross-section 24 in a pipe 30 or a vessel intersects an axis AP of a pipe or a vessel with an angle α (as shown in the example of the pipe 30 in FIG. 3). Referring to FIG. 3 and FIG. 12, the selected cross-section 24 may be perpendicular or substantially perpendicular to the axis AP of the pipe 30 or the vessel, that is the angle α is substantially 90 degrees. In some embodiments, the angle α may be greater or less than 90 degrees. In some embodiments, the acoustic wave sensors 22 may be installed around the cross section 24 of the peripheral wall 26 evenly. In some embodiments, the acoustic wave sensors may be installed unevenly around the peripheral wall of the measured cross-section according to the characteristics of the axial flow profile such that more sound wave paths are obtained for the flow field with relatively complex variations of axial velocity.

(32) It should be appreciated that the shape of the cross-section of the measured object, as well as certain variations of the reconstruction formulas for the axial velocity distribution would not influence the effect of the measurement method and system of the present disclosure.

(33) The measurement methods of the present disclosure are suitable for the pipes and vessels with different cross-section shapes such as circular, elliptical, and rectangular etc. Therefore, reconstructing the axial velocity distribution for a circular cross-sectional pipe is used as an example to validate the acoustic measurement method of this invention.

(34) As shown in FIG. 3 and FIG. 12, a cross-section 24 within a stable-flow segment in a circular cross-sectional pipe is selected, and six (6) acoustic wave sensors 22 are disposed around a peripheral wall 26 of the pipe 30 at the measured cross-section 24 of the pipe, Other conditions include: the diameter of the pipe is one (1) meter, the medium is air, and the temperature is room temperature. The preset axial velocity distribution is reconstructed by the method of this invention using a simulated acoustic measurement.

(35) Taking the geometric center of the cross-section of the pipe as the coordinate origin, the preset simulation axial velocity field of the round cross section in the pipe with a center of the flow field at the geometric center is constructed as follows:
u.sub.0(x,y)=5e.sup.−(x.sup.2.sup.+y.sup.2.sup.) (4)

(36) u.sub.0 (x, y) is illustrated in a coordinate system as shown in FIG. 4.

(37) Substituting formula (4) into formula (1) and integrating along sound wave paths as shown in FIG. 12, all sound wave travelling times on the sound wave paths are obtained. The sound wave paths shown in FIG. 12 and corresponding sound travelling times (ms) in the simulated pipe are listed in Table 1 in FIG. 13, wherein 6 digits after decimal point are kept.

(38) Based on the acoustic method 100 to determine the axial velocity distribution and flow rate, and by using a 3 order precision fitting by Taylor series expansion to calculate the axial velocity distribution u(x,y) and the reconstruction formula (1), the results of reconstructed flow field u(x,y) with the acoustic simulation measurement is obtained. The relative reconstruction errors of u(x,y) using the method 100 to the preset simulation axial velocity field are shown in FIG. 5. FIG. 5 shows the reconstruction relative errors of an axial velocity distribution in a flow field having a center at a geometric round center, using six of acoustic wave sensors to simulate measurement in a round cross section pipe with 3 order fitting precision by Taylor series expansion. As can be seen in FIG. 5, the relative error of the reconstructed u(x,y) using the method 100 to the preset simulation axial velocity field is approximately within a range of −0.5463 to 0.2872.

(39) Similarly, the reconstructed axial velocity distribution u(x,y) with a 4 order precision Taylor series expansion is conducted, and the relative reconstruction error is shown in FIG. 6. As can be seen in FIG. 6, the relative error of the reconstructed u(x,y) using the method 100 to the preset simulation axial velocity field is approximately in a range of −0.0088 to 0.0025.

(40) As can be seen, the simulated acoustic measurement is sufficiently accurate and reliable. And when the fitting precision with Taylor series expansion is increased, the simulated measurement precision is also increased correspondingly.

(41) Similarly, in another example, a preset simulation field may be built with the flow center deviating from the geometric center in a round cross section of a pipe, and a center of the round cross section is used as an origin of a coordinate. The preset axial velocity distribution in the pipe is determined using the following formula:
u.sub.0(x,y)=5e.sup.−(1.5(x−0.1).sup.2.sup.+(y+0.15).sup.2.sup.) (5)

(42) The preset simulation field of the formula (5) is illustrated in FIG. 7.

(43) By substituting formula (5) into formula (1) and integrating along sound wave paths shown in FIG. 12, all sound wave travelling times (ms) are obtained and listed in Table 2 in FIG. 14, wherein 6 digits are remained behind a decimal point.

(44) Based on the method 100 to measure the axial velocity distribution and flow rate, and using a precision fitting of 3 order (Formula (2)) and 4 order of Taylor series, respectively, to calculate the axial velocity distribution u(x,y), the reconstructed flow field result u(x,y) with the acoustic simulation measurement using method 100 is obtained. The reconstruction errors are given in FIGS. 8 and 9, respectively. As can be seen in FIG. 8, the relative error of the reconstructed u(x,y) using the method 100 (3 order of Taylor series) to the preset simulation axial velocity field is approximately in a range of −1.5924 to 0.7029. As can be seen in FIG. 9, the relative error of the reconstructed u(x,y) using the method 100 (4 order of Taylor series) to the preset simulation axial velocity field is approximately in a range of −0.1030 to 0.3008.

(45) Again, it can be seen that the simulated acoustic measurement is sufficiently accurate and reliable. And when the fitting precision with Taylor series expansion is increased, the simulated measurement precision by the acoustic measurement is increased correspondingly. That is, the precision of acoustic measurement is affected by the fitting order/precision of Taylor series expansion.

(46) To further determine the influence factors to the precision of the acoustic measurement, the simulated acoustic measurements in a round cross section of a pipe was conducted with eight (8) acoustic wave sensors. The preset axial velocity distribution in pipe is determined using the formula (5). That is, the simulated flow field is built with the flow center deviating from the geometric center in a round cross section of a pipe, and a center of the round cross section is used as an origin of a coordinate. Taylor series expansion fitting precision of 3 order and Taylor series expansion fitting precision 4 order are used for reconstruction. The relative reconstruction errors of the axial velocity distribution for the acoustic measurement using Taylor series expansion fitting precision of 3 order is illustrated in FIG. 10. As can be seen in FIG. 10 for fitting with 3 order of Taylor series, the relative error of the reconstructed u(x,y) using the method 100 to the preset simulation axial velocity field is approximately in a range of −2.4172 to 0.5419.

(47) The relative reconstruction errors of the axial velocity distribution for the acoustic measurement using Taylor series expansion fitting precision of 4 order is illustrated in FIG. 11. As can be seen in FIG. 11 for fitting with 4 order of Taylor series, the relative error of the reconstructed u(x,y) using the method 100 to the preset simulation axial velocity field is approximately in a range of −0.119 to 0.2895.

(48) As it can be seen, the precision of the acoustic wave measurement is influenced by the fitting order/precision with Taylor series expansion. The use of eight of acoustic wave sensors does not show significant improvement to the simulated measurement results compared with using 6 of acoustic wave sensors.

(49) The above two typical exemplary scenarios of simulated acoustic measurements of the axial velocity distribution show that the acoustic measurement of axial velocity distribution in a pipe is feasible and also reliable, and the acoustic measurement is based on the mechanism of sound wave path bending caused by a flow perpendicular or having angle to the sound propagation direction. Further, the precision of the acoustic measurement is affected by the precision of order of fitting series and irrelevant to the number of acoustic wave sensors used when a certain precision/accuracy has been achieved. In other words, when the number of the acoustic wave sensors is greater than a certain number, the effect of the number of the acoustic wave sensors to the measurement accuracy is small. Such a certain number is determined by the complexity of the spatial variation of the axial flow field, with a larger value for a more complex axial flow field. Nevertheless, as the fitting order/precision of series is increased, the number of acoustic wave sensors needs to be increased correspondingly to solve the equation because of the increase of number of coefficients needed in the fitting formula for the axial velocity distribution.

(50) According to another aspect, the present disclosure provides a system 20 for acoustic measurement of the axial velocity distribution and the flow rate in a pipe 30 or a vessel. Referring to FIG. 12, the system 20 comprises a plurality of acoustic wave sensors 22 installed on the peripheral wall 26 of a measured cross-section 24 in the pipe 30, a digital to analog conversion card (AD card) 40, an analog to digital conversion card (DA card) 50, and a measuring computer 60. The measuring computer 60 includes a preset measurement software to measure the axial velocity distribution and the flow rate as described above. For example, the measuring computer 60 may include one or more processor, a memory and a plurality of instructions stored in the memory executable by the one or more processor to implement one or more steps of the method 100 or the routines to determine the axial velocity distribution and flowrate in the pipe or the vessel.

(51) The digital to analog conversion card 40 is connected with the measuring computer 60 and the plurality of the acoustic wave sensors 22, respectively, to transfer the digital signals coded with the measurement software or controlled by the one or more processors into analog acoustic signals to be emitted by the acoustic wave sensors 22.

(52) The analog to digital conversion card 50 is connected with the measuring computer 60 and the acoustic wave sensors 22, respectively, to transfer the collected acoustic signals into digital signals and input into the measuring computer 60.

(53) The measuring computer 60 controls, via the measurement software or the one or processors and the memory, the acoustic wave sensors 22 to emit sound waves and measure the sound wave travelling time on the sound wave paths from each acoustic wave sensor 22 to all other acoustic wave sensors 22, and substitute the sound wave travelling times obtained from the acoustic wave sensors 22 into a reconstruction formula for the axial velocity distribution to obtain the axial velocity distribution. In some embodiments, the reconstruction formula (1) for an axial velocity distribution in a pipe or a vessel is used to obtain the axial flow distribution u(x,y):

(54) $\begin{matrix} \int_{l_{i}} u (x, y) \sqrt{1 + {(y^{'})}^{2}} d x = L_{i} \sqrt{c^{2} - {(\frac{L_{i}}{Δ t_{i}})}^{2}}, i = 1, 2, .Math., N & (1) \end{matrix}$
wherein, I.sub.i is the ith sound wave path, L.sub.i is a distance between the two acoustic sensors at the two ends of the ith path, Δt.sub.i is the average of the two sound wave travelling times of the sound waves along the ith sound path in opposite directions, N is the number of effective wave sound paths, c is the static sound speed at the measuring medium in the pipe or the vessel at the measuring conditions;

(55) The measuring computer, via the preset measuring software or one or more processors and memory, integrates the obtained axial velocity distribution u(x,y) along the cross-section, so that the flow rate is obtained.

(56) The method and system of the present disclosure are suitable to different media in a pipe, such as gas, liquid, two-phase or multiphase media for effective measurements to the axial velocity distribution and the flow rate in the pipe. Furthermore, the method and system of the present disclosure can be used to measure flow fields in a combustion chamber, a fluidized bed, a chemical reactor, and an open jet flow, etc.

Selected Embodiments and Claim Concepts

(57) This section describes additional aspects and features of the method and system to determine the axial velocity distribution and flow rate presented without limitation as a series of paragraphs, some or all of which may be alphanumerically designated for clarity and efficiency. Each of these paragraphs can be combined with one or more other paragraphs, and/or with disclosure from elsewhere in this application, in any suitable manner. Some of the paragraphs below may expressly refer to and further limit other paragraphs, providing without limitation examples of some of the suitable combinations.

(58) A0. A method to measure an axial velocity distribution and a flow rate in a pipe, characterized, comprising:

(59) selecting a single cross-section located at a stable-flow segment or an interested segment in the pipe or a vessel, installing a plurality of acoustic wave sensors along a peripheral wall of the pipe or the vessel at the cross-section, wherein a plurality of effective sound wave paths are formed between the acoustic wave sensors and only one sound wave path is counted as an effective sound wave path between one pair of the acoustic wave sensors;

(60) measuring a sound wave travelling time along each sound wave path, and substituting the sound wave traveling time into the following equation so as to reconstruct the axial velocity distribution u(x,y) on the cross section of the pipe or the vessel;

(61) $\int_{l_{i}} u (x, y) \sqrt{1 + {(y^{'})}^{2}} d x = L_{i} \sqrt{c^{2} - {(\frac{L_{i}}{Δ t_{i}})}^{2}}, i = 1, 2, .Math., N$

(62) wherein, I.sub.i is an ith sound wave path, L.sub.i is a distance between the two acoustic wave sensors at two ends of an ith sound wave path, Δt.sub.i is the sound traveling time along the ith sound wave path, N is a number of effective sound wave paths, c is a static sound speed at measuring physical conditions of a medium;

(63) integrating the axial velocity distribution u(x,y) obtained from the above step along the cross-section to obtain flow rate in the pipe.

(64) A1. The method of paragraph A0, wherein measuring a sound wave travelling time along each sound wave path including measuring sound wave travelling times along each sound wave path in opposite directions, respectively and an average of the sound travelling times along the each sound wave path in opposite directions is a sound wave travelling time in an effective sound wave path.

(65) A2. The method of any one of paragraph A0 to A1, wherein each acoustic wave sensor emits sound waves in turn, and when one acoustic wave sensor emits a sound wave, the rest of the acoustic wave sensors record the sound wave and wherein all of the acoustic wave sensors emit sound waves with the same frequency.

(66) A3. The method of any one of paragraph A0 to A1, wherein every two to three acoustic wave sensors emit sound waves simultaneously with different frequencies for each, and wherein the sound waves are identifiable by filtering.

(67) A4. The method of any one of paragraph A0 to A1, wherein the acoustic wave sensor has an integrated function of both emitting sound wave and receiving acoustic signals.

(68) A5. The method of any one of paragraph A0 to A1, wherein the acoustic wave sensor is a combination of a sound wave emitter and a sound wave receiver.

(69) A6. The method of any one of paragraph A0 to A1, wherein the axial velocity distribution u(x,y) is calculated by fitting with a polynomial of Taylor series expansion.

(70) A7. The method of any one of paragraph A0 to A1, wherein, when the axial velocity distribution is characterized as a free jet flow, the axial velocity distribution along the cross-section is calculated with proximate fitting:

(71) $u (x, y) = U e^{- \frac{{(x - x_{0})}^{2} + {(y - y_{0})}^{2}}{2 σ^{2}}} + u_{0}$

(72) where (x.sub.0, y.sub.0) are coordinates of a point of a maximum velocity, U is the maximum velocity magnitude of the jet, u.sub.0 is a velocity at a far edge of the jet flow, and a is an expansion width.

(73) A8. The method of any one of paragraph A0 to A1, wherein the measured cross-section in the pipe is perpendicular or approximately perpendicular to an axis of the pipe.

(74) B0. A system based on one of paragraph A0 to paragraph A9, the system comprises: several acoustic wave sensors installed on a peripheral wall around a measured cross-section in a pipe or a vessel, a digital to analog conversion card, an analog to digital conversion card, and a measuring computer for measurement; a measuring software is preset in the measuring computer including a measurement software to measure an axial velocity distribution in the pipe, and the flow rate according to the previous described method;

(75) the digital to analog conversion card is connected to the measuring computer and the acoustic wave sensors, respectively to transfer the digital acoustic signal, coded by the measurement software to analog acoustic signals, and emits by acoustic wave sensors;

(76) the analog to digital card is, connected to the acoustic wave sensors and the measuring computer, respectively, to transfer the measurement acoustic signals collected by the acoustic wave sensors into digital signals, and input into the measuring computer; wherein

(77) via the measuring software, the measurement computer controls all acoustic wave sensors to emit sound waves, measures each sound wave travelling time in each sound wave path from one acoustic wave sensor to all the other acoustic wave sensors, and substitutes all sound wave traveling time along sound wave paths into a reconstruction formulas to reconstruct the axial velocity distribution u(x,y) along the cross-section in the pipe or vessel,

(78) $\int_{l_{i}} u (x, y) \sqrt{1 + {(y^{'})}^{2}} d x = L_{i} \sqrt{c^{2} - {(\frac{L_{i}}{Δ t_{i}})}^{2}}, i = 1, 2, .Math., N$

(79) where, I.sub.i is an ith sound wave path, L.sub.i is a distance between two acoustic wave sensors at two ends of an ith path, Δt.sub.i is an average of two sound wave traveling time of along the ith sound wave path, N is the number of effective sound wave paths, c is the sound speed at a medium in the pipe at measuring physical conditions; and

(80) wherein the measuring software of the measurement computer integrates the axial velocity u(x,y) along the cross-section to obtain the flow rate in the pipe.

(81) It should be noted that the above exemplary embodiments are used to describe the technical solutions of the present invention, and are not restricted to the exemplary embodiments. Although the invention is described in detail with reference to the preferred embodiments, the technicians in the related arts should understand that: the exemplary embodiments of the present invention can be modified or the technical characteristic can be substituted without depart from the spirit of the present invention, and which should be covered in the scope of the technical solutions of the present invention.

Method and system of acoustic wave measurement of axial velocity distribution and flow rate

Assignee

Inventors

Cpc classification

Classification Explorer

G01F1/667

PHYSICS

Classification Explorer

G01P5/242

PHYSICS

Classification Explorer

G01P5/245

PHYSICS

Classification Explorer

G01H5/00

PHYSICS

Classification Explorer

G01P5/001

PHYSICS

Classification Explorer

G01F1/662

PHYSICS

Classification Explorer

G01P5/24

PHYSICS

Classification Explorer

G01F7/00

PHYSICS

International classification

Classification Explorer

G01P5/24

PHYSICS

Classification Explorer

G01F1/66

PHYSICS

Classification Explorer

G01F7/00

PHYSICS

Classification Explorer

G01H5/00

PHYSICS

Abstract

Claims

Description