Carbon source feeding amount adjustment method, device, product and medium for denitrification process
By acquiring real-time operating parameters of the denitrification process, calculating the sludge denitrification rate and influent load, and combining the deviation of the terminal nitrate nitrogen concentration, the carbon source dosage is dynamically adjusted, solving the problem of mismatch in carbon source dosage in traditional methods, achieving precise matching between carbon source dosage and system requirements, and improving the denitrification effect.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- YANGTZE ECOLOGY & ENVIRONMENT CO LTD
- Filing Date
- 2026-02-25
- Publication Date
- 2026-06-05
AI Technical Summary
Traditional carbon source addition methods cannot accurately reflect the actual activity changes of microorganisms in the denitrification system, resulting in a mismatch between the amount of carbon source added and the actual needs, leading to problems of over- or under-addition.
By acquiring real-time operating parameters during the denitrification process, the sludge denitrification rate and influent load are calculated. Combined with the deviation of the final nitrate nitrogen concentration, the carbon source dosage is dynamically adjusted to ensure that it is precisely matched with the system requirements.
It achieves precise matching between carbon source dosage and actual needs, avoiding over- or under-dosing caused by changes in microbial activity in traditional methods, and improves the stability and efficiency of denitrification.
Smart Images

Figure CN122144894A_ABST
Abstract
Description
Technical Field
[0001] This application relates to the field of wastewater treatment technology, specifically to a method, equipment, product, and medium for adjusting the carbon source dosage in a denitrification process. Background Technology
[0002] With the increasing severity of water pollution, wastewater treatment technology has received widespread attention. In biological wastewater treatment, denitrification is a crucial step in removing nitrate nitrogen from water, playing a vital role in preventing eutrophication. Denitrification refers to the biochemical process under anaerobic conditions where denitrifying bacteria utilize organic carbon sources as electron donors to reduce nitrate nitrogen to nitrogen gas.
[0003] In actual denitrification processes, carbon source addition is a key factor affecting denitrification efficiency. Currently, engineering practices widely employ carbon source addition methods based on experience or fixed ratios. For example, CN115373256A discloses a wastewater treatment dosage control method. This method uses a multi-stage dosing approach to adjust a single PID control process into multiple PID control processes to achieve dosage control. It determines the corresponding control parameters based on the precipitated ion concentration and controls the opening and closing of the actuator to adjust the dosage.
[0004] However, since the activity of microorganisms in the denitrification system is dynamically affected by various environmental factors such as temperature, pH, and dissolved oxygen, and traditional addition methods are often based on fixed empirical parameters or simple feedback control, they cannot accurately reflect the actual denitrification capacity of the microorganisms in the current system, resulting in the carbon source addition amount not being precisely matched with the actual needs. Summary of the Invention
[0005] This application provides a method, equipment, product, and medium for adjusting the carbon source dosage in a denitrification process, which ensures a precise match between the carbon source dosage and actual requirements.
[0006] The first aspect of this application provides a method for adjusting the carbon source dosage in a denitrification process, specifically including: The real-time operating parameters of the denitrification process and the volume of the reaction tank in which the denitrification process is located are obtained. The real-time operating parameters include the influent flow rate, influent nitrate nitrogen concentration, terminal nitrate nitrogen concentration and mixed liquor suspended solids concentration. Based on real-time operating parameters and reactor volume, the sludge denitrification rate, which characterizes the current activity of denitrifying microorganisms, is obtained. The influent load is determined based on the influent flow rate and influent nitrate nitrogen concentration. The baseline carbon source dosage is determined based on the sludge denitrification rate and influent load. The carbon source feedback adjustment amount is obtained based on the deviation between the terminal nitrate nitrogen concentration and the preset nitrate nitrogen target value. The carbon source baseline dosage is then corrected based on the carbon source feedback adjustment amount to obtain the carbon source dosage.
[0007] By adopting the above technical solution, real-time operating parameters and the volume of the reaction tank in the denitrification process are obtained. These real-time operating parameters include influent flow rate, influent nitrate nitrogen concentration, terminal nitrate nitrogen concentration, and mixed liquor suspended solids concentration. Since the denitrification process is dynamically affected by environmental factors such as temperature and pH, real-time acquisition of these parameters can capture real-time changes in the system state. Based on the real-time operating parameters and the reaction tank volume, the sludge denitrification rate, which characterizes the current activity of denitrifying microorganisms, is obtained. Because the sludge denitrification rate can accurately reflect the actual nitrogen removal capacity of microorganisms under current environmental conditions, it solves the problem that traditional methods cannot accurately assess microbial activity. The influent load is determined first, as it directly determines the amount of nitrogen the system needs to remove. Accurate calculation of the influent load clarifies the actual carbon source requirement. The baseline carbon source dosage is determined based on the sludge denitrification rate and the influent load. By combining the actual treatment capacity of the microorganisms with the treatment requirements, it ensures that the carbon source dosage matches the actual needs of the system, avoiding over- or under-dosing problems caused by a fixed dosage ratio. The carbon source feedback adjustment amount is obtained based on the deviation between the terminal nitrate nitrogen concentration and the preset nitrate nitrogen target value. The baseline carbon source dosage is then corrected based on the carbon source feedback adjustment amount to obtain the carbon source dosage. Since the feedback adjustment can be corrected in real time according to the actual treatment effect, it ensures a precise match between the carbon source dosage and the actual requirements.
[0008] Optionally, obtaining the sludge denitrification rate, which characterizes the current activity of denitrifying microorganisms, based on real-time operating parameters, specifically includes: The actual nitrate removal rate is obtained by multiplying the difference between the influent nitrate nitrogen concentration and the terminal nitrate nitrogen concentration by the influent flow rate. The total amount of denitrifying sludge is obtained based on the volume of the reaction tank and the concentration of suspended solids in the mixed liquor. The ratio of the actual amount of nitrate nitrogen removed to the total amount of denitrifying sludge is used as the sludge denitrification rate.
[0009] By adopting the above technical solution, the difference between the influent nitrate nitrogen concentration and the terminal nitrate nitrogen concentration is multiplied by the influent flow rate to obtain the actual nitrate nitrogen removal amount. Since the difference between the influent nitrate nitrogen concentration and the terminal nitrate nitrogen concentration reflects the actual nitrate nitrogen concentration removed by the system, multiplying it by the influent flow rate can accurately calculate the total amount of nitrate nitrogen actually removed by the system in the current time period. The total amount of denitrifying sludge is obtained based on the reaction tank volume and the mixed liquor suspended solids concentration. Since the mixed liquor suspended solids concentration represents the amount of microorganisms per unit volume, multiplying it by the reaction tank volume can accurately calculate the total amount of microorganisms participating in the denitrification reaction in the entire reaction tank. The ratio of the actual nitrate nitrogen removal amount to the total amount of denitrifying sludge is used as the sludge denitrification rate. Since the sludge denitrification rate represents the nitrate nitrogen removal capacity per unit microorganism per unit time, the actual activity level of microorganisms under the current environmental conditions can be accurately quantified by calculating the ratio.
[0010] Optionally, determining the influent load based on the influent flow rate and influent nitrate nitrogen concentration specifically includes: Calculate the product of influent flow rate and influent nitrate nitrogen concentration to obtain the nitrate nitrogen flux entering the denitrification process per unit time with the wastewater; The homogenized concentration is obtained by arithmetically averaging the influent nitrate nitrogen concentration and the terminal nitrate nitrogen concentration. Multiplying the homogenized concentration by the influent flow rate yields the equivalent flux characterizing the overall migration of internal pollutant stock. Based on the volume of the reaction tank and the influent flow rate of the denitrification process, the hydraulic retention time is determined, and the buffer coefficient is determined based on the hydraulic retention time. The external shock load is obtained based on the nitrate nitrogen flux and the buffer coefficient; the internal inertial load is obtained based on the equivalent flux and the buffer coefficient. The influent load is determined based on external impact load and internal inertial load.
[0011] By employing the above technical solution, the product of the influent flow rate and the influent nitrate nitrogen concentration is calculated to obtain the nitrate nitrogen flux entering the denitrification process per unit time. The nitrate nitrogen flux directly reflects the nitrogen load intensity input to the system and can quantify the impact of external shocks on the system. The homogenized concentration is obtained by arithmetically averaging the influent and terminal nitrate nitrogen concentrations. During denitrification, the nitrate nitrogen concentration exhibits a gradient change from influent to terminal; the homogenized concentration represents the average concentration level throughout the entire reaction process. Multiplying the homogenized concentration by the influent flow rate yields the equivalent flux, which characterizes the overall migration of internal pollutant stock. The equivalent flux reflects the migration intensity of pollutants within the system. The method quantifies the impact of internal loads on the system. Based on the reactor volume and influent flow rate of the denitrification process, the hydraulic retention time is determined, and a buffer coefficient is determined accordingly. The hydraulic retention time determines the residence time of pollutants within the system, while the buffer coefficient reflects the system's ability to buffer load fluctuations. External shock loads are obtained based on nitrate nitrogen flux and the buffer coefficient, and internal inertial loads are obtained based on equivalent flux and the buffer coefficient. These external shock loads and internal inertial loads quantify the combined impact of external inputs and internal migrations on the system, respectively. The influent load is determined based on these two loads, comprehensively reflecting the overall nitrogen load pressure faced by the system. Optionally, determining the influent load based on external impact load and internal inertial load specifically includes: The ratio of external impact load to internal inertial load is taken as the impact inertia ratio. When the impact inertia ratio is not greater than the preset lower boundary threshold, the external impact load is squared to obtain the impact effect product. The impact effect product is divided by the internal inertia load to obtain the disturbance compensation amount. The disturbance compensation amount is added to the internal inertia load to obtain the influent load. When the impact inertia ratio is greater than the preset lower boundary threshold and less than the preset upper boundary threshold, the external impact load and the internal inertia load are geometrically averaged to obtain a fusion reference value; the external impact load and the internal inertia load are arithmetically averaged to obtain a mean reference value; the average value of the fusion reference value and the mean reference value is taken as the influent load. When the impact inertia ratio is not less than the preset upper boundary threshold, the impact inertia ratio is used as the forward amplification factor; the forward amplification factor is multiplied by the external impact load to obtain the risk response amount; the risk response amount is summed with the external impact load to obtain the influent load.
[0012] By adopting the above technical solution, the ratio of external impact load to internal inertial load is used as the impact-inertia ratio. This allows for the determination of the current load impact level faced by the system. When the impact-inertia ratio is not greater than a preset lower threshold, it indicates that the external impact is relatively weak. The impact effect product is obtained by squaring the external impact load. Squaring amplifies weak signals and prevents them from being ignored. The impact effect product is divided by the internal inertial load to obtain the disturbance compensation amount, ensuring that even weak impacts receive an appropriate response. The disturbance compensation amount is added to the internal inertial load to obtain the inflow load, avoiding insufficient response due to an excessively small impact. When the impact-inertia ratio is greater than a preset lower threshold but less than a preset upper threshold, it indicates that the external impact and internal inertia are at a comparable level. The external impact load and the internal inertial load are then geometrically averaged. By obtaining a fusion baseline value, the geometric mean avoids the influence of extreme values on the balance of weights. The arithmetic mean of the external impact load and the internal inertial load is calculated to obtain the mean reference value. The arithmetic mean provides a linear trade-off reference. The average of the fusion baseline value and the mean reference value is used as the influent load, avoiding the limitations of a single averaging method. When the impact-inertia ratio is not less than the preset upper boundary threshold, it indicates that the external impact is much greater than the internal inertia. The impact-inertia ratio is used as the look-ahead amplification factor, and the response is amplified using the impact intensity itself. The look-ahead amplification factor is multiplied by the external impact load to obtain the risk response amount, providing an additional safety margin. The risk response amount is summed with the external impact load to obtain the influent load, avoiding the risk of system overload due to strong impacts. Precise response strategies are adopted for different impact levels to ensure the accuracy of load calculation.
[0013] Optionally, determining the baseline carbon source dosage based on the sludge denitrification rate and influent load specifically includes: The activity correction coefficient is determined based on the sludge denitrification rate and the preset baseline denitrification rate. The required carbon-nitrogen ratio is determined based on the preset theoretical carbon-nitrogen ratio and the activity correction coefficient. The product of the influent load and the required carbon-nitrogen ratio is used as the baseline carbon source dosage.
[0014] By adopting the above technical solution, an activity correction coefficient is determined based on the sludge denitrification rate and a preset baseline denitrification rate. The activity correction coefficient can quantify the degree of deviation of the current microbial activity from the standard state. The required carbon-nitrogen ratio is determined based on the preset theoretical carbon-nitrogen ratio and the activity correction coefficient. The theoretical carbon-nitrogen ratio reflects the carbon source demand ratio under standard conditions. After adjustment by the activity correction coefficient, an actual carbon-nitrogen ratio that adapts to the current microbial activity state can be obtained, avoiding the problem that a fixed carbon-nitrogen ratio cannot adapt to changes in microbial activity. The product of the influent load and the required carbon-nitrogen ratio is used as the baseline carbon source dosage. The influent load determines the amount of nitrogen that needs to be treated. Multiplying it by the adjusted carbon-nitrogen ratio can accurately calculate the carbon source dosage that matches the current microbial activity and treatment requirements. This avoids the problem of over- or under-dosing caused by changes in microbial activity in traditional methods, ensuring an accurate match between the carbon source dosage and the actual demand.
[0015] Optionally, the step of obtaining the carbon source feedback adjustment amount based on the deviation between the terminal nitrate nitrogen concentration and the preset nitrate nitrogen target value specifically includes: The deviation value is obtained by calculating the difference between the terminal nitrate nitrogen concentration and the preset target value of nitrate nitrogen. Multiply the deviation value by the preset proportional response coefficient to obtain the direct adjustment component; Multiply the deviation value by the influent flow rate to obtain the compensation flux; Multiply the compensation flux by the preset load compensation coefficient to obtain the compensation adjustment component; The carbon source feedback regulation is obtained by summing the direct regulation component and the compensation regulation component.
[0016] By adopting the above technical solution, the difference between the terminal nitrate nitrogen concentration and the preset nitrate nitrogen target value is calculated to obtain the deviation value. The deviation value directly reflects the gap between the current treatment effect and the target requirement. The deviation value is multiplied by the preset proportional response coefficient to obtain the direct adjustment component. The proportional response can linearly adjust according to the magnitude of the deviation; the larger the deviation, the stronger the adjustment. The deviation value is multiplied by the influent flow rate to obtain the compensation flux. Changes in influent flow rate affect the treatment effect, and the compensation flux can quantify the amplification effect of flow rate changes on the deviation. The compensation flux is multiplied by the preset load compensation coefficient to obtain the compensation adjustment component. The load compensation coefficient can convert flux changes into specific carbon source adjustment amounts, avoiding the problem of only considering the concentration deviation while ignoring the flow rate influence. The direct adjustment component and the compensation adjustment component are summed to obtain the carbon source feedback adjustment amount. The direct adjustment component responds to the concentration deviation, while the compensation adjustment component compensates for the flow rate influence. The combination of the two can comprehensively consider the key factors affecting the treatment effect, ensuring the accuracy and comprehensiveness of the feedback adjustment.
[0017] Optionally, the step of correcting the carbon source baseline dosage based on the carbon source feedback adjustment to obtain the carbon source dosage specifically includes: Calculate the absolute value of the deviation between the terminal nitrate nitrogen concentration and the preset target value of nitrate nitrogen, and use the ratio of the absolute value of the deviation to the preset target value of nitrate nitrogen as the adjustment intensity factor characterizing the degree of deviation of the current system from steady state; Multiply the carbon source feedback adjustment amount by the adjustment intensity factor to obtain the dynamic correction amount; The carbon source dosage is obtained by algebraically summing the baseline carbon source dosage and the dynamic correction dosage.
[0018] By adopting the above technical solution, the absolute value of the deviation between the terminal nitrate nitrogen concentration and the preset nitrate nitrogen target value is calculated. The ratio of the absolute value of the deviation to the preset nitrate nitrogen target value is used as the adjustment intensity factor to characterize the degree of deviation of the current system from the steady state. The adjustment intensity factor can quantify the relative degree of the system's deviation from the target state, avoiding the problem that the absolute deviation cannot reflect the severity of the relative deviation. The carbon source feedback adjustment amount is multiplied by the adjustment intensity factor to obtain the dynamic correction amount. The adjustment intensity factor amplifies or reduces the feedback adjustment amount. The greater the degree of system deviation, the stronger the correction force. The correction force is smaller when the system is close to the target, avoiding the over-adjustment or under-adjustment problems that may be caused by fixed intensity adjustment. The carbon source baseline dosage and the dynamic correction amount are algebraically summed to obtain the carbon source dosage. The carbon source baseline dosage provides the basic dosage, and the dynamic correction amount is finely adjusted according to the actual degree of deviation. The combination of the two can ensure the basic treatment effect while making accurate corrections according to the actual situation, ensuring that the final carbon source dosage meets the basic treatment requirements and can accurately respond to changes in system state.
[0019] In a second aspect, this application provides an electronic device for adjusting the carbon source dosage in a denitrification process. The electronic device includes one or more processors and a memory. The memory is coupled to the one or more processors and is used to store computer program code, which includes computer instructions. The one or more processors invoke the computer instructions to cause the electronic device for adjusting the carbon source dosage in the denitrification process to perform the method described in the first aspect and any possible implementation thereof.
[0020] Thirdly, this application provides a computer program product containing instructions that, when run on an electronic device for adjusting the carbon source dosage in a denitrification process, cause the electronic device to perform the method described in the first aspect and any possible implementation thereof.
[0021] Fourthly, this application provides a computer-readable storage medium including instructions that, when executed on a device for adjusting the carbon source dosage in a denitrification process, cause the electronic device to perform the method described in the first aspect and any possible implementation thereof. Attached Figure Description
[0022] Figure 1 This is a schematic diagram of the architecture of a carbon source dosage adjustment system for a denitrification process provided in an embodiment of this application; Figure 2 This is a schematic flowchart of a method for adjusting the carbon source dosage in a denitrification process provided in an embodiment of this application; Figure 3 This is a schematic diagram of an influent load segmentation calculation strategy provided in an embodiment of this application; Figure 4 This is an exemplary hardware structure diagram of an electronic device for adjusting the carbon source dosage in a denitrification process, provided in an embodiment of this application. Detailed Implementation
[0023] Figure 1 An exemplary system architecture for a carbon source dosage adjustment system in a denitrification process is shown.
[0024] like Figure 1 As shown, the system architecture may include electronic device 11, network 12, and data acquisition device 13. Network 12 serves as the medium for providing a communication link between electronic device 11 and data acquisition device 13. Network 12 may include various connection types, such as wired, wireless communication links, or industrial Ethernet.
[0025] Operators can use electronic equipment 11 to interact with data acquisition device 13 via network 12 to receive or send process parameters, etc. Various process control applications, such as parameter monitoring applications and carbon source dosing control applications, can be installed on electronic equipment 11.
[0026] Electronic device 11 is hardware and can be various industrial control devices with displays, including but not limited to industrial PCs, touch screens, process controllers, and remote terminals.
[0027] The data acquisition device 13 can be a sensor system that provides various detection services, such as an online instrument that monitors parameters like influent flow rate, nitrate nitrogen concentration, and mixed liquor suspended solids concentration. The sensor system can analyze and process the acquired data and feed back the processing results (carbon source dosage adjustment instructions) to the electronic equipment.
[0028] The following detailed explanation uses the electronic device side as an example.
[0029] This embodiment provides a method for adjusting the carbon source dosage in a denitrification process. Figure 2 This is a schematic flowchart of a method for adjusting the carbon source dosage in a denitrification process provided in an embodiment of this application, as shown below. Figure 2 As shown, the method includes steps S101 to S105: S101: Obtain the real-time operating parameters of the denitrification process and the volume of the reaction tank in which the denitrification process is located. The real-time operating parameters include the influent flow rate, influent nitrate nitrogen concentration, terminal nitrate nitrogen concentration and mixed liquor suspended solids concentration.
[0030] In the embodiments of this application, denitrification refers to the biochemical process in which microorganisms use an external carbon source to convert nitrate nitrogen into nitrogen gas under anaerobic conditions, and it is widely used in wastewater treatment systems to remove nitrogen pollutants.
[0031] Real-time operating parameters represent various dynamic indicators during the denitrification process, reflecting the current system's working status and treatment efficiency. These parameters include influent flow rate, influent nitrate nitrogen concentration, final nitrate nitrogen concentration, and mixed liquor suspended solids concentration. Influent flow rate represents the amount of wastewater entering the denitrification reactor per unit time, typically expressed in cubic meters per hour. Influent nitrate nitrogen concentration refers to the nitrate nitrogen content in the wastewater entering the denitrification reactor. Final nitrate nitrogen concentration represents the remaining nitrate nitrogen content in the effluent after denitrification treatment. Mixed liquor suspended solids concentration represents the activated sludge concentration in the denitrification reactor.
[0032] Specifically, the electronic equipment collects all real-time operating data from the denitrification system through a sensor network in the data acquisition device, forming a complete parameter acquisition system. The influent flow rate is obtained through a flow meter installed on the influent pipe of the denitrification reactor, reflecting the hydraulic load of the system. The influent nitrate nitrogen concentration in the wastewater entering the denitrification reactor is obtained through an online nitrate nitrogen analyzer installed on the influent pipe, reflecting the nitrogen load entering the system. The final nitrate nitrogen concentration after treatment is obtained through an online nitrate nitrogen analyzer installed at the effluent end of the denitrification reactor, reflecting the denitrification treatment effect. The mixed liquor suspended solids concentration is obtained through a solids sensor installed in the reactor, reflecting the amount of microorganisms in the system. The reactor volume is calculated from the system design parameters or through a water level sensor.
[0033] For example, in the denitrification tank of a municipal wastewater treatment plant, electronic equipment collects comprehensive operating parameters through various online monitoring devices: electromagnetic flow meters measure the wastewater flow rate entering the denitrification tank in real time, recording the hourly wastewater volume; an online nitrate nitrogen monitor installed on the influent pipe obtains the influent nitrate nitrogen concentration, showing a concentration of tens of milligrams per liter; an online analyzer at the effluent monitoring point obtains the final nitrate nitrogen concentration after denitrification treatment, showing a single-digit milligram concentration per liter; a mixed liquor suspended solids (MSS) meter installed inside the tank obtains the activated sludge concentration, showing a concentration of several thousand milligrams per liter; and the effective volume of the tank is determined by inputting its geometric dimensions or by using a water level sensor. The electronic equipment transmits all collected parameters to the computing unit, providing necessary parameters for calculating the sludge denitrification rate and determining the carbon source dosage.
[0034] S102: Based on real-time operating parameters and reaction tank volume, the sludge denitrification rate, which characterizes the current activity of denitrifying microorganisms, is obtained.
[0035] In this embodiment, the sludge denitrification rate is a key indicator reflecting the activity of microorganisms during denitrification. It represents the ability of a unit mass of denitrifying microorganisms to remove nitrate nitrogen per unit time, directly determining the denitrification efficiency and the required carbon source. It is the theoretical basis for precise carbon source dosing and is usually expressed as milligrams of nitrate nitrogen per gram of mixed liquor suspended solids per hour. The sludge denitrification rate is not a fixed value but changes dynamically with environmental conditions, microbial community structure, and activity. Real-time calculation of the sludge denitrification rate can accurately reflect the current microbial activity state and avoid blind adjustments to the carbon source dosage.
[0036] Specifically, the electronic equipment calculates the current sludge denitrification rate using the mass balance principle. The calculation first determines the total amount of nitrate nitrogen removed from the reactor, which is calculated by multiplying the influent flow rate by the difference between the influent nitrate nitrogen concentration and the final nitrate nitrogen concentration. Then, it determines the total number of microorganisms participating in denitrification in the reactor, which is calculated by multiplying the reactor volume by the mixed liquor suspended solids concentration. Finally, the total amount of nitrate nitrogen removed is divided by the total number of microorganisms and the reaction time to obtain the nitrate nitrogen removal capacity per unit of microorganism per unit time, i.e., the sludge denitrification rate. This calculation process considers changes in flow rate, load fluctuations, and microbial concentration fluctuations under actual operating conditions, accurately reflecting the current activity state of the microorganisms.
[0037] Based on the above embodiments, as an optional embodiment, the step of obtaining the sludge denitrification rate characterizing the current denitrifying microbial activity based on real-time operating parameters and reaction tank volume may include steps S201 to S203: S201: Multiply the difference between the influent nitrate nitrogen concentration and the terminal nitrate nitrogen concentration by the influent flow rate to obtain the actual nitrate nitrogen removal rate.
[0038] In the embodiments of this application, the influent nitrate nitrogen concentration represents the content of nitrate nitrogen in the wastewater entering the denitrification reactor, reflecting the nitrogen load intensity that needs to be treated. For example, when the influent nitrate nitrogen concentration is 30 mg / L, it indicates that 30 mg of nitrate nitrogen per liter of wastewater needs to be removed.
[0039] Specifically, the electronic equipment calculates the concentration difference between the influent nitrate nitrogen concentration and the terminal nitrate nitrogen concentration through subtraction. This concentration difference represents the nitrate nitrogen removal capacity per unit volume of the denitrification process. Subsequently, the electronic equipment acquires the influent flow rate and multiplies the concentration difference by the influent flow rate. The result of this multiplication is the actual nitrate nitrogen removal amount. This calculation process achieves a conversion from a concentration-based dimension to a mass-based dimension; multiplying the concentration difference by the flow rate yields the total mass of nitrate nitrogen actually removed by the system per unit time.
[0040] S202: The total amount of denitrifying sludge is obtained based on the volume of the reaction tank and the concentration of suspended solids in the mixed liquor.
[0041] In this embodiment, the total amount of denitrifying sludge is used to represent the total mass of microbial sludge participating in the denitrification reaction in the entire reaction tank, and is a key parameter for evaluating the biochemical treatment capacity of the system.
[0042] Specifically, the electronic equipment multiplies the reactor volume by the mixed liquor suspended solids concentration, and the result of this multiplication is the total amount of denitrifying sludge. This calculation process combines the spatial information of the reactor with the sludge density information. The reactor volume provides the three-dimensional spatial range of sludge distribution, while the mixed liquor suspended solids concentration provides the mass information of sludge per unit volume. The total amount of denitrifying sludge obtained by multiplying the two represents the total amount of microorganisms available for denitrification within the entire reactor, laying the foundation for assessing the system's biochemical treatment potential and calculating the denitrification efficiency per unit of sludge.
[0043] S203: The ratio of the actual amount of nitrate nitrogen removed to the total amount of denitrifying sludge is used as the sludge denitrification rate.
[0044] Specifically, the electronic device uses the actual nitrate nitrogen removal as the dividend and the total amount of denitrifying sludge as the divisor in a division operation. The result of this division is the sludge denitrification rate. This calculation process converts the overall system performance into unit sludge efficiency. The actual nitrate nitrogen removal provides information on the total amount of nitrogen removed by the system under current operating conditions, while the total amount of denitrifying sludge provides information on the basic number of microorganisms participating in the reaction. The sludge denitrification rate obtained by dividing the two eliminates the influence of the total sludge amount and accurately reflects the intrinsic activity level of microorganisms under current conditions.
[0045] S103: Determine the influent load based on the influent flow rate and influent nitrate nitrogen concentration.
[0046] In this embodiment, the influent load represents the total load intensity that the denitrification system needs to handle. It takes into account the effects of external influent impact and internal pollutant migration, and reflects the overall treatment pressure faced by the system. For example, when the influent load is high, it indicates that the total amount of nitrogen that the system needs to handle is large and the treatment difficulty increases.
[0047] Specifically, the electronic equipment first calculates the nitrate nitrogen flux by multiplying the influent flow rate and the influent nitrate nitrogen concentration. The nitrate nitrogen flux represents the total amount of nitrate nitrogen entering the denitrification process per unit time. Then, the electronic equipment arithmetically averages the influent nitrate nitrogen concentration and the terminal nitrate nitrogen concentration to obtain the homogenized concentration. Multiplying the homogenized concentration by the influent flow rate yields the equivalent flux, which characterizes the overall migration capacity of the internal pollutant stock. Next, the electronic equipment calculates the hydraulic retention time by the ratio of the reactor volume to the current influent flow rate and converts it into a buffer coefficient to quantify the system's buffering capacity. Then, the nitrate nitrogen flux is combined with the buffer coefficient to obtain the external shock load, and the equivalent flux is combined with the buffer coefficient to obtain the internal inertial load. Finally, by analyzing the relative relationship between the external shock load and the internal inertial load, different fusion strategies are used to determine the final influent load.
[0048] Based on the above embodiments, as an optional embodiment, the step of determining the influent load based on the influent flow rate and the influent nitrate nitrogen concentration may include steps S301 to S306: S301: Calculate the product of the influent flow rate and the influent nitrate nitrogen concentration to obtain the nitrate nitrogen flux entering the denitrification process per unit time.
[0049] In the embodiments of this application, nitrate nitrogen flux represents the total mass of nitrate nitrogen flowing into the denitrification system per unit time, reflecting the nitrogen load intensity brought to the system by external wastewater input. For example, when the nitrate nitrogen flux is 25 kg / h, it indicates that 25 kg of nitrate nitrogen enters the denitrification process with the wastewater every hour.
[0050] Specifically, the electronic device uses the influent flow rate as the first multiplier and the influent nitrate nitrogen concentration as the second multiplier in a multiplication operation. The result of the multiplication operation is the nitrate nitrogen flux. This calculation process combines flow rate information in the volume dimension with pollutant information in the concentration dimension. The influent flow rate indicates the volume of wastewater entering the system per unit time, and the influent nitrate nitrogen concentration indicates the mass of nitrate nitrogen contained in each unit volume of wastewater. Multiplying the two eliminates the influence of the volume factor and directly yields the total mass of nitrate nitrogen entering the system per unit time.
[0051] S302: The homogenized concentration is obtained by arithmetically averaging the influent nitrate nitrogen concentration and the terminal nitrate nitrogen concentration.
[0052] In the embodiments of this application, the homogenization concentration represents the average value of the influent nitrate nitrogen concentration and the terminal nitrate nitrogen concentration, reflecting the overall level of nitrate nitrogen concentration during denitrification and used to characterize the balance state of pollutants within the system. For example, when the homogenization concentration is 15 mg / L, it indicates that the average level of nitrate nitrogen concentration during denitrification is 15 mg / L.
[0053] Specifically, the electronic equipment adds the influent nitrate nitrogen concentration and the terminal nitrate nitrogen concentration to obtain a total concentration, and then divides the total concentration by two to obtain the homogenized concentration. This calculation process realizes the conversion from input-output concentration to process average concentration. The influent nitrate nitrogen concentration provides the concentration information at the front end of the system, and the terminal nitrate nitrogen concentration provides the concentration information at the back end of the system. The homogenized concentration obtained by adding the two and dividing by two eliminates the influence of the concentration difference between the front and back ends and represents the typical concentration level of the entire denitrification process.
[0054] S303: Multiply the homogenized concentration by the influent flow rate to obtain the equivalent flux characterizing the overall migration of internal pollutant stock.
[0055] In the embodiments of this application, the equivalent flux represents the pollutant migration rate calculated based on the homogenized concentration. It is used to quantify the overall migration capacity of the internal pollutant stock during the denitrification process and reflects the dynamic flow characteristics of pollutants within the system. For example, when the equivalent flux is 16 kg / h, it indicates that the pollutant stock within the system migrates at a rate of 16 kg / h.
[0056] Specifically, the electronic device uses the homogenized concentration as the first multiplier and the influent flow rate as the second multiplier in a multiplication operation. The result of the multiplication operation is the equivalent flux. This calculation process combines internal equilibrium concentration information with system flow rate information. The homogenized concentration provides the typical distribution density of pollutants during denitrification, and the influent flow rate provides the system's volumetric processing capacity. The equivalent flux obtained by multiplying the two represents the migration rate of the internal pollutant stock under the current flow rate conditions, reflecting the intensity of mass transfer driven by the concentration gradient within the system.
[0057] S304: Based on the volume of the reaction tank and the influent flow rate of the denitrification process, determine the hydraulic retention time, and determine the buffer coefficient based on the hydraulic retention time.
[0058] In the embodiments of this application, the hydraulic retention time represents the average residence time of wastewater in the reaction tank, reflecting the length of time that wastewater comes into contact with microorganisms and reacts, and determining the sufficiency of the denitrification reaction. For example, when the hydraulic retention time is 2 hours, it means that the wastewater stays in the reaction tank for an average of 2 hours for denitrification treatment. The buffer coefficient is used to represent the system's ability to mitigate external shocks, quantifies the strength of the buffering effect of the reaction tank on fluctuations in the influent load, and reflects the system's shock resistance performance.
[0059] Specifically, the electronic equipment acquires two parameters: the reactor volume and the current influent flow rate. It performs a division operation using the reactor volume as the divisor and the current influent flow rate as the divisor to obtain the hydraulic retention time. This time parameter directly reflects the average residence time of wastewater in the reactor. The electronic equipment then calculates the time ratio by comparing the hydraulic retention time with a preset baseline retention time. When the time ratio is greater than 1, it indicates that the actual retention time exceeds the standard retention time, and the system has a strong buffering capacity. The electronic equipment subtracts 1 from the time ratio to obtain the excess time proportion, then multiplies the excess time proportion by the buffer gain coefficient and adds a preset basic buffer coefficient to obtain the buffer coefficient. When the time ratio is less than 1, it indicates that the actual retention time is lower than the standard retention time, and the system's buffering capacity is weak. The electronic equipment directly multiplies the time ratio by the preset basic buffer coefficient to obtain the buffer coefficient. This calculation process converts the length of the retention time into a quantitative assessment of the system's buffering performance. The time ratio reflects the comparison between actual operation and standard conditions, while the buffer gain coefficient and the basic buffer coefficient, as inherent system parameters, embody the buffering characteristics of the denitrification process. This piecewise calculation method can accurately reflect the buffering characteristics of the denitrification system under different operating conditions. When the residence time is sufficient, the system can effectively alleviate load fluctuations. The buffer coefficient increases linearly with the excess time, reflecting the cumulative effect of volume buffering. When the residence time is insufficient, the system's buffering capacity is limited. The buffer coefficient is proportional to the time ratio, reflecting the direct impact of insufficient time on buffering performance. The piecewise function avoids the limitation that a single linear relationship cannot accurately describe different operating states.
[0060] S305: Based on the nitrate nitrogen flux and buffer coefficient, the external shock load is obtained; based on the equivalent flux and buffer coefficient, the internal inertial load is obtained.
[0061] In this embodiment, the external shock load represents the instantaneous impact intensity of the external influent on the denitrification system, reflecting the actual influence of the influent nitrate nitrogen load under the buffering effect of the system. For example, when the external shock load is 18 kg / h, it indicates that the external influent generates an effective impact of 18 kg / h on the denitrification process after being buffered by the system. The internal inertial load is used to represent the load effect generated by the migration of pollutants in the system, reflecting the influence intensity of mass transfer driven by the internal concentration gradient on the system operation.
[0062] Specifically, the electronic equipment multiplies the nitrate nitrogen flux with the buffer coefficient to obtain the external shock load. The external shock load characterizes the actual shock intensity generated by the external influent under the adjustment of the system's buffer capacity. Subsequently, the electronic equipment multiplies the equivalent flux with the buffer coefficient to obtain the internal inertial load. The internal inertial load characterizes the inertial influence intensity of internal pollutant migration under the adjustment of the system's buffer capacity. The multiplication process realizes the conversion from the original flux to the effective load. The buffer coefficient, as an adjustment factor, can differentiate the load from different sources according to the current buffer capacity of the system. When the buffer coefficient is large, it indicates that the system's buffer capacity is strong, and the original flux is amplified to reflect the load carrying capacity under sufficient buffering. When the buffer coefficient is small, it indicates that the system's buffer capacity is weak, and the original flux is suppressed to reflect the load limiting effect under limited buffering.
[0063] S306: Determine the influent load based on external impact load and internal inertial load.
[0064] In the embodiments of this application, the influent load represents the comprehensive load intensity that the denitrification system needs to handle, taking into account the combined effects of external influent impact and internal pollutant migration, reflecting the overall treatment pressure and operational difficulty faced by the system.
[0065] Specifically, the electronic equipment analyzes the ratio of external impact load to internal inertial load to determine the dominant load type of the system. When the external impact load is significantly greater than the internal inertial load, it indicates that the system is dominated by external impact. The electronic equipment adopts the external-dominant mode, using the external impact load as the base value, multiplying the internal inertial load by a coordination coefficient, and adding it to the external impact load to obtain the influent load. When the internal inertial load is significantly greater than the external impact load, it indicates that the system is dominated by internal migration. The electronic equipment adopts the internal-dominant mode, using the internal inertial load as the base value, multiplying the external impact load by a coordination coefficient, and adding it to the internal inertial load to obtain the influent load. When the two loads are similar, it indicates that the system is in a balanced state. The electronic equipment adopts the balanced mode, weighting the external impact load and the internal inertial load to obtain the influent load. This classification method can accurately reflect the dominant characteristics of the system load under different operating conditions, avoiding the overestimation or underestimation of load that may be caused by simple addition. The introduction of coordination coefficients and weighting coefficients ensures the reasonable integration between different load components.
[0066] Those skilled in the art will understand that the total load of a denitrification system comes not only from newly entering wastewater (external shock load) but also from the existing stock of pollutants in the reaction tank that are migrating and transforming (internal inertial load). This internal inertial load plays a crucial buffering and stabilizing role in system operation, preventing drastic fluctuations in system state due to instantaneous external changes. Simply adding the two together ignores their differences in dominant roles under different operating conditions. Therefore, this step analyzes the relative relationship between the two and employs different fusion strategies to more accurately reflect the overall treatment pressure currently faced by the system. The above-mentioned method for calculating the buffer coefficient is only a preferred embodiment, and its core idea is to establish a functional relationship between hydraulic retention time and system buffering capacity. Other nonlinear or piecewise functions that can reflect this relationship, such as models built based on exponential or logarithmic functions, should also fall within the scope of protection of this invention. Parameters such as the baseline retention time and buffer gain coefficient can be calibrated through step response experiments or historical data analysis of a specific wastewater treatment plant.
[0067] Based on the above embodiments, as an optional embodiment, the step of determining the influent load based on external impact load and internal inertial load may include steps S401 to S404: S401: The ratio of external impact load to internal inertial load is used as the impact inertia ratio.
[0068] In the embodiments of this application, the impact inertia ratio represents the relative relationship between the external impact load and the internal inertial load, quantifies the degree of comparison between the external impact intensity and the internal migration intensity of the system, and reflects the load structure characteristics of the denitrification system.
[0069] Specifically, the electronic device acquires the two parameters calculated in the aforementioned steps: the external impact load and the internal inertial load. It then performs a division operation, using the external impact load as the dividend and the internal inertial load as the divisor, to obtain the impact-inertia ratio. This ratio directly reflects the intensity ratio of the external impact to the internal inertia. The division operation transforms the absolute load value into a relative proportional relationship. As a dimensionless parameter, the impact-inertia ratio intuitively expresses the dominant relationship between the two loads. When the ratio is greater than 1, it indicates that the external impact is dominant; when the ratio is less than 1, it indicates that the internal inertia is dominant; and when the ratio is close to 1, it indicates that the two loads are roughly equal, and the system is in a relatively balanced state.
[0070] S402: When the impact inertia ratio is not greater than the preset lower boundary threshold, the external impact load is squared to obtain the impact effect product. The impact effect product is divided by the internal inertia load to obtain the disturbance compensation amount. The disturbance compensation amount is added to the internal inertia load to obtain the influent load.
[0071] In the embodiments of this application, the impact effect product represents the cumulative effect value of the external impact load after nonlinear amplification, reflecting the intensity of the secondary influence of the external impact under the condition of internal inertia dominance. For example, when the impact effect product is 500 (kg / h)², it indicates the cumulative effect of the external impact after squaring. The disturbance compensation amount is used to represent the correction increment to the internal inertial load, quantifying the contribution of the external impact to the total system load under the internal dominance mode, and reflecting the compensation effect under weak impact conditions.
[0072] Specifically, the electronic device compares the impact inertia ratio with a preset lower boundary threshold. When the impact inertia ratio is not greater than the lower boundary threshold, it indicates that the internal inertial load is dominant. The electronic device then acquires the external impact load, performs a square operation, and multiplies the external impact load by itself to obtain the impact effect product. The impact effect product reflects the nonlinear amplification effect of the external impact on the internal stable system. Subsequently, the electronic device uses the impact effect product as the dividend and the internal inertial load as the divisor to perform a division operation to obtain the disturbance compensation amount. The division operation realizes the unit conversion from the effect product to the load increment, ensuring the consistency of the dimensions of the compensation amount and the original load. Finally, the electronic device adds the disturbance compensation amount to the internal inertial load to obtain the influent load. This calculation method is based on the physical mechanism of the denitrification system. When internal migration dominates the system, although the external impact is relatively weak, it will produce a continuous disturbance effect in the stable internal environment. The square operation simulates the cumulative amplification process of this disturbance.
[0073] Those skilled in the art will understand that when the impact inertia ratio is very low, the system operates primarily through internal circulation and is stable, with external impacts representing weak disturbances. In this case, the impact is not a simple linear superposition; using the square of the external impact load can effectively amplify this weak signal, preventing it from being drowned out by system noise. Dividing by the internal inertia load normalizes this amplified effect, matching its dimensions with the dominant internal load, thus obtaining a reasonable disturbance compensation amount that can be responded to by the system. The preset lower boundary threshold is not arbitrarily set; it can be determined through statistical analysis of historical operating data from the wastewater treatment plant. For example, the 25th percentile of the long-term distribution of the impact inertia ratio or an empirical safety value (such as 0.5) can be set as the lower boundary threshold.
[0074] S403: When the impact inertia ratio is greater than the preset lower boundary threshold and less than the preset upper boundary threshold, the external impact load and the internal inertia load are geometrically averaged to obtain a fusion reference value; the external impact load and the internal inertia load are arithmetically averaged to obtain a mean reference value; the average value of the fusion reference value and the mean reference value is used as the influent load.
[0075] In this embodiment, the fusion reference value represents the comprehensive reference value obtained by geometrically averaging the external impact load and the internal inertial load, reflecting the geometric center trend of the two loads under equilibrium conditions and having the characteristic of suppressing the influence of extreme values; the mean reference value is used to represent the arithmetic mean of the external impact load and the internal inertial load, reflecting the direct average level of the two loads and embodying the basic characteristics of linear fusion.
[0076] Specifically, the electronic device compares the impact inertia ratio with preset lower and upper boundary thresholds. When the impact inertia ratio is greater than the lower boundary threshold and less than the upper boundary threshold, it indicates that the system is in equilibrium. The electronic device acquires the external impact load and the internal inertia load, performs a geometric mean calculation, multiplies the two load values, and takes the square root to obtain the fused reference value. The fused reference value reflects the geometric center level of the two loads. Subsequently, the electronic device adds the external impact load and the internal inertia load and divides by 2 to obtain the mean reference value. The mean reference value reflects the arithmetic center level of the two loads. Finally, the electronic device adds the fused reference value and the mean reference value and divides by 2 to obtain the final influent load. This dual-average calculation method fully utilizes the different characteristics of geometric and arithmetic means. The geometric mean weakens the influence of numerical differences, while the arithmetic mean maintains the linear relationship of the values. The average value of the two achieves a balance between robustness and accuracy.
[0077] Those skilled in the art will understand that within the preset lower and upper boundary thresholds, external shocks and internal inertia are roughly equivalent, and both influence the system. This scheme employs a dual averaging strategy to achieve robust integration: the geometric mean (integration benchmark value) effectively suppresses the impact of extreme load values, resulting in a more stable outcome; while the arithmetic mean (mean reference value) directly reflects the linear contribution of both. Averaging these two averages again aims to balance the robustness of the calculation with the sensitivity to numerical changes, yielding a stable and accurate comprehensive load assessment. The preset upper boundary threshold is not arbitrarily set; it can be determined through statistical analysis of historical operating data from the wastewater treatment plant. For example, the 75th percentile of the long-term distribution of the shock-inertia ratio or an empirical value (1.5) can be set as the upper boundary threshold.
[0078] S404: When the impact inertia ratio is not less than the preset upper boundary threshold, the impact inertia ratio is used as the forward amplification factor; the forward amplification factor is multiplied by the external impact load to obtain the risk response amount; the risk response amount is summed with the external impact load to obtain the influent load.
[0079] In the embodiments of this application, the look-ahead amplification factor represents the risk amplification factor determined based on the impact inertia ratio, reflecting the degree of potential risk faced by the system under the condition of external impact dominance, and quantifying the threat level of impact intensity to system stability; the risk response amount is used to represent the increased load margin in response to external impact risks, reflecting the safety margin requirements of the system under high impact conditions.
[0080] Specifically, the electronic equipment compares the impact inertia ratio with a preset upper boundary threshold. When the impact inertia ratio is not less than the upper boundary threshold, it indicates that the system is dominated by external impacts. The electronic equipment directly uses the impact inertia ratio as the forward amplification factor, which reflects the strength of the external impact relative to the internal inertia. Subsequently, the electronic equipment multiplies the forward amplification factor with the external impact load to obtain the risk response amount. The risk response amount represents the additional load handling capacity that needs to be reserved to cope with the impact risk. Finally, the electronic equipment adds the risk response amount with the external impact load to obtain the inflow load. This calculation method is based on the forward-looking principle of risk management. When external impacts dominate, the system is prone to load fluctuations and decreased processing efficiency. Increasing the risk response amount can improve the system's shock resistance. The selection of the forward amplification factor ensures the rationality and effectiveness of the risk response.
[0081] Those skilled in the art will understand that when the impact inertia ratio is very high, the system faces significant external impact risks, which may lead to a deterioration in processing effectiveness or even system collapse. In this case, the core of the control strategy is risk anticipation and proactive response. This scheme uses the impact inertia ratio itself as a forward amplification factor, meaning that the stronger the impact, the stronger the response. By multiplying the external impact load by this factor to obtain the "risk response amount," it is equivalent to adding a safety margin proportional to the risk level on top of the base load, thereby enabling the system to prepare in advance and effectively resist strong impacts.
[0082] To more intuitively understand the strategy described above for determining influent load based on external impact load and internal inertial load, please refer to [link to relevant documentation]. Figure 3 .
[0083] Figure 3 This is a schematic diagram of an influent load segmentation calculation strategy provided in an embodiment of this application. For example... Figure 3 As shown, this strategy uses the impact inertia ratio as the core criterion. The numerical range of this index is divided into three intervals by two preset thresholds (lower boundary threshold and upper boundary threshold), corresponding to three different system operating states and calculation strategies: Internal inertia-dominated region (weak impact region): When the impact inertia ratio is not greater than the lower boundary threshold, it indicates that the system is operating stably and the external impact is a weak disturbance. At this time, a "disturbance compensation strategy" is adopted, which performs nonlinear amplification processing on the external impact load to ensure that weak signals can be effectively responded to and prevent control lag.
[0084] Balanced Transition Zone: When the impact inertia ratio is between the lower and upper boundary thresholds, it indicates that the external impact and internal inertia are evenly matched. At this point, a "robust fusion strategy" is adopted, which combines geometric and arithmetic means to smoothly and robustly fuse the two loads to obtain a stable and accurate comprehensive load assessment.
[0085] External shock dominance zone (strong shock zone): When the shock inertia ratio is not less than the upper boundary threshold, it indicates that the system faces significant external shock risk. At this time, a "proactive response strategy" is adopted, using the shock inertia ratio itself as a risk amplification factor, adding a safety margin proportional to the risk level above the base load, enabling the system to prepare in advance and effectively resist strong shocks.
[0086] Through this segmented intelligent computing strategy, the present invention can adaptively adjust the load calculation method according to the different operating conditions of the system, thereby achieving an accurate characterization of the actual processing pressure of the system.
[0087] S104: Determine the baseline carbon source dosage based on the sludge denitrification rate and influent load.
[0088] In the embodiments of this application, the carbon source baseline dosage represents the basic amount of carbon source added to meet the denitrification process, reflecting the minimum carbon source requirement to maintain normal denitrification under the current sludge activity and influent load conditions.
[0089] Electronic equipment analyzes the relationship between influent load and sludge denitrification rate to determine the system load status. When the influent load is greater than the sludge denitrification rate, it indicates that the system is in an overload state. The electronic equipment subtracts the sludge denitrification rate from the influent load to obtain the load difference. Then, it divides the load difference by the sludge denitrification rate to obtain the load excess ratio. Subsequently, it multiplies the load excess ratio by a preset carbon source dosage coefficient to obtain the basic dosage. Then, it multiplies the basic dosage by a safety margin coefficient to obtain the carbon source baseline dosage. When the influent load is less than or equal to the sludge denitrification rate, it indicates that the system is in a normal or light load state. The electronic equipment divides the influent load by the sludge denitrification rate to obtain the load sufficiency ratio. Then, it multiplies the load sufficiency ratio by a preset maintenance dosage coefficient to obtain the carbon source baseline dosage. The advantage of this classification calculation method is that it achieves dynamic matching and precise addition. The overload ratio directly reflects the degree of insufficient processing capacity, while the sufficient load ratio directly reflects the degree of surplus processing capacity. Under overload conditions, the insufficient processing capacity is compensated by adding more carbon source, while under normal load conditions, the basic processing effect is guaranteed by maintenance addition. At the same time, the classification strategy avoids the waste of carbon source caused by excessive addition and the deterioration of processing effect caused by insufficient addition.
[0090] Based on the above embodiments, as an optional embodiment, the step of determining the baseline carbon source dosage according to the sludge denitrification rate and influent load may include steps S501 to S503: S501: Determine the activity correction coefficient based on the sludge denitrification rate and the preset baseline denitrification rate.
[0091] In the embodiments of this application, the activity correction coefficient represents the degree of deviation of the current sludge activity from the standard activity level, reflects the ratio of the actual denitrification capacity to the ideal denitrification capacity, and is used to quantify the strength of sludge biological activity.
[0092] Specifically, the electronic device uses the sludge denitrification rate as the dividend and the baseline denitrification rate as the divisor to perform a division operation to obtain an initial ratio. When the initial ratio is within a preset normal range, the electronic device directly uses the initial ratio as the activity correction coefficient. When the initial ratio exceeds the normal range, the electronic device smooths the initial ratio, multiplying the excess by a decay coefficient and adding it to the normal range boundary value to obtain the activity correction coefficient. When the initial ratio is below the normal range, the electronic device enhances the initial ratio, multiplying the deficiency by an amplification coefficient and subtracting it from the normal range boundary value to obtain the activity correction coefficient. The advantage of this segmented processing method is that it effectively prevents system oscillations caused by abnormal activity. Maintaining the original ratio within the normal range accurately reflects the actual activity state. When the ratio exceeds the range, decay processing avoids excessive carbon source reduction caused by excessively high activity. When the ratio is below the range, enhancement processing avoids slow response caused by excessively low activity. The segmented strategy achieves differentiated processing and smooth transition under different activity states.
[0093] S502: Determine the required carbon-nitrogen ratio based on the preset theoretical carbon-nitrogen ratio and the activity correction coefficient.
[0094] In the embodiments of this application, the required carbon-nitrogen ratio represents the actual carbon-nitrogen ratio required for the denitrification process under the current sludge activity conditions. It reflects the optimal ratio of carbon source to nitrogen load and is used to guide the precise control of carbon source addition.
[0095] Specifically, the electronic device uses the theoretical carbon-to-nitrogen ratio as the dividend and the activity correction coefficient as the divisor to perform a division operation to obtain the required carbon-to-nitrogen ratio. This calculation process realizes the conversion from theoretical ratio to actual ratio. The division operation establishes an inverse proportional relationship between the activity level and carbon source requirement. When the activity correction coefficient is greater than 1, it indicates that the sludge activity is higher than the standard level, and the division operation makes the required carbon-to-nitrogen ratio less than the theoretical carbon-to-nitrogen ratio, reflecting the high-efficiency utilization of carbon sources by highly active sludge. When the activity correction coefficient is less than 1, it indicates that the sludge activity is lower than the standard level, and the division operation makes the required carbon-to-nitrogen ratio greater than the theoretical carbon-to-nitrogen ratio, reflecting the compensation requirement of low-active sludge for additional carbon sources. The advantage of this inverse proportional adjustment method is that it establishes a scientific correspondence between the activity level and carbon source requirement. High-active sludge can make fuller use of carbon sources for denitrification, so the carbon source addition ratio can be appropriately reduced. Low-active sludge has a weaker metabolic capacity and needs more carbon sources to maintain the same denitrification effect. The inverse proportional calculation ensures that optimal treatment effect and economic benefits can be achieved under different activity conditions.
[0096] S503: Use the product of the influent load and the required carbon-nitrogen ratio as the baseline carbon source dosage.
[0097] Specifically, the electronic equipment uses the influent load as the first multiplier to represent the total amount of nitrogen to be treated, and the required carbon-to-nitrogen ratio as the second multiplier to represent the carbon source requirement per unit of nitrogen removal. A multiplication operation is then performed to obtain the baseline carbon source dosage. This multiplication process directly achieves a quantitative conversion from nitrogen load to carbon source requirement. The influent load provides quantitative information about the treatment scale, while the required carbon-to-nitrogen ratio provides proportional information about the treatment efficiency. The result of multiplying these two values accurately reflects the total carbon source requirement of the system under current operating conditions. The advantage of this direct multiplication calculation method is that it establishes a linear correspondence between load and dosage. When the load increases, the carbon source dosage increases proportionally to ensure treatment effectiveness; when the load decreases, the carbon source dosage decreases proportionally to achieve cost savings. Simultaneously, the introduction of the required carbon-to-nitrogen ratio ensures automatic adjustment of the dosage under different activity conditions, achieving dynamic optimization and precise control.
[0098] S105: Based on the deviation between the terminal nitrate nitrogen concentration and the preset nitrate nitrogen target value, obtain the carbon source feedback adjustment amount, and correct the carbon source baseline dosage based on the carbon source feedback adjustment amount to obtain the carbon source dosage.
[0099] Specifically, the electronic device acquires two parameters: the terminal nitrate nitrogen concentration detected by the system and the preset target nitrate nitrogen value. It calculates the difference between the terminal nitrate nitrogen concentration and the target value to obtain the concentration deviation. When the concentration deviation is positive, it indicates that the actual treatment effect is not up to standard. The electronic device multiplies the concentration deviation by a preset feedback gain coefficient to obtain a positive adjustment amount as the carbon source feedback adjustment amount. When the concentration deviation is negative, it indicates that the actual treatment effect exceeds the standard. The electronic device multiplies the absolute value of the concentration deviation by a preset saving coefficient to obtain a negative adjustment amount as the carbon source feedback adjustment amount. Subsequently, the electronic device performs an algebraic operation between the carbon source feedback adjustment amount and the carbon source baseline dosage to obtain the final carbon source dosage. The advantage of this feedback adjustment method is that it achieves precise control and dynamic optimization of the treatment effect. Positive adjustment ensures timely compensation when the treatment is not up to standard, negative adjustment achieves cost savings when the treatment exceeds the standard, the setting of the feedback gain coefficient and the saving coefficient ensures the rationality of the adjustment and the stability of the system, and the algebraic operation achieves the organic combination of baseline dosage and feedback adjustment.
[0100] In this embodiment, the carbon source feedback adjustment amount represents the carbon source addition correction value determined based on the deviation between the system effluent effect and the target effect, reflecting the amount of carbon source compensation or reduction required to achieve the expected treatment effect. Based on the above embodiments, as an optional embodiment, the step of determining the baseline carbon source dosage according to the sludge denitrification rate and influent load may include steps S601 to S605: S601: Calculate the difference between the terminal nitrate nitrogen concentration and the preset target nitrate nitrogen value to obtain the deviation value.
[0101] In the embodiments of this application, the deviation value represents the numerical difference between the terminal nitrate nitrogen concentration and the preset nitrate nitrogen target value, reflecting the degree of deviation between the actual treatment effect and the expected treatment effect of the system.
[0102] Specifically, the electronic device acquires two parameters: the terminal nitrate nitrogen concentration monitored in real time and the system's preset target nitrate nitrogen value. The terminal nitrate nitrogen concentration is used as the minuend, and the target nitrate nitrogen value is used as the subtrahend to perform a subtraction operation, yielding a deviation value. This deviation value directly reflects the magnitude and direction of the difference between the current treatment effect and the target effect. The subtraction process converts absolute concentration into relative deviation. As an unbiased parameter, the deviation value objectively assesses system performance. A positive deviation value indicates that the terminal nitrate nitrogen concentration is higher than the target value, meaning the system's treatment effect is substandard and its processing capacity needs to be increased. A negative deviation value indicates that the terminal nitrate nitrogen concentration is lower than the target value, meaning the system's treatment effect is exceeding the standard and its processing intensity can be appropriately reduced. A deviation value of zero indicates that the terminal nitrate nitrogen concentration equals the target value, meaning the system's treatment effect just meets the standard.
[0103] S602: Multiply the deviation value by the preset proportional response coefficient to obtain the direct adjustment component.
[0104] In this embodiment, the direct adjustment component represents the instantaneous carbon source adjustment amount calculated based on the deviation value, reflecting the direct response adjustment intensity required for the current treatment effect deviation, and is used to achieve rapid correction of the deviation; the proportional response coefficient is used to represent the proportional relationship parameter of the conversion of the deviation value into the adjustment amount, reflecting the system's response sensitivity to the deviation.
[0105] Specifically, the electronic device acquires two parameters: the deviation value calculated in the aforementioned steps and the system's preset proportional response coefficient. The deviation value is used as the first multiplier to represent the degree of deviation requiring correction, and the proportional response coefficient is used as the second multiplier to represent the system's response intensity. A multiplication operation is then performed to obtain the direct adjustment component. This multiplication process achieves a quantitative conversion from concentration deviation to carbon source adjustment amount. The deviation value provides quantitative information on the adjustment direction and adjustment demand, while the proportional response coefficient provides proportional information on the adjustment intensity. The result of multiplying the two accurately reflects the immediate adjustment amount required by the system under the current deviation conditions. The advantage of this proportional adjustment method lies in its rapid response and accurate adjustment. It can provide an adjustment response of corresponding intensity according to the actual magnitude of the deviation, avoiding the limitations of fixed adjustment mode in adapting to different degrees of deviation. Simultaneously, the linear proportional relationship ensures the predictability of the adjustment and the stability of the control.
[0106] S603: Multiply the deviation value by the influent flow rate to obtain the compensation flux.
[0107] In the embodiments of this application, the compensation flux represents the comprehensive compensation parameter calculated based on the deviation value and the influent flow rate, reflecting the pollutant treatment flux that needs to be compensated under the current flow conditions, and is used to quantify the compensation requirements for the impact of flow rate changes on the deviation.
[0108] Specifically, the electronic equipment uses the deviation value as the first multiplier to represent the degree of deviation at the concentration level and the influent flow rate as the second multiplier to represent the treatment scale at the flow rate level, performing a multiplication operation to obtain the compensation flux. This multiplication process transforms a single concentration deviation into a comprehensive flux deviation. The deviation value provides information on the deviation in mass concentration, while the influent flow rate provides information on the scale of volumetric flow rate. The compensation flux obtained by multiplying the two comprehensively reflects the absolute deviation under the current flow rate conditions. The advantage of this flow-weighted calculation method is that it achieves flow rate normalization in deviation assessment, ensuring the fairness and accuracy of deviation response under different flow rate conditions. At high flow rates, the amplified impact of deviation reflects the severity of deviation under large-scale treatment, while at low flow rates, the reduced impact of deviation reflects the relativity of deviation under small-scale treatment.
[0109] S604: Multiply the compensation flux by the preset load compensation coefficient to obtain the compensation adjustment component.
[0110] In this embodiment, the compensation adjustment component represents the flow-weighted carbon source adjustment amount calculated based on the compensation flux, reflecting the amount of carbon source compensation required after considering the influence of influent flow, and is used to achieve accurate deviation correction under flow change conditions; the load compensation coefficient is used to represent the proportional parameter for the conversion of compensation flux into carbon source adjustment amount, reflecting the intensity setting of flux compensation.
[0111] Specifically, the electronic equipment uses the compensation flux as the first multiplier to represent the comprehensive deviation that needs to be compensated, and the load compensation coefficient as the second multiplier to represent the system's compensation response intensity. Multiplication is then performed to obtain the compensation adjustment component. This multiplication process converts the compensation flux into units of carbon source adjustment. The compensation flux provides quantitative information about the flow-weighted deviation, while the load compensation coefficient provides proportional information about the compensation intensity. The result of multiplying the two accurately reflects the amount of compensatory carbon source adjustment required by the system under the current flow conditions. The advantage of this compensation adjustment method lies in its full consideration of the impact of flow rate changes on treatment effectiveness. At high flow rates, the compensation response is enhanced to address dilution effects and load shocks; at low flow rates, the compensation response is weakened to avoid over-adjustment and resource waste. The introduction of the load compensation coefficient enables precise control and optimized configuration of the compensation intensity.
[0112] S605: Summing the direct adjustment component and the compensation adjustment component yields the carbon source feedback adjustment amount.
[0113] In the embodiments of this application, the carbon source feedback adjustment amount represents the overall carbon source adjustment amount determined after comprehensively considering the instantaneous deviation response and flow compensation effect, reflecting the complete carbon source addition correction value required based on the system feedback signal.
[0114] Specifically, the electronic device uses the direct adjustment component as the first addend to represent the basic adjustment requirement based on the deviation, and the compensation adjustment component as the second addend to represent the compensation adjustment requirement based on the flow rate. The carbon source feedback adjustment amount is obtained through addition. This addition process realizes the transformation from separate adjustment to comprehensive adjustment. The direct adjustment component provides the basic adjustment amount for deviation correction, and the compensation adjustment component provides the compensation adjustment amount for flow rate influence. The result of adding the two fully reflects the overall feedback adjustment intensity required by the system under current operating conditions. The advantage of this superimposed adjustment method is that it fully considers the multiple causes of deviations during denitrification. Concentration deviations require rapid correction through direct adjustment, while flow rate changes require adaptive correction through compensation adjustment. Superimposed adjustment ensures the comprehensiveness and effectiveness of the control strategy, achieving precise control and stable operation under complex conditions.
[0115] Based on the above embodiments, as an optional embodiment, the step of correcting the carbon source baseline dosage based on the carbon source feedback adjustment amount to obtain the carbon source dosage may include steps S701 to S703: S701: Calculate the absolute value of the deviation between the terminal nitrate nitrogen concentration and the preset nitrate nitrogen target value, and use the ratio of the absolute value of the deviation to the preset nitrate nitrogen target value as the adjustment intensity factor characterizing the degree of deviation of the current system from steady state.
[0116] In the embodiments of this application, the adjustment intensity factor represents a dimensionless index of the degree to which the current system deviates from the steady state, reflecting the severity of the relative deviation between the actual processing effect and the target effect, and is used to quantify the stability level of the system's operating state.
[0117] Specifically, the electronic device calculates the difference between the terminal nitrate nitrogen concentration and the target nitrate nitrogen value. Then, it performs an absolute value operation on this difference to obtain the absolute deviation value. This absolute value eliminates the influence of the deviation direction, retaining only the magnitude of the deviation. Subsequently, the electronic device uses the absolute deviation value as the dividend and the target nitrate nitrogen value as the divisor to perform a division operation to obtain the adjustment intensity factor. This division process converts absolute deviation into relative deviation. The absolute deviation value provides quantitative information about the degree of deviation, while the target nitrate nitrogen value provides reference information for the evaluation benchmark. The ratio obtained by dividing the two accurately reflects the severity of the deviation relative to the target value. The advantage of this relativistic evaluation method is that it establishes a unified deviation evaluation standard. Under high target value conditions, the same absolute deviation corresponds to a smaller relative deviation, while under low target value conditions, the same absolute deviation corresponds to a larger relative deviation. The relative deviation more accurately reflects the actual difficulty of system control and the urgency of the adjustment needs.
[0118] S702: Multiply the carbon source feedback adjustment amount by the adjustment intensity factor to obtain the dynamic correction amount.
[0119] In the embodiments of this application, the dynamic correction amount represents the carbon source correction dosage after dynamic adjustment based on the degree of system deviation. It reflects the actual carbon source adjustment intensity required after considering the severity of the system steady-state deviation, and is used to achieve precise control of adaptive intensity.
[0120] Specifically, the electronic device uses the carbon source feedback adjustment as the first multiplier to represent the basic adjustment requirement, and the adjustment intensity factor as the second multiplier to represent the adjustment coefficient indicating the severity of deviation. Multiplication is then performed to obtain the dynamic correction. This multiplication process transforms the system from fixed adjustment to adaptive adjustment. The carbon source feedback adjustment provides the basic intensity information, while the adjustment intensity factor provides the dynamic scaling information. The result of multiplying these two factors accurately reflects the optimal adjustment intensity required by the system at the current level of deviation. The advantage of this dynamic scaling method lies in establishing a scientific correspondence between the degree of deviation and the adjustment intensity. The more severe the deviation, the larger the adjustment intensity factor, and the dynamic correction is close to or equal to the original adjustment, ensuring strong correction. Conversely, the less severe the deviation, the smaller the adjustment intensity factor, and the dynamic correction is significantly smaller than the original adjustment, achieving gentle adjustment. Dynamic scaling ensures the accuracy and economy of the adjustment strategy.
[0121] S703: The carbon source dosage is obtained by algebraically summing the baseline carbon source dosage and the dynamic correction amount.
[0122] In the embodiments of this application, the carbon source addition amount represents the final carbon source addition amount determined by the denitrification system. It is a complete carbon source addition scheme that comprehensively considers the baseline requirements and dynamic corrections, reflecting the optimal carbon source requirements of the system under the current operating conditions.
[0123] Specifically, the electronic equipment analyzes the sign of the dynamic correction value to determine the calculation direction. When the dynamic correction value is positive, it indicates that the carbon source dosage needs to be increased. The electronic equipment performs a positive addition operation with the dynamic correction value to obtain the carbon source dosage. When the dynamic correction value is negative, it indicates that the carbon source dosage needs to be reduced. The electronic equipment performs a negative addition operation with the dynamic correction value to obtain the carbon source dosage. The algebraic calculation process realizes the transformation from separate control to integrated control. The carbon source dosage provides the main framework of control, the dynamic correction value provides fine adjustment of control, and algebraic summation ensures the accurate implementation of positive and negative corrections. The advantage of this dual control method is that it fully leverages the respective advantages of feedforward control and feedback control. Feedforward control provides rapid response and basic guarantee based on load prediction, while feedback control provides precise correction and quality improvement based on effect detection. The dual control mechanism effectively solves the limitations of a single control mode, realizing high-precision control and high-efficiency operation of the denitrification process.
[0124] To more clearly illustrate the technical solution provided in this application, a specific numerical example will be used below to simulate how electronic device 11 executes the carbon source dosage adjustment method described in this application at a certain operating moment.
[0125] I. Initial Parameter Settings Suppose that at a certain point in time, the real-time operating parameters and inherent system parameters acquired by the data acquisition device 13 are as follows: Real-time operating parameters: Influent flow rate (Q_in): 1000 m³ / h Influent nitrate nitrogen concentration (N_in): 30 mg / L Terminal nitrate nitrogen concentration (N_eff): 4.5 mg / L Mixed liquor suspended solids concentration (MLSS): 3500 mg / L (i.e., 3.5 kg / m³ or 3.5 g / L) Reaction tank volume (V): 4000 m³ Preset parameters: (Note: These parameters can be determined through statistical analysis of historical data from the wastewater treatment plant, system identification, or on-site commissioning and adjustment.) Baseline denitrification rate (SDNR_base): 2.5 mg-N / (g-MLSS·h) Theoretical carbon-to-nitrogen ratio (C / N_theory): 3.0 (assuming the external carbon source is sodium acetate) Baseline hydraulic residence time (HRT_base): 2.5 h Buffer gain coefficient: 0.5 Base buffer factor: 1.0 Impact inertia ratio lower threshold: 0.8 Impact inertia ratio upper threshold: 1.5 Target value for nitrate nitrogen (N_target): 4.0 mg / L Proportional response coefficient (Kp): 0.2 (kg / h) / (mg / L) Load compensation factor (Kc): 0.1 (dimensionless) II. Calculation Process Step 1: Calculate the sludge denitrification rate (corresponding to claim 2) Calculate the actual nitrate nitrogen removal rate: Actual nitrate nitrogen removal rate = (N_in - N_eff) × Q_in = (30 mg / L - 4.5 mg / L) × 1000 m³ / h = 25.5 mg / L × 1000 m³ / h = 25,500 g / h = 25.5 kg / h Calculate the total amount of denitrification sludge: Total denitrification sludge volume = MLSS × V = 3.5 kg / m³ × 4000 m³ = 14,000 kg Calculate the sludge denitrification rate (SDNR): SDNR = (Actual nitrate nitrogen removal rate × 1000 g / kg) / (Total denitrification sludge × 1000 g / kg) Note that the unit conversion is mg / (g·h) SDNR = (25.5 kg / h × 1,000,000 mg / kg) / (14,000 kg × 1000 g / kg) =25,500,000 / 14,000,000 ≈ 1.82 mg-N / (g-MLSS·h) Analysis: The actual denitrification activity of the microorganisms is 1.82, which is lower than the baseline value of 2.5, indicating that the microbial activity is low.
[0126] Step 2: Determine the influent load Calculate nitrate nitrogen flux and equivalent flux: Nitrate nitrogen flux = N_in × Q_in = 30 mg / L × 1000 m³ / h = 30 kg / h Homogenized concentration = (N_in + N_eff) / 2 = (30 + 4.5) / 2 = 17.25 mg / L Equivalent flux = Homogenized concentration × Q_in = 17.25 mg / L × 1000 m³ / h = 17.25 kg / h Calculate the buffer coefficient: Hydraulic retention time (HRT) = V / Q_in = 4000 m³ / 1000 m³ / h = 4.0 h Since HRT (4.0 h) > HRT_base (2.5 h), the buffer coefficient = ((HRT / HRT_base - 1) × buffer gain coefficient) + base buffer coefficient = ((4.0 / 2.5 - 1) × 0.5) + 1.0 = (1.6 - 1) × 0.5 + 1.0 = 0.3 + 1.0 = 1.3 Calculate external impact load and internal inertial load: External shock load = Nitrate nitrogen flux × Buffer coefficient = 30 kg / h × 1.3 = 39 kg / h Internal inertial load = Equivalent flux × Buffer coefficient = 17.25 kg / h × 1.3 ≈ 22.43 kg / h The influent load is determined based on the segmented impact inertia ratio: Impact inertia ratio = External impact load / Internal inertia load = 39 / 22.43 ≈ 1.74 Since the impact inertia ratio (1.74) > the preset upper boundary threshold (1.5), the system is in a state dominated by external strong impact, and the third case of claim 4 is executed.
[0127] Forward magnification factor = Impact inertia ratio = 1.74 Risk response capacity = External shock load × Forward amplification factor = 39 kg / h × 1.74 ≈ 67.86 kg / h Inflow load = External shock load + Risk response amount = 39 kg / h + 67.86 kg / h = 106.86 kg / h Analysis: The comprehensive inflow load calculated by this model is 106.86 kg / h, which is much greater than the simple nitrate nitrogen flux of 30 kg / h, reflecting a proactive response to strong shock risks.
[0128] Step 3: Determine the baseline carbon source dosage Determine the activity correction factor and required carbon-nitrogen ratio: Activity correction factor = SDNR / SDNR_base = 1.82 / 2.5 = 0.728 Required C / N ratio = C / N_theory / Activity correction factor = 3.0 / 0.728 ≈ 4.12 Analysis: Due to the low microbial activity, the required C / N ratio calculated by the model is increased from the theoretical value of 3.0 to 4.12 to compensate for the insufficient activity.
[0129] Calculate the baseline carbon source dosage: Carbon source baseline dosage = Influent load × Required carbon-to-nitrogen ratio = 106.86 kg / h × 4.12 ≈ 440.26 kg / h Step 4: Calculate the final carbon source dosage Calculate the carbon source feedback adjustment: Deviation (e) = N_eff - N_target = 4.5 mg / L - 4.0 mg / L = 0.5 mg / L Direct adjustment amount = e × Kp = 0.5 mg / L × 0.2 (kg / h) / (mg / L) = 0.1 kg / h Compensation flux = e × Q_in = 0.5 mg / L × 1000 m³ / h = 500 g / h = 0.5 kg / h Compensation adjustment component = Compensation flux × Kc = 0.5 kg / h × 0.1 = 0.05 kg / h Carbon source feedback adjustment amount = Direct adjustment component + Compensation adjustment component = 0.1 kg / h + 0.05 kg / h = 0.15 kg / h The baseline dosage is adjusted based on feedback adjustment. Adjustment intensity factor = |e| / N_target = |0.5| / 4.0 = 0.125 Dynamic correction amount = Carbon source feedback adjustment amount × Adjustment intensity factor = 0.15 kg / h × 0.125 = 0.01875 kg / h Analysis: Since the system deviates from the target value not far, the adjustment intensity factor is small. The feedback adjustment amount is scaled to avoid over-adjustment.
[0130] Calculate the final carbon source dosage: Final carbon source dosage = baseline carbon source dosage + dynamic correction dosage = 440.26 kg / h + 0.01875 kg / h ≈ 440.28 kg / h.
[0131] The following describes an electronic device for adjusting the carbon source dosage in an exemplary denitrification process provided in an embodiment of this application. Figure 4This is an exemplary hardware structure diagram of an electronic device for adjusting the carbon source dosage in a denitrification process, provided in an embodiment of this application.
[0132] In some embodiments, the electronic device for adjusting the carbon source dosage in the denitrification process is a computer device, or the electronic device for adjusting the carbon source dosage in the denitrification process includes a computer device. The computer device includes a processor, memory, and a network interface connected via a system bus. The processor of the computer device provides computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and internal memory. The non-volatile storage medium stores an operating system, computer programs, and a database. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage medium. The database of the computer device stores data. The network interface of the computer device is used to communicate with other external terminals or servers via a network connection. In some embodiments, the network interface can be a wired network interface; in some embodiments, the network interface can also be a wireless network interface. When the computer program is executed by the processor, it implements the methods in the embodiments of this application.
[0133] Those skilled in the art will understand that Figure 4 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.
Claims
1. A method for adjusting the carbon source dosage in a denitrification process, characterized in that, The method includes: The real-time operating parameters of the denitrification process and the volume of the reaction tank in which the denitrification process is located are obtained. The real-time operating parameters include the influent flow rate, influent nitrate nitrogen concentration, terminal nitrate nitrogen concentration and mixed liquor suspended solids concentration. Based on real-time operating parameters and reactor volume, the sludge denitrification rate, which characterizes the current activity of denitrifying microorganisms, is obtained. The influent load is determined based on the influent flow rate and influent nitrate nitrogen concentration. The baseline carbon source dosage is determined based on the sludge denitrification rate and influent load. The carbon source feedback adjustment amount is obtained based on the deviation between the terminal nitrate nitrogen concentration and the preset nitrate nitrogen target value. The carbon source baseline dosage is then corrected based on the carbon source feedback adjustment amount to obtain the carbon source dosage.
2. The method for adjusting the carbon source dosage in the denitrification process according to claim 1, characterized in that, The sludge denitrification rate, which characterizes the current activity of denitrifying microorganisms, is obtained based on real-time operating parameters, specifically including: The actual nitrate removal rate is obtained by multiplying the difference between the influent nitrate nitrogen concentration and the terminal nitrate nitrogen concentration by the influent flow rate. The total amount of denitrifying sludge is obtained based on the volume of the reaction tank and the concentration of suspended solids in the mixed liquor. The ratio of the actual amount of nitrate nitrogen removed to the total amount of denitrifying sludge is used as the sludge denitrification rate.
3. The method for adjusting the carbon source dosage in the denitrification process according to claim 1, characterized in that, The determination of the influent load based on the influent flow rate and influent nitrate nitrogen concentration specifically includes: Calculate the product of influent flow rate and influent nitrate nitrogen concentration to obtain the nitrate nitrogen flux entering the denitrification process per unit time with the wastewater; The homogenized concentration is obtained by arithmetically averaging the influent nitrate nitrogen concentration and the terminal nitrate nitrogen concentration. Multiplying the homogenized concentration by the influent flow rate yields the equivalent flux characterizing the overall migration of internal pollutant stock. Based on the volume of the reaction tank and the influent flow rate of the denitrification process, the hydraulic retention time is determined, and the buffer coefficient is determined based on the hydraulic retention time. The external shock load is obtained based on the nitrate nitrogen flux and the buffer coefficient; the internal inertial load is obtained based on the equivalent flux and the buffer coefficient. The influent load is determined based on external impact load and internal inertial load.
4. The method for adjusting the carbon source dosage in the denitrification process according to claim 3, characterized in that, The determination of the influent load based on external impact load and internal inertial load specifically includes: The ratio of external impact load to internal inertial load is taken as the impact inertia ratio. When the impact inertia ratio is not greater than the preset lower boundary threshold, the external impact load is squared to obtain the impact effect product. The impact effect product is divided by the internal inertia load to obtain the disturbance compensation amount. The disturbance compensation amount is added to the internal inertia load to obtain the influent load. When the impact inertia ratio is greater than the preset lower boundary threshold and less than the preset upper boundary threshold, the external impact load and the internal inertia load are geometrically averaged to obtain a fusion reference value; the external impact load and the internal inertia load are arithmetically averaged to obtain a mean reference value; the average value of the fusion reference value and the mean reference value is taken as the influent load. When the impact inertia ratio is not less than the preset upper boundary threshold, the impact inertia ratio is used as the forward amplification factor; the forward amplification factor is multiplied by the external impact load to obtain the risk response amount; the risk response amount is summed with the external impact load to obtain the influent load.
5. The method for adjusting the carbon source dosage in the denitrification process according to claim 1, characterized in that, The determination of the baseline carbon source dosage based on the sludge denitrification rate and influent load specifically includes: The activity correction coefficient is determined based on the sludge denitrification rate and the preset baseline denitrification rate. The required carbon-nitrogen ratio is determined based on the preset theoretical carbon-nitrogen ratio and the activity correction coefficient. The product of the influent load and the required carbon-nitrogen ratio is used as the baseline carbon source dosage.
6. The method for adjusting the carbon source dosage in the denitrification process according to claim 1, characterized in that, The step of obtaining the carbon source feedback adjustment amount based on the deviation between the terminal nitrate nitrogen concentration and the preset nitrate nitrogen target value specifically includes: The deviation value is obtained by calculating the difference between the terminal nitrate nitrogen concentration and the preset target value of nitrate nitrogen. Multiply the deviation value by the preset proportional response coefficient to obtain the direct adjustment component; Multiply the deviation value by the influent flow rate to obtain the compensation flux; Multiply the compensation flux by the preset load compensation coefficient to obtain the compensation adjustment component; The carbon source feedback regulation is obtained by summing the direct regulation component and the compensation regulation component.
7. The method for adjusting the carbon source dosage in the denitrification process according to claim 1, characterized in that, The process of correcting the carbon source baseline dosage based on the carbon source feedback adjustment to obtain the carbon source dosage specifically includes: Calculate the absolute value of the deviation between the terminal nitrate nitrogen concentration and the preset target value of nitrate nitrogen, and use the ratio of the absolute value of the deviation to the preset target value of nitrate nitrogen as the adjustment intensity factor characterizing the degree of deviation of the current system from steady state; Multiply the carbon source feedback adjustment amount by the adjustment intensity factor to obtain the dynamic correction amount; The carbon source dosage is obtained by algebraically summing the baseline carbon source dosage and the dynamic correction dosage.
8. An electronic device for adjusting the carbon source dosage in a denitrification process, characterized in that, The electronic device includes: one or more processors and a memory; the memory is coupled to the one or more processors, the memory being used to store computer program code, the computer program code including computer instructions, and the one or more processors invoking the computer instructions to cause the electronic device to perform the method as described in any one of claims 1-7.
9. A computer program product containing instructions, characterized in that, When the computer program product is run on an electronic device for adjusting the carbon source dosage in a denitrification process, the electronic device causes the electronic device to perform the method as described in any one of claims 1-7.
10. A computer-readable storage medium comprising instructions, characterized in that, When the instruction is executed on an electronic device for adjusting the carbon source dosage in a denitrification process, the electronic device causes the electronic device to perform the method as described in any one of claims 1-7.