To perform a similar type of surface-subsidence analysis, consider the following recommendations:
Use a hex-dominant mesh.
If possible, use sweep meshing in the reservoir and in the domain above and below it.
Because the global stiffness matrix is unsymmetric, use the Newton-Raphson solver with the unsymmetric option (NROPT,UNSYMM) to avoid convergence issues.