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

книги из ГПНТБ / Пак, В. В. Шахтные вентиляционные установки местного проветривания

.pdf
Скачиваний:
10
Добавлен:
21.10.2023
Размер:
10.03 Mб
Скачать

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

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

а — характеристики давления; б — схема расположения штуцеров

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

При работе без воздушного сепаратора (см. рис. 54) на периферии (в пределах лопаток) имеется явно выраженная область разрежения, распространяющаяся со стороны передних кромок лопаток на рас­ стояние около 35—40% величины их хорды. Эта область сохраняется в диапазоне от максимальной производительности до режима макси­ мального давления. По мере дальнейшего уменьшения производи­ тельности на аэродинамической характеристике вентилятора по­ является впадина, а в указанной области давление становится выше атмосферного. Такое повышение давления можно объяснить появле­ нием вращающегося срыва и выбрасыванием части возмущенного потока за пределы входных кромок рабочих лопаток.

При работе вентилятора с воздушным сепаратором (см. рис. 55) в точке над лопаткой, расположенной от входной кромки на рас­ стоянии около 40% длины хорды, разрежение при максимальном расходе существенно выше по сравнению с работой без сепаратора.

101

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

=±45°.

Проведенная серия испытаний показала, что для колес с мери­

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

Расход воздуха, проходящего через сепаратор, определяли с по­ мощью шести осредняющих расходомерных зондов, установленных равномерно по окружности в направляющем аппарате. По этим зон­ дам определяли расход воздуха через направляющий аппарат ()НА, а расход воздуха, подаваемого вентилятором в нагнетательный трубопровод, Qa замеряли по трубе Вентури. Таким образом, расход воздуха через канал сепаратора составлял

