In this article, we construct a new combined characteristic mixed finite element procedure to simulate the incompressible wormhole propagation. In this procedure, we use the classical mixed finite element method to solve the pressure equation and a modified mass-preserving characteristic finite element method for the solute transport equation, and solve the porosity function straightly by the given concentration. This combined method not only keeps mass balance globally but also preserves maximum principle for the porosity. We considered the corresponding convergence and derive the optimal L2-norm error estimate. Finally, we present some numerical examples to confirm theoretical analysis.