Scientific journal
Advances in current natural sciences
ISSN 1681-7494
"Перечень" ВАК
ИФ РИНЦ = 0,775

NOISE MODELING IN HIGH-RATE REALTIME GNSS TIMESERIES

Pupatenko V.V. 1 Shestakov N.V. 2, 3
1 Yu.A. Kosygin Institute of Tectonics and Geophysics Far Eastern Branch
2 Far Eastern Federal University
3 Institute of Applied Mathematics Far Eastern Branch
In this paper the results of noise modeling in high-rate GNSS timeseries are presented. GNSS data processing was performed using Precise Point Positioning (PPP) approach in realtime mode. Initial data for noise modeling includes 10–13 days of 1 s observations from 128 global distributed continuous GNSS stations belong to IGS and CORS networks. Modified sidereal filtering was applied to coordinate timeseries for reducing multipath error. The LRTGNM (Low Real-Time GNSS Noise Models), MRTGNM (Median Real-Time GNSS Noise Models) and HRTGNM (High Real-Time GNSS Noise Models) was developed based on statistical distribution of GNSS timeseries noise. The models was compared with Peterson’s seismic noise models NHNM (New High Noise Model) и NLNM (New Low Noise Model) as well as with high-quality accelerometers noise models AHNM (accelerometric high noise model) and ALNM (accelerometric low noise model). It is shown that for periods 300 s and less the both seismic instruments are more sensitive than the GNSS receiver.
GPS
GNSS
noise models
Precise Point Positioning
real-time GNSS processing

Данные глобальных навигационных спутниковых систем (ГНСС), в первую очередь американской системы GPS, уже более 25 лет успешно используются в геодинамических исследованиях, наблюдениях за вулканической активностью и в сейсмологии [14].

В результате так называемой кинематической обработки ГНСС-измерений можно формировать временные серии (ряды) координат геодезических пунктов с малым шагом дискретизации 1 с и менее [2, 4], т.е. практически мгновенных положений ГНСС-станций. С их помощью получают записи сейсмических волн, генерируемых сильными землетрясениями. Таким образом, ГНСС-оборудование может быть использовано в качестве сейсмометра [1].

Основные преимущества ГНСС-технологий перед классическими сейсмическими инструментами заключаются в возможности регистрации сколь угодно сильных колебаний земной поверхности любых периодов без искажения фазы и амплитуды [1], основной недостаток – высокий уровень шума, обусловленный сантиметровой точностью определения мгновенных координат [4, 8].

В координатных рядах, полученных обработкой ГНСС-измерений, можно выделить две составляющие. Одна из них – полезный сигнал – обусловлена перемещениями антенны вследствие, например, природных или техногенных подвижек земной поверхности или прохождения сейсмических волн. Другая составляющая – шум – совместное действие случайных, систематических и грубых ошибок на результаты определения координат.

Исследование характеристик шума временных рядов ГНСС-координат имеет большое практическое значение, так как может дать информацию о его причинах и позволяет судить о надежности и точности выделения «полезного сигнала» из исходного координатного ряда. Наиболее просто и удобно определять характеристики шума, используя ряды, которые полезного сигнала не содержат.

Обзор более 15 работ, посвященных моделированию шума в рядах координат, получаемых обработкой данных ГНСС-измерений, можно найти в [11]. Во всех этих работах обработка данных ГНСС-измерений выполнена методом относительного позиционирования, т.е. определялись взаимные перемещения пары спутниковых приемников, установленных в стабильной точке с известными координатами (база) и определяемой подвижной точке (ровер). Вектор, компонентами которого являются разности координат базы и ровера, называется базовой линией, а расстояние между базой и ровером, определяемое модулем этого вектора, называется длиной базовой линии.

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

В настоящей работе выполнено моделирование шума рядов координат, получаемых при обработке ГНСС-данных методом Precise Point Positioning (PPP) [15] в режиме реального времени. Этот метод позволяет получать мгновенные значения абсолютных координат без привязки к базовым станциям и является предпочтительным для получения волновых форм землетрясений и основным при обработке данных ГНСС-измерений в режиме реального времени [8]. Построенные модели шума рядов координат сопоставлены со стандартными моделями сейсмического шума [5, 13], что позволяет делать выводы о возможности применения ГНСС-аппаратуры для решения сейсмологических задач и/или совместного использования с сейсмическими инструментами.

Материалы и методы исследования

В качестве исходных данных были использованы наблюдения с шагом дискретизации 1 с, полученные на 128 станциях (рис. 1) ГНСС-сетей IGS [10] и CORS [7]. Продолжительность записей составила от 10 до 13 суток. Данные измерений на всех станциях не содержат пропусков более 15 минут. В целом выбранные станции имеют достаточно равномерное пространственное распределение, со сгущением на территории США и с некоторым недостатком станций на Азиатском континенте и в Антарктиде (т.к. данные многих станций там содержали значительные пропуски или имели шаг дискретизации более 1 с). Данные ГНСС-измерений станций сетей IGS и CORS доступны на официальных сайтах [12, 7].

pup1.tif

Рис. 1. ГНСС-станции, использованные для построения моделей шума

Для упрощения работы с данными реальных ГНСС-наблюдений, инструментальные данные за выбранный временной интервал и необходимая сопутствующая информация сначала накапливались в специальном архиве, после чего проводилась их обработка путем имитации вычислений в режиме реального времени. Обработка была осуществлена в известной программе BKG Ntrip Client [3] с использованием поправок часов и орбит спутников CLK91, предоставляемых сервисом BKG.

Для уменьшения влияния ошибок, связанных с переотражением спутникового сигнала от препятствий и подстилающей поверхности, к получаемым рядам координат применялась сидерическая фильтрация [6]. Фильтрованный сигнал (по каждой компоненте) рассчитывался по следующей формуле:

put01.wmf

где X(t) – исходный сигнал; Y(t) – фильтрованный сигнал; оператор filt[ ] обозначает линейную полосовую фильтрацию сигнала; T – период повторения взаимного расположения навигационных спутников и приемной антенны, установленной в определяемом пункте.

Длительность сигнала после применения сидерической фильтрации уменьшается на 2·T. Период T был определен автокорреляцией рядов координат и составил для большинства станций 86154 или 86155 с.

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

Методика определения спектральной плотности мощности (СПМ) основана на алгоритме Уэлча и была максимально приближена к той, которая использовалась в работе [5]. Весь ряд координат по каждой компоненте каждой станции разбивался на сегменты длиной 43200 отсчетов. Соседние сегменты были сдвинуты друг относительно друга на 8640 отсчетов (на 20 %). Сегменты с более чем 120 пропущенными отсчетами отбрасывались, меньшие пропуски заполнялись линейной интерполяцией.

В каждом сегменте исключался линейный тренд и применялась оконная функция Ханна, после чего рассчитывалась СПМ. Полученные кривые СПМ усреднялись скользящим средним с шириной в одну октаву, т.е. для некоторого периода T усреднялись все значения СМП, попадающие в диапазон от T·2-½ до T·2½. В качестве итоговых значений были взяты восемь точек на октаву, т.е. периоды двух соседних точек отличались в 21/8 раз. Максимальный период был ограничен значением 10000 с.

Таким образом, для каждого сегмента были определены 95 значений СПМ, соответствующие средним значениям периодов от 2,83 до 10000 с. Общее количество полученных сегментов со всех станций составило 9059, таким образом, для каждого периода была получена выборка из 9059 значений СПМ. В этих выборках определялись значения СПМ, соответствующие уровням вероятности 5 %, 50 % и 95 %, аппроксимацией которых были построены следующие модели шума в рядах координат, полученных обработкой данных ГНСС-измерений в режиме реального времени методом PPP:

1. LRTGNM (Low Real-Time GNSS Noise Models) – модель низкого (5 %) шума.

2. MRTGNM (Median Real-Time GNSS Noise Models) – модель медианного (50 %) шума.

3. HRTGNM (High Real-Time GNSS Noise Models) – модель высокого (95 %) шума.

Каждая из указанных моделей построена отдельно для горизонтальных компонент (к названию добавляется нижний индекс NE, например, MRTGNMNE) и отдельно для вертикальной компоненты (нижний индекс UD, например MRTGNMUD).

Результаты исследования и их обсуждение

Построенные модели аналитически выражаются через набор опорных точек, для которых заданы значения периода и СПМ, приведенные в таблице.

Полученные модели низкого (LRTGNM), медианного (MRTGNM) и высокого (HRTGNM) шума

Период, с

Уровень шума, дБ относительно 1 м2/Гц

Горизонтальные компоненты

Вертикальная компонента

LRTGNMNE

MRTGNMNE

HRTGNMNE

LRTGNMUD

MRTGNMUD

HRTGNMUD

2,83

– 57

– 53,5

– 47,5

– 48,5

– 45,5

– 39

6

– 51,5

– 49

– 45

– 43,5

– 41,5

8

– 52

– 48,5

– 44,5

– 44

– 41

– 36,5

20

– 32

30

– 46

– 42,5

– 37,5

– 38

– 35,5

– 28,5

70

– 25

400

– 32,5

– 27,5

– 22

– 25

– 20

– 14,5

2000

– 17

– 11,5

– 4

3000

– 22

– 16

– 9

10000

– 13

– 4,5

2,5

– 5

4

11

Для произвольного периода T значение СПМ P может быть определено линейной интерполяцией при логарифмическом масштабе периодов:

put02.wmf

где Ti и Ti+1 – периоды соседних опорных точек (Тi ≤ Т ≤ Тi+1); Pi и Pi+1 – соответствующие им значения СПМ.

При необходимости модели шума могут также задаваться в скоростях и ускорениях:

put03.wmf put04.wmf

Полученные результаты моделирования шума рядов ГНСС-координат могут быть использованы для сравнения точности регистрации колебаний земной поверхности при помощи ГНСС и сейсмических приборов (широкополосными сейсмометрами и акселерометрами). На рис. 2 приведено сравнение построенных моделей шума с моделями сейсмического шума NHNM (New High Noise Model) и NLNM (New Low Noise Model) для широкополосных сейсмометров [13] и с моделями шума высокоточных акселерометров AHNM (accelerometric high noise model) и ALNM (accelerometric low noise model) [5].

В диапазоне периодов от 3 до 1000 с модели MRTGNM значительно выше моделей сейсмического шума NLNM и NHNM. При дальнейшем увеличении периода разница между ними уменьшается, т.е. чувствительность ГНСС-приемника и широкополосных сейсмометров становится сопоставимой. Модели шума для акселерометров ALNM и AHNM построены только для периодов до 150 с, они также значительно ниже моделей MRTGNM. Если провести экстраполяцию моделей ALNM и AHNM на более длинные периоды, то на периодах 300–2000 с эти модели были бы сопоставимы с моделями MRTGNM, а на периодах более 2000 с, по-видимому, чувствительность ГНСС-приемника значительно превосходит чувствительность высокоточных акселерометров.

pup2.tif

Рис. 2. Сравнение построенных моделей и моделей сейсмического шума: 1 и 2 – модели высокого (NHNM) и низкого (NLNM) сейсмического шума [13]; 3 и 4 – модели высокого (AHNM) и низкого (ALNM) шума высокоточных акселерометров [5]; 5 и 6 – построенные модели MRTGNMNE и MRTGNMUD

Таким образом, полученные результаты доказывают нецелесообразность замены «стандартных» сейсмических приборов (широкополосных сейсмических станций и высокоточных акселерометров) ГНСС-приемниками, вопреки ожиданиям некоторых исследователей [9]. Однако для получения более полного спектра сейсмических колебаний объединение записей сейсмического и ГНСС-оборудования [8] вполне оправдано, т.к. только спутниковые приемники в состоянии зарегистрировать сейсмические волны с периодами выше тысячи секунд и статические смещения.

Построенные модели шума рядов ГНСС-координат также позволяют оценить применимость ГНСС для регистрации сейсмических колебаний в заданном районе от определенного типа землетрясений. Для этого необходимо рассчитать реальные или теоретические (прогнозные) спектры сейсмических колебаний и сравнить их с моделями шума.

Исследования выполнены в рамках государственного задания Института тектоники и геофизики им. Ю.А. Косыгина ДВО РАН и при частичном финансировании РФФИ (проекты 16-05-00097a (Пупатенко В.В.), 17-55-53110 (Шестаков Н.В.), 17-55-50030 (Шестаков Н.В.).