南黄海苏北浅滩浒苔体内的农药污染与食用风险评价

杨丰源 ，和 庆 ，陈以芹 ，蔡春尔 ，陈彦百 ，李娟英 *

（1.上海海洋大学水产动物遗传育种中心协同创新中心，上海 201306；2.上海海洋大学海洋生态与环境学院，上海 201306）

自2007年以来，我国南黄海近海海域连续10年爆发浒苔形成的大规模绿潮，覆盖海域面积最大超过1000 km2，规模堪称世界之最[1]。绿潮爆发对黄海近海海域生态环境和功能造成严重破坏，也造成巨大经济损失[2]。因此，解决绿潮爆发带来的危害，探索浒苔资源化利用的有效途径，对系统解决海洋绿潮带来的环境、生态和经济问题具有重要意义。研究表明，浒苔是重要的食品、饲料、医药和新能源材料[3-8]。其中，食用既能解决浒苔的爆发问题，又可以提供有营养的海洋藻类。浒苔是一种高蛋白（26.46%～29.01%）、低脂肪（0.76%～0.91%）和富含纤维（4.87%～6.58%）的海藻食品，并且含有丰富的氨基酸（17种）、脂肪酸（13种）和多种矿物质（Ca、Fe和Zn等）[5-6]，其中铁含量在中国食物营养成分表上记载为我国食物之最。如今，在日本及我国福建南部、江苏和浙江等地浒苔已成为一种大受欢迎的食品[9]。目前所使用的农药种类大致可分为有机磷类、有机氯类、拟除虫菊酯类和氨基甲酸酯类等。虽然我国从1983年开始禁止使用有机氯农药（OCPs），很多地区环境中OCPs污染也逐年减轻，但仍有部分残留。有机磷（OPPs）和拟除虫菊酯类农药（PEs）是继OCPs之后两类大量生产和使用的农药[10]，调查研究发现，南黄海沉积物中的OCPs范围是0.5～7.6 ng·g-1[11]；浒苔爆发的另一区域山东省近海胶州湾海域，水体中 OPPs的含量在 0.2～77.6 ng·L-1[12]，沉积物中OCPs的含量为7.7～213 ng·g-1[13]；该区域内的条斑紫菜、坛紫菜、海带、裙带菜和龙须菜等5种海藻体内的 OCPs浓度为 18.13～60.3 ng·g-1[14]。由此可见，海域环境中的农药可能导致浒苔体内某种程度的残留，从而影响浒苔的质量和食用安全，并给后续消费浒苔的人群带来一定的健康风险。

因此本文针对南黄海苏北浅滩地区的浒苔进行了农药残留的检测，选择在浒苔爆发严重的南黄海苏北浅滩地区如东、大丰和吕四等地进行研究。该地区的浒苔种类包括肠浒苔、曲浒苔、缘管浒苔、盘浒苔和浒苔属浒苔等5种。采用微波萃取-气相色谱质谱联用法，对浒苔中8种OCPs、4种PEs和9种OPPs农药进行了检测，并在此基础上对浒苔食用的健康风险进行评价，结果不仅为南黄海地区浒苔能否食用提供了数据支撑，同时也可以为该地区农药的科学使用提供理论指导。

1 材料与方法

1.1 试剂

无水硫酸钠、Florisil硅土（100～200目）、层析硅胶（80～100目）和氧化铝（均购自国药集团，使用前分别于马弗炉中在450、650、650℃和400℃下干燥6 h，其中层析硅胶和氧化铝在冷却一段时间后加3%蒸馏水混匀脱活，所有试剂在干燥室中冷却至室温，储存于干燥器内待用）；色谱纯试剂：甲醇（赛默飞世尔有限公司）、乙酸乙酯（国药集团）、丙酮（默克股份两合公司）、二氯甲烷（国药集团）和正己烷（上海Sigma公司）等有机溶剂；分析纯试剂：无水硫酸钠、硅胶、Florisil硅土和铜片等（均购自国药集团）。8种OCPs标液包括 α-HCH、β-HCH、γ-HCH、δ-HCH、p，p'-DDD、p，p'-DDE、六氯苯和环氧七氯，购自德国Dr.Ehrenstorfer公司；4种PEs标液包括氯氟氰菊酯、联苯菊酯、氯菊酯和甲氰菊酯，购自国家标准物质研究中心；9种OPPs标液包括乐果、敌敌畏、马拉硫磷、甲基对硫磷、杀螟硫磷、三唑磷、喹硫磷、伏杀硫磷和毒死蜱，购自安普实验科技股份有限公司，且标液均在有效期内。

1.2 研究区域

苏北浅滩养殖区是早期爆发绿潮的重要区域[1]，本研究区域（图1）主要位于我国江苏省沿海大丰、如东和吕四境内，该区域北临黄海，东接东海，属于苏北浅滩紫菜养殖区域，养殖面积可达上万公顷，根据紫菜养殖区面积大小进行细化，其中，如东0.82万hm2，大丰大约0.67万hm2，吕四大约0.07万hm2，因此，将如东和大丰的采样点分成RD1、RD2和RD3以及DF1、DF2、DF3和DF4多个样点。浒苔是附着生长在紫菜养殖筏架上的一种大型绿藻，大量养殖筏架为苏北浅滩区的浒苔提供了重要的附着基，故该区域是浒苔的主要生长区，浒苔在冬季附着在紫菜养殖筏架上，随着温度的上升，浒苔开始大量繁殖，可达至数千吨[15]。6、7月份是浒苔大量繁殖的季节。

图1 采样点示意图Figure1 Sampling sites schematic

1.3 样品采集

于2016年1—8月分别在江苏如东、大丰、吕四等紫菜养殖筏架进行浒苔的采集。所采浒苔样品由科研人员从养殖筏架上混合采集，并置于加入冰袋的保温箱中运回实验室，储存在-80℃冰箱中保存；实验同时采集了养殖筏架下的沉积物样品。每个养殖筏架按对角线布点法用彼得森采泥器（DXCN1/40）采集多个分点的表层沉积物，混匀装入铝罐中，放置于加入冰袋的保温箱中运回实验室后于-20℃冰箱保存；48 h内完成浒苔样品的前处理，7 d内完成样品分析。

1.4 样品前处理

生物和沉积物样品中OCPs、PEs和OPPs的萃取、净化及测定方法参见参考文献[16-18]，简述如下，首先，分别称取冷冻干燥后的3 g沉积物和1 g浒苔样品，每个样品2个平行，用二氯甲烷溶液对样品进行微波萃取，之后合并萃取液，向萃取液中加入铜片去除硫化物，然后将萃取液全部转移至旋转蒸发瓶，浓缩至1～2 mL后用石油醚进行脱脂。最后，浓缩液通过预淋洗过的净化柱净化，再用洗脱液洗脱，OCPs、PEs和OPPs的洗脱液分别是正己烷、正己烷/二氯甲烷（体积比7∶3）和丙酮/二氯甲烷（体积比1∶1），最后利用旋转蒸发仪将收集到的溶液蒸发至近干，加入正己烷定容至 500 μL，并抽取 100 μL到内插管中待GC-MS分析。

1.5 气质分析条件

采用 GC-MS（Agilent 7890A/5975C）和 DB-5MS（30 m×0.25 mm×0.25 μm）色谱柱进行分离分析；以流速为1.2 mL·min-1的氦气（99.99%）为载气；进样方式为脉冲不分流进样，进样量为1 μL。气相条件：进样口温度为280℃；色谱柱升温程序：初温50℃，以25℃·min-1升到100℃，再以5℃·min-1升到300℃，保留5 min；质谱条件：离子源为EI源，能量为70 eV，离子源温度为230℃，传输线温度为300℃，质量分析器为四极杆，四极杆温度为150℃，扫描方式为选择离子扫描（SIM），根据标准色谱图中各组分的保留时间，对样品中的各组分进行定性分析，采用单点外标法峰面积计算法进行定量分析，农药中目标化合物的定性、定量离子见表1。

1.6 质量控制

实验前对浒苔及沉积物样品进行了OCPs、PEs和OPPs的回收率实验，其回收率在85%～105%之间，相对标准偏差在6%～10%，浒苔方法检出限分别为 1.5、2.04 μg·kg-1和 1.43 μg·kg-1，沉积物方法检出限分别为 1.5、0.40 μg·kg-1和 0.50 μg·kg-1。为保证实验结果准确性和科学性，浒苔和沉积物中农药含量测定实验设定了1个空白样和2个平行样。

表1 GC-MS分析农药中目标化合物的定性、定量离子Table1 Quantification and identification ions for the GC-MS analysis of the target compounds

1.7 食用风险评价

1.7.1 农药日摄入量（Estimated Daily Intake，EDI）

市民因食用浒苔摄入的三类农药含量的计算表达式如下：

式中：EDI代表农药日摄入量，ng·kg-1·d-1（dw）；c代表浒苔中农药残留，ng·g-1（dw）；Con代表市民每日平均浒苔消耗量，本文Con以0.182 kg·周-1计算，即0.026 kg·d-1[19]（由于人们日常食用浒苔较少，其消费量难以确定，但明显低于其他水产品的消费量，为了保守评估，此处参照鱼虾等其他食用水产品的摄入量0.026 kg·d-1，作为浒苔的日摄入量）；Bw代表市民平均体重，约为70 kg[20]。

1.7.2 目标危险系数（Target Hazard Quotient，THQ）

THQ是以测定的人体摄入剂量与参考剂量的比值为评价标准，适用于单一污染物评价，假设人体摄入剂量等于吸收剂量。

式中：ADI（Acceptable Daily Intake）代表每日允许摄入量，ng·kg-1·d-1；参考 GB 2763—2014《食品中农药最大残留限量》进行计算，其中，α-六六六、δ-六六六、联苯菊酯、甲氰菊酯、高效氯氟氰菊酯、氯菊酯、乐果、甲基对硫磷、杀螟硫磷、马拉硫磷、毒死蜱、喹硫磷和三唑磷分别为 5×103、5×103、1×104、3×104、4×104、5×104、2×103、3×103、6×103、3×105、1×104、5×102ng·kg-1·d-1和 1×103ng·kg-1·d-1。

若 THQ＜1，则无显著健康风险；若 THQ＞1，暴露人群有明显健康风险。THQ值越大，相应的风险越大。

1.7.3 风险指数（Hazard Index，HI）

HI用于评价复合污染的健康风险，公式如下：

式中：HI代表风险指数，若HI＜1，表明复合农药污染无显著健康风险；若HI＞1，则存在明显健康风险。HI值越大，相应的风险越大。

1.7.4 风险商

根据农业部颁发的《食品中农药残留风险评估指南》对苏北南黄海浅滩地区的浒苔中农药残留进行风险评估，结果依据GB 2763—2014《食品中农药最大残留限量》进行判定，其中最大残留限量值（MRL）的指标按国标判定，无MRL值的结果不判定。农药残留风险评估分短期（急性）和长期（慢性）摄入，其中短期摄入量（NESTI）和长期摄入量（NEDI）计算依据联合国粮食及农业组织（FAO）和世界卫生组织（WHO）联合管理的三个专家咨询机构之一的农药残留联席会议（JMPR）建议方法[21]，引入短期参考剂量（ARfD）和长期参考剂量（ADI）值判定风险商。

本研究中采集到的浒苔样品个体质量3～5 g之间，中位值3.24 g，即浒苔个体质量采用3.24 g（＜25 g）；LP是摄入量；Bw是消费者平均体重；Cm是农药残留量的中值或平均值，本研究平均值高于中位值，参照最严风险评估原则，计算中选取平均值[22]。借鉴欧盟（EU）推荐的评估方法和技术，HQ代表风险商，表示其风险程度，HQ＜1表示摄入食品是风险可接受，风险较低，HQ值越大风险越高[23]。

2 结果与讨论

2.1 南黄海苏北浅滩浒苔中农药类污染物的残留水平

根据南黄海苏北浅滩浒苔体内农药的检测结果（表 2）可知，ΣOCPs浓度范围为 15.86～56.93 ng·g-1，8种有机氯农药中α-六六六、六氯苯和δ-六六六均有检出，浓度高低顺序为α-六六六＞六氯苯＞δ-六六六；ΣPEs的浓度检出范围为 42.03～133.49 ng·g-1，4 种拟除虫菊酯类农药均有检出，其检出率均较高，其中联苯菊酯的浓度最高（64.05 ng·g-1）；ΣOPPs浓度范围为9.42～13.36 ng·g-1，9种有机磷中被检出的有乐果、甲基对硫磷、毒死蜱、杀螟硫磷、马拉硫磷、喹硫磷和三唑磷，其中甲基对硫磷，杀螟硫磷、马拉硫磷、喹硫磷和三唑磷的检出率较高，三唑磷浓度最高，达6.32 ng·g-1，其余每种 OPPs的浓度均小于 5 ng·g-1。由上述讨论可知，三类农药在浒苔体内的浓度以PEs的残留量为最高，OCPs次之，含量最少的是OPPs。这种差异可能与农药的环境行为、归趋密切相关，比如在不同环境条件中的半衰期、稳定性、光解和水解性质等以及农药的化学结构和理化性质，其中农药的正辛醇/水分配系数（lgKow）最为重要[24]，三类农药中，PEs的lgKow在 6～7 之间[25]，高于 OCPs的 lgKow4 左右[26-27]，而OPPs的lgKow值最小，基本在4以下[28]，由于分配系数与农药在浒苔体内的富集能力呈正相关[29]，因此浒苔体内三类农药的浓度顺序为PEs＞OCPs＞OPPs。

