一种基于内源释放的大型浅水湖泊模型及其构建方法技术

本发明专利技术公开了一种基于内源释放的大型浅水湖泊模型及其构建方法，构建方法包括以下步骤：资料收集及实地调查监测；分区取样测定沉积物特性；构建水动力模块并进行边界条件和参数敏感性分析；构建水质模块并进行参数敏感性分析；构建沉积物模块；构建沉积物成岩模块并进行参数敏感性分析；判定沉积物再悬浮的底部切应力条件，耦合沉积物和水质模块。本发明专利技术将实地调查、地理信息技术及数值模拟有机结合，考虑水动力、水质和生态要素，合理分析大型浅水湖泊在不同水动力、边界作用下湖区及底泥营养盐变化，为浅水湖泊保护、利用、技术咨询、管理提供依据。

A Large Shallow Lake Model Based on Endogenous Release and Its Construction Method

【技术实现步骤摘要】
一种基于内源释放的大型浅水湖泊模型及其构建方法
本专利技术涉及水环境数值模拟领域，具体地涉及一种基于内源释放的大型浅水湖泊模型及其构建方法。
技术介绍
国内外对于河湖沉积物氮磷内源释放的基本研究方法通常为：基于野外观测或室内实验，研究沉积物的静/动态释放规律或某种影响因子对沉积物释放强度的影响，归纳和分析氮磷释放与水质及富营养化问题之间的关系，或估算氮磷释放通量。然而实验室环境与野外实际环境差别很大。实验室所获得的数据引入了较大的不确定性；难以动态确定沉积物氮磷的释放规律，过于简化的估算模型所计算出的内源负荷与实际存在较大的误差；沉积物再悬浮过程可能被不同程度地简化或忽略；此外，从沉积物中有机氮磷的衰变到向上覆水释放，这一过程的机理研究少于实验研究，包含机理的模型研究远少于估算或根据实验的回归结论。这些研究现状所存在的问题，对研究的手段和方法提出了更高的要求。计算有机物在沉积物种的衰减变化以及向上覆水释放的全过程的模型被称为沉积成岩模型。由于野外监测条件的客观限制、水体本身的特异性和模型参数不确定性等因素，沉积成岩模型在湖泊水环境模拟中的应用案例极少，沉积物营养盐释放通量，通常依然采用参数估计和参数调节的方法进行率定。。再悬浮沉积物在过去的几十年里受到越来越多的关注，它的研究是非常重要的，作为内源污染物的重要来源，沉积物再悬浮对水体环境的影响有时可以超过长期的外荷载。悬浮物会影响湖泊营养物的循环，间接导致了水体的富营养化，悬浮物也可以影响光衰减、有机污染物，甚至增加水体中有毒重金属的含量。因此，有必要了解沉积物再悬浮的动力机制。
技术实现思路
本专利技术为了提供一种可靠的耦合模型以研究沉积物再悬浮过程中营养盐释放的动力机制，提供了一种基于内源释放的大型浅水湖泊模型及其构建方法。本专利技术所采取的技术方案为：一种基于内源释放的大型浅水湖泊模型构建方法，包括以下步骤S01构建水动力模块：对水动力边界条件敏感性和水动力参数敏感性进行分析量化水动力边界条件和水动力参数的敏感性，据敏感性对参数进行调整优化；S02构建水质模块：根据大气沉降系数、湖泊边界出入流水质数据和湖泊内部的水质数据为所述水质模块提供水质边界和初始水质条件，在水质参数范围内进行抽样，量化水质参数敏感性，选定最佳水质参数取值；S03构建沉积物模块：包括设置黏性沉积物层和非黏性泥砂层，获取沉积物边界条件、初始条件和沉积物参数构建所述沉积物模块；S04构建沉积物成岩模块：对成岩参数进行抽样，利用区域敏感性分析方法量化成岩参数敏感性，根据结果选定最佳成岩参数取值；S05判定沉积物再悬浮的底部切应力条件，耦合沉积物和水质模块。进一步的，所述步骤S03中沉积物参数包括底部摩擦系数、垂向分层数、临界沉降系数、临界冲刷系数、表面侵蚀率、沉降速度以及泥床形态和固化特性。进一步的，所述步骤S02中湖泊边界出入流水质数据和湖泊内部的水质数据通过YSI水质监测仪或实验分析获取，所述水质监测仪获取数据包括DO、Chl-a、NH4+-N，所述实验分析获取数据包括TN、NH4+-N、NO3--N、NO2--N、TP、PO43-。进一步的，所述步骤S05具体包括：a模拟不同切应力条件下沉积物再悬浮情景，绘制切应力和水体中悬浮物的浓度的关系，得到相应的拟合函数；b构建边界层的沉积物通量，并分为再悬浮通量和沉降通量；c设置不同的营养盐吸附/解吸系数，将底泥营养盐释放与底泥起悬沉降过程联系起来，耦合沉积物模块和水质模块。进一步的，所述成岩参数包括沉降到沉积物中藻类的Gi分类份数、沉积物通量区域常数参数、衰减速率、温度影响系数、沉积物初始温度和温度扩散系数、底栖生物量的底部累积应力。进一步的，采用拉丁超立方抽样方法在水质参数范围内进行抽样，利用区域敏感性分析方法量化水质参数敏感性，选定最佳水质参数取值；采用拉丁超立方抽样方法在确定的成岩参数范围内进行抽样，利用区域敏感性分析方法量化成岩参数敏感性，选定最佳成岩参数取值。进一步的，在所述步骤S01前包括分区取样测定沉积物特性，具体包括：划分湖区、选定研究点、分层取水样和泥样、测定水样和泥样；所述划分湖区根据湖区大小、边界特征等将浅水湖泊划分不同的湖区，便于选定研究点；所述选定研究点主要在不同湖区内藻华爆发、水体黑臭、底泥淤积的区域选择典型的研究点，为了进一步地测定研究；所述分层取水样和泥样包括根据水深将水层分成2～3层进行分别取水样，泥样取出后按照取样深度分成4～5层；所述测定水样和泥样包括测定水样的性质和泥样的特性。本专利技术所产生的有益效果包括：本专利技术将模块构建和参数敏感性研究结合起来优化参数取值，为模块耦合打好基础，通过耦合后的模型研究营养盐释放；本专利技术结合野外观测、室内实验、模型构建过程，构建了一种运用数值模拟技术反演湖泊沉积物再悬浮营养盐释放过程的方法，克服了以往并没有相关方法的缺陷。本专利技术方法简单，实施方便，为科学合理构建的模型提供依据。附图说明图1本专利技术中方法的流程图；图2野外观测平台示意图；图3流量边界示意图；图4风场数据示意图；图5太湖分区示意图；图6分层泥样示意图；图7太湖模型网格划分图；图8太湖底泥分布示意图；图9切应力和水体中悬浮物的浓度关系图。具体实施方式下面结合具体实施例对本专利技术做进一步详细的解释说明，但应当理解为本专利技术的保护范围并不受具体实施例的限制。如图1，本专利技术中一种基于内源释放的大型浅水湖泊模型构建方法，包括如下步骤：(1)资料收集及实地调查监测；(2)分区取样测定沉积物特性；(3)构建水动力模块并进行边界条件和参数敏感性分析；(4)构建水质模块并进行参数敏感性分析；(5)构建沉积物模块；(6)构建沉积物成岩模块并进行参数敏感性分析；(7)判定沉积物再悬浮的底部切应力条件，耦合沉积物和水质模块。步骤(1)中资料收集包括查阅文献、利用地图软件、遥感技术等获取所研究区域水域范围及底泥分布的初步判定，实地调查监测包括咨询湖区周边气象站、水文站、水质监测站获取风场、湖区水位、出入湖流量、出入湖水质、地形等信息，搭建小型监测站搭建野外观测平台，利用垂向布置的多个自动风速风向仪和能够观测湖体表层流场的ADV(声学多普勒流速仪)或其他流速测量设备获取高频实时三维流速数据。步骤(2)中分区取样测定沉积物特性包括划分湖区、选定研究点、分层取水样和泥样、测定水样和泥样。划分湖区根据湖区大小、边界特征等将浅水湖泊划分不同的湖区，便于选定研究点。选定研究点主要在不同湖区内藻华爆发、水体黑臭、底泥淤积的区域选择典型的研究点，为了进一步地测定研究。分层取水样和泥样包括根据水深将水层分成2～3层进行分别取水样，泥样取出后垂直分成4～5层。测定水样和泥样包括测定水样的性质(pH值、电导率、氧化还原电位、温度、溶解氧浓度、悬浮物浓度、磷酸盐浓度、总磷浓度、氨氮浓度、总氮浓度、叶绿素a浓度等)和泥样的特性(泥厚、粒径、容重、密度、含水率等)，为模块构建提供边界条件、初始条件、参数设置、率定过程等需要的数据。步骤(3)中构建水动力模块并进行边界条件和参数敏感性分析包括利用地图软件、遥感技术等获取的湖泊岸线数据和地形数据进行网格划分，根据获取的气象数据为模型提供大气边界，根据获取的湖区周边的风场数据进行空间差值为模型提供风场边界，根据收集的出入本文档来自技高网...

【技术保护点】
1.一种基于内源释放的大型浅水湖泊模型的构建方法，其特征在于：包括以下步骤S01 构建水动力模块：对水动力边界条件敏感性和水动力参数敏感性进行分析量化水动力边界条件和水动力参数的敏感性；S02 构建水质模块：根据大气沉降系数、湖泊边界出入流水质数据和湖泊内部的水质数据为所述水质模块提供水质边界和初始水质条件，在水质参数范围内进行抽样，量化水质参数敏感性，选定最佳水质参数取值；S03 构建沉积物模块：包括设置黏性沉积物层和非黏性泥砂层，获取沉积物边界条件、初始条件和沉积物参数构建所述沉积物模块；S04 构建沉积物成岩模块：对成岩参数进行抽样，利用区域敏感性分析方法量化成岩参数敏感性，根据结果选定最佳成岩参数取值；S05 判定沉积物再悬浮的底部切应力条件，耦合沉积物和水质模块；步骤S03和步骤S04可置换顺序。

【技术特征摘要】
1.一种基于内源释放的大型浅水湖泊模型的构建方法，其特征在于：包括以下步骤S01构建水动力模块：对水动力边界条件敏感性和水动力参数敏感性进行分析量化水动力边界条件和水动力参数的敏感性；S02构建水质模块：根据大气沉降系数、湖泊边界出入流水质数据和湖泊内部的水质数据为所述水质模块提供水质边界和初始水质条件，在水质参数范围内进行抽样，量化水质参数敏感性，选定最佳水质参数取值；S03构建沉积物模块：包括设置黏性沉积物层和非黏性泥砂层，获取沉积物边界条件、初始条件和沉积物参数构建所述沉积物模块；S04构建沉积物成岩模块：对成岩参数进行抽样，利用区域敏感性分析方法量化成岩参数敏感性，根据结果选定最佳成岩参数取值；S05判定沉积物再悬浮的底部切应力条件，耦合沉积物和水质模块；步骤S03和步骤S04可置换顺序。2.根据权利要求1所述的基于内源释放的大型浅水湖泊模型的构建方法，其特征在于：所述步骤S03中沉积物参数包括底部摩擦系数、垂向分层数、临界沉降系数、临界冲刷系数、表面侵蚀率、沉降速度以及泥床形态和固化特性。3.根据权利要求1所述的基于内源释放的大型浅水湖泊模型的构建方法，其特征在于：所述步骤S02中湖泊边界出入流水质数据和湖泊内部的水质数据通过YSI水质监测仪或实验分析获取，所述水质监测仪获取数据包括DO、Chl-a、NH4+-N，所述实验分析获取数据包括TN、NH4+-N、NO3--N、NO2--N、TP、PO43。4.根据权利要求1所述的基于内源释放的大型浅水湖泊模型的构建方法，其特征在于：所述步骤S05具体包括：a模拟不同切应力条件下沉积物再...

【专利技术属性】
技术研发人员：李一平，施媛媛，程一鑫，程月，罗凡，徐芸蔚，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏,32