Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Шкинкис, Ц. Н. Проблемы гидрологии дренажа

.pdf
Скачиваний:
5
Добавлен:
22.10.2023
Размер:
13.48 Mб
Скачать

Таблица 18

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

 

 

Глубина

Средняя

Скорость

Средняя движущая

Почва

фронта

влажность

инфильтрации,

сила, создаваемая

смачивания,

переходной

градиентом

 

 

зоны,

см/с

потенциала,

 

 

см

%

дин/г

Песчаный

сугли-

ю л

31,6

23-10-4

5,48-103

нок

 

20.5

31.0

18,3-10-4

4,40-103

 

 

30.6

30.0

12,7-10-4

3,52-103

Пылеватый

су-

15,1

36,5

3,12-10-4

4,16-103

ГЛИНОК

 

19,3

36,3

2,59-10-4

3,92-103

 

 

29,6

35,7

1,92-10-4

2,87-103

По А. Н. Костикову, гидравлический уклон i на глубине I от по­ верхности определяется по формуле

•_ е + I Ч- Я) —L

/оо\

где е — высота слоя воды на поверхности почвы; I — толщина слоя почвы, в которой происходит инфильтрация; Н0— капиллярный на­ пор; L — избыточное давление воздуха в почве (сверх атмосфер­ ного) .

В формуле (33) члены е, I, Н0 соответствуют членам рgh, рgl и

2 — формулы (31). Таким образом, Костяков вводит в уравнение

Пуазейля

дополнительный член — L — противодавление воздуха.

А. Н.

Костяков отмечает также, что скорость инфильтрации

со временем приближается к коэффициенту фильтрации: v-*-k, т. е.

i >-1,0.

На самом же деле в природе гидравлический уклон убы­

вает более резко, чем это следует из формулы (33)

Из приведенных

зависимостей совершенно ясно, что на полях

мелкого

дренажа,

где путь просачивания воды от поверхности

земли до поверхности грунтовой воды значительно короче, чем на полях глубокого дренажа (H i— Я2), инфильтрационные потоки быстрее достигнут этой поверхности и начнется повышение зеркала грунтовой воды. С повышением уровня грунтовой воды путь инфильтрационного потока все больше сокращается и, следовательно, скорость подъема уровня все более возрастает.

Кривая депрессии на поле мелкого дренажа особенно быстро поднимается при наличии повторных паводков, так как влажность1

1 Гидравлический уклон уменьшается быстрее, чем следует из формулы (33), главным образом потому, что по мере смыкания менисков просачивающейся воды с менисками почвенной влаги уменьшается капиллярный напор, т. е. член Я 0 очень быстро стремится к нулю. (Прим, ред.)

50

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

С увеличением ординат у кривой депрессии

возрастает пьезо-

„ .

Аг/2 Аг/1

метрический уклон этой кривой вблизи дрен, т. е. —— > —— (см.

рис. 12). Принимая, что движение грунтовых вод на дренированном поле является потенциальным, потенциал же скорости удовлетво­ ряет уравнению Лапласа

V2? = 0

(34)

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

Согласно формуле (17), при увеличении напора и, следова­ тельно, пьезометрического уклона кривой депрессии вблизи дрен удельный дренажный сток должен возрастать. Если в начале па-

водкового цикла при одинаковом пьезометрическом уклоне Ау

на расстоянии х от дрены величина дренажного стока на полях мелкого и глубокого дренажа одинаковая, т. е. q' — q" (рис. 12), то после некоторого промежутка времени просачивания талых вод (АТ) вследствие подъема кривой депрессии на величину Д при мелком дренаже сток будет увеличиваться на Дq'. В результате q' + A q '^q " . Это и является основной причиной резкого увеличения стока в системах мелкого дренажа в начале паводкового цикла.

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

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

Продолжительность дренажного стока. По продолжительности дренажного стока Тс можно судить о роли дренажа при регулиро­ вании водного режима почв. Величина Тс в первую очередь зависит от конкретных гидрометеорологических условий данного года и от года к году меняется в широких пределах. Так, за период с 1957 по

4*

51

1971 г. в средней части Латвийской ССР (Кокнесский опытный уча­ сток) при среднеглубоком дренаже (7= 1,2 м и £ = 20 м) продолжи­ тельность стока менялась от 48 до 223 суток в году. Средняя за 14 лет Тс здесь составляет 118 суток, или 32% времени (табл. 19). Наименьшая средняя продолжительность стока наблюдается в июне, наибольшая — в апреле.

q л/(с-га)

Рис. 13. Кривые обеспечения средних суточных модулей дре­ нажного стока для среднеглубокого (1,2 м) дренажа, имеющего различные расстояния между дренами Е. 1961-62 г. Кокнесе.

1 — £=14 м; 2 — Е - 30 м.

Расстояние между дренами Е на продолжительность стока су­ щественного влияния не оказывает. Однако в большинстве случаев наблюдается тенденция к уменьшению Тс с увеличением Е (см.

таб л .19).

Продолжительность дренажного стока Тс значительно возра­ стает при увеличении глубины закладки дрен t. Так, на Кокнесском опытном участке при £ = 14 м с увеличением t от 0,9 до 1,2 м вели­ чина Тс возросла на 23%• На Кандавском участке при £ = 20 м увеличение t от 1,2 до 1,5 м дало прирост Тс на 44% (табл. 20). На Римейкском опытном участке при глубоком систематическом дренаже продолжительность стока в среднем в 1,7 раза больше, чем при соответствующем мелком дренаже (табл. 20). Еще боль­ шая разница в продолжительности стока получается при глубоком

52

Таблица 19

Средняя продолжительность дренажного стока (в сутках и °/о) на легких минеральных почвах при глубине закладки дрен 1,2 м

Расстояние

 

 

 

 

 

 

 

 

 

 

 

 

Всего

между

X

XI

XII

I

II

in

IV

V

VI

VII

 

IX

VIII

за год

дренами,

м

 

 

 

 

 

 

 

 

 

 

 

 

(X -IX )

 

 

Кокнесский опытный участок (1957—1971 гг.)

 

 

14

8

20

20

9

9

14

27

14

1

2

5

5

134

 

2,2

5,5

5,5

2,5

2,5

3,8

7,3

3,8

0,3

0,5

1,4

1,4

36,7

20

8

17

17

8

7

13

25

12

_

1

5

5

118

 

2,2

4,6

4,6

2,2

1,9

3,6

6,8

3,3

0,3

1,4

1,4

32,3

 

 

Кандавский опытный участок (1968—1971 гг.)

 

 

14

20

26

31

28

22

25

28

31

22

10

10

14

267

 

5,5

7,1

8,5

7,7

6,0

7,0

7,7

8,5

6,0

2,7

2,7

3,8

73,2

40

11

22

26

22

20

20

27

24

8

2

5

9

196

 

3,0

6,0

7,1

6,0

5,5

5,5

7,4

6,6

2,2

0,5

1,4

2,5

53,7

Таблица 20

Продолжительность дренажного стока (сутки) в зависимости от глубины закладки дрен t

Расстояние

Глубина

между

закладки

дренами,

дрен,

м

м

 

 

 

Гидрологические годы

 

 

 

1958-59

1959-60

1960-61

1961-62

1962-63

1963-64

1964-65

1965-66

1966-67

1967-68

В среднем за 10 лет

сутки %

 

 

 

Римейкский опытный участок

 

 

 

20

1,5

268

225

365

365

242

232

219

319

303

285

282

20

0,9

129

101

257

267

ПО

83

146

162

199

162

162

36

1,5

292

314

365

365

277

280

243

365

349

318

317

36

0,9

68

49

160

200

109

32

69

93

157

104

104

 

 

 

Кандавский опытный участок

 

 

 

20

1,5

290

257

349

365

365

284

292

239

157

274

273

20

1,2

228

77

191

297

138

168

150

93

143

257

189

и мелком разреженном дренаже (£ = 36 м). В отдельные многовод­ ные годы в системах глубокого дренажа сток происходит непре­ рывно. Это свидетельствует о большом значении глубокого дре­ нажа при регулировании водного режима легких суглинистых и супесчаных почв.

Обеспеченность дренажного стока. При выяснении влияния сте­ пени дренирования (£ и t) на режим дренажного стока и определе­ нии расчетного дренажного стока для каждого варианта дренажа составлялись кривые обеспеченности средних суточных модулей дренажного стока (см. рис. 13). Из этих кривых, а также табл. 21 и 22 следует, что величина модулей дренажного стока q равной обеспеченности значительно возрастает при уменьшении расстоя­ ния между дренами £. Так, в 1961-62 г. при обеспеченности Р = = 2,0% в варианте частого систематического дренажа (£=12 м) ве-

Таблица 21

Модули дренажного стока (л /(с-га )) при различных значениях обеспеченности в зависимости от расстояния между дренами. 1961-62 г. Римейкас

Расстояние

 

 

 

 

Обеспеченность Р%

 

 

 

 

дренами,

0,3

0,6

1,0

2,0

3,0

5,0

7,0

10,0

15,0

20,0

м

12

5,00

3,05

2,32

1,70

1,29

0,90

0,80

0,61

0,42

0,32

20

2,90

2,06

1,58

1,15

0,96

0,75

0,64

0,53

0,38

0,27

36

1,44

1,22

1,04

0,77

0,65

0,52

0,45

0,35

0,28

0,14

Таблица 22

 

Влияние глубины закладки дрен t и расстояний между дренами Е

 

на обеспеченность модулей дренажного стока

(л /(с-га )).

1966-67

г. Кандава

 

 

 

 

 

Обеспеченность Р%

 

 

 

Е м

t М

0,3

0,6

1,0

2,0

3,0

5,0

7,0

10,0

 

 

20

1,5

0,68

0,61

0,57

0,48

0,41

0,30

0,24

0,19

20

0,9

1,05

0,91

0,78

0,58

0,43

0,15

0,09

0,06

14

1,2

1,14

1,22

1,07

0,80

0,55

0,35

0,25

0,18

20

1,2

0,99

0,91

0,84

0,73

0,58

0,35

0,26

0,20

30

1,2

0,74

0,64

0,55

0,42

0,34

0,24

0,18

0,13

40

1,2

0,58

0,48

0,44

0,36

0,26

0,19

0,15

0,10

личина

9=1,70 л /(с-га), а

в варианте

разреженного

дренажа

(£ = 36

м) 9 = 0,77 л/(с-га).

Еще большая

разница в

величине

имеет место при £ < 2,0% . Это значит, что

расстояние между дре­

нами наиболее сильно влияет на величину

наибольшего дренаж­

ного стока. По мере увеличения значений £

разница в стоке между

частым и разреженным дренажем сокращается.

В этом можно убедиться из графика связи между равнообеспе­ ченными модулями дренажного стока и густотой дренажной сети

(рис. 14).

Благодаря большей мощности зоны аэрации и соответственно большей водовместимости почвы пик дренажного паводкового стока в системах глубокого дренажа менее высокий, чем в системах мел­

54

кого дренажа. Поэтому в первом случае, по сравнению со вторым, модули дренажного стока при малых значениях обеспеченности (Р < 1,0%) обычно меньше. Исключение бывает только тогда, когда почва перед большими паводками имеет влажность, близкую к наи­ меньшей (полевой) влагоемкости. В этом случае дренажный сток при мелком дренаже может быть меньше, чем при глубоком, неза­ висимо от величины обеспеченности (рис. 15).

При обеспеченности Р > 10% глубокий дренаж, как правило, дает значительно больший сток, чем мелкий (табл. 23).

Таблица 23

