The results of numerical simulation of the nonlinear evolution of an acoustic pulse propagating above an earthquake epicenter with preset energy parameters are presented. The interaction with charged plasma particles in the lower ionosphere is taken into account. The interaction between neutral and charged components is described by the diffusion mechanism of the disturbance of the density of charged ionospheric plasma components. The simulation was carried out for three values of earthquake magnitude M = 5; 6; and 7. It is ascertained that the numerical simulation results agree with analytical estimates made during the study of main regularities of seismic events that induce weak nonlinear effects in the atmosphere.