Hydraulic Injection‐Induced Velocity Changes Revealed by Surface Wave Coda and Polarization Data at a Shale Play Site in Southwest China

2020 ◽  
Vol 125 (10) ◽  
Author(s):  
Yan Zhang ◽  
Fenglin Niu ◽  
Kai Tao ◽  
Jieyuan Ning ◽  
Haichao Chen ◽  
...  
Author(s):  
N. H. Ngo ◽  
H. Tran ◽  
R. R. Gamache ◽  
J. M. Hartmann

A short overview of recent results on the effects of pressure (collisions) regarding the shape of isolated infrared lines of water vapour is presented. The first part of this study considers the basic collisional quantities, which are the pressure-broadening and -shifting coefficients, central parameters of the Lorentzian (and Voigt) profile and thus of any sophisticated line-shape model. Through comparisons of measured values with semi-classical calculations, the influences of the molecular states (both rotational and vibrational) involved and of the temperature are analysed. This shows the relatively unusual behaviour of H 2 O broadening, with evidence of a significant vibrational dependence and the fact that the broadening coefficient (in cm −1 atm −1 ) of some lines increases with temperature. In the second part of this study, line shapes beyond the Voigt model are considered, thus now taking ‘velocity effects’ into account. These include both the influence of collisionally induced velocity changes that lead to the so-called Dicke narrowing and the influence of the dependence of collisional parameters on the speed of the radiating molecule. Experimental evidence of deviations from the Voigt shape is presented and analysed. The interest of classical molecular dynamics simulations, to model velocity changes, together with semi-classical calculations of the speed-dependent collisional parameters for line-shape predictions from ‘first principles’, are discussed.


1993 ◽  
Vol 247 ◽  
pp. 589-601 ◽  
Author(s):  
Yehuda Agnon

Nonlinear interaction of surface wave trains is studied. Zakharov's kernel is extended to include the vicinity of trio resonance. The forced wave amplitude and the wave velocity changes are then first order rather than second order. The model is applied to remove near-resonance singularities in expressions for the change of speed of one wave train in the presence of another. New results for Wilton ripples and the drift current and setdown in shallow water waves are readily derived. The ideas are applied to the derivation of forced waves in the vicinity of quartet and quintet resonance in an evolving wave field.


Geophysics ◽  
2009 ◽  
Vol 74 (1) ◽  
pp. D25-D36 ◽  
Author(s):  
Rodrigo Felício Fuck ◽  
Andrey Bakulin ◽  
Ilya Tsvankin

Time-lapse traveltime shifts of reflection events recorded above hydrocarbon reservoirs can be used to monitor production-related compaction and pore-pressure changes. Existing methodology, however, is limited to zero-offset rays and cannot be applied to traveltime shifts measured on prestack seismic data. We give an analytic 3D description of stress-related traveltime shifts for rays propagating along arbitrary trajectories in heterogeneous anisotropic media. The nonlinear theory of elasticity helps to express the velocity changes in and around the reservoir through the excess stresses associated with reservoir compaction. Because this stress-induced velocity field is both heterogeneous and anisotropic, it should be studied using prestack traveltimes or amplitudes. Then we obtain the traveltime shifts by first-order perturbation of traveltimes that accounts not only for the velocity changes but also for 3D deformation of reflectors. The resulting closed-form expression can be used efficiently for numerical modeling of traveltime shifts and, ultimately, for reconstructing the stress distribution around compacting reservoirs. The analytic results are applied to a 2D model of a compacting rectangular reservoir embedded in an initially homogeneous and isotropic medium. The computed velocity changes around the reservoir are caused primarily by deviatoric stresses and produce a transversely isotropic medium with a variable orientation of the symmetry axis and substantial values of the Thomsen parameters [Formula: see text] and [Formula: see text]. The offset dependence of the traveltime shifts should play a crucial role in estimating the anisotropy parameters and compaction-related deviatoric stress components.


