Predton of old Paraffn Preptaton Usng old Phase Equaton of tate Proeedngs of European Congress of Chemal Engneerng (ECCE-6) Copenhagen, 16- eptember 7 Predton of old Paraffn Preptaton Usng old Phase Equaton of tate Mohsen Vafae eft, a Hamd Mehdzadeh, a Al Mousav, b a Department of Chemal Engneerng, Tarbat Modares Unversty, Tehran, Iran b Iran's Researh Insttute of Petroleum Industry, Pazhooheshgah, Blvd.Kharabad Junton, Old Qom Road, Tehran, Iran Abstrat In the tradtonal methods for wax preptaton predton, the fugaty of sold s alulated usng a hypothetal thermodynam yle. In ths study, a new method based on TT sold equaton of state s proposed to alulate the fugates of sold phases. The mult-sold approah s used for desrpton of sold phases. Obtaned parameters for ths equaton of state s useful for future works. Keywords: wax preptaton, sold equaton of state, mult-sold 1. Introduton Wax deposton from gas and ol produton faltes and ppelnes s undesrable. The flow-lnes and proess equpments may be plugged by wax deposton. Dfferent physal and hemal methods have been proposed to remove deposted solds, whh nreases operatng osts. A relable model for wax preptaton alulatons s hghly valued for desgn and operaton of flowlnes. ne 199s, so many efforts have been done to predt ondtons under whh the waxes an preptate and the amount of preptated wax n funtons of pressure, temperature and omposton. At frst, the alulatons were based on the sold-lqud equlbrum assumpton. Later on, the gas phase was nluded n the alulatons [1, ]. There are two learly defned assumptons for the determnaton of the thermodynam equlbrum wax lqud n establshed mult-omponent hydroarbons systems: sold soluton, and the formaton of multple sold phases. In the former ase, dfferent methods were proposed based on the atvty oeffent model assumng the non-dealty of lqud and sold phases [3, 4]. old phase transton and vapor phase were then onsdered n other works [5-7]. The non-dealty was defned usng Wlson or UNIQUAC equatons. Lra-Galeana et al. [8] developed the mult-sold approah n 1996. In ths model, t s assumed that the sold wax onssts of several pure sold phases, n whh the number and nature of them wll be obtaned from phase stablty analyss. Coutnho showed that the sold phase s a mult-sold soluton n nature and supported by expermental data [9].
M. Vafae eft et al. In ths work, the wax-preptaton model based on sold phase equaton of state wll be presented. The mult-sold approah s used beause of ts wde aeptablty and lmtaton of usng sold phase equaton of state. The parameters of the equaton of state are obtaned n the ase that vapor, lqud and sold phases are nluded n the system.. old Phase Equaton of tate The equatons of state that an be appled to a system ontanng sold, lqud and vapor phases smultaneously are sare[1,11]. In ths projet, TT 1 equaton of state s used [11]. The general form of the equaton s as below RT a P = (1) v b ( v + ub)( v + wb) In ths equaton, u and w are 3 and -.5 respetvely. Also, a =.4757R T / P b =.7474RT / P Z =.9696 a = α a (5) Alpha funton for lqud and vapor phases ould be used n onventonal polynomal form or n exponental form. In ths equaton, a new alpha funton s ntrodued for the sold phase, whh wll be dsussed later. In order to alulate fugaty of eah omponent n pure sold state, the followng equaton s used [1]. ΔH v (6) d ln f = dt + dp RT RT Where, f: fugaty Δ H : enthalpy hange n result of hange n system temperature v: partal molar volume By ntegratng the above equaton from the trple pont pressure to the system pressure for lqud and sold phases and dvdng the two equatons, the followng f relaton wll be obtaned. In ths equaton, ΔH and Δ v are supposed to be ndependent of pressure and temperature. f f ΔH T ΔvP ln = 1 L (7) f RT Tr RT () (3) (4) 1 Twu-m-Tassone
Predton of old Paraffn Preptaton Usng old Phase Equaton of tate uppose that the pressure system s equal to zero, we wll have f f ΔH T ln = 1 L f RT Tr (8) Ths equaton s smlar to the equaton that was proposed by Prausntz et al. [13] for the alulaton of sold phase fugaty. If the fugaty s alulated at zero pressure, * * f * * 1 a v + w ln = 1 ln b ln( v 1) ln * * (9) P ( w u) b v + u Where, 1/ * * * * 1 a a a v = u w u + w 4 uw + * * * (1) b b b * Pa a = (11) RT * Pb b = (1) RT And, u and w are the parameters of the equaton of state. By equalzng the two equatons, we wll get * * L f + Δ = H * * 1 a v w f T 1 ln b ln( v 1) ln ln 1 * * + t ( w u) b v u P RT (13) T The only unknown varable n the above equaton s. Therefore, the equaton s solved to obtan a * for dfferent temperatures. The parameters of predefned sold alpha funton are alulated by orrelatng the data to the followng equaton.5 n α ( T ) = 1+ l (1 T ) + m (1 T ) (.7 T ) (14) In whh α = a / a a a = * s r α an be alulated usng R T P r r In order to alulate the optmum parameters, the followng goal funton was used to mnmze the dfferene between the presented data and the alulated data from eq. 9. 1 (1.5 ) (1 ) n (.7 ) (17) f x = + l T + m T T α ( ) ( ) s r, r, r, The smplex-nelder-mead algorthm was utlzed to obtan the optmum parameters, whh mnmzes the goal funton. Due to the nonlnearty of the funton, the results a * (15) (16)
M. Vafae eft et al. wll drastally depend on the ntal guess for the optmal parameters. To avod ths problem, the optmzaton problem s run for dfferent startng ponts. 3. Wax Preptatng Model The followng expressons are used to desrbe the system n whh vapor, lqud and pure solds are n equlbrum state. Mass balane for preptatng omponents: v V L L f F y n + x n + n z n = (18) Where = 1,..., N, j = 1,..., N P N P : Number of preptated solds N C : Number of omponents Mass balane for non-preptatng omponents: v V L L f F y n + x n z n = Equalty of fugates n the lqud and gas phases for eah omponent: V L f f = Equalty of fugates n the lqud and sold phases for preptatng omponents: f f = L ummaton of mole fratons n lqud and gas phases are equal to unty N x 1= = 1 N = 1 y 1= (19) () (1) () (3) All the equatons above onsttute a system of equatons, whh should be solved n order to defne the equlbrum system ompletely. An error funton s ntrodued to hek the onvergene of the system of equatons. NC + N P + (4) f ( δ ) = δ = 1 Where δ 's are the rght hand expressons n equatons 13-17. 4. Results and Dsussons The omposton of ol samples and some synthet mxtures, whh were used n ths projet are shown n tables 1-4. ne the referenes are dfferent, the omponents onsdered for eah sample are not the same.
Predton of old Paraffn Preptaton Usng old Phase Equaton of tate Table 1. Mole fratons for two synthet mxtures [7] Component Mxture C Mxture B n-c 1.5876.511 n-c 18.513.819 n-c 19.486.694 n-c.463.59 n-c 1.44.56 n-c.418.433 n-c 3.397.373 n-c 4.378.319 n-c 5.359.74 n-c 6.34.36 n-c 7.37. n-c 8.176 n-c 9.148 n-c 3.17 Table. Mole fratons for a synthet mxture [14] Component Bm 13 Component Bm 13 C 1 8.1 C 34.61 C 18 7.9 C 35.53 C 19 6.9 C 36.45 5. C Table 3. Mole fratons of heavy ol fratons [5] Ol 5 Peusoomponent Mole perent MW P-C 1+ 4.467 167. N-C 1+ 6.487 16. A-C 1+ 15.16 16. P-C 15+.996 37. N-C 15+ 3.867 33. A-C 15+ 8.9664 33. P-C + 1.546 37. N-C +.1514 3. A-C + 5.199 3. P-C 5+.7856 375. N-C 5+ 1.389 37. A-C 5+ 3.49 37. P-C 3+.358 449. N-C 3+ 1.4348 44. A-C 3+ 1.4348 44. P-C 35+.1377 511. N-C 35+ 1.5694 51. A-C 35+.174 51. P-C 4+.648 59. N-C 4+ 1.1964 587. A-C 4+.491 587. P-C 46+.59 713. N-CP1 +.3143 74. A-CP1 + 1.885 74. N-CP +.57 91. A-CP + 1.3396 91.
M. Vafae eft et al. Table 4. Mole fratons of heavy ol fratons [5] Ol 6 Pseudoomponent Mole perent MW P-CP1 3.59 157. N-CP1 4.771 157. A-CP1 4.771 157. P-CP.7858 1. N-CP 4.5495 1. A-CP 4.5495 1. P-CP3 1.855 5. N-CP3.989 5. A-CP3 4.4744 5. P-CP4 1.38 3. N-CP4.918 3. A-CP4 4.357 3. P-CP5.3674 563. N-CP5.5116 563. A-CP5 5.1937 563. P-CP6.581 654. N-CP6 1.319 654. A-CP6 1.677 654. P-CP7.736 666. N-CP7 1.99 666. A-CP7.8634 666. N-CP8.8611 744. A-CP8.177 744. The parameters of the equaton of state (l s, m s, n s ) are alulated as dsussed n the prevous seton. The optmal values for the mentoned parameters for the ol samples are n tables 5-9. The ntal guess for the system of equlbrum equatons s gven from the results of a two-phase flash alulaton. Then, the dogleg method [15] s appled to hek the onvergene rtera,.e. the rght-hand sde expresson n equaton 19 should be less than 1e-7. If the rteron s not met, the program wll shft to smplex algorthm that uses the results of the prevous step as the ntal ponts. There s a normalzng step, whh flters the nomng physally unaeptable data. The physal propertes data are derved from databooks or are estmated from proposed methods n the lterature. Table 5. Parameters of EO for mxture B Component l m n n-c 1-9.44 13.39 -.1876 n-c 18 7.1395 -.44956-4.3144 n-c 19-8.993 15.535 -.35543 n-c 7.9884 -.5493-4.596 n-c 1-9.735 16.518 -.1681 n-c 8.616 -.57798-4.59 n-c 3-9.518 17.335 -.37988 n-c 4 9.41 -.64356-4.3 n-c 5-9.765 18.17.6465 n-c 6 1.11 -.748-4.9 n-c 7-9.7534 18.749.19831 n-c 8 1.816 -.76167-4.198
Predton of old Paraffn Preptaton Usng old Phase Equaton of tate Table 5. Parameters of EO for mxture B (ontd.) n-c 9-9.995 19.34.441 n-c 3 11.444 -.81451-4.1968 Table 6. Parameters of EO for mxture C Component l m n n-c 1-9.44 13.39 -.1876 n-c 18 7.1395 -.44956-4.3144 n-c 19-8.993 15.535 -.35543 n-c 7.9884 -.5493-4.596 n-c 1-9.735 16.518 -.1681 n-c 8.616 -.57798-4.59 n-c 3-9.518 17.335 -.3799 n-c 4 9.41 -.64356-4.3 n-c 5-9.765 18.17.6465 n-c 6 1.11 -.748-4.9 n-c 7-9.7534 18.749.19831 Table 7. Parameters of EO for bm13 Component l m n C 1-1.833 14.648 -.581 C 18 7.355 -.47341-4.847 C 19 6.4639 -.9946-4.563 C 8.3 -.53895-4.5 C 34 1.96-1.438-4.151 C 35-16.46 8.979.11 13.63-1.1167-4.143 C 36 Table 8. Parameters of EO for Ol 5 Pseudoomponent l m n P-C 1+ -8.61 1.786 -.1739 N-C 1+ 3.15 -.8398-5.1386 A-C 1+ 3.67 -.194-4.9696 P-C 15+ 5.9567 -.761-4.538 N-C 15+ -5.7139 9.9654 -.17938 A-C 15+ -6.577 9.9896 -.1857 P-C + 7.49 -.398-4.47 N-C + -6.3593 11.679 -.15433 A-C + -5.837 1.416 -.1188 P-C 5+ 8.9575 -.5176-4.374 N-C 5+ -7.447 13.844 -.153 A-C 5+ -6.94 11.417 -.147 P-C 3+ 1.561 -.678-4.3345 N-C 3+ -8.654 16.1 -.9867 A-C 3+ -6.588 1.693 -.1949 P-C 35+ 11.91 -.73866-4.353 N-C 35+ -1.19 18.873 -.79 A-C 35+ -7.483 14.39 -.187 P-C 4+ 13.681 -.8873-4.681 N-C 4+ -11.777 1.77 -.559 A-C 4+ -8.138 15.977 -.19855 P-C 46+ 16.49-1.1337-4.9 N-CP1 + -14.97 7.17 -.1868 A-CP1 + -9.4777 19.316 -.1736 N-CP + -19.61 34.39.9795 A-CP + -13.495 3.969 -.113
M. Vafae eft et al. Table 9. Parameters of EO for Ol 6 Pseudoomponent l m P-CP1-7.671 11.668 N-CP1-4.749 7.8714 A-CP1 3.1885 -.13346 P-CP -7.9844 13.48 N-CP 3.8369 -.164 A-CP 3.6483 -.11553 P-CP3 6.4993 -.37657 N-CP3-4.3969 9.15 A-CP3-4.3753 8.6468 P-CP4 7.6419 -.4964 N-CP4-4.6555 1.188 A-CP4-4.311 8.8143 P-CP5 13.934-1.69 N-CP5-7.8951 17.94 A-CP5-4.135 1.168 P-CP6 16.199-1.4885 N-CP6-9.3345.954 A-CP6-4.411 13.697 P-CP7 16.5-1.569 N-CP7-9.589 1.355 A-CP7-4.4491 13.94 N-CP8-1.87 3.979 A-CP8-4.6938 15.65 n -.15 -.1643-4.751 -.4753-4.84-4.83-4.186 -.9495 -.11886-4.155 -.64131 -.137-3.8195.63654 -.78385-3.754.9889 -.5575-3.7459.9417 -.5763.11 -.33513 α s ndependent of pressure and t an be used for the sold volume predton [11]. The error funton for lghter omponents, lke C 1, are greater than heavy fratons. The reason s that the alulatons are done n a temperature far from the meltng pont. Preptated Wax Wt% 1 9 8 7 Ths Work 6 Expt. 5 4 3 1 7 9 31 33 35 Temperature, K Fgure 1. Expermental and alulated amount of preptated wax for Ol 5
Predton of old Paraffn Preptaton Usng old Phase Equaton of tate Preptated Wax Wt% 16 14 1 1 8 6 4 Ths Work Expt. 7 8 9 3 31 3 Temperature (K) Fgure. Expermental and alulated amount of preptated wax for Ol 6 3 5 Preptated Wax Wt% 15 1 5 Ths Work Expt. 5 6 7 8 9 3 31 3 Temperature (K) Fgure 3. Expermental and alulated amount of preptated wax for bm13 6 Preptated Wax Wt% 5 4 3 1 Ths Work Expt. 65 75 85 95 35 315 Temperature (K) Fgure 4. Expermental and alulated amount of preptated wax for mxture B
M. Vafae eft et al. 6 mxture C Preptated Wax Wt% 5 4 3 1 Ths Work Expt. 65 75 85 95 35 Temperature (K) Fgure 5. Expermental and alulated amount of preptated wax for mxture C 5. Conluson Complex behavor of sold phase n ol mxture and wde range of ts applaton n sold preptaton and deposton petroleum fluds (wax, asphaltene, ) need to be modeled va applable and effent methods. Here applaton of a sold EO for desrpton of sold phase was tested for wax preptaton n petroleum mxtures. In ths work, TT sold equaton of state s used for desrbng wax preptaton phenomena n some synthet and real ol mxtures. Ths sold equaton of state s based on an alpha funton. Usng thermodynam yle for pure sold fugaty from pure lqud fugaty, the TT parameters were tuned before ts applaton for wax preptaton predton. The multsold phase approah s used for determnaton of the nature and number of sold phases. As t an be seen n the prevous setons, the obtaned results n ths work are n good agreements wth the expermental data. Referenes [1] Won K. V., Flud Phase Equlbra, 53, 377 (1989). [] Won K. V., Flud Phase Equlbra, 3, 65 (1986). [3] Dardon J.-L., Pauly J., Coutnho J. A. P., and Montel F., Energy and Fuels, 15, 73-735 (1). [4] Coutnho J. A. P. and Dardon J.-L., Energy and Fuels, 15, 1454 (1). [5] Vafae-eft M., Mousav-Dehghan. A., and Mohammad-Zadeh Bahar M., Flud Phase Equlbra, 173, 65 (). [6] Nhta D. V., Gousl L., and Froozabad A., PE 74686 (1). [7] Feyz F. and Dalrsefat R., Fuel (7), do:1.116/j.fuel.6.11.34, Artle n Press. [8] Lra-Galeana C., Froozabad A., and Prausntz J. M., AIChE, 4, 39 (1996).
Predton of old Paraffn Preptaton Usng old Phase Equaton of tate [9] Coutnho J. A. P., Edmonds B., Moorwood T., zzepansk R., and Zhang X., PE 7834 (). [1] Yokozek A., Internatonal Journal of Thermophyss, 4, 589 (3). [11] Twu C. H., Tassone V., and m W. D., AIChEJ, 49, 957 (3). [1] Ness H. C. V. and Abbott M. H., Classal Thermodynams of Non-Eletrolyte olutons wth Applaton tp Phase Equlbra. New York: M-Graw Hll (198). [13] Prausntz J., Lhtenthaler M., R. N., and de Azevedo E. G., Moleular Thermodynams of Flud-Phase Equlbra. Upper addle Rver, New Jersey: Prnte Hall (1999). [14] Esobar-Remolna J. C. M., Flud Phase Equlbra, 4, 197 (6). [15] Powell, M. J. D., Numeral Methods for Nonlnear Algebra Equatons, P. Rabnowtz (197)