中国有色金属学报 2004,(04),681-685 DOI:10.19476/j.ysxb.1004.0609.2004.04.028
预焙铝电解槽电流效率与阳极电流分布的数学模型
张秋萍
北方工业大学现场总线与自动化重点实验室,北方工业大学现场总线与自动化重点实验室 北京100041 ,北京100041
摘 要:
在复杂的综合数学模型的基础上, 利用正交多元回归法研究了预焙电解槽中电流效率与阳极电流分布的关系, 得到一个代数方程式, 同时, 用这个代数方程式分析了阳极电流分布与电流效率的关系。分析结果表明:阳极电流分布与电流效率都随时间和空间而改变;阳极电流分布的改变引起电流效率的改变;在KuhnTucker理论的基础上, 还讨论了系列电流不变时的最佳电流效率, 由于电解槽中磁场、流场分布不均匀等原因, 并非严格均匀的阳极电流分布才能得到最高的电流效率。
关键词:
铝电解 ;正交多元回归 ;电流效率 ;电流分布 ;
中图分类号: TF351
收稿日期: 2003-08-07
基金: 国家自然科学基金资助项目 (59974001); 北京市现场总线与自动化重点实验室资助项目;
Mathematical model of current efficiency with anode current distribution in aluminum production
Abstract:
On the bases of the comprehensive mathematical model, the relation between current efficiency and anode current distribution in aluminum production was examined through using orthoganal regression design method. An algebraic equation was developed under some reasonable simplification, and based on it, the variation of current efficiency and anode current distribution was analyzed. The results show that both current efficiency and current distribution change with time and position. By use of the Kuhn Tucker theory, the optional current efficiency under given total serial current was discussed. Not strictly smooth anode current distribution results in optional current efficiency because of the uneven physical fields distribution.
Keyword:
current efficiency; aluminum electrolysis; current distribution; orthogonal regression design;
Received: 2003-08-07
铝电解槽中电流分布瞬息万变, 在不同的时间和空间有不同的值, 是电解过程中最难掌握的参数。 电流效率是铝电解过程最重要的技术经济指标, 在铝生产中有很多因素影响电流效率, 例如: 温度, 电流强度, 电流密度, 金属流动, 电解质组成, 极间距离, 槽龄等。 通过分析可以发现几乎所有与电流效率有关的因素都与电流分布有关, 或者直接或间接地影响电流分布
[1 ,2 ]
, 金属流动与电流分布有关, 温度及其分布通过改变炉膛形状来改变电流分布, 电解质成分也可以由电流分布的变化来反映, 极间的距离和铝水平也影响电流分布。 由于在任何情况下, 外部电流分布及其变化比其他任何参数更容易测定。 因此, 研究电流分布与电流效率之间的关系非常有意义, 也是电流效率连续监测的基础。
国内外这一方面的工作报道较少, 但研究电流分布和研究电流效率的报道很多, 这些研究对本文研究具有重要的参考意义
[3 ,4 ,5 ,6 ]
。
预焙电解槽物理场的特征在许多文章中有描述, 针对某厂160 kA预焙电解槽, 我们也做过一系列的工作
[7 ,8 ,9 ,10 ]
。 结论如下: 现行的母线配置下垂直磁场较大且严重不均匀, 某一电解槽区域垂直磁场达到4×10-3 T, 金属流动速度较快且不均匀, 平均流动速度为8 cm/s, 最大流动速度大于20cm/s。 电解槽内的电流分布受炉帮和单个阳极电流的影响, 水平电流随炉帮的变化而变化, 如果炉帮延伸到阳极底掌投影下, 则水平电流增加, 引起铝液明显上下波动, 这就导致电流效率很低。
1电流效率的计算模型
1.1 电流效率的综合数学模型
有关电流效率的机理模型和经验模型有许多报道, 但都不能满足分析电流分布与电流效率关系的需要。 在Lillebuen模型基础上
[11 ]
, 充分考虑整个传质过程中的所有因素, 通过重新定义及重新计算模型中的一些变量, 推导了一个铝电解过程中电流效率的综合数学模型。 下面仅给出方程式, 模型的推导和参数的物理意义见文献
[
12 ]
。
{ η = 1 0 0 - 2 1 9 A Ι - 1 D 0 . 6 7 m e μ - 0 . 5 ? ? ? U - 0 . 8 3 e d - 0 . 1 7 ρ 1 . 5 c % w ( A l ) ( 1 - f ) f = c m / c * m D m e = D m ( σ r / σ ) 0 . 5 ? ? ? ( 1 )
进一步研究模型中的变量, 可知它们由下述工作参数决定, 即温度 (T ) , 电流密度 (J c ) , 熔体流动速度 (u ) , 阴阳极间距离 (d ) , 界面张力 (σ ) , 金属饱和度分数 (f ) 。 已知这些参数就可以通过方程 (1) 计算出电流效率, 如果电解槽炉帮形状变化不是很快的话, 这些参数可以由电流分布直接计算得到或者由半经验公式粗略决定。 熔体流动速度 (u ) 主要取决于电磁力和重力, 随着时间变化的电磁力只取决于磁场和电流分布, 因此, 通过解关于热场、 电场、 磁场和流场的方程组, 即可求得熔体流动速度。 电流密度d c 由直接解电场方程组决定, T , d , σ 和f 等参数由半经验公式 (2) 计算得到
[13 ]
:
{ Τ = Τ 0 + B J e ( J e - G ) J i = C i / Ι i σ = C 0 - C 2 J c f = f 0 - C 4 J a ? ? ? ( 2 )
这里T 0 , B , G , C i , C 0 , C 2 , f 0 , C 4 是实验结果计算得到的常量。 J e , J c , J a 分别是电解质电流密度, 阴极电流密度和阳极电流密度, 即电流分布。 所以方程式 (1) 、 (2) 给出了电流效率和电流分布的关系, 如果已知电流分布, 就可以计算出电流效率。 由于在大型电解槽中电流分布不均匀, 因此电流效率也不一致。 采用区域参数对应的是区域电流效率, 采用平均参数对应的就是平均电流效率。 本文作者采用二氧化碳分析法和铜稀释法, 仅仅用整个电解槽的平均电流效率验证模型。 实验结果表明, 一周时间内电流效率的计算值和测量值相对误差为0.5%。 因此, 方程式 (1) 、 (2) 可以分析电流分布与电流效率的关系, 但从方程组知道电流分布与电流效率是通过很复杂的物理场计算模型联系起来, 计算工作量大而且通过表达式很难分析它们之间的关系。 为了解决这个问题, 我们采用正交多元回归法简化模型。
1.2 电流效率的代数方程模型
以160 kA预焙电解槽为例, 此方法可以推广到其他电解槽。 将160 kA预焙电解槽划分成9个区域
[11 ]
, 其中8个区域有垂直电流, 而另一个没有垂直电流, 除第9区域外, 其余每个区域由3个阳极组成, 区域划分见图1。 当然区域可以任意划分, 但若在同一区域内有一致的物理场, 分析会更方便。 如果区域划分更细, 结果更精确, 但计算却更复杂。
图1 160 kA预焙电解槽的分区
Fig.1 Zone distribution of160 kA prebake cell
从区域的划分得知, 共有8个区有电流进入, 方程式中应包含8个变量。 可用8元正交回归方法解决这个问题, 建立代数方程直接显示每区电流对平均电流效率的影响。 正交多元回归法是设计实验、 建立模型的有效方法。 首先选点计算基于复杂模型的一系列电流效率值, 然后用这些值得到一个代数方程式。
建立回归模型所需的电流效率值的数目为
[14 ]
:
表1 正交回归的表头设计
Table 1 Head design for orthogonal regression table
x 0
x 1
x 2
x 3
x 4
x 5
x 6
x 7
x 8
x 2 1
x 2 2
x 2 3
x 2 4
x 2 5
x 2 6
x 2 7
x 2 8
x 1 x 2
x 1 x 3
x 1 x 4
x 1 x 5
x 1 x 6
x 1 x 7
x 1 x 8
x 2 x 3
x 2 x 4
x 2 x 5
x 2 x 6
x 2 x 7
x 2 x 8
x 3 x 4
x 3 x 5
x 3 x 6
x 3 x 7
x 3 x 8
x 4 x 5
x 4 x 6
x 4 x 7
x 4 x 8
x 5 x 6
x 5 x 7
x 5 x 8
x 6 x 7
x 6 x 8
x 7 x 8
x 1 x 2 x 3
表2 因子水平编码表
Table 2 Level code for factors
x i
+2.045
+1
0
-1
-2.045
I
24 135
21 000
18 000
15 000
11 865
N =2p +2p +1=28 +2×8+1=273
中心常数为:
K = (2p /N ) 1/2 =0.968 4
星号臂为:
r = (2p -2N ) 1/2 -2p -1=2.045
上述方程式中, N 是需要计算的点数目; p 是因素的数目; k 是中心常数; r 是星号臂。
如果全盘考虑所有项, 即包括所有因素的相互作用, 则方程式中会出现很多项, 那样太复杂而且没必要。 从实际应用和分析中, 我们可以忽略大多数高次幂项。 这里我们只采用一次项和二次项, 且只用一个三次项比较其重要性。 稍后检查出的重要性表明所有三次项不严重影响电流效率。
电流与因子的转换关系为:
x i = (I -18 000) /3 000
式中 I 为电流强度, 0水平为18 kA, +1水平为21 kA。 实验期间, 电解铝厂的系列电流大约在150 kA左右。 当然, 建立的电流效率和电流分布的关系可以延伸到更宽的范围。
电流效率和区域电流的代数表达式如下:
η =B 0 +∑a i x i +∑b i x
2 i
+
∑ (x i ∑C ij x j ) (3)
因为很难测量区域电流效率来验证得到的数学模型, 我们仅考虑整个电解槽的平均电流效率。 用于计算平均电流效率的B 0 , a i , b i 和c ij 在表3中给出。 实验期间, 电解质成分的分子比为2.4, Al2 O3 3%, MgF2 3%, CaF2 3%左右, 系列电流为约150 kA。 表3的系数可能会随操作环境的不同而不同。
为了测试正交多元回归模型的精确度, 任选一些点分别用正交回归模型和综合模型计算电流效率。 结果如图2所示。 从图2的结果来看, 我们认为正交回归模型可代替复杂综合模型来分析电流分布对电流效率的影响。
2电流分布对电流效率的影响
由式 (3) 很容易计算出平均电流效率和区域电流效率, 图3所示为各个区域电流对平均电流效率的影响, 曲线上的数 (如1, 2等) 表示相关区域。
图3中, 整个电流被固定在153 kA。 我们使区域电流从10 kA变化到30 kA, 这是符合铝电解过程阳极作业的特点。 在实验期间, 实际区域电流变化在13~24 kA范围。 在图3的分析中, 只让一个区域的电流变化, 总电流保持连续, 剩下的电流平均分布在其他区域。 图3结果表明, 对大多数区域, 在考察的范围内存在最大电流效率, 但有些区域, 例如3和5, 电流效率在给定范围内一直增加 (曲线5) 或一直减少 (曲线3) 。 我们认为各区域的电流分布对电流效率有不同的影响, 多数区域存在最佳电流分布, 这是一个普遍存在的规律。 当然, 各种电解槽的定量关系不相同, 这与电解槽构造有关, 与电解槽物理场分布密切相关。
表3 用回归方程 (3) 计算平均电流效率的系数表
Table 3 Coefficients for computation of even current efficiency by Eqn.3
Coefficient
Value
Coefficient
Value
Coefficient
Value
Coefficient
Value
Coefficient
Value
a1 a2 a3 a4 a5 a6 a7 a8 b1 b2
0.256 0.183 -0.074 0.076 0.191 0.032 0.069 -0.062 -0.051 -0.069
b3 b4 b5 b6 b7 b8 c12 c13 c14 c15
0.021 -0.045 0.003 -0.030 -0.034 -0.042 -0.024 -0.004 0.025 0.000
c16 c17 c18 c23 c24 c25 c26 c27 c28 c34
0.052 0.041 0.017 0.034 -0.025 -0.001 0.007 0.049 0.059 -0.027
c35 c36 c37 c38 c45 c46 c47 c48 c56 c57
0.008 0.048 -0.006 -0.040 0.036 0.017 0.067 -0.101 -0.031 0.024
c58 c67 c68 c78 c123 B0
0.081 -0.047 0.005 0.025 0.001 0.863
图2 回归模型和综合模型的计算结果比较
Fig.2 Calculation results for bothcomprehension model and regression model1—Results for comprehension model;2—Results for regression model
图3 电流分布对电流效率的关系图
Fig.3 Relation between currentdistribution and current efficiency
由方程式 (3) 出发, 也很容易作出二元图和多元图。 下面就最佳的电流效率作一简要描述。
3 最佳电流效率
什么是最佳的电流效率, 怎样的电流分布能得到最高的电流效率?这些问题很难精确解答, 很可能是一个多门学科的综合而归结为一个数学难题。 在表3中, 最高电流效率为91.78%, 这不是在整个范围和多因素综合考虑的最大值, 只是在特定条件下的最佳值。 虽然很难得到整个区域中的最高电流效率, 但易知对所有阳极均匀分布的电流得不到最大电流效率。 如果所有阳极电流是相等的, 则电流效率等于91.12%。
基于Kuhn Tucker 理论
[15 ]
, 能得到当∑x i =3 (总电流153 kA) 时电流效率的区域最大值, 过程如下: 引进函数φ (x i , λ ) =-η (x i ) +λ (∑x j -3)
令
{ ? φ / ? λ = 0 ? φ / ? x i = 0
所以
{ ∑ x j - 3 = 0 - a i - 2 b i x i - ∑ C i j x j + λ = 0
解上述方程得x 1 =0.788; x 2 =0.653; x 3 =1.885; x 4 =-1.245; x 5 =1.284; x 6 =0.845; x 7 =-1.241; x 8 =0.31 (注: I =3 000× x i +18 000; 例如: x 4 =-1.245, I =14 265A) 。 在这种情况下, 平均电流效率为91.17%, 比电流平均分布在各个区域时的电流效率大0.05%。 而在图3中, 有许多点的电流效率都大于91.12%。 这似乎与传统的观点有点矛盾。 由于电解槽中存在不平衡的磁场分布, 不均匀的电流分布部分补偿了不均匀的磁场分布, 减少了电磁力对熔体的作用, 从而减低了金属流动速度, 提高了电流效率。 如何才能达到最佳的电流效率?由于电流分布明显地影响电流效率, 因此从电流分布的观点, 应该尽量接近最佳电流分布。 新换的阳极不能象其他阳极一样导电, 而要在10 h以后才能正常工作。 如果新阳极能设置在最佳电流分布很低的地方, 而且在它正常导电时移到别处, 那么电流效率就能增加。 这样操作可能会增加人力成本, 但若每区的最佳电流分布差别很大的话, 这是值得的。 另外, 组装阳极后, 按组装质量分级, 在安装阳极前, 把导电性好的阳极安装在最佳电流分布较高的区域。 如果条件许可, 电解槽上安装调节装置, 也是一种有效的方法。
4 结论
以160 kA预焙电解槽为例, 应用铝电解过程电流效率的综合模型通过正交多元回归方法对电流效率的计算, 得到一个正交多元回归模型, 此模型直接描述了预焙铝电解槽电流效率与电流分布间的关系。 同时, 应用简单的正交多元回归方程式分析了电流分布对电流效率的影响。 结论为: 在某些限定条件下, 对大多数阳极的电流分布有一个最佳电流效率值; 但在不同条件下最佳电流分布值不一样, 这与电解槽结构和物理场分布相关; 且在电解槽投入生产后, 对特定的电解槽, 建立特定的电流分布对电流效率的关系方程, 对于指导电解槽作业, 提高生产效率有参考意义。
参考文献
[1] 邱竹贤.预焙槽炼铝[M].北京:冶金工业出版社, 1988.6.QIUZhu xian.AluminumMetallurgywithPrebakeCell[M].Beijing:MetallurgicalIndustryPress, 1988.
[2] GrjotheimK, KrohnC, MalinovskyM, etal.AluminumElectrolysis[M]. (2nded) .Dusseldorf:Aluminum Verlag, 1982.
[3] AureliuP, AugustinM, GheorgheD, etal.Analysisofthephenomenaatthestart upofaluminumelectrolysiscellbythevisualizationofthecurrentdistribution[A].ProceedingsoftheTechnicalSessionsatthe131stAnnualMeeting[C].184ThornHillRoadWarrendale, Pennsylvania, Washington:APublicationofTMS, 2002.389396.
[4] ArkhipovGV, PinginVV.Investigatingthermoelectricfieldsandcathodebottomintegrityduringcellpreheating, start upandinitialoperatingperiod[A].ProceedingsoftheTechnicalSessionsatthe131stAnnualMeeting[C].184ThornHillRoadWarrendale, Pennsylvania, Washington:APublicationofTMS, 2002.347354.
[5] ZHOUPing, ZHOUNai jun, MEIChi, etal.NumericalcalculationandindustrialmeasurementsofmetalpadvelocitiesinHall Heroultcells[J].TransNonferrousMetSocChina, 2003, 13 (1) :208212.
[6] 吴建康, 黄 珉, 黄 俊, 等.铝电解槽电解质铝液流动及铝液表面变形计算[J].中国有色金属学报, 2003, 13 (1) :241244.WUJian kang, HUANGMin, HUANGJun, etal.Computationofflowfieldofelectrolyte aluminiumliquidandsurfacedistortionofaluminumliquidinreductioncell[J].TheChineseJournalofNonferrousMetals, 2003, 13 (1) :241244.
[7] 曾水平, 彩祺凤, 梅 炽, 等.铝电解槽内磁场的三维数值分析[J].中南工业大学学报, 1995, 26 (5) :618622.ZENGShui ping, CAIQi feng, MEIChi, etal.Calculationandanalysisofthemagneticfieldinaluminareductioncell[J].JournalofCentralSouthUniversityofTechnology, 1995, 26 (5) :618622.
[8] 曾水平, 刘业翔, 梅 炽.铝电解槽内三维电流分布数值模拟[J].有色金属, 1996, 48 (3) :8892.ZENGShui ping, LIUYe xiang, MEIChi.3 Dcurrentdistributioninaluminareductioncell[J].NonferrousMatals, 1996, 48 (3) :8892.
[9] 刘业翔, 梅 炽, 曾水平.铝电解槽内三维力场的数值模拟[J].中国有色金属学报, 1996, 6 (1) :2731.LIUYe xiang, MEIChi, ZENGShui ping.Computersimulationofelectromagneticforcefieldonthemeltinaluminareductioncell[J].TheChineseJournalofNonferrousMetals, 1996, 6 (1) :2731.
[10] ZENGShui ping.Effectofcurrentdistributiononcurrentefficiencyin160kAprebakecells[A].ProceedingsoftheTechnicalSessionsatthe131stAnnualMeeting[C].184ThornHillRoadWarrendale, Pennsylvania, Washington:APublicationofTMS, 2002.503510.
[11] HaupinWE.ProductionofAluminiumandAlumina[M].London:LondonPress, 1987.
[12] ZENGShui ping, LIUYe xiang, MEIChi.Mathematicalmodelforcontinuousdetectionofcurrentefficiencyinaluminumproduction[J].TransNonferrousMetSocChina, 1998, 8 (4) :683687.
[13] 曾水平, 刘业翔.160kA预焙铝电解槽区域电流效率[J].中国有色金属学报, 2000, 10 (2) :274277.ZENGShui ping, LIUYe xiang.Zonecurrentefficiencyin160kAprebakealuminumelectrolysiscell[J].TheChineseJournalofNonferrousMetals, 2000, 10 (2) :274277.
[14] 汪锡孝.实验研究方法[M].长沙:湖南科技出版社, 1989.WANGXi xiao.ExperimentalMethods[M].Changsha:HunanScienceTechnologyPress, 1989.
[15] FoxRL.工程设计最优方法[M].北京:科学出版社, 1984.FoxRL.OptimuminEngineeringDesign[M].Beijing:SciencePress, 1984.