姜睿 呂清林 霍穩(wěn)周 包洪洲 張霞

摘 ?????要：運(yùn)用Aspen軟件建立閃蒸分離模塊采用9種不同的物性方程計(jì)算氫氣在1，4丁二醇中的溶解度。計(jì)算中的溫度范圍為298.15～373.15 K，壓力為1.542～9.8 MPa。對(duì)比模擬值和文獻(xiàn)中實(shí)驗(yàn)值，分析不同物性方程對(duì)于氫氣在1，4丁二醇中的溶解度計(jì)算的使用范圍。考察了進(jìn)料組分流量比對(duì)計(jì)算過程影響，結(jié)果表明，Aspen軟件計(jì)算過程有較高的穩(wěn)定性。優(yōu)選方程對(duì)氫氣在1，4丁二醇內(nèi)溶解度進(jìn)行計(jì)算，分析結(jié)果對(duì)1，4丁二醇液相加氫精制過程單元設(shè)計(jì)及工藝條件選擇有意義。

關(guān) ?鍵 ?詞：Aspen模擬;氫氣溶解度;1，4丁二醇;加氫精制

中圖分類號(hào)：TQ 015?????????文獻(xiàn)標(biāo)識(shí)碼：?A ???????文章編號(hào)： 1671-0460（2019）11-2691-04

Calculation of Hydrogen Solubility in 1，4-Butanediol Based on Aspen Plus

JIANG Rui， LV Qing-lin， HUO Wen-zhou， BAO Hong-zhou， ZHANG Xia

（Sinopec Dalian Research Institute of Petroleum and Petrochemicals， Liaoning Dalian116041， China）

Abstract： A simulation model of flash separation unit was established by Aspen Plus to calculate the solubility of hydrogen in 1，4-butanediol with nine different physical properties equations. The solubilities were calculated at the temperature of 298.15～373.15K and the pressure of?1.542～9.8?MPa. The calculated data were?compared with the measured?value in the literature， and the application range of different physical properties equations for the calculation of the solubility of hydrogen in 1，4-butanediol was analyzed. The influence of feed component flow ratio on the calculation process was investigated. The results?showed?that the Aspen software had?high stability in the calculation process. The optimum method?was?used to calculate the solubility of hydrogen in 1，4-butanediol.?The analysis results have guiding significance for the design and process condition selection of 1，4-butanediol hydrofining.

Key?words： Aspen simulation; Solubility of hydrogen; 1，4-butanediol; Hydrofining

1，4丁二醇（BDO）是一種重要的有機(jī)化工產(chǎn)品，在聚氨酯、聚對(duì)苯甲二酸丁二醇酯、四氫呋喃等重要化工原料生產(chǎn)中有重要作用，廣泛應(yīng)用在化工、紡織、醫(yī)藥、造紙及汽車等多個(gè)領(lǐng)域中。近年來，隨著國內(nèi)化工產(chǎn)業(yè)發(fā)展，1，4丁二醇生產(chǎn)及需求也在逐年增加，特別是高純度，高質(zhì)量的1，4丁二醇，受到國內(nèi)外市場(chǎng)的青睞^[1]。

目前國內(nèi)在1，4丁二醇生產(chǎn)中，炔醛法與順酐法為主要的生產(chǎn)工藝。炔醛法反應(yīng)過程為乙炔與甲醇反應(yīng)生成1，4-丁炔二醇，丁炔二醇加氫生產(chǎn)1，4-丁二醇。順酐法生產(chǎn)1，4-丁二醇反應(yīng)過程為順酐經(jīng)過酯化生成馬來酸二甲酯，馬來酸二甲酯加氫生成1，4-丁二醇。然而這兩種方法生產(chǎn)過程中都會(huì)生成環(huán)狀縮醛雜質(zhì)，縮醛雜質(zhì)與1，4丁二醇沸點(diǎn)相近，不容易分離。同時(shí)，縮醛為后續(xù)加工中的成色物質(zhì)前體，對(duì)產(chǎn)品質(zhì)量影響很大，在1，4丁二醇精制生產(chǎn)中需要加以去除，單純的循環(huán)精餾精制1，4丁二醇會(huì)消耗大量能源及造成大量1，4丁二醇損耗，較難獲得國標(biāo)優(yōu)級(jí)品的1，4-丁二醇產(chǎn)品。

