CN109272246A - A kind of river clustering method based on algal bloom - Google Patents
A kind of river clustering method based on algal bloom Download PDFInfo
- Publication number
- CN109272246A CN109272246A CN201811165966.5A CN201811165966A CN109272246A CN 109272246 A CN109272246 A CN 109272246A CN 201811165966 A CN201811165966 A CN 201811165966A CN 109272246 A CN109272246 A CN 109272246A
- Authority
- CN
- China
- Prior art keywords
- index
- river
- water body
- distance
- water
- 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
- 238000000034 method Methods 0.000 title claims abstract description 66
- 239000005422 algal bloom Substances 0.000 title claims abstract description 43
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 117
- 235000016709 nutrition Nutrition 0.000 claims abstract description 11
- 230000035764 nutrition Effects 0.000 claims abstract description 11
- 239000000126 substance Substances 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 22
- 230000003287 optical effect Effects 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 9
- 230000000877 morphologic effect Effects 0.000 claims description 8
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 6
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 4
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 claims description 4
- YUWBVKYVJWNVLE-UHFFFAOYSA-N [N].[P] Chemical compound [N].[P] YUWBVKYVJWNVLE-UHFFFAOYSA-N 0.000 claims description 4
- UMVBXBACMIOFDO-UHFFFAOYSA-N [N].[Si] Chemical compound [N].[Si] UMVBXBACMIOFDO-UHFFFAOYSA-N 0.000 claims description 4
- HIVGXUNKSAJJDN-UHFFFAOYSA-N [Si].[P] Chemical compound [Si].[P] HIVGXUNKSAJJDN-UHFFFAOYSA-N 0.000 claims description 4
- 229930002875 chlorophyll Natural products 0.000 claims description 4
- 235000019804 chlorophyll Nutrition 0.000 claims description 4
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 claims description 4
- 230000033116 oxidation-reduction process Effects 0.000 claims description 4
- 229910052698 phosphorus Inorganic materials 0.000 claims description 4
- 239000011574 phosphorus Substances 0.000 claims description 4
- 229910052710 silicon Inorganic materials 0.000 claims description 4
- 239000010703 silicon Substances 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 3
- 238000005452 bending Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 3
- 230000007423 decrease Effects 0.000 claims description 3
- 238000011144 upstream manufacturing Methods 0.000 claims description 3
- 229910052757 nitrogen Inorganic materials 0.000 claims description 2
- 229910017435 S2 In Inorganic materials 0.000 claims 1
- 239000002352 surface water Substances 0.000 claims 1
- 230000002265 prevention Effects 0.000 abstract description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012851 eutrophication Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- BHEPBYXIRTUNPN-UHFFFAOYSA-N hydridophosphorus(.) (triplet) Chemical compound [PH] BHEPBYXIRTUNPN-UHFFFAOYSA-N 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 230000000050 nutritive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Educational Administration (AREA)
- Economics (AREA)
- Development Economics (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Primary Health Care (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of river clustering method based on algal bloom, comprising the following steps: establish river algal bloom index system;Determine system index weights;Parameter score;River cluster.The present invention comprehensively considers river algal bloom index, including river form, water body power, water body physical and chemical factor, water body optics, water nutrition index, establish river algal bloom index system, and each index weights are determined using subjective and objective combination weights method, finally river cluster is carried out using hierarchical clustering method, to which the river of different zones, different characteristic be sorted out, the universal law for understanding river algal bloom, provides a kind of effective technical management tool for wawter bloom risk profile and wawter bloom prevention and control.
Description
Technical field
The present invention relates to a kind of clustering methods, and in particular to a kind of river clustering method based on algal bloom.
Background technique
With social development, mankind's activity is frequent, and great variety occurs for river channel ecology environment, and different degrees of wawter bloom causes society
It can pay close attention to.Algal bloom is response of the aquatic ecosystem to water eutrophication, and algal bloom directly destroys ecology of water landscape,
The survival and development for inhabiting biology in water area ecological environment are endangered, or even threaten human health.It is quick-fried for research river wawter bloom
Hair mechanism, environment scholar carry out from river form, river water body hydrodynamic characterisitic, water temperature, nutritive salt characteristic, sociales etc.
Research finds algal bloom key factor are as follows: slow-flowing stream fluidised form;It is suitable for illumination, water temperature condition;Sufficient nitrogen and phosphorous nutrient.
Currently, for river algal bloom Mechanism Study mainly for single river, for different zones different characteristic
River, if there is cluster, if with similar wawter bloom risk and mechanism even prevention and control measure, lack research.
For the universal law for understanding the river algal bloom of different zones different characteristic, the present invention proposes a kind of quick-fried based on wawter bloom
The river New Clustering of hair.
Summary of the invention
To solve the deficiencies in the prior art, it is an object of that present invention to provide one kind comprehensively consider river form, water body power,
Water body physical and chemical factor, water body optics, water nutrition index, using the entitled mode of subjective and objective combination weights method, science, reasonably
The method for carrying out river cluster.
The technical solution adopted by the present invention are as follows:
A kind of river clustering method based on algal bloom, comprising the following steps:
S1, river algal bloom index system, including first class index layer and two-level index layer are established;
The first class index layer is made of several first class index, including river morphological index, water body dynamic index, water body reason
Change level of factor, water body optical index, water nutrition index;
The two-level index layer is made of several two-level index,
River morphological index includes length index, area index, wriggles and spend index,
Water body dynamic index includes water level index, velocity parameters, the vertical coefficient of stability index of water body,
Water body physical and chemical factor index includes temperature index, pH value index, conductivity indices, oxidation-reduction potential index, DO
Index,
Water body optical index includes transparency index, turbidity index, light attenuation coefficient index, euphotic zone depth index, mixes
Layer depth index is closed,
Water nutrition index includes total nitrogen (TN) index, total phosphorus (TP) index, dissolved silicon hydrochlorate (D-Si) index, nitrogen phosphorus
(TN/TP) than index, silicon nitrogen ratio (D-Si/TN) index, silicon phosphorus ratio (D-Si/TP) index, chlorophyll (Chla) index;
S2, using subjective and objective combination weights method, determine each two-level index weight in each river;
S3, each first class index score in each river is determined using weighted sum method based on each two-level index weight in each river;
S4, using hierarchical clustering method, carry out river cluster.
Subjective and objective combination weights method in above-mentioned steps S2 includes subjective weighting method, objective weighted model;
Determine system index weights, specifically:
Assuming that two-level index number is J:
wj=α aj+(1-α)bjJ=1,2 ..., J, 0≤α≤1 (1)
In formula (1), wjFor the comprehensive weight of j-th of index, ajAnd bjThe subjective weight of respectively j-th index and objective
Weight, α is preference coefficient, according to policymaker to the preference of different enabling legislations;
The subjective weighting method uses expert survey, and n experts understand river algal bloom index system, and anonymity is finger
It marks j and assigns power, ajIt is averaged;
The objective weighted model uses Information Entropy, and steps are as follows:
A1: nondimensionalization processing
Nondimensionalization processing is carried out to index system data using Maximum Approach, it is assumed that there is N river to participate in clustering,
Assuming that two-level index number is J, matrix [x is constitutedij']N×J, matrix [x is obtained after nondimensionalization processingij]N×J:
In formula (2), max xj' it is maximum value under same index;
A2: parameter entropy
K=(lnN)-1 (4)
In formula (3), fijThe specific gravity of the index is accounted for for lower i-th river of j-th of index;
In formula (4), K is constant;
In formula (5), HjFor the entropy of j-th of index;
A3: parameter entropy weight
Weighted sum method is used in above-mentioned steps S3, determines each first class index score in each river, specifically:
Assuming that first class index number is P,
In formula (7), SipFor i-th p-th of river first class index score, l is that the 1st second level refers under p-th of first class index
Corresponding j value is marked, m is two-level index number under p-th of first class index.
River cluster is carried out using hierarchical clustering method in above-mentioned steps S4, steps are as follows:
B1: distance calculates
Assuming that there is N river to participate in clustering, N class is established,It calculates, is obtained apart from square through distance
Battle array D(0),(0)To cluster original state;
The method that the distance calculates, including knearest neighbour method, longest distance method, group average distance method;
The knearest neighbour method, it is assumed that N1 and N2 is two rivers, N1 the and N2 shortest distance are as follows:
DN1,N2=min (du,v) (8)
In formula (8), DN1,N2The shortest distance between all first class index in all first class index in the river N1 and the river N2,
U is p-th of the river N1 first class index score, and v is p-th of the river N2 first class index score, du,vFor the distance between u and v;
The longest distance method, N1 and N2 longest distance are as follows:
DN1,N2=max (du,v) (9)
The group average distance method, N1 and N2 group average distance are as follows:
In formula (10), P is first class index number;
B2: new classification is established
Assuming that min0For Distance matrix D(0)Minimum value:
min0=min (D(0)) (11)
Assuming that min0It isWithThe distance between two classes,WithMerge into one kindOther classes not merged establish new classification:
B3: new classifying distance calculates
It calculatesWithThe distance between, obtain Distance matrix D(1);
B4: it computes repeatedly
Step B2-B3 is repeated, until N river cluster is 1, then stops clustering;
B5: cluster
Dendrogram is drawn, determines clusters number.
Above-mentioned sinuous degree is river bending degree:
In formula (12), β is the degree that wriggles, L1For Talweg length, L2For river upstream and downstream point-to-point transmission linear distance.
The above-mentioned vertical coefficient of stability of water body:
In formula (13), χ is the vertical coefficient of stability of water body, and g is acceleration of gravity, ρHFor bottom water body density, ρ0For superficial water
Volume density, ρavgFor the vertical averag density of water body, H is the river depth of water;
The water body density is calculated by the corresponding water body density of water temperature and water body silt content:
In formula (14), ρ is water body density, ρTFor the corresponding water body density of water temperature, ρsFor silt bulk density, δ is water body containing sand
Amount;
The corresponding water body density of the water temperature:
In formula (15), T is water temperature.
Above-mentioned light attenuation coefficient:
In formula (16), ε is light attenuation coefficient, and z is from the river water surface to measured place depth, and E (z) is to utilize underwater light quantum
Instrument measures the irradiation level of depth z, and E (0) is water surface irradiation level;
The euphotic zone depth:
In formula (17), zeuFor euphotic zone depth;
The layer depth takes river surface temperature to decline 0.5 DEG C of corresponding depth of water.
The invention has the beneficial effects that:
A kind of river clustering method based on algal bloom of the invention, comprehensively considers river algal bloom important indicator,
Including river form, water body power, water body physical and chemical factor, water body optics, water nutrition index, river algal bloom index is established
System, and each index weights are determined using subjective and objective combination weights method, river cluster finally is carried out using hierarchical clustering method, thus
The river of different zones, different characteristic is sorted out, understand river algal bloom universal law, be wawter bloom risk profile and
Wawter bloom prevention and control provide a kind of effective technical management tool.
Detailed description of the invention
Fig. 1 is a kind of flow chart of river clustering method based on algal bloom of the invention;
Fig. 2 is the structural schematic diagram of algal bloom index system in river of the invention;
Fig. 3 is river cluster result of the invention.
Specific embodiment
Specific introduce is made to the present invention below in conjunction with the drawings and specific embodiments.
As shown in Figure 1, the invention discloses a kind of river clustering method based on algal bloom, comprising the following steps:
S1, river algal bloom index system is established;
S2, system index weights are determined;
S3, parameter score;
S4, river cluster.
Algal bloom index system in river in step S1, including first class index layer and two-level index layer;
First class index layer includes river morphological index, water body dynamic index, water body physical and chemical factor index, water body optical index
With water nutrition index;
Two-level index layer is made of several two-level index;River morphological index includes length index, area index and wriggles
Index is spent, water body dynamic index includes that water level index, velocity parameters and the vertical coefficient of stability index of water body, water body physical and chemical factor refer to
Mark includes temperature index, pH value index, conductivity indices, oxidation-reduction potential index and DO index, and water body optical index includes
Transparency index, turbidity index, light attenuation coefficient index, euphotic zone depth index and layer depth index, water nutrition refer to
Mark includes total nitrogen (TN) index, total phosphorus (TP) index, dissolved silicon hydrochlorate (D-Si) index, nitrogen phosphorus (TN/TP) than index, silicon nitrogen
Than (D-Si/TN) index, silicon phosphorus ratio (D-Si/TP) index and chlorophyll (Chla) index;
The degree that wriggles is river bending degree:
In formula (12), β is the degree that wriggles, L1For Talweg length, L2For river upstream and downstream point-to-point transmission linear distance;
The vertical coefficient of stability of water body:
In formula (13), χ is the vertical coefficient of stability of water body, and g is acceleration of gravity, ρHFor bottom water body density, ρ0For superficial water
Volume density, ρavgFor the vertical averag density of water body, H is the river depth of water;
Water body density is calculated by the corresponding water body density of water temperature and water body silt content:
In formula (14), ρ is water body density, ρTFor the corresponding water body density of water temperature, ρsFor silt bulk density, δ is water body containing sand
Amount;
The corresponding water body density of water temperature:
In formula (15), T is water temperature;
Light attenuation coefficient:
In formula (16), ε is light attenuation coefficient, and z is from the river water surface to measured place depth, and E (z) is to utilize underwater light quantum
Instrument measures the irradiation level of depth z, and E (0) is water surface irradiation level;
Euphotic zone depth:
In formula (17), zeuFor euphotic zone depth;
Layer depth takes river surface temperature to decline 0.5 DEG C of corresponding depth of water.
System index weights are determined in step S2, it is assumed that two-level index number is J, using subjective and objective combination weights method:
wj=α aj+(1-α)bjJ=1,2 ..., J, 0≤α≤1 (1)
In formula (1), wjFor the comprehensive weight of j-th of index, ajAnd bjThe subjective weight of respectively j-th index and objective
Weight, α is preference coefficient, according to policymaker to the preference of different enabling legislations;
Subjective weighting method uses expert survey, and n experts understand river algal bloom index system, and anonymity is index j
Assign power, ajIt is averaged;
Objective weighted model uses Information Entropy, and steps are as follows:
A1: nondimensionalization processing
Nondimensionalization processing is carried out to index system data using Maximum Approach, it is assumed that there is N river to participate in clustering,
Assuming that two-level index number is J, matrix [x is constitutedij']N×J, matrix [x is obtained after nondimensionalization processingij]N×J:
In formula (2), maxxj' it is maximum value under same index;
A2: parameter entropy
K=(ln N)-1 (4)
In formula (3), fijThe specific gravity of the index is accounted for for lower i-th river of j-th of index;
In formula (4), K is constant;
In formula (5), HjFor the entropy of j-th of index;
A3: parameter entropy weight
Parameter score in step S3, it is assumed that first class index number is P, determines that each river is each using weighted sum method
First class index score:
In formula (7), SipFor i-th p-th of river first class index score, l is that the 1st second level refers under p-th of first class index
Corresponding j value is marked, m is two-level index number under p-th of first class index.
River clusters in step S4, carries out river cluster using hierarchical clustering method, steps are as follows:
B1: distance calculates
Assuming that there is N river to participate in clustering, N class is established,It calculates, is obtained apart from square through distance
Battle array D(0),(0)To cluster original state;Knearest neighbour method, longest distance method, group average distance may be selected in the method that distance calculates
Method;
Knearest neighbour method, it is assumed that N1 and N2 is two rivers, N1 the and N2 shortest distance are as follows:
DN1,N2=min (du,v) (8)
In formula (8), DN1,N2The shortest distance between all first class index in all first class index in the river N1 and the river N2,
U is p-th of the river N1 first class index score, and v is p-th of the river N2 first class index score, du,vFor the distance between u and v;
Longest distance method, N1 and N2 longest distance are as follows:
DN1,N2=max (du,v) (9)
Group average distance method, N1 and N2 group average distance are as follows:
In formula (10), P is first class index number;
B2: new classification is established
Assuming that min0For Distance matrix D(0)Minimum value:
min0=min (D(0)) (11)
Assuming that min0It isWithThe distance between two classes,WithMerge into one kindOther classes not merged establish new classification:
B3: new classifying distance calculates
It calculatesWithThe distance between, obtain Distance matrix D(1);
B4: it computes repeatedly
Return step B2, computes repeatedly, until N river cluster is 1, then stops clustering;
B5: cluster
Dendrogram is drawn, determines clusters number.
Embodiment 1:
By taking the numerous tributaries of a certain large-scale river-like reservoir as an example, the river cluster based on algal bloom is carried out, steps are as follows:
1. establishing river algal bloom index system
Algal bloom index system in river is as shown in Fig. 2, include first class index layer and two-level index layer;
First class index layer includes river morphological index, water body dynamic index, water body physical and chemical factor index, water body optical index
With water nutrition index, amount to 5 first class index;
Two-level index layer is made of two-level index, and river morphological index includes that length index, area index and the degree that wriggles refer to
Mark, water body dynamic index includes water level index, velocity parameters and the vertical coefficient of stability index of water body, water body physical and chemical factor index packet
Temperature index, pH value index, conductivity indices, oxidation-reduction potential index and DO index are included, water body optical index includes transparent
Spend index, turbidity index, light attenuation coefficient index, euphotic zone depth index and layer depth index, water nutrition index packet
Total nitrogen (TN) index, total phosphorus (TP) index, dissolved silicon hydrochlorate (D-Si) index, nitrogen phosphorus (TN/TP) are included than index, silicon nitrogen ratio
(D-Si/TN) index, silicon phosphorus ratio (D-Si/TP) index and chlorophyll (Chla) index amount to 23 two-level index.
2. determining system index weights
It shares 10 rivers and participates in clustering, by taking the 1st article of river as an example, river algal bloom index value such as 1 institute of table
Show;
Each index weights of expert survey, each index weights of Information Entropy and each index weights of subjective and objective combination weights method such as table
Shown in 2, preference factor alpha takes 0.5.
1 river algal bloom index value of table (by taking the 1st article of river as an example)
2 index weights of table
3. parameter score
It shares 10 rivers and participates in clustering, 23 two-level index, river algal bloom first class index score such as 3 institute of table
Show.
3 river algal bloom first class index score of table
4. river clusters
It is carried out using knearest neighbour method apart from calculating, Distance matrix D(0)As shown in table 4.
4 Distance matrix D of table(0)
Distance matrix D(0)In, 0.02 is minimum value, therefore one kind is merged into river 2, river 5 and river 6, it is assumed that close
And 1, other classes not merged establish new classification, Distance matrix D(1)As shown in table 5.
It is carried out using knearest neighbour method apart from calculating, by river 1 for merging 1 (river 2, river 5 and river 6), point
Not Ji Suan river 1 and river 2, river 5 and river 6 difference, be minimized.
5 Distance matrix D of table(1)
Distance matrix D(1)In, 0.03 is minimum value, therefore one kind is merged into river 2, river 5, river 6 and river 10, false
It is set as merging 2, other classes not merged establish new classification, Distance matrix D(2)As shown in table 6.
6 Distance matrix D of table(2)
Distance matrix D(2)In, 0.05 is minimum value, therefore river 2, river 5, river 6, river 7 and river 10 are merged into
One kind, other classes not merged establish new classification.And so on, when 10 rivers finally being merged into 1 (cluster is 1),
End of clustering.
As shown in figure 3, the river cluster result based on algal bloom,
When river is divided into two classes, river 3, river 8 are one kind, remaining river is one kind;
When river is divided into three classes, river 1 is one kind, and river 3, river 8 are one kind, remaining river is one kind.
Claims (7)
1. a kind of river clustering method based on algal bloom, which comprises the following steps:
S1, river algal bloom index system, including first class index layer and two-level index layer are established;
The first class index layer is made of several first class index, including river morphological index, water body dynamic index, water body physics and chemistry because
Sub- index, water body optical index, water nutrition index;
The two-level index layer is made of several two-level index,
River morphological index includes length index, area index, wriggles and spend index,
Water body dynamic index includes water level index, velocity parameters, the vertical coefficient of stability index of water body,
Water body physical and chemical factor index includes temperature index, pH value index, conductivity indices, oxidation-reduction potential index, DO index,
Water body optical index includes transparency index, turbidity index, light attenuation coefficient index, euphotic zone depth index, mixed layer
Depth index,
Water nutrition index includes total nitrogen (TN) index, total phosphorus (TP) index, dissolved silicon hydrochlorate (D-Si) index, nitrogen phosphorus (TN/
TP) than index, silicon nitrogen ratio (D-Si/TN) index, silicon phosphorus ratio (D-Si/TP) index, chlorophyll (Chla) index;
S2, using subjective and objective combination weights method, determine each two-level index weight in each river;
S3, each first class index score in each river is determined using weighted sum method based on each two-level index weight in each river;
S4, using hierarchical clustering method, carry out river cluster.
2. a kind of river clustering method based on algal bloom according to claim 1, which is characterized in that the step S2
In subjective and objective combination weights method include subjective weighting method, objective weighted model;
Determine system index weights, specifically:
Assuming that two-level index number is J:
wj=α aj+(1-α)bjJ=1,2 ..., J, 0≤α≤1 (1)
In formula (1), wjFor the comprehensive weight of j-th of index, ajAnd bjThe subjective weight and objective weight of respectively j-th index,
α is preference coefficient, according to policymaker to the preference of different enabling legislations;
The subjective weighting method uses expert survey, and n experts understand river algal bloom index system, and anonymity is index j
Assign power, ajIt is averaged;
The objective weighted model uses Information Entropy, and steps are as follows:
A1: nondimensionalization processing
Nondimensionalization processing is carried out to index system data using Maximum Approach, it is assumed that there is N river to participate in clustering, it is assumed that
Two-level index number is J, constitutes matrix [xij']N×J, matrix [x is obtained after nondimensionalization processingij]N×J:
In formula (2), maxxj' it is maximum value under same index;
A2: parameter entropy
K=(lnN)-1 (4)
In formula (3), fijThe specific gravity of the index is accounted for for lower i-th river of j-th of index;
In formula (4), K is constant;
In formula (5), HjFor the entropy of j-th of index;
A3: parameter entropy weight
3. a kind of river clustering method based on algal bloom according to claim 1, which is characterized in that the step S3
It is middle to use weighted sum method, determine each first class index score in each river, specifically:
Assuming that first class index number is P,
In formula (7), SipFor i-th p-th of river first class index score, l is the 1st two-level index pair under p-th of first class index
The j value answered, m are two-level index number under p-th of first class index.
4. a kind of river clustering method based on algal bloom according to claim 1, which is characterized in that the step S4
Middle to carry out river cluster using hierarchical clustering method, steps are as follows:
B1: distance calculates
Assuming that there is N river to participate in clustering, N class is established,It is calculated through distance, obtains Distance matrix D(0), (0) is cluster original state;
The distance calculates, including knearest neighbour method, longest distance method, group average distance method;
The knearest neighbour method, it is assumed that N1 and N2 is two rivers, N1 the and N2 shortest distance are as follows:
DN1,N2=min (du,v) (8)
In formula (8), DN1,N2The shortest distance between all first class index in all first class index in the river N1 and the river N2, u are
P-th of the river N1 first class index score, v are p-th of the river N2 first class index score, du,vFor the distance between u and v;
The longest distance method, N1 and N2 longest distance are as follows:
DN1,N2=max (du,v) (9)
The group average distance method, N1 and N2 group average distance are as follows:
In formula (10), P is first class index number;
B2: new classification is established
Assuming that min0For Distance matrix D(0)Minimum value:
min0=min (D(0)) (11)
Assuming that min0 isWithThe distance between (i=1,2 ..., N, r=1,2 ..., N, i ≠ r) two classes,With
Merge into one kindOther classes not merged establish new classification:Q=N-2;
B3: new classifying distance calculates
It calculatesWithThe distance between, obtain Distance matrix D(1);
B4: it computes repeatedly
Step B2-B3 is repeated, until N river cluster is 1, then stops clustering;
B5: cluster
Dendrogram is drawn, determines clusters number.
5. a kind of river clustering method based on algal bloom according to claim 1, which is characterized in that described to wriggle
Degree is river bending degree:
In formula (12), β is the degree that wriggles, L1For Talweg length, L2For river upstream and downstream point-to-point transmission linear distance.
6. a kind of river clustering method based on algal bloom according to claim 1, which is characterized in that the water body hangs down
To the coefficient of stability:
In formula (13), χ is the vertical coefficient of stability of water body, and g is acceleration of gravity, ρHFor bottom water body density, ρ0It is close for surface water
Degree, ρavgFor the vertical averag density of water body, H is the river depth of water;
The water body density is calculated by the corresponding water body density of water temperature and water body silt content:
In formula (14), ρ is water body density, ρTFor the corresponding water body density of water temperature, ρsFor silt bulk density, δ is water body silt content;
The corresponding water body density of the water temperature:
In formula (15), T is water temperature.
7. a kind of river clustering method based on algal bloom according to claim 1, which is characterized in that the optical attenuation
Coefficient:
In formula (16), ε is light attenuation coefficient, and z is from the river water surface to measured place depth, and E (z) is to be surveyed using underwater quanta meter
The irradiation level of depth z is obtained, E (0) is water surface irradiation level;
The euphotic zone depth:
In formula (17), zeuFor euphotic zone depth;
The layer depth takes river surface temperature to decline 0.5 DEG C of corresponding depth of water.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811165966.5A CN109272246B (en) | 2018-10-08 | 2018-10-08 | River channel clustering method based on bloom outbreak |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811165966.5A CN109272246B (en) | 2018-10-08 | 2018-10-08 | River channel clustering method based on bloom outbreak |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109272246A true CN109272246A (en) | 2019-01-25 |
CN109272246B CN109272246B (en) | 2021-11-30 |
Family
ID=65195689
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811165966.5A Active CN109272246B (en) | 2018-10-08 | 2018-10-08 | River channel clustering method based on bloom outbreak |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109272246B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112802069A (en) * | 2020-12-30 | 2021-05-14 | 北京师范大学 | Method for judging dam bank of shallow lake water-blocking dam |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100710445B1 (en) * | 2005-07-21 | 2007-04-24 | 배연재 | Stream water quality assessment technique using benthic macroinvertebrates |
CN106202960A (en) * | 2016-07-21 | 2016-12-07 | 沈阳环境科学研究院 | A kind of health assessment method based on Lake Water ecosystem |
CN107122927A (en) * | 2017-06-27 | 2017-09-01 | 河海大学 | A kind of water transfer drainage water environment improvement integrated evaluating method |
-
2018
- 2018-10-08 CN CN201811165966.5A patent/CN109272246B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100710445B1 (en) * | 2005-07-21 | 2007-04-24 | 배연재 | Stream water quality assessment technique using benthic macroinvertebrates |
CN106202960A (en) * | 2016-07-21 | 2016-12-07 | 沈阳环境科学研究院 | A kind of health assessment method based on Lake Water ecosystem |
CN107122927A (en) * | 2017-06-27 | 2017-09-01 | 河海大学 | A kind of water transfer drainage water environment improvement integrated evaluating method |
Non-Patent Citations (1)
Title |
---|
彭进平等: "水动力条件对湖泊水体磷素质量浓度的影响", 《生态环境》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112802069A (en) * | 2020-12-30 | 2021-05-14 | 北京师范大学 | Method for judging dam bank of shallow lake water-blocking dam |
CN112802069B (en) * | 2020-12-30 | 2023-07-04 | 北京师范大学 | Determination method for water-blocking dike cofferdam of shallow lake |
Also Published As
Publication number | Publication date |
---|---|
CN109272246B (en) | 2021-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chang et al. | A novel index for assessing the water quality of urban landscape lakes based on water transparency | |
Hohenegger et al. | Remarks on west Pacific Nummulitidae (foraminifera) | |
Cavender-Bares et al. | Microbial size spectra from natural and nutrient enriched ecosystems | |
Lande et al. | Suspension times of particles in the upper ocean | |
Hu et al. | Comparing biological classifications of freshwater phytoplankton: a case study from South China | |
Lee et al. | Eutrophication dynamics of Tolo Harbour, Hong Kong | |
Bowman et al. | Tidal stirring and the distribution of phytoplankton in Long Island and Block Island Sounds | |
CN109117861B (en) | Point set multi-level aggregative property analysis method considering spatial position | |
CN110969346A (en) | Drainage basin water ecological function regional treatment demand evaluation method based on index screening | |
Baird et al. | The use of mechanistic descriptions of algal growth and zooplankton grazing in an estuarine eutrophication model | |
La Ferla et al. | Are prokaryotic cell shape and size suitable to ecosystem characterization? | |
Udovič et al. | Co-occurrence of functional groups in phytoplankton assemblages dominated by diatoms, chrysophytes and dinoflagellates | |
CN109272246A (en) | A kind of river clustering method based on algal bloom | |
Centis et al. | Structure of the diatom community of the River Adige (North-Eastern Italy) along a hydrological gradient | |
Sui et al. | Seasonal water quality changes and the eutrophication of lake Yilong in Southwest China | |
Jiang et al. | Analysis of Spatial and Temporal Evolution of Regional Water Resources Carrying Capacity and Influencing Factors—Anhui Province as an Example | |
Liu et al. | Spatial and seasonal dynamics of water environmental capacity in mountainous rivers of the Southeastern Coast, China | |
Wang et al. | Monitored and modeled correlations of sediment and nutrients with Chesapeake Bay water clarity | |
CN106813645B (en) | A kind of section integrated evaluating method of healthy stream | |
Yang et al. | Environmental factors controlling the spatial distribution of subfossil Chironomidae in surface sediments of Lake Dongping, a warm temperate lake in North China | |
Liu et al. | Effect of nutrient supply on cell size evolution of marine phytoplankton | |
Ge et al. | Development and testing of a planktonic index of biotic integrity (P-IBI) for Lake Fuxian, China | |
Shen et al. | Distribution and interaction characteristics of water quality at the stratified confluence reservoirs | |
Luo et al. | Evaluation of the suitability of human settlements in Qingpu district, Shanghai | |
Gin | Microbal size spectra from diverse marine ecosystems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |