Numerical models for the simulation of bleed-off control on post-injection seismicity in ‎enhanced geothermal systems

This dataset comprises the input files for numerical simulation of stimulating an enhanced geothermal system and the associated risk of induced seismicity. Numerical models explore the hydromechanical mechanisms of post-injection induced seismicity and the potential impact of bleed-off to control seismicity. The models are 2D plane strain including a homogeneous and isotropic granitic reservoir rock intersected by an interconnected fracture network and a distant fault to the wellbore. The fault in our numerical models is either hydraulically isolated from the fracture network, i.e., our “basecase” scenario, or intersected by them, a scenario referred to as “connected”. Post-injection protocols to control seismicity include well shut-in or bleed-off. During bleed-off, the wellhead is opened to impose hydrostatic pressure conditions along the well over a period of either 24 h or 1 h. The ‎hydrostatic pressure profile is maintained afterwards. Different injection periods are considered to inspect the effect of this parameter on the effectiveness of bleed-off in controlling post-injection seismicity. Following this brief explanation, the dataset includes seven models provided in seven distinct folders: “Basecase_Inj250h_shutin”, basecase scenario of 250 h injection followed by well shut-in, “Basecase_Inj250h_bleedoff_1h”, basecase scenario of 250 h injection followed by well bleed-off over 1 h, “Basecase_Inj250h_bleedoff_24h”, basecase scenario of 250 h injection followed by well bleed-off over 24 h, “Connected_Inj84h_shutin”, connected scenario of 84 h injection followed by well shut-in, “Connected_Inj84h_bleedoff_1h”, connected scenario of 84 h injection followed by well bleed-off over 1 h, “Connected_Inj200h_shutin”, connected scenario of 200 h injection followed by well shut-in, “Connected_Inj200h_bleedoff_1h”, connected scenario of 200 h injection followed by well bleed-off over 1 h. In each folder, there is a file with the name of the folder ending as “_gen.dat” which contains the input data of the model, including material properties, initial and boundary conditions and the time intervals. There is also a file ending as “_gri.dat” that includes the information on the mesh. The file “root.dat” includes the name of the model. To run the simulation, execute the Code_Bright executables “Cb_v9_3.exe” or “Cb_v9_1_3par4”, i.e., the parallel version of the code running on 4 CPU cores, in each folder.