Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model reduction”

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,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Multi-dimensional simulations of the expanding supernova remnant of SN 1987A

The expanding remnant from SN 1987A is an excellent laboratory for investigating the physics of supernovae explosions. There is still a large number of outstanding questions, such as the reason for the asymmetric radio morphology, the structure of the pre-supernova environment, and the efficiency of particle acceleration at the supernova shock. We explore these questions using three-dimensional simulations of the expanding remnant between days 820 and 10,000 after the supernova. We combine a hydrodynamical simulation with semi-analytic treatments of diffusive shock acceleration and magnetic field amplification to derive radio emission as part of an inverse problem. Simulations show that an asymmetric explosion, combined with magnetic field amplification at the expanding shock, is able to replicate the persistent one-sided radio morphology of the remnant. We use an asymmetric Truelove and McKee progenitor with an envelope mass of 10 M {sub ☉} and an energy of 1.5 × 10{sup 44} J. A termination shock in the progenitor's stellar wind at a distance of 0.''43-0.''51 provides a good fit to the turn on of radio emission around day 1200. For the H II region, a minimum distance of 0.''63 ± 0.''01 and maximum particle number density of (7.11 ± 1.78) × 10{sup 7} m{sup –3} produces a good fit to the evolving average radius and velocity of the expanding shocks from day 2000 to day 7000 after explosion. The model predicts a noticeable reduction, and possibly a temporary reversal, in the asymmetric radio morphology of the remnant after day 7000, when the forward shock left the eastern lobe of the equatorial ring.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Electron cyclotron heating can drastically alter reversed shear Alfvén eigenmode activity in DIII-D through finite pressure effects

Here, a recent DIII-D experiment investigating the impact of electron cyclotron heating (ECH) on neutral beam driven reversed shear Alfvén eigenmode (RSAE) activity is presented. The experiment includes variations of ECH injection location and timing, current ramp rate, beam injection geometry (on/off-axis), and neutral beam power. Essentially all variations carried out in this experiment were observed to change the impact of ECH on AE activity significantly. In some cases, RSAEs were observed to be enhanced with ECH near the off-axis minimum in magnetic safety factor (${{q}_{\min}}$ ), in contrast to the original DIII-D experiments where the modes were absent when ECH was deposited near ${{q}_{\min}}$ . It is found that during intervals when the geodesic acoustic mode (GAM) frequency at ${{q}_{\min}}$ is elevated and the calculated RSAE minimum frequency, including contributions from thermal plasma gradients, is very near or above the nominal TAE frequency (f TAE ), RSAE activity is not observed or RSAEs with a much reduced frequency sweep range are found. This condition is primarily brought about by ECH modification of the local electron temperature (T e ) which can raise both the local T e at ${{q}_{\min}}$ as well as its gradient. A q-evolution model that incorporates this reduction in RSAE frequency sweep range is in agreement with the observed spectra and appears to capture the relative balance of TAE or RSAE-like modes throughout the current ramp phase of over 38 DIII-D discharges. Detailed ideal MHD calculations using the NOVA code show both modification of plasma pressure and pressure gradient at ${{q}_{\min}}$ play an important role in modifying the RSAE activity. Analysis of the ECH injection near the ${{q}_{\min}}$ case where no frequency sweeping RSAEs are observed shows the typical RSAE is no longer an eigenmode of the system. What remains is an eigenmode with poloidal harmonic content reminiscent of the standard RSAE, but absent of the typical frequency sweeping behavior. The remaining eigenmode is also often strongly coupled to gap TAEs. Analysis with the non-perturbative gyro fluid code TAEFL confirms this change in RSAE activity and also shows a large drop in the resultant mode growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

MHz gravitational waves from short-term anisotropic inflation

We reveal the universality of short-term anisotropic inflation. As a demonstration, we study inflation with an exponential type gauge kinetic function which is ubiquitous in models obtained by dimensional reduction from higher dimensional fundamental theory. It turns out that an anisotropic inflation universally takes place in the later stage of conventional inflation. Remarkably, we find that primordial gravitational waves with a peak amplitude around 10{sup −26}∼10{sup −27} are copiously produced in high-frequency bands 10 MHz∼100 MHz. If we could detect such gravitational waves in future, we would be able to probe higher dimensional fundamental theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The completed SDSS-IV extended baryon oscillation spectroscopic survey: geometry and growth from the anisotropic void–galaxy correlation function in the luminous red galaxy sample

ABSTRACT We present an analysis of the anisotropic redshift-space void–galaxy correlation in configuration space using the Sloan Digital Sky Survey extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 luminous red galaxy (LRG) sample. This sample consists of LRGs between redshifts 0.6 and 1.0, combined with the high redshift z > 0.6 tail of the Baryon Oscillation Spectroscopic Survey Data Release 12 CMASS sample. We use a reconstruction method to undo redshift-space distortion (RSD) effects from the galaxy field before applying a watershed void-finding algorithm to remove bias from the void selection. We then perform a joint fit to the multipole moments of the correlation function for the growth rate fσ8 and the geometrical distance ratio DM/DH, finding $f\sigma _8(z_\rm {eff})=0.356\pm 0.079$ and $D_M/D_H(z_\rm {eff})=0.868\pm 0.017$ at the effective redshift $z_\rm {eff}=0.69$ of the sample. The posterior parameter degeneracies are orthogonal to those from galaxy clustering analyses applied to the same data, and the constraint achieved on DM/DH is significantly tighter. In combination with the consensus galaxy BAO and full-shape analyses of the same sample, we obtain fσ8 = 0.447 ± 0.039, DM/rd = 17.48 ± 0.23, and DH/rd = 20.10 ± 0.34. These values are in good agreement with the ΛCDM model predictions and represent reductions in the uncertainties of $13{{\ \rm per\ cent}}$, $23{{\ \rm per\ cent}}$, and $28{{\ \rm per\ cent}}$, respectively, compared to the combined results from galaxy clustering, or an overall reduction of 55 per cent in the allowed volume of parameter space.

Nadathur, Seshadri↗
DOE OSTI2025

Quantum error mitigation for Fourier moment computation

Hamiltonian moments in Fourier space—expectation values of the unitary evolution operator under a Hamiltonian at different times—provide a convenient framework to understand quantum systems. They offer insights into the energy distribution, higher-order dynamics, response functions, correlation information, and physical properties. This paper focuses on the computation of Fourier moments within the context of a nuclear effective field theory on superconducting quantum hardware. The study integrates echo verification and noise renormalization into Hadamard tests using control reversal gates. These techniques, combined with purification and error suppression methods, effectively address quantum hardware decoherence. The analysis, conducted using noise models, reveals a significant reduction in noise strength by two orders of magnitude. Moreover, quantum circuits involving up to 266 gates over five qubits demonstrate high accuracy under these methodologies when run on IBM superconducting quantum devices. Published by the American Physical Society 2025

Kiss, Oriel (ORCID:0000000174613342)↗
DOE OSTI2013

Search for gamma-ray spectral lines with the Fermi Large Area Telescope and dark matter implications

Weakly interacting massive particles (WIMPs) are a theoretical class of particles that are excellent dark matter candidates. WIMP annihilation or decay may produce essentially monochromatic γ rays detectable by the Fermi Large Area Telescope (LAT) against the astrophysical γ -ray emission of the Galaxy. We have searched for spectral lines in the energy range 5–300 GeV using 3.7 years of data, reprocessed with updated instrument calibrations and an improved energy dispersion model compared to the previous Fermi-LAT Collaboration line searches. We searched in five regions selected to optimize sensitivity to different theoretically motivated dark matter density distributions. We did not find any globally significant lines in our a priori search regions and present 95% confidence limits for annihilation cross sections of self-conjugate WIMPs and decay lifetimes. Our most significant fit occurred at 133 GeV in our smallest search region and had a local significance of 3.3 standard deviations, which translates to a global significance of 1.5 standard deviations. We discuss potential systematic effects in this search, and examine the feature at 133 GeV in detail. We find that the use both of reprocessed data and of additional information in the energy dispersion model contributes to the reduction in significance of the linelike feature near 130 GeV relative to significances reported in other works. We also find that the feature is narrower than the LAT energy resolution at the level of 2 to 3 standard deviations, which somewhat disfavors the interpretation of the 133 GeV feature as a real WIMP signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Universal Collaboration Strategies for Signal Detection: A Sparse Learning Approach

Here, this paper considers the problem of high-dimensional signal detection in a large distributed network whose nodes can collaborate with their one-hop neighboring nodes (spatial collaboration). We assume that only a small subset of nodes communicate with the fusion center (FC). We design optimal collaboration strategies which are universal for a class of deterministic signals. By establishing the equivalence between the collaboration strategy design problem and sparse principal component analysis (PCA), we solve the problem efficiently and evaluate the impact of collaboration on detection performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Optimal Energy Dispatch of Distributed PVs for the Next Generation of Distribution Management Systems

Advanced Distribution Management Systems (ADMS) are being widely adopted by electric utilities for managing and optimizing the operations of their distribution systems. Distributed photovoltaic (DPV) systems with smart inverters can be controlled to adjust active power and reactive power outputs, and they are envisioned to become a part of (centrally or distributed) controllable assets managed by the ADMS for optimizing grid operations. This paper proposes an optimal energy dispatch strategy controlling DPV systems for regulating distribution voltages and achieving conservation voltage reduction. A convex optimization model is proposed with the use of linearized power flow, and the gradient projection algorithm is used to solve the optimal active power and reactive power outputs of smart inverters. The proposed optimal energy dispatch is implemented using an open-source ADMS platform, and simulation results have demonstrated the effectiveness of the proposed approach on improving distribution grid operations.

42 ENGINEERING↗
DOE OSTI2018

A new family of transcriptional regulators of tungstoenzymes and molybdate/tungstate transport

Bacterial genes for molybdenum-containing and tungsten-containing enzymes are often differentially regulated depending on the metal availability in the environment. Here in this paper, we describe a new family of transcription factors with an unusual DNA-binding domain related to excisionases of bacteriophages. These transcription factors are associated with genes for various molybdate and tungstate-specific transporting systems as well as molybdo/tungsto-enzymes in a wide range of bacterial genomes. We used a combination of computational and experimental techniques to study a member of the TF family, named TaoR (for tungsten-containing aldehyde oxidoreductase regulator). In Desulfovibrio vulgaris Hildenborough, a model bacterium for sulfate reduction studies, TaoR activates expression of aldehyde oxidoreductase aor and represses tungsten-specific ABC-type transporter tupABC genes under tungsten-replete conditions. TaoR binding sites at aor promoter were identified by electrophoretic mobility shift assay and DNase I footprinting. We also reconstructed TaoR regulons in forty five Deltaproteobacteria by comparative genomics approach and predicted target genes for TaoR family members in other Proteobacteria and Firmicutes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

WE-H-BRA-02: Radiobiological Modeling of Tumor Control Probability (TCP) and Radiation-Induced Pneumonitis (RP) for Lung Cancer Patients Treated with Monte Carlo-Based Lung SBRT

Purpose: To present radiobiological modeling of TCP using tumor size-adjusted BED(s-BED)and PTV(D99) to lung SBRT patients treated with X-ray Voxel Monte Carlo(XVMC) algorithm, apply parameterized Lyman-NTCP model to predict grade-2 RP and subsequently, compare with clinical outcomes/observations. Methods: Dosimetric parameters and clinical follow-up for XVMC-based lung-SBRT patients were retrospectively evaluated. Patients were treated at Novalis-TX with hybrid(2 non-coplanar partial-arcs plus 3–6 static-beams)plan using HD-MLC/6MV-SRS-beam.For TCP,s-BED modelling was utilized: TCP=EXP[sBED-TCD50]/k/(1.0+EXP[sBED-TCD50]/k), where k=31Gy corresponding to TCD50=0Gy and s-BED was defined as BED10 minus 10 times the tumor diameter(in centimeters)by Ohri et al.(IJROBP,2012). For 2-yr local-control, we used more-realistic MC-computed PTVD99 as a predictive parameter, s-BED(D99).Due to relatively shorter median follow-up interval(12-months),Kaplan-Meier curves were generated to estimate 2-yr observed local-control and compared to predicted-rate by TCP modeling. For NTCP, we employed parameterized Lyman-NTCP model utilizing normal-lung DVH and α/β=3Gy fitted to predict grade-2 RP after lung-SBRT. Results: Total 108 patients (137 tumors) treated for 35–70Gy in 3–5 fractions, either primary-lung(n=74)or metastatic-lung(n=53)tumors were included.F or the given prescription dose with MC-computed MUs, 2-yr local-control rates with s-BED(D99) was 87±8%. Kaplan-Meier generated observed local-control rate at 2-yr was 87.5%,suggesting that PTV(D99) could be a potential predictor (p-value=0.38). Observed vs predicted TCP for primary-lung tumors and metastatic tumors were 97% vs 88±7% and 94% vs 86±9%.NTCP model predicted well for symptomatic-RP with predicted vs observed (3±5% vs 2%). Radiographic and clinically significant RP was observed in 13% and 2% of patients. Higher rates of radiographic change were observed in patients who received >50Gy compared to ≤50Gy(24% vs 10%). Conclusion: Utilizing MC-computed PTVD99, our TCP results were well correlated with clinical outcome. The predicted grade-2 RP rate was comparable to clinical observations. Clinical application of these radiobiological models may potentially allow for target dose escalation and/or lung-toxicity reduction. Further validation of these radiobiological models with longer follow up interval for large cohorts of lung-SBRT patients is anticipated.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Explainable and trustworthy artificial intelligence for correctable modeling in chemical sciences

Data science has primarily focused on big data, but for many physics, chemistry, and engineering applications, data are often small, correlated and, thus, low dimensional, and sourced from both computations and experiments with various levels of noise. Typical statistics and machine learning methods do not work for these cases. Expert knowledge is essential, but a systematic framework for incorporating it into physics-based models under uncertainty is lacking. Here, we develop a mathematical and computational framework for probabilistic artificial intelligence (AI)–based predictive modeling combining data, expert knowledge, multiscale models, and information theory through uncertainty quantification and probabilistic graphical models (PGMs). We apply PGMs to chemistry specifically and develop predictive guarantees for PGMs generally. Our proposed framework, combining AI and uncertainty quantification, provides explainable results leading to correctable and, eventually, trustworthy models. The proposed framework is demonstrated on a microkinetic model of the oxygen reduction reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Editors’ Choice—Natural Convection Boundary Layer Thickness at Elevated Chloride Concentrations and Temperatures and the Effects on a Galvanic Couple

The natural convection boundary layer ( δ n c ) and its influence on cathodic current in a galvanic couple under varying electrolytes as a function of concentration (1 − 5.3 M NaCl) and temperature (25 °C−45 °C) were understood. Polarization scans were obtained under quiescent conditions and at defined boundary layer thicknesses using a rotating disk electrode on platinum and stainless steel 304L (SS304L); these were combined to determine δ n c . With increasing chloride concentration and temperature, δ n c decreased. Increased mass transport (Sherwood number) results in a decrease in δ n c , providing a means to predict this important boundary. Using Finite Element Modeling, the cathodic current was calculated for an aluminum alloy/SS304L galvanic couple as a function of water layer ( WL ) thickness and cathode length. Electrolyte domains were delineated, describing (i) dominance of ohmic resistance over mass transport under thin WL , (ii) the transition from thin film to bulk conditions at δ n c , and (iii) dominance of mass transport under thick WL . With increasing chloride concentration, cathodic current decreased due to decreases in mass transport. With increasing temperature, increased cathodic current was related to increases in mass transport and solution conductivity. This study has implications for sample sizing and corrosion prediction under changing environments.

25 ENERGY STORAGE↗
DOE OSTI2017

Grand challenges in understanding the interplay of climate and land changes

Half of the Earth s land surface has been altered by human activities, creating various consequences on the climate and weather systems at local to global scales, which in turn affects a myriad of land surface processes and our adaptation behaviors. We here review the status and major knowledge gaps of studying the interactions of land and atmospheric changes and present eleven grand challenge areas for scientific research and adaptation communities in the coming decade: (1) collective and separate impacts of major land changes and the interactions with non-land-change factors such as atmospheric CO2 increase, (2) carbon and other biogeochemical cycles, (3) climatically relevant biospheric emissions such as aerosols, (4) water cycle, (5) agriculture, (6) urbanization, (7) gradual acclimation of plants, communities, and ecosystems to climate and environmental changes, (8) plant migration, (9) land use projections, (10) reduction of uncertainties in models and data, and finally (11) adaptation strategies. We conclude that we need to create and maintain a close cross-disciplinary coordination between measurements and process representation in models to analyze complex multi-facet interrelated perturbations and feedbacks between land and climate changes. Along with major scientific research thrusts, land-use and land cover change mitigation and adaptation assessments should be strengthened to identify barriers that need to be overcome, evaluate and prioritize opportunities, and examine how decision making processes work in specific contexts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Optimizing the Costs of Solid Sorbent-Based CO 2 Capture Process Through Heat Integration

The focus of this project was the ADAsorb™ CO 2 Capture Process, a temperature-swing adsorption process that incorporates a three-stage fluidized bed as the adsorber and a single-stage fluidized bed as the regenerator. ADAsorb™ system was designed, fabricated, and tested under DOE award DEFE0004343. Two amine-based sorbents were evaluated in conjunction with the ADAsorb™ process: “BN”, an ion-exchange resin; and “OJ”, a metal organic framework (MOF) sorbent. Two cross heat exchanger designs were evaluated for use between the adsorber and regenerator: moving bed and fluidized bed. The fluidized bed approach was rejected fairly early in the project because the additional electrical load to power blowers or fans to overcome the pressure drop required for fluidization was estimated to be nominally three times the electrical power that could be generated from the steam saved through the use of the cross heat exchanger. The Energy Research Center at Lehigh University built and utilized a process model of the ADAsorb™ capture process and integrated this model into an existing model of a supercritical PC power plant. The Lehigh models verified that, for the ADAsorb™ system, the largest contributor to parasitic power was lost electrical generation, which was primarily electric power which the host plant could not generate due to the extraction of low pressure (LP) steam for sorbent heating, followed by power for the CO 2 compressor and the blower or fan power required to fluidize the adsorber and regenerator. Sorbent characteristics such as the impacts of moisture uptake, optimized adsorption and regeneration temperature, and sensitivity to changes in pressure were also included in the modeling study. Results indicate that sorbents which adsorb more than 1-2% moisture by weight are unlikely to be cost competitive unless they have an extremely high CO 2 working capacity that well exceeds 15% by weight. Modeling also revealed that reductions in adsorber pressure drop could negatively affect the CO 2 adsorption characteristics for sorbents with certain isobar adsorption characteristics like sorbent BN. Thus, reductions in pressure drop do not provide the efficiency benefits expected. A techno-economic assessment conducted during the project revealed that without heat integration, the a metal organic framework (MOF) sorbent used in conjunction with the ADAsorb™ process provided the opportunity for improved performance over the benchmark MEA process. While the addition of a cross heat exchanger and heat integration was found to significantly improve net unit heat rate, the additional equipment costs required to realize these improvements almost always outweighed the improvement in performance. The exception to this was for a supported amine sorbent and the addition of a moving bed cross heat exchanger alone or in conjunction with waste heat from the compressor used for supplemental regenerator heating. Perhaps one of the most important points to be drawn from the work conducted during this project is the significant influence of sorbent characteristics alone on the projected COE and LCOE associated with the ADAsorb™ process, and the implications associated with future improvements to solid sorbent CO 2 capture. The results from this project suggest that solid sorbent CO 2 capture will continue to see performance gains and lower system costs as further sorbent improvements are realized.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2016

Predictive Engineering Tools for Injection-Molded Long-Carbon-Fiber Thermoplastic Composites. Topical Report

This project aimed to integrate, optimize, and validate the fiber orientation and length distribution models previously developed and implemented in the Autodesk® Simulation Moldflow® Insight (ASMI) software package for injection-molded long-carbon-fiber (LCF) thermoplastic composite structures. The project was organized into two phases. Phase 1 demonstrated the ability of the advanced ASMI package to predict fiber orientation and length distributions in LCF/polypropylene (PP) and LCF/polyamide-6, 6 (PA66) plaques within 15% of experimental results. Phase 2 validated the advanced ASMI package by predicting fiber orientation and length distributions within 15% of experimental results for a complex three-dimensional (3D) Toyota automotive part injection-molded from LCF/PP and LCF/PA66 materials. Work under Phase 2 also included estimate of weight savings and cost impacts for a vehicle system using ASMI and structural analyses of the complex part. The present report summarizes the completion of Phases 1 and 2 work activities and accomplishments achieved by the team comprising Pacific Northwest National Laboratory (PNNL); Purdue University (Purdue); Virginia Polytechnic Institute and State University (Virginia Tech); Autodesk, Inc. (Autodesk); PlastiComp, Inc. (PlastiComp); Toyota Research Institute North America (Toyota); Magna Exteriors and Interiors Corp. (Magna); and University of Illinois. Figure 1 illustrates the technical approach adopted in this project that progressed from compounding LCF/PP and LCF/PA66 materials, to process model improvement and implementation, to molding and modeling LCF/PP and LCF/PA66 plaques. The lessons learned from the plaque study and the successful validation of improved process models for fiber orientation and length distributions for these plaques enabled the project to go to Phase 2 to mold, model, and optimize the 3D complex part.

42 ENGINEERING↗
DOE OSTI2017

Small Scale Field Test Demonstrating CO 2 Sequestration In Arbuckle Saline Aquifer And By CO 2 -Eor At Wellington Field, Sumner County, Kansas

The objectives of this project are to understand the processes that occur when a maximum of 70,000 metric tonnes of CO 2 are injected into two different formations to evaluate the response in different lithofacies and depositional environments. The evaluation will be accomplished through the use of both in situ and indirect MVA (monitoring, verification, and accounting) technologies. The project will optimize for carbon storage accounting for 99% of the CO 2 using lab and field testing and comprehensive characterization and modeling techniques. Site characterization and CO 2 injection should demonstrate state-of-the-art MVA tools and techniques to monitor and visualize the injected CO 2 plume and to refine geomodels developed using nearly continuous core, exhaustive wireline logs, and well tests and a multi-component 3-D seismic survey. Reservoir simulation studies will map the injected CO 2 plume and estimate tonnage of CO 2 stored in solution, as residual gas, and by mineralization and integrate MVA results and reservoir models shall be used to evaluate CO 2 leakage. A rapid-response mitigation plan was developed to minimize CO 2 leakage and provide a comprehensive risk management strategy. The CO 2 was intended to be supplied from a reliable facility and have an adequate delivery and quality of CO 2 . However, several unforeseen circumstances complicated this plan: (1) the initially negotiated CO 2 supply facility went offline and contracts associated with CO 2 supply had to be renegotiated, (2) a UIC Class VI permit proved to be difficult to obtain due to the experimental nature of the project. Both subjects are detailed in separate deliverables attached to this report. The CO 2 enhanced oil recovery (EOR) and geologic storage in Mississippian carbonate reservoir was sucessully deployed. Approximately 20,000 metric tons of CO 2 was injected in the upper part of the Mississippian reservoir to verify CO 2 EOR viability in carbonate reservoirs and evaluate a potential of transitioning to geologic CO 2 storage through EOR. A total of 1,101 truckloads, 19,803 metric tons—an average of 120 tonnes per day—were delivered over the course of injection that lasted from January 9 to June 21, 2016. After cessation of CO 2 injection, the KGS 2-32 well was converted to water injector and continues to operate. CO 2 EOR progression in the field was monitored weekly with fluid level, temperature, and production recording and formation fluid composition sampling. It is important to note that normally, CO 2 EOR pilots are less efficient than commercial operations due to lack of directional and precise well control, lack of surface facilities for CO 2 recycling, and other factors. As a result of this pilot CO 2 injection, the observed incremental average oil production increase was ~68% with only ~18% of injected CO 2 produced back. Decline curve analysis forecasts of additional cumulative oil produced were 32.44M STB to the end of 2027. Wellington Mississippian pilot efficiency by the end of forecast calculations is 11 MCF per barrel of produced oil. Using 32M STB oil production and $1,964,063 cost of CO 2 , CO 2 EOR cost per barrel of oil production is ~$60. Simple but robust monitoring technologies proved to be very efficient in detecting and locating CO 2 . High CO 2 reservoir retentions with low yields within an actively producing field could help to estimate real-world risks of CO 2 geological storage for future projects. The Wellington Field CO 2 EOR was executed in a controlled environment with high efficiency. This case study proves that CO 2 EOR could be successfully applied in Kansas carbonate reservoirs if CO 2 sources and associated infrastructure are available. Recent developments in unconventional resources development in Mid-Continent USA and associated large volume disposal of backflow water and the resulting seismic activity have brought more focus and attention to the Arbuckle Group in southern Kansas. Despite the commercial interest, limited essential information about reservoir properties and structural elements has impeded the management and regulation of disposal, an issue brought to the forefront by recent seismicity in and near areas of large volumes and rates of brine disposal. The Kansas Geological Survey (KGS) collected, compiled, and analyzed available data, including well logs, core data, step rate tests, drill stem tests, 2-D and 3-D seismic data, water level measurements, and others types of data. Several exploratory wells were drilled and core was collected and modern suites of logs were analyzed. Reservoir properties were populated into several site-specific geological models. The geological models illustrate the highly heterogeneous nature of the Arbuckle Group. Vertical and horizontal variability results in several distinct hydro-stratigraphic units that are the result of both depositional and diagenetic processes. During the course of this project, it has been demonstrated that advanced seismic interpretation methods can be used successfully for characterization of the Mississippian reservoir and Arbuckle saline aquifer. Analysis of post-stack 3-D seismic data at the Mississippian reservoir showed the response of a gradational velocity transition. Pre-stack gather analysis showed that porosity zones of the Mississippian and Arbuckle reservoirs exhibit characteristic amplitude versus offset (AVO) response. Simultaneous AVO inversion estimated P- and S-impedances. The 3-D survey gather azimuthal anisotropy analysis (AVAZ) provided information about the fault and fracture network and showed good agreement to the regional stress field and well data. Mississippian reservoir porosity and fracture predictions agreed well with the observed mobility of injected CO 2 in KGS well 2-32. Fluid substitution modeling predicted acoustic impedance reduction in the Mississippian carbonate reservoir introduced by the presence of CO 2 . Seismicity in the United States midcontinent has increased by orders of magnitude over the past decade. Spatiotemporal correlations of seismicity to wastewater injection operations have suggested that injection-related pore fluid pressure increases are inducing the earthquakes. In this investigation, we examine earthquake occurrence in southern Kansas and northern Oklahoma and its relation to the change in pore pressure. The main source of data comes from the Wellington Array in the Wellington oil field, in Sumner County, Kansas, which has monitored for earthquakes in central Sumner County, Kansas, since early 2015. The seismometer array was established to monitor CO 2 injection operations at Wellington Field. Although no seismicity was detected in association with the spring 2016 Mississippian CO 2 injection, the array has recorded more than 2,500 earthquakes in the region and is providing valuable understanding to induced seismicity. A catalog of earthquakes was built from this data and was analyzed for spatial and temporal changes, stress information, and anisotropy information. The region of seismic concern has been shown to be expanding through use of the Wellington earthquake catalog, which has revealed a northward progression of earthquake activity reaching the metropolitan area of Wichita. The stress orientation was also calculated from this earthquake catalog through focal mechanism inversion. The calculated stress orientation was confirmed through comparison to other stress measurements from well data and previous earthquake studies in the region. With this knowledge of the stress orientation, the anisotropy in the basement could be understood. This allowed for the anisotropy measurements to be correlated to pore pressure increases. The increase in pore pressure was monitored through time-lapse shear-wave anisotropy analysis. Since the onset of the observation period in 2010, the orientation of the fast shear wave has rotated 90°, indicating a change associated with critical pore pressure build up. The time delay between fast and slow shear wave arrivals has increased, indicating a corresponding increase in anisotropy induced by pore pressure rise. In-situ near-basement fluid pressure measurements corroborate the continuous pore pressure increase revealed by the shear-wave anisotropy analysis over the earthquake monitoring period. This research is the first to identify a change in pore fluid pressure in the basement using seismological data and it was recently published in the AAAS journal Science Advances (Nolte et al., 2017). The shear-wave splitting analysis is a novel application of the technique, which can be used in other regions to identify an increase in pore pressure. This increasing pore fluid pressure has become more regionally extensive as earthquakes are occurring in southern Kansas, where they previously were absent. These monitoring techniques and analyses provide new insight into mitigating induced seismicity’s impact on society.

02 PETROLEUM↗