5.4. Гомологические модели функциональных связей

advertisement
5.4. Гомологические модели функциональных связей
Существование гомологических рядов данных, конгруэнтных закономерностях их
связи дают возможность решать специальные задачи, основываясь на полученных
выводах. Одна из таких задач – восстановление функций влияния реакции природных
систем на воздействия различных факторов.
5.4.1. Свойственные функции. При планировании хозяйственной деятельности
необходимо знать экстремальные особенности природной ситуации, в частности, уметь
рассчитывать вероятность появления негативных явлений и оценивать возможный ущерб
от их воздействия на конкретной территории. Эти проблемы рассматриваются в рамках
теории анализа и управления рисками, но, несмотря на большие усилия, методология и
методика анализа природных рисков мало разработаны [Рагозин, 1999]. Такое положение,
прежде всего, связано с трудностями понимания существа решаемых задач и отсутствием
моделей, адекватных природе явления.
Расчеты вероятностей проявления негативных последствий основаны на методах
исчисления рисков, в основу которых положены представления о свойственных функциях
геосистем (ландшафтов) и факторных системах, формирующихся из взаимодействующих
факторов среды и влияющих на реализацию этих функций [Черкашин, 2005] (см. п. 3.6.6).
Выявление функций, свойственных каждой геосистеме, отражающих закономерности
воздействия факторов среды на реакцию системы, изучаются разными методами с
использованием разнообразных моделей [Гуламов, Логофет, 1997, Пых, Малкина-Пых, 1997].
Свойственные функции выделяются по критериям полезности (или опасности)
геосистем для воспроизводства окружающей среды, природных ресурсов и обеспечения
хозяйственной деятельности. В этом смысле они близки по содержанию к представлениям
об экологической нише в понимании Ю. Одума [Одум, 1975]: экологическая ниша включает
не только физическое пространство, занимаемое организмом, но и функциональную роль
организма в сообществе и его положение относительно градиентов внешних факторов.
Представления о свойственных местоположению функциях объединяет как модели
типа «воздействие-реакция», «влияние-ущерб», так и модели оценки вероятности
возникновения критических значений факторов среды типа кривых обеспеченности
явлений (надежности, живучести). Во всех случаях речь идет о
функциях многих
переменных:
y = F ( x, u , n )
(1)
- непосредственно учитываемых факторов воздействия x; совокупности всех неучтенных в
расчетах
факторов
–
бонитета
u
(см.
п.
2.5.3);
набора
числовых
значений,
характеризующих классификационный тип природной среды – код ситуации n.
Функции такого вида рассматриваются в факторальной географии и факториальной
экологии, но до сих пор о них известно мало в силу их разнообразия, многофакторности
влияния
и
сложности
преимущественно
изучения,
массового
поэтому
характера,
любая
полезна
дополнительная
для
понимания
информация,
особенностей
свойственных функций, связанных с типологией и классификацией природных ситуаций.
Подобную информацию предоставляет космическая съемка в разных каналах спектра,
комплексно отражающая различные свойства земной поверхности, изменяющиеся от
места к месту. Растровые изображения удобны для сравнительного количественного
анализа с использованием ГИС-технологий и математических моделей. Показатели
яркости элементов изображений выше или ниже пороговых несут информацию о
критических ситуациях, например, о пожароопасности территории или, напротив,
избыточной увлажненности грунтов.
В
качестве
теоретической
основы
создания
функциональных
моделей
экологических ниш предлагается использовать методы теории конгруэнции [Черкашин,
1997, 2000], приводящие к уравнениям, адекватность которых оценивается по материалам
космической съемки для разных геосистем горно-таежных ландшафтов Лено-Ангарского
плато.
Выделяя на снимке однородные по функциональным характеристикам ареалы,
соответствующие определенным типам биогеоценозов или геосистем, на следующем
этапе появляется возможность дать их количественное описание, например, в виде
частотного распределения пикселов растрового изображения по яркости в границах
выделов. Для решения этой задачи проводится ординация
объектов в пространстве
показателей яркости некоррелированных каналов, т.е. строится модель ниши объекта в
гиперобъеме факторов, отраженных на снимке, что обычно осуществляется по методу
параллелепипедной классификации и сегментации снимка.
Объектом
исследования
стала
осваиваемая
горно-таежная
территория
Ковыктинского газоконденсатного месторождения (Жигаловский район Иркутской
области), где проводилась оценка воздействия геолого-разведочных работ на природные
экосистемы. Использовался космический снимок местности с аппарата LANDSAT (август,
2000г., разрешение 30 м). Осуществлена автоматическая обработка изображения с
выделением природных границ и его векторизация для создания ГИС. В ходе наземных
маршрутных исследований и вертолетных облетов территории осуществлялась заверка
обработанных снимков и составлена карта растительного покрова (биогеоценозов,
экосистем), представленная коренными лесами с преобладанием светлохвойных и
темнохвойных пород, гарями и разными стадиями восстановления коренных сообществ.
Каждый выдел карты растительности соответствует пространственной нише
оптимального (надежного) существования определенного набора явлений, свойственных
местной геосистеме. Структура этой ниши оценивается по космоснимку с помощью
распределения пикселов по яркости в разных каналах (рис 5.16). В качестве основного
использовали инфракрасный канал, который обычно несет информацию о критических
состояниях геосистем. Кривые распределения имеют колоколообразную форму и по
аналогии с термином «плотность приспособленности», принятого для видов, могут
называться «плотностью явления», имея в виду, что за характеристиками яркости снимка
стоят конкретные эколого-географические процессы и явления.
Рассматривается
модель плотного заполнения нишами пространства факторов
(ординации) с перекрытием (рис.5.17). Для каждой ниши выделяется центральная точка
оптимума и множества допустимых (свойственных) отклонений от нее, за пределами
(критическими явлениями) которой находятся несвойственные (рискованные) состояния
данного местоположения. Структура пространства ординации устроена так, что
критические явления одной ниши примыкают к точке оптимума другой: событие,
выходящее за пределы одной ниши попадают в сферу влияния другой. Ниши не
перекрываются
только
в
точках
оптимума.
Обычно
одной
точке
двухмерного
ординационного пространства соответствует три пересекающиеся области ниш.
Эта картина с учетом формы кривых плотности явлений проиллюстрирована на
рис. 5.18. Колоколообразные кривые выделяют допустимые области ниш в плоскости,
линии которой условно названы экобазой. Критические точки и точки оптимума
практически совпадают, поэтому каждая ниша имеет две базы – нижнюю пессимальную
критическую (экобазу) и верхнюю оптимальную базу. Отсюда следует, что теоретическое
моделирование кривых плотности явлений в виде свойственных функций должно дать
объяснение закономерностей строения (топологии) пространства ординации.
Эта задача решается в рамках теории функциональных систем [Черкашин, 1997,
2000], описывающей причинно-следственные связи (см. п. 3.6.6). В данном случае
исходим из гипотезы, что все функции (1) воспроизводятся на основе принципа
конгруэнции – существования пучка прямых с центом ( k 0 , k )
ln y − k 0 = A(ln x, u )(ln x − k ) ,
(2)
где k 0 = ln y 0 , k = ln x 0 - индивидуальные константы, связанные с естественной шкалой
x 0 , y 0 измерения переменных x, y - в данном случае значения яркости x и частоты его
встречаемости y (в долях от единицы):
ln( y y 0 ) = A ln( x x 0 ) , y = Bx A , ln B = ln y 0 + A ln x 0 = k 0 + Ak .
(3)
Согласно (2) линии-векторы, исходящие из центра пучка, при изменении
A(ln x, u ) вырисовывают сложные кривые (см. рис.5.18), моделирующие свойственные
функции, т.е. эти функции являются годографом вектора пучка. Центр пучка (узел
конгруэнции) лежит на функциональной кривой и рассматривается как особая
критическая точка распределения. Согласно принятой модели пространства ординации
критическая точка соседнего распределения совпадает
с оптимумом данного. Таким
образом, формальный смысл критических точек заключается в их свойстве быть узлами
конгруэнции векторов годографа кривой плотности явлений. Наличие такого свойства
говорит о гомологичности натурных рядов переменных y и x и возможности калибровки
их связей. Системой калибровки становится модель ниши (2).
В формуле (3) при постоянном значении А зависимость
y ( x ) - степенная
(аллометрическая) функция, аппроксимирующая сложные кривые в окрестностях узла
конгруэнции. Появление конгруэнции для линеаризованных степенных функций связано с
существованием линейной зависимости коэффициентов А и ln B [Lumer, 1939, Кофман,
1982]. Это свойство рассматривается как фундаментальная закономерность, позволяющая
сводить функциональные связи к простейшему виду, и, как следует из (3), сохраняющаяся
для случая, когда А становится функцией нескольких переменных. Именно поиск вида
этой функции А, управляющей ориентацией образующих векторов, является основной
задачей в данном случае.
5.4.2. Подобие функциональных ниш. Для обоснования гипотез оценки
параметров используется метод сравнения распределений путем совмещения их точек
оптимума. Формально это достигается
изменением значений ln x
∆ k = k 2 − k1 . Таким образом, два распределения y1
на величину
и y2 оказываются сравнимы в
логарифмических координатах:
ln y1 − k 01 = A1 (ln x, u )(ln x − k1 ), ln y2 − k02 = A2 (ln x + ∆ k , u )(ln x + ∆ k − k 2 ) .
Теперь разделив почленно второе уравнение на первое с учетом k1 = k 2 − ∆ k ,
получим:
ln y 2 − k 02 = [ A2 (ln x + ∆ k , u ) / A1 (ln x, u )] ⋅ (ln y1 − k 01 ) .
(4)
Сопоставление распределений в широкой окрестности их оптимума дают
практически
прямые
линии
(рис.5.19),
что
означает:
отношение
A2 (ln x + ∆ k , u ) / A1 (ln x, u ) = К 12 - константа. Из (4) следует
ln y 2 = K 12 ln y1 + С , C = − k 01 K 12 + k 02 ,
и
свободный
член
C
здесь
(5)
является
линейной
функцией
К 12 :
C = − 2.93K 12 − 2.94 . В этом случае все линии вида (5) пересекаются в одной точке, т.е.
образуют конгруэнцию, в которой в силу (5) и условности выбора и y1 и y2, справедливо
k 01 ≡ k 02 - это значение становится экобазой для всех распределений плотности явлений.
Расчеты подтверждают эту гипотезу: k 01 = − 2,93, k 02 = − 2,94 . Таким образом,
критический уровень явлений для всех природных условий одинаковый, а выявленные
зависимости относятся к одной фундаментальной закономерности (конгруэнции). Кроме
того, все распределения образуют гомологический ряд подобия свойственных функций,
допускающих калибровку с учетом смещения центров.
Основной вопрос – что представляет собой функция A(ln x, u ) , сохраняющая
подобие
при
смещении
значения
аргумента
(свойство
автомодельности):
К 12 A(ln x, u ) = A(ln x + ∆ k , u ) . Она с точностью до коэффициента пропорциональности
повторяется через постоянный интервал значений аргумента ∆ k . Этому условию
соответствуют постоянные, периодические и экспоненциальные функции. Объединяющей
является формула экспоненты с комплексным коэффициентом а
A(ln x, u ) = A0 (u ) exp(a ln x ), a = α + ω i, i =
− 1,
(6)
Где, согласно уравнению Эйлера, exp(a ln x) = eα ln x [cos(ω ln x) + i sin(ω ln x)] . При α=0
для
действительной
части
выражения,
стоящего
в
скобках,
получаем
A(ln x, u ) = A0 (u ) cos(ω ln x ) . По пересечению линии экобазы k 0 ≈ − 2,93 с кривыми
плотности явлений, построенных в логарифмической шкале, определяется положение узла
каждой конгруэнции k = k n и из формулы (2) рассчитывается функция A(ln x, u ) . Она
действительно представляет собой фрагмент периодической функции с периодом
T = 2π / ω ≈ 1 и частотой ω = 2π / Т ≈ 6,28 . Следовательно, значения
отличаются на 1
kn
примерно
( k n ≈ n ). Это позволяет говорить о дискретности (квантованности)
пространства ординации (см. рис.5.21). При T ≈ 1 положение максимумов соседних
колоколообразных кривых различаются по значениям яркости (х) примерно в
exp(1) = 2,718 раза, т.е исходные кривые вытягиваются вдоль оси х. Как видим, в этих
закономерностях одновременно и естественно проявляются фундаментальные константы
1, π , e
.
Дополнительная экспоненциальная составляющая e α ln x в уравнении (4) дает
возможность учесть затухание колебаний и несколько повысить точность расчетов,
особенно на концах распределения. Таким образом, общий вид расчетной формулы,
моделирующей структуру экологической ниши, следующий
ln y n − k 0 = A0 (u )eα ln x cos(ω ln x)[ln x − Tn] ,
(7)
где n – порядковый номер функциональной кривой в ряду (см. рис.5.21). Оценке подлежит
только функция от бонитета A0(u), которая определяется влиянием неучтенных факторов.
Значение u является своеобразным гомотопическим параметром, переводящим одну
функцию связи в другую. Таким образом, зависимость (4.27) раскрывается как сложная
функция трех параметров x, u и n:
 x 
