![](/user_photo/_userpic.png)
книги из ГПНТБ / Смирнов, Д. Н. Автоматическое регулирование процессов очистки сточных и природных вод
.pdfродинамического |
подобия, |
или в реальном |
смесителе. |
|||
Первый способ точнее, но |
осуществить |
его |
на стадии |
|||
проектирования |
сооружений |
далеко |
не всегда |
возмож |
||
но. Остановимся |
на втором |
способе |
как |
наиболее до |
||
ступном в практике. |
|
|
|
|
|
|
Значение 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 |
' ——r—i |
'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 — автоматический |
|||||
|
|
потенциометр |
рН-метра |
с позиционным |
контактным устройством |
|