Skip to content

Commit

Permalink
Merge branch 'well-drilling-cost-correlation-consolidation'
Browse files Browse the repository at this point in the history
  • Loading branch information
@softwareengineerprogrammer
softwareengineerprogrammer committed Jun 7, 2024
2 parents 4ba73e8 + ce3fd0c commit 25553d4
Show file tree
Hide file tree
Showing 3 changed files with 50 additions and 137 deletions.
90 changes: 8 additions & 82 deletions src/geophires_x/Economics.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -95,43 +95,13 @@ def calculate_cost_of_one_vertical_well(model: Model, depth_m: float, well_corre
)

if well_correlation is WellDrillingCostCorrelation.SIMPLE:
cost_of_one_well = vertical_drilling_cost_per_m * depth_m * 1E-6

elif well_correlation is WellDrillingCostCorrelation.VERTICAL_SMALL:
cost_of_one_well = (0.30212 * depth_m ** 2 + 584.91124 * depth_m + 751368.47270) * 1E-6
elif well_correlation is WellDrillingCostCorrelation.VERTICAL_LARGE:
cost_of_one_well = (0.28180 * depth_m ** 2 + 1275.52130 * depth_m + 632315.12640) * 1E-6
elif well_correlation is WellDrillingCostCorrelation.DEVIATED_SMALL:
cost_of_one_well = (0.28977 * depth_m ** 2 + 882.15067 * depth_m + 680562.50150) * 1E-6
elif well_correlation is WellDrillingCostCorrelation.DEVIATED_LARGE:
cost_of_one_well = (0.25528 * depth_m ** 2 + 1716.71568 * depth_m + 500866.89110) * 1E-6
elif well_correlation is WellDrillingCostCorrelation.VERTICAL_SMALL_INT1:
cost_of_one_well = (0.13710 * depth_m ** 2 + 129.61033 * depth_m + 1205587.57100) * 1E-6
elif well_correlation is WellDrillingCostCorrelation.VERTICAL_LARGE_INT1:
cost_of_one_well = (0.18927 * depth_m ** 2 + 293.45174 * depth_m + 1326526.31300) * 1E-6
elif well_correlation is WellDrillingCostCorrelation.DEVIATED_SMALL_INT1:
cost_of_one_well = (0.15340 * depth_m ** 2 + 120.31700 * depth_m + 1431801.54400) * 1E-6
elif well_correlation is WellDrillingCostCorrelation.DEVIATED_LARGE_INT1:
cost_of_one_well = (0.19950 * depth_m ** 2 + 296.13011 * depth_m + 1697867.70900) * 1E-6
elif well_correlation is WellDrillingCostCorrelation.VERTICAL_SMALL_INT2:
cost_of_one_well = (0.00804 * depth_m ** 2 + 455.60507 * depth_m + 921007.68680) * 1E-6
elif well_correlation is WellDrillingCostCorrelation.VERTICAL_LARGE_INT2:
cost_of_one_well = (0.00315 * depth_m ** 2 + 782.69676 * depth_m + 983620.25270) * 1E-6
elif well_correlation is WellDrillingCostCorrelation.DEVIATED_SMALL_INT2:
cost_of_one_well = (0.00854 * depth_m ** 2 + 506.08357 * depth_m + 1057330.39000) * 1E-6
elif well_correlation is WellDrillingCostCorrelation.DEVIATED_LARGE_INT2:
cost_of_one_well = (0.00380 * depth_m ** 2 + 838.90249 * depth_m + 1181947.04400) * 1E-6
elif well_correlation is WellDrillingCostCorrelation.VERTICAL_SMALL_IDEAL:
cost_of_one_well = (0.00252 * depth_m ** 2 + 439.44503 * depth_m + 590611.90110) * 1E-6
elif well_correlation is WellDrillingCostCorrelation.VERTICAL_LARGE_IDEAL:
cost_of_one_well = (-0.00240 * depth_m ** 2 + 752.93946 * depth_m + 524337.65380) * 1E-6
elif well_correlation is WellDrillingCostCorrelation.DEVIATED_SMALL_IDEAL:
cost_of_one_well = (0.00719 * depth_m ** 2 + 455.85233 * depth_m + 753377.73080) * 1E-6
elif well_correlation is WellDrillingCostCorrelation.DEVIATED_LARGE_IDEAL:
cost_of_one_well = (0.00376 * depth_m ** 2 + 762.52696 * depth_m + 765103.07690) * 1E-6
cost_of_one_well = vertical_drilling_cost_per_m * depth_m * 1E-6
else:
cost_of_one_well = well_correlation.calculate_cost_MUSD(depth_m)

# account for adjustment factor
cost_of_one_well = well_cost_adjustment_factor * cost_of_one_well

return cost_of_one_well


Expand Down Expand Up @@ -195,60 +165,16 @@ def calculate_cost_of_non_vertical_section(model: Model, length_m: float, well_c
if not NonverticalsCased:
# assume that casing & cementing costs 50% of drilling costs
casing_factor = 0.5

if model.economics.Nonvertical_drilling_cost_per_m.Provided or well_correlation is WellDrillingCostCorrelation.SIMPLE:
cost_of_non_vertical_section = casing_factor * ((num_nonvertical_sections * nonvertical_drilling_cost_per_m * length_per_section_m)) * 1E-6
elif well_correlation is WellDrillingCostCorrelation.VERTICAL_SMALL:
cost_per_section = (0.30212 * length_per_section_m ** 2 + 584.91124 * length_per_section_m + 751368.47270) * 1E-6
cost_of_non_vertical_section = casing_factor * num_nonvertical_sections * cost_per_section
elif well_correlation is WellDrillingCostCorrelation.VERTICAL_LARGE:
cost_per_section = (0.28180 * length_per_section_m ** 2 + 1275.52130 * length_per_section_m + 632315.12640) * 1E-6
cost_of_non_vertical_section = casing_factor * num_nonvertical_sections * cost_per_section
elif well_correlation is WellDrillingCostCorrelation.DEVIATED_SMALL:
cost_per_section = (0.28977 * length_per_section_m ** 2 + 882.15067 * length_per_section_m + 680562.50150) * 1E-6
cost_of_non_vertical_section = casing_factor * num_nonvertical_sections * cost_per_section
elif well_correlation is WellDrillingCostCorrelation.DEVIATED_LARGE:
cost_per_section = (0.25528 * length_per_section_m ** 2 + 1716.71568 * length_per_section_m + 500866.89110) * 1E-6
cost_of_non_vertical_section = casing_factor * num_nonvertical_sections * cost_per_section
elif well_correlation is WellDrillingCostCorrelation.VERTICAL_SMALL_INT1:
cost_per_section = (0.13710 * length_per_section_m ** 2 + 129.61033 * length_per_section_m + 1205587.57100) * 1E-6
cost_of_non_vertical_section = casing_factor * num_nonvertical_sections * cost_per_section
elif well_correlation is WellDrillingCostCorrelation.VERTICAL_LARGE_INT1:
cost_per_section = (0.18927 * length_per_section_m ** 2 + 293.45174 * length_per_section_m + 1326526.31300) * 1E-6
cost_of_non_vertical_section = casing_factor * num_nonvertical_sections * cost_per_section
elif well_correlation is WellDrillingCostCorrelation.DEVIATED_SMALL_INT1:
cost_per_section = (0.15340 * length_per_section_m ** 2 + 120.31700 * length_per_section_m + 1431801.54400) * 1E-6
cost_of_non_vertical_section = casing_factor * num_nonvertical_sections * cost_per_section
elif well_correlation is WellDrillingCostCorrelation.DEVIATED_LARGE_INT1:
cost_per_section = (0.19950 * length_per_section_m ** 2 + 296.13011 * length_per_section_m + 1697867.70900) * 1E-6
cost_of_non_vertical_section = casing_factor * num_nonvertical_sections * cost_per_section
elif well_correlation is WellDrillingCostCorrelation.VERTICAL_SMALL_INT2:
cost_per_section = (0.00804 * length_per_section_m ** 2 + 455.60507 * length_per_section_m + 921007.68680) * 1E-6
cost_of_non_vertical_section = casing_factor * num_nonvertical_sections * cost_per_section
elif well_correlation is WellDrillingCostCorrelation.VERTICAL_LARGE_INT2:
cost_per_section = (0.00315 * length_per_section_m ** 2 + 782.69676 * length_per_section_m + 983620.25270) * 1E-6
cost_of_non_vertical_section = casing_factor * num_nonvertical_sections * cost_per_section
elif well_correlation is WellDrillingCostCorrelation.DEVIATED_SMALL_INT2:
cost_per_section = (0.00854 * length_per_section_m ** 2 + 506.08357 * length_per_section_m + 1057330.39000) * 1E-6
cost_of_non_vertical_section = casing_factor * num_nonvertical_sections * cost_per_section
elif well_correlation is WellDrillingCostCorrelation.DEVIATED_LARGE_INT2:
cost_per_section = (0.00380 * length_per_section_m ** 2 + 838.90249 * length_per_section_m + 1181947.04400) * 1E-6
cost_of_non_vertical_section = casing_factor * num_nonvertical_sections * cost_per_section
elif well_correlation is WellDrillingCostCorrelation.VERTICAL_SMALL_IDEAL:
cost_per_section = (0.00252 * length_per_section_m ** 2 + 439.44503 * length_per_section_m + 590611.90110) * 1E-6
cost_of_non_vertical_section = casing_factor * num_nonvertical_sections * cost_per_section
elif well_correlation is WellDrillingCostCorrelation.VERTICAL_LARGE_IDEAL:
cost_per_section = (-0.00240 * length_per_section_m ** 2 + 752.93946 * length_per_section_m + 524337.65380) * 1E-6
cost_of_non_vertical_section = casing_factor * num_nonvertical_sections * cost_per_section
elif well_correlation is WellDrillingCostCorrelation.DEVIATED_SMALL_IDEAL:
cost_per_section = (0.00719 * length_per_section_m ** 2 + 455.85233 * length_per_section_m + 753377.73080) * 1E-6
cost_of_non_vertical_section = casing_factor * num_nonvertical_sections * cost_per_section
elif well_correlation is WellDrillingCostCorrelation.DEVIATED_LARGE_IDEAL:
cost_per_section = (0.00376 * length_per_section_m ** 2 + 762.52696 * length_per_section_m + 765103.07690) * 1E-6
else:
cost_per_section = well_correlation.calculate_cost_MUSD(length_per_section_m)
cost_of_non_vertical_section = casing_factor * num_nonvertical_sections * cost_per_section

# account for adjustment factor

cost_of_non_vertical_section = well_cost_adjustment_factor * cost_of_non_vertical_section

return cost_of_non_vertical_section


Expand Down
57 changes: 36 additions & 21 deletions src/geophires_x/OptionList.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -51,27 +51,42 @@ class ReservoirVolume(str, Enum):


class WellDrillingCostCorrelation(str, Enum):
"""Note: order must be retained since input is read as an int"""

VERTICAL_SMALL = "vertical small diameter, baseline"
DEVIATED_SMALL = "deviated small diameter, baseline"
VERTICAL_LARGE = "vertical large diameter, baseline"
DEVIATED_LARGE = "deviated large diameter, baseline"

SIMPLE = "Simple"

VERTICAL_SMALL_INT1 = "vertical small diameter, intermediate1"
VERTICAL_SMALL_INT2 = "vertical small diameter, intermediate2"
DEVIATED_SMALL_INT1 = "deviated small diameter, intermediate1"
DEVIATED_SMALL_INT2 = "deviated small diameter, intermediate2"
VERTICAL_LARGE_INT1 = "vertical large diameter, intermediate1"
VERTICAL_LARGE_INT2 = "vertical large diameter, intermediate2"
DEVIATED_LARGE_INT1 = "deviated large diameter, intermediate1"
DEVIATED_LARGE_INT2 = "deviated large diameter, intermediate2"
VERTICAL_SMALL_IDEAL = "vertical open-hole, small diameter, ideal"
DEVIATED_SMALL_IDEAL = "deviated liner, small diameter, ideal"
VERTICAL_LARGE_IDEAL = "vertical open-hole, large diameter, ideal"
DEVIATED_LARGE_IDEAL = "deviated liner, large diameter, ideal"
"""Note: order must be retained since input is read as an int; first int arg is duplicative of order"""

VERTICAL_SMALL = 1, "vertical small diameter, baseline", 0.30212, 584.91124, 751368.47270
DEVIATED_SMALL = 2, "deviated small diameter, baseline", 0.28977, 882.15067, 680562.50150
VERTICAL_LARGE = 3, "vertical large diameter, baseline", 0.28180, 1275.52130, 632315.12640
DEVIATED_LARGE = 4, "deviated large diameter, baseline", 0.25528, 1716.71568, 500866.89110

SIMPLE = 5, "Simple", 0, 1846*1E6, 0 # Based on Fervo Project Cape cost per meter (~$1846/m)

VERTICAL_SMALL_INT1 = 6, "vertical small diameter, intermediate1", 0.13710, 129.61033, 1205587.57100
VERTICAL_SMALL_INT2 = 7, "vertical small diameter, intermediate2", 0.00804, 455.60507, 921007.68680
DEVIATED_SMALL_INT1 = 8, "deviated small diameter, intermediate1", 0.15340, 120.31700, 1431801.54400
DEVIATED_SMALL_INT2 = 9, "deviated small diameter, intermediate2", 0.00854, 506.08357, 1057330.39000
VERTICAL_LARGE_INT1 = 10, "vertical large diameter, intermediate1", 0.18927, 293.45174, 1326526.31300
VERTICAL_LARGE_INT2 = 11, "vertical large diameter, intermediate2", 0.00315, 782.69676, 983620.25270
DEVIATED_LARGE_INT1 = 12, "deviated large diameter, intermediate1", 0.19950, 296.13011, 1697867.70900
DEVIATED_LARGE_INT2 = 13, "deviated large diameter, intermediate2", 0.00380, 838.90249, 1181947.04400
VERTICAL_SMALL_IDEAL = 14, "vertical open-hole, small diameter, ideal", 0.00252, 439.44503, 590611.90110
DEVIATED_SMALL_IDEAL = 15, "deviated liner, small diameter, ideal", 0.00719, 455.85233, 753377.73080
VERTICAL_LARGE_IDEAL = 16, "vertical open-hole, large diameter, ideal", -0.00240, 752.93946, 524337.65380
DEVIATED_LARGE_IDEAL = 17, "deviated liner, large diameter, ideal", 0.00376, 762.52696, 765103.07690


def __new__(cls, *args, **kwds):
obj = str.__new__(cls)
obj._value_ = args[1]
return obj

# ignore the first param since it's already set by __new__
def __init__(self, idx: int, _: str, c2:float, c1:float, c0:float):
self._c2 = c2
self._c1 = c1
self._c0 = c0

def calculate_cost_MUSD(self, meters) -> float:
return (self._c2 * meters ** 2 + self._c1 * meters + self._c0) * 1E-6



Expand Down
40 changes: 6 additions & 34 deletions src/geophires_x/SUTRAEconomics.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -424,48 +424,20 @@ def Calculate(self, model: Model) -> None:

# CAPEX
# Drilling

self.C1well = 0
if self.ccwellfixed.Valid:
self.C1well = self.ccwellfixed.value
self.Cwell.value = self.C1well * (model.wellbores.nprod.value + model.wellbores.ninj.value)
else:
if model.reserv.depth.value > 7000.0 or model.reserv.depth.value < 500:
print("Warning: simple user-specified cost per meter used for drilling depth < 500 or > 7000 m")
print('Warning: simple user-specified cost per meter used for drilling depth < 500 or > 7000 m')
model.logger.warning(
"Warning: simple user-specified cost per meter used for drilling depth < 500 or > 7000 m"
'Warning: simple user-specified cost per meter used for drilling depth < 500 or > 7000 m'
)
if self.wellcorrelation.value == WellDrillingCostCorrelation.VERTICAL_SMALL:
self.C1well = (0.30212 * model.reserv.depth.value ** 2 + 584.91124 * model.reserv.depth.value + 751368.47270) * 1E-6
elif self.wellcorrelation.value == WellDrillingCostCorrelation.VERTICAL_LARGE:
self.C1well = (0.28180 * model.reserv.depth.value ** 2 + 1275.52130 * model.reserv.depth.value + 632315.12640) * 1E-6
elif self.wellcorrelation.value == WellDrillingCostCorrelation.DEVIATED_SMALL:
self.C1well = (0.28977 * model.reserv.depth.value ** 2 + 882.15067 * model.reserv.depth.value + 680562.50150) * 1E-6
elif self.wellcorrelation.value == WellDrillingCostCorrelation.DEVIATED_LARGE:
self.C1well = (0.25528 * model.reserv.depth.value ** 2 + 1716.71568 * model.reserv.depth.value + 500866.89110) * 1E-6
elif self.wellcorrelation.value == WellDrillingCostCorrelation.VERTICAL_SMALL_INT1:
self.C1well = (0.13710 * model.reserv.depth.value ** 2 + 129.61033 * model.reserv.depth.value + 1205587.57100) * 1E-6
elif self.wellcorrelation.value == WellDrillingCostCorrelation.VERTICAL_LARGE_INT1:
self.C1well = (0.18927 * model.reserv.depth.value ** 2 + 293.45174 * model.reserv.depth.value + 1326526.31300) * 1E-6
elif self.wellcorrelation.value == WellDrillingCostCorrelation.DEVIATED_SMALL_INT1:
self.C1well = (0.15340 * model.reserv.depth.value ** 2 + 120.31700 * model.reserv.depth.value + 1431801.54400) * 1E-6
elif self.wellcorrelation.value == WellDrillingCostCorrelation.DEVIATED_LARGE_INT1:
self.C1well = (0.19950 * model.reserv.depth.value ** 2 + 296.13011 * model.reserv.depth.value + 1697867.70900) * 1E-6
elif self.wellcorrelation.value == WellDrillingCostCorrelation.VERTICAL_SMALL_INT2:
self.C1well = (0.00804 * model.reserv.depth.value ** 2 + 455.60507 * model.reserv.depth.value + 921007.68680) * 1E-6
elif self.wellcorrelation.value == WellDrillingCostCorrelation.VERTICAL_LARGE_INT2:
self.C1well = (0.00315 * model.reserv.depth.value ** 2 + 782.69676 * model.reserv.depth.value + 983620.25270) * 1E-6
elif self.wellcorrelation.value == WellDrillingCostCorrelation.DEVIATED_SMALL_INT2:
self.C1well = (0.00854 * model.reserv.depth.value ** 2 + 506.08357 * model.reserv.depth.value + 1057330.39000) * 1E-6
elif self.wellcorrelation.value == WellDrillingCostCorrelation.DEVIATED_LARGE_INT2:
self.C1well = (0.00380 * model.reserv.depth.value ** 2 + 838.90249 * model.reserv.depth.value + 1181947.04400) * 1E-6
elif self.wellcorrelation.value == WellDrillingCostCorrelation.VERTICAL_SMALL_IDEAL:
self.C1well = (0.00252 * model.reserv.depth.value ** 2 + 439.44503 * model.reserv.depth.value + 590611.90110) * 1E-6
elif self.wellcorrelation.value == WellDrillingCostCorrelation.VERTICAL_LARGE_IDEAL:
self.C1well = (-0.00240 * model.reserv.depth.value ** 2 + 752.93946 * model.reserv.depth.value + 524337.65380) * 1E-6
elif self.wellcorrelation.value == WellDrillingCostCorrelation.DEVIATED_SMALL_IDEAL:
self.C1well = (0.00719 * model.reserv.depth.value ** 2 + 455.85233 * model.reserv.depth.value + 753377.73080) * 1E-6
elif self.wellcorrelation.value == WellDrillingCostCorrelation.DEVIATED_LARGE_IDEAL:
self.C1well = (0.00376 * model.reserv.depth.value ** 2 + 762.52696 * model.reserv.depth.value + 765103.07690) * 1E-6

if self.wellcorrelation.value == WellDrillingCostCorrelation.SIMPLE:
self.C1well = self.wellcorrelation.value.calculate_cost_MUSD(model.reserv.depth.value)

self.C1well = self.C1well * self.ccwelladjfactor.value
self.Cwell.value = self.C1well * (model.wellbores.nprod.value + model.wellbores.ninj.value)
Expand Down

0 comments on commit 25553d4

Please sign in to comment.