This talk presents results from our recent efforts for reviving the 2-field (fluid pressure and solid displacement) approach for numerically solving poroelasticity problems on quadrilateral meshes. The Darcy equation is solved for fluid pressure by the weak Galerkin finite element methods, which establish the discrete weak gradient and numerical velocity in the Arbogast-Correa spaces. The elasticity equation is solved for solid displacement by the enriched Lagrangian elements, which were motivated by the Bernardi-Raugel elements for Stokes flow. These two types of solvers are coupled through the implicit Euler temporal discretization to solve poroelasticity. Numerical experiments on two widely tested benchmarks will be presented to demonstrate the new solvers are locking-free.  We discuss also extension to 3-dim and implementation in deal.II. This is based on a series of joint work with several collaborators.