Summary: Two-dimensional hyperbolic heat conduction problems are investigated by using the hybrid numerical scheme. The thermal wave of such problems propagates with a finite velocity. Thus numerical oscillations in the vicinity of the thermal wave front can be observed, and a hybrid numerical method is presented, to reduce these oscillatory magnitudes. This method is that the time-dependent terms in the governing differential equations are removed by using the Laplace transform technique, and then the control volume method is used to discretize the space domain in the transform domain. The key of the present method is the selection of the shape functions. Various examples with the irregular geometry are given.