CN105844073A - Source/sink landscape recognition method and apparatus of non-point source pollution - Google Patents
Source/sink landscape recognition method and apparatus of non-point source pollution Download PDFInfo
- Publication number
- CN105844073A CN105844073A CN201610036702.4A CN201610036702A CN105844073A CN 105844073 A CN105844073 A CN 105844073A CN 201610036702 A CN201610036702 A CN 201610036702A CN 105844073 A CN105844073 A CN 105844073A
- Authority
- CN
- China
- Prior art keywords
- landscape cell
- landscape
- cell
- cost
- resistance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a source/sink landscape recognition method and apparatus of non-point source pollution. The method includes steps of determining resistance coefficients and dynamic coefficients of a current landscape unit and related landscape units thereof by employing resistance coefficients and dynamic coefficients of the current landscape unit and the related landscape units influencing a non-point source pollution process; calculating the resistance cost and the power cost of the current landscape unit by employing the resistance coefficients and the dynamic coefficients of the current landscape unit and the related landscape units thereof and the distances from the current landscape unit and the related landscape units thereof to sub-basin water outlets; and recognizing source/sink landscapes of the non-point source pollution according to the resistance cost and the power cost of the current landscape unit. According to the method and apparatus, the dynamic process of the formation of non-point source pollution from the landscape units to the sub-basin water outlets can be recognized and monitored, and an important significance is attached to further deep revealing of the coupling relationship between large-scale landscape patterns and the ecological process.
Description
Technical field
The present invention relates to source/remittance view recognition methods and the device of a kind of pollution of area source.
Background technology
Pollution of area source generating process is by local weather, the hydrology, topography and geomorphology, Land_use change, soil and plants
The capped effect waiting factor of influence, makes pollution of area source there is obvious spatial diversity on region, therefore, it is difficult to
Monitoring and the impact on water environment of the prevention and control pollution of area source.
Therefore, how from preferably, more fully in angle recognition and monitored area pollution of area source develops,
There is certain realistic meaning.
Summary of the invention
It is an object of the invention to provide source/remittance view recognition methods and the device of a kind of pollution of area source, can be more preferably
Ground solves identification and the problem of monitoring of pollution of area source in region.
According to an aspect of the invention, it is provided the source of a kind of pollution of area source/remittance view recognition methods, including:
Utilize and affect the current landscape cell of pollution of area source process and the drag parameter of relevant landscape cell thereof and move
Force parameter, determines described current landscape cell and the resistance coefficient of relevant landscape cell thereof and the coefficient of impact；
Utilize described current landscape cell and the resistance coefficient of relevant landscape cell thereof and the coefficient of impact, described currently
Landscape cell and relevant landscape cell thereof, to the distance of sub basin outlet, calculate the resistance of described current landscape cell
Power cost and power cost；
Resistance cost according to described current landscape cell and power cost, enter the source/remittance view of pollution of area source
Row identifies.
Preferably, the resistance system of described current landscape cell and relevant landscape cell thereof is determined by following steps
Number:
Obtain affect the current landscape cell of pollution of area source process and relevant landscape cell thereof include vegetative coverage
The drag parameter of rate, surface roughness and Manning roughness coefficient；
To any one landscape cell in described current landscape cell and relevant landscape cell thereof, by by described
After the vegetation coverage of landscape cell is added with surface roughness, with the Manning roughness coefficient phase of described landscape cell
Take advantage of, obtain the resistance coefficient of described landscape cell.
Preferably, described current landscape cell and the dynamical system of relevant landscape cell thereof are determined by following steps
Number:
Obtain affect the current landscape cell of pollution of area source process and relevant landscape cell thereof include average rainfall
The kinetic parameter of amount, slope factor, the erosion durability of soil and runoff curve number；
To any one landscape cell in described current landscape cell and relevant landscape cell thereof, by by described
After the average rainfall of landscape cell, slope factor and the erosion durability of soil are added, with the runoff of described landscape cell
Curve number is multiplied, and obtains the coefficient of impact of described landscape cell.
Preferably, the resistance cost of described current landscape cell is calculated by below equation:
Wherein, MCRR is the resistance cost of described current landscape cell, and n is to arrive subflow from current landscape cell
On a certain path of territory outlet j the quantity of landscape cell of process, m is sub basin outlet number, DijFor
The distance of sub basin outlet j, R is arrived from landscape cell iiResistance coefficient for landscape cell i.
Preferably, exist between described current landscape cell and relevant landscape cell thereof to sub basin outlet many
During paths, calculate the resistance cost corresponding to every paths respectively, and choose the resistance cost of minimum as institute
State the resistance cost of current landscape cell.
Preferably, the power cost of described current landscape cell is calculated by below equation:
Wherein, MCRE is the power cost of described current landscape cell, and n is to arrive subflow from current landscape cell
On a certain path of territory outlet j the quantity of landscape cell of process, m is sub basin outlet number, DijFor
The distance of sub basin outlet j, E is arrived from landscape cell iiThe coefficient of impact for landscape cell i.
Preferably, exist between described current landscape cell and relevant landscape cell thereof to sub basin outlet many
During paths, calculate the power cost corresponding to every paths respectively, and choose the power cost of minimum as institute
State the power cost of current landscape cell.
Preferably, by following steps identification source/remittance view:
Compare resistance cost and the power cost of described current landscape cell；
If described resistance cost is less than power cost, then described current landscape cell is identified as pollution of area source
" source " landscape cell；
If described resistance cost is more than power cost, then described current landscape cell is identified as pollution of area source
" converge " landscape cell；
If described resistance cost is equal to power cost, then described current landscape cell is identified as pollution of area source
" the balance of source and sink " landscape cell.
According to a further aspect in the invention, it is provided that the source of a kind of pollution of area source/remittance view identification device, including:
Coefficient determination module, affects the current landscape cell of pollution of area source process and relevant view list thereof for utilizing
The drag parameter of unit and kinetic parameter, determine described current landscape cell and the resistance coefficient of relevant landscape cell thereof
And the coefficient of impact；
Cost calculation module, for utilize described current landscape cell and relevant landscape cell thereof resistance coefficient and
The coefficient of impact, described current landscape cell and relevant landscape cell thereof, to the distance of sub basin outlet, calculate institute
State resistance cost and the power cost of current landscape cell；
Contamination Assessment module, for the resistance cost according to described current landscape cell and power cost, source, opposite
Source/remittance the view polluted is identified.
Preferably, the resistance cost of the more described current landscape cell of described Contamination Assessment module and power cost,
If described resistance cost is less than power cost, then described current landscape cell is identified as " source " of pollution of area source
Landscape cell；If described resistance cost is more than power cost, then described current landscape cell is identified as source, face dirty
" remittance " landscape cell of dye；If described resistance cost is equal to power cost, then described current landscape cell is known
Wei " the balance of source and sink " landscape cell of pollution of area source.
Compared with prior art, the beneficial effects of the present invention is:
The present invention is capable of identify that source remittance view, and monitors the landscape cell pollution of area source formation to sub basin outlet
Dynamic process, to the most deeply disclose large scale landscape pattern have important with the coupled relation of ecological process
Meaning.
Accompanying drawing explanation
Fig. 1 is the source/remittance view recognition methods block diagram of the pollution of area source that the embodiment of the present invention provides；
Fig. 2 is the source/remittance view identification device block diagram of the pollution of area source that the embodiment of the present invention provides；
Fig. 3 a is the cost calculation schematic diagram between the unit that the embodiment of the present invention provides；
Fig. 3 b is calculating with " source " Cummulative distance costs model as starting point of providing of the embodiment of the present invention
The source figure layer schematic diagram of principle；
Fig. 3 c is calculating with " source " Cummulative distance costs model as starting point of providing of the embodiment of the present invention
The cost figure layer schematic diagram of principle；
Fig. 3 d is calculating with " source " Cummulative distance costs model as starting point of providing of the embodiment of the present invention
The Cost Distance analysis result schematic diagram of principle；
Fig. 4 a be the embodiment of the present invention provide with " remittance " as starting point Cummulative distance costs model calculating
The source figure layer schematic diagram of principle；
Fig. 4 b be the embodiment of the present invention provide with " remittance " as starting point Cummulative distance costs model calculating
The cost figure layer schematic diagram of principle；
Fig. 4 c be the embodiment of the present invention provide with " remittance " as starting point Cummulative distance costs model calculating
The Cost Distance analysis result schematic diagram of principle；
Fig. 5 is " source/remittance " Plan Recognize flow chart affecting pollution of area source that the embodiment of the present invention provides.
Detailed description of the invention
Below in conjunction with accompanying drawing to a preferred embodiment of the present invention will be described in detail, it will be appreciated that described below
Preferred embodiment be merely to illustrate and explain the present invention, be not intended to limit the present invention.
Fig. 1 is the source/remittance view recognition methods block diagram of the pollution of area source that the embodiment of the present invention provides, such as Fig. 1 institute
Showing, step includes:
Step S101: utilize current landscape cell and the resistance of relevant landscape cell thereof affecting pollution of area source process
Force parameter and kinetic parameter, determine described current landscape cell and the resistance coefficient of relevant landscape cell thereof and power
Coefficient.
Drag parameter includes: current landscape cell and the vegetation coverage of relevant landscape cell, surface roughness
And Manning roughness coefficient.
Kinetic parameter includes: current landscape cell and the average rainfall of relevant landscape cell thereof, slope factor,
The erosion durability of soil and runoff curve number.
The resistance coefficient of each landscape cell is determined by following steps: the vegetation first obtaining this landscape cell is covered
Lid rate, surface roughness and Manning roughness coefficient, secondly by obtained vegetation coverage and surface roughness phase
Add, be then multiplied with Manning roughness coefficient, obtain the resistance coefficient of this landscape cell.
The coefficient of impact of each landscape cell is determined by following steps: obtain the average rainfall of this landscape cell
Amount, slope factor, the erosion durability of soil and runoff curve number, secondly by obtained average rainfall, the gradient because of
Son and the erosion durability of soil are added, and are then multiplied with runoff curve number, obtain the coefficient of impact of this landscape cell.
Described relevant landscape cell includes from current landscape cell to a certain path of corresponding sub basin outlet
The landscape cell of process.
Step S102: utilize described current landscape cell and the resistance coefficient of relevant landscape cell thereof and dynamical system
Current landscape cell several, described and relevant landscape cell thereof to the distance of sub basin outlet, calculate described currently
The resistance cost of landscape cell and power cost.
1, the resistance cost of described current landscape cell is calculated by below equation:
Wherein, MCRR is the resistance cost of described current landscape cell, and n is to arrive subflow from current landscape cell
On a certain path of territory outlet j the quantity of landscape cell of process, m is sub basin outlet number, DijFor
The distance of sub basin outlet j, R is arrived from landscape cell iiResistance coefficient for landscape cell i.
It should be noted that between described current landscape cell and relevant landscape cell thereof to sub basin outlet
When there is mulitpath, calculate the resistance cost corresponding to every paths respectively, and choose the resistance cost of minimum
Resistance cost as described current landscape cell.
2, the power cost of described current landscape cell is calculated by below equation:
Wherein, MCRE is the power cost of described current landscape cell, and n is to arrive subflow from current landscape cell
On a certain path of territory outlet j the quantity of landscape cell of process, m is sub basin outlet number, DijFor
The distance of sub basin outlet j, E is arrived from landscape cell iiThe coefficient of impact for landscape cell i.
It should be noted that between described current landscape cell and relevant landscape cell thereof to sub basin outlet
When there is mulitpath, calculate the power cost corresponding to every paths respectively, and choose the power cost of minimum
Power cost as described current landscape cell.
Step S103: according to resistance cost and the power cost of described current landscape cell, to pollution of area source
Source/remittance view is identified.
Specifically include: compare resistance cost and the power cost of described current landscape cell, if described resistance cost
Less than power cost, then described current landscape cell is identified as " source " landscape cell of pollution of area source；If institute
State resistance cost and be more than power cost, then described current landscape cell is identified as " remittance " view of pollution of area source
Unit；If described resistance cost is equal to power cost, then described current landscape cell is identified as pollution of area source
" the balance of source and sink " landscape cell.
The present invention is capable of identify that " source/remittance " landscape pattern of large scale pollution of area source process, has following a few side
The technical meaning in face:
1, the present invention has taken into full account the heterogeneity of background of regional environment, dirty evaluating source, view resistance/power opposite
The potential impact of dye process rises to true impact between the two, can compare in different areas and answer
With, greatly extend the range of application of source remittance landscape theories.
The present invention can be applied to the identification of " source/remittance " landscape pattern of Large-scale areas pollution of area source process,
It is no longer limited to " source/remittance " general layout formed from the mosaic mode that landscape types is intrinsic to consider pollution of area source,
But from affect view stream pollution of area source forming process must be obstructed/restriction of kinetic factor identifies " source/remittance "
General layout, has therefore taken into full account the heterogeneity of background of regional environment.All of sub basin water in being determined by region
Mouth locus, by the potential impact of source, view resistance/power opposite pollution course, rises to and arrives for landscape cell
The pollution of area source of sub basin outlet forms the impact of dynamic process, has more real realistic meaning, and energy
It is applied to area, different basin.Therefore, the present invention greatly extends the range of application of source remittance landscape theories, more
In depth the identification of " source/remittance " landscape pattern is applied in concrete environmental background.
2, affect " source/remittance " landscape pattern of pollution of area source process by identifying, for identify rapidly and accurately with
Delimit pollution of area source key area and effective means be provided.
The present invention is capable of identify that " source/remittance " landscape pattern affecting this ecological process of pollution of area source, wherein " source "
General layout refers to promote the landscape types of pollution of area source process development, is the starting point that outwards spreads of pollution of area source material
And base, under such general layout state, pollution of area source is the most incident, the emphasis ground occurred for pollution of area source
District；General layout of " converging " refers to the place that a pollution of area source process disappears, and is to stop to delay pollution of area source process to send out
The landscape types of exhibition.Therefore, by identifying " source/remittance " landscape pattern affecting pollution of area source process, it is possible to
Judge the general layout situation in " source " and " remittance " that affect pollution of area source on big regional scale, further determine that
The key area that in this region, pollution of area source occurs, it is possible to provide reference into planning and prevention and control pollution of area source.
3, by identifying " source/remittance " landscape pattern affecting pollution of area source process, it is possible to formed for pollution of area source
Danger important directive significance is provided.
Due to " source ", " remittance " effect existing during pollution of area source, therefore " source/remittance " landscape pattern
Whether reasonable layout influences whether the risk degree that pollution of area source is formed.If " source " " remittance " scape in region
See in spatial distribution, reach poised state, form rational Spatial Distribution Pattern, it will produce less non-dots
Source pollutants；Whereas if " source " " remittance " landscape pattern irrational distribution, and there is more " source " scape
Sight is distributed near exit, basin, it will have more area source pollutants to produce, and the danger of pollution of area source will
Significantly increase.Therefore by identifying " source/remittance " general layout affecting pollution of area source process, it is possible to for pollution of area source
The danger formed provides important directive significance.
Fig. 2 is the source/remittance view identification device block diagram of the pollution of area source that the embodiment of the present invention provides, such as Fig. 2 institute
Show, including coefficient determination module, cost calculation module and Contamination Assessment module.
Coefficient determination module affects the current landscape cell of pollution of area source process and relevant view list thereof for utilizing
The drag parameter of unit and kinetic parameter, determine described current landscape cell and the resistance coefficient of relevant landscape cell thereof
And the coefficient of impact.
Cost calculation module for utilize described current landscape cell and relevant landscape cell thereof resistance coefficient and
The coefficient of impact, described current landscape cell and relevant landscape cell thereof, to the distance of sub basin outlet, calculate institute
State resistance cost and the power cost of current landscape cell.
Contamination Assessment module is for the resistance cost according to described current landscape cell and power cost, and source, opposite is dirty
Source/remittance the view of dye is identified.Specifically, the resistance of the more described current landscape cell of Contamination Assessment module
Power cost and power cost, if described resistance cost is less than power cost, then by described current landscape cell identification
" source " landscape cell for pollution of area source；If described resistance cost is more than power cost, then work as prospect by described
See unit and be identified as " remittance " landscape cell of pollution of area source；If described resistance cost is equal to power cost, then will
Described current landscape cell is identified as " the balance of source and sink " landscape cell of pollution of area source..
The present invention, from the view resistance/power affecting large scale pollution of area source forming process, merges each view and wants
Element (such as Land_use change, soil, the hydrology, topography and geomorphology, vegetative coverage etc.) builds affects pollution of area source shape
View resistance/the coefficient of impact become, simultaneously, it is considered to landscape cell and the distance factor of sub basin outlet, finally
Obtain based on Large-scale areas pollution of area source process " source/remittance " landscape pattern and identify model.
1, below using the outlet of reservoir area of Three Gorges sub basin as the Rendezvous Point of area source pollutants, and it is applied to
As a example by little accumulation resistance model, illustrate.
Minimum principle designed by accumulation resistance model is to calculate from source through spent by the view of different resistances
Expense or overcome resistance institute work." source " therein refers in general layout with process study, those energy
Promote the landscape types of ecological process development, be things or the event starting point to external diffusion and base, there is inside
Homogeneity and the ability to surrounding extension, therefore, the minimum accumulation resistance model of application be necessary for determining for
" source " of specific ecological process.But from the point of view of this ecological process of pollution of area source, affect pollution of area source
Whether " source " have the feature that is hidden, it is difficult to directly go to identify and judge from specific landscape types
It is " source ", especially cannot be determined by experiment the method for quantization in the range of large-scale dimension to determine source, face
" source " of pollution course, so it will be unpractical for directly applying traditional minimum accumulation resistance model.Thus
" source/remittance " landscape pattern calculated in pollution of area source process to the cumulative cost of different landscape unit from source identifies
Existing defects in technology.The present invention is reservoir area of Three Gorges sub basin in view of targeted spatial domain, and basin conduct
Gathering ground, its outlet is the place that pollution of area source material finally collects through different landscape unit, i.e. source, face is dirty
Dye process is the final locus terminated in basin.Therefore the present invention is directed to the former of minimum accumulation resistance model
Reason, finally collects part using the outlet of sub basin as area source pollutants matter in basin, calculates landscape cell and arrives
The Cummulative distance costs of sub basin outlet.So, one side avoids and cannot determine pollution of area source process
" source " and the problem of None-identified " source/remittance " landscape pattern, on the other hand also further expanded consuming
The application of distance costs distance model, i.e. inversely applies it according to the feature of pollution of area source process.
2, create with the kinetic theory of view stream and affect the view resistance/coefficient of impact of pollution of area source process, and
As the view resistance and the power face that affect pollution of area source process in minimum accumulation resistance model.
The most crucial process affecting pollution of area source formation includes that rainfall runoff, Soil erosion and contaminant transportation three are big
Aspect, this three broad aspect is directed to the effect of view stream, is mainly reflected in rainwash or interflow subsurface drainage, soil
Earth corrodes stream and pollutant mass flow, and pollutant mass flow therein is accompanied by the generation of rainfall runoff and Soil erosion
And advance, therefore pollution of area source material moves to " remittance " this process inherently by various views from " source "
The dynamic action of stream and drag effect.Wherein, it is provided that the main source of view mobilization force include gravity, rainfall,
Runoff and Soil structure etc., provide the main source of view flow resistance then to include vegetation cover, earth's surface
Degree of roughness etc..Therefore, the present invention chooses affects the scenic powered index of pollution of area source process (i.e. kinetic parameter)
Including the slope factor in average rainfall, Soil erosion and the erosion durability of soil, and runoff yield ability can be reflected
Runoff curve number, the leading indicator (i.e. drag parameter) choosing the view resistance affecting pollution of area source process includes
Vegetation cover, roughness of ground surface and Manning roughness coefficient, the most as shown in table 1.
The present invention passes through resistance/coefficient of impact, sets up the view resistance/power face affecting pollution of area source process.Specifically
Say, give resistance coefficient and the coefficient of impact to landscape cell, and by the view on these landscape cell Special compositions
Resistance area and power face, and the view resistance/coefficient of impact determines the controllability of pollution of area source forming process, if scape
See dynamic action and be more than drag effect, show as pollution of area source process and be easier to occur, be the most then susceptible to.
The mode that embodies of the view resistance/coefficient of impact is to use linear mode generally to change, specifically such as formula (1) and
Shown in formula (2).
R=(VFC+M) × N (1)
E=(P+S+K) × CN (2)
The index representated by each parameter in above-mentioned formula (1) and formula (2) refers to table 1.
View resistance/dynamic index the parameter of table 1 modelling
3, on the basis of minimum accumulation resistance model (Minimum Cumulative Resistance, MCR)
On construct the view resistance/power cost model affecting pollution of area source process.
The essence of MCR is the Integrative expression of Cost Distance (cost distance), is that the one to reality is abstract
Express, refer to bring " cost " difference due to the change of landscape, Cost Distance be determined by material,
Energy calculates at the cost coefficient of different surfaces, shown in concrete computing formula such as formula (3).
In formula, MCR is minimum accumulation Resistance Value, RiRepresent the resistance coefficient that certain species is moved by landscape cell i,
∑ represents and passes through the distance of all unit and the accumulation of resistance between landscape cell i and source j, and min represents and commented
The speckle of valency takes accumulation resistance minima for different sources；Minimum accumulation Resistance Value reflects the latent of ecological process
In probability and trend, can determine whether this unit and source unit by the size of the minimum accumulation Resistance Value of landscape cell
Connectedness and similarity.
The key that Cost Distance calculates is accumulative, and Fig. 3 a to Fig. 3 d provides one with " source " tiring out as starting point
Meter expends distance costs model Computing Principle schematic diagram, and Fig. 4 a to Fig. 4 c provides one with " remittance " as starting point
Cummulative distance costs model Computing Principle schematic diagram, as shown in Fig. 3 a to Fig. 4 c.Each unit passes through
The calculating of cost, had both considered the distance between grid cell, it is also considered that the current cost of each unit, cumulative cost
Being that the current cost from the cumulative cost of a upper unit and this unit is added, computing formula is: accum_cost
=al+ ((cost_a+cost_b)/2) × D, wherein, accum_cost is the tired of a certain unit b
The current cost of meter；Al is the accumulative current cost of a upper adjacent cells a；Cost_a is the cost of unit a；
Cost_b is the cost of unit b；D is the distance between adjacent cells a, b, by the central point of two unit
Calculate.If a length of d of element sides, the distance between adjacent cells is exactly d the most up and down；Diagonal
Distance between adjacent cells is
Due to above-mentioned Cummulative Cost Distance model only for be the view of different resistance, only take into account " source "
The view resistance that run into of expansion, but for this ecological process of pollution of area source, owing to there is the work of view stream
With, both there is the effect of view resistance, there is also the effect of scenic powered, therefore at resistance area and power simultaneously
In the case of face is simultaneous, need according to the above view resistance/coefficient of impact, the meter of minimum cumulative resistance
Calculate formula and principle sets up corresponding view resistance cost model and scenic powered cost model, the most simultaneously
Shown in formula (4) and formula (5).
In formula: MCRR and MCRE is a certain landscape cell Study of Platycladus orientalis-Robinia pseudoacacia resistance and the minimum cumulative cost of power；N is
From a certain path that current landscape cell arrives sub basin outlet j the quantity of landscape cell of process；m
For sub basin outlet number；DijDistance for i-th landscape cell to corresponding basin outlet j；VFCi、MiWith
NiIt is respectively the vegetation coverage of i-th landscape cell, roughness of ground surface and Manning roughness coefficient；Pi、Si、KiWith
CNiIt is respectively the average precipitation of i-th landscape cell, slope factor, soil erodibility and runoff curve number.
4, application ArcGIS software raster symbol-base device identification judges to affect " source/remittance " scape of pollution of area source process
See general layout.
By view resistance cost MCRR and scenic powered cost MCRE, by ArcGIS raster symbol-base device
(MCRR>MCRE, 1, Con (MCRR<MCRE, 2,0), in formula, MCRR, MCRE are respectively for conditional function: Con
Representing view resistance/power cost grid map, 1 and 2 represent " source " and " remittance " respectively, and circulation identifies impact
" source/remittance " space cell that pollution of area source is formed.If MCRR > MCRE, i.e. this grid space unit is to right
Answering the view resistance cost of sub basin outlet more than power cost, this unit shows as " remittance " of pollution of area source,
If MCRE > MCRR, i.e. this grid space unit are less than dynamic to the view resistance cost of corresponding sub basin outlet
Power cost, shows as " source "；If MCRR=MCRE, then this landscape cell affects the scape that pollution of area source is formed
Seeing resistance cost and be equal to power cost, " source " and " remittance " general layout reaches balance.
Fig. 5 is " source/remittance " Plan Recognize flow chart affecting pollution of area source that the embodiment of the present invention provides, as
Shown in Fig. 5, step includes:
The first step: obtain the vegetation coverage VFC of each relevant landscape cell, surface roughness M, graceful peaceful rough system
Number N, average rainfall P, slope factor S, soil erodibility K, runoff curve number CN.
Second step: utilize vegetation coverage VFC that formula (1) and the first step obtain, surface roughness M, graceful
Peaceful rough coefficient N, calculates the resistance coefficient R of relevant landscape cell respectively；Formula (2) and the first step is utilized to obtain
Average rainfall P, slope factor S, soil erodibility K, runoff curve number CN, calculate relevant scape respectively
See the coefficient of impact of unit.
3rd step: utilize formula (4) and relevant landscape cell to the distance of each sub basin outlet, relevant scape
See the resistance coefficient R of unit, calculate the view resistance cost MCRR of each landscape cell；Utilize formula (5)
With the distance of relevant landscape cell to each sub basin outlet, the coefficient of impact E of relevant landscape cell, calculate every
The scenic powered cost MCRE of individual landscape cell.
4th step: for certain landscape cell, compares its view resistance cost MCRR and scenic powered cost
MCRE, if MCRR < MCRE, is then defined as " source " landscape cell by this landscape cell；If MCRR > MCRE,
Then this landscape cell is defined as " remittance " landscape cell.
The present invention can be that pollution of area source hazard assessment provides technological means and important references, with the face of reservoir area of Three Gorges
As a example by source is polluted, a feature of reservoir area of Three Gorges maximum is that water environment pollution faces the situation of sternness, especially
The pollution of area source impact on the water body such as river, lake, its critical process is the loss of nutrient.The present invention is permissible
Using reservoir area of Three Gorges as a Large-scale areas entirety, in source remittance landscape theories, view flow resistance is theoretical and tires out
Carry out affecting " source remittance " Plan Recognize of pollution of area source process on the basis of meter resistance Cost Theory.Due to " source "
Being pollution of area source process and the starting point of expansion, the source of pollution of area source material nutrient, " remittance " has and receives undertaking
The function that pollution of area source material migrates, is therefore producing on pollution of area source " source " landscape pattern than " remittance " scape
See general layout and have more danger, if controlled at spatial and temporal scales by spatially combining different types of view
On the balance of nutrient loss, the danger of pollution of area source will be reduced, therefore the present invention is dangerous at pollution of area source
Property evaluation is provided that important technological means and reference value.
The present invention can be the basin for the purpose of prevention and control pollution of area source or Ecological planning of regional landscape offer technology hands
Section.Landscape Ecological Planning for the purpose of prevention and control pollution of area source is with regional landscape ecological system global optimization as base
This target, by research landscape pattern and pollution of area source process and mankind's activity and the interaction of view, builds
The space structure of vertical regional landscape ecological system optimization utilization and pattern, make the views such as gallery, speckle, substrate want
The quantity of element and spatial distribution thereof are reasonable, make flow of information, material stream lead to energy smoothness.It can thus be seen that
Landscape Ecological Planning to be implemented is necessary for " source/remittance " effect wanting landscape pattern in the pollution course of source, opposite to be played,
And the empty heterogeneity at " source/remittance " general layout place, only recognize these effect and relations, it becomes possible to clear
The key area of pollution of area source process, deathtrap and safety zone is affected in the range of Large-scale areas, thus
Go to reasonable disposition " source/remittance " landscape pattern with having purpose and planning, make pollution of area source material in heterogeneous landscape
Redistribute, thus reach the purpose of Control the Plan Pollution Sources.Therefore the present invention can be by identifying Large-scale areas
" source/remittance " view of pollution of area source process provides technological means for Ecological planning of regional landscape.
In sum, the present invention is capable of identify that " source/remittance " general layout affecting pollution of area source process, concrete
The hazard assessment of pollution of area source, basin or Ecological planning of regional landscape, prevention and control it are mainly reflected on practical use
The aspects such as pollution of area source.
Although above the present invention being described in detail, but the invention is not restricted to this, the art technology
Personnel can carry out various amendment according to the principle of the present invention.Therefore, all amendments made according to the principle of the invention,
All should be understood to fall into protection scope of the present invention.
Claims (10)
1. the source of a pollution of area source/remittance view recognition methods, it is characterised in that including:
Utilize and affect the current landscape cell of pollution of area source process and the drag parameter of relevant landscape cell thereof and move
Force parameter, determines described current landscape cell and the resistance coefficient of relevant landscape cell thereof and the coefficient of impact；
Utilize described current landscape cell and the resistance coefficient of relevant landscape cell thereof and the coefficient of impact, described currently
Landscape cell and relevant landscape cell thereof, to the distance of sub basin outlet, calculate the resistance of described current landscape cell
Power cost and power cost；
Resistance cost according to described current landscape cell and power cost, enter the source/remittance view of pollution of area source
Row identifies.
Method the most according to claim 1, it is characterised in that by following steps determine described currently
Landscape cell and the resistance coefficient of relevant landscape cell thereof:
Obtain affect the current landscape cell of pollution of area source process and relevant landscape cell thereof include vegetative coverage
The drag parameter of rate, surface roughness and Manning roughness coefficient；
To any one landscape cell in described current landscape cell and relevant landscape cell thereof, by by described
After the vegetation coverage of landscape cell is added with surface roughness, with the Manning roughness coefficient phase of described landscape cell
Take advantage of, obtain the resistance coefficient of described landscape cell.
Method the most according to claim 1, it is characterised in that by following steps determine described currently
Landscape cell and the coefficient of impact of relevant landscape cell thereof:
Obtain affect the current landscape cell of pollution of area source process and relevant landscape cell thereof include average rainfall
The kinetic parameter of amount, slope factor, the erosion durability of soil and runoff curve number；
To any one landscape cell in described current landscape cell and relevant landscape cell thereof, by by described
After the average rainfall of landscape cell, slope factor and the erosion durability of soil are added, with the runoff of described landscape cell
Curve number is multiplied, and obtains the coefficient of impact of described landscape cell.
Method the most according to claim 1, it is characterised in that by below equation calculate described currently
The resistance cost of landscape cell:
Wherein, MCRR is the resistance cost of described current landscape cell, and n is to arrive subflow from current landscape cell
On a certain path of territory outlet j the quantity of landscape cell of process, m is sub basin outlet number, DijFor
The distance of sub basin outlet j, R is arrived from landscape cell iiResistance coefficient for landscape cell i.
Method the most according to claim 4, it is characterised in that when described current landscape cell and phase thereof
Close landscape cell to when there is mulitpath between sub basin outlet, calculate the resistance corresponding to every paths respectively
Power cost, and choose the resistance cost resistance cost as described current landscape cell of minimum.
Method the most according to claim 1, it is characterised in that by below equation calculate described currently
The power cost of landscape cell:
Wherein, MCRE is the power cost of described current landscape cell, and n is to arrive subflow from current landscape cell
On a certain path of territory outlet j the quantity of landscape cell of process, m is sub basin outlet number, DijFor
The distance of sub basin outlet j, E is arrived from landscape cell iiThe coefficient of impact for landscape cell i.
Method the most according to claim 6, it is characterised in that when described current landscape cell and phase thereof
Close landscape cell to when there is mulitpath between sub basin outlet, calculate moving corresponding to every paths respectively
Power cost, and choose the power cost power cost as described current landscape cell of minimum.
Method the most according to claim 1, it is characterised in that by following steps identification source/remittance view:
Compare resistance cost and the power cost of described current landscape cell；
If described resistance cost is less than power cost, then described current landscape cell is identified as pollution of area source
" source " landscape cell；
If described resistance cost is more than power cost, then described current landscape cell is identified as pollution of area source
" converge " landscape cell；
If described resistance cost is equal to power cost, then described current landscape cell is identified as pollution of area source
" the balance of source and sink " landscape cell.
9. the source of a pollution of area source/remittance view identification device, it is characterised in that including:
Coefficient determination module, affects the current landscape cell of pollution of area source process and relevant view list thereof for utilizing
The drag parameter of unit and kinetic parameter, determine described current landscape cell and the resistance coefficient of relevant landscape cell thereof
And the coefficient of impact；
Cost calculation module, for utilize described current landscape cell and relevant landscape cell thereof resistance coefficient and
The coefficient of impact, described current landscape cell and relevant landscape cell thereof, to the distance of sub basin outlet, calculate institute
State resistance cost and the power cost of current landscape cell；
Contamination Assessment module, for the resistance cost according to described current landscape cell and power cost, source, opposite
Source/remittance the view polluted is identified.
Device the most according to claim 9, it is characterised in that described Contamination Assessment module is relatively more described
The resistance cost of current landscape cell and power cost, if described resistance cost is less than power cost, then by described
Current landscape cell is identified as " source " landscape cell of pollution of area source；If described resistance cost is more than power cost,
Then described current landscape cell is identified as " remittance " landscape cell of pollution of area source；If described resistance cost is equal to
Power cost, then be identified as " the balance of source and sink " landscape cell of pollution of area source by described current landscape cell.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610036702.4A CN105844073B (en) | 2016-01-20 | 2016-01-20 | A kind of source of pollution of area source/remittance landscape recognition methods and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610036702.4A CN105844073B (en) | 2016-01-20 | 2016-01-20 | A kind of source of pollution of area source/remittance landscape recognition methods and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105844073A true CN105844073A (en) | 2016-08-10 |
| CN105844073B CN105844073B (en) | 2019-03-26 |
Family
ID=56580554
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610036702.4A Active CN105844073B (en) | 2016-01-20 | 2016-01-20 | A kind of source of pollution of area source/remittance landscape recognition methods and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105844073B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107563803A (en) * | 2017-08-24 | 2018-01-09 | 北京工商大学 | A kind of market area partition method based on cost grid |
| CN111898896A (en) * | 2020-07-24 | 2020-11-06 | 中国科学院城市环境研究所 | Watershed non-point source pollution loss intensity evaluation method considering soil attributes |
| CN113553789A (en) * | 2021-07-30 | 2021-10-26 | 中国科学院城市环境研究所 | Resistance distance calculation method and system for ecological effect of green land landscape on building energy consumption |
| CN114662930A (en) * | 2022-03-24 | 2022-06-24 | 中国科学院地理科学与资源研究所 | Agricultural non-point source pollution risk identification method and electronic equipment |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103077307A (en) * | 2012-12-31 | 2013-05-01 | 清华大学 | Drought and water-starved river ecological restoration-based hydrodynamic modeling method |
| CN104091052A (en) * | 2014-06-25 | 2014-10-08 | 安徽省环境科学研究院 | Non-point source pollution spatial distributed model |
| CN104299161A (en) * | 2013-07-19 | 2015-01-21 | 重庆师范大学 | Method and device for obtaining graphic data of county-scale abandoned land |
| CN105000748A (en) * | 2015-04-30 | 2015-10-28 | 天津大学 | Typical farmland non-point source pollutant 'point-line-face' comprehensive intercepting and control system |
-
2016
- 2016-01-20 CN CN201610036702.4A patent/CN105844073B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103077307A (en) * | 2012-12-31 | 2013-05-01 | 清华大学 | Drought and water-starved river ecological restoration-based hydrodynamic modeling method |
| CN104299161A (en) * | 2013-07-19 | 2015-01-21 | 重庆师范大学 | Method and device for obtaining graphic data of county-scale abandoned land |
| CN104091052A (en) * | 2014-06-25 | 2014-10-08 | 安徽省环境科学研究院 | Non-point source pollution spatial distributed model |
| CN105000748A (en) * | 2015-04-30 | 2015-10-28 | 天津大学 | Typical farmland non-point source pollutant 'point-line-face' comprehensive intercepting and control system |
Non-Patent Citations (6)
| Title |
|---|
| JONATHAN A.GREENBERG等: "Low-cost distance analysis for spatial interpolation", 《COMPUTER&GEOSCIENCES》 * |
| 叶玉瑶等: "生态阻力面模型构建及其在城市扩展模拟中的应用", 《地理学报》 * |
| 张军: "GIS支持下淮河流域淮南市中段农业非点源污染研究", 《中国优秀硕士学位论文全文数据库 工程科技I辑》 * |
| 潘竟虎等: "基于空间主成分和最小累积阻力模型的内陆河景观生态安全评价与格局优化-以张掖市甘州区为例", 《应用生态学报》 * |
| 郝芳华等: "大尺度区域非点源污染负荷计算方法", 《环境科学学报》 * |
| 陈利顶等: "基于"源-汇"生态过程的景观格局识别方法-景观空间负荷对比指数", 《生态学报》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107563803A (en) * | 2017-08-24 | 2018-01-09 | 北京工商大学 | A kind of market area partition method based on cost grid |
| CN107563803B (en) * | 2017-08-24 | 2020-12-08 | 北京工商大学 | Quotient circle dividing method based on distance consumption grid |
| CN111898896A (en) * | 2020-07-24 | 2020-11-06 | 中国科学院城市环境研究所 | Watershed non-point source pollution loss intensity evaluation method considering soil attributes |
| CN113553789A (en) * | 2021-07-30 | 2021-10-26 | 中国科学院城市环境研究所 | Resistance distance calculation method and system for ecological effect of green land landscape on building energy consumption |
| CN114662930A (en) * | 2022-03-24 | 2022-06-24 | 中国科学院地理科学与资源研究所 | Agricultural non-point source pollution risk identification method and electronic equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105844073B (en) | 2019-03-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kurczyn et al. | Mesoscale eddies in the northeastern Pacific tropical‐subtropical transition zone: Statistical characterization from satellite altimetry | |
| Hang et al. | Wind conditions in idealized building clusters: macroscopic simulations using a porous turbulence model | |
| Rui et al. | GIS-based emergency response system for sudden water pollution accidents | |
| CN104898183A (en) | Modeling evaluation method for urban heavy rain inundation | |
| Erpicum et al. | Detailed inundation modelling using high resolution DEMs | |
| CN105844073A (en) | Source/sink landscape recognition method and apparatus of non-point source pollution | |
| CN104933899A (en) | Ship seaworthiness early-warning system and method | |
| Wu et al. | Integrated hydrologic and hydrodynamic modeling to assess water exchange in a data-scarce reservoir | |
| El Kasmi et al. | Turbulence modeling of atmospheric boundary layer flow over complex terrain: a comparison of models at wind tunnel and full scale | |
| CN105761302A (en) | Three-dimensional digital navigation channel system and application thereof | |
| Li et al. | A visual analytics approach for flood risk analysis and decision-making in cultural heritage | |
| Dinçer et al. | Enhancing wind turbine site selection through a novel wake penalty criterion | |
| Ju et al. | Planning for the change: mapping sea level rise and storm inundation in Sherman Island using 3Di hydrodynamic model and LiDAR | |
| Zhao et al. | Tsunami waveforms and runup of undular bores in coastal waters | |
| CN205068800U (en) | Boats and ships seaworthiness early warning system | |
| Goeury et al. | Finding good solutions to telemac optimization problems with a metaheuristic | |
| Choi et al. | Estimates of tsunami risk zones on the coasts adjacent to the East (Japan) Sea based on the synthetic catalogue | |
| Devkota et al. | Response characteristics of the Perdido and Wolf Bay system to inflows and sea level rise. | |
| Dewals et al. | Semi-explicit modelling of watersheds with urban drainage systems | |
| Koh et al. | Tsunami-driven debris hazard assessment at a coastal community: Focusing on shipping container debris hazards at Honolulu Harbor, Hawaii | |
| Sitzenfrei | Stochastic generation of urban water systems for case study analysis | |
| Chen et al. | Wargame Mapping and Implementation for Emergency Evacuation of Residents in Urban Waterlogging Disaster | |
| Tian et al. | Influence of model geometrical fitting and turbulence parameterization on phytoplankton simulation in the Gulf of Maine | |
| KR20160061521A (en) | Decision method and system of debris blocking facilities site on river using 3-dimentional hydrodynamic models, and the system | |
| Dewals et al. | Integrated assessment of flood protection measures in the context of climate change: hydraulic modelling and economic approach |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |