Numerical description of shale gas desorption stages