Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Well Log”

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.

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

At least 451 records · Page 25

High pressure/high temperature multiphase simulations of dodecane injection to nitrogen: Application on ECN Spray-A

The present work investigates the complex phenomena associated with pressure/high temperature dodecane injection for the Engine Combustion Network (ECN) Spray-A case, employing more elaborate thermodynamic closures, to avoid well known deficiencies concerning density and speed of sound prediction using traditional cubic models. A tabulated thermodynamic approach is proposed here, based on log 10 (p)-T tables, providing very high accuracy across a large range of pressures, spanning from 0 to 2500 bar, with only a small number of interpolation points. The tabulation approach is directly extensible to any thermodynamic model, existing or to be developed in the future. Here NIST REFPROP properties are used, combined with PC-SAFT Vapor-Liquid-Equilibrium to identify the liquid in mixtures penetration, hence avoiding the use of an arbitrary threshold for mass fraction. Identified liquid and vapour penetration are compared against experimental data from the ECN database showing a good agreement, within approximately 3–8% for axial penetration of liquid, 2% for vapor axial penetration and within experimental uncertainty for radial distribution of mass fraction. Analysis of the vortex evolution indicates that driving mechanisms behind the jet break-up are vortex tilting/stretching, then baroclinic torque, leading to Rayleigh-Taylor instabilities, closely followed by vortex dilation and finally viscous effects.

42 ENGINEERING↗

Efficient Estimation of Pauli Observables by Derandomization

In this work, we consider the problem of jointly estimating expectation values of many Pauli observables, a crucial subroutine in variational quantum algorithms. Starting with randomized measurements, we propose an efficient derandomization procedure that iteratively replaces random single-qubit measurements by fixed Pauli measurements; the resulting deterministic measurement procedure is guaranteed to perform at least as well as the randomized one. In particular, for estimating any L low-weight Pauli observables, a deterministic measurement on only of order log(L) copies of a quantum state suffices. In some cases, for example, when some of the Pauli observables have high weight, the derandomized procedure is substantially better than the randomized one. Specifically, numerical experiments highlight the advantages of our derandomized protocol over various previous methods for estimating the ground-state energies of small molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Constraining AGN Torus Sizes with Optical and Mid-infrared Ensemble Structure Functions

We propose a new method to constrain the size of the dusty torus in broad-line active galactic nuclei (AGNs) using optical and mid-infrared (MIR) ensemble structure functions (SFs). Because of the geometric dilution of the torus, the MIR response to optical continuum variations has suppressed variability with respect to the optical that depends on the geometry (e.g., size, orientation, opening angle) of the torus. More extended tori have steeper MIR SFs with respect to the optical SFs. We demonstrate the feasibility of this SF approach using simulated AGN light curves and a geometric torus model. While it is difficult to use SFs to constrain the orientation and opening angle due to the insensitivity of the SF on these parameters, the size of the torus can be well determined. Applying this method to the ensemble SFs measured for 587 SDSS quasars, we measure a torus R–L relation of $\mathrm{log}\,{R}_{\mathrm{eff}}(\mathrm{pc})={0.51}_{-0.04}^{+0.04}\times \mathrm{log}({{L}}_{\mathrm{bol}}/{10}^{46}\,\mathrm{erg}\ {{\rm{s}}}^{-1})-{0.38}_{-0.01}^{+0.01}$ in the WISE W1 band and sizes ~1.4 times larger in the W2 band, which are in good agreement with dust reverberation mapping measurements. Compared with the reverberation mapping technique, the SF method is much less demanding in data quality and can be applied to any optical+MIR light curves for which a lag measurement may not be possible, as long as the variability process and torus structure are stationary. While this SF method does not extract all information contained in the light curves (i.e., the transfer function), it provides an intuitive interpretation for the observed trends of AGN MIR SFs compared with optical SFs.

79 ASTRONOMY AND ASTROPHYSICS↗

Far-ultraviolet energy distributions of the metal-poor A stars HD 109995 and HD 161817

Low-resolution IUE spectra at wavelengths between 1300 and 3400 A of the metal-poor stars HD 109995 (A1p) and HD 161817 (A4p) have been compared with model-atmosphere energy distributions computed by Kurucz (1979). Good overall agreement is found. Effective temperatures, metal abundances, and angular diameters could be determined. Assuming an absolute visual magnitude of 0.7, the previously determined gravity log = 3 yields masses of 0.5 solar masses for both stars. It is found that the theoretical UBV colors calculated earlier agree reaonably well with the ones observed for these stars.

Boehm-Vitense, E.↗

ROSAT observations of the Chamaeleon star forming cloud

Two-soft X-ray images of the Chamaeleon I star forming cloud obtained with the ROSAT Position Sensitive Proportional Counter (PSPC) are presented. Seventy reliable, and perhaps 19 additional, X-ray sources are found. Up to Ninety percent of these sources are certainly or probably identified with T Tauri stars formed in the cloud. Twenty to 35 are probably previously unrecognized 'weak' T Tauri (WTT) stars. T Tauri X-ray luminosities range from log bar-L(sub x) less than 28.8 to 31.1 erg/s, or 10(exp 2) - 10(exp 4) times solar levels, with mean value of log bar-L(sub x) = 10(exp 29.2) erg/s in the 0.2-2.5 keV band. The X-ray luminosities of well-studied Chamaeleon cloud members are correlated with a complex of four stellar properties: effective temperature, mass, radius and bolometric luminosity. The spatial distribution, Hertzsprung-Russell (H-R) diagram locations of the stars indicate WTT and CTT are coeval. The total premain sequence population of the cloud is likely to be greater than 100 stars, with WTT stars outnumbering 'classical' T Tauri (CTT) stars by approximately 2:1.

Feigelson, E. D.↗

Trapped Melt in IIIAB Irons: Solid/Liquid Elemental Partitioning During the Fractionation of the IIIAB Magma

Group IIIAB, the largest iron-meteorite group, shows compositional trends (including a three-order-of-magnitude It concentration range) indicating that it formed by fractional crystallization of a metallic magma. Because about 200 irons are available, and all degrees of crystallization are well represented, IIIAB offers an excellent set of samples for the study of crystallization at all depths of the asteroidal core. On log-log Ir-Au, and Ir-As diagrams IIIAB forms a broad band; the breadth represents real meteorite-to-meteorite variations, far outside experimental or sampling uncertainties. A successful model must explain the width of this band; I suggest that it mainly resulted from the trapping of parental magma within the crystallizing solid. Because S is essentially insoluble in metal, the abundance of FeS is a measure of the fraction of trapped liquid. The trapped-melt model is supported by the observation that irons having higher S contents plot closer to the inferred composition of the magmatic parental liquid. The lowest S values are found in the irons occupying the left envelope of the IIIAB Ir-Au or Ir-As compositional fields, thus it is this set of irons that should be interpreted as the solid products of a fractionating magma. This simplifies the modeling of the crystallization process and allows inferences regarding the distribution ratios for other elements in the evolved IIIAB system. The large (multiton) Cape York irons show wide variations in their trapped-melt fractions; their compositions seem best understood in terms of a low initial S content of the IIIAB magma, about 20 mg/g. The inferred initial IIIAB distribution coefficient for Ir, 4.6, is much higher than published values based on laboratory studies of low-S systems; I suggest that low-S (and low-P) partition-ratio measurements tend to err in the direction of unity. In IIIAB distribution coefficients for Au, As, and Ni were still < 1 when the most evolved IIIAB irons formed, another indication of a low initial S content.

Wasson, John T.↗

Transcriptomics-based Machine Learning (ML) Analysis Predicts Space-Exposed Murine Livers

NASA has employed high-throughput molecular assays to identify sub-cellular changes impacting human physiology during spaceflight. Machine learning (ML) methods hold the promise to improve our ability to identify important signals within highly dimensional molecular data. However, the inherent limitation of study subject numbers within a spaceflight mission minimizes the utility of ML approaches. To overcome the sample power limitations, data from multiple spaceflight missions must be aggregated while appropriately addressing intra- and inter-study variabilities. Here we describe an approach to log transform, scale and normalize data from six heterogeneous, mouse liver derived transcriptomics datasets (ntotal=137) which enabled ML-methods to perform well (AUC ≥ 0.87) in classifying spaceflown vs ground control animals rather than mission-of-origin. Concordance was found between liver-specific biological processes identified from harmonized ML-based analysis and study-by-study classical omics analysis. This work demonstrates the feasibility of applying ML methods on integrated, heterogeneous datasets of small sample size.

