Бесплатный автореферат и диссертация по наукам о земле на тему

Разработка метода определения поверхности скольжения оползня по данным геодезического мониторинга

ВАК РФ 25.00.32, Геодезия

Автореферат диссертации по теме "Разработка метода определения поверхности скольжения оползня по данным геодезического мониторинга"

005050268

На правах рукописи ф]/

Кузнецов Александр Игоревич

РАЗРАБОТКА МЕТОДА ОПРЕДЕЛЕНИЯ ПОВЕРХНОСТИ СКОЛЬЖЕНИЯ ОПОЛЗНЯ ПО ДАННЫМ ГЕОДЕЗИЧЕСКОГО МОНИТОРИНГА

Специальность 25.00.32 - Геодезия

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

6 МАР 2013

Москва-2013

005050268

Работа выполнена в ФГБОУ ВПО «Государственный университет по землеустройству».

Научный руководитель: Баранов Владимир Николаевич

доктор технических наук, профессор.

Официальные оппоненты: Клюшин Евгений Борисович

доктор технических наук, профессор Московского государственного университета геодезии и картографии;

Писаренко Владимир Кондратьевич

кандидат технических наук, доцент, зав. кафедрой Московской государственной академии коммунального хозяйства и строительства

Ведущая организация:

ФГУТТ Центральный научно-исследовательский институт геодезии, аэросъёмки и картографии им. Ф.Н. Красовского.

Защита состоится 28 марта 2013 года в 10.00 на заседании диссертационного совета Д 212.143.03 при Московском государственном университете геодезии и картографии (МИИГАиК) по адресу: 105064, г. Москва, Гороховский переулок, д. 4, зал заседаний Учёного совета.

С диссертацией можно ознакомиться в библиотеке МИИГАиК.

Автореферат диссертации разослан и размещён на сайте МИИГАиК www.miigaik.ru 26 февраля 2013 г.

Ученый секретарь ^ Климков

диссертационного совета «гР/о-и —— Юрий Михайлович

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы.

Развитие городов и сельских населённых пунктов сопровождается хозяйственным освоением новых территорий, включая площади, поражённые оползневыми процессами. Эффективное использование таких земель требует изучения геометрии и динамики оползней в рамках инженерных изысканий.

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

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

Объект и предмет исследования - объектом является оползневое тело, предметом - поверхность скольжения оползня (ПСО).

Цель работы - разработка метода определения высот ПСО по данным геодезического мониторинга.

Основные задачи исследования:

1) анализ проблемы на основе опубликованных материалов;

2) разработка теории метода определения высот ПСО по данным геодезического мониторинга;

3) разработка методики косвенного измерения высот ПСО.

Методы исследований: математический анализ, аналитическая геометрия, теория ошибок измерений, статистический анализ, натурный эксперимент.

Научная новизна работы. Разработан и внедрён метод определения высот ПСО по данным наземного геодезического мониторинга.

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

На защиту выносятся следующие научные результаты:

1) расчётный метод определения высот ПСО по результатам геодезического мониторинга земной поверхности: исходные постулаты, алгоритмы обработки данных, оценка точности входных и выходных данных;

2) методика измерения высот ПСО на основе геодезического мониторинга. Публикации н апробация работы. Результаты диссертационных исследований доложены на общероссийских и международных конференциях: СЕОРОЯМ (Москва, 2009), ГУЗ (Москва, 2009 [3], 2010 [5]), МИИГАиК (Москва, 2010), ПНИИИС (Москва, 2009 [2], 2010 [4], 2011 [9]).

По оползневой тематике опубликовано 9 статей, из них три [6-8] - в рецензируемых журналах ВАК. Результаты внедрены в производство ООО «Экостроймониторинг» [9] и использованы в учебном процессе на кафедре Геодезии и геоинформатики ГУЗ.

Объём и структура работы. Диссертация состоит из введения, пяти глав, заключения, приложений и библиографического списка. Общий объём диссертации составляет 168 страниц, в том числе 34 рисунка, 3 таблицы и 5 приложений (35 страниц). Библиографический список содержит 71 наименование.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Первая глава «Основання для изучения поверхности скольжения оползня—(ПСО) по ве»сторам смешений деформационных пунктов»

посвящена постановке задачи диссертационного исследования.

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

Предметом диссертационного исследования является геометрия ПСО -нижней границы между смещающимися и неподвижными породами. В работе приводится обзор возможностей существующих методов изучения ПСО и даны сведения об их точности. Установлено, что в обычной практике инженерных изысканий достаточная точность определения высот ПСО составляет 1 м.

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

На базе этих постулатов советскими учёными было предложено несколько алгоритмов построения ПСО по данным геодезического мониторинга: Г.Л. Фи-сенко (1956), И.О. Тихвинский (1977), Г.И. Лебедь (1980). Эти алгоритмы, используя приближённый графический путь решения, не удовлетворяли метрологическому требованию Единства измерений, т.к. не предусматривали полноценной оценки точности результатов (отсутствовала связь между точностью

входных и выходных данных). Этот недостаток стал препятствием к широкому распространению метода изучения ПСО по векторам смещений пунктов.

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

Во второй главе «Основы теории метода определения высот ПСО по векторам смещений пунктов» даны основные понятия и выведены формулы для оценки случайной погрешности метода.

При определении высот ПСО по векторам смещений пунктов реализуется тригонометрический принцип - превышение находят через угол наклона

1/^,,.И1 хорды к ~{к +1) и расстояние между точками (рисунок 1).

Пункт К

в точке к+1

Рисунок 1 - Вертикальный разрез оползня по линии К-(К+1)

Способ вычисления тангенса угла наклона хорды зависит от исходной модели линии скольжения. Основанием для выбора модели является тип оползня и

качество используемой геодезической информации. Если принять модель ПСО в виде дуги окружности, предполагая линейное изменение наклона ПСО (асеквентные оползни сдвига или низкая достоверность измерений), получим:

Л"Ч> — с .„„ ^»-нм „ "*-(»+!) — шп---= 0,

у„ соэса^, - а;) + у,н со.ч(огмм> - а™,)

. (1)

гсоБ (а,-О 2соз(ог,+1-<+|)

где _ дирекционный угол и горизонтальное проложение линии;

. _ Углы наклона векторов смещения в соответствующих точках; °Г. ам ~ дирекционные углы максимального наклона ПСО;

- угол наклона проекции вектора смещения К на створ К-(К+1)\ ~ Угол наклона проекции вектора смещения К+1 на створ {К+1)-К (рисунок 1).

Если считать линию скольжения параболой (оползни выдавливания), то:

г,«*.* _о ^нм:1»1'«.!^

-,)нм| :;--—г—

2 2

• (2)

соз(а4 -О со-а£,)

где 'м - уклоны векторов смещения.

Необходимость учёта от" возникает в связи с проектированием векторов

смещений на направление нивелирования. Параметр ат соотносится с

направлением падения мезорельефа оползня и определяется по горизонталям топографического плана.

Если в направлении смещения построить профиль по нескольким деформационным пунктам и знать высоту ПСО в одном из них Нто можно найти высоты других точек поверхности скольжения:

И, (3)

где нк - высота определяемой точки к на поверхности скольжения;

[й] - сумма превышений по ПСО от исходной точки до определяемой. Метод получения высот ПСО по результатам геодинамических измерений (векторам смещений пунктов) обладает следующей особенностью - поверх-

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

На основе выражения (2) в диссертации выведена формула для суммы

где к, г - порядковый номер и число пунктов в ходе; ДЛ"("_1Н1и1) - проекция линии (Ы)-(Л+1) на направление максимального наклона ПСО в точке к, С - максимальный уклон ПСО в точке к; кк - превышение на ПСО (рисунок 2: И, ЛЛ)-

Рисунок 2 - Ход динамического нивелирования в разрезе

На рисунке 2 проекции векторов смещений 1, 2 и 3 на линию хода параллельны соответствующим хордам 1-А, А-Б, Б-3, а слагаемые Ик суммы (4) соответствуют превышениям на поверхности скольжения /г,_, , /гь , . Найдём СКП функции (4):

превышений [к] в ходе динамического нивелирования:

А Ах;:

(4)

>

3

В диссертации выведены формулы, выражающие СКП максимального уклона ПСО т. „ через группы параметров: Л\- ,Ду ,Дг или /' ,е , а™. Во втором случае имеем:

т. +

е'*' } ) + ~) -- «Г))>'„- ^

: cos2 к-<) ' {6)

где ак,а" - дирекционные углы в точке к вектора смещения и направления

максимального наклона ПСО; а{к_,н*+1)- дирекционный угол линии (К-1)-(К+1);

Ь =1ап ^ и ек = ц> ак - уклон и тангенс дирекционного угла вектора смещения;

"Ч ''Пск~ СКП соответствующих параметров.

Применительно к продольным профилям, когда линия (К-1)-(К+1) совпадает с направлением смещения пунктов, = а, = ам , формула (5) примет вид:

»,,=Jt (7)

где S(i_,Mt.„ - расстояние между соответствующими пунктами профиля; 5,ср - средний шаг профиля в районе пункта К; горизонтальное проложение вектора смещения на пункте К; т^ _ СКП уклона вектора смещения пункта К;

m&x'mi* - СКП измеренных плановых и высотных смещений.

Для априорной оценки точности профильных сетей динамического нивелирования упростим формулу (7). Допуская равенство погрешностей одного превышения S?т, =mh -const , получим:

(8)

где ть - СКП одного превышения; г - число пунктов в ходе (профиле).

Для априорной оценки точности пяогцадных сетей динамического нивелирования (пункты произвольно распределены по территории оползня) из формулы (5), приняв а, = а", получим СКП превышения в ходе /и,,,:

ДХ-

~Ы>,ср-Гг, (9)

где , т1С[1 - СКП уклона вектора смещения к и среднее на объекте; ЛХ™ - среднее приращение абсцисс между смежными пунктами сети в оползневой системе координат (ось X" совпадает с направлением смещения оползня); остальные обозначения — см. формулу (4).

Для постановки динамического нивелирования необходимо определить плотность сети на разных участках оползня и продолжительность мониторинга.

Задавшись погрешностью одного превышения т„ и погрешностью уклона ПСО т.

'к *

из формулы (7) получим средний шаг продольного профиля в районе пункта К: с«.-"». . т^Тсозу,

где - допустимая погрешностью определения высот ПСО; £ы -

числовое значение длины оползня в метрах; Т - продолжительность мониторинга; - ожидаемая скорость смещения оползня; /•"'' - проектное число пунктов в профиле; остальные обозначения - см. формулы (6) и (7). Из формулы (10) можно найти расчётную продолжительность наблюдений Т:

'»ыТ,,^,,, ' (П)

где - средний шаг профиля; у р - предварительная оценка среднего наклона ПСО; г - количество пунктов в профиле; остальные обозначения -см. формулу (10).

В заключение второй главы проведён анализ влияния различных источников случайных погрешностей на результат динамического нивелирования и выработаны методические рекомендации для организации мониторинга:

• координаты X, К, 2 деформационных пунктов следует определять с точностью: тх г <5 мм, т2 <1 мм;

• зарегистрированное смещение оползня должно превышать 0,1 м;

• для уменьшения влияния погрешности т^ пункты оползневой сети следует

закладывать вдоль линий, ориентированных по направлению предполагаемого смещения оползня, и на удалении более 20 метров от его бортов;

• для уменьшения случайных погрешностей метода следует создавать деформационные сети с расстоянием по направлению смещения: в головной и подошвенной частях оползня ~ 5м, в зоне выноса Б?' <20м .

5—третьей_главе «Методические погрешности динамического

нивелирования» критически рассмотрены исходные постулаты метода.

Главное внимание уделено изменению кривизны ПСО. При использовании кусочно-круговой и кусочно-параболической моделей ПСО кривизна модели изменяется скачкообразно. Из-за несоответствия формы реальной ПСО и её теоретической модели, посчитанные по формуле (1) или (2) превышения будут содержать систематическую погрешность е4.

Для расчёта 0(, предлагается использовать высотные невязки трёхсторонних полигонов (на рисунке 3 - невязки , /2 , /", ):

\Л\ = Кг+Кг+К^ГГ, (12)

где й|_2Л-зЛ-|- превышения на ПСО между соответствующими точками хода;

/• дОП

ь - невязка трехстороннего полигона; -допустимая высотная невязка.

Направление смешения -►

'........ ' л ч

Я...................................... .

Трещина отрыва оползня .................................

Рисунок 3 - Схема вытянутого хода из пята пунктов, включающего три трехсторонних полигона

Преобладание в ходе невязок одного знака указывает на присутствие в сумме превышений монотонно изменяющейся систематической погрешности ®(",'Г ■ Если исключены грубые погрешности измерений, что обеспечивается на этапе предварительной обработки векторов смещений (глава 4), то недопустимая невязка сигнализирует о наличие в сумме превышений переменной систематической погрешности ЭЦ исходной теоретической модели ПСО.

Анализ данных мониторинга крупного оползня, расположенного в Коломенском музее-заповеднике г. Москвы, показывает, что размер систематической погрешности превышения существенно зависит от расположения пунктов на склоне и максимален в местах резких перегибов ПСО.

Ниже представлены формулы для расчёта монотонной и переменной составляющих систематической погрешности в ходе, выведенные для случая расположения пунктов на одном профиле (рисунок 3):

(13) (14)

где /к - невязки перекрывающихся трёхсторонних полигонов, входящих в профиль (на рис. 3 полигоны: 1-2-3-1, 2-3-4-2, 3-4-5-3); Л/= (¡Л | - ./Г" )>0 - часть

невязки /,, , превышающая в допуск /,,"""; N - общее количество невязок.

Сравнение количественных оценок монотонной и переменной составляющих систематической погрешности позволяет выбрать оптимальный способ обработки данных динамического нивелирования (глава 4).

Для уменьшения влияния систематической погрешности моделирования ПСО 01(1| рекомендуется:

- создавать деформационные оползневые сети профильной структуры;

- ограничивать максимальные расстояния между смежными деформационными пунктами 20 метрами;

-сгущать профили в местах резких перегибов ПСО, где |/А | >/;""';

- оптимизировать теоретические модели ПСО (см. главу 4).

Приведём итоги анализа других факторов, приводящих к методическим погрешностям динамического нивелирования:

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

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

• О форме ПСО мы судим по измеренным смещениям деформационных пунктов, отождествляемых с телом оползня. Однако в зарегистрированном смещении пункта присутствует составляющая, прямо не связанная с движением оползня.

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

ции по закладке деформационных оползневых пунктов, обеспечивающие их максимальную представительность, сохранность и удобство использования.

В четвёртой главе «Приёмы обработки данных динамического нивелирования» рассмотрены способы вычислительной обработки данных, уменьшающие влияние погрешностей на результаты динамического нивелирования.

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

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

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

Предобработка измеренных смещений пункта включает следующие этапы:

1. Корреляционный анализ плановых смещений Л\-, Ду, а также полных плановых Ал и высотных смещений &2.

Прямолинейный характер движения должен подтверждаться высокими

линейными корреляциями и (или) 'дг/ду. Низкий коэффициент корреля-

ции (меньше 0,5) может объясняться недостаточной точностью определения одного из параметров корреляции или её отсутствием. Если обеих корреляций не установлено, то пункт исключается из дальнейшей обработки.

2. Определение вероятнейших значений дирекционного угла и угла накчона вектора смещения с помощью уравнения регрессии.

Рассчитываем коэффициенты при уравнениях линейной регрессии е, I и их погрешности тс, т1 :

МДу) '

т, =-т—(Л-—; '"<

(т(Ду)У и-2

<г(Аг) [Г^г, <г(Аз)

,2

А?/Аг

(15)

(16)

где а, у - вероятнейшие значения дирекционного угла и угла наклона вектора смещения; '-Д(,Д1, - коэффициенты линейной корреляции между смещениями; о(Ах)М&у)М^),о(&*) - СКО ряда соответствующих смещений.

Найденные параметры используются для расчёта превышений и оценки точности динамического нивелирования.

3. Регрессионная оценка точности векторов смещений применяется при необходимости проведения независимой оценки точности геодезической информации, полученной из полевых измерений.

а) Регрессионная оценка плановых координат пункта проводится по отклонениям фактической траектории движения пункта от линейного тренда. Приняв = т,"" = (круговая область погрешностей - рисунок 4), получим:

|£(ГС05Г)3

(17)

где '"Гк - СКП регрессии (оценка плановых координат); у - отклонения ординат от линии тренда; г - румб вектора смещения; п - число циклов наблюдений.

А 1-

Линия тренда У= е + С/ \

Область \ погрешностей

Рисунок 4 - К расчёту СКП плановых координат по отклонениям от

линеинои регрессии

Регрессионная погрешность координат тх'г не идентична погрешности тх, полученной из уравнивания. Они соотносятся следующим образом:

. (18) где тх- СКП определения положения оползня, вызванная нестабильностью направления смещения пункта в горизонтальной плоскости.

б) Регрессионная оценка точности высотных измерений для каждого пункта проводится по данным его плановых и высотных смещений. Так как точность определения смещения в плане и по высоте может существенно различаться, то область погрешностей ограничена эллипсом с полуосями: - СКП планового смещения, и т'™ - СКП высоты (рисунок 5).

Высота /

Линия тренда Ж-! *Д? +С,

Рисунок 5 - К расчёту СКП высот по отклонениям от линейной регрессии

По аналогии с (17), рассчитаем характеристику разброса результатов относительно тренда, моделирующего траекторию движения пункта:

1

л — 2

(19)

где т'о' — СКП регрессии; г - отклонения по высоте от линии тренда (рисунок 5); V- наклон вектора смещения; и - количество циклов наблюдений. СКП высоты т^' получим из канонического уравнения эллипса:

11т1 ~(то' и«1")'

тЛзтО С°5 У

(20)

Далее по формуле (18), заменив индексы Х,У на 2, можно найти нестабильность направления смещения пункта в вертикальной плоскости т2.

Анализ экспериментальных данных, собранных автором в Коломенском музее-заповеднике, показывает, что «шумовая» компонента движения пункта в горизонтальной плоскости тхх в большинстве случаев пренебрежимо мала. В тоже время высотные статистические оценки те" заметно отличались от точности нивелирования тг. Несмотря на качественную закладку, сезонная «шумовая» компонента движения пунктов по вертикали составила тг => 1-2 мм.

Таким образом, опытные данные показывают, что в расчётах по формулам (7), (Ю), (11) погрешности смещений т^ и, особенно, т^ следует выражать через регрессионные оценки координат пунктов и т'7" , полнее

характеризующие точность определения смещений оползня.

Далее в работе излагаются два способа обработки профильных сетей, позволяющие уменьшить систематическую погрешность динамического нивелирования ©¡ф вызванную несовершенством кусочно-круговой и кусочно-параболической моделей ПСО.

Первый регрессионный способ предполагает выбор математической модели линии скольжения в процессе статистической обработки всего ряда данных. Общий вид модели находится при условии максимального приближения её к обрабатываемому ряду данных и соответствия модели типу оползня.

Задача решается в два этапа. Сначала находится статистическая зависимость уклона линии скольжения в точке от её местоположения на профиле . Потом уравнение регрессии интегрируется по аргументу 5 для расчёта высот отдельных точек. В итоге получаем массив высот Нк вспомогательных точек профиля, который можно представить в виде таблицы или графика.

Регрессионный способ можно рекомендовать для моделирования линий скольжения типичных оползней выдавливания, имеющих правильную форму (парабола, клотоида и т.п.) при большой плотности оползневой сети. У таких оползней преобладают монотонные погрешности: > 0^]'.

Второй сплайн-способ вычисления высот следует применять в том случае, если в невязках превалируют переменные систематические погрешности ©¡^ > (сложная ПСО с резкими перегибами). Тогда для интерполяции уклонов целесообразно применение кубических сплайнов ¿с„„„„, = /(5), заданных на каждом интервале области определения отдельной формулой. Графически такой сплайн представляет собой кривую с плавно изменяющейся кривизной. Для вычисления высот вспомогательных точек профиля полученный сплайн икгегрируюг по аргументу 5.

Моделирование ПСО с помощью сплайнов является наиболее универсальным способом обработки данных динамического нивелирования.

В пятой главе «Разработка методики динамического нивелирования» создаётся алгоритм динамического нивелирования с достаточной для практики точностью определения высот ПСО и проводится апробация методики.

Структурная схема методики представлена в таблице 1.

Приведём основные положения методики динамического нивелирования, согласованные с ГОСТ Р 8.563-2009 ГСИ. Методики (методы) измерений.

• Область применения - изучение высот ПСО глубоких (до 50 метров) одноярусных оползней сдвига и выдавливания, находящихся в стадии подготовки основных смещений или повторных смещений (скорость 0,1-1 метр в год), путём продолжительного (1 год и более) геодезического мониторинга.

• СКП определения высот ПСО:

- до 1 метра (при пересечении профилями только верхних границ оползня);

- до 0,5 метра (профили пересекают нижние и верхние границы оползня).

• Геодезическая сеть - в процессе проектирования предпочтение следует отдавать профильной структуре сети с преобладанием продольных створов.

1. Подготовка измерений:

сбор информации об объекте;

проектирование сети и разработка программы измерений (8-11); закладка геодезической сети.

2. Выполнение измерений:

- выполнение ряда циклов полевых измерений с соблюдением допусков;

- определение координат и высот пунктов деформационной сета, расчёт смещений;

- оценка точности координат и высот; оценка точности смещений (если цикла два);'

- измерение по плану направлений максимального наклона ПСО.

3. Предварительная статистическая обработка (только для рядов смещений):

- контроль векторов смещений на основе корреляционного анализа;

- определение усреднённых угловых характеристик векторов смещений и их погрешностей по формулам (15, 16);

- регрессионная оценка точности положения пунктов (17-20) (при необходимости).

4. Оценка точности динамического нивелирования (внутренняя):

- расчёт случайной погрешности определения высот ПСО (5-7);

- расчет монотонной и переменной составляющих систематической погрешности определения высот ПСО (13-14);

5. Определение высот ПСО:

Профильная сеть

- расчёт превышений:

а} юсочно-круговой или кусочно-параболический способ: -расчёт превышений через наклоны илиуклоны векторов смещений смежных пунктов (I, 2);

- последовательное суммирование превышений (4);

б) сплайн- способ•

- интерполяция уклонов в профиле к}'-бическими сплайнами / я = ,'

- интегрирование функции по аргументу Б;

в) регрессионный способ:

- аппроксимация уклонов в профиле уравнением регрессии / = f(S); ■ интегрирование функции по аргументу 5;

- высотная привязка (3)\

- переход к ПСО интерполированием высот между профилями.

Площадная сеть

- расчёт превышений:

а) триангуляционный способ:

- разбиение сети на треугольники:

- расчёт превышений по сторонам треугольников через уклоны векторов смещений смежных пунктов (2); -расчёт весов превышений: Р~1/т;';

- уравнивание превышений;

б) способ квадратов ■

- построение планов продольного и поперечного уклонов ПСО;

- построение на планах сетки квадратов с шагом, соответствующим расстояниям л!ежду пунктами;

- расчёт превышений между узлами сетки квадратов (2);

- расчёт весов превышений: Р~ 1/т,';

- уравнивание превышений;

- высотная привязка сети;

- построение ПСО

6. Уравнивание высот ПСО

(при наличии избыточных данных: инклинометрии, тензометрии и т.п.)

7. Окончательная оценка точности (внешняя) модели ПСО:

- количественная (по результатам сравнения с данными прямых измерений высот

ПСО);

- качественная (по результатам экспертной оценки).

Количество продольных профилей определяется из расчёта один профиль на 100-150 м ширины оползня. Продольные профили не рекомендуется размешать ближе 20 м к бортам оползня. Расстояния между пунктами продольных профилей: в головной части оползня и зоне выдавливания - 5-10 м, в зоне выноса - до 20 м. В продольные профили желательно включать инклинометрические скважины.

Опорные пункты располагаются по периметру оползневого массива. Требования к закладке опорных и деформационных пунктов определены в Своде правил по инженерно-геодезическим изысканиям для строительства.

• Периодичность и продолжительность мониторинга - зависит от скорости движения оползня. Рекомендуемая продолжительность не меньше 1 года, при этом за время мониторинга зарегистрированное смещение должно превысить 0,1 м. Обшее количество циклов должно быть не меньше четырёх.

• Точность измерений - при медленных смещениях (около 0,1 м в год) должна обеспечивать определение плановых координат с погрешностью до 5 мм, высот - до 1 мм.

Направление максимального наклона ПСО, соответствующее направлению падению мезорельефа оползня, определяется по топографическому плану с максимально возможной точностью.

• Алгоритм обработки данных - зависит от плотности пунктов и конфигурации деформационной сети. В основных этапах он представлен в таблице 1 (графы 3-7). Процедура включает априорную оценку точности динамического нивелирования (при необходимости), контроль измеренных смещений, расчёт высот ПСО, уравнивание высот (при наличии избыточных данных), контроль и апостериорную оценку точности результатов.

В конце главы проводится апробация методики динамического нивелирования в части обработки результатов мониторинга. Объектом исследования послужил оползень выдавливания, расположенный в районе Хорошёвского спрямления реки Москвы. Деформационная сеть на объекте включает 6 продольных профилей, из которых 4 пересекают тело активного оползня.

Средний шаг профиля 13 м, полное смещение оползня за время наблюдений 0,1 м, число контрольных инклинометрических скважин - 5 шт.

На рисунке 6 представлен графический результат динамического нивелирования - план ПСО с горизонталями. СКП определения высот ПСО, рассчитанная по формуле (7), составила 1,6 м.

Контроль динамического нивелирования проводился по отклонениям найденных высот от данных скважинной инклинометрии. Расхождения в высотах точек, полученных двумя методами, не превысили 2 метров, что лежит в пределах допустимого расхождения, рассчитанного с использованием формулы (7).

Профиль 5

Рисунок 6 - План ПСО с горизонталями

ЗАКЛЮЧЕНИЕ

Ниже приведены основные итоги диссертационного исследования.

1. Проведённый анализ состояния дел при инженерных изысканиях опасных склоновых процессов выявил актуальность разработки метода исследования геометрии ПСО на основе геодезического мониторинга.

2. Показано, что главным препятствием на пути изучения ПСО по векторам смещений деформационных пунктов является недостаточная теоретическая проработка метода и, как следствие, отсутствие аттестованных методик выполнения измерений.

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

4. Даны методические рекомендации для уменьшения отдельных погрешностей метода: СКП определения направления максимального наклона, погрешности теоретической модели ПСО, погрешности за счёт неустойчивости деформационных пунктов относительно оползневого тела и других погрешностей. Предложена оптимальная структура оползневой сети - профильная с преобладанием линий, вытянутых по направлению смещения.

5. Разработана методика расчёта предельных расстояний между смежными пунктами в разных частях оползня и установлен критерий для дополнительного сгущения сети. Выведена формула для определения продолжительности мониторинга.

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

7. Разработаны способы моделирования ПСО (регрессионный, сплайн-способ), позволяющие уменьшить влияние непостоянства кривизны ПСО на результат динамического нивелирования. Даны рекомендации по их применению.

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

9. Методика успешно апробирована на оползнях выдавливания в г. Москве (Коломенское, Хорошёво).

10. Результаты проведённых исследований позволяют рекомендовать предлагаемый геодезический метод изучения геометрии ПСО для широкого применения при стационарном изучении блоковых оползней.

ПУБЛИКАЦИИ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ

1. Применение ГНСС при геодезических наблюдениях оползневого склона реки Москвы в районе Чертановского коллектора; монография «Информационное обеспечение кадастров пространственными данными» / А.И. Кузнецов. -М.:ГУЗ,2006 -с. 169-176.

2. Деформационные пункт: альтернатива традициям / АЛ. Кузнецов, С.А. Моисеенко //Журнал«Инженерныеизыскания»2/2009 -с.38-43.

3. Моделирование оползневых процессов по данным геодезических наблюдений на территории Коломенского парка г. Москвы / АЛ. Кузнецов // Материалы международной научно-практической конференции, посвященной 230-летию со дня основания ГУЗ.- М.: ГУЗ,2009 -с. 160-167.

4. Изучение поверхности скольжения оползня традиционными геодезическими методами / А.И. Кузнецов, С.А. Моисеенко, В.А. Волков // Материалы шестой общероссийской конференции изыскательских организаций - 2010,- М.: ПНИИИС, 2011 -с. 238-240.

5. Определение поверхности скольжения оползня по результатам геодезического мониторинга / АЛ. Кузнецов // Материалы шестой международной научно-практической конференции молодых учёных и специалистов - 2010, -М.: ГУЗ, 2011 - с. 156-161.

6. Опыт использования данных геодезического мониторинга для построения поверхности скольжения оползня. / АЛ. Кузнецов, С.А. Моисеенко, В.А. Волков // Журнал «Инженерные изыскания» 2/2011 -с. 56-59.

7. К вопросу организации геодезического мониторинга оползневых процессов / АЛ. Кузнецов //Журнал «Инженерные изыскания» 10/2011 - с. 30-33.

8. Оценка точности результатов геодезического мониторинга оползневых процессов с использованием регрессионного анализа / А.И. Кузнецов // Журнал «Геодезия и картография» 10/2011-е. 8-13.

9. К вопросу о структуре оползня на Карамышевской набережной р. Москвы / В.А. Волков, A.B. Тихонов, A.B. Калинина, С.М. Аммосов, АИ. Кузнецов // Материалы седьмой общероссийской конференции изыскательских организаций. - М: ПНИИИС, 2011.-е. 21-24.

Редакционно-издательский отдел ГУЗ Сдано в производство 15.02.2013. Подписано в печать 13.02.2013. Формат 60х 84'/16. Объем 1,0 п.л. Бумага офсетная. Тираж 130 экз. Заказ №

Отдел оперативной полиграфии ГУЗ 105064, г. Москва, ул. Казакова, 15