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

книги из ГПНТБ / Местер, И. М. Автоматизация контроля и регулирования рудничного проветривания

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

воздействия и соответствует сложному дифференциальному звену. Наличие в числителе оператора р свидетельствует о том, что со­ ставляющая переходного процесса, обусловленная этим членом, имеет место только при скорости изменения во времени входного воздействия Q (t), отличающейся от нуля. Знак коэффициента Къ

положительный, следовательно, приращению Q (t) соответствует

 

 

 

 

приращение С (t) того же

 

 

 

 

знака.

 

 

выражения

 

 

 

 

 

 

Второй член

 

 

 

 

(1.10)

описывает приближен­

 

 

 

 

но объект без учета эффек­

 

 

 

 

та

выработанного

простран­

 

 

 

 

ства и соответствует зави­

 

 

 

 

симости

(1.3).

отрицательного

 

 

 

 

 

 

Наличие

 

 

 

 

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

 

 

 

 

ной функции (1.10) свиде­

 

 

 

 

тельствует о том, что, по

 

 

 

 

крайней мере, один ее нуль

 

 

 

 

может оказаться

и

в

правой

 

 

 

 

полуплоскости

поэтому

 

 

 

 

объект можно считать неми­

 

 

 

 

нимально-фазовым.

 

 

 

 

 

 

 

 

Реакция объекта, описы­

 

 

 

 

ваемого

выражением

(1.10),

 

 

 

 

на

 

скачкообразное

возмуще­

 

 

 

 

ние

по

управляющему воз­

 

 

 

 

действию ДQ (t) при AG (t)п =

 

 

 

 

-

0

 

 

 

 

 

 

 

 

 

 

 

AC(t) = ACE(t) + ACx (t),

 

Т Ц

Р ЛЛГГО TIТІ ПТАТГПІП ПЛППѴ ГіТІТГГіТ» гК ч ттп л

 

 

 

 

, ч

 

 

 

( 1- 11)

ТІЧЛ

 

 

Л

 

 

 

 

 

jTiiL/.

і . и .

и с и е л и д д и и

ш ѵ и л -

^де

 

составляющая,

 

 

ции объекта

 

 

 

\і)

 

 

 

 

обусловленная

выработан­

ным пространством; АС% (t) — составляющая, обусловленная лавой. Переходя от изображений к оригиналам и полагая без нарушения

общности,

что Г 2 =£ Т х, получаем

 

 

= > 2- ВУ1 ( е -^ ° - е - ;/г .);

(1.12)

 

AC(th = Kw ( i - e t/Ti).

(1.13)

На рис.

1.5 показан примерный вид переходных функций АС (t);

AC, (t) и АС (t).

 

В некоторых работах [1, 15, 16] отмечается, что скачкообразное возмущение по Q (t) не всегда вызывает характерный всплеск С (t). При принятой нами математической модели объекта всплеск дей-

20

ствителыю не наблюдается, но только при определенных соотноше­ ниях параметров объекта 7\, Т 2, Тх, Къ и KXq.

Приравнивая нулю при t = 0 импульсную переходную функцию объекта, находим искомые соотношения

Кв

K XQ

(1.14)

T i T o

т х

 

При плавном нарастании

управляющего воздействия АQ (t)

и начальных нулевых условиях получаем следующие выражения для переходных процессов:

для

составляющей АС (t)B

 

 

АC(t), = K,Q ( l + Ч [ ю Н - < / Ч ) ] - Ч І ° » Р ( - < / Ч )1 у'

(І16)

для

составляющей АС (t)x

 

 

AC(t)x = K XQQ ( t -

(1.16)

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

строго с учетом реальных

физических процессов, протекающих

в объекте и, в частности, в

выработанном пространстве [14].

По данным исследований,

длительность переходных газодинами­

ческих процессов, обусловленных значительным ступенчатым воз­ мущением по воздуху, колеблется в широких пределах (от десятков минут до нескольких часов и даже суток), что полностью подтвержда­ ет высказанное в работе [4] предположение о существенной разноинерционности аэродинамических и газодинамических процессов на добычных участках газовых шахт. Установлено также, что параметры, определяющие динамические свойства объекта [в случае переда­ точной функции (1.10) — это значения Тх, Тг, Г2, Къ и K Xq ], могут изменяться в широких пределах от объекта к объекту и, кроме того, изменяться с течением времени у одного и того же объекта, что значительно усложняет задачу построения системы автоматической стабилизации концентрации на исходящей добычных участков.

Выполненные нами исследования в условиях ряда шахт Кара­ гандинского угольного бассейна, а также аналогичные работы для условий шахт Донбасса [17] показали, что динамические свойства выработанного пространства не проявляются при идентификации объекта методами статистической динамики по данным пассивного эксперимента (рис. 1.6). Это объясняется, по-видимому, тем, что в процессе нормальной эксплуатации резкие возмущения по воздуху возникают очень редко, как правило, лишь в аварийных случаях. Как видно из рисунка, максимумы оценок нормированных взаимокорреляционных функций р (t)cq и р (т)cg процессов С (t) Q (t) и С (t)

21

ß t * h

G (t) почти совпадают с осью ординат (постоянная времени 7\. значительно меньше интервала корреляции, равного 4 ч для кривой а и 10 ч для кривой б, при этом р (0 )Gq > 0 , а р (0 )со < 0 , что соот­ ветствует основным статическим свойствам объекта.

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

стотными

 

методами.

На

а

АС

 

 

 

 

 

 

 

рис. 1.7 приведены обоб­

 

Айк*,а

 

 

 

 

 

 

 

щенные переходные и нор­

 

 

 

 

 

 

 

 

 

 

мированные

амплитудно-

 

 

 

 

 

 

 

 

 

 

частотные характеристики

 

 

 

 

 

 

 

 

 

 

объекта,

полученные

по

 

 

 

 

 

 

 

 

 

 

данным

активных

экспе­

 

 

 

 

 

 

 

 

 

 

риментов

 

[1,

4].

 

 

 

 

 

 

 

 

 

 

 

 

 

Кривые

1 характерны

 

 

 

 

 

 

 

 

 

 

для

добычных

участков

 

 

 

 

 

 

 

 

 

 

шахт Донецкого бассейна,

 

 

 

 

 

 

 

 

 

 

кривые 2 — для

наиболее

 

 

 

 

 

 

 

 

 

 

тяжелых

(с точки

зрения

 

 

 

 

 

 

 

 

 

 

динамики переходных про­

 

 

 

 

 

 

 

 

 

 

цессов) условий некоторых

 

 

 

 

 

 

 

 

 

 

участков

шахт

Караган­

 

 

 

 

 

 

 

 

 

 

динского

бассейна.

 

что

 

 

 

 

 

 

 

 

 

 

Из рисунка видно,

 

 

 

 

 

 

 

 

 

 

резонансный пик АЧХ на­

 

 

 

 

 

 

 

 

 

 

ступает

в

первом случае

 

 

 

 

 

 

 

 

 

 

при

(ор =

0,01 рад/мин,

 

 

 

 

 

 

 

 

 

 

во

втором — при

сор

=

 

 

 

 

 

 

 

 

 

 

= 0,1 рад/мин.

 

считать

 

 

 

 

 

 

 

 

 

 

Объект

можно

 

 

 

 

 

 

 

 

 

 

практически безынерцион­

 

 

 

 

 

 

 

 

 

 

ным

(I т] (со) I

— I г)

(0) |) С

 

 

О

0,2

0,4- 0,6 0,8

1,0 1,2 ш, рад/мин

</ 0,05 при

полосе про­

Рис.

1.7. Характеристики реальных

газоди­

пускания: в первом слу­

 

 

 

намических объектов:

 

 

 

чае

До =

0,017 рад/мин,

а — переходные характеристики: 1

=

10

мин;.

во втором

случае

Аса =

Т,

=

5

мин;

Г2 = 15 мин;

Кв =

35

мин-1 ;

2

= 0,0025 рад/мин.

 

 

 

=

40

мин;

Г, = 50 мин;

Т2 =

130

мин; Кв =

 

 

=

350 мин- ‘; б — нормированные АЧХ (кривые 1

и 2),

Таким образом, рас­

соответствующие переходным

характеристикам 1 и 2,

сматриваемый объект даже

 

 

и разностная АЧХ (кривая з

— I

Др (ю) |)

 

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

Если допустить, что в процессе нормальной эксплуатации газовыделение может изменяться в 2,5—3 раза, то в не меньшей степени должно изменяться и управляющее воздействие Q 0, а следовательно,

23-

в таком же диапазоне будут изменяться постоянная времени объекта (1.4) и его коэффициент усиления по обоим каналам [(1.5) и (1.6)]. В результате объект даже при значительной его идеализации описы­ вается дифференциальным уравнением с переменными, случайно изменяющимися во времени коэффициентами.

Рассмотрим подготовительные участки угольных шахт. Подго­ товительным участкам как объектам контроля и регулирования по фактору проветривания уделяется все большее внимание, особенно в условиях метанообилышх шахт восточных районов страны. Доля общешахтного расхода воздуха, необходимого для проветривания подготовительных участков, непрерывно возрастает и доходит, например, по шахтам Карагандинского угольного бассейна до

40-50% .

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

1. Основными причинами нарушений проветривания подгото­ вительных забоев являются отключения электроэнергии (48,7%), плановые остановки подготовительного забоя (25,3%) и вентилятора главного проветривания (21,4%). Среднее время нормального функ­ ционирования подготовительного забоя (между нарушениями про­ ветривания) 125 ч. Среднее время нарушения проветривания под­ готовительного забоя 3,76 ч. Среднее количество нарушений про­ ветривания подготовительного забоя составляет 68 раз в год, при

•этом суммарное время нарушения проветривания 255 ч.

2.Продолжительность остановки ВМП в 85% случаев превышает 30 мин, что в большинстве случаев достаточно для загазирования подготовительного забоя.

3.В 70% общего количества простоев ВМП их остановка вызвана причинами, не зависящими от условий работы подготовительных забоев, при этом средняя длительность простоя ВМП составляет 2,5 ч. Таким образом, резервирование электроснабжения ВМП имеет ограниченную эффективность.

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

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

Для автоматизации этого опасного и сложного процесса необхо­ димо прежде всего выявить специфические особенности объекта управления.

24

Обобщенную модель процесса (рис. 1.8) разгазирования неза­ висимо от метода можно представить как взаимодействие выталки­ ваемого из тупика длиной L потока метана GT с потоком чистого воздуха Q, подаваемого на разбавление в условную зону смешения длиной X, через которую протекает поток смешения QCM.В зависимости от конкретных условий углы а и ß могут изменяться в широких пределах. Измеряемая на расстоянии х от начала условной зоны смешения концентрация метана является регулируемым параметром. Часть Qn общего потока Q, подаваемая ВМП для выталкивания метанового облака из тупика, представляет собой регулирующее

воздействие.

Gr

и

Q

следует

 

 

 

Параметры

 

 

 

рассматривать как внешние слу­

 

 

 

чайно изменяющиеся

в

процессе

 

 

 

разгазирования возмущения.

 

 

 

При

неизменном

значении Q

 

 

 

параметр

GT является функцией

 

 

 

общего количества метана G0, вы­

 

 

 

делившегося за время Ѳо наруше­

 

 

 

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

 

 

 

ния последнего по длине хтв тупи­

 

 

 

ковой выработке

 

 

 

 

 

 

 

 

ѳ L

 

 

 

 

Рис.

1.8. Обобщенная модель процес­

G0(Ѳ,

L) — j" j

q (Ѳ,

xT)dxrdQ,(1.17)

са

смешения при

разгазировании

 

 

о о

 

 

 

 

 

тупиковой выработки

где

q (Ѳ, жт) — погонное

абсолютное

газовыделение

как функ­

ция

длины

тупиковой

выработки

и

длительности

загазирова-

ния,

м3/м-сек.

 

 

 

 

 

 

Внутренний интеграл в выражении (1.17) представляет собой общее количество газа, выделяющегося в единицу времени, как функция Ѳ (м3/сек)

Т\

 

/ (Ѳ) = J q (Ѳ, X T) dxT.

(1.18)

о

 

В большинстве случаев можно считать / (Ѳ) невозрастающей функцией времени.

Общее количество метана, выделившегося за время Ѳ0, с учетом

(1.17) и (1.18)

е„

 

G(0) = j f(Q)dQ.

(1.19)

о

 

Если с момента Ѳ = Ѳ0 начать разгазирование и предположить, что концентрация метана С будет при этом поддерживаться по­ стоянной и равной допустимому значению Сдоп (относительные

25

единицы), то длительность разгазирования тр может быть определена

из уравнения [26]

Ѳ„+Тр

Сдоп<?тр-

J /(Ѳ)<Ю =

0.

(1.20)

 

о

 

 

Если предположить, что за

время разгазирования / (Ѳ) = /0 =

= const, то выражение (1.20)

упростится

и при

этом

_

/обо

 

( 1. 21)

Р

Сд0ПѲ —/о

 

 

 

Из выражений (1.20) и (1.21) следует, что при Q = const расчет­ ное значение тр представляет собой оптимальную (минимально возможную) длительность процесса разгазирования.

Горизонтальная выраіат ка

БремсБерг

Уклон

Рис. 1.9. Распределение концентрации метана по длине в условиях различных загазированных подготовительных выработок

Если при прочих равных условиях осуществить разгазирование за время, меньшее, чем тр, то концентрация на исходящей струе будет превышать значение 6'доп, что недопустимо. Поддержание же на исходящей струе более низкой концентрации, чем Сдоп, приводит к резкому возрастанию длительности разгазирования, что также нежелательно.

Из приведенных зависимостей следует, что количество GT метана, подаваемого в зону смешения, зависит от трех переменных пара­ метров: Ѳ, хТ и Qn.

26

При неизменном значении Qn Qu 0 зависимость GT (т), отсчи­ тываемая с момента начала разгазирования, воспроизводит в опре­ деленном масштабе, зависящем от значения Qn 0, функцию

Я (^т)ѳ=Ѳо*

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

На рис. 1.9 показаны экспериментальные кривые Ст— / (хт), полученные в результате обработки данных [191 и показывающие распределение концентрации метана по длине тупиковой выработки при различной длительности Ѳ загазироваиия.

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

в

работах [19, 201. Наиболее высокий градиент изменения концен­

трации имеет

место в

восстающих выработках, наименьший —

в

наклонных.

Отдельные

всплески Ст свидетельствуют о наличии

местных усиленных выделений метана в районе того сечения выра­

ботки, которому соответствует максимум всплеска

кривой Ст —

=-■ / О т ­

что параметр

данные работы [20] свидетельствуют также о том,

GT является переменным во времени, а максимальная скорость его изменения во времени при прочих равных условиях тем больше, чем больше длительность загазироваиия Ѳ0. В связи с тем что концентра­ ция нарастает от устья к забою тупиковой выработки, при разгазировании с постоянной скоростью выталкивания (постоянная произ­ водительность ВМП) параметр Ст сначала возрастает, достигая своего максимального значения, соответствующего однократному обмену воздуха в объеме тупиковой выработки, затем спадает из-за разбавления заднего фронта метанового облака выталкивающим потоком воздуха до установившегося значения, соответствующего режиму нормального проветривания.

При автоматической стабилизации концентрации С на исходящей струе величина (?п должна изменяться во времени по закону, обрат­ ному зависимости Ст= / (жт), что необходимо для обеспечения постоянства величины Gr в течение процесса разгазирования и, сле­ довательно, минимальной его длительности (рассматривается наи­ более тяжелый случай, когда Ст> С*,*» — Сдоп).

Ввиду невозможности достаточно точного определения функции Ст= / (жт), а также из-за случайного характера изменения во времени количества воздуха С стабилизацию величины целесообразно осуществлять с помощью замкнутой САР, при этом объект управле­ ния представляет собой два последовательно включенных звена: тупиковую выработку (ТВ) и зону смешения (ЗС).

Входной величиной звена ТВ является количество подаваемого на выталкивание воздуха Qn, выходной — количество выталки­ ваемого газа GT. Расчеты, выполненные на основании данных,

27

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

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

Для математического описания 3G могут быть использованы различные по сложности модели: идеального выталкивания, идеаль­ ного смешения, одно- и двухпараметрической диффузии [13].

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

- + АС -

К $ (Т) АQn+

к № (Т) АСт,

(1.22)

где

V

 

 

Т (х) ~

 

(1.23)

<?(0) (т) + <2<°> (т)

 

iHQ) /ѵ\

Ст0) (т) —Сдоп

(1.24)

Лоб 1 4 -

С;(0)(т) + (?(0){х) »

iAG) / ч

<?п0) (т)

 

(1.25)

обВ 4

()'«>(г) + (?Тй) '

 

Так как в процессе разгазирования параметры ЗД°) (т),

(Е°) (т)

и С(.°) (т) могут изменяться (с учетом аварийных режимов) в широких пределах (от нуля до максимального конечного значения), то как постоянная времени объекта Т (т), так и его коэффициенты усиления по каналам С Qn и С Ст могут при этом теоретически изме­ няться от минимального, конечного значения до бесконечности. Практически же диапазон изменения этих параметров, хотя и ко­ нечен, весьма широк, что в первую очередь следует учитывать при разработке системы автоматизированного управления процессом разгазирования опасных по газу тупиковых выработок.

Так, например, по данным работы [20] отношение максимального газовыделения к минимальному из тупиковой выработки составляет около 8. Если при этом допустить, что во время разгазирования по непредвиденным причинам в два раза уменьшилась величина @(0) (т), то коэффициент усиления объекта К^> (т) может измениться в про­ цессе разгазирования в 16 раз.

Возможность отнесения рассматриваемого объекта к классу систем с медленно изменяющимися параметрами зависит от характера изменения во времени усредненного значения концентрации Ст

28

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

При принятых допущениях

Ст (т) = /(;гт)-

(1-2(5)

Если пренебречь диффузионными явлениями в процессе разга­

зирования, то можно считать, что

 

=

j <?п('с)йт + 4 0)-

(1-27)

где S r — сечение тупиковой

выработки; хт10) — путь,

пройденный

газовым облаком с начала разгазирования.

Если предположить, что САР процесса разгазирования функ­ ционирует нормально и обеспечивает качественную стабилизацию концентрации на уровне Сдоп за счет поддержания соответствующего значения (7£0), то при СТ % С в любой момент времени выполняется условие

Qn (т) С'т (т) = = const. (1.28)

Учитывая (1.27) и (1.28), получим

Ч-^)г7Ы+*“’}-с-м-0- (Ш>

Решение уравнения (Е29) с целью получения в явном виде функции Ст(т) весьма затруднительно в общем виде, однако в частном случае, когда /Ст является линейной функцией хт, т. е.

/ст (Хт) = а о) + ахт,

(1.30)

оно имеет следующий вид:

Ст(т) = ] / - ^ - т + (С П 2.

(1-31)

где а — градиент изменения Ст(хт) вдоль длины тупиковой выра­ ботки.

Можно показать, что скорость изменения концентрации С* (т) в процессе разгазирования при постоянной максимальной скорости выталкивания Qn = (?Птах всегда больше, чем при соблюдении условия (1.28), при этом зависимость

с ; « = / ( * г,|»п_вп_ = / ( | ^ )

(1.32)

можно рассматривать как оценку сверху истинного значения функции Ст (т), соответствующей решению уравнения (1.31).

Скорость изменения во времени коэффициента усиления объекта по каналу «газовыделение — концентрация», как видно из (1.24), равна скорости изменения во времени параметра Ст, являющегося, следовательно, важным параметром идентификации объекта.

29

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