海洋污染物运动轨迹预测及滨海电厂取水安全预警方法与流程

1.本发明属于滨海电厂保障技术领域，尤其涉及一种海洋污染物运动轨迹预测及滨海电厂取水安全预警方法。


背景技术：

2.近年海洋生物或异物影响滨海电厂厂取水安全的事件国内外曾多次发生，且有增多趋势，冷源安全已成为影响滨海电厂厂安全的重要因素。根据世界核营运者协会(wano)分析，大约20％此类事件对安全相关系统有直接影响。海生物大量聚集、污染物突然涌入等极端情况下，导致拦污结构的破坏、机组停堆等事故严重影响滨海电厂的取水安全。
3.由于污染物种类繁多，且囿于科学的预测手段，目前应对污染物的主要方式为被动清理的方式，即当在潮流及冷源取水作用下，污染物堆积至取水口门的各道拦污网之上，再由人工进行打捞。对于日常运营时期，经常性维护能够保证电厂的安全运行。但对于极端状况下，大量污染物突然涌入取水口门后，由于短时间内人力物力不能满足污染物清理要求，导致污染物大量聚集致使拦污网结构破坏或严重堵塞拦污网而影响取水，从而进一步引起电厂的停堆，影响电厂的运营效率。
4.海上污染物的运动受潮流、风力及冷源取水作用综合影响，对于特定地区，潮流具有周期性的特征，因此在掌握地区潮流及取水特性情况下，可对污染物的运动轨迹进行分析及预测，由此结合现场的监测，可计算污染物到达取水口门的时间，为滨海电厂运营单位提供有效的预警信息，从而为滨海电厂的防灾减灾工作提供有效的技术支撑。


技术实现要素：

5.针对现有技术存在的问题，本发明提供了一种海洋污染物运动轨迹预测及滨海电厂取水安全预警方法。
6.本发明是这样实现的，一种海洋污染物运动轨迹预测及滨海电厂取水安全预警方法，其特征在于：包括如下步骤：
7.(1)建立潮流数学模型：首先根据工程区域的特征，建立大范围潮流数学模型，外海边界距离工程位置50～100km；数学模型基于流体不可压缩性及navier～stokes方程，并服从boussinesq假定和静水压力假定；利用潮汐预报软件chinatide提供潮位边界资料，通过对水文实测资料进行验证，保证该工程海域下数学模型的各项参数有效性，确保潮流数学模型的准确性；
8.(2)建立粒子追踪模型：将呈现离散状、无自游泳能力或弱自游泳能力进入明渠的漂浮类污染物在数学模型中将其概化为粒子；粒子追踪模型采用langevin方程来描述粒子的迁移运动；
9.(3)确定粒子源释放点：以工程取水口为中心，将其周围海域根据距离及方位分别设置若干个粒子源释放点；
10.(4)统计粒子的通量及分布情况：根据潮流数学模型和粒子追踪模型计算不同点
在不同时刻、不同潮型释放污染物后，粒子随潮流及明渠取水作用下的运动轨迹，并统计粒子通过取水口拦污网位置的通量及分布情况；此时完成污染物运动轨迹预测；
11.(5)制作概率云图：根据统计结果，制作污染物进入明渠的概率云图，并评估污染物对取水安全的影响，通过概率云图可判定污染物来向的危险程度，重点加强对该方位的监测，通过提供观测设备的性能，监测范围及距离越远，预留应对措施的时间；
12.(6)实时监测及标记：根据现场对污染物开展实时监测，对已进入工程区域观测范围的污染物进行及时标记，掌握其距取水口位置、方位及时间；
13.(7)预测已标记粒子的后续运动轨迹：将已标记粒子信息输入到粒子数学模型中，通过粒子数学模型计算，预测已标记粒子的后续运动轨迹并计算至取水口的时间，为工程运营单位提出预警，以及时采取应急措施，包括提前调配人员及设备以提高清理能力，或根据污染物的类型选择不同的打捞方案。
14.优选的，步骤1中的实测水文资料主要包括潮位、流速、流向资料。
15.优选的，步骤2中的无自游泳能力或弱自游泳能力的物体包括泡沫、秸秆、树枝、塑料及水草中的一种或多种。
16.优选的，步骤3中依工程取水口为中心，以距离10km范围内、每间隔1km设置一个粒子释放点，且以22.5
°
为一个方向。
17.优选的，步骤4中计算不同点位在涨急时刻、涨停时刻、落急时刻、落停时刻释放污染物后的运动情况。
18.本发明具有的优点和技术效果：本发明通过建立滨海电厂海域的污染物运动数值模型，通过统计污染物进入取水口的概率，提出污染物重点监测范围及方位，为现场监控提供指导。通过结合现场监测及数值模拟的方法，可快速对滨海电厂周围海域已发现的污染物的运动轨迹进行计算，并可预测污染物运动至取水口门的时间，解决了目前被动清理垃圾、无法有效应对污染物突然涌入污染物的难题，为电厂的运营部门提供预警信息。
19.由于本发明采用上述技术方案，电厂运营部门可加强对工程海域重点区域的监测，在发现大量污染物进入重点区域后，可提前迅速启动应急预案，调配应急资源，有效解决在应对污染物突然涌入取水口、无法及时清理垃圾导致机组停堆的难题，提高电厂的运营效率。
附图说明
20.图1是本发明海洋污染物运动轨迹预测及滨海电厂取水安全预警流程图。
21.图2本发明应用实施例工程范围地形及网格图；
22.图3是潮位验证曲线图；
23.图4是流速流向验证曲线图；
24.图5是污染物释放位置示意图；
25.图6是wsw方向距堤头1km位置、涨急时刻释放污染物运动轨迹图；
26.图7是涨急时刻不同位置污染物进入明渠通量概率分布图。
具体实施方式
27.为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明
进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
28.请参阅图1，一种海洋污染物运动轨迹预测及滨海电厂取水安全预警方法，其特征在于：包括如下步骤：
29.(1)建立潮流数学模型：首先根据工程区域的特征，建立大范围潮流数学模型，外海边界距离工程位置50～100km；数学模型基于流体不可压缩性及navier～stokes方程，并服从boussinesq假定和静水压力假定；利用潮汐预报软件chinatide提供潮位边界资料，通过对水文实测资料进行验证，保证该工程海域下数学模型的各项参数有效性，确保潮流数学模型的准确性；
30.(2)建立粒子追踪模型：将呈现离散状、无自游泳能力或弱自游泳能力进入明渠的漂浮类污染物在数学模型中将其概化为粒子；粒子追踪模型采用langevin方程来描述粒子的迁移运动；
31.(3)确定粒子源释放点：以工程取水口为中心，将其周围海域根据距离及方位分别设置若干个粒子源释放点；
32.(4)统计粒子的通量及分布情况：根据潮流数学模型和粒子追踪模型计算不同点在不同时刻、不同潮型释放污染物后，粒子随潮流及明渠取水作用下的运动轨迹，并统计粒子通过取水口拦污网位置的通量及分布情况；此时完成污染物运动轨迹预测；
33.(5)制作概率云图：根据统计结果，制作污染物进入明渠的概率云图，并评估污染物对取水安全的影响，通过概率云图可判定污染物来向的危险程度，重点加强对该方位的监测，通过提供观测设备的性能，监测范围及距离越远，预留应对措施的时间；
34.(6)实时监测及标记：根据现场对污染物开展实时监测，对已进入工程区域观测范围的污染物进行及时标记，掌握其距取水口位置、方位及时间；
35.(7)预测已标记粒子的后续运动轨迹：将已标记粒子信息输入到粒子数学模型中，通过粒子数学模型计算，预测已标记粒子的后续运动轨迹并计算至取水口的时间，为工程运营单位提出预警，以及时采取应急措施，包括提前调配人员及设备以提高清理能力，或根据污染物的类型选择不同的打捞方案。
36.本发明进一步采取的技术方案如下：
37.优选的，步骤1中的实测水文资料主要包括潮位、流速、流向资料。之所以选用上述潮位、流速、流向水文资料，其目的在于通过以上参数的验证，可确保数学模型能真实模拟工程区域的潮流特征，可进一步为粒子的运动模拟提供可靠的数学模型基础。
38.优选的，步骤2中的污染物指无自游泳能力或弱自游泳能力的物体，包括泡沫、秸秆、树枝、塑料及水草中的一种或多种。根据已发生的取水安全事故，该类型污染物种类较多，且相较于具有自主游泳能力的生物类型污染物，其水流跟随性良好，因此可通过数学模型对此类型污染物进行模拟。
39.优选的，步骤3中依滨海电厂取水口为中心，以距离10km范围内、每间隔1km设置一个粒子释放点，且以22.5
°
为一个方向。采用该技术方案具有的优点在于将工程海域空间位置离散化、精确化，为制作污染物进入明渠的概率云图提供基础；另外根据不同海域的潮流特征，10km范围外的污染物在一个潮流周期内进入明渠的概率非常小，因此对10km范围内的粒子运动进行模拟能够满足预警的需求。
40.优选的，步骤4中计算不同点位在涨急时刻、涨停时刻、落急时刻、落停时刻释放污染物后的运动情况。由于潮流流速存在周期变化特征，采用四个时刻目的在于将时间离散化，对不同时刻下粒子的运动进行模拟。
41.具体工程实施例：
42.(1)、首先建立工程区域的潮流数学模型，并根据实测水文资料验证模型的准确性；
43.(1.1)针对某电厂工程首先建立潮流数学模型，计算网格见图1，计算区域南北向长约54km，东西宽约76km，外边界至-50米等深线；
44.(1.2)对潮流数学模型进行验证；
45.根据实测水文资料主要包括潮位、流速、流向资料，验证结果见图2和图3，可以看出潮位、流速、流向的计算结果与实测资料符合良好，满足规范(《水运工程模拟试验技术规范》jts/t 231-2021)要求，因此模型可以用于工程区域潮流计算。
46.(2)建立粒子追踪模型：进入明渠的漂浮类污染物主要是塑料、泡沫以及木头等，其总体呈现离散状，在数学模型中将其概化为粒子。粒子追踪模型采用langevin方程来描述粒子的迁移运动。
47.(3)
48.3.1)首先以工程取水口为中心，将周围海域分为16个方向，即e、ene、
…i…
ene(每22.5
°
为一个方向，其中以e为起点，逆时针旋转，ene为终点，对应i从1至16进行编号)，作为污染物出现的方向；
49.3.2)在每个方位距离口门不同位置设置若干污染物来源点，如1km、2km、
…j…
nkm；
50.由于污染物的来源只能为水体中，假设工程周围nkm范围内均为水体，则以上污染物的来源位置共有16n个。
51.3.3)以上各位置可表示为；
52.假定取水口门位置为(x0,y0)，各点位置为:
[0053][0054][0055]
分别计算每个位置处的污染物在潮流及取水作用下的运动轨迹。并根据以下方法确定该处污染物进入明渠的概率：
[0056]
(4)根据潮流数学模型和粒子追踪模型计算不同点在不同时刻、不同潮型释放污染物后，粒子随潮流及明渠取水作用下的运动轨迹，并统计粒子通过取水口拦污网位置的通量及分布情况；此时完成污染物运动轨迹预测；
[0057]
监测污染物通过取水口拦污网的污染物通量；其中污染物释放持续1个小时，以污染物通量百分比pp评估各工况下污染物的堆积，其中：
[0058]
假定第k个工况下污染物位于x(i,j)，y(i,j)位置，该污染点持续释放1个小时的污染量，每秒释放污染物量为pt，共释放总量为tp，计算该工况下，污染物的运动轨迹，并进入明渠口的污染物量ip，则该工况下，污染物进入明渠的比例(概率)为
[0059]
pp＝ip/tp
[0060]
pp——通过某断面污染物通量百分比
[0061]
ip——通过某断面的污染物通量
[0062]
tp：污染物释放总量
[0063]
其中污染物通量统计时间为从污染物释放时刻后24h以内。污染物释放位置示意见图5，本实施例中通过模型计算可以直观看到污染物在出现后，随潮流及明渠取水作用下的运动轨迹。以wsw方向、距离堤头1km处在涨急时刻释放粒子后的运动轨迹，进行说明；具体参见图5，图中虚线表示轨迹，空心圆点表示污染物初始释放位置，实心圆点表示24小时后的最后位置。
[0064]
由图可以观测到，污染物出现后，主要随涨潮潮流向西运动，至落潮期间，随落潮流向东运动，至口门附近，由于取水的作用，污染物向明渠内运动，并在靠西堤侧向明渠内运动，因此可以判断西堤侧为污染物堆积区域，因此应重点在西堤侧部署垃圾清运装置，并在设计时应予考虑。
[0065]
(5)根据各工况下统计结果，制作污染物进入明渠的概率云图，并评估污染物对取水安全的影响；本实施例中根据以上计算的方法，重复计算不同工况下的污染物运动情况，统计得到如图7的概率云图。
[0066]
由概率云图可知，进入明渠的污染物主要以wsw～w向为主，不同时刻进入明渠的通量大小有较大差别，总体来看，当污染物距离堤头超过6km，在一个潮周期内，几乎没有污染物进入明渠。因此在现场观测时，应重点加强对wsw～w向内6km范围内的污染物进行监测，根据现场监控设备的性能，可适当扩大监测范围。
[0067]
(6)实时监测及标记：根据现场对污染物开展实时监测，对已进入工程区域观测范围的污染物进行及时标记，掌握其距取水口位置、方位及时间；
[0068]
(7)根据现场对污染物的实时监测，对已进入工程区域观测范围的污染物进行及时记录；将污染物信息输入到粒子数学模型中，预测该污染物的后续运动轨迹并计算至取水口的时间，为滨海电厂运营单位提出预警，以及时采取应急措施。本实施例中具体操作为；通过观测现场污染的情况，将污染物位置、时刻输入数学模型，就可以得到其运动轨迹及到达明渠的时间，为滨海电厂运营单位提出预警，可以提前为应急管理预留一定的时间，以及时采取应急措施，包括提前调配人员及设备以提高清理能力，或根据污染物的类型选择不同的打捞方案。
[0069]
通过上述技术方案通过现场观测的方式对滨海电厂取水口附近海域的漂浮类污染物进行标记，并对其后续运动轨迹作出预测，计算其进入明渠的概率及到达时间，直接为滨海电厂冷源安全部门提出预警，降低滨海电厂因污染物堆积导致机组停堆的风险，提高滨海电厂运营效率。
[0070]
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。