Search NASA⌕ Search

SEARCH · Search NASA

Results for “well fractures”

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 307 records · Page 17

Coupling Thermal-Hydraulics with Reactor Pressure Vessel Fracture Models

Because of the importance of maintaining the integrity of the reactor pressure vessel (RPV) in light-water reactors, significant effort has been devoted over the last several decades to understand the potential failure mechanisms, characterize how material properties evolve due to environmental conditions, and develop tools to assess the probability of fracture in aged RPVs during potential transient events. One important aspect of this problem that has not yet been fully addressed is the effects of nonuniform cooling on the RPV inner wall due to the way low-temperature water injected into the primary coolant system is distributed within the RPV during potential loss-of-coolant accident transients. Colder temperatures are expected in the regions near the inlets, causing a “cold-plume” effect. Rapid cooling of the RPV wall increases the likelihood of fracture initiation, and localized cooling due to cold-plume effects could potentially exacerbate this problem. The Grizzly code can perform a three-dimensional simulation of the RPV thermomechanical response as well as the probabilistic assessment of fracture for a population of pre-existing flaws, which is necessary to account for spatial variations in coolant temperature. This report documents a first effort to couple Grizzly RPV fracture models with computational fluid dynamics models performed using the Cardinal code, which allows for a realistic representation of the temperature distribution. This report demonstrates this coupling approach on two test cases: a simplified proof-of-concept scenario and a realistic small-break loss-of-coolant accident scenario. In both cases, the effects of the spatially varying temperature are evident; although, because the small-break loss-of-coolant accident scenario was far from challenging the RPV’s margins, it is still inconclusive how important cold-plume effects are in realistic accident scenarios. However, this effort represents an important step toward performing realistic analyses to determine to what extent spatially nonuniform cooling might affect RPV integrity under more aggressive accident scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The Influence of an Interlayer on Dual Hydraulic Fractures Propagation

Multi-cluster hydraulic fracturing of long-range horizontal wells is an approach for enhancing the productivity of low-permeability shale reservoirs. In this study, RFPA-Petrol (rock failure process analysis on petroleum problems) is applied for modeling hydraulic fracture propagation in multilayered formations. RFPA-Petrol based on coupled hydraulic-mechanical-damage (HMD) modeling was first tested by modeling a laboratory scale experiment on a physical (cement) model with a single completion. The modeling demonstrated the capability of RFPA-Petrol for simulating hydraulic fracture propagation. Then, we used RFPA-Petrol to investigate how the difference in material properties between oil-bearing layers and interlayers and the fracturing fluid properties influence the propagation of dual fractures in multilayered laboratory-scale models. In this case, the models with geological discontinuities in the vertical direction are strongly heterogeneous and RFPA-Petrol simulations successfully modeled the fracture configurations.

58 GEOSCIENCES↗

NiAl alloys for structural uses

Alloys based on the intermetallic compound NiAl are of technological interest as high temperature structural alloys. These alloys possess a relatively low density, high melting temperature, good thermal conductivity, and (usually) good oxidation resistance. However, NiAl and NiAl-base alloys suffer from poor fracture resistance at low temperatures as well as inadequate creep strength at elevated temperatures. This research program explored macroalloying additions to NiAl-base alloys in order to identify possible alloying and processing routes which promote both low temperature fracture toughness and high temperature strength. Initial results from the study examined the additions of Fe, Co, and Hf on the microstructure, deformation, and fracture resistance of NiAl-based alloys. Of significance were the observations that the presence of the gamma-prime phase, based on Ni3Al, could enhance the fracture resistance if the gamma-prime were present as a continuous grain boundary film or 'necklace'; and the Ni-35Al-20Fe alloy was ductile in ribbon form despite a microstructure consisting solely of the B2 beta phase based on NiAl. The ductility inherent in the Ni-35Al-20Fe alloy was explored further in subsequent studies. Those results confirm the presence of ductility in the Ni-35Al-20Fe alloy after rapid cooling from 750 - 1000 C. However exposure at 550 C caused embrittlement; this was associated with an age-hardening reaction caused by the formation of Fe-rich precipitates. In contrast, to the Ni-35Al-20Fe alloy, exploratory research indicated that compositions in the range of Ni-35Al-12Fe retain the ordered B2 structure of NiAl, are ductile, and do not age-harden or embrittle after thermal exposure. Thus, our recent efforts have focused on the behavior of the Ni-35Al-12Fe alloy. A second parallel effort initiated in this program was to use an alternate processing technique, mechanical alloying, to improve the properties of NiAl-alloys. Mechanical alloying in the conventional sense requires ductile powder particles which, through a cold welding and fracture process, can be dispersion strengthened by submicron-sized oxide particles. Using both the Ni-35Al-Fe alloys to contain approx. 1 v/o Y2O3. Preliminary results indicate that mechanically alloyed and extruded NiAl-Fe + Y2O3 alloys when heat treated to a grain-coarsened condition, exhibit improved creep resistance at 1000 C when compared to NiAl; oxidation resistance comparable to NiAl; and fracture toughness values a factor of three better than NiAl. As a result of the research initiated on this NASA program, a subsequent project with support from Inco Alloys International is underway.

Koss, D. A.↗

Developing Upscaling Approach for Swarming Hydraulic Fractures Observed at Hydraulic Fracturing Test Site through Multiscale Simulations

We show that this work aims to address a challenge posed by recent observations of tightly spaced hydraulic fractures in core samples from the hydraulic fracturing test site (HFTS) in the Middle Wolfcamp Formation. Many fractures in retrieved cores have subfoot spacing, which is at odds with conventional models in which usually one hydraulic fracture is initiated per cluster. Models assuming a single fracture at each cluster, although a common practice, often predict excessive fracture propagation that is inconsistent with microseismic observation. Here, we aim to develop a numerical approach to effectively account for densely spaced hydraulic fractures in field-scale simulations. Because it is impractical to explicitly model all aforementioned fractures, we develop a new upscaling law that enables existing simulation tools to predict reservoir response to fracture swarms. The upscaling law is derived based on an energy equivalence argument and validated through multiscale simulations using a high-fidelity code, GEOS. The swarming fractures are first modeled with a spacing that is much smaller than the cluster spacing; these fractures are then approximated by an upscaled, single fracture based on the proposed upscaling law. The upscaled fracture is shown to successfully match the energy input rate and produce the total fracture aperture and average propagation length of the explicitly simulated swarm. Afterward, the upscaling approach is further implemented in 3D field-scale simulations and validated against the HFTS microseismic data of a horizontal well. Our results show that hydraulic fracture swarming can significantly affect fracture propagation behaviors compared with the propagation of single fractures as assumed by conventional modeling approaches. Under the considered situations, the conventional treatment yields fast propagation speed that far exceeds that indicated by the microseismic data. We also illustrate that this discrepancy can be reduced readily through the implementation of the upscaling law. Our results demonstrate the importance of accounting for the fracture swarming effect in field-scale simulations and the efficacy of this approach to enable realistic predictions of reservoir responses to fracture swarms, without the need to model tightly spaced fractures individually.

02 PETROLEUM↗

Investigation of Post-Flight Solid Rocket Booster Thermal Protection System

After every Shuttle mission, the Solid Rocket Boosters (SRBs) are recovered and observed for missing material. Most of the SRB is covered with a cork-based thermal protection material (MCC-l). After the most recent shuttle mission, STS-114, the forward section of the booster appeared to have been impacted during flight. The darkened fracture surfaces indicated that this might have occurred early in flight. The scope of the analysis included microscopic observations to assess the degree of heat effects and locate evidence of the impact source as well as chemical analysis of the fracture surfaces and recovered foreign material using Fourier Transform Infrared Spectroscopy and Scanning Electron Microscopy/Energy Dispersive Spectroscopy. The amount of heat effects and presence of soot products on the fracture surface indicated that the material was impacted prior to SRB re-entry into the atmosphere. Fragments of graphite fibers found on these fracture surfaces were traced to slag inside the Solid Rocket Motor (SRM) that forms during flight as the propellant is spent and is ejected throughout the descent of the SRB after separation. The direction of the impact mark matches with the likely trajectory of SRBs tumbling prior to re-entry.

Nelson, Linda A.↗

Key Learnings from Hydraulic Fracturing Test Site - 2 (HFTS-2), Delaware Basin

Hydraulic fractures (HF) play a dominant role in oil and gas production from unconventional reservoirs, fracture characteristics, geometry, and extent are generally not well known and are typically inferred from various monitoring methods or subsurface models. The Hydraulic Fracturing Test Site-2 (HFTS-2) is a field-based research experiment performed in the Wolfcamp formation of the Permian (Delaware) Basin. HFTS-2, a unique joint industry program (JIP) that brings together government, academia, and industry, focused on field data acquisition and interpretation. The project has built a comprehensive HF research dataset that includes microseismic (MS), 1,500 ft of core, diagnostic formation injection tests (DFIT), fluid characterization, advanced logs, geochemistry, downhole gauges, varying completions, and permanent fiber optic (FO) cable in vertical and horizontal wells. These data were integrated and used to characterize HF geometry/dimensions and resulting depletion profiles, gain insights for completion design, and achieve a better understanding of parent well depletion effects on child well performance. The analysis and integration of the comprehensive HFTS-2 diagnostics dataset provide insights for optimal well placement and completions and improves understanding about the effects of parent well depletion. Furthermore, this dataset can be used to blind test and develop subsurface workflows for unconventional resources. Key learnings from HFTS-2 and the new technologies tested can easily be transferred to other unconventional plays.

58 GEOSCIENCES↗

Oil Recovery Prediction for Polymer Flood Field Test of Heavy Oil on Alaska North Slope Via Machine Assisted Reservoir Simulation

The first ever polymer flood field pilot to enhance the recovery of heavy oils on the Alaska North Slope is ongoing. This study constructs and calibrates a reservoir simulation model to predict the oil recovery performance of the pilot through machine-assisted reservoir simulation techniques. To replicate the early water breakthrough observed during waterflooding, transmissibility contrasts are introduced into the simulation model, forcing viscous fingering effects. In the ensuing polymer flood, these transmissibility contrasts are reduced to replicate the restoration of injection conformance during polymer flooding, as indicated by a significant decrease in water cut. Later, transmissibility contrasts are reinstated to replicate a water surge event observed in one of the producing wells during polymer flooding. This event may represent decreased injection conformance from fracture overextension; its anticipated occurrence in the other production well is included in the final forecast. The definition of polymer retention in the simulator incorporates the tailing effect reported in laboratory studies; this tailing effect is useful to the simultaneous history match of producing water cut and produced polymer concentration. The top 24 best-matched simulation models produced at each stage of the history matching process are used to forecast oil recovery. The final forecast clearly demonstrates that polymer flooding significantly increases the heavy oil production for this field pilot compared to waterflooding alone. This exercise displays that a simulation model is only valid for prediction if flow behavior in the reservoir remains consistent with that observed during the history matched period. Critically, this means that a simulation model calibrated for waterflooding may not fully capture the benefits of an enhanced oil recovery process such as polymer flooding. Therefore, caution is recommended in using basic waterflood simulation models to scope potential enhanced oil recovery projects.

Keith, Cody Douglas↗

Coupled Thermo–Hydro–Mechanical Modeling of Hydro-Shearing Stimulation in an Enhanced Geothermal System in the Raft River Geothermal Field, USA

The Raft River Enhanced Geothermal System (EGS) demonstration project aims to improve the geothermal production by enhancing the reservoir permeability via shear stimulation in well RRG-9. In this paper, we performed a series of 3D thermo– hydro–mechanical (THM) simulations to investigate the influence of hydraulic and thermal effects on the development of this EGS. The model includes synthetic fracture populations based on borehole televiewer images and in-situ stress measurements from well RRG-9. Fracture permeability evolution is determined using an empirical permeability law developed from laboratory experiments. The model was calibrated by comparing the hydraulic response to field observations, including wellhead pressure, injection rate, and well injectivity. Particularly, we analyzed the enhancement of reservoir permeability and the spatial extent of the stimulation zone for the given injection schedule. Our results indicate that the permeability enhancement of fractured geothermal reservoir is caused by the combined effects of injection-induced cooling and fluid pressure increase. The decrease of temperature plays a dominant role in reactivation of natural fractures under the hydro-shearing mechanism, while the higher injection pressure promotes shear failure and enlarges the stimulation zone. Finally, for the specified extensional stress state, the model favors greatest permeability enhancement along the maximum horizontal principal stress, a moderate vertical enhancement, and a smaller gain along the minimum horizontal principal stress.

58 GEOSCIENCES↗

Chemical and isotopic evolution of flowback fluids from the Utica Gas Shale Play, Eastern Ohio USA

For this work, hydraulic fracturing flowback fluids were collected from two Utica/Point Pleasant well pads in eastern Ohio. One site was in the wet gas zone (UPPW), while the site ~50 km south consisted of four wells on the same pad in the dry gas zone (UPPS). Samples of input fluids also were collected before and during hydraulic fracturing. Flowback fluids are Na-Ca-Cl brines with total dissolved salt (TDS) concentrations that increase over several months from ~100 to 200 g/L. The slightly higher TDS of the dry gas fluids are in part due to recycled flowback used as input fluids, and in part due to lower volume of water used in hydraulic fracturing. Concentrations of most major ions (Ca 2+ , Mg 2+ , Na + , Sr2 + , Fe, Mn, Cl – , Br – ) are similar for the five wells sampled, although small but systematic changes occur in the major element ratios over time. Most notably, an increase in the Sr/Cl ratio corresponds to an increase in the 87 Sr/ 86 Sr ratio in the fluids, suggesting interactions with a more radiogenic Sr source in the subsurface. Dissolved Ba concentrations and Ra activities were different between the two sites, reflecting a high SO 4 2– fluid used at the UPPW site, and water-rock reactions occurring during hydraulic fracturing at the UPPS site. Water oxygen (δ 18 O) and hydrogen (δD) isotopes for input fluids used in the UPPW4 well and fresh water used for the UPPS wells fall on the Global Meteoric waterline (GMWL). Flowback fluids from both sites are relatively enriched in δ 18 O and δD compared to the input, but do not appear to follow a simple mixing trend, suggesting reaction and isotopic exchange with carbonates and fractionation due to imbibition in the rock. Cl isotopes, δ 37 Cl in the FP fluids varied from ~ –0.43 to +0.13‰, the largest variation was observed in the earlier stages of flowback, while in the later stages δ 37 Cl exhibited a small but systematic increase over time, suggesting diffusion control of isotopic composition. Some of the trace species measured (dissolved Fe, Mn, thiosulfate and organic acids) do not follow the same trends as the major ions, suggesting contributions from input fluids, or microbially mediated reactions are exerting control on their concentrations. Input water chemistry exerts an important control on the Sr and Ba concentrations in flowback water. High SO 4 2– in the input fluids used for hydraulically fracturing the wet gas well leads to precipitation of barite-celestite in flowback fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Inference of Induced Fracture Geometries Using Fiber-Optic Distributed Strain Sensing in Hydraulic Fracture Test Site 2

Fiber-optic Distributed Strain Sensor observations in a vertical well allow rapid inference of the time evolution of nearby hydraulic fracture heights and their vertical aperture distributions during stimulation and depressurization. We describe a computationally efficient deconvolutional model based on an elastostatic kernel for fracture opening within an isotropic, homogeneous, linear elastic medium. Furthermore, an application of this model to observations from the Hydraulic Fracture Test Site 2 reveals that ultimate fracture heights are typically greater than those indicated by microseismic monitoring, and some hydraulic fractures continue their upward growth for several hours after high-pressure injections ceased.

58 GEOSCIENCES↗

Structural Design Parameters for Germanium

The fracture toughness and slow crack growth parameters of germanium supplied as single crystal beams and coarse grain disks were measured. Although germanium is anisotropic (A* 1.7), it is not as anisotropic as SiC, NiAl, or Cu. Thus the fracture toughness was similar on the 100, 110, and 111 planes, however, measurements associated with randomly oriented grinding cracks were 6 to 30 higher. Crack extension in ring loaded disks occurred on the 111 planes due to both the lower fracture energy and the higher stresses on stiff 111 planes. Germanium exhibits a Weibull scale effect, but does not exhibit significant slow crack growth in distilled water. (n 100), implying that design for quasi static loading can be performed with scaled strength statistics. Practical values for engineering design are a fracture toughness of 0.69 0.02 MPam (megapascals per square root meter) and a Weibull modulus of m 6 2. For well ground and reasonable handled coupons, average fracture strength should be greater than 40 megapascals. Aggregate, polycrystalline elastic constants are Epoly 131 gigapascals, vpoly 0.22.

fracture toughness↗

Utah FORGE Project 3-2417: Simulations for Distributed Acoustic Sensing Strain Signatures as an Indicator of Fracture Connectivity

This dataset encompasses simulations of strain signatures from both hydraulically connected and "near-miss" fractures in enhanced geothermal systems (EGS). The files and results are presented from the perspective of digital acoustic sensing's (DAS) potential to differentiate the two fracture types. This dataset was acquired by the FOGMORE R&D project (Fiber Optic Geophysical MOnitoring of Reservoir Evolution), under Utah FORGE R&D Project 3-2417. Included are simulation and results via MatLab and COMSOL files, as well as a thesis and paper summarizing the results. Some stimulated fractures may be incomplete, approaching but not intersecting the production well. These "near-miss" fractures can be addressed in future stimulation stages or re-stimulated to complete the connection. We propose the use of fiber optic distributed acoustic sensing (DAS) as a method by which near-miss stimulated fractures may be identified and distinguished from hydraulically connected fractures. The low-frequency sub-nanostrain signatures of both complete and near-miss fractures in DAS data are simulated in this study using a hydrogeomechanical discrete fracture network model. The spatial distribution of strain was found to be an accurate indicator. However, this indicator must be evaluated in the context of DAS gauge length and spatial sampling. These simulations are a precursor to tests conducted at FORGE in 2023.

15 GEOTHERMAL ENERGY↗

Source apportionment of airborne volatile organic compounds near unconventional oil and gas development

Oil and natural gas (ONG) extraction emits volatile organic compounds (VOCs). Certain VOCs are identified as hazardous air pollutants (HAPS) while others contribute to ozone formation. This study examines the impact of ONG operations on VOC levels during the development of multi-well ONG pads in suburban Broomfield, Colorado. From October 2018 to December 2020, weekly VOC measurements were taken at 18 sites across the area. These included spots near well pads, in adjacent neighborhoods, and at a background site, covering various stages of well pad development including drilling, hydraulic fracturing, flowback, and production. Analysis using Positive Matrix Factorization (PMF) identified six factors, including combustion, background/biogenic sources, light and complex alkanes, drilling activities, and ONG acetylene. Factors linked to local ONG activities exhibited clear temporal and spatial correlations with Broomfield well development. Benzene source analysis revealed distinct contribution gradients, with ONG-related sources notably influencing areas near the well pads, particularly in pre-production. ONG-related weekly benzene contributions varied from 9% to 63% at a community background site and 18% to 89% in a neighborhood close to a well pad.

54 ENVIRONMENTAL SCIENCES↗

Titanium-Oxygen Reactivity Study

A program has been conducted at Astronautics to investigate the likelihood of occurrence of the catastrophic oxidation of titanium alloy sheet under conditions which simulate certain cases of accidental failure of the metal while it is in contact with liquid or gaseous oxygen. Three methods of fracturing the metal were used; they consisted of mechanical puncture, tensile fracture of welded joints, and perforation by very high velocity particles. The results of the tests which have been conducted provide further evidence of the reactivity of titanium with liquid and gaseous oxygen. The evidence indicates that the rapid fracturing of titanium sheet while it is in contact with oxygen initiates the catastrophic oxidation reaction. Initiation occurred when the speed of the fracture was some few feet per second, as in both the drop-weight puncture tests and the static tensile fracture tests of welded joints, as well as when the speed was several thousand feet per second, as in the simulated micrometeoroid penetration tests. The slow propagation of a crack, however, did not initiate the reaction. It may logically be concluded that the localized frictional heat of rapid fracture and/or spontaneous oxidation (exothermic) of minute particles emanating from the fracture cause initiation of the reaction. Under conditions of slow fracture, however, the small heat generated may be adequately dissipated and the reaction is not initiated. A portion of the study conducted consisted of investigating various means by which the reaction might be retarded or prevented. Providing a "barrier" at the titanium-oxygen interface consisting of either aluminum metal or a coating of a petroleum base corrosion inhibitor appeared to be only partially effective in retarding the reaction. The accidental puncturing or similar rupturing of thin-walled pressurized oxygen tanks on missiles and space vehicle will usually constitute loss of function, and may sometimes cause their catastrophic destruction by explosive decompression regardless of the type of material used for their construction. In the case of tanks constructed of titanium alloys the added risk is incurred of catastrophic burning of the tanks. In view of this it is recommended that thin-walled tanks constructed of titanium alloys should not be used to contain liquid or gaseous oxygen.

Chafey, J. E.↗

Direct Injection of Biomineralizing Agents to Restore Injectivity and Wellbore Integrity

Summary In this manuscript, we describe the second of two field demonstrations of microbially induced calcium carbonate precipitation (MICP) performed in a failed waterflood injection well in Indiana. In 2012, fracture-related flow pathways developed in the wellbore cement, causing injection water to bypass the oil-bearing formation and enter a high-permeability sandstone thief zone, thereby substantially decreasing injection pressure. In the first field demonstration, our study team characterized the well's mode of failure and successfully applied MICP to decrease flow through the defective cement. However, because the MICP treatment was conducted using a bailer delivery system, the degree of permeability reduction achievable was not adequate to fully restore the historic injection pressure of 1,400 psi at 1 gal/min. For the second field demonstration (reported herein), a direct injection system was developed that substantially increased the injection volume of MICP-promoting fluids. Two strategies were implemented to produce more ureolytic microbes: resuspending concentrated frozen cells immediately before injection and scaling up the bioreactor growth capacity. Multiple pulses of microbes and urea-calcium media were pumped into a string of 1-in.-diameter tubing separated by brine spacers and injected continuously at a flow rate of 3.4 to 1.4 gal/min. During the third day of injection, an injection pressure of 1,384 psi at a flow rate of 1.4 gal/min was achieved, and the experiment was terminated. This study demonstrates that MICP can be successfully used in large-volume applications where the time frame for the delivery of reactants is limited. This finding has significant relevance for commercialization of the MICP biotechnology in the oil and gas industry.

Engineering↗

Combined Structural Analysis of Core and Image Log of TGH 76-31 South East of Mt Baker, Washington State

Despite active volcanism, few geothermal energy resources have been developed in the Cascades Range. Temperature Gradient Hole 76-31 was drilled to ~440 m measured depth to probe for zones where fractures provide fluid conduits that transport deep volcanic heat to shallow depths that could support baseload, carbon neutral electrical generation. These zones were predicted by a Play Fairway assessment (PFA) of resource potential along a zone 11 km west-southwest of the summit of Mount Baker Volcano. Rock core, temperature logs, and an acoustic image log were obtained. By comparison to outcrops, the core has been initially interpreted as the Chilliwack group comprised of partially metamorphosed basaltic to andesitic volcanics. Core mapping reveals complex, steeply dipping networks of fractures and brecciation along slickensided strike slip faults; clay alteration is common in many of these structures. The majority of fractures are thoroughly healed by layers of chlorite and calcite; preservation of pyrite indicates anoxic conditions. The majority of fracture porosity resides in very dense fractures a few centimeters or less in length. The image log provides good insights into attitude of fractures that fully transect the core, but generally underestimates fracture density. The combination of complex, non-planar fracture zones containing many short fractures and healing promote misinterpretation of natural fracture attitude and density in the image log. The healing and anoxic minerals are consistent with the conductive temperature gradient measured in the well below a shallow isothermal zone, although, several fractures are open or only partially healed and resulted in fluid entries into the well. Here, the equilibrated measured temperature gradient of 64ºC/km and calculated heat flow of 145 mW/m2 is more than twice the regional average, indicating local influence of the Mt Baker magmatic system at the Little Park Creek TGH site.

Cascades↗

A comprehensive review of proppant embedment in shale reservoirs: Experimentation, modeling and future prospects

This paper provides a comprehensive review on the application of proppants to maintain fracture permeability over the lifetime of a well based on published observations from experiments and modeling. The review identifies and describes important processes occurring during proppant embedment, during hydraulic fracturing, laboratory testing of fracture conductivity, proppant embedment and modeling of proppant embedment. Lastly, this paper identifies the challenges and knowledge gaps that also provide future avenues of research and opportunities for collaborative technological development which requires an interdisciplinary approach of science, engineering in academia, government, and private sector.

03 NATURAL GAS↗