A dynamic multi-objective optimization control method for wastewater treatment process

Abstract

The present application relates to sewage treatment technical field, disclose a kind of dynamic multi-objective optimization control method of sewage treatment process, comprising: step 1, by double-chamber microbial fuel cell to the effluent of regulating tank is carried out biological energy collection, and output voltage is supplied to multi-parameter sensor and edge controller by conversion circuit;Step 2, multi-parameter sensor is with 5min as period to collect sewage COD value, ammonia nitrogen concentration, dissolved oxygen content and oxidation-reduction potential data, transmission to edge controller;Step 3, edge controller is divided into low load, normal load, impact load three kinds of working condition mode according to real-time inflow and design flow ratio value.This application adopts working condition mode dynamic division and reinforcement learning control instruction generation technical scheme, reach the technical effect that aeration intensity and carbon source dosing amount are self-adaptive adjustment with inflow load, compared with fixed cycle aeration and artificial experience dosing mode in prior art, solve the deficiency of control response lag and energy consumption drug consumption virtual high.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of sewage treatment, in particular to a dynamic multi-objective optimization control method for a sewage treatment process. BACKGROUND

[0002] Due to the development of China's economic society and the improvement of people's living standards, water resources consumption and sewage discharge are increasing. In order to reduce the impact of water pollution, the growth of the national economy and the enhancement of people's environmental awareness have brought new opportunities for the development of sewage treatment process automation technology. However, the sewage treatment process is a complex dynamic system, with long biochemical reaction period, complex pollutant composition, real-time changes in water inflow and water composition, aeration energy consumption, pumping energy consumption, and mutual coupling and influence between water quality.

[0003] For example, the Chinese invention application with publication number CN119717523A discloses a dynamic multi-objective optimization control method for a sewage treatment process based on probability and dimension joint prediction, relating to the technical field of sewage treatment. It includes: obtaining process data of the sewage treatment process, wherein; constructing an optimization objective function according to a support vector regression-based predictor and the process data; and using a dynamic multi-objective optimization algorithm based on probability and dimension joint prediction and a fast elitist multi-objective genetic algorithm to solve the objective function to obtain a set of Pareto optimal solutions.

[0004] For example, the Chinese invention application with publication number CN118778522A discloses a dynamic multi-objective optimization control method and system for a sewage treatment process. First, real-time data in the sewage treatment process are collected, and a self-adaptive fuzzy neural network is used to establish an optimization objective model for water quality and energy consumption; then, a multi-objective optimization method is used to optimize and solve the optimization solution of nitrate nitrogen concentration and dissolved oxygen concentration, and an environmental detection operator is used to detect whether the environment has changed.

[0005] The above-mentioned patents have the following shortcomings:

[0006] 1. The existing control method uses fixed cycle aeration and manual experience dosing mode, and cannot establish a dynamic coupling mechanism with water inflow fluctuation and water quality concentration change, resulting in excessive aeration energy consumption and reagent dosage under low load conditions, and out-of-specification water indicators due to control response lag during high load periods, which cannot adapt to the characteristics of significant day-night differences and significant seasonal fluctuations in rural sewage discharge.

[0007] 2. The traditional equipment lacks remote state monitoring and intelligent fault tolerance mechanism, and cannot automatically switch to a safe operation mode and trigger a remote alarm when the aeration device is abnormal or the sensor fails suddenly, which increases the risk of unplanned system shutdown, and the fault recovery period is too long due to the response time of rural area maintenance personnel, which easily causes secondary pollution.

[0008] 3. The prior art does not construct a synergistic optimization model of energy consumption and drug consumption. In the long-term operation, due to the lack of process parameter self-learning ability, the actual operation energy consumption is significantly higher than the design value, the drug consumption increases unreasonably due to equipment aging, and the contradiction between the limited operation and maintenance budget of the township sewage treatment facility is increasingly prominent.

[0009] To this end, the present application provides a sewage treatment process dynamic multi-objective optimization control method to solve the problems mentioned above. SUMMARY

[0010] In view of the deficiencies of the prior art, the present application provides a sewage treatment process dynamic multi-objective optimization control method to solve the problems mentioned in the background art.

[0011] To achieve the above purpose, the present application is implemented by the following technical scheme: a sewage treatment process dynamic multi-objective optimization control method, comprising:

[0012] Step 1, biological energy is collected from the effluent of the adjusting tank by the double-chamber microbial fuel cell, and the output voltage is supplied to the multi-parameter sensor and the edge controller through the conversion circuit;

[0013] Step 2, the multi-parameter sensor collects the COD value, ammonia nitrogen concentration, dissolved oxygen content and oxidation-reduction potential data of the sewage with a period of 5 minutes, and transmits them to the edge controller;

[0014] Step 3, the edge controller divides the low load, normal load and impact load three working condition modes according to the ratio of the real-time inflow to the design flow;

[0015] Step 4, the edge controller generates a set of collaborative control instructions of aeration intensity, carbon source dosage rate and sludge return flow based on the improved Q-learning algorithm according to the current working condition mode;

[0016] Step 5, the aeration intensity is adjusted by the variable frequency micro-aeration fan, and the collaborative control instructions are executed by synchronously controlling the opening degree of the sheet-shaped slow-release bin and the return period of the peristaltic pump;

[0017] Step 6, when the data change rate of the multi-parameter sensor exceeds the set threshold, the ring-shaped ultrasonic assembly is started to clean the sensor probe and the historical control parameters are called;

[0018] Step 7, the edge controller packs and uploads the device running state data to the remote monitoring platform through the LoRa communication module;

[0019] Step 8, the experience replay pool parameters of the Q-learning algorithm are iteratively updated based on 24h running data.

[0020] Preferably, in step 1, the anode chamber of the dual-chamber microbial fuel cell is filled with a carbon felt biofilm carrier, and the cathode chamber is provided with a polytetrafluoroethylene air diffusion layer.

[0021] The operating mode classification criteria in step 3 are as follows: a real-time flow ratio of less than 0.3 is considered a low load, between 0.3 and 1.2 is considered a normal load, and greater than 1.2 is considered an impact load.

[0022] Preferably, in step 4, the improved Q-learning algorithm constructs a three-dimensional state space including flow ratio, pollutant concentration gradient, and equipment health, and defines a reward function with processing efficiency, energy consumption cost, and equipment wear as dimensions.

[0023] Preferably, step 1, the bioenergy harvesting from the effluent of the equalization tank using a dual-chamber microbial fuel cell, further includes:

[0024] Sub-step 1.1: The anode cavity is filled with a carbon felt biofilm carrier. As wastewater flows, an organic oxidation reaction occurs, generating electrons and protons. The open-circuit voltage V... oc :

[0025] V oc =E cathode -E anode ,

[0026] Among them, E cathode E is the cathode oxygen reduction potential. anode This represents the oxidation potential of the organic substrate at the anode.

[0027] Sub-step 1.2: A polytetrafluoroethylene air diffusion layer is installed in the cathode cavity. Protons pass through the proton exchange membrane and undergo a reduction reaction with oxygen, resulting in an output current density i:

[0028]

[0029] Among them, V operate R is the operating voltage. internal This refers to the battery's internal resistance.

[0030] Sub-step 1.3, when V is detected operate When the voltage drops below 3.6V, switch to supercapacitor power supply mode. Switching criteria:

[0031] V operate <V threshold and

[0032] Among them, V threshold For the preset voltage threshold, The rate of change of voltage;

[0033] Sub-step 1.4, output voltage is stabilized to 5V±0.2V by DC-DC converter, conversion efficiency η:

[0034]

[0035] Where, P out is output power, P in is input power, P loss is circuit loss power.

[0036] Preferably, in step 2, the multi-parameter sensor collects sewage parameters at a cycle of 5min further comprises:

[0037] Sub-step 2.1, start sensor self-calibration mode, calculate zero point drift compensation value of each channel:

[0038] ΔC i = C raw,i -C ref,i ,

[0039] Where, ΔC i is the zero point drift compensation value of the i channel, C raw,i is the original reading of the i channel, C ref,i is the reference value of the standard solution;

[0040] Sub-step 2.2, after the solid-state ion sieve membrane captures the target ions, calculate the ammonia nitrogen concentration C NH3 :

[0041] C NH3 = K (E m -E0) / S,

[0042] Where, K is the membrane screening coefficient, E m is the measured potential, E0 is the reference potential, and S is the Nernst slope;

[0043] Sub-step 2.3, microelectrode array detects dissolved oxygen content, and temperature compensation is carried out:

[0044] DO adj = DO raw / [1+α(T-T ref )],

[0045] Where, DO adj is the corrected value of dissolved oxygen concentration, DO raw is the original value of dissolved oxygen, α is the temperature coefficient, T is the real-time water temperature, and T ref is the standard temperature;

[0046] Sub-step 2.4, when the parameter mutation rate δ is detected to be greater than 10%, trigger abnormal processing:

[0047] δ = |(C t -C t-1 )| / (Δt·C t-1 ),

[0048] If δ exceeds the threshold, discard the current data,

[0049] wherein, δ is the data mutation rate, C t is the pollutant concentration measurement value of the current sampling period, C t-1 is the pollutant concentration history value of the last sampling period, and Δt is the time interval of adjacent sampling periods.

[0050] Sub-step 2.5, encapsulate the data by CRC-16 check and format transmission:

[0051] Data packet = [Header] [C NH3 ] [D adj ] [ORP] [CRC],

[0052] wherein, Data packet is the data frame structure, Header is the data frame header information, ORP is the oxidation-reduction potential measurement value, and CRC is the cyclic redundancy check code.

[0053] Preferably, in step 3, the edge controller divides the working condition mode further comprises:

[0054] Sub-step 3.1, calculate the real-time flow ratio η:

[0055] η = Q actual / Q design ,

[0056] wherein, Q actual is the current inflow, and Q design is the design treatment flow;

[0057] Sub-step 3.2, sliding average filtering processing is performed on η to update the effective flow ratio η eff :

[0058] η eff = (η t + η t-1 + η t-2 ) / 3,

[0059] When |η t - η t-1 | > 0.5, trigger the mutation flag:

[0060] Flag abrupt = 1,

[0061] wherein, η effη t is the real-time flow ratio of the current sampling period t-1 is the historical flow ratio of the last sampling period t-2 is the historical flow ratio of the last-but-one sampling period abrupt is the flow mutation flag

[0062] Sub-step 3.3, dividing the working condition mode:

[0063]

[0064] wherein Mode is the working condition mode classification result, η low is the low load determination threshold, η high is the impact load determination threshold

[0065] Sub-step 3.4, when Flag abrupt = 1, revising the working condition mode:

[0066] Mode = argmin(|η eff -η hist |),

[0067] wherein argmin(·) is a mathematical operator, η hist is the historical same-period flow ratio mean value

[0068] Sub-step 3.5, generating the working condition identifier code:

[0069] CaseID = [Mode][Flag abrupt ][Q actual / Q design × 100],

[0070] wherein CaseID is the working condition code

[0071] Preferably, in the step 4, the edge controller runs the improved Q-learning algorithm to generate the control instruction set further comprises:

[0072] Sub-step 4.1, constructing a three-dimensional state space M = (η eff , ΔC, H), calculating the pollutant concentration gradient ΔC:

[0073] ΔC = (C t - C t-1 ) / ΔT,

[0074] wherein C t is the current pollutant concentration, C t-1 is the last-period concentration, and ΔT is the sampling interval

[0075] Sub-step 4.2, define the device health H:

[0076] H = 1 -∑(β k ·F k ) / F max ,

[0077] where β k is the kth fault weight, F k is the number of failures, F max is the maximum allowed number of failures;

[0078] Sub-step 4.3, design the action space A = (Δf, Δα, Δq),

[0079] where Δf is the change in aeration frequency, Δα is the increase in slow-release bin opening, and Δq is the change in return flow rate;

[0080] Sub-step 4.4, construct the reward function R:

[0081] R = ω1(1 - C out / C in ) + ω2(1 - P / P max ) + ω3·H,

[0082] where ω1, ω2, ω3 are dynamic weights, C out is the effluent concentration, C in is the influent pollutant concentration, P is the real-time power, and P max is the rated power;

[0083] Sub-step 4.5, use the double-ε-greedy strategy to select actions:

[0084] ε = ε0·exp(-k·N episode ),

[0085] When rand(·) < ε, explore randomly, otherwise select the action with the maximum Q value, ε0 is the initial exploration rate, k is the decay coefficient, N episode is the cumulative number of algorithm training periods, ε is the exploration rate parameter, and rand(·) is an arbitrary number generation function;

[0086] Sub-step 4.6, apply the constraint condition to modify the action:

[0087] f current + Δf > f max , Δf = f max - f current ,

[0088] where f current is the current actual operating frequency of the aeration fan, and f max is the maximum safe frequency of the aeration fan.

[0089] Preferably, in step 5, adjusting aeration intensity by a variable frequency micro-aeration blower and executing control instructions further comprises:

[0090] Sub-step 5.1, parsing the cooperative control instruction set, calculating the target value f of aeration frequency target :

[0091] f target =f current +Δf·sgn(C out -C std ),

[0092] wherein f current is the current frequency, Δf is the aeration frequency change amount, C out is the effluent concentration, C std is the discharge standard limit value, and sgn is the sign function;

[0093] Sub-step 5.2, controlling the opening degree α of the sheet-shaped slow-release bin according to the following formula:

[0094] α=α0+k p ·(H target -H actual ),

[0095] wherein α0 is the reference opening degree, k p is the proportional coefficient, H target is the target carbon-nitrogen ratio, and H actual is the measured carbon-nitrogen ratio;

[0096] Sub-step 5.3, adjusting the peristaltic pump backflow period according to the following formula:

[0097] T pump =T base ·exp(-λ·ΔC),

[0098] wherein T pump is the peristaltic pump working period, T base is the reference period, λ is the attenuation coefficient, and ΔC is the pollutant concentration gradient;

[0099] Sub-step 5.4, the actuator linkage constraint condition is determined according to the following formula:

[0100] f target >0.8f max , α>50°,

[0101] α=50°-0.2(f target -0.8f max ),

[0102] to prevent high-frequency aeration and high-opening degree carbon injection from occurring simultaneously, wherein ftarget f is the aeration frequency target value max f is the maximum safe frequency of the aeration fan, a is the opening degree of the sheet-shaped slow-release bin;

[0103] Sub-step 5.5, update the execution state flag according to the following formula:

[0104] Flag exec = round(f target / f max ×100)∣∣round(α)×1000∣∣round(T pump ),

[0105] Wherein, Flag exec is the execution mechanism state combination flag, and round is the rounding operation.

[0106] Preferably, the step 6, when the multi-parameter sensor data change rate exceeds the set threshold, the abnormal processing further comprises:

[0107] Sub-step 6.1, calculate the data mutation rate δ according to the following formula:

[0108]

[0109] Wherein, δ is the data mutation rate, C t is the current pollutant concentration, C t-1 is the last period concentration, and ΔT is the sampling interval;

[0110] Sub-step 6.2, when δ> δ threshold , start ultrasonic cleaning according to the following formula:

[0111] T clean = min(120,30+15·[δ / 5]),

[0112] Wherein, δ threshold is the mutation threshold, and T clean is the cleaning time;

[0113] Sub-step 6.3, call the historical control parameter according to the following formula:

[0114]

[0115] Wherein, H param is the historical optimal parameter index value, C t-i is the pollutant concentration measurement value of i sampling periods before the current time t, argmin k is the mathematical operator, and C hist,k-i is the water quality data corresponding to the k-i group of parameters in the historical database;

[0116] Sub-step 6.4, data recovery is performed according to the following formula:

[0117]

[0118] wherein, C valid is the effective data output value, C hist,Hparam is the H param th parameter corresponding to the pollutant concentration data set;

[0119] Sub-step 6.5, the fault code encoding is generated according to the following formula:

[0120]

[0121] wherein, FaultCode is the fault code, round is the rounding operation, sensorID is the sensor identity code, clean count is the ultrasonic cleaning cumulative number counter.

[0122] Preferably, in step 7, the edge controller further comprises uploading device running state data through the communication module.

[0123] Sub-step 7.1, the device state matrix L is constructed according to the following formula:

[0124]

[0125] wherein, f current is the current frequency, a is the opening degree of the sheet-shaped slow-release bin, T pump is the working period of the peristaltic pump, δ is the data mutation rate, C out is the effluent concentration, C in is the influent pollutant concentration, H is the device health, FaultCode is the fault code, and CaseID is the working condition characteristic code identifier.

[0126] Sub-step 7.2, the data packet is compressed according to the following formula:

[0127]

[0128] wherein, Payload is the actual transmission data packet structure, Zigzag(L) is the device state matrix compression encoding method, is the XOR encryption, and CRC32 is the check polynomial.

[0129] Sub-step 7.3, the transmission power is dynamically adjusted according to the following formula:

[0130] P TX = min(P max , P base + 3·[RSSI / 10]),

[0131] wherein P TX is the actual transmit power control amount, P max is the rated power, P base is the reference transmit power start value, and RSSI is the received signal strength;

[0132] Sub-step 7.4, retransmit control according to the following formula:

[0133] N retry is the data packet retransmission number control, and ACK is the data packet confirmation flag;

[0134] When the confirmation frame ACK is not received, retransmit adaptively according to the signal strength;

[0135] Sub-step 7.5, update the communication state flag according to the following formula:

[0136]

[0137] wherein CommFlag is the communication state composite flag, and round is the rounding operation.

[0138] The present application provides a sewage treatment process dynamic multi-objective optimization control method. It has the following beneficial effects:

[0139] 1. The present application adopts the working condition mode dynamic division and reinforcement learning control instruction generation technical scheme, achieves the technical effect that the aeration intensity and carbon source dosage adaptively adjust with the influent load, compared with the fixed cycle aeration and artificial experience dosing mode in the prior art, solves the problems of control response lag and high energy consumption and drug consumption.

[0140] 2. The present application adopts the multi-parameter abnormality detection and remote state monitoring linkage technical scheme, achieves the technical effect of equipment fault self-diagnosis and intelligent switching of safety mode, compared with the lack of remote fault tolerance mechanism in the prior art, solves the problems of system shutdown and secondary pollution caused by sudden failure.

[0141] 3. The present application adopts the multi-objective optimization reward function and experience replay parameter iteration technical scheme, achieves the technical effect of energy consumption and drug consumption collaborative optimization and long-term running adaptation, compared with the lack of process parameter self-learning ability in the prior art, solves the problems of irrational growth of operating cost and contradiction between operation and maintenance budget. BRIEF DESCRIPTION OF DRAWINGS

[0142] Figure 1 is the flowchart of the present application. DETAILED DESCRIPTION

[0143] For those skilled in the technical field, the technical solutions in the embodiments of the present application will be described clearly and completely below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, other embodiments obtained by those of ordinary skill in the art without creative labor should be within the scope of protection of the present application.

[0144] The present application will be described in detail below in conjunction with the drawings:

[0145] Embodiments:

[0146] Please refer to the accompanying Figure 1 , the embodiment of the present application provides a sewage treatment process dynamic multi-objective optimization control method, comprising:

[0147] Step 1, the output voltage is collected by the double-chamber microbial fuel cell to the regulating pool effluent, and the output voltage is supplied to the multi-parameter sensor and the edge controller by the conversion circuit;

[0148] Substep 1.1, fill the carbon felt biofilm carrier in the anode chamber, and when the sewage flows, the organic matter oxidation reaction produces electrons and protons, and the open circuit voltage V oc :

[0149] V oc =E cathode -E anode ,

[0150] Wherein, E cathode is the cathode oxygen reduction potential, and E anode is the anode organic substrate oxidation potential;

[0151] Substep 1.2, the cathode chamber is provided with a polytetrafluoroethylene air diffusion layer, and the protons pass through the proton exchange membrane and react with oxygen to output the current density i:

[0152]

[0153] Wherein, V operate is the working voltage, and R internal is the internal resistance of the battery;

[0154] Substep 1.3, when V operate is lower than 3.6V, switch to super capacitor power supply mode, and the switching judgment condition is:

[0155] V operate <V threshold and

[0156] Wherein, V threshold is a preset voltage threshold, For voltage rate of change;

[0157] Sub-step 1.4, the output voltage is stabilized to 5V±0.2V by the DC-DC converter, and the conversion efficiency η:

[0158]

[0159] Where, P out is the output power, P in is the input power, P loss is the circuit loss power;

[0160] Step 2, the multi-parameter sensor collects the sewage COD value, ammonia nitrogen concentration, dissolved oxygen content and oxidation-reduction potential data with a period of 5 minutes, and transmits them to the edge controller;

[0161] Sub-step 2.1, start the sensor self-calibration mode, and calculate the zero point drift compensation value of each channel:

[0162] ΔC i =C raw,i -C ref,i ,

[0163] Where, ΔC i is the zero point drift compensation value of the i-th channel, C raw,i is the original reading of the i-th channel, and C ref,i is the reference value of the standard solution;

[0164] Sub-step 2.2, after the solid-state ion sieve membrane captures the target ions, calculate the ammonia nitrogen concentration C NH3 :

[0165] C NH3 =K(E m -E0) / S,

[0166] Where, K is the membrane screening coefficient, E m is the measured potential, E0 is the reference potential, and S is the Nernst slope;

[0167] Sub-step 2.3, the micro-electrode array detects the dissolved oxygen content and performs temperature compensation:

[0168] DO adj =DO raw / [1+α(T-T ref ]),

[0169] Where, DO adj is the dissolved oxygen concentration correction value, DO raw is the original value of dissolved oxygen, α is the temperature coefficient, T is the real-time water temperature, and T ref is the standard temperature;

[0170] Sub-step 2.4, when the parameter mutation rate δ>10% is detected, triggering abnormal processing:

[0171] δ = |(C t -C t-1 )| / (Δt·C t-1 ),

[0172] If δ exceeds the threshold, discard the current data,

[0173] wherein δ is the data mutation rate, C t is the pollutant concentration measurement value of the current sampling period, C t-1 is the pollutant concentration history value of the last sampling period, and Δt is the time interval of adjacent sampling periods.

[0174] Sub-step 2.5, encapsulate the data by CRC-16 check and format transmission:

[0175] Data packet = [Header] [C NH3 ] [D adj ] [ORP] [CRC],

[0176] wherein Data packet is the data frame structure, Header is the data frame header information, ORP is the oxidation-reduction potential measurement value, and CRC is the cyclic redundancy check code.

[0177] Step 3, the edge controller divides the low load, normal load and impact load three working condition modes according to the real-time water inflow and design flow ratio;

[0178] Sub-step 3.1, calculate the real-time flow ratio η:

[0179] η = Q actual / Q design ,

[0180] wherein Q actual is the current water inflow, and Q design is the design treatment flow;

[0181] Sub-step 3.2, sliding average filtering processing is performed on η, and the effective flow ratio η eff is updated:

[0182] η eff = (η t + η t-1 + η t-2 ) / 3,

[0183] When |η t -η t-1 |>0.5, trigger mutation flag:

[0184] Flag abrupt = 1,

[0185] wherein η eff is the effective flow rate ratio, η t is the real-time flow rate ratio of the current sampling period, η t-1 is the historical flow rate ratio of the last sampling period, η t-2 is the historical flow rate ratio of the second last sampling period, and Flag abrupt is the flow rate mutation flag;

[0186] Sub-step 3.3, dividing the working condition mode:

[0187]

[0188] wherein Mode is the working condition mode classification result, η low is the low load determination threshold, and η high is the impact load determination threshold;

[0189] Sub-step 3.4, when Flag abrupt = 1, correcting the working condition mode:

[0190] Mode = argmin(|η eff - η hist |),

[0191] wherein argmin(·) is a mathematical operator, and η hist is the historical same period flow rate ratio mean value;

[0192] Sub-step 3.5, generating the working condition identifier code:

[0193] CaseID = [Mode][Flag abrupt ][Q actual / Q design x 100],

[0194] wherein CaseID is the working condition code;

[0195] Step 4, the edge controller runs the improved Q-learning algorithm based on the current working condition mode to generate the aeration intensity, carbon source dosage rate, and sludge return flow collaborative control instruction set;

[0196] Sub-step 4.1, constructing a three-dimensional state space M = (η eff , ΔC, H) and calculating the pollutant concentration gradient ΔC:

[0197] ΔC = (C t - C t-1 ) / ΔT,

[0198] wherein Ct C is the current pollutant concentration t-1 C is the concentration of the last cycle, ΔT is the sampling interval

[0199] Sub-step 4.2, define the device health H:

[0200] H = 1 -∑(β k ·F k ) / F max ,

[0201] Where β k is the kth fault weight, F k is the number of failures, F max is the maximum allowed number of failures

[0202] Sub-step 4.3, design the action space A = (Δf, Δα, Δq),

[0203] Where Δf is the change in aeration frequency, Δα is the increase in slow-release bin opening, and Δq is the change in return flow rate

[0204] Sub-step 4.4, construct the reward function R:

[0205] R = ω1(1-C out / C in )+ω2(1-P / P max )+ω3·H,

[0206] Where ω1, ω2, ω3 are dynamic weights, C out is the effluent concentration, C in is the influent pollutant concentration, P is the real-time power, and P max is the rated power

[0207] Sub-step 4.5, use the double-ε-greedy strategy to select actions:

[0208] ε = ε0·exp(-k·N episode ),

[0209] When rand(·) < ε, explore randomly, otherwise select the action with the maximum Q value, ε0 is the initial exploration rate, k is the decay coefficient, N episode is the cumulative number of algorithm training cycles, ε is the exploration rate parameter, and rand(·) is an arbitrary number generation function

[0210] Sub-step 4.6, apply the constraint condition to modify the action:

[0211] f current + Δf > f max , Δf = f max - f current ,

[0212] wherein f current is the current actual operating frequency of the aeration blower, f max is the maximum safe frequency of the aeration blower;

[0213] Step 5, adjust the aeration intensity by the variable frequency micro-aeration blower, and execute the cooperative control instruction synchronously by controlling the opening degree of the sheet-shaped slow-release bin and the reflux cycle of the peristaltic pump;

[0214] Sub-step 5.1, analyze the cooperative control instruction set, and calculate the target value f target of the aeration frequency:

[0215] f target = f current + Δf· sgn(C out -C std ),

[0216] wherein f current is the current frequency, Δf is the aeration frequency change amount, C out is the effluent concentration, C std is the discharge standard limit value, and sgn is the sign function;

[0217] Sub-step 5.2, control the opening degree α of the sheet-shaped slow-release bin according to the following formula:

[0218] α = α0+ k p ·(H target -H actual ),

[0219] wherein α0is the reference opening degree, k p is the proportional coefficient, H target is the target carbon-nitrogen ratio, and H actual is the measured carbon-nitrogen ratio;

[0220] Sub-step 5.3, adjust the reflux cycle of the peristaltic pump according to the following formula:

[0221] T pump = T base · exp(-λ· ΔC),

[0222] wherein T pump is the working cycle of the peristaltic pump, T base is the reference cycle, λ is the attenuation coefficient, and ΔC is the pollutant concentration gradient;

[0223] Sub-step 5.4, the linkage constraint condition of the actuator is determined according to the following formula:

[0224] f target > 0.8f max , α > 50°,

[0225] α = 50°- 0.2(ftarget -0.8f max ),

[0226] Prevent high-frequency aeration and high-opening carbon injection from occurring simultaneously, where f target is the target value of aeration frequency, f max is the maximum safe frequency of the aeration fan, and a is the opening of the sheet-shaped slow-release bin;

[0227] Sub-step 5.5, update the execution state flag according to the following formula:

[0228] Flag exec = round(f target / f max ×100)∣∣round(α)×1000∣∣round(T pump ),

[0229] where Flag exec is the execution mechanism state combination flag, and round is the rounding operation;

[0230] Step 6, when the multi-parameter sensor data change rate exceeds the set threshold, start the ring ultrasonic assembly cleaning sensor probe and call the historical control parameters;

[0231] Sub-step 6.1, calculate the data mutation rate δ according to the following formula:

[0232]

[0233] where δ is the data mutation rate, C t is the current pollutant concentration, C t-1 is the previous period concentration, and ΔT is the sampling interval;

[0234] Sub-step 6.2, when δ > δ threshold , start ultrasonic cleaning according to the following formula:

[0235] T clean = min(120, 30 + 15·[δ / 5]),

[0236] where δ threshold is the mutation threshold, and T clean is the cleaning duration;

[0237] Sub-step 6.3, call the historical control parameters according to the following formula:

[0238]

[0239] where H param is the historical optimal parameter index value, and C t-iCi is the pollutant concentration measurement value of the current time t, i sampling period, argmin k is a mathematical operator, C hist,k-i is the water quality data corresponding to the k-i group of parameters in the historical database;

[0240] Sub-step 6.4, data recovery is performed as follows:

[0241]

[0242] Wherein, C valid is the effective data output value, is the H param th parameter corresponding to the pollutant concentration data set;

[0243] Sub-step 6.5, the fault code encoding is generated as follows:

[0244]

[0245] Wherein, FaultCode is the fault code, round is the rounding operation, sensorID is the sensor identity code, clean count is the ultrasonic cleaning cumulative number counter;

[0246] Step 7, the edge controller packs and uploads the device running state data to the remote monitoring platform through the LoRa communication module;

[0247] Sub-step 7.1, the device state matrix L is constructed as follows:

[0248]

[0249] Wherein, f current is the current frequency, α is the opening degree of the sheet-shaped slow-release bin, T pump is the working period of the peristaltic pump, δ is the data mutation rate, C out is the outlet concentration, C in is the inlet pollutant concentration, H is the device health, FaultCode is the fault code, and CaseID is the working condition characteristic code identifier;

[0250] Sub-step 7.2, the data packet is compressed as follows:

[0251]

[0252] Wherein, Payload is the actual transmission data packet structure, Zigzag(L) is the device state matrix compression encoding method, is the exclusive or encryption, and CRC32 is the check polynomial;

[0253] Sub-step 7.3, dynamically adjust the transmission power as follows:

[0254] P TX = min(P max ,P base + 3·[RSSI / 10]),

[0255] Wherein, P TX is the actual transmission power control, P max is the rated power, P base is the reference transmission power starting value, RSSI is the received signal strength;

[0256] Sub-step 7.4, retransmission control as follows:

[0257] N retry is the number of data packet retransmission control, ACK is the data packet confirmation flag;

[0258] When the confirmation frame ACK is not received, the signal strength is adaptively retransmitted;

[0259] Sub-step 7.5, update the communication status flag as follows:

[0260]

[0261] Wherein, CommFlag is the communication state composite flag, round is the rounding operation;

[0262] Step 8, based on 24h operation data iterative update Q-learning algorithm experience playback pool parameters.

[0263] In step 1, the anode cavity of the double-cavity microbial fuel cell is filled with carbon felt biofilm carrier, and the cathode cavity is provided with polytetrafluoroethylene air diffusion layer.

[0264] In step 3, the working condition mode division standard is: the real-time flow ratio is less than 0.3, which is determined as low load, 0.3 to 1.2 is the conventional load, and greater than 1.2 is determined as impact load.

[0265] In step 4, the improved Q-learning algorithm constructs a three-dimensional state space containing flow ratio, pollutant concentration gradient and equipment health degree, and defines a reward function with treatment efficiency, energy consumption cost and equipment loss as dimensions.

[0266] Through the bioenergy self-supply system to realize the equipment power supply closed loop, using microbial fuel cell to directly convert organic chemical energy into electrical energy, break through the limitation of traditional sewage treatment system relying on external power supply, significantly improve the system self-operation ability when deployed in remote areas, at the same time through intelligent power supply switching mechanism to ensure the continuous and stable operation of key equipment.

[0267] A multi-dimensional water quality perception network is constructed, dynamic calibration and abnormal data filtering technology are adopted, sensor drift and transient interference problems are effectively overcome, and accurate continuous monitoring of water quality parameters is realized, providing a high-reliability data foundation for subsequent intelligent decision-making, and greatly reducing the risk of control misalignment caused by measurement errors.

[0268] An adaptive operating condition recognition model based on flow characteristics is established, transient flow fluctuation interference is eliminated through sliding filtering and historical data comparison, and the operating load level is accurately divided, providing a classification benchmark for process parameter optimization under different inflow conditions, and solving the control response disorder problem caused by the lag of traditional method operating condition judgment.

[0269] Innovative design of multi-objective reinforcement learning control architecture, the health status of equipment is considered in decision-making, the dynamic reward mechanism is used to balance processing efficiency and operating cost, and the collaborative optimization of aeration, dosing, reflux and other execution units is realized, breaking through the limitations of isolated adjustment of each process parameter under the control mode of artificial experience.

[0270] Step 5: Develop intelligent linkage control strategy for actuator, introduce physical constraint condition to prevent equipment overload operation, through frequency adjustment and proportional opening control guided by sign function, ensure that process adjustment is quickly and accurately implemented, while maintaining treatment effect, and minimize mechanical wear and energy waste.

[0271] Build an abnormal operating condition self-recovery system, use mutation rate threshold to trigger ultrasonic cleaning and historical parameter calling mechanism, effectively deal with sudden conditions such as sensor contamination failure, and ensure continuous and stable operation of the system through multi-source data fusion repair technology, significantly reduce the ecological risk caused by unplanned shutdown.

[0272] Design an extensible Internet of Things communication architecture, realize efficient transmission in low-bandwidth environment through matrix data organization and adaptive power control, establish a panoramic visualization channel for equipment operating status, provide real-time data support for remote intelligent operation and maintenance, and solve the problem of facility supervision in rural areas.

[0273] Create an experience knowledge continuous evolution model, update algorithm parameters through historical data driving, make the control system have process adaptability and equipment aging compensation ability, ensure the continuous optimization of treatment efficiency and cost control indicators in long-term operation, and break through the technical bottleneck of performance degradation of static control model.

[0274] Although embodiments of the present application have been shown and described, it will be understood by those of ordinary skill in the art that various changes, modifications, substitutions and alterations can be made without departing from the principles and spirit of the present application, and the scope of the present application is defined by the appended claims and their equivalents.

Claims

1. A dynamic multi-objective optimization control method for wastewater treatment processes, characterized in that, The utility model relates to a kind of biological energy collection system and method for adjusting pool effluent, comprising: Step one, biological energy is collected to the effluent of adjusting pool by double-chamber microbial fuel cell, and output voltage is supplied to multi-parameter sensor and edge controller by conversion circuit; The anode chamber of the double-chamber microbial fuel cell is filled with carbon felt biofilm carrier, and the cathode chamber is provided with polytetrafluoroethylene air diffusion layer; Step two, multi-parameter sensor with Wastewater is collected periodically. Values, ammonia nitrogen concentration, dissolved oxygen content, and redox potential data are transmitted to the edge controller; Step three, edge controller is divided into three kinds of working condition modes according to the ratio of real-time water inflow and design flow, namely low load, normal load and impact load; The working condition mode division criterion is that real-time flow ratio is less than is determined as low load, to is regular load, and greater than is determined as impact load; Step four, edge controller based on the current operating mode to run the improvement The algorithm generates aeration intensity, carbon source dosing rate, and sludge return flow coordination control instruction set; The improvements The algorithm constructs a three-dimensional state space containing flow ratio, pollutant concentration gradient, and equipment health, and defines a reward function with treatment efficiency, energy consumption cost, and equipment wear as dimensions. Step five, the aeration intensity is adjusted by variable frequency micro-aeration fan, and the opening degree of sheet-shaped slow-release bin and the return period of peristaltic pump are synchronously controlled to execute cooperative control instruction, further comprising: Sub-step , parsing the cooperative control instruction set, calculating the aeration frequency target value : ， wherein is the current frequency, is the aeration frequency change amount, is the effluent concentration, is the discharge standard limit, is the symbol function; Sub-step The opening of the tablet release compartment is controlled according to the following formula : ， wherein, is the reference opening degree, is the proportional coefficient, is the target carbon-nitrogen ratio, is the measured carbon-nitrogen ratio; Sub-step The peristaltic pump backflow period is adjusted according to the following formula: ， wherein, is a period of operation of the peristaltic pump, is a reference period, is a decay coefficient, is a contaminant concentration gradient; Sub-step The actuator linkage constraint condition is determined according to the following formula: ， ， ， Prevent high frequency aeration and high opening carbon from happening at the same time, wherein, is the aeration frequency target value, is the maximum safe frequency of the aeration blower, is the sheet-shaped slow-release bin opening; sub-step the execution status flag is updated according to the following equation: ， wherein is a flag for the actuator state combination, is a rounding operation; Step six, when the data change rate of multi-parameter sensor exceeds the set threshold, the sensor probe is cleaned by starting annular ultrasonic assembly, and historical control parameters are called; Step seven, the edge controller packs the device running state data through the communication module and uploads it to the remote monitoring platform. Step seven, the edge controller packs the device running state data through the communication module and uploads it to the remote monitoring platform. Step eight, based on Running data iteration update Experience replay pool parameters for the algorithm.

2. The dynamic multi-objective optimization control method of wastewater treatment process according to claim 1, characterized in that, In the step one, further comprising: Sub-step In the anode chamber, carbon felt biofilm carrier is filled, and when sewage flows, organic matter oxidation reaction occurs to produce electrons and protons, and open circuit voltage : ， wherein, is the cathode oxygen reduction potential, is the anode organic substrate oxidation potential; Sub-step , the cathode cavity is provided with a polytetrafluoroethylene air diffusion layer, and protons pass through a proton exchange membrane to react with oxygen to output current density : ， wherein, V is the operating voltage, R is the internal resistance of the battery; Sub-step When detecting Below , switch to super capacitor power supply mode, switch determination condition: ， wherein, is a predetermined voltage threshold, is a voltage change rate; Sub-step , by the converter stabilizes the output voltage to , conversion efficiency : , wherein, Pout is the output power, Pin is the input power, Pc is the circuit loss power.

3. The dynamic multi-objective optimization control method of wastewater treatment process according to claim 1, characterized in that, In the step two, further comprising: Sub-step , start sensor self-calibration mode, calculate each channel zero-point drift compensation value: ， wherein, is the first channel zero drift compensation value, is the first channel raw reading, is the standard solution reference value; Sub-step After the target ions are captured by the solid-state ion sieve membrane, the ammonia nitrogen concentration is calculated : ， wherein, is a membrane screening coefficient, is a measured potential, is a reference potential, is a Nernst slope; Sub-step The microelectrode array detects the dissolved oxygen content and performs temperature compensation: ， wherein, is a dissolved oxygen concentration correction value, is a dissolved oxygen raw value, is a temperature coefficient, is a real-time water temperature, is a standard temperature; Sub-step When the parameter mutation rate is detected, an exception handling is triggered:​ ， If Discarding current data beyond threshold, wherein, is a data mutation rate, is a pollutant concentration measurement value of a current sampling period, is a pollutant concentration history value of a previous sampling period, is a time interval of adjacent sampling periods; Sub-step , by checking will be encapsulated, format transmission: ， wherein, is a data frame structure, is a data frame header information, is an oxidation-reduction potential measurement value, is a cyclic redundancy check code.

4. The dynamic multi-objective optimization control method of wastewater treatment process according to claim 1, characterized in that, In the step three, further comprising: Sub-step , calculating real-time flow ratio : ， wherein, Qp = current inflow flow rate, Qd = design treatment flow rate; Sub-step , the effective flow rate ratio is updated by performing a moving average filter process :​ ， When a mutation signature is triggered: ， wherein, is an effective flow rate ratio, is a real-time flow rate ratio of a current sampling period, is a history flow rate ratio of a previous sampling period, is a history flow rate ratio of a previous previous sampling period, is a flow rate mutation flag; Sub-step , dividing the working condition mode: ， wherein, is the operating mode classification result, is the low load determination threshold, is the impact load determination threshold; Sub-step When the operating mode is corrected: ， wherein, is a mathematical operator, is the historical same period flow ratio mean; sub-step generating the operating condition identifier code: ， wherein, is a working condition code.

5. The dynamic multi-objective optimization control method of wastewater treatment process according to claim 1, characterized in that, In the step four, further comprising: Sub-step , constructing a three-dimensional state space , calculating a pollutant concentration gradient : ， wherein, is the current pollutant concentration, is the concentration of the previous cycle, is the sampling interval; Sub-step , defining device health : ， wherein, is the first class failure weight, is the number of failures, is the maximum allowed number of failures; Sub-step , design action space , wherein, is the aeration frequency change amount, is the slow-release bin opening amount increase, is the backflow change rate; Sub-step , constructing a reward function : ， wherein, , , is a dynamic weight, is an effluent concentration, is an influent pollutant concentration, is a real-time power, is a rated power; Sub-step , employing a double strategy screening action: ， When < arbitrary exploration, otherwise screening value maximum action, is the initial exploration rate, is the decay coefficient, is the algorithm training period cumulative number, is the exploration rate parameter, is an arbitrary number generation function; sub-step , apply the constraint condition to correct the action: ＞ ， ， wherein, is the current actual operating frequency of the aeration blower, is the maximum safe frequency of the aeration blower.

6. The dynamic multi-objective optimization control method of wastewater treatment process according to claim 1, characterized in that, In the step six, further comprising: Sub-step The data mutation rate is calculated as follows : ， wherein, is the data mutation rate, is the current pollutant concentration, is the concentration of the previous cycle, is the sampling interval; sub-step When ultrasonic cleaning is initiated according to the following equation: ， wherein, is a mutation threshold, is a cleaning duration; Sub-step The historical control parameters are invoked as follows: ， wherein, is a historical optimal parameter index value, is a current time before a pollutant concentration measurement value of a sampling period, is a mathematical operator, is water quality data corresponding to the group parameters in the historical database is water quality data corresponding to the group parameters in the historical database Sub-step Data recovery is performed according to the following equation: ， wherein, is an effective data output value, is the first a parameter corresponding to the pollutant concentration data set; Sub-step The fault code encoding is generated according to the following equation: ， wherein, is a fault code, is a rounding operation, is a sensor identification code, is an ultrasonic cleaning cumulative number counter.

7. The dynamic multi-objective optimization control method of wastewater treatment process according to claim 1, characterized in that, In the step seven, further comprising: Sub-step , the device state matrix is constructed as follows : ， wherein, is the current frequency, is the tablet release compartment opening, is the peristaltic pump duty cycle, is the data mutation rate, is the effluent concentration, is the influent contaminant concentration, is the equipment health, is the fault code, is the operating condition signature code identifier; Sub-step compressing the data packet as follows: ， wherein, is a real transmission data packet structure, is a device state matrix compression encoding method, is an XOR encryption, is a check polynomial; Sub-step The transmit power is dynamically adjusted as follows: ， wherein, is an actual transmit power control amount, is a nominal power, is a reference transmit power start value, is a received signal strength; Sub-step , retransmission control is performed as follows: , for the number of data packet retransmissions, for the data packet acknowledgement flag; When the acknowledgement frame is not received, the retransmission is adapted to the signal strength. sub-step the communication status flag is updated according to the following equation: ， wherein is a communication status composite flag, is a rounding operation.

Images