СПОСОБ ПРИГОТОВЛЕНИЯ БИОКОМПОСТА ИЗ ОРГАНИЧЕСКИХ ОТХОДОВ И УСТРОЙСТВО ДЛЯ ЕГО РЕАЛИЗАЦИИ

1. Способ приготовления биокомпоста из органических отходов, включающий в себя смешивание птичьего помета или навоза с органическим наполнителем до получения сыпучей массы, компостирование полученной смеси в ферментере при принудительном вентилировании объемами воздуха, снижающимися по мере ускорения аэробно-анаэробного процесса ферментации с одновременным контролем температуры и влажности, отличающийся тем, что в нем осуществляют два режима вентиляции, в первом режиме, с переменным расходом, предназначенном для разогрева ферментируемой биомассы, воздух из вытяжной системы ферментера подают на воздуховоздушный рекуператор, через который пропускают и приточный воздух, во втором режиме с постоянным максимальным расходом, приточный воздух подают напрямую в ферментируемую биомассу, а воздух из вытяжной системы подают напрямую в окружающую среду, при этом переход из первого режима во второй осуществляют по достижению температуры ферментируемой биомассы 65-70°C, периодически, через промежутки времени, равные 1-2 часа, измеряют температуру, влажность и содержание кислорода в приточном воздуха и в вытяжном окне системы вентиляции, посредством математической модели оценивают состояние ферментируемой биомассы для каждого момента времени измерения, по полученным оценкам формируют оптимальные программы изменения расхода воздуха, обеспечивающего минимальное время разогрева ферментируемой биомассы.
2. Способ по п.1, отличающийся тем, что математическая модель, посредством которой оценивают состояния ферментируемой биомассы, имеет следующий вид:
<mrow><mover accent="true"><mi>X</mi><mo>˙</mo></mover><mo>=</mo><mi>Ф</mi><mrow><mo>(</mo><mrow><mi>X</mi><mo>,</mo><mtext> </mtext><mtext>Fu</mtext><mrow><mo>(</mo><mtext>t</mtext><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mo>,</mo><mtext>    (1)</mtext></mrow>
<mrow><mover accent="true"><mi>X</mi><mo>˙</mo></mover><mo>=</mo><mfrac bevelled="true"><mrow><mo>∂</mo><mi>X</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac></mrow>
где
- производная вектора состояния по времени,
Ф(X,Fu(t)) - оператор модели, имеющий следующий вид
<mrow><mi>Ф</mi><mrow><mo>(</mo><mrow><mi>X</mi><mo>,</mo><mi>F</mi><mo>,</mo><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mo>=</mo><mrow><mo>[</mo><mtable columnalign="left"><mtr><mtd><msub><mi>a</mi><mrow><mn>11</mn></mrow></msub><msub><mi>x</mi><mn>1</mn></msub><mo>+</mo><msub><mi>a</mi><mrow><mn>12</mn></mrow></msub><msub><mi>x</mi><mn>2</mn></msub><mo>+</mo><msub><mi>a</mi><mrow><mn>13</mn></mrow></msub><msub><mi>x</mi><mn>3</mn></msub><mo>+</mo><msub><mi>a</mi><mrow><mn>15</mn></mrow></msub><msub><mi>x</mi><mn>5</mn></msub><mo>+</mo><msub><mi>b</mi><mrow><mn>11</mn></mrow></msub><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mrow><mo>(</mo><mrow><msub><mi>f</mi><mrow><mn>12</mn></mrow></msub><mo>−</mo><msub><mi>x</mi><mn>1</mn></msub></mrow><mo>)</mo></mrow><mo>+</mo><mi>ϕ</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>a</mi><mrow><mn>21</mn></mrow></msub><msub><mi>x</mi><mn>1</mn></msub><mo>+</mo><msub><mi>a</mi><mrow><mn>22</mn></mrow></msub><msub><mi>x</mi><mn>2</mn></msub><mo>+</mo><msub><mi>a</mi><mrow><mn>24</mn></mrow></msub><msub><mi>x</mi><mn>4</mn></msub><mo>+</mo><msub><mi>c</mi><mrow><mn>21</mn></mrow></msub><mrow><mo>(</mo><mrow><msub><mi>x</mi><mn>2</mn></msub><mo>−</mo><msub><mi>f</mi><mn>1</mn></msub></mrow><mo>)</mo></mrow><mo>+</mo><msub><mi>b</mi><mrow><mn>21</mn></mrow></msub><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mrow><mo>(</mo><mrow><msub><mi>x</mi><mn>1</mn></msub><mo>−</mo><msub><mi>f</mi><mrow><mn>12</mn></mrow></msub></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>a</mi><mrow><mn>31</mn></mrow></msub><msub><mi>x</mi><mn>1</mn></msub><mo>+</mo><msub><mi>a</mi><mrow><mn>33</mn></mrow></msub><msub><mi>x</mi><mn>3</mn></msub><mo>+</mo><msub><mi>a</mi><mrow><mn>34</mn></mrow></msub><msub><mi>x</mi><mn>4</mn></msub><mo>+</mo><msub><mi>b</mi><mrow><mn>31</mn></mrow></msub><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mrow><mo>(</mo><mrow><msub><mi>x</mi><mn>3</mn></msub><mo>−</mo><msub><mi>f</mi><mn>2</mn></msub></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>a</mi><mrow><mn>42</mn></mrow></msub><msub><mi>x</mi><mn>2</mn></msub><mo>+</mo><msub><mi>a</mi><mrow><mn>43</mn></mrow></msub><msub><mi>x</mi><mn>3</mn></msub><mo>+</mo><msub><mi>a</mi><mrow><mn>44</mn></mrow></msub><msub><mi>x</mi><mn>4</mn></msub><mo>+</mo><msub><mi>a</mi><mrow><mn>45</mn></mrow></msub><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mrow><mo>(</mo><mrow><msub><mi>x</mi><mn>3</mn></msub><mo>−</mo><msub><mi>x</mi><mn>4</mn></msub></mrow><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mrow><mn>41</mn></mrow></msub><mrow><mo>(</mo><mrow><msub><mi>x</mi><mn>4</mn></msub><mo>−</mo><msub><mi>f</mi><mn>2</mn></msub></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>a</mi><mrow><mn>51</mn></mrow></msub><msub><mi>x</mi><mn>1</mn></msub><mo>+</mo><msub><mi>a</mi><mrow><mn>55</mn></mrow></msub><msub><mi>x</mi><mn>5</mn></msub><mo>+</mo><msub><mi>b</mi><mrow><mn>51</mn></mrow></msub><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mrow><mo>(</mo><mrow><msub><mi>f</mi><mn>3</mn></msub><mo>−</mo><msub><mi>x</mi><mn>5</mn></msub></mrow><mo>)</mo></mrow></mtd></mtr></mtable><mo>]</mo></mrow><mo>,</mo><mtext>    (2)</mtext></mrow>
где:
<mrow><msub><mover accent="true"><mi>x</mi><mo>˙</mo></mover><mn>1</mn></msub></mrow>
1
2
3
5
12
1
1
2
2
в первой строке оператора модели отражена взаимосвязь скорости изменения температуры ферментируемой массы
от самой температуры x от температуры воздушной среды в ферментере x, от влажности ферментируемой биомассы x, от содержания кислорода в воздушной среде ферментера x, от потери теплоты с вентиляционным воздухом в режиме нагрева, пропорциональной произведению расхода вентиляционного воздуха u(t) на разность между температурой приточного воздуха после рекуператора f и температурой наружного воздуха f, от теплоты реакции ферментации, задаваемой моделью φ(t)-dt+dt;
<mrow><msub><mover accent="true"><mi>x</mi><mo>˙</mo></mover><mn>2</mn></msub></mrow>
2
1
4
2
1
1
12
во второй строке оператора модели отражена зависимость скорости изменения температуры воздушной среды в ферментере
от самой температуры x, от температуры ферментируемой биомассы x от относительной влажности воздушной среды в ферментере x, от потерь теплоты сквозь ограждающую конструкцию ферментера, пропорциональную разности между температурой воздушной среды в ферментере x и температурой наружного воздуха f, от притока теплоты за счет реакции ферментации, пропорциональной произведению расхода вентиляционного воздуха u(t) на разность между температурой ферментируемой биомассы x и температурой приточного воздуха после рекуператора f;
<mrow><msub><mover accent="true"><mi>x</mi><mo>˙</mo></mover><mn>3</mn></msub></mrow>
3
1
4
3
2
в третьей строке оператора модели отражена зависимость скорость изменения влажности ферментируемой биомассы
от самой влажности x, от ее температуры x от относительной влажности воздушной среды в ферментере x и от потери влаги за счет вентилирования, пропорциональной произведению расхода вентиляционного воздуха u(t) на разность между влажностью ферментируемой биомассы x и относительной влажностью наружного воздуха f;
<mrow><msub><mover accent="true"><mi>x</mi><mo>˙</mo></mover><mn>4</mn></msub></mrow>
4
2
3
3
4
4
2
в четвертой строке оператора модели отражена зависимость скорости изменения относительной влажности воздушной среды в ферментере
от самой относительной влажности x, от температуры воздушной среды в ферментере x, от влажности ферментируемой биомассы x, от притока влаги от ферментируемой биомассы, пропорционального произведению расхода вентиляционного воздуха u(t) на разность между влажностью ферментируемой биомассы x и относительной влажностью воздушной среды в ферментере x, от потерь влаги, пропорциональной произведению расхода вентиляционного воздуха u(t) на разность между относительной влажностью воздушной среды в ферментере x и относительной влажностью наружного воздуха f;
<mrow><msub><mover accent="true"><mi>x</mi><mo>˙</mo></mover><mn>5</mn></msub></mrow>
5
1
3
5
в пятой сроке оператора модели отражена зависимость скорости изменения содержания кислорода в воздушной среде ферментера
от самого содержания кислорода x, от температуры компостируемой биомассы x и от притока кислорода, пропорционального произведению расхода вентиляционного воздуха u(t) на разность между содержанием кислорода в окружающем воздухе f и в воздушной среде ферментера x;
Tr
-1
1
2
3
4
5
1
2
3
4
5
X=[x x x x x] - вектор состояния ферментируемой биомассы, представленный в транспонированной форме, компонентами которого являются: температура компостируемой массы - x, °C; температура воздушной среды в компостере - x, °C; влажность компостируемой массы - x, %; влажность воздушной среды в компостере - x, %; содержание кислорода в воздушной среде ферментера - x, г кг; Tr - индекс транспонирования векторов и матриц;
Tr
-1
1
2
3
1
2
3
F=[f f f] - вектор параметров окружающей среды, компонентами которого являются: температура наружного воздуха - f, °C; относительная влажность наружного воздуха - f, %; содержание кислорода в наружном воздухе - f, г кг;
3
-1
u(t)- расход вентиляционного воздуха, м час;
t - переменная времени, час;
1
2
2
-1
φ(t)-dt+dt - функция выделения теплоты при ферментации, (°С);
12
1
2
f=f+k(u)x - функция рекуперации теплоты, °C;
11-
55
11-
51
41
1
2-
-1
a a, b b, c, d, d параметры модели, имеющие размерность (час) и объединяемые в вектор P, который уточняют после периодических измерений состояния ферментируемой биомассы.
3. Способ по п.2, отличающийся тем, что оценку состояния ферментируемой биомассы формируют по следующей последовательности операций, выполняемых в моменты измерений t=1, 2, 3…
<mrow><mover accent="true"><mover accent="true"><mi>X</mi><mo>⌢</mo></mover><mo>˙</mo></mover><mrow><mo>(</mo><mrow><mi>t</mi><mrow><mo>|</mo><mrow><mi>t</mi><mo>−</mo><mn>1</mn></mrow></mrow></mrow><mo>)</mo></mrow><mo>=</mo><mi>Ф</mi><mrow><mo>(</mo><mrow><mover accent="true"><mi>X</mi><mo>⌢</mo></mover><mo>,</mo><mi>F</mi><mo>,</mo><mi>u</mi><mo>,</mo><mi>t</mi></mrow><mo>)</mo></mrow><mo>,</mo></mrow>
<mrow><mover accent="true"><mi>P</mi><mo>˙</mo></mover><mrow><mo>(</mo><mrow><mi>t</mi><mrow><mo>|</mo><mrow><mi>t</mi><mo>−</mo><mn>1</mn></mrow></mrow></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>∂</mo><mi>Ф</mi><mrow><mo>(</mo><mrow><mover accent="true"><mi>X</mi><mo>⌢</mo></mover><mo>,</mo><mi>F</mi><mo>,</mo><mi>u</mi><mo>,</mo><mi>t</mi></mrow><mo>)</mo></mrow></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mi>P</mi><mrow><mo>(</mo><mrow><mi>t</mi><mo>−</mo><mn>1</mn><mrow><mo>|</mo><mrow><mi>t</mi><mo>−</mo><mn>1</mn></mrow></mrow></mrow><mo>)</mo></mrow><mo>+</mo><mi>P</mi><mrow><mo>(</mo><mrow><mi>t</mi><mo>−</mo><mn>1</mn><mrow><mo>|</mo><mrow><mi>t</mi><mo>−</mo><mn>1</mn></mrow></mrow></mrow><mo>)</mo></mrow><msup><mrow><mo stretchy="false">(</mo><mfrac><mrow><mo>∂</mo><mi>Ф</mi><mrow><mo>(</mo><mrow><mover accent="true"><mi>X</mi><mo>⌢</mo></mover><mo>,</mo><mi>F</mi><mo>,</mo><mi>u</mi><mo>,</mo><mi>t</mi></mrow><mo>)</mo></mrow></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo stretchy="false">)</mo></mrow><mrow><mi>T</mi><mi>r</mi></mrow></msup><mo>,</mo></mrow>
<mrow><mover accent="true"><mi>X</mi><mo>⌢</mo></mover><mrow><mo>(</mo><mrow><mi>t</mi><mrow><mo>|</mo><mi>t</mi></mrow></mrow><mo>)</mo></mrow><mo>=</mo><mover accent="true"><mi>X</mi><mo>⌢</mo></mover><mrow><mo>(</mo><mrow><mi>t</mi><mrow><mo>|</mo><mrow><mi>t</mi><mo>−</mo><mn>1</mn></mrow></mrow></mrow><mo>)</mo></mrow><mo>+</mo><mi>K</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mrow><mo>[</mo><mrow><mi>Y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>−</mo><mi>H</mi><mover accent="true"><mi>X</mi><mo>⌢</mo></mover><mrow><mo>(</mo><mrow><mi>t</mi><mrow><mo>|</mo><mrow><mi>t</mi><mo>−</mo><mn>1</mn></mrow></mrow></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow><mo>,</mo></mrow>
<mrow><mi>P</mi><mrow><mo>(</mo><mrow><mi>t</mi><mrow><mo>|</mo><mi>t</mi></mrow></mrow><mo>)</mo></mrow><mo>=</mo><mi>P</mi><mrow><mo>(</mo><mrow><mi>t</mi><mrow><mo>|</mo><mrow><mi>t</mi><mo>−</mo><mn>1</mn></mrow></mrow></mrow><mo>)</mo></mrow><mo>−</mo><mi>K</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msup><mi>H</mi><mrow><mi>T</mi><mi>r</mi></mrow></msup><mi>P</mi><mrow><mo>(</mo><mrow><mi>t</mi><mrow><mo>|</mo><mrow><mi>t</mi><mo>−</mo><mn>1</mn></mrow></mrow></mrow><mo>)</mo></mrow><mo>,</mo><mtext>                                                         (3)</mtext></mrow>
<mrow><mi>K</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>P</mi><mrow><mo>(</mo><mrow><mi>t</mi><mrow><mo>|</mo><mrow><mi>t</mi><mo>−</mo><mn>1</mn></mrow></mrow></mrow><mo>)</mo></mrow><mi>H</mi><msup><mrow><mrow><mo>[</mo><mrow><mi>H</mi><mi>P</mi><mrow><mo>(</mo><mrow><mi>t</mi><mo>−</mo><mn>1</mn></mrow><mo>)</mo></mrow><msup><mi>H</mi><mrow><mi>T</mi><mi>r</mi></mrow></msup><mo>+</mo><msup><mi>Q</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow><mo>]</mo></mrow></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo>,</mo></mrow>
Tr
2
4
5
2
4
5
<mrow><mover accent="true"><mi>X</mi><mo>⌢</mo></mover><mrow><mo>(</mo><mrow><mi>t</mi><mrow><mo>|</mo><mrow><mi>t</mi><mo>−</mo><mn>1</mn></mrow></mrow></mrow><mo>)</mo></mrow><mo>,</mo><mtext> </mtext><mover accent="true"><mi>X</mi><mo>⌢</mo></mover><mrow><mo>(</mo><mrow><mi>t</mi><mrow><mo>|</mo><mi>t</mi></mrow></mrow><mo>)</mo></mrow></mrow>
где: Y=[x x x] - вектор измеряемых состояний, включающий в себя температуру - x, влажность - x и содержание кислорода - x в воздухе перед вытяжным окном; (t|t) - слева от вертикальной черты обозначено прогнозное время, справа - время измерения и оценивания;
- прогнозная и уточненная оценки вектора состояния; К(t)-[5×3] -матрица коэффициентов сглаживания; H=[0 1 0 1 1] - матрица структуры измерения; Р(t|t), P(t|t-1)-[5×5] - прогнозная и уточненная матрица ошибок оценивания; Q-[3×3] - матрица ошибок измерения.
4. Способ по п.2, отличающийся тем, что оптимальную программу изменения расхода вентиляционного воздуха минимизирующую время разогрева биомассы формируют в моменты измерений по следующей последовательности операций:
<mrow><msub><mi>u</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>u</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>−</mo><msub><mi>Δ</mi><mi>i</mi></msub><msubsup><mi>λ</mi><mi>i</mi><mrow><mi>T</mi><mi>r</mi></mrow></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mfrac><mrow><mo>∂</mo><msub><mi>Φ</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mrow><mo>(</mo><mrow><mi>Χ</mi><mo>,</mo><mi>F</mi><mo>,</mo><msub><mi>u</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow><mrow><mo>∂</mo><mi>u</mi></mrow></mfrac><mo>,</mo><mtext>   (4)</mtext></mrow>
<mrow><mover accent="true"><mi>λ</mi><mo>˙</mo></mover><mo>=</mo><mo>−</mo><msup><mi>λ</mi><mrow><mi>T</mi><mi>r</mi></mrow></msup><mfrac><mrow><mo>∂</mo><mrow><mo>{</mo><mrow><msub><mi>Φ</mi><mn>1</mn></msub><mrow><mo>(</mo><mrow><mi>Χ</mi><mo>,</mo><mi>F</mi><mo>,</mo><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow><mo>}</mo></mrow></mrow><mrow><mo>∂</mo><mi>Χ</mi></mrow></mfrac><mo>,</mo><mtext>                           (5)</mtext></mrow>
<mrow><msub><mi>T</mi><mrow><mn>1</mn><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>T</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>−</mo><msub><mi>Δ</mi><mi>i</mi></msub><mfrac><mrow><mo>∂</mo><msub><mi>Φ</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mrow><mo>(</mo><mrow><mi>Χ</mi><mo>,</mo><mi>F</mi><mo>,</mo><msub><mi>u</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>,</mo><mtext>            (6)</mtext></mrow>
11
1
11
1
g=1, если (x-x)>0, g=0, если (x-x)<0,
11
1
i
где: i - индекс итераций формирования программы; λ - вектор вспомогательных переменных одинаковой размерности с вектором состояния биомассы X; x - заданное значение температуры разогрева биомассы (65-70°C); T - время разогрева, являющееся варьируемой переменной (фактором оптимизации); Δ - шаг оптимизации, g=0,6-0,75 - управляющий параметр алгоритма.
5. Способ по п.2, отличающийся тем, что вектор параметров P математической модели уточняют в моменты измерений по следующей последовательности операций:
<mrow><msub><mi>P</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>P</mi><mn>1</mn></msub><msub><mrow /><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>−</mo><msub><mi>δ</mi><mi>i</mi></msub><mfrac><mn>1</mn><mi>t</mi></mfrac><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mn>0</mn><mi>t</mi></munderover><mrow><mrow><mo>(</mo><mrow><mi>Y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>−</mo><mover accent="true"><mi>Y</mi><mo>⌢</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mfrac><mrow><mo>∂</mo><msub><mi>Φ</mi><mi>i</mi></msub><mrow><mo>(</mo><mrow><mi>Χ</mi><mo>,</mo><mi>F</mi><mo>,</mo><msub><mi>u</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow><mrow><mo>∂</mo><mi>P</mi></mrow></mfrac><mi>d</mi><mi>t</mi><mo>,</mo><mtext>          (7)</mtext></mrow></mrow></mstyle></mrow>
i
<mrow><mover accent="true"><mi>Y</mi><mo>⌢</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow>
где: i - индекс операции, δ - рабочий шаг операции,
- моделируемое значение вектора измеряемых состояний ферментируемой биомассы.
6. Устройство для реализации способа по п.1, содержащее ферментер с перфорированным днищем, систему вентиляции с нагнетательным вентилятором для подачи воздуха снизу вверх через перфорированное днище во внутреннее пространство ферментера, систему управления работой вентилятора с термодатчиком, отличающееся тем, что система вентиляции содержит воздуховоздушный рекуператор с перепускным тройником на входе и управляемой заслонкой, изменяющей направления потока, обеспечивающей пропуск воздуха из внутреннего пространства ферментера через рекуператор и в обход его, нагнетательный вентилятор переменной производительности, установленный между рекуператором и перфорированным днищем ферментера, датчики температуры, влажности и содержания кислорода, установленные на входе и выходе системы вентиляции, и автоматическое управляющее устройство, на вход которого подключены датчики температуры, влажности и содержания кислорода, а выход которого соединен с управляемым приводом нагнетательного вентилятора и приводом управляемой заслонки.