目前有效的精制BDO方法是利用加氫的反應(yīng)方式將粗BDO內(nèi)的縮醛進(jìn)行加氫精制，生成1，4-丁二醇，去除雜質(zhì)，得到高品質(zhì)BDO產(chǎn)品^[2，3]。其加氫過程為液相加氫，在此過程中，氫氣在BDO內(nèi)溶解度為BDO加氫精制的關(guān)鍵參數(shù)。然而，目前對(duì)有機(jī)物內(nèi)氫氣溶解度的模擬及計(jì)算大多集中于苯類，正庚烷及汽柴油等烴類中^[4-6]，對(duì)煉油工業(yè)中加氫過程或液體儲(chǔ)氫過程進(jìn)行指導(dǎo)^[7，8]。然而，對(duì)于氫氣在醇類內(nèi)溶解度，特別是二醇類內(nèi)溶解度的模擬及計(jì)算研究較少^[9]，文獻(xiàn)^[10]中對(duì)氫氣在BDO內(nèi)溶解度實(shí)驗(yàn)測(cè)量值并不多，在較大的溫度及壓力范圍內(nèi)覆蓋數(shù)據(jù)點(diǎn)較少，其對(duì)BDO精制過程單元設(shè)計(jì)指導(dǎo)性不高，從物性模型計(jì)算手段對(duì)氫氣在BDO內(nèi)溶解度研究有重要意義。

在本文研究中，利用Aspen Plus軟件，對(duì)氫氣在BDO內(nèi)溶解過程建立模擬流程，并對(duì)溶解度進(jìn)行了模擬計(jì)算。模擬過程中，選用了PR、PENG-ROB、SRK、PSRK、NRTL及UNQUAC等9種物性方程對(duì)氫氣在BDO內(nèi)溶解度進(jìn)行了計(jì)算，對(duì)比模擬值與文獻(xiàn)中實(shí)驗(yàn)值，總結(jié)Aspen軟件計(jì)算范圍，分析氫氣在BDO內(nèi)隨溫度及壓力變化趨勢(shì)，選取不同溫度壓力下下適用的模擬計(jì)算方法，為BDO精制工藝分析及裝置設(shè)計(jì)提供指導(dǎo)。

1 ?模擬流程的建立

利用Aspen Plus V8.4 建立模擬流程，采用分離過程單元操作，以閃蒸分離器作為操作單元進(jìn)行氣液相平衡分離計(jì)算。模擬步驟如下：

（1）啟動(dòng)Aspen Plus V8.4，選擇模板General with Metric Units。

（2）進(jìn)入Components選項(xiàng)中，輸入組分，在Method模塊中，選擇物性方程。

（3）在Simulation模塊內(nèi)，建立如圖1內(nèi)所示流程，建立閃蒸流程，輸入物流參數(shù)，運(yùn)行模擬流程，在液相物料計(jì)算結(jié)果中可以得到氫氣在BDO內(nèi)的溶解摩爾分率。

2 ?結(jié)果與討論

2.1 ?不同模型計(jì)算結(jié)果對(duì)比

在Aspen Plus的模擬計(jì)算中，不同物性方法的選取對(duì)于計(jì)算結(jié)果有明顯區(qū)別。在石油化工工業(yè)VLE體系中，Aspen軟件推薦采用PR，SRK及PSRK方法進(jìn)行計(jì)算。在煉油工業(yè)中富氫系統(tǒng)及加氫精制計(jì)算中，Aspen軟件推薦采用PENG-ROB，GRAYSON及RK-SOAVE方法進(jìn)行計(jì)算。然而煉油工業(yè)內(nèi)加氫工藝多針對(duì)烴類物質(zhì)，而對(duì)于醇類反應(yīng)而言，Aspen軟件推薦采用WILSON，NRTL及UNIQUAC等方法進(jìn)行計(jì)算?？梢钥闯觯瑲錃庠诖碱悆?nèi)溶解度，特別針對(duì)BDO來說，并沒有一個(gè)固定的推薦方法。為此，本文分別采用以上9種方法，分別對(duì)氫氣在BDO內(nèi)溶解度進(jìn)行了計(jì)算。

在Aspen Plus中設(shè)定具體操作條件如下，以氫氣及1，4丁二醇為主要組分，進(jìn)料溫度選擇298.15～373.15 K，壓力選擇1.542～9.89 MPa，氫氣進(jìn)料量為1 kmol/h，1，4丁二醇進(jìn)料量為1 kmol/h。閃蒸器內(nèi)操作條件設(shè)定為與物料溫度壓力一致的條件，進(jìn)行分離計(jì)算。其不同物性模型方法計(jì)算結(jié)果及其與文獻(xiàn)中實(shí)驗(yàn)值比較如圖2-4所示。

不同的物性方法得到的模擬計(jì)算值之間存在很大差異，其值與實(shí)際測(cè)量值之間偏差有明顯不同。從圖1中看出，在298.15 K時(shí)，NRTL方法與WILSON方法計(jì)算結(jié)果相同，其結(jié)果與實(shí)驗(yàn)測(cè)定值偏差最小，平均偏差-19.34%。從圖2中看出，當(dāng)溫度增加到323.15 K時(shí)，NRTL與WILSON方法計(jì)算結(jié)果與實(shí)驗(yàn)測(cè)定值偏差增加至-29.2%。而PR-BM方法計(jì)算結(jié)果與實(shí)驗(yàn)測(cè)定值偏差最小，其值為-15.5%。當(dāng)溫度為373.15 K時(shí)，依舊為PR-BM方法計(jì)算結(jié)果與實(shí)驗(yàn)測(cè)定值偏差最小，其值為14.6%。

綜上，NRTL，WILSON及PR-BM三種方法對(duì)于氫氣在1，4丁二醇內(nèi)溶解度計(jì)算較為合適，其中NRTL及WILSON適用于相對(duì)低溫（298.15 K）下氫氣在1，4丁二醇內(nèi)溶解度計(jì)算，PR-BM方法適用于相對(duì)高溫（323.15～373.15 K）下氫氣在1，4丁二醇內(nèi)溶解度計(jì)算。

2.2 ?氫氣及1，4丁二醇流率比例對(duì)模擬結(jié)果影響

在之前學(xué)者研究中^[7，11]，計(jì)算氫氣及1，4丁二醇流股的比例偏差過大，軟件計(jì)算結(jié)果容易造成偏差，其在選擇好計(jì)算方法基礎(chǔ)上，對(duì)物性模擬方法計(jì)算比例也應(yīng)該進(jìn)行確定。

為此，在之前利用Aspen軟件建立的閃蒸模擬流程中，以氫氣與1，4丁二醇兩種組分作為閃蒸單元的進(jìn)料組分，兩種物料進(jìn)料溫度為373.15?K、壓力為5.25?MPa，閃蒸器內(nèi)溫度與壓力同進(jìn)料溫度壓力一致。物性計(jì)算方法選用PR-BM，兩種組分比率不同，在氫氣與BD進(jìn)料摩爾比率為3000∶1～1∶

1?000條件下，對(duì)模擬過程進(jìn)行計(jì)算，具體結(jié)果如表1所示。

根據(jù)計(jì)算結(jié)果可以發(fā)現(xiàn)，在計(jì)算范圍內(nèi)，液相物料內(nèi)的計(jì)算結(jié)果并沒有發(fā)生明顯變化，其中當(dāng)氫氣與1，4丁二醇摩爾比大于129時(shí)，氫氣完全溶解在液相物料中，計(jì)算結(jié)果也很好的反應(yīng)了該過程。同之前ChemCAD軟件的計(jì)算結(jié)果^[11]相比，Aspen軟件在物性數(shù)據(jù)計(jì)算中并沒有出現(xiàn)因進(jìn)料組分比例過大引起的計(jì)算誤差，其計(jì)算結(jié)果有較高的穩(wěn)定精度。通過計(jì)算結(jié)果，也可以說明Aspen 軟件適用于1，4丁二醇精制過程中計(jì)算。

2.3 ?氫氣在1.4丁二醇內(nèi)溶解度變化規(guī)律

由于1，4丁二醇熔點(diǎn)為293.35 K，一般將其加熱以保障1，4丁二醇在裝置內(nèi)流動(dòng)性，為此，本文結(jié)合實(shí)際工況，在298.15 K條件下，利用NRTL方法計(jì)算，壓力范圍2～9 MPa下氫氣在1，4丁二醇內(nèi)溶解度。在323.15～373.15 K條件下，利用PR-BM方法計(jì)算壓力范圍2～9 MPa下氫氣在1，4丁二醇內(nèi)溶解度。其中氫氣進(jìn)料量與1， 4丁二醇進(jìn)料量比值為1∶1，閃蒸器內(nèi)操作條件設(shè)定為與物料溫度壓力一致的條件，進(jìn)行分離計(jì)算。計(jì)算結(jié)果如表2所示。

