A semi-analytical model for horizontal-well productivity in shale gas reservoirs: Coupling of multi-scale seepage and matrix shrinkage