On two-phase flow in karstic geometry: modeling, analysis and numerical simulations