Sensitivity of time lapse seismic data to the compliance of hydraulic fractures

We study the sensitivity of seismic waves to changes in the fracture normal and tangential compliances by analyzing the fracture sensitivity wave equation, which is derived by differentiating the elastic wave equation with respect to the fracture compliance. The sources for the sensitivity wavefield are the sensitivity moments, which are functions of fracture compliance, background elastic properties and the stress acting on the fracture surface. Based on the analysis of the fracture sensitivity wave equation, we give the condition for the weak scattering approximation to be valid for fracture scattering. Under the weak scattering approximation, we find that the percentage change of fracture compliance in hydraulic fracturing is equal to the percentage change of the recorded time-lapse seismic data. This could provide a means for monitoring the opening/closing of fractures in hydraulic fracturing through time-lapse seismic surveys.Eni-MIT Energy Initiative Founding Member Progra

Detectionof gas hydrates infaults using azimuthal seismic velocity analysis,Vestnesa Ridge, W-Svalbard Margin

Accepted for publication in Journal of Geophysical Research. Solid Earth. Copyright 2020 American Geophysical Union. Further reproduction or electronic distribution is not permitted.Joint analysis of electrical resistivity and seismic velocity data is primarily used to detect the presence of gas hydrate‐filled faults and fractures. In this study, we present a novel approach to infer the occurrence of structurally‐controlled gas hydrate accumulations using azimuthal seismic velocity analysis. We perform this analysis using ocean‐bottom seismic (OBS) data at two sites on Vestnesa Ridge, W‐Svalbard Margin. Previous geophysical studies inferred the presence of gas hydrates at shallow depths (up to ~190‐195 m below the seafloor) in marine sediments of Vestnesa Ridge. We analyze azimuthal P‐wave seismic velocities in relation with steeply‐dipping near surface faults to study structural controls on gas hydrate distribution. This unique analysis documents directional changes in seismic velocities along and across faults. P‐wave velocities are elevated and reduced by ~0.06‐0.08 km/s in azimuths where the raypath plane lies along the fault plane in the gas hydrate stability zone (GHSZ) and below the base of the GHSZ, respectively. The resulting velocities can be explained with the presence of gas hydrate‐ and free gas‐filled faults above and below the base of the GHSZ, respectively. Moreover, the occurrence of elevated and reduced (>0.05 km/s) seismic velocities in groups of azimuths bounded by faults, suggests compartmentalization of gas hydrates and free gas by fault planes. Results from gas hydrate saturation modelling suggest that these observed changes in seismic velocities with azimuth can be due to gas hydrate saturated faults of thickness greater than 20 cm and considerably smaller than 300 cm

Seismic attenuation in fractured porous media: insights from a hybrid numerical and analytical model

Seismic attenuation in fluid-saturated porous rocks can occur by geometric spreading, wave scattering or the internal dissipation of energy, most likely due to the squirt-flow mechanism. In principle, the pattern of seismic attenuation recorded on an array of sensors contains information about the medium, in terms of material heterogeneity and anisotropy, as well as material properties such as porosity, crack density, and pore-fluid composition and mobility. In practice, this inverse problem is challenging. Here we provide some insights into the effects of internal dissipation by analysing synthetic data produced by a hybrid numerical and analytical model for seismic wave propagation in a fractured medium embedded within a layered geological structure. The model is made up of one anisotropic and three isotropic horizontal layers. The anisotropic layer consists of a porous, fluid-saturated material containing vertically aligned inclusions representing a set of fractures. This combination allows squirt-flow to occur between the pores in the matrix and the model fractures. Our results show that the fluid mobility and the associated relaxation time of the fluid-pressure gradient control the frequency range over which attenuation occurs. This induced attenuation increases with incidence angle and azimuth away from the fracture strike-direction. Azimuthal variations in the induced attenuation are elliptical allowing the fracture orientations to be obtained from the axes of the ellipse. These observations hold out the potential of using seismic attenuation as an additional diagnostic in the characterisation of rock formations for a variety of applications including hydrocarbon exploration and production, subsurface storage of CO2, and geothermal energy extraction