CN104457918B - 一种回水流量检测偏差修正方法 - Google Patents
一种回水流量检测偏差修正方法 Download PDFInfo
- Publication number
- CN104457918B CN104457918B CN201410728525.7A CN201410728525A CN104457918B CN 104457918 B CN104457918 B CN 104457918B CN 201410728525 A CN201410728525 A CN 201410728525A CN 104457918 B CN104457918 B CN 104457918B
- Authority
- CN
- China
- Prior art keywords
- pipeline
- data record
- wall thickness
- flowmeter
- caliber
- 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.)
- Active
Links
Landscapes
- Measuring Volume Flow (AREA)
Abstract
本发明属于矿山、选矿及化工流程生产过程中工业循环利用回水的流量检测和修正,主要用于解决由于流量计和流体管道内壁慢时变黏结所产生的流量检测值反向漂移问题,并提出一种工业回水供水管道内壁黏结所造成慢时变缩径问题的快速修正解决办法,提高回水生产流程的流量检测精度,给出管道壁厚在线测量方法,根据不同介质声波传导率的不同采用便携式测厚仪测量管道同一截面的壁厚,由各点壁厚推算回水管道的有效内径,提出一种回水检测流量慢时变漂移快速修正算法,为复杂条件下回水流量的自动控制成为可能,该算法方便简单易实现。
Description
技术领域
本发明涉及矿山、选矿及化工流程生产过程中工业循环利用回水的流量检测和修正,主要用于解决由于流量计和流体管道内壁慢时变黏结所产生的流量检测值反向漂移问题。
背景技术
为实现资源循环利用,减少新水消耗量,大多数矿山、选矿和化工等高能耗企业将生产过程排出废水在净化池中进行中和、沉淀后,净化池底流经渣浆泵打回堆存库或其它处理流程,回水(净化池中的上清液)则返回生产流程继续循环利用。为降低生产成本,工业生产回水一般不进行化学处理,只进行简单的絮凝、沉淀过程,经常出现回水固体杂质含量高,甚至浑浊现象,尤其是回水分界面附近最为明显。
为减少检测仪表的维护量、提高使用寿命,固液两相流体的流量检测一般选用电磁流量计。测量介质的流通面大是电磁流量计的特点,但同时电磁流量计的正、负和参比电极长期与回水接触,尤其是低压回水管路易发生黏结现象,致使测量电极之间的电导率发生缓慢变化,无法稳定可靠地准确计量流量。为消除电磁流量计电极黏结所产生的流量检测偏差,一般采用自清洗电磁流量计,即按固定周期产生高频振荡,依靠流量计内腔测量介质的振动冲刷去除正、负电极和参比电极表面的黏结物,保证在测量介质特性不变的条件下,电极电导率的相对稳定。
在电磁流量计的允许测量范围内,根据液固两相流体流量与流速和管径截面成正比的测量原理,在测量介质和流速一定(电极电导率不变)的情况下,流量大小与管径成正比。由于带有腐蚀性的工业回水杂质含量高,黏度大,杂质极易在流量计内壁沉积结垢,这种慢时变结垢的直接影响就是流量计及其供水管线的管径变细。随着时间推移,流量检测数据与实际值相比偏差越来越大,最终导致无法为生产作业提供指导,更谈不上控制回路的自动投运。
发明内容
为克服上述现有回水流量检测技术的不足,实现对工业生产过程回水流量的准确检测,方便自动控制回路的投运,本发明提出一种回水流量检测快速修正方法。回水流量快速修正方法包括自清洗电磁流量计、便携式测厚仪、上位监控机、修正模型、显示界面等四部分组成。回水流量快速修正方法要解决的问题如下:
(1)为防止计算精度对不同类数据造成的影响,消除大数吃小数现象,对辨识模型所需的定长数据进行归一化处理;
(2)利用便携式测厚仪定期对流量计同一截面上的上、下、左、右四点进行厚度测量,折算后转换为流量计的有效管径,形成由时间尺度、有效管径值所构成的二维定长队列,应用最小二乘非线性拟合算法,在线辨识管径记忆因子等参数;
(3)开发一种回水流量计内壁慢时变快速修正模型,根据流量计过程检测值、管径修正值、记忆因子、投运时间等参数,在线快速计算流量计管径;
(4)利用组态软件的上位脚本工具,编制数据输入窗口、数据处理、参数辨识和流量修正算法模型;
(5)监控界面显示流量修正,为回水流量自动调节提供条件。
为实现上述功能,本发明所采用的技术方案如下:
(1)在不同时间段内用便携式测厚仪对流量计或回水管道同一截面上、下、左、右四点的壁厚进行多点测量,将测量时间、测量壁厚形成一条数据记录,推入定长队列中,并对每条数据记录中的四点壁厚值进行处理,最终形成由测量时间尺度、流量计管径所构成的二维数据记录集; (1)
(2)
(3)
(4)
(5)
式中,、是第条壁厚测量数据记录相对时间,单位为秒/s;
是第条管道上部壁厚测量值,单位毫米/;
是第条管道下部壁厚测量值,单位毫米/;
是第条管道左侧壁厚测量值,单位毫米/;
是第条管道右侧壁厚测量值,单位毫米/;
是管道测量壁厚数据记录集;
是管道测量壁厚数据记录元素;
是变换后管道测量壁厚数据记录集;
或变换后管道壁厚元素集;
是第条壁厚记录测量间隔时间,单位为天;
是第条管道等效壁厚,单位毫米/;
是定长队列的数据记录长度;
是管道有效内径数据记录集;
或变换后管道内径元素集;
是第条壁厚记录对应的管道有效流通内径,单位毫米/;
是管道起始流通内径,单位毫米/;
是管道起始壁厚,单位毫米/。
(2)对二维定长队列数据记录集进行归一化处理,依据数据集内流量计管径测量值符合独立正态分布条件和3σ误差限原理剔除异常数据,即剔除大于和小于的记录,消除异常数据对辨识模型的影响;
; (6)
(7)
(8)
(9)
(10)
(11)
根据上式(7)、(8)、(9)、(10)、(11)即可得出;
同理, (12)
根据上式(7)、(8)、(9)、(10)、(12)即可得出;
由数据记录剔除限、可得,模型辨识数据记录集。
(13)
(3)建立流量计管径参数的最小二乘辨识模型,根据数据记录集计算记忆因子;
(14)
是归一化管道有效内径数据记录集;
或归一化管道有效内径数据记录元素集;
是第条管道有效流通内径归一化标准值;
是第条管道有效流通换算内径间隔时间,单位为天;
是第条壁厚记录对应的管道有效流通内径,单位毫米/;
是管道有效流通内径最小值,单位毫米/;
是管道有效流通内径标称值,单位毫米/;
是条归一化管道有效内径对应的平均值;
是条归一化管道有效内径标准差;
是归一化管道有效内径的下限;
是归一化管道有效内径的上限;
和是单因素独立正态分布概率密度函数;
是归一化管道有效内径数据记录集中小于下限的概率;
是归一化管道有效内径数据记录集中大于上限的概率;
是大于下限,小于上限的归一化管道有效内径数据记录集;
是集合内的元素记录数;
是管道标称内管径;
是管道归一化后管径值;
是流量计慢时变缩径记忆因子;
是辨识模型的系数;
=2.71828,自然数;
是流量计相邻测量间隔时间(天)。
对式(14)进行变换可得,
(15)
令,式(15)两边取对数
(16)
和是中间参数;
是自然对数。
对式(16)中参数引入变量,,则
(17)
(18)
将模型辨识数据记录集带入式(16)即可得到系数,。
,、和是中间参数;
是集合内的元素记录数;
是对求和;
是对求和;
是对平方求和;
是对的乘积求和。
(4)反变换系数,,即模型系数和记忆因子,,带入式(14)即可得到辨识模型,经反归一化变换后,在线修正由于回水黏结所导致的管径慢时变缩径;
(19)
是流量计当前管道内径值;
是流量计归一化后管道内径值;
流量计折算管径记录集内最大值;
流量计折算管径记录集内最小值。
(5)根据流量计垂直流通面积与流量成正比原理,将辨识模型管径修正值带入修正模型,实现流量在线修正和集中显示;
(20)
是修正后流量值;
是流通管道标称管径值;
是流量计修正前实际流量值;
与现有技术相比,本发明的有益效果具体如下:
(1)提出一种工业回水供水管道内壁黏结所造成慢时变缩径问题的快速修正解决办法,提高回水生产流程的流量检测精度;
(2)给出管道壁厚在线测量方法,根据不同介质声波传导率的不同采用便携式测厚仪测量管道同一截面的壁厚,由各点壁厚推算回水管道的有效内径;
(3)提出一种回水检测流量慢时变漂移快速修正算法,为复杂条件下回水流量的自动控制成为可能,该算法方便简单易实现。
附图说明:
附图1为回水流量快速修正方法信号传输示意图;
附图2为数据处理流程。
图中:自清洗电磁流量调节阀1、流量计2、IO模块3、PLC处理器4、阀位调节信号5、流量检测信号6、调节阀阀位反馈信号7、上位监控机8、便携式超声波测厚仪9、回水输送管道10。
具体实施方式
下面结合附图对本发明进一步说明。
如附图1所示,其中自清洗电磁流量调节阀1和自清洗电磁流量计2安装在回水输送管道10上,其中自清洗电磁流量调节阀1的调节阀阀位反馈信号7和自清洗电磁流量计2的过程检测流量通过IO模块3将信号接入PLC系统4,通过通讯总线在上位监控机8中进行集中监控。
为实现慢时变缩径所导致的回水流量偏差,采用便携式超声波测厚仪9对流量计的同一垂直截面的上、下、左、右壁厚进行检测。为保证测量准确性,首先对管道表面的涂层进行打磨,然后均匀喷涂耦合剂,消除气泡对测量精度的影响。测量数据通过便携式超声波测厚仪9专用USB数据线与上位监控机8进行连接,导入到上位监控机8中手动删除异常或错误数据,形成如附图2左侧所示的数据记录,再对同一截面的4条测量记录求平均值,进行修正后得到如附图2所示的数据集合,并以数据文件的形式进行存储,文件名即是年月日时分秒。
首先从最新数据文件中读取辨识模型的数据记录,进行从式(1)至式(11)所示的数据记录变换、归一化和数据处理,形成模型辨识数据集。采用变换方法,通过式(12)至式(17)计算,将流量计慢时变缩径非线性模型转化为线性模型,通过线性最小二乘算法辨识出模型系数和记忆因子,最后再反归一化计算出适合在线快速修正的流量计内径计算模型。完成内径修正后,根据式(18)实时修正检测流量值,在监控界面中显示,并可以作为流量、调节阀开度调节控制回路的过程值。
Claims (6)
1.一种回水流量检测偏差修正方法,其特征在于:所述的该方法步骤如下:
1)在不同时间段内用便携式测厚仪对流量计或回水管道同一截面上、下、左、右四点的壁厚进行多点测量,将测量时间、测量壁厚形成一条数据记录,推入定长队列中:
,并对每条数据记录中的四点壁厚值进行处理,最终形成由测量时间尺度、流量计管径所构成的二维数据记录集;
2)根据步骤1)得到的二维数据记录集,对二维数据记录集进行归一化处理,依据二维数据记录集内流量计管径测量值符合独立正态分布条件和3σ误差限原理剔除异常数据,剔除异常数据为大于PU和小于PL的记录,消除异常数据对辨识模型的影响,并得到模型辨识数据记录集,其中PU是归一化管道有效内径数据记录集中大于上限的概率,PL是归一化管道有效内径数据记录集中小于下限的概率;
3)根据步骤2)中消除异常数据的二维数据记录集,建立流量计管径参数的最小二乘辨识模型,根据数据记录集计算记忆因子λ;
4)根据步骤3),计算得到辨识模型,经反归一化变换后,在线修正由于回水黏结所导致的管径慢时变缩径,得到辨识模型管径修正值D;
5)根据流量计管径与流量成正比原理,将步骤4)中得到的辨识模型管径修正值D带入修正模型,实现流量在线修正和集中显示。
2.如权利要求1所述的一种回水流量检测偏差修正方法,其特征在于:所述步骤1)中得到由测量时间尺度、流量计管径所构成的二维数据记录集过程如下:
; (1)
; (2)
; (3)
; (4)
; (5)
式中,、是第条壁厚测量数据记录相对时间,单位为秒/s;
是第条管道上部壁厚测量值,单位毫米/;
是第条管道下部壁厚测量值,单位毫米/;
是第条管道左侧壁厚测量值,单位毫米/;
是第条管道右侧壁厚测量值,单位毫米/;
是管道测量壁厚数据记录集;
是管道测量壁厚数据记录元素;
是变换后管道测量壁厚数据记录集;
或变换后管道壁厚元素集;
是第条壁厚记录测量间隔时间,单位为天;
是第条管道等效壁厚,单位毫米/;
是定长队列的数据记录长度;
是管道有效内径数据记录集;
或变换后管道内径元素集;
是第条壁厚记录对应的管道有效流通内径,单位毫米/;
是管道起始流通内径,单位毫米/;
是管道起始壁厚,单位毫米/。
3.如权利要求1所述的一种回水流量检测偏差修正方法,其特征在于:所述步骤2)中得到的模型辨识数据记录集过程如下:
; (6)
; (7)
; (8)
; (9)
; (10)
; (11)
根据上式(5)、(6)、(7)、(8)即可得出DL;
同理,; (12)
根据上式(5、(6)、(7)、(9)即可得出DU;
由数据记录剔除限DL、DU可得,模型辨识数据记录集,
; (13)
; (14)
是归一化管道有效内径数据记录集;
或归一化管道有效内径数据记录元素集;
是第条管道有效流通内径归一化标准值;
是第条管道有效流通换算内径间隔时间,单位为天;
是第条壁厚记录对应的管道有效流通内径,单位毫米/;
是管道有效流通内径最小值,单位毫米/;
是管道有效流通内径标称值,单位毫米/;
是条归一化管道有效内径对应的平均值;
是条归一化管道有效内径标准差;
是归一化管道有效内径的下限;
是归一化管道有效内径的上限;
和是单因素独立正态分布概率密度函数;
是归一化管道有效内径数据记录集中小于下限的概率;
是归一化管道有效内径数据记录集中大于上限的概率;
是大于下限,小于上限的归一化管道有效内径数据记录集;
是集合内的元素记录数;
是管道标称内管径;
是管道归一化后管径值;
是流量计慢时变缩径记忆因子;
是辨识模型的系数;
=2.71828,自然数;
是流量计相邻测量间隔时间(天)。
4.如权利要求1所述的一种回水流量检测偏差修正方法,其特征在于:所述步骤3)中得到最小二乘辨识模型,根据模型辨识数据记录集计算记忆因子λ的过程如下:
对式(14)进行变换可得,
; (15)
令,式(15)两边取对数
; (16)
和是中间参数;
是自然对数;
对式(16)中参数引入变量,,则
; (17)
; (18)
将模型辨识数据记录集带入式(16)即可得到反变换系数,;
,、反变换系数和是中间参数;
是集合内的元素记录数;
是对求和;
是对求和;
是对平方求和;
是对的乘积求和。
5.如权利要求4所述的一种回水流量检测偏差修正方法,其特征在于:所述的反变换系数,,即模型系数和记忆因子,,带入式(14)即可得到辨识模型,经反归一化变换后,在线修正由于回水黏结所导致的管径慢时变缩径;
; (19)
是流量计当前管道内径值;
是流量计归一化后管道内径值;
流量计折算管径记录集内最大值;
流量计折算管径记录集内最小值。
6.如权利要求5所述的一种回水流量检测偏差修正方法,其特征在于:所述的根据流量计垂直流通面积与流量成正比原理,将辨识模型管径修正值带入修正模型,实现流量在线修正和集中显示;
; (20)
是修正后流量值;
是流通管道标称管径值;
是流量计修正前实际流量值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410728525.7A CN104457918B (zh) | 2014-12-03 | 2014-12-03 | 一种回水流量检测偏差修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410728525.7A CN104457918B (zh) | 2014-12-03 | 2014-12-03 | 一种回水流量检测偏差修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104457918A CN104457918A (zh) | 2015-03-25 |
CN104457918B true CN104457918B (zh) | 2017-12-01 |
Family
ID=52904342
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410728525.7A Active CN104457918B (zh) | 2014-12-03 | 2014-12-03 | 一种回水流量检测偏差修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104457918B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109781215B (zh) * | 2019-01-17 | 2020-07-17 | 合肥工业大学 | 一种电磁流量传感器多场耦合感应电动势的数值仿真方法 |
CN112597431A (zh) * | 2020-12-24 | 2021-04-02 | 金川集团股份有限公司 | 一种基于积分扩散因子的流量插值积算方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000504411A (ja) * | 1996-01-18 | 2000-04-11 | イエダ リサーチ アンド デベロツプメント カンパニー リミテツド | 流体が通過するシステムをモニタする装置 |
WO2005088262A2 (en) * | 2004-03-10 | 2005-09-22 | Cidra Corporation | Method and apparatus for measuring parameters of a stratified flow |
CN102818659A (zh) * | 2011-06-09 | 2012-12-12 | 沈阳航发热计量技术有限公司 | 一种供热失水测量系统及测量方法 |
CN103074873A (zh) * | 2013-02-07 | 2013-05-01 | 天津大学 | 一种非理想流场下的多声道超声流量计声道布置方法 |
-
2014
- 2014-12-03 CN CN201410728525.7A patent/CN104457918B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000504411A (ja) * | 1996-01-18 | 2000-04-11 | イエダ リサーチ アンド デベロツプメント カンパニー リミテツド | 流体が通過するシステムをモニタする装置 |
WO2005088262A2 (en) * | 2004-03-10 | 2005-09-22 | Cidra Corporation | Method and apparatus for measuring parameters of a stratified flow |
CN102818659A (zh) * | 2011-06-09 | 2012-12-12 | 沈阳航发热计量技术有限公司 | 一种供热失水测量系统及测量方法 |
CN103074873A (zh) * | 2013-02-07 | 2013-05-01 | 天津大学 | 一种非理想流场下的多声道超声流量计声道布置方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104457918A (zh) | 2015-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104793638B (zh) | 打叶风分智能追踪控制方法 | |
CN105403265A (zh) | 一种自动校正零点漂移的超声水表及其校正方法 | |
CN108151835A (zh) | 一种基于多参数检测的泵站信息智能监测终端的率定方法 | |
CN104457918B (zh) | 一种回水流量检测偏差修正方法 | |
CN104502161B (zh) | 粉尘采样器检定装置 | |
CN102538914A (zh) | 一种具有校正功能的电子显示水表 | |
CN109931801B (zh) | 一种火电厂循环冷却水塔池自动补水系统及补水方法 | |
CN205785444U (zh) | 一种超声波液位计校验系统 | |
CN104034377A (zh) | 一种空气压缩机流量测量装置及其标定系统和标定方法 | |
CN110231069A (zh) | 一种冷却塔水量损失的精确测量方法 | |
CN112580183B (zh) | 一种在线学习水泵模型实时流量精确控制方法 | |
CN108709613A (zh) | 一种可实现多流量计同时校准的系统和方法 | |
CN106404091A (zh) | 一种基于体积法的泵站流量测法 | |
CN108106679A (zh) | 一种电站磨煤机入口风量的测量方法及系统 | |
CN105571661B (zh) | 基于示值误差拟合的电磁流量计仪表特征系数的计算方法 | |
CN106444931A (zh) | 白酒酒精度在线监测的方法及监测设备 | |
CN211061457U (zh) | 一种高压直流输电换流阀冷却水系统集中式水质监测装置 | |
CN207717197U (zh) | 流量计的校准装置 | |
CN106679773A (zh) | 基于分段线性的电磁流量计仪表特征系数计算方法 | |
CN207082035U (zh) | 一种凝结水精处理控制系统 | |
CN101482429A (zh) | 一种检定弯管流量计的装置 | |
CN203502327U (zh) | 一种火电机组热力系统腐蚀在线监测装置及其系统 | |
CN216997726U (zh) | 烧结球团生产用水水质调节系统 | |
CN217359031U (zh) | 一种电热水器检测系统 | |
CN205940718U (zh) | 一种智能活塞体积管检定装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |