磁化除垢影响因素作用规律及其数学模型研究

王宇斌，王望泊*，戚晓宇，王 刚，王 妍

(1.西安建筑科技大学 资源工程学院，西安 710055； 2.三门峡万象实业有限公司，河南 三门峡 472000)

工业生产中存在设备温度偏高的问题，通常采用水冷壁降温，但因其管径较小，分布密集[1]，且管路内壁常附着有溶解度低的碳酸钙垢[2]，导致管道堵塞和散热困难，影响工业循环冷却水的安全运行，并对生产效率和企业的经济效益造成不利影响[3-4]。为保证企业的正常生产并降低其成本，现场经常对其进行除垢处理。目前多采用化学法除垢，该法的药剂消耗量大、成本较高[5-6]，且容易造成设备腐蚀和环境污染[7-8]。与化学法相比，电磁除垢法作为一种新型的物理方法，能促进碳酸钙在水中的溶解或造成碳酸钙晶型的转变，从而抑制碳酸钙沉淀的产生[9-10]，该方法绿色高效、节能环保，且适用范围广[11-12]。此外，已有磁化除垢研究普遍使用的单因素试验法，难以揭示不同因素对磁化除垢的影响规律，也造成其除垢效果不够稳定[13]。课题组的前期单因素探索试验发现，磁化时间和电流频率等因素对溶液电导率有一定影响。为了解交变磁场作用下不同因素对除垢效果的影响大小，以揭示磁化除垢的作用规律，研究采用正交试验的方法对饱和碳酸钙溶液进行磁化处理，并对试验结果进行正交多项式回归分析。在此基础上建立数学模型，通过模型进一步对磁化除垢的条件进行优化，研究可为工业生产中的除垢问题提供理论参考。

1 试验部分

1.1 试剂与仪器

碳酸氢钠、无水氯化钙，分析纯，天津市河东区红岩试剂厂；一次蒸馏水；Φ0.6 mm铜漆包线；1 000 mL烧杯。

MHS-5200A型信号发生器，郑州明禾电子科技有限公司；TDL-80-2B型低速离心机，上海安亭科学仪器厂；mp515-01型电导仪，上海三信仪表厂。

1.2 试验流程及检测方法

取物质的量之比1∶1的无水氯化钙和碳酸氢钠固体溶于1 000 mL的蒸馏水中，配制成钙离子浓度为400、600、800和1 000 mg/L碳酸钙溶液，分别在磁化时间为0.5～3.0 h、电流频率为0.6～1.5 kHz、线圈匝数为100～190圈和4种波形等不同条件下对溶液进行磁化处理，再对磁化溶液进行离心分离，并取上清液检测其电导率。试验过程中交变磁场的磁感应强度可根据赫姆霍兹定理和麦克斯韦方程计算[14]。烧杯缠绕区中点的磁感应强度与电流大小和线圈匝数成正比。

根据相关参考文献[15-16]可知，当碳酸钙沉淀溶解时，溶液中钙离子和碳酸根离子浓度增大，导电能力增强，即可用电导率反映溶液离子浓度的变化。故研究采用操作简单且精确度较高的溶液电导率来表征溶液钙离子浓度的变化以反映除垢效果。检测时取100 mL的磁化溶液，在3 000 r/min的转速下离心处理30 min，最后使用电导仪检测上清液的电导率，并以5次检测的平均值作为最终结果。

2 正交试验设计

参考已有研究[17]并结合理论分析，选取磁化时间(A)、电流频率(B)、线圈匝数(C)、c0(Ca2+) (D)及波形(E)等5个主要因素进行试验。设计试验时采用5因素4水平L16(45)的正交方法，并将5个因素的水平选在适当范围内，正交因素及其对应水平如表1所示，正交试验的结果如表2所示。

表1 正交试验各因素水平用量Table 1 Different levels of factor and the amount of orthogonal experiment

