一种基于水生态分区理论的河口适宜生态流量评估方法与流程

一种基于水生态分区理论的河口适宜生态流量评估方法
1.本发明属于环境保护技术领域，具体为一种基于水生态分区理论的河口适宜生态流量评估方法。


背景技术：

2.河口是河流与海洋交汇处，在河口和海洋之间起到物质循环、能量流动、信息传递等的生态服务功能，在维持生态系统持续发展等方面起到重要的作用(zhang et al.，2017)。近年来，在人类开发活动和自然因素的综合影响下，河口生态系统退化及生态功能显著下降的问题日益突出。在人类开发活动中，如水利工程、围填海工程、河口采砂、森林砍伐、工农业污水等，这些活动都影响入海径流的水量、水质和时空变化。入海径流的改变导致河口生态系统退化及生态功能显著下降的问题日益突出。因此，开展河口生态流量的研究，有助于指导生态系统的修复和保护，对其可持续发展有重大意义。
3.生态流量的研究开始于20世纪40年代，目前已有将近50个国家开展了相关研究，主要研究对象有湖泊、河流和湿地等淡水生态系统。迄今为止，根据不同的流量定义，生态流量的计算方法超过200种，主要可以分为：水文学法、水力学法、生境模拟法和整体分析法等(王中根等，2020)。而针对淡水输入对河口生态系统影响的研究从20世纪80年代开始兴起，近30年来不断发展和完善。河口生态系统受到海陆间的交互作用，水动力过程复杂，影响因素众多，决定了其生态流量机理和计算方法的复杂性，相应造成了淡水生态系统的生态流量研究成果直接应用于河口生态系统中存在适用性问题。由于径流和潮汐交互作用，水环境各类指标在不同河口甚至同一河口具有明显的区域特征(刘静等，2017)。河口上、中、下游段的水体理化指标的差异最终反映生物群落分布。传统的水生态系统监测管理、保护和管理大多只基于政治行政单位和水文分区单元，导致其分散而片面，未能从生态系统的整体角度上揭示水生生态系统状况，进行合适的分区管理。目前，国际上淡水资源管理的基本空间单元主要有3类，即行政分区单元，流域水文单元和生态分区单元。其中，水生态分区是根据水生生物与其承受的环境压力之间的响应关系，从多个尺度上划分出具有相同生态潜力的区域，在生态治理中更具有优势(孙小银等，2010)。杨爱民等(2008)的中国生态水文分区研究，将生态流量因素与水生态分区结合，但并没有直接考虑到水生生物等重要因素，亦同样属于水文分区范畴。淡水输入对河口水环境过程的影响已成为河口生物种群分布的变化的主要因素。美国旧金山海湾实地调查资料分析认为，河口淡水输入及相应盐度变化度对生物群落分布及生物量具有显著影响(powell et al.,2002)。生物体能够整合环境的时间尺度上压力效应信息，可能比简单的物理化学指标更能反应生态系统的时间变化特征，能够从生物学角度对环境状况进行监测和评价。浮游植物作为水生态系统中重要的初级生产者，其群落结构特征能够敏感地反应水生态系统对水环境污染和生态破坏的响应。reynolds等和padis
á
k等依据浮游植物对特定的环境因子的敏感性和耐受性，将在某一已知环境中共同占优势的不同种群的浮游植物归纳同一功能类群(reynolds et al.,2002,padis
á
k et al.,2009)。很大程度弥补了浮游植物个体组成随机性大、种类适应生境互相重叠等不足。目前，浮游植物功能类群方法已经广泛运用到在湖泊、河流、水库及湿地
等水域生态系统中(贾鹏等，2021)。根据水生生物与承受的环境压力之间的响应关系，选择生物类群(如：浮游植物类群)为分区指标，从多个维度划分出具有相同生态性质的区域，根据各分区单元的情况，制定合适生态流量方案，确定优先监测和保护的区域，推断未检测区域生态走向和环境状况等，是其他分区方法所不具备的，也是当前环境管理所迫切需要的。
4.总体来说，河口独特地理环境下，不同区域的某些物理和化学因素以及生物群落结构差异明显，在实际的水资源管理中，现有大多数生态流量方法较少考虑针对生物类群以及不同生态区域进行分区管理，不利于河口生态系统的环境保护，生态恢复，以及水资源的科学规划管理。鉴于上述原因，我们提出了一种以生物类群划分的水域水生态分区，为河口水生态系的研究、评价、修复和管理提供合适的空间单元，解决实际应用问题的生态流量评估方法。


技术实现要素：

5.本发明的目的在于基于水生态分区理论，提出以相似浮游植物类群的河口水生态分区，根据影响浮游植物类群组成的关键环境因子空间分布情况表征水生态分区空间变化，建立与生态流量之间响应关系的分析方法。
6.为了实现上述目的，本发明采用了如下技术方案：
7.一种基于水生态分区理论的河口适宜生态流量评估方法，具体包括以下步骤：
8.s1、确定研究区域，设计研究方案，收集研究区域的水文数据，同时采集研究区域现场的浮游植物数据和环境数据信息；
9.s2、结合s1中所得的浮游植物数据信息，通过定性分析方法和定量分析方法完成浮游植物的物种鉴定；
10.s3、通过s2中所述的定性分析数据采用聚类分析方法确定研究区域浮游植物的群落类型和水生态分区的划分；
11.s4、通过s2中所述的定量分析数据采用约束排序分析方法确定研究区域的关键环境因子和浮游植物类群的代表物种；
12.s5、结合s1、s3与s4中所得信息，构建不同流量情景下的水动力-水质二维模型；
13.s6、根据s5中所得的水动力-水质二维模型的结果，分析确定研究区域的适宜生态流量。
14.优选地，所述s1中提到的水文数据收集，包括研究区域不同时期的径流量以及潮位信息。同时采集浮游植物数据和环境数据。
15.优选地，所述s3中提到确定研究区域浮游植物的群落类型和水生态分区的划分，具体包括以下步骤：
16.a1、依据浮游植物定性数据进行聚类分析，确定不同类型的浮游植物群落；
17.a2、根据不同浮游植物类群所在采样点划分不同水生态分区。
18.优选地，所述s4中提到的确定研究区域的关键环境因子和浮游植物类群的代表物种，具体包括以下步骤：
19.b1、进行排序分析前，使用浮游植物丰度数据进行去趋势对应分析(dca)来确定是采用线性排序还是单峰排序；若dca四个轴中第一长轴的数值小于3，则选用冗余分析(rda)；若数值大于4，则选用典范对应分析(cca)；若数值在3-4之间，则两者都合适；
20.b2、数据分析时，将所有物种数据均通过lg(x 1)进行转化；
21.b3、去除方差膨胀系统(vif)大于20的环境变量，避免环境变量之间的多重共线性造成模型的不稳定；
22.b4、关键环境因子的确定：关键环境因子通过多元回归变量筛选，选择forward.sel函数进行变量的前向选择；
23.b5、通过浮游植物马格列夫指数和关键环境因子作广义可加模型，得到浮游植物的马格列夫指数和关键环境因子的响应关系，选择马格列夫指数大于其均值的数值对应的环境因子数值作为适宜值；
24.b6、浮游植物类群的代表物种选择：基于物种与约束排序轴的得分大小，确定合适的数量的代表物种。
25.优选地，所述b4中提到的变量的前向选择，具体包括以下步骤：
26.c1、依次分别运行每个环境变量与响应变量的约束排序分析；
27.c2、基于设置的标准选择“最好”的显著的解释变量；
28.c3、寻找模型中第二(第三、第四)解释变量；
29.c4、重复c1-c3所述操作，直至无显著性的解释性变量为止。
30.优选地，所述s5中提到的不同流量情景下的水动力-水质二维模型的构建，具体包括以下步骤：
31.d1、模型数据收集和前处理：构建研究研究区域的网格地形文件；
32.d2、搭建模型：设置内河口上边界和外河口下边界，模型输入主要条件：模拟物质，扩撒系统，降解系数，初始条件；
33.d3、模型率定和验证：通过保守物种对扩散系数进行率定，通过模拟区域不同监测点的降解物质浓度对降解速率进行率定；
34.d4、方案设计：设置不同流量情景下的水生态分区时空变化，通过查阅资料，确定合理的等梯度的流量变化以及流量范围。
35.优选地，所述s6中提到的研究区域的适宜生态流量的确定，具体包括以下步骤：
36.e1、计算各生态分区的面积：根据不同流量情景模拟下的结果，提取每个网格中各个关键环境因子模拟值，按照各分区的关键环境因子适宜值划分，然后将对应分区不同种类关键环境因子划分网格进行对照，取它们之间网格重合部分定义为该区生态类型的区域，不重合部分划分为其他生态类型区域，最后计算网格的面积得出各生态类群区域的面积与空间位置，该研究区各生态分区的面积以下公式表示：
37.s
总
＝s1 s2 s3
…
sn38.其中：s
总
表示研究区域的总面积；
39.s1表示ⅰ类型生态分区；
40.s2表示ⅱ类型生态分区；
41.s3表示ⅲ类型生态分区；
42.……
43.sn表示其他类型的生态分区；
44.e2、通过改变流量输入情景，可以得到一系列梯度的不同流量输入下，水生态分区的面积大小和空间位置变化响应关系，根据水生态分区的时空变化，依据生态系统健康和
水资源利用的要求，选择合适的生态流量模式。
45.与现有技术相比，本发明提供了一种基于水生态分区理论的河口适宜生态流量评估方法，具备以下有益效果：
46.(1)本发明考虑了不同浮游植物类群划分的水生态分区的生态流量分析技术，能较细致的了解河口各区域生态的情况，综合考虑各水生态分区的生态问题，避免生态流量计算结果的片面性。
47.(2)本发明结合水生态分区与生态流量的响应关系，可以为研究区确定优先监测和保护的区域，推断未监测区域生态走向与环境状态等，提高河口生态治理的有效性和可行性。
附图说明
48.图1为本发明提出的一种基于水生态分区理论的河口适宜生态流量评估方法的方法流程示意图；
49.图2为本发明提出的一种基于水生态分区理论的河口适宜生态流量评估方法的实施例3中浮游植物ward最小方差聚类分析结果示意图；
50.图3为本发明提出的一种基于水生态分区理论的河口适宜生态流量评估方法的实施例3中基于相似浮游植物类群的水生态分区划分示意图；
51.图4为本发明提出的一种基于水生态分区理论的河口适宜生态流量评估方法的实施例3中浮游植物与监测环境因子的rda排序示意图；
52.图5为本发明提出的一种基于水生态分区理论的河口适宜生态流量评估方法的实施例3中浮游植物与关键环境因子的rda排序示意图；
53.图6为本发明提出的一种基于水生态分区理论的河口适宜生态流量评估方法的实施例3中浮游植物丰富度与关键环境因子的响应关系示意图；
54.图7为本发明提出的一种基于水生态分区理论的河口适宜生态流量评估方法的实施例3中模型的网格地形图；
55.图8为本发明提出的一种基于水生态分区理论的河口适宜生态流量评估方法的实施例3中模型边界设置示意图；
56.图9为本发明提出的一种基于水生态分区理论的河口适宜生态流量评估方法的实施例3中同时期的40％、80％、100％、120％、160％的流量情景盐度场变化示意图；
57.图10为本发明提出的一种基于水生态分区理论的河口适宜生态流量评估方法的实施例3中各水生态分区在不同流量情景下的响应曲线。
具体实施方式
58.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例只是本发明一部分实施例，而不是全部的实施例。
59.实施例1：
60.请参阅图1，一种基于水生态分区理论的河口适宜生态流量评估方法，具体包括以下步骤：
61.s1、确定研究区域，设计研究方案，收集研究区域的水文数据，同时采集研究区域
现场的浮游植物数据和环境数据信息；
62.s1中提到的水文数据收集，包括研究区域不同时期的径流量以及潮位信息；
63.s2、结合s1中所得的浮游植物数据信息，通过定性分析方法和定量分析完成浮游植物的物种鉴定；
64.s3、通过s2中的定性分析方法确定研究区域浮游植物的群落类型和水生态分区的划分；
65.s3中提到确定研究区域浮游植物的群落类型和水生态分区的划分，具体包括以下步骤：
66.a1、依据浮游植物定性数据进行聚类分析，确定不同类型的浮游植物群落；
67.a2、根据不同浮游植物类群所在采样点划分不同水生态分区；
68.s4、通过s2中的定量分析方法确定研究区域的关键环境因子和浮游植物类群的代表物种；
69.s4中提到的确定研究区域的关键环境因子和浮游植物类群的代表物种，具体包括以下步骤：
70.b1、进行排序分析前，使用浮游植物丰度数据进行去趋势对应分析(dca)来确定是采用线性排序还是单峰排序；若dca四个轴中第一长轴的数值小于3，则选用冗余分析(rda)；若数值大于4，则选用典范对应分析(cca)；若数值在3-4之间，则两者都合适；
71.b2、数据分析时，将所有物种数据均通过lg(x 1)进行转化；
72.b3、去除方差膨胀系统(vif)大于20的环境变量，避免环境变量之间的多重共线性造成模型的不稳定；
73.b4、关键环境因子的确定：关键环境因子通过多元回归变量筛选，选择forward.sel函数进行变量的前向选择；
74.b4中提到的变量的前向选择，具体包括以下步骤：
75.c1、依次分别运行每个环境变量与响应变量的约束排序分析；
76.c2、基于设置的标准选择“最好”的显著的解释变量；
77.c3、寻找模型中第二(第三、第四)解释变量；
78.c4、重复c1-c3操作，直至无显著性的解释性变量为止；
79.b5、通过浮游植物马格列夫指数和关键环境因子作广义可加模型，得到浮游植物的马格列夫指数和关键环境因子的响应关系，选择马格列夫指数大于其均值的数值对应的环境因子数值作为适宜值；
80.b6、浮游植物类群的代表物种选择：基于物种与约束排序轴的得分大小，确定合适的数量的代表物种；比如选择前两轴的累计拟合度超过0.6(这是个随意的值)的物种，拟合优度越高，物种与相对应的轴越匹配，也就是该物种能被所采集的环境因子解释效果越好；
81.s5、结合s3与s4中所得信息，构建不同流量情景下的水动力-水质二维模型；
82.s5中提到的不同流量情景下的水动力-水质二维模型的构建，具体包括以下步骤：
83.d1、模型数据收集和前处理：构建研究研究区域的网格地形文件；
84.d2、搭建模型：设置内河口上边界和外河口下边界，模型输入主要条件：模拟物质，扩撒系统，降解系数，初始条件；
85.d3、模型率定和验证：通过保守物种对扩散系数进行率定，通过模拟区域不同监测
点的降解物质浓度对降解速率进行率定；
86.d4、方案设计：设置不同流量情景下的水生态分区时空变化，通过查阅资料，确定合理的等梯度的流量变化以及流量范围；
87.s6、根据s5中所得的水动力-水质二维模型，分析确定研究区域的适宜生态流量；
88.s6中提到的研究区域的适宜生态流量的确定，具体包括以下步骤：
89.e1、计算各生态分区的面积：根据不同流量情景模拟下的结果，提取每个网格中各个关键环境因子模拟值，按照各分区的关键环境因子适宜值划分，然后将对应分区不同种类关键环境因子划分网格进行对照，取它们之间网格重合部分定义为该区生态类型的区域，不重合部分划分为其他生态类型区域，最后计算网格的面积得出各生态类群区域的面积与空间位置，该研究区各生态分区的面积以下公式表示：
90.s
总
＝s1 s2 s3
…
sn91.其中：s
总
表示研究区域的总面积；
92.s1表示ⅰ类型生态分区；
93.s2表示ⅱ类型生态分区；
94.s3表示ⅲ类型生态分区；
95.……
96.sn表示其他类型的生态分区；
97.e2、通过改变流量输入情景，可以得到一系列梯度的不同流量输入下，水生态分区的面积大小和空间位置变化响应关系，根据水生态分区的时空变化，依据生态系统健康和水资源利用的要求，选择合适的生态流量模式。
98.本发明针对目前较少采用生物类群为生物指标的水生态分区与河口生态流量的计算方法，基于水生态分区理论，提出以相似浮游植物类群的河口水生态分区，根据影响浮游植物类群组成的关键环境因子空间分布情况表征水生态分区空间变化，建立与生态流量之间响应关系的分析技术。与现有技术相比本发明有效效果在于：1、考虑不同浮游植物类群划分的水生态分区的生态流量分析技术，能较细致的了解河口各区域生态的情况，综合考虑各水生态分区的生态问题，避免生态流量计算结果的片面性；2、结合水生态分区与生态流量的响应关系，可以为研究区确定优先监测和保护的区域，推断未监测区域生态走向与环境状态等，提高河口生态治理的有效性和可行性。
99.实施例2：
100.请参阅图1，基于实施例1但有所不同之处在于，
101.1.研究区域的数据收集
102.主要分为两种类型：
103.(1)水文数据收集，包括不同时期径流量、潮位等水文数据；
104.(2)现场浮游植物数据和环境数据的采集：
105.①
通过现场调查，设置合理的采样点位置。如采样点的布设尽可能涵盖到每个出海口，依据河口的上、中、下游段布设合理梯度的采样点；
106.②
环境调查采样和检测方法按《海洋监测规范》(gb17378-2007)、《海洋调查规范》(gb 12763-2007)、《水与废水监测分析方法(第四版增补版)》与《近岸海域环境监测规》(hj442-2008)规范执行。
107.2.浮游植物物种鉴定
108.浮游植物鉴定分为定性和定量分析。
109.(1)定性数据根据一个物种是否存在于一个样方中，存在记为1，不存在记为0，就构成了二元生态数据，这种二元数据有着重要的生态意义，因为种出现与否与环境密切关系。同时避免了使用物种丰度数据出现极端的数据干扰。定性数据用于聚类分析，可以区分不同浮游植类群。
110.(2)定量数据是计算每个样点的单位体积中不同物种浮游植物的生物体个数，即浮游植物的丰度，而浮游植物丰度大小可以作为生物响应变量，结合环境因子(解释变量)约束排序分析，可以研究物种丰度与环境因子的之间的关系。
111.3.浮游植物群落类型与水生态分区的划分
112.聚类分析的目的是识别在环境中不连续的对象子集，它可以解决哪些对象有足够的相似性能够被归于一组，并且确定组与组之间的差异或分离程度。聚类方法因其算法而不同，本发明基于一种最小二乘法线性模型准则的聚类方法，分组的依据是使组内的平方和(即方差分析的方差最小化)，具体方法如下：
113.(1)依据浮游植物定性数据进行聚类分析，确定不同类型的浮游植物群落；
114.(2)根据不同浮游植物类群所在采样点划分不同水生态分区。
115.4.关键环境因子和浮游植物类群的代表物种的确定
116.为了探究不同环境因子与浮游植物群落结构之间的关系，使用约束排序，该方法是多元回归与传统排序的组合，明确地探索响应变量矩阵(物种数据)和解释变量矩阵(环境数据)两个矩阵之间的关系，在动植物群落与环境因子关系种的研究得到广泛应用；本发明用于排序的物种要求满足下面两个条件：该物种在各样点出现的频度＞12.5％,该物种在至少一个样点的相对密度≥1％，物种矩阵经过log(x 1)转换，具体操作如下：
117.(1)在进行排序分析前。先使用浮游植物丰度数据进行去趋势对应分析(dca)来确定是采用线性排序还是单峰排序。dca四个轴中第一长轴的数值小于3，则选用冗余分析(rda),数值大于4，则选用典范对应分析(cca)，数值在3-4之间，则两者都合适。
118.(2)数据分析时，所有物种数据均经过了lg(x 1)的转化。
119.(3)为了避免环境变量之间的多重共线性造成模型的不稳定。去除了方差膨胀系统(vif)大于20的环境变量。
120.(4)关键环境因子的确定：关键通过多元回归变量筛选，通常有三种模型，前向、后向和逐步变量选择；
121.本发明选择数据分析工具r语言软件中adespatial包中的函数forward.sel进行变量的前向选择，选择过程如下：
122.·
依次分别运行每个环境变量与响应变量的约束排序分析；
123.·
基于设置的标准选择“最好”的显著的解释变量；
124.·
接下来寻找模型中第二个(第三、第四)解释变量；
125.·
这个过程一直持续，直到无显著性的解释变量为止。
126.(5)通过浮游植物马格列夫指数和关键环境因子作广义可加模型，得到浮游植物的马格列夫指数和关键环境因子的响应关系，选择马格列夫指数大于其均值的数值对应的环境因子数值作为适宜值；
127.(6)浮游植物类群的代表物种选择：基于物种与约束排序轴的得分大小，确定合适的数量的代表物种。比如选择前两轴的累计拟合度超过0.6(这是个随意的值)的物种，拟合优度越高，物种与相对应的轴越匹配，也就是该物种能被所采集的环境因子解释效果越好。
128.5.基于不同流量情景的水动力-水质模型的模拟
129.河口区域环境因子的变化主要是由于入海口径流和海洋潮汐的影响，根据历史流量为参考值设置河口淡水输入，设置不同流量情景，流量设置可以分为水量和水质两部分组合在同一潮位下，不同流量大小下模拟相同时长，计算出不同流量情景下的水生态分区面积变化，具体内容如下：
130.(1)模型数据收集和前处理，构建研究研究区域的网格地形文件；
131.(2)搭建模型，设置内河口上边界和外河口下边界，模型输入主要参数条件：模拟物质，扩撒系统，降解系数，初始条件等；
132.(3)模型率定和验证。其中扩散系数和降解系数是属于率定系数，扩散系数可以通过保守物种，如盐度和示踪剂的进行率定。降解速率可以通过模拟区域不同监测点的降解物质浓度(如cod)进行率定，也可以使用相关文献中类似的区域降解物质的降解速率；
133.(4)方案设计：设置不同流量情景下的水生态分区时空变化。通过查阅资料，确定合理的等梯度的流量变化以及流量范围。
134.6.水生态分区面积与流量之间的响应关系
135.计算各生态分区的面积。根据不同流量情景模拟下的结果，提取每个网格中各个关键环境因子模拟值，按照各分区的关键环境因子适宜值划分，然后将对应分区不同种类关键环境因子划分网格进行对照，取它们之间网格重合部分定义为该区生态类型的区域，不重合部分划分为其他生态类型区域，最后计算网格的面积得出各生态类群区域的面积与空间位置，该研究区各生态分区的面积以下公式表示：
136.s
总
＝s1 s2 s3
…
sn137.其中：s
总
表示研究区域的总面积；
138.s1表示ⅰ类型生态分区；
139.s2表示ⅱ类型生态分区；
140.s3表示ⅲ类型生态分区；
141.……
142.sn表示其他类型生态分区；
143.通过改变流量输入情景，可以得到一系列梯度的不同流量输入下，水生态分区的面积大小和空间位置变化响应关系，根据水生态分区的时空变化，依据生态系统健康和水资源利用的要求，选择合适的生态流量模式。
144.实施例3：
145.请参阅图2-10，基于实施例1-2但有所不同之处在于，
146.本发明以2021年1月珠江东四口门作为实例分析，确定了三个浮游植物类群以及盐度、ph两个关键环境因子，以相似浮游植物类群划分为三个水生态分区，使用影响浮游植物类群组成的关键环境因子空间分布情况表征水生态分区空间分布，根据水动力-水质模模拟不同情景下流量与水生态分区的响应关系，结果如附图10所示，模拟了同时期的40％、60％、80％、100％、120％、140％和160％的径流量，在相同潮位边界情景下，不同流量下各
分区的之间响应关系。
147.本发明具体操作流程如下：
148.(1)基于浮游植物类群进行水生态分区划分
149.使用各采样点的浮游植物定性数据进行ward聚类，图2为聚类距离绘制的聚类树，通过采样点间聚类分组进行水生态分区划分，分区空间分布如图3所示，即采样点s8所在的区域为微咸水生态分区，采样点s5所在的区域为河海交汇水生态分区，采样点s21所在的区域为近岸海水生态分区。
150.(2)基于冗余分析(rda)的关键因子和代表物种的确定
151.将浮游植物数据与环境数据经过上述流程中的数据转化方法进行处理，得到生物数据(响应变量)和环境数据(解释变量)的两个数据矩阵，使用数据分析工具r语言软件中vegan包中rda分析得到两矩阵的之间关系，提取、解读和绘制rda的排序结果(如图4所示)，通过前向选择可以确定关键环境因子，绘制关键环境因子与浮游植物的简约型排序图(如图5所示)。
152.(3)关键环境因子适宜数值范围的选择
153.环境因子的适宜范围通过生物指标和环境因子的响应关系来确定。生物多样性指数是指用简单的数值表示群落内种类多样性的程度，是用来判断群落或生态系统的稳定性指标。其中马格列夫指数是考虑群落的物种数量和总个体数，将一定大小的样本中的物种数量定义为多样性指数，马格列夫指数越大，物种丰富度越大。相比于单一物种数或生物量更加综合体现浮游植物群落与环境因子的关系，考虑生物与环境之间关系不是简单线性加和，所以本发明采用广义可加模型(gams)模型来模拟生物指数与环境因子的响应关系，确定关键环境因子的适宜范围。如图6所示，分别为三个水生态分区的马格列夫指数与盐度和ph两个因子的响应曲线，选择各分区内的马格列夫指数大于均值部分对应的环境因子数值作为适宜值，各分区的关键环境因子的适宜范围如下表所示：
154.表1各生态分区关键环境因子的适宜范围
[0155][0156][0157]
(4)基于不同流量情景的水动力-水质模型的模拟
[0158]
将收集的地形数据通过mike21 mesh generator网格生成器生态网格地形图(如图7所示)，分别设置上下游的开边界，上游开边界设置有虎门、洪奇门、蕉门和横门四个入海的河道边界，下游开边界设置以澳门-桂山岛-香港的外海边界(如图8所示)。设置好模型各参数后，进行模型的率定和验证。基于率定好模型，在其他参数不变的条件下，模拟同时期的40％、60％、80％、100％、120％、140％和160％生态流量情景，得到不同流量情景条件下的关键环境因子变化结果，图9所示为同时期的40％、80％、100％、120％、160％的流量情景盐度场变化图。
[0159]
(5)水生态分区面积与流量之间的响应关系
[0160]
根据不同流量情节条件下的关键环境因子空间变化结果，统计各分区的关键环境因子的适宜值范围的网格面积，根据各水生态分区的盐度和ph适宜值范围，确定各分区的
不同流量下的分区面积(如表2和图10所示)。
[0161]
表2各分区面积与流量之间的响应关系
[0162][0163]
注：q为日均径流量、面积单位为km2。
[0164]
从图10中可以得出ⅰ类型分区(微咸水生态类型)的面积在40％-120％日均流量的情景下，其面积随流量变化明显，而超过日均流量120％后，其面积随流量变化幅度减小；ⅱ类型分区(河海交汇水生态类型)的面积随流量变不明显，但模型模拟空间结果看，空间位置随流量的增大向外海移动；ⅲ类型分区(近海岸水生态类型)，则是随流量的增大而减小，其中是因为分区区域因流量增大，导致分区空间位置向外海移动，移出了研究区域所导致面积减小。
[0165]
以上，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。