Edith Cowan University Research Online
Theses: Doctorates and Masters Theses
2019
Wettability of quartz surfaces under carbon dioxide geo- sequestration conditions. A theoretical study
Aleksandr Abramov Edith Cowan University
Follow this and additional works at: https://ro.ecu.edu.au/theses
Part of the Engineering Commons
Recommended Citation Abramov, A. (2019). Wettability of quartz surfaces under carbon dioxide geo-sequestration conditions. A theoretical study. https://ro.ecu.edu.au/theses/2232
This Thesis is posted at Research Online. https://ro.ecu.edu.au/theses/2232 Edith Cowan University
Copyright Warning
You may print or download ONE copy of this document for the purpose of your own research or study.
The University does not authorize you to copy, communicate or otherwise make available electronically to any other person any copyright material contained on this site.
You are reminded of the following:
Copyright owners are entitled to take legal action against persons who infringe their copyright.
A reproduction of material that is protected by copyright may be a copyright infringement. Where the reproduction of such material is done without attribution of authorship, with false attribution of authorship or the authorship is treated in a derogatory manner, this may be a breach of the author’s moral rights contained in Part IX of the Copyright Act 1968 (Cth).
Courts have the power to impose a wide range of civil and criminal sanctions for infringement of copyright, infringement of moral rights and other offences under the Copyright Act 1968 (Cth). Higher penalties may apply, and higher damages may be awarded, for offences and infringements involving the conversion of material into digital or electronic form. Wettability of Quartz Surfaces under Carbon Dioxide Geo-Sequestration Conditions. A Theoretical Study
Aleksandr Abramov
Submitted for the Degree of Doctor of Philosophy
Edith Cowan University
School of Engineering
2019 Abstract
The wettability of rocks under reservoir conditions is important to ensure and secure long term underground storage of carbon dioxide. The composition of those rocks vary significantly and are influenced by the fact that quartz is the second most abundant mineral in the earth's continental crust. Thus, the CO2 wettability of quartz dominates the overall CO2 trapping performance of storage and cap rocks. If depleted oil or gas reservoirs are used for storage of
CO2 quartz surfaces of rocks in reservoirs which have been previously exposed to hydrocarbons might be covered with chemisorpt hydrocarbon molecules. The CO2 wettability of these chemically modified quartz is studied in this work with molecular dynamics.
To model quartz surfaces with chemisorpt hydrocarbons both CLAYFF and DREIDING force fields are coupled at atomic site charge level using the density functional theory and the Bader charge analysis. Augmented with modified charges of the OC bond, CLAYFF and the DREIDING force fields are applied to solve the practical problem of calculating the contact angle of a water droplet on alkylated quartz surfaces in a CO2 environment. A systematic computational study of wettability of fully hydroxylated and alkylated (001) -quartz surface in carbon dioxide atmosphere with respect to surface concentration of pentyl groups is performed. Alkylated quartz surfaces have been shown to be extremely hydrophobic even when the surface density of hydroxyl groups is close to the highest naturally observed. The study also verifies that a comprehensive description of wettability of alkylated quartz surface requires three parameters: the theoretical contact angle, the apparent contact angle and the hidden contact angle. These contact angles are determined at the tip level of pentyl groups and the level of the quartz surface. The hidden contact angle is calculated as the angle of a water "skirt", which is formed between the level of the quartz surface and the tip level of pentyl groups.
Additionally, the concept and the method of how to determine computational contact angles of a liquid droplet resting on a solid surface from individual snapshots of molecular dynamics simulations have been formulated, implemented and analysed in this work. Spherical coordinates to circumscribe a sphere around given configuration of water molecules form the basis of the method, which is thus natural and consistent with the droplet's geometric computational framework.
2 Declaration I certify that this thesis does not, to the best of my knowledge and belief:
i. incorporate without acknowledgment any material previously submitted for a degree or diploma in any institution of higher education; ii. contain any material previously published or written by another person except where due reference is made in the text of this thesis; or iii. contain any defamatory material.
Aleksandr Abramov 27/07/2019
3 Acknowledgments This work was supported by Edith Cowan University Petroleum Engineering PhD Scholarship and by resources provided by the Pawsey Supercomputing Centre with funding from the Australian Government and the Government of Western Australia.
4 Table of Contents List of tables ...... 8
List of figures ...... 9
1. Introduction...... 22
1.1. Scope and objective ...... 22
1.2. Methods, procedures and process ...... 22
1.3. Results, observations and conclusions ...... 22
1.4. Novel and additive information ...... 22
1.5. Thesis organization ...... 23
1.6. Publications ...... 23
2. Literature review ...... 25
2.1. Global warming ...... 25
2.1.1. Overview ...... 25
2.1.2. Response strategies ...... 29
2.2. Carbon capture and storage ...... 30
2.2.1. Overview ...... 30
2.2.2. Carbon geo-sequestration projects ...... 32
2.3. Wettability and carbon dioxide storage capacity of rocks ...... 34
2.4. Properties of carbon dioxide in geological conditions ...... 37
2.5. Quartz wettability in presence of carbon dioxide ...... 40
2.5.1. Experimental studies ...... 41
2.5.2. Theoretical studies ...... 47
2.6. Quartz crystal ...... 49
2.7. Quartz surface ...... 51
2.8. Conclusions ...... 56
3. Computational methods ...... 57
5 3.1. Overview ...... 57
3.2. Force fields ...... 57
3.3. Molecular dynamics ...... 63
3.4. Density functional theory ...... 67
3.5. Conclusions ...... 69
4. The suitability of spheroidal constructions to estimate the contact angle using snapshots of molecular dynamics simulations ...... 71
4.1. Annotation ...... 71
4.2. Introduction and motivation ...... 71
4.3. Method development ...... 73
4.4. Computational details ...... 79
4.4.1. Force fields...... 79
4.4.2. Reconstruction of the pristine quartz surface ...... 81
4.4.3. Simulation setups ...... 82
4.5. Method application and discussion ...... 85
4.6. Conclusions ...... 90
5. Application of the CLAYFF and the DREIDING force fields for modelling of alkylated quartz surfaces ...... 91
5.1. Annotation ...... 91
5.2. Introduction ...... 91
5.3. Computational methodology ...... 93
5.4. Coupling of the CLAYFF and the DREIDING force fields ...... 97
5.5. Application of modified force fields to alkylated quartz surfaces ...... 102
5.6. Conclusions ...... 109
6. Wettability of fully hydroxylated and alkylated (001) alpha-quartz surface in carbon dioxide atmosphere ...... 111
6 6.1. Annotation ...... 111
6.2. Introduction ...... 111
6.3. Computational methodology ...... 114
6.4. Force fields ...... 115
6.5. Modelling of the oxygen-carbon bond ...... 118
6.6. Preparation of alkylated quartz surfaces ...... 120
6.7. Simulation setups ...... 123
6.8. Determination of the contact angle ...... 125
6.9. Qualitative analysis of the results ...... 127
6.10. Quantitative analysis of the results ...... 132
6.11. Conclusions...... 138
6.12. Supplementary material ...... 140
6.12.1. Surface hydroxyl group concentration vs surface pentyl group concentration 140
6.12.2. Water iso-density charts for surfaces 1 through 11 ...... 140
6.12.3. Simulation snapshots for surfaces 1 through 11 ...... 151
6.12.4. Simulation snapshots with removed pentyl groups for surfaces 1 through 11 162
7. Conclusions...... 168
8. References ...... 170
7 List of tables
Table 1. The Buckingham potential parameters used for the pristine quartz surface (van Beest et al., 1990)...... 80
Table 2. The Lennard-Jones potential parameters used to model CO2-H2O-SiO2 systems (Abascal & Vega, 2005; Harris & Yung, 1995; S. Iglauer et al., 2012; van Beest et al., 1990). . 80
Table 3. Summary of computational details...... 84
Table 4. Contact angle computational results...... 85
Table 5. Non-bond potential parameters used in the simulations (Cygan et al., 2004; Mayo et al., 1990)...... 94
Table 6. Bond stretch parameters used in the simulations (Cygan et al., 2004; Mayo et al., 1990)...... 96
Table 7. Non-bond potential parameters used in the simulations (Cygan et al., 2004; Mayo et al., 1990)...... 116
Table 8. Bond stretch parameters used in the simulations (Cygan et al., 2004; Mayo et al., 1990)...... 117
Table 9. Wettability zones according to the theoretical and the apparent contact angles... 135
8 List of figures
Figure 1. Emissions of carbon dioxide due to industrial and agricultural activities (1 Petagram of carbon = 1 PgC = 1015 grams of carbon = 1 Gigatonne of carbon = 1 GtC). Image source: (IPCC, 2013)...... 26
Figure 2. Atmospheric CO2 concentration, CO2 partial pressure (pCO2, left axis) and pH of oceanic surface. MLO: Mauna Loa Observatory, Hawaii; SPO: South Pole; HOT: Hawaii Ocean Time-Series station. Image source: (IPCC, 2013)...... 27
Figure 3. Climate change indicators. Image source: (IPCC, 2013)...... 28
Figure 4. Schematic illustration of carbon dioxide capture and storage technology. Geological storage formations may vary and can be selected from specially dedicated for storage formations, depleted oil or gas reservoirs, deep aquifers, coal beds or salt caverns...... 31
Figure 5. A liquid droplet on a solid surface and variables used in the Young equation...... 34
Figure 6. Capillary rise of a liquid illustrating parameters used in the Young-Laplace equation...... 35
Figure 7. Schematic illustration of CO2 structural trapping...... 36
Figure 8. Phase diagram of carbon dioxide. Image source: ChemicaLogic Corporation, Copyright (c) 1999 ChemicaLogic Corporation, USA. All rights reserved...... 37
Figure 9. CO2 density as a function of temperature and pressure. Image source: (Bachu, 2003; IPCC, 2005)...... 39
Figure 10. Possible geological pressure and temperature conditions in "cold" and "warm" basins worldwide superimposed onto carbon dioxide phase diagram. Thickened lines highlight conditions at optimal depths 800-1000 m for "cold" and 1500-2000 m for "warm" basins, respectively. Interpolation between optimal "cold" and "warm" conditions (gradient coloured lines) shows range of pressures and temperatures expected at CO2 storage depths. Image is constructed on basis of (Bachu, 2003)...... 40
Figure 11. Liquid droplet on a rough and inhomogeneous surface. The Wenzel (left) and the Cassie-Baxter (right) models...... 43
9 Figure 12. Advancing contact angle of a deionized water droplet on rough quartz surface (RMS 560 nm) cleaned with air plasma in carbon dioxide atmosphere at different conditions of pressure (in MPa) and temperature (in Kelvin). Source of data: (Ahmed Z. Al-Yaseri et al., 2016)...... 45
Figure 13. Advancing contact angle of a water droplet on rough quartz surface (RMS 560 nm) cleaned with air plasma with varying concentration of NaCl (in wt%) at 10 MPa of carbon dioxide pressure at two temperatures (in Kelvin). Source of data: (Ahmed Z. Al-Yaseri et al., 2016)...... 46
Figure 14. Approximate phase diagram of silica. Image source: (serc.carleton.edu), drawn by Dexter Perkins and John Brady on basis of (Swamy et al., 1994)...... 50
Figure 15. The primitive unit cell of α-quartz. Yellow balls are silicon atoms; red balls are oxygen atoms...... 51
Figure 16. Top view of hydroxylated (001) quartz surface. Only top SiO4H2 groups are shown. Yellow balls are silicon atoms, red balls are oxygen atoms, white balls are hydrogen atoms. Short (strong) hydrogen bonds are shown with straight dashed lines. To highlight the zigzag structure of alternating strong and weak hydrogen bonds three long (weak) hydrogen bonds are shown with dashed blue ellipses along a horizontal line...... 53
Figure 17. Side views (top: along x direction; bottom: along y direction) of a four-layer slab (four heights of the primitive quartz unit cell) of hydroxylated quartz. Atoms are shown with bonds, yellow - silicon atoms, red - oxygen atoms, white - hydrogen atoms...... 54
Figure 18. Top view of pristine (001) quartz surface. Only top layer atoms are shown. Yellow balls are silicon atoms and red balls are oxygen atoms...... 55
Figure 19 Side views (top: along x direction; bottom: along y direction) of a four-layer slab (four heights of the primitive quartz unit cell) of pristine quartz. Atoms are shown with bonds, yellow - silicon atoms, red - oxygen atoms...... 56
Figure 20. Point charges and auxiliary charge distributions in the Ewald sum...... 62
Figure 21. Simulation cell with schematic illustration of hydrophobic surface, water droplet and steps of the scanning algorithm...... 74
10 Figure 22. Simulation cell with schematic illustration of a hydrophilic surface, showing the water droplet, steps and elements of the scanning algorithm...... 75
Figure 23. Calculation of water density in spherical coordinates for hydrophobic and hydrophilic surfaces...... 77
Figure 24. Illustration of contact angle calculations...... 78
Figure 25. Quartz unit cells used in this work. Left: primitive hexagonal unit cell of α-quartz, isometric view. Right: reconstruction of the orthorhombic cell (solid bold rectangle) from 2x2x1 super cell, top view (dashed parallelogram shows the primitive cell). Red balls - oxygen atoms, blue balls - silicon atoms...... 81
Figure 26. Four-layer quartz slab used in this work. Left: view of the super cell 3x2x1 along z direction with only top atoms shown for clarity (the orthorhombic unit cell is highlighted with dashed rectangle). Right: slab views along x (right top) and y (right bottom) directions. Red balls - oxygen atoms, blue balls - silicon atoms...... 82
Figure 27. Initial simulation setups for systems investigated at two CO2 pressures - 4 MPa (left) and 10 MPa (right). The hemispherical water droplets above the quartz surfaces surrounded by CO2 molecules are clearly visible...... 84
Figure 28. Water density in spherical coordinates for the system at 4 MPa CO2 pressure. Horizontal line shows half of water's normal density...... 87
Figure 29. Water density in spherical coordinates for the system at 10 MPa CO2 pressure. Horizontal line shows half of water's normal density...... 87
Figure 30. Water droplet and its circumscribed sphere on the quartz surface at 4 MPa CO2 pressure, isometric view (left) and top view (right). Orange balls are oxygen atoms and white balls are hydrogen atoms of the TIP4P/2005 water...... 88
Figure 31. Water droplet and its circumscribed sphere on the quartz surface at 10 MPa CO2 pressure, isometric view (left) and top view (right). Orange balls are oxygen atoms and white balls are hydrogen atoms of the TIP4P/2005 water...... 88
Figure 32. Iso-density chart constructed in cylindrical coordinates for a water droplet on a pristine quartz surface at 4 MPa CO2 pressure. Blue line shows half of water normal density. Red dots depict data points used to fit the tangential line, the straight line (5, 7 and 9 points
11 were used for the averaging, the case for 7 points is shown). The dashed line shows the circumscribed sphere contour found in the spherical coordinates...... 89
Figure 33. Iso-density chart constructed in cylindrical coordinates for a water droplet on a pristine quartz surface at 10 MPa CO2 pressure. Blue line shows half of water normal density. Red dots depict data points used to fit the tangential line, the straight line (5, 7 and 9 points were used for the averaging, the case for 7 points is shown). The dashed line shows the circumscribed sphere contour found in the spherical coordinates...... 89
Figure 34. Conceptual illustration of the force fields coupling problem. White balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms, yellow balls - silicon atoms. Top and side views of the quartz surface/crystal are shown on the left, a side view of the alkylated quartz is shown on the right. The polar OC bond is highlighted with dashed ellipses. FF stands for Force Field...... 98
Figure 35. Optimized at the DFT level of theory quartz slab, side view and views along the x and y axis. White balls - hydrogen atoms, light blue balls - carbon atoms, red balls - oxygen atoms, yellow balls - silicon atoms. Blue lines demark the unit cell of the simulated system, images of the unit cell are shown with lines instead of the balls...... 99
Figure 36. Charge density distribution of the three-layer quartz slab computed at the DFT level of theory, views along x (left) and y (right) axises. White balls - hydrogen atoms, light blue balls - carbon atoms, red balls - oxygen atoms, yellow balls - silicon atoms. Grey halo shows the iso- surface of the charge density, the iso-value was chosen such that the charge clouds around the oxygen and silicon atoms just start to overlap...... 100
Figure 37. Absolute charge difference between modified force field and DFT results for the two atoms of the OC bond. FF stands for Force Field...... 101
Figure 38. Top layer of the surface unit cell of hydroxylated quartz with a pentyl surface group
2 concentration of 1.45 C5H11/nm . White balls - hydrogen atoms, light blue balls - carbon atoms, red balls - oxygen atoms, yellow balls - silicon atoms...... 103
Figure 39. Top layer of the surface unit cell of hydroxylated quartz with a pentyl surface group
2 concentration of 3.18 C5H11/nm . White balls - hydrogen atoms, light blue balls - carbon atoms, red balls - oxygen atoms, yellow balls - silicon atoms...... 104
12 Figure 40. Side views of the initial simulation setups for quartz surfaces with C5H11 density 1.45 (left) and 3.18 (right) groups per square nanometre. White balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms, yellow balls - silicon atoms. CO2 molecules are shown with lines to improve visibility of the water droplet...... 106
Figure 41. Circle fitted to the water iso-density for the hydroxylated quartz surface with a
2 pentyl concentration of 1.45 C5H11/nm . Dashed line shows the tip level of pentyl groups. 107
Figure 42. Circle fitted to the water iso-density for the hydroxylated quartz surface with a
2 pentyl concentration of 3.18 C5H11/nm . Dashed line shows the tip level of pentyl groups. 108
Figure 43. Simulation snapshot of the quartz surface with a pentyl concentration of 1.45 groups per square nm. View along x axis (left) and along y axis (right). The sphere fitted to the iso-density chart is illustrated in light purple color. Water molecules and pentyl groups are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms. Atoms of the CO2 molecules and of the quartz slab are removed for clarity...... 108
Figure 44. Simulation snapshot of the quartz surface with a pentyl concentration of 3.18 groups per square nm. View along x axis (left) and along y axis (right). The sphere fitted to the iso-density chart is illustrated in light purple color. Water molecules and pentyl groups are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms. Atoms of the CO2 molecules and of the quartz slab are removed for clarity...... 109
Figure 45. Optimized at the DFT level of theory three-layer quartz slab with integrated pentyl group. White balls - hydrogen atoms, light blue balls - carbon atoms, red balls - oxygen atoms, yellow balls - silicon atoms. Blue lines demark the unit cell of the simulated system, images of which are shown with lines instead of the balls...... 119
Figure 46. Top view of the unit cell of the top layer of the fully hydroxylated quartz surface.
White balls - hydrogen atoms, red balls - oxygen atoms, yellow balls - silicon atoms. Top SiO4H2 groups are shown with balls, lower atoms are shown with lines of the same colour. Dashed parallelogram demarks the 2x1x1 super cell of the orthorhombic unit cell. Numbers illustrate the cumulative alkylation sequence...... 121
13 Figure 47. Top views of the top layer unit cells. 16 quartz surfaces with surface pentyl concentrations ranging from 0.289 to 4.629 groups/nm2 (surfaces 1 to 16) are shown, unit cells are ordered from left to right and from bottom to top. White balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms, yellow balls - silicon atoms...... 122
Figure 48. Top view (left) and views along x (right top) and y (right bottom) directions of fully pentylated quartz surface (surface 16). White balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms, yellow balls - silicon atoms. Top SiO4HC5H11 groups are shown with balls, lower atoms are shown with lines of the same colour. Only top layer unit cell of four-layer slab is demonstrated...... 123
Figure 49. Space available to CO2 in the hydroxylated (left) and fully pentylated (right) systems.
Dashed red line demarks volume available to CO2. White balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms, yellow balls - silicon atoms...... 124
Figure 50. Initial simulation setup for the quartz surface with four C5H11 groups in the unit cell (surface 4). Views along x (left) and y (right) axis. White balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms, yellow balls - silicon atoms. To make the water droplet more visible CO2 molecules are shown with lines...... 125
Figure 51. Schematic illustration of calculations of the water iso-density...... 126
Figure 52. Water iso-density chart for quartz surface (surface 4) with concentration of pentyl groups 1.157 C5H11 per square nm. The iso-density data points are shown with blue dots, fitted circle is shown with blue line, the tip level of the pentyl groups is shown with the dashed line, the dotted line illustrates the planar projection of the water "skirt"...... 128
Figure 53. Averaged contours of water droplets for surfaces 1 (S1) through 11 (S11) with corresponding concentrations of pentyl groups per square nm (ascending water droplet). Dashed horizontal lines represent the lowest and the highest tip levels of the pentyl groups...... 129
Figure 54. Simulation snapshot of the quartz surface with pentyl concentration 2.893 groups per square nm (surface 10). View along x axis: top left; view along y axis: top right; view along z axis: bottom left; isometric view: bottom right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface is illustrated with square in light yellow color.
14 Water molecules and pentyl groups are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms...... 130
Figure 55. Side views of a water droplet on the quartz surface with removed pentyl groups for pentyl concentration 1.157 groups per square nm (surface 4). View along x axis: left; view along y axis: right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface position is indicated by the lowest laying water molecules. Water molecules are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms. .... 131
Figure 56. Parameters used to quantify the wettability of alkylated quartz surface. Water droplet and its "skirt" are shown in blue, C5H11 groups are illustrated with zigzags...... 132
Figure 57. Dependence of the circumscribed sphere (circle) radius on concentration of pentyl groups...... 132
Figure 58. Dependence of the vertical position of the circumscribed sphere (circle) centre on concentration of pentyl groups. Position of the quartz surface is at zero...... 133
Figure 59. Circumscribed sphere (circle) or spherical droplet contact with the quartz surface
(height zero) and the tip level of pentyl groups as a function of C5H11 concentration...... 134
Figure 60. Dependence of the theoretical contact angle on concentration of pentyl groups...... 134
Figure 61. Dependence of the apparent contact angle on concentration of pentyl groups. 135
Figure 62. Dependence of the hidden contact angle on concentration of pentyl groups. .... 136
Figure 63. Surface are per one C5H11 group as a function of its surface concentration. Dotted line averages values at the plateau...... 138
Figure 64. Surface hydroxyl group concentration as a function of surface pentyl group concentration...... 140
Figure 65. Water iso-density chart for quartz surface (surface 1) with concentration of pentyl groups 0.289 C5H11 per square nm. The iso-density data points are shown with blue dots, fitted circle is shown with blue line, the tip level of the pentyl groups is shown with the dashed line, the dotted line illustrates the planar projection of the water "skirt"...... 140
15 Figure 66. Water iso-density chart for quartz surface (surface 2) with concentration of pentyl groups 0.579 C5H11 per square nm. The iso-density data points are shown with blue dots, fitted circle is shown with blue line, the tip level of the pentyl groups is shown with the dashed line, the dotted line illustrates the planar projection of the water "skirt"...... 141
Figure 67. Water iso-density chart for quartz surface (surface 3) with concentration of pentyl groups 0.868 C5H11 per square nm. The iso-density data points are shown with blue dots, fitted circle is shown with blue line, the tip level of the pentyl groups is shown with the dashed line, the dotted line illustrates the planar projection of the water "skirt"...... 142
Figure 68. Water iso-density chart for quartz surface (surface 4) with concentration of pentyl groups 1.157 C5H11 per square nm. The iso-density data points are shown with blue dots, fitted circle is shown with blue line, the tip level of the pentyl groups is shown with the dashed line, the dotted line illustrates the planar projection of the water "skirt"...... 143
Figure 69. Water iso-density chart for quartz surface (surface 5) with concentration of pentyl groups 1.447 C5H11 per square nm. The iso-density data points are shown with blue dots, fitted circle is shown with blue line, the tip level of the pentyl groups is shown with the dashed line, the dotted line illustrates the planar projection of the water "skirt"...... 144
Figure 70. Water iso-density chart for quartz surface (surface 6) with concentration of pentyl groups 1.736 C5H11 per square nm. The iso-density data points are shown with blue dots, fitted circle is shown with blue line, the tip level of the pentyl groups is shown with the dashed line, the dotted line illustrates the planar projection of the water "skirt"...... 145
Figure 71. Water iso-density chart for quartz surface (surface 7) with concentration of pentyl groups 2.025 C5H11 per square nm. The iso-density data points are shown with blue dots, fitted circle is shown with blue line, the tip level of the pentyl groups is shown with the dashed line, the dotted line illustrates the planar projection of the water "skirt"...... 146
Figure 72. Water iso-density chart for quartz surface (surface 8) with concentration of pentyl groups 2.314 C5H11 per square nm. The iso-density data points are shown with blue dots, fitted circle is shown with blue line, the tip level of the pentyl groups is shown with the dashed line, the dotted line illustrates the planar projection of the water "skirt"...... 147
Figure 73. Water iso-density chart for quartz surface (surface 9) with concentration of pentyl groups 2.604 C5H11 per square nm. The iso-density data points are shown with blue dots, fitted
16 circle is shown with blue line, the tip level of the pentyl groups is shown with the dashed line, the dotted line illustrates the planar projection of the water "skirt"...... 148
Figure 74. Water iso-density chart for quartz surface (surface 10) with concentration of pentyl groups 2.893 C5H11 per square nm. The iso-density data points are shown with blue dots, fitted circle is shown with blue line, the tip level of the pentyl groups is shown with the dashed line...... 149
Figure 75. Water iso-density chart for quartz surface (surface 11) with concentration of pentyl groups 3.182 C5H11 per square nm. The iso-density data points are shown with blue dots, fitted circle is shown with blue line, the tip level of the pentyl groups is shown with the dashed line...... 150
Figure 76. Simulation snapshot of the quartz surface with pentyl concentration 0.289 groups per square nm (surface 1). View along x axis: top left; view along y axis: top right; view along z axis: bottom left; isometric view: bottom right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface is illustrated with square in light yellow color. Water molecules and pentyl groups are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms...... 151
Figure 77. Simulation snapshot of the quartz surface with pentyl concentration 0.579 groups per square nm (surface 2). View along x axis: top left; view along y axis: top right; view along z axis: bottom left; isometric view: bottom right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface is illustrated with square in light yellow color. Water molecules and pentyl groups are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms...... 152
Figure 78. Simulation snapshot of the quartz surface with pentyl concentration 0.868 groups per square nm (surface 3). View along x axis: top left; view along y axis: top right; view along z axis: bottom left; isometric view: bottom right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface is illustrated with square in light yellow color. Water molecules and pentyl groups are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms...... 153
Figure 79. Simulation snapshot of the quartz surface with pentyl concentration 1.157 groups per square nm (surface 4). View along x axis: top left; view along y axis: top right; view along
17 z axis: bottom left; isometric view: bottom right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface is illustrated with square in light yellow color. Water molecules and pentyl groups are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms...... 154
Figure 80. Simulation snapshot of the quartz surface with pentyl concentration 1.447 groups per square nm (surface 5). View along x axis: top left; view along y axis: top right; view along z axis: bottom left; isometric view: bottom right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface is illustrated with square in light yellow color. Water molecules and pentyl groups are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms...... 155
Figure 81. Simulation snapshot of the quartz surface with pentyl concentration 1.736 groups per square nm (surface 6). View along x axis: top left; view along y axis: top right; view along z axis: bottom left; isometric view: bottom right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface is illustrated with square in light yellow color. Water molecules and pentyl groups are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms...... 156
Figure 82. Simulation snapshot of the quartz surface with pentyl concentration 2.025 groups per square nm (surface 7). View along x axis: top left; view along y axis: top right; view along z axis: bottom left; isometric view: bottom right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface is illustrated with square in light yellow color. Water molecules and pentyl groups are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms...... 157
Figure 83. Simulation snapshot of the quartz surface with pentyl concentration 2.314 groups per square nm (surface 8). View along x axis: top left; view along y axis: top right; view along z axis: bottom left; isometric view: bottom right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface is illustrated with square in light yellow color. Water molecules and pentyl groups are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms...... 158
Figure 84. Simulation snapshot of the quartz surface with pentyl concentration 2.604 groups per square nm (surface 9). View along x axis: top left; view along y axis: top right; view along
18 z axis: bottom left; isometric view: bottom right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface is illustrated with square in light yellow color. Water molecules and pentyl groups are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms...... 159
Figure 85. Simulation snapshot of the quartz surface with pentyl concentration 2.893 groups per square nm (surface 10). View along x axis: top left; view along y axis: top right; view along z axis: bottom left; isometric view: bottom right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface is illustrated with square in light yellow color. Water molecules and pentyl groups are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms...... 160
Figure 86. Simulation snapshot of the quartz surface with pentyl concentration 3.182 groups per square nm (surface 11). View along x axis: top left; view along y axis: top right; view along z axis: bottom left; isometric view: bottom right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface is illustrated with square in light yellow color. Water molecules and pentyl groups are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms, light blue balls - carbon atoms...... 161
Figure 87. Side views of a water droplet on the quartz surface with removed pentyl groups for pentyl concentration 0.289 groups per square nm (surface 1). View along x axis: left; view along y axis: right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface position is indicated by the lowest laying water molecules. Water molecules are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms. .... 162
Figure 88. Side views of a water droplet on the quartz surface with removed pentyl groups for pentyl concentration 0.579 groups per square nm (surface 2). View along x axis: left; view along y axis: right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface position is indicated by the lowest laying water molecules. Water molecules are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms. .... 162
Figure 89. Side views of a water droplet on the quartz surface with removed pentyl groups for pentyl concentration 0.868 groups per square nm (surface 3). View along x axis: left; view along y axis: right. The sphere fitted to the iso-density chart is illustrated in light purple color.
19 The surface position is indicated by the lowest laying water molecules. Water molecules are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms. .... 163
Figure 90. Side views of a water droplet on the quartz surface with removed pentyl groups for pentyl concentration 1.157 groups per square nm (surface 4). View along x axis: left; view along y axis: right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface position is indicated by the lowest laying water molecules. Water molecules are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms. .... 163
Figure 91. Side views of a water droplet on the quartz surface with removed pentyl groups for pentyl concentration 1.447 groups per square nm (surface 5). View along x axis: left; view along y axis: right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface position is indicated by the lowest laying water molecules. Water molecules are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms. .... 164
Figure 92. Side views of a water droplet on the quartz surface with removed pentyl groups for pentyl concentration 1.736 groups per square nm (surface 6). View along x axis: left; view along y axis: right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface position is indicated by the lowest laying water molecules. Water molecules are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms. .... 164
Figure 93. Side views of a water droplet on the quartz surface with removed pentyl groups for pentyl concentration 2.025 groups per square nm (surface 7). View along x axis: left; view along y axis: right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface position is indicated by the lowest laying water molecules. Water molecules are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms. .... 165
Figure 94. Side views of a water droplet on the quartz surface with removed pentyl groups for pentyl concentration 2.314 groups per square nm (surface 8). View along x axis: left; view along y axis: right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface position is indicated by the lowest laying water molecules. Water molecules are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms. .... 165
Figure 95. Side views of a water droplet on the quartz surface with removed pentyl groups for pentyl concentration 2.604 groups per square nm (surface 9). View along x axis: left; view along y axis: right. The sphere fitted to the iso-density chart is illustrated in light purple color.
20 The surface position is indicated by the lowest laying water molecules. Water molecules are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms. .... 166
Figure 96. Side views of a water droplet on the quartz surface with removed pentyl groups for pentyl concentration 2.893 groups per square nm (surface 10). View along x axis: left; view along y axis: right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface position is indicated by the lowest laying water molecules. Water molecules are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms. .... 166
Figure 97. Side views of a water droplet on the quartz surface with removed pentyl groups for pentyl concentration 3.182 groups per square nm (surface 11). View along x axis: left; view along y axis: right. The sphere fitted to the iso-density chart is illustrated in light purple color. The surface position is indicated by the lowest laying water molecules. Water molecules are shown in VDW representation, white balls - hydrogen atoms, red balls - oxygen atoms. .... 166
21 1. Introduction
1.1. Scope and objective Carbon geo-sequestration is a viable mitigation strategy to tackle consequences of excessive anthropogenic emissions of carbon dioxide (IPCC, 2005). The strategy, among other options, relies on depleted oil or gas reservoirs as possible storage formations. After being exposed to hydrocarbons, cap rocks of these formations may undergo chemical transformations, which change their wettability (Ali et al., 2019) and thus capacity, reliability and security of the structural and residual trapping (Stefan Iglauer, Pentland, & Busch, 2014). The aim of this work is to understand influence of these chemical transformations on the wettability of quartz surfaces in presence of carbon dioxide and to quantify this influence.
1.2. Methods, procedures and process To study the effects of chemisorpt hydrocarbons on the wettability of the hydroxylated quartz a range of methods of computational chemistry, materials science and solid state physics is used. Known and reliable force fields (Cygan, Liang, & Kalinichev, 2004; Mayo, Olafson, & Goddard, 1990; van Beest, Kramer, & van Santen, 1990) are adapted to the new chemical environments using the density functional theory (Hohenberg & Kohn, 1964; Kohn & Sham, 1965). Investigation on how the surface density of alkyl groups affects quartz's wettability in carbon dioxide atmosphere is then performed using classical molecular dynamics (Gonzalez, 2011). A possibility to use spheroidal geometric constructions to calculate the contact angle using individual snapshots of molecular dynamics simulations is studied with the most stable surface of the pristine quartz.
1.3. Results, observations and conclusions Insights gained at the molecular level indicate that the wettability of alkylated quartz is a complex phenomenon complete characterisation of which requires a number of parameters. Reported contact angles obtained at different surface levels are compared with the contact angle that is most likely to be measured in an experiment. Quantitative confirmation of the hydrophobicity of the hydroxylated and alkylated quartz is obtained.
1.4. Novel and additive information This research contributes to the scientific and engineering community in two ways. Firstly, there is a fundamental contribution in the form of modified and adapted for the modelling of
22 the alkylated quartz force fields. Secondly, there is a practical contribution in the form of a deeper understanding of the wetting processes of chemically modified with hydrocarbons quartz surfaces. Reported quantification of the wettability of the latter surfaces with respect to the concentration of alkyl groups is of particular importance for the engineering applications and development of effective carbon geo-sequestration projects.
1.5. Thesis organization The dissertation is organized as follows. In chapter two author presents an overview of information in existing literature on the subject of climate change, carbon capture and storage technology, parameters affecting applicability and efficiency of the technology, specifically the contact angle in quartz-water-carbon dioxide systems. The chapter ends with a summary of identified knowledge gaps. In chapter three an overview of existing computational methods utilized to eliminate the knowledge gaps is given. The following chapters are organized as self- sufficient and potentially publishable summaries of research results obtained as independent and yet connected to each other parts of this project. In chapter four a possibility of using spheroidal geometric constructions for calculation of the contact angle is explored. The most stable pristine quartz surface is used as a proving ground for proposed in the chapter method. In chapter five an oxygen-carbon bond connecting quartz crystal to alkyl rest is modelled using the density functional theory. Charges on both atoms of the bond are estimated using the Bader charge analysis. In chapter six a systematic study of wettability of hydroxylated and alkylated quartz with respect to the surface pentyl density is performed. From chapter four through chapter six, the narrative thus naturally proceeds from development of a method to calculate the contact angle of a spherical droplet (chapter four), to a solution of encountered scientific problem of adapting known force fields to chemical environments of studied alkylated quartz surfaces (chapter five), and finally to application of modified contact angle calculation method and modified force fields to systematic study of wettability of alkylated quartz (chapter six). The work is concluded in chapter seven.
1.6. Publications A. Abramov, S. Iglauer, "Analysis of individual molecular dynamics snapshots simulating wetting of surfaces using spheroidal geometric constructions", accepted for publication in The Journal of Chemical Physics.
23 A. Abramov, S. Iglauer, "Application of the CLAYFF and the DREIDING force fields for modelling of alkylated quartz surfaces", Langmuir, DOI: 10.1021/acs.langmuir.9b00527, URL: https://pubs.acs.org/doi/10.1021/acs.langmuir.9b00527.
A. Abramov, A. Keshavarz, S. Iglauer, "Wettability of fully hydroxylated and alkylated (001) alpha-quartz surface in carbon dioxide atmosphere", Journal of Physical Chemistry C, DOI: 10.1021/acs.jpcc.9b00263, URL: https://pubs.acs.org/doi/10.1021/acs.jpcc.9b00263.
24 2. Literature review
2.1. Global warming
2.1.1. Overview Operating for over than 30 years (since 1988) an intergovernmental scientific body - The Intergovernmental Panel on Climate Change (IPCC) unambiguously links human industrial activity (Figure 1) and changes in the climate and ecosystems observed globally (Figure 2, Figure 3). The IPCC uses recorded from 1850 instrumental measurements and data derived from paleoclimate archives to set a long-term context (IPCC, 2013). As Figure 1 shows industrial activity of the civilization manifests itself in exponentially growing emissions of carbon dioxide, which reach 10 GtC per annum. Action of natural CO2 sinks significantly mitigated the impact, but from 1800 measured atmospheric growth rate of CO2 have been steadily increasing and reached more than 3 GtC per annum, Figure 1. These anthropogenic emissions result in 25% increase in CO2 concentration over a period of 50 years, from 1960 to 2010, Figure 2. Increased carbon dioxide concentrations are consistently observed in the Northern and Southern Hemispheres. Higher content of the acid gas in the atmosphere also leads to increased ocean acidity; reduced pH levels have been recorded for several decades in a row (IPCC, 2013), Figure 2.
Numerous indicators demonstrate that Earth's climate is changing (warming), and from 1950 trends of the change are steady and the rate is alarming, Figure 3: land surface temperature increased by 1C, sea surface temperature increased by 0.3C, marine air temperature increased by 0.4C, sea level increase by more than 100 mm, summer arctic sea-ice extent decreased by 5106 km2. Increasing trends are also demonstrated by tropospheric temperature, ocean heat content and specific humidity; decreasing trends are observed for Northern Hemisphere snow cover and glacier mass balance, see Figure 3. Thus, collected and analysed by the intergovernmental body paleoclimate data shows clear connection between
CO2 concentration and global temperature. With more than 40% CO2 concentration increase since preindustrial times, atmospheric CO2 content became higher than it has ever been over the past 800 000 years (more than 410 ppm in 2018 (scripps.ucsd.edu)). Every of the past three decades had been warmer than the previous decades since 1950, with the first decade of the twenty first century been the warmest of all (IPCC, 2013).
25
Figure 1. Emissions of carbon dioxide due to industrial and agricultural activities (1 Petagram of carbon = 1 PgC = 1015 grams of carbon = 1 Gigatonne of carbon = 1 GtC). Image source: (IPCC, 2013).
26
Figure 2. Atmospheric CO2 concentration, CO2 partial pressure (pCO2, left axis) and pH of oceanic surface. MLO: Mauna Loa Observatory, Hawaii; SPO: South Pole; HOT: Hawaii Ocean Time-Series station. Image source: (IPCC, 2013).
This warming is explained by increased heat uptake from the sun by the climate system. The difference between absorbed and reflected energy, i.e. the radiative forcing, has been increasing with increase in CO2 concentration due to anthropogenic industrial activity. In fact, the increase in the atmospheric concentration of CO2 since 1750 is the largest contributor to the total radiative forcing. As a result of warming caused by anthropogenic radiative forcing ice sheets and glaciers are losing mass. Arctic sea ice area has been declining since 1980, as well as the Greenland ice sheet. With smaller snow and ice area more of sun's heat is absorbed by the land and oceans, which creates an alarming positive feedback influencing the climate system.
Existing climate system trends, see Figure 3, led the IPCC to the following conclusions (IPCC, 2013): (i) warming in the climate system is undeniable; (ii) human influence on the climate system is clear; (iii) continued greenhouse gases emission causes further climate change and long term negative consequences. Limiting the climate change requires significant and permanent reduction of greenhouse gases emission.
27
Figure 3. Climate change indicators. Image source: (IPCC, 2013).
Greenhouse gases emission causes the primary effect of the climate change - the warming. There are number of secondary effects among which are the following, see Figure 2 and Figure 3: rise of sea levels, increase in oceans acidity, change of weather and precipitation patterns, spread of tropical diseases, to list some. Mentioned global processes are intrinsically negative, although some local effects may be considered as to say "positive", such as appearance of new sea routes and accessibility of otherwise ice covered areas, warmer winters and reduced energy consumption for heating. The net worldwide outcome of these developments is believed to be extremely negative and catastrophic.
In October 2018 the IPCC issued a special report "Global warming of 1.5C" (IPCC, 2018), where the authors argue that increase in global temperature above 1.5C, which is estimated
28 to happen around 2030 if current trends of CO2 emission are preserved, could render the planet unsuitable for living in long term due to irreversible changes, such as the loss of ecosystems. A Greenpeace representative summarized the report in three words written in capital letters: "ACT NOW, IDIOTS" (www.bbc.com).
2.1.2. Response strategies Apart from a passive strategy of doing nothing regarding the climate change it is possible to single out several active strategies which can be used on their own or in combination (IPCC, 2014): (i) adaptation; (ii) mitigation and (iii) climate engineering.
As the word "adaptation" suggests the corresponding strategy implies changing the living environment so it better fits warmer conditions. These changes may include usage of heat resistant and less water demanding plants, flood protective structures, higher foundations and street levels, light coloured roofs, weather control techniques, to name a few.
The mitigation strategies can be developed on premises that emission of greenhouse gases originates from energy production satisfying certain energy consumption, improvements can be implemented on both sides. Limiting energy consumption measures include altering behaviour, lifestyle and diet of the population (mode and demand for mobility, households' energy use, choice of longer-lasting products, reduction in food wastes), as well as energy conservation and energy efficiency. According to the IPCC these are the key mitigation strategies on the energy demand side to limit CO2 concentration to about 450 ppm by 2100 (IPCC, 2014). On the energy supply side reducing the carbon intensity (decarbonizing) of electricity generation is crucial to achieve cost-effective mitigation of negative consequences of the climate change. This can be accomplished by switching the energy production to renewable sources and specifically to nuclear energy, which share of world electricity generation has been in decline since 1993 (IPCC, 2014). If such switch in energy production is not feasible in certain locations, then large CO2 emitters can be equipped with carbon capture and storage facilities which remove carbon dioxide from the exhaust fumes. In addition to these measures natural or artificial reservoirs absorbing carbon dioxide, or carbon sinks, can be preserved, enhanced or created. Oceans (at the expense of acidification) and forests are among the largest carbon sinks.
29 Special case of a mitigation strategy is the climate engineering or geoengineering. The strategy is largely represented by two subcategories: solar radiation management and removal of CO2 from the atmosphere. The former method aims to reflect part of incoming sun light back to space. The letter method is analogous to the carbon capture and storage although it deals with very low CO2 concentrations due to which it is much less efficient.
2.2. Carbon capture and storage
2.2.1. Overview
Carbon dioxide Capture and Storage (CCS) technology can reduce CO2 emissions of fossil fuel power plants, where the gas is captured, transported to a place of storage and isolated from the atmosphere deep underground (Bachu, 2003; GCCSI, 2018; IPCC, 2005). The technology is appealing as it can be used on already operating power stations, as well as on prospective ones. The capturing is the most efficient when it is done at the emission source where three different types of capture systems may be applied: (i) post-combustion; (ii) pre-combustion and (iii) oxyfuel combustion (IPCC, 2005). In the first system CO2 removed after combustion of organic or fossil fuels. In the second system carbon based fuel is partially oxidised to form hydrogen and carbon monoxide, which then reacts with water to form CO2 and hydrogen. CO2 is next removed to leave hydrogen for further processing or combustion. In the oxyfuel (or oxy-firing) combustion fuel is burnt in pure oxygen which produces only CO2 and water. The first two systems are used by the industry for purposes different from the CCS, the third system is in demonstration stage.
To separate CO2 from fume gases absorption, adsorption or membrane technologies can be applied. Out of three methods absorption via amine scrubbing is the most mature and widely used by the industry technology, for example, in natural gas processing. After separation CO2 is transported in supercritical state using existing or specially constructed pipeline networks to a place of storage where it is then injected into a geological formation at depth below 800 m (IPCC, 2005). This formation may be an oil reservoir where CO2 is used for the enhanced oil recovery or it may be specially dedicated for CO2 storage formation, see Figure 4.
30 CO2 Injection well Injection 800 m or or more m 800
Geological storage formation
Figure 4. Schematic illustration of carbon dioxide capture and storage technology. Geological storage formations may vary and can be selected from specially dedicated for storage formations, depleted oil or gas reservoirs, deep aquifers, coal beds or salt caverns.
The storage formation must be determined carefully (Bachu, 2003) so that acting trapping mechanisms prevent CO2 from leaking back to the surface because of its buoyancy. Four trapping mechanisms have been identified and are being actively studied: (i) structural trapping; (ii) residual trapping; (iii) solubility (or dissolution) trapping and (iv) mineral trapping (Stefan Iglauer et al., 2014).
The structural trapping is well understood mechanism preventing upward hydrocarbon migration. Over geological timescale oil, gas or both accumulate and separate from water beneath impenetrable cap rocks. As with hydrocarbons these natural seals stop low density supercritical CO2 from moving toward the surface retaining it underground. It should be noted that actual seal configuration does not have to be dome like, overlaying along geological faults seal formations can play similar role.
31 In the residual trapping micrometre sized bubbles of CO2 disconnected from the main body of the gas are immobilised in the pore space of the storage formation. CO2 mobility is halted by the capillary forces which balance or exceed the buoyancy forces.
In the solubility or dissolution trapping some portion of CO2 dissolves into water present in the pore space of surrounding rocks. This causes increase in water density and its corresponding downward motion.
The mineral trapping refers to formation of minerals which are produced as a result of chemical reactions between CO2, water, ions present in water and rocks of the storage reservoir. Reacted CO2 becomes permanently captured in the solid phase.
Among the major difficulties reducing confidence in the carbon capture and storage technology is limited ability to robustly predict storage security of the first two trapping processes.
2.2.2. Carbon geo-sequestration projects The Global CCS Institute (GCCSI) headquartered in Melbourne, Australia
(www.globalccsinstitute.com), in its CO2RE database identifies 23 large-scale CCS projects already operating or under construction, with additional 28 pilot and demonstration projects conceived or functioning (GCCSI, 2018). These facilities capture 43 Million Tonnes Per Annum
(MTPA) of CO2.
The current level of development in CCS technology was reached in more than 45 years. Enhanced oil recovery (EOR) pioneered practical implementations of CCS in 1972 with Val Verde EOR facility in the Sharon Ridge oilfield in Texas. Governmental regulations introduced in Norway in 1991 motivated Statoil to start the Sleipner CO2 storage facility in 1996, which became the first in the world project where CO2 was injected into a dedicated geological formation. Injection of CO2 on the Snohvit gas field, which started in 2008, was required by a license condition also imposed by the Norwegian State. In Weyburn and Midale oil fields in
Canada CO2 have been injected for EOR from 2000 and 2005, respectively. The In Salah CO2 storage facility in Algeria injected carbon dioxide into the Krechba formation, a depleted gas reservoir, from 2004 to 2010. Over the next five years geochemical, geophysical and production techniques were applied to observe the CO2 storage process. Obtained data was then used to improve modern monitoring, modelling and verification methodologies.
32 The second decade of the twenty-first century was the most fruitful in terms of CCS applications in diverse industrial areas and reached sequestration capacity. The technology is now in use in coal-fired power plants, iron and steel factories, chemical and hydrogen production plants.
The first large-scale carbon capture facility in power generation was commenced in 2014 with reconstructed Unit 3 of Boundary Dam Power Station in Saskatchewan province in western Canada (1 MTPA). The first commercial application of CCS in the oil sands industry was also achieved in Canada by Shell in cooperation with the Canadian Federal Government and the Provincial Government of Alberta in 2015 in the Quest facility (1 MTPA).
After 2009 six large-scale CO2 sequestration facilities began operation in the US. Launched in 2017 the Illinois industrial CCS facility uses dedicated storage (1 MTPA), other five capture and storage units operate in EOR mode: Century natural gas processing plant (8.4 MTPA), Air Products hydrogen production steam methane reformer (1 MTPA), Coffeyville gasification plant producing fertilisers (1 MTPA), Lost Cabin natural gas processing plant (0.9 MTPA), Petra Nova coal-fired power plant with carbon capture (1.4 MTPA).
In Brazil Petrobras started operating CCS facilities in Santos Basin in 2015. In 2017 ten floating production storage and offloading units at the Lula, Sapinhoa and Lapa fields were injecting
CO2 for EOR (2.5 MTPA).
In Saudi Arabia the Uthmaniyah facility was functioning since 2015. CO2 from a natural gas processing plant is used for EOR on the largest in the world Ghawar field (0.8 MTPA).
The first large-scale CCS facility in the iron and steel industry started operation in 2016 in the United Arab Emirates. The Emirates Steel Industries factory in Mussafah captures and transports CO2 to Abu Dhabi National Oil Company for EOR (0.8 MTPA).
The EOR project in Jilin oilfield in China was in development since 1990. This first large-scale CCS facility operating in China entered Phase III in 2018 (0.6 MTPA).
The Gorgon liquefied natural gas project in Australia involving the development of the Greater
Gorgon gas fields, in some of which the natural gas contains high concentrations of CO2, is expected to capture and reinject up to 4 MTPA of CO2 upon completion.
33 2.3. Wettability and carbon dioxide storage capacity of rocks Reliability and long term security of the structural and the residual trapping mechanisms are strongly influenced by the wettability of storage formations in CO2 atmosphere and related capillary effects. The wettability or wetting may intuitively be understood as an ability of a liquid to spread over a solid surface. The degree of spreading is determined by a balance between adhesive and cohesive forces. The macroscopic phenomenon of wetting essentially originates from intermolecular interactions occurring when the liquid and the solid are brought in contact with each other. It is implied that apart from contacting liquid and solid a third party is in play - the gas in which wetting is observed. In such a setup wettability of the solid may be understood as its preference to be in contact with or to be wet by either of two fluids the liquid or the gas.
Wettability can quantitatively be characterised by the contact angle given by the Young equation, which expresses thermodynamic equilibrium between three phases - the solid phase, the liquid phase and the gas phase (Atkins & De Paula, 2006):
훾 − 훾 푐표푠휃 = , 훾 where is the contact angle, is the interfacial tension and S, L and G stand for the solid, the liquid and the gas, respectively, see Figure 5 and Figure 6.
LG
Gas Liquid
Solid SL SG
Figure 5. A liquid droplet on a solid surface and variables used in the Young equation.
The contact angle can take values between 0 and 180 degrees. The value 90 degrees separates wetting (less than 90 degrees) and non-wetting (more than 90 degrees) regimes. Complete
34 surface wetting is identified by the zero contact angle and completely non-wetting state is characterised by the contact angle 180 degrees.
The contact angle is related to the capillary pressure, which is the pressure difference across the interface which separates two fluids (Stefan Iglauer et al., 2014):
푃 =∆푃 = 푃 − 푃 , where Pnw and Pw are pressures of the non-wetting (carbon dioxide) and wetting (water) phases.
If a formation pore is modelled as a capillary tube, then it can be shown that the contact angle is connected to the capillary pressure through the Young-Laplace equation (Dake, 1978):
∆푃 =2훾퐶, where is the interfacial tension of the fluid pair, C is the mean curvature of the interface.
SG
rc
LG SL
Figure 6. Capillary rise of a liquid illustrating parameters used in the Young-Laplace equation.
35 For a capillary tube shown in Figure 6:
푃 =2훾/푅,
푅 = 푟 /푐표푠휃,
푃 =2훾푐표푠휃/푟 , where R is the principle radius of curvature, and rc is the radius of the capillary tube (or pore radius).
Pc
Cap rock Pb h
CO2+residual H2O
H2O+residual CO2
Figure 7. Schematic illustration of CO2 structural trapping.
Capillary pressure which prevents non-wetting phase (carbon dioxide in the context of carbon capture and storage) from entering the pores is counterbalanced by the buoyancy of the non- wetting phase:
푃 =∆휌푔ℎ, where is the density difference between CO2 and water, g is acceleration due to gravity, h is the height of the CO2 plume, see Figure 7.
Equating Pc and Pb an estimation of the CO2 storage capacity is obtained (Stefan Iglauer et al., 2014):
2훾푐표푠휃 ℎ= . ∆휌푔푟
There are four variables in obtained equation: water-CO2 interfacial tension, the contact angle, the density difference between CO2 and water, and the radius of the capillary tube (the
36 pore radius). The latter one is the most difficult to estimate due to complex pore morphology of real rocks. Because of the same reason it is also less representative of real pore environments. On the other hand, the first two variables have extensively been studied in experimental and theoretical works. The third variable, the density difference, predominantly depends on properties of CO2 which in turn strongly depend on pressure and temperature conditions.
2.4. Properties of carbon dioxide in geological conditions
At standard conditions of pressure and temperature CO2 is a stable gas with the density 1.976 kg/m3 (IPCC, 2005). Physical state of carbon dioxide varies with temperature and pressure according to its phase diagram, Figure 8 (www.chemicalogic.com). Three phases solid, liquid and gas coexist at -56.5C and 5.18 bar. CO2 becomes supercritical at 31.1C and
73.9 bar. It is important to note that at supercritical conditions CO2 fills available volume as a gas, but can have densities comparable with densities of its liquid state from 150 to more than 800 kg/m3, depending on pressure and temperature, see Figure 9 (Bachu, 2003; IPCC, 2005).
Figure 8. Phase diagram of carbon dioxide. Image source: ChemicaLogic Corporation, Copyright (c) 1999 ChemicaLogic Corporation, USA. All rights reserved.
37 Density of carbon dioxide is crucial for the storage efficiency and safety. Larger densities increase amount of stored CO2 in terms of its mass and decrease its upward acting buoyancy force. As Figure 9 illustrates CO2 density strongly linked to pressure and temperature and thus to depth at which it is stored. The greater the depth the greater the pressure and the greater the temperature, but their effect on CO2 density is opposite. In the beginning CO2 density increases rapidly with pressure, then temperature counteracts pressure which stabilises the density. Considering interplay of pressure and temperature with depth it was show that there exist an optimal depth range from 800 to 1000 m for "cold" geological basins and from 1500 to 2000 m for "warm" geological basins, where CO2 density reaches plateau and almost ceases increasing or even decreases (Bachu, 2003). Expected storage pressure and temperature conditions and physical state of CO2 at optimal depths can be determined from the hydrostatic and the geothermal gradients.
Variations of pressure P and temperature T with depth are given by formulae:
푃 = 휌 푔푧,
푇 = 푇 + 퐺푧, where ρw is water density, g is the acceleration due to gravity, z is the depth, Ts is the surface temperature, G is the geothermal gradient.
Eliminating the depth, the pressure and the temperature can be combined in one equation and superimposed onto CO2 phase diagram (Bachu, 2003):
휌 푔 푃 = (푇 − 푇 ). 퐺
38
Figure 9. CO2 density as a function of temperature and pressure. Image source: (Bachu, 2003; IPCC, 2005).
Overlay of geological conditions and physical states of CO2 in coordinates pressure versus temperature is depicted in Figure 10. The "cold" and the "warm" basin conditions shown in Figure 10 are defined by the surface temperatures -2 and 30C, respectively, and by the geothermal gradients 20 and 60C/km, respectively (Bachu, 2003). The CO2 saturation line is constructed with the Span and Wagner equation of state (Span & Wagner, 1996) using a web based tool (www.energy.psu.edu). The diagram in Figure 10 illustrates that if stored at the optimal depths, CO2 could be either in supercritical or liquid states with pressures ranging from 7.85 to 19.62 MPa and temperatures covering interval from 14 to 150C. In "cold" basins
3 as depth increases from 800 to 1000 m CO2 density decreases from 874 kg/m (at 14C and 7.85 MPa) to 868.2 kg/m3 (at 18C and 9.81 MPa), respectively. In "warm" basins as depth
3 increases from 1500 to 2000 m CO2 density increases from 273.5 kg/m (at 120C and 14.7 MPa) to 320 kg/m3 (at 150C and 19.62 MPa), respectively.
Figure 10 also depicts one (out of many) possible "critical path", the path in the pressure- temperature space which crosses the critical point when CO2 moves upwards from the storage formation. Apart from this unique route characteristic phase transformations in warmer
39 conditions are from the supercritical state to the gas state and in the colder conditions from the liquid state to the gas state (with uptake of the latent heat). In colder conditions more general metamorphoses from the supercritical to the liquid and then to the gas states are possible as well.
20 18 16 14 12 10 8 Supercritical fluid Liquid Pressure, MPa 6 Gas 4 2 0 0 20 40 60 80 100 120 140 160 Temperature, C Saturation line Critical point Supercritical region "Warm" basin "Cold" basin Critical path
Figure 10. Possible geological pressure and temperature conditions in "cold" and "warm" basins worldwide superimposed onto carbon dioxide phase diagram. Thickened lines highlight conditions at optimal depths 800- 1000 m for "cold" and 1500-2000 m for "warm" basins, respectively. Interpolation between optimal "cold" and "warm" conditions (gradient coloured lines) shows range of pressures and temperatures expected at CO2 storage depths. Image is constructed on basis of (Bachu, 2003). 2.5. Quartz wettability in presence of carbon dioxide Quartz is the second most abundant mineral in the Earth's continental crust (about 12% by weight) (Anderson & Anderson, 2010; Marshall & Fairbridge, 1999), it is behind feldspars which compose about half of the crust. In addition to its global natural abundance quartz is present in large concentrations in cap rocks, more than 50% by weight (Stefan Iglauer, Al- Yaseri, Rezaee, & Lebedev, 2015). Quartz wettability is therefore crucial for characterisation of overall wettability of real multi-mineral cap rocks. To understand and quantify wettability of quartz surface it has extensively been studied not only in numerous experiments but also with advanced methods of theoretical chemistry and materials science.
40 2.5.1. Experimental studies Experimental wettability studies explored different gaseous atmospheres including air (Mazurek, Pogorzelski, & Boniewicz-Szmyt, 2009; Roshan, Al-Yaseri, Sarmadivaleh, & Iglauer, 2016), water vapour saturated air (Lamb & Furlong, 1982), hydrogen sulphide and carbon dioxide (Broseta, Tonnet, & Shah, 2012; Espinoza & Santamarina, 2010). Variety of water solutions have been tested ranging from distilled water (Sklodowska & Matlakowska, 1997) to brines of different salinity (Ahmed Z. Al-Yaseri, Lebedev, Barifcani, & Iglauer, 2016; Ahmed Zarzor Al-Yaseri, Roshan, Lebedev, Barifcani, & Iglauer, 2016; Broseta et al., 2012). Diversity of combinations surface-liquid-gas examined in the experiments is enriched by multiplicity of contact angle measurement methods. The are many approaches to measure the contact angle experimentally which are described in specialized textbooks, for example, Chapter 1 "Contact Angle and Wetting Properties" in (Bracco & Holst, 2013) offers nine methodologies for assessing the contact angle experimentally. Some common techniques to measure the contact angle relatively widely seen in research literature are: the sessile drop method (Sarmadivaleh, Al-Yaseri, & Iglauer, 2015), the captive bubble method (Kaveh, Rudolph, van Hemert, Rossen, & Wolf, 2014; Saraji, Goual, Piri, & Plancher, 2013), the capillary rise method (Koh, Hao, Smith, Chau, & Bruckard, 2009) and the direct measurement with the atomic force microscopy (Deng, Xu, Lu, Wang, & Shi, 2018).
In the sessile drop method, a water droplet is placed on the surface of a mineral in atmosphere of the less dense phase, i.e. air, CO2 or other gases. The captive bubble method is an analogous approach where the less dense phase is injected underneath the surface of a mineral submerged into denser phase, i.e. water. In the capillary rise method, the contact angle is calculated using the values of experimentally measured capillary rise and radius (Bracco & Holst, 2013). In contrast to other methods, where droplet sizes are in the range of millimetres, the atomic force microscopy is used to measure the contact angle of micrometre-sized droplets. In the technique, the force between the sharp tip and a surface as a function of their separation is detected and used to reproduce the shape of the sample. Very high resolution of less than nanometre is attainable with the atomic force microscopy (Bracco & Holst, 2013).
At the micrometre level of the atomic force microscopy, size of the water droplet (or bubble, in a more general case) has negligible effect on the contact angle measured in air atmosphere (Deng et al., 2018). For larger droplets, however, size matters and affects outcome of
41 measurements. The gravitational influence becomes significant when radius of the droplet exceeds the threshold value given by the capillary length. This threshold condition is fulfilled when the gravitational forces become comparable with the surface tension forces. Relative importance of one type of forces over the other is formally expressed by a dimensionless number, the Bond number (also known as the Eötvös number) (Hager, 2012):
∆휌푔푟 퐵 = , 훾 where is the density difference between the liquid and vapour phases, g is acceleration due to gravity, r is the droplet radius, is the interfacial tension of the fluid pair.
The droplet size therefore becomes an important parameter affecting measurements of the contact angle when:
퐵 ≥ 1.
Experimental assessments of the contact angle for all droplet sizes are influenced by the surface roughness and heterogeneity (Drelich, Miller, & Good, 1996). It was demonstrated in (Deng et al., 2018) that even for a micrometre sized water droplet the contact angles measured along the triple-phase contact line may have almost twofold difference, ranging from 27.8 to 48.9. To quantify effects of surface's roughness and heterogeneity on the contact angle and to link contact angles obtained for ideal smooth and homogeneous surfaces with real contact angles on rough and heterogeneous surfaces two models have been developed (Cassie & Baxter, 1944; Wenzel, 1936). Effect of surface roughness on the contact angle is expressed by an equation proposed by Robert Wenzel (Wenzel, 1936):
푐표푠휃 = 푅 푐표푠휃, where W and are the Wenzel contact angle and the contact angle obtained on the smooth and homogeneous surface, respectively; Rf is the roughness factor which is defined as the ratio of the area of the rough surface to the area of projected flat surface (in practice the roughness factor is always larger than one).
Wenzel's theory assumes that the roughness grooves are completely filled with liquid for which the contact angle is measured, the assumption which is not required for the Cassie- Baxter model (Cassie & Baxter, 1944):
42 푐표푠휃 = 푓 푐표푠휃 + 푓 푐표푠휃 , where CB is the Cassie-Baxter contact angle; f1 and f2 are the area fractions for surface components 1 and 2, for which the contact angles are 1 and 2, respectively (in principle, the model can be extended to any number of surface components).
Gas Gas
Liquid Liquid
Surface Surface
Figure 11. Liquid droplet on a rough and inhomogeneous surface. The Wenzel (left) and the Cassie-Baxter (right) models.
The Cassie-Baxter model can handle situations where the roughness grooves are filled with a gas (Figure 11), which becomes one of the surface components, let us say the second one. In a water-air system this means that 2=180, and with f2=1-f1 one arrives at:
푐표푠휃 = 푓 (푐표푠휃 +1) − 1.
It is interesting to note that for a given surface material the contact angle can be tuned if necessary just by altering f1. In the context of CO2 storage, just because of surface roughness and the Cassie-Baxter wetting regime, the contact angles of water may be heightened and CO2 storage capacity and security may be reduced.
As may be inferred from presented in this section material the contact angle is a sensitive to number of parameters characteristic. Among those parameters are surface chemistry (hydrophilicity and hydrophobicity (Stefan Iglauer, 2017)) and surface conditions (roughness and homogeneity, pre-treatment and crystal face (Deng et al., 2018)), size of the droplet (or bubble), method used to measure the contact angle (Drelich et al., 1996). In addition, in the context of carbon geo-sequestration, attention should be paid to salinity of the denser phase and ionic strength of the solution, experimental atmosphere including type of gas and its
43 pressure and temperature. Increase in salinity, valency of dissolved ions, and in carbon dioxide pressure causes increase in CO2-wettability, while effect of temperature is less clear (Stefan Iglauer, 2017). Understanding of influence of these latter four factors (pressure and temperature, salinity, type of ions) on water/CO2 wetting behaviour of quartz is crucial for correct estimation of CO2 storage capacity and security of existing and prospective carbon geo-sequestration projects.
An attempt (Roshan et al., 2016) to understand effects of pressure and temperature on the contact angle was made using an analytical expression derived in (Garcia, Osborne, & Subashi, 2008, 2009):
∆휌 푐표푠휃 =−1− 푉 (푧)푑푧, 훾 where is the density difference between the liquid and vapour phases, is the interfacial tension of the fluid pair, Vs is the potential energy of the adsorbate molecule due to the substrate (the substrate potential), zmin is the position of the minimum of the potential near the surface.
Because of high CO2 compressibility there is one strongly dominating factor affecting the contact angle with varying pressure - carbon dioxide density, which upon increase enhances
CO2-mineral interactions (Stefan Iglauer, 2017). In terms of presented equation, the effect of pressure may be understood as follows. With increase in pressure the density difference between the liquid (water) and vapour (CO2) phases decreases, see Figure 9. It was shown that the ratio of the van der Waals potential integral over the interfacial tension for given temperature, brine composition and mineral substrate is constant (Ahmed Zarzor Al-Yaseri et al., 2016). Reduced due to increased pressure density difference in the second term on the right hand side of the equation is thus the main factor responsible for increase of the contact angle.
For the sake of completeness effect of pressure on the interfacial tension should be characterised too. Pressure increase from ambient conditions causes rapid decrease in water-
CO2 interfacial tension from about 70 to 30 mN/m at temperatures ranging from 278 to 333 K (Hebach et al., 2002), after pressure reaches the saturation or supercritical values the interfacial tension remains almost unchangeable (Hebach et al., 2002). It is necessary to note
44 that in relative terms the effect of pressure on water-CO2 interfacial tension is much smaller than its effect on CO2 density, Figure 9.
With variations in temperature there is no sole dominating factor, the density difference and the interfacial tension of the fluid pair are both affected; as well as the kinetic energy of the adsorbate molecules, which in turn results in variations in explored by the adsorbate molecules energies of the potential well. That is why the effect of temperature on the contact angle is more intricate. For example, the contact angle decreases with increase in temperature for coal and increases with increase in temperature for quartz (Stefan Iglauer, 2017). Qualitative picture of the advancing contact angle dependence on pressure and temperature of a deionized water droplet on quartz surface cleaned with air plasma in carbon dioxide atmosphere is shown in Figure 12.
Figure 12. Advancing contact angle of a deionized water droplet on rough quartz surface (RMS 560 nm) cleaned with air plasma in carbon dioxide atmosphere at different conditions of pressure (in MPa) and temperature (in Kelvin). Source of data: (Ahmed Z. Al-Yaseri et al., 2016).
It is clear from Figure 12 that increase in temperature and pressure has the same qualitative effect on the contact angle of a water droplet on quartz surface in CO2 atmosphere, both parameters increase the contact angle. It is also clear that at the geological sequestration conditions the contact angle is high, up to 50, and the quartz surface is only weakly water- wet (Ahmed Z. Al-Yaseri et al., 2016).
Experimentally measured advancing contact angles a should be expected to be a good approximation of the static contact angles defined by the Young's equation (Bracco & Holst,
45 2013). In general, the wetting phenomena is dynamic, the contact angles formed by expanding liquid are called "advancing" - a, and the contact angles formed by contracting liquid are called "receding" - r. The difference between the advancing and the receding angles originates from surface roughness and heterogeneity and is referred to as the contact angle hysteresis:
퐻=휃 −휃 .
After the PT-conditions, water salinity and type of dissolved ions are known to significantly affect quartz wettability in CO2 atmosphere, Figure 13. The effect of salinity is explained by screening of surface charges by salt ions of opposite charge which from double layer reducing surface electrostatic potential and thus reducing overall surface polarity and hydrophilicity (Stefan Iglauer, 2017). The higher the salt concentration the higher the contact angel. The effect increases with higher valency (i.e. higher charge concentration per unit of volume) and with higher ionic strength: