CN105207205A - Distributed energy system energy optimization regulation and control method fusing demand side response - Google Patents

Abstract

Translated fromChinese

本发明公开了一种融合需求侧响应的分布式能源系统能量优化调控方法，包括建立分布式能源系统优化调控模型，通过引入用户制冷设备的响应能力来参与分布式能源系统优化运行；通过考虑系统内建筑的热平衡方程，来定量描述获得用户室内温度与制冷设备出力之间的数学关系。本发明在满足用户用能舒适性的基础上，减少对电池储能系统的依赖性，有效提高分布式能源系统运行的经济性；建立了分布式能源系统优化调控模型，通过引入用户侧制冷设备的响应能力来参与分布式能源系统优化运行，同时通过考虑系统内建筑的热平衡方程，来定量描述获得用户室内温度与制冷设备出力之间的数学关系，保证用户的舒适性。结果表明引入用户需求响应可有效节约用能成本。

The invention discloses an energy optimization control method for a distributed energy system that integrates demand-side response, including establishing an optimal control model for a distributed energy system, and participating in the optimal operation of the distributed energy system by introducing the response capability of user refrigeration equipment; The heat balance equation of the internal building is used to quantitatively describe the mathematical relationship between the user's indoor temperature and the output of the refrigeration equipment. On the basis of satisfying the user's energy comfort, the present invention reduces the dependence on the battery energy storage system and effectively improves the economy of the operation of the distributed energy system; establishes an optimal control model for the distributed energy system, and introduces user-side refrigeration equipment To participate in the optimal operation of the distributed energy system, and by considering the heat balance equation of the building in the system, quantitatively describe the mathematical relationship between the indoor temperature of the user and the output of the refrigeration equipment to ensure the comfort of the user. The results show that the introduction of user demand response can effectively save energy costs.

Description

Translated fromChinese

一种融合需求侧响应的分布式能源系统能量优化调控方法A distributed energy system energy optimization control method integrating demand side response

技术领域technical field

本发明属于分布式能源技术领域，尤其涉及一种融合需求侧响应的分布式能源系统能量优化调控方法。The invention belongs to the technical field of distributed energy sources, and in particular relates to an energy optimization control method of a distributed energy system integrating demand side response.

背景技术Background technique

分布式能源系统是指分布在用户端的能源综合利用系统，可通过对风/光/储等分布式能源的综合管理和优化调度实现用户冷/热/电多种能源的梯级利用，以有效提高综合能源利用效率，降低用户用能成本。然而，目前分布式能源系统中高比例间歇性分布式电源的消纳需要配置大容量的电池储能系统，以保障用户供电质量与可靠性。当前储能设备成本高，大容量投入在经济上和技术上都比较难以实现，给分布式能源系统的安全高效运行带来了一定困难。如何在保障运行经济性的同时实现对高比例间歇性分布式电源的充分消纳，是当前分布式能源系统的推广与应用面临的主要技术瓶颈。针对此问题，部分学者开展了深入的研究：由于目前技术的欠缺，电池储能系统的成本依旧居高不下；同时，分布式电源的功率输出具有随机不确定性，由此会导致电池储能系统的频繁冲放电，加快了蓄电池的老化过程。The distributed energy system refers to the comprehensive energy utilization system distributed at the user end. Through the comprehensive management and optimal scheduling of distributed energy such as wind/light/storage, the cascade utilization of various energy sources such as cold/heat/electricity can be realized to effectively improve Comprehensive energy utilization efficiency reduces energy consumption costs for users. However, the consumption of a high proportion of intermittent distributed power in the current distributed energy system requires the configuration of a large-capacity battery energy storage system to ensure the quality and reliability of power supply for users. At present, the cost of energy storage equipment is high, and large-capacity investment is relatively difficult to achieve economically and technically, which brings certain difficulties to the safe and efficient operation of distributed energy systems. How to fully accommodate the high-proportion intermittent distributed power while ensuring the economical operation is the main technical bottleneck facing the promotion and application of the current distributed energy system. In response to this problem, some scholars have carried out in-depth research: due to the lack of current technology, the cost of battery energy storage systems is still high; at the same time, the power output of distributed power sources has random uncertainty, which will lead to battery energy storage. The frequent charging and discharging of the system accelerates the aging process of the battery.

当前电池储能设备成本高，大容量投入在经济上和技术上都比较难以实现，给分布式能源系统的安全高效运行带来了一定困难。At present, the cost of battery energy storage equipment is high, and large-capacity investment is relatively difficult to achieve economically and technically, which brings certain difficulties to the safe and efficient operation of distributed energy systems.

发明内容Contents of the invention

本发明的目的在于提供一种融合需求侧响应的分布式能源系统能量优化调控方法，旨在解决当前电池储能设备成本高，大容量投入在经济上和技术上都比较难以实现，给分布式能源系统的安全高效运行带来困难的问题。The purpose of the present invention is to provide an energy optimization and control method for distributed energy systems that integrates demand-side response, aiming at solving the problem of high cost of current battery energy storage equipment, large-capacity investment is difficult to achieve economically and technically, and provides distributed The safe and efficient operation of energy systems poses difficult issues.

本发明是这样实现的，一种融合需求侧响应的分布式能源系统能量优化调控方法，所述融合需求侧响应的分布式能源系统能量优化调控方法建立分布式能源系统优化调控模型，通过引入用户制冷设备的响应能力来参与分布式能源系统优化运行；同时通过考虑系统内建筑的热平衡方程，来定量描述获得用户室内温度与制冷设备出力之间的数学关系。The present invention is achieved in this way, an energy optimization control method for a distributed energy system that integrates demand-side response. The energy optimization control method for a distributed energy system that integrates demand-side response establishes an optimal control model for a distributed energy system. The responsiveness of the refrigeration equipment is used to participate in the optimal operation of the distributed energy system; at the same time, by considering the heat balance equation of the building in the system, it is quantitatively described to obtain the mathematical relationship between the user's indoor temperature and the output of the refrigeration equipment.

进一步，所述融合需求侧响应的分布式能源系统能量优化调控方法具体包括：Further, the energy optimization control method of the distributed energy system that integrates the demand side response specifically includes:

步骤一，基于分布式能源系统所接入的主网电力市场环境、分布式能源系统自身的网络构成、需求侧响应资源、分布式电源、储能设备的特性，建立分布式能源系统的能量管理优化调度模型、约束及目标，获得相关参数信息；Step 1. Based on the main grid power market environment connected to the distributed energy system, the network composition of the distributed energy system itself, demand-side response resources, distributed power sources, and the characteristics of energy storage equipment, the energy management of the distributed energy system is established. Optimize the scheduling model, constraints and objectives, and obtain relevant parameter information;

步骤二，对于分布式能源系统内的分布式电源和负荷，所处的外界环境进行能量管理调度周期内的预测；Step 2. For the distributed power sources and loads in the distributed energy system, the external environment where they are located is predicted within the energy management scheduling cycle;

步骤三，基于所建立目标分布式能源系统的能量管理优化调度模型，进行能量管理优化调度求解；Step 3, based on the established energy management optimal scheduling model of the target distributed energy system, the energy management optimal scheduling solution is performed;

步骤四，根据优化模型，通过内点法获得运算周期内分布式能源系统能量管理每个阶段的运行方案一天之内的出力调度方案；Step 4, according to the optimization model, obtain the output scheduling plan within one day of the operation plan of each stage of the energy management of the distributed energy system within the operation cycle through the interior point method;

步骤五基于优化所得调度方案，利用能量管理系统向各受控对象下发控制命令。Step 5: Based on the optimized scheduling scheme, the energy management system is used to issue control commands to each controlled object.

进一步，所述构建优化调控模型包括用户舒适度与制冷设备出力之间数学模型，具体构建方法如下：P_ex为分布式能源系统与外电网的交换功率，当为正时，分布式能源系统从电网获得电量，为负时，分布式能源系统向外电网售电；P_wt、P_pv分别表示风力发电和光伏发电设备的出力；P_el为分布式能源系统中的电负荷；P_bt表示分布式能源系统中蓄电池的出力，当其为正时，蓄电池充电，为负时，蓄电池放电；P_ec为电制冷机的输入电功率；Q_ec为电制冷机的输出制冷功率；Q_cl为分布式能源系统中的冷负荷，全部为夏季降低室温所用；Further, the construction optimization control model includes a mathematical model between user comfort and refrigeration equipment output. The specific construction method is as follows: P_ex is the exchange power between the distributed energy system and the external power grid. When it is positive, the distributed energy system starts from When the electricity obtained by the grid is negative, the distributed energy system sells electricity to the external grid; P_wt and P_pv respectively represent the output of wind power and photovoltaic power generation equipment; P_el is the electric load in the distributed energy system; P_bt represents the distribution The output of the battery in the type energy system, when it is positive, the battery is charged, and when it is negative, the battery is discharged; P_ec is the input electric power of the electric refrigerator; Q_ec is the output cooling power of the electric refrigerator; Q_cl is the distributed The cooling load in the energy system is all used to lower the room temperature in summer;

步骤一，获得用户室内温度与制冷设备出力之间的数学关系，依据能量守恒的原则得到建筑的热平衡方程，如式(1)所示：Step 1. Obtain the mathematical relationship between the user’s indoor temperature and the output of the refrigeration equipment, and obtain the heat balance equation of the building according to the principle of energy conservation, as shown in formula (1):

$\mathrm{<span><span>\&Delta;</span>\&Delta;</span>}\mathrm{<span><span>Q</span>Q</span>}<span><span>=</span>=</span>\mathrm{<span><span>\&rho;</span>\&rho;</span>}<span><span>\&CenterDot;</span>\&CenterDot;</span>\mathrm{<span><span>C</span>C</span>}<span><span>\&CenterDot;</span>\&CenterDot;</span>\mathrm{<span><span>V</span>V</span>}<span><span>\&CenterDot;</span>\&CenterDot;</span>\frac{{\mathrm{<span><span>dT</span>dT</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}}{\mathrm{<span><span>d</span>d</span>}\mathrm{<span><span>\&tau;</span>\&tau;</span>}}<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>1</span>1</span>}<span><span>)</span>)</span>$

