Ответ на вопрос
Кратко и по делу — принципы, алгоритмы, критичные источники ошибок и способы их учета.
1) Наблюдения и базовая математика
- Основные датчики: GNSS (GPS/GLONASS/Galileo), InSAR, вспомогательно VLBI/SLR/DORIS для привязки системы.
- Наблюдательные уравнения (псевдодальность и фазовая):
\[
P = \rho + c\,(dT - d\tau) + I + T + \varepsilon_P
\]
\[
\Phi = \frac{\rho}{\lambda} + N + \frac{c\,(dT - d\tau)}{\lambda} - \frac{I}{\lambda} + \frac{T}{\lambda} + \varepsilon_\Phi
\]
где \(\rho\) — геометрический векторный диапазон, \(c\) — скорость света, \(dT,d\tau\) — ошибки часов приёма/спутника, \(I,T\) — ионосферный и тропосферный задержки, \(N\) — целые неоднозначности, \(\lambda\) — длина волны, \(\varepsilon\) — шум/мультипуть.
2) Алгоритмы обработки GNSS для мм‑точности
- Схемы:
- Дифференциальная обработка (двойные разности, базели) — высокая точность для локальных сетей.
- PPP (Precise Point Positioning) — позиционирование на одиночной станции с точными эфемеридами и часами.
- Сетевое решение (network adjustment) с глобальной моделью ошибок.
- Оценка параметров: вектор положения, часы, тропосферные параметры (ZTD и градиенты), плавающие/целые неоднозначности; метод — взвешенный НМНК / Калмановский фильтр (реальное‑время) с возможностью последующей сглаживающей процедуры.
- Разрешение неоднозначностей: метод LAMBDA (integer least squares) — этапы: float‑оценка → декорреляция → поиск целых решений → валидация (ratio‑test, порог обычно \(\ge 3\)).
- Стохастическое моделирование: вес по углу возвышения, корреляции по времени; оценка компонент дисперсии (variance component estimation) для реалистичных погрешностей.
- Привязка к референсной системе: реализация ITRF через IGS решения; оценка EOP/масштаба/сдвигов станций.
3) InSAR и его комбинация с GNSS
- InSAR даёт относительную карту фазовых сдвигов; методы PS / SBAS / SqueeSAR для миллиметровой точности в линии видимости.
- Обработки: орбитальная ретракция, атмосферная коррекция (DEM, атмосферная модель или GNSS), фазовое развёртывание, фильтрация постоянных отражателей.
- Комбинация GNSS+InSAR даёт 3D‑векторы движений и коррекцию атмосферных и орбитальных ошибок InSAR.
4) Критичные источники неопределённости и их учёт (с типичными порядками величин)
- Фазовый шум измерения: \(\sim\! \) \(\) \(1\!-\!5\ \mathrm{mm}\) на одно измерение — учёт через весовые матрицы и фильтрацию.
- Неоднозначности \(N\): неверная фиксация даёт ошибки \(\sim\! \mathrm{cm}\) — надёжное разрешение LAMBDA + валидация.
- Орбиты и часы спутников: прецизионные эфемериды IGS ухватывают ошибки до порядка \(\sim\! 2\!-\!5\ \mathrm{cm}\); при PPP невязка орбит/часов трансформируется в позицию — использовать IGS ORB/CLOCK и, по возможности, поправки (UBX, CLK‑CORR).
- Тропосфера:
- Гидростатическая часть моделируется с точностью \(\sim\! 1\ \mathrm{mm}\) вертикально; влажная часть («ZWD») изменчива \(\sim\! 5\!-\!50\ \mathrm{cm}\) и даёт остатки вертикально \(\sim\! \mathrm{mm}\!-\!\mathrm{cm}\). Меры: оценивать ZTD как параметр (в каждый эпизод), применять продвинутые mapping functions (VMF1/GMF), моделирование градиентов.
- Ионосфера: при двойной частоте первый порядок устраняется (iono‑free LC), остаются higher‑order эффекты \(\ll 1\ \mathrm{mm}\) — при необходимости учитывать через модели/рефракции.
- Антенны (PCV/offset): необработанные фазовые центры дают см‑уровень ошибок; использовать калибровки IGS ANTEX и локальные измерения, корректно учитывать высоту/ориентацию.
- Мультипуть: может давать от \(\mathrm{mm}\) до \(\mathrm{dm}\) в зависимости от сайта; смягчение: правильная установка антенны, antena pattern, sidereal filter, сцепленные весы/робастные методы.
- Земные и атмосферные нагрузки:
- Твердая Земля (солидный прилив): вертикальный эффект до \(\sim\! 0.3\ \mathrm{m}\) — обязателен (IERS конвенции).
- Океаническое приливное нагружение: до \(\sim\! 40\ \mathrm{mm}\) локально — применять модели (FES/GOT).
- Атмосферное давление: вариации до \(\sim\! 10\ \mathrm{mm}\) — применять глобальные поле‑нагрузки (ECMWF/NCEP) или estimate/undo.
- Гидрологическая загрузка: сезонные \(\sim\! 10\!-\!30\ \mathrm{mm}\) — если требуется, моделировать/учитывать.
- Earth orientation parameters (EOP): ошибки EOP дают глобальные смещения мм‑уровня — использовать оперативные IERS/IGS EOP.
- Референсная сеть/связь антенн (tie vectors): невязки антенно‑вндерных измерений и локальных связей дают систематические смещения — контролировать через локальные измерения и комбинировать техники.
5) Учёт временных корреляций и оценка погрешностей времени‑серий
- Сигналы времени‑серий не белый шум: присутствуют flicker‑noise и random‑walk. Оценить PSD/коэфф. спектра/ML‑оценкой (maximum likelihood) и применять GLS вместо простого OLS для корректных доверительных интервалов.
- Моделировать и оценивать: линейный тренд (скорости), сезонные гармоники, скачки (co‑seismic), постсейсмические релаксации (экспоненциальные/логарифмические), трансляционные смещения.
- Валидация: независимые решения (сетевое vs PPP), бенчмаркинг по IGS/локальным точкам, синхронизация со сейсмологией/InSAR.
6) Практический рабочий процесс (рекомендуемый)
- Предобработка: синхронизация часов, флагирование/фильтрация наблюдений, применение ANTEX.
- Модельные поправки: IERS приливы, ocean load, pole tide, relativistic correction, Sagnac, solid Earth tides.
- Выбор стратегии: сетевое решение для малых‑средних сетей; PPP+ambiguity fixing для автономных станций.
- Шаг оценки: float‑решение → декорреляция и LAMBDA → фиксация → реоценка параметров и ковариаций.
- Постобработка: корректировки нагрузок, выравнивание к ITRF, построение временных рядов и оценка шума (ML/GLS), комбинирование с InSAR/SLR/VLBI для устранения длинноволновых ошибок.
- Валидация и QC: анализ редуцированных величин, остатков, повторные решения с разными моделями.
7) Итоговые рекомендации
- Для реальной мм‑точности: использовать точные IGS эфемериды/часы, калиброванные PCV, оценивать тропосферный ZTD и градиенты, применять LAMBDA‑фиксацию, моделировать все нагрузки по IERS и океаническим моделям, правильно учитывать стохастику времени‑серий (GLS/ML). Комбинировать GNSS и InSAR для полноты 3D‑вектора и для коррекции атмосферы/орбиты.
Если нужно — могу дать компактную схему входных данных, программного пакета (например GAMIT/GLOBK, Bernese, GAMIT, GipsyX, TEQC, GMTSAR, SNAP) и последовательность команд/настроек для реализации этого рабочего процесса.
Еще