CN107655533A - A kind of Ultrasonic Wave Flowmeter signal processing method and system based on backward energy integration - Google Patents
A kind of Ultrasonic Wave Flowmeter signal processing method and system based on backward energy integration Download PDFInfo
- Publication number
- CN107655533A CN107655533A CN201711043275.3A CN201711043275A CN107655533A CN 107655533 A CN107655533 A CN 107655533A CN 201711043275 A CN201711043275 A CN 201711043275A CN 107655533 A CN107655533 A CN 107655533A
- Authority
- CN
- China
- Prior art keywords
- signal
- point
- area
- echo
- peak point
- 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.)
- Granted
Links
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 Ultrasonic Wave Flowmeter signal processing method based on backward energy integration is proposed, according to the ultrasonic echo signal of sampling, passes through the maximum positive peak point for comparing to obtain backward energy signal；Then section starting point is determined forward, with 0.1 times more than echo-signal normalization amplitude for Rule of judgment, finds integration starting point；From integration starting point to maximum positive peak point, a peak point is often found, will be normalized after the squared magnitude, calculated the trapezoidal area surrounded between adjacent peak point and add up；Judge the size of cumulative area value and preset area value；It is cumulative to stop integration after cumulative area value exceedes preset area value, preset area is subtracted by cumulative area value and is worth to area difference；The ordinate of known area difference and adjacent peak point, the intersection point of small trapezoidal upper bottom edge and X-axis is solved, as characteristic point；Following current, the adverse current propagation time of ultrasonic echo signal are determined, and calculates gas flow.
Description
Technical field
The present invention relates to field of flow measurement, is a kind of Ultrasonic Wave Flowmeter transmitter, and particularly one kind is based on back
The Ultrasonic Wave Flowmeter signal processing method and system of wave energy integration.Using gas ultrasound wave current amount proposed by the present invention
Signal processing method processing ultrasonic echo signal is counted, realizes the measurement of gas flow.
Background technology
Ultrasonic Wave Flowmeter is widely used in day due to the advantages that its measurement accuracy is high, range ratio is big, without crushing
The gas flow measurement field of the large diameter pipelines such as right gas.In measuring method, because propagation time difference method is caused by temperature change
Measurement error it is smaller, be the most frequently used, maximally effective measuring method.When measuring flow using propagation time difference method, it is necessary first to
Drive transmitting transducer to produce ultrasonic signal, obtained ultrasonic signal is then converted into echo letter by receive transducer again
Number；Corresponding following current, adverse current propagation time are calculated by measuring the time difference between transmitting and reception signal, obtains pipeline
The face velocity of cross section, and then calculate gas flow.The initial time of ultrasonic signal is known, and during echo-signal arrival
Quarter is unknown, it is necessary to be determined according to the characteristic point of a certain stabilization in echo-signal.Because ultrasonic signal is in gas
During propagation, its energy attenuation is serious, and can produce corresponding path offset with the increase of gas flow rate, further weakens energy
Amount so that echo-signal amplitude is faint, easily by noise jamming.And the feature directly determined by straightforward procedures such as maximal peak points
Point, larger fluctuation can be all produced, influences measurement accuracy.
Domestic and foreign scholars are studied gas ultrasonic flowmeter signal processing method, mainly there is following several method.
(1) method based on fixed threshold
From the angle of echo-signal amplitude, the measuring method of fixed threshold be by setting a fixed threshold value,
When the amplitude of ultrasonic echo signal is more than this value, just using this moment as during the arrival of ultrasonic echo signal
Between, and try to achieve with this propagation time of ultrasonic echo signal.Tamotsu Kobayashi et al. (Kobayashi T.,
Shigeta K.,Fujii T.Transit-time difference type ultrasonic flowmeter[P].U.S.:
20020062690A1,2002-05-30.) describe fixed threshold method hardware realize, that is, utilize comparator and with reference to electricity
Pressure value determines the arrival time of ultrasonic echo signal.R.Bates et al. (Bates R., et al.Development of
a custom on-line ultrasonic vapour analyzer/flowmeter for the ATLAS inner
detector,with application to gaseous tracking and Cherenkov detectors[J]
.Journal of Instrumentation,2013,8:C01002. the software realization of above-mentioned fixed threshold method) is refer to,
By starting clock while ultrasonic signal is sent, and stop clock after echo-signal amplitude exceedes preset value to count
Calculate the propagation time of ultrasonic signal.However, service life of noise pollution, drift or ultrasonic transducer etc. can all cause
The change of ultrasonic echo signal amplitude, once ultrasonic echo signal amplitude fluctuations are larger, obtained according to fixed threshold method
The echo signal reception time will malfunction.
(2) method based on energy gradient
A kind of letter based on energy gradient is it is also proposed from the angle, Shen Ziwen et al. of echo-signal energy variation
Number processing method (Shen Ziwen, Xu Kejun, Fang Min, waits Ultrasonic Wave Flowmeter signal processing methods of the based on energy gradient
[J] Chinese journal of scientific instrument, 2015,36 (9):2138-2144.), the pretreatment such as it is filtered to the echo-signal of sampling first,
Then each peak point of echo-signal is found, and to the progress square of its amplitude, the energy (energy of signal as the peak point
Squared magnitude with signal is in proportionate relationship, FANGM, XU K J, ZHU W J, et al.Energy transfer model
and its applications of ultrasonic gas flow-meter under static and dynamic
flow rates.[J].Review of Scientific Instruments,2016,87(1):015107.), obtain each
The energy gradient of peak point, and envelope fitting is carried out to it, in order to reduce amount of calculation, ensure the real-time of signal transacting,
The envelope of ultrasound echo signal energy is asked for using the method for linear interpolation.Energy curves in the case where obtaining each flow
Afterwards, the threshold value for distinguishing characteristic ripple is set according to energy curves.When actually measuring, according to threshold value set in advance and
The energy curves of echo-signal, characteristic point corresponding to differentiation are asked in real time.Because this method is needed before flow is measured
Threshold parameter selection is carried out, operation is relative complex.
(3) method based on energy jump
Daniel companies of the U.S. using detection energy jump method search characteristic point method (William Freund,
Winsor Letton,James Mc-Clellan,Baocang Jia,Anni Wey,Wen Chang.Method and
apparatus for measuring the time of flight of a signal,US patent NO.5983730,
Nov.16,1999).Because the energy of echo-signal is during being entirely activated to propagation and arriving reception again, from weak to strong, then from
Arrive by force weak, therefore backward energy rate of change will first increase and reduce again, the threshold variations point according to energy gradient, you can it is determined that super
Acoustic transit time.Specifically, the squared magnitude of echo-signal each point is asked for first；Then, asked for often using moving average
The average energy value of a bit, then draw out the change curve of the ratio between forward and backward 2 average energies, i.e. energy gradient curve；Most
Afterwards, the critical point of energy gradient is found, determines ultrasonic propagation time.But the patent does not disclose the crucial skill of method
Art, as the time interval of sliding window is chosen, and ask for time interval during energy gradient between forward and backward 2 points etc..
The content of the invention
It is difficult to determine stable characteristic point for Ultrasonic Wave Flowmeter echo-signal, when judging the arrival of echo-signal
Carve, and the problem of instrument is required compared with high real-time, the present invention proposes the gas ultrasound ripple integrated based on backward energy
Flow meter signal processing method.Because the energy of signal and the squared magnitude of signal are in proportionate relationship, therefore, in order to simplify meter
Calculate, the present invention uses the energy of square sign echo-signal of echo-signal amplitude.In the hardware using DSP and FPGA as double-core
(Xu Kejun, Fang Min, Wang Wei, red legend is handsome, a kind of Ultrasonic Wave Flowmeters based on FPGA and DSP of Shen Ziwen, Shen on platform
Please Chinese invention patent, publication number CN104697593A, applying date 2015.06.10), write software program, real-time implementation base
In the signal processing method of backward energy integration.
The calculation procedure of Ultrasonic Wave Flowmeter signal processing method based on backward energy integration is：Basis is adopted first
The ultrasonic echo signal of sample, pass through the maximum positive peak point for comparing to obtain backward energy signal；Then determine that section is risen forward
Initial point, with 0.1 times more than echo-signal normalization amplitude for Rule of judgment, find integration starting point；From integration starting point (product
Point starting point is also a peak point) arrive between maximum positive peak point, a peak point is often found, will be returned after the squared magnitude
One changes, and calculates the trapezoidal area surrounded between adjacent peak point and adds up；Judge the size of cumulative area value and preset area value；When
Stopping integration is cumulative after cumulative area value exceedes preset area value, subtracts preset area by cumulative area value and is worth to difference in areas
It is worth (small trapezoidal area)；The ordinate of known area difference (small trapezoidal area) and adjacent peak point (small trapezoidal upper bottom
Value), the intersection point (trapezoidal height) of the small trapezoidal upper bottom edge of solution and X-axis, as characteristic point；Determine the suitable of ultrasonic echo signal
Stream, adverse current propagation time, and calculate gas flow.
Brief description of the drawings
Fig. 1 is the signal processing method flow chart based on backward energy integration.
Fig. 2 is 0m3The backward energy signal and its bag under backward energy signal and its envelope and different flow under/h
Winding thread figure.
Fig. 3 is to determine integrating range starting point schematic diagram.
Fig. 4 is to determine selection integration starting point schematic diagram.
Fig. 5 is the integral area schematic diagram based on backward energy integration method.
Fig. 6 is to determine preset area value schematic diagram.
Fig. 7 is the feature point calculating method based on backward energy integration.
Fig. 8 is the hardware block diagram of Ultrasonic Wave Flowmeter.
Fig. 9 is gas ultrasonic flowmeter system master chip dsp software flow chart.
Embodiment
Fig. 1 is the signal processing method flow chart based on backward energy integration.After echo-signal is filtered etc. into processing,
Just start to calculate characteristic point.
(1) echo-signal maximum positive peak point and minimal negative peak point are found.In the echo signal data of collection, pass through
Point-by-point comparison is met two peak points of condition.
(2) will be compared after two peak points square.It is in proportionate relationship according to the squared magnitude of the energy of signal and signal,
By maximum positive peak point and minimal negative peak point respectively square after compare size.To intactly represent backward energy signal, then
Need all to carry out the amplitude of the echo-signal every bit of collection square, so processing can increase the operand of algorithm, influence journey
The whole efficiency of sequence.The present invention is carried out using the ascent stage each peak point for finding echo-signal to the amplitude of each peak point
Square processing, can both make full use of the information entrained by backward energy signal, the amount of calculation of whole algorithm can be reduced again.And
And after a peak point is searched out, just carry out follow-up area integral and calculate, (cumulative area value is big when reaching preparatory condition
In preset area value SA) after can complete to calculate, further reduce calculating process repetitive cycling.Reduced using such flow
The operands of data, it is favorably improved the real-time of system.
(3) backward energy signal maximum positive peak point is obtained.To determine that integration starting point is prepared, backward energy signal is most
The amplitude (ordinate) of big positive peak point uses AmaxRepresent.
(4) integration starting point, squared magnitude normalization are determined.In the case of known maximum positive peak point, find first
The section starting point of integration, then according to the requirement (0.1 times that is more than echo-signal normalization amplitude) of setting, find and integrate
Initial point, area is integrated since being integrated starting point.The requirement specifically set is described in detail in Fig. 3.
(5) quadrature area.Start reference area value after finding integration starting point, specifically include 6 steps：
Quadrature the first step of area：First peak point after searching integration starting point.
Quadrature the second step of area：To being normalized after the peak point squared magnitude.Peak point amplitude AnRepresent, return
Amplitude after one change is all usedRepresent, the backward energy signal amplitude of the peak point is its square, uses An 2Represent, n is represented from product
Divide starting point to each peak point subscript between maximum positive peak point.Normalizing formula is
Quadrature the 3rd step of area：Trapezoidal area between two neighboring peak point is calculated, and is added up.It is as shown in figure 5, adjacent
Two peak points can be surrounded with axis of abscissas one it is small trapezoidal, obtain the energy of echo-signal by the way that trapezoidal area is cumulative
Value.Two neighboring peak point (amplitude normalization) coordinate is used respectively(Xn+1,) represent, n and n+1 represent from
Starting point is integrated to each peak point subscript between maximum positive peak point.Small trapezoidal area SnCalculation formula is
All small trapezoidal areas are added up, obtain cumulative area value.Cumulative area S calculation formula is
Quadrature the 4th step of area：Whether preset area value S is more thanAJudge cumulative area value and preset area value
Size：
If cumulative area value is not up to preset area value, normalized amplitude is (under trapezoidal after preserving the peak point square
Bottom), find next peak point backward, then repeat second and third, four steps, continue to determine whether to be more than preset area value SA；
If cumulative area value has reached preset area value, next step is proceeded by.
Quadrature the 5th step of area：Reference area difference S'.Cumulative area value subtracts preset area and is worth to difference in areas
Value S'.
Quadrature the 6th step of area：Special sign point Transverse coordinates are asked by area difference.The step describes in detail in the figure 7.
(6) value of characteristic point is obtained.Following current, adverse current propagation time can then be obtained by the value of characteristic point, obtain pipeline
The face velocity of cross section, and then calculate gas flow.With a certain sound channel, downstream propagation times tdExemplified by, calculation formula is：
td=xtIn/f (4) formula, xtFor the characteristic point abscissa value being calculated, f is the sample frequency of echo-signal；
Gas flow Q calculation formula is：
In formula,For the average speed in sound channel, S is the cross-sectional area of pipeline, and α represents one and number of channels, sound channel point
The relevant correction factor of the factor such as cloth and Reynolds number, L represent sound channel length, and θ represents angle between transducer and gas pipeline,
td、tuUltrasonic wave following current, the adverse current propagation time of same sound channel are represented respectively.
Fig. 2 is 0m3The backward energy signal and its bag under backward energy signal and its envelope and different flow under/h
Winding thread figure.The figure of top half is 0m3Backward energy signal and its envelope under/h, this is by by 50 groups of 0m3Returning under/h
The pretreatments such as ripple signal is filtered, squared magnitude, averaging and normalization, using curve-fitting tool in MATLAB to echo
All peak points of energy signal carry out envelope and are fitted what is obtained.The figure of the latter half is the backward energy letter under different flow
Number and its envelope, by 0m3/h、200m3/h、400m3/h、600m3/h、800m3/h、1000m3Each 50 groups of echo-signals under/h
Carry out what same processing fitting obtained.It is can be seen that by different flow backward energy signal envelope when flow increases,
The backward energy signal decline stage can be distorted, and (there are some researches show when detected fluid flow velocity increases to certain limit, return
Ripple signal can produce serious distortion.SABATINI A M.A digital-signal-processing technique for
ultrasonic signal modeling and classification[J].IEEE Transactions on
Instrumentation Measurement,2001,50(1):15-21.), the amplitude of echo-signal, which can undergo, starts from scratch soon
Speed increase, the process of slow-decay to zero after maximum is reached, i.e., during releasing energy, be slow and jiggly release
Put.But in the ascent stage of backward energy signal, because the ultrasonic wave of excited target transducer transmitting directly acts on, substantially not
Exist reflection sound wave etc. influence, energy concentrate and stably, thus, it can be seen that each represented by 1,2,3,4,5,6 in fig. 2
The waveform profiles of ascent stage under individual flow are basically unchanged, and as the increase of flow, the waveform profiles of ascent stage occur
Translation, but waveform profiles are substantially parallel under each flow.Because the energy entrained by the ultrasonic wave launched every time is identical
(frequency, amplitude, the phase of pumping signal are all identical), so, the energy entrained by backward energy signal ascent stage is same
It is stable under flow, is integrated by the area surrounded to the ascent stage, the characteristic point that can be stablized, and then try to achieve
The flow of tested gas.
Fig. 3 is to determine that integrating range originates point diagram.When the flow velocity of tested gas exceedes certain scope, i.e., in big flow
Under, the noise signal adulterated in echo-signal substantially becomes big, it can be seen that the noise signal amplitude after being normalized in Fig. 3 is
Through more than 0.1, in fig. 2, the energy for having been described above the same flow lower ascent stage is stable, so the section end point of integration
◇ is assured that as maximal peak point, then, determine that the section starting point that area starts to integrate is just very crucial.If section
Starting point is chosen at first point of sampling, then noise signal will be taken as echo-signal to calculate into area, not only increase
Add amount of calculation, interference is introduced and calculated, and the real information of echo-signal can not be reflected, reduce measurement accuracy.And by most
Big peak point determines forward section starting point, although the position (Transverse coordinates of peak-peak) it may fluctuate, can have
Effect removes noise jamming, can reduce the influence of positional fluctuation by algorithm process.Take into account big flow and low discharge ascent stage
Peak value, choose maximal peak point 206 to 220 points (each 2048 points of echo signal sample, each 25 points of cycle) forward
Position, as integrating range starting point, as shown in the o in Fig. 2.Integrate starting point ● selection describe in detail in Fig. 4.
Fig. 4 is to determine selection integration starting point diagram.The present invention, by many experiments, is finding area with exhaustive method
Between after starting point, normalized from section starting point to section end point (maximal peak point), setting with echo-signal amplitude
0.1 on the basis of (after square represent backward energy signal, then amplitude be normalized to 0.01, the calculation process according to Fig. 1,
First find the peak point of echo-signal, then square be the peak point of backward energy signal), when normalized peak point amplitude expires
It is exactly required peak point, these peak points are as follow-up areal calculation best results when foot is more than or equal to 0.1.Such as Fig. 4 institutes
Show, in integration starting point ● the peak value between the end point ◇ of section is all the peak point for the condition that meets.Can from Fig. 4
Go out, when backward energy signal normalization amplitude is 0.01, the area very little that is covered before the starting point o of section with X-axis, i.e.,
The noise signal of signal leading portion disturbs very little, can be ignored.
Fig. 5 is the integral area schematic diagram based on backward energy integration method.It is determined that area integral section starting point,
And then after finding integration starting point, it can both carry out specific areal calculation.In figure ● representative be exactly we
Find and meet the integration starting point that imposes a condition and ascent stage each peak point, including maximal peak point b shown in figure, secondary big peak
It is worth the peak point c of point a and interlude.The method for using linear fit at 2 points between each peak point is adjacent, structure are into a line
Section, then the trapezoidal base (solid line in figure) adjacent with representated by this 2 points ordinate, axis of abscissas collectively constitute one
It is individual small trapezoidal, each small trapezoidal area is calculated and added up, while precision is ensured, simplifies calculating process, avoids
To amount of calculation caused by signal node-by-node algorithm.To for example above-mentioned formula (2) of each small trapezoidal areal calculation Suo Shu.
Fig. 6 is to determine preset area value schematic diagram.After integration starting point and ascent stage each peak point required for finding, I
Need set an area preset value, when cumulative area value exceedes preset area value, just stopping accumulation calculating area, open
The abscissa of the characteristic point to begin needed for solution.So, this area preset value is equally very crucial.When the area value that we choose
When smaller, corresponding peak value points are just few, i.e. for the covered area of institute at the peak point c of interlude, it is exactly our solutions to put 3
Obtained characteristic point；When the area value that we choose is larger, corresponding peak value points are just more, i.e. the covered area of institute is secondary
At big peak point a, point 1 is exactly that we solve obtained characteristic point.When backward energy signal waveform fluctuates, backward energy
Signal fluctuation causes peak point amplitude to change, and then the cumulative area value being calculated can change, and fluctuation causes to add up
The part that area value increases or decreases with the trapezoidal area of last part in cumulative area value, it is necessary to be compensated, i.e., ladder in figure
Shape S1 or trapezoidal S3 area value (two parts shaded area in figure) compensates, because area preset value is certain,
The area value for needing to compensate part is equal, i.e., trapezoidal S1 area value is equal to trapezoidal S3 area values, and because of trapezoidal S3's
Upper bottom is significantly less than trapezoidal S1 upper bottom, so, trapezoidal S3 height (the peak point c of interlude abscissa and spy
Sign point the distance between 3 abscissas) it is higher (between a of second largest peak value abscissa and the abscissa of characteristic point 1 than trapezoidal S1
Distance) it is big more, i.e., required by the obtained fluctuation of the Transverse coordinates of characteristic point 1 of solution to be significantly less than the fluctuation of the abscissa of characteristic point 3, that
, the stability of the abscissa of characteristic point 1 is higher than the stability of the abscissa of characteristic point 3.The stability of characteristic point is to ensure gas
The key of flow measurement precision.So area preset value should be chosen at second largest peak value point a and maximal peak point b (shown in Fig. 5)
Between, and no more than maximal peak point b, area discrepancy caused by the peak value of big flow and low discharge ascent stage is taken into account,
The area value after amplitude calculates after normalizing is chosen, determines preset area value SADuring between 60 to 63, best results.
Fig. 7 is the feature point calculating method based on backward energy integration.Area preset value be chosen at second largest peak value point a and
Maximal peak point b, when cumulative area value is more than preset area value SAWhen, stop cumulative, it is known that difference in areas (trapezoidal area) and phase
The ordinate (upper bottom value) of adjacent peak point, it is possible to solve trapezoidal height.Second largest peak value point a and maximal peak point b are (such as
Shown in Fig. 5) coordinate be respectively (x1, y1), (x2, y2), wherein, y1 and y2 are the width of normalized backward energy signal
Value, then, the slope k of linear fit is between this 2 points
Cumulative area value and preset area value SABetween area difference S' be
The distance between the abscissa of characteristic point 1 and maximal peak point b abscissas are x, i.e., trapezoidal S' height, then, S''s
Area value is
Unknown number x (an only solution is eligible) can be obtained according to formula (6)：
Then the value of final the solved abscissa of characteristic point 1 is (x2-x).
Fig. 8 is the hardware block diagram of Ultrasonic Wave Flowmeter.It is minimum mainly to include primary instrument, DSP minimum systems, FPGA
System, drive signal generation and amplifying circuit, transmitting/receiving signal channel switching circuit, echo-signal conditioning and Acquisition Circuit,
Power management module.
The primary instrument is made up of four transducers and flowtube, and four transducers are respectively and fixedly installed to flowtube
On；Four transducers are respectively transducer 1, transducer 2, transducer 3 and transducer 4, set up two-channel direct-injection type structure；Transducing
Device is transceiver type, is controlled by transmitting/receiving signal channel switching circuit, and it is sharp for launching to can not only be used for transmitting transducer
Signal is encouraged, is used for receives echo-signal but also as receive transducer.
The FPGA minimum systems are made up of fpga chip, clock, reset circuit and FLASH (flash memory), wherein, FPGA cores
Piece model EP2C8Q208C8N.Include DAC (digital analog converter) control module, ROM (read-only storage), frequency dividing inside FPGA
Module, RAM_2SPORT.
The DSP minimum systems are by dsp chip, clock, outer watchdog, FRAM (ferroelectric memory), man-machine interface, string
Port communications form, wherein, dsp chip model TMS320F28335.
The power management module for system provide needed for amplitude digital power and analog power.
Fig. 9 is gas ultrasonic flowmeter system master chip dsp software flow chart.After system electrification, DSP carries out initial
Change, complete the distribution to GPIO (universal input/output) mouth, timer interruption initialization, liquid crystal display initialization, flowmeter show
Show the initialization tasks such as initialization.DSP sends excitation initial signal to FPGA, and switches transmitting, receiving channel.Detect timer
Whether time tag position is 0, if 0, is then directly entered DSP and receives unloading initial signal wait state；If 1, then liquid is updated
Crystalline substance display, communicates with host computer, and by timer time mark position 0.After if DSP receives unloading initial signal, start to read
The echo-signal deposited in FPGA is taken, integral area is calculated, tries to achieve characteristic point abscissa, obtain the propagation time, seek each sound channel
Flow, and then export final flow rate Q.
Claims (4)
1. a kind of Ultrasonic Wave Flowmeter signal processing method and system based on backward energy integration, including primary instrument,
DSP minimum systems, FPGA minimum systems, drive signal generation and amplifying circuit, transmitting/receiving signal channel switching circuit, return
Ripple signal condition and Acquisition Circuit, the power management module and Ultrasonic Wave Flowmeter signal transacting based on backward energy integration
Software；After system electrification, DSP is initialized, and completes the distribution to GPIO mouths, timer interruption initializes, at the beginning of liquid crystal display
The initialization tasks such as beginningization, flowmeter display initialization；DSP sends excitation initial signal to FPGA, and switches transmitting, receives and lead to
Road；After DSP receives unloading initial signal, then start to read the echo-signal deposited in FPGA；Integral area is calculated, tries to achieve feature
Point abscissa；The propagation time is obtained, seeks the flow of each sound channel, and then export final flow rate Q；It is characterized in that：Surpass for gas
Acoustic wave flow meter echo-signal is difficult to determine stable characteristic point, judges the due in of echo-signal, and instrument requires have
There is the problem of compared with high real-time, it is proposed that the Ultrasonic Wave Flowmeter signal processing method based on backward energy integration, because
The energy of signal and the squared magnitude of signal are in proportionate relationship, therefore, are calculated to simplify, use square of echo-signal amplitude
Characterize the energy of echo-signal.
2. a kind of Ultrasonic Wave Flowmeter signal processing method based on backward energy integration as claimed in claim 1 and it is
System, it is characterised in that：After echo-signal is filtered etc. into processing, just start to calculate characteristic point；
(1) echo-signal maximum positive peak point and minimal negative peak point are found；
(2) will be compared after two peak points square；
(3) backward energy signal maximum positive peak point is obtained；
(4) integration starting point, squared magnitude normalization are determined；
(5) quadrature area.Start reference area value after finding integration starting point, specifically include 6 steps：
Quadrature the first step of area：First peak point after searching integration starting point；
Quadrature the second step of area：To being normalized after the peak point squared magnitude；
Quadrature the 3rd step of area：Trapezoidal area between two neighboring peak point is calculated, and is added up；
Quadrature the 4th step of area：Whether preset area value S is more thanAJudge the size of cumulative area value and preset area value；
Quadrature the 5th step of area：Reference area difference S'；
Quadrature the 6th step of area：Special sign point Transverse coordinates are asked by area difference；
(6) value of characteristic point is obtained；Following current, adverse current propagation time can then be obtained by the value of characteristic point, obtain conduit cross-sectional
The face velocity in face, and then calculate gas flow.
3. a kind of Ultrasonic Wave Flowmeter signal processing method based on backward energy integration as claimed in claim 1 and it is
System, it is characterised in that：In the case of known maximum positive peak point, choosing maximal peak point, 206 to 220 points are (each to return forward
2048 points of ripple signal sampling, each 25 points of cycle) position, as integrating range starting point；Finding section starting point
Afterwards, from section starting point to the end point of section, setting on the basis of echo-signal amplitude normalized 0.1, normalization is worked as
Peak point amplitude meet be more than or equal to 0.1 when, be exactly required peak point.
4. a kind of Ultrasonic Wave Flowmeter signal processing method based on backward energy integration as claimed in claim 1 and it is
System, it is characterised in that：Area preset value should be chosen between second largest peak value point a and maximal peak point b, and no more than maximum peak
It is worth point b, chooses the area value after amplitude calculates after normalizing, determine preset area value SABetween 60 to 63.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711043275.3A CN107655533B (en) | 2017-10-31 | 2017-10-31 | A kind of Ultrasonic Wave Flowmeter signal processing method and system based on backward energy integral |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711043275.3A CN107655533B (en) | 2017-10-31 | 2017-10-31 | A kind of Ultrasonic Wave Flowmeter signal processing method and system based on backward energy integral |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107655533A true CN107655533A (en) | 2018-02-02 |
| CN107655533B CN107655533B (en) | 2019-09-10 |
Family
ID=61095941
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711043275.3A Active CN107655533B (en) | 2017-10-31 | 2017-10-31 | A kind of Ultrasonic Wave Flowmeter signal processing method and system based on backward energy integral |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107655533B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110554378A (en) * | 2019-09-10 | 2019-12-10 | 深圳锐越微技术有限公司 | Single-channel Doppler radar radial motion direction identification method and device |
| CN112890859A (en) * | 2021-01-19 | 2021-06-04 | 广州多浦乐电子科技股份有限公司 | Data processing method and system based on ultrasonic signal feature vector |
| CN113272518A (en) * | 2018-11-29 | 2021-08-17 | Bp探索操作有限公司 | DAS data processing to identify fluid inflow location and fluid type |
| CN115078909A (en) * | 2022-06-16 | 2022-09-20 | 国网江苏省电力有限公司宿迁供电分公司 | Non-ideal pulse leading edge detection method and device based on energy inflection point |
| CN115655396A (en) * | 2022-12-26 | 2023-01-31 | 水发航宇星物联科技(辽宁)有限公司 | Ultrasonic wave fault wave prevention method and ultrasonic wave fluid detection device |
| CN116211453A (en) * | 2023-05-09 | 2023-06-06 | 浙江伽奈维医疗科技有限公司 | Temperature control method and system for varicose vein radio frequency ablation catheter |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4670238B2 (en) * | 2003-12-10 | 2011-04-13 | パナソニック株式会社 | Gas shut-off device |
| CN102113899A (en) * | 2009-12-30 | 2011-07-06 | 东软飞利浦医疗设备系统有限责任公司 | Method for extracting two-dimensional dynamic blood vessel information based on ultrasonic echo in real time |
| CN102323592A (en) * | 2011-08-12 | 2012-01-18 | 中国科学院光电技术研究所 | Normalization method for target echo signal |
| CN106643940A (en) * | 2017-02-20 | 2017-05-10 | 重庆川仪自动化股份有限公司 | Method for computing propagation time of ultrasonic flow meters on basis of echo energy |
| CN106768109A (en) * | 2017-02-21 | 2017-05-31 | 合肥工业大学 | Ultrasonic Wave Flowmeter signal processing method based on echo ascent stage peak fitting and based on backward energy point location |
| CN106871980A (en) * | 2017-02-13 | 2017-06-20 | 合肥工业大学 | The excitation of gas ultrasonic flowmeter batch (-type) and signal processing method and system based on adjacent peak maximum difference |
-
2017
- 2017-10-31 CN CN201711043275.3A patent/CN107655533B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4670238B2 (en) * | 2003-12-10 | 2011-04-13 | パナソニック株式会社 | Gas shut-off device |
| CN102113899A (en) * | 2009-12-30 | 2011-07-06 | 东软飞利浦医疗设备系统有限责任公司 | Method for extracting two-dimensional dynamic blood vessel information based on ultrasonic echo in real time |
| CN102323592A (en) * | 2011-08-12 | 2012-01-18 | 中国科学院光电技术研究所 | Normalization method for target echo signal |
| CN106871980A (en) * | 2017-02-13 | 2017-06-20 | 合肥工业大学 | The excitation of gas ultrasonic flowmeter batch (-type) and signal processing method and system based on adjacent peak maximum difference |
| CN106643940A (en) * | 2017-02-20 | 2017-05-10 | 重庆川仪自动化股份有限公司 | Method for computing propagation time of ultrasonic flow meters on basis of echo energy |
| CN106768109A (en) * | 2017-02-21 | 2017-05-31 | 合肥工业大学 | Ultrasonic Wave Flowmeter signal processing method based on echo ascent stage peak fitting and based on backward energy point location |
Non-Patent Citations (4)
| Title |
|---|
| JIE YANG: "Design of Multi-channel Ultrasonic Flowmeter Based on ARM", 《INTERNATIONAL CONFERENCE ON ELECTRICAL & CONTROL ENGINEERING》 * |
| 汪伟: "一种气体超声波流量计信号处理方法研究", 《电子测量与仪器学报》 * |
| 沈子文: "基于能量变化率的气体超声波流量计信号处理方法", 《仪器仪表学报》 * |
| 田雷: "基于回波峰值拟合的气体超声流量计信号处理", 《电子测量与仪器学报》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113272518A (en) * | 2018-11-29 | 2021-08-17 | Bp探索操作有限公司 | DAS data processing to identify fluid inflow location and fluid type |
| CN110554378A (en) * | 2019-09-10 | 2019-12-10 | 深圳锐越微技术有限公司 | Single-channel Doppler radar radial motion direction identification method and device |
| CN110554378B (en) * | 2019-09-10 | 2021-05-11 | 深圳锐越微技术有限公司 | Single-channel Doppler radar radial motion direction identification method and device |
| CN112890859A (en) * | 2021-01-19 | 2021-06-04 | 广州多浦乐电子科技股份有限公司 | Data processing method and system based on ultrasonic signal feature vector |
| CN112890859B (en) * | 2021-01-19 | 2024-01-09 | 广州多浦乐电子科技股份有限公司 | Data processing method and system based on ultrasonic signal feature vector |
| CN115078909A (en) * | 2022-06-16 | 2022-09-20 | 国网江苏省电力有限公司宿迁供电分公司 | Non-ideal pulse leading edge detection method and device based on energy inflection point |
| CN115078909B (en) * | 2022-06-16 | 2024-03-15 | 国网江苏省电力有限公司宿迁供电分公司 | Non-ideal pulse leading edge detection method and device based on energy inflection point |
| CN115655396A (en) * | 2022-12-26 | 2023-01-31 | 水发航宇星物联科技(辽宁)有限公司 | Ultrasonic wave fault wave prevention method and ultrasonic wave fluid detection device |
| CN116211453A (en) * | 2023-05-09 | 2023-06-06 | 浙江伽奈维医疗科技有限公司 | Temperature control method and system for varicose vein radio frequency ablation catheter |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107655533B (en) | 2019-09-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107655533B (en) | A kind of Ultrasonic Wave Flowmeter signal processing method and system based on backward energy integral | |
| CN106768109B (en) | Based on echo ascent stage peak fitting and based on the Ultrasonic Wave Flowmeter signal processing method of backward energy point location | |
| CN103344288B (en) | A kind of transit-time ultrasonic flow meter measuring method analyzed based on zero point | |
| CN106643939B (en) | The method for calculating ultrasonic propagation time for ultrasonic flowmeter | |
| CN103776516B (en) | Ultrasonic flow metering system with upstream pressure transducer | |
| JP2003344131A (en) | Ultrasonic flow velocity distribution and flowmeter | |
| CN107860430B (en) | Time difference measurement method of ultrasonic gas flowmeter based on time difference method | |
| CN110186522B (en) | Moisture overreading compensation and flow measurement method combining vortex street amplitude characteristic | |
| CN107843305B (en) | A kind of Ultrasonic Wave Flowmeter signal processing method based on echo signal envelope fitting | |
| CN107328447A (en) | Well head gas-liquid biphase flowmeter amount method and device | |
| CN110646042A (en) | Cross-correlation interpolation method for calculating flight time difference of low-power-consumption ultrasonic flowmeter | |
| Tian et al. | Energy peak fitting of echo based signal processing method for ultrasonic gas flow meter | |
| CN105181793A (en) | Method for measurement of two-phase flow gas holdup based on ultrasonic two-frequency signal | |
| CN111044110A (en) | Gas ultrasonic flowmeter signal processing method based on similarity analysis | |
| CN112304376B (en) | Ultrasonic flowmeter flow measuring method based on data fusion | |
| US10416013B2 (en) | Ultrasonic flow measuring apparatus detecting propagation time period from a zero cross point when a peak waveform exceeds a trigger level | |
| CN207036192U (en) | Gas liquid two purpose ultrasonic flowmeter | |
| US9962067B2 (en) | Real time detection of aspiration short shots using pressure signal | |
| CN113375737A (en) | Flow velocity metering method of time difference type ultrasonic gas flowmeter | |
| KR101764870B1 (en) | Signal processing system for ultrasonic floemeter | |
| CN116558587A (en) | Ultrasonic fluid flow measurement method and system based on flow prediction | |
| Ma et al. | Signal processing method based on connection fitting of echo peak point with a large slope for ultrasonic gas flow meter | |
| CN110131591B (en) | Method, device and equipment for positioning pipeline leakage | |
| Zheng et al. | Study on the measurement method of wet gas flow velocity by ultrasonic flow meter | |
| CN104713631B (en) | The detection method of bulk sound velocity in a kind of oil well |
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 |