1.一种海底混输管道腐蚀速率预测方法，其特征在于，包括以下步骤：(0) (0) (0) (0)

1)不等时间间距测量海底混输管道的腐蚀速率检测值序列X ，X ＝{x (1),x(0) (0)(2),...,x (n)}，x (k)表示序号为k时管道腐蚀速率的检测值；

(0)

2)根据步骤1)得到的海底混输管道的腐蚀速率检测值序列X 构建用于海底混输管道腐蚀速率预测的ARIMA模型及NEGM(1,1,τ)模型；

3)将ARIMA模型和NEGM(1,1,τ)模型分别利用等权重法、灰色关联度准则及误差平方和最小准则建立三个组合预测模型；

4)根据均方根误差RMSE、平均相对误差AARD及平均误差MAPE对步骤3)建立的三个组合预测模型进行性能评估，并选择性能最优的组合预测模型，然后利用性能最优的组合预测模型进行海底混输管道腐蚀速率的预测。

2.根据权利要求1所述的海底混输管道腐蚀速率预测方法，其特征在于，建立用于海底混输管道腐蚀速率预测的NEGM(1,1,τ)模型的具体过程为：(0) (1) (1)

11)由海底混输管道的腐蚀速率检测值序列X 计算一阶累加1‑AGO序列X ，其中，X(1) (1) (1)＝{x (1),x (2),...,x (n)}， tj表示序号为j时对应的检测时间；

(0)

12)由海底混输管道的腐蚀速率检测值序列X 和一阶累加1‑AGO序列计算海底混输管(0) (1)道的腐蚀速率检测值序列X 的光滑比ρ(k)及级比σ (k) ，其中，(1)

根据计算得到的光滑比ρ(k)及级比σ (k)判断海(0)

底混输管道的腐蚀速率检测值序列X 是否满足光滑性条件，当海底混输管道的腐蚀速率(0)检测值序列X 满足光滑性条件时，则转至步骤13)中；当海底混输管道的腐蚀速率检测值(0)序列X 不满足光滑性条件时，则转至步骤1)中；

(1) (0) (1) (1) (1)

13)根据一阶累加1‑AGO序列X 计算得序列X 紧邻均值生成序列Z ，其中Z ＝(z(1) (1) (1) (1) (1)(2),z (3),…,z (n))，z (k)＝0.5(x (k)+x (k‑1)),k＝2,3,…,n；

(1)

14)根据一阶累加1‑AGO序列X 计算NEGM(1,1,τ)模型的建模初始值 其中，τ为 的时间参数，τk表示海底混输管道腐蚀速率的1‑AGO序列各分量的权重系数，(1)

15)设 为白化方程参数列，根据一阶累加1‑AGO序列Z 构建灰微分方程a及μ均为常数；

16)求 解灰微分方 程 得时间响应 函数 其中 ，并以此得到海底混输管道的腐

(0)

蚀速率检测值序列X 的还原值

则预测得到的海底混

输管道腐蚀速率预测值序列

3.根据权利要求1所述的海底混输管道腐蚀速率预测方法，其特征在于，基于残差平方和最小原则建立 优化模型，令 再代入的表达式，得S对τ的一阶导数，令S对τ的一阶导数为零，得 的时间参数τ。

4.根据权利要求1所述的海底混输管道腐蚀速率预测方法，其特征在于，建立用于海底混输管道腐蚀速率预测的ARIMA模型的具体操作过程为：(0) (0)

21)将海底混输管道的腐蚀速率检测值序列X 进行平稳化处理，得平稳序列Y ；

(0)

22)绘制平稳序列Y 的自相关函数及偏自相关函数图，根据偏自相关函数图的性质以及ARIMA模型的定阶准则确定最优的模型；

23)对步骤22)中的ARIMA模型进行参数估计，得用于海底混输管道腐蚀速率预测的ARIMA模型。

5.根据权利要求1所述的海底混输管道腐蚀速率预测方法，其特征在于，设N为海底混输管道腐蚀速率检测次数，yk为序号为k时腐蚀速率检测值， 表示序号为k时的预测值，则有