A method for dynamically defining responsibility for cross-border pollution in a large-scale complex river network region
By identifying transboundary pollution responsibility areas using a one-dimensional steady-state water quality model and river hydrodynamic characteristics, and calculating dynamic allowable discharges and contribution weights, the ambiguity in defining transboundary pollution responsibility in large-scale complex river network areas has been resolved, enabling accurate determination of pollution responsibility and dynamic control.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- NINGBO UNIV
- Filing Date
- 2026-04-09
- Publication Date
- 2026-07-24
AI Technical Summary
In the management of cross-border pollution in large-scale complex river network areas, traditional methods cannot accurately identify the specific contribution of pollution emissions from each administrative or functional unit to cross-border pollution, resulting in a vague management scope, difficulty in defining responsibilities, and neglect of the impact of river network characteristics and location distribution on pollution transport.
By employing a one-dimensional steady-state water quality model combined with river hydrodynamic characteristics, effective responsibility delineation areas for cross-border pollution are identified. Through iterative calculations, the dynamic allowable discharge and contribution weight of pollution control spatial units are determined, thereby achieving refined management and control level classification.
It has enabled precise identification of responsibility and dynamic control over cross-border pollution, overcome the rigidity of traditional methods, ensured the objectivity and precision of management, and provided operable technical support for the governance of water environment in large-scale river network areas.
Smart Images
Figure CN122452912A_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of water environment management and pollution control technology, specifically a method for dynamic liability delineation of cross-border pollution in large-scale complex river network areas. Background Technology
[0002] The management of transboundary pollution and the delineation of pollution responsibility in large-scale plain river network areas have always been challenging issues for river basin authorities and local government environmental protection departments. These issues directly impact water resource protection, water security, and even harmonious development within these areas. Due to the extensive coverage of plain river network areas and their involvement of multiple administrative regions, characterized by dense and widespread pollution discharges, intricate river systems, and complex and unpredictable pollution transport processes, the relationship between dispersed pollution discharges and transboundary water body receptors is often ambiguous, making it difficult to accurately identify the specific contribution of pollution discharges from each administrative or functional unit to transboundary pollution. This results in significant problems in river basin and cross-administrative region management, including unclear management scope and mutual shirking of responsibility among different regions. Against the backdrop of my country's ecological civilization construction, the impact of this problem on pollution control and water ecological protection in the southeastern coastal river network areas has become increasingly prominent.
[0003] First, traditional methods for defining liability for transboundary river pollution are mostly based on clear upstream and downstream responsibilities. They involve upstream and downstream consultation and decision-making through cross-sectional pollution flux monitoring. However, in large-scale complex river network areas with multiple polluting entities and unclear upstream and downstream relationships in multi-directional connected river systems, the previous work on defining and controlling transboundary pollution liability cannot be carried out smoothly. Secondly, existing large-scale pollution control and pollution liability determination work in river basins and administrative regions mostly rely on statistics and calculations of the amount of pollution discharged in each administrative region as the principle for determining pollution liability. This ignores the different pollution carrying capacity or pollution tolerance of each responsible area due to the different characteristics of their river network systems, and thus cannot achieve accurate and objective determination of pollution liability. Finally, since pollutants often diffuse and attenuate during their transport to the transboundary river after entering the transboundary water body, the magnitude of the pollution transport flux will vary depending on the spatial distribution of the pollution source. Therefore, the impact on the water quality of the transboundary section will also vary depending on the location of the pollution source. Traditional pollution liability determination and pollution control cannot reflect the spatial distance relationship caused by the location distribution of each administrative unit.
[0004] Therefore, there is an urgent need to construct a set of methods for defining and intelligently managing cross-border pollution responsibility applicable to large-scale river network areas. This method should take the identification of effective impact areas of cross-border pollution as a prerequisite, the dynamic over-pollution capacity of each spatial unit as the judgment criterion, the amount of excessive emissions as the definition basis, and consider the process attenuation factor for correction and verification. Summary of the Invention
[0005] To address the shortcomings of existing technologies, this invention provides a method for dynamic liability delineation of cross-border pollution in large-scale complex river network areas. This method solves important technical problems in large-scale plain river network areas, such as unclear boundaries for cross-border pollution control and difficulties in delineating liability in each area due to the complex topology of the water system, the dispersed nature of pollution sources, and dynamic emissions.
[0006] To achieve the above objectives, the present invention provides the following technical solution: a method for dynamic liability delineation of transboundary pollution in large-scale complex river network areas, comprising the following steps: Based on the hydrodynamic characteristics of important transboundary rivers in the plain river network area, the distribution and emission characteristics of pollution sources along the river, and through calculation using a one-dimensional steady-state water quality model along the river, a trial algorithm was used to determine the district and county-level administrative regions whose pollution impact weights on national and provincial control transboundary sections are within a set range as effective responsibility delineation areas for transboundary pollution. The identified effective responsibility delineation area for cross-border pollution is used as the system boundary. The county and district-level administrative units within the effective responsibility delineation area for cross-border pollution are subdivided into several pollution control spatial units of different scales within a set range. For the multiple pollution control spatial units obtained by division, based on the river distribution, river hydrodynamic dynamics and water quality control requirements within each pollution control spatial unit, the dynamic allowable discharge of each pollution control spatial unit that can reflect the dynamic pollution carrying capacity of the river is calculated. Based on the dynamic allowable emissions of each control space unit, the dynamic pollution exceedance emissions of each pollution control space unit that exceed the benchmark are determined as the benchmark for liability determination. The influence coefficient of pollution discharge on the transport process of water quality at cross-boundary sections within each pollution control spatial unit is calculated using a one-dimensional steady-state water quality model. This coefficient is then used as a correction coefficient to correct the impact of dynamic pollution exceeding standards in each pollution control spatial unit and to determine the effective dynamic contribution weight of each pollution control spatial unit to cross-boundary pollution. This allows for a refined classification of the control levels of each pollution control spatial unit.
[0007] Furthermore, the iterative one-dimensional steady-state water quality model along the river is determined based on the river length and calculation accuracy, and takes the following form: ; ; ...; ; in, This refers to the concentration of the upstream water in the first section; (This refers to the upstream water flow rate in the first section). This refers to the concentration of the downstream effluent from the first stage. This refers to the first stage of water flow rate; The concentration of pollutants in the first tributary or generalized discharge outlet; This represents the flow rate of the first tributary or generalized outlet; the process is iterated and repeated in this way. The pollutant concentration at the tributary or generalized discharge outlet N in the N+1th segment; This refers to the flow rate into a tributary or a generalized discharge outlet. The pollutant concentration at the tributary or generalized discharge outlet N in the N+1th segment; The flow rate into a tributary or a generalized discharge outlet; The concentration of pollutants effluent in the (N+1)th segment; Outflow rate; The degradation coefficient of the target pollutant; The average flow velocity of the river channel; The distance is the flow distance of the N+1 segment of the river channel.
[0008] Furthermore, the calculation process for the pollution impact weight is as follows: Based on a one-dimensional steady-state water quality model along the river, the water quality concentration of major transboundary rivers at the inflow boundaries of the trial county and district administrative regions is used as the inflow pollution concentration. Based on the hydrodynamic characteristics of the transboundary river and the pollution input results of all tributaries and generalized discharge outlets along the river, the water quality concentration at the national and provincial control transboundary sections under the influence of pollution sources along the river is calculated iteratively. ; Simultaneously, the water quality concentration of the blank control group was calculated in the iterative calculation, ignoring the impact of pollution along the route. ; Obtain pollution impact weight expression: .
[0009] Furthermore, the process for calculating the dynamic allowable discharge of each pollution control spatial unit, which reflects the dynamic pollution carrying capacity of the river channel, is as follows: The characteristic parameters of each pollution control spatial unit are obtained, including flow rate, flow velocity, actual pollutant concentration and target pollutant concentration; The characteristic parameters are input into the pollution holding capacity calculation formula to calculate the pollution holding capacity, where the pollution holding capacity calculation formula is in the form of: ; in, The dynamic contamination-carrying capacity of space unit i on day k. The number of sewage-receiving river sections within a spatial unit that have a water environment function positioning; Indicates the first unit within the cell The first section of the river The average flow rate on the k-th day of the day; The first one determined according to the water environment function positioning The first effective control unit Target concentration of pollutants at the inflow section of each river segment; For the first The first effective control unit Target concentration of pollutants at the outflow section of each river segment; Indicates the first The first effective control unit The first-level comprehensive degradation coefficient of each river section; For the first The first effective control unit The length of each river section; For the first The first effective control unit The average daily flow velocity on day k of each river segment.
[0010] Furthermore, the process for obtaining dynamic pollution exceedance data for each pollution control spatial unit is as follows: Determine the number of point sources, land use area, and pollution loss intensity within each effective control unit, and calculate the total pollution load: ; in, For the first Total pollution emission load of each spatial unit on day k; For the first Within the effective spatial unit, the first The average daily emission load of individual point sources is mostly average. and The first The area of land use type h within each effective spatial unit and its pollution loss intensity on day k; Calculate the pollution overload of each spatial unit. : .
[0011] Furthermore, the process of revising the impact of dynamic pollution exceeding emission standards in each pollution control spatial unit and determining the effective dynamic contribution weight of each pollution control spatial unit to transboundary pollution is as follows: Calculate the sum of pollution overloads for all effective control units; Calculate the correction factor for pollution emission process in spatial units. The calculation formula is as follows: ; In the formula, k j To control the degradation coefficient of transboundary river j, Q j To control the flow of the transboundary river j; j v represents the distance that the transboundary river j flows from its entry point to the transboundary section within the control unit; jThe average flow velocity of the transboundary river j from its entry point to the transboundary section; The pollution excess load of each effective control unit is normalized based on the sum of the pollution excess loads, and the dynamic pollution responsibility contribution rate is obtained by combining the pollution transport process correction coefficients calculated for each unit. : .
[0012] Furthermore, the process of classifying the control levels of each pollution control spatial unit into more refined categories is as follows: If the dynamic control index is not less than the first threshold, it is classified as a key control area; If the dynamic control index is less than the first threshold and not less than the second threshold, it is classified as a general control zone. If the dynamic control index is less than the second threshold, it is classified as a light control zone; If the dynamic control index is less than 0, it is classified as an area with no impact.
[0013] The present invention has the following beneficial effects: This invention first uses a one-dimensional steady-state water quality model to screen out control units that have an effective impact on the water quality of the assessment section, ensuring the relevance of subsequent analyses. Then, it utilizes a coupled model to reconstruct the dynamic process of river network hydrodynamics and pollutant migration and transformation, avoiding distortion in pollution carrying capacity assessment caused by static calculations. By accurately calculating the over-discharge load of point and non-point sources, pollution liability determination focuses on the actual exceeding limits. Combining multi-dimensional factors to calculate dynamic correction coefficients and control indices, the control level is dynamically adjusted according to hydrological conditions and pollution load. This overcomes the rigidity of traditional one-size-fits-all control and, through the coordinated use of various technical means, ensures the objectivity of pollution liability determination and the accuracy of control measures, providing operable and traceable technical support for the refined management of the water environment in plain river network areas.
[0014] Of course, any product implementing this invention does not necessarily need to achieve all of the advantages described above at the same time. Attached Figure Description
[0015] Figure 1 This is the overall flowchart of the present invention.
[0016] Figure 2 A schematic diagram illustrating the process of identifying the scope of cross-boundary responsibility in a complex river network area.
[0017] Figure 3 A schematic diagram showing the results of the responsibility delineation in a complex river network area.
[0018] Figure 4 Map showing the distribution of spatial units within the defined area and the amount of emissions exceeding standards in each unit (day n).
[0019] Figure 5 A diagram showing the corrected cross-border pollution contribution rate of each unit within the defined area (day n). Detailed Implementation
[0020] Please refer to the flowchart. Figure 1 This invention provides a technical solution: a method for dynamic liability delineation of cross-border pollution in large-scale complex river network areas, comprising the following steps: The first step is to identify the scope of responsibility and divide spatial units: Based on the hydrodynamic characteristics of important transboundary rivers in the plain river network area, the distribution and emission characteristics of pollution sources along the river, and through calculations using an iterative one-dimensional steady-state water quality model along the river, a trial-and-error algorithm is used to determine that the district and county-level administrative regions with an influence weight of 60%-80% on national and provincial control transboundary sections are the effective responsibility delineation areas for transboundary pollution. This area will serve as the subsequent scope for transboundary pollution responsibility delineation and control, referring to case studies. Figure 2 , Figure 3 Using the identified effective responsibility demarcation zone for cross-border pollution as the system boundary, the county and district-level administrative units within the scope are further subdivided into several scales (100-500 km). 2 The pollution control spatial unit serves as the smallest control spatial unit for subsequent responsibility definition and management, as referenced in the case study. Figure 4 .
[0021] The second step is the calculation of dynamic pollution exceedance emissions for each spatial unit: For the multiple pollution control spatial units identified, based on the internal river distribution, hydrodynamic characteristics, and water quality control requirements, the dynamic allowable emissions for each control spatial unit are calculated to reflect the river's dynamic pollution carrying capacity. Based on these calculations, the dynamic pollution exceedance emissions exceeding the baseline for each pollution control spatial unit are determined as the benchmark for liability determination, referencing a case study. Figure 4 .
[0022] The third step is to define the dynamic responsibility for cross-boundary pollution in each spatial unit: This involves calculating the impact coefficient of pollution emissions within each pollution control spatial unit on the transport process of water quality at the cross-boundary section as a correction coefficient. This corrects for the impact of dynamic pollution exceeding standards in each unit and determines the effective dynamic contribution weight of each spatial unit to cross-boundary pollution, enabling refined management and control level classification for each spatial unit. (Refer to the case study.) Figure 5 .
[0023] Furthermore, the iterative one-dimensional steady-state water quality model along the river is determined based on the river length and calculation accuracy, and takes the following form: ; ; ...; ; in, The concentration of the upstream water in the first section (mg / L); The upstream inflow rate of the first section (m³) 3 / s); The concentration of the downstream effluent from the first stage (mg / L); The first stage of outflow rate (m³) 3 / s); The pollutant concentration (mg / L) at the first tributary or generalized discharge outlet. The flow rate (m³) of the first tributary or generalized outlet. 3 / s); and so on, iterating in sequence. The pollutant concentration (mg / L) at the tributary or generalized discharge outlet N in the N+1 segment. The flow rate (m³) of a tributary or generalized outlet. 3 / s); The pollutant concentration (mg / L) at the tributary or generalized discharge outlet N in the N+1 segment. The flow rate (m³) of a tributary or generalized outlet. 3 / s); The effluent pollutant concentration (mg / L) for the N+1th segment. Outflow rate; The degradation coefficient (d) of the target pollutant -1 ); The average flow velocity in the river channel is (m / s). The distance (m) for the N+1 segment of the river channel.
[0024] This method, through spatial node iteration, can accurately analyze the spatial distribution and cumulative contribution of multiple pollution sources along the river, avoiding the distortion of water quality response caused by the simplification and merging of discharge outlets in conventional methods. Its advantage lies in truly reflecting the gradual degradation and superposition process of pollutants in the river network, significantly improving the accuracy and spatial refinement of identifying the causal relationship between cross-boundary section water quality and upstream responsibility areas.
[0025] Furthermore, the calculation process for the pollution impact weight is as follows: Based on a one-dimensional steady-state water quality model along the river, the water quality concentration of major transboundary rivers at the inflow boundaries of the trial county and district administrative regions is used as the inflow pollution concentration. Based on the hydrodynamic characteristics of the transboundary river and the pollution input results of all tributaries and generalized discharge outlets along the river, the water quality concentration at the national and provincial control transboundary sections under the influence of pollution sources along the river is calculated iteratively. ; Simultaneously, the water quality concentration of the blank control group was calculated in the iterative calculation, ignoring the impact of pollution along the route. ; Obtain pollution impact weight expression: .
[0026] In identifying effective responsibility zones for cross-border pollution, calculating the pollution impact weight of each spatial unit plays a crucial quantitative supporting role. By employing a trial-and-error algorithm to determine the impact coefficient of the target area, the proportion of the target area's contribution to the exceedance of water quality standards at the cross-section is quantified. Based on the trial calculation results, it is determined whether the weight value reaches a preset threshold (e.g., 60%-80%), thus scientifically judging whether a certain area should be included in the effective responsibility zone. This process ensures the objectivity, quantifiability, and verifiability of responsibility zone boundary delineation, avoiding subjective assumptions.
[0027] Furthermore, the process for calculating the dynamic allowable discharge of each pollution control spatial unit, which reflects the dynamic pollution carrying capacity of the river channel, is as follows: The characteristic parameters of each pollution control spatial unit are obtained, including flow rate, flow velocity, actual pollutant concentration and target pollutant concentration; The characteristic parameters are input into the pollution holding capacity calculation formula to calculate the pollution holding capacity, where the pollution holding capacity calculation formula is in the form of: ; in, Let represent the dynamic pollution carrying capacity (g / d) of spatial unit i on day k. The number of sewage-receiving river sections within a spatial unit that have a water environment function positioning; Indicates the first unit within the cell The first section of the river The average flow rate on the k-th day of the day (m³) 3 / s); The first one determined according to the water environment function positioning The first effective control unit Target concentration of pollutants (mg / L) at the inflow section of each river segment; For the first The first effective control unit Target concentration of pollutants (mg / L) at the outflow section of each river segment. Indicates the first The first effective control unit The first-order comprehensive degradation coefficient (d) of each river section -1 ); For the first The first effective control unit The length of each river segment (m); For the first The first effective control unit The average daily flow velocity (m / s) on day k of each river segment.
[0028] The core advantage of calculating the permissible discharge capacity of a unit based on the dynamic pollution carrying capacity of a river lies in resolving the mismatch between static total emission control and dynamic water environment capacity. This method accurately defines the compliant discharge threshold for each unit under dynamic environments by verifying the maximum permissible discharge under different hydrodynamic conditions in real time. This provides a scientific benchmark for determining "excessive discharge" that aligns with objective environmental capacity, avoiding misjudgments of responsibility due to sudden drops in environmental capacity during low water levels. It shifts the determination of responsibility for cross-border pollution from static management based on fixed indicators to dynamic and precise control in response to natural fluctuations.
[0029] Furthermore, the process for obtaining dynamic pollution exceedance data for each pollution control spatial unit is as follows: Determine the number of point sources, land use area, and pollution loss intensity within each effective control unit, and calculate the total pollution load: ; in, For the first Total pollution load (g / d) of each spatial unit on day k. For the first Within the effective spatial unit, the first The average daily emission load (g / d) of each point source is mostly average; and The first The area (m²) of land use type h within each effective spatial unit. 2 ) and its pollution loss intensity on day k (g / d) m 2 This intensity is related to daily rainfall; Calculate the pollution overload of each spatial unit. : .
[0030] Acquiring dynamic data on excessive emissions from each pollution control unit enables precise identification of pollution stakeholders and quantification of their "excessive" contributions. Compared to static total emissions, this method can separate background emissions from compliant emissions, capturing in real-time the impact of instantaneous events such as illegal discharges and rainy season overflows on downstream sections. This provides a crucial basis for scientifically defining the dynamic pollution responsibility weights among units in complex river networks and avoiding the unfair phenomenon of "compliant entities bearing the blame," significantly improving the scientific rigor and timeliness of cross-boundary compensation mechanisms.
[0031] Furthermore, the process of revising the impact of dynamic pollution exceeding emission standards in each pollution control spatial unit and determining the effective dynamic contribution weight of each pollution control spatial unit to transboundary pollution is as follows: Calculate the sum of pollution overloads for all effective control units; Calculate the correction factor for pollution emission process in spatial units. The calculation formula is as follows: ; In the formula, k j To control the degradation coefficient (d) of transboundary river j -1 ), Q j To control the flow rate (m) of the transboundary river j 3 / s); l j The distance (m) that transboundary river j flows from its entry point to the transboundary section within the control unit; v j The average flow velocity (m / s) of the transboundary river j from its entry point to the transboundary section. The pollution excess load of each effective control unit is normalized based on the sum of the pollution excess loads, and the dynamic pollution responsibility contribution rate is obtained by combining the pollution transport process correction coefficients calculated for each unit. : .
[0032] This method introduces spatial weights reflecting transport distance, degradation coefficient, and hydraulic conditions to correct the instantaneous excess emissions of each unit to an effective equivalent impact on downstream cross-border sections. This allows pollution contributions from different spatial locations to be compared and summed on a uniform and equitable scale, thereby accurately quantifying the true dynamic responsibility weight of each unit for cross-border pollution.
[0033] Furthermore, the process of classifying the control levels of each pollution control spatial unit into more refined categories is as follows: If the dynamic control index is not less than the first threshold, it is classified as a key control area; If the dynamic control index is less than the first threshold and not less than the second threshold, it is classified as a general control zone. If the dynamic control index is less than the second threshold, it is classified as a light control zone; If the dynamic control index is less than 0, it is classified as an area with no impact.
[0034] The core advantage of implementing refined management and control levels for each pollution control spatial unit lies in achieving optimized allocation of management resources and differentiated policy implementation. Classifying units based on indicators such as pollution load contribution, remaining environmental capacity, and dynamic exceedance risk (e.g., key control, general control) can change the previous "one-size-fits-all" management model. Stricter emission standards and higher-frequency monitoring can be implemented for higher-level units, while risk prevention is emphasized for lower-level units.
[0035] To address the shortcomings of existing technologies, this invention provides a method for dynamic liability delineation of cross-border pollution in large-scale complex river network areas. This method solves important technical problems in large-scale plain river network areas, such as unclear boundaries of cross-border pollution control and difficulties in delineating liability in each area due to the complex topology of the water system, the dispersed nature of pollution sources, and dynamic emissions.
[0036] In a specific embodiment, taking a typical plain river network area as an example: Identification of the scope of responsibility for cross-boundary sections: To ensure the effectiveness and targeting of pollution control within the Dachi River network area, the main approach to delineating the comprehensive water environment management area of the cross-boundary river network area is as follows: By combining the one-dimensional water quality transport model of cross-boundary rivers and the cross-boundary section influence weight model, and through continuous trial calculations, the influence range within 12km above the cross-boundary section (administrative boundary) is determined as the influence range with an influence weight of 60%-80%; by superimposing the main influence ranges of water quality from multiple cross-boundary river systems, the effective management scope of this cross-boundary pollution responsibility delineation is obtained.
[0037] Spatial unit division and dynamic calculation of excessive emissions: Based on the water system connectivity, administrative boundaries and the results of the existing local pollution control unit division, the effective responsibility delineation area is refined into 31 control units. The daily dynamic pollution carrying capacity of each unit is calculated, which serves as the daily pollution emission allowance benchmark. Based on this, combined with dynamic pollution discharge monitoring and statistical data, the annual pollution carrying capacity and emission load are calculated on a daily basis, and the daily COD, NH3-N and TP exceedances of each unit are identified.
[0038] Pollution contribution correction and liability determination: Based on the hydrodynamic and pollution transport characteristics of important transboundary rivers within each unit, the daily excess emissions of each unit are corrected and normalized. The effective daily contribution weight of each unit to transboundary pollution COD, NH3-N and TP is calculated, and the liability is classified and determined to carry out targeted pollution control.
[0039] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.
[0040] This invention is described with reference to flowchart illustrations and / or block diagrams of systems, 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 1 A device that provides the functions specified in one or more boxes.
[0041] 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.
[0042] 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.
[0043] Although preferred embodiments of the invention have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including both the preferred embodiments and all changes and modifications falling within the scope of the invention.
[0044] Obviously, those skilled in the art can make various modifications and variations to this invention without departing from its spirit and scope. Therefore, if these modifications and variations fall within the scope of the claims of this invention and their equivalents, this invention also intends to include these modifications and variations.
Claims
1. A method for dynamic liability delineation of transboundary pollution in large-scale complex river network areas, characterized in that, Includes the following steps: Based on the hydrodynamic characteristics of important transboundary rivers in the plain river network area, the distribution and emission characteristics of pollution sources along the river, and through calculation using a one-dimensional steady-state water quality model along the river, a trial algorithm was used to determine the district and county-level administrative regions whose pollution impact weights on national and provincial control transboundary sections are within a set range as effective responsibility delineation areas for transboundary pollution. The identified effective responsibility delineation area for cross-border pollution is used as the system boundary. The county and district-level administrative units within the effective responsibility delineation area for cross-border pollution are subdivided into several pollution control spatial units of different scales within a set range. For the multiple pollution control spatial units obtained by division, based on the river distribution, river hydrodynamic dynamics and water quality control requirements within each pollution control spatial unit, the dynamic allowable discharge of each pollution control spatial unit that can reflect the dynamic pollution carrying capacity of the river is calculated. Based on the dynamic allowable emissions of each control space unit, the dynamic pollution exceedance emissions of each pollution control space unit that exceed the benchmark are determined as the benchmark for liability determination. The influence coefficient of pollution discharge on the transport process of water quality at cross-boundary sections within each pollution control spatial unit is calculated using a one-dimensional steady-state water quality model. This coefficient is then used as a correction coefficient to correct the impact of dynamic pollution exceeding standards in each pollution control spatial unit and to determine the effective dynamic contribution weight of each pollution control spatial unit to cross-boundary pollution. This allows for a refined classification of the control levels of each pollution control spatial unit.
2. The method for quantitative control and liability determination of transboundary pollution in complex river network areas according to claim 1, characterized in that, The iterative one-dimensional steady-state water quality model along the river is determined based on the river length and calculation accuracy, and takes the following form: ； ； ...； ； in, This refers to the concentration of the upstream water in the first section; This refers to the upstream water flow rate of the first section; This refers to the concentration of the downstream effluent from the first stage. This refers to the first stage of water flow rate; The concentration of pollutants in the first tributary or generalized discharge outlet; This represents the flow rate of the first tributary or generalized outlet; the process is iterated and repeated in this way. The pollutant concentration at the tributary or generalized discharge outlet N in the N+1th segment; This refers to the flow rate into a tributary or a generalized discharge outlet. The pollutant concentration at the tributary or generalized discharge outlet N in the N+1th segment; The flow rate into a tributary or a generalized discharge outlet; The concentration of pollutants effluent in the (N+1)th segment; Outflow rate; The degradation coefficient of the target pollutant; The average flow velocity of the river channel; The distance is the flow distance of the N+1 segment of the river channel.
3. The method for quantitative control and liability determination of transboundary pollution in complex river network areas according to claim 2, characterized in that, The calculation process for the pollution impact weight is as follows: Based on a one-dimensional steady-state water quality model along the river, the water quality concentration of major transboundary rivers at the inflow boundaries of the trial county and district administrative regions is used as the inflow pollution concentration. Based on the hydrodynamic characteristics of the transboundary river and the pollution input results of all tributaries and generalized discharge outlets along the river, the water quality concentration at the national and provincial control transboundary sections under the influence of pollution sources along the river is calculated iteratively. ; Simultaneously, the water quality concentration of the blank control group was calculated in the iterative calculation, ignoring the impact of pollution along the route. ; Obtain pollution impact weight expression: 。 4. The method for quantitative control and liability determination of transboundary pollution in complex river network areas according to claim 1, characterized in that, The process for calculating the dynamic allowable discharge of each pollution control spatial unit that reflects the dynamic pollution carrying capacity of the river is as follows: The characteristic parameters of each pollution control spatial unit are obtained, including flow rate, flow velocity, actual pollutant concentration and target pollutant concentration; The characteristic parameters are input into the pollution holding capacity calculation formula to calculate the pollution holding capacity, where the pollution holding capacity calculation formula is in the form of: ； in, The dynamic contamination-carrying capacity of space unit i on day k. The number of sewage-receiving river sections within a spatial unit that have a water environment function positioning; Indicates the first unit within the cell The first section of the river The average flow rate on the k-th day of the day; The first one determined according to the water environment function positioning The first effective control unit Target concentration of pollutants at the inflow section of each river segment; For the first The first effective control unit Target concentration of pollutants at the outflow section of each river segment; Indicates the first The first effective control unit The first-level comprehensive degradation coefficient of each river section; For the first The first effective control unit The length of each river section; For the first The first effective control unit The average daily flow velocity on day k of each river segment.
5. The method for quantitative control and liability determination of transboundary pollution in complex river network areas according to claim 1, characterized in that, The process for obtaining dynamic pollution exceedance data for each pollution control spatial unit is as follows: Determine the number of point sources, land use area, and pollution loss intensity within each effective control unit, and calculate the total pollution load: ； in, For the first Total pollution emission load of each spatial unit on day k; For the first Within the effective spatial unit, the first The average daily emission load of individual point sources is mostly average. and The first The area of land use type h within each effective spatial unit and its pollution loss intensity on day k; Calculate the pollution overload of each spatial unit. : 。 6. The method for quantitative control and liability determination of transboundary pollution in complex river network areas according to claim 1, characterized in that, The process of correcting the impact of dynamic pollution exceeding emission standards in each pollution control spatial unit and determining the effective dynamic contribution weight of each pollution control spatial unit to transboundary pollution is as follows: Calculate the sum of pollution overloads for all effective control units; Calculate the correction factor for pollution emission process in spatial units. The calculation formula is as follows: ； In the formula, k j To control the degradation coefficient of transboundary river j, Q j To control the flow of the transboundary river j; j v represents the distance that the transboundary river j flows from its entry point to the transboundary section within the control unit; j The average flow velocity of the transboundary river j from its entry point to the transboundary section; The pollution excess load of each effective control unit is normalized based on the sum of the pollution excess loads, and the dynamic pollution responsibility contribution rate is obtained by combining the pollution transport process correction coefficients calculated for each unit. : 。 7. The method for responsibility delineation and adaptive management of dynamic pollution carrying capacity calculation in plain river network areas according to claim 6, characterized in that, The process of classifying the control levels of each pollution control spatial unit in a refined manner is as follows: If the dynamic control index is not less than the first threshold, it is classified as a key control area; If the dynamic control index is less than the first threshold and not less than the second threshold, it is classified as a general control zone. If the dynamic control index is less than the second threshold, it is classified as a light control zone; If the dynamic control index is less than 0, it is classified as an area with no impact.