Exemplo n.º 1
0
    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
Exemplo n.º 2
0
  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
Exemplo n.º 3
0
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) 
        
         

        '''
Exemplo n.º 4
0
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)