A 2-dimensional heterogeneous reactor model was developed and simulated for a tube reactor of packed bed where the steam-$CO_2$ combined reforming reaction of natural gas proceeded to produce synthesis gas. Under the reactor feeding rate, 45 $Nm^3$/h, of the reactant gas stream, the 2-dimensional heterogeneous reactor model showed the similar results to those from the ASPEN simulator although there were some discrepancies between the two in the temperature and the $H_2$/CO ratio of the reformed gas at the reactor exit. The calculated enthalpy difference between the reformed gas at the reactor exit and the reactant gas fed to the reactor was closely correspondent to the total amount of heat transferred to the reactor interior from the furnace. This supports that the 2-dimensional heterogeneous reactor model was reasonably established and the numerical solution was properly obtained.