Модули дренажного стока (л/с -га)) в зависимости от глубины закладки дрен —20 м). Римейкскии опытный участок

Гидроло­

Глубина

 

 

 

Обеспеченность Р%

 

 

 

закладки

 

 

 

 

 

 

 

 

 

 

гический

 

 

 

 

 

 

 

 

 

 

дрен,

0,3

0,6

1,0

2,0

3,0

5,0

7,0

10,0

15,0

20,0

год

м

1958-59

0,9

1,19

0,96

0,78

0,50

0,35 0,22 0,15

0,09

 

1,5

0,80

0,74

0,69 0,58 0,46 0,32 0,27

0,22

1959-60

0,9

3,20

2,28

1,06

0,21 0,12 0,10 0,07

0,04

_

_

 

1,5

0,50

0,40

0,32 0,21 0,17 0,13 0,10

0,08

1960-61

0,9

1,30

0,76

0,62

0,48

0,40 0,30

0,25

0,21

0,16

0,12

 

1,5

1,00

0,57

0,47

0,37

0,32

0,27

0,24

0,21

0,18

0,15

1961-62

0,9

3,20

2,25

1,60

1,12

0,91

0,69

0,58

0,47

0,32

0,21

 

1,5

2,00

1,68

1,23

0,90

0,76

0,62

0,53

0,43

0,33

0,24

1962-63

0,9

0,70

0,59

0,50

0,37

0,30

0,20

0,16

0,10

0,06

0,04

 

1,5

0,62

0,59

0,52

0,35

0,28

0,23

0,20

0,15

0,14

0,12

1963-64

0,9

1,28

0,86

0,56

0,20

0,09

0,04

0,01

0,0

0,0

0,0

 

1,5

0,94

0,82

0,63

0,38

0,29

0,21

0,16

0,13

0,09

0,07

1964-65

0,9

1,26

0,94

0,78

0,57

0,43

0,26

0,19

0,13

0,07

0,04

 

1,5

1,44

1,00

0,73

0,46

0,37

0,28

0,24

0,18

0,14

0,10

1965-66

0,9

1,55

0,80

0,52

0,36

0,30

0,23

0,18

0,12

0,06

0,03

 

1,5

2,34

1,46

0,97

0,66

0,50

0,31

0,23

0,16

0,12

0,10

1966-67

0,9

2,00

1,67

1,33

0,80

0,61

0,44

0,33

0,23

0,15

0,09

 

1,5

1,69

1,53

1,35

0,86

0,57

0,43

0,35 0,27

0,20

0,15

1967-68

0,9

1,49

1,10

0,82

0,54

0,43

0,30

0,23

0,15

0,06

0,00

 

1,5

1,42

1,17

0,95

0,67

0,54

0,40

0,34 0,37

0,18

0,12

Среднее

0,9

1,72

1,22

0,86

0,52

0,39

0,28

0,22

0,15

0,11

0,06

 

1,5

1,27

0,99

0,79

0,54

0,43

0,32

0,27

0,22

0,17

0,13

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

55

q л/(с-га)

Рис. 14. Связь между равнообеспечен­

Рис. 15. Связь между годовым дренажным стоком

Q для средне-

ными модулями дренажного

стока

и

глубокого дренажа (1,2 м), имеющего различные расстояния между

густотой дренажной сети |.

1961-62

г.

дренами

Е. 1957— 1968 гг. Кокнесе.

 

Римейкас.

 

 

 

 

 

 

 

 

Корреляционная связь между годовым дренажным стоком для

мелкого

(0,9

м)

и

глубокого

(1,5

м)

систематического дренажа

(£ = 20

м)

Римейкского опытного участка

видна

на рис.

16. При­

мерно такая же корреля­

 

 

 

 

 

 

 

 

 

ционная связь существует

 

 

 

 

 

 

 

 

 

между годовым дренаж­ 600Г

 

 

 

 

 

 

 

 

ным

стоком

для

мелкого

 

 

 

 

 

 

 

 

 

