A volume averaging theory (VAT) established in the field of fluid-saturated porous media has been successfully exploited to derive a general set of bioheat transfer equations for blood flows and its surrounding biological tissue. A closed set of macroscopic governing equations for both velocity and temperature fields in intra- and extravascular phases has been established, for the first time, using the theory of anisotropic porous media. Firstly, two individual macroscopic energy equations are derived for the blood flow and its surrounding tissue under the thermal non-equilibrium condition. The blood perfusion term is identified and modeled in consideration of the transvascular flow in the extravascular region, while the dispersion and interfacial heat transfer terms are modeled according to conventional porous media treatments. It is shown that the resulting two-energy equation model reduces to Pennes model, Wulff model and their modifications, under appropriate conditions. Subsequently, the two-energy equation model has been extended to the three-energy equation version, in order to account for the countercurrent heat transfer between closely spaced arteries and veins in the circulatory system and its effect on the peripheral heat transfer. This general form of three-energy equation model naturally reduces to the energy equations for the tissue, proposed by Chato, Keller and Seiler. Controversial issues on blood perfusion, dispersion and interfacial heat transfer coefficient are discussed in a rigorous mathematical manner. (c) 2007 Published by Elsevier Ltd.