Beispiel #1
1
import sgems
import math

seed = 14071789
num_realizations = 10000
nrow, ncol = 158, 93
delr, delc = 1000, 1000

offset = (503800.0, 2778800.0)
realization_path = "D:/Users/jwhite/Projects/Broward/Geostats/SGEMS/realizations/"
realization_prefix = "layer1_thk_omni_"

try:
    sgems.execute("DeleteObjects Layer_Thk::1000_grid")
except:
    pass
load_command = "LoadObjectFromFile  D:/Users/jwhite/Projects/Broward/Geostats/SGEMS/Layer1_cdf.sgems::s-gems"
sgems.execute(load_command)

newGrid_command = (
    "NewCartesianGrid  1000_grid::"
    + str(ncol)
    + "::"
    + str(nrow)
    + "::1::"
    + str(delc)
    + "::"
    + str(delr)
    + "::1.0::"
    + str(offset[0])
    + "::"
Beispiel #2
0
    def execute(self):
	print self.params
	# ESTABELECA OS PARAMETROS INICIAIS
	propertie = self.params['prop']['property']
	grid = self.params['prop']['grid']
	
	measured = float(self.params['Measured']['value'])
	Indicated = float(self.params['Indicated']['value'])
	Inferred = float(self.params['Inferred']['value'])
	dx = float(self.params['dx']['value'])
	dy = float(self.params['dy']['value'])
	dz = float(self.params['dz']['value'])
	
	grid = self.params['prop']['grid'] 
	x = sgems.get_property(grid, "_X_")
	y = sgems.get_property(grid, "_Y_")
	z = sgems.get_property(grid, "_Z_")
	distance = []
	
	

	numx = int(round((max(x) - min(x))/dx,0)) + 1
	numy = int(round((max(y) - min(y))/dy,0)) + 1
	numz = int(round((max(z) - min(z))/dz,0)) + 1

		

	cmd="NewCartesianGrid rescat::"+str(numx)+"::"+str(numy)+"::"+str(numz)+"::"+str(dx)+"::"+str(dy)+"::"+str(dz)+"::"+str(min(x))+"::"+str(min(y))+"::"+str(min(z))+"::0"
        print cmd
	sgems.execute(cmd)

	xh = sgems.get_property("rescat","_X_")
	yh = sgems.get_property("rescat","_Y_")
	zh = sgems.get_property("rescat","_Z_")	

	

	
	v = []
	for i, j, k in zip(xh,yh,zh):
		valores = 0
		menor_valor = 1000
		for xind, yind, zind in zip(x,y,z):
			distance = (math.sqrt((xind-i)**2+(yind-j)**2+(zind-k)**2))
			if (distance < measured):
				valores = 0
			elif( distance < Indicated and distance > measured):
				valores = 1
			elif(distance < Inferred and distance > Indicated):
				valores = 2
			else:
				valores = 3
			if (valores < menor_valor):
				menor_valor = valores		
		v.append(menor_valor)
	print (v)
	sgems.set_property("rescat","R",v)
	
	
    	return True 
Beispiel #3
0
    def execute(self):
      
        #aqui vai o codigo
        #getting variables
        point_set_name = self.params['gridselectorbasic']['value']
        new_grid_name = self.params['lineEdit']['value']
        buffer = float(self.params['doubleSpinBox_4']['value'])
        sx, sy, sz = float(self.params['doubleSpinBox']['value']), float(self.params['doubleSpinBox_2']['value']), float(self.params['doubleSpinBox_3']['value'])

        x, y, z = np.array(sgems.get_X(point_set_name)), np.array(sgems.get_Y(point_set_name)), np.array(sgems.get_Z(point_set_name))
        grid_dic = autogrid(x, y, z, sx, sy, sz, buffer)

        sgems.execute('NewCartesianGrid  {}::{}::{}::{}::{}::{}::{}::{}::{}::{}::0,00'.format(new_grid_name, grid_dic['nx'], grid_dic['ny'], grid_dic['nz'], grid_dic['sx'], grid_dic['sy'], grid_dic['sz'], grid_dic['ox'], grid_dic['oy'], grid_dic['oz']))

        return True
Beispiel #4
0
import sgems

# 1. The output path of where to store the generated facies models
path = '/home/lewisli/Desktop/WCA/Generated Models/'

# Load project files
sgems.execute("LoadProject \
/home/lewisli/Desktop/WCA/InitialData/InitialData.prj/")

sgems.execute(
    "ClearPropertyValueIf ProbacubeLess::ProbacubeLess_PCAChanProb::-996660000::-10000"
)
sgems.execute(
    "ClearPropertyValueIf ProbacubeLess::ProbacubeLess_PCALevProb::-996660000::-10000"
)
sgems.execute(
    "ClearPropertyValueIf ProbacubeLess::ProbacubeLess_PCAPoorLevProb::-996660000::-10000"
)
sgems.execute(
    "ClearPropertyValueIf ProbacubeLess::ProbacubeLess_PCAShaleProb::-996660000::-10000"
)

sgems.execute(
    "ClearPropertyValueIf ProbacubeMore::ProbacubeMore_PCAChanProb::-996660000::-10000"
)
sgems.execute(
    "ClearPropertyValueIf ProbacubeMore::ProbacubeMore_PCALevProb::-996660000::-10000"
)
sgems.execute(
    "ClearPropertyValueIf ProbacubeMore::ProbacubeMore_PCAPoorLevProb::-996660000::-10000"
)
Beispiel #5
0
import sgems

# 1. The output path of where to store the generated facies models
path = '/home/lewisli/Desktop/WCA/Generated Models/'

# Load project files
sgems.execute("LoadProject \
/home/lewisli/Desktop/WCA/InitialData/InitialData.prj/")

sgems.execute("ClearPropertyValueIf ProbacubeLess::ProbacubeLess_PCAChanProb::-996660000::-10000")
sgems.execute("ClearPropertyValueIf ProbacubeLess::ProbacubeLess_PCALevProb::-996660000::-10000")
sgems.execute("ClearPropertyValueIf ProbacubeLess::ProbacubeLess_PCAPoorLevProb::-996660000::-10000")
sgems.execute("ClearPropertyValueIf ProbacubeLess::ProbacubeLess_PCAShaleProb::-996660000::-10000")

sgems.execute("ClearPropertyValueIf ProbacubeMore::ProbacubeMore_PCAChanProb::-996660000::-10000")
sgems.execute("ClearPropertyValueIf ProbacubeMore::ProbacubeMore_PCALevProb::-996660000::-10000")
sgems.execute("ClearPropertyValueIf ProbacubeMore::ProbacubeMore_PCAPoorLevProb::-996660000::-10000")
sgems.execute("ClearPropertyValueIf ProbacubeMore::ProbacubeMore_PCAShaleProb::-996660000::-10000")

sgems.execute("ClearPropertyValueIf Probacubes::Probacubes_PCAChanProb::-996660000::-10000")
sgems.execute("ClearPropertyValueIf Probacubes::Probacubes_PCALevProb::-996660000::-10000")
sgems.execute("ClearPropertyValueIf Probacubes::Probacubes_PCAPoorLevProb::-996660000::-10000")
sgems.execute("ClearPropertyValueIf Probacubes::Probacubes_PCAShaleProb::-996660000::-10000")

# Part 1: Run Snesim to generate facies realizations
#TI_for_Study = ["TI1","TI2","TI3"]  # TI used for the current study
TI_for_Study = ["TI1","TI2","TI3"]  # TI used for the current study
seed = 125622
facies_name=["NTG"]
Pcube_name = ["Probacubes","ProbacubeLess",
"ProbacubeMore"]  # name of the probability cube to be employed (1st one should be used by default)
Beispiel #6
0
import sgems
sgems.execute("ReloadPythonPlugins")
Beispiel #7
0
import sgems
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/velocity_grid::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/well_data::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray1::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray2::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray3::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray4::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray5::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray6::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray7::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray8::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray9::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray10::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray11::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray12::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray13::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray14::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray15::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray16::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray17::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray18::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray19::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray20::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray21::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray22::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray23::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray24::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level6_2ray25::All')
sgems.execute('SaveGeostatGrid  velocity_grid (50x40x1)::C:/Users/noe/Documents/MATLAB/Research/DataSets/SixIterations.prj/velocity_grid (50x40x1)::s-gems')
sgems.execute('SaveGeostatGrid  well_data::C:/Users/noe/Documents/MATLAB/Research/DataSets/SixIterations.prj/well_data::s-gems')
sgems.execute('SaveGeostatGrid  Level6_2ray1::C:/Users/noe/Documents/MATLAB/Research/DataSets/SixIterations.prj/Level6_2ray1::s-gems')
Beispiel #8
0
def ar2gasprop_to_ar2gems(grid, grid_name, prop, prop_name):
    sgems.execute('NewCartesianGrid  {}::{}::{}::{}::{}::{}::{}::{}::{}::{}::0,00'.format(grid_name,
                                                                                          grid.dim()[0], grid.dim()[1], grid.dim()[2],
                                                                                          grid.cell_size()[0], grid.cell_size()[1], grid.cell_size()[2],
                                                                                          grid.origin()[0], grid.origin()[1], grid.origin()[2]))
    sgems.set_property(grid_name, prop_name, prop)
Beispiel #9
0
import sgems
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/LargerVariogram.prj/velocity_grid (50x40x1)::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/LayerCake/well_data::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray1::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray2::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray3::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray4::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray5::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray6::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray7::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray8::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray9::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray10::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray11::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray12::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray13::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray14::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray15::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray16::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray17::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray18::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray19::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray20::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray21::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray22::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray23::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray24::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Google Drive/SGeMS/DataSets/RayRepo/ThirdLevel6ray25::All')
sgems.execute('SaveGeostatGrid  velocity_grid (50x40x1)::C:/Users/noe/Google Drive/SGeMS/DataSets/LargerVariogram.prj/velocity_grid (50x40x1)::s-gems')
sgems.execute('SaveGeostatGrid  well_data::C:/Users/noe/Google Drive/SGeMS/DataSets/LargerVariogram.prj/well_data::s-gems')
sgems.execute('SaveGeostatGrid  ThirdLevel6ray1::C:/Users/noe/Google Drive/SGeMS/DataSets/LargerVariogram.prj/ThirdLevel6ray1::s-gems')
Beispiel #10
0
from os import listdir
from os.path import isfile, join
import shutil, errno
import sgems


baseProjectName = '/home/lewisli/Desktop/WCA/generatedModels/FaciesModel/FaciesModel.prj'
sgems.execute('LoadProject ' + baseProjectName)

numFaciesModelperGrid = 10
grids = ['TI1']
numCopiesPerFacieModel = 10

for k in range(0,len(grids)):
    # First make copies of facies models
    for i in range(0,numFaciesModelperGrid):
        propertyName = grids[k] + '::NTG__real' + str(i)
        tempName = grids[k] + '::temp__real' + str(i)
        sgems.execute('CopyProperty ' + propertyName + '::' + tempName + '::0::0')
        
        # Delete original
        sgems.execute('DeleteObjectProperties '+ propertyName)
    
    for i in range(0,numFaciesModelperGrid):  
        tempName = grids[k] + '::temp__real' + str(i)  
        for j in range(0,numCopiesPerFacieModel):
            newID = k*numFaciesModelperGrid*numCopiesPerFacieModel + i*numCopiesPerFacieModel + j
            newPropertyName = grids[k] + '::NTG__real' + str(newID)
            
            print newPropertyName + ' from ' + tempName
Beispiel #11
0
numRealMaxMin = 0
varioMax = []
# Generate Variogram range
for i in range(0, numReal):
    varioMax.append(triRand(30, 40, 35))
    print varioMax[i]

# Number of histograms
numHisto = 10

TI_for_Study = [
    "TI1", "TI2", "TI3", "TI4", "TI5", "TI8", "TI9", "TI10", "TI11", "TI12",
    "TI13"
]  # TI used for the current study

sgems.execute("LoadObjectFromFile " + initalDataDir +
              "WellsDataPerFacies/WellData_facies0_Sgems.out::All")
sgems.execute("LoadObjectFromFile " + initalDataDir +
              "WellsDataPerFacies/WellData_facies1_Sgems.out::All")
sgems.execute("LoadObjectFromFile " + initalDataDir +
              "WellsDataPerFacies/WellData_facies2_Sgems.out::All")
sgems.execute("LoadObjectFromFile " + initalDataDir +
              "WellsDataPerFacies/WellData_facies3_Sgems.out::All")

# We have to clear well data of NaN
sgems.execute("ClearPropertyValueIf WellFacies0::Porosity::-996660000::-10000")
sgems.execute("ClearPropertyValueIf WellFacies1::Porosity::-996660000::-10000")
sgems.execute("ClearPropertyValueIf WellFacies2::Porosity::-996660000::-10000")
sgems.execute("ClearPropertyValueIf WellFacies3::Porosity::-996660000::-10000")

sgems.execute('NewCartesianGrid  poro::' + nx + '::' + ny + '::' + nz +
              '::1.0::1.0::1.0::0::0::0')
Beispiel #12
0
import sgems
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/ComOn.prj/velocity_grid (40x40x1)::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/well_data::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/ComOn.prj/All_StraightRays::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray1::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray2::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray3::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray4::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray5::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray6::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray7::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray8::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray9::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray10::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray11::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray12::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray13::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray14::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray15::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray16::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray17::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray18::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray19::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray20::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray21::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray22::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray23::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray24::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray25::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray26::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/Level2_1ray27::All')
Beispiel #13
0
import sgems
sgems.execute('LoadObjectFromFile  D:/Users/jwhite/Projects/Broward/Geostats/SGEMS/Layer1_thk.sgems::s-gems')
data = sgems.get_property('Layer_Thk','Layer_1_thk NGVD_meters')
print data
import sgems
from sgemsbook import path

sgems.execute('LoadProject ' + path + '/DataSet/stanford5_filtersim.prj')

sgems.execute('RunGeostatAlgorithm  filtersim_std::/GeostatParamUtils/XML::<parameters>  <algorithm name="filtersim_std" />     <GridSelector_Sim  value="grid_2d_region"  />     <Property_Name_Sim  value="filtersim_aff_r0" />     <Nb_Realizations  value="1" />     <Seed  value="211175" />     <PropertySelector_Training  grid="TI_2D_large"   property="thick"  />     <Scan_Template  value="11 11 5" />     <Patch_Template_ADVANCED  value="7 7 3" />     <Is_Contv  value="0"  />     <Is_Catv  value="1"  />     <Treat_Cate_As_Cont  value="0"  />     <Nb_Facies  value="3" />     <Marginal_Cpdf  value="0.63 0.33 0.04" />     <Hard_Data  grid="well"   property="facies"  />     <Use_SoftField  value="0"  />     <SoftData_properties count="0"   value=""  />     <TauModelObject  value="1 1" />     <Use_Region  value="1"  />     <Region_Indicator_Prop  value="affinity"  />     <Use_Previous_Simulation  value="0"  />     <Previous_Simulation_Prop  value="affinity"  />     <Active_Region_Code  value="0" />     <Cmin_Replicates  value="10 10 10" />     <Nb_Multigrids_ADVANCED  value="3" />     <Constraint_Marginal_ADVANCED  value="0.5" />     <Filter_Default  value="1"  />     <Filter_User_Define  value="0"  />     <User_Def_Filter_File  value="" />     <CrossPartition  value="1"  />     <KMeanPartition  value="0"  />     <Nb_Bins_ADVANCED  value="4" />     <Nb_Bins_ADVANCED2  value="2" />     <Use_Normal_Dist  value="1"  />     <Debug_Level  value="0" />     <Data_Weights  value="0.5 0.3 0.2" />   </parameters>   ')

sgems.execute('RunGeostatAlgorithm  filtersim_std::/GeostatParamUtils/XML::<parameters>  <algorithm name="filtersim_std" />     <GridSelector_Sim  value="grid_2d_region"  />     <Property_Name_Sim  value="filtersim_aff_r1" />     <Nb_Realizations  value="1" />     <Seed  value="211175" />     <PropertySelector_Training  grid="TI_2D_large"   property="facies"  />     <Scan_Template  value="11 11 5" />     <Patch_Template_ADVANCED  value="7 7 3" />     <Is_Contv  value="0"  />     <Is_Catv  value="1"  />     <Treat_Cate_As_Cont  value="0"  />     <Nb_Facies  value="3" />     <Marginal_Cpdf  value="0.63 0.33 0.04" />     <Hard_Data  grid="well"   property="facies"  />     <Use_SoftField  value="0"  />     <SoftData_properties count="0"   value=""  />     <TauModelObject  value="1 1" />     <Use_Region  value="1"  />     <Region_Indicator_Prop  value="affinity"  />     <Use_Previous_Simulation  value="1"  />     <Previous_Simulation_Prop  value="filtersim_aff_r0__real0"  />     <Active_Region_Code  value="1" />     <Cmin_Replicates  value="10 10 10" />     <Nb_Multigrids_ADVANCED  value="3" />     <Constraint_Marginal_ADVANCED  value="0.5" />     <Filter_Default  value="1"  />     <Filter_User_Define  value="0"  />     <User_Def_Filter_File  value="" />     <CrossPartition  value="1"  />     <KMeanPartition  value="0"  />     <Nb_Bins_ADVANCED  value="4" />     <Nb_Bins_ADVANCED2  value="2" />     <Use_Normal_Dist  value="1"  />     <Debug_Level  value="0" />     <Data_Weights  value="0.5 0.3 0.2" />   </parameters>   ')

sgems.execute('RunGeostatAlgorithm  filtersim_std::/GeostatParamUtils/XML::<parameters>  <algorithm name="filtersim_std" />     <GridSelector_Sim  value="grid_2d_region"  />     <Property_Name_Sim  value="filtersim_aff_r2" />     <Nb_Realizations  value="1" />     <Seed  value="211175" />     <PropertySelector_Training  grid="TI_2D_large"   property="thin"  />     <Scan_Template  value="11 11 5" />     <Patch_Template_ADVANCED  value="7 7 3" />     <Is_Contv  value="0"  />     <Is_Catv  value="1"  />     <Treat_Cate_As_Cont  value="0"  />     <Nb_Facies  value="3" />     <Marginal_Cpdf  value="0.63 0.33 0.04" />     <Hard_Data  grid="well"   property="facies"  />     <Use_SoftField  value="0"  />     <SoftData_properties count="0"   value=""  />     <TauModelObject  value="1 1" />     <Use_Region  value="1"  />     <Region_Indicator_Prop  value="affinity"  />     <Use_Previous_Simulation  value="1"  />     <Previous_Simulation_Prop  value="filtersim_aff_r1__real0"  />     <Active_Region_Code  value="2" />     <Cmin_Replicates  value="10 10 10" />     <Nb_Multigrids_ADVANCED  value="3" />     <Constraint_Marginal_ADVANCED  value="0.5" />     <Filter_Default  value="1"  />     <Filter_User_Define  value="0"  />     <User_Def_Filter_File  value="" />     <CrossPartition  value="1"  />     <KMeanPartition  value="0"  />     <Nb_Bins_ADVANCED  value="4" />     <Nb_Bins_ADVANCED2  value="2" />     <Use_Normal_Dist  value="1"  />     <Debug_Level  value="0" />     <Data_Weights  value="0.5 0.3 0.2" />   </parameters>   ')

Beispiel #15
0
import sgems
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/FirstSlow.prj/velocity_grid (40x40x1)::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/well_data::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays1::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays2::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays3::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays4::All')
sgems.execute('SaveGeostatGrid  velocity_grid (40x40x1)::C:/Users/noe/Documents/Research/DataSets/FirstSlow.prj/velocity_grid (40x40x1)::s-gems')
sgems.execute('SaveGeostatGrid  well_data::C:/Users/noe/Documents/Research/DataSets/FirstSlow.prj/well_data::s-gems')
sgems.execute('SaveGeostatGrid  FirstLevelRays1::C:/Users/noe/Documents/Research/DataSets/FirstSlow.prj/FirstLevelRays1::s-gems')
sgems.execute('SaveGeostatGrid  FirstLevelRays2::C:/Users/noe/Documents/Research/DataSets/FirstSlow.prj/FirstLevelRays2::s-gems')
sgems.execute('SaveGeostatGrid  FirstLevelRays3::C:/Users/noe/Documents/Research/DataSets/FirstSlow.prj/FirstLevelRays3::s-gems')
sgems.execute('SaveGeostatGrid  FirstLevelRays4::C:/Users/noe/Documents/Research/DataSets/FirstSlow.prj/FirstLevelRays4::s-gems')
sgems.execute('RunGeostatAlgorithm  bssim::/GeostatParamUtils/XML::<parameters> <algorithm name="bssim" />     <Grid_Name  value="velocity_grid (40x40x1)"  /> <Property_Name  value="Level1_1" />     <Nb_Realizations  value="2" />     <Seed  value="40807263" />    <Kriging_Type  value="Simple Kriging (SK)"  /> <SK_mean  value="5" /> <Block_Cov_Approach  value="FFT with Covariance-Table"  /> <Simulation_Path  value="Block First"  /> <Check_Block_Reproduction  value="1"  /> <Generate_Etype  value="1"  /> <Debug_Info  value="None (level 0)"  /> <Hard_Data  grid="well_data"   property="velocities"  />   <Assign_Hard_Data  value="1"  /> <Max_Conditioning_Data_Point  value="12" />  <Search_Ellipsoid_Point  value="0.5 0.5 0.25  90 0 0" /> <Max_Conditioning_Data_Block  value="12" />  <Search_Points_In_Blocks_Method  value="Include All"  /> <Max_Cond_Points_In_Blocks  value="12" /> <Block_From_File  value="0"  /> <Block_From_Pset  value="1"  />  <Block_Data_File  value="D:/DataSet/tomography_blocks.dat" /> <Number_of_Blocks  value="4" /> <Block_Average_Values  value="5.0334 5.1435 5.3713 4.4778 " />     <Consider_Block_Error  value="0"  />  <Block_Grid_0 value="FirstLevelRays1"  /> <Block_Grid_1 value="FirstLevelRays2"  /> <Block_Grid_2 value="FirstLevelRays3"  /> <Block_Grid_3 value="FirstLevelRays4"  /> <Variogram_Cov  nugget="0.1" structures_count="1"  > <structure_1  contribution="0.9"  type="Spherical"   > <ranges max="0.5"  medium="0.5"  min="0.25"   /> <angles x="90"  y="0"  z="0"   />    </structure_1>  </Variogram_Cov>  <cdf_type  value="Soares"  />     <Gaussian_mean  value="1" /> <Gaussian_variance  value="1" />     <LN_mean  value="1" />  <LN_variance  value="1" />     <U_min  value="0" />     <U_max  value="1" /> <nonParamCdf  ref_on_file ="0"  ref_on_grid ="1" break_ties ="0" filename =""   grid ="All_StraightRays"  property ="velocities"> <LTI_type  function ="Power"  extreme ="3.7"  omega ="3" /> <UTI_type  function ="Power"  extreme ="7.1"  omega ="0.333" />  </nonParamCdf>    <is_local_correction  value="1"  />   </parameters>') 
sgems.execute('SaveGeostatGrid velocity_grid (40x40x1)::C:/Users/noe/Documents/Research/DataSets/FirstSlow.prj/velocity_grid (40x40x1)::s-gems')
sgems.execute('SaveGeostatGrid velocity_grid (40x40x1)::C:/Users/noe/Documents/Research/DataSets/FirstSlow.prj/velocity_grid (40x40x1)Level1__1::s-gems::0::Level1_1__real0')
sgems.execute('SaveGeostatGrid velocity_grid (40x40x1)::C:/Users/noe/Documents/Research/DataSets/FirstSlow.prj/velocity_grid (40x40x1)Level1__2::s-gems::0::Level1_1__real1')
Beispiel #16
0
import sgems
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/DoubleWidthMoreVariogramY.prj/velocity_grid (50x80x1)::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/well_data::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray1::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray2::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray3::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray4::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray5::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray6::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray7::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray8::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray9::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray10::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray11::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray12::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray13::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray14::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray15::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray16::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray17::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray18::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray19::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray20::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray21::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray22::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray23::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray24::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray25::All')
sgems.execute('SaveGeostatGrid  velocity_grid (50x80x1)::C:/Users/noe/Documents/MATLAB/Research/DataSets/DoubleWidthMoreVariogramY.prj/velocity_grid (50x80x1)::s-gems')
sgems.execute('SaveGeostatGrid  well_data::C:/Users/noe/Documents/MATLAB/Research/DataSets/DoubleWidthMoreVariogramY.prj/well_data::s-gems')
sgems.execute('SaveGeostatGrid  Level2_1ray1::C:/Users/noe/Documents/MATLAB/Research/DataSets/DoubleWidthMoreVariogramY.prj/Level2_1ray1::s-gems')
Beispiel #17
0
import sgems
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/FourItrWOGauss.prj/velocity_grid (50x80x1)::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/well_data::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray1::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray2::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray3::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray4::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray5::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray6::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray7::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray8::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray9::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray10::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray11::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray12::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray13::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray14::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray15::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray16::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray17::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray18::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray19::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray20::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray21::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray22::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray23::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray24::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray25::All')
sgems.execute('SaveGeostatGrid  velocity_grid (50x80x1)::C:/Users/noe/Documents/MATLAB/Research/DataSets/FourItrWOGauss.prj/velocity_grid (50x80x1)::s-gems')
sgems.execute('SaveGeostatGrid  well_data::C:/Users/noe/Documents/MATLAB/Research/DataSets/FourItrWOGauss.prj/well_data::s-gems')
sgems.execute('SaveGeostatGrid  Level4_4ray1::C:/Users/noe/Documents/MATLAB/Research/DataSets/FourItrWOGauss.prj/Level4_4ray1::s-gems')
Beispiel #18
0
import sgems
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/TOC_All_Blocks::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/TOC_Cond_Data::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block1::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block2::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block3::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block4::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block5::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block6::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block7::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block8::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block9::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block10::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block11::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block12::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block13::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block14::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block15::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block16::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block17::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block18::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block19::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block20::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block21::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block22::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block23::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block24::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block25::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block26::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block27::All')
sgems.execute('LoadObjectFromFile C:/Users/ahino/OneDrive/Documents/School/Research/Code/Toc Problem/TOC_Project.prj/block28::All')
# Path where initial data is stored
initalDataDir = '/home/lewisli/Desktop/WCA/initialData/'


# Use sgsim for each facies for vshale
# Since we are not conditioning

def triRand(a, c, b):
    c = float(c)
    t = (c-a)/(b-a)
    y = sqrt(random())
    d = a if random() < t else b
    return d + (c-d) * y


sgems.execute('LoadProject  /home/lewisli/Desktop/WCA/generatedModels/Porosity/PoroProjectFull.prj')

print 'here'

nbfacies = 4
SaveRealName='VShalesity.prj'
nx = '78'
ny = '59'
nz = '116'
seed = 14071789

MinValueVShale = [0.20736,0.040977,0.032918,4.0396e-06]
MaxValueVShale = [1.0001,0.97252,0.947399,0.882]

MinValuePoro = [1.4e-07,0.002,0.001,0.04]
MaxValuePoro = [0.2199,0.304,0.334,0.3604]
import sgems

for i in range(1,19):
	sgems.execute('RunGeostatAlgorithm  kriging::/GeostatParamUtils/XML::<parameters>  <algorithm name="kriging" />     <Grid_Name  value="kriging grid"  />     <Property_Name  value="titi" />     <Kriging_Type  type="Simple Kriging (SK)" >    <parameters mean="0.27" />  </Kriging_Type>    <Hard_Data  grid="sample data"   property="porosity"  />     <Search_Ellipsoid  value="80 80 80  0 0 0" />    <Min_Conditioning_Data  value="0" />     <Max_Conditioning_Data  value="' + str(10+5*i) +'" />     <Variogram  nugget="0" structures_count="1"  >    <structure_1  contribution="0.003"  type="Spherical"   >      <ranges max="38"  medium="38"  min="38"   />      <angles x="0"  y="0"  z="0"   />    </structure_1>  </Variogram>  </parameters>')  
Beispiel #21
0
import sgems

sgems.execute(
    'LoadObjectFromFile  D:/Users/jwhite/Projects/Broward/Geostats/SGEMS/Layer1_thk.sgems::s-gems'
)
data = sgems.get_property('Layer_Thk', 'Layer_1_thk NGVD_meters')
print data
Beispiel #22
0
import sgems
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/MoreSources.prj/velocity_grid (40x40x1)::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/well_data::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/MoreSources.prj/All_StraightRays::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays1::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays2::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays3::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays4::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays5::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays6::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays7::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays8::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays9::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays10::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays11::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays12::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays13::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays14::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays15::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays16::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays17::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays18::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays19::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays20::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays21::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays22::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays23::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays24::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays25::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays26::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays27::All')
Beispiel #23
0
import sgems
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/velocity_grid::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/well_data::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray1::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray2::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray3::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray4::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray5::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray6::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray7::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray8::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray9::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray10::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray11::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray12::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray13::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray14::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray15::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray16::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray17::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray18::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray19::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray20::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray21::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray22::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray23::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray24::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray25::All')
sgems.execute('SaveGeostatGrid  velocity_grid (50x45x1)::C:/Users/noe/Documents/MATLAB/Research/DataSets/TwoIterations.prj/velocity_grid (50x45x1)::s-gems')
sgems.execute('SaveGeostatGrid  well_data::C:/Users/noe/Documents/MATLAB/Research/DataSets/TwoIterations.prj/well_data::s-gems')
sgems.execute('SaveGeostatGrid  Level2_1ray1::C:/Users/noe/Documents/MATLAB/Research/DataSets/TwoIterations.prj/Level2_1ray1::s-gems')
Beispiel #24
0
import sgems
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/LargerSearchEllipsoid.prj/velocity_grid (50x40x1)::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/well_data::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray1::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray2::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray3::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray4::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray5::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray6::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray7::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray8::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray9::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray10::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray11::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray12::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray13::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray14::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray15::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray16::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray17::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray18::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray19::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray20::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray21::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray22::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray23::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray24::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/RayRepo/ThirdLevel6ray25::All')
sgems.execute('SaveGeostatGrid  velocity_grid (50x40x1)::C:/Users/noe/Documents/MATLAB/Research/LargerSearchEllipsoid.prj/velocity_grid (50x40x1)::s-gems')
sgems.execute('SaveGeostatGrid  well_data::C:/Users/noe/Documents/MATLAB/Research/LargerSearchEllipsoid.prj/well_data::s-gems')
sgems.execute('SaveGeostatGrid  ThirdLevel6ray1::C:/Users/noe/Documents/MATLAB/Research/LargerSearchEllipsoid.prj/ThirdLevel6ray1::s-gems')
import sgems
import math

seed = 14071789
num_realizations = 5000
nrow,ncol = 158,93
delr,delc = 1000,1000 

offset = (503800.0,2778800.0)
realization_path = 'D:/Users/jwhite/Projects/Broward/Geostats/SGEMS/l1_ds_reals/'
realization_prefix = 'layer1_thk_omni_ds'

prop_names = ['H_thk','Q1_thk','Q2_thk','Q3_thk','Q4_thk','Q5_thk','T1_thk','T2_thk','T3_thk',]

try:
    sgems.execute('DeleteObjects tbl29::halfmile')
except:
    pass
load_command = 'LoadObjectFromFile  D:/Users/jwhite/Projects/Broward/Geostats/12_layer/Table_29.sgems::s-gems' 
sgems.execute(load_command)

newGrid_command = 'NewCartesianGrid  halfmile::116::197::1::2650::2650::1.0::668350::288415::0'
sgems.execute(newGrid_command)
#--for OK 
#dssim_command_aniso = 'RunGeostatAlgorithm  dssim::/GeostatParamUtils/XML::<parameters>  '+\
#                '<algorithm name="dssim" />     <Grid_Name  value="1000_grid"  />     '+\
#                '<Property_Name  value="DsSim" />     <Nb_Realizations  value="'+str(num_realizations)+'" />     '+\
#                '<Seed  value="14071789" />     <Kriging_Type  value="Ordinary Kriging (OK)"  />     '+\
#                '<Hard_Data  grid="Layer_Thk"   property="Layer_1_thk"  />     <Assign_Hard_Data  value="1"  />     '+\
#                '<Max_Conditioning_Data  value="12" />     <Search_Ellipsoid  value="500000 55000 0  10 0 0" />    '+\
#                '<cdf_type  value="Uniform"  />     <LN_mean  value="1" />     <LN_variance  value="1" />     '+\
Beispiel #26
0
MaxValuePoro = [0.23558,0.30323,0.33307,0.36032]
numReal = 1
numRealMaxMin = 0
varioMax = []
# Generate Variogram range
for i in range (0,numReal):
    varioMax.append(triRand(30,40,35))
    print varioMax[i]
    

# Number of histograms
numHisto = 10;

TI_for_Study = ["TI1","TI2","TI3","TI4","TI5","TI8","TI9","TI10","TI11","TI12","TI13"]  # TI used for the current study

sgems.execute("LoadObjectFromFile " + initalDataDir + "WellsDataPerFacies/WellData_facies0_Sgems.out::All")
sgems.execute("LoadObjectFromFile " + initalDataDir + "WellsDataPerFacies/WellData_facies1_Sgems.out::All")
sgems.execute("LoadObjectFromFile " + initalDataDir + "WellsDataPerFacies/WellData_facies2_Sgems.out::All")
sgems.execute("LoadObjectFromFile " + initalDataDir + "WellsDataPerFacies/WellData_facies3_Sgems.out::All")

# We have to clear well data of NaN
sgems.execute("ClearPropertyValueIf WellFacies0::Porosity::-996660000::-10000")
sgems.execute("ClearPropertyValueIf WellFacies1::Porosity::-996660000::-10000")
sgems.execute("ClearPropertyValueIf WellFacies2::Porosity::-996660000::-10000")
sgems.execute("ClearPropertyValueIf WellFacies3::Porosity::-996660000::-10000")

sgems.execute('NewCartesianGrid  poro::'+nx+'::'+ny+'::'+nz+'::1.0::1.0::1.0::0::0::0')
sgems.execute('NewCartesianGrid  temp::'+nx+'::'+ny+'::'+nz+'::1.0::1.0::1.0::0::0::0')

# Run algorithm for all facies
Beispiel #27
0
import sgems
import math

seed = 14071789
num_realizations = 5000
nrow,ncol = 158,93
delr,delc = 1000,1000 

offset = (503800.0,2778800.0)
realization_path = 'D:/Users/jwhite/Projects/Broward/Geostats/SGEMS/l1_ds_reals/'
realization_prefix = 'layer1_thk_omni_ds'

try:
    sgems.execute('DeleteObjects Layer_Thk::1000_grid')
except:
    pass
load_command = 'LoadObjectFromFile  D:/Users/jwhite/Projects/Broward/Geostats/SGEMS/Layer1_thk.sgems::s-gems' 
sgems.execute(load_command)

newGrid_command = 'NewCartesianGrid  1000_grid::'+str(ncol)+'::'+str(nrow)+'::1::'+str(delc)+'::'+str(delr)+'::1.0::'+str(offset[0])+'::'+str(offset[1])+'::0'
sgems.execute(newGrid_command)
#--for OK 
dssim_command_aniso = 'RunGeostatAlgorithm  dssim::/GeostatParamUtils/XML::<parameters>  '+\
                '<algorithm name="dssim" />     <Grid_Name  value="1000_grid"  />     '+\
                '<Property_Name  value="DsSim" />     <Nb_Realizations  value="'+str(num_realizations)+'" />     '+\
                '<Seed  value="14071789" />     <Kriging_Type  value="Ordinary Kriging (OK)"  />     '+\
                '<Hard_Data  grid="Layer_Thk"   property="Layer_1_thk"  />     <Assign_Hard_Data  value="1"  />     '+\
                '<Max_Conditioning_Data  value="12" />     <Search_Ellipsoid  value="500000 55000 0  10 0 0" />    '+\
                '<cdf_type  value="Uniform"  />     <LN_mean  value="1" />     <LN_variance  value="1" />     '+\
                '<U_min  value="0" />     <U_max  value="30" />     <nonParamCdf  ref_on_file ="0"  ref_on_grid ="1"  '+\
                'break_ties ="0" filename =""   grid =""  property ="">  <LTI_type  function ="Power"  extreme ="0"  '+\
Beispiel #28
0
import sgems
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/DoubleWidth.prj/velocity_grid (50x80x1)::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/well_data::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray1::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray2::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray3::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray4::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray5::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray6::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray7::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray8::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray9::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray10::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray11::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray12::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray13::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray14::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray15::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray16::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray17::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray18::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray19::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray20::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray21::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray22::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray23::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray24::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level2_1ray25::All')
sgems.execute('SaveGeostatGrid  velocity_grid (50x80x1)::C:/Users/noe/Documents/MATLAB/Research/DataSets/DoubleWidth.prj/velocity_grid (50x80x1)::s-gems')
sgems.execute('SaveGeostatGrid  well_data::C:/Users/noe/Documents/MATLAB/Research/DataSets/DoubleWidth.prj/well_data::s-gems')
sgems.execute('SaveGeostatGrid  Level2_1ray1::C:/Users/noe/Documents/MATLAB/Research/DataSets/DoubleWidth.prj/Level2_1ray1::s-gems')
Beispiel #29
0
import sgems
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/velocity_grid::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/well_data::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level5_2ray1::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level5_2ray2::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level5_2ray3::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level5_2ray4::All')
sgems.execute('SaveGeostatGrid  velocity_grid (50x40x1)::C:/Users/noe/Documents/MATLAB/Research/DataSets/FiveIterations.prj/velocity_grid (50x40x1)::s-gems')
sgems.execute('SaveGeostatGrid  well_data::C:/Users/noe/Documents/MATLAB/Research/DataSets/FiveIterations.prj/well_data::s-gems')
sgems.execute('SaveGeostatGrid  Level5_2ray1::C:/Users/noe/Documents/MATLAB/Research/DataSets/FiveIterations.prj/Level5_2ray1::s-gems')
sgems.execute('SaveGeostatGrid  Level5_2ray2::C:/Users/noe/Documents/MATLAB/Research/DataSets/FiveIterations.prj/Level5_2ray2::s-gems')
sgems.execute('SaveGeostatGrid  Level5_2ray3::C:/Users/noe/Documents/MATLAB/Research/DataSets/FiveIterations.prj/Level5_2ray3::s-gems')
sgems.execute('SaveGeostatGrid  Level5_2ray4::C:/Users/noe/Documents/MATLAB/Research/DataSets/FiveIterations.prj/Level5_2ray4::s-gems')
sgems.execute('RunGeostatAlgorithm  bssim::/GeostatParamUtils/XML::<parameters> <algorithm name="bssim" />     <Grid_Name  value="velocity_grid (50x40x1)"  /> <Property_Name  value="Level5_2" />     <Nb_Realizations  value="1" />     <Seed  value="44615721" />    <Kriging_Type  value="Simple Kriging (SK)"  /> <SK_mean  value="2276.7" /> <Block_Cov_Approach  value="FFT with Covariance-Table"  /> <Simulation_Path  value="Block First"  /> <Check_Block_Reproduction  value="1"  /> <Generate_Etype  value="1"  /> <Debug_Info  value="None (level 0)"  /> <Hard_Data  grid="well_data"   property="velocities"  />   <Assign_Hard_Data  value="1"  /> <Max_Conditioning_Data_Point  value="12" />  <Search_Ellipsoid_Point  value="3 1 1  90 0 0" /> <Max_Conditioning_Data_Block  value="12" />  <Search_Points_In_Blocks_Method  value="Not Include"  /> <Max_Cond_Points_In_Blocks  value="12" /> <Block_From_File  value="0"  /> <Block_From_Pset  value="1"  />  <Block_Data_File  value="D:/DataSet/tomography_blocks.dat" /> <Number_of_Blocks  value="4" /> <Block_Average_Values  value="1788.8282 2657.8604 2511.8348 4105.2632 " />     <Consider_Block_Error  value="0"  />  <Block_Grid_0 value="Level5_2ray1"  /> <Block_Grid_1 value="Level5_2ray2"  /> <Block_Grid_2 value="Level5_2ray3"  /> <Block_Grid_3 value="Level5_2ray4"  /> <Variogram_Cov  nugget="0.1" structures_count="1"  > <structure_1  contribution="0.9"  type="Spherical"   > <ranges max="10"  medium="2"  min="1"   /> <angles x="90"  y="0"  z="0"   />    </structure_1>  </Variogram_Cov>  <cdf_type  value="Soares"  />     <Gaussian_mean  value="1" /> <Gaussian_variance  value="1" />     <LN_mean  value="1" />  <LN_variance  value="1" />     <U_min  value="0" />     <U_max  value="1" /> <nonParamCdf  ref_on_file ="0"  ref_on_grid ="1" break_ties ="0" filename =""   grid ="well_data"  property ="velocities"> <LTI_type  function ="Power"  extreme ="1000"  omega ="3" /> <UTI_type  function ="Power"  extreme ="5000"  omega ="0.333" />  </nonParamCdf>    <is_local_correction  value="1"  />   </parameters>') 
sgems.execute('SaveGeostatGrid velocity_grid (50x40x1)::C:/Users/noe/Documents/MATLAB/Research/DataSets/FiveIterations.prj/velocity_grid (50x40x1)::s-gems')
sgems.execute('SaveGeostatGrid velocity_grid (50x40x1)::C:/Users/noe/Documents/MATLAB/Research/DataSets/FiveIterations.prj/velocity_grid (50x40x1)Level5__2::s-gems::0::Level5_2__real0')
Beispiel #30
0
def ar2gasgrid_to_ar2gems(grid_name, grid):
    sgems.execute('NewCartesianGrid  {}::{}::{}::{}::{}::{}::{}::{}::{}::{}::0,00'.format(grid_name,
                                                                                          grid.dim()[0], grid.dim()[1], grid.dim()[2],
                                                                                          grid.cell_size()[0], grid.cell_size()[1], grid.cell_size()[2],
                                                                                          grid.origin()[0], grid.origin()[1], grid.origin()[2]))
Beispiel #31
0
		#print inner_list
        # in alternative, if you need to use the file content as numbers
        # inner_list = [int(elt.strip()) for elt in line.split(',')]
        list_of_lists.append(inner_list)
print list_of_lists[0][1]


#loop through each row of table
Nbreal=1  # Number of model generated by each combination of the parameter

nx = '200'
ny = '200'
nz = '1'

# import the parameters from the design table
sgems.execute('NewCartesianGrid  temp::'+nx+'::'+ny+'::'+nz+'::1.0::1.0::1.0::0::0::0')
for i in range(len(list_of_lists)):
	sgems.execute('RunGeostatAlgorithm  TetrisTiGen::/GeostatParamUtils/XML::<parameters>  <algorithm name="TetrisTiGen" />     <Grid value="temp" region=""  />  \
	   <Seed  value="'+str(list_of_lists[i][11])+'" />   \
	   <Nb_Realizations  value="'+str(Nbreal)+'" />   \
	   <Property  value="barrier_flow'+str(i)+'" />   \
	   <XmlFile>   <Tetris>   <Objects numObjects="1">    <ObjectManager numElements="1" numOperations="0" type="Element" name="fault">     <Element value="1" name="Cuboid">      <ShapeManager type="Cuboid" name="Cuboid" code="1">       <Dimension>        <dx>         <CDF name="Constant">          <mean>      \
	   <Params value="'+str(list_of_lists[i][8])+'" name="Constant"/>      \
	   </mean>         </CDF>        </dx>        <dy>         <CDF name="Constant">          <mean>           <Params value="1" name="Constant"/>          </mean>         </CDF>        </dy>        <dz>         <CDF name="Constant">          <mean>           <Params value="1" name="Constant"/>          </mean>         </CDF>        </dz>       </Dimension>       <Rotation>        <ax>         <CDF name="Constant">          <mean>           <Params value="0" name="Constant"/>          </mean>         </CDF>        </ax>        <ay>         <CDF name="Constant">          <mean>           <Params value="0" name="Constant"/>          </mean>         </CDF>        </ay>        <az>         <CDF name="Constant">          <mean>          \
	   <Params value="'+str(list_of_lists[i][9])+'" name="Constant"/>      \
	   </mean>         </CDF>        </az>       </Rotation>      </ShapeManager>     </Element>    </ObjectManager>   </Objects>   <Simulations numSimulations="1">    <Simulation name="Simulation 1">     <ObjectSimulated type="Single Object" name="fault"/>   \
	   <Stop target="'+str(list_of_lists[i][10])+'" type="proportion"/>    \
	   <positioning type="random"/>    </Simulation>   </Simulations>  </Tetris>  </XmlFile>   </parameters>   ')
	filename='R'+str(i)+'_barrier.out::gslib'
	sgems.execute('SaveGeostatGrid  temp::'+path+filename+'::0::barrier_flow'+str(i)+'__real0')
import os

import sgems

nodata = -9966699

os.chdir("RES_DIR")
sgems.execute("DeleteObjects computation_grid")
sgems.execute("DeleteObjects PROJECT_NAME")
sgems.execute("DeleteObjects finished")

for file in OBJECT_FILES:
    sgems.execute("LoadObjectFromFile  {}::All".format(file))

sgems.execute("NewCartesianGrid  computation_grid::GRID")

#~sgems.execute('RunGeostatAlgorithm  ALGORITHM_NAME::/GeostatParamUtils/XML::ALGORITHM_XML')
#~sgems.execute('SaveGeostatGrid  computation_grid::FEATURE_OUTPUT.grid::gslib::0::OUTPUT_LIST')
#~sgems.execute('SaveGeostatGrid  computation_grid::FEATURE_OUTPUT(var).grid::gslib::1::OUTPUT_LIST_krig_var')
Beispiel #33
0
    def execute(self):
        '''# Execute the funtion read_params
        read_params(self.params)
        print self.params'''

        #Get the grid and rock type propery
        grid = self.params['propertyselectornoregion']['grid']
        prop = self.params['propertyselectornoregion']['property']

        #Get the X, Y and Z coordinates and RT property
        X = sgems.get_property(grid, '_X_')
        Y = sgems.get_property(grid, '_Y_')
        Z = sgems.get_property(grid, '_Z_')
        RT_data = sgems.get_property(grid, prop)

        # Getting properties
        grid_krig = self.params['gridselectorbasic_2']['value']
        grid_var = self.params['gridselectorbasic']['value']
        props = (self.params['orderedpropertyselector']['value']).split(';')
        n_var = int(self.params['indicator_regionalization_input']
                    ['number_of_indicator_group'])
        n_prop = int(self.params['orderedpropertyselector']['count'])
        min_cond = self.params['spinBox_2']['value']
        max_cond = self.params['spinBox']['value']

        # Error messages
        if len(grid_var) == 0 or len(grid_krig) == 0:
            print 'Select the variables'
            return False

        if n_var != n_prop:
            print 'Number of variables and number of variograms models are diferent.'
            return False

        #Creating an empty list to store the interpolated distances
        SG_OK_list = []

        # Loop in every variable
        for i in xrange(0, n_var):

            # Getting variables
            prop_HD = props[i]
            prop_name = "Interpolated_" + str(prop_HD)
            prop_name_var = "Interpolated_" + str(prop_HD) + ' krig_var'
            var_str = ''
            indicator_group = "Indicator_group_" + str(i + 1)
            elipsoide = self.params['ellipsoidinput']['value']
            n_struct = int(
                self.params['indicator_regionalization_input'][indicator_group]
                ['Covariance_input']['structures_count'])

            # Error message
            if n_struct == 0:
                print 'Variogram have no structures'
                return False

            # Loop in every variogram structure
            for j in xrange(0, n_struct):
                # Getting variogram parameters
                Structure = "Structure_" + str(j + 1)

                cov_type = self.params['indicator_regionalization_input'][
                    indicator_group]['Covariance_input'][Structure][
                        'Two_point_model']['type']

                cont = self.params['indicator_regionalization_input'][
                    indicator_group]['Covariance_input'][Structure][
                        'Two_point_model']['contribution']

                if cov_type == 'Nugget Covariance':
                    #Writing variogram parameters on a variable in nugget effect case
                    var_str = var_str + '<{} type="{}">  <Two_point_model  contribution="{}"  type="{}"   >    </Two_point_model>    </Structure_1> '.format(
                        Structure, 'Covariance', cont, cov_type, Structure)

                else:
                    range1 = self.params['indicator_regionalization_input'][
                        indicator_group]['Covariance_input'][Structure][
                            'Two_point_model']['ranges']['range1']
                    range2 = self.params['indicator_regionalization_input'][
                        indicator_group]['Covariance_input'][Structure][
                            'Two_point_model']['ranges']['range2']
                    range3 = self.params['indicator_regionalization_input'][
                        indicator_group]['Covariance_input'][Structure][
                            'Two_point_model']['ranges']['range3']

                    rake = self.params['indicator_regionalization_input'][
                        indicator_group]['Covariance_input'][Structure][
                            'Two_point_model']['angles']['rake']
                    dip = self.params['indicator_regionalization_input'][
                        indicator_group]['Covariance_input'][Structure][
                            'Two_point_model']['angles']['dip']
                    azimuth = self.params['indicator_regionalization_input'][
                        indicator_group]['Covariance_input'][Structure][
                            'Two_point_model']['angles']['azimuth']

                    # Writing variogram parameters on a variable in other cases
                    var_str = var_str + '<{} type="{}">  <Two_point_model  contribution="{}"  type="{}"   >      <ranges range1="{}"  range2="{}"  range3="{}"   />      <angles azimuth="{}"  dip="{}"  rake="{}"   />    </Two_point_model>    </{}> '.format(
                        Structure, 'Covariance', cont, cov_type, range1,
                        range2, range3, azimuth, dip, rake, Structure)

            # Calling ordinary kriging for each variable, using the variograms parameters above
            sgems.execute(
                'RunGeostatAlgorithm  kriging::/GeostatParamUtils/XML::<parameters>  <algorithm name="kriging" />     <Variogram  structures_count="{}" >    {}  </Variogram>    <ouput_kriging_variance  value="1"  />     <output_n_samples_  value="0"  />     <output_average_distance  value="0"  />     <output_sum_weights  value="0"  />     <output_sum_positive_weights  value="0"  />     <output_lagrangian  value="0"  />     <Nb_processors  value="-2"  />    <Grid_Name value="{}" region=""  />     <Property_Name  value="{}" />     <Hard_Data  grid="{}"   property="{}"   region=""  />     <Kriging_Type  type="Ordinary Kriging (OK)" >    <parameters />  </Kriging_Type>    <do_block_kriging  value="1"  />     <npoints_x  value="5" />     <npoints_y  value="5" />     <npoints_z  value="5" />     <Min_Conditioning_Data  value="{}" />     <Max_Conditioning_Data  value="{}" />     <Search_Ellipsoid  value="{}" />    <AdvancedSearch  use_advanced_search="0"></AdvancedSearch>  </parameters>'
                .format(n_struct, var_str, grid_krig, prop_name, grid_var,
                        prop_HD, min_cond, max_cond, elipsoide))

            SG_OK_list.append(sgems.get_property(grid_krig, prop_name))

            #Deleting kriged distances
            sgems.execute('DeleteObjectProperties  {}::{}'.format(
                grid_krig, prop_name))
            sgems.execute('DeleteObjectProperties  {}::{}'.format(
                grid_krig, prop_name_var))

        RT = (self.params['orderedpropertyselector']['value']).split(';')

        #Determinig geomodel based on minimum estimed signed distance function
        GeoModel = SG_OK_list[0][:]

        t = 0
        for i in range(len(SG_OK_list[0])):
            sgmin = 10e21
            for j in range(len(SG_OK_list)):
                if SG_OK_list[j][i] < sgmin:
                    sgmin = SG_OK_list[j][i]
                    t = j
            if math.isnan(SG_OK_list[j][i]):
                GeoModel[i] = float('nan')
            else:
                GeoModel[i] = (int(RT[t].split('RT_')[-1]))

        #Creating GeoModel property
        lst_props_grid = sgems.get_property_list(grid_krig)
        prop_final_data_name = 'Geologic_Model'

        if (prop_final_data_name in lst_props_grid):
            flag = 0
            i = 1
            while (flag == 0):
                test_name = prop_final_data_name + '-' + str(i)
                if (test_name not in lst_props_grid):
                    flag = 1
                    prop_final_data_name = test_name
                i = i + 1

        #Assign conditioning data to grid node
        for i in range(len(RT_data)):
            if not math.isnan(RT_data[i]):
                closest_node = sgems.get_closest_nodeid(
                    grid_krig, X[i], Y[i], Z[i])
                GeoModel[closest_node] = RT_data[i]

        sgems.set_property(grid_krig, prop_final_data_name, GeoModel)

        #Operating softmax transformation
        if self.params['softmax_check']['value'] == '1':

            gamma = float(self.params['Gamma']['value'])
            Prob_list = SG_OK_list[:]

            for i in range(len(SG_OK_list[0])):
                soma = 0
                for j in range(len(SG_OK_list)):
                    soma = soma + math.exp(-SG_OK_list[j][i] / gamma)
                for j in range(len(SG_OK_list)):
                    Prob_list[j][i] = math.exp(
                        -SG_OK_list[j][i] / gamma) / soma

            #Creating probabilities propreties
            for k in range(len(Prob_list)):
                prop_final_data_name = 'Probability_RT' + str(
                    RT[k].split('RT_')[-1])

                if (prop_final_data_name in lst_props_grid):
                    flag = 0
                    i = 1
                    while (flag == 0):
                        test_name = prop_final_data_name + '-' + str(i)
                        if (test_name not in lst_props_grid):
                            flag = 1
                            prop_final_data_name = test_name
                        i = i + 1

                sgems.set_property(grid_krig, prop_final_data_name,
                                   Prob_list[k])

            #Operating servo-system
            if self.params['servo_check']['value'] == '1':
                var_rt_grid = self.params['targe_prop']['grid']
                var_rt_st = self.params['targe_prop']['property']
                var_rt_region = self.params['targe_prop']['region']
                if len(var_rt_grid) == 0 or len(var_rt_st) == 0:
                    print 'Select the target proportion property'
                    return False

                #Getting variables
                var_rt = sgems.get_property(var_rt_grid, var_rt_st)

                #Getting parameters
                lambda1 = float(self.params['Lambda']['value'])
                mi = lambda1 / (1 - lambda1)

                #Checking if a region exist
                if len(var_rt_region) == 0:
                    #Variable without a region
                    var_region = var_rt

                else:
                    region_rt = sgems.get_region(var_rt_grid, var_rt_region)
                    #Geting the variable inside the region
                    var_region = []
                    for i in range(len(var_rt)):
                        if region_rt[i] == 1:
                            var_region.append(var_rt[i])

                #Getting the target proportion
                target_prop = proportion(var_region, RT)

                #Getting the random path
                ran_path = random_path(Prob_list[0])

                #Removing the blocks outside the region from randon path
                if len(var_rt_region) != 0:
                    for i in range(len(region_rt)):
                        if region_rt[i] == 0:
                            ran_path.remove(i)

                #servo system
                p = 0
                GeoModel_corrected = GeoModel[:]

                visited_rts = []
                for j in ran_path:
                    visited_rts.append(GeoModel[j])
                    instant_proportions = proportion(visited_rts, RT)

                    sgmax = 10e-21
                    for i in range(len(Prob_list)):
                        Prob_list[i][j] = Prob_list[i][j] + (
                            mi * (target_prop[i] - instant_proportions[i]))
                        if Prob_list[i][j] > sgmax:
                            sgmax = Prob_list[i][j]
                            p = i

                    GeoModel_corrected[j] = int(RT[p][-1])
                    visited_rts[-1] = int(RT[p].split('RT_')[-1])

                #Correcting servo servo-system by the biggest proportion on a neighborhood
                GeoModel_corrected_servo_prop = GeoModel_corrected[:]
                ran_path_servo_correction = random_path(
                    GeoModel_corrected_servo_prop)
                for i in ran_path_servo_correction:
                    vizinhanca = neighb(grid_krig, i)

                    blk_geo_model_corrected_servo = []
                    for j in vizinhanca:
                        blk_geo_model_corrected_servo.append(
                            GeoModel_corrected_servo_prop[j])

                    proportions_servo = proportion(
                        blk_geo_model_corrected_servo, RT)
                    indice_max_prop = proportions_servo.index(
                        max(proportions_servo))

                    GeoModel_corrected_servo_prop[i] = int(
                        RT[indice_max_prop].split('RT_')[-1])

                #Creating Geologic_Model_Servo_System property
                prop_final_data_name = 'Geologic_Model_Servo_System'

                if (prop_final_data_name in lst_props_grid):
                    flag = 0
                    i = 1
                    while (flag == 0):
                        test_name = prop_final_data_name + '-' + str(i)
                        if (test_name not in lst_props_grid):
                            flag = 1
                            prop_final_data_name = test_name
                        i = i + 1

                #Creating Geologic_Model_Corrected property
                prop_final_data_name1 = 'Geologic_Model_Corrected'

                if (prop_final_data_name1 in lst_props_grid):
                    flag = 0
                    i = 1
                    while (flag == 0):
                        test_name1 = prop_final_data_name1 + '-' + str(i)
                        if (test_name1 not in lst_props_grid):
                            flag = 1
                            prop_final_data_name1 = test_name1
                        i = i + 1

                #Assign conditioning data to grid node
                for i in range(len(RT_data)):
                    if not math.isnan(RT_data[i]):
                        closest_node = sgems.get_closest_nodeid(
                            grid_krig, X[i], Y[i], Z[i])
                        GeoModel_corrected[closest_node] = RT_data[i]
                        GeoModel_corrected_servo_prop[closest_node] = RT_data[
                            i]

                #Setting properties
                sgems.set_property(grid_krig, prop_final_data_name,
                                   GeoModel_corrected)
                sgems.set_property(grid_krig, prop_final_data_name1,
                                   GeoModel_corrected_servo_prop)

        return True
Beispiel #34
0
import sgems
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/FourItrMoreVariY.prj/velocity_grid (50x80x1)::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/well_data::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray1::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray2::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray3::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray4::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray5::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray6::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray7::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray8::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray9::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray10::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray11::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray12::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray13::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray14::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray15::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray16::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray17::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray18::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray19::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray20::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray21::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray22::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray23::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray24::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/MATLAB/Research/DataSets/RayRepo/Level4_4ray25::All')
sgems.execute('SaveGeostatGrid  velocity_grid (50x80x1)::C:/Users/noe/Documents/MATLAB/Research/DataSets/FourItrMoreVariY.prj/velocity_grid (50x80x1)::s-gems')
sgems.execute('SaveGeostatGrid  well_data::C:/Users/noe/Documents/MATLAB/Research/DataSets/FourItrMoreVariY.prj/well_data::s-gems')
sgems.execute('SaveGeostatGrid  Level4_4ray1::C:/Users/noe/Documents/MATLAB/Research/DataSets/FourItrMoreVariY.prj/Level4_4ray1::s-gems')
Beispiel #35
0
import sgems
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/LimitedBlock.prj/velocity_grid (40x40x1)::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/well_data::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/LimitedBlock.prj/All_StraightRays::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays1::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays2::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays3::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays4::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays5::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays6::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays7::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays8::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays9::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays10::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays11::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays12::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays13::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays14::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays15::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays16::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays17::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays18::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays19::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays20::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays21::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays22::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays23::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays24::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays25::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays26::All')
sgems.execute('LoadObjectFromFile C:/Users/noe/Documents/Research/DataSets/RayRepo/FirstLevelRays27::All')