CN110335172B - Watershed water environment capacity distribution method and device, electronic equipment and storage medium - Google Patents

Abstract

The application discloses a watershed water environment capacity distribution method and device, electronic equipment and a storage medium, and belongs to the technical field of water environment monitoring. The method comprises the following steps: acquiring a watershed to be distributed, which needs to perform water environment capacity distribution; acquiring each control unit related to the watershed to be distributed and river reach corresponding to each control unit one by one; acquiring the maximum water environment capacity allowed to be discharged by each river reach; determining the water environment capacity allowed by the control unit corresponding to each river reach based on the maximum water environment capacity allowed to be discharged by each river reach; and correcting the water environment capacity allowed by the control unit corresponding to each river reach. In the embodiment of the application, the control units in the distribution domain are divided, the catchment range of each control unit is determined, then the maximum water environment capacity allowed to be discharged of each river reach is obtained, and then the water environment capacity allowed by the control unit corresponding to each river reach is determined and corrected, so that the objectivity and fairness of the distribution result are guaranteed.

Description

Watershed water environment capacity distribution method and device, electronic equipment and storage medium

Technical Field

The application belongs to the technical field of water environment monitoring, and particularly relates to a watershed water environment capacity distribution method and device, electronic equipment and a storage medium.

Background

At present, the idea of controlling the total amount of water pollutants in China is implemented, and the reasonable distribution of water environment capacity (particularly the maximum load of water pollutants under the condition of meeting the requirement of water environment quality, and is also called as the water load or pollutant carrying capacity) is the key of controlling the total amount of water pollutants. The pollutant discharge amount reduction mainly comprises the steps of distributing quantitative discharge limits to different control units according to a certain method, and then reducing the pollutant discharge amount according to the current pollution discharge situation in the control units.

At present, the distribution method of the water environment capacity mainly comprises an equal proportion distribution method, a fair interval distribution method, a linear programming method and the like. The equal proportion distribution method is simple to operate, distribution indexes are single, and unfairness of multi-row and multi-division is easily caused due to lack of consideration on regional development. The fair interval allocation method is the optimization of the unfairness of the equal proportion allocation method, but still does not depart from the defect of equal proportion allocation of a single allocation index. The linear programming method always reflects the characteristic difference of different distribution objects from the aspect of economic benefit, but the unfairness phenomenon of distribution results is easily caused by objectively emphasizing the optimization of the economic benefit.

Disclosure of Invention

In view of this, an object of the present application is to provide a watershed water environment capacity allocation method, device, electronic device and storage medium, so as to improve the unfairness of the allocation result easily caused by the allocation method, and further ensure objectivity and fairness of the allocation result.

The embodiment of the application is realized as follows:

in a first aspect, an embodiment of the present application provides a watershed water environment capacity allocation method, including acquiring a watershed to be allocated, where water environment capacity allocation is required; acquiring each control unit related to the watershed to be distributed and river reach corresponding to each control unit one by one; acquiring the maximum water environment capacity allowed to be discharged by each river reach; determining the water environment capacity allowed by the control unit corresponding to each river reach based on the maximum water environment capacity allowed to be discharged by each river reach; and correcting the water environment capacity allowed by the control unit corresponding to each river reach. In the embodiment of the application, the control units in the distribution domain are divided, the catchment range of each control unit is determined, then the maximum water environment capacity allowed to be discharged by each river reach corresponding to the control units one by one is obtained, then the water environment capacity allowed by the control unit corresponding to each river reach is determined, the determined water environment capacity is further corrected, equal proportion distribution is not performed, and the objectivity and fairness of the distribution result are guaranteed.

With reference to one possible implementation manner of the embodiment of the first aspect, the modifying the allowable water environment capacity of the control unit corresponding to each river reach includes: determining that the value of the water environment capacity allowed by the control unit corresponding to the upstream river reach is unchanged; acquiring the retention coefficients of pollutants of other river sections except the uppermost stream section; determining that the value of the water environment capacity allowed by a control unit of the river reach with the pollutant retention coefficient not greater than a preset threshold value is not changed; and correcting the water environment capacity allowed by the control unit of the target river reach with the pollutant retention coefficient larger than the preset threshold value. In the embodiment of the application, during correction, the influence of the pollutants of the upstream river reach on the downstream river reach is taken into account, the water environment capacity allowed by the control unit of the target river reach of which the pollutant retention coefficient is larger than the preset threshold value in other river reach except the upstream river reach is corrected, and the influence of the pollutants of the upstream river reach on the downstream river reach is taken into account, so that the accuracy and the objectivity of the distribution result are further ensured.

With reference to a possible implementation manner of the embodiment of the first aspect, the modifying the water environmental capacity allowed by the control unit of the target river reach with the pollutant retention coefficient greater than the preset threshold value includes: acquiring all upstream river reach positioned at the upstream of the target river reach; acquiring the pollutant contribution rate of each upstream river reach to the target river reach; obtaining the sum of the contribution of all the upstream river reach to the water environment capacity of the target river reach based on the maximum water environment capacity of each upstream river reach and the pollutant contribution rate of each upstream river reach to the target river reach; and correcting the water environment capacity allowed by the control unit corresponding to the target river reach as the contribution sum of the water environment capacity. In the embodiment of the application, when the water environment capacity allowed by the control unit corresponding to the target river reach is corrected, the contribution rate of all upstream river reaches positioned at the upstream of the target river reach to the pollutants of the target river reach is obtained, then the water environment capacity contribution sum of all upstream river reaches to the target river reach is obtained based on the maximum water environment capacity of each upstream river reach and the pollutant contribution rate of each upstream river reach to the target river reach, finally the water environment capacity allowed by the control unit corresponding to the target river reach is corrected into the water environment capacity contribution sum, the conversion of the pollutants in the migration process is considered in distribution, and the scientificity of the accounting result of the contribution amount is ensured.

With reference to one possible implementation manner of the embodiment of the first aspect, before modifying the sum of the contribution of the aquatic environment capacity allowed by the control unit corresponding to the target river reach to the aquatic environment capacity, the method further includes: determining the contribution of the aquatic environment capacity is less than the aquatic environment capacity allowed by the target river reach. In the embodiment of the application, before the water environment capacity allowed by the control unit corresponding to the corrected target river reach is the sum of the water environment capacity contributions, the sum of the water environment capacity contributions is determined to be smaller than the water environment capacity allowed by the target river reach, so that the accuracy of the corrected result is ensured.

With reference to one possible implementation manner of the embodiment of the first aspect, before obtaining the pollutant contribution rate of each upstream river reach to the target river reach, the method further includes: the rate of contribution of pollutants from each upstream stretch to each downstream stretch located downstream from itself is determined. In the embodiment of the application, the pollutant contribution rate of each upstream river reach to each downstream river reach downstream of the upstream river reach is determined in advance, so that corresponding data can be directly acquired when the demand exists, and the cost is saved and the processing efficiency is accelerated.

With reference to one possible implementation of the embodiment of the first aspect, determining a pollutant contribution rate of each upstream stretch to each downstream stretch located downstream of the upstream stretch comprises: establishing a transition probability matrix according to the upstream and downstream relations of the river reach corresponding to each control unit; establishing a retention coefficient matrix according to the pollutant inflow amount and the pollutant output amount of the river reach corresponding to each control unit and the upstream-downstream relation of each river reach; establishing a pollutant transfer coefficient matrix according to the transfer probability matrix and the retention coefficient matrix; and determining the pollutant contribution rate of each upstream river reach to each downstream river reach positioned downstream of the upstream river reach based on the pollutant transfer coefficient matrix. In the embodiment of the application, each control unit is set to be a discrete state variable by introducing a Markov chain method, upstream and downstream river networks corresponding to the control units are represented by using a topological relation, a migration and transportation path of river pollutants in an area is definitely researched, retention coefficients of the pollutants in each river reach are definitely determined, a retention coefficient matrix is established, then the contribution rate of each upstream river reach to the pollutants in each downstream river reach located at the downstream of each upstream river reach is determined based on the pollutant transfer coefficient matrix established by the transfer probability matrix and the retention coefficient matrix, and the scientificity and the accuracy of the accounting result of the contribution amount are ensured.

With reference to a possible implementation manner of the embodiment of the first aspect, establishing a retention coefficient matrix according to the inflow amount of pollutants into the river reach, the output amount of pollutants of the river reach corresponding to each control unit, and the upstream-downstream relationship of each river reach includes: calculating pollutant retention coefficients corresponding to all the river reach according to the pollutant inflow amount and the pollutant output amount of the river reach corresponding to all the control units and the upstream and downstream relations of all the river reach; and establishing a retention coefficient matrix based on the respective corresponding pollutant retention coefficients of the river reach and the upstream and downstream relations of the river reach.

With reference to one possible implementation of the embodiment of the first aspect, before obtaining the contaminant retention coefficients of the river reach other than the most upstream river reach, the method further includes: and calculating the pollutant retention coefficient corresponding to each river reach. In the embodiment of the application, the pollutant retention coefficient corresponding to each river reach is calculated in advance, so that corresponding data can be directly acquired when the requirement exists, and the cost is saved and the processing efficiency is accelerated.

With reference to one possible implementation manner of the embodiment of the first aspect, calculating the pollutant retention coefficient corresponding to each river reach includes: calculating the pollutant retention coefficient corresponding to each river reach according to the pollutant inflow amount, the pollutant output amount and the position of each river reach, wherein the calculation formula of the pollutant retention coefficient is r_i＝(Load_in-Load_out)/Load_in，r_iIs the corresponding pollutant retention coefficient of river reach i, Load_inThe pollutant Load capacity corresponding to the river reach i is named as Loadin plus Load (i-1) out, subi is the pollutant inflow capacity corresponding to the river reach i, Load (i-1) out is the pollutant output capacity corresponding to the directly upstream i-1 of the river reach i, and Loadout is the pollutant output capacity corresponding to the river reach i. In the embodiment of the application, when the retention coefficient of the pollutants is calculated, the influence of the directly upstream on the downstream river reach is fully considered, and the reasonability of the calculation result is ensured.

With reference to a possible implementation manner of the embodiment of the first aspect, the obtaining each control unit related to the watershed to be distributed and a river reach corresponding to each control unit one to one includes: and dividing the control units in the to-be-distributed flow domain based on a home management principle, and determining the river reach corresponding to each control unit. In the embodiment of the application, the control units in the to-be-distributed flow domain are divided based on the home management principle, and the river reach corresponding to each control unit is determined, so that the operation is simple, and environment management decision makers can conveniently make decisions.

With reference to one possible implementation manner of the embodiment of the first aspect, after the correcting the allowable water environment capacity of the control unit corresponding to each river reach, the method further includes: and determining the water environment capacity of the target pollution source corresponding to each control unit according to the corrected water environment capacity allowed by each control unit. In the embodiment of the application, the corrected water environment capacity of the target pollution source in each control unit is further distributed to control the specific water pollutant discharge amount of various pollution sources.

With reference to one possible implementation manner of the embodiment of the first aspect, determining, according to the corrected allowable water environment capacity of each control unit, the water environment capacity of the target pollution source corresponding to each control unit includes: acquiring the pollutant river entering amount of each preset pollution source in each control unit; and determining the water environment capacity of the target pollution source corresponding to each control unit according to the corrected water environment capacity allowed by each control unit and the acquired river inflow amount of the pollutants of each preset pollution source corresponding to each control unit. In the embodiment of the application, when the water environment capacity of the target pollution source in a certain control unit is further distributed, the water environment capacity of the target pollution source is determined according to the water environment capacity allowed by the control unit after the water environment capacity of each preset pollution source in the control unit is obtained and the water environment capacity of each preset pollution source in the control unit is corrected, so that the contribution of each sewage discharge unit is fully considered, and the rationality of the distribution result is ensured.

With reference to a possible implementation manner of the embodiment of the first aspect, the acquiring the pollutant river entering amount of each preset pollution source in each control unit includes: acquiring pollutant discharge amount of each preset pollution source in each control unit and river entering coefficients corresponding to each preset pollution source; and obtaining the river entering amount of the pollutants of each preset pollution source in each control unit according to the pollutant discharge amount of each preset pollution source in each control unit and the river entering coefficient corresponding to each preset pollution source. In the embodiment of the application, the pollutant discharge amount of each preset pollution source in the control unit and the river entering coefficient corresponding to each preset pollution source are obtained, so that the corresponding pollutant river entering amount of the preset pollution sources is obtained, and the controllability and the reasonability of the pollutant river entering amount are guaranteed.

With reference to one possible implementation manner of the embodiment of the first aspect, the target pollution source is at least one of a rural direct-discharge living source, an agricultural planting source, a livestock breeding source, an industrial point source and a direct-discharge urban living source. In the embodiment of the application, when the water environment capacity of the target pollution source in the control unit is further distributed, a main stream pollution source with a large influence result is considered, and the reasonability of the distribution result is further ensured.

In a second aspect, an embodiment of the present application further provides a watershed water environment capacity distribution device, including: the device comprises a first acquisition module, a second acquisition module, a third acquisition module, a first determination module and a correction module; the system comprises a first acquisition module, a first distribution module and a second acquisition module, wherein the first acquisition module is used for acquiring a watershed to be distributed, which needs to distribute water environment capacity; the second acquisition module is used for acquiring each control unit related to the watershed to be distributed and river reach corresponding to each control unit one by one; the third acquisition module is used for acquiring the maximum water environment capacity allowed to be discharged by each river reach; the first determining module is used for determining the water environment capacity allowed by the control unit corresponding to each river reach based on the maximum water environment capacity allowed to be discharged by each river reach; and the correction module is used for correcting the water environment capacity allowed by the control unit corresponding to each river reach.

With reference to a possible implementation manner of the embodiment of the second aspect, the modification module is specifically configured to: determining that the value of the water environment capacity allowed by the control unit corresponding to the upstream river reach is unchanged; acquiring the retention coefficients of pollutants of other river sections except the uppermost stream section; determining that the value of the water environment capacity allowed by a control unit of the river reach with the pollutant retention coefficient not greater than a preset threshold value is not changed; and correcting the water environment capacity allowed by the control unit of the target river reach with the pollutant retention coefficient larger than the preset threshold value.

With reference to a possible implementation manner of the embodiment of the second aspect, the modification module is specifically configured to: acquiring all upstream river reach positioned at the upstream of the target river reach; acquiring the pollutant contribution rate of each upstream river reach to the target river reach; obtaining the sum of the contribution of all the upstream river reach to the water environment capacity of the target river reach based on the maximum water environment capacity of each upstream river reach and the pollutant contribution rate of each upstream river reach to the target river reach; and correcting the water environment capacity allowed by the control unit corresponding to the target river reach as the contribution sum of the water environment capacity.

With reference to one possible implementation manner of the embodiment of the second aspect, the method further includes: and the second determination module is used for determining that the contribution sum of the water environment capacity is smaller than the water environment capacity allowed by the target river reach before the correction module is used for correcting the water environment capacity allowed by the control unit corresponding to the target river reach to be the contribution sum of the water environment capacity.

With reference to one possible implementation manner of the embodiment of the second aspect, the method further includes: and the third determining module is used for determining the contribution sum of the water environmental capacity to be less than the allowable water environmental capacity of the target river reach before the correcting module is used for acquiring the pollutant contribution rate of each upstream river reach to the target river reach.

With reference to a possible implementation manner of the embodiment of the second aspect, the third determining module is specifically configured to: establishing a transition probability matrix according to the upstream and downstream relations of the river reach corresponding to each control unit; establishing a retention coefficient matrix according to the pollutant inflow amount and the pollutant output amount of the river reach corresponding to each control unit and the upstream-downstream relation of each river reach; establishing a pollutant transfer coefficient matrix according to the transfer probability matrix and the retention coefficient matrix; and determining the pollutant contribution rate of each upstream river reach to each downstream river reach positioned downstream of the upstream river reach based on the pollutant transfer coefficient matrix.

With reference to a possible implementation manner of the embodiment of the second aspect, the third determining module is specifically configured to: calculating pollutant retention coefficients corresponding to all the river reach according to the pollutant inflow amount and the pollutant output amount of the river reach corresponding to all the control units and the upstream and downstream relations of all the river reach; and establishing a retention coefficient matrix based on the respective corresponding pollutant retention coefficients of the river reach and the upstream and downstream relations of the river reach.

With reference to one possible implementation manner of the embodiment of the second aspect, the method further includes: and the calculation module is used for calculating the pollutant retention coefficient corresponding to each river reach before the correction module is used for acquiring the pollutant retention coefficients of other river reaches except the most upstream river reach.

With reference to a possible implementation manner of the embodiment of the second aspect, the calculating module is specifically configured to: calculating the pollutant retention coefficient corresponding to each river reach according to the pollutant inflow amount, the pollutant output amount and the position of each river reach, wherein the calculation formula of the pollutant retention coefficient is r_i＝(Load_in-Load_out)/Load_in，r_iIs the corresponding pollutant retention coefficient of river reach i, Load_inThe pollutant Load capacity corresponding to the river reach i is named as Loadin plus Load (i-1) out, subi is the pollutant inflow capacity corresponding to the river reach i, Load (i-1) out is the pollutant output capacity corresponding to the directly upstream i-1 of the river reach i, and Loadout is the pollutant output capacity corresponding to the river reach i.

With reference to a possible implementation manner of the embodiment of the second aspect, the second obtaining module is specifically configured to: and dividing the control units in the to-be-distributed flow domain based on a home management principle, and determining the river reach corresponding to each control unit.

With reference to one possible implementation manner of the embodiment of the second aspect, the method further includes: and the fourth determining module is used for determining the water environment capacity of the target pollution source corresponding to each control unit according to the corrected water environment capacity allowed by each control unit after the correction determining module corrects the water environment capacity allowed by the control unit corresponding to each river reach.

With reference to a possible implementation manner of the embodiment of the second aspect, the fourth determining module is specifically configured to: acquiring the pollutant river entering amount of each preset pollution source in each control unit; and determining the water environment capacity of the target pollution source corresponding to each control unit according to the water environment capacity allowed by each control unit and the acquired river inflow amount of the pollutants of each preset pollution source corresponding to each control unit.

With reference to a possible implementation manner of the embodiment of the second aspect, the fourth determining module is specifically configured to: acquiring pollutant discharge amount of each preset pollution source in each control unit and river entering coefficients corresponding to each preset pollution source; and obtaining the river entering amount of the pollutants of each preset pollution source in each control unit according to the pollutant discharge amount of each preset pollution source in each control unit and the river entering coefficient corresponding to each preset pollution source.

In combination with one possible implementation manner of the embodiment of the second aspect, the target pollution source is at least one of a rural direct-discharge living source, an agricultural planting source, a livestock breeding source, an industrial point source and a direct-discharge urban living source.

In a third aspect, an embodiment of the present application further provides an electronic device, including: a memory and a processor, the memory and the processor connected; the memory is used for storing programs; the processor is configured to invoke a program stored in the memory to perform the method of the first aspect embodiment and/or any possible implementation manner of the first aspect embodiment.

In a fourth aspect, embodiments of the present application further provide a storage medium, on which a computer program is stored, where the computer program is executed by a computer to perform the method provided in the foregoing first aspect and/or any possible implementation manner of the first aspect.

Additional features and advantages of the application will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by the practice of the embodiments of the application. The objectives and other advantages of the application may be realized and attained by the structure particularly pointed out in the written description and drawings.

Drawings

In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and it is obvious for those skilled in the art to obtain other drawings without creative efforts. The foregoing and other objects, features and advantages of the application will be apparent from the accompanying drawings. Like reference numerals refer to like parts throughout the drawings. The drawings are not intended to be to scale as practical, emphasis instead being placed upon illustrating the subject matter of the present application.

Fig. 1 shows a flow chart of a watershed water environment capacity allocation method provided by an embodiment of the present application.

Fig. 2 shows a schematic diagram of a catchment space topological relation of each control unit in a flow domain to be distributed according to an embodiment of the present application.

Fig. 3 is a schematic flow chart of another watershed water environment capacity allocation method provided by an embodiment of the application.

Fig. 4 shows a schematic flow chart of another watershed water environment capacity allocation method provided by the embodiment of the application.

Fig. 5 is a schematic flow chart of another watershed water environment capacity allocation method provided by an embodiment of the application.

Fig. 6 shows a schematic block diagram of a watershed water environment capacity distribution device provided by an embodiment of the application.

Fig. 7 shows a schematic structural diagram of an electronic device provided in an embodiment of the present application.

Detailed Description

The technical solutions in the embodiments of the present application will be described below with reference to the drawings in the embodiments of the present application.

It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures. Meanwhile, relational terms such as "first," "second," and the like may be used solely in the description herein to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other identical elements in a process, method, article, or apparatus that comprises the element.

Further, the term "and/or" in the present application is only one kind of association relationship describing the associated object, and means that three kinds of relationships may exist, for example, a and/or B may mean: a exists alone, A and B exist simultaneously, and B exists alone.

First embodiment

Referring to fig. 1, steps included in a method for allocating a watershed water environment capacity according to an embodiment of the present application will be described with reference to fig. 1.

Step S101: and acquiring the watershed to be distributed, which needs to distribute the water environment capacity.

Step S102: and acquiring each control unit related to the watershed to be distributed and river reach corresponding to each control unit one by one.

And after the watershed to be distributed is obtained, obtaining each control unit related to the watershed to be distributed and river reach corresponding to each control unit one by one. For example, the control units in the to-be-distributed watershed are divided based on the home management principle, and the river reach corresponding to each control unit is determined, so that each control unit related to the to-be-distributed watershed and the river reach corresponding to each control unit one to one can be obtained. For example, according to the home management principle, taking a street town as each administrative subject, dividing the control units in the domain to be distributed, and comprehensively considering the administrative boundary and the river catchment range to determine the river reach corresponding to each control unit.

As an embodiment, the obtained control unit and the corresponding river reach may be obtained from a database, that is, stored in the database in advance, and corresponding data may be obtained from the database when there is a need. For example, the control units and the corresponding river reach in the hot watershed are determined to be divided in advance, and the obtained data is stored in the database, that is, the database stores the control units and the corresponding river reach data corresponding to each of the plurality of hot watersheds, for example, the control units and the corresponding river reach data of the divided hot watershed a, the control units and the corresponding river reach data of the divided hot watershed B, and the like.

When the control units in the watershed to be distributed are divided and the river reach corresponding to each control unit is determined, hydrologic and water quality data of the watershed to be distributed need to be collected and analyzed, a basic database is established, a Geographical Information System (GIS) is used as a support to extract and research a river network water system diagram of the watershed and distribute the river network water system diagram to each control unit, and the river reach corresponding to each control unit is determined.

Step S103: and acquiring the maximum water environment capacity allowed to be discharged in each river reach.

And after the control unit of the river reach to be distributed and the corresponding river reach are determined, acquiring the maximum water environment capacity allowed to be discharged by each river reach, and accounting the maximum water environment capacity of the river reach in the river reach based on the existing hydrological water quality data. According to a differential equation of a one-dimensional steady-state attenuation law of river pollutants:

and is composed of

Obtaining by solution: c ═ C₀e-Kxu (1), wherein C is the in-process contaminant concentration mg/L; u is the average flow velocity m/s of the river section; x is the distance km along the way; k is the comprehensive degradation coefficient d-1 of pollutants; c0 is the concentration of water pollutant in mg/L at the upstream section.

Considering the influence of river drain pollutant discharge to river water quality, bring the drain pollutant discharge into the calculation to the assumed complete required on-way distance of mixing of the pollutant that discharges into the river course is less than river course length far away, and the pollutant just can be in the river course cross section misce bene in short time, then section pollutant concentration is under the river reach:

let C_L＝C_SThe water environment capacity of each river reach is obtained by the formula (1) and the formula (2):

W_j＝[C_s-C₀exp(-kx/u)]exp(kx_i/u)Q (3)

wherein, C_SThe target water quality mg/L of the cross section is controlled for researching the river reach, and can be regarded as a fixed value and is the water quality expected to reach the standard C_LIs the concentration mg/L of pollutants on the lower section of the river reach; c₀The concentration of the pollutants on the upper section of the river reach is mg/L; w_jRepresenting the maximum water environment capacity t allowed to be discharged by the river reach j; q represents the generalized sewage discharge amount m of the sewage discharge outlet³. Wherein, C_SCan be regarded as a constant value, C₀And Q varies from subject to subject.

The research basin control section can be generally arranged at the positions of a water quality monitoring port, a water intake, an upstream and downstream key node or a basin outlet and the like in the basin. And (4) substituting the monitored data into the formula (3) to obtain the maximum water environment capacity allowed to be discharged in each river reach. Alternatively, the data required to calculate the maximum allowable aquatic capacity for each river reach may be obtained by field measurements.

Wherein, when the drain is generalized, a plurality of adjacent drains in the same control unit are generalized into a centralized drain. The calculation method of the pollutant concentration of the generalized sewage draining exit can be that the pollutant concentrations of a plurality of sewage draining exits in the same generalized sewage draining exit are accumulated to obtain the pollutant concentration of the generalized sewage draining exit. For example, 3 drain are generalized to be a drain, and then the pollutant concentration of this generalized drain equals the accumulation of the pollutant concentration of these 3 drains, obtains the pollutant concentration of generalized drain.

Step S104: and determining the water environment capacity allowed by the control unit corresponding to each river reach based on the maximum water environment capacity allowed to be discharged by each river reach.

And (4) determining the water environment capacity allowed by the control unit corresponding to each river reach based on the maximum water environment capacity allowed to be discharged by each river reach after the maximum water environment capacity allowed to be discharged by each river reach is obtained based on the formula (3).

For convenience of understanding, the following description will be given by taking the control unit of a certain basin to be distributed shown in fig. 2 as an example, and assuming that the control units of a certain basin to be distributed are divided into 6, the catchment ranges of the 6 control units are analyzed, so as to obtain the catchment spatial topological relation shown in fig. 2. Wherein the arrows in the figure indicate the direction of the bus-flow. The principle of calculating the maximum allowable discharge water environment capacity of each river reach based on the above equation (3) is the same, and for easy understanding, the following description will be given by taking W1 as an example, and the standard control section target water quality C will be described_sAnd the average flow velocity u of the river section of the river reach 1, the distance x along the way, the sewage discharge quantity Q of the generalized sewage discharge outlet and the pollutant concentration C of the section on the river reach are measured₀The maximum water environment capacity allowed to be discharged by W1 can be obtained by substituting the formula. Wherein, because the river reach 1, the river reach 2 and the river reach 6 are positioned at the most upstream, the concentration C of the pollutant on the section of the river reach₀I.e. its own concentration of contaminants, and for the downstream sections 3, 4, 5, C₀I.e. the concentration of the contaminant at the outlet of the upstream section (i.e. at the inlet itself). Obtaining the maximum water environment capacity W allowed to be discharged at each river section_jThen, the water environment capacity allowed by the control unit corresponding to each river reach, for example, the water environment capacity W allowed by the control unit corresponding to the river reach 1, is determined₀May be equal to W1 multiplied by a factor, such as W₀T × W1, where the coefficient t ranges from 0.9 to 1, for example, when t is 1, the allowable water environment capacity of the control unit is equal to the maximum allowable water environment capacity of the corresponding river reach, for example, the allowable water environment capacity of the control unit 1 is equal to W1.

Second embodiment

The present embodiment is different from the first embodiment in that, as shown in fig. 3, after step S104, the method further includes: and correcting the water environment capacity allowed by each control unit. In the correction, it is determined that the value of the water environment capacity allowed by the control unit corresponding to the most upstream river reach is not changed, for example, the values of the water environment capacity allowed by the control units corresponding to the river reach 1, the river reach 2, and the river reach 6 in fig. 2 are maintained. When the water environment capacity allowed by the control unit corresponding to other river reach except the most upstream river reach is corrected, the process may be as follows: acquiring the retention coefficients of pollutants of other river sections except the uppermost stream section; the water environment capacity allowed by the control unit of the target river reach with a pollutant retention coefficient greater than a preset threshold (as an embodiment, the preset threshold may be 0, or a positive number close to 0, such as 0.1) is corrected. With reference to fig. 2, pollutant retention coefficients of the river reach 3, the river reach 4 and the river reach 5 are obtained, water environment capacity allowed by a control unit of a target river reach of which the pollutant retention coefficients are larger than a preset threshold value is corrected, and if the pollutant retention coefficients of the river reach 3, the river reach 4 and the river reach 5 are all larger than the preset threshold value, the river reach 3, the river reach 4 and the river reach 5 are all target river reach; similarly, if it is assumed that the pollutant retention coefficient of only the river reach 5 is greater than the preset threshold, only the river reach 5 is the target river reach. And the value of the water environment capacity allowed by the control unit of the river reach except the upstream river reach, of which the pollutant retention coefficient is not more than the preset threshold value, is unchanged.

When the water environment capacity allowed by the control unit of the target river reach with the pollutant retention coefficient larger than the preset threshold is corrected, the water environment capacity may be: acquiring all upstream river reach positioned at the upstream of the target river reach; acquiring the pollutant contribution rate of each upstream river reach to the target river reach; obtaining the sum of the contribution of all the upstream river reach to the water environment capacity of the target river reach based on the maximum water environment capacity of each upstream river reach and the pollutant contribution rate of each upstream river reach to the target river reach; and correcting the water environment capacity allowed by the control unit corresponding to the target river reach as the contribution sum of the water environment capacity. For convenience of understanding, the description will be given with reference to fig. 2, taking the target river reach as the river reach 3 as an example, and the river reaches located upstream of the target river reach 3 are the river reach 1 and the river reach 2, the pollutant contribution rate of the river reach 1 to the river reach 3 is obtained, and the pollutant contribution rate of the river reach 2 to the river reach 3 is obtained, if G13, if G23, based on the maximum aquatic environment capacity W1 of the river reach 1 and the pollutant contribution rate G13 to the target river reach 3, and the maximum water environment capacity W2 of the river reach 2 and the pollutant contribution rate G23 to the target river reach 3, the sum of the water environment capacity contribution to the target river reach 3 is obtained as W1G 13+ W2G 23, the water environment capacity allowed by the control unit corresponding to the target river reach 3 is corrected to be the sum of the water environment capacity contribution, that is, the allowable water environment capacity of the target river reach 3 before the correction is W3, and after the correction, it is W1 × G13+ W2 × G23.

Taking the target river reach as the river reach 4 as an example, the river reach located at the upstream of the target river reach 4 is the river reach 1, the river reach 2, the river reach 3 and the river reach 6, the pollutant contribution rate of the river reach 1 to the river reach 4 is obtained, the pollutant contribution rate of the river reach 2 to the river reach 4 is obtained under the assumption of G14, the pollutant contribution rate of the river reach 3 to the river reach 4 is obtained under the assumption of G24, the pollutant contribution rate of the river reach 3 to the river reach 4 is obtained under the assumption of G34, the pollutant contribution rate of the river reach 6 to the river reach 4 is obtained under the assumption of G64, the pollutant contribution rate G14 to the target river reach 4, the maximum water environment capacity W2 to the river reach 2, the pollutant contribution rate G38 to the target river reach 4, the maximum water environment capacity W3 to the target river reach 3 and the pollutant contribution rate G34 to the target river reach 4 are obtained on the basis of the maximum water environment capacity W1 of the river reach 1, the water environment capacity W6329 of the river reach 6, the pollutant contribution rate G6338 to the target river reach 4, the pollutant contribution rate G2 to the water environment G638 + 6 + 9, and the water environment capacity G639 + 6 to the target river reach 4 and the allowable water environment capacity of the control unit corresponding to the target river reach 4 is corrected to be the sum of the contributions of the water environment capacities, that is, before correction, the allowable water environment capacity of the target river reach 4 is W4, and after correction, the allowable water environment capacity is W1 × G14+ W2 × G24+ W3 × G34+ W6 × G64. When the target river reach is river reach 5, the correction process is similar to target river reach 3 or target river reach 4, and is not burdensome.

As an embodiment, when the allowable water environment capacity of the control unit corresponding to the target river reach is corrected to the sum of the water environment capacity contributions, it may be determined whether the allowable water environment capacity of the control unit corresponding to the target river reach is greater than the sum of the water environment capacity contributions, and when the sum of the water environment capacity contributions is determined to be less than the allowable water environment capacity of the target river reach, the correction is performed, for example, when the target river reach 3 is corrected, it is determined whether W1G 13+ W2G 23 is greater than W3, only when W1G 13+ W2G 23 is less than W3, the correction is performed, otherwise, the allowable water environment capacity of the target river reach 3 remains unchanged to W3, and the process of the other target river reach is similar to that.

As an embodiment, the above-mentioned obtaining of the pollutant retention coefficients of the other river segments except for the most upstream river segment may be obtained from a database, that is, the pollutant retention coefficients of the river segments under a plurality of hot watersheds are stored in the database in advance, and corresponding data may be obtained when the database is used. For example, the pollutant retention coefficients of the individual river segments shown in fig. 2 are stored. In this embodiment, before the pollutant retention coefficients of the other sections except the uppermost stream section are obtained, the pollutant retention coefficients corresponding to the respective sections need to be calculated.

When the pollutant retention coefficient corresponding to each river reach is calculated, the pollutant retention coefficient corresponding to each river reach can be calculated according to the inflow amount of the pollutant, the output amount of the pollutant and the position of the pollutant corresponding to each river reach. Wherein the calculation formula of the retention coefficient of the pollutants is r_i＝(Load_in-Load_out)/Load_in，r_iIs the corresponding pollutant retention coefficient of river reach i, Load_inThe pollutant Load capacity corresponding to the river reach i is named as Loadin plus Load (i-1) out, subi is the pollutant inflow capacity corresponding to the river reach i, Load (i-1) out is the pollutant output capacity corresponding to the directly upstream i-1 of the river reach i, and Loadout is the pollutant output capacity corresponding to the river reach i. For example, for river reach 1, the corresponding pollutant Load is equal to the corresponding pollutant inflow, i.e. Load, and for river reach 3, the pollutant Load is equal to the pollutant inflow of river reach 3+ the pollutant output of river reach 2 + the pollutant output of river reach 1, i.e. Load sub3+ Load2out + Load1 out.

Wherein the pollutant output quantity of each river reach can be obtained by measuring the pollutant output quantity of the outlet end face of each river reach.

The pollutant inflow amount of each river reach can be obtained by setting the pollutant discharge amount of the pollution source in the corresponding control unit, wherein the pollutant inflow amount is calculated by the following formula: the method comprises the following steps of (1) determining the river entering quantity of a direct-discharge living source pollutant in a rural area by using ei ═ Ns-i + Nz-i + Nx-i + Cs-i + Gs-i, wherein Ns-i is the river entering quantity of the direct-discharge living source pollutant in the rural area, Nz-i is the river entering quantity of the agricultural planting source pollutant, Nx-i is the river entering quantity of the livestock and poultry breeding source pollutant, Gs-i is the river entering quantity of the industrial point source pollutant, and Cs-i is the river entering quantity of the direct-discharge urban living source pollutant; and wherein Ns-i-Ns-out a; nz-i ═ Nz-out ═ b; nx-i ═ Nx-out ═ c; Cs-i-Cs-out d; gs-i is Gs-out, wherein Ns-out is the discharge amount of pollutants of rural direct-discharge living sources, Nz-out is the discharge amount of pollutants of agricultural planting sources, Nx-out is the discharge amount of pollutants of livestock and poultry breeding sources, Cs-out is the discharge amount of pollutants of direct-discharge urban living sources, Gs-out is the discharge amount of pollutants of industrial point sources, and a is the river coefficient of rural direct-discharge living sources, and is 0.1 for example; b is an agricultural planting source river coefficient, for example, 0.05, c is a livestock and poultry breeding source river coefficient, for example, 0.01, d is a direct-discharge town living source, for example, 1.0, and e is an industrial point source river coefficient, for example, 1.0. The river entry coefficients of the pollution sources can be obtained by referring to hydrological data, and the pollutant discharge amount is the discharge amount corresponding to the optimization index.

The above-mentioned method for obtaining the pollutant contribution rate of each upstream river reach to the target river reach may be obtained from a database, that is, the database stores the pollutant contribution rate of each upstream river reach under the hot watershed to each downstream river reach downstream of the database in advance, and when the database is used, the corresponding data is obtained. For example, the pollutant contribution rates of the river reach 1 and the river reach 2 to the river reach 3, the river reach 4 and the river reach 5 are stored; the contribution rates of the river reach 3 and the river reach 6 to the pollutants of the river reach 4 and the river reach 5 respectively; river reach 4 contributes to the pollutant of river reach 5. In this embodiment, before obtaining the pollutant contribution rate of each upstream river reach to the target river reach, the pollutant contribution rate of each upstream river reach to each downstream river reach located downstream of the upstream river reach needs to be determined.

Alternatively, in determining the rate of contribution of pollutants from each upstream stretch to each downstream stretch downstream from itself, the process may be: establishing a transition probability matrix according to the upstream and downstream relations of the river reach corresponding to each control unit; establishing a retention coefficient matrix according to the pollutant inflow amount and the pollutant output amount of the river reach corresponding to each control unit and the upstream-downstream relation of each river reach; establishing a pollutant transfer coefficient matrix according to the transfer probability matrix and the retention coefficient matrix; and determining the pollutant contribution rate of each upstream river reach to each downstream river reach positioned downstream of the upstream river reach based on the pollutant transfer coefficient matrix.

When the transition probability matrix is established according to the upstream and downstream relations of the river reach corresponding to each control unit, each control unit is set as a discrete state variable, and the upstream and downstream river networks corresponding to the control units are represented by using the topological relation, so that the upstream and downstream relations of each control unit in the drainage basin can be represented by using the schematic diagram shown in fig. 2. The upstream and downstream relationship of the control unit can be quantified by the following formula H (i, j).

Then a transition probability matrix can be obtained based on the above formula:

wherein H is a transition probability matrix, and referring to fig. 2, the river reach 1 is directly upstream of the river reach 3, so the third column of the first row in the transition probability matrix is 1; river 2 is directly upstream of river 3, so the third column in the second row is also 1, river 3 is directly upstream of river 4, so the fourth column in the third row is 1, river 4 is directly upstream of river 5, so the fifth column in the third row is 1, river 6 is directly upstream of river 4, so the fourth column in the sixth row is 1.

When the retention coefficient matrix is established according to the pollutant inflow amount and the pollutant output amount of each river reach corresponding to each control unit and the upstream and downstream relationship of each river reach, the process is combined with the graph 2, namely, the pollutant inflow amount and the pollutant output amount of each river reach 1, 2, 3, 4, 5 and 6 are respectively obtained, and after the pollutant inflow amount and the pollutant output amount of each river reach are obtained, the retention coefficient matrix can be established by combining the upstream and downstream relationship of each river reach. Further, after the pollutant inflow amount and the pollutant output amount of each river reach are obtained, the pollutant retention coefficient corresponding to each river reach can be calculated by combining the upstream and downstream relations of each river reach; after calculating the pollutant retention coefficients corresponding to the river reach, establishing a retention coefficient matrix by combining the upstream and downstream relations of the river reach:

wherein r in the retention coefficient matrix₁Indicates the retention coefficient of the pollutant, r, corresponding to the river reach 1₂Denotes the retention coefficient of the contaminant, r, corresponding to the river reach 2_nRepresenting the retention coefficient of the contaminant corresponding to the river reach n, in the case of fig. 2, a 6 x 6 retention coefficient matrix may be obtained, for example

After obtaining a transition probability matrix H and a retention coefficient matrix R, a pollutant transition coefficient matrix can be established according to the transition probability matrix and the retention coefficient matrix,

wherein the content of the first and second substances,

is a pollutant transfer coefficient matrix, H is a transfer probability matrix, R is a retention coefficient matrix, and I is an identity matrix.

Then I-R are:

the resulting pollutant transfer coefficient moment is calculated according to a formula

Comprises the following steps:

obtaining a pollutant transfer coefficient matrix

Based on

The contaminant contribution rate Gij of each upstream stretch to each downstream stretch downstream from itself can be determined, where Gij represents the contaminant contribution rate of the upstream stretch i to the downstream stretch j. For example, a matrix

The coefficient of 0.8 in the first row and the third column in (1) is the pollutant contribution rate G13 of the river reach 1 to the river reach 3, and the matrix

The coefficient of 0.8 in the second row and the third column in (1) is the pollutant contribution rate G23 of the river reach 2 to the river reach 3, and the matrix

The coefficient of 0.75 in the third row and the fourth column in the middle is the pollutant contribution rate G34 of the river reach 3 to the river reach 4, and the matrix

The coefficient of 0.82 of the fourth row and the fifth column in the middle is the pollutant contribution rate G45 of the river reach 4 to the river reach 5, and the matrix is

The coefficient of 0.75 in the sixth row and the fourth column in (1) is the pollutant contribution rate G64 of the river reach 6 to the river reach 4.

As can be seen from the topological relationship shown in fig. 2, the river reach 1 and the river reach 2 are both the straight upstream of the river reach 3, the river reach 3 and the river reach 6 are both the straight upstream of the river reach 4, the river reach 4 is the straight upstream of the river reach 5, the river reach 1 and the river reach 2 span the river reach 3 in the middle of the river reach 4, and are the second-level upstream, the river reach 3 and the river reach 6 span the river reach 4 in the middle of the river reach 5, and are the second-level upstream and the third-level upstream, the river reach 1 and the river reach 2 span the river reach 4 and the river reach 5 in the middle of the river reach 5, and are the third-level upstream. For immediate upstream, based on the contaminant transfer coefficient matrix

The contribution Gij of each upstream stretch to the pollutants in its immediate downstream stretch, e.g. G13, G23, G34, G64, G45, is directly obtained. For the secondary upstream, in order to obtain the pollutant contribution rate Gij of the secondary upstream to the downstream thereof, two iterations are required to obtain the pollutant contribution rate Gij of the secondary upstream to the downstream thereof, such as G14, G24, G35, and G65. Namely:

wherein, the matrix

The coefficient of 0.6 in the first row and the fourth column in the middle is the pollutant contribution rate G14 of the river reach 1 to the river reach 4, and the matrix is

The coefficient of 0.6 in the second row and the fourth column in (1) is the pollutant contribution rate G24 of the river reach 2 to the river reach 4, and the matrix

The coefficient of 0.615 in the third row and the fifth column is the pollutant contribution rate G35 of the river reach 3 to the river reach 5, and the matrix is

The coefficient of 0.615 in the sixth row and the fifth column in (1) is the pollutant contribution rate G65 of the river reach 6 to the river reach 5.

For the third-level upstream, three iterations are required to obtain the pollutant contribution rate Gij of the third-level upstream to the downstream thereof, so that the pollutant contribution rate Gij of the third-level upstream to the downstream thereof can be obtained, such as G15 and G25. Namely:

wherein, the matrix

The coefficient of 0.492 of the first row and the fifth column in the matrix is the pollutant contribution rate G15 of the river reach 1 to the river reach 5

The coefficient of 0.492 in the second row and the fifth column in (1) is the pollutant contribution rate G25 of the river reach 2 to the river reach 5.

Obtaining a pollutant transfer coefficient matrix

And then, iterating for different times to finally obtain the pollutant contribution rate Gij of each upstream river reach to each downstream river reach located at the downstream of the upstream river reach. It should be noted that, the above iteration number is only for the topological relation shown in fig. 2, and different topological relations need the same iteration number, for example, if it is desired to calculate the pollutant contribution rate Gij of the upstream of the fourth level to the downstream thereof, it is necessary to iterate four times, for example, if it is assumed that the river reach 7 located downstream of the river reach 5 is also included on the basis of fig. 2, it is necessary to iterate four times when calculating the pollutant contribution rate G17 of the upstream river reach 1 to the river reach 7 of the fourth level and calculating the pollutant contribution rate G27 of the upstream river reach 2 to the river reach 7.

The pollutant contribution Gij of each upstream stretch to each downstream stretch located downstream from itself can be obtained in the manner described above.

Third embodiment

This embodiment is different from the first embodiment in that: as shown in fig. 4, after step S104, the method further includes: and determining the water environment capacity of the target pollution source corresponding to each control unit according to the water environment capacity allowed by each control unit. That is, after the water environment capacity allowed by each control unit is allocated, the water environment capacity can be further allocated secondarily to the target pollution source in each control unit. Of course, the secondary distribution of the water environment capacity to the target pollution source in each control unit may be performed after the correction of the water environment capacity allowed by each distributed control unit is completed, that is, the secondary distribution of the water environment capacity to the target pollution source in the control unit is performed based on the corrected water environment capacity allowed by each control unit, as shown in fig. 5.

When the target pollution source in each control unit is subjected to secondary distribution of water environment capacity, the secondary distribution may be: acquiring the pollutant river entering amount of each preset pollution source in each control unit; and determining the water environment capacity of the target pollution source corresponding to each control unit according to the water environment capacity allowed by each control unit and the acquired river inflow amount of the pollutants of each preset pollution source corresponding to each control unit. For example, when the water environment capacity of a target pollution source in the control unit 1 corresponding to the river reach 1 is secondarily distributed, the river inflow amount Ns-i of rural directly-discharged living source pollutants in the control unit 1, the river inflow amount Nz-i of agricultural planting source pollutants, the river inflow amount Nx-i of livestock and poultry breeding source pollutants, the river inflow amount Gs-i of industrial point source pollutants and the river inflow amount Cs-i of directly-discharged urban living source pollutants are obtained; then based on the water environment capacity W allowed by the control unit 1₀(for the control unit 1, W₀T × W1, the value range of the coefficient t is 0.9-1), and the obtained Ns-i, Nz-i, Nx-i, Gs-i, and Cs-i can determine the water environment capacity of the corresponding target pollution source in the control unit 1, for example, if the target pollution source is a rural directly-discharged living source, the water environment capacity W of the rural directly-discharged living source is determined_s＝W₀(Ns-i/ei) and similarly, when the target pollution source is an agricultural planting source, the water environment capacity W of the agricultural planting source_s＝W₀(Nz-i/ei), similarly, when the target pollution source is the direct-discharge urban living source, the water environment capacity W of the direct-discharge urban living source_s＝W₀*(Cs-i/ei)。

When the control unit 3 secondarily allocates the water environment capacity of the target pollution source in the control unit 3, if the water environment capacity is based on the water environment capacity W allowed by the control unit 3 before the correction₀W3, for example, when assigning the water environmental capacity of a target pollution source, such as an agricultural planting source, within the control unit 3_sT W3 Nz-i/ei); if the water environment capacity W is allowed by the corrected control unit 3₀When the distribution is performed at W1G 13+ W2G 23, and the water environment capacity of the agricultural planting source in the control unit 3 is distributed, W is used_s(W1 × G13+ W2 × G23) × (Nz-i/ei). For another example, when the control unit 4 allocates the water environmental capacity of the target pollution source such as the livestock and poultry raising source in the control unit 4, in one embodiment, W may be_sIn another embodiment, W4 (Nx-i/ei) may be W_s＝(W1*G14+W2*G24+W3*G34+W6*G64)*(Nx-i/ei)。

It should be noted that, the above only illustrates the process of distributing the water environmental capacity of a part of target pollution sources, such as livestock and poultry breeding sources, and the principle of the distribution process of different target pollution sources is the same. In addition, when the water environment capacity of the target pollution source corresponding to each control unit is determined according to the water environment capacity allowed by each control unit, the water environment capacity of the target pollution source in the control unit may be allocated based on the water environment capacity allowed by the control unit, in an embodiment, as shown in fig. 5, at this time, the water environment capacity of the target pollution source in the control unit is allocated based on the water environment capacity allowed by the modified control unit.

Wherein ei is Ns-i + Nz-i + Nx-i + Cs-i + Gs-i. It should be noted that, in the above description, the preset pollution source is a rural directly-discharging living source, an agricultural planting source, a livestock and poultry breeding source, an industrial point source, and a directly-discharging urban living source, and the corresponding target pollution source is at least one of the rural directly-discharging living source, the agricultural planting source, the livestock and poultry breeding source, the industrial point source, and the directly-discharging urban living source.

Wherein, when acquiring the pollutant river entering amount of each preset pollution source in each control unit, the pollutant river entering amount may be: acquiring pollutant discharge amount of each preset pollution source in each control unit and river entering coefficients corresponding to each preset pollution source; and obtaining the river entering amount of the pollutants of each preset pollution source in each control unit according to the pollutant discharge amount of each preset pollution source in each control unit and the river entering coefficient corresponding to each preset pollution source. For example, for the control unit 2, assuming that the preset pollution sources are still 5 kinds as described above, the pollutant emission amount Ns-out and the river entering coefficient a of the rural in-line living source are obtained, assuming 0.1, the pollutant emission amount Nz-out and the river entering coefficient b of the agricultural planting source are obtained, assuming 0.05, the pollutant emission amount Nx-out and the river entering coefficient c of the livestock and poultry breeding source are obtained, assuming 0.01, the pollutant emission amount Cs-ou and the river entering coefficient d of the in-line town living source are obtained, assuming 1.0, the pollutant emission amount Gs-out and the river entering coefficient e of the industrial point source are obtained, assuming 1.0. Then the river inflow of pollutants of each preset pollution source in the control unit 2 can be obtained, for example Ns-i is Ns-out 0.1; nz-i ═ Nz-out 0.05; nx-i ═ Nx-out 0.01; Cs-i-Cs-out 1; gs-i ═ Gs-out 1.

Fourth embodiment

The embodiment of the present application further provides a watershed water environment capacity distribution device 100, as shown in fig. 6. The watershed water environment capacity distribution device 100 includes: a first obtaining module 110, a second obtaining module 120, a third obtaining module 130, and a first determining module 140.

The first acquiring module 110 is configured to acquire a watershed to be allocated, where water environment capacity allocation is required.

The second obtaining module 120 is configured to obtain each control unit related to the drainage basin to be distributed, and river reach corresponding to each control unit one to one. Optionally, the second obtaining module 120 is specifically configured to divide the control units in the domain to be allocated based on the home management principle, and determine the river reach corresponding to each control unit.

A third acquisition module 130 for acquiring a maximum aquatic environment capacity allowed to be discharged per river reach.

The first determination module 140 is configured to determine the allowable aquatic environment capacity of the control unit corresponding to each river reach based on the maximum aquatic environment capacity allowed to be discharged by each river reach.

Optionally, in an embodiment, the watershed water environment capacity distribution device 100 further includes: and the correction module is used for correcting the water environment capacity allowed by the control unit corresponding to each river reach. Optionally, the modification module is specifically configured to: determining that the value of the water environment capacity allowed by the control unit corresponding to the upstream river reach is unchanged; acquiring the retention coefficients of pollutants of other river sections except the uppermost stream section; determining that the value of the water environment capacity allowed by a control unit of the river reach with the pollutant retention coefficient not greater than a preset threshold value is not changed; and correcting the water environment capacity allowed by the control unit of the target river reach with the pollutant retention coefficient larger than the preset threshold value. Optionally, the modification module is specifically configured to: acquiring all upstream river reach positioned at the upstream of the target river reach; acquiring the pollutant contribution rate of each upstream river reach to the target river reach; obtaining the sum of the contribution of all the upstream river reach to the water environment capacity of the target river reach based on the maximum water environment capacity of each upstream river reach and the pollutant contribution rate of each upstream river reach to the target river reach; and correcting the water environment capacity allowed by the control unit corresponding to the target river reach as the contribution sum of the water environment capacity.

And the watershed water environment capacity distribution device 100 further includes: and the second determination module is used for determining that the contribution sum of the water environment capacity is smaller than the water environment capacity allowed by the target river reach before the correction module is used for correcting the water environment capacity allowed by the control unit corresponding to the target river reach to be the contribution sum of the water environment capacity.

And the watershed water environment capacity distribution device 100 further includes: and the third determining module is used for determining the contribution sum of the water environmental capacity to be less than the allowable water environmental capacity of the target river reach before the correcting module is used for acquiring the pollutant contribution rate of each upstream river reach to the target river reach. Optionally, the third determining module is specifically configured to: establishing a transition probability matrix according to the upstream and downstream relations of the river reach corresponding to each control unit; establishing a retention coefficient matrix according to the pollutant inflow amount and the pollutant output amount of the river reach corresponding to each control unit and the upstream-downstream relation of each river reach; establishing a pollutant transfer coefficient matrix according to the transfer probability matrix and the retention coefficient matrix; and determining the pollutant contribution rate of each upstream river reach to each downstream river reach positioned downstream of the upstream river reach based on the pollutant transfer coefficient matrix. Optionally, the third determining module is specifically configured to: calculating pollutant retention coefficients corresponding to all the river reach according to the pollutant inflow amount and the pollutant output amount of the river reach corresponding to all the control units and the upstream and downstream relations of all the river reach; and establishing a retention coefficient matrix based on the respective corresponding pollutant retention coefficients of the river reach and the upstream and downstream relations of the river reach.

And the watershed water environment capacity distribution device 100 further includes: and the calculation module is used for calculating the pollutant retention coefficient corresponding to each river reach before the correction module is used for acquiring the pollutant retention coefficients of other river reaches except the most upstream river reach. Optionally, the calculation module is specifically configured to: calculating the pollutant retention coefficient corresponding to each river reach according to the pollutant inflow amount, the pollutant output amount and the position of each river reach, wherein the calculation formula of the pollutant retention coefficient is r_i＝(Load_in-Load_out)/Load_in，r_iIs the corresponding pollutant retention coefficient of river reach i, Load_inThe pollutant Load capacity corresponding to the river reach i is named as Loadin plus Load (i-1) out, subi is the pollutant inflow capacity corresponding to the river reach i, Load (i-1) out is the pollutant output capacity corresponding to the directly upstream i-1 of the river reach i, and Loadout is the pollutant output capacity corresponding to the river reach i.

In one embodiment, the watershed water environment capacity distribution device 100 further comprises: and the fourth determining module is used for determining the water environment capacity of the target pollution source corresponding to each control unit according to the water environment capacity allowed by each control unit after the first determining module is used for determining the water environment capacity allowed by the control unit corresponding to each river reach. Optionally, the fourth determining module is specifically configured to: acquiring the pollutant river entering amount of each preset pollution source in each control unit; and determining the water environment capacity of the target pollution source corresponding to each control unit according to the water environment capacity allowed by each control unit and the acquired river inflow amount of the pollutants of each preset pollution source corresponding to each control unit. Optionally, the fourth determining module is specifically configured to: acquiring pollutant discharge amount of each preset pollution source in each control unit and river entering coefficients corresponding to each preset pollution source; and obtaining the river entering amount of the pollutants of each preset pollution source in each control unit according to the pollutant discharge amount of each preset pollution source in each control unit and the river entering coefficient corresponding to each preset pollution source.

The target pollution source is at least one of a rural direct-discharge living source, an agricultural planting source, a livestock breeding source, an industrial point source and a direct-discharge urban living source.

It should be noted that, in the present specification, the embodiments are all described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments may be referred to each other.

The device 100 for distributing the volume of the watershed water environment according to the embodiment of the present application has the same implementation principle and the same technical effects as those of the foregoing method embodiments, and for the sake of brief description, reference may be made to the corresponding contents in the foregoing method embodiments for the parts of the device embodiments that are not mentioned.

Fifth embodiment

The embodiment of the present application further provides an electronic device 200 capable of executing the watershed water environment capacity allocation method, as shown in fig. 7. The electronic device 200 includes: a transceiver 210, a memory 220, a communication bus 230, and a processor 240.

The elements of the transceiver 210, the memory 220, and the processor 240 are electrically connected to each other directly or indirectly to achieve data transmission or interaction. For example, the components may be electrically coupled to each other via one or more communication buses 230 or signal lines. The transceiver 210 is used for transceiving data. The memory 220 is used for storing a computer program, such as a software functional module shown in fig. 6, that is, the watershed water environment capacity distribution device 100. The watershed water environment capacity distribution device 100 includes at least one software functional module, which can be stored in the memory 220 in the form of software or firmware (firmware) or solidified in an Operating System (OS) of the electronic device 200. The processor 240 is configured to execute an executable module stored in the memory 220, such as a software functional module or a computer program included in the watershed water environment capacity distribution device 100. For example, the processor 240 is configured to obtain a watershed to be allocated, where water environmental capacity allocation is required; the distribution system is also used for acquiring each control unit related to the watershed to be distributed and river reach corresponding to each control unit one by one; and also for obtaining the maximum water environmental capacity allowed to be discharged per river reach; and determining the water environment capacity allowed by the control unit corresponding to each river reach based on the maximum water environment capacity allowed to be discharged by each river reach.

The Memory 220 may be, but is not limited to, a Random Access Memory (RAM), a Read Only Memory (ROM), a Programmable Read-Only Memory (PROM), an Erasable Read-Only Memory (EPROM), an electrically Erasable Read-Only Memory (EEPROM), and the like.

The processor 240 may be an integrated circuit chip having signal processing capabilities. The processor may be a general-purpose processor, including a Central Processing Unit (CPU), a Network Processor (NP), and the like; but also Digital Signal Processors (DSPs), Application Specific Integrated Circuits (ASICs), Field Programmable Gate Arrays (FPGAs) or other programmable logic devices, discrete Gate or transistor logic devices, discrete hardware components. The various methods, steps, and logic blocks disclosed in the embodiments of the present application may be implemented or performed. A general purpose processor may be a microprocessor or the processor 240 may be any conventional processor or the like.

The electronic device 200 includes, but is not limited to, a computer, a server, and the like.

It should be noted that, in the present specification, the embodiments are all described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments may be referred to each other.

Sixth embodiment

The present embodiment also provides a non-volatile computer-readable storage medium (hereinafter, referred to as a storage medium), where the storage medium stores a computer program, and when the computer program is executed by a computer, the computer program is executed by the electronic device 200 to perform the steps included in the watershed water environment capacity allocation method provided by the foregoing method embodiment.

In the embodiments provided in the present application, it should be understood that the disclosed apparatus and method can be implemented in other ways. The apparatus embodiments described above are merely illustrative, and for example, the flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of apparatus, methods and computer program products according to various embodiments of the present application. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.

In addition, functional modules in the embodiments of the present application may be integrated together to form an independent part, or each module may exist separately, or two or more modules may be integrated to form an independent part.

The functions, if implemented in the form of software functional modules and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present application or portions thereof that substantially contribute to the prior art may be embodied in the form of a software product stored in a storage medium and including instructions for causing a computer device (which may be a personal computer, a notebook computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: various media capable of storing program codes, such as a usb disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk.

The above description is only for the specific embodiments of the present application, but the scope of the present application is not limited thereto, and any person skilled in the art can easily conceive of the changes or substitutions within the technical scope of the present application, and shall be covered by the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.

Claims (26)

1. A watershed water environment capacity allocation method is characterized by comprising the following steps:

acquiring a watershed to be distributed, which needs to perform water environment capacity distribution;

acquiring each control unit related to the watershed to be distributed and river reach corresponding to each control unit one by one;

acquiring the maximum water environment capacity allowed to be discharged by each river reach;

determining the water environment capacity allowed by the control unit corresponding to each river reach based on the maximum water environment capacity allowed to be discharged by each river reach;

correcting the water environment capacity allowed by the control unit corresponding to each river reach, wherein the correction of the water environment capacity allowed by the control unit corresponding to each river reach comprises the following steps:

determining that the value of the water environment capacity allowed by the control unit corresponding to the upstream river reach is unchanged;

acquiring the retention coefficients of pollutants of other river sections except the uppermost stream section;

determining that the value of the water environment capacity allowed by a control unit of the river reach with the pollutant retention coefficient not greater than a preset threshold value is not changed;

correcting the water environment capacity allowed by the control unit of the target river reach with the pollutant retention coefficient larger than the preset threshold, wherein the correction of the water environment capacity allowed by the control unit of the target river reach with the pollutant retention coefficient larger than the preset threshold comprises the following steps:

acquiring all upstream river reach positioned at the upstream of the target river reach;

acquiring the pollutant contribution rate of each upstream river reach to the target river reach;

obtaining the sum of the contribution of all the upstream river reach to the water environment capacity of the target river reach based on the maximum water environment capacity of each upstream river reach and the pollutant contribution rate of each upstream river reach to the target river reach;

and correcting the water environment capacity allowed by the control unit corresponding to the target river reach as the contribution sum of the water environment capacity.

2. The method of claim 1, wherein before modifying the sum of the contribution of the aquatic environment capacity to the aquatic environment capacity allowed by the control unit corresponding to the target river reach, the method further comprises:

determining the contribution of the aquatic environment capacity is less than the aquatic environment capacity allowed by the target river reach.

3. The method of claim 1, wherein prior to obtaining the contaminant contribution rate of each upstream leg to the target leg, the method further comprises:

the rate of contribution of pollutants from each upstream stretch to each downstream stretch located downstream from itself is determined.

4. The method of claim 3, wherein determining the contaminant contribution rate of each upstream stretch to each downstream stretch located downstream from the upstream stretch comprises:

establishing a transition probability matrix according to the upstream and downstream relations of the river reach corresponding to each control unit;

establishing a retention coefficient matrix according to the pollutant inflow amount and the pollutant output amount of the river reach corresponding to each control unit and the upstream-downstream relation of each river reach;

establishing a pollutant transfer coefficient matrix according to the transfer probability matrix and the retention coefficient matrix;

and determining the pollutant contribution rate of each upstream river reach to each downstream river reach positioned downstream of the upstream river reach based on the pollutant transfer coefficient matrix.

5. The method of claim 4, wherein establishing a retention coefficient matrix according to the pollutant inflow amount and the pollutant output amount of the river reach corresponding to each control unit and the upstream and downstream relationship of each river reach comprises:

calculating pollutant retention coefficients corresponding to all the river reach according to the pollutant inflow amount and the pollutant output amount of the river reach corresponding to all the control units and the upstream and downstream relations of all the river reach;

and establishing a retention coefficient matrix based on the respective corresponding pollutant retention coefficients of the river reach and the upstream and downstream relations of the river reach.

6. The method of claim 1, wherein prior to obtaining the contaminant retention factor for the sections other than the most upstream section, the method further comprises: and calculating the pollutant retention coefficient corresponding to each river reach.

7. The method of claim 6, wherein calculating the contaminant retention coefficient for each of the river reach segments comprises:

calculating the pollutant retention coefficient corresponding to each river reach according to the pollutant inflow amount, the pollutant output amount and the position of each river reach, wherein the calculation formula of the pollutant retention coefficient is r_i＝(Load_in-Load_out)/Load_in，r_iIs the corresponding pollutant retention coefficient of river reach i, Load_inThe pollutant Load capacity corresponding to the river reach i is named as Loadin plus Load (i-1) out, subi is the pollutant inflow capacity corresponding to the river reach i, Load (i-1) out is the pollutant output capacity corresponding to the directly upstream i-1 of the river reach i, and Loadout is the pollutant output capacity corresponding to the river reach i.

8. The method according to claim 1, wherein acquiring each control unit related to the watershed to be distributed and river reach corresponding to each control unit in a one-to-one manner comprises:

and dividing the control units in the to-be-distributed flow domain based on a home management principle, and determining the river reach corresponding to each control unit.

9. The method according to any one of claims 1-8, wherein after determining the water environmental capacity allowed by the control unit corresponding to each river reach, the method further comprises:

and determining the water environment capacity of the target pollution source corresponding to each control unit according to the water environment capacity allowed by each control unit.

10. The method of claim 9, wherein determining the aquatic environment capacity of the target pollution source corresponding to each control unit according to the aquatic environment capacity allowed by each control unit comprises:

acquiring the pollutant river entering amount of each preset pollution source in each control unit;

and determining the water environment capacity of the target pollution source corresponding to each control unit according to the water environment capacity allowed by each control unit and the acquired river inflow amount of the pollutants of each preset pollution source corresponding to each control unit.

11. The method of claim 10, wherein obtaining the pollutant river entry volume for each predetermined pollution source in each control unit comprises:

acquiring pollutant discharge amount of each preset pollution source in each control unit and river entering coefficients corresponding to each preset pollution source;

and obtaining the river entering amount of the pollutants of each preset pollution source in each control unit according to the pollutant discharge amount of each preset pollution source in each control unit and the river entering coefficient corresponding to each preset pollution source.

12. The method of claim 10, wherein the target pollution source is at least one of a rural in-line living source, an agricultural planting source, a livestock breeding source, an industrial point source, and an in-line town living source.

13. A watershed water environmental capacity distribution device, comprising:

the system comprises a first acquisition module, a first distribution module and a second acquisition module, wherein the first acquisition module is used for acquiring a watershed to be distributed, which needs to distribute water environment capacity;

the second acquisition module is used for acquiring each control unit related to the watershed to be distributed and river reach corresponding to each control unit one by one;

the third acquisition module is used for acquiring the maximum water environment capacity allowed to be discharged by each river reach;

the first determining module is used for determining the water environment capacity allowed by the control unit corresponding to each river reach based on the maximum water environment capacity allowed to be discharged by each river reach;

the fourth acquisition module is used for acquiring the pollutant retention coefficients of other river sections except the uppermost stream river section;

the correction module is used for correcting the water environment capacity allowed by the control unit of the target river reach with the pollutant retention coefficient larger than a preset threshold, wherein the correction module is specifically used for: acquiring all upstream river reach positioned at the upstream of the target river reach; acquiring the pollutant contribution rate of each upstream river reach to the target river reach; obtaining the sum of the contribution of all the upstream river reach to the water environment capacity of the target river reach based on the maximum water environment capacity of each upstream river reach and the pollutant contribution rate of each upstream river reach to the target river reach; and correcting the water environment capacity allowed by the control unit corresponding to the target river reach as the contribution sum of the water environment capacity.

14. The apparatus of claim 13, further comprising: and the second determination module is used for determining that the contribution sum of the water environment capacity is smaller than the water environment capacity allowed by the target river reach before the correction module is used for correcting the water environment capacity allowed by the control unit corresponding to the target river reach to be the contribution sum of the water environment capacity.

15. The apparatus of claim 13, further comprising: and the third determining module is used for determining the contribution sum of the water environmental capacity to be less than the allowable water environmental capacity of the target river reach before the correcting module is used for acquiring the pollutant contribution rate of each upstream river reach to the target river reach.

16. The apparatus of claim 15, wherein the third determining module is specifically configured to:

establishing a transition probability matrix according to the upstream and downstream relations of the river reach corresponding to each control unit; establishing a retention coefficient matrix according to the pollutant inflow amount and the pollutant output amount of the river reach corresponding to each control unit and the upstream-downstream relation of each river reach; establishing a pollutant transfer coefficient matrix according to the transfer probability matrix and the retention coefficient matrix; and determining the pollutant contribution rate of each upstream river reach to each downstream river reach positioned downstream of the upstream river reach based on the pollutant transfer coefficient matrix.

17. The apparatus of claim 16, wherein the third determining module is specifically configured to: calculating pollutant retention coefficients corresponding to all the river reach according to the pollutant inflow amount and the pollutant output amount of the river reach corresponding to all the control units and the upstream and downstream relations of all the river reach; and establishing a retention coefficient matrix based on the respective corresponding pollutant retention coefficients of the river reach and the upstream and downstream relations of the river reach.

18. The apparatus of claim 13, further comprising: and the calculation module is used for calculating the pollutant retention coefficient corresponding to each river reach before the fourth acquisition module is used for acquiring the pollutant retention coefficients of other river reaches except the most upstream river reach.

19. The apparatus of claim 18, wherein the computing module is specifically configured to: calculating the pollutant retention coefficient corresponding to each river reach according to the pollutant inflow amount, the pollutant output amount and the position of each river reach, wherein the calculation formula of the pollutant retention coefficient is r_i＝(Load_in-Load_out)/Load_in，r_iIs the corresponding pollutant retention coefficient of river reach i, Load_inThe pollutant Load capacity corresponding to the river reach i is named as Loadin plus Load (i-1) out, subi is the pollutant inflow capacity corresponding to the river reach i, Load (i-1) out is the pollutant output capacity corresponding to the directly upstream i-1 of the river reach i, and Loadout is the pollutant output capacity corresponding to the river reach i.

20. The apparatus of claim 13, wherein the second obtaining module is specifically configured to: and dividing the control units in the to-be-distributed flow domain based on a home management principle, and determining the river reach corresponding to each control unit.

21. The apparatus of any one of claims 13-20, further comprising: and the fourth determining module is used for determining the water environment capacity of the target pollution source corresponding to each control unit according to the water environment capacity allowed by each control unit after the first determining module is used for determining the water environment capacity allowed by the control unit corresponding to each river reach.

22. The apparatus of claim 21, wherein the fourth determining module is specifically configured to: acquiring the pollutant river entering amount of each preset pollution source in each control unit; and determining the water environment capacity of the target pollution source corresponding to each control unit according to the water environment capacity allowed by each control unit and the acquired river inflow amount of the pollutants of each preset pollution source corresponding to each control unit.

23. The apparatus of claim 22, wherein the fourth determining module is specifically configured to: acquiring pollutant discharge amount of each preset pollution source in each control unit and river entering coefficients corresponding to each preset pollution source; and obtaining the river entering amount of the pollutants of each preset pollution source in each control unit according to the pollutant discharge amount of each preset pollution source in each control unit and the river entering coefficient corresponding to each preset pollution source.

24. The apparatus of claim 22, wherein the target pollution source is at least one of a rural in-line living source, an agricultural planting source, a livestock breeding source, an industrial point source, and an in-line town living source.

25. An electronic device, comprising: a memory and a processor, the memory and the processor connected;

the memory is used for storing programs;

the processor is to invoke a program stored in the memory to perform the method of any of claims 1-12.

26. A storage medium having stored thereon a computer program which, when executed by a computer, performs the method of any one of claims 1-12.

Images