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)