Algoritmos Petrofisica

GT = (216 - 80) / (9997 / 100)
Tf[] = (GT * (DEPTH[] / 100)) + 80
If(RHOB[] < 2.1 and DT[] > 100 and NPHI[] > .51) Then CARB[] = 1 Else CARB[] = 0
Rmf[] = 2.15 * ((75 + 6.77) / (Tf[]+ 6.77))
Rmfe[] = Rmf[] * 0.85
Kc[] = 61 + (0.133 * Tf[])
RWE[] = Rmfe[] / (10 ^ (SP_ED[] / -Kc[]))
; Density Porosity
PHID[] = (RhoM - RHOB[]) / (RhoM - RhoF)
; DIFE[] = NPHI[] - PHID[]Vshl[] = min(1, max(0, (GR[]-GRcln) / (GRshl-GRcln) ))
RT[] = LLD_ED[]
Rxo[] = MSFL_ED[]
PHIDC[] = PHID[] - (Vshl[] * PHIDSH)
NPHIC[] = NPHI[] - (Vshl[] * HIshl)
PHIE[] = (PHIDC[] +NPHIC[]) / 2
If(PHIE[] < 0) Then PHIE[]=0 Else PHIE[]=PHIE[]
If(Vshl[] < .75) Then PHIE[] = PHIE[] Else PHIE[] = 0
; Indonesian SwI; Requires, PhiE[], Vshl[] & RT[]
SwI[] = sqrt(1/RT[]) /( Vshl[]^(1- 0.5*Vshl[]) /sqrt(Rshl) + sqrt( PHIE[]^m / (a*Rw)) )
If(CARB[] = 0) Then PHIE[] = PHIE[] Else PHIE[] = 0
If(CARB[] = 0) Then SwI[] = SwI[] Else SwI[] = 1
If (SwI[] <= 0.75)Then SwI[] = SwI[] Else SwI[] = 1
If (PHIE[] >= .06) Then SwI[] = SwI[] Else SwI[] = 1
If (Vshl[] <= .75) Then SwI[] = SwI[] Else SwI[] = 1

; Indonesian SxoI; Requires, PhiE[], Vshl[]& Rxo[]
SxoI[] = sqrt(1/Rxo[]) / ( Vshl[]^(1- 0.5*Vshl[]) /sqrt(Rshl) + sqrt( PHIE[]^m / (a*Rmf[])) )
If (SxoI[] <= 1) Then SxoI[] = SxoI[] Else SxoI[] = 1
If (PHIE[] >= .06) ThenSxoI[] = SxoI[] Else SxoI[] = 1
If (Vshl[] <= .75) Then SxoI[] = SxoI[] Else SxoI[] = 1
If(SxoI[] >= SwI[]) Then SxoI[] = SxoI[] Else SxoI[] = SwI[]
BVW[] = PHIE[] * SwI[]
BVXO[] =PHIE[] * SxoI[]

; “Saturacion de PetroleoMovil”
So_Mov[] = SxoI[] - SwI[]
BVW[] = PHIE[] * SwI[]
BVXO[] = PHIE[] * SxoI[]

; Timur Permeability
KTim[] = (100 * PHIE[]^2.25/SwI[])^2If(Vshl[] < 0.5 and PHIE[] > 0.08) Then AN[] = 1 Else AN[] = 0
; “Calculo de NETPAY”
If(Vshl[] < 0.5 and PHIE[] > 0.08 and SwI[] < 0.5) Then NETPAY[] = 0.5 Else NETPAY[] = 0