Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “horizontal”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Characterization of elastic mechanical properties of Tuscaloosa Marine Shale from well logs using the vertical transversely isotropic model

To avoid steep declines in the Tuscaloosa Marine Shale (TMS) production, wells are fracture-stimulated to release the hydrocarbons trapped in the matrix of the formation. An accurate estimation of Young’s modulus and Poisson’s ratio is essential for hydraulic fracture propagation. In addition, ignoring the highly heterogeneous and anisotropic character of TMS can lead to erroneous stress values, which subsequently affect hydraulic fracture width estimates and the overall hydraulic fracturing process. We have developed an empirical 1D geomechanical model that takes into account VTI anisotropy, and it is used to characterize the elastic mechanical properties of TMS in two wells. In the analyzed formation, the vertical Poisson’s ratio is less than the horizontal Poisson’s ratio, which suggests the necessity of an alternative to the ANNIE equations. The stiffness coefficients [Formula: see text] and [Formula: see text] were estimated using the relationships developed from the ultrasonic core data available for the two TMS. Further, correlations between the static and dynamic properties from laboratory tests were used to improve the minimum horizontal stress calculation. We compare VTI Young’s moduli, Poisson’s ratios, and minimum horizontal stress with the isotropic solution. VTI modeling improves the estimation of the elastic mechanical properties. The isotropic solution underestimates the minimum horizontal stress in the formation. Moreover, it was shown that the 20 ft shale interval below the TMS base is characterized by a low Young’s modulus (the vertical Young’s modulus is equal to 20 GPa, whereas the horizontal Young’s modulus is equal to 40 GPa) and may be a frac barrier.

Geochemistry & Geophysics↗
DOE OSTI2013

Slip and Dilation Tendency Analysis of the Tuscarora Geothermal Area

Critically stressed fault segments have a relatively high likelihood of acting as fluid flow conduits (Sibson, 1994). As such, the tendency of a fault segment to slip (slip tendency; Ts; Morris et al., 1996) or to dilate (dilation tendency; Td; Ferrill et al., 1999) provides an indication of which faults or fault segments within a geothermal system are critically stressed and therefore likely to transmit geothermal fluids. The slip tendency of a surface is defined by the ratio of shear stress to normal stress on that surface: Ts = T / on (Morris et al., 1996). Dilation tendency is defined by the stress acting normal to a given surface: Td = (o1-on) / (o1-o3) (Ferrill et al., 1999). Slip and dilation were calculated using 3DStress (Southwest Research Institute). Slip and dilation tendency are both unitless ratios of the resolved stresses applied to the fault plane by ambient stress conditions. Values range from a maximum of 1, a fault plane ideally oriented to slip or dilate under ambient stress conditions to zero, a fault plane with no potential to slip or dilate. Slip and dilation tendency values were calculated for each fault in the focus study areas at, McGinness Hills, Neal Hot Springs, Patua, Salt Wells, San Emidio, and Tuscarora on fault traces. As dip is not well constrained or unknown for many faults mapped in within these we made these calculations using the dip for each fault that would yield the maximum slip tendency or dilation tendency. As such, these results should be viewed as maximum tendency of each fault to slip or dilate. The resulting along-fault and fault-to-fault variation in slip or dilation potential is a proxy for along fault and fault-to-fault variation in fluid flow conduit potential. Stress Magnitudes and directions Stress field variation within each focus area was approximated based on regional published data and the world stress database (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2010; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012; Moeck et al., 2010; Moos and Ronne, 2010 and Reinecker et al., 2005) as well as local stress information if applicable. For faults within these focus systems we applied either a normal faulting stress regime where the vertical stress (sv) is larger than the maximum horizontal stress (shmax) which is larger than the minimum horizontal stress (sv>shmax>shmin) or strike-slip faulting stress regime where the maximum horizontal stress (shmax) is larger than the vertical stress (sv) which is larger than the minimum horizontal stress (shmax >sv>shmin) depending on the general tectonic province of the system. Based on visual inspection of the limited stress magnitude data in the Great Basin we used magnitudes such that shmin/shmax = .527 and shmin/sv= .46, which are consistent with complete and partial stress field determinations from Desert Peak, Coso, the Fallon area and Dixie valley (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2011; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012). Slip and dilation tendency for the Tuscarora geothermal field was calculated based on the faults mapped Tuscarora area (Dering, 2013). The Tuscarora area lies in the Basin and Range Province, as such we applied a normal faulting stress regime to the Tuscarora area faults, with a minimum horizontal stress direction oriented 115, based on inspection of local and regional stress determinations, as explained above. Under these stress conditions north-northeast striking, steeply dipping fault segments have the highest dilation tendency, while north-northeast striking 60 degrees dipping fault segments have the highest tendency to slip. Tuscarora is defined by a left-step in a major north- to-north northeast striking, west-dipping range-bounding normal fault system. Faults within the broad step define an anticlinal accommodation ...

15 GEOTHERMAL ENERGY↗
DOE OSTI2013

Slip and Dilation Tendency Analysis of the San Emidio Geothermal Area

Critically stressed fault segments have a relatively high likelihood of acting as fluid flow conduits (Sibson, 1994). As such, the tendency of a fault segment to slip (slip tendency; Ts; Morris et al., 1996) or to dilate (dilation tendency; Td; Ferrill et al., 1999) provides an indication of which faults or fault segments within a geothermal system are critically stressed and therefore likely to transmit geothermal fluids. The slip tendency of a surface is defined by the ratio of shear stress to normal stress on that surface: Ts = T / on (Morris et al., 1996). Dilation tendency is defined by the stress acting normal to a given surface: Td = (o1-on) / (o1-o3) (Ferrill et al., 1999). Slip and dilation were calculated using 3DStress (Southwest Research Institute). Slip and dilation tendency are both unitless ratios of the resolved stresses applied to the fault plane by ambient stress conditions. Values range from a maximum of 1, a fault plane ideally oriented to slip or dilate under ambient stress conditions to zero, a fault plane with no potential to slip or dilate. Slip and dilation tendency values were calculated for each fault in the focus study areas at, McGinness Hills, Neal Hot Springs, Patua, Salt Wells, San Emidio, and Tuscarora on fault traces. As dip is not well constrained or unknown for many faults mapped in within these we made these calculations using the dip for each fault that would yield the maximum slip tendency or dilation tendency. As such, these results should be viewed as maximum tendency of each fault to slip or dilate. The resulting along-fault and fault-to-fault variation in slip or dilation potential is a proxy for along fault and fault-to-fault variation in fluid flow conduit potential. Stress Magnitudes and directions Stress field variation within each focus area was approximated based on regional published data and the world stress database (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2010; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012; Moeck et al., 2010; Moos and Ronne, 2010 and Reinecker et al., 2005) as well as local stress information if applicable. For faults within these focus systems we applied either a normal faulting stress regime where the vertical stress (sv) is larger than the maximum horizontal stress (shmax) which is larger than the minimum horizontal stress (sv>shmax>shmin) or strike-slip faulting stress regime where the maximum horizontal stress (shmax) is larger than the vertical stress (sv) which is larger than the minimum horizontal stress (shmax >sv>shmin) depending on the general tectonic province of the system. Based on visual inspection of the limited stress magnitude data in the Great Basin we used magnitudes such that shmin/shmax = .527 and shmin/sv= .46, which are consistent with complete and partial stress field determinations from Desert Peak, Coso, the Fallon area and Dixie valley (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2011; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012). Slip and dilation tendency for the San Emidio geothermal field was calculated based on the faults mapped Tuscarora area (Rhodes, 2011). The San Emidio area lies in the Basin and Range Province, as such we applied a normal faulting stress regime to the San Emidio area faults, with a minimum horizontal stress direction oriented 115, based on inspection of local and regional stress determinations, as explained above. This is consistent with the shmin determined through inversion of fault data by Rhodes (2011). Under these stress conditions north-northeast striking, steeply dipping fault segments have the highest dilation tendency, while north-northeast striking 60 degrees dipping fault segments have the highest tendency to slip. Interesting, the San Emidio geothermal field lies in an area of primarily north striking faults, w...

15 GEOTHERMAL ENERGY↗
DOE OSTI2013

Slip and Dilation Tendency Analysis of the Salt Wells Geothermal Area

Critically stressed fault segments have a relatively high likelihood of acting as fluid flow conduits (Sibson, 1994). As such, the tendency of a fault segment to slip (slip tendency; Ts; Morris et al., 1996) or to dilate (dilation tendency; Td; Ferrill et al., 1999) provides an indication of which faults or fault segments within a geothermal system are critically stressed and therefore likely to transmit geothermal fluids. The slip tendency of a surface is defined by the ratio of shear stress to normal stress on that surface: Ts = T / on (Morris et al., 1996). Dilation tendency is defined by the stress acting normal to a given surface: Td = (o1-on) / (o1-o3) (Ferrill et al., 1999). Slip and dilation were calculated using 3DStress (Southwest Research Institute). Slip and dilation tendency are both unitless ratios of the resolved stresses applied to the fault plane by ambient stress conditions. Values range from a maximum of 1, a fault plane ideally oriented to slip or dilate under ambient stress conditions to zero, a fault plane with no potential to slip or dilate. Slip and dilation tendency values were calculated for each fault in the focus study areas at, McGinness Hills, Neal Hot Springs, Patua, Salt Wells, San Emidio, and Tuscarora on fault traces. As dip is not well constrained or unknown for many faults mapped in within these we made these calculations using the dip for each fault that would yield the maximum slip tendency or dilation tendency. As such, these results should be viewed as maximum tendency of each fault to slip or dilate. The resulting along-fault and fault-to-fault variation in slip or dilation potential is a proxy for along fault and fault-to-fault variation in fluid flow conduit potential. Stress Magnitudes and directions Stress field variation within each focus area was approximated based on regional published data and the world stress database (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2010; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012; Moeck et al., 2010; Moos and Ronne, 2010 and Reinecker et al., 2005) as well as local stress information if applicable. For faults within these focus systems we applied either a normal faulting stress regime where the vertical stress (sv) is larger than the maximum horizontal stress (shmax) which is larger than the minimum horizontal stress (sv>shmax>shmin) or strike-slip faulting stress regime where the maximum horizontal stress (shmax) is larger than the vertical stress (sv) which is larger than the minimum horizontal stress (shmax >sv>shmin) depending on the general tectonic province of the system. Based on visual inspection of the limited stress magnitude data in the Great Basin we used magnitudes such that shmin/shmax = .527 and shmin/sv= .46, which are consistent with complete and partial stress field determinations from Desert Peak, Coso, the Fallon area and Dixie valley (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2011; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012). Slip and dilation tendency for the Salt Wells geothermal field was calculated based on the faults mapped in the Bunejug Mountains quadrangle (Hinz et al., 2011). The Salt Wells area lies in the Basin and Range Province (N. Hinz personal comm.) As such we applied a normal faulting stress regime to the Salt Wells area faults, with a minimum horizontal stress direction oriented 105, based on inspection of local and regional stress determinations. Under these stress conditions north-northeast striking, steeply dipping fault segments have the highest dilation tendency, while north-northeast striking 60 degrees dipping fault segments have the highest tendency to slip. Several such faults intersect in high density in the core of the accommodation zone in the Bunejug Mountains and local to the Salt Wells geothermal .

15 GEOTHERMAL ENERGY↗
DOE OSTI2013

Slip and Dilation Tendency Analysis of McGinness Hills Geothermal Area

Slip and Dilation Tendency in focus areas Critically stressed fault segments have a relatively high likelihood of acting as fluid flow conduits (Sibson, 1994). As such, the tendency of a fault segment to slip (slip tendency; Ts; Morris et al., 1996) or to dilate (dilation tendency; Td; Ferrill et al., 1999) provides an indication of which faults or fault segments within a geothermal system are critically stressed and therefore likely to transmit geothermal fluids. The slip tendency of a surface is defined by the ratio of shear stress to normal stress on that surface: Ts = T / on (Morris et al., 1996). Dilation tendency is defined by the stress acting normal to a given surface: Td = (o1-on) / (o1-o3) (Ferrill et al., 1999). Slip and dilation were calculated using 3DStress (Southwest Research Institute). Slip and dilation tendency are both unitless ratios of the resolved stresses applied to the fault plane by ambient stress conditions. Values range from a maximum of 1, a fault plane ideally oriented to slip or dilate under ambient stress conditions to zero, a fault plane with no potential to slip or dilate. Slip and dilation tendency values were calculated for each fault in the focus study areas at, McGinness Hills, Neal Hot Springs, Patua, Salt Wells, San Emidio, and Tuscarora on fault traces. As dip is not well constrained or unknown for many faults mapped in within these we made these calculations using the dip for each fault that would yield the maximum slip tendency or dilation tendency. As such, these results should be viewed as maximum tendency of each fault to slip or dilate. The resulting along-fault and fault-to-fault variation in slip or dilation potential is a proxy for along fault and fault-to-fault variation in fluid flow conduit potential. Stress Magnitudes and directions Stress field variation within each focus area was approximated based on regional published data and the world stress database (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2010; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012; Moeck et al., 2010; Moos and Ronne, 2010 and Reinecker et al., 2005) as well as local stress information if applicable. For faults within these focus systems we applied either a normal faulting stress regime where the vertical stress (sv) is larger than the maximum horizontal stress (shmax) which is larger than the minimum horizontal stress (sv>shmax>shmin) or strike-slip faulting stress regime where the maximum horizontal stress (shmax) is larger than the vertical stress (sv) which is larger than the minimum horizontal stress (shmax >sv>shmin) depending on the general tectonic province of the system. Based on visual inspection of the limited stress magnitude data in the Great Basin we used magnitudes such that shmin/shmax = .527 and shmin/sv= .46, which are consistent with complete and partial stress field determinations from Desert Peak, Coso, the Fallon area and Dixie valley (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2011; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012). Slip and dilation tendency for the McGinness Hills geothermal field was calculated based on the faults mapped McGinness Hills area (Siler 2012, unpublished). The McGinness Hills area lies in the Basin and Range Province, as such we applied a normal faulting stress regime to the McGinness area faults, with a minimum horizontal stress direction oriented 115, based on inspection of local and regional stress determinations, as explained above. Under these stress conditions north-northeast striking, steeply dipping fault segments have the highest dilation tendency, while north-northeast striking 60 degrees dipping fault segments have the highest tendency to slip. The McGinness Hills geothermal system is characterized by a left-step in a north-northeast striking west-dipping fault system...

15 GEOTHERMAL ENERGY↗
DOE OSTI2019

Microseismic Focal Mechanisms and Implications for Changes in Stress during the 2014 Newberry EGS Stimulation

We adapt the relative polarity method from Shelly et al. (2016) to compute focal mechanisms for microearthquakes associated with the 2014 hydroshearing stimulation at the Newberry volcano geothermal site. We focus the analysis on events relocated by Aguiar and Myers (2018), who report that six event clusters predominantly comprise the 2014 sequence. Data quality allows focal mechanism analysis for four of the six event clusters. We use Hardebeck and Shearer (2002, 2003; hereafter HASH) to compute focal mechanisms based on first–motion polarities and S/P amplitude ratios. We manually determine P– and S–wave polarities for a well–recorded reference event in each cluster, then use waveform cross correlation to determine whether recordings of other events in the cluster are the same or reversed polarity at each network station. Here, most waveform polarities are consistent with the affiliated reference event, indicating similar focal mechanisms within each cluster. The deeper clusters are east–west–striking normal faults, whereas the shallower clusters, close to the top of the open–hole section of the borehole, are strike slip with east–west motion. Regional studies and prestimulation borehole breakouts find the maximum stress direction is vertical and maximum horizontal stress is approximately north–south. Fault geometry and focal mechanisms of microseismicity during the stimulation suggest that increased pressure from fluid injection predominantly caused changes in horizontal stress, consistent with predictions from numerical studies of stress change caused by fluid injection. At shallow depths, where previous studies suggest the difference between vertical and horizontal stress is lowest, injection appears to have rotated the direction of maximum stress from vertical to horizontal, resulting in strike–slip motion. At greater depth, vertical stress continued to be the dominant direction during the stimulation, but fault orientation indicates either reactivation of pre–existing fractures or rotation of the direction of maximum horizontal stress from approximately north–south to east–west.

15 GEOTHERMAL ENERGY↗
DOE OSTI2017

Brine Extraction and Treatment Strategies to Enhance Pressure Management and Control of CO 2 Plumes in Deep Geologic Formations

The overall goal of the this project is to develop and validate pressure management and carbon dioxide (CO 2 ) plume control strategies that can address technical and economic barriers to commercial deployment of CO 2 storage technologies, based on computational and field demonstration work at the Archer Daniels Midland Company (ADM) facility where the Illinois Basin–Decatur Project (IBDP) and the Illinois-Industrial Carbon Capture and Storage (IL-ICCS) projects are located. To accomplish the overall goal, the ISGS designed a brine extraction storage test (BEST) that could be completed in two phases. The goal of BEST Phase I was to evaluate the feasibilities of extraction well(s) placement, the brine extraction to CO 2 injection rate ratio, extraction well completion, and brine treatment and handling. The goal of BEST Phase II would be to validate the brine extraction and treatment options deemed feasible in Phase I by (1) demonstrating the efficacy of brine extraction (BE) in managing pressure (i.e., formation) and the CO 2 plume, and (2) demonstrating treatment of extracted brine with high total dissolved solids (TDS; >200,000 mg/L) using multiple advanced treatment technologies. This report details work done in Phase I. Several brine extraction and treatment scenarios were tested, simulated, and analyzed for their effectiveness in extracting brine. Initially a vertical well was studied; however, geologic modeling, reservoir modeling, and the existing facility and wellbore infrastructure dictated that the location of a vertical brine extraction well was limited to an area with no existing monitoring wells and where the well would be in relative proximity to an existing CO 2 plume. Consequently, a vertical well was excluded, and a horizontal brine extraction well placed above the existing CO 2 plume near two existing wells was studied. The horizontal well option allows the project to leverage the availability of cased-hole logs and cross-well tomography to monitor CO 2 saturation and plume distribution, respectively. Because of the proximity of the horizontal well option to two existing wells, no additional monitoring well (or caprock penetration) is required. The recommended brine extraction pilot design options are (1) a horizontal extraction well at the base of the Middle Mt. Simon, which is 350–520 ft (107–158 m) above the CO 2 plume at CCS#1 and VW#1; or (2) a vertical extraction well 0.5 mi (0.8 km) from CCS#2 in a direction approximately southeast of CCS#2, perpendicular to the direction of high hydraulic connectivity. A horizontal extraction well has advantages over a vertical extraction well, including less risk of drilling into an existing CO 2 plume and it can be located between two other wells that can be used for monitoring. Thus, because the two existing wells can serve as monitoring wells, it eliminates the need for a third verification well and allows for a lower extraction rate to control the CO 2 plume and pressure. Managing pressure and the CO 2 plume distribution via brine extraction creates the obvious and important challenge of handling and treating the extracted brine. There were three options for brine disposal: (1) underground injection control (UIC) disposal well, (2) brine treatment and industrial use, and (3) brine pretreatment and discharge into municipal wastewater system. The primary design elements were budget and permitting requirements. The disposal well would be a vertical well drilled and completed into the Potosi Dolomite. For the range of extraction rates anticipated, the cost of this well is relatively constant. The cost of brine treatment is highly depends on the extraction rate, which depends on the well orientation. If relatively high rates are required, the vertical disposal well option is more favorable; for relatively lower rates, the two brine treatment options have lower costs. Life-cycle-analysis studies on extracted brine handling options suggest that a UIC well has a lower environmental impact than brine treatment. Both brine disposal options using brine treatment require removal of suspended solids from the extracted brine. The most suitable commercially available technology and the most promising emerging and innovative technology are recommended for implementation in Phase II. Though the challenges of this project are written specific to Decatur, every CO 2 storage site considering the use of brine extraction integrated with CO 2 storage will have similar, if not identical, technical and logistical challenges.

58 GEOSCIENCES↗
DOE OSTI2022

Using Natural Gas Liquids to Recover Unconventional Oil and Gas Resources (Final Report)

This document presents final technical findings for the project Using Natural Gas Liquids to Recover Unconventional Oil and Gas Resources (FE0031782). The project is part of the U.S. Department of Energy Oil and Gas Program to develop and advance technologies that can significantly improve the recovery efficiencies of unconventional oil and gas resources. The overall objective of this project is to improve the ultimate recovery from unconventional oil and gas (UOG) resources in the United States by developing a method for using unrefined natural gas liquids (NGLs) as treatment fluids to improve hydrocarbon production. Horizontal extended-lateral drilling coupled with high volume hydraulic fracturing has significantly increased production from UOG resources in the U.S. However, the recovery efficiency is low compared to the estimated oil and gas in place. Recent data indicate that less than 10% of the oil in the liquid-rich UOG reservoirs is produced. Alternative completion methods using NGLs could increase production (Battelle, 2016; Wan et al., 2013; Wan, 2013; Downey et al, 2021); however, field validation tests are needed to develop an approach that is economical, efficient, and compatible in the UOG setting to advance towards commercial deployment. This project aims to develop and field test a method to improve recovery of oil resources in UOG shale plays by using Y-Grade NGLs, or a similar combination of NGLs, as treatment fluids. Refined NGLs have been used as a hydraulic treatment fluid in UOG plays for decades and are shown to be particularly effective because their miscibility with oil allows oil to flow more freely; however, the use of Y-Grade (unrefined) NGLs has not been studied. The use of Y-Grade NGLs would be advantageous over refined NGLs because Y-Grade NGLs do not require infrastructure or investment in refining and are already being produced from many UOG reservoirs. The concept was tested and monitored in the field at a commercial well site owned by project partner Hopco, Ltd. The project team, which consists of multiple oil and gas operators, Linde Gas North America LLC (Linde) and the Ohio Division of Geological Survey (ODGS), has extensive experience with oil and gas production in the Appalachian basin and the ability to work together quickly to solve technical issues and research needs. A key part of the proposed work was the use of existing wells for field testing and monitoring. A total of four wells (three vertical and one horizontal) were available for this project. One of the vertical wells was utilized as the test well for the NGL treatment test. A nearby vertical well was used for microseismic monitoring. The remaining vertical well and the horizontal well provided a baseline for typical UOG production in the oil window of the Utica/Point Pleasant (UPP). Major technical tasks of the project include characterization of the geotechnical properties of the UPP with an emphasis on the field site; design and planning for the NGLs testing; field testing and monitoring; analysis and integration of field data; and economic and resource/reserve assessment. The Shoman monitoring well and Doughty NGL treatment well were successfully plugged back during September-October 2020 in preparation for treatment and monitoring. A nitrogen diagnostic fracture injection test (DFIT) was completed on the Doughty well on July 22, 2021, consisting of 133,000 scf (91 Bbl.) of nitrogen. A nitrogen foam frac was completed in the Utica-Point Pleasant interval on August 17, 2021 with funding from outside sources. A microseismic monitoring array was installed in the Shoman well and monitored microseismic activity during the Doughty well frac job. Y-Grade NGL injection commenced on August 26, 2021. A total of 215 Bbl. was injected but the job was shut down due to a small leak on the suction hose on the pump truck. The Y-Grade treatment resumed on August 27, 2021 and an additional 726 Bbl. of Y-Grade NGL was injected at a well head pressure of 3850 psi. Total volume of injected Y-Grade over the two days of injection was 941 Bbl. The well was shut-in for 17 days following injection to allow the Y-Grade NGLs to soak on the formation. Y-Grade treatment flow back commenced on 9/13/2021 on a weekly basis. Production data, including surface pressures, oil, nitrogen, natural gas, and flow times was measured and recorded. Periodic gas samples were collected and analyzed to determine composition of flowback gas. As of July 2022, the treatment well had produced 726 Bbl of oil and 2,888 mcf gas. In August 2022, tubing and packer in the well were removed and a pump was installed to enhance oil recovery. Currently, the operator is producing the well about 2 days a week for a few hours. The performance of the NGL treatment test was evaluated based on reservoir simulations of the treatment process, processing of well testing data, analysis of micro seismic monitoring data, and production data analysis. This analysis suggested that oil production in the small test would continue through 2025. Upscaling the treatment to a horizontal Utica Point-Pleasant well would allow more oil recovery, but the process would involve more investment, services, and operational support. An economic analysis was conducted for scenarios aimed at upscaling the NGL treatment process for more typical horizontal UPP wells in the Appalachian Basin.

02 PETROLEUM↗
DOE OSTI2024

Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring. We have deployed strainmeters and then conducted injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters (30 to 40m depth) was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 n/d over a few days. The signal at two strainmeters at shallow depth were consistent, although the magnitude of the horizontal strains were different reflecting the different radial directions from the well. The signal at a deep strainmeter deployed at reservoir depth was much different, with tensile vertical strains and compressive horizontal strains. These data can be interpreted by inverting poroelastic forward models developed using numerical and analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. This model is fast and can be inverted to estimate reservoir stiffness and geometry. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. The proxy model is periodically updated and refined using the finite element model to ensure accuracy. This approach significantly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. We have shown that the strain tensor in the caprock is sensitive to pressure in the reservoir, boundaries in the reservoir, and pressure in the caprock caused by leaks. These results indicate that coupling strain tensor data with inversion has the potential to help evaluate reservoirs during initial characterization, and to monitor them during the CO2 injection and storage process.

Murdoch, Larry↗
DOE OSTI2019

Impact of model improvements on 80 m wind speeds during the second Wind Forecast Improvement Project (WFIP2)

Abstract. During the second Wind Forecast Improvement Project (WFIP2; October 2015–March 2017, held in the Columbia River Gorge and Basin area of eastern Washington and Oregon states), several improvements to the parameterizations used in the High Resolution Rapid Refresh (HRRR – 3 km horizontal grid spacing) and the High Resolution Rapid Refresh Nest (HRRRNEST – 750 m horizontal grid spacing) numerical weather prediction (NWP) models were tested during four 6-week reforecast periods (one for each season). For these tests the models were run in control (CNT) and experimental (EXP) configurations, with the EXP configuration including all the improved parameterizations. The impacts of the experimental parameterizations on the forecast of 80 m wind speeds (wind turbine hub height) from the HRRR and HRRRNEST models are assessed, using observations collected by 19 sodars and three profiling lidars for comparison. Improvements due to the experimental physics (EXP vs. CNT runs) and those due to finer horizontal grid spacing (HRRRNEST vs. HRRR) and the combination of the two are compared, using standard bulk statistics such as mean absolute error (MAE) and mean bias error (bias). On average, the HRRR 80 m wind speed MAE is reduced by 3 %–4 % due to the experimental physics. The impact of the finer horizontal grid spacing in the CNT runs also shows a positive improvement of 5 % on MAE, which is particularly large at nighttime and during the morning transition. Lastly, the combined impact of the experimental physics and finer horizontal grid spacing produces larger improvements in the 80 m wind speed MAE, up to 7 %–8 %. The improvements are evaluated as a function of the model's initialization time, forecast horizon, time of the day, season of the year, site elevation, and meteorological phenomena. Causes of model weaknesses are identified. Finally, bias correction methods are applied to the 80 m wind speed model outputs to measure their impact on the improvements due to the removal of the systematic component of the errors.

