Soil dynamic parameters are crucial calculation parameters for the seismic response of soil. However, there is currently insufficient theoretical research on soil dynamic characteristics in Heilongjiang Province. To address this, we analyzed data on the dynamic shear modulus ratio and damping ratio of 286 groups of clay, clayey soil, and sandy soil at different burial depths in the region. We derived the characteristic parameters of each group of soil dynamic parameters and statistically obtained a fitting formula for the variation of these parameters with soil burial depth. Using this formula, we calculated the dynamic shear modulus ratio and damping ratio data of different soil types at different burial depths and verified the rationality of the fitting formula. We established six soil seismic response models using soil types with abundant data and input typical seismic motions in Heilongjiang Province. We calculated the fitting parameters a₁, a₂, E(λ_max), E(M), etc. of the fitting formula and studied the influence of the variability of characteristic parameters on peak ground acceleration, characteristic period, and plateau value of ground response spectrum by scaling the fitting parameters. Results show that with the increase of a₁ and E(λ_max) the damping ratio increases, resulting in a decrease in the site amplification factor F_PGA and an increase in the characteristic period T_g. The increase of a₂ and E(M) will lead to a decrease in the damping ratio, resulting in an increase in F_PGA and a decrease in T_g. The plateau value β has a large variability, fluctuating within a certain range without a definite increasing or decreasing trend. Variability has the greatest impact on F_PGA, followed by its impact on T_g, and then on β. The soil seismic response is the most sensitive to the variability of a₁. As the intensity of the input seismic motion increases, the error increases. Overall, this study provides necessary supplements to the research on soil dynamic parameters in Heilongjiang Province.