SEARCH · Geology Reports

Results for “Inverse Problems”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7Linked to original sources

USGS2012

Approaches in highly parameterized inversion - PEST++, a Parameter ESTimation code optimized for large environmental models

An object-oriented parameter estimation code was developed to incorporate benefits of object-oriented programming techniques for solving large parameter estimation modeling problems. The code is written in C++ and is a formulation and expansion of the algorithms included in PEST, a widely used parameter estimation code written in Fortran. The new code is called PEST++ and is designed to lower the barriers of entry for users and developers while providing efficient algorithms that can accommodate large, highly parameterized problems. This effort has focused on (1) implementing the most popular features of PEST in a fashion that is easy for novice or experienced modelers to use and (2) creating a software design that is easy to extend; that is, this effort provides a documented object-oriented framework designed from the ground up to be modular and extensible. In addition, all PEST++ source code and its associated libraries, as well as the general run manager source code, have been integrated in the Microsoft Visual Studio® 2010 integrated development environment. The PEST++ code is designed to provide a foundation for an open-source development environment capable of producing robust and efficient parameter estimation tools for the environmental modeling community into the future.

Techniques and Methods↗

USGS2008

Geophysical characterization of the American River levees, Sacramento, California, with electromagnetics, capacitively coupled resistivity, and DC resistivity

A geophysical characterization of a portion of American River levees in Sacramento, California was conducted in May, 2007. Targets of interest included the distribution and thickness of sand lenses that underlie the levees and the depth to a clay unit that underlies the sand. The concern is that the erosion of these sand lenses can lead to levee failure in highly populated areas of Sacramento. DC resistivity (Geometric’s OhmMapper and Advanced Geosciences, Inc.’s SuperSting R8 systems) and electromagnetic surveys (Geophex’s GEM-2) were conducted over a 6 mile length of the levee on roads and bicycle and horse trails. 2-D inversions were conducted on all the geophysical data. The OhmMapper and SuperSting surveys produced consistent inversion results that delineated potential sand and clay units. GEM-2 apparent resistivity data were consistent with the DC inversion results. However, the GEM-2 data could not be inverted due to low electromagnetic response levels, high ambient electromagnetic noise, and large system drifts. While this would not be as large a problem in conductive terrains, it is a problem for a small induction number electromagnetic profiling system such as the GEM-2 in a resistive terrain (the sand lenses). An integrated interpretation of the geophysical data acquired in this investigation is presented in this report that includes delineation of those areas consisting of predominantly sand and those areas consisting predominantly of clay. In general, along most of this part of the American River levee system, sand lenses are located closest to the river and clay deposits are located further away from the river. The interpreted thicknesses of the detected sand deposits are variable and range from 10 ft up to 60 ft. Thus, despite issues with the GEM-2 inversion, this geophysical investigation successfully delineated sand lenses and clay deposits along the American River levee system and the approximate depths to underlying clay zones. The results of this geophysical investigation should help the USACE to maintain the current levee system while also assisting the designers and planners of levee enhancements with the knowledge of what is to be expected from the near-surface geology and where zones of concern may be located.

California↗

USGS2024

The importance of nodal plane orientation diversity for earthquake focal mechanism stress inversions

Inversions of earthquake focal mechanisms are among the most accessible and reliable methods for determining crustal stress. However, the use of this method varies widely, and assumptions that underpin it are often violated, potentially compromising stress estimates. We investigate the consequences of violating the little-studied assumption that the focal mechanisms have diverse orientations. Our approach is to employ data-informed synthetic mechanisms, with nodal plane orientations defined by recent earthquake lineaments in the Midland Basin, western Texas, and rakes consistent with slip in the mapped stress field. Using both the traditional stress inversion method that assumes constant shear stress magnitudes on the causative faults as well as a recently published variable shear stress method, we show that low fault plane diversity can cause maximum horizontal stress ( S Hmax ) orientation and relative principal stress magnitude (faulting regime) estimates to differ markedly from the true values. This problem is compounded for catalogs with even modest amounts of noise (≤15°) or few (e.g., 20) mechanisms. Significantly, traditional approaches for quantifying uncertainty such as the bootstrap can severely underestimate the true uncertainty under these circumstances. To remedy this, we provide simple tools to quantify nodal plane orientation diversity and stress inversion reliability.

Geological Society of London Special Publications↗

USGS1998

Singularity and Nonnormality in the Classification of Compositional Data

Geologists may want to classify compositional data and express the classification as a map. Regionalized classification is a tool that can be used for this purpose, but it incorporates discriminant analysis, which requires the computation and inversion of a covariance matrix. Covariance matrices of compositional data always will be singular (noninvertible) because of the unit-sum constraint. Fortunately, discriminant analyses can be calculated using a pseudo-inverse of the singular covariance matrix; this is done automatically by some statistical packages such as SAS. Granulometric data from the Darss Sill region of the Baltic Sea is used to explore how the pseudo-inversion procedure influences discriminant analysis results, comparing the algorithm used by SAS to the more conventional Moore-Penrose algorithm. Logratio transforms have been recommended to overcome problems associated with analysis of compositional data, including singularity. A regionalized classification of the Darss Sill data after logratio transformation is different only slightly from one based on raw granulometric data, suggesting that closure problems do not influence severely regionalized classification of compositional data.

Mathematical Geology↗

USGS2013

Effects of error covariance structure on estimation of model averaging weights and predictive performance

When conducting model averaging for assessing groundwater conceptual model uncertainty, the averaging weights are often evaluated using model selection criteria such as AIC, AICc, BIC, and KIC (Akaike Information Criterion, Corrected Akaike Information Criterion, Bayesian Information Criterion, and Kashyap Information Criterion, respectively). However, this method often leads to an unrealistic situation in which the best model receives overwhelmingly large averaging weight (close to 100%), which cannot be justified by available data and knowledge. It was found in this study that this problem was caused by using the covariance matrix, C E , of measurement errors for estimating the negative log likelihood function common to all the model selection criteria. This problem can be resolved by using the covariance matrix, C ek , of total errors (including model errors and measurement errors) to account for the correlation between the total errors. An iterative two-stage method was developed in the context of maximum likelihood inverse modeling to iteratively infer the unknown C ek from the residuals during model calibration. The inferred C ek was then used in the evaluation of model selection criteria and model averaging weights. While this method was limited to serial data using time series techniques in this study, it can be extended to spatial data using geostatistical techniques. The method was first evaluated in a synthetic study and then applied to an experimental study, in which alternative surface complexation models were developed to simulate column experiments of uranium reactive transport. It was found that the total errors of the alternative models were temporally correlated due to the model errors. The iterative two-stage method using C ek resolved the problem that the best model receives 100% model averaging weight, and the resulting model averaging weights were supported by the calibration results and physical understanding of the alternative models. Using C ek obtained from the iterative two-stage method also improved predictive performance of the individual models and model averaging in both synthetic and experimental studies.

Water Resources Research↗

USGS2010

Time-dependent seismic tomography

Of methods for measuring temporal changes in seismic-wave speeds in the Earth, seismic tomography is among those that offer the highest spatial resolution. 3-D tomographic methods are commonly applied in this context by inverting seismic wave arrival time data sets from different epochs independently and assuming that differences in the derived structures represent real temporal variations. This assumption is dangerous because the results of independent inversions would differ even if the structure in the Earth did not change, due to observational errors and differences in the seismic ray distributions. The latter effect may be especially severe when data sets include earthquake swarms or aftershock sequences, and may produce the appearance of correlation between structural changes and seismicity when the wave speeds are actually temporally invariant. A better approach, which makes it possible to assess what changes are truly required by the data, is to invert multiple data sets simultaneously, minimizing the difference between models for different epochs as well as the rms arrival-time residuals. This problem leads, in the case of two epochs, to a system of normal equations whose order is twice as great as for a single epoch. The direct solution of this system would require twice as much memory and four times as much computational effort as would independent inversions. We present an algorithm, tomo4d, that takes advantage of the structure and sparseness of the system to obtain the solution with essentially no more effort than independent inversions require. No claim to original US government works Journal compilation ?? 2010 RAS.

Geophysical Journal International↗

USGS1979

Monitoring massive fracture growth at 2-km depths using surface tiltmeter arrays

Tilt due to massive hydraulic fractures induced in sedimentary rocks at depths of up to 2.2 km have been recorded by surface tiltmeters. Injection of fluid volumes up to 4 · 10 5 liters and masses of propping agent up to 5 · 10 5 kg is designed to produce fractures approximately 1 km long, 50–100 m high and about 1 cm wide. The surface tilt data adequately fit a dislocation model of a tensional fault in a half-space. Theoretical and observational results indicate that maximum tilt occurs at a distance off the strike of the fracture equivalent to 0.4 of the depth to the fracture. Azimuth and extent of the fracture deduced from the geometry of the tilt field agree with other kinds of geophysical measurements. Detailed correlation of the tilt signatures with pumping parameters (pressure, rate, volume, mass) have provided details on asymmetry in geometry and growth rate. Whereas amplitude variations in tilt vary inversely with the square of the depth, changes in flow rate or pressure gradient can produce a cubic change in width. These studies offer a large-scale experimental approach to the study of problems involving fracturing, mass transport, and dilatancy processes.

Tectonophysics↗

USGS2005

Constraints on fault slip rates of the southern California plate boundary from GPS velocity and stress inversions

We use Global Positioning System (GPS) velocities and stress orientations inferred from seismicity to invert for the distribution of slip on faults in the southern California plate-boundary region. Of particular interest is how long-term slip rates are partitioned between the Indio segment of the San Andreas fault (SAF), the San Jacinto fault (SJF) and the San Bernardino segment of the SAE We use two new sets of constraints to address this problem. The first is geodetic velocities from the Southern California Earthquake Center's (SCEC) Crustal Motion Map (version 3 by Shen et al.), which includes significantly more data than previous models. The second is a regional model of stress-field orientations at seismogenic depths, as determined from earthquake focal mechanisms. While GPS data have been used in similar studies before, this is the first application of stress-field observations to this problem. We construct a simplified model of the southern California fault system, and estimate the interseismic surface velocities using a backslip approach with purely elastic strain accumulation, following Meade et al. In addition, we model the stress orientations at seismogenic depths, assuming that crustal stress results from the loading of active faults. The geodetically derived stressing rates are found to be aligned with the stress orientations from seismicity. We therefore proceed to invert simultaneously GPS and stress observations for slip rates of the faults in our network. We find that the regional patterns of crustal deformation as imaged by both data sets can be explained by our model, and that joint inversions lead to better constrained slip rates. In our preferred model, the SJF accommodates ???15 mm yr-1 and the Indio segment of the SAF ???23 mm yr-1 of right-lateral motion, accompanied by a low slip rate on the San Bernardino segment of the SAF 'Anomalous' fault segments such as around the 1992 Mw = 7.3 Landers surface rupture can be detected. There, observed stresses deviate strongly from the long-term loading as predicted by our simple model. Evaluation of model misfits together with information from palaeoseismology may provide further insights into the time dependence of strain accumulation along the San Andreas system. ?? 2004 RAS.

Geophysical Journal International↗

USGS1982

Earthquake location in island arcs

A comprehensive data set of selected teleseismic P-wave arrivals and local-network P- and S-wave arrivals from large earthquakes occurring at all depths within a small section of the central Aleutians is used to examine the general problem of earthquake location in island arcs. Reference hypocenters for this special data set are determined for shallow earthquakes from local-network data and for deep earthquakes from combined local and teleseismic data by joint inversion for structure and location. The high-velocity lithospheric slab beneath the central Aleutians may displace hypocenters that are located using spherically symmetric Earth models; the amount of displacement depends on the position of the earthquakes with respect to the slab and on whether local or teleseismic data are used to locate the earthquakes. Hypocenters for trench and intermediate-depth events appear to be minimally biased by the effects of slab structure on rays to teleseismic stations. However, locations of intermediate-depth events based on only local data are systematically displaced southwards, the magnitude of the displacement being proportional to depth. Shallow-focus events along the main thrust zone, although well located using only local-network data, are severely shifted northwards and deeper, with displacements as large as 50 km, by slab effects on teleseismic travel times. Hypocenters determined by a method that utilizes seismic ray tracing through a three-dimensional velocity model of the subduction zone, derived by thermal modeling, are compared to results obtained by the method of joint hypocenter determination (JHD) that formally assumes a laterally homogeneous velocity model over the source region and treats all raypath anomalies as constant station corrections to the travel-time curve. The ray-tracing method has the theoretical advantage that it accounts for variations in travel-time anomalies within a group of events distributed over a sizable region of a dipping, high-velocity lithospheric slab. In application, JHD has the practical advantage that it does not require the specification of a theoretical velocity model for the slab. Considering earthquakes within a 260 km long by 60 km wide section of the Aleutian main thrust zone, our results suggest that the theoretical velocity structure of the slab is presently not sufficiently well known that accurate locations can be obtained independently of locally recorded data. Using a locally recorded earthquake as a calibration event, JHD gave excellent results over the entire section of the main thrust zone here studied, without showing a strong effect that might be attributed to spatially varying source-station anomalies. We also calibrated the ray-tracing method using locally recorded data and obtained results generally similar to those obtained by JHD. ?? 1982.

Physics of the Earth and Planetary Interiors↗

USGS2011

Mechanical analysis of fault slip data: Implications for paleostress analysis

Stress inversions are a useful and popular tool for structural geologist and seismologist alike. These methods were first introduced by Wallace (1951) and Bott (1959) and subsequent studies continue to be based on their assumptions: the remote stress tensor is spatially uniform for the rock mass containing the faults and temporally constant over the history of faulting in that region, and the slip on each fault surface has the same direction and sense as the maximum shear stress resolved on that surface from the remote stress tensor. Furthermore, successful implementation requires that slip accumulates on faults of diverse orientation. Many studies employ these methods on isolated faults or on fault systems with limited ranges of orientations, which can lead to erroneous results. We propose a new method that incorporates the effects of mechanical interaction of the entire fault or fault system, and solves the complete mechanical problem rather than employing empirical relationships between slip and stress or strain (or strain rate). The method requires knowledge of the fault geometry and information on at least one slip vector component along portions of the known fault geometry. For example, if throw is known, the strike-slip component can be solved for. We test the method using a single synthetic fault with anisotropic roughness similar to that measured at fault outcrops. While the orientation of remote stress may be determined precisely, the lack of diverse fault orientations introduces a systematic error in the remote stress ratio. We further test the effect of diversity of fault orientations and find that Wallace–Bott type inversions do not perform as well for limited ranges of orientations when compared to the proposed method. Finally, we use published data from the 1999 Chi-Chi, Taiwan, earthquake, and find that the method using surface data only, and surface data with subsurface focal mechanisms, produce similar results. The resulting stress orientations are in good agreement with results from Wallace–Bott inversions. Furthermore, the slip distribution is in general agreement with kinematic slip inversions using coseismic surface deformation. Stress inversion methods using fault slip data can thus be improved upon, significantly in some cases, by solving a mechanical boundary value problem that takes into account the geometry of faults or fault systems. As a bonus, the solution provides the stress, strain, and displacement fields throughout the region and the slip distributions on the faults.

Journal of Structural Geology↗

USGS2009

Obtaining parsimonious hydraulic conductivity fields using head and transport observations: A Bayesian geostatistical parameter estimation approach

Flow path delineation is a valuable tool for interpreting the subsurface hydrogeochemical environment. Different types of data, such as groundwater flow and transport, inform different aspects of hydrogeologic parameter values (hydraulic conductivity in this case) which, in turn, determine flow paths. This work combines flow and transport information to estimate a unified set of hydrogeologic parameters using the Bayesian geostatistical inverse approach. Parameter flexibility is allowed by using a highly parameterized approach with the level of complexity informed by the data. Despite the effort to adhere to the ideal of minimal a priori structure imposed on the problem, extreme contrasts in parameters can result in the need to censor correlation across hydrostratigraphic bounding surfaces. These partitions segregate parameters into facies associations. With an iterative approach in which partitions are based on inspection of initial estimates, flow path interpretation is progressively refined through the inclusion of more types of data. Head observations, stable oxygen isotopes ( 18 O/ 16 O ratios), and tritium are all used to progressively refine flow path delineation on an isthmus between two lakes in the Trout Lake watershed, northern Wisconsin, United States. Despite allowing significant parameter freedom by estimating many distributed parameter values, a smooth field is obtained.

Water Resources Research↗

USGS2012

Reflectance spectroscopy: Quantitative analysis techniques for remote sensing applications

Several methods for the analysis of remotely sensed reflectance data are compared, including empirical methods and scattering theories, both of which are important for solving remote sensing problems. The concept of the photon mean optical path length and the implications for use in modeling reflectance spectra are presented. It is shown that the mean optical path length in a particulate surface is in rough inverse proportion to the square root of the absorption coefficient. Thus, the stronger absorber a material is, the less photons will penetrate into the surface. The concept of apparent absorbance (‐In reflectance) is presented, and it is shown that absorption bands, which are Gaussian in shape when plotted as absorption coefficient (true absorbance) versus photon energy, are also Gaussians in apparent absorbance. However, the Gaussians in apparent absorbance have a smaller intensity and a width which is a factor of √2 larger. An apparent continuum in a reflectance spectrum is modeled as a mathematical function used to isolate a particular absorption feature for analysis. It is shown that a continuum should be removed by dividing it into the reflectance spectrum or subtracting it from the apparent absorbance and that the fitting of Gaussians to absorption features should be done using apparent absorbance versus photon energy. Kubelka‐Munk theory is only valid for materials with small total absorption and for bihemispherical reflectance, which are rarely encountered in geologic remote sensing. It is shown that the recently advocated bidirectional reflectance theories have the potential for use in deriving mineral abundance from a reflectance spectrum.

Journal of Geophysical Research↗

USGS2013

Assessing grain-size correspondence between flow and deposits of controlled floods in the Colorado River, USA

Flood-deposited sediment has been used to decipher environmental parameters such as variability in watershed sediment supply, paleoflood hydrology, and channel morphology. It is not well known, however, how accurately the deposits reflect sedimentary processes within the flow, and hence what sampling intensity is needed to decipher records of recent or long-past conditions. We examine these problems using deposits from dam-regulated floods in the Colorado River corridor through Marble Canyon–Grand Canyon, Arizona, U.S.A., in which steady-peaked floods represent a simple end-member case. For these simple floods, most deposits show inverse grading that reflects coarsening suspended sediment (a result of fine-sediment-supply limitation), but there is enough eddy-scale variability that some profiles show normal grading that did not reflect grain-size evolution in the flow as a whole. To infer systemwide grain-size evolution in modern or ancient depositional systems requires sampling enough deposit profiles that the standard error of the mean of grain-size-change measurements becomes small relative to the magnitude of observed changes. For simple, steady-peaked floods, 5–10 profiles or fewer may suffice to characterize grain-size trends robustly, but many more samples may be needed from deposits with greater variability in their grain-size evolution.

Journal of Sedimentary Research↗

USGS2015

Inverse modeling with RZWQM2 to predict water quality