17 WIND ENERGY↗
NASA NTRS1961

Theory of Electrostatic Fields in the Ionosphere at Equatorial Latitudes

The properties of the elongated electrostatic fields that are required to provide the coupling mechanism in the dynamo-motor concept of the E and F regions of the ionosphere are examined theoretically for the conditions that prevail in equatorial latitudes. The analysis is developed for an electrostatic field of arbitrary horizontal scale in a horizontally stratified partly ionized gas subject to an imposed magnetic field having the form of a parabolic arch over the equator. The anisotropic character and continuous variation with height of the conductivity are retained throughout , and numerical solutions are determined for the attenuation of the electric field with distance along the field line. The results are similar qualitatively to those found previously upon analysis of the corresponding problem for middle latitudes, but the attenuation of the electrostatic field with height is considerably greater. It develops, in particular, that the coupling between E and F regions is very small for fields having horizontal wavelengths of a few kilometers. It follows that the dynamo-motor concept could not be used to account for the presence of irregularities of this scale. On the other hand, it is found that almost all the attenuation occurs at heights near that of the E region. If an electrostatic field having a horizontal wavelength of a few kilometers could be produced at heights of, say 200 km or greater, it follows that it would be very elongated and extend with little change in amplitude from hemisphere to hemisphere. The results also indicate that effective coupling could be achieved for fields having horizontal wavelengths of the order of tens of kilometers, but the assumptions introduced to simplify the analysis may impair the quantitative reliability of the results for fields of this scale.

Briggs, Benjamin R.↗
NASA NTRS1972

Malfunction of Television Camera Monitor on Command Module

The command module television camera monitor exhibited loss of horizontal synchronization during the initial usage. This condition cleared and performance of the monitor was normal until the press conference telecast during the transearth coast phase. At that time, the monitor had the same horizontal synchronization problem reported during the initial usage. The horizontal hold control adjustment would not correct the condition. The monitor was turned off for approximately 5 minutes, then turned back on, after which the monitor's picture was normal. The loss of horizontal synchronization was most likely caused by either an intermittent condition in the precision voltage regulator circuit assembly in the low voltage power supply, or the shift of the stable range of the horizontal potentiometer setting when warming up after turn-on.

Source record↗
NASA NTRS1976

PECHCV, PECHFV, PEFHCV and PEFHFV: A set of atmospheric, primitive equation forecast models for the Northern Hemisphere, volume 3

As part of the SEASAT program of NASA, a set of four hemispheric, atmospheric prediction models were developed. The models, which use a polar stereographic grid in the horizontal and a sigma coordinate in the vertical, are: (1) PECHCV - five sigma layers and a 63 x 63 horizontal grid, (2) PECHFV - ten sigma layers and a 63 x 63 horizontal grid, (3) PEFHCV - five sigma layers and a 187 x 187 horizontal grid, and (4) PEFHFV - ten sigma layers and a 187 x 187 horizontal grid. The models and associated computer programs are described.

Wellck, R. E.↗
NASA NTRS1984

Longwall shearer tracking system

A tracking system for measuring and recording the movements of a longwall shearer vehicle includes an optical tracking assembly carried at one end of a desired vehicle path and a retroreflector assembly carried by the vehicle. Continuous horizontal and vertical light beams are alternately transmitted by means of a rotating Dove prism to the reflector assembly. A vertically reciprocating reflector interrupts the continuous light beams and converts these to discrete horizontal and vertical light beam images transmitted at spaced intervals along the path. A second rotating Dove prism rotates the vertical images to convert them to a second series of horizontal images while the first mentioned horizontal images are left unrotated and horizontal. The images are recorded on a film.

Poulsen, P. D.↗
NASA NTRS1986

Use of the Sun to determine pointing of ST radar beams

Verification of the beam pointing direction for ST (stratosphere troposphere) radars is a technically difficult problem. Consequently it is not usually done. For measurement of horizontal wind, the lack of precise knowledge of the beam pointing direction is usually of little consequence as any errors cause only a small uncertainty in the measured velocity. However, instantaneous vertical velocities are typically more than an order of magnitude less than horizontal velocities and average vertical velocities are more than two orders of magnitude less than average horizontal velocities. Hence small pointing errors for vertical beams can result in large errors due to contamination by horizontal winds. Experimental confirmation of pointing accuracy using the measured winds is difficult but has been achieved where horizontal and vertical winds are measured at the same site. Using the techniques described by Balsley and Riddle, 1984 and Riddle and Balsley, 1985, pointing accuracy was measured at Ponape, Federated States of Micronesia. The experiment is described in detail.

Riddle, A. C.↗
NASA NTRS1991

The 3-D description of vertical current sheets with application to solar flares

Following a brief review of the processes which have been suggested for explaining the occurrence of solar flares we suggest a new scenario which builds on the achievements of the previous suggestion that the current sheets, which develop naturally in 3-D cases with gravity from impacting independent magnetic structures (i.e., approaching current systems), do not consist of horizontal currents but are instead predominantly vertical current systems. This suggestion is based on the fact that as the subphotospheric sources of the magnetic field displace the upper photosphere and lower chromosphere regions, where plasma beta is near unity, will experience predominantly horizontal mass motions which will lead to a distorted 3-D configurations of the magnetic field having stored free energy. In our scenario, a vertically flowing current sheet separates the plasma regions associated with either of the subphotospheric sources. This reflects the balanced tension of the two stressed fields which twist around each other. This leads naturally to a metastable or unstable situation as the twisted field emerges into a low beta region where vertical motions are not inhibited by gravity. In our flare scenario the impulsive energy release occurs, initially, not by reconnection but mainly by the rapid change of the magnetic field which has become unstable. During the impulsive phase the field lines contort in such way as to realign the electric current sheet into a minimum energy horizontal flow. This contortion produces very large electric fields which will accelerate particles. As the current evolves to a horizontal configuration the magnetic field expands vertically, which can be accompanied by eruptions of material. The instability of a horizontal current is well known and causes the magnetic field to undergo a rapid outward expansion. In our scenario, fast reconnection is not necessary to trigger the flare, however, slow reconnection would occur continuously in the current layer at the locations of potential flaring. During the initial rearrangement of the field strong plasma turbulence develops. Following the impulsive phase, the final current sheet will experience faster reconnection which we believe responsible for the gradual phase of the flare. The reconnection will dissipate part of the current and will produce sustained and extended heating in the flare region and in the postflare loops.

Fontenla, Juan M.↗
NASA NTRS1992

Finite amplitude gravity waves in the Venus atmosphere generated by surface topography

A two-dimensional, fully nonlinear, nonhydrostatic, gravity wave model is used to study the evolution of gravity waves generated near the surface of Venus. The model extends from near the surface to well above the cloud layers. Waves are forced by applying a vertical wind at the bottom boundary. The boundary vertical wind is determined by the product of the horizontal wind and the gradient of the surface height. When wave amplitudes are small, the near-surface horizontal wind is the zonally averaged basic-state zonal wind, and the length scales of the forcing that results are characteristic of the surface height variation. When the forcing becomes larger and wave amplitudes affect the near-surface horizontal wind field, the forcing spectrum becomes more complicated, and a spectrum of waves is generated that is not a direct reflection of the spectrum of the surface height variation. Model spatial resolution required depends on the amplitude of forcing; for very nonlinear cases considered, vertical resolution was 250 m, and horizontal resolution was slightly greater than 1 km. For smaller forcing amplitudes, spatial resolution was much coarser, being 1 km in the vertical and about 10 km in the horizontal. Background static stability and mean wind are typical of those observed in the Venus atmosphere.

Young, R. E.↗
NASA NTRS1999

Thermal Infrared Spectroscopy from Mars Landers and Rovers: A New Angle on Remote Sensing

The MINUTES instrument of the Athena Precursor Experiment (APEX) on the Mars Surveyor 2001 lander mission will perform the first thermal infrared remote sensing observations from the surface of another planet. Experience gained from this experiment will be used to guide observations from identical instruments mounted on the Athena rovers, to be launched in 2003 and 2005. The utility of infrared spectrometers in determining the mineralogic composition of geologic surfaces from airborne and spaceborne platforms has been amply demonstrated. However, relatively little experience exists in using functionally similar instruments on the ground in the context of planetary science. What work has been done on this problem has mostly utilized field spectrometers that are designed to look down on nearby target rocks. While many Mini-TES observations will be made with this type of geometry, it is likely that other observations will be made looking horizontally at the more vertically-oriented facets of rock targets, to avoid spectral contamination from dust mantles. On rover missions, the Mini-TES may also be pointed horizontally at rocks several meters away, to determine if they are worthy of approaching for in situ observations and possible sample cacheing. While these observations will undoubtedly prove useful, there are important, and perhaps unappreciated, differences between horizontal-viewing, surface-based spectroscopy and the more traditional nadir-viewing, orbit or aircraft-based observations. Plans also exist to step the Mini-TES in a rastering motion to build hyperspectral scenes. Horizontal viewing hyperspectral cubes also possess unique qualities that call for innovative analysis techniques. The effect of viewing geometry: In thermal emission spectroscopy, regardless of whether an instrument is looking down on or horizontally at a target, the same basic equation governs the radiance reaching the sensor .

Moersch, J.↗