Search NASA⌕ Search

SEARCH · Search NASA

Results for “VERTICAL MIGRATION”

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 55 records · Page 3

Assisted migration in a warmer and drier climate: less climate buffering capacity, less facilitation and more fires at temperate latitudes?

Assisted tree migration has been proposed as a conceptual solution to mitigate lags in biotic responses to anthropogenic climate change. The rationale behind this concept is that tree species currently growing under warmer and drier climates will be more resistant and resilient to the new climatic conditions than tree species naturally growing in currently wetter and colder climates. However, we hypothesize that, by being more stress‐tolerant to warmer and drier conditions, translocated species should exhibit different functional attributes, which could induce important ecological and societal costs and overcome the desired benefits of maintaining wood production and other ecosystem services. We used principal component analysis (PCA) to analyze variation in seven traits of 106 tree and tall shrub species from contrasting latitudinal distributions in western North America and Europe to predict the potential functional changes of forest ecosystems due to the translocation of tree species from low to high latitudes. We show that species from both continents differed primarily by their position on the leaf economy spectrum (LES) and their size traits. Even though, in Europe, differences in LES were significantly correlated to species southern latitudinal positions, in both continents differences in size traits were significantly correlated to latitude. These results suggest that assisted migration by translocating more conservative species of shorter stature in currently cooler climates should decrease the buffering capacity of forest canopies, decrease facilitation for understory species, and increase wildfire risks, whose effects have the potential to accelerate climate warming through negative atmospheric feedback processes. As an alternative solution to assisted migration that may accelerate rather than mitigate climate change, we recommend that foresters gradually diversify the vertical structure and layering of the existing forest canopy to maintain a sustainable water cycle and energy balance between the soil, the tree and the atmosphere without increasing the wildfire risk.

Environmental Sciences & Ecology↗

Engineering of Grain Boundaries in CeO 2 Enabling Tailorable Resistive Switching Properties

Abstract Defect engineering in valence change memories aimed at tuning the concentration and transport of oxygen vacancies are studied extensively, however mostly focusing on contribution from individual extended defects such as single dislocations and grain boundaries. In this work, the impact of engineering large numbers of grain boundaries on resistive switching mechanisms and performances is investigated. Three different grain morphologies, that is, “random network,” “columnar scaffold,” and “island‐like,” are realized in CeO 2 thin films. The devices with the three grain morphologies demonstrate vastly different resistive switching behaviors. The best overall resistive switching performance is shown in the devices with “columnar scaffold” morphology, where the vertical grain boundaries extending through the film facilitate the generation of oxygen vacancies as well as their migration under external bias. The observation of both interfacial and filamentary switching modes only in the devices with a “columnar scaffold” morphology further confirms the contribution from grain boundaries. In contrast, the “random network” or “island‐like” structures result in excessive or insufficient oxygen vacancy concentration migration paths. The research provides design guidelines for grain boundary engineering of oxide‐based resistive switching materials to tune the resistive switching performances for memory and neuromorphic computing applications.

36 MATERIALS SCIENCE↗

Analysis of Hydrology and Interim Measure Performance in the Chromium Plume at Los Alamos National Laboratory - 20498

Hexavalent chromium, Cr(VI), is present in the regional aquifer beneath Los Alamos National Laboratory (LANL) at concentrations greater than the New Mexico groundwater standard (50 ppb). The source of the Cr is blowdown from a power plant where potassium dichromate was used as a corrosion inhibitor from 1956 to 1972. Blowdown comingled with effluent from another source and migrated approximately 3 miles down a canyon, primarily as surface flow, until reaching an alluvial infiltration zone. The infiltrating water moved vertically through a thick (approximately 900 ft), complex stratigraphic sequence of unsaturated tuffs, basalt, and unconsolidated alluvial deposits. Horizontal transport also took place within spatially limited perched groundwater horizons in the vadose zone (VZ). Breakthrough at the regional aquifer water table is thought to have occurred in multiple hydraulic windows or 'drip points' located approximately 3 miles from the initial release site at the power plant. The regional aquifer is highly heterogeneous and anisotropic, and displays unconfined behavior in an upper zone and confined behavior at depth. The Cr(VI) contamination has so far remained confined to the upper approximately 25 m of the aquifer. The ambient hydraulic gradient in the regional aquifer is toward the east/southeast, but the water table is extremely flat in the Cr plume area compared to surrounding areas. There is some local variability in flow direction, potentially related to effects of seasonal pumping from nearby Los Alamos county water-supply wells, which draw from the regional aquifer several hundred feet below the water table. An interim measure (IM) is currently operating at the site to achieve hydraulic control of plume migration using a series of extraction wells, ion-exchange treatment, and injection along the downgradient portion of the plume. The migration of Cr contamination in the aquifer has been previously evaluated using a calibrated probabilistic numerical model, which has also been used to design the spatial configuration of extraction and injection wells for the IM, and to identify locations for monitoring wells. The model is built using the Finite Element Heat and Mass (FEHM) transfer code, and the Model Analysis and Decision Support (MADS) software is used for analyses. To date, the portion of the IM system that has operated in a continuous manner has been very effective at hydraulic control. Modifications to the flow field associated with the influence of IM pumping and injection are evaluated here, with modeling constraints to be provided by recent observations of tracers injected into both monitoring wells and IM injection wells, as well as opportunistic observations of geochemical signatures of treated ion-exchange effluent injected into IM injection wells. The tracer and geochemical data provide insights into modified flow directions within the Cr plume area that cannot be otherwise obtained. Results of the modeling work are being used to inform adaptive management of IM operations and will be a valuable tool for development and refinement of long-term remediation strategies. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Vertical Transport of Water/Contaminant Due to Intermittent Water Influx on the Ground Surface of Vadose Zone - 20088

For movement of moisture and transport of contaminants in the vertical direction, two equations, the Richards equation and the advection-dispersion equation, are solved in a vadose zone for prescribed water and contaminant influx conditions at the ground surface. Two different boundary conditions are considered, one with a single peak of sinusoidal variation and four peaks of finite duration, each of which is sinusoidal in shape. The moisture migration is solved and then the contaminant transport, which is coupled to moisture movement, is solved. In the case of a single influx on the ground surface, at first the hydraulic head and moisture content show a peak. Then above the peak and below the ground surface, the peaks tend to die out, showing a tendency to return to the initial state. The vertical profile of concentration also shows a single peak, which soon becomes stagnant because downward water migration is faster than concentration profile movement. In the case of four intermittent array of peaks, the hydraulic head and moisture content profiles tend to merge and combine together. The concentration profiles show appearance of four peaks clearly. However, although the peak value of concentration remains quite constant, the values may vary with the sorption coefficient because of repeated influx of water and contaminant through the top boundary, the ground surface. The penetration depth decreases with an increasing sorption coefficient. (authors)

54 ENVIRONMENTAL SCIENCES↗

Influence of Fault Architecture on Induced Earthquake Sequence Evolution Revealed by High-Resolution Focal Mechanism Solutions

The increasing seismicity and improved seismic observation network in recent years provide an opportunity to explore factors that influence the triggering processes, spatiotemporal evolution, and maximum magnitude of induced sequences. We map the fault architecture and stress state of four induced sequences in Oklahoma to determine their influence on the seismicity. We systematically relocate the earthquakes and compute hundreds of focal mechanisms of small to medium events (1.0 < M < 5.1) using various techniques, including machine learning, for the Guthrie, Woodward, Cushing, and Fairview sequences. The detailed fault geometry and spatiotemporal evolution of seismicity and stress states reveal different dominant driving forces for each sequence. In Cushing and Fairview (largest event ≥M5.0), the main fault structures are near-vertical narrow strike-slip faults, with most of the small earthquake fault planes optimally oriented. The two sequences exhibit discontinuous temporal migration but strong earthquake self-driven rupture growth. In Guthrie and Woodward (largest event <M5.0), the two sequences show more complex diffuse fault structures with varying dipping angles along depth. The inverted focal mechanisms show a mix of strike-slip faulting and normal faulting in both sequences, and the normal faulting events are less optimally oriented than strike-slip events. The two sequences are dominated by continuous diffusive migration in time driven by pore pressure propagation. The above results suggest that fault architecture and stress state influence sequence evolution, major driving forces, and possibly maximum magnitude.

58 GEOSCIENCES↗

Seismic Monitoring at the Farnsworth CO2-EOR Field Using Time-Lapse Elastic-Waveform Inversion of 3D-3C VSP Data

During the Development Phase of the U.S. Southwest Regional Partnership on Carbon Sequestration, supercritical CO2 was continuously injected into the deep oil-bearing Morrow B formation of the Farnsworth Unit in Texas for Enhanced Oil Recovery (EOR). The project injected approximately 94 kilotons of CO2 to study geologic carbon storage during CO2-EOR. A three-dimensional (3D) surface seismic dataset was acquired in 2013 to characterize the subsurface structures of the Farnsworth site. Following this data acquisition, the baseline and three time-lapse three-dimensional three-component (3D-3C) vertical seismic profiling (VSP) data were acquired at a narrower surface area surrounding the CO2 injection and oil/gas production wells between 2014 and 2017 for monitoring CO2 injection and migration. With these VSP datasets, we inverted for subsurface velocity models to quantitatively monitor the CO2 plume within the Morrow B formation. We first built 1D initial P-wave (Vp) and S-wave (Vs) velocity models by upscaling the sonic logs. We improved the deep region of the Vp and Vs models by incorporating the deep part of a migration velocity model derived from the 3D surface seismic data. We improved the shallow region of 3D Vp and Vs models using 3D traveltime tomography of first arrivals of VSP downgoing waves. We further improved the 3D baseline velocity models using elastic-waveform inversion (EWI) of the 3D baseline VSP upgoing data. Our advanced EWI method employs alternative tomographic and conventional gradients and total-variation-based regularization to ensure the high-fidelity updates of the 3D baseline Vp and Vs models. We then sequentially applied our 3D EWI method to the three time-lapse datasets to invert for spatiotemporal changes of Vp and Vs in the reservoir. Our inversion results reveal the volumetric changes of the time-lapse Vp and Vs models and show the evolution of the CO2 plume from the CO2 injection well to the oil/gas production wells.

42 ENGINEERING↗

Multi-Scale Seismic Measurements for Site Characterization and CO2 Monitoring in an Enhanced Oil Recovery/Carbon Capture, Utilization, and Sequestration Project, Farnsworth Field, Texas

To address the challenges of climate change, significantly more geologic carbon sequestration projects are beginning. The characterization of the subsurface and the migration of the plume of supercritical carbon dioxide are two elements of carbon sequestration that can be addressed through the use of the available seismic methods in the oil and gas industry. In an enhanced oil recovery site in Farnsworth, TX, we employed three separate seismic techniques. The three-dimensional (3D) surface seismic survey required significant planning, design, and processing, but produces both a better understanding of the subsurface structure and a three-dimensional velocity model, which is essential for the second technique, a timelapse vertical seismic profile, and the third technique, cross-well seismic tomography. The timelapse 3D Vertical Seismic Profile (3D VSP) revealed both significant changes in the reservoir between the second and third surveys and geo-bodies that may represent the extent of the underground carbon dioxide. The asymmetry of the primary geo-body may indicate the preferential migration of the carbon dioxide. The third technique, cross-well seismic tomography, suggested a strong correlation between the well logs and the tomographic velocities, but did not observe changes in the injection interval.

Energy & Fuels↗

Time-Lapse Gravity Monitoring of CO2 Migration Based on Numerical Modeling of a Faulted Storage Complex

In this study, the performance of both surface and borehole time-lapse gravity monitoring to detect CO2 leakage from a carbon storage site is evaluated. Several hypothetical scenarios of CO2 migration in a leaky fault, and thief zones at different depths at the Kimberlina site (California, USA) constitute the basis of the approach. The CO2 displacement is simulated using the TOUGH2 simulator applied to a detailed geological model of the site. The gravity responses to these CO2 plumes are simulated using forward modeling with sensors at ground surface and in vertical boreholes. Results of inversion on one scenario are also presented. The surface-based gravity responses obtained for the different leakage scenarios demonstrate that leakage can be detected at the surface in all the scenarios but the time to detection is highly variable (10 to 40 years) and dependent on the detection threshold considered. Borehole measurements of the vertical component of gravity provide excellent constrains in depth when they are located in proximity of the density anomaly associated with the presence of CO2, thus discriminating multiple leaks in different thief zones. Joint inversion of surface and borehole data can bring valuable information of the occurrence of leakages and their importance by providing a reasonable estimate of mass of displaced fluids. This study demonstrates the importance of combining multiphase flow simulations with gravity modeling in order to define if and when gravity monitoring would be applicable at a given storage site.

Time-lapse gravity monitoring, leak detection, CCS↗

Strain age cracking in a simulated heat affected zone of Inconel 740H laser-powder bed fusion components

Inconel® 740H components produced via laser-powder-bed fusion (L-PBF) additive manufacturing are ideal for supercritical CO2 primary heat exchangers. Arc welding is often needed, and subsequent post-weld heat treatment aging at 790–840 °C is required to improve strength via gamma prime (γ’) precipitation; however, strain age cracking (SAC) can occur in the heat affected zone (HAZ) during this process. This study uses stress relaxation testing at 800 °C on simulated HAZ specimens from vertically and horizontally built L-PBF IN740H to assess SAC susceptibility across heating rates of 40–3480 °C/h and weld-induced strains of 3–10 %. Vertical builds require greater mechanical energy input and exhibit longer times to fracture than horizontal builds, and time to fracture occurs sooner at slower heating rates. Creep voids were observed in γ’-denuded regions along grain boundaries, and cracking propagated along migrated grain boundaries and at interfaces between elongated secondary phases (γ’ or carbides) and the γ-matrix.

14 SOLAR ENERGY↗

Geophysical monitoring using active seismic techniques at the Citronelle Alabama CO 2 storage demonstration site

Between August 2012 and September 2014, about 114,000 metric tonnes of CO 2 was captured from the coal-fired Plant Barry Power Station at Bucks Alabama and injected into the Paluxy Formation above the oil pool in the southeast unit of the Citronelle Oilfield. Various monitoring methods were deployed at land surface and in project wells to measure system performance, comply with permit requirements and test new and innovative monitoring tools. The monitoring program relied heavily on active seismic methods for subsurface imaging of geologic structure and time-lapse seismic techniques to track the CO 2 migration in the injection interval. Both conventional geophone/hydrophone and fiber-optic based Distributed Acoustic Sensing (DAS) arrays were deployed and tested, allowing a side by side comparison of the equipment and techniques. Geophysical imaging of the subsurface was successful using DAS in the offset vertical seismic profile (OVSP) survey configuration. A high resolution OVSP image of the subsurface was obtained in 2014 with DAS, which exceeded project expectations in comparison to a lower resolution image obtained in 2012 using a conventional 80-level geophone array. A time-lapse image of the redistribution of CO 2 after injection ended in September 2014 was obtained with two DAS OVSP surveys from June 2014 and December 2015, thus successfully demonstrating its proof-of-concept. Unfortunately, a pre-injection baseline survey with DAS, which was in its initial stage of technology development in 2012, did not have sufficient quality for use, making it difficult to interpret the acquired DAS time-lapse difference. Additional research in this area has since demonstrated the utility of time-lapse DAS OVSP. DAS data were also acquired during a cross-well seismic survey conducted in 2014. Unfortunately, the DAS technique was not success in the cross-well survey configuration because the system noise level was too high in the crosswell frequency output range (100–1200 Hz) of the piezoelectric source (increasing by a factor of ten compared to VSP frequency band). Additionally, the cross-well geometry causes sub-horizontal (broadside) incidence on the vertical DAS fiber cable, which is known to be problematic. Current research is focused on improving the DAS cable response to broadside acoustic energy. Time-lapse seismic surveys using commercially available conventional arrays were also acquired. In contrast to the DAS acquired data, the cross-well seismic results obtained with the conventional array was highly successful and clearly showed the CO 2 remained in zone at the end of injection. Time-lapse differencing of the OSVP surveys acquired with the conventional arrays proved to be inconclusive. Finally, changes in wellbore conditions between surveys and unavoidable changes in equipment (the array used for the baseline survey was retired) affected data quality, making it difficult to interpret the OVSP results.

58 GEOSCIENCES↗

Enhancing cathode composites with conductive alignment synergy for solid-state batteries