и глубокого

разреженно­

 

 

 

 

 

 

 

 

 

го дренажа

при £ = 36

м.

500

 

 

 

 

 

> / \ /

 

Установлено,

что

для

 

 

 

 

 

 

 

дренажных

систем,

име­

 

 

 

 

 

<

f / A

5

 

ющих различную степень

 

 

 

 

 

V

/ &

 

 

Ш

 

 

 

 

Л

 

 

 

дренирования

и /),

но

 

 

 

 

 

 

 

 

одинаковые

почвенные

и

 

 

 

 

 

Л

 

 

 

гидрологические

условия,

 

 

 

 

 

i

 

 

 

300

 

 

 

 

 

 

 

 

связь между годовым сто­

 

 

 

У

 

 

 

 

ком тесная,

коэффициент

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

корреляции г>0,9.

 

 

 

 

 

 

 

 

 

 

 

Хорошая корреляцион­

200

 

;

/

 

 

 

 

 

ная связь существует так­

 

 

V

 

 

 

 

 

 

же между средними су­

 

 

 

 

 

 

 

 

 

точными

модулями

дре­

100.

/

о----------

300

400 С£=й9ммм

нажного стока для дре­

 

 

100

 

200

нажных систем,

имеющих

Рис.

16. Связь между годовым дренажным сто­

различные расстояния ме­

ком Q для глубокого (1,5

м) и мелкого (0,9

м)

жду

дренами

(рис.

17),

систематического

дренажа

при

£= 20

м.

а также для

 

систем,

име­

 

 

1958—1968 гг. Римейкас.

 

 

ющих различные глубины закладки дрен. В обоих случаях величина коэффициента корреля­ ции больше 0,9 (табл. 24).

Таблица 24

Характеристика корреляционной связи между средними суточными модулями дренажного стока при различных Е и t

Опытный

Е

м

участок

t

м

Ю

О'

Число член ряда

!

Средняя Уравнение регрессии ошибка уравне­

ний регрессии

Коэффициент

корреляции

Е'

по Е"

Е " по Е'

 

 

и t '

по /"

и t" по V

со

 

 

 

со

Кандава

Е '

=

14

м

80

0,908+0,014

ЧЕ' =

qE„ =

±0,46 ±0,08

 

t

=

1,2м

 

=

