Prediction of ground motion by earthquakes is still the main issue for seismologists, geo-technicians, engineers, and researchers. In this paper k-means algorithm has been used to disaggregate seismic parameters to evaluate their inter-correlations. A goal is to quantify the weights and effects of each parameter with respect other ones. From database, about 4900.0 data, divided in 22 categories, have been collected. The main divisions regard the wave components in horizontal and vertical axis, and the soil characteristics. The studied seismic zone is the “Norpirenaica oriental”, placed at Pyrenees area between Spain and France, classified as a very high seismic hazard. Numerical and analytical analyses have been carried out to implement the algorithm. These analyses could quantify the role of the sand horizontal stratigraphy, the non-linear response, and the elasticity of the soil. Results are plotted in stochastic distributions and elastic spectra accelerations. Rigorously, results are valid only for the studied seismogenic zone under predefined constrictions and range.