CN116640886A - Fermentation control method, device, equipment and readable storage medium - Google Patents
Fermentation control method, device, equipment and readable storage medium Download PDFInfo
- Publication number
- CN116640886A CN116640886A CN202310658351.0A CN202310658351A CN116640886A CN 116640886 A CN116640886 A CN 116640886A CN 202310658351 A CN202310658351 A CN 202310658351A CN 116640886 A CN116640886 A CN 116640886A
- Authority
- CN
- China
- Prior art keywords
- fermentation
- parameter
- feeding
- theoretical
- concentration
- 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.)
- Pending
Links
- 238000000855 fermentation Methods 0.000 title claims abstract description 251
- 230000004151 fermentation Effects 0.000 title claims abstract description 251
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000003860 storage Methods 0.000 title claims abstract description 13
- 230000008859 change Effects 0.000 claims abstract description 76
- 239000000758 substrate Substances 0.000 claims abstract description 58
- 239000003513 alkali Substances 0.000 claims abstract description 50
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 48
- 239000007788 liquid Substances 0.000 claims abstract description 31
- 238000012840 feeding operation Methods 0.000 claims abstract description 29
- YPFDHNVEDLHUCE-UHFFFAOYSA-N 1,3-propanediol Substances OCCCO YPFDHNVEDLHUCE-UHFFFAOYSA-N 0.000 claims description 41
- 239000000047 product Substances 0.000 claims description 31
- DNIAPMSPPWPWGF-VKHMYHEASA-N (+)-propylene glycol Chemical compound C[C@H](O)CO DNIAPMSPPWPWGF-VKHMYHEASA-N 0.000 claims description 27
- 229940035437 1,3-propanediol Drugs 0.000 claims description 27
- 229920000166 polytrimethylene carbonate Polymers 0.000 claims description 27
- 230000008569 process Effects 0.000 claims description 19
- 230000006870 function Effects 0.000 claims description 17
- 238000012545 processing Methods 0.000 claims description 13
- 230000000875 corresponding effect Effects 0.000 claims description 11
- 230000001276 controlling effect Effects 0.000 claims description 10
- 238000004088 simulation Methods 0.000 claims description 10
- 239000002028 Biomass Substances 0.000 claims description 9
- 230000002378 acidificating effect Effects 0.000 claims description 9
- 239000006052 feed supplement Substances 0.000 claims description 5
- 239000000126 substance Substances 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000002596 correlated effect Effects 0.000 claims description 3
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 239000002585 base Substances 0.000 claims description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 18
- 238000010586 diagram Methods 0.000 description 12
- 238000004519 manufacturing process Methods 0.000 description 6
- 238000012986 modification Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 238000002360 preparation method Methods 0.000 description 5
- 230000001133 acceleration Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 239000012467 final product Substances 0.000 description 4
- 244000005700 microbiome Species 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 3
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 2
- 230000004075 alteration Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 241000894006 Bacteria Species 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000006227 byproduct Substances 0.000 description 1
- 230000005574 cross-species transmission Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004066 metabolic change Effects 0.000 description 1
- 230000002503 metabolic effect Effects 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 230000000813 microbial effect Effects 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 241001148471 unidentified anaerobic bacterium Species 0.000 description 1
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q3/00—Condition responsive control processes
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12M—APPARATUS FOR ENZYMOLOGY OR MICROBIOLOGY; APPARATUS FOR CULTURING MICROORGANISMS FOR PRODUCING BIOMASS, FOR GROWING CELLS OR FOR OBTAINING FERMENTATION OR METABOLIC PRODUCTS, i.e. BIOREACTORS OR FERMENTERS
- C12M41/00—Means for regulation, monitoring, measurement or control, e.g. flow regulation
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12P—FERMENTATION OR ENZYME-USING PROCESSES TO SYNTHESISE A DESIRED CHEMICAL COMPOUND OR COMPOSITION OR TO SEPARATE OPTICAL ISOMERS FROM A RACEMIC MIXTURE
- C12P7/00—Preparation of oxygen-containing organic compounds
- C12P7/02—Preparation of oxygen-containing organic compounds containing a hydroxy group
- C12P7/04—Preparation of oxygen-containing organic compounds containing a hydroxy group acyclic
- C12P7/18—Preparation of oxygen-containing organic compounds containing a hydroxy group acyclic polyhydric
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Abstract
The embodiment of the invention provides a fermentation control method, a device, equipment and a readable storage medium, wherein the basic feeding rate of a substrate of a fermentation broth is determined according to an alkali adding change parameter when the alkali adding change parameter meets a first preset condition by continuously acquiring the alkali adding change parameter of the fermentation broth, and feeding operation is performed based on the basic feeding rate. When the fermentation time of the fermentation liquid is within a preset period, a first fermentation state parameter of the fermentation liquid in the current preset period is obtained, and the first fermentation state parameter is processed by utilizing a preset predictive control algorithm so as to obtain a theoretical feeding change parameter of a fermentation substrate in the current preset period. And finally, controlling the feeding operation of the fermentation liquid in the next preset period by utilizing the theoretical feeding change parameter. On the basis of the basic feeding rate, the feeding operation of the fermentation broth in the next preset period is controlled by using a predictive control algorithm, so that the accuracy of feeding rate control can be further improved.
Description
Technical Field
The present invention relates to the field of bioengineering technology, and more particularly, to a fermentation control method, apparatus, device, and readable storage medium.
Background
The 1, 3-propylene glycol fermentation system is a batch fermentation process which takes anaerobic bacteria as a fermentation strain and crude glycerol as a substrate, wherein the fermentation bacteria produce acidic substances in the fermentation process, the pH value is maintained at 7 by corresponding control equipment, and the temperature is maintained at about 37 ℃ by the corresponding control equipment. Wherein, but in fermentation process on-line measuring value includes: current feed rate, accumulated feed amount, current alkali rate, accumulated alkali amount, pH and temperature.
If the optimal control algorithm is adopted to control the fermentation process in real time, information such as biomass concentration, substrate concentration, product concentration and the like is needed, however, the existing fermentation control system can not provide corresponding biological fermentation process information, and only the online detection value can be provided. The traditional off-line track fermentation measurement lacks real-time feedback control, a large amount of experiments are needed to find a fermentation track with higher product concentration, and if the amount of fed-batch is too small, the process cannot reach the maximum product concentration. Overfeeding can have negative effects on fermentation, leading to metabolic spillover and thus the formation of byproducts.
In order to increase the final product concentration, a suitable feed rate must be found to saturate the rate of product formation. The optimal feed rate is constantly changing with strong nonlinearity due to the exponential growth rate, metabolic changes, volume dynamics and the remaining disturbances of the microbial growth. The nonlinear dynamics of the microorganism and the disturbance existing in the actual process make it desirable to obtain a higher product concentration, and many links for improvement are needed for the control of the fermentation process.
Disclosure of Invention
The embodiment of the invention solves the technical problem of low preparation efficiency in the process of preparing 1, 3-propanediol by fermentation in the prior art by providing a fermentation control method, a device, equipment and a readable storage medium.
In a first aspect, the present invention provides a fermentation control method applied to a process for preparing 1, 3-propanediol from a fermentation broth; the fermentation control method comprises the following steps: continuously acquiring an alkali adding change parameter for the fermentation broth, determining a basic feeding rate for a substrate of the fermentation broth according to the alkali adding change parameter when the alkali adding change parameter meets a first preset condition, and carrying out feeding operation based on the basic feeding rate; when the fermentation time of the fermentation liquid is within a preset period, acquiring a first fermentation state parameter of the fermentation liquid within the current preset period, and processing the first fermentation state parameter by utilizing a preset predictive control algorithm to obtain a theoretical feeding change parameter of the fermentation substrate within the current preset period; and controlling the feeding operation of the fermentation liquor in the next preset period by utilizing the theoretical feeding variation parameters.
As an optional implementation manner, the alkali adding variation parameter meets a first preset condition, and includes: and the alkali adding change parameter is monotonically decreased along with the time change within the preset time.
As an alternative embodiment, the determining the basic flow rate for the fermentation broth substrate according to the alkali addition variation parameter includes: determining the basic flow rate based on the product of the alkali addition change parameter and a first preset coefficient; wherein the magnitude of the first preset coefficient is positively correlated with the volume of the fermentation broth.
As an alternative embodiment, the first fermentation status parameter includes: the concentration of biomass in the fermentation broth, the concentration of substrate in the fermentation broth, the concentration of 1, 3-propanediol in the fermentation broth, the concentration of acidic species in the fermentation broth, and the volume of the fermentation broth;
the processing the first fermentation state parameter by using a preset predictive control algorithm to obtain a theoretical feeding change parameter of the fermentation substrate in a current preset period comprises the following steps: taking the concentration of 1, 3-propylene glycol in the fermentation liquor as an objective function, and taking the biomass concentration in the fermentation liquor, the concentration of a substrate in the fermentation liquor, the concentration of an acidic substance in the fermentation liquor, the volume of the fermentation liquor and the theoretical feeding rate in the current preset period as constraint conditions of the objective function; and determining a corresponding theoretical feeding rate in the process of calculating the extremum of the objective function so as to obtain the theoretical feeding change parameter based on a plurality of theoretical feeding rates.
As an alternative embodiment, after the controlling the feeding operation of the fermentation broth in the next preset period by using the theoretical feeding variation parameter, the method further includes: acquiring a second fermentation state parameter of the fermentation liquid in the next preset period, and processing the second fermentation state parameter and the theoretical feed supplement change parameter by utilizing the predictive control algorithm to obtain a theoretical concentration change parameter of 1, 3-propanediol in the fermentation liquid in the next preset period; and correcting the predictive control algorithm by using the actual concentration variation parameter of the 1, 3-propanediol in the first fermentation state parameter and the theoretical concentration variation parameter.
As an alternative embodiment, the correcting the predictive control algorithm by using the actual concentration variation parameter of the 1, 3-propanediol in the first fermentation status parameter and the theoretical concentration variation parameter includes: correspondingly modifying a second preset coefficient in the predictive control algorithm based on the difference value between the actual concentration variation parameter and the theoretical concentration variation parameter; and if the difference value between the actual concentration variation parameter and the theoretical concentration variation parameter is larger, the second preset coefficient is smaller.
As an alternative embodiment, after said modifying the predictive control algorithm, the method further comprises: circularly executing the step of correspondingly correcting a second preset coefficient in the predictive control algorithm based on the difference value between the actual concentration variation parameter and the theoretical concentration variation parameter; until the number of cycles is greater than a first preset threshold, or until the difference between the actual concentration variation parameter and the theoretical concentration variation parameter is less than a second preset threshold.
In a second aspect, the present invention provides a fermentation control apparatus for use in a process for preparing 1, 3-propanediol from a fermentation broth; the fermentation control device comprises:
the data calculation unit is used for continuously acquiring an alkali adding change parameter aiming at the fermentation broth, determining a basic feeding rate aiming at a substrate of the fermentation broth according to the alkali adding change parameter when the alkali adding change parameter meets a first preset condition, and carrying out feeding operation based on the basic feeding rate;
the simulation prediction unit is used for acquiring a first fermentation state parameter of the fermentation broth in a current preset period when the fermentation time of the fermentation broth is in the preset period, and processing the first fermentation state parameter by utilizing a preset prediction control algorithm to acquire a theoretical feeding change parameter of the fermentation substrate in the current preset period;
and the feeding control unit is used for controlling the feeding operation of the fermentation liquid in the next preset period by utilizing the theoretical feeding change parameter.
In a third aspect, the present invention provides, by an embodiment of the present invention, a fermentation control apparatus comprising a memory, a processor and code stored on the memory and executable on the processor, the processor implementing any of the embodiments of the first aspect when executing the code.
In a fourth aspect, the present invention provides, by way of example of the present invention, a readable storage medium having stored thereon a program which when executed by a processor implements any of the embodiments of the first aspect.
One or more technical solutions provided in the embodiments of the present invention at least have the following technical effects or advantages:
firstly, continuously obtaining an alkali adding change parameter for fermentation liquor, determining a basic feeding rate for a substrate of the fermentation liquor according to the alkali adding change parameter when the alkali adding change parameter meets a first preset condition, and carrying out feeding operation based on the basic feeding rate.
Secondly, when the fermentation time of the fermentation liquid is within a preset period, a first fermentation state parameter of the fermentation liquid in the current preset period is obtained, and the first fermentation state parameter is processed by utilizing a preset predictive control algorithm so as to obtain a theoretical feeding change parameter of a fermentation substrate in the current preset period.
And finally, controlling the feeding operation of the fermentation liquid in the next preset period by utilizing the theoretical feeding change parameter. The basic feeding rate can ensure the basic accuracy of the feeding rate, and based on the basic feeding rate, the feeding operation of the fermentation broth in the next preset period is controlled by utilizing the theoretical feeding change parameter obtained by the predictive control algorithm, so that the accuracy of feeding rate control can be further improved, and the technical problem of low preparation efficiency in the process of preparing 1, 3-propanediol by fermentation in the prior art is effectively solved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for the description of the embodiments will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a fermentation control method in an embodiment of the invention;
FIG. 2 is a graph showing the theoretical concentration variation parameters of the target product without feed control in the embodiment of the invention;
FIG. 3 is a schematic diagram of the feed rate control based on the basic stream acceleration rate in an embodiment of the present invention;
FIG. 4 shows the theoretical concentration variation parameters of the target product under the feeding control based on the basic stream acceleration rate in the embodiment of the invention;
FIG. 5 is a schematic diagram of the feed rate control based on the basic stream acceleration rate and the predictive control algorithm in an embodiment of the invention;
FIG. 6 shows the theoretical concentration variation parameters of the target product under the feed supplement control based on the basic stream acceleration rate and the predictive control algorithm in the embodiment of the invention;
FIG. 7 is a diagram showing the variation of the difference between the theoretical concentration variation parameter and the actual concentration variation parameter according to the embodiment of the present invention;
FIG. 8 is a schematic diagram of a fermentation control apparatus according to an embodiment of the present invention;
FIG. 9 is a schematic diagram showing the structure of a fermentation control apparatus according to an embodiment of the present invention;
FIG. 10 is a schematic diagram of a readable storage medium structure according to an embodiment of the present invention.
Detailed Description
The embodiment of the invention solves the technical problem of low preparation efficiency in the process of preparing 1, 3-propanediol by fermentation in the prior art by providing a fermentation control method, a device, equipment and a readable storage medium.
The technical scheme provided by the embodiment of the invention aims to solve the technical problems, and the overall thought is as follows:
firstly, continuously obtaining an alkali adding change parameter for fermentation liquor, determining a basic feeding rate for a substrate of the fermentation liquor according to the alkali adding change parameter when the alkali adding change parameter meets a first preset condition, and carrying out feeding operation based on the basic feeding rate.
Secondly, when the fermentation time of the fermentation liquid is within a preset period, a first fermentation state parameter of the fermentation liquid in the current preset period is obtained, and the first fermentation state parameter is processed by utilizing a preset predictive control algorithm so as to obtain a theoretical feeding change parameter of a fermentation substrate in the current preset period.
And finally, controlling the feeding operation of the fermentation liquid in the next preset period by utilizing the theoretical feeding change parameter.
In order to better understand the above technical solutions, the following detailed description will refer to the accompanying drawings and specific embodiments.
First, the term "and/or" appearing herein is merely an association relationship describing associated objects, meaning that there may be three relationships, e.g., a and/or B, may represent: a exists alone, A and B exist together, and B exists alone. In addition, the character "/" herein generally indicates that the front and rear associated objects are an "or" relationship.
It should be noted that the terms "first," "second," and the like in the description and the claims of the present invention and the above figures are used for distinguishing between similar objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used may be interchanged where appropriate such that the embodiments of the invention described herein may be capable of operation in sequences other than those illustrated or otherwise described.
In a first aspect, the present invention provides a fermentation control method, which can be applied to a process for preparing 1, 3-propanediol from a fermentation broth, and can be applied to a process for preparing other fermentation products, such as ethanol.
Referring to fig. 1, the fermentation control method may include the following steps S101 to S103:
step S101: continuously obtaining alkali adding change parameters for the fermentation liquor, determining a basic feeding rate for a substrate of the fermentation liquor according to the alkali adding change parameters when the alkali adding change parameters meet a first preset condition, and carrying out feeding operation based on the basic feeding rate.
Specifically, the alkali addition change parameter meeting the first preset condition includes: the alkali adding change parameter is monotonically decreased along with the change of time within the preset duration.
In the specific implementation process, the alkali adding change parameter in the preset time period characterizes the change condition of the alkali adding rate in the preset time period, and the preset time period can be set according to actual application requirements. As long as the maximum value exists in the alkali adding rate within the preset duration, and after the corresponding time of the maximum value, the alkali adding rate is reduced along with the time, the shortage of the substrate of the fermentation broth is represented, and the feeding operation can be started from the corresponding time of the maximum value of the alkali adding rate.
For how to determine the basic flow rate for the substrate of the fermentation broth, in particular, the basic flow rate may be determined based on the product of the alkali addition variation parameter and the first preset coefficient. Wherein the magnitude of the first preset coefficient is positively correlated with the volume of the fermentation broth.
In one embodiment, the base stream rate may be calculated using the following formula:
F 1 ＝e+fF 2
wherein F is 1 For basic flow rate of substrate of fermentation broth, F 2 The alkali liquor adding rate (alkali adding rate) is represented, f is a first preset coefficient, and e is a third preset coefficient.
In an alternative embodiment, the first preset factor may be 2.8 and the third preset factor may be 0.01.
It should be noted that, the first preset coefficient and the third preset coefficient may be determined by historical experimental data, where the historical experimental data includes: the relation between the substrate concentration of the fermentation broth and time when no feed is fed, the relation between the substrate concentration of the fermentation broth and time when the feed is fed, the relation between the substrate concentration of the fermentation broth and the alkali adding rate when the feed is not fed, and the like.
Second, the first and third preset coefficients are related to the volume of the fermentation broth, and if the volume of the fermentation broth is large, the values of the first and third preset coefficients also need to be increased accordingly.
Step S102: when the fermentation time of the fermentation liquid is within a preset period, a first fermentation state parameter of the fermentation liquid in the current preset period is obtained, and the first fermentation state parameter is processed by utilizing a preset predictive control algorithm so as to obtain a theoretical feeding change parameter of a fermentation substrate in the current preset period.
The preset period may be a certain period of time after the feeding operation based on the basic stream rate, and may be set according to an actual application scenario or application experience.
Specifically, the first fermentation status parameter may include: the biomass concentration in the fermentation broth, the substrate concentration in the fermentation broth, the 1, 3-propanediol concentration in the fermentation broth, the acidic material concentration in the fermentation broth, and the volume of the fermentation broth.
For how to obtain the theoretical feeding change parameter of the fermentation substrate in the current preset period, specifically, the concentration of 1, 3-propanediol in the fermentation liquid can be taken as an objective function, and the biomass concentration in the fermentation liquid, the concentration of the substrate in the fermentation liquid, the concentration of acidic substances in the fermentation liquid, the volume of the fermentation liquid and the theoretical feeding rate can be taken as constraint conditions of the objective function in the current preset period. And determining the corresponding theoretical feeding rate in the process of calculating the extremum of the objective function, so as to obtain theoretical feeding change parameters based on a plurality of theoretical feeding rates.
In the specific implementation process, the substrate feeding rate is used as a main control quantity, the pH value of the fermentation liquor is used as an auxiliary control quantity, and the actual feeding rate is controlled based on the optimal substrate feeding rate by calculating the optimal substrate feeding rate track so as to maximize the yield of the 1, 3-propanediol at the end of fermentation.
In this regard, in an alternative embodiment, the problem of optimizing the substrate feed rate may be described as follows:
the objective function J can be expressed as:
the constraints of the objective function described above may be expressed as follows:
P b (t+H)＝f(t,X(t),S(t),P a (t),Q(t),V(t))；
Q min ≤Q(t)≤Q max ；
V 0 ≤V(t)≤V max ；
X(t)＝X i ；
S(t)＝S i ；
P a (t)＝P ai ；
in the constraint condition, t is a preset period, t f For the preset period end time, P b The concentration of the 1, 3-propylene glycol is H, X (t) is the biomass concentration in the fermentation broth at the moment t, S (t) is the substrate concentration in the fermentation broth at the moment t and P are the prediction period of a prediction control algorithm a (t) is the concentration of acidic substances in the fermentation broth at time t, V (t) is the volume of the fermentation broth at time t, and Q (t) is the feed rate in the fermentation broth at time t.
Wherein Q is min At the feed rate ofLimit value, Q max Is the upper limit value of the feeding rate, V 0 Is the lower limit of the volume of fermentation liquor, V max An upper limit of the volume of the fermentation broth.
Substituting the first fermentation state parameter into the constraint condition, and determining a corresponding theoretical feeding rate in the process of calculating the objective function, so as to obtain a theoretical feeding change parameter based on a plurality of theoretical feeding rates.
Step S103: and controlling the feeding operation of the fermentation broth in the next preset period by utilizing the theoretical feeding variation parameters.
Because the duration of the current preset period is consistent with that of the next preset period, after the theoretical feeding change parameter is obtained, the feeding operation of the fermentation liquid in the next preset period can be controlled by utilizing the theoretical feeding change parameter.
As an alternative embodiment, after the feeding operation of the fermentation broth in the next preset period is controlled by using the theoretical feeding variation parameter, a second fermentation state parameter of the fermentation broth in the next preset period may be further obtained, and the second fermentation state parameter and the theoretical feeding variation parameter are processed by using a predictive control algorithm, so as to obtain a theoretical concentration variation parameter of 1, 3-propanediol in the fermentation broth in the next preset period.
And then, correcting the predictive control algorithm by using the actual concentration variation parameter and the theoretical concentration variation parameter of the 1, 3-propanediol in the first fermentation state parameters.
Specifically, the second preset coefficient in the predictive control algorithm may be modified correspondingly based on the difference between the actual concentration variation parameter and the theoretical concentration variation parameter. And if the difference value between the actual concentration variation parameter and the theoretical concentration variation parameter is larger, the second preset coefficient is smaller.
In one embodiment, the second preset coefficients may be modified accordingly using the following formula:
u＝0.01+2.8*dcdt+0.003*e
wherein u is a second preset coefficient, dcdt is an alkali addition rate, and e is a third preset coefficient.
As an alternative embodiment, after the correction of the predictive control algorithm, the step of performing corresponding correction on the second preset coefficient in the predictive control algorithm based on the difference between the actual concentration variation parameter and the theoretical concentration variation parameter may also be performed in a cyclic manner. Until the number of cycles is greater than a first preset threshold, or until the difference between the actual concentration variation parameter and the theoretical concentration variation parameter is less than a second preset threshold.
In the specific implementation process, the theoretical concentration change parameter characterizes a theoretical concentration value of the target product changing along with time in a preset period. Correspondingly, the actual concentration variation parameter represents the actual concentration value of the target product varying along with time in a preset period, and can be obtained by detecting the actual concentration value of the target product on line.
The first preset threshold may be set according to an actual application scenario, and in an alternative embodiment, the first preset threshold may be 100. The second preset threshold may be set according to an actual application scenario, and in an alternative embodiment, the second preset threshold may be 10 g/l.
In order to better show the application effect of the fermentation control method, simulation experiments are carried out according to the following conditions: the fermentation substrate is glycerol, the volume of the fermentation tank is 5 liters, the concentration of the matrix glycerol in the fermentation liquid is 500 g/liter, the lower limit value of the volume of the fermentation liquid is 3 liters, the upper limit value of the volume of the fermentation liquid is 4 liters according to actual fermentation site data, the prediction period of a prediction control algorithm is 1 hour, the fermentation termination time is 30 hours, the lower limit value of the feeding rate of the matrix glycerol is 0, and the upper limit value of the feeding rate of the matrix glycerol is 0.2 g/liter.
The fermentation simulation experiment can be based on MATLAB software, and the optimization problem in the predictive control algorithm is added with a punishment item of the accumulated amount of the fed materials on the basis of the optimization problem of the previous purpose, so that the highest concentration of the target product can be obtained by using the minimum glycerol consumption. The objective function is shown as follows:
wherein y (k) represents the biomass concentration in the fermentation broth, the substrate concentration in the fermentation broth, the concentration of 1, 3-propanediol, the concentration of acidic species in the fermentation broth, and the volume of the fermentation broth at time k. u (k) represents the feed rate at time k, and Q and R are both weight matrices.
Converting the optimization problem into an NLP (non-linear programming) problem through a single-shooter method, solving an optimal action sequence at the moment by utilizing an optimization open source algorithm library CasADI, taking a first action as a current action value, putting the first action into fermentation simulation at the moment, and circularly repeating to obtain a series of state tracks of the concentration of the target product (1, 3-propanediol) in a preset period, wherein the series of state tracks are theoretical concentration change parameters.
And then in different fermentation stages, correspondingly modifying a second preset coefficient in the predictive control algorithm based on the difference value of the actual concentration variation parameter and the theoretical concentration variation parameter. For example, if the actual concentration value obtained at a certain sampling time is greater than the theoretical concentration value, the first preset coefficient is appropriately adjusted before the next sampling time, so that the subsequent basic flow rate is reduced to approach the optimal theoretical concentration value of the target product at the next sampling time.
In the first simulation experiment, if the fermentation broth does not have corresponding feed control, theoretical concentration change parameters of the target product in a preset period are shown in fig. 2, and the final product concentration of the target product (1, 3-propanediol) is about 29g/L when no feed control is performed can be seen from fig. 2.
In the simulation experiment II, if the substrate feeding operation is performed only based on the basic feeding rate shown in fig. 3, the theoretical concentration change parameter of the target product in the preset period is shown in fig. 4, and the final product concentration of the target product (1, 3-propylene glycol) is about 56g/L without feeding control can be seen from fig. 4. Wherein the starting time of the feeding is determined by the derivative of the alkali adding rate, and the feeding is started when the alkali adding rate reaches the peak to start to decrease. The reason for setting the feeding time is that when the alkali adding rate reaches the highest value, the microorganism propagation rate and the metabolism rate are always in the peak period, the substrate consumption rate of the microorganism is the fastest, and the substrate can be fed to the microorganism just at the moment. Second, empirically, the substrate concentration at this time is often not too high to produce excessive substrate inhibition.
Simulation experiment three, if the feed control shown in fig. 5 is adopted: on the basis of feeding operation based on the basic flow rate, a preset predictive control algorithm is utilized to process the first fermentation state parameter, and the theoretical feeding change parameter of the fermentation substrate in the current preset period is obtained. And finally, controlling the feeding operation of the fermentation broth in the next preset period by utilizing the theoretical feeding variation parameters.
Correspondingly, the theoretical concentration change parameters of the target product in the preset period are shown in fig. 6, and the final product concentration of the target product (1, 3-propylene glycol) is about 70g/L in the no-feed control process can be seen from fig. 6.
With continued Reference to fig. 6, the MPC Reference curve was characterized: the current moment prediction control algorithm utilizes the theoretical concentration change track (theoretical concentration change parameter) of the substrate in the next preset period, which is obtained by utilizing the theoretical feed supplement change parameter of the previous preset period. If the preset period is 1 hour, the running time of the predictive control algorithm starts from 20 hours, the fermentation substrate is sampled for the first time at 20 hours, and the actual fermentation state parameter of the hour is obtained by sampling at 21 hours.
If the predictive control algorithm is based on the basic flow rate and the actual fermentation state parameters obtained at the moment of 20 hours, the theoretical concentration change parameters of the substrate at 20-21 hours are calculated, and meanwhile, the actual fermentation state parameters at the hour are sampled once again and obtained at 22 hours. And comparing the actual fermentation state parameters (including the actual change parameters of the substrate concentration) obtained at the moment 21h with the theoretical concentration change parameters of the 21h, and if the substrate concentration value at the moment 21h is higher than the actual substrate concentration value at the moment 21h, increasing the substrate feeding rate in the next hour. And the like, and stopping feeding until the fermentation time reaches 30 hours.
As shown in fig. 7, the error between the theoretical concentration variation parameter and the MPC Reference is smaller and smaller under the adjustment control of the predictive control algorithm, which characterizes that the actual concentration variation of the substrate approaches to the direction optimized by the predictive control algorithm, and the theoretical concentration variation parameter obtained by using the predictive control algorithm is closer to the actual application scene. Therefore, the concentration of the target product (1, 3-propylene glycol) can be improved through the theoretical feed supplement change parameter obtained by the predictive control algorithm.
In summary, the three simulation experiment results show that: the final fermentation product concentration of the fermentation broth is only 29g/L without feed control; the final fermentation product concentration of the fermentation broth under feed control with a basal feed rate was 56g/L. The final fermentation product concentration of the fermentation liquid can reach 70g/L under the basic flow rate and the feed control of a predictive control algorithm, and is greatly improved compared with the former two.
In a second aspect, the present invention provides a fermentation control apparatus, which can be applied to a process for preparing 1, 3-propanediol from a fermentation broth, based on the same inventive concept. Referring to fig. 8, the fermentation control apparatus may include:
a data calculation unit 201, configured to continuously obtain an alkali addition variation parameter for the fermentation broth, determine a basic feeding rate for a substrate of the fermentation broth according to the alkali addition variation parameter when the alkali addition variation parameter meets a first preset condition, and perform a feeding operation based on the basic feeding rate;
the simulation prediction unit 202 is configured to obtain a first fermentation state parameter of the fermentation broth in a current preset period when the fermentation time of the fermentation broth is in the preset period, and process the first fermentation state parameter by using a preset prediction control algorithm to obtain a theoretical feeding variation parameter of the fermentation substrate in the current preset period;
and the feeding control unit 203 is configured to control feeding operation of the fermentation broth in a next preset period by using the theoretical feeding variation parameter.
Since the fermentation control apparatus described in this embodiment is an electronic device for implementing the fermentation control method in this embodiment, based on the fermentation control method described in this embodiment, those skilled in the art can understand the specific implementation of the electronic device and various modifications thereof, so how the electronic device implements the method in this embodiment will not be described in detail herein. Any electronic device used by those skilled in the art to implement the fermentation control method according to the embodiments of the present invention falls within the scope of the present invention.
In a third aspect, embodiments of the present invention provide a fermentation control apparatus, which may be applied to the preparation of 1, 3-propanediol, based on the same inventive concept.
Referring to fig. 9, a fermentation control apparatus according to an embodiment of the present invention includes: memory 901, processor 902, and code stored on the memory and executable on processor 902, processor 902 implementing any one of the embodiments of the fermentation control methods described above when the code is executed.
Where in FIG. 9, a bus architecture (represented by bus 900), the bus 900 may include any number of interconnected buses and bridges, with the bus 900 linking together various circuits, including one or more processors, represented by processor 902, and memory, represented by memory 901. Bus 900 may also link together various other circuits such as peripheral devices, voltage regulators, power management circuits, etc., as are well known in the art and, therefore, will not be described further herein. The bus interface 905 provides an interface between the bus 900 and the receiver 903 and the transmitter 904. The receiver 903 and the transmitter 904 may be the same element, i.e. a transceiver, providing a unit for communicating with various other apparatus over a transmission medium. The processor 902 is responsible for managing the bus 900 and general processing, while the memory 901 may be used to store data used by the processor 902 in performing operations.
In a fourth aspect, as shown in fig. 10, the present invention provides, by way of an embodiment of the present invention, a readable storage medium 1000 having stored thereon a program 1001, which program 1001, when executed by a processor, implements any of the embodiments of the foregoing fermentation control method.
The technical scheme provided by the embodiment of the invention at least has the following technical effects or advantages:
the basic feeding rate can ensure the basic accuracy of the feeding rate, and based on the basic feeding rate, the feeding operation of the fermentation broth in the next preset period is controlled by utilizing the theoretical feeding change parameter obtained by the predictive control algorithm, so that the accuracy of feeding rate control can be further improved, and the technical problem of low preparation efficiency in the process of preparing 1, 3-propanediol by fermentation in the prior art is effectively solved.
It will be appreciated by those skilled in the art that embodiments of the invention may be provided as a method, system, or computer product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the invention may take the form of a computer product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer instructions. These computer instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
While preferred embodiments of the present invention have been described, additional variations and modifications in those embodiments may occur to those skilled in the art once they learn of the basic inventive concepts. It is therefore intended that the following claims be interpreted as including the preferred embodiments and all such alterations and modifications as fall within the scope of the invention.
It will be apparent to those skilled in the art that various modifications and variations can be made to the present invention without departing from the spirit or scope of the invention. Thus, it is intended that the present invention also include such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.
Claims (10)
1. A fermentation control method is characterized by being applied to a process of preparing 1, 3-propanediol from fermentation liquor; the fermentation control method comprises the following steps:
continuously acquiring an alkali adding change parameter for the fermentation broth, determining a basic feeding rate for a substrate of the fermentation broth according to the alkali adding change parameter when the alkali adding change parameter meets a first preset condition, and carrying out feeding operation based on the basic feeding rate;
when the fermentation time of the fermentation liquid is within a preset period, acquiring a first fermentation state parameter of the fermentation liquid within the current preset period, and processing the first fermentation state parameter by utilizing a preset predictive control algorithm to obtain a theoretical feeding change parameter of the fermentation substrate within the current preset period;
and controlling the feeding operation of the fermentation liquor in the next preset period by utilizing the theoretical feeding variation parameters.
2. The method of claim 1, wherein the alkalization variation parameter satisfies a first preset condition, comprising:
and the alkali adding change parameter is monotonically decreased along with the time change within the preset time.
3. The method of claim 2, wherein determining a base flow rate for the broth substrate based on the alkalization variation parameter comprises:
determining the basic flow rate based on the product of the alkali addition change parameter and a first preset coefficient; wherein the magnitude of the first preset coefficient is positively correlated with the volume of the fermentation broth.
4. The method of claim 3, wherein the first fermentation status parameter comprises: the concentration of biomass in the fermentation broth, the concentration of substrate in the fermentation broth, the concentration of 1, 3-propanediol in the fermentation broth, the concentration of acidic species in the fermentation broth, and the volume of the fermentation broth;
the processing the first fermentation state parameter by using a preset predictive control algorithm to obtain a theoretical feeding change parameter of the fermentation substrate in a current preset period comprises the following steps:
taking the concentration of 1, 3-propylene glycol in the fermentation liquor as an objective function, and taking the biomass concentration in the fermentation liquor, the concentration of a substrate in the fermentation liquor, the concentration of an acidic substance in the fermentation liquor, the volume of the fermentation liquor and the theoretical feeding rate in the current preset period as constraint conditions of the objective function;
and determining a corresponding theoretical feeding rate in the process of calculating the extremum of the objective function so as to obtain the theoretical feeding change parameter based on a plurality of theoretical feeding rates.
5. The method of claim 4, further comprising, after said controlling the feeding of the fermentation broth during a next preset period using the theoretical feeding variation parameter:
acquiring a second fermentation state parameter of the fermentation liquid in the next preset period, and processing the second fermentation state parameter and the theoretical feed supplement change parameter by utilizing the predictive control algorithm to obtain a theoretical concentration change parameter of 1, 3-propanediol in the fermentation liquid in the next preset period;
and correcting the predictive control algorithm by using the actual concentration variation parameter of the 1, 3-propanediol in the first fermentation state parameter and the theoretical concentration variation parameter.
6. The method of claim 5, wherein modifying the predictive control algorithm using the actual concentration variation parameter of 1, 3-propanediol in the first fermentation status parameter and the theoretical concentration variation parameter comprises:
correspondingly modifying a second preset coefficient in the predictive control algorithm based on the difference value between the actual concentration variation parameter and the theoretical concentration variation parameter;
and if the difference value between the actual concentration variation parameter and the theoretical concentration variation parameter is larger, the second preset coefficient is smaller.
7. The method of claim 5, further comprising, after said modifying said predictive control algorithm:
circularly executing the step of correspondingly correcting a second preset coefficient in the predictive control algorithm based on the difference value between the actual concentration variation parameter and the theoretical concentration variation parameter;
until the number of cycles is greater than a first preset threshold, or until the difference between the actual concentration variation parameter and the theoretical concentration variation parameter is less than a second preset threshold.
8. A fermentation control device is characterized by being applied to a process of preparing 1, 3-propanediol from fermentation liquor; the fermentation control device comprises:
the data calculation unit is used for continuously acquiring an alkali adding change parameter aiming at the fermentation broth, determining a basic feeding rate aiming at a substrate of the fermentation broth according to the alkali adding change parameter when the alkali adding change parameter meets a first preset condition, and carrying out feeding operation based on the basic feeding rate;
the simulation prediction unit is used for acquiring a first fermentation state parameter of the fermentation broth in a current preset period when the fermentation time of the fermentation broth is in the preset period, and processing the first fermentation state parameter by utilizing a preset prediction control algorithm to acquire a theoretical feeding change parameter of the fermentation substrate in the current preset period;
and the feeding control unit is used for controlling the feeding operation of the fermentation liquid in the next preset period by utilizing the theoretical feeding change parameter.
9. A fermentation control apparatus comprising a memory, a processor and code stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1-7 when executing the code.
10. A readable storage medium having stored thereon a program, which when executed by a processor, implements the method of any of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310658351.0A CN116640886A (en) | 2023-06-05 | 2023-06-05 | Fermentation control method, device, equipment and readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310658351.0A CN116640886A (en) | 2023-06-05 | 2023-06-05 | Fermentation control method, device, equipment and readable storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116640886A true CN116640886A (en) | 2023-08-25 |
Family
ID=87615196
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310658351.0A Pending CN116640886A (en) | 2023-06-05 | 2023-06-05 | Fermentation control method, device, equipment and readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116640886A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117448147A (en) * | 2023-10-19 | 2024-01-26 | 威海紫光优健科技股份有限公司 | Process control system for converting red ginseng into Rg3 based on microorganisms and preparation method |
-
2023
- 2023-06-05 CN CN202310658351.0A patent/CN116640886A/en active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117448147A (en) * | 2023-10-19 | 2024-01-26 | 威海紫光优健科技股份有限公司 | Process control system for converting red ginseng into Rg3 based on microorganisms and preparation method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110054274B (en) | Water purification flocculation precipitation dosing control method | |
| CN108153145B (en) | PID and filtering algorithm based train speed automatic control method | |
| CN116640886A (en) | Fermentation control method, device, equipment and readable storage medium | |
| CN110101106B (en) | Moisture control method and system for dampening and humidifying process based on fuzzy feedforward feedback algorithm | |
| JP7059989B2 (en) | Control system and control method | |
| US20090117647A1 (en) | Method for on-line optimization of a fed-batch fermentation unit to maximize the product yield | |
| WO2007085880A1 (en) | A method for on-line prediction of future performance of a fermentation unit. | |
| CN101825870A (en) | Method and system for controlling supply quantity of water-treatment flocculating agent | |
| CN109960149B (en) | Automatic disturbance rejection control parameter setting method for pH neutralization process | |
| De Franco et al. | Bayesian learning for the Markowitz portfolio selection problem | |
| CN106605179B (en) | Predicted value shaping system and method, and control system and control method | |
| CN103618486A (en) | Fuzzy-control direct-current motor speed control method | |
| CN101995845A (en) | Field programmable gate array (FPGA)-based auto-leveling control system and method | |
| Li et al. | Next day load forecasting using SVM | |
| Vargas et al. | On-line maximization of biogas production in an anaerobic reactor using a pseudo-super-twisting controller | |
| CN114660924A (en) | Servo motor energy-saving control system and method | |
| CN113721468B (en) | Self-healing control method and system for additive faults of sewage treatment system actuator | |
| Gorrini et al. | Reaction rate estimation and model predictive control of hybridoma cell cultures | |
| CN109870986B (en) | On-line control method of stirred tank reactor based on neural network and data drive | |
| JPH0654678A (en) | Method for automatic feed culture of microorganism | |
| CN102953966B (en) | Self-adaptive iterative control method of variable-frequency pumps in pump station group | |
| Kavanagh et al. | Productivity improvement of recombinant Escherichia coli fermentation via robust optimization | |
| CN105199973B (en) | A kind of Yeast Cultivation online adaptive control method based on differential evolution algorithm | |
| Hartwich et al. | Adaptive control vector parameterization for nonlinear model‐predictive control | |
| CN201867815U (en) | Water quality control system with accidental modeling and data drive |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |