CN110427733B - Method for obtaining algae concentration based on phosphorus cycle - Google Patents
Method for obtaining algae concentration based on phosphorus cycle Download PDFInfo
- Publication number
- CN110427733B CN110427733B CN201910857739.7A CN201910857739A CN110427733B CN 110427733 B CN110427733 B CN 110427733B CN 201910857739 A CN201910857739 A CN 201910857739A CN 110427733 B CN110427733 B CN 110427733B
- Authority
- CN
- China
- Prior art keywords
- algae
- rate
- above formula
- follows
- growth
- 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
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
Abstract
The invention discloses an algae concentration obtaining method based on phosphorus cycle. The method can effectively simulate the growth and death process of the algae, so that the method has the capability of simulating the change of the algae, and effectively reduces the number of unknown biochemical parameters by using the variables of phosphorus and key nutrient salts. The method can effectively simulate the change process of the algae in the water body, not only does not reduce the simulation precision, but also greatly reduces the model parameters and increases the calculation efficiency because phosphorus is selected as key nutrient-limiting salt.
Description
Technical Field
The invention relates to the technical field of algae simulation, in particular to an algae concentration acquisition method based on phosphorus circulation.
Background
At present, excessive discharge of nutrient salts causes the eutrophication degree of lakes and reservoirs to be increased, and water bloom is frequent. Numerical simulation is one of effective methods for analyzing the time-space change of the water bloom in the water body. The existing algae calculation models are numerous and can meet the requirements of bloom simulation, but the models have the following defects:
a. the simulated nutrient salts are numerous and cannot represent the core limit of key nutrient salts;
b. the simulation of the biochemical process is excessive, so that a plurality of ecological parameters which cannot be directly obtained are generated, and the precision of a simulation result is seriously influenced;
c. the operation efficiency is low.
Disclosure of Invention
Aiming at the defects in the prior art, the method for obtaining the algae concentration based on the phosphorus cycle solves the problem of inaccurate analog calculation of the algae concentration.
In order to achieve the purpose of the invention, the invention adopts the technical scheme that: a method for obtaining algae concentration based on phosphorus cycle comprises the following steps:
s1, establishing a dynamic model of the biological debris;
s2, establishing a kinetic model of the nutrient salt according to the kinetic model of the biological debris;
s3, establishing a dynamic model of algae according to the dynamic model of the biological debris and the dynamic model of the nutrient salt;
and S4, calculating according to the dynamic model of the algae to obtain the algae concentration.
Further, the method comprises the following steps: the dynamic model of the biological debris in the step S1 is as follows:
in the above formula, DET is the biomass of the chips, t is the time,to be the rate of change of the concentration of debris,to account for the rate of change of debris due to algae death,snk (DET) is the rate of debris sedimentation;
wherein the rate of mineralised hydrolysis of the finesThe calculation formula of (2) is as follows:
in the above formula, r DET The crumb based hydrolysis rate.
Further: the kinetic model of the nutritive salt in the step S2 is as follows:
in the above formula, PO 4 Is the nutrient salt biomass, t is the time,in order to determine the rate of change of the nutrient salt concentration,the rate of change of nutrient salts caused by the growth of algae,the rate of change of nutrient salts caused by hydrolysis of the crumb;
wherein the rate of change of nutrient salts caused by algae growthThe calculation formula of (c) is:
rate of change of nutrient salts by hydrolysis of debrisThe calculation formula of (2) is as follows:
in the above-mentioned formula, the compound has the following structure,the rate of mineralized hydrolysis of the fines.
Further: the dynamic model of the algae in the step S3 is as follows:
in the above formula, P is the biomass of algae, t is the time,is a model of the dynamics of the algae,is the rate at which the algae grow and,is the rate of algal death.
Further, the method comprises the following steps: the rate of growth of the algaeThe calculation formula of (2) is as follows:
in the above formula, r P P is the algal biomass, the algal basal growth rate.
Further: the basal growth rate r of the algae P The calculation formula of (2) is as follows:
in the above formula, r 0 Is the maximum growth rate of the algae,in order to be a function of the temperature limit,in order to be a function of the light confinement,as a nutrient salt limiting function;
in the above formula, T is water temperature, T 1 Lower temperature limit, T, for proper algae growth 2 The upper limit of temperature suitable for the growth of algae;
in the above formula, I Z Solar radiation at a depth Z below the water surface, I opt Optimum illumination intensity for algae growth, wherein the depth under water surface is the solar radiation I at Z Z The calculation formula of (c) is:
I Z =I 0 e -(KW+KC·P)D
in the above formula, I 0 The water surface solar radiation is adopted, KW is the extinction coefficient of all substances except algae, KC is an algae extinction coefficient factor, and D is the water depth;
in the above formula, K P At a phosphate half-saturation concentration, PO 4 Is the concentration of nutrient salt.
in the above formula,. Mu. 1 First order death rate for algae, μ 2 Second order death rate for algae.
Further: the calculation formula of the algae concentration in the step S4 is:
in the above formula, P n+1 The next time point is the algae concentration, P n At is the algae concentration at the current time, and Δ t is the calculation step length.
The invention has the beneficial effects that: the method of the invention effectively reduces the number of unknown biochemical parameters by using the variable of phosphorus and key nutritive salt, and can accurately calculate the concentration of the algae.
The method can effectively simulate the change process of the algae in the water body, not only does not reduce the simulation precision, but also greatly reduces the model parameters and increases the calculation efficiency because phosphorus is selected as key nutrient limiting salt, and can effectively simulate the growth and death processes of the algae.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a diagram of the interaction process of nutritive salt-algae-bio-debris in an embodiment of the present invention;
FIG. 3 is a schematic representation of algae change over time in an embodiment of the present invention.
Detailed Description
The following description of the embodiments of the present invention is provided to facilitate the understanding of the present invention by those skilled in the art, but it should be understood that the present invention is not limited to the scope of the embodiments, and it will be apparent to those skilled in the art that various changes may be made without departing from the spirit and scope of the invention as defined and defined in the appended claims, and all matters produced by the invention using the inventive concept are protected.
Nutrient salt-algae-bio-debris three-variable circulation mode, as shown in FIG. 2, exogenous input to nutrient salt (phosphate PO) 4 ) The algae (phytoplankton P) ingests nutrient salts (phosphate PO) 4 ) The algae (phytoplankton P) die to form biological debris, and the biological debris DET decomposes to form nutrient salt (phosphate PO) 4 )。
As shown in fig. 1, a method for obtaining algae concentration based on phosphorus cycle includes the following steps:
s1, establishing a dynamic model of the biological debris;
debris is primarily organisms of phytoplankton residues, the kinetics of which are affected by algal death and settlement, so the kinetic model for biological debris is:
in the above formula, DET is the biomass of debris, t is the time,to be the rate of change of the biomass of the crumb,for the rate of change of debris due to algae death,snk (DET) is the rate of debris sedimentation;
wherein the rate of mineralised hydrolysis of the finesThe calculation formula of (2) is as follows:
in the above formula, r DET Is the crumb based hydrolysis rate.
S2, establishing a kinetic model of the nutrient salt according to the kinetic model of the biological debris;
phosphorus plays an important role in the growth of algae. The circulation of phosphorus is mainly based on phytoplankton absorption and debris hydrolysis, so the kinetic model of nutritive salt is as follows:
in the above formula, PO 4 Is the nutrient salt biomass, t is the time,in order to determine the rate of change of the nutrient salt concentration,the rate of change of nutrient salts caused by the growth of algae,the rate of change of nutrient salts caused by hydrolysis of the crumb;
algae absorb phosphorus for growth through photosynthesis, so the rate of nutrient salt change caused by algae growthThe calculation formula of (2) is as follows:
the algae debris is converted into phosphate through hydrolysis and mineralization processes, so that the change rate of nutrient salt caused by the hydrolysis of the debrisThe calculation formula of (c) is:
in the above formula, the first and second carbon atoms are,the rate of mineralized hydrolysis of the fines.
S3, establishing a dynamic model of the algae according to the dynamic model of the biological debris and the dynamic model of the nutrient salt;
water quality models rarely model specific species of algae. Algae may be grouped into one or several groups. The kinetic model of algae is influenced by growth and death, and is as follows:
in the above formula, P is the biomass of algae, t is the time,is a model of the dynamics of the algae,is the rate of growth of algae,Is the rate of algae death.
The growth of algae is influenced by water temperature, nutrient salts, light and basal growth rate, so the growth rate of algaeThe calculation formula of (2) is as follows:
in the above formula, r P P is the algal biomass, the algal basal growth rate.
Basal growth rate r of algae P The calculation formula of (2) is as follows:
in the above formula, r 0 The maximum growth rate of the algae is set as,in order to be a function of the temperature limit,in order to be a function of the light confinement,as a nutrient salt limiting function;
the growth of the algae is influenced by the temperature, and the algae grows fastest in a certain proper temperature interval. Above or below this temperature, the growth rate of the algae is affected. Temperature limit functionThe calculation formula of (c) is:
in the above formula, T is water temperature, T 1 Lower temperature limit, T, for proper algae growth 2 The upper limit of temperature suitable for the growth of algae;
the growth of algae is also influenced by light (solar radiation). Light limiting functionThe calculation formula of (2) is as follows:
in the above formula, I Z Solar radiation at a depth Z below the water surface, I opt Optimum illumination intensity for algae growth, wherein the depth under water surface is the solar radiation I at Z Z The calculation formula of (2) is as follows:
I Z =I 0 e -(KW+KC·P)D
in the above formula, I 0 The water surface solar radiation is adopted, KW is the extinction coefficient of all substances except algae, KC is an algae extinction coefficient factor, and D is the water depth;
the limiting nutrient salts that algae growth is subject to are mainly nitrogen, phosphorus and silicon. The mass ratio of nitrogen and phosphorus in domestic water is generally more than 15, and the water is a typical phosphorus-limited water body. Therefore, phosphorus in the water body is key nutrient salt for limiting the growth of algae, and a Michaelis-Menton model is selected to describe the relationship between the growth of algae and the phosphorus concentration and a nutrient salt limiting functionThe calculation formula of (c) is:
in the above formula, K P At a phosphate half-saturation concentration, PO 4 Is the nutrient salt concentration.
Death of algae is a nonlinear process, so the rate of algae deathThe calculation formula of (2) is as follows:
in the above formula,. Mu. 1 First order death rate for algae, μ 2 Second order death rate for algae.
And S4, calculating according to the dynamic model of the algae to obtain the algae concentration.
The calculation formula of the algae concentration is as follows:
in the above formula, P n+1 The next time point is the algae concentration, P n At the current time, the algae concentration is calculated as Δ t.
The values of the parameters in the present invention are shown in table 1.
TABLE 1 ecological parameters and suggested values
As can be seen from the data in FIG. 3, the method of the present invention effectively reduces the number of unknown biochemical parameters with the variables of phosphorus and key nutrient salts, and allows accurate calculation of algae concentration.
The method can effectively simulate the change process of the algae in the water body, not only does not reduce the simulation precision, but also greatly reduces the model parameters and increases the calculation efficiency because phosphorus is selected as key nutrient limiting salt, and can effectively simulate the growth and death processes of the algae.
Claims (4)
1. A method for obtaining algae concentration based on phosphorus cycle is characterized by comprising the following steps:
s1, establishing a dynamic model of the biological debris;
s2, establishing a kinetic model of the nutrient salt according to the kinetic model of the biological debris;
s3, establishing a dynamic model of algae according to the dynamic model of the biological debris and the dynamic model of the nutrient salt;
s4, calculating according to the dynamic model of the algae to obtain the concentration of the algae;
the dynamic model of the biological debris in the step S1 is as follows:
in the above formula, DET is the biomass of the chips, t is the time,to be the rate of change of the biomass of the crumb,to account for the rate of change of debris due to algae death,snk (DET) is the rate of debris sedimentation;
wherein the rate of mineralisation and hydrolysis of the debrisThe calculation formula of (2) is as follows:
in the above formula, r DET Is the crumb based hydrolysis rate;
the kinetic model of the nutritive salt in the step S2 is as follows:
in the above formula, PO 4 Is the nutrient salt biomass, t is the time,the rate of change of the nutrient salt concentration is,the rate of change of nutrient salts caused by the growth of algae,the rate of change of nutrient salts caused by hydrolysis of the crumb;
wherein the rate of change of nutrient salts caused by algae growthThe calculation formula of (2) is as follows:
in the above-mentioned formula, the compound has the following structure,is the rate of mineralized hydrolysis of the fines;
the dynamic model of the algae in the step S3 is as follows:
in the above formula, P is the biomass of algae, t is the time,is a model of the dynamics of the algae,is the rate at which the algae grow and,is the rate of algal death;
the calculation formula of the algae concentration in the step S4 is as follows:
in the above formula, P n+1 The next time point is the algae concentration, P n At is the algae concentration at the current time, and Δ t is the calculation step length.
3. The method of claim 2, wherein the basal growth rate r of the algae is set according to the algae concentration P The calculation formula of (2) is as follows:
in the above formula, r 0 The maximum growth rate of the algae is set as,in order to be a function of the temperature limit,in order to be a function of the light confinement,as a nutrient salt limiting function;
in the above formula, T is water temperature, T 1 Lower temperature limit, T, for proper algae growth 2 The upper limit of temperature suitable for the growth of algae;
in the above formula, I Z Solar radiation at a depth Z below the water surface, I opt Optimum illumination intensity for algae growth, wherein the depth under water surface is the solar radiation I at Z Z The calculation formula of (2) is as follows:
I Z =I 0 e -(KW+KC·P)D
in the above formula, I 0 The water surface solar radiation is adopted, KW is the extinction coefficient of all substances except algae, KC is an algae extinction coefficient factor, and D is the water depth;
in the above formula, K P At a phosphate half-saturation concentration, PO 4 Is the nutrient salt concentration.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910857739.7A CN110427733B (en) | 2019-09-09 | 2019-09-09 | Method for obtaining algae concentration based on phosphorus cycle |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910857739.7A CN110427733B (en) | 2019-09-09 | 2019-09-09 | Method for obtaining algae concentration based on phosphorus cycle |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110427733A CN110427733A (en) | 2019-11-08 |
CN110427733B true CN110427733B (en) | 2022-11-29 |
Family
ID=68418973
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910857739.7A Active CN110427733B (en) | 2019-09-09 | 2019-09-09 | Method for obtaining algae concentration based on phosphorus cycle |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110427733B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112385536B (en) * | 2020-10-30 | 2022-05-06 | 南京信息工程大学 | Method for measuring and calculating light-limited depth of large floating algae |
CN115114867B (en) * | 2022-07-11 | 2023-07-11 | 中国水利水电科学研究院 | Emergent aquatic plant biomass prediction method |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5863434A (en) * | 1994-12-14 | 1999-01-26 | University Of Ottawa/Universite D'ottawa | Psychrophilic anaerobic treatment of waste in a sequencing semibatch/batch bioreactor |
CN103942358A (en) * | 2014-03-07 | 2014-07-23 | 同济大学 | Optimized design method of photobioreactor for cultivating microalgae in waste water |
CN106345400A (en) * | 2016-10-27 | 2017-01-25 | 中国科学院城市环境研究所 | Porous phosphorus removal adsorbent based on hydrated iron oxide and preparation method thereof |
CN107292436A (en) * | 2017-06-16 | 2017-10-24 | 北京工商大学 | Blue-green alga bloom Forecasting Methodology based on nonlinear kinetics temporal model |
CN109255168A (en) * | 2018-08-27 | 2019-01-22 | 广州资源环保科技股份有限公司 | A kind of method of predictions and simulations Shallow Lake Ecosystems total phosphorus circulation |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5011675B2 (en) * | 2005-08-09 | 2012-08-29 | 味の素株式会社 | Method for simulating material production process |
US20120058248A1 (en) * | 2008-12-17 | 2012-03-08 | LiveFuels, Inc. | Systems and methods for reducing algal biomass |
GB201001957D0 (en) * | 2010-02-05 | 2010-03-24 | Biotarisk As | Pollution monitoring |
WO2013170031A1 (en) * | 2012-05-09 | 2013-11-14 | The Regents Of The University Of California | Method for in silico modeling of gene product expression and metabolism |
FR2999299B1 (en) * | 2012-12-12 | 2021-05-07 | Ifp Energies Now | METHOD OF EXPLOITATION OF A SEDIMENTARY BASIN BY MEANS OF A STRATIGRAPHIC SIMULATION COUPLED WITH A MODEL OF PRODUCTION AND DEGRADATION OF ORGANIC MATTER |
-
2019
- 2019-09-09 CN CN201910857739.7A patent/CN110427733B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5863434A (en) * | 1994-12-14 | 1999-01-26 | University Of Ottawa/Universite D'ottawa | Psychrophilic anaerobic treatment of waste in a sequencing semibatch/batch bioreactor |
CN103942358A (en) * | 2014-03-07 | 2014-07-23 | 同济大学 | Optimized design method of photobioreactor for cultivating microalgae in waste water |
CN106345400A (en) * | 2016-10-27 | 2017-01-25 | 中国科学院城市环境研究所 | Porous phosphorus removal adsorbent based on hydrated iron oxide and preparation method thereof |
CN107292436A (en) * | 2017-06-16 | 2017-10-24 | 北京工商大学 | Blue-green alga bloom Forecasting Methodology based on nonlinear kinetics temporal model |
CN109255168A (en) * | 2018-08-27 | 2019-01-22 | 广州资源环保科技股份有限公司 | A kind of method of predictions and simulations Shallow Lake Ecosystems total phosphorus circulation |
Non-Patent Citations (2)
Title |
---|
Modeling of P-dynamics and algal growth in a stratified reservoir—mechanisms of P-cycle in water and interaction between overlying water and sediment;Eiji Komatsu et al;《ecological modelling 》;20061231;第331-349页 * |
长江口多藻类生态动力学模型研究;王彪 等;《环境生态学》;20190630;第1卷(第2期);第12-16页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110427733A (en) | 2019-11-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Marsullo et al. | Dynamic modeling of the microalgae cultivation phase for energy production in open raceway ponds and flat panel photobioreactors | |
Malapascua et al. | Photosynthesis monitoring to optimize growth of microalgal mass cultures: application of chlorophyll fluorescence techniques | |
Salgueiro et al. | Bioremediation of wastewater using Chlorella vulgaris microalgae: Phosphorus and organic matter | |
CN110427733B (en) | Method for obtaining algae concentration based on phosphorus cycle | |
Viruela et al. | Water resource recovery by means of microalgae cultivation in outdoor photobioreactors using the effluent from an anaerobic membrane bioreactor fed with pre-treated sewage | |
Kesaano et al. | Dissolved inorganic carbon enhanced growth, nutrient uptake, and lipid accumulation in wastewater grown microalgal biofilms | |
Wu et al. | An integrated microalgal growth model and its application to optimize the biomass production of Scenedesmus sp. LX1 in open pond under the nutrient level of domestic secondary effluent | |
Buono et al. | Productivity and biochemical composition of Tetradesmus obliquus and Phaeodactylum tricornutum: effects of different cultivation approaches | |
del Rio‐Chanona et al. | Dynamic modeling of green algae cultivation in a photobioreactor for sustainable biodiesel production | |
Znad et al. | CO2 biomitigation and biofuel production using microalgae: photobioreactors developments and future directions | |
CN108546648A (en) | A kind of environment-friendly type microalgae culture method | |
Kawaroe et al. | Comparison of density, specific growth rate, biomass weight, and doubling time of microalgae Nannochloropsis sp. cultivated in Open Raceway Pond and Photobioreactor | |
Quinn et al. | Scale‐Up of flat plate photobioreactors considering diffuse and direct light characteristics | |
CN109097283A (en) | A kind of method of microalgae alkalinity flocculation harvest and Cyclic culture | |
Naderi et al. | Modified photobioreactor for biofixation of carbon dioxide by Chlorella vulgaris at different light intensities | |
Concas et al. | Nannochloris eucaryotum growth: Kinetic analysis and use of 100% CO2 | |
Kroon et al. | Modelling microalgal productivity in a high rate algal pond based on wavelength dependent optical properties | |
DE69216003D1 (en) | METHOD FOR PRODUCING TAXOL BY MEANS OF TAXUS SPECIES CELL CULTURE | |
US20220220527A1 (en) | Method for measuring the activity of a culture of microalgae | |
Yang et al. | A mathematical model of neutral lipid content in terms of initial nitrogen concentration and validation in Coelastrum sp. ha-1 and application in Chlorella sorokiniana | |
Graneli et al. | Can microbenthic photosynthesis influence below-halocline oxygen conditions in the Kattegat? | |
CN103757064B (en) | A kind of method of pH regulation and control microalgae grease Rapid Accumulation | |
Romagnoli et al. | Growth kinetic model for microalgae cultivation in open raceway ponds: a system dynamics tool | |
Esteves et al. | Comparative assessment of microalgal growth kinetic models based on light intensity and biomass concentration | |
Wang et al. | Effects of Different Conditions on the Growth of Microcystis aeruginosa. |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |