Диагностирование и прогноз экологического состояния природно-технической гидрогеологической системы
На мой взгляд, есть несколько способов повышения эффективности модели, а именно: — замена показателей, не полностью связанных с функцией отклика и не сильно влияющих на ее изменение. В нашем случае, такими признаками, на мой взгляд, являются защищенность и вертикальный градиент, т.к. он имеет важную роль в местах перетекания, а таких мест на карте не много, поэтому его роль незначительна… Читать ещё >
Диагностирование и прогноз экологического состояния природно-технической гидрогеологической системы (реферат, курсовая, диплом, контрольная)
1. Формирование исходной модели
1.1 Стратиграфия и литология изучаемой территории
1.2 Гидрогеологические условия
1.3 Техногенные условия
2. Формирование информационной модели
3. Тестирование модели Заключение Список литературы
Целью данной работы является проведение диагностирования и составления прогноза состояния природно-технической гидрогеологической системы подольско-мячковского горизонта по содержанию хлора в его водах.
Основными задачами являются: анализ территории по набору гидрогеологических, геологических и техногенных условий состояния изучаемой территории; выяснение влияния имеющихся условий на загрязнение подольско-мячковского горизонта хлором; формирование информационно-диагностической модели; тестирование сформированной модели для прогноза состояния подольско-мячковского горизонта.
1. Формирование исходной модели
1.1 Стратиграфия и литология изучаемой территории
Исследуемая территория расположена в восточной части Московской области. Стратиграфия и литология изучается на основе гидрогеологического разреза. В геологическом строении данного района представлены следующие стратиграфические единицы: отложения среднего и верхнего карбона, средней юры, нижнего мела и четвертичные отложения.
Средний отдел карбона представлен глинистыми известняками подольско-мячковского горизонта и глинами. Мощность известняков по разрезу изменяется от 34 до 39 м., а глин >2 м. Верхний отдел карбона представлен известняками касимовского водоносного горизонта и пестроцветными глинами. Мощность касимовского горизонта около 50 м. Мощность глинистых отложений составляет 8−12 м.
Средний отдел юры на изучаемой территории представлен келловейским ярусом, который сложен черными глинами мощностью на западе 15 м., на востоке 9 м.
Нижний мел на данной территории представлен аптским ярусом, сложенным песками мощностью 11 м. Пески обводнены.
Отложения четвертичной системы имеют широкое распространение и представлены пролювиальными, флювиогляциальными и ледниковыми осадками. Среднечетвертичные ледниковые отложения представлены: моренными суглинками мощностью на западе 21 м., на востоке 11 м., флювиогляциальными песками и супесью. Пролювиальные отложения представлены суглинками мощностью от 2 до 5 м.
1.2 Гидрогеологические условия
На изучаемой территории развито два водоносных комплекса: касимовский и подольско-мячковский.
Водовмещающими породами касимовского горизонта являются известняки с прослоями глинистых известняков мощностью 48−49 м. Водоносный горизонт напорный, верхним водоупором служат юрские черные глины, а нижним пестроцветные глины. Пьезометрический уровень расположен на глубине в западной части разреза 42 м, в восточной — 31 м.
Водовмещающими породами подольско-мячковского горизонта также являются известняки мощностью 34−39 м. Водоупор представлен глинами. Пьезометрический уровень этого горизонта расположен на глубине 45 м, то есть ниже уровня касимовского горизонта.
Пьезометрический уровень подольско-мячковского горизонта расположен ниже уровня касимовского горизонта. Таким образом, возможен переток воды через водоупор из касимовского горизонта в подольско-мячковский, а следовательно, подольско-мячковский водоносный горизонт может быть подвержен загрязнению. Используя данные по скважинам был построен геолого-гидрогеологический разрез по линии 1−13.
В качестве исследуемых показателей были приняты следующие: защищенность подольско-мячковского горизонта; суммарная защищенность подольско-мячковского горизонта; положение пьезометрического уровня касимовского водоносного горизонта; положение пьезометрического уровня подольско-мячковского водоносного горизонта; вертикальный градиент, определяющий направление и интенсивность перетекания между касимовским и подольско-мячковским водоносными горизонтами; водопроводимость подольско-мячковского водоносного горизонта; функциональный тип территории; уровень концентрации хлоридов в подземных водах подольско-мячковского водоносного горизонта; По первым четырем показателям были построены карты на основе данных по скважинам. Также были построены карты водопроводимости, типа территории и концентрации хлоридов в воде.
1.3 Техногенные условия
Анализ рисунков 6, 7 и 8 показал, что для территорий с/х назначения загрязнение хлоридами подольско-мячковского горизонта наибольшее и превышает 3 ПДК, но есть небольшая территория на северо-западе района, где содержание хлоридов < 1 ПДК. Что, видимо, связано с хорошей защищенностью на данной территории водоносного горизонта.
Территория индустриального назначения характеризуется содержанием хлоридов < 1 ПДК и 1−3 ПДК.
Селитебным районам соответствует содержание хлоридов в подольско-мячковском горизонте < 1 ПДК и 1−3 ПДК, но есть небольшая зона, где это содержание превышает 3 ПДК. Видимо, это связано с тем, что здесь сосредоточены большие участки с/х и индустриального использования.
Таким образом, можно сделать вывод, что максимальное загрязнение хлоридами водоносного горизонта происходит в основном в области распространения с/х угодий.
Для формирования информационной модели имеющиеся карты значений показателей разбили на 96 блоков, размер каждого составил 1×1 км, и вычислили средние значения этих показателей по блоку.
Вертикальный градиент вычислялся как разность уровней касимовского и подольско-мячковского водоносных горизонтов, отнесенная к мощности кривякинского водоупора: Подсчитанные средние значения показателей представлены в таблице 1. Таблица 1. Средние значения показателей по блокам № блока | Суммарная защищенность | Защищенность | Уровень | Уровень | Вертикальный градиент | |
25,6 | 38,3 | — 5,3 | ||||
24,5 | 3,9 | 34,1 | 32,4 | 0,435 897 | ||
28,9 | 7,5 | 32,8 | 31,9 | 0,12 | ||
32,5 | 11,4 | 32,1 | 33,2 | — 0,9 649 | ||
33,7 | 14,8 | 31,1 | 34,2 | — 0,20 946 | ||
33,1 | 17,0 | 31,4 | 34,7 | — 0,19 412 | ||
33,0 | 14,6 | 31,2 | 36,3 | — 0,34 932 | ||
32,2 | 12,6 | 30,1 | 38,9 | — 0,69 841 | ||
33,3 | 14,4 | 33,3 | 38,1 | — 0,33 333 | ||
37,0 | 19,1 | 36,6 | 37,1 | — 0,2 618 | ||
39,0 | 22,5 | 40,5 | 38,4 | 0,93 333 | ||
39,2 | 24,0 | 0,41 667 | ||||
25,2 | 34,3 | 39,1 | — 4,8 | |||
23,8 | 3,3 | 32,8 | 33,4 | — 0,18 182 | ||
26,7 | 6,8 | 31,1 | 29,8 | 0,191 176 | ||
32,0 | 10,6 | 29,9 | 31,0 | — 0,10 377 | ||
37,3 | 14,0 | 28,4 | 32,3 | — 0,27 857 | ||
38,3 | 15,5 | 28,1 | 33,6 | — 0,35 484 | ||
34,9 | 13,5 | 29,6 | 37,8 | — 0,60 741 | ||
32,7 | 11,7 | 29,8 | 40,6 | — 0,92 308 | ||
36,3 | 13,6 | 31,7 | 40,3 | — 0,63 235 | ||
39,1 | 18,2 | 37,1 | — 0,11 538 | |||
42,5 | 22,3 | 39,1 | 38,1 | 0,44 843 | ||
42,6 | 23,6 | 42,5 | 40,3 | 0,9 322 | ||
25,4 | 1,5 | 33,9 | 41,2 | — 4,86 667 | ||
23,5 | 4,2 | 32,2 | 33,8 | — 0,38 095 | ||
25,7 | 6,6 | 30,8 | 32,1 | — 0,19 697 | ||
31,0 | 10,8 | 29,3 | — 0,43 519 | |||
34,9 | 14,3 | 28,1 | 36,1 | — 0,55 944 | ||
33,6 | 14,2 | 32.8 | — 0,43 662 | |||
32,5 | 12,5 | 36.2 | 42,4 | — 0,496 | ||
31,4 | 10,5 | 33.9 | 43,7 | — 0,93 333 | ||
34,4 | 12,4 | 31.6 | 39,8 | — 0,66 129 | ||
37,9 | 16,9 | 33.3 | 36,4 | — 0,18 343 | ||
40,1 | 20,1 | 37.1 | 36,4 | 0,34 826 | ||
41,1 | 22,3 | 40.6 | 36,1 | 0,201 794 | ||
25,1 | 2,3 | 33.3 | 42,1 | — 3,82 609 | ||
23,2 | 5,5 | 32.1 | 39,1 | — 1,27 273 | ||
25,0 | 9,9 | 31.4 | 38,4 | — 0,70 707 | ||
28,5 | 11,6 | 29.7 | 39,7 | — 0,86 207 | ||
27,5 | 12,5 | 29.5 | 41,3 | — 0,944 | ||
26,2 | 11,9 | 33.3 | 43,1 | — 0,82 353 | ||
26,5 | 10,6 | 41.2 | 46,5 | — 0,5 | ||
30,2 | 9,3 | 39.1 | 45,0 | — 0,63 441 | ||
32,5 | 11,1 | 33.8 | 39,4 | — 0,5045 | ||
34,5 | 14,7 | 31.5 | 36,5 | — 0,34 014 | ||
36,6 | 17,7 | 35.3 | 35,7 | — 0,0226 | ||
38,2 | 20,1 | 40.5 | 34,5 | 0,298 507 | ||
25,0 | 3,3 | 34.1 | 44,0 | — 3 | ||
24,4 | 7,0 | 32.3 | 43,9 | — 1,65 714 | ||
22,5 | 11,3 | 31.2 | 43,8 | — 1,11 504 | ||
24,3 | 10,9 | 30.8 | 43,9 | — 1,20 183 | ||
24,9 | 10,0 | 31.2 | 44,3 | — 1,31 | ||
26,9 | 9,1 | 31.3 | 45,0 | — 1,50 549 | ||
27,2 | 8,1 | 36.9 | 43,6 | — 0,82 716 | ||
26,8 | 7,3 | 42.9 | 42,8 | 0,13 699 | ||
27,8 | 9,7 | 37.5 | 42,1 | — 0,47 423 | ||
30,3 | 13,3 | 34.7 | 38,8 | — 0,30 827 | ||
33,1 | 15,7 | 35.4 | 36,1 | — 0,4 459 | ||
35,1 | 18,2 | 38.8 | 34,3 | 0,247 253 | ||
26,1 | 3,9 | 34.7 | 43,6 | — 2,28 205 | ||
24,1 | 8,7 | 33.2 | 43,4 | — 1,17 241 | ||
22,1 | 9,7 | 32.5 | 42,8 | — 1,6 186 | ||
23,8 | 9,1 | 32.6 | 42,9 | — 1,13 187 | ||
26,9 | 7,4 | 32.3 | 43,5 | — 1,51 351 | ||
28,4 | 6,6 | 32.5 | 44,5 | — 1,81 818 | ||
27,4 | 6,4 | 32.2 | 43,3 | — 1,73 438 | ||
25,4 | 8,2 | 38.1 | 43,3 | — 0,63 415 | ||
24,8 | 10,7 | 37.8 | 42,4 | — 0,42 991 | ||
27,9 | 13,7 | 36.5 | 38,6 | — 0,15 328 | ||
31,5 | 15,8 | 36.3 | 36,1 | 0,12 658 | ||
33,1 | 18,1 | 38.2 | 34,1 | 0,226 519 | ||
26,6 | 4,9 | 35.6 | 42,1 | — 1,32 653 | ||
24,8 | 8,5 | 34.2 | 41,2 | — 0,82 353 | ||
22,7 | 8,2 | 34.5 | 40,5 | — 0,73 171 | ||
23,6 | 6,6 | 34.7 | 40,4 | — 0,86 364 | ||
26,3 | 4,7 | 34.3 | 41,5 | — 1,53 191 | ||
28,4 | 7,3 | 34.1 | 42,4 | — 1,13 699 | ||
27,8 | 9,8 | 34.6 | 43,4 | — 0,89 796 | ||
26,4 | 11,0 | 34.8 | 44,5 | — 0,88 182 | ||
25,5 | 12,1 | 35.6 | 41,2 | — 0,46 281 | ||
26,3 | 15,2 | 36.3 | 37,2 | — 0,5 921 | ||
29,5 | 17,2 | 37.1 | 35,4 | 0,98 837 | ||
31,5 | 19,6 | 37.2 | 34,1 | 0,158 163 | ||
26,7 | 5,4 | 35.7 | 39,7 | — 0,74 074 | ||
25,4 | 6,1 | 35.9 | 39,1 | — 0,52 459 | ||
23,5 | 5,1 | 36.3 | 38,2 | — 0,37 255 | ||
22,7 | 5,6 | 36.7 | 37,4 | — 0,125 | ||
25,1 | 5,8 | 36.4 | 38,4 | — 0,34 483 | ||
27,1 | 8,8 | 37.4 | 39,2 | — 0,20 455 | ||
28,1 | 11,7 | 37.5 | 41,1 | — 0,30 769 | ||
29,6 | 14,9 | 37.9 | 42,1 | — 0,28 188 | ||
28,9 | 16,8 | 37.8 | 40,9 | — 0,18 452 | ||
28,3 | 18,3 | 37.4 | 37,6 | — 0,1 093 | ||
29,5 | 19,6 | 37.5 | 34,3 | 0,163 265 | ||
30,9 | 21,1 | 37.2 | 34,1 | 0,146 919 | ||
max | 42,6 | 46,5 | 0,4359 | |||
min | 22,1 | 28,1 | 29,8 | — 5,3 | ||
· Далее было произведено разделение территории на обучающую и тестируемую. В нашем случае все 96 блоков на карте являются обучающими. А значения тестируемой части карты были даны изначально. Для оценки взаимной информативности имеющиеся средние значения показателей необходимо ранжировать. Для этого все имеющиеся 96 значений поделили на 3 ранга.
Важным является то, что число значений относящихся к тому или иному рангу не должно быть меньше 10. Иначе необходимо изменять шаги ранжирования. В нашем случае этим обстоятельством осложнилось ранжирование вертикального градиента, поэтому пришлось изменить шаги ранжирования, которые оказались неравнозначными. Ранжирование показателей приведено в таблице 2.
Таблица 2. Ранжирование показателей № п/п | Показатели | Значения показателей по рангам | |||
Защищенность | [0 — 8) | [8 — 16) | [16 — 24] | ||
Суммарная защищенность | [22,1 — 28,9) | [28,9 — 35,7) | [35,7- 42,6] | ||
Уровень | [28,1 — 33,1) | [33,1 — 38,1) | [38,1 — 43,0] | ||
Уровень | [29,8 — 35,4) | [35,4 — 40,9) | [40,9 — 46,5] | ||
Вертикальный градиент | [-5,3: -1,0) | [-1,0: -0,3) | [-0,3 0,4539) | ||
Водопроводимость | < 100 | 100−200 | > 200 | ||
Тип территории | индустриал | селитебная | с/х | ||
Содержание хлоридов в | < 1 ПДК | 1−3 ПДК | > 3 ПДК | ||
Проранжировав показатели, произведем переход от средних значений показателей к рангам. Ранжированные значения показателей по блокам приведены в таблице 3.
Таблица 3. Ранжированные значения показателей по блокам № блока | Суммарная защищенность | Защищенность | Уровень | Уровень | Вертикальный градиент | Водопроводимость T | Тип территории | Содержание хлоридов в | |
Число значений каждого из показателей, попавших в тот или иной ранг, приведены в таблице 4.
Таблица 4. Число значений попавших в ранги № п/п | Показатели | Число значений показателей входящих в ранг | |||
Защищенность | |||||
Суммарная защищенность | |||||
Уровень | |||||
Уровень | |||||
Вертикальный градиент | |||||
Водопроводимость | |||||
Тип территории | |||||
Содержание хлоридов в | |||||
Аналогичные формулы для: сложных двоичных признаков: сложных троичных признаков. В нашем случае в качестве функции отклика является содержание хлоридов в водоносном горизонте. Таким образом, имея ввиду наличие 3-х рангов по хлору полная информативность равняется: начения логарифмов приведены в приложении 1. Узнав значения логарифмов получаем полную информативность: бит.
В дальнейшем при расчете взаимной информативности признаков 145,9 бит будут приняты за 100% информативности, по отношению к которой и будет выясняться наибольшая информативность признаков.
Производим расчет информативности простых признаков.
1. Суммарная защищенность.
2.
Cl | |||||
145,9−143=2,9 бит | |||||
В дальнейшем при расчете информативности в скобочках сразу будут приведены значения .
2. Защищенность.
Cl | |||||
27=39,31 | |||||
49=74,3 | |||||
20=30,8 | |||||
145,9−144,41=1,5 бит | |||||
1. Уровень касимовского горизонта
Cl | |||||
35=51,7 | |||||
49=70,7 | |||||
12=18,7 | |||||
145,9−141, 1=4,8 бит | |||||
2. Уровень подольско-мячковского горизонта
Cl | |||||
19=29,1 | |||||
40=60,9 | |||||
37=55,8 | |||||
145,9−145,8=0,1 бит | |||||
3. Вертикальный градиент
Cl | |||||
20=30,4 | |||||
39=55,2 | |||||
37=57,9 | |||||
145,9−143,5=2,4 бит | |||||
4. Водопроводимость
Cl Т | |||||
; | 16=15,3 | ||||
40=61,2 | |||||
; | 40=37,4 | ||||
145,9−113,9=32 бит | |||||
Тип территории Cl ТТ | |||||
; | 19=17,5 | ||||
42=53,1 | |||||
35=53,5 | |||||
145,9−124,1=21,8 бит | |||||
Все произведенные выше расчеты информативности простых признаков приведены в таблице 5.
Таблица 5. Информативность простых признаков Показатели | Информативность | ||
бит | % | ||
Суммарная защищенность | 2,9 | 1,9 | |
Защищенность | 1,5 | 1,0 | |
Уровень касимовского горизонта | 4,8 | 3,3 | |
Уровень подольско-мячковского горизонта | 0,1 | 0,06 | |
Вертикальный градиент | 2,4 | 1,6 | |
Водопроводимость | 21,9 | ||
Тип территории | 21,8 | 14,9 | |
Уровень концентрации хлоридов | 145,9 | ||
Анализируя полученные результаты можно сказать, что половина признаков оказались практически не информативными, и ни один из признаков не показал минимальной информативности в 30% для дальнейшего их использования для диагностики и прогнозирования экологического состояния ПТГГС.
В сложившейся ситуации необходимо сформировать сложные двоичные и троичные признаки, чтобы повысить информативность имеющихся признаков с помощью способа перебора. Для этого выбираем 4 наиболее информативных признака таких как: суммарная защищенность; уровень касимовского горизонта; водопроводимость; тип территории ().
Нам необходимо выбрать 3 наиболее информативных двоичных признака и 1 троичный.
Вначале сформируем двоичные признаки и выберем наиболее информативные из них.
Здесь и в дальнейшем будут указаны в отдельном столбце сразу значения относительной информативности.
1. Суммарная защищенность и уровень касимовского горизонта ().
Cl | Сумма | Относительная информативность | ||||
45,3 | ||||||
20,3 | ||||||
18,9 | ||||||
; | 4,9 | |||||
; | 5,5 | |||||
27,8 | ||||||
3,75 | ||||||
; | ; | |||||
; | 2,8 | |||||
сумма | 130,25 | |||||
145,9−130,25=15,7 бит | ||||||
3. Суммарная защищенность и водопроводимость ().
Cl | Сумма | Относительная информативность | ||||
; | 25,7 | |||||
; | 8,8 | |||||
23,6 | ||||||
; | 3,2 | |||||
; | 3,7 | |||||
; | 3,6 | |||||
25,3 | ||||||
; | ||||||
; | ||||||
сумма | 105,9 | |||||
145,9−105,9=40 бит | ||||||
4. Суммарная защищенность и тип территории ()
Cl | Сумма | Относительная информативность | ||||
24,2 | ||||||
18,8 | ||||||
20,1 | ||||||
; | ; | |||||
; | 3,8 | |||||
; | 3,9 | |||||
25,9 | ||||||
; | ||||||
; | 13,1 | |||||
сумма | 111,8 | |||||
145,9−111,8=34,1 бит | ||||||
5. Уровень касимовского горизонта и водопроводимость ()
Cl | Сумма | Относительная информативность | ||||
; | 21,5 | |||||
27,5 | ||||||
; | 20,1 | |||||
25,2 | ||||||
; | 9,7 | |||||
; | 4,8 | |||||
; | ; | |||||
; | 2,8 | |||||
; | ; | |||||
сумма | 104,1 | |||||
145,9−104,1=41,8 бит | ||||||
6. Уровень касимовского горизонта и тип территории ()
Cl | Сумма | Относительная информативность | ||||
23,5 | ||||||
; | ||||||
28,1 | ||||||
; | 13,8 | |||||
; | ||||||
; | 5,5 | |||||
; | 3,6 | |||||
сумма | 116,5 | |||||
145,9−116,5=29,4 бит | ||||||
7. Водопроводимость и тип территории
Cl | Сумма | Относительная информативность | ||||
; | ; | |||||
; | 25,9 | |||||
; | ||||||
30,3 | ||||||
; | ; | |||||
; | 4,3 | |||||
; | 7,6 | |||||
; | ; | |||||
; | 3,2 | |||||
сумма | 82,3 | |||||
145,9−82,3=63,6 бит | ||||||
Анализируя имеющиеся результаты расчетов, выбираем 3 двоичных признака с наибольшей информативностью:, , .
Теперь необходимо сформировать троичные признаки.
1. Суммарная защищенность, уровень касимовского горизонта и водопроводимость
Cl | Сумма | Относительная информативность | ||||
; | 15,8 | |||||
10,9 | ||||||
; | 8,3 | |||||
; | ; | |||||
; | ; | |||||
; | 7,6 | |||||
; | ||||||
10,1 | ||||||
; | ; | |||||
12,5 | ||||||
; | 3,2 | |||||
; | 3,6 | |||||
12,9 | ||||||
; | 2,8 | |||||
; | ; | |||||
; | ; | |||||
; | ||||||
; | ; | |||||
; | ; | |||||
; | ; | |||||
сумма | 91,7 | |||||
145,9−91,7=54,2 бит | ||||||
2. Уровень касимовского горизонта, водопроводимость и тип территории
Cl | Сумма | Относительная информативность | ||||
; | ; | |||||
; | 9,7 | |||||
; | 5,5 | |||||
; | 8,3 | |||||
; | 13,2 | |||||
; | ; | |||||
; | 2,8 | |||||
; | ; | |||||
; | ; | |||||
; | 12,5 | |||||
13,8 | ||||||
; | 2,8 | |||||
; | ; | |||||
; | 2,8 | |||||
; | ; | |||||
; | ||||||
; | 3,2 | |||||
; | ; | |||||
; | 4,8 | |||||
; | ||||||
сумма | 71,7 | |||||
145,9−71,7=71,7 бит | ||||||
3. Суммарная защищенность, уровень касимовского горизонта и тип территории
Cl | Сумма | Относительная информативность | ||||
9,2 | ||||||
; | 5,5 | |||||
; | 2,8 | |||||
; | 8,9 | |||||
15,9 | ||||||
; | ; | |||||
7,6 | ||||||
; | ; | |||||
; | ; | |||||
13,8 | ||||||
; | ; | |||||
; | 8,9 | |||||
; | ||||||
; | ; | |||||
; | ; | |||||
; | ; | |||||
; | ||||||
; | ; | |||||
; | ; | |||||
; | ; | |||||
сумма | 86,3 | |||||
145,9−86,3=59,6 бит | ||||||
4. Суммарная защищенность, водопроводимость и тип территории
Cl | Сумма | Относительная информативность | ||||
; | ; | |||||
; | 8,3 | |||||
; | ||||||
18,7 | ||||||
; | ; | |||||
; | ; | |||||
; | ; | |||||
; | 13,8 | |||||
; | ||||||
; | ||||||
; | ||||||
; | 7,6 | |||||
; | 3,2 | |||||
; | ; | |||||
; | ; | |||||
; | ; | |||||
; | 2,8 | |||||
; | ; | |||||
; | ; | |||||
сумма | 70,3 | |||||
145,9−70,3=75,6 бит | ||||||
В результате проведения расчетов информативности методом перебора, необходимо выбрать наиболее информативных три двоичных и один троичный признаков для того, чтобы произвести тестирование модели. Итоговые значения наиболее информативных сложных двоичных и троичных признаков сведены в таблице 6.
Таблица 6. Наиболее информативные сложные признаки
Показатели | Информативность | ||
Бит | % | ||
40,0 | 27,4 | ||
41,8 | 28,6 | ||
63,6 | 43,6 | ||
75,6 | 51,9 | ||
Уровень концентрации хлоридов | 145,9 | ||
3. Тестирование модели
После определения более информативных сложных признаков необходимо произвести тестирование сформированной модели для выяснения её пригодности для дальнейшего прогноза. Для этого используется контрольная модель, с помощью которой, с использованием матриц взаимных переходов, считается сумма голосов. В контрольной части смотрят сочетание значений признаков и по матрице взаимных переходов определяют количество попаданий этих сочетаний в 1, 2 и 3 ранг по хлору и в итоге суммируют их. Тот ранг, в который попадет наибольшее значение сочетаний, будет считаться прогнозным. Этот ранг будет сравниваться с фактическим значением ранга по хлору. Чем больше будет совпадений прогнозных значений с фактическими, тем более пригодной для прогноза окажется наша модель.
Результаты диагностирования модели приведены в таблице 7
Таблица 7. Диагностирование модели № блока | Ранги | Содержание хлора | ||||
фактическое | прогнозное | |||||
Анализируя таблицу стоит обратить внимание на так называемый «отказ», который свидетельствует об недостаточности значений показателей входящих в ранг. В нашей ситуации при нашей модели, разбитой на 96 блоков, идеальным количеством значений входящих в 3 ранга составляло бы 32. К сожалению, в нашей модели такое условие не выполняется. В некоторых случаях число значений входящих в ранг не превышает и 12. Это снижает эффективность нашей модели.
На основании составленной таблицы необходимо посчитать ошибку по формуле:, где — количество совпадений прогнозных значений и фактических.
Такая высокая ошибка свидетельствует об неэффективности сформированной нами информационной модели.
Если следовать методике диагностирования и прогноза экологического состояния ПТГГС, то после проведения тестирования имеющейся информационной модели, нам необходимо произвести прогноз состояния ПТГГС. В нашем случае, т.к. модель оказалась неэффективной, прогноз, с точки зрения здравого смысла, производить не стоит, потому что его результаты не будут корректными и достоверными. Но т.к. главной нашей целью является освоение методики проведения диагностирования и прогноза экологического состояния ПТГГС, прогноз необходимо произвести.
Результаты прогноза приведены в таблице 8.
Таблица 8. Результаты прогноза
№ блока | Ранги | Содержание хлора | |||
прогнозное | |||||
Как уже было сказано выше, данный прогноз является некорректным и не представляет никакой ценности.
Заключение
В процессе проведения данной работы нами были оценены гидрогеологические, геологические и техногенные условия изучаемой территории, выяснено влияния имеющихся условий на загрязнение подольско-мячковского горизонта хлором, была сформирована информационно-диагностическая модель, подвергавшаяся тестированию. Выяснилось, что имеющаяся модель не эффективна, потому как имеет большую ошибку. В таком случае необходимо повышать эффективность модели, теми или иными способами.
На мой взгляд, есть несколько способов повышения эффективности модели, а именно: — замена показателей, не полностью связанных с функцией отклика и не сильно влияющих на ее изменение. В нашем случае, такими признаками, на мой взгляд, являются защищенность и вертикальный градиент, т.к. он имеет важную роль в местах перетекания, а таких мест на карте не много, поэтому его роль незначительна; - введение дополнительного показателя, такого как минерализация вод подольско-мячковского горизонта; - разбиение имеющихся значений не на 3 ранга, а, предположим, на 2 или 4 ранга; - изменить способ разбиения на обучающую и контрольную части; - уточнение достоверности исходных данных; - разбиение на ранги таким образом, чтобы в каждый ранг входило одинаковое количество значений показателя.
1. Учебное пособие Природно-технические гидрогеологические системы. М.: 2006 г.
2. Белоусова А. П., Гавич И. К., Лисенков А. Б. Попов Е.В., Экологическая гидрогеология