1999 ◽  
Vol 138 (1) ◽  
pp. 205-220 ◽  
Author(s):  
Kazunori Yoshizawa ◽  
Kiyoshi Yomogida ◽  
Seiji Tsuboi

2021 ◽  
Author(s):  
◽  
Rachel Heckels

<p>Ambient seismic noise is used to examine the spatial and temporal surface wave velocity structures and ambient seismic noise fields in the vicinity of different fault zone environments. This study focuses on two distinct regions of central South Island, New Zealand. The Canterbury Plains is a sedimentary basin with many minor faults, which was considered to have low seismic hazard prior to the 2010 – 2011 Canterbury earthquake sequence. We focus on the time period immediately following the 2010 Darfield earthquake, which ruptured the previously unmapped Greendale Fault. The second region of interest is the central Southern Alps. The locked portion of the Alpine Fault currently poses one of the largest seismic hazards for New Zealand. The wealth of data from both permanent and temporary seismic deployments in these regions make them ideal areas in which to assess the effectiveness of ambient noise for velocity modelling in regions surrounding faults at different stages of their seismic cycles.  Temporal velocity changes are measured following the Mw 7.1 Darfield earthquake of 4 September 2010 in the Canterbury Plains. Nine-component cross-correlations are computed from temporary and permanent seismic stations lying on and surrounding the Greendale Fault. Using the Moving-Window Cross-Spectral method, surface wave velocity changes are calculated for the four months immediately following the earthquake until 10 January 2011, for 0.1 — 1.0 Hz. An average increase in seismic velocity of 0.14 ± 0.04 % is determined throughout the region, providing the first such estimate of postseismic relaxation rates in Canterbury. Depth analyses further showed that velocity changes are confined to the uppermost 5 km of the subsurface and we attribute this to postseismic relaxation via crack-healing of the Greendale Fault and throughout the surrounding region.  Rayleigh and Love wave dispersion is examined throughout the Canterbury region. Multi-component cross-correlation functions are analysed for group and phase dispersion curves. These are inverted using frequency-time analysis for 2-D phase and group velocity maps of Rayleigh and Love waves. A high-velocity zone to the southeast of the region coincides with volcanic rocks of Banks Peninsula. Dispersion curves generated from the surface wave tomography are further inverted for one-dimensional shear velocity profiles. These models show a thin, low-velocity near surface layer consistent with the basin sediments, which thins towards the foothills of the Southern Alps. A near-surface damage zone is identified along the length of the Greendale Fault, with consistent reduced Vs velocities to depth of up to 5 km.  Surface and shear wave velocity maps are computed for the central Southern Alps to image the seismic structure of the region. Tomographic surface maps at periods of 5 – 12 s are produced from dispersion measurements of three-component cross-correlation functions. At periods of 5 – 8 s a strong NE-SW trending velocity contrast highlights the Alpine Fault. One-dimensional shear velocity models, computed from the surface wave maps, are in agreement with previous models produced by other conventional methods. An analysis of surface wave amplitudes through signal-to-noise ratios of cross-correlations reveals strong directional effects. Calculated signal-to-noise ratios are up to eight times higher for surface waves travelling north-west than for waves travelling to the south or east. We attribute this to a combination of more energetic ocean wave signals from the Southern Ocean compared to the Tasman Sea.</p>


2021 ◽  
Author(s):  
◽  
Francesco Civilini

