We report on a refined macroscopic model for slightly compressible gas slip flow in porous media developed by upscaling the pore-scale boundary value problem. The macroscopic model is validated by comparisons with an analytic solution on a two-dimensional (2-D) ordered model structure and with direct numerical simulations on random microscale structures. The symmetry properties of the apparent slip-corrected permeability tensor in the macroscale momentum equation are analysed. Slip correction at the macroscopic scale is more accurately described if an expansion in the Knudsen number, beyond the first order considered so far, is employed at the closure level. Corrective terms beyond the first order are a signature of the curvature of solid–fluid interfaces at the pore scale that is incompletely captured by the classical first-order correction at the macroscale. With this expansion, the apparent slip-corrected permeability is shown to be the sum of the classical intrinsic permeability tensor and tensorial slip corrections at the successive orders of the Knudsen number. All the tensorial effective coefficients can be determined from intrinsic and coupled but easy-to-solve closure problems. It is further shown that the complete form of the slip boundary condition at the microscale must be considered and an important general feature of this slip condition at the different orders in the Knudsen number is highlighted. It justifies the importance of slip-flow correction terms beyond the first order in the Knudsen number in the macroscopic model and sheds more light on the physics of slip flow in the general case, especially for large porosity values. Nevertheless, this new nonlinear dependence of the apparent permeability with the Knudsen number should be further verified experimentally.< Réduire