Machine Learning↗

Utah FORGE 2-2439v2: Characterizing In-Situ Stress with Laboratory Modelling and Field Measurements - 2024 Annual Workshop Presentation

This is a presentation on A Multi-Component Approach to Characterizing In-Situ Stress at the Utah FORGE Site: Laboratory Modelling and Field Measurements project by The University of Pittsburgh, presented by Andrew Bunger. The project characterizes the stress in the Utah FORGE EGS reservoir using three methods: Method 1: Demonstrate complimentary laboratory rock-core stress estimation combined with Machine Learning approach for measuring in-situ stress from field sonic log data; Method 2: Complete field based in-situ measurement (mini-frac); and Method 3: Develop a mechanics-based method for connection near wellbore stress measurements to stresses away from the well-bore. This presentation was featured in the Utah FORGE R&D Annual Workshop on August 14, 2024.

15 GEOTHERMAL ENERGY↗

TESS Delivers Five New Hot Giant Planets Orbiting Bright Stars from the Full-frame Images

We present the discovery and characterization of five hot and warm Jupiters—TOI-628 b (TIC 281408474; HD 288842), TOI-640 b (TIC 147977348), TOI-1333 b (TIC 395171208, BD+47 3521A), TOI-1478 b (TIC 409794137), and TOI-1601 b (TIC 139375960)—based on data from NASA’s Transiting Exoplanet Survey Satellite (TESS). The five planets were identified from the full-frame images and were confirmed through a series of photometric and spectroscopic follow-up observations by the TESS Follow-up Observing Program Working Group. The planets are all Jovian size (R {sub P} = 1.01–1.77 R {sub J}) and have masses that range from 0.85 to 6.33 M {sub J}. The host stars of these systems have F and G spectral types (5595 ≤ T {sub eff} ≤ 6460 K) and are all relatively bright (9.5 < V < 10.8, 8.2 < K < 9.3), making them well suited for future detailed characterization efforts. Three of the systems in our sample (TOI-640 b, TOI-1333 b, and TOI-1601 b) orbit subgiant host stars (log g < 4.1). TOI-640 b is one of only three known hot Jupiters to have a highly inflated radius (R {sub P} > 1.7 R {sub J}, possibly a result of its host star’s evolution) and resides on an orbit with a period longer than 5 days. TOI-628 b is the most massive, hot Jupiter discovered to date by TESS with a measured mass of 6.31{sub −0.30}{sup +0.28} M {sub J} and a statistically significant, nonzero orbital eccentricity of e = 0.074{sub −0.022}{sup +0.021}. This planet would not have had enough time to circularize through tidal forces from our analysis, suggesting that it might be remnant eccentricity from its migration. The longest-period planet in this sample, TOI-1478 b (P = 10.18 days), is a warm Jupiter in a circular orbit around a near-solar analog. NASA’s TESS mission is continuing to increase the sample of well-characterized hot and warm Jupiters, complementing its primary mission goals.

79 ASTRONOMY AND ASTROPHYSICS↗

Groundwater Variability Across Temporal and Spatial Scales in the Central and Northeastern U.S.

Depth-to-water measurements from 181 monitoring wells in unconfined or semi-confined aquifers in nine regions of the central and northeastern U.S. were analyzed. Groundwater storage exhibited strong seasonal variations in all regions, with peaks in spring and lows in autumn, and its interannual variability was nearly unbounded, such that the impacts of droughts, floods, and excessive pumping could persist for many years. We found that the spatial variability of groundwater storage anomalies (deviations from the long term mean) increases as a power function of extent scale (square root of area). That relationship, which is linear on a log-log graph, is common to other hydrological variables but had never before been shown with groundwater data. We describe how the derived power function can be used to determine the number of wells needed to estimate regional mean groundwater storage anomalies with a desired level of accuracy, or to assess uncertainty in regional mean estimates from a set number of observations. We found that the spatial variability of groundwater storage anomalies within a region often increases with the absolute value of the regional mean anomaly, the opposite of the relationship between soil moisture spatial variability and mean. Recharge (drainage from the lowest model soil layer) simulated by the Variable Infiltration Capacity (VIC) model was compatible with observed monthly groundwater storage anomalies and month-to-month changes in groundwater storage.

Groundwater↗

Groundwater Variability Across Temporal and Spatial Scales in the Central and Northeastern U.S.

Depth-to-water measurements from 181 monitoring wells in unconfined or semi-confined aquifers in nine regions of the central and northeastern U.S. were analyzed. Groundwater storage exhibited strong seasonal variations in all regions, with peaks in spring and lows in autumn, and its interannual variability was nearly unbounded, such that the impacts of droughts, floods, and excessive pumping could persist for many years. We found that the spatial variability of groundwater storage anomalies (deviations from the long term mean) increases as a power function of extent scale (square root of area). That relationship, which is linear on a log-log graph, is common to other hydrological variables but had never before been shown with groundwater data. We describe how the derived power function can be used to determine the number of wells needed to estimate regional mean groundwater storage anomalies with a desired level of accuracy, or to assess uncertainty in regional mean estimates from a set number of observations. We found that the spatial variability of groundwater storage anomalies within a region often increases with the absolute value of the regional mean anomaly, the opposite of the relationship between soil moisture spatial variability and mean. Recharge (drainage from the lowest model soil layer) simulated by the Variable Infiltration Capacity (VIC) model was compatible with observed monthly groundwater storage anomalies and month-to-month changes in groundwater storage.

water cycle↗

OLCF Summit Supercomputer GPU Snapshots During Double-Bit Errors and Normal Operations

As we move into the exascale era, the power and energy footprints of high-performance computing (HPC) systems have grown significantly larger. Due to the harsh power and thermal conditions the system, components are exposed to extreme operating conditions. Operation of such modern HPC systems requires deep insights into long term system behavior to maintain its efficiency as well as its longevity. To help the HPC community to gain such insights, we provide double-bit errors using system telemetry data and logs collected from the Summit supercomputer, equipped with 27,648 Tesla V100 GPUs with 2nd-generation high-bandwidth memory (HBM2). The dataset relies on Nvidia XID records internally collected by GPU firmware at the time of failure occurrence, on the reboot-time logs of each Summit node, on node-level job scheduler records collected after each job termination, and on a 1Hz data rate from the baseboard management controllers (BMCs) of each Summit compute node using the OpenBMC event subscription protocol. Technical details can be found in the paper Oles et. al “Understanding GPU Memory Corruption at Extreme Scale: The Summit Case Study” ICS’24 (https://doi.org/10.1145/3650200.3656615).

97 MATHEMATICS AND COMPUTING↗

The effect of Cr, Co, Al, Mo and Ta on a series of cast Ni-base superalloys on the stability of an aluminide coating during cyclic oxidation in Mach 0.3 burner rig

The influence of varying the content of Co, Cr, Mo, Ta, and Al in a series of cast Ni-based gamma/gamma'superalloys on the behavior of aluminide coatings was studied in burner rig cyclic oxidation tests at 1100 C. The alloys had nominally fixed levels of Ti, W, Cb, Zr, C, and B. The alloy compositions were based on a full 2(sup 5)-fractional statistical design supplemented by 10 star point alloys and a center point alloy. This full central composite design of 43 alloys plus two additional alloys with extreme Al levels allowed a complete second degree estimating equation to be derived from the 5-compositional variables. The weight change/time data for the coated samples fitted well to the paralinear oxidation model and enabled a modified oxidation attack parameter, K'(sub a) to be derived to rank the alloys and log K' (sub a ) to be used as the dependent variable in the estimating equation to determine the oxidation resistance of the coating as a function of the underlying alloy content. The most protective aluminide coatings are associated with the highest possible base ally contents of CR and Al and at a 4 percent Ta level. The Mo and Co effects interact but at fixed levels of 0, 5, or 10% Co. A 4% Mo level is optimum.

Zaplatynsky, I.↗

Electron-H2 Collisions Studied Using the Finite Element Z-Matrix Method

We have applied the Z-matrix method, using a mixed basis of finite elements and Gaussians, to study e-H2 elastic and inelastic collisions. Special attention is paid to the quality of the basis set and the treatment of electron correlation. The calculated cross sections are invariant, to machine accuracy, with respect to the choice of parameters a, b, d, e as long as they satisfy Equation (3). However, the log derivative approach, i.e., the choice a = -e = 1, b = d = 0 appears to converge slightly faster than other choices. The cross sections agree well with previous theoretical results. Comparison will be made with available experimental data.

Huo, Winifred M.↗

Collecting and Processing Earth Science Data Metrics at NASA ESDIS

Since the launch of Terra satellite in 1999, the number of Earth Science remote sensing data products created and distributed by NASA's Earth Observing System (EOS) Data and Information System (EOSDIS) has increased from a few hundred to nearly ten thousand. NASA's Earth Science Data and Information System (ESDIS) Metrics System (EMS) collects metrics on data ingest, archive, and distribution by its Distributed Active Archive Centers (DAACs) and the Science Investigator-led Systems (SIPS), known as Data Providers. These metrics are critical in helping NASA management as well as data producers in resource planning and gaining a wide range of knowledge of data users and data usage.EMS receives flat files, or log files of data archive, ingest, and distribution either in their raw format, such as Apache web logs, or text files of log records formatted by the Data Providers. Tens of millions of records are processed each day to extract metrics on data products, user information, distribution protocols and services, and so on. The metrics are then made available to designated parties.This presentation provides an overview of the EMS processing workflow and improvement efforts made in recent years to handle ever-increasing number of data records and new metrics requirements, discusses several key steps including mapping log records to data products and identifying user communities along with geo-distribution, and demonstrates typical metrics capabilities produced by the EMS system. Challenges and potential approaches to improve the system are also discussed.

Pan, Jianfu↗

Oilwell Conversion (Well API 121913310501) to Geothermal Heat Storage Well for Flexible Electricity Storage

Geothermal growth is limited by a lack of geographically dispersed high-temperature thermal resources and high initial upfront investment in characterization and well construction. This project intended to address the challenges of energy supply intermittency and enhance grid resilience, reliability, and energy security by storing energy provided from currently available renewable resources in the subsurface to harvest it a later time during at-peak energy demand. This project intended to improve geothermal adoption, reduce initial project risk, and improve price competitiveness through utilizing existing oil and gas infrastructure such as non-productive wells, non-economic fields, dry holes, and orphaned wells. The project also intended to address the lack of geographically dispersed thermal resources and enhance grid resilience, reliability, and energy security by introducing an economical method for storing energy from currently available renewable resources in the subsurface for usage during at-peak energy demand. During this research, the project furthered the understanding of the feasibility of utilizing abandoned oil and gas wells as geothermal heat storage wells. The project team investigated the heat storage and hydrogeological characteristics of subsurface reservoirs in the Illinois Basin to evaluate their response to heat injection for determining the evolution of temperature profiles and heat losses over time using existing and available data sets. The project team then performed modeling and simulation to evaluate the heat losses of returning fluids during heat extraction. The outputs were used to select an optimal candidate reservoir and location in Southern Illinois. The team designed and performed a small-scale field test in an existing oil well to refine the model and to demonstrate the permitting and regulatory pathways necessary for the conversion of oil and gas assets to geothermal use. The field test also serves as a proof of concept and can guide the procedures for future research and implementation. Additionally, the project team, conducted initial market research and customer discovery to develop a go to market strategy for an Advanced Geothermal Energy Storage (AGES) system. The project team in this research also identified the parameters to be refined in future research, to improve the current go to market strategy economic model. To this end several subject matter experts were also identified to assist in future research with geothermal infrastructure setup, energy storage policy and law, energy storage market demand, potential siting based on demand etc. Future research will involve further sophistication of the site commercial modeling, implementing a larger-scale test, and further refinement of the thermodynamic modeling/simulation process. The output will be lifecycle costs and economics suitable for comparison to alternative approaches from a validated full-scale demonstration for venture capital investment into this technology. The project successfully demonstrated the ability to leverage existing oilfield infrastructure, permits, and land access and leasing agreements, to enable geothermal storage projects to come online faster and cheaper than a greenfield development could. This technology could allow for greater energy independence and security through long-term energy storage solutions. The longer duration allows for greater storage for renewables currently limited by hours-long storage durations of lithium-ion. The AGES system would support the growth of renewable energy farms, and provide greater opportunities for a cleaner energy infrastructure.

15 GEOTHERMAL ENERGY↗

Clutter Distributions for Tomographic Image Standardization in Ground-Penetrating Radar

Multistatic ground-penetrating radar (GPR) signals can be imaged tomographically to produce 3-D distributions of image intensities. In the absence of objects of interest, these intensities can be considered to be estimates of clutter. These clutter intensities spatially vary over several orders of magnitude and vary across different arrays, which makes a direct comparison of these raw intensities difficult. However, by gathering statistics on these intensities and their spatial variation, a variety of metrics can be determined. In this study, the clutter distribution is found to fit better to a two-parameter Weibull distribution than Gaussian or log-normal distributions. Based on the spatial variation of the two Weibull parameters, scale and shape, more information may be gleaned from these data. How well the GPR array is illuminating various parts of the ground, in depth and cross track, may be determined from the spatial variation of the Weibull scale parameter, which may in turn be used to estimate an effective attenuation coefficient in the soil. The transition in depth from clutter- to noise-limited conditions (which is one possible definition of GPR penetration depth) can be estimated from the spatial variation of the Weibull shape parameter. Lastly, the underlying clutter distributions also provide an opportunity to standardize image intensities to determine when a statistically significant deviation from background (clutter) has occurred, which is convenient for buried threat detection algorithm development that needs to be robust across multiple different arrays.

42 ENGINEERING↗

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects by Incorporating Machine Learning Techniques

The overall objective of this project is to develop a framework for reliable characterization and prediction of the state of stress in the overburden and underburden (including the basement) in CO 2 storage reservoirs using machine learning and integrated geomechanics and geophysical methods. Specifically, we propose to develop workflow encompassing of technologies and/or methods to predict stress and pressure changes due to CO 2 injection in an active tertiary recovery site and their impacts on subtle fault activation, fractures and occurrence of microseismic events and compare responses to field observations. In this project, we anticipate using dataset from the Farnsworth field Unit (FWU) which is operated by Purdure Petroleum. A novel elastic-waveform VSP inversion technique will be used to estimate high-resolution spatial and temporal changes of elastic moduli in CO 2 storage reservoirs, which will be combined with velocity-stress relationship derived from laboratory tests to obtain subsurface pressure and stress. Clustered microseismic data will be jointly inverted for improved focal mechanisms. Least-squares reverse-time migration of microseismic waveform data will be performed to directly image fracture/fault zones. Additionally, a deep neural network machine learning technique with convolutional and recurrent layers will be used for learning the spectro-temporal structures in microseismic waveforms. The results of this geotechnical data analysis will be integrated to develop a high-resolution 3D mechanical earth model extending from the overburden sealing formations to the underburden including the basement. Mechanical properties will be derived through integration of mechanical logs, tests, available results from chemo-mechanical laboratory tests, and elastic inversion of seismic data using a combination of Bayesian and stochastic methods as well as machine learning technique. Failure features (faults/fractures) will be represented and/or modeled based on seismic and core data analysis. A transient hydrodynamic-geomechanical model will be developed through coupling with the calibrated FWU reservoir simulation model. The full physics coupled model will be used to train a reduced order proxy model using machine learning algorithm for estimating stress which will then be used with appropriate constitutive relationships and forward seismological models to simulate pressure changes and induced microseismicity. An advanced optimization framework will be developed to perform a history match to minimize error between field observations and simulated. The history matched proxy model will be verified against the full-physics equivalent. The field observations that will be used in the coupled model calibration process include pressure/stress inverted from VSP, moment magnitude from microseismic analysis, real time downhole pressure measurements, production and injection data. Parameter sensitivity and uncertainty analysis will be performed to characterize the impact of model parameter uncertainty on stress estimates. The proposed project will have significant impact on future field implementation of the proposed technology. Because the project field site is an ongoing CO 2 EOR development, the value of the new technology will be demonstrated in an operational context and evaluated as a viable risk mitigation strategy. Cost/benefit will be evaluated together with the various commercial incentives for CO 2 sequestration available to oil and gas operators. The extensive available dataset and ongoing data acquisition under the SWP Phase III work plan provides flexibility for investigation of multiple approaches and reduces technical risk.

58 GEOSCIENCES↗