def test_volume_kw(self): grid = EclGrid(self.egrid_file()) vol = grid.createVolumeKeyword( ) self.assertEqual( len(vol) , grid.getNumActive()) for active_index , volume in enumerate(vol): self.assertEqual( volume , grid.cell_volume( active_index = active_index )) vol = grid.createVolumeKeyword( active_size = False ) self.assertEqual( len(vol) , grid.getGlobalSize()) for global_index , volume in enumerate(vol): self.assertEqual( volume , grid.cell_volume( global_index = global_index ))
class RegionDefinitionTest(ExtendedTestCase): def setUp(self): self.grid = EclGrid( self.createTestPath("Statoil/ECLIPSE/Mariner/MARINER.EGRID")) with open( self.createTestPath("Statoil/ECLIPSE/Mariner/faultblock.grdecl") ) as fileH: self.kw = EclKW.read_grdecl( fileH , "FAULTBLK" , ecl_type = EclTypeEnum.ECL_INT_TYPE ) self.faults = FaultCollection( self.grid , self.createTestPath("Statoil/ECLIPSE/Mariner/faults.grdecl")) self.poly_file1 = self.createTestPath("Statoil/ECLIPSE/Mariner/pol1.xyz") self.poly_file2 = self.createTestPath("Statoil/ECLIPSE/Mariner/pol2.xyz") self.poly_file3 = self.createTestPath("Statoil/ECLIPSE/Mariner/pol3.xyz") self.poly_file4 = self.createTestPath("Statoil/ECLIPSE/Mariner/pol4.xyz") self.poly_file5 = self.createTestPath("Statoil/ECLIPSE/Mariner/pol5.xyz") self.poly_file6 = self.createTestPath("Statoil/ECLIPSE/Mariner/pol6.xyz") self.poly_file7 = self.createTestPath("Statoil/ECLIPSE/Mariner/pol7.xyz") self.poly_file8 = self.createTestPath("Statoil/ECLIPSE/Mariner/pol8.xyz") self.poly_file9 = self.createTestPath("Statoil/ECLIPSE/Mariner/pol9.xyz") self.poly_file10 = self.createTestPath("Statoil/ECLIPSE/Mariner/pol10.xyz") self.poly_file11 = self.createTestPath("Statoil/ECLIPSE/Mariner/pol11.xyz") self.fault_blocks = [] for k in range(self.grid.getNZ()): blocks = FaultBlockLayer( self.grid , k) blocks.scanKeyword( self.kw ) self.fault_blocks.append( blocks ) def test_create(self): with self.assertRaises(TypeError): regionDef = RegionDefinition( "GG" ) with self.assertRaises(ValueError): regionDef = RegionDefinition.create( 1 , self.faults , ["Test" , 1]) regionDef = RegionDefinition.create( 1 , self.faults , ["DF28_LC" , self.poly_file1 , self.poly_file2 ]) def test_create(self): defRegion = RegionDefinition.create( 5 , self.faults , ["DF4_MC","DF15_MC","DF43_MC","DF25_LC","DF21_C","DF1_C",self.poly_file6, "DF26_MC","DF34_MCS","DF32_MC",self.poly_file5]) region_kw = EclKW.create( "REGIONS" , self.grid.getGlobalSize() , EclTypeEnum.ECL_INT_TYPE ) region_kw.assign( 0 ) for k in range(self.grid.getNZ()): with self.assertRaises(NotImplementedError): block_list = defRegion.findInternalBlocks( self.grid , self.fault_blocks[k] )
class RegionDefinitionTest(ExtendedTestCase): def setUp(self): self.grid = EclGrid( self.createTestPath("Statoil/ECLIPSE/Mariner/MARINER.EGRID")) with open( self.createTestPath( "Statoil/ECLIPSE/Mariner/faultblock.grdecl")) as fileH: self.kw = EclKW.read_grdecl(fileH, "FAULTBLK", ecl_type=EclTypeEnum.ECL_INT_TYPE) self.faults = FaultCollection( self.grid, self.createTestPath("Statoil/ECLIPSE/Mariner/faults.grdecl")) self.poly_file1 = self.createTestPath( "Statoil/ECLIPSE/Mariner/pol1.xyz") self.poly_file2 = self.createTestPath( "Statoil/ECLIPSE/Mariner/pol2.xyz") self.poly_file3 = self.createTestPath( "Statoil/ECLIPSE/Mariner/pol3.xyz") self.poly_file4 = self.createTestPath( "Statoil/ECLIPSE/Mariner/pol4.xyz") self.poly_file5 = self.createTestPath( "Statoil/ECLIPSE/Mariner/pol5.xyz") self.poly_file6 = self.createTestPath( "Statoil/ECLIPSE/Mariner/pol6.xyz") self.poly_file7 = self.createTestPath( "Statoil/ECLIPSE/Mariner/pol7.xyz") self.poly_file8 = self.createTestPath( "Statoil/ECLIPSE/Mariner/pol8.xyz") self.poly_file9 = self.createTestPath( "Statoil/ECLIPSE/Mariner/pol9.xyz") self.poly_file10 = self.createTestPath( "Statoil/ECLIPSE/Mariner/pol10.xyz") self.poly_file11 = self.createTestPath( "Statoil/ECLIPSE/Mariner/pol11.xyz") self.fault_blocks = [] for k in range(self.grid.getNZ()): blocks = FaultBlockLayer(self.grid, k) blocks.scanKeyword(self.kw) self.fault_blocks.append(blocks) def test_create(self): with self.assertRaises(TypeError): regionDef = RegionDefinition("GG") with self.assertRaises(ValueError): regionDef = RegionDefinition.create(1, self.faults, ["Test", 1]) regionDef = RegionDefinition.create( 1, self.faults, ["DF28_LC", self.poly_file1, self.poly_file2]) def test_create(self): defRegion = RegionDefinition.create(5, self.faults, [ "DF4_MC", "DF15_MC", "DF43_MC", "DF25_LC", "DF21_C", "DF1_C", self.poly_file6, "DF26_MC", "DF34_MCS", "DF32_MC", self.poly_file5 ]) region_kw = EclKW.create("REGIONS", self.grid.getGlobalSize(), EclTypeEnum.ECL_INT_TYPE) region_kw.assign(0) for k in range(self.grid.getNZ()): with self.assertRaises(NotImplementedError): block_list = defRegion.findInternalBlocks( self.grid, self.fault_blocks[k])