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])
+ "::"
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
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
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"
)
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)
import sgems
sgems.execute("ReloadPythonPlugins")
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')
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)
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')
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
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')
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')
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> ')
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')
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')
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')
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>')
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
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')
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')
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" /> '+\
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
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" '+\
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')
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')
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]))
#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')
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
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')
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')