The possibility is here explored to use an ‘equivalent’ homogeneous configuration to simulate 1D seismic response of heterogeneous engineering-geological bodies when relatively weak seismic impedance contrasts internal to the bodies (it was assumed a shear wave velocity variation between the alternating layers equal to 150 m/s) only exist above the seismic bedrock. This equivalent configuration is obtained by considering an equivalent Vs value the harmonic average of the actual Vs values and a linear combination of G/G0 and D curves relative to the lithotechnical components present in the actual configuration. To evaluate feasibility of this approach, a wide set of numerical simulations was carried out by randomly generating subsoil layering including sequences of alternating thin layers of geotechnical units (e.g. sands and clays) each characterized by a characteristic nonlinear curve. Outcomes of these simulations are compared with those provided by considering a single homogeneous layer characterized by equivalent nonlinear curves obtained as a weighted average of the original curves. By comparing the heterogeneous and the homogeneous columns seismic response in terms of amplification factors and fundamental period, the results confirm the possibility to model a 1D column characterized by a generic lithostratigraphic succession with an equivalent one without introducing significative errors that, at least for the studied cases, do not exceed the 6%. This conclusion is substantially confirmed by extending the comparison to a real case, i.e., the 113 m-thick heterogeneous soil profile at Mirandola site (Northern Italy), presented in the last part.