El método de diferencias finitas se utiliza ampliamente en la simulación numérica de ondas sísmicas, la imagen y la inversión de formas de onda. En el método de diferencias finitas, el operador de diferencia finita se utiliza para reemplazar aproximadamente al operador diferencial, que se puede obtener truncando la serie de convolución espacial. Las propiedades de la función de ventana truncada, como los lóbulos principales y laterales de la respuesta de amplitud de la función de ventana, determinan la precisión de las diferencias finitas, lo que afecta significativamente los resultados de imagen e inversión sísmica. Aunque la dispersión numérica es inevitable en este proceso, se puede suprimir de manera más efectiva utilizando operadores de diferencias finitas de mayor precisión. En este artículo, utilizamos el algoritmo del rebaño de krill, en contraste con el PSO estándar y el CDPSO (una variante de PSO), para optimizar el operador de diferencias finitas. Los resultados de la simulación numérica verifican que el algoritmo del rebaño de krill tiene un buen rendimiento en la mejora de la precisión del operador diferencial.