The integrity of wellbores at the interbed between the caprock and salt is a serious concern in the Big Hill site. For the remediation and life extension of wellbores, more accurate predictions from the global model are needed. The Big Hill global model is improved using the M-D viscoplastic contact surface model and the mesh containing the interbed layer with contact surfaces between the salt and caprock layers, and fault blocks in overburden and caprock layers. The model calibration has been performed based on the cavern volumetric closures obtained from the Caveman calculations. The results agree well from 1991 to the early 2000s. The difference starts to widen after that, it might be because of frequent fluid movement and raw water injection. Therefore, the predictions from this improved model could be used to examine the structural integrity of caverns in Big Hill salt dome.
BC-4 is an abandoned brining cavern situated in the middle of the site. Its presence poses a concern for several reasons: 1) the cavern was leached up into the caprock; 2) it is similar to BC-7, a brining cavern on the northwest corner of the dome that collapsed in 1954 and now is the home to Cavern Lake; 3) a similar collapse of BC-4 would have catastrophic consequences for the future operation of the site. There exists a previously mapped fault feature in the caprock and thought to extend into the salt dome than runs in close proximity to BC-4. There are uncertainties about the true extent of the fault, and no explicit analysis has been performed to predict the effects of the fault on BC-4 stability. Additional knowledge of the fault and its effects is becoming more crucial as an enhanced monitoring program is developed and installed.
It has been recognized that as cavern operations become more frequent due to oil sales, field conditions may arise which require a faster turnaround time of analysis to address potential cavern impacts. This letter describes attempts to implement a strategy of transferring an intermediate solution of a Big Hill (BH) geomechanical model from a previous finite element mesh with a specified cavern geometry, to a new mesh with a new cavern geometry created by leaching from an oil sale operation.
The Department of Energy maintains an up-to-date documentation of the number of available full drawdowns of each of the caverns owned by the Strategic Petroleum Reserve (SPR). This information is important for assessing the SPR's ability to deliver oil to domestic oil companies expeditiously if national or world events dictate a rapid sale and deployment of the oil reserves. Sandia was directed to develop and implement a process to continuously assess and report the evolution of drawdown capacity, the subject of this report. A cavern has an available drawdown if after that drawdown, the long-term stability of the cavern, the cavern field, or the oil quality are not compromised. Thus, determining the number of available drawdowns requires the consideration of several factors regarding cavern and wellbore integrity and stability, including stress states caused by cavern geometry and operations, salt damage caused by dilatant and tensile stresses, the effect of enhanced creep on wellbore integrity, and the sympathetic stress effect of operations on neighboring caverns. A consensus has now been built regarding the assessment of drawdown capabilities and risks for the SPR caverns. The process involves an initial assessment of the pillar-to-diameter (P/D) ratio for each cavern with respect to neighboring caverns. A large pillar thickness between adjacent caverns should be strong enough to withstand the stresses induced by closure of the caverns due to salt creep. The first evaluation of P/D includes a calculation of the evolution of P/D after a number of full cavern drawdowns. The most common storage industry standard is to keep this value greater than 1.0, which should ensure a pillar thick enough to prevent loss of fluids to the surrounding rock mass. However, many of the SPR caverns currently have a P/D less than 1.0 or will likely have a low P/D after one or two full drawdowns. For these caverns, it is important to examine the structural integrity with more detail using geomechanical models. Finite-element geomechanical models have been used to determine the stress states in the pillars following successive drawdowns. By computing the tensile and dilatant stresses in the salt, areas of potential structural instability can be identified that may represent "red flags" for additional drawdowns. These analyses have found that many caverns will maintain structural integrity even when grown via drawdowns to dimensions resulting in a P/D of less than 1.0. The analyses have also confirmed that certain caverns should only be completely drawn down one time. As the SPR caverns are utilized and partial drawdowns are performed to remove oil from the caverns (e.g., for occasional oil sales authorized by the Congress or the President), the changes to the cavern caused by these procedures must be tracked and accounted for so that an ongoing assessment of the cavern's drawdown capacity may be continued. A proposed methodology for assessing and tracking the available drawdowns for each cavern was presented in Sobolik et al. (2018). This report includes an update to the baseline drawdowns for each cavern, and provides an initial assessment of the evolution of drawdown expenditure for several caverns
A finite element numerical analysis model , that consists of a realistic mesh capturing the geometries of Big Hill (BH) Strategic Petroleum Reserve (SPR) site using the multi - mechanism deformation (M - D) salt constitutive model and including data taken daily of the wellhead pressure and level of the oil - brine interface , has been upgraded . The upgraded model contains the shear zone to examine the interbed behavior in a realistic manner. The salt creep rate is not uniform in the salt dome, and creep test data for BH salt is limited. Therefore, a model calibration is necessary to simulate the geomechanical behavior of the salt dome. Cavern volumetric closures of SPR caverns calculated from sonar survey reports are used for the field baseline measurement. The structure factor, A 2 , and transient strain lim it factor, K 0 , in the M - D constitutive model are used for model calibration. An A 2 value obtained experimentally from the BH salt and K 0 value of WIPP salt are used as the baseline values. To adjust the magnitude of A 2 and K 0 , multiplication factors A2F an d K0F are defined, respectively. The A2F and K0F values of the salt dome and salt drawdown layer of elements surrounding each SPR cavern have been determined through a number of back fitting analyses. The trendlines of the predictions and sonar data match up well for BH 101, 103, 104, 106, 110, 111, 112, and 113. The prediction curves are close to the sonar data for BH 102 and 114. However, the prediction curves for BH 105, 107, 108, and 109 are not close to the sonar data. An inconsistency was found in the sonar data, i.e. the volume measured later is larger than that before in some time intervals , even if the leached volume is taken into account , for BH 101, 104, 106, 107, and 112. Project discussions are needed to determine possibilities on how to resolve the issues and determine the best path forward for future computer modeling attempts.
The Department of Energy maintains an up-to-date documentation of the number of available full drawdowns of each of the caverns owned by the Strategic Petroleum Reserve (SPR). This information is important for assessing the SPR's ability to deliver oil to domestic oil companies expeditiously if national or world events dictate a rapid sale and deployment of the oil reserves. What factors go into assessing available drawdowns? Determining the number of drawdowns requires the consideration of several factors regarding cavern and wellbore integrity and stability, including stress states caused by cavern geometry and operations, salt damage caused by dilatant and tensile stresses, the effect of enhanced creep on wellbore integrity, and the sympathetic stress effect of operations on neighboring caverns. A consensus has now been built regarding the assessment of drawdown capabilities and risks for the SPR caverns. The process involves an initial assessment of the pillar-to-diameter (P/D) ratio for each cavern with respect to neighboring caverns. Ideally, it is desired to keep this value greater than 1.0, which is in line with most industry design standards and should ensure cavern integrity and prevent loss of fluids to the surrounding rock mass. However, many of the SPR caverns currently have a P/D less than 1.0, or will likely have a low P/D after one or two full drawdowns. For these caverns, it is important to examine the structural integrity with more detail using geomechanical models. Finite-element geomechanical models have been used to determine the stress states in the pillars following successive drawdowns. By computing the tensile and dilatant stresses in the salt, areas of potential structural instability can be identified that may represent "red flags" for additional drawdowns. These analyses have found that many caverns will maintain structural integrity even when grown via drawdowns to dimensions resulting in a P/D of less than 1.0. The analyses have also confirmed that certain caverns should only be completely drawn down one time. As the SPR caverns are utilized and partial drawdowns are performed to remove oil from the caverns (e.g., for occasional oil sales authorized by the Congress or the President), the changes to the cavern volumes casused by these procedures must be tracked and accounted for so that an ongoing assessment of the cavern's drawdown capacity may be continued. A proposed methodology for assessing and tracking the available drawdowns for each cavern is presented in this report.
Geotechnical concerns arise due to the close proximity of the some of the caverns to each other (e.g., Caverns 15 and 17) or to the edge of the salt dome (e.g., Cavern 20). There are nine abandoned caverns, one of which collapsed (Cavern 7) in 1954 and another (Cavern 4) which is believed to be in a quasi-stable condition. This report provides explanations for these geotechnical concerns. The structural integrity of the pillar between BC-15 and 17 is examined. No salt fall is expected through 2045. However, the dilatant damaged area increases with time, especially, at the chimney area of BC-17. One drawdown leach for both caverns could be allowed if they are normally operated as a gallery, depressurized simultaneously. The possibility of a loss in integrity of BC-20 is examined in the salt between the dome edge and the cavern. The edge pillar is predicted to have experienced tensile stress since September 1999, but the small tensile stressed area is predicted to disappear in 2018 because BC-20 is filled fully with brine rather than oil since 2/7/2013. Even though BC-20 is no longer used as an SPR cavern, we need to continue monitoring the cavern integrity. BC-4 is also currently filled with brine and will not hold pressure at the wellhead. The cavern extends upward into the caprock and has no effective salt roof The results indicate that any sort of caprock roof collapse for BC-4 is not imminent but salt falls will likely occur from the near-roof portions of the cavern. The uncertainty due to salt falls illustrates the importance of continued monitoring of the area around BC-4 for behavior such as subsidence and tilt which may indicate a change in the cavern's integrity status.
The study described in this report involves heated and unheated pressurized slot testing to determine thermo-mechanical properties of the Tptpll (Tertiary, Paintbrush, Topopah Spring Tuff Formation, crystal poor, lower lithophysal) and Tptpul (upper lithophysal) lithostratigraphic units at Yucca Mountain, Nevada. A large volume fraction of the proposed repository at Yucca Mountain may reside in the Tptpll lithostratigraphic unit. This unit is characterized by voids, or lithophysae, which range in size from centimeters to meters, making a field program an effective method of measuring bulk thermal-mechanical rock properties (thermal expansion, rock mass modulus, compressive strength, time-dependent deformation) over a range of temperature and rock conditions. The field tests outlined in this report provide data for the determination of thermo-mechanical properties of this unit. Rock-mass response data collected during this field test will reduce the uncertainty in key thermal-mechanical modeling parameters (rock-mass modulus, strength and thermal expansion) for the Tptpll lithostratigraphic unit, and provide a basis for understanding thermal-mechanical behavior of this unit. The measurements will be used to evaluate numerical models of the thermal-mechanical response of the repository. These numerical models are then used to predict pre- and post-closure repository response. ACKNOWLEDGEMENTS The authors would like to thank David Bronowski, Ronnie Taylor, Ray E. Finley, Cliff Howard, Michael Schuhen (all SNL) and Fred Homuth (LANL) for their work in the planning and implementation of the tests described in this report. This is a reprint of SAND2004-2703, which was originally printed in July 2004. At that time, it was printed for a restricted audience. It has now been approved for unlimited release.
A finite element numerical analysis model has been constructed that consists of a mesh that effectively captures the geometries of Bayou Choctaw (BC) Strategic Petroleum Reserve (SPR) site and multimechanism deformation (M-D) salt constitutive model using the daily data of actual wellhead pressure and oil–brine interface location. The salt creep rate is not uniform in the salt dome, and the creep test data for BC salt are limited. Therefore, the model calibration is necessary to simulate the geomechanical behavior of the salt dome. The cavern volumetric closures of SPR caverns calculated from CAVEMAN are used as the field baseline measurement. The structure factor, A2, and transient strain limit factor, K0, in the M-D constitutive model are used for the calibration. The value of A2, obtained experimentally from BC salt, and the value of K0, obtained from Waste Isolation Pilot Plant salt, are used for the baseline values. To adjust the magnitude of A2 and K0, multiplication factors A2F and K0F are defined, respectively. The A2F and K0F values of the salt dome and salt drawdown skins surrounding each SPR cavern have been determined through a number of back analyses. The cavern volumetric closures calculated from this model correspond to the predictions from CAVEMAN for six SPR caverns. Therefore, this model is able to predict behaviors of the salt dome, caverns, caprock, and interbed layers. The geotechnical concerns associated with the BC site from this analysis will be explained in a follow-up paper.
The Department of Energy, in response to requests from the U.S. Congress, wishes to maintain an up-to-date table documenting the number of available full drawdowns of each of the caverns owned by the Strategic Petroleum Reserve. This information is important for assessing the SPR's ability to deliver oil to domestic oil companies expeditiously if national or world events dictate a rapid sale and deployment of the oil reserves. The evaluation of drawdown risks require the consideration of several factors regarding cavern and wellbore integrity and stability, including stress states caused by cavern geometry and operations, salt damage caused by dilatant and tensile stresses, the effect on enhanced creep on wellbore integrity, the sympathetic stress effect of operations on neighboring caverns. Based on the work over the past several years, a consensus has been built regarding the assessment of drawdown capabilities and risks for the SPR caverns. This report draws upon the recently Bayou Choctaw model upgrade and analyses to reevaluate and update the available drawdowns for each of those caverns. BC-18, 19, 101 and 102 are predicted to have conditional five available drawdowns remaining BC-15 and 17 have only one remaining drawdowns due to their proximity.
The three-dimensional finite element mesh capturing realistic geometries of the Bayou Choctaw site has been constructed using the sonar and seismic survey data obtained from the field. The mesh consists of hexahedral elements because the salt constitutive model is coded using hexahedral elements. Various ideas and techniques to construct finite element mesh capturing artificially and naturally formed geometries are provided. The techniques to reduce the number of elements as much as possible to save on computer run time while maintaining the computational accuracy is also introduced. The steps and methodologies could be applied to construct the meshes of Big Hill, Bryan Mound, and West Hackberry strategic petroleum reserve sites. The methodology could be applied to the complicated shape masses for various civil and geological structures.
A finite element numerical analysis model has been constructed that consists of a realistic mesh capturing the geometries of Bayou Choctaw (BC) Strategic Petroleum Reserve (SPR) site and multi - mechanism deformation ( M - D ) salt constitutive model using the daily data of actual wellhead pressure and oil - brine interface. The salt creep rate is not uniform in the salt dome, and the creep test data for BC salt is limited. Therefore, the model calibration is necessary to simulate the geomechanical behavior of the salt dome. The cavern volumetric closures of SPR caverns calculated from CAVEMAN are used for the field baseline measurement. The structure factor, A 2 , and transient strain limit factor, K 0 , in the M - D constitutive model are used for the calibration. The A 2 value obtained experimentally from the BC salt and K 0 value of Waste Isolation Pilot Plant (WIPP) salt are used for the baseline values. T o adjust the magnitude of A 2 and K 0 , multiplication factors A2F and K0F are defined, respectively. The A2F and K0F values of the salt dome and salt drawdown skins surrounding each SPR cavern have been determined through a number of back fitting analyses. The cavern volumetric closures calculated from this model correspond to the predictions from CAVEMAN for six SPR caverns. Therefore, this model is able to predict past and future geomechanical behaviors of the salt dome, caverns, caprock , and interbed layers. The geological concerns issued in the BC site will be explained from this model in a follow - up report .
This letter report provides a Bayou Choctaw (BC) Strategic Petroleum Reserve (SPR) well grading system based on the geomechanical simulation. The analyses described in this letter were used to evaluate the caverns’ geomechanical effect on wellbore integrity, which is an important component in the well integrity grading system recently developed by Roberts et al. [2015]. Using these analyses, the wellbores for caverns BC-17 and 20 are expected to be significantly impacted by cavern geomechanics, BC-18 and 19 are expected to be medium impacted; and the other caverns are expected to be less impacted.
This report summarizes the work performed in the prioritization of cavern access wells for remediation and monitoring at the Bayou Choctaw Strategic Petroleum Reserve site. The grading included consideration of all 15 wells at the Bayou Choctaw site, with each active well receiving a separate grade for remediation and monitoring. Numerous factors affecting well integrity were incorporated into the grading including casing survey results, cavern pressure history, results from geomechanical simulations, and site geologic factors. The factors and grading framework used here are the same as those used in developing similar well remediation and monitoring priorities at the Big Hill, Bryan Mound, and West Hackberry Strategic Petroleum Reserve Sites.
The three-dimensional finite element mesh capturing realistic geometries of Bayou Choctaw site has been constructed using the sonar and seismic survey data obtained from the field. The mesh is consisting of hexahedral elements because the salt constitutive model is coded using hexahedral elements. Various ideas and techniques to construct finite element mesh capturing artificially and naturally formed geometries are provided. The techniques to reduce the number of elements as much as possible to save on computer run time with maintaining the computational accuracy is also introduced. The steps and methodologies could be applied to construct the meshes of Big Hill, Bryan Mound, and West Hackberry strategic petroleum reserve sites. The methodology could be applied to the complicated shape masses for not only various civil and geological structures but also biological applications such as artificial limbs.
Oil leaks were found in wellbores of Caverns 105 and 109 at the Big Hill Strategic Petroleum Reserve site. According to the field observations, two instances of casing damage occurred at the depth of the interbed between the caprock bottom and salt top. A three dimensional finite element model, which contains wellbore element blocks and allows each cavern to be configured individually, is constructed to investigate the wellbore damage mechanism. The model also contains element blocks to represent interface between each lithology and a shear zone to examine the interbed behavior in a realistic manner. The causes of the damaged casing segments are a result of vertical and horizontal movements of the interbed between the caprock and salt dome. The salt top subsides because the volume of caverns below the salt top decrease with time due to salt creep closure, while the caprock subsides at a slower rate because the caprock is thick and stiffer. This discrepancy produces a deformation of the well. The deformed wellbore may fail at some time. An oil leak occurs when the wellbore fails. A possible oil leak date of each well is determined using an equivalent plastic strain failure criterion. A well grading system for a remediation plan is developed based on the predicted leak dates of each wellbore.
Oil leaks were found in wellbores of Caverns 105 and 109 at the Big Hill Strategic Petroleum Reserve site. According to the field observations, two instances of casing damage occurred at the depth of the interbed between the caprock bottom and salt top. A three dimensional finite element model, which contains wellbore element blocks and allows each cavern to be configured individually, is constructed to investigate the wellbore damage mechanism. The model also contains element blocks to represent interface between each lithology and a shear zone to examine the interbed behavior in a realistic manner. The causes of the damaged casing segments are a result of vertical and horizontal movements of the interbed between the caprock and salt dome. The salt top subsides because the volume of caverns below the salt top decrease with time due to salt creep closure, while the caprock subsides at a slower rate because the caprock is thick and stiffer. This discrepancy yields a deformation of the well. The deformed wellbore may fail at some time. An oil leak occurs when the wellbore fails. A possible oil leak date of each well is determined using the equivalent plastic strain failure criterion. A well grading system for a remediation plan is developed based on the predicted leak dates of each wellbore.
This report summarizes the work performed in developing a framework for the prioritization of cavern access wells for remediation and monitoring at the Big Hill Strategic Petroleum Reserve site. This framework was then applied to all 28 wells at the Big Hill site with each well receiving a grade for remediation and monitoring. Numerous factors affecting well integrity were incorporated into the grading framework including casing survey results, cavern pressure history, results from geomechanical simulations, and site geologic factors. The framework was developed in a way as to be applicable to all four of the Strategic Petroleum Reserve sites.
JAS3D, a three dimensional iterative solid mechanics code, has been used for structural analyses for the Strategic Petroleum Reserve system since the 1990s. JAS3D is no longer supported by Sandia National Laboratories, and has been replaced by Adagio. To validate the transition from JAS3D to Adagio, the existing JAS3D input decks and user subroutines for Bayou Choctaw and Big Hill models were converted for use with Adagio. The calculation results from the Adagio runs are compared to the JAS3D. Since the Adagio results are very similar to the JAS3D results, Adagio is judged to be performing satisfactorily.
This report addresses recent well integrity issues related to cavern 114 at the Big Hill Strategic Petroleum Reserve site. DM Petroleum Operations, M&O contractor for the U.S. Strategic Petroleum Reserve, recognized an apparent leak in Big Hill cavern well 114A in late summer, 2012, and provided written notice to the State of Texas as required by law. DM has since isolated the leak in well A with a temporary plug, and is planning on remediating both 114 A- and B-wells with liners. In this report Sandia provides an analysis of the apparent leak that includes: (i) estimated leak volume, (ii) recommendation for operating pressure to maintain in the cavern between temporary and permanent fixes for the well integrity issues, and (iii) identification of other caverns or wells at Big Hill that should be monitored closely in light of the sequence of failures there in the last several years.
The finite-element shock hydrodynamics code ALEGRA has recently been upgraded to include an X-FEM implementation in 2D for simulating impact, sliding, and release between materials in the Eulerian frame. For validation testing purposes, the problem of long-rod penetration in semi-infinite targets is considered in this report, at velocities of 500 to 3000 m/s. We describe testing simulations done using ALEGRA with and without the X-FEM capability, in order to verify its adequacy by showing X-FEM recovers the good results found with the standard ALEGRA formulation. The X-FEM results for depth of penetration differ from previously measured experimental data by less than 2%, and from the standard formulation results by less than 1%. They converge monotonically under mesh refinement at first order. Sensitivities to domain size and rear boundary condition are investigated and shown to be small. Aside from some simulation stability issues, X-FEM is found to produce good results for this classical impact and penetration problem.
Oil leaks were found in well casings of Caverns 105 and 109 at the Big Hill Strategic Petroleum Reserve site. According to the field observations, two instances of casing damage occurred at the depth of the interface between the caprock and top of salt. This damage could be caused by interface movement induced by cavern volume closure due to salt creep. A three dimensional finite element model, which allows each cavern to be configured individually, was constructed to investigate shear and vertical displacements across each interface. The model contains interfaces between each lithology and a shear zone to examine the interface behavior in a realistic manner. This analysis results indicate that the casings of Caverns 105 and 109 failed by shear stress that exceeded shear strength due to the horizontal movement of the top of salt relative to the caprock, and tensile stress due to the downward movement of the top of salt from the caprock, respectively. The casings of Caverns 101, 110, 111 and 114, located at the far ends of the field, are predicted to be failed by shear stress in the near future. The casings of inmost Caverns 107 and 108 are predicted to be failed by tensile stress in the near future.