In this paper, we intend to mitigate absolute accelerations and displacements in the low-frequency regime of multiple-degrees-of-freedom fuel storage tanks subjected to stochastic seismic excitations. Therefore, we propose to optimize a finite locally resonant metafoundation equipped with massive resonators and fully nonlinear hysteretic devices. The optimization process takes into account the stochastic nature of seismic records in the stationary frequency domain; the records are modelled with the power spectral density S0 and modified with a Kanai-Tajimi filter. Moreover, the massive superstructure of a fuel storage tank is also considered in the optimization procedure. To optimize the nonlinear behaviour of dampers, we use a Bouc-Wen hysteretic model; the relevant nonlinear differential equations are reduced to a system of linear equations through the stochastic equivalent linearization technique. The optimized system is successively verified against natural seismic records by means of nonlinear transient time history analyses. Finally, we determine the dispersion relations for the relevant periodic metafoundation.