【水文模型】评价指标

WW、forever

已于 2023-12-04 16:40:22 修改

阅读量1.1w

点赞数 23

分类专栏： # 水文模型文章标签：水文模型

于 2023-02-12 21:31:46 首次发布

本文链接：https://blog.csdn.net/qq_44246618/article/details/128999913

版权

水文模型专栏收录该内容

35 篇文章 181 订阅

订阅专栏

模型模拟结果的评价主要通过有些评价指标结合时间序列曲线图来进行。常用的评价指标包括均方根差（Root Mean Square Error, RMSE）、均方误差（Mean Square Error, RSE）、平均绝对误差（Mean Abosolute Error, MSE）、相关系数（Correlation Coefficient, CC）、纳什效率系数（Nash-Sutcliffe efficiency coefficient, NSE）和偏差百分比（Percent bias, Pbias）等，用于检测模型的系统误差和随机误差以及实测值与模拟值之间的关系或偏差程度。

1 皮尔逊相关系数（Pearson’s correlation coefficient, PCC）

1.1 简介

在这里插入图片描述

1.2 MATLAB实现代码

MATLAB实现代码如下：

function PCC = GetPCC(Qs, Qo)
% 输入变量
% Qs    模拟径流
% Qo   实测径流
% 输出变量
% PCC       皮尔逊相关系数（Pearson’s correlation coefficient, PCC）

if length(Qs)==length(Qo)
    COV = cov(Qs, Qo);
    PCC = COV(1,2)/std(Qs)/std(Qo);
else
    error("实测径流和模拟径流长度不等");
end

end

2 百分比偏差（Percent bias, Pbias）

2.1 简介

百分比偏差Pbias（%），又称相对误差，表示测量模拟的水流的百分比偏差大于或小于相应的推断的自然水流，越接近0模型效果越好：
在这里插入图片描述
说明： 体现模拟值和实测值的累积偏差，当模拟水文过程与实际趋势贴合情况良好时，PBIAS可更加精确地评估模型总水量平衡的效果。最优值为0，表示模型模拟准确，正值表示模型低估观测值的偏差，负值表示模型高估观测值的偏差。

2.2 MATLAB实现代码

MATLAB实现代码如下：

function Pbias = GetPbias(Qs, Qo)
% 输入变量
% Qs    模拟径流
% Qo   实测径流
% 输出变量
% Pbias     百分比偏差/相对误差（Percent bias, Pbias）

if length(Qs)==length(Qo)
    Qs_avg = mean( Qs );
    Qo_avg = mean( Qo );
    Pbias = (Qs_avg -Qo_avg )/Qo_avg *100;
else
    error("实测径流和模拟径流长度不等");
end

end

3 纳什效率系数（the Nash-Sutcliffe efficiency coefficient, NSE）

3.1 简介

在这里插入图片描述
说明： 体现模拟径流和实测径流的拟合程度，但会赋予洪峰段更高的计算权重，容易忽视平水期或枯水期的拟合情况

3.2 MATLAB实现代码

MATLAB实现代码如下：

function NSE = GetNSE(Qs, Qo)
% 输入变量
% Qs    模拟径流
% Qo   实测径流
% 输出变量
% NSE       纳什效率系数（Nash-Sutcliffe efficiency coefficient, NSE）

if length(Qs)==length(Qo)
    QoAve = mean(Qo);
    NSE = 1- sum( (Qs-Qo).^2 )/sum( (Qo-QoAve).^2 );
else
    error("实测径流和模拟径流长度不等");
end

end

4 克林-古普塔效率系数（Kling-Gupta efficiency coefficient, KGE）

4.1 简介

KGE测量的是一个点到最优点的欧几里得距离，相关系数（correlation coefficient, PCC）、偏置比（bias ration, BR）和相对变化率（relative variablity, RV）等于1：
取值范围：（-∞，1）
在这里插入图片描述

4.2 MATLAB实现代码

MATLAB实现代码如下：

function KGE = GetKGE(Qs, Qo)
% 输入变量
% Qs    模拟径流
% Qo   实测径流
% 输出变量
% KGE       克林-古普塔效率系数（Kling-Gupta efficiency coefficient, KGE）

if length(Qs)==length(Qo)
    QsAve = mean(Qs);
    QoAve = mean(Qo);
    COV = cov(Qs, Qo);
    CC = COV(1,2)/std(Qs)/std(Qo);
    BR = QsAve/QoAve;
    RV = ( std(Qs)/QsAve )/ ( std(Qo)/QoAve );
    KGE = 1- ( (CC-1)^2+(BR-1)^2+(RV-1)^2 )^0.5;
else
    error("实测径流和模拟径流长度不等");
end

end

5 决定系数R^2(Coefficient of Determination)

5.1 简介

在这里插入图片描述

说明： 模拟值与实测值的线性相关程度，弊端是R2对模拟值整体偏高或偏低的偏差响应不明显

5.2 MATLAB实现代码

MATLAB实现代码如下：

function R2 = GetR2(Qs, Qo)
% 输入变量
% Qs    模拟径流
% Qo   实测径流
% 输出变量
% R2         决定/相关系数R2（Coefficient of Determination）

if length(Qs)==length(Qo)
    QsAve = mean(Qs);
    QoAve = mean(Qo);
    tempX = sum( (Qo-QoAve).*(Qs-QsAve));
    tempY = sum( (Qo-QoAve).^2 )^0.5 * sum( (Qs-QsAve).^2)^0.5 ;
    R2 = tempX/tempY;
else
    error("实测径流和模拟径流长度不等");
end

end

6 归一化均方根误差(NRMSE)

6.1 简介

6.2 MATLAB实现代码

7 均方根误差与观测值标准差的比率（RSR）

7.1 简介

在这里插入图片描述
说明： RSR从最优值0变化到较大的正值，该值越小，表示模型性能越好。

7.2 MATLAB实现代码

MATLAB实现代码如下：

function RSR = GetRSR(Qs, Qo)
% 输入变量
% Qs    模拟径流
% Qo   实测径流
% 输出变量
% RSR       均方根误差与观测值标准差的比率（RSR）

if length(Qs)==length(Qo)
    QoAve = mean(Qo);
    RSR = ( sum(Qo-Qs).^2 )^0.5/( sum(Qo-QoAve).^2 )^0.5;
else
    error("实测径流和模拟径流长度不等");
end

end