CN110335172A

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

Info

Publication number: CN110335172A
Application number: CN201910623640.0A
Authority: CN
Inventors: 张维; 郭健; 刘晓东; 李燃; 闫佩; 罗彦鹤; 高郁杰; 唐丽丽; 李怀明; 张雷波; 王兴; 郭洪鹏
Original assignee: Tianjin Huanke Environmental Planning Technology Development Co ltd
Current assignee: Tianjin Huanke Environmental Planning Technology Development Co ltd
Priority date: 2019-07-11
Filing date: 2019-07-11
Publication date: 2019-10-15
Anticipated expiration: 2039-07-11
Also published as: CN110335172B

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

Basin water environment capacity allocation method, device, electronic equipment and storage medium

Technical field

The application belongs to monitoring water environment technical field, and in particular to a kind of basin water environment capacity allocation method, device, Electronic equipment and storage medium.

Background technique

China implements the thinking of water pollutant overall control at present, and environmental capacity of water (refers in particular to meeting quality of water environment Under it is required that, water body hold pollutant maximal workload, therefore also referred to as water body load or pollutant carrying capacity) reasonable distribution be water The key of contaminant transport model.Pollutant discharge amount reduction distributes to quantitative discharge limit not essentially according to certain method Then same control unit carries out pollutant discharge amount reduction according to the blowdown status in control unit.

Currently, the distribution method of environmental capacity of water mainly has proportional allocations method, equitable interval distribution method, linear programming technique Deng.Wherein, proportional allocations method is easy to operate, and indicator of distribution is single, due to the considerations of lacking to regional development, easily causes " more Arrange more points " unjust phenomenon.And equitable interval distribution method is the optimization of reciprocity proportionate allocation unjustness, but still without departing from The defect of proportional allocations single assigned index.Linear programming technique often from economic benefit, embodies different distribution objects Feature difference, but easily cause because of the objective optimization for emphasizing economic benefit the unjust phenomenon of allocation result.

Summary of the invention

In consideration of it, the application be designed to provide a kind of basin water environment capacity allocation method, device, electronic equipment and Storage medium easily causes the unjust phenomenon of allocation result to improve above-mentioned distribution method, and then ensure that the visitor of allocation result The property seen and fairness.

Embodiments herein is achieved in that

In a first aspect, the embodiment of the present application provides a kind of basin water environment capacity allocation method, including obtain need into The basin to be allocated of row environmental capacity of water distribution；Obtain each control unit relevant to the basin to be allocated, and with it is each A one-to-one section of control unit；Obtain the maximum environmental capacity of water that each section allows to discharge；Permitted based on each section Perhaps the maximum environmental capacity of water discharged determines the environmental capacity of water that the corresponding control unit in each section allows；To each section The environmental capacity of water that corresponding control unit allows is modified.In the embodiment of the present application, by the control in basin to be allocated Unit processed carry out divide and clearly each control unit charge for remittance range, then obtain each with the one-to-one section of control unit Allow the maximum environmental capacity of water discharged, so determine again environmental capacity of water that the corresponding control unit in each section allows and Also determining environmental capacity of water is further corrected, is no longer proportional allocations, ensure that allocation result objectivity and Fairness.

A kind of possible embodiment of embodiment with reference to first aspect allows the corresponding control unit in each section Environmental capacity of water is modified, comprising: determines that the value for the environmental capacity of water that the corresponding control unit in most upstream section allows is constant； The pollutant for obtaining other sections in addition to the section of most upstream is detained coefficient；Determine that pollutant is detained coefficient no more than preset threshold Section control unit allow environmental capacity of water value it is constant；The mesh that coefficient is greater than the preset threshold is detained to pollutant The environmental capacity of water that the control unit of mark section allows is modified.In the embodiment of the present application, in amendment, by upper reach Influence of the pollutant to lower reache is taken into account, and it is big to be detained coefficient to the pollutant in other sections in addition to the section of most upstream The environmental capacity of water that control unit in the target section of preset threshold allows is modified, by the pollutant of upper reach under The influence of trip section is taken into account, to be further ensured that the accuracy and objectivity of allocation result.

A kind of possible embodiment of embodiment with reference to first aspect is detained coefficient to pollutant and is greater than the default threshold The environmental capacity of water that the control unit of the target section of value allows is modified, comprising: obtains the institute for being located at target section upstream There is upper reach；Each upper reach is obtained to the pollutant contribution rate of the target section；Maximum based on each upper reach Environmental capacity of water and each upper reach obtain all upper reaches pair to the pollutant contribution rate of the target section The target section environmental capacity of water contribution and；The environmental capacity of water that correcting the corresponding control unit in the target section allows is institute State environmental capacity of water contribution and.In the embodiment of the present application, in the environmental capacity of water allowed the corresponding control unit in target section When being modified, by obtaining pollutant contribution rate of all upper reaches for being located at target section upstream to the target section, The maximum environmental capacity of water and each upper reach for being then based on each upper reach contribute the pollutant of the target section Rate, obtain all upper reaches to the environmental capacity of water contribution of the target section and, finally by the corresponding control in target section The environmental capacity of water that unit processed allows is modified to the environmental capacity of water and contributes and consider pollutant in distribution and migrating Conversion in journey ensure that the science of the accounting result of contribution amount.

A kind of possible embodiment of embodiment with reference to first aspect is correcting the corresponding control unit in target section The environmental capacity of water of permission be the environmental capacity of water contribute and before, the method also includes: determine the environmental capacity of water Contribution and the environmental capacity of water allowed less than the target section.In the embodiment of the present application, in the corresponding control in amendment target section The environmental capacity of water that unit processed allows is that the environmental capacity of water is contributed and before, first determines environmental capacity of water contribution and be small In the environmental capacity of water that the target section allows, to guarantee the accuracy of correction result.

A kind of possible embodiment of embodiment with reference to first aspect is obtaining each upper reach to the target section Pollutant contribution rate before, the method also includes: determine each upper reach to be located at itself downstream each downstream The pollutant contribution rate of section.In the embodiment of the present application, by the way that each upper reach is determined in advance to positioned at the every of itself downstream The pollutant contribution rate of one lower reache, to directly acquire corresponding data when there is demand, with save the cost and Accelerate treatment effeciency.

A kind of possible embodiment of embodiment with reference to first aspect determines each upper reach to positioned at itself downstream Each lower reache pollutant contribution rate, comprising: according to the upstream and downstream of the corresponding section of each control unit close System, establishes transition probability matrix；Enter river amount, pollutant output quantity according to the pollutant of the corresponding section of each control unit And the upstream-downstream relationship of each section, it establishes and is detained coefficient matrix；According to the transition probability matrix, the delay coefficient square Battle array establishes pollutant transfer-coefficient matrix；Based on the pollutant transfer-coefficient matrix, determine each upper reach to positioned at certainly The pollutant contribution rate of each lower reache in body downstream.It, will by introducing Markov Chain method in the embodiment of the present application Each control unit is set as discrete state variable, and carries out table with topological relation to the corresponding upstream and downstream network of waterways of control unit Sign, specifies the delay system of the migration transport path and clear pollutant of river pollutant in survey region in each section Number is established and is detained coefficient matrix, is then based on transition probability matrix, is detained the pollutant transfer ratio square that coefficient matrix is established Battle array determines pollutant contribution rate of each upper reach to each lower reache for being located at itself downstream, ensure that contribution amount Accounting result science and accuracy.

A kind of possible embodiment of embodiment with reference to first aspect, according to the corresponding section of each control unit Pollutant enter the upstream-downstream relationship of river amount, pollutant output quantity and each section, establish and be detained coefficient matrix, comprising: root Enter river amount, pollutant output quantity and the upstream and downstream of each section according to the pollutant of the corresponding section of each control unit to close System calculates the corresponding pollutant in each section and is detained coefficient；It is detained coefficient based on the corresponding pollutant in each section And the upstream-downstream relationship of each section, it establishes and is detained coefficient matrix.