This chapter presents guidelines for autocalibration of the Root Zone Water Quality Model (RZWQM2) by inverse modeling using PEST parameter estimation software ( Doherty, 2010 ). Two sites with diverse climate and management were considered for simulation of N losses by leaching and in drain flow: an almond [ Prunus dulcis (Mill.) D.A. Webb] orchard in the San Joaquin Valley, California and the Walnut Creek watershed in central Iowa, which is predominantly in corn ( Zea mays L.)–soybean [ Glycine max (L.) Merr.] rotation. Inverse modeling provides an objective statistical basis for calibration that involves simultaneous adjustment of model parameters and yields parameter confidence intervals and sensitivities. We describe operation of PEST in both parameter estimation and predictive analysis modes. The goal of parameter estimation is to identify a unique set of parameters that minimize a weighted least squares objective function, and the goal of predictive analysis is to construct a nonlinear confidence interval for a prediction of interest by finding a set of parameters that maximizes or minimizes the prediction while maintaining the model in a calibrated state. We also describe PEST utilities (PAR2PAR, TSPROC) for maintaining ordered relations among model parameters (e.g., soil root growth factor) and for post-processing of RZWQM2 outputs representing different cropping practices at the Iowa site. Inverse modeling provided reasonable fits to observed water and N fluxes and directly benefitted the modeling through: (i) simultaneous adjustment of multiple parameters versus one-at-a-time adjustment in manual approaches; (ii) clear indication by convergence criteria of when calibration is complete; (iii) straightforward detection of nonunique and insensitive parameters, which can affect the stability of PEST and RZWQM2; and (iv) generation of confidence intervals for uncertainty analysis of parameters and model predictions. Composite scaled sensitivities, which reflect the total information provided by the observations for a parameter, indicated that most of the RZWQM2 parameters at the California study site (CA) and Iowa study site (IA) could be reliably estimated by regression. Correlations obtained in the CA case indicated that all model parameters could be uniquely estimated by inverse modeling. Although water content at field capacity was highly correlated with bulk density (−0.94), the correlation is less than the threshold for nonuniqueness (0.95, absolute value basis). Additionally, we used truncated singular value decomposition (SVD) at CA to mitigate potential problems with highly correlated and insensitive parameters. Singular value decomposition estimates linear combinations (eigenvectors) of the original process-model parameters. Parameter confidence intervals (CIs) at CA indicated that parameters were reliably estimated with the possible exception of an organic pool transfer coefficient (R45), which had a comparatively wide CI. However, the 95% confidence interval for R45 (0.03–0.35) is mostly within the range of values reported for this parameter. Predictive analysis at CA generated confidence intervals that were compared with independently measured annual water flux (groundwater recharge) and median nitrate concentration in a collocated monitoring well as part of model evaluation. Both the observed recharge (42.3 cm yr −1 ) and nitrate concentration (24.3 mg L −1 ) were within their respective 90% confidence intervals, indicating that overall model error was within acceptable limits.

Book chapter↗

USGS2007

Physically based estimation of soil water retention from textural data: General framework, new models, and streamlined existing models

Numerous models are in widespread use for the estimation of soil water retention from more easily measured textural data. Improved models are needed for better prediction and wider applicability. We developed a basic framework from which new and existing models can be derived to facilitate improvements. Starting from the assumption that every particle has a characteristic dimension R associated uniquely with a matric pressure ψ and that the form of the ψ– R relation is the defining characteristic of each model, this framework leads to particular models by specification of geometric relationships between pores and particles. Typical assumptions are that particles are spheres, pores are cylinders with volume equal to the associated particle volume times the void ratio, and that the capillary inverse proportionality between radius and matric pressure is valid. Examples include fixed-pore-shape and fixed-pore-length models. We also developed alternative versions of the model of Arya and Paris that eliminate its interval-size dependence and other problems. The alternative models are calculable by direct application of algebraic formulas rather than manipulation of data tables and intermediate results, and they easily combine with other models (e.g., incorporating structural effects) that are formulated on a continuous basis. Additionally, we developed a family of models based on the same pore geometry as the widely used unsaturated hydraulic conductivity model of Mualem. Predictions of measurements for different suitable media show that some of the models provide consistently good results and can be chosen based on ease of calculations and other factors.

Vadose Zone Journal↗

USGS2010

Double-porosity models for a fissured groundwater reservoir with fracture skin

Theories of flow to a well in a double-porosity groundwater reservoir are modified to incorporate effects of a thin layer of low-permeability material or fracture skin that may be present at fracture-block interfaces as a result of mineral deposition or alteration. The commonly used theory for flow in double- porosity formations that is based upon the assumption of pseudo–steady state block-to-fissure flow is shown to be a special case of the theory presented in this paper. The latter is based on the assumption of transient block-to-fissure flow with fracture skin. Under conditions where fracture skin has a hydraulic conductivity that is less than that of the matrix rock, it may be assumed to impede the interchange of fluid between the fissures and blocks. Resistance to flow at fracture-block interfaces tends to reduce spatial variation of hydraulic head gradients within the blocks. This provides theoretical justification for neglecting the divergence of flow in the blocks as required by the pseudo–steady state flow model. Coupled boundary value problems for flow to a well discharging at a constant rate were solved in the Laplace domain. Both slab-shaped and sphere-shaped blocks were considered, as were effects of well bore storage and well bore skin. Results obtained by numerical inversion were used to construct dimensionless-type curves that were applied to well test data, for a pumped well and for an observation well, from the fractured volcanic rock terrane of the Nevada Test Site.

