用于供能系统的控制方法、系统、设备及存储介质与流程

1.本技术涉及控制系统及方法技术领域，尤其涉及一种用于供能系统的控制方法、系统、设备及存储介质。


背景技术：

2.目前，由于储能装置的出现，用户在拥有储能装置后进行存储控制(即改变负荷曲线)，可以在满足其用能需求的情况下实现降低其用能损失量的目的，或者用户通过改变其用能行为直接改变其负荷曲线来实现降低其用能损失量的目的。这样的情况下，能源系统(例如电力系统、热能系统等)中多个用户实施上述行为时，会造成能源浪费的现象，并且不利于整个能源系统的发展和平稳运行。
3.所以，考虑用户可能改变其负荷曲线的情形，如何稳定能源系统中各用户的能耗，令其合理化，且避免能源浪费以及保证能源系统的稳定运行，已成为业界亟待解决的技术问题。


技术实现要素：

4.鉴于以上所述相关技术的缺点，本技术的目的在于提供一种用于供能系统的控制方法、系统、设备及存储介质，用于提供可靠的供能参数以稳定用户的能耗，解决现有技术的问题。
5.为实现上述目的及其他相关目的，本技术的第一方面提供一种用于供能系统的控制方法，所述供能系统向多个用能对象供能；所述供能系统的控制方法包括：根据每一用能对象在每一时段中各时间间隔内的初始供能需求、以及每一时段中各时间间隔内的供能参数，得到预期影响量矩阵；其中，所述预期影响量矩阵中包括每一用能对象在每一所述时段中对于所述供能系统的预期供能参数影响量；在初始影响量矩阵上叠加和初始供能需求相关的噪声矩阵为所述预期影响量矩阵；根据所述预期影响量矩阵、以及矩阵位置集合，在最小化影响量矩阵的秩的条件下构造优化影响量矩阵；其中，所述优化影响量矩阵中包括每一用能对象在每一时段中对于所述供能系统的优化供能参数影响量；所述影响量矩阵和预期影响量矩阵是同型矩阵，且所述矩阵位置集合中包括所述预期影响量矩阵及所述影响量矩阵中部分的相应元素位置；根据所述优化影响量矩阵，确定所述供能系统的实际供能参数，以用于控制所述供能系统按该实际供能参数对各所述用能对象供能。
6.在本技术的第一方面的某些实施方式中，所述根据每一用能对象在每一时段中各时间间隔内的初始供能需求、以及每一时段中各时间间隔内的供能参数，得到预期影响量矩阵，包括：根据每一用能对象在每一时段中各时间间隔内的初始供能需求、以及每一时段中各时间间隔内的供能参数，得到每一用能对象在每一时段中各时间间隔内的预期供能需求；根据每一用能对象在每一时段中各时间间隔内的预期供能需求，确定每一用能对象在每一时段中对于所述供能系统的预期供能参数影响量，以得到所述预期影响量矩阵。
7.在本技术的第一方面的某些实施方式中，每一用能对象在每一时段中对于所述供
能系统的供能参数影响量是根据每一用能对象在每一时段中各时间间隔内的供能需求和每一时段中各时间间隔内的供能参数所确定的。
8.在本技术的第一方面的某些实施方式中，单个用能对象在单个时段中各时间间隔的初始供能需求的总量与预期供能需求的总量相等；单个用能对象在单个时段中各时间间隔内的初始供能需求同各时间间隔内的预期供能需求的差异之和满足预设比例；所述根据每一用能对象在每一时段中各时间间隔内的初始供能需求、以及每一时段中各时间间隔内的供能参数，得到每一用能对象在每一时段中各时间间隔内的预期供能需求，包括：根据单个用能对象在单个时段中各时间间隔内的初始供能需求、以及单个时段各时间间隔内的供能参数，确定使单个用能对象在单个时段中的用能损失量最小所对应的单个用能对象在单个时段中各时间间隔内的预期供能需求。
9.在本技术的第一方面的某些实施方式中，根据所述预期影响量矩阵、以及矩阵位置集合，在最小化影响量矩阵的秩的条件下构造优化影响量矩阵，包括：利用矩阵补全算法，在最小化影响量矩阵的秩的条件下和预设误差约束下构造所述优化影响量矩阵；其中，所述预设误差约束为根据所述矩阵位置集合对影响量矩阵和预期影响量矩阵的差进行投影操作得到的矩阵的范数满足预设误差。
10.在本技术的第一方面的某些实施方式中，所述矩阵补全算法包括：不动点延续算法、加速近端梯度算法、半定规划法、奇异值阈值算法、增广拉格朗日乘子法等。
11.在本技术的第一方面的某些实施方式中，根据所述预期影响量矩阵相对应的奇异值矩阵，确定所述不动点延续算法中影响量矩阵的初始矩阵。
12.在本技术的第一方面的某些实施方式中，每一用能对象的实际供能损失参数为所述优化影响量矩阵中每一用能对象各时段的优化供能参数影响量的均值。
13.为实现上述目的及其他相关目的，本技术的第二方面提供一种用于供能系统的控制系统，所述供能系统向多个用能对象供能；所述供能系统的控制系统包括：预期影响量矩阵确定模块，用于根据每一用能对象在每一时段中各时间间隔内的初始供能需求、以及每一时段中各时间间隔内的供能参数，得到预期影响量矩阵；其中，所述预期影响量矩阵中包括每一用能对象在每一所述时段中对于所述供能系统的预期供能参数影响量；在初始影响量矩阵上叠加和初始供能需求相关的噪声矩阵为所述预期影响量矩阵；优化影响量矩阵确定模块，用于根据所述预期影响量矩阵、以及矩阵位置集合，在最小化影响量矩阵的秩的条件下构造优化影响量矩阵；其中，所述优化影响量矩阵中包括每一用能对象在每一时段中对于所述供能系统的优化供能参数影响量；所述影响量矩阵和预期影响量矩阵是同型矩阵，且所述矩阵位置集合中包括所述预期影响量矩阵及所述影响量矩阵中部分的相应元素位置；供能参数确定模块，用于根据所述优化影响量矩阵，确定所述供能系统的实际供能参数，以用于控制所述供能系统按该实际供能参数对各所述用能对象供能。
14.在本技术的第二方面的某些实施方式中，所述预期影响量矩阵确定模块根据每一用能对象在每一时段中各时间间隔内的初始供能需求、以及每一时段中各时间间隔内的供能参数，得到预期影响量矩阵，包括：根据每一用能对象在每一时段中各时间间隔内的初始供能需求、以及每一时段中各时间间隔内的供能参数，得到每一用能对象在每一时段中各时间间隔内的预期供能需求；根据每一用能对象在每一时段中各时间间隔内的预期供能需求，确定每一用能对象在每一时段中对于所述供能系统的预期供能参数影响量，以得到所
述预期影响量矩阵。
15.在本技术的第二方面的某些实施方式中，每一用能对象在每一时段中对于所述供能系统的供能参数影响量是根据每一用能对象在每一时段中各时间间隔内的供能需求和每一时段中各时间间隔内的供能参数所确定的。
16.在本技术的第二方面的某些实施方式中，单个用能对象在单个时段中各时间间隔的初始供能需求的总量与预期供能需求的总量相等；单个用能对象在单个时段中各时间间隔内的初始供能需求同各时间间隔内的预期供能需求的差异之和满足预设比例；所述根据每一用能对象在每一时段中各时间间隔内的初始供能需求、以及每一时段中各时间间隔内的供能参数，得到每一用能对象在每一时段中各时间间隔内的预期供能需求，包括：根据单个用能对象在单个时段中各时间间隔内的初始供能需求、以及单个时段各时间间隔内的供能参数，确定使单个用能对象在单个时段中的用能损失量最小所对应的单个用能对象在单个时段中各时间间隔内的预期供能需求。
17.在本技术的第二方面的某些实施方式中，所述优化影响量矩阵确定模块根据所述预期影响量矩阵、以及矩阵位置集合，在最小化影响量矩阵的秩的条件下构造优化影响量矩阵，包括：利用矩阵补全算法，在最小化影响量矩阵的秩的条件下和预设误差约束下构造所述优化影响量矩阵；其中，所述预设误差约束为根据所述矩阵位置集合对影响量矩阵和预期影响量矩阵的差进行投影操作得到的矩阵的范数满足预设误差。
18.在本技术的第二方面的某些实施方式中，所述矩阵补全算法包括：不动点延续算法、加速近端梯度算法、半定规划法、奇异值阈值算法、增广拉格朗日乘子法等。
19.在本技术的第二方面的某些实施方式中，根据所述预期影响量矩阵相对应的奇异值矩阵，确定所述不动点延续算法中影响量矩阵的初始矩阵。
20.在本技术的第二方面的某些实施方式中，每一用能对象的实际供能损失参数为所述优化影响量矩阵中每一用能对象各时段的优化供能参数影响量的均值。
21.为实现上述目的及其他相关目的，本技术的第三方面提供一种计算机设备，包括：存储装置，用于存储至少一程序；处理装置，用于调用所述至少一程序，以执行如本技术第一方面中任一项所述的用于供能系统的控制方法。
22.在本技术的第三方面的某些实施方式中，所述计算机设备还包括通信装置，用于与供能对象通信，以输出所述实际供能参数。
23.为实现上述目的及其他相关目的，本技术的第四方面提供一种计算机可读存储介质，存储有至少一程序，所述至少一程序在被调用时执行并实现如本技术第一方面中任一项所述的用于供能系统的控制方法。
24.综上所述，本技术的公开的用于供能系统的控制方法、系统、设备及存储介质，根据每一用能对象在每一时段中各时间间隔内的初始供能需求、以及每一时段中各时间间隔内的供能参数，得到预期影响量矩阵；根据所述预期影响量矩阵、以及矩阵位置集合，在最小化影响量矩阵的秩的条件下构造优化影响量矩阵；根据所述优化影响量矩阵，确定所述供能系统的实际供能参数，以用于控制所述供能系统按该实际供能参数对各所述用能对象供能。本技术可以在考虑用户可能改变负荷曲线的情况下，为电力系统的每一用户定制可以满足其实际用能需求的长期的供能参数。使得供能系统通过长期的供能参数进行供能时，既可以满足用户实际的用能需求又可以使得用户不会主动改变负荷曲线，进而改善了
能源浪费的现象并且保证了电力系统的平稳运行。
附图说明
25.本技术所涉及的发明的具体特征如所附权利要求书所显示。通过参考下文中详细描述的示例性实施方式和附图能够更好地理解本技术所涉及发明的特点和优势。对附图简要说明书如下：
26.图1显示为本技术的应用场景在一实施例中的结构示意图。
27.图2显示为本技术用于供能系统的控制方法在一实施例中的流程示意图。
28.图3显示为根据贪心算法确定单个时段中的预期供能需求的方法在一实施例中的流程示意图。
29.图4显示为本技术用于供能系统的控制系统在一实施例中的结构示意图。
30.图5显示为本技术的计算机设备在一实施例中的结构示意图。
31.图6显示为本技术的计算机设备在另一实施例中的结构示意图。
具体实施方式
32.以下由特定的具体实施例说明本技术的实施方式，熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本技术的其他优点及功效。
33.在下述描述中，参考附图，附图描述了本技术的若干实施例。应当理解，还可使用其他实施例，并且可以在不背离本公开的精神和范围的情况下进行模块或单元组成、电气以及操作上的改变。下面的详细描述不应该被认为是限制性的，并且本技术的实施例的范围仅由公布的专利的权利要求所限定。这里使用的术语仅是为了描述特定实施例，而并非旨在限制本技术。
34.如同在本文中所使用的，单数形式“一”、“一个”和“该”旨在也包括复数形式，除非上下文中有相反的指示。应当进一步理解，术语“包含”、“包括”表明存在所述的特征、步骤、操作、元件、组件、项目、种类、和/或组，但不排除一个或多个其他特征、步骤、操作、元件、组件、项目、种类、和/或组的存在、出现或添加。此处使用的术语“或”和“和/或”被解释为包括性的，或意味着任一个或任何组合。因此，“a、b或c”或者“a、b和/或c”意味着“以下任一个：a；b；c；a和b；a和c；b和c；a、b和c”。仅当元件、功能、步骤或操作的组合在某些方式下内在地互相排斥时，才会出现该定义的例外。
35.电力系统是由发电厂、送变电线路、供配电所和用电等环节组成的电能生产与消费系统。它的功能是将自然界的一次能源通过发电动力装置转化成电能，再经输电、变电和配电将电能供应到各用户。为实现这一功能，电力系统还需要包括控制系统，以对电能的生产和分配的过程进行测量、调节、控制、保护、通信、以及调度等。
36.电力系统的主体结构有电源、变电所、输电/配电线路、以及负荷中心。所述电源包括但不限于水电站、火电厂、以及核电站等发电厂。根据能源类型的不同，电力系统的类型包括但不限于：基于传统燃料发电的电网、基于新能源发电的微电网等。所述新能源包括但不限于：煤炭、天然气、风能、水能、光伏、生物质能、海洋能、以及其他形式等中的一种或多种方式。
37.电能的生产与供应使用应当保持平衡。现有技术通常基于用户的负荷曲线来获取
未来一段时间内的负荷预测，从而根据预测结果为不同的用户制定不同的供能策略以满足用户的负荷。然而，用户的能源消耗模式、即负荷曲线往往会随着时间的变化而变化。例如，用户在拥有储能装置后会进行存储控制，这样可以在满足其用能需求的情况下实现降低其用能损失量的目的，或者用户直接通过改变其用能行为来实现降低其用能损失量的目的。这样的情况下，会使预测的供能策略与用户的实际用能不匹配，从而造成能源浪费的现象并且不利于整个电力系统的发展和平稳运行。
38.所以，为了更高效地控制供能系统的运行，在考虑用户可能改变负荷曲线的情况下有针对性地为不同用户确定稳定的供能参数(例如供能量、供能量对应的供能指令、或者供能损失参数)，能够避免套利者获利，一方面提高对供能系统运行的控制效率，另一方面增加供能系统的可靠性和稳定性。
39.基于此，本技术提供一种用于供能系统的控制系统，以获得优化的实际供能参数。
40.请参阅图1，图1显示为本技术的应用场景在一实施例中的结构示意图。
41.如图1所示，所述供能系统100向多个用能对象200供能。其中，所述供能系统可供应如电能、热能等能源；所述供能系统可以被应用在家庭、企业、商场、学校等任何用能场所。
42.其中，所述用能对象200指一个或多个需要用能的设备，例如工业设备，例如制造/加工设备、流水线设备等；又例如家用电器设备，例如电视、冰箱、空调等，又例如照明系统等的灯具。在一些示例中，所述用能对象可以是接入同一计量装置的至少一个用能设备的总称。
43.所述供能系统的类型以及向用能对象供能的传输介质可根据能源的类型而确定。所述能源的类型包括但不限于为电能、热能等。例如，当所述能源为电能时，所述供能系统包括但不限于为电网，所述传输介质包括但不限于为电线。
44.本技术的应用场景可以是传统的火力发电电网，也可以是基于分布式能源的微电网等，此处不作限制。
45.在一些实施例中，用能对象一侧还可以包括储能装置，所述储能装置用于储备一定容量的能源。所述储能装置包括但不限于为蓄电池、超级电容器等。其中，所述蓄电池可为铅酸电池、镍镉电池、镍氢电池、钠硫电池、锂电池或燃料电池等。
46.可以理解的是，由于并非限制能源类型为电能，因此对其表述亦不限制于电能，在本文中通过“能源”、以及“能量”来表达；相应的，电网系统、电力系统等，也通过“供能系统”来替代性地表达。
47.需说明的是，所述供能系统只是一种概括性的表达，而非限制其必须为一独立的实体装置。所述供能系统和用能对象之间的耦合关系包括直接或间接的连接关系，只要能令用能对象从供能系统获得能量即可，而并不限制其具体实现形式。
48.请参阅图2，图2显示为本技术用于供能系统的控制方法在一实施例中的流程示意图。所述控制方法用于在考虑用能对象的实际负荷曲线相对于历史负荷曲线可能发生改变的情况下，有针对性地为不同的用能对象确定长期稳定的供能参数。所述供能参数包括但不限于：每一用能对象的供能量、供能量对应的供能指令、或者供能损失参数，用于稳定能源系统中各用户的用能行为，令其合理化，且避免能源浪费。在一些实施例中，所述供能损失参数可以是表示对供能系统进行供能所产生损失进行补偿的参数。
49.所述用于供能系统的控制方法可由本技术所述的计算机设备执行，或者其他可执行本技术的控制方法的计算机设备执行。
50.如图2所示，所述控制方法包括步骤s110、步骤s120、步骤s130。
51.步骤s110：所述计算机设备根据每一用能对象在每一时段中各时间间隔内的初始供能需求、以及每一时段中各时间间隔内的供能参数，得到预期影响量矩阵；其中，所述预期影响量矩阵中包括每一用能对象在每一所述时段中对于所述供能系统的预期供能参数影响量；在初始影响量矩阵上叠加和初始供能需求相关的噪声矩阵为所述预期影响量矩阵。
52.在以下实施例中，根据所述供能损失参数进行举例。
53.在一些实施例中，所述供能参数影响量等数据是在一时段内进行衡量、评估、拟合、或计算的。其中，所述时段是指用能对象的用能时长；所述用能时长包括但不限于：一或多日、一或多周、一或多月、或者一或多年等。
54.其中，每个所述时段包括一个或具有不同时间特征信息的多个周期性的时间间隔。举例来说，设所述时段为一天，所述时间间隔为一天中的每个小时。或者，在一天中按供能系统的供能损失参数的不同可以分为多个时间间隔，即：0am～7am时间间隔、7am～11pm时间间隔、11pm～5pm时间间隔、5pm～7pm时间间隔、7pm～12pm时间间隔。
55.需要说明的是，所述时段内各时间间隔的时序分配仅为举例，并不以此为限。
56.在本技术的示例中，以t表示所述时段中时间间隔的总数，以t表示所述时间间隔；其中，t的可行域为[1，t]。例如，所述时间间隔为每个小时，则所述时段可以为一天，即t可以为24。
[0057]
在实际场景中，用能对象i之间的用能量各有不同，其中，为一时段中每个时间间隔t内用能对象i的用能量。所述供能参数影响量能够评估用能对象的用能量对于供能系统的供能损失量的边际影响。
[0058]
在一些实施例中，每一用能对象在每一时段中对于所述供能系统的供能参数影响量是根据每一用能对象在每一时段中各时间间隔内的供能需求和每一时段中各时间间隔内的供能参数所确定的。
[0059]
公式(1)给出了每个用能对象i的供能参数影响量mci
i
的定义：
[0060][0061]
在上述公式中，为在该时段内该用能对象i的总用能量；其中，为在该时段内该用能对象i的总用能量；其中，是根据在该时段内每个时间间隔t内用能对象i的用能量计算获得的。
[0062]
需说明的是，在其它实施例中，所述用能量和总用能量的计算方式也可以加以变化，例如，每一时段内总用能量可以是各时间间隔的用能量的加权和，每一时间间隔的用能量可以是各时刻的用能量的加权和等，并非以此为限。
[0063]
在一示例中，供能损失量和用能量的数学关系由公式(2)所示的损失量函数表示：
[0064][0065]
其中，n表示供能系统所供能的用能对象的数量。所述损失量函数用于表示在时间间隔t内n个用能对象的用能量l
t
和供能系统的供能损失量c(l
t
)之间的关系；为在时间间隔t内用能对象i的用能量，即供能系统在时间间隔t内对用能对象i的供能量；a为损失量函数的二次项系数、b为损失量函数的一次项系数、c为损失量函数的常数项。
[0066]
当然，虽然本实施例中通过二次函数来表示所述n个用能对象的用能量l
t
和供能系统的供能损失量c(l
t
)之间的数学关系，但是并不限制于此，在其它实施例中可以通过其它多项式形式加以替代。
[0067]
在一些实施例中，所述供能损失参数举例为在各时间间隔内供能系统的供能价格p(t)。
[0068]
在一示例中，所述各时间间隔内的供能损失参数p(t)可以由公式(3)表示：
[0069][0070]
在一示例中，综合公式(1)至公式(3)，所述供能参数影响量mci
i
可由如下公式表示：
[0071][0072]
在公式(4)中，p(t)为每个时间间隔内的供能损失参数，为每个时间间隔t内用能对象i的用能量，为的1-范数(1-norm)。
[0073]
在一些实施例中，在一时段中各时间间隔内用能对象i的供能需求可以直接用其用能量来表示，其中，为该时段中每个时间间隔内t该用能对象i的用能量。在另一些实施例中，在一时段中各时间间隔内用能对象i的供能需求也可以用归一化的用能
量来表示，其中，来表示，其中，为用能量的归一化表达。
[0074]
基于公式(4)，可根据一时段中各时间间隔内用能对象i的供能需求和各时间间隔内的供能损失参数p
t
确定所述用能对象i的供能参数影响量mci
i
。
[0075]
需说明的是，在其它实施例中，所述供能参数影响量mci
i
的计算方式也可以加以变化，例如在公式(4)的基础上增加其它权重系数等，并非以此为限。
[0076]
基于上述理解，所述步骤s110包括步骤s111和步骤s112。
[0077]
步骤s111：所述计算机设备根据每一用能对象在每一时段中各时间间隔内的初始供能需求、以及每一时段中各时间间隔内的供能参数，得到每一用能对象在每一时段中各时间间隔内的预期供能需求。
[0078]
所述初始供能需求为真实用能量的历史数据所对应的供能需求。所述各时间间隔内的供能参数举例为各时间间隔内的供能损失参数。所述预期供能需求包括预计可能的供能需求，例如，在实际的用能过程中，为了使用能对象在一时段中的总用能量相同、以及在所述供能损失参数的作用下降低该时段中的用能损失量，会对在该时段的每一时间间隔内的用能量进行调整。所述用能损失量举例为用能成本。所述用能损失量用于补偿所述供能系统的供能损失量。
[0079]
具体地，基于所述初始供能需求以及所述各时间间隔内的供能损失参数可以对每一时间间隔内的供能需求进行预期调整，以使得在所述供能损失参数的作用下降低在该时段中的用能损失量。需要说明的是，将用能对象的用能量作归一化处理后作为所述初始供能需求，可以提升所述用于供能系统的控制方法的计算效率。实际应用中，也可以采用各时间间隔t内该用能对象i的用能量作为初始供能需求来获取预期供能需求。
[0080]
在一些实施例中，所述初始供能需求与预期供能需求满足如下公式：
[0081][0082][0083]
公式(5)表示的是单个用能对象在单个时段中各时间间隔的初始供能需求的总量与预期供能需求的总量相等。公式(6)表示的是单个用能对象在单个时段中各时间间隔内的初始供能需求同各时间间隔内的预期供能需求的差异之和满足预设比例θ。其中，所述预期供能需求为
[0084]
所述初始供能需求与预期供能需求在满足公式(5)和公式(6)所述的约束的情况下，根据单个用能对象在单个时段中各时间间隔内的初始供能需求、以及单个时段各时间间隔内的供能参数，确定使单个用能对象在单个时段中的用能损失量最小所对应的单个用能对象在单个时段中各时间间隔内的预期供能需求。
[0085]
以所述供能参数为供能损失参数p
t
为例，通过在公式(5)和公式(6)的约束下，确定使得公式(7)所示的用能损失量最小所对应的各时间间隔内的预期供能需求d
i
'
t
。
[0086][0087]
需要说明的是，也可以为公式(5)和公式(6)设置期限。例如，根据实际需求，在一时段中部分时间间隔内实施所述用能量改变的行为。即，求解问题(7)时可以采取分段方式求解，或者只在一时段中连续的几个时间间隔内进行求解。
[0088]
在另一些实施例中，也可以采用各时间间隔t内一用能对象i的用能量作为所述初始供能需求来求解预期供能需求公式(5)、公式(6)、以及公式(7)可以分别作出适应性改变，转化成如公式(8)、公式(9)、以及公式(10)所示的优化问题：
[0089][0090][0091][0092]
其中，在公式(8)和公式(9)的约束下，确定使得公式(10)所示的用能损失量最小所对应的各时间间隔内的用能量可以直接将得到的各时间间隔内的用能量作为所述预期供能需求来确定所述预期供能参数影响量。
[0093]
在一些示例中，可以根据单个用能对象在单个时段中各时间间隔内的初始供能需求、以及该时段中各时间间隔内的供能损失参数，应用贪心算法来求解公式(7)和公式(10)所对应的优化问题，以确定所述预期供能需求。
[0094]
具体地，以用能对象i在单个时段中归一化的用能量作为所述初始供能需求为例。请参阅图3，图3显示为根据贪心算法确定单个时段中的预期供能需求的方法在一实施例中的流程示意图。图3所示的方法包括：步骤s121、步骤s122、步骤s123、步骤s124、步骤s125。
[0095]
步骤s121：对所述单个时段中的每一时间间隔内的供能损失参数进行排序，以得到升序排列的损失参数排序p(1)＜p(2)＜.....p(t)，并将当前供能改变量δd的初始值设为0，将当前供能损失参数顺序j的初始值设为1，以及确定供能改变量阈值。在本实施例中，所述供能改变量阈值为θ
·
sum1，其中，其中在其他实施例中，如果所述初始供能需求直接用一时段的用能量来表示，则所述供能该变量阈值所涉及的sum1可以直接用前文所述的来替代。
[0096]
步骤s122：判断当前供能改变量δd是否小于等于所述供能改变量阈值θ
·
sum1，如果是，则执行步骤s123，如果否，则结束。其中，所述当前供能改变量δd会根据步骤s125
进行更新。
[0097]
步骤s123：判断与所述当前供能损失参数顺序j相对应的高供能损失参数顺序的初始供能需求d
i(t+1-j)
是否在所述供能改变量阈值与当前供能改变量的差值θ
·
sum
1-δd以上：
[0098]
如果是，执行步骤s124，如果否，则执行步骤s125。其中，所述当前供能损失参数顺序j会根据步骤s125进行更新。
[0099]
步骤s124：将当前供能损失参数顺序j下的初始供能需求d
i(j)
加上所述供能改变量阈值同所述当前供能改变量的差值θ
·
sum
1-δd后作为当前供能损失参数顺序下的预期供能需求d
i
'
(j)
，将与所述当前供能损失参数顺序相对应的高供能损失参数顺序的初始供能需求d
i(t+1-j)
减去所述供能改变量阈值与所述当前供能改变量的差值θ
·
sum
1-δd后作为当前高供能损失参数顺序下的所述预期供能需求d
i
'
(t+1-j)
。然后结束获得该用能对象在该时段的预期供能需求的操作。
[0100]
例如，一时段中的初始供能需求为该时段中六个时间间隔所对应的供能损失参数分别为(1，3，6，5，4，2)。对该时段中的供能损失参数按损失参数升序的到的损失参数排序为：1<2<3<4<5<6，设θ值为0.2，则供能损失参数最低的时间间隔所对应的初始供能需求为0.05。第一次进行判断时，当前供能损失参数顺序j＝1相对应的高供能损失参数顺序j＝6的初始供能需求为0.4。由于0.4>(0.2-0)，所以，将供能损失参数顺序j＝6所对应的供能需求从0.4变成0.2，以0.2为供能损失参数顺序j＝6所对应的时间间隔内的预期供能需求。相应地，将供能损失参数顺序j＝1所对应的供能需求从0.05变成0.25，以0.25为供能损失参数顺序j＝1所对应的时间间隔内的预期供能需求。将(0.25，0.25，0.2，0.15，0.1，0.05)作为该时段的预期供能需求
[0101]
步骤s125：将与所述当前供能损失参数相对应的高供能损失参数下的初始供能需求d
i(t+1-j)
设为0后作为当前高供能损失参数下的预期供能需求d
i
'
(t+1-j)
，将与所述当前供能损失参数下的初始供能需求d
i(j)
加上与所述当前供能损失参数相对应的高供能损失参数下的初始供能需求d
i(t+1-j)
后作为当前供能损失参数下的预期供能需求d
i
'
(j)
，并将当前供能损失参数顺序j加1，将当前供能改变量同与所述当前供能损失参数相对应的高供能损失参数下的初始供能需求之和d
i(t+1-j)
+δd作为下一当前供能改变量。将更新后的当前供能改变量δd和当前供能损失参数顺序j分别作为下一次执行步骤s122时的当前供能改变量δd和下一次执行步骤s123时的当前供能损失参数顺序j。
[0102]
例如，一时段中的初始供能需求为该时段中十个时间间隔所对应的供能损失参数分别为(1，8，7，11，6，9，5，4，3，2)。对该时段中的供能损失参数按损失参数升序得到的损失参数排序为：1<2<3<4<5<6<7<8<9<11，设θ值为0.3，则供能损失参数最低的时间间隔所对应的初始供能需求为0，第一次进行判断时，当前供能损失参数顺序j＝1相对应的高供能损失参数顺序j＝10的供能损失参数为11，当前供能损失参数顺序j＝1相对应的高供能损失参数顺序j＝10的初始供能需求为0.25。由于0.25<(0.3-0)，所以，将供能损失参数顺序j＝10所对应的供能需求从0.25变成0，以0为供能损失参数顺序j＝10所对应的预期供能需求。相应地，将供能损失参数顺序j＝1所对应的供能需求从0变成(0.25+0)，以0.25为供能
损失参数顺序j＝1所对应的预期供能需求。然后，将供能损失参数顺序j加1，此时所述当前供能损失参数顺序j为2。将供能改变量δd加上供能损失参数顺序j＝10所对应的供能需求，此时，当前供能改变量δd为0.25。由于当前供能损失参数顺序j＝2相对应的高供能损失参数顺序j＝9的初始供能需求为0.25>(0.3-0.25)，所以，将供能损失参数顺序j＝2所对应的供能需求从0变成((0.3-0.25)+0)，以0.05为供能损失参数顺序j＝2所对应的预期供能需求。将供能损失参数顺序j＝9所对应的供能需求从0.25变成(0.25-(0.3-0.25))，以0.2为供能损失参数顺序j＝9所对应的预期供能需求。将(0.25，0.05，0.1，0，0.1，0.2，0.15，0.05，0.05，0.05)作为该时段的预期供能需求
[0103]
其中，所述当前供能改变量初始时为所述初始值0，当前供能损失参数顺序j初始时为1。所述当前供能损失参数初始时所对应的供能损失参数p(1)最低。其中，所述当前供能改变量δd和当前供能损失参数顺序j会根据步骤s125进行更新，以作用至步骤s122和步骤s123。
[0104]
在确定每一用能对象i在每一时段中的预期供能需求后，可执行步骤s112。
[0105]
步骤s112：根据每一用能对象在每一时段中各时间间隔内的预期供能需求，确定每一用能对象在每一时段中对于所述供能系统的预期供能参数影响量，以得到所述预期影响量矩阵。
[0106]
具体地，将多个用能对象在多个时段中的预期供能需求带入公式(4)，获得多个用能对象在多个时段中的预期供能参数影响量，以得到所述预期影响量矩阵y。其中，所述预期影响量矩阵y中的每一个元素代表每一用能对象在每一所述时段中对于所述供能系统的预期供能参数影响量。
[0107]
其中，所述预期影响量矩阵中的每一个元素代表每一用能对象在每一所述时段中对于所述供能系统的预期供能参数影响量y
i,j
。所述预期供能参数影响量y
i,j
满足如下公式：
[0108][0109]
其中，为每一用能对象在每一时段的所述初始供能需求所对应的供能参数影响量，ε
i,j
为每一用能对象在每一时段的噪声所对应的供能参数影响量。m为用能对象的个数，n为观察期间中所述时段的个数。其中，所述噪声是所述预期供能需求与所述初始供能需求的差异。
[0110]
所以，在初始影响量矩阵上叠加和初始供能需求相关的噪声矩阵为所述预期影响量矩阵，所述预期影响量矩阵y可由如下公式表示：
[0111]
y＝θ
*
+z
ꢀꢀꢀꢀꢀꢀ
(12)
[0112]
其中，θ
*
为初始影响量矩阵，所述初始影响量矩阵θ
*
中的每一元素为每一用能对象在每一时段的所述初始供能需求所对应的供能参数影响量z为和初始供能需求相关的噪声矩阵，所述噪声矩阵z中的每一元素为每一用能对象在每一时段的噪声所对应的供能参数影响量ε
i,j
。
[0113]
需要说明的是，实际上用户不会在整个观察期间的每个时段中均进行供能需求的改变。所以，所述在实际中，预期影响量矩阵中只有部分的供能参数影响量是存在噪声的，
即而预期影响量矩阵中其他供能参数影响量中噪声所对应的供能参数影响量ε
i,j
为0。在这种情况下，初始影响量矩阵θ
*
是低秩性，预期影响量矩阵y也是低秩的。所以，我们可以采用前文求解得到的所述预期影响量矩阵y中的部分元素来构造优化影响量矩阵，使得所述优化影响量矩阵所对应的实际供能参数更接近用户的真实供能需求。所以，供能系统采用基于该优化影响量矩阵得到的长期稳定的实际供能参数进行供能时，用户不会为了降低用能损失量而去改变自身的用能需求。
[0114]
步骤s120可以根据前文求解得到的所述预期影响量矩阵y中的部分元素来构造优化影响量矩阵。
[0115]
步骤s120：所述计算机设备根据所述预期影响量矩阵、以及矩阵位置集合，在最小化影响量矩阵的秩的条件下构造优化影响量矩阵；其中，所述优化影响量矩阵中包括每一用能对象在每一时段中对于所述供能系统的优化供能参数影响量；所述影响量矩阵和预期影响量矩阵是同型矩阵，且所述矩阵位置集合中包括所述预期影响量矩阵及所述影响量矩阵中部分的相应元素位置。
[0116]
在一些实施例中，利用矩阵补全算法，在最小化影响量矩阵的秩的条件下和预设误差约束下构造所述优化影响量矩阵。所述误差约束可由如下公式表示：
[0117]
||p
ω
(θ-y)||
f
＜δ
ꢀꢀꢀꢀꢀꢀ
(13)
[0118]
其中，δ为预设误差，p
ω
为投影操作。公式(13)表示的是根据矩阵位置集合对影响量矩阵θ和预期影响量矩阵y的差进行投影操作得到的矩阵的范数满足预设误差。所述投影操作p
ω
可由公式(14)来解释：
[0119][0120]
其中，ω为矩阵位置集合，所述影响量矩阵θ和预期影响量矩阵y是同型矩阵，且所述矩阵位置集合ω中包括所述预期影响量矩阵y及所述影响量矩阵θ中部分的相应元素位置。例如，矩阵位置集合ω中包括m
×
n的预期影响量矩阵y的对角线位置，由于所述影响量矩阵θ和预期影响量矩阵y是同型矩阵，所以，所述矩阵位置集合ω中还包括所述影响量矩阵θ的对角线位置。
[0121]
所述矩阵位置集合ω是通过抽样算法和预计抽样的矩阵位置个数n确定的。例如，所述抽样算法为简单随机抽样，所述预计抽样的矩阵位置个数为n，基于所述抽样概率抽取矩阵中的n个位置。其中，m与n的乘积为所述预期影响量矩阵中矩阵位置的总个数。需要说明的是，抽样的个数太多，会使构造的矩阵更接近y；抽样的个数太少，会使得无法准确构造所述优化影响量矩阵。
[0122]
在所述预设误差约束下和最小化影响量矩阵的秩的条件下，构造优化影响量矩阵的问题可以转化为公式(15)所对应的秩最小化的优化问题。
[0123][0124]
利用矩阵补全算法可以求解上述优化问题以构造所述优化影响量矩阵。其中，所述优化影响量矩阵中包括每一用能对象在每一时段中对于所述供能系统的优化供能参数
影响量。
[0125]
所述矩阵补全算法包括但不限于不动点延续算法、加速近端梯度算法、半定规划法、奇异值阈值算法、增广拉格朗日乘子法等。
[0126]
所述公式(15)所对应的非凸优化问题可以转化成如下公式所示的最小值问题：
[0127][0128]
其中，||θ||
*
为所述影响量矩阵θ的奇异值之和，||.||
f
代表矩阵的f范数，矩阵的f范数是每一个矩阵位置的元素值的平方求和后的次幂，其中，m为所述用能对象的个数，n为观察期间中所述时段的个数，p为确定所述矩阵位置集合ω的抽样概率，σ是所述噪声矩阵z的标准差。所述噪声矩阵可以根据前文求解得到的所述预期影响量矩阵y与所述为初始影响量矩阵θ
*
之间的差值得到。
[0129]
利用不动点延续算法求解公式(16)的最小值问题时，需要获得如下微分方程的解：
[0130][0131]
其中，表示的是关于所述影响量矩阵θ的次梯度。不动点延续算法构造矩阵a、和b，并将b矩阵的转置与a矩阵的乘积ab
t
作为所述优化影响量矩阵。
[0132]
其中，所述矩阵a和矩阵b是基于初始矩阵a0和初始矩阵b0，利用所述不动点算法不断迭代得到的。需要说明的是，初始矩阵a0和初始矩阵b0的设置与所述不动点算法的收敛速度有关。
[0133]
在一些实施例中，所述初始矩阵a0和初始矩阵b0与前文求解得到的所述预期影响量矩阵y相关。
[0134]
具体地，对所述预期影响量矩阵y进行奇异值分解可以得到所述预期影响量矩阵y的奇异值矩阵σ
m
×
n
，奇异值矩阵σ
m
×
n
是m
×
n的对角矩阵，其中，通过简化后的奇异值矩阵可以构造所述初始矩阵a0和初始矩阵b0。
[0135]
例如，所述奇异值矩阵中的各对角元素为(σ1，σ2,...,σn)。将所述奇异值矩阵σ
m
×
n
中的各对角元素按照大小排序，得到排序后的奇异值矩阵，排序后的奇异值矩阵中的各对角元素为(σ1'，σ2',...,σn')，其中，(σ1'＞σ2'＞...＞σn')。当所述奇异值矩阵中的对角元素小于奇异值阈值时则将对应的对角元素设为0，以得到简化的奇异值矩阵所述奇异值阈值可以是所述对角线元素最大值的十分之一、百分之一、或者千分之一等，但并不以此为限。所述简化后的奇异值矩阵可由如下公式表示：
[0136][0137]
其中，k为简化后的奇异值矩阵中非0的对角元素的个数。是对角矩阵，每一对角元素为奇异值矩阵σ
m
×
n
中在奇异值阈值以上的对角元素。0
m
×
n
代表的是m
×
n维的
0矩阵。
[0138]
根据简化的奇异值矩阵得到的所述初始矩阵a0和初始矩阵b0可由如下公式表示：
[0139][0140]
其中，是对所述对角矩阵中的每个元素开根号得到的。为所述优化影响量矩阵的初始矩阵。
[0141]
在其他实施例中，所述优化影响量矩阵的初始矩阵可以为m
×
n维的0矩阵，也可以是所述预期影响量矩阵y，但不以此为限。
[0142]
需要说明的是，在步骤s120中，以用能对象的个数m为行并且以观察期间中所述时段的个数n为列。在其他实施方式中，也可用观察期间中所述时段的个数n为行并且以用能对象的个数m为列。
[0143]
根据步骤s120中构造的优化影响量矩阵可执行步骤s130。
[0144]
步骤s130：根据所述优化影响量矩阵，确定所述供能系统的实际供能参数，以用于控制所述供能系统按该实际供能参数对各所述用能对象供能。
[0145]
所述实际供能参数包括但不限于：每一用能对象的实际供能量、实际供能量对应的实际供能指令、或者实际供能损失参数。所述实际供能参数用于稳定能源系统中各用户的用能行为，令其合理化，且避免能源浪费。
[0146]
在一些实施例中，所述实际供能损失参数为所述优化影响量矩阵中每一用能对象各时段的优化供能参数影响量的均值。例如，所述优化影响量矩阵中一用能对象在n个时段中的优化供能参数影响量为(a1,a2,....,a
n
)，对各优化供能参数影响量求和后与n的比值作为所述稳定的实际供能损失参数。
[0147]
在一些实施例中，基于得到的所述实际供能损失参数，确定所述每一用能对象的实际供能量、或实际供能量对应的实际供能指令。例如，将与所述实际供能损失参数相对应的供能量作为长期的实际供能量，以控制不包括储能装置的供能系统按照实际供能量进行产能，以供给所述各用能对象。又如，将与所述实际供能损失参数相对应的供能量作为稳定的实际供能量后确定所述实际供能指令，以供包括储能装置的供能系统按照所述实际供能指令进行产能、储能、以及释能。其中，所述实际供能指令中包括：每一时间间隔的实际产能量、每一时间间隔的实际储能量、每一时间间隔的实际释能量。
[0148]
需要说明的是，供能系统通过所述实际供能参数进行供能时，既可以满足用户实际的用能需求又可以使得用户不会轻易改变其负荷曲线，进而改善了能源浪费的现象并且保证了电力系统的平稳运行。
[0149]
本技术还提供一种用于供能系统的控制系统，用于在考虑用能对象的实际负荷曲线相对于历史负荷曲线可能发生改变的情况下，有针对性地为不同的用能对象确定长期稳定的供能参数。所述供能参数包括但不限于：每一用能对象的供能量、供能量对应的供能指令、或者供能损失参数，用于稳定能源系统中各用户的用能行为，令其合理化，且避免能源
浪费。在一些实施例中，所述供能损失参数可以是表示对供能系统进行供能所产生损失进行补偿的参数。
[0150]
所述控制系统可以应用在图5或图6所示的计算机系统中。例如，以软件的形式搭载在图5或图6所示的计算机系统中以实现相应功能。
[0151]
请参阅图4，图4显示为本技术用于供能系统的控制系统在一实施例中的结构示意图。如图所4示，所述用于供能系统的控制系统40包括：预期影响量矩阵确定模块41、优化影响量矩阵确定模块42、供能参数确定模块43。
[0152]
所述预期影响量矩阵确定模块41用于根据每一用能对象在每一时段中各时间间隔内的初始供能需求、以及每一时段中各时间间隔内的供能参数，得到预期影响量矩阵；其中，所述预期影响量矩阵中包括每一用能对象在每一所述时段中对于所述供能系统的预期供能参数影响量；在初始影响量矩阵上叠加和初始供能需求相关的噪声矩阵为所述预期影响量矩阵。
[0153]
在以下实施例中，根据所述供能损失参数进行举例。
[0154]
在实际场景中，用能对象i之间的用能量各有不同，其中，为一时段中每个时间间隔t内用能对象i的用能量。所述供能参数影响量能够评估用能对象的用能量对于供能系统的供能损失量的边际影响。
[0155]
在一些实施例中，每一用能对象在每一时段中对于所述供能系统的供能参数影响量是根据每一用能对象在每一时段中各时间间隔内的供能需求和每一时段中各时间间隔内的供能参数所确定的。
[0156]
公式(1)给出了每个用能对象i的供能参数影响量mci
i
的定义：
[0157][0158]
在上述公式中，为在该时段内该用能对象i的总用能量；其中，为在该时段内该用能对象i的总用能量；其中，是根据在该时段内每个时间间隔t内用能对象i的用能量计算获得的。
[0159]
需说明的是，在其它实施例中，所述用能量和总用能量的计算方式也可以加以变化，例如，每一时段内总用能量可以是各时间间隔的用能量的加权和，每一时间间隔的用能量可以是各时刻的用能量的加权和等，并非以此为限。
[0160]
在一示例中，供能损失量和用能量的数学关系由公式(2)所示的损失量函数表示：
[0161][0162]
其中，n表示供能系统所供能的用能对象的数量。所述损失量函数用于表示在时间间隔t内n个用能对象的用能量l
t
和供能系统的供能损失量c(l
t
)之间的关系；为在时间间隔t内用能对象i的用能量，即供能系统在时间间隔t内对用能对象i的供能量；a
为损失量函数的二次项系数、b为损失量函数的一次项系数、c为损失量函数的常数项。
[0163]
当然，虽然本实施例中通过二次函数来表示所述n个用能对象的用能量l
t
和供能系统的供能损失量c(l
t
)之间的数学关系，但是并不限制于此，在其它实施例中可以通过其它多项式形式加以替代。
[0164]
在一些实施例中，所述供能损失参数举例为在各时间间隔内供能系统的供能价格p(t)。在一示例中，所述各时间间隔内的供能损失参数p(t)可以由公式(3)表示：
[0165][0166]
在一示例中，综合公式(1)至公式(3)，所述供能参数影响量mci
i
可由如下公式表示：
[0167][0168]
在公式(4)中，p(t)为每个时间间隔内的供能损失参数，为每个时间间隔t内用能对象i的用能量，为的1-范数(1-norm)。
[0169]
在一些实施例中，在一时段中各时间间隔内用能对象i的供能需求可以直接用其用能量来表示，其中，为该时段中每个时间间隔内t该用能对象i的用能量。在另一些实施例中，在一时段中各时间间隔内用能对象i的供能需求也可以用归一化的用能量来表示，其中，来表示，其中，为用能量的归一化表达。
[0170]
基于公式(4)，可根据一时段中各时间间隔内用能对象i的供能需求和各时间间隔内的供能损失参数p
t
确定所述用能对象i的供能参数影响量mci
i
。
[0171]
需说明的是，在其它实施例中，所述供能参数影响量mci
i
的计算方式也可以加以变化，例如在公式(4)的基础上增加其它权重系数等，并非以此为限。
[0172]
基于上述理解，首先，所述预期影响量矩阵确定模块41根据每一用能对象在每一时段中各时间间隔内的初始供能需求、以及每一时段中各时间间隔内的供能参数，得到每一用能对象在每一时段中各时间间隔内的预期供能需求。
[0173]
所述初始供能需求为真实用能量的历史数据所对应的供能需求。所述各时间间隔内的供能参数举例为各时间间隔内的供能损失参数。所述预期供能需求包括预计可能的供能需求，例如，在实际的用能过程中，为了使用能对象在一时段中的总用能量相同、以及在所述供能损失参数的作用下降低该时段中的用能损失量，会对在该时段的每一时间间隔内的用能量进行调整。所述用能损失量举例为总用能成本。所述用能损失量用于补偿所述供能系统的供能损失量。
[0174]
具体地，基于所述初始供能需求以及所述各时间间隔内的供能损失参数可以对每一时间间隔内的供能需求进行预期调整，以使得在所述供能损失参数的作用下降低在该时段中的用能损失量。需要说明的是，将用能对象的用能量作归一化处理后作为所述初始供能需求，可以提升计算效率。实际应用中，也可以采用各时间间隔t内该用能对象i的用能量作为初始供能需求来获取预期供能需求。
[0175]
在一些实施例中，所述初始供能需求与预期供能需求满足如下公式：
[0176][0177][0178]
公式(5)表示的是单个用能对象在单个时段中各时间间隔的初始供能需求的总量与预期供能需求的总量相等。公式(6)表示的是单个用能对象在单个时段中各时间间隔内的初始供能需求同各时间间隔内的预期供能需求的差异之和满足预设比例θ。其中，所述预期供能需求为
[0179]
所述初始供能需求与预期供能需求在满足公式(5)和公式(6)所述的约束的情况下，根据单个用能对象在单个时段中各时间间隔内的初始供能需求、以及单个时段各时间间隔内的供能参数，确定使单个用能对象在单个时段中的用能损失量最小所对应的单个用能对象在单个时段中各时间间隔内的预期供能需求。
[0180]
以所述供能参数为供能损失参数p
t
为例，通过在公式(5)和公式(6)的约束下，确定使得公式(7)所示的用能损失量最小所对应的各时间间隔内的预期供能需求d
i
'
t
。
[0181][0182]
需要说明的是，也可以为公式(5)和公式(6)设置期限。例如，根据实际需求，在一时段中部分时间间隔内实施所述用能量改变的行为。即，求解问题(7)时可以采取分段方式求解，或者只在一时段中连续的几个时间间隔内进行求解。
[0183]
在另一些实施例中，也可以采用各时间间隔t内一用能对象i的用能量作为所述初始供能需求来求解预期供能需求公式(5)、公式(6)、以及公式(7)可以分别作出适应性改变，转化成如公式(8)、公式(9)、以及公式(10)所示的优化问题：
[0184][0185][0186][0187]
其中，在公式(8)和公式(9)的约束下，确定使得公式(10)所示的用能损失量最小所对应的各时间间隔内的用能量可以直接将得到的各时间间隔内的用能量作为所述预期供能需求来确定所述预期供能参数影响量。
[0188]
在一些示例中，可以根据单个用能对象在单个时段中各时间间隔内的初始供能需求、以及该时段中各时间间隔内的供能损失参数，应用贪心算法来求解公式(7)和公式(10)所对应的优化问题，以确定所述预期供能需求。
[0189]
具体地，以用能对象i在单个时段中归一化的用能量作为所述初始供能需求为例。所述预期影响量矩阵确定模块41根据贪心算法确定单个时段中的预期供能需求。
[0190]
需要说明的是，根据贪心算法确定单个时段中的预期供能需求的方式与前文所述的相同或相似，在此不再详述。
[0191]
在确定每一用能对象i在每一时段中的预期供能需求后，所述预期影响量矩阵确定模块41根据每一用能对象在每一时段中各时间间隔内的预期供能需求，确定每一用能对象在每一时段中对于所述供能系统的预期供能参数影响量，以得到所述预期影响量矩阵。
[0192]
将多个用能对象在多个时段中的预期供能需求带入公式(4)，获得多个用能对象在多个时段中的预期供能参数影响量，以得到所述预期影响量矩阵y。其中，所述预期影响量矩阵y中的每一个元素代表每一用能对象在每一所述时段中对于所述供能系统的预期供能参数影响量。
[0193]
其中，所述预期影响量矩阵中的每一个元素代表每一用能对象在每一所述时段中对于所述供能系统的预期供能参数影响量y
i,j
。所述预期供能参数影响量y
i,j
满足如下公式：
[0194][0195]
其中，为每一用能对象在每一时段的所述初始供能需求所对应的供能参数影响量，ε
i,j
为每一用能对象在每一时段的噪声所对应的供能参数影响量。m为所述用能对象的个数，n为观察期间中所述时段的个数。其中，所述噪声是所述预期供能需求与所述初始供能需求的差异。
[0196]
所以，在初始影响量矩阵上叠加和初始供能需求相关的噪声矩阵为所述预期影响量矩阵，所述预期影响量矩阵y可由如下公式表示：
[0197]
y＝θ
*
+z
ꢀꢀꢀꢀꢀꢀ
(12)
[0198]
其中，θ
*
为初始影响量矩阵，所述初始影响量矩阵θ
*
中的每一元素为每一用能对象在每一时段的所述初始供能需求所对应的供能参数影响量z为和初始供能需求相
关的噪声矩阵，所述噪声矩阵z中的每一元素为每一用能对象在每一时段的噪声所对应的供能参数影响量ε
i,j
。
[0199]
需要说明的是，实际上用户不会在整个观察期间的每个时段中均进行供能需求的改变。所以，所述在实际中，预期影响量矩阵中只有部分的供能参数影响量是存在噪声的，即而预期影响量矩阵中其他供能参数影响量中噪声所对应的供能参数影响量ε
i,j
为0。在这种情况下，初始影响量矩阵θ
*
是低秩性，预期影响量矩阵y也是低秩的。所以，我们可以采用前文求解得到的所述预期影响量矩阵y中的部分元素来构造优化影响量矩阵，使得所述优化影响量矩阵所对应的实际供能参数更接近用户的真实供能需求。所以，供能系统采用基于该优化影响量矩阵得到的长期稳定的实际供能参数进行供能时，用户不会为了降低用能损失量而去改变自身的用能需求。
[0200]
所述优化影响量矩阵确定模块42可以根据前文求解得到的所述预期影响量矩阵y中的部分元素来构造优化影响量矩阵。
[0201]
优化影响量矩阵确定模块42用于根据所述预期影响量矩阵、以及矩阵位置集合，在最小化影响量矩阵的秩的条件下构造优化影响量矩阵；其中，所述优化影响量矩阵中包括每一用能对象在每一时段中对于所述供能系统的优化供能参数影响量；所述影响量矩阵和预期影响量矩阵是同型矩阵，且所述矩阵位置集合中包括所述预期影响量矩阵及所述影响量矩阵中部分的相应元素位置。
[0202]
在一些实施例中，利用矩阵补全算法，在最小化影响量矩阵的秩的条件下和预设误差约束下构造所述优化影响量矩阵。所述误差约束可由如下公式表示：
[0203]
||p
ω
(θ-y)||
f
＜δ
ꢀꢀꢀꢀꢀꢀ
(13)
[0204]
其中，δ为预设误差，p
ω
为投影操作。公式(13)表示的是根据矩阵位置集合对影响量矩阵θ和预期影响量矩阵y的差进行投影操作得到的矩阵的范数满足预设误差。所述投影操作p
ω
可由公式(14)来解释：
[0205][0206]
其中，ω为矩阵位置集合，所述影响量矩阵θ和预期影响量矩阵y是同型矩阵，且所述矩阵位置集合ω中包括所述预期影响量矩阵y及所述影响量矩阵θ中部分的相应元素位置。例如，矩阵位置集合ω中包括m
×
n的预期影响量矩阵y的对角线位置，由于所述影响量矩阵θ和预期影响量矩阵y是同型矩阵，所以，所述矩阵位置集合ω中还包括所述影响量矩阵θ的对角线位置。
[0207]
所述矩阵位置集合ω是通过抽样算法和预计抽样的矩阵位置个数n确定的。例如，所述抽样算法为简单随机抽样，所述预计抽样的矩阵位置个数为n，基于所述抽样概率抽取矩阵中的n个位置。其中，m与n的乘积为所述预期影响量矩阵中矩阵位置的总个数。需要说明的是，抽样的个数太多，会使构造的矩阵更接近y；抽样的个数太少，会使得无法准确构造所述优化影响量矩阵。
[0208]
在所述预设误差约束下和最小化影响量矩阵的秩的条件下，构造优化影响量矩阵的问题可以转化为公式(15)所对应的秩最小化的优化问题。
[0209][0210]
利用矩阵补全算法可以求解上述优化问题以构造所述优化影响量矩阵。其中，所述优化影响量矩阵中包括每一用能对象在每一时段中对于所述供能系统的优化供能参数影响量。
[0211]
所述矩阵补全算法包括但不限于不动点延续算法、加速近端梯度算法、半定规划法、奇异值阈值算法、增广拉格朗日乘子法等。
[0212]
所述公式(15)所对应的非凸优化问题可以转化成如下公式所示的最小值问题：
[0213][0214]
其中，||θ||
*
为所述影响量矩阵θ的奇异值之和，||.||
f
代表矩阵的f范数，矩阵的f范数是每一个矩阵位置的元素值的平方求和后的次幂，其中，m为所述用能对象的个数，n为观察期间中所述时段的个数，p为确定所述矩阵位置集合ω的抽样概率，σ是所述噪声矩阵z的标准差。所述噪声矩阵可以根据前文求解得到的所述预期影响量矩阵y与所述为初始影响量矩阵θ
*
之间的差值得到。
[0215]
利用不动点延续算法求解公式(16)的最小值问题时，需要获得如下微分方程的解：
[0216][0217]
其中，表示的是关于所述影响量矩阵θ的次梯度。不动点延续算法构造矩阵a、和b，并将b矩阵的转置与a矩阵的乘积ab
t
作为所述优化影响量矩阵。
[0218]
其中，所述矩阵a和矩阵b是基于初始矩阵a0和初始矩阵b0，利用所述不动点算法不断迭代得到的。需要说明的是，初始矩阵a0和初始矩阵b0的设置与所述不动点算法的收敛速度有关。
[0219]
在一些实施例中，所述初始矩阵a0和初始矩阵b0与前文求解得到的所述预期影响量矩阵y相关。
[0220]
具体地，对所述预期影响量矩阵y进行奇异值分解可以得到所述预期影响量矩阵y的奇异值矩阵σ
m
×
n
，奇异值矩阵σ
m
×
n
是m
×
n的对角矩阵，其中，通过简化后的奇异值矩阵可以构造所述初始矩阵a0和初始矩阵b0。
[0221]
例如，所述奇异值矩阵中的各对角元素为(σ1，σ2,...,σn)。将所述奇异值矩阵σ
m
×
n
中的各对角元素按照大小排序，得到排序后的奇异值矩阵，排序后的奇异值矩阵中的各对角元素为(σ1'，σ2',...,σn')，其中，(σ1'＞σ2'＞...＞σn')。当所述奇异值矩阵中的对角元素小于奇异值阈值时则将对应的对角元素设为0，以得到简化的奇异值矩阵所述奇异值阈值可以是所述对角线元素最大值的十分之一、百分之一、或者千分之一等，但并不以此为限。所述简化后的奇异值矩阵可由如下公式表示：
[0222][0223]
其中，k为简化后的奇异值矩阵中非0的对角元素的个数。是对角矩阵，每一对角元素为奇异值矩阵σ
m
×
n
中在奇异值阈值以上的对角元素。0
m
×
n
代表的是m
×
n维的0矩阵。
[0224]
根据简化的奇异值矩阵得到的所述初始矩阵a0和初始矩阵b0可由如下公式表示：
[0225][0226]
其中，是对所述对角矩阵中的每个元素开根号得到的。为所述优化影响量矩阵的初始矩阵。
[0227]
在其他实施例中，所述优化影响量矩阵的初始矩阵可以为m
×
n维的0矩阵，也可以是所述预期影响量矩阵y，但不以此为限。需要说明的是，在本技术中以用能对象的个数m为行并且以观察期间中所述时段的个数n为列。在其他实施方式中，也可用观察期间中所述时段的个数n为行并且以用能对象的个数m为列。
[0228]
所述供能参数确定模块43用于根据所述优化影响量矩阵，确定所述供能系统的实际供能参数，以用于控制所述供能系统按该实际供能参数对各所述用能对象供能。
[0229]
所述实际供能参数包括但不限于：每一用能对象的实际供能量、实际供能量对应的实际供能指令、或者实际供能损失参数。所述实际供能参数用于稳定能源系统中各用户的用能行为，令其合理化，且避免能源浪费。
[0230]
在一些实施例中，所述实际供能损失参数为所述优化影响量矩阵中每一用能对象各时段的优化供能参数影响量的均值。例如，所述优化影响量矩阵中一用能对象在n个时段中的优化供能参数影响量为(a1,a2,....,a
n
)，对各优化供能参数影响量求和后与n的比值作为所述稳定的实际供能损失参数。
[0231]
在一些实施例中，基于得到的所述实际供能损失参数，确定所述每一用能对象的实际供能量、或实际供能量对应的实际供能指令。例如，将与所述实际供能损失参数相对应的供能量作为长期的实际供能量，以控制不包括储能装置的供能系统按照实际供能量进行产能，以供给所述各用能对象。又如，将与所述实际供能损失参数相对应的供能量作为稳定的实际供能量后确定所述实际供能指令，以供包括储能装置的供能系统按照所述实际供能指令进行产能、储能、以及释能。其中，所述实际供能指令中包括：每一时间间隔的实际产能量、每一时间间隔的实际储能量、每一时间间隔的实际释能量。
[0232]
需要说明的是，供能系统通过所述实际供能参数进行供能时，既可以满足用户实际的用能需求又可以使得用户不会轻易改变其负荷曲线，进而改善了能源浪费的现象并且保证了电力系统的平稳运行。
[0233]
请参阅图5，图5显示为本技术的计算机设备在一实施例中的结构示意图，如图5所
示，所述计算机系统50包括存储装置51、以及处理装置52。
[0234]
所述存储装置51用于存储至少一程序，在一些实施例中，所述存储装置51包括一个或多个存储器，用以在执行时执行基于本技术技术思想而示例的用于供能系统的控制方法。在实施例中，所述存储器可包括高速随机存取存储器，并且还可包括非易失性存储器，例如一个或多个磁盘存储设备、闪存设备或其他非易失性固态存储设备。在某些实施例中，存储器还可以包括远离一个或多个处理器的存储器，例如经由rf电路或外部端口以及通信网络访问的网络附加存储器，其中所述通信网络可以是因特网、一个或多个内部网、局域网、广域网、存储局域网等，或其适当组合。存储器控制器可控制设备的诸如cpu和外设接口之类的其他组件对存储器的访问。
[0235]
所述处理装置52与所述存储装置51连接，在一些实施例中，所述处理装置52包括一个或多个处理器，所述处理装置52用于调所述至少一程序以执行本技术所述的用于供能系统的控制方法。在实施例中，所述处理器可操作地与存储器和/或非易失性存储设备耦接。更具体地，处理器可执行在存储器和/或非易失性存储设备中存储的指令以在所述计算机设备中执行操作，诸如生成实际供能参数，以用于控制所述供能系统按该实际供能参数对各所述用能对象供能。如此，处理器可包括一个或多个通用微处理器、一个或多个专用处理器、一个或多个现场可编程逻辑阵列、或它们的任何组合。
[0236]
请参阅图6，图6显示为本技术的计算机设备在另一实施例中的结构示意图，如图6所示，所述计算机设备60包括通信装置61、存储装置62、以及处理装置63。所述存储装置62、以及处理装置63与图5所示的存储装置51、以及处理装置52相同或相似，在此不再详述。
[0237]
所述通信装置61用于与供能对象通信，以输出所述实际供能参数。通过所述通信装置61，所述计算机设备61能够与多个供能对象进行数据交互或通信传输。例如，所述计算机设备将实际供能参数通过所述通信装置61发送至各供能对象，以供所述各供能对象根据所述实际供能参数支付用能费用。其中，所述通信装置61可以包括：有线通信电路和/或无线通信电路，有线通信电路如有线网卡，无线通信电路如2g-5g、wifi、314/433射频电路等等，从而能接入通信网络，通信网络可以包括因特网、移动网络、局域网(lan)、广域网(wan)、存储局域网(san)、或者一个或多个内部网等，或任何其它适当通信网络中的一个或多个。
[0238]
在一实施例中，图5或图6所示的计算机设备例如为装载有app应用程序或具备网页/网站访问性能的电子设备，例如台式电脑、智能电视、服务器、移动终端(如手机、平板电脑、笔记本电脑等)，该些电子设备还可以包含显示单元/音频单元等输入输出单元。所述显示单元可包括触摸部件，该触摸部件通过检测对象触摸其屏幕(例如，电子显示器的表面)的发生和/或位置来促进用户输入；例如，供能系统侧的操作者通过触摸所述电子设备的屏幕以输入所述奇异值阈值等。
[0239]
另外，本技术还公开一种计算机可读存储介质，存储有至少一计算机程序，所述计算机程序被执行时实现上述针对用于供能系统的控制方法所描述的至少一种实施例。
[0240]
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本技术的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个
人计算机，服务器，或者网络设备等)执行本技术各个实施例所述方法的全部或部分步骤。
[0241]
于本技术提供的实施例中，所述计算机可读写存储介质可以包括只读存储器、随机存取存储器、eeprom、cd-rom或其它光盘存储装置、磁盘存储装置或其它磁存储设备、闪存、u盘、移动硬盘、或者能够用于存储具有指令或数据结构形式的期望的程序代码并能够由计算机进行存取的任何其它介质。另外，任何连接都可以适当地称为计算机可读介质。例如，如果指令是使用同轴电缆、光纤光缆、双绞线、数字订户线(dsl)或者诸如红外线、无线电和微波之类的无线技术，从网站、服务器或其它远程源发送的，则所述同轴电缆、光纤光缆、双绞线、dsl或者诸如红外线、无线电和微波之类的无线技术包括在所述介质的定义中。然而，应当理解的是，计算机可读写存储介质和数据存储介质不包括连接、载波、信号或者其它暂时性介质，而是旨在针对于非暂时性、有形的存储介质。如申请中所使用的磁盘和光盘包括压缩光盘(cd)、激光光盘、光盘、数字多功能光盘(dvd)、软盘和蓝光光盘，其中，磁盘通常磁性地复制数据，而光盘则用激光来光学地复制数据。
[0242]
在一个或多个示例性方面，本技术所述方法的计算机程序所描述的功能可以用硬件、软件、固件或其任意组合的方式来实现。当用软件实现时，可以将这些功能作为一个或多个指令或代码存储或传送到计算机可读介质上。本技术所公开的方法或算法的步骤可以用处理器可执行软件模块来体现，其中处理器可执行软件模块可以位于有形、非临时性计算机可读写存储介质上。有形、非临时性计算机可读写存储介质可以是计算机能够存取的任何可用介质。
[0243]
本技术上述的附图中的框图，图示了按照本技术各种实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。基于此，框图中的每个方框可以代表一个模块、程序段、或代码的一部分，该模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个接连地表示的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这根据所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以通过执行规定的功能或操作的专用的基于硬件的系统来实现，或者可以通过专用硬件与计算机指令的组合来实现。
[0244]
上述实施例仅例示性说明本技术的原理及其功效，而非用于限制本技术。任何熟悉此技术的人士皆可在不违背本技术的精神及范畴下，对上述实施例进行修饰或改变。因此，举凡所属技术领域中具有通常知识者在未脱离本技术所揭示的精神与技术思想下所完成的一切等效修饰或改变，仍应由本技术的权利要求所涵盖。