A kind of possible embodiment of embodiment with reference to first aspect is obtaining other sections in addition to the section of most upstream Pollutant be detained before coefficient, the method also includes: calculate the corresponding pollutant in each section and be detained coefficient.This Shen Coefficient please be detained by the corresponding pollutant in each section of calculated in advance, directly to obtain when there is demand in embodiment Corresponding data are taken, with save the cost and accelerate treatment effeciency.

A kind of possible embodiment of embodiment with reference to first aspect, it is stagnant to calculate the corresponding pollutant in each section Stay coefficient, comprising: river amount, pollutant output quantity and the position locating for itself are entered according to the corresponding pollutant in each section It sets, calculates the corresponding pollutant in each section and be detained coefficient, wherein the calculation formula that pollutant is detained coefficient is r_i= (Load_in-Load_out)/Load_in, r_iIt is detained coefficient, Load for the corresponding pollutant of section i_inFor the corresponding pollutant of section i Load, Loadin=subi+Load (i-1) out, subi are that i corresponding pollutant in section enters river amount, and Load (i-1) out is Section i directly under the corresponding pollutant output quantity of upstream i-1, Loadout is the corresponding pollutant output quantity of section i.The application In embodiment, when calculating pollutant delay coefficient, the influence directly under upstream to lower reache is fully considered, ensure that calculating As a result reasonability.

A kind of possible embodiment of embodiment with reference to first aspect obtains relevant each to the basin to be allocated Control unit, and with the one-to-one section of each control unit, comprising: based on localized management principle to the stream to be allocated Control unit in domain is divided, and determines each corresponding section of control unit.In the embodiment of the present application, it is based on possession Management principle divides the control unit in basin to be allocated, and determines each corresponding section of control unit, operation Simply, convenient for decision of environment management person because gesture applies plan.

A kind of possible embodiment of embodiment with reference to first aspect allows to the corresponding control unit in each section Environmental capacity of water be modified after, the method also includes the water environments that are allowed according to revised each control unit Capacity determines the environmental capacity of water in the corresponding target stains source of each control unit.In the embodiment of the present application, after to amendment Each control unit in the environmental capacity of water in target stains source further distributed, it is dirty with the specific water for controlling all kinds of pollution sources Contaminate object discharge amount.

A kind of possible embodiment of embodiment with reference to first aspect allows according to revised each control unit Environmental capacity of water determines the environmental capacity of water in the corresponding target stains source of each control unit, comprising: obtain each control unit The pollutant of interior each default pollution sources enters river amount；According to the environmental capacity of water of revised each control unit permission and obtain The pollutant of the corresponding each default pollution sources of each control unit got enters river amount, determines the corresponding mesh of each control unit Mark the environmental capacity of water of pollution sources.In the embodiment of the present application, the environmental capacity of water in target stains source in some control unit When further being distributed, the pollutant by obtaining each default pollution sources in the control unit enters river amount, then basis The environmental capacity of water that the pollutant of each default pollution sources enters river amount and the revised control unit and allows determines target dirt The environmental capacity of water in dye source, it has fully considered the contribution of each sewage unit, ensure that the reasonability of allocation result.

A kind of possible embodiment of embodiment with reference to first aspect obtains each default dirt in each control unit The pollutant in dye source enter river amount, comprising: obtain each default pollution sources in each control unit pollutant discharge amount and Each default corresponding river pollutant sources of pollution sources；According to the pollutant discharge amount of each default pollution sources in each control unit And each the default corresponding river pollutant sources of pollution sources, the pollutant for obtaining each default pollution sources in the control unit enter River amount.In the embodiment of the present application, by obtaining the pollutant discharge amount of each default pollution sources in control unit and each The default corresponding river pollutant sources of pollution sources are entered river amount with this to obtain the pollutant of corresponding default pollution sources, ensure that with this Pollutant enters the controllability and reasonability of river amount.

A kind of possible embodiment of embodiment with reference to first aspect, the target stains source are the in line life in rural area At least one of source, agricultural planting source, livestock and poultry cultivation source, industrial point source, in line urban life source.In the embodiment of the present application, When the environmental capacity of water in target stains source is further distributed in control unit, it is contemplated that the influence biggish master of result Pollution sources are flowed, and then ensure that the reasonability of allocation result.

Second aspect, the embodiment of the present application also provides a kind of basin water environment capacity allocation devices, comprising: first obtains Module, second obtain module, third acquisition module, the first determining module and correction module；First obtains module, for obtaining Need to carry out the basin to be allocated of environmental capacity of water distribution；Second obtains module, related to the basin to be allocated for obtaining Each control unit, and with the one-to-one section of each control unit；Third obtains module, for obtaining each section Allow the maximum environmental capacity of water discharged；First determining module, the maximum water environment for allowing to discharge based on each section are held Amount determines the environmental capacity of water that the corresponding control unit in each section allows；Correction module, for the corresponding control in each section The environmental capacity of water that unit processed allows is modified.

In conjunction with a kind of possible embodiment of second aspect embodiment, the correction module is specifically used for: determination is most upper The value for swimming the environmental capacity of water that the corresponding control unit in section allows is constant；Obtain the dirt of other sections in addition to the section of most upstream It contaminates object and is detained coefficient；Determine that pollutant is detained coefficient no more than the environmental capacity of water that the control unit of the section of preset threshold allows Value it is constant；The environmental capacity of water that the control unit that pollutant is detained the target section that coefficient is greater than the preset threshold is allowed It is modified.

In conjunction with a kind of possible embodiment of second aspect embodiment, the correction module is specifically used for: acquisition is located at All upper reaches of target section upstream；Each upper reach is obtained to the pollutant contribution rate of the target section；Based on every The maximum environmental capacity of water and each upper reach of a upper reach obtain institute to the pollutant contribution rate of the target section Have the upper reach to the contribution of the environmental capacity of water of the target section and；Correcting the corresponding control unit in the target section allows Environmental capacity of water be the environmental capacity of water contribution and.

In conjunction with a kind of possible embodiment of second aspect embodiment, further includes: the second determining module, for described Correction module is used to correct the environmental capacity of water that the corresponding control unit in the target section allows as environmental capacity of water contribution Before, environmental capacity of water environmental capacity of water contribution and allowed less than the target section is determined.

In conjunction with a kind of possible embodiment of second aspect embodiment, further includes: third determining module, for described Correction module is for determining the environmental capacity of water before obtaining each upper reach to the pollutant contribution rate of the target section Contribution and the environmental capacity of water allowed less than the target section.

In conjunction with a kind of possible embodiment of second aspect embodiment, the third determining module is specifically used for: according to The upstream-downstream relationship of each corresponding section of control unit, establishes transition probability matrix；Respectively according to each control unit The pollutant of corresponding section enters river amount, pollutant output quantity and the upstream-downstream relationship of each section, establishes and is detained coefficient square Battle array；Pollutant transfer-coefficient matrix is established according to the transition probability matrix, the delay coefficient matrix；Based on the pollutant Transfer-coefficient matrix determines pollutant contribution rate of each upper reach to each lower reache for being located at itself downstream.

In conjunction with a kind of possible embodiment of second aspect embodiment, the third determining module is specifically used for: according to The pollutant of each corresponding section of control unit enters river amount, pollutant output quantity and the upstream and downstream of each section and closes System calculates the corresponding pollutant in each section and is detained coefficient；It is detained coefficient based on the corresponding pollutant in each section And the upstream-downstream relationship of each section, it establishes and is detained coefficient matrix.

In conjunction with a kind of possible embodiment of second aspect embodiment, further includes: computing module, in the amendment The pollutant that module is used to obtain other sections in addition to the section of most upstream is detained before coefficient, calculates each section and respectively corresponds to Pollutant be detained coefficient.

In conjunction with a kind of possible embodiment of second aspect embodiment, the computing module is specifically used for: according to each The corresponding pollutant in section enters river amount, pollutant output quantity and itself location, and it is respectively right to calculate each section The pollutant answered is detained coefficient, wherein the calculation formula that pollutant is detained coefficient is r_i=(Load_in-Load_out)/Load_in, r_i It is detained coefficient, Load for the corresponding pollutant of section i_inFor the corresponding pollutant burden amount of section i, Loadin=subi+Load (i-1) out, subi are that i corresponding pollutant in section enters river amount, and Load (i-1) out is the corresponding directly under upstream i-1 of section i Pollutant output quantity, Loadout are the corresponding pollutant output quantity of section i.

In conjunction with a kind of possible embodiment of second aspect embodiment, described second obtains module, is specifically used for: being based on Localized management principle divides the control unit in the basin to be allocated, and determines each corresponding river of control unit Section.

In conjunction with a kind of possible embodiment of second aspect embodiment, further includes: the 4th determining module, for described After amendment determining module is modified the environmental capacity of water that the corresponding control unit in each section allows, according to revised The environmental capacity of water that each control unit allows, determines the environmental capacity of water in the corresponding target stains source of each control unit.

In conjunction with a kind of possible embodiment of second aspect embodiment, the 4th determining module is specifically used for: obtaining The pollutant of each default pollution sources in each control unit enters river amount；The environmental capacity of water allowed according to each control unit The pollutant of each default pollution sources corresponding with each control unit got enters river amount, determines that each control unit is corresponding Target stains source environmental capacity of water.

In conjunction with a kind of possible embodiment of second aspect embodiment, the 4th determining module is specifically used for: obtaining The pollutant discharge amount of each default pollution sources in each control unit and each default corresponding river pollutant sources of pollution sources； According to the pollutant discharge amount of each default pollution sources in each control unit and each default pollution sources it is corresponding enter river Coefficient, the pollutant for obtaining each default pollution sources in the control unit enter river amount.

In conjunction with a kind of possible embodiment of second aspect embodiment, the target stains source is the in line life in rural area At least one of source, agricultural planting source, livestock and poultry cultivation source, industrial point source, in line urban life source.

The third aspect, the embodiment of the present application also provides a kind of electronic equipment, comprising: memory and processor, it is described to deposit Reservoir is connected with the processor；The memory is for storing program；The processor is stored in the storage for calling Program in device is to execute above-mentioned first aspect embodiment and/or with reference to first aspect any possible embodiment party of embodiment The method that formula provides.

Fourth aspect, the embodiment of the present application also provides a kind of storage mediums, are stored thereon with computer program, the meter Calculation machine program executes above-mentioned first aspect embodiment when being run by computer and/or any of embodiment can with reference to first aspect The method that the embodiment of energy provides.

Other feature and advantage of the application will be illustrated in subsequent specification, also, partly be become from specification It is clear that being understood and implementing the embodiment of the present application.The purpose of the application and other advantages can be by written Specifically noted structure is achieved and obtained in specification and attached drawing.

Detailed description of the invention

In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, below will be to institute in embodiment Attached drawing to be used is needed to be briefly described, it should be apparent that, the accompanying drawings in the following description is only some implementations of the application Example, for those of ordinary skill in the art, without creative efforts, can also obtain according to these attached drawings Obtain other attached drawings.By the way that shown in attached drawing, above and other purpose, the feature and advantage of the application will be more clear.In whole Identical appended drawing reference indicates identical part in attached drawing.Attached drawing, emphasis deliberately are not drawn by actual size equal proportion scaling It is that the purport of the application is shown.

Fig. 1 shows a kind of flow diagram of basin water environment capacity allocation method provided by the embodiments of the present application.

Fig. 2 shows the charge for remittance space topological of each control unit in basin to be allocated provided by the embodiments of the present application passes The schematic diagram of system.

Fig. 3 shows the flow diagram of another basin water environment capacity allocation method provided by the embodiments of the present application.

Fig. 4 shows the flow diagram of another basin water environment capacity allocation method provided by the embodiments of the present application.

Fig. 5 shows the flow diagram of another basin water environment capacity allocation method provided by the embodiments of the present application.

Fig. 6 shows a kind of module diagram of basin water environment capacity allocation device provided by the embodiments of the present application.

Fig. 7 shows the structural schematic diagram of a kind of electronic equipment provided by the embodiments of the present application.

Specific embodiment

Below in conjunction with the attached drawing in the embodiment of the present application, technical solutions in the embodiments of the present application is described.

It should also be noted that similar label and letter indicate similar terms in following attached drawing, therefore, once a certain Xiang Yi It is defined in a attached drawing, does not then need that it is further defined and explained in subsequent attached drawing.Meanwhile the application's The relational terms of such as " first ", " second " or the like are used merely to an entity or operation and another entity in description Or operation distinguishes, without necessarily requiring or implying between these entities or operation there are any this actual relationship or Person's sequence.Moreover, the terms "include", "comprise" or any other variant thereof is intended to cover non-exclusive inclusion, to make Obtaining the process, method, article or equipment including a series of elements not only includes those elements, but also including not arranging clearly Other element out, or further include for elements inherent to such a process, method, article, or device.Not more In the case where limitation, the element that is limited by sentence "including a ...", it is not excluded that including process, the side of the element There is also other identical elements in method, article or equipment.

Furthermore term "and/or" in the application, only a kind of incidence relation for describing affiliated partner, expression can deposit In three kinds of relationships, for example, A and/or B, can indicate: individualism A exists simultaneously A and B, these three situations of individualism B.

First embodiment

Referring to Fig. 1, being a kind of basin water environment capacity allocation method provided by the embodiments of the present application, below in conjunction with figure 1 pair of its step of including, is illustrated.

Step S101: the basin to be allocated for needing to carry out environmental capacity of water distribution is obtained.

Step S102: obtaining each control unit relevant to the basin to be allocated, and with each control unit one One corresponding section.

After getting basin to be allocated, obtain relevant to basin to be allocated each control unit and with each control The one-to-one section of unit.For example, the control unit in the basin to be allocated is divided based on localized management principle, And determine each corresponding section of control unit, and then available each control relevant with the basin to be allocated is singly Member, and with the one-to-one section of each control unit.For example, the town Yi Jie is each Subject of Administration according to localized management principle, Control unit in the basin to be allocated is divided, Administrative boundaries is comprehensively considered and river charge for remittance range determines each control The corresponding section of unit.

As an implementation, the control unit of acquisition and corresponding section can be obtains from database, It that is to say and be previously stored in the database, obtain corresponding data from database when having demand.For example, right in advance Control unit and corresponding section in popular basin carry out division confirmation, and in the database by the storage of obtained data, Namely multiple popular corresponding control units in basin and corresponding section data are stored in database, for example, storage There are the control unit and corresponding section data of the popular basin A of division, is stored with the control unit of the popular basin B of division And corresponding section data etc..

When being divided to the control unit in basin to be allocated, and determining the corresponding section of each control unit, It needs the data such as the hydrology to basin to be allocated, water quality to be collected and analyze, establishes basic database, be with basic database The network of waterways water in research basin is extracted in support using GIS-Geographic Information System (Geographic Information Systems, GIS) System's figure, and distribute to each control unit, specify the corresponding section of each control unit.

Step S103: the maximum environmental capacity of water that each section allows to discharge is obtained.

After the control unit in basin clearly to be allocated and corresponding section, the maximum that each section allows to discharge is obtained Environmental capacity of water, the maximum environmental capacity of water based on section in existing hydrology-water quality data watershed are calculated.According to river Flow the differential equation of pollutant one-dimensional stable attenuation law:

And byIt solves:, wherein C is along journey pollutant concentration mg/L；u For cross section of river mean flow rate m/s；X is along journey distance km；K is the comprehensive degradation coefficient d-1 of pollutant；C0 is upstream section water Pollutant concentration mg/L.

Sewage draining exit pollutant emission is included in meter by the influence in view of the pollutant emission of Drain contamination for river channel mouth to water quality of river It calculates, and it is required along journey distance far smaller than channel length, and pollutant to assume that the pollutant for being discharged into river is thoroughly mixed It can be uniformly mixed on river channel cross section within a short period of time, then section pollutant concentration under section are as follows:

