水体表面漂浮物迁移、堆积和清漂过程的数值模拟方法
【专利摘要】本发明提供一种水体表面漂浮物迁移、堆积和清漂过程的数值模拟方法，其特征在于，包括以下步骤：步骤1. 漂浮物迁移过程模拟：基于漂浮物集合F中各个漂浮物在水体运动过程中的受力情况，将每个漂浮物都分别概化为具有一个质心的实体，求取每个漂浮物实体的六个属性的变化的过程；步骤2. 漂浮物堆积过程模拟：将具有相同离散网格编号Qij的漂浮物合并为一个新的漂浮物，重新求取该新的漂浮物的质量、运动速度和质心坐标，并将新的漂浮物纳入漂浮物集合中；步骤3. 清漂过程模拟：根据漂浮物清理计划，将漂浮物进行剔除。本方法能够快速、准确地模拟漂浮物在水体中的迁移规律，从而能够提高海水体漂浮物的综合管理效率。
【专利说明】
水体表面漂浮物迁移、堆积和清漂过程的数值模拟方法
技术领域
[0001] 本发明属于水力学、环境管理和环境模拟领域，具体涉及一种水体表面漂浮物迀 移、堆积和清漂过程的数值模拟方法。 技术背景
[0002] 漂浮物是依赖于浮力漂浮在水体表面的固体废弃物，包括原木、树枝、稻草、桔杆、 塑料制品和人畜尸体等。漂浮物的存在，对水体水质，水面景观，供水，水产，航运等构成不 利影响，不仅破坏生态环境，还会威胁饮水安全。尤其是在具有发电功能的水库区域，漂浮 物顺流而下，易于聚集在河道凹岸、拦河坝前，不仅影响发电水头，减小水电枢纽的发电效 益，对枢纽运行安全构成威胁。漂浮物的存在对水体正常功能的发挥造成了极大的影响，开 展漂浮物的监测、清理和排除，是减少漂浮物危害的重要手段。而开展漂浮物的治理，需要 对水体中漂浮物的运动规律有全面的把握。漂浮物在水体表面的迀移模拟技术，能提供漂 浮物运动的实时信息，有助于漂浮物治理方案的构建。基于漂浮物模拟技术提供的实时信 息，在合适的位置和时间开展清漂工作，进行漂浮物监控，有助于快速有效的清理流域漂浮 物、降低排漂清漂的费用，节省人力物力。
[0003] 近年来，随着漂浮物对水质、水生态、水工建筑的不利影响逐渐受到重视，针对漂 浮物的研究不断发展。目前，主要在海洋区域针对漂浮物的相关研究较多，国外学者通过拉 格朗日和空间建模技术，构建了漂浮物预测模型，模型能够模拟海洋中漂浮物的运移路径 和轨迹，但对漂浮物的堆积以及质量缺乏量化和考虑，且缺少对河流、水库等水体漂浮物的 迀移运动过程的模拟研究。我国学者在对漂浮物的研究中，主要采用采样和实验分析方法， 包括研究漂浮物的丰度、粒径、空间分布及其组成成分，研究漂浮物对水工建筑的影响以及 漂浮物的治理措施等。关于漂浮物的定量研究较少，如通过分析漂浮物的水流特性，提出漂 浮物模拟的相似条件，应用正态模型或小变率的变态模型来模拟漂浮物的堆积形态;通过 假定漂浮物与水流质点运动的轨迹一致，基于平面二维水流数值模型用水流速度模拟漂浮 物运动，分析漂浮物的堆积问题。目前已有的研究在一定程度上为漂浮物数值模拟的进一 步开展提供了技术借鉴，但现有研究仍然不能快速、准确的模拟漂浮物从产生到堆积或清 理的过程。
[0004] 漂浮物在水体中受水和风的双重动力作用，单纯用水的流速代替漂浮物的速度进 行模拟，考虑不够全面，且水流数学模型采用离散化的网格形式，模拟的是网格的速度，格 点的速度与漂浮物质心的速度存在一定的差别；而在漂浮物的堆积问题中，清漂措施会极 大的影响漂浮物的总量和运移情况。因此，在现有研究的基础上，通过合理概化漂浮物受水 流、风速的影响，综合考虑漂浮物的入汇、堆积以及清漂过程，模拟漂浮物的运移过程。为了 更加有效的进行漂浮物治理，建立漂浮物模型，快速、准确地模拟漂浮物在水体中的迀移规 律，对大尺度流域的产漂量进行定量研究，对漂浮物的去除，提高海水体漂浮物的综合管理 有重要意义。
【发明内容】
[0005] 本发明是为了解决上述问题而进行的，目的在于提供一种能够快速、准确地模拟 漂浮物在水体中的迀移规律的体表面漂浮物迀移、堆积和清漂过程的数值模拟方法。
[0006] 本发明为了实现上述目的，采用了以下方案：
[0007] 本发明提供一种水体表面漂浮物迀移、堆积和清漂过程的数值模拟方法，其特征 在于，包括以下步骤：
[0008] 步骤1.漂浮物迀移过程模拟：
[0009] 基于漂浮物集合F中各个漂浮物在水体运动过程中的受力情况，将每个漂浮物都 分别概化为具有一个质心的实体，求取每个漂浮物实体的六个属性的变化的过程，这六个 属性分别为质量Qf，正东方向流速uf，正北方向流速vf，质心在正东和正北方向的大地坐标 Xf和yf，以及质心所在的离散网格编号Qij;
[0010] 步骤2.漂浮物堆积过程模拟:将具有相同离散网格编号Qu的漂浮物合并为一个新 的漂浮物，重新求取该新的漂浮物的质量、运动速度和质心坐标，并将新的漂浮物纳入漂浮 物集合中；
[0011] 步骤3.清漂过程模拟：
[0012] 根据水体的漂浮物清理计划，将处于清理计划所在的清漂时间段和清理区域的漂 浮物进行剔除:求取当前时刻下漂浮物的六个属性，对比漂浮物所在的网格是否处于清漂 区域，若是，则进一步对比当前时刻是否处于清漂时间段内，若是，则将漂浮物清理出漂浮 物集合F。
[0013] 进一步地，本发明提供的水体表面漂浮物迀移、堆积和清漂过程的数值模拟方法 中，步骤1.漂浮物迀移过程模拟，还包括以下子步骤：
[0014] 步骤1-1.漂浮物入汇条件给定:漂浮物入汇条件是漂浮物计算的前提，给定在时 亥|Jt所汇入的漂浮物的质量Qf 0，排放位置xf 0和yf 0，流速uf 0 = 0和vf 0 = 0，以及排放处的网格 编号Q_，设漂浮物实体为float，漂浮物组成的漂浮物集合为F={float| float为漂浮物实 体}，则每排入一个漂浮物，F中的漂浮物的数量增加1个；
[0015] 步骤1-2.水动力模型构建:基于二维规则网格构建技术，构建用于漂浮物模拟的 二维规则网格体系G[ij]及其水动力模型，该水动力模型基于二维圣维南方程，如下：
[0021]式中：h为水深，u为正东方向的流速，v为正北方向的流速，X为正东方向的投影坐 标，y为正北方向的投影坐标，t为时间，g为重力加速度，p为水体密度，c为谢才系数，f为柯 氏力常数，/ = 2flsin供，?>为炜度，Ω为地转角速度，约为2π/(24Χ3600)弧度每秒，ξχ、Cy为 分别为X、y方向上的涡动粘滞系数；；τ χ、τ y为分别为X、y方向上的风切应力， 其表达形式为：
[0024]式中：Ca为风阻力系数;pa为空气密度;Wx，W y为分别为x、y方向上的风速；
[0025]通过二维水动力模型，模拟任意时刻每个网格上的正东方向和正北方向的表面流 速，Uij 和 Vij;
[0026]步骤1-3.气象插值:构建二维规则网格内及其周边的气象站集合W，W中包括有每 个气象站的位置(xw，yw)，该气象站t时刻监测到的正东方向的风速wu和正北方向的风速 wv，通过反距离权重法，将集合W中所有气象站的风速插值到每个二维规则网格，得到t时刻 每个网格上的正东方向和正北方向的风速为WUij，WVi j;反距离权重插值的方程如下：
[0030] 式中，Μ为第Si个气象站所占权重，wu(Si)，wv( Si)为第Si个气象站实测的正东和正 北方向的风速值，di〇为网格ij与其第si个气象站之间的距离，N为气象站个数；
[0031] 步骤1-4.漂浮物运移速度计算:漂浮物在水体中受水和风的作用力而发生变化， 将漂浮物所受的力概化成流速和风速对漂浮物运动速度的影响，漂浮物质心所在位置的运 动速度uf和vf通过下式求取：
[0032] uf = cw X wub+Ca X ub,
[0033] vf = cw X wvb+Ca X vb,
[0034] 式中，cw为漂浮物运动的水阻系数，ca为漂浮物运动的风阻系数， WUb，WVb，Ub，vb* 别为漂浮物质心所在位置正东和正北方向的风速和流速;在所述步骤1-2和所述步骤1-3 中，求取到的速度为网格的流速(Uij和Vij)与风速(WUij和WVij) #1113，'\￥￥13，1113，￥13通过漂浮物 所在网格1及其相邻网格2~9的流速和风速，利用反距离权重法插值获得：
[0036]其中，i为与漂浮物相关的网格编号，1为漂浮物质心所在的网格，2,3,4,5分别为1 的西、东、南、北方向的网格，6,7为4号网格的西、东方向的网格、8,9为5号网格的西、东方向 的网格，cU沩漂浮物的质心(xf，yf)到网格质心ij(x，y)的距离；
[0037]步骤1-5.漂浮物位置计算:漂浮物在时刻t所处的位置为xf，yf，则其位置满足：
[0038] xft = xf(t-i)+ufX At
[0039] yft = yf(t-i)+vfX At,
[0040] 其中，xft，yft为漂浮物当前时刻t的位置，xft-i，yft-i为上一时刻的位置，Δ t为当 前时刻与上一时刻的时间差，等于一个计算时间步长;若t-i时刻为漂浮物排入的时刻， xft-i - xfo 7 yft-i - yfo ；
[0041] 步骤1-6,漂浮物所在网格计算：
[0042] 在漂浮物模拟计算中，一个计算时间步长内，漂浮物运移的速度可能超过一个网 格的范围，采用模式分离法，通过动步长模式，计算漂浮物所在的网格;设定当前漂浮物所 在的网格为i j，网格的X方向的长度为dx，y方向的长度为dy，模型模拟步长为dt，则用于计 算漂浮物运移的步长为dtf=max(dt/(dx/uf)，dt/(dy/vf));
[0043] 在一个漂浮物模拟步长后，漂浮物的移动距离为dxf和dyf，其当前的位置为xft = xft-i+dxf，yft = yft-i+dyf，设运动前漂浮物所在的网格号为i j (i，j)，在i j网格的东西南北 方向外扩至少1层网格进行搜索；
[0044] 若漂浮物质心所在的位置Q与网格ijN四个顶点(A，B，C，D)所组成的向量及网格四 条边的向量，满足关系式：
[0045] ABXAQ>0
[0046] BCXBQ>0
[0047] CDXCQ>0
[0048] DAXDQ>0，
[0049] 则Q在网格i jN(AB⑶)的范围内，且有漂浮物所在的网格编号Qij = i jN。
[0050] 进一步地，本发明提供的水体表面漂浮物迀移、堆积和清漂过程的数值模拟方法 中，步骤2.漂浮物堆积过程模拟，包括以下子步骤：
[0051] 步骤2-1.判断漂浮物集合中是否存在具有相同编号的漂浮物：
[0052]在完成一个步长dt的计算后，检索漂浮物集合F中所有漂浮物所在的网格编号，在 集合F中，从第1个漂浮物向第N个漂浮物依次扫描，这里N为F中，漂浮物的总数，以当前编号 为ne[l，N]的漂浮物的所在网格编号Qij(n)为判断依据，从编号为第N的漂浮物向第n+1个 漂浮物，编号nk依次减1进行扫描，判断Qi j (η)与Qi j (nk)的大小；
[0053]若不存在网格编号相同的漂浮物，则继续下一个步长的计算；
[0054]若存在具有相同编号的漂浮物，则转入步骤2-2;
[0055]步骤2-2.合并相同网格内的漂浮物：
[0056] 在漂浮物体系F中，当如」（1〇=〇〇_(111〇时，合并漂浮物11与1^为新的漂浮物，并将 该新的漂浮物的编号设为η，它的坐标与速度满足：
[0058] 新的漂浮物的质量为两个合并的漂浮物质量之和:〇以11)=〇以11)+〇以111〇;
[0059] 步骤2-3.对漂浮物集合F进行重新排列：
[0060] 在将漂浮物η和漂浮物nk合并后，重新排列漂浮物体系，漂浮物1~nk-ι保持序号 不变，漂浮物nk+Ι~N序号依次减1，得到合并后总数为N-1的漂浮物集合F;
[0061] 重复执行步骤2-1~步骤2-3,直到具有相同ij号的漂浮物都被合并并且新的漂浮 物都被重新排列完成为止。
[0062] 进一步地，本发明提供的水体表面漂浮物迀移、堆积和清漂过程的数值模拟方法 中，步骤3.清漂过程模拟，还包括以下子步骤：
[0063]步骤3-1.定义清漂方案：按照水体中清漂的规划方案，定义清漂的时间段为TC [七1，〖2]，清漂位置集为?0=作1_|?1_为清漂的网格编号};
[0064]步骤3-2.清理满足条件的漂浮物：
[0065]在模拟过程中，若满足条件：当前所模拟的时间t满足且当前漂浮物所 在的网格Qi j (n) e FC，则将编号为η的漂浮物清理出漂浮物集合F;
[0066] 并将漂浮物1~η保持不变，漂浮物η+1~Ν序号依次减1，得到清出漂浮物编号为η 后的新的漂浮物集合F;
[0067] 重复执行步骤3-1至步骤3-2,直到所有满足条件的漂浮物都被清理完为止。
[0068] 进一步地，本发明提供的水体表面漂浮物迀移、堆积和清漂过程的数值模拟方法 还可以具有这样的特征:在步骤1-6中，是在网格的东西南北方向外扩3层网格进行搜索。 [0069]发明的作用与效果
[0070] 本发明与现有技术相比，具有以下优点和效果：
[0071] (1)将漂浮物概化为具有6个属性的实体，漂浮物移动、堆积和清理的过程，即是对 漂浮物6个属性的定量求解的过程。本技术相比于现有漂浮物的相关技术，具有明显的创新 性，技术可行，计算速度快。
[0072] (2)基于水力学领域成熟的平面二维水动力模拟技术，以离散网格的形式获得水 体表面的流速。分析漂浮物在水体所受的水流作用力和风作用力，通过流速和风速的概化 模式，基于漂浮物的质心与网格体系的距离关系，利用反距离权重法，获得漂浮物的运移速 度。通过该方法的处理，避免了对漂浮物的形态，在水体中所受复杂的力学结构的解析，具 有更好的可行性。
[0073] (3)利用模式分裂法，通过构建模型模拟步长和漂浮物运移计算步长的方式，进行 漂浮物移动的模拟，确保漂浮物在移动的过程中，每次计算步长内其运移的距离不会超过 一个网格的边长，从而确保对漂浮物所在网格的搜索更准确和可靠。
[0074] (4)基于构建的漂浮物迀移模拟技术，能对任何水体中浮在水面的漂浮物进行模 拟和预测，确定漂浮物在任意时刻所处的位置以及水体中任意位置所含有的漂浮物质量， 分析漂浮物的运移轨迹，为漂浮物的监控、清理、处置提供数据支撑。
[0075] 综上所述，本发明的体表面漂浮物迀移、堆积和清漂过程的数值模拟方法，能够快 速、准确地模拟漂浮物在水体中的迀移规律，对大尺度流域的产漂量进行定量研究，对漂浮 物的去除，从而能够切实提高海水体漂浮物的综合管理效率。
【附图说明】
[0076] 图1为本发明实施例的模拟范围网格体系的示意图；
[0077] 图2为本发明实施例的漂浮物入汇位置和过程的示意图；
[0078] 图3为本发明实施例中通过气象站插值获得网格上的风速的示意图；
[0079] 图4为本发明实施例中漂浮物质心与周边9个网格的反距离权重插值的示意图；
[0080] 图5为本发明实施例中漂浮物运移后搜索其所在网格ij号的示意图；
[0081] 图6为本发明实施例中合并具有相同网格号的漂浮物的流程图；
[0082] 图7为本发明实施例中模拟2天后漂浮物运移状态的效果图。
【具体实施方式】
[0083] 以下结合附图对本发明涉及的水体表面漂浮物迀移、堆积和清漂过程的数值模拟 方法的具体实施方案进行详细地说明。
[0084] 实施例
[0085] 本实施例的水体表面漂浮物迀移、堆积和清漂过程的数值模拟方法包括以下步 骤：
[0086]步骤1.漂浮物迀移过程模拟
[0087] 1-1.漂浮物水动力模型构建:利用规则网格绘制软件，如delft3d等软件绘制要模 拟区域的二维规则网格体系G[ij]，如图1所示。基于二维圣维南方程，构建水动力模型，水 动力模型基于二维圣维南方程，如下：
[0088] 连续方程：
[0090] 动量方程：
[0093]式中：式中:h为水深;u为正东方向的流速;v为正北方向的流速，X为正东方向的投 影坐标，y为正北方向的投影坐标;g为重力加速度;P为水体密度;C为谢才系数;f为柯氏力 常数，/__:=_2?.Ηη鲁,P为炜度，Ω为地转角速度，约为2π/(24 X 3600)弧度每秒；ξχ、ξγ为分别 为X、y方向上的涡动粘滞系数
;τχ、"^为分别为X、y方向上的风切应力，其表 达形式为：
[0094] Tr^CapaWx{W； + W；f2 (4)
[0095] r, =Capjrv(fr；+fFl2)h2 (5)
[0096] 式中：Ca为风阻力系数;pa为空气密度;Wx，Wy为分别为x、y方向上的风速。通过二维 水动力模型，模拟任意时刻每个网格上的正东方向和正北方向的表面流速，uu和 Vlj。在本 实施例中，设定模拟时间为30天，模拟步长为360秒，模拟的第1天的流速如附图1所示。
[0097] 1-2.给定入汇条件:在漂浮物的计算中，漂浮物入汇条件是漂浮物计算的前提，在 本实施例中给定从时刻〇到第30天的漂浮物的入汇量。如图2中，为30天内在入汇点SI，S2处 的漂浮物的入汇位置及其入汇过程线。
[0098] 1-3.气象插值:构建二维规则网格内及其周边的气象站集合W，在本实施例中W包 括有2个气象站，分别W1和W2。其中，W1的位置为(x wl，ywl)，W2的位置为(xw2，yw2)。每个气象站 t时刻监测到的正东方向的风速(wu)和正北方向的风速(wv)。通过反距离权重法，将集合W 中所有气象站的风速插值到每个二维规则网格，得到t时刻每个网格上的正东方向和正北 方向的风速为WUij，WVij。反距离权重插值的方程如下：
[0101] 式中，Ai为第Si个气象站所占权重，Xij，yij为网格ij的x，y值;wu(Si)，wv(Si)为第Si 个气象站实测的正东和正北方向的风速值;dlQ为网格ij与其第Sl个气象站之间的距离。插 值获得的第1天的网格上的风速如图3所示。
[0102] 1-4.漂浮物运移速度计算:漂浮物在水体中受水和风的作用力而发生变化。将漂 浮物所受的力概化成流速和风速对漂浮物运动速度的影响。漂浮物质心所在位置的运动速 度uf，vf通过下式求取： u\ -c x wuh +c xuh
[0103] ' . '' . , (8) V1 =Γ xuvb +c, XW
[0104] 其中，cw为漂浮物运动的水阻系数，ca为漂浮物运动的风阻系数;在本例中，水阻系 数为〇. 8，风阻系数为0.2。wub，wvb，ub，vb分别为漂浮物质心所在位置正东和正北方向的风 速和流速。利用漂浮物所在的位置与周边9个网格ijl~ij9的距离，通过反距离权重法，获 得漂浮物质心到9个网格的权重（如图4所示），将9个网格的流速(Uij和vij)与风速(wuij和 wVlj)插值到漂浮物所在的位置，通过公式(8)求取到漂浮物所在位置的速度uf，wf。
[0105] 1-5.漂浮物位置计算:漂浮物在时刻t所处的位置为Xf，yf，则其位置满足： xf, =xfu_h+ufxAi
[0106] f f (9) = +vfxAt
[0107]其中，xft，yft为漂浮物当前时刻t的位置，xft-i，yft-i为上一时刻的位置，Δ t为当 前时刻与上一时刻的时间差，等于一个计算时间步长。若t-Ι时刻为漂浮物排入的时刻，则 xft-i - xfo 7 yft-i - yfo 〇
[0108] 1-6.漂浮物所在网格计算:在漂浮物模拟计算中，一个时间步长60秒内，漂浮物运 移的速度可能超过一个网格的范围，采用模式分离法，通过动步长模式，计算漂浮物所在的 网格。设定当前漂浮物所在的网格为i j，网格的X方向的长度为dx，y方向的长度为dy，模型 模拟步长为dt，贝lj用于计算漂浮物运移的步长为dtf=max(dt/(dx/uf)，(11:/((17八;〇)。在一 个漂浮物模拟步长后，漂浮物的移动距离为dxf和dyf，其当前的位置为xf t = xft-i+dxf^yft =yft-i+dyf。设运动前漂浮物所在的网格号为i j (i，j)，在i j网格的东西南北方向外扩3层 网格，g卩[i-3，i+3]与[j-3，j+3]的范围内搜索。若漂浮物质心所在的位置Q与网格ij N四个 顶点(A，B，C，D)所组成的向量及网格四条边的向量，满足关系式：
[0110] 贝1JQ在网格i jN( ABCD)的范围内，且有漂浮物所在的网格编号Qi j = i jN。漂浮物由i j 处运移到Q处过程中，网格搜索的范围和确定漂浮物所在网格的ij如图5所示。
[0111] 步骤2.漂浮物堆积过程模拟
[0112] 2-1.漂浮物集合中，判断所有漂浮物所在网格:在完成一个步长dt的计算后，检索 漂浮物集合F中，所有漂浮物所在的网格编号。在集合F中，从第1个漂浮物向第N个漂浮物依 次扫描(N为F中，漂浮物的总数），以当前编号为ne[l，N]的漂浮物的所在网格编号Qij(n) 为判断依据，从编号为第N的漂浮物向第n+1个漂浮物，编号nk依次减1进行扫描，判断Qij (η)与Qij(nk)的大小。若不存在若不存在网格编号相同的漂浮物，则继续下一个步长的计 算。若存在具有相同编号的漂浮物，则转入步骤2-2。
[0113] 2-2.合并相同网格内的漂浮物：当漂浮物体系F中，Qi j (n) = Qi j (nk)时，合并漂浮 物n与nk为新的漂浮物，其编号为η，其坐标与速度满足：
[0115] 漂浮物质量为两个合并的漂浮物质量之和，Qf(n)=Qf(n)+Qf(nk)
[0116] 2-3.漂浮物体系F重新排列：在将漂浮物η和漂浮物nk合并后，重新排列漂浮物体 系，漂浮物1~nk-Ι保持序号不变，漂浮物nk+Ι~N序号依次减1，得到合并后总数为N-1的漂 浮物集合F。重复执行步骤2-1~步骤2-3,直到具有所有相同ij号的漂浮物合并完成为止。 漂浮物堆积的计算流程如图6所示。
[0117]步骤3.漂浮物清理过程模拟
[0118] 3-1.定义清漂方案:按照水体中清漂的规划方案，定义清漂的时间段为TC[tl，t2] 和清漂位置集FC = {Fi j | Fi j为清漂的网格编号}。
[0119] 3-2.清理满足条件的漂浮物:在模拟过程中，若当前所模拟的时间t满足 t2,且当前漂浮物所在的网格Qij(n)eFC，则将编号为η的漂浮物清理出F。将漂浮物1~η保 持不变，漂浮物η+1~Ν序号依次减1，则得到清出漂浮物编号为η后的漂浮物体系F。重复执 行，直到所有满足条件的漂浮物被清理完为止。经过计算，模拟的2天后计算范围内漂浮物 的状态如图7所示，从图中可以看出，漂浮物的移动轨迹与水流流场的轨迹相似。
[0120] 实施例的作用与效果：
[0121] (1)将漂浮物概化为具有6个属性的实体，漂浮物移动、堆积和清理的过程，即是对 漂浮物6个属性的定量求解的过程，技术可行，计算速度快。
[0122] (2)基于水力学领域成熟的平面二维水动力模拟技术，以离散网格的形式获得水 体表面的流速。分析漂浮物在水体所受的水流作用力和风作用力，通过流速和风速的概化 模式，基于漂浮物的质心与网格体系的距离关系，利用反距离权重法，获得漂浮物的运移速 度。通过该方法的处理，避免了对漂浮物的形态，在水体中所受复杂的力学结构的解析，具 有更好的可行性。
[0123] (3)利用模式分裂法，通过构建模型模拟步长和漂浮物运移计算步长的方式，进行 漂浮物移动的模拟，确保漂浮物在移动的过程中，每次计算步长内其运移的距离不会超过 一个网格的边长，从而确保对漂浮物所在网格的搜索更准确和可靠。
[0124] (4)基于构建的漂浮物迀移模拟技术，能对任何水体中浮在水面的漂浮物进行模 拟和预测，确定漂浮物在任意时刻所处的位置以及水体中任意位置所含有的漂浮物质量， 分析漂浮物的运移轨迹，为漂浮物的监控、清理、处置提供数据支撑。
[0125]以上实施例仅仅是对本发明技术方案所做的举例说明。本发明所涉及的水体表面 漂浮物迀移、堆积和清漂过程的数值模拟方法并不仅仅限定于在以上实施例中所描述的内 容，而是以权利要求所限定的范围为准。本发明所属领域技术人员在该实施例的基础上所 做的任何修改或补充或等效替换，都在本发明的权利要求所要求保护的范围内。
【主权项】
1. 一种水体表面漂浮物迁移、堆积和清漂过程的数值模拟方法，其特征在于，包括w下 步骤： 步骤1.漂浮物迁移过程模拟： 基于漂浮物集合F中各个漂浮物在水体运动过程中的受力情况，将每个漂浮物都分别 概化为具有一个质屯、的实体，求取每个漂浮物实体的六个属性的变化的过程，运六个属性 分别为质量Qf，正东方向流速uf，正北方向流速vf，质屯、在正东和正北方向的大地坐标xf和 yf，W及质屯、所在的离散网格编号化j; 步骤2.漂浮物堆积过程模拟： 将具有相同离散网格编号化北勺漂浮物合并为一个新的漂浮物，重新求取该新的漂浮物 的质量、运动速度和质屯、坐标，并将新的漂浮物纳入漂浮物集合中； 步骤3.清漂过程模拟： 根据水体的漂浮物清理计划，将处于清理计划所在的清漂时间段和清理区域的漂浮物 进行剔除:求取当前时刻下漂浮物的六个属性，对比漂浮物所在的网格是否处于清漂区域， 若是，则进一步对比当前时刻是否处于清漂时间段内，若是，则将漂浮物清理出漂浮物集合 F。2. 根据权利要求1所述的水体表面漂浮物迁移、堆积和清漂过程的数值模拟方法，其特 征在于： 步骤1.漂浮物迁移过程模拟，包括W下子步骤： 步骤1-1.漂浮物入汇条件给定:给定在时刻t所汇入的漂浮物的质量Qfo,排放位置xfo 和yfo,流速uf〇 = 0和vf日=0, W及排放处的网格编号化j日，设漂浮物实体为f loat，漂浮物组 成的漂浮物集合为F= {float I float为漂浮物实体}，则每排入一个漂浮物，F中的漂浮物的 数量增加1个； 步骤1-2.水动力模型构建:基于二维规则网格构建技术，构建用于漂浮物模拟的二维 规则网格体系G[ij]及其水动力模型，该水动力模型基于二维圣维南方程，如下： 连续方程：式中：h为水深，U为正东方向的流速，V为正北方向的流速，X为正东方向的投影坐标，y 为正北方向的投影坐标，g为重力加速度，t为时间，P为水体密度，C为谢才系数，f为柯氏力 常数，/ = 2nsinp，口为缔度，Ω为地转角速度，约为2V(24X3600)弧度每秒，ξχ、ξy为分别 为X、y方向上的满动粘滞系数：Tx、Ty为分别为x、y方向上的风切应力，其表 达形式为：式中：Ca为风阻力系数;Pa为空气密度;Wx，Wy为分别为X、y方向上的风速； 通过二维水动力模型，模拟任意时刻每个网格上的正东方向和正北方向的表面流速， IHj和 VU; 步骤1-3.气象插值:构建二维规则网格内及其周边的气象站集合W，W中包括有每个气 象站的位置(xw，yw)，该气象站t时刻监测到的正东方向的风速WU和正北方向的风速WV，通 过反距离权重法，将集合W中所有气象站的风速插值到每个二维规则网格，得到t时刻每个 网格上的正东方向和正北方向的风速为WUij，WVij，反距离权重插值的方程如下：式中，λ功第Si个气象站所占权重，WU(Si)，WV(Si)为第Si个气象站实测的正东和正北方 向的风速值，dio为网格ij与其第Si个气象站之间的距离，N为气象站个数； 步骤1-4.漂浮物运移速度计算:漂浮物在水体中受水和风的作用力而发生变化，将漂 浮物所受的力概化成流速和风速对漂浮物运动速度的影响，漂浮物质屯、所在位置的运动速 度uf和vf通过下式求取： uf = Cw X WUb+Ca X ub, vf = Cw X Wvb+Ca X vb， 式中，Cw为漂浮物运动的水阻系数，Ca为漂浮物运动的风阻系数，^16，群13，油，乂6分别为 漂浮物质屯、所在位置正东和正北方向的风速和流速;在所述步骤1-2和所述步骤1-3中，求 取到的速度为网格的流速(Uij和Vij)与风速(WUij和WVij) ;*址)，群13，1113，￥13通过漂浮物所在网 格1及其相邻网格2~9的流速和风速，利用反距离权重法插值获得：其中，i为与漂浮物相关的网格编号，1为漂浮物质屯、所在的网格，2,3,4,5分别为1的 西、东、南、北方向的网格，6，7为4号网格的西、东方向的网格、8，9为5号网格的西、东方向的 网格，di功漂浮物的质屯、(xf，yf倒网格质屯、U (X，y)的距离； 步骤1-5.漂浮物位置计算:漂浮物在时刻t所处的位置为xf，yf，则其位置满足： xft = xf(t-i)+ufX At yft = yf(t-i)+vf X Δ t， 其中，为漂浮物当前时刻t的位置，为上一时刻的位置，Δ t为当前时 刻与上一时刻的时间差，等于一个计算时间步长，若t-1时刻为漂浮物排入的时刻，xft-i = xfo，yft-i = yfo; 步骤1-6，漂浮物所在网格计算： 在漂浮物模拟计算中，一个计算时间步长内，漂浮物运移的速度可能超过一个网格的 范围，采用模式分离法，通过动步长模式，计算漂浮物所在的网格;设定当前漂浮物所在的 网格为ij，网格的X方向的长度为dx，y方向的长度为dy，模型模拟步长为dt，则用于计算漂 浮物运移的步长为化f=max(化/(dx/uf)，化/(dy/vf)); 在一个漂浮物模拟步长后，漂浮物的移动距离为dxf和dyf，其当前的位置为xft = xft-i+ dxf，yf t =州*-1+(1州，设运动前漂浮物所在的网格号为ij(ij)，在ij网格的东西南北方向 外扩至少1层网格进行捜索； 若漂浮物质屯、所在的位置Q与网格ijN四个顶点(A，B，C，D)所组成的向量及网格四条边 的向量，满足关系式： ABXAQ>0 BCXBQ>0 CDXCQ>0 DAXDQ>0， 贝化在网格Un (ABCD)的范围内，且有漂浮物所在的网格编号Qi j = Un。3.根据权利要求1所述的水体表面漂浮物迁移、堆积和清漂过程的数值模拟方法，其特 征在于： 步骤2.漂浮物堆积过程模拟，包括W下子步骤： 步骤2-1.判断漂浮物集合中是否存在具有相同编号的漂浮物： 在完成一个步长dt的计算后，检索漂浮物集合F中所有漂浮物所在的网格编号，在集合 F中，从第1个漂浮物向第N个漂浮物依次扫描，运里N为F中，漂浮物的总数，W当前编号为η e[l，N]的漂浮物的所在网格编号Qij(n)为判断依据，从编号为第Ν的漂浮物向第η+1个漂 浮物，编号nk依次减1进行扫描，判断Qi j (η)与Qi j (nk)的大小； 若不存在网格编号相同的漂浮物，则继续下一个步长的计算； 若存在具有相同编号的漂浮物，则转入步骤2-2; 步骤2-2.合并相同网格内的漂浮物： 在漂浮物体系F中，当Qij(n)=Qij(nk)时，合并漂浮物η与nk为新的漂浮物，并将该新 的漂浮物的编号设为η，它的坐标与速度满足：新的漂浮物的质量为两个合并的漂浮物质量之和:Qf(n)=Qf(n)+Qf(nk); 步骤2-3.对漂浮物集合F进行重新排列： 在将漂浮物η和漂浮物nk合并后，重新排列漂浮物体系，漂浮物1~nk-1保持序号不变， 漂浮物nk+1~N序号依次减1，得到合并后总数为N-1的漂浮物集合F; 重复执行步骤2-1~步骤2-3,直到具有相同ij号的漂浮物都被合并并且新的漂浮物都 被重新排列完成为止。4. 根据权利要求1所述的水体表面漂浮物迁移、堆积和清漂过程的数值模拟方法，其特 征在于： 步骤3.清漂过程模拟，包括W下子步骤： 步骤3-1.定义清漂方案：按照水体中清漂的规划方案，定义清漂的时间段为TC[tl， *2]，清漂位置集为。〔={。。'巾。'为清漂的网格编号}; 步骤3-2.清理满足条件的漂浮物： 在模拟过程中，若满足条件：当前所模拟的时间t满足且当前漂浮物所在的 网格Qij(n)eFC，则将编号为η的漂浮物清理出漂浮物集合F; 并将漂浮物1~η保持不变，漂浮物η+1~Ν序号依次减1，得到清出漂浮物编号为η后的 新的漂浮物集合F; 重复执行步骤3-1至步骤3-2,直到所有满足条件的漂浮物都被清理完为止。5. 根据权利要求2所述的水体表面漂浮物迁移、堆积和清漂过程的数值模拟方法，其特 征在于： 其中，在步骤1-6中，是在网格的东西南北方向外扩3层网格进行捜索。
【文档编号】G06F17/50GK106096212SQ201610570763
【公开日】2016年11月9日
【申请日】2016年7月19日
【发明人】张万顺, 王永桂, 彭虹
【申请人】武汉大学