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.

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

ENH: Add calculation of water zone volumes #925

Merged
merged 1 commit into from
Dec 19, 2024
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
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
42 changes: 39 additions & 3 deletions src/fmu/dataio/export/rms/inplace_volumes.py
Original file line number Diff line number Diff line change
Expand Up @@ -124,6 +124,43 @@ def _convert_table_from_rms_to_legacy_format(table: pd.DataFrame) -> pd.DataFram
columns="REAL", errors="ignore"
)

@staticmethod
def _compute_water_zone_volumes_from_totals(table: pd.DataFrame) -> pd.DataFrame:
"""
Calculate 'water' zone volumes by subtracting HC-zone volumes from 'Total'
volumes which represents the entire zone. Total volumes are removed after
'water' zone volumes have been added to the table.
"""
_logger.debug("Computing water volumes from Totals...")

total_suffix = "_TOTAL"
total_columns = [col for col in table.columns if col.endswith(total_suffix)]

if not total_columns:
raise RuntimeError(
"Found no 'Totals' volumes in the table. Please ensure 'Totals' "
"are reported and rerun the volumetric job before export."
)
Comment on lines +139 to +143
Copy link
Collaborator

Choose a reason for hiding this comment

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

This can also be captured more upfront by evaluation the job settings which already are read

(this example from Drogon)

myjob = rmsjobs.Job.get_job(
    owner=["Grid models", "Geogrid", "Grid"],
    type="Volumetrics",
    name="geogrid_volumes",
)

print("USE TOTALS", mydict["Output"][0]["UseTotals"])

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 know, did initially check several things on the job settings e.g. required calculations.. but found it difficult to write tests, had to mock so much. I'll keep it here for now and we can move later if we want 👍


for total_col in total_columns:
volumetric_col = total_col.replace(total_suffix, "")

water_zone_col = f"{volumetric_col}_WATER"
oil_zone_col = f"{volumetric_col}_OIL"
gas_zone_col = f"{volumetric_col}_GAS"

# first set water zone data equal to the Total
# then subtract data from the oil/gas zone
table[water_zone_col] = table[total_col]

if oil_zone_col in table:
table[water_zone_col] -= table[oil_zone_col]

if gas_zone_col in table:
table[water_zone_col] -= table[gas_zone_col]

return table.drop(columns=total_columns)

