Application of high order seismic modelling methods to 2D elastic random media

A seismic modeling method of high accuracy is required to quantify partly subtle effects of random media wave propagation. The modeling method must also allow coarse sampling to include as much statistics as possible in discrete data. Both requirements favour high order methods. But can we use high order difterentiators in random media? And if yes, which specific method is best?