一种淡水水生动物急性水质基准的推导方法技术

本发明专利技术公开了一种淡水水生动物急性水质基准的推导方法，包括以下步骤：a、急性毒性数据的收集和筛选；步骤b、依次分别求出筛选出的急性毒性数据的属平均急性值及对应的累积概率；步骤c、以步骤b中求得的属平均急性值的自然对数为横坐标X，以其对应的累积概率P为纵坐标Y构造三次样条函数S(x)，对Y赋值，对应求得一个X值，对X取反对数，得短期危险浓度STHC100Y，即为针对目标污染物的急性水质基准值。本发明专利技术产生的有益效果是数据可信度高、计算过程简单、客观性强和实效精确性好，同时，克服了现有水质基准计算模型的不足，丰富了水质基准的计算方法，为科学地计算水质基准提供了技术支持。

Derivation of an acute water quality standard for freshwater aquatic animals

The present invention discloses a method for deriving an acute water quality reference for freshwater aquatic animals, including the following steps: A, collection and screening of acute toxicity data; step B, in order to determine the average acute value of the selected acute toxicity data and the corresponding cumulative probability respectively; step C, an average acute value obtained in step b. The natural logarithm of the natural logarithm is the horizontal coordinate X, and the corresponding cumulative probability P constructs the three cubic spline function S (x) for the longitudinal coordinates Y, assigns a value to the Y, corresponds to a X value, and gets the opposition of X, and obtains the short-term dangerous concentration STHC100Y, that is, the acute water quality reference value for the target pollutants. The beneficial effect of the invention is high reliability of data, simple calculation process, strong objectivity and good accuracy, at the same time, it overcomes the shortage of the existing water quality standard calculation model, enriches the calculation method of water quality datum, and provides technical support for the scientific calculation of water quality datum.