除了测定浒苔体内农药的残留浓度外，本文也将浒苔与其他文献中的海洋藻类进行了比较，位于福建省沿海地区的坛紫菜中ΣOCPs（α-HCH、β-HCH、γ-HCH、δ-HCH、p，p′-DDT、o，p′-DDT、p，p′-DDE、p，p′-DDD 和六氯苯）浓度为 6.59～23.49 ng·g-1，同位于福建省沿海地区的海带，其ΣOCPs浓度为3.68～10.32 ng·g-1[30]，大连的裙带菜体内ΣOCPs浓度为 1.37～2.81 ng·g-1[31]，江苏省连云港的条斑紫菜体内ΣOCPs浓度为 10.5～11 ng·g-1[32]，本文南黄海地区浒苔体内的Σ OCPs浓度稍高于其他地区海藻中的浓度，但仍属于同一浓度水平。福建省紫菜的ΣPEs（联苯菊酯、甲氰菊酯、高效氯氟氰菊酯和氯菊酯）浓度为13.1～24.7 ng·g-1，福建省海带的ΣPEs浓度为 12.75～24.18 ng·g-1[33]，本文的浒苔中ΣPEs浓度与其相似，其中联苯菊酯含量最高且检出率高达100%，在其他海藻中联苯菊酯的检出率也是最高的且浓度最大，是由于联苯菊酯的lgKow为6，易富集在生物体内，所以联苯菊酯的浓度较高。山东省青岛市的紫菜中ΣPEs含量较低，且只有三唑磷被检出，其浓度是 44～184 ng·g-1[34]，而本文浒苔的∑OPPs浓度范围为未检出至20.87 ng·g-1，比紫菜的浓度低，且都是三唑磷浓度最高，目前三唑磷已成为长江流域防治水稻螟虫的主打药品[35]，故江苏沿海地区使用量较大，导致三唑磷残留浓度较高。

表2 浒苔体内3类农药的组成及浓度（ng·g-1dw）Table2 Composition and concentrations of 3 types of pesticides in Enteromorpha（ng·g-1dw）

2.2 南黄海苏北浅滩浒苔中农药的富集效应

为了进一步了解三类农药在浒苔体内的富集效应，用生物-沉积物富集因子（BSAF=浒苔体内某农药类污染物的浓度/沉积物中某农药类污染物的浓度）[36]（图2）来描述浒苔对沉积物中农药的累积程度。沉积物中，∑OCPs、∑PEs和∑OPPs浓度范围分别为191.83～314.04、1.10～1.84 ng·g-1和5.74～7.63 ng·g-1（表3）。根据BSAF的计算结果可知，OCPs的BSAF值均小于1，OPPs的BSAF值为1～10之间，而PEs的 BSAF值最大，在10～100范围内。由于浒苔体内∑PEs浓度最大，而沉积物中∑PEs浓度又最小，所以PEs的BSAF值最大，无论是浒苔还是沉积物中的∑OPPs浓度都是三类农药中最小的且数值相近，因此，OPPs的BSAF值接近于1，而∑OCPs在沉积物中含量最大，且高于∑OCPs在浒苔体内的含量，故OCPs的BSAF值均小于1。由三类农药的BSAF值大小可知，其与沉积物中农药的含量相关，而沉积物中农药的含量与农药在沉积物有机碳与上覆水之间的分配系数Koc有关，这和浒苔中农药含量与lgKow有关相似。其中，OCPs的 lgKoc值均大于 3[29]，OPPs的lgKoc值大多小于3[37]，故沉积物中 OCPs的浓度大于OPPs，PEs例外，由于PEs在水中能快速降解[38]，因此，沉积物中PEs含量最低。

图2 浒苔-沉积物富集因子（BSAF）Figure2 The BSAF of Enteromorpha

2.3 南黄海苏北浅滩浒苔中农药类污染物的食用风险评价

本文所研究的浒苔体内有3种OCPs、4种PEs及7种OPPs农药被检出。其中，六氯苯的ADI值于20世纪70～80年代撤销，无法进行健康风险评价，故进行相关健康风险评价的农药有2种OCPs、4种PEs和7种OPPs（表4），maxTHQ值为每种农药在所有样点中的最大目标危险系数，顺序如下：α-六六六＞喹硫磷＞三唑磷＞联苯菊酯＞甲氰菊酯＞乐果＞甲基对硫磷＞高效氯氟氰菊酯＞氯菊酯＞δ-六六六＞杀螟硫磷＞毒死蜱＞马拉硫磷，所测的maxTHQ值均小于相应的ADI值，并且THQ值均远小于1，证明浒苔中每种农药单体对人体健康都没有潜在风险。aveTHQ值为每种农药在22个样点中的算数平均目标危险系数，α-六六六在OCPs中的贡献率为99.54%，联苯菊酯在PEs中的贡献率为80.51%，喹硫磷在OPPs中的贡献率为58.74%，表明从单一农药的平均含量看，α-六六六在OCPs中贡献最大，联苯菊酯在PEs中贡献最大，喹硫磷在OPPs中贡献最大。同时，每种农药的aveTHQ值都远小于1，对人体健康没有危害；HI值的分析结果显示：OPPs＞PEs＞OCPs，三类农药中 OPPs贡献最大。

表3 沉积物体内3类农药的组成及浓度（ng·g-1dw）Table3 Composition and concentrations of 3 types of pesticides in sediment（ng·g-1dw）

表4 食用浒苔3类农药摄入量与相关的健康风险评价Table4 Estimated daily intake and health risk of 3 types of pesticides in Enteromorpha

浒苔中三类农药于样点S1～S22的复合污染HI值如图3所示，其中，由于浒苔样品在3至5月为爆发中期，故以该时段所采集的全部浒苔样品为研究对象，可发现OCPs、PEs和OPPs HI范围分别为8.74×10-10～6.58×10-9、1.03×10-9～4.10×10-9和 0～8.48×10-9，三类农药的总 HI值为 2.62×10-9～1.59×10-8，远小于1，说明南黄海苏北浅滩浒苔体内被检测的三类农药浓度远低于能够危害人体健康的浓度水平。三类农药中，OPPs的贡献最大，在9种OPPs中检测出7种，是三类农药中被检测出数量最多的，并且在这7种农药中，喹硫磷的 ADI值是 5×102ng·kg-1·d-1，是三类农药中最小的，同时其ADI值远小于其他农药。因此，喹硫磷的THQ最大，故喹硫磷是OPPs中贡献最大的，所以，OPPs的HI值最高，与其他农药相比，OPPs对人类的健康风险是最大的。尽管PEs在风险监控中略低于OPPs，但是，PEs在浒苔体内的残留浓度和富集程度均是三类农药中最高的，其HI值较低，原因在于PEs检测的种类少于OPPs，故PEs的复合污染HI值低于OPPs，因此，仍需加强对PEs的控制。

根据《食品中农药残留风险评估指南》评估结果，从长期（慢性）摄入量可看出，与上述THQ值所得结论一致，从短期（急性）摄入量可看出，OPPs中三唑磷的健康风险最大，其后依次是乐果、甲基对硫磷、杀螟硫磷、毒死蜱和马拉硫磷，PEs中联苯菊酯的风险最大，高效氯氟氰菊酯的风险最小。短期评估与长期评估结果一致。三唑磷（ARfD=0.001 mg·kg-1bw）、乐果（ARfD=0.02 mg·kg-1bw）、联苯菊酯（ARfD=0.01 mg·kg-1bw）和高效氯氟氰菊酯（ARfD=0.02 mg·kg-1bw）急性毒性大，ARfD值小，则风险商值大，但风险程度最大的三唑磷残留风险商为3.64%，且远小于1，说明其健康风险相对较低。

图3 浒苔中复合农药污染的健康风险指数（HI）Figure3 Health hazard index（HI）of pesticides mixture in Enteromorpha

3 结论

（1）通过对浒苔体内三类农药含量的基础调查研究，可知浒苔体内农药含量为PEs＞OCPs＞OPPs，且PEs农药含量远高于其他两类农药。浒苔与其他可食用海藻相比，浒苔中OCPs和PEs的含量比其他地区的藻类高，而浒苔中OPPs的残留量略低于其他海藻。

（2）通过对目标物的生物-沉积物富集因子（BSAF）分析发现，浒苔对PEs的富集能力最强，OPPs的富集能力次之，而对OCPs的富集能力最弱，农药的富集能力可能与农药的环境行为、归趋以及农药的化学结构和理化性质有关，主要因素是农药的Koc和降解程度。

（3）对浒苔食用风险评价的结果表明，从单一农药的平均含量看，α-六六六、联苯菊酯和喹硫磷在三类农药中贡献最大，由于农药ADI值小，则风险商值大，但每种农药的aveTHQ值都远小于1，对人体健康没有危害；从复合污染HI值分析，三类农药的总HI值为 2.62×10-9～1.59×10-8，且远小于 1，说明南黄海苏北浅滩的浒苔从农药角度分析食用是安全的。

（4）根据《食品中农药残留风险评估指南》评估结果显示，从长期摄入量可得出与食用风险评价一致的结果，从短期摄入量可得出，有机磷农药中三唑磷的健康风险最大，拟除虫菊酯类农药中联苯菊酯的风险最大。由于农药ARfD值小，则风险商值大，但风险程度最大的三唑磷残留风险商为3.64%，且远小于1，说明南黄海苏北浅滩地区浒苔农药中健康风险相对较低。

[1]于仁成,刘东艳.我国近海藻华灾害现状、演变趋势与应对策略[J].中国科学院院刊,2016,31（10）：1167-1174.

YU Ren-cheng,LI Dong-yan.Harmful algal blooms in the coastal waters of China：Current situation,long-term changes and prevention strategies[J].Bulletin of Chinese Academy of Sciences,2016,31（10）：1167-1174.

[2]张林慧,张建恒,赵 升,等.2014年青岛海域消亡漂浮浒苔生理特征研究[J].上海海洋大学学报,2016,25（4）：591-598.

ZHANG Lin-hui,ZHANG Jian-heng,ZHAO Sheng,et al.Physiological characteristics of the floating ulva macroaglae along Qingdao coast in 2014[J].Journal of Shanghai Ocean Unversity,2016,25（4）：591-598.

[3]何 清,胡晓波,周峙苗,等.东海绿藻缘管浒苔营养成分分析及评价[J].海洋科学,2006,30（1）：34-38.

HE Qing,HU Xiao-bo,ZHOU Zhi-miao,et al.Evaluation on nutrition components of Enteromorpha linza[J].Marine Sciences,2006,30（1）：34-38.

[4]王文娟,赵 宏,米 锴,等.大型绿藻浒苔属植物研究进展[J].湖南农业科学,2009（8）：1-4.

WANG Wen-juan,ZHAO Hong,MI Kai,et al.Research advances on Enteromorpha：Large green algae[J].Hunan Agricultural Sciences,2009（8）：1-4.

[5]林英庭,朱风华,徐 坤,等.青岛海域浒苔营养成分分析与评价[J].饲料工业,2009,30（3）：46-49.

LIN Ying-ting,ZHU Feng-hua,XU Kun,et al.The nutrition analysis and evaluation on Enteromorpha prolifera at Qingdao sea area[J].Feed Industry Magazine,2009,30（3）：46-49.

[6]赵 明,陈建美,蔡 葵,等.浒苔堆肥化处理及对大白菜产量和品质的影响[J].中国土壤与肥料,2010（2）：66-70.ZHAO Ming,CHEN Jian-mei,CAI Kui,et al.Effects of Enteromorpha prolifera composting on yield and quality of Chinese cabbage[J].Soil and Fertilizer Sciences in China,2010（2）：66-70.

[7]孙 昊,王万章,王秀仑.利用孔石莼类海藻制作生物质板材的工艺[J].农业工程学报,2009,25（4）：309-312.

SUN Hao,WANG Wan-zhang,WANG Xiu-lun.Technical process for producing the biomass boards using sea alga Ulva pertusa kjellman species[J].Transactions of the Chinese of Agricultural Engineering,2009,25（4）：309-312.

[8]佚 名.海洋藻类或可成为新型能源[J].江西煤炭科技,2008（3）：4.

YI Ming.Marine algae or can be a new energy source[J].Jiangxi Coal Science&Technology,2008（3）：4.

[9]于 波,汤国民,刘少青.浒苔绿潮的发生、危害及防治对策[J].山东农业科学,2012,44（3）：102-104.

YU Bo,TANG Guo-min,LIU Shao-qing.Occurrence,harm and prevention of Entermorpha prolifera tidal[J].Shandong Agricultural Science,2012,44（3）：102-104.

[10]黄群腾.水环境中36种农药残留的同时分析方法及其应用[D].厦门：厦门大学,2008.

HUANG Qun-teng.Simultaneously determination method for 36 pesticides in aquatic environment and its application[D].Xiamen：Xiamen University,2008.

[11]王教凯.生物标志物指示南黄海沉积物中有机氯农药的分布及来源的可行性研究[D].青岛：中国海洋大学,2014.

WANG Jiao-kai.Feasibility study on biomarkers indicating the distribution and source of organochloride pestcides in sediment from Southern Yellow Sea[D].Qingdao：Ocean University of China,2014.

[12]王 凌.海洋环境中典型有机磷污染物分析及其生态效应研究[D].青岛：中国海洋大学,2006.

WANG Ling.Study on analytical methods of typical organophophorus pesticide pollutents in marine environment and their ecologic effects[D].Qingdao：Ocean Unversity of China,2006.

[13]李 星,曹彦忠,张进杰,等.ASE-SPE/GC-MS/MS同时测定海洋沉积物中71种农药残留[J].分析测试学报,2013,32（10）：1180-1186.

LI Xing,CAO Yan-zhong,ZHANG Jin-jie,et al.Simultaneous determination of 71 pesticides residues in marine sediments by ASE-SPE/GC-MS/MS[J].Journal of Instrumental Analysis,2013,32（10）：1180-1186.

[14]彭全材,宋金明,李 琛,等.胶州湾5种海藻中的多不饱和脂肪酸与有机氯农药共摄入风险的评估研究[J].海洋与湖沼,2014,45（1）：80-87.

PENG Quan-cai,SONG Jin-ming,LI Chen,et al.Risk assessment of ingesting polyunsa turated fatty acids and organochlorine pesticide residues in five seaweeds from Jiaozhou Bay[J].Oceanologia at Limnologia Sinica,2014,45（1）：80-87.

[15]蔡子豪,杜 晶,孙 彬,等.南黄海绿潮藻的分子鉴定及营养价值初探[J].浙江农业学报,2016,28（7）：1206-1215.

CAI Zi-hao,DU Jing,SUN Bin,et al.Molecular identification of green tide algae in the South Yellow Sea and preliminary study on its nutritivevalue[J].Acta Agricultures Zhejiangensis,2016,28（7）：1206-1215.

[16]陈美娜.基于脱附动力学评价养殖底泥中疏水性有机污染物生物有效性的研究[D].上海：上海海洋大学,2016.

CHEN Mei-na.Research on bioavailability of sediment-associated of aquaculture zone hydrophobic organic contaminants based on their desorption kinetics[D].Shanghai：Shanghai Ocean Unversity,2016.

[17]李振华.固相微萃取用于评价养殖底泥中菊酯类农药生物有效性和预测生物累积的研究[D].上海：上海海洋大学,2015.

LI Zhen-hua.Research on bioavailability of sediment-associated pyrethroids and prediction of biological accumulation based on solid phase microextraction[D].Shanghai：Shanghai Ocean Unversity,2015.

[18]徐佳艳,彭自然,和 庆,等.长三角地区池塘养殖水产品体内农药类污染与食用风险评价[J]生态毒理学报,2017,12（3）：485-495.

XU Jia-yan,PENG Zi-ran,HE Qing,et al.Assessment pesticide residues in pond aquaculture products in Yangtze River Delta and the health hazard for human consumption[J].Asian Journal of Ecotoxicology,2017,12（3）：485-495.

[19]高俊全,李筱薇,赵京玲.2000年中国总膳食研究：膳食铅、镉摄入量[J].卫生研究,2006,35（6）：750-754.

GAO Jun-quan,LI Xiao-wei,ZHAO Jing-ling.2000 Chinese total diet study：The dietary lead and cadmium intakes[J].Journal of Hygiene Research,2006,35（6）：750-754.

[20]Zhuang P,Mcbride M B,Xia H P,et al.Health risk from heavy metals via consumption of food crops in the vicinity of Dabaoshan mine,south China[J].Science of the Total Environment,2009,407（5）：1551-1561.

[21]World Health Organization.Principles and methods for the risk assessment of chemicals in food-Chapter 6：Dietary exposure assessment of chemicals in food[M].Geneva：World Health Organization,2010：112-115.

[22]张 磊,刘兆平.食品化学物风险评估中一些重要参数的选择和使用[J].中国食品卫生杂志,2015,27（3）：308-311.

ZHANG Lei,LIU Zhao-ping.The choices and usage of some important parameters in risk assessment of food chemicals[J].Chinese Journal of Food Hygiene,2015,27（3）：308-311.

[23]Authority E F S.Scientific opinion on risk assessment for a selected group of pesticides from the triazole group to test possible,methodologies to assess cumulative effects from exposure through food from these pesticides on human health[J].EFSA Journal,2009,7（9）：1167.

[24]Mailhot H.Prediction of algal bioaccumulation and uptake rate of nine organic compounds by ten physicochemical properties[J].Environmental Science&Technology,1987,21（10）：1009-1013.

[25]Feo M L,Ginebreda A,Eljarrat E,et al.Presence of pyrethroid pesticides in water and sediments of Ebro River Delta[J].Journal of Hydrology,2010,393（3）：156-162.

[26]毕峻奇.六六六、滴滴涕和多环芳烃在茶树体内迁移的模拟研究[D].福州：福建农林大学,2013.

BI Jun-qi.Simulation study on transfer of HCHs,DDTs and PAHs in tea tree[D].Fuzhou：Fujian Agriculture and Forestry University,2013.

[27]许国飞,张治国,孟 俊,等.我国大气中有机氯杀虫剂的研究进展[J].安徽农业科学,2009,37（6）：2690-2693.

XU Guo-fei,ZHANG Zhi-guo,MENG Jun,et al.Research progress on organochlorine pesticides in atmosphere in China[J].Journal of Anhui Agricultural Sciences,2009,37（6）：2690-2693.

[28]Benfenati E,Gini G,Piclin N,et al.Predicting logP of pesticides using different software[J].Chemosphere,2003,53（9）：1155-1164.

[29]沈国兴,严国安,彭金良,等.农药对藻类的生态毒理学研究Ⅱ：毒性机理及其富集和降解[J].环境科学进展,1999（6）：131-140.

SHEN Guo-xing,YAN Guo-an,PENG Jin-liang,et al.Study on ecotoxicology for pesticides to algaeⅡ：Toxic mechanism and accumulation,degradation[J].Advances in Environmental Science,1999（6）：131-140.

[30]叶江雷,弓振斌,金贵娥,等.福建沿海地区可食性海洋贝、藻类体中 HCH、DDT 的分布特征[J].福建水产,2012,34（6）：455-460.

YE Jiang-lei,GONG Zhen-bin,JIN Gui-e,et al.Distribution characteristics of HCH,DDT in the edible marine shellfish and algae of Fujian coast[J].Journal of Fujian Fisheries,2012,34（6）：455-460.

[31]吴金浩,李 爱,刘桂英,等.微波辅助萃取-气相色谱法检测裙带菜中的多氯联苯和有机氯农药残留[J].环境化学,2013,32（1）：156-161.

WU Jin-hao,LI Ai,LIU Gui-ying,et al.Determination of polychlorinated biphenyls and organochlorine pesticide residues in wakame by gas chromatography with electron capture detection using microwaveassisted extraction[J].Environmental Chemistry,2013,32（1）：156-161.

[32]周明莹,曲克明,陈碧鹃,等.紫菜中六六六、滴滴涕和扑草净的检测[J].渔业科学进展,2014,35（6）：128-132.

ZHOU Ming-ying,QU Ke-ming,CHEN Bi-juan,et al.Residual determination method of HCHs,DDTs and prometryn in porphyra[J].Progress in Fishery Sciences,2014,35（6）：128-132.

[33]谢湘云.海洋环境样品中痕量有机氯和拟除虫菊酯农药的测定方法研究[D].厦门：厦门大学,2007.

XIE Xiang-yun.Determination of trace organochlorine and pyrethroid pesticides in marine environmental samples[D].Xiamen：Xiamen University,2007.

[34]孙晓杰,郭萌萌,王苏玥,等.分散固相萃取-在线凝胶色谱-气相色谱-质谱联用法快速检测紫菜中的农药多残留[J].色谱,2014,32（10）：1124-1130.

SUN Xiao-jie,GUO Meng-meng,WANG Su-yue,et al.Rapid determination of pesticide multiresidues in porphyra by dispersive solidphase extraction coupled with online gel permeation chromatographygas chromato graphy mass spectrometry[J].Chinese Journal of Chromatography,2014,32（10）：1124-1130.

[35]杨英利,朱小燕,李爱民,等.3种淡水藻对三唑磷的降解研究[J].环境科学研究,2007,20（1）：85-89.

YANG Ying-li,ZHU Xiao-yan,LI Ai-min,et al.Study on biodegradation of triazophos by three species of freshwater algae[J].Research of Environmental Science,2007,20（1）：85-89.

[36]苏秋克,祁士华,吴辰熙,等.洪湖特色水产品对湖水及沉积物中有机氯农药的积累模式[J].地质科技情报,2007（4）：85-90.

SU Qiu-ke,QI Shi-hua,WU Chen-xi,et al.Organism productions accumulating model of organochlorine pesticides in water and sediments from Honghu Lake,Hubei Province,China[J].Geological Science and Technology Information,2007（4）：85-90.

[37]Daré J K,Silva C F,Freitas M P.Revealing chemophoric sites in organophosphorus insecticides through the MIA-QSPR modeling of soil sorption data[J].Ecotoxicology and Environmental Safety,2017,144：560-563.

[38]岳文洁.氯氰菊酯农药对海洋微藻和大型海藻的毒性效应研究[D].广州：暨南大学,2009.

YUE Wen-jie.Studied the toxic effects of cypermethrin on marine microalga and seaweed[D].Guangzhou：Jinan University,2009.