Enhancing transport and chemomechanical properties in cathode composites is crucial for the performance of solid-state batteries. Our study introduces the filler-aligned structured thick (FAST) electrode, which notably improves mechanical strength and ionic/electronic conductivity in solid composite cathodes. The FAST electrode incorporates vertically aligned nanoconducting carbon nanotubes within an ion-conducting polymer electrolyte, creating a low-tortuosity electron/ion transport path while strengthening the electrode’s structure. This design not only mitigates recrystallization of the polymer electrolyte but also establishes a densified local electric field distribution and accelerates the migration of lithium ions. The FAST electrode showcases outstanding electrochemical performance with lithium iron phosphate as the active material, achieving a high capacity of 148.2 milliampere hours per gram at 0.2 C over 100 cycles with substantial material loading (49.3 milligrams per square centimeter). This innovative electrode design marks a remarkable stride in addressing the challenges of solid-state lithium metal batteries.

Science & Technology - Other Topics↗

Unraveling the Crystallization Kinetics of 2D Perovskites with Sandwich-Type Structure for High-Performance Photovoltaics

2D perovskite solar cells with high stability and high efficiency have attracted significant attention. A systematical static and dynamic structure investigation is carried out to show the details of 2D morphology evolution. A dual additive approach is used, where the synergy between an alkali metal cation and a polar solvent leads to high-quality 2D perovskite films with sandwich-type structures and vertical phase segregation. Such novel structure can induce high-quality 2D slab growth and reduce internal and surface defects, resulting in a high device efficiency of 16.48% with enhanced continuous illumination stability and improved moisture (55-60%) and thermal (85°C) tolerances. Transient absorption spectra reveal the carrier migration from low n to high n species with different kinetics. An [PbI 6 ] 4- octagon coalescence transformation mechanism coupled with metal and organic cations wrapped is proposed. By solvent vapor annealing, a recrystallization and reorientation of the 2D perovskite slabs occurs to form an ideal structure with improved device performance and stability.

36 MATERIALS SCIENCE↗

Interfacial intermixing and anti-phase boundaries in GaP/Si(0 0 1) heterostructures

Epitaxial GaP/Si heterostructures grown by molecular beam epitaxy (MBE) and migration-enhanced epitaxy (MEE) have been studied, primarily using aberration-corrected scanning transmission electron microscopy (AC-STEM). Atomically-resolved structure images, which are sensitive to atomic-column intensity revealed the detailed geometry of antiphase boundaries that were both parallel and also inclined to the growth direction. The polar/non-polar GaP/Si interfaces were neither atomically flat laterally nor abrupt vertically. Measurements of intensity profiles using both bright-field and dark-field AC-STEM images, as well as chemically-sensitive g 002 dark-field imaging, indicated substantial interfacial intermixing, which increased significantly from ~1.3 nm (MEE growth at 440 °C) to ~2.1 nm (MBE growth at 600 °C). The finite interface width will impact theoretical predictions of charge imbalance and strong electric fields across the heterointerface.

36 MATERIALS SCIENCE↗

Photoelectron Spectroscopy and Theoretical Study on Monosolvated Cyanate Analogue Clusters ECX-· Sol (ECX- = NCSe-, AsCSe- and AsCS-; Sol = H2O, CH3CN)

Six monosolvated cyanate analogue clusters ECX-?Sol (ECX– = NCSe-, AsCSe-, and AsCS-; Sol = H2O and CH3CN) were investigated using negative ion photoelectron spectroscopy (NIPES). NIPES experiments show that these clusters possess overall similar spectra compared to their respective isolated ECX- anions but shift to higher electron binding energy with solvent CH3CN stabilizing the excess electrons slightly more than H2O. For the ECX-·H2O series, vertical detachment energies and their increments relative to the bare species are measured to be 3.700 / 0.370, 3.085 / 0.425, 3.085 / 0.430 eV for NCSe-, AsCSe- and AsCS-, respectively, while the corresponding values in ECX-·CH3CN series are 3.835 / 0.505, 3.145 / 0.475, and 3.135 / 0.480 eV. Ab initio electronic structure calculations indicate the excess charges located at the terminal N and Se atoms in NCSe- and migrated to the central C atom in AsCSe- and AsCS-. For NCSe-, the solvation is driven by the interactions with the two negatively charged terminal ends, while for AsCSe- and AsCS-, the solvation revolves around the interactions with the central C atom, where all the excess negative charge is concentrated. Two nearly degenerate isomers for NCSe-·H2O are identified, one forming a single strong N…H-O hydrogen bond (HB) and the other featuring a bidentate HB with two hydroxyl H atoms pointing to N and Se ends. In contrast, the negative central C atom in AsCSe-/AsCS- allows the formation of a bifurcated HB with H2O. Similar effects are observed for the acetonitrile case, in which the three H atoms of methyl group interact with the two negatively charged terminal ends in NCSe-, while prefer binding to the central negative carbon atom in AsCSe-/AsCS-. The different binding motifs derived in this work may suggest distinctly different solvation properties in NCSe- versus AsCSe-/AsCS-, with the former anion leading to asymmetric solvation at the N end of the solute, while the latter species creating more ‘isotropic’ solvation around the central C equatorial plane. This work was supported by U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences, and performed using EMSL, a national scientific user facility sponsored by DOE’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory, which is operated by Battelle Memorial Institute for the DOE. The theoretical calculations were conducted on EMSL’s “Cascade” Supercomputer.

Yuan, Qinqin↗

Experiments and Modeling of Proppant Embedment and Fracture Conductivity for the Caney Shale, Oklahoma, USA

ABSTRACT: The ultimate aim of hydraulic fracturing is to have a long and conductive flow path that extends from the wellbore into the formation. The effective fracture length is part of a hydraulically propped fracture which contributes to production. The difficulty in achieving economical production targets from shale reservoirs is at the forefront in many exploration companies. Fracture conductivity loss is related to; proppant embedment under depletion, proppant crushing, damage as a result of fracturing fluid, fines migration and proppant-pack permeability-damage are some of the factors that contribute to production decline after hydraulic fracturing in shale reservoirs. The Caney Shale is a calcareous organic-rich mudrock. Various studies have investigated the effect that clay on shale well productivity, however, there is currently no literature on the Caney shale in relation to horizontal wells; all available literature exists in vertical wells as well as on formations of the Caney that are shallow in comparison to an emerging play which is twice the depth. In this paper we investigate stress-dependent fracture conductivity of proppant-filled fractures and proppant embedment in Caney shale through laboratory and modeling studies. API fracture conductivity tests were conducted using 2% KCl on five locations within the Caney shale that consisted of selecting three brittle(reservoir) zones and two ductile zones. Confining pressures range from 1,000 psi to 12,000 psi at 210°F. Conductivity, permeability as well as embedment were measured during the test. Our experimental results have confirmed that improved fracture conductivity is attributed to; proppant size, the increase in porosity of the proppant pack, closure pressure changes and the reduction in fracture conductivity are a function of many factors such as fracture closure stress. The findings from this study could help the stimulation design by providing new insights into the critical factors that are to be determined to facilitate the choice of proppants as well as fracturing fluids for long term production and recovery from shale reservoirs. 1 INTRODUCTION The development of low permeability formations, like shales, has been aided by hydraulic fracturing of horizontal wells (Radonjic et al., 2020). Hydraulic fracturing fluid is injected at a high pressure to induce tensile fractures that can link to and stimulate natural fractures (Katende et al., 2021a,b). Preserving adequate conductivity in hydraulic fractures over the life of the wells is required for economic production; nevertheless, conserving such conductivity can be difficult in some circumstances, particularly in soft, clay-rich formations (Wang et al., 2021). Proppant particles help to keep the fractures open when the pumping stops and the fracturing fluid returns to the wellbore, producing one or more propped hydraulic fractures of varying length, breadth, and height (Katende et al., 2021a). The proppant pack within the hydraulic fracture boosts well output by providing a greater permeability flowpath for hydrocarbons (Duenckel et al., 2016). Proppant in the fracture is under complicated stress conditions, and the interplay between the rock formation and the proppant pack has a significant impact on proppant-pack permeability (Karazincir et al., 2019). Proppant may be embedded (Katende et al., 2021a) in the rock or crushed into small pieces if the proppant size and strength characteristics are not specified appropriately, resulting in a loss in proppant-pack permeability and fracture aperture, and consequently a fall in well output.

Katende, A.↗

Physical Properties and Gas Hydrate at a Near-Seafloor Thrust Fault, Hikurangi Margin, New Zealand

The Papaku Fault Zone, drilled at International Ocean Discovery Program (IODP) Site U1518, is an active splay fault in the frontal accretionary wedge of the Hikurangi Margin. In logging-while-drilling data, the 33-m-thick fault zone exhibits mixed modes of deformation associated with a trend of downward decreasing density, P-wave velocity, and resistivity. Methane hydrate is observed from ~30 to 585 m below seafloor (mbsf), including within and surrounding the fault zone. Hydrate accumulations are vertically discontinuous and occur throughout the entire logged section at low to moderate saturation in silty and sandy centimeter-thick layers. In this paper, we argue that the hydrate distribution implies that the methane is not sourced from fluid flow along the fault but instead by local diffusion. This, combined with geophysical observations and geochemical measurements from Site U1518, suggests that the fault is not a focused migration pathway for deeply sourced fluids and that the near-seafloor Papaku Fault Zone has little to no active fluid flow.

58 GEOSCIENCES↗

A Multi-variable, Multi-value Hydrodynamic Modeling Approach to Support Marine Energy and Coastal Resilience Applications

Funded by the Marine Energy Seedling Laboratory Program, the study reported on here aimed to develop a multi-variable, multi-value hydrodynamic modeling approach based on the Finite-Volume Community Ocean Model (FVCOM) to support a wide range of marine renewable energy and coastal resilience applications. We migrated the marine hydrokinetic energy module to simulate tidal energy extraction by tidal turbines under wave-current interactions. We explored the particle tracking module in FVCOM and applied it to Sequim Bay, Washington, to illustrate the surface current field and mimic the pollutant transport near the surface. A tracer-water age module was added to FVCOM to examine water exchange through the entrance. Lastly, a user-specified vertical coordinate option was added to FVCOM to improve simulation of the effect of floating turbines on surface currents and scalar transport. The above-described model development has been tested using a simplified model domain and demonstrated in real-world coastal systems. The outcome of this project will provide a useful modeling code and approach for marine renewable energy and coastal resilience research and applications.

16 TIDAL AND WAVE POWER↗

Wellbore Stability and Mud Loss Management in Geothermal Drilling: Optimizing Mud Weight to Mitigate Tensile Wellbore Fracturing at The Geysers, California

As part of a U.S. Department of Energy (DOE) Geothermal Technologies Office-funded initiative, Geysers Power Company, LLC, a subsidiary of Calpine Corporation, has been working to enhance drilling performance at the world’s largest geothermal field, The Geysers, in northern California. In a recent drilling operation of the GDC-36 well, excessive mud losses were encountered, initially addressed through repeated but largely ineffective cement plugging. Ultimately, the most effective strategy was to drill blind through the loss zones, made feasible by the high rate of penetration (ROP) achieved with PDC bits, allowing significant progress before the mud tanks were depleted and water-sensitive argillic formation layers could collapse. In response to these challenges, the project team explored alternative methods to minimize downtime and risks associated with cement plugging and continuous mud loss and to contemplate the driving mechanisms for the losses. Wellbore imaging using Formation MicroImager (FMI) and Ultrasonic Borehole Imager (UBI) tools revealed longitudinal tensile fractures, which were attributed to mud weights exceeding the minimum circumferential stress resulting from the native stress field and formation pressure. This study examines the mud losses encountered and leverages wellbore imaging data to understand the mechanisms behind mud induced tensile fracturing in specific rock facies. Understanding fracture behavior across different lithologies is crucial, as fractures within the reservoir can enhance steam migration throughout the system. The reservoir at The Geysers lies within the Mesozoic Franciscan Assemblage, a tectonic mélange formed by subduction. It consists of metamorphosed turbidite sandstone (greywacke) and mudstone (argillite), oceanic upper crust (including greenstone and chert), and serpentinized ultramafic rocks - each exhibiting distinct geomechanical fracturing properties. The structural fabric of the Franciscan Assemblage was shaped by low-angle Mesozoic thrust faulting and later overprinted by sub-vertical strike-slip structures related to the Pacific-North American plate boundary. A wellbore stability model was developed using core measurements and logs to simulate fracturing scenarios during drilling under varying stress conditions. These simulations guided the development of an optimized mud weight management strategy that should enable adaptive adjustments during drilling, reducing the likelihood of tensile fracturing and mud losses, ultimately improving operational efficiency.

15 GEOTHERMAL ENERGY↗