5, \2(fgr/ — =

0 ,16?£, +

 

 

£ " = 4 0

м

 

 

- 0 ,1 0 7

+ 0,051

 

 

t

=

1,2 м

 

 

 

 

 

 

 

 

 

Римейкас

Е =

0,9

м

85

0,931+0,010

Яг =

qt„=

±0,15 ±0,11

 

£

= 2 0

м

 

=

1 .1 8 ? ,,- =

0 ,73?,,+

 

 

Г =

1,5 м

 

 

- 0,071

+0,108

 

 

£

= 2 0

м

 

 

 

 

 

 

 

 

 

57

Чв'=и, л/Сс-га) Из приведенных данных следует, что показатели дренажного стока одной дре­ нажной системы могут быть использо­ ваны при прогнозировании стока в си­ стеме с другими Е и t.

2. Максимальный дренажный сток

Максимальный дренажный сток в зависимости от степени дренирова­ ния. Максимальный дренажный сток ха­ рактеризует способность дренажа отво­ дить избыточную воду в период наиболь­ шего притока. По величине максималь­ ного стока можно судить о водопроницае­ мости и водоотдаче почвы. Однако дан­ ных о максимальном стоке в литературе довольно мало.

На иловато-суглинистых почвах опыт­ ной станции Мидовланд (США) после интенсивного и продолжительного дождя (70 мм) И. Ниел [102] наблюдал макси­ мальный дренажный сток величиной 4,7 л/(с>га). По данным Флодквист, в Швеции наблюдался максимальный дре­ нажный сток величиной до 6,0 л/(с*га).

Особенно большие значения модулей дренажного стока, наблюдавшиеся в 1931 г. на опытном поле Фридрихштейна, приводит И. Роте [279]: 28/IV — 3,38 л/(с • га), 7/1V — 5,19 и 9,13 л/(с • га), 28/IV — 12,10 л/(с-га).

Рис. 17. Связь между средними суточными моду­ лями дренажного стока q для среднеглубокого дренажа (1,2 м), имеющего различные расстояния между дренами (£ '= 14 м и £"= 40 м). 1961-62 г.

Кандава.

По данным А. А. Ксензова

[80],

на тяжелых

почвах Карелии

срочные модули дренажного

стока

достигают

10—12 л /(с-га),

а максимальные средние суточные — 5,75 л/ (с • га). Исследователями отмечается, что большой дренажный сток яв­

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

58

Характерно, что величина

максимального дренажного

стока

<7макс меньше в тех областях,

где меньше значение дренажного

стока в водном балансе почв.

По данным А. Н. Леушева [84],

в Бе­

лорусской ССР <7макс наблюдался в 1962 г., когда его величина при

f= 1 0

м

была 1,566 л/(с-га), при £ = 1 5 м — 1,200 л/(с-га), а при

£ = 20

м — 1,008 л/(с • га).

Данные полевых наблюдений, полученные на легких суглинис­ тых и супесчаных почвах Латвийской ССР, показывают, что вели­ чина наибольших модулей дренажного стока в условиях атмосфер­ ного питания может превысить 5 л /(с • га), а в условиях смешанного питания даже 8 л/(с>га) (табл. 25; 26; 27). Отмечается, что наи­ больший дренажный сток, как правило, формируется весной во время снеготаяния, т. е. перед оттаиванием почвы.

 

 

 

 

 

 

 

 

Таблица 25

 

Характерные наибольшие модули дренажного стока на Кокнесском

 

 

 

 

 

опытном участке

 

 

 

 

Расстояние Гдубина

 

q л/(с*га)

 

Глубина

Глубина

 

 

 

1958

 

I960

промерзания

снегового

между

 

закладки,

 

почвы,

см

покрова,

см

дренами,

 

 

 

 

 

 

 

м

 

 

 

 

 

 

 

м

 

15/1V

16/IV

12/IV

1958

I960

1958

1960

 

 

Атмосфер ное питан ие

14

1,2

> 5 , 0

> 5 , 0

> 5 , 0

 

 

20

1,2

> 5 , 0

> 5 , 0

5 .2

 

 

20

0 ,9

> 5 , 0

> 5 , 0

1.2

5 - 1 0

5 0 - 6 0 3 6 - 4 0 3 0 - 3 2

 

 

 

Смешанное питание

 

 

14

0 ,9

> 1 0 , 0

> ю , о

2,4

 

 

30

0,9

3,5

4,3

1,°

 

 

25

1,05

4,1

9,5

 

 

 

За период наших исследований наиболее значительный весен­ ний дренажный сток с характерными резкими суточными его коле­ баниями, вызванными суточным колебанием температуры воздуха, наблюдался весной 1958 и 1960 гг. (рис. 18). При этом в 1958 г. сток проходил при незначительной, а в 1960 г. — при значительной глубине промерзания почвы. Так, на Кокнесском участке в первом случае глубина промерзания была 5—10 см, а во втором — 50— 60 см. Однако в обоих случаях в период наиболее интенсивного снеготаяния отводящие дрены железобетонных контрольных колод­ цев на опытных системах оказались не в состоянии отвести все притекающие дренажные воды (см. рис. 18), несмотря на то, что при проектировании этих систем был принят завышенный в не­ сколько раз расчетный дренажный сток. Надо подчеркнуть, что даже при глубине промерзания почвы более 0,5 м величина макси­ мального дренажного стока превысила 5,0 л/(с-га). Из этого сле­ дует, что мерзлый слой дренированных легких суглинистых и

59

Соседние файлы в папке книги из ГПНТБ