Modeling Horizontal Salt Cavern Leaching in Bedded Salt Formations

Abstract

Summary The leaching of a salt cavern will trigger a series of rock-fluid interactions, including salt rock dissolution, cavity expansion, and brine transport caused by convection, turbulence, and diffusion effects. These interactions have influences on one another. The primary objectives of this study include developing a 3D multiphysical coupled model for horizontal salt cavern leaching and quantifying these interactions. The species transport equation and standard κ-ε equation were combined to describe the brine transport dynamics within the cavity. Based on the velocity and concentration distribution characteristics predicted, the interface movement equation implemented with mesh deformation techniques was applied to describe the cavity expansion. Next, the Volgograd cavern monitored data were collected for model validation. The predicted results are consistent with the field data. The average relative errors are 11.0% for brine displacing concentration and 4.5% for cavity volume. The results suggest that the cavity can be divided into three regions, including the main flow region, circulation region, and reflux region. The results also suggest that the brine concentration distribution is relatively uniform. With the dissolution threshold angle and anisotropic dissolution rates considered, the resultant cavity cross section is crown top and cone bottom. The results also show that the cavity can be divided into dissolution and erosion sections according to its position relative to the injection point.