Можно вместо точного оптимального профиля рассчитать квазиоптимальный, в котором программа изменения во времени режимного параметра (температуры) заменена ступенчатой (рис. 10.7). В этой программе в начальный период ферментации параметр поддерживается оптимальным для роста биомассы, а затем происходит переключение (момент tп) на значение, оптимальное для биосинтеза, которое и поддерживается до конца ферментации tк:
(10 4)
Рис. 10.7. Примерный вид оптимального профиля квазиоптимальной ступенчатой программы
Время tп переключения с одного участка на другой в такой программе можно определить, если известны значения μ:
μ = μX при T = TX (максимальная скорость роста); (10.5)
μ = μP при T= TP (медленный рост до конца процесса). (10.6)
Кроме того, необходимо знать ограничение Xmax, определяемое характеристиками аппарата. Очевидно, что при t= tK
X ≤ Xmax, (10.7)
В первом периоде ферментации (при t < tn)
X = X0eμXt (10.8)
В точке переключения (при t = tп)
X = X0eμXtgп (10.9)
В дальнейшем (при t > tп) рост биомассы подчиняется выражению
X = [ X0eµXtп]eµ(t–tп) (10.10)
При t = tK
X = Xmax (10.11)
Xmax = (X0eµXtп)eµP(tK–tп) (10.12)
откуда
(10.13)
Аналогичным образом можно рассчитать оптимальную ступенчатую программу для величины рН, для концентрации субстрата и других параметров. Надо только знать комплекс значений всех факторов, оптимальный для роста, и другой комплекс – оптимальный для биосинтеза, и соответствующие удельные скорости роста μX и μp>.
Приведем практический пример расчета квазиоптимальной программы управления по величине рН. На рис. 10.8 представлена зависимость удельной скорости роста μ и удельной скорости биосинтеза антибиотика qP от рН.
На рис. 10.8 можно определить, что оптимальные значения pH для роста биомассы рНX = 4,8, а для биосинтеза антибиотика рНР = 7,0. При этом в условиях, оптимальных для роста биомассы (рН 4,8), удельная скорость роста μX = 0,03 ч–1, а в условиях, оптимальных для биосинтеза антибиотика (рН 7,0), μΡ = 0,01 ч–1. Начальная концентрация биомассы в аппарате X0 = 2 г/л, максимально допустимая для данного аппарата Xmax = 40 г/л.
Рис. 10.8. Пример зависимости роста микроорганизмов и биосинтеза антибиотика от величины рН
Задано также время ферментации tK = 200 ч. Цифровая ось указана только для μ. Ось для величины qp не представлена, поскольку цифровые значения qP не играют роли для определения параметров ступенчатой программы управления.
Необходимо определить ступенчатую квазиоптимальную программу управления.
Рис. 10.9. Расчетная ступенчатая программа управления процессом ферментации для примера, рассмотренного в разделе 10.3
Имея эти значения, можно найти точку переключения квазиоптимальной программы tп:
Величина tK максимальна при tп = 0:
Если с начала до конца рост идет при рН = рНX = 4,8, то tK = 100 ч, а если при рН = рНP = 7,0, то tK = 300 ч.
При заданном времени ферментации tK = 200 ч профиль рН соответствует показанному на рис. 10.9.