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')