Carbon dioxide disposal into deep aquifer has been an important venue to trap excess gas emission which causes global warming. In the CO2 sequestration process, CO2 is captured from the point source and injected into the saline aquifer deep enough where CO2 reaches the supercritical state and it has a very high density compared to gaseous state. This process is described by the two-phase two-component model, which involves two nonlinear time dependent advection-diffusion equations. The difficulty lies in the injection phase when the advection terms highly dominate over the diffusion terms. Discontinuous Galerkin (DG) methods, which are famous for the properties of high order accuracy, locality and locally mass conservation, have proved to be promising for advection dominated transport equations. I develop a new fully implicit fully coupled DG method called “partial upwind” method, to discretize the equations. For time discretization, it uses the backward Euler method to allow large time steps. For space discretization, it uses the usual interior penalty DG discretization for the elliptic terms and the upwind for part of the advection terms. The other part of the advection terms are handled specially for stabilization purpose. Numerical simulations show that the new DG method works well for the CO2 sequestration problems in homogenous porous media and has shown great potential in heterogenous porous media. I also compare the new method with the primal interior penalty DG method and show that the new method is superior to the usual DG method for some subsurface fluid flow problems. Though DG methods perform well for the CO2 sequestration problem, they are indeed more costly than traditional numerical methods. The first order finite volume (FV) method, on the other hand, is very efficient. A new coupled finite volume and discontinuous Galerkin method, which uses the accuracy of DG methods on some parts of the domain and the efficiency of FV methods everywhere else to reduce the computational cost, is also studied for the time-dependent advection-diffusion equations. Theoretical and numerical results show that the new coupled method converges and can be both accurate and efficient at the same time for some typical examples. We want to apply the new coupled method to the CO2 sequestration problems in the future.