In this paper, we systematically investigated the geometric quantum discord of Heisenberg spin chain models in a Non-Markovian environment using a Non-Markovian quantum state diffusion method proposed by Diosi and Gisin in 1998. Geometrical quantum discord was used to describe quantum correlations as proposed by Dakic et al. In the first place the reduced density matrix of the system was calculated by the Non-Markovian quantum state diffusion master equation. Afterward, the reduced density matrix was brought into the quantum discord formula, thus achieving a numerical precise simulation of the geometric quantum discord of the Heisenberg spin chain. At length the evolution of the geometric quantum discord of Heisenberg spin chain models with various parameters vs time was discussed in this paper, using the maximally entangled state |ψAB〉=12|11〉+00e as the initial state of Heisenberg spin chain model of the system. According to the numerical simulation results, the environmental correlation coefficient γ, parameters J, parameters a and parameters η can all influence the evolution of the system geometric quantum discord dynamics to different degrees. When the environmental correlation coefficient γ is very small, the geometric quantum discord shows a significant upward trend, showing that non-Markovian environmentsplay a positive role on the geometric quantum discord of the system. Also large parameters J, a and η also play a positive role on the geometric quantum discord of the system. The present study has a certain role and significance for increasing the geometric quantum discord of the system, providing a theoretical basis for the experimental researchers in practice.