@staticmethod
def _add_missing_columns_to_table(table: pd.DataFrame) -> pd.DataFrame:
"""Add columns with nan values if not present in table."""
Expand All @@ -149,14 +186,12 @@ def _transform_and_add_fluid_column_to_table(
are renamed into 'BULK' and 'PORV' columns. To separate the data an additional
FLUID column is added that indicates the type of fluid the row represents.
"""
table_index = [
col for col in _enums.InplaceVolumes.index_columns() if col in table
]
Comment on lines -152 to -154
Copy link
Collaborator Author

Choose a reason for hiding this comment

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

This was also given as input argument


tables = []
for fluid in (
_enums.InplaceVolumes.Fluid.gas.value,
_enums.InplaceVolumes.Fluid.oil.value,
_enums.InplaceVolumes.Fluid.water.value,
):
fluid_suffix = fluid.upper()
fluid_columns = [
Expand Down Expand Up @@ -188,6 +223,7 @@ def _convert_table_from_legacy_to_standard_format(
table_index = [
col for col in _enums.InplaceVolumes.index_columns() if col in table
]
table = self._compute_water_zone_volumes_from_totals(table)
table = self._transform_and_add_fluid_column_to_table(table, table_index)
table = self._add_missing_columns_to_table(table)
return self._set_table_column_order(table)
Expand Down
168 changes: 168 additions & 0 deletions tests/data/drogon/tabular/volumes/geogrid.csv
Original file line number Diff line number Diff line change
Expand Up @@ -335,3 +335,171 @@ oil,Volon,EastLowland,Calcite,,86379.79939277154,,0.0,0.0,0.0,,0.0,
oil,Volon,EastLowland,Offshore,,0.0,,0.0,0.0,0.0,,0.0,
oil,Volon,EastLowland,Lowershoreface,,0.0,,0.0,0.0,0.0,,0.0,
oil,Volon,EastLowland,Uppershoreface,,0.0,,0.0,0.0,0.0,,0.0,
water,Valysar,WestLowland,Floodplain,,199635209.06231964,,19889504.858085632,,,,,
water,Valysar,WestLowland,Channel,,188904690.15197682,,49852644.8797493,,,,,
water,Valysar,WestLowland,Crevasse,,66498186.74765427,,13614114.21556282,,,,,
water,Valysar,WestLowland,Coal,,47005819.028200254,,0.0,,,,,
water,Valysar,WestLowland,Calcite,,0.0,,0.0,,,,,
water,Valysar,WestLowland,Offshore,,0.0,,0.0,,,,,
water,Valysar,WestLowland,Lowershoreface,,0.0,,0.0,,,,,
water,Valysar,WestLowland,Uppershoreface,,0.0,,0.0,,,,,
water,Valysar,CentralSouth,Floodplain,,28174636.855486907,,2711703.3576335907,,,,,
water,Valysar,CentralSouth,Channel,,22735443.45432435,,6156691.954032898,,,,,
water,Valysar,CentralSouth,Crevasse,,9887763.792533554,,2015381.9916992192,,,,,
water,Valysar,CentralSouth,Coal,,4107955.788754794,,0.0,,,,,
water,Valysar,CentralSouth,Calcite,,0.0,,0.0,,,,,
water,Valysar,CentralSouth,Offshore,,0.0,,0.0,,,,,
water,Valysar,CentralSouth,Lowershoreface,,0.0,,0.0,,,,,
water,Valysar,CentralSouth,Uppershoreface,,0.0,,0.0,,,,,
water,Valysar,CentralNorth,Floodplain,,38746747.15721468,,3257848.738916397,,,,,
water,Valysar,CentralNorth,Channel,,20313542.091495283,,5319292.084007263,,,,,
water,Valysar,CentralNorth,Crevasse,,12964778.222036486,,2584124.102832794,,,,,
water,Valysar,CentralNorth,Coal,,14076425.203160832,,0.0,,,,,
water,Valysar,CentralNorth,Calcite,,0.0,,0.0,,,,,
water,Valysar,CentralNorth,Offshore,,0.0,,0.0,,,,,
water,Valysar,CentralNorth,Lowershoreface,,0.0,,0.0,,,,,
water,Valysar,CentralNorth,Uppershoreface,,0.0,,0.0,,,,,
water,Valysar,NorthHorst,Floodplain,,8210484.753254818,,770498.3396244047,,,,,
water,Valysar,NorthHorst,Channel,,2578974.2676892485,,687172.2579650879,,,,,
water,Valysar,NorthHorst,Crevasse,,761616.8469049805,,157481.72686386108,,,,,
water,Valysar,NorthHorst,Coal,,2037701.0165579505,,0.0,,,,,
water,Valysar,NorthHorst,Calcite,,0.0,,0.0,,,,,
water,Valysar,NorthHorst,Offshore,,0.0,,0.0,,,,,
water,Valysar,NorthHorst,Lowershoreface,,0.0,,0.0,,,,,
water,Valysar,NorthHorst,Uppershoreface,,0.0,,0.0,,,,,
water,Valysar,CentralRamp,Floodplain,,1061413.6835421966,,82375.53680419922,,,,,
water,Valysar,CentralRamp,Channel,,689661.5301341028,,178625.55529785156,,,,,
water,Valysar,CentralRamp,Crevasse,,68555.0179801248,,13122.741760253906,,,,,
water,Valysar,CentralRamp,Coal,,133394.72849460552,,0.0,,,,,
water,Valysar,CentralRamp,Calcite,,0.0,,0.0,,,,,
water,Valysar,CentralRamp,Offshore,,0.0,,0.0,,,,,
water,Valysar,CentralRamp,Lowershoreface,,0.0,,0.0,,,,,
water,Valysar,CentralRamp,Uppershoreface,,0.0,,0.0,,,,,
water,Valysar,CentralHorst,Floodplain,,8923441.539376646,,737789.9258155823,,,,,
water,Valysar,CentralHorst,Channel,,3874852.635276627,,1028251.5404396057,,,,,
water,Valysar,CentralHorst,Crevasse,,1046442.6939407401,,212083.40999984718,,,,,
water,Valysar,CentralHorst,Coal,,4122777.8013486303,,0.0,,,,,
water,Valysar,CentralHorst,Calcite,,0.0,,0.0,,,,,
water,Valysar,CentralHorst,Offshore,,0.0,,0.0,,,,,
water,Valysar,CentralHorst,Lowershoreface,,0.0,,0.0,,,,,
water,Valysar,CentralHorst,Uppershoreface,,0.0,,0.0,,,,,
water,Valysar,EastLowland,Floodplain,,138781101.69286424,,10602248.395273207,,,,,
water,Valysar,EastLowland,Channel,,57556967.579790935,,14855282.31137085,,,,,
water,Valysar,EastLowland,Crevasse,,32904003.154506125,,6366257.547336578,,,,,
water,Valysar,EastLowland,Coal,,33782575.87539426,,0.0,,,,,
water,Valysar,EastLowland,Calcite,,0.0,,0.0,,,,,
water,Valysar,EastLowland,Offshore,,0.0,,0.0,,,,,
water,Valysar,EastLowland,Lowershoreface,,0.0,,0.0,,,,,
water,Valysar,EastLowland,Uppershoreface,,0.0,,0.0,,,,,
water,Therys,WestLowland,Floodplain,,0.0,,0.0,,,,,
water,Therys,WestLowland,Channel,,0.0,,0.0,,,,,
water,Therys,WestLowland,Crevasse,,0.0,,0.0,,,,,
water,Therys,WestLowland,Coal,,0.0,,0.0,,,,,
water,Therys,WestLowland,Calcite,,53179530.480960734,,0.0,,,,,
water,Therys,WestLowland,Offshore,,224308300.77785882,,27138859.65850281,,,,,
water,Therys,WestLowland,Lowershoreface,,157251330.06353703,,40617617.03001963,,,,,
water,Therys,WestLowland,Uppershoreface,,111674468.20940188,,34788357.914711,,,,,
water,Therys,CentralSouth,Floodplain,,0.0,,0.0,,,,,
water,Therys,CentralSouth,Channel,,0.0,,0.0,,,,,
water,Therys,CentralSouth,Crevasse,,0.0,,0.0,,,,,
water,Therys,CentralSouth,Coal,,0.0,,0.0,,,,,
water,Therys,CentralSouth,Calcite,,1570953.4599857288,,0.0,,,,,
water,Therys,CentralSouth,Offshore,,5148527.537121184,,480000.0074880719,,,,,
water,Therys,CentralSouth,Lowershoreface,,5091024.065522364,,1325732.81913805,,,,,
water,Therys,CentralSouth,Uppershoreface,,3061684.296614888,,994067.8604154126,,,,,
water,Therys,CentralNorth,Floodplain,,0.0,,0.0,,,,,
water,Therys,CentralNorth,Channel,,0.0,,0.0,,,,,
water,Therys,CentralNorth,Crevasse,,0.0,,0.0,,,,,
water,Therys,CentralNorth,Coal,,0.0,,0.0,,,,,
water,Therys,CentralNorth,Calcite,,8246046.31052004,,0.0,,,,,
water,Therys,CentralNorth,Offshore,,30136718.68757325,,2900918.633241721,,,,,
water,Therys,CentralNorth,Lowershoreface,,20880188.429421045,,5023761.500228882,,,,,
water,Therys,CentralNorth,Uppershoreface,,12440213.544097727,,4028126.994163514,,,,,
water,Therys,NorthHorst,Floodplain,,0.0,,0.0,,,,,
water,Therys,NorthHorst,Channel,,0.0,,0.0,,,,,
water,Therys,NorthHorst,Crevasse,,0.0,,0.0,,,,,
water,Therys,NorthHorst,Coal,,0.0,,0.0,,,,,
water,Therys,NorthHorst,Calcite,,3910879.9040033277,,0.0,,,,,
water,Therys,NorthHorst,Offshore,,8525384.040342128,,898389.3357550651,,,,,
water,Therys,NorthHorst,Lowershoreface,,1813954.4037235645,,448523.21085357666,,,,,
water,Therys,NorthHorst,Uppershoreface,,2230344.651985049,,726788.3869018554,,,,,
water,Therys,CentralRamp,Floodplain,,0.0,,0.0,,,,,
water,Therys,CentralRamp,Channel,,0.0,,0.0,,,,,
water,Therys,CentralRamp,Crevasse,,0.0,,0.0,,,,,
water,Therys,CentralRamp,Coal,,0.0,,0.0,,,,,
water,Therys,CentralRamp,Calcite,,14269.20473185298,,0.0,,,,,
water,Therys,CentralRamp,Offshore,,570056.6082884455,,42336.00001788139,,,,,
water,Therys,CentralRamp,Lowershoreface,,573928.800745388,,133559.80555534363,,,,,
water,Therys,CentralRamp,Uppershoreface,,283942.140126402,,93575.89129638672,,,,,
water,Therys,CentralHorst,Floodplain,,0.0,,0.0,,,,,
water,Therys,CentralHorst,Channel,,0.0,,0.0,,,,,
water,Therys,CentralHorst,Crevasse,,0.0,,0.0,,,,,
water,Therys,CentralHorst,Coal,,0.0,,0.0,,,,,
water,Therys,CentralHorst,Calcite,,1873912.9772147671,,0.0,,,,,
water,Therys,CentralHorst,Offshore,,4164290.322972057,,413992.2520942539,,,,,
water,Therys,CentralHorst,Lowershoreface,,7019327.266154291,,1655530.3278546333,,,,,
water,Therys,CentralHorst,Uppershoreface,,4653060.394801656,,1524199.657714963,,,,,
water,Therys,EastLowland,Floodplain,,0.0,,0.0,,,,,
water,Therys,EastLowland,Channel,,0.0,,0.0,,,,,
water,Therys,EastLowland,Crevasse,,0.0,,0.0,,,,,
water,Therys,EastLowland,Coal,,0.0,,0.0,,,,,
water,Therys,EastLowland,Calcite,,7565198.396226491,,0.0,,,,,
water,Therys,EastLowland,Offshore,,32971054.08148795,,2403078.7679725736,,,,,
water,Therys,EastLowland,Lowershoreface,,42226056.55160797,,9751004.178666055,,,,,
water,Therys,EastLowland,Uppershoreface,,30247661.497023392,,9434800.083554983,,,,,
water,Volon,WestLowland,Floodplain,,151553410.404796,,13848884.900240175,,,,,
water,Volon,WestLowland,Channel,,314039375.983029,,76394654.75561242,,,,,
water,Volon,WestLowland,Crevasse,,0.0,,0.0,,,,,
water,Volon,WestLowland,Coal,,0.0,,0.0,,,,,
water,Volon,WestLowland,Calcite,,42194832.373705655,,0.0,,,,,
water,Volon,WestLowland,Offshore,,0.0,,0.0,,,,,
water,Volon,WestLowland,Lowershoreface,,0.0,,0.0,,,,,
water,Volon,WestLowland,Uppershoreface,,0.0,,0.0,,,,,
water,Volon,CentralSouth,Floodplain,,50301143.083524644,,4752792.43396759,,,,,
water,Volon,CentralSouth,Channel,,111281241.46517287,,27250731.639530182,,,,,
water,Volon,CentralSouth,Crevasse,,0.0,,0.0,,,,,
water,Volon,CentralSouth,Coal,,0.0,,0.0,,,,,
water,Volon,CentralSouth,Calcite,,13374219.834380861,,0.0,,,,,
water,Volon,CentralSouth,Offshore,,0.0,,0.0,,,,,
water,Volon,CentralSouth,Lowershoreface,,0.0,,0.0,,,,,
water,Volon,CentralSouth,Uppershoreface,,0.0,,0.0,,,,,
water,Volon,CentralNorth,Floodplain,,13437653.734824093,,1219808.1165828705,,,,,
water,Volon,CentralNorth,Channel,,18540833.941328164,,4579542.983924866,,,,,
water,Volon,CentralNorth,Crevasse,,0.0,,0.0,,,,,
water,Volon,CentralNorth,Coal,,0.0,,0.0,,,,,
water,Volon,CentralNorth,Calcite,,2912524.102400574,,0.0,,,,,
water,Volon,CentralNorth,Offshore,,0.0,,0.0,,,,,
water,Volon,CentralNorth,Lowershoreface,,0.0,,0.0,,,,,
water,Volon,CentralNorth,Uppershoreface,,0.0,,0.0,,,,,
water,Volon,NorthHorst,Floodplain,,2485998.1940900474,,220445.77733993536,,,,,
water,Volon,NorthHorst,Channel,,5144868.677435892,,1311122.0843620305,,,,,
water,Volon,NorthHorst,Crevasse,,0.0,,0.0,,,,,
water,Volon,NorthHorst,Coal,,0.0,,0.0,,,,,
water,Volon,NorthHorst,Calcite,,755851.579945738,,0.0,,,,,
water,Volon,NorthHorst,Offshore,,0.0,,0.0,,,,,
water,Volon,NorthHorst,Lowershoreface,,0.0,,0.0,,,,,
water,Volon,NorthHorst,Uppershoreface,,0.0,,0.0,,,,,
water,Volon,CentralRamp,Floodplain,,1002797.0335718896,,106327.63387680054,,,,,
water,Volon,CentralRamp,Channel,,2912622.8234018395,,755271.9564285278,,,,,
water,Volon,CentralRamp,Crevasse,,0.0,,0.0,,,,,
water,Volon,CentralRamp,Coal,,0.0,,0.0,,,,,
water,Volon,CentralRamp,Calcite,,183954.78222594503,,0.0,,,,,
water,Volon,CentralRamp,Offshore,,0.0,,0.0,,,,,
water,Volon,CentralRamp,Lowershoreface,,0.0,,0.0,,,,,
water,Volon,CentralRamp,Uppershoreface,,0.0,,0.0,,,,,
water,Volon,CentralHorst,Floodplain,,14479252.611059882,,1352125.556652069,,,,,
water,Volon,CentralHorst,Channel,,21087337.78153704,,5291756.395442961,,,,,
water,Volon,CentralHorst,Crevasse,,0.0,,0.0,,,,,
water,Volon,CentralHorst,Coal,,0.0,,0.0,,,,,
water,Volon,CentralHorst,Calcite,,2670829.178936201,,0.0,,,,,
water,Volon,CentralHorst,Offshore,,0.0,,0.0,,,,,
water,Volon,CentralHorst,Lowershoreface,,0.0,,0.0,,,,,
water,Volon,CentralHorst,Uppershoreface,,0.0,,0.0,,,,,
water,Volon,EastLowland,Floodplain,,44700188.399142474,,3830575.4531288147,,,,,
water,Volon,EastLowland,Channel,,112391110.29544623,,26596352.009635925,,,,,
water,Volon,EastLowland,Crevasse,,0.0,,0.0,,,,,
water,Volon,EastLowland,Coal,,0.0,,0.0,,,,,
water,Volon,EastLowland,Calcite,,14198059.664785597,,0.0,,,,,
water,Volon,EastLowland,Offshore,,0.0,,0.0,,,,,
water,Volon,EastLowland,Lowershoreface,,0.0,,0.0,,,,,
water,Volon,EastLowland,Uppershoreface,,0.0,,0.0,,,,,
Loading
Loading