1.一种城市综合能源网统一平台潮流计算方法，其特征在于，具体步骤包括：
s1：了解综合能源系统中各类异质能源并建立各类异质能源耦合至系统中的相应模型，所述模型包括热电联产、燃气轮机组、换热器、电蓄热、电制氢和电动汽车等，其中各类异质能源以节点分布式电源或负荷的形式接入系统中；
所述热电联产机组的模型根据热电比是否为一个确定值，可将其分为如下两种形式：
式中：pchp为热电联产机组的电出力，hchp为热电联产机组的热出力，cm为定热电比，cz为变热电比，ηe为热电联产机组的冷凝效率，fin为燃料输入速率；这里需要注意的是cz虽然是一个变化的值，但是可以根据实际情况，在某个时段内会保持不变，在通常情况下，热电联产机组以定热电比运行；所述热电联产机组在电热耦合的节点处电功率以及热功率的分配，实现了两种不同种类能源电能和热能的耦合；
所述燃气轮机就是一种常见的热电联产单元，燃汽轮机组通过将天然气管道末端节点的天然气转化为电能实现了电气耦合，并为电力系统提供有功与无功；所述燃气轮机可以视作固定热电比的情况如方程(1)所示，所述燃气轮机的工作特性可用方程(3)描述：
其中l为天然气流量，q为天然气热值，单位为mj/m3，η为转换效率；
建立换热器模型，所述换热器在传热过程中，其特征温差δtave和换热量q之间满足类似欧姆定律的关系，即：
其中特征温差δtave为换热器两端温度差，可视为换热器两端冷热流体的进出口平均温度的差值；r为换热器的耗散热阻，一般采用逆流换热器时，r的表达式为：
式中：gu为用户端循环水的热容量流，是循环水的质量流量mu和热容cp的乘积；au＝ka/gu；gh为储热装置中循环水的热容量流，是储热装置中循环水质量流量mh和热容ch的乘积；ah＝ka/gh，k和a分别为换热器的换热系数和换热面积；
建立电蓄热与电制氢模型，所述电蓄热与电制氢模型都是储能装置，电蓄热模型的热储能过程将电功率转换为热功率，一般处于电节点末端等效为一个储热的电负荷，电蓄热模型处于放热过程时位于电热耦合节点处，产生的热能为其他热节点供能；电蓄热模型的储热和放热过程可以类比做电容的充放电过程：
对于电容的充放电过程来说：
式中：i为电流，du/dt为电压对时间的变化率，c为电容；
那么对于电蓄热模型的热储能过程来说：
其中，m为储热装置中热水的总质量，cp为水的比热容，q为换热量，电蓄热的放热过程中可用电热转换系数ceh来描述：
式中：peh和heh分别为电蓄热的放热过程中的电功率和热功率；
电制氢模型中储氢环节可近似看作电解水的过程，近似视为“电容”或者恒功率源；氢储能中用氢环节可以视为注入的功率源；
建立电动车模型，所述单个电动车充电过程可简化视为恒功率模型，即忽略暂态过程，将电动车模型视为pq节点处的固定负荷；
s2：将所述各类异质能源模型统一纳入电力系统的分析决策体系，即考虑所述各类异质能源模型在电力系统潮流分析计算中所扮演的角色，建立综合能源网统一至电平台的同质化等效模型；
对配电网各母线节点而言，所述热电联产机组、电蓄热的放热过程和电制氢的用氢环节等电、气、热过程均可视为注入功率源，电动汽车、换热器、电蓄热的储热过程、电制氢的储氢环节等各类储能过程可以视为电负荷；
s3：对于断面潮流而言，忽略电动汽车储能的暂态过程，可以将所述电动汽车视为容量已知的电负荷，因此可以根据功率源等效的不同将恒定有功与无功功率的节点视为pq节点；将已知有功功率和电压幅值的节点视为pv节点，即所述燃气轮机可根据分轴结构燃气轮机中对励磁方式的控制方式不同可作不同的等效，当控制方式为电压控制时可等效为pv节点，当控制方式为功率因数控制时可等效为pq节点；电制氢的用氢环节，做燃料电池时被等效为pv节点；所述热电联产机组，一般通过燃气轮机组实现并用控制功率因数的方式，因此视为pq节点；
s4：对于某一确定节点各种异质能源的具体情况分析计算时，可将换热器视为“电阻”，储热装置视为“电容”，从而实现模型简化进而可以构建系统整体的能量平衡方程，简单电热互联系统的整体约束方程如下：
ptp(t)+pchp(t)+pw(t)＝pel(t)+peh(t)(9)
式中：ptp(t)为火电机组电出力，pchp(t)为热电联产机组电出力，pw(t)为风机电出力，pel(t)为系统原有电负荷，peh(t)为引入电加热装置而新增的电负荷，公式(10)中m为储热装置中热水的总质量，cp为水的比热容；tu,i为热用户的回水温度，td,o(t)为储热装置水温，r为换热器等效耗散热阻；
s5：对所述在简单电热互联系统的整体约束方程下的综合能源网统一至电平台的同质化等效模型采用前推回代法进行分析计算，传统的前推回代法只适用于配电网pq节点的潮流计算，对于模型简化而带来的pv节点问题，需要对算法进行改进，这里采用影响因子矩阵的修正方法，由于pv节点有功功率与电压幅值已知，因此可通过无功补偿改变无功功率进而使电压相角满足要求，可以发现改变某一节点无功补偿，各节点的电压变化量与无功补偿量呈近似于线性关系：
对上式进行简化：
其中节点2为pv节点，v1和v2分别为两个节点的电压，δq2为节点2处注入的无功补偿，q2为节点2初始无功功率，p2为节点2初始有功功率，补偿后电压幅值相应变为v2′，δv2为补偿后的节点2的电压差，r12+jx12为节点间的转移阻抗；可以看出，无功补偿量与电压变化量的比值为一常数，以此为基础构建影响因子矩阵i，其中：
iij为影响因子矩阵对应第i行第j列元素，通过对所有pv节点无功功率进行修正，获得各节点电压修正量，可以写成：
δu＝iδq(15)
δu为pv节点电压不平衡向量；δq为pv节点无功补偿量向量；i表示影响因子矩阵，阶数即为网络中pv节点个数；
s6：判断所有的电压是否满足收敛判定条件，若不满足则对注入无功功率进行修正，pv节点的收敛判定条件为：
式中
s7：迭代至各节点均满足精度要求，计算结束，输出潮流计算的相关结果；否则返回至s5重新进行计算直到各节点均满足精度要求。
2.根据权利要求1所述城市综合能源网统一平台潮流计算方法，其特征在于，对于所述前推回代法中的产生的弱环网问题，采用解环的处理方法，对于负荷侧按固定比例分配的环网，可以按照这个比例分配负荷进行直接解环；若这个比例系数无法确定，即在环网选择一个节点，在此节点断开使网络变为辐射型，然后采用前面所述的前推回代法进行潮流计算；
显然对于断开的节点处：
式中：p1与p2为解环后新生成的两个节点分别注入的有功功率，q1与q2为解环后新生成的两个节点分别注入的无功功率，p和q是原环网支路上流动的有功与无功。
3.根据权利要求1所述城市综合能源网统一平台潮流计算方法，其特征在于，所述前推回代法中对影响因子矩阵进行拟合时，拟合曲线虽然拟合度较高，但是无功变化量对影响因子矩阵较为敏感，因此需要尽可能保证拟合精度，采用最小二乘法对对无功在500-1000kvar附近取点拟合，所得斜率即为影响因子矩阵对应数据；
最小二乘法计算公式如下：其中x，y为自变量与因变量，k为拟合斜率。
4.根据权利要求1所述城市综合能源网统一平台潮流计算方法，其特征在于，前推回代法包括前推和回代两个过程，首先各节点电压赋初值，然后进行功率前推：从末端节点开始功率前推，对每条支路，由受端节点的注入功率和支路阻抗数据计算送端前一节点支路的功率，直到完成首节点功率的计算时停止，功率变化量如下：
其中p、q为分别为节点的注入有功、无功功率；u为节点电压；r+jx为支路阻抗；
电压回代：从首节点开始，对每条支路，由送端节点的节点电压和支路功率计算受端节点的节点电压，直到完成所有末端节点电压的计算时停止；电压变化量实部虚部如下：