CN111157064B - Multi-frequency information fusion continuous wave flow measuring method and device and electronic equipment - Google Patents
Multi-frequency information fusion continuous wave flow measuring method and device and electronic equipment Download PDFInfo
- Publication number
- CN111157064B CN111157064B CN201911403014.7A CN201911403014A CN111157064B CN 111157064 B CN111157064 B CN 111157064B CN 201911403014 A CN201911403014 A CN 201911403014A CN 111157064 B CN111157064 B CN 111157064B
- Authority
- CN
- China
- Prior art keywords
- frequency
- measurement
- side tone
- frequencies
- minimum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/66—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters
Abstract
The invention discloses a method, a device and electronic equipment for measuring continuous wave flow by multi-frequency information fusion, wherein the method comprises the following steps: according to the estimated fluid flow, the minimum lateral sound frequency is obtained when the difference between the phase change value of the ultrasonic wave propagating in the fluid along the countercurrent direction and the phase change value of the ultrasonic wave propagating in the fluid along the downstream direction is not more than 180 degrees; selecting a plurality of groups of measuring side tone frequencies and calculating side tone frequencies based on the minimum side tone frequency; according to the measured phase difference value under each group of the measured side tone frequency, the actual phase difference value corresponding to each group of the measured side tone frequency and the calculated side tone frequency is obtained; respectively calculating the fluid flow corresponding to each group of side-tone frequencies for measurement according to each group of side-tone frequencies for measurement and the corresponding actual phase difference value; and averaging the fluid flow corresponding to the lateral sound frequency for each group of measurement to obtain the final fluid flow. The invention can improve the measurement precision and ensure that the measurement sensitivity is not reduced.
Description
Technical Field
The present invention relates to the field of ultrasonic flow measurement, and in particular, to a method and an apparatus for continuous wave flow measurement with multi-frequency information fusion, and an electronic device.
Background
The ultrasonic flowmeter obtains the average flow velocity information of the pipeline by processing ultrasonic signals by utilizing the obvious difference of forward and backward wave propagation of sound waves in the flowing of the pipeline, and then predicts the flowing flow of the pipeline. The ultrasonic flowmeter has the advantages of no invasion of a measured fluid, no moving part, no influence on fluid flow and the like, thereby being widely applied to various industrial fields.
In the invention patent with patent number ZL2015102172276 and invention name "ultrasonic flow measuring method and device based on lateral sound phase measurement", the author proposes an improved lateral sound phase measurement method aiming at the problem of complex system design in the traditional lateral sound, thereby reducing the system complexity, however, the method only uses the phase value of a single frequency to obtain the flow, thereby neglecting the phase values of other frequencies, the data utilization is not sufficient, and the obtained flow precision is still to be improved.
Disclosure of Invention
The invention provides a continuous wave flow measuring method and device with multi-frequency information fusion and electronic equipment, and aims to solve the technical problems that data is not sufficiently utilized and the flow precision needs to be improved in the conventional ultrasonic flow measuring method.
The technical scheme adopted by the invention is as follows:
a multi-frequency information fusion continuous wave flow measuring method comprises the following steps:
according to the estimated fluid flow, the minimum lateral sound frequency is obtained when the difference between the phase change value of the ultrasonic wave propagating in the fluid along the countercurrent direction and the phase change value of the ultrasonic wave propagating in the fluid along the downstream direction is not more than 180 degrees;
selecting a plurality of groups of measuring side tone frequencies and calculating side tone frequencies based on the minimum side tone frequency, wherein each group of calculating side tone frequencies is equal to the difference between the maximum measuring side tone frequency and the rest measuring side tone frequencies, and the minimum calculating side tone frequency does not exceed the minimum side tone frequency;
according to the measured phase difference value under each group of the measured side tone frequency, the actual phase difference value corresponding to each group of the measured side tone frequency and the calculated side tone frequency is obtained;
respectively calculating the fluid flow corresponding to each group of side-tone frequencies for measurement according to each group of side-tone frequencies for measurement and the corresponding actual phase difference value;
and averaging the fluid flow corresponding to the lateral sound frequency for each group of measurement to obtain the final fluid flow.
In an embodiment, the selecting a plurality of sets of side tone frequencies for measurement and side tone frequencies for calculation based on the minimum side tone frequency, wherein each set of side tone frequencies for calculation is equal to a difference between a maximum side tone frequency for measurement and the remaining side tone frequencies for measurement, and the minimum side tone frequency for calculation does not exceed the minimum side tone frequency specifically includes the step;
selecting a group of sidetone frequencies not exceeding the minimum sidetone frequency as a minimum side tone frequency for calculation;
setting a preset maximum set sidetone frequency greater than the minimum sidetone frequency as a maximum measurement sidetone frequency;
determining the upper limit value of progressive multiples between adjacent groups of side tone frequencies according to the absolute error of ultrasonic phase measurement;
determining other side tone frequencies for calculation between the minimum side tone frequency for calculation and the maximum side tone frequency for measurement according to the upper limit value of the progressive multiple between adjacent groups of side tone frequencies;
and determining the frequencies of the rest groups of side tones for measurement according to the difference between the frequency of the maximum side tone for measurement and the frequencies of all the side tones for calculation.
In one embodiment, the progressive multiple k between adjacent sets of sidetone frequencies satisfies:
where δ represents the absolute error of the ultrasonic phase measurement.
In one embodiment, the number of sets N of sidetone frequencies for measurement satisfies:
wherein, TmaxFor single measurement of maximum time constraint, UmaxMeasuring the maximum range of flow rate for the flowmeter, fmaxThe maximum frequency of the flowmeter, L is the length of a pipeline in the flowmeter, and C is the propagation speed of the ultrasonic wave in the static fluid.
In an embodiment, calculating the fluid flow rate corresponding to each group of measuring side-tone frequencies according to each group of measuring side-tone frequencies and the corresponding actual phase difference values thereof specifically includes:
wherein, ViRepresenting the sidetone frequency f for measurementiCorresponding fluid flow, L is the length of the pipeline in the flowmeter, R is the radius of the pipeline, C is the propagation speed of the ultrasonic wave in the static fluid, phidiff(fi) For measuring the frequency f of the sidetoneiThe corresponding actual phase difference value.
In an embodiment, the averaging of the fluid flows corresponding to the side-tone frequencies for each group of measurements to obtain the final fluid flow V specifically includes:
according to another aspect of the present invention, there is also provided a continuous wave flow measuring apparatus with multi-frequency information fusion, including:
the minimum side-tone frequency determining module is used for solving the minimum side-tone frequency when the difference value between the phase change value of the ultrasonic wave in the fluid in the upstream direction and the phase change value of the ultrasonic wave in the fluid in the downstream direction is not more than 180 degrees according to the estimated fluid flow value;
a measurement side tone frequency and calculation side tone frequency determination module for selecting a plurality of sets of measurement side tone frequencies and calculation side tone frequencies based on the minimum side tone frequency, wherein each set of calculation side tone frequencies is equal to a difference between a maximum measurement side tone frequency and the remaining measurement side tone frequencies, and the minimum calculation side tone frequency does not exceed the minimum side tone frequency;
the actual phase difference value calculating module is used for calculating the actual phase difference value corresponding to each group of side tone frequency for measurement and the side tone frequency for calculation according to the measurement phase difference value under each group of side tone frequency for measurement;
each group of fluid flow solving modules respectively calculate the fluid flow corresponding to each group of side tone frequency according to each group of side tone frequency for measurement and the corresponding actual phase difference value;
and the final fluid flow calculating module is used for calculating the average of the fluid flows corresponding to the side tone frequencies for each group of measurement to obtain the final fluid flow.
In one embodiment, the module for determining the frequency of the side tone for measurement and the frequency of the side tone for calculation includes:
a minimum calculation sidetone frequency determination module for selecting a group of sidetone frequencies not exceeding the minimum sidetone frequency as a minimum calculation sidetone frequency;
a maximum measurement side tone frequency determination module for setting a preset maximum set side tone frequency larger than the minimum side tone frequency as a maximum measurement side tone frequency;
the progressive multiple determining module is used for determining the upper limit value of the progressive multiple between adjacent groups of side tone frequencies according to the absolute error of ultrasonic phase measurement;
the other groups of side tone frequency determination modules for calculation determine the other groups of side tone frequencies for calculation between the minimum side tone frequency for calculation and the maximum side tone frequency for measurement according to the upper limit value of the progressive multiple between adjacent groups of side tone frequencies;
and the other groups of measuring side tone frequency determining modules determine the other groups of measuring side tone frequencies according to the difference value between the maximum measuring side tone frequency and each calculating side tone frequency.
According to another aspect of the present invention, there is also provided a storage medium including a stored program, which, when executed, controls an apparatus in which the storage medium is located to perform the continuous wave flow measurement method of multi-frequency information fusion.
According to another aspect of the present invention, there is also provided an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor implementing the multi-frequency information fusion continuous wave flow measurement method when executing the program.
The invention has the following beneficial effects:
the invention selects a plurality of groups of side tone frequencies for measurement and side tone frequencies for calculation on the basis of the minimum side tone frequency when the difference value is not more than 180 degrees, and finally obtains the average value of each fluid flow as the final fluid flow measurement result through the actual phase difference value corresponding to the side tone frequencies for measurement and the side tone frequencies for calculation, the obtained measurement error is superior to the measurement error obtained by a single frequency, the measurement precision of the flow is improved, and simultaneously, because the number of the side tone frequency for measurement is not changed, a new side tone frequency point for measurement is not additionally added, and only the calculation amount is slightly increased in the data processing by increasing the side tone frequency for calculation, because the newly increased calculation amount belongs to a low-speed phase signal, the increase of the calculation amount is negligible, therefore, the invention can improve the measurement precision and ensure that the measurement sensitivity (response time) is not reduced.
In addition to the objects, features and advantages described above, other objects, features and advantages of the present invention are also provided. The present invention will be described in further detail below with reference to the accompanying drawings.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this application, illustrate embodiments of the invention and, together with the description, serve to explain the invention and not to limit the invention. In the drawings:
fig. 1 is a flow chart of a continuous wave flow measurement method of multi-frequency information fusion according to a preferred embodiment of the present invention.
Fig. 2 is a flow chart of the selection process of the side-tone frequency for measurement and the side-tone frequency for calculation according to the preferred embodiment of the present invention.
Fig. 3(a) is a schematic diagram of the error of the flow measurement corresponding to a single frequency of 4.0MHz in the preferred embodiment of the present invention.
Fig. 3(b) is a schematic diagram of the error of the flow measurement corresponding to a single frequency of 4.5MHz in the preferred embodiment of the present invention.
Fig. 3(c) is a schematic diagram of the error of the flow measurement corresponding to a single frequency of 4.9MHz in the preferred embodiment of the present invention.
Fig. 3(d) is a schematic diagram of the error of the flow measurement value corresponding to a single frequency of 5MHz according to the preferred embodiment of the present invention.
Fig. 4 is a schematic diagram of flow measurement error after fusion of multi-frequency measurement information according to a preferred embodiment of the present invention.
Fig. 5(a) is a schematic diagram of the error of the flow measurement corresponding to a single frequency of 1.0MHz according to another preferred embodiment of the present invention.
Fig. 5(b) is a schematic flow chart of the error of the flow measurement corresponding to a single frequency of 1.2MHz according to another preferred embodiment of the present invention.
Fig. 6 is a schematic diagram of flow measurement error after fusion of multi-frequency measurement information according to another preferred embodiment of the present invention.
Fig. 7 is a schematic diagram of a module of a continuous wave flow measuring device with multi-frequency information fusion according to a preferred embodiment of the invention.
Fig. 8 is a detailed schematic diagram of the measuring side-tone frequency and calculating side-tone frequency determination module in accordance with the preferred embodiment of the present invention.
Detailed Description
It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict. The present invention will be described in detail below with reference to the embodiments with reference to the attached drawings.
As shown in fig. 1, a method for measuring a continuous wave flow by multi-frequency information fusion includes the steps of:
s100, according to the estimated fluid flow, obtaining the minimum sidetone frequency when the difference between the phase change value of the ultrasonic wave propagating in the fluid along the countercurrent direction and the phase change value of the ultrasonic wave propagating in the fluid along the downstream direction is not more than 180 degrees;
s200, selecting a plurality of groups of measuring side tone frequencies and calculating side tone frequencies based on the minimum side tone frequency, wherein each group of calculating side tone frequencies is equal to the difference between the maximum measuring side tone frequency and the rest measuring side tone frequencies, and the minimum calculating side tone frequency does not exceed the minimum side tone frequency;
s300, solving actual phase difference values corresponding to the side tone frequencies for measurement and the side tone frequencies for calculation according to the measurement phase difference values under the side tone frequencies for measurement of each group;
s400, respectively calculating the fluid flow corresponding to each group of side tone frequencies for measurement according to each group of side tone frequencies for measurement and the corresponding actual phase difference value;
s500, averaging the fluid flow corresponding to the side tone frequency for each group of measurement to obtain the final fluid flow.
As shown in fig. 2, in a preferred embodiment, selecting a plurality of sets of side tone frequencies for measurement and side tone frequencies for calculation based on the minimum side tone frequency, wherein each set of side tone frequencies for calculation is equal to the difference between the maximum side tone frequency for measurement and the remaining side tone frequencies for measurement, and the minimum side tone frequency for calculation does not exceed the minimum side tone frequency specifically includes the steps;
s201, selecting a group of side tone frequencies not exceeding the minimum side tone frequency as the minimum side tone frequency for calculation;
s202, taking a preset maximum set side tone frequency larger than the minimum side tone frequency as a maximum side tone frequency for measurement;
s203, determining the upper limit value of progressive multiples between adjacent groups of lateral tone frequencies according to the absolute error of ultrasonic phase measurement;
s204, determining other side tone frequencies for calculation between the minimum side tone frequency for calculation and the maximum side tone frequency for measurement according to the upper limit value of the progressive multiple between adjacent groups of side tone frequencies;
and S205, determining the frequencies of the rest groups of measuring side tones according to the difference between the maximum frequency of the measuring side tone and the frequencies of all the calculating side tones.
In this embodiment, the upper limit value of the progressive multiple is only a limit on the maximum value of the progressive multiple, and other progressive multiples not exceeding the maximum value may be selected between adjacent side tone frequencies, except that the maximum value is used as the progressive multiple.
Optionally, the progressive multiple k between adjacent sets of sidetone frequencies satisfies:
where δ represents the absolute error of the ultrasonic phase measurement.
Optionally, calculating the fluid flow rate corresponding to each group of side-tone frequencies for measurement according to each group of side-tone frequencies for measurement and the corresponding actual phase difference value thereof specifically includes:
wherein, ViRepresenting the sidetone frequency f for measurementiCorresponding fluid flow, L is the length of the pipeline in the flowmeter, R is the radius of the pipeline, C is the propagation speed of the ultrasonic wave in the static fluid, phidiff(fi) For measuring the frequency f of the sidetoneiThe corresponding actual phase difference value.
Optionally, the averaging the fluid flows corresponding to the side tone frequencies for each group of measurement to obtain the final fluid flow V specifically includes:
where N is the number of groups of sidetone frequencies for measurement.
In the above embodiment, the actual phase difference values at all frequencies are obtained by a method of ambiguity resolution, the flow measurement values corresponding to different frequencies are obtained by a formula, and then the final flow measurement value is obtained by averaging, assuming that the phase measurement errors at different frequencies are the same and can be represented as δ, the flow measurement error in the above embodiment can be theoretically represented as:
the pipe length L of the flowmeter is 0.2m, the propagation speed C of the ultrasonic wave in the static fluid is 1500m/s, the pipe radius R is 0.004m, and the estimated fluid flow U is 10 m/s. Thus, the phase difference Φdiff(f) Absence of fuzzy number
Provided that f < 281 KHz. For the situation, if the maximum frequency of the ultrasonic wave is 5MHz, and the progressive multiple k between two adjacent groups of lateral sound frequencies is less than or equal to 5, the minimum lateral sound frequency for calculation is 0.1MHz, the maximum lateral sound frequency for measurement is 5MHz, and the maximum lateral sound frequency for measurement is selected before passing throughBy the method, the rest of the side tone frequencies for calculation are 0.5MHz and 1MHz, and then the side tone frequencies for calculation are respectively subtracted by utilizing the maximum side tone frequency for measurement of 5MHz, so that all the side tone frequencies for measurement can be obtained: 5MHz, 4.9MHz, 4MHz, 4.5 MHz.
It will be understood by those skilled in the art that the actual phase difference corresponding to the high frequency sound wave can be obtained by the following formula:
wherein [ ·]0.5Characterization of rounding off,. phifrac(f1) For the side tone frequency of f1Number of measurable phases in one cycle, phidiff(f2) For the side tone frequency of f2Actual phase difference of time (no ambiguity number), Φdiff(f1) For the side tone frequency of f1The actual phase difference (no ambiguity number) of time, k is the progressive multiple between two adjacent sets of side tone frequencies.
For the value of k, there is an assumption that the phase fluctuation does not bring about a change in the value of the above formula, i.e., satisfies
For the ultrasonic probe, it is difficult to achieve the frequency bandwidth distribution required for the measurement. For this reason, the processing needs to be performed by using a folding sound method. The third frequency f can be obtained from equation (7)3＝f1-f2Phase difference of (2):
with selective side frequency f3＝f1-f2So that the phase difference phi is obtained at the frequencydiff(f3) Within a period, the problem of fuzzy of the whole period, namely the formula (7)Can be expressed as:
Φdiff(f3)＝Φfrac(f1)-Φfrac(f2). (8)
in combination with equation (5), the phase difference without the whole period ambiguity corresponding to the high frequency sound wave can be obtained, as in the present embodiment, f1＝5MHz；f2＝4.9MHz；f3＝f1-f2The phase difference without ambiguity number at 0.1MHz is used for 0.1MHz at 5MHz to 4.9MHz, and the phase difference without ambiguity number corresponding to different high frequency signals can be obtained by the formulas (5) and (8) because the phase difference without ambiguity number at 0.1MHz is used for 0.1MHz, and the phase difference in one period obtained by measurement at each measurement side tone frequency is given in table 1.
TABLE 1 measured (with full cycle ambiguity) phase differences
The actual (full cycle ambiguity free) phase difference and flow values in the above table can be obtained using equations (5) and (3), see table 2.
TABLE 2 actual (without full cycle ambiguity) phase difference and flow values
| Side tone frequency f (MHz) | Φdiff(°) | Flow rate value: ml/s (formula 3) |
| 0.1 | 64 | 139.62 |
| 0.5 | 320 | 139.62 |
| 1 | 640 | 139.62 |
| 4 | 2560 | 139.62 |
| 4.5 | 2880 | 139.62 |
| 4.9 | 3136 | 139.62 |
| 5 | 3200 | 139.62 |
The actual (without full cycle ambiguity) phase differences in table 2 are specifically calculated as follows:
for the case of 0.1MHz, since it is less than 281KHz, there is no ambiguity in the corresponding phase difference, and 0.1MHz is 5MHz to 4.9MHz, so that the following formula (8) can be obtained:
Φdiff(0.1MHz)＝Φfrac(5MHz)-Φfrac(4.9MHz)＝-40°-(-104°)＝64°；
further, for the actual phase difference of each subsequent sidetone frequency, the iterative calculation process is as follows:
in table 2 above, the phase difference measurement is assumed to be error free. Assuming that the phase difference fluctuates randomly within ± 1 ° at each frequency, the measured values fluctuate as shown in fig. 3(a) -3(d) at a single frequency.
After multi-frequency fusion, the obtained flow measurement error is shown in fig. 4, and table 3 lists the standard deviation corresponding to the corresponding measurement data:
TABLE 3 measurement of corresponding standard deviation
| Frequency (MHz) | 4 | 4.5 | 4.9 | 5 | After fusion |
| Standard deviation (ml/s) | 0.0321 | 0.0286 | 0.0253 | 0.0253 | 0.0134 |
It is obvious that the measurement error obtained from the information fusion is significantly better than that obtained from a single frequency, and at the same time, in the above embodiment, the number of the measurement side-tone frequencies is not additionally increased, and therefore, the measurement sensitivity (response time) is not reduced.
Similarly, the measurement method of the present invention can be applied to the patent "ultrasonic flow measurement method and apparatus based on lateral tone phase measurement", since the highest frequency in the patent is set to 1MHz, which is different from the above-mentioned setting condition, but the data fusion scheme of the present invention can also be adopted. Specifically, the sidetone frequencies are required to be 1MHz and 200KHz, and the sidetone frequencies for measurement are 1MHz and 1.2MHz, or 1MHz and 0.8MHz, by using the folded tone method. In the following analysis, only the cases of 1MHz and 1.2MHz are considered, and the analysis of the other case is similar and will not be described again, and table 4 gives the phase value in one period measured under the corresponding frequency.
TABLE 4 phase differences with full period ambiguity
The actual (full-cycle ambiguity free) phase difference and flow values in the above table can be obtained using equations (5) and (3), as shown in table 5:
TABLE 5 actual (without full cycle ambiguity) phase difference and flow values
| Side tone frequency f (MHz) | Φdiff(°) | Flow rate value: ml/s (formula 3) |
| 0.2 | 128 | 139.62 |
| 1 | 640 | 139.62 |
| 1.2 | 768 | 139.62 |
The specific calculation is as follows: for the case of 0.1MHz, since it is less than 281KHz, the corresponding phase difference has no ambiguity problem, and for the phase difference of the subsequent frequency, the iterative calculation process is as follows:
similar to the above-described conditions, assuming that the phase difference fluctuates randomly within ± 1 ° at each frequency, the measured values fluctuate as shown in fig. 5(a) -5(b) at a single frequency. The flow measurement error obtained after the multi-frequency fusion is performed is shown in fig. 6, table 6 lists the standard deviation corresponding to the corresponding measurement data, it can be seen that the measurement error obtained according to the information fusion is significantly better than the measurement error obtained by a single frequency, and meanwhile, in the above embodiment, the number of the side-tone frequencies for measurement is not additionally increased, and the measurement sensitivity (response time) is not reduced.
Table 6. measurement of corresponding standard deviation:
| frequency (MHz) | 1 | 1.2 | Fusion |
| Standard deviation (ml/s) | 0.1262 | 0.1042 | 0.0825 |
Optionally, the number N of groups of the measurement sidetone frequencies satisfies:
wherein, TmaxFor single measurement of maximum time constraint, UmaxMeasuring the maximum range of flow rate for the flowmeter, fmaxThe maximum frequency of the flowmeter, N is the group number of the side-tone frequencies for measurement, k is the progressive multiple between two adjacent groups of the side-tone frequencies, L is the length of a pipeline in the flowmeter, R represents the radius of the pipeline, and C is the ultrasonic wave in the static fluidThe propagation speed of (c). For data fusion, the higher the number of groups N of measurement side tone frequencies, the better the measurement accuracy. However, the number N of groups of the side tone frequencies for measurement is constrained by the measurement sensitivity (response time), and although the number is large, the measurement accuracy is better, but the measurement sensitivity (response time) is affected, so that the number must be constrained, the equation (17) can well realize the perfect balance between the accuracy improvement and the measurement sensitivity, and the related parameters need to be optimally set according to the measurement requirements in engineering application.
According to the corresponding formula:
(v) error in the same measured phase difference +diffThe higher the measurement frequency, the better the measurement accuracy, and this phenomenon can be compared with fig. 3(a) -3(d) and fig. 5(a) -5 (b). Therefore, in the case of information fusion, it is desirable that the higher the frequency corresponding to the fused phase, the better. For the previous embodiment, if the frequency of fusion is higher, such as 5MHz and 4.99MHz, the measurement accuracy is theoretically higher than the fusion of 5MHz and 4.9 MHz. In addition, the phase difference information of 0.01MHz can be obtained by measuring the phase difference signals of 5MHz and 4.99MHz, so that a wider measuring range can be obtained, and the phase difference information of 0.1MHz can be obtained by measuring the phase difference signals of 5MHz and 4.9 MHz. When the lowest frequency of 0.01MHz is adopted, the number of groups of the side tone frequency to be calculated is 4 according to the example that the frequency multiple is less than or equal to 5: 0.01MHz, 0.05MHz, 0.25MHz, 1.25MHz, the number of sets of side tone frequencies for measurement was 5: 5MHz, 4.99MHz, 4.95MHz, 4.75MHz, 3.75 MHz. In contrast, for the lowest frequency of 0.1MHz, the number of sets of side tone frequencies to be calculated is 3: 0.1MHz, 0.5MHz, 1MHz, the number of sets of side tone frequencies for measurement required is 4: 5MHz, 4.9MHz, 4.5MHz, 4 MHz. The more the frequency groups are, the longer the measurement time is, and the slower the response is, so that the optimization selection problem exists, which is specifically as follows:
constraint conditions are as follows: assuming a single measurement with a maximum time constraint of TmaxThe maximum range of the flow rate measured by the flowmeter is UmaxThe flowmeter is the mostHigh frequency of fmax. Then the minimum requirement for obtaining the side tone frequency is satisfied:
since each frequency requires one measurement of forward flow and reverse flow during the measurement process, the time taken for each frequency measurement is:
in the above formula, L/(C + U) and L/(C-U) represent the time taken for forward and reverse flow measurements, respectively. The approximation of equation (20) is based on the assumption that the flow velocity (generally measured flow velocity is 0-10 m/s) is much smaller than the sound velocity (sound wave propagation velocity in water at normal temperature is 1500m/s), and therefore, the number of usable measurement side-tone frequency groups satisfies:
assuming that the frequency multiples are consistent and set to k (1 < k ≦ 6 in engineering applications), the lowest frequency can be expressed as:
substituting equation (19) into the above equation yields:
meanwhile, the constraint condition shown in the formula (17) can be obtained by combining the formula (21), and as long as the group number N of the side tone frequencies for measurement meets the corresponding constraint, the response speed can not be slowed down while the accuracy of the measurement method is improved.
As shown in fig. 7, another embodiment of the present invention further provides a continuous wave flow measuring apparatus with multi-frequency information fusion, including:
the minimum side-tone frequency determining module is used for solving the minimum side-tone frequency when the difference value between the phase change value of the ultrasonic wave in the fluid in the upstream direction and the phase change value of the ultrasonic wave in the fluid in the downstream direction is not more than 180 degrees according to the estimated fluid flow value;
a measurement side tone frequency and calculation side tone frequency determination module for selecting a plurality of sets of measurement side tone frequencies and calculation side tone frequencies based on the minimum side tone frequency, wherein each set of calculation side tone frequencies is equal to a difference between a maximum measurement side tone frequency and the remaining measurement side tone frequencies, and the minimum calculation side tone frequency does not exceed the minimum side tone frequency;
the actual phase difference value calculating module is used for calculating the actual phase difference value corresponding to each group of side tone frequency for measurement and the side tone frequency for calculation according to the measurement phase difference value under each group of side tone frequency for measurement;
each group of fluid flow solving modules respectively calculate the fluid flow corresponding to each group of measuring side tone frequency according to each group of measuring side tone frequency and the corresponding actual phase difference value;
and the average fluid flow calculating module is used for calculating the average of the fluid flows corresponding to the side tone frequencies for each group of measurement to obtain the final fluid flow.
The ultrasonic flow measuring device based on the side-tone phase measurement of the invention obtains the minimum side-tone frequency that the difference value of the phase change value of ultrasonic wave propagating along the countercurrent direction and the phase change value propagating along the downstream direction in the fluid does not exceed 180 degrees according to the estimated value of the fluid flow, then selects a plurality of groups of side-tone frequencies for measurement and side-tone frequencies for calculation based on the minimum side-tone frequency, wherein, each group of side-tone frequencies for calculation are respectively equal to the difference between the maximum side-tone frequency for measurement and the rest side-tone frequencies for measurement, and the minimum side-tone frequency for calculation does not exceed the minimum side-tone frequency, then obtains the actual phase difference value corresponding to each group of side-tone frequencies for measurement and the side-tone frequencies for calculation according to the measurement phase difference value of the ultrasonic wave propagating along the countercurrent direction and the phase change value along the downstream direction in the fluid under each group of side-tone frequencies for measurement, and finally, averaging the fluid flow corresponding to each group of side-tone frequencies for measurement to obtain the final fluid flow, so that the technical problems that the data is not sufficiently utilized and the flow accuracy needs to be improved in the conventional ultrasonic flow measurement method are solved, and the response speed of the flowmeter cannot be reduced while a high-accuracy measurement result is obtained.
Optionally, as shown in fig. 8, the module for determining the frequency of the side tone for measurement and the frequency of the side tone for calculation includes:
a minimum calculation sidetone frequency determination module for selecting a group of sidetone frequencies not exceeding the minimum sidetone frequency as a minimum calculation sidetone frequency;
a maximum measurement side tone frequency determination module for setting a preset maximum set side tone frequency larger than the minimum side tone frequency as a maximum measurement side tone frequency;
the progressive multiple determining module is used for determining the upper limit value of the progressive multiple between adjacent groups of side tone frequencies according to the absolute error of ultrasonic phase measurement;
the other groups of side tone frequency determination modules for calculation determine the other groups of side tone frequencies for calculation between the minimum side tone frequency for calculation and the maximum side tone frequency for measurement according to the upper limit value of the progressive multiple between adjacent groups of side tone frequencies;
and the other groups of measuring side tone frequency determining modules determine the other groups of measuring side tone frequencies according to the difference value between the maximum measuring side tone frequency and each calculating side tone frequency.
Another embodiment of the present invention provides a storage medium including a stored program, wherein an apparatus in which the storage medium is located is controlled to perform the continuous wave flow measurement method as the multi-frequency information fusion when the program is executed.
Another embodiment of the present invention provides an electronic device, which includes a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor implements the multi-frequency information fusion continuous wave flow measurement method when executing the program.
It should be noted that the steps illustrated in the flowcharts of the figures may be performed in a computer system such as a set of computer-executable instructions and that, although a logical order is illustrated in the flowcharts, in some cases, the steps illustrated or described may be performed in an order different than presented herein.
The functions of the method of the present embodiment, if implemented in the form of software functional units and sold or used as independent products, may be stored in one or more storage media readable by a computing device. Based on such understanding, part of the contribution of the embodiments of the present invention to the prior art or part of the technical solution may be embodied in the form of a software product stored in a storage medium and including instructions for causing a computing device (which may be a personal computer, a server, a mobile computing device, a network device, or the like) to execute all or part of the steps of the method described in the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
The embodiments are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same or similar parts among the embodiments are referred to each other.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (8)
1. A multi-frequency information fusion continuous wave flow measuring method is characterized by comprising the following steps:
according to the estimated fluid flow, the minimum lateral sound frequency is obtained when the difference between the phase change value of the ultrasonic wave propagating in the fluid along the countercurrent direction and the phase change value of the ultrasonic wave propagating in the fluid along the downstream direction is not more than 180 degrees;
selecting a plurality of groups of measuring side tone frequencies and calculating side tone frequencies based on the minimum side tone frequency, wherein each group of calculating side tone frequencies is equal to the difference between the maximum measuring side tone frequency and the rest measuring side tone frequencies, and the minimum calculating side tone frequency does not exceed the minimum side tone frequency;
according to the measured phase difference value under each group of the measured side tone frequency, the actual phase difference value corresponding to each group of the measured side tone frequency and the calculated side tone frequency is obtained;
respectively calculating the fluid flow corresponding to each group of side-tone frequencies for measurement according to each group of side-tone frequencies for measurement and the corresponding actual phase difference value;
averaging the fluid flow corresponding to each group of measurement side-tone frequency to obtain the final fluid flow;
selecting a plurality of sets of side tone frequencies for measurement and side tone frequencies for calculation based on the minimum side tone frequency, wherein each set of side tone frequencies for calculation is equal to a difference between the maximum side tone frequency for measurement and the remaining side tone frequencies for measurement, and the minimum side tone frequency for calculation does not exceed the minimum side tone frequency;
selecting a group of sidetone frequencies not exceeding the minimum sidetone frequency as a minimum side tone frequency for calculation;
setting a preset maximum set sidetone frequency greater than the minimum sidetone frequency as a maximum measurement sidetone frequency;
determining the upper limit value of progressive multiples between adjacent groups of side tone frequencies according to the absolute error of ultrasonic phase measurement;
determining the other groups of side tone frequencies for calculation between the minimum side tone frequency for calculation and the maximum side tone frequency for measurement according to the upper limit value of the progressive multiple between adjacent groups of side tone frequencies;
and determining the frequencies of the rest groups of side tones for measurement according to the difference between the frequency of the maximum side tone for measurement and the frequencies of all the side tones for calculation.
2. The method of claim 1, wherein the progressive multiple k between adjacent groups of sidetone frequencies satisfies:
where δ represents the absolute error of the ultrasonic phase measurement.
3. The method of claim 2, wherein the multi-frequency information fusion continuous wave flow measurement is performed by a plurality of frequency sensors,
the number N of groups of the measurement sidetone frequencies satisfies:
wherein, TmaxFor single measurement of maximum time constraint, UmaxMeasuring the maximum range of flow rate for the flowmeter, fmaxThe maximum frequency of the flowmeter, L is the length of a pipeline in the flowmeter, and C is the propagation speed of the ultrasonic wave in the static fluid.
4. The method of claim 3, wherein the calculating the fluid flow rate corresponding to each set of measuring side-tone frequencies according to each set of measuring side-tone frequencies and the corresponding actual phase difference values comprises:
wherein, ViRepresenting the sidetone frequency f for measurementiCorresponding fluid flow, L is the length of the pipeline in the flowmeter, R is the radius of the pipeline, C is the propagation speed of the ultrasonic wave in the static fluid, phidiff(fi) For measuring side audioRate fiThe corresponding actual phase difference value.
5. The method of claim 4, wherein the averaging of the fluid flow rates corresponding to the side-tone frequencies for each measurement set to obtain the final fluid flow rate V is specifically:
6. a multi-frequency information fused continuous wave flow measuring device is characterized by comprising:
the minimum side-tone frequency determining module is used for solving the minimum side-tone frequency when the difference value between the phase change value of the ultrasonic wave in the fluid in the upstream direction and the phase change value of the ultrasonic wave in the fluid in the downstream direction is not more than 180 degrees according to the estimated fluid flow value;
a measurement side tone frequency and calculation side tone frequency determination module for selecting a plurality of sets of measurement side tone frequencies and calculation side tone frequencies based on the minimum side tone frequency, wherein each set of calculation side tone frequencies is equal to a difference between a maximum measurement side tone frequency and the remaining measurement side tone frequencies, and the minimum calculation side tone frequency does not exceed the minimum side tone frequency;
the actual phase difference value calculating module is used for calculating the actual phase difference value corresponding to each group of side tone frequency for measurement and the side tone frequency for calculation according to the measurement phase difference value under each group of side tone frequency for measurement;
each group of fluid flow solving modules respectively calculate the fluid flow corresponding to each group of measuring side tone frequency according to each group of measuring side tone frequency and the corresponding actual phase difference value;
the average fluid flow calculating module is used for calculating the average of the fluid flows corresponding to the side tone frequencies for each group of measurement to obtain the final fluid flow;
the measurement side tone frequency and calculation side tone frequency determination module includes:
a minimum calculation sidetone frequency determination module for selecting a group of sidetone frequencies not exceeding the minimum sidetone frequency as a minimum calculation sidetone frequency;
a maximum measurement side tone frequency determination module for setting a preset maximum set side tone frequency larger than the minimum side tone frequency as a maximum measurement side tone frequency;
the progressive multiple determining module is used for determining the upper limit value of the progressive multiple between adjacent groups of side tone frequencies according to the absolute error of ultrasonic phase measurement;
the other groups of side tone frequency determination modules for calculation determine the other groups of side tone frequencies for calculation between the minimum side tone frequency for calculation and the maximum side tone frequency for measurement according to progressive multiples between adjacent groups of side tone frequencies;
and the other groups of measuring side tone frequency determining modules determine the other groups of measuring side tone frequencies according to the difference value between the maximum measuring side tone frequency and each calculating side tone frequency.
7. A storage medium comprising a stored program, characterized in that, when the program is executed, the storage medium is controlled by an apparatus in which the storage medium is located to perform the multi-frequency information fusion continuous wave flow measurement method according to any one of claims 1 to 5.
8. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor implements the multi-frequency information fusion continuous wave flow measurement method according to any one of claims 1 to 5 when executing the program.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911403014.7A CN111157064B (en) | 2019-12-31 | 2019-12-31 | Multi-frequency information fusion continuous wave flow measuring method and device and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911403014.7A CN111157064B (en) | 2019-12-31 | 2019-12-31 | Multi-frequency information fusion continuous wave flow measuring method and device and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111157064A CN111157064A (en) | 2020-05-15 |
| CN111157064B true CN111157064B (en) | 2021-06-25 |
Family
ID=70559497
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911403014.7A Active CN111157064B (en) | 2019-12-31 | 2019-12-31 | Multi-frequency information fusion continuous wave flow measuring method and device and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111157064B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111829601B (en) * | 2020-07-21 | 2022-03-18 | 中国人民解放军国防科技大学 | Multi-state parameter synchronous measurement method and device of fluid, electronic equipment and medium |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS63101713A (en) * | 1986-10-17 | 1988-05-06 | Tokyo Keiki Co Ltd | Continuous wave phase difference type ultrasonic flowmeter |
| CN101672693A (en) * | 2009-09-30 | 2010-03-17 | 西华大学 | Non-fuzziness measuring method on ultrasonic transmission time |
| CN102639970A (en) * | 2009-11-30 | 2012-08-15 | 独立行政法人产业技术综合研究所 | Flow rate measuring device |
| CN105115552A (en) * | 2015-04-30 | 2015-12-02 | 中国人民解放军国防科学技术大学 | Ultrasonic wave flow measuring method and device based on sidetone phase measurement |
| RU2597666C1 (en) * | 2015-07-31 | 2016-09-20 | Федеральное государственное бюджетное учреждение науки Институт проблем управления им. В.А. Трапезникова Российской академии наук | Method of measuring mass flow rate of liquid media |
| CN107290564A (en) * | 2017-07-28 | 2017-10-24 | 电子科技大学 | A kind of ultrasonic flow-velocity measuring method based on phase difference |
| CN107741382A (en) * | 2017-08-29 | 2018-02-27 | 中国人民解放军国防科技大学 | Liquid density measurement method and system based on continuous sound wave propagation |
-
2019
- 2019-12-31 CN CN201911403014.7A patent/CN111157064B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS63101713A (en) * | 1986-10-17 | 1988-05-06 | Tokyo Keiki Co Ltd | Continuous wave phase difference type ultrasonic flowmeter |
| CN101672693A (en) * | 2009-09-30 | 2010-03-17 | 西华大学 | Non-fuzziness measuring method on ultrasonic transmission time |
| CN102639970A (en) * | 2009-11-30 | 2012-08-15 | 独立行政法人产业技术综合研究所 | Flow rate measuring device |
| CN105115552A (en) * | 2015-04-30 | 2015-12-02 | 中国人民解放军国防科学技术大学 | Ultrasonic wave flow measuring method and device based on sidetone phase measurement |
| RU2597666C1 (en) * | 2015-07-31 | 2016-09-20 | Федеральное государственное бюджетное учреждение науки Институт проблем управления им. В.А. Трапезникова Российской академии наук | Method of measuring mass flow rate of liquid media |
| CN107290564A (en) * | 2017-07-28 | 2017-10-24 | 电子科技大学 | A kind of ultrasonic flow-velocity measuring method based on phase difference |
| CN107741382A (en) * | 2017-08-29 | 2018-02-27 | 中国人民解放军国防科技大学 | Liquid density measurement method and system based on continuous sound wave propagation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111157064A (en) | 2020-05-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7581453B2 (en) | Ultrasonic flow meter system | |
| AU2009302284B2 (en) | Viscous fluid flow measurement using a differential pressure measurement and a sonar measured velocity | |
| CN111157064B (en) | Multi-frequency information fusion continuous wave flow measuring method and device and electronic equipment | |
| EP2275787A2 (en) | Coriolis flowmeter | |
| EP0017475A1 (en) | Acoustic flowmeter with Reynolds number compensation | |
| US11549841B2 (en) | Ultrasonic meter employing two or more dissimilar chordal multipath integration methods in one body | |
| CN105115552B (en) | The ultrasonic flow measuring method and device of phase are surveyed based on sidetone | |
| CN108431554B (en) | Fluid measuring device | |
| EP3370074B1 (en) | Method for detecting background noise of sensor, and device thereof | |
| EP2423650A2 (en) | Coriolis mass flowmeter | |
| CN103649691A (en) | Method for determining the flow rate of fluids using the ultrasonic transit-time method | |
| Willatzen | Ultrasonic flowmeters: temperature gradients and transducer geometry effects | |
| CN115854945A (en) | Method and device for determining sediment deposition thickness | |
| Willatzen | Sound propagation in a moving fluid confined by cylindrical walls—a comparison between an exact analysis and the local-plane-wave approximation | |
| Ma et al. | Signal processing method based on connection fitting of echo peak point with a large slope for ultrasonic gas flow meter | |
| US3097526A (en) | fischbacher | |
| JP2001249039A (en) | Ultrasonic gas flow-velocity measuring method | |
| JPH0454892B2 (en) | ||
| CN113687345B (en) | Doppler radar water flow velocity measurement method and device | |
| JP7174574B2 (en) | Ultrasonic flowmeter and determination method of zero-crossing time in ultrasonic flowmeter | |
| JP4548485B2 (en) | Pulse doppler type ultrasonic flow meter, its program | |
| CN105890684A (en) | Novel setting method for determining sound channel positions by adopting Gauss-Jacobi polynomial | |
| Tougher et al. | On the approach to the critical solution in leading order thin-film coating and rimming flow | |
| CN117740090A (en) | Ultrasonic clutter elimination method, system and flowmeter | |
| CN115900850A (en) | Detection method of electromagnetic flowmeter and electromagnetic flowmeter |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |