Krg

In [1]:

#krow

#krow

In [2]:

from krpy import Corey, Krog, sg_denormalize, sg_normalize
import pandas as pd
import numpy as np

from krpy import Corey, Krog, sg_denormalize, sg_normalize import pandas as pd import numpy as np

In [3]:

cor = Corey(no=2, ng=1.7, krg_end=0.2,kro_end=0.98, pcg_end=0,npc=1.0, pd=2.0)
cor

cor = Corey(no=2, ng=1.7, krg_end=0.2,kro_end=0.98, pcg_end=0,npc=1.0, pd=2.0) cor

Out[3]:

Corey(nw=2.0, no=2.0, ng=1.7, npcwo=2.0, npcog=2.0, krw_end=1.0, kro_end=0.98, krg_end=0.2, pco_end=0.0, pcg_end=0.0, pdwo=0.0, pdog=0.0)

In [4]:

krog = Krog.from_corey(cor, swir=0.22, sor=0.1,sgc=0.1, n=10)

krog = Krog.from_corey(cor, swir=0.22, sor=0.1,sgc=0.1, n=10)

/Users/scuervo91/Documents/dev/krpy/krpy/kr.py:209: RuntimeWarning: divide by zero encountered in true_divide
  kr_ratio = krg / kro

In [5]:

df = krog.df()
print(df)

df = krog.df() print(df)

                 krg       kro  pcog       sgn   kr_ratio
saturation                                               
0.100000    0.000000  0.980000   0.0  0.000000   0.000000
0.164444    0.004773  0.774321   0.0  0.111111   0.006164
0.228889    0.015508  0.592840   0.0  0.222222   0.026160
0.293333    0.030898  0.435556   0.0  0.333333   0.070938
0.357778    0.050387  0.302469   0.0  0.444444   0.166586
0.422222    0.073632  0.193580   0.0  0.555556   0.380370
0.486667    0.100386  0.108889   0.0  0.666667   0.921916
0.551111    0.130462  0.048395   0.0  0.777778   2.695773
0.615556    0.163708  0.012099   0.0  0.888889  13.530992
0.680000    0.200000  0.000000   0.0  1.000000        inf

In [6]:

print(krog.to_ecl())

print(krog.to_ecl())

SGOF
0.100000 0.000000 0.980000 0.0
0.164444 0.004773 0.774321 0.0
0.228889 0.015508 0.592840 0.0
0.293333 0.030898 0.435556 0.0
0.357778 0.050387 0.302469 0.0
0.422222 0.073632 0.193580 0.0
0.486667 0.100386 0.108889 0.0
0.551111 0.130462 0.048395 0.0
0.615556 0.163708 0.012099 0.0
0.680000 0.200000 0.000000 0.0/

In [7]:

krog.interpolate(0.5)

krog.interpolate(0.5)

Out[7]:

	krg	kro	pcog	sgn	kr_ratio
saturation
0.5	0.106609	0.096373	0.0	0.689655	1.288921

In [8]:

krog.plot(norm=True)

krog.plot(norm=True)

Out[8]:

[<AxesSubplot:xlabel='Gas Saturation []', ylabel='Kr []'>]