本发明提供了一种基于三角综合评估的新型能源供给系统优化配置方法，其首先选取负荷失电率、负荷失热率、系统年化成本和能量过剩率作为新型能源供给系统电热供应稳定性、经济性和新能源利用情况三方面的评估指标，建立三角综合评估模型，并设定三角综合评估模型的面积为优化配置方法的目标函数；然后以三角综合评估模型的面积最大为目标，服从约束条件，建立新型能源供给系统的优化数学模型；最后采用分层迭代算法，求解出新型能源供给系统的最优配置方案，并通过实际算例验证了该优化方法的正确性和可行性，优化配置方案可以满足偏远地区电力热力的刚性需求，并有利于能源供给系统自身的经济运行和可再生能源的高效利用。
1.一种基于三角综合评估的新型能源供给系统优化配置方法，所述新型能源供给系统分为电力系统层和热力系统层，电力系统层包含光伏发电单元、变流器单元以及储电单元，热力系统层包含储热单元以及电热转换单元；其特征在于，所述新型能源供给系统优化配置方法包括如下步骤：(1)选取负荷失电率或负荷失热率作为新型能源供给系统电热供应稳定性的评估指标，选取系统年化成本作为新型能源供给系统经济性的评估指标，选取能量过剩率作为新型能源供给系统新能源利用情况的评估指标，建立三角综合评估模型并计算模型面积；
(2)以三角综合评估模型面积最大为目标，在服从运行约束条件的前提下，建立新型能源供给系统的优化数学模型如下：其中：Atri为三角综合评估模型的面积，anorm为标准化后系统的负荷失电率或负荷失热率，bnorm为标准化后的系统年化成本，cnorm为标准化后系统的能量过剩率；
(3)采用分层迭代算法遍历所有配置方案求解上述优化数学模型，得到包含各单元配置容量的系统最优配置方案，以使系统在电热供应稳定性、经济性和新能源利用情况三方面达到综合最优。
2.根据权利要求1所述的新型能源供给系统优化配置方法，其特征在于：所述负荷失电率通过以下公式计算得到：其中：LPSP为系统的负荷失电率，Psue(t)为t时刻系统的供电功率，Peload(t)为t时刻系统的电负荷功率，t为自然数，T为给定的时间长度，p(Psue(t)＜Peload(t))＝0或1，即当Psue(t)＜Peload(t)时p(Psue(t)＜Peload(t))＝1，否则p(Psue(t)＜Peload(t))＝0。
3.根据权利要求1所述的新型能源供给系统优化配置方法，其特征在于：所述负荷失热率通过以下公式计算得到：其中：LHSP为系统的负荷失热率，Psuh(t)为t时刻系统的供热功率，Phload(t)为t时刻系统的热负荷功率，t为自然数，T为给定的时间长度，p(Psuh(t)＜Phload(t))＝0或1，即当Psuh(t)＜Phload(t)时p(Psuh(t)＜Phload(t))＝1，否则p(Psuh(t)＜Phload(t))＝0。
4.根据权利要求1所述的新型能源供给系统优化配置方法，其特征在于：所述系统年化成本通过以下公式计算得到：其中：ASC为系统年化成本，Ci(j)为系统中第j个单元的初始年化投资成本，Cm(j)为系统第j个单元的年化运行维护成本，j为自然数。
5.根据权利要求1所述的新型能源供给系统优化配置方法，其特征在于：所述能量过剩率通过以下公式计算得到：其中：EXR为系统的能量过剩率，Pexr(t)为t时刻系统过剩的能量功率，Pload(t)为t时刻系统的负荷功率，t为自然数，T为给定的时间长度。
6.根据权利要求1所述的新型能源供给系统优化配置方法，其特征在于：所述步骤(1)中建立三角综合评估模型的具体实现过程如下：
1.1从坐标轴原点向不同方向延伸出三条数轴，使得这三条数轴将以原点为圆心的圆等分成三部分，各数轴之间的夹角为120°；
1.2对负荷失电率或负荷失热率、系统年化成本和能量过剩率这三个评估指标进行标准化处理；
1.3将标准化处理后的三个评估指标分别在三条数轴上标定对应的点，并将各点依次连接构成三角形，进而计算该三角形的面积。
7.根据权利要求1所述的新型能源供给系统优化配置方法，其特征在于：所述步骤(2)中的运行约束条件如下：系统功率平衡约束：
Peload(t)+Phload(t)＝Ppv(t)+Pbat(t)+Psto(t)储电单元的容量和功率交换能力约束：
SOCmin≤SOC(t)≤SOCmax 0≤|Pbat(t)|≤Pbat_max储热单元的容量和功率交换能力约束：
Qmin≤Q(t)≤Qmax 0≤|Psto(t)|≤Psto_max
SOC(t)＝SOC(t-1)+α·Pbat(t)Δt
Psto(t)＝ηetPin(t) Pout(t)＝1/ηeo·Psto(t)其中：Peload(t)为t时刻系统的电负荷功率，Phload(t)为t时刻系统的热负荷功率，PPV(t)为t时刻光伏发电单元的总输出功率，Pbat(t)为t时刻储电单元的储电或放电功率，储电时为负值，放电时为正值，Psto(t)为t时刻储热单元的储热或放热功率，储热时为正值，放热时为负值；SOC(t)和SOC(t-1)分别为t时刻和t-1时刻储电单元的荷电量，SOCmin和SOCmax分别为荷电量的下限值和上限值，Pbat_max为储电单元储电或放电功率的上限值；Q(t)为t时刻储热单元的储热量，Qmin和Qmax分别为储热量的下限值和上限值，Psto_max为储热单元储热或放热功率的上限值；t为自然数，Δt为时间间隔，ηc为储电单元的充电效率，ηd为储电单元的放电效率；Pin(t)为t时刻电热转换单元的输出功率，Pout(t)为t时刻储热单元给热负荷供热的实际输出功率，ηet为储热单元的吸热效率，ηeo为储热单元的放热效率，α为中间变量。
8.根据权利要求1所述的新型能源供给系统优化配置方法，其特征在于：所述步骤(3)中的分层迭代算法具体过程如下：
3.1获得目标地区给定时间段内各时刻系统的电负荷功率、热负荷功率和光伏发电单元的单位输出功率，进而确定光伏发电单元总功率、储电单元容量、储热单元容量、变流器单元功率和电热转换单元功率的最大可执行范围；
3.2将系统划分成电力系统层和热力系统层，并设定热力系统层所需热量由光伏发电单元产生后直接经电热转换单元转换至热力子系统，而不由储电单元供给；
3.3选取负荷失电率、系统年化成本、能源过剩率作为电力系统层的三个评估指标，选取负荷失热率、系统年化成本、能源过剩率作为热力系统层的三个评估指标；
3.4对于任一电力系统层容量配置方案，根据该方案对电力系统层进行迭代操作如下：
当PPV_e(t)＞Peload(t)情况下，判断SOC(t)与SOCmax的大小：若SOC(t)≥SOCmax，则不进行充电操作并令SOC(t)＝SOCmax；若SOC(t)＜SOCmax，则根据下式确定储电功率Pbat(t)大小并完成储电单元的充电操作；
当PPV_e(t)≤Peload(t)情况下，判断SOC(t)与SOCmin的大小：若SOC(t)≤SOCmin，则不进行放电操作并令SOC(t)＝SOCmin；若SOC(t)＞SOCmin，则根据下式确定放电功率Pbat(t)大小并完成储电单元的放电操作；
其中：Pbat(t)为t时刻储电单元的储电或放电功率，储电时为负值，放电时为正值，Peload(t)为t时刻系统的电负荷功率，SOC(t)为t时刻储电单元的荷电量，PPV_e(t)为t时刻电力系统层光伏发电单元的输出功率，SOCmin和SOCmax分别为荷电量的下限值和上限值，Pinverter为变流器单元的额定功率；
遍历求得给定时间段内各时刻的SOC(t)和Pbat(t)，计算该方案下电力系统层三个的评估指标，并建立三角综合评估模型并计算其面积；在电力系统层最大可执行范围内，按照设定步长改变电力系统层的容量配置方案，并重新根据步骤3.4进行迭代操作；
3.5对于任一热力系统层容量配置方案，根据该方案对热力系统层进行迭代操作如下：
当PPV_h(t)＞Phload(t)情况下，判断Q(t)与Qmax的大小：若Q(t)≥Qmax，则不进行储热操作并令Q(t)＝Qmax；若Q(t)＜Qmax，则根据下式确定储热功率Psto(t)大小并完成储热单元的吸热操作；
当PPV_h(t)≤Phload(t)情况下，判断Q(t)与Qmin的大小：若Q(t)≤Qmin，则不进行放热操作并令Q(t)＝Qmin；若Q(t)＞Qmin，则根据下式确定放热功率Psto(t)大小并完成储热单元的放热操作；
其中：Psto(t)为t时刻储热单元的储热或放热功率，储热时为正值，放热时为负值，Phload(t)为t时刻系统的热负荷功率，Q(t)为t时刻储热单元的储热量，PPV_h(t)为t时刻热力系统层光伏发电单元经电热转换单元的输出功率，Qmin和Qmax分别为储热量的下限值和上限值，Pwater_cycle为水循环装置的额定功率；
遍历求得给定时间段内各时刻的Q(t)和Psto(t)，计算该方案下热力系统层三个的评估指标，并建立三角综合评估模型并计算其面积；在热力系统层最大可执行范围内，按照设定步长改变热力系统层的容量配置方案，并重新根据步骤3.5进行迭代操作；
3.6确定电力系统层和热力系统层的最优容量配置方案，即得到包含各单元配置容量的系统最优配置方案，以使系统在电热供应稳定性、经济性和新能源利用情况三方面达到综合最优。
一种基于三角综合评估的新型能源供给系统优化配置方法
技术领域
[0001] 本发明属于能源供给技术领域，具体涉及一种基于三角综合评估的新型能源供给系统优化配置方法。
背景技术
[0002] 目前，我国部分偏远地区由于自然环境恶劣、居民居住分散，存在用电难、用热难的问题；考虑到实际供电工程建设成本极高，基于可再生能源的开发利用，以解决偏远地区的供电供热问题，则成为现阶段能源供给领域的研究热点。因此，新型能源供给系统应运而生，新型能源供给系统以光伏发电为主，综合利用新能源分布式发电、储能、储热的先进技术，构成光、储、热一体化的能源供应，可以实现对分布式可再生能源的充分消纳以及电力、热力等多种形式的能量供给，有助于满足偏远地区居民用户的用电供热需求。
[0003] 新型能源供给系统的结构如图1所示，由光伏发电单元、变流器单元、储电单元、储热单元及电热转换单元组成。其中，光伏发电单元、变流器单元和储电单元共同构成了新型能源供给系统的电力系统层，主要实现太阳能到电能的能量转换，用于为偏远地区的交流电负荷提供电力供应；储热单元和电热转换单元共同构成了新型能源供给系统的热力系统层，主要实现电能到热能的能量转换，用于为偏远地区的热负荷提供热力供应。此外，在热力系统层中，通常需要采用额外的水循环装置，用于实现储热单元、电热转换单元和热负荷三者之间的热能传输。
[0004] 为了实现新型能源供给系统的合理运行，需要对所含单元进行合理的容量配置。现有针对交(直)流微网系统的优化配置方法，一般只涉及到电力供配的情况，根据交(直)流微网系统对应的负荷数据，选取一或多个性能指标以及合适的数值优化算法，进而优化得到在该负荷数据下性能指标较优的优化配置结果。其中，性能指标的选取需要对应交(直)流微网系统的关键特性，通常从经济性和稳定性两方面考虑；常用的数值优化算法包括迭代算法、神经网络算法、遗传算法等先进算法，通常根据优化问题的计算难度和收敛精度来选择性能匹配的优化算法。而新型能源供给系统涉及到光伏、电力和热力等多种能量形式的供配调控，极大地复杂了系统配置的优化过程；因此，针对新型能源供给系统，采用多目标优化算法有利于实现系统性能的综合最优。由于新型能源供给系统具有多种能量形式，且需要对偏远地区提供电力热力供应，性能指标的选择需要综合考虑电力系统层和热力系统层的可靠性情况、整个能源供给系统的经济性情况以及对于太阳能的利用情况。此外，现有的多目标算法一般是将多个目标进行线性加和或基于模糊隶属度函数将多目标问题转变为非线性单目标优化问题；但这两种方法中，线性加权系数依赖经验设计，具有极大的不确定性，而模糊隶属度函数的选取则十分复杂，如何针对新型能源供给系统的多个性能指标进行合理的综合综合评估，则决定了优化算法的执行过程和优化结果的满意程度。
发明内容
[0005] 鉴于上述，本发明提供了一种基于三角综合评估的新型能源供给系统优化配置方法，以实现新型能源供给系统在电热供应稳定性、经济性和新能源利用情况三方面综合最优为目标，优化配置方案可以满足偏远地区电力热力的刚性需求，并有利于能源供给系统自身的经济运行和可再生能源的高效利用。
[0006] 一种基于三角综合评估的新型能源供给系统优化配置方法，所述新型能源供给系统分为电力系统层和热力系统层，电力系统层包含光伏发电单元、变流器单元以及储电单元，热力系统层包含储热单元以及电热转换单元；所述新型能源供给系统优化配置方法包括如下步骤：
[0007] (1)选取负荷失电率或负荷失热率、系统年化成本和能量过剩率分别作为新型能源供给系统电热供应稳定性、经济性和新能源利用情况三方面的评估指标，建立三角综合评估模型并计算模型面积；
[0008] (2)以三角综合评估模型面积最大为目标，在服从运行约束条件的前提下，建立新型能源供给系统的优化数学模型如下：
[0009]
[0010] 其中：Atri为三角综合评估模型的面积，anorm为标准化后系统的负荷失电率或负荷失热率，bnorm为标准化后的系统年化成本，cnorm为标准化后系统的能量过剩率；
[0011] (3)采用分层迭代算法遍历所有配置方案求解上述优化数学模型，得到包含各单元配置容量的系统最优配置方案，以使系统在电热供应稳定性、经济性和新能源利用情况三方面达到综合最优。
[0012] 进一步地，所述负荷失电率通过以下公式计算得到：
[0013]
[0014] 其中：LPSP为系统的负荷失电率，Psue(t)为t时刻系统的供电功率，Peload(t)为t时刻系统的电负荷功率，t为自然数，T为给定的时间长度，p(Psue(t)＜Peload(t))＝0或1，即当Psue(t)＜Peload(t)时p(Psue(t)＜Peload(t))＝1，否则p(Psue(t)＜Peload(t))＝0。
[0015] 进一步地，所述负荷失热率通过以下公式计算得到：
[0016]
[0017] 其中：LHSP为系统的负荷失热率，Psuh(t)为t时刻系统的供热功率，Phload(t)为t时刻系统的热负荷功率，t为自然数，T为给定的时间长度，p(Psuh(t)＜Phload(t))＝0或1，即当Psuh(t)＜Phload(t)时p(Psuh(t)＜Phload(t))＝1，否则p(Psuh(t)＜Phload(t))＝0。
[0018] 进一步地，所述系统年化成本通过以下公式计算得到：
[0019]
[0020] 其中：ASC为系统年化成本，Ci(j)为系统中第j个单元的初始年化投资成本，Cm(j)为系统第j个单元的年化运行维护成本，j为自然数。
[0021] 进一步地，所述能量过剩率通过以下公式计算得到：
[0022]
[0023] 其中：EXR为系统的能量过剩率，Pexr(t)为t时刻系统过剩的能量功率，Pload(t)为t时刻系统的负荷功率，t为自然数，T为给定的时间长度。
[0024] 进一步地，所述步骤(1)中建立三角综合评估模型的具体实现过程如下：
[0025] 1.1从坐标轴原点向不同方向延伸出三条数轴，使得这三条数轴将以原点为圆心的圆等分成三部分，各数轴之间的夹角为120°；
[0026] 1.2对负荷失电率或负荷失热率、系统年化成本和能量过剩率这三个评估指标进行标准化处理；
[0027] 1.3将标准化处理后的三个评估指标分别在三条数轴上标定对应的点，并将各点依次连接构成三角形，进而计算该三角形的面积。
[0028] 进一步地，所述步骤(2)中的运行约束条件如下：
[0029] 系统功率平衡约束：
[0030] Peload(t)+Phload(t)＝Ppv(t)+Pbat(t)+Psto(t)
[0031] 储电单元的容量和功率交换能力约束：
[0032] SOCmin≤SOC(t)≤SOCmax 0≤|Pbat(t)|≤Pbat_max
[0033] 储热单元的容量和功率交换能力约束：
[0034] Qmin≤Q(t)≤Qmax 0≤|Psto(t)|≤Psto_max
[0035] 储电单元的效率约束：
[0036] SOC(t)＝SOC(t-1)+α·Pbat(t)Δt
[0037]
[0038] 储热单元的效率约束：
[0039] Psto(t)＝ηetPin(t) Pout(t)＝1/ηeo·Psto(t)
[0040] 其中：Peload(t)为t时刻系统的电负荷功率，Phload(t)为t时刻系统的热负荷功率，PPV(t)为t时刻光伏发电单元的总输出功率，Pbat(t)为t时刻储电单元的储电或放电功率，储电时为负值，放电时为正值，Psto(t)为t时刻储热单元的储热或放热功率，储热时为正值，放热时为负值；SOC(t)和SOC(t-1)分别为t时刻和t-1时刻储电单元的荷电量，SOCmin和SOCmax分别为荷电量的下限值和上限值，Pbat_max为储电单元储电或放电功率的上限值；Q(t)为t时刻储热单元的储热量，Qmin和Qmax分别为储热量的下限值和上限值，Psto_max为储热单元储热或放热功率的上限值；t为自然数，Δt为时间间隔，ηc为储电单元的充电效率，ηd为储电单元的放电效率；Pin(t)为t时刻电热转换单元的输出功率，Pout(t)为t时刻储热单元给热负荷供热的实际输出功率，ηet为储热单元的吸热效率，ηeo为储热单元的放热效率。
[0041] 进一步地，所述步骤(3)中的分层迭代算法具体过程如下：
[0042] 3.1获得目标地区给定时间段内各时刻系统的电负荷功率、热负荷功率和光伏发电单元的单位输出功率，进而确定光伏发电单元总功率、储电单元容量、储热单元容量、变流器单元功率和电热转换单元功率的最大可执行范围；
[0043] 3.2将系统划分成电力系统层和热力系统层，并设定热力系统层所需热量由光伏发电单元产生后直接经电热转换单元转换至热力子系统，而不由储电单元供给；
[0044] 3.3选取负荷失电率、系统年化成本、能源过剩率作为电力系统层的三个评估指标，选取负荷失热率、系统年化成本、能源过剩率作为热力系统层的三个评估指标；
[0045] 3.4对于任一电力系统层容量配置方案，根据该方案对电力系统层进行迭代操作如下：
[0046] 当PPV_e(t)＞Peload(t)情况下，判断SOC(t)与SOCmax的大小：若SOC(t)≥SOCmax，则不进行充电操作并令SOC(t)＝SOCmax；若SOC(t)＜SOCmax，则根据下式确定储电功率Pbat(t)大小并完成储电单元的充电操作；
[0047]
[0048] 当PPV_e(t)≤Peload(t)情况下，判断SOC(t)与SOCmin的大小：若SOC(t)≤SOCmin，则不进行放电操作并令SOC(t)＝SOCmin；若SOC(t)＞SOCmin，则根据下式确定放电功率Pbat(t)大小并完成储电单元的放电操作；
[0049]
[0050] 其中：Pbat(t)为t时刻储电单元的储电或放电功率，储电时为负值，放电时为正值，Peload(t)为t时刻系统的电负荷功率，SOC(t)为t时刻储电单元的荷电量，PPV_e(t)为t时刻电力系统层光伏发电单元的输出功率，SOCmin和SOCmax分别为荷电量的下限值和上限值，Pinverter为变流器单元的额定功率；
[0051] 遍历求得给定时间段内各时刻的SOC(t)和Pbat(t)，计算该方案下电力系统层三个的评估指标，并建立三角综合评估模型并计算其面积；在电力系统层最大可执行范围内，按照设定步长改变电力系统层的容量配置方案，并重新根据步骤3.4进行迭代操作；
[0052] 3.5对于任一热力系统层容量配置方案，根据该方案对热力系统层进行迭代操作如下：
[0053] 当PPV_h(t)＞Phload(t)情况下，判断Q(t)与Qmax的大小：若Q(t)≥Qmax，则不进行储热操作并令Q(t)＝Qmax；若Q(t)＜Qmax，则根据下式确定储热功率Psto(t)大小并完成储热单元的吸热操作；
[0054]
[0055] 当PPV_h(t)≤Phload(t)情况下，判断Q(t)与Qmin的大小：若Q(t)≤Qmin，则不进行放热操作并令Q(t)＝Qmin；若Q(t)＞Qmin，则根据下式确定放热功率Psto(t)大小并完成储热单元的放热操作；
[0056]
[0057] 其中：Psto(t)为t时刻储热单元的储热或放热功率，储热时为正值，放热时为负值，Phload(t)为t时刻系统的热负荷功率，Q(t)为t时刻储热单元的储热量，PPV_h(t)为t时刻热力系统层光伏发电单元经电热转换单元的输出功率，Qmin和Qmax分别为储热量的下限值和上限值，Pwater_cycle为水循环装置的额定功率；
[0058] 遍历求得给定时间段内各时刻的Q(t)和Psto(t)，计算该方案下热力系统层三个的评估指标，并建立三角综合评估模型并计算其面积；在热力系统层最大可执行范围内，按照设定步长改变热力系统层的容量配置方案，并重新根据步骤3.5进行迭代操作；
[0059] 3.6确定电力系统层和热力系统层的最优容量配置方案，即得到包含各单元配置容量的系统最优配置方案，以使系统在电热供应稳定性、经济性和新能源利用情况三方面达到综合最优。
[0060] 基于上述技术方案，本发明方法具有以下有益技术效果：
[0061] 1.本发明选取负荷失电率、负荷失热率、系统年化成本和能量过剩率等多个评估指标，建立三角模型进行综合评估，有利于实现新型能源供给系统在电热供应稳定性、经济性和新能源利用情况三方面的综合最优。
[0062] 2.本发明建立三角综合评估模型，并设定三角综合评估模型的面积作为优化目标函数，实现多目标的平衡化处理，避免了传统方法将多个指标加和作为优化目标的加权系数设定问题。
[0063] 3.本发明采用分层迭代方法分别对电力系统层和热力系统层进行迭代处理，简化了整个新型能源供给系统的优化过程。
附图说明
[0064] 图1为新型能源供给系统的结构示意图。
[0065] 图2为某偏远地区电负荷、夏季典型日和冬季典型日1MW光伏发电单元输出功率和热负荷的预测数据示意图。
[0066] 图3为新型能源供给系统电力系统层最优容量配置下各单元的功率(容量)曲线示意图。
[0067] 图4为新型能源供给系统热力系统层最优容量配置下各单元的功率(容量)曲线示意图。
[0068] 图5为新型能源供给系统最优容量配置的成本分布示意图。
具体实施方式
[0069] 为了更为具体地描述本发明，下面结合附图及具体实施方式对本发明的技术方案进行详细说明。
[0070] 本发明基于三角综合评估的新型能源供给系统优化配置方法，包括如下步骤：
[0071] (1)选取负荷失电率、负荷失热率、系统年化成本和能量过剩率作为新型能源供给系统电热供应稳定性、经济性和新能源利用情况三方面的评估指标，建立三角综合评估模型，并设定三角综合评估模型的面积为所述优化配置方法的目标函数。
[0072] 负荷失电率通过以下公式计算得到：
[0073]
[0074] 其中：Psue(t)为t时刻系统的供电功率，Peload(t)为t时刻的电负荷功率，T为总时段数，LPSP为系统的负荷失电率，等于供电功率小于电负荷功率的时间概率。
[0075] 负荷失热率通过以下公式计算得到：
[0076]
[0077] 其中：Psuh(t)为t时刻系统的供热功率，Phload(t)为t时刻的热负荷功率，T为总时段数，LHSP为系统的负荷失热率，等于供热功率小于热负荷功率的时间概率。
[0078] 系统年化成本通过以下公式计算得到：
[0079]
[0080] 其中：ASC为系统年化成本，Ci(j)为系统第j个单元的初始年化投资成本，Cm(j)为系统第j个单元的年化运行维护成本。
[0081] 能源过剩率通过以下公式计算得到：
[0082]
[0083] 其中：EXR为系统的能量过剩率，Pexr(t)为t时刻系统过剩的能量功率，Pload(t)为t时刻系统的负荷功率。
[0084] 建立三角综合评估模型包括如下步骤：
[0085] 1.1以原点为圆心延伸出3条数轴，将圆等分为3个部分，各数轴之间的夹角为120°。
[0086] 1.2选取3个评估指标，并通过以下公式进行标准化处理：
[0087]
[0088] 其中：LPSPdesire为期望的负荷失电率，max(0,LPSP-LPSPdesire)等于0和LPSP-LPSPdesire中较大的值，LPSPmax为max(0,LPSP-LPSPdesire)可能的最大值，LPSPmin为max(0,LPSP-LPSPdesire)可能的最小值。
[0089]
[0090] 其中：LHSPdesire为期望的负荷失热率，max(0,LHSP-LHSPdesire)等于0和LHSP-LHSPdesire中较大的值，LHSPmax为max(0,LHSP-LHSPdesire)可能的最大值，LHSPmin为max(0,LHSP-LHSPdesire)可能的最小值。
[0091]
[0092] 其中：ASCmax为ASC可能的最大值，ASCmin为ASC可能的最小值。
[0093]
[0094] 其中：EXRmax为EXR可能的最大值，EXRmin为EXR可能的最小值。
[0095] 1.3将标准化处理后的三个指标依次在各数轴上标定对应的点，并将各点依次连接，构成一个三角形。
[0096] 1.4计算三角形的面积，作为目标函数，目标函数通过以下公式计算得到：
[0097]
[0098] 其中：Atri为三角形的面积，anorm、bnorm、cnorm为从LPSPnorm、LHSPnorm、ASCnorm、EXRnorm中选取的三个评估指标。
[0099] (2)以三角综合评估模型的面积最大为目标，服从以下约束条件，建立新型能源供给系统的优化数学模型。
[0100] 系统功率平衡约束：
[0101] Peload(t)+Phload(t)＝Ppv(t)+Pbat(t)+Psto(t) (10)
[0102] 其中：PPV(t)为t时刻光伏发电单元的输出功率，Pbat(t)为t时刻储电单元的储电或发电功率，储电单元储电时为负值，发电时为正值，Psto(t)为t时刻储热单元的储热或放热功率，储热单元储热时为正值，发热时为负值。
[0103] 储电单元的容量和功率交换能力约束：
[0104] SOCmin≤SOC(t)≤SOCmax (11)
[0105] 其中：SOC(t)为t时刻储电单元的荷电量，SOCmin为最小荷电量，SOCmax为最大荷电量。
[0106] 0≤|Pbat(t)|≤Pbat_max (12)
[0107] 其中：Pbat_max为储电单元的最大储电或放电功率，受到变流器单元功率的Pinverter限制。
[0108] 储热单元的容量和功率交换能力约束：
[0109] Qmin≤Q(t)≤Qmax (13)
[0110] 其中：Q(t)为t时刻储热单元的储热量，Qmin为最小储热量，Qmax为最大储热量。
[0111] 0≤|Psto(t)|≤Psto_max (14)
[0112] 其中：Psto_max为储热单元的最大储热或放热功率。
[0113] 储电单元的效率约束：
[0114] SOC(t)＝SOC(t-1)+α·Pbat(t)Δt (15)
[0115] 其中：SOC(t-1)为t-1时刻储电单元的荷电量，Δt为时间间隔，ηc为储电单元的充电效率，ηd为储电单元的放电效率。
[0116] 储热单元的效率约束：
[0117] Psto(t)＝ηetPin(t) (16)
[0118] Pout(t)＝1/ηeo·Psto(t) (17)
[0119] 式中：Pin(t)为电热转换单元的输出功率，Pout(t)为储热单元给热负荷供热的实际输出功率,ηet为储热单元的吸热效率，ηeo为储热单元的放热效率。
[0120] 新型能源供给系统的优化数学模型如下：
[0121]
[0122] (3)采用分层迭代算法，求解出新型能源供给系统的最优配置方案，该方案包含了各单元的配置容量，以使得系统在电热供应稳定性、经济性和新能源利用情况三方面达到综合最优；分层迭代算法包括如下步骤：
[0123] 3.1获得目标地区典型m天1～24h各时刻的电负荷数据Peload(t)、热负荷数据Phload(t)和单位光伏单元的输出功率PPV_one(t)，由PPV_one(t)、Peload(t)、Phload(t)，确定光伏发电单元总功率、储电单元容量、储热单元容量、变流器单元功率和电热转换单元功率的最大可执行范围。
[0124] 3.2将新型能源供给系统划分成电力系统层和热力系统层，并设定热力系统层所需热量是由光伏发电单元产生后，直接经电热转换单元转换至热力子系统，而不由储电单元供给。
[0125] 3.3选取负荷失电率、系统年化成本、能源过剩率作为电力系统层的三个评估指标，选取负荷失热率、系统年化成本、能源过剩率作为热力系统层的三个评估指标。
[0126] 3.4针对某一个电力系统层容量配置方案，对电力系统层进行迭代操作如下：
[0127] 当t时刻PPV_e(t)>Peload(t)时，判断储电单元荷电量SOC(t)是否小于储电单元总额定容量即SOC(t)
[0128]
[0129] 其中：PPV_e(t)为电力系统层光伏发电单元的输出功率。
[0130] 当t时刻PPV_e(t)≤Peload(t)时，判断储电单元荷电量SOC(t)是否大于SOCmin即SOC(t)>SOCmin。若SOC(t)小于SOCmin，则不进行放电操作，并令SOC(t)＝SOCmin；若SOC(t)大于SOCmin，则根据下式确定放电功率大小，并完成储电单元的放电操作。
[0131]
[0132] 依次求得m天1～24h各时刻的SOC(t)和Pbat(t)，计算该容量配置方案下电力系统层的评估指标，并建立三角评估模型，计算对应的目标函数数值。
[0133] 在电力系统层的最大执行范围内，按照设定步长，改变电力系统层的容量配置，并重新进行上述的迭代操作；其中，步长由电力系统层包含的各单元单位功率或容量决定。
[0134] 3.5针对某一个热力系统层的容量配置情方案，对热力系统层进行迭代操作如下：
[0135] 当t时刻光伏发电单元经电热转换单元的输出功率PPV_h(t)>Phload(t)，判断储热单元储热量Q(t)是否小于Qmax即Q(t)
[0136]
[0137] 其中：PPV_h(t)为热力系统层光伏发电单元经电热转换单元的输出功率，Pwater_cycle为水循环装置的额定功率。
[0138] 当t时刻PPV_h(t)≤Phload(t)时，判断Q(t)是否大于Qmin即Q(t)>Qmin。若Q(t)小于或等于Qmin，则不进行储热操作，并令Q(t)＝Qmin；若Q(t)大于Qmin，则根据下式确定储热单元放热功率，并完成放热操作。
[0139]
[0140] 依次求得m天1～24h各时刻Q(t)和Psto(t)，计算该容量配置方案下热力系统层的评估指标，并建立三角评估模型，计算对应的目标函数数值。
[0141] 在热力系统层的最大执行范围内，按照设定步长，改变热力系统层的容量配置，并重新进行上述迭代操作。
[0142] 3.6确定电力系统层和热力系统层最优的容量配置方案，即得到新型能源供给系统的配置方案，该方案可使得系统在电热供应稳定性、经济性和新能源利用情况三方面达到综合最优。
[0143] 下面根据实际算例对本发明技术效果进行说明。
[0144] 新型能源供给系统典型结构如图1所示，包括光伏发电单元、储电单元、储热单元、变流器单元和电热转换单元。其中，储电单元采用Li蓄电池，电热转换单元采用量子电锅炉，能效转换比为2，储热单元需配备水循环单元来完成吸热和放热。图2给出了某偏远地区电负荷、夏季典型日和冬季典型日1MW光伏发电单元输出功率和热负荷的预测数据；其中，电负荷最大功率为746kW，夏季典型日热负荷最大功率为564kW，冬季典型日热负荷最大功率为1128kW。表1给出了新型能源供给系统各单元的成本参数；其中，储电单元最大荷电量为额定荷电量的0.9，最小荷电量为额定荷电量的0.3，储热单元最大储热量为额定储热量的0.9，最小储热量为额定储热量的0.4。设定储电单元的充电效率ηc和放电效率ηd都为1，储热单元的吸热效率ηet和放热效率ηeo都为1，优化周期为1年，分层迭代步长为1小时，即总时段数T为8760。
[0145] 表1
[0146]
[0147]
[0148] 设定期望的负荷失电率LPSPdesire为0.01，使用期限为20年，由图2所示的电负荷和1MW光伏发电单元输出功率的预测数据，经分层迭代算法得到电力系统层的最优容量配置方案为：光伏发电单元2.3MW，变流器单元1.0MW，储电单元6890组。该方案下，电力系统层的负荷失电率LPSP为0，年化成本ASC为170.2万元/年，能量过剩率EXR单日最大为0.3360，最小为0.0496；图3为新型能源供给系统电力系统层最优容量配置下各单元的功率(容量)曲线。
[0149] 设定期望的负荷失热率LHSPdesire为0.01，使用期限为20年，由图2所示的热负荷和1MW光伏发电单元输出功率的预测数据，经分层迭代算法得到热力系统层的最优容量配置方案为：光伏发电单元2.1MW，电热转换单元(电锅炉)2.1MW，储热单元水循环系统20套，储热单元22574.9kWh。该方案下，热力系统层的负荷失热率LHSP为0，年华成本为315.4万元/年，能量过剩率单日最大为1.0691，最小为0.04257；图4为新型能源供给系统热力系统层最优容量配置下各单元的功率(容量)曲线。
[0150] 由此确定，针对图2的预测数据，新型能源供给系统综合最优的容量配置方案如下：光伏发电单元4.4MW，变流器单元1.0MW，储电单元6890组，电热转换单元(电锅炉)2.1MW，储热单元水循环系统20套，储热单元22574.9kWh，总年华成本为485.6万元/年；图5为新型能源供给系统最优容量配置的成本分布图。
[0151] 上述对实施例的描述是为便于本技术领域的普通技术人员能理解和应用本发明。熟悉本领域技术的人员显然可以容易地对上述实施例做出各种修改，并把在此说明的一般原理应用到其他实施例中而不必经过创造性的劳动。因此，本发明不限于上述实施例，本领域技术人员根据本发明的揭示，对于本发明做出的改进和修改都应该在本发明的保护范围之内。
