Ultrasonic flow velocity analysis method and computer readable storage medium

Abstract

The present invention belongs to the field of analysis and detection technology, and relates to an ultrasonic flow velocity analysis method and a computer-readable storage medium, specifically comprising: an upstream ultrasonic transducer and a downstream ultrasonic transducer alternately transmit ultrasonic waves at a preset frequency, collect ultrasonic signals received by the downstream ultrasonic transducer and the upstream ultrasonic transducer after the ultrasonic waves pass through the flue gas in real time, and also collect the flue gas temperature and pressure in the flue gas; respectively process the ultrasonic signals received by the upstream ultrasonic transducer and the downstream ultrasonic transducer to obtain the calculation of the upstream target time and the downstream target time; calculate the flue gas flow velocity according to the upstream target time, the downstream target time, the distance between the upstream ultrasonic transducer and the downstream ultrasonic transducer, and the angle between the connecting line of the upstream ultrasonic transducer and the downstream ultrasonic transducer and the flue gas radial direction. The ultrasonic flow velocity analysis method of the present invention effectively improves the measurement accuracy and stability of the ultrasonic flowmeter.

Description

Technical Field

[0001] The present invention belongs to the technical field of analysis and detection, and in particular relates to an ultrasonic flow velocity analysis method and a computer-readable storage medium. Background Art

[0002] As the requirements for the measurement accuracy of flue gas velocity are getting higher and higher, the traditional Pitot tube method based on the differential pressure principle has problems such as single-point monitoring cannot reflect the average velocity of the pipeline cross section, is not suitable for occasions with high humidity, high dust and corrosive flue gas, and inaccurate low flow rate measurement. New technical methods are needed to replace them, and ultrasonic flow velocity monitoring technology has become an ideal alternative.

[0003] The traditional ultrasonic flow velocity is based on the principle of flight time difference. A pair of ultrasonic transducers are installed on the upstream and downstream sides of the pipeline, facing each other; the upstream transducer transmits ultrasonic waves, and the downstream transducer receives ultrasonic signals to obtain the downstream flight time; the downstream transducer transmits ultrasonic waves, and the upstream transducer receives ultrasonic signals to obtain the upstream flight time; the downstream flight time and the upstream flight time contain the flue gas flow rate information, and the flow rate can be calculated; the accuracy of the flight time directly determines the accuracy of the flow rate measurement. At present, the flight time extraction methods include the threshold method and the cross-correlation method; however, in actual application sites, due to factors such as temperature changes, aging or contamination of core drive devices (transducers, high-voltage relays, etc.), and rapid changes in environmental pressure fields or flow fields, performance degradation or deterioration may occur, as follows:

[0004] (1) Ultrasonic signals are aged or contaminated by devices. Although there is a post-amplification circuit, the noise is also amplified accordingly while the post-amplification circuit amplifies the signal. The signal-to-noise ratio is not improved, and the analysis accuracy is reduced.

[0005] (2) Temperature changes, environmental pressure fields, or rapid changes in flow fields cause the echo to be delayed or to exceed the signal analysis area in advance, resulting in the inability to analyze effective signals and poor measurement accuracy;

[0006] (3) The phenomenon of wave jumping occurs, that is, the effective signal is shifted forward or backward by one or more cycles, resulting in periodic jumps in the flow velocity measured in different cycles. Summary of the invention

[0007] Based on the above-mentioned shortcomings and deficiencies in the prior art, one of the objects of the present invention is to at least solve one or more of the above-mentioned problems in the prior art. In other words, one of the objects of the present invention is to provide an ultrasonic flow velocity analysis method and a computer-readable storage medium that meet one or more of the above-mentioned needs.

[0008] In order to achieve the above-mentioned object of the invention, the present invention adopts the following technical solutions:

[0009] An ultrasonic flow velocity analysis method is provided, wherein an upstream ultrasonic transducer and a downstream ultrasonic transducer are respectively installed on both sides of the flue, and the upstream ultrasonic transducer and the downstream ultrasonic transducer are mutually irradiated. The ultrasonic flow velocity analysis method comprises the following steps:

[0010] S1, the upstream ultrasonic transducer and the downstream ultrasonic transducer emit ultrasonic waves alternately at a preset frequency, and collect ultrasonic signals received by the downstream ultrasonic transducer and the upstream ultrasonic transducer after the ultrasonic waves pass through the flue gas in real time;

[0011] The flue gas temperature and pressure in the flue are also collected;

[0012] S2. Process the ultrasonic signals received by the upstream ultrasonic transducer and the downstream ultrasonic transducer as follows:

[0013] Extracting the effective analysis interval of the ultrasonic signal according to the flue gas temperature, pressure and the distance between the upstream ultrasonic transducer and the downstream ultrasonic transducer to obtain the target interval of the ultrasonic signal;

[0014] Performing Hilbert transform on the target interval of the ultrasonic signal, and extracting the real part and the imaginary part after the Hilbert transform;

[0015] Poll the real part after Hilbert transform to determine whether there are S consecutive peak time points whose corresponding real parts are all greater than the set threshold; if so, enter the normal measurement mode, and extract the S consecutive peak time points after the first time greater than the set threshold as the first target time set; wherein the real parts corresponding to the S consecutive peak time points are all greater than the set threshold, and S is an integer greater than 2;

[0016] An inverse tangent function is constructed based on the real part and the imaginary part after the Hilbert transform, and Q consecutive zero-crossing time points with a value of zero for the inverse tangent function are extracted as a second target time set; wherein the interval between adjacent zero-crossing time points is within a preset target range, and Q is an integer greater than 3;

[0017] Extracting the minimum value of the intersection of the first target time set and the second target time set as the target time, so as to obtain the upstream target time and the downstream target time respectively;

[0018] S3. Calculate the flue gas flow rate according to the upstream target time, the downstream target time, the distance between the upstream ultrasonic transducer and the downstream ultrasonic transducer, and the angle between the line connecting the upstream ultrasonic transducer and the downstream ultrasonic transducer and the flue radial direction.

[0019] As a preferred solution, in step S2, the effective analysis interval of the ultrasonic signal is:

[0020] ;

[0021] Where, F is the ultrasonic frequency of the ultrasonic transducer, T is the flue gas temperature, P is the flue gas pressure, L is the distance between the upstream ultrasonic transducer and the downstream ultrasonic transducer, C 0 is the speed of sound at 273.15K.

[0022] As a preferred solution, in step S2, for the target interval of the ultrasonic signal received by the upstream ultrasonic transducer , and its corresponding real and imaginary parts after Hilbert transformation are and ;

[0023] Among them, the Hilbert transform H is:

[0024] ;

[0025] For the target interval of the ultrasonic signal received by the downstream ultrasonic transducer , and its corresponding real and imaginary parts after Hilbert transformation are and ;

[0026] Among them, the Hilbert transform H is:

[0027] .

[0028] As a preferred solution, in step S2, constructing an inverse tangent function based on the real part and the imaginary part after the Hilbert transform includes:

[0029] Based on the real part and the imaginary part The inverse tangent function constructed is ;

[0030] Based on the real part and the imaginary part The inverse tangent function constructed is .

[0031] As a preferred solution, in step S2, adjacent zero-crossing time points t i+1 With t i The following conditions are met:

[0032] ; Where i∈[1,Q-1], f is the sampling frequency.

[0033] As a preferred solution, in step S5, the calculation formula for the flue gas flow rate is:

[0034] ;

[0035] Where, L is the distance between the upstream ultrasonic transducer and the downstream ultrasonic transducer, θ is the angle between the line connecting the upstream ultrasonic transducer and the downstream ultrasonic transducer and the flue radial direction, t r and t v are the upstream target time and the downstream target time respectively.

[0036] As a preferred solution, the step S2 further includes:

[0037] If there are no S consecutive wave peak time points whose corresponding real parts are all greater than the set threshold, the high-frequency pulse width vibration maintenance mode is started;

[0038] Among them, the pulse width of the ultrasonic wave emitted in the high-frequency pulse width vibration maintenance mode is 2 to 3 times the pulse width of the ultrasonic wave emitted in the normal measurement mode.

[0039] As a preferred solution, after starting the high-frequency pulse width vibration maintenance mode, repeat steps S1 to S2 to determine whether there are S consecutive peak time points whose real parts corresponding to them are all greater than the set threshold; if so, enter the normal measurement mode; if not, enter the excitation gain amplification mode to adjust the excitation gain of each ultrasonic transducer.

[0040] As a preferred solution, the adjustment process of the excitation gain includes the following steps:

[0041] (1) According to the highest peak amplitude of the real part after Hilbert transform Calculate the gain amplification factor β:

[0042] ; Among them, V max is the peak intensity of the ultrasonic signal at the maximum excitation power of the ultrasonic transducer;

