Transient shale gas flow model