毛維東，孫 邃，馬 賽，鄒 湘，郭中權，呂順之，肖 艷

(1.中煤科工集團杭州研究院有限公司，浙江 杭州 311201；2.江蘇新宜中澳環(huán)境技術有限公司，江蘇 宜興 214205)

0 引 言

我國是世界上最大的煤炭生產國和消費國，煤炭資源分布“西多東少、北富南貧”，與水資源“東多西少、南富北貧”的分布特點相反。西部煤炭資源豐富，占據全國煤炭資源總量的70%，但西部地區(qū)處于干旱—半干旱帶，水資源匱乏，地表植被稀疏，生態(tài)環(huán)境脆弱[1-2]。

高礦化度礦井水是指含鹽量大于1 000 mg/L的礦井水[3]。西北地區(qū)的礦井水以高礦化度礦井水為主，占比達50%以上[4]。隨著煤炭生產開發(fā)布局的優(yōu)化，國家能源供應戰(zhàn)略西移[5]，高礦化度礦井水在全國煤礦礦井水的占比大幅度增加，成為制約西北煤礦發(fā)展最為突出的問題。

現有高礦化度礦井水除鹽工藝主要包括化學藥劑法、離子交換法、蒸餾法、電滲析法和反滲透法等[3, 6]，其中反滲透法在高礦化度礦井水處理中的應用越來越廣[7-8]。但各大膜廠家提供的工藝設計軟件，包括ROSA、WAVE、WINFLOWS等，均是針對單級反滲透系統的苦咸水、海水處理需求，缺少針對多級膜系統的過程設計方法。在面對礦井水水質的多變性以及無機鹽成分的復雜性時，現有反滲透系統在實際運行過程中，雖然出水水質基本可以得到保證，但膜元件易堵塞、易結垢，更換頻繁，導致運行成本較高[9]。

為降低工程投入和運行費用，提高礦井水回收率，本文以某煤礦高礦化度礦井水零排放工程的多級反滲透系統為例，結合CFD模擬技術，對原有多級反滲透系統進行了工藝優(yōu)化。

1 工程概況

鄂爾多斯某煤礦位于內蒙古自治區(qū)鄂爾多斯市伊金霍洛旗境內，屬于缺水地區(qū)，水資源的缺乏直接影響了煤礦的正常經營。根據環(huán)保部門要求，礦井水需全部回用于生產生活，實現零排放。該煤礦高礦化度礦井水零排放工程，設計處理能力為600 m3/h，工藝流程分為凈化處理、深度處理、濃縮處理和蒸發(fā)結晶4個單元。其中深度處理和濃縮處理單元采用了多級反滲透系統進行除鹽處理，該系統設計處理能力為600 m3/h，工藝流程圖如圖1所示。

圖1 礦井水凈化處理、深度處理、濃縮處理工藝流程[8]Fig.1 The process flow of purification, advanced and concentration treatment for coal mine water[8]

整套脫鹽系統中包含的三個反滲透系統在圖1中以灰色框表示：其中一級膜濃縮的設計水回收率為75%；二級膜濃縮的設計水回收率為65%；三級膜濃縮的設計水回收率為60%。這其中一級、二級膜濃縮的膜組件均為海德能PROC10卷式膜元件。三級膜濃縮則采用德國ROCHEM水處理公司的DTRO膜片。

實際運行過程中，多級反滲透系統出現了進水指標超出原設計范圍的問題，如進水濁度、鈣、鎂、TOC含量等，造成膜污染與結垢，影響系統運行效果和膜原件壽命[10-11]。

2 CFD模擬礦井水多級反滲透處理系統

2.1 卷式RO膜元件模擬

反滲透膜過濾脫鹽過程的流體力學模擬涉及到溶質(無機鹽)和溶劑(水)分別在進水流道和產水流道的物料傳質模擬。水流流動的控制方程為不可壓縮流體的質量守恒和動量守恒方程(納維—斯托克斯方程)。

無機鹽在整個流場內的運動通過物料守恒方程表示：

其中，Cj表示對應無機鹽j的濃度，ui和xi分別代表沿著第i個坐標方向的流速和距離，D為對應無機鹽的擴散系數，Ssj為溶質的質量源匯項。

卷式RO膜元件采用陶氏BW30-400，由10片長2.05 m、寬0.9 m的膜頁卷制而成。濃水通道厚度為濃水網厚度，產水通道厚度由膜袋之間的清水網厚度決定[12]。運行中，進水從一側濃水通道進入，濃水從另一側濃水通道出口排出，而產水從中心集水管中排出[13]。

單個膜袋與相鄰膜袋可以認為是并聯運行，有一定相互作用(共享同一濃水通道)，但獨立性更顯著，所以在模擬過程中采用幾何建模軟件ANSYS DesignModeler對單膜片、單濃水通道進行建模，如圖2所示，采用ANSYS CFX 19.1自配的Meshing模塊對模擬幾何結構劃分網格，其中，邊界層網格設置5層與膜片平行的六面體網格，第一層厚度設置為0.02 mm，逐層增長率為1.5倍，來精確模擬邊界層無機鹽成分的濃差極化現象。其余區(qū)域采用平均網格尺寸為0.05 mm的四面體網格建模，經網格獨立性測試后，總網格數近1 300萬。

圖2 卷式反滲透單進水流道模型Fig.2 The single inlet model for a spiral woundRO element

膜元件設計進水8.6 t/h，對應進水流道內平均雷諾數Re為2 600，水流處于層流紊流過渡區(qū)，適用不完全紊流的k-ω紊流模型模擬。濃水網表面設置為無滑移邊界條件。采用高精度中心差分法解穩(wěn)態(tài)水流控制方程，所有模擬流速、壓強收斂于10-4，無機鹽濃度收斂于10-6，其他采用默認值。

卷式RO元件內復合無機鹽濃度沿軸向分布如圖3所示，在卷式RO膜元件運行過程中，復合無機鹽濃度沿軸向逐漸增高，在濃水通道靠近膜面附近，出現明顯的濃差極化現象。

2.2 DTRO膜元件模擬

對DTRO的導流盤進行3D掃描，得到導流盤的尺寸及突起結構的排列分布位置，導流盤半徑為96 mm，突起形狀為1.5 mm的半球體，如圖4所示。

圖4 導流盤結構組成及剖面圖顯示Fig.4 Composition of DTRO and display of cutaway view

DTRO導流盤具有對稱性，對幾何結構進行對稱簡化后再進行網格劃分，網格劃分步驟與卷式反滲透膜方法一致，經過獨立性驗證后網格總數1 050萬。進口設計為流體速度進口，進口速度設置為0.1 m/s；出口設計為壓力出口，壓力出口設置為80 bar。膜片過膜通量與無機鹽通量設置方法與上述研究卷式反滲透膜設置一致。

圖5 DTRO流速分布云圖Fig.5 Distribution of velocity in the partial DTRO element

經過CFD模擬，得到DTRO流速分布云圖(如圖5所示)和DTRO鹽度分布云圖(如圖6所示)。通過DTRO內部的流速分布云圖及鹽濃度分布云圖，可以看出，由于導流盤上存在圓柱凸起，在其后方產生“卡門渦街”，增強了流體的湍流動能，其自清洗效果優(yōu)異，減少了膜污染，降低了濃差極化現象。

圖6 DTRO鹽度分布云圖Fig.6 Distribution of salinity in the partial DTRO element

通過DTRO整體內的鹽度分布云圖(圖7)，我們可以看出，高鹽度主要出現在膜片中心附近，分布較為均勻，無明顯的濃差極化現象產生；從污染物的分布云圖(圖8)觀測到，導流盤突起附近出現了明顯的污染物堆積，與實際的污染膜情況基本一致。

圖7 整體鹽度分布云圖Fig.7 Distribution of salinity in the whole DTRO element

3 工藝優(yōu)化

鄂爾多斯某煤礦高礦化度礦井水零排放工程中的三級膜濃縮單元采用了多級耦合膜濃縮系統進行除鹽處理，該系統設計處理能力為600 m3/h。

一級膜濃縮采用苦咸水反滲透膜單元(BWRO)，該套脫鹽系統中的BWRO單元采用兩段式設計，其中第一段包含18個并聯的膜組件，第二段包含9個并聯的膜組件，每個膜組件里包含6個膜元件PROC10，具體如圖9所示。優(yōu)化后，一級膜濃縮單元依然采用兩段式設計，其中第一段膜組件變?yōu)?7個并聯，第二段膜組件依然為9個并聯，每個膜組件里包含6個膜元件PROC10。

圖9 一級膜濃縮組件排布Fig.9 The arrangement of the first stage RO membrane elements

二級膜濃縮處理采用苦咸水反滲透膜單元(BWRO)，三級膜濃縮采用碟管式反滲透膜單元(DTRO)。BWRO單元也采用兩段式設計，其中第一段包含6個并聯的膜組件，第二段包含3個并聯的膜組件，每個膜組件里包含6個膜元件PROC10，具體如圖10所示。DTRO單元包含多個互相疊加的圓環(huán)形膜片，如圖11所示。

優(yōu)化后，二級膜濃縮單元也依然采用兩段式設計，其中第一段膜組件變?yōu)?個并聯，第二段膜組件依然為3個并聯，每個膜組件里包含6個膜元件。DTRO單元不變。

4 技術性分析

