賈愛迪，高丙朋，徐媛媛

(新疆大學 電氣工程學院，烏魯木齊 830047)

生物氧化預處理氧化槽內氣液混合相的密度機理模型

賈愛迪，高丙朋*，徐媛媛

(新疆大學 電氣工程學院，烏魯木齊 830047)

在新疆高寒高海拔地區(qū)生物氧化提金預處理的研究中，用機理建模法建立了氧化槽內氣液混合相的密度機理數(shù)學模型。經過對氧化槽內氣液兩相的混合密度進行數(shù)值計算，擬合出礦區(qū)環(huán)境溫度、礦區(qū)大氣壓強以及礦漿質量濃度對氧化槽內氣液兩相的混合密度的影響曲線，提出了一種對氧化槽設備的改進方案。結果表明，氧化槽內氣液兩相的混合密度與礦區(qū)環(huán)境溫度以及礦漿質量濃度成正比例關系，而與礦區(qū)大氣壓強成反比例關系；改進后的氧化槽設備可提高生物氧化反應的速度和效率。

生物氧化；金礦；氣液混合相；混合密度；密度機理模型

生物氧化提金工藝受到多種因素的影響和制約，是涉及熱力學、動力學、生物化學以及電化學等多門學科的復雜工藝[1-2]，生物氧化預處理工藝是有效提高生物氧化提金效率的前期處理工藝技術，其影響因素[3-5](如溫度、礦漿pH值、氧化還原電位(ORP)、空氣量、磨礦細度以及礦漿濃度等)的分析研究在不斷的深入，已經得到了很多最優(yōu)參數(shù)范圍，例如，生物氧化預處理工藝中礦漿質量濃度的最優(yōu)范圍為15%～20%[6]；最有利于氧化槽內生物氧化反應的背景溫度為42℃等[7]。目前，對生物氧化預處理氧化槽內溫度場的建模與數(shù)值模擬的研究在漸漸深入，也有學者進行了一些細菌氧化機理、動力學過程的研究[8-9]，而對于生物氧化預處理氧化槽內密度場的建模與數(shù)值模擬仍然有待完善，因此，必須研究生物氧化預處理過程中氧化槽內氣液混合相的密度，從而為槽內密度場的建立與分析奠定基礎。

本文通過建立氧化槽內氣液兩相的密度機理模型，從而得到高寒環(huán)境、晝夜溫差較大的工況[3,10]下氣液兩相的混合密度的影響參數(shù)以及變化趨勢，為生物氧化預處理氧化槽內密度場的數(shù)值模擬奠定了基礎；同時，該模型也為研究生物氧化預處理過程中氧化還原電位(ORP)與氧化槽內氣液混合相密度之間的影響關系提供了理論支撐，為生物氧化提金預處理工藝控制系統(tǒng)參數(shù)的優(yōu)化提供理論基礎。

1 生物氧化提金預處理工藝

生物氧化預處理工藝是利用細菌作催化劑對包裹目的元素的礦石進行氧化，從而使得非目標元素被氧化成離子進入溶液中被處理回收，以便于得到暴露在氧化渣中的目標元素進行下一步處理，大大提高了金的浸出率。

圖1所示為生物氧化預處理氧化槽，槽體直徑10 m，高10 m，一般液位控制在9 m(溢漿口高度)。生物氧化預處理過程中空壓機不斷的向氧化槽內的液相礦漿通入空氣，由電動調節(jié)閥調節(jié)進氣量，空氣從氧化槽底部進入，經攪拌電機攪拌后均勻分布于礦漿中，如此，在氧化槽內部呈現(xiàn)復雜的氣液兩相混合反應狀態(tài)[11-12]。

圖1 生物氧化預處理過程氧化槽工藝Fig.1 Oxidation tank artwork of bio-oxidation pretreatment

2 氧化槽內氣液混合相密度機理數(shù)學模型

生物氧化預處理氧化槽內是礦漿溶液與所通入空氣同時存在的氣液兩相流，由于槽內有攪拌電機在進行攪拌，可以視槽內氣液兩相混合液是均勻的，因此，其氣液混合相的密度為：

其矩陣形式為：

式(2)中，ρ為氧化槽內氣液混合相的混合密度；ρg為氧化槽內氣相密度；ρl為氧化槽內礦漿溶液密度；α為氧化槽內氣相含率。

2.1 解析氣相密度ρg

氣相密度定義為氧化槽內氣相質量與氧化槽內氣液兩相體積之比。在生物氧化預處理過程中，向氧化槽內通入一定質量空氣的同時，槽內生物化學反應消耗了一定質量的氧氣，由此，可以得到理論氣相密度：

式(3)中，Q為槽內單位時間通入空氣的質量；Qg為槽內生物化學反應單位時間所消耗的氧氣質量；H為槽內液面高度；v為槽內氣相的氣升速度；V為槽內氣液兩相體積。

2.2 解析液相密度ρl

生物氧化預處理過程中氧化槽內的液相為調配好的一定質量濃度的礦漿，所占體積仍然是氣液兩相體積，因此，氧化槽內液相礦漿的密度為：

式(4)中，ms為氧化槽內含砷難處理金礦石的質量；Cl為經配礦調漿后的礦漿的質量濃度。

2.3 解析氣相含率α

氧化槽內氣液兩相在攪拌電機作用下均勻混合，同時存在，理論上氣液混合相中氣相含率：

式(5)中，Vg為氧化槽內所含氣相體積；P為當?shù)卮髿鈮毫?；Pb為生物化學反應背景壓力；T為當?shù)丨h(huán)境溫度(開氏溫度)；Tb為生物化學反應背景溫度(開氏溫度)。

將解析公式(3)、(4)、(5)均代入式(1)進行合并整理，從而得到了生物氧化預處理氧化槽內氣液混合兩相的密度機理模型：

通過分析生物氧化預處理過程中的氧化槽內氣液兩相的混合密度機理模型，從式(6)可以看出該混合密度存在多方面影響因素，基于新疆高寒高海拔地區(qū)金礦的工程實際特性，選取礦區(qū)環(huán)境溫度T、礦區(qū)大氣壓強P以及礦漿質量濃度Cl作為目標影響因素進行分析研究，從而得到礦區(qū)海拔較高，壓強較低以及晝夜溫差大等特性對氧化槽內氣液混合相密度的影響。

3 變量參數(shù)對氣液兩相混合密度的影響

利用工程實際中的參數(shù)，可以擬合出生物氧化預處理工藝的礦區(qū)環(huán)境溫度、礦區(qū)大氣壓強以及礦漿質量濃度等影響參數(shù)變化對氧化槽內氣液兩相的混合密度的影響曲線。以新疆某高寒地區(qū)的含砷難處理金礦為工程背景，該礦區(qū)緯度較高，地形屬于干旱半干旱北溫帶戈壁緩丘，地勢平坦，海拔高度約1300 m；該礦區(qū)屬于典型的大陸性氣候，晝夜溫差較大，年平均氣溫在8.5～10.5℃，最高溫度為42℃，而最低溫度為-35℃，其生物氧化預處理工藝參數(shù)如表1所示。

表1 生物氧化預處理參數(shù)Tab.1 Biological oxidation pretreatment parameters

生物氧化預處理工藝氧化槽內安裝有冷水盤管與熱水盤管，用來調節(jié)氧化槽內背景溫度，溫度低于理想值時開啟熱水盤管加熱，而當溫度高于理想值時開啟冷水盤管降溫，從而使得反應在最適宜溫度下進行，得到較高的生物氧化效率。雖然氧化槽內背景溫度可以基本維持恒定，但是由于氧化槽體積較大，均在室外放置，因此環(huán)境溫度對氧化槽內氣液兩相的密度有較大的影響，而該礦區(qū)晝夜溫差較大，使得氧化槽內氣液混合相的氣相含率受進氣溫度影響，同時也影響了氣液兩相的混合密度。根據已知參數(shù)以及密度機理模型(式(2)～(5))擬合出氧化槽內礦漿濃度分別為15%、18%及20%時對應氣液兩相的混合密度隨礦區(qū)環(huán)境溫度T變化的曲線，如圖2所示。

圖2 密度隨礦區(qū)環(huán)境溫度的變化曲線Fig.2 The curve between density and environment temperature

由圖2可知，生物氧化預處理過程中，氧化槽內礦漿質量濃度為15%時對應的密度大，而20%時對應的密度卻??；且氣液兩相的密度隨溫度的升高而減小。

中國是擁有生物氧化提金礦廠最多的國家，各廠所在地理位置對應的海拔高度不同，對應的各礦區(qū)大氣壓強也不同(海拔越高，壓強越低)。因此，生物氧化預處理氧化槽內氣液兩相的混合密度也會發(fā)生改變。假設廠區(qū)位置分布于海平面(0 m)到海拔高度(1300 m)之間，其對應的大氣壓強則從101325 Pa減小到84854 Pa，根據已知參數(shù)以及密度機理模型(式(2)～(5))，可得到密度隨不同海拔高度的礦區(qū)大氣壓強的變化曲線，如圖3所示。

圖3 密度隨礦區(qū)大氣壓強的變化曲線Fig.3 The curve between density and atmospheric pressure

由圖3可知，氧化槽內氣液兩相的密度隨氣壓的增大而增大；當氣壓變化18 kPa，即海拔變化1300 m，對應的密度變化了僅僅0.001 kg/m3，可見，礦區(qū)大氣壓強對氧化槽內氣液混合相的密度影響非常小，因此，礦區(qū)所在位置的海拔因素影響變得很小，在分析時考慮權重也可以變小。

工程實際中對難處理金精礦的生物氧化預處理，需要將礦漿質量濃度Cl調配在理論最優(yōu)范圍內(即15%～20%)[9]，可見，當?shù)V漿質量濃度不同時對應氧化槽內氣液兩相的混合密度也會有所不同。根據已知參數(shù)以及密度機理模型(式(2)～(5))擬合出礦區(qū)環(huán)境溫度分別為最低溫238 K和最高溫315 K時氧化槽內氣液兩相的混合密度受礦漿質量濃度變化的影響曲線，如圖4所示。

圖4 密度ρ隨礦漿質量濃度Cl的變化曲線Fig.4 The curve between density and mass concentration

由圖4可知，礦區(qū)環(huán)境溫度較低時混合相密度較溫度較高時的大，且氧化槽內氣液兩相的密度隨礦漿濃度的增大而減??；當?shù)V漿濃度變化 5%，對應的密度變化400 kg/m3，由變化幅值可見，礦漿質量濃度對氧化槽槽內氣液混合相的密度影響較大，必須要作為其影響因素在工程實際中加以考慮。

4 對氧化槽設備的改進

基于氧化槽內氣液混合相的密度影響因素的分析，為減小晝夜溫差極大的礦區(qū)環(huán)境溫度對槽內氣液混合相密度的影響，可以在空氣進入氧化槽之前對其進行預熱保溫處理，使其保持工藝要求的最適宜溫度(一般為42 )℃。在進氣管外嵌套一個恒溫水管，并加上保溫層進行保溫，可達到預熱進氣溫度的目的；為使得氧化槽內進入的礦漿質量濃度盡量保持在工藝參數(shù)最優(yōu)值(一般為 18%)參加反應，可以將氧化槽設備原本的頂部進漿改為底部進漿。如此，最優(yōu)濃度的礦漿進入氧化槽后便可及時在富氧區(qū)活躍細菌的催化作用下開始反應，在攪拌作用下反應更加充分快速，則同時提高了細菌生物氧化反應的速度和效率，也縮短了反應時間，可提高工藝經濟性。改進后的氧化槽設備如圖5所示(注：只更改氧化槽設備，并不改變生物氧化預處理工藝流程)。與圖1進行比較，改進部分為給進氣管外加了恒溫熱水管，將進漿管由氧化槽頂部更改為底部。

圖5 改進后的氧化槽設備Fig.5 The improved oxidation tank equipment

對氧化槽設備進行改進后，可降低礦區(qū)環(huán)境溫度以及礦漿質量濃度變化對槽內氣液混合相密度的影響，進而使該混合密度盡量作為一個穩(wěn)定的物性參數(shù)在生物氧化反應中體現(xiàn)，故混合密度的變化對細菌生物氧化反應的干擾系數(shù)會越小，對反應速度與效率的影響也就越小。

5 結語

本文通過對新疆高寒高海拔地區(qū)生物氧化提金預處理工藝過程的研究，得到了生物氧化預處理過程中氧化槽內氣液混合相的密度機理數(shù)學模型。結合氣液兩相流的混合密度公式分析獲得了礦區(qū)環(huán)境溫度T、礦區(qū)大氣壓強P以及礦漿質量濃度Cl對氧化槽內氣液兩相的混合密度ρ變化的影響關系，并據此對氧化槽設備進行了改進。這些工作為氧化槽內密度場的建立以及數(shù)值模擬研究提供了理論支持，在后續(xù)的研究中將以該密度模型為基礎，分析生物氧化預處理過程中氧化還原電位(ORP)與氧化槽內氣液混合相的密度之間的關系，從而進一步揭示高寒高海拔地區(qū)生物氧化還原反應的內在規(guī)律和特點，促進生物氧化提金技術的發(fā)展和完善。

[1] 朱長亮, 楊洪英, 湯興光, 等. 含砷難處理金礦的細菌氧化預處理研究現(xiàn)狀[J]. 貴金屬, 2010, 31(1): 48-52. ZHU C L, YANG H Y, TANG X G, et al. Current status of studies on bacterial pre-oxidation and leaching of refractory gold ores with As[J]. Precious metals, 2010, 31(1): 48-52.

