[go: up one dir, main page]

CN111204867B - Membrane Bioreactor-MBR Membrane Fouling Intelligent Decision-Making Method - Google Patents

Membrane Bioreactor-MBR Membrane Fouling Intelligent Decision-Making Method Download PDF

Info

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
Application number
CN202010154288.3A
Other languages
Chinese (zh)
Other versions
CN111204867A (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.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
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 Beijing University of Technology filed Critical Beijing University of Technology
Publication of CN111204867A publication Critical patent/CN111204867A/en
Application granted granted Critical
Publication of CN111204867B publication Critical patent/CN111204867B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F3/00Biological treatment of water, waste water, or sewage
    • C02F3/30Aerobic and anaerobic processes
    • C02F3/301Aerobic and anaerobic treatment in the same reactor
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2203/00Apparatus and plants for the biological treatment of water, waste water or sewage
    • C02F2203/006Apparatus and plants for the biological treatment of water, waste water or sewage details of construction, e.g. specially adapted seals, modules, connections
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2209/00Controlling or monitoring parameters in water treatment
    • C02F2209/005Processes using a programmable logic controller [PLC]
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2209/00Controlling or monitoring parameters in water treatment
    • C02F2209/03Pressure
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2209/00Controlling or monitoring parameters in water treatment
    • C02F2209/04Oxidation reduction potential [ORP]
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2209/00Controlling or monitoring parameters in water treatment
    • C02F2209/06Controlling or monitoring parameters in water treatment pH
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2209/00Controlling or monitoring parameters in water treatment
    • C02F2209/08Chemical Oxygen Demand [COD]; Biological Oxygen Demand [BOD]
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2209/00Controlling or monitoring parameters in water treatment
    • C02F2209/11Turbidity
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2209/00Controlling or monitoring parameters in water treatment
    • C02F2209/15N03-N
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2209/00Controlling or monitoring parameters in water treatment
    • C02F2209/22O2
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2209/00Controlling or monitoring parameters in water treatment
    • C02F2209/40Liquid flow rate
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2303/00Specific treatment goals
    • C02F2303/14Maintenance of water treatment installations
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00Water conservation; Efficient water supply; Efficient water use
    • Y02A20/152Water 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

Membrane bioreactor-MBR membrane pollution intelligent decision-making method
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)
② first hidden layer:
Figure BDA0002403529950000031
third, the second hidden layer:
Figure BDA0002403529950000032
output layer:
Figure BDA0002403529950000033
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:
Figure BDA0002403529950000034
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:
Figure BDA0002403529950000041
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
Figure BDA0002403529950000042
Figure BDA0002403529950000051
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:
Figure BDA0002403529950000052
wherein
Figure BDA0002403529950000053
Is the control limit for the SPE statistic,
Figure BDA0002403529950000054
λ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:
Figure BDA0002403529950000061
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,
Figure BDA0002403529950000062
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:
Figure BDA0002403529950000063
updating parameters:
Figure BDA0002403529950000064
Figure BDA0002403529950000065
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
Figure RE-GDA0002456115500000081
Figure RE-GDA0002456115500000091
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
Figure BDA0002403529950000092
Figure BDA0002403529950000101
Figure BDA0002403529950000111

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)
② first hidden layer:
Figure FDA0003313492920000021
third, the second hidden layer:
Figure FDA0003313492920000022
output layer:
Figure FDA0003313492920000023
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:
Figure FDA0003313492920000031
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:
Figure FDA0003313492920000032
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
Figure FDA0003313492920000041
Figure FDA0003313492920000042
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:
Figure FDA0003313492920000051
wherein
Figure FDA0003313492920000052
Is the control limit for the SPE statistic,
Figure FDA0003313492920000053
λ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:
Figure FDA0003313492920000061
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,
Figure FDA0003313492920000062
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:
Figure FDA0003313492920000063
updating parameters:
Figure FDA0003313492920000064
Figure FDA0003313492920000065
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.
CN202010154288.3A 2019-06-24 2020-03-07 Membrane Bioreactor-MBR Membrane Fouling Intelligent Decision-Making Method Active CN111204867B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

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