Solutions of the deep bed filtration equations are sought for situations involving spatially varying initial porosites and temporally varying superficial velocities and inlet concentrations. Results of a general nature are obtained for both a quasistatic model and a diffusionless model. In the latter case it is shown that one numerical solution can sometimes be used to infer information about a class of problems.