Skip to content

Commit

Permalink
update cost
Browse files Browse the repository at this point in the history
  • Loading branch information
@Zhuoran29
Zhuoran29 committed Jan 16, 2025
1 parent 0134f23 commit dce9c9d
Show file tree
Hide file tree
Showing 2 changed files with 38 additions and 33 deletions.
35 changes: 20 additions & 15 deletions src/watertap_contrib/reflo/analysis/case_studies/permian/permian_fo_cryst.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -561,7 +561,7 @@ def run_permian_FO(operating_condition,
permian_fo_config = {
"feed_vol_flow": 0.22, # initial value for fo model setup
"feed_TDS_mass": 0.119, # mass fraction, 0.119 is about 130 g/L
"recovery_ratio": 0.36,
"recovery_ratio": 0.41,
"RO_recovery_ratio":1, # RO recovery ratio
"NF_recovery_ratio":0.8, # Nanofiltration recovery ratio
"feed_temperature":25,
Expand Down Expand Up @@ -683,16 +683,16 @@ def run_permian_FO(operating_condition,
LCOWs =[]
failed= []

rr = [ 0.30, 0.32, 0.36, 0.38, 0.4,0.43, 0.45]
rr = [ 0.30, 0.32, 0.36, 0.38,0.43, 0.45]
strong_draw_mass = [i*0.03 + 0.80 for i in range(6)]
yields = [0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95]
nacl_prices = [0, -0.01, -0.02, -0.03,-0.04,-0.05,-0.06,-0.07,-0.08]

for v in nacl_prices:
for v in rr:
permian_fo_config = {
"feed_vol_flow": 0.22, # initial value for fo model setup
"feed_TDS_mass": 0.119, # mass fraction, 0.119 is about 130 g/L
"recovery_ratio": 0.45,
"recovery_ratio": v,
"RO_recovery_ratio":1, # RO recovery ratio
"NF_recovery_ratio":0.8, # Nanofiltration recovery ratio
"feed_temperature":25,
Expand All @@ -711,7 +711,7 @@ def run_permian_FO(operating_condition,
"nacl_yield": 0.9, # Yield
}
permian_cost_config = {
"nacl_recovery_price": v
"nacl_recovery_price": 0
}

try:
Expand Down Expand Up @@ -741,11 +741,11 @@ def run_permian_FO(operating_condition,

# fo_elec_cost = 0.07 * value(m.fs.treatment.FO.fs.fo.costing.electricity_flow) * 10290.711324821756
fo_heat_cost = 0.02 * value(m.fs.treatment.FO.fs.fo.costing.thermal_energy_flow) * 8766
chem_elec_cost = 0.07 * 10290.711324821756 * value(m.fs.treatment.costing._registered_flows["electricity"][0])
ec_elec_cost = 0.07 * 10290.711324821756 * value(m.fs.treatment.costing._registered_flows["electricity"][1])
filt_elec_cost = 0.07 * 10290.711324821756 * value(m.fs.treatment.costing._registered_flows["electricity"][2])
fo_elec_cost = 0.07 * 10290.711324821756 * value(m.fs.treatment.costing._registered_flows["electricity"][3])
cryst_elec_cost = 0.07 * sum([value(m.fs.treatment.costing._registered_flows['electricity'][i]) for i in [4,5,6,7]]) * 10290.711324821756
chem_elec_cost = 0.07 * 8766 * value(m.fs.treatment.costing._registered_flows["electricity"][0])
ec_elec_cost = 0.07 * 8766 * value(m.fs.treatment.costing._registered_flows["electricity"][1])
filt_elec_cost = 0.07 * 8766 * value(m.fs.treatment.costing._registered_flows["electricity"][2])
fo_elec_cost = 0.07 * 8766 * value(m.fs.treatment.costing._registered_flows["electricity"][3])
cryst_elec_cost = 0.07 * sum([value(m.fs.treatment.costing._registered_flows['electricity'][i]) for i in [4,5,6,7]]) * 8766
cryst_heat_cost = 0.02 * sum([value(m.fs.treatment.costing._registered_flows['heat'][i]) for i in [1]]) * 8766
cryst_nacl_revenue = -0.07 * value(m.fs.treatment.costing._registered_flows['NaCl_recovered'][0]) * 8766

Expand All @@ -767,6 +767,11 @@ def run_permian_FO(operating_condition,
fo_capexs.append(fo_capex*capital_recovery_rate/flow_vol)
cryst_capexs.append(cryst_capex*capital_recovery_rate/flow_vol)

# chem_opexs.append((chem_opex ) /flow_vol)
# ec_opexs.append( (ec_opex ) /flow_vol)
# filt_opexs.append((filt_opex )/flow_vol)
# fo_opexs.append( (fo_opex ) /flow_vol)
# cryst_opexs.append( (cryst_opex ) /flow_vol)
chem_opexs.append((chem_opex + h2o2_cost +chem_elec_cost) /flow_vol)
ec_opexs.append( (ec_opex + alum_cost +ec_elec_cost) /flow_vol)
filt_opexs.append((filt_opex + filt_elec_cost)/flow_vol)
Expand All @@ -783,8 +788,8 @@ def run_permian_FO(operating_condition,

#%% Make plots
import matplotlib.pyplot as plt

plt.stackplot(nacl_prices, NaCls,
# ec_opexs = [i*0.1 for i in ec_opexs]
plt.stackplot(rr, NaCls,
chem_capexs, chem_opexs,
ec_capexs, ec_opexs,
filt_capexs, filt_opexs,
Expand All @@ -796,7 +801,7 @@ def run_permian_FO(operating_condition,
'Cart filt CAPEX', 'Cart filt OPEX',
'FO CAPEX', 'FO OPEX',
'Cryst CAPEX', 'Cryst OPEX',
# 'Elec', 'Heat','Aluminum','H2O2'
'Elec', 'Heat','Aluminum','H2O2'
],
hatch =['','', '\\\\',
'', '\\\\',
Expand All @@ -818,10 +823,10 @@ def run_permian_FO(operating_condition,
plt.rcParams['figure.dpi']=300

# Show the legend
plt.legend(loc='upper center', bbox_to_anchor=(0.5, 1.15), ncol =2,prop={'size': 6})
plt.legend(loc='upper center', bbox_to_anchor=(0.5, 1.34), ncol =2,prop={'size': 8})

plt.ylabel('LCOW ($/m3)')
plt.xlabel('NaCl prices ($/kg)')
plt.xlabel('FO recovery rate')
plt.title('')
# Display the chart
plt.show()
Expand Down
36 changes: 18 additions & 18 deletions src/watertap_contrib/reflo/analysis/case_studies/permian/permian_fo_dwi.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -390,11 +390,11 @@ def run_permian_FO(operating_condition,

# scaling (based on grid participation), setup order
# deactivate constraints,
m.fs.costing = REFLOSystemCosting()
m.fs.costing.cost_process()
m.fs.costing.add_annual_water_production(flow_vol)
m.fs.costing.add_LCOW(flow_vol)
m.fs.costing.initialize()
# m.fs.costing = REFLOSystemCosting()
# m.fs.costing.cost_process()
# m.fs.costing.add_annual_water_production(flow_vol)
# m.fs.costing.add_LCOW(flow_vol)
# m.fs.costing.initialize()


print(f"DOF after add costing: {degrees_of_freedom(m)}")
Expand Down Expand Up @@ -455,12 +455,12 @@ def run_permian_FO(operating_condition,
h2o2_cost = value(m.fs.treatment.costing.aggregate_flow_costs["hydrogen_peroxide"])


chem_elec_cost = 0.07 * value(m.fs.treatment.chem_addition.unit.electricity[0]) * 10290.711324821756
ec_elec_cost = 0.07 * value(m.fs.treatment.ec.unit.costing.electricity_flow) * 10290.711324821756
filt_elec_cost = 0.07 * value(m.fs.treatment.cart_filt.unit.electricity[0]) * 10290.711324821756
fo_elec_cost = 0.07 * value(m.fs.treatment.FO.fs.fo.costing.electricity_flow) * 10290.711324821756
chem_elec_cost = 0.07 * value(m.fs.treatment.chem_addition.unit.electricity[0]) * 8766
ec_elec_cost = 0.07 * value(m.fs.treatment.ec.unit.costing.electricity_flow) * 8766
filt_elec_cost = 0.07 * value(m.fs.treatment.cart_filt.unit.electricity[0]) * 8766
fo_elec_cost = 0.07 * value(m.fs.treatment.FO.fs.fo.costing.electricity_flow) * 8766
fo_heat_cost = 0.02 * value(m.fs.treatment.FO.fs.fo.costing.thermal_energy_flow) * 8766
dwi_elec_cost = 0.07 * value(m.fs.treatment.DWI.costing.pumping_power_required) * 10290.711324821756
dwi_elec_cost = 0.07 * value(m.fs.treatment.DWI.costing.pumping_power_required) * 8766

print('LCOW ', lcow)

Expand Down Expand Up @@ -515,7 +515,7 @@ def run_permian_FO(operating_condition,
LCOWs =[]
failed= []

rr = [ 0.2,0.24,0.28,0.32,0.35,0.4,0.44, 0.45]
rr = [ 0.2,0.24,0.28,0.3,0.35,0.4,0.44, 0.45]
strong_draw_mass = [i*0.03 + 0.80 for i in range(6)]

for v in rr:
Expand All @@ -527,7 +527,7 @@ def run_permian_FO(operating_condition,
"NF_recovery_ratio":0.8, # Nanofiltration recovery ratio
"feed_temperature":25,
"strong_draw_temp":25, # Strong draw solution inlet temperature (C)
"strong_draw_mass_frac": 0.85, # Strong draw solution mass fraction
"strong_draw_mass_frac": 0.9, # Strong draw solution mass fraction
"product_draw_mass_frac": 0.01, # FO product draw solution mass fraction
"HX1_cold_out_temp": 78 + 273.15, # HX1 coldside outlet temperature
"HX1_hot_out_temp": 32 + 273.15, # HX1 hotside outlet temperature
Expand Down Expand Up @@ -561,11 +561,11 @@ def run_permian_FO(operating_condition,

# fo_elec_cost = 0.07 * value(m.fs.treatment.FO.fs.fo.costing.electricity_flow) * 10290.711324821756
fo_heat_cost = 0.02 * value(m.fs.treatment.FO.fs.fo.costing.thermal_energy_flow) * 8766
chem_elec_cost = 0.07 * 10290.711324821756 * value(m.fs.treatment.costing._registered_flows["electricity"][0])
ec_elec_cost = 0.07 * 10290.711324821756 * value(m.fs.treatment.costing._registered_flows["electricity"][1])
filt_elec_cost = 0.07 * 10290.711324821756 * value(m.fs.treatment.costing._registered_flows["electricity"][2])
fo_elec_cost = 0.07 * 10290.711324821756 * value(m.fs.treatment.costing._registered_flows["electricity"][3])
dwi_elec_cost = 0.07 * 10290.711324821756 * value(m.fs.treatment.costing._registered_flows["electricity"][4])
chem_elec_cost = 0.07 * 8766 * value(m.fs.treatment.costing._registered_flows["electricity"][0])
ec_elec_cost = 0.07 * 8766 * value(m.fs.treatment.costing._registered_flows["electricity"][1])
filt_elec_cost = 0.07 * 8766 * value(m.fs.treatment.costing._registered_flows["electricity"][2])
fo_elec_cost = 0.07 * 8766 * value(m.fs.treatment.costing._registered_flows["electricity"][3])
dwi_elec_cost = 0.07 * 8766 * value(m.fs.treatment.costing._registered_flows["electricity"][4])


var_opex_total = value(m.fs.treatment.costing.total_variable_operating_cost)
Expand Down Expand Up @@ -635,7 +635,7 @@ def run_permian_FO(operating_condition,
plt.rcParams['figure.dpi']=300

# Show the legend
plt.legend(loc='upper right', ncol =2,prop={'size': 6})
plt.legend(loc='upper right', bbox_to_anchor=(0.5, 1.34), ncol =2,prop={'size': 8})

plt.ylabel('LCOW ($/m3)')
plt.xlabel('FO recovery rate')
Expand Down

0 comments on commit dce9c9d

Please sign in to comment.