[0043] (2) Amplify the excitation gain of the ultrasonic transducer according to the gain amplification factor β, then repeat steps S1 to S2 to poll the real part after the Hilbert transform, obtain the highest peak amplitude of the real part after the Hilbert transform, and judge whether the highest peak amplitude meets the target condition; if so, enter the normal measurement mode; if not, loop to step (1) and use the highest peak amplitude Replace the highest peak amplitude Recalculate the gain amplification factor;

[0044] The target conditions are: ;in, is the highest peak amplitude of the real part after Hilbert transform corresponding to the jth amplification of the excitation gain of the ultrasonic transducer, j∈[1,N], N is the number of times the excitation gain is amplified.

[0045] The present invention also provides a computer-readable storage medium, wherein instructions are stored in the computer-readable storage medium. When the instructions are executed on a computer, the computer executes the ultrasonic flow velocity analysis method as described in any of the above schemes.

[0046] Compared with the prior art, the present invention has the following beneficial effects:

[0047] (1) The ultrasonic flow velocity analysis method of the present invention effectively improves the measurement accuracy and stability of the ultrasonic flow meter;

[0048] (2) The ultrasonic flow velocity analysis method of the present invention uses high-frequency pulse width vibration maintenance to eliminate the attenuation of ultrasonic signals emitted by the ultrasonic transducer caused by dirt attached to the ultrasonic transducer, and uses a step-by-step amplitude excitation gain amplification mode to prevent signal mutations and the occurrence of wrong waves, while increasing the effective use period of the equipment;

[0049] (3) The computer-readable storage medium of the present invention can realize intelligent ultrasonic flow velocity analysis. BRIEF DESCRIPTION OF THE DRAWINGS

[0050] Figure 1 is a schematic diagram of the layout structure of the ultrasonic transducer in the flue of Example 1 of the present invention;

[0051] Figure 2 is a flow chart of an ultrasonic flow velocity analysis method according to Embodiment 1 of the present invention;

[0052] Figure 3 is a graph of the real part after Hilbert transformation of Example 1 of the present invention;

[0053] Figure 4 is a graph of the inverse tangent function of Example 1 of the present invention;

[0054] Figure 5 It is a graph showing the change of different measured flow velocities over time when the wave skipping phenomenon occurs in the prior art;

[0055] Figure 6 is a graph showing the change of different measured flow velocities over time obtained by the ultrasonic flow velocity analysis method of Example 1 of the present invention;

[0056] Figure 7 It is a flow chart of the ultrasonic flow velocity analysis method of Example 1 of the present invention after starting the high-frequency pulse width vibration maintenance mode. DETAILED DESCRIPTION

[0057] In order to more clearly illustrate the embodiments of the present invention, the specific implementation methods of the present invention will be described below with reference to the accompanying drawings. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those of ordinary skill in the art, other accompanying drawings and other implementation methods can be obtained based on these accompanying drawings without creative work.

[0058] Embodiment 1:

[0059] like Figure 1 As shown, ultrasonic transducers are installed on the upstream and downstream sides of the flue I, namely, the upstream ultrasonic transducer II and the downstream ultrasonic transducer III. The upstream ultrasonic transducer II and the downstream ultrasonic transducer III are mutually reflected. The distance between the upstream ultrasonic transducer II and the downstream ultrasonic transducer III is L, and the angle between the line connecting the upstream ultrasonic transducer and the downstream ultrasonic transducer and the radial direction of the flue is θ.

[0060] like Figure 2 As shown, the ultrasonic flow velocity analysis method of this embodiment includes the following steps:

[0061] (1) The upstream ultrasonic transducer and the downstream ultrasonic transducer emit ultrasonic waves alternately at a preset frequency F, and collect in real time the ultrasonic signals (hereinafter referred to as echo signals) received by the downstream ultrasonic transducer and the upstream ultrasonic transducer after the ultrasonic waves pass through the flue gas, and obtain the echo signal X(t) received by the upstream ultrasonic transducer and the echo signal Y(t) received by the downstream ultrasonic transducer, respectively, where t is time;

[0062] In addition, the flue gas temperature T and pressure P are collected through sensors installed in the flue;

[0063] Specifically, the performance parameters of the upstream ultrasonic transducer and the downstream ultrasonic transducer of this embodiment are the same, the power is measured by the highest peak amplitude of the echo signal, the transducer excitation gain is set to level 4, and the peak intensity of the echo signal under the maximum excitation power of the transducer is set to V max ; Usually 50%V max The echo signal strength is used as the factory excitation power of the transducer.

