水-丙二醇單甲醚-1，2-丙二醇三元物系分離模擬

蔡 乾，曹貴平，呂 慧

（華東理工大學 化工學院 聯合化學反應工程研究所 化學工程聯合國家重點實驗室，上海 200237）

在過氧化氫氧化丙烯合成環氧丙烷的過程中，產物脫除甲醇等組分后的廢水中含有低濃度的丙二醇單甲醚（MME）和1，2-丙二醇（PG）［1-2］。MME是一種低毒性的優良溶劑，被廣泛用于工業和消費品中［1-2］，同時也是燃料抗凍劑、洗滌劑和萃取劑等的重要原料［3-4］。PG是生產化妝品和不飽和聚氨酯的原料［5-6］，在除冰溶劑、個人護理品以及食品的防腐劑中發揮著不容忽視的作用［7-8］，因此，對MME和PG的回收利用具有重要的價值。體系中的水（WT）與MME會形成最低恒沸點共沸物［3，9］，如采用精餾方法分離，存在流程長、能耗高的缺點。Zhao等［10］通過萃取與精餾的耦合方式分離了WT-MME二元體系，但MME-PG-WT三元體系含有少量的PG，萃取時PG在萃取相和萃余相存在不同的分布，影響后續分離效率，因此萃取劑的選擇十分重要。

本工作通過三元相圖研究了苯（BE）、甲苯（MB）、乙苯（EB）、二甲苯（DB）和異丙苯（IB）萃取劑對MME和PG的萃取效果。利用單級萃取模擬計算，篩選了適合WT-MME-PG體系的最佳萃取劑和最佳萃取溫度。通過多級萃取模擬，分析了萃取的最佳操作條件，為后續的精餾過程提供基礎數據，對工業回收化學品及處理有機廢水提供指導。

1 熱力學方程

合適的熱力學方程是成功模擬的基礎。前人在研究WT-MME體系的分離時發現NRTL熱力學方程可以準確地預測體系中各組分的性質以及二元交互參數［9-12］，胡松等［13］采用NRTL熱力學方程對環氧丙烷-水-甲醇-PG四元體系進行了模擬計算。本工作在前人研究的基礎上，采用NRTL熱力學方程對本體系進行模擬。在Aspen Plus計算時，使用了如式（1）所示的NRTL方程式［14-15］，該方程定義了液相組分中組分i的活度系數（γi）：

式中，x為液相中各組分的摩爾分數；G為吉布斯自由能，kJ/mol；τ為關于組分i和j的二元相互作用的參數，與溫度有關；G和τ可通過式（2）～（4）計算：

式中，Δgi j表示組分分子對之間的相互作用能；R為氣體常數，J/（mol·K）；T為絕對溫度，K；a，b，c為非隨機參數，與分子間相互作用能量有關。

NRTL熱力學方程下體系中各組分的參數見表1～4，其中S為萃取劑。

表1 NRTL熱力學方程下WT-MME-PG的參數Table 1 WT-MME-PG parameters under NRTL thermodynamic equation

表2 NRTL熱力學方程下WT-S的參數Table 2 WT-S parameters under NRTL thermodynamic equation

表3 NRTL熱力學方程下MME-S的參數Table 3 MME-S parameters under NRTL thermodynamic equation

表4 NRTL熱力學方程下PG-S的參數Table 4 PG-S parameters under NRTL thermodynamic equation

2 計算基礎

2.1 指標

萃取劑的選擇直接關系到萃取操作的可行性、經濟性及后續分離操作的難易程度。理想的萃取劑不僅能對目標組分有優異的萃取能力，而且可減小后續的分離難度、降低分離成本。分配系數（k）和萃取率（η）被廣泛用作衡量萃取劑的性能優劣的指標［15］。ki=yi/xi，其中，i表示組分；yi和xi分別表示組分i在萃取相與萃余相中的質量分數。ki越大，表明組分i在萃取相中的含量越多，即萃取效果越明顯。萃取率是量化萃取劑效果的重要指標，經萃取后，組分i的萃取率為：ηi=mEyi/（mEyi+mRxi），其中，mE為組分i在萃取相中的質量流量，mR為組分i的總質量流量。萃取劑在萃余相中的殘余量（RS）也是評價萃取劑優劣的一個重要指標，RS越小，后續對萃余相中的關鍵組分的提純就越方便，同時也可以減少萃取劑的損耗。

2.2 模擬基礎

針對WT-MME-PG體系，期望萃取劑對MME的萃取效果較好，而對PG和WT不具萃取效果，理想情況是經萃取后，全部的MME和萃取劑作為萃取相流出，全部的PG和WT作為萃余相流出，這樣可以很大程度上簡化后續的分離過程。設定萃取壓力為常壓，萃取溫度為25 ℃，原料進料量為1 000 kg/h，原料含有96.88%（w）WT、1.69%（w）MME和1.43%（w）PG，萃取劑與原料質量進料比為1。

2.3 三元相圖

三元相圖被廣泛用作選擇合適萃取劑的分析和解釋圖形工具［16-17］。典型的單級萃取三元相圖如圖1所示，三角形的3個頂點表示3個純組分。溶解度曲線將三角形內部分為兩個區域，陰影部分為兩相區域，陰影部分以外區域為單相區，萃取操作需要在兩相區進行，兩相區域越大，表示萃取劑和溶劑互溶度越小，越有利于萃取。

圖1 單級萃取圖解法Fig.1 Graphic method of single stage extraction.A：pure solvent A；B：solute B；S：pure extractant；F：raw mixture of A and B；E：extract phase；R：raffinate phase；E’：extract phase after extractant removal completely；R’：raffinate phase after extractant removal completely.

萃取操作時，含有組分A和B的原料液（F點）和萃取劑（S點）的混合組成位于兩相區內（M點），萃取劑與原料液充分混合，靜置分層為萃取相（E點）和萃余相（R點），各液相的質量關系可由杠桿規則描述。分別連接S點和E點、S點和R點并延長，交AB邊于E’點和R’點，E’點和R’點分別表示當萃取劑完全脫除后萃取液和萃余液組成?？梢钥闯?，相比F點，E’點中的A組分含量增加，B組分含量減少；R’點中的A組分含量減少，B組分含量有一定程度的增加。因此，經單級萃取后，萃取劑對原料液中的溶質A具有萃取作用，實際操作時，可通過改變萃取條件滿足設定要求。

3 結果與討論

3.1 三元相圖的繪制

采用Aspen Plus軟件的Ternary Diag功能繪制三元相圖。WT-MME-BE的三元相圖見圖2。對圖2下部紅色方框區域進行放大處理，可清晰看出萃取相和萃余相中各組分的組成，其余四種溶劑作萃取劑時的三元相圖放大圖見圖3。WT-PG-S三元相圖放大圖見圖4。WT-MME-PG和MMEPG-S體系為均相，三元相圖不再繪出。

圖2 WT-MME-BE三元相圖Fig.2 Ternary phase diagram of WT-MME-BE.

圖3 WT-MME-S三元相圖放大圖Fig.3 Enlarged ternary phase diagram of WT-MME-S.

圖4 WT-PG-S三元相圖放大Fig.4 Enlarged ternary phase diagram of WT-PG-S.

3.2 三元相圖分析

根據杠桿原理可求得各體系中代表萃取相和萃余相的E點和R點的組成，進而求得各組分的分配系數，WT-MME-S和WT-PG-S體系中各組分的分配系數見表5。

由表5可看出，在WT-MME-S體系中，MME的分配系數在BE和IB作為萃取劑時分別為0.417和0.154；在MB，EB，DB作萃取劑時均為1.125，表明BE和IB對MME的萃取效果較差，不予考慮。MB作萃取劑時WT的分配系數為0.001，小于EB和DB作萃取劑時的0.008和0.003，表明MB作萃取劑時萃取相的WT含量相對較少。EB，MB，DB作萃取劑時它們的分配系數分別為491.500，991.000，988.000，說明EB作萃取劑時在萃余相中有較多的殘留。在WT-PG-S體系中，采用EB、DB和MB作萃取劑時，PG的分配系數分別為0.250，0.167，0；WT的分配系數在EB，MB，DB三種萃取劑下接近于WT-MME-S體系。綜合比較可知，5種萃取劑中BE和IB對MME萃取效果較差，MB，EB，DB對MME的萃取效果接近，且萃取相中WT和PG含量均較少。選擇MB，EB，DB作為備選萃取劑，采用Aspen Plus軟件分離模擬得到最優萃取劑，并確定最佳操作參數。

表5 WT-MME-S和WT-PG-S體系中各組分分配系數Table 5 Partition coefficients(k) of components in WT-MME-S and WT-PG-S systems

3.3 單級萃取

