A rock-physics investigation of unconsolidated saline permafrost: P-wave properties from laboratory ultrasonic measurements

Geophysics ◽  
2016 ◽  
Vol 81 (1) ◽  
pp. WA233-WA245 ◽  
Author(s):  
Shan Dou ◽  
Seiji Nakagawa ◽  
Douglas Dreger ◽  
Jonathan Ajo-Franklin

Saline permafrost is sensitive to thermal disturbances and is prone to subsidence, which renders it a major source of geohazard in Arctic coastal environments. Seismic methods could be used to map and monitor saline permafrost at scales of geotechnical interests because of the ice-content dependencies of seismic properties. We have developed a comprehensive study of the ultrasonic P-wave properties (i.e., velocity and attenuation) of synthetic saline permafrost samples for a range of salinities and temperatures, and measurements conducted on a fine-grained permafrost core obtained from Barrow, Alaska. The resulting data consist of P-wave properties presented as functions of temperature and salinity. Notable observations include the following: P-wave velocities showed marked reductions in the presence of dissolved salts and complex variations resulting from the water-to-ice phase transitions; strong P-wave attenuation was present in the temperature intervals in which the samples were partially frozen. When presented as functions of ice saturation, the data sets lead us to two key findings: (1) neither a purely cementing nor a purely pore-filling model of the pore-scale distributions of ice could adequately fit the observed velocity data and (2) although the velocities increase monotonically with increasing ice saturations, P-wave attenuation reaches a maximum at intermediate ice saturations—contrary to the ordinary expectation of decreasing attenuation with increasing velocities. The observed ice-content dependencies of P-wave properties, along with the implications on the probable pore-scale distributions of ice, provide a valuable basis for rock-physics modeling, which in turn could facilitate seismic characterizations of saline permafrost.

Geophysics ◽  
2020 ◽  
pp. 1-26
Author(s):  
Xiaomin Zhao ◽  
Mark E. Willis ◽  
Tanya Inks ◽  
Glenn A. Wilson

Several recent studies have advanced the use of time-lapse distributed acoustic sensing (DAS) vertical seismic profile (VSP) data in horizontal wells for determining hydraulically stimulated fracture properties. Hydraulic fracturing in a horizontal well typically generates vertical fractures in the rock medium around each stage. We model the hydraulically stimulated formation with vertical fracture sets about the lateral wellbore as a horizontally transverse isotropic (HTI) medium. Rock physics modeling is used to relate the anisotropy parameters to fracture properties. This modeling was used to develop an inversion for P-wave time delay to fracture height and density of each stage. Field data from two horizontal wells were analyzed, and fracture height evaluated using this technique agreed with microseismic analysis.


Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. D625-D641 ◽  
Author(s):  
Dario Grana

The estimation of rock and fluid properties from seismic attributes is an inverse problem. Rock-physics modeling provides physical relations to link elastic and petrophysical variables. Most of these models are nonlinear; therefore, the inversion generally requires complex iterative optimization algorithms to estimate the reservoir model of petrophysical properties. We have developed a new approach based on the linearization of the rock-physics forward model using first-order Taylor series approximations. The mathematical method adopted for the inversion is the Bayesian approach previously applied successfully to amplitude variation with offset linearized inversion. We developed the analytical formulation of the linearized rock-physics relations for three different models: empirical, granular media, and inclusion models, and we derived the formulation of the Bayesian rock-physics inversion under Gaussian assumptions for the prior distribution of the model. The application of the inversion to real data sets delivers accurate results. The main advantage of this method is the small computational cost due to the analytical solution given by the linearization and the Bayesian Gaussian approach.


Geophysics ◽  
2015 ◽  
Vol 80 (2) ◽  
pp. E83-E95
Author(s):  
Richard T. Houck ◽  
Adrian Ciucivara ◽  
Scott Hornbostel

Unconstrained 3D inversion of marine controlled source electromagnetic data (CSEM) data sets produces resistivity volumes that have an uncertain relationship to the true subsurface resistivity at the scale of typical hydrocarbon reservoirs. Furthermore, CSEM-scale resistivity is an ambiguous indicator of hydrocarbon presence; not all resistivity anomalies are caused by hydrocarbon reservoirs, and not all hydrocarbon reservoirs produce a distinct resistivity anomaly. We have developed a method for quantifying the effectiveness of resistivities from CSEM inversion in detecting hydrocarbon reservoirs. Our approach uses probabilistic rock-physics modeling to update information from a preexisting prospect assessment, based on uncertain resistivities from CSEM. The result is an estimate the probability of hydrocarbon presence that accounts for uncertainty in the resistivity and in rock properties. Examples using synthetic and real CSEM data sets demonstrate that the effectiveness of CSEM inversion in identifying hydrocarbon reservoirs depends on the interaction between the uncertainty associated with the inversion-derived resistivity and the range of rock and fluid properties that were expected for the targeted prospect. Resistivity uncertainty that has a small effect on hydrocarbon probability for one set of rock property distributions may have a large effect for a different set of rock properties. Depending on the consequences of this interaction, resistivities from CSEM inversion might reduce the risk associated with predictions of hydrocarbon presence, but they cannot be expected to guarantee a specific well outcome.


Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. C177-C191 ◽  
Author(s):  
Yunyue Li ◽  
Biondo Biondi ◽  
Robert Clapp ◽  
Dave Nichols

Seismic anisotropy plays an important role in structural imaging and lithologic interpretation. However, anisotropic model building is a challenging underdetermined inverse problem. It is well-understood that single component pressure wave seismic data recorded on the upper surface are insufficient to resolve a unique solution for velocity and anisotropy parameters. To overcome the limitations of seismic data, we have developed an integrated model building scheme based on Bayesian inference to consider seismic data, geologic information, and rock-physics knowledge simultaneously. We have performed the prestack seismic inversion using wave-equation migration velocity analysis (WEMVA) for vertical transverse isotropic (VTI) models. This image-space method enabled automatic geologic interpretation. We have integrated the geologic information as spatial model correlations, applied on each parameter individually. We integrate the rock-physics information as lithologic model correlations, bringing additional information, so that the parameters weakly constrained by seismic are updated as well as the strongly constrained parameters. The constraints provided by the additional information help the inversion converge faster, mitigate the ambiguities among the parameters, and yield VTI models that were consistent with the underlying geologic and lithologic assumptions. We have developed the theoretical framework for the proposed integrated WEMVA for VTI models and determined the added information contained in the regularization terms, especially the rock-physics constraints.


Sign in / Sign up

Export Citation Format

Share Document