The nonlinear propagation of self-gravitational shock structures (SGSSs) in a self-gravitating, super-dense, degenerate quantum plasma system (containing non-degenerate extremely heavy nuclei and ultra-relativistic degenerate electrons) has been investigated. The well-known reductive perturbation technique, which is valid in the small but finite amplitude limit, has been used to examine the nonlinear propagation of these SGSSs in such degenerate quantum plasma systems. The nonlinear dynamics of these SGSSs has been found to be governed by the Burgers equation, which is derived analytically and solved numerically in planar coordinates. These SGSSs in such plasma systems are shown to be formed due to the presence of a viscous force (which is the source of dissipation) acting on inertial extremely heavy nuclei species of the plasma system. The fundamental properties (viz., amplitude, width, speed, etc.) of these SGSSs are influenced in the ultra-relativistic limit. Our considered plasma model and the numerical analysis of the Burgers equation can be applied to astrophysical compact objects like neutron stars.