CN116468275B - 考虑致灾生物运动特征及节律行为的取水风险评估方法 - Google Patents
考虑致灾生物运动特征及节律行为的取水风险评估方法 Download PDFInfo
- Publication number
- CN116468275B CN116468275B CN202310430759.2A CN202310430759A CN116468275B CN 116468275 B CN116468275 B CN 116468275B CN 202310430759 A CN202310430759 A CN 202310430759A CN 116468275 B CN116468275 B CN 116468275B
- Authority
- CN
- China
- Prior art keywords
- water
- disaster
- water intake
- causing
- flow field
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 280
- 230000033001 locomotion Effects 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000012502 risk assessment Methods 0.000 title claims abstract description 21
- 230000006399 behavior Effects 0.000 title claims abstract description 15
- 230000033764 rhythmic process Effects 0.000 title claims abstract description 14
- 238000013508 migration Methods 0.000 claims abstract description 56
- 230000005012 migration Effects 0.000 claims abstract description 56
- 238000013461 design Methods 0.000 claims abstract description 11
- 238000010168 coupling process Methods 0.000 claims abstract description 7
- 230000027387 rhythmic behavior Effects 0.000 claims abstract description 7
- 230000008878 coupling Effects 0.000 claims abstract description 6
- 238000005859 coupling reaction Methods 0.000 claims abstract description 6
- 230000008859 change Effects 0.000 claims description 27
- 230000008569 process Effects 0.000 claims description 24
- 238000004088 simulation Methods 0.000 claims description 24
- 230000009182 swimming Effects 0.000 claims description 23
- 238000004364 calculation method Methods 0.000 claims description 20
- 238000011156 evaluation Methods 0.000 claims description 12
- 238000011160 research Methods 0.000 claims description 9
- 238000002474 experimental method Methods 0.000 claims description 7
- 238000013178 mathematical model Methods 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 230000002706 hydrostatic effect Effects 0.000 claims description 6
- 238000012876 topography Methods 0.000 claims description 6
- 238000012795 verification Methods 0.000 claims description 6
- 230000009471 action Effects 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 241000238557 Decapoda Species 0.000 description 38
- 241000238552 Penaeus monodon Species 0.000 description 13
- 241000238553 Litopenaeus vannamei Species 0.000 description 12
- 238000005259 measurement Methods 0.000 description 10
- 230000027288 circadian rhythm Effects 0.000 description 4
- 239000000498 cooling water Substances 0.000 description 4
- 230000000630 rising effect Effects 0.000 description 3
- 241000995704 Fenneropenaeus chinensis Species 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000009545 invasion Effects 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 239000013535 sea water Substances 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 235000003826 Artemisia Nutrition 0.000 description 1
- 235000003261 Artemisia vulgaris Nutrition 0.000 description 1
- 241000238421 Arthropoda Species 0.000 description 1
- 239000002028 Biomass Substances 0.000 description 1
- 241000195493 Cryptophyta Species 0.000 description 1
- 206010013554 Diverticulum Diseases 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 241000927735 Penaeus Species 0.000 description 1
- 235000014001 Prunus serrulata Nutrition 0.000 description 1
- 241000392970 Prunus serrulata Species 0.000 description 1
- 241000242583 Scyphozoa Species 0.000 description 1
- 241001125843 Trichiuridae Species 0.000 description 1
- 230000003872 anastomosis Effects 0.000 description 1
- 244000030166 artemisia Species 0.000 description 1
- 235000009052 artemisia Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000002060 circadian Effects 0.000 description 1
- 230000005895 circadian behavior Effects 0.000 description 1
- 239000003653 coastal water Substances 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000026058 directional locomotion Effects 0.000 description 1
- 238000004836 empirical method Methods 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000011010 flushing procedure Methods 0.000 description 1
- 238000010030 laminating Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000010865 sewage Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Water Supply & Treatment (AREA)
- Game Theory and Decision Science (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Public Health (AREA)
- Quality & Reliability (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Primary Health Care (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种考虑致灾生物运动特征及节律行为的取水风险评估方法，包括：构建包括电厂取水工程在内的取水海域流场模型；流场模型结果验证及流场分析；构建典型取水致灾生物迁移模型；典型取水致灾生物迁移模型与取水海域流场模型耦合；取水风险评估。本发明针对核电取水风险评估方法无法反应各电厂取水设计、海域特点以及各致灾生物差异性的局限性，提出一种利用流场模型和生物迁移模型耦合的方式，充分考虑致灾生物运动特征、节律行为、取水工程设计特征及取水流场变化的取水风险评估方法，使得取水风险评估的结果能够更客观、准确的反应实际的各电厂取水差异性情况，为核电厂取水风险提供科学有效的评估结果，为保障核电安全提供技术支持。
Description
技术领域
本发明涉及一种考虑致灾生物运动特征及节律行为的取水风险评估方法，是一种水工工程运行安全性的评估方法，是一种电厂冷却取水工程运行安全风险的预测或反演的方法。
背景技术
为了减少近岸生物灾害对核电取水安全的影响，需充分明晰和评估核电取水工程运行的安全性。然而，目前采用的评估方法主要基于专家经验法，提取核电取水系统中的关键指标，如拦污网设置、监测系统配备及所处海域致灾生物密度等进行评估，存在一定的局限性，没有充分考虑由于致灾生物运动特点、生物节律特征以及核电取水汇流引起的局部海域流场变化等综合造成的取水风险。例如，同处于一海湾中的两个毗邻电厂(电厂A、电厂B)取水，均采用明渠取水。由于两个电厂毗邻，因此，海水过滤系统设计及布置基本一致，海流、海洋水文气象条件一致。以传统的基于经验法的取水评估方法，两电厂面临致灾生物应当一致，其取水风险也应当一致。然而，历次该海湾中致灾生物毛虾暴发时，电厂B冷源安全问题频发，而电厂A受到影响较小，以此证明现有的评估方法存在很大问题。如何综合考虑多项因素，提高取水风险评估的准确性，是一种亟需解决的问题。
发明内容
为了克服现有技术的问题，本发明提出了一种考虑致灾生物运动特征及节律行为的取水风险评估方法。所述的方法综合考虑致灾生物随水流运动的特点及其在关键生活史阶段体现出的生物节律行为、取水流场分布规律等多重影响因素的核电厂取水风险，为避免海洋生态灾害所造成的核电取水冷源丧失事件提出了科学的依据。
本发明的目的是这样实现的：一种考虑致灾生物运动特征及节律行为的取水风险评估方法，所述方法的步骤如下：
步骤1，构建包括电厂取水工程在内的取水海域流场模型：获取取水海域地形、潮流条件、岸线特征，以及电厂取水工程的设计尺寸、取水流速、取水水深等，建立取水海域流场数学模型，确定模拟范围，划分计算网格，确定流场模拟的边界条件、初始条件并选择合适的计算参数，选择计算控制潮型，模拟包含核电取水工程，以及取水影响海域在内的海域流场模拟；流场模型以水深积分形式的浅水方程为基础，模型采用基于Boussinesq和垂向静水压力假定并采用曲线坐标，垂向上使用σ垂向坐标系统：
其中：d为参考平面以下的水深值，H为全水深，z为物理空间垂向坐标，ζ为自由水面高程；
模型控制方程包括连续方程、动量方程及密度状态方程，并通过湍流理论使方程闭合，具体控制方程如下：
水深平均的连续方程：
其中，ζ为自由水面高程；t为时间；ξ和η为曲线坐标系下水平方向的坐标；Gξξ、Gηη为曲线坐标系转换为直角坐标系的转换系数；U、V分别为ξ和η方向的水深平均流速；Q表示源汇项，如取排水等；
模型在水平方向ξ和η的动量方程：
其中：u、v为曲线正交坐标系下ξ和η方向的流速分量，f为Coriolis参数，f＝2Ωsinφ，Ω为地球自转角速度，φ为纬度；Fξ和Fη分别为ξ和η方向的紊动动量通量，反应雷诺应力；Pξ和Pη为ξ和η方向的水压力梯度；Mξ和Mη为ξ和η方向上动量的源汇项；
由于基于浅水方程，模型基于静水压力：二维模型不考虑水体密度ρ0变化。
床底阻力采用二次型形式：
其中：谢才系数C2D根据谢才公式计算，n为糙率系数，可依据实测资料率定，或根据经验取值，范围可在0.01～0.03之间；
步骤2，流场模型结果验证及流场分析：采用实测数据对流场模拟结果进行潮位及流场的验证，以确保计算的潮位过程与实测过程基本一致或整体上计算潮位与实测值符合；
步骤3，构建典型取水致灾生物迁移模型：根据致灾生物的运动方向、运动速度，随昼夜时间在水体中的位置变化，及其运动特征和节律行为提出致灾生物迁移模型；
模拟方程如下：
(1)致灾生物运动方程：
或
ΔS＝(ulive+uflow)×Δt，
其中：ΔS为致灾生物运动位移；Δt为致灾生物运动时间步长；ulive为致灾生物游动速度矢量；uflow为该时间步长下致灾生物所在位置水流速度矢量；
(2)致灾生物节律行为运动方程：
Height＝f(Time)
其中：Height为致灾生物所在位置水深，f(Time)为致灾生物在水体中所处的垂向位置随时刻变化的相关关系，该函数关系通过对致灾生物昼夜垂直分布的相关研究或实验获得；
步骤4，典型取水致灾生物迁移模型与取水海域流场模型耦合：取水海域流场模型的输出结果将作为致灾生物迁移模型的输入条件，驱动致灾生物在流场变化下进行迁移运动，水力信息单向传递给致灾生物迁移模型，根据生物的迁移运动规律获取生物运动状态、生物所处位置随时间的变化，进而获得致灾生物在取水环境流场作用下的迁移过程；
步骤5，取水风险评估：以一个完整的典型潮作为计算时间总长度，以进入取水明渠的致灾生物占初始投放数量的百分比作为评价取水风险的指标，基于取水海域流场变化规律、致灾生物迁移模型计算结果，分析得到对应环境水流和取水工程结构下，取水运行时受到取水汇流作用影响进入明渠的致灾生物占初始投放致灾生物数量的百分比评价取水风险，100％为取水风险评估最大值，即初始投放致灾生物全部进入取水明渠，0％为取水风险评估最小值，即初始投放致灾生物均没有进入取水明渠。
本发明的优点和有益效果是：本发明针对目前核电取水风险评估方法无法反应各电厂取水设计、取水海域特点及各致灾生物差异性的局限性，提出了一种利用流场模型和生物迁移模型耦合的方式，充分考虑致灾生物运动特征、节律行为及取水工程设计特点、取水流场变化特征的取水风险评估方法，使得取水风险评估的结果能够更客观、准确的反应实际的各电厂取水风险差异性情况，为核电厂取水风险提供科学有效的评估结果，为保障核电安全提供技术支持。
附图说明
下面结合附图和实施例对本发明作进一步说明。
图1是本发明实施例所述方法的流程图；
图2是本发明实施例所述σ垂向坐标示意图。
具体实施方式
实施例：
本实施例是一种考虑致灾生物运动特征及节律行为的取水风险评估方法。所述的方法充分考虑各电厂海域水文、地形、潮流条件的不同，致灾生物运动特征、节律行为的特点及取水工程设计的特征，将这些重要的因素均纳入取水风险评估中，体现各电厂面临取水环境、致灾生物及自身设计特点的差异性所导致的取水风险不同，为电厂取水风险防控，取水安全运行提供依据。为便于说明，本实施例以上述两个毗邻的核电厂(电厂A、电厂B)取水为应用实例，对取水风险评估的过程进行具体描述。
电厂A和电厂B两电厂冷却水均采用直流冷却系统，以海水作为循环冷却水源，温排水排入厂址附近海域。
电厂A和电厂B所处海湾的东、北、西三面为低山丘陵环抱。湾内南北长约35km，东西向最宽处约25km，最窄处约14km。湾内海区等深线分布大致与岸线平行，湾口处的海床最深，海湾面积约600km2。
两电厂附近潮汐观测站的统计结果表明，所处海湾的潮汐性质属于不正规半日混合潮，涨落潮流以南北向往复流为主，涨潮流流向以北向为主，落潮流以南向为主。根据实测资料显示，夏季最大涨潮流速为105cm/s，流向为322°；最大落潮流速为66cm/s，流向为178°。冬季最大涨潮流速为57cm/s，流向为334°；最大落潮流速为56cm/s，流向为142°。
威胁电厂取水冷源安全的主要生物为毛虾(中国毛虾，Acetes chinesis，隶属节肢动物门、甲壳纲、十足目、樱虾科、毛虾属)。如2015年1月10日，电厂B因毛虾涌入取水口导致反应堆降功率运行，期间打捞毛虾总量超20吨；同年11月-12月，电厂B又遭遇毛虾数次入侵，取水系统反冲洗排水渠收集毛虾量超1.2吨/小时；2016年1月8日-9日，电厂B再次因毛虾入侵导致机组停堆。电厂B冷源安全问题频发，但距其取水位置仅800米的电厂A未受到显著威胁。本实施例所述考虑致灾生物运动特征及节律行为的取水风险评估方法的步骤如下，流程如图1所示：
步骤1，构建包括电厂取水工程在内的取水海域流场模型：获取取水海域地形、潮流条件、岸线特征，以及电厂取水工程的设计尺寸、取水流速、取水水深等，建立取水海域流场数学模型，确定模拟范围，划分计算网格，确定流场模拟的边界条件、初始条件并选择合适的计算参数，如海床糙率系数、涡黏系数等，选择计算控制潮型，模拟包含核电取水工程，以及取水可能影响海域在内的海域流场模拟。流场模型以水深积分形式的浅水方程为基础，模型采用基于Boussinesq和垂向静水压力假定并采用曲线坐标，垂向上可以使用σ垂向坐标系统(见图2)，即将物理空间垂向坐标z转换为σ，采用σ垂向坐标，可将自由水面和不规则地形海床变成坐标系中的σ＝0和σ＝-1的坐标平面，使得计算水域垂向上具有相同的网格数，提高了浅水部分垂向分辨率，地形贴合性较好：
其中：d为参考平面以下的水深值，H是全水深，z为物理空间垂向坐标，ζ为自由水面高程(参考平面z＝0以上的数值)。在海床床底σ＝-1，自由水面σ＝0。
图2中网格中垂向分层线的水深是不同的，但是距离自由水面(倒三角表示的是自由水面)的相对水深是一致的，阴影部分表达的海底。
模型控制方程包括连续方程、动量方程及密度状态方程，并通过湍流理论使方程闭合，具体控制方程如下：
水深平均的连续方程：
其中，d为参考平面以下的水深值，ζ为自由水面高程(参考平面z＝0以上的数值)；t为时间；ξ和η为曲线坐标系下水平方向的坐标；Gξξ、Gηη为曲线坐标系转换为直角坐标系的转换系数；U、V分别为ξ和η方向的水深平均流速；Q表示源汇项，如取排水、降雨和蒸发等。
模型在水平方向ξ和η的动量方程：
其中：u、v为曲线正交坐标系下ξ和η方向的流速分量，f为Coriolis参数，f＝2Ωsinφ，Ω为地球自转角速度，φ为纬度；Fξ和Fη分别为ξ和η方向的紊动动量通量，反应雷诺应力；Pξ和Pη为ξ和η方向的水压力梯度；Mξ和Mη为ξ和η方向上动量的源汇项。
由于基于浅水方程，模型基于静水压力：二维模型不考虑水体密度ρ0变化。
床底阻力采用二次型形式：
其中：谢才系数C2D可根据谢才公式计算，n为糙率系数，可依据实测资料率定，或根据经验取值，范围可在0.01～0.03之间。
应用实例中电厂A和电厂B毗邻，距离很近(两个电厂的取水明渠直线距离不超过约800m)，考虑到电厂所处海湾的湾潮流涨落潮特征、取水可能影响的海域范围以及致灾生物来源情况，取水海域流场模型构建范围以A电厂和B电厂为重点关注区域，包括整个海湾，及其外围部分区域，大约1000km2的海域。重点模拟区域为核电取水所在的B电厂海域。模型采用非结构化网格，网格数量为57851，节点数量为30046个，对靠近核电厂取水重点关注区域的网格进行局部加密，网格尺寸约30-80m，其余模拟区域网格尺寸约为80-1000m。
初始条件：采用静流条件起算，即从流速为零的静止流场开始计算。
边界条件：岸边界为固定壁面，以不透水设计，即水体与岸边无水流交换；开边界采用潮位边界，潮位由TPXO全球海潮模型给定。
取排水条件设计：电厂A和电厂B取水流量分别为：90m3/s、225m3/s。
海床糙率系数确定：依据实测资料率定，选取糙率系数n取值0.018-0.025之间。
涡黏系数确定：天然海域的水流均为湍流。数学模型需要特定的湍流模型闭合数值方程。湍流求解中一般通过引入广义涡黏系数ν求解湍流，单位m2/s。水平黏性系数取值多数跟一般跟局部流速、水深有关，且应与网格尺度的大小相协调。二维模拟中广义水平黏度νH包含二部分：第一部分为分子黏度νm，在海域模拟中相对而言为极小的值；第二部分为水平湍动黏度，可称其为数值黏度；本实例中水平湍动黏度采用Smagorinsky公式(包含局部水深、流速和网格尺度信息)计算，并增加率定的背景值根据多次全潮水文测验结果率定后，水平湍动系数背景值/>取为1.0m2/s。
计算控制潮型确定：根据电厂A、电厂B的海洋生物进入明渠的观测结果，夏季小潮时核电厂明渠内生物量较高，即表明此潮型下取水对海洋生物的影响相对较强，故此，本应用实例中选取夏季小潮为计算控制潮型。
步骤2，流场模型结果验证及流场分析：采用实测数据对流场模拟结果进行潮位及流场的验证，以确保计算的潮位过程与实测过程基本一致或整体上计算潮位与实测值符合。
流场模型计算的整体流场特性与观测结果应尽量吻合，需进行定点潮位、流速与流向过程验证。通过调整模型糙率、开边界水力参数等方式，达到模拟准确度的要求，以保证流场的相似性。具体包括：潮位，高低潮时间的相位允许偏差为±0.5h，最高最低潮位值允许偏差为±0.1m；流速，憩流时间和最大流速出现的时间允许偏差为±0.5h，流速过程线的形态基本一致，涨、落潮段平均流速允许偏差为±10％；流向，往复流时测点主流流向允许偏差为±10°，平均流向允许偏差为±10°；旋转流时测点流向允许偏差为±15°；流场中各固定点模拟计算的潮位、流速、流向与实测过程中固定点的实测潮位、流速、流向应基本吻合，即构建的取水海域流场模型计算的流场能够反映取水海域本身海流变化的基本特征以及取水工程运行下取水汇流流场分布变化特征，可以此为基础，开展致灾生物对取水安全风险影响的评估。
接上述A、B电厂应用实例：
应用实例中采用实测资料对模型流场结果进行了验证，潮位验证结果表明，计算的潮位过程与同一时间实测得到的潮位过程进行对比，个别高低潮时刻潮位略有偏差，偏差小于0.1m，总体吻合较好，潮位偏差符合相关规程和规范的要求。夏季小潮定点测量的流速、流向计算值与同一时间实测得到的流速、流向值进行对比，总体基本吻合。验证总体结果表明：夏季小潮期间，模型计算整体流场分布特性与实测结果符合良好，各固定点模拟计算的流速、流向与实测过程也基本吻合，即本应用实例的研究模型计算流场能够反映工程海区海流的基本特征，可据此开展致灾生物对取水安全风险影响的评估。
步骤3，构建典型取水致灾生物迁移模型：根据致灾生物的运动方向、运动速度，随昼夜时间在水体中的位置变化，及其运动特征和节律行为提出水致灾生物迁移模型。
取水海域中的生物和杂物容易受到取水汇流作用的影响进入取水明渠，进而影响取水安全，从已经出现的取水冷源丧失事件中可以观察到，中国沿海影响取水的主要威胁生物多为毛虾、水母、藻类，以及没有主动游动能力或具有较弱游动能力的生物等。典型取水致灾生物迁移模型的构建，充分考虑致灾生物的运动特征及其在运动迁移过程中体现出的节律行为特征等，包括：致灾生物的运动方向、运动速度，随昼夜时间在水体中的位置变化节律等，根据其运动特征及节律行为提出致灾生物迁移模型，模拟方程如下：
(1)致灾生物运动方程：
或
ΔS＝(ulive+uflow)×Δt，
其中：ΔS为致灾生物运动位移；Δt为致灾生物运动时间步长；ulive为致灾生物运动速度矢量，与水流流向相同，通过对致灾生物运动的相关研究或实验可以获得；uflow为该时间步长下致灾生物所在位置水流速度矢量；
(2)致灾生物节律行为运动方程：
生物在水体中的运动，除了表现出随流运动的趋向性特征外，还随着昼夜时间的变化在水深方向上呈现出一定的垂向迁移规律，即生物昼夜节律行为，即致灾生物在水体中的垂向位置是与时间相关的函数，控制方程如下：
Height＝f(Time)
其中：Height为致灾生物所在位置水深，f(Time)为致灾生物在水体中所处的垂向位置随时刻变化的相关关系，该函数关系通过对致灾生物昼夜垂直分布的相关研究或实验获得。
应用实例中电厂A和电厂B的取水主要致灾生物为毛虾，根据相关研究及实验结果显示，毛虾体型小，体长约为1～2cm，因游动能力较弱不作长距离洄游，呈现明显的季节性定向的浅、深水之间的移动，随潮流推移游动于近岸区域，游动方向与水流方向一致。电厂所处海湾的海域毛虾随昼夜时间在水体中的位置变化表现为白天主要栖息于中下层或近底层，夜晚向底层迁移。毛虾迁移过程由水流流速、毛虾自身游动速度及其随昼夜在水体中垂向迁移能够到达的位置决定，其中水流流速和流向通过大范围取水流场数学模型获得，毛虾游动速度和垂向迁移可以通过实验或国内外相关研究获得。
虾类的游动速度可通过实验测定，相关研究表明，虾的游动速度与其体长密切相关。实验和观测结果表明：虾类游动速度约为其体长的2-3倍，即毛虾相对体长游动速度为2.5±0.5BL/s，其中BL(body length)即为虾体长。因此，在案例中，毛虾的游动能力参考游动速度，即ulive，设定范围为：0-6cm/s。在迁移模拟过程中，每一个时间步长开始时，毛虾个体都将从速度范围中随机获得一个游动速度，进行该时间步长内的迁移运动。毛虾自身游动体现为顺流的行为习性，其迁移方向始终跟随水流方向的变化而改变。
毛虾顺水流方向运动，运动位移ΔS为：
ΔS＝(ulive+uflow)×Δt，
其中：ΔS为顺流方向上毛虾的位移；ulive为毛虾游动速度矢量，与水流流向相同；uflow为该时间步长下毛虾所在位置水流速度矢量；Δt为时间步长。
毛虾垂向迁移主要由其自身的昼夜迁移节律决定，根据电厂所处海湾的海域毛虾调查研究结果，毛虾大部分时间主要分布在5.2～7.2m的水深位置，夜间(22：50-00：10)出现明显的垂向迁移，主要分布在7.2～8.2m的水深位置。因此，在本应用实例中，毛虾的垂向所处位置时随时间变化的分段阈值，即在迁移模拟过程中，每一个时间步长开始时，判断时间若处于00：10-22：50之间，则毛虾个体从水深5.2～7.2m范围中随机获得一个水深，进行后续的迁移运动；判断时间若处于22：50-00：10之间，则毛虾个体从水深7.2～8.2m范围中随机获得一个水深，进行该时间步长内的迁移运动。
毛虾迁移体现在其自身游动和受昼夜节律影响的垂向迁移，其中：
受昼夜节律影响的垂向迁移可以表达为：
其中：Height为毛虾所处水深，Time为模拟对应的时刻，f(Time)为致灾生物在水体中所处的垂向位置随时刻变化的相关关系。
步骤4，典型取水致灾生物迁移模型与取水海域流场模型耦合：取水海域流场模型的输出结果将作为致灾生物迁移模型的输入条件，驱动致灾生物在流场变化下进行迁移运动，水力信息(各点流速、流向、水深)单向传递给致灾生物迁移模型，根据生物的迁移运动规律获取生物运动状态、生物所位置随时间的变化，进而获得致灾生物在取水环境流场作用下的迁移过程。
本应用实例中，以毛虾群体中的单个个体为例，说明毛虾个体运动迁移模型与流场模拟结果耦合过程如下，毛虾群体的运动即是单个个体运动的总和：
1)本时刻t0，根据毛虾个体所处位置(x0，y0，z0)和取水海域流场模型获得初始时刻毛虾所在位置水流的流速uflow＝(ux，uy，uz)，其中z0为根据t0时刻在H＝f(t0)的取值，获得毛虾在垂向所处水深z0，此时刻，毛虾个体的游动速度为ulive为0-6cm/s速度范围中随机产生的一个游动速度，速度的方向与uflow方向相同。
2)计算时间步长为Δt，毛虾顺水流方向运动，每个时间步长下，毛虾的位移变动ΔS为：
将毛虾游动运动的速度矢量分解为沿x轴的速度分量ulive-x，沿y轴的速度分量ulive-y，沿z轴的速度分量ulive-z，同样的，水流速度矢量分解为沿x轴的速度分量uflow-x，沿y轴的速度分量uflow-y，沿z轴的速度分量uflow-z，并以水流速度矢量与y轴正方向的夹角α，与x轴正方向的夹角β，与z轴正方向的夹角γ为该时间步长下毛虾游动运动方向的标记。则Δt时间步长内，根据ΔS＝(ulive+uflow)×Δt，毛虾个体位置(x0，y0，z0)更新为(xt，yt，zt)，其中：
xt＝x0+ΔSx＝x0+(ulive-x+uflow-x)×Δt＝x0+(ulive+uflow)×sinα×Δt+(ulive+uflow)×cosβ×Δt
zt＝z0+ΔSz＝z0+(ulive-z+uflow-z)×Δt＝z0+(ulive+uflow)×sinβ×Δt+(ulive+uflow)×cosγ×Δt
并根据t+Δt所对应的时刻Time，带入Height＝f(Time)公式，判断zt是否满足毛虾昼夜行为节律，即zt是否Height＝f(Time)的合理取值范围，若满足重复过程1)获得该位置处水流的流速uflow，进行下一个时间步长该毛虾个体的迁移运动；
3)若不满足，则毛虾个体从游动速度范围0-6cm/s中再次随机获得一个游动速度，重复过程2)进行迁移运动；
4)若重复过程2)n次，仍无法满足毛虾垂向位置在合理范围内，ulive＝0，进入下一时间步长。
5)重复过程1)～4)，进行下一个毛虾个体的运动迁移，直至毛虾群体迁移结束。
毛虾迁移模型与流场模拟结果耦合遵循毛虾在迁移过程中，其在水深方向的分布始终满足生物的昼夜行为节律，并始终顺水流方向迁移。
步骤5，取水风险评估：以一个完整的典型潮(即计算控制潮型)作为计算时间总长度，以进入取水明渠口门的致灾生物占初始投放数量的百分比作为评价取水风险的指标，基于取水海域流场变化规律、致灾生物迁移数学模型计算结果，分析得到对应环境水流和取水工程结构下，电厂取水运行时受到取水汇流作用影响进入取水系统的致灾生物占初始投放致灾生物数量的百分比作为评价取水风险，100％为取水风险评估最大值，即初始投放致灾生物全部进入取水明渠，0％为取水风险评估最小值，即初始投放致灾生物均没有进入取水明渠。
应用实例中，评估计算结果显示：
(1)致灾生物随潮运移，其中大部分致灾生物表现为随潮向外部海域运移，少部分致灾生物随潮进入取排水口前缘水域以及电厂B海域。
(2)受核电取水口前缘海域流场影响，来自外海的致灾生物呈现自外向内逐渐堆积的分布特征，具体表现为：第1个潮周结束时，致灾生物在电厂B排水明渠南侧附近；第2潮周结束时，致灾生物刚刚开始进入电厂B取水明渠；第4个潮周结束时，致灾生物少量分别进入电厂A和和电厂B取水明渠；第6个潮周结束时，致灾生物大量进入电厂A和电厂B取水明渠。
以进入取水明渠口门的致灾生物占初始投放数量的百分比(即取水卷吸百分数)作为评价取水风险的指标，基于取水海域流场变化规律、致灾生物迁移数学模型计算结果表明：各取水口取水卷吸百分数：第4个潮周结束时刻，电厂A取水卷吸百分数特征值为0.4％，电厂B取水卷吸百分数特征值为2.0％；第6个潮周结束时刻，电厂A取水卷吸百分数特征值为1.9％，电厂B取水卷吸百分数特征值为6.9％；第10个潮周结束时刻，电厂A取水卷吸百分数特征值为4.7％，电厂B取水卷吸百分数特征值为14.0％；第15个潮周结束时刻，即计算控制潮型，夏季典型小潮的最后一个潮周计算结束，电厂A取水卷吸百分数特征值为6.5％，电厂B取水卷吸百分数特征值为19.5％。
即电厂A取水风险评估值为6.5，电厂B取水风险评估值为19.5，两个核电厂取水风险差异显著，电厂A取水风险约为电厂B取水风险的3倍，这一结果与目前掌握的A、B电厂的实际取水风险数据基本吻合。
最后应说明的是，以上仅用以说明本发明的技术方案而非限制，尽管参照较佳布置方案对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案(比如电厂周边地形环境、各种公式的运用、步骤的先后顺序等)进行修改或者等同替换，而不脱离本发明技术方案的精神和范围。
Claims (1)
1.一种考虑致灾生物运动特征及节律行为的取水风险评估方法，其特征在于，所述方法的步骤如下：
步骤1，构建包括电厂取水工程在内的取水海域流场模型：获取取水海域地形、潮流条件、岸线特征，以及电厂取水工程的设计尺寸、取水流速、取水水深，建立取水海域流场数学模型，确定模拟范围，划分计算网格，确定流场模拟的边界条件、初始条件并选择合适的计算参数，计算参数包括：海床糙率系数、涡黏系数；选择计算控制潮型，模拟包含核电取水工程，以及取水可能影响海域在内的海域流场模拟，流场模型以水深积分形式的浅水方程为基础，模型采用基于Boussinesq和垂向静水压力假定并采用曲线坐标，垂向上使用σ垂向坐标系统：
其中：d为参考平面以下的水深值，H为全水深，z为物理空间垂向坐标，ζ为自由水面高程，在海床底部σ＝-1，在自由水面σ＝0；
模型控制方程包括连续方程、动量方程及密度状态方程，并通过湍流理论使方程闭合，具体控制方程如下：
水深平均的连续方程：
其中，ζ为自由水面高程；t为时间；ξ和η为曲线坐标系下水平方向的坐标；Gξξ、Gηη为曲线坐标系转换为直角坐标系的转换系数；U、V分别为ξ和η方向的水深平均流速；Q表示源汇项；
模型在水平方向ξ和η的动量方程：
其中：u、v为曲线正交坐标系下ξ和η方向的流速分量，f为Coriolis参数，f＝2Ωsinφ，Ω为地球自转角速度，φ为纬度；Fξ和Fη分别为ξ和η方向的紊动动量通量，反应雷诺应力；Pξ和Pη为ξ和η方向的水压力梯度；Mξ和Mη为ξ和η方向上动量的源汇项；
由于基于浅水方程，模型基于静水压力：二维模型不考虑水体密度ρ0变化；
床底阻力采用二次型形式：
其中谢才系数C2D，根据谢才公式计算，n为糙率系数；
步骤2，流场模型结果验证及流场分析：采用实测数据对流场模拟结果进行潮位及流场的验证，以确保计算的潮位过程与实测过程基本一致或整体上计算潮位与实测值符合；
步骤3，构建典型取水致灾生物迁移模型：根据致灾生物的运动方向、运动速度，随昼夜时间在水体中的位置变化，及其运动特征和节律行为提出致灾生物迁移模型；
模拟方程如下：
(1)致灾生物运动方程：
或
ΔS＝(ulive+uflow)×Δt，
其中：ΔS为致灾生物运动位移；Δt为致灾生物运动时间步长；ulive为致灾生物游动速度矢量；uflow为该时间步长下致灾生物所在位置水流速度矢量；
(2)致灾生物节律行为运动方程：
Height＝f(Time)其中：Height为致灾生物所在位置水深，f(Time)为致灾生物在水体中所处的垂向位置随时刻变化的相关关系，该函数关系通过对致灾生物昼夜垂直分布的相关研究或实验获得；
步骤4，典型取水致灾生物迁移模型与取水海域流场模型耦合：取水海域流场模型的输出结果将作为致灾生物迁移模型的输入条件，驱动致灾生物在流场变化下进行迁移运动，水力信息单向传递给致灾生物迁移模型，根据生物的迁移运动规律获取生物运动状态、生物所处位置随时间的变化，进而获得致灾生物在取水环境流场作用下的迁移过程；
步骤5，取水风险评估：以一个完整的典型潮作为计算时间总长度，以进入取水明渠口门的致灾生物占初始投放数量的百分比作为评价取水风险的指标，基于取水海域流场变化规律、致灾生物迁移模型计算结果，分析得到对应环境水流和取水工程结构下，取水运行时受到取水汇流作用影响进入取水系统的致灾生物占初始投放致灾生物数量的百分比作为评价取水风险，100％为取水风险评估最大值，0％为取水风险评估最小值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310430759.2A CN116468275B (zh) | 2023-04-21 | 2023-04-21 | 考虑致灾生物运动特征及节律行为的取水风险评估方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310430759.2A CN116468275B (zh) | 2023-04-21 | 2023-04-21 | 考虑致灾生物运动特征及节律行为的取水风险评估方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116468275A CN116468275A (zh) | 2023-07-21 |
| CN116468275B true CN116468275B (zh) | 2023-09-29 |
Family
ID=87176734
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310430759.2A Active CN116468275B (zh) | 2023-04-21 | 2023-04-21 | 考虑致灾生物运动特征及节律行为的取水风险评估方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116468275B (zh) |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20100030351A (ko) * | 2008-09-10 | 2010-03-18 | 한국원자력연구원 | 원자력발전소의 인적유발사건 영향 평가 시스템 및 그 방법 |
| CN109190718A (zh) * | 2018-11-07 | 2019-01-11 | 苏州热工研究院有限公司 | 一种核电厂海生物探测预警多源信息融合算法 |
| CN112613239A (zh) * | 2020-11-24 | 2021-04-06 | 国家海洋局北海环境监测中心(中国海监北海区检验鉴定中心） | 一种水母灾害预警模型的构建方法 |
| CN112801518A (zh) * | 2021-02-02 | 2021-05-14 | 徐粱钰 | 一种水母灾害风险风险评估方法 |
| WO2021253902A1 (zh) * | 2020-06-18 | 2021-12-23 | 中国长江三峡集团有限公司 | 一种计算多点取水口水头损失系数及支管流量分配的方法 |
| CN114113314A (zh) * | 2021-11-25 | 2022-03-01 | 厦门大学 | 一种用于滨海核电站的海洋致灾生物水声监测系统及方法 |
| CN114492973A (zh) * | 2022-01-18 | 2022-05-13 | 苏州热工研究院有限公司 | 一种核电厂周边海域海洋生物预测方法 |
| CN114519495A (zh) * | 2021-12-29 | 2022-05-20 | 中交天津港湾工程研究院有限公司 | 海洋污染物运动轨迹预测及滨海电厂取水安全预警方法 |
| CN115408887A (zh) * | 2022-10-31 | 2022-11-29 | 交通运输部天津水运工程科学研究所 | 一种近海电厂取水卷吸机制数值模拟与评估方法 |
| KR20220166929A (ko) * | 2021-06-11 | 2022-12-20 | 부경대학교 산학협력단 | 취수구 유입 해양이물질 차단 및 제거 시스템 |
| CN115587662A (zh) * | 2022-10-31 | 2023-01-10 | 广东核电合营有限公司 | 一种核电厂周边海域海洋生物分布预测方法 |
| CN115618651A (zh) * | 2022-11-22 | 2023-01-17 | 交通运输部天津水运工程科学研究所 | 一种近海电厂取水卷吸效应研究的物理模拟方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20220324737A1 (en) * | 2021-04-03 | 2022-10-13 | Power and Water Utility Company for Jubail and Yanbu (MARAFIQ) | Methods for monitoring marine waterbody sustainability at discharge sites |
-
2023
- 2023-04-21 CN CN202310430759.2A patent/CN116468275B/zh active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20100030351A (ko) * | 2008-09-10 | 2010-03-18 | 한국원자력연구원 | 원자력발전소의 인적유발사건 영향 평가 시스템 및 그 방법 |
| CN109190718A (zh) * | 2018-11-07 | 2019-01-11 | 苏州热工研究院有限公司 | 一种核电厂海生物探测预警多源信息融合算法 |
| WO2021253902A1 (zh) * | 2020-06-18 | 2021-12-23 | 中国长江三峡集团有限公司 | 一种计算多点取水口水头损失系数及支管流量分配的方法 |
| CN112613239A (zh) * | 2020-11-24 | 2021-04-06 | 国家海洋局北海环境监测中心(中国海监北海区检验鉴定中心） | 一种水母灾害预警模型的构建方法 |
| CN112801518A (zh) * | 2021-02-02 | 2021-05-14 | 徐粱钰 | 一种水母灾害风险风险评估方法 |
| KR20220166929A (ko) * | 2021-06-11 | 2022-12-20 | 부경대학교 산학협력단 | 취수구 유입 해양이물질 차단 및 제거 시스템 |
| CN114113314A (zh) * | 2021-11-25 | 2022-03-01 | 厦门大学 | 一种用于滨海核电站的海洋致灾生物水声监测系统及方法 |
| CN114519495A (zh) * | 2021-12-29 | 2022-05-20 | 中交天津港湾工程研究院有限公司 | 海洋污染物运动轨迹预测及滨海电厂取水安全预警方法 |
| CN114492973A (zh) * | 2022-01-18 | 2022-05-13 | 苏州热工研究院有限公司 | 一种核电厂周边海域海洋生物预测方法 |
| CN115408887A (zh) * | 2022-10-31 | 2022-11-29 | 交通运输部天津水运工程科学研究所 | 一种近海电厂取水卷吸机制数值模拟与评估方法 |
| CN115587662A (zh) * | 2022-10-31 | 2023-01-10 | 广东核电合营有限公司 | 一种核电厂周边海域海洋生物分布预测方法 |
| CN115618651A (zh) * | 2022-11-22 | 2023-01-17 | 交通运输部天津水运工程科学研究所 | 一种近海电厂取水卷吸效应研究的物理模拟方法 |
Non-Patent Citations (2)
| Title |
|---|
| 核电厂冷源安全的挑战及对策;吴侨军;;电力安全技术(11);全文 * |
| 潜在影响防城港核电冷源系统的藻类暴发特点及其监测防控技术;贺立燕;宋秀贤;於凡;王凯;宋书群;俞志明;;海洋与湖沼(03);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116468275A (zh) | 2023-07-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Buttolph et al. | Two-dimensional depth-averaged circulation model CMS-M2D: Version 3.0, Report 2: Sediment transport and morphology change | |
| Frey et al. | GenCade version 1 model theory and user’s guide | |
| Foreman et al. | The ability of hydrodynamic models to inform decisions on the siting and management of aquaculture facilities in British Columbia | |
| Yang et al. | Integrated modeling of flood flows and tidal hydrodynamics over a coastal floodplain | |
| MacWilliams et al. | An overview of multi-dimensional models of the Sacramento–San Joaquin Delta | |
| CN116468275B (zh) | 考虑致灾生物运动特征及节律行为的取水风险评估方法 | |
| Howarth | Hydrography of the Irish Sea | |
| Zhang et al. | Simulation and analysis of Back siltation in a navigation channel using MIKE 21 | |
| Li et al. | Salinity calculations in the coastal modeling system | |
| Fu et al. | Calculation of flow field and analysis of spawning sites for Chinese sturgeon in the downstream of Gezhouba Dam | |
| MacWilliams et al. | Hydrodynamic simulation of circulation and residence time in Clifton Court Forebay | |
| Kamarainen et al. | Estimates of phosphorus entrainment in Lake Mendota: a comparison of one‐dimensional and three‐dimensional approaches | |
| Xiang et al. | Simulation of Remote-Sensed Chlorophyll Concentration with a Coupling Model Based on Numerical Method and CA-SVM in Bohai Bay, China | |
| Schoenbaechler et al. | TxBLEND model calibration and validation for the Nueces estuary | |
| Mo et al. | A long-term numerical model of morphodynamic evolution and its application to the Modaomen Estuary | |
| Ribeiro | Coupled Modelling of the Tagus and Sado Estuaries and Their Associated Mesoscale Patterns | |
| Yang et al. | Numerical study of circulation and temperature-salinity distributions in the Bras d’Or Lakes | |
| Christensen | Hydrodynamic modeling of newly emergent coastal deltaic floodplains | |
| Wood | Dependence of flow and transport through the Williamson River delta, upper Klamath Lake, Oregon, on wind, river inflow, and lake elevation | |
| Kazhyken et al. | Application of a dispersive wave hydro-sediment-morphodynamic model in the Ria Formosa lagoon | |
| Groza | Modeling changes in hydrodynamics and sediment transport due to inlet modification at Little Lagoon Pass, Alabama | |
| Dalledonne | Analysis and assessment of uncertainties in coastal morphodynamic modelling | |
| Yan et al. | Simulation of different sea-crossing traffic route structures’ effects on the marine environment for the Dalian large-scale offshore airport island | |
| Li et al. | Analysis of the Invasion of Acetes into the Water Intake of the Daya Bay Nuclear Power Base | |
| Ahn et al. | Using ECHO as a decision support tool to manage coastal systems subject to physical changes |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |