Implementation of automatic heating control systems allows us to reduce heat consumption by 10% in residential areas and 40% in office and educational buildings. Currently, there are heating control systems, however, they are applied only to a single-level two-pipe heating system. Development of an imitation model of heat flows redistribution is necessary to synthesize the system of interconnected dynamic heating control of a building. Unlike existing solutions, this research work considers the problem of unbalanced heat flow in a multi-level hierarchical heating system. Calculation of convective heat transfer in the room assumes that the air temperature at any given time is the same throughout the entire room. When we determine heat transfer through walling, it is assumed that the walling or its part has the same temperature of the planes perpendicular to the direction of air flow. In this case, the heat transfer process is described by a one-dimensional heat transfer equation. The developed model of heating control systems allows us to connect the automatic control modules, change control algorithms at the compilation stage and in the system state during the simulation process. In comparison with possible analogue models based on AnyLogic or ANSYS modeling systems, the presented model is the model of controlled object. It is easily combined with models of automatic control units and considers the problem of imbalance of heat flows. An example of the functional scheme of the local temperature control system around one battery is considered. Implementation of developed imitation model makes it possible to ensure a new level of quality control of technological processes of production and consumption of power energy resources by using modern information technologies and synthesizing a system of interconnected dynamic heating control. Possibilities of such modeling are focused on development of the uninterrupted and high-quality heat supply system, maintaining energy-efficient operating modes, as well as actual economic effect. The model under consideration allows us to simulate redistribution of heat flows in different operating modes of the heating system.