通過計算流體力學分別模擬原始方案與優(yōu)化方案，從而得到膜內離子污染分布，進而分析工藝優(yōu)化前后進水條件的差異，如表1所示。由表1可以看出，優(yōu)化方案與原始方案相比，優(yōu)化方案進口壓力比原始方案進口壓力大約增加0.136 MPa，優(yōu)化方案進水流量比原始方案進水流量增加0.49 m3/h。

原始方案和優(yōu)化方案多級反滲透系統(第一段BWRO)中Ca2+整體的運行表現如圖12所示，從圖中可以看出膜元件內的Ca2+濃度沿程逐漸增高，比較優(yōu)化方案與原始方案，膜元件內Ca2+濃度沿程分布基本沒有差別并且Ca2+都在進口處發(fā)生沉積。

圖12 優(yōu)化方案與原始方案第一級第一段膜內Ca2+分布Fig.12 The Ca2+ distribution in the first section membrane element of the first stage for both optimized and original schemes

圖13 優(yōu)化方案與原始方案第一級第一段膜內分布Fig.13 The distribution in the first section membrane element of the first stage for both optimized and original schemes

原始方案和優(yōu)化方案多級反滲透系統(第一段BWRO)中Mg2+整體的運行表現如圖14所示，從圖中可以看出膜元件內的Mg2+濃度沿程逐漸增高，比較優(yōu)化方案與原始方案，膜元件內Mg2+濃度沿程分布基本沒有差別并且Mg2+在整個過程中不發(fā)生沉積。

圖14 優(yōu)化方案與原始方案第一級第一段膜內Mg2+分布Fig.14 The Mg2+ distribution in the first section membrane element of the first stage for both optimized and original schemes

原始方案和優(yōu)化方案多級反滲透系統(第一段BWRO)中Na+整體的運行表現如圖15所示，從圖中可以看出膜元件內的Na+濃度沿程逐漸增高，比較優(yōu)化方案與原始方案，膜元件內Na+濃度沿程分布基本沒有差別并且Na+在整個過程中不發(fā)生沉積。

圖15 優(yōu)化方案與原始方案第一級第一段膜內Na+分布Fig.15 The Na+ distribution in the first section membrane element of the first stage for both optimized and original schemes

圖16 優(yōu)化方案與原始方案第一級第一段膜內分布Fig.16 The distribution in the first section membrane element of the first stage for both optimized and original schemes

綜上所述，優(yōu)化方案多級反滲透系統膜內離子污染分布與原始方案基本保持一致，由此可以得出結論：優(yōu)化方案與原始方案之間投入阻垢劑用量基本保持一致。

5 經濟性分析

本次工藝優(yōu)化對一級膜濃縮單元、二級膜濃縮單元進行了優(yōu)化，三級膜濃縮單元維持原設計不變，對一級膜濃縮單元及二級膜濃縮單元進行經濟性分析對比。

如表2所示，某煤礦礦井水多級膜濃縮單元處理規(guī)模為600 m3/h，優(yōu)化后多級反滲透脫鹽系統的回收率由96.5%提高到97%，由于回收率的提高，最終濃水量由21 m3/h下降到18 m3/h；整個系統一級膜濃縮和二級膜濃縮的元件數量從216支下降到204支，減少了12支，約占原始系統的5%，按照平均膜價4 000元/支計算，由于優(yōu)化方案減少了12支膜，所以初始膜投入也相應減少48 000元。優(yōu)化后的電力消耗僅比原始方案多了0.02 kW·h/m3。

表2 優(yōu)化方案與原始方案主要參數對比

如表3所示，按照該煤礦總處理量為600 t/h，電價為0.86元/(kW·h)計算，優(yōu)化方案比原始方案每小時需要多投入10.32元，優(yōu)化方案的噸水電力消耗高0.017元/m3。在蒸發(fā)結晶處理費用上，按照成本79.09元/t濃水計算[8]，優(yōu)化方案比原始方案要節(jié)省237.27元/h的投入，優(yōu)化方案的噸水蒸發(fā)結晶成本低0.395元/m3。此時，與蒸發(fā)結晶處理節(jié)省的費用相比，電力所增加的投入可以忽略不計。

表3 優(yōu)化方案與原始方案主要投入對比

6 結 論

(1)CFD技術可以有效模擬卷式RO膜元件和DTRO膜元件的工作狀態(tài)，指導膜系統設計和運行維護。

(2)與原始方案相比，優(yōu)化方案的多級反滲透系統膜內離子污染分布基本同原來一致，所以阻垢劑用量也不變，藥劑費沒有增加。

(3)優(yōu)化后，BWRO單元和SWRO單元各減少6支膜，整個多級反滲透系統從216支下降到204支，減少了12支(約5%)，膜元件投資降低5%左右。

(4)與原始方案相比，礦井水回收率由96.5%增加到97%，噸水電力消耗升高了0.017元/m3，但噸水蒸發(fā)結晶費用降低了0.395元/m3，節(jié)省運行費用0.378元/m3，經濟效益明顯。