Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Add drilling cost curves to geothermal cost cmod #1191

Merged
merged 7 commits into from
Aug 1, 2024
Merged

Add drilling cost curves to geothermal cost cmod #1191

Show file tree
Hide file tree
Changes from 5 commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
1b916e8
Check for calculation of negative prodcution wells in well requirement
mjprilliman Jun 24, 2024
042f6b1
Add drilling cost curves to geothermal cost cmod, report production and
mjprilliman Jun 25, 2024
17eacb8
Merge branch 'patch' into drilling-cost-ssc-equation
mjprilliman Jul 22, 2024
fa33583
Add drilling cost calculations as optional inputs to cmod geothermal
mjprilliman Jul 22, 2024
369aeff
Merge branch 'patch' into drilling-cost-ssc-equation
mjprilliman Jul 26, 2024
184f478
Update vartable dependencies, add section breaks to switch cases
mjprilliman Jul 31, 2024
406acf0
Merge branch 'patch' into drilling-cost-ssc-equation
mjprilliman Jul 31, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
1 change: 1 addition & 0 deletions shared/lib_geothermal.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1567,6 +1567,7 @@ double CGeothermalAnalyzer::GetNumberOfWells(void)
}

if (mp_geo_out->md_NumberOfWellsInj < 0) mp_geo_out->md_NumberOfWellsInj = 0;
if (mp_geo_out->md_NumberOfWells < 0) mp_geo_out->md_NumberOfWells = 0;
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Both injection and regular wells drilled - good catch!

}

return mp_geo_out->md_NumberOfWells;
Expand Down
168 changes: 166 additions & 2 deletions ssc/cmod_geothermal_costs.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,9 +75,26 @@ static var_info _cm_vtab_geothermal_costs[] = {
{ SSC_INPUT, SSC_NUMBER, "hp_flash_pressure", "HP Flash Pressure", "psia", "", "GeoHourly", "conversion_type=1", "", "" },
{ SSC_INPUT, SSC_NUMBER, "lp_flash_pressure", "LP Flash Pressure", "psia", "", "GeoHourly", "conversion_type=1", "", "" },
{ SSC_INPUT, SSC_NUMBER, "flash_count", "Flash Count", "(1 -2)", "", "GeoHourly", "conversion_type=1", "", "" },
// Outputs

{ SSC_INPUT, SSC_NUMBER, "geotherm.cost.inj_cost_curve_welltype", "Injection well type", "0/1", "", "GeoHourly", "", "", "" },
{ SSC_INPUT, SSC_NUMBER, "geotherm.cost.prod_cost_curve_welltype", "Production well type", "0/1", "", "GeoHourly", "", "", "" },
{ SSC_INPUT, SSC_NUMBER, "geotherm.cost.inj_cost_curve_welldiam", "Injection well diameter type", "0/1", "0=LargerDiameter,1=SmallerDiameter", "GeoHourly", "", "", "" },
{ SSC_INPUT, SSC_NUMBER, "geotherm.cost.prod_cost_curve_welldiam", "Production well diameter type", "0/1", "0=LargerDiameter,1=SmallerDiameter", "GeoHourly", "", "", "" },
{ SSC_INPUT, SSC_NUMBER, "geotherm.cost.inj_cost_curve", "Injection well diameter type", "0/1", "0=LargerDiameter,1=SmallerDiameter", "GeoHourly", "", "", "" },
{ SSC_INPUT, SSC_NUMBER, "geotherm.cost.prod_cost_curve", "Production well diameter type", "0/1", "0=LargerDiameter,1=SmallerDiameter", "GeoHourly", "", "", "" },
{ SSC_INPUT, SSC_NUMBER, "resource_depth", "Resource Depth", "m", "", "GeoHourly", "", "", "" },
{ SSC_INPUT, SSC_NUMBER, "geotherm.cost.prod_wells_drilled", "Number of drilled production wells", "0/1", "0=LargerDiameter,1=SmallerDiameter", "GeoHourly", "", "", "" },
{ SSC_INPUT, SSC_NUMBER, "geotherm.cost.inj_wells_drilled", "Number of drilled injection wells", "0/1", "0=LargerDiameter,1=SmallerDiameter", "GeoHourly", "", "", "" },
{ SSC_INPUT, SSC_NUMBER, "calc_drill_costs", "Calculate drill costs", "0/1", "0=LargerDiameter,1=SmallerDiameter", "GeoHourly", "?=0", "", "" },
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

curious why you are using "." instead of "_"?

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

We've just been working off what was there before. It's probably worth changing this in the future.



// Outputs

{ SSC_OUTPUT, SSC_NUMBER, "baseline_cost", "Baseline Cost", "$/kW", "", "GeoHourly", "?", "", "" },
{ SSC_OUTPUT, SSC_NUMBER, "inj_total_cost", "Total Injection well cost", "$/kW", "", "GeoHourly", "?", "", "" },
{ SSC_OUTPUT, SSC_NUMBER, "prod_total_cost", "Total Production well cost", "$/kW", "", "GeoHourly", "?", "", "" },


var_info_invalid };


