CN108981857A - Acoustic resonance level measuring method based on iteration Weighted Fusion - Google Patents
Acoustic resonance level measuring method based on iteration Weighted Fusion Download PDFInfo
- Publication number
- CN108981857A CN108981857A CN201810901827.8A CN201810901827A CN108981857A CN 108981857 A CN108981857 A CN 108981857A CN 201810901827 A CN201810901827 A CN 201810901827A CN 108981857 A CN108981857 A CN 108981857A
- Authority
- CN
- China
- Prior art keywords
- frequency
- sequence
- group
- point
- leak source
- 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
- G01F23/00—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm
- G01F23/22—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water
- G01F23/28—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water by measuring the variations of parameters of electromagnetic or acoustic waves applied directly to the liquid or fluent solid material
- G01F23/296—Acoustic waves
- G01F23/2966—Acoustic waves making use of acoustical resonance or standing waves
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Electromagnetism (AREA)
- Thermal Sciences (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a kind of acoustic resonance level measuring methods based on iteration Weighted Fusion.It is able to solve the decline problem of measurement accuracy caused by existing measurement method resonant frequency point missing inspection.The present invention under same liquid level, successively acquires multi-group data using hardware system, and respectively this several groups of data are carried out with the extraction of Fast Fourier Transform and resonant frequency point.The average value of difference between first group of adjacent resonant frequency point is found out, leak source confidence interval is set, finds out leak source position and leak source numerical value, obtains having mended the new frequency sequence formed after leak source.It is successively iterated with the average value of difference between new frequency sequence consecutive points with next group of mean value and merges to obtain fused resonant frequency point sequence.Using between fused resonant frequency point sequence consecutive points difference and liquid level reduction formula obtain multiple level measurements, this multiple level measurements segmentation plus compensation are averaged as final liquid level height value.Accurate level gauging result can be obtained in the present invention.
Description
Technical field
The invention belongs to measuring techniques, are related to a kind of acoustic resonance level measuring method based on iteration Weighted Fusion.
Background technique
Level measuring method is widely used in many industrial circles such as petroleum, chemical industry, sewage treatment.But presently, there are
Level measuring method has had more than 20, for example utilizes hand-dipping, glass tube, float-type, capacitance method, ultrasonic wave, radar, light
Fibre etc. measures liquid level.Wherein, since Ultrasonic Liquid Level Measurement can realize contactless, high reliability measurement, therefore its
That applies is the most extensive.But in actual application, since the surface of fluid to be measured usually can be there is some foams, residual
The barriers such as slag.When ultrasonic wave encounters the foreign matter that these swim on liquid level, it is easy to happen parasitic reflection phenomenon, that is, is changed former
Come the direction of propagation of sound wave, to influence the degradation of liquid level measurement precision.
And because its wavelength is longer diffraction can occur when encountering foreign matter for low-frequency sound wave, i.e. sound wave can continue around foreign matter
It propagates forward, so as to avoid the generation of parasitic reflection.In document " the level gauging side based on sound wave on fixed frequency range resonance principle
Method " in, propose the level measuring method using sound wave on fixed frequency range resonance principle.This method is by the collected sound of microphone
Wave time-domain signal is transformed into frequency-region signal by FFT (Fast Fourier Transform (FFT)), and it is solid to obtain it by fast frequency detection method
Determine the resonant frequency point in frequency range, and is converted using the equal difference relationship between adjacent resonant frequency point and the liquid level based on resonant frequency
Formula obtains multiple groups liquid level height value, is then averaged to obtain final liquid level height value.But in actual application
Cheng Zhong, due to microphone sensitivity the problem of and ADC (analog-digital converter) module in signal conversion process can exist one
The presence of some uncertain factors such as fixed deviation and ambient temperature, noise, therefore collected acoustic signals exist
There can be the decaying of noise or signal to a certain extent.To which when extracting resonance point, partial resonance point possibly can not be got,
So that the accuracy of the level measuring method of resonance frequency of sound wave on fixed frequency range will be greatly reduced.
Summary of the invention
The purpose of the present invention is in view of the deficienciess of the prior art, propose a kind of acoustic resonance based on iteration Weighted Fusion
Level measuring method.
The present invention successively extracts being total to for the synthetic waveform that multiple groups are generated under same liquid level with sound wave in fixed frequency range
Then vibration frequency point is merged multiple groups resonant frequency by the method for iteration Weighted Fusion, and fusion results are passed through liquid
Position reduction formula is converted to liquid level information, specifically includes the following steps:
Phonocatheter is inserted into the bottom of testing liquid, microphone, loudspeaking by step (1) perpendicular to the liquid level of testing liquid
Device and temperature sensor are placed side by side on the other end of phonocatheter, and guarantee that loudspeaker and microphone are right against the liquid of testing liquid
Face;
It is one group of D=[A, B] Hz even variation at any time that step (2), which generates a band frequency range by controller,
Sine wave swept-frequency signal, frequency sweep time are Δ T, and Δ T ∈ [1s, 5s], the swept-frequency signal is by loudspeaker perpendicular to liquid level to be measured hair
It propagates out and in phonocatheter, the value of swept frequency range [A, B] Hz should meet A > > 20Hz and B-A > 3F0, wherein F0It is first
Resonant frequency under the conditions of beginning;When the sound wave that loudspeaker issues propagates to testing liquid liquid level through phonocatheter, sound wave can be sent out
Raw to reflect to form back wave, the frequency sweep sound wave that back wave and loudspeaker issue, which is overlapped and will form composite wave, propagates to Mike
At wind;Microphone collects synthetic waveform, then converts digital signal for acoustic signals by ADC (analog-digital converter);
Collected digital signal is made it by time domain wave by step (3) by the FFT (Fast Fourier Transform) of controller
Shape is converted to frequency-domain waveform;
Step (4) utilizes patent " level measuring method based on resonance frequency of sound wave on fixed frequency range " (CN
101852638) the fast frequency detection method in extracts the resonant frequency point in frequency domain;
In the case where same liquid level, repetition step (2)-(4) successively acquire N and are combined into wave train step (5),
When due to extracting resonance point by step (4), it may occur that the case where partial resonance point is omitted, that is, first group of collected frequency
Count n1Not equal to second group collected frequency points n2It is not equal to ... frequency points n collected not equal to N groupN, because
This demand takes the average value of the difference between its side frequency point, (is referred to as leaked using the resonance point that the mean value is sought being missed
Point), its step are as follows:
Step (5.1) seeks the average value of the difference between its adjacent two o'clock to resonant frequency point: being adopted by step (4)
Collect the resonant frequency point sequence that N is combined into wave, successively remembers that first group of resonance point sequence isSecond group
ForThird group is denoted as..., N group isf
Indicate frequency values, the first digit of lower target of f indicates that group number, lower target second indicate Frequency point institute position in the sequence
It sets, n1,n2,......,nNIt is the total number of resonant frequency point in each group sequence respectively, when being extracted due to step (4) resonant frequency
It will appear leak source, then have n1≠n2≠n3≠......≠nN, to P1Find out the average value of the difference between its adjacent two dot frequencyWherein
Leak source confidence interval is arranged in step (5.2), finds out leak source position and finds out the leak source: can obtain in step (5.1)
To first group of mean frequency valueFind out its frequency sequence withReliability distribution series when for benchmark, i.e.,Leak source confidence interval I is set1=[0.6,1.55), I2=[1.56,2.3), I3=[2.3,
3.0),I4=[3.0,4.0), I5=[4.0,4.9), I6=[4.9,5.59), I7=[5.59,6.5);Work as R1,i∈I1When, g (i)=
0;Work as R1,i∈I2When, g (i)=1;Work as R1,i∈I3When, g (i)=2;Work as R1,i∈I4When, g (i)=3;Work as R1,i∈I5When, g (i)=
4;Work as R1,i∈I6When, g (i)=5;Work as R1,i∈I7When, g (i)=6;Wherein g (i) is the number of i-th of frequency sequence section leak source,
G (i)=0,1,2 ... .., i=1,2 ... n1-1;When finding frequency sequence leak source section, when finding out leak source number, with the section
The equal part of difference between frequency takes the frequency values of leak source, T=g (1)+g (2)+... ..+g (i), i=1,2 ... .n1, t is all leakages in frequency sequence
The total number of point, the frequency sequence P after then obtaining first group of complement point1'(f1,1,f1,2..., f1,n1,f1,1',...,f1,t), and
To P1' resequenced to obtain new frequency sequence P1"={ f1(k), k=1,2 ... .., n1+ t }, k is after rearranging
Frequency sequence number, wherein f1(n1+t)>.....>f1(2)>f1(1);
Step (6) is iterated Weighted Fusion to second group of data that step (5) obtains: by obtaining in step (5.2)
Complement point after frequency sequence P1", the average value of the difference between adjacent two dot frequency is found out, new mean frequency value is obtained
Wherein k=1,2 ... .., n1+ t, n1For the total number of collected frequency sequence before first group of data complement point, t is
The total number of all leak sources;By obtaining second group of frequency sequence mean value in step (5.1)WithIt obtains new
Mean frequency value, and withFor the benchmark of second group of complement point, repeats step (5.2) and obtain the frequency sequence of second group of complete leak source of benefit
P2'(f2,1,f2,2..., f2,n2,f2,1',...,f2,t), P is obtained after rearrangement2"={ f2(k), k=1,2 ... .., n2+
T }, k is to rearrange rear frequency sequence number, wherein f2(n2+t)>.....>f2(2)>f2(1);By repeating step (5), (6)
Frequency sequence after available next group of complement point;
Step (7) is by obtaining the new frequency sequence P of N group after step (5), the fusion of (6) iterationN"={ fN(k), k=
1,2,…,nN+ t }, using the frequency equal difference relationship of adjacent resonance point, seek the intervening sequence e (j) of adjacent point-to-point transmission, e (j)=
fN(k+1)-fN(k) j=k, k=1,2 ... ..nN+ t-1, and it is right to calculate using following liquid level reduction formula sequence institute
The liquid level height value answered, available total nN+ t-1 liquid level L (j):
Wherein T is the temperature of environment, and t is the total number of leak source, and Δ L is compensation term.
In practical application, when liquid level is larger, sound wave experienced the phase offset (pressure of zero degree in every secondary reflection
Wave property), therefore there is certain " phase offset ", influence measurement result;That is, need to compensate error caused by this factor, it adopts
With segmented compensation function, liquid level is divided into three sections, each section of measurement height increases certain compensation rate Δ L
L (j) summation of acquirement is averaged, final liquid level L is obtained
Controller in this method, using the controller of the ARM framework with stronger processing power, such as the life of ST company
The development board of the STM32F7 system class of production;Audio processing chip used above only need to be performed acoustic signals generation conveying and
Acquisition, such as the audio decoder chips W M8994 carried on STM32F746G-DISCO development board.
The key technology of the above method is: successively acquiring multiple groups sound wave by reflecting superimposed composite wave signal, leads to
The resonant frequency point that quick frequency detecting method successively extracts each group composite wave is crossed, by the method for iteration Weighted Fusion by multiple groups
There are the resonance dot frequencies of the sound wave of leak source to be merged when probably due to extracting, and obtains one group of new resonance point frequency sequence,
To reduce because of the affected measurement results of uncertain factor some in the hardware systems such as microphone and ambient enviroment
Error, and then improve the accuracy of level gauging.
It can be effectively reduced using the method that the present invention mentions and extract the resonance frequency of composite wave in fast frequency detection method
Error caused by some uncertain factors when rate point in the presence of existing error and ambient enviroment, so that it is surveyed
Accuracy of measurement can reach 1%, better than the measurement accuracy of most of ultrasonic liquid levellers on the market, while also reduce measurement cost.
Detailed description of the invention
The hardware structure diagram of Fig. 1 the method for the present invention;
The flow chart of Fig. 2 the method for the present invention;
The time domain waveform of the collected sound wave of Fig. 3 microphone;
Fig. 4 is frequency domain figure corresponding to Fig. 3 time domain waveform;
The resonant frequency point for the sound wave that Fig. 5 is extracted;
Fig. 6 iteration Weighted Fusion method seeks leak source flow chart.
Specific embodiment
The present invention is further detailed with reference to the accompanying drawing.
The hardware configuration of acoustic resonance level measuring method based on iteration Weighted Fusion as shown in Figure 1, include controller 1,
Loudspeaker 2, microphone 3, temperature sensor 4, phonocatheter 5 and snorkel 6.Controller therein selects STM32F746G-DISCO
The onboard STM32F746G of development board, audio decoder chip is using carrying on STM32F746G-DISCO development board
WM8994, loudspeaker select model F10 dynamic speaker, and microphone selects TCM-340 electret mike, temperature sensor
DS18B20 is selected, phonocatheter selects the pvc pipe of diameter 7.5cm.Controller passes through loudspeaker for one group of fixed frequency range D=[A, B]
The sine wave swept-frequency signal of Hz exports in the form of a sound wave, and the acoustic signals are even variations.The acoustic signals are through leading sound
Pipe can reflect to form back wave after propagating to testing liquid surface, and the frequency sweep sound wave that back wave is issued with loudspeaker is folded
Add and to form composite wave and collected by microphone, digital signal is then converted for composite wave signal by ADC (analog-digital converter)
It is transmitted to controller.
Based on the acoustic resonance level measuring method process of iteration Weighted Fusion as shown in Fig. 2, core are as follows: in same liquid
Multiple groups composite wave signal is successively acquired under the height of position, obtains its resonant frequency point sequence 23 using fast frequency detection method.Meter
The average value 24 of the difference between first group of side frequency point is calculated, leak source confidence interval is set, finds out leak source position and finds out the leakage
Point 25 calculates the mean frequency value 26 of frequency sequence new after having mended leak source.With the average out to base of 24,26 resulting mean frequency value sums
Standard obtains the position and numerical value 27 of next group of leak source, acquires the difference 28 between the adjacent resonant frequency point of its new frequency sequence.It will
The interval of 28 obtained adjacent resonance points and the temperature value obtained according to temperature sensor, calculate the height 30 of liquid level,
Then segmented compensation 31 appropriate is carried out, final liquid level 32 is found out, carries out the output 32 of result.
Below according to the process of this method, and actual measurement environment and level gauging example are combined, detailed introduction should
Each step of method.
The setting of 1 test environment and parameter
Parameters in hardware configuration Fig. 1 are set: maximum measurement height w=10.6m, measurable maximum distance
hmax=10m, the length dz=0.6m in dead zone, the diameter d=0.07m of phonocatheter (pvc pipe), wherein the frequency sweep time is 2s, frequency sweep
Frequency range D=[1000,2000] Hz.Within the time of Δ t=2s, the swept frequency of sound wave is uniform from 1000Hz to 2000Hz
Variation.Acoustic signals are issued by loudspeaker, are reflected after reaching testing liquid surface by phonocatheter, then by microphone
It collects, this process can have deviation, therefore choosing frequency range is the spectrogram in [1000Hz, 2000Hz].Fig. 3, Fig. 4 difference
It is the time domain waveform and frequency-domain waveform of the collected acoustic signals of L=3.8m microphone.
2 are iterated fusion to wherein two groups of resonant frequency point sequences using iteration Weighted Fusion method
The flow chart of the step as shown in fig. 6, specific steps are as follows:
(1) mean value of resonant frequency point sequence is sought.As can be seen that being detected using fast frequency from the spectrogram of Fig. 5
When method extracts resonance point, it is possible that omitting some resonant frequency points, the accuracy rate of surveyed liquid level is caused significantly
It reduces.In order to reduce error caused by this reason, the resonant frequency point of omission is found out simultaneously using the method that iteration weights
It fills.It needs that data are normalized before fusion, such as two groups of resonant frequency sequences in Fig. 5WithThe number of the resonance point of this two groups of sound waves has n respectively1And n2It is a, and n1
≠n2, wherein f indicates frequency values.The average value of the difference between the adjacent resonant frequency point of this two groups of data is found out respectively
(2) leak source confidence interval is set, finds out leak source position and finds out leak source
By step (1) obtain the adjacent point-to-point transmission of two groups of resonant frequency point every mean valueWithWithFor the group
Benchmark find out the reliability distribution series of its frequency sequence,Leak source confidence interval is set
I1=[0.6,1.55) I2=[1.56,2.3), I3=[2.3,3.0), I4=[3.0,4.0), I5=[4.0,4.9), I6=
[4.9,5.59),I7=[5.59,6.5);Work as R1,i∈I1When, g (i)=0;Work as R1,i∈I2When, g (i)=1;Work as R1,i∈I3When,
G (i)=2;Work as R1,i∈I4When, g (i)=3;Work as R1,i∈I5When, g (i)=4;Work as R1,i∈I6When, g (i)=5;Work as R1,i∈I7
When, g (i)=6;Wherein number of the g (i) for i-th of frequency sequence section leak source, g (i)=0,1,2 ... .., i=1,2,
...n1-1;When finding frequency sequence leak source section, when finding out leak source number, taken with the equal part of the difference between the interval frequency
The frequency values of leak source, i.e.,
T=g (1)+g (2)+... ..+g (i), i=1,2 ... .n1, t is the total number of all leak sources in frequency sequence;Acquire first
Frequency sequence P after group complement point1'(f1,1,f1,2..., f1,n1,f1,1',...,f1,t), and carry out the P that resequences to obtain1"={ f1
(k), k=1,2 ... .., n1+ t }, k is to rearrange rear frequency sequence number, wherein f1(n1+t)>.....>f1(2)>f1(1);
Seek P1" side frequency point interval average value
(3) frequency sequence after seeking second group of complete leak source of benefit with the method for iteration
There are the influence of the uncertain factors such as error when influence and hardware system due to ambient enviroment measure, cause
To resonant frequency sequence itself there are certain uncertainties.In order to reduce error brought by this uncertainty, with the
The mean value at the new frequency sequence interval after one group of complement pointWith being averaged for the mean value sum before second group of complement pointTo reduce its survey
Error is measured, i.e.,WithTo seek second group of frequency sequence P2'(f2,1,f2,2..., f2,n2,f2,1',...,
f2,t), P is obtained after rearrangement2"={ f2(k), k=1,2 ... .., n2+t};Similarly, it can be with the method for iteration Weighted Fusion
Obtain the frequency sequence P of N groupN"={ fN(k), k=1,2 ..., nN+t}。
(4) difference between side frequency point is sought.Frequency point sequence after obtaining iteration Weighted Fusion according to above method
PN"={ fN(k), k=1,2 ..., nN+ t }, taking the intervening sequence between two adjacent Frequency points is e (j), wherein the tool of e (j)
Gymnastics is made as shown in formula (6).
E (j)=fN(k+1)-fN(k) j=k, k=1,2 ... ..nN+t-1 (6)
(5) height of liquid level is sought.Intervening sequence e (j) after being merged by above method between side frequency point,
Shown in the formula such as formula (7) for acquiring liquid level according to each sequence:
Wherein, T is the temperature of environment, and t is the total number of leak source, and Δ L is compensation term;In practical application, work as liquid level
When larger, sound wave experienced the phase offset (pressure wave characteristics) of zero degree in every secondary reflection, therefore " phase is inclined in the presence of certain
Move ", influence measurement result;That is, need to compensate error caused by this factor, using segmented compensation function, liquid level is divided into
Three sections, each section of measurement height increases certain compensation rate Δ L
L (j) summation of acquirement is averaged, final liquid level L is obtained
Claims (2)
1. the acoustic resonance level measuring method based on iteration Weighted Fusion, it is characterised in that method includes the following steps:
One end of phonocatheter is inserted into the bottom of testing liquid by step (1) perpendicular to the liquid level of testing liquid, and microphone is raised
Sound device and temperature sensor are placed side by side on the other end of phonocatheter, and guarantee that loudspeaker and microphone are right against testing liquid
Liquid level;
Step (2) generates the sine that a band frequency range is one group of D=[A, B] Hz even variation at any time by controller
Wave swept-frequency signal, the frequency sweep time be Δ T, Δ T ∈ [1s, 5s], the swept-frequency signal by loudspeaker perpendicular to liquid level to be measured issue simultaneously
It is propagated in phonocatheter, the value of swept frequency range [A, B] Hz should meet A > > 20Hz and B-A > 3F0, wherein F0For initial strip
Resonant frequency under part;When the sound wave that loudspeaker issues propagates to testing liquid liquid level through phonocatheter, sound wave can occur instead
It penetrates to form back wave, the frequency sweep sound wave that back wave and loudspeaker issue, which is overlapped and will form composite wave, propagates to microphone
Place;Microphone collects synthetic waveform, then converts digital signal for acoustic signals by analog-digital converter;
Collected digital signal is converted to it by time domain waveform by the Fast Fourier Transform of controller by step (3)
Frequency-domain waveform;
Step (4) is using in patent " level measuring method based on resonance frequency of sound wave on fixed frequency range " (CN 101852638)
Fast frequency detection method extract frequency domain in resonant frequency point;
In the case where same liquid level, repetition step (2)-(4) successively acquire N and are combined into wave train step (5), seek
The average value of difference between side frequency point, the resonance point being missed is sought using the average value, the specific steps are as follows:
Step (5.1) seeks the average value of the difference between its adjacent two o'clock to resonant frequency point: collecting N by step (4)
It is combined into the resonant frequency point sequence of wave, successively remembers that first group of resonance point sequence isSecond group isThird group is denoted as..., N group isF table
Show frequency values, the first digit of lower target of f indicates that group number, lower target second indicate Frequency point institute position in the sequence,
n1,n2,......,nNIt is the total number of resonant frequency point in each group sequence respectively, due to meeting when step (4) resonant frequency extracts
There is leak source, then has n1≠n2≠n3≠......≠nN, to P1Find out the average value of the difference between its adjacent two dot frequency
Wherein
Leak source confidence interval is arranged in step (5.2), finds out leak source position and finds out the leak source: the can be obtained in step (5.1)
One group of mean frequency valueFind out its frequency sequence withReliability distribution series when for benchmark, i.e.,Setting
Leak source confidence interval I1=[0.6,1.55), I2=[1.56,2.3), I3=[2.3,3.0), I4=[3.0,4.0), I5=[4.0,
4.9),I6=[4.9,5.59), I7=[5.59,6.5);Work as R1,i∈I1When, g (i)=0;Work as R1,i∈I2When, g (i)=1;Work as R1,i
∈I3When, g (i)=2;Work as R1,i∈I4When, g (i)=3;Work as R1,i∈I5When, g (i)=4;Work as R1,i∈I6When, g (i)=5;Work as R1,i
∈I7When, g (i)=6;Wherein number of the g (i) for i-th of frequency sequence section leak source, g (i)=0,1,2 ... ..;When find frequency
Rate sequence leak source section, when finding out leak source number, the frequency values of leak source is taken with the equal part of the difference between the interval frequency, under
The formula in face calculates:T=g (1)
+ g (2)+... ..+g (i), i=1,2 ... .n1, t is the total number of all leak sources in frequency sequence, then obtains first group of benefit
Frequency sequence P after point1'(f1,1,f1,2..., f1,n1,f1,1',...,f1,t), and to P1' resequenced to obtain new frequency
Rate sequence P1"={ f1(k), k=1,2 ... .., n1+ t }, k is to rearrange rear frequency sequence number, wherein f1(n1+t)>.....
>f1(2)>f1(1);
Step (6) is iterated Weighted Fusion to second group of data that step (5) obtains: by mending obtained in step (5.2)
Frequency sequence P after point1", the average value of the difference between adjacent two dot frequency is found out, new mean frequency value is obtained
Wherein n1For the total number of collected frequency sequence before first group of data complement point, t is the total number of all leak sources;Pass through
Second group of frequency sequence mean value is obtained in step (5.1)WithObtain new mean frequency value, and withIt is
The benchmark of two groups of complement points repeats step (5.2) and obtains the frequency sequence P of second group of complete leak source of benefit2'(f2,1,f2,2..., f2,n2,
f2,1',...,f2,t), P is obtained after rearrangement2"={ f2(k), k=1,2 ... .., n2+ t }, k is to rearrange rear frequency sequence
Columns, wherein f2(n2+t)>.....>f2(2)>f2(1);Frequency after obtaining next group of complement point by repetition step (5), (6)
Rate sequence;
Step (7) is by obtaining the new frequency sequence P of N group after step (5), the fusion of (6) iterationN"={ fN(k), k=1,
2,…,nN+ t }, using the frequency equal difference relationship of adjacent resonance point, seek the intervening sequence e (j) of adjacent point-to-point transmission, e (j)=fN
(k+1)-fN(k), liquid level height value corresponding to the sequence and using following liquid level reduction formula is calculated, obtains total nN+t-
1 liquid level L (j):
Wherein T is the temperature of environment, and t is the total number of leak source, and Δ L is compensation term.
2. the acoustic resonance level measuring method according to claim 1 based on iteration Weighted Fusion, it is characterised in that: work as liquid
When position height is larger, sound wave experienced the phase offset of zero degree in every secondary reflection, therefore there is certain " phase offset ", shadow
Ring measurement result;That is, need to compensate error caused by this factor, using segmented compensation function, liquid level is divided into three sections,
Each section of measurement height increases certain compensation rate Δ L
L (j) summation of acquirement is averaged, final liquid level L is obtained
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810901827.8A CN108981857B (en) | 2018-08-09 | 2018-08-09 | Acoustic resonance liquid level measurement method based on iterative weighted fusion |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810901827.8A CN108981857B (en) | 2018-08-09 | 2018-08-09 | Acoustic resonance liquid level measurement method based on iterative weighted fusion |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108981857A true CN108981857A (en) | 2018-12-11 |
CN108981857B CN108981857B (en) | 2020-03-10 |
Family
ID=64556271
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810901827.8A Active CN108981857B (en) | 2018-08-09 | 2018-08-09 | Acoustic resonance liquid level measurement method based on iterative weighted fusion |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108981857B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110260948A (en) * | 2019-06-20 | 2019-09-20 | 杭州电子科技大学 | Level measuring method based on acoustic resonance frequency nonlinear filtering |
CN113432674A (en) * | 2021-06-30 | 2021-09-24 | 杭州艾科赛德电子科技有限公司 | Container liquid level detection system based on sonar |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB776526A (en) * | 1954-12-10 | 1957-06-05 | Jack Kritz | Acoustic flowmeter |
CN101852638A (en) * | 2010-05-18 | 2010-10-06 | 杭州电子科技大学 | Liquid level measurement method based on resonance frequency of sound wave on fixed frequency range |
CN104154955A (en) * | 2014-05-19 | 2014-11-19 | 北京理工大学 | Method and system for dynamic measurement of liquid surface morphology and dosage of liquid propellant in storage tank |
CN106643982A (en) * | 2016-11-30 | 2017-05-10 | 清华大学 | Liquid level measuring method based on acoustic wave resonance frequency |
CN107741263A (en) * | 2017-10-13 | 2018-02-27 | 杭州电子科技大学 | Liquid level method of estimation based on frequency domain information fusion |
CN108304617A (en) * | 2018-01-05 | 2018-07-20 | 哈尔滨工程大学 | Buoyant raft structure wideband line spectrum vibration noise Fast Prediction method |
-
2018
- 2018-08-09 CN CN201810901827.8A patent/CN108981857B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB776526A (en) * | 1954-12-10 | 1957-06-05 | Jack Kritz | Acoustic flowmeter |
CN101852638A (en) * | 2010-05-18 | 2010-10-06 | 杭州电子科技大学 | Liquid level measurement method based on resonance frequency of sound wave on fixed frequency range |
CN104154955A (en) * | 2014-05-19 | 2014-11-19 | 北京理工大学 | Method and system for dynamic measurement of liquid surface morphology and dosage of liquid propellant in storage tank |
CN106643982A (en) * | 2016-11-30 | 2017-05-10 | 清华大学 | Liquid level measuring method based on acoustic wave resonance frequency |
CN107741263A (en) * | 2017-10-13 | 2018-02-27 | 杭州电子科技大学 | Liquid level method of estimation based on frequency domain information fusion |
CN108304617A (en) * | 2018-01-05 | 2018-07-20 | 哈尔滨工程大学 | Buoyant raft structure wideband line spectrum vibration noise Fast Prediction method |
Non-Patent Citations (1)
Title |
---|
徐晓滨等: "基于声波共振的液位测量实验系统设计与测试", 《实验技术与管理》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110260948A (en) * | 2019-06-20 | 2019-09-20 | 杭州电子科技大学 | Level measuring method based on acoustic resonance frequency nonlinear filtering |
CN113432674A (en) * | 2021-06-30 | 2021-09-24 | 杭州艾科赛德电子科技有限公司 | Container liquid level detection system based on sonar |
CN113432674B (en) * | 2021-06-30 | 2023-09-26 | 杭州艾科赛德电子科技有限公司 | Container liquid level detection system based on sonar |
Also Published As
Publication number | Publication date |
---|---|
CN108981857B (en) | 2020-03-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2352930C1 (en) | Method for reduction of digitised data in probe emat - "salamander" | |
CN101852638B (en) | Liquid level measurement method based on resonance frequency of sound wave on fixed frequency range | |
CN102563361A (en) | Device and method for detecting and positioning leakage of gas transmission pipeline based on conventional data and sound wave signals | |
CN111980676B (en) | Method and processing device for evaluating well cementation quality by array acoustic logging | |
CN104897780A (en) | Method for positioning acoustic emission source by using acoustic emission signal energy | |
CN108981857A (en) | Acoustic resonance level measuring method based on iteration Weighted Fusion | |
RU2017139770A (en) | INSTALLATION AND METHOD OF REMOTE MEASUREMENT OF GEOMETRIC PARAMETERS OF THE PIPELINE AT THE STAGE OF THE REMOVAL BY MEANS OF SOUND WAVES IN THE MODE OF REAL TIME | |
CN110469782A (en) | A kind of pipe leakage positioning device based on sef-adapting filter | |
CN103321633A (en) | Method and device for detecting working fluid level depth of oil well | |
CN109991590A (en) | A kind of System and method for of the testing transducer low frequencies characteristic in confined space pressurized tank | |
CN106770668B (en) | Method for detecting quality of single-hole foundation pile by acoustic transmission method | |
CN105180853A (en) | Electromagnetic ultrasonic metal material thickness measurement method | |
CA Oliveira et al. | Three-dimensional global scale underwater sound modeling: The T-phase wave propagation of a Southern Mid-Atlantic Ridge earthquake | |
CN205352653U (en) | Accurate ranging system of supersound in low speed wind tunnel | |
CN107741263B (en) | Liquid level estimation method based on frequency domain information fusion | |
CN103389153A (en) | Ultrasonic wave transit time measuring circuit by secondary reflection waves | |
CN113533504A (en) | Sub-surface crack quantitative measurement method based on laser ultrasonic surface wave frequency domain parameters | |
CN110987318B (en) | Automatic detection device and detection method for gas leakage of high-pressure pipeline | |
CN109270172A (en) | Verify the method and device of ultrasonic water immersion piezoelectric probe | |
JP2004061361A (en) | Piping breakage investigating apparatus | |
CN112362732A (en) | Free-form surface defect detection method and system based on laser ultrasonic elliptical imaging | |
CN104501909B (en) | A kind of small-range liquid level emasuring device and measuring method based on ultrasonic wave | |
CN107064995B (en) | A kind of acoustics method for positioning underground pipeline based on DSP | |
CN110260948A (en) | Level measuring method based on acoustic resonance frequency nonlinear filtering | |
CN106195648B (en) | A kind of experimental test procedures of the equivalent pipe range of reducer pipe |
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 |