1.一种多能耦合系统的潮流计算方法，其特征在于，包括：
构建天然气网络模型、供热网络模型、电力网络模型和耦合设备模型；
根据所述天然气网络模型、所述供热网络模型、所述电力网络模型和所述耦合设备模型，得到多能耦合系统的稳态模型；
根据所述多能耦合系统的所述稳态模型，得到所述多能耦合系统的非线性代数方程；
根据所述非线性代数方程和连续牛顿法，得到所述多能耦合系统的五阶精度解和四阶精度解；
根据所述五阶精度解和所述四阶精度解，得到所述多能耦合系统的潮流；
所述根据所述多能耦合系统的所述稳态模型，得到所述多能耦合系统的非线性代数方程，包括：
根据所述多能耦合系统的所述稳态模型，得到所述多能耦合系统的相关变量；
根据所述相关变量，得到所述非线性代数方程；
所述非线性代数方程为：  ；
其中， ，为所述相关变量，所述相关变量包括所述多能耦合系统温度、流量、功率以及电压；
所述根据所述非线性代数方程和连续牛顿法，得到所述多能耦合系统的五阶精度解和四阶精度解包括：
根据所述非线性代数方程，得到非线性函数；
所述非线性函数为： ；
根据所述非线性函数，通过所述连续牛顿法，得到所述多能耦合系统的所述五阶精度解和所述四阶精度解；
其中，所述五阶精度解为：
；
其中， 为所述五阶精度解， 为中间变量， ， 为系数， 为第 次迭代
的解， 为迭代步长；
所述四阶精度解为：
；
其中， 为所述四阶精度解， 为第 次迭代的解， 为所述迭代步长， 为系数， 为所述中间变量；
所述根据所述五阶精度解和所述四阶精度解，得到所述多能耦合系统的潮流包括：
根据所述五阶精度解与所述四阶精度解的差值，调整所述迭代步长；
通过所述迭代步长，得到所述多能耦合系统的所述潮流；
所述根据所述五阶精度解与所述四阶精度解的差值，调整所述迭代步长，包括：
对所述五阶精度解与所述四阶精度解的差值进行归一化处理，得到所述多能耦合系统的潮流计算误差表达式：
；
其中， 代表 范数， 为第 次迭代的解， 为迭代步长， 为所述五阶精度解与所述四阶精度解的差值；
当 小于预设第二误差限 时，则减小所述迭代步长，得到减小后的所述迭代步长；
所述减小后的所述迭代步长为：
；
其中， 为所述减小后的所述迭代步长， 为所述迭代步长的最小值， 为所述迭代步长；
并重复根据所述非线性代数方程和连续牛顿法，得到所述多能耦合系统的五阶精度解和四阶精度解；
当 大于所述预设第二误差限 时，则增大所述迭代步长，得到增大后的所述迭代步长；
所述增大后的所述迭代步长为：
；
其中， 为所述增大后的所述迭代步长， 为所述迭代步长；
当增大所述迭代步长的次数大于预设次数时，得到所述多能耦合系统的所述潮流。
2.根据权利要求1所述的多能耦合系统的潮流计算方法，其特征在于，所述构建天然气网络模型、供热网络模型、电力网络模型和耦合设备模型，包括：
根据气流方程、节点质量流量守恒方程以及压缩机模型，得到所述天然气网络模型；
所述气流方程为：
；
其中， 为所述多能耦合系统的节点 与节点 之间的管道内的气体流量； 为管道系数； 为所述节点 的压强； 为所述节点 的压强； 为符号函数，所述符号函数满足 ； 为所述管道的长度； 为所述管道的直径； 为所述
管道的摩擦系数； 为气体温度； 为气体压缩系数； 为天然气相对密度；
所述节点质量流量守恒方程为：
；
其中， 为节点 的注入气体质量流量， 为所述管道 的气体质量流量， 为流入所述节点 的管道集合；
所述压缩机模型为：
；
其中， 为所述多能耦合系统的压缩机的压缩比， 为所述压缩机的高压侧的压强， 为所述压缩机的低压侧的压强。
3.根据权利要求2所述的多能耦合系统的潮流计算方法，其特征在于，所述构建天然气网络模型、供热网络模型、电力网络模型和耦合设备模型，还包括：
根据水力模型和热力模型，得到所述供热网络模型；
所述热力模型包括管道温度传输方程、节点温度混合方程以及节点换热方程，所述水力模型包括节点质量流量连续方程和水头损失方程；
其中，所述管道温度传输方程为：
；
其中， 为所述管道的入口温度， 为所述管道的出口温度， 为环境温度， 为水的比热容， 为所述管道的质量流量，为所述管道的热阻系数， 为管道长度；
所述节点温度混合方程为：
；
；
其中， 为所述节点 的温度， 为所述管道 的质量流量， 为所述管道 的所述出口温度， 为所述管道 的所述入口温度， 为流入所述管道的节点 的管道集合，为流出所述管道的节点 的管道集合；
所述节点换热方程为：
；
其中， 为所述节点产生/消耗的热功率， 为节点供水温度， 为节点回水温度，为所述水的比热容；
所述节点质量流量连续方程为：
；
其中， 为节点 的注入质量流量， 为管道 的质量流量， 为流入所述节点 的管道集合；
所述水头损失方程为：
；
其中， 为所述管道的压降系数， 为所述管道 的质量流量， 为环路 中的所述管道集合。
4.根据权利要求1所述的多能耦合系统的潮流计算方法，其特征在于，所述构建天然气网络模型、供热网络模型、电力网络模型和耦合设备模型，包括：
根据节点功率注入方程和节点电压方程，得到所述电力网络模型；
所述节点功率注入方程为：
；
；
其中， 为所述节点 的注入有功功率， 为所述节点 的注入无功功率， 为所述节点的电压幅值， 为所述节点 的电压幅值， 为所述节点 和所述节点 之间的导纳的实部， 为所述节点 和 之间的导纳的虚部； 为所述节点 和所述节点 之间的电压相角差，的范围为1，2，3……n；
所述节点电压方程为：
；
；
其中， 为所述节点 的电压幅值， 为所述电压幅值的给定值， 为平衡节点的电压相角， 为所述电压相角的给定值， 为电源电压与所述平衡节点的集合。
5.根据权利要求2所述的多能耦合系统的潮流计算方法，其特征在于，所述构建天然气网络模型、供热网络模型、电力网络模型和耦合设备模型，还包括：
根据燃机轮机模型、锅炉模型、热电联产机组模型和循环泵模型，得到所述耦合设备模型；
所述燃机轮机模型为：
；
其中， 为燃气轮机消耗的天燃气流量， ， 与 为所述燃气轮机的能耗系数， 为天然气的净热值， 为所述燃气轮机的输出电功率；
所述锅炉模型为：
；
其中， 为锅炉消耗的天燃气流量， 为所述锅炉输出的热功率， 为所述锅炉的功率转换效率， 为所述天然气的净热值；
所述热电联产机组模型为：
；
其中， 为热电联产机组消耗的天然气流量， 为天然气的净热值， 为所述热电联产机组输出的热功率， 为所述热电联产机组输出的电功率， 为所述热电联产机组的功率转换效率；
所述循环泵模型为：
；
其中， 为循环泵的消耗的电功率， 为所述循环泵的入口质量流量， 为重力加速度， 为水头压力， 为循环泵效率。
6.一种多能耦合系统的潮流计算装置，其特征在于，包括：
稳态模型构建模块，用于构建天然气网络模型、供热网络模型、电力网络模型和耦合设备模型；根据所述天然气网络模型、所述供热网络模型、所述电力网络模型和所述耦合设备模型，得到多能耦合系统的稳态模型；
潮流计算模块，用于根据所述多能耦合系统的所述稳态模型，得到所述多能耦合系统的非线性代数方程；根据所述非线性代数方程和连续牛顿法，得到所述多能耦合系统的五阶精度解和四阶精度解；根据所述五阶精度解和所述四阶精度解，得到所述多能耦合系统的潮流；
所述根据所述多能耦合系统的所述稳态模型，得到所述多能耦合系统的非线性代数方程，包括：
根据所述多能耦合系统的所述稳态模型，得到所述多能耦合系统的相关变量；
根据所述相关变量，得到所述非线性代数方程；
所述非线性代数方程为：  ；
其中， ，为所述相关变量，所述相关变量包括所述多能耦合系统温度、流量、功率以及电压；
所述根据所述非线性代数方程和连续牛顿法，得到所述多能耦合系统的五阶精度解和四阶精度解包括：
根据所述非线性代数方程，得到非线性函数；
所述非线性函数为： ；
根据所述非线性函数，通过所述连续牛顿法，得到所述多能耦合系统的所述五阶精度解和所述四阶精度解；
其中，所述五阶精度解为：
；
其中， 为所述五阶精度解， 为中间变量， ， 为系数， 为第 次迭代
的解， 为迭代步长；
所述四阶精度解为：
；
其中， 为所述四阶精度解， 为第 次迭代的解， 为所述迭代步长， 为系数， 为所述中间变量；
所述根据所述五阶精度解和所述四阶精度解，得到所述多能耦合系统的潮流包括：
根据所述五阶精度解与所述四阶精度解的差值，调整所述迭代步长；
通过所述迭代步长，得到所述多能耦合系统的所述潮流；
所述根据所述五阶精度解与所述四阶精度解的差值，调整所述迭代步长，包括：
对所述五阶精度解与所述四阶精度解的差值进行归一化处理，得到所述多能耦合系统的潮流计算误差表达式：
；
其中， 代表 范数， 为第 次迭代的解， 为迭代步长， 为所述五阶精度解与所述四阶精度解的差值；
当 小于预设第二误差限 时，则减小所述迭代步长，得到减小后的所述迭代步长；
所述减小后的所述迭代步长为：
；
其中， 为所述减小后的所述迭代步长， 为所述迭代步长的最小值， 为所述迭代步长；
并重复根据所述非线性代数方程和连续牛顿法，得到所述多能耦合系统的五阶精度解和四阶精度解；
当 大于所述预设第二误差限 时，则增大所述迭代步长，得到增大后的所述迭代步长；
所述增大后的所述迭代步长为：
；
其中， 为所述增大后的所述迭代步长， 为所述迭代步长；
当增大所述迭代步长的次数大于预设次数时，得到所述多能耦合系统的所述潮流。
7.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5中任一项所述的多能耦合系统的潮流计算方法。