CN103810391B - 运动发酵单胞菌生产乙醇中振荡过程的调控方法 - Google Patents
运动发酵单胞菌生产乙醇中振荡过程的调控方法 Download PDFInfo
- Publication number
- CN103810391B CN103810391B CN201410059227.3A CN201410059227A CN103810391B CN 103810391 B CN103810391 B CN 103810391B CN 201410059227 A CN201410059227 A CN 201410059227A CN 103810391 B CN103810391 B CN 103810391B
- Authority
- CN
- China
- Prior art keywords
- oscillatory
- zymomonas mobilis
- data
- regulation
- production
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 73
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 title claims abstract description 50
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 49
- 241000588902 Zymomonas mobilis Species 0.000 title claims abstract description 25
- 230000001105 regulatory effect Effects 0.000 title claims abstract description 12
- 230000001276 controlling effect Effects 0.000 title abstract description 4
- 235000019441 ethanol Nutrition 0.000 title abstract 4
- 238000000855 fermentation Methods 0.000 claims abstract description 26
- 230000004151 fermentation Effects 0.000 claims abstract description 26
- 230000033228 biological regulation Effects 0.000 claims abstract description 21
- 230000003534 oscillatory effect Effects 0.000 claims description 50
- 238000010790 dilution Methods 0.000 claims description 15
- 239000012895 dilution Substances 0.000 claims description 15
- 238000004088 simulation Methods 0.000 claims description 12
- 239000002028 Biomass Substances 0.000 claims description 8
- 239000000758 substrate Substances 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 abstract description 2
- 230000010355 oscillation Effects 0.000 description 7
- 238000005265 energy consumption Methods 0.000 description 6
- 240000004808 Saccharomyces cerevisiae Species 0.000 description 3
- 235000014680 Saccharomyces cerevisiae Nutrition 0.000 description 3
- 239000012528 membrane Substances 0.000 description 3
- 230000008016 vaporization Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000014759 maintenance of location Effects 0.000 description 2
- 230000001172 regenerating effect Effects 0.000 description 2
- 241000894006 Bacteria Species 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002551 biofuel Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000009931 harmful effect Effects 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 238000011031 large-scale manufacturing process Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000004540 process dynamic Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000005654 stationary process Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Preparation Of Compounds By Using Micro-Organisms (AREA)
- Micro-Organisms Or Cultivation Processes Thereof (AREA)
Abstract
运动发酵单胞菌生产乙醇中振荡过程的调控方法,采用实际生产数据建立运动发酵单胞菌连续发酵生产乙醇的动态模型,通过动态模型模拟计算获得运动发酵单胞菌乙醇连续发酵中的振荡过程与操作变量之间的定量关系,利用振荡过程与操作变量之间的关系来调控实际生产过程中的振荡现象来趋利避害,分析比较调控后的实际生产的数据与数值模拟是否相符,如果数值相符,继续采用该振荡过程与操作变量之间的关系调控;如果实际生产的数值与模拟数值不符,修改动态模型的参数,重复操作,直到最终调控后的实际生产数据与模拟预测的数据相符为止,本发明为提高发酵过程的稳定性,提高产物转化率,降低生产成本提供依据。
Description
技术领域
本发明属于生物乙醇的生产技术领域,特别涉及运动发酵单胞菌生产乙醇中振荡过程的调控方法。
背景技术
随着后化石经济时代的到来,用可再生的生物质资源替代化石能源已成必然。我国是石油资源相对贫乏的国家,据测算我国石油稳定供给不超过20年,很可能在我国实现“全面小康”的2020年,迎来石油供应丧失平衡的“拐点年”。为缓解能源压力,我国政府在不断推进可再生能源产业的相关研究,国家发改委于2007年制定的中长期可再生能源开发计划中制定目标,到2020年采用生物燃料替代1000万吨石油产品,占到全部交通燃料消费量的15%左右。生物乙醇需求量巨大,但是生产成本高制约了生物乙醇的大规模生产。目前生物乙醇的生产成本构成中,能耗成本约占30%,仅次于原料成本(约占60%),因此降低能耗是降低生物乙醇生产成本的主要方向。提高发酵产物中乙醇的含量,减少分离过程的能耗,是生物乙醇生产领域迫切需要解决的技术问题之一。
目前工业连续发酵生产生物乙醇中使用的菌种主要是酿酒酵母和运动发酵单胞菌。运动发酵单胞菌的乙醇产率高达97%,酿酒酵母只有90%,产物浓度高有利于降低能耗,除此之外,相比酿酒酵母,运动发酵单胞菌还具有生长和发酵能耗低、耐酒精度高、能在较高糖浓度中生长发酵、不需要控制氧浓度等优点,这些优点可以进一步降低能耗。因此,使用运动发酵单胞菌连续发酵生产生物乙醇过程具有更大的成本优势。
但是运动发酵单胞菌生产过程中存在复杂的振荡现象,在实验中观测到生物质和产物浓度持续大幅度的振荡行为。与稳态操作相比,振荡行为更难以预测和控制,从而影响生产过程的平稳,因此需要研究引发振荡现象的特定条件,开发相应的调控方法来减弱振荡过程带来的不良影响。为了充分发挥运动发酵单胞菌的优势,提高产物的转化率,必须解决稳定性的问题,因此,需要深入研究运动发酵单胞菌乙醇连续发酵过程的特性并提出振荡过程的定量调控方法。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供运动发酵单胞菌生产乙醇中振荡过程的调控方法,为提高发酵过程的稳定性,提高产物转化率,降低生产成本提供依据。
为了实现上述目的,本发明采用的技术方案是:
运动发酵单胞菌生产乙醇中振荡过程的调控方法,包括以下步骤:
第一步,采用实际生产数据建立运动发酵单胞菌连续发酵生产乙醇的动态模型其中x是该发酵过程模型的状态变量,α是该发酵过程模型的操作变量;所述的状态变量包括产物浓度、底物浓度和生物质浓度;所述的操作变量包括稀释速率和初始底物浓度;
第二步,通过动态模型模拟计算获得运动发酵单胞菌乙醇连续发酵中的振荡过程与操作变量之间的定量关系,具体为:
步骤2.1,通过改变操作变量,求取该发酵过程动态模型中对应的Hopf奇异点;
步骤2.2,通过模拟Hopf奇异点在这些操作点上的振荡过程,获取该过程的振荡过程数据,所述的振荡过程数据包括振幅和周期;
步骤2.3,通过分析振荡过程数据与操作变量的数值,建立振荡过程与操作变量之间的定量关系;
第三步,利用振荡过程与操作变量之间的关系来调控实际生产过程中的振荡现象来趋利避害;调控振荡现象包括调控振荡的振幅的大小和周期的长短;
第四步,比较调控后的实际生产的数据与数值模拟结果是否相符,如果数值相符,继续采用该振荡过程与操作变量之间的关系来调控振荡过程;如果实际生产的数值与模拟数值结果不符,修改动态模型的参数,重复第一步至第三步,直到最终调控后的实际生产数据与模拟预测的数据相符为止。
与现有技术相比,本发明的有益效果是:
对于发酵过程这样的慢过程,本发明相比传统的实验方法,能够克服实验过程由于周期较长、实验不充分而只有有限的实验数据;有些动态特性有可能被忽略掉的不足。使用本发明的方法可以得到完整的动态特性的数据,建立精确的调控模型,进而对连续发酵中的振荡过程进行有效调控。
附图说明
图1是本发明通过改变操作变量来调控运动发酵单胞菌连续生产乙醇中振荡过程的流程图。
图2是运动发酵单胞菌连续生产乙醇的流程示意图。
图3是进料稀释速率为0.1h-1时系统振荡过程。
图4是进料稀释速率为0.2h-1时系统振荡过程。
图5是预测振荡过程振幅为10kg/m3时系统的振荡过程。
图6是预测振荡过程振幅为20kg/m3时系统的振荡过程。
图7是进料稀释速率为0.052h-1时系统振荡逐渐衰减过程。
图8是进料稀释速率为0.050h-1时系统振荡迅速衰减过程。
具体实施方式
下面结合附图和实施例详细说明本发明。
如图1所示,运动发酵单胞菌生产乙醇中振荡过程的调控方法,包括以下步骤:
第一步,如图2所示,采用实际生产数据建立运动发酵单胞菌连续发酵生产乙醇的动态模型其中x是该发酵过程模型的状态变量,α是该发酵过程模型的操作变量,动态模型具体形式如下,其中状态变量Cx是生物质浓度,Cs是底物浓度,Ce是关键组分浓度,Cp是产物浓度,CpM是膜侧的产物浓度;操作变量为进料稀释速率Din、初始底物浓度Cs0;Dout是采出的稀释速率,DMout是膜侧采出的稀释速率;
发酵过程模型其中的数值如下:
表1动态模型中的参数值
其中,AM是渗透汽化膜的有效面积,Ce0是进料中的关键组分浓度,Cp0是进料中的产物浓度,CpM0是渗透汽化膜侧进料中的产物浓度,Cs0是初始底物浓度,Cx0是进料中的生物质浓度,DMin是膜侧稀释速率,Din是进料稀释速率,Ks是发酵反应的米氏常数,Vf是发酵有效体积,VM是渗透汽化膜侧的有效体积,P是渗透汽化速率,Ypx是产物的生物质收率因子,Ysx是底物的生物质收率因子,k1,k2,k3是反应常数,mp是产物保留常数,ms是底物保留常数,ρ是反应混合液的密度;
第二步,通过动态模型模拟计算获得运动发酵单胞菌乙醇连续发酵中的振荡过程与操作变量之间的定量关系,具体为:
步骤2.1,通过改变操作变量,例如,进料稀释速率、初始底物浓度等,求取该发酵过程动态模型中对应的Hopf奇异点;
发酵过程的Hopf奇异点数值如下表所示,在该点操作条件下,发酵过程的振荡现象随操作变量进料稀释速率的变化而改变,
表2发酵过程的Hopf奇异点的数值
步骤2.2,通过模拟Hopf奇异点在这些操作点上的振荡过程,获取该过程的振荡过程数据;所述的振荡过程数据包括振幅和周期;
改变稀释速率Din,模拟振荡现象,获取振荡过程的周期T和振幅Amp。Din=0.1h-1,Din=0.2h-1时,模拟计算的结果如图3、图4所示。
从计算结果可以得到Din=0.1h-1和Din=0.2h-1时,振荡过程中产物浓度的振幅与周期,同样的方法,可以得到Din=0.3h-1至Din=0.9h-1时对应的振幅和周期,如下表所示。
表3不同稀释速率下振荡过程的振幅与周期
步骤2.3,通过分析振荡过程数据与操作变量的数值,建立振荡过程与操作变量之间的定量关系,例如,振幅、周期与操作变量之间的关系。
Din=0.000037Amp 3-0.00012Amp 2-0.00028Amp+0.054
T=3.0252Din -0.702,或者Din=4.8421T-1.424
第三步,利用振荡过程与操作变量之间的关系来调控实际生产过程中的振荡现象来趋利避害,调控振荡现象包括调控振荡的振幅的大小和周期的长短;
比如,减小振荡的振幅,增大振荡的周期等;
预测当振幅Amp=10.00kg/m3时,Din=0.076h-1,T=18.42hr
模拟值Amp=9.30kg/m3,T=17.8hr,如图5所示。预测结果与模拟结果相符。
预测当振幅Amp=20.00kg/m3时,Din=0.297h-1,T=7.10hr
模拟值Amp=20.06kg/m3,T=7.1hr,如图6所示。预测结果与模拟结果相一致。
从预测公式可以得,当Din=0.54h-1时,Amp接近零,实际上在Din=0.50h-1到Din=0.54h-1之间,振荡过程的振幅数值接近于零,而且在有的数值下,振荡过程逐渐衰减,最终保持平稳过程。当Din=0.52h-1时,如图7所示。从中可以看到,经过一段时间的发展,振荡过程逐渐衰减并消失。当Din=0.50h-1时,振荡过程迅速衰减,如图8所示。
第四步,比较调控后的实际生产的数据与数值模拟结果是否相符,如果数值相符,继续采用该振荡过程与操作变量之间的关系来调控振荡过程;如果实际生产的数值与模拟数值结果不符,修改动态模型的参数,重复第一步至第三步,直到最终调控后的实际生产数据与模拟预测的数据相符为止。
Claims (1)
1.运动发酵单胞菌生产乙醇中振荡过程的调控方法,其特征在于,包括以下步骤:
第一步,采用实际生产数据建立运动发酵单胞菌连续发酵生产乙醇的动态模型其中x是该发酵过程模型的状态变量,α是该发酵过程模型的操作变量;所述的状态变量包括产物浓度、底物浓度和生物质浓度;所述的操作变量包括稀释速率和初始底物浓度;
第二步,通过动态模型模拟计算获得运动发酵单胞菌乙醇连续发酵中的振荡过程与操作变量之间的定量关系,具体为:
步骤2.1,通过改变操作变量,求取该发酵过程动态模型中对应的Hopf奇异点;
步骤2.2,通过模拟Hopf奇异点的振荡过程,获取该过程的振荡过程数据,所述的振荡过程数据包括振幅和周期;
步骤2.3,通过分析振荡过程数据与操作变量的数值,建立振荡过程与单一操作变量之间的定量关系;
第三步,利用振荡过程与操作变量之间的定量关系来调控实际生产过程中的振荡现象来趋利避害;调控振荡现象包括调控振荡的振幅的大小和周期的长短;
第四步,比较调控后的实际生产的数据与数值模拟结果是否相符,如果数值相符,继续采用该振荡过程与操作变量之间的关系来调控振荡过程;如果实际生产的数值与模拟数值结果不符,修改动态模型的参数,重复第一步至第三步,直到最终调控后的实际生产数据与模拟预测的数据相符为止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410059227.3A CN103810391B (zh) | 2014-02-21 | 2014-02-21 | 运动发酵单胞菌生产乙醇中振荡过程的调控方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410059227.3A CN103810391B (zh) | 2014-02-21 | 2014-02-21 | 运动发酵单胞菌生产乙醇中振荡过程的调控方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103810391A CN103810391A (zh) | 2014-05-21 |
CN103810391B true CN103810391B (zh) | 2017-05-17 |
Family
ID=50707151
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410059227.3A Active CN103810391B (zh) | 2014-02-21 | 2014-02-21 | 运动发酵单胞菌生产乙醇中振荡过程的调控方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103810391B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101565706A (zh) * | 2008-04-25 | 2009-10-28 | 武汉市星站生物技术有限公司 | 构建基因重组运动发酵单胞菌应用于乙醇发酵 |
CN102337304A (zh) * | 2010-07-27 | 2012-02-01 | 中国科学院成都生物研究所 | 运动发酵单胞菌快速生产乙醇的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8476048B2 (en) * | 2010-12-17 | 2013-07-02 | E I Du Pont De Nemours And Company | Xylose utilizing zymomonas mobilis with improved ethanol production in biomass hydrolysate medium |
-
2014
- 2014-02-21 CN CN201410059227.3A patent/CN103810391B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101565706A (zh) * | 2008-04-25 | 2009-10-28 | 武汉市星站生物技术有限公司 | 构建基因重组运动发酵单胞菌应用于乙醇发酵 |
CN102337304A (zh) * | 2010-07-27 | 2012-02-01 | 中国科学院成都生物研究所 | 运动发酵单胞菌快速生产乙醇的方法 |
Non-Patent Citations (2)
Title |
---|
Analysis of Hopf Points for a Zymomonas mobilis Continuous Fermentation Process Producing Ethanol;Hangzhou wang等;《Industrial & Engineering Chemistry Research》;20121218;第1645-1654页 * |
化工过程稳定稳态点抗扰动能力定量分析方法;王杭州等;《化工学报》;20120131;第63卷(第1期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103810391A (zh) | 2014-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111914414B (zh) | 质子交换膜电解池三维全电池模型的建立方法 | |
Petre et al. | Advanced nonlinear control strategies for a fermentation bioreactor used for ethanol production | |
CN105091590A (zh) | 烧结主抽风机智能节电的方法及系统 | |
CN103810391B (zh) | 运动发酵单胞菌生产乙醇中振荡过程的调控方法 | |
Singh et al. | Frequency regulation of isolated hybrid wind/diesel, power generation with fuel cell system | |
Liu et al. | Kinetic modeling for redox potential-controlled repeated batch ethanol fermentation using flocculating yeast | |
CN105787241B (zh) | 一种析氢电极材料性能的优化方法 | |
Converti et al. | Simplified modeling of fed‐batch alcoholic fermentation of sugarcane blackstrap molasses | |
CN108388132A (zh) | 一种聚合酶链反应过程系统控制模型建模方法 | |
Dong et al. | Development and optimization of an unstructured kinetic model for sodium gluconate fermentation process | |
Lin et al. | Computer simulation of bioprocess | |
CN111783315B (zh) | 基于数据驱动的燃料乙醇发酵罐早停止方法 | |
CN103617316A (zh) | 一种基于多核支持向量机的发酵过程动态预测方法 | |
Ciesielski et al. | The dumping of oscillatory phenomena in the process of bioethanol production by continuous fermentation | |
CN113110628A (zh) | 基于pso的压水堆除氧器水位控制方法 | |
Ljunggren | Kinetic analysis and modeling of enzymatic hydrolysis and SSF | |
CN107194129B (zh) | 一种应用乳酸修正的adm1模型分析餐厨垃圾产甲烷的方法 | |
Mahar et al. | Parametric effect on ethanol production though use of Monod model | |
WO2019114402A1 (zh) | 一种秸秆直燃循环流化床锅炉变负荷蒸汽压跟踪控制方法 | |
CN101718495A (zh) | 烧结机点火炉的复合控制装置 | |
Carvajal et al. | Simulation of the batch fermentation stage in the process to obtain ethanol from final molasse | |
Muñoz-Tamayo et al. | Optimization of lipid production by oleaginous yeast in continuous culture | |
TW201530277A (zh) | 高爐鐵水溫度預測方法 | |
Sridhar | Global optimization of continuous fermentation involving Zymomonas mobilis | |
Gomez et al. | Mathematical modeling of a raceway pond system for biofuels production |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |