Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
учебники / Клюшин Е.Б., Михелев Д.Ш., Барков Д.П. - Практикум по прикладной геодезии - 1993.pdf
Скачиваний:
379
Добавлен:
26.04.2015
Размер:
13.91 Mб
Скачать

Решив эту систему по формуле (285), получим уравнение осадки, аппроксимирующее результаты наблюдений;

St= 37,0448,09011+7,751512 - 0,552813.

(288)

Спрогнозировав осадки на

1989 г., на восьмой цикл наблюдений

и подставив в формулу (288)

значение /= 8 , можно

записать:

Sts= 37,04-48,0901 • 8 + 7,7515 • 64-0,5528-512= — 134,6 мм.

Аналогично можно вести прогноз и на девятый цикл. Далее прогнозировать нецелесообразно из-за возможных значительных оши­ бок аппроксимации.

§ 39. Исследование устойчивости реперов исходной геодезической основы

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

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

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

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

1$еР1>'И >/& ~.

(289)

где Scp — среднее значение смещений репера;

t — нормированный

множитель (критерий предельных ошибок), для подобных случаев принимаемый равным трем; ц— средняя квадратическая ошибка единицы веса; QSc — обратный вес определения величины Scp.

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

Для способов второй группы устойчивость оценивается по измене­ нию высотного положения реперов относительно средней отметки, вычисленной в текущем циклеj по отметкам п реперов исходной основы:

(290)

В этом случае задача сводится к выявлению изменения высотного положения неустойчивых реперов относительно средней отметки при условии, что менее 50% всех реперов сместились одновременно в одну сторону и на одну и ту же величину. При этом смещение относительно средней отметки можно считать значимым, если со­

блюдается условие,

аналогичное (289). Решение задачи рассмотрим

на примере способа, предложенного Ю. Е. Федосеевым.

Предположим,

что на площадке имеется п реперов и между

ними равноточно измерены превышения в нескольких циклах. По результатам измерений вычислены отметки всех реперов Ни (i— номер репера, j — номер цикла). При этом за исходный может быть принят любой репер, так как системой отсчета отметок в цикле в дальнейшем будет служить средняя отметка из отметок всех реперов, вычисленная

по формуле

(290). Отметки

реперов в этой новой системе

 

 

 

л

 

Н l j = Hj j — H cpj

= Hi j

= 1

(291)

 

 

 

 

 

n

Если в

течение

времени,

прошедшего между циклами j и j + 1,

все реперы оставались стабильными, то их отклонения от средней плоскости в пределах ошибок измерений сохранят свои значения.

Если же

к

моменту

проведения j +1

цикла изменилось положение

репера с

номером К

на величину Ак,

то

 

 

 

 

 

(292)

Н о

тогда

 

изменятся значение средней отметки на величину Ак/п

и

отметки

всех реперов относительно

средней плоскости:

Из уравнений (293) следует, что в результате изменения высоты одного из реперов остальные тоже изменили свое положение от-

А* ^ носительно средней отметки на величину---. Эта величина названа

п

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

Если

потеряли устойчивость два репера К и г, то реакция

системы

будет равна

Л= Д*+Д,

аотметки реперов К и г относительно средней отметки изменятся соответственно на величины

ДК- А?± А' ; Дг_ АЛ ± ^ ,

Следовательно, по изменению отметок реперов относительно средней отметки можно судить о количественных характеристиках устойчивости реперов. Таким образом, задача анализа сводится к выявлению изменения положения неустойчивых реперов относитель­ но средней отметки (средней плоскости) на фоне реакции всей

системы (рис. 6 6 ).

При анализе

длинных

рядов наблюдений кривые,

соответствующие

уклонениям

могут

настолько удалиться друг

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

относительно

средней отметки в смежных циклах,

т. е.

8 ,-,j+

, —H'ij.

(294)

Рис. 66. Схема графического изображения изменения отметки реперов от­ носительно средней плоскости

В случае изменения положения репера К на величину А* это изменение проявится и в разностях 5, т. е.

(295)

