Probabilistic analysis and comparison of stress-dependent rock physics models

A rock physics model attempts to account for the nonlinear stress dependence of seismic velocity by relating changes in stress and strain to changes in seismic velocity and anisotropy. Understanding and being able to model this relationship is crucial for any time-lapse geophysical or geohazard modelling scenario. In this study, we take a number of commonly used rock physics models and assess their behaviour and stability when applied to stress versus velocity measurements of a large (dry) core data set of different lithologies. We invert and calibrate each model and present a database of models for over 400 core samples. The results of which provide a useful tool for setting a priori parameter constraints for future model inversions. We observe that some models assume an increase in VP/VS ratio (hence Poisson’s ratio) with stress. A trait not seen for every sample in our data set. We demonstrate that most model parameters are well constrained. However, third-order elasticity models become ill-posed when their equations are simplified for an isotropic rock. We also find that third-order elasticity models are limited by their approximation of an exponential relationship via functions that lack an exponential term. We also argue that all models are difficult to parametrize without the availability of core data. Therefore, we derive simple relationships between model parameters, core porosity and clay content. We observe that these relationship are suitable for estimating seismic velocities of rock but poor when comes to predicting changes related to effective stress. The findings of this study emphasize the need for improvement to models if quantitatively accurate predictions of time-lapse velocity and anisotropy are to be made. Certain models appear to better fit velocity depth log data than velocity–stress core data. Thus, there is evidence to suggest a limitation in core data as a representation of the stress dependence of the subsurface. The differences in the stress dependence of the subsurface compared to that measured under laboratory conditions could potentially be significant. Although potentially difficult to investigate, its importance is of great significance if we wish to accurately interpret the stress dependence of subsurface seismic velocities

Reservoir stress path and induced seismic anisotropy: Results from linking coupled fluid-flow/geomechanical simulation with seismic modelling

We present a workflow linking coupled fluid-flow and geomechanical simulation with seismic modelling to predict seismic anisotropy induced by nonhydrostatic stress changes. We generate seismic models from coupled simulations to examine the relationship between reservoir geometry, stress path and seismic anisotropy. The results indicate that geometry influences the evolution of stress, which leads to stress-induced seismic anisotropy. Although stress anisotropy is high for the small reservoir, the effect of stress arching and the ability of the side-burden to support the excess load limit the overall change in effective stress and hence seismic anisotropy. For the extensive reservoir, stress anisotropy and induced seismic anisotropy are high. The extensive and elongate reservoirs experience significant compaction, where the inefficiency of the developed stress arching in the side-burden cannot support the excess load. The elongate reservoir displays significant stress asymmetry, with seismic anisotropy developing predominantly along the long-edge of the reservoir. We show that the link between stress path parameters and seismic anisotropy is complex, where the anisotropic symmetry is controlled not only by model geometry but also the nonlinear rock physics model used. Nevertheless, a workflow has been developed to model seismic anisotropy induced by non-hydrostatic stress changes, allowing field observations of anisotropy to be linked with geomechanical models

Can good fat boost your fitness level?

Fish oil supplements may seem like a relatively recent health fad but they have actually been produced in the UK on a large scale since 1935 by the company Seven Seas Ltd. Since then, the fish oil supplement market has continued to grow, with many beneficial effects claimed for health

Complementary hydro-mechanical coupled finite/discrete element and microseismic modelling to predict hydraulic fracture propagation in tight shale reservoirs

This paper presents a novel approach to predict the propagation of hydraulic fractures in tight shale reservoirs. Many hydraulic fracture modelling schemes assume that the fracture direction is pre-seeded in the problem domain discretization. This is a severe limitation as the reservoir often contains large numbers of pre-existing fractures that strongly influence the direction of the propagating fracture. To circumvent these shortcomings a new fracture modelling treatment is proposed where the introduction of discrete fracture surfaces is based on new and dynamically updated geometrical entities rather than the topology of the underlying spatial discretization. Hydraulic fracturing is an inherently coupled engineering problem with interactions between fluid flow and fracturing when the stress state of the reservoir rock attains a failure criterion. This work follows a staggered hydro-mechanical coupled finite/discrete element approach to capture the key interplay between fluid pressure and fracture growth. In field practice the fracture growth is hidden from the design engineer and microseismicity is often used to infer hydraulic fracture lengths and directions. Microsesimic output can also be computed from changes of the effective stress in the geomechanical model and compared against field microseismicity. A number of hydraulic fracture numerical examples are presented to illustrate the new technology

The host galaxies and explosion sites of long-duration gamma-ray bursts: Hubble Space Telescope near-infrared imaging

We present the results of a Hubble Space Telescope WFC3/F160WSnapshot survey of the host galaxies of 39 long-duration gamma-ray bursts (LGRBs) at z < 3. We have non-detections of hosts at the locations of four bursts. Sufficient accuracy to astrometrically align optical afterglowimages and determine the location of the LGRB within its hostwas possible for 31/35 detected hosts. In agreement with other work, we find the luminosity distribution of LGRB hosts is significantly fainter than that of a star formation rate-weighted field galaxy sample over the same redshift range, indicating LGRBs are not unbiasedly tracing the star formation rate. Morphologically, the sample of LGRB hosts is dominated by spiral-like or irregular galaxies. We find evidence for evolution of the population of LGRB hosts towards lower luminosity, higher concentrated hosts at lower redshifts. Their half-light radii are consistent with other LGRB host samples where measurements were made on rest-frame UV observations. In agreement with recent work, we find their 80 per cent enclosed flux radii distribution to be more extended than previously thought, making them intermediate between core-collapse supernova (CCSN) and superluminous supernova (SLSN) hosts. The galactocentric projectedoffset distribution confirms LGRBs as centrally concentrated, much more so than CCSNe and similar to SLSNe. LGRBs are strongly biased towards the brighter regions in their host light distributions, regardless of their offset. We find a correlation between the luminosity of the LGRB explosion site and the intrinsic column density, NH, towards the burst. © 2017 The Authors

A multimethod Global Sensitivity Analysis to aid the calibration of geomechanical models via time-lapse seismic data

Time-lapse seismic attributes are used extensively in the history matching of production simulator models. However, although proven to contain information regarding production induced stress change, it is typically only loosely (i.e. qualitatively) used to calibrate geomechanical models. In this study we conduct a multimethod Global Sensitivity Analysis (GSA) to assess the feasibility and aid the quantitative calibration of geomechanical models via near-offset time-lapse seismic data. Specifically, the calibration of mechanical properties of the overburden. Via the GSA, we analyse the near-offset overburden seismic traveltimes from over 4000 perturbations of a Finite Element (FE) geomechanical model of a typical High Pressure High Temperature (HPHT) reservoir in the North Sea. We find that, out of an initially large set of material properties, the near-offset overburden traveltimes are primarily affected by Young's modulus and the effective stress (i.e. Biot) coefficient. The unexpected significance of the Biot coefficient highlights the importance of modelling fluid flow and pore pressure outside of the reservoir. The FE model is complex and highly nonlinear. Multiple combinations of model parameters can yield equally possible model realizations. Consequently, numerical calibration via a large number of random model perturbations is unfeasible. However, the significant differences in traveltime results suggest that more sophisticated calibration methods could potentially be feasible for finding numerous suitable solutions. The results of the time-varying GSA demonstrate how acquiring multiple vintages of time-lapse seismic data can be advantageous. However, they also suggest that significant overburden near-offset seismic time-shifts, useful for model calibration, may take up to 3 yrs after the start of production to manifest. Due to the nonlinearity of the model behaviour, similar uncertainty in the reservoir mechanical properties appears to influence overburden traveltime to a much greater extent. Therefore, reservoir properties must be known to a suitable degree of accuracy before the calibration of the overburden can be considered

Architecture of the chromatin remodeler RSC and insights into its nucleosome engagement.

Eukaryotic DNA is packaged into nucleosome arrays, which are repositioned by chromatin remodeling complexes to control DNA accessibility. The Saccharomyces cerevisiae RSC (Remodeling the Structure of Chromatin) complex, a member of the SWI/SNF chromatin remodeler family, plays critical roles in genome maintenance, transcription, and DNA repair. Here, we report cryo-electron microscopy (cryo-EM) and crosslinking mass spectrometry (CLMS) studies of yeast RSC complex and show that RSC is composed of a rigid tripartite core and two flexible lobes. The core structure is scaffolded by an asymmetric Rsc8 dimer and built with the evolutionarily conserved subunits Sfh1, Rsc6, Rsc9 and Sth1. The flexible ATPase lobe, composed of helicase subunit Sth1, Arp7, Arp9 and Rtt102, is anchored to this core by the N-terminus of Sth1. Our cryo-EM analysis of RSC bound to a nucleosome core particle shows that in addition to the expected nucleosome-Sth1 interactions, RSC engages histones and nucleosomal DNA through one arm of the core structure, composed of the Rsc8 SWIRM domains, Sfh1 and Npl6. Our findings provide structural insights into the conserved assembly process for all members of the SWI/SNF family of remodelers, and illustrate how RSC selects, engages, and remodels nucleosomes