由表2可知，试验编号为9条件下的磁化溶液的电导率最高，达到了1.436 mS/cm，说明此时溶液中钙离子的浓度最高，除垢效果最明显。试验编号为14的磁化溶液的电导率最低为0.631 mS/cm。试验编号为3、4、5、6、10、15和16的磁化溶液的电导率相对较高，在1.084～1.413 mS/cm之间，试验编号为1、2、7、8、11、12和13的磁化溶液的电导率相对较低，在0.691～0.881 mS/cm之间。

结合他人研究可知，在交变磁场的作用下，钙离子的电化学特性和物理特性如黏附特性会发生一定改变, 并使溶液中的碳酸钙晶核数目减少且难以长大[19]。此外，碳酸钙晶体的稳定性受离子的极化程度影响，离子的极化程度愈高, 碳酸钙晶体的稳定性愈低[20]。交变磁场可促使碳酸钙晶体发生极化，造成碳酸钙晶体的稳定性降低并溶解于溶液，所以溶液中的离子数量增多且电导率增大，也即通过交变磁场的适当磁化处理可以促进碳酸钙晶体的溶解因而达到除垢的效果。

表2 磁化除垢正交试验结果Table 2 Result table of orthogonal experiment of magnetic disposing incrustation

3 正交试验结果的回归分析

根据相关研究[18]，可利用正交多项式处理正交试验结果，进而得到溶液电导率与各因素之间的定量关系。方差分析的各项结果见表3及表4。

当试验水平为4时，函数关系为三次多项式，设三次多项式的回归方程为：

(1)

式(1)中y为响应值，b为回归系数，φ1(x)、φ2(x)和φ3(x)分别为x的一次、二次和三次多项式。

由表4可知，因子A1、D1和D3为显著影响因素，其余因子为不显著影响因素。含有影响因子A1、D1和D3的回归方程式为：

(2)

相关系数代入式(2)可得出方程：

(3)

表3 方差的各项结果Table 3 Variance results

注：Bi=∑(cη)；Di=snr；di=Bi/Di；c为多项式系数；η为电导率；n为试验重复次数，r为重复使用S的次数；λi为调整系数；Q为各个因子的平方和。

表4 磁化溶液电导率的各个水平方差分析结果 Table 4 Variance analysis of various factors on the impact of the conductivity of the magnetized solution

式(3)为表述磁化溶液的电导率与A磁化时间和D溶液初始钙离子浓度间相关关系的线性回归方程，即磁化除垢的数学模型。根据该方程可计算溶液电导率残差的标准差，结果见表5。

表5 溶液电导率的残差估计Table 5 Residual estimate of the conductivity of the magnetized solution

由表5可计算出实验值与预测值的误差，计算过程如下：

(4)

(5)

由式(5)可知，m为测试次数，S为试验值与预测值的误差。此处磁化碳酸钙溶液电导率的误差为0.334 3。通过残差估计可以在一定范围内计算出的不同条件下溶液的电导率，且该计算值与试验值的误差在(±0.334 3)%之间。由此可见，通过回归分析建立的数学模型的准确性较高，可在工业生产中普遍应用。

4 结论

1)磁化除垢的优化条件为：磁化时间为0.5 h，电流频率为0.6 kHz，线圈匝数为100圈，c0(Ca2+)为1 000 mg/L，波形为方波。在此条件下获得磁化溶液的电导率的平均值为1.474 mS/cm，比相同条件下的未磁化溶液增加了0.147 mS/cm，此条件下磁化除垢的效果最好。

2)磁化时间和溶液初始钙离子浓度是影响溶液电导率的显著因素，而电流频率、线圈匝数和波形为不显著因素。研究建立了磁化时间和初始钙离子浓度与溶液电导率之间相关关系的数学模型，该模型的准确性和可靠性较高，可为揭示磁化除垢的作用规律提供理论参考，也可为工业生产中除垢问题的解决提供新的思路。