The most popular methods used for solving transient heat conduction problems, like FEM and BEM require discretization of the domain or the boundary. The discretization problem escalates for unsteady issues, because an iterative process is required to solve them. An alternative to avoid the mentioned problem are parametric integral equations systems (PIES), which do not require classical discretization of the boundary and the domain, while being numerically solved. PIES have been previously used with success to solve steady-state problems. Moreover, they have been recently tested also with success for transient heat conduction problems, without internal heat sources. The purpose of this paper is to generalize PIES based on analytical modification of classical BIE for transient heat conduction with internal heat source and NURBS for boundary modeling. The obtained generalization of PIES is tested on examples, mostly with defined exact solution.