∫ ∫ ∫
It is convenient to build a THM formulation by Galerkin weighted residual method for modeling the behavior of reservoir formation under elevated temperature. With the Onsager reciprocal principle, it is generally possible to have a finite element formulation for the THM problem with a symmetrical coefficient matrix. However, this will require many more parameters and will result in many additional computation tasks. The asymmetrical coefficient matrix has resulted from the simplification described in this chapter. However, there are currently several numerical techniques, such as mathematical programming, that can be effectively used for solving the set of equations that have an asymmetrical coefficient matrix. The computer implementation of the model presented is an important aspect for the application of the model, and it will be the task for the next stage of this work.