To suppress the fluctuation from both load side and intermittent renewable energy (IRE), nuclear power plants (NPPS) should be operated in load-following mode to improve economic competitiveness. The modular high temperature gas-cooled reactor (MHTGR) belongs to the category of small nuclear reactor (SMRs) and is suitable for load-following by the virtue of online refueling ability and inherent safety. To realize economies of scale for MHTGR, multi-modular scheme that multiple nuclear steam supply system (NSSS) modules are connected in parallel providing superheated steam for common turbine is recommended to achieve desired power ratings. However, because of the large heat capacity in the pebble-bed of MHTGR and thermal coupling of different NSSSs through common secondary loop fluid network, the current control strategy which suppresses the nuclear power, coolant temperatures measurement from their set-points without considering thermal dynamic of NSSS itself, may not favorable for heat transfer in the NSSS. To improve the load-following ability, a multivariable dynamic matrix control (DMC) is constituted to dynamic compensate the thermal energy variation of NSSS. The implementation of the DMC has a typical cascade structure, where DMC revises the set-points of NSSS module in outer loop and the existing PID control law is adopt for stabilization in inner loop. Numerical results show that this cascade DMC can improve the transient of thermal power under power maneuvering, and can also attenuate the nuclear power, helium flowrate set-points and feed-water temperature step disturbance.