钨矿强化碱煮搅拌釜的放大设计研究.pdf

5 0 有色金属 冶炼部分 h t t p ／／y s y l ．b g r i m m ．c n 2 0 1 5 年第7 期 d o i 1 0 ．3 9 6 9 I ／i ．i s s n ．1 0 0 7 - 7 5 4 5 ．2 0 1 5 ．0 7 ．0 1 4 钨矿强化碱煮搅拌釜的放大设计研究 刘飞飞，王兵，林飞，罗穆峰 江西理工大学，江西赣州3 4 1 0 0 0 摘要以钨矿机械活化碱煮搅拌釜试验模型为对象，建立了几何相似放大的搅拌釜模型，并基于F l u e n t 软件对模型进行模拟分析。结果表明，放大的搅拌釜在顶部两端区域易形成流动死区，不利于固液混 合。优化该放大模型的搅拌釜结构和搅拌桨安装方式后能避免产生流动死区。改变搅拌釜模型的容积 和搅拌桨转速，对比分析优化前后搅拌釜内流体的混合特性和搅拌功率，建立了混合特性评价体系，同 等条件下，优化后的搅拌釜能缩短混合时间，提高混合效率。 关键词钨矿；机械活化；设计；反应釜；放大；数值模拟 中图分类号T F 8 4 1 ．1 ；T F 8 0 3 ．2文献标志码A文章编号1 0 0 7 7 5 4 5 2 0 1 5 0 7 0 0 5 0 0 5 S c a l e - u pD e s i g no fM e c h a n i c a lA c t i v a t i o nA l k a l iB o i l i n gR e a c t o rf o rT u n g s t e nO r e L I UF e i f e i ，W A N GB i n g ，L I NF e i ，L U OM u f e n g J i a n g x iU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g y ，G a n z h o u3 4 1 0 0 0 ，J i a n g x i ，C h i n a A b s t r a c t B a s e do nm o d e lo fm e c h a n i c a la c t i v a t i o na l k a l ib o i l i n gr e a c t o rf o rt u n g s t e no r e ，t h er e a c t o rm o d e l w a se s t a b l i s h e d b ya d o p t i n gg e o m e t r i cs i m i l a r i t ya m p l i f i c a t i o nm e t h o da n da n a l y z e d b y n u m e r i c a l s i m u l a t i o nw i t hF l u e n ts o f t w a r e ．T h er e s u l t ss h o wt h a ta m p l i f i e dm o d e lt e n d st op r o d u c ef l o w “d e a d z o n e s ”o nb o t hs i d e sa tt h et o po fr e a c t o r ，w h i c hi sn o tc o n d u c i v et os o l i d l i q u i dm i x t u r e ．F l o w “d e a d z o n e s ”d i s a p p e a r sf o l l o w i n go p t i m i z a t i o no fr e a c t o r ’Ss t r u c t u r ea n di m p e l l e ri n s t a l l a t i o n ．T h em i x i n g c h a r a c t e r i s t i c se v a l u a t i o ns y s t e mw a se s t a b l i s h e dw i t ha d j u s t m e n to fs t i r r i n gs p e e da n dv o l u m eo fr e a c t o r ， a n dc o m p a r a t i v ea n a l y s i so nm i x i n gc h a r a c t e r i s t i co ff l u i da n ds t i r r i n gp o w e ro fr e a c t o r ．T h er e s u l t ss h o w t h a to p t i m i z e dr e a c t o rc a ns h o r t e nm i x i n gt i m ea n di m p r o v em i x i n ge f f i c i e n c y ． K e yw o r d s t u n g s t e no r e ；m e c h a n i c a la c t i v a t i o n ；d e s i g n ；r e a c t o r ；s c a l e u p ；n u m e r i c a ls i m u l a t i o n 由于传统低品位、高杂质含量的钨矿原料浸出 工艺存在浸出率低、能耗高、效率低等缺点，难以满 足市场对金属钨的需求，改善钨冶炼碱煮浸出工艺 迫在眉睫，为此，李洪桂等[ 1 3 研究了钨的机械化活化 强化浸出工艺，有效解决了上述问题，但是该工艺的 主要设备 热球磨机 磨损快，限制了该设备的推广 应用，作者[ 2 。33 与赣州市某公司合作研发设计了集磨 矿与浸出为一体的钨矿机械活化碱煮设备试验样 机，该样机能有效避免磨损，且工艺具有高效、节能、 收稿日期2 0 1 5 一O l 一1 6 基金项目国家自然科学基金资助项目 6 1 3 6 4 0 1 4 作者简介刘飞飞 1 9 6 2 一 ．男，江西人，博士，教授． 流程短等优点，不足的是，该样机所设计的容积较 小，不能满足生产需求，为进一步投入生产应用，需 对该设备进行中试和扩大试验。 近年来随着计算机技术的发展，运用C F D 方法 模拟分析搅拌釜内流场已成为一种非常有效的手 段，如运用数值模拟方法对搅拌釜内固液悬浮流场 进行分析[ 4 坩] ，刘丽艳等[ 7 ] 应用C F D 方法模拟计算 了粉煤灰氢氧化钠桨液槽内各项性能参数的变化 情况，并计算得到一种适合于粉煤灰氢氧化钠桨 万方数据 2 0 1 5 年第7 期有色金属 冶炼部分 h t t p ／／y s y l ．b g r i m m ．c n 5 1 液槽固液悬浮的放大判据。O s h i n o w o [ 8 3 用C F D 方 法分析了反应槽内的速度分布，J a w o r s k i [ 9 ] 利用F l u e n t 软件模拟计算了双层涡轮桨的混合过程。因此， 本文采用C F D 方法模拟计算几何相似放大的钨矿机 械活化碱煮反应釜，运用F l u e n t 软件模拟分析该放大 反应釜内的流体特性，并对其结构进行优化。 1C F D 方法 运用软件F l u e n t 对搅拌设备釜内流场进行仿 真分析，通常可采用的数值模拟方法包括“黑箱”模 型法、动量源法、内外迭代法、多重参考系法、滑移网 格法[ 1o ] 等。为了解决釜内静止挡板与旋转搅拌桨 之间的交互作用，通常采用内外迭代法、多重参考系 法和滑移网格法。三者的共同点在于搅拌釜内均 分为两个区域，即搅拌桨桨叶区域和静止挡板与釜 体组成的桨外区域。三者不同之处在于相比于内 迭代法和滑移网格法，多重参考系法主要用于模拟 两区域之间具有微弱的相互作用，且两区域之间的 交互作用采用时均解法，其次，该方法计算量小，求 解计算耗时短，计算结果较精确，应用比较广泛，因 此，本文采用多重参考系法对搅拌釜内流体进行模 拟分析。 2 数值模拟及分析 2 ．1 搅拌釜模型的建立和物性参数 建立了几何相似放大的搅拌釜模型，结构简图 如图1 所示，搅拌釜主要由釜体、导流筒、挡板、强化 装置、J H 型轴流搅拌桨[ 1 2 ] 等组成，釜内均布四块挡 板，搅拌桨转速为2 4 0r ／m i n 。搅拌釜尺寸依照试验 样机实际尺寸2 倍建立，其结构尺寸如下 m m 釜 径D 一6 4 0 、桨径d 一2 6 8 、槽体圆柱高度H 一7 6 0 、桨 离釜底高度h 一2 0 8 、挡板宽度W e 一6 0 、槽体椭球短 轴直径d 1 3 2 0 。 两相物系为氢氧化钠溶液和白钨矿 小颗粒 状 的混合物，固液比为1 1 ．1 ，固体颗粒直径为 7 5 “m 。 2 ．2 网格划分 由于搅拌釜内结构复杂，故选取四面体非结构 化网格进行网格化处理，为了减少耗时，得到较精确 的计算结果，网格划分采用不同的网格尺度，搅拌桨 旋转部分生成较密的四面体网格，图2 为搅拌釜网 格示意图。 2 ．3 边界条件和模型求解设置 模型采用多重参考系法求解计算，需将搅拌釜 图1 搅拌浸出槽结构简图 F i g ．1D i a g r a mo fr e a c t o rs t r u c t u r e 图2 搅拌浸出槽网格示意图 F i g ．2 G r i dd i a g r a mo fr e a c t o r 划分为动区域和静区域，定义搅拌桨桨叶附近区域 为动区域，其他区域为静区域，动区域和静区域之间 产生交界面”i n t e r f a c e ”。模型采用稳态非耦合求解 器，多相流模型为混合模型，相数为2 ，湍流模型采 用标准肛e 模型。残差精度设为1 1 0 ～，计算时监 视力矩趋于稳定值即为计算收敛。 2 ．4 搅拌槽内流场分析 图3 为z 一0 平面上对应的速度矢量图。 从图3 可看出，搅拌釜内的流体形成了轴流循 环，在搅拌桨桨叶端附近，流体轴向速度较大，沿着 桨叶至釜壁的方向，流体速度逐渐减小，撞击釜壁 后，沿着釜壁向上流动，再流经强化装置，实现流体 的循环流。不足的是，由于搅拌桨的“抽吸”作用力 过大，导流筒外流体受到搅拌桨的强吸作用，大部分 流体未流向釜顶就转向流至导流筒内，因此，在釜顶 万方数据 5 2 有色金属 冶炼部分 h t t p ／] y s y l ．b g r i m m ．c n 2 0 1 5 年第7 期 ≯繇 ≯ 工、 图3 搅拌釜z 0 截面速度矢量图 F i g ．3V e l o c i t yv e c t o rd i a g r a mo f r e a c t o ra tz20s e c t i o n 两端的流体流动不明显，易形成“流动死区”，且处在 该区域的流体由于重力的影响，有向下流动的趋势， 与沿釜壁向上流动的流体形成对流，减弱了向上流 动的流体轴向速度，且在靠近釜壁、搅拌桨上方区域 产生漩涡流，上述现象均不利于物料的混合，在很大 程度上延长了混合时间，降低了混合效率。 出现上述现象的原因可能是采用几何相似放大 的模型，依靠相似原理[ 1 引，即外形结构参数对应成 比例，但体积成幂次增加，在相同转速条件下，搅拌 桨单位体积功率降低，转矩动力偏小，随着搅拌桨的 尺寸增大，搅拌桨旋转运动时产生的抽吸力越大，导 流筒内压力变小，而流体往往从高压区流向低压区， 故流体过早地流向导流筒。依此可以推测，随着釜 径D 的增大，搅拌釜内“流动死区”区域容积增大， 混合效率下降，搅拌功耗增加，同时增加了制造成 本，对工业应用搅拌釜而言，显然是不经济的。 3 搅拌釜结构的优化 3 ．1 模型结构参数 从上述分析可知，几何相似法直接放大的搅拌 釜存在一些不足，因此，对搅拌釜的结构作如下优化 设计 1 保持搅拌釜体积不变，降低搅拌釜的高度； 2 釜内安装2 个桨径不大的搅拌桨。根据搅拌 釜体积不变，可计算出优化后的搅拌釜模型的釜径 D 和槽体椭球体短轴直径d 。。其结构尺寸如下 r a m D 一9 0 0 、d 。一2 2 8 、槽体圆柱高度H 一3 8 0 、 桨离釜底高度h 一1 2 0 、挡板宽度W b 一3 0 、两桨间 距h 1 3 7 0 。 3 ．2 仿真结果分析 在相同转速条件下，运用F l u e n t 软件对搅拌釜再 次数值模拟，图4 为搅拌釜z 0 截面的速度矢量图。 4 7 e 0 I 耵J e “l l 3 2 e { } L 1 8 e “儿 1 0 e “ 1 0 3 e l l l 5 7 e “抑 8 3 e “M l O e { H } 3 6 e H l 6 3 e “HJ 8 9 e “且l 1 6 e lx l 4 2 e “m 8 6 e H l O S e O 1 2 l e “K J 4 8 e I H l c e 一 ” 3 6 e 一 1 j 图4 优化后搅拌釜z O 截面速度矢量图 F i g ．4V e l o c i t yv e c t o rd i a g r a mo fr e a c t o ra t z - - - - - 0s e c t i o na f t e ro p t i m i z a t i o n 由图4 可知，流体亦产生了循环流，在搅拌桨下 方和上方附近的流体流动速度较大，流体经搅拌作 用，同时甩向反应釜中心区域和釜壁，在两搅拌桨的 作用下，流体的轴向速度变大，向上流动时速度衰减 慢，能流向更高的高度，避免釜顶产生“流动死区”。 且搅拌桨的尺寸小、轻便、易于安装。 3 ．3 不同转速下的混合性能分析 为了得出优化前后搅拌釜内流体混合性能的优 劣性，在不同转速条件下对搅拌釜进行混合时间的 模拟测定，分析优化前后搅拌釜的单位体积混合能 的变化。根据功的计算公式，可推导出单位体积混 合能的表达式为 W 。一P ，t 1 式中，w 。为单位体积混合能 J ／m 3 ；P ，为单 位体积功率 W ／m 3 ；t 为混合时间 s ，可通过混合 时间模拟计算得出。 P 。计算表达式[ 1 4 1 如下 P 。一P ／V 一2 n n M ／V 2 式中，P 为搅拌功率 w ；竹为搅拌桨转速 r ／ m i n ；M 为搅拌扭矩 N m ；V 为釜体体积 m 3 。 运用F l u e n t 软件计算出搅拌釜的转速分别为 2 1 0 、2 4 0 、2 7 0 、3 0 0r ／m i n 时，优化前的迭代扭矩值 分别为0 ．9 4 、1 ．2 3 、1 ．5 6 、1 ．9 3N m ；优化后的迭 代扭矩值分别为1 ．1 1 、1 ．4 4 、1 ．6 8 、2 ．2 5N m 。对 浏豳固豳圜阑阀翻阀㈤㈣黼湖矧川川二] 叫 叫 一 一 一 一 一 一 一 ■■_纛熏燃篷澄■■羹鍪■_■● 万方数据 2 0 1 5 年第7 期 有色金属 冶炼部分 h t t p f ／y s y l ．b g r i m m ．c n 5 3 应上述转速下优化前的混合时间分别为 s 8 6 ．5 、 8 1 ．4 、7 5 ．2 、7 2 ．9 ；优化后的混合时间分别为 S 7 1 ．3 、6 7 ．6 、6 1 ．5 、5 6 ．1 。 搅拌釜总体积V 约为1 7 0L ，将上述数据代入 公式 2 中，即可求出单位体积功率P ，。结果如下 搅拌釜的转速分别为2 1 0 、2 4 0 、2 7 0 、3 0 0r ／r a i n 时，优 化前的单位体积功率分别为1 1 7 、1 8 2 、2 5 9 、3 5 6w ／ m 3 ；优化后分别为1 4 3 、2 1 3 、3 0 0 、4 1 6w ／m 3 。 可知，随着转速的增大，搅拌釜内物料混合时间 减小，转速增大有利于物料的混合。单位体积功率 随着转速的增大而增加，从而增加了功率消耗。从 上述分析可知，转速增大虽然能缩短混合时间，但同 时增大了功率消耗，因此，转速不宜过大。相比于优 化前的搅拌釜，优化后的搅拌釜在双搅拌桨的搅拌 作用下，搅拌力矩大，虽然单位体积功率大，但混合 时间短。 将上述数据带入公式 1 可求出单位体积混合 能W 。，由此得到了如图5 所示的单位体积混合能 W ，与转速7 1 的关系曲线。 乒 g ● 晕 色 _ _ 罐 姐 赠 器 递 g 暑 图5 转速与单位体积混合能的关系曲线 F i g ．5 V a r i a t i o nc u r v eb e t w e e nr o t a t i n g s p e e da n du n i tm i x i n ge n e r g y 由图5 可知，随着搅拌桨转速的增大，流体的单 位体积混合能增大，在相同转速条件下，相对于优化 前的搅拌釜，优化后的搅拌釜单位体积混合能降低。 综合上述分析可知，相对于优化前的搅拌釜，优 化后的搅拌釜缩短了物料混合时间，降低了搅拌釜 的总功耗，混合效率更高。 3 ．4 不同体积下的搅拌功率 搅拌功率是衡量搅拌釜经济性能的重要依据， 在转速行一2 4 0r ／m i n 时，不同的搅拌釜体积条件 下，采用F l u e n t 软件求解出搅拌釜体积分别为2 0 、 1 7 0 、5 7 0 、13 5 0L 时的扭矩值，并分析优化前后搅拌 釜的搅拌功率的变化规律。已知扭矩M 和转速翘， 由公式 2 可求出不同体积下的搅拌功率P ，结果 如图6 所示。从图6 可以看出，随着搅拌釜容积的 增大，搅拌功率增大，相对于优化前的搅拌釜，优化 后的搅拌釜功耗更低，且搅拌功率增长速度慢。 芝 碍 罨 被 藉 图6 体积一搅拌功率图 F i g ．6 V a r i a t i o nc u r v eb e t w e e nv o l u m e a n ds t i r r i n gp o w e r 4结论 1 通过运用F l u e n t 模拟分析钨矿机械活化碱 煮搅拌釜模型，发现了优化前的搅拌釜在釜顶两端 易形成“流动死区”，对搅拌釜的结构和搅拌桨的安 装方式进行优化设计后，避免了在釜顶产生“流动死 区”，流体轴向速度较大。 2 优化后的搅拌釜混合时间和单位体积混合能 低于优化前，有效提高了混合效率。 3 搅拌功率随搅拌釜容积的增大而增加，在同 等条件下，优化后的搅拌釜功耗低于优化前。 参考文献 [ 1 ] 李洪桂，刘茂盛，思泽金，等．白钨精矿与黑白钨混合矿 碱分解方法及设备中国，8 5 1 0 0 3 5 0 [ P ] ．1 9 8 6 0 8 2 7 ． E 2 ] 刘飞飞，贺强．带双层桨的钨矿碱煮浸出槽内流场仿真 分析[ J ] ，湿法冶金，2 0 1 3 ，3 2 2 1 0 1 1 0 4 ． [ 3 ] Z H A N GK e - w e i ，L I UF e i f e i ．T h eS t r u c t u r eO p t i m i z a t i o nD e s i g no ft h eR e a c t i o nT a n kw i t hA g i t a t o r [ c ] ／／ T h e3 r dI n t e r n a t i o n a lC o n f e r e n c eo nM e c h a n i cA u t o m a t i o na n dC o n t r o lE n g i n e e r i n g ，B a o t o u ，2 0 1 2 1 9 8 6 1 9 8 9 ． E 4 ] 李永纲，黄雄斌．立式圆槽内多轴搅拌器固一液悬浮性 能[ J ] ．过程工程学报，2 0 1 2 ，1 2 2 1 8 l 一1 8 6 ． [ 5 ] 贾海洋，唐克伦，唐永亮，等．搅拌槽内部固一悬浮流场 的数值模拟及实验研究[ J ] ．机械设计与制造，2 0 1 3 ，1 2 万方数据 5 4 有色金属 冶炼部分 h t t p ／／y s y l ．b g r i m m ．c n 2 0 1 5 年第7 期 0 5 1 。f 0 1 j 2 ． [ 6 ] 王振松，李良超，黄雄斌．固一液搅拌槽内槽底流场的 C F D 模拟[ J ] ．北京化工大学学报自然科学版，2 0 0 5 ， 3 2 4 5 - 9 ． [ 7 ] 刘丽艳，李晨辰，朱国瑞，等．粉煤灰一氢氧化钠浆液槽 放大准则数值模拟研究[ J ] ．现代化工，2 0 1 4 ，3 4 2 1 4 8 1 5 1 ． [ 8 ] O s h i n o w oL ，J a w o r s k iz ，D y s t e rKN ，e ta 1 ．P r e d i c t i n g t h et a n g e n t i a lv e l o c i t yf i e l di ns t i r r e dt a n k su s i n gt h e m u l t i p l er e f e r e n c ef r a m e M R F m o d e lw i t hv a l i d a t i o n b yL D Am e a s u r e m e n t s [ C ] ／／P r o c e e d i n go f1 0 t hE u r o p e a nC o n f e r e n c eo nM i x i n g ，D e l f t ，2 0 0 0 2 8 1 - 2 8 8 ． [ 9 ] J a w o r s k iZ ，B u j a l s k iW ，O t o m oN ，e ta 1 ．C F DS t u d yo f H o m o g e n i z a t i o nw i t hD u a lR U S H T O NT u r b i n e sC o r n p a r i s o nw i t hE x p e r i m e n t a lR e s u l t s [ J ] ．T r a n sI C h e m E ， 2 0 0 0 ，7 8 A 3 2 7 3 3 3 ． [ 1 0 ] L u oJY ，G o s m a nAD ，I s s aRI ，e ta 1 ．F u l lf l o wf i e l d c o m p u t a t i o no fm i x i n gi nb a f f l e ds t i r r e dr e a c t o r s [ J ] ． T r a n sI C h e m E ，1 9 9 3 ，7 1 3 4 2 3 4 4 ． [ 1 1 ] P e r i c l e o u sKA ，P a t e lM K ．T h es o u r c e - s i n ka p p r o a c h i nt h em o d e l i n go fs t i r r e dr e a c t o r s E J ] ．P h y s i c oC h e m i c a lH y d r o d y n a m i c s ，1 9 8 7 ，9 2 7 9 2 9 7 ． [ 1 2 ] 张锁龙，沈惠平，张国忠．J H 型轴流式搅拌桨流场分 析及设计口] ．化学工程，1 9 9 9 ，2 7 5 2 6 2 9 ． [ 1 3 ] 晋日亚，贺增第，朱志伟．关于化工设备放大方法的讨 论[ J ] ．化学工程师，1 9 9 9 ，3 4 6 2 5 2 7 ． [ 1 4 ] 何珊珊，李振花，万茂荣，等．搅拌槽中固液悬浮及其 放大的研究[ J ] ．化学工业与工程，1 9 9 5 ，1 2 3 7 - 1 4 ， 3 5 ． 上接第4 9 页 4结论 采用细菌氧化一氰化提金工艺处理含砷含碳难 处理金矿的最佳条件为矿石粒度一0 ．0 7 4r n m 占 8 1 ％、氧化矿浆温度3 0 ℃、矿浆p H 1 ．6 左右、矿 浆浓度2 0 ％、细菌氧化时间4 天，金浸出率从直接 氰化的不足5 ％提高到9 3 ．8 1 ％。 参考文献 E l l 陈宁，姚英杰，程蓉，等．微生物分解难处理金矿的动力 学研究[ J ] ．有色金属 冶炼部分 ，2 0 0 8 1 3 8 4 1 ． [ 2 ] 王仍坚，蒋开喜，王海北，等．合成硫化铜的细菌浸出研 究[ J ] ．有色金属 冶炼部分 ，2 0 1 2 1 1 1 6 1 8 ． [ 3 1 郭持皓，李云，王云，等．从难处理含砷金矿中湿法回收 砷[ J ] ．有色金属 冶炼部分 ，2 0 1 2 5 8 - 1 0 ． [ 4 ] 杨松荣，邱冠周，胡岳华，等．含砷难处理金矿石生物氧 化工艺及应用[ M ] ．北京冶金工业出版社，2 0 0 6 ． [ 5 ] 张建涛，程矗，陈帅，等．脱除某难处理金矿中砷和硫的 研究[ J ] ．有色金属 冶炼部分 ，2 0 1 4 5 1 2 1 6 ． [ 6 ] R a w l i n g sDE ，J o h n s o nDB ．B i o m i n i n g ．S p r i n g e r ， 2 0 0 7 ． r 7 ] F r a n z m a n nP ，H a d d a dC7H a w k e sR ，e ta 1 ．E f f e c t so f t e m p e r a t u r eo nt h er a t e so fi r o na n ds u l f u ro x i d a t i o nb y s e l e c t e db i o l e a c h i n gB a c t e r i aa n dA r c h a e a A p p l i c a t i o n o ft h eR a t k o w s k ye q u a t i o n [ J ] ．M i n e r a l sE n g i n e e r i n g ， 2 0 0 5 ，1 8 1 3 1 3 0 4 1 3 1 4 ． [ 8 ] D i x o nDG ．A n a l y s i so fh e a tc o n s e r v a t i o nd u r i n gc o p p e r s u l p h i d eh e a pl e a c h i n g [ J ] ．H y d r o m e t a l l u r g y ，2 0 0 0 ，5 8 i 2 7 4 1 ． [ 9 ] Y a n gC ，L iY ，Z h o uY ，e ta t ．B a c t e r i a lc o m m u n i t yd y n a m i c sd u r i n gab l o o mc a u s e db yA k a s h i w os a n g u i n ei n t h eX i a m e ns e aa r e a ，C h i n a [ J ] ．H a r m f u la l g a e ，2 0 1 2 ， 2 0 1 3 2 1 4 1 ． [ 1 0 ] R a w l i n g sDE ．C h a r a c t e r i s t i c sa n da d a p t a b i l i t yo fi r o n - a n ds u l f u r o x i d i z i n gm i c r o o r g a n i s m su s e df o rt h er e c o v e r yo fm e t a l sf r o mm i n e r a l sa n dt h e i rc o n c e n t r a t e s [ J ] ． M i c r o b i a lc e l lf a c t o r i e s ，2 0 0 5 ，4 1 1 3 ． [ 1 1 ] F o u c h e rS ，B a t t a g l i a - B r u n e tF ，d ’H u g u e sP ，e ta 1 ．E v o l u t i o no ft h e b a c t e r i a lp o p u l a t i o nd u r i n gt h eb a t c hb i o l e a c h i n go fac o b a h i f e r o u sp y r i t ei nas u s p e n d e d - s o l i d s b u b b l ec o l u m na n dc o m p a r i s o nw i t ham e c h a n i c a l l ya g i t a t e dr e a c t o r [ J ] ．H y d r o m e t a l l u r g y ，2 0 0 3 ，7 1 1 5 - 1 2 ． [ 1 2 ] Z h o uQG ，1 3 0F ，B oZH ，e ta 1 ．I s o l a t i o no fas t r a i no f A c i d i t h i o b a c i l l u sc a l d u sa n di t sr o l ei nb i o l e a e h i n go f e h a l c o p y r i t e [ J ] ．W o r l dJ o u r n a lo fM i c r o b i o l o g ya n dB i o t e c h n o l o g y ，2 0 0 7 ，2 3 9 1 2 1 7 1 2 2 5 ． 万方数据