【技术实现步骤摘要】
一种淡水水生动物急性水质基准的推导方法
本专利技术属于水质基准计算
，具体涉及一种淡水水生动物急性水质基准的推导方法。
技术介绍
水环境质量基准(以下简称“水质基准”)，是指水环境中的污染物质或有害因素对人体健康、水生态系统与使用功能不产生有害效应的最大剂量或水平。淡水水生生物水质基准(freshwaterqualitycriteriafortheprotectionofaquaticorganisms)是指能够保护淡水水生生物及其生态功能的水质基准，包括短期水质基准和长期水质基准。水质基准的推导技术路线主要为：(1)收集评价数据，如动物急性毒性、动物慢性毒性、植物毒性、生物累积或其他数据，(2)根据毒性数据分别推导出水生动物的最终急性毒性值(FAV)、最终慢性毒性值(FCV)、最终残留毒性值(FRV)、最终植物毒性值(FPV)或最低重要效应值，(3)计算基准最大浓度CMC和基准持续浓度CCC用来表示水生生物水质基准值，(4)对结果和数据评估。目前淡水水生生物水质基准的推导方法主要有评价因子法、毒性百分数排序法及物种敏感度分布曲线法三种。评价因子法的优点是简单易行，所需数据量少且适用范围广，但由于评价因子法推导出的基准值仅取决于最敏感物种的毒性数据且评价因子的确定通常是从经验出发，因此结果的准确性存在较大的不确定因素。毒性百分数排序法计算基准时，仅仅考虑最敏感的四个属种，但当对生态系统具有重要生态价值或对人类具有经济、娱乐价值的属种比以上四个属种更敏感时，即有可能造成基准值偏高而造成“欠保护”。并且该方法计算较为复杂，如以下基准最大浓度CMC的计算公式为例即可说明其复杂性。此外，对于大多数污染物来说，现有的毒性数据无法满足毒性百分数排序法的数据要求，要使用该方法还需要开展大量的毒理学实验研究，而这是一项耗时耗资耗力的工作。物种敏感度分布曲线法利用不同数学模型对所有毒性数据进行拟合，根据相关系数选择较优的拟合方程，但是实际工作中不可能穷尽所有的计算公式，不排除存在更好的数学模型，因此，该方法受研究人员数理知识背景等因素影响大，具有一定的不确定性。近年来，我国不断完善各项环境标准，尤其在水环境领域更是密集，不少标准都在制定或修订中。然而，作为标准的制定依据，水质基准还缺乏植根于我国国情的持续研究，难以为标准制定提供有力支撑，导致现行水质标准的科学准确性和实效精确性都可能存在问题。CN103336889A公开了一种水生生物水质基准推导方法，采用毒性百分数排序法计算基准，受试生物采用水生动物和水生植物，采用了大量毒性数据，但最终仅选择其中4个属的数据进行分析，造成“欠保护”问题，同时，由于在目前急性毒性数据的收集中，水生动物的数据更容易收集，数据更多且筛选出的数据可信度更高；而大部分水生植物对污染物的灵敏度较低，所需毒性实验的时间更久，并且数据较少，进行对比后，难将异值点排除，容易对基准的推导造成偏差，造成水质基准的精确性较差。CN107392780A公开了一种用于淡水水生生物水质基准预测的系统，采用物种敏感度分布曲线法利用不同的数学模型对所有毒性数据进行拟合。然而，当毒性数据少于5个时则不能采用本系统；此外，当来自相同分类学类群的多个生物均为受试生物时，物种的相似性可能导致存在多个相近的种平均急性值，进而使多个相邻的数据点在分布曲线中近乎垂直，这种分布情况可能会很难找到一个合适的拟合度满足要求的分布模型。样条插值是一种工业设计中常用的、得到平滑曲线的一种插值方法，三次样条又是其中用的较为广泛的一种。三次样条函数插值是通过构造多项式来形成一条把所有主干点连接起来的平滑曲线，然后进行插值的方法。三次样条插值既可保持分段低次插值多项式，又可提高插值函数光滑性，但至今仍无人将其用于水质基准的推导计算。
技术实现思路
针对以上技术问题，本专利技术提供了一种数据可信度高、计算过程简单、客观性强和实效精确性好的淡水水生动物急性水质基准的推导方法，同时，克服了现有水质基准计算模型的不足，丰富了水质基准的计算方法，为科学地计算水质基准提供了技术支持。为了解决以上技术问题，本专利技术提供的一种淡水水生动物急性水质基准的推动方法，包括以下步骤：步骤a、急性毒性数据的收集和筛选；步骤b、数据处理：依次分别求出筛选出的急性毒性数据的属平均急性值GMAV及对应的累积概率P；步骤c、急性基准值的计算：以步骤b中求得的属平均急性值GMAV的自然对数为横坐标X，以其对应的累积概率P为纵坐标Y，构造三次样条函数S(x)，对Y赋值，其中Y＝0～1，对应求得一个X值，对X取反对数，得短期危险浓度STHC100Y，即为针对目标污染物的急性水质基准值。优选地，所述步骤a中急性毒性数据的受试物种为淡水水生动物，包括本地物种和在当地大量繁殖的外来物种，符合三门八科的物种筛选原则。优选地，所述急性毒性数据的筛选遵循以下原则：毒性试验必须设置对照试验，对照组和试验组的理化条件、试验生物和试验步骤应完全相同，不能使用物种表现出疾病症状或者大量死亡的对照组试验结果；不能使用未有明确测试终点、暴露时间的试验数据；水蚤或其他枝角类动物所做实验应使用龄期小于24h的生物，其急性毒性值采用48hEC50或LC50表示；其他水生动物的急性毒性值采用96hEC50或LC50表示，若暴露时间不足96h，采用48hEC50或LC50代替。优选地，所述步骤b的具体步骤如下：(1)利用筛选出的急性毒性数据，计算出目标污染物对生物的属平均急性值GMAV；(2)将所述的属平均急性值GMAV按从小到大的顺序并分别赋予从1至N的排号R，其中，R为自然数，排号1为最小属平均急性毒性值，排号N为最大属平均急性毒性值；(3)利用公式P＝R/(N+1)，计算每个属的属平均急性值GMAV的累积概率P，其中，P代表水生生物被危害的比例。进一步地，所述步骤(1)中同一物种有多个有效毒性终点LC50或EC50，则以多个毒性终点LC50或EC50的几何平均值作为该种平均急性值SMAV，计算公式如下：其中：EC50n代表同一物种中，第n个有效毒性终点；n代表同一物种中有效毒性终点的个数；EC50可用LC50替代。进一步地，所述步骤(1)中同一个属的一个以上物种的急性毒性数据存在时，属平均急性值GMAV为种平均急性值SMAV的几何平均值，计算公式如下：其中：SMAVn代表第n种物种的种平均急性值；n代表物种个数。进一步地，所述步骤(2)中有多个属平均急性值GMAV是相等的，任意将它们排成连续的等级。优选地，所述步骤c的具体步骤如下：对所述的属平均急性值GMAV取自然对数后作为横坐标X，记为X＝(x1,x2,....,xN)，以其对应的累积概率P为纵坐标Y，记为Y＝(y1,y2,.....,yN)；在区间[x1,xN]上拟合求得三次样条函数S(x)，使它满足S(xn)＝yn，n＝1,2….N，同时满足边界条件S”(x1)＝S”(xN)＝0；基于风险安全评价方法对Y赋值，根据该Y值确定其对应的X的区间[xi,xi+1]，进而确定在区间[xi,xi+1]上的三次样条函数，将Y值代入该三次样条函数，对应求得一个X值，将X取反对数，得到目标污染物的短期危险浓度STHC100Y，其中，Y的取值范围为0～1。优选地，所述本文档来自技高网...

【技术保护点】
一种淡水水生动物急性水质基准的推导方法，其特征在于，包括以下步骤：步骤a、急性毒性数据的收集和筛选：所述的急性毒性数据的受试物种为淡水水生动物，包括本地物种和在当地大量繁殖的外来物种，符合三门八科的物种筛选原则；步骤b、数据处理：依次分别求出筛选出的急性毒性数据的属平均急性值GMAV及对应的累积概率P；步骤c、急性基准值的计算：以步骤b中求得的属平均急性值GMAV的自然对数为横坐标X，以其对应的累积概率P为纵坐标Y，构造三次样条函数S(x)，对Y赋值，其中Y＝0～1，对应求得一个X值，对X取反对数，得短期危险浓度STHC100Y，即为针对目标污染物的急性水质基准值。

【技术特征摘要】
1.一种淡水水生动物急性水质基准的推导方法，其特征在于，包括以下步骤：步骤a、急性毒性数据的收集和筛选：所述的急性毒性数据的受试物种为淡水水生动物，包括本地物种和在当地大量繁殖的外来物种，符合三门八科的物种筛选原则；步骤b、数据处理：依次分别求出筛选出的急性毒性数据的属平均急性值GMAV及对应的累积概率P；步骤c、急性基准值的计算：以步骤b中求得的属平均急性值GMAV的自然对数为横坐标X，以其对应的累积概率P为纵坐标Y，构造三次样条函数S(x)，对Y赋值，其中Y＝0～1，对应求得一个X值，对X取反对数，得短期危险浓度STHC100Y，即为针对目标污染物的急性水质基准值。2.根据权利要求1所述的一种淡水水生动物急性水质基准的推导方法，其特征在于，毒性试验必须设置对照试验，对照组和试验组的理化条件、试验生物和试验步骤应完全相同，不能使用物种表现出疾病症状或者大量死亡的对照组试验结果；不能使用未有明确测试终点、暴露时间的试验数据；水蚤或其他枝角类动物所做实验应使用龄期小于24h的生物，其急性毒性值采用48hEC50或LC50表示；其他水生动物的急性毒性值采用96hEC50或LC50表示，若暴露时间不足96h，采用48hEC50或LC50代替。3.根据权利要求1所述的一种淡水水生动物急性水质基准的推导方法，其特征在于，所述的步骤b的具体步骤如下：(1)利用筛选出的急性毒性数据，计算出目标污染物对生物的属平均急性值GMAV；(2)将所述的属平均急性值GMAV按从小到大的顺序并分别赋予从1至N的排号R，其中，R为自然数，排号1为最小属平均急性毒性值，排号N为最大属平均急性毒性值；(3)利用公式P＝R/(N+1)，计算每个属的属平均急性值GMAV的累积概率P，其中，P代表水生生物被危害的比例。4.根据权利要求3所述的一种淡水水生动物急性水质基准的推导方法，其特征在于，所述步骤(1)中同一物种有多个有效毒性终点LC50或EC50，则以多个毒性终点LC50或EC50的几何平均值作为该种平均急性值SMAV，计算公式如下：

【专利技术属性】
技术研发人员：陈曲，宋凡浩，白英臣，吴丰昌，赵天慧，栗婷婷，王珺瑜，刘大庆，陈季康，
申请(专利权)人：中国环境科学研究院，
类型：发明
国别省市：北京,11