Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “rock failure”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Impact of Electrolyte Additives on the Lifetime of High Voltage NMC Lithium-Ion Pouch Cells

This work involves improving the lifetime of lithium-ion cells during high voltage cycling using electrolyte additives. Three generations of electrolyte additives were investigated and screened in NMC442/graphite pouch cells using a 24 h voltage-hold protocol at 40 degrees C to accelerate oxidative reactions occurring at 4.4 V. Once promising additives and combinations were identified, they were then tested in cobalt-free NMC640/graphite cells for long-term cycling to upper cutoff voltages of 4.3, 4.4, and 4.5 V at temperatures of 20, 40, and 55 degrees C. Degradation mechanisms were probed using dV/dQ analysis, micro-X-ray fluorescence spectroscopy, and electrochemical impedance spectroscopy. The primary failure mode of cells held at high voltages is due to increase in cell impedance, which is correlated to the dissolution of transition metals, specifically manganese, originating from the positive electrode. We believe this dissolution is presumably due to the formation of a high impedance rock salt surface layer on the NMC positive electrode particles. Such deleterious outcomes can be limited by selecting an appropriate electrolyte additive package. It is hoped that this paper can provide a starting point for developing NMC Li-ion cells that can operate to voltages as high as 4.4 V and still display long lifetimes.

Azam, Saad↗
DOE OSTI2025

Use of Rig Parameter Data in Bit Constraint Models for Improved Drilling Performance at The Geysers

Surface parameter measurements are routinely used during deep well construction to monitor and guide drilling conditions for improved performance and reduced costs. However, these measurements are of reduced value without a standard to aid in evaluation and decision making. A method is demonstrated whereby drill bit constraint models are used to interpret drilling response parameters. Drill rig parameter data for well GDC-36 at the Geysers Geothermal Field Power were acquired by Geysers Power Company and drilling contractor Kenai Drilling using Pason US DataHub and evaluated. Drilling parameters are evaluated using laboratory-validated rock reduction models for predicting the phenomenological response of drag bits (Detournay and Defourny, 1992) along with other model constraints in computational algorithms. The method is used to evaluate overall bit performance, monitor bit integrity, and detect the presence of drillstring vibrations and other conditions contributing to bit failure; comparisons are made to observations of bit wear and damage. The method will be applied in real-time to improve decision-making on subsequent wells and has applicability to development of advanced analytics on future geothermal wells using real-time electronic drilling recorder (EDR) data for improved performance and reduced drilling costs.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Development of Leaching Processes for the Recovery of Rare Earth Elements from Acid Mine Drainage Precipitates

Acid mine drainage (AMD) has been a challenge for mine operators to address and has had significant environmental impacts when there is a failure to address it. Fortunately, there is a regulation in the US that requires the treatment of AMD before water can be discharged into the environment. AMD is generated when sulfide minerals are exposed to water and air from mining. The acidic water then leaches metals from the surrounding rock, creating the potential for environmental contamination of this acidic water containing dissolved metals. AMD can contain high concentrations of metals like iron, aluminum, and manganese and also have been shown to contain trace concentrations of critical minerals, including rare earth elements (REEs). This environmental waste stream is now being researched as a potential source for REEs. There is a patented process for processing and extracting REEs from AMD which produces rare earth oxide preconcentrate (REOP) from AMD treatment precipitates and a final stage mixed rare earth oxide (MREO) product. This body of research examines a selective leaching process for each of these products and determines the activation energy associated with the leaching processes. Acid leaching with HCl was examined to selectively leach REEs from REOP while contaminant metals remained in the solid residue. The maximum leaching recovery of the REEs from the REOP was approximately 90% and an activation energy of 6.4 kJ/mol at pH 3.0. Ammonium chloride leaching was examined to selectively remove contaminant metals from a MREO product. The ammonium chloride process successfully leached major contaminant metals in excess of 85% recovery and had a maximum activation energy of 40.0 kJ/mol.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Analyzing the Validity of Brazilian Testing Using Digital Image Correlation and Numerical Simulation Techniques

Characterizing the mechanical behavior of rocks plays a crucial role to optimize the fracturing process in unconventional reservoirs. However, due to the intrinsic anisotropy and heterogeneity in unconventional resources, fracture process prediction remains the most significant challenge for sustainable and economic hydrocarbon production. During the deformation tracking under compression, deploying conventional methods (strain gauge, extensometer, etc.) is insufficient to measure the deformation since the physical attachment of the device is restricted to the size of the sample, monitoring limited point-wise deformation, producing difficulties in data retrieval, and a tendency to lose track in failure points, etc. Where conventional methods are limited, the application of digital image correlation (DIC) provides detailed and additional information of strain evolution and fracture patterns under loading. DIC is an image-based optical method that records an object with a camera and monitors the random contrast speckle pattern painted on the facing surface of the specimen. To overcome the existing limitations, this paper presents numerical modeling of Brazilian disc tests under quasi-static conditions to understand the full-field deformation behaviors and finally, it is validated by DIC. As the direct tensile test has limitations in sample preparation and test execution, the Brazilian testing principle is commonly used to evaluate indirectly the tensile strength of rocks. The two-dimensional numerical model was built to predict the stress distribution and full-field deformation on Brazilian disc under compression based on the assumptions of a homogenous, isotropic and linear elastic material. The uniaxial compression test was conducted using the DIC technique to determine the elastic properties of Spider Berea sandstone, which were used as inputs for the simulation model. The model was verified by the analytical solution and compared with the digital image correlation. The numerical simulation results showed that the solutions matched reasonably with the analytical solutions where the maximum deviation of stress distribution was obtained as 14.59%. The strain evolution (normal and shear strains) and displacements along the central horizontal and vertical planes were investigated in three distinguishable percentages of peak loads (20%, 40%, and 90%) to understand the deformation behaviors in rock. The simulation results demonstrated that the strain evolution contours consistently matched with DIC generated contours with a reasonable agreement. The changes in displacement along the central horizontal and vertical planes showed that numerical simulation and DIC generated experimental results were repeatable and matched closely. In terms of validation, Brazilian testing to measure the indirect tensile strength of rocks is still an issue of debate. The numerical model of fracture propagation supported by digital image correlation from this study can be used to explain the fracturing process in the homogeneous material and can be extended to non-homogeneous cases by incorporating heterogeneity, which is essential for rock mechanics field applications.

58 GEOSCIENCES↗
DOE OSTI2023

Discrete Measurements of the Least Horizontal Principal Stress from Core Data: An Application of Viscoelastic Stress Relaxation

Summary The emerging Paradox Oil Play in southeastern Utah is among the most significant unconventional plays in the western USA. The mean total undiscovered oil resources within just the Pennsylvanian Cane Creek interval of the Paradox Basin are believed to exceed 215 million barrels. However, to date, less than 5% (~9 million barrels) of the total Cane Creek resource has been produced from fewer than 40 wells, and only approximately one-half of those are horizontal wells. More than 95% of production is from the central Cane Creek Unit (CCU). Natural fractures are a key feature of many production wells, but stimulation by induced hydraulic fractures is not consistently successful. We hypothesize that more effective production in this play will rely on better fundamental characterization, especially on better quantification of the state of stress. Approximately 110 ft of core, well logs, and a diagnostic fracture injection test (DFIT) were acquired from the State 16-2 well within the CCU. With these data, we applied two methods to constrain and clarify the state of stress. The first technique, the Simpson’s coefficient method, provides lower bounds on the two horizontal principal stresses and relies on only limited data. Alternatively, the viscoelastic stress relaxation (VSR) method is used to estimate the least horizontal principal stress, building on observations that principal stresses become more isotropic as the viscous behavior of a rock is more pronounced. Results of these two methods support the hypothesis that the state of stress in the CCU of the Paradox Basin is nearly lithostatic and isotropic. Other factors consistent with this hypothesis include high formation pore pressure, which tends to reduce the possible stress states by changing the frictional failure equilibrium; lack of induced fractures in the core, which should be present in the case of stress anisotropy; and interbedded halite layers, which given their high degree of ductility, probably lead to greater VSR for the entire sedimentary package.

Energy & Fuels↗
DOE OSTI2021

Fixed Bit Rotary Drilling Failure Criteria Effect on Drilling Vibration

Polycrystalline diamond compact (PDC) bit use has been expanded to geothermal drilling applications due to their improved performance in drilling compared to roller cone bits. Multiple vibration models have been developed yet and there are limited use of field data to study the ability of numerical vibration model to predict recorded vibration in the field. The objective of this paper is to evaluate the capabilities of two different rock-bit interaction models with two different ROP formulations for predicting the actual behavior seen in the field. The first model (Model-I) considers the rock-bit interaction, incorporating the effect of subsurface formation characteristics. The second model (Model-II) is a PDC bit model that uses calculated threshold weight-on-cutter (WOCt) to study the interaction of the bit with the rock surface. The models, which couple the drill stem axial and torsional motions, are solved numerically using the fourth-order Runge-Kutta method and compared with recorded downhole drilling data from the Chocolate Mountain Aerial Gunnery Range (CMAGR) geothermal well located in California. An 8–½" PDC bit with 8 blades was modeled and vibration simulation was performed to investigate stick-slip vibrations at 1400 ft. depth. Simulation results using field data indicate that Model-II encounters a low range of stick-slip severity (SSV) as compared to Model-I, which is in the moderate range. Comparison with the field value for SSV shows Model-II is the best representative of the field case for this depth. The sensitivity analysis shows that results are sensitive to model constants making prediction of SSV from models uncertain.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Modeling Validation of Mechanical and Thermal Borehole Breakouts From Polyaxial Laboratory Tests to Determine Maximum Horizontal in Situ Stress

To accurately assess rock behavior in the deep subsurface, it is necessary to measure the in situ stress directions and magnitudes. The current methods for measuring the in situ stress state in the deep subsurface primarily include hydraulic fracturing tests (i.e., the minimum horizontal stress) or occasionally observing existing compressive borehole breakouts (i.e., the maximum horizontal stress) that have occurred naturally from drilling. If there are no existing compressive breakouts in a borehole, the maximum horizontal in situ stress cannot be estimated with much confidence. In response to this data gap, a new thermal breakout technology is being developed that will provide a method for thermally inducing borehole breakouts and obtaining consistent measurements of the maximum horizontal stress magnitude. This thermal breakout technology involves heating the borehole and increasing the thermoelastic compressive stress in the rock until a breakout develops, which is directly correlated to the maximum horizontal stress magnitude. In support of developing the thermal breakout technology, polyaxial laboratory tests have been performed on small-scale boreholes within rock blocks where mechanically- and thermally-induced borehole breakouts have been created. Numerical models along with the principal of superposition were created and used to analyze the polyaxial laboratory tests to predict the maximum horizontal stress, given the same data that would be obtained in an actual subsurface borehole thermal breakout test. Multiple failure criteria were used to evaluate the best prediction of the breakout onset and maximum horizontal stress. The maximum horizontal stress predictions were compared to the actual maximum horizontal stress applied at breakout using acoustic events recorded from emission sensors in the polyaxial tests. The results showed consistent results that could be used to refine the modeling approach and failure criterion that are used to make the maximum horizontal stress predictions. This study provided insight and validation for the thermal breakout stress measurement concept.

Jones, Matt↗
DOE OSTI2021

Probing the Damage Zone at Parkfield

