An integrated model with stable numerical methods for fractured underground gas storage. (20th March 2023)