This paper deals with a time stepping method for the finite element simulation of two phase flow hydraulic and mechanical (H2M) coupled processes in porous media, which is a common phenomenon in geological applications such as CO₂ storage facilities. The computation task arising from the numerical modeling of H2M coupled processes in such real geological application is intensive. Therefore, the high performance computing is of interest to the corresponding researchers. In the present study, we present a time stepping method with PI (proportional and integral feedback) automatic control to improve the computation efficiency of the modeling of H2M coupled processes. We apply the PI control to solve the nonlinear coupled partial differential equations with a first order finite difference scheme for time discretization and the Picard method for linearization. The efficiency of the present method is demonstrated by applying it to a CO  storage benchmark.