"BHP" : 75*U.psi,
"schedule" : [0.*U.yr, 2*U.yr]
}
}
# print input
print(("<%s> Execution started."%time.asctime()))
DIM=2
domain=Rectangle(N_X, N_Y, l0=L_X, l1=L_Y)
N=1000
for I in wellspecs:
domain.setTagMap(I, N)
N+=1
domain.addDiracPoint(wellspecs[I]["X0"][:DIM], I)
print(("<%s> Well %s introduced to domain."%(time.asctime(), I)))
#domain=Brick(N_X, N_Y,N_Z,l0=L_X, l1=L_Y,l2=L_Z)
print(("<%s> Domain has been generated."%time.asctime()))
print("length x-direction = %f km"%(sup(domain.getX()[0])/U.km))
print("cell size in x direction = %f m"%(CELL_X/U.m))
print("length y-direction = %f km"%(sup(domain.getX()[1])/U.km))
print("cell size in y direction = %f m"%(CELL_Y/U.m))
print("fracture permeability in x direction= %f mD"%(PERM_F_X/(U.mDarcy)))
print("fracture permeability in y direction= %f mD"%(PERM_F_Y/(U.mDarcy)))
print("fracture permeability in z direction= %f mD"%(PERM_F_Z/(U.mDarcy)))
"Q": [0., 2000 * U.Barrel / U.day],
"BHP": 75 * U.psi,
"schedule": [0. * U.yr, 2 * U.yr]
}
}
# print input
print(("<%s> Execution started." % time.asctime()))
DIM = 2
domain = Rectangle(N_X, N_Y, l0=L_X, l1=L_Y)
N = 1000
for I in wellspecs:
domain.setTagMap(I, N)
N += 1
domain.addDiracPoint(wellspecs[I]["X0"][:DIM], I)
print(("<%s> Well %s introduced to domain." % (time.asctime(), I)))
#domain=Brick(N_X, N_Y,N_Z,l0=L_X, l1=L_Y,l2=L_Z)
print(("<%s> Domain has been generated." % time.asctime()))
print("length x-direction = %f km" % (sup(domain.getX()[0]) / U.km))
print("cell size in x direction = %f m" % (CELL_X / U.m))
print("length y-direction = %f km" % (sup(domain.getX()[1]) / U.km))
print("cell size in y direction = %f m" % (CELL_Y / U.m))
print("fracture permeability in x direction= %f mD" % (PERM_F_X / (U.mDarcy)))
print("fracture permeability in y direction= %f mD" % (PERM_F_Y / (U.mDarcy)))
print("fracture permeability in z direction= %f mD" % (PERM_F_Z / (U.mDarcy)))
mkDir(OUTPUT_DIR)