[0064] (2) The echo signal X(t) received by the upstream ultrasonic transducer and the echo signal Y(t) received by the downstream ultrasonic transducer are processed as follows:

[0065] Firstly, the effective analysis interval of the echo signal is extracted according to the flue gas temperature T, pressure P and the distance L between the upstream ultrasonic transducer and the downstream ultrasonic transducer to obtain the target interval of the echo signal;

[0066] Specifically, the effective analysis interval of the echo signal in this embodiment is:

[0067] ;

[0068] Wherein, F is the ultrasonic frequency of the ultrasonic transducer, that is, the above-mentioned preset frequency; C 0 is the speed of sound at 273.15K.

[0069] By extracting the effective analysis interval of the echo signal, the target interval segments of the echo signal X(t) are obtained respectively. And the target interval of the echo signal Y(t) , t * is the time within the target interval.

[0070] Next, the target interval segment of the echo signal is subjected to Hilbert transformation, and the real part and the imaginary part after the Hilbert transformation are extracted;

[0071] Specifically, for the target interval of the echo signal X(t) received by the upstream ultrasonic transducer , and its corresponding real and imaginary parts after Hilbert transformation are and ;

[0072] Among them, the Hilbert transform H is:

[0073] ;

[0074] For the target interval of the echo signal Y(t) received by the downstream ultrasonic transducer , and its corresponding real and imaginary parts after Hilbert transformation are and ;

[0075] Among them, the Hilbert transform H is:

[0076] ;

[0077] After obtaining the real part and imaginary part after the above Hilbert transform, poll the real part after the Hilbert transform to determine whether there are three consecutive peak time points whose real parts R are all greater than the set threshold k; if not, start the high-frequency pulse width vibration maintenance mode; if so, enter the normal measurement mode, extract the three consecutive peak time points after the first time greater than the set threshold as the first target time set; among which, the real parts corresponding to the three consecutive peak time points are all greater than the set threshold; if Figure 3 As shown, the time points corresponding to the three consecutive peaks in the dotted box can be obtained, that is, the three consecutive peak time points, and the real parts R corresponding to the peak time points are all greater than the set threshold k. Among them, the number of consecutive peak time points is not limited to 3, and the specific number can be adjusted according to actual application needs.

[0078] In addition, this embodiment also constructs an inverse tangent function Z based on the real part and the imaginary part after the Hilbert transform, and extracts 8 consecutive zero-crossing time points with an inverse tangent function value of zero as the second target time set; wherein the interval between adjacent zero-crossing time points is within a preset target range; the number of the above zero-crossing time points is not limited to 8, and the specific number can be adjusted according to actual application requirements.

[0079] This embodiment constructs an inverse tangent function based on the real part and the imaginary part after the Hilbert transform, including:

[0080] Based on the real part and the imaginary part The inverse tangent function constructed is ;

[0081] Based on the real part and the imaginary part The inverse tangent function constructed is ;

[0082] In addition, the interval between adjacent zero-crossing time points in this embodiment is within the preset target range, specifically, the interval between adjacent zero-crossing time points t i+1 With t i The following conditions are met:

[0083] ; Where i∈[1,7] is an integer, and f is the sampling frequency;

[0084] like Figure 4 As shown, eight consecutive zero-crossing time points within the dashed box can be obtained.

[0085] In addition, the above process of extracting 8 consecutive zero-crossing time points with an inverse tangent function value of zero as the second target time set specifically includes: extracting the zero-crossing time points of the inverse tangent function and forming an array, polling the zero-crossing time points in the array, and determining whether there are 8 consecutive zero-crossing time points that satisfy the interval between adjacent zero-crossing time points within a preset target range; if so, then 8 consecutive zero-crossing time points are used to obtain the second target time set; if not, starting the high-frequency pulse width vibration maintenance mode.

[0086] Finally, the minimum value of the intersection of the first target time set and the second target time set is extracted as the target time;

[0087] Therefore, according to the above processing of the echo signal X(t) and the echo signal Y(t), the upstream target time t is obtained accordingly. r and downstream target time t v .

[0088] (3) According to the upstream target time t r , downstream target time t vThe flue gas velocity is calculated based on the distance L between the upstream ultrasonic transducer and the downstream ultrasonic transducer and the angle θ between the line connecting the upstream ultrasonic transducer and the downstream ultrasonic transducer and the flue radial direction.

[0089] Specifically, the calculation formula of the flue gas flow rate V in this embodiment is:

[0090] .

[0091] like Figure 5 , a wave-jump phenomenon occurs, resulting in periodic jumps in the flow velocity measured in different periods; Figure 6 As shown, after the ultrasonic flow velocity analysis method of this embodiment is used to analyze the flue gas flow velocity, there is no periodic jump in the flow velocity measured in different periods, which effectively improves the measurement accuracy and stability of the ultrasonic flowmeter.

[0092] like Figure 7 As shown, the pulse width of the ultrasonic wave emitted in the high-frequency pulse width vibration maintenance mode of this embodiment is 2 to 3 times the pulse width of the ultrasonic wave emitted in the normal measurement mode. As an example, the value is 2 times, with an interval of 3s per cycle, for a total of 3 to 5 cycles. After starting the high-frequency pulse width vibration maintenance mode, single measurement and signal acquisition are entered, that is, the above steps (1) to (2) are repeated to determine whether the real parts corresponding to 3 consecutive peak time points are all greater than the set threshold; if so, the normal measurement mode is entered; if not, the excitation gain amplification mode is entered to adjust the excitation gain of each ultrasonic transducer to increase the signal strength of the emitted ultrasonic wave.

[0093] The adjustment process of the excitation gain in this embodiment specifically includes the following steps:

[0094] (a) According to the highest peak amplitude of the real part after the current Hilbert transform Calculate the gain amplification factor β:

[0095] ; Among them, V max is the peak intensity of the echo signal at the maximum excitation power of the ultrasonic transducer;

[0096] (b) amplifying the excitation gain of the ultrasonic transducer according to the gain amplification factor β, and then entering single measurement and signal acquisition, that is, repeating the above steps (1) to step (2) to poll the real part after the Hilbert transform, obtaining the highest peak amplitude of the real part after the Hilbert transform, and judging whether the highest peak amplitude of the real part after the Hilbert transform meets the target condition; if so, entering the normal measurement mode; if not, looping to step (a) and using the highest peak amplitude Replace the highest peak amplitude Recalculate the gain amplification factor;

[0097] The target conditions are: ;in, It is the maximum peak amplitude of the real part after Hilbert transform corresponding to the j-th amplification of the excitation gain of the ultrasonic transducer, j∈[1,N], N is the number of times the excitation gain is amplified; as an example, in this embodiment, N is taken as 3. When the maximum peak amplitude of the real part after Hilbert transform corresponding to the excitation gain of the ultrasonic transducer after 3 amplifications still does not meet the target condition, a fault is prompted.

[0098] The computer-readable storage medium of this embodiment stores instructions. When the instructions are executed on a computer, the computer executes the ultrasonic flow velocity analysis method to realize intelligent ultrasonic flow velocity analysis.

[0099] The above description is only a detailed description of the preferred embodiments and principles of the present invention. For ordinary technicians in this field, according to the ideas provided by the present invention, there will be changes in the specific implementation methods, and these changes should also be regarded as the protection scope of the present invention.

Claims

1. An ultrasonic flow velocity analysis method, wherein an upstream ultrasonic transducer and a downstream ultrasonic transducer are respectively installed on the upstream and downstream sides of the flue, and the upstream ultrasonic transducer and the downstream ultrasonic transducer are mutually reflected, characterized in that: The ultrasonic flow velocity analysis method comprises the following steps: S1, the upstream ultrasonic transducer and the downstream ultrasonic transducer emit ultrasonic waves alternately at a preset frequency, and collect ultrasonic signals received by the downstream ultrasonic transducer and the upstream ultrasonic transducer after the ultrasonic waves pass through the flue gas in real time; The flue gas temperature and pressure in the flue are also collected; S2. Process the ultrasonic signals received by the upstream ultrasonic transducer and the downstream ultrasonic transducer as follows: Extracting the effective analysis interval of the ultrasonic signal according to the flue gas temperature, pressure and the distance between the upstream ultrasonic transducer and the downstream ultrasonic transducer to obtain the target interval of the ultrasonic signal; Performing Hilbert transform on the target interval of the ultrasonic signal, and extracting the real part and the imaginary part after the Hilbert transform; Poll the real part after Hilbert transform to determine whether there are S consecutive peak time points whose corresponding real parts are all greater than the set threshold; if so, enter the normal measurement mode, and extract the S consecutive peak time points after the first time greater than the set threshold as the first target time set; wherein the real parts corresponding to the S consecutive peak time points are all greater than the set threshold, and S is an integer greater than 2; An inverse tangent function is constructed based on the real part and the imaginary part after the Hilbert transform, and Q consecutive zero-crossing time points with a value of zero for the inverse tangent function are extracted as a second target time set; wherein the interval between adjacent zero-crossing time points is within a preset target range, and Q is an integer greater than 3; Extracting the minimum value of the intersection of the first target time set and the second target time set as the target time, so as to obtain the upstream target time and the downstream target time respectively; S3, calculating the flue gas flow rate according to the upstream target time, the downstream target time, the distance between the upstream ultrasonic transducer and the downstream ultrasonic transducer, and the angle between the connecting line of the upstream ultrasonic transducer and the downstream ultrasonic transducer and the flue radial direction; In step S2, the effective analysis interval of the ultrasonic signal is: Wherein, F is the transmitting ultrasonic frequency of the ultrasonic transducer, T is the flue gas temperature, P is the flue gas pressure, L is the distance between the upstream ultrasonic transducer and the downstream ultrasonic transducer, and C0 is the speed of sound at 273.15K; The step S2 further comprises: If there are no S consecutive wave peak time points whose corresponding real parts are all greater than the set threshold, the high-frequency pulse width vibration maintenance mode is started; Among them, the pulse width of the ultrasonic wave emitted in the high-frequency pulse width vibration maintenance mode is 2 to 3 times the pulse width of the ultrasonic wave emitted in the normal measurement mode; After starting the high-frequency pulse width vibration maintenance mode, repeat steps S1 to S2 to determine whether there are P consecutive peak time points whose real parts are all greater than the set threshold; if so, enter the normal measurement mode; if not, enter the excitation gain amplification mode to adjust the excitation gain of each ultrasonic transducer.

2. The ultrasonic flow velocity analysis method according to claim 1, characterized in that: In step S2, for the target interval X(t * ), the corresponding real and imaginary parts after Hilbert transformation are R top =Re(H[X(t * )]) and I top =Im(H[X(t * )]); Among them, the Hilbert transform H is: For the target interval Y(t * ), the corresponding real and imaginary parts after Hilbert transformation are R down =Re(H[Y(t * )]) and I down =Im(H[Y(t * )]); Among them, the Hilbert transform H is:

3. The ultrasonic flow velocity analysis method according to claim 2, characterized in that: In step S2, constructing an inverse tangent function based on the real part and the imaginary part after the Hilbert transform includes: Based on the real part R top and the imaginary part I top The inverse tangent function constructed is Based on the real part R down and the imaginary part I down The inverse tangent function constructed is 4. The ultrasonic flow velocity analysis method according to claim 3, characterized in that: In step S2, adjacent zero-crossing time points t i+1 With t i The following conditions are met: Where i∈[1,Q-1] and f is the sampling frequency.

5. The ultrasonic flow velocity analysis method according to any one of claims 1 to 4, characterized in that: In step S3, the calculation formula of the flue gas flow rate is: Where, L is the distance between the upstream ultrasonic transducer and the downstream ultrasonic transducer, θ is the angle between the line connecting the upstream ultrasonic transducer and the downstream ultrasonic transducer and the flue radial direction, t r and t v are the upstream target time and the downstream target time respectively.

6. The ultrasonic flow velocity analysis method according to claim 1, characterized in that: The adjustment process of the excitation gain comprises the following steps: (1) According to the highest peak amplitude of the real part after Hilbert transform Calculate the gain amplification factor β: Among them, V max is the peak intensity of the ultrasonic signal at the maximum excitation power of the ultrasonic transducer; (2) Amplify the excitation gain of the ultrasonic transducer according to the gain amplification factor β, then repeat steps S1 to S2 to poll the real part after the Hilbert transform, obtain the highest peak amplitude of the real part after the Hilbert transform, and judge whether the highest peak amplitude meets the target condition; if so, enter the normal measurement mode; if not, loop to step (1) and use the highest peak amplitude Replace the highest peak amplitude Recalculate the gain amplification factor; The target conditions are: in, is the highest peak amplitude of the real part after Hilbert transform corresponding to the jth amplification of the excitation gain of the ultrasonic transducer, j∈[1,N], N is the number of times the excitation gain is amplified.

7. A computer-readable storage medium, wherein instructions are stored in the computer-readable storage medium, characterized in that: When the instructions are executed on a computer, the computer is enabled to execute the ultrasonic flow velocity analysis method as described in any one of claims 1 to 6.

Images