論文題目內熱整合型蒸餾塔之動態模擬.DOC

論文題目：內熱整合型蒸餾塔之動態模擬 Dynamic Simulation for Internally Heat-Integrated Distillation Column 研究生：史詔仰 (Shih, Chao-Yang) 指導教授：黃琦聰 博士(Dr. Huang, Chi-Tsung) 文獻上有關內部熱整合型蒸餾塔的探討大部份均著重於其之模式化、分析與設計；有關其控制方面，日本學者指出精餾段之壓力控制十分不穩定。本研究動態模擬之架構乃基於DuPont公司Franks (1972) 之模組化(modular)結構，並以FORTRAN語言撰寫，自行發展出內熱整合型蒸餾塔之動態模擬軟體。熱力學模式則採用SRK與Modified Raoult’s law兩種；而活性係數則是利用NRTL來計算。另一方面，本研究亦提出一種HIDiC之控制策略，並考慮精餾段各板壓力之動態來進行HIDiC之模擬。於低壓下，採用Modified Raoult’s law，針對苯、甲苯、二甲苯之三成份系統，利用此軟體來進行動態模擬。於中壓下之動態模擬則採用SRK之熱力學模式，考慮汽體滯留量之變化，並以丙烷與丙烯之二成份系統為例題。由電腦模擬的結果發現，本研究所提出之控制策略，對於進料組成之變化及進料流量的變化干擾，均能有相當不錯的控制效果。 簡介 HIDiC的主要原理主要是利用蒸餾塔的汽提段(stripping section)與精餾(rectifying section)段進行熱轉移，將精餾段所移出的熱，來供給汽提段加熱之用，如圖1所示為傳統之蒸餾塔，圖2 所示為HIDiC。理論上在傳統的蒸餾操作下，精餾段的溫度比汽提段低，並無法傳熱給汽提段，因此在HIDiC的設計上就必須提高精餾段的操作壓力，使整段的溫度分布高於汽提段，且其間必須裝設有蒸汽再壓縮機將汽提段的蒸汽傳送到精餾段，以完成整個蒸餾操作。利用McCabe-Thiele圖解方式，也可以簡單說明HIDiC的能源需求低於傳統的蒸餾塔。如圖3為傳統蒸餾塔之平衡曲線與操作線之關係；然而，由於理想的HIDiC並沒有外部回流，其蒸餾操作曲線是呈弧形狀態，如圖4所示。換言之，HIDiC的內部回流是隨著熱傳量的大小而變化；其變化的趨勢從塔頂到塔底來觀察，其在精餾段是由小變大、汽提段則是由大變小，且呈連續變化。另一方面，這樣的操作亦可減少反向質傳的現象出現，而且因沒有物質的直接接觸，也避免了反混(remixing)的現象。 圖1 傳統型蒸餾塔 圖2 HIDiC 圖 3 McCabe -Thiele圖 傳統型蒸餾塔 圖 4 McCabe -Thiele圖 內熱整合型蒸餾塔 低壓之平衡板計算 在進行低壓(1~5atm)之動態模擬時，本研究將塔內板視為理想板以平衡板(equilibrium stage)的概念來處理，其平衡板模型如圖5所示。為簡化模擬的複雜性及考慮實際的情況而作了下列的假設與限制： (1) 氣相與液相的流動均是當氣液組成達到平衡後才離開該板。 (2) 假設於每板中的溫度、壓力與組成分佈皆為常數，也就是板上每點的條件皆相同。 (3) 每板中汽體滯留量（vapor holdup）極少于以忽略不計，即假設每板中汽體流動時無夾帶液體。 (4) 每板的板效率皆視為100%；氣相與液相均處於熱力學平衡，即VLE (vapor-liquid-equilibrium)。 (5) 每板的液體之停留與流動的關係均以Francis weir formula表示之。 (6) 板上的滯留液體不會有分相(phase split)的行為。 而在其模式化上，則用質能平衡MESH方程式來描述。 低壓之熱力學模式 處理VLE計算時，則使用Modified 圖5 低壓平衡板簡圖 Raoult’s Law，即 其中為i成份的液相活性係數；P與Pisat分別為總壓力和i成份的飽和蒸汽壓，其單位是bar；Ki則是i成份的相平衡常數。可為NRTL model ,UNIQUAC model或Wilson model。 中高壓之平衡板計算 若蒸餾塔操作壓力大於5-10 atm時，為中高壓蒸餾，此時汽體密度增加，在平衡板的考量方面，除了必須考慮每板汽體滯留量影響之外，其餘假設與低壓部分相同，其中高壓平衡板模擬簡圖如圖6。而在模式化上亦是利用質能平衡MESH方程式來描述。 中壓之熱力學模式 對於真實氣體混合物於中高壓時在氣液平衡(VLE)的計算上，則使用狀態方程式(equation of state) SRK (Soave-Redlich-Kwong; Soave, 1971)來處理。其為