К вычислению ковариационной функции интенсивности световой волны в среде с крупномасштабными случайными неоднородностями. Кумулянтный подход