Quasi-Static Poroelastic Equations As A Symmetric Positive System And Its Numerical Approximation