Expand Down Expand Up @@ -295,6 +312,146 @@ class cm_geothermal_costs : public compute_module
void exec() override
{
SGeothermal_Inputs geo_inputs;

int calc_drill_costs = as_integer("calc_drill_costs");
if (calc_drill_costs == 1) {
int inj_cost_curve_welldiam = as_integer("geotherm.cost.inj_cost_curve_welldiam");
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

shouldn't these inputs be required if "calc_drill_costs"==1?

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Yes, I will update this.

int inj_cost_curve_welltype = as_integer("geotherm.cost.inj_cost_curve_welltype");
int inj_cost_curve = as_integer("geotherm.cost.inj_cost_curve");
int prod_cost_curve_welldiam = as_integer("geotherm.cost.prod_cost_curve_welldiam");
int prod_cost_curve_welltype = as_integer("geotherm.cost.prod_cost_curve_welltype");
int prod_cost_curve = as_integer("geotherm.cost.prod_cost_curve");

//Drilling cost calculations
double resource_depth = as_double("resource_depth");
double inj_well_cost = 0;
if (inj_cost_curve_welldiam == 0) {
if (inj_cost_curve_welltype == 0) {
switch (inj_cost_curve) {
case 0:
inj_well_cost = 0.281801107 * pow(resource_depth, 2) + 1275.521301 * resource_depth + 632315.1264;
case 1:
inj_well_cost = 0.189267288 * pow(resource_depth, 2) + 293.4517365 * resource_depth + 1326526.313;
case 2:
inj_well_cost = 0.003145418 * pow(resource_depth, 2) + 782.70 * resource_depth + 983620.25;
case 3:
inj_well_cost = -0.002397497 * pow(resource_depth, 2) + 752.94 * resource_depth + 524337.65;
}
}
else {
switch (inj_cost_curve) {
case 0:
inj_well_cost = 0.2528 * pow(resource_depth, 2) + 1716.72 * resource_depth + 500866.89;
case 1:
inj_well_cost = 0.19950 * pow(resource_depth, 2) + 296.13 * resource_depth + 1697867.71;
case 2:
inj_well_cost = 0.0038019 * pow(resource_depth, 2) + 838.90 * resource_depth + 1181947.04;
case 3:
inj_well_cost = 0.0037570 * pow(resource_depth, 2) + 762.53 * resource_depth + 765103.08;
}

}
}
else {
if (inj_cost_curve_welltype == 0) {
switch (inj_cost_curve) {
case 0:
inj_well_cost = 0.30212 * pow(resource_depth, 2) + 584.91 * resource_depth + 751368.47;
case 1:
inj_well_cost = 0.13710 * pow(resource_depth, 2) + 129.61 * resource_depth + 1205587.57;
case 2:
inj_well_cost = 0.0080395 * pow(resource_depth, 2) + 455.61 * resource_depth + 921007.69;
case 3:
inj_well_cost = 0.0025212 * pow(resource_depth, 2) + 439.45 * resource_depth + 590611.90;
}

}
else {
switch (inj_cost_curve) {
case 0:
inj_well_cost = 0.28977 * pow(resource_depth, 2) + 882.15 * resource_depth + 680562.50;
case 1:
inj_well_cost = 0.15340 * pow(resource_depth, 2) + 120.32 * resource_depth + 1431801.54;
case 2:
inj_well_cost = 0.0085389 * pow(resource_depth, 2) + 506.08 * resource_depth + 1057330.39;
case 3:
inj_well_cost = 0.0071869 * pow(resource_depth, 2) + 455.85 * resource_depth + 753377.73;
}

}
}
double inj_wells_drilled = as_double("geotherm.cost.inj_wells_drilled");
double inj_total_cost = inj_wells_drilled * inj_well_cost;
assign("inj_total_cost", inj_total_cost);

double prod_well_cost = 0;
if (prod_cost_curve_welldiam == 0) {
if (prod_cost_curve_welltype == 0) {
switch (prod_cost_curve) {
case 0:
prod_well_cost = 0.281801107 * pow(resource_depth, 2) + 1275.521301 * resource_depth + 632315.1264;
case 1:
prod_well_cost = 0.189267288 * pow(resource_depth, 2) + 293.4517365 * resource_depth + 1326526.313;
case 2:
prod_well_cost = 0.003145418 * pow(resource_depth, 2) + 782.70 * resource_depth + 983620.25;
case 3:
prod_well_cost = -0.002397497 * pow(resource_depth, 2) + 752.94 * resource_depth + 524337.65;
}
}
else {
switch (prod_cost_curve) {
case 0:
prod_well_cost = 0.2528 * pow(resource_depth, 2) + 1716.72 * resource_depth + 500866.89;
case 1:
prod_well_cost = 0.19950 * pow(resource_depth, 2) + 296.13 * resource_depth + 1697867.71;
case 2:
prod_well_cost = 0.0038019 * pow(resource_depth, 2) + 838.90 * resource_depth + 1181947.04;
case 3:
prod_well_cost = 0.0037570 * pow(resource_depth, 2) + 762.53 * resource_depth + 765103.08;
}

}
}
else {
if (prod_cost_curve_welltype == 0) {
switch (prod_cost_curve) {
case 0:
prod_well_cost = 0.30212 * pow(resource_depth, 2) + 584.91 * resource_depth + 751368.47;
case 1:
prod_well_cost = 0.13710 * pow(resource_depth, 2) + 129.61 * resource_depth + 1205587.57;
case 2:
prod_well_cost = 0.0080395 * pow(resource_depth, 2) + 455.61 * resource_depth + 921007.69;
case 3:
prod_well_cost = 0.0025212 * pow(resource_depth, 2) + 439.45 * resource_depth + 590611.90;
}

}
else {
switch (prod_cost_curve) {
case 0:
prod_well_cost = 0.28977 * pow(resource_depth, 2) + 882.15 * resource_depth + 680562.50;
case 1:
prod_well_cost = 0.15340 * pow(resource_depth, 2) + 120.32 * resource_depth + 1431801.54;
case 2:
prod_well_cost = 0.0085389 * pow(resource_depth, 2) + 506.08 * resource_depth + 1057330.39;
case 3:
prod_well_cost = 0.0071869 * pow(resource_depth, 2) + 455.85 * resource_depth + 753377.73;
}

}
}

double prod_wells_drilled = as_double("geotherm.cost.prod_wells_drilled");
double prod_total_cost = prod_wells_drilled * prod_well_cost;
assign("prod_total_cost", prod_total_cost);

if (!is_assigned("conversion_type")) {
return;
}
else {
//keep going
}
}
int conversion_type = as_integer("conversion_type");

int ppi_base_year = as_integer("ppi_base_year");
Expand Down Expand Up @@ -569,6 +726,12 @@ class cm_geothermal_costs : public compute_module

}







//OM Cost calculations
/*
double unit_plant = as_double("gross_output");
Expand Down Expand Up @@ -603,12 +766,13 @@ class cm_geothermal_costs : public compute_module
double pump_om_cost = total_annual_pump_cost + oil_downhole_pump_cost;
double field_om = well_om + surface_equip_om + chemical_cost + makeup_water_cost + pump_om_cost + field_labor_om;


//Annual Tax and Insurance
double total_capital_cost = as_double("total_installed_cost"); //todo move OM to separate cmod, pass in total installed cost
double tax_insurance_rate = 0.0075;
double annual_tax_insurance_cost = total_capital_cost * tax_insurance_rate;
*/

};

};
Expand Down