Реферат: Диагностирование и прогноз экологического состояния природно-технической гидрогеологической системы - Refy.ru - Сайт рефератов, докладов, сочинений, дипломных и курсовых работ

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

Рефераты по геологии » Диагностирование и прогноз экологического состояния природно-технической гидрогеологической системы

Содержание

Введение

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 блоков, размер каждого составил 1x1 км, и вычислили средние значения этих показателей по блоку.

Вертикальный градиент вычислялся как разность уровней касимовского и подольско-мячковского водоносных горизонтов, отнесенная к мощности кривякинского водоупора: Подсчитанные средние значения показателей представлены в таблице 1. Таблица 1. Средние значения показателей по блокам № блока Суммарная защищенность Защищенность

Уровень

Уровень

Вертикальный градиент

1 25,6 0 33 38,3 -5,3
2 24,5 3,9 34,1 32,4 0,435897
3 28,9 7,5 32,8 31,9 0,12
4 32,5 11,4 32,1 33,2 -0,09649
5 33,7 14,8 31,1 34,2 -0,20946
6 33,1 17,0 31,4 34,7 -0,19412
7 33,0 14,6 31,2 36,3 -0,34932
8 32,2 12,6 30,1 38,9 -0,69841
9 33,3 14,4 33,3 38,1 -0,33333
10 37,0 19,1 36,6 37,1 -0,02618
11 39,0 22,5 40,5 38,4 0,093333
12 39,2 24,0 43 42 0,041667
13 25,2 0 34,3 39,1 -4,8
14 23,8 3,3 32,8 33,4 -0,18182
15 26,7 6,8 31,1 29,8 0,191176
16 32,0 10,6 29,9 31,0 -0,10377
17 37,3 14,0 28,4 32,3 -0,27857
18 38,3 15,5 28,1 33,6 -0,35484
19 34,9 13,5 29,6 37,8 -0,60741
20 32,7 11,7 29,8 40,6 -0,92308
21 36,3 13,6 31,7 40,3 -0,63235
22 39,1 18,2 35 37,1 -0,11538
23 42,5 22,3 39,1 38,1 0,044843
24 42,6 23,6 42,5 40,3 0,09322
25 25,4 1,5 33,9 41,2 -4,86667
26 23,5 4,2 32,2 33,8 -0,38095
27 25,7 6,6 30,8 32,1 -0,19697
28 31,0 10,8 29,3 34 -0,43519
29 34,9 14,3 28,1 36,1 -0,55944
30 33,6 14,2 32.8 39 -0,43662
31 32,5 12,5 36.2 42,4 -0,496
32 31,4 10,5 33.9 43,7 -0,93333
33 34,4 12,4 31.6 39,8 -0,66129
34 37,9 16,9 33.3 36,4 -0,18343
35 40,1 20,1 37.1 36,4 0,034826
36 41,1 22,3 40.6 36,1 0,201794
37 25,1 2,3 33.3 42,1 -3,82609
38 23,2 5,5 32.1 39,1 -1,27273
39 25,0 9,9 31.4 38,4 -0,70707
40 28,5 11,6 29.7 39,7 -0,86207
41 27,5 12,5 29.5 41,3 -0,944
42 26,2 11,9 33.3 43,1 -0,82353
43 26,5 10,6 41.2 46,5 -0,5
44 30,2 9,3 39.1 45,0 -0,63441
45 32,5 11,1 33.8 39,4 -0,5045
46 34,5 14,7 31.5 36,5 -0,34014
47 36,6 17,7 35.3 35,7 -0,0226
48 38,2 20,1 40.5 34,5 0,298507
49 25,0 3,3 34.1 44,0 -3
50 24,4 7,0 32.3 43,9 -1,65714
51 22,5 11,3 31.2 43,8 -1,11504
52 24,3 10,9 30.8 43,9 -1,20183
53 24,9 10,0 31.2 44,3 -1,31
54 26,9 9,1 31.3 45,0 -1,50549
55 27,2 8,1 36.9 43,6 -0,82716
56 26,8 7,3 42.9 42,8 0,013699
57 27,8 9,7 37.5 42,1 -0,47423
58 30,3 13,3 34.7 38,8 -0,30827
59 33,1 15,7 35.4 36,1 -0,04459
60 35,1 18,2 38.8 34,3 0,247253
61 26,1 3,9 34.7 43,6 -2,28205
62 24,1 8,7 33.2 43,4 -1,17241
63 22,1 9,7 32.5 42,8 -1,06186
64 23,8 9,1 32.6 42,9 -1,13187
65 26,9 7,4 32.3 43,5 -1,51351
66 28,4 6,6 32.5 44,5 -1,81818
67 27,4 6,4 32.2 43,3 -1,73438
68 25,4 8,2 38.1 43,3 -0,63415
69 24,8 10,7 37.8 42,4 -0,42991
70 27,9 13,7 36.5 38,6 -0,15328
71 31,5 15,8 36.3 36,1 0,012658
72 33,1 18,1 38.2 34,1 0,226519
73 26,6 4,9 35.6 42,1 -1,32653
74 24,8 8,5 34.2 41,2 -0,82353
75 22,7 8,2 34.5 40,5 -0,73171
76 23,6 6,6 34.7 40,4 -0,86364
77 26,3 4,7 34.3 41,5 -1,53191
78 28,4 7,3 34.1 42,4 -1,13699
79 27,8 9,8 34.6 43,4 -0,89796
80 26,4 11,0 34.8 44,5 -0,88182
81 25,5 12,1 35.6 41,2 -0,46281
82 26,3 15,2 36.3 37,2 -0,05921
83 29,5 17,2 37.1 35,4 0,098837
84 31,5 19,6 37.2 34,1 0,158163
85 26,7 5,4 35.7 39,7 -0,74074
86 25,4 6,1 35.9 39,1 -0,52459
87 23,5 5,1 36.3 38,2 -0,37255
88 22,7 5,6 36.7 37,4 -0,125
89 25,1 5,8 36.4 38,4 -0,34483
90 27,1 8,8 37.4 39,2 -0,20455
91 28,1 11,7 37.5 41,1 -0,30769
92 29,6 14,9 37.9 42,1 -0,28188
93 28,9 16,8 37.8 40,9 -0,18452
94 28,3 18,3 37.4 37,6 -0,01093
95 29,5 19,6 37.5 34,3 0,163265
96 30,9 21,1 37.2 34,1 0,146919
max 42,6 24 43 46,5 0,4359
min 22,1 1 28,1 29,8 -5,3

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

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

 Таблица 2. Ранжирование показателей № п/п

Показатели Значения показателей по рангам
1 2 3
1 Защищенность [0 – 8) [8 – 16) [16 – 24]
2 Суммарная защищенность [22,1 – 28,9) [28,9 – 35,7) [35,7- 42,6]
3

Уровень

[28,1 – 33,1) [33,1 – 38,1) [38,1 – 43,0]
4

Уровень

[29,8 - 35,4) [35,4 – 40,9) [40,9 – 46,5]
5 Вертикальный градиент [-5,3: -1,0) [-1,0: -0,3) [-0,3 0,4539)
6 Водопроводимость < 100 100-200 > 200
7 Тип территории индустриал селитебная с/х
8

Содержание хлоридов в

< 1 ПДК 1-3 ПДК > 3 ПДК

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

Таблица 3. Ранжированные значения показателей по блокам № блока Суммарная защищенность Защищенность

Уровень

Уровень

Вертикальный градиент

Водопроводимость T

Тип территории

Содержание хлоридов в

1 1 1 2 2 1 3 3

1

2 1 1 2 1 3 3 3

1

3 2 1 1 1 3 3 2

2

4 2 2 1 1 3 2 2

2

5 2 2 1 1 3 2 3

3

6 2 3 1 1 3 2 3

2

7 2 2 1 2 2 2 2

1

8 2 2 1 2 2 3 2

1

9 2 2 2 2 2 3 2

1

10 3 3 2 2 3 3 1

1

11 3 3 3 2 3 2 2

2

12 3 3 3 3 3 1 3

3

13 1 1 2 2 1 3 3

1

14 1 1 1 1 3 3 3

1

15 1 1 1 1 3 3 2

2

16 2 2 1 1 3 2 2

2

17 3 2 1 1 3 2 3

3

18 3 2 1 1 2 2 3

2

19 2 2 1 2 2 2 2

1

20 2 2 1 2 2 3 2

1

21 3 2 1 2 2 3 2

1

22 3 3 2 2 3 3 1

1

23 3 3 3 2 3 2 2

2

24 3 3 3 2 3 1 3

3

25 1 1 2 3 1 3 2

2

26 1 1 1 1 2 3 3

1

27 1 1 1 1 3 3 3

1

28 2 2 1 1 2 2 2

2

29 2 2 1 2 2 2 3

3

30 2 2 1 2 2 2 3

2

31 2 2 2 3 2 2 2

1

32 2 2 2 3 2 3 2

2

33 2 2 1 2 2 3 2

1

34 3 3 2 2 3 3 2

2

35 3 3 2 2 3 2 2

2

36 3 3 3 2 3 1 3

3

37 1 1 2 3 1 3 3

1

38 1 1 1 2 1 3 2

2

39 1 2 1 2 2 3 2

2

40 1 2 1 2 2 2 2

2

41 1 2 1 3 2 2 3

3

42 1 2 2 3 2 2 3

2

43 1 2 3 3 2 2 2

1

44 2 2 3 3 2 3 2

2

45 2 2 2 2 2 3 2

1

46 2 2 1 2 2 3 1

1

47 3 3 2 2 3 2 2

2

48 3 3 3 1 3 1 3

3

49 1 1 2 3 1 3 2

2

50 1 1 1 3 1 2 1

2

51 1 2 1 3 1 2 1

1

52 1 2 1 3 1 3 2

1

53 1 2 1 3 2 2 3

3

54 1 2 1 3 1 1 3

3

55 1 2 2 3 2 3 1

1

56 1 1 3 3 3 3 2

1

57 1 2 2 3 2 1 1

2

58 2 2 2 2 2 1 2

3

59 2 2 2 2 3 2 3

2

60 2 3 3 1 3 2 3

1

61 1 1 2 3 1 3 2

2

62 1 2 2 3 1 2 1

2

63 1 2 1 3 1 2 1

1

64 1 2 1 3 1 3 2

1

65 1 1 1 3 1 2 3

3

66 1 1 1 3 1 1 3

3

67 1 1 1 3 1 3 1

1

68 1 2 3 3 2 3 2

1

69 1 2 2 3 2 1 1

2

70 1 2 2 2 3 1 2

3

71 2 2 2 2 3 2 3

2

72 2 3 3 1 3 2 3

1

73 1 1 2 3 1 3 2

2

74 1 2 2 3 2 2 1

1

75 1 2 2 2 2 2 1

1

76 1 1 2 2 2 3 2

2

77 1 1 2 3 1 2 3

3

78 1 1 2 3 1 1 3

2

79 1 2 2 3 2 3 1

1

80 1 2 2 3 2 3 2

2

81 1 2 2 3 2 1 1

2

82 1 2 2 2 3 1 2

3

83 2 3 2 2 3 2 3

2

84 2 3 2 1 3 2 3

1

85 1 1 2 2 2 3 2

2

86 1 1 2 2 2 2 1

2

87 1 1 2 2 2 2 1

1

88 1 1 2 2 3 3 2

1

89 1 1 2 2 2 2 3

3

90 1 2 2 2 3 1 3

3

91 1 2 2 3 2 3 1

1

92 2 2 2 3 3 3 2

1

93 2 3 2 3 3 1 1

2

94 1 3 2 2 3 1 2

2

95 2 3 2 1 3 2 3

3

96

2

3

2

1

3

2

3

1

Число значений каждого из показателей, попавших в тот или иной ранг, приведены в таблице 4.

 Таблица 4. Число значений попавших в ранги № п/п  Показатели Число значений показателей входящих в ранг
1 2 3
1 Защищенность 27 49 20
2 Суммарная защищенность 22 30 14
3

Уровень

35 49 12
4

Уровень

19 40 37
5 Вертикальный градиент 20 39 37
6 Водопроводимость 16 40 40
7 Тип территории 19 42 35
8

Содержание хлоридов в

39 36 21

Аналогичные формулы для: сложных двоичных признаков: сложных троичных признаков. В нашем случае в качестве функции отклика является содержание хлоридов в водоносном горизонте. Таким образом, имея ввиду наличие 3-х рангов по хлору полная информативность равняется: начения логарифмов приведены в приложении 1. Узнав значения логарифмов получаем полную информативность:  бит.

В дальнейшем при расчете взаимной информативности признаков 145,9 бит будут приняты за 100% информативности, по отношению к которой и будет выясняться наибольшая информативность признаков.

Производим расчет информативности простых признаков.

1.  Суммарная защищенность.

2. 

 Cl

1 2 3
1 22 20 10

 2 14 12 4

3 3 6 5

39 38 19 145,9-143=2,9 бит

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

2. Защищенность.

 Cl

1 2 3
1 11 12 4 27=39,31
 2 22 17 10 49=74,3
3 6 9 5 20=30,8

39 38 19 145,9-144,41=1,5 бит

1.  Уровень касимовского горизонта

 Cl

1 2 3
1 15 12 8 35=51,7
 2 19 23 7 49=70,7
3 5 3 4 12=18,7

39 38 19 145,9-141, 1=4,8 бит

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

 Cl

1 2 3
1 8 7 4 19=29,1
 2 16 16 8 40=60,9
3 15 15 7 37=55,8

39 38 19 145,9-145,8=0,1 бит

3.  Вертикальный градиент

 Cl

1 2 3
1 8 8 4 20=30,4
 2 19 15 5 39=55,2
3 12 15 10 37=57,9

39 38 19 145,9-143,5=2,4 бит

4.  Водопроводимость

 Cl Т

1 2 3
1 - 6 10 16=15,3
 2 13 18 9 40=61,2
3 26 14 - 40=37,4

39 38 19 145,9-113,9=32 бит

 Тип территории Cl ТТ

1 2 3
1 12 7 - 19=17,5
 2 16 23 3 42=53,1
3 11 8 16 35=53,5

39 38 19 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
Водопроводимость 32 21,9
Тип территории 21,8 14,9

Уровень концентрации хлоридов

145,9

100

Анализируя полученные результаты можно сказать, что половина признаков оказались практически не информативными, и ни один из признаков не показал минимальной информативности в 30% для дальнейшего их использования для диагностики и прогнозирования экологического состояния ПТГГС.

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

Нам необходимо выбрать 3 наиболее информативных двоичных признака и 1 троичный.

Вначале сформируем двоичные признаки и выберем наиболее информативные из них.

Здесь и в дальнейшем будут указаны в отдельном столбце сразу значения относительной информативности.

1. Суммарная защищенность и уровень касимовского горизонта ().

 Cl

1 2 3 Сумма

Относительная информативность

12 11 15 5 31 45,3
 21 6 6 2 14 20,3
22 6 5 2 13 18,9
32 2 3 - 5 4,9
33 - 2 4 6 5,5
11 8 5 5 18 27,8
31 1 1 1 3 3,75
13 3 - - 3 0
23 2 1 - 3 2,8
сумма 96 130,25
145,9-130,25=15,7 бит

3.   Суммарная защищенность и водопроводимость ().

 Cl

1 2 3 Сумма

Относительная информативность

13 17 10 - 27 25,7
23 7 3 - 10 8,8
22 5 8 3 16 23,6
33 3 1 - 4 3,2
32 - 5 1 6 3,7
31 1 - 4 5 3,6
12 6 5 5 16 25,3
11 - 5 5 10 10
21 - 1 1 2 2
сумма 96 105,9
145,9-105,9=40 бит

4.  Суммарная защищенность и тип территории ()

Cl

1 2 3 Сумма

Относительная информативность

13 7 2 8 17 24,2
22 8 6 1 15 18,8
23 5 5 3 13 20,1
31 2 - - 2 0
32 1 5 - 6 3,8
33 - 1 5 6 3,9
12 6 12 2 20 25,9
21 1 1 - 2 2
11 9 6 - 15 13,1
сумма 96 111,8
145,9-111,8=34,1 бит

5.  Уровень касимовского горизонта и водопроводимость ()

 Cl

1 2 3 Сумма

Относительная информативность

23 13 9 - 22 21,5
12 4 8 6 18 27,5
13 11 4 - 15 20,1
22 6 8 3 17 25,2
21 - 6 4 10 9,7
32 3 2 - 5 4,8
31 - - 4 4 0
33 2 1 - 3 2,8
11 - - 2 2 0
сумма 96 104,1
145,9-104,1=41,8 бит

6.  Уровень касимовского горизонта и тип территории ()

 Cl

1 2 3 Сумма

Относительная информативность

23 6 5 4 15 23,5
12 8 8 - 16 16
13 3 3 8 14 20
22 6 11 3 20 28,1
21 8 6 - 14 13,8
32 3 3 - 6 6
33 2 - 4 6 5,5
11 4 1 - 5 3,6
сумма 96 116,5
145,9-116,5=29,4 бит

7.  Водопроводимость и тип территории

 Cl

1 2 3 Сумма

Относительная информативность

33 7 - - 7 0
32 12 14 - 26 25,9
22 4 8 - 12 11
23 4 7 9 20 30,3
31 7 - - 14 0
13 - 1 7 8 4,3
21 5 3 - 8 7,6
11 - 4 - 4 0
12 - 1 3 4 3,2
сумма 96 82,3
145,9-82,3=63,6 бит

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

Теперь необходимо сформировать троичные признаки.

1. Суммарная защищенность, уровень касимовского горизонта и водопроводимость

 Cl

1 2 3 Сумма

Относительная информативность

123 9 7 - 16 15,8
112 2 2 3 7 10,9
113 6 3 - 9 8,3
111 - - 2 2 0
133 2 - - 2 0
121 - 5 3 8 7,6
221 - 1 1 2 2
222 3 3 1 7 10,1
232 2 - - 2 0
122 3 3 2 8 12,5
223 3 1 - 4 3,2
213 4 1 - 5 3,6
212 2 5 2 9 12,9
323 2 1 - 3 2,8
332 - 2 - 2 0
331 - - 4 4 0
312 - 1 1 2 2
313 1 - - 1 0
322 - 2 - 2 0
233 - 1 - 1 0
сумма 96 91,7
145,9-91,7=54,2 бит

2. Уровень касимовского горизонта, водопроводимость и тип территории

 Cl

1 2 3 Сумма

Относительная информативность

233 4 - - 4 0
132 6 4 - 10 9,7
122 2 4 - 6 5,5
123 - 3 6 9 8,3
232 5 9 - 14 13,2
231 4 - - 4 0
322 1 2 - 3 2,8
313 - - 4 4 0
133 3 - - 3 0
222 1 2 - 8 12,5
223 2 4 3 9 13,8
332 2 1 - 3 2,8
131 2 - - 2 0
121 2 1 - 3 2,8
113 - - 2 2 0
211 - 4 1 4 0
212 - 1 3 4 3,2
323 2 - - 2 0
221 3 2 - 5 4,8
213 - 1 1 2 2
сумма 96 71,7
145,9-71,7=71,7 бит

3. Суммарная защищенность, уровень касимовского горизонта и тип территории

 Cl

1 2 3 Сумма

Относительная информативность

122 1 8 1 10 9,2
112 2 4 - 6 5,5
111 2 1 - 3 2,8
113 4 - 5 9 8,9
121 6 5 1 12 15,9
132 3 - - 3 0
223 2 2 1 5 7,6
233 2 - - 2 0
232 - 1 - 1 0
123 4 2 3 9 13,8
221 - 1 - 1 0
212 5 4 - 9 8,9
213 - 2 2 4 4
321 2 - - 2 0
332 - 2 - 2 0
333 - - 4 4 0
313 - 1 1 2 2
312 1 - - 1 0
322 - 2 - 2 0
211 1 - - 1 0
сумма 96 86,3
145,9-86,3=59,6 бит

4. Суммарная защищенность, водопроводимость и тип территории

Cl

1 2 3 Сумма

Относительная информативность

133 7 - - 7 0
232 6 3 - 9 8,3
222 3 3 - 6 6
223 4 5 3 12 18,7
331 2 - - 2 0
322 - 4 - 4 0
313 - - 4 4 0
132 5 10 - 15 13,8
323 - 1 1 2 2
332 1 1 - 2 2
122 1 1 - 2 2
121 5 3 - 8 7,6
113 - 1 3 4 3,2
131 4 - - 4 0
111 - 3 - 3 0
212 - - 1 1 0
112 - 1 2 3 2,8
211 - 1 - 1 0
231 1 - - 1 0
сумма 96 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

100

3. Тестирование модели

После определения более информативных сложных признаков необходимо произвести тестирование сформированной модели для выяснения её пригодности для дальнейшего прогноза. Для этого используется контрольная модель, с помощью которой, с использованием матриц взаимных переходов, считается сумма голосов. В контрольной части смотрят сочетание значений признаков и по матрице взаимных переходов определяют количество попаданий этих сочетаний в 1, 2 и 3 ранг по хлору и в итоге суммируют их. Тот ранг, в который попадет наибольшее значение сочетаний, будет считаться прогнозным. Этот ранг будет сравниваться с фактическим значением ранга по хлору. Чем больше будет совпадений прогнозных значений с фактическими, тем более пригодной для прогноза окажется наша модель.

Результаты диагностирования модели приведены в таблице 7

Таблица 7. Диагностирование модели № блока Ранги Содержание хлора
1 2 3 фактическое прогнозное
97

2 1
98

1 2
99

3 1
100

2 3
101

2 2
102

2 3
103

3 2
104

2 2
105

2 1
106

1 3
107

2 2
108

1 3
109

2 1
110

1 2
111

3 1
112

2 3
113

2 2
114

2 3
115

3 1
116

2 2
117

2 1
118

1 3
119

2 2
120

1 3

Анализируя таблицу стоит обратить внимание на так называемый «отказ», который свидетельствует об недостаточности значений показателей входящих в ранг. В нашей ситуации при нашей модели, разбитой на 96 блоков, идеальным количеством значений входящих в 3 ранга составляло бы 32.. К сожалению, в нашей модели такое условие не выполняется. В некоторых случаях число значений входящих в ранг не превышает и 12. Это снижает эффективность нашей модели.

На основании составленной таблицы необходимо посчитать ошибку по формуле: , где - количество совпадений прогнозных значений и фактических.

Такая высокая ошибка свидетельствует об неэффективности сформированной нами информационной модели.

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

Результаты прогноза приведены в таблице 8.

Таблица 8. Результаты прогноза

№ блока Ранги Содержание хлора
1 2 3 прогнозное
121

1
122

2
123

1
124

3
125

2
126

3
127

1
128

2
129

1
130

3
131

2
132

3

Как уже было сказано выше, данный прогноз является некорректным и не представляет никакой ценности.


Заключение

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

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


Список литературы

1.  Учебное пособие Природно-технические гидрогеологические системы. М.: 2006 г.

2.  Белоусова А.П., Гавич И.К., Лисенков А.Б. Попов Е.В., Экологическая гидрогеология