Refine
Document Type
- Article (1)
- Doctoral Thesis (1)
Language
- English (2) (remove)
Is part of the Bibliography
- yes (2) (remove)
Keywords
- Fluid (2) (remove)
Institute
- Institut für Geowissenschaften (2) (remove)
In this study we report the structure of supercritical H2O-SiO2 fluid composed of 50 mol% H2O and 50 mol% SiO2 at 3000 K and 2400 K. investigated by means of ab initio molecular dynamics of models comprising 192 and 96 atoms. The density is set constant to 138 g/cm(3), which yields a pressure of 4.3 GPa at 3000 K and 3.6 GPa at 2400 K. Throughout the trajec[ories, water molecules are formed and dissociated via the network modifying reaction 2 SiOH = SiOSi + H2O The calculation of the reaction constant K- [OH](2)/[H2O][O2-] is carried out on the basis of the experimentally relevant Q ' species notation and agrees well with an extrapolation of experimental data to 3000 K. After quench from 3000 K to 2400 K, the degree of polymerization of the silicate network in the 192-atom models increases noticeably within several tens of picoseconds, accompanied by release of molecular H2O. An unexpected opposite trend is observed in smaller 96-atom models, due to a finite size effect, as several uncorrelated models of 192 and 96 atoms indicate. The temperature-dependent slowing down of the H2O-silica interaction dynamics is described on the basis of the bond autocorrelation function. (C) 2016 Elsevier B.V. All rights reserved.
Fluids in the Earth's crust can move by creating and flowing through fractures, in a process called `hydraulic fracturing’. The tip-line of such fluid-filled fractures grows at locations where stress is larger than the strength of the rock. Where the tip stress vanishes, the fracture closes and the fluid-front retreats. If stress gradients exist on the fracture's walls, induced by fluid/rock density contrasts or topographic stresses, this results in an asymmetric shape and growth of the fracture, allowing for the contained batch of fluid to propagate through the crust.
The state-of-the-art analytical and numerical methods to simulate fluid-filled fracture propagation are two-dimensional (2D). In this work I extend these to three dimensions (3D). In my analytical method, I approximate the propagating 3D fracture as a penny-shaped crack that is influenced by both an internal pressure and stress gradients. In addition, I develop a numerical method to model propagation where curved fractures can be simulated as a mesh of triangular dislocations, with the displacement of faces computed using the displacement discontinuity method. I devise a rapid technique to approximate stress intensity and use this to calculate the advance of the tip-line. My 3D models can be applied to arbitrary stresses, topographic and crack shapes, whilst retaining short computation times.
I cross-validate my analytical and numerical methods and apply them to various natural and man-made settings, to gain additional insights into the movements of hydraulic fractures such as magmatic dikes and fluid injections in rock. In particular, I calculate the `volumetric tipping point’, which once exceeded allows a fluid-filled fracture to propagate in a `self-sustaining’ manner. I discuss implications this has for hydro-fracturing in industrial operations. I also present two studies combining physical models that define fluid-filled fracture trajectories and Bayesian statistical techniques. In these studies I show that the stress history of the volcanic edifice defines the location of eruptive vents at volcanoes. Retrieval of the ratio between topographic to remote stresses allows for forecasting of probable future vent locations. Finally, I address the mechanics of 3D propagating dykes and sills in volcanic regions. I focus on Sierra Negra volcano in the Gal\'apagos islands, where in 2018, a large sill propagated with an extremely curved trajectory. Using a 3D analysis, I find that shallow horizontal intrusions are highly sensitive to topographic and buoyancy stress gradients, as well as the effects of the free surface.