1,721,004 research outputs found
Comparison of two- and three-dimensional full waveform inversion imaging using wide-angle seismic data from the Deep Galicia Margin
Full text linkFull waveform inversion (FWI) is a data-fitting technique capable of generating high-resolution velocity models with a resolution down to half the seismic wavelength. FWI is applied typically to densely sampled seismic data. In this study, we applied FWI to 3D wide-angle seismic data acquired using sparsely spaced ocean bottom seismometers (OBSs) from the Deep Galicia Margin west of Iberia. Our dataset samples the S-reflector, a low-angle detachment present in this area. Here we highlight differences between 2D, 2.5D and 3D-FWI performances using a real sparsely spaced dataset. We performed 3D FWI in the time domain and compared the results with 2D and 2.5D FWI results from a profile through the 3D model. When overlaid on multichannel seismic images, the 3D FWI results constrain better the complex faulting within the pre- and syn-rift sediments and crystalline crust compared to the 2D result. Furthermore, we estimate variable serpentinisation of the upper mantle below the S-reflector along the profile using 3D FWI, reaching a maximum of 45 per cent. Differences in the data residuals of the 2D, 2.5D and 3D inversions suggest that 2D inversion can be prone to overfitting when using a sparse dataset. To validate our results, we performed tests to recover the anomalies introduced by the inversions in the final models using synthetic datasets. Based on our comparison of the velocity models, we conclude that the use of 3D data can partially mitigate the problem of receiver sparsity in FWI
Integrated geophysical characterization of crustal domains in the eastern Black Sea
Full text linkRifting may lead ultimately to continental breakup, but the identification and characterization of the resulting crustal distribution remains challenging. Also, spatial and temporal changes in breakup magmatism may affect the geophysical character of the newly formed oceanic crust, resulting in contrasting interpretations of crustal composition and distribution. In the Eastern Black Sea Basin (EBSB), the evolution from rifting to breakup has been long debated, with several interpretations for the distribution of stretched continental and oceanic crust. We interpret basement morphological variations from long-offset seismic reflection profiles, highlighting a northwest-to-southeast transition from faulted and tilted continental blocks to a rough and then smoother basement. We model magnetic anomalies to further constrain the various basement domains, and infer the presence of a weakly magnetized, stretched continental crust in the northwest, and a 0.4-3.8 A/m layer coinciding with the smooth basement in the central and southeastern area. We conclude that the EBSB oceanic crust extends farther to the northwest than was suggested previously from an abrupt change in crustal thickness and lower-crustal velocity. The apparent discrepancy between these different types of geophysical evidence may result from changes in magma supply during breakup, affecting the thickness and velocity structure of the resulting oceanic crust.</p
The effect of heterogeneities in hydrate saturation on gas production from natural systems
Full text linkUnderstanding the rate and time evolution of gas release from natural gas hydrate systems is important when evaluating the potential of gas hydrate as a future energy source, or the impact of gas from hydrate on climate. The release of gas from hydrate is heavily influenced by a number of factors, many of which vary through the hydrate system. The fundamental heterogeneity of natural gas hydrate systems is often poorly represented in models. Here we simulate depressurisation-induced gas production from a single vertical well in 34 models with heterogeneous 2D distributions of hydrate that include layered, columnar or random configurations and comparable models with homogenous saturation distributions. We found that the temporal evolution of gas production rate follows a consistent trend for all models, but at any time the gas production rate across the models varied by up to ±35% in the first year of production, and by up to ±25% thereafter. The primary control on the gas production rate is the overall amount of hydrate in the system, but local variations in hydrate saturation cause significant fluctuations in the time evolution of production. These hydrate variations can cause changes in the gas flow path through the system and associated drops in gas production rate continuing for multiple years. Overall, our results suggest that small levels of heterogeneity in hydrate systems can cause variations in the gas production rate similar in scale to much larger variations in homogenous systems. Our work provides an error margin for previously modelled gas production rates, and a note of caution for potential commercial development of gas hydrate
Seismic characterisation of multiple BSRs in the Eastern Black Sea Basin
Full text linkLong offset seismic reflection data reveal the presence of four Bottom Simulating Reflectors (BSR0-3) within folded sediments of the Tuapse Trough, along the NE margin of the Eastern Black Sea Basin (EBSB). Multiple BSRs are observed in other sites worldwide, however, their origin and formation mechanisms are still debated. Here, we investigate the formation mechanisms of the EBSB multiple BSRs based on their seismic character and on their physical properties derived from reflected and refracted arrival seismic velocities. Seismic reflection data are downward continued to enhance refracted arrivals. A 2D travel-time velocity model of the sub-seabed, using combined travel-times from non-downward-continued reflected and downward-continued refracted signals, shows variations in the physical properties at the BSRs and nearby sediments. The P-wave velocity (VP) increase of 1.55–1.72 km/s between the seafloor and BSR0 (258 mbsf) reflects normal compaction trends in sediments, whereas the VP of 1.75–1.83 km/s between BSR0 and BSR1 (360 mbsf) is higher than that expected for sediments at that depth. Beneath BSR1, a VP decrease from 1.83 km/s to 1.61 km/s occurs within a 70-80 m-thick layer including BSR2 (395 mbsf) and extending to BSR3 (438 mbsf). Beneath BSR3, VP increases. Based on an analytical model linking seismic velocity to physical properties, these VP trends can be explained by a gas hydrate saturation from 0 to 2% between the seafloor and BSR0, reaching 4 ± 2% just above BSR1. A free gas saturation of up to 20–25% is estimated within the low-velocity zone between BSR1 and BSR3. BSR1 likely represents the present-day base of the gas hydrate stability zone (BGHSZ), which aligns with the theoretical BGHSZ assuming a geothermal gradient of 26–30 °C/km. Based on seismic polarities and results from travel-time analysis and rock physics modelling, we suggest that hydrate dissociation and recycling processes may explain the negative polarity of BSR2 and BSR3, which are still visible due to the presence of relict gas, and inferred higher gas hydrate saturations close to the present-day base of the stability zone at BSR1. Also, structural and stratigraphic controls seem to have favoured focused free gas flow and hydrate formation at the top of an anticlinal structure, thus likely controlling multiple BSR generation in the EBSB
A social, environmental and economic evaluation protocol for potential gas hydrate exploitation projects
Full text linkThere is increasing global interest in the potential commercial development of methane gas hydrate as a widespread and abundant unconventional source of natural gas. Previous work has focussed on understanding the nature and distribution of the resource, and potential recovery technology, neglecting assessment of the associated social, economic and environmental consequences. This gap needs to be addressed for any commercial gas hydrate development business case to succeed. Here we develop a multi-criteria decision analysis (MCDA) protocol of gas hydrate development using the ELECTRE III method. Our protocol proposes criteria that evaluate the social, environmental and economic impacts of gas hydrate development proposals, which are weighted to represent the priorities of six identified stakeholder groups. We have tested the protocol on potential commercial gas hydrate development in Alaska through a series of interviews. Our results show that there is no universal preference structure, even within stakeholder groups, indicating that buy-in from all groups is a complex compromise. However, there are two fundamentally opposing groups, one composed of individuals from governmental and industry backgrounds who prioritise economic criteria, and another represented by members of the local community and environmental advocates who prioritise social and environmental criteria. The protocol concludes that gas hydrate development in Alaska is unlikely to be supported under present-day conditions. This work provides the first structured foundation for comprehensive assessment of future development proposals of gas hydrate or other natural resources.</p
Laboratory observations of frequency-dependent ultrasonic P-wave velocity and attenuation during methane hydrate formation in Berea sandstone
Full text linkKnowledge of the effect of methane hydrate saturation and morphology on elastic wave attenuation could help reduce ambiguity in seafloor hydrate content estimates. These are needed for seafloor resource and geohazard assessment, as well as to improve predictions of greenhouse gas fluxes into the water column. At low hydrate saturations, measuring attenuation can be particularly useful as the seismic velocity of hydrate-bearing sediments is relatively insensitive to hydrate content. Here, we present laboratory ultrasonic (448–782 kHz) measurements of P-wave velocity and attenuation for successive cycles of methane hydrate formation (maximum hydrate saturation of 26 per cent) in Berea sandstone. We observed systematic and repeatable changes in the velocity and attenuation frequency spectra with hydrate saturation. Attenuation generally increases with hydrate saturation, and with measurement frequency at hydrate saturations below 6 per cent. For hydrate saturations greater than 6 per cent, attenuation decreases with frequency. The results support earlier experimental observations of frequency-dependent attenuation peaks at specific hydrate saturations. We used an effective medium rock-physics model which considers attenuation from gas bubble resonance, inertial fluid flow and squirt flow from both fluid inclusions in hydrate and different aspect ratio pores created during hydrate formation. Using this model, we linked the measured attenuation spectral changes to a decrease in coexisting methane gas bubble radius, and creation of different aspect ratio pores during hydrate formation
An anisotropic model for the electrical resistivity of two-phase geologic materials
Full text linkElectrical and electromagnetic surveys of the seafloor provide valuable information about the macro and microscopic properties of subseafloor sediments. Sediment resistivity is highly variable and governed by a wide range of properties including pore-fluid salinity, pore-fluid saturation, porosity, pore geometry, and temperature. A new anisotropic, twophase, effective medium model describes the electrical resistivity of porous rocks and sediments. The only input parameters required are the resistivities of the solid and fluid components, their volume fractions and grain shape. The approach makes use of the increase in path length taken by an electrical current through an idealized granular medium comprising of aligned ellipsoidal grains. The model permits both solid and fluid phases to have a finite conductivity useful for dealing with surface charge conduction effects associated with clay minerals and gives results independent of grain size hence, valid for a wide range of sediment types. Furthermore, the model can be used to investigate the effects of grain aspect ratio and alignment on electrical resistivity anisotropy. Good agreement was found between the model predictions and laboratory measurements of resistivity and porosity on artificial sediments with known physical properties
Joint elastic-electrical effective medium models of reservoir sandstones
No full textImprovements in the joint inversion of seismic and marine controlled source electromagnetic data sets will require better constrained models of the joint elastic-electrical properties of reservoir rocks. Various effective medium models were compared to a novel laboratory data set of elastic velocity and electrical resistivity (obtained on 67 reservoir sandstone samples saturated with 35 g/l brine at a differential pressure of 8 MPa) with mixed results. Hence, we developed a new three-phase effective medium model for sandstones with pore-filling clay minerals based on the combined self-consistent approximation and differential effective medium model. We found that using a critical porosity of 0.5 and an aspect ratio of 1 for all three components, the proposed model gave accurate model predictions of the observed magnitudes of P-wave velocity and electrical resistivity and of the divergent trends of clean and clay-rich sandstones at higher porosities. Using only a few well-constrained input parameters, the new model offers a practical way to predict?in situ?porosity and clay content in brine saturated sandstones from co-located P-wave velocity and electrical resistivity data sets
An effective medium inversion algorithm for gas hydrate quantification and its application to laboratory and borehole measurements of gas hydrate-bearing sediments
Full text linkThe presence of gas hydrate in marine sediments alters their physical properties. In some circumstances, gas hydrate may cement sediment grains together and dramatically increase the seismic P- and S-wave velocities of the composite medium. Hydrate may also form a load-bearing structure within the sediment microstructure, but with different seismic wave attenuation characteristics, changing the attenuation behaviour of the composite. Here we introduce an inversion algorithm based on effective medium modelling to infer hydrate saturations from velocity and attenuation measurements on hydrate-bearing sediments. The velocity increase is modelled as extra binding developed by gas hydrate that strengthens the sediment microstructure. The attenuation increase is modelled through a difference in fluid flow properties caused by different permeabilities in the sediment and hydrate microstructures. We relate velocity and attenuation increases in hydrate-bearing sediments to their hydrate content, using an effective medium inversion algorithm based on the self-consistent approximation (SCA), differential effective medium (DEM) theory, and Biot and squirt flow mechanisms of fluid flow. The inversion algorithm is able to convert observations in compressional and shear wave velocities and attenuations to hydrate saturation in the sediment pore space. We applied our algorithm to a data set from the Mallik 2L–38 well, Mackenzie delta, Canada, and to data from laboratory measurements on gas-rich and water-saturated sand samples. Predictions using our algorithm match the borehole data and water-saturated laboratory data if the proportion of hydrate contributing to the load-bearing structure increases with hydrate saturation. The predictions match the gas-rich laboratory data if that proportion decreases with hydrate saturation. We attribute this difference to differences in hydrate formation mechanisms between the two environments
- …