其中dT_in/dτ为室内温度的变化率；ρ·V为室内空气的质量；C为比热容；ΔQ为室内热量的变化量；Among them, dT_in /dτ is the rate of change of indoor temperature; ρ V is the quality of indoor air; C is the specific heat capacity; ΔQ is the change of indoor heat;

步骤二，式(1)进一步转化为式(2)：Step 2, formula (1) is further transformed into formula (2):

${\mathrm{<span><span>k</span>k</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>a</span>a</span>}\mathrm{<span><span>l</span>l</span>}\mathrm{<span><span>l</span>l</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>a</span>a</span>}\mathrm{<span><span>l</span>l</span>}\mathrm{<span><span>l</span>l</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span><span><span>(</span>(</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>u</span>u</span>}\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>)</span>)</span><span><span>+</span>+</span>{\mathrm{<span><span>k</span>k</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span><span><span>(</span>(</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>u</span>u</span>}\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>)</span>)</span><span><span>+</span>+</span>\stackrel{<span><span>\&CenterDot;</span>\&CenterDot;</span>}{\mathrm{<span><span>I</span>I</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>S</span>S</span>}\mathrm{<span><span>C</span>C</span>}<span><span>+</span>+</span>\mathrm{<span><span>Q</span>Q</span>}<span><span>-</span>-</span>{\mathrm{<span><span>Q</span>Q</span>}}_{\mathrm{<span><span>c</span>c</span>}\mathrm{<span><span>l</span>l</span>}}<span><span>=</span>=</span>\mathrm{<span><span>\&rho;</span>\&rho;</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>C</span>C</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>V</span>V</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\frac{{\mathrm{<span><span>dT</span>dT</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}}{\mathrm{<span><span>d</span>d</span>}\mathrm{<span><span>\&tau;</span>\&tau;</span>}}<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>2</span>2</span>}<span><span>)</span>)</span>$

其中，k_wall·F_wall·(T_out-T_in)为建筑外墙与室外传递的热量；k_wall为建筑外墙的传热系数；F_wall为建筑外墙面积；(T_out-T_in)为室内外温度差；k_win·F_win·(T_out-T_in)为建筑外窗与室外传递的热量；k_win为建筑外窗的传热系数；F_win为建筑外窗的面积；表示太阳热辐射传递的热量；为太阳辐射功率，表示与光照垂直时每平方米每秒接受的热量；SC为遮阳系数，与是否有遮阳板、玻璃材质有关；Q为室内热源的发热功率；Q_cl为制冷设备的制冷功率；Among them, k_wall F_wall (T_out -T_in ) is the heat transfer between the building exterior wall and the outside; k_wall is the heat transfer coefficient of the building exterior wall; F_wall is the area of the building exterior wall; (T_out -T_in ) is the indoor and outdoor temperature difference; k_win F_win (T_out -T_in ) is the heat transfer between the building’s exterior windows and the outside; k_win is the heat transfer coefficient of the building’s exterior windows; F_win is the area of the building’s exterior windows; Indicates the heat transferred by solar thermal radiation; is the solar radiation power, indicating the heat received per square meter per second when it is perpendicular to the light; SC is the shading coefficient, which is related to whether there is a sun visor or glass material; Q is the heating power of the indoor heat source; Q_cl is the cooling power of the refrigeration equipment ;

步骤三，约束条件与优化决策变量：Step 3, constraints and optimization decision variables:

电母线的能量平衡方程如下：The energy balance equation of the electric bus is as follows:

P_ex+P_WT+P_PV＝P_el+P_bt+P_ec(3)P_ex +P_WT +P_PV ＝P_el +P_bt +P_ec (3)

冷母线能量平衡方程：Cold bus energy balance equation:

Q_ec＝Q_cl(4)Q_ec = Q_cl (4)

电制冷设备能量转换方程，其中EER为电制冷设备的制冷能效比：Energy conversion equation of electric refrigeration equipment, where EER is the cooling energy efficiency ratio of electric refrigeration equipment:

Q_ec＝P_ec·EER(5)Q_ec =P_ec ·EER(5)

式子(4)和(5)合为式(6)：Formulas (4) and (5) are combined into formula (6):

Q_cl＝P_ec·EER(6)Q_cl ＝P_ec ·EER(6)

式(1)和(6)即为该分布式能源系统利用母线式结构所得到的能量平衡方程；Equations (1) and (6) are the energy balance equations obtained by using the busbar structure of the distributed energy system;

步骤四，将式(6)带入热平衡方程(2)，得到以电制冷设备功率表达的热平衡方程：Step 4, bring equation (6) into heat balance equation (2), and obtain the heat balance equation expressed in terms of the power of the electric refrigeration equipment:

$\begin{array}{c}{\mathrm{<span><span>k</span>k</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>a</span>a</span>}\mathrm{<span><span>l</span>l</span>}\mathrm{<span><span>l</span>l</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>a</span>a</span>}\mathrm{<span><span>l</span>l</span>}\mathrm{<span><span>l</span>l</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span><span><span>(</span>(</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>u</span>u</span>}\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>)</span>)</span><span><span>+</span>+</span>{\mathrm{<span><span>k</span>k</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span><span><span>(</span>(</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>u</span>u</span>}\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>)</span>)</span><span><span>+</span>+</span>\stackrel{<span><span>\&CenterDot;</span>\&CenterDot;</span>}{\mathrm{<span><span>I</span>I</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>\mathrm{<span><span>S</span>S</span>}\mathrm{<span><span>C</span>C</span>}<span><span>+</span>+</span>\mathrm{<span><span>Q</span>Q</span>}\\ <span><span>-</span>-</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>\mathrm{<span><span>E</span>E.</span>}\mathrm{<span><span>E</span>E.</span>}\mathrm{<span><span>R</span>R</span>}<span><span>=</span>=</span>\mathrm{<span><span>\&rho;</span>\&rho;</span>}<span><span>\&CenterDot;</span>\&CenterDot;</span>\mathrm{<span><span>C</span>C</span>}<span><span>\&CenterDot;</span>\&CenterDot;</span>\mathrm{<span><span>V</span>V</span>}<span><span>\&CenterDot;</span>\&CenterDot;</span>\frac{{\mathrm{<span><span>dT</span>dT</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}}{\mathrm{<span><span>d</span>d</span>}\mathrm{<span><span>\&tau;</span>\&tau;</span>}}\end{array}<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>7</span>7</span>}<span><span>)</span>)</span>$

分布式能源系统中各设备的功率上下限约束如式(8)所示：The power upper and lower limits of each device in the distributed energy system are shown in formula (8):

$\left\{\begin{array}{c}\underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}}}<span><span><</span><</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}}}\\ \underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}<span><span><</span><</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}\\ \underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}}}<span><span><</span><</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}}}\end{array}\right.<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>8</span>8</span>}<span><span>)</span>)</span>$

步骤五，蓄电池荷电状态的约束，如下式(9)所示：Step 5, the constraint on the state of charge of the battery, as shown in the following formula (9):

$\underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}<span><span><</span><</span>{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>9</span>9</span>}<span><span>)</span>)</span>$

其中，W_bt表示蓄电池某一时刻的电量，例如时刻t时：Among them, W_bt represents the power of the battery at a certain moment, for example, at time t:

$\underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}<span><span><</span><</span>{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>(</span>(</span>\mathrm{<span><span>t</span>t</span>}<span><span>)</span>)</span>}<span><span>=</span>=</span>{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>(</span>(</span>\mathrm{<span><span>0</span>0</span>}<span><span>)</span>)</span>}<span><span>+</span>+</span>{<span><span>\&Sigma;</span>\&Sigma;</span>}_{\mathrm{<span><span>1</span>1</span>}}^{\mathrm{<span><span>t</span>t</span>}}{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>(</span>(</span>\mathrm{<span><span>t</span>t</span>}<span><span>)</span>)</span>}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>10</span>10</span>}<span><span>)</span>)</span>$

其中，W_bt(0)为蓄电池的初始电量；Among them, W_bt(0) is the initial power of the battery;

步骤六，分布式能源系统优化调控目标函数：Step six, the distributed energy system optimizes the control objective function:

