一种基于生物絮凝的水沙动力学模拟方法技术

本发明专利技术提供了一种基于生物絮凝的水沙动力学模拟方法，所述方法包括如下步骤：A)给定初始河道断面资料，利用上游边界条件、下游边界条件及河道糙率计算出各个断面的水力要素，进而得到整个河道的水力要素分布；B)利用所述水力要素计算出水体紊动剪切速率，然后利用所述水体紊动剪切速率和测定的水质参数计算出悬移泥沙生物絮凝体的特征数据；C)利用上游来沙过程线、所述各个断面的水力要素以及所述悬移泥沙生物絮凝体的特征计算出相应断面泥沙生物絮凝体的冲刷淤积及运动输移，给出冲淤之后的河道断面情况；D)在相应的时间节点上，重复步骤A)-C)直到完成整个时段的计算，模拟河道冲淤的全过程，给出河道的演变规律。

【技术实现步骤摘要】
一种基于生物絮凝的水沙动力学模拟方法
本专利技术涉及水利工程
，特别涉及生物絮凝的水沙动力学领域。
技术介绍
我国河道含沙量高、输沙量大，泥沙输移对河床演变、水利工程建设等均产生重要影响，因此泥沙研究在我国一直受到高度重视，并已形成一套相对成熟的研究体系。但传统的泥沙研究，以干净泥沙为研究对象，干净泥沙主要由岩石风化产生，其表面特性、颗粒特性和群体特性均以岩石颗粒为基础，忽略了微生物等的影响。近年来随着社会经济的快速发展，生产和生活用水的大量排放，点源和面源污染导致河道水体污染日趋严重、营养化程度增高，河流中的干净泥沙颗粒作为水环境中的重要载体，吸附了大量污染物，为微生物在其表面的吸附成膜提供了更为有利的条件，从而进一步加剧了水环境中泥沙颗粒、生物膜生长及污染物吸附的相互影响问题。这一问题传统的泥沙研究无法回答。目前已有的水沙动力学模型，只能模拟干净泥沙颗粒的絮凝及输移，以及由此所导致的河床冲淤演变及流量、水位、糙率等水沙动力特性变化。微生物是自然界中个体最小、数目最多、分布范围最广的生命体，在生长环境中大多是以附着状态而非游离状态存在，其通过代谢活动分泌称之为生物膜的胞外聚合物质(EPS)粘附于固体基底表面。在河道、湖泊水体中，有大量的微生物。观测发现，微生态系统对河流、湖泊、潮间带地区等水沙环境具有重要作用，微生物通过新陈代谢活动分泌生物膜等作用于泥沙颗粒，并形成新的泥沙生物絮凝体，又称之为“生物絮凝泥沙”或“生物絮凝体”。与传统模型研究中的普通泥沙颗粒相比，生物絮凝泥沙的理化性质均发生很大变化，生物絮凝泥沙的输移规律也发生改变。微生物吸附及生物膜的形成会强烈地改变泥沙颗粒本身的物理化学特征和运动输移规律等，并且导致泥沙颗粒的新的环境效应问题。但是对于生物絮凝影响下的水沙特性，由于其作用机理不同，影响因子由于侧重于具有生物活性的微生态系统而更显复杂化，单纯的物理化学理论已无法准确全面的描述其实际情况。因此，为了更加准确全面的描述泥沙生物絮凝体的输移状况，亟需建立一种耦合泥沙颗粒和微生态系统之间相关关系的模型系统。
技术实现思路
本专利技术主要针对实际工程需要，建立了一种基于泥沙生物絮凝的水沙动力学模型系统，该系统涵盖水位流量预报、河道演进、泥沙生物絮凝变化等模拟功能，具有功能齐全，方便实用的特点。据此，本专利技术提供了一种基于生物絮凝的水沙动力学模拟方法，所述方法包括如下步骤：A)给定初始河道断面资料，利用上游边界条件、下游边界条件及河道糙率计算出各个断面的水力要素，进而得到整个河道的水力要素分布；B)利用所述水力要素计算出水体紊动剪切速率，然后利用所述水体紊动剪切速率和测定的水质参数计算出悬移泥沙生物絮凝体的特征数据；C)利用上游来沙过程线、所述各个断面的水力要素以及所述悬移泥沙生物絮凝体的特征计算相应断面泥沙生物絮凝体的冲刷淤积及运动输移，给出冲淤之后的河道断面情况；D)在相应的时间节点上，重复步骤A)-C)直到完成整个时段的计算，模拟河道冲淤的全过程，给出河道的演变规律。其中，步骤A)对应于水动力学模块，步骤B)对应于泥沙生物絮凝模块，步骤C)对应于泥沙动力学模块。由于传统的水沙动力学模型系统忽略生物絮凝对水沙动力学的影响，因而不能真实的反应河道中实际情况。然而，专利技术人通过研究发现，生物絮凝对水沙动力学具有一定的影响，其不可作为一个忽略因素来用于实际的水沙动力学评价当中。因此，将生物絮凝的因素考虑到水沙动力学系统当中，无疑使河道相关测报的数据更加科学、全面，从而更精准的指导实际现场。在一个具体实施方式中，所述水力要素包括流量、流速、水位和水深。在一个具体实施方式中，所述水质参数可以包括水体营养浓度、生物膜最大指定生长速率、营养物质的半饱和浓度、悬移泥沙生物絮凝体携带生物膜的能力系数、生物絮凝计算系数和破坏计算系数。在一个具体实施方式中，所述悬移泥沙生物絮凝体的特征可以包括悬移泥沙生物絮凝体的粒径；优选所述悬移泥沙生物絮凝体的特征数据还包括悬移泥沙生物絮凝体的密度；特别优选所述悬移泥沙生物絮凝体的特征数据还包括悬移泥沙生物絮凝体的组成和/或其结构。其中，悬移泥沙生物絮凝体的组成主要包括悬移泥沙和附着于悬移泥沙上的生物膜两部分。附着于悬移泥沙上的生物膜可以简称为悬移生物膜；悬移泥沙生物絮凝体的结构主要指悬移泥沙生物絮凝体的孔隙率等。在一个优选的实施方式中，所述悬移泥沙生物絮凝体的特征包括悬移泥沙生物絮凝体的粒径和密度。为了更精准地将生物絮凝应用于泥沙动力学模拟方法中，并且为了使具体的操作更具可行性，在本专利技术的一个实施例中，将所述悬移泥沙生物絮凝体分为悬移泥沙和悬移生物膜两部分，并引入分形维数理论。由此可根据絮体的分形几何特征，结合悬移泥沙生物絮凝体的体积来推导求出其粒径。在一个具体实施方式中，根据方程求出所述悬移泥沙生物絮凝体的粒径L，其中，V为所述悬移泥沙生物絮凝体的总体积，d为所述悬移泥沙生物絮凝体的分形维数，LP为所述悬移泥沙生物絮凝体中的所述悬移泥沙的泥沙颗粒的粒径，是所述悬移泥沙生物絮凝体中的所述悬移泥沙的泥沙颗粒的体积；以及任选地，ρf＝ζρB+(1-ζ)ρS，其中，ρf为悬移泥沙生物絮凝体的密度，ρB为所述悬移泥沙生物絮凝体中的所述悬移生物膜的密度，ρS为所述悬移泥沙生物絮凝体中的所述悬移泥沙的密度。在一个具体实施方式中，根据方程V＝VS+VB＝(1-ζ)V+ζV或求出所述悬移泥沙生物絮凝体的总体积V；其中，VS为所述悬移泥沙生物絮凝体中的所述悬移泥沙的总体积，VB为所述悬移泥沙生物絮凝体中的所述悬移生物膜的总体积，ζ表示所述悬移泥沙生物絮凝体中的所述悬移生物膜的体积含量，且ζ＝VB/V，(1-ζ)表示所述悬移泥沙生物絮凝体中的所述悬移泥沙的体积含量，t为时间。在此求出的V可用于方程在一个具体实施方式中，用于方程中的通过方程求得，其中，Cs是悬移泥沙生物絮凝体中的所述悬移泥沙的浓度，且CS＝(1-Ω)C，其中，Ω是总悬移泥沙生物絮凝体中的所述悬移生物膜的质量含量，C是悬移泥沙生物絮凝体的浓度；(1-ζ)是悬移泥沙生物絮凝体中的所述悬移泥沙的体积含量；G为水体紊动剪切速率；ka为悬移泥沙生物絮凝体的絮凝速率，kb为悬移泥沙生物絮凝体的破坏速率，且ka和kb计算方程分别为和其中，k'a为无量纲悬移泥沙生物絮凝体的絮凝计算系数，kb'为无量纲悬移泥沙生物絮凝体的破坏计算系数，μ为水的动力粘度系数，Fy为悬移泥沙生物絮凝体的强度。在此求出的可用于方程在一个具体实施方式中，用于方程中的通过方程求得，其中，CB是所述悬移泥沙生物絮凝体中的所述悬移生物膜的浓度，且CB＝ΩC；η为所述悬移泥沙生物絮凝体中的所述悬移生物膜的生长速率，且ηmax为所述悬移泥沙生物絮凝体中的所述悬移生物膜最大指定生长速率，N为水体的营养浓度，Km为营养物质的半饱和浓度；K为悬移泥沙生物絮凝体携带生物膜的能力，K＝βVp＝β(L3-V)，β为悬移泥沙生物絮凝体携带生物膜的能力系数，Vp为悬移泥沙生物絮凝体中的孔隙体积。在本专利技术中，最后得到的是河道冲淤的全过程，这其中包括河道的河床高程、河道的断面形状、河道断面的平均挟沙量、河道断面的平均含沙量以及泥沙的沉降速度等。在一个优选的实施方式中，所述河道冲淤全过程包括河道的河本文档来自技高网...

【技术保护点】
一种基于生物絮凝的水沙动力学模拟方法，所述方法包括如下步骤：A)给定初始河道断面资料，利用上游边界条件、下游边界条件及河道糙率计算出各个断面的水力要素，进而得到整个河道的水力要素分布；B)利用所述水力要素计算出水体紊动剪切速率，然后利用所述水体紊动剪切速率和测定的水质参数计算出悬移泥沙生物絮凝体的特征数据；C)利用上游来沙过程线、所述各个断面的水力要素以及所述悬移泥沙生物絮凝体的特征计算出相应断面泥沙生物絮凝体的冲刷淤积及运动输移，给出冲淤之后的河道断面情况；D)在相应的时间节点上，重复步骤A)‑C)直到完成整个时段的计算，模拟河道冲淤的全过程，给出河道的演变规律。

【技术特征摘要】
1.一种基于生物絮凝的水沙动力学模拟方法，所述方法包括如下步骤：A)给定初始河道断面资料，利用上游边界条件、下游边界条件及河道糙率计算出各个断面的水力要素，进而得到整个河道的水力要素分布；B)利用所述水力要素计算出水体紊动剪切速率，然后利用所述水体紊动剪切速率和测定的水质参数计算出悬移泥沙生物絮凝体的特征数据；C)利用上游来沙过程线、所述各个断面的水力要素以及所述悬移泥沙生物絮凝体的特征计算出相应断面泥沙生物絮凝体的冲刷淤积及运动输移，给出冲淤之后的河道断面情况；D)在相应的时间节点上，重复步骤A)-C)直到完成整个时段的计算，模拟河道冲淤的全过程，给出河道的演变规律。2.根据权利要求1所述的方法，其特征在于，所述水力要素包括流量、流速、水位和水深。3.根据权利要求1所述的方法，其特征在于，所述水质参数包括水体营养浓度、生物膜最大指定生长速率、营养物质的半饱和浓度、悬移泥沙生物絮凝体携带生物膜的能力系数、生物絮凝计算系数和破坏计算系数。4.根据权利要求1所述的方法，其特征在于，所述悬移泥沙生物絮凝体的特征包括悬移泥沙生物絮凝体的粒径。5.根据权利要求4所述的方法，其特征在于，所述悬移泥沙生物絮凝体的特征数据还包括悬移泥沙生物絮凝体的密度。6.根据权利要求5所述的方法，其特征在于，所述悬移泥沙生物絮凝体的特征数据还包括悬移泥沙生物絮凝体的组成和/或其结构。7.根据权利要求1-6任意一项所述的方法，其特征在于，将所述悬移泥沙生物絮凝体分为悬移泥沙和悬移生物膜两部分，并引入分形维数理论。8.根据权利要求7所述的方法，其特征在于，根据方程求出所述悬移泥沙生物絮凝体的粒径L，其中，V为所述悬移泥沙生物絮凝体的总体积，d为所述悬移泥沙生物絮凝体的分形维数，LP为所述悬移泥沙生物絮凝体中的所述悬移泥沙的泥沙颗粒的粒径，是所述悬移泥沙生物絮凝体中的所述悬移泥沙的泥沙颗粒的体积。9.根据权利要求8所述的方法，其特征在于，根据方程V＝VS+VB＝(1-ζ)V+ζV或求出所述悬移泥沙生物絮凝体的总体积V；其中，VS为所述悬移泥沙生物絮凝体中的所述悬移泥沙的总体积，...

【专利技术属性】
技术研发人员：赵慧明，曹文洪，汤立群，王崇浩，王玉海，刘大滨，郭传胜，秦俊桃，
申请(专利权)人：中国水利水电科学研究院，
类型：发明
国别省市：北京;11