<p>We present three projects that use different bandwidths of the ambient noise spectrum to solve geophysical problems. Specifically, we use signals within the noise field to determine surface and shear wave velocities, image the shallow and deep crust, and monitor time-dependent deformation resulting from geothermal fluid injection and extraction.  Harrat Al-Madinah, a Cenozoic bimodal alkaline volcanic field in west-central Saudi Arabia, is imaged using shear-velocities obtained from natural ambient seismic noise. To our knowledge, this project is the first analysis of Saudi Arabia structure using ambient noise methods. Surface wave arrivals are extracted from a year's worth of station-pair cross-correlations, which are approximations of the empirical Green's function of the interstation path. We determine group and phase velocity surface wave dispersion maps with a 0.1 decimal degree resolution and resolve a zone of slow surface wave velocity south-east of the city of Medina, which is spatially correlated with the most recent historical eruption (the 1256 CE Medina eruption). Dispersion curves are calculated at each grid-point of the surface-wave velocity maps and inverted to obtain measurements of shear-velocity with depth. The 1D velocity models are then used to produce average shear-velocity models for the volcanic field. A shear-velocity increase ranging from 0.5 to 1.0 km/s, suggesting a layer interface, is detected at approximately 20 km depth and compared to P-wave measurement from a previous refraction study. We compute cross-section profiles by interpolating the inversions into a pseudo-3D model and resolve a zone of slow shear-velocity below the 1256 CE eruption location. These areas are also spatially correlated with low values of Bouguer gravity. We hypothesize that the low shear-velocity and gravity measurements are caused by fluids and fractures created from prior volcanic eruptions.   We use the coda of cross-correlations extracted from ambient noise to determine shear-velocity changes at Rotokawa and Ngatamariki, two electricity producing geothermal fields located in the North Island of New Zealand. Stacks of cross correlations between stations prior to the onset of production are compared to cross correlations of moving stacks in time periods of well stimulation and the onset of electricity production using the Moving Window Cross Spectral technique. An increase between 0.05% to 0.1% of shear-velocity is detected at Rotokawa coinciding with an increase of injection. The shear-velocity subsequently decreases by approximately 0.1% when the rate of production surpasses the rate of injection. A similar amplitude shear-velocity increase is detected at Ngatamariki during the beginning of injection. After the initial increase, the shear-velocity at Ngatamariki fluctuates in response to differences in injection and production rates. A straight-ray pseudo-tomography analysis is conducted at the geothermal fields, which reveals that localized positive velocity changes are co-located with injection wells.  Lastly, we use ambient noise and active sources at the Ngatamariki geothermal field to determine the structure of the top 200 meters using the Refraction Microtremor technique. We deployed a linear 72-channel array of vertical geophones with ten meter spacing at two locations of the geothermal field and determine average 1D and 2D shear-velocity profiles. We were able to image depths between 57 to 93 meters for 2D profiles and up to 165 meters for 1D profiles. A shear-velocity anomaly was detected across one of the lines that coincided with the inferred location of a fault determined from nearby well logs. This suggests that the method can be used to cheaply and quickly constrain near-surface geology at geothermal fields, where ambient noise is abundant and typical reflection and refraction surveys require large inputs of energy and are hindered by attenuation and scattering in near-surface layers.</p>


2021 ◽  
Vol 11 (4) ◽  
pp. 1-18
Author(s):  
Q. Bai ◽  
H. Konietzky

This contribution proposes a numerical microstructural modeling approach to investigate stress-induced seismic velocity changes on anisotropic rocks. By introducing pre-existing cracks with preferential orientations in bonded-particle assemblies, the transverse isotropic structure of the Whitby Mudstone is simulated. Using power-law distributed aperture and calibrated micro-properties, we successfully reproduce stress-dependent velocity changes on Whitby Mudstones with different anisotropic angles in relation to the applied loads. The proposed model also duplicates the directional dependence of wave speed with respect to the bedding plane as expected theoretically. The numerical models show that velocity increase results from the closure of pre-existing cracks due to load increase. Direct relations are established between velocity changes and opened crack density (or crack closure), which displays a similar tendency compared with theoretical predictions. This relation can be used to quantify the micromechanisms behind the velocity changes. The proposed model provides the ability to directly examine the micro-processes underlying velocity changes.


2021 ◽  
Author(s):  
Luc Illien ◽  
Christoph Sens-Schönfelder ◽  
Christoff Andermann ◽  
Odin Marc ◽  
Kristen Cook ◽  
...  

&lt;p&gt;Following the passage of seismic waves, most geomaterials experience non-linear mesoscopic elasticity (&lt;em&gt;NLME&lt;/em&gt;). This is described by a drop in elastic moduli that precedes a subsequent recovery of physical properties over a relaxation timescale. Thanks to the development of seismic interferometry techniques that allows for the continuous monitoring of relative seismic velocity changes &lt;em&gt;&amp;#948;v&lt;/em&gt; in the subsurface, observations of &lt;em&gt;NLME&lt;/em&gt; (&lt;em&gt;&amp;#948;v&lt;/em&gt;&lt;sub&gt;&lt;em&gt;NLME&lt;/em&gt;&lt;/sub&gt;) in the field are now numerous. In parallel, a growing community uses seismic interferometry to monitor velocity changes induced by seasonal hydrological variations (&lt;em&gt;&amp;#948;v&lt;sub&gt;hydro&lt;/sub&gt;&lt;/em&gt;). Monitoring of these variations are often independently done and a linear superposition of both effects is mostly assumed when decomposing the observed &lt;em&gt;&amp;#948;v&lt;/em&gt; signal (&lt;em&gt;&amp;#948;v&lt;/em&gt; =&amp;#160; &lt;em&gt;&amp;#948;v&lt;sub&gt;NLME&lt;/sub&gt;&lt;/em&gt; + &lt;em&gt;&amp;#948;v&lt;sub&gt;hydro&lt;/sub&gt;&lt;/em&gt;). However, transient hydrological behaviour following co-seismic ground shaking has been widely reported in boreholes measurements and streamflow, which suggests that&amp;#160; &lt;em&gt;&amp;#948;v&lt;sub&gt;hydro&lt;/sub&gt;&lt;/em&gt; may be impacted by the transient variation of material properties caused by &lt;em&gt;NLME&lt;/em&gt;. In this presentation, we attempt to characterize the relative seismic velocity variations &lt;em&gt;&amp;#948;v&lt;/em&gt; retrieved from a small dense seismic array in Nepal that was deployed in the aftermath of the &amp;#160;2015 Mw 7.8 Gorkha earthquake and that is prone to highly variable hydrological conditions. We first investigated the effect of aftershocks in computing &lt;em&gt;&amp;#948;v&lt;/em&gt; at a 10-minute resolution centered around significant ground shaking events. After correcting &lt;em&gt;&amp;#948;v&lt;/em&gt; for &lt;em&gt;NLME&lt;/em&gt; caused by the Gorkha earthquake and its subsequent aftershocks, we test whether the corresponding residuals are in agreement with the background hydrological behaviour which we inferred from a calibrated hydrological model. This is not the case and we find that transient hydrological properties improve the data description in the early phase after the mainshock. We report three distinct relaxation time scales that are relevant for the recovery of seismic velocity at our field site:&amp;#160; &lt;strong&gt;1.&lt;/strong&gt; A long time scale activated by the main shock of the Gorkha earthquake (~1 year) &lt;strong&gt;2.&lt;/strong&gt; A relatively short timescale (1-3 days) that occurs after moderate aftershocks. &lt;strong&gt;3.&lt;/strong&gt; An intermediate timescale (4-6 months) during the 2015 monsoon season that corresponds to the recovery of the hydrological system. This timescale could correspond to an enhanced permeability caused by Gorkha ground shaking. Our study demonstrates the capability of seismic interferometry to monitor transient hydrological properties after earthquakes at a spatial scale that is not available with classical hydrological measurements. This investigation demands calibrated hydrological models and a framework in which the different forcing of &lt;em&gt;&amp;#948;v&lt;/em&gt; are coupled.&lt;/p&gt;


Sign in / Sign up

Export Citation Format

Share Document