Numerical simulation of ground surface subsidence due to coal-bed methane extraction