CN105973349B - A kind of on-line prediction method of Coriolis flowmeter intrinsic frequency - Google Patents
A kind of on-line prediction method of Coriolis flowmeter intrinsic frequency Download PDFInfo
- Publication number
- CN105973349B CN105973349B CN201610267912.4A CN201610267912A CN105973349B CN 105973349 B CN105973349 B CN 105973349B CN 201610267912 A CN201610267912 A CN 201610267912A CN 105973349 B CN105973349 B CN 105973349B
- Authority
- CN
- China
- Prior art keywords
- intrinsic frequency
- frequency
- signal
- coriolis flowmeter
- amplitude
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 13
- 238000005259 measurement Methods 0.000 claims abstract description 20
- 238000005070 sampling Methods 0.000 claims abstract description 11
- 230000009466 transformation Effects 0.000 claims abstract description 11
- 238000001228 spectrum Methods 0.000 claims description 4
- 230000007704 transition Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims 2
- 230000006641 stabilisation Effects 0.000 claims 1
- 238000011105 stabilization Methods 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 4
- 230000005526 G1 to G0 transition Effects 0.000 abstract description 3
- 238000000691 measurement method Methods 0.000 abstract description 2
- 238000012545 processing Methods 0.000 abstract description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 238000013016 damping Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012857 repacking Methods 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F25/00—Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
The present invention relates to a kind of on-line prediction methods of Coriolis flowmeter intrinsic frequency.Mainly solve the problems, such as that existing Coriolis flowmeter intrinsic frequency measurement method requires to stop just being capable of measuring by flowmeter.It is characterized by: detecting electromagnet output signal in real-time sampling Coriolis flowmeter measurement pipe, when signal amplitude changes more than 1/10, changes the stationary phase data in preceding ten periods using amplitude, calculate original natural frequency value f0, intrinsic frequency variation delta f is calculated, and then pass through f=f based on amplitude variation phase signal characteristic parameter identification using the sampled data in amplitude variation phase in continuous two ten periods0+ Δ f calculates the current intrinsic frequency f of measurement pipe.The variation of measurement pipe signal frequency is found out by the data processing of Fourier transformation by the detection signal for sampling multiple periods, goes to calculate current intrinsic frequency in real time, the normal metering without influencing Coriolis flowmeter.
Description
Technical field
The present invention relates to a kind of prediction techniques of flowmeter intrinsic frequency, and in particular to a kind of Coriolis flowmeter intrinsic frequency
On-line prediction method.
Background technique
Coriolis flowmeter CMF (hereinafter referred to as Coriolis flowmeter) is a kind of high-precision mass flowmenter, is answered extensively
For the measurement of such as oil product, chemical medium, papermaking black liquid and slurry medium in process industrial field, using fluid in vibration
The generated Coriolis force directly proportional to mass flow when flowing in pipeline detects measurement pipe or more two-way sinusoidal signal
Phase difference directly calculates mass flow.Wherein, the amplitude stability of vibration control system is before guaranteeing flow measurement accuracy of measurement
Propose condition.In general, the driving frequency of measurement pipe should be made to match with intrinsic frequency, if can know in real time to make its amplitude stability
The design parameter of intrinsic frequency, the vibration control system to Coriolis flowmeter is very helpful.
In digital Coriolis flowmeter, to guarantee amplitude stability, real-time measurement is needed to track intrinsic frequency, it is existing solid
Having frequency test method, including free damping damped method etc. to be required to Coriolis flowmeter stopping can just measure, and cannot achieve
The intrinsic frequency on-line prediction of flowmeter.
Summary of the invention
In order to overcome the shortcomings of that background technique, the present invention provide a kind of on-line prediction side of Coriolis flowmeter intrinsic frequency
Method mainly solves existing Coriolis flowmeter intrinsic frequency measurement method and requires to ask what flowmeter stopping being just capable of measuring
Topic.
The technical scheme adopted by the invention is that: a kind of on-line prediction method of Coriolis flowmeter intrinsic frequency is adopted in real time
Electromagnet output signal is detected in sample Coriolis flowmeter measurement pipe, when signal amplitude changes more than 1/10, is changed using amplitude
The stationary phase data in preceding ten periods calculate original natural frequency value f0, utilize amplitude variation phase in continuous two ten periods
Sampled data calculates intrinsic frequency variation delta f, and then pass through f=f based on amplitude variation phase signal characteristic parameter identification0+
Δ f calculates the current intrinsic frequency f of measurement pipe.
Number of cycles sampling is carried out to amplitude stability phase signal, it, will be corresponding when power spectrum maximum by spectrum analysis
Frequency, regard as the fundamental frequency f of signal0, i.e., former intrinsic frequency f0。
Number of cycles twice is carried out to magnitude transition phase signal to sample, and reflection is obtained by Fourier transformation and is sampled for the first time
The real part a1 and imaginary part b1 parameter of signal characteristic and real part a2 and imaginary part the b2 parameter of second of sampled signal feature, in turn
Pass through
Calculate intrinsic frequency variation delta f.
Signal characteristic parameter real part a1, a2 and imaginary part b1, b2 are obtained by Fourier transformation.
The beneficial effects of the present invention are: by the detection signal for sampling multiple periods, at the data of Fourier transformation
Reason, finds out the variation of measurement pipe signal frequency, goes to calculate current intrinsic frequency in real time, without influencing Coriolis flowmeter just
Often metering.
Specific embodiment
In the vibration control system of Coriolis flowmeter, using a kind of Frequency Measurement Algorithm based on Fourier, when frequency changes
When, the variation of frequency is found out according to the phase angle difference between two adjacent data windows, by sampling the detection signal in multiple periods, is passed through
The data processing of Fourier transformation finds out the variation of measurement pipe signal frequency, goes to calculate current intrinsic frequency in real time, without
Influence the normal metering of Coriolis flowmeter.
Number of cycles sampling is carried out to amplitude stability phase signal, it, will be corresponding when power spectrum maximum by spectrum analysis
Frequency, regard as the fundamental frequency f of signal0, i.e., former intrinsic frequency f0。
Intrinsic frequency f0, then its corresponding calculating angular frequency is denoted as ω0.Meanwhile the true intrinsic frequency of measurement pipe being denoted as
F, corresponding true angular frequency are denoted as ω.Due to f0Although be not with f it is stringent equal, its value is not much different, therefore, can be with
Make following reasonable hypothesis:
1, intrinsic frequency f=f0+ Δ f, Δ f indicate intrinsic frequency and calculate the difference of frequency, and Δ f value very little；
2, the calculating angular frequency being calculated0It is fluctuated near true angular frequency；
The derived object of following mathematical formulae exactly acquires the value of Δ f:
The signal that intrinsic frequency indicates are as follows:
F (t)=Fsin (2 π f0t+2πΔft+α0) (1)
Wherein F and α0Indicate amplitude and initial phase.
The resulting sample series of a cycle will be sampled and carry out Fourier transformation, obtained:
Due to ω ≈ ω0, so formula (2) is approximately equal to following formula:
(1) is substituted into (3), and is further spread out, available:
Equally, the sample series in second period carry out same Fourier transformation, obtain:
(1) is substituted into (5), same expansion is done, obtains:
It is available by (4) and (6):
It eliminates in (7) and (8), it obtains:
It can thus be concluded that:
Electromagnet output signal is detected in real-time sampling Coriolis flowmeter measurement pipe, when signal amplitude changes more than 1/10
When, change the stationary phase data in preceding ten periods using amplitude, calculates original natural frequency value f0, connected using the amplitude variation phase
The sampled data in continuous two ten periods is reflected based on amplitude variation phase signal characteristic parameter identification by Fourier transformation
The real part a1 and imaginary part b1 parameter of first time sampled signal feature and the real part a2 and imaginary part b2 of second of sampled signal feature
Parameter, and then pass through
Calculate intrinsic frequency variation delta f.Then pass through f=f0+ Δ F calculates the current intrinsic frequency f of measurement pipe.
The beneficial effects of the present invention are: by the detection signal for sampling multiple periods, at the data of Fourier transformation
Reason, finds out the variation of measurement pipe signal frequency, goes to calculate current intrinsic frequency in real time, without influencing Coriolis flowmeter just
Often metering.
Every technical staff's notice: of the invention although the present invention is described according to above-mentioned specific embodiment
Invention thought be not limited in the invention, any repacking with inventive concept will all be included in this patent protection of the patent right
In range.
Claims (2)
1. a kind of on-line prediction method of Coriolis flowmeter intrinsic frequency, it is characterised in that: the measurement of real-time sampling Coriolis flowmeter
Electromagnet output signal is detected on pipe, when signal amplitude changes more than 1/10, changes the stabilization in preceding ten periods using amplitude
Issue evidence calculates original natural frequency value f0, using the sampled data in amplitude variation phase in continuous two ten periods, it is based on width
It is worth the identification of variation phase signal characteristic parameter, calculates intrinsic frequency variation delta f, and then pass through f=f0+ Δ f calculates measurement pipe
Current intrinsic frequency f,
Original natural frequency value f0Calculation method: the sampling of number of cycles is carried out to amplitude stability phase signal, passes through frequency spectrum point
Analysis, corresponding frequency when by power spectrum maximum regard as the fundamental frequency f of signal0, i.e., former intrinsic frequency f0。
The calculation method of intrinsic frequency variation delta f: number of cycles twice is carried out to magnitude transition phase signal and is sampled, by Fu
Leaf transformation obtains the real part a of reflection first time sampled signal feature1With imaginary part b1Parameter and second sampled signal feature
Real part a2With imaginary part b2Parameter, and then pass through
Intrinsic frequency variation delta f is calculated, wherein T is the sampling period.
2. the on-line prediction method of Coriolis flowmeter intrinsic frequency according to claim 1, it is characterised in that: by Fu
Leaf transformation obtains signal characteristic parameter real part a1, a2With imaginary part b1, b2。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610267912.4A CN105973349B (en) | 2016-04-19 | 2016-04-19 | A kind of on-line prediction method of Coriolis flowmeter intrinsic frequency |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610267912.4A CN105973349B (en) | 2016-04-19 | 2016-04-19 | A kind of on-line prediction method of Coriolis flowmeter intrinsic frequency |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105973349A CN105973349A (en) | 2016-09-28 |
| CN105973349B true CN105973349B (en) | 2019-04-09 |
Family
ID=56993136
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610267912.4A Expired - Fee Related CN105973349B (en) | 2016-04-19 | 2016-04-19 | A kind of on-line prediction method of Coriolis flowmeter intrinsic frequency |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105973349B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113988144B (en) * | 2021-12-27 | 2022-03-08 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | Coriolis flowmeter natural frequency obtaining method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102095430A (en) * | 2010-11-18 | 2011-06-15 | 合肥工业大学 | Sensor dynamic error frequency-domain correction technology based on step response |
| CN102435849A (en) * | 2011-10-26 | 2012-05-02 | 深圳市科陆电子科技股份有限公司 | Method for improving frequency measurement precision by using digital filter |
| CN102506951A (en) * | 2011-10-28 | 2012-06-20 | 合肥工业大学 | Method and system of digital driving following for Coriolis mass flowmeter |
| CN103983849A (en) * | 2014-05-07 | 2014-08-13 | 江苏天浩达科技有限公司 | Real-time high-accuracy power harmonic analysis method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7343253B2 (en) * | 2005-07-11 | 2008-03-11 | Invensys Systems, Inc. | Coriolis mode processing techniques |
-
2016
- 2016-04-19 CN CN201610267912.4A patent/CN105973349B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102095430A (en) * | 2010-11-18 | 2011-06-15 | 合肥工业大学 | Sensor dynamic error frequency-domain correction technology based on step response |
| CN102435849A (en) * | 2011-10-26 | 2012-05-02 | 深圳市科陆电子科技股份有限公司 | Method for improving frequency measurement precision by using digital filter |
| CN102506951A (en) * | 2011-10-28 | 2012-06-20 | 合肥工业大学 | Method and system of digital driving following for Coriolis mass flowmeter |
| CN103983849A (en) * | 2014-05-07 | 2014-08-13 | 江苏天浩达科技有限公司 | Real-time high-accuracy power harmonic analysis method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105973349A (en) | 2016-09-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2683385C (en) | Wet gas measurement | |
| US10139259B2 (en) | System and method for metering gas based on amplitude and/or temporal characteristics of an electrical signal | |
| WO2014039873A1 (en) | Ultrasonic flow metering using compensated computed temperature | |
| WO2013002656A1 (en) | Flow rate determination method and apparatus | |
| CN207317905U (en) | Flowmeter calibrating installation | |
| Ghaoud et al. | Modelling and tracking a vortex flow-meter signal | |
| CN106092524A (en) | A kind of method using vibration signal accurately to extract tach signal | |
| CN106324340B (en) | A kind of method of synchronized phasor and frequency measurement dynamic property | |
| CN108333426A (en) | Power system frequency measurement method based on fourier algorithm | |
| CN103063913A (en) | Frequency tracking method for Fourier transform | |
| CN105973349B (en) | A kind of on-line prediction method of Coriolis flowmeter intrinsic frequency | |
| CN101241151B (en) | Real time frequency measurement method based on amplitude linear change model | |
| Zhang et al. | Mathematical model of time difference for Coriolis flow sensor output signals under gas-liquid two-phase flow | |
| CN112212951A (en) | Nonlinear correction method for in-situ calibration of permanent magnet sodium flowmeter based on signal frequency band selection | |
| Świsulski et al. | Methods of measurement signal acquisition from the rotational flow meter for frequency analysis | |
| Chen et al. | Research on a transit-time liquid ultrasonic flowmeter under unstable flow fields | |
| CN204514403U (en) | A kind of differential pressure mass flowmeter for vortex street | |
| Liu et al. | Signal processing of Coriolis mass flowmeters under gas-liquid two-phase flow conditions | |
| EP3312573B1 (en) | System and method for fluid metering | |
| CN104483549A (en) | Phase difference measuring method based on data extension and Hilbert transform | |
| CN114814275A (en) | Dynamic calculation method and device for multiphase fluid flow velocity | |
| EP3322962B1 (en) | Method for the measurement of an unsteady liquid flow rate, in particular of a high pressure liquid flow | |
| Xin et al. | Analysis and optimization of detecton curve of BLTD cementing mud flowmeter | |
| Song et al. | Study on detection and location of mobile oil pipeline leakage based on pressure gradient method | |
| Hu et al. | Research and implementation of amplitude ratio calculating system for coriolis mass flowmeter |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 325000 Wenzhou City National University Science Park incubator, No. 38 Dongfang South Road, Ouhai District, Wenzhou, Zhejiang Applicant after: Wenzhou Vocational College of Science and Technology Address before: 325000 No. 1000 Hongqiao Road, Ouhai District, Wenzhou, Zhejiang, 1000 Applicant before: Wenzhou Vocational College of Science and Technology |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190409 |