SEARCH · Engineering Papers

Results for “geostatistics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2023

DEEPEN 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

DEEPEN: Final 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

Part of the DEEPEN (DE-risking Exploration of geothermal Plays in magmatic ENvironments) project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology, detailed in the journal article below, is based on the method proposed by Poux & O'brien (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., (2022). For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty.

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗

DOE OSTI2021

Assessing the Consistency of Estimated Ground Cover Fractions between the BLM AIM Method and Optical Remote Sensing Method

Since 2012, Argonne National Laboratory (Argonne) has supported the Bureau of Land Management (BLM) in developing remote sensing methodologies for long-term environmental monitoring of Palo Verde Mesa in eastern Riverside County, California, including methods for: detailed mapping of ephemeral streams, estimating fractional cover of desert-land surface components (e.g., trees, shrubs, litters, and bare ground), evaluating erosion risk or land stability, and characterizing vegetation alliances using spatial structure and geostatistical approaches. These studies showed the promise of remote sensing for monitoring changes in desert landscapes by providing information that would be difficult to obtain through field surveys. During this time the BLM has also worked to establish long-term monitoring protocols and compiled field-observation data collected using standardized protocols from the Assessment, Inventory, and Monitoring (AIM) strategy. The AIM data can be compared to data derived from remote sensing methods to evaluate their relative operational utility in monitoring landscape change. If ground cover estimated using remotely sensed imagery is comparable to AIM ground cover estimates, then remote sensing can be used to monitor whether any land cover change in desert landscapes may be related to solar energy development. Therefore, the goal of this study was to determine the consistency in ground cover estimates between AIM data and those derived from publicly- available remotely sensed imagery, such as that available through the U.S. Department of Agriculture, National Agricultural Imagery Program (NAIP), to examine feasibility of a remote sensing method for complementing AIM monitoring. Based on the image analysis in this study, we also provide recommendations for how small unmanned aerial system (sUAS) data may be used to complement BLM’s AIM data and NAIP imagery for future vegetation monitoring. The ground cover types we originally planned to investigate were trees, shrubs, and bare ground. However, the small sample size and a limited range of cover fraction of trees and shrubs in the AIM dataset (e.g., 33 samples with a maximum shrub cover of 18%, 14 samples with a maximum tree cover of 17%) did not allow for performing a meaningful evaluation for the remote sensing approach. Therefore, we conducted the study focusing on bare ground, foliar, and rock cover, all of which are indicators reported in the AIM remote sensing dataset.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

Applying 3D Geologic Modeling Workflows to the Argillite Reference Case (Rev. 1)

The objective of this short report is to document the application of our 3D geologic modeling workflow to an argillite (shale) host rock. Over the past four years, our team at Los Alamos National Laboratory has developed a geologic modeling workflow that can be applied to generic alluvial basins such as those found in the western United States. In “frontier” or “exploratory” basins where data are sparse, the first steps are to collect, evaluate and integrate available subsurface data into conceptual geologic models. Those models form the basis for constructing the geologic framework model, a 3D geocellular model ideally constrained by seismic and borehole data. To date we have constructed our models using “synthetic” well data derived from conceptual models, without the prospect of validating our workflow using “real” subsurface data. We were tasked to investigate whether our workflow designed for alluvial basin sediments could be applied to other potential repository host rocks. This task also provided the opportunity to work with high-quality subsurface data collected specifically for siting and evaluating a nuclear waste repository. Nagra, the Swiss governmental agency responsible for the disposal of the nation’s radioactive waste, generously provided us with data from two deep boreholes drilled through their argillaceous target formation. The aim of our proof-of-concept demonstration is to evaluate whether geostatistical methods offer a viable approach to property modeling in argillaceous rocks. Nagra provided us with the well data on the condition that we maintain confidentiality with all transferred information and results. Fortunately, Nagra posts numerous technical reports on its public website that describe the subsurface geology in great detail. All of the information and illustrations in this report related to the Swiss repository enterprise are taken from the Nagra public website.

58 GEOSCIENCES↗

DOE OSTI2024

Generative network-based approaches to generate stochastic realizations

Generative Adversarial Network (GAN) – based models have been successfully applied in generating different geological models in the literature. However, it is still challenging to use GAN to generate geological realizations with extremely sparse conditioning data (e.g. several well data), which may be regarded as local noise by GAN during the training process. In this work, we propose a novel conditional Generative Adversarial Neural Operator (cGANO) to tackle this challenge. In cGANO, the mapping between conditioning data and output is established through the U-shaped neural operators (UNO), which better preserves local information. Another advantage of using UNO comes from its grid-independent property, which makes the generation of downscaling stochastic geologic realizations possible. We tested the model performance on the IBDP geostatistical dataset with 100 realizations.

58 GEOSCIENCES↗

DOE OSTI2022

Alternative Conceptual Models of the Subsurface at the Hanford Site

This report describes the use of a multipoint geostatistics (MPS) framework with multiple data types for developing alternative conceptual-mathematical models of the subsurface at the Hanford Site. Developing multiple conceptual-mathematical models of the site is important for assessing potential uncertainties related to the features, events, and processes that control the subsurface flow of groundwater and the transport of contaminants of concern (COCs). Evaluating the effects of these uncertainties can result in better site characterization, monitoring, and remedial decisions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Initial Uncertainty Analysis of Carbon Tetrachloride Contamination and Remediation in the Ringold A and Lower Mud Units at the Central Plateau

The long-term effectiveness of groundwater cleanup at the Hanford Site Central Plateau depends on predictive models that can capture key uncertainties in contaminant fate and transport. Carbon tetrachloride (CCl 4 ), a persistent and toxic compound, presents particular challenges due to variability in degradation rates, uncertainty in initial plume distribution, and subsurface heterogeneity. These uncertainties directly influence plume persistence, migration pathways, and remedy performance, and thus must be systematically evaluated to support long-term remediation planning. To address these gaps, a large-scale Monte Carlo analysis was conducted using the Plateau to River (P2R) model framework. The modeling approach parameterized three primary uncertainty factors: (1) degradation rate, (2) initial plume distribution, and (3) hydraulic conductivity. Degradation was represented as a first-order process, with half-lives ranging from 70 to 700 years. Initial plume distributions were created using a geostatistical simulation method (sgsim), which generates many equally plausible versions of how contaminants might be distributed underground. From this, 100 different scenarios were mapped onto the P2R grid. Variability in hydraulic conductivity was represented in a similar way, with 100 scenarios each for the Ringold Lower Mud and Ringold A units (layers 6 and 7), based on fitted exponential variograms and conditioned to well data. In total, more than 1000 realizations were simulated to assess plume behavior under uncertainty. Results demonstrate that degradation kinetics exert the strongest control over plume persistence: Shorter half-lives produced rapid mass reduction, while longer half-lives yielded persistent plumes with limited attenuation. A nonlinear response was observed, with steep mass reductions at half-lives greater than 200 years and near-linear declines beyond this threshold, reflecting interactions between degradation and pumping. The initial plume distribution strongly influenced early transport patterns, with broader sources generating larger plume footprints, although pump-and-treat operations constrained plume migration to managed areas. By comparison, hydraulic conductivity variability in the Ringold units had only a secondary influence, modifying spreading behavior without altering the dominant migration pathways governed by source configuration and hydraulic controls. Overall, the analysis highlights that uncertainty in degradation rate and initial plume configuration are the primary drivers of variability in plume predictions, while conductivity heterogeneity plays a limited role. These findings underscore the need for improved site-specific data on degradation processes and source characterization to enhance the reliability of long-term performance assessments and to better inform remedial decision-making at the Central Plateau.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1988

Estimating the signal-to-noise ratio of AVIRIS data

To make the best use of narrowband airborne visible/infrared imaging spectrometer (AVIRIS) data, an investigator needs to know the ratio of signal to random variability or noise (signal-to-noise ratio or SNR). The signal is land cover dependent and varies with both wavelength and atmospheric absorption; random noise comprises sensor noise and intrapixel variability (i.e., variability within a pixel). The three existing methods for estimating the SNR are inadequate, since typical laboratory methods inflate while dark current and image methods deflate the SNR. A new procedure is proposed called the geostatistical method. It is based on the removal of periodic noise by notch filtering in the frequency domain and the isolation of sensor noise and intrapixel variability using the semi-variogram. This procedure was applied easily and successfully to five sets of AVIRIS data from the 1987 flying season and could be applied to remotely sensed data from broadband sensors.

Curran, Paul J.↗

NASA NTRS1988

Zones of information in the AVIRIS spectra

To make the best use of Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data an investigator needs to know the ratio of signal to random variability or noise (S/N ratio). The signal is land-cover dependent and decreases with both wavelength and atmospheric absorption and random noise comprises sensor noise and intra-pixel variability. The three existing methods for estimating the S/N ratio are inadequate as typical laboratory methods inflate, while dark current and image methods deflate the S/N ratio. We propose a new procedure called the geostatistical method. It is based on the removal of periodic noise by notch filtering in the frequency domain and the isolation of sensor noise and intra-pixel variability using the semi-variogram. This procedure was applied easily and successfully to five sets of AVIRIS data from the 1987 flying season.

Curran, Paul J.↗

NASA NTRS1989

Autocorrelation and regularization in digital images. II - Simple image models

The variogram function used in geostatistical analysis is a useful statistic in the analysis of remotely sensed images. Using the results derived by Jupp et al. (1988), the basic second-order, or covariance, properties of scenes modeled by simple disks of varying size and spacing after imaging into disk-shaped pixels are analyzed to explore the relationship betwee image variograms and discrete object scene structure. The models provide insight into the nature of real images of the earth's surface and the tools for a complete analysis of the more complex case of three-dimensional illuminated discrete-object images.

Jupp, David L. B.↗

NASA NTRS1989

Estimation of signal-to-noise - A new procedure applied to AVIRIS data

Curran, Paul J.↗

NASA NTRS1991

Estimating The SNR Of AVIRIS Data

Combination of established statistical and filtering techniques yields improved estimates of signal-to-noise ratio (SNR) of multispectral imagery produced by Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). Method yields more accurate estimate of random variability in AVIRIS imagery. Combination of techniques called "geostatistical" method. In first step of method, periodic noise removed by notch filtering in frequency domain. In subsequent steps, sensor noise and variability within picture elements isolated by use of semi-variance of picture elements.

Curran, Paul J.↗

NASA NTRS1989

Decorrelation distance of snow in the Colorado River Basin

The problem of estimating areal averages from point measurement has been extensively studied by mining engineers and hydrologists. Its application to satellite measurements has recently been introduced. The semivariaogram has been used in many geostatistical applications to estimate spatial structures of observed properties, such as mineral distributions. An examination is made of snow variations in Colorado from daily snow data collected in 11 SNOTEL stations. The associated semivariogram is estimated. The objective is to estimate the spatial structure of the snow field so that the point data can be used for comparison with, and validation for, satellite measurements.

Chang, A. T. C.↗

NASA NTRS1991

Environmental analysis using integrated GIS and remotely sensed data - Some research needs and priorities

This paper discusses some basic scientific issues and research needs in the joint processing of remotely sensed and GIS data for environmental analysis. Two general topics are treated in detail: (1) scale dependence of geographic data and the analysis of multiscale remotely sensed and GIS data, and (2) data transformations and information flow during data processing. The discussion of scale dependence focuses on the theory and applications of spatial autocorrelation, geostatistics, and fractals for characterizing and modeling spatial variation. Data transformations during processing are described within the larger framework of geographical analysis, encompassing sampling, cartography, remote sensing, and GIS. Development of better user interfaces between image processing, GIS, database management, and statistical software is needed to expedite research on these and other impediments to integrated analysis of remotely sensed and GIS data.

Davis, Frank W.↗

NASA NTRS1993

Functional patterns in an annual grassland during an AVIRIS overflight

This study relates Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) imagery to ground measurements of vegetation distribution, physiology, and productivity at Stanford University's Jasper Ridge Biological Preserve. Primary efforts focused on a 9-ha region of annual grassland where we completed a detailed ground-based study in conjunction with a 15 May 1991 AVIRIS overflight. Spectral mixture analysis and the normalized difference vegetation index (NDVI) calculated from AVIRIS data were used to evaluate spatial patterns of vegetation type, productivity, and potential physiological activity. Concurrent ground sampling revealed a high degree of correlation between NDVI and estimates of canopy chemistry, structure, productivity, and CO2 flux, supporting the use of imaging spectrometry to estimate spatial and temporal trends in vegetation physiology and productivity in this relatively simple grassland ecosystem. Geostatistical analyses of both ground and AVIRIS data supported the conclusion that the AVIRIS pixel size was suitable for describing the influence of major landscape features in this grassland and that spatial detail would be lost at slightly larger pixel sizes typical of other imaging spectrometers.

Gamon, John A.↗

NASA NTRS1996

New results from mapping Antarctica at high resolution from radar altimeter data

The limits of the altimeter data evaluation were analyzed using geostatistical methods. It is often argued that satellite radar altimeter data over ice can not be used to map ice surfaces with a slope exceeding 0.5 deg. The maps presented, obtained from the Geosat satellite geodetic mission, concern the Antarctica north of 72.1 deg South. The grids constructed from other satellite observations facilitate the analysis of the changes in the Antarctic ice stream/ice shelf systems. The evaluation was applied to the Lambert Glacier/Amery Ice Shelf system.

Matassa, Michael↗

NASA NTRS1997

Evaluation of GIS Technology in Assessing and Modeling Land Management Practices

There is an increasing concern of land owners to protect and maintain healthy and sustainable agroecosystems through the implementation of best management practices (BMP). The objectives of this study were: (1) To develop and evaluate the use of a Geographic Information System (GIS) technology for enhancing field-scale management practices; (2) evaluate the use of 2-dimensional displays of the landscape and (3) define spatial classes of variables from interpretation of geostatistical parameters. Soil samples were collected to a depth of 2 m at 15 cm increments. Existing data from topographic, land use, and soil survey maps of the Winfred Thomas Agricultural Research Station were converted to digital format. Additional soils data which included texture, pH, and organic matter were also generated. The digitized parameters were used to create a multilayered field-scale GIS. Two dimensional (2-D) displays of the parameters were generated using the ARC/INFO software. The spatial distribution of the parameters evaluated in both fields were similar which could be attributed to the similarity in vegetation and surface elevation. The ratio of the nugget to total semivariance, expressed as a percentage, was used to assess the degree of spatial variability. The results indicated that most of the parameters were moderate spatially dependent Biophysical constraint maps were generated from the database layers, and used in multiple combination to visualize results of the BMP. Understanding the spatial relationships of physical and chemical parameters that exists within a field should enable land managers to more effectively implement BMP to ensure a safe and sustainable environment.

Archer, F.↗