考察不同萃取劑條件下的單級萃取。選擇Aspen Plus軟件中的Decanter模塊進行單級萃取模擬，設定溶劑比（萃取劑與原料液質量流量比）為1?？疾觳煌腿┫聹囟葘Ω鹘M分萃取率及RS的影響，結果見圖5。由圖5a可知，MME的萃取率隨溫度的升高而增大，當萃取劑分別為MB，EB，DB時，MME的萃取率在5 ℃時分別為0.425，0.401，0.403，升至30 ℃時則分別為0.459，0.441，0.440，因此，與EB和DB相比，MB對MME的萃取效果更為顯著。由圖5b可看出，PG的萃取率隨溫度的升高而增大，當萃取劑為MB時，PG的萃取率在5 ℃時為0，隨著溫度升至30 ℃，PG的萃取率達到0.001；PG的萃取率在EB和DB作萃取劑時一致，在5 ℃時為0.141，升至30 ℃時為0.150。從圖5c可知，WT的萃取率隨溫度的升高而增大，在5 ℃時，WT的萃取率在MB，EB，DB作萃取劑時分別為0.2×10-4，2.2×10-4，2.8×10-4，在30 ℃時 則 分 別 為0.5×10-4，7.1×10-4，6.8×10-4，在MB為萃取劑時，萃取相中的WT含量低于EB和DB作萃取劑時的含量。從圖5d可看出，采用MB、EB和DB萃取劑時，隨溫度的升高，RS增大，其中，MB作萃取劑時的RS略高于EB和DB作萃取劑時的RS。

圖5 不同萃取劑下溫度對各組分萃取率及RS的影響Fig.5 Effect of temperature on extraction rate(η) of component and extractant residue in raffinate phase(RS) with different extractants.

綜合而言，MB對MME的萃取效果優于EB和DB，而對PG和WT的萃取效果不如EB和DB，因此滿足MME作為萃取相流出、PG和WT作為萃余相流出的目標，盡管MB作萃取劑時的RS略高，但由于殘余量較少，此影響可以忽略，故對于WT-MME-PG體系，MB為最佳萃取劑。

升高溫度可增加MME的萃取率，但會導致PG和WT萃取率的增大和RS的增加，而實際目標希望RS越小越好、萃取相中的PG和WT含量越低越好，且由于MME的萃取率可通過增大溶劑比和萃取級數來提高，因此實際萃取溫度越低越好，故取萃取溫度為5 ℃。

3.4 多級萃取

為了確定最佳溶劑比和萃取級數，進一步進行了多級萃取模擬研究。選用Aspen Plus軟件中的Extract模塊，MB作萃取劑，萃取溫度5 ℃，考察溶劑比和萃取級數對各組分的萃取率及RS的影響，如圖6所示。由圖6a可知，MME的萃取率隨溶劑比和萃取級數的增大而增大。當溶劑比取1.8時，在萃取級數為19時MME的萃取率達到1。由圖6b～c可看出，PG和WT的萃取率均隨著溶劑比的增加而增大，隨萃取級數的增加先增大然后基本保持不變。從圖6d可看出，在不同溶劑比下，RS隨萃取級數的增大而減小，在萃取級數大于17后趨于穩定。在考慮設備成本和溶劑成本需要減小溶劑比和萃取級數的前提下，為了獲得最大MME萃取率和最小RS，因此，取溶劑比和萃取級數分別為1.8和19。

圖6 不同溶劑比下萃取級數對各組分萃取率以及RS的影響Fig.6 Effects of extraction stages(NT) on η of component and RS under different solvent ratios(δ).Conditions：MB extractant，5 ℃.δ：● 1.6；■ 1.8；▲ 2.0

以MB為萃取劑，在萃取溫度5 ℃、溶劑比1.8、萃取級數19的條件下進行模擬計算，帶詳細物流信息的萃取效果見圖7。

圖7 最優條件下的萃取效果Fig.7 Extraction effect diagram under optimal conditions.Conditions：MB extractant，5 ℃，δ=1.8，NT=19. FT：total flow rate.

從圖7可看出，萃取相為MME和MB，含微量的PG和WT，萃余相為PG和WT，含微量的MME和MB。

4 結論

1）采用Aspen Plus軟件模擬WT-MME-PG的初步萃取分離，選擇5種萃取劑，通過繪制三元相圖分析不同萃取劑對各組分分配系數的影響，發現MB，EB，DB對MME有較好的萃取效果。

2）單級萃取下，比較了MB，EB，DB作萃取劑時MME和PG的萃取率及萃余相中的殘余量，實驗結果表明，MB為適宜的萃取劑，5 ℃為適宜的萃取溫度。

3）通過多級萃取模擬研究，發現MB作萃取劑，5 ℃下的最優操作條件為溶劑比1.8、萃取級數19。在該條件下，MME和MB作為萃取相流出，WT和PG作為萃余相流出，極大地簡化了后續分離工作。本工作對萃取劑的選擇、工業回收有機組分及處理有機廢水提供了一定的參考。