def calculation(self, parameters):
par = {}
n = len(parameters)
n1 = len(self._phi)
n2 = len(self._co2)
n3 = len(self._p)
for i in range(n):
par[parameters[i]] = np.zeros(((n3, n2, n1)))
vpBrineCo2 = ep.vp(self._ksat_l, self._gdry_l, self._rhosat_l,
self._phi, self._co2, self._p)
vsBrineCo2 = ep.vs(self._gdry_l, self._rhosat_l, self._phi, self._co2,
self._p)
ipBrineCo2 = ep.Ip(vpBrineCo2, self._rhosat_l, self._phi, self._co2,
self._p)
isBrineCo2 = ep.Ip(vsBrineCo2, self._rhosat_l, self._phi, self._co2,
self._p)
rtBrineCo2 = ep.vp_vs(vpBrineCo2, vsBrineCo2, self._phi, self._co2,
self._p)
vpBrineOil = ep.vp(self._ksat_loil, self._gdry_l, self._rhosat_loil,
self._phi, self._co2, self._p)
vsBrineOil = ep.vs(self._gdry_l, self._rhosat_loil, self._phi,
self._co2, self._p)
ipBrineOil = ep.Ip(vpBrineOil, self._rhosat_loil, self._phi, self._co2,
self._p)
isBrineOil = ep.Ip(vsBrineOil, self._rhosat_loil, self._phi, self._co2,
self._p)
rtBrineOil = ep.vp_vs(vpBrineOil, vsBrineOil, self._phi, self._co2,
self._p)
par['vpBrineCo2'] = vpBrineCo2
par['vsBrineCo2'] = vsBrineCo2
par['ipBrineCo2'] = ipBrineCo2
par['isBrineCo2'] = isBrineCo2
par['rtBrineCo2'] = rtBrineCo2
par['vpBrineOil'] = vpBrineOil
par['vsBrineOil'] = vsBrineOil
par['ipBrineOil'] = ipBrineOil
par['isBrineOil'] = isBrineOil
par['rtBrineOil'] = rtBrineOil
return par
def calculation (self,parameters):
par = {}
n = len(parameters)
n1 = len(self._phi)
n2 = len(self._co2)
n3 = len(self._p)
for i in range(n):
par[parameters[i]] = np.zeros(((n3,n2,n1)))
vpBrineCo2 = ep.vp(self._ksat_l,self._gdry_l,self._rhosat_l,self._phi,self._co2,self._p)
vsBrineCo2 = ep.vs(self._gdry_l,self._rhosat_l,self._phi,self._co2,self._p)
ipBrineCo2 = ep.Ip(vpBrineCo2,self._rhosat_l,self._phi,self._co2,self._p)
isBrineCo2 = ep.Ip(vsBrineCo2,self._rhosat_l,self._phi,self._co2,self._p)
rtBrineCo2 = ep.vp_vs(vpBrineCo2,vsBrineCo2,self._phi,self._co2,self._p)
vpBrineOil = ep.vp(self._ksat_loil,self._gdry_l,self._rhosat_loil,self._phi,self._co2,self._p)
vsBrineOil = ep.vs(self._gdry_l,self._rhosat_loil,self._phi,self._co2,self._p)
ipBrineOil = ep.Ip(vpBrineOil,self._rhosat_loil,self._phi,self._co2,self._p)
isBrineOil = ep.Ip(vsBrineOil,self._rhosat_loil,self._phi,self._co2,self._p)
rtBrineOil = ep.vp_vs(vpBrineOil,vsBrineOil,self._phi,self._co2,self._p)
par['vpBrineCo2'] = vpBrineCo2
par['vsBrineCo2'] = vsBrineCo2
par['ipBrineCo2'] = ipBrineCo2
par['isBrineCo2'] = isBrineCo2
par['rtBrineCo2'] = rtBrineCo2
par['vpBrineOil'] = vpBrineOil
par['vsBrineOil'] = vsBrineOil
par['ipBrineOil'] = ipBrineOil
par['isBrineOil'] = isBrineOil
par['rtBrineOil'] = rtBrineOil
return par
Vp_Vs_sat_l_oil = np.zeros(((NP,nCO2,Np)))
Ip_sat_l_oil = np.zeros(((NP,nCO2,Np)))
Is_sat_l_oil = np.zeros(((NP,nCO2,Np)))
rho_sat_l = np.zeros(((NP,nCO2,Np)))
rho_sat_l_oil = np.zeros(((NP,nCO2,Np)))
for i in range(NP):
for k in range(nCO2):
Vssat_l1= ep.vs(Gdry_constant[i,:],rhosat_l_contact[i,:,k],phi)
Vpsat_l1= ep.vp(Ksat_l_contact[i,:,k],Gdry_constant[i,:],rhosat_l_contact[i,:,k],phi)
Vp_Vs_sat_l1= ep.vp_vs(Vpsat_l1,Vssat_l1,phi)
Ip_sat_l1=ep.Ip(Vpsat_l1,rhosat_l_contact[i,:,k],phi)
Is_sat_l1=ep.Ip(Vssat_l1,rhosat_l_contact[i,:,k],phi)
Vssat_l3= ep.vs(Gdry_constant[i,:],rhosat_l_contact_oil[i,:,k],phi)
Vpsat_l3= ep.vp(Ksat_l_contact_oil[i,:,k],Gdry_constant[i,:],rhosat_l_contact_oil[i,:,k],phi)
Vp_Vs_sat_l3= ep.vp_vs(Vpsat_l3,Vssat_l3,phi)
Ip_sat_l3=ep.Ip(Vpsat_l3,rhosat_l_contact_oil[i,:,k],phi)
Is_sat_l3=ep.Ip(Vssat_l3,rhosat_l_contact_oil[i,:,k],phi)
'''
Vssat_l_oil = np.zeros(((NP, nCO2, Np)))
Vp_Vs_sat_l_oil = np.zeros(((NP, nCO2, Np)))
Ip_sat_l_oil = np.zeros(((NP, nCO2, Np)))
Is_sat_l_oil = np.zeros(((NP, nCO2, Np)))
rho_sat_l = np.zeros(((NP, nCO2, Np)))
rho_sat_l_oil = np.zeros(((NP, nCO2, Np)))
for i in range(NP):
for k in range(nCO2):
Vssat_l1 = ep.vs(Gdry_constant[i, :], rhosat_l_contact[i, :, k], phi)
Vpsat_l1 = ep.vp(Ksat_l_contact[i, :, k], Gdry_constant[i, :],
rhosat_l_contact[i, :, k], phi)
Vp_Vs_sat_l1 = ep.vp_vs(Vpsat_l1, Vssat_l1, phi)
Ip_sat_l1 = ep.Ip(Vpsat_l1, rhosat_l_contact[i, :, k], phi)
Is_sat_l1 = ep.Ip(Vssat_l1, rhosat_l_contact[i, :, k], phi)
Vssat_l3 = ep.vs(Gdry_constant[i, :], rhosat_l_contact_oil[i, :, k],
phi)
Vpsat_l3 = ep.vp(Ksat_l_contact_oil[i, :, k], Gdry_constant[i, :],
rhosat_l_contact_oil[i, :, k], phi)
Vp_Vs_sat_l3 = ep.vp_vs(Vpsat_l3, Vssat_l3, phi)
Ip_sat_l3 = ep.Ip(Vpsat_l3, rhosat_l_contact_oil[i, :, k], phi)
Is_sat_l3 = ep.Ip(Vssat_l3, rhosat_l_contact_oil[i, :, k], phi)
'''
Vssat_l1= ep.vs(Gdry_l[i,:],rhosat_l_contact[i,:,k],phi)