УДК 551.4 (470.621)

ББК 26.823 (2 Рос.Ады)

С 87

Варшанина Т.П., Плисенко О.А.,

Коробков В.Н., Солодухин А.А.

Центр интеллектуальных геоинформационных технологий НИИ комплексных проблем АГУ

 

Структурно подобная модель рельефа для точечного прогнозирования геоморфологических процессов

Аннотация. Разрабатывается структурно подобная геоинформационная модель рельефа на основе метода структурной маски энергетических географических полей и информационно-математической структурной модели трёхмерной геометрической поверхности рельефа. Модель предназначена для моделирования динамических свойств геоморфосистем, точечного прогнозирования геопространственных процессов.

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

 

Varshanina T.P., Plisenko O.A.,

Korobkov V.N., Solodukhin A.A.

 

STRUCTURALLY SIMILAR MODEL OF THE TOPOGRAPHY FOR POINT FORECASTING OF GEOMORPHOLOGICAL PROCESSES

 

Abstract. A structurally similar geoinformation model of the topography is elaborated using a method of a structural mask of power geographical fields and the information and mathematical structural model of a three-dimensional geometrical surface of the topography. The model is intended for modeling of the dynamic properties of geomorphic systems and for point forecasting of geospatial processes.

Keywords. Structurally similar model of topography, method of a structural mask of power geographical fields, information and mathematical structural model of a three-dimensional geometrical surface of the topography, modeling of the dynamic properties of geomorphic systems, point forecasting of geospatial processes.

 

Рельеф физической земной поверхности, формирующийся в результате динамических процессов происходящих в теле Земли в поле космопланетарных энергетических воздействий, является одним из главных факторов дифференциации объектов географического пространства. Вследствие этого создание структурно подобной модели рельефа удовлетворительно соответствующей иерархии эмпирически объективных сущностей геопространства является решающим шагом на пути осуществления сущностно-логической структуры данных в глобальном проекте «Электронная Земля» и успешного прогнозирования широкого спектра геопространственных процессов.

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

К условиям, ограничивающим область понятий сущностно-логической модели рельефа, как и любой модели реального мира, относятся следующие [1]:

- иерархия понятий, отображающих структуру рельефа, должна соответствовать иерархии эмпирически объективных географических сущностей;

- модель рельефа (как сложного самоорганизующегося объекта) должна быть свернута в модель, отражающую на каждом структурном уровне параметры порядка самоорганизации (структурообразующие параметры);

- структурная модель рельефа по классификации подобий может быть только практически приближенно подобной моделью.

В онтологии знаний о рельефе, как природном образовании в системе геопространства, базовым является понятие геоморфосистема.

Общепринятого определения понятия «геоморфосистема» нет. А.Н. Ласточкин [7, 8] трактует это понятие как форма земной поверхности, описываемая системой элементарных поверхностей. На наш взгляд, геоморфосистема – это участок земной поверхности, описываемый системой элементарных поверхностей и квазиоднородный по специфике и интенсивности рельефообразующих процессов.

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

В настоящее время отсутствуют способы регистрации в фокусе земной поверхности интегральных геофизических полей, вследствие этого в данной работе предлагается вычислять их структуру и динамику методом структурной маски [2]. Метод основан на следующих логических построениях. Так как структура каждой геоморфосистемы формируется порождающим её полем энергии, структурообразующие элементы геоморфосистемы – параметры поля высот, являются мерой напряженности реального геофизического поля, т.е. мерой параметра порядка его структурной организации. Таким образом, структура географических энергетических полей определяется градиентами структурообразующих иерархию географических систем величин, трактуемых в качестве меры их параметров порядка.

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

Алгоритм вычисления структурной маски поля эндогенной энергии Пространственно-временная структура поля эндогенной энергии в фокусе земной поверхности воссоздается через наибольшие градиенты мощности слоя земной коры (разностный слой) между полибазисными поверхностями рельефа, фиксирующими положение поверхностей выравнивания соответствующего геологического времени [3, Философов]. Динамика положения полибазисных поверхностей соответствует ритмам энергетического режима тектонических течений. Мощность разностного слоя четвертичного возраста указывает на скорость вертикальных движений в каждой точке земной поверхности за четвертичное время. Наибольшие градиенты мощности разностного слоя в каждой точке земной поверхности воссоздают структуру современного поля эндогенной энергии в фокусе земной поверхности.

Алгоритм вычисления структурной маски поля эндогенной энергии в фокусе земной поверхности реализуется комплексом методов с помощью геоинформационных инструментов на примере территории Республики Адыгея. Мощность разностного слоя между топографической поверхностью и базисной поверхностью первого порядка (четвертичного возраста) отображается в изолиниях одинаковой мощности разностного слоя – изогипсопахитах. Строится равномерная решетка с постоянным шагом, в каждом узле которой вычисляется наибольший градиент мощности разностного слоя (рис.1) структурно подобный напряженности эндогенного энергетического поля.

Поле эндогенной энергии создает структурно подобную ему систему тектонических элементов территории – систему блоков и ограничивающих их разломов. Эти элементы на земной поверхности достоверно идентифицируются по комплексу индикационных признаков [11] и визуально по матрице высот (рис. 2). Вычисление результирующего вектора наибольшего градиента мощности разностного слоя для каждого тектонического элемента, где величина относительной скорости равна величине градиента, а направление – вектору (рис.3), показало высокую степень соответствия полученной геодинамической модели территории её известным тектоническим и тектонофизическим моделям [10]. Верификация геодинамической модели произведена также путём сравнения соответствия модельных и инструментальных данных.

Доказательством структурного подобия модельных данных, вычисленных с помощью разработанной геодинамической модели, и инструментальных данных, представленных на картосхемах пространственной дифференциации скоростей вертикальных тектонических движений на территории исследования по данным повторного нивелирования за период 1925-1992гг [9]. является высокий и достоверный коэффициент корреляции между ними (табл. 1).

Таблица 1

Коэффициент корреляции между модельными и инструментальными данными скорости вертикальных тектонических движений на территории Адыгеи

Системно-структурная геодинамическая модель

Инструментальные данные

средние значения

максимальные значения

средние значения

0,78

t=8,12%*

0,74

t=2,35%*

максимальные значения

0,81

t=9,30%*

0,70

t=8,63%*

* При tтабл = 2,576

 

Скорость вертикальных движений тектонических блоков по инструментальным данным определялась статистически по Grid-поверхностям, разработанным по соответствующим оцифрованным картосхемам. Модельные данные скоростей вертикальных движений тектонических блоков в мм/год были рассчитаны с учетом продолжительности периода, которому принадлежит базисная поверхность. Вычисление коэффициента корреляции производилось между максимальными и средними значениями скоростей вертикальных движений тектонических блоков в мм/год, полученными по модельным и инструментальным данным.

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

Разработка приближено подобной трёхмерной геометрической поверхности рельефа возможна по классификации А.Н. Ласточкина [5]. На основании этой классификации разработаны алгоритмы построения цифровой информационно-математической модели, позволяющие выделять и классифицировать полную группу структурных элементов поверхности рельефа (структурные линии: гребневые (L1), килевые(L2), максимальных (L3) и минимальных (L4) уклонов, выпуклых (L5) и вогнутых (L6) перегибов склонов) и элементарных поверхностей (структурных единиц геоморфосистем).

Автоматизированная технология выделения геоморфосистем основана как на использовании строгих вычислительных методов, так и методов искусственного интеллекта для решения слабоформализованных задач. Обеспечено автоматизированное определение положения характерных точек поверхности рельефа (рис. 1), по которым строятся структурные линии (рис. 2). Классификация выделенных элементарных поверхностей производится с помощью классификационного конструктора (рис. 3). Комплекс геоморфологических параметров каждой поверхности хранится в объектно-ориентированной базе данных, обеспечивающей сопряженность с данными о параметрах эндогенного поля энергии, характере литологического основания и морфоклиматической принадлежности иерархии объектов поверхности рельефа.

 1

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

2 

Рис.2. Автоматизированное выделение и типизация характерных и особых точек, характерных линий по изолинейным картам.

32

Рис.3. Основные формы классификационного конструктора

 

Морфоклиматическое районирование территории исследования произведено по схеме, предложенной В.И. Кружалиным [6]. В Адыгее выделено четыре морфоклиматических области с характерными для них рельефообразующими факторами, обусловливающими тип и интенсивность действующих геоморфологических экзогенных процессов [1, стр. 229].

Алгоритм вычисления структурной маски поля экзогенной энергии.

 Информационно-математическая модель трёхмерной геометрической поверхности рельефа идеальна для расчёта структурно подобного поля приходящей к нему потенциальной экзогенной энергии, которое является продуктом интерференции энергетических потоков как проходящих через атмосферу извне, так и формирующихся в ней. Структуру поля приходящей экзогенной энергии определяет прямая радиация, так как рассеянная поступает со всего небосвода. Перерасчёт в рельефе прямой радиации производится в авторской программе «San» [3]. Суточная сумма прямой радиации, поступающей на склоны с заданными уклоном и экспозицией, определяется с интервалом 15 минут. На каждом шаге учитывается высота и азимут Солнца, количество прямой радиации, падающее на площадку, перпендикулярную солнечным лучам. Вычисление производится на основе аппроксимированной функции определяющей рассеяние света в атмосфере, таблично заданных функций распределения солнечной энергии в спектре солнечного излучения на внешней границе атмосферы, а также спектра поглощения атмосферы. На основе полученных результатов рассчитывается средняя суточная сумма прямой радиации на склонах основных экспозиций и ее отношение к количеству радиации, поступающей на горизонтальную поверхность.

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

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

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

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

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

- границ иерархии геоморфосистем;

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

- границ структурных элементов геологической среды и почв;

- экологических ниш биоты и ландшафтов;

- точечного прогнозирования экзогенных геоморфологических процессов по данным мониторинга.

Построение и анализ пространственно-временного поля эндогенной энергии территории исследования

С помощью геоинформационного модуля вычислена пространственно-временная структура поля эндогенной энергии на территорию Адыгеи [4]. Относительно автоматизировано ранжированной гидросети построены полибазисные поверхности 2,3,4 и 5 порядков, отражающие динамику поверхности рельефа за плиоцен-четвертичное время.

Карты остаточного рельефа построены по методике В.П. Философова. С некоторой долей условности базисная поверхность 4-го порядка отнесена к началу, 3-го – к концу плиоцена; 2-го – к началу, 1-го – к концу четвертичного периода. Построение поверхностей остаточного рельефа соответствующих порядков показало, что эндогенное поле энергии в реологических условиях геологической среды уже в плиоцене сформировало основные тектонические элементы территории – различного ранга тектонические блоки и линеаменты (рис. 4). От плиоцена к четвертичному времени увеличивалось дробление территории, расширялись области опусканий, приуроченные в основном к северной части республики (рис. 5, 6, 7).

 4

Рис.4. Остаточный рельеф относительно базисной поверхности 4-го порядка.

 5

Рис.5. Остаточный рельеф относительно базисной поверхности 3-го порядка.

6 

 

Рис.6. Остаточный рельеф относительно базисной поверхности 2-го порядка.

 7

Рис.7. Остаточный рельеф относительно базисной поверхности 1-го порядка.

 

Структура эндогенного поля энергии построена на основании исследования дифференциации тектонических движений за плиоцен-четвертичное время.

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

С этой позиции проанализированы поля наибольших градиентов мощности «разностного слоя» между базисными поверхностями 2, 3, 4 порядков, спроецированные на современное блоковое строение территории. Градиент мощности разностного слоя измерялся в метрах на единицу расстояния в один метр.

Векторное поле направления и относительной скорости горизонтального перемещения блоков БП 4-3 (по данным мощности «разностного слоя» между базисными поверхностями БП-4 и БП-3, рис. 8) показало, что в плиоцене блоки в северо-западной части исследуемой территории отличались малой относительной скоростью горизонтального движения, которое было направлено к югу, в сторону осевой зоны Большого Кавказа. Эта тенденция направления горизонтального движения блоков сохранялась в западной части территории до долгоживущего коро-мантийного Ахтырского разлома. На остальной части территории картина поля тектонических движений все более усложнялась в направлении Главного Кавказского хребта, свидетельствуя о сложной системе стесненных деформаций.

8 

Рис.8. БП 4-3: направление/скорость движения блоков.

 

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

Векторное поле тектонических напряжений БП 3-2 (рис. 9) свидетельствует о том, что на территории от Азово-Кубанского прогиба до Ахтырского разлома и Северо-Юрской депрессии, если и есть некоторые изменения по сравнению с БП 4-3 в направлении движения блоков, то в целом сохраняется уровень интенсивности тектонических взаимодействий.

9 

Рис.9. БП 3-2: направление/скорость движения блоков.

 

К югу от Северо-Юрской депрессии по сравнению с БП 4-3 интенсивность тектонических взаимодействий уменьшается. Нагорье Лагонаки продолжает выступать в качестве концентратора тектонических напряжений.

В конце четвертичного периода (БП 2-1, рис. 10) интенсивность тектонических взаимодействий в среднегорно-высокогорной части республики продолжает снижаться. Нагорье Лагонаки остается концентратором тектонических напряжений.

10 

Рис.10. БП 2-1: направление/скорость движения блоков.

 

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

Пространственно-временная дифференциация тектонического напряжения в пределах тектонических элементов отображается в результате построения градиентного поля относительно изогипсопахит (изолиний равной мощности) разностного слоя в узлах решетки с длиной стороны ячейки 250 м (рис. 11,12).

11 

Рис.11. Градиенты наибольшей мощности «разностного слоя», построенные относительно его изогипсопахит.

 

12 

 

Рис.12. Дифференциация тектонических напряжений в МС блоках и линеаментах, отображаемое величиной и направлением наибольшего градиента разностного слоя: а – плиоцен, б – плиоцен-четвертичное время, в – четвертичный период.

 

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

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

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

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

 13

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

 

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

«Хорошо структурированная» системная геодинамическая модель территории необходимое основание для применения методов и инструментария теории топологии:

- для поиска количественных соотношений «структура–свойство» и «структура-активность»;

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

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

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

- для поиска критических точек неустойчивости («разветвления»);

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

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

Разработан инструментарий для поэлементного анализа и прогнозирования характеристик саморазвивающихся геоморфосистем (рис.14, 15).

 

14 

Рис.14. Обобщенная схема классов модулей построения и анализа поведения геосистемы на примере подсистемы «Геодинамическая модель».

 

15 

Рис.15. Главное окно программы Моделирование динамических свойств геосистем.

 

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

– для исследования пространственно-временных закономерностей полей напряжений, формируемых в земной коре по иерархии блоков актуотектонической модели, и выявления концентраторов напряжении;

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

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

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

 

Литература:

  1. (Варшанина Т.П., Плисенко О.А. Интегрированная ГИС региона (на примере Республики Адыгея). Под ред. Б.И. Кочурова. – Москва – Майкоп: Издательский дом «Камертон», 2011. – 360 с.)
  2. Варшанина Т.П. Разработка хорошо структурированной модели геопространства на основе метода структурной маски энергетических полей. Вестник Адыгейского государственного университета. Серия «Естественно-математические и технические науки». Майкоп: изд-во АГУ, 2012. – Вып. 4 (110). С 176-180.
  3. Варшанина Т.П., Митусов Д.В. Климатические ресурсы ландшафтов Республики Адыгея. – Майкоп: Изд-во Адыгейского госуниверситета, 2005. – 237 с.
  4. Варшанина Т.П., Плисенко О.А., Солодухин А.А., Коробков В.Н. Структурно-подобная геодинамическая модель Краснодарского края и Республики Адыгея: Под ред. Б.И. Кочурова. – Москва – Майкоп: Издательский дом «Камертон», 2011. – 128 с.
  5. Зинченко А.Г., Ласточкин А.Н. Методика геоморфологического картографирования шельфа и континентального склона Российской Федерации. М.: ЗАО «Геоинформмарк», 2001. 38 с.
  6. Кружалин В.И. Экологическая геоморфология суши. / В.И. Кружалин. – М.: Научный мир, 2001. – 176 с.
  7. Ласточкин А.Н. Морфодинамическая концепция общей геоморфологии. / А.Н. Ласточкин. – Л.: Издательство Ленинградского университета, 1991. – 220 с.
  8. Ласточкин А.Н. Системно-морфологическое основание наук о Земле. / А.Н. Ласточкин. // Геотопология, структурная география и общая теория геосистем. - СПб.: Изд-во НИИХ СПбГУ, 2002. – 762 с.
  9. Лилиенберг Д.А. Картографо-геодинамические модели современных движений морфоструктур Кавказа и Закавказья для разных эпох. / Д.А. Лилиенберг. // Геоморфология. 1997. - Вып.4. С. 63-75.
  10. Милановский Е.Е. Новейшая тектоника Кавказа. / Е.Е. Милановский. М.: Недра, 1968. - 482 с.
  11. Ранцман Е.Я., Гласко М.П. Морфоструктурные узлы. – места экстремальных природных явлений. – М.: Медиа-Пресс, 2004. - 224 с.

___________________________________________________________________________

 

Варшанина Татьяна Павловна, доцент, кандидат биологических наук, заведующая Центром интеллектуальных геоинформационных технологий НИИ комплексных проблем АГУ, Федеральное государственное образовательное учреждение высшего профессионального образования Адыгейский государственный университет, 385000,г. Майкоп, ул. Первомайская 208, 8 (8772) 59-39-44, Этот адрес электронной почты защищен от спам-ботов. У вас должен быть включен JavaScript для просмотра. .

Плисенко Ольга Анатольевна, зав. сектором Центра интеллектуальных геоинформационных технологий НИИ комплексных проблем АГУ,  385000,г. Майкоп, ул. Первомайская 208, 8 (8772) 59-39-44, Этот адрес электронной почты защищен от спам-ботов. У вас должен быть включен JavaScript для просмотра. .

Коробков Виктор Николаевич, программист Центра интеллектуальных геоинформационных технологий НИИ комплексных проблем АГУ, 385000,г. Майкоп, ул. Первомайская 208, 8 (8772) 59-39-44, Этот адрес электронной почты защищен от спам-ботов. У вас должен быть включен JavaScript для просмотра. .

Солодухин Александр Александрович, программист Центра интеллектуальных геоинформационных технологий НИИ комплексных проблем АГУ, 385000,г. Майкоп, ул. Первомайская 208, 8 (8772) 59-39-44, Этот адрес электронной почты защищен от спам-ботов. У вас должен быть включен JavaScript для просмотра. .

Центр интеллектуальных геоинформационных технологий НИИ комплексных проблем Адыгейского государственного университета

Яндекс.Метрика
© Адыгейский государственный университет. НИИ комплексных проблем.