Приведены интегральные представления для трехмерного поля излучения в плоскослоистой и двухслойной средах. Для двухслойной среды исследовано поведение поля в окрестности разрыва диэлектрической проницаемости, где общее интегральное представление неприменимо, и получены хорошо сходящиеся представления поля в этой области. Приведены выражения для тангенциальной составляющей поверхностной волны через аналитическую функцию двух переменных при любых параметрах разрыва. Получено асимптотическое разложение поля на малых расстояниях от антенны, расположенной на разрыве сред, что позволяет устранить сингулярность при численном моделировании поля. Таким образом, сделаны необходимые шаги к разработке модели тока, учитывающей влияние подстилающей поверхности, что, в свою очередь, является новым этапом в развитии алгоритма интерпретации данных георадиолокации методом обратной задачи.