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

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

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

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

подобия,

или в реальном

смесителе.

Первый способ точнее, но

осуществить

его

на стадии

проектирования

сооружений

далеко

не всегда

возмож­

но. Остановимся

на втором

способе

как

наиболее до­

ступном в практике.

 

 

 

 

 

Значение D определяется решением обратной задачи

из уравнения (19).

 

 

 

 

 

В реальном

ершовом смесителе

экспериментально

получают несколько значений установившихся

концент­

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

нении (19), должны быть стабилизированы. По

этим

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

распределения

концентрации вещества по длине смесителя.

 

Записав уравнение (19) в логарифмической форме,

получим уравнение прямой:

 

 

lnC(oo,*) = lnC1 +

ftx,

(20)

где

 

 

ь = - У + + у > к - ^ ^

2D

( 2 1 )

v ;

Отсюда

D =

£ ± ± - .

 

(22)

Эта прямая строится по экспериментальным значе­

ниям In С(оо, х). Для х=0

In С(оо, х)=\п

Сь

Величина

Ь равна тангенсу угла наклона прямой к оси х.

 

Подставив значение b

в выражение

(21),

находим

величину D.

 

 

 

Степень точности данного метода определения коэф­ фициента турбулентной диффузии позволяет пользо­ ваться им при расчете устойчивости САР.

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

Для вывода передаточной функции смесителя-реак­ тора задаемся длиной х\ и получаем однозначное реше­ ние уравнения (15):

_

Уо» + 4Р-(К + р)

—у

 

С в ы х ( р 1 , х 1 ) = С 1 е

2 0

\

(23)

40

Передаточная функция на участке от входа в смеси­ тель-реактор до выхода из него

W(p1,x1)

= c * i a

(Pl,Xl>

 

= е

w

. (24)

Рассмотренные

в этом

разделе

расчетные

способы

получения

статических

и

динамических характеристик

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

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

5. Расчет настроечных параметров и устойчивости систем автоматического регулирования

При выборе типа регулятора приходится искать ком­

промиссное решение между

стремлением упростить

и, следовательно, удешевить

систему регулирования

и необходимостью обеспечить требуемое качество регу­ лирования.

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

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

41

Яд = °р/°о,

(25)

где 0О — отклонение параметра при том же

возмуще­

нии и отключенном регуляторе.

 

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

саны в работе

[28].

 

 

Общей динамической характеристикой объектов ре­

гулирования,

как сказано выше,

служит

отношение

х / Т 0 . Тип регулятора определяется

значением

этого от­

ношения:

 

 

 

Позиционный

< 0 , 2

Непрерывный

0,2—1,5

 

Импульсный

> 1,5

При небольших величинах т область применения не­ прерывных регуляторов расширяется.

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

= (Р макс i °макс f+г) (26)

О

макс i

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

Показателем качества регулирования в установив­ шемся режиме является величина статической неравно­ мерности. Остающееся отклонение параметра свойствен­

но П-регуляторам

и в меньшей

степени

ПИ-регулято-

рам. В последних

возможно

длительное

отклонение

параметра регулирования, вызванное характером возму­ щений (см. п. 2 данной главы).

Параметрами настройки непрерывных регуляторов служат: степень (скорость) обратной связи б для П-ре- гулятора, степень обратной связи и время изодрома Г и для ПИ-регулятора; степень обратной связи, время изо­

дрома

и время предварения

Гпр для

ПИД-регулятора;

степень

обратной

связи

и

время

предварения для

ПД-регулятора.

 

 

 

 

 

В табл. 2 приводятся

рекомендуемые

Всесоюзным

теплотехническим

институтом

им. Ф. Э.

Дзержинского

42

Расчетные формулы настройки регуляторов

Т/Г.

п

|

ПИ

 

 

6

 

б

 

<0, 2

ет

 

1,1ет

 

2,6

— —0,08

2,6

— —0,08

0,2—1,5

Г„

 

Г 0

Р

т

 

Р

т

 

+ 0 , 7

 

 

т

 

- ^ + 0 , 6

 

 

 

 

 

 

 

0

 

' 0

Регулятор

 

|

п и д

Г

б

и

 

3,3т

0,8ет

0,8

т

3.7

Т Г - ° ' 1 3

т

Р

t

 

"т^Г

 

 

• ^ - + 1 . 5

 

 

 

>1,5

2

2

0,6т

1,7

Р

Р

 

 

 

Р

Т а б л и ц а 2

тт

ипр

2,5т (0,75-=-2)т

1

(0,75-=-1,5)т

— т

т

 

" Т о "

0,7т (0,75^-2) т

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

рывных

регуляторов по экспериментально полученным

параметрам объектов

регулирования.

 

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

торов составлены для

получения переходного процесса

со степенью затухания яр = 0,75.

 

Анализ работы САР

в переходном

режиме и опреде­

ление ее

устойчивости

осуществляются

с помощью час­

тотных характеристик регулируемого объекта и всей системы. Частотные характеристики отражают поведе­

ние системы при

действии

на

вход в нее

возмущения

в форме непрерывных

гармонических

колебаний.

 

Отношение амплитуд

выходных

и

входных

колеба­

ний

в зависимости

от частоты

Л(ю)

называют

ампли­

тудно-частотной

характеристикой:

 

 

 

 

 

 

 

А((о) = Ааьа.

 

 

 

 

(27)

Зависимость разности фаз этих колебаний от часто­

ты ф(со) называют

фазо-частотной характеристикой:

 

 

 

Ф Н

= фв ы х — фв х -

 

 

(28)

 

Комплексную

 

функцию

частоты

W(ia)

называют

амплитудно-фазовой

характеристикой

системы

(АФХ):

 

 

W (tea)

= А (со) е 1 ф ( ш )

,

 

 

(29)

где

Л(со)— модуль

комплексной

функции;

изображает­

 

ся в

виде

вектора длиной,

равной

отноше­

 

нию амплитуд.

 

 

 

 

 

 

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

На кривой разгона (см. рис. 8) по оси ординат отло­ жены значения х в процентах от предельного значения в пересчете на процентное возмущение. Для построения АФХ кривая разгона аппроксимируется несколькими вертикальными прямолинейными отрезками. Каждый отрезок Хг можно рассматривать как кривую разгона элементарного звена, обладающего только транспорт­ ным запаздыванием, равным 30 сек. Весь объект пред­ ставляет тогда систему, состоящую из параллельных элементарных звеньев.

44

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

дого звена будут наблюдаться

гармонические колеба­

ния со сдвигом фаз Atpi = Атт ,

Амплитуда и частота

1т М

Рис. 9. Амплитудно-фазовая характеристика ершового смеси­ теля и САР

этих колебаний не отличаются от входных, так как кро­ ме запаздывания звенья не имеют других свойств. Таким образом, выходную функцию каждого звена можно изобразить на комплексной плоскости вектором ri(xt, Афг). Выходная функция всего объекта изобразится суммой этих векторов 2>j. Так получают одну точку го­

дографа вектора

R((ui)

для частоты юь. Чтобы

построить

всю

АФХ объекта,

необходимо

нанести

не

менее

15—20 точек.

 

 

 

 

 

 

На

рис. 9 показана

построенная

изложенным

спосо­

бом АФХ ершового смесителя

(см. рис. 8). Для частоты

ю = 0,015 рад/сек

дан

порядок

суммирования

векторов

г,. Каждый последующий вектор откладывается от кон­

ца

предыдущего,

начиная с положительного

направле­

ния

оси абсцисс,

под

углом аи — 57° Atcah- Здесь

а =

= 57-30-0,015 = 25,6°.

Необходимо отметить,

что

точ­

ность определения модуля вектора изложенным спосо­ бом падает с увеличением частоты. Ошибка может до­ стигнуть 30%- Слева от годографа объекта нанесены АФХ системы регулирования с ПИ-регулятором для разных значений времени изодрома Ти. АФХ системы

45

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

Исследования показали, что для расчета САР про­ цессов реагентной очистки сточных вод применим кри­ терий Найквиста—Михайлова, позволяющий опреде­ лять устойчивость замкнутой системы регулирования, которая обеспечивается устойчивостью разомкнутой системы (без обратной связи) при условии, что ее АФХ не охватывает в комплексной плоскости точку с коорди­

натами

(—1, 0).

 

 

 

 

 

 

АФХ

одноконтурной

САР выражается

произведени­

ем АФХ объекта и регулятора:

 

 

 

 

 

 

 

W(m) =

Wo6(i«>)Wp(m).

 

 

 

(30)

 

Построение области устойчивости такой САР рас­

смотрим на примере ершового смесителя

(см. рис. 9).

АФХ объекта

пересекает

вещественную

полуось в

точке

с

координатой —2. Подключение к

такому

объекту

П-регулятора

изменяет

длину вектора

в kv

раз

(где

kv

— коэффициент передачи регулятора). Для

обеспече­

ния устойчивости САР в данном случае

kv,

очевидно,

должен быть не более —1/—2=0,5. В

результате

воз­

действия

астатической

(интегральной)

составляющей

каждый вектор АФХ поворачивается на 90° и изменяет­ ся в ер/со раз, что приближает годограф к критической точке.

Построим на этом же графике АФХ разомкнутой системы с ПИ-регулятором. Для этого выберем на го­

дографе объекта

частоту

со^ и к концу соответствующе­

го ей вектора A0k

под углом 90° прибавим вектор ААи-

Длина его подсчитывается

по формуле

 

 

Ш£ I и

Соединив начало координат с концом дополнитель­ ного вектора, получим вектор АФХ разомкнутой систе­ мы Ak для данной частоты щ. Таким же образом строим АФХ разомкнутой системы и для других частот, раз­ личных при Ти и КР=1. Частоты при этом нужно выби­ рать так, чтобы получать отрезки годографов вблизи отрицательной вещественной полуоси (интерес пред­ ставляют лишь точки их пересечения с последней).

Значение предельного (критического) коэффициента

46

передачи регулятора

для

каждого Ги определяется

по

формуле

 

 

 

 

 

 

 

*р.кр=1'#*.

(32)

где Rk—отрезок

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

вещественной полуоси,

отсекаемый

АФХ

разомкнутой системы.

 

Полученные данные позволяют графически предста­

вить область устойчивости

САР

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

настройки регулятора

Та и /гр .к р

(рис. 1 0 ) .

 

'П11

1

[

- х -

i-

-—i

~ J

1

' ——ri

'20

w

so

su

wo

/го

m

tso

us TU)cen

 

 

Рис.

10. Область устойчивости

CAP

 

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

АЛ, = - ^

KpA0i^Ttt,

(33)

СО,- I и

где knp— отношение времени предварения Тир к времени изодрома Ги .

В плоскости параметров настройки граница устойчи­ вости с ПИД-регулятором располагается выше кривой, нанесенной на рис. 10, вследствие уменьшения угла по­ ворота вектора АФХ под действием дифференциальной составляющей.

47

Г л а в а I I I

АВТОМАТИЧЕСКОЕ Р Е Г У Л И Р О В А Н И Е ПРОЦЕССОВ

НЕ Й Т Р А Л И З А Ц И И СТОЧНЫХ ВОД

1.Регулирование очистки сточных вод, содержащих преимущественно кислоты

При автоматическом регулировании процессов ней­ трализации, очистки от металлов и ряда других процес­ сов химической очистки сточных вод для лучшей стаби­ лизации регулируемого параметра большое значение имеет предварительное усреднение стоков по составу. САР на очистных сооружениях с предварительным ус­ реднением стоков не сложны по структуре и строятся по принципу регулирования параметра рН. Такие САР до­ статочно полно описаны в ряде работ [60, 61]. Здесь мы рассмотрим наиболее сложные случаи регулирования величины рН.

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

рованному принципу

(рис. 11);

 

 

Регулятор системы

получает

на

вход три сигнала:

два непрерывных

(от рН-метра

7 исходной воды и от до­

затора реагента)

и один дискретный

(от рН-метра 10 на

выходе из смесителя), сигнализирующий лишь о пере­ ходе величины рН через заданное значение. При паде­ нии величины рН в конце смесителя ниже этого значе­ ния открывается регулирующий орган (нож дозатора) на заранее заданную величину; при увеличении значе­ ния рН нож дозатора возвращается в исходное положе­ ние. Ход ножа (размах) составляет в зависимости от настройки 30—60% полного хода.

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

48

Рис. 11. Принципиальная

технологическая схема

пропорционально-дискретной

САР

/ — автоматический потенциометр контрольного рН-метра;

2 — высокоомные преобразователи

рН-мет-

ра;

3 —датчики

рН-метра; 4 — ершовый смеситель-реактор; 5 — исполнительный механизм с реостат­

ным

датчиком

обратной связи;

6 дозатор реагента; 7 автоматический потенциометр рН-метра с

функциональным вторичным датчиком; 8

— регулятор; 9 — резисторы настройки; 10 автоматический

 

 

потенциометр

рН-метра

с позиционным

контактным устройством

 

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