资源描述:
2 0 1 5 年第3 期 有色金属 冶炼郏分 h t t p //y s y l .b g r i m m .c n 2 3 d o i 1 0 .3 9 6 9 /j .i 蟠n .1 0 0 7 7 5 4 5 .2 0 1 5 .0 3 .0 0 7 大型铝电解槽电流效率数值建模研究 邱泽晶1 ,张红亮2 ,徐宇杰2 ,王振宇1 1 .南瑞 武汉 电气设备与工程能效测评中心,武汉4 3 0 0 7 4 ; 2 .中南大学冶金与环境学院,长沙4 1 0 0 8 3 摘要对铝电解槽内电流损失机理和现代大型槽 4 0 0k A 的强湍流特征进行剖析,采用界面更新理论 描述铝液电解质界面附近的传质行为,建立了电流损失计算方程、质量传递系数方程及电流效率全域 计算方程,进而结合多相流体力学及熔盐电化学相关理论,在C F X 平台上建立了关联熔体运动的大型 铝电解槽电流效率全域数值计算模型。应用该模型对某4 0 0k A 级大型铝电解槽的电流效率进行了瞬 态实时计算,求得全槽时均电流效率约为9 2 .5 7 %。 关键词铝电解槽;电流效率;数值模拟;多相流 中图分类号T F 8 2 1文献标志码A文章编号1 0 0 7 7 5 4 5 2 0 1 5 0 3 0 0 2 3 一0 7 N u m e r i c a lM o d e l i n go fC u r r e n tE f f i c i e n c yi nL a r g e - s c a l e A l u m i n u mR e d u c t i o nC e U s Q I UZ e j i n 9 1 ,Z H A N GH o n g l i a n 9 2 ,X UY u j i e 2 ,W A N GZ h e n y u l 1 .N A R I W u h a n E l e c t r i c a lE q u i p m e n ta n dE n g i n e e “n gE f f i c i e n c y E v a l u a t i o nC e n t e r ,W u h a n4 3 0 0 7 4 。C h i n a 2 .S c h 0 0 1o fM e t a l l u r g ya n dE n v i r o n m e n t ,C e n t r a lS o u t hU n i v e r s i t y ,C h a n g s h a4 1 0 0 8 3 ,C h i n a A b s t r a c t T h em e c h a n i s mo fc u r r e n t1 0 s sa n ds t r o n gt u r b u l e n c ec h a r a c t e r i s t i c si nl a r g e s c a l ea l u m i n u mc e l l s ≥ 4 0 0k A w e r ea n a l y z e d .M a s st r a n s f e rn e a r b ym e t a l - b a t hi n t e r f a c ew a sd e s c r i b e dw i t hs u r f a c er e n e w a l m o d e l .T h ec u r r e n tl o s sc a l c u l a t i o ne q u a t i o n ,m a s st r a n s f e rc o e f f i c i e n te q u a t i o na n dc u r r e n te f f i c i e n c y C E c a l c u l a t i o ne q u a t i o nw e r eb u i l t .W i t hc o m b i n a t i o no fm u l t i p h a s ef l u i dm e c h a n i c sa n dm 0 1 t e ns a l t e l e c t r o c h e m i s t r yt h e o “e s , ag l o b a ln u m e r i c a lc a l c u l a t i o nm o d e lo fC Ea s s o c i a t e dw i t hm e l tf l o ww a s e s t a b l i s h e d .T h i sm o d e lw a se f f e c t i v e l ya p p l i e dt oa n a l y z eg l o b a la n dl o c a lC Ei n4 0 0k A g r a d ea l u m i n u m r e d u c t i o nc e l l s .T h ec a l c u l a t e dt i m e a v e r a g e dC Ei s9 2 .5 7 %. K e yw o r d s a l u m i n u mr e d u c t i o nc e l l ;c u r r e n te f f i c i e n c y ;n u m e r i c a ls i m u l a t i o n ;m u l t ip h a s ef l o w 当前,我国已建和在建大型铝电解槽 4 0 0 k A 超过2o o o 万t ,所消耗的电能超过全国总发电 量的5 %,节能降耗已成为维持我国铝电解工业可 持续发展的基础问题。为了实现大幅节能、提高劳 动生产率及降低单位产能的投资成本,近年来,铝电 解槽正向低电压、大容量方向迅猛发展,4 0 0k A 级 以上电解槽在我国铝工业界快速兴起,以实现低电 压为目标的各种新型结构电解槽及临界稳态生产工 艺相继面世5 。。这一发展方向使我国铝工业实现 进一步节能降耗的目标成为可能,但同时也带来了 收稿日期2 0 1 4 一l O 0 8 基金项目国家自然科学基金青年科学基金资助项目 5 1 1 0 4 1 8 7 ;国家自然科学基金面上项目 5 1 2 7 4 2 4 1 ;国家自然科学基金 创新群体资助项目 6 1 3 2 1 0 0 3 作者简介邱泽晶 1 9 8 4 一 。男,湖北鄂州人,硕士,高级工程师;通信作者张红亮 1 9 7 9 一 ,男,湖南临湘人,博士,副教授. 万方数据 2 4 有色金属 冶炼部分 h t t p //y s y l .b g r i m m .c n 2 0 1 5 年第3 期 新的挑战,即如何兼顾低电压与高电流效率。为此, 建立大型铝电解槽的全域电流效率的数学模型,全 面掌握电流效率的分布及变化情况对当前铝电解工 业、尤其是以低电压为代表的节能型铝电解工艺来 说,具有重要的意义。 电流效率是铝电解生产最重要的技术经济指标 之一,直接关系电耗指标的高低,目前已有较多的相 关研究报道,具体到电解槽电流效率的计算方面,也 涌现出不同的计算模型,总体而言,可分为两种一 种是经验计算模型[ 6 ] ,即根据大量现场实测数据建 立相关工艺参数与电流效率之间的经验函数关系; 另一种是理论计算模型[ 7 ] ,即根据电流损失的机理, 基于电化学原理建立其数学描述。这两种方法皆有 优缺点,为此,有学者[ 8 。9 1 尝试结合这两方面的研究, 在假设熔体为静止的前提下,给出了电流效率计算 模型。此外,还有学者[ 1 叩将铝液一电解质界面视为 固定壁面,从而在电流效率的计算中可以考虑流场 影响。由此可见,针对铝电解槽内电流效率的计算 模型,尽管已有一些理论与经验模型,但皆不完善, 尤其是未能充分考虑电解槽流场的影响,因而有较 大的缺陷,无法为现代大型低电压节能型铝电解槽 的电流效率预测提供可靠的支持。 本文针对大型铝电解槽内熔体运动剧烈、电流 效率影响因素复杂的特点,结合电流损失机理及相 关多相传质理论,建立了基于瞬态多相流的电流效 率全域计算模型,并应用所建模型对某4 0 0k A 级 系列铝电解槽的电流效率进行了分析,获得该槽内 电流效率的具体分布。 1 铝电解槽电流效率全域计算模型的 建立 铝电解槽内电流流经阴极而没有产出铝,就造 成了电流的损失。电流损失主要是铝的二次反应损 失和钠的二次反应损失,具体而言,包括铝的二次 反应损失、钠的二次反应损失、A l 。C 。的生成和氧 化、电子导电性、杂质引起的损失、阴极和阳极之间 的瞬时短路、铝或钠渗透进电解槽的内衬材料而引 起的电流损失。 在已知铝电解槽内阴极电流密度分布的情况 下,槽内每一区域的电流效率 ∞ 可按下式计算 刀一堑≥1 0 0 % 1 ‘∞f 式中,i 。,为该区域阴极炭块电流密度 A /m 2 ; i 如,,为该区域损失电流密度 A /m 2 。 显然,要计算电解槽内铝液一电解质界面上的 电流效率分布,必须知道该界面上的阴极电流密度 及损失的电流密度,前者可以通过电一磁一流场耦 合计算得到,后者则需要应用相关经验公式得到,而 为了模型与实际一致,需考虑多个影响因素,其中主 要为流场的某些变量。 通过对铝电解槽内电流损失的机理分析,电解 过程中电流损失的电流密度可近似应用式 2 计算 求得L 8 J i 妇一F 忌N 。c N 。[ 一a e x p 一i ≮u ] 2 式中,F 为法拉第常数 9 64 8 5C /m 0 1 ;是N 。为铝 液一电解质界面处钠单质的质量传递系数;T 为熔 体温度 12 2 3K ;己,为浓差极化压降 V ;f 。。为密 闭体系内钠的平衡浓度 m o l /L ;口为钠在电解质 体系内相对于c 。。值的饱和度。 闭体系内钠的平衡浓度,可根据实验数据及相 关经验关系式得到,而浓差极化电势U 可基于如下 方程[ 1 u 建立 u o .0 0 1 5 R 一与箬等 3 式中R 为分子比,计算时取2 .3 。 该电流效率计算理论模型中,除了钠的质量传 递系数走N 。外,其他变量均为已知,考虑到大型工业 铝电解槽中钠的传质行为受到熔体运动的显著影 响,造成了走N 。值在铝液一电解质界面上的不同区 域和不同时刻都存在差异,故只要计算得到是。。的 局部分布,便可得到铝液一电解质界面全域空间上 电流效率的分布。 为了计算界面处忌。。的值,本文应用界面更新理 论[ 1 2 。胡来描述铝液一电解质界面处钠的传质行为, 并推导得到钠的时均局部质量传递系数忌N 。 ‰一 等存 i 4 式中,e 为湍动能耗散率 m 2 s - 3 ;D 为扩散 系数,计算中取1 1 0 _ 8m 2 s;u 为运动黏度 m 2 /s ; 以多相磁流体计算为基础,结合式 1 至式 4 可求得大型槽内电流效率 田 的空间与时间的分布。 有关多相流的计算的理论及建模过程可参见前期相 关研究的结果[ 1 4 。18 。。 2 大型铝电解槽电流效率算例与结果 分析 以某4 0 0k A 级铝电解槽为算例,应用前文所 万方数据 2 0 1 5 年第3 期有色金属 冶炼部分 h t t p //y s y l .b g r i m m .c n 2 5 述理论建立电流效率全域数值计算模型,结果如图 1 所示,并进行电流效率大小及分布进行分析。 图l某4 0 0k A 铝电解槽多相流及 电流效率模型计算网格 F i g .1 C a I c u I a t i o nm e s ho fm u l t ip h a s e sa n d c u r r e n te f f i c i e n c ym o d e lf o r4 0 0k Ac e 儿 该4 0 0k A 级槽采用单阳极设计,共4 8 块阳极, 2 4 块阴极炭块,并用非通长双钢棒出电,大面加工 距离为2 8 0m m ,小面加工距离为3 9 0m m ,主要结 构参数电流强度4 2 0k A 、电解槽尺寸1 77 8 0m m 43 2 0m m 14 2 2m m 、槽膛尺寸1 75 4 0m m 41 4 0m m 5 6 0m m 、阴极炭块尺寸36 8 0m m 6 6 5 m m 4 8 5m m 、阴极钢棒尺寸45 2 0m m 2 0 0m m x1 0 0m m 、阳极尺寸17 0 0m m 6 6 0m m 6 3 5 m m 、钢壳厚度1 6m m 、阴极钢棒个数2 4 2 。此外, 根据实际生产信息,取铝水平2 7 0m m 、电解质水平 2 0 0m m 、极距4 5m m 。计算中所涉及到的熔体主要 物性参数为铝液密度22 7 0k g /m 3 、黏度1 .1 8 1 0 | ’P a s ;电解质密度21 3 0k g /m 3 、黏度2 .5 1 1 0 卅P a s ;气泡密度o .3 9 8k g /m 3 、黏度5 .0 6 1 0 4P a s 。 通过应用本课题组前期所开发的铝电解槽多物 理场模型与方法,首先计算得到稳态4 0 0k A 级铝 电解槽的电磁场分布结果,其阴极电流密度匕分布 如图2 所示 铝液上表面 。 、,n 1 二j jl ‘J I一 f } “ I 一, 1 q i 一5 j 卅 一l ‘”朋 一} I ~}一.{ ‘ I一‘{ { - 228 { J { 一22 曩 图24 0 0k A 级铝电解槽阴极电流密度分布 F i g .2 C a t h o d ec u r r e n td e n s i t yd i s t r i b u t i o no f4 0 0k Ac e U 由图2 可看出,槽内阴极电流密度的分布呈现 一定的规律,其中阳极中缝、边缝及间缝对应区域 阴极电流密度相对较低、而阳极底掌对应区域则 较高。 应用所开发的三维瞬态三相磁流体模型计算得 到电解槽瞬时的界面波动与湍流耗散率的结果,为 方便起见,本文仅选取流场处于相对稳定的1 0 0s 的时间周期,在该时间内,铝液一电解质界面上的平 均湍动能耗散率的变化情况如图3 所示。 由图3 可知,由于铝液一电解质界面的瞬时波 动及流场的不停运动,处于铝液一电解质界面处的平 均湍动能耗散率亦随着时间在o .0 0 9 ~o .0 1 5m 2 s - 3 的区间内波动,该时间内的总均值为1 .1 3 6 1 0 _ 2 m 2 s ~。 在得到电流密度与湍动能耗散率之后,便可应 用式 4 计算该1 0 0s 周期内铝液一电解质界面处 1 2 1 厂 z l ⋯j }8 | ⋯Ⅲ,㈨L ⋯”⋯1 ■了_ _ 百i i _ 一百~百t , } ;t I { | 7 、 图3 不同时刻算例中铝电解槽铝液一电解质 界面上平均湍动能耗散率 F i g .3A V e r a g et u r b u l e n c ee d d yd i s s i p a t i o n a td i f f e r e n tt i m e 钠的传质系数随时间与空间的变化,结果如图4 所示。 万方数据 .2 6 有色金属 冶炼部分 h t t p //y s y l .b g r i m m .c n 2 0 1 5 年第3 期 一一 结合式 1 ~ 3 计算得到该4 0 0k A 级铝电解槽 l a l c 在不同时刻电流效率的空间分布,结果如图5 所示。 、/ m m ’㈣ 嘶 } 6 1 8 f 124 H2 a 1 0s ; b 2 0s ; c 3 0s ; d 4 0s ; e 5 0s ; f 6 0s ; g 7 0s ; h 8 0s ; i 9 0s ; j 1 0 0s 图4 算例中铝液一电解质界面钠传质系数瞬态空间分布结果 F i g .4 S o d i u mm a s st r a n s f e rc o e f f i c i e n t0 fm e t a 卜b a t hi n t e r f a c ea td i f f e r e n tt i m e 万方数据 2 0 1 5 年第3 期 有色金属 冶炼部分 h t t p //y s y l .b g r i m m c n 2 7 C E l l H jH ,∽H 8 8 l lH t Jl l q J .q lq 2l l q {∽q 4 q I 1 %I J 町 a 1 0s ; b 2 0s ; c 3 0s ; d 4 0 s ; e 5 0s ; f 6 0s ; g 7 0s ; h 8 0s ; i 9 0s ; j l o o s 图5 算例中电流效率瞬态空间分布 F i g .5 D i s t r i b u t i o no fc u r r e n te f f i c i e n c y o fm e t a 卜b a t hi n t e r f a c ea td i f f e r e n tt i m e 万方数据 2 8 有色金属 冶炼部分 h t t p //y s y l .b g r i m m .c n 2 0 1 5 年第3 期 由图4 可知,铝液~电解质界面上钠的传质系 数的平均值在一定范围内扰动,且其具体分布也存 在一定的时域变化,但其总体分布形态则基本保持 不变。其规律可总结如下首先,传质系数的变化基 本与阳极底掌的投影相重合,这表明了气泡对于传 质系数的重大影响;其次,在靠近进电侧中间部分与 出电侧的两个端部界面传质系数相对较大,造成这 一现象的原因与槽内熔体的运动状态相关。 由图5 可见,槽内电流效率在空间分布上同样 会呈现一定的规律,在阳极底掌、阳极间缝与中缝区 域,由于钠的传质系数存在较大的差异,这些区域的 电流效率亦存在较明显的差异,说明电解槽内熔体 的运动对于电解槽内电流效率具有较大影响。 为了得到某时刻全槽的电流效率,同样在1 0 0s 时间间隔内对电流效率进行空间域的平均,得到全 槽平均电流效率随时间的变化规律如图6 所示。 时间,s 图6 算例铝电解槽在不同时刻的 全槽电流效率 F i g .6A V e r a g ef u c e l lc u r r e n te f f i c i e n c y a td i f f e r e n tt i m e 由图6 可知,全槽电流效率的平均值是随时间 的变化而变化,这一规律与实际生产高度一致,整体 波动范围在9 2 .4 %~9 2 .7 %,属于幅度较小的扰 动。通过对电流效率计算结果的时域的平均,可得 到该4 0 0k A 级电解槽在所取时间段内的电流效率 的平均值为9 2 .5 7 %。 通过上述4 0 0k A 级计算实例可知,大型铝电 解槽的电流效率大小及分布受槽内磁流体行为的影 响显著,在当前铝电解槽普遍实行低电压节能技术 的前提下,可通过适当改变影响磁流体行为的各种 因素 包括母线配置、阴极形状、阳极形状、槽内衬结 构与工艺条件等 来提高电解槽的整体电流效率。 3结论 1 建立了基于槽内熔体瞬态运动的大型铝电解 槽电流效率全域数值计算模型。即基于电解槽内铝 液一电解质界面附近瞬态湍动能耗散率可计算得到 电流损失,进而得到铝液一电解质界面处的瞬态电 流效率的分布结果。 2 应用所建电流效率全域模型对某4 0 0k A 级 系列铝电解槽的电流效率进行了分析,获得了该槽 内电流效率分布的具体特征,最后求得全槽时均电 流效率约为9 2 .5 7 %。 参考文献 [ 1 ] 杨小华,李贤.2 4 0k A 铝电解槽内衬结构优化[ J ] .有色 金属 冶炼部分 ,2 0 1 4 9 2 6 2 9 。 [ 2 ] 赵应彬,王平,陈谦.3 0 0k A 新式异型阴极双钢棒铝电 解槽生产实践[ J ] .有色金属 冶炼部分 ,2 0 1 4 2 2 8 3 0 . [ 3 ] 王轩,赵仁涛,王永良,等.3 5 0 k A 铝电解槽焙烧自动分 流装置的设计与实现[ J ] .有色金属 冶炼部分 ,2 0 1 4 4 1 4 1 7 . [ 4 ] 马琼,刘进县.5 0 0k A 预焙阳极电解槽针振和摆动现象 的处理[ J ] .有色金属 冶炼部分 ,2 0 1 4 1 2 1 1 一1 4 . [ 5 ] 秦卫中,王有来,李勇,等.大型铝电解槽低电压物料与 能量双平衡控制技术[ J ] .有色金属 冶炼部分 ,2 0 1 3 4 二2 7 3 1 . [ 6 ] D e w i n gEw . L o s so fc u r r e n te f f i c i e n c yi na l u m i n u m e l e c t r o l y s i sc e i l s [ J ] .M e t a l l u r g i c a lT r a n s a c t i o n sB , 1 9 9 1 ,2 2 B 1 7 7 1 8 2 . [ 7 ] P e t e r s o nRD ,w a n gX . I n f l u e n c eo fd i s s o l v e dm e t a l si n c r y o l i t i cm e l t so nH a l lc e l lc u r r e n ti n e f f i c 涎n c y [ C ] //R o o yEL .L i g h tM e t a l s1 9 9 1 .N e w0 r l e a n s ,L A T M S , 1 9 9 1 3 3 1 3 3 7 . [ 8 ] s t e r t e nA ,S o l l iPA . A ne l e c t r o c h e m i c a lc u r r e n te f f i c i e n c ym o d e lf o ra l u m i n i u me l e c t r 0 1 y s i sc e l l s [ j ] .J o u n a l o fA p p l i e dE l e c t r o c h e m i s t r y ,1 9 9 6 ,2 6 1 8 7 1 9 3 . [ 9 ] S t e r t e nA . C u r r e n te f f i c i e n c yj na l u m i n i u mr e d u c t i o n c e l l s [ J ] ,J o u n a lo fA p p l i e dE l e c t r o c h e m i s t r y ,1 9 8 8 ,1 8 4 7 3 4 8 3 . [ 1 0 ] H a a r b e r gT ,s 0 1 h e i mA ,J o h a n s e nsT ,e ta 1 .E f f e c to f a n o d i cg a sr e l e a s eo nc u r r e n te f f i c i e n c yi nH a I I H e r o u l t c e l l s [ c ] //w e l c hBJ .L i g h tM e t a l s1 9 9 8 . S a nA n t o n i o ,T X T M S ,19 9 8 4 7 5 4 8 1 . [ 1 1 ] H a a r b e r gGM , s e nKs ,T h o n s t a dJ ,e ta 1 .M e a s u r e m e n to fe I e c t r o n i cc o n d u c t i o ni n c r y o l i t ea l u m i n am e l t s a n de s t i m a t i o no fi t se f f e c to nc u r r e n te f f i c i e n c y [ J ] . M e t a l l u r g i c aJT r a n s a c t i o n sB ,1 9 9 3 ,2 4 B 7 2 9 7 3 5 . 0 O O 0 0 0 醉鞍毽西 万方数据 2 0 1 5 年第3 期有色金属 冶炼部分 h t t p ;//y s y Lb g r i m m .c n 2 9 [ 1 2 ] Ja j u e eB ,M a r g a r i t i sA ,K a r a m a n e vD ,e ta 1 .A p p l i c a t i o n o fs u r f a c er e n e w a I s t r e t c hm o d e lf o ri n t e f f a c em a s s t r a n s f e r [ J ] . c h e m i c a lE n g i n e e r i n gs c i e n c e ,2 0 0 6 ,6 l 3 9 1 7 . [ 1 3 ] c y r 订c o u r t e s s o l e ,J a c q u e l i n eE t a y . F I o w sa n dm a s s t r a n s f e r si nt w os u p e r i m p o s e di i q u i d1 a y e r si na ni n d u c t i o nf u r n a c e [ J ] .I n t e r n a t i o n a lJ o u r n a lo fH e a ta n d M a s sT r a n s f e r ,2 0 1 3 ,6 5 1 0 8 9 3 . [ 1 4 ] L IJ i e ,X uY u j i e ,Z H A N GH 。n g _ 1 i a n g ,e ta 1 .A ni n h o m o g e n e o u st h r e e _ p h a s em o d e If D rt h ef l o wi na l u m i n i u m r e d u c t i o nc e l l s [ J ] .I n t e r n a t i o n a lJ o u r n a lo fM u l t i p h a s e F l o w ,2 0 1 1 ,3 7 4 6 5 4 . [ 1 5 ] Z H A N GH o n g l i a n g 。L IJ i e ,W A N Gz h i g a n g ,e ta 1 . T h en u m e r i c a lm o d e l i n go fm e l tf l o w a n dM H Di n s t a b i l i t i e si na na l u m i n u mr e d u c t i o nc e I l [ J ] .J O M ,2 0 1 0 ,6 2 1 1 2 6 3 1 . [ 1 6 ] Z H A N GH o n g l i a n g ,Z H O Uc HENNQ ,w uB i n g ,e t a 1 .Av i r t u a la l u m i n u mr e d u c t i o nc e l l [ J ] .J O M ,2 0 1 3 , 6 5 1 1 1 4 5 2 1 4 5 8 . [ 1 7 ] X uY u - j i e ,Z H A N GH o n g l i a n g ,L IJ i e ,e ta 1 .An o n l i n e a rs h a n o w w a t e rm o d e lc o m b i n e dw i t hg a sb u b b l e e f f e c tf o rm e l tf l o w sa n di n t e r f a c ei n s t a b i l i t vi na l u m i n u mr e d u c t i o nc e l l s [ J ] .J O M ,2 0 1 3 ,6 5 1 1 1 4 5 9 一 1 4 6 6 。 [ 1 8 ] z H A N GH o n g .1 i a n g ,Y A N Gs h u a i ,Z H A N GH } h u i ,e t a 1 . N u m e r i c a ls i m u l a t i o no fa l u m i n a m i x i n gp r o c e s s w i t ham u I t i c o m p o n e n tf l o wm o d e lc o u p l e dw i t he l e c t r D m a g n e t i cF o r c e si na l u m i n u mr e d u c t i o nc e l l s [ J ] . J O M ,2 0 1 4 ,6 6 7 1 2 1 0 一1 2 1 7 . 上接第2 2 页 [ 6 ] 周洪武,铅阳极泥冶炼技术简评和电热连续熔炼的可 行性[ J ] .有色冶炼,2 0 0 2 ,3 1 4 7 1 1 . [ 7 ] 胡建辉,张传福.铜阳极泥预处理脱铜工艺优化[ J ] .贵 金属,2 0 0 2 ,2 3 4 1 5 . [ 8 ] 王日.铜阳极泥处理工艺优化[ J ] ,南昌水专学报, 2 0 0 4 ,2 3 4 7 5 7 8 . [ 9 ] 王爱荣,陈志刚,涂百乐.卡尔多炉处理铜阳极泥的生 产实践[ J ] .有色金属 冶炼部分 ,2 0 1 4 8 1 8 2 1 , 万方数据
展开阅读全文