1.一种铝水反应制氢装置的故障检测方法，包括：步骤S1：建立铝水反应制氢装置结构模型；

步骤S2：建立铝水反应制氢装置的水滴蒸发和反应模型；

步骤S3：基于水滴蒸发和反应模型建立铝水反应制氢装置的传热模型；

步骤S4：基于传热模型和反应器结构模型仿真计算铝水制氢过程中铝水反应制氢装置管壁温度分布变化；

步骤S5：基于管壁温度分布变化和铝水反应制氢装置结构参数预测铝水反应制氢装置故障点；

其中，步骤S3具体包括建立铝水反应制氢装置参考坐标系，根据水滴蒸发和反应模型中水滴的上升高度对铝水反应制氢装置进行分区，获得内热源区、内热源衰减区和非内热源区；步骤S5之后还包括统计历史制氢过程中各时刻的铝水反应制氢装置底部杂质历史高度和历史产氢量，基于铝水反应制氢装置底部杂质历史高度和历史产氢量预测杂质超量故障。

2.根据权利要求1所述的铝水反应制氢装置的故障检测方法，其特征在于：所述基于铝水反应制氢装置底部杂质历史高度和历史产氢量预测杂质超量故障包括；

拟合各时刻的铝水反应制氢装置底部杂质历史高度和历史产氢量之间的关系；

获取当前制氢过程的氢产量需求，基于各时刻的铝水反应制氢装置底部杂质历史高度和历史产氢量之间的关系预测当前制氢过程的杂质预测高度；

若杂质预测高度大于杂质高度阈值，则预测当前制氢过程会发生杂质超量故障。

3.根据权利要求1所述的铝水反应制氢装置的故障检测方法，其特征在于：所述步骤S2具体包括，统计历史铝水制氢反应过程中反应水进水管的温度变化，确定制氢反应时间阈值 ，根据铝水制氢反应时间 计算单位体积的内热源散发的热量：其中， 为铝水制氢反应时间 的内热源单位体积散发的热量， 为制氢反应时间阈值， 为单个水滴在第一时间段内产生的热量， 为单个水滴在第二时间段内产生的热量， 、 分别为第一时间段和第二时间段进水管阀门开关系数， 为进水管的体积，为单个水滴的体积， 为铝水反应制氢装置圆柱体底面半径， 为第一时间段内水滴能够到达的高度， 为水滴上升高度衰减系数。

4.根据权利要求1所述的铝水反应制氢装置的故障检测方法，其特征在于：所述对铝水反应制氢装置进行分区，具体包括：内热源区圆柱体下底面圆心坐标为 ，底面半径为 ，上顶面圆心坐标为，其中， ， ， ， 为内热源区半径系数，内热源区的温度为 ， 、b分别为高度分区阈值， ；

内热源衰减区圆柱体下底面圆心坐标为 ，底面半径为 ，上顶面圆心坐标为 ，其中 ， ， ， 为内热源衰减区半径系数，内热源衰减区的温度为 ， ， ， 为内热源衰减区第i点的半径；hmin为单个水滴上升最小高度，hmax为单个水滴上升最大高度；R为圆柱体底面半径；

非内热源区为铝水反应制氢装置中除内热源区、内热源衰减区的其他区域。

5.根据权利要求1所述的铝水反应制氢装置的故障检测方法，其特征在于：建立铝水反应制氢装置的传热模型具体包括：根据水滴蒸发和反应模型计算内热源区的热量 、内热源衰减区释放的热量 ；

根据内热源区的热量、内热源衰减区释放的热量计算传递至铝水反应制氢装置管壁的总热量 ，其中 为非内热源区吸收的热量；

根据铝水反应制氢装置壁面温度分布计算铝水反应制氢装置管壁的总热量，得到传热模型：

其中， 传递至铝水反应制氢装置管壁的总热量， 、 分别为内热源区和内热源衰减区的平均比定压热容， 、 、 分别为内热源区、内热源衰减区和管壁的温度， 为热量传递损耗。

6.根据权利要求1所述的铝水反应制氢装置的故障检测方法，其特征在于：所述步骤S4具体包括：

导入步骤S3建立的传热模型和步骤S1建立的反应器结构模型，读入初始条件参数和温度场初值；计算内热源区、内热源衰减区的热源放热功率以获得内热源区和内热源衰减区的温度；对反应器结构模型进行有限元剖分，基于内热源区和内热源衰减区的温度迭代计算铝水反应制氢装置管壁温度分布变化。

7.根据权利要求6所述的铝水反应制氢装置的故障检测方法，其特征在于：对反应器结构模型进行有限元剖分具体包括：铝水反应制氢装置模型整体为第一区域，铝水反应制氢装置内放置反应物的区域为第二区域，所述第一区域包括所述第二区域；

对第一区域进行有限元网格剖分，根据第一约束条件对第二区域进行有限元网格剖分优化；

对优化后的第二区域的网格进行检查，若通过检查，则完成铝水反应制氢装置制氢的建模，否则，返回根据第一约束条件对第二区域进行有限元网格剖分优化；

所述第一约束条件为网格尺寸大小与网格距离水滴位置远近的函数；三维空间内函数在y轴方向上的自变量变化对网格尺寸大小的影响大于在x轴方向上的自变量及在z轴方向上的自变量的变化对网格尺寸大小的影响。

8.根据权利要求7所述的铝水反应制氢装置的故障检测方法，其特征在于：三维空间内第一约束条件为：

其

中，Smax为设定的面积最大值，Smin为设定的面积最小值，S(i,j,k)为面积约束公式，i0为水滴位置在x轴方向上的像素坐标值，j0为水滴位置在y轴方向上的像素坐标值，k0为水滴位置在z轴方向上的像素坐标值，i,j,k为节点坐标，n为指数系数，a1、a2、a3为匹配系数， 为范围系数。

9.根据权利要求1所述的铝水反应制氢装置的故障检测方法，其特征在于：所述步骤S5具体包括：

基于铝水反应制氢装置结构参数计算铝水反应制氢装置管壁各个位置能够承受的最大温升 ；

基于铝水反应制氢装置管壁温度分布变化计算铝水制氢反应过程管壁各个位置的实际温升 ；

对比管壁各个位置能够承受的最大温升 和实际温升 ，若 ，则位置k为预测故障点，其中 为铝水反应制氢装置管壁位置k处能够承受的最大温升，为铝水制氢反应过程管壁位置k处的实际温升， 为温升承受阈值。

10.一种基于铝水制氢过程仿真的故障预测系统，包括：模型建立模块，用于建立铝水反应制氢装置结构模型，建立铝水反应制氢装置的水滴蒸发和反应模型，并基于水滴蒸发和反应模型建立铝水反应制氢装置的传热模型；

仿真计算模块：用于基于传热模型和反应器结构模型仿真计算铝水制氢过程中铝水反应制氢装置管壁温度分布变化；

故障预测模块，用于基于管壁温度分布变化和铝水反应制氢装置结构参数预测铝水反应制氢装置故障点；

其中，建立铝水反应制氢装置的传热模型包括建立铝水反应制氢装置参考坐标系，根据水滴蒸发和反应模型中水滴的上升高度对铝水反应制氢装置进行分区，获得内热源区、内热源衰减区和非内热源区；故障预测模块还被配置为用于统计历史制氢过程中各时刻的铝水反应制氢装置底部杂质历史高度和历史产氢量，基于铝水反应制氢装置底部杂质历史高度和历史产氢量预测杂质超量故障。