通過計(jì)算數(shù)據(jù)可以看出，在相同溫度下，隨著壓力增加，氫氣在1，4丁二醇內(nèi)溶解度也隨之增加。在相同壓力下，隨著溫度增加，氫氣在1，4丁二醇內(nèi)溶解度同樣增加。氫氣在1，4丁二醇內(nèi)溶解度內(nèi)溶解趨勢(shì)與在烴類內(nèi)溶解趨勢(shì)相同，均隨溫度壓力升高而增加。對(duì)于1，4丁二醇精制流程來說，先將1，4丁二醇內(nèi)進(jìn)行預(yù)熱到反應(yīng)溫度，再進(jìn)行充氫溶解，對(duì)于提高氫氣在1，4丁二醇內(nèi)溶解度，進(jìn)而提高縮醛的催化轉(zhuǎn)化，提高1，4丁二醇產(chǎn)品純度是非常有利的。

3??結(jié)論

（1）利用Aspen軟件建立閃蒸分離模塊采用9種物性方程來計(jì)算氫氣在1，4丁二醇中的溶解度，并對(duì)計(jì)算結(jié)果與文獻(xiàn)實(shí)驗(yàn)值進(jìn)行對(duì)比。發(fā)現(xiàn)在298.15?K時(shí)，NRTL方法與WILSON方法計(jì)算結(jié)果相同，其結(jié)果與實(shí)驗(yàn)測(cè)定值偏差最小。而在325.15～373.15?K區(qū)間，PR-BM方法計(jì)算結(jié)果與實(shí)驗(yàn)測(cè)定值偏差最小。

（2）對(duì)物料進(jìn)料組分流率對(duì)模擬計(jì)算結(jié)果影響進(jìn)行分析，Aspen軟件在較大物料進(jìn)料比下同樣保持很高的計(jì)算穩(wěn)定精度，完全可以滿足1，4丁二醇的精制流程計(jì)算。

（3）利用優(yōu)選的模擬方法對(duì)不同溫度及壓力條件下，氫氣在1，4丁二醇內(nèi)溶解度進(jìn)行了計(jì)算，發(fā)現(xiàn)隨著溫度及壓力的增加，氫氣在1，4丁二醇內(nèi)溶解度也隨之升高，對(duì)1，4丁二醇液相加氫精制過程提出工藝優(yōu)化指導(dǎo)。

參考文獻(xiàn)：

[1]李玉芳. 1，4-丁二醇的生產(chǎn)技術(shù)及國內(nèi)外市場(chǎng)分析[J]. 化學(xué)工業(yè)， 2012， 30（10）： 39-44.

[2]付秋紅， 張寶國， 張艷俠， 等. 高活性醛加氫催化劑研究[J]. 當(dāng)代化工， 2015， 11： 2569-2571.

[3]伍德， 威萊特， 克雷， 等. 1，4丁二醇的純化方法： 中國， CN 1216973A[P]. 1999-05-19.

[4]Aslam R， Mu?Ller K， Mu?Ller M， et al. Measurement of hydrogen solubility in potential Liquid Organic Hydrogen Carriers[J]. Journal of chemical & engineering data， 2015， 61（1）： 643-649.

[5]Ronze D， Fongarland P， Pitault I， et al. Hydrogen solubility in straight run gasoil[J]. Chemical Engineering Science， 2002， 57（4）： 547-553.

[6]羅化峰， 郭劍虹， 凌開成， 等. 氫氣在烴類混合溶劑中高壓溶解度的測(cè)定[J]. 煤炭轉(zhuǎn)化， 2011， 34（2）： 55-58.

[7]韓興華， 董婷， 白昊， 等. 基于 Aspen Plus 的氫氣在烴類中溶解度的計(jì)算[J]. 中北大學(xué)學(xué)報(bào) （自然科學(xué)版）， 2016， 5： 16.

[8]劉晨光， 闕國和. 氫氣在石油餾分中溶解度的測(cè)定[J]. 煉油設(shè)計(jì)， 1999， 29（5）： 33-36.

[9]Safamirzaei M， Modarress H， Mohsen-Nia M. Modeling the hydrogen solubility in methanol， ethanol， 1-propanol and 1-butanol[J]. Fluid Phase Equilibria， 2010， 289（1）： 32-39.

[10]Brunner E. Solubility of hydrogen in diols and their ethers[J]. The Journal of Chemical Thermodynamics， 1980， 12（10）： 993-1002.

[11]王世麗， 劉永剛， 張瑞芹. 運(yùn)用 ChemCAD 計(jì)算氫氣在烴類中的溶解度[J]. 現(xiàn)代化工， 2013， 5： 114-118.