[2] 馬紅周, 燕超, 王耀寧, 等. 氫氧化鈉浸出含砷金礦中砷[J]. 貴金屬, 2015, 36(1): 14-16. MA H Z, YAN C, WANG Y N. Leaching the arsenic in the gold ore by sodium hydroxide[J]. Precious metals, 2015, 36(1): 14-16.

[3] 南新元, 陳飛, 孔軍. 高寒地區(qū)某金礦生物氧化預處理過程溫度控制研究[J]. 貴金屬, 2014, 35(2): 38-42. NAN X Y, CHEN F, KONG J. The research on temperature control in biological oxidation pretreatment process of a gold mine in alpine region[J]. Precious metals, 2014, 35(2): 38-42.

[4] ARRASCUE M E L. Biooxidation of arsenopyrite concentrate using BIOX process: Industrial experience in the Tamboraque, Peru[J]. Hydrometallurgy, 2006, 83(1): 90-96.

[5] 伍喜慶, 程征, 楊平偉. 礦漿pH對表面磁化和菱鐵礦-赤鐵礦協(xié)同磁化的效應[J]. 中南大學學報: 自然科學版, 2014, 45(3): 678-683. WU X Q, CHEN Z, YANG P W. Effects of pulp pH on surface magnetization and siderite hematite synergic magnetization[J]. Journal of Central South University: Science and technology, 2014, 45(3): 678-683.

[6] 楊松榮. 含砷難處理金礦石生物氧化提金基礎與工程化研究[D]. 長沙: 中南大學, 2004. YANG S R. Fundamental and industrialization investigation on biooxidation of arsenic-bearing refractory gold ore[D]. Changsha: Central South University, 2004.

[7] 高丙朋, 南新元. 生物氧化預處理工藝控制系統(tǒng)研究與應用[J]. 黃金, 2011, 32(1): 50-52. GAO B P, NAN X Y. Research and application of control system in biological oxidation pretreatment artwork process[J]. Gold, 2011, 32(1): 50-52.

[8] LI Q, LI D, QIAN F. Pre-oxidation of high-sulfur and high-arsenic refractory gold concentrate by ozone and ferric ions in acidic media[J]. Hydrometallurgy, 2009, 97(1): 61-66.

[9] 高丙朋, 南新元, 魏霞. 基于迭代LS-SVM生物氧化提金預處理工藝參數(shù)優(yōu)化算法的研究[J]. 貴金屬, 2012, 33(2): 40-43. GAO B P, NAN X Y, WEI X. Study on optimization algorithm of process parameters for gold extraction pretreatment based on iterative LS-SVM[J]. Precious metals, 2012, 33(2): 40-43.

[10] 鄭曄. 高寒環(huán)境下生物氧化提金工藝研究與應用[J].中國礦業(yè), 2013, 22(6): 70-74. ZHENG Y. Research and application of biological oxidation technology on gold extraction in gold area[J]. China mining magazine, 2013, 22(6): 70-74.

[11] MISHIMA K. Visualization and measurement of gasliquid metal two-phase flow with large density difference using thermal neutrons as microscopic probes[J]. Nuclear instruments and methods in physics research, 2011, 424: 229-234.

[12] DORAO C A. Effect of inlet pressure and temperature on density wave oscillations in a horizontal channel[J]. Chemical engineering science, 2015, 134: 767-773.

Density Mechanism Model on Gas-liquid Mixing Phase in the Oxidation Tank of Biological Oxidation Pretreatment

JIA Aidi, GAO Bingpeng*, XU Yuanyuan
(Eelectrical Engineer School, Xinjiang University, Urumqi 830047, China )

In the research of biological oxidation gold pretreatment to recovere gold from the cold and high areas in Xinjiang, a gas-liquid mixing phase density model by the mechanism modeling method was established. By the calculation of the gas-liquid mixing phase density in the oxidation tank, the fitting curve that the influence of environment temperature, atmospheric pressure and mass concentration to mixture density was resulted, and the oxidation tank equipment improvement was developed. The results show that there are direct proportion relationship between mixture density with environment temperature and mass concentration. But there are inverse proportion relationship between mixture density and atmospheric pressure. This modified oxidation tank equipment can improve the reaction speed and efficiency of biological oxidation reaction.

biological oxidation; gold mine; gas-liquid mixing phase; mixture density; density mechanism model

TF831

：A

：1004-0676(2016)03-0046-05

2015-11-15

新疆自然科學基金項目(2012211A004)。

賈愛迪，女，碩士研究生，研究方向：軟測量建模及工業(yè)控制優(yōu)化。E-mail: 654615992@qq.com

*通訊作者：高丙朋，男，博士，副教授，研究方向：智能控制與系統(tǒng)開發(fā)及PLC應用。E-mail: 155258792@qq.com