Dynamic behavior of the propagation of incoherent optical spatial solitons in a nonlocal nonlinear medium is investigated. By means of the Hirota method, we obtain the soliton solutions, based on which we find that |ψ| is positively related to n 0 , but inversely related to ω and κ, whereas n is independent of them, with ψ as the slowly varying amplitude of the beam, n as the refractive index change, n 0 , ω, and κ as the unperturbed refractive index, frequency of the propagating beam, and beam intensity distribution, respectively. Head-on and bound-state soliton interactions are both given, and one interaction period decreases with ω and κ increasing in the bound state one. With the external perturbations taken into consideration, the associated chaotic motions of the perturbed model are studied, and the corresponding power spectra and phase projections are obtained. Both the weak and developed chaotic states are investigated, and the difference between them roots in the relative magnitude of nonlinearities and perturbations. Such chaotic motions can be weakened via increasing ω and κ or decreasing n 0 . Periodic motion is obtained with the nonlinearities and perturbations balanced.