Rocks are heterogeneous materials that exhibit nonlinear elastic (anelastic) behavior at scales ranging from the laboratory to Earth. In the laboratory, typical, complex relationships exist between stress and strain that include hysteresis, finite relaxation times, strain rate, and history dependence. These behaviors are linked to important characteristics such as stress, porosity, permeability, material integrity, and material failure. In this work, we adopted a “pump-probe” type experiment common in laboratory studies, using solid earth tides as the low-frequency pump and empirical Green's function as the high-frequency probe. By probing the velocity at different points in the pump cycle, we constrained important information about the strain-modulus relationship. Near the San Andreas Fault, we observed strongly nonlinear elastic behavior that characterizes the damage zone. We also constrained important aspects of hysteretic behavior that are related to damage properties and possibly pore pressure. Away from the fault, the nonlinear behavior is diminished.

58 GEOSCIENCES↗
DOE OSTI2025

Experimental data for damage mechanics simulation challenge

While there are many computational approaches for simulating damage in rock and other materials, few have been ground truth tested with either known experimental data or with blind data sets. Here, in this work, we present a bench-mark laboratory data set for a damage mechanics challenge to compare computational approaches on damage evolution in brittle-ductile materials. The samples were fabricated through additive manufacturing to produce repeatable specimens designed to fail in controlled ways. The failure was induced in the samples using a 3-point bending test to produce different Modes such as Mode I and mixed Modes including I-II, I-III and I-II-III Modes to generate a calibration data set and a blind challenge data set. Data collected included spatial and temporal measurements from traditional digital load–displacement sensors, 2D digital image correlation measurement to map surface deformations, 3D X-ray microscopy to ground-truth the crack-failure geometry, and laser profilometry to capture surface roughness. The data sets are available, on a data repository, to the community to advance computational models to improve our ability to predict damage in brittle-ductile materials.

3-point bending↗
DOE OSTI2020

Macro- and micro-compression testing of shales

Here, macro- and micro-geomechanical properties of shales were investigated at millimeter- and micrometer-scale using an example of the Woodford Shale. Geomechanical properties, such as uniaxial (unconfined) compressive strength and Young's modulus, were quantified at the millimeter-scale and compared with results from micro-compression testing experiments of FIB-SEM-nanofabricated micro-pillars. Size-scale and compositional/structural heterogeneity effect on elasto-plastic deformation and failure behavior of shales were investigated. Also, relationship between elemental/mineral composition and unconfined compressive strength and Young's modulus, at a micrometer-scale, was discussed. It was shown that geomechanical properties of shales are scale-dependent and are strongly affected by compositional/structural anisotropy. In particular, micro-compression testing experiments showed non-uniformly distributed compressive strength and Young's modulus within the investigated Woodford Shale rock samples. It was demonstrated, that geomechanical properties of shales, investigated at the micrometer-scale, tend to be affected by different minerals and/or pores, while the same properties, investigated at the millimeter-scale, are governed by different micro-facies (and/or micro-fractures) present within the rock. This shows that geomechanical properties of shales can be dramatically different depending on the scale of investigation and compositional/structural heterogeneity of these rocks, and therefore are not easily transferable across the scales.

58 GEOSCIENCES↗
DOE OSTI2022

Precursory Off-Fault Deformation in Restraining and Releasing Step Overs: Insights From Discrete Element Method Models

Accelerating geophysical activity is detected preceding some, but not all, large earthquakes. This observation may indicate that no precursors occur before some earthquakes, or that the instrumentation lacks the required sensitivity. Here, to aid crustal monitoring efforts, we use discrete element method models to identify the locations and styles of deformation that may provide useful information about approaching fault reactivation. We model the reactivation of two healed rough faults in a variety of step over configurations, embedded in a host rock with varying amounts of damage subject to shear velocity loading parallel to the faults. Both the fault geometry and ratio of fault to host rock strength control the amount of off-fault deformation. Consistent with field observations, models with larger steps and more preexisting host rock damage produce higher amounts of off-fault deformation. We assess the size of the continuous regions of high velocities and strains to compare the value of the precursory information of each velocity and strain component. Comparing the three components of the velocity vector suggests that the fault-parallel velocity produces the largest and most temporally continuous regions of elevated velocity. The size of these regions increases toward failure, indicating the usefulness of tracking this component. Comparing the volumetric and shear components of the three-dimensional strain tensor suggests that during most of the interseismic period, the shear strain provides more information about approaching fault slip than the volumetric strain. However, in the days and months preceding fault reactivation, both the shear and volumetric strains provide similarly valuable information.

58 GEOSCIENCES↗
DOE OSTI2026

Lattice-Oxygen-Driven Selective Oxidation Strategy for Stable Argyrodite Solid-State Lithium Metal Batteries

All-solid-state lithium metal batteries (ASSLMBs) with Li6PS5Cl argyrodite electrolytes and high-voltage LiNi0.8Mn0.1Co0.1O2 (NMC811) cathodes offer high energy density but suffer from rapid capacity fading due to the layered-to-rock-salt transition of NMC811 and structural degradation of Li6PS5Cl from parasitic interfacial reactions. Here, we demonstrate a catholyte engineering strategy using a Li2S scavenging additive to suppress interfacial reactivity and preserve the structural and electrochemical stability of both NMC811 and Li6PS5Cl. Incorporating 0.10 wt.% Li2S enables exceptional cycling stability, achieving 76% capacity retention after 550 cycles at C/10 and 88% retention after 800 cycles at C/3 at 60 degrees C, compared with rapid failure in pristine cells. Spectroscopic, electrochemical, and morphological analyses confirm that Li2S maintains electrode integrity by sustaining particle contact and suppressing phase decomposition. This work elucidates interfacial degradation pathways in NMC811/argyrodite systems and introduces a low-cost, scalable strategy to stabilize nickel-rich oxide cathodes in ASSLMBs, advancing their practical viability.

25 ENERGY STORAGE↗
DOE OSTI2020

Thermally-Induced Breakouts: Insights from True-Triaxial Tests with Acoustic Emission Monitoring

Borehole breakouts are used to constrain the magnitude of maximum horizontal stress. However, when the borehole wall strength is higher than the in situ tangential stress, borehole wall failure does not develop. Additional compressive stress can be induced by heating borehole walls. To validate this concept experimentally, we conducted room-temperature and elevated temperature true-triaxial tests on Berea sandstone and Niagaran dolomite samples. We used acoustic emission sensors to capture the onset of breakout development, and we measured the temperature close to borehole wall to assess the magnitude of induced thermal hoop stress. The test results show that within a specific rock type, the breakouts develop in similar manner in room-temperature and elevated-temperature tests. Therefore, the maximum horizontal stress can be constrained from the following dataset: critical tangential stress at which breakout develops, minimum horizontal stress, elastic and thermal properties, and temperature change at the borehole wall.

Trzeciak, Maciej↗
DOE OSTI2021

Study on Foam Fracturing of Granite for the Development of Enhanced Geothermal Systems

Foam fracturing has been considered as a potential technical approach to address water concerns with hydraulic stimulation in the development of enhanced geothermal systems (EGS). This paper reports the recent progress of a project sponsored by the U.S. DOE GTO Waterless Stimulation Initiative. The project investigates the foam fracturing as an alternative approach to the hydraulic fracturing. The paper describes the development of foam fracturing testing system at Oak Ridge National Laboratory (ORNL). The system can be used to perform the hydraulic fracturing of rock-like materials with both water and foamed fluids under pressure to 6,000 psi (41.4 MPa). The system features the sections of foam generation and foam injection and has the capability of testing materials in both monotonic and cyclic injection modes.Experimental results of foam fracturing on Charcoal Black granite will be reported. Cylindrical granite specimens with a blind hole were studied along with water and nitrogen-gas-in-water foam as the fracturing fluids. The effects of fracturing fluids and injection modes on the breakdown pressure and failure responses of the material were examined. The observations from the experimental work and the implications to the EGS application will be presented and discussed.

Wang, Hong↗
DOE OSTI2020

Estimating perturbed stress from 3-D borehole displacements induced by fluid injection in fractured or faulted shales

SUMMARY Hydrofracturing stress measurements in fractured and anisotropic shales are notoriously difficult, because opening of existing geological features tends to prevent the creation of a pure hydraulic fracture perpendicular to the least compressive principal stress. Here we show how adding 3-D borehole-displacement measurements while conducting the hydraulic injection test helps to better constrain the principal stress orientations and magnitudes. We developed a 3-D fully coupled hydromechanical numerical model to analyse the displacement, fluid pressure and injection flow-rate data measured during an injection pressure-step-rate test conducted to activate a faulted borehole interval in the Mont Terri Opalinus Clay (Switzerland). We find that injected fluids can only penetrate the fault when it is at or above the Coulomb failure pressure. Borehole displacement orientations are sensitive to a ∼15° variation in the stress–tensor orientation and a 1 MPa stress magnitude variation. Although some dispersion occurs while rupture is propagating along the fault plane ∼4 m away from the borehole, the maximum density of displacement orientations consistently informs about the stress orientation. Thus, an extended injection step-rate approach coupled with an accurate in situ measurement of the borehole wall displacements can be used to better constrain the local stress field perturbations in fractured shales and in heterogeneous rock in general.

Guglielmi, Yves↗
DOE OSTI2020

Experimental Evaluation of Geopolymer, Nano-Modified, and Neat Class H Cement by Using Diametrically Compressive Tests

Abstract Cement is the prime barrier material to maintain zonal isolation during wellbore life-cycle. Different types of the cement have different performances on sealability. The primary objective of this study is evaluation of three types of cement systems, neat Class H, nano-synthetic graphite, and rock-based geopolymer cements, on interfacial bond integrity. Diametrical compressive tests are used on casing–cement–formation systems. Cement performances are tested in a conventional configuration (Berea Sandstone). A universal testing system is utilized to analyze the relationship applied load–axial strain and the digital image correlation (DIC) technique is used to generate the strain field maps. Based on the relationship and DIC results, bond strengths in different cement systems are evaluated. Additionally, microscope images provide insight explanations of the bond from micro-view. Based on the experimental results, failures of all cases initiate and propagate from the casing–cement interface. Then, radial cracks develop through the formation. This is because the mechanical interlocking bonds are formed only at the cement–sandstone interface to enhance the bond strength. The comparison of geopolymer and nano-modified cements with the conventional cement (neat Class H) shows that the unconventional cements have a lower maximum load and higher axial strain before failure of the samples. It indicates that both the two kinds of cements are more ductile than conventional cement and can support more deformation before the initiation of cracks.

Energy & Fuels↗
DOE OSTI2022

Operation Upshot-Knothole

A thundering boom and sheet of fire announced a dramatic, once in a lifetime event – the firing of a nuclear projectile from a 280 mm cannon. Travelling seven miles downrange, the projectile detonated 500 feet over the Nevada Test Site’s Frenchman Flat with a force of fifteen kilotons. Although the test was successful, the massive cannon, designed for tactical battlefield use, had a critical flaw. It was too large and cumbersome for mobile combat. The cannon test, codenamed Grable, was part of Operation Upshot-Knothole, the fourth test series conducted in Nevada. Beginning with the Annie test on March 17, 1953, and concluding with the Climax test on June 4th, the Operation consisted of three airdrops, seven tower detonations, and the one artillery shot. As with previous Nevada tests, the demarcation between the weapon development series, Upshot, and the DOD effects tests, Knothole, was blurred. All eleven tests generated data in support of weapon development as well as tactical battlefield applications. During Upshot-Knothole, the University of California Radiation Laboratory (UCRL), conducted its first ever tests, Ruth and Ray. Although neither test was successful, these failures, along with another a year later at Operation Castle, provided the new laboratory with design data leading to future successful tests. The Federal Civil Defense Administration (FCDA) participated in a test series with the general goal of preparing for societal recovery from a nuclear attack. Known as Operation Doorstep, the FCDA’s experiments evaluated the effects of a nuclear attack on food supplies, housing and urban structures, and even a copse of trees transplanted from nearby Mt. Charleston. A contingent of the press witnessed test Annie from an outcropping of rocks, subsequently named News Nob. The presence of the news media, along with the FCDA’s participation, marked a fundamental change in the Nevada test program - a change that that now included the general populace, whose survival was, of course, the rationale for testing. The Department of Defense continued and expanded its Desert Rock exercises of troop exercises, helicopter operations, and evaluation of blast damage to military aircraft, tanks, and guns. Annie, Nancy, Badger, Simon, Encore, and Grable involved larger numbers of personnel than the other five shots. Below is an edited and condensed account of Desert Rock activities compiled by the Defense Threat Reduction Agency.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Real-Time Drilling Optimization System for Improved Overall Rate of Penetration and Reduced Cost Per Foot in Geothermal Drilling

The key to success in geothermal drilling is economic feasibility, and a major cost in the development of geothermal resources is the actual drilling of the wells. In this project, a real-time drilling optimization system for geothermal drilling was developed. The system couples three individual components while drilling. The first component is a drill stem vibration analysis model, the second is Mechanical Specific Energy (MSE) analyses, and the third is a detailed PDC Rate of Penetration (ROP) drill bit model for optimum RPM and WOB combinations. The benefit of the coupled system is that the range of WOB and RPM could be selected to avoid drill stem vibrations. Secondly, MSE is used as an efficiency measure and the detailed PDC drill bit model ensures the drill bit does not endure temperatures that exceed the temperature at which the PDC cutters experience accelerated wear. The new detailed PDC bit model is based on rock/bit interaction that physically tracks the PDC cutter wear flats as the bit drills ahead giving the capability to calculate the temperature being generated underneath the worn cutters to better advise on operational parameters to avoid accelerated cutter wear and failure and to ensure that operational parameters are applied so that overall ROP is maximized. By combining the drill stem vibrations and the detailed PDC bit cutter wear and “safe” non-accelerated cutter wear temperature and optimum ranges of operating parameters, it results in higher ROP and lower cost drilling. Single cutter PDC testing performed in different lithologies at Sandia was utilized to verify the PDC cutter forces and depth of cut for new and worn cutters. Based on single cutter PDC temperature modeling, verification using single cutter data from the testing done by National Oilwell Varco (NOV) was performed. Sandia’s Hard-Rock Drilling Facility (HRDF) was utilized to test different drill bit configurations with different cutter designs and wear status with different induced modes of vibration to obtain the critical bit RPM/WOB ranges resulting in ineffective drilling and low ROP. The collected test data were further used to verify and calibrate the full hole PDC ROP model that was developed based on single cutter interaction data. A full coupled drill stem vibration model was formulated and verified with geothermal field data from the Chocolate Mountain Aerial Gunnery Range (CMAGR). A graphical user interface (GUI) was developed using Tkinter library in the computer programming language Python, which integrates all the developed models in one system. The developed system consists mainly of the PDC ROP model, PDC bit wear model, PDC cutter temperature model, Mechanical Specific Energy (MSE) model, and drillstring vibration model integrated into one system. The developed system can be used for both, post well analysis and real-time optimization using different criteria such as ROP maximization or MSE minimization. The software uses Differential Evolution Algorithm (DEA) to find optimum values for operational parameters based on last foot drilled while avoiding the drillstring vibration and cutter temperature critical operating parameters.

15 GEOTHERMAL ENERGY↗