1.一种面向园区多能源系统的潮流优化计算方法，其特征在于：包括以下步骤：
步骤1、从多能源系统运行成本、可再生能源发电量和系统网损三个维度，构建园区多能源系统潮流优化目标函数；
其中，园区多能源系统的运行成本最小的目标函数为：
min{cePe.l+cgPg.l}
园区多能源系统内可再生能源发电量最大的目标函数为：
max{Ppv+Pwt}
园区多能源系统的网损最小的目标函数为：
min{cePe.loss+cgPg.loss}
上式中，Pe.l和Pg.l分别表示园区多能源系统与电力网络和天然气网络联络线交换功率；ce和cg表示电力和天然气价格；Ppv和Pwt分别表示园区多能源系统内光伏与风力的发电量，且有Pr＝Ppv+Pwt；Pe.loss和Pg.loss分别表示园区多能源系统内电力网络和天然气网络的网损。
步骤2、在所述园区多能源系统潮流优化目标函数的基础上，建立园区多能源系统潮流约束条件；
步骤3、在所述园区多能源系统潮流约束条件下，通过多目标优化算法求解满足园区多能源系统潮流优化目标的电气和天然气网络潮流。
2.根据权利要求1所述的一种面向园区多能源系统的潮流优化计算方法，其特征在于：所述步骤2的园区多能源系统潮流约束条件由等式约束和不等式约束构成；其中，等式约束包括三相电气潮流方程约束和天然气潮流方程约束；不等式约束包括电气网络安全约束和天然气网络安全约束；
其中，三相电气潮流方程约束为：fe(P,Q,V,θ)＝0；
天然气潮流方程约束为：fg(M,p)＝0；
电气网络安全约束：fe(V)≤0；
天然气网络安全约束：fg(p)≤0；
上式中，P，Q是电气系统节点有功与无功功率；V，θ代表电气系统节点电压幅值和相角；M代表天然气流量；p代表天然气压力。