We study a class of preconditioners to solve large-scale linear systems arising from fully implicit reservoir simulation.These methods are discussed in the framework of the auxiliary space preconditioning method for generality.Unlike in the case of classical algebraic preconditioning methods,we take several analytical and physical considerations into account.In addition,we choose appropriate auxiliary problems to design the robust solvers herein.More importantly,our methods are user-friendly and general enough to be easily ported to existing petroleum reservoir simulators.We test the efciency and robustness of the proposed method by applying them to a couple of benchmark problems and real-world reservoir problems.The numerical results show that our methods are both efcient and robust for large reservoir models. We study a class of preconditioners to solve large-scale linear systems arising from fully implicit reservoir simulation. These methods are discussed in the framework of the auxiliary space preconditioning method for generality .Unlike in the case of classical algebraic preconditioning methods, we take several analytical and physical considerations into account. In addition, we choose appropriate auxiliary problems to design the robust solvers herein. More importantly, our methods are user-friendly and general enough to be easily ported to existing petroleum reservoir simulators. We test the efciency and robustness of the proposed method by applying them to a couple of benchmark problems and real-world reservoir problems. numerical results show that our methods are both efcient and robust for large reservoir models.