In this Letter, we report the results of constraining the dynamical dark energy with a divergence-free parameterization, w(z) = w(0) + w(a)(In(2+z)/1+z - In2), in the presence of spatial curvature and massive neutrinos, with the 7-yr WMAP temperature and polarization data, the power spectrum of LRGs derived from SDSS DR7, the Type la supernova data from Union2 sample, and the new measurements of H-0 from HST, by using a MCMC global fit method. Our focus is on the determinations of the spatial curvature, Omega(k), and the total mass of neutrinos, Sigma m(v) in such a dynamical dark energy scenario, and the influence of these factors to the constraints on the dark energy parameters, w(0) and w(a). We show that Omega(k) and Sigma m(v) can be well constrained in this model; the 95% CL limits are: -0.0153 < Omega(k) < 0.0167 and Sigma m(v) < 0.56 eV. Comparing to the case in a flat universe, we find that the error in w(0) is amplified by 25.51%, and the error in wa is amplified by 0.14%; comparing to the case with a zero neutrino mass, we find that the error in w(0) is amplified by 12.24%, and the error in w(a) is amplified by 1.63%. (C) 2012 Elsevier B.V. All rights reserved.