Вывод об изменении положения реперов возможно сделать лишь

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

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

Примем следующую классификацию смещений:

скачкообразные смещения Af, превосходящие доверительный ин­ тервал;

медленные смещения r|f, которые не превосходят для смежных циклов доверительный интервал, но постепенно накапливаются и мо­

гут

превзойти

его.

 

 

 

 

 

Величины

разностей

могут

также

содержать

случайные ошибки

г,, как следствие ошибок

измерений.

 

 

Ък j

В общем случае для сместившегося репера с номером К разность

можно

представить

в следующем

виде:

 

 

 

 

 

 

 

 

(296)

 

Второй

член в выражении

(296)

содержит

также собственные

смещения репера К. Чтобы их отделить от реакции системы, прибавим и вычтем величину --(Дкч]+ т\кtj + ткj). После приведения подобных

членов для 8К j в смежных циклах будем иметь

Накопление смещений к моменту проведения цикла будет выражаться суммой, т. е.

+ %.; + 4 j) + -

I

K,J+ T\K.J+ 4 , J)

(298)

п ;=

1

 

J

 

 

 

1

 

 

 

Так как т— случайные

величины,

матеметическое

ожидание которых

равно нулю, то

можно

считать,

что

 

В этом случае выражение (298) примет вид

X

X (АК., + ЛК.Л-Г^

X X (Ai,y+Л1,j)+

j =1

п 7=1

L

J=1*=1

1 т

~|

 

 

+- X (АК.>+ЛК.Л •

 

(30°)

У= 1

J

 

 

В формуле (300) первый член характеризует сумму собственных смещений репера К, а выражение в квадратных скобках — реакцию системы на скачки и медленные смещения всех реперов.

Для определения численных величин смещений репера К необ­ ходимо их выразить через вычисленные значения 5. Из формулы (300) с некоторым допущением можно вывести

j= 1

 

j= 1

 

Я- Ь « 1

 

,

т

л

т

 

 

+^ГТ .X

52.,+ ... + ^ Т

.X

5».г

(301)

 

Выражение (301) — рабочая

формула для вычисления смещений

любого

из

исследуемой системы реперов.

 

Рассмотрим последовательность вычислений на примере анализа устойчивости куста опорных реперов сети (рис. 67).

Пусть выполнено 12 циклов равноточных измерений превышений h между реперами (табл. 35). Вычисления выполняются в следующей

последовательности (табл. 36).

 

 

 

 

Т а б лица 35. Ведомость

наблюдаемых

 

превышений

 

 

 

 

Цикл

/ll _ 2 , ММ

hi -з, мм

А2-з>мм

 

1

+0,90

-0,34

-1,24

 

2

+0,82

-0,50

-1,32

 

3

+0,76

-0,58

-1,34

 

4

+0,62

-0,78

-1,40

 

5

+0,24

-0,86

-0,62

 

6

-0,40

-0,96

-0,56

 

7

-0,52

-1,12

-0,60

 

8

-0,68

-1,22

-0,54

 

9

-0,84

-1,28

-0,46

 

10

-1,00

-1,48

-0,48

 

11

-1,04

-1,98

-0,94

Рис. 67. Схема сети опорных

12

-1,08

-2,00

-0,92

реперов

1. Определяют отметки всех

реперов относительно одного из

них,

например,

первого

(графы 2 ,

3 , 4

).

2.

Находят

среднюю

отметку

HcpJ

(графа 5).

36. Ведомость вычисления изменений отметок реперов, мм

о

о

о

о

о

 

сч

 

rN

 

s

о

 

о

 

о

 

 

o'

о"

о"

о

о

о"

 

o'

 

s

 

 

 

5

 

+

 

 

 

+

I

 

+

 

 

 

 

 

 

 

>*

 

 

 

 

 

 

 

 

 

 

к

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

п

о

 

о\

(N(п (N

40

(N

 

г-

CN

н

 

 

к

о о о

(NО О О О о

 

о

ГО

п

 

о" о" о" О О

о"

o'

о" о*

о"

 

о

 

 

 

+ W 1

1

1 +

1

1

X

 

 

 

о

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

00

 

 

 

 

 

 

 

 

 

 

8.

 

 

 

 

 

 

 

оо Tfr

/—ч

в

 

 

 

4fr

 

(N

00

к

о

о о

оо

< 4 (N

 

О

о

о

о

ОО

Я

о

 

 

 

 

 

 

 

 

 

 

o'

о

о" о"о"

О

o'

о

о

о" о о

я

 

<D

 

I

I

 

 

1

1

1

1

+ 1 1

в

 

 

 

 

 

 

 

 

 

'w '

S

о

о

vr>

—>

<41 0 0

о

оо г-

ГО

00

о

ml оо

o'

о о

о

о"

o'

о

о

о1оо

о" о"

+

+

 

+

+ +I + + + + + + +

оо о о о о о о о

оо о о" о о" о о" о" I I 1 + 1 + 1 1 +

 

 

3

го

о\

 

г-

ON

00

40

 

о

г-

 

о"

о

ю

 

 

о*4

o'

о

o'

о*

о

 

I

I

 

 

I

I

 

I

I

 

 

о

г-

го ю

го 1Г>

ГО

г-

 

 

Г"

40

 

о

о

о

 

 

о"

о

о"

о

о

o'

о"

о

о

о

 

 

 

+

+

+

 

I

I

 

ON

—>

чо ю

г-

1Г>

«Л>

ГО

о

го

 

о

о

го TJ-

 

чо г-

оо

о

о"

о"

o'

о*4

o'

o'

о

o'

о"

1

1

1

+

+

+

+

+

+

+

+6Г0

+

90‘0+

0,05-

0,37-

0,45-

0,55-

0,63-

0,70-

0,83-

 

о

 

 

 

 

 

 

 

 

 

о

оо 00

NO

40

rN

(N

00

оо

го m

in

 

оо ON

*—•

(N

(N

 

о

o'

о

о"

о

о"

 

 

г

 

 

I

I

 

I

I

 

I

I

 

о

<N

40

<N

Tf

§

cs

oo

(N

о

ON

оо г-

40

<N

 

40

00

o^

o'

o'

o'

о*

o'

o'

о

o"

o'

 

+

+

 

 

 

 

 

I

I

 

 

 

 

 

о

8

 

о

о

о

 

 

 

 

 

о

о

о

 

 

 

 

о"

о*

 

о"

o'

o'

-ч

 

Л

 

н

(N(N

ас

о о

 

оо

 

 

+

 

аг

 

 

 

 

п

 

 

о

 

 

с

о

о

 

o'

о"

 

+

+

 

г-

г-

 

ач

ON

 

о

О*4

 

I

I

 

го

 

Л

о

о

XX

о"

о

0)

I

I

сок

 

 

л

 

 

а

 

 

с

1,01-

 

fit

 

S

 

 

а>

 

 

СО

 

 

а

оо

 

X

ON

 

 

л

 

 

us

 

 

ю

 

 

о

"Г*

 

н

 

о

8

8

о

о

Таблица

3. Вычисляют

отметки

реперов Н\ в системе высот

средней

отметки (графы

6 , 7, 8 ). Контролем вычислений служит

равенство

(в пределах ошибок округлений) суммы отметок

нулю

(графа 9).

4. Вычисляют

разности

р характеризующие

изменение

отметок

реперов H lj от цикла к циклу (графы 10, 11, 12).

Контролем

служит

равенство нулю суммы разностей по строкам

(графа

13). Для

наглядности по величинам

5itj строится график

(рис. 6 8 ).

 

 

Н,мм

Рис. 68. График изменения

отметок реперов относительно средней

плоскости

5.

Анализируют

численные значения разностей. Из

материалов

по оценке точности высотных измерений средняя квадратическая

ошибка определения превышения между реперами

равна

0,05 мм,

тогда предельная ошибка будет составлять 0,15 мм.

Если

рассмат­

ривать вычисленные разности, то в основном их значения не превосходят предельной ошибки измерений. Исключение составляют значения разностей между четвертым и пятым, а также между десятым и одиннадцатым циклами. Очевидно, что в эти периоды произошли скачкообразные изменения положения реперов: Рп 2 за время, прошедшее между четвертым и пятым циклами, опустился; опустился также Рп 3 за время между десятым и одиннадцатым циклами. Основанием для такого предположения служат в обоих случаях величины и знаки разностей, соответствующие выводам из формул (295), т. е. реакция системы в (п—1 ) раз (для данного примера в два раза) меньше величины смещения репера и проти­ воположна ему по знаку.

Дальнейшие вычисления проводятся по формуле (301) для опре­ деления смещений реперов. Поскольку эти вычисления для каждого

из реперов однотипны, рассмотрим

их на примере Рп 2. Так как

число реперов равно трем (см. рис. 6 8

), формула (301) будет иметь вид

Обозначив в формуле (302) каждый из ее членов через Ti j9 можно определить их численные значения (табл. 37). В графе 2 этой

таблицы даны

численные значения

собственных смещений

Рп 2,

а в графах 3

и 4— соответственно

влияние смещений Рп 7

и Рп 3

на численные значения смещений Рп 2. Вычисления следует выполнять с использованием значений Ьк} и 5itj из табл. 36, суммируя для каждого цикла все предыдущие значения и умножая эту сумму на

соответствующий коэффициент.

Например,

для цикла

5

и Рп 2

суммируются четыре значения и

Г2 5 будет

равно ~ (0

,0 0

0 , 0 1

0,013— 0,54)= —0,87 мм. В этом же цикле для Рп 7 и Рп 3 необходимо из значений 5 исключить реакцию системы на скачкообразное

смещение

Рп 2. Для этого

значение 52 5

в

графе

И

табл. 36 надо

разделить

на п 1

(в данном случае на 2

)

и

полученную величину

со своим

знаком

сложить

со значениями

 

б1 > 5

и

8 3 > 5

в графах 1 0

и 12 соответственно. Полученные величины, свободные от влияния

скачкообразного изменения

положения

Рп 2,

показаны

в

скобках.

Аналогичным образом выполнены вычисления для цикла

1 1

. Однако

здесь скачкообразно

изменилось положение Рп 3, поэтому,

учитывая

его влияние на другие, собственное его значение 6

в

этом цикле

должно

быть принято равным нулю.

 

 

 

 

 

В соответствии с формулой (302) искомое смещение Рп 2 в каждом

цикле

определяется

суммированием величин

Tt у

В табл. 37 эта

величина обозначена

знаком

Е (графа

5).

 

 

 

 

Полученные численные значения смещения Рп 2 вводятся в качестве поправки в его первоначальную отметку (графа 6 ). В данном примере отметка Рп 2 в начальном цикле была принята за нулевую. Для наглядности изменения отметки репера под влиянием его смещений изображаются на графике (рис. 69).

Т а б л и ц а 37.

Ведомость

вычисления

собственных

смещений реперов, мм

Цикл

Т и

Т и

T yj

I

я 2испр

1

2

3

4

5

6

1

0,00

0,00

0,00

0,00

0,00

2

0,00

+ 0,04

-0,04

0,00

0,00

3

-0,02

+ 0,06

-0,06

-0,02

-0,02

4

-0,06

+ 0Л2

-0,10

-0,04

-0,04

5

-0,87

+ 0,14

-0,12

-0,85

-0,85

6

-0,99

+ 0,18

-0,14

-0,95

-0,95

7

-1,02

+ 0,24

-0,16

-1,00

-0,94

8

-1,14

+ 0,28

-0,17

-1,04

-1,04

9

-1,26

+ 0,31

-0,18

-1,12

-1,12

10

-1,32

+ 0,38

-0,20

-1,14

-1,14

И

-1,35

+ 0,38

-0,20

-1,17

-1,17

12

-1,38

+ 0,40

-0,20

-1,18

-1,18

Соседние файлы в папке учебники