the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Numerical Analysis of the Effect of Heterogeneity on CO2 Dissolution Enhanced by Gravity Driven Convection
Abstract. Dissolution trapping of CO2 in brine can mitigate the risk of supercritical CO2 leakage during long-term Geological Carbon Sequestration (GCS). The dissolution of overlying supercritical CO2 into brine increases the density of brine in its upper portion, which causes Gravity-Driven Convection (GDC) and thus significantly increases the rate of CO2 dissolution. To date, most studies on GDC-driven dissolution are based on homogeneous media and only few studies exist on the effect of heterogeneity on GDC-driven dissolution. Here, we study the effect of heterogeneity and anisotropy on GDC-driven dissolution rate by using numerical simulations with randomly obtained permeability fields. Dissolution rates calculated by these simulations are related to properties of the permeability field by using least-squares regression. We obtained two empirical formulas for predicting the asymptotic GDC-driven dissolution rate. In the first formula the dissolution rate is almost linearly proportional to the dimensionless equivalent vertical permeability. In the second one the dissolution rate is linearly proportional to a dimensionless vertical finger-tip velocity. This indicates that the GDC-driven dissolution can be predicted using either the equivalent vertical permeability or the vertical finger-tip velocity. Furthermore, both formulas demonstrate that higher anisotropy results in lower dissolution rates, suggesting that pronounced horizontal stratification can inhibit the dissolution of CO2.
Status: closed
-
RC1: 'Comment on egusphere-2024-3511', Ming Yang, 19 Dec 2024
Review of paper “Numerical Analysis of the Effect of Heterogeneity on CO2 Dissolution Enhanced by Gravity Driven Convection” submitted to HESS
General comments:
It is a pleasure procedure to review this interesting manuscript, which deals with the dissolution trapping of CO2 during geological carbon sequestration, an important means to reduce carbon emissions. In the process of sequestering CO2 in the deep saline aquifer, a concern is that the CO2 is less dense than the saline water, which gives rise to the possibility of upward CO2 leak. It is noticed that dissolution of CO2 in the brine is a crucial mechanism that reduces the risk of leak by transforming the supercritical CO2 into aqueous CO2. The significant challenge lies in the fact that the dissolution process is susceptible to instabilities driven by gravity and permeability heterogeneity, both of which are ubiquitous and involve big uncertainties. Therefore, it is important to characterize the dissolution rate of CO2 in the saline water, which is the objective of this manuscript.
Through reading the manuscript, I can feel that the authors have spent a big effort in preparing this work. They carefully selected the parameters for their numerical model based on the data from a deep literature review. They have conducted a large quantity of numerical simulations to build a data pool for systematic analysis. Finally, they give fruitful analysis and discussions of the results. Based on my own reading of the manuscript, I give my support of publication of this work on HESS. Several revisions are needed before it is accepted for publication.
Major comments:
(1)This manuscript improves the current predictor for enhanced CO2 dissolution due to gravity driven convection. Most of current predictors for enhanced CO2 dissolution address the homogeneous cases, which may have limitations for the real heterogeneous problem. Two representative works, as shown in Table 5 and Figure 6, have offered the preliminary predictor for the CO2 dissolution in heterogeneous fields. However, while one of them simplifies the predictor by neglecting the anisotropic effect, the other may not properly incorporate the ansotropic effect. The authors provided a new explanation of the anisotropic effect on the density driven dissolution. Their numerical results show that for a fixed vertical equivalent permeability increasing the horizontal equivalent permeability can reduce the dissolution efficiency, because increased horizontal permeability can increase horizontal mass exchange, make the horizontal mass distribution more uniform, and thus reduce the instability. My question is: I realize that the gamma_1=0.08 obtained by the data regression is quite similar to that for the 0.09 listed in Table 1. Do they have any relations? This is based on my observation that the alpha_1=1.1 is very close to 1.0 for the homogeneous.
(2)Furthermore, the authors introduced a new predictor using finger velocity, which is particularly intriguing because finger velocity can be measured using optical fiber technology. Given the rapid advancements and expanding applications of optical fibers in the field of geosciences, this novel formula could serve as a valuable tool for monitoring the trapping of CO2 through dissolution. My question is: The regression value for gamma_2=0.34 is quite different from 0.08 or 0.09, could you please explain why the gamma_2 is so different from gamma_1? I am wondering if it is simply a result of the data regression or maybe it have physical explanations.
(3)In the study of density driven instability, the fingers are usually irregular, as shown in the Figure 4 of your manuscript. Could you please explain why the fingers in Figure 3 are quite uniform?
(4)I can understand that the authors and many other researchers use two-dimensional model in both laboratory and numerical simulation, because of the high cost of three-dimensional model. I do agree that the two-dimensional study is very useful, but it would be nice if the authors can provide some discussions on the three-dimensional effect on GDC. I know it is hard to design a new simulation of new parameters, so I would appreciate it if the authors could find some related three-dimensional studies and give a short comparison of the difference of three dimensional simulations and two dimensional simulations. This may give more confidence for readers using the results from this work.
Minor comments:
1-The supplementary materials provided important and comprehensive information about the numerical modeling, but it is not fully referred in the paper. Please, give more clear reference of the supplementary materials in the paper.
2-Line 195 and 205, the appendix is missing.
3-Line 294, it would be better if we say in panel (a) of Figure 4 rather than in the first panel of Figure 4.
4-Line 299, I am wondering if you added white randomness (‘white noise’ maybe) in the heterogeneous fields?
5-Line 306, please remove the comma in ‘dissolution process, shows’.
6-Line 388, please remove the period before the references.
7-The overall structure of this manuscript is very nice, dividing the whole article into logical sections of proper titles. However, it may be better if the authors can reorganize the section ‘6 Conclusions’. We can see that in the section ‘5 Results and Discussion’. The authors first describe the general impact of heterogeneity on the development of instability, and then perform log-linear regressions of the simulation results. However, in the conclusion section the authors do not organize these results in the same order. Moreover, it would be better if the first paragraph is split so that we have a short summary of this work before writing the conclusions. This does not affect the comprehending of this work, but it would be nicer if the authors can reorganize the conclusions.
Ming Yang,
Tsinghua University-
AC1: 'Reply on RC1', yufei Wang, 25 Feb 2025
Dear Dr. Yang,
We sincerely appreciate your expertise and attention in evaluating our manuscript. Your comprehensive analysis and constructive feedback have enhanced the rigor and clarity of our work.Attached to this letter, you will find a point-by-point response to each of your comments in a separate PDF file. Should any aspect of our response require further clarification or adjustment, please do not hesitate to contact us.
Respectfully,
Yufei Wang
On behalf of all authors
-
AC1: 'Reply on RC1', yufei Wang, 25 Feb 2025
-
RC2: 'Comment on egusphere-2024-3511', Anonymous Referee #2, 23 Jan 2025
General comments
The article addresses the dissolution trapping of CO2 during geological carbon sequestration in deep saline aquifers, a critical process that serves as a key indicator for evaluating the safety of long-term CO2 storage. By trapping CO2 in the aqueous phase, dissolution trapping helps prevent CO2 leakage. However, the dissolution process can result in gravity-driven instabilities, such as Gravity-Driven Convection (GDC), and is influenced by formation properties like permeability and porosity, which are subject to uncertainties due to heterogeneities. The article aims to quantitatively analyze the explicit relationships between GDC-driven dissolution rates and permeability heterogeneities reported in the literature through numerical experiments.
The article is well-written and easy to understand. It is evident that the authors have carefully guided readers step by step through their work, supported by an in-depth literature review. Additionally, they have conducted a substantial number of numerical simulations across various setups to provide robust evidence. Based on my assessment of the manuscript, I believe this work is highly deserving of publication, subject to revisions.
Major comments
These are primarily comments from me and do not involve major revisions to the article.
- The 2D numerical setup is well-chosen, in my opinion, and all choices are clearly explained. Simulation data on a 100x100 grid over an extended period are obtained from 554 realizations of the permeability distributions. Additionally, the time-step size is governed by CFL conditions, which can result in a significant number of time steps. I assume this process is computationally intensive, and I was wondering how challenging it would be to extend this work to 3D, as doing so could provide valuable insights. Perhaps providing an approximate value of the CPU time per simulation could serve as a good starting point and indicator.
- The authors clearly state the software used for the numerical experiments and provide open access to the code and data. The user guide is extremely valuable, as it includes all the numerical details of the solver. However, I believe that for full reproducibility of this specific work, it would benefit from an additional short note, which I may have overlooked.
- Multiple realizations of permeability distributions are provided, and Figure 4 effectively illustrates the differences in numerical simulation results between them. I was wondering whether the upscaling of highly heterogeneous permeability fields could yield similar results for the predictors or if a relationship could be established between the upscaled and original permeability fields.
Minor comments
- Table 3: Could you comment on the potential effects of introducing a relationship between porosity and permeability realizations? Specifically, how might such a relationship influence the dissolution rates.
- Line 231: 553 realizations and line 354: 554 realizations.
- In Figure 4, the range of permeability values appears to be relatively narrow (from -14 to -10 on a logarithmic scale). Could you comment on the potential impact of using a wider range of permeability values?
Citation: https://doi.org/10.5194/egusphere-2024-3511-RC2 -
AC2: 'Reply on RC2', yufei Wang, 25 Feb 2025
Dear Reviewer,
We sincerely appreciate your expertise and attention in evaluating our manuscript. Your comprehensive analysis and constructive feedback have enhanced the rigor and clarity of our work.
Attached to this letter, you will find a point-by-point response to each of your comments in a separate PDF file. Should any aspect of our response require further clarification or adjustment, please do not hesitate to contact us.
Respectfully,
Yufei Wang
On behalf of all authors
-
CC1: 'Comment on egusphere-2024-3511', Giacomo Medici, 28 Jan 2025
General comments
Very good research on geological carbon storage. Please, take into account my comments to improve the manuscript.
Specific comments
Lines 15-87. Please, more emphasis on the fact that your research has a double focus: anisotropy and heterogeneity. My suggestions aim to increase the impact of your research in different sub-fields of applied geoscience.
Lines 21-22. “By being sealed under the low permeability caprock”. Consider inserting recent papers on low permeability sedimentary layers that incorporate a discussion on applications for CO2 storage.
- English, K.L., English, J.M., Moscardini, R., Haughton, P.D., Raine, R.J. and Cooper, M., 2024. Review of Triassic Sherwood Sandstone Group reservoirs of Ireland and Great Britain and their future role in geoenergy applications. Geoenergy, 2(1), pp.geoenergy2023-042.
- Medici, G., Munn, J.D., Parker, B.L., 2024. Delineating aquitard characteristics within a Silurian dolostone aquifer using high-density hydraulic head and fracture datasets. Hydrogeology Journal, 32(6), 1663-1691.
- Newell, A.J. and Shariatipour, S.M., 2016. Linking outcrop analogue with flow simulation to reduce uncertainty in sub-surface carbon capture and storage: an example from the Sherwood Sandstone Group of the Wessex Basin, UK. Geological Society, London, Special Publications, 436(1), 231-246.
Line 87. Please, disclose the specific objectives of your research by using numbers (e.g., i, ii, and iii).
Line 133. Be more specific when you mention anisotropy. “Flow anisotropy”? We are not talking about a mechanical anisotropy here due to the presence of rock discontinuities.
Line 295. Too soon for introducing references in this paragraph.
Line 367. “Geological Carbon Sequestration (GCS)”. You have already defined the acronym.
Line 424. Integrate the two relevant papers on the effect of low permeability layers on deep aquifers/reservoirs.
Figures and tables
Figure 1a. You made bigger part of the reservoir with a rectangle. The spatial scale is unclear.
Figure 2. Spatial scale unclear also here.
Figure 4. Make the letters larger.
Figure 4. Make the letters and numbers larger.
Citation: https://doi.org/10.5194/egusphere-2024-3511-CC1 -
AC3: 'Reply on CC1', yufei Wang, 25 Feb 2025
Dear Dr. Medici,
We sincerely appreciate your expertise and attention in evaluating our manuscript. Your comprehensive analysis and constructive feedback have enhanced the rigor and clarity of our work.
Attached to this letter, you will find a point-by-point response to each of your comments in a separate PDF file. Should any aspect of our response require further clarification or adjustment, please do not hesitate to contact us.
Respectfully,
Yufei Wang
On behalf of all authors
-
AC3: 'Reply on CC1', yufei Wang, 25 Feb 2025
Status: closed
-
RC1: 'Comment on egusphere-2024-3511', Ming Yang, 19 Dec 2024
Review of paper “Numerical Analysis of the Effect of Heterogeneity on CO2 Dissolution Enhanced by Gravity Driven Convection” submitted to HESS
General comments:
It is a pleasure procedure to review this interesting manuscript, which deals with the dissolution trapping of CO2 during geological carbon sequestration, an important means to reduce carbon emissions. In the process of sequestering CO2 in the deep saline aquifer, a concern is that the CO2 is less dense than the saline water, which gives rise to the possibility of upward CO2 leak. It is noticed that dissolution of CO2 in the brine is a crucial mechanism that reduces the risk of leak by transforming the supercritical CO2 into aqueous CO2. The significant challenge lies in the fact that the dissolution process is susceptible to instabilities driven by gravity and permeability heterogeneity, both of which are ubiquitous and involve big uncertainties. Therefore, it is important to characterize the dissolution rate of CO2 in the saline water, which is the objective of this manuscript.
Through reading the manuscript, I can feel that the authors have spent a big effort in preparing this work. They carefully selected the parameters for their numerical model based on the data from a deep literature review. They have conducted a large quantity of numerical simulations to build a data pool for systematic analysis. Finally, they give fruitful analysis and discussions of the results. Based on my own reading of the manuscript, I give my support of publication of this work on HESS. Several revisions are needed before it is accepted for publication.
Major comments:
(1)This manuscript improves the current predictor for enhanced CO2 dissolution due to gravity driven convection. Most of current predictors for enhanced CO2 dissolution address the homogeneous cases, which may have limitations for the real heterogeneous problem. Two representative works, as shown in Table 5 and Figure 6, have offered the preliminary predictor for the CO2 dissolution in heterogeneous fields. However, while one of them simplifies the predictor by neglecting the anisotropic effect, the other may not properly incorporate the ansotropic effect. The authors provided a new explanation of the anisotropic effect on the density driven dissolution. Their numerical results show that for a fixed vertical equivalent permeability increasing the horizontal equivalent permeability can reduce the dissolution efficiency, because increased horizontal permeability can increase horizontal mass exchange, make the horizontal mass distribution more uniform, and thus reduce the instability. My question is: I realize that the gamma_1=0.08 obtained by the data regression is quite similar to that for the 0.09 listed in Table 1. Do they have any relations? This is based on my observation that the alpha_1=1.1 is very close to 1.0 for the homogeneous.
(2)Furthermore, the authors introduced a new predictor using finger velocity, which is particularly intriguing because finger velocity can be measured using optical fiber technology. Given the rapid advancements and expanding applications of optical fibers in the field of geosciences, this novel formula could serve as a valuable tool for monitoring the trapping of CO2 through dissolution. My question is: The regression value for gamma_2=0.34 is quite different from 0.08 or 0.09, could you please explain why the gamma_2 is so different from gamma_1? I am wondering if it is simply a result of the data regression or maybe it have physical explanations.
(3)In the study of density driven instability, the fingers are usually irregular, as shown in the Figure 4 of your manuscript. Could you please explain why the fingers in Figure 3 are quite uniform?
(4)I can understand that the authors and many other researchers use two-dimensional model in both laboratory and numerical simulation, because of the high cost of three-dimensional model. I do agree that the two-dimensional study is very useful, but it would be nice if the authors can provide some discussions on the three-dimensional effect on GDC. I know it is hard to design a new simulation of new parameters, so I would appreciate it if the authors could find some related three-dimensional studies and give a short comparison of the difference of three dimensional simulations and two dimensional simulations. This may give more confidence for readers using the results from this work.
Minor comments:
1-The supplementary materials provided important and comprehensive information about the numerical modeling, but it is not fully referred in the paper. Please, give more clear reference of the supplementary materials in the paper.
2-Line 195 and 205, the appendix is missing.
3-Line 294, it would be better if we say in panel (a) of Figure 4 rather than in the first panel of Figure 4.
4-Line 299, I am wondering if you added white randomness (‘white noise’ maybe) in the heterogeneous fields?
5-Line 306, please remove the comma in ‘dissolution process, shows’.
6-Line 388, please remove the period before the references.
7-The overall structure of this manuscript is very nice, dividing the whole article into logical sections of proper titles. However, it may be better if the authors can reorganize the section ‘6 Conclusions’. We can see that in the section ‘5 Results and Discussion’. The authors first describe the general impact of heterogeneity on the development of instability, and then perform log-linear regressions of the simulation results. However, in the conclusion section the authors do not organize these results in the same order. Moreover, it would be better if the first paragraph is split so that we have a short summary of this work before writing the conclusions. This does not affect the comprehending of this work, but it would be nicer if the authors can reorganize the conclusions.
Ming Yang,
Tsinghua University-
AC1: 'Reply on RC1', yufei Wang, 25 Feb 2025
Dear Dr. Yang,
We sincerely appreciate your expertise and attention in evaluating our manuscript. Your comprehensive analysis and constructive feedback have enhanced the rigor and clarity of our work.Attached to this letter, you will find a point-by-point response to each of your comments in a separate PDF file. Should any aspect of our response require further clarification or adjustment, please do not hesitate to contact us.
Respectfully,
Yufei Wang
On behalf of all authors
-
AC1: 'Reply on RC1', yufei Wang, 25 Feb 2025
-
RC2: 'Comment on egusphere-2024-3511', Anonymous Referee #2, 23 Jan 2025
General comments
The article addresses the dissolution trapping of CO2 during geological carbon sequestration in deep saline aquifers, a critical process that serves as a key indicator for evaluating the safety of long-term CO2 storage. By trapping CO2 in the aqueous phase, dissolution trapping helps prevent CO2 leakage. However, the dissolution process can result in gravity-driven instabilities, such as Gravity-Driven Convection (GDC), and is influenced by formation properties like permeability and porosity, which are subject to uncertainties due to heterogeneities. The article aims to quantitatively analyze the explicit relationships between GDC-driven dissolution rates and permeability heterogeneities reported in the literature through numerical experiments.
The article is well-written and easy to understand. It is evident that the authors have carefully guided readers step by step through their work, supported by an in-depth literature review. Additionally, they have conducted a substantial number of numerical simulations across various setups to provide robust evidence. Based on my assessment of the manuscript, I believe this work is highly deserving of publication, subject to revisions.
Major comments
These are primarily comments from me and do not involve major revisions to the article.
- The 2D numerical setup is well-chosen, in my opinion, and all choices are clearly explained. Simulation data on a 100x100 grid over an extended period are obtained from 554 realizations of the permeability distributions. Additionally, the time-step size is governed by CFL conditions, which can result in a significant number of time steps. I assume this process is computationally intensive, and I was wondering how challenging it would be to extend this work to 3D, as doing so could provide valuable insights. Perhaps providing an approximate value of the CPU time per simulation could serve as a good starting point and indicator.
- The authors clearly state the software used for the numerical experiments and provide open access to the code and data. The user guide is extremely valuable, as it includes all the numerical details of the solver. However, I believe that for full reproducibility of this specific work, it would benefit from an additional short note, which I may have overlooked.
- Multiple realizations of permeability distributions are provided, and Figure 4 effectively illustrates the differences in numerical simulation results between them. I was wondering whether the upscaling of highly heterogeneous permeability fields could yield similar results for the predictors or if a relationship could be established between the upscaled and original permeability fields.
Minor comments
- Table 3: Could you comment on the potential effects of introducing a relationship between porosity and permeability realizations? Specifically, how might such a relationship influence the dissolution rates.
- Line 231: 553 realizations and line 354: 554 realizations.
- In Figure 4, the range of permeability values appears to be relatively narrow (from -14 to -10 on a logarithmic scale). Could you comment on the potential impact of using a wider range of permeability values?
Citation: https://doi.org/10.5194/egusphere-2024-3511-RC2 -
AC2: 'Reply on RC2', yufei Wang, 25 Feb 2025
Dear Reviewer,
We sincerely appreciate your expertise and attention in evaluating our manuscript. Your comprehensive analysis and constructive feedback have enhanced the rigor and clarity of our work.
Attached to this letter, you will find a point-by-point response to each of your comments in a separate PDF file. Should any aspect of our response require further clarification or adjustment, please do not hesitate to contact us.
Respectfully,
Yufei Wang
On behalf of all authors
-
CC1: 'Comment on egusphere-2024-3511', Giacomo Medici, 28 Jan 2025
General comments
Very good research on geological carbon storage. Please, take into account my comments to improve the manuscript.
Specific comments
Lines 15-87. Please, more emphasis on the fact that your research has a double focus: anisotropy and heterogeneity. My suggestions aim to increase the impact of your research in different sub-fields of applied geoscience.
Lines 21-22. “By being sealed under the low permeability caprock”. Consider inserting recent papers on low permeability sedimentary layers that incorporate a discussion on applications for CO2 storage.
- English, K.L., English, J.M., Moscardini, R., Haughton, P.D., Raine, R.J. and Cooper, M., 2024. Review of Triassic Sherwood Sandstone Group reservoirs of Ireland and Great Britain and their future role in geoenergy applications. Geoenergy, 2(1), pp.geoenergy2023-042.
- Medici, G., Munn, J.D., Parker, B.L., 2024. Delineating aquitard characteristics within a Silurian dolostone aquifer using high-density hydraulic head and fracture datasets. Hydrogeology Journal, 32(6), 1663-1691.
- Newell, A.J. and Shariatipour, S.M., 2016. Linking outcrop analogue with flow simulation to reduce uncertainty in sub-surface carbon capture and storage: an example from the Sherwood Sandstone Group of the Wessex Basin, UK. Geological Society, London, Special Publications, 436(1), 231-246.
Line 87. Please, disclose the specific objectives of your research by using numbers (e.g., i, ii, and iii).
Line 133. Be more specific when you mention anisotropy. “Flow anisotropy”? We are not talking about a mechanical anisotropy here due to the presence of rock discontinuities.
Line 295. Too soon for introducing references in this paragraph.
Line 367. “Geological Carbon Sequestration (GCS)”. You have already defined the acronym.
Line 424. Integrate the two relevant papers on the effect of low permeability layers on deep aquifers/reservoirs.
Figures and tables
Figure 1a. You made bigger part of the reservoir with a rectangle. The spatial scale is unclear.
Figure 2. Spatial scale unclear also here.
Figure 4. Make the letters larger.
Figure 4. Make the letters and numbers larger.
Citation: https://doi.org/10.5194/egusphere-2024-3511-CC1 -
AC3: 'Reply on CC1', yufei Wang, 25 Feb 2025
Dear Dr. Medici,
We sincerely appreciate your expertise and attention in evaluating our manuscript. Your comprehensive analysis and constructive feedback have enhanced the rigor and clarity of our work.
Attached to this letter, you will find a point-by-point response to each of your comments in a separate PDF file. Should any aspect of our response require further clarification or adjustment, please do not hesitate to contact us.
Respectfully,
Yufei Wang
On behalf of all authors
-
AC3: 'Reply on CC1', yufei Wang, 25 Feb 2025
Viewed
Since the preprint corresponding to this journal article was posted outside of Copernicus Publications, the preprint-related metrics are limited to HTML views.
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
181 | 0 | 0 | 181 | 0 | 0 |
- HTML: 181
- PDF: 0
- XML: 0
- Total: 181
- BibTeX: 0
- EndNote: 0
Viewed (geographical distribution)
Since the preprint corresponding to this journal article was posted outside of Copernicus Publications, the preprint-related metrics are limited to HTML views.
Country | # | Views | % |
---|---|---|---|
United States of America | 1 | 80 | 45 |
China | 2 | 36 | 20 |
Israel | 3 | 7 | 3 |
Norway | 4 | 7 | 3 |
Brazil | 5 | 5 | 2 |
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1
- 80