We report the numerical and experimental studies of the two-dimensional Brillouin gain spectrum (BGS) distribution deformation induced by the self-phase modulation in the Brillouin optical time domain reflectometry (BOTDR) with a 20.6 km sensing distance. The BGS distribution deformation is investigated by analyzing the evolution of the point spread function along the fiber in the two-dimensional model of the BOTDR. In the simulation and experimental results, the specific deformation degree of the BGS distribution induced by the self-phase modulation is related to the pump pulse profile, pump pulse peak power, BGS demodulation method, and detected scattered light component. By comprehensively analyzing the evolution of the point spread function induced by the self-phase modulation and using the image deconvolution, a typical BOTDR sensor with a 25 ns pump pulse reaches the 20 cm spatial resolution over the 20.6 km sensing fiber.