TCSES 141-2024 页岩气开发地块特征污染物土壤生态阈值确定技术指南

ICS13.020.10CCSZ04团 体 标 准T/CSES141-2024页岩气开发地块特征污染物土壤生态阈值确定技术指南Technicalguidanceforderivingsoilecologicalthresholdsforcharacteristiccontaminantsinshalegasdevelopmentsites（发布稿）2024-04-16发布 2024-04-16实施中国环境科学学会发布T/CSES141T/CSES141—2024目 次前 言 II范围 1规范性引用文件 1术语和定义 1土壤生态阈值确定程序 3土壤环境暴露场景分析 4生态受体及生态过程的筛选 4特征污染物生态毒性数据获取及筛选 5土壤生态阈值的推导 6土壤生态阈值推导的不确定性分析 8土壤生态阈值的审核 8附录A（资料性）页岩气开发地块特征污染物土壤环境暴露场景分析示例 10附录B（资料性）常用的毒性数据库 11附录C（资料性）中国西南地区页岩气开发区域土壤性质参数 12附录D（资料性）物种敏感性分布曲线拟合函数公式 16附录E（资料性）拟合优度评价模型 17附录F（资料性）页岩气开发地块特征污染物土壤生态阈值确定报告编制大纲 18参考文献 19I前 言GB/T1请注意本文件的某些内容可能涉及专利。本文件的发布机构不承担识别专利的责任。本文件由生态环境部华南环境科学研究所和中国热带农业科学院环境与植物保护研究所提出。本文件由中国环境科学学会归口。本文件起草单位：生态环境部华南环境科学研究所，中国热带农业科学院环境与植物保护研究所，中国矿业大学（北京本文件主要起草人：林必桂，任明忠，王钰钰，陈启宇，崔凡，赫云兰，周上博，张六一，刘志磊。II范围本文件规定了页岩气开发地块土壤生态阈值确定的程序、方法与技术要求。本文件适用于我国页岩气开发地块及其周边区域自然保护地及农用地土壤中特征污染物的生态阈值制定。本文件不适用于致病性生物、放射性物质以及其他物理性污染（如噪声、热污染）的土壤生态阈值的确定。规范性引用文件下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。其中，注日期的引用（包括所有的修改单适用于本文件。GB/T41611页岩气术语和定义HJ25.1建设用地土壤污染状况调查技术导则HJ/T166土壤环境监测技术规范HJ682建设用地土壤污染风险管控和修复术语HJ831淡水生物水质基准推导技术指南HJ1185区域性土壤环境背景含量统计技术导则（试行）HJ1231土壤环境词汇NY/T395农田土壤环境质量监测技术规范术语和定义下列术语和定义适用于本文件。页岩气shalegas以游离态、吸附态为主，少量溶解态，赋存于富有机质页岩层段中的天然气。注：一般具有自生自储、大面积连续分布、储层低孔低渗、单井无自然产能或低产，需通过增产改造才能获得工业气流等特点。[来源：GB/T41611—2022,3.1]页岩气开发地块shalegasdevelopmentsites页岩气开发矿权区域，包括页岩气生产平台在内的矿权所属地块。注3~65气开发地块属于第二种情况。地块生态风险评估ecologicalriskassessmentforlandparcels1对地块土壤中的污染物危害动物、植物、微生物和其他生态系统过程与功能的概率或水平与程度进行评估的过程。[来源修改于：HJ682-2019,2.4.4]土壤生态阈值soilecologicalthreshold,SET保证页岩气开发场地及周边生态受体安全的土壤中污染物的最大允许含量。注：本文件土壤生态阈值包含外源生态安全阈值和土壤环境背景含量两部分。土壤环境背景含量environmentalbackgroundconcentrationofsoil,Cb一定时间条件下，仅受地球化学过程和非点源输入影响的土壤中元素或化合物的含量。[来源：HJ1185—2021,3.1]外源土壤生态阈值exoticecologicalsafetythreshold定环境条件和生态系统类型，推导出的土壤中外来化学物质浓度或暴露水平的定量标准。生态受体ecologicalreceptor一般指地块及其周边环境中可能受到污染物影响的人群或生物类群，也可泛指地块周边受影响的功能水体和自然及人文景观等。[来源：HJ1231—2022,5.7]敏感受体sensitivereceptor受地块污染物影响的潜在生物类群中，对污染物反应最敏感的受体。[来源：HJ1231—2022,5.7]ecologicalprocess是指生物群落与非生物环境之间，以及生物群落内各种生物之间相互作用的过程，包括基于微生物的土壤生态过程，如硝化作用、有机质矿化、酶活性等。3.10无观察效应浓度noobservedeffectconcentration,NOEC外源污染物在规定暴露条件下，与对照组相比，对受试生物未产生显著效应（P＞0.05）的最高浓度。[来源：HJ831—2022,3.11,有修改]3.11风险商riskquotient,RQ通常用于对某个单一化合物进行危害风险评估，其计算方式是通过实际检测或者利用模型预测出的环境中该化合物的浓度与表明此物质胁迫程度的毒理数值相比，得到风险商值（R2[来源：T/CSES23—2021,3.20,有修改]3.12物种敏感度分布speciessensitivitydistribution,SSD描述不同物种由于生活史、生理构造、行为特征和地理分布等的不同，对某一污染物的敏感性差异遵循的概率分布规律。[来源：HJ831—2022,3.7]3.13不确定性分析uncertaintyanalysis对风险评估过程的不确定性因素进行综合分析评价。[来源：HJ1231—2022,5.21,有修改]土壤生态阈值确定程序页岩气开发地块特征污染物土壤生态阈值的确定主要包括以下步骤（1分析页岩气开发地块特征污染物的土壤环境暴露场景；确定需要保护的生态受体和生态过程及其暴露途径；有效毒性数据的获取及筛选以及土壤污染物背景含量获取；阈值推导方法的比较选择及确定；阈值的推导；阈值的不确定性分析；阈值的审核。3数据归一化数据归一化页岩气开发地块特征污染物暴露场景分析页岩气开发地块特征污染物暴露场景分析SSD法外推排序法外推评估因子法外推满足SSD法数据要求满足SSD法数据要求否是满足排序法数据要求是否满足评估因子法数据要求是否外源土壤生态阈值(SETSoi_Ex)外推土壤污染物背景值(Cb)有效毒性数据的获取及筛选文献数据生态毒性数据库确定需要保护的生态受体和生态过程土壤生态阈值(SETSoil)的确定及不确定度分析生态毒性实验收集生态毒性数据确定页岩气开发地块特征污染物土壤生态阈值审核通过确定页岩气开发地块特征污染物土壤生态阈值图1页岩气开发地块特征污染物土壤生态阈值确定程序土壤环境暴露场景分析页岩气开发地块特征污染物生态安全阈值的推导首先应开展土壤环境暴露场景构建。土壤环境暴露场景构建主要分析识别以下几个方面（：与土壤污染相关的环境介质及其在页岩气开发过程中的产生环节；土壤高关注污染物的识别及特征污染物的筛选；生态受体及生态过程的筛选以保护页岩气开发地块土壤环境及周边自然保护地和农田生态系统为主，其关键受体包括陆生（生产者（消费者（分解者态阈值需包含但不限于如下类型的重要土壤生态受体和生态功能：43个；优先选择植物生长指标，包括生物量、出苗率、芽长、根伸长和产量等指标。2长、生物量和生长率等等指标。3个。特征污染物生态毒性数据获取及筛选毒性数据来源物种毒性参数的数据可通过土壤生态毒性实验以及各种毒性数据库和公开发表的文献等来源获取。常用的毒性数据库参考附录B。（8章应的生物毒理学实验以补充相关毒性数据；补充试验应充分考虑本土关键物种为受试对象。毒性数据的筛选原则从毒性数据库或者文献资料中获取的毒性数据的筛选应遵循可靠性、丰富性、实用性的原则：（GB/T21809GB/T32O1269、O1698、O1863应能根据文献资料确定土壤污染物的暴露时间和毒性终点（如生长、繁殖，并可根据剂量-ECx；pH等；体系的毒性数据；陆生植物、无脊椎动物和土壤生态过程相关毒性效应数据应分别筛选归类，单独评估。毒性数据的归一化处理通用要求SSD数据归一化到统一的土壤条件下（我国西南地区页岩气开发区域土壤理化性质参数见附录C16.0me/100g土壤（ECx）EC10。自然保护地土壤特征污染物数据归一化自然保护地土壤特征污染物数据归一化采用种间外推归一化法，即以生物分类学相近的物种建立的回归模型为基础，数据按下列公式进行归一化：�𝑎𝑎∙�𝑝𝑝𝑝𝑝

−𝑝𝑝𝑝𝑝+𝑏𝑏∙𝑙𝑙𝑙𝑙𝑙�𝑂𝑂𝑂𝑠𝑠𝑠�+𝑐𝑐∙𝑙𝑙𝑙𝑙𝑙�𝑂𝐶𝑂𝑠𝑠𝑠��𝑥𝑥𝑠𝑠𝑠

=𝑥𝑝∙

𝑠𝑠𝑠𝑠𝑠𝑠

𝑝𝑝

𝑂𝑂𝑂𝑂𝑝𝑝

𝑂𝑂𝐶𝐶𝑂𝑂𝑝𝑝

（式1）式中：ECx_std——标准土壤条件下的x%效应浓度；5pHstd——标准土壤条件下的土壤pH值；OCstd——标准土壤条件下的土壤有机质含量；CECstd——标准土壤条件下的土壤阳离子交换量；ECx_p——实验土壤条件下的x%效应浓度；pHp——实验土壤条件下的土壤pH；OCp——实验土壤条件下的土壤有机质含量；CECp——实验土壤条件下的土壤阳离子交换量；a、b、c——相近物种的log型归一化模型中pH、OC、CEC的回归模型参数。注：若无法从数据库和文献资料中获得足够数量的实验土壤条件下ECx_std毒性数据(如EC10)，则可利用具有明显剂量效应关系的原始数据（如EC50）通过合适方法推导出ECx_std。农用地土壤特征污染物数据归一化1质相对应的物种毒性数值且有足够数据建立自身的回归模型时，也可利用多元线性归一化法将毒性参数值归一化到统一的土壤条件下，如公式（2）所示：lo(𝑥𝑥)=𝑎𝑎∙lo(𝑝𝑝𝑝𝑝+𝑏𝑏∙lo(𝑂𝐸𝐸+𝑐𝑐∙lo(𝐸𝐸)𝑑𝑑 （式中，ECx——x%效应浓度；pH——供试土壤pH值；OC——供试土壤有机质含量；CEC——供试土壤阳离子交换量；a、b、c、d——回归模型参数。土壤特征污染物数据归一化取值原则对采用同一物种和同一供试土壤获得的同一测试终点的毒性数据，取几何平均值。（（1（2通过土壤性质对各污染物的生态毒性数据进行（R2）p值大（E（12CCEC土壤生态阈值的推导推导方法本标准规定的页岩气开发地块特征污染物土壤生态阈值的推导采用外源生态阈值结合土壤环境背景含量方法，如式（3）所示：SETSoil=SETSoil_ExCb （式中，SETSoil——页岩气开发地块特征污染物土壤生态阈值；SETSoi_Ex——推导的外源土壤生态阈值；Cb——土壤环境背景含量。6（8.28.4节SES_xSSDSSDCb（HJ要求执行0。95%SSDHC55%分布预测值作为危害浓度。SSD法的外源土壤生态阈值推导SSDSSD法的外源生态安全阈值推导的毒性数据需要满足以下要求：2个营养级及以上，如生产者、消费者、分解者。8个以上（10-15个及以上为优。SSDESTSoi_Ex正态分布、Bu、eibll、aa等分布函数（D。拟合分布模（rootmeansquareerror，RME）确定最佳的拟合函数（。基于排序分布法的外源土壤生态阈值推导SSDECx（有不同响应水平时，优先选择或转化为EC10）分别按照从小到大的顺序进行排序，得到毒性数据序列，按公式（4）计算毒性数据在序列中的百分位数。𝑗𝑗=𝑖𝑖𝑛𝑛+1

⸱100 （4）式中：j——毒性数据序列中某一数据的百分位数；i——毒性数据序列中某一数据的顺序数；n——毒性数据总数。（j）ECx（EC10）SSD5%ESTSoil_Ex值。基于评估因子法的外源土壤生态阈值推导若无法获得足够有关土壤环境中生产者、消费者或分解者的生态毒理学数据，可选择毒性数据的最低值，采用评估因子法进行外推，评估因子值（AF）选择参考表1。𝐸𝑠𝑙𝑙𝑖𝑖𝑙𝑙𝐸𝐸𝑥𝑥=𝑥𝑥/𝐴𝐴𝐴𝐴 5）表1 评估因子AF取值数据类型评估因子生物（如植物、蚯蚓或微生物）短期毒性试验的L（E）C501000一种生物（如植物）长期毒性试验的NOEC值100来自3个营养级的生物长期毒性试验的NOEC值507来自3个营养级的3种以上生物长期毒性试验的NOEC值10已知物种敏感性分布曲线（SSD法）1~5，根据具体情况调整野外试验数据或生态模型试验数据根据具体情况决定取舍土壤生态阈值推导的不确定性分析不确定性分析要求（95%置信区间环境背景含量（Cb）的不确定度范围。推导的外源土壤生态阈值（SETSoil_Ex）的不确定性分析①SSD法推导阈值（HC5）0.05显著水平上的、可拟合度低95%95%SSD95%信区间的数值。土壤环境背景含量（Cb）的不确定性分析（如部分重金属元素（（𝑥𝑥算术标准差（、几何平均值（M、几何标准差（）J1852217..2条款9%（正态分布采用𝑥𝑥2𝐸𝑀/𝐷2~𝑀𝐷2，2.5%~97.5%分位数。土壤生态阈值（SETSoil）的不确定性分析Cb0时，SETSoilSETSoil_ExCb0时，SETSoilSETSoil_ExCb两个输入参数不确定度的平方和的平方根。土壤生态阈值的审核自审核项目土壤生态阈值的最终确定需要仔细审核阈值推导所用数据以及推导步骤，以确保阈值合理可靠。自审项目如下：67节规定的数据质量要求。GB/T3359、GB/T4883和GB/T8056等标准的规定，分析并判别是否存在明显异常数据。是否遗漏其它重要数据。专家审核项目土壤生态阈值推导数据是否可靠。8土壤生态阈值推导过程是否符合技术标准。土壤生态阈值的确定是否合理。土壤生态阈值的不确定性分析是否合理。9T/CSES141T/CSES141—2024附 录 A（资料性）页岩气开发地块特征污染物土壤环境暴露场景分析示例页岩气开发各工程阶段（井场建设、钻井、压裂试气、采气、集输等）产生的污染介质及其所含污染物类型非常复杂，因此应从污染介质特征和化学物质生态毒性效应等方面筛选优先关注的特征污染物，并进一步分析其土壤环境暴露场景。洗井废水、压裂返排液、废水处理站的尾水以及凝析油等；固体废弃物主要包括钻井岩屑、水基岩屑、油基岩屑、清管废渣和污水处理站污泥等。（一开至三开钻井（洗井（射孔压裂（冒滴漏、池底底部出现裂缝）作业事故引发的溢流井喷等事故工况则会给土壤环境带来很大的生态风险。页岩气开发使用的钻井液、压裂返排液等添加种类繁多的化学药剂（（也会被钻井液或压裂返排液带出地表。目前，受关注的污染物种类主要有以下几种：季铵盐类化合物、非离子表面活性剂、芳烃化合物、重金属、高含量无机盐。结合污染物生态毒理效应研究结果，优先关注的特征污染物包括多环芳烃、重金属（Ba,,b,Tl,d土壤环境污染物的关键受体根据暴露途径分成土壤直接接触、土壤和食物摄入、土壤地3类。土壤直接接触途径是土壤生态阈值确定的主要暴露途径，并以陆生植物（生产者无脊椎动物（消费者）和土壤微生物（分解者）及其主导的生态功能作为该途径的关键受体。（尤其是采用本土关键物种时10附 录 B（资料性）常用的毒性数据库常用的毒性数据库如下所示：欧洲化学品管理局（ECA）的国际统一化学品信息数据库（CLD：oeECHA(europa.eu)，https://iuclid6.echa.europa.eu；（RMtp:/.e-；（NLM）TOXLINE毒理学数据库：/newtoxnet/toxline.htm；ECOTOX毒性数据库：ECOTOX|Search()，/ecotox；Elsevier公司的“ECOTOΧ-CD”毒性数据库：；（6）中国科学院南京土壤研究所建立的场地土壤污染物环境与毒性基础数据库：:8995/#/login等。11附 录 C（资料性）中国西南地区页岩气开发区域土壤性质参数（htt:/d3oilcb.n/xtn/j/ntduton（不同土纲或土种等C1（-长宁地区及周边区域pH6.502.0%16.0me/100g土壤。表C1中国不同典型类型土壤有机质含量状况区域土种土壤pH有机质含量有效阳离子交换量（/）（/）（%）（me/100g土壤）（me/100g土壤）四川长宁-昭通区域暗红泥田6.136.0—6.52.89%10暗紫黄泥土7.006.6—7.420.518—23暗紫鸡粪土田6.535.9—7.42.72%18暗紫砂泥土6.506.6—7.41.63%20.518—23暗紫石骨土6.906.6—7.41.55%18.516—21白鳝泥田5.134.6—6.62.25%白眼砂土7.957.7—8.51.16%6北碚黄泥土4.885.0—6.45.76%1210—15扁砂黄泥土4.804.4—6.04.96%10大足黄泥土4.804.0-6.3/4.54.77%10浮泥7.406.9—7.66.60%39.536—43钙紫石骨土8.157.6—8.61.06%16海漂土5.305.5—6.012.02%汉昌潮砂泥田7.006.2—7.72.33%黑汤大土6.005.8—6.317.55%2928—30红泥5.504.0—5.8/5.52.23%10红紫砂泥土5.455.0—6.22.56％6红紫砂土4.704.5—6.31.30%红棕紫泥土8.407.8—8.61.28%23.520—27红棕紫砂泥土8.307.8— 1.02%22.520—25厚褐红土6.576.4—6.82.00%16厚毡洼炭土7.036.8—7.219.40%华坪灰泡砂土5.855.5—6.27.65%黄大土5.756.0左右4.84%12黄红紫黄泥土7.877.7—8.31.18%19黄泡土7.877.4—8.05.50%2015—25黄羊肝土5.835.3—6.42.45%1212黄紫泥5.375.2—5.91.91%10黄紫泥田6.686.0— 1.83%1813—23灰红泥田6.836.3—7.54.82%12灰泡大泥5.005.0左右5.60%19.518—21灰泡大土5.505.5—6.05.63%22灰紫砂泥田7.207.0左右3.27%灰棕潮泥土8.077.5—8.01.78%20灰棕潮砂泥土7.937.1—8.41.07%12灰棕黄紫泥土6.356.7—7.31.23%2218—26灰棕紫砂泥土6.556.6—7.41.16%20.519—22夹白鳝黄泥田6.605.0—7.32.32%15夹黄紫泥田7.906.7—8.51.80%2521—29假白鳝紫泥田5.904.5—6.32.26%1513—17剑阁紫泥土7.406.6—7.51.53%18.516—21靖外红土6.005.3—6.43.04%砾质褐红大土6.906.9左右1.83%9.6面黄泥土5.801.47%6.56—7名山白鳝泥土4.905.0左右1.59%8.5内家紫砂泥田7.907.0—8.12.04%1712—22彭山灰棕砂土8.057.5—8.61.16%7蓬安死黄泥田5.905.7—7.02.22%1513—17蒲江黄泥田5.935.5—7.12.80%11.58—15黔江黄泥土5.055.0左右1.89%16砂灰泡土5.775.8左右6.15%20砂灰汤泥6.456.1—6.87.01%石骨子夹泥田6.956.5—7.61.14%18石灰红泥土7.406.8—8.22.99%1816— 20石柱黄泥土6.335.0—6.42.43%水富黄泥6.256.0左右5.46%12酸紫黄泥土5.874.9—6.21.27%1916—22酸紫泥田5.774.7—6.41.61%18酸紫砂土4.754.5—6.01.22%11.510—13天全泡土5.575.0— 6.70%129—15铜梁大泥田7.036.5—7.82.36%2321—25万源泡土6.735.3—6.93.95%13西昌红泥土5.375.0-5.54.08%10西昌黄红泥土5.505.0—6.42.49%18下湿紫泥田5.405.6—7.52.70%22.521—24祥云青泥田7.807.1—7.84.99%16小土黄泥田6.105.0—6.92.14%12.59—16新津紫砂土8.107.5—8.61.02%9.58—11窑泥土5.505.0左右/5.22.38%1513仪陇黄砂土5.805.0-6.21.12%10越州红泥5.305.0—6.52.19%118—14昭通黄泥6.506.5左右1.80%1312—14砖红紫石骨土8.057.8—8.40.98%1613—19紫潮泥土7.737.6—8.11.38%18紫潮砂泥土7.931.77%12.511—14紫砂土7.136.6—7.41.94%10棕黑汤泥6.706.0-7.023.98%35棕黑汤土7.477.0-7.513.95%30棕灰泡泥6.035.4— 6.14.34%20棕灰泡土6.175.5—6.54.60%2622—30棕灰汤土5.655.8左右4.44%17棕紫泥土8.037.8—8.51.15%20.518—23棕紫砂泥土8.177.7—8.31.00%20棕紫砂土8.107.8— 8.50.74%15棕紫石骨土7.807.7—8.50.83%1817—19域（川东）暗紫砂泥土6.706.6—7.41.63％20.518—23暗紫石骨土6.906.6—7.41.55%18.516—21北碚黄泥土4.885.0—6.45.76%12.510—15扁石黄砂土5.375.2—6.42.58%1918—20大足黄泥土4.464.0—6.34.77%10钙紫大泥土7.877.6—8.51.45%2120—22钙紫石骨土8.157.6—8.61.06%16红棕石骨土8.207.8—8.50.90%20.519—22红棕紫泥土8.407.8—8.61.28%23.520—27红棕紫砂泥土8.307.8— 8.41.02%黄潮砂泥土6.876.1—7.32.32%15.514—17黄紫泥田6.686.0— 7.41.83%1813—23灰棕潮泥土8.077.5—8.01.78%20灰棕潮砂泥土7.937.1—8.41.07%12灰棕黄紫泥土6.506.7—7.31.23%2218—26灰棕紫砂泥土6.556.6—7.41.16%20.519—22剑阁紫泥土7.576.6—7.51.53%18.516—21麻砂黄棕泥5.705.5—6.57.72%21.8名山白鳝泥土4.975.0左右1.59%8.5内家紫砂泥田7.907.0—8.12.04%1712—22彭山灰棕砂土8.057.5—8.61.16%7蓬安死黄泥田5.905.7—7.02.22%1513—17破黄泥田5.404.8— 6.41.85%16黔江黄泥土5.055.0左右1.89%16石骨子田8.007.8—8.31.44%17石灰黄泥土7.736.9—8.12.26%17.512—23石柱黄泥土6.345.0—6.42.43%14酸紫黄泥土5.874.9—6.21.27%1916—22酸紫砂土4.754.5—6.01.22%11.510—13铜梁大泥田7.036.5—7.82.36%2321—25下湿黄泥田5.935.3—6.53.35%24新津紫砂土8.107.5—8.61.02%9.58—11鸭屎黄泥田7.656.6—8.03.59%2218—26仪陇黄砂土5.805.0-6.21.12%10紫潮砂泥土7.931.77%12.511—14统计结果几何均值（GM）6.50/2.0%16.0/几何标准偏差（GSD）1.11/3.3%5.8/15T/CSES141—2024T/CSES141—2024附 录 D（资料性）物种敏感性分布曲线拟合函数公式常用的物种敏感性分布曲线拟合函数公式如下：BurrIII型函数：y1 bckx 式中：y—累积概率，%；x—毒性值，mg/kg；b、c、k为函数的三个参数。Log-normal型函数：ylnx 式中：y—累积概率，%；x—毒性值，mg/kg；μ和σ为函数参数。Log-logistic型函数：y1x1x 式中：y—累积概率，%；x—毒性值，mg/kg；α、β、γ为函数参数。型函数：y1e式中：y—累积概率，%；x—毒性值，mg/kg；α和β为函数参数。Gamma型函数：x/y式中：y—累积概率，%；x—毒性值，mg/kg；α和β为函数参数。16附 录 E（资料性）拟合优度评价模型决定系数（coefficientofdetermination，R2）∑𝑛𝑛

(𝑦𝑖𝑖−𝚤𝚤)2𝑅𝑅2=1− 𝑖𝑖=1 ∑𝑛𝑛

𝑦𝑦2∑𝑛𝑛∑𝑖𝑖=1

�𝑦𝑦−𝑖𝑖=1�𝑖𝑖𝑛𝑛𝑖𝑖式中：R2为决定系数，取值范围是[0,1]，R2越接近1，说明拟合优度越大，模型拟合越精准。𝑦𝑖𝑖为实测值，𝚤𝚤为预测值，n为样本数量。均方根（roommeansquareerrorsRMSE）∑𝑛𝑛

(𝑦𝑖𝑖−𝚤𝚤)2𝑅𝑅𝑀𝑀𝐸𝐸𝐸𝐸=�𝑖𝑖=1 𝑛𝑛RM