Water Resources Research↗

USGS2018

Gas hydrate quantification using full-waveform inversion of sparse ocean-bottom seismic data: A case study from Green Canyon Block 955, Gulf of Mexico

We present a case study of gas hydrate quantification using dense short-offset multichannel seismic (MCS) and sparse long-offset ocean-bottom-seismometer (OBS) data in lease block Green Canyon 955 (GC955), Gulf of Mexico (GOM), where the presence of gas hydrate was interpreted using logging while drilling (LWD) data acquired by the GOM Gas Hydrate Joint Industry Project Leg II expedition. We use frequency-domain full-waveform inversion (FWI) of seven OBS gathers to invert for a P-wave velocity model of an approximately 7 km long MCS profile connecting two LWD sites, GC955-H and GC955-Q. We build the starting model for FWI using traveltime inversion (TI) of the MCS and OBS data. In addition, we use the TI model for depth migrating the MCS stack. At the LWD sites, we constrain the hydrate saturation (Sgh) using sonic and resistivity logs. Unfortunately, as is typical of seismic quantification problems, the FWI model resolution is not sufficient to extrapolate the LWD-based Sgh. Therefore, we apply Backus averaging to the sonic log, at 60 m wavelength, bringing it within approximately 8% of the FWI model and make the assumption that averaging the sonic log is same as redistributing the gas hydrate within the Backus wavelength. In this manner, instead of Sgh, the FWI model is able to estimate the total gas hydrate volume. In the end, we use the FWI model and the migrated stack to constrain the locations and bulk volumes of free gas and gas hydrate. Our results demonstrate that with careful processing, reasonable estimates on locations and bulk volumes of submarine gas hydrate accumulations can be achieved even with sparse seismic data that are not adequate for amplitude-based assessments.

Geophysics↗

USGS2023

Using DC resistivity ring array surveys to resolve conductive structures around tunnels or mine-workings

In underground environments, conventional direct current (DC) resistivity surveys with a single linear array of electrodes produce fundamentally non-unique inversions. These non-uniqueness and model resolution issues stem from limitations placed on the location of transmitters (TXs) and receivers (RXs) by the geometry of existing tunnels and boreholes. Poor excitation and/or sampling of the region of interest (ROI) can create artifacts and reduce the resolution of the recovered model. To address these problems we propose the use of an ensemble of ring arrays, which are created by placing one or more electrodes in each face (sidewalls, floor, and ceiling) of the tunnel to form a ring of electrodes at each along-tunnel location. Using a series of increasingly complex synthetic models, we assess the benefits of ring arrays and show that they can be used to better constrain the location and shape of anomalous bodies around the tunnel. Although ring arrays significantly improve the resolution of the recovered model, the size of the comprehensive ring array survey increases rapidly with the number of electrodes used. To balance model resolution and survey size, we developed a physics-based survey design methodology. In this methodology, TXs are selected based upon secondary charge accumulations on a test block that is moved throughout the ROI. Although this survey design methodology does not produce a strictly optimal survey, it balances model resolution and survey size in a practical and computationally efficient manner. Since the ring array more accurately estimates the around-tunnel location of targets and ensures that targets on all sides of the tunnel are detected, it is ideally suited to tunnel-based environments. Our results show that only about 6 % "> 6% of the possible TXs and 0.5 % "> 0.5% of the RXs in the comprehensive ring array survey are needed to retain the improvements in resolution. Therefore, economical ring array surveys can be designed for both reconnaissance and target characterization. Following the inversion of the reconnaissance dataset, additional rings can be added to reduce the inter-ring spacing or off-tunnel boreholes can be added to the region around identified anomalies to increase resolution as required.

Journal of Applied Geophysics↗