Труды Кольского научного ценра РАН. № 6, вып.1. 2019 г.

измеренных компонент поля, умноженным на геометрические коэффициенты, с учетом значений силы тока в источнике выполнялись вычисления значений кажущегося сопротивления для измеряемых компонент поля. Углубленная обработка предполагает определение погрешности измерений сигнала путем вычисления статистических оценок (среднего значения, дисперсии и среднеквадратичного отклонения) по нескольким фрагментам записи на данной частоте. Для этого рассматривались синхронные записи силы тока в источнике и измеряемых компонент поля для каждого из используемых фрагментов записи. В каждом из фрагментов выполнялся спектральный анализ авто- и взаимных корреляционных функций наблюдаемых компонент поля. Погрешность измерений амплитуд наблюдаемых компонент силы тока в источнике и фазовых сдвигов между компонентами поля определялась по среднеквадратичному отклонению (СКВО) от среднего значения. Среднее значение амплитуды и СКВО рассчитывались по трем фрагментам записи поля на данной частоте колебаний тока в источнике. Начало и конец фрагмента длительностью в 1/3 длительности записи сигнала — от 15 мин (для частоты сигнала 0,381 Гц) до 6 мин (для частоты сигнала 194,2 Гц) — на каждой частоте определялись с точностью до 1 с путем синхронизации записей силы тока в источнике и компонент поля в измерительной точке. Фрагменты записей тока и сигналов поля прореживались до одной частоте дискретизации, равной: 1 кГц — для записей сигналов на частотах 94,2-194,2 Гц; 500 Гц — для частот 38,2-64,2 Гц; 250 Гц — для частот 1,942-19,42 Гц; 62,5 Гц — для частот 0,382-0,942 Гц. Децимация выполнялась путем замены участка записи из ^-отсчетов средним значением на участке, где к — шаг децимации. Определение амплитуд в каждом из трех фрагментов выполнено по Фурье-спектру автокорреляционной функции с окном Блэкмана — Харриса с 50-процентным перекрытием. Фазовые сдвиги между компонентами поля определялись по Фурье-спектрам взаимнокорреляционных функций с тем же окном. Пример спектральной оценки амплитуды тока на частоте 9,422 Гц для каждого источника поля (BJ1 J1-401 и J1-153) приведены на рис. 4 и 5. Спектры рассчитаны методом быстрого преобразования Фурье (БПФ) [14]. Таким образом, сила тока в линиях определялась с относительной погрешностью не более 0,5 %. Обработка записей сигналов от измерительных датчиков проводилась по тем же фрагментам, с той же частотой дискретизации и теми же окнами, что и обработка записей тока. Синхронизация фрагментов записей сигналов и тока выполнялась по GPS-меткам с погрешностью не более чем 0,5 с. Поскольку удаление от центров питающих линий L1 и L2 для ряда точек наблюдений не превышает пяти длин самих линий (рис. 2), то нормальное поле источника в дальней зоне вычислялось с учетом конечных размеров и формы источника. Нормальное поле длинной линии, лежащей на плоской границе земля — воздух, определялось путем интегрирования поля горизонтального электрического диполя (ГЭД) в дальней зоне по контуру питающей линии. Поле ГЭД в дальней зоне, записанное в локальной системе координат источника, определяется [15] ортом x(s), параллельным касательному вектору ds, ортом z(s). направленным по вертикали, и ортом j(s ), являющимся векторным произведением единичных векторов z(s) и х(.ѵ): J'(.v) = z ( s ) х .V(.V). 26