$\begin{array}{c}\begin{array}{cc}\mathrm{<span><span>min</span>min</span>}& \mathrm{<span><span>C</span>C</span>}<span><span>=</span>=</span>\underset{\mathrm{<span><span>i</span>i</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}\end{array}<span><span>\&lsqb;</span>\&lsqb;</span><span><span>(</span>(</span>\frac{{\mathrm{<span><span>C</span>C</span>}}_{\mathrm{<span><span>p</span>p</span>}\mathrm{<span><span>h</span>h</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>+</span>+</span>{\mathrm{<span><span>C</span>C</span>}}_{\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>e</span>e</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}}{\mathrm{<span><span>2</span>2</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>+</span>+</span>\frac{{\mathrm{<span><span>C</span>C</span>}}_{\mathrm{<span><span>p</span>p</span>}\mathrm{<span><span>h</span>h</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>C</span>C</span>}}_{\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>e</span>e</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}}{\mathrm{<span><span>2</span>2</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span><span><span>|</span>|</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>|</span>|</span><span><span>)</span>)</span>\\ <span><span>+</span>+</span><span><span>(</span>(</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>W</span>W</span>}\mathrm{<span><span>T</span>T</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>C</span>C</span>}}_{{\mathrm{<span><span>WT</span>WT</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>+</span>+</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>P</span>P</span>}\mathrm{<span><span>V</span>V</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>C</span>C</span>}}_{{\mathrm{<span><span>PV</span>PV</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>+</span>+</span><span><span>|</span>|</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>|</span>|</span><span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>C</span>C</span>}}_{{\mathrm{<span><span>bt</span>bt</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>+</span>+</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>C</span>C</span>}}_{{\mathrm{<span><span>ec</span>ec</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>)</span>)</span>\\ <span><span>+</span>+</span>\mathrm{<span><span>\&gamma;</span>\&gamma;</span>}<span><span>\&CenterDot;</span>\&CenterDot;</span><span><span>|</span>|</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>t</span>t</span>}}<span><span>|</span>|</span><span><span>\&rsqb;</span>\&rsqb;</span>\end{array}<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>13</span>13</span>}<span><span>)</span>)</span>$

式中第一项为分布式能源系统与电网能量交换带来的净支出；C_ph,i为i时刻从电网购电的电价；C_se,i为i时刻向电网售电的电价；The first item in the formula is the net expenditure brought by the energy exchange between the distributed energy system and the grid; C_ph,i is the price of electricity purchased from the grid at time i; C_se,i is the price of electricity sold to the grid at time i;

式中第二项代表分布式能源系统中各设备的使用维护成本。分别代表风机、光伏、蓄电池和电制冷机单位时间段、单位功率的使用维护成本；The second item in the formula Represents the use and maintenance cost of each device in the distributed energy system. Represent the use and maintenance costs per unit time period and unit power of fans, photovoltaics, batteries and electric refrigerators, respectively;

式中第三项r·｜T_in,i-T_set｜为影响用户舒适性而设的罚函数项，γ为罚因子，为用户对舒适性的敏感程度，单位为元/℃；γ根据不同的用户敏感性来选择，称为用户敏感系数；In the formula, the third item r·|T_in,i- T_set ｜is a penalty function item set up to affect user comfort, γ is a penalty factor, and is the user’s sensitivity to comfort, and the unit is yuan/°C; γ is based on Different user sensitivities are selected, called user sensitivity coefficients;

步骤七，约束和目标函数共同组成了一个混合整数线性规划模型，以每15分钟为一个时间点，全天共96个时刻，分布式能源系统的融合需求侧响应优化调控模型如(14)所示：Step 7: Constraints and objective functions together form a mixed integer linear programming model, with every 15 minutes as a time point, a total of 96 moments throughout the day, the integrated demand response optimization control model of the distributed energy system is shown in (14) Show:

$\left\{\begin{array}{cc}\mathrm{<span><span>min</span>min</span>}& \mathrm{<span><span>C</span>C</span>}<span><span>=</span>=</span>\underset{\mathrm{<span><span>i</span>i</span>}<span><span>=</span>=</span>\mathrm{<span><span>1</span>1</span>}}{\overset{\mathrm{<span><span>96</span>96</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}}<span><span>\&lsqb;</span>\&lsqb;</span><span><span>(</span>(</span>\frac{{\mathrm{<span><span>C</span>C</span>}}_{\mathrm{<span><span>p</span>p</span>}\mathrm{<span><span>h</span>h</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>+</span>+</span>{\mathrm{<span><span>C</span>C</span>}}_{\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>e</span>e</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}}{\mathrm{<span><span>2</span>2</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}\frac{{\mathrm{<span><span>C</span>C</span>}}_{\mathrm{<span><span>p</span>p</span>}\mathrm{<span><span>h</span>h</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>C</span>C</span>}}_{\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>e</span>e</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}}{\mathrm{<span><span>2</span>2</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span><span><span>|</span>|</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>|</span>|</span><span><span>)</span>)</span>\\ & <span><span>+</span>+</span><span><span>(</span>(</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>W</span>W</span>}\mathrm{<span><span>T</span>T</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>C</span>C</span>}}_{{\mathrm{<span><span>WT</span>WT</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>+</span>+</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>P</span>P</span>}\mathrm{<span><span>V</span>V</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>C</span>C</span>}}_{{\mathrm{<span><span>PV</span>PV</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>+</span>+</span><span><span>|</span>|</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>|</span>|</span><span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>C</span>C</span>}}_{{\mathrm{<span><span>bt</span>bt</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>+</span>+</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>C</span>C</span>}}_{{\mathrm{<span><span>ec</span>ec</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>)</span>)</span>\\ & <span><span>+</span>+</span>\mathrm{<span><span>\&gamma;</span>\&gamma;</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span><span><span>|</span>|</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>t</span>t</span>}}<span><span>|</span>|</span><span><span>\&rsqb;</span>\&rsqb;</span>\\ \mathrm{<span><span>s</span>the\; s</span>}<span><span>.</span>.</span>\mathrm{<span><span>t</span>t</span>}<span><span>.</span>.</span>& {\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>+</span>+</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>W</span>W</span>}\mathrm{<span><span>T</span>T</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>+</span>+</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>P</span>P</span>}\mathrm{<span><span>V</span>V</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>l</span>l</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>=</span>=</span>\mathrm{<span><span>0</span>0</span>}\\ & {\mathrm{<span><span>k</span>k</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>a</span>a</span>}\mathrm{<span><span>l</span>l</span>}\mathrm{<span><span>l</span>l</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>a</span>a</span>}\mathrm{<span><span>l</span>l</span>}\mathrm{<span><span>l</span>l</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span><span><span>(</span>(</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>u</span>u</span>}\mathrm{<span><span>t</span>t</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>)</span>)</span><span><span>+</span>+</span>{\mathrm{<span><span>k</span>k</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span><span><span>(</span>(</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>u</span>u</span>}\mathrm{<span><span>t</span>t</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>)</span>)</span>\\ & <span><span>+</span>+</span>{\stackrel{<span><span>\&CenterDot;</span>\&Center\; Dot;</span>}{\mathrm{<span><span>I</span>I</span>}}}_{\mathrm{<span><span>t</span>t</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>S</span>S</span>}\mathrm{<span><span>C</span>C</span>}<span><span>+</span>+</span>{\mathrm{<span><span>Q</span>Q</span>}}_{\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>E</span>E.</span>}\mathrm{<span><span>E</span>E.</span>}\mathrm{<span><span>R</span>R</span>}<span><span>-</span>-</span>\mathrm{<span><span>\&rho;</span>\&rho;</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>C</span>C</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>V</span>V</span>}<span><span>(</span>(</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}<span><span>+</span>+</span>\mathrm{<span><span>1</span>1</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>)</span>)</span><span><span>=</span>=</span>\mathrm{<span><span>0</span>0</span>}\\ & \underset{\mathrm{<span><span>i</span>i</span>}<span><span>=</span>=</span>\mathrm{<span><span>1</span>1</span>}}{\overset{\mathrm{<span><span>96</span>96</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}}{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>=</span>=</span>\mathrm{<span><span>0</span>0</span>}\\ & \underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}}}<span><span><</span><</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}}}\\ & \underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}<span><span><</span><</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}\\ & \underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}}}<span><span><</span><</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}}}\\ & \underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}<span><span><</span><</span>{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>(</span>(</span>\mathrm{<span><span>t</span>t</span>}<span><span>)</span>)</span>}<span><span>=</span>=</span>{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>(</span>(</span>\mathrm{<span><span>0</span>0</span>}<span><span>)</span>)</span>}<span><span>+</span>+</span>\underset{\mathrm{<span><span>1</span>1</span>}}{\overset{\mathrm{<span><span>t</span>t</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}}{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>(</span>(</span>\mathrm{<span><span>t</span>t</span>}<span><span>)</span>)</span>}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}\\ & \underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span><</span><</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span><</span><</span>\frac{\mathrm{<span><span>1</span>1</span>}}{{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}}}\end{array}\right.<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>14</span>14</span>}<span><span>)</span>)</span>$

其中t＝1,2,…,96。where t=1,2,...,96.

进一步，优化模型中日结束时刻与起始时刻的电池蓄电量相等，即P_bt还需满足下式约束：Furthermore, in the optimization model, the battery storage capacity at the end of the day is equal to that at the beginning of the day, that is, P_bt also needs to satisfy the following constraints:

∑P_bt＝0(11)ΣP_bt =0(11)

最后，需要考虑建筑室内温度上下限约束：Finally, the upper and lower limits of the indoor temperature of the building need to be considered:

$\underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}}<span><span>-</span>-</span><span><span><</span><</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}}<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>12</span>12</span>}<span><span>)</span>)</span>$

该优化模型中P_wt、P_pv、P_el和均为已知预测量，优化过程中的变量包括P_ex、P_bt、P_ec和T_in，看出T_in的值可根据热平衡方程式(2)由上述已知预测量和P_ec求得，因此T_in为非独立变量，独立决策变量只有三个：P_ex、P_bt、P_ec。In this optimization model, P_wt , P_pv , P_el and Both are known predictive quantities. The variables in the optimization process include P_ex , P_bt , P_ec and T_in . It can be seen that the value of T_in can be obtained from the above known predictive quantities and P_ec according to the heat balance equation (2). Therefore, T_in is a dependent variable, and there are only three independent decision variables: P_ex , P_bt , and P_ec .

本发明提供的融合需求侧响应的分布式能源系统能量优化调控方法，将用户侧电冷热转换装置(电制冷机或热泵等设备)纳入分布式能源系统的调控过程之中，以分布式能源系统综合调控成本最小为目标，建立了分布式能源系统调控方法的优化模型，并进行优化求解，可在满足用户用能舒适性的基础上，减少对电池储能系统的依赖性，有效提高分布式能源系统运行的经济性；建立了分布式能源系统优化调控模型，通过引入用户侧制冷设备的响应能力来参与分布式能源系统优化运行，同时通过考虑系统内建筑的热平衡方程，来定量描述获得用户室内温度与制冷设备出力之间的数学关系，保证用户舒适性。结果表明引入用户需求响应可有效节约用能成本。本发明可通过在温度舒适度范围内对楼宇室温进行调节，实现需求侧资源的优化响应管理，从而降低微网的运行成本。同时，本发明通过对需求侧资源的优化响应管理可降低储能系统的容量配置，在一定程度上降低微网的投资成本。另外，由于用户对舒适性的要求不同，因此本发明引入了用户敏感系数，并在微网调度目标中考虑了因影响用户舒适性而加入的惩罚项，使得优化调度能够对用户的不同敏感系数作出相应调整，满足用户个性化用能需求。The distributed energy system energy optimization control method integrated with the demand side response provided by the present invention incorporates the user-side electric cooling and heat conversion device (electric refrigerator or heat pump and other equipment) into the control process of the distributed energy system, and uses the distributed energy With the goal of minimizing the cost of system comprehensive control, an optimization model of the distributed energy system control method is established and optimized to solve the problem, which can reduce the dependence on the battery energy storage system and effectively improve the energy distribution while satisfying the user's energy comfort. The economy of the operation of the distributed energy system; the optimal control model of the distributed energy system is established, and the response ability of the user-side refrigeration equipment is introduced to participate in the optimal operation of the distributed energy system. At the same time, by considering the heat balance equation of the building in the system, quantitative description is obtained. The mathematical relationship between the user's indoor temperature and the output of the refrigeration equipment ensures user comfort. The results show that the introduction of user demand response can effectively save energy costs. The invention can adjust the room temperature of the building within the range of temperature comfort to realize the optimal response management of resources on the demand side, thereby reducing the operating cost of the microgrid. At the same time, the present invention can reduce the capacity configuration of the energy storage system through the optimized response management of the resources on the demand side, and reduce the investment cost of the microgrid to a certain extent. In addition, since users have different requirements for comfort, the present invention introduces user sensitivity coefficients, and considers the penalty item that affects user comfort in the micro-grid scheduling target, so that the optimal scheduling can respond to different sensitivity coefficients of users Make corresponding adjustments to meet the individual energy consumption needs of users.

