Sensitivity analysis of flow in unsaturated heterogeneous porous media: Theory, numerical model, and its verification