From c4df93b127757c288c7f9f8b18a3732d1ebdc082 Mon Sep 17 00:00:00 2001
From: aheinri5 <Anna@netzkritzler.de>
Date: Tue, 23 Feb 2021 20:36:33 +0100
Subject: [PATCH] [docs] added plot with logarithmic storage
---
vignettes/ensemble_workflow_vignette.Rmd | 70 +++++++++++++++++++++++-
1 file changed, 68 insertions(+), 2 deletions(-)
diff --git a/vignettes/ensemble_workflow_vignette.Rmd b/vignettes/ensemble_workflow_vignette.Rmd
index 235b2b3..f26dfcb 100644
--- a/vignettes/ensemble_workflow_vignette.Rmd
+++ b/vignettes/ensemble_workflow_vignette.Rmd
@@ -213,11 +213,77 @@ ggplot(per_por_sto_df,
geom_point(aes(color = as.factor(perc))) +
xlab("x point coordinate") +
labs(color = "%") +
- facet_grid(cols = vars(name))
+ facet_grid(cols = vars(name),
+ labeller = as_labeller(c(per = "permeability",
+ por = "porosity",
+ sto = "storage")))
```
-Ta-Daa! We can see `permeability` has the most influence on the pressure. Though they may seem suspicious, `porosity` and `storage` are being plotted correctly - the points are just being placed right on top of each other. Since `porosity` can't go over the value `1` which was the original value, our value vector only went from -50% to 0% which is why the line colors of `porosity` and `storage` differ.
+Ta-Daa! We can see `permeability` has the most influence on the pressure. Though they may seem suspicious, `porosity` and `storage` are being plotted correctly - the points are just being placed right on top of each other. Since `porosity` can't go over the value `1` which was the original value, our value vector only went from -50% to 0% which is why the line colors of `porosity` and `storage` differ. Maybe we want to try and use a logarithmic approach for `storage`. This won't work with the built-in functionality of `OGS6_Ensemble` so we'll set up our Ensemble a little differently.
+
+```{r}
+# Calculate log value
+log_val <- log(as.numeric(
+ ogs6_obj$media[[1]]$properties[[4]]$value),
+ base = 10)
+
+# Apply changes to log value
+log_vals <- vapply(percentages, function(x){
+ log_val + (log_val * (x / 100))
+}, FUN.VALUE = numeric(1))
+
+# Transfer back to non-logarithmic scale
+back_transf_vals <- 10^log_vals
+
+# Change sim_path to fit your system
+ogs6_obj$sim_path <- "D:/ogs6_theis/axisym_theis_sim_log_storage"
+# Set up new ensemble
+ogs6_ens_sto <-
+ OGS6_Ensemble$new(
+ ogs6_obj = ogs6_obj,
+ parameters =
+ list(
+ sto = list(
+ ogs6_obj$media[[1]]$properties[[4]]$value,
+ values = back_transf_vals)
+ ),
+ percentages_mode = FALSE,
+ sequential_mode = TRUE
+ )
+```
+
+As before, we can run the simulation right away.
+
+```{r, eval = FALSE}
+ogs6_ens_sto$run_simulation()
+```
+
+```{r, include = FALSE}
+lapply(ogs6_ens_sto$ensemble, ogs6_read_output_files)
+```
+
+Let's check if we can observe any influence of `storage` on `pressure` now.
+
+```{r}
+# Get combined dataframe
+sto_df <-
+ ogs6_ens_sto$get_point_data(
+ keys = c("pressure"),
+ start_at_timestep = ogs6_ens_sto$ensemble[[1]]$pvds[[1]]$last_timestep)
+
+# Supply percentages manually since we couldn't use `percentages_mode`
+percs <- vapply(sto_df$sim_id,
+ function(x){percentages[[x]]},
+ FUN.VALUE = numeric(1))
+
+ggplot(sto_df,
+ aes(x = x,
+ y = pressure)) +
+ geom_point(aes(color = as.factor(percs))) +
+ xlab("x point coordinate") +
+ labs(color = "%")
+```
## Theis solution for well pumping
--
GitLab