[go: up one dir, main page]

CN112485399B - Method for Estimating Toxin Concentration in Wastewater - Google Patents

Method for Estimating Toxin Concentration in Wastewater Download PDF

Info

Publication number
CN112485399B
CN112485399B CN202011506021.2A CN202011506021A CN112485399B CN 112485399 B CN112485399 B CN 112485399B CN 202011506021 A CN202011506021 A CN 202011506021A CN 112485399 B CN112485399 B CN 112485399B
Authority
CN
China
Prior art keywords
sewage
concentration
toxins
toxin
estimating
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
Application number
CN202011506021.2A
Other languages
Chinese (zh)
Other versions
CN112485399A (en
Inventor
赵顺毅
赵珂
王劲夫
刘飞
栾小丽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Saida Environmental Technology Co ltd
Original Assignee
Jiangnan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jiangnan University filed Critical Jiangnan University
Priority to CN202011506021.2A priority Critical patent/CN112485399B/en
Publication of CN112485399A publication Critical patent/CN112485399A/en
Application granted granted Critical
Publication of CN112485399B publication Critical patent/CN112485399B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/18Water
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING 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
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/02Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving viable microorganisms
    • C12Q1/04Determining presence or kind of microorganism; Use of selective media for testing antibiotics or bacteriocides; Compositions containing a chemical indicator therefor
    • C12Q1/06Quantitative determination
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Organic Chemistry (AREA)
  • Wood Science & Technology (AREA)
  • Biochemistry (AREA)
  • Theoretical Computer Science (AREA)
  • Immunology (AREA)
  • Computing Systems (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Analytical Chemistry (AREA)
  • Zoology (AREA)
  • Biophysics (AREA)
  • Biotechnology (AREA)
  • General Engineering & Computer Science (AREA)
  • Toxicology (AREA)
  • Genetics & Genomics (AREA)
  • Microbiology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Molecular Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

本发明公开了一种污水毒素浓度估计方法,将具有一定细胞指数CI的培养细胞暴露在污水样本中,连续采集污水样本内培养细胞的细胞数量测量值;根据细胞对毒素的摄取机理,以及毒素所引起的细胞数量变化创建细胞毒素动态模型;使用粒子滤波算法对细胞毒素动态模型和细胞数量测量值进行处理,估计出污水样本中的毒素浓度。本发明采用粒子滤波算法根据暴露在污水样本中的培养细胞的细胞数量测量值和细胞毒素动态模型对污水中的毒素浓度做连续在线估计,并准确估计出污水中毒素浓度,而采用的粒子滤波算法对于存在大的过程噪声、观测噪声或建模不确定的非线性动态系统可以取得良好的估计效果,保证毒素浓度估计的准确性。

Figure 202011506021

The invention discloses a method for estimating the concentration of toxins in sewage. The cultured cells with a certain cell index CI are exposed to sewage samples, and the measured value of the cell number of the cultured cells in the sewage samples is continuously collected; The resulting change in cell number creates a cytotoxin dynamic model; the cytotoxin dynamic model and cell number measurements are processed using a particle filter algorithm to estimate the toxin concentration in the sewage sample. The present invention adopts the particle filter algorithm to continuously estimate the toxin concentration in the sewage on-line according to the cell number measurement value of the cultured cells exposed in the sewage sample and the cytotoxin dynamic model, and accurately estimates the toxin concentration in the sewage. The algorithm can achieve good estimation results for nonlinear dynamic systems with large process noise, observation noise or modeling uncertainty, and ensure the accuracy of toxin concentration estimation.

Figure 202011506021

Description

Sewage toxin concentration estimation method
Technical Field
The invention relates to a toxin concentration determination method, in particular to a sewage toxin concentration estimation method.
Background
Sewage treatment is a process of separating pollutants contained in sewage or converting the pollutants into harmless substances by using various different methods to purify the sewage. Generally, urban and industrial water treatment units combine primary, secondary and tertiary treatment processes to efficiently treat various pollutants in water, i.e., inorganic pollutants, organic pollutants and biological pollutants. The first-stage treatment mainly adopts a physical method to remove suspended matters and partial organic matters; the secondary treatment is to remove colloid and soluble substances in the sewage by adopting a biological method, a chemical method or a physical and chemical method on the basis of the primary treatment; and the third-stage treatment further removes refractory substances, nitrogen, phosphorus, heavy metals and the like. Due to the complexity of wastewater, it is often difficult to determine the specific poison that causes the performance degradation of a wastewater treatment plant, and therefore, timely measurement of toxicity is critical to the smooth operation of a wastewater treatment plant.
The basic premise of a successful early warning system is to continuously collect accurate data describing the risk of toxic contamination, but most of the key quality indicators need to be obtained by off-line sample analysis, which is usually expensive and requires frequent and costly maintenance, and the lack of appropriate key variable information can seriously affect the monitoring of the toxin concentration in the sewage treatment process. It is therefore important to use real-time measurements and established models to achieve continuous online estimation of key variables.
In general, the response of cells to small doses of toxicants is limited (stable), while the response of large doses of toxicants is unlimited (unstable), a phenomenon that clearly indicates the non-linearity of cytotoxicity. Common estimation methods for nonlinear systems, such as Extended Kalman Filtering (EKF), Unscented Kalman Filtering (UKF), ensemble kalman filtering (EnKF), rolling time domain estimation (MHE), are suitable for nonlinear dynamic systems with small uncertainty, but are difficult to achieve good effects in nonlinear dynamic systems with large process noise, observation noise, or modeling uncertainty.
Disclosure of Invention
The invention aims to provide a sewage toxin concentration estimation method, which can accurately estimate the concentration of toxins in sewage.
In order to solve the above technical problems, the present invention provides a method for estimating the concentration of toxins in wastewater, comprising the steps of,
exposing cultured cells with a certain cell index CI in a sewage sample, and continuously collecting cell number measured values of the cultured cells in the sewage sample;
creating a cytotoxin dynamic model according to the uptake mechanism of the toxin by the cells and the change of the cell number caused by the toxin;
and processing the cytotoxin dynamic model and the cell number measured value by using a particle filtering algorithm to estimate the toxin concentration in the sewage sample.
In a preferred embodiment of the present invention, further comprising, creating the cytotoxin kinetic model comprises,
creating a kinetic model of the toxin uptake process by the cells:
Figure BDA0002844944300000021
creating an apoptotic mechanism for intracellular endotoxin concentration:
Figure BDA0002844944300000022
cicharacterizing the intracellular endotoxin concentration; c. CeCharacterizing the extracellular toxin concentration; n represents the cell number, and N is approximately equal to CI; k is a radical of1,k2,k3,KiIs a relevant parameter; kciCharacterizing the relative cell killing rate in the presence of the toxin; k is a radical ofsThe relative cell proliferation rate was characterized in the absence of toxin.
In a preferred embodiment of the present invention, further comprising, creating the cytotoxin kinetic model comprises,
exposing the cultured cells with a certain cell index CI to toxins with different concentrations, and continuously collecting cell number measured values of the cultured cells under the action of the toxins;
substituting the cell number measured value of the cultured cells under the action of the toxin into the (formula 1) and the (formula 2), identifying the unknown parameters in the kinetic model and the cell endotoxin concentration apoptosis mechanism by adopting a nonlinear regression method, and calculating the related parameter k1,k2,k3,Ki
In a preferred embodiment of the present invention, the cultured cells further comprise a cell index CI of 1.0 to 1.2 for cells exposed to different concentrations of the toxin.
In a preferred embodiment of the present invention, the cultured cells exposed to the wastewater sample have a cell index CI of 1.0 to 1.2.
In a preferred embodiment of the present invention, the cultured cell is further cultured in the cell line NIH 3T 3.
In a preferred embodiment of the present invention, the initial cell number of cultured cells exposed to the contaminated water sample and cultured cells exposed to different concentrations of toxin each comprises 10000 cells per well.
In a preferred embodiment of the present invention, the method for estimating the concentration of the toxin in wastewater further comprises discretizing the dynamic model of the toxin in wastewater to obtain a discretized dynamic model of the toxin in wastewater
xk=f(xk-1)+wk
yk=g(xk)+vk
Figure BDA0002844944300000031
Is a state vector;
Figure BDA0002844944300000032
is a measurement vector;
Figure BDA0002844944300000033
to relate to xkA non-linear function of (d);
Figure BDA0002844944300000034
is an observation function;
Figure BDA0002844944300000035
is process noise and is subject to distribution
Figure BDA0002844944300000036
Figure BDA0002844944300000037
To observe noise and obey distribution
Figure BDA0002844944300000038
And the particle filtering algorithm is combined with the discretized cytotoxin dynamic model to process the measured value of the cell number in the sewage sample.
In a preferred embodiment of the present invention, the particle filter algorithm further comprises a function of defining a state variable x ═ c for processing the measured cell number in the wastewater samplei,ce,N]Defining a measurement vector y ═ N, based on a discretized dynamic model of the cytotoxin, and estimating the true intracellular endotoxin concentration c during the cytotoxic process from the noisy measurements of the number of cellsiExtracellular toxin concentration ceAnd the number of cells N.
In a preferred embodiment of the present invention, the discretization method further comprises one of an euler method, a longge-kutta method, and a linear multi-step method.
The invention has the beneficial effects that:
the method for estimating the concentration of the toxin in the sewage adopts a particle filter algorithm to continuously estimate the concentration of the toxin in the sewage on line according to the cell number measured value of the cultured cells exposed in the sewage sample and a cytotoxin dynamic model, and accurately estimates the concentration of the toxin in the sewage. The adopted particle filtering algorithm can obtain a good estimation effect on a nonlinear dynamic system with large process noise, observation noise or uncertain modeling, and the accuracy of toxin concentration estimation is ensured.
Drawings
FIG. 1 is a flow chart of a method for estimating the concentration of toxins in wastewater according to a preferred embodiment of the present invention;
FIG. 2 shows the estimation of toxin concentration in wastewater using the method of the present invention at ceWhen equal to 0.62A plot of estimated toxin chromium (VI) concentration estimates versus true values;
FIG. 3 shows the estimation of toxin concentration in wastewater using the method of the present invention at ceAn error plot of estimated toxin chromium (VI) concentration at 0.62;
FIG. 4 shows the estimation of toxin concentration in wastewater using the method of the present invention at ceA plot of estimated toxin mercury (II) chloride concentration versus true at 22.35;
FIG. 5 shows the estimation of toxin concentration in wastewater using the method of the present invention at ceError plot of estimated toxin mercuric (II) chloride concentration at 48.3.
Detailed Description
The present invention is further described below in conjunction with the following figures and specific examples so that those skilled in the art may better understand the present invention and practice it, but the examples are not intended to limit the present invention.
Examples
It should be noted that (1) CI is a cell index value, which is an indirect indicator of toxic pollutants and biological consequences in the aquatic environment, and a dynamic model of the cellular toxicity response established based on the CI value can be used for online estimation of cytotoxicity.
(2) NIH 3T3 is a mouse embryo fibroblast cell line established by the national institutes of health (NIH for short), and is characterized in that the cell line is passaged once every 3 days and inoculated with 3 x 10^5cells/ml each time, and the cell line is commonly used for the research of transfection and gene expression in laboratories.
The embodiment of the invention discloses a sewage toxin concentration estimation method, which is used for estimating the concentration of toxins in sewage on line in the sewage treatment process, and is shown by referring to FIG. 1, and comprises the following steps,
in the first step, a dynamic model of the cytotoxin is created based on the mechanism of uptake of the toxin by the cells, and the change in the number of cells caused by the toxin.
(1) And obtaining the cultured cells with a certain cell index CI.
Selecting cultured cell with cell line NIH 3T3, placing the cultured cell in CO2Culturing in an incubator until the cells are culturedThe index value CI is 1.0-1.2, and the purpose of culture is to ensure that the initial cell number of the cultured cells put into use reaches 10000 cells per hole.
(2) Exposing cultured cells with a cell index CI of 1.0-1.2 in toxins with different concentrations, continuously and automatically monitoring the reaction of the cells to the toxic substances by using an RT-CES system every hour, continuously recording a CI value every hour for 24 hours, and continuously collecting cell number measurement values of the cultured cells in the toxins with different concentrations.
(3) Creating a kinetic model of the toxin uptake process by the cells:
Figure BDA0002844944300000061
and creating an apoptotic mechanism for intracellular endotoxin concentration:
Figure BDA0002844944300000062
cicharacterizing the intracellular endotoxin concentration; c. CeCharacterizing the extracellular toxin concentration; n represents the cell number, and N is approximately equal to CI; k is a radical of1,k2,k3,KiIs a relevant parameter; kciCharacterizing the relative cell killing rate in the presence of the toxin; k is a radical ofsThe relative cell proliferation rate was characterized in the absence of toxin.
k1ceRepresents a linear diffusion component;
Figure BDA0002844944300000063
representing the saturated carrier-mediated component.
Substituting the cell number measured value of the cultured cells under the action of the toxic elements obtained in the step (2) into the (formula 1) and the (formula 2), identifying unknown parameters in the kinetic model and the cell endotoxin concentration apoptosis mechanism by adopting a nonlinear regression method, and calculating the related parameter k1,k2,k3,KiAnd determining the dynamic model of the cytotoxin.
A second step of discretizing the obtained dynamic model of the cytotoxin, wherein the discretization method comprises one of an Euler method, a Runge-Kutta method and a linear multi-step method; obtaining a discretized cytotoxin dynamic model,
xk=f(xk-1)+wk
yk=g(xk)+vk
Figure BDA0002844944300000064
is a state vector;
Figure BDA0002844944300000065
is a measurement vector;
Figure BDA0002844944300000066
to relate to xkA non-linear function of (d);
Figure BDA0002844944300000067
is an observation function;
Figure BDA0002844944300000068
is process noise and is subject to distribution
Figure BDA0002844944300000069
Figure BDA00028449443000000610
To observe noise and obey distribution
Figure BDA00028449443000000611
And step three, collecting a water body sample (hereinafter referred to as a sewage sample) after sewage treatment, exposing the same cultured cells obtained in the step (1) in the sewage sample, continuously and automatically monitoring the reaction of the cells to the toxic substances every hour by using an RT-CES system, continuously recording a CI value every 24 hours, and continuously collecting the cell number measured value of the cultured cells in the sewage sample.
And step four, processing the discretized cytotoxin dynamic model and the cell number measured value obtained in the step three by using a particle filtering algorithm, and estimating the toxin concentration in the sewage sample, wherein the specific method comprises the following steps:
define state variable x ═ ci,ce,N]A measurement vector y is defined as N.
S01, initializing: from the prior distribution p (x)0) Generating initial particles
Figure BDA0002844944300000071
S02, importance sampling: sampling distributions from importance
Figure BDA0002844944300000072
To generate a priori particles
Figure BDA0002844944300000073
S03, calculating the weight value: (3) calculating the weight of each particle
Figure BDA0002844944300000074
Normalizing the weights according to the newly obtained measurements
Figure BDA0002844944300000075
S04, resampling: if it is not
Figure BDA0002844944300000076
From discrete distributions
Figure BDA0002844944300000077
In the random sampling of N samples
Figure BDA0002844944300000078
And all weights are set to a constant value
Figure BDA0002844944300000079
Obtaining an approximate posterior probability density function
Figure BDA00028449443000000710
Wherein N isthrFor the set threshold, δ (·) is a dirac δ function;
s05, outputting: obtaining a state estimate
Figure BDA00028449443000000711
That is, the true intracellular endotoxin concentration c during cytotoxicity was estimated from the measurement of the number of cells containing noiseiExtracellular toxin concentration ceAnd the number of cells N.
Referring to FIG. 2, the method for estimating the toxin concentration in wastewater using the present invention is shown in step cePlot of estimated toxin chromium (VI) concentration versus true at 0.62. Wherein the first plot of FIG. 2 is a comparison of the estimated intracellular endotoxin concentration value with the actual value; the second graph is a comparison graph of the estimated value of the concentration of the extracellular toxin and the real value; the third graph is a comparison of the estimated cell number value and the actual value.
Referring to FIG. 3, the method for estimating the toxin concentration in wastewater using the present invention is shown in step ceError plot of estimated toxin chromium (VI) concentration at 0.62. Wherein the first graph of figure 3 is the intracellular endotoxin concentration estimation error; the second graph is the extracellular toxin concentration estimation error; the third plot is the cell number estimation error.
Referring to FIG. 4, the method for estimating the toxin concentration in wastewater using the present invention is shown in step cePlot of estimated toxin mercury (II) chloride concentration versus true at 22.35. Wherein the first plot of FIG. 4 is a comparison of the estimated intracellular endotoxin concentration value with the actual value; the second graph is a comparison graph of the estimated value of the concentration of the extracellular toxin and the real value; the third graph is a comparison of the estimated cell number value and the actual value.
Referring to FIG. 5, the estimation of toxin concentration in wastewater using the present invention is shownCounting method in ceError plot of estimated toxin mercuric (II) chloride concentration at 48.3. Wherein the first graph of figure 5 is the intracellular endotoxin concentration estimation error; the second graph is the extracellular toxin concentration estimation error; the third plot is the cell number estimation error.
According to the graph 2, the real value can be tracked by the sewage toxin concentration estimation method, and the error map of the graph 3 is combined to show that the estimation errors of the sewage toxin concentration estimation method on the cell endotoxin concentration, the cell exotoxin concentration and the cell number are small, so that the accurate estimation on the cell endotoxin concentration, the cell exotoxin concentration and the cell number can be realized.
Similarly, it can be seen from fig. 4 that the estimation and prevention of the concentration of the toxins in wastewater of the present invention can track the true value, and the error map of fig. 5 is combined to determine that the estimation error of the method for estimating the concentration of the toxins in the cells to the concentration of the toxins in the cells, and the number of the cells is very small, so that the accurate estimation of the concentration of the toxins in the cells, and the number of the cells can be realized.
The above-mentioned embodiments are merely preferred embodiments for fully illustrating the present invention, and the scope of the present invention is not limited thereto. The equivalent substitution or change made by the technical personnel in the technical field on the basis of the invention is all within the protection scope of the invention. The protection scope of the invention is subject to the claims.

Claims (9)

1.一种污水毒素浓度估计方法,其特征在于:包括以下步骤,1. a sewage toxin concentration estimation method is characterized in that: comprise the following steps, 将具有一定细胞指数CI的培养细胞暴露在污水样本中,连续采集污水样本内培养细胞的细胞数量测量值;Expose the cultured cells with a certain cell index CI to the sewage sample, and continuously collect the measured value of the cell number of the cultured cells in the sewage sample; 根据细胞对毒素的摄取机理,以及毒素所引起的细胞数量变化创建细胞毒素动态模型;Create a cytotoxin dynamic model according to the uptake mechanism of toxins by cells and the changes in cell number caused by toxins; 使用粒子滤波算法对所述细胞毒素动态模型和所述细胞数量测量值进行处理,估计出污水样本中的毒素浓度;Using a particle filter algorithm to process the cytotoxin dynamic model and the cell number measurements to estimate the toxin concentration in the sewage sample; 创建所述细胞毒素动态模型包括,Creating the cytotoxic dynamic model includes, 创建细胞对毒素摄取过程的动力学模型:Create a kinetic model of the cellular uptake of toxins:
Figure FDA0003170429020000011
Figure FDA0003170429020000011
创建细胞内毒素浓度凋亡机制:Create apoptotic mechanism of intracellular toxin concentration:
Figure FDA0003170429020000012
Figure FDA0003170429020000012
ci表征细胞内毒素浓度;ce表征细胞外毒素浓度;N表征细胞数量,N≈CI;k1,k2,k3,Ki为相关参数;kci表征存在毒素时的相对细胞杀伤率;ks表征无毒素时的相对细胞增殖率。c i is the intracellular toxin concentration; c e is the extracellular toxin concentration; N is the number of cells, N≈CI; k 1 , k 2 , k 3 , K i are the relevant parameters; kc i is the relative cell killing in the presence of toxin rate; k s represents the relative cell proliferation rate in the absence of toxin.
2.如权利要求1所述的污水毒素浓度估计方法,其特征在于:创建所述细胞毒素动态模型包括,2. The method for estimating sewage toxin concentration according to claim 1, wherein: creating the cytotoxin dynamic model comprises: 将具有一定细胞指数CI的培养细胞暴露在不同浓度的毒素中,连续采集有毒素作用下培养细胞的细胞数量测量值;Expose cultured cells with a certain cell index CI to different concentrations of toxins, and continuously collect measurements of the number of cells in cultured cells under the action of toxins; 将有毒素作用下培养细胞的细胞数量测量值代入所述(式1)和(式2)中,并采用非线性回归方法对所述动力学模型和细胞内毒素浓度凋亡机制中的未知参数进行辨识,计算出所述相关参数k1,k2,k3,KiSubstitute the cell number measurements of cells cultured under the action of toxins into the (Equation 1) and (Equation 2), and use a nonlinear regression method for the kinetic model and the unknown parameters of the intracellular toxin concentration in the apoptotic mechanism. Carry out identification, and calculate the relevant parameters k 1 , k 2 , k 3 , and K i . 3.如权利要求2所述的污水毒素浓度估计方法,其特征在于:暴露在不同浓度毒素中的培养细胞的细胞指数CI为1.0~1.2。3 . The method for estimating the concentration of toxins in sewage according to claim 2 , wherein the cell index CI of the cultured cells exposed to different concentrations of toxins is 1.0-1.2. 4 . 4.如权利要求1所述的污水毒素浓度估计方法,其特征在于:暴露在污水样本中的培养细胞的细胞指数CI为1.0~1.2。4 . The method for estimating the concentration of sewage toxins according to claim 1 , wherein the cell index CI of the cultured cells exposed to the sewage sample is 1.0-1.2. 5 . 5.如权利要求1所述的污水毒素浓度估计方法,其特征在于:所述培养细胞的细胞系为NIH 3T3。5. The method for estimating the concentration of sewage toxins according to claim 1, wherein the cell line of the cultured cells is NIH 3T3. 6.如权利要求5所述的污水毒素浓度估计方法,其特征在于:暴露在污水样本中培养细胞、暴露在不同浓度毒素中的培养细胞的初始细胞数均为每孔10000个细胞。6 . The method for estimating the concentration of toxins in sewage according to claim 5 , wherein the initial cell number of the cultured cells exposed to the sewage sample and the cultured cells exposed to different concentrations of toxins is 10,000 cells per well. 7 . 7.如权利要求1~6任一项所述的污水毒素浓度估计方法,其特征在于:污水毒素浓度估计方法还包括,7. The method for estimating sewage toxin concentration according to any one of claims 1 to 6, wherein the method for estimating sewage toxin concentration further comprises: 对所述细胞毒素动态模型进行离散化,得到离散化后的细胞毒素动态模型Discretize the cytotoxin dynamic model to obtain a discretized cytotoxin dynamic model xk=f(xk-1)+wkx k =f(x k-1 )+w k ; yk=g(xk)+vky k =g(x k )+v k ;
Figure FDA0003170429020000031
为状态向量;
Figure FDA0003170429020000032
为测量向量;
Figure FDA0003170429020000033
为关于xk的非线性函数;
Figure FDA0003170429020000034
为观测函数;
Figure FDA0003170429020000035
为过程噪声,且服从分布
Figure FDA0003170429020000036
Figure FDA0003170429020000037
为观测噪声,且服从分布
Figure FDA0003170429020000038
Figure FDA0003170429020000031
is the state vector;
Figure FDA0003170429020000032
is the measurement vector;
Figure FDA0003170429020000033
is a nonlinear function about x k ;
Figure FDA0003170429020000034
is the observation function;
Figure FDA0003170429020000035
is process noise and obeys the distribution
Figure FDA0003170429020000036
Figure FDA0003170429020000037
is observation noise and obeys the distribution
Figure FDA0003170429020000038
使用所述粒子滤波算法结合离散化后的细胞毒素动态模型对污水样本中的细胞数量测量值进行处理。Cell number measurements in sewage samples were processed using the particle filter algorithm in combination with a discretized cytotoxin dynamic model.
8.如权利要求7所述的污水毒素浓度估计方法,其特征在于:使用所述粒子滤波算法处理污水样本中细胞数量测量值包括,定义状态变量x=[ci,ce,N],定义测量向量y=N,基于离散化后的细胞毒素动态模型,并根据含有噪声的细胞数量测量值估计出细胞毒性过程中真实的细胞内毒素浓度ci、细胞外毒素浓度ce和细胞数量N。8 . The method for estimating the concentration of sewage toxins according to claim 7 , wherein: using the particle filter algorithm to process the measured value of the number of cells in the sewage sample comprises: defining a state variable x=[ ci , c e , N], 8 . Define the measurement vector y=N, based on the discretized cytotoxin dynamic model, and estimate the real intracellular toxin concentration c i , extracellular toxin concentration c e and cell number in the process of cytotoxicity according to the measured value of the number of cells containing noise N. 9.如权利要求7所述的污水毒素浓度估计方法,其特征在于:所述离散化的方法包括欧拉方法、龙格-库塔方法、线性多步方法中的一种。9 . The method for estimating sewage toxin concentration according to claim 7 , wherein the discretization method comprises one of Euler method, Runge-Kutta method, and linear multi-step method. 10 .
CN202011506021.2A 2020-12-18 2020-12-18 Method for Estimating Toxin Concentration in Wastewater Active CN112485399B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011506021.2A CN112485399B (en) 2020-12-18 2020-12-18 Method for Estimating Toxin Concentration in Wastewater

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011506021.2A CN112485399B (en) 2020-12-18 2020-12-18 Method for Estimating Toxin Concentration in Wastewater

Publications (2)

Publication Number Publication Date
CN112485399A CN112485399A (en) 2021-03-12
CN112485399B true CN112485399B (en) 2021-09-07

Family

ID=74914827

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011506021.2A Active CN112485399B (en) 2020-12-18 2020-12-18 Method for Estimating Toxin Concentration in Wastewater

Country Status (1)

Country Link
CN (1) CN112485399B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1482887A1 (en) * 1987-06-10 1989-05-30 Украинский Научно-Исследовательский Институт Разведения И Искусственного Осеменения Крупного Рогатого Скота Method of assessing toxic effect of chemicals contained in aqueous medium
JP2001133452A (en) * 1999-11-01 2001-05-18 Nippon Steel Corp Evaluation method for toxic substances
CN1664579A (en) * 2004-12-31 2005-09-07 湖南大学 Application of Spirogyra in Testing the Toxicity of Contaminated Wastewater and Its Test Method
CN104155287A (en) * 2014-09-01 2014-11-19 湖南大学 Method for quantitatively analyzing content of cadmium ion in environmental water body sample
CN108120812A (en) * 2017-12-14 2018-06-05 常州大学 In a kind of dyeing waste water in groups biological toxicity tests and differentiate appraisal procedure
CN109061091A (en) * 2018-09-14 2018-12-21 中国科学院东北地理与农业生态研究所 A method of estimation lake dissolved organic carbon concentration
CN110672804A (en) * 2019-09-30 2020-01-10 华南师范大学 An intelligent detection method for typical organic pollutants in urban river water bodies
CN110878331A (en) * 2019-12-06 2020-03-13 西安外事学院 A kind of microorganism-based nitroaromatic pollutant detection method
CN111122815A (en) * 2019-12-23 2020-05-08 深圳大学 A method for evaluating road rainwater toxicity

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1482887A1 (en) * 1987-06-10 1989-05-30 Украинский Научно-Исследовательский Институт Разведения И Искусственного Осеменения Крупного Рогатого Скота Method of assessing toxic effect of chemicals contained in aqueous medium
JP2001133452A (en) * 1999-11-01 2001-05-18 Nippon Steel Corp Evaluation method for toxic substances
CN1664579A (en) * 2004-12-31 2005-09-07 湖南大学 Application of Spirogyra in Testing the Toxicity of Contaminated Wastewater and Its Test Method
CN104155287A (en) * 2014-09-01 2014-11-19 湖南大学 Method for quantitatively analyzing content of cadmium ion in environmental water body sample
CN108120812A (en) * 2017-12-14 2018-06-05 常州大学 In a kind of dyeing waste water in groups biological toxicity tests and differentiate appraisal procedure
CN109061091A (en) * 2018-09-14 2018-12-21 中国科学院东北地理与农业生态研究所 A method of estimation lake dissolved organic carbon concentration
CN110672804A (en) * 2019-09-30 2020-01-10 华南师范大学 An intelligent detection method for typical organic pollutants in urban river water bodies
CN110878331A (en) * 2019-12-06 2020-03-13 西安外事学院 A kind of microorganism-based nitroaromatic pollutant detection method
CN111122815A (en) * 2019-12-23 2020-05-08 深圳大学 A method for evaluating road rainwater toxicity

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Dynamic Monitoring of Cytotoxicity on Microelectronic Sensors;James Zan Xing et al;《Chem. Res. Toxicol》;20050228;第18卷(第2期);第154-161页 *
实时细胞电子分析技术的应用研究进展;严国俊 等;《中国药学杂志》;20140228;第49卷(第3期);第169-173页 *
模型不确定非线性Markov跳变系统的滤波算法;赵顺毅 等;《控制与决策》;20121130;第27卷(第11期);第1616-1620页 *

Also Published As

Publication number Publication date
CN112485399A (en) 2021-03-12

Similar Documents

Publication Publication Date Title
Ahmadisharaf et al. Calibration and validation of watershed models and advances in uncertainty analysis in TMDL studies
Antanasijević et al. Modelling of dissolved oxygen in the Danube River using artificial neural networks and Monte Carlo Simulation uncertainty analysis
Willacker et al. Reservoirs and water management influence fish mercury concentrations in the western United States and Canada
WO2016179864A1 (en) Fresh water acute standard prediction method based on metal quantitative structure-activity relationship
Miskewitz et al. In-stream dissolved oxygen impacts and sediment oxygen demand resulting from combined sewer overflow discharges
Qian et al. Chlorophyll a as an indicator of microcystin: Short-term forecasting and risk assessment in Lake Erie
Qian et al. Implications of Stein’s paradox for environmental standard compliance assessment
Grayman History of water quality modeling in distribution systems:(013)
Meador et al. Relative importance of water‐quality stressors in predicting fish community responses in midwestern streams
CN112485399B (en) Method for Estimating Toxin Concentration in Wastewater
CN115408646A (en) River pollutant monitoring method and system based on big data
Zhulidov et al. Critical analysis of water quality monitoring in the Russian Federation and former Soviet Union
Malek et al. The effect of chemical parameters on water quality index in machine learning studies: a meta-analysis
Wang et al. No detected toxic concentrations in in situ algal growth inhibition tests—a convenient approach to aquatic ecotoxicology
Albani et al. Bayesian inference and wind field statistical modeling applied to multiple source estimation
CN116976709B (en) Vulnerability assessment method based on river basin ecosystem
Khan et al. Predicting dissolved oxygen concentration in urban watersheds: A comparison of fuzzy number based and Bayesian data-driven approaches
Yağanoğlu et al. Determination of spatial and temporal changes in surface water quality of Filyos River (Turkey) using principal component analysis and cluster analysis
Wang et al. Predicting cytotoxicity of binary pollutants towards a human cell panel in environmental water by experimentation and deep learning methods
Suhartanto et al. Estimation of Rainfall from Climatology Data Using Artificial Neural Networks in Palembang City South Sumatera
Wei et al. Improved 1-D hydrodynamic simulation of reservoir release based on coarse resolution topography data using ISS-SRC-assisted ensemble Kalman filter
Anyanwu et al. Experimental investigations and mathematical modelling of corrosion growth rate on carbon steel under the influence of soil pH and resistivity
CN106682271A (en) Method for determining SWMM water quality scouring model parameter
Douglas Bioaerosol emissions from open windrow composting facilities: Emission characterisation and dispersion modelling improvements
Clerc et al. Exploratory study on a statistical method to analyse time resolved data obtained during nanomaterial exposure measurements

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
TR01 Transfer of patent right

Effective date of registration: 20240514

Address after: Room 401, No. 1118 Jiangbin East Avenue, Dongzhou Street, Fuyang District, Hangzhou City, Zhejiang Province, 310000

Patentee after: Hangzhou Saida Environmental Technology Co.,Ltd.

Country or region after: China

Address before: 1800 No. 214122 Jiangsu city of Wuxi Province Li Lake Avenue

Patentee before: Jiangnan University

Country or region before: China

TR01 Transfer of patent right