A method for predicting operational risks of wastewater treatment plants

Abstract

This invention discloses a method for predicting operational risks of wastewater treatment devices, relating to the field of wastewater treatment technology. The method includes: establishing a fixed correspondence between aerators and corresponding local aeration monitoring areas, forming a standard gas release reference state; extracting continuous high-value areas of local gas-liquid mixing ratio within the target local aeration monitoring area, and connecting these continuously corresponding high-value areas of local gas-liquid mixing ratio at lateral positions in different vertical layers to form a local gas release channel; then extracting and connecting the high-value positions of local upflow velocity to form a local upflow high-value axis; and identifying the lateral contraction state of the local gas release channel. This invention establishes a fixed correspondence between aerators and corresponding local aeration monitoring areas to form a standard gas release reference state; subsequently, it analyzes the continuous high-value areas of local gas-liquid mixing ratio, the local gas release channel, and the local upflow high-value axis.

Description

Technical Field

[0001] This invention relates to the field of wastewater treatment technology, specifically a method for predicting operational risks of wastewater treatment equipment. Background Technology

[0002] During the operation of a wastewater treatment plant, the aeration unit plays a crucial role in supplying oxygen, promoting water mixing, and maintaining the conditions for biochemical reactions. In actual operation, there is usually a correlation between the gas release state, liquid surface response state, and local oxygen transfer state of the aerator, which together affect the operational stability of the wastewater treatment plant. In order to ensure treatment effect and operational continuity, existing technologies usually monitor aeration intensity, liquid surface disturbance, dissolved oxygen changes, or related operating parameters, and judge the operating status of the plant accordingly. As wastewater treatment plants become more complex in terms of scale and operating conditions, monitoring results of single parameters or local conditions are sometimes insufficient to fully reflect the operational risks caused by local anomalies in aerators. In particular, when local gas release conditions change, the liquid level response and local oxygen transfer conditions may not always maintain a synchronous relationship. Existing technologies have proposed analytical approaches that utilize liquid level fluctuations, dissolved oxygen values, or local flow conditions. These solutions play a positive role in equipment operation monitoring, condition analysis, and anomaly identification. However, from the perspective of predicting the overall operational risks of wastewater treatment devices, if only the gas release state, liquid surface response state, or local oxygen transfer state are obtained separately, without continuous analysis of the correlation between the three, it is difficult to form a bottom-up abnormal chain judgment result in some operating scenarios. Especially when it is necessary to conduct a comprehensive analysis of changes in the lower gas release structure, changes in the liquid surface response, and changes in local oxygen transfer recovery, how to establish a risk prediction method with a continuous correspondence is still a direction that can be further studied and improved in this field. Therefore, this invention proposes a method for predicting the operational risks of wastewater treatment devices. Summary of the Invention

[0003] The purpose of this invention is to provide a method for predicting operational risks of wastewater treatment devices, in order to solve the problems mentioned in the background art.

[0004] This invention can be achieved through the following technical solution: a method for predicting operational risks of wastewater treatment equipment, comprising: Step 1: Establish a fixed correspondence between aerators and corresponding local aeration monitoring areas, and form a standard gas release reference state; Step 2: Extract continuous high-value areas of local gas-liquid mixing ratio within the target local aeration monitoring area, and connect the continuously corresponding continuous high-value areas of local gas-liquid mixing ratio at the lateral positions in different vertical layers to form local gas release channels. Then, extract the high-value positions of local upflow velocity and connect them to form a local upflow high-value axis. Identify the lateral contraction state of the local gas release channels, as well as the continuous correspondence between the local upflow high-value axis and the local gas release channels. When both the lateral contraction state and the continuous corresponding state are satisfied, the gas release morphology shift feature is extracted. Step 3: Extend the local gas release channel vertically into the liquid surface area, extract the liquid surface response center, and extract the positional correspondence between the liquid surface response center and the liquid surface position corresponding to the local gas release channel, the local concentration state of liquid surface disturbance, and the abrupt transition state of the liquid surface boundary to form surface response distortion characteristics. Step 4: Extract the local dissolved oxygen recovery state within the local oxygen transfer response range corresponding to the liquid surface response center, and when the gas release morphology shift characteristics and surface response distortion characteristics have been formed, form the reverse deviation characteristics of oxygen transfer results based on the state that the local dissolved oxygen recovery rate has not increased synchronously and the local dissolved oxygen recovery amplitude is lower than the corresponding level under the standard gas release reference state. Step 5: Verify the abnormal chain closure of the gas release mode shift characteristics, surface response distortion characteristics, and oxygen transfer result reverse deviation characteristics, and output the prediction results of the operation risk of the sewage treatment device.

[0005] A further technical improvement of the present invention is that the method for forming a standard gas release reference state in step one includes: When the corresponding aerator is in the preset stable aeration condition, the local gas-liquid mixing ratio distribution, local upflow velocity distribution, liquid surface response and local dissolved oxygen recovery status of the corresponding local aeration monitoring area are obtained respectively. Among them, the local dissolved oxygen recovery state is the local dissolved oxygen recovery rate and local dissolved oxygen recovery amplitude extracted within the preset recovery time interval, with the moment when the aerator enters the preset stable aeration condition as the recovery starting point; The corresponding distribution and recovery relationships of the local gas-liquid mixing ratio, local upflow velocity, liquid surface response, and local dissolved oxygen recovery within the same local aeration monitoring area are then determined as the standard gas release reference state for that local aeration monitoring area.

[0006] A further technical improvement of the present invention is that the method for forming a local gas release channel in step two includes: The horizontal expansion range of the local gas-liquid mixing ratio of each vertical layer in the local aeration monitoring area under the standard gas release reference state is used as the benchmark gas release space. Under the current working conditions, the local gas-liquid mixing ratio values ​​of multiple lateral acquisition positions in each vertical layer are obtained respectively; in each vertical layer, continuous acquisition positions are selected along the lateral direction that are located inside the reference gas release space and whose local gas-liquid mixing ratio values ​​are higher than the local gas-liquid mixing ratio values ​​of the corresponding lateral acquisition positions under the standard gas release reference state, and the area corresponding to the continuous acquisition positions is determined as the continuous high value area of ​​local gas-liquid mixing ratio. Calculate the lateral distribution range of the region with continuously high local gas-liquid mixing ratio; The region with a continuous high value of local gas-liquid mixing ratio that is smaller in lateral distribution range than the corresponding reference gas release space lateral expansion range is identified as the region in the current vertical layer that participates in the formation of local gas release channels. Extract the lateral center position of the region participating in the formation of local gas release channels in each vertical layer; calculate the offset of the lateral center position between adjacent vertical layers in order from bottom to top; Adjacent areas whose lateral center position offset does not exceed a preset offset are connected sequentially to form a local gas release channel.

[0007] A further technical improvement of the present invention is that the method for forming a local upflow high-value axis in step two includes: Local upflow velocity values ​​were obtained at multiple lateral acquisition locations within each vertical layer. Within each vertical layer, the lateral acquisition position where the local upflow velocity value is higher than the local upflow velocity values ​​of other lateral acquisition positions in the same vertical layer is selected as the high value position of the local upflow velocity in that vertical layer. Then, following the order from bottom to top, the locations of the highest local upflow velocities corresponding to the horizontal positions in adjacent vertical layers are connected sequentially to form a local upflow high value axis.

[0008] A further technical improvement of the present invention lies in: the method for extracting the gas release pattern shift features in step two includes: The lateral distribution range of the regions involved in the formation of local gas release channels in each vertical layer was extracted, and the lateral expansion range of the local gas-liquid mixing ratio of the corresponding vertical layer under the standard gas release reference state was extracted. When the lateral distribution range of the regions participating in the formation of local gas release channels in multiple adjacent vertical layers is smaller than the lateral expansion range of the corresponding local gas-liquid mixing ratio, and the lateral distribution range of the regions participating in the formation of local gas release channels in the upper vertical layer is smaller than the lateral distribution range of the regions participating in the formation of local gas release channels in the lower vertical layer, this state is regarded as the lateral contraction state of the local gas release channels. When the local upwelling high value axis is located within the region participating in the formation of the local gas release channel in multiple adjacent vertical layers, this state is regarded as the continuous correspondence between the local upwelling high value axis and the local gas release channel. When the lateral contraction state of the local gas release channel and the continuous correspondence between the local upflow high value axis and the local gas release channel are simultaneously established, the section formed by the continuous vertical extension of the local upflow high value axis within the region participating in the formation of the local gas release channel is defined as the gas release morphology shift feature.

[0009] A further technical improvement of the present invention is that the method for extracting the liquid surface response center in step three includes: Multiple liquid level acquisition points are set at horizontal intervals within the liquid surface area; the liquid level height change value of each liquid level acquisition point under the current working conditions is continuously acquired, and the liquid level disturbance value of each liquid level acquisition point is calculated based on the liquid level height change value. Under standard gas release reference conditions, the liquid level height change values ​​at each liquid level sampling location are continuously acquired, and the standard liquid level disturbance value at each liquid level sampling location is calculated based on the liquid level height change values. Extend the local gas release channel vertically to the liquid surface area to determine the liquid surface sampling location where the local gas release channel reaches the liquid surface; Starting from the liquid surface sampling position where the local gas release channel reaches the liquid surface, read the liquid surface disturbance value and the standard liquid surface disturbance value of adjacent liquid surface sampling positions one by one along both sides of the horizontal direction, and calculate the liquid surface disturbance offset value of each liquid surface sampling position. Compare the liquid surface disturbance offset values ​​of adjacent liquid surface sampling positions on both sides in the direction away from the starting point; retain the liquid surface sampling positions whose liquid surface disturbance offset values ​​are less than the previous liquid surface sampling position's liquid surface disturbance offset value, and stop further comparison in that direction when the liquid surface disturbance offset value is greater than or equal to the previous liquid surface sampling position's liquid surface disturbance offset value. The area enclosed by the liquid surface sampling positions on both sides of the horizontal axis and the liquid surface sampling position at the starting point is defined as the liquid surface response range; Within the liquid surface response range, the liquid surface sampling location with the largest liquid surface disturbance offset value is selected and determined as the liquid surface response center.

[0010] A further technical improvement of the present invention is that the method for forming surface response distortion characteristics in step three includes: Extend the local gas release channel vertically to the liquid surface area and determine the liquid surface sampling location corresponding to the local gas release channel; Extract the lateral position difference between the liquid surface response center and the liquid surface acquisition position corresponding to the local gas release channel. When the lateral position difference does not exceed the preset position deviation value, this state is taken as the positional correspondence between the liquid surface response center and the liquid surface position corresponding to the local gas release channel. The sum of liquid surface disturbance offset values ​​at each liquid surface acquisition position within the liquid surface response range is calculated. When the ratio of the sum of liquid surface disturbance offset values ​​at each liquid surface acquisition position within the liquid surface response range to the sum of liquid surface disturbance offset values ​​at each liquid surface acquisition position within the liquid surface range is greater than a preset concentration ratio, this state is regarded as a local concentration state of liquid surface disturbance. Extract the difference in liquid surface disturbance offset between the two boundary positions of the liquid surface response range and the adjacent liquid surface acquisition position outside it, and extract the difference in liquid surface disturbance offset between the adjacent liquid surface acquisition positions inside the liquid surface response range. When the difference in liquid surface disturbance offset between the two boundary positions of the liquid surface response range and the adjacent liquid surface acquisition position outside it is greater than the difference in liquid surface disturbance offset between the adjacent liquid surface acquisition positions inside the liquid surface response range, this state is regarded as the abrupt transition state of the liquid surface boundary. When the positional correspondence between the liquid surface response center and the corresponding liquid surface position of the local gas release channel, the local concentration state of liquid surface disturbance, and the abrupt transition state of the liquid surface boundary are all established simultaneously, surface response distortion characteristics are formed.

[0011] A further technical improvement of the present invention is that the method for forming the reverse deviation feature of the oxygen transfer result in step four includes: Multiple dissolved oxygen collection points are set within the local oxygen transfer response range corresponding to the liquid surface response center; the dissolved oxygen values ​​at each collection point are continuously acquired under the current operating conditions. Within the preset recovery time interval, the rate of change and the magnitude of change of dissolved oxygen values ​​at each dissolved oxygen collection location were calculated. Calculate the average rate of change of dissolved oxygen values ​​at each dissolved oxygen sampling location, and determine the average rate as the local dissolved oxygen recovery rate; calculate the average magnitude of change of dissolved oxygen values ​​at each dissolved oxygen sampling location, and determine the average magnitude as the local dissolved oxygen recovery magnitude. Extract the standard local dissolved oxygen recovery rate and standard local dissolved oxygen recovery amplitude within the corresponding local oxygen transfer response range under the standard gas release reference state; When the gas release morphology shift characteristic and surface response distortion characteristic have been formed, and the local dissolved oxygen recovery rate is less than the standard local dissolved oxygen recovery rate, and the local dissolved oxygen recovery amplitude is less than the standard local dissolved oxygen recovery amplitude, the oxygen transfer result reverse deviation characteristic is formed.

[0012] Compared with the prior art, the present invention has the following beneficial effects: This invention establishes a fixed correspondence between aerators and corresponding local aeration monitoring areas to form a standard gas release reference state. Subsequently, it analyzes the continuous high-value areas of local gas-liquid mixing ratio, local gas release channels, and local upflow high-value axes to further identify the lateral contraction state of local gas release channels and the continuous correspondence between local upflow high-value axes and local gas release channels, thereby extracting gas release morphology shift characteristics. This invention can implement changes in gas release structure onto local spatial objects and their continuous changes, providing a unified basis for subsequent liquid surface response analysis and local oxygen transfer analysis. Furthermore, during the liquid surface analysis stage, this invention extends the local gas release channel vertically into the liquid surface area, extracts the liquid surface response center, and further extracts the positional correspondence between the liquid surface response center and the corresponding liquid surface position of the local gas release channel, the local concentration state of liquid surface disturbance, and the abrupt transition state of the liquid surface boundary, thereby forming surface response distortion characteristics. By analyzing the correspondence between the lower gas release structure and the liquid surface response state, the judgment of the liquid surface response is not limited to the result of a single liquid surface fluctuation, but is based on the surface response connected to the local gas release channel, which is beneficial to improving the pertinence of surface anomaly identification. On the other hand, in the local oxygen transfer analysis stage, this invention extracts the local dissolved oxygen recovery state within the local oxygen transfer response range corresponding to the liquid surface response center. When the gas release pattern shift feature and surface response distortion feature have been formed, based on the state that the local dissolved oxygen recovery rate has not increased synchronously and the local dissolved oxygen recovery amplitude is lower than the corresponding level under the standard gas release reference state, a reverse deviation feature of the oxygen transfer result is formed. Finally, the anomaly chain closure verification of the gas release pattern shift feature, surface response distortion feature and reverse deviation feature of the oxygen transfer result is performed, and the operation risk prediction result of the sewage treatment device is output. It can include the gas release state, liquid surface response state and local oxygen transfer state into the same anomaly chain for continuous judgment, which is conducive to forming a more complete operation risk prediction result. Attached Figure Description

[0013] To facilitate understanding by those skilled in the art, the present invention will be further described below with reference to the accompanying drawings.

[0014] Figure 1 This is a flowchart of the method of the present invention. Detailed Implementation

[0015] To further illustrate the technical means and effects of the present invention in achieving its intended purpose, the following detailed description of the specific implementation methods, structures, features, and effects of the present invention, in conjunction with the accompanying drawings and preferred embodiments, is provided.

[0016] Please see Figure 1 As shown, the present invention provides a method for predicting operational risks of wastewater treatment equipment, including: Step 1: Establish a fixed correspondence between aerators and corresponding local aeration monitoring areas, and form a standard gas release reference state; Specifically, a fixed correspondence is established between aerators and corresponding local aeration monitoring areas, including: In this embodiment, the installation position of each aerator at the bottom of the aeration tank is first determined. Using the installation position of a single aerator as a reference, multiple gas-liquid mixing ratio sampling positions are set horizontally at intervals above it, and multiple sampling layers are set vertically. Then, under the condition that the corresponding aerator continuously outputs gas, the gas-liquid mixing ratio data of each horizontal sampling position within each sampling layer is continuously acquired. The range between the outermost horizontal sampling positions where the gas-liquid mixing ratio value is continuously higher than the preset lower limit value is determined as the coverage boundary range of the aerator's gas release effect in the horizontal direction. Simultaneously, the gas-liquid mixing ratio data of each sampling layer is continuously acquired vertically, and the highest sampling layer where the gas-liquid mixing ratio value is still higher than the preset lower limit value is determined as the propagation boundary position of the aerator's gas release effect in the vertical direction. Then, based on the horizontal coverage boundary range and the vertical propagation boundary position, local aeration monitoring areas corresponding one-to-one with each aerator are delineated. Subsequently, each local aeration monitoring area is numbered, and the number is fixedly bound to the corresponding aerator. Thus, each aerator corresponds to a unique local aeration monitoring area.

[0017] When the corresponding aerator is in a preset stable aeration condition, the local gas-liquid mixing ratio distribution, local upflow velocity distribution, liquid surface response, and local dissolved oxygen recovery status are acquired within the corresponding local aeration monitoring area, including: In this embodiment, the output process of the aerator is continuously recorded, and the gas-liquid mixing ratio distribution data and upflow velocity distribution data of each sampling layer within the corresponding local aeration monitoring area are acquired simultaneously. The range of sampling locations where the gas-liquid mixing ratio value in each sampling layer is continuously higher than the preset lower limit of the gas-liquid mixing ratio is defined as the gas release propagation profile. When the aerator output remains unchanged within a preset stable observation period, and the change in the boundary position of the gas release propagation profile of each sampling layer does not exceed the preset profile change limit within two consecutive sampling cycles, this period is defined as the preset stable aeration condition. Subsequently, under the preset stable aeration condition, multiple sampling layers are set up vertically, and multiple horizontal sampling positions are set up within each sampling layer. The gas-liquid mixing ratio data and upflow velocity data of each horizontal sampling position are acquired to form a local gas-liquid mixing ratio distribution state and a local upflow velocity distribution state; wherein, the upflow velocity is expressed as the ratio of length to time. Simultaneously, multiple liquid level sampling points are set at horizontal intervals within the liquid surface area to continuously acquire the liquid level height change values ​​at each sampling point. The liquid level height change values ​​are expressed in units of length, and the amplitude of the liquid level height change value within a unit sampling time is determined as the liquid surface disturbance value to form the liquid surface response state. Then, multiple dissolved oxygen sampling points are set within the corresponding local aeration monitoring area to continuously acquire the dissolved oxygen value change process at each sampling point to form the raw data of the local dissolved oxygen recovery state. The dissolved oxygen value change rate is expressed as the ratio of concentration to time, and the dissolved oxygen value change amplitude is expressed in units of concentration.

[0018] Taking the moment when the aerator enters the preset stable aeration condition as the recovery starting point, the local dissolved oxygen recovery rate and local dissolved oxygen recovery magnitude are extracted within the preset recovery time interval, including: In this embodiment, when the aerator transitions from an unstable output state to a preset stable aeration state, this transition time is defined as the recovery starting point. From this starting point, dissolved oxygen values ​​at each sampling location are continuously acquired within a preset recovery time interval. Subsequently, the rate of change and magnitude of dissolved oxygen value change at each sampling location within the preset recovery time interval are calculated. Then, the average value of the rate of change of dissolved oxygen value at each sampling location is calculated and determined as the local dissolved oxygen recovery rate. Finally, the average value of the magnitude of change of dissolved oxygen value at each sampling location is calculated and determined as the local dissolved oxygen recovery magnitude. This method ensures that the statistical caliber of the local dissolved oxygen recovery rate and local dissolved oxygen recovery magnitude is consistent with the average value caliber used in subsequent claims.

[0019] The corresponding distribution and recovery relationships of the local gas-liquid mixing ratio, local upflow velocity, liquid surface response, and local dissolved oxygen recovery within the same local aeration monitoring area are determined as the standard gas release reference state, including: In this embodiment, firstly, within the same local aeration monitoring area, the local gas-liquid mixing ratio data and upflow velocity data are paired point-by-point according to the same sampling layer and the same lateral sampling position, forming a corresponding distribution relationship between the local gas-liquid mixing ratio distribution and the local upflow velocity distribution. Then, the paired data within each sampling layer are extended vertically to the liquid surface area, and paired with the liquid surface disturbance values ​​at the liquid surface sampling position in a one-to-one lateral correspondence manner, forming a corresponding distribution relationship between the lower gas release propagation distribution and the liquid surface response distribution. Subsequently, the corresponding distribution relationship is uniformly paired with the local dissolved oxygen recovery rate and local dissolved oxygen recovery amplitude already extracted within the same local aeration monitoring area, forming a corresponding recovery relationship. Finally, the corresponding distribution relationship and the corresponding recovery relationship are jointly determined as the standard gas release reference state for this local aeration monitoring area.

[0020] Step 2: Extract continuous high-value areas of local gas-liquid mixing ratio within the target local aeration monitoring area, and connect the continuously corresponding continuous high-value areas of local gas-liquid mixing ratio at the lateral positions in different vertical layers to form local gas release channels. Then, extract the high-value positions of local upflow velocity and connect them to form a local upflow high-value axis. Identify the lateral contraction state of the local gas release channels, as well as the continuous correspondence between the local upflow high-value axis and the local gas release channels. When both the lateral contraction state and the continuous corresponding state are satisfied, the gas release morphology shift feature is extracted. Specifically, the horizontal expansion range of the local gas-liquid mixing ratio at each vertical layer of the local aeration monitoring area under the standard gas release reference state is used as the benchmark gas release space. Under the current operating conditions, continuous high-value areas of the local gas-liquid mixing ratio in each vertical layer are extracted, including: In this embodiment, the standard gas release reference state established in step one is first retrieved, and the lateral expansion range of the local gas-liquid mixing ratio for each vertical layer in the corresponding local aeration monitoring area is extracted. The lateral expansion range in each vertical layer is determined as the reference gas release space for that layer. Subsequently, under the current operating conditions, multiple vertical layers are set up along the vertical direction, and multiple lateral acquisition positions are set up along the lateral direction within each vertical layer. The local gas-liquid mixing ratio value at each lateral acquisition position is acquired point by point. Then, within each vertical layer, the local gas-liquid mixing ratio value at each lateral acquisition position is compared sequentially with the local gas-liquid mixing ratio value at the same lateral acquisition position under the standard gas release reference state. Continuous acquisition positions located within the reference gas release space and with local gas-liquid mixing ratio values ​​higher than those at the corresponding lateral acquisition positions under the standard gas release reference state are retained, and the area corresponding to these continuous acquisition positions is determined as the continuous high-value area of ​​the local gas-liquid mixing ratio. In this way, the continuous high-value areas of local gas-liquid mixing ratio in each vertical layer are uniformly established based on the comparison with the standard gas release reference state, rather than based on the judgment of the absolute high value of a single layer.

[0021] Based on the continuous high value region of local gas-liquid mixing ratio, a local gas release channel is formed, and a local upflow high value axis is formed, including: In this embodiment, the lateral distribution range of the continuous high-value region of local gas-liquid mixing ratio within each vertical stratum is first calculated. Regions with a lateral distribution range smaller than the lateral expansion range of the corresponding reference gas release space are identified as areas participating in the formation of local gas release channels within the current vertical stratum. Next, the lateral center position of the region participating in the formation of local gas release channels in each vertical stratum is extracted. The offset of the lateral center position between adjacent vertical strata is calculated in a bottom-up order. Adjacent regions with lateral center position offsets not exceeding a preset offset are sequentially connected to form local gas release channels. Subsequently, under the same vertical stratum and lateral acquisition position arrangement, the local upflow velocity values ​​at each lateral acquisition position are obtained. Within each vertical stratum, the lateral acquisition position with a local upflow velocity value higher than the local upflow velocity values ​​at other lateral acquisition positions within the same vertical stratum is identified as the high-value position of the local upflow velocity for that vertical stratum. Then, in a bottom-up order, the corresponding high-value positions of the lateral positions in adjacent vertical strata are sequentially connected to form a local upflow high-value axis. This step allows for the formation of anomalous structural objects that can extend continuously in the vertical direction, based on both the local gas-liquid mixing ratio distribution and the local upflow velocity distribution.

[0022] Identifying the lateral contraction state of local gas release channels and the continuous correspondence between the local upflow high-value axis and the local gas release channels, including: In this embodiment, the lateral distribution range of the regions involved in the formation of local gas release channels in each vertical stratum is first extracted, and the lateral expansion range of the local gas-liquid mixing ratio of the corresponding vertical stratum under the standard gas release reference state is also extracted. Then, strata are compared in ascending order. When the lateral distribution range of the regions involved in the formation of local gas release channels in multiple adjacent vertical strata is smaller than the lateral expansion range of the corresponding local gas-liquid mixing ratio, and the lateral distribution range of the regions involved in the formation of local gas release channels in the upper vertical stratum is smaller than the lateral distribution range of the regions involved in the formation of local gas release channels in the lower vertical stratum, this state is defined as a lateral contraction state of the local gas release channels. Next, it is determined whether the local upflow high-value axis is located within the regions involved in the formation of local gas release channels in multiple adjacent vertical strata. When the local upflow high-value axis is located within the regions involved in the formation of local gas release channels in multiple adjacent vertical strata, this state is defined as a continuous correspondence between the local upflow high-value axis and the local gas release channels.

[0023] When both the lateral contraction state of the local release channel and the continuous correspondence between the local upflow high-value axis and the local release channel are simultaneously established, the release morphology shift feature is extracted: In this embodiment, when both of the above states are simultaneously established, a segment formed by the continuous vertical extension of the local upflow high-value axis within the region participating in the formation of the local gas release channel is extracted, and this segment is identified as the gas release pattern shift feature. Specifically, firstly, the position segment where the local upflow high-value axis falls within the region participating in the formation of the local gas release channel is determined in multiple adjacent vertical layers. Then, the position segments that continuously fall within the region in each vertical layer are sequentially connected along the vertical direction to form a continuous vertical segment. This continuous vertical segment is output as the gas release pattern shift feature.

[0024] Step 3: Extend the local gas release channel vertically into the liquid surface area, extract the liquid surface response center, and extract the positional correspondence between the liquid surface response center and the liquid surface position corresponding to the local gas release channel, the local concentration state of liquid surface disturbance, and the abrupt transition state of the liquid surface boundary to form surface response distortion characteristics. Specifically, the liquid surface response center is extracted around the receiving position of the local gas release channel on the liquid surface, including: In this embodiment, multiple liquid surface sampling points are first set at lateral intervals within the liquid surface area. The lateral spacing between each sampling point is determined according to the ratio of the lateral width of the target local aeration monitoring area on the liquid surface to the preset number of sampling points. Subsequently, the liquid surface height change value of each sampling point is continuously acquired at a preset sampling period, and the liquid surface height change value is expressed in length units. Within each sampling period, the absolute value of the difference between the liquid surface height value at the current sampling time and the liquid surface height value at the previous sampling time is taken as the liquid surface disturbance value of that sampling point within that sampling period. At the same time, under the standard gas release reference state, the liquid surface height change value of each sampling point is continuously acquired according to the same liquid surface sampling point arrangement and the same preset sampling period, and the standard liquid surface disturbance value is obtained using the same calculation caliber.

[0025] Next, the localized gas release channel formed in step two is extended vertically into the liquid surface area. The liquid surface sampling position where the localized gas release channel reaches the liquid surface is determined. Starting from this sampling position, the liquid surface disturbance values ​​and standard liquid surface disturbance values ​​of adjacent sampling positions are extracted one by one along both sides. The absolute value of the difference between the liquid surface disturbance value and the standard liquid surface disturbance value is determined as the liquid surface disturbance offset value. Subsequently, the liquid surface disturbance offset values ​​of adjacent sampling positions on both sides are compared in the direction away from the starting point. When the liquid surface disturbance offset value of the subsequent sampling position is less than that of the previous sampling position, the sampling position is retained. When the liquid surface disturbance offset value of the subsequent sampling position is greater than or equal to that of the previous sampling position, the extraction in that direction is stopped. The liquid surface response range is defined by the liquid surface sampling positions retained on both sides of the horizontal axis and the liquid surface sampling position at the starting point. Within the liquid surface response range, the liquid surface sampling position with the largest liquid surface disturbance offset value is selected as the liquid surface response center. Among them, the outermost liquid surface sampling position among the liquid surface sampling positions retained on both sides of the horizontal axis is determined as the left boundary position and the right boundary position of the liquid surface response range, respectively.

[0026] Extract the positional correspondence between the liquid surface response center and the corresponding liquid surface position in the local gas release channel, and extract the local concentrated state of liquid surface disturbance, including: In this embodiment, based on the already determined liquid surface response center, the local gas release channel is extended vertically into the liquid surface area to determine the corresponding liquid surface sampling position. Then, the lateral position difference between the liquid surface response center and the liquid surface sampling position corresponding to the local gas release channel is calculated. This lateral position difference is expressed in units of length and is specifically calculated as the absolute value of the difference between their lateral coordinates. When the lateral position difference does not exceed a preset position deviation value, it is determined that the liquid surface response center and the liquid surface position corresponding to the local gas release channel have a positional correspondence. The preset position deviation value is determined based on the statistical results of the lateral position difference between the liquid surface response center and the corresponding liquid surface position of the local gas release channel under the standard gas release reference state. Specifically, it can be the average value of the statistical results plus a preset correction value, or the maximum allowable value of the statistical results.

[0027] Next, the sum of liquid surface disturbance offset values ​​at each liquid surface sampling location within the liquid surface response range is calculated, and the sum of liquid surface disturbance offset values ​​at all liquid surface sampling locations within the liquid surface range is also calculated. The former is then compared to the latter, and the ratio is a dimensionless proportion. When this proportion is greater than a preset concentration proportion, it is determined that the liquid surface disturbance has formed a local concentration state within the liquid surface response range. The preset concentration proportion is determined based on the statistical results of the proportion of the liquid surface disturbance offset value within the liquid surface response range to the overall liquid surface disturbance offset value under the standard gas release reference state. Specifically, it can be the average value of the statistical results plus a preset correction proportion, or the upper limit of the statistical results.

[0028] Extracting the abrupt transition state at the liquid surface boundary, including: In this embodiment, based on the determined liquid surface response range, the outermost liquid surface sampling positions on the left and right sides of the liquid surface response range are determined respectively, and these are taken as the left and right boundary positions of the liquid surface response range, respectively. Then, the liquid surface disturbance offset difference between the left boundary position and its adjacent left-outer liquid surface sampling position, and the liquid surface disturbance offset difference between the right boundary position and its adjacent right-outer liquid surface sampling position are calculated. Simultaneously, the liquid surface disturbance offset differences between all adjacent liquid surface sampling positions within the liquid surface response range are extracted, and the average value of all adjacent liquid surface sampling differences is calculated. This average value is used as the comparison caliber for the liquid surface disturbance offset differences between adjacent liquid surface sampling positions within the liquid surface response range. Then, the difference in liquid surface disturbance offset between the left boundary position of the liquid surface response range and its adjacent left outer liquid surface sampling position, and the difference in liquid surface disturbance offset between the right boundary position of the liquid surface response range and its adjacent right outer liquid surface sampling position are compared with the average value respectively. When the difference in liquid surface disturbance offset between the two boundary positions and their adjacent outer liquid surface sampling positions is greater than the average value, the liquid surface boundary abrupt transition state is determined to be established.

[0029] When the positional correspondence, the localized concentration of liquid surface disturbance, and the abrupt transition of the liquid surface boundary all occur simultaneously, surface response distortion characteristics are formed, including: In this embodiment, when the positional correspondence between the liquid surface response center and the liquid surface position corresponding to the local gas release channel is established, the local concentration state of liquid surface disturbance is established, and the abrupt transition state of the liquid surface boundary is established, the judgment results of the above three states are jointly verified; when the three states simultaneously meet the conditions, the corresponding liquid surface response result is determined as the surface response distortion feature. Specifically, the comparison results of the lateral position difference, the calculation results of the proportion of liquid surface disturbance offset value, and the comparison results between the difference of liquid surface boundary position offset value and the difference of the average offset value within the liquid surface response range are first uniformly summarized, and then the three results are judged simultaneously; when all three results meet the requirements, the surface response distortion feature is output.

[0030] Step 4: Extract the local dissolved oxygen recovery state within the local oxygen transfer response range corresponding to the liquid surface response center, and when the gas release morphology shift characteristics and surface response distortion characteristics have been formed, form the reverse deviation characteristics of oxygen transfer results based on the state that the local dissolved oxygen recovery rate has not increased synchronously and the local dissolved oxygen recovery amplitude is lower than the corresponding level under the standard gas release reference state. Specifically, the local oxygen transfer response range is determined around the liquid surface response center, and multiple dissolved oxygen collection points are deployed, including: In this embodiment, the liquid surface response center already formed in step three is first invoked, and its vertical projection is used as the center of the local oxygen transfer response range. Then, using the center as a reference, the horizontal width of the liquid surface response range is taken as the horizontal boundary of the local oxygen transfer response range, and the vertical boundary of the water body between a preset depth below the liquid surface and a preset depth below the liquid surface response center is taken as the vertical boundary of the local oxygen transfer response range. The preset depth is determined by the vertical height of the target local aeration monitoring area and the width of the liquid surface response range already formed in step three, ensuring that the local oxygen transfer response range covers the main oxygen transfer change area below the liquid surface response center. Next, multiple dissolved oxygen sampling points are arranged at intervals along the horizontal and vertical sides within the local oxygen transfer response range. The spacing between adjacent horizontally adjacent dissolved oxygen sampling points is determined by the ratio of the horizontal boundary of the local oxygen transfer response range to the preset number of horizontal sampling points, and the spacing between adjacent vertically adjacent dissolved oxygen sampling points is determined by the ratio of the vertical boundary of the local oxygen transfer response range to the preset number of vertical sampling points. Through this step, the local oxygen transfer response range is defined as the local water body analysis range extending downwards from the liquid surface response center, rather than an abstract affected area.

[0031] Continuously acquire dissolved oxygen values ​​at each sampling location under the current operating conditions, and extract the local dissolved oxygen recovery rate and local dissolved oxygen recovery magnitude, including: In this embodiment, dissolved oxygen values ​​at each sampling location are first continuously acquired at a preset sampling period. The preset sampling period is determined based on the frequency of dissolved oxygen changes within the target local aeration monitoring area, specifically taking a period less than one-fifth to one-tenth of a complete dissolved oxygen fluctuation cycle within the local oxygen transfer response range corresponding to the liquid surface response center. Subsequently, within a preset recovery time interval, the rate of change and amplitude of dissolved oxygen value changes at each sampling location are calculated. The preset recovery time interval is defined as the moment when the corresponding aerator enters the preset stable aeration condition under the current operating condition, and is determined according to the same time length as the standard gas release reference state in step one. Then, the average value of the rate of change of dissolved oxygen value at each sampling location is calculated, and this average value is determined as the local dissolved oxygen recovery rate. Next, the average value of the amplitude of change of dissolved oxygen value at each sampling location is calculated, and this average value is determined as the local dissolved oxygen recovery amplitude. Thus, the recovery results of multiple dissolved oxygen sampling locations within the local oxygen transfer response range are uniformly implemented as two recovery quantities: local dissolved oxygen recovery rate and local dissolved oxygen recovery amplitude.

[0032] Extract the standard local dissolved oxygen recovery rate and standard local dissolved oxygen recovery range under the standard gas release reference state, and complete the corresponding comparisons, including: In this embodiment, the standard gas release reference state established in step one is invoked. Within the same horizontal boundary range, vertical boundary range, preset sampling period, and preset recovery time interval as the current local oxygen transfer response range, the standard local dissolved oxygen recovery rate and standard local dissolved oxygen recovery amplitude under the standard gas release reference state are extracted. Subsequently, the local dissolved oxygen recovery rate extracted under the current operating condition is compared with the standard local dissolved oxygen recovery rate, and the local dissolved oxygen recovery amplitude extracted under the current operating condition is compared with the standard local dissolved oxygen recovery amplitude. When the local dissolved oxygen recovery rate is less than the standard local dissolved oxygen recovery rate, and the local dissolved oxygen recovery amplitude is less than the standard local dissolved oxygen recovery amplitude, it is determined that the local oxygen transfer recovery result shows a weakening recovery relative to the standard gas release reference state. Through this step, the recovery results within the local oxygen transfer response range are uniformly included in the same spatial and temporal scope as the standard gas release reference state for comparison.

[0033] Given that the characteristics of gas release morphology shift and surface response distortion have been established, joint verification is performed in a unified order to generate reverse deviation characteristics of oxygen transfer results, including: In this embodiment, the gas release pattern shift feature formed in step two and the surface response distortion feature formed in step three are first invoked, and it is determined whether both are valid. Once both features are valid, it is then determined whether the local dissolved oxygen recovery rate is less than the standard local dissolved oxygen recovery rate, and whether the local dissolved oxygen recovery amplitude is less than the standard local dissolved oxygen recovery amplitude. Only when the gas release pattern shift feature is valid, the surface response distortion feature is valid, the local dissolved oxygen recovery rate is less than the standard local dissolved oxygen recovery rate, and the local dissolved oxygen recovery amplitude is less than the standard local dissolved oxygen recovery amplitude, is the result determined as a reverse deviation feature of the oxygen transfer result.

[0034] Step 5: Verify the abnormal chain closure of the gas release mode shift characteristics, surface response distortion characteristics, and oxygen transfer result reverse deviation characteristics, and output the operation risk prediction results of the sewage treatment device; Specifically, the release morphology shift characteristics, surface response distortion characteristics, and oxygen transfer result reverse deviation characteristics are uniformly extracted, and the correspondence between abnormal results in the same local aeration monitoring area is established, including: In this embodiment, the gas release morphology shift feature already formed in step two is first invoked to extract its corresponding local aeration monitoring area number, corresponding vertical segment range, and formation time. Then, the surface response distortion feature already formed in step three is invoked to extract its corresponding local aeration monitoring area number, corresponding liquid surface response range, and formation time. Subsequently, the oxygen transfer result reverse deviation feature already formed in step four is invoked to extract its corresponding local aeration monitoring area number, corresponding local oxygen transfer response range, and formation time. Afterward, the three types of feature results are merged according to the local aeration monitoring area number, and the gas release morphology shift feature, surface response distortion feature, and oxygen transfer result reverse deviation feature with the same number are grouped into the same anomaly chain verification unit. In this embodiment, the three types of feature results grouped into the same anomaly chain verification unit also correspond to the same aerator number, ensuring that the anomaly chain verification unit simultaneously possesses local aeration monitoring area consistency and aerator consistency.

[0035] Temporal closure verification was performed on the formation order of the three types of feature results, including: In this embodiment, within the same anomaly chain verification unit, the formation time of the gas release pattern shift feature is first compared with the formation time of the surface response distortion feature. When the formation time of the gas release pattern shift feature is earlier than or equal to the formation time of the surface response distortion feature, the temporal relationship between the two features is retained. Then, the formation time of the surface response distortion feature is compared with the formation time of the oxygen transfer result reverse deviation feature. When the formation time of the surface response distortion feature is earlier than or equal to the formation time of the oxygen transfer result reverse deviation feature, the temporal relationship between the two features is retained. Afterwards, it is determined whether the sequence of "gas release pattern shift feature forms first, surface response distortion feature forms later, and oxygen transfer result reverse deviation feature forms last" is simultaneously satisfied within the same anomaly chain verification unit. When the three satisfy the above sequence relationship, the temporal closure verification is confirmed. Through this step, the formation order between the three types of feature results can be implemented as a temporal closure relationship where downstream results inherit from upstream results, rather than the simultaneous occurrence of three independent anomaly results.

[0036] Positional closure verification is performed on the spatial continuity relationship of the three types of feature results, including: In this embodiment, within the same anomaly chain verification unit, the extension position of the vertical segment corresponding to the gas release pattern offset feature on the liquid surface is first extracted, and the liquid surface response center and liquid surface response range corresponding to the surface response distortion feature are extracted. Then, it is determined whether the liquid surface response center is located within the liquid surface extension range of the vertical segment corresponding to the gas release pattern offset feature, and whether the liquid surface response range covers the liquid surface acquisition position with the largest liquid surface disturbance offset value within the liquid surface extension range. Next, the local oxygen transfer response range corresponding to the oxygen transfer result reverse deviation feature is extracted, and it is determined whether the local oxygen transfer response range is located within the vertical projection range of the liquid surface response center. When both of the aforementioned conditions are met simultaneously, the spatial continuity relationship is established. Through this step, the positional relationship between the three elements—"gas release pattern offset feature—surface response distortion feature—oxygen transfer result reverse deviation feature"—can be established as a continuous continuity relationship from bottom to top, forming a unified spatial chain from the vertical anomaly propagation segment in step two, the liquid surface anomaly response range in step three, and the local oxygen transfer response range in step four.

[0037] The temporal closure verification results and spatial continuity verification results are jointly verified, and the operational risk prediction results of the wastewater treatment device are output, including: In this embodiment, the system first checks whether the temporal closure verification is valid within the same anomaly chain verification unit, and then checks whether the spatial continuity relationship is valid. When both the temporal closure verification and the spatial continuity relationship are valid, the system further checks whether the gas release pattern shift characteristic, surface response distortion characteristic, and oxygen transfer result reverse deviation characteristic are valid at the current verification time. When all three characteristics are valid at the current verification time, the anomaly chain verification unit is determined to be an anomaly chain closure valid, and the corresponding wastewater treatment device operation risk prediction result is output. Specifically, the wastewater treatment device operation risk prediction result includes at least: the corresponding aerator number, the corresponding local aeration monitoring area number, the time when the anomaly chain closure is valid, and the current operation risk status of the local aeration monitoring area.

[0038] The above description is merely a preferred embodiment of the present invention and is not intended to limit the present invention in any way. Although the present invention has been disclosed above with reference to preferred embodiments, it is not intended to limit the present invention. Any person skilled in the art can make some modifications or alterations to the above-disclosed technical content to create equivalent embodiments without departing from the scope of the present invention. Any simple modifications, equivalent changes and alterations made to the above embodiments based on the technical essence of the present invention without departing from the scope of the present invention shall still fall within the scope of the present invention.

Claims

1. A method for predicting operational risks of wastewater treatment plants, characterized in that, include: Step 1: Establish a fixed correspondence between aerators and corresponding local aeration monitoring areas, and form a standard gas release reference state; Step 2: Extract continuous high-value areas of local gas-liquid mixing ratio within the target local aeration monitoring area, and connect the continuously corresponding continuous high-value areas of local gas-liquid mixing ratio at the lateral positions in different vertical layers to form local gas release channels. Then, extract the high-value positions of local upflow velocity and connect them to form a local upflow high-value axis. Identify the lateral contraction state of the local gas release channels, as well as the continuous correspondence between the local upflow high-value axis and the local gas release channels. When both the lateral contraction state and the continuous corresponding state are satisfied, the gas release morphology shift feature is extracted. Step 3: Extend the local gas release channel vertically into the liquid surface area, extract the liquid surface response center, and extract the positional correspondence between the liquid surface response center and the liquid surface position corresponding to the local gas release channel, the local concentration state of liquid surface disturbance, and the abrupt transition state of the liquid surface boundary to form surface response distortion characteristics. Step 4: Extract the local dissolved oxygen recovery state within the local oxygen transfer response range corresponding to the liquid surface response center, and when the gas release morphology shift characteristics and surface response distortion characteristics have been formed, form the reverse deviation characteristics of oxygen transfer results based on the state that the local dissolved oxygen recovery rate has not increased synchronously and the local dissolved oxygen recovery amplitude is lower than the corresponding level under the standard gas release reference state. Step 5: Verify the abnormal chain closure of the gas release mode shift characteristics, surface response distortion characteristics, and oxygen transfer result reverse deviation characteristics, and output the prediction results of the operation risk of the sewage treatment device.

2. The method for predicting operational risks of a wastewater treatment plant according to claim 1, characterized in that, The methods for establishing a standard gas release reference state in step one include: When the corresponding aerator is in the preset stable aeration condition, the local gas-liquid mixing ratio distribution, local upflow velocity distribution, liquid surface response and local dissolved oxygen recovery status of the corresponding local aeration monitoring area are obtained respectively. Among them, the local dissolved oxygen recovery state is the local dissolved oxygen recovery rate and local dissolved oxygen recovery amplitude extracted within the preset recovery time interval, with the moment when the aerator enters the preset stable aeration condition as the recovery starting point; The corresponding distribution and recovery relationships of the local gas-liquid mixing ratio, local upflow velocity, liquid surface response, and local dissolved oxygen recovery within the same local aeration monitoring area are then determined as the standard gas release reference state for that local aeration monitoring area.

3. The method for predicting operational risks of a wastewater treatment plant according to claim 1, characterized in that, The methods for forming localized gas release channels in step two include: The horizontal expansion range of the local gas-liquid mixing ratio of each vertical layer in the local aeration monitoring area under the standard gas release reference state is used as the benchmark gas release space. Under the current working conditions, the local gas-liquid mixing ratio values ​​of multiple lateral acquisition positions in each vertical layer are obtained respectively; in each vertical layer, continuous acquisition positions are selected along the lateral direction that are located inside the reference gas release space and whose local gas-liquid mixing ratio values ​​are higher than the local gas-liquid mixing ratio values ​​of the corresponding lateral acquisition positions under the standard gas release reference state, and the area corresponding to the continuous acquisition positions is determined as the continuous high value area of ​​local gas-liquid mixing ratio. Calculate the lateral distribution range of the region with continuously high local gas-liquid mixing ratio; The region with a continuous high value of local gas-liquid mixing ratio that is smaller in lateral distribution range than the corresponding reference gas release space lateral expansion range is identified as the region in the current vertical layer that participates in the formation of local gas release channels. Extract the lateral center position of the region participating in the formation of local gas release channels in each vertical layer; calculate the offset of the lateral center position between adjacent vertical layers in order from bottom to top; Adjacent areas whose lateral center position offset does not exceed a preset offset are connected sequentially to form a local gas release channel.

4. The method for predicting operational risks of a wastewater treatment plant according to claim 3, characterized in that, Step two involves methods for forming a local upwelling high-value axis, including: Local upflow velocity values ​​were obtained at multiple lateral acquisition locations within each vertical layer. Within each vertical layer, the lateral acquisition position where the local upflow velocity value is higher than the local upflow velocity values ​​of other lateral acquisition positions in the same vertical layer is selected as the high value position of the local upflow velocity in that vertical layer. Then, following the order from bottom to top, the locations of the highest local upflow velocities corresponding to the horizontal positions in adjacent vertical layers are connected sequentially to form a local upflow high value axis.

5. The method for predicting operational risks of a wastewater treatment plant according to claim 4, characterized in that, Step two involves extracting the gas release pattern shift features, including: The lateral distribution range of the regions involved in the formation of local gas release channels in each vertical layer was extracted, and the lateral expansion range of the local gas-liquid mixing ratio of the corresponding vertical layer under the standard gas release reference state was extracted. When the lateral distribution range of the regions participating in the formation of local gas release channels in multiple adjacent vertical layers is smaller than the lateral expansion range of the corresponding local gas-liquid mixing ratio, and the lateral distribution range of the regions participating in the formation of local gas release channels in the upper vertical layer is smaller than the lateral distribution range of the regions participating in the formation of local gas release channels in the lower vertical layer, this state is regarded as the lateral contraction state of the local gas release channels. When the local upwelling high value axis is located within the region participating in the formation of the local gas release channel in multiple adjacent vertical layers, this state is regarded as the continuous correspondence between the local upwelling high value axis and the local gas release channel. When the lateral contraction state of the local gas release channel and the continuous correspondence between the local upflow high value axis and the local gas release channel are simultaneously established, the section formed by the continuous vertical extension of the local upflow high value axis within the region participating in the formation of the local gas release channel is defined as the gas release morphology shift feature.

6. The method for predicting operational risks of a wastewater treatment plant according to claim 1, characterized in that, Step 3 involves extracting the liquid surface response center, including: Multiple liquid level acquisition points are set at horizontal intervals within the liquid surface area; the liquid level height change value of each liquid level acquisition point under the current working conditions is continuously acquired, and the liquid level disturbance value of each liquid level acquisition point is calculated based on the liquid level height change value. Under standard gas release reference conditions, the liquid level height change values ​​at each liquid level sampling location are continuously acquired, and the standard liquid level disturbance value at each liquid level sampling location is calculated based on the liquid level height change values. Extend the local gas release channel vertically to the liquid surface area to determine the liquid surface sampling location where the local gas release channel reaches the liquid surface; Starting from the liquid surface sampling position where the local gas release channel reaches the liquid surface, read the liquid surface disturbance value and the standard liquid surface disturbance value of adjacent liquid surface sampling positions one by one along both sides of the horizontal direction, and calculate the liquid surface disturbance offset value of each liquid surface sampling position. Compare the liquid surface disturbance offset values ​​of adjacent liquid surface sampling positions on both sides in the direction away from the starting point; retain the liquid surface sampling positions whose liquid surface disturbance offset values ​​are less than the previous liquid surface sampling position's liquid surface disturbance offset value, and stop further comparison in that direction when the liquid surface disturbance offset value is greater than or equal to the previous liquid surface sampling position's liquid surface disturbance offset value. The area enclosed by the liquid surface sampling positions on both sides of the horizontal axis and the liquid surface sampling position at the starting point is defined as the liquid surface response range; Within the liquid surface response range, the liquid surface sampling location with the largest liquid surface disturbance offset value is selected and determined as the liquid surface response center.

7. The method for predicting operational risks of a wastewater treatment plant according to claim 6, characterized in that, The methods for forming surface response distortion characteristics in step three include: Extend the local gas release channel vertically to the liquid surface area and determine the liquid surface sampling location corresponding to the local gas release channel; Extract the lateral position difference between the liquid surface response center and the liquid surface acquisition position corresponding to the local gas release channel. When the lateral position difference does not exceed the preset position deviation value, this state is taken as the positional correspondence between the liquid surface response center and the liquid surface position corresponding to the local gas release channel. The sum of liquid surface disturbance offset values ​​at each liquid surface acquisition position within the liquid surface response range is calculated. When the ratio of the sum of liquid surface disturbance offset values ​​at each liquid surface acquisition position within the liquid surface response range to the sum of liquid surface disturbance offset values ​​at each liquid surface acquisition position within the liquid surface range is greater than a preset concentration ratio, this state is regarded as a local concentration state of liquid surface disturbance. Extract the difference in liquid surface disturbance offset between the two boundary positions of the liquid surface response range and the adjacent liquid surface acquisition position outside it, and extract the difference in liquid surface disturbance offset between the adjacent liquid surface acquisition positions inside the liquid surface response range. When the difference in liquid surface disturbance offset between the two boundary positions of the liquid surface response range and the adjacent liquid surface acquisition position outside it is greater than the difference in liquid surface disturbance offset between the adjacent liquid surface acquisition positions inside the liquid surface response range, this state is regarded as the abrupt transition state of the liquid surface boundary. When the positional correspondence between the liquid surface response center and the corresponding liquid surface position of the local gas release channel, the local concentration state of liquid surface disturbance, and the abrupt transition state of the liquid surface boundary are all established simultaneously, surface response distortion characteristics are formed.

8. The method for predicting operational risks of a wastewater treatment plant according to claim 1, characterized in that, The methods for generating the reverse deviation characteristic of oxygen transfer results in step four include: Multiple dissolved oxygen collection points are set within the local oxygen transfer response range corresponding to the liquid surface response center; the dissolved oxygen values ​​at each collection point are continuously acquired under the current operating conditions. Within the preset recovery time interval, the rate of change and the magnitude of change of dissolved oxygen values ​​at each dissolved oxygen collection location were calculated. Calculate the average rate of change of dissolved oxygen values ​​at each dissolved oxygen sampling location, and determine the average rate as the local dissolved oxygen recovery rate; calculate the average magnitude of change of dissolved oxygen values ​​at each dissolved oxygen sampling location, and determine the average magnitude as the local dissolved oxygen recovery magnitude. Extract the standard local dissolved oxygen recovery rate and standard local dissolved oxygen recovery amplitude within the corresponding local oxygen transfer response range under the standard gas release reference state; When the gas release morphology shift characteristic and surface response distortion characteristic have been formed, and the local dissolved oxygen recovery rate is less than the standard local dissolved oxygen recovery rate, and the local dissolved oxygen recovery amplitude is less than the standard local dissolved oxygen recovery amplitude, the oxygen transfer result reverse deviation characteristic is formed.

Images