CN111204867B - Membrane Bioreactor-MBR Membrane Fouling Intelligent Decision-Making Method - Google Patents
Membrane Bioreactor-MBR Membrane Fouling Intelligent Decision-Making Method Download PDFInfo
- Publication number
- CN111204867B CN111204867B CN202010154288.3A CN202010154288A CN111204867B CN 111204867 B CN111204867 B CN 111204867B CN 202010154288 A CN202010154288 A CN 202010154288A CN 111204867 B CN111204867 B CN 111204867B
- Authority
- CN
- China
- Prior art keywords
- data
- membrane
- hidden layer
- value
- water
- 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
- 239000012528 membrane Substances 0.000 title claims abstract description 113
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000009285 membrane fouling Methods 0.000 title claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 91
- 230000008569 process Effects 0.000 claims abstract description 41
- 230000035699 permeability Effects 0.000 claims abstract description 33
- 239000010865 sewage Substances 0.000 claims abstract description 21
- 238000011156 evaluation Methods 0.000 claims abstract description 10
- 239000013598 vector Substances 0.000 claims description 34
- 238000004519 manufacturing process Methods 0.000 claims description 26
- 239000011159 matrix material Substances 0.000 claims description 26
- 230000006870 function Effects 0.000 claims description 17
- 238000012549 training Methods 0.000 claims description 17
- 238000005273 aeration Methods 0.000 claims description 14
- 239000010802 sludge Substances 0.000 claims description 14
- 238000005259 measurement Methods 0.000 claims description 11
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 10
- 239000001301 oxygen Substances 0.000 claims description 10
- 229910052760 oxygen Inorganic materials 0.000 claims description 10
- 238000004891 communication Methods 0.000 claims description 7
- 210000002569 neuron Anatomy 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 238000011160 research Methods 0.000 claims description 7
- 238000004140 cleaning Methods 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000013210 evaluation model Methods 0.000 claims description 5
- 230000001965 increasing effect Effects 0.000 claims description 5
- 238000012544 monitoring process Methods 0.000 claims description 5
- 230000033116 oxidation-reduction process Effects 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 229910002651 NO3 Inorganic materials 0.000 claims description 2
- NHNBFGGVMKEFGY-UHFFFAOYSA-N Nitrate Chemical compound [O-][N+]([O-])=O NHNBFGGVMKEFGY-UHFFFAOYSA-N 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000005457 optimization Methods 0.000 claims description 2
- 230000002087 whitening effect Effects 0.000 claims description 2
- 230000004927 fusion Effects 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 7
- 241000196324 Embryophyta Species 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000004065 wastewater treatment Methods 0.000 description 2
- 241000170793 Phalaris canariensis Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000008234 soft water Substances 0.000 description 1
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F3/00—Biological treatment of water, waste water, or sewage
- C02F3/30—Aerobic and anaerobic processes
- C02F3/301—Aerobic and anaerobic treatment in the same reactor
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- C—CHEMISTRY; METALLURGY
- C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F2203/00—Apparatus and plants for the biological treatment of water, waste water or sewage
- C02F2203/006—Apparatus and plants for the biological treatment of water, waste water or sewage details of construction, e.g. specially adapted seals, modules, connections
-
- C—CHEMISTRY; METALLURGY
- C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F2209/00—Controlling or monitoring parameters in water treatment
- C02F2209/005—Processes using a programmable logic controller [PLC]
-
- C—CHEMISTRY; METALLURGY
- C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F2209/00—Controlling or monitoring parameters in water treatment
- C02F2209/03—Pressure
-
- C—CHEMISTRY; METALLURGY
- C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F2209/00—Controlling or monitoring parameters in water treatment
- C02F2209/04—Oxidation reduction potential [ORP]
-
- C—CHEMISTRY; METALLURGY
- C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F2209/00—Controlling or monitoring parameters in water treatment
- C02F2209/06—Controlling or monitoring parameters in water treatment pH
-
- C—CHEMISTRY; METALLURGY
- C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F2209/00—Controlling or monitoring parameters in water treatment
- C02F2209/08—Chemical Oxygen Demand [COD]; Biological Oxygen Demand [BOD]
-
- C—CHEMISTRY; METALLURGY
- C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F2209/00—Controlling or monitoring parameters in water treatment
- C02F2209/11—Turbidity
-
- C—CHEMISTRY; METALLURGY
- C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F2209/00—Controlling or monitoring parameters in water treatment
- C02F2209/15—N03-N
-
- C—CHEMISTRY; METALLURGY
- C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F2209/00—Controlling or monitoring parameters in water treatment
- C02F2209/22—O2
-
- C—CHEMISTRY; METALLURGY
- C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F2209/00—Controlling or monitoring parameters in water treatment
- C02F2209/40—Liquid flow rate
-
- C—CHEMISTRY; METALLURGY
- C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F2303/00—Specific treatment goals
- C02F2303/14—Maintenance of water treatment installations
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/152—Water filtration
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Water Supply & Treatment (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Marketing (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Chemical & Material Sciences (AREA)
- Organic Chemistry (AREA)
- Entrepreneurship & Innovation (AREA)
- Biodiversity & Conservation Biology (AREA)
- Microbiology (AREA)
- Hydrology & Water Resources (AREA)
- Operations Research (AREA)
- Environmental & Geological Engineering (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Separation Using Semi-Permeable Membranes (AREA)
- Activated Sludge Processes (AREA)
Abstract
An intelligent decision-making method for membrane bioreactor-MBR membrane pollution belongs to the field of online early warning of water quality parameters in sewage treatment. Firstly, deeply analyzing a membrane pollution mechanism, constructing a membrane pollution prediction model based on a deep belief network, and realizing accurate prediction of membrane water permeability; secondly, fusion of a water permeability predicted value and related parameter variables is utilized, and intelligent early warning of membrane pollution is achieved based on comprehensive evaluation of membrane pollution; secondly, an intelligent decision-making model based on data and knowledge is established based on main factors of membrane pollution, decision support is provided for operators, harm caused by membrane pollution is reduced, the occurrence rate of membrane pollution is reduced, the safety of MBR sewage treatment process is guaranteed, and the efficient and stable operation of MBR sewage treatment plants is promoted.
Description
Technical Field
The invention belongs to the field of online early warning of water quality parameters in sewage treatment, and provides an intelligent early warning system for membrane bioreactor-MBR membrane pollution for the first time. On the basis of real operation data of an MBR membrane sewage treatment process, extracting characteristic variables of MBR membrane water permeability by a characteristic analysis method, and establishing a soft measurement model by using a deep belief network to predict the membrane water permeability which is difficult to directly measure in the MBR membrane sewage treatment process; establishing a membrane pollution comprehensive evaluation model by utilizing the water permeability prediction value and combining other collectable process variables of a water plant, determining main factors of membrane pollution and realizing intelligent early warning of the membrane pollution; based on data and knowledge drive, an intelligent decision model is established, accurate identification of membrane pollution is realized, the incidence rate of membrane pollution is reduced, and the sewage treatment performance of the membrane is improved.
Background
MBR is a novel sewage treatment process combining membrane separation technology and biological treatment technology, has the obvious advantages of good solid-liquid separation effect, low sludge load, small occupied area and the like, and has great development potential. However, membrane fouling is a bottleneck problem in MBR wastewater treatment processes, which can lead to reduced effluent quality and increased operating costs, and even to a breakdown of the wastewater treatment process. Currently, the main approach to the membrane fouling problem is to periodically clean and replace the membranes. In the actual sewage treatment process, the cleaning and the replacement of the membrane do not have a strict, objective and quantitative standard and are mainly carried out according to human experience. However, frequent cleaning can lead to rupture and corrosion of membrane filaments, which can reduce membrane life, increase production energy consumption and operation cost, and greatly restrict popularization and application of MBR. Therefore, the method can accurately identify the membrane pollution, reduce the incidence rate of the membrane pollution, and become a key point for ensuring the stable operation of the MBR and popularizing the MBR technology. However, MBR sewage treatment process is complex and difficult to directly model, and the monitoring of pollution condition is a difficult problem in the current control field. Currently, an effective early warning decision system is not available in a membrane sewage treatment plant which is built and put into operation so as to realize intelligent early warning decision in the membrane sewage treatment process. Therefore, the research of new decision-making technology for solving the membrane pollution problem in the sewage treatment process becomes an important subject of research in the field of sewage control and has important practical significance.
The invention relates to an intelligent decision-making method for membrane bioreactor-MBR membrane pollution, which utilizes a characteristic analysis method to extract characteristic variables and establishes a soft measurement model of membrane water permeability based on a deep belief network to realize accurate prediction of the water permeability in the process of membrane sewage treatment, utilizes a water permeability prediction value to establish a membrane pollution comprehensive evaluation model in combination with other collectable process variables of a water plant, but an intelligent decision-making system for membrane pollution at home and abroad does not form a complete theoretical system, builds an intelligent decision-making method for MBR membrane pollution including soft and hard platforms based on an intelligent method, and has high development and application values in the aspects of filling technical blanks at home and abroad, integrating a sewage treatment industrial chain and the like.
Disclosure of Invention
The membrane bioreactor-MBR membrane pollution intelligent decision method comprises the steps of operation process data acquisition, operation process data preprocessing, membrane pollution intelligent prediction and membrane pollution intelligent decision, and specifically comprises the following steps:
(1) acquiring data of an operation process: the MBR membrane processing system is used as a research object, and the data of the operation process is collected through a collection instrument arranged on a process field, and the MBR membrane processing system comprises the following steps: the method comprises the following steps of (1) realizing data acquisition by using inlet water chemical oxygen demand, inlet water pH value, inlet water biological oxygen demand, outlet water chemical oxygen demand, outlet water pH value, outlet water biological oxygen demand, anaerobic zone oxidation-reduction potential, anoxic zone oxidation-reduction potential, aerobic zone nitrate, aerobic zone dissolved oxygen, water production pressure, water production turbidity, water production flow, sludge concentration and aeration quantity; the data collected by the instrument is transmitted to the programmable logic controller through a communication protocol, the programmable logic controller transmits the operation process data to the upper computer through the communication protocol, and the data in the upper computer is transmitted to the data processing server through the local area network;
(2) preprocessing operation process data: the operating data of the membrane pool is taken as a research object, a characteristic analysis model is established by utilizing a partial least square method, and 5 principal component variables are obtained, wherein the principal component variables are respectively as follows: the water production pressure, the water production turbidity, the water production flow, the sludge concentration and the aeration rate, and 5 main component variables are used as input variables of the membrane pollution intelligent prediction model block; the water permeability is used as an output variable of the membrane pollution intelligent prediction module;
(3) intelligent prediction of membrane pollution: the water permeability prediction is realized by intelligent prediction of membrane pollution, wherein: the water permeability is obtained by predicting a deep belief network-DBN, the DBN consists of 1 input layer, 2 hidden layers and 1 output layer, the number of neurons in the input layers is 5, the number of neurons in each hidden layer is M, M is a positive integer which is more than 2 and less than 30, and the number of neurons in the output layers is 1, namely the connection mode is 5-M-M-1; n groups of data are used as training samples of the soft measurement model; DBN input at time t is x (t) ═ x1(t),…,x5(t)],x1(t) value x representing water pressure at time t2(t) value x representing turbidity of water produced at time t3(t) value x representing the water production flow at time t4(t) value x representing sludge concentration at time t5(t) value representing aeration amount at time t, soft water permeability prediction based on DBNThe calculation mode of the measurement model is as follows in sequence:
inputting a layer: h is0(t)=x(t)
(1)
wherein h is0(t) represents the output vector of the input layer at time t, h1(t) represents the output vector of the first hidden layer at time t, h2(t) represents the output vector of the second hidden layer at time t, w0,1(t) represents the weight vector between the input layer and the first hidden layer at time t, w1,2(t) represents the weight vector between the first hidden layer and the second hidden layer at time t, w2,3(t) represents the weight vector between the second hidden layer and the output layer at time t, b1(t) denotes the bias vector of the first hidden layer at time t, b2(t) represents the bias vector of the second hidden layer at time t, and y (t) represents the actual output of the DBN at time t;
the DBN training is divided into two processes: unsupervised pre-training and supervised weight fine tuning; setting the iteration number of each layer of pre-training as 100, the iteration number of a back propagation algorithm as 1000, the expected error as 0.01, and the initial weight and bias as 0.01; the specific training steps are as follows:
unsupervised pre-training: the update rule of the parameters obtained by the contrast divergence algorithm is as follows:
wherein, w0,1(t +1) is the input layer and the first at time t +1Weight vectors between hidden layers, w1,2(t +1) is the weight vector between the first hidden layer and the second hidden layer at time t +1, b1(t +1) is the bias vector for the first hidden layer at time t +1, b2(t +1) is the bias vector for the second hidden layer at time t +1, EdataIs a expectation of training data, EmodelIs an expectation of model definition, μw1E (0,0.02) is the learning rate of the weights of the input layer and the first hidden layer, muw2E (0,0.02) is the learning rate of the weights of the first hidden layer and the second hidden layer, mub1E (0,0.01) is the learning rate of the first hidden layer, μb2E (0,0.01) is the learning rate of the second hidden layer;
secondly, adjusting the weight value by a back propagation algorithm: training layer by layer to obtain the initial value of the DBN parameter, then fine-tuning the weight value through a back propagation algorithm to obtain a better model effect, and adjusting the weight value by adopting an error back propagation method as follows:
wherein, w2,3(t +1) is the weight vector between the second hidden layer and the output layer at time t +1, yd(t) is the desired output of DBN at time t, ηoutE (0,3) is the weight learning rate between the second hidden layer and the output layer, η2E (0,3) is the weight learning rate between the first hidden layer and the second hidden layer, η1E (0,3) is a weight learning rate between the input layer and the first hidden layer;
(4) intelligent decision making of membrane pollution: the membrane pollution intelligent decision is to establish a membrane pollution comprehensive evaluation model by utilizing a water permeability prediction value and combining other membrane pollution related variables, and provide a decision aiming at membrane pollution, and the process is as follows:
1) determining an evaluation index of membrane fouling, u ═ u { [ u ]1,u2,u3,u4,u5,u6},u1Is the value of the water production flow, u2Is the value of the water production pressure u3Is the value of turbidity of the produced water, u4Is the value of the sludge concentration u5Is aeratedValue sum u6Is a predicted value of the water permeability;
2) establishing monitoring statistics, wherein the process variable is u epsilon R6×nWherein 6 is a variable dimension, and n is the number of samples; and (3) carrying out independent component and principal component decomposition on the membrane pollution evaluation matrix u:
u=As+PtT+fT (7)
wherein A ∈ R6×rFor the mixing matrix, s ∈ Rr×nFor independent matrices, t ∈ Rn×kAs principal component scoring matrix, load matrix P is formed by R6×k,f∈Rn×6The final residual matrix is obtained, r represents the number of independent elements, and k is the number of selected principal elements; to estimate A and s, a unmixing matrix W is required to obtain a reconstructed signal Is an estimated value of a source signal, and performs whitening processing on u to obtain a score vector:
Z=Λ-1/2HTu=Bs (8)
wherein Λ is a diagonal matrix containing all eigenvalues, H is a corresponding eigenvector matrix and the unit is orthogonal, and B is a unit orthogonal matrix; the unmixing matrix is:
W=BTΛ-1/2HT (9)
decomposing the original data by independent components and principal components to obtain independent component matrix and principal component scoring matrix, and establishing I2And T2Statistics, namely establishing SPE statistics for the residual information; independent component vector is sl=[s1l, s2l,…,srl]T∈Rr×1And a pivot score vector tl=[tl1,tl2,…,tlk]T∈Rk×1Residual vector fl=[fl1, fl2,…,fl6]T∈ R 6×11,2, …, n (n is the number of samples); establishment of I2、T2And the SPE statistics are as follows:
whereinIs the control limit for the SPE statistic,λj(j ═ 1, 2., k) is the eigenvalue corresponding to each principal component covariance matrix, k is the number of principal components, Fk.(n-k),αThe upper limit value of F distribution with confidence coefficient alpha and degree of freedom k and (n-k); c. CαCorresponding to the lower limit value of the normal distribution under the confidence coefficient alpha, determining I by adopting a nuclear density estimation method2Control limit of statistic, determining T based on F distribution2A control limit for the statistic;
3) intelligent decision making: obtaining the membrane pollution evaluation index u (t) at the time t on line, and monitoring the statistic I at the time t2、T2And SPE is beyond the control limit, and the occurrence of membrane pollution is determined; in order to determine the main factors of membrane pollution, measures are taken aiming at the determination factors, and a multi-classifier based on a kernel function is established for distinguishing the membrane pollution factors, and the method mainly comprises the following steps: flow of produced water reaches peak value of 460m3H, the pressure of produced water is lower than 20kp, and the aeration quantity is lower than 2400m3(ii) the sludge concentration is higher than 13000mg/l, (ii) the water production pressure is lower than 20kp, the water permeability is lower than 30LMH/bar, (iii) the water production pressure is lower than 20kp, the water permeability is lower than 60LMH/bar, and the turbidity of the produced water is higher than 5 NUT; combining 6 binary classifiers into a multi-classifier, wherein for a first binary classifier, the data labels not belonging to the first class are-1, the data labels belonging to the first class are +1, for a second binary classifier, the data labels not belonging to the second class are-1, the data labels belonging to the second class are +1, for a third binary classifier, the data labels not belonging to the third class are-1, the data labels belonging to the third class are +1, and for a fourth binary classifier, the data labels not belonging to the third class are +1The classifier labels-1 for data not belonging to the fourth class, +1 for data belonging to the fourth class, -1 for data not belonging to the fifth class, +1 for data belonging to the fifth class, -1 for data not belonging to the sixth class, and +1 for data belonging to the sixth class, and the q-th classifier (q ═ 1,2, …,6) optimizes the objective function as:
wherein C ∈ (0,5) is a penalty coefficient, wqIs the connection weight of the q-th classifier, T is the transpose symbol, ξp qRelaxation variables for the p-th data sample of the q-th classifier, bqFor the bias of the q-th classifier,is a Gaussian kernel function, σ ∈ (0,1) is the width of the kernel function, up qP-th data sample, u, of the q-th classifier, respectivelyv qIs the v-th data sample of the q-th classifier (v ═ 1,2, …, n), ypIs the category of the p sample; introducing a Lagrange multiplier, and converting the optimization problem into:
updating parameters:
wherein, yvIs the class of the v-th sample, αpThe glanz multiplier for the p-th data sample,αvfor the evaluation index u, the lagrangian multiplier for the v-th data sample is as follows:
f(u)=sgn(αpypK(u,up q)+bq) (15)
wherein sgn represents a sign function, the input is positive, the output is 1, and the output is-1 when the input is negative; the output value can be obtained by inputting the sample to be detected into the decision function, and the category corresponding to the function with the maximum value is the category to which the sample to be detected belongs; the 6-type corresponding operation suggestions comprise: running for 4 hours in the state is not suitable, and the water yield is reduced to 260m3Lower than h, reducing the water yield of the membrane pool to 260m3Per hour, controlling transmembrane pressure difference to be less than 40kPa, increasing aeration to 3000m3Fourthly, controlling the sludge concentration of the membrane tank to be 8000- & lt12000 mg/L, & ltv & gt, adjusting the operation parameters, and reducing the water yield to 260m3H, or increasing aeration amount to 5000m3Over per hour, performing online physical cleaning and online chemical cleaning within 24 hours; and applying the trained decision model to a sewage treatment process, carrying out characteristic matching on fault data by the model, outputting fault categories, carrying out operation suggestions corresponding to the fault categories, and providing decision support for the production process.
Drawings
FIG. 1 is a diagram of a deep belief network architecture;
FIG. 2 is a diagram of 10-step prediction results of a water outlet water permeability soft measurement model, wherein black line circles are actual water permeability values, and black line points are predicted values of a deep belief network soft measurement model;
FIG. 3 is a diagram showing the MBR membrane contamination early warning result, wherein (a) is a non-Gaussian statistic I2Changing result, wherein the black line points are statistics change curve, the black line is control limit, and the graph (b) is Gaussian statistic T2Changing results, wherein black line band points are a statistic change curve, black lines are control limits, and a graph (c) is a non-Gaussian statistic SPE changing result, wherein the black line band points are the statistic change curve, and the black lines are the control limits;
FIG. 4 is a plot of MBR membrane fouling decision results, wherein black circles are the decision model outputs and black stars are the actual decision suggestion categories;
FIG. 5 is a data flow direction indicating diagram of the MBR membrane pollution intelligent early warning system;
Detailed Description
(1) Specific implementation of membrane pollution intelligent decision system design and software and hardware function integration
The hardware platform environment built in an actual sewage treatment plant is shown in fig. 2. Collecting operation process data through a collection instrument installed in a process site; data collected by the instrument is transmitted to the PLC through a Modbus communication protocol, the PLC transmits operation process data to the upper computer through an RS232 communication protocol, and the data in the upper computer is transmitted to the data processing server through the local area network; issuing operation process data to a water plant work manager through a Web Server in a Browser/Server mode, and displaying the prediction of the water permeability and the early warning result of membrane pollution in a Client/Server mode; the main functions of the developed MBR membrane pollution intelligent early warning system can be realized by firstly inquiring the operation parameters of the membrane pool, secondly predicting the water permeability on line, thirdly early warning the membrane pollution and fourthly carrying out intelligent decision on the membrane pollution.
The invention adopts component technology in software industry to package the membrane pollution data preprocessing module, the membrane pollution intelligent prediction module, the membrane pollution intelligent early warning module and the membrane pollution intelligent decision-making module into functional modules, thereby enhancing the reusability of the model and making up the blank of popularization of MBR membrane pollution intelligent decision-making technology to human-computer interaction interface in actual system operation at home and abroad; the NET platform is adopted for software development, so that an ActiveX control is conveniently created, and the usable environment range of software is expanded; a field bus technology is adopted to establish a full-flow system communication network to realize information transmission among modules; meanwhile, the MBR membrane pollution intelligent decision-making system provided by the invention realizes the connection of a central control room and each on-site data acquisition point, forms a centralized management early warning system, is easy to expand, has independent functions of each part, can add software and hardware modules according to actual prediction requirements and is fused with other systems, can realize the stability and reliability of the system and ensures the early warning precision of membrane pollution.
(2) Specific implementation of film pollution intelligent decision method research
The invention obtains an MBR membrane pollution intelligent decision-making method; the method is characterized in that characteristic variables of MBR membrane water permeability are obtained through characteristic analysis, a soft measurement model of MBR membrane water permeability is established by utilizing a deep belief network, intelligent detection of MBR membrane water permeability is realized, a membrane pollution comprehensive evaluation model is established by combining a predicted value of the membrane water permeability with other collectable process variables of a water plant, identification of membrane pollution and identification of main factors are realized, the intelligent level of the sewage treatment plant is improved, and normal operation of a sewage treatment process is guaranteed.
Firstly, input variables are collected through an online detection instrument arranged on a process field, the variables to be collected comprise 5 types, and parameter information and collection positions are shown in a table 1.
TABLE 1 Process variable types collected
Secondly, establishing a soft measurement model by using a deep belief network, and realizing 510-step prediction of the membrane water permeability based on a multi-step prediction strategy. And training and testing the deep belief network by adopting the data acquired in real time. 80 sets of data were selected for testing. The data collected are shown in table 2.
And thirdly, performing comprehensive evaluation on membrane pollution by using the predicted value of the water permeability and current values of other related collected variables (water production flow, water production pressure, aeration amount and sludge concentration) so as to judge whether membrane pollution occurs.
And fourthly, making a membrane pollution operation decision according to the relevant variables of the membrane pollution.
TABLE 2 Soft measurement model test data
Claims (1)
1. The membrane bioreactor-MBR membrane pollution intelligent decision-making method is characterized by comprising the following steps:
(1) acquiring data of an operation process: the MBR membrane treatment system is used as a research object, and the data of the operation process is collected through a collection instrument arranged on a process field, and the MBR membrane treatment system comprises the following steps: the method comprises the following steps of (1) realizing data acquisition by using inlet water chemical oxygen demand, inlet water pH value, inlet water biological oxygen demand, outlet water chemical oxygen demand, outlet water pH value, outlet water biological oxygen demand, anaerobic zone oxidation-reduction potential, anoxic zone oxidation-reduction potential, aerobic zone nitrate, aerobic zone dissolved oxygen, water production pressure, water production turbidity, water production flow, sludge concentration and aeration quantity; the data collected by the instrument is transmitted to the programmable logic controller through a communication protocol, the programmable logic controller transmits the operation process data to the upper computer through the communication protocol, and the data in the upper computer is transmitted to the data processing server through the local area network;
(2) preprocessing operation process data: the operating data of the membrane pool is taken as a research object, a characteristic analysis model is established by utilizing a partial least square method, and 5 principal component variables are obtained, wherein the principal component variables are respectively as follows: the water production pressure, the water production turbidity, the water production flow, the sludge concentration and the aeration rate, and 5 main component variables are used as input variables of the membrane pollution intelligent prediction model block; the water permeability is used as an output variable of the membrane pollution intelligent prediction module;
(3) intelligent prediction of membrane pollution: the water permeability prediction is realized by intelligent prediction of membrane pollution, wherein: the water permeability is obtained by predicting a deep belief network-DBN, the DBN consists of 1 input layer, 2 hidden layers and 1 output layer, the number of neurons of the input layer is 5, the number of neurons of each hidden layer is M, M is a positive integer larger than 2 and smaller than 30, and the number of neurons of the output layer is less than or equal to the number of neurons of each hidden layerThe number of warp elements is 1, namely the connection mode is 5-M-M-1; n groups of data are used as training samples of the soft measurement model; DBN input at time t is x (t) ═ x1(t),…,x5(t)],x1(t) value x representing water pressure at time t2(t) value x representing turbidity of water produced at time t3(t) value x representing the water production flow at time t4(t) value x representing sludge concentration at time t5(t) represents the value of aeration at the time t, and the calculation modes of the soft measurement model for predicting the water permeability based on the DBN sequentially comprise the following steps:
inputting a layer: h is0(t)=x(t) (1)
wherein h is0(t) represents the output vector of the input layer at time t, h1(t) represents the output vector of the first hidden layer at time t, h2(t) represents the output vector of the second hidden layer at time t, w0,1(t) represents the weight vector between the input layer and the first hidden layer at time t, w1,2(t) represents the weight vector between the first hidden layer and the second hidden layer at time t, w2,3(t) represents the weight vector between the second hidden layer and the output layer at time t, b1(t) denotes the bias vector of the first hidden layer at time t, b2(t) represents the bias vector of the second hidden layer at time t, and y (t) represents the actual output of the DBN at time t;
the DBN training is divided into two processes: unsupervised pre-training and supervised weight fine tuning; setting the iteration number of each layer of pre-training as 100, the iteration number of a back propagation algorithm as 1000, the expected error as 0.01, and the initial weight and bias as 0.01; the specific training steps are as follows:
unsupervised pre-training: the update rule of the parameters obtained by the contrast divergence algorithm is as follows:
wherein, w0,1(t +1) is the weight vector between the input layer and the first hidden layer at time t +1, w1,2(t +1) is the weight vector between the first hidden layer and the second hidden layer at time t +1, b1(t +1) is the bias vector for the first hidden layer at time t +1, b2(t +1) is the bias vector for the second hidden layer at time t +1, EdataIs a expectation of training data, EmodelIs an expectation of model definition, μw1E (0,0.02) is the learning rate of the weights of the input layer and the first hidden layer, muw2E (0,0.02) is the learning rate of the weights of the first hidden layer and the second hidden layer, mub1E (0,0.01) is the learning rate of the first hidden layer, μb2E (0,0.01) is the learning rate of the second hidden layer;
secondly, adjusting the weight value by a back propagation algorithm: training layer by layer to obtain the initial value of the DBN parameter, then finely adjusting the weight value through a back propagation algorithm to obtain a better model effect, and adjusting the weight value by adopting an error back propagation method to be as follows:
wherein, w2,3(t +1) is the weight vector between the second hidden layer and the output layer at time t +1, yd(t) is the desired output of DBN at time t, ηoutE (0,3) is the weight learning rate between the second hidden layer and the output layer, η2E (0,3) is the weight learning rate between the first hidden layer and the second hidden layer, η1E (0,3) is a weight learning rate between the input layer and the first hidden layer;
(4) intelligent decision making of membrane pollution: the membrane pollution intelligent decision is to establish a membrane pollution comprehensive evaluation model by utilizing a water permeability prediction value and combining other membrane pollution related variables, and provide a decision aiming at membrane pollution, and the process is as follows:
1) determining an evaluation index of membrane fouling, u ═ u { [ u ]1,u2,u3,u4,u5,u6},u1Is the value of the water production flow, u2Is the value of the water production pressure u3Is the value of turbidity of the produced water, u4Is the value of the sludge concentration u5Is the value of aeration amount and u6Is a predicted value of the water permeability;
2) establishing monitoring statistics, wherein the process variable is u epsilon R6×nWherein 6 is a variable dimension, and n is the number of samples; and (3) carrying out independent component and principal component decomposition on the membrane pollution evaluation matrix u:
u=As+PtT+fT (7)
wherein A ∈ R6×rFor the mixing matrix, s ∈ Rr×nFor independent matrices, t ∈ Rn×kFor the principal component scoring matrix, the load matrix P is the same as R6 ×k,f∈Rn×6The final residual matrix is r, the number of independent elements is represented by r, and k is the number of selected principal elements; to estimate A and s, a unmixing matrix W is required to obtain a reconstructed signal Is an estimated value of a source signal, and performs whitening processing on u to obtain a score vector:
Z=Λ-1/2HTu=Bs (8)
wherein Λ is a diagonal matrix containing all eigenvalues, H is a corresponding eigenvector matrix and the unit is orthogonal, and B is a unit orthogonal matrix; the unmixing matrix is:
W=BTΛ-1/2HT (9)
raw data is passed through independent components and principal componentsAfter decomposition, obtaining independent component matrix and principal component scoring matrix, and establishing I2And T2Statistics, namely establishing SPE statistics for the residual information; independent component vector is sl=[s1l,s2l,…,srl]T∈Rr×1And a pivot score vector tl=[tl1,tl2,…,tlk]T∈Rk×1Residual vector fl=[fl1,fl2,…,fl6]T∈R6×11,2, …, n (n is the number of samples); establishment of I2、T2And the SPE statistics are as follows:
whereinIs the control limit for the SPE statistic,λjeigenvalues corresponding to the covariance matrix of each principal component, where j is 1, 2. k is the number of pivot elements, Fk.(n-k),αIs the upper limit value of F distribution with confidence coefficient alpha and freedom degrees k and (n-k); c. CαCorresponding to the lower limit value of normal distribution under the confidence coefficient alpha, determining I by adopting a nuclear density estimation method2Control limit of statistic, determining T based on F distribution2A control limit for the statistic;
3) intelligent decision making: obtaining the membrane pollution evaluation index u (t) at the time t on line, and monitoring the statistic I at the time t2、T2And SPE is beyond the control limit, and the occurrence of membrane pollution is determined; in order to determine the main factors of membrane pollution, measures are taken aiming at the determination factors, and a multi-classifier based on a kernel function is established for distinguishing the membrane pollution factors, and the method mainly comprises the following steps: flow of produced water reaches peak value of 460m3H, the water production pressure is lower than 20kPa, and the aeration quantity is lower than 2400m3High sludge concentrationAt 13000mg/l, water pressure lower than 20kPa and water permeability lower than 30LMH/bar, etc.; combining 6 binary classifiers into a multi-classifier, wherein for a first binary classifier, data not belonging to the first class is labeled with-1, data belonging to the first class is labeled with +1, for a second binary classifier, data not belonging to the second class is labeled with-1, data belonging to the second class is labeled with +1, for a third binary classifier, data not belonging to the third class is labeled with-1, data belonging to the third class is labeled with +1, for a fourth binary classifier, data not belonging to the fourth class is labeled with-1, data belonging to the fourth class is labeled with +1, for a fifth binary classifier, data not belonging to the fifth class is labeled with-1, data belonging to the fifth class is labeled with +1, for a sixth binary classifier, data not belonging to the sixth class is labeled with-1, the data belonging to the sixth class is labeled +1, and the q-th classifier (q 1,2, …,6) optimizes the objective function as:
wherein C ∈ (0,5) is a penalty coefficient, wqIs the connection weight of the q-th classifier, T is the transposed symbol, xip qRelaxation variables for the p-th data sample of the q-th classifier, bqFor the bias of the q-th classifier,is a Gaussian kernel function, σ ∈ (0,1) is the width of the kernel function, up qP-th data sample, u, of the q-th classifier, respectivelyv qIs the v-th data sample of the q-th classifier (v ═ 1,2, …, n), ypIs the category of the p sample; introducing a Lagrange multiplier, and converting the optimization problem into:
updating parameters:
wherein, yvIs the class of the v-th sample, αpIs the Grenarian multiplier, alpha, of the p-th data samplevFor the evaluation index u, the lagrangian multiplier for the v-th data sample is as follows:
f(u)=sgn(αpypK(u,up q)+bq) (15)
wherein sgn represents a sign function, the input is positive, the output is 1, and the output is-1 when the input is negative; the output value can be obtained by inputting the sample to be detected into the decision function, and the category corresponding to the function with the maximum value is the category to which the sample to be detected belongs; the 6-type corresponding operation suggestions comprise: running in this state should not exceed 4h, and the water yield is reduced to 260m3Lower than h, reducing the water yield of the membrane pool to 260m3Per hour, controlling transmembrane pressure difference to be less than 40kPa, increasing aeration to 3000m3Fourthly, controlling the sludge concentration of the membrane tank to be 8000- & lt12000 mg/L, & ltv & gt, adjusting the operation parameters, and reducing the water yield to 260m3H, or increasing aeration amount to 5000m3Over per hour, performing online physical cleaning and online chemical cleaning within 24 hours; and applying the trained decision model to a sewage treatment process, performing feature matching on fault data by using the model, outputting fault categories, performing operation suggestions corresponding to the fault categories, and providing decision support for the production process.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910549549 | 2019-06-24 | ||
CN2019105495499 | 2019-06-24 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111204867A CN111204867A (en) | 2020-05-29 |
CN111204867B true CN111204867B (en) | 2021-12-10 |
Family
ID=70781469
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010154288.3A Active CN111204867B (en) | 2019-06-24 | 2020-03-07 | Membrane Bioreactor-MBR Membrane Fouling Intelligent Decision-Making Method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111204867B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113283481A (en) * | 2021-05-14 | 2021-08-20 | 群智未来人工智能科技研究院(无锡)有限公司 | Intelligent membrane pollution decision-making method based on knowledge type-two fuzzy |
CN116738296B (en) * | 2023-08-14 | 2024-04-02 | 大有期货有限公司 | Comprehensive intelligent monitoring system for machine room conditions |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5795268B2 (en) * | 2009-02-27 | 2015-10-14 | オーソ−クリニカル・ダイアグノスティックス・インコーポレイテッドOrtho−Clinical Diagnostics, Inc. | Method for detecting an impending analytical failure of a networked diagnostic clinical analyzer |
CN101944275B (en) * | 2010-08-26 | 2012-07-04 | 天津市环境保护科学研究院 | Membrane pollution diagnosis and early warning decision making system of hollow fiber device |
US9177550B2 (en) * | 2013-03-06 | 2015-11-03 | Microsoft Technology Licensing, Llc | Conservatively adapting a deep neural network in a recognition system |
CN104478172B (en) * | 2014-12-17 | 2017-02-22 | 中国科学院生态环境研究中心 | Method and device for treating organic wastewater with two-phase anaerobic tubular membrane bioreactor |
US9984201B2 (en) * | 2015-01-18 | 2018-05-29 | Youhealth Biotech, Limited | Method and system for determining cancer status |
CN105373110B (en) * | 2015-12-16 | 2018-06-05 | 浙江中烟工业有限责任公司 | Cigarette ultrahigh speed film wrapping machine multi-state production process is monitored on-line and method for diagnosing faults |
CN107909084B (en) * | 2017-11-15 | 2021-07-13 | 电子科技大学 | A convolution-linear regression network-based prediction method for haze concentration |
CN109133351A (en) * | 2018-08-29 | 2019-01-04 | 北京工业大学 | Membrane bioreactor-MBR fouling membrane intelligent early-warning method |
CN109325683A (en) * | 2018-09-19 | 2019-02-12 | 北京城市排水集团有限责任公司 | A kind of film properties integrated evaluating method based on analytic hierarchy process (AHP) |
CN109062196B (en) * | 2018-10-31 | 2020-12-15 | 东北大学 | A blast furnace process monitoring and fault diagnosis method integrating PCA-ICA |
CN109473182B (en) * | 2018-11-12 | 2021-10-19 | 北京北排科技有限公司 | MBR (Membrane biological reactor) membrane water permeability intelligent detection method based on deep belief network |
US11551136B2 (en) * | 2018-11-14 | 2023-01-10 | Tencent America LLC | N-best softmax smoothing for minimum bayes risk training of attention based sequence-to-sequence models |
CN112101402B (en) * | 2020-07-22 | 2024-04-09 | 北京工业大学 | Membrane pollution early warning method based on knowledge fuzzy learning |
-
2020
- 2020-03-07 CN CN202010154288.3A patent/CN111204867B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN111204867A (en) | 2020-05-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Nourani et al. | Artificial intelligence based ensemble modeling of wastewater treatment plant using jittered data | |
CN117170221B (en) | An artificial intelligence control system for sewage treatment | |
CN109133351A (en) | Membrane bioreactor-MBR fouling membrane intelligent early-warning method | |
CN109492265B (en) | Wastewater effluent index prediction method based on dynamic nonlinear PLS soft measurement method | |
CN102262147A (en) | Soft measurement method and system for effluent chemical oxygen demand (COD) of waste water treatment system | |
CN202075272U (en) | Soft measurement system applied to sewage treatment | |
CN112101402B (en) | Membrane pollution early warning method based on knowledge fuzzy learning | |
CN110320335B (en) | Multi-element robust soft measurement method for effluent quality index of sewage treatment | |
CN106682316A (en) | Real-time effluent total-phosphorus monitoring system based on peak radial basis function neural network | |
CN105574326A (en) | Self-organizing fuzzy neural network-based soft measurement method for effluent ammonia-nitrogen concentration | |
CN108375534A (en) | MBR fouling membrane intelligent early-warning methods | |
US20220267169A1 (en) | Intelligent Detection System of Effluent Total Nitrogen based on Fuzzy Transfer Learning Algorithm | |
CN112417765B (en) | Sewage treatment process fault detection method based on improved teacher-student network model | |
Wu et al. | Data-driven intelligent warning method for membrane fouling | |
CN111204867B (en) | Membrane Bioreactor-MBR Membrane Fouling Intelligent Decision-Making Method | |
WO2021114320A1 (en) | Wastewater treatment process fault monitoring method using oica-rnn fusion model | |
CN108595892A (en) | Soft-measuring modeling method based on time difference model | |
CN106706491B (en) | Intelligent detection method for membrane bioreactor-MBR water permeability | |
CN117192063A (en) | Water quality prediction method and system based on coupled Kalman filtering data assimilation | |
Wang et al. | Prediction of effluent ammonia nitrogen in wastewater treatment plant based on self-organizing hybrid neural network | |
Ren et al. | Multivariable identification of membrane fouling based on compacted cascade neural network | |
CN201330211Y (en) | Working parameter self-optimizing simulation system for sewage treatment plant | |
Zhang et al. | Effluent Quality Prediction of Wastewater Treatment System Based on Small-world ANN. | |
CN103675005B (en) | The industrial melt index soft measurement instrument of optimum FUZZY NETWORK and method | |
Wang et al. | Neurodynamics-driven prediction model for state evolution of coastal water-quality |
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 |