On Application of the Weak Galerkin Finite Element Method to a Two-Phase Model for Subsurface Flow