附图说明Description of drawings

图1是本发明实施例提供的融合需求侧响应的分布式能源系统能量优化调控方法流程图。Fig. 1 is a flowchart of an energy optimization control method for a distributed energy system integrating demand side response provided by an embodiment of the present invention.

图2是本发明实施例提供的分布式能源系统母线式结构示意图。Fig. 2 is a schematic diagram of a busbar structure of a distributed energy system provided by an embodiment of the present invention.

图3是本发明实施例提供的分布式能源系统日用电量曲线示意图。Fig. 3 is a schematic diagram of a daily power consumption curve of a distributed energy system provided by an embodiment of the present invention.

图4是本发明实施例提供的建筑内热源发热曲线示意图。Fig. 4 is a schematic diagram of a heat generation curve of a heat source in a building provided by an embodiment of the present invention.

图5是本发明实施例提供的实时电价曲线示意图。Fig. 5 is a schematic diagram of a real-time electricity price curve provided by an embodiment of the present invention.

图6是本发明实施例提供的环境温度曲线示意图。Fig. 6 is a schematic diagram of an environmental temperature curve provided by an embodiment of the present invention.

图7是本发明实施例提供的太阳辐射强度曲线示意图。Fig. 7 is a schematic diagram of a solar radiation intensity curve provided by an embodiment of the present invention.

图8是本发明实施例提供的风机出力曲线示意图。Fig. 8 is a schematic diagram of a fan output curve provided by an embodiment of the present invention.

图9是本发明实施例提供的光伏出力曲线示意图。Fig. 9 is a schematic diagram of a photovoltaic output curve provided by an embodiment of the present invention.

图10是本发明实施例提供的联络线交换功率曲线示意图。Fig. 10 is a schematic diagram of a tie line exchange power curve provided by an embodiment of the present invention.

图11是本发明实施例提供的蓄电池充放电曲线示意图。Fig. 11 is a schematic diagram of the charging and discharging curve of the storage battery provided by the embodiment of the present invention.

图12是本发明实施例提供的电制冷机出力曲线示意图。Fig. 12 is a schematic diagram of an output curve of an electric refrigerator provided by an embodiment of the present invention.

图13是本发明实施例提供的室内温度变化曲线示意图。Fig. 13 is a schematic diagram of an indoor temperature change curve provided by an embodiment of the present invention.

图14是本发明实施例提供的电制冷机出力对比示意图。Fig. 14 is a schematic diagram of output comparison of electric refrigerators provided by an embodiment of the present invention.

图15是本发明实施例提供的需求侧响应效果示意图。Fig. 15 is a schematic diagram of the demand side response effect provided by the embodiment of the present invention.

图16是本发明实施例提供的实施方式流程图。Fig. 16 is an implementation flow chart provided by an embodiment of the present invention.

具体实施方式Detailed ways

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。In order to make the object, technical solution and advantages of the present invention more clear, the present invention will be further described in detail below in conjunction with the examples. It should be understood that the specific embodiments described here are only used to explain the present invention, not to limit the present invention.

本发明提出一种融合需求侧响应的分布式能源系统能量优化调控策略，建立了分布式能源系统优化调控模型，通过引入用户制冷设备的响应能力来参与分布式能源系统优化运行，同时通过考虑系统内建筑的热平衡方程，来定量描述获得用户室内温度与制冷设备出力之间的数学关系，保证用户舒适性。结果表明引入用户需求响应可有效节约用能成本。本发明还可应用在冷热电联供机组、电制冷机、吸收式制冷机、电储能的楼宇供能系统中，通过在用户舒适度范围内对建筑物室内温度进行调节，实现需求侧资源的优化响应管理，从而降低楼宇在整个调度周期中的用能成本。The present invention proposes a distributed energy system energy optimization control strategy that integrates demand-side response, establishes a distributed energy system optimization control model, and participates in the distributed energy system optimal operation by introducing the response capability of the user's refrigeration equipment, and at the same time by considering the system The heat balance equation of the internal building is used to quantitatively describe the mathematical relationship between the indoor temperature of the user and the output of the refrigeration equipment to ensure the comfort of the user. The results show that the introduction of user demand response can effectively save energy costs. The present invention can also be applied to building energy supply systems of combined cooling, heating and power units, electric refrigerators, absorption refrigerators, and electric energy storage. Optimized response management of resources, thereby reducing the energy cost of buildings in the entire dispatch cycle.

下面结合附图对本发明的应用原理作详细的描述。The application principle of the present invention will be described in detail below in conjunction with the accompanying drawings.

S101：基于分布式能源系统所接入的主网电力市场环境、分布式能源系统自身的网络构成、需求侧响应资源、分布式电源、储能设备的特性，建立分布式能源系统的能量管理优化调度模型、约束及目标，获得相关参数信息；S101: Based on the power market environment of the main grid connected to the distributed energy system, the network composition of the distributed energy system itself, demand-side response resources, distributed power sources, and the characteristics of energy storage equipment, establish energy management optimization for distributed energy systems Scheduling models, constraints and goals, and obtaining relevant parameter information;

S102：对于分布式能源系统内的分布式电源和负荷，所处的外界环境进行能量管理调度周期内的预测量，本发明优化模型中P_wt、P_pv、P_el和均为已知预测量；S102: For the distributed power sources and loads in the distributed energy system, the external environment in which they are located is used to predict the energy management scheduling cycle. In the optimization model of the present invention, P_wt , P_pv , P_el and Both are known predictive quantities;

S103：基于所建立目标分布式能源系统的能量管理优化调度模型，以P_ex、P_bt、P_ec作为决策变量，进行能量管理优化调度求解；S103: Based on the established energy management optimal scheduling model of the target distributed energy system, use P_ex , P_bt , and P_ec as decision variables to solve energy management optimal scheduling;

S104：根据优化模型，通过内点法获得运算周期内分布式能源系统能量管理每个阶段的P_ex、P_bt、P_ec运行方案一天之内的出力调度方案；S104: According to the optimization model, obtain the output scheduling plan within one day of the P_ex , P_bt , and P_ec operation plans of each stage of the energy management of the distributed energy system within the operation cycle through the interior point method;

S105：基于优化所得调度方案，利用能量管理系统向各受控对象下发控制命令。S105: Based on the optimized scheduling scheme, use the energy management system to issue control commands to each controlled object.

下面结合附图16及具体的实施例对本发明的应用原理作进一步的描述。The application principle of the present invention will be further described below in conjunction with FIG. 16 and specific embodiments.

以图2所示一种典型的分布式能源系统为例，构建其优化调控模型。系统中分布式电源包含风机和光伏发电的分布式能源系统母线式结构如图2所示。P_ex为分布式能源系统与外电网的交换功率，当其为正时，分布式能源系统从电网获得电量，为负时，分布式能源系统向外电网售电；P_wt、P_pv分别表示风力发电和光伏发电设备的出力；P_el为分布式能源系统中的一般电负荷；P_bt表示分布式能源系统中蓄电池的出力，当其为正时，蓄电池充电，为负时，蓄电池放电；P_ec为电制冷机的输入电功率；Q_ec为电制冷机的输出制冷功率；Q_cl为分布式能源系统中的冷负荷，全部为夏季降低室温所用。电制冷机可理解为中央空调主机，为冷水机组，由压缩机、蒸发器、冷凝器、膨胀阀等部分构成，可由分布式能源系统能量管理系统直接控制其功率大小连续变化，以实现参与系统优化运行的功能。Taking a typical distributed energy system shown in Figure 2 as an example, its optimal control model is constructed. The busbar structure of the distributed energy system including wind turbines and photovoltaic power generation in the system is shown in Figure 2. P_ex is the exchange power between the distributed energy system and the external grid. When it is positive, the distributed energy system obtains electricity from the grid; when it is negative, the distributed energy system sells electricity to the external grid; P_wt and P_pv represent The output of wind power and photovoltaic power generation equipment; P_el is the general electrical load in the distributed energy system; P_bt represents the output of the battery in the distributed energy system, when it is positive, the battery is charged, and when it is negative, the battery is discharged; P_ec is the input electric power of the electric refrigerator; Q_ec is the output cooling power of the electric refrigerator; Q_cl is the cooling load in the distributed energy system, all of which are used to lower the room temperature in summer. The electric refrigerator can be understood as the main engine of the central air conditioner, which is a water chiller, which is composed of a compressor, an evaporator, a condenser, an expansion valve, etc., and its power can be directly controlled by the energy management system of the distributed energy system to continuously change, so as to realize the participation in the system Optimized function.

用户舒适度与制冷设备出力之间数学模型Mathematical model between user comfort and refrigeration equipment output

利用用户制冷设备参与系统优化调控运行，必须建立在保证用户舒适性的基础上。对于用户制冷设备而言，衡量用户舒适性的指标是室内温度。因此，首先要获得用户室内温度与制冷设备出力之间的数学关系。The use of user refrigeration equipment to participate in system optimization and regulation operation must be based on ensuring user comfort. For user refrigeration equipment, the index to measure user comfort is the indoor temperature. Therefore, it is first necessary to obtain the mathematical relationship between the user's indoor temperature and the output of the refrigeration equipment.

要精确模拟建筑内温度的变化情况，必须确定影响温度的各个因素。建筑物墙体的传热过程较为复杂，且要考虑对流等因素的影响。为此在保证其准确性在可接受范围内的情况下，这里假设建筑物墙体稳态传热，且不考虑建筑的围护结构蓄热和对流因素的影响，依据能量守恒的原则得到建筑的热平衡方程，如式(1)所示：To accurately simulate temperature changes in a building, it is necessary to determine the factors that affect temperature. The heat transfer process of building walls is relatively complicated, and the influence of factors such as convection must be considered. Therefore, under the condition of ensuring its accuracy within an acceptable range, it is assumed here that the building wall conducts heat in a steady state, and regardless of the influence of the heat storage and convection factors of the building envelope, the building is obtained according to the principle of energy conservation. The heat balance equation, as shown in formula (1):

$\mathrm{<span><span>\&Delta;</span>\&Delta;</span>}\mathrm{<span><span>Q</span>Q</span>}<span><span>=</span>=</span>\mathrm{<span><span>\&rho;</span>\&rho;</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>C</span>C</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>V</span>V</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\frac{{\mathrm{<span><span>dT</span>dT</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}}{\mathrm{<span><span>d</span>d</span>}\mathrm{<span><span>\&tau;</span>\&tau;</span>}}<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>1</span>1</span>}<span><span>)</span>)</span>$

其中dT_in/dτ为室内温度的变化率；ρ·V为室内空气的质量；C为比热容；ΔQ为室内热量的变化量。Among them, dT_in /dτ is the rate of change of the indoor temperature; ρ·V is the quality of the indoor air; C is the specific heat capacity; ΔQ is the variation of the indoor heat.

影响建筑内部热量的主要因素有：室内外温差造成的冷/热耗散，太阳热辐射，建筑内人体及设备发热，及制冷/热设备的效果。以夏季制冷为例，式(1)可进一步转化为式(2)：The main factors affecting the internal heat of a building are: cold/heat dissipation caused by the temperature difference between indoor and outdoor, solar heat radiation, heating of human body and equipment in the building, and the effect of cooling/heating equipment. Taking summer cooling as an example, formula (1) can be further transformed into formula (2):

${\mathrm{<span><span>k</span>k</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>a</span>a</span>}\mathrm{<span><span>l</span>l</span>}\mathrm{<span><span>l</span>l</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>a</span>a</span>}\mathrm{<span><span>l</span>l</span>}\mathrm{<span><span>l</span>l</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span><span><span>(</span>(</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>u</span>u</span>}\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>)</span>)</span><span><span>+</span>+</span>{\mathrm{<span><span>k</span>k</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span><span><span>(</span>(</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>u</span>u</span>}\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>)</span>)</span><span><span>+</span>+</span>\stackrel{<span><span>\&CenterDot;</span>\&Center\; Dot;</span>}{\mathrm{<span><span>I</span>I</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>S</span>S</span>}\mathrm{<span><span>C</span>C</span>}<span><span>+</span>+</span>\mathrm{<span><span>Q</span>Q</span>}<span><span>-</span>-</span>{\mathrm{<span><span>Q</span>Q</span>}}_{\mathrm{<span><span>c</span>c</span>}\mathrm{<span><span>l</span>l</span>}}<span><span>=</span>=</span>\mathrm{<span><span>\&rho;</span>\&rho;</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>C</span>C</span>}<span><span>\&CenterDot;</span>\&CenterDot;</span>\mathrm{<span><span>V</span>V</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\frac{{\mathrm{<span><span>dT</span>dT</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}}{\mathrm{<span><span>d</span>d</span>}\mathrm{<span><span>\&tau;</span>\&tau;</span>}}<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>2</span>2</span>}<span><span>)</span>)</span>$

其中，k_wall·F_wall·(T_out-T_in)为建筑外墙与室外传递的热量；k_wall为建筑外墙的传热系数；F_wall为建筑外墙面积；(T_out-T_in)为室内外温度差；k_win·F_win·（T_out-T_in）为建筑外窗与室外传递的热量；k_win为建筑外窗的传热系数；F_win为建筑外窗的面积；表示太阳热辐射传递的热量；为太阳辐射功率，表示与光照垂直时每平方米每秒接受的热量；SC为遮阳系数，与是否有遮阳板、玻璃材质等有关；Q为室内热源的发热功率，如人体及用电设备的发热；Q_cl为制冷设备的制冷功率。Among them, k_wall F_wall (T_out -T_in ) is the heat transfer between the building exterior wall and the outside; k_wall is the heat transfer coefficient of the building exterior wall; F_wall is the area of the building exterior wall; (T_out -T_in ) is the indoor and outdoor temperature difference; k_win F_win (T_out -T_in ) is the heat transfer between the building’s exterior windows and the outside; k_win is the heat transfer coefficient of the building’s exterior windows; F_win is the area of the building’s exterior windows; Indicates the heat transferred by solar thermal radiation; is the solar radiation power, which means the heat received per square meter per second when it is perpendicular to the light; SC is the shading coefficient, which is related to whether there are sun visors, glass materials, etc.; Q is the heating power of indoor heat sources, such as the human body and electrical equipment Heat generation; Q_cl is the cooling power of the refrigeration equipment.

约束条件与优化决策变量Constraints and Optimization Decision Variables

电母线的能量平衡方程如下：The energy balance equation of the electric bus is as follows:

P_ex+P_WT+P_PV＝P_el+P_bt+P_ec(3)P_ex +P_WT +P_PV ＝P_el +P_bt +P_ec (3)

冷母线能量平衡方程：Cold bus energy balance equation:

Q_ec＝Q_cl(4)Q_ec = Q_cl (4)

电制冷设备能量转换方程，其中EER为电制冷设备的制冷能效比：Energy conversion equation of electric refrigeration equipment, where EER is the cooling energy efficiency ratio of electric refrigeration equipment:

Q_ec＝P_ec·EER(5)Q_ec =P_ec ·EER(5)

式子(4)和(5)可合为式(6)：Formulas (4) and (5) can be combined into formula (6):

Q_cl＝P_ec·EER(6)Q_cl ＝P_ec ·EER(6)

式(1)和(6)即为该分布式能源系统利用母线式结构所得到的能量平衡方程。Equations (1) and (6) are the energy balance equations obtained by using the busbar structure of the distributed energy system.

将式(6)带入热平衡方程(2)，得到以电制冷设备功率表达的热平衡方程：Put equation (6) into heat balance equation (2), and get the heat balance equation expressed in terms of power of electric refrigeration equipment:

$\begin{array}{c}{\mathrm{<span><span>k</span>k</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>a</span>a</span>}\mathrm{<span><span>l</span>l</span>}\mathrm{<span><span>l</span>l</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>a</span>a</span>}\mathrm{<span><span>l</span>l</span>}\mathrm{<span><span>l</span>l</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span><span><span>(</span>(</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>u</span>u</span>}\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>)</span>)</span><span><span>+</span>+</span>{\mathrm{<span><span>k</span>k</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span><span><span>(</span>(</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>u</span>u</span>}\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>)</span>)</span><span><span>+</span>+</span>\stackrel{<span><span>\&CenterDot;</span>\&Center\; Dot;</span>}{\mathrm{<span><span>I</span>I</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>S</span>S</span>}\mathrm{<span><span>C</span>C</span>}<span><span>+</span>+</span>\mathrm{<span><span>Q</span>Q</span>}\\ <span><span>-</span>-</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>E</span>E.</span>}\mathrm{<span><span>E</span>E.</span>}\mathrm{<span><span>R</span>R</span>}<span><span>=</span>=</span>\mathrm{<span><span>\&rho;</span>\&rho;</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>C</span>C</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>V</span>V</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\frac{{\mathrm{<span><span>dT</span>dT</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}}{\mathrm{<span><span>d</span>d</span>}\mathrm{<span><span>\&tau;</span>\&tau;</span>}}\end{array}<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>7</span>7</span>}<span><span>)</span>)</span>$

分布式能源系统中各设备的功率上下限约束如式(8)所示：The power upper and lower limits of each device in the distributed energy system are shown in formula (8):

$\left\{\begin{array}{c}\underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}}}<span><span><</span><</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}}}\\ \underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}<span><span><</span><</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}\\ \underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}}}<span><span><</span><</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}}}\end{array}\right.<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>8</span>8</span>}<span><span>)</span>)</span>$

此外，蓄电池在工作过程中除了需要考虑其最大充放电功率的约束外，还需要考虑蓄电池荷电状态的约束，如下式(9)所示：In addition, in the working process of the battery, in addition to the constraints of its maximum charging and discharging power, it is also necessary to consider the constraints of the state of charge of the battery, as shown in the following equation (9):

$\underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}<span><span><</span><</span>{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>9</span>9</span>}<span><span>)</span>)</span>$

其中，W_bt表示蓄电池某一时刻的电量，例如时刻t时，Among them, W_bt represents the battery power at a certain moment, for example, at time t,

$\underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}<span><span><</span><</span>{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>(</span>(</span>\mathrm{<span><span>t</span>t</span>}<span><span>)</span>)</span>}<span><span>=</span>=</span>{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>(</span>(</span>\mathrm{<span><span>0</span>0</span>}<span><span>)</span>)</span>}<span><span>+</span>+</span>{<span><span>\&Sigma;</span>\&Sigma;</span>}_{\mathrm{<span><span>1</span>1</span>}}^{\mathrm{<span><span>t</span>t</span>}}{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>(</span>(</span>\mathrm{<span><span>t</span>t</span>}<span><span>)</span>)</span>}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>10</span>10</span>}<span><span>)</span>)</span>$

其中，W_bt(0)为蓄电池的初始电量。Among them, W_bt(0) is the initial power of the battery.

由于该优化模型是针对全天的情况进行计算的，为保证连续性，需要另该日结束时刻与起始时刻的电池蓄电量相等，即P_bt还需满足下式约束：Since the optimization model is calculated for the whole day, in order to ensure continuity, the battery storage capacity at the end of the day is equal to the start of the day, that is, P_bt also needs to satisfy the following constraints:

