CN106203688A - A kind of economy and the optimal screening method of point-source pollution abatement technology under environmental constraints - Google Patents

Abstract

The invention discloses a kind of economy and the optimal screening method of point-source pollution abatement technical measures under environmental constraints, including: carry out point-source pollution load prediction under conventional development pattern；Water Functional Zone is divided to propose water environmental model constraints；Identify the point-source pollution decreasing measure set of technical feasibility；Set up region social economy constraints；Determine the cost function that point-source pollution is cut down；Determine the cost function that point-source pollution is cut down.The present invention enters river amount dynamically based on point source pollutant, consider the double constraints condition such as economic development total amount, water correction target, compared with cost function by the system identification of point-source pollution abatement technology, it is achieved that the optimal screening of multiple point-source pollution technical combinations；The point-source pollution abatement engineering construction carrying out Technological Economy feasible for basin provides technical support.

Description

A kind of economy and the optimal screening method of point-source pollution abatement technology under environmental constraints

Technical field

The invention belongs to point-source pollution technical field, particularly relate to a kind of economy and point-source pollution abatement skill under environmental constraints The optimal screening method of art.

Background technology

Traditional point-source pollution abatement, generally according to reducing discharging proportion requirement, proposes multinomial point-source pollution abatement technology (such as dirt Water processes, cleans production etc.), there are following two aspect shortcomings: (1) these abatement technology are generally not set up with quality of water environment Quantitative response relation, mainly based on the control of qualified discharge, can be real to the point-source pollution abatement technology planned or have been carried out Existing water quality objective, it is impossible to the effect assessment of quantitative；(2) point-source pollution of the lower economic optimum of environmental goals constraint cannot be realized Abatement technology screening, cause point source cut down technological investment high cost, the disengagement zone socio-economic development stage and cannot be effective Implement.

Summary of the invention

Present invention aim at providing a kind of economy and the optimal screening method of point-source pollution abatement technology under environmental constraints, Aim to solve the problem that the point-source pollution abatement technology screening that cannot realize the lower economic optimum of environmental goals constraint, cause point-source pollution to be cut down Technological investment high cost, disengagement zone socio-economic development stage and the problem that cannot effectively implement.

For reaching above-mentioned purpose, the present invention adopts the following technical scheme that

A kind of economy and the optimal screening method of point-source pollution abatement technology under environmental constraints, including step:

Step 1, under the conditions of carrying out point-source pollution load prediction under the conditions of present situation abatement, i.e. quantitative forecast present situation abatement, respectively What in the industrial wastewater of Water Functional Zone and sanitary sewage, each pollutant were inscribed when each enters river amount；

Step 2, divides a Water Functional Zone unit to establish water quality objective, and Con trolling index requirement is measured in the river that enters of point source pollutant of retrodicting, I.e. water quality objective constraints, this step particularly as follows:

Water function area dividing according to different rivers and water quality objective, water quality condition because of the present circumstance, build quality of water environment mould Type；According to quality of water environment model, point Water Functional Zone retrodict reach the constraints of water quality objective, i.e. pollutant enter river amount control Index request WG processed_i,k,t；Water Functional Zone to river, measures Con trolling index with COD and ammonia nitrogen for entering river；To lake and The Water Functional Zone of reservoir, with COD, ammonia nitrogen, total nitrogen and total phosphorus for entering river amount Con trolling index；

Step 3, the river that enters that river is measured and step 2 is retrodicted that enters predicted based on step 1 is measured Con trolling index requirement, is identified each water The point-source pollution abatement technology set of the technical feasibility of functional areas；Described point-source pollution abatement technology set includes contamination sources Head emission-reduction technology, pollutant process interrupter technique and/or pollutant end-of-pipe control technology, and the pollution that each abatement technology is corresponding Thing clearance parameter, economic characteristics, applicable elements；

Step 4, sets up the investment and recovery condition of Water Functional Zone unit spot source pollutants abatement, this step particularly as follows:

Based on statistical analysis and Correlative plan, identify the present situation of social economic development of region and following scale, determine point source Pollute the maximum possible scale of investment of abatement, set up scale of investment and the ratio of regional economy overall size of point-source pollution abatement Example, thus obtain point source pollutant write-down investment constraint KK_i,t, the i.e. Water Functional Zone unit of described regional economy overall size GDP value；

Step 5, determines the unit cost functions of point-source pollution abatement technology set Point Source every pollution abatement technology, point Source is polluted the unit cost functions of abatement technology and is cut down capital expenditure and the operation and maintenance cost structure of technology by point-source pollution Become；

Step 6, binding site source pollutants write-down investment constraint KK_i,tUnit cost functions with point-source pollution abatement technology C_i,k,t, from point-source pollution abatement technology set filter out optimum point-source pollution abatement technology, this step particularly as follows:

Retrain based on unit water quality objective constraints and unit spot source pollutants write-down investment and build constraints, from respectively In the point-source pollution abatement technology set of Water Functional Zone, optimal screening goes out the optimum point-source pollution abatement technology of each Water Functional Zone.

In step 1, under the conditions of present situation abatement in the industrial wastewater of each Water Functional Zone and sanitary sewage each pollutant when each The river amount that enters inscribed uses equation below prediction, wherein unit i.e. Water Functional Zone unit:

M_i,k,t=D_i,k,t×K_i,k,t (1)

D_i,k,t=WP_i,k,t×SP_i,k,t (2)

WP_i,k,t=(WDP_i,t-RDP_i,t)×CW_i,k,t×10^-2+(WIP_i,t-RIP_i,t)×CW_i,k,t×10^-2 (3)

${\mathrm{CW}}_{i,k,t}=min({\mathrm{CS}}_k,\frac{({\mathrm{WDP}}_{i,t}\×{\mathrm{PDP}}_{i,k,t}+{\mathrm{WIP}}_{i,t}\×{\mathrm{PIP}}_{i,k,t})(1-L_k\×{10}^{-2})}{{\mathrm{WDP}}_{i,t}+{\mathrm{WIP}}_{i,t}})---\left(4\right)$

SP_i,k,t=(QDP_i,t-WDP_i,t)×PDP_i,k,t+(QIP_i,t-WIP_i,t)×PIP_i,k,t (5)

In formula (1)～(5):

M_i,k,tRepresent i unit k pollutant t enters river amount, unit: t/a；

D_i,k,tRepresent the discharge capacity of i unit k pollutant t, unit: t/a；

K_i,k,tRepresent the river pollutant sources of i unit k pollutant t；

WP_i,k,tRepresent the discharge capacity of the wastewater reuse approach factory k pollutant of i unit t, unit: t/a；

SP_i,k,tRepresent the in line amount of k pollutant, unit: t/a in the polluter of i unit t；

WDP_i,tRepresent the wastewater reuse approach factory sanitary sewage disposal amount of i unit t, unit: ten thousand m³/a；

RDP_i,tRepresent reuse amount after the wastewater reuse approach factory sanitary sewage disposal of i unit t, unit: ten thousand m³/a；

CW_i,k,tRepresent the Tail water reuse concentration of the wastewater reuse approach factory k pollutant of i unit t, unit: mg/l；

WIP_i,tRepresent the wastewater reuse approach factory Industrial Wastewater Treatment amount of i unit t, unit: ten thousand m³/a；

RIP_i,tRepresent reuse amount after the wastewater reuse approach factory Industrial Wastewater Treatment of i unit t, unit: ten thousand m³/a；

PDP_i,k,tRepresent the generation concentration of k pollutant, unit: mg/l in i unit t sanitary sewage；

PIP_i,k,tRepresent the generation concentration of k pollutant, unit: mg/l in i unit t industrial wastewater；

L_kRepresent the clearance of wastewater reuse approach factory k pollutant, unit: %；

CS_kRepresent the discharge standard of sewage treatment plant's k pollutant, unit: mg/l；

QDP_i,tRepresent the sanitary sewage discharge capacity of i unit t, unit: ten thousand m³/a；

QIP_i,tRepresent the discharged volume of industrial waste water of i unit t, unit: ten thousand m³/a。

In step 2, described quality of water environment model is zero dimension, one or more dimensions form.

In step 2, pollutant enter river overall control index request WG_i,k,t=Min (M_i,k,t,CPT2_i,k,t), wherein:

WG_i,k,tCon trolling index requirement, unit: t/a are measured in the river that enters for i unit t k pollutant；

M_i,k,tThe river that enters for i unit k pollutant t is measured, unit: t/a；

CPT2_i,k,tFor the Water Functional Zone water environment capacity of i unit t k pollutant, unit: t/a.

In step 4, point source pollutant write-down investment constraint KK_i,t=GDP_i,t·GK_i,t, wherein:

KK_i,tPollutant write-down investment for i unit t retrains, unit: Wan Yuan；

GDP_i,tFor the regional economy overall size of i unit t, the i.e. GDP of Water Functional Zone unit, unit: Wan Yuan；

GK_i,tThe scale of investment cut down for i unit t point-source pollution and the ratio of regional economy overall size are immeasurable Guiding principle.

In step 4, the scale of investment of point-source pollution abatement is 0.5%-3.0% with the ratio of regional economy overall size.

In step 5, the unit cost functions of point-source pollution abatement technology:

$C_{i,k,t}=\frac{Dr\×{\mathrm{Inv}}_{i,k,t}\×R{(1+R)}^{Lf}}{{(1+R)}^{Lf}-1}+{\mathrm{Ope}}_{i,k,t}$

Wherein:

C_i,k,tFor the totle drilling cost of i unit t abatement unit quantity k pollutant, unit: unit/year；

Inv_i,k,tFor the cost of investment of i unit t abatement unit quantity k pollutant, unit: unit/year；

Ope_i,k,tFor the operating cost of i unit t abatement unit quantity k pollutant, unit: unit/year；

Dr is cost of investment coefficient of depreciation, unit: 1/ year；

R is depreciation interest rate；

Lf is that point-source pollution abatement is arranged service life of technology, unit: year.

In step 6, described object function is Max (W1_i,k,t+W2_i,k,t+W3_i,k,t), wherein:

W1_i,k,tThe k pollutant abatement amount realized for i unit t enterprise qualified discharge, unit: t/a；

W2_i,k,tThe k pollutant abatement amount realized for i unit t wastewater reuse approach factory, unit: t/a；

W3_i,k,tThe k pollutant abatement amount of realization, unit: t/a is cut down the most further for i unit t.

In step 6, described constraints includes:

C1_i,k,t×W1_i,k,t+C2_i,k,t×W2_i,k,t+C3_i,k,t×W3_i,k,t≤KK_i,t

$W{1}_{i,k,t}+W{2}_{i,k,t}+W{3}_{i,k,t}\≤{\mathrm{PG}}_{i,k,t}-\frac{{\mathrm{WG}}_{i,k,t}}{K_{i,k,t}}$

W1_{I, k, t}≤LW1_{I, k, t}

W2_{I, k, t}≤LW2_{I, k, t}

W1_{I, k, t}≥0

W2_{I, k, t}≥0

W3_{I, k, t}≥0

Wherein:

C1_i,k,tFor the totle drilling cost of the abatement unit quantity k pollutant that i unit t wastewater reuse approach mode realizes, list Position: ten thousand yuan/t；

C2_i,k,tAssembly for the abatement unit quantity k pollutant that i unit t polluter dispersion up to standard administration way realizes This, unit: ten thousand yuan/t；

C3_i,k,tFor the totle drilling cost of the abatement unit quantity k pollutant that i unit t realizes the most further, list Position: ten thousand yuan/t；

W1_i,k,tThe k pollutant abatement amount realized for i unit t enterprise qualified discharge, unit: t/a；

W2_i,k,tThe k pollutant abatement amount realized for i unit t wastewater reuse approach factory, unit: t/a；

W3_i,k,tThe k pollutant abatement amount of realization, unit: t/a is cut down the most further for i unit t；

KK_i,tPollutant write-down investment for i unit t retrains, unit: Wan Yuan；

PG_i,k,tFor the generation amount of i unit t k pollutant, unit: t/a；

WG_i,k,tCon trolling index requirement, unit: t/a are measured in the river that enters for i unit t k pollutant；

K_i,k,tRiver pollutant sources for i unit t k pollutant；

LW1_i,k,tFor the i unit t enterprise maximum attainable k pollutant abatement amount of qualified discharge, unit: t/a；

LW2_i,k,tFor the maximum attainable k pollutant abatement amount of i unit t wastewater reuse approach factory, unit: t/a.

Compared to the prior art, the invention have the advantages that and beneficial effect:

(1) comprehensive by force, considered economic development total amount and water correction target double constraints condition, passed through point source The system identification polluting abatement technology is compared with cost function, it is achieved that the optimization sieve of multiple point-source pollution abatement technical combinations Choosing.

(2) workable, the point-source pollution abatement engineering construction carrying out Technological Economy feasible for basin provides decision-making to prop up Hold and technical support.

Accompanying drawing explanation

Fig. 1 is the idiographic flow schematic diagram of the inventive method；

Fig. 2 is the embodiment representative basin COD pollutant water environment capacity provided and the schematic diagram entering river amount；

Fig. 3 is the one-tenth of representative basin different units each point-source pollution abatement technical units COD abatement amount that embodiment provides This schematic diagram compared.

Detailed description of the invention

In order to make the purpose of the present invention, technical scheme and advantage clearer, below in conjunction with embodiment, to the present invention It is further elaborated.Should be appreciated that specific embodiment described herein, only in order to explain the present invention, is not used to Limit the present invention.

Below in conjunction with the accompanying drawings and the application principle of the present invention is further described by specific embodiment.

As it is shown in figure 1, the economy of the embodiment of the present invention and the optimal screening side of point-source pollution abatement technology under environmental constraints Method comprises the following steps:

S101: carry out point-source pollution load prediction under the conditions of present situation abatement；

S102: a point Water Functional Zone unit establishes water quality objective, point source pollutant of retrodicting enters river amount Con trolling index；

S103: according to step S101 and S102 result, identify the point-source pollution abatement technology set of technical feasibility；

S104: set up the investment and recovery condition of computing unit point source pollutant abatement；

S105: according to step S103 result, determine the cost function of point source every abatement technology；

S106: integrating step S104 and S105, propose the economic optimization screening scheme of point-source pollution abatement technology.

The optimal screening method of the present invention be described more detail below:

Step one, carries out point-source pollution load prediction under the conditions of present situation abatement.

Abatement condition refers to sewage treatment facility and waste water treating and reutilizing amount condition.In this step, consider the dirt of present situation Water processing establishment, waste water treating and reutilizing amount and the impact of different river pollutant sources, the industry of the following each Water Functional Zone of quantitative forecast What in waste water and sanitary sewage, each pollutant were inscribed when each enters river amount.During prediction, the most only need to consider the sewage of present situation Treatment facility and waste water treating and reutilizing amount, provide basis for cutting down the laying of technical scheme future.River pollutant sources size is by multiple Factor affects, such as discharging modes (such as pipeline or channel), discharge environment (such as temperature) etc., and usually 0.3～0.9, need to be according to reality Border situation carries out measuring.

The fundamental formular of point-source pollution load prediction is as follows, and wherein i unit i.e. represents i-th Water Functional Zone unit:

M_i,k,t=D_i,k,t×K_i,k,t (1)

D_i,k,t=WP_i,k,t×SP_i,k,t (2)

WP_i,k,t=(WDP_i,t-RDP_i,t)×CW_i,k,t×10^-2+(WIP_i,t-RIP_i,t)×CW_i,k,t×10^-2 (3)

${\mathrm{CW}}_{i,k,t}=min({\mathrm{CS}}_k,\frac{({\mathrm{WDP}}_{i,t}\×{\mathrm{PDP}}_{i,k,t}+{\mathrm{WIP}}_{i,t}\×{\mathrm{PIP}}_{i,k,t})(1-L_k\×{10}^{-2})}{{\mathrm{WDP}}_{i,t}+{\mathrm{WIP}}_{i,t}})---\left(4\right)$

SP_i,k,t=(QDP_i,t-WDP_i,t)×PDP_i,k,t+(QIP_i,t-WIP_i,t)×PIP_i,k,t (5)

In formula (1)～(5):

M_i,k,tRepresent i unit k pollutant t enters river amount, unit: t/a；

D_i,k,tRepresent the discharge capacity of i unit k pollutant t, unit: t/a；

K_i,k,tRepresent the river pollutant sources of i unit k pollutant t；

WP_i,k,tRepresent the discharge capacity of the wastewater reuse approach factory k pollutant of i unit t, unit: t/a；

SP_i,k,tRepresent the in line amount of k pollutant, unit: t/a in the polluter of i unit t；

WDP_i,tRepresent the wastewater reuse approach factory sanitary sewage disposal amount of i unit t, unit: ten thousand m³/a；

RDP_i,tRepresent reuse amount after the wastewater reuse approach factory sanitary sewage disposal of i unit t, unit: ten thousand m³/a；

CW_i,k,tRepresent the Tail water reuse concentration of the wastewater reuse approach factory k pollutant of i unit t, unit: mg/l；

WIP_i,tRepresent the wastewater reuse approach factory Industrial Wastewater Treatment amount of i unit t, unit: ten thousand m³/a；

RIP_i,tRepresent reuse amount after the wastewater reuse approach factory Industrial Wastewater Treatment of i unit t, unit: ten thousand m³/a；

PDP_i,k,tRepresent the generation concentration of k pollutant, unit: mg/l in i unit t sanitary sewage；

PIP_i,k,tRepresent the generation concentration of k pollutant, unit: mg/l in i unit t industrial wastewater；

L_kRepresent the clearance of wastewater reuse approach factory k pollutant, unit: %；

CS_kRepresent the discharge standard of sewage treatment plant's k pollutant, unit: mg/l；

QDP_i,tRepresent the sanitary sewage discharge capacity of i unit t, unit: ten thousand m³/a；

QIP_i,tRepresent the discharged volume of industrial waste water of i unit t, unit: ten thousand m³/a。

Step 2, according to Water Functional Zone water quality objective retrodict pollutant enter river amount Con trolling index requirement.

This step is according to different river water function area dividing and present situation water quality condition, and point Water Functional Zone proposes water quality objective constraint Condition, the river amount Con trolling index that enters of the most each Water Functional Zone pollutant (such as COD and ammonia nitrogen) requires WG_i,k,t。

Enter river amount Con trolling index requirement WG_i,k,tBe calculated as follows:

WG_i,k,t=Min (M_i,k,t,CPT2_i,k,t) (6)

In formula (6):

WG_i,k,tCon trolling index requirement, unit: t/a are measured in the river that enters for i unit t k pollutant；

M_i,k,tThe river that enters for i unit k pollutant t is measured, unit: t/a, M_i,k,tCalculated by formula (1) and obtain；

CPT2_i,k,tFor the Water Functional Zone water environment capacity of i unit t k pollutant, unit: t/a.

This step particularly as follows:

Water function area dividing according to different rivers and water quality objective, water quality condition because of the present circumstance, build quality of water environment mould Type；According to quality of water environment model, retrodict and reach the constraints of water quality objective in point Water Functional Zone, the most each Water Functional Zone pollutant Con trolling index requirement is measured in the river that enters of (such as COD and ammonia nitrogen).Water Functional Zone to river, measures for entering river with COD and ammonia nitrogen Con trolling index；To lake and the Water Functional Zone of reservoir, with COD, ammonia nitrogen, total nitrogen and total phosphorus for entering river amount Con trolling index. In the present invention, according to reliability and the level of detail of data, quality of water environment model can use zero dimension, one or more dimensions form.

Fig. 2 is representative basin COD (COD) water environment capacity and the schematic diagram entering river amount.

Step 3, identifies the point-source pollution abatement technology set of technical feasibility.

The river that enters of each Water Functional Zone pollutant based on step one prediction is measured, and each Water Functional Zone that step 2 is retrodicted is dirty Con trolling index requirement is measured in the river that enters of dye thing, identifies the point-source pollution abatement technology set of the technical feasibility of each Water Functional Zone.

This set should consider the source emission-reduction technology of pollutant, and the process of pollutant to be considered blocks and end treatment Technology；Source emission-reduction technology mainly considers clearer production technology popularization, industry restructuring, recycling economy foundation etc., and process hinders Disconnected technology mainly considers perfect, the construction of sewage treatment plant of sewage collecting official website, the structure of regenerated water utilizing system, and end is controlled Reason technology mainly considers wetland construction, sludge handling, water ecology reparation etc..Need to set up abatement skill for various abatement technology Art data base, stores the pollutants removal rate parameter of various abatement technology, economic characteristics, applicable elements etc., for technology below Condition is preferably provided.

Step 4, sets up the investment and recovery condition of Water Functional Zone unit spot source pollutants abatement.

Based on statistical analysis and Correlative plan, identify the present situation of social economic development of region and following scale, determine point source Pollute the maximum possible scale of investment of abatement, set up point-source pollution abatement scale of investment and regional economy overall size (as GDP) ratio, thus obtain point source pollutant write-down investment constraint KK_i,t。

KK_i,tBe calculated as follows:

KK_i,t=GDP_i,t·GK_i,t (7)

In formula (7):

KK_i,tPollutant write-down investment for i unit t retrains, unit: Wan Yuan；

GDP_i,tFor the regional economy overall size of i unit t, the i.e. GDP of Water Functional Zone unit, unit: Wan Yuan；

GK_i,tThe scale of investment cut down for i unit t point-source pollution and the ratio of regional economy overall size are immeasurable Guiding principle.

Generally, point-source pollution abatement, as a part for environmental infrastructure investment, generally accounts for this regional economy scale The 0.5%-3.0% of total amount, i.e. GK_i,tGenerally take 0.5%-3.0%；Socio-economic development present situation according to discrete cell, location With trend, this ratio GK_i,tCan suitably adjust and change.

Step 5, the point-source pollution abatement technology set obtained according to step 3, determine that point source is every and pollute abatement skill The unit cost functions of art.

The unit cost of point-source pollution abatement technology set Point Source every pollution abatement technology is built based on investigation and analysis Function.Assume the cost of point source every pollution abatement technology and abatement nothing impact, the unit cost letter each other on pollutant Number had both needed to consider year depreciable cost, it is also contemplated that annual operating and maintenance cost.Concrete, the cost of point-source pollution abatement technology is mainly by base This construction cost and operation and maintenance cost two parts are constituted.Capital expenditure concentrates expenditure the most in the short term, such as engineering Reconnoitre, survey and draw, design, civil engineering and installation cost etc.；Operation and maintenance cost generally year-by-year payment within the project life cycle, such as people Work, power consumption and plant maintenance expense etc..Totle drilling cost is equal to cost of investment depreciation and operation and maintenance cost sum.For considering time valency Being worth, the cost of investment depreciation of water appliance is generally converted according to initial outlay cost single-candidate present worth in its life cycle, and can Can be affected by the water saving device subsidy for capital expenditure policy of government.Typically, higher interest rate and shorter abatement technology implementation week Phase, the totle drilling cost causing abatement technology is increased.

Seeing Fig. 3, the cost that representative basin point source difference pollutes abatement technical units COD abatement amount is as follows:

$C_{i,k,t}=\frac{Dr\×{\mathrm{Inv}}_{i,k,t}\×R{(1+R)}^{Lf}}{{(1+R)}^{Lf}-1}+{\mathrm{Ope}}_{i,k,t}---\left(8\right)$

In formula (8):

C_i,k,tFor the totle drilling cost of i unit t abatement unit quantity k pollutant, unit: unit/year；

Inv_i,k,tFor the cost of investment of i unit t abatement unit quantity k pollutant, unit: unit/year；

Ope_i,k,tFor the operating cost of i unit t abatement unit quantity k pollutant, unit: unit/year；

Dr is cost of investment coefficient of depreciation, unit: 1/ year；

R is depreciation interest rate；

Lf is the service life of abatement technology, unit: year.

Step 6, binding site source pollutants write-down investment constraint KK_i,tUnit cost functions with point-source pollution abatement technology C_i,k,t, treatment in accordance with local conditions optimal screening goes out the point-source pollution abatement technology that each Water Functional Zone is optimum, thus the point source obtaining optimum is dirty Dye abatement technical combinations.

The cost absorbing and benefit of quantitative analysis point-source pollution abatement technology, is calculated by optimal screening under investment and recovery The point-source pollution abatement technical combinations minimizing cost meeting economy with environmental goals is proposed.River amount Con trolling index is entered in difference On the basis of being separately optimized, take each optimization unit investment enveloping outer enclosure, expand to whole basin aspect, propose basin aspect different The recommendation abatement technical scheme of unit.

In this step, it is object function to the maximum, based on unit with the pollutant abatement amount that point-source pollution abatement technology realizes Water quality objective constraints and the constraint of unit spot source pollutants write-down investment build constraints, and the point source from each Water Functional Zone is dirty In dye abatement technology set, optimal screening goes out the optimum point-source pollution abatement technology of each Water Functional Zone.

Object function is shown in formula (9):

Max(W1_i,k,t+W2_i,k,t+W3_i,k,t) (9)

Constraints is shown in formula (10)～(16):

C1_i,k,t×W1_i,k,t+C2_i,k,t×W2_i,k,t+C3_i,k,t×W3_i,k,t≤KK_i,t (10)

$W{1}_{i,k,t}+W{2}_{i,k,t}+W{3}_{i,k,t}\≤{\mathrm{PG}}_{i,k,t}-\frac{{\mathrm{WG}}_{i,k,t}}{K_{i,k,t}}---\left(11\right)$

W1_{I, k, t}≤LW1_{I, k, t} (12)

W2_{I, k, t}≤LW2_{I, k, t} (13)

W1_{I, k, t}≥0 (14)

W2_{I, k, t}≥0 (15)

W3_{I, k, t}≥0 (16)

In formula (9)～(16):

W1_i,k,tThe k pollutant abatement amount realized for i unit t enterprise qualified discharge, unit: t/a；

W2_i,k,tThe k pollutant abatement amount realized for i unit t wastewater reuse approach factory, unit: t/a；

W3_i,k,tThe k pollutant abatement amount of realization, unit: t/a is cut down the most further for i unit t；

C1_i,k,tFor the totle drilling cost of the abatement unit quantity k pollutant that i unit t wastewater reuse approach mode realizes, list Position: ten thousand yuan/t；

C2_i,k,tAssembly for the abatement unit quantity k pollutant that i unit t polluter dispersion up to standard administration way realizes This, unit: ten thousand yuan/t；

C3_i,k,tFor the totle drilling cost of the abatement unit quantity k pollutant that i unit t realizes the most further, list Position: ten thousand yuan/t；

PG_i,k,tFor the generation amount of i unit t k pollutant, unit: t/a；

K_i,k,tRiver pollutant sources for i unit t k pollutant；

LW1_i,k,tFor the i unit t enterprise maximum attainable k pollutant abatement amount of qualified discharge, unit: t/a；

LW2_i,k,tFor the maximum attainable k pollutant abatement amount of i unit t wastewater reuse approach factory, unit: t/a；

C1_i,k,t、C2_i,k,t、C3_i,k,tUse formula (8) to calculate to obtain.

The inventive method is comprehensive by force, and system considers the double constraints condition such as economic development total amount, water correction target；Logical The system identification crossing point-source pollution abatement technology is compared with cost function, it is achieved the optimization sieve of multiple point-source pollution technical combinations Choosing；Workable, can be that the point-source pollution abatement engineering construction that Technological Economy is carried out feasible in basin provides technical support etc. excellent Point, it is adaptable in the range of the large scale of basin, Soil erosion class pollution of area source enters the calculating of river amount；

The foregoing is only presently preferred embodiments of the present invention, not in order to limit the present invention, all essences in the present invention Any amendment, equivalent and the improvement etc. made within god and principle, should be included within the scope of the present invention.

Claims (9)

1. economy and an optimal screening method for point-source pollution abatement technology under environmental constraints, is characterized in that, including:

Step 1, under the conditions of carrying out point-source pollution load prediction under the conditions of present situation abatement, i.e. quantitative forecast present situation abatement, each water merit What in the industrial wastewater in energy district and sanitary sewage, each pollutant were inscribed when each enters river amount；

Step 2, point Water Functional Zone unit establishment water quality objective, Con trolling index requirement, i.e. water are measured in the river that enters of point source pollutant of retrodicting Matter targeting constraints, this step particularly as follows:

Water function area dividing according to different rivers and water quality objective, water quality condition because of the present circumstance, build quality of water environment model；Root According to quality of water environment model, point Water Functional Zone retrodict reach the constraints of water quality objective, i.e. pollutant enter river amount control refer to Mark requires WG_i,k,t；Water Functional Zone to river, measures Con trolling index with COD and ammonia nitrogen for entering river；To lake and reservoir Water Functional Zone, with COD, ammonia nitrogen, total nitrogen and total phosphorus for enter river amount Con trolling index；

Step 3, the river that enters that river is measured and step 2 is retrodicted that enters predicted based on step 1 is measured Con trolling index requirement, is identified each water function The point-source pollution abatement technology set of the technical feasibility in district；Described point-source pollution abatement technology set includes that pollutant source subtracts Drainage technique, pollutant process interrupter technique and/or pollutant end-of-pipe control technology, and pollutant corresponding to each abatement technology go Except rate parameter, economic characteristics, applicable elements；

Step 4, sets up the investment and recovery condition of Water Functional Zone unit spot source pollutants abatement, this step particularly as follows:

Based on statistical analysis and Correlative plan, identify the present situation of social economic development of region and following scale, determine point-source pollution The maximum possible scale of investment of abatement, sets up the scale of investment of point-source pollution abatement and the ratio of regional economy overall size, from And obtain point source pollutant write-down investment constraint KK_i,t, the described regional economy overall size i.e. GDP value of Water Functional Zone unit；

Step 5, determines the unit cost functions of point-source pollution abatement technology set Point Source every pollution abatement technology, and point source is dirty The unit cost functions of dye abatement technology is made up of capital expenditure and the operation and maintenance cost of point-source pollution abatement technology；

Step 6, binding site source pollutants write-down investment constraint KK_i,tUnit cost functions C with point-source pollution abatement technology_i,k,t, From point-source pollution abatement technology set filter out optimum point-source pollution abatement technology, this step particularly as follows:

Retrain based on unit water quality objective constraints and unit spot source pollutants write-down investment and build constraints, from each water merit In the point-source pollution abatement technology set in energy district, optimal screening goes out the optimum point-source pollution abatement technology of each Water Functional Zone.

2. economy as claimed in claim 1 and the optimal screening method of point-source pollution abatement technology, its feature under environmental constraints It is:

In step 1, under the conditions of present situation abatement, in the industrial wastewater of each Water Functional Zone and sanitary sewage, each pollutant are inscribed when each Enter river amount use equation below prediction, wherein unit i.e. Water Functional Zone unit:

M_i,k,t=D_i,k,t×K_i,k,t (1)

D_i,k,t=WP_i,k,t×SP_i,k,t (2)

WP_i,k,t=(WDP_i,t-RDP_i,t)×CW_i,k,t×10^-2+(WIP_i,t-RIP_i,t)×CW_i,k,t×10^-2 (3)

${\mathrm{CW}}_{i,k,t}=min({\mathrm{CS}}_k,\frac{({\mathrm{WDP}}_{i,t}\×{\mathrm{PDP}}_{i,k,t}+{\mathrm{WIP}}_{i,t}\×{\mathrm{PIP}}_{i,k,t})(1-L_k\×{10}^{-2})}{{\mathrm{WDP}}_{i,t}+{\mathrm{WIP}}_{i,t}})---\left(4\right)$

SP_i,k,t=(QDP_i,t-WDP_i,t)×PDP_i,k,t+(QIP_i,t-WIP_i,t)×PIP_i,k,t (5)

In formula (1)～(5):

M_i,k,tRepresent i unit k pollutant t enters river amount, unit: t/a；

D_i,k,tRepresent the discharge capacity of i unit k pollutant t, unit: t/a；

K_i,k,tRepresent the river pollutant sources of i unit k pollutant t；

WP_i,k,tRepresent the discharge capacity of the wastewater reuse approach factory k pollutant of i unit t, unit: t/a；

SP_i,k,tRepresent the in line amount of k pollutant, unit: t/a in the polluter of i unit t；

WDP_i,tRepresent the wastewater reuse approach factory sanitary sewage disposal amount of i unit t, unit: ten thousand m³/a；

RDP_i,tRepresent reuse amount after the wastewater reuse approach factory sanitary sewage disposal of i unit t, unit: ten thousand m³/a；

CW_i,k,tRepresent the Tail water reuse concentration of the wastewater reuse approach factory k pollutant of i unit t, unit: mg/l；

WIP_i,tRepresent the wastewater reuse approach factory Industrial Wastewater Treatment amount of i unit t, unit: ten thousand m³/a；

RIP_i,tRepresent reuse amount after the wastewater reuse approach factory Industrial Wastewater Treatment of i unit t, unit: ten thousand m³/a；

PDP_i,k,tRepresent the generation concentration of k pollutant, unit: mg/l in i unit t sanitary sewage；

PIP_i,k,tRepresent the generation concentration of k pollutant, unit: mg/l in i unit t industrial wastewater；

L_kRepresent the clearance of wastewater reuse approach factory k pollutant, unit: %；

CS_kRepresent the discharge standard of sewage treatment plant's k pollutant, unit: mg/l；

QDP_i,tRepresent the sanitary sewage discharge capacity of i unit t, unit: ten thousand m³/a；

QIP_i,tRepresent the discharged volume of industrial waste water of i unit t, unit: ten thousand m³/a。

3. economy as claimed in claim 1 and the optimal screening method of point-source pollution abatement technology, its feature under environmental constraints It is:

In step 2, described quality of water environment model is zero dimension, one or more dimensions form.

4. economy as claimed in claim 1 and the optimal screening method of point-source pollution abatement technology, its feature under environmental constraints It is:

In step 2, pollutant enter river overall control index request WG_i,k,t=Min (M_i,k,t,CPT2_i,k,t), wherein:

WG_i,k,tCon trolling index requirement, unit: t/a are measured in the river that enters for i unit t k pollutant；

M_i,k,tThe river that enters for i unit k pollutant t is measured, unit: t/a；

CPT2_i,k,tFor the Water Functional Zone water environment capacity of i unit t k pollutant, unit: t/a.

5. economy as claimed in claim 1 and the optimal screening method of point-source pollution abatement technology, its feature under environmental constraints It is:

In step 4, point source pollutant write-down investment constraint KK_i,t=GDP_i,t·GK_i,t, wherein:

KK_i,tPollutant write-down investment for i unit t retrains, unit: Wan Yuan；

GDP_i,tFor the regional economy overall size of i unit t, the i.e. GDP of Water Functional Zone unit, unit: Wan Yuan；

GK_i,tScale of investment and the ratio of regional economy overall size, dimensionless for i unit t point-source pollution abatement.

6. economy as claimed in claim 1 and the optimal screening method of point-source pollution abatement technology, its feature under environmental constraints It is:

In step 4, the scale of investment of point-source pollution abatement is 0.5%-3.0% with the ratio of regional economy overall size.

7. economy as claimed in claim 1 and the optimal screening method of point-source pollution abatement technology, its feature under environmental constraints It is:

In step 5, the unit cost functions of point-source pollution abatement technology:

$C_{i,k,t}=\frac{Dr\×{\mathrm{Inv}}_{i,k,t}\×R{(1+R)}^{Lf}}{{(1+R)}^{Lf}-1}+{\mathrm{Ope}}_{i,k,t}$

Wherein:

C_i,k,tFor the totle drilling cost of i unit t abatement unit quantity k pollutant, unit: unit/year；

Inv_i,k,tFor the cost of investment of i unit t abatement unit quantity k pollutant, unit: unit/year；

Ope_i,k,tFor the operating cost of i unit t abatement unit quantity k pollutant, unit: unit/year；

Dr is cost of investment coefficient of depreciation, unit: 1/ year；

R is depreciation interest rate；

Lf is that point-source pollution abatement is arranged service life of technology, unit: year.

8. economy as claimed in claim 1 and the optimal screening method of point-source pollution abatement technology, its feature under environmental constraints It is:

In step 6, described object function is Max (W1_i,k,t+W2_i,k,t+W3_i,k,t), wherein:

W1_i,k,tThe k pollutant abatement amount realized for i unit t enterprise qualified discharge, unit: t/a；

W2_i,k,tThe k pollutant abatement amount realized for i unit t wastewater reuse approach factory, unit: t/a；

W3_i,k,tThe k pollutant abatement amount of realization, unit: t/a is cut down the most further for i unit t.

9. economy as claimed in claim 1 and the optimal screening method of point-source pollution abatement technology, its feature under environmental constraints It is:

In step 6, described constraints includes:

C1_i,k,t×W1_i,k,t+C2_i,k,t×W2_i,k,t+C3_i,k,t×W3_i,k,t≤KK_i,t

$W{1}_{i,k,t}+W{2}_{i,k,t}+W{3}_{i,k,t}\≤{\mathrm{PG}}_{i,k,t}-\frac{{\mathrm{WG}}_{i,k,t}}{K_{i,k,t}}$

W1_{I, k, t}≤LW1_{I, k, t}

W2_{I, k, t}≤LW2_{I, k, t}

W1_{I, k, t}≥0

W2_{I, k, t}≥0

W3_{I, k, t}≥0

Wherein:

C1_i,k,tFor the totle drilling cost of the abatement unit quantity k pollutant that i unit t wastewater reuse approach mode realizes, unit: ten thousand Unit/t；

C2_i,k,tFor the totle drilling cost of the abatement unit quantity k pollutant that i unit t polluter dispersion up to standard administration way realizes, list Position: ten thousand yuan/t；

C3_i,k,tFor the totle drilling cost of the abatement unit quantity k pollutant that i unit t realizes the most further, unit: ten thousand Unit/t；

W1_i,k,tThe k pollutant abatement amount realized for i unit t enterprise qualified discharge, unit: t/a；

W2_i,k,tThe k pollutant abatement amount realized for i unit t wastewater reuse approach factory, unit: t/a；

W3_i,k,tThe k pollutant abatement amount of realization, unit: t/a is cut down the most further for i unit t；

KK_i,tPollutant write-down investment for i unit t retrains, unit: Wan Yuan；

PG_i,k,tFor the generation amount of i unit t k pollutant, unit: t/a；

WG_i,k,tCon trolling index requirement, unit: t/a are measured in the river that enters for i unit t k pollutant；

K_i,k,tRiver pollutant sources for i unit t k pollutant；

LW1_i,k,tFor the i unit t enterprise maximum attainable k pollutant abatement amount of qualified discharge, unit: t/a；

LW2_i,k,tFor the maximum attainable k pollutant abatement amount of i unit t wastewater reuse approach factory, unit: t/a.