Analysis of impact of typhoon track and intensity on cage aquaculturearea in sea area of Maniao Bay
-
摘要:
台风过境产生的巨浪和强流对网箱养殖区具有严重影响。为评估台风路径和强度对海南马袅湾内大型网箱养殖区的影响程度，基于Holland风场模型和Delft3D软件构建了波流耦合模型，并利用观测数据对模型进行了验证。依据历史台风路径及登陆位置，将台风路径分为4种类型，分别选取“威马逊”、“圆规”、“黄蜂”和“泰德”作为代表性台风，通过平移和改变台风强度重构了72场台风过程，模拟了台风期间波流场的变化，统计分析了不同路径和强度的台风过境过程中马袅湾养殖区有效波高及流速的最大值，将其作为评估网箱养殖区受影响强弱的因素，探讨了台风路径和强度对网箱养殖区的影响，并提出了不同强度台风过境时的预警路径范围。结果表明，第I类台风路径的影响范围最广、程度最大，其中“超强台风”等级时，R4路径下有效波高最大值超过5.5 m，流速最大值超过1.5 m·s−1，养殖区受严重影响；第II、第III类台风多数仅造成中等以下影响；第IV类台风路径在所有强度下对网箱养殖区基本无影响。
Abstract:Extreme waves and currents generated by the passage of typhoons have an important impact on the cage culture area. In order to assess the degree of hazard of typhoons with different paths and intensities on the large-scale aquaculture areas with in-cell nets in Maniao Bay, Hainan, we constructed a wave-current coupling model based on the Holland wind field model and the Delft3D software, and the model was validated by using observational data. The typhoon paths were categorized into four types based on the paths and landfall locations of the historical typhoons, and "Rammasun", "Kompasu", "Vongfong" and "Ted" were selected as the representative typhoons of the four types of paths. Seventy-two typhoon processes were reconstructed by translating and changing the typhoon intensity, and we simulated the process of change in the wave current field during the typhoon, obtained the statistics of different paths and intensity of the typhoon transit process of the maximum effective wave height and current velocity in the Maniao Bay aquaculture area, which were used as an assessment of the strength of the influence of the Maniao Bay nets aquaculture area, the strength of the impacts of the factors. Moreover, we analyzed theimpacts of different paths and typhoon intensities on the net box aquaculture area, and proposed the range of warning paths during the transit of typhoons of different intensities. The results show that among the four types of typhoon paths, the first type of typhoon paths had the widest and greatest impacts on the net-pen aquaculture area, in which the maximum effective wave height in the net-pen aquaculture area under the R4 path was higher than 5.5 m, and the maximum current velocity was more than 1.5 m·s−1 under the 'super typhoon' category, which severely affected the aquaculture area; the second and third types of typhoon paths had the greatest impacts on the net-pen aquaculture area. Most of the paths of Category II and III typhoons had less moderate impacts on the net-pen culture area; and the paths of Category IV typhoons had basically no impacts on the net-pen culture area under the intensity of each typhoon.
-
Keywords:
- Typhoon wave /
- Aquaculture cages /
- Early warning assessment /
- Delft3D /
- Wave-current coupling model
-
氯霉素(chloramphenicol，CAP)是一种具有旋光活性的高效广谱抗菌药，农业生产上可以有效治疗家畜家禽及水产生物所患的由致病性革兰阴性菌引起的疾病[1-3]。然而，氯霉素可以引起一系列严重的“三致”毒性和临床不良反应[4]。氯霉素对实验动物的致癌性证据充分，而对人体的致癌性证据有限。因此世界卫生组织(World Health Organization，WHO)中的国际癌症研究机构(International Agency for Research on Cancer，IARC)在1990年将其归为2A类致癌物质(对人很可能致癌)[1,5]。此外，氯霉素可通过短期可逆性的骨髓抑制，引起各类血细胞的减少，如血小板减少性紫癜、粒细胞缺乏症等。严重时氯霉素可引起不可逆的再生障碍性贫血，患者造血功能极少恢复，虽然发生较晚但死亡率高，而且存活者也容易发生粒细胞性白血病。氯霉素的另一个严重不良反应是会引起灰婴综合征(gray baby syndrome，GBS)，摄入大剂量氯霉素后，血药浓度异常增高引起婴儿的循环衰竭[6]。
因此，我国农业部235号公告明确规定，氯霉素为禁止用于所有食品动物的药物，在所有食品动物中不得检出[7]。近年来，有媒体报道称在贝类产品中检出氯霉素，原因可能是运输者为了防止贝类在运输途中滋生细菌，人为非法添加氯霉素。针对该情况，笔者对2015—2017年市场、渔获码头贝类进行了跟踪抽样调查，采用GB/T 20756—2006液相色谱-串联质谱法，对贝类样品中氯霉素残留状况进行了检测，并采用国际食品法典委员会(CAC)推荐的遗传毒性致癌物危险性评估技术——暴露边界比（margin of exposure，MOE）法[8]，通过@risk软件对贝类中非法添加氯霉素可能引起的食品安全风险进行了评估。
目前，MOE法在风险评估中应用广泛[9-10]。MOE是对人或实验动物产生很小但可测量作用的剂量，与估计的人的暴露量之间的比率。遗传毒性的致癌物的MOE值通过该类物质的基准剂量最低限值(lower confidence limit on the benchmark dose，BMDL) 与估计人群暴露量的比值计算得到，BMDL是描述遗传毒性致癌物危害特征的起点。MOE值越小，说明该物质对人类的致癌危险性越大，反之就越小。王君和吴永宁[8]运用MOE法评价了遗传毒性致癌物丙烯酰胺的危险性；Thompson等[11]综合了机械学和药代动力学信息，得出了六价铬的参考剂量和MOE值；Okauru等[12]运用MOE法比较了多环芳烃和热水对食管癌发病的风险。本研究采用MOE法评价氯霉素的暴露量，可以很好地衡量相关风险。
@risk软件是一款基于蒙特卡罗模拟技术[13]加载到Microsoft Excel上进行风险评估和决策分析的专业软件[14]，该软件为Excel增添了高级模型和风险分析功能，允许在建立模型时应用各种概率分布函数，对各种可能出现的结果进行模拟，得出构成风险的各种事件的发生概率，并能够以各种图表展示分析结果[13,15-16]。
本研究首次采用蒙特卡罗模拟技术，基于@risk 7.5.2版本，通过对检测到的市售贝类产品中氯霉素质量分数数据进行分布拟合，并模拟抽样10 000次[13,16-17]，进而得到膳食暴露量的概率分布曲线及其中位数(P50)和第95百分位数(P95)的值，通过计算膳食摄入量和人群估计暴露量，对所得结果进行暴露边界比评价。属于定量的风险评估方法，得到的数据直观，可与其他类似风险进行横向比较。该方法可以较好地评估市售贝类产品中氯霉素的暴露风险，对该风险曾经发生、现在已基本消除的情况进行量化描述，以期说明通过加大执法力度进行风险管理可以消除该风险。
1. 材料与方法
1.1 实验材料
1.1.1 样品采集
2015—2017年，根据《水产品抽样规范》(GB/T 30891—2014)的抽样方法，在市场采集贝类样品，每个档口只采集1个贝类品种且只算作1个样品，样品取净肉400 g，共300个样品，包括东风螺(Babylonia areolata)、翡翠贻贝(Perna viridis)、文蛤(Meretrix meretrix)、菲律宾蛤仔(Ruditapes philippinarum)、杂色蛤仔(R. variegata)、波纹巴非蛤(Paphia undulata)等。具体采样情况见表1。所采样品包括蛤类、螺类、贻贝等，均为市场常见品种，能够一定程度上反映当地贝类的消费情况。
表 1 2015—2017年市售贝类产品采样情况Table 1. Sample collection of commercially available shellfish in 2015−2017年份
year样品种类
species数量
quantity2015 文蛤 Meretrix meretrix 23 菲律宾蛤仔 Ruditapes philippinarum 23 杂色蛤仔 R. variegata 33 波纹巴非蛤 Paphia undulata 10 东风螺 Babylonia areolata 11 小计 subtotal 100 2016 文蛤 M. meretrix 26 波纹巴非蛤 P. undulata 11 菲律宾蛤仔 R. philippinarum 8 翡翠贻贝 P. viridis 55 小计 subtotal 100 2017 文蛤 M. meretrix 15 杂色蛤仔 R. variegata 6 波纹巴非蛤 P. undulata 12 翡翠贻贝 P. viridis 67 小计 subtotal 100 总计 total 300 1.1.2 试剂耗材
氯霉素(Dr.Ehrenstorfer，99.9%)、氯霉素-d5溶液(Dr.Ehrenstorfer，100.0 μg·mL–1)、正己烷(Sigma，≥95%)、乙酸乙酯(Sigma，≥99.7%)、甲醇(Fisher Scientific，99.9%)、氢氧化铵溶液(Macklin，25%~28%)、无水硫酸钠(广州化学试剂厂，≥99.0%)
1.1.3 仪器
API 3000-Agilent 1200 Series液相色谱-串联质谱仪、METTLER TOLEDO XS204电子天平、Thermo Scientific BIOFUGE STRATOS离心机、IKA MS3 basic涡旋振荡器、IKA T25 digital组织匀浆机、BUCHI R-215旋转蒸发仪。
1.2 样品分析
根据《可食动物肌肉、肝脏和水产品中氯霉素、甲砜霉素、氟苯尼考残留量的测定 液相色谱-串联质谱法》(GB/T 20756—2006)的相关规定进行检测。测定中采用空白实验、平行双样、质控样品测定、阳性复测等进行内部质量控制。
1.3 数据分析
1.3.1 检测结果
采用标准曲线法得到供试液中氯霉素的体积分数，贝类产品中氯霉素的残留量X=cV/m，c为供试液中氯霉素的体积分数，V为溶解残余物所用体积，m为贝类样品的质量。
1.3.2 暴露评估
1) 贝类产品中氯霉素质量分数数据的分布拟合。采用@risk软件，将贝类产品中氯霉素质量分数的检测值进行分布拟合。参照世界卫生组织《食品中痕量污染物的可靠评价》[18]和美国环境保护署(USEPA)建议的数据处理方法[19-20]，对于未检出样品按方法最低检出浓度的1/2进行处理取值，本研究中氯霉素的方法检出限为0.1 μg·kg–1，故对于未检出样品按照0.05 μg·kg–1进行处理取值。将得到的不同拟合曲线的拟合度运用卡方检验、安德森-达林检验、柯尔莫可洛夫-斯米洛夫检验这3种统计检验方法进行检验，并综合考虑3种方法的结果，最终确定最佳拟合分布[16-17,21-22]。
2) 膳食摄入量。根据GEMS/Food在2012年发布的数据，我国所在的G09区(东亚及东南亚地区)，软体动物的平均摄入量为12.2 g·d–1 [23]，采用最保守的估计值10 g·d–1为贝类产品的膳食摄入量。本研究在软件中设定膳食摄入量为(10±3) g·d–1，μ=10 g·d–1，σ=1的随机变化正态分布函数。
3) 每人每日膳食平均暴露量 [以下简称膳食暴露量(estimated dietary exposure，EXP)]。计算公式为：
$${\rm{EX}}{{\rm{P}}_{\rm{D}}}{{ = C}} \times I_D$$ (1) 其中EXPD为日膳食暴露量，C为贝类产品中氯霉素的质量分数，ID为日膳食摄入量。根据以上公式，将相关参数(包括软件拟合的贝类产品中氯霉素质量分数的分布函数)输入@risk软件，随机从贝类产品中氯霉素质量分数分布中抽取数值计算膳食暴露量，每次模拟过程循环10 000次，得到膳食暴露量P50和P95的值。
1.3.3 风险特征描述
暴露边界比评估。采用CAC推荐的遗传毒性致癌物危险性评估技术——MOE法[8,10,24]对市售贝类产品中氯霉素残留进行风险特征描述。根据公式MOE=BMDL/EXP，结合动物及临床实验数据，将BMDL定为每千克体质量1 mg·d–1。
EXP计算公式为EXP=EXPD/mb，其中mb为标准人体质量，按照60 kg计算，膳食暴露量按照公式(1)得到的膳食暴露量的P50和P95取值。
通过得到的MOE的P50和P95值进行评价，参照加拿大卫生部对于具有遗传毒性致癌物的MOE结果判定方法进行判定，即MOE值<5 000、5 000~500 000以及>500 000时分别表示被评估的危害物属于高、中、低危害程度[7]。
2. 结果
2.1 检测结果
所有样品均采用国标方法检测，经仪器分析的数据见表2。其中2015年的检出率为33%，2016年的检出率为3%，2017年的全部100个样品均未检出氯霉素残留，故在本研究中不继续计算2017年氯霉素的膳食暴露量。
表 2 2015—2017年市售贝类产品中氯霉素残留检测结果Table 2. Result of determination of chloramphenicol residues in commercially available shellfish in 2015−2017年份
year采集样品数
quantity of collected samples阳性样品数
quantity of positive samples检出值/μg·kg–1
range of concentration2015 100 33 0.36~4 050 2016 100 3 50.4~820 2017 100 0 – 总计 total 300 36 0.36~4 050 2.2 膳食暴露量
2.2.1 贝类产品中氯霉素质量分数数据的分布拟合
根据2.1的检测结果，在软件中分别进行2015年、2016年数据的分布拟合，经3种检验综合评价，数据更符合正态分布，得到最优的拟合曲线，分别为：
2015年，RiskNormal (121.43, 524.71, RiskTruncate(0, ))
2016年，RiskNormal (9.412 5, 82.293, RiskTruncate(0, ))
2.2.2 暴露评估结果
根据公式 (1)，将得到的@risk相关参数(表3)输入软件中。
表 3 氯霉素的膳食暴露量分析@risk软件计算参数Table 3. @risk parameter to calculate dietary exposure of chloramphenicol描述
description分布
distribution2015年氯霉素残留量 chloramphenicol residue, 2015 RiskNormal (121.43, 524.71, RiskTruncate(0, )) 2015年膳食摄入量 daily intake, 2015 RiskNormal (0.01, 1, RiskTruncate(0.007, 0.013)) 2015年膳食暴露量 dietary exposure, 2015 RiskOutput()+B1*B2 2016年氯霉素残留量 chloramphenicol residue, 2016 RiskNormal (9.4125, 82.293, RiskTruncate(0, )) 2016年膳食摄入量 daily intake, 2016 RiskNormal (0.01, 1, RiskTruncate(0.007, 0.013)) 2016年膳食暴露量 dietary exposure, 2016 RiskOutput()+B4*B5 在@risk软件中进行蒙特卡罗模拟，每次模拟过程循环10 000次，得到的膳食暴露量的概率分布曲线见图1。
![]()
图 1 氯霉素的膳食暴露量概率分布曲线Figure 1. Probability distribution curve of dietary exposure of chloramphenicol根据得到的氯霉素的膳食暴露量概率分布曲线，分别得到对应年份膳食暴露量的P50和P95 (表4)。
表 4 氯霉素的膳食暴露量Table 4. Dietary exposure of chloramphenicol年份
year膳食暴露量/μg·kg–1 dietary exposure 概率分布中位百分数
P50概率分布第95百分位数
P952015 3.91 11.5 2016 0.579 1.73 2017 – – 2.3 MOE
本研究采用单一食品选择研究法，定量估计的分析方式，通过计算得到的暴露量以及MOE值见表5。2015年居民消费贝类产品所致氯霉素的平均暴露量为每千克体质量0.065 1 μg·d–1，计算得到的MOE值为15 361，介于5 000~500 000，按照加拿大卫生部的标准，属于中危害水平。高暴露量人群的暴露量为每千克体质量0.192 μg·d–1，计算得到的MOE值为5 208，介于5 000~500 000，按照加拿大卫生部的标准，属于中危害水平。2016年居民消费贝类产品所致氯霉素的平均暴露量为每千克体质量0.009 65 μg·d–1，计算得到的MOE值为103 627，介于5 000~500 000，按照加拿大卫生部的标准，属于中危害水平。高暴露量人群的暴露量为每千克体质量0.028 8 μg·d–1，计算得到的MOE值为34 722，介于5 000~500 000，按照加拿大卫生部的标准，属于中危害水平。2017年由于在市售贝类产品中未检出氯霉素，属于无危害水平。
表 5 氯霉素的人群估计暴露量及暴露边界比Table 5. Estimated dietary exposure and margin of exposure of chloramphenicol年份
year人群估计暴露量/μg·d–1
EXP暴露边界比
MOE概率分布中位百分数
P50概率分布第95百分位数
P95概率分布中位百分数
P50概率分布第95百分位数
P952015 0.065 1 0.192 15 361 5 208 2016 0.009 65 0.028 8 103 627 34 722 2017 – – – – 3. 讨论
3.1 氯霉素BMDL的取值
目前的致癌理论认为，遗传毒性致癌物是一类通过诱发体细胞基因突变从而激活致癌基因或者灭活抑癌基因的表达从而诱发癌变的致癌物，在体细胞基因中几个分子甚至一个分子的基因突变就有可能诱发癌症[8]。氯霉素的致癌性也有类似特征，给其危害的评估带来一定困难。对市售贝类产品中氯霉素残留进行暴露评估，是一种行之有效的方法。
在致畸性研究中，氯霉素未表现出对大鼠、兔的致畸作用，但即使小剂量给药，均可引起高死胎率。几乎所有氯霉素致突变研究，均可得到阳性结果。氯霉素在离体细胞实验中可引发DNA断裂、染色体畸变并增加姐妹染色单体互换[25]。对于具有遗传毒性和致癌性的物质，传统假设它们没有阈值，并在任意暴露水平都会存在一定风险。故氯霉素无每日允许摄入量(acceptable daily intake，ADI)值，在世界粮农组织/世界卫生组织的报告中同样认为其无ADI值[3]。
BMDL的取值通常根据实验动物及临床实验数据得到，根据WHO的国际化学品安全规划(International Programme on Chemical Safety，IPCS)提出的“化学品危险性评估的剂量-反应分析建模的原则”[9]，风险评估时可通过实验数据的拟合分析得到基准剂量(benchmark dose，BMD)，进而得到BMDL值。由于氯霉素尚缺少BMDL值，需要根据各方面的实验数据综合考虑。研究表明，小鼠单次口服每千克体质量50 μg·d–1和100 μg·d–1的氯霉素可以引起染色体畸变[8]。成人或较大儿童应用氯霉素剂量大于每千克体质量100 μg·d–1，早产儿或新生儿大于每千克体质量25 μg·d–1可能引起灰婴综合征[2]，氯霉素注射液的临床用量为成人每日0.5~1 g[6]，按60 kg标准人体质量折算后为每千克体质量8~17 μg·d–1。以上剂量均为短期给药剂量，综合考虑长期食用的影响，本研究采用的BMDL为每千克体质量0.1 μg·d–1。该取值是否合理，仍需要更多数据进行评估。
3.2 氯霉素的限量值及最低执行限量要求
我国农业部235号公告明确规定氯霉素为禁用药物[3]。早在1969年，世界卫生组织(WHO)食品添加剂联合专家委员会(JECFA)和联合国粮农组织(FAO)第12届会议就建议禁止使用氯霉素[5,22]。在欧盟规定(EEC) No.2377/90以及以后的各版本限量标准中，氯霉素被列在不制定残留限量标准(禁用)的物质中[11]。氯霉素在各个国家均无最高残留限量值(maximum residue limit，MRL)。
我国相关水产品检验检测标准中，由于采用的方法不同，氯霉素的检出限在0.1~0.3 μg·kg–1。欧盟委员会2003/181/EC决议的官方公报[26]以及全球水产养殖联盟(Global Aquaculture Alliance，GAA)制定的最佳水产养殖规范(Best Aquaculture Practice，BAP)中关于水产品中氯霉素残留的最低执行限量要求(minimum required performance limits，MRPL)为0.3 μg·kg–1。
本研究氯霉素的检出限为0.1 μg·kg–1，符合相关规定和检测要求；同时对于未检出的样品采用检出限的1/2进行处理取值符合相关指南的要求。
3.3 膳食摄入量
相对其他海产品，我国消费者人群贝类的膳食摄入量数据在文献报道中较少涉及。根据GEMS/Food在2012年发布的数据，我国所在的G09区，软体动物(包括贝类和头足类)的平均摄入量为12.2 g·d–1，该区淡水鱼的平均摄入量为27.0 g·d–1，海水鱼的平均摄入量为1.29 g·d–1，甲壳类的平均摄入量为5.29 g·d–1 [23]。《2002年中国居民营养与状况调查报告》中只统计了鱼类及虾蟹类的数据，数据显示我国居民平均每标准人日鱼类摄入量为24.8 g、虾蟹类的摄入量为4.8 g[27]。《广东省居民膳食营养与健康状况十年变化分析》(2002—2012)中只统计了鱼虾类的数据，每标准人日鱼虾类摄入量为56.8 g[28]。《中国居民营养与健康状况监测2010—2013年综合报告》的数据也只统计了鱼虾类的数据，城乡居民平均每标准人日鱼虾类摄入量为23.7 g，与2002年相比下降了5.9 g[29]。唐洪磊等[30]2006年的调查显示，广东省沿海城市，虾蟹贝类居民日均摄入量中值为24.8 g·d–1，均值为66.6 g·d–1。综合以上数据，本研究将10 g·d–1定为估计的我国居民贝类产品的膳食摄入量。
3.4 MOE法及其安全边界值
MOE分析技术是欧洲食品安全局在2005年FAO/WHO JECFA第64次会议上提出的，是评估遗传毒性致癌物危险性的有效方法。目前，虽然对于MOE值的分界值界定，国际社会尚没有一个公认值，不同国家及其机构的判定标准不一，但均可作为判定物质安全性的参考。加拿大卫生部对于MOE值分别界定了5 000和500 000两个分界值，分别界定出被评估危害物质危害程度的高、中、低[8]。英国癌症委员会(COC)的一个草案认为，MOE值达到10 000时可能已经足够安全。欧盟也认为MOE值达到10 000或以上，风险即较低[8,10,31]。本研究综合以上标准，采用较严格的加拿大卫生部标准，可以较好地区分不同年度市售贝类产品中氯霉素残留风险的严重程度。
3.5 贝类产品中氯霉素残留对消费者健康的风险评估
在2015年之前，市售贝类产品中人为违规添加氯霉素的情况均无人报道，其严重程度亦无从考量。自从2015年有媒体报道相关问题以来，政府部门便采取了一系列行动，以期能有效地控制该问题。本研究在2017年的采样检测时发现，市售贝类产品中氯霉素残留量已低于检测限，证明该问题已基本得到有效控制，对消费者的健康危害低。同时，监管部门仍应持续保持对此类风险问题的关注，以防死灰复燃。
-
![]()
图 1 海南省马袅湾网箱位置示意图
Figure 1. Geographic location, distribution of cages in Maniao Bay, Hainan Province
![]()
图 5 “威马逊”“启德”和“海鸥”台风有效波高结果比较
Figure 5. Verification of significant wave height results for typhoons Rammasun, Kaitai and Kalmaegi
![]()
图 7 不同台风路径下各等级台风发生次数
Figure 7. Number of various levels of typhoons under different tracks
![]()
图 8 “威马逊”台风下养殖区有效波高最大值分布及最大流速时刻流场
注：a. 养殖区有效波高最大值分布；b. 养殖区内最大流速时刻流场。
Figure 8. Distribution of maximum significant wave height and maximum flow velocity moment flow field in aquaculture area under typhoon Rammasun
Note: a. Distribution of maximum significant wave height in the aquaculture area; b. Maximum flow velocity moment flow field in the aquaculture area.
![]()
图 10 马袅湾养殖区内有效波高及流速最大值
Figure 10. Maximum significant wave height and current velocity in aquaculture area of Maniao Bay
![]()
图 11 “圆规”台风下养殖区有效波高最大值分布及最大流速时刻流场
注：a. 养殖区有效波高最大值分布；b. 养殖区内最大流速时刻流场。
Figure 11. Distribution of maximum significant wave height and maximum flow velocity moment flow field in aquaculture area under typhoon Kompasu
Note: a. Distribution of maximum significant wave height in the aquaculture area; b. Maximum flow velocity moment flow field in the aquaculture area.
![]()
图 13 马袅湾养殖区内有效波高及流速最大值
Figure 13. Maximum significant wave height and current velocity in aquaculture area of Maniao Ba
![]()
图 14 “黄蜂”台风下养殖区有效波高最大值分布及最大流速时刻流场
注：a. 养殖区有效波高最大值分布；b. 养殖区内最大流速时刻流场。
Figure 14. Distribution of maximum significant wave height and maximum flow velocity moment flow field in aquaculture area under typhoon Vongfong
Note: a. Distribution of maximum significant wave height in the aquaculture area; b. Maximum flow velocity moment flow field in the aquaculture area.
![]()
图 16 马袅湾养殖区内有效波高及流速最大值
Figure 16. Maximum significant wave height and currentvelocity in aquaculture area of Maniao Bay
![]()
图 17 “泰德”台风下养殖区有效波高最大值分布及最大流速时刻流场
注：a. 养殖区有效波高最大值分布；b. 养殖区内最大流速时刻流场。
Figure 17. Distribution of maximum significant wave height and maximum flow velocity moment flow field in aquaculture area under typhoon Ted
Note: a. Distribution of maximum significant wave height in the aquaculture area; b. Maximum flow velocity moment flow field in the aquaculture area.
![]()
图 19 马袅湾养殖区内有效波高及流速最大值
Figure 19. Maximum significant wave height and currentvelocity in aquaculture area of Maniao Bay
表 1 天文潮模拟精度
Table 1 Simulation accuracy of astronomical tide
观测站
Observation station平均绝对误差
MAE/m均方根误差
RMSE/m相关系数
CC秀英 Xiuying 0.134 0.166 0.978 湛江 Zhanjiang 0.109 0.145 0.990 
下载: 导出CSV
表 2 台风浪模拟精度
Table 2 Significant wave height accuracy
观测站
Observation station平均绝对
误差
MAE/m均方根
误差
RMSE/m平均绝对
百分比误差
MAPE/%相关系数
Correlation
coefficient秀英 Xiuying 0.179 0.255 7.4 0.992 浮标测点Buoy2 0.425 0.520 18.3 0.980 东方 Dongfang 0.238 0.335 14.6 0.934
下载: 导出CSV
表 3 台风等级相关参数信息
Table 3 Information on parameters related to typhoon levels
台风等级
Typhoon intensity
level最大风速
Maximum wind
speed/(m·s−1)中心气压
Central
pressure/hpa超强台风 Super Typhoon 51.0 944 强台风 Severe Typhoon 41.5 962 台风 Typhoon 32.7 977
下载: 导出CSV
-
[1] 赵广生, 牛小静. 海南岛历史风暴潮模拟和灾害风险评估[J]. 水动力学研究与进展A辑, 2022, 37(6): 831-836. [2] SHI X W, ZHENG S X, LIU Q, et al. Research on numerical simulation of typhoon waves with different return periods in nearshore areas: case study of Guishan island waters in Guangdong province, China[J]. Stoch Environ Res Risk Assess, 2021, 35: 1771-1781. doi: 10.1007/s00477-020-01960-4
[3] 姜雨青, 张俊波, 国志兴, 等. 我国沿海渔业养殖设施空间分布与台风影响关联性研究[J]. 海洋湖沼通报, 2023, 45(1): 183-192. [4] 石振强. 深水网箱浮圈水动力特性理论及试验研究[D]. 哈尔滨: 哈尔滨工程大学, 2015: 1-2. [5] 黄小华, 郭根喜, 陶启友, 等. HDPE圆形重力式网箱受力变形特性的数值模拟[J]. 南方水产科学, 2013, 9(5): 126-131. [6] 刘春宏. 海水网箱养殖区周围流场特性研究[D]. 天津: 天津大学, 2020: 34-35. [7] 姜雨青. 台风对我国沿海设施渔业影响的关联性及灾害评估[D]. 上海: 上海海洋大学, 2021: 53-55. [8] 黄滨, 关长涛, 崔勇, 等. 台风“米雷”对山东网箱养殖业灾害性影响的调查与技术解析[J]. 渔业现代化, 2011, 38(4): 17-21. [9] 张润川. 台风灾害链网络构建及其演化特征研究[D]. 福州: 福建师范大学, 2023: 34-44. [10] 吴亚男, 左军成, 吕晓凤, 等. 舟山海域风暴潮增水数值分析: 台风“灿鸿”为例[J]. 海洋湖沼通报, 2017(6): 38-45. [11] 原楠, 陈新平, 陈学恩, 等. 罗源湾海域台风风暴潮数值模拟研究[J]. 海洋通报, 2019, 38(1): 20-30. [12] 罗志发, 黄本胜, 邱静, 等. 粤港澳大湾区风暴潮时空分布特征及影响因素[J]. 水资源保护, 2022, 38(3): 72-79,153. [13] BI C, ZHAO Y, SUN X, et al. An efficient artificial neural network model to predict the structural failure of high-density polyethylene offshore net cages in typhoon waves[J]. Ocean Engin, 2020, 196: 106793. doi: 10.1016/j.oceaneng.2019.106793
[14] ZHANG Y, GUO H, LIU S, et al. Fatigue vulnerability of sea cage to storm wave loads[J]. J Mar Sci Technol, 2022, 28(1): 153-164.
[15] LI L, FU S, XU Y. Nonlinear hydro-elastic analysis of an aquaculture fish cage in irregular waves[J]. Mar Struc, 2013, 34: 56-73. doi: 10.1016/j.marstruc.2013.08.002
[16] YIN C, HUANG H, WANG D, et al. The characteristics of storm wave behavior and its effect on cage cultureusing the ADCIRC+SWAN model in Houshui Bay, China[J]. J Ocean Univ China, 2020, 19(2): 307-319. doi: 10.1007/s11802-020-3941-3
[17] 孙熊雄. 基于机器学习的重力式深水网箱受灾风险评估分析研究[D]. 大连: 大连理工大学, 2019: 26-31. [18] 刘圣聪, 黄六一, 倪益, 等. 水下圆形重力式网箱缆绳张力对波、流的响应特性[J]. 水产学报, 2019, 43(6): 1518-1526. [19] 孟范兵. 波流作用下深水重力式网箱水动力特性研究[D]. 大连理工大学, 2014: 44-53. [20] 黄婉茹, 郭敬, 谭骏, 等. 风暴潮灾害承灾体风险预警技术研究[J]. 自然灾害学报, 2023, 32(4): 85-93. [21] 于松辰. 波流作用下重力式网箱动力响应研究[D]. 哈尔滨: 哈尔滨工程大学, 2019: 1-3. [22] 朱磊, 刘会欣. 基于Delft3D模型的风暴潮增减水模拟研究: “9711”号台风为例[J]. 海洋湖沼通报, 2018(5): 1-10. [23] 孙志林, 王辰, 钟汕虹, 等. 浪潮耦合的舟山渔港台风暴潮数值模拟[J]. 海洋通报, 2019, 38(2): 150-158. doi: 10.11840/j.issn.1001-6392.2019.02.004 [24] HU K, MESELHE E, RHODE R, et al. The impact of levee openings on storm surge: a numerical analysis in Coastal Louisiana[J]. Appl Sci, 2022, 12(21): 10884. doi: 10.3390/app122110884
[25] MCLAUGHLIN C, LAW B, MULLIGAN R. Modeling surface waves and tide-surge interactions leading to enhanced total water levels in a macrotidal bay[J]. Coast Engin J, 2022, 64(1): 24-41. doi: 10.1080/21664250.2021.1965417
[26] FUJITA T. Pressure distribution within typhoon[J]. Geophys Mag, 1952, 23: 437-451.
[27] MYERS V A. Maximum hurricane winds[J]. Bull Amer Meteorol Soc, 1957, 38: 227-228.
[28] JELESNIANSKI C P. A numerical calculation of storm tidesinduced by a tropical storm impinging on a continental shelf[J]. Mon Weather Rev, 1965, 93(6): 343-358. doi: 10.1175/1520-0493(1993)093<0343:ANCOS>2.3.CO;2
[29] HOLLAND G J. An analytic model of the wind and pressure profiles in hurricanes[J]. Mon Weather Rev, 1980, 108(8): 1212-1218. doi: 10.1175/1520-0493(1980)108<1212:AAMOTW>2.0.CO;2
[30] VICKERY P J, SKERLJ P F, STECKLEY A C, et al. Hurricane wind field model for use in hurricane simulations[J]. J Struct Eng, 2000, 126(10): 1203-1221. doi: 10.1061/(ASCE)0733-9445(2000)126:10(1203)
[31] WILLOUGHBY H E, RAHN M E. Parametric representation of the primary hurricane vortex. part I: observations and evaluation of the Holland(1980) model[J]. Mon Weather Rev, 2004, 132(12): 3033-3048. doi: 10.1175/MWR2831.1
[32] MING Y, ZHANG W, YU H, et al. An overview of the China meteorological administration tropical cyclone database[J]. J Atmos Ocean Technol, 2014, 31(2): 287-301. doi: 10.1175/JTECH-D-12-00119.1
[33] LU X, YU H, YING M, et al. Western north pacific tropical cyclone database created by the China meteorological administration[J]. Adv Atmos Sci, 2021, 38(4): 1-10.
[34] PHILLIPS N A. A coordinate system having some special advantages for numerical forecasting[J]. J Meteorol Res, 1957, 14: 184-185. doi: 10.1175/1520-0469(1957)014<0184:ACSHSS>2.0.CO;2
[35] BOOIJ N, RIS R C, HOLTHUIJSEN L H. A third-generation wave model for coastal regions 1. Model description and validation[J]. J Geophy Res Atmos, 1999, 104(4): 7649-7666.
[36] WU G, SHI F, KIRBY J T, et al. Modeling wave effects on storm surge and coastal inundation[J]. Coast Engin, 2018, 140: 371-382. doi: 10.1016/j.coastaleng.2018.08.011
[37] 徐亚男. 风暴潮与波浪耦合数值预报模型的研究[D]. 天津: 天津大学, 2012: 26-27. [38] PINDSOO K, SOOMERE T. Contribution of wave set-up into the total water level in the Tallinn area[J]. Proc Est Acad Sci, 2015, 64(3): 338-348. doi: 10.3176/proc.2015.3S.03
[39] 薛鸿超, 顾家龙, 任汝述. 海岸动力学[M]. 北京: 人民交通出版社, 1980: 129-167. [40] BROWN M J, WOLF J. Coupled wave and surge modelling for the eastern Irish Sea and implications for model wind-stress[J]. Cont Shelf Res, 2009, 29(10): 1329-1342. doi: 10.1016/j.csr.2009.03.004
[41] 黄山. 琼州海峡风暴潮数值模拟与不同台风路径对增水的影响研究[D]. 广州: 华南理工大学, 2016: 31-32. [42] 石青. 近岸台风浪及植被水流环境中的风浪数值模拟研究[D]. 大连: 大连理工大学, 2020: 39-42. [43] 尹超. 海南岛近海热带气旋引起的台风浪风险评估研究[D]. 青岛: 中国科学院大学(中国科学院海洋研究所), 2020: 43-45. -
期刊类型引用(6)
1. 罗兰，王甫，黄承鹏，秦威振，任燕，林启辉. 深圳市市售动物源性食品中氯霉素残留量暴露评估. 中国食品卫生杂志. 2023(04): 581-586 . 
百度学术
2. 孙蕾. 市售水产品氯霉素残留污染状况及健康风险评估. 现代食品. 2022(22): 171-173 . 百度学术
3. 田甜，巫剑，文金华，曾祥林. 广西北部湾鲜活水产品中氯霉素残留的膳食暴露风险评估. 现代食品科技. 2021(05): 310-318+286 . 百度学术
4. 王胜，孙言凤，肖永华，黄常刚. 2011—2020年W市市售水产品中4种违禁渔药污染状况及暴露风险评估. 食品安全质量检测学报. 2021(24): 9612-9618 . 百度学术
5. 翟毓秀，郭萌萌，江艳华，姚琳，赵艳芳，吴海燕，李风铃，谭志军. 贝类产品质量安全风险分析. 中国渔业质量与标准. 2020(04): 1-25 . 百度学术
6. 黄思聪，龙敏仪. QuEChERS-液相色谱-串联质谱法测定牛蛙肉中酰胺醇类抗生素及其代谢物. 肉类研究. 2019(06): 39-43 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 0
- HTML全文浏览量: 0
- PDF下载量: 0
- 被引次数: 7
粤公网安备 44010502001741号
