def testInjectionWellCO2SmallWellR(self):
###
# Testing SemiAnalyticalWell for horizontal injection well settings
###
# load global physical properties
gpp = SimulationParameters(working_fluid = 'co2')
gpp.time_years = 10.
gpp.dT_dz = 0.06
gpp.well_radius = 0.02
gpp.m_dot_IP = 0.1
vertical_well = SemiAnalyticalWell(params = gpp,
dz_total = -3500.,
T_e_initial = 15.)
# initial state
initial_state = FluidState.getStateFromPT(1.e6, 25., gpp.working_fluid)
vertical_well_results = vertical_well.solve(initial_state)
self.assertMessages(gpp.working_fluid,
(vertical_well_results.state.P_Pa, 1.1431e6),
(vertical_well_results.state.T_C, 222.2246),
(vertical_well_results.state.h_Jkg, 6.8717e5))
def testInjectionWellCO2(self):
###
# Testing SemiAnalyticalWell for horizontal injection well settings
###
# load global physical properties
gpp = SimulationParameters(working_fluid = 'co2')
gpp.time_years = 10.
gpp.dT_dz = 0.06
gpp.well_radius = 0.279
gpp.m_dot_IP = 5.
vertical_well = SemiAnalyticalWell(params = gpp,
dz_total = -3500.,
T_e_initial = 15.)
# initial state
initial_state = FluidState.getStateFromPT(1.e6, 25., gpp.working_fluid)
vertical_well_results = vertical_well.solve(initial_state)
self.assertMessages(gpp.working_fluid,
(vertical_well_results.state.P_Pa, 1.7245e6),
(vertical_well_results.state.T_C, 156.08),
(vertical_well_results.state.h_Jkg, 6.1802e5))
horizontal_well = SemiAnalyticalWell(params = gpp,
dr_total = 3000.,
T_e_initial = vertical_well.T_e_initial + gpp.dT_dz * abs(vertical_well.dz_total))
# initial state
initial_state = FluidState.getStateFromPT(vertical_well_results.state.P_Pa, vertical_well_results.state.T_C, gpp.working_fluid)
horizontal_well_results = horizontal_well.solve(initial_state)
self.assertMessages(gpp.working_fluid,
(horizontal_well_results.state.P_Pa, 1.7235e6),
(horizontal_well_results.state.T_C, 212.746),
(horizontal_well_results.state.h_Jkg, 6.755e5))
def testInjectionWellWater(self):
###
# Testing SemiAnalyticalWell for vertical and horizontal injection well settings
###
# load parameters
gpp = SimulationParameters(working_fluid = 'water')
gpp.time_years = 10.
gpp.dT_dz = 0.06
gpp.well_radius = 0.279
gpp.m_dot_IP = 5.
vertical_well = SemiAnalyticalWell(params = gpp,
dz_total = -3500.,
T_e_initial = 15.)
# initial state
initial_state = FluidState.getStateFromPT(1.e6, 25., gpp.working_fluid)
vertical_well_results = vertical_well.solve(initial_state)
self.assertMessages(gpp.working_fluid,
(vertical_well_results.state.P_Pa, 3.533e7),
(vertical_well_results.state.T_C, 67.03),
(vertical_well_results.state.h_Jkg, 3.0963e5))
horizontal_well = SemiAnalyticalWell(params = gpp,
dr_total = 3000.,
T_e_initial = vertical_well.T_e_initial + gpp.dT_dz * abs(vertical_well.dz_total))
# initial state
initial_state = FluidState.getStateFromPT(vertical_well_results.state.P_Pa, vertical_well_results.state.T_C, gpp.working_fluid)
horizontal_well_results = horizontal_well.solve(initial_state)
self.assertMessages(gpp.working_fluid,
(horizontal_well_results.state.P_Pa, 3.533e7),
(horizontal_well_results.state.T_C, 121.99),
(horizontal_well_results.state.h_Jkg, 5.3712e5))