TITLE "JK1 EOS MODEL" TITLE "ADVANCED PVT & EOS FLUID CHARACTERIZATION COURSE" TITLE "VIENNA, AUSTRIA, NOVEMBER 24-29, 2013 TITLE "Your Name" STAB ; initialize K-values with stability test TEST1 ; test single-phase solution for 2nd phase TEST2 ; test two-phase solution for 3rd phase -- Define a fluid characterization corresponding to the -- measured compositions, to simplify their input. CHAR LAB-INPUT EOS SRK NAME MW H2S N2 CO2 C1 C2 C3 I-C4 N-C4 NEO-C5 ISO-C5 N-C5 22DM-C4 C-C5 23DM-C4 2M-C5 3M-C5 C6* 86 N-C6 MC-C5 Benzene C-C6 2M-C6 3M-C6 224TM-C5 C7* 99 N-C7 MC-C6 Toluene C8* 110 N-C8 E-Benzene M-Xylene O-Xylene C9* 126 N-C9 C10* 141 N-C10 C11 147 C12 161 C13 175 C14 190 C15 206 C16 222 C17 237 C18 251 C19 263 C20 275 C21 291 C22 305 C23 318 C24 331 C25 345 C26 359 C27 374 C28 388 C29 402 C30 416 C31+ 485 END -- Because lab actually measures mass fractions, these should be input to PVT program. -- If mass fractions not reported, then mole fractions can be input together with MWs -- to calculate mass fractions. Once calcualted, you can "conserve" them. CONSERVE MASS MIX SEPGAS MOLES 0.000 0.484 2.037 82.531 7.266 5.155 0.779 1.047 0.000 0.228 0.187 0.004 0.022 0.000 0.037 0.019 0.000 0.039 0.025 0.005 0.027 0.008 0.006 0.000 0.014 0.012 0.025 0.007 0.017 0.003 0.001 0.003 0.001 0.006 0.001 0.003 0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 MIX SEPOIL MOLES 0.000 0.072 0.847 17.472 6.165 11.852 3.689 6.843 0.066 3.026 3.169 0.087 0.000 0.473 1.125 0.679 0.000 1.614 1.252 0.264 1.801 0.477 0.481 0.000 1.256 1.298 3.052 1.035 2.931 1.175 0.145 0.932 0.190 2.304 0.854 2.622 0.724 2.486 2.075 2.037 1.805 1.884 1.318 0.998 1.025 0.806 0.669 0.595 0.518 0.453 0.416 0.375 0.341 0.315 0.290 0.272 0.221 1.129 MIX WELLSTREAM MOLES 0.000 0.393 1.776 68.250 7.024 6.625 1.418 2.319 0.014 0.842 0.842 0.022 0.017 0.104 0.276 0.164 0.000 0.385 0.294 0.062 0.416 0.111 0.110 0.000 0.287 0.294 0.689 0.233 0.657 0.260 0.033 0.207 0.042 0.510 0.188 0.578 0.160 0.546 0.455 0.447 0.396 0.413 0.289 0.219 0.225 0.177 0.147 0.131 0.114 0.100 0.091 0.082 0.075 0.069 0.064 0.060 0.049 0.248 -- Define lumping rules for single carbon numbers 6 through 10. -- These rules will be used for composition conversions. LUMP C6 22DM-C4, C-C5, 23DM-C4, 2M-C5, 3M-C5, C6*, N-C6 LUMP C7 MC-C5, Benzene, C-C6, 2M-C6, 3M-C6, 224TM-C5, C7*, N-C7 LUMP C8 MC-C6, Toluene, C8*, N-C8 LUMP C9 E-Benzene, M-Xylene, O-Xylene, C9*, N-C9 LUMP C10 C10*, N-C10 LUMP I-C5 NEO-C5, ISO-C5 -- Soreide correlation used to generate relation between specific gravity and MW of -- undefined hydrocarbons (C7, C8, ..., CN+). -- Command below prints the Soreide equation and default parameters (changeable). SOREIDE REPLACE SOREIDE FACTOR BY 0.298 NOTE: "SETUP SINGLE CARBON NUMBER (SCN) EOS MODEL" NOTE: "Use GAMMA model to get MWS of C11+ for characterization, also C31+" GAMMA SPLIT C11 SHAPE 1 ; exponential BOUND 143 ; (14*11-4)-7 (or use LMW of C11 found in CHAR below) AVERAGE 240.3 ; use lab C11+ MW (or best-fit C11+ MW) END MIX FEED WELLSTREAM ; used for averageing of lumped properties. CHAR EOSSCN BASE LAB-INPUT ; use MWs from the LAB-INPUT characterization as basis for calculating the MWs of the lumped SCN in the EOSSCN characterization. EOS SRK NAME LMW MW TC(R) PC(PSIA) AF VSHIFT SG TB(R) N2 CO2 C1 C2 C3 I-C4 N-C4 I-C5 N-C5 C6 C7 C8 C9 C10 C11 C12 C13 C14 C15 C16 C17 C18 C19 C20 C21 C22 C23 C24 C25 C26 C27 C28 C29 C30 C31+ END BINARY N2 CO2 CO2 0.000 C1 0.020 0.120 C2 0.060 0.150 REPLACE BINARY OF N2 WITH C3 TO C31+ BY 0.08 REPLACE BINARY OF CO2 WITH C3 TO C31+ BY 0.15 REPLACE MW OF C11 TO C31+ by -99 ; replace the MWs of the C11+ fraction from the LAB-INPUT chararacterization by a non-physical value. CORRELATE ; calculate the MWs of C11+ components using the Gamma distribution model defined. LUMP C7+ C7 25*1 -- Convert the previously specified, measured compositions (LAB-INPUT characterization), -- to the "EOSSCN" characterization : GAMMA SPLIT NONE ; don't use gamma splitting MIX WELLSTREAM-SCN WELLSTREAM 1 MOLE ; MIX SEPGAS-SCN SEPGAS 1 MOLE ; MIX SEPOIL-SCN SEPOIL 1 MOLE ; GAMMA SPLIT C11 FIT SEPOIL-SCN C11 AMOUNTS IGNORE C31+ END Note: "Saturation pressure calculation of Wellstream-SCN" MIX FEED WELLSTREAM-SCN 1 MOLE PRES PSIA TEMP 263 F PSAT Note: "Saturation pressure calculation of SepOil-SCN" MIX FEED SEPOIL-SCN 1 MOLE PRES PSIA TEMP 94 F PSAT MIX FEED SEPOIL-SCN 1 MOLE TEMP 94 F PSAT TEMP 77 F PSAT DLE ID "SEP TEST OF SEP OIL" STAGE TEMP PRES LDEN F PSIA G/CC 1 94 665 0.7071 2 80 91 3 120 45 4 75 14.7 0.7888 5 60 14.7 0.7961 END MIX FEED WELLSTREAM-SCN 1 MOLE TEMP 263 F PSAT DLE ID "SEP TEST OF WELLSTREAM-SCN MIXTURE" BASIS 14.60 LBMOL STAGE TEMP PRES LDEN GOR MREM LVOL GSG CUMOGR Y-N2 Y-CO2 Y-C1 Y-C2 Y-C3 Y-I-C4 Y-N-C4 Y-I-C5 Y-N-C5 Y-C6 Y-C7+ X-N2 X-CO2 X-C1 X-C2 X-C3 X-I-C4 X-N-C4 X-I-C5 X-N-C5 X-C6 X-C7+ F PSIA G/CC SCF/BBL SCF BBL BBL/MMSCF % % % % % % % % % % % % % % % % % % % % % % 1 94 665 0.7071 3436 4320 1.2573 0.705 0.484 2.037 82.531 7.266 5.155 0.779 1.047 0.228 0.187 0.121 0.165 0.072 0.847 17.472 6.165 11.852 3.689 6.843 3.092 3.169 3.978 42.822 2 80 91 4583 0.926 3 120 45 4630 1.264 4 75 14.7 0.788 4734 1.0079 1.758 5 60 14.7 0.7961 1.0 END MIX FEED WELLSTREAM-SCN 1 MOLE TEMP 264 F PRES 5535 PSIG CCE ID "LAB CCE TEST" SATP 5535 PSIG PRES RVOL DEN LRVOL GZ PSIG G/CC % 11000 0.82736 0.47758 1.77731 10440 0.83715 0.47199 1.70691 10000 0.84486 0.46769 1.65013 9000 0.86686 0.45581 1.52404 8000 0.89274 0.4426 1.39542 7000 0.92595 0.42673 1.26675 6500 0.94639 0.41751 1.20243 6000 0.97085 0.40699 1.13883 5535 1.00000 0.39513 0.00 1.08234 5178 1.03101 4.26 4774 1.07253 15.84 4280 1.13887 21.07 3898 1.20738 22.84 3456 1.3127 23.83 3070 1.43810 24.61 2446 1.74876 24.16 1827 2.30634 22.97 1383 3.05722 21.93 938 4.56767 20.78 END MIX FEED WELLSTREAM-SCN 1 MOLE TEMP 264 F PRES 5535 PSIG CVD ID "LAB CVD TEST WITHOUT BLEEDING EMULATION" TEMP PRES RMREM LRVOL GZ Y-C1 Y-C7+ LAPI F PSIG % % % % 263 5535 0 0 1.0823 68.250 9.529 4500 7.976 17.000 0.9202 71.913 6.200 3500 21.980 22.842 0.8515 74.339 3.900 2500 39.965 23.832 0.8369 75.538 2.600 1700 56.259 22.309 0.8568 76.049 2.000 900 73.951 19.063 0.9069 74.723 2.000 263 500 263 400 263 300 263 250 263 200 263 150 263 100 263 50 263 25 263 0 14.967 60 0 40.78 END