
книги из ГПНТБ / Снежные лавины (прогноз и защита) [сб. ст
.].pdfП р и л о ж е н и е
Результаты промежуточных вычислений по различным программам при опознавании лавиноопасных ситуаций
( по данным табл. 3 ;
* |
П D |
o n a u |
u a . |
У |
|
I |
П |
in |
1У |
У1 |
|
"Альтерна |
"Образ-3" "Поиок-А" |
"Поиск-Б" |
"Поиск-В |
"Поиок" |
|
тива" |
|
|
|
|
|
I 13,2
210,11
36,21
43,92
512,89
65,10
7 6,24
8- 0,63
91,45
10 |
0,25 |
I I |
1,72 |
12- 1,12
131,76
140,04
15- 0,46
16 0,64 Г ? 0,29
183,35
190,72
20- 0,63
210,73
223,60
234,30
242,50
250,18
26- 4,22
27- 0,67
28- 0,99
29- 1,48
301,93
10/0 |
0,528 |
0,637 |
0,999 |
0,992 |
|
0 /0 |
0,594 |
0,224 |
0,993 |
0,908 |
|
20/8 |
0,63 |
0,735 |
0,877 |
0,727 |
|
150/5 |
0,745 |
0,586 |
0,616 |
0,521 |
|
10/8 |
0,420 |
0,885 |
1,000 |
1,000 |
|
20/0 |
0,532 |
0,460 |
0,831 |
0,714 |
|
170/8 |
1,000 |
1,000 |
0,732 |
0,688 |
|
0/33 |
0,222 |
0,153 |
Ю“ 5 |
0,083 |
|
10/33 |
0,465 |
0,444 |
0,226 |
0 . 4J 5 |
|
0/48 |
ЗЛО-4 |
0,305 |
I Q - Z 0 |
I 0"6 |
|
ю-14 |
|||||
10/42 |
0,082 |
0,115 |
Ю-4 |
||
10/50 |
0,282 |
0,202 |
ю-15 |
Ю "4 |
|
10/50 |
0,295 |
0,194 |
Ю -20 |
Ю-4 |
|
0/8 |
0,507 |
0,472 |
0,279 |
0,360 |
|
0/8 |
0,234 |
0,027 |
Ю-7 |
0,080 |
|
20/58 |
0,035 |
0,161 |
10-5 |
0,015 |
|
0/133 |
0,163 |
0,200 |
Ю -7 |
0,003 |
|
130/32 |
0 . 4СГ7 |
0,452 |
0,272 |
0,304 |
|
10/100 |
0,110 |
0,242 |
0,062 |
0,203 |
|
0/133 |
ЗЛО-4 |
0,497 |
Ю -20 |
Ю -6 |
|
0/42 |
0,430 |
0,390 |
0,084 |
0,27 |
|
40/0 |
0,460 |
0,972 |
0,234 |
0,435 |
|
0/18 |
0,395 |
0,487 |
0,065 |
0,187 |
|
0/8 |
0,425 |
0,354 |
0,016 |
0,159 |
|
40/17 |
0,622 |
0,427 |
0,672 |
0,401 |
|
0/ Г 7 |
0,630 |
Ю -20 |
ю-7 |
Ю -20 |
|
0/ I I 7 |
0,039 |
0,104 |
i0 -20 |
Ю “ 5 |
|
40/17 |
0,182 |
0,172 |
I 0~i5 |
Ю " 4 |
|
0/42 |
0,467 |
0,353 |
0,004 |
0,116 |
|
0/25 |
0,152 |
0,363 |
0,170 |
0,282 |
20
Пояснения. 1) во второй колонке "Приложения"шро1раыма I )
приведены значения линейной решающей функции, вычисленной по программе "Альтернатива". Пороговое значение этой функции
=й,Ь2. 2) В третьей колонке ^программа Щ приведены соотно шения "гол осов ", вычисленные по программе "О браз-3". Число над чертой обозначает количество голосов за то , что ситуация
является лавиноопасной, под чертой - не лавиноопасной. 3 ; В ко лонках 4-7 ( программы Ш -Ш приведены значения адекватности, вычисленные по различным модификациям программы "Поиск". В программах Ши ьУ, предназначенных для ручного счета,корреля ционные связи между определяющими параметрами игнорируются, в
программах У и У1 - |
учитываются. Близость в программе Шопре |
||||||
деляется как |
1 / |
5 02 |
, в программе |
1У - как |
ъ х р ( - р о2) , |
||
в программе У - |
как |
ехр (~ Р 2) |
и в программе |
У1 - |
как |
||
e x p (-j3 n2) |
.Здесь j ° 0 |
- расстояние, |
вычисленное |
без |
учета корреляционных связей между определяющими параметрами;
Р - расстояние,вычисленное |
с учетом этих связей и ft - об |
|||||||
щее чиоло |
определяющих параметров. |
|
|
|||||
|
|
|
|
|
Литература |
|
|
|
Ф л я й г |
В. |
Внимание, |
лавины ! ИЛ,I960. |
|
|
|||
Л а в и н о о п а с н ы е |
|
р а й о н ы |
С о в е т с к о г о |
|||||
С о ю з а . |
Под ред. |
Г.К.Тушинского. |
Изд-во МГУ,1970. |
|
||||
Р у к о в о д с т в о |
|
п о |
с н е г о л а в и н н ы м |
р а |
||||
б о т |
а м . М ., |
Гидрометеоиздат,1965. |
|
|
||||
А к к у р а т о в |
В . Н . |
Прогноз наступления лавинной |
||||||
опасности по величине метелевого переноса и температур |
||||||||
ного сжатия снега. Сб. "Вопросы использования снега и |
|
|||||||
борьбы со снежными заносами и лавинами". М., изд. АН СССР, |
||||||||
1956. |
|
|
|
|
|
|
|
|
Б е л е н ь к и й Б. М. |
|
Климатические |
условия района |
|
||||
карьеров комбината "Апатит" имени С.М.Кирова. Сб."Снег |
|
|||||||
и лавины Хибин".М ., изд. МГУ,1967, ротапринт. |
|
|||||||
М и л л е р |
Р . |
Л. |
, К а н |
Д. С. - Статистический анализ |
|
вгеологических науках. М ., "Мир", 1965.
Бо н г а р д М. М. Проблема узнавания. М., " Наука", 1967.
2 -1159 |
21 |
|
В.Ф.Гракович
ОБ ИСПОЛЬЗОВАНИИ МЕТОДА РАСПОЗНАВАНИЯ ОБРАЗОВ ДЛЯ ОЦЕНКИ ЛАВИННОЙ СИТУАЦИИ
ПО КОМПЛЕКСУ МЕТЕОРОЛОГИЧЕСКОЙ ИНФОРМАЦИИ
В научной практике возникает 1фуг задач, связанных с распознаванием и предсказанием лавинных ситуаций. Вероятность прогноза схода снежной лавины зависит от большого числа инфор мационных параметров, а также от характера их комбинаций. Дан ная работа - лишь попытка выяснить корреляционные зависимости времени схода снежных лавин от изменений метеорологических факторов, которые приводят к срыву лавины. Анализируется воз можность применения метода распознавания образов для прогноза таких природных явлений,как сход снежных лавин.
Причина срыва снежных лавин - комплекс условий,при кото рых формируется снежный покров, а также диагенез снежной тол щи, изменяющий внутренние связи. Разнообразие географических условий сильно усложняет работу исследователя. Поэтому многие авторы /Аккуратов,I960; Тушинский,1962/ подчеркивают мысль о непригодности универсальных стандартных критериев для определе ния времени схода лавины и невозможности выработать единый ме тод прогнозирования лавинной опасности. Для каждого района при ходится искать свои способы решения этой проблемы. Большинство лавин, опасных для жизни и имущества людей, сходит во время сильных метелей или сразу после них. В тех случаях, когда схо дят лавины сублимационного диафтореза,необходимо обратить внима ние на анализ условий погоды,предшествующих событию. Понимание основ физических процессов,происходящих внутри снежной толщи, необходимо для правильной интерпретации натурных наблюдений и решения комплекса проблем,связанных с предсказанием времени падения лавин с долгопериодными процессами формирования. В этом случае вывод об устойчивости различных слоев снега на склоне сформулировать достаточно сложно, так как они подвергаются различным метеорологическим воздействиям и внутренним струк турным изменениям. Проблема прогноза схода снежных лавин под-
22
разделяется |
на две части: I) |
распознавание |
лавинной |
ситуации |
в отличие от |
нелавинной и 2) |
непосредственно |
прогноз |
с задан |
ной заблаговременностью.
Известны различные подходы к решению задачи прогнозиро вания снежных лавин. Статистический анализ факторов - о д а из методов оценки лавинной опасности. Метод применим, в первую очередь, к лавинам из метелсвого снега, которые сходят во время метели или сразу после ее окончания. Для определенных районов такие прогнозы дают хорошие результаты /Аккуратов, 1360/. Опасности от лавин сублимационного диафтореза предви деть чрезвычайно трудно, так как оценить каждый фактор, дейст вующий на сход лавины, практически невозможно без создания мощного математического аппарата и обработки статистического материала на ЭВМ. Проводимый метод прогнозирования по комплек су нескольких факторов /Справочник по прогнозированию,1964/мо жет оказать помощь при оценке лавинной опасности, однако,уве ренного ответа на вопрос о возможности схода лавины он не дает.
Совершенно ясно, что существующие методы должны быть точными
иобъективными.
Визложенных исследованиях вопроса схода снежных лавин
использовались комплексы метеорологических станций Торс кол
и Эльбрусской географического (факультета МГУ.
Для упрощения анализа синоптические данные (скорость ветра,
количество осадков, температура воздуха}образующие многомерный случайный вектор Е .подвергались предварительному сглажи ванию. Для выявления совместных зависимостей (осадки, ско рость ветра, температура воздуха) составлена программа "Графи ка ", где с помощью ЭВМ “0дра-1013" строятся сглаженные графи
ки, определяющие синоптическую ситуацию на каждый день (р и сЛ ). В дальнейшем графики метеорологических данных анализирова
лись и по ним строились функции распределения относительно времени схода лавин. Анализ графиков за IS67-I972 г г . показал,
что влияние температуры на возникновение лавины более много сторонне, чем любого другого фактора.Это подтверждается иссле дованием физических процессов,происходящих в снежной толше
/Тушинский,1962/. При теплой погоде,тогда температура близка к нулю, неустойчивость снега сильно увеличивается. Низкие тем
пературы |
с сильным ветром |
приводят к возникновению снежных |
"д осок ". |
Анализ количества осадков позволяет установить зави |
|
симость |
между количеством |
выпавшего снега, интенсивностью его |
23
г-' (s
Р и с.1 , сглаженный график ыетеоданных (апрель 1968)
(Г- |
количество |
осадков; V - скорость ветра; |
Т °- |
температура |
воздуха. |
выпадения и образованием снежных пластов.имеющих потенциаль ную возможность создать лавину. Оказывается, что существует некоторая критическая интенсивность выпадения -снега,вызываю- щая сход лавины.
Задачу распознавания лавинной ситуации оудем решать методом.предложенным М.М.Бонгардом /1 9 6 7 /. Она формулируется следующим образом. Имеется некоторое начальное множество объек тов , каждый из которых описывается конечным наоором извест ных параметров. Предполагается, что существует решающее пра вило , которое всегда разбивает множество объектов на непересекающиеся подмножества 1классы;„ В случае разбивания на два класса математическая часть решения задачи значительно упрощается.хотя идея остается такой же.
24
Разбивание множества объектов порождает и множество их
описаний.Следовательно, нужно построить систему, которая по описанию произвольного объекта устанавливала бы его принадлеж ность к соответствующему классу. При этом считаем, что сама "узнающая" система вырабатывает необходимое правило классифи кации объектов на классы в процессе обучения."Узнающая" сис тема с помощью стандартных математических операций,применя емых к объекту, и заложенной в программе логики определенный объект будет отнесен ЭВМ к конкретному классу.
В силу специфики работы ЭВМ возникает необходимость квантования параметров,осуществляемого на основании опыта ис следования и эвристического в том .смысле, что в его основе ле жат те или иные дополнительные интуитивные предположения,среди
которых наиболее распространенно предположение о "компактно сти" /Цыпкин,1970/. Вопросы, связанные с исследованием свойств параметров и их квантованием, в данной работе не рассматри ваются.
После представления параметров в удобном для машины ви
де каждый из а параметров ( , . . . <2^ ) будет являться набором двоичных чисел. Вместо указанных параметров
получаются новые параметры I 6 1. , i 2. , . . . 6ГП ) . То, что некоторые наборы булевских переменных изображают один фактиче ский параметр, учитывается при их кодировке. Таким образом,
каждый объект из пространства представляет собой набор двоич ных чиоел.или булевский вектор. Задача сводится к определению логических функций в качестве искомых правил для построения таксона. Булевские векторы.являющиеся отображением объектов,
на которых будет проводиться |
обучение,составляют матрицу Е • |
|
I , |
, 2 |
/ |
е ; |
£ |
|
i ; |
с |
|
prop |
го |
с и |
Необходимо найти характеристические функции X - .равные единице на Aj и нулю на Ag, где А( , к г - подмножества фа зового пространства,соответствующие объектам
*i { О ф } м .
Всего логических функций от S |
переменных существует |
25
2r 2SC k |
и перебрать все функции немыслимо, так как при не |
||
больших |
S получаются астрономические цифры. |
Но любую логи |
|
ческую функцию можно представить через другие |
логические^ |
^ |
|
функции. Тогда число функций X ^ будет уменьшено до 2 |
С к . |
В качестве такой логической функции выберем конъюнкцию. Эту функцию будем применять к любому сочетанию трех векторов из матрицы Е . Перебрав все конъюнкции трех векторов, получим совокупность трех признаков. Теперь остается выбрать из них самые надежные.Это можно сделать так: выбрав определенное чис ло признаков с максимальными "коэффициентами надежности" и задав "порог" lv , выбрать те признаки, для которых превыша ется "порог" К .Отобранные признаки и будут таксонами, с помощью которых можно отнести любой объект к тому или иному классу. На этом заканчивается один из этапов процесса обуче ния. Может случиться, что по каким-либо причинам забравован нужный признак. Тогда возьмем объекты,мало покрытые признака ми, и проводим доучивание, уменьшив "порог" 1C .Выберем до полнительные таксоны и включим их в число уже отобранных. Пос ле того , как весь класс объектов покрыт достаточно плотно так сонами, можно производить экзамен, т . е . вводить данные об очередном объекте для отнесения его к соответствующему классу.
Рассмотрим реализацию этого алгоритма. Для распознавания
лавинной ситуации |
составили |
программу |
"Лавина- I " |
для ЭВМ |
|
"Ыинск-22". Объектами класса |
Д |
были сошедшие |
лавины,за |
||
регистрированные |
по времени на данном |
склоне,класса А - |
нелавинная ситуация в момент назревания лавинном ситуации,
наждыи период характеризовался набором значений метеоданных:
<о -сумма осадков; Y - скорость |
ветра; Т |
-температура |
воздуха; Т -продолжительность |
периода. |
На основании вве |
денных в машину описании нескольких лавинных и нелавинных ситуа ций программа должна построить таксоны, с помощью которых она будет различать описание новых событий.
Поскольку знания метеорологических условии недостаточно для предсказания схода лавины, то необходимо учесть всю предисторию создания снежной толши, накапливаемой в лавиносборе.
Для этого , описание любого объекта (множество А ) формирова лось как совокупность описании всех слоев снежной толши. Цен ность такого метода заключается в том, что характеристика сне говой обстановки на каждый, указанный заранее день задается описанием всех слоев снега, которые существовали до этого дня
26
после схода последней лавины, а не только последнего снего пада. .Условно были выделены следующие подсобытия, или слои: непогода,хорошая погода,длинный период, последние два дня перед сходом лавины.Смена погодных условий, влияющих на раз витие снежной толщи и возникновение лавин показана на р и с.2.
Рис.2 Схема погодных условий
Схема описания события для ввода в |
ЭВМ приведена в табл .i , |
||||||||
цийры в ней |
указывают разрядность параметра. |
|
|
||||||
|
|
|
|
|
|
Т а б л и ц а |
1 |
||
нодсооытия |
|
репогода |
..орошая |
погода |
долгимI последние |
||||
количество |
|
1 |
|
1 |
|
период!2 |
дня |
|
|
|
|
|
|
|
|
||||
двоичных |
4 |
12 3 , 2 |
I |
i |
2 |
4 2 3 .3 |
2 3 |
1 I |
|
разрядов |
з : |
||||||||
Параметры |
(Г V. ,т.°;t |
V. т .с \Х. |
ь; |
|
d |
||||
X Vi Н t. <1 \ |
|||||||||
|
е . Ч Ч V |
1 |
L 1 i |
27
Здесь : -количество осадков (гумма); \ -скорость ветра (средняя); Т 0 - средняя температура воздуха за период;
-продолжительность периода; Jс -интенсивность снегопада в последние 2 дня; 2 -признак лавины (сходила или не сходила;. ЭВМ обраоатывает булевские векторы, которые получа ются из исходных метеорологических данных путем их квантова
ния и кодирования. |
|
|
|
|
||
С помощью программы "Лавина- i " |
реализован алгоритм обу |
|||||
чения и доучивания. |
Входная информация для программы подго |
|||||
тавливается вспомогательной |
программой |
'‘М етео-Г', которая пре |
||||
образует |
метеорологические |
данные в |
конечное число матриц |
|||
классов |
и |
А■а . |
Матрицы класса А, |
содержат данные об ус |
||
ловиях схода |
лавин, |
иш сса |
А £ - |
об |
условиях нормирования |
снежной толки до критического момента. Вид матрицы класса А , соответствующий вектору Е ( 1 ) , приведен в таол. 2.
наличие нескольких матриц класса А, ж. А2 объясняет
ся следующей причиной. Размерность матрицы одного класса рав
на |
М |
N1 |
,где |
М |
- суммарные |
числа битов |
кодов всех парамет |
|
ров, |
а |
N sn .m |
|
,гдс |
П. |
-количество |
объектов данного |
|
класса, |
а |
гп |
- |
число |
подсобытий одного |
объекта.Разрядная |
сетка ЗЕш"Минск-22" ограничивает размерность числом 33, поэто
му продолжение матрицы, начиная со |
столбца ( в +36), где |
|||||
I = |
1 ,2 ,3 . . . . |
размешаются в следующих матрицах. Максималь |
||||
ная граница |
N |
= 72 • 36 т .е . матрица |
каждого класса,мо |
|||
жет быть не более 72. Размерность |
Is! , |
как правило, не огра- |
||||
кичивается. В результате булевский |
вектор представляет собой |
|||||
i -ю строку |
всех |
матриц данного класса. |
|
|||
|
Операторная схема программы "Лавина-Г’ имеет вид: |
|||||
где |
П, |
- |
ввод |
исходных данных; |
ф,, - начальное формирова |
|
ние; |
А( |
-обработка матриц лавин; |
О Д |
-организация цик |
||
ла по перебору сочетании по первой |
переменной ; P2Q ) -орга |
|||||
низация по перебору сочетаний по второй |
переменной ; |
Ф3(0 формирование цикла для перебора по третьей переменной
р ( I) - |
организация цикла по |
перебору |
сочетаний по третьей |
||
переменной; |
П2 - печать |
результатов; |
(С -конец работы; |
||
A j -обработка матрицы |
; |
А3 - построение |
решетки |
||
признаков матриц А, и А г |
. Б |
результате |
работы |
программы |
формируются и печатаются таксонометрические решетки указан ных матриц.
28
Достоинством программы "Лавина- I " по сравнению с аналогич ными программам /Бонгард,1967/ является то , что она не накладывает ограничений на число параметров,используемых для описания подсобытий, интервал квантования параметров,кодирова ние и число событий, на которых происходит обучение.
На основе экспериментального материала исследовалось 10 лавинных и 17 нелавинных случаев в зимние периоды (1967-
1972 г г . ) .
Для цроведения экзамена использовано девять событий, отме ченных на одном лавиносборе. Анализ проводился о нескольких позиций: исследовалась применимость алгоритма распознавания для прогнозирования схода лавин, а также оценивалась информатив ность параметров, правильность описания события, способ коди рования информации. Обучение проводилось на небольшом коли честве лавинных событий,поэтому трудно дать однозначный ответ о применимости метода распознавания для прогноза схода,хотя результат достаточно обнадеживающий.
Проверка применяемого метода на указанном материале пока зала, что общая его оправдываемость составляет 93$. Это позво ляет с уверенностью считать, что первую часть задачи, т .е .
распознавание лавинной и нелавинной ситуаций по комплексу до ступных измерению характеристик,программа выполняет успешно. Вторую чаоть задачи, которая предполагает непосредственно крат косрочный прогноз схода снежных лавин,можно решить при наличии получаемой синоптической ситуации, сопутствующей образованию снегонакопления на горных склонах. Необходимо отметить, что при регулярной обработке на ЭВМ поступающих данных с различных лавиносборов и анализе предыдущих случаев, вероятность правиль ных ответов будет, безусловно, повышаться.Анализ результатов, полученных с помощью программ " Лавина- I " .показал;
1)можно построить признаки, позволяющие разделять лавин ные и нелавинные события по косвенным характеристикам;
2)выбор параметров для описания события был правильным;
3)представление события с помощью слоиотой структуры позволяет выявить накапливаемые изменения в стратиграфии снега;
4)программа составлена достаточно гибко и позволяет про
водить обучение на любом количестве исходных .данных (до 100 событий). Небольшое изменение программы позволяет онять
это ограничение;
29