A numerical method for the gas production process in gas hydrate reservoirs