A Simulation and Optimization Control Method for Social Water Cycle in Industrial Enterprises

Abstract

This invention discloses a method for simulating and optimizing the social water cycle in industrial enterprises. Based on the various technological processes of the industrial enterprise and the stages of the social water cycle ("water intake-water use-water consumption-water discharge"), a simplified network diagram of the industrial enterprise's water system is drawn. Based on the physical mechanisms, theoretical relationships, logical relationships, and correlations of each water-using process stage, theoretical equations, empirical formulas, and correlation formulas for each water-using process stage of the industrial enterprise's water system are proposed. Based on system dynamics theory, a simulation model of the social water cycle in the industrial enterprise is constructed. This method ultimately achieves integrated optimization and efficient utilization of the entire plant's water system, achieving the comprehensive goals of water conservation and emission reduction. This breakthrough in the research of the social water cycle in industrial enterprises improves the fundamental theory of the social water cycle and promotes the development of water resources science.

Description

Technical Field

[0001] This invention relates to the field of hydrology, specifically to a method for simulating and optimizing the social water cycle in industrial enterprises. Background Technology

[0002] The social water cycle is a side cycle coupled with the natural water cycle process driven by the water demand of the economy and society. It is not only closely related to the development and evolution of the economic and social system, but also follows the basic laws of natural and social processes.

[0003] Industry is a typical unit with the most complete and intense water cycle in society. It generally has significant characteristics such as continuous water intake, continuous water supply, cascaded water use, water recycling, water reuse, reclaimed water reuse, and wastewater discharge. Moreover, it is greatly affected by market economic regulation and policy system, and is a microcosm of the social water cycle system.

[0004] Industry remains a major water consumer in my country. Despite decades of industrial development and water-saving renovations, overall water efficiency still lags significantly behind international advanced levels. Previous research on industrial water cycles has primarily focused on indicative calculations and evaluations, leaning towards macro-level studies without delving into specific industrial sectors and enterprises. The research lacks analysis of the principles of industrial water cycles and micro-level research. This study aims to comprehensively consider the entire industrial enterprise's water system as an organic whole. It will comprehensively examine various processes, equipment, and technologies involved in the water cycle, as well as all possible scenarios for wastewater reuse, regeneration, and recycling. Using the principles and technologies of process system integration, it will systematically analyze and optimize the dynamic mechanism of the water cycle pathway. Water will be used progressively according to quality requirements to improve the water system's reuse rate, thereby simultaneously reducing both fresh water consumption and wastewater discharge. Ultimately, it aims to achieve integrated optimization and efficient utilization of the entire plant's water system, achieving the comprehensive goals of water conservation and emission reduction. This breakthrough in industrial enterprise-based social water cycle research will improve the fundamental theory of social water cycles and promote the development of water resources science. Summary of the Invention

[0005] To address the aforementioned shortcomings in existing technologies, this invention provides a method for simulating and optimizing the social water cycle in industrial enterprises.

[0006] To achieve the above-mentioned objectives, the technical solution adopted by this invention is as follows:

[0007] A method for simulating and optimizing the social water cycle in industrial enterprises includes the following steps:

[0008] S1. Draw a generalized diagram of the industrial enterprise's water system network model based on the data of each process flow of the industrial enterprise and each link of the social water cycle "water intake-water use-water consumption-water discharge".

[0009] S2. Extract the system dynamics unit mathematical model of each water use system based on the physical mechanism of each water use process;

[0010] S3. Based on system dynamics theory, construct a simulation model of the social water cycle of industrial enterprises according to the relationship between system dynamics units, and optimize and regulate the water cycle according to the simulation model.

[0011] Furthermore, the simplified diagram of the industrial enterprise water system network model in S1 is specifically as follows:

[0012] Water usage will be restricted to a fixed flow rate of chemically demineralized water.

[0013] Ultrafiltration, reverse osmosis, and ion exchange water treatment processes all produce effluents of different concentrations. The impurity enrichment units of the chemical water treatment system and the integrated water treatment system are divided, and the units are connected by the water production rate α.

[0014] Furthermore, the water-using process in S2 is an analysis of the plant's impurity control strategy, including the concentrations of chloride ions, total suspended solids, and sulfate ions.

[0015] Total suspended solids control the water used in physical forward and reverse filtration processes, as well as the makeup water for domestic water systems;

[0016] Chloride ion concentration controls the concentration rate of the circulating cooling tower and the amount of water used for desulfurization.

[0017] Sulfate ions control the desulfurization water and wastewater processes.

[0018] Furthermore, step S3 specifically includes the following steps:

[0019] S31. Based on system dynamics theory, add additional constraints to reduce the dimensionality of the system dynamics unit mathematical model constructed in step S2.

[0020] S32. Perform nonlinear constraint transformation on the dimensionality-reduced system dynamics mathematical model;

[0021] S33. Determine the initial optimization point, solve the system dynamics model after nonlinear constraint transformation, and obtain the optimal water use regulation scheme.

[0022] Furthermore, the additional constraints in S31 include:

[0023] The water-using unit does not directly reuse its own water;

[0024] The ash and slag system and other water systems are terminal units, which do not drain water or supply water to the outside.

[0025] Chemical water treatment systems and integrated water treatment systems' enrichment units only accept water supplied by the treatment units; boiler water systems only accept demineralized water; and domestic water systems only accept demineralized water and fresh water.

[0026] Furthermore, step S32 specifically includes the following steps:

[0027] S321. Using the inlet and outlet flow rates of the water-using unit as known quantities, the impurity conservation constraint of the water-using unit is quantitatively represented.

[0028] S322. Add the inlet and outlet impurity concentration constraint as the penalty function to the objective function, and add a connection number optimization component;

[0029] S323. Set up a linear equation system for the unit outlet impurity concentration to correct for the negative outlet concentration caused by impurity load.

[0030] Furthermore, the specific method for quantifying the water-using unit impurity conservation constraint in S321 is as follows:

[0031]

[0032] Among them, F i,j For water-using unit A i Water supply unit A j The flow rate of impurities in the water supply; Let n be the concentration of impurities in the wastewater of the nth water-using unit, where n is the number of water sources. For the water consumption of the water-using unit, [M s ] is the loading matrix of impurity s in all cells, F w Fresh water supply volume The concentration of impurities s in fresh water.

[0033] Furthermore, the specific representation of adding the inlet and outlet impurity concentration constraint as a penalty function to the objective function in S322, and adding a connection number optimization component, is as follows:

[0034]

[0035]

[0036]

[0037] Among them, g + Let sgn be the auxiliary constraint function, μ be the sign function, μ be the penalty coefficient, and f′ be the penalized objective function. Each water source supplies water to water user unit A. i The inlet concentration and maximum inlet concentration of impurities s after water mixing. Water-using unit A iThe outlet concentration and maximum outlet concentration after enriching a certain amount of impurities s.

[0038] Furthermore, the specific method for correcting the negative outlet concentration caused by impurity load in S323 is as follows:

[0039]

[0040] Among them, Q k Let b be the impurity load at the k-th outlet. k C represents the maximum throughput of the k-th outlet. k This is the negative impurity concentration correction value for the k-th outlet.

[0041] The present invention has the following beneficial effects:

[0042] 1) The treatment unit is described by negative impurity load, and the impurity enrichment unit is described by positive impurity load. The treatment unit is set with a maximum impurity removal load. If the total impurity mass at the inlet is less than the maximum treatment capacity, the outflow is considered to be with the best water quality. Otherwise, the impurity enrichment process is used as an equivalent process.

[0043] (2) The special water-using subsystems in thermal power plants are equivalent to a fixed water production rate water treatment unit, a fixed water flow rate unit, and a terminal water-using unit, and are respectively treated as follows: the fixed water production rate water treatment unit is set with a fixed water production rate and is divided into a water treatment process and an impurity enrichment process, and the flow rate is allocated according to the set water production rate; the fixed flow rate unit meets the actual demand by setting a fixed inlet flow rate; the terminal water-using unit makes its water supply equal to the set water consumption.

[0044] (3) Introduce an adjacency matrix to describe the water supply relationship between units. Assume that each unit can supply water to other units and receive water from other units, and extend it to a water network model that can describe the "supply-use-consumption-discharge" relationship of the water-using unit network.

[0045] (4) Using the adjacency matrix as the optimization dependent variable and the water use mechanism of each unit as the constraint, the water use units are spliced ​​together to construct an industrial enterprise social water cycle optimization and control model, and the "water saving-energy saving" target can be used to solve the model. Attached Figure Description

[0046] Figure 1 This is a schematic diagram of the process for the simulation and optimization control method of social water cycle in industrial enterprises according to the present invention.

[0047] Figure 2 This is a schematic diagram of the water use process and water system division in a thermal power plant according to an embodiment of the present invention.

[0048] Figure 3 This is a schematic diagram of a typical power company water use relationship model according to an embodiment of the present invention.

[0049] Figure 4 This is a schematic diagram showing the equivalent breakdown of the water treatment process in an embodiment of the present invention.

[0050] Figure 5 This is the optimized water network diagram for this embodiment. Detailed Implementation

[0051] The specific embodiments of the present invention are described below to enable those skilled in the art to understand the present invention. However, it should be understood that the present invention is not limited to the scope of the specific embodiments. For those skilled in the art, various changes are obvious as long as they are within the spirit and scope of the present invention as defined and determined by the appended claims. All inventions utilizing the concept of the present invention are protected.

[0052] A method for simulating and optimizing the social water cycle in industrial enterprises, such as... Figure 1 As shown, it includes the following steps:

[0053] S1. Draw a generalized diagram of the industrial enterprise's water system network model based on the data of each process flow of the industrial enterprise and each link of the social water cycle "water intake-water use-water consumption-water discharge".

[0054] In this embodiment, based on the various technological processes of a typical industrial enterprise and the various stages of the social water cycle ("water intake-water use-water consumption-water discharge"), a simplified network diagram of the water system of a typical thermal power plant industrial enterprise is drawn, such as... Figure 2 As shown.

[0055] The boiler water system (water-steam circulation system) in the plant area is regulated by the plant's power generation load and can only use demineralized water for replenishment, making optimized design impossible. Therefore, boiler water use is limited to a fixed flow rate of chemically demineralized water. The cooling water system is not a fully functional impurity enrichment unit; circulating water consumption depends on unit output and local weather conditions, and drainage is controlled by the water tower concentration ratio. Flue gas desulfurization water consumption and impurity load are related to coal consumption and coal type. Ignoring the influence of external factors on desulfurization evaporation and open-loop cooling tower evaporation, the measured impurity load is used as the calculated value. Ultrafiltration, reverse osmosis, and ion exchange water treatment processes all produce effluents of varying concentrations. Therefore, separate impurity enrichment units are established for the chemical water treatment system and the integrated water treatment system, connected by a water production rate α to satisfy the mass transfer process model. The equivalent decomposition of the water treatment process is shown in [link to model]. Figure 4 .

[0056] S2. Extract the system dynamics unit mathematical model of each water use system based on the physical mechanism of each water use process;

[0057] Based on the physical mechanisms of each water-using process, the mathematical models of the system dynamics units of each water-using system are compiled, as shown in Table 1:

[0058] Table 1 Mathematical Model of System Dynamics Model Unit

[0059]

[0060]

[0061] The case study company's water usage data comes from the company's water balance test results and actual measurement data in the plant area. Analysis of the plant's main impurity control strategies, particularly chloride ions (Cl...) - ), total suspended solids (SS) and sulfate ions The concentration of SS (suspended solids) best reflects the water use control process in the water network of a thermal power plant. SS mainly controls the water used in physical forward and reverse filtration processes; water with high SS levels cannot be used as makeup water for chemical demineralization systems or domestic water systems. Cl (Cl-) - Concentration mainly affects the concentration ratio of the circulating cooling tower and the amount of desulfurization water used, making it an important control indicator; Concentration affects the desulfurization water and drainage process; other ions need to meet the boiler water requirements, but do not dominate in the chemical desalination system, as shown in Table 2.

[0062] Table 2 Optimization data for each water use unit

[0063]

[0064] S3. Based on system dynamics theory, construct a simulation model of the social water cycle of industrial enterprises according to the relationship between system dynamics units, and optimize and regulate the water cycle according to the simulation model.

[0065] In this embodiment, based on system dynamics theory, and according to the theoretical, logical, and interrelationships between various system dynamics units, a simulation model of the social water cycle in an industrial enterprise is constructed, such as... Figure 3 As shown.

[0066] We selected typical thermal power plant industrial enterprises and conducted practical research on the integrated optimization and regulation of social water cycle with the goal of water conservation and emission reduction.

[0067] The plant's water usage processes are divided into nine systems: domestic and fire-fighting water system, phase I circulating cooling water system, phase II circulating cooling water system, boiler water system, desulfurization water system, ash and slag system, integrated water treatment system, chemical water treatment system, and other water systems. Equivalent treatment of each unit and selection of key impurity indicators are then performed to construct an optimization model. The specific steps include the following:

[0068] S31. Based on system dynamics theory, add additional constraints to reduce the dimensionality of the system dynamics unit mathematical model constructed in step S2.

[0069] Additional constraints: Adding extra constraints can make the model results more reasonable and can also reduce the search space, thus simplifying the model. Dimensionality reduction is achieved through the following assumptions:

[0070] A. The unit does not directly reuse itself;

[0071] B. The ash and slag system and other water systems are the terminal units, which do not discharge water or supply water to the outside.

[0072] C. Chemical water treatment systems and integrated water treatment systems' enrichment units only accept water supplied by the treatment units; boiler water systems only accept demineralized water; and domestic water systems only accept demineralized water and fresh water.

[0073] S32. Perform nonlinear constraint transformation on the dimensionality-reduced system dynamics mathematical model;

[0074] This embodiment includes the following steps:

[0075] S321. Using the inlet and outlet flow rates of the water-using unit as known quantities, the impurity conservation constraint of the water-using unit is quantitatively represented.

[0076] Since the inlet and outlet concentrations of each unit need to be within the limits set for each impurity inlet and outlet, the calculation of the impurity concentration at the unit inlet depends on the outlet concentration of each water supply unit, while unit A... i Export concentration The calculations require the conservation of impurities within each unit. Adding an outlet concentration optimization variable to the control variables for each unit would not only create additional optimization variables but also introduce nonlinear equality constraints (impurity conservation constraints), making the model difficult to solve. Therefore, we treat the inlet and outlet flow rates of each unit as known quantities and express the unit impurity conservation constraints in matrix form.

[0077]

[0078] Among them, F i,j For water-using unit A i Water supply unit A j The flow rate of impurities in the water supply; Let n be the concentration of impurities in the wastewater of the nth water-using unit, where n is the number of water sources. For the water consumption of the water-using unit, [M s ] is the loading matrix of impurity s in all cells, F w Fresh water supply volume The concentration of impurities s in fresh water.

[0079] S322. Add the inlet and outlet impurity concentration constraint as the penalty function to the objective function, and add a connection number optimization component;

[0080] Given the unit flow rates, if the entire water network satisfies the flow conservation relationship, then a solution must exist. Inlet and outlet impurity concentration constraints are added to the objective function as penalty functions, and a connectivity optimization component is added, expressed as:

[0081]

[0082]

[0083]

[0084] Among them, g + Let sgn be the auxiliary constraint function, μ be the sign function, μ be the penalty coefficient, and f′ be the penalized objective function. Each water source supplies water to water user unit A. i The inlet concentration and maximum inlet concentration of impurities s after water mixing. Water-using unit A i The outlet concentration and maximum outlet concentration after enriching a certain amount of impurities s.

[0085] S323. Set up a linear equation system for the unit outlet impurity concentration to correct for the negative outlet concentration caused by impurity load.

[0086] For processing units, it is necessary to correct for negative outlet concentrations caused by negative impurity loads. When the amount of impurities removed is less than the maximum removal amount, the outlet concentration of the impurity is considered to be 0. Let Q·C=b be a system of linear equations for solving the outlet impurity concentration of the unit. If the impurity load of the unit is negative, it is judged to be a processing unit, and if the outlet concentration is negative, it is judged to be that the maximum processing capacity has not been reached. Then, a correction calculation is performed:

[0087] C = Q -1 ·b

[0088] If there are k values ​​that need to be corrected, and C is known but b is unknown, then the linear equations corresponding to the concentration units that need to be corrected can be extracted from the quantitative expression of the impurity conservation constraint of the water unit, forming a system of k linear equations. (Since Q is full rank, Q) -1 If the rank is full, then there must be a solution, and the corresponding b can be calculated. k The modified C is calculated using the quantification expression of the impurity conservation constraint for the rewater unit, thus concluding the modification process. Integrating the inlet and outlet concentration calculations into the objective function removes nonlinear equality constraints while reducing the dimensionality of optimization variables. However, the objective function contains concentration calculations and penalties added to concentration and connectivity, making it nonlinear and not guaranteed to be smooth, which significantly complicates the solution.

[0089] The negative concentration correction method is as follows:

[0090] The correction will cause the inlet and outlet concentrations of all units except the correction unit to increase or remain unchanged. Furthermore, there will be an interaction between multiple negative concentrations that need to be corrected. This process is extremely complex when considering the water network flow state. To solve this problem, for the units that need to be corrected, only two states are considered: reaching the maximum processing capacity and the minimum outlet concentration. First, the correction is performed according to the minimum outlet concentration. If a unit exceeds the maximum processing capacity, the unit is corrected according to the maximum processing capacity. If 0 < i < j < k are the concentrations of 3 negative impurities that need to be corrected, and their concentration correction values are C′ i , C′ j , C′ k , and Q·C′ = b′, then b′ can be obtained by solving the following system of equations i , b′ j , b′ k and the corrected C′:

[0091]

[0092] If there is an element in b′ that exceeds the maximum value and this value is in the p-th row, then at the corresponding position p there is and C′ p changes from the known minimum value to an unknown value for solution. Corresponding to the above formula, only C′ p is used as the position that does not need to be corrected for re-solving. This method is also applicable to correcting multiple positions simultaneously. Finally, the obtained C′ and b′ are the final results.

[0093] S33. Determine the optimization initial point, solve the system dynamics model after nonlinear constraint transformation, and obtain the optimal water use regulation plan.

[0094] Determination of the initial point: The solution of the model highly depends on the initial point x0. For any x0 (the row vector obtained by splicing the rows of the unit adjacency matrix), the equality constraint Aeq·x0 = beq after Gaussian elimination can be used. Using the values of [x r+1 , …, x n to correct the values of [x1, …, x r to make it satisfy the equality constraint. Assign equal values to x0 in sequence, and then solve this nonlinear optimization problem (NLP) with these initial points respectively, and obtain the final non-dominated solution set. In this paper, the "Trust-RegionInterior-Point Method" nonlinear optimization algorithm in MATLAB is used to solve the model, and the solution results are as shown in Figure 5 and Table 3.

[0095] Table 3 Comparison table of optimization effects

[0096]

[0097]

[0098] After optimization, the fresh water consumption decreased by 220.96 t / h, a 6.94% reduction compared to the original water consumption, and direct sewage discharge decreased by 65%. The optimized integrated water treatment system reduced its treatment capacity by 11.34 t / h, while the chemical water treatment system increased its capacity by 45.16 t / h, and desulfurization water consumption decreased by 18 t / h. The increased water consumption in the domestic and fire-fighting water systems can be considered as treated domestic sewage mixed with this portion for reuse. The concentrated water from both the integrated and chemical water treatment systems is directly returned to the inlet to mix with other water and then re-enter the treatment system for further treatment, achieving direct reuse. Under the three impurity indicators, the optimized inlet and outlet concentrations meet the limits for each inlet and outlet, and each has a certain impurity concentration close to the maximum outlet concentration, consistent with the maximum mass transfer theory. The integrated water treatment, chemical water treatment, and FGD calculations all use the maximum mass transfer process of the system for calculating the inlet and outlet concentrations. The integrated water treatment can handle total suspended solids (SS) and sulfate. The treatment is effective, with SS removal approaching its maximum removal capacity; the chemical water treatment system effectively removes Cl. - The removal is close to the maximum removal capacity; the desulfurization water system only removes water near its maximum capacity. It has the removal capacity and has reached its maximum removal capacity. The inlet and outlet impurity concentrations of each unit after optimization are shown in Table 4, and the impurity removal status of integrated water, chemical water, and desulfurization are shown in Table 5.

[0099] Table 4. Optimized unit inlet and outlet water quality (mg / L)

[0100]

[0101] Table 5 Utilization Rate of Impurity Removal Capacity of Water Treatment System (M / M) max )

[0102]

[0103] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1A device that provides the functions specified in one or more boxes.

[0104] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0105] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0106] Specific embodiments have been used to illustrate the principles and implementation methods of this invention. The descriptions of the embodiments above are only for the purpose of helping to understand the method and core ideas of this invention. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of this invention. Therefore, the content of this specification should not be construed as a limitation of this invention.

[0107] Those skilled in the art will recognize that the embodiments described herein are intended to help the reader understand the principles of the invention, and should be understood that the scope of protection of the invention is not limited to such specific statements and embodiments. Those skilled in the art can make various other specific modifications and combinations based on the technical teachings disclosed in this invention without departing from the spirit of the invention, and these modifications and combinations are still within the scope of protection of this invention.

Claims

1. An industrial enterprise social water cycle simulation, emulation and optimization control method, characterized in that, Comprise the following steps: S1, according to the industrial enterprise each technological process and social water cycle "water intake-water consumption-water consumption-drainage" each link data draw industrial enterprise water system network model generalization figure; S2, according to the physical mechanism of each water use process link, the system dynamics unit mathematical model of each water use system is extracted; S3, based on the theory of system dynamics, according to the relationship between system dynamics units, an industrial enterprise social water cycle simulation model is constructed, and water cycle optimization control is carried out according to the simulation model, specifically including the following steps: S31, based on the theory of system dynamics, additional constraint conditions are added to the system dynamics unit mathematical model constructed in step S2 to reduce dimension; S32, the nonlinear constraint transformation is carried out on the reduced system dynamics mathematical model, specifically including the following steps: S321, the inlet and outlet flow of the water unit is used as a known quantity, and the impurity conservation constraint of the water unit is quantitatively represented, specifically as follows: wherein, is the water using unit is the water using unit is the impurity flow rate of the water supply; is the effluent impurity concentration of the nth water using unit, n is the number of water sources, is the water consumption of the water using unit, is the load matrix of impurity s in all units, is the fresh water supply water volume, is the concentration of impurity s contained in the fresh water; S322, the inlet and outlet impurity concentration constraint is added to the objective function as a penalty function, and a connection number optimization component is added, specifically represented as follows: in, For auxiliary constraint functions, For symbolic functions, The penalty coefficient, Let the objective function be the one after the penalty is applied. , Each water source is directed to the water-using unit. The inlet concentration and maximum inlet concentration of impurities s after water mixing. , Each is a water-using unit The outlet concentration and maximum outlet concentration after enriching a certain amount of impurity s; S323, a linear equation set of unit outlet impurity concentration is set, and the negative outlet concentration caused by impurity load is corrected, specifically as follows: wherein, is the impurity load for the k th outlet, is the maximum throughput for the k th outlet, is the negative impurity concentration correction value for the k th outlet; S33, the optimization initial point is determined, the system dynamics model after nonlinear constraint transformation is solved, and the optimal water control scheme is obtained.

2. The industrial enterprise sociohydrologic cycle simulation and optimization method of claim 1, wherein, The water system network model generalization figure in S1 is specifically: The water use is limited to only use fixed flow chemical desalted water; The ultrafiltration, reverse osmosis, ion exchange water treatment process all exist different concentration of effluent, divide out chemical water treatment system impurity enrichment unit, comprehensive water treatment system impurity enrichment unit, with water production rate Carry out the connection between units.

3. The industrial enterprise sociohydrologic cycle simulation and optimization method of claim 1, wherein, The water use process link in S2 is to analyze the impurity control strategy of the plant area, including: the concentration of chloride ions, total suspended solids and sulfate ions, wherein, The total suspended solids control the water use of the physical positive filter and reverse filter process and the water supplement of the domestic water system; The chloride ion concentration controls the concentration ratio of the circulating cooling tower and the desulfurization water consumption; The sulfate ion controls the desulfurization water and the drainage process.

4. The industrial enterprise sociohydrologic cycle simulation and optimization method of claim 1, wherein, The additional constraint conditions in S31 include: The water unit is not directly reused itself; The ash system and other water use systems are terminal units, which do not drain and do not supply water to the outside; The chemical water treatment system and the enrichment unit of the comprehensive water treatment system only accept the water supply of the treatment unit;The boiler water system only accepts desalted water;The domestic water system only accepts desalted water and fresh water.

Images