( ? с = <?в - < ? н а -

(137)

Осредняющие расходомерные зонды в направляющем аппарате предварительно были протарированы по трубе Вентури (при этом канал сепаратора закрыли). _

На рис. 56 показано изменение значения Qc в зависимости от режима работы вентилятора К-9Г- 2. Аналогичные измерения были проведены и для вентилятора К-10Г-16. Несмотря на значительное различие оптимальных параметров этих вентиляторов (см. табл. 20) качественно картина течения в сепараторе была подобной.

Испытания позволили проследить характерные зависимости тече­ ния в воздушном сепараторе:

расход воздуха через сепаратор практически мало изменяется от угла установки лопаток направляющего аппарата;

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

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

ствует Qu;

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

тора в количестве, примерно равном 0,6@н.

102

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

зонда

 

диаметром

10 мм

 

 

 

 

 

 

 

 

(рис. 57). Испытания про­

 

 

 

 

 

 

 

 

водили на модели К-9Г-12

 

 

 

 

 

 

 

 

при угле установки за­

 

 

 

 

 

 

 

 

крылков

Она = 0° на че­

 

 

 

 

 

 

 

 

тырех

режимах:

 

 

 

 

 

 

 

 

 

Qx= 0,314 —режим мак­

 

 

 

 

 

 

 

 

симальной производитель­

 

 

 

 

 

 

 

 

ности;

 

 

 

 

 

 

 

 

 

 

 

 

 

~QU= 0,236 — режим оп­

 

 

 

 

 

 

 

 

тимальной производитель­

 

 

 

 

 

 

 

 

ности;

 

 

 

 

 

 

 

 

 

 

 

 

 

Q1U=

0,11 — режим, со­

 

 

 

 

 

 

 

 

ответствующий

 

нижней

 

 

 

 

 

 

 

 

точке впадины на харак­

 

 

 

 

 

 

 

 

теристике давления;

 

 

 

 

 

 

 

 

Q1V ='0,055 —

режим,

 

 

 

 

 

 

 

 

близкий неполному дроссе­

 

 

 

 

 

 

 

 

лированию,

 

 

 

 

 

 

 

 

 

 

 

Поля скоростей и углов

 

 

 

 

 

 

 

 

потока

в

направляющем

 

 

 

 

 

 

 

 

аппарате и воздушном се­

 

 

 

 

 

 

 

 

параторе

позволили уста­

 

 

 

 

 

 

 

 

новить,

что:

 

 

 

 

 

 

 

 

 

 

 

значения

с 1а

осевой

0 ,0 5

 

 

 

 

 

+J2

Z .

скорости в межлопаточном

 

 

 

 

 

 

 

 

 

 

 

 

 

 

канале

 

направляющего

0,0

 

 

 

 

 

 

°

аппарата

на соответствен­

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

ных режимах

 

мало изме- . qq$

 

 

/

/

 

 

 

няются вдоль радиуса;

 

>

 

 

 

 

 

 

/ /

 

 

 

 

 

при

работе

на режи­

 

у

/ /

 

 

 

 

 

мах

Q1U

и

Qiy

выброса

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

потока со стороны рабо­

 

 

 

 

 

 

-

 

 

 

 

 

 

 

 

 

чего колеса в межлопаточ- ~ 0 ,1 5

0 ,1

0 ,1 5

0 ,2 0

0 ,2 5

0 ,3

 

ном канале не обнару­

0 ,0 5

 

 

 

 

 

 

 

 

 

жено;

канале

сепаратора

Рис. 56.

Расход воздуха через канал сепа­

в

 

ратора

при 0НА:

 

 

между обечайкой и корпу­

 

------------ 0е ;

---------------- (-40°

 

сом на

режиме

QT имеет

 

 

 

 

 

 

 

 

 

 

место

подсос

потока со скоростью

с 1а с,

примерно

в

3—4 раза

меньше скорости с1а на этом режиме;

на оптимальном режиме @п подсос потока через канал сепаратора не наблюдается;

103

на режимах Qlu и Qlv в канале сепаратора имеет место выброс возмущенного потока из рабочего колеса в направлении, противо­ положном движению основного потока в межлопаточном канале направляющего аппарата; скорость обратного возмущенного по­

- 4 0

Канал Воздушного сепаратора,

Оfa S Сtа с

_^ - - — — _

тока с 1ас на режиме Qlv примерно в два-три раза

больше скорости с'la на этом режиме;

выбрасываемый в канал сепаратора поток на ре­

жимах Qnl

и Q1V сильно

закручен в

направлении

вращения

рабочего ко­

леса; степень закручива­ ния потока на этих режи­ мах примерно одинакова и характеризуется по от­ ношению к оси вентиля­ тора углом потока ас = = 130°.

ОЛ

 

 

 

 

N .

 

Измеренные двумя опи­

 

 

 

 

 

 

 

санными способами

значе­

0,2

 

 

 

 

 

 

ния

расхода

воздуха,

 

 

 

 

 

 

^ -----V,

проходящего

 

через

ре­

0

 

 

 

 

 

 

шетку направляющего ап­

 

 

 

 

 

 

парата и канал сепара­

 

 

 

 

 

 

 

тора,

удовлетворительно

- 0,2

 

 

 

 

 

 

согласуются между собой,

 

 

 

 

 

 

ч

что свидетельствует о до­

- 0 , 4

 

 

 

 

 

стоверности

полученных

 

 

 

 

 

 

 

результатов. Характер те­

0 ,5

0 ,6

0 ,7

0 ,8

0 ,9

1 ,0

1,1 F

чения в воздушном сепара­

торе вентилятора К-9Г-12

 

 

 

 

 

 

 

Рис. 57. Поля скоростей чо

и с1 ао и углов

хорошо согласуется

с ре­

зультатами

исследований

выхода потока ас

в

воздушном

сепараторе

вентилятора К-9Г-12 при

Q:

 

лопаточного

 

сепаратора

---------- 0,3 1 4 ;------------- с 0,236;

 

 

0,055

осевого вентилятора К-06

 

 

 

 

 

 

 

[55].

Это дает

основание

ная скорость

потока

с 1а,

 

 

считать, что максималь-

выбрасываемого

при

срыве

через

канал сепаратора на режимах, близких к полному дросселированию, зависит от окружной скорости колеса на периферии uR и мало зависит от угла установки лопаток направляющего аппарата (или рабочего' колеса), а также от типа осевого вентилятора. Так, например, для вентилятора К-06 указанная скорость составляла с 1ас — 0,8uR, а для вентилятора К-9Г-12 — составила с1ас = 0,75uR.

Используя полученные количественные зависимости, можно опре-

104

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

= Д,

(138)

где # — наружный диаметр рабочего колеса.

Пренебрегая толщиной обечайки, которая обычно не более 0,005#,

определим площадь кольцевого

канала сепаратора

/ — 0,6(?н че­

 

 

 

 

 

 

0,75и,

рез внутренний диаметр

корпуса

направляющего

аппарата Дик

и диаметр#. Эта площадь выражается как

U,/OUr

=

Ц- (Дик — # а),

 

 

 

 

 

откуда приближенно # На =

[ / —

+ # .

 

 

 

 

 

Г Up

 

 

 

 

Выразим далее QHчерез

безразмерный коэффициент производи­

тельности, после чего будем иметь

 

 

 

 

# н а

= # 1 / 1 + 0,785<?я.

 

(139)

Выражение (139) получено по результатам экспериментального исследования потоков в канале сепаратора и является приближен­ ным. Однако, как показал опыт, оно вполне приемлемо для инженер­ ных расчетов при определении конструктивных размеров противосрывных устройств вентиляторов с меридиональным ускорением.

Дальнейшие экспериментальные исследования воздушного сепа­ ратора позволили определить влияние размеров его отдельных элементов (см. рис. 53) на аэродинамические характеристики венти­ лятора с меридиональным ускорением потока.

Выступание лопаток

в канал сепаратора исследовано в диапа­

зоне =

(0,01 -f- 0,1) # .

Оптимальным оказался размер 1К —

= (0,035

0,045) Д , что в долях хорды лопатки на периферии колеса

составляет (0,3 -f- 0,35) Ъ.

Расстояние от обечайки сепаратора до лопаток колеса в диапа­ зоне б = (0,0065 0,045) # мало влияет на характеристики, макси­ мальный к. п. д. вентилятора при этом изменяется на 1,0—1,5%. Рекомендуемые значения этой величины составляют б = (0,015 —

-f- 0,02) # .

Угол перехода от корпуса направляющего аппарата к корпусу вентилятора исследован в диапазоне ак = 20 -f- 90°. Изменение угла перехода в этом диапазоне слабо влияет на характеристики, раци­ ональными можно считать значения ак = 60 4- 90°.

Длина обечайки канала сепаратора 10 влияет на эффективность раскручивания потока в нем. Оптимальная длина обечайки при

безлопаточном канале

сепаратора составляет 10 = (0,35

0,45) # .

В тех случаях, где

по конструктивным соображениям

длина

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

105

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

спрямляющая решетка из 12—14 лопаток при густоте тс =

1,2 1,3,

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

около 45°

и осевой выход из его канала при срывных режимах.

 

Ширина канала на выходе при работе с входным коллектором (как, например, у вентиляторов ВМ-8м и ВМ-12м) должна составлять

к =

(0,06 -f- 0,085) Д ,

причем меньшее значение соответствует Q =

Чу

 

 

 

 

=

0,2 -f- 0,25,

а

большее

 

 

 

Схема /

О

 

0,3 0,35.

 

 

0,8

 

 

 

 

Q =

 

 

 

 

 

 

 

При работе без кол­

0,6

 

/

\

 

лектора

 

(вентиляторы

 

 

//

 

ВМ-5м и ВМ-бм) расстоя­

 

 

<

*\

 

ние

от

фланца

входного

ИУ

 

 

 

 

патрубка до обечайки се­

 

 

 

Схема 2

паратора

увеличивается

ОД

 

 

 

 

до

 

к — (0,15 -f- 0,18) Д.

 

~~---

 

 

 

В

последнем

случае на

 

 

 

 

входе в вентилятор

обра­

0,3

 

/

\

 

зуется патрубок с кольце­

\

 

вым уступом, который,

 

Л _

/

\

Схема J

согласно

исследованиям

 

\

/

\

 

[59], обеспечивает сравни­

од

 

•j

\

 

тельно

малые потери при

 

 

входе в вентилятор. Коэф­

 

 

 

ГЛ\

 

фициент

 

потерь

в

этом

 

 

 

т

 

случае

составляет

около

 

 

 

 

 

£с =

0,23,

а

при

уста­

 

0,1 0,2

0,3 ОД О

новке на входе коллектора

 

с

радиусом закругления

 

 

 

 

 

Рис. 58. Характеристики вентилятора К-9Г-12

г =

 

0,1Д

соответственно

при работе с сепаратором и без него по схеме:

£к =

0,11

[60].

Если вход

 

 

г ; -------

вентилятора

выполнен

 

 

 

 

 

просто,

в

виде

 

патрубка

с фланцем, как, например, у серийных вентиляторов местного проветривания типа СВМ, коэффициент потерь при входе соста­ вляет около £п = 0,57.

На рис. 58 показаны характеристики вентилятора К-9Г-12 (0на=

=0°) при работе с входным коллектором без сепаратора (схема 7),

свходным коллектором и воздушным сепаратором (схема 2), с коль­ цевым уступом и воздушным сепаратором (схема 3).

_ Сравнение характеристик показывает, что в диапазоне режимов

Q = 0,21 0,36 применение сепаратора с входным коллектором приводит к снижению давления и к. п. д. примерно на 1% (по срав­

нению со схемой 1). При Q < 0 ,2 сепаратор обеспечивает устранение впадины и значительное повышение к. п. д. Напорные характери­ стики при работе вентилятора по схемам 3 ж 2 практически совпа­

106

дают, а к. п. д. в области оптимума снижается на 1—2% . Регулиро­ вочные характеристики вентилятора К-9Г-12 с сепаратором по схеме 2 (при 0на = ±50°) показаны на рис. 45.

§ 4. Пример аэродинамического расчета рабочего колеса с меридиональным ускорением потока

Последовательность аэродинамического расчета проследим на примере вентилятора со встроенным электродвигателем по схеме НА + К + СА, который предназначен для местного проветривания и должен обеспечивать на оптималь­ ном режиме следующие параметры установки:

 

полное давление Ну,

 

к гс/м 2 ....................................................

 

185

 

производительность Q, м * /се к ................................................

 

3,1

 

полный к. п. д. т)у.....................................................................

 

 

 

 

0,75

 

скорость вращения колеса п, об/мин......................................

 

2950

 

диаметр присоединительных патрубков dT, м ........................

 

0,5

 

Определим дополнительные расчетные параметры вентилятора. Диаметр

рабочего колеса принимаем Д =

0,495 м. Исходя из габаритов встраиваемого

.электродвигателя относительный диаметр втулки на выходе d2 =

0,7 и d2 —

= 0,345

м. Выбираем степень ме-

1

 

ридпонального

ускорения

(см.

____------ -

 

рис.

44)

щ = 1,35, тогда относи­

 

 

 

тельный диаметр втулки на входе

 

 

 

Jx =

V 1 — Tit (1 — dt) =

0,555;

 

 

 

dx =

0,275.

 

 

 

при­

 

 

 

 

Число рабочих лопаток

 

 

 

нимаем

гл =

7.

Принимаем

 

 

 

пг =

0

в связи с

малым

значе­

 

 

 

нием расчетного давления {Ну<=&

0,3). Угловая скорость ш =

ТТп

=3Q = 308,8 1/сек.

Осевая скорость потока

на

 

 

входе в колесо и на выходе из

 

 

него составляет:

 

 

 

 

4Q

=23,2 м/сек;

 

 

я(Д»-«Ч)

 

 

 

 

 

 

4Q

= 31,2 м/сек.

 

 

я(Д*-«г*)

Рис. 59. Построение профиля сечения

Плотность

перемещаемого

 

лопатки рабочего колеса

воздуха р = 0,122

кгс сек2/м4.

сечений

z = 5, значения гидравлических

Принимаем число расчетных

радиусов на входе и выходе пз сечений определяем по формуле

 

= |/~Д2 — (2ч0 — 1)

,

где R — наружный радиус колеса; п0 — порядковый номер расчетного се­

чения, начиная с периферийного;

г — радиус втулки рабочего колеса соответ­

ственно на входе и выходе;

z — число расчетных сечений.

По заданным и дополнительным расчетным параметрам определяем коорди­ наты расчетных сечений, углы треугольников скоростей на входе и выходе, яагруженность тсрСж . Результаты расчетов приведены в табл. 22 (схему тре­

угольников скоростей и расположения сечений см. на рис. 40).

107

о

00

Расчетная формула

Г1 '— j / ”R2—(2п0 — 1)

га =

"|/ № — (2л0 — 1)

.

 

1

 

rcp=-2~ (ri+ r2)

 

 

2ягСр'

 

KX= 0 ) r i .................................

.....................

и2 = <ог2 .......................................................

 

с

 

 

2“

P“ 2l]y ......................................................

 

■di= arctg«i/cl a ........................................

 

0 2 = arctg — — .............................

с2а

у2 = arctg с2ц/с2 а ........................................

Ц>1=УГс?д+“! • • • ..............................

»2 = УГcja+ (lt2—С2ц)2 ..........................

^ ^ a r c t g - ^ T - d - ^ i t g ^

Лр

А____c2a tg ^гидНCjatg 'в‘1

0mM~ arctg---------------------------- • ■

- п _4 (cXa tg ^ 1 —c2a tg ©2ИД) (cos ид

•'Ср'-'Ж—

,--------------------- -

«2в+ с1в

 

 

 

 

Т а б л и ц а

22

Размерность

Значения параметров в расчетных сечениях

 

1 -1

II—II

111-Ш

IV -IV

V—V

 

 

 

мм

339

220

200

178

152

 

мм

241

228

214

198

182

 

мм

240

224

207

188

167

 

мм

215

201

186

169

150

 

м/сек

73.7

68,0

г61,8

54,9

47.0

 

м/сек

74,4

70.3

66,6

61,1

56.0

 

м/сек

26.8

28.4

30.0

32.6

35,6

 

градус—мин

72—30

71-10

69-25

67-5

63-45

градус—мин 56-45

53-20

49-35

42—25

33-10

градус—мин

40-40

42-15 43-50

46-15

48—

45

м/сек

80,0

71,9

66.0

59.6

52,4

 

м/сек

52,0

52,3

48,1

41,1

37,3

 

градус—мин

54-5

50—20

46—15

38-15

28—10

градус—мин

65

62-45

60

55-35

4 9 -

30

 

0,933

0,945

0,967

1,08

1,18

 

I

 

Расчетная формула

тср =

"ТсрСж

~~п ...............................................................

 

Ь Ж

до_

(^i— tg (23+ 9тСр)+ 'в'гид

Рр

1,365+0,28тср+0,05с ...................

j v

( ^ - к?) + tg [44,6+ (0,3+ 0,08с) тср] - Раид

°3,17+0,45тср+ 0,016с

е°=90°— р ? ........................................................

6 = Тернер ...........................................

/ = Ь

1 %

 

 

10 0

 

 

100

 

 

Ъ*

Ь0,5/

Re ~ "g/

<р—arcsin

Ь

 

 

2Rc

 

 

 

 

Т а б л и ц а 23

Размерность

Значения параметров в расчетных сечениях

1-1

1 1 -и

1 II-IH

IV—IV

V—V

к

0,883

0,898

0,915

0,947

1,00

градус—мин

59-35

56-35

53-20

47-25

40

%

5,2

5,9

6,35

8,07

9.83

градус—мин

30-25

33-25

36-40

42—35

50

М М

. 190

180

170

160

150

мм

9,9

10,7

11,0

13,0

14,7

мм

9,5

11,0

12,4

13,7

15,0

мм

462

390

345

257

198

градус—мин

11—55

13-25

14-20

18-10

22—20

мм

193

183

173

163

154

Принимая на периферии в сечении I —I коэффициент силы Жуковского Сж = 1,0 (см. табл. 19), а у втулки (учитывая сравнительно малые аэродина­

мические нагрузки) в сечении V—V величину Сж = 1,2, определим густоту тср и длину хорды профиля на этих сечениях. Зададимся линейным изменением хорды и максимальной толщины профиля с по высоте лопатой, принимая в се­

чениях I —I и V—V значение с соответственно равным 5 и 10% . Входной угол атаки выбираем постоянным для всех сечений, т. е. а х = 0°. Далее вычислим параметры, характеризующие геометрическпе размеры профилен по сечениям, результаты которых приведены в табл. 23.

На полученные для каждого сечения средние линии накладываем симметрич­ ный профиль С-4, координаты которого приведены в табл. 24.

 

 

 

 

 

Т а б л и ц а 24

x — x j l ......................................

0

0,05

0,10

0,20

0,30

0,40

0,50

У = У/1......................................

0

0,308

0,402

0,483

0,500

0,489

0,457

х х/1......................................

0,60

0,70

0,80

0,90

1,0

гilc

гг/с

v = y / i......................................

0,405

0,337

0,254

0,160

0

0,12

0,06

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

откладывают на соответствующем от входной точки расстоянии х — xl значение

у — ус, согласно табл. 24. Полученные точки контура профиля соединяют плавной линией с соответствующими радиусами на входе и выходе (рис. 59). Углы 0Копределяют графически по вычерченным сечениям. Один из углов 0К, например на периферии в сечении I —I (см. рис. 40), принимают в качестве угла установки лопатки колеса. Полученные таким образом профили изгибают на конических поверхностях соответствующих сечений, устанавливая их под углами 0С. Затем определяют положение центра тяжести профиля известным способом [44]. Совместив центры тяжести на одном радиусе, получают основные контуры лопатки.

Г л а в а V

РАЗРАБОТКА РЯДА ОСЕВЫХ ЭЛЕКТРИЧЕСКИХ ВЕНТИЛЯТОРОВ

Широкий диапазон вентиляционных режимов и условий про­ ветривания тупиковых выработок вызывает необходимость разра­ ботки типоразмерного ряда осевых вентиляторов местного проветри­ вания. Правилами безопасности предусматриваются вентиляторы местного проветривания с электрическим и пневматическим при­ водом. Действующий в настоящее время отечественный стандарт на осевые ВМП (ГОСТ 6625—64 *) предусматривает пять типоразмеров электрических и три типоразмера пневматических вентиляторов. В последние годы в дополнение к стандарту были разработаны мощ­ ные пневматические вентиляторы ВМП-бм диаметром рабочего колеса 600 мм и вентиляторы ВМ-12м диаметром рабочего колеса 1200 мм для горнорудной промышленности.

110

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