Organic contaminated soil oxidation treatment control system with multi-parameter coupling control

Abstract

The invention relates to the technical field of treatment control, in particular to a multi-parameter coupling control organic contaminated soil oxidation treatment control system which comprises a pollution distribution gradient zoning module, an oxidizing agent block usage amount inference module, an oxidizing agent adding dynamic regulation and control module, an oxidation reaction process dynamic monitoring module and an oxidizing agent adding stop judgment module. According to the method, the pollution distribution gradient is dynamically divided, the oxidant usage amount distribution is adjusted based on the spatial distribution characteristics of the organic pollutants, and the real-time comparison of the total organic carbon content degradation rate and the oxidant residual rate in the treatment stage is combined to dynamically correct the oxidant adding proportion, so that oxidant waste and reaction imbalance are effectively avoided; through interactive analysis of the multi-layer soil oxidant concentration and the total organic carbon residual concentration change trend, diffusion limitation and degradation stagnation stages are accurately identified, based on the dynamic state of the oxidation reaction, oxidant adding is adjusted or stopped in time, and the organic pollutant removal efficiency is improved.

Description

Technical Field

[0001] The present invention relates to the technical field of process control, and particularly to an oxidation treatment control system for organically contaminated soil with multi-parameter coupling control. Background Art

[0002] The technical field of process control mainly focuses on the precise regulation and coordinated control of multiple process variables in systems such as industrial processes, environmental governance, and energy conversion, aiming to improve the overall performance and treatment effect of the system. In the direction of environmental governance, process control technology is widely applied in fields such as wastewater treatment, waste gas purification, solid waste management, and contaminated soil remediation. This technical field emphasizes establishing mathematical models, designing optimal control algorithms, and integrating various sensors and actuators to dynamically adjust parameters such as temperature, pressure, flow rate, and chemical agent dosage in the treatment process, achieving real-time and efficient control of the pollutant degradation, transformation, or removal process. Typical application schemes include algorithm systems based on PID control, adaptive control, model predictive control, fuzzy control, etc., cooperating with automation equipment and data acquisition systems to ensure the stability, reliability, and sustainability of the treatment process.

[0003] Among them, the oxidation treatment control system for organically contaminated soil with multi-parameter coupling control is a treatment control scheme for repairing organically contaminated soil. Specifically, by coupling and regulating multiple key control parameters in the oxidation reaction process, it realizes the efficient oxidative degradation of organic pollutants in the contaminated soil. The purpose of this system is to improve the removal efficiency of organic pollutants in the soil remediation process, shorten the remediation cycle, reduce energy and reagent consumption, and at the same time enhance the controllability and repeatability of the remediation process. This scheme is widely applied in scenarios such as industrial pollution sites, chemical leakage accident areas, and agricultural land remediation, and is particularly suitable for complex contaminated soil treatment projects that require precise control of the dynamic process of oxidation reaction.

[0004] Traditional control systems lack precise response to the dynamic characteristics of pollution distribution and the real-time changes in the reaction process. The control strategy mainly realizes global regulation through PID or model predictive control, and cannot be refined to local pollution intensity differences. During the treatment process, the fixed oxidant dosing strategy is likely to cause excess or deficiency of oxidant in some blocks, and there are stage imbalance problems in the oxidation reaction. It fails to identify the state of degradation stagnation or diffusion limitation in real time, resulting in increased energy consumption and reagent dosage, extended treatment cycle, large fluctuations in pollutant residues, and limited treatment effects for complex contaminated soil. Typically, the oxidation reaction in high-pollution areas decays prematurely, and the oxidant utilization rate in low-pollution areas is insufficient, with poor overall repair consistency. Summary of the Invention

[0005] The purpose of the present invention is to solve the shortcomings existing in the prior art, and propose an oxidation treatment control system for organically contaminated soil with multi-parameter coupling control.

[0006] To achieve the above object, the present invention adopts the following technical solutions: An organic contaminated soil oxidation treatment control system with multi-parameter coupling control, the system includes: The pollution distribution gradient zoning module obtains the measured values of the total organic carbon content of multiple sampling layers in the soil depth, collects the maximum increasing value and the minimum decreasing value in the depth direction, divides the soil into multiple pollution intensity blocks, and generates pollution intensity block division information; The oxidant block usage inference module, based on the pollution intensity block division information, calculates the oxidant demand of each block according to the total organic carbon content of the block, establishes the oxidant dosage for each block, and generates the oxidant usage distribution information for each block; The oxidant dosing dynamic regulation module, based on the oxidant usage distribution information for each block, obtains the time series data of the oxidant residual concentration, determines whether the oxidant consumption rate exceeds the degradation requirement of the total organic carbon content, adjusts the oxidant dosing ratio for the next round, and generates the oxidant dosing adjustment information for each block; The oxidation reaction process dynamic monitoring module, based on the oxidant dosing adjustment information for each block, compares the change trend of the oxidant concentration diffusion rate with the decline trend of the remaining concentration of the total organic carbon content, detects whether the oxidation reaction of each block enters the diffusion-limited stage or the degradation stagnation stage, and obtains the dynamic state information of the block oxidation reaction.

[0007] On the other hand, the pollution intensity block division information includes block number, block boundary position, and block pollution level. The oxidant usage distribution information for each block includes block number, required oxidant dosage for the block, and the range of oxidant penetration depth for the block. The oxidant dosing adjustment information for each block includes block number, oxidant dosing adjustment ratio, and oxidant dosing adjustment period. The dynamic state information of the block oxidation reaction includes block number, oxidant diffusion trend state, and pollutant degradation trend state.

[0008] On the other hand, the pollution distribution gradient zoning module includes: The depth concentration acquisition sub-module obtains the measured values of the total organic carbon content of multiple sampling layers in the soil depth, monitors the original values of the total organic carbon content of the sampling layers, establishes the organic carbon distribution sequence in the depth direction, and generates the soil organic carbon distribution sequence; The pollution gradient determination sub-module, according to the soil organic carbon distribution sequence, calls the maximum increasing value and the minimum decreasing value, calculates the pollution gradient factor value, determines whether the pollution gradient factor value is within the pollution area division interval, and generates the pollution intensity division basis; The block boundary demarcation sub-module, based on the pollution intensity division basis, determines whether the pollution gradients of consecutive layers continuously meet the partition threshold conditions, demarcates the upper and lower boundary positions of the block, numbers each demarcated block, and generates the pollution intensity block division information.

[0009] On the other hand, the oxidant block usage inference module includes: The block parameter acquisition sub-module collects the water content, particle size distribution, total organic carbon content, and effective porosity of each block based on the pollution intensity block division information, classifies and organizes them according to the block index, and establishes a block physical parameter set. The oxidant penetration ability determination sub-module calls the water content, particle size distribution, and effective porosity based on the block physical parameter set to determine whether the oxidant permeability coefficient is greater than the set penetration threshold, screens the blocks that meet the conditions, and obtains the oxidant penetration ability screening result. The oxidant demand calculation sub-module calculates the oxidant demand of each block by calling the total organic carbon content and the effective reaction interface area of the corresponding block according to the oxidant penetration ability screening result, summarizes the oxidant demand of each block, and establishes the block oxidant usage distribution information.

[0010] On the other hand, the oxidant dosing dynamic regulation module includes: The degradation rate acquisition sub-module obtains the block oxidant usage distribution information, collects the degradation rate sequences of the main components of the total organic carbon content in each block at the current treatment stage, detects the difference in the component degradation rate between adjacent time points, judges the trend direction of the difference, calculates the average change in the degradation rate of each block, and establishes the block degradation rate trend value. The residual rate calculation sub-module collects the oxidant residual concentration time series of each block according to the block degradation rate trend value, calculates the change difference in the residual concentration between adjacent time points, judges the trend of the change rate, calculates the average change in the residual concentration of each block, and establishes the block oxidant residual rate value. The dosing ratio adjustment sub-module calls the block degradation rate trend value based on the block oxidant residual rate value, compares the change trend of the oxidant residual rate with the degradation rate trend of the total organic carbon content, calculates the oxidant dosing adjustment coefficient for each block, compares the oxidant dosing adjustment coefficient with the adjustment threshold. If the dosing adjustment coefficient is greater than the adjustment threshold, increase the oxidant dosing ratio; otherwise, maintain the current oxidant dosing ratio, and establish the block oxidant dosing adjustment information.

[0011] On the other hand, the dynamic monitoring module of the oxidation reaction process includes: The oxidant residual trend extraction sub-module collects the measured values of the oxidant residual concentration at multiple depths of the soil based on the block oxidant dosing adjustment information, calculates the difference in the oxidant concentration between adjacent depth layers, judges the change direction of the concentration difference, calculates the diffusion rate of the oxidant in each layer within the time series, screens the change rate of the oxidant concentration diffusion trend within the block, and establishes the oxidant diffusion trend change rate information. The organic carbon concentration trend extraction sub-module collects the remaining concentration change value of the total organic carbon content in the corresponding soil layer according to the oxidant diffusion trend change rate information, calculates the difference in the decrease of the organic carbon concentration in adjacent time periods, judges the stability of the decreasing trend, calculates the decreasing rate of the total organic carbon concentration in the time series, and establishes the information on the decreasing rate of the organic carbon concentration. Based on the information on the decreasing rate of the organic carbon concentration, the oxidation reaction state determination sub-module calls the oxidant diffusion trend change rate, compares the difference in trends between the two, calculates the oxidation reaction state change factor, judges whether the change factor exceeds the diffusion limitation determination threshold or the degradation stagnation determination threshold, obtains the state of the oxidation reaction stage of each block, and establishes the dynamic state information of the block oxidation reaction.

[0012] On the other hand, the system further includes: Based on the dynamic state information of the block oxidation reaction, the oxidant dosing stop determination module detects whether the change in the remaining concentration of the total organic carbon content is lower than the degradation rate threshold of the total organic carbon content. If both conditions are met, it determines to stop dosing for the corresponding block, establishes the timing strategy for stopping dosing for the entire block, and generates the block oxidant dosing stop instruction set. The block oxidant dosing stop instruction set includes the block number, the stop dosing time point, and the stop dosing confirmation flag.

[0013] On the other hand, the oxidant dosing stop determination module includes: Based on the dynamic state information of the block oxidation reaction, the diffusion ability determination sub-module obtains the current oxidant residual concentration gradient of each block, collects the change value of the block oxidant diffusion rate, calculates the difference in the diffusion rate in adjacent time periods, judges the change trend of the diffusion ability, compares the trend change rate with the diffusion ability threshold, screens the blocks with insufficient diffusion ability, and establishes a list of blocks with insufficient diffusion ability. Based on the list of blocks with insufficient diffusion ability, the organic carbon degradation state determination sub-module collects the change value of the remaining concentration of the total organic carbon content in the corresponding block, calculates the degradation rate of the total organic carbon content in adjacent time periods, judges whether the degradation rate is lower than the degradation rate threshold of the total organic carbon content, screens the blocks that meet the degradation stagnation state, and establishes a list of blocks in the degradation stagnation state. Based on the list of blocks in the degradation stagnation state, the dosing stop strategy formulation sub-module calls the list of blocks with insufficient diffusion ability, compares the block numbers in the two lists, screens the blocks that simultaneously meet the conditions of insufficient diffusion ability and degradation stagnation state, assigns the stop dosing time point, generates the stop dosing order of the blocks, and establishes the block oxidant dosing stop instruction set.

[0014] Compared with the prior art, the advantages and positive effects of the present invention are as follows: In the present invention, by dynamically dividing the pollution distribution gradient, adjusting the distribution of the oxidant usage amount based on the spatial distribution characteristics of organic pollutants, and combining the real-time comparison of the degradation rate of the total organic carbon content and the residual rate of the oxidant in the treatment stage, the dosing ratio of the oxidant is dynamically corrected, effectively avoiding the waste of the oxidant and the reaction imbalance. Through the cross-analysis of the change trends of the oxidant concentration in multiple soil layers and the remaining concentration of the total organic carbon, the diffusion-limited and degradation-stagnant stages are accurately identified. Based on the dynamic state of the oxidation reaction, the dosing of the oxidant is adjusted or stopped in a timely manner, ensuring the high efficiency and controllability of the oxidation process, improving the removal efficiency of organic pollutants, shortening the treatment cycle, reducing the consumption of energy and reagents, enhancing the repeatability and stability of the repair process, and meeting the dynamic treatment requirements of complex contaminated soil. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] Figure 1 It is a system flowchart of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0016] In order to make the objectives, technical solutions and advantages of the present invention clearer and more understandable, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0017] In the description of the present invention, it should be understood that the orientation or positional relationship indicated by the terms "length", "width", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", etc. is based on the orientation or positional relationship shown in the drawings, and is only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore should not be construed as a limitation of the present invention. In addition, in the description of the present invention, the meaning of "a plurality of" is two or more unless otherwise specifically defined.

[0018] Please refer to Figure 1 , the present invention provides a technical solution: an oxidation treatment control system for organic contaminated soil with multi-parameter coupling control, the system includes a pollution distribution gradient zoning module, an oxidant block usage amount inference module, an oxidant dosing dynamic regulation module, a dynamic monitoring module for the oxidation reaction process, and an oxidant dosing stop determination module; The pollution distribution gradient zoning module obtains the measured values of the total organic carbon content of multiple sampling levels in the soil depth, detects the change range of the total organic carbon content between each sampling level, collects the maximum increasing value and the minimum decreasing value in the depth direction, calculates the ratio of the two as a determination factor, determines whether the determination factor is greater than the pollution distribution zoning threshold, and if the condition is met, divides the soil into multiple pollution intensity blocks and generates pollution intensity block division information. Total organic carbon content is an international standard soil pollution detection index, commonly used to measure the overall level of soil organic pollutants; Based on the pollution intensity block division information, the oxidant block usage inference module collects the water content, particle size distribution, total organic carbon content, and effective porosity of each block, determines the oxidant permeability coefficient and the effective reaction interface, calculates the oxidant demand of the block according to the total organic carbon content of the block, establishes the oxidant dosage of the block, and generates the oxidant usage distribution information of the block; The oxidant permeability coefficient can be obtained through the joint test of oxidant type and soil physical properties, reflecting the effective migration rate of the oxidant under a specific soil structure; Based on the oxidant usage distribution information of the block, the oxidant dosing dynamic regulation module collects the degradation rate sequences of the main components of the total organic carbon content of each block in the current treatment stage, obtains the time series data of the oxidant residual concentration, calculates the change rate of the oxidant residual concentration within the time window, obtains the change trend of the total organic carbon content degradation rate sequence, compares the oxidant residual rate with the total organic carbon content degradation rate, and determines whether the oxidant consumption rate exceeds the degradation demand of the total organic carbon content. If it exceeds, the oxidant dosing ratio for the next round is reduced; if it does not exceed, the current dosing ratio is maintained, and the oxidant dosing adjustment information for the block is generated; Based on the oxidant dosing adjustment information of the block, the dynamic monitoring module of the oxidation reaction process collects the measured values of the oxidant residual concentration at multiple depths of the soil, obtains the change value of the remaining concentration of the total organic carbon content in the corresponding soil layer, compares the change trend of the oxidant concentration diffusion rate with the downward trend of the remaining concentration of the total organic carbon content, detects whether the oxidation reaction of each block enters the diffusion-limited stage or the degradation stagnation stage, and obtains the dynamic state information of the block oxidation reaction; Based on the dynamic state information of the block oxidation reaction, the oxidant dosing stop determination module determines whether the diffusion ability of persulfate in the block is insufficient to maintain the effective degradation of the total organic carbon content, and detects whether the change in the remaining concentration of the total organic carbon content is lower than the degradation rate threshold of the total organic carbon content. If both conditions are met, it is determined that the dosing of the corresponding block is stopped, the full-block stop dosing time sequence strategy is established, and the oxidant stop dosing instruction set for the block is generated; The degradation rate threshold of the total organic carbon content is the pollutant treatment control standard set for the project, generally evaluated based on the repair target value and the dynamic degradation trend, and often determined in combination with the national soil pollution risk control standard; The pollution intensity block division information includes block number, block boundary position, and block pollution level. The block oxidant usage distribution information includes block number, required oxidant dosage for the block, and the range of oxidant penetration depth for the block. The block oxidant dosing adjustment information includes block number, oxidant dosing adjustment ratio, and oxidant dosing adjustment period. The block oxidation reaction dynamic state information includes block number, oxidant diffusion trend state, and pollutant degradation trend state. The block oxidant stop dosing instruction set includes block number, stop dosing time point, and stop dosing confirmation flag.

[0019] The pollution distribution gradient zoning module includes: The depth concentration acquisition sub-module obtains the measured values of the total organic carbon content at multiple sampling levels in the soil depth, monitors the original values of the total organic carbon content at the sampling levels, establishes the organic carbon distribution sequence in the depth direction, and generates the soil organic carbon distribution sequence; Obtain the measured values of the total organic carbon content at multiple sampling levels in the soil depth, monitor the original values of the total organic carbon content at each sampling level, and establish the organic carbon distribution sequence in the depth direction. The specific process is as follows: First, set the soil profile depth range from 0 cm to 100 cm, set the depth interval of the sampling levels to 10 cm, and sequentially collect 10 layers from 0 - 10 cm, 10 - 20 cm, 20 - 30 cm... up to 90 - 100 cm. Five replicate points are arranged at each sampling level, and the samples at each sampling level are tested for content using the high-temperature combustion method to obtain the original content values. During the monitoring process, abnormal points are excluded, and values with a relative deviation exceeding 10% in repeated detections are removed. The mean value is calculated using the remaining valid values to obtain the mean value for each layer. For example, the five measured values for the 0 - 10 cm layer are 12.1 g / kg, 12.4 g / kg, 11.9 g / kg, 12.2 g / kg, and 12.0 g / kg, and the calculated mean value is 12.12 g / kg. The mean values for the remaining layers are calculated in the same way, and the depth direction content sequence {12.12, 11.80, 10.75, 10.50, 9.80, 8.90, 7.65, 6.80, 5.95, 5.30} g / kg is sorted and summarized. Since no thresholds, reference values, weights, or coefficients appear in this paragraph, no additional supplements are required, and the soil organic carbon distribution sequence is established.

[0020] The pollution gradient determination sub-module calls the maximum increasing value and the minimum decreasing value according to the soil organic carbon distribution sequence, and uses the formula: ; operates to obtain the pollution gradient factor value, determines whether the pollution gradient factor value is within the pollution area division interval, and generates the basis for pollution intensity division; in, represents the pollution gradient factor value, is the maximum incremental difference between all adjacent layers in the depth direction, is the minimum decreasing difference, is the normalized value of the sampling layer spacing, reflecting the impact of spatial sampling layer distribution on gradient determination. is the normalized value of the average total organic carbon content in the profile section, is the number of profile sampling layers, It is the normalized value of the coefficient of variation of soil density, reflecting the impact of soil structure on pollution distribution and enhancing calculation flexibility; According to the soil organic carbon distribution sequence, the maximum increment and minimum decrement values are called, and the intervals between adjacent layers are first calculated. Difference, using the difference calculation formula , calculate 9 sets of difference sequences, and filter the maximum increment value in the difference and minimum decrement value , example sequence calculation, , , normalized value of sampling layer spacing , mean , the number of sampling layers , normalized value of soil density variation coefficient , substitute into the formula: ; in, Represents the pollution gradient factor value, which is the main evaluation index of this step. is the maximum incremental difference between all adjacent layers in the depth direction, is the minimum decreasing difference, is the normalized value of the sampling layer spacing, reflecting the uniformity of the profile space sampling density. is the normalized value of the average total organic carbon content in the profile section, is the number of profile sampling layers, It is the normalized value of soil density variation coefficient, reflecting the influence of soil structure on pollution distribution and enhancing calculation flexibility. Used to quantify the difference in pollution increase and decrease in profile depth, multiplication Used to enhance the uniformity of spatial sampling, the denominator By normalizing the pollution load and the influence of the number of sampling layers through exponentiation and square root, The regulatory effect of nonlinear soil structure variation on the stability of pollution distribution is mapped logarithmically, and the overall formula uses division to balance the difference between the numerator and denominator to construct a gradient sensitivity index.

[0021] The advantage of the formula is that by introducing and two normalization terms, it breaks through the simple difference evaluation mode of the traditional TOC gradient, can dynamically reflect the actual impact of the profile sampling density and soil density variation on the pollution gradient, and improve the accuracy and adaptability of pollution zoning determination.

[0022] The denominator part of the formula calculation , the numerator part , calculate . The zoning interval is set to , this interval refers to the recommended interval of the light pollution section gradient of the organic carbon change grading standard of the soil profile, and is set in combination with the average gradient data of the actual site profile test. The threshold interval does not depend on historical experience, but is set according to the recommended section of the standardized profile measurement standard, and is controlled by the standard deviation of the depth distribution of the profile. When the standard deviation fluctuation is greater than 5%, the interval threshold is linearly adjusted with the change of the standard deviation. In this example, the standard deviation of the profile TOC is 1.45 g / kg, which is in the normal section, and the interval maintains the standard value. This result indicates being within the interval, there is a pollution gradient in the profile, and the subsequent pollution block division stage can be entered based on this.

[0023] The block boundary delimitation sub-module, based on the pollution intensity division basis, judges whether the pollution gradients of consecutive horizons continuously meet the zoning threshold conditions, delimits the upper and lower boundary positions of the block, numbers each delimited block, and generates pollution intensity block division information; Based on the pollution intensity division basis, judge whether the pollution gradients of each consecutive horizon continuously meet the zoning threshold conditions, and set the continuous determination threshold of the pollution gradient to , this value refers to the lateral profile uniformity standard of the profile, and the recommended value of 0.15 - 0.25 g / kg is used as the section determination band in the soil pollution detection standard. The middle value of 0.2 g / kg is taken in this section. Calculate the absolute value of the difference between adjacent horizons, judge whether it is lower than this threshold, and screen out the continuous segments. For example, the differences between the 1st - 3rd layers are -0.32, -1.05, -0.25 g / kg. The 1st - 2nd layers do not meet the requirements, and the 2nd - 3rd layers meet the requirements. Continue to judge downward. The difference between the 3rd - 4th layers of -0.25 meets the requirements. Judge in this way to form continuous segment blocks, set the block boundaries as the upper and lower boundary horizons of the continuous segments, number the blocks, and the block numbers are {Block 1: 20 - 50 cm, Block 2: 50 - 100 cm}. This threshold of 0.2 g / kg is a fixed value and does not change with the number of sampling layers, but can be dynamically adjusted according to the differences in the physical properties of the profile. If the overall coefficient of variation of the profile > 0.25, the threshold is recommended to be adjusted to 0.25 g / kg to improve the zoning sensitivity.

[0024] Establish pollution intensity block division information, as shown in Table 1.

[0025] Table 1 Information Table for Pollutant Intensity Block Division

[0026] As shown in Table 1, two pollutant intensity blocks are divided in this section. The result shows that there are obvious regional differences in pollutant intensity in the section, and the information of each subsequent block is used to guide the division of oxidant dosing areas.

[0027] The oxidant block usage inference module includes: Based on the pollutant intensity block division information, the block parameter acquisition sub-module collects the water content, particle size distribution, total organic carbon content, and effective porosity of each block, classifies and organizes them according to the block index, and establishes a block physical parameter set; Obtain the pollutant intensity block division information, collect the water content, particle size distribution, total organic carbon content, and effective porosity of each block, classify and organize them according to the block index. First, read the pollutant intensity block division information to determine the total number of blocks , assuming that a total of 3 blocks, namely Block 1, Block 2, and Block 3, are divided in the example, and each block corresponds to a different section depth range. For each block, a water content sensor is arranged on-site for real-time collection, and the particle size distribution is detected using a laser particle size analyzer. The normalized value of the uniformity coefficient of the block particle size distribution is calculated as D60 / D10 and the normalized value is taken. The high-temperature combustion method is used for detection, and the block mean value is normalized. The effective porosity is calculated by combining the permeability method and the specific gravity method. The physical parameters of each block are sorted out and classified according to the block number to generate a block physical parameter set, as shown in Table 2 for example; Table 2 Block Physical Parameter Table

[0028] As shown in Table 2, a block physical parameter set is formed for use in subsequent determination steps.

[0029] Based on the block physical parameter set, the oxidant permeability determination sub-module calls the water content, particle size distribution, and effective porosity to determine whether the oxidant permeability coefficient is greater than the set permeability threshold, screens the blocks that meet the conditions, and obtains the oxidant permeability screening result; Based on the block physical parameter set, the water content, particle size distribution, and effective porosity are called to determine whether the oxidant permeability coefficient is greater than the set permeability threshold. First, the permeability coefficient is calculated using the measured porosity and particle size distribution, and the formula is used, where is an empirical coefficient, which can be set to 0.15, is the effective porosity percentage / 100. The permeability coefficient of each block is calculated in turn. For example, for Block 1: , the seepage threshold is set to 0.040, which is set according to the recommended value in the in-situ chemical oxidation construction guide for groundwater contaminated sites and is adjusted according to the average particle size of the soil layer. When D50 > 0.15mm, it takes 0.04cm / s. The permeability coefficient of Block 1 is greater than the threshold, so it is determined to meet the conditions. Blocks 2 and 3 are calculated accordingly, and the blocks that meet the requirements are screened out to establish the screening results of the oxidant seepage capacity.

[0030] According to the screening results of the oxidant seepage capacity, the total organic carbon content and the effective reaction interface area of the corresponding block are called by the oxidant demand calculation sub-module, and the formula is used: ; The oxidant demand of each block is obtained through operation, and the oxidant demand of each block is summarized to establish the distribution information of the oxidant usage amount of each block; Among them, represents the oxidant demand of the block, represents the normalized value of the total organic carbon content of the th block, represents the normalized value of the particle size distribution uniformity coefficient of the th block, represents the , represents the normalized value of the effective seepage volume of the oxidant, indicating the ratio of the actual permeable volume of the oxidant in the block to the standard reference volume, represents the normalized value of the overall volume of the block, indicating the ratio of the total volume of the block to the standard reference volume, represents the total number of blocks, which depends on the number of blocks in the pollution intensity block division information; According to the screening results of the oxidant seepage capacity, the and the effective reaction interface area of the corresponding block are called, and the formula is used: ; Among them, represents the oxidant demand of the block, is the normalized value of the total organic carbon content of the th block, is the normalized value of the particle size distribution uniformity coefficient of the th block, is the stability constant and can be set to , is the normalized value of the effective seepage volume of the oxidant, is the normalized value of the overall volume of the block, is the total number of blocks. Summation is used to accumulate the contributions of the oxidant requirements of each block, and multiplication reflects the combined effect of organic carbon load and soil layer structure, division balances the reaction interface effect, and square root processing smooths the contribution of the reaction area, corrects the amplification effect of the oxidant penetration volume on the demand. The advantage of the formula is that through , , these three dynamic terms, the comprehensive response ability of the demand calculation to soil structure heterogeneity, oxidant migration effectiveness, and reaction contact efficiency is enhanced.

[0031] Example calculation, parameters of block 1 , , , , calculate: ; Multiply the whole by , the demand of block 1 in units. Calculate block 2 and block 3 in turn, and summarize the total of the three blocks is 4.7289 units. This result shows that the oxidant demand can dynamically reflect the differences in organic carbon load and oxidant migration efficiency of different blocks through this calculation method, and form the distribution information of the oxidant usage amount of each block.

[0032] The dynamic regulation module for oxidant dosing includes: The degradation rate acquisition sub-module obtains the distribution information of the oxidant usage amount of each block, collects the degradation rate sequences of the main components of the total organic carbon content of each block in the current treatment stage, detects the difference in the component degradation rate between adjacent time points, judges the trend direction of the difference, calculates the average value of the change in the degradation rate of each block, and establishes the trend value of the degradation rate of each block; Obtain the distribution information of the oxidant usage amount of each block, collect the degradation rate sequences of the main components of the total organic carbon content of each block in the current treatment stage. First, read the distribution information of the oxidant usage amount of each block according to the block number, determine the monitoring block range, and set the total number of blocks , block 1, block 2, block 3. The selected time span for collection is from day 0 to day 30 of the treatment stage, the monitoring frequency is once every 5 days, and the collected index is the residual amount of the main components (benzene series, phenols, aliphatic compounds) in the block, calculate the component degradation rate between adjacent time points, and use , where is the TOC content at time point, = 5 days, calculate the benzene series in example block 1 , , then . Determine the trend direction of the difference. If is monotonically decreasing, the trend is positive. If the volatility is greater than the set threshold of 5% (recommended value of fluctuation tolerance according to the soil organic pollutant remediation monitoring guide), it is regarded as an unstable trend. The trend direction determination is added to the block degradation rate trend value, and the average value of the block degradation rate change is calculated. For example, the average degradation rate of three monitoring points in block 1 , and finally the block degradation rate trend value is established.

[0033] The residual rate calculation sub-module, based on the block degradation rate trend value, collects the time series of the oxidant residual concentration in each block, calculates the difference in the residual concentration change between adjacent time points, determines the change rate trend, calculates the average value of the residual concentration change in each block, and establishes the block oxidant residual rate value; According to the block degradation rate trend value, collect the time series of the oxidant residual concentration in each block. First, select monitoring blocks 1, 2, and 3. The collection time period is the same as that for TOC monitoring, from day 0 to day 30, once every 5 days. The oxidant is persulfate, and collect the persulfate concentration in the pore water of the block , calculate the difference in the residual concentration change between adjacent time points , calculate the oxidant residual rate , for example, block 2 , , then . Determine the change rate trend. If it continuously decreases, the trend is positive. If the fluctuation amplitude exceeds 5%, it is regarded as a trend fluctuation. Calculate the average value of the residual concentration change within the time series according to the trend result. For block 2 , and finally establish the block oxidant residual rate value.

[0034] The dosing ratio adjustment sub-module, based on the block oxidant residual rate value, calls the block degradation rate trend value, compares the change trend of the oxidant residual rate with the degradation rate trend of the total organic carbon content, and uses the formula: ; Calculate to obtain the oxidant dosing adjustment coefficient for each block. Compare the oxidant dosing adjustment coefficient with the adjustment threshold. If the dosing adjustment coefficient is greater than the adjustment threshold, increase the oxidant dosing ratio. Otherwise, maintain the current oxidant dosing ratio and establish the block oxidant dosing adjustment information; Among them, represents the oxidant dosing adjustment coefficient of the block, represents the oxidant residual rate of the block, represents the degradation rate of the total organic carbon content of the block, represents a small positive value, usually taking the value of , represents the normalized value of the oxidant dosing duration in the represents the normalized value of the corresponding degradation stage duration in the represents the normalized value of the oxidant dosing adjustment decision threshold, represents the block number.

[0035] Based on the block oxidant residue rate value, call the block degradation rate trend value, compare the change trend of the oxidant residue rate with the degradation rate trend of the total organic carbon content, and use the formula: ; Among them, is the oxidant dosing adjustment coefficient for the block, is the oxidant residue rate for the block, is the degradation rate of the total organic carbon content for the block, is a small positive value, which can be set to , is the normalized value of the oxidant dosing duration in the block, is the normalized value of the corresponding degradation stage duration in the block, is the normalized value of the oxidant dosing adjustment decision threshold. The subtraction in the formula is used to measure the rate mismatch between the two, and adding to the denominator prevents division by zero, and the square root term measures the amplification effect of the dosing persistence difference on the adjustment coefficient. The advantage of the formula is that by dynamically integrating the time normalization term into the adjustment calculation, the adjustment strategy not only considers the rate difference but also takes into account the progress difference of different block treatment stages, enhancing the sensitivity of the adjustment decision.

[0036] Example calculation, for block 3 , , , , calculate: ; Threshold , this value is optimized according to the suggestions in the dynamic optimization guidelines for the chemical oxidation process. The typical threshold is in the range of 0.3 - 0.4, and it is optimized on-site according to the economy of the oxidant and the degradation matching. Currently, taking 0.35 corresponds to the operating condition balance point. This result shows <s , currently maintain the original oxidant dosing ratio without adjustment, and establish the information on the adjustment of oxidant dosing in the block.

[0037] The dynamic monitoring module for the oxidation reaction process includes: Based on the information on the adjustment of oxidant dosing in the block, the oxidant residue trend extraction sub-module collects the measured values of the oxidant residue concentration at multiple depths in the soil, calculates the difference in oxidant concentration between adjacent depth layers, determines the change direction of the concentration difference, calculates the diffusion rate of the oxidant in each layer within the time series, screens the change rate of the oxidant concentration diffusion trend within the block, and establishes the information on the change rate of the oxidant diffusion trend. Based on the information on the adjustment of oxidant dosing in the block, collect the measured values of the oxidant residue concentration at multiple depths in the soil. First, determine the block numbers to be monitored according to the adjustment information. The example monitoring blocks are 1, 2, and 3. The monitoring depth range of the soil profile is set to 0 - 100 cm, and sampling layers are set every 10 cm, with a total of 10 layers vertically. The residual concentration of persulfate in the pore water is collected for each layer. , The selected sampling time period is from the 0th day to the 30th day of the treatment stage, and sampling is carried out every 5 days. A concentration time-depth matrix is constructed, and the concentration difference between adjacent depth layers is calculated. , for example, the concentration of the 10 - 20 cm layer in block 2. , , then . Determine the change direction of the difference. If remaining negative indicates that the diffusion trend is normal. If it alternates between positive and negative and the fluctuation is greater than 10%, it is determined that the diffusion trend is unstable. This threshold of 10% is set according to the in-situ chemical oxidation engineering monitoring specification and reflects the reference standard for the vertical diffusion uniformity of the profile. Calculate the diffusion rate of the oxidant in each layer within the time series. , for example, a certain layer in block 2. . Screen out the representative sections in each block and calculate the change rate of the oxidant concentration diffusion trend. , using the average absolute value trend calculation, establish the information on the change rate of the oxidant diffusion trend.

[0038] Based on the information on the change rate of the oxidant diffusion trend, the total organic carbon concentration trend extraction sub-module collects the change values of the remaining concentration of the total organic carbon content in the corresponding soil layer, calculates the difference in the decrease of the organic carbon concentration between adjacent time periods, determines the stability of the decreasing trend, calculates the decreasing rate of the total organic carbon concentration within the time series, and establishes the information on the decreasing rate of the organic carbon concentration.

[0039] According to the information on the change rate of the oxidant diffusion trend, collect the change values of the remaining concentration of the total organic carbon content in the corresponding soil layer. The monitoring profile is stratified corresponding to the oxidant diffusion monitoring layer. The total organic carbon is monitored for each layer. , The sampling time is synchronized with the oxidant. From the 0th day to the 30th day, sampling is carried out every 5 days. Calculate the , for example, a certain layer in block 3. 、 If . To judge the stability of the downward trend, if it shows continuous decline and the fluctuation of the decline difference between adjacent stages is < 5%, the trend is determined to be stable. If it is ≥ 5%, it is an unstable trend. This 5% threshold is recommended according to the evaluation standard for the remediation effect of organically contaminated soil and reflects the standard for the stability of TOC changes in the profile. Calculate the decline rate of the total organic carbon concentration within the time series , for a certain layer in Example Block 3 . Integrate the data of each layer to establish the information on the decline rate of organic carbon concentration at the block level.

[0040] The oxidation reaction state determination sub-module, based on the information on the decline rate of organic carbon concentration, calls the change rate of the oxidant diffusion trend, compares the trend difference between the two, calculates the change factor of the oxidation reaction state, determines whether the change factor exceeds the diffusion-limited determination threshold or the degradation stagnation determination threshold, obtains the oxidation reaction stage state of each block, and establishes the dynamic state information of the block oxidation reaction; Based on the information on the decline rate of organic carbon concentration, call the change rate of the oxidant diffusion trend, compare the trend difference between the two, calculate the change factor of the oxidation reaction state, and use the formula: ; Wherein, is the change factor of the oxidation reaction state of the th block, is the change rate of the oxidant diffusion trend of the th block, is the decline rate of the organic carbon concentration of the th block, is the stability constant. The absolute value of the subtraction in the formula quantifies the difference between the oxidant diffusion trend and the organic carbon degradation trend. Adding to the denominator prevents division by zero and normalizes the degree of difference. The advantage of the formula is that it can dynamically capture the changes in the oxidation reaction kinetic state through trend comparison, avoid the lag risk brought by static rate judgment, and improve the sensitivity of reaction control.

[0041] Example calculation for Block 1 、 , calculate: ; Diffusion-limited determination threshold and degradation stagnation determination threshold . These thresholds are recommended according to the technical guidelines for the regulation of in-situ chemical oxidation processes and usually take values in the range of 0.05 - 0.1. They can be optimized on-site according to the diffusion-degradation coupling performance. In this example, select 、 . This result shows that Slightly higher than the degradation stagnation threshold, it is determined that the current block is in the critical state of slow degradation, and the dynamic state information of the block oxidation reaction is obtained.

[0042] The oxidant dosing stop determination module includes: Based on the dynamic state information of the block oxidation reaction, the diffusion ability determination sub-module obtains the current oxidant residual concentration gradient of each block, collects the change value of the oxidant diffusion rate of the block, calculates the difference in diffusion rate between adjacent time periods, judges the change trend of the diffusion ability, compares the trend change rate with the diffusion ability threshold, screens out the blocks with insufficient diffusion ability, and establishes a list of blocks with insufficient diffusion ability; Based on the dynamic state information of the block oxidation reaction, the current persulfate residual concentration gradient of each block is obtained. First, read the persulfate concentrations at different depths in the cross-section of each block point , the monitoring time period is from day 0 to day 30, and it is collected every 5 days. The depth range is 0 - 100 cm, with 10 cm per layer and a total of 10 layers, to construct a cross-section residual concentration matrix. Calculate the concentration gradient of each layer , for example, in block 1, the 20 - 30 cm layer 、 , then . Collect the change value of the oxidant diffusion rate of each block , = 5 days, example . Calculate the difference in diffusion rate between adjacent time periods , judge the change trend of the diffusion ability. If is continuously negative more than 3 times, or the fluctuation amplitude exceeds the set threshold of 2`0%, it is determined that the diffusion ability is insufficient. This 20% threshold is recommended according to the in-situ oxidation reaction process monitoring guidelines and refers to the cross-section diffusion ability fluctuation tolerance. Screen out the blocks that meet the conditions and establish a list of blocks with insufficient diffusion ability. For example, block 1 and block 3 are included in the list, and block 2 remains normal.

[0043] The organic carbon degradation state determination sub-module collects the change value of the remaining concentration of the total organic carbon in the corresponding block according to the list of blocks with insufficient diffusion ability, calculates the degradation rate of the total organic carbon content in adjacent time periods, judges whether the degradation rate is lower than the total organic carbon content degradation rate threshold, screens out the blocks that meet the degradation stagnation state, and establishes a list of blocks in the degradation stagnation state; According to the list of blocks with insufficient diffusion ability, collect the change value of the remaining concentration of the total organic carbon in the corresponding block , the monitoring period is synchronized with the diffusion monitoring, from day 0 to day 30, and it is collected every 5 days. Calculate the adjacent time periods , calculate the degradation rate , for example, block 1 、 , ]. Determine whether the degradation rate is lower than the total organic carbon content degradation rate threshold , threshold , set according to the recommended interval (2.0 - 3.0) of the mild reaction retardation threshold in the organic pollutant remediation assessment standard for contaminated soil, and select the median value. If , filter into the list of blocks in the degradation stagnation state. Example block 1 meets the criteria, block 3 meets the criteria, block 2 is not included in the list, and a list of blocks in the degradation stagnation state is established.

[0044] The dosing stop strategy formulation sub-module is based on the list of blocks in the degradation stagnation state, calls the list of blocks with insufficient diffusion ability, compares the block numbers in the two lists, filters the blocks that simultaneously meet the conditions of insufficient diffusion ability and degradation stagnation state, assigns the dosing stop time points, generates the dosing stop order of the blocks, and establishes the dosing stop instruction set for the blocks of the oxidant; Based on the list of blocks in the degradation stagnation state, call the list of blocks with insufficient diffusion ability, compare the block numbers in the two lists, and filter the blocks that simultaneously meet the conditions of insufficient diffusion ability and degradation stagnation state. Example blocks 1 and 3 meet the conditions at the same time. To ensure the orderly process of stopping dosing, assign the dosing stop time points, set a minimum interval of 5 days strategy to prevent excessive disturbance of the site flow field caused by large-area sudden stops at the same time. According to the current dosing cycle of the oxidant , according to , for example, if block 1 is currently on the 30th day, then days, block 3 days, generate the dosing stop order for each block. Establish the dosing stop instruction set for the blocks of the oxidant, as shown in Table 3.

[0045] Table 3 Dosing stop instruction set for the blocks of the oxidant

[0046] As shown in Table 3, the dosing stop time sequence has been planned. The result shows that blocks 1 and 3 will stop dosing in batches according to the established strategy, and the dosing stop instruction set for the blocks of the oxidant is generated.

[0047] The above is only a preferred embodiment of the present invention, and does not limit the present invention in other forms. Any person skilled in the art may use the disclosed technical content to make changes or modifications into equivalent embodiments with equivalent changes and apply them to other fields. However, any simple modification, equivalent change, and modification made to the above embodiments based on the technical essence of the present invention without departing from the technical solution content of the present invention still fall within the protection scope of the technical solution of the present invention.

Claims

1. The oxidation treatment control system for organic contaminated soil with multi-parameter coupling control is characterized in that The system includes: The pollution distribution gradient zoning module obtains the measured values of the total organic carbon content in multiple sampling horizons in the soil depth, collects the maximum increasing value and the minimum decreasing value in the depth direction, divides the soil into multiple pollution intensity blocks, and generates pollution intensity block division information; The oxidant block usage inference module, based on the pollution intensity block division information, calculates the oxidant demand of each block according to the total organic carbon content of the block, establishes the oxidant dosage of the block, and generates the oxidant usage distribution information of the block; The oxidant dosing dynamic regulation module, based on the oxidant usage distribution information of the block, obtains the time series data of the oxidant residual concentration, determines whether the oxidant consumption rate exceeds the degradation demand of the total organic carbon content, adjusts the oxidant dosing ratio in the next round, and generates the oxidant dosing adjustment information of the block; The oxidation reaction process dynamic monitoring module, based on the oxidant dosing adjustment information of the block, compares the change trend of the oxidant concentration diffusion rate with the decreasing trend of the remaining concentration of the total organic carbon content, detects whether the oxidation reaction in each block enters the diffusion-limited stage or the degradation stagnation stage, and obtains the dynamic state information of the block oxidation reaction.

2. The control system for oxidizing and treating organic contaminated soil with multi-parameter coupling control according to claim 1, wherein The pollution intensity block division information includes block number, block boundary position, and block pollution level. The oxidant usage distribution information of the block includes block number, required oxidant dosage of the block, and the range of oxidant penetration depth of the block. The oxidant dosing adjustment information of the block includes block number, oxidant dosing adjustment ratio, and oxidant dosing adjustment period. The dynamic state information of the block oxidation reaction includes block number, oxidant diffusion trend state, and pollutant degradation trend state.

3. The multi-parameter coupled control organic contaminated soil oxidation treatment control system according to claim 2, wherein The pollution distribution gradient zoning module includes: The depth concentration acquisition sub-module obtains the measured values of the total organic carbon content in multiple sampling horizons in the soil depth, monitors the original values of the total organic carbon content of the sampling horizons, establishes the organic carbon distribution sequence in the depth direction, and generates the soil organic carbon distribution sequence; The pollution gradient determination sub-module, according to the soil organic carbon distribution sequence, calls the maximum increasing value and the minimum decreasing value, calculates and obtains the pollution gradient factor value, determines whether the pollution gradient factor value is within the pollution area division interval, and generates the pollution intensity division basis; The block boundary demarcation sub-module, based on the pollution intensity division basis, determines whether the pollution gradients of consecutive horizons continuously meet the zoning threshold conditions, demarcates the upper and lower boundary positions of the block, numbers each demarcated block, and generates the pollution intensity block division information.

4. The multi-parameter coupled control organic contaminated soil oxidation treatment control system according to claim 3, wherein The oxidant block usage inference module includes: The block parameter acquisition sub-module, based on the pollution intensity block division information, collects the moisture content, particle size distribution, total organic carbon content, and effective porosity of each block, classifies and organizes them according to the block index, and establishes a block physical parameter set; The oxidant penetration ability determination sub-module, based on the block physical parameter set, calls the moisture content, particle size distribution, and effective porosity, determines whether the oxidant penetration coefficient is greater than the set penetration threshold, screens the blocks that meet the conditions, and obtains the oxidant penetration ability screening result; Based on the screening results of the oxidant penetration ability, the oxidant demand calculation sub-module calls the total organic carbon content and the effective reaction interface area of the corresponding block, calculates the oxidant demand for each block, summarizes the oxidant demand for each block, and establishes the distribution information of the oxidant usage amount for each block.

5. The multi-parameter coupled control organic contaminated soil oxidation treatment control system according to claim 4, characterized in that, The oxidant dosing dynamic regulation module includes: The degradation rate acquisition sub-module obtains the distribution information of the oxidant usage amount for each block, collects the degradation rate sequences of the main components of the total organic carbon content in each block at the current treatment stage, detects the difference in the component degradation rate between adjacent time points, judges the trend direction of the difference, calculates the average change in the degradation rate for each block, and establishes the degradation rate trend value for each block; The residual rate calculation sub-module, based on the degradation rate trend value of each block, collects the time series of the oxidant residual concentration in each block, calculates the change difference in the residual concentration between adjacent time points, judges the trend of the change rate, calculates the average change in the residual concentration for each block, and establishes the oxidant residual rate value for each block; The dosing ratio adjustment sub-module, based on the oxidant residual rate value of each block, calls the degradation rate trend value of the block, compares the change trend of the oxidant residual rate with the degradation rate trend of the total organic carbon content, calculates the oxidant dosing adjustment coefficient for each block, compares the oxidant dosing adjustment coefficient with the adjustment threshold. If the dosing adjustment coefficient is greater than the adjustment threshold, increase the oxidant dosing ratio; otherwise, maintain the current oxidant dosing ratio, and establish the oxidant dosing adjustment information for each block.

6. The multi-parameter coupled control organic contaminated soil oxidation treatment control system according to claim 5, characterized in that, The dynamic monitoring module for the oxidation reaction process includes: The oxidant residual trend extraction sub-module, based on the oxidant dosing adjustment information of each block, collects the measured values of the oxidant residual concentration at multiple depths in the soil, calculates the difference in the oxidant concentration between adjacent depth layers, judges the change direction of the concentration difference, calculates the diffusion rate of the oxidant for each layer in the time series, screens the change rate of the oxidant concentration diffusion trend within the block, and establishes the change rate information of the oxidant diffusion trend; The organic carbon concentration trend extraction sub-module, according to the change rate information of the oxidant diffusion trend, collects the change value of the remaining concentration of the total organic carbon content in the corresponding soil layer, calculates the difference in the decrease of the organic carbon concentration between adjacent time periods, judges the stability of the decreasing trend, calculates the decreasing rate of the total organic carbon concentration in the time series, and establishes the decreasing rate information of the organic carbon concentration; The oxidation reaction state determination sub-module, based on the decreasing rate information of the organic carbon concentration, calls the change rate of the oxidant diffusion trend, compares the difference between the two trends, calculates the change factor of the oxidation reaction state, judges whether the change factor exceeds the diffusion limitation determination threshold or the degradation stagnation determination threshold, obtains the state of the oxidation reaction stage for each block, and establishes the dynamic state information of the oxidation reaction for each block.

7. The multi-parameter coupled control organic contaminated soil oxidation treatment control system according to claim 6, characterized in that, The system further includes: The oxidant dosing stop determination module, based on the dynamic state information of the oxidation reaction for each block, detects whether the change in the remaining concentration of the total organic carbon content is lower than the degradation rate threshold of the total organic carbon content. If both conditions are met, it is determined that the dosing for the corresponding block is stopped, the timing strategy for stopping the dosing for the entire block is established, and the oxidant stop dosing instruction set for the block is generated; The oxidant stop dosing instruction set for the block includes the block number, the stop dosing time point, and the stop dosing confirmation flag.

8. The multi-parameter coupled control organic contaminated soil oxidation treatment control system according to claim 7, characterized in that, The oxidant dosing stop determination module includes: The diffusion ability determination sub-module, based on the dynamic state information of the block oxidation reaction, obtains the current oxidant residual concentration gradient of each block, collects the change value of the oxidant diffusion rate in the block, calculates the difference in the diffusion rate in adjacent time periods, judges the change trend of the diffusion ability, compares the trend change rate with the diffusion ability threshold, screens out the blocks with insufficient diffusion ability, and establishes a list of blocks with insufficient diffusion ability; The organic carbon degradation state determination sub-module, according to the list of blocks with insufficient diffusion ability, collects the change value of the remaining concentration of the total organic carbon content in the corresponding block, calculates the degradation rate of the total organic carbon content in adjacent time periods, judges whether the degradation rate is lower than the total organic carbon content degradation rate threshold, screens out the blocks that meet the degradation stagnation state, and establishes a list of blocks in the degradation stagnation state; The dosing stop strategy formulation sub-module, based on the list of blocks in the degradation stagnation state, calls the list of blocks with insufficient diffusion ability, compares the block numbers in the two lists, screens out the blocks that simultaneously meet the conditions of insufficient diffusion ability and degradation stagnation state, assigns the dosing stop time points, generates the block dosing stop sequence, and establishes a block oxidant dosing stop instruction set.

Images