2010年10月Oct 2010第33卷第5期Large Scale Nitrogenous Fertilizer IndustryVol, 33 No5基于非平衡級(jí)速率模型的低溫甲醇洗吸收塔模擬研究楊亭1王東巖2關(guān)威1張雙雙1常虹張述偉1(1.大連理工大學(xué)化工學(xué)院遼寧大連,116012;2陜西延長(zhǎng)石油集團(tuán)榆林煤化公司陜西榆林,719000)摘要針對(duì)低溫甲醇洗吸收塔的特點(diǎn)建立了采用新物性計(jì)算方法的非平衡級(jí)速率模型,與傳統(tǒng)采用PSRK( Predictive soave-Rich- Kwong)物性計(jì)算方法的平衡級(jí)模型模擬結(jié)果和實(shí)際生產(chǎn)數(shù)據(jù)進(jìn)行對(duì)比,對(duì)低溫甲醇洗吸收塔進(jìn)行模擬研究。研究表明:通過(guò)采用新的物性計(jì)算方法并優(yōu)化物性組分參數(shù),非平衡級(jí)速率模型較好地反映了實(shí)際工況,證明該模型可以應(yīng)用于低溫甲醇洗系統(tǒng)的工藝計(jì)算、實(shí)際生產(chǎn)的模擬分析以及操作條件的優(yōu)化,為此裝置設(shè)計(jì)提供可靠技術(shù)支持。關(guān)鍵詞低溫甲醇洗吸收塔非平衡級(jí)速率模型參數(shù)優(yōu)化近年來(lái),隨著我國(guó)煤化工熱潮的掀起,工業(yè)上有收斂速度快適應(yīng)強(qiáng)非理想系物系的特點(diǎn),但對(duì)脫硫脫碳的方法很多,其中冷法以低溫甲醇洗迭代初值有嚴(yán)格的要求,且計(jì)算量大。非平衡級(jí)Rectisol為代表,熱法以聚乙二醇二甲醚法 Selexol速率模型將分離過(guò)程視為一個(gè)速率過(guò)程,取代了工藝最為著名低溫甲醇洗技術(shù)是20世紀(jì)50年平衡過(guò)程,基本假設(shè)是在塔板上氣液相主體并不代由德國(guó)林德公司和魯奇公司共同開(kāi)發(fā)的一種氣處于平衡狀態(tài),相平衡僅存在于氣液相界面處。該體凈化方法,目前全世界共有低溫甲醇洗裝置幾模型直接用傳質(zhì)、傳熱速率方程表征傳遞過(guò)程,從十佘套。低溫甲醇洗工藝凈化方法是一種典型的內(nèi)在的傳遞過(guò)程機(jī)理來(lái)揭示精餾特性,避免使用物理吸收酸性組分的方法,其優(yōu)點(diǎn)如下:吸收能級(jí)效率和等板高度(HETP)等難以確定的量,可更力大,可以同時(shí)脫除HS、CO2等多種有害組分;凈準(zhǔn)確預(yù)測(cè)塔內(nèi)的濃度、流量和溫度分布,為精餾分化度高凈化氣體中含硫量小于0000%,C02含離過(guò)程的設(shè)計(jì)和模擬開(kāi)辟了新的道路們。本文以某量小于000%;吸收選擇性好,低溫時(shí)HS和CO2廠低溫甲醇洗吸收塔實(shí)際生產(chǎn)數(shù)據(jù)為依據(jù),采用在甲醇中的溶解度大,吸收溶液的循環(huán)量小;甲醇非平衡級(jí)速率模型,并采用新的物性計(jì)算方法及的熱穩(wěn)定性和化學(xué)穩(wěn)定性很好,吸收酸性氣體后修正物性組分參數(shù),準(zhǔn)確地模擬了低溫甲醇洗工不發(fā)生降解;在低溫下甲醇黏度小,具有良好的傳藝系統(tǒng)的吸收塔,為該系統(tǒng)的工藝模擬提供了重?zé)岷蛡髻|(zhì)性能。低溫甲醇洗工藝凈化法被廣泛應(yīng)要的技術(shù)支持。用于以渣油或煤為原料的制氨以及城市煤氣化等工業(yè)上,對(duì)該工藝的準(zhǔn)確模擬對(duì)于實(shí)際工廠操作1非平衡級(jí)模型和設(shè)計(jì)至關(guān)重要。速率平衡過(guò)程視為質(zhì)量和熱量的速率傳遞過(guò)長(zhǎng)期以來(lái),化工分離設(shè)備的設(shè)計(jì)與模擬仍以程,雙膜理論被考慮在氣液相之間,平衡條件僅僅平衡級(jí)模型為主,通過(guò)板效率和等板高度對(duì)計(jì)算假設(shè)在它的界面處,膜中的質(zhì)量和熱量計(jì)算需考結(jié)果進(jìn)行校正。而板效率和等板高度這些因素往慮板式塔或填料塔的特點(diǎn)、氣液量負(fù)荷、不同的傳往不能由理論方法確定,使得其用于實(shí)際操作和設(shè)計(jì)時(shí)受到很大限制4。1985年 Kirshnarllurthy等收稿日期:201003-10收到修改稿日期:2010-10-15作者簡(jiǎn)介:楊亭,女,1986年出生,大連理工大學(xué)化T學(xué)院化學(xué)建立了在傳質(zhì)速率方程基礎(chǔ)上的非平衡級(jí)速率模T藝專(zhuān)業(yè)在讀堿士,聯(lián)系電話(huà)011389390m96型,并使用 Newon- RaphSon法進(jìn)行求解,該方法具 newcomen2010年第33卷遞性質(zhì)等}。非平衡級(jí)速率模型假設(shè):①穩(wěn)態(tài)操作;于計(jì)算傳質(zhì)傳熱系數(shù)和相間面積關(guān)系式的計(jì)算方②氣相到達(dá)上一級(jí)前完全混合,液相為全混流;③法正確與否,是決定該塔模擬計(jì)算是否正確的關(guān)在相界面上氣液兩相達(dá)到相平衡;④在每一級(jí)塔鍵。模型中界面面積因子、熱傳遞因子以及滯留比板上傳質(zhì)速率相等。其非平衡級(jí)速率模型建立如例因子需要不斷地調(diào)整。本文對(duì)該吸收塔精餾過(guò)圖1所示。程中的傳質(zhì)系數(shù)計(jì)算采用ACHE(1958)關(guān)聯(lián)式傳熱系數(shù)計(jì)算采用 Chiltion and Colbun0關(guān)聯(lián)式相界面積計(jì)算采用 Scheffe and Weiland關(guān)聯(lián)式r2實(shí)際工況模擬液相21工藝流程典型的低溫甲醇洗流程包括吸收塔、閃蒸罐、汽提塔和甲醇再生塔。其中吸收塔作為低溫甲醇洗系統(tǒng)的核心設(shè)備,承擔(dān)著脫除工藝氣體中H2CO2的任務(wù)。變換氣依次通過(guò)吸收塔的脫硫段和脫碳段,最終H2S和CO2的摩爾分率可滿(mǎn)足如下設(shè)圖1非平衡級(jí)速率模型計(jì)控制指標(biāo):HS含量小于0.0000%,C02含量小非平衡級(jí)速率模型對(duì)于板式塔或填料塔,認(rèn)于000%。從吸收塔塔頂出來(lái)的凈化氣大部分去為它由一序列非平衡段組成。對(duì)每一非平衡段可甲醇合成或者氨合成,部分去城市煤氣系統(tǒng)以寫(xiě)出模型方程組如下。吸收塔分為4段,最下段(B1)為脫硫段,其余物料平衡方程3段(B2、B3、B4)為脫碳段。在脫硫B1段原料氣經(jīng)氣相F+V+1y-VD,+N,=0(1)富含CO2的甲醇液(11號(hào)流股)洗滌,脫除HS液相F+L-Lx,+M=0(2)COS和部分CO2等組分后進(jìn)入脫碳B2段脫除傳質(zhì)速率方程CO2,進(jìn)入脫碳段的氣體不含硫,接著繼續(xù)進(jìn)入脫碳B3段進(jìn)一步脫除CO2組分。從脫碳B3段頂部氣相N∑.m4(m-m)-ym∑N=0(3)出來(lái)仍含有部分CO2的氣體再進(jìn)入脫碳B4段,在脫碳B4段塔頂用貧甲醇液(4號(hào)流股)洗滌,凈化液相N-∑4m,(一貼,)N0(4)氣(3號(hào)流股)由塔頂引出。吸收塔中間設(shè)有兩個(gè)中膜界面汽液平衡方程:y-Kx=0(間冷卻器(E1和E2),分別冷卻富甲醇液5號(hào)和號(hào)流股。該吸收塔如圖2所示。液膜能量平衡方程:Q-Q=0(6)氣膜能量平衡方程:Q-Q=0(7)界面摩爾分率歸一化方程:與傳統(tǒng)平衡級(jí)模型不同,非平衡級(jí)速率模型將分離過(guò)程視為質(zhì)量和熱量的速率傳遞過(guò)程,并非平衡過(guò)程。該方法包含氣液相傳質(zhì)傳熱平衡模型、用于確定相間傳遞速率的氣液相傳質(zhì)傳熱速率模型、表面氣液平衡關(guān)系式和用于計(jì)算傳質(zhì)傳熱系數(shù)和相間面積的關(guān)系式。2吸收塔流程在采用非平衡級(jí)速率模型的計(jì)算中,氣液相1~14—流股;B—塔段;E—換熱器傳質(zhì)傳熱平衡模型、用于確定相間傳遞速率的氣22原料氣組成液相傳質(zhì)傳熱速率模型、表面氣液平衡關(guān)系式、用進(jìn)人該系統(tǒng)的原料氣2號(hào)流股中,除含有N2和第5期楊亭等,基于非平衡級(jí)速率模型的低溫甲醇洗吸收塔模擬研究H2之外,還含有CO、H0、HS和CO2等雜質(zhì)。原料氣組成見(jiàn)表1。表1原料氣組成摩爾流量/組成(摩爾分?jǐn)?shù)),%溫度/K壓力/MPakmolh- CH,OH CO, H,7.6176152.60.0862.420.1023物性計(jì)算方法選擇經(jīng)驗(yàn)進(jìn)行估測(cè)。非平衡級(jí)速率模型的模擬避免了物性計(jì)算方法是否恰當(dāng),直接影響到計(jì)算結(jié)這些不確定的因素,直接計(jì)算相流動(dòng)接觸間的質(zhì)果的精確,進(jìn)而影響模擬結(jié)果的準(zhǔn)確性。針對(duì)本量、熱量傳遞。該模型的模擬遠(yuǎn)比平衡級(jí)模型要復(fù)文低溫甲醇洗系統(tǒng)吸收塔的特點(diǎn),發(fā)現(xiàn)傳統(tǒng)平衡雜,需要提供相應(yīng)塔的結(jié)構(gòu)數(shù)據(jù),并且可以同時(shí)進(jìn)級(jí)模型采用PSRK物性方法計(jì)算的其吸收塔體系行塔的水力學(xué)計(jì)算。其流程建立的單元操作模型有諸多誤差,比如液相密度、焓以及相平衡。本文見(jiàn)表3。在對(duì)吸收塔采用非平衡級(jí)速率模型的同時(shí),采用表3非平衡級(jí)速率模型中用到的單元操作模塊新的物性計(jì)算方法并通過(guò)修正其物性組分參數(shù)131塔段來(lái)計(jì)算其蒸汽壓、液相密度和相平衡。然后以工廠塔板數(shù)提供的實(shí)際生產(chǎn)數(shù)據(jù)(見(jiàn)表2)為依據(jù),將塔底1號(hào)塔板類(lèi)型扁平泡罩扁平泡罩扁平泡罩齒形泡塔徑/2450流股的模擬結(jié)果與實(shí)際生產(chǎn)數(shù)據(jù)做對(duì)比,從而確泡罩直徑/mm80定吸收塔模擬的物性計(jì)算方法。塔板型式單流道雙流道表21號(hào)流股的模擬結(jié)果與實(shí)際數(shù)據(jù)對(duì)比溢流堰高/mm1號(hào)流股計(jì)算值物性計(jì)算方法實(shí)際值2.5模擬結(jié)果PSRK新物性溫度/K通常教硏室在對(duì)低溫甲醇洗系統(tǒng)的模擬研究采壓力/MPa76177.6177617用SRK物性方法及傳統(tǒng)平衡級(jí)模型,與其他方法密度/kg:m3相比,模型結(jié)果與實(shí)際值較吻合121,但也存在黏度/gcm-20.006定的誤差。因此,在教研室原來(lái)模擬研究的基礎(chǔ)熱力學(xué)導(dǎo)熱系數(shù)/W·mk-10.3940.147上,本文繼續(xù)研究探討,能否采用其新物性方法和表面張力/mN·m1非平衡級(jí)速率模型模擬,使其模擬結(jié)果與實(shí)際生從表2可以看出,采用新物性計(jì)算方法,塔底產(chǎn)數(shù)據(jù)的誤差更小。對(duì)吸收塔(圖2)的模擬研究1號(hào)流股的溫度偏差008%,密度偏差3.31%,黏中,本文建立了兩種模型,模型Ⅰ為吸收塔B1、度偏差333%,熱力學(xué)導(dǎo)熱系數(shù)偏差2.04%;而采B2、B3段采用新物性計(jì)算方法和非平衡級(jí)速率模用PSRK物性方法的計(jì)算值與實(shí)際數(shù)據(jù)相比,其結(jié)型,B4段由于其貧甲醇溶液(4號(hào)流股)的進(jìn)料溫果相差很大,如溫度偏差0.19%,密度偏差度為-57.7℃,而該塔的塔壓為75.5MPa,因此屬于14.38%,黏度偏差83.33%,熱力學(xué)導(dǎo)熱系數(shù)偏差非常典型的高壓低溫體系,所以該段仍采用PSRKl68.03%。因此,從采用新物性計(jì)算方法和PSRK物性方法和非平衡級(jí)速率模型;模型Ⅱ?yàn)槲账镄苑謩e計(jì)算出的流股結(jié)果與實(shí)際值的比較可以全流程(包括B1B2、B3、B4段)全采用PSRK物性看出,采用新物性計(jì)算方法得出的流股密度、黏方法和傳統(tǒng)平衡級(jí)模型進(jìn)行模擬研究。度、導(dǎo)熱系數(shù)以及表面張力值同實(shí)際值非常吻合,通過(guò)采用新的物性計(jì)算方法并優(yōu)化物性組分說(shuō)明該新物性計(jì)算方法相對(duì)于FSRK物性方法更參數(shù)以及模型參數(shù),將該模型Ⅰ應(yīng)用到吸收塔適合低溫甲醇洗體系計(jì)算。的模擬中,得出各個(gè)流股和該塔的模擬計(jì)算結(jié)果24非平衡級(jí)速率模型的單元操作模塊將其模擬結(jié)果與模型Ⅱ的計(jì)算結(jié)果及實(shí)際生產(chǎn)數(shù)在傳統(tǒng)平衡級(jí)模型的模擬計(jì)算中,常需要應(yīng)據(jù)對(duì)比,結(jié)果如表4所示。用板效率和HETP等因素來(lái)評(píng)價(jià)板式塔和填料塔從表4可見(jiàn),對(duì)于流股中組分N2、CO、H2O的性能。而對(duì)于板效率和理論級(jí)當(dāng)量高度這些因HS的組成含量而言,模型Ⅰ和模型Ⅱ的計(jì)算結(jié)果素來(lái)說(shuō),往往沒(méi)有理論的方法來(lái)確定,只能靠實(shí)際非常相近,并且兩種計(jì)算結(jié)果與實(shí)際值都能很好2010年第33卷地吻合。而對(duì)于流股的溫度、流量和流股中其他組為0.08%、溫度偏差為0.2K、CO2純度偏差為分含量來(lái)說(shuō),兩種模型的計(jì)算結(jié)果存在一定的差0.15%;塔頂氣相流率(3號(hào)物流)偏差為0.09%、溫距。而模型Ⅰ的計(jì)算結(jié)果同實(shí)際生產(chǎn)數(shù)據(jù)非常吻度偏差為0.IK、H2純度偏差為0.08%;側(cè)線1l號(hào)合,模型Ⅱ的計(jì)算結(jié)果同實(shí)際值相差較大。采用模物流流率偏差為0.22%、溫度偏差為O.3K、CO2純型Ⅱ計(jì)算時(shí),除塔底物流1號(hào)流股和塔頂凈化氣3度偏差為0.52%;側(cè)線7號(hào)物流流率偏差為號(hào)流股的計(jì)算值與實(shí)際值很接近外,其余流股的0.02%、溫度偏差為36K、CO2純度偏差為132%計(jì)算值與實(shí)際值都有一定的差距,如塔段循環(huán)物側(cè)線5號(hào)物流流率偏差為0.08%、溫度偏差為3.9流5號(hào)和7號(hào)流股數(shù)據(jù)與實(shí)際值相差較大,特別K、CO2純度偏差為3.25%。并且,由表2可以得知是流股中組分CO2的含量采用模型Ⅰ對(duì)流股的密度、黏度、導(dǎo)熱系數(shù)等物性采用模型Ⅰ對(duì)低溫甲醇洗的吸收塔進(jìn)行模擬值進(jìn)行模擬計(jì)算時(shí),其結(jié)果與實(shí)際值也都符合良計(jì)算,各流股的溫度、壓力、流量以及關(guān)鍵組分的好。這說(shuō)明模型Ⅰ充分反應(yīng)了該實(shí)際工況,并證明摩爾組成值都與實(shí)際值非常吻合。例如,由模型Ⅰ模型Ⅰ對(duì)低溫甲醇洗的模擬是可行的,為后續(xù)低得計(jì)算結(jié)果,吸收塔塔底液相流率(1號(hào)物流)偏差溫甲醇洗系統(tǒng)流程模擬奠定了基礎(chǔ)。表4模擬結(jié)果與實(shí)際值對(duì)比流股模型溫度K流量組分組成(摩爾分率),%kmolII CH,OHH2模型I266.03111.60.6510.3實(shí)際值26623109.10.6500.00l模型Ⅱ26573110.70.6430.341001100010.0010.0020.000模型I216.1397040.0000.00009530.0380000.0020.002實(shí)際值2162396660.0000.0000.9530.0370.0000.0020.0020.0000.006模型Ⅱ216.03953.20.0000.0000.95300380.0000.0020.0020.0000模型I249.3482970.8610.12200110.0010.0000.00000000.0030001實(shí)際值2530.8610.1260.00800.0000.0000.0000.000模型Ⅱ模型I:際值250模型Ⅱ25305104.20.8150.160.0000.0030.0009模型I267.06314.10.6590.3280.0090.0010.0000.0000.0000實(shí)際值266.76312906600.3270.0090.0010.0000.0000.0000020.000模型Ⅱ263.8627590.6640.32100110.00100000.0000.0000.002000模型I26703068.60.6590.3280.0090.0010.0000.0000.0000.0020.000實(shí)際值266.730610.000000.0020.000模型Ⅱ263.83012.50.6640.3210.0110.0010.000.0000.0000.0020.0003結(jié)論的借鑒意義。通過(guò)以上模擬分析,可得出以下結(jié)論3)和傳統(tǒng)平衡模型模擬相比,采用非平衡級(jí)1)在流程模擬中,選擇合適的物性方法是決速率模型模擬的計(jì)算求解和收斂難度增大。非平定模擬結(jié)果準(zhǔn)確的關(guān)鍵。本文通過(guò)采用新的物性衡級(jí)速率模型的模擬避免了使用平衡級(jí)模型中難計(jì)算方法并優(yōu)化物性組分參數(shù),使其液相密度、黏以準(zhǔn)確計(jì)算的板效率和HETP,其模擬計(jì)算的精確度、導(dǎo)熱系數(shù)以及相平衡的計(jì)算準(zhǔn)確,從而建立了度和可靠性增加,調(diào)整參數(shù)比平衡級(jí)復(fù)雜,更接近低溫甲醇洗系統(tǒng)的熱力學(xué)模型。實(shí)際工廠數(shù)據(jù)。2)實(shí)際工況的模擬表明:應(yīng)用新的物性計(jì)算符號(hào)說(shuō)明方法,采用優(yōu)化參數(shù)后的非平衡級(jí)速率模型,低溫a—界面面積,m2甲醇洗系統(tǒng)吸收塔的模擬結(jié)果與實(shí)際值符合良F一進(jìn)料流量, kmol/好,對(duì)低溫甲醇洗后續(xù)系統(tǒng)的模擬和優(yōu)化有很好k-傳質(zhì)系數(shù),kmol/(mh)第5期楊亭等,基于非平衡級(jí)速率模型的低溫甲醇洗吸收塔模擬研究L一液相流量,kmol/hMulticomponent Separation Processes. AIChE J, 1985, 31(3)N一傳質(zhì)通量,kmol/h449~464K一相平衡常數(shù)6王峰,趙寧,李軍平等,采用非平衡級(jí)模型模擬催化精餾過(guò)Q—一熱負(fù)荷,J/h程的研究,石油化工,2007,36(11):1128-1133V一氣相流量,kmol/h7徐偉齊,湯志剛,周長(zhǎng)城等.非平衡級(jí)模型用于SO2化學(xué)吸收過(guò)x一液相摩爾分?jǐn)?shù)程的模擬清華大學(xué)學(xué)報(bào),2003,43(12):1616-1619y一氣相摩爾分?jǐn)?shù)8 Al-Baghli NA, Pruess SA, Yesavage VF, et al. A Rate-based Model for The Design of Gas Absorbers for The Removal of CO, and上標(biāo)H S Using Aqueous Solutions of MEA and DEA. Fluid Phase E-L一液相quilibria,2001,185(1-2):31-43V一氣相9 AIChE Bubble Tray Design Manual; Prediction of Fractionation相界面Efficieney. American Institute of Chemical Engineers, New York下標(biāo)i,m一組分10 Taylor R, Krishna R Multicomponent Mass Transfer. Wiley-Inter-j—塔板數(shù)參考文獻(xiàn)11 Schelte rd, Weiland RH. Mass Transfer Characteristics of valve1楊誼,張述偉,陳曉峰.“煤代油”低溫甲醇洗工藝流程改造方Trays. Ind. Eng. Chem. Res, 1987, 26(2): 228-236案的研究,大氮肥,2004,27(5):192~29512關(guān)威,張述偉,管風(fēng)寶等.城市煤氣低溫甲醇洗系統(tǒng)凈化工藝模擬與改造研究.化肥設(shè)計(jì),2006,44(6):172張述偉,胡乃平,王長(zhǎng)英等.用低溫甲醇洗凈化不含硫原料氣的13 Poling B, Prausnitz J. O'Connell J. The Properties of Gases and工藝流程簡(jiǎn)化.高?；瘜W(xué)工程學(xué)報(bào),1999,13(6):5Liquids. New York: Me Graw-Hill Professional Publishing, 20003張述偉,陸明亮,徐志武.低溫甲醉洗系統(tǒng)模擬與分析氮肥設(shè)14劉健,張述偉,孫道青.低溫甲醇洗凈化天然氣工藝流程的計(jì),1994,32(1):25-31研究天然氣化工,2007,32(45):47-54魏海國(guó),費(fèi)維揚(yáng)非平衡級(jí)模型及其在化工分離過(guò)程中的應(yīng)用.15 Wilson GM Vapor-Liquid Equilibriun.x. A New Expression for石油化工,2001,30(6):486-490the Excess Free Energy of Mixing. Journal of the American Chem-5 Kirshnamurthy R, Taylor R. A Non-equilibrium Stage Model ofal society,1964,86(2):127-130SIMULATION STUDY ON RECTISOL ABSORBER BASEDON NON-EQUILIBRIUM RATE MODELYang Ting, Guan Wei, Zhang Shuangshuang, Chang Hong, Zhang ShuweiDepartment of Material and Chemical engineering, Dalian University of Technology, Dalian, 116012)Wang dongyan(Shaanxi y anchang Petroleum Group Yulin Coal Chemical Industry Corp., Yulin, 719000)Abstract With the absorber features in Rectisol unit a non-equilibrium stage and rate-based modelemploying a new method of physical property calculation is established to compare with both the equilibriumstage model simulated results by employing a traditional PSRK(Predictive Soave Redlich-Kwong method ofphysical property calculation and actual operating data and to carry out a simulation study on the absorber inRectisol The study shows that by employing the new method of physical property calculation and optimizingthe parameters of physical property components the non -equilibrium stage and rate-based model better reflects the actual operating conditions, certifies that the model may be applied to the process calculations, simulation analysis of actual operation and optimization of operating conditions of the Rectisol system; thus a reliable technical support is provided for the unit designKey words: rectisol, absorber, non-equilibrium stage and rate-based model, parameters optimization,,,,,,,,,,,,,,歡迎投稿、歡迎訂閌、歡迎刊登廣告!