This UKCCSRC (UK Carbon Capture and Storage Research Centre) Call 1 project involved the development, testing and validation of a two-fluid transient flow model for simulating outflow following the failure of high pressure CO2 pipelines is presented. The project made use of experimental data and used experimental data available from other UK/EC funded projects. The model developed accounts for thermal and mechanical non-equilibrium effects during depressurisation by utilising simple constitutive relations describing inter-phase mass, heat and momentum transfer in terms of relaxation to equilibrium. Pipe wall/fluid heat exchange on the other hand is modelled by coupling the fluid model with a finite difference transient heat conduction model. This paper describes the model, the details of its numerical solution and its validation as well as parametric analysis of relevant parameters. http://www.sciencedirect.com/science/article/pii/S1750583614002394, DOI: 10.1016/j.ijggc.2014.08.013. UKCCSRC grant UKCCSRC-C1-07.
The partitioning coefficients of water between iron and silicate melts at 20, 50, 90 and 135 gigapascals (corresponding to 2800, 3500, 3900 and 4200 kelvin) were calculated by using ab initio molecular dynamics and thermodynamic integration techniques. The Gibbs free energy of a series of iron and silicate melts with different concentrations of H2/H2O were calculated. Then the chemical potentials of H2/H2O were derived from the concentration dependent Gibbs free energies at each pressure temperature. The partitioning coefficients can be calculated by equating the chemical potential of H2/H2O in iron and silicate melts. The Weeks-Chandler-Andersen (WCA) system with established thermodynamics was used as the reference.
Synchrotron X-radiography (images) and diffraction data collected to measure anelasticity of zinc. NERC grant NE/H016309/1 - Experimental determination of mantle rheology. NERC grant NE/L006898/1 - The strength of the lower mantle.
Synchrotron X-radiography (images) and diffraction data collected to measure rheology of Quartz coesite and stishovite.
Synthetic seismic waveforms computed using the spectral element method for 1-D and 3-D Earth models for a variety of scenarios of structures in the deep Earth's interior.
This poster on the UKCCSRC (UK Carbon Capture and Storage Research Centre) Call 1 project, Multi-Phase Flow Modelling for Hazardous Assessment, was presented at the Cambridge Biannual, 02.04.14. Grant number: UKCCSRC-C1-07.
The dataset contains 15 plots and data for time-dependent pressures and temperatures at various locations along a 2582-m-long well and at various simulation times. The realistic scenarios taken into considerations are applied to the Goldeneye depleted reservoir in the North Sea. Pure CO2 is injected into the well and then discharged in the Goldeneye reservoir. Six different scenarios are considered: three different injection durations (linear ramp-up of the inlet mass flow rate from 0 to 33.5 kg/s over 5 minutes, 30 minutes, and 2 hours) and two different upstream temperatures (278.15 K and 283.15 K). Data is currently restricted until publication.
Processed SAR interferograms for the Wells, Nevada earthquake. Grant abstract: How do earthquakes happen? Understanding the nature of earthquakes is a key fundamental question in Geociences that holds many implications for society. Earthquakes are typically associated with a sudden release of energy that has slowly accumulated over hundreds to thousands of years, being strongly controlled by friction in faults buried several kilometers beneath our feet under quite extreme conditions. For example, the amount of heat produced in just a few seconds is such that it can dramatically change the nature of the fault zone near the sliding surface. Moreover, there is abundant evidence of substantial frictional weakening of faults (i.e., fault strength weakens with increasing slip or slip rate) during earthquakes. However, there are still many open questions related to earthquake source processes: How similar are earthquakes in different temperature-pressure conditions? What is the earthquake's energy budget, which controls the intensity of ground motions? What are the physical mechanisms responsible for fault weakening? Recent progress in seismological imaging methods, theoretical fracture mechanics and rupture dynamics simulations can help solve these questions. Huge volumes of freely available seismic and geodetic data from around the world now allow the routine calculation of earthquake models where earthquakes are typically described as single space-time points. Time is now ripe for systematically building robust, more detailed seismic models bearing information on earthquake's physics by using recently developed sophisticated modelling tools along with high-quality images of the 3-D Earth's interior structure enabled by high performance computing facilities. Moreover, it is now possible to model ruptures theoretically in detail using both analytical fracture mechanics calculations and numerical rupture dynamics simulations, and, for example, estimate the fault temperature during the rupture process, which is the most direct way to quantify friction. However, systematic quantitative links between these calculations and seismological observations are still lacking. This project addresses these issues through a coordinated effort involving seismology and rock mechanics aiming at estimating fault temperature rise during earthquakes from new macroscopic seismic source models. We will use advanced seismic source imaging methods to build a new set of robust kinematic, static and dynamic earthquake source parameters for a large selected set of global earthquakes (e.g., average fault length, width, rupture speed and time history, stress drop, radiated and fracture energy). These solutions will then be used as input parameters to estimate fault temperature using analytical and numerical rupture dynamics calculations. This will lead to an improved understanding of how local fault processes occurring at scales from few microns to tens of centimetres translate into macroscopic seismological properties, how energy is partitioned during earthquakes and which are the mechanisms responsible for fault weakening. Ultimately this project will shed new light on many basic questions in earthquake science such as the similarity of earthquakes in different P-T conditions and the potential geological record left by ruptures (e.g., melt). More broadly, this project will benefit hazard models and any studies relying on accurate earthquake source parameters such as studies in seismic tomography, active tectonics and microseismicity (e.g., associated with hydraulic fracturing).
Data generated at UCL on a conventional triaxial apparatus. The data includes the 'raw' data, as well as semi-processed data used to plot figures and draw main conclusions. This dataset is used and fully described/interpreted in the paper: M. Jefferd, N. Brantut, P.G. Meredith and T.M. Mitchell, Compactive Deformation of Sandstone under Crustal Pressure and Temperature Conditionsserpentinite, submitted to J. Geophys. Res. And in the UCL PhD Thesis M.Jefferd, Sandstone under Crustal Pressure and Temperature
Detrital zircon age data, details of Expedition 362 samples . For more information see published report, https://doi.org/10.1016/j.epsl.2017.07.019 IODP Sites U1480–U1481, located on the Indian oceanic plate, east of the NinetyEast Ridge and west of the north Sumatran subduction margin Site U1480 ~ 3°2.0447'N 91°36.3481'E 4147.5 Site U1481 ~ 2°45.261'N 91°45.5771'E