1.一种计及用户不满意度的综合能源系统优化方法，其特征在于包括以下步骤：(101)获得综合能源系统机组运行参数、状态信息和用户实际负荷需求信息；

(102)根据所述综合能源系统机组运行参数和状态信息，建立综合能源系统机组设备模型；所述综合能源系统机组设备模型包括冷热电联产机组装置模型、热泵装置模型、光伏发电装置模型、电制冷机装置模型、暖通空调系统模型、储能装置模型；所述步骤(102)根据综合能源系统机组运行参数和状态信息，建立综合能源系统机组设备模型，具体包括：所述冷热电联产机组装置模型为：

其中， 分别表示t时刻冷热电联产机组的电输出量、热输出量、冷输出量， 分别表示冷热电联产机组的电转换效率、热转换效率、冷转换效率， 表示t时刻冷热电联产机组的天然气消耗量；

所述热泵装置模型为：

其中， 表示热泵装置t时刻的热输出量，ηHP表示热泵装置的热转换效率， 表示热泵装置在t时刻消耗的电量；

所述光伏发电装置模型为：

其中， 表示光伏发电装置在t时刻的发电量，ηPV为光伏发电装置的发电效率，S为光伏发电装置的面积， 为光伏发电装置单位面积在t时刻的光照强度；

所述电制冷机装置模型为：

其中， 表示电制冷机装置t时刻的冷输出量，ηER表示电制冷机装置的冷转换效率，表示电制冷机装置在t时刻消耗的电量；

所述暖通空调系统模型为：

其中， 表示暖通空调系统t时刻的冷输出量，ηHVAC表示暖通空调系统的冷转换效率， 表示暖通空调系统在t时刻消耗的电量，T0表示当前时刻下室外环境温度值，T表示当前时刻下HVAC系统的温度设定值，b1、b2、b3为HVAC系统的能耗参数；

所述储能装置模型为：

其中，S(t)表示储能装置在t时刻的储能状态，Ω表示储能装置的总容量，σESS表示储能装置的自损耗系数，Pch(t)、Pdis(t)分别表示储能装置在t时刻的充、放能功率，ηch、ηdis分别表示储能装置的充能、放能效率，σ表示储能装置的充能、放能状态变量，σ＝1表示充能，σ＝

0表示放能，即保证储能装置不能同时进行充放能；

(103)根据所述综合能源系统机组设备模型建立综合能源系统多能源流动模型；所述综合能源系统多能源流动模型包括电能流、热能流、冷能流及气能流模型；所述步骤(103)根据所述综合能源系统机组设备模型建立综合能源系统多能源流动模型，具体包括：所述综合能源系统多能源流动模型中，

电能流：

热能流：

冷能流：

气能流：

其中，Pt表示综合能源系统在t时刻从外部电网购买的电量， 表示t时刻电储能系统的电输出量，且有 分别表示电储能装置在t时刻的放电、充电功率，表示能源用户的用电设备i在t时刻的电能消耗量，Ht表示综合能源系统在t时刻从外部热网购买的热量， 表示t时刻热储能系统的热输出量，且有 分别表示热储能装置在t时刻的放热、充热功率， 表示能源用户的用热设备i在t时刻的热能消耗量， 表示t时刻冷储能系统的冷输出量，且有 分别表示热冷能装置在t时刻的放冷、充冷功率， 表示能源用户的用冷设备i在t时刻的冷能消耗量， 表示能源用户的生产车间i在t时刻的环境温度平衡冷能消耗量，Gt表示综合能源系统在t时刻从外部气网购买的天然气量， 表示冷热电联产机组i在t时刻消耗的的天然气量；

(104)根据所述用户实际负荷需求信息，建立用户不满意度模型；所述用户不满意度模型包括能源用户设备生产效率不满意度模型和车间工人热不舒适度模型；所述步骤(104)根据用户实际负荷需求信息，建立用户不满意度模型，具体包括：所述能源用户设备生产效率不满意度模型为：

其中， 表示生产设备i在t时刻的不满意度，θ表示不满意度系数， 表示生产设备i在t时刻的实际负荷功率， 表示生产设备i在t时刻的期望负荷功率；

所述车间工人热不舒适度模型为：

4 2

PPD(Ti)＝100‑95exp(‑0.03353PMV(Ti)‑0.2179PMV(Ti))(19)其中，Mi表示人的代谢速率，Wi表示人体的机械功率， 表示人的衣服表面温度， 表示衣服的表面与裸露体表面积之比， 表示穿着衣服的等效热阻， 表示人体周围空气蒸汽分压，表示室内平均辐射温度， 表示对流换热系数， 表示相对空气流速；

(105)根据所述综合能源系统机组设备模型、用户实际负荷需求信息以及用户不满意度模型，以系统总运行成本最低为目标，建立优化目标指标问题；

(106)根据所述优化目标指标问题对综合能源系统和用户负荷进行优化，获得用户负荷需求、HVAC系统温度设定值和综合能源系统各设备出力。

2.根据权利要求1所述的一种计及用户不满意度的综合能源系统优化方法，其特征在于：所述步骤(105)根据综合能源系统机组设备模型、用户实际负荷需求信息以及用户不满意度模型，以系统总运行成本最低为目标，建立优化目标指标问题，所述系统总运行成本包括用户不满意度成本、综合能源系统能源成本及综合能源系统运行维护成本，具体包括：所述用户不满意度成本函数为：

其中，CDis表示用户的不满意度总费用，CIE和CPPD分别表示用户生产效率不满意度费用和车间工人热不舒适度费用，K1、K2表示不满意度成本转换系数， 和 分别表示用户生产效率不满意度和车间工人热不舒适度；

所述以系统总运行成本最低为目标，建立优化目标指标问题中，综合能源系统日常运行费用：

CDOC＝CEnegy+COM(21)

综合能源系统能源费用：

综合能源系统运行维护费用：

综合能源系统优化问题：

其中，CEnegy、COM分别表示综合能源系统的能源成本和运行维护成本，CE、CH、CG分别表示电能成本、热能成本、天然气成本，Pt、Ht、Ft分别表示综合能源系统在t时刻的电、热、气购买量， 分别表示t时刻的电、热、气单价， 统一表示为设备i在t时刻的运行功率， 表示设备i的单位运行功率的运行维护费用，τ为综合能源系统运行成本和用户不满意度成本之间的折衷因子， 分别表示设备i的最小、最大运行功率，Smin、Smax分别表示储能装置的最小、最大储能状态，SL、ST分别表示一个运行周期内储能装置的初始、最终储能状态，Pch,max、Pdis,max分别表示储能装置的最大充能、放能功率。

3.一种权利要求1或2任一权利要求所述的计及用户不满意度的综合能源系统优化方法的系统，其特征在于包括：数据获取模块(201)，用于获得综合能源系统机组运行参数、状态信息和用户实际负荷需求信息；

综合能源系统各能源设备模块(202)，用于建立综合能源系统机组设备模型；

多能流模型建立模块(203)，用于建立综合能源系统多能源流动模型；

用户不满意度模型建立模块(204)，用于建立用户不满意度模型；

优化问题模块(205)，用于建立优化目标指标问题；

优化模块(206)，用于求解优化问题，获得用户负荷需求、HVAC系统温度设定值和综合能源系统各设备出力。