Merge pull request #7 from slimgroup/abhinav

Abhinav

GahlotLi2024SEG/paper.bib

@@ -0,0 +1,48 @@
+@conference {gahlot2023NIPSWSifp,
+	title = {Inference of CO2 flow patterns {\textendash} a feasibility study},
+	booktitle = {Neural Information Processing Systems (NeurIPS)},
+	year = {2023},
+	note = {(NeurIPS 2023 Workshop - Tackling Climate Change with Machine Learning (Spotlight))},
+	month = {10},
+	abstract = {},
+	keywords = {Amortized Variational Inference, Bayesian inference, CCS, conditional normalizing flows, geological carbon storage, monitoring, NIPS, Summary Statistics},
+	doi = {10.48550/arXiv.2311.00290},
+	url = {https://slim.gatech.edu/Publications/Public/Conferences/NIPS/2023/gahlot2023NIPSWSifp/paper.pdf},
+	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/NIPS/2023/gahlot2023NIPSWSifp/poster.pdf},
+	author = {Abhinav Prakash Gahlot and Huseyin Tuna Erdinc and Rafael Orozco and Ziyi Yin and Felix J. Herrmann}
+}
+
+@article{herrmann2023president,
+  title={President's Page: Digital twins in the era of generative AI},
+  author={Herrmann, Felix J},
+  journal={The Leading Edge},
+  volume={42},
+  number={11},
+  pages={730--732},
+  year={2023},
+  publisher={Society of Exploration Geophysicists}
+}
+
+@article{thelen2022comprehensivea,
+  title={A comprehensive review of digital twin—part 1: modeling and twinning enabling technologies},
+  author={Thelen, Adam and Zhang, Xiaoge and Fink, Olga and Lu, Yan and Ghosh, Sayan and Youn, Byeng D and Todd, Michael D and Mahadevan, Sankaran and Hu, Chao and Hu, Zhen},
+  journal={Structural and Multidisciplinary Optimization},
+  volume={65},
+  number={12},
+  pages={354},
+  year={2022},
+  publisher={Springer}
+}
+
+@conference {yu2023IMAGEmsc,
+	title = {Monitoring Subsurface CO2 Plumes with Sequential Bayesian Inference},
+	booktitle = {International Meeting for Applied Geoscience and Energy},
+	year = {2023},
+	note = {(IMAGE, Houston)},
+	month = {08},
+	abstract = {To monitor and predict CO2 plume dynamics during geological carbon storage, reservoir engineers usually perform two-phase flow simulations. While these simulations may provide useful insights, their usefulness is limited due to numerous complicating factors including uncertainty in the dynamics of the plume itself. To study this phenomenon, we consider stochasticity in the dynamic caused by unknown random changes in the injection rate. By conditioning the CO2 plume predictions on seismic observations, we correct the CO2 plume predictions and quantify uncertainty with machine learning.},
+	keywords = {Bayesian inference, CCS, deep learning, Imaging, monitoring, SEG, Uncertainty quantification},
+	url = {https://slimgroup.github.io/IMAGE2023/SequentialBayes/abstract.html},
+	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2023/yu2023IMAGEmsc},
+	author = {Ting-ying Yu and Rafael Orozco and Ziyi Yin and Mathias Louboutin and Felix J. Herrmann}
+}

GahlotLi2024SEG/paper.qmd

@@ -0,0 +1,44 @@
+---
+title: "A Digital Twin for Geological Carbon Storage with Controlled Injectivity"
+author:
+  - name:
+      Abhinav Prakash Gahlot^1\*^, Haoyun Li,^1\*^, Ziyi Yin^1^, Rafael Orozco^1^, and Felix J. Herrmann^1^ \
+
+      ^1^ Georgia Institute of Technology, ^\*^ first two authors contributed equally
+bibliography: paper.bib
+---
+
+## Introduction
+
+Digital Twins refer to dynamic virtual replicas of subsurface systems, integrating real-time data and employing advanced generative Artificial Intelligence (genAI) methodologies, such as neural posterior density estimation via simulation-based inference and sequential Bayesian inference. Thanks to combination of these advanced Bayesian inference techniques, our approach is capable of addressing challenges of monitoring and controlling CO~2~ storage projects. These challenges include dealing with the subsurface's complexity and heterogeneity (seismic and fluid-flow properties), operations optimization, and risk mitigation, e.g. via injection rate control. Because our Digital Twin is capable of handling diverse monitoring data, consisting of time-lapse seismic and data collected at (monitoring) wells, it entails a technology that serves as a platform to integrate seemingly disparate and siloed fields, e.g. geophysics and reservoir engineering. In addition, recent breakthroughs in genAI, allow Digital Twins to capture uncertainty in a principled way [@yu2023IMAGEmsc;@herrmann2023president;@gahlot2023NIPSWSifp]. By employing training and inference recursively, the Digital Twin trains its neural networks on samples of the simulated current state---i.e., the CO~2~ saturation/pressure, paired with simulated imaged seismic and/or data collected at (monitoring) wells. These training pairs of the simulated state and simulated observations are obtained by sampling the posterior distribution, $\mathbf{x}_{k-1}\sim p(\mathbf{x}_{k-1}\vert \mathbf{y}^\mathrm{o}_{1:k-1})$, at the previous timestep, $k-1$, conditioned on field data, $\mathbf{y}^\mathrm{o}_{1:k-1}$, collected over all previous timesteps, $1:k-1$, followed by advancing the state to the current timestep, followed by simulating (seismic/well) observations associated with that state. Given these simulated state-observation pairs, the Digital Twin's networks are trained, so they are current and ready to produce samples of the posterior when the new field data comes in---i.e. $\mathbf{x}_{k}\sim p(\mathbf{x}_{k}\vert \mathbf{y}^\mathrm{o}_{1:k})$. While this new neural approach to data assimilation for CO~2~ storage projects provides what is called an uncertainty-informed *Digital Shadow*, it lacks decision making and control, which would make it a Digital Twin [@thelen2022comprehensivea], capable of optimizing storage operations while mitigating risks including the risk of fracturing the cap rock by exceeding the fracture pressure. The latter risk is illustrated in Figure 1, where the first row contains simulated samples of the differential pressure at timestep, $k=4$, without control. These samples for the simulated state exceed the fracture pressure and are denoted by the red areas. During this talk, we will demonstrate how the Digital Twin can make informed decisions to avoid exceeding the fracture pressure.
+
+
+## Methodology
+
+To make uncertainty informed decisions on adapting the injection rate, 128 samples of the state (see second row Figure 1), $\{\mathbf{x}^{(m)}_{3}\}_{m=1}^{128}$, and permeability, $\{K^{(m)}\}_{m=1}^{128}$, are drawn at time-step $k=3$ from the posterior distribution, $\mathbf{x}_{3}\sim p(\mathbf{x}_{3}\vert \mathbf{y}^\mathrm{o}_{1:3})$, for the state conditioned on the observed time-lapse data, and from the distribution for the permeability, $K\sim p(K)$. To find the optimized injection rate, we first calculate for each sample, $K$ and $\mathbf{x}^{(m)}_{3}$, the optimized injectivity by $\max_{q_3} q_3\Delta t \ \ \text{subject to} \ \ \mathbf{x}_{4}['p']<\mathbf{p}_{\max}$ where $\mathbf{p}_{\max}$ is the depth-dependent fracture pressure, and $\mathbf{x}_{4}=\mathcal{M}_3(\mathbf{x}_{3}, \mathbf{K}; q_3)$ the state's time-advancement denoted by the symbol $\mathcal{M}_3$. For the fluid-flow simulations, the open-source Julia package [JUDI.jl](https://github.com/slimgroup/JUDI.jl) and [JutulDarcy.jl](https://github.com/sintefmath/JutulDarcy.jl) are used. Results of these optimizations are included in Figure 2(a), which contains a histogram of the empirical fracture frequency as a function of injection rates at time $k=3$. Given this histogram, our task is to maximize the injection rate given a pre-defined confidence interval (e.g. 95 %), so that the fracture probability remains below a certain percentage e.g. 1%.  With these simulations, and the fact that fracture/no-fracture occurrences are distributed according to the Bernoulli distribution, we will demonstrate that we are able to select an injection rate that limits fracture occurrence to the prescribed probability with a prescribed confidence interval. For instance, we can compute $Pr([\mathbf{x}_4['p']>p_{\max}]<0.01)<1-0.025$, which corresponds to selecting an injection rate that leads to fracture rate of $<1\%$ with $97.5\%$ confidence. The confidence interval is halved, because only conservative (left) injection rates will be selected (see Figure 2(b)).
+
+
+## Results
+
+To calculate injection rates that mitigate the risk of exceeding the fracture pressure, we proceed as follows. First, because the optimized injection rates are close to the fracture pressure, we consider these optimization as approximations to the injection rates where the fracture pressure are exceeded. Next, Kernel Density Estimation (KDE) is applied to produce the smooth red curve in Figure 2(a). This smoothed probability function is used to calculate the Cumulative Density Function (CDF), plotted in Figure 2(b). Using the fact that non-fracture/fracture occurrence entails a Bernoulli distribution, confidence intervals can be calculated, $\pm Z_{\frac{\alpha}{2}} \sqrt{\frac{\hat{p}(1 - \hat{p})}{128}}$ where $\hat{p}$ represents the CDF (blue line) and $Z_{\frac{\alpha}{2}}=1.96$ with $\alpha=0.05$. From the CDF and confidence intervals (denoted by the grey areas), the following conclusions can be drawn: First, if the initial injection rate of $q_3 = 0.0500 m^3/s$ is kept, the fracture probability lies between 24.47 -- 40.71% (vertical dashed line) and has a maximum likelihood of 32.59%, which are all way too high. Second, if we want to limit the fracture occurrence rate to 1% (red dashed line), then we need to lower the injection rate to $q_3=0.0387\mathrm{m^3/s}$. To ensure the low fracture occurrence rate of 1%, the reduced injection rate is chosen as the smallest injection rate within the confidence interval. As can be observed from Figure 2(b), lowering the injection rate avoids exceeding the fracture pressure at the expense of injecting less CO~2~. Out of 128 samples, 43 samples are fractured with the initial injection rate, while only one sample is fractured with the controlled injection rate. 
+
+## Conclusion and discussions
+
+The above example illustrates how Digital Twins can be used to mitigate risks associated with CO~2~ storage projects. Specifically, we used the Digital Twin's capability to produce samples of its state (pressure), conditioned on observed seismic and/or well data. Using these samples, in conjunction with samples from the permeability distribution, we were able to capture statistics on the fracture occurrence frequency as a function of the injection rate. Given these statistics, we were in a position to set a fracture frequency and choose the corresponding injection rate as a function of the confidence interval. By following this procedure, exceeding the fracture pressure was avoided by lowering the injection rate. The decision to lower the injection rate, and by which amount, was informed by the Digital Twin, which uses seismic and/or well data to capture reservoir's state including its uncertainty.
+
+::: {#fig-flow}
+![](./figs/Figure1.png){width="90%"}
+
+Pressure state at time-step $k=3$ and a comparative analysis of pressure outputs from the digital twin at time-step $k=4$
+:::
+
+::: {#fig-fracture}
+![](./figs/Figure2.png){width="90%"}
+
+Injection rate samples and KDE curve(a), CDF curve of fracture probability versus injection rate with 95% confidence interval(b)
+:::
+
+## References
+
+::: {#refs}
+:::

_quarto.yml

@@ -16,8 +16,8 @@ website:
       contents:
       - section: "Abstracts"
         contents:
-        - file: DigitalTwin/abstract.qmd
-          text: "Digital twin with control"
+        - file: GahlotLi2024SEG/abstract.qmd
+          text: "A Digital Twin for Geological Carbon Storage with Controlled Injectivity"
         - file: yin2024SEG/paper.qmd
           text: "WISER: full-Waveform variational Inference via Subsurface Extensions with Refinements"
   page-footer:

index.qmd

@@ -6,8 +6,8 @@ This is a Quarto website.
 
 All submissions to the Image24 conference with additional figures, references, ...
 
-List of abstrac:
+List of abstracts:
 
-- [Digital twin with control](DigitalTwin/abstract.qmd) An uncertainty-aware digital twin for geological carbon storage
+- [Digital twin with control](GahlotLi2024SEG/abstract.qmd) A Digital Twin for Geological Carbon Storage with Controlled Injectivity
 
 - [WISER](yin2024SEG/paper.qmd): full-Waveform variational Inference via Subsurface Extensions with Refinements