1.一种工业园区综合能源系统与工厂生产计划协同优化方法，其特征在于：包括如下步骤：步骤S1、针对工业园区内综合能源系统及工业用户的结构，在原有能量供应、转换、存储设备模型的基础上，建立了综合能源系统氢气制取模型、二氧化碳回收利用及存储模型、生产原料存储区模型；

能量的供应设备模型包括天然气网、电网、光伏发电机组、太阳能热水器；能量的转换设备模型包括燃气锅炉、燃气轮机、P2G单元、吸收式制冷机、电制冷机；能量的存储设备模型包括电能存储设备；综合能源系统氢气制取模型包括电解水单元；二氧化碳回收利用及存储模型包括二氧化碳捕获设备；生产原料存储区模型包括初始原料仓库、车间1中间产品存储仓库、车间2中间产品存储仓库、车间z中间产品存储仓库；

所述P2G单元及电解水单元模型为：

上式中： 表示电转气单元产生天然气的速率；δ表示消耗氢气的速率与生产天然气速率之间比例系数； 表示水的电解槽产生氢气的速率；

上式中： 表示电转气单元消耗二氧化碳的速率；λ表示消耗氢气的速率与生产天然气消耗二氧化碳速率之间比例系数；

上式中： 表示电转气单元消耗电能的功率；表示消耗氢气的速率与耗电功率之间转换系数；

上式中： 表示电转气单元消耗电能的最大功率；

上式中： 表示电解水槽产出氧气的速率；τ表示产出氢气的速率与产出氧气的速率之间比例系数；

上式中： 表示电转气单元产出天然气的速率；θ表示天然气的产出速率转换为天然气功率的系数；

综合能源系统二氧化碳回收利用模型为：

上式中： 表示捕获由热电联产机组及燃气锅炉产出的二氧化碳速率； 表示由于捕获技术限制被排放到空气中二氧化碳速率；ρ表示转换系数； 表示燃气锅炉产热功率； 表示燃气锅炉产热效率； 分别表示热电联产机组的产电及产热功率；

分别表示热电联产机组的产电及产热效率；

上式中： 表示回收的二氧化碳用于P2G单元消耗的速率； 表示存储到二氧化碳罐的速率；

综合能源系统二氧化碳存储模型为：

上式中：CO2s,t表示二氧化碳存储量； 表示二氧化碳初始存储量； 表示二氧化碳出售速率；

上式表示二氧化碳初末时刻二氧化碳存储量保持不变；

0≤CO2s,t≤CO2s,max，

上式中：CO2s,max表示二氧化碳最大存储容量；

上式中： 表示二氧化碳最大存储速率；

步骤S2、园区综合能源管理中心接收到工业用户的生产任务、车间能耗类型、生产约束，构建工业用户不同车间能耗的多能耦合模型及工厂生产约束模型；车间能耗类型包括纯电能车间、纯冷能车间、纯热能车间、电热耦合车间、电冷耦合车间、热冷耦合车间；

工业用户不同车间能耗的多能耦合模型为：

上式中：PE,t表示工厂消耗电能的总功率；z表示车间标号；E'表示消耗纯电能车间集合；x为二进制，表示车间运行状态，取1时代表车间工作，反之不工作；PE,z表示z车间消耗的电功率；EH'表示热点耦合车间集合； 表示电热耦合车间z消耗的电功率；EC'表示电冷耦合车间集合； 表示电冷耦合车间z消耗的电功率；

上式中：PH,t表示工厂消耗热能的总功率；H'表示消耗纯热能车间集合；PH,z表示z车间消耗的热功率； 表示电热耦合车间消耗的总热功率；HC'表示热冷耦合车间集合； 表示热冷耦合车间z消耗的热功率；

上式中：PC,t表示工厂消耗冷能的总功率；C'表示消耗纯冷能车间集合； 表示电冷耦合车间消耗的总冷功率； 表示热冷耦合车间消耗的总冷功率；

其中：

上式中： 表示z车间电热耦合系数； 表示z车间电冷耦合系数； 表示z车间热冷耦合系数；

生产约束包括不同车间的生产顺序、生产速度、运行状态，以及存储仓库的初始数量、剩余容量；

P(EHC),z表示车间z各种耦能耗功率；υ＝{υ1…υZ}表示不同车间的生产速率； 表示车间z的中间产品初始存储数量； 表示车间z中间产品存储容量；

具体的，中间产品存储仓库模型：

上式中： 表示车间z中间产品存储数量； 表示车间z中间产品初始存储量；υz,t表示车间z不同时刻加工速率；

上式表示初末时刻中间产品存储量保持不变；

上式中： 表示车间z中间产品最大存储容量；

υz,t＝xz,t·υz，

上式中：υz表示车间z加工速率；步骤S3、园区综合能源管理中心将接收到工业用户信息、综合能源系统的设备信息、外部气网及电网的价格、功率信息进行整合，在完成工厂的日生产任务前提下，以园区内综合能源系统所有设备的运行效率、工业用户不同车间的多能耦合模型及工厂生产流程为约束条件，以工业园区一天能源消耗总费为目标函数，并构建综合能源系统与工业用户运行的协同优化模型；

步骤S4、园区综合能源管理中心对综合能源系统与工业用户运行的协同优化模型求解，得到园区综合能源系统各设备出力计划及工厂生产调度安排；

步骤S5、园区综合能源管理中心按照求解结果控制综合能源系统的运行，并将求得的结果发送给工厂用户，使其按相应的安排从事生产活动。

2.根据权利要求1所述的一种工业园区综合能源系统与工厂生产计划协同优化方法，其特征在于：所述步骤S3的具体过程如下：完成工厂的日生产任务前提下，以工业园区一天能源消耗总费为目标函数，具体包括向电网购电成本、天然气站购气成本、二氧化碳排放成本、捕获二氧化碳售出的收入、电解水制氢副产品氧气售出收入；

具体的，工业园区一天能源消耗总费的目标函数表达式为：上式中：Ctotal表示工业园区一天能源消耗总费用；Ct表示工业园区一天中某时刻能源消耗总费用；

其中：

上式中： 表示电网购电功率；pE,t表示电网购电价格； 表示天然气站购气功率；

pg,t表示天然气站购气价格； 表示二氧化碳排放速率； 表示二氧化碳排放价格；

表示二氧化碳购买速率； 表示二氧化碳购买价格； 表示电解水制氢产出氧气的速率； 表示出售氧气的价格； 表示出售二氧化碳的速率； 表示出售二氧化碳的价格；

约束条件包括：电热冷功率平衡约束、电网及气网交互功率约束、各能源设备运行约束、储能约束、工厂生产约束；

综合能源系统生产电、热、冷之间的关系表达式为：

上式中： 表示园区综合能源系统提供的电、热、冷功率； 表示蓄电池充放电功率； 表示光伏产电功率； 表示电转气耗电功率； 表示电转气产气功率； 表示太阳能热水器产热功率；

其中：

上式中：ω表示百分数； 表示热电联产机组产热效率和产电效率； 表示燃EC AC气锅炉效率；COP 表示电制冷机转换系数；COP 表示吸收式制冷机转换系数；

具体的，园区内综合能源系统与工业用户电、热、冷功率平衡约束为：具体的，电网及气网交互功率约束为：

上式中： 表示电网购电及气网购气最大功率。