适应风电消纳的铝电解阳极气泡行为建模.pdf

适应风电消纳的铝电解阳极气泡行为建模 ① 何桂雄１， 熊 敏１， 张红亮２， 唐艳梅１， 江 南２ （１．中国电力科学研究院，北京 １００１９２； ２．中南大学 冶金与环境学院，湖南 长沙 ４１００８３） 摘 要 以适应风电消纳的异型阳极为研究对象，采用计算流体力学（ＣＦＤ）中的多气泡族（ＭＵＳＩＧ）模型对异型阳极表面生成的气 泡在熔融电解质中的运动行为进行计算。 结果显示，在阳极气泡的尺寸分布上，普通阳极底掌被所有尺寸气泡覆盖，而倾角与倒角 阳极底掌则以小尺寸气泡为主，且随着倾角或倒角增大，大尺寸气泡在异型阳极底掌覆盖范围大幅减小；在气泡逸出效果上，由于 倒角阳极的气泡逸出路径要比普通或倾角阳极短，不易形成局部大尺寸气泡群，因而倒角阳极的排气效果要比倾角阳极更好。 关键词 新能源消纳； 阳极气泡； 多气泡族（ＭＵＳＩＧ）模型； 倾角阳极； 倒角阳极 中图分类号 ＴＦ８２１文献标识码 Ａｄｏｉ１０．３９６９／ ｊ．ｉｓｓｎ．０２５３－６０９９．２０１６．０６．０１９ 文章编号 ０２５３－６０９９（２０１６）０６－００７５－０４ Ｓｔｕｄｙ ｏｎ Ｂｕｂｂｌｅ Ｂｅｈａｖｉｏｒ ｏｆ Ａｂｎｏｒｍａｌ Ａｎｏｄｅ ｄｕｒｉｎｇ Ｅｌｅｃｔｒｏｌｙｓｉｓ ｏｆ Ａｌｕｍｉｎａ ＨＥ Ｇｕｉ⁃ｘｉｏｎｇ１， ＸＩＯＮＧ Ｍｉｎ１， ＺＨＡＮＧ Ｈｏｎｇ⁃ｌｉａｎｇ２， ＴＡＮＧ Ｙａｎ⁃ｍｅｉ１， ＪＩＡＮＧ Ｎａｎ２ （１．Ｃｈｉｎａ Ｅｌｅｃｔｒｉｃ Ｐｏｗｅｒ Ｒｅｓｅａｒｃｈ Ｉｎｓｔｉｔｕｔｅ， Ｂｅｉｊｉｎｇ １００１９２， Ｃｈｉｎａ； ２．Ｓｃｈｏｏｌ ｏｆ Ｍｅｔａｌｌｕｒｇｙ ａｎｄ Ｅｎｖｉｒｏｎｍｅｎｔ， Ｃｅｎｔｒａｌ Ｓｏｕｔｈ Ｕｎｉｖｅｒｓｉｔｙ， Ｃｈａｎｇｓｈａ ４１００８３， Ｈｕｎａｎ， Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ Ｔｈｅ ｍｏｖｅｍｅｎｔ ｏｆ ｂｕｂｂｌｅ ｇｅｎｅｒａｔｅｄ ｏｎ ａｎ ａｂｎｏｒｍａｌ ａｎｏｄｅ ｓｕｒｆａｃｅ ｉｎ ｔｈｅ ｍｏｌｔｅｎ ｅｌｅｃｔｒｏｌｙｔｅｓ ｗａｓ ｃａｌｃｕｌａｔｅｄ ｂｙ ｕｓｉｎｇ ａ ＣＦＤ⁃ｍｏｄｅｌ ＭＵＳＩＧ． Ｔｈｅ ｒｅｓｅａｒｃｈ ｏｆ ｔｈｅ ｓｉｚｅ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ａｎｏｄｅ ｂｕｂｂｌｅｓ ｓｈｏｗｓ ｔｈａｔ ｔｈｅ ｂｏｔｔｏｍ ｏｆ ｎｏｒｍａｌ ａｎｏｄｅ ｉｓ ｃｏｖｅｒｅｄ ｂｙ ａｌｌ ｓｉｚｅｓ ｏｆ ｂｕｂｂｌｅｓ， ａｎｄ ｉｎｃｌｉｎｅｄ ａｎｏｄｅ ａｎｄ ｃｈａｍｆｅｒｅｄ ａｎｏｄｅ ａｒｅ ｃｏｖｅｒｅｄ ｍａｉｎｌｙ ｂｙ ｆｉｎｅ⁃ｓｉｚｅ ｂｕｂｂｌｅ ｗｉｔｈ ｔｈｅ ｃｏｖｅｒｉｎｇ ａｒｅａ ｄｅｃｒｅａｓｉｎｇ ｓｉｇｎｉｆｉｃａｎｔｌｙ ｗｉｔｈ ａｎ ｉｎｃｒｅａｓｅ ｉｎ ｉｎｃｌｉｎｅｄ ｏｒ ｃｈａｍｆｅｒｅｄ ａｎｇｌｅ． Ｔｈｅ ｒｅｓｅａｒｃｈ ａｂｏｕｔ ｔｈｅ ｂｕｂｂｌｅ ｅｓｃａｐｉｎｇ ｅｆｆｅｃｔ ｓｈｏｗｓ ｔｈａｔ ｔｈｅ ｃｈａｍｆｅｒｅｄ ａｎｏｄｅ ｈａｓ ａ ｓｈｏｒｔｅｒ ｅｓｃａｐｉｎｇ ｐａｓｓａｇｅ ｔｈａｎ ｔｈａｔ ｏｆ ｎｏｒｍａｌ ｏｒ ｉｎｃｌｉｎｅｄ ａｎｏｄｅ， ｒｅｓｕｌｔｉｎｇ ｉｎ ｆｅｗｅｒ ｌａｒｇｅ⁃ｓｉｚｅ ｂｕｂｂｌｅ ｇｒｏｕｐｓ ｆｏｒｍｅｄ ｉｎ ｔｈｅ ｌｏｃａｌ ａｒｅａ． Ｉｔ ｉｓ ｃｏｎｃｌｕｄｅｄ， ｔｈｅ ｃｈａｍｆｅｒｅｄ ａｎｏｄｅ ｉｓ ｂｅｔｔｅｒ ｉｎ ｂｕｂｂｌｅ ｅｓｃａｐｉｎｇ ｃｏｍｐａｒｅｄ ｔｏ ｔｈｅ ｉｎｃｌｉｎｅｄ ａｎｏｄｅ． Ｋｅｙ ｗｏｒｄｓ ｎｅｗ ｅｎｅｒｇｙ ａｃｃｏｍｍｏｄａｔｉｏｎ； ａｎｏｄｅ ｂｕｂｂｌｅ； ＭＵＳＩＧ ｍｏｄｅｌ； ｉｎｃｌｉｎｅｄ ａｎｏｄｅ； ｃｈａｍｆｅｒｅｄ ａｎｏｄｅ 我国新能源产业发展迅猛，但却面临着无法有效 消纳的难题。 以风电为例，目前我国风电资源丰富的 西北地区，发展速度极快，但这些地区的电网建设相对 落后，总体可再生能源接入能力不高，导致大范围弃风 现象。 通过来自 ２０１５ 年上半年风电并网运行报告，弃 风最严重的地区包括甘肃、新疆、蒙西等地区，有接近 ３０％的弃风幅度［１］。 与此同时，我国铝电解工业是有 名的耗电大户，且我国铝电解工业发展迅猛，２０１５ 年 中国原铝产量达到 ３ ２００ 万吨，约占全球原铝总产量 的 ５６．７％［２］，消耗了达 ４ ３７８ 亿千瓦时的电能，占全社 会用电总量的 ７．９％［３］。 更严峻的是，我国铝电解行业 平均电力成本约占原铝总生产成本的 ４５％ ～５５％，电 力供应与电价已成为铝电解企业盈亏的关键因素。 考虑到我国电解铝高产能地区恰好与风电资源丰 富的地区基本重叠，通过电网与铝电解相关技术的研 发与设备升级，配合电价补偿机制，既可以解决铝电解 的电耗问题，又可以实现风电等新能源的消纳，对于铝 电解和风电两个产业皆具有极其重要的意义。 应用铝电解消纳风电，首先要考虑由此带来的电 流波动以及电解槽稳定问题，在工业过程中即要求解 决铝电解槽的“心脏” 阳极的问题。 目前工业所用 阳极多为平底炭阳极，此类阳极与电解质的润湿性较 差，且底掌面积较大，不利于原铝生产中产生的阳极气 体（主要为 ＣＯ２和 ＣＯ 的混合气体）逸出。 该缺陷在电 ①收稿日期 ２０１６－０６－０９ 基金项目 国家科技支撑项目（２０１５ＢＡＡ０１Ｂ０４） 作者简介 何桂雄（１９８４－），男，北京人，工程师，硕士研究生，研究方向为电力需求侧管理、荷网源协同控制。 通讯作者 张红亮（１９７９－），男，湖南临湘人，副教授，工学博士，研究方向为冶金过程数值仿真与系统节能。 第 ３６ 卷第 ６ 期 ２０１６ 年 １２ 月 矿矿 冶冶 工工 程程 ＭＩＮＩＮＧ ＡＮＤ ＭＥＴＡＬＬＵＲＧＩＣＡＬ ＥＮＧＩＮＥＥＲＩＮＧ Ｖｏｌ．３６ №６ Ｄｅｃｅｍｂｅｒ ２０１６ 万方数据 解槽处于电流不稳定状态时表现尤为突出，例如，当电 解槽由于消纳风电而产生系列电流波动时，电解槽内 气泡的产生将更加杂乱无章，部分阳极气泡会积聚在 阳极底掌下的电解质中而无法释放，一方面大幅减小 实际导电接触面积，使得极间电压降增大至 ０．１５ ～ ０．３５ Ｖ［４］之间；另一方面可能引起极间温度过高，甚至 导致阳极效应的频繁发生，消耗大量额外电能等［５－６］。 因此，研究铝电解槽内阳极气体的快速排出成了铝电 解消纳新能源的关键基础问题之一，当前可应用的工 业方案包括倾角阳极和倒角阳极，均能一定程度加速 气泡排出，但其内部气泡尺寸分布和逸出行为尚未知， 故未能使其效果达到最优化。 铝电解阳极气泡的研究一直是铝电解过程流场的 研究热点，国内外有大量学者开展了相关研究［７－１０］。 目前有少量对于倒角阳极的研究，但大部分集中在平 底阳极上，且都将阳极气泡的尺寸假设为均一大小，这 显然与实际存在一定的偏差。 此外，尽管异形结构对 阳极气泡的逸出有促进作用，但缺乏系统研究，尤其是 在不同尺寸的气泡行为上，基本鲜有报道。 本文在某 ３００ ｋＡ 预焙铝电解槽消纳风电的背景 下，以消纳风电的阳极气泡作为具体研究对象，采用计 算流体力学（ＣＦＤ）方法对阳极表面生成的气泡在熔融 电解质中的运动行为进行计算，并引入 ＭＵＳＩＧ 模型对 阳极气泡进行分组，运用大型商用软件 ＡＮＳＹＳ 建立异 形阳极结构模型，着重分析不同阳极结构（包括平底、 倾角与倒角阳极）的电解槽中阳极气泡在阳极底掌的 尺寸分布及逸出行为，以期探索阳极气泡在异型阳极 中的行为规律，为铝电解槽消纳风电过程中阳极的选 用提供理论依据。 １ 模型建立 １．１ 铝电解槽多相流体的控制方程 铝电解槽内为一个多场⁃多相复杂体系，主要由气 （电解过程的阳极气体）⁃液（熔融电解质和铝水）⁃固 （扩散的氧化铝颗粒和炭粒等）等构成。 当前，流体力 学理论体系中，针对多相问题的计算方法主要有欧拉⁃ 欧拉法和欧拉⁃拉格朗日法。 考虑到本文所建立的模型 仅仅包含液（熔融电解质）⁃气（阳极气泡）两相，因此，主 要采用欧拉⁃欧拉法，其连续方程与动量方程如下 ∂ ∂ｔ （ｒαρα） ＋ ∇（ｒαραＵα） ＝ ０（１） ∂ ∂ｔ （ｒαραＵα） ＋ ∇［ｒα（ ραＵα Ｕ α）］ ＝ － ｒα∇ｐα＋ ∇{ｒαμαｅｆｆ［∇Ｕα＋ （∇Ｕα）Ｔ］}＋ ＳＭα ＋ Ｍ α （２） 式中 Ｕα为第 α 相平均速度，ｍ／ ｓ；ｒα为第 α 相体积分 数； ρα为第 α 相密度，ｋｇ／ ｍ３；ＳＭα为外部作用力，Ｎ； μαｅｆｆ为黏度，Ｐａｓ； ｐα为压力，Ｐａ；Ｍα为流体间动量， ｋｇｍ／ ｓ。 １．２ 预焙阳极的多气泡族（ＭＵＳＩＧ）模型 在铝业界，通常认为阳极气泡在阳极底掌经历以 下过程产生→长大→合并→逸出。 在通有直流电的 炭阳极⁃电解质界面上所发生的电化学反应会一直产 生阳极气泡，初始属于小体积，但随着气泡的聚集与气 泡合并作用，气泡会不断长大，并在部分区域合成大尺 寸气泡，而其他尺寸依然有着小尺寸的存在，随着气泡 不断逸出⁃产生的循环，电解槽的阳极底掌下面的气泡 实际是具备多个尺寸分布的。 在计算流体力学软件 ＣＦＸ 平台中，多气泡族模型 （ｍｕｌｔｉｐｌｅ ｓｉｚｅｄ⁃ｇｒｏｕｐ， ＭＵＳＩＧ）的基本原理是将气泡的 作用间接处理，将其离散成多个小组，每个小组按照气 泡的直径分组，通过对动量方程进行简化获得流速的 解，因而，应用该模型可将气泡的作用更加真实地刻画 出来。 在当前已有的包含气泡的模型中，基本都是将 气泡尺寸视为同等大小，这显然与实际存在较大差异。 本文所建立的 ＭＵＳＩＧ 中，气泡的组分如下 ∂ ∂ｔ ｆｉ＋ ∂ ∂Ｘｊ μｊｆｊ ＝ Ｂ Ｂ － Ｄ Ｂ ＋ Ｂ Ｃ － Ｄ Ｃ （３） 式中 ｉ、ｊ 分别为第 ｉ 组和第 ｊ 组气体的分组；ＢＢ为大尺 寸气泡生成率，ｍ３／ ｓ；ＤＢ为大尺寸气泡消失率，ｍ３／ ｓ； ＢＣ为小尺寸气泡生成率，ｍ３／ ｓ；ＤＣ为小尺寸气泡消失 率，ｍ３／ ｓ。 对于气泡聚并生成率 Ｃｉｊ，采用 Ｐｒｉｎｃｅ ＆ Ｂｌａｎｃｈ 模 型加以描述，其聚并速率用下式表示 Ｃｉｊ ＝ Ｆ Ｃ θＴ ｉｊ ＋ θ Ｂ ｉｊ ＋ θ ＬＳ ｉｊ η ｉｊ （４） 式中 ＦＣ为修正系数，ＦＣ＝１； θＴ ｉｊ为湍流作用力，Ｎ；θ Ｂ ｉｊ为 气泡浮力，Ｎ；θＬＳ ｉｊ为剪切力，Ｎ；ηｉｊ为气泡聚并概率。 ＢＣ＝ １ ２ ∑ ｉ ｊ ＝ １ ∑ ｉ ｋ ＝ １ ＣｊｋＮｉＮｊ（５） 式中 ｋ 为序号；Ｎｉ和 Ｎｊ分别为第 ｉ 组和第 ｊ 组气泡的 密度，ｋｇ／ ｍ３。 根据铝电解试验及经验，本文将所研究的虚拟阳 极气泡分为 ５ 组不同大小，假定最小直径 ｄｍｉｎ为 １ ｍｍ、 最大直径为 ５ ｍｍ，其余气泡直径 ｄｉ用下式计算 ｄｉ ＝ ｄ ｍｉｎ ＋ Δｄ（ｉ － ０．５）（６） Δｄ ＝ ｄｍａｘ － ｄ ｍｉｎ Ｎ （７） 式中 Δｄ 为气泡的区域均值，ｍｍ。 则计算得到 ｄ１～ｄ５ 分别为 １．４ ｍｍ、２．２ ｍｍ、３．０ ｍｍ、３．８ ｍｍ 和 ４．６ ｍｍ。 ６７矿 冶 工 程第 ３６ 卷 万方数据 １．３ 计算模型和计算流程 本文研究的异形阳极结构如图 １ 所示，采用的计 算流程如图 ２ 所示，选用 ３００ ｋＡ 铝电解槽的阳极与电 解质体系，假定的尺寸参数为阳极尺寸 １ ６００ ｍｍ ６６０ ｍｍ ５５０ ｍｍ、 中缝 １５０ ｍｍ、间缝 ４０ ｍｍ、极距为 ４５ ｍｍ、大面加工距离 ３００ ｍｍ、电解质水平 ２２０ ｍｍ。 此外，计算过程未考虑磁场的作用。 在 ＡＮＳＹＳ 中建立 其物理模型和网格划模型如图 ３ 所示。 本文所建立模 型采用的物性参数及相关假设参见文献［１１－１２］。 θ ab c gg 图 １ 异型阳极结构示意 （ａ） 普通阳极炭块； （ｂ） 倾角阳极炭块； （ｃ） 倒角阳极炭块 1451 443 6 14514 4;6 [ANSYS] ,*, 6 [MUSIG] 14B446,A CFD;D0 140 [CFX] 图 ２ 仿真计算流程 图 ３ 模型网格划分图 ２ 异型阳极结构的气泡行为 应用前文建立的电解质⁃阳极多气泡族气泡两相 稳态模型，以单个阳极单元为研究范围，在倾斜角度 θ 分别为 １、２，倒角边长 ｇ 分别为 ３ ｃｍ、４ ｃｍ 的槽内， 分析了气泡行为。 ２．１ 倾角阳极结构的气泡行为 对比了普通阳极和倾角分别为 １、２的异型阳极 结构的气泡行为，如图 ４～６ 所示。 在垂直面上电解质 层气泡的平均直径如图 ７ 所示，各尺寸气泡分布如 表 １ 所示。 由图 ４～７ 及表 １ 可以看出，无论阳极是否 有倾角，其底部气泡尺寸的分布规律均为尺寸最小的 气泡布满整个底掌，且距离底掌距离最近，而它们又被 下面一层较厚的中大尺寸的气泡所包裹，其位置略微 图 ４ 距离普通阳极底掌 １ ｍｍ 处电解质中各尺寸气泡分布 （ａ） ｄ１＝１．４ ｍｍ； （ｂ） ｄ２＝２．２ ｍｍ； （ｃ） ｄ３＝３．０ ｍｍ； （ｄ） ｄ４＝３．８ ｍｍ； （ｅ） ｄ５＝４．６ ｍｍ 图５ 距离倾角为１的阳极底掌１ ｍｍ处电解质中各尺寸气泡分布 （ａ） ｄ１＝１．４ ｍｍ； （ｂ） ｄ２＝２．２ ｍｍ； （ｃ） ｄ３＝３．０ ｍｍ； （ｄ） ｄ４＝３．８ ｍｍ； （ｅ） ｄ５＝４．６ ｍｍ ７７第 ６ 期何桂雄等 适应风电消纳的铝电解阳极气泡行为建模 万方数据 图６ 距离倾角为２的阳极底掌１ ｍｍ电解质层中各尺寸气泡分布 （ａ） ｄ１＝１．４ ｍｍ； （ｂ） ｄ２＝２．２ ｍｍ； （ｃ） ｄ３＝３．０ ｍｍ； （ｄ） ｄ４＝３．８ ｍｍ； （ｅ） ｄ５＝４．６ ｍｍ 图 ７ 倾角阳极的电解质层气泡平均直径分布 （ａ） 普通阳极； （ｂ） 倾角 １； （ｃ） 倾角 ２ 表１ 不同倾角的阳极底掌下１ ｍｍ处电解质层中各尺寸气泡分布 倾斜角度 ／ （） 不同气泡直径（ｍｍ）下的气泡分布／ ％ ｄ１＝１．４ｄ２＝２．２ｄ３＝３．０ｄ４＝３．８ｄ５＝４．６ 普通阳极１７．８５６．５４６．２６６．３１６３．０４ １１０．８２２．０９１．８８１．８６８３．３５ ２１１．１３１．７２１．５８１．５７８４．００ 靠近阳极底掌的中部位置。 主要原因是靠近底掌区域 为反应集中区域，故气泡多为小尺寸，而随着气泡运动 与合并，则逐渐在其下部增大。 由表 １ 可以看出，随着倾角增大，中等尺寸气泡开 始增多，但大尺寸气泡反而有下降趋势，其主要原因为 气泡的逸出减少了大尺寸气泡的数量。 因此，可以说 阳极有倾角能够加速气泡的逸出，有利于减少气膜电 阻，但同时，阳极有倾角时会加速小气泡的聚并，并促 进中型气泡的累积。 ２．２ 倒角阳极结构的气泡行为 应用多尺寸两相流模型，计算了不同倒角长度的 气泡行为。 不同倒角长度阳极气泡平均直径分布如 图 ８ 所示，气泡直径统计结果如表 ２ 所示。 图 ８ 不同倒角阳极电解质层中气泡的平均直径分布 （ａ） 普通阳极； （ｂ） ａ＝３ ｃｍ； （ｃ） ａ＝４ ｃｍ 表 ２ 不同倒角阳极的气泡直径分布 倒角边长 ／ ｃｍ 气泡最大直径 ／ ｍｍ 气泡最小直径 ／ ｍｍ 气泡平均直径 ／ ｍｍ 普通阳极４．６２．０３．５ ３４．５２．４３．０ ４４．５２．０４．０ 从图 ８ 和表 ２ 可以看出，普通阳极底掌处有一层 大尺寸气膜，但在倒角长度为 ３ ｃｍ 的阳极底掌出现了 一层中等尺寸的气膜，但随着倾角长度增大，该层中等 尺寸的气泡平均直径却略有增大。 实验结果验证了阳 极倒角对气泡的排出有一定促进作用，因为倾角的出 现直接导致气泡在阳极底掌的水平距离缩短，气泡更 容易从阳极端部逸出。 但该倾角的作用有一个限度， 即存在一个最佳值，需要通过反复计算来确定。 ３ 结 论 １） 建立了一种可以刻画各种尺寸的气泡在阳极 底掌行为的多相流模型，完成了模型的理论推导，建立 了相应的网格及计算流程。 ２） 在阳极气泡的尺寸分布上，普通阳极底掌为所 有尺寸的气泡覆盖，而倾角与倒角阳极底掌则是小尺 寸为主，且随着倾角或倒角增大，大尺寸气泡在异型阳 极底掌覆盖范围大幅减小；在气泡逸出效果上，由于倒 角阳极的气泡逸出路径要比普通或倾角阳极短，不易 形成局部的大尺寸气泡群，因而倒角阳极的排气效果 要比倾角阳极更好，但倒角长度存在一个最佳值，需要 通过计算优化得到。 参考文献 ［１］ 国家能源局． ２０１５ 年上半年全国风电并网运行情况［ＥＢ／ ＯＬ］． ２０１５－７－２７． ｈｔｔｐ∥ｗｗｗ．ｎｅａ．ｇｏｖ．ｃｎ／２０１５－０７／２７／ ｃ １３４４５１６７８．ｈｔｍ． （下转第 ８３ 页） ８７矿 冶 工 程第 ３６ 卷 万方数据 ２） 粗硫酸镍的生产不仅可以降低电解液中镍含 量，还可以脱除其他杂质，对电解液的净化有一定的 益处。 参考文献 ［１］ 柯 浪，彭映林，郑雅杰． 铜电解液电积脱铜制备高纯阴极铜［Ｊ］． 矿冶工程， ２０１３，３３（１）７４－７８． ［２］ Ｗｅｉｚｈｉ Ｚｅｎｇ， Ｓｈｉｊｉｅ Ｗａｎｇ， Ｍｉｃｈａｅｌ Ｌ Ｆｒｅｅ． Ｅｘｐｅｒｉｍｅｎｔａｌ ａｎｄ Ｓｉｍｕ⁃ ｌａｔｉｏｎ Ｓｔｕｄｉｅｓ ｏｆ Ｅｌｅｃｔｒｏｌｙｔｅ Ｆｌｏｗ ａｎｄ Ｓｌｉｍｅ Ｐａｒｔｉｃｌｅ Ｔｒａｎｓｐｏｒｔ ｉｎ ａ Ｐｉｌｏｔ Ｓｃａｌｅ Ｃｏｐｐｅｒ Ｅｌｅｃｔｒｏｒｅｆｉｎｉｎｇ Ｃｅｌｌ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ ｔｈｅ Ｅｌｅｃｔｒｏ⁃ ｃｈｅｍｉｃａｌ Ｓｏｃｉｅｔｙ， ２０１６，１６３（５）Ｅ１１１－Ｅ１２２． ［３］ Ａｎａ Ｓｏｆｉａ Ｖａｒｅｌａ， Ｍａｔｔｈｉａｓ Ｋｒｏｓｃｈｅｌ， Ｔｏｂｉａｓ Ｒｅｉｅｒ， ｅｔ ａｌ． Ｃｏｎｔｒｏｌｌｉｎｇ ｔｈｅ ｓｅｌｅｃｔｉｖｉｔｙ ｏｆ ＣＯ２ｅｌｅｃｔｒｏｒｅｄｕｃｔｉｏｎ ｏｎ ｃｏｐｐｅｒ Ｔｈｅ ｅｆｆｅｃｔ ｏｆ ｔｈｅ ｅｌｅｃｔｒｏｌｙｔｅ ｃｏｎｃｅｎｔｒａｔｉｏｎ ａｎｄ ｔｈｅ ｉｍｐｏｒｔａｎｃｅ ｏｆ ｔｈｅ ｌｏｃａｌ ｐＨ ［ Ｊ］． Ｃａｔａｌｙｓｉｓ Ｔｏｄａｙ， ２０１６，２６０８－１３． ［４］ 翟秀静． 中国金属冶金学［Ｍ］． 北京冶金工业出版社， ２０１４． ［５］ 朱祖泽，贺家齐． 现代铜冶金学［Ｍ］． 北京科学出版社， ２００３． ［６］ 吴 展，李 伟，陈志华，等． 高效萃取剂 ＡＤ１００ 从粗硫酸镍溶液 中回收铜的研究［Ｊ］． 矿冶工程， ２０１３，３３（２）１０５－１０７． ［７］ Ｆｏｒｓｅｎ Ｏ， Ｈｅｔｔｕｌａ Ｅ， Ｌｉｌｌｉｕｓ Ｋ． Ｂｅｈａｖｉｏｕｒ ｏｆ Ａｒｓｅｎｉｃ， Ａｎｔｉｍｏｎｙ ａｎｄ Ｂｉｓｍｕｔｈ ｉｎ ｔｈｅ Ｓｏｌｉｄｉｆｉｃａｔｉｏｎ ａｎｄ Ｅｌｅｃｔｒｏｌｙｓｉｓ ｏｆ Ｎｉｃｋｅｌ⁃Ｏｘｙｇｅｎ⁃Ｂｅａｒ⁃ ｉｎｇ Ｃｏｐｐｅｒ Ａｎｏｄｅｓ［Ｃ］∥Ｐｈｙｓｉｃａｌ Ｃｈｅｍｉｓｔｒｙ ｏｆ Ｅｘｔｒａｃｔｉｖｅ Ｍｅｔａｌｌｕｒｇｙ． Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ａｎ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｓｙｍｐｏｓｉｕｍ Ｈｅｌｄ ａｔ ｔｈｅ ＡＩＭＥ Ａｎｎｕａｌ Ｍｅｅｔｉｎｇ， １９８５３５３－３７７． ［８］ 李志坚． 通电解废液中提取硫酸镍［Ｊ］． 江西有色金属， ２００７，２１ （４）２４－２６． ［９］ 金美珊． 膜法电解硫酸镍溶液的实验研究［Ｄ］． 兰州兰州交通大 学环境与市政工程学院， ２０１４． ［１０］ 郭学益，石文堂，李 栋，等． 采用旋流电积技术从电镀污泥中回 收铜和镍［Ｊ］． 中国有色金属学报， ２０１２，２０（１２）２４２５－２４３０． ［１１］ 李大平． 电解法制备硫酸镍工艺研究［Ｊ］． 河南化工， １９９９（８） １０－１２． ［１２］ 欧阳红英． 浓缩冷冻结晶法处理废酸液的应用［Ｊ］． 环境与开 发， ２０００，１５（３）４１－４２． ［１３］ 姚素平． 电热浓缩法生产粗硫酸镍的工艺与实践［Ｊ］． 有色冶金 设计与研究， １９９５，１６（４）２１－２７． ［１４］ 余智艳． 蒸汽间接加热浓缩生产粗硫酸镍工艺应用［Ｊ］． 有色冶 金设计与研究， １９９９，２０（１）１８－２２． ［１５］ 杨立新，黄 凯，郭学益． 蒸发氧化法制备高品质氧化镉粉末 ［Ｊ］． 矿冶工程， ２００３，２３（５）６３－６５． ［１６］ 张均杰，葛晓明． 铜电解脱铜后液提取粗硫酸镍的生产实践［Ｊ］． 中国有色冶金， ２０１１（５）１５－１８． ［１７］ 苏 峰，李敬忠，李俊标，等． 粗硫酸镍提取工艺及生产实践［Ｊ］． 铜业工程， ２０１４（２）２０－２２． ［１８］ 唐谟堂，曹 刿． 湿法冶金设备［Ｍ］． 长沙 中南大学出版社， ２００９． ［１９］ 王 润． 硫酸镍结晶过程的生产控制［Ｊ］． 中国有色冶金， ２０１１ （４）２５－２７． ［２０］ 李洪贵． 冶金原理［Ｍ］． 北京 科学出版社， ２００８． ［２１］ 郑平友，余劲松，张淑萍，等． 蒸发结晶系统传热传质规律的研究 ［Ｊ］． 科学技术与工程， ２００６，６（８）１００２－１００６． ［２２］ 崔忠圻，刘北兴． 金属学与热处理原料［Ｍ］． 哈尔滨 哈尔滨工 业大学出版社， ２００８． 􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍 （上接第 ７８ 页） ［２］ 生意社． ２０１５ 年全球原铝 １２ 月产量统计［ＥＢ／ ＯＬ］． ２０１６－０２－１４． ｈｔｔｐ∥ｗｗｗ．１００ｐｐｉ．ｃｏｍ． ［３］ 国家能源局． 国家能源局发布 ２０１５ 年全社会用电量［ＥＢ／ ＯＬ］． ２０１６－１－１９． ｈｔｔｐ∥ｗｗｗ．ｈｚｇｄ．ｃｏｍ．ｃｎ／ ｚｘ／ ｘｗ／２０１６０１１９／２６８６．ｈｔｍｌ． ［４］ 杨少华，谢宝如，杨凤丽，等． 铝电解槽阳极气泡行为研究［Ｊ］． 有 色金属科学与工程， ２０１３（３）２０－２４． ［５］ 王志刚，赖延清，刘 伟，等． 铝电解阳极气泡行为研究进展［Ｊ］． 轻金属， ２００７（５）２７－３１． ［６］ 李 劼，丁凤其，李民军，等． 预焙铝电解槽阳极效应的智能预报 方法［Ｊ］． 中南工业大学学报（自然科学版）， ２００１，３２（１）２９－ ２２． ［７］ Ｐｅｒｒｏｎ Ａ， Ｋｉｓｓ Ｌ Ｉ， Ｐｏｎｃｓｋ Ｓ． Ａｎ ｅｘｐｅｒｉｍｅｎｔａｌ ｉｎｖｅｓｔｉｇａｔｉｏｎ ｏｆ ｔｈｅ ｍｏｔｉｏｎ ｏｆ ｓｉｎｇｌｅ ｂｕｂｂｌｅｓ ｕｎｄｅｒ ａ ｓｌｉｇｈｔｌｙ ｉｎｃｌｉｎｅｄ ｓｕｒｆａｃｅ［Ｊ］． Ｉｎｔｅｒｎａ⁃ ｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｍｕｌｔｉｐｈａｓｅ Ｆｌｏｗ， ２００６，３２（５）６０６－６２２． ［８］ 铁 军，邹建成． 铝电解阳极气泡长大过程的计算机模拟［Ｊ］． 计 算机与应用化学， ２００２，１９（１）１７２－１７６． ［９］ 曾 健． 阳极炭块倒角降低阳极毛耗的实践［Ｊ］． 有色冶金节能， ２０１２，２（１）２２－２４． ［１０］ 秦晓明，孙海生，崔瑞军． “船形”阳极炭块在预焙电解槽中的应 用［Ｊ］． 轻金属， ２００９（９）４５－４６． ［１１］ 刘 伟． 铝电解槽多物理场数学建模及应用研究［Ｄ］． 长沙中 南大学冶金科学与工程学院， ２００８． ［１２］ 江 南． 若干新型电极结构的铝电解槽物理场仿真研究［Ｄ］． 长 沙中南大学冶金与环境学院， ２０１５． 􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍􀪍 􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝 􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝 􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝 􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀤝􀦝 􀦝 􀦝 􀦝 中国知网版权声明 本刊已许可中国学术期刊（光盘版）电子杂志社在中国知网及其系列数据库产品中以数字化方式 复制、汇编、发行、信息网络传播本刊全文。 该社著作权使用费与本刊稿酬一并支付。 作者向本刊提交 文章发表的行为即视为同意本刊上述声明。 矿冶工程编辑部 ２０１６ 年 １２ 月 ３８第 ６ 期潘 霞等 铜电解脱铜后液除 Ｎｉ 并制取粗硫酸镍 万方数据