∑P_bt＝0(11)ΣP_bt =0(11)

最后，需要考虑建筑室内温度上下限约束：Finally, the upper and lower limits of the indoor temperature of the building need to be considered:

$\underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}}<span><span>-</span>-</span><span><span><</span><</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}}<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>12</span>12</span>}<span><span>)</span>)</span>$

该优化模型中P_wt、P_pv、P_el和均为已知预测量。优化过程中的变量包括P_ex、P_bt、P_ec和T_in，可以看出T_in的值可根据热平衡方程式(2)由上述已知预测量和P_ec求得，因此T_in为非独立变量，独立决策变量只有三个：P_ex、P_bt、P_ec。In this optimization model, P_wt , P_pv , P_el and All are known predictors. The variables in the optimization process include P_ex , P_bt , P_ec and T_in . It can be seen that the value of T_in can be obtained from the above-mentioned known predictive measurements and P_ec according to the heat balance equation (2), so T_in is not independent There are only three independent decision variables: P_ex , P_bt , and P_ec .

分布式能源系统优化调控目标函数Optimal control objective function of distributed energy system

$\begin{array}{c}\begin{array}{cc}\mathrm{<span><span>min</span>min</span>}& \mathrm{<span><span>C</span>C</span>}<span><span>=</span>=</span>\underset{\mathrm{<span><span>i</span>i</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}\end{array}<span><span>\&lsqb;</span>\&lsqb;</span><span><span>(</span>(</span>\frac{{\mathrm{<span><span>C</span>C</span>}}_{\mathrm{<span><span>p</span>p</span>}\mathrm{<span><span>h</span>h</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>+</span>+</span>{\mathrm{<span><span>C</span>C</span>}}_{\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>e</span>e</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}}{\mathrm{<span><span>2</span>2</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>+</span>+</span>\frac{{\mathrm{<span><span>C</span>C</span>}}_{\mathrm{<span><span>p</span>p</span>}\mathrm{<span><span>h</span>h</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>C</span>C</span>}}_{\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>e</span>e</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}}{\mathrm{<span><span>2</span>2</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span><span><span>|</span>|</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>|</span>|</span><span><span>)</span>)</span>\\ <span><span>+</span>+</span><span><span>(</span>(</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>W</span>W</span>}\mathrm{<span><span>T</span>T</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>C</span>C</span>}}_{{\mathrm{<span><span>WT</span>WT</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>+</span>+</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>P</span>P</span>}\mathrm{<span><span>V</span>V</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>C</span>C</span>}}_{{\mathrm{<span><span>PV</span>PV</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>+</span>+</span><span><span>|</span>|</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>|</span>|</span><span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>C</span>C</span>}}_{{\mathrm{<span><span>bt</span>bt</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>+</span>+</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>C</span>C</span>}}_{{\mathrm{<span><span>ec</span>ec</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>)</span>)</span>\\ <span><span>+</span>+</span>\mathrm{<span><span>\&gamma;</span>\&gamma;</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span><span><span>|</span>|</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>t</span>t</span>}}<span><span>|</span>|</span><span><span>\&rsqb;</span>\&rsqb;</span>\end{array}<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>13</span>13</span>}<span><span>)</span>)</span>$

优化调控的主要目标是在保证用户舒适性的基础上，最小化分布式能源系统的综合运行成本。因此其目标函数应有两部分组成，一是经济成本，二是用户因舒适性未被满足而带来的惩罚，其中经济成本又包括购电成本和各设备的使用维护成本。图2所示分布式能源系统融合用户需求侧响应后的优化调控模型的目标函数如式(13)所示。The main goal of optimal regulation is to minimize the comprehensive operating cost of distributed energy systems on the basis of ensuring user comfort. Therefore, its objective function should consist of two parts, one is the economic cost, and the other is the penalty caused by the user’s unsatisfied comfort. The economic cost includes the cost of electricity purchase and the cost of using and maintaining each device. The objective function of the optimal control model of the distributed energy system shown in Figure 2 after integrating user demand side response is shown in formula (13).

式中第一项为分布式能源系统与电网能量交换带来的净支出；C_ph,i为i时刻从电网购电的电价；C_se,i为i时刻向电网售电的电价。The first item in the formula is the net expenditure brought about by the energy exchange between the distributed energy system and the grid; C_ph,i is the price of electricity purchased from the grid at time i; C_se,i is the price of electricity sold to the grid at time i.

式中第二项代表分布式能源系统中各设备的使用维护成本。分别代表风机、光伏、蓄电池和电制冷机单位时间段、单位功率的使用维护成本。The second item in the formula Represents the use and maintenance cost of each device in the distributed energy system. Represent the use and maintenance costs per unit time period and unit power of fans, photovoltaics, batteries and electric refrigerators, respectively.

式中第三项r·｜T_in,i-T_set｜为影响用户舒适性而设的罚函数项，γ为罚因子，可理解为用户对舒适性的敏感程度，单位为元/℃。γ可以根据不同的用户敏感性来选择，称其为用户敏感系数。可以看出γ越大，需求侧响应带来的惩罚将会越大，反之需求侧响应带来的惩罚较小。The third item in the formula, r·|T_in,i -T_set |, is a penalty function item that affects user comfort, and γ is a penalty factor, which can be understood as the user's sensitivity to comfort, and the unit is yuan/°C. γ can be selected according to different user sensitivities, which is called user sensitivity coefficient. It can be seen that the larger γ is, the greater the penalty brought by demand-side response will be, and vice versa, the penalty brought by demand-side response will be smaller.

上述一系列约束和目标函数共同组成了一个混合整数线性规划模型(MILP)。以每15分钟为一个时间点，全天共96个时刻，图2所示分布式能源系统的融合需求侧响应优化调控模型如(15)所示。The above series of constraints and the objective function together constitute a mixed integer linear programming model (MILP). Taking every 15 minutes as a time point, a total of 96 moments throughout the day, the integrated demand-side response optimization control model of the distributed energy system shown in Figure 2 is shown in (15).

$\left\{\begin{array}{cc}\mathrm{<span><span>min</span>min</span>}& \mathrm{<span><span>C</span>C</span>}<span><span>=</span>=</span>\underset{\mathrm{<span><span>i</span>i</span>}<span><span>=</span>=</span>\mathrm{<span><span>1</span>1</span>}}{\overset{\mathrm{<span><span>96</span>96</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}}<span><span>\&lsqb;</span>\&lsqb;</span><span><span>(</span>(</span>\frac{{\mathrm{<span><span>C</span>C</span>}}_{\mathrm{<span><span>p</span>p</span>}\mathrm{<span><span>h</span>h</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>+</span>+</span>{\mathrm{<span><span>C</span>C</span>}}_{\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>e</span>e</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}}{\mathrm{<span><span>2</span>2</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}\frac{{\mathrm{<span><span>C</span>C</span>}}_{\mathrm{<span><span>p</span>p</span>}\mathrm{<span><span>h</span>h</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>C</span>C</span>}}_{\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>e</span>e</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}}{\mathrm{<span><span>2</span>2</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span><span><span>|</span>|</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>|</span>|</span><span><span>)</span>)</span>\\ & <span><span>+</span>+</span><span><span>(</span>(</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>W</span>W</span>}\mathrm{<span><span>T</span>T</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>C</span>C</span>}}_{{\mathrm{<span><span>WT</span>WT</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>+</span>+</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>P</span>P</span>}\mathrm{<span><span>V</span>V</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>C</span>C</span>}}_{{\mathrm{<span><span>PV</span>PV</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>+</span>+</span><span><span>|</span>|</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>|</span>|</span><span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>C</span>C</span>}}_{{\mathrm{<span><span>bt</span>bt</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>+</span>+</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>C</span>C</span>}}_{{\mathrm{<span><span>ec</span>ec</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>)</span>)</span>\\ & <span><span>+</span>+</span>\mathrm{<span><span>\&gamma;</span>\&gamma;</span>}<span><span>\&CenterDot;</span>\&CenterDot;</span><span><span>|</span>|</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>i</span>i</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>t</span>t</span>}}<span><span>|</span>|</span><span><span>\&rsqb;</span>\&rsqb;</span>\\ \mathrm{<span><span>s</span>the\; s</span>}<span><span>.</span>.</span>\mathrm{<span><span>t</span>t</span>}<span><span>.</span>.</span>& {\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>+</span>+</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>W</span>W</span>}\mathrm{<span><span>T</span>T</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>+</span>+</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>P</span>P</span>}\mathrm{<span><span>V</span>V</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>l</span>l</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>=</span>=</span>\mathrm{<span><span>0</span>0</span>}\\ & {\mathrm{<span><span>k</span>k</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>a</span>a</span>}\mathrm{<span><span>l</span>l</span>}\mathrm{<span><span>l</span>l</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>a</span>a</span>}\mathrm{<span><span>l</span>l</span>}\mathrm{<span><span>l</span>l</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span><span><span>(</span>(</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>u</span>u</span>}\mathrm{<span><span>t</span>t</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>)</span>)</span><span><span>+</span>+</span>{\mathrm{<span><span>k</span>k</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span><span><span>(</span>(</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>o</span>o</span>}\mathrm{<span><span>u</span>u</span>}\mathrm{<span><span>t</span>t</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>)</span>)</span>\\ & <span><span>+</span>+</span>{\stackrel{<span><span>\&CenterDot;</span>\&CenterDot;</span>}{\mathrm{<span><span>I</span>I</span>}}}_{\mathrm{<span><span>t</span>t</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>{\mathrm{<span><span>F</span>f</span>}}_{\mathrm{<span><span>w</span>w</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>\&CenterDot;</span>\&CenterDot;</span>\mathrm{<span><span>S</span>S</span>}\mathrm{<span><span>C</span>C</span>}<span><span>+</span>+</span>{\mathrm{<span><span>Q</span>Q</span>}}_{\mathrm{<span><span>t</span>t</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>E</span>E.</span>}\mathrm{<span><span>E</span>E.</span>}\mathrm{<span><span>R</span>R</span>}<span><span>-</span>-</span>\mathrm{<span><span>\&rho;</span>\&rho;</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>C</span>C</span>}<span><span>\&CenterDot;</span>\&Center\; Dot;</span>\mathrm{<span><span>V</span>V</span>}<span><span>(</span>(</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}<span><span>+</span>+</span>\mathrm{<span><span>1</span>1</span>}}<span><span>-</span>-</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>)</span>)</span><span><span>=</span>=</span>\mathrm{<span><span>0</span>0</span>}\\ & \underset{\mathrm{<span><span>i</span>i</span>}<span><span>=</span>=</span>\mathrm{<span><span>1</span>1</span>}}{\overset{\mathrm{<span><span>96</span>96</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}}{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span>=</span>=</span>\mathrm{<span><span>0</span>0</span>}\\ & \underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}}}<span><span><</span><</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>x</span>x</span>}}}\\ & \underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}<span><span><</span><</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}\\ & \underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}}}<span><span><</span><</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>e</span>e</span>}\mathrm{<span><span>c</span>c</span>}}}\\ & \underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}<span><span><</span><</span>{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>(</span>(</span>\mathrm{<span><span>t</span>t</span>}<span><span>)</span>)</span>}<span><span>=</span>=</span>{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>(</span>(</span>\mathrm{<span><span>0</span>0</span>}<span><span>)</span>)</span>}<span><span>+</span>+</span>\underset{\mathrm{<span><span>1</span>1</span>}}{\overset{\mathrm{<span><span>t</span>t</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}}{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}<span><span>(</span>(</span>\mathrm{<span><span>t</span>t</span>}<span><span>)</span>)</span>}<span><span><</span><</span>\stackrel{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>W</span>W</span>}}_{\mathrm{<span><span>b</span>b</span>}\mathrm{<span><span>t</span>t</span>}}}\\ & \underset{<span><span>\&OverBar;</span>\&OverBar;</span>}{{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span><</span><</span>{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}<span><span><</span><</span>\frac{\mathrm{<span><span>1</span>1</span>}}{{\mathrm{<span><span>T</span>T</span>}}_{\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}<span><span>,</span>,</span>\mathrm{<span><span>t</span>t</span>}}}}\end{array}\right.<span><span>-</span>-</span><span><span>-</span>-</span><span><span>-</span>-</span><span><span>(</span>(</span>\mathrm{<span><span>14</span>14</span>}<span><span>)</span>)</span>$

其中t＝1,2,…,96。where t=1,2,...,96.

通过以下的仿真对本发明的应用效果作进一步的说明。The application effect of the present invention will be further described through the following simulation.

通过图2所示的分布式能源系统来验证所提融合需求侧响应的分布式能源系统能量优化调控策略的有效性。系统中参与需求侧响应的含制冷机系统建筑设定为小型独栋办公建筑，长30m，宽20m，层高3m，共三层。建筑外墙采用190mm单排孔砌砖，内外25mm绝热砂浆；窗户为PVC材质塑料窗，玻璃为普通中空玻璃，并设窗户面积占侧面外墙面积的50％。相关参数见表1。The distributed energy system shown in Figure 2 is used to verify the effectiveness of the energy optimization control strategy of the proposed distributed energy system that integrates demand side response. The building with refrigeration system participating in the demand side response in the system is set as a small single-family office building with a length of 30m, a width of 20m, and a floor height of 3m, with three floors in total. The exterior wall of the building adopts 190mm single-row hole bricklaying, and 25mm thermal insulation mortar inside and outside; the windows are made of PVC plastic windows, and the glass is ordinary insulating glass, and the window area accounts for 50% of the side exterior wall area. The relevant parameters are shown in Table 1.

表1参与需求侧响应的含制冷机系统建筑参数信息表Table 1 Information table of building parameters with chiller system participating in demand response

设建筑办公时间为8:00到20:00，该建筑日常规用电(不含制冷用电)曲线如图3所示。建筑内热源发热主要由设备和人体发热两部分组成。设备发热可以近似认为与其用电成正比。办公时间为8:00到20:00，在此时间段内再加入人体发热。可以得到内热源发热曲线，如图4所示。电价采用美国纽约州夏季某典型日的电价折算成人民币的电价。如图5所示。上图所示的是从电网购电的价格，售电时，以该价格乘以某一系数为售电价格，取该系数为0.8。以华北地区8月某典型日为例，室外温度变化如图6。Assuming that the office hours of the building are from 8:00 to 20:00, the daily power consumption curve of the building (excluding cooling power consumption) is shown in Figure 3. The heating of heat sources in buildings is mainly composed of equipment and human body heating. The heating of equipment can be approximately considered to be directly proportional to its power consumption. Office hours are from 8:00 to 20:00, and human body heat is added during this time period. The heating curve of the internal heat source can be obtained, as shown in Figure 4. The electricity price adopts the electricity price of a typical day in summer in New York State of the United States and converts it into RMB electricity price. As shown in Figure 5. The figure above shows the price of electricity purchased from the grid. When selling electricity, the price is multiplied by a certain coefficient for the electricity sales price, and the coefficient is taken as 0.8. Taking a typical day in August in North China as an example, the outdoor temperature changes are shown in Figure 6.

从相关文献可以得到我国北方夏季典型日的太阳辐射强度曲线，如图7所示。该值为与太阳直射方向垂直时所接受的太阳辐射强度，考虑到太阳直射方向与建筑外窗的角度关系、部分外窗背阳以及玻璃的遮阳系数等因素，近似取为0.45*F_win*I_t。空气密度ρ和空气比热容C分别取1.2kg/m³和1000j/(kg·℃)。风力发电和太阳能发电出力与也天气情况有关，预测其出力曲线如下图8，9所示。选用蓄电池容量为300kWh，但蓄电池过充过放会对蓄电池寿命造成极大影响，因此设蓄电池蓄电量不超过240kWh，不低于60kWh。最大充放电功率均为10kW。设蓄电池初始时刻电量为100kWh。分布式能源系统与外电网联络点交换功率的限制双向均为400kW。电制冷设备功率上下限为0到120kWh，制冷能效比EER取3。风机、光伏、蓄电池、电制冷设备使用维护成本分别取0.11、0.08、0.02、0.01元/kWh。加入需求侧响应，设用户在工作时间可以接受温度在设定温度的正负2.5℃的范围内波动，用户的设定温度依然是22.5度。用户的敏感系数γ＝0.1。其它条件不变，优化结果如图10～13所示。The solar radiation intensity curve of a typical summer day in northern my country can be obtained from relevant literature, as shown in Figure 7. This value is the received solar radiation intensity when it is perpendicular to the direction of direct sunlight. Considering the angle relationship between the direction of direct sunlight and the exterior windows of the building, some of the exterior windows are backed by the sun, and the shading coefficient of the glass, it is approximately taken as is 0.45*F_win *I_t . Air density ρ and air specific heat capacity C are taken as 1.2kg/m³ and 1000j/(kg·℃) respectively. The output of wind power and solar power is also related to the weather conditions, and the predicted output curves are shown in Figures 8 and 9 below. The selected battery capacity is 300kWh, but overcharge and overdischarge of the battery will have a great impact on the life of the battery, so the storage capacity of the battery is not more than 240kWh, not less than 60kWh. The maximum charging and discharging power is 10kW. Assume that the initial battery capacity is 100kWh. The exchange power limit between the distributed energy system and the contact point of the external grid is 400kW in both directions. The upper and lower limits of the power of electric refrigeration equipment are 0 to 120kWh, and the refrigeration energy efficiency ratio EER is 3. Use and maintenance costs of fans, photovoltaics, batteries, and electric refrigeration equipment Take 0.11, 0.08, 0.02, and 0.01 yuan/kWh respectively. Add demand-side response, assume that the user can accept temperature fluctuations within the range of plus or minus 2.5°C of the set temperature during working hours, and the user's set temperature is still 22.5°C. User's sensitivity coefficient γ=0.1. Other conditions remain unchanged, the optimization results are shown in Figures 10-13.

分布式能源系统与外电网的功率交换和蓄电池的工作情况，有需求侧响应的情况下二者没有明显变化。电制冷设备的工作情况以及室内温度在工作时间(8:00到20:00)有明显不同，出现了明显波动。该情况下，总运行成本为912.9元。The power exchange between the distributed energy system and the external power grid and the working conditions of the storage battery do not change significantly when there is a demand-side response. The working conditions of the electric refrigeration equipment and the indoor temperature are obviously different during the working hours (8:00 to 20:00), and there are obvious fluctuations. In this case, the total operating cost is 912.9 yuan.

首先对比引入需求侧响应先后，电制冷设备出力曲线，如图14所示。可以看出加入需求侧响应后，电制冷设备的出力以无需求侧响应时电制冷机的出力曲线为基准上下波动。其出力高出基准的部分为蓄冷，即“充电”，低于基准的部分为放冷，即“放电”。将制冷设备出力在这两种情况下的差值绘制成曲线的形式将与蓄电池工作情况类似。以电制冷机基准出力减去加入需求侧响应后的出力，得到曲线如图15所示。First, compare the introduction of demand-side response and the output curve of electric refrigeration equipment, as shown in Figure 14. It can be seen that after the demand side response is added, the output of the electric refrigeration equipment fluctuates up and down based on the output curve of the electric refrigerator without the demand side response. The part whose output is higher than the standard is cold storage, that is, "charging", and the part that is lower than the standard is cooling, that is, "discharge". Plotting the difference in cooling plant output under these two conditions will be in the form of a curve similar to battery operation. The base output of the electric refrigerator is subtracted from the output after adding the demand side response, and the obtained curve is shown in Figure 15.

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。The above descriptions are only preferred embodiments of the present invention, and are not intended to limit the present invention. Any modifications, equivalent replacements and improvements made within the spirit and principles of the present invention should be included in the protection of the present invention. within range.

Claims (8)

1. one kind merges the energy-optimised regulate and control method of distributed energy resource system of Demand Side Response, it is characterized in that, the energy-optimised regulate and control method of distributed energy resource system of described fusion Demand Side Response sets up distributed energy resource system Optimum Regulation model, participates in distributed energy resource system optimizing operation by the responding ability introducing user's refrigeration plant; Simultaneously by the equation of heat balance of building in consideration system, quantitative description obtain user indoor temperature and refrigeration plant exert oneself between mathematical relationship.

2. the energy-optimised regulate and control method of distributed energy resource system merging Demand Side Response as claimed in claim 1, it is characterized in that, the energy-optimised regulate and control method of distributed energy resource system of described fusion Demand Side Response specifically comprises:

Step one, based on major network Power Market, the network formation of distributed energy resource system self, the characteristic of Demand Side Response resource, distributed power source, energy storage device that distributed energy resource system accesses, set up the energy management Optimal Operation Model of distributed energy resource system, obtain relevant parameter information;

Step 2, for the distributed power source in distributed energy resource system and load, residing external environment carries out the prediction in energy management dispatching cycle;

Step 3, based on the energy management Optimal Operation Model of set up target distribution formula energy resource system, carries out energy management Optimized Operation and solves;

Step 4, according to Optimized model, the scheduling scheme of exerting oneself within the operating scheme one day being obtained distributed energy resource system energy management each stage in execution cycle by interior point method.

3. the as claimed in claim 2 energy-optimised regulate and control method of distributed energy resource system merging Demand Side Response, is characterized in that, described structure Optimum Regulation model comprise users'comfort and refrigeration plant exert oneself between Mathematical Modeling, concrete construction method is as follows: P_exfor the exchange power of distributed energy resource system and external power grid, when being timing, distributed energy resource system obtains electricity from electrical network, and for time negative, distributed energy resource system is to external power grid sale of electricity; P_wt, P_pvrepresent exerting oneself of wind power generation and photovoltaic power generation equipment respectively; P_elfor the electric loading in distributed energy resource system; P_btrepresent exerting oneself of storage battery in distributed energy resource system, when it is timing, charge in batteries, for time negative, battery discharging; P_ecfor the input electric power of electric refrigerating machine; Q_ecfor the output refrigeration work consumption of electric refrigerating machine; Q_clfor the refrigeration duty in distributed energy resource system, be all that to reduce room temperature summer used;

Step one, obtain user indoor temperature and refrigeration plant exert oneself between mathematical relationship, keep according to energy

Permanent principle obtains the equation of heat balance built, shown in (1):

$\mathrm{\&Delta;}Q=\mathrm{\&rho;}\&CenterDot;C\&CenterDot;V\&CenterDot;\frac{{\mathrm{dT}}_{in}}{d\mathrm{\&tau;}}---\left(1\right)$

Wherein dT_in/ d τ is the rate of change of indoor temperature; ρ V is the quality of room air; C is specific heat capacity; Δ Q is the variable quantity of indoor heat;

Step 2, the cooling/heating power stage of human body and equipment heating and cooling/heating equipment in building, formula (1) is further converted to formula (2):

$k_{wall}\&CenterDot;F_{wall}\&CenterDot;(T_{out}-T_{in})+k_{win}\&CenterDot;F_{win}\&CenterDot;(T_{out}-T_{in})+\stackrel{\&CenterDot;}{I}\&CenterDot;F_{win}\&CenterDot;SC+Q-Q_{cl}=\mathrm{\&rho;}\&CenterDot;C\&CenterDot;V\&CenterDot;\frac{{\mathrm{dT}}_{in}}{d\mathrm{\&tau;}}---\left(2\right)$

Wherein, k_wallf_wall(T_out-T_in) be the heat of external wall and outdoor transmission; k_wallfor the thermal transmission coefficient of external wall; F_wallfor external wall area; (T_out-T_in) for indoor and outdoor temperature poor; k_winf_win(T_out-T_in) be the heat of external window of building and outdoor transmission; k_winfor the thermal transmission coefficient of external window of building; F_winfor the area of external window of building;represent the heat that sun heat radiation transmits;for solar radiation power, the heat of every square metre of acceptance per second when representing vertical with illumination; Whether SC is shading coefficient, and have sunshading board, glass material relevant; Q is the heating power of indoor airflow; Q_clfor the refrigeration work consumption of refrigeration plant;

Step 3, constraints and Optimal Decision-making variable:

The energy-balance equation of goddess of lightning's line is as follows:

P_ex+P_WT+P_PV＝P_el+P_bt+P_ec(3)

Cold bus energy-balance equation:

Q_ec＝Q_cl(4)

Electricity refrigeration plant power conversion equation, wherein EER is the refrigeration efficiency ratio of electric refrigeration plant:

Q_ec＝P_ec·EER(5)

Formula (4) and (5) are combined into formula (6):

Q_cl＝P_ec·EER(6)

Formula (1) and (6) are the energy-balance equation that this distributed energy resource system utilizes bus-type structure to obtain;

Step 4, brings formula (6) into equation of heat balance (2), obtains the equation of heat balance of expressing with electric refrigeration plant power:

$k_{wall}\&CenterDot;F_{wall}\&CenterDot;(T_{out}-T_{in})+k_{win}\&CenterDot;F_{win}\&CenterDot;(T_{out}-T_{in})+\stackrel{\&CenterDot;}{I}\&CenterDot;F_{win}\&CenterDot;SC+Q$

$-P_{ec}\&CenterDot;EER=\mathrm{\&rho;}\&CenterDot;C\&CenterDot;V\&CenterDot;\frac{{\mathrm{dT}}_{in}}{d\mathrm{\&tau;}}---\left(7\right)$

In distributed energy resource system, the power bound constraint of each equipment is such as formula shown in (8):

$\left\{\begin{array}{c}\underset{\&OverBar;}{P_{\mathrm{ex}}}<P_{\mathrm{ex}}<\stackrel{\&OverBar;}{P_{\mathrm{ex}}}\\ \underset{\&OverBar;}{P_{\mathrm{bt}}}<P_{\mathrm{bt}}<\stackrel{\&OverBar;}{P_{\mathrm{bt}}}\\ \underset{\&OverBar;}{P_{ec}}<P_{ec}<\stackrel{\&OverBar;}{P_{ec}}\end{array}\right.---\left(8\right)$

Step 5, the constraint of storage battery charge state, shown in (9):

$\underset{\&OverBar;}{W_{bt}}<W_{bt}<\stackrel{\&OverBar;}{W_{bt}}---\left(9\right)$

Wherein, W_btrepresent the electricity in storage battery a certain moment, such as, during moment t:

$\underset{\&OverBar;}{W_{bt}}<W_{bt\left(t\right)}=W_{bt\left(0\right)}+{\mathrm{\&Sigma;}}_1^t{P}_{bt\left(t\right)}<\stackrel{\&OverBar;}{W_{bt}}---\left(10\right)$

Wherein, W_{bt (0)}for the initial quantity of electricity of storage battery;

Step 6, distributed energy resource system Optimum Regulation target function:

$\begin{array}{c}\min C=\underset{i}{\mathrm{\&Sigma;}}\&lsqb;(\frac{C_{ph,i}+C_{se,i}}{2}\&CenterDot;P_{ex,i}+\frac{C_{ph,i}-C_{se,i}}{2}\&CenterDot;|P_{ex,i}\left|\right)\\ +(P_{WT,i}\&CenterDot;C_{{\mathrm{WT}}_{om}}+P_{PV,i}\&CenterDot;C_{{\mathrm{PV}}_{om}}+|P_{bt,i}|\&CenterDot;C_{{\mathrm{bt}}_{om}}+P_{ec,i}\&CenterDot;C_{{\mathrm{ec}}_{om}})\\ +\mathrm{\&gamma;}\&CenterDot;|T_{in,i}-T_{set}|\&rsqb;\end{array}---\left(13\right)$

Section 1 in formulafor distributed energy resource system and power grid energy exchange the net disbursement brought; C_{ph, i}for the i moment is from the electricity price of electrical network power purchase; C_{se, i}for the i moment is to the electricity price of electrical network sale of electricity;

Section 2 in formularepresent the working service cost of each equipment in distributed energy resource system;represent the working service cost of blower fan, photovoltaic, storage battery and electric refrigerating machine unit interval section, unit power respectively;

Section 3 γ in formula | T_{in, i}-T_set| for affecting the penalty function item that user's comfortableness is established, γ is penalty factor, for user is to the sensitivity of comfortableness, unit be unit/DEG C; γ selects according to different user's sensitiveness, is called user's coefficient of sensitivity;

Step 7, constraint and target function constitute a MILP model jointly, and be a time point with every 15 minutes, whole day is totally 96 moment, and the fusion Demand Side Response Optimum Regulation model of distributed energy resource system is as shown in (14):

Wherein t=1,2 ..., 96.

4. the energy-optimised regulate and control method of distributed energy resource system merging Demand Side Response as claimed in claim 3, it is characterized in that, in Optimized model, the end of day moment is equal with the battery charge capacity of initial time, i.e. P_btalso need to meet following formula constraint:

∑P_bt＝0(11)

Finally, need to consider the constraint of architecture indoor temperature bound:

$\underset{\&OverBar;}{T_{in}}<T_{in}<\stackrel{\&OverBar;}{T_{in}}---\left(12\right)$

P in this Optimized model_wt, P_pv, P_elwithbe known premeasuring, the variable in optimizing process comprises P_ex, P_bt, P_ecand T_in, find out T_invalue according to heat balance equation (2) by above-mentioned known premeasuring and P_ectry to achieve, therefore T_infor non-independent variable, independent decision-making variable only has three: P_ex, P_bt, P_ec.

5. one kind uses the supply of cooling, heating and electrical powers unit of the energy-optimised regulate and control method of the distributed energy resource system merging Demand Side Response described in claim 1-4 any one.

6. one kind uses the electric refrigerating machine of the energy-optimised regulate and control method of the distributed energy resource system merging Demand Side Response described in claim 1-4 any one.

7. one kind uses the Absorption Refrigerator of the energy-optimised regulate and control method of the distributed energy resource system merging Demand Side Response described in claim 1-4 any one.

8. one kind uses the building energy supplying system of the electric energy storage of the energy-optimised regulate and control method of the distributed energy resource system merging Demand Side Response described in claim 1-4 any one.