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    
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

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+
 FIT SEPOIL-SCN C11 MW      
END