CN105973349B - 一种科氏流量计固有频率的在线预测方法 - Google Patents
一种科氏流量计固有频率的在线预测方法 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
本发明涉及一种科氏流量计固有频率的在线预测方法。主要解决了现有的科氏流量计固有频率测量方法都需要将流量计停止才能测量的问题。其特征在于:实时采样科氏流量计测量管上检测电磁铁输出信号,当信号幅值变化超过1/10时,利用幅值变化前十个周期的稳定期数据,计算出原有固有频率值f0,利用幅值变化期连续两个十周期的采样数据,基于幅值变化期信号特征参数识别,计算固有频率变化量Δf,进而通过f=f0+Δf计算出测量管当前固有频率f。通过采样多个周期的检测信号,经过傅里叶变换的数据处理,求出测量管信号频率的变化,去实时计算出当前的固有频率,而不影响科氏流量计的正常计量。
Description
技术领域
本发明涉及一种流量计固有频率的预测方法,具体涉及一种科氏流量计固有频率的在线预测方法。
背景技术
科里奥利流量计CMF(以下简称科氏流量计)是一种高精度的质量流量计,广泛应用于过程工业领域中如油品、化工介质、造纸墨液和浆体等介质的测量,利用流体在振动的管道中流动时所产生的与质量流量成正比的科里奥利力,检测测量管上下两路正弦信号的相位差,直接地计算质量流量。其中,振动控制***的振幅稳定是保证流量计测量精度的前提条件。通常,为使其振幅稳定,应使测量管的驱动频率与固有频率相吻合,若能实时知道固有频率的具体参数,将对科氏流量计的振动控制***非常有帮助。
在数字式科氏流量计中,为保证振幅稳定,需要实时测量跟踪固有频率,现有的固有频率测试方法,包括自由阻尼衰减法等均要求科氏流量计停止才能进行测量,无法实现流量计的固有频率在线预测。
发明内容
为了克服背景技术的不足,本发明提供一种科氏流量计固有频率的在线预测方法,主要解决了现有的科氏流量计固有频率测量方法都需要将流量计停止才能测量的问题。
本发明所采用的技术方案是:一种科氏流量计固有频率的在线预测方法,实时采样科氏流量计测量管上检测电磁铁输出信号,当信号幅值变化超过1/10时,利用幅值变化前十个周期的稳定期数据,计算出原有固有频率值f0,利用幅值变化期连续两个十周期的采样数据,基于幅值变化期信号特征参数识别,计算固有频率变化量Δf,进而通过f=f0+Δf计算出测量管当前固有频率f。
对振幅稳定期信号进行一次整数周期采样,通过频谱分析,将功率谱最大时对应的频率,认定为信号的基频f0,即原固有频率f0。
对振幅过渡期信号进行两次整数周期采样,通过傅里叶变换获得反映第一次采样信号特征的实部a1与虚部b1参数,以及第二次采样信号特征的实部a2与虚部b2参数,进而通过
计算出固有频率变化量Δf。
通过傅里叶变换获取信号特征参数实部a1,a2与虚部b1,b2。
本发明的有益效果是:通过采样多个周期的检测信号,经过傅里叶变换的数据处理,求出测量管信号频率的变化,去实时计算出当前的固有频率,而不影响科氏流量计的正常计量。
具体实施方式
在科氏流量计的振动控制***中,采用一种基于傅里叶的测频算法,当频率变化时,根据两个相邻数据窗间的相角差求出频率的变化,通过采样多个周期的检测信号,经过傅里叶变换的数据处理,求出测量管信号频率的变化,去实时计算出当前的固有频率,而不影响科氏流量计的正常计量。
对振幅稳定期信号进行一次整数周期采样,通过频谱分析,将功率谱最大时对应的频率,认定为信号的基频f0,即原固有频率f0。
固有频率f0,则其对应的计算角频率记为ω0。同时,将测量管真实的固有频率记为f,其对应的真实角频率记为ω。由于f0与f虽然不是严格相等,但其值相差不多,因此,可以做出如下合理的假设:
1、固有频率f=f0+Δf,Δf表示固有频率与计算频率的差值,并且Δf值很小;
2、计算得到的计算角频率ω0在真实角频率ω附近波动;
以下数学公式的推导目标,就是求得Δf的值:
固有频率表示的信号为:
f(t)=Fsin(2πf0t+2πΔft+α0) (1)
其中F和α0表示幅值和初始相位。
将采样一个周期所得的采样系列进行傅里叶变换,得:
由于ω≈ω0,所以公式(2)约等于以下公式:
将(1)代入(3),并进一步展开,可以得到:
同样,第二个的周期的采样系列进行同样的傅里叶变换,得:
将(1)代入(5),做同样的展开,得:
由(4)和(6)可以得到:
消除(7)和(8)中的项,得:
由此可得:
实时采样科氏流量计测量管上检测电磁铁输出信号,当信号幅值变化超过1/10时,利用幅值变化前十个周期的稳定期数据,计算出原有固有频率值f0,利用幅值变化期连续两个十周期的采样数据,基于幅值变化期信号特征参数识别,通过傅里叶变换获得反映第一次采样信号特征的实部a1与虚部b1参数,以及第二次采样信号特征的实部a2与虚部b2参数,进而通过
计算出固有频率变化量Δf。然后通过f=f0+ΔF计算出测量管当前固有频率f。
本发明的有益效果是:通过采样多个周期的检测信号,经过傅里叶变换的数据处理,求出测量管信号频率的变化,去实时计算出当前的固有频率,而不影响科氏流量计的正常计量。
各位技术人员须知:虽然本发明已按照上述具体实施方式做了描述,但是本发明的发明思想并不仅限于此发明,任何运用本发明思想的改装,都将纳入本专利专利权保护范围内。
Claims (2)
1.一种科氏流量计固有频率的在线预测方法,其特征在于:实时采样科氏流量计测量管上检测电磁铁输出信号,当信号幅值变化超过1/10时,利用幅值变化前十个周期的稳定期数据,计算出原有固有频率值f 0,利用幅值变化期连续两个十周期的采样数据,基于幅值变化期信号特征参数识别,计算固有频率变化量Δf,进而通过f=f0+Δf计算出测量管当前固有频率f,
原有固有频率值f0的计算方法:对振幅稳定期信号进行一次整数周期采样,通过频谱分析,将功率谱最大时对应的频率,认定为信号的基频f0,即原固有频率f0。
固有频率变化量Δf的计算方法:对振幅过渡期信号进行两次整数周期采样,通过傅里叶变换获得反映第一次采样信号特征的实部a1与虚部b1参数,以及第二次采样信号特征的实部a2与虚部b2参数,进而通过
计算出固有频率变化量Δf,其中T为采样周期。
2.根据权利要求1所述的科氏流量计固有频率的在线预测方法,其特征在于:通过傅里叶变换获取信号特征参数实部a1,a2与虚部b1,b2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610267912.4A CN105973349B (zh) | 2016-04-19 | 2016-04-19 | 一种科氏流量计固有频率的在线预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610267912.4A CN105973349B (zh) | 2016-04-19 | 2016-04-19 | 一种科氏流量计固有频率的在线预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105973349A CN105973349A (zh) | 2016-09-28 |
CN105973349B true CN105973349B (zh) | 2019-04-09 |
Family
ID=56993136
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610267912.4A Expired - Fee Related CN105973349B (zh) | 2016-04-19 | 2016-04-19 | 一种科氏流量计固有频率的在线预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105973349B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113988144B (zh) * | 2021-12-27 | 2022-03-08 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | 一种科氏流量计固有频率获得方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102095430A (zh) * | 2010-11-18 | 2011-06-15 | 合肥工业大学 | 基于阶跃响应的传感器动态误差频域修正技术 |
CN102435849A (zh) * | 2011-10-26 | 2012-05-02 | 深圳市科陆电子科技股份有限公司 | 使用数字滤波器提高频率测量精度的方法 |
CN102506951A (zh) * | 2011-10-28 | 2012-06-20 | 合肥工业大学 | 一种科氏质量流量计的数字驱动跟踪方法和*** |
CN103983849A (zh) * | 2014-05-07 | 2014-08-13 | 江苏天浩达科技有限公司 | 一种实时高精度的电力谐波分析方法 |
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/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102095430A (zh) * | 2010-11-18 | 2011-06-15 | 合肥工业大学 | 基于阶跃响应的传感器动态误差频域修正技术 |
CN102435849A (zh) * | 2011-10-26 | 2012-05-02 | 深圳市科陆电子科技股份有限公司 | 使用数字滤波器提高频率测量精度的方法 |
CN102506951A (zh) * | 2011-10-28 | 2012-06-20 | 合肥工业大学 | 一种科氏质量流量计的数字驱动跟踪方法和*** |
CN103983849A (zh) * | 2014-05-07 | 2014-08-13 | 江苏天浩达科技有限公司 | 一种实时高精度的电力谐波分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105973349A (zh) | 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 | |
EP2726827A1 (en) | Flow rate determination method and apparatus | |
Ghaoud et al. | Modelling and tracking a vortex flow-meter signal | |
CN106092524A (zh) | 一种使用振动信号精确提取转速信号的方法 | |
CN108333426A (zh) | 基于傅氏算法的电力***频率测量方法 | |
CN103063913A (zh) | 用于傅里叶变换的频率跟踪方法 | |
CN105973349B (zh) | 一种科氏流量计固有频率的在线预测方法 | |
Ravula et al. | Experimental validation of leak and water-ingression detection in low-pressure gas pipeline using pressure and flow measurements | |
CN101241151B (zh) | 基于幅值线性变化模型的实时频率测量方法 | |
Zhang et al. | Mathematical model of time difference for Coriolis flow sensor output signals under gas-liquid two-phase flow | |
Ś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 (zh) | 一种差压涡街质量流量计 | |
EP3312573B1 (en) | System and method for fluid metering | |
CN104483549A (zh) | 基于数据延拓和Hilbert变换的相位差测量方法 | |
EP3322962B1 (en) | Method for the measurement of an unsteady liquid flow rate, in particular of a high pressure liquid flow | |
Liu et al. | Signal processing of Coriolis mass flowmeters under gas-liquid two-phase flow conditions | |
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 | |
Zhu et al. | A Method of Detecting Measurement Uncertainty of Voltage Unbalance | |
CN204730893U (zh) | 一种集成迭代流量计 | |
Shen et al. | A novel method for CMF signal processing based on the revised sliding recursive DTFT algorithm | |
JP6318505B2 (ja) | 磁気式酸素分析方法及び磁気式酸素分析計 |
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 |