An integrated model with stable numerical methods for fractured underground gas storage