精餾塔塔板液位的密度修正算法

張煒

(中國天辰工程有限公司，天津 300400)

1 概 述

精餾塔是化工生產中的重要設備，利用不同介質的揮發度不同，在塔板上使得氣相和液相充分混合傳質，實現輕組分上升重組分下降，從而將不同揮發度的物質分離開來[1]。其中塔板液位是精餾塔的重要參數。如果塔板液位過高會造成泡沫夾帶增多從而降低分離質量，甚至造成淹塔使得上層塔板的液體直接進入了下層塔板，破壞了精餾塔的工作過程；如果塔板液位過低會使得下層塔板的氣相物質未與液相物質接觸而直接進入上層塔板，造成串塔同樣破壞了精餾塔的工作過程。

由于長期以來塔板液位的測量比較困難，塔板上的液體呈氣液兩相且劇烈“沸騰”的狀態，使傳統的差壓液位計無法應用，且塔板間空間狹小又難以安裝頂裝式液位計。一般只測量塔底液位和塔頂冷凝罐的液位，并根據進料和出料流量的變化和塔內溫度的變化估計塔板液位的高低。

2 密度修正算法

筆者在某項目中使用差壓液位變送器結合密度修正算法較好地解決了塔板液位測量的問題。該項目中采用泡罩塔結構，可以認為塔板上方的液體總質量是基本固定的。當液體質量過大時，液體會溢流到下層塔板；當液體質量過小時，冷凝液就會積存在塔板上。塔內上升氣流豐富，在塔板上形成均勻的氣液兩相混合物。在這2個假設的基礎上建立液位密度修正算法： 塔板上液體的總質量不變；氣液兩相充分混合，即其密度是均勻的[2]。

泡罩塔在靜態時有一個靜態液位，即沒有上升氣流時塔板的最高液位；在靜態液位的高度開一個取壓口，即高壓側取壓口；在預計的最高塔板液位上方再開一個取壓口，即低壓側取壓口。采用差壓液位變送器來測量這兩個取壓口的差壓[3]。

在塔正常工作時，如圖1所示，上升氣流使得塔板液位上升，塔板上方形成氣液兩相物質。

圖1 沸騰狀態的塔板液位示意

靜態時塔板上方液體的總質量m為

m=ρ1L1S

(1)

式中：ρ1——靜態液位時液體的密度，kg/m3；L1——塔板上方的靜態液位，即沒有上升氣流時塔板的最高液位，m；S——塔板面積，m2。

塔正常工作時液位上升到L，根據假設，m不變，則氣液兩相物質的密度ρ2：

ρ2=m/(LS)=ρ1L1/L

(2)

式中：L——塔正常工作時的液位，即氣液兩相物質頂部液位，m。

塔正常工作時，上下取壓口之間的差壓Δp：

Δp=ρ2g(L-L1)

(3)

式中：g——重力加速度，9.8m/s2。

將式(2)代入式(3)得到：

Δp=[ρ1L1g(L-L1)]/L

(4)

根據式(4)計算得塔正常工作時塔板液位高度L：

(5)

在式(5)中ρ1，L1是已知量，再根據差壓變送器測量得到的Δp即可計算出塔板液位L。

假設ρ1=1000kg/m3；L1=0.5m。代入式(5)得到：

L=2450/(4900-Δp)

(6)

以Δp為橫坐標，L為縱坐標得到圖2所示的差壓液位曲線。

圖2 差壓液位曲線

從圖2中可以看出，隨著Δp的增加，L的上升速度越來越快，當Δp趨近于4900Pa時，L趨向于上限，因此Δp或L應作大信號切除。

3 改進的密度修正算法

第2節密度修正算法的推導是建立在塔板上方液體總質量不變的基礎上的，但在實際生產過程中總質量會隨著上下塔板間差壓的變化而變化。當差壓變大時塔板能托舉的液體質量就增加，反之則減少。在理想狀態下，塔板上增加的液體產生的壓力與板間差壓應相等。

如圖3所示，在塔板上下氣相側增加2個取壓點，用來測量板間差壓。

圖3 差壓補償的液位計算方法

設塔板下方的氣相壓力減去塔板上方的氣相壓力為板間差壓ΔpB，塔板上增加的液體質量為Δm，則：

Δm=ΔpBS/g

(7)

結合式(1)得到塔工作時塔板上方的液體總質量：

m=ρ1L1S+ΔpBS/g

(8)

采用與第2節相同的推導過程可以得出：

(9)

式(9)中L1，ρ1是已知量，再通過差壓變送器測量得到Δp和ΔpB，然后可以計算出塔板液位L。

下面假設ρ1=1000kg/m3；L1=0.5m；ΔpB分別取500，1000，2000，3000Pa。以Δp為橫坐標，L為縱坐標得到圖4所示的差壓液位曲線。

圖4 差壓補償下的差壓液位曲線

從圖4可以看出隨著ΔpB的增加，相同Δp計算得到的L在下降。

4 結束語

筆者介紹的密度修正算法考慮了液位升高后氣液兩相混合物密度降低的情況，通過假設塔板上方液體總質量不變和氣液兩相混合物密度均勻得出了密度修正計算公式。該算法已在國內外多家工廠應用，取得了不錯的使用效果，并在密度修正算法的基礎上加入板間差壓補償，使得計算結果更加接近實際。

參考文獻：

[1]侯穎.三元物系共沸精餾塔建模與仿真[J].計算機仿真，2014(01)： 248-252.

[2]馬星若，阮文娟.差壓式水位計冷凝水密度修正對汽包水位測量的影響[J].上海電力，2003(06)： 531-533.

[3]王泉.精餾塔塔板阻力波動原因分析及預防措施[J].深冷技術，2012(06)： 9-11.

[4]何離慶.過程控制系統與裝置[M].重慶： 重慶大學出版社，1995.

[5]黃科林.精餾塔的最優設計[J].化工設計，1991(04)： 21-27.

[6]劉道德.化工設備的選擇與工藝設計[M].長沙： 中南工業大學出版社，1992.

[7]吳兆亮，邸進申.板式塔塔板主要結構尺寸的新設計方法——負荷性能圖法[J].化學工程，1997，25(01)： 23-26.

[8]張靜.MATLAB在控制系統中的應用[M].北京： 電子工業出版社，2007.

[9]張曉華.控制系統數字仿真與CAD[M].北京： 機械工業出版社，2004.

[10]王樹清.工業過程控制工程[M].北京： 化學工業出版社，2003.

[11]潘立登，潘仰東.系統辨識與建模[M].北京： 化學工業出版社，2004.