Seismic attenuation and velocity dispersion in fractured rocks: The role played by fracture contact areas

The presence of fractures in fluid‐saturated porous rocks is usually associated with strong seismic P‐wave attenuation and velocity dispersion. This energy dissipation can be caused by oscillatory wave‐induced fluid pressure diffusion between the fractures and the host rock, an intrinsic attenuation mechanism generally referred to as wave‐induced fluid flow. Geological observations suggest that fracture surfaces are highly irregular at the millimetre and sub‐millimetre scale, which finds its expression in geometrical and mechanical complexities of the contact area between the fracture faces. It is well known that contact areas strongly affect the overall mechanical fracture properties. However, existing models for seismic attenuation and velocity dispersion in fractured rocks neglect this complexity. In this work, we explore the effects of fracture contact areas on seismic P‐wave attenuation and velocity dispersion using oscillatory relaxation simulations based on quasi‐static poroelastic equations. We verify that the geometrical and mechanical details of fracture contact areas have a strong impact on seismic signatures. In addition, our numerical approach allows us to quantify the vertical solid displacement jump across fractures, the key quantity in the linear slip theory. We find that the displacement jump is strongly affected by the geometrical details of the fracture contact area and, due to the oscillatory fluid pressure diffusion process, is complex‐valued and frequency‐dependent. By using laboratory measurements of stress‐induced changes in the fracture contact area, we relate seismic attenuation and dispersion to the effective stress. The corresponding results do indeed indicate that seismic attenuation and phase velocity may constitute useful attributes to constrain the effective stress. Alternatively, knowledge of the effective stress may help to identify the regions in which wave induced fluid flow is expected to be the dominant attenuation mechanism.     

Virtual signals combination – phase analysis for 2D and 3D data representation (acoustic medium)

We formulate the Kirchhoff‐Helmholtz representation theory for the combination of seismic interferometry signals synthesized by cross‐correlation and by cross‐convolution in acoustic media. The approach estimates the phase of the virtual reflections from the boundary encompassing a volume of interest and subtracts these virtual reflections from the total seismic‐interferometry wavefield. The reliability of the combination result, relevant for seismic exploration, depends on the stationary‐phase and local completeness in partial coverage regions. The analysis shows the differences in the phase of the corresponding seismic interferometry (by cross‐correlation) and virtual reflector (by cross‐convolution) signals obtained by 2D and 3D formulations, with synthetic examples performed to remove water layer multiples in ocean bottom seismic (OBS) acoustic data.     

Modelling the wave phenomena in acoustic and elastic media with sharp variations of physical properties using the grid‐characteristic method

This paper introduces a novel method of modelling acoustic and elastic wave propagation in inhomogeneous media with sharp variations of physical properties based on the recently developed grid‐characteristic method which considers different types of waves generated in inhomogeneous linear‐elastic media (e.g., longitudinal, transverse, Stoneley, Rayleigh, scattered PP‐, SS‐waves, and converted PS‐ and SP‐waves). In the framework of this method, the problem of solving acoustic or elastic wave equations is reduced to the interpolation of the solutions, determined at earlier time, thus avoiding a direct solution of the large systems of linear equations required by the FD or FE methods. We apply the grid‐characteristic method to compare wave phenomena computed using the acoustic and elastic wave equations in geological medium containing a hydrocarbon reservoir or a fracture zone. The results of this study demonstrate that the developed algorithm can be used as an effective technique for modelling wave phenomena in the models containing hydrocarbon reservoir and/or the fracture zones, which are important targets of seismic exploration.     

Three-dimensional Marchenko internal multiple attenuation on narrow azimuth streamer data of the Santos Basin,Brazil

In recent years, a variety of Marchenko methods for the attenuation of internal multiples has been developed. These methods have been extensively tested on two-dimensional synthetic data and applied to two-dimensional field data, but only little is known about their behaviour on three-dimensional synthetic data and three-dimensional field data. Particularly, it is not known whether Marchenko methods are sufficiently robust for sparse acquisition geometries that are found in practice. Therefore, we start by performing a series of synthetic tests to identify the key acquisition parameters and limitations that affect the result of three-dimensional Marchenko internal multiple prediction and subtraction using an adaptive double-focusing method. Based on these tests, we define an interpolation strategy and use it for the field data application. Starting from a wide azimuth dense grid of sources and receivers, a series of decimation tests are performed until a narrow azimuth streamer geometry remains. We evaluate the effect of the removal of sail lines, near offsets, far offsets and outer cables on the result of the adaptive double-focusing method. These tests show that our method is most sensitive to the limited aperture in the crossline direction and the sail line spacing when applying it to synthetic narrow azimuth streamer data. The sail line spacing can be interpolated, but the aperture in the crossline direction is a limitation of the acquisition. Next, we apply the adaptive Marchenko double-focusing method to the narrow azimuth streamer field data from the Santos Basin, Brazil. Internal multiples are predicted and adaptively subtracted, thereby improving the geological interpretation of the target area. These results imply that our adaptive double-focusing method is sufficiently robust for the application to three-dimensional field data, although the key acquisition parameters and limitations will naturally differ in other geological settings and for other types of acquisition.     

Comparison of acoustic and elastic full-waveform inversion of 2D towed-streamer data in the presence of salt

Over the last years, full-waveform inversion has become an important tool in the list of processing and imaging technologies available to the industry. For marine towed-streamer data, full-waveform inversion is typically applied using an acoustic approximation because S-waves do not propagate in water and elastic effects in recorded data are generally assumed to be small. We compare acoustic and elastic modelling and full-waveform inversion for a field data set acquired offshore Angola over sediments containing a salt body with significant topology. Forward modelling tests reveal that such geological structures lead to significant mode conversions at interfaces and, consequently, to significant relative amplitude differences when elastically and acoustically modelled traces are compared. Using an acoustic approach for modelling in full-waveform inversion therefore leads to problems matching the synthetic data with the field data, even for recorded pressure data and with trace normalization applied. Full-waveform inversion is unable to find consistent model updates. Applying elastic full-waveform inversion leads to more consistent and reliable model updates with less artefacts, at the expense of additional computation cost. Although two-dimensional marine towed-streamer data are least favourable for the application of full-waveform inversion compared to three-dimensional data or ocean-bottom data, it is recommended to check on the existence of elastic effects before deciding on the final processing and imaging approach.     

Research note: Seismic attenuation due to wave‐induced fluid flow at microscopic and mesoscopic scales

Wave‐induced fluid flow at microscopic and mesoscopic scales arguably constitutes the major cause of intrinsic seismic attenuation throughout the exploration seismic and sonic frequency ranges. The quantitative analysis of these phenomena is, however, complicated by the fact that the governing physical processes may be dependent. The reason for this is that the presence of microscopic heterogeneities, such as micro‐cracks or broken grain contacts, causes the stiffness of the so‐called modified dry frame to be complex‐valued and frequency‐dependent, which in turn may affect the viscoelastic behaviour in response to fluid flow at mesoscopic scales. In this work, we propose a simple but effective procedure to estimate the seismic attenuation and velocity dispersion behaviour associated with wave‐induced fluid flow due to both microscopic and mesoscopic heterogeneities and discuss the results obtained for a range of pertinent scenarios.     

Estimating turbulent dissipation rates from acoustic backscatter

Recent work has suggested the possibility of using high-frequency (10 kHz to 1 MHz) acoustics to remotely sense turbulent microstructure in the ocean. Once developed, this technique will provide a powerful tool for examining the spatial and temporal distributions of ocean mixing in a way not possible using traditional microstructure methods. In this paper turbulent dissipation rates are estimated through an inversion of 307 kHz acoustic scattering data collected in the lee of a sill in Knight Inlet, British Columbia. These data have been shown previously to be strongly correlated with temperature and shear microstructure measurements. Here we show that inversion methods can be used to get turbulent dissipation rates from acoustic backscatter, provided that independent measurements of temperature and salinity stratification are available. The temperature–salinity characteristics of the environment, however, can place limitations on the inversion technique. The strong negative salinity gradient in Knight Inlet decreases the slope of the functional relationship between dissipation rate and scattering cross-section for high dissipation rates and increases the uncertainty of the inversion. This limitation on the inversion technique is not an issue throughout most of the world's oceans (where dT/dS>0) and, in places where it could be a problem, the limitation can be overcome by using multi-frequency techniques.     

Analytical approximations of bulk and shear moduli for dry rock based on the differential effective medium theory

Differential effective medium theory has been applied to determine the elastic properties of porous media. The ordinary differential equations for bulk and shear moduli are coupled and it is more difficult to obtain accurate analytical formulae about the moduli of dry porous rock. In this paper, in order to decouple these equations we first substitute an analytical approximation for the dry‐rock modulus ratio into the differential equation and derive analytical solutions of the bulk and shear moduli for dry rock with three specific pore shapes: spherical pores, needle‐shaped pores and penny‐shaped cracks. Then, the validity of the analytical approximations is tested by integrating the full differential effective medium equation numerically. The analytical formulae give good estimates of the numerical results over the whole porosity range for the cases of the three given pore shapes. These analytical formulae can be further simplified under the assumption of small porosity. The simplified formulae for spherical pores are the same as Mackenzie's equations. The analytical formulae are relatively easy to analyse the relationship between the elastic moduli and porosity or pore shapes and can be used to invert some rock parameters such as porosity or pore aspect ratio. The predictions of the analytical formulae for experimental data show that the formulae for penny‐shaped cracks are suitable to estimate the elastic properties of micro‐crack rock such as granite, they can be used to estimate the crack aspect ratio while the crack porosity is known and also to estimate the crack porosity evolution with pressure if the crack aspect ratio is given.     

Deterring effects of 8-45 kHz tone pulses on harbour seals (Phoca vitulina) in a large pool

The marine aquaculture industry suffers losses due to pinniped attacks which damage net enclosures and fish stocks. Acoustic harassment devices (AHDs) emit loud sounds which are intended to deter pinnipeds from approaching aquaculture enclosures. At present, many AHDs emit sounds in the 8-20 kHz frequency range. It is not known whether sounds of higher frequencies have a deterrent effect on seals. Therefore five captive harbour seals (Phoca vitulina) were subjected to four series of tone pulses together spanning a broad frequency range (8, 16, 32 and 45 kHz). Pulse duration was 250 ms and pulse interval was 5s. Each of the four sounds was made deterrent by increasing the amplitude. The seals reacted by swimming away from the sounds. The displacement effect of each sound was judged by comparing the animals' surface positions, and number of surfacings, during ten 45 min baseline periods with ten 45 min test periods per frequency (one frequency per day in rotation, 40 sessions in total). The seals were displaced by all four frequencies throughout the 40 trial days. The seals came to the surface more often when the test tones were produced than in the baseline periods. The initial displacement distances did not change over the 40 test days. This suggests that operating AHDs for only short periods will be more effective and less likely to result in habituation by the seals than operating them continuously. The discomfort threshold sound pressure level (SPL) was established for each of the four pulse frequencies. The acoustic discomfort threshold SPL is defined as the boundary SPL between the area that the animals generally occupied during the transmission of the sounds and the area that they generally did not enter during sound transmission. The discomfort threshold SPL may depend on the context.