2. 三门峡万象实业有限公司, 河南 三门峡 472000
2. Sanmenxia Vientiane Industrual Co., Ltd., Henan Sanmenxia 472000, China
工业生产中存在设备温度偏高的问题,通常采用水冷壁降温,但因其管径较小,分布密集[1],且管路内壁常附着有溶解度低的碳酸钙垢[2],导致管道堵塞和散热困难,影响工业循环冷却水的安全运行,并对生产效率和企业的经济效益造成不利影响[3-4]。为保证企业的正常生产并降低其成本,现场经常对其进行除垢处理。目前多采用化学法除垢,该法的药剂消耗量大、成本较高[5-6],且容易造成设备腐蚀和环境污染[7-8]。与化学法相比,电磁除垢法作为一种新型的物理方法,能促进碳酸钙在水中的溶解或造成碳酸钙晶型的转变,从而抑制碳酸钙沉淀的产生[9-10],该方法绿色高效、节能环保,且适用范围广[11-12]。此外,已有磁化除垢研究普遍使用的单因素试验法,难以揭示不同因素对磁化除垢的影响规律,也造成其除垢效果不够稳定[13]。课题组的前期单因素探索试验发现,磁化时间和电流频率等因素对溶液电导率有一定影响。为了解交变磁场作用下不同因素对除垢效果的影响大小,以揭示磁化除垢的作用规律,研究采用正交试验的方法对饱和碳酸钙溶液进行磁化处理,并对试验结果进行正交多项式回归分析。在此基础上建立数学模型,通过模型进一步对磁化除垢的条件进行优化,研究可为工业生产中的除垢问题提供理论参考。
1 试验部分 1.1 试剂与仪器碳酸氢钠、无水氯化钙,分析纯,天津市河东区红岩试剂厂;一次蒸馏水;Φ0.6mm铜漆包线;1000mL烧杯。
MHS-5200A型信号发生器,郑州明禾电子科技有限公司;TDL-80-2B型低速离心机,上海安亭科学仪器厂;mp515-01型电导仪,上海三信仪表厂。
1.2 试验流程及检测方法取物质的量之比1:1的无水氯化钙和碳酸氢钠固体溶于1000mL的蒸馏水中,配制成钙离子浓度为400、600、800和1000mg/L碳酸钙溶液,分别在磁化时间为0.5~3.0h、电流频率为0.6~1.5kHz、线圈匝数为100~190圈和4种波形等不同条件下对溶液进行磁化处理,再对磁化溶液进行离心分离,并取上清液检测其电导率。试验过程中交变磁场的磁感应强度可根据赫姆霍兹定理和麦克斯韦方程计算[14]。烧杯缠绕区中点的磁感应强度与电流大小和线圈匝数成正比。
根据相关参考文献[15-16]可知,当碳酸钙沉淀溶解时,溶液中钙离子和碳酸根离子浓度增大,导电能力增强,即可用电导率反映溶液离子浓度的变化。故研究采用操作简单且精确度较高的溶液电导率来表征溶液钙离子浓度的变化以反映除垢效果。检测时取100mL的磁化溶液,在3000r/min的转速下离心处理30min,最后使用电导仪检测上清液的电导率,并以5次检测的平均值作为最终结果。
2 正交试验设计参考已有研究[17]并结合理论分析,选取磁化时间(A)、电流频率(B)、线圈匝数(C)、c0(Ca2+)(D)及波形(E)等5个主要因素进行试验。设计试验时采用5因素4水平L16(45)的正交方法,并将5个因素的水平选在适当范围内,正交因素及其对应水平如表 1所示,正交试验的结果如表 2所示。
水平 | A 磁化 时间/h |
B 电流 频率/kHz |
C 线圈 匝数/圈 |
D c0(Ca2+)/ (mg·L-1) |
E 波形 |
1 | 0.5 | 0.6 | 100 | 400 | 三角波 |
2 | 1.0 | 0.9 | 130 | 600 | 正弦波 |
3 | 2.0 | 1.2 | 160 | 800 | 方波 |
4 | 3.0 | 1.5 | 190 | 1000 | 升锯齿波 |
试验 编号 |
A | B | C | D | E | 电导率/ (mS·cm-1) |
1 | 1 | 1 | 1 | 1 | 1 | 0.776 |
2 | 1 | 2 | 2 | 2 | 2 | 0.768 |
3 | 1 | 3 | 3 | 3 | 3 | 1.216 |
4 | 1 | 4 | 4 | 4 | 4 | 1.413 |
5 | 2 | 1 | 2 | 3 | 4 | 1.267 |
6 | 2 | 2 | 1 | 4 | 3 | 1.375 |
7 | 2 | 3 | 4 | 1 | 2 | 0.881 |
8 | 2 | 4 | 3 | 2 | 1 | 0.797 |
9 | 3 | 1 | 3 | 4 | 2 | 1.436 |
10 | 3 | 2 | 4 | 3 | 1 | 1.258 |
11 | 3 | 3 | 1 | 2 | 4 | 0.702 |
12 | 3 | 4 | 2 | 1 | 3 | 0.691 |
13 | 4 | 1 | 4 | 2 | 3 | 0.749 |
14 | 4 | 2 | 3 | 1 | 4 | 0.631 |
15 | 4 | 3 | 2 | 4 | 1 | 1.341 |
16 | 4 | 4 | 1 | 3 | 2 | 1.084 |
EⅠ | 4.173 | 4.228 | 3.937 | 2.979 | 4.172 | ET=16.385 |
EⅡ | 4.320 | 4.032 | 4.067 | 3.016 | 4.169 | |
EⅢ | 4.087 | 4.140 | 4.080 | 4.825 | 4.031 | |
EⅣ | 3.805 | 3.985 | 4.301 | 5.565 | 4.013 | |
S | 0.035 | 0.009 | 0.017 | 1.276 | 0.006 | CT=16.779 |
注:EⅠ、EⅡ、EⅢ、EⅣ为各因子在同一水平所对应的数据之和;ET为电导率的数据总和;CT为校正项,CT=E2T/16;S为变差平方和。 |
由表 2可知,试验编号为9条件下的磁化溶液的电导率最高,达到了1.436mS/cm,说明此时溶液中钙离子的浓度最高,除垢效果最明显。试验编号为14的磁化溶液的电导率最低为0.631mS/cm。试验编号为3、4、5、6、10、15和16的磁化溶液的电导率相对较高,在1.084~1.413mS/cm之间,试验编号为1、2、7、8、11、12和13的磁化溶液的电导率相对较低,在0.691~0.881mS/cm之间。
结合他人研究可知,在交变磁场的作用下,钙离子的电化学特性和物理特性如黏附特性会发生一定改变, 并使溶液中的碳酸钙晶核数目减少且难以长大[19]。此外,碳酸钙晶体的稳定性受离子的极化程度影响,离子的极化程度愈高, 碳酸钙晶体的稳定性愈低[20]。交变磁场可促使碳酸钙晶体发生极化,造成碳酸钙晶体的稳定性降低并溶解于溶液,所以溶液中的离子数量增多且电导率增大,也即通过交变磁场的适当磁化处理可以促进碳酸钙晶体的溶解因而达到除垢的效果。
3 正交试验结果的回归分析根据相关研究[18],可利用正交多项式处理正交试验结果,进而得到溶液电导率与各因素之间的定量关系。方差分析的各项结果见表 3及表 4。
编号 | φ1(A) | φ2(A) | φ3(A) | φ1(B) | φ2(B) | φ3(B) | φ1(C) | φ2(C) | φ3(C) | φ1(D) | φ2(D) | φ3(D) | φ1(E) | φ2(E) | φ3(E) |
1 | -3 | 1 | -1 | -3 | 1 | -1 | -3 | 1 | -1 | -3 | 1 | -1 | -3 | 1 | -1 |
2 | -3 | 1 | -1 | -1 | -1 | 3 | -1 | -1 | 3 | -1 | -1 | 3 | -1 | -1 | 3 |
3 | -3 | 1 | -1 | 1 | -1 | -3 | 1 | -1 | -3 | 1 | -1 | -3 | 1 | -1 | -3 |
4 | -3 | 1 | -1 | 3 | 1 | 1 | 3 | 1 | 1 | 3 | 1 | 1 | 3 | 1 | 1 |
5 | -1 | -1 | 3 | -3 | 1 | -1 | -1 | -1 | 3 | 1 | -1 | -3 | 3 | 1 | 1 |
6 | -1 | -1 | 3 | -1 | -1 | 3 | -3 | 1 | -1 | 3 | 1 | 1 | 1 | -1 | -3 |
7 | -1 | -1 | 3 | 1 | -1 | -3 | 3 | 1 | 1 | -3 | 1 | -1 | -1 | -1 | 3 |
8 | -1 | -1 | 3 | 3 | 1 | 1 | 1 | -1 | -3 | -1 | -1 | 3 | -3 | 1 | -1 |
9 | 1 | -1 | -3 | -3 | 1 | -1 | 1 | -1 | -3 | 3 | 1 | 1 | -1 | -1 | 3 |
10 | 1 | -1 | -3 | -1 | -1 | 3 | 3 | 1 | 1 | 1 | -1 | -3 | -3 | 1 | -1 |
11 | 1 | -1 | -3 | 1 | -1 | -3 | -3 | 1 | -1 | -1 | -1 | 3 | 3 | 1 | 1 |
12 | 1 | -1 | -3 | 3 | 1 | 1 | -1 | -1 | 3 | -3 | 1 | -1 | 1 | -1 | -3 |
13 | 3 | 1 | 1 | -3 | 1 | -1 | 3 | 1 | 1 | -1 | -1 | 3 | 1 | -1 | -3 |
14 | 3 | 1 | 1 | -1 | -1 | 3 | 1 | -1 | -3 | -3 | 1 | -1 | 3 | 1 | 1 |
15 | 3 | 1 | 1 | 1 | -1 | -3 | -1 | -1 | 3 | 3 | 1 | 1 | -3 | 1 | -1 |
16 | 3 | 1 | 1 | 3 | 1 | 1 | -3 | 1 | -1 | 1 | -1 | -3 | -1 | -1 | 3 |
Bi | -1.34 | -0.43 | 0.33 | -0.62 | 0.04 | -0.57 | 1.11 | 0.09 | 0.33 | 9.57 | 0.70 | -2.84 | -0.62 | -0.02 | 0.03 |
Di | 160 | 32 | 160 | 160 | 32 | 160 | 160 | 32 | 160 | 160 | 32 | 160 | 160 | 32 | 160 |
di×10-3 | -8.38 | -13.44 | 2.06 | -3.88 | 1.25 | -3.56 | 6.94 | 2.81 | 2.06 | 59.81 | 21.88 | -17.75 | -3.88 | -0.63 | 0.19 |
λi | 2 | 1 | 3/10 | 2 | 1 | 3/10 | 2 | 1 | 3/10 | 2 | 1 | 3/10 | 2 | 1 | 3/10 |
Q×10-3 | 11.22 | 5.78 | 0.68 | 2.40 | 0.05 | 2.03 | 7.70 | 0.25 | 0.68 | 572.41 | 15.31 | 50.41 | 2.40 | 0.01 | 0.01 |
效应项 | A1 | A2 | A3 | B1 | B2 | B3 | C1 | C2 | C3 | D1 | D2 | D3 | E1 | E2 | E3 |
注:Bi=∑(cη);Di=snr;di=Bi/Di;c为多项式系数;η为电导率;n为试验重复次数,r为重复使用S的次数;λi为调整系数;Q为各个因子的平方和。 |
方差来源 | 离差×10-3 | f | 均方离差×10-3 | F | 显著性 | F表值 |
A/h | 17.68 | 3 | 5.89 | 7.27 | * | F0.10(3, 3)=5.39 |
A1 | 11.22 | 1 | 11.22 | 13.85 | (*) | F0.10(1, 3)=5.54 |
A2 | 5.78 | 1 | 5.78 | 7.14 | F0.01(1, 3)=34.12 | |
A3 | 0.68 | 1 | 0.68 | 0.84 | F0.01(1, 3)=34.12 | |
B/kHz | 4.48 | 3 | 1.49 | 1.84 | F0.10(3, 3)=5.39 | |
B1 | 2.40 | 1 | 2.40 | 2.96 | F0.10(1, 3)=5.54 | |
B2 | 0.05 | 1 | 0.05 | 0.06 | F0.01(1, 3)=34.12 | |
B3 | 2.03 | 1 | 2.03 | 2.51 | F0.01(1, 3)=34.12 | |
C/N | 8.63 | 3 | 2.87 | 3.54 | F0.10(3, 3)=5.39 | |
C1 | 7.70 | 1 | 7.70 | 9.51 | F0.10(1, 3)=5.54 | |
C2 | 0.25 | 1 | 0.25 | 0.31 | F0.01(1, 3)=34.12 | |
C3 | 0.68 | 1 | 0.68 | 0.84 | F0.01(1, 3)=34.12 | |
D/(mg·L-1) | 638.13 | 3 | 212.71 | 262.60 | ** | F0.10(3, 3)=5.39 |
D1 | 572.41 | 1 | 572.41 | 706.68 | (**) | F0.10(1, 3)=5.54 |
D2 | 15.31 | 1 | 15.31 | 18.90 | F0.01(1, 3)=34.12 | |
D3 | 50.41 | 1 | 50.41 | 62.23 | (*) | F0.01(1, 3)=34.12 |
E | 2.42 | 3 | 0.81 | |||
E1 | 2.40 | |||||
E2 | 0.01 | |||||
E3 | 0.01 | |||||
误差 | 2.42 | |||||
总和 | 671.34 |
当试验水平为4时,函数关系为三次多项式,设三次多项式的回归方程为:
$ \hat y = {b_0} + {b_1}{\phi _1}\left( x \right) + b{\phi _2}\left( x \right) + {b_3}{\varphi _3}\left( x \right) $ | (1) |
式(1)中y为响应值,b为回归系数,ϕ1(x)、ϕ2(x)和ϕ3(x)分别为x的一次、二次和三次多项式。
由表 4可知,因子A1、D1和D3为显著影响因素,其余因子为不显著影响因素。含有影响因子A1、D1和D3的回归方程式为:
$ \begin{array}{l} \hat y = \bar y + {b_{A1}}{\lambda _1}\left( {A' - \bar A'} \right) + {b_{D3}}{\lambda _3}\left( {D' - \bar D'} \right) + \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;{b_{D1}}{\lambda _1}\left[ {{{\left( {D' - \bar D'} \right)}^2} - \frac{{{N^2} - 1}}{{12}}} \right] \end{array} $ | (2) |
相关系数代入式(2)可得出方程:
$ \begin{array}{l} \hat y = - 33.424 \times {10^{ - 3}}{X_A} - 4.213 \times \\ {10^{ - 3}}{X_D} + 2.990 \times {10^{ - 6}}X_D^2 + 2.319 \end{array} $ | (3) |
式(3)为表述磁化溶液的电导率与A磁化时间和D溶液初始钙离子浓度间相关关系的线性回归方程,即磁化除垢的数学模型。根据该方程可计算溶液电导率残差的标准差,结果见表 5。
A/ h |
D/ (mg·L-1) |
ŷ | |||
0.5 | 400 | 0.776 | 1.089 | -0.319 | 0.102 |
0.5 | 600 | 0.768 | 0.851 | -0.083 | 0.007 |
0.5 | 800 | 1.216 | 0.845 | 0.371 | 0.137 |
0.5 | 1000 | 1.413 | 1.079 | 0.334 | 0.111 |
1.0 | 800 | 1.267 | 0.829 | 0.438 | 0.192 |
1.0 | 1000 | 1.375 | 1.063 | 0.312 | 0.098 |
1.0 | 400 | 0.881 | 1.079 | -0.198 | 0.039 |
1.0 | 600 | 0.797 | 0.834 | -0.037 | 0.001 |
2.0 | 1000 | 1.436 | 1.029 | 0.407 | 0.166 |
2.0 | 600 | 1.258 | 0.801 | 0.457 | 0.209 |
2.0 | 800 | 0.702 | 0.795 | -0.093 | 0.009 |
2.0 | 400 | 0.691 | 1.045 | -0.354 | 0.126 |
3.0 | 600 | 0.749 | 0.767 | -0.018 | 0.000 |
3.0 | 400 | 0.631 | 1.012 | -0.381 | 0.145 |
3.0 | 1000 | 1.341 | 0.996 | 0.345 | 0.119 |
3.0 | 800 | 1.084 | 0.762 | 0.322 | 0.104 |
注: |
由表 5可计算出实验值与预测值的误差,计算过程如下:
$ \sum {{{\left( {\bar y - \hat y} \right)}^2}} = 1.565 $ | (4) |
$ S = \sqrt {\left[ {\frac{{\sum {{{\left( {\bar y' - \hat y} \right)}^2}} }}{{m - 2}}} \right]} = \sqrt {\frac{{1.565}}{{16 - 2}}} = 0.3343 $ | (5) |
由式(5)可知,m为测试次数,S为试验值与预测值的误差。此处磁化碳酸钙溶液电导率的误差为0.3343。通过残差估计可以在一定范围内计算出的不同条件下溶液的电导率,且该计算值与试验值的误差在(±0.3343)%之间。由此可见,通过回归分析建立的数学模型的准确性较高,可在工业生产中普遍应用。
4 结论1) 磁化除垢的优化条件为:磁化时间为0.5h,电流频率为0.6kHz,线圈匝数为100圈,c0(Ca2+)为1000mg/L,波形为方波。在此条件下获得磁化溶液的电导率的平均值为1.474mS/cm,比相同条件下的未磁化溶液增加了0.147mS/cm,此条件下磁化除垢的效果最好。
2) 磁化时间和溶液初始钙离子浓度是影响溶液电导率的显著因素,而电流频率、线圈匝数和波形为不显著因素。研究建立了磁化时间和初始钙离子浓度与溶液电导率之间相关关系的数学模型,该模型的准确性和可靠性较高,可为揭示磁化除垢的作用规律提供理论参考,也可为工业生产中除垢问题的解决提供新的思路。
[1] |
刘秀菊, 王玲. 冶金工业循环冷却水中微生物的危害及控制[J]. 黑龙江冶金, 2005, 25(3): 47-48. Liu Xiuju, Wang Ling. Harm and control of microorganism in circulating cooling water in metallurgical industry[J]. Metallurgy and Materials, 2005, 25(3): 47-48. (in Chinese) |
[2] |
张玉玲, 李伟, 李旭东, 等. 天冬氨酸-柠檬酸共聚物对碳酸钙的抑制性能[J]. 精细化工, 2018, 35(11): 1877-1884. Zhang Yuling, Li Wei, Li Xudong, et al. Scale inhibition performance of copoly(aspartic acid-citric acid) on calcium carbonate[J]. Fine Chemicals, 2018, 35(11): 1877-1884. (in Chinese) |
[3] |
蒋文斌, 陈华德, 李青, 等. 用于除垢的新型电磁场发生装置设计[J]. 中国计量学院学报, 2011, 22(3): 263-267. Jiang Wenbin, Chen Huade, Li Qing, et al. Design of a new electromagnetic field generator for anti-scaling[J]. Journal of China University of Metrology, 2011, 22(3): 263-267. (in Chinese) |
[4] |
韩良敏, 冯毅. 磁场抑垢的影响因素与机理研究[J]. 给水排水, 2010, 46(S1): 378-380. Han Liangmin, Feng Yi. Influencing factors and mechanism research of magnetic anti-scaling[J]. Water & Wastewater Engineering, 2010, 46(S1): 378-380. (in Chinese) |
[5] |
陈璟. 机械除垢与化学除垢技术在现场的应用[J]. 吐哈油气, 2010, 15(4): 386-388. Chen Jing. Field applications of mechanical cleaning and chemical cleaning technique[J]. Tuha Oil & Gas, 2010, 15(4): 386-388. (in Chinese) |
[6] |
王兵, 李长俊, 朱伟, 等. 结垢及除垢技术在管道中的应用研究[J]. 石油化工腐蚀与防护, 2008, 25(1): 28-30, 36. Wang Bing, Li Changjun, Zhu Wei, et al. Appfication of de-fouling technology in pipelines[J]. Corrosion & Protection in Petrochemical Industry, 2008, 25(1): 28-30, 36. (in Chinese) |
[7] |
肖曾利, 蒲春生. 磁防垢技术及应用发展现状[J]. 断块油气田, 2010, 17(1): 121-125. Xiao Zengli, Pu Chunsheng. Review of technology and application on magnetized water treatment for anti-scaling[J]. Fault-Block Oil & Gas Field, 2010, 17(1): 121-125. (in Chinese) |
[8] |
陆希希, 郭越, 林雪松, 等. 超声波除碳酸钙垢的实验研究[J]. 当代化工, 2017, 46(7): 1318-1321. Lu Xixi, Guo Yue, Lin Xuesong, et al. Study on calcium carbonate scale removal with ultrasonic wave in laboratory[J]. Contemporary Chemical Industry, 2017, 46(7): 1318-1321. (in Chinese) |
[9] |
费继友, 李玉泉, 白鑫. 水的电磁变频除垢防垢技术和实验研究[J]. 化工自动化及仪表, 2011, 38(2): 157-161. Fei Jiyou, Li Yuquan, Bai Xin. Research on electromagnetic antiscaling and descaling in water[J]. Control and Instruments in Chemical Industry, 2011, 38(2): 157-161. (in Chinese) |
[10] |
孔垂鑫.交变电磁场抑垢机理实验研究[D].济南: 山东大学, 2016 Kong Chuixin. Experimental study on the anti-scaling mechanism of the alternating electromagnetic field[D]. Jinan: Shandong University, 2016(in Chinese) http://cdmd.cnki.com.cn/Article/CDMD-10422-1016164434.htm |
[11] |
Cho Y I, Lane J, Kim W. Pulsed-Power treatment for physical water treatment[J]. International Communications in Heat and Mass Transfer, 2005, 32(7): 861-871. |
[12] |
Yacoby M J, Yacoby M J. The effect of magnetic and constant electric field antiscaling devices on the mechanism of CaCO3 scale formation[J]1995, doi:http://hdl.handle.net/10210/9091
|
[13] |
赵丹亚, 郭平. 磁化处理水防垢除垢机理探讨[J]. 云南工学院学报, 1994(1): 83-88. Zhao Danya, Guo Ping. Approaching of the machanism of scale removing and preventing of magnetized water[J]. Journal of Yunnan Institute of Technology, 1994(1): 83-88. (in Chinese) |
[14] |
姜德宁, 周昆. 基于交变磁场的水管防垢试验研究[J]. 天津工程师范学院学报, 2008, 18(1): 21-25. Jiang Dening, Zhou Kun. Experimental study on antiscale of pipes based on alternating magnetic field[J]. Journal of Tianjing University of Technology and Education, 2008, 18(1): 21-25. (in Chinese) |
[15] |
杨海涛. 电导率法探究电磁场强度对碳酸钙结晶的影响[J]. 电子制作, 2014(5): 30-30. Yang Haitao. Exploration of the effect of electromagnetic field intensity on crystallization of calcium carbonate by conductivity method[J]. Practical Electronics, 2014(5): 30-30. (in Chinese) |
[16] |
张小霓, 于萍, 罗运柏. 溶液电导率法对碳酸钙结晶动力学的研究[J]. 应用化学, 2004, 21(2): 187-191. Zhang Xiaoni, Yu Ping, Luo Yunbai. Study on kinetics of formation of CaCO3 crystal by solution conductivity measurement[J]. Chinese Journal of Applied Chemistry, 2004, 21(2): 187-191. (in Chinese) |
[17] |
王宇斌, 彭祥玉, 张小波, 等. 基于正交实验的红土镍矿酸浸条件优化[J]. 无机盐工业, 2017, 49(6): 29-32. Wang Yubin, Peng Xiangyu, Zhang Xiaobo, et al. Optimization of nickel laterite acid leaching conditions based on orthogonal experiment[J]. Inorganic Chemicals Industry, 2017, 49(6): 29-32. (in Chinese) |
[18] |
王宇斌, 文堪, 彭祥玉, 等. 镁质红土镍矿酸浸结果影响因素的数学模型[J]. 中国有色冶金, 2018, 47(3): 76-81. Wang Yubin, Wen Kan, Peng Xiangyu, et al. Mathematical model of factors affecting the acid leaching results of magnesian nickel laterite[J]. China Nonferrous Metallurgy, 2018, 47(3): 76-81. (in Chinese) |
[19] |
顾晓.新型电子防垢技术在大港油田注水系统的应用[C]//中国石油学会全国油田注水技术研讨会, 2005 http://www.wanfangdata.com.cn/details/detail.do?_type=conference&id=6129172
|
[20] |
刘军营, 张明, 李素玲. 水垢的预测与磁化处理技术[J]. 节能技术, 2002, 20(5): 27-29. Liu Junying, Zhang Ming, Li Suling. Forecasting and magnetic disposing technology for incrustation[J]. Energy Conservation Technology, 2002, 20(5): 27-29. (in Chinese) |