yn = e  Tn 
e 
k0
A0 ( u ) xα cos(ω ln x )
(8)
Таким образом, основываясь на свойствах конгруэнции функциональных зависимостей
плотности явлений от факторов среды, предложен и проверен метод анализа свойственных
функций, основанный на математическом моделировании структуры природных ниш с
применением данных распределения яркостных характеристик космических изображений.
Доказано существование инвариантного уровня критических изменений - экобазы и
дискретного размещения вдоль линии экобазы узлов конгруэнции, являющихся началом
вектора годографа формирования кривых распределений. Обоснован экспоненциальный
комплексный вид функций управления ориентацией образующих векторов (функций
чувствительности). На этой основе разрабатывается методика оценки уровня критических
изменений (на уровне экобазы), переводящего систему одного типа в систему другого (см.
п. 6.3).
Рис. 5.16.
Рис.5.17.
Рис. 5.18.
Рис. 5.19.
10
12
10 12
12
30
11
12
10
12
12
18
19
12
10
5
12
18
16
27
31
1
27
11
16
22
10
12
27
18
18
19
25
18
18
18
18
1
25
22
23 18
25
17
4
5
1
18
1
23
19
17
1
1
9
10
12
11
27
12
18 11
1
23
181
18
1
18
5
23
23
1
1 23
24
27
Рис.5.21.
27
12
27
30
19
31
11
12
9
4
1
63 1
30 18
6
30
20
31
21
20
31
31
21
31
30
21
31
11
7
21
30
21
1
20
1
31
30 21
27
27
27
18
1
10
30
27
11
25 18
18
21
16
27
30
23
27
18
12
27
31
27
20
30
18
1
31
21
23
31
22
30
6
30
9
18
21
1
13
30
31
30
8
.
7
7
7
8
8
6
4
4
6
6
13
7
8
7
.
2
6
6
8
2
3
11
8
7
7
2
6
8
6
8
6
2
2
8
8
14
9
31
11
8
8
4
3
7
6
4
3
3
7
28
30
31
27 22
5
6
7
2
3
6
1
2
6
1 31
5
6
2
6
11
1
6
14
13
3
248
22
2
24
20
2
23
28
31
21
5
4
24
12
30 29
9
18
30
31
23
18
10
22
25
27
2 2 2 72 9
27
27
10
20
31
1
19
24
27
27
18
12
31
6
30
22
27
27
31
20
12
11
28 11
2
1
31
27
31
31
27
22
27 24
23
30
24
6
3
3
2
5
1
24
31
5
1
12
22
27
1
11
22
16
5 1
27
22
11 23
28
30
19
2
3
2
5
2
2
1
1
31
2 52 3 3 0
22
20
27
30 23
24
9
1
18
11
27
12
27
1
22
24
10 5
27
1
2
12
27
12
27
23
25
18
16
1
23
15
27
16
23
23
30
30
24
30
24
31
15
23
23
27
23
27
22
11
27
23
27
27
5
22
27
27
18
18 27
18
27
2
5
5 3
2
1
2
1
2
2
22
2
5
2
27
27 23
23
27 31
19
31
1
23
27
31
16
23
25
22
3
1
4
3
27
1
27
27
27
18
27
23
1
17
1
27
27
1
5
1
1
4
27
31
30
1
27
27
1
5
23
25
6
2
3
2
23
24
29
27
18
23
24
19
17
16
18
24
27 22
22
18
30 23 12
18
17
18
27
27 23 27
27
23
27
24
1 72 2
24
23
19
19
27
1 81 7
17
23
27
11
12
27
18
19
17
22
27
18
15
25
27
2
1
3
8
22
1
5
2
4
22
3
11
3
30
1
28
23
23
18
11
1
3
5
3
3
1
25
27
23
23
27
3
1
3
11
1
27
22
27
19
23
27
18
16
22
10
18
17
19
17
15
19
22
23 27 23
22
3 01 2
22
5
5
1
22
22
15
1
2
5
27
27
27
27
16
18
23
22
15
11
15
22
27
31
18
1
27
27
23
1 62 3
27
1
19
17
31
17
11
22
11
18
18
23
1
23
23
19
19
22
27
27
23
19
31
27
18 22
22
23
23
12
25
18
27
23
15
27
19
18
12
1
22
15
27
18
18
12
12
12
23
18
-
Download