专利名称：流动体系的粘度检测方法
技术领域：
本发明涉及声波检测领域，尤其涉及搅拌釜、管道和环管中流体粘度的声波检测
方法。
背景技术：
粘度是流体的重要物理性质和技术指标之一，流体的粘性是流体微团间发生相对
滑移时产生切向阻力的性质，是流体抵抗剪切变形的能力或属性。粘度的测定在许多工业
部门和科学研究领域中都具有重要意义，特别是在石油化工、医药、冶金、食品等行业中。搅
拌釜、管道和环管在众多领域都有广泛的应用，对搅拌釜、管道和环管中流体粘度的测量是
控制生产流程、保证安全生产、控制与评定产品质量及科学研究的重要手段。 搅拌釜、管道和环管是工业生产中的核心设备，伴随着混合、溶解、反应等过程的
进行，不断发生着各种物理化学变化，从而使流体的粘度也随之而变。若能实现流体粘度
的实时在线检测，将能大大提高产品质量，并降低生产消耗。然而传统的流体粘度检测方
法，如旋转剪切速率降速式粘度测量法(ZL96109032. 4)、强迫共振法的流体粘度测量方法
(ZL200710030380. 3)等，都只能用于离线检测，无法应用于在线检测。而Brookf ield公司
和Cambridge公司的粘度检测仪(US5， 531， 102)，具有较高的检测精度，测量结果能够基
本满足很多工业过程的要求，然而，这些粘度检测仪及相关配套器件价格昂贵、维护保养复
杂，同时检测时还需要将检测仪插入到流体中，从而改变流体的流场分布，不利于流程的稳
定运行。 专利号为ZL01822778. 3的中国发明专利中提出了 "声探测器"的概念，通过将探 测器安装在与聚合反应器伸出的管子相联的线管外壁并用来监测聚合物的流变性质。然 而，探测器需要定位在与聚合反应器伸出的管子相联的线管外壁，由于线管包括结构细部， 这就产生了诸如聚合物堵塞线管以及维护保养复杂等困难。为了克服以上困难，我们采用 了高频频率采样，直接检测流体撞击壁面的声发射信号。通过实验发现将声波传感器直接 定位在管路的外壁，拾取流体撞击壁面的声波信号，通过分析能准确得到流体的粘度。该方 法不仅适用于管路中流体粘度的测定，同时能准确的测定搅拌釜和环管中流体的粘度。
发明内容
本发明提出了一种具有准确度高、安全环保、简易快捷、价格合理的流体粘度声发 射检测技术，通过提取声发射信号中的特征变量，并结合特征变量与流体粘度的回归模型， 进而实现流体粘度的在线检测。 —种流动体系的粘度的检测方法，包括以下步骤 1)在流动体系所在的装置外壁设置声发射信号接收装置； 2)接收流体撞击壁面所产生的声发射信号； 3)选取声发射信号中的平均能量或特征频段的能量作为特征变量； 4)将3)中选取的特征变量代入预先标定的特征变量与流体粘度间的粘度预测模型，计算得到流动体系的粘度。 所述的流动体系包括液体、气液体系或气液固体系。 所述的流动体系所在的装置指的是包括内部存在是流动体系的搅拌釜、管道和环管。 所述的声发射信号接收装置的安装在流动体系所在的装置外壁，不会影响装置内 部多相流体的运动或内部的化学反应。 所述的声发射信号中的平均能量，是通过对采集到的声发射信号，使用快速傅立 叶变进行频谱分析后，直接积分得到的。 所述的声发射信号中的特征频段的能量，是通过对采集到的声发射信号，使用快 速傅立叶变进行频谱分析，得到流体撞击壁面产生声发射信号的特征频段后，对特征频段 的进行积分得到的。 所述的选取声发射信号中特征频段判断指的是分析接收装置采集到的声发射信 号，利用快速傅立叶变换做频谱分析，将符合以下规律的区域作为特征频段当流体粘度增 加时，壁面附近的流体的速度降低，撞击壁面变得平缓，声发射信号的能量降低；当流体粘 度降低时，壁面附近的流体的速度增加，撞击壁面变得剧烈，声发射信号的能量增加。
所述的特征变量与流体粘度间的预测模型，是所述的特征变量与流体粘度成如下 的对应关系 u = A一exp (E/t》+y。 其中ii为流体的粘度；E为特征变量；A、^、y。为系数，通过预先标定的已知粘度
的流动系统的特征变量与流体粘度间的标准曲线确定。 所述的流动体系的声波信号的频率一般在0 lMHz之间。 所述的声发射传感器的接收频率范围为lkHz 2MHz，接收频率范围为lkHz 2腿z的声发射信号，接收频率范围以30kHz 500kHz为佳。
本发明与现有的方法相比具有如下一些优点 1)声发射接收装置是非插入式的，安装简易方便，不会影响装置内部多相流体的 运动或内部的化学反应； 2)不需要发射源。振动信号是流体在运动过程中自身产生的，安全环保； 3)对测量条件要求低，能在比较恶劣的环境下全天候工作，即使在高温高压等苛
刻环境下仍能正常工作； 4)反应灵敏，测量误差小，适用面广。
图1是搅拌釜冷模装置示意图； 其中，1 一搅拌釜；2_传感器；3_前置放大器；4_放大器；5_带采集卡的电脑。 图2是实施例1冷模装置中采用平均能量作为特征变量标定实验结果。 图3是实施例2的声发射信号频谱分析； 其中直线a、 b之间为特征频段，即积分区域。 图4是实施例2冷模装置中采用特征频段的能量作为特征变量标定实验结果。 图5是实施例3直径560mm搅拌釜在1. 33、 1. 50、 1. 67、 1. 83、2. 00r/s五个不同转速下采用平均能量作为特征变量的标准曲线。 图6是实施例3中实验测定的粘度数据和通过公式计算得到的粘度数据的比较。
图7是实施例4直径980mm的搅拌釜在1. 33、 1. 50、 1. 67、 1. 83 、2. OOr/s五个不同 转速下采用特征频段的能量作为特征变量的标准曲线。 图8是实施实例4中实验测定的粘度数据和通过公式计算得到的粘度数据的比 较。 图9是实施例5的管道在0. 5、1. 0、1. 5、2. 0、2. 5m/s五个不同流速下采用平均能 量作为特征变量的标准曲线。 图10是实施例6中实验测定的粘度数据和通过公式计算得到的粘度数据的比较。
图11是环管冷模装置示意图； 其中，l，-环管反应器2-传感器3-前置放大器4-放大器5-带采集卡的 电脑。 图12是实施例6的环管在0. 5、1. 0、1. 5、2. 0、2. 5m/s五个不同流速下采用平均能 量作为特征变量的标准曲线。 图13是实施例6中实验测定的粘度数据和通过公式计算得到的粘度数据的比较。
具体实施例方式
本发明将声发射信号接收装置放置于搅拌釜、管道和环管外壁面位置，通过采集 搅拌釜、管道和环管中流体撞击器壁产生的声发射信号，分析并提取声发射信号中平均能 量或特征频段的能量为特征变量E，借助预先标定建立的E与流体粘度i！之间的对应关系， 由此可以实现流体粘度的在线检测。根据检测结果，可以对搅拌釜、管道和环管进行优化设 计，达到指导生产、提高生产效率的目的。 搅拌釜、管道和环管内部的动态声发射信号通过设置在外壁处的声发射接收装置
进入放大装置进行信号的放大，以保证信号的长距离输送，然后进入声信号采集装置进行
信号的A/D转换，最后进入声发射信号处理装置(计算机)进行处理和分析。 通过置于搅拌釜、管道和环管外壁处的声发射传感器，采集声发射信号，分析计算
获得特征变量E，当流体粘度增加时，壁面附近的流体的速度降低，撞击壁面变得平缓，声
发射信号的能量降低；当搅流体粘度降低时，壁面附近的流体的速度增加，撞击壁面变得剧
烈，声发射信号的能量增加。 为实施本发明方法所设计的一套专用检测装置，包括声发射信号的接收装置，信 号采集装置以及信号处理装置。其中声发射信号的接收装置为一个或多个振动换能器；信 号采集装置为一个或多个信号采集卡(A/D转化器)；信号处理装置为带处理软件的处理 器。 所述的振动信号的接收装置的信号输出端与信号放大装置的输入端连接，信号放 大装置的输出端与信号采集装置的输入端连接，所述的信号放大装置为一个或多个信号放 大器。该放大装置可以根据实际需要选择是否使用。 振动信号接收装置的接收频率范围为lkHz 2MHz，其中接收频率范围以30kHz 500kHz为佳。
实施例1
采用如图1所示的冷模装置，搅拌釜1为高150mm、直径110mm的有机玻璃容器。 搅拌叶轮为六叶盘式涡轮，叶轮直径55mm。搅拌转速通过电子恒速搅拌器测定(杭州仪表 电机厂，额定电压220V，量程80-1000r *min—、精确到±10r *min—"。实验中采用的采样 频率为500kHz，采样时间为10s。声波测量仪器包括声波接收装置2、前置放大器3、放大器 4、信号采集装置和信号处理装置5(带采集卡的电脑)。标定实验所用物料为水-蔗糖溶 液，并通过落球法测定流体的粘度。 分析声信号接收装置采集到的声发射信号，首先对信号进行快速傅立叶变换 (FFT)，将时间域内信号转换成频域信号，得到声发射信号的频谱图，通过积分得到平均能 通过标定实验建立流体粘度y与特征变量E之间的关系，具体步骤如下将已知
粘度的流体加入搅拌釜，通过采集声发射信号计算不同粘度的平均能量作为特征变量E;
将流体粘度P做横坐标，平均能量做纵坐标，通过拟合得到指数拟合曲线图2; 通过拟合曲线得到流体粘度y与特征变量E的关系 E = 2. 57E8exp (- ii /0. 15) +1105. 56 实施例2 采用实施实例1所使用的冷模装置。标定实验所用物料为水_蔗糖溶液，并通过 落球法测定流体的粘度，同时在搅拌釜底通入气体。 分析声信号接收装置采集到的声发射信号，对信号进行快速傅立叶变换做频谱分 析，将符合以下规律的区域作为特征频段当搅拌釜中流体粘度增加时，由搅拌桨带动的流 体转速降低，撞击搅拌釜壁面变得平缓，声发射信号的特征能量降低；当搅拌釜中流体粘度 降低时，由搅拌桨带动的流体转速增加，撞击搅拌釜壁面变得剧烈，声发射信号的特征能量 增加。将流体撞击壁面产生声发射信号的特征频段(1 5kHz)积分得到能量，作为特征变 通过标定实验建立流体粘度P与特征变量之间的关系，具体步骤如下将已知粘 度的流体加入搅拌釜，通过采集声发射信号计算不同粘度的特征频段的能量E'(图3);将 流体粘度P做横坐标，特征频段的能量E'做纵坐标，通过拟合得到指数拟合曲线图4;
通过拟合曲线得到流体粘度y与特征频段的能量E'的关系
E' =3. 58E13exp (- ii /0. 10) +468278. 92
实施例3 利用声发射技术对直径560mm的搅拌釜的流体粘度进行了检测。搅拌釜高 1000mm、直径560mm的，材质为有机玻璃，转速调节范围50 250r/min。声发射信号采用声 波测量仪器进行测量。其中声波测量仪器包括声波接收装置、放大装置、信号采集装置和信 号处理装置，声传感器置于搅拌釜壁面上，采样频率lMHz，采样时间5s。标定实验所用物料 为水_蔗糖溶液，并通过落球法测定流体的粘度。通过加入不同数量的蔗糖来改变流体的 粘度，同时采集声发射信号。标定实验在五个转速(1.33、1.50、1.67、1.83、2.00r/s)下进 行，对采集到的声发射信号进行FFT变换，积分后得到平均能量作为特征变量E，与流体粘 度P作图，具体结果如图5所示。 通过拟合曲线得到转速为1. 33r/s流体粘度y与E的关系
E = 3. 97E14exp(-ii /0. 07)+1101. 17
通过拟合曲线得到转速为1. 50r/s流体粘度y与E的关系
E = 6. 82E10exp (- ii /0. 10) +1096. 55 通过拟合曲线得到转速为1. 67r/s流体粘度y与E的关系
E = 2. 42E8exp (- ii /0. 14) +1095. 09 通过拟合曲线得到转速为1. 83r/s流体粘度y与E的关系
E = 6. 52E9exp (- ii /0. 12) +1173. 27 通过拟合曲线得到转速为2. OOr/s流体粘度y与E的关系
E = 3. 58E7exp (- ii /0. 18) +1100. 58 将实验测定的粘度数据和通过公式计算得到的粘度数据比较，如图6所示。
实施例4 利用声发射技术对直径980mm的搅拌釜的流体粘度进行了检测。搅拌釜高 1600mm、直径980mm的，材质为有机玻璃，转速调节范围10 250r/min。声发射信号采用 声波测量仪器进行测量。其中声波测量仪器包括声波接收装置、放大装置、信号采集装置和 信号处理装置，声传感器置于搅拌釜壁面上，采样频率lMHz，采样时间5s。标定实验所用物 料为水-蔗糖溶液，并通过落球法测定流体的粘度。实验中向不同粘度的蔗糖溶液通入气 体，气体量维持在80m3/hr不变，同时采集声发射信号。标定实验在五个转速(1. 33、1. 50、 1. 67、 1. 83、2. 00r/s)下进行，对采集到的声发射信号进行FFT变换，积分后得到特征频段 (2 4kHz)的能量E'，与流体粘度ii作图，具体结果如图7所示。
通过拟合曲线得到转速为1. 33r/s流体粘度y与E'的关系
E' =6. 56E16exp (- ii /0. 07) +473072. 59 通过拟合曲线得到转速为1. 50r/s流体粘度y与E'的关系
E' = 3. 03E13exp(-ii/0. 10)+471431. 14
通过拟合曲线得到转速为1. 67r/s流体粘度y与E'的关系
E' = 1. 04E14exp(-ii/0. 09)+474686. 69 通过拟合曲线得到转速为1. 83r/s流体粘度y与E'的关系
E' = 4. 69Ellexp(-ii/0. 13)+488146. 24
通过拟合曲线得到转速为2. 00r/s流体粘度y与E'的关系
E' =5. 16E10exp (- ii /0. 16) +488799. 17 将实验测定的粘度数据和通过公式计算得到的粘度数据比较，如图8所示。
实施例5 利用声发射技术对的管道中的粘度进行了检测。管路直径100mm，长度5000mm。声 发射信号采用声波测量仪器进行测量。其中声波测量仪器包括声波接收装置、放大装置、信 号采集装置和信号处理装置，声传感器置于管路壁面上，采样频率lMHz，采样时间5s。标定 实验所用物料为水_蔗糖溶液，并通过落球法测定流体的粘度。通过加入不同数量的蔗糖 来改变流体的粘度，同时采集声发射信号。标定实验在五个流速(0. 5、1. 0、1. 5、2. 0、2. 5m/ s)下进行，对采集到的声发射信号进行FFT变换，积分后得到平均能量作为特征变量E，与 流体粘度P作图，具体结果如图9所示。 通过拟合曲线得到流速为0. 5m/s时流体粘度y与E的关系
E = 2. 18E8exp (- ii /0. 15) +1957. 94
通过拟合曲线得到流速为1. 0m/s时流体粘度y与E的关系
E = 6. 15E8exp (- ii /0. 15) +2240. 12 通过拟合曲线得到流速为1. 5m/s时流体粘度y与E的关系
E = 2. 95E7exp (- ii /0. 18) +2572. 79 通过拟合曲线得到流速为2. Om/s时流体粘度y与E的关系
E = 5. 19E8exp (- ii /0. 14) +2683. 10 通过拟合曲线得到流速为2. 5m/s时流体粘度y与E的关系
E = 1. 10E9exp (- ii /0. 14) +2823. 79 将实验测定的粘度数据和通过公式计算得到的粘度数据比较，如图10所示。
实施例6 采用如图ll所用的环管装置。声发射信号采用声波测量仪器进行测量。其中声波 测量仪器包括传感器、前置放大器3、放大器4、信号采集装置和信号处理装置5(带采集卡 的电脑)，声传感器置于环管反应器l的管路外壁面上，采样频率lMHz，采样时间5s。标定 实验所用物料为水_蔗糖溶液，并通过落球法测定流体的粘度。通过加入不同数量的蔗糖 来改变流体的粘度，同时采集声发射信号。标定实验在五个流速(0. 5、1. 0、1. 5、2. 0、2. 5m/ s)下进行，对采集到的声发射信号进行FFT变换，积分后得到平均能量作为特征变量E，与 流体粘度P作图，具体结果如图12所示。 通过拟合曲线得到流速为0. 5m/s时流体粘度y与E的关系
E = 1. 48E14exp (- ii /0. 07) +4180. 99 通过拟合曲线得到流速为1. Om/s时流体粘度y与E的关系
E = 7. 09E6exp (- ii /0. 22) +4967. 08 通过拟合曲线得到流速为1. 5m/s时流体粘度y与E的关系
E = 1. 19E9exp (- ii /0. 14) +5697. 25 通过拟合曲线得到流速为2. Om/s时流体粘度y与E的关系
E = 9. 85E8exp (_ ii /0. 14)+6131. 98通过拟合曲线得到流速为2. 5m/s时流体粘度y与E的关系
E = 1. 70E8exp (- ii /0. 16) +6770. 43 将实验测定的粘度数据和通过公式计算得到的粘度数据比较，如图13所示。
权利要求
一种流动体系的粘度检测方法，包括以下步骤1)在流动体系所在的装置外壁设置声发射信号接收装置；2)接收流体撞击壁面所产生的声发射信号；3)选取声发射信号中的平均能量或特征频段的能量作为特征变量；4)将3)中选取的特征变量代入预先标定的特征变量与流体粘度间的粘度预测模型，计算得到流动体系的粘度。
2. 根据权利要求l所述的粘度检测方法，其特征在于所述的流动体系包括液体、气液体系或气液固体系。
3. 根据权利要求1所述的粘度检测方法，其特征在于所述的流动体系所在的装置指内部存在流动体系的搅拌釜、管道和环管。
4. 根据权利要求1所述的粘度检测方法，其特征在于所述的平均能量通过对采集到的声发射信号，使用快速傅立叶变进行频谱分析后，直接积分得到。
5. 根据权利要求1所述的粘度检测方法，其特征在于所述的特征频段的能量通过对采集到的声发射信号，使用快速傅立叶变进行频谱分析，得到流体撞击壁面产生声发射信号的特征频段后，对特征频段的进行积分得到。
6. 根据权利要求5所述的粘度检测方法，其特征在于所述的特征频段指的是分析采集到的声发射信号，利用快速傅立叶变换做频谱分析，将符合以下规律的区域作为特征频段当流体粘度增加时，壁面附近的流体的速度降低，撞击壁面变得平缓，声发射信号的能量降低；当流体粘度降低时，壁面附近的流体的速度增加，撞击壁面变得剧烈，声发射信号的能量增加。
7. 根据权利要求1所述的粘度检测方法，其特征在于所述的特征变量与流体粘度间的预测模型，是所述的特征变量与流体粘度成如下的对应关系ii = A一exp (E/t》+y0其中P为流体的粘度，E为特征变量，Aptpy。为系数。
8. 根据权利要求l所述的粘度检测方法，其特征在于接收频率范围为lkHz 2MHz的声发射信号。
9. 根据权利要求8所述的粘度检测方法，其特征在于接收频率范围为30kHz 500kHz的声发射信号。
全文摘要
本发明公开了一种流动体系的粘度的测定方法，包括以下步骤1)在流动体系所在的装置外壁设置声发射信号接收装置；2)接收流体撞击壁面所产生的声发射信号；3)选取声发射信号中的平均能量或特征频段的能量作为特征变量；4)将3)中选取的特征变量代入预先标定的特征变量与流体粘度间的粘度预测模型，计算得到流动体系的粘度。本发明使用非侵入式的声发射装置进行检测，操作简单、安装方便，具有灵敏度高、测量误差小的特点。
文档编号G01N11/16GK101706403SQ20091015260
公开日2010年5月12日 申请日期2009年9月10日 优先权日2009年9月10日
发明者任聪静, 汪兵, 王靖岱, 阳永荣, 陈杰勋, 黄正梁 申请人:浙江大学