Enable C_L=C_SThe environmental capacity of water of each section is obtained by above-mentioned formula (1) and formula (2) are as follows:

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

Wherein, C_STo study section control section target water quality mg/L, and definite value can be regarded as, it is expected water up to standard Matter, C_LFor section pollutant concentration mg/L under section；C₀For section pollutant concentration mg/L on section；W_jIndicate that section j allows to arrange The maximum environmental capacity of water t put；Q indicates generalization sewage draining exit quantity of wastewater effluent m³.Wherein, C_SDefinite value, C can be regarded as₀And Q It is different and different according to research object.

Research watershed control section generally may be provided at water quality monitoring mouth, intake, upstream and downstream key node in basin Or the positions such as basin outlet.Each each section, which can be acquired, by data substitution formula (3) for obtaining monitoring allows to discharge Maximum environmental capacity of water.Optionally, data needed for calculating the maximum environmental capacity of water that each section allows can pass through reality Ground measurement obtains.

It wherein, is that will be located at the multiple sewage draining exits closed in the same control unit generally to change when progress sewage draining exit is generally changed The sewage draining exit concentrated for one.The calculation method of the pollutant concentration of generalization sewage draining exit can be, to the same generalization sewage draining exit The pollutant concentration of interior multiple sewage draining exits adds up, and obtains the pollutant concentration of generalization sewage draining exit.For example, by 3 blowdowns Mouth is generalized as a sewage draining exit, then the pollutant concentration of the generalization sewage draining exit is equal to the tired of the pollutant concentration of this 3 sewage draining exits Add, obtains the pollutant concentration of generalization sewage draining exit.

Step S104: allow the maximum environmental capacity of water discharged based on each section, determine the corresponding control in each section The environmental capacity of water that unit allows.

Based on above-mentioned formula (3) after obtaining the maximum environmental capacity of water that each section allows to discharge, it is based on each section Allow the maximum environmental capacity of water discharged, determines the environmental capacity of water that the corresponding control unit in each section allows.

In order to make it easy to understand, being illustrated by taking the control unit in some basin to be allocated shown in Fig. 2 as an example below, it is assumed that It is 6 that the control unit in some basin to be allocated, which delimited, analyzes the charge for remittance range of this 6 control units, obtains Fig. 2 Shown in charge for remittance spatial topotaxy.Wherein, the arrow in figure indicates confluence direction.Each river is calculated based on above-mentioned formula (3) The principle for allowing the maximum environmental capacity of water discharged of section is the same, in order to make it easy to understand, being said below by taking W1 as an example It is bright, by the control section target water quality C of standard_s, and measure the cross section of river mean flow rate u of obtained section 1, along journey distance X, generalization sewage draining exit quantity of wastewater effluent Q, section pollutant concentration C on section₀It substitutes into formula among the above and W1 permission can be obtained The maximum environmental capacity of water of discharge.Wherein, since section 1, section 2, section 6 are in most upstream, section pollutes on section Object concentration C₀The as pollutant concentration of itself, and for the section in downstream 3, section 4, section 5, C₀As upstream river The pollutant concentration in section exit (namely itself inlet).Allow the maximum environmental capacity of water W of discharge obtaining each section_j Afterwards, the environmental capacity of water that the corresponding control unit in each section allows is determined, for example, what the corresponding control unit in section 1 allowed Environmental capacity of water W₀W1 can be equal to multiplied by a coefficient, such as W₀=t*W1, wherein the value range of coefficient t is 0.9-1, such as works as t When taking 1, the environmental capacity of water that control unit allows is equal to the maximum environmental capacity of water for allowing to discharge of corresponding section, such as controls The environmental capacity of water that unit 1 processed allows is equal to W1.

Second embodiment

The present embodiment compared with first embodiment, the difference is that, as shown in figure 3, after step s 104, this method Further include: the environmental capacity of water allowed each control unit is modified.In amendment, determine that the section of most upstream is corresponding The value for the environmental capacity of water that control unit allows is constant, such as section 1 in Fig. 2, section 2, the corresponding control unit in section 6 permit Perhaps the value of environmental capacity of water remains unchanged.The water ring that the corresponding control unit in other sections in addition to the section of most upstream is allowed When border capacity is modified, process may is that the pollutant for obtaining other sections in addition to the section of most upstream is detained coefficient；It is right Pollutant is detained coefficient and is greater than preset threshold (as an implementation, which can be 0, or close to 0 positive number Environmental capacity of water as the control unit of target section 0.1) allows is modified.For Fig. 2, section 3, section are obtained 4, the respective pollutant in section 5 is detained coefficient, is detained the control unit that coefficient is greater than the target section of preset threshold to pollutant The environmental capacity of water of permission is modified, it is assumed that wherein section 3, section 4, section 5 pollutant be detained coefficient be all larger than it is default Threshold value, then section 3, section 4, section 5 are target section；Similarly, if assuming, the pollutant for wherein there was only section 5 is detained coefficient Greater than preset threshold, then only section 5 is target section.Pollutant is detained coefficient and removes most upstream section no more than preset threshold The value for the environmental capacity of water that the control unit of outer section allows is constant.

Wherein, hold in the water environment that the control unit for being detained the target section that coefficient is greater than preset threshold to pollutant allows When amount is modified, it may is that all upper reaches for obtaining and being located at target section upstream；Each upper reach is obtained to the mesh Mark the pollutant contribution rate of section；Maximum environmental capacity of water and each upper reach based on each upper reach are to the mesh Mark the pollutant contribution rate of section, obtain all upper reaches to the environmental capacity of water contribution of the target section and；Amendment The environmental capacity of water that the corresponding control unit in the target section allows be the environmental capacity of water contribute and.In order to make it easy to understand, It is illustrated in conjunction with Fig. 2, by taking target section is section 3 as an example, then the section positioned at 3 upstream of target section is section 1 and section 2, section 1 is obtained to the pollutant contribution rate of section 3, it is assumed that is G13, is obtained section 2 to the pollutant contribution rate of section 3, vacation It is set as G23, then the maximum environmental capacity of water W1 based on the section 1 and pollutant contribution rate G13 to the target section 3, and The maximum environmental capacity of water W2 of section 2 and pollutant contribution rate G23 to the target section 3 is obtained to the target section 3 Environmental capacity of water is contributed and is W1*G13+W2*G23, and the environmental capacity of water of the corresponding control unit in target section 3 permission is corrected And, namely before amendment, the environmental capacity of water of the permission of target section 3 is W3, is W1* after amendment for environmental capacity of water contribution G13+W2*G23。

By taking target section is section 4 as an example, then the section positioned at 4 upstream of target section is section 1, section 2, section 3, river Section 6 obtains section 1 to the pollutant contribution rate of section 4, it is assumed that and it is G14, obtains section 2 to the pollutant contribution rate of section 4, It is assumed to be G24, obtains section 3 to the pollutant contribution rate of section 4, it is assumed that is G34, obtains pollutant tribute of the section 6 to section 4 Offer rate, it is assumed that be G64, the then maximum environmental capacity of water W1 based on section 1 and the pollutant contribution rate to the target section 4 G14, the maximum environmental capacity of water W2 of section 2 and the maximum water ring to the pollutant contribution rate G24 of the target section 4, section 3 Border capacity W3 and to the pollutant contribution rate G34 of the target section 4 and maximum environmental capacity of water W6 of section 6 and right The pollutant contribution rate G64 of the target section 4, acquisition contribute the environmental capacity of water of the target section 4 and are W1*G14+W2* G24+W3*G34+W6*G64 corrects the environmental capacity of water of the corresponding control unit in target section 4 permission as water environment appearance Amount contribution and, namely before amendment, the environmental capacity of water of the permission of target section 4 is W4, is W1*G14+W2*G24+W3* after amendment G34+W6*G64.When target section is section 5, makeover process is similar with target section 3 or target section 4, herein no longer Burden.

Wherein, as an implementation, in the environmental capacity of water for correcting the corresponding control unit permission in the target section For environmental capacity of water contribution and when, can first judge environmental capacity of water that the corresponding control unit in target section allows is whether It is contributed greater than the environmental capacity of water and is determining water environment environmental capacity of water contribution and allowed less than the target section When capacity, just it is modified, for example, first judge whether W1*G13+W2*G23 is greater than W3 when being modified to target section 3, It only when W1*G13+W2*G23 is less than W3, is just modified, the environmental capacity of water that otherwise target section 3 allows remains unchanged It is still W3, the process of remaining target section is similar therewith.

Wherein, the pollutant of other sections of the above-mentioned acquisition in addition to the section of most upstream is detained coefficient, as a kind of implementation Mode, acquisition modes can be each under multiple popular basins from being previously stored in obtain in database namely database The pollutant of a section is detained coefficient, obtains corresponding data when in use.Such as it is stored with each section shown in Fig. 2 Pollutant be detained coefficient.Under this kind of embodiment, it is detained system in the pollutant for obtaining other sections in addition to the section of most upstream Before number, needs to calculate the corresponding pollutant in each section and be detained coefficient.

When in each section of calculating, corresponding pollutant is detained coefficient, it can be corresponding according to each section Pollutant enters river amount, pollutant output quantity and itself location, calculates the corresponding pollutant in each section and is detained Coefficient.Wherein, it is r that pollutant, which is detained the calculation formula of coefficient,_i=(Load_in-Load_out)/Load_in, r_iIt is corresponding for section i Pollutant is detained coefficient, Load_inFor the corresponding pollutant burden amount of section i, Loadin=subi+Load (i-1) out, subi Entering river amount for the corresponding pollutant of section i, Load (i-1) out is section i directly under the corresponding pollutant output quantity of upstream i-1, Loadout is the corresponding pollutant output quantity of section i.For example, for section 1, corresponding pollutant burden amount is equal to pair The pollutant answered enters river amount namely Loadin=subi, and for section 3, pollutant burden amount is equal to the pollution of section 3 Object enters river amount+section 2 pollutant output quantity+section 1 pollutant output quantity namely Loadin=sub3+Load2out+ Load1out。

Wherein, the respective pollutant output quantity in above-mentioned each section can pass through the dirt of the exit end face of the respective section of measurement It contaminates object output quantity and obtains.

Wherein, the above-mentioned respective pollutant in each section enters river amount, can be by setting the dirt in corresponding control unit The pollutant discharge amount in dye source and obtain, wherein pollutant enter river amount calculation formula: ei=Ns-i+Nz-i+Nx-i+Cs-i+ Gs-i, wherein Ns-i is that the in line life source pollutants in rural area enter river amount, and Nz-i is that agricultural planting source pollutants enter river amount, Nx-i Enter river amount for livestock and poultry cultivation source pollutants, Gs-i is that industrial point source pollutant enters river amount, and Cs-i is the pollution of in line urban life source Object enters river amount；And wherein, Ns-i=Ns-out*a；Nz-i=Nz-out*b；Nx-i=Nx-out*c；Cs-i=Cs-out*d； Gs-i=Gs-out*e, wherein Ns-out is the in line life source pollutants discharge amount in rural area, and Nz-out is the pollution of agricultural planting source Object discharge amount, Nx-out are livestock and poultry cultivation source pollutants discharge amount, and Cs-out is in line urban life source discharge amount of pollution, Gs- Out is industrial point source pollutant discharge amount, and above-mentioned a is the in line life source river pollutant sources in rural area, such as takes 0.1；Above-mentioned b is Agricultural planting source river pollutant sources, such as 0.05 is taken, above-mentioned c is livestock and poultry cultivation source river pollutant sources, such as takes 0.01, above-mentioned d For in line urban life source, such as 1.0 are taken, above-mentioned e is industrial point source river pollutant sources, such as takes 1.0.Wherein, above-mentioned each The river pollutant sources of a pollution sources can be obtained refering to hydrological data, and above-mentioned pollutant discharge amount is the corresponding row of optimizing index High-volume.

Wherein, the above-mentioned each upper reach of acquisition can be the pollutant contribution rate of the target section, acquisition modes From each upper reach being previously stored in obtained in database namely database under popular basin to positioned at itself downstream Each lower reache pollutant contribution rate, obtain corresponding data when in use.Such as it is stored with section 1, river Section 2 is respectively to the pollutant contribution rate of section 3, section 4, section 5；Section 3, section 6 are respectively to the pollutant of section 4, section 5 Contribution rate；Pollutant contribution rate of the section 4 to section 5.Under this kind of embodiment, each upper reach is being obtained to the target river Before the pollutant contribution rate of section, it is thus necessary to determine that pollution of each upper reach to each lower reache for being located at itself downstream Object contribution rate.

Optionally, pollutant contribution rate of each upper reach to each lower reache for being located at itself downstream is being determined When, process may is that according to the upstream-downstream relationship of the corresponding section of each control unit, establish transition probability matrix； Enter river amount, pollutant output quantity and the upstream and downstream of each section according to the pollutant of the corresponding section of each control unit Relationship is established and is detained coefficient matrix；Pollutant transfer ratio is established according to the transition probability matrix, the delay coefficient matrix Matrix；Based on the pollutant transfer-coefficient matrix, determine each upper reach to each the downstream river for being located at itself downstream The pollutant contribution rate of section.

It, will be each when establishing transition probability matrix according to the upstream-downstream relationship of the corresponding section of each control unit Control unit is set as discrete state variable, and is characterized to the corresponding upstream and downstream network of waterways of control unit with topological relation, Then each control unit upstream-downstream relationship in basin can be characterized with schematic diagram shown in Fig. 2.Can with following formula H (i, J) upstream-downstream relationship of control unit quantified.

Then it is based on The available transition probability matrix of above-mentioned formula:

Wherein, H is transition probability matrix, and for Fig. 2, section 1 is section 3 directly under upstream, therefore transition probability The third of the first row in matrix is classified as 1；Section 2 is section 3 directly under upstream, therefore the third column of the second row are also 1, section 3 be section 4 directly under upstream, therefore the 4th of the third line is classified as 1, and section 4 is section 5 directly under upstream, therefore the third line 5th is classified as 1, and section 6 is section 4 directly under upstream, therefore the 4th of the 6th row is classified as 1.

Wherein, river amount, pollutant output quantity and each are entered according to the pollutant of the corresponding section of each control unit The upstream-downstream relationship of a section when establishing delay coefficient matrix, for process combination Fig. 2, i.e., obtains section 1, section respectively 2, section 3, section 4, section 5, the respective pollutant in section 6 enter river amount, pollutant output quantity, are getting each section respectively Pollutant enter river amount, pollutant output quantity and then combine each section upstream-downstream relationship, so that it may establish be detained coefficient Matrix.Further, after getting the respective pollutant in each section and entering river amount, pollutant output quantity, in conjunction with each The upstream-downstream relationship of section can calculate the corresponding pollutant in each section and be detained coefficient；Calculating each section Corresponding pollutant is detained coefficient and then combines the upstream-downstream relationship of each section, establishes and is detained coefficient matrix:

Wherein, the r being detained in coefficient matrix₁Indicate that the corresponding pollutant in section 1 is detained coefficient, r₂Indicate that section 2 is corresponding Pollutant is detained coefficient, r_nIndicate that n corresponding pollutant in section is detained coefficient, for Fig. 2, an available 6*6's is stagnant Coefficient matrix is stayed, such as is obtained

After obtaining transition probability matrix H and being detained coefficient matrix R, according to the transition probability matrix, the delay Coefficient matrix can establish pollutant transfer-coefficient matrix,Wherein,For pollutant transfer ratio square Battle array, H are transition probability matrix, and R is to be detained coefficient matrix, and I is unit matrix.

Then I-R are as follows:

Resulting pollutant transfer ratio square is then calculated according to formulaAre as follows:

Obtaining pollutant transfer-coefficient matrixAfterwards, being based on shouldEach upper reach can be determined to positioned at itself The pollutant contribution rate Gij of each lower reache in downstream, wherein Gij indicates upper reach i to positioned at itself lower reache j Pollutant contribution rate.For example, matrixIn the tertial coefficient 0.8 of the first row be pollutant tribute of the section 1 to section 3 Offer rate G13, matrixIn the tertial coefficient 0.8 of the second row be pollutant contribution rate G23 of the section 2 to section 3, matrixIn the third line the 4th arrange coefficient 0.75 be pollutant contribution rate G34 of the section 3 to section 4, matrixIn the 4th The coefficient 0.82 that row the 5th arranges is pollutant contribution rate G45 of the section 4 to section 5, matrixIn the 6th row the 4th column Coefficient 0.75 is pollutant contribution rate G64 of the section 6 to section 4.

By topological relation shown in Fig. 2 it is found that section 1, section 2 are section 3 directly under upstream, section 3 and section 6 It is section 4 directly under upstream, section 4 is section 5 directly under upstream, and section 1 and the section 2 middle span for section 4 Section 3 has been got over, it is second level upstream that being that second level upstream, similarly section 3 and section 6 are intermediate for section 5, which spans section 4, Similarly, section 1, section 2 correspond to for section 5, and centre, which is seen, spans section 4 and section 5, are three-level upstream.For directly under upper For trip, it is based on pollutant transfer-coefficient matrixCan directly obtain each upper reach to be located at itself directly under lower reache Pollutant contribution rate Gij, such as G13, G23, G34, G64, G45.And for second level upstream, to obtain second level upstream to it The pollutant contribution rate Gij in downstream, needs to carry out second iteration, so that it may obtain second level upstream and contribute pollutant downstream Rate Gij, such as obtains G14, G24, G35, G65.That is:

Wherein, matrixIn the coefficient 0.6 that arranges of the first row the 4th be section 1 to the pollutant contribution rate of section 4 G14, matrixIn the second row the 4th column coefficient 0.6 be pollutant contribution rate G24 of the section 2 to section 4, matrixIn the third line the 5th arrange coefficient 0.615 be pollutant contribution rate G35 of the section 3 to section 5, matrixIn The coefficients 0.615 of the 6th row the 5th column be section 6 to the pollutant contribution rate G65 of section 5.

And for three-level upstream, to obtain three-level upstream to pollutant contribution rate Gij downstream, need to carry out Iteration three times, so that it may obtain three-level upstream to pollutant contribution rate Gij downstream, such as obtain G15, G25.That is:

Wherein, matrixIn the first row the 5th arrange coefficient 0.492 be pollutant tribute of the section 1 to section 5 Offer rate G15, matrixIn the coefficients 0.492 of the second row the 5th column be section 2 to the pollutant contribution rate of section 5 G25。

Obtaining pollutant transfer-coefficient matrixAfterwards, it by iteration difference number, can finally obtain Pollutant contribution rate Gij to each upper reach to each lower reache for being located at itself downstream.Wherein, it needs to illustrate Be, above-mentioned the number of iterations only for topological relation shown in Fig. 2 for, different topological relations, the iteration needed Number is not only identical, such as, it is desirable to level Four upstream is calculated to pollutant contribution rate Gij downstream, it is necessary to iteration four times, example Such as, it is assumed that further include then calculating level Four upper reach 1 to section 7 positioned at the section 7 in 5 downstream of section on the basis of Fig. 2 Pollutant contribution rate G17, and when computationally trip section 2 is to the pollutant contribution rate G27 of section 7, it is necessary to iteration four times.

Dirt of each upper reach to each lower reache for being located at itself downstream can be obtained through the above way Contaminate object contribution rate Gij.

3rd embodiment

The present embodiment compared with first embodiment, the difference is that: it is described as shown in figure 4, after step s 104 Method further include: according to the environmental capacity of water that each control unit allows, determine the corresponding target stains source of each control unit Environmental capacity of water.It, can also be further right namely after completing to be allocated the environmental capacity of water that each control unit allows Target stains source in each control unit carries out the secondary distribution of environmental capacity of water.Certainly to the target stains in each control unit Source carries out the secondary distribution of environmental capacity of water, be also possible to the environmental capacity of water for completing to allow each control unit of distribution into Carried out after row amendment, namely the environmental capacity of water allowed based on revised each control unit to the target stains source in itself into The secondary distribution of row environmental capacity of water, as shown in Figure 5.

When carrying out the secondary distribution of environmental capacity of water to the target stains source in each control unit, acquisition may is that The pollutant of each default pollution sources in each control unit enters river amount；The environmental capacity of water allowed according to each control unit The pollutant of each default pollution sources corresponding with each control unit got enters river amount, determines that each control unit is corresponding Target stains source environmental capacity of water.For example, the water environment to the target stains source in the corresponding control unit 1 in section 1 is held When amount carries out secondary distribution, first obtains the in line life source pollutants in rural area in unit 1 processed and enter river amount Ns-i, obtain agricultural planting Source pollutants enter river amount Nz-i, obtain livestock and poultry cultivation source pollutants and enter river amount Nx-i, obtain industrial point source pollutant and enter river amount Gs- I obtains in line urban life source pollutants and enters river amount Cs-i；It is then based on the environmental capacity of water W of the permission of control unit 1₀(for For control unit 1, W₀The value range of=t*W1, coefficient t are 0.9-1), and Ns-i, Nz-i, Nx-i, Gs- for getting I, Cs-i can determine the environmental capacity of water in corresponding target stains source in control unit 1, for example, target stains source is rural area In line life source, then the environmental capacity of water W in the in line life source in rural area_s=W₀* (Ns-i/ei), similarly, when target stains source is agriculture Industry plants source, then the environmental capacity of water W in agricultural planting source_s=W₀* (Nz-i/ei), similarly, when target stains source is in line cities and towns Life source, then the environmental capacity of water W in line urban life source_s=W₀*(Cs-i/ei)。

For control unit 3, secondary distribution is carried out in the environmental capacity of water to the target stains source in control unit 3 When, if it is the environmental capacity of water W allowed based on the control unit 3 before amendment₀=t*W3 is allocated, such as to control unit When the environmental capacity of water in the target stains source such as agricultural planting source in 3 is allocated, W_s=t*W3*Nz-i/ei)；If it is base In the environmental capacity of water W that revised control unit 3 allows₀=W1*G13+W2*G23 is allocated, then in control unit 3 The environmental capacity of water in agricultural planting source when being allocated, W_s=(W1*G13+W2*G23) * (Nz-i/ei).In another example for For control unit 4, when being allocated to the environmental capacity of water in the target stains source such as livestock and poultry cultivation source in control unit 4, one Under kind embodiment, it can be, W_s=t*W4* (Nx-i/ei) under another embodiment, can be W_s=(W1*G14+W2* G24+W3*G34+W6*G64)*(Nx-i/ei)。

Wherein, it should be noted that above-mentioned only example holds the water environment in partial target pollution sources such as livestock and poultry cultivation sources The process being allocated is measured, the principle of the assigning process in different target stains sources is the same.In addition, according to each control The environmental capacity of water that unit allows can be base when determining the environmental capacity of water in the corresponding target stains source of each control unit The environmental capacity of water in the target stains source in it is allocated in the environmental capacity of water that control unit allows, a kind of embodiment Under, under mode as shown in Figure 5, at this point, being the environmental capacity of water based on the permission of revised control unit to the target in it The environmental capacity of water of pollution sources is allocated.

Wherein, ei=Ns-i+Nz-i+Nx-i+Cs-i+Gs-i.It should be noted that above-mentioned is using default pollution sources as agriculture The in line life source in village, agricultural planting source, livestock and poultry cultivation source, industrial point source, in line urban life source, corresponding target is dirty Dye source is the in line life source in rural area, agricultural planting source, livestock and poultry cultivation source, industrial point source, at least one in line urban life source Kind.

Wherein, when the pollutant for obtaining each default pollution sources in each control unit enters river and measures, it may is that acquisition The pollutant discharge amount of each default pollution sources in each control unit and each default corresponding river pollutant sources of pollution sources； According to the pollutant discharge amount of each default pollution sources in each control unit and each default pollution sources it is corresponding enter river Coefficient, the pollutant for obtaining each default pollution sources in the control unit enter river amount.For example, for control unit 2, Assuming that default pollution sources remain as above-mentioned 5 kinds, then the pollutant discharge amount Ns-out in the in line life source in rural area is obtained, and enter River system number a, it is assumed that be 0.1, obtain the pollutant discharge amount Nz-out and river pollutant sources b in agricultural planting source, it is assumed that it is 0.05, Obtain the pollutant discharge amount Nx-out and river pollutant sources c in livestock and poultry cultivation source, it is assumed that be 0.01, obtain in line urban life source Pollutant discharge amount Cs-ou and river pollutant sources d, it is assumed that be 1.0, obtain the pollutant discharge amount Gs-out of industrial point source with And river pollutant sources e, it is assumed that be 1.0.The pollutant that each default pollution sources in control unit 2 can then be obtained enters river Amount, such as Ns-i=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 also provides a kind of basin water environment capacity allocation devices 100, as shown in Figure 6.The basin water Environmental carrying capacity distributor 100 includes: that the first acquisition module 110, second obtains module 120, third obtains module 130 and the One determining module 140.

First obtains module 110, needs to carry out the basin to be allocated of environmental capacity of water distribution for obtaining.

Second obtains module 120, for obtaining each control unit relevant to the basin to be allocated, and with it is each The one-to-one section of control unit.Optionally, second module 120 is obtained, be specifically used for based on localized management principle to described Control unit in basin to be allocated is divided, and determines each corresponding section of control unit.

Third obtains module 130, the maximum environmental capacity of water for allowing to discharge for obtaining each section.

First determining module 140, the maximum environmental capacity of water for allowing to discharge based on each section, determines each section The environmental capacity of water that corresponding control unit allows.

Optionally, under a kind of embodiment, basin water environment capacity allocation device 100, further includes: correction module is used for The environmental capacity of water allowed the corresponding control unit in each section is modified.Optionally, correction module is specifically used for: really The value for determining the environmental capacity of water that the corresponding control unit in most upstream section allows is constant；Obtain other rivers in addition to the section of most upstream The pollutant of section is detained coefficient；Determine that pollutant is detained coefficient no more than the water ring that the control unit of the section of preset threshold allows The value of border capacity is constant；The water ring that the control unit that pollutant is detained the target section that coefficient is greater than the preset threshold is allowed Border capacity is modified.Optionally, the correction module, is specifically used for: obtaining all upstream rivers for being located at target section upstream Section；Each upper reach is obtained to the pollutant contribution rate of the target section；Maximum water environment based on each upper reach is held Amount and each upper reach obtain all upper reaches to the target river to the pollutant contribution rate of the target section Section environmental capacity of water contribution and；The environmental capacity of water that correcting the corresponding control unit in the target section allows is the water environment Capacity contribution and.

And basin water environment capacity allocation device 100 further include: the second determining module, for being used in the correction module It contributes for the environmental capacity of water with before, determines in correcting the environmental capacity of water that the corresponding control unit in the target section allows Environmental capacity of water environmental capacity of water contribution and allowed less than the target section.

And basin water environment capacity allocation device 100 further include: third determining module, for being used in the correction module Before obtaining each upper reach to the pollutant contribution rate of the target section, determines the environmental capacity of water contribution and be less than The environmental capacity of water that the target section allows.Optionally, the third determining module, is specifically used for: single according to each control The upstream-downstream relationship of the corresponding section of member, establishes transition probability matrix；According to the corresponding section of each control unit Pollutant enter the upstream-downstream relationship of river amount, pollutant output quantity and each section, establish and be detained coefficient matrix；According to described Transition probability matrix, the delay coefficient matrix establish pollutant transfer-coefficient matrix；Based on the pollutant transfer ratio square Battle array determines pollutant contribution rate of each upper reach to each lower reache for being located at itself downstream.Optionally, described Three determining modules, are specifically used for: entering river amount, pollutant output quantity according to the pollutant of the corresponding section of each control unit And the upstream-downstream relationship of each section, it calculates the corresponding pollutant in each section and is detained coefficient；It is each based on each section Self-corresponding pollutant is detained the upstream-downstream relationship of coefficient and each section, establishes and is detained coefficient matrix.

And basin water environment capacity allocation device 100 further include: computing module, for being used to obtain in the correction module The pollutant for removing other sections outside the section of most upstream is detained before coefficient, and it is stagnant to calculate the corresponding pollutant in each section Stay coefficient.Optionally, the computing module, is specifically used for: entering river amount, pollutant according to the corresponding pollutant in each section Output quantity and itself location calculate the corresponding pollutant in each section and are detained coefficient, wherein pollutant is detained The calculation formula of coefficient is r_i=(Load_in-Load_out)/Load_in, r_iIt is detained coefficient, Load for the corresponding pollutant of section i_in For the corresponding pollutant burden amount of section i, Loadin=subi+Load (i-1) out, subi are that i corresponding pollutant in section enters River amount, Load (i-1) out are section i directly under the corresponding pollutant output quantity of upstream i-1, and Loadout is that section i is corresponding Pollutant output quantity.

Under a kind of embodiment, optionally, basin water environment capacity allocation device 100 further include: the 4th determining module is used In after the environmental capacity of water that first determining module is used to determine that the corresponding control unit in each section to allow, according to every The environmental capacity of water that a control unit allows, determines the environmental capacity of water in the corresponding target stains source of each control unit.It is optional Ground, the 4th determining module, is specifically used for: the pollutant for obtaining each default pollution sources in each control unit enters river Amount；The environmental capacity of water each default pollution sources corresponding with each control unit got allowed according to each control unit Pollutant enter river amount, determine the environmental capacity of water in the corresponding target stains source of each control unit.Optionally, the described 4th really Cover half block, is specifically used for: obtaining the pollutant discharge amount of each default pollution sources in each control unit and each presets The corresponding river pollutant sources of pollution sources；According to the pollutant discharge amount of each default pollution sources in each control unit and each The default corresponding river pollutant sources of pollution sources, the pollutant for obtaining each default pollution sources in the control unit enter river amount.

Wherein, the target stains source be the in line life source in rural area, it is agricultural planting source, livestock and poultry cultivation source, industrial point source, straight Arrange at least one of urban life source.

It should be noted that all the embodiments in this specification are described in a progressive manner, each embodiment weight Point explanation is the difference from other embodiments, and the same or similar parts between the embodiments can be referred to each other.

The technology of basin water environment capacity allocation device 100 provided by the embodiment of the present application, realization principle and generation Effect is identical with preceding method embodiment, and to briefly describe, Installation practice part does not refer to place, can refer to preceding method reality Apply corresponding contents in example.

5th embodiment

The embodiment of the present application also provides a kind of electronic equipments that above-mentioned basin water environment capacity allocation method can be performed 200, as shown in Figure 7.The electronic equipment 200 includes: transceiver 210, memory 220, communication bus 230 and processor 240.

The transceiver 210, the memory 220, each element of processor 240 directly or indirectly electrically connect between each other It connects, to realize the transmission or interaction of data.For example, these elements between each other can by one or more communication bus 230 or Signal wire, which is realized, to be electrically connected.Wherein, transceiver 210 is used for sending and receiving data.Memory 220 is for storing computer program, such as It is stored with software function module shown in Fig. 6, i.e. basin water environment capacity allocation device 100.Wherein, basin water environment capacity Distributor 100 include at least one can be stored in the form of software or firmware (firmware) in the memory 220 or The software function module being solidificated in the operating system (operating system, OS) of the electronic equipment 200.The processing Device 240, for executing the executable module stored in memory 220, such as basin water environment capacity allocation device 100 includes Software function module or computer program.For example, processor 240, for obtain need to carry out environmental capacity of water distribution to point With basin；And be also used to obtain each control unit relevant to the basin to be allocated, and with each control unit one One corresponding section；And it is also used to obtain the maximum environmental capacity of water that each section allows to discharge；And it is based on each section Allow the maximum environmental capacity of water discharged, determines the environmental capacity of water that the corresponding control unit in each section allows.

Wherein, memory 220 may be, but not limited to, random access memory (Random Access Memory, RAM), read-only memory (Read Only Memory, ROM), programmable read only memory (Programmable Read-Only Memory, PROM), erasable read-only memory (Erasable Programmable Read-Only Memory, EPROM), Electricallyerasable ROM (EEROM) (Electric Erasable Programmable Read-Only Memory, EEPROM) etc..

Processor 240 may be a kind of IC chip, the processing capacity with signal.Above-mentioned processor can be General processor, including central processing unit (Central Processing Unit, CPU), network processing unit (Network Processor, NP) etc.；It can also be digital signal processor (Digital Signal Processor, DSP), dedicated integrated Circuit (Application Specific Integrated Circuit, ASIC), field programmable gate array (Field Programmable Gate Array, FPGA) either other programmable logic device, discrete gate or transistor logic, Discrete hardware components.It may be implemented or execute disclosed each method, step and the logic diagram in the embodiment of the present application.It is general Processor can be microprocessor or the processor 240 is also possible to any conventional processor etc..

Wherein, above-mentioned electronic equipment 200, including but not limited to computer, server etc..

Sixth embodiment

The embodiment of the present application also provides a kind of non-volatile computer read/write memory medium, (hereinafter referred to as storage is situated between Matter), computer program is stored on the storage medium, such as by above-mentioned electronic equipment when which is run by computer The step of basin water environment capacity allocation method that execution above method embodiment provides when 200 operation is included.

In several embodiments provided herein, it should be understood that disclosed device and method can also pass through Other modes are realized.The apparatus embodiments described above are merely exemplary, for example, flow chart and block diagram in attached drawing Show the device of multiple embodiments according to the application, the architectural framework in the cards of method and computer program product, Function and operation.In this regard, each box in flowchart or block diagram can represent the one of a module, section or code Part, a part of the module, section or code, which includes that one or more is for implementing the specified logical function, to be held Row instruction.It should also be noted that function marked in the box can also be to be different from some implementations as replacement The sequence marked in attached drawing occurs.For example, two continuous boxes can actually be basically executed in parallel, they are sometimes It can execute in the opposite order, this depends on the function involved.It is also noted that every in block diagram and or flow chart The combination of box in a box and block diagram and or flow chart can use the dedicated base for executing defined function or movement It realizes, or can realize using a combination of dedicated hardware and computer instructions in the system of hardware.

In addition, each functional module in each embodiment of the application can integrate one independent portion of formation together Point, it is also possible to modules individualism, an independent part can also be integrated to form with two or more modules.

It, can be with if the function is realized and when sold or used as an independent product in the form of software function module It is stored in a computer readable storage medium.Based on this understanding, the technical solution of the application is substantially in other words The part of the part that contributes to existing technology or the technical solution can be embodied in the form of software products, the meter Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be a People's computer, laptop, server or network equipment etc.) execute each embodiment the method for the application whole Or part steps.And storage medium above-mentioned include: USB flash disk, mobile hard disk, read-only memory (Read-Only Memory, ROM), random access memory (Random Access Memory, RAM), magnetic or disk etc. are various can store program The medium of code.

The above, the only specific embodiment of the application, but the protection scope of the application is not limited thereto, it is any Those familiar with the art within the technical scope of the present application, can easily think of the change or the replacement, and should all contain Lid is within the scope of protection of this application.Therefore, the protection scope of the application shall be subject to the protection scope of the claim.

Claims

1. a kind of basin water environment capacity allocation method characterized by comprising

Obtain the basin to be allocated for needing to carry out environmental capacity of water distribution；

Obtain each control unit relevant to the basin to be allocated, and with the one-to-one section of each control unit；

Obtain the maximum environmental capacity of water that each section allows to discharge；

Allow the maximum environmental capacity of water discharged based on each section, determines the water ring that the corresponding control unit in each section allows Border capacity；

The environmental capacity of water allowed the corresponding control unit in each section is modified.

2. the method according to claim 1, wherein the water environment allowed the corresponding control unit in each section Capacity is modified, comprising:

Determine that the value for the environmental capacity of water that the corresponding control unit in most upstream section allows is constant；

The pollutant for obtaining other sections in addition to the section of most upstream is detained coefficient；

It is constant no more than the value for the environmental capacity of water that the control unit of the section of preset threshold allows to determine that pollutant is detained coefficient；

The environmental capacity of water allowed the control unit that pollutant is detained the target section that coefficient is greater than the preset threshold carries out Amendment.

3. according to the method described in claim 2, it is characterized in that, being detained the mesh that coefficient is greater than the preset threshold to pollutant The environmental capacity of water that the control unit of mark section allows is modified, comprising:

Obtain all upper reaches for being located at target section upstream；

Each upper reach is obtained to the pollutant contribution rate of the target section；

Maximum environmental capacity of water and each upper reach based on each upper reach contribute the pollutant of the target section Rate, obtain all upper reaches to the environmental capacity of water contribution of the target section and；

Correct the environmental capacity of water that the corresponding control unit in the target section allows be the environmental capacity of water contribute and.

4. according to the method described in claim 3, it is characterized in that, correcting what the corresponding control unit in the target section allowed Environmental capacity of water be the environmental capacity of water contribute and before, the method also includes:

Determine environmental capacity of water environmental capacity of water contribution and allowed less than the target section.

5. according to the method described in claim 3, it is characterized in that, obtaining pollution of each upper reach to the target section Before object contribution rate, the method also includes:

Determine pollutant contribution rate of each upper reach to each lower reache for being located at itself downstream.

6. according to the method described in claim 5, it is characterized in that, determining each upper reach to positioned at each of itself downstream The pollutant contribution rate of a lower reache, comprising:

According to the upstream-downstream relationship of the corresponding section of each control unit, transition probability matrix is established；

Enter the upper of river amount, pollutant output quantity and each section according to the pollutant of the corresponding section of each control unit Downstream relationship is established and is detained coefficient matrix；

Pollutant transfer-coefficient matrix is established according to the transition probability matrix, the delay coefficient matrix；

Based on the pollutant transfer-coefficient matrix, determine each upper reach to each lower reache for being located at itself downstream Pollutant contribution rate.

7. according to the method described in claim 6, it is characterized in that, according to the pollution of the corresponding section of each control unit Object enters river amount, pollutant output quantity and the upstream-downstream relationship of each section, establishes and is detained coefficient matrix, comprising:

Enter the upper of river amount, pollutant output quantity and each section according to the pollutant of the corresponding section of each control unit Downstream relationship calculates the corresponding pollutant in each section and is detained coefficient；

It is detained the upstream-downstream relationship of coefficient and each section based on the corresponding pollutant in each section, establishes and be detained coefficient Matrix.

8. a kind of basin water environment capacity allocation device characterized by comprising

First obtains module, needs to carry out the basin to be allocated of environmental capacity of water distribution for obtaining；

Second obtains module, for obtaining each control unit relevant with the basin to be allocated, and singly with each control The one-to-one section of member；

Third obtains module, the maximum environmental capacity of water for allowing to discharge for obtaining each section；

First determining module, the maximum environmental capacity of water for allowing to discharge based on each section, determines that each section is corresponding The environmental capacity of water that control unit allows；

Correction module, the environmental capacity of water for allowing the corresponding control unit in each section are modified.

9. a kind of electronic equipment characterized by comprising memory and processor, the memory are connected with the processor；

The memory is for storing program；

The processor is for calling the program being stored in the memory to execute as described in any one of claim 1-7 Method.

10. a kind of storage medium, which is characterized in that be stored thereon with computer program, the computer program is transported by computer Such as method of any of claims 1-7 is executed when row.