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

курсовая работа

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

Делись добром ;)