CN110323740A - 计及光伏出力和空调负荷动态相关性的经济调峰方法 - Google Patents
计及光伏出力和空调负荷动态相关性的经济调峰方法 Download PDFInfo
- Publication number
- CN110323740A CN110323740A CN201910551461.0A CN201910551461A CN110323740A CN 110323740 A CN110323740 A CN 110323740A CN 201910551461 A CN201910551461 A CN 201910551461A CN 110323740 A CN110323740 A CN 110323740A
- Authority
- CN
- China
- Prior art keywords
- user
- air conditioner
- air
- photovoltaic
- power
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000004378 air conditioning Methods 0.000 claims abstract description 76
- 238000005315 distribution function Methods 0.000 claims abstract description 52
- 230000001186 cumulative effect Effects 0.000 claims abstract description 38
- 230000005611 electricity Effects 0.000 claims abstract description 31
- 230000005855 radiation Effects 0.000 claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims description 48
- 241000039077 Copula Species 0.000 claims description 24
- 230000007704 transition Effects 0.000 claims description 18
- 238000012360 testing method Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 240000002853 Nelumbo nucifera Species 0.000 claims description 6
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 6
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 6
- 238000005265 energy consumption Methods 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 238000005057 refrigeration Methods 0.000 description 8
- 238000013480 data collection Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 4
- 241000196324 Embryophyta Species 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 101000604223 Homo sapiens Nocturnin Proteins 0.000 description 1
- 102100038815 Nocturnin Human genes 0.000 description 1
- 230000002301 combined effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/008—Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
-
- H02J3/383—
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2310/00—The network for supplying or distributing electric power characterised by its spatial reach or by the load
- H02J2310/10—The network having a local or delimited stationary reach
- H02J2310/12—The local stationary network supplying a household or a building
- H02J2310/14—The load or loads being home appliances
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B10/00—Integration of renewable energy sources in buildings
- Y02B10/10—Photovoltaic [PV]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Air Conditioning Control Device (AREA)
Abstract
本发明公开了一种计及光伏出力和空调负荷动态相关性的经济调峰方法,包括:获取用户空调用电行为的控制温度和舒适温度区间、预设期间的室外环境温度和太阳光辐射度;计算光伏电源的电池温度、光伏出力情况、空调负荷有功功率、光伏出力情况和空调负荷有功功率的边缘累积概率分布函数及光伏有功出力样本值和用户空调负荷功率样本值,建立目标函数,由目标函数和约束条件组成规划问题,并进行迭代求解,获取目标函数的最小值为结果,根据结果得到用户空调用电行为模式调度方案。该方法充分考虑了光伏分布式电源和空调系统运行方式之间的动态相关性、调峰幅度大、经济性高和能源综合利用效率高。
Description
技术领域
本发明涉及新能源综合利用技术领域,特别涉及一种计及光伏出力和空调负荷动态相关性的经济调峰方法。
背景技术
气象系统是一个复杂的系统,光照、湿度、温度与风速等气象变量间具有一定的互异性、相似性、相关性和耦合性等。气象变量之间的相关性决定了分布式光伏电源和冷负荷之间的相关性。随着配电网光伏供电技术的快速发展,光伏电站已具备一定的无功输出和调节能力,但同时也带来了电压波动、越限和损耗增加等电能质量问题。对于光伏电站的电力系统,空调负荷作为柔性负荷存在着电力系统调控的巨大潜力,是优质的需求响应调控资源。现有的对于含光伏电站和空调负荷的同一地区电力系统的综合资源调度,仍然是分别只考虑光伏发电出力的参数和只考虑空调冷负荷的参数,这种方式没有考虑光伏出力和空调冷负荷之间的动态相关性对电力系统稳定运行的影响,从而导致最终的调度方案准确度和能源利用效率较低以及经济成本较高。
发明内容
本发明旨在至少在一定程度上解决相关技术中的技术问题之一。
为此,本发明的目的在于提出一种计及光伏出力和空调负荷动态相关性的经济调峰方法,该方法可推广至更多分布式电源和不同类型负荷的配电系统,有效促进能源综合配置和协调规划,提高分布式电源的渗透率和发电效率。
为达到上述目的,本发明实施例提出了一种计及光伏出力和空调负荷动态相关性的经济调峰方法,包括:
S1,获取用户空调用电行为的控制温度和舒适温度区间、预设期间的室外环境温度和太阳光辐射度;
S2,根据所述预设期间的室外环境温度和太阳光辐射度计算光伏电源的电池温度;
S3,获取电力系统的光伏电源额定功率,根据所述太阳光辐射度、所述光伏电源的电池温度和所述电力系统的光伏电源额定功率求得光伏出力情况;
S4,根据所述用户空调用电行为的控制温度得到用户空调系统的开关状态,根据所述预设期间的室外环境温度和所述用户空调系统的开关状态计算空调负荷有功功率;
S5,计算所述光伏出力情况和所述空调负荷有功功率的边缘累积概率分布函数,再根据所述边缘累积概率分布函数计算Copula函数求得联合概率分布函数,并转换为矩阵形式以得到初步多维变量矩阵;
S6,利用Markov转移概率矩阵对所述初步多维变量矩阵进行时序重构计算新的多维变量矩阵以得到新的联合概率分布函数;
S7,根据所述新的联合概率分布函数随机生成预设数量的所述光伏出力情况和所述空调负荷有功功率的边缘累积概率分布值,将所述光伏出力情况和所述空调负荷有功功率的边缘累积概率分布值代入对应的所述边缘累积概率分布函数得到光伏有功出力样本值和用户空调负荷功率样本值;
S8,建立目标函数,由所述目标函数和约束条件组成规划问题,将所述光伏有功出力样本值作为电力系统中首节点的平衡功率输入,所述用户空调负荷功率样本值作为电力系统中PV节点的注入功率,根据所述舒适温度区间、所述光伏有功出力样本值和所述用户空调负荷功率样本值对所述规划问题进行迭代求解,获取所述目标函数的最小值为结果,根据所述结果得到用户空调用电行为模式调度方案。
本发明实施例的计及光伏出力和空调负荷动态相关性的经济调峰方法,通过将Copula函数多维马尔科夫链的转移概率矩阵相结合,获取了分布式光伏电源系统中的有功出力和用户空调负荷功率分布两部分变量的动态联合概率分布,继而得到两部分变量的边缘累积概率分布样本,作为电力系统节点的输入,将该方法与分布式光伏电源系统和空调冷负荷系统的传统分开计算方法进行比较,充分考虑了分布式光伏电源系统和空调冷负荷系统之间的动态相关性和时序性,因而可以准确描绘出两个系统通过耦合对源-荷的联合影响,大大提高电力系统随机生产仿真的精度。
另外,根据本发明上述实施例的计及光伏出力和空调负荷动态相关性的经济调峰方法还可以具有以下附加的技术特征:
进一步地,在本发明的一个实施例中,所述根据所述用户空调用电行为的控制温度得到用户空调系统的开关状态,包括:
若用户室内温度大于所述用户空调用电行为的控制温度,则所述用户空调系统为开启状态;
若用户室内温度小于所述用户空调用电行为的控制温度,则所述用户空调系统的为关闭状态。
进一步地,在本发明的一个实施例中,所述空调负荷有功功率的计算公式为:
其中,PAC为空调负荷有功功率,SAC为用户空调系统的开关状态,1表示空调开启,0表示空调关闭;COP为空调能耗比,QAC为等效热导率。
进一步地,在本发明的一个实施例中,所述S5包括:
S51,利用矩估计方法分别对所述光伏出力情况和所述空调负荷功率计算所述边缘累积概率分布函数;
S52,对步骤S51中得到的所述边缘累积概率分布函数的分布类型和参数性质进行t检验,并判断假设检验是否通过,若是则进入步骤S53,若否则返回步骤S51;
S53,根据所述光伏出力情况和所述空调负荷功率的边缘累积概率分布函数,计算所述光伏出力情况和所述用户空调负荷有功功率的Copula函数表达式作为所述联合概率分布函数,并转换为矩阵形式以得到所述初步多维变量矩阵。
进一步地,在本发明的一个实施例中,计算所述光伏出力情况和所述用户空调负荷有功功率的Copula函数表达式类型包括:正态Copula函数、t-Copula函数、Gumbel-Copula、Clayton-Copula和Frank-Copula函数,并采用欧式平方距离最小化原则确定欧式平方距离中最小值所对应的Copula函数以作为所述联合概率分布函数。
进一步地,在本发明的一个实施例中,所述S6具体包括:
S61,在所述初步多维变量矩阵任意选取一组列向量作为所述新的多维变量矩阵的初始列,继而计算下一时刻多维变量的状态转移概率,依据最大状态转移概率原则选取下一时刻变量所属区域;
S62,对所述初步多维变量矩阵中从初始列开始按列依次判断所属区域,并寻找符合所属区域的多维列向量作为所述新的多维变量矩阵的下一列,得到所述新的多维变量矩阵,继而得到所述新的联合概率分布函数。
进一步地,在本发明的一个实施例中,所述Markov转移概率矩阵为:
Pij=P(Xn+1=si|Xn=sj),si,sj∈s
其中,Pij为状态转移概率;S为变量状态序列。
进一步地,在本发明的一个实施例中,所述目标函数为:
f=max|Pmax-Pmax’|
其中,Pmax,Pmax’分别为未计及和已计及光伏出力和空调负荷相关性的有功功率峰值负荷。
进一步地,在本发明的一个实施例中,所述约束条件包括:
潮流方程等式约束:
P1=PPV’
Pi=PAC’,i∈{2,…,m}
其中,Ui、Pi和Qi分别为每个节点的电压、有功和无功功率;Gij和Bij为每个节点的电导和电纳;n表示电力系统的节点个数;m表示电力系统的PV节点个数;P1为首节点的有功功率,PPV’为光伏有功出力样本值,PAC’为空调负荷有功功率样本值,ΔPi为有功功率不平衡量,δij为节点i和节点j电压的相位差,ΔQi为无功功率不平衡量;
控制变量不等式约束:
为室内温度最小值,Tindoor为室内温度,为室内温度最大值;
电压上下限约束:
为节点i电压最小值,为节点i电压最大值,β为节点电压置信水平。
进一步地,在本发明的一个实施例中,所述目标函数为最大化电力系统中峰值负荷削减量,所述目标函数的控制变量为用户室内温度,状态变量为节点电压。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明一个实施例的计及光伏出力和空调负荷动态相关性的经济调峰方法流程图;
图2为根据本发明一个实施例的计及光伏出力和空调负荷动态相关性的经济调峰方法流程框图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
下面参照附图描述根据本发明实施例提出的计及光伏出力和空调负荷动态相关性的经济调峰方法。
图1为根据本发明一个实施例的计及光伏出力和空调负荷动态相关性的经济调峰方法流程图。
如图1所示,该计及光伏出力和空调负荷动态相关性的经济调峰方法包括以下步骤:
步骤S1,获取用户空调用电行为的控制温度和舒适温度区间、预设期间的室外环境温度和太阳光辐射度。
具体地,用户空调用电行为的控制温度和舒适温度区间数据通过调查海量用户并收集问卷采集;从气象数据采集系统(包括各个传感器,智能气象数据采集仪,GPRS DTU通讯模块)中获取整个夏季期间每日24点室外环境温度数据;从太阳辐射测量仪采集夏季期间每日24点的太阳辐射度;所述的气象数据采集系统、太阳辐射测量仪为一个系统,每小时采集的室外环境温度和太阳辐射量称为第i个变量,整个夏季期间共i×24×61有个变量。
可以理解的是,预设期间恶可以进行设置,在上述实施例中,预设期间为整个夏季期间每日24点。
步骤S2,根据预设期间的室外环境温度和太阳光辐射度计算光伏电源的电池温度。
具体地,通过采集的预设期间的室外环境温度和太阳光辐射度计算光伏电源的电池温度,通过如下公式计算:
式中,NOCT为工作单元温度,表示在入射辐射为0.8kW/m2、环境温度为20℃和无负载操作时的电池温度;GT是撞击PV阵列的太阳辐射[kW/m2];ηmp,STC是标准测试条件下的最大功率点效率[%];αp是功率的温度系数[%/℃],;TC,STC是在标准测试条件下的电池温度[25℃];τα为0.9;Ta是环境温度[℃]。
步骤S3,获取电力系统的光伏电源额定功率,根据所述太阳光辐射度、所述光伏电源的电池温度和所述电力系统的光伏电源额定功率求得光伏出力情况。
具体地,通过如下公式计算光伏出力情况:
YPV为获取接入配电网的光伏电源额定容量[kW];fPV是光伏电源系统的功率降额因数;GT,STC是在标准测试条件下的入射辐射[kW/m2];TC为光伏电源的电池温度[℃]。
可以理解的是,通过获取电力系统的光伏电源额定功率,并结合采集的太阳光辐射度和计算出的光伏电源的电池温度,来求得光伏出力情况。其中,不同的电力系统有不同的光伏电源额定功率。
步骤S4,根据用户空调用电行为的控制温度得到用户空调系统的开关状态,根据预设期间的室外环境温度和用户空调系统的开关状态计算空调负荷有功功率。
进一步地,在本发明的一个实施例中,根据用户空调用电行为的控制温度得到用户空调系统的开关状态,包括:
若用户室内温度大于用户空调用电行为的控制温度,则用户空调系统为开启状态;
若用户室内温度小于用户空调用电行为的控制温度,则用户空调系统的为关闭状态。
具体地,在夏季,空调的模式为制冷模式,以降低室内温度,用户空调用电行为的控制温度为用户空调显示的温度,若用户室内温度大于用户空调用电行为的控制温度,则用户空调系统为开启状态。若用户室内温度小于用户空调用电行为的控制温度,则用户空调系统的为关闭状态。
进一步地,根据预设期间的室外环境温度和用户空调系统的开关状态计算空调负荷有功功率,其中,空调负荷有功功率的计算公式来源为DesT软件,具体为:
其中,PAC为空调负荷有功功率,SAC为用户空调系统的开关状态,1表示空调开启,0表示空调关闭;COP为空调能耗比,QAC为等效热导率。
步骤S5,计算光伏出力情况和空调负荷有功功率的边缘累积概率分布函数,再根据边缘累积概率分布函数计算Copula函数求得联合概率分布函数,并转换为矩阵形式以得到初步多维变量矩阵。
进一步地,步骤S5进一步包括:
S51,利用矩估计方法分别对光伏出力情况和空调负荷功率计算边缘累积概率分布函数;
S52,对步骤S51中得到的边缘累积概率分布函数的分布类型和参数性质进行t检验,并判断假设检验是否通过,若是则进入步骤S53,若否则返回步骤S51;
S53,根据光伏出力情况和空调负荷功率的边缘累积概率分布函数,计算光伏出力情况和用户空调负荷有功功率的Copula函数表达式作为联合概率分布函数,并转换为矩阵形式以得到初步多维变量矩阵。
进一步地,计算光伏出力情况和用户空调负荷有功功率的Copula函数表达式类型包括:正态Copula函数、t-Copula函数、Gumbel-Copula、Clayton-Copula和Frank-Copula函数,并采用欧式平方距离最小化原则确定欧式平方距离中最小值所对应的Copula函数以作为联合概率分布函数。
具体地,利用矩估计方法分别计算光伏出力情况和空调负荷有功功率的边缘累积概率分布函数,得到光伏出力情况的边缘累积概率分布函数和空调负荷有功功率的边缘累积概率分布函数,根据各自的边缘累积概率分布函数计算Copula函数以得到联合概率分布函数。
可以理解的是,对各自的边缘累积概率分布函数计算Copula函数的类型包括:正态Copula函数、t-Copula函数、Gumbel-Copula、Clayton-Copula和Frank-Copula函数,在计算出Copula函数后通过欧式平方距离最小化原则确定出一个欧式平方距离最小的所对应的Copula函数作为联合概率分布函数,将联合概率分布函数转换为矩阵形式得到初步多维变量矩阵。
步骤S6,利用Markov转移概率矩阵对初步多维变量矩阵进行时序重构计算新的多维变量矩阵以得到新的联合概率分布函数。
进一步地,步骤S6进一步包括:
S61,在初步多维变量矩阵任意选取一组列向量作为新的多维变量矩阵的初始列,继而计算下一时刻多维变量的状态转移概率,依据最大状态转移概率原则选取下一时刻变量所属区域;
S62,对初步多维变量矩阵中从初始列开始按列依次判断所属区域,并寻找符合所属区域的多维列向量作为新的多维变量矩阵的下一列,得到新的多维变量矩阵,继而得到新的联合概率分布函数。
具体地,变量区域计算方法为:
Xt+1=Xmin+ε1,t(Xmax-Xmin)
其中,ε1,t为服从均匀分布的随机数,ε1,t∈[0,1];[Xmin,Xmax]为变量上下限。
Markov转移概率矩阵为:
Pij=P(Xn+1=si|Xn=sj),si,sj∈s
其中,Pij为状态转移概率;S为变量状态序列。
步骤S7,根据新的联合概率分布函数随机生成预设数量的光伏出力情况和空调负荷有功功率的边缘累积概率分布值,将光伏出力情况和空调负荷有功功率的边缘累积概率分布值代入对应的边缘累积概率分布函数得到光伏有功出力样本值和用户空调负荷功率样本值。
具体地,根据上述产生的新的联合概率分布函数随机生成预设数量的光伏出力情况的边缘累积概率分布值和空调负荷有功功率的边缘累积概率分布值,将光伏出力情况和空调负荷有功功率的边缘累积概率分布值分别代入各自对应的边缘累积概率分布函数得到光伏有功出力样本值和用户空调负荷功率样本值。
步骤S8,建立目标函数,由所述目标函数和约束条件组成规划问题,将所述光伏有功出力样本值作为电力系统中首节点的平衡功率输入,所述用户空调负荷功率样本值作为电力系统中PV节点的注入功率,根据所述舒适温度区间、所述光伏有功出力样本值和所述用户空调负荷功率样本值对所述规划问题进行迭代求解,获取所述目标函数的最小值为结果,根据所述结果得到用户空调用电行为模式调度方案。
可以理解的是,将所述光伏有功出力样本值作为电力系统中首节点的平衡功率输入,所述用户空调负荷功率样本值作为电力系统中PV节点的注入功率,来辅助求解目标函数和约束条件组成的规划问题。其中,目标函数为最大化电力系统中峰值负荷削减量,所述目标函数的控制变量为用户室内温度,状态变量为节点电压。
进一步地,目标函数为:
f=max|Pmax-Pmax’|
其中,Pmax,Pmax’分别为未计及和已计及光伏出力和空调负荷相关性的有功功率峰值负荷。
约束条件包括:
潮流方程等式约束:
P1=PPV’
Pi=PAC’,i∈{2,…,m}
其中,Ui、Pi和Qi分别为每个节点的电压、有功和无功功率;Gij和Bij为每个节点的电导和电纳;n表示电力系统的节点个数;m表示电力系统的PV节点个数;P1为首节点的有功功率,PPV’为光伏有功出力样本值,PAC’为空调负荷有功功率样本值,ΔPi为有功功率不平衡量,δij为节点i和节点j电压的相位差,ΔQi为无功功率不平衡量;
控制变量不等式约束:
为室内温度最小值,Tindoor为室内温度,为室内温度最大值;
电压上下限约束:
为节点i电压最小值,为节点i电压最大值,β为节点电压置信水平。
具体地,对由目标函数和约束条件组成规划问题进行迭代求解,以获取目标函数最小值为结果,获取用户空调用电24小时行为模式调度方案。
综上所述,通过进行具体的源荷一体化调度计算,来对本发明的方法进行详细说明,如图2所示。
(s1)用户空调用电行为的控制温度和舒适温度区间数据通过调查海量用户并收集问卷采集;从气象数据采集系统(包括各个传感器,智能气象数据采集仪,GPRS DTU通讯模块)中获取整个夏季期间每日24点室外环境温度数据;从太阳辐射测量仪采集夏季期间每日24点的太阳辐射度;所述的气象数据采集系统、太阳辐射测量仪为一个系统,每小时采集的室外环境温度和太阳辐射量称为第i个变量,整个夏季期间共i×24×61有个变量。
(s2)根据上述步骤(s1)采集的室外环境温度或太阳光辐射度计算第i个变量的光伏电池温度。
(s3)根据上述步骤(s1)采集的太阳辐射度和上述步骤(s2)计算的光伏电池温度并通过获取的接入配电网的光伏电源额定容量和功率降额因数,计算第i个变量的光伏出力情况。
(s4)根据上述步骤(s1)采集的用户空调用电行为的控制温度确定空调开关状态,若用户室内温度高于控制温度,则空调模式为开启状态,反之,则为关闭状态。
(s5)根据上述步骤(s1)采集的室外环境温度和步骤(s4)确定的空调开关状态,获取用户空调的等效热容、等效热阻、等效热耗率并确定仿真时间间隔,可以计算用户室内温度随时间变化轨迹。
(s6)根据上述步骤(s4)确定的空调开关状态,获取用户的空调能效比,计算用户空调负荷有功功率。
(s7)根据上述步骤(s3)计算的光伏有功出力和步骤(s6)的用户空调负荷有功功率,采用矩估计方法分别计算第i个变量的边缘累积概率分布函数。
(s8)对上述步骤(s7)计算的边缘累积概率分布的分布类型和参数性质进行t检验。通过假设检验的继续步骤(s8),未通过假设检验的返回步骤(s7)。
(s9)根据上述步骤(s1)采集的有功出力或热负荷功率和上述步骤(s7)计算的边缘累积概率分布函数并通过上述步骤(s8)检验的边缘概率分布公式,计算5个类型Copula函数,计算欧式平方距离。
(s10)取上述步骤(s9)计算的欧式平方距离中最小所对应的Copula函数作为联合概率分布函数,得到初步多维变量矩阵。
(s11)从上述步骤(s10)得到的初步多维变量矩阵任意选取一组列向量作为新的多维变量矩阵的初始列,继而计算下一时刻多维变量的状态转移概率,依据最大状态转移概率原则选取下一时刻变量所属的区域。
(s12)对初步多维变量矩阵中从初始列开始按列依次判断所属区域,并寻找符合所属区域的多维列向量作为新的多维变量矩阵的下一列,得到新的多维变量矩阵,继而得到新的联合概率分布函数。
(s13)根据上述步骤(s12)得到新的联合概率分布函数,随机生产各个变量的边缘累积概率分布值,代入至上述步骤(s7)的边缘累积概率分布函数中,得到光伏有功出力和用户空调负荷功率样本值;
(s14)将上述步骤(s13)计算的光伏有功出力样本值作为配电网(电力系统)首节点的平衡功率,用户空调负荷功率样本值作为配电网中PV节点的注入功率,结合上述步骤(s1)采集的用户空调用电行为的舒适温度区间数据,由目标函数和约束条件组成规划问题进行迭代求解,以获取目标函数最小值为结果,得到用户空调用电行为模式调度方案。目标函数是:
f=max|Pmax-Pmax’|
约束条件表达式如下:
P1=PPV’
Pi=PAC’,i∈{2,…,m}
根据本发明实施例提出的计及光伏出力和空调负荷动态相关性的经济调峰方法,通过将Copula函数多维马尔科夫链的转移概率矩阵相结合,获取了分布式光伏电源系统中的有功出力和用户空调负荷功率分布两部分变量的动态联合概率分布,继而得到两部分变量的边缘累积概率分布样本,作为电力系统节点的输入,将该方法与分布式光伏电源系统和空调冷负荷系统的传统分开计算方法进行比较,充分考虑了分布式光伏电源系统和空调冷负荷系统之间的动态相关性和时序性,因而可以准确描绘出两个系统通过耦合对源-荷的联合影响,大大提高电力系统随机生产仿真的精度。
该方法将分布式光伏电源系统和空调冷负荷系统看作一个整体进行分析,彰显不同类型能源供应商与用户协同合作的重要性,将需求侧与供给侧深度融合、统筹规划,促成能源运营管理更加高效。本方法中,计及了光伏发电出力和空调冷负荷动态相关性,显著提升空调负荷有功调节能力,缓解了夏季用电高峰期间的供需不平衡矛盾,降低最高峰值拉闸限电的可能性,调峰能力强,可靠性和经济性高。并且本方法可推广至更多分布式电源和不同类型负荷的配电系统,有效促进能源综合配置和协调规划,提高分布式电源的渗透率和发电效率。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (10)
1.一种计及光伏出力和空调负荷动态相关性的经济调峰方法,其特征在于,包括以下步骤:
S1,获取用户空调用电行为的控制温度和舒适温度区间、预设期间的室外环境温度和太阳光辐射度;
S2,根据所述预设期间的室外环境温度和太阳光辐射度计算光伏电源的电池温度;
S3,获取电力系统的光伏电源额定功率,根据所述太阳光辐射度、所述光伏电源的电池温度和所述电力系统的光伏电源额定功率求得光伏出力情况;
S4,根据所述用户空调用电行为的控制温度得到用户空调系统的开关状态,根据所述预设期间的室外环境温度和所述用户空调系统的开关状态计算空调负荷有功功率;
S5,计算所述光伏出力情况和所述空调负荷有功功率的边缘累积概率分布函数,再根据所述边缘累积概率分布函数计算Copula函数求得联合概率分布函数,并转换为矩阵形式以得到初步多维变量矩阵;
S6,利用Markov转移概率矩阵对所述初步多维变量矩阵进行时序重构计算新的多维变量矩阵以得到新的联合概率分布函数;
S7,根据所述新的联合概率分布函数随机生成预设数量的所述光伏出力情况和所述空调负荷有功功率的边缘累积概率分布值,将所述光伏出力情况和所述空调负荷有功功率的边缘累积概率分布值代入对应的边缘累积概率分布函数得到光伏有功出力样本值和用户空调负荷功率样本值;
S8,建立目标函数,由所述目标函数和约束条件组成规划问题,将所述光伏有功出力样本值作为电力系统中首节点的平衡功率输入,所述用户空调负荷功率样本值作为电力系统中PV节点的注入功率,根据所述舒适温度区间、所述光伏有功出力样本值和所述用户空调负荷功率样本值对所述规划问题进行迭代求解,获取所述目标函数的最小值为结果,根据所述结果得到用户空调用电行为模式调度方案。
2.根据权利要求1所述的方法,其特征在于,所述根据所述用户空调用电行为的控制温度得到用户空调系统的开关状态,包括:
若用户室内温度大于所述用户空调用电行为的控制温度,则所述用户空调系统为开启状态;
若用户室内温度小于所述用户空调用电行为的控制温度,则所述用户空调系统的为关闭状态。
3.根据权利要求1所述的方法,其特征在于,所述空调负荷有功功率的计算公式为:
其中,PAC为空调负荷有功功率,SAC为用户空调系统的开关状态,1表示空调开启,0表示空调关闭;COP为空调能耗比,QAC为等效热导率。
4.根据权利要求1所述的方法,其特征在于,所述S5包括:
S51,利用矩估计方法分别对所述光伏出力情况和所述空调负荷功率计算所述边缘累积概率分布函数;
S52,对步骤S51中得到的所述边缘累积概率分布函数的分布类型和参数性质进行t检验,并判断假设检验是否通过,若是则进入步骤S53,若否则返回步骤S51;
S53,根据所述光伏出力情况和所述空调负荷功率的边缘累积概率分布函数,计算所述光伏出力情况和所述用户空调负荷有功功率的Copula函数表达式作为所述联合概率分布函数,并转换为矩阵形式以得到所述初步多维变量矩阵。
5.根据权利要求4所述的方法,其特征在于,
计算所述光伏出力情况和所述用户空调负荷有功功率的Copula函数表达式类型包括:正态Copula函数、t-Copula函数、Gumbel-Copula、Clayton-Copula和Frank-Copula函数,并采用欧式平方距离最小化原则确定欧式平方距离中最小值所对应的Copula函数以作为所述联合概率分布函数。
6.根据权利要求1所述的方法,其特征在于,所述S6具体包括:
S61,在所述初步多维变量矩阵任意选取一组列向量作为所述新的多维变量矩阵的初始列,继而计算下一时刻多维变量的状态转移概率,依据最大状态转移概率原则选取下一时刻变量所属区域;
S62,对所述初步多维变量矩阵中从初始列开始按列依次判断所属区域,并寻找符合所属区域的多维列向量作为所述新的多维变量矩阵的下一列,得到所述新的多维变量矩阵,继而得到所述新的联合概率分布函数。
7.根据权利要求1所述的方法,其特征在于,所述Markov转移概率矩阵为:
Pij=P(Xn+1=si|Xn=sj),si,sj∈s
其中,Pij为状态转移概率;S为变量状态序列。
8.根据权利要求1所述的方法,其特征在于,所述目标函数为:
f=max|Pmax-Pmax’|
其中,Pmax,Pmax’分别为未计及和已计及光伏出力和空调负荷相关性的有功功率峰值负荷。
9.根据权利要求1所述的方法,其特征在于,所述约束条件包括:
潮流方程等式约束:
P1=PPV’
Pi=PAC’,i∈{2,…,m}
其中,Ui、Pi和Qi分别为每个节点的电压、有功和无功功率;Gij和Bij为每个节点的电导和电纳;n表示电力系统的节点个数;m表示电力系统的PV节点个数;P1为首节点的有功功率,PPV’为光伏有功出力样本值,PAC’为空调负荷有功功率样本值,ΔPi为有功功率不平衡量,δij为节点i和节点j电压的相位差,ΔQi为无功功率不平衡量;
控制变量不等式约束:
为室内温度最小值,Tindoor为室内温度,为室内温度最大值;
电压上下限约束:
为节点i电压最小值,为节点i电压最大值,β为节点电压置信水平。
10.根据权利要求1所述的方法,其特征在于,
所述目标函数为最大化电力系统中峰值负荷削减量,所述目标函数的控制变量为用户室内温度,状态变量为节点电压。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910551461.0A CN110323740B (zh) | 2019-06-24 | 2019-06-24 | 计及光伏出力和空调负荷动态相关性的经济调峰方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910551461.0A CN110323740B (zh) | 2019-06-24 | 2019-06-24 | 计及光伏出力和空调负荷动态相关性的经济调峰方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110323740A true CN110323740A (zh) | 2019-10-11 |
CN110323740B CN110323740B (zh) | 2020-08-21 |
Family
ID=68120153
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910551461.0A Active CN110323740B (zh) | 2019-06-24 | 2019-06-24 | 计及光伏出力和空调负荷动态相关性的经济调峰方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110323740B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111313403A (zh) * | 2020-02-26 | 2020-06-19 | 上海电力大学 | 一种基于马尔可夫随机场的低压配电系统网络拓扑识别方法 |
CN112290536A (zh) * | 2020-09-23 | 2021-01-29 | 电子科技大学 | 基于近端策略优化的电-热综合能源系统在线调度方法 |
CN113408924A (zh) * | 2021-06-30 | 2021-09-17 | 中国农业大学 | 基于统计机器学习的园区综合能源系统的规划方法 |
CN115657545A (zh) * | 2022-10-24 | 2023-01-31 | 海宁云多科技有限公司 | 一种基于空调与光伏联合运行的用电调控系统和方法 |
CN116454924A (zh) * | 2023-06-16 | 2023-07-18 | 华东交通大学 | 一种基于相关系数的柔性负荷功率调控方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102510079A (zh) * | 2011-10-23 | 2012-06-20 | 重庆市电力公司 | 利用太阳能发电的水源热泵冷电联合调度系统及调度方法 |
CN106709610A (zh) * | 2017-01-12 | 2017-05-24 | 浙江大学 | 一种微电网电储能、冰蓄冷联合优化调度方法 |
CN106992523A (zh) * | 2017-04-25 | 2017-07-28 | 国网上海市电力公司 | 一种用于含光伏和热负荷的电力系统的潮流计算方法 |
CN107294087A (zh) * | 2017-06-23 | 2017-10-24 | 清华大学 | 一种含气象能源的综合能源系统典型场景集生成方法 |
CN107294212A (zh) * | 2017-07-25 | 2017-10-24 | 山东大学 | 考虑不同空调负荷特性的微网双层优化调度方法及系统 |
CN108054750A (zh) * | 2017-12-08 | 2018-05-18 | 西南大学 | 一种计及夏季空调负荷差异特性的电网调峰控制系统与方法 |
US20190052120A1 (en) * | 2016-02-03 | 2019-02-14 | Gree Electric Appliances, Inc. Of Zhuhai | Photovoltaic Energy Storage Air Conditioner and Control Method Thereof |
CN109461091A (zh) * | 2018-05-25 | 2019-03-12 | 中国农业大学 | 考虑光伏和冷负荷相关性的用电负荷计算方法及信息系统 |
CN109659972A (zh) * | 2018-11-30 | 2019-04-19 | 国网浙江省电力有限公司经济技术研究院 | 基于Pair-copula理论的多点光伏出力概率预测方法和系统 |
CN109783841A (zh) * | 2018-11-28 | 2019-05-21 | 河海大学 | 一种基于多场景状态转移矩阵与条件概率抽样的光伏出力时间序列模拟方法 |
-
2019
- 2019-06-24 CN CN201910551461.0A patent/CN110323740B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102510079A (zh) * | 2011-10-23 | 2012-06-20 | 重庆市电力公司 | 利用太阳能发电的水源热泵冷电联合调度系统及调度方法 |
US20190052120A1 (en) * | 2016-02-03 | 2019-02-14 | Gree Electric Appliances, Inc. Of Zhuhai | Photovoltaic Energy Storage Air Conditioner and Control Method Thereof |
CN106709610A (zh) * | 2017-01-12 | 2017-05-24 | 浙江大学 | 一种微电网电储能、冰蓄冷联合优化调度方法 |
CN106992523A (zh) * | 2017-04-25 | 2017-07-28 | 国网上海市电力公司 | 一种用于含光伏和热负荷的电力系统的潮流计算方法 |
CN107294087A (zh) * | 2017-06-23 | 2017-10-24 | 清华大学 | 一种含气象能源的综合能源系统典型场景集生成方法 |
CN107294212A (zh) * | 2017-07-25 | 2017-10-24 | 山东大学 | 考虑不同空调负荷特性的微网双层优化调度方法及系统 |
CN108054750A (zh) * | 2017-12-08 | 2018-05-18 | 西南大学 | 一种计及夏季空调负荷差异特性的电网调峰控制系统与方法 |
CN109461091A (zh) * | 2018-05-25 | 2019-03-12 | 中国农业大学 | 考虑光伏和冷负荷相关性的用电负荷计算方法及信息系统 |
CN109783841A (zh) * | 2018-11-28 | 2019-05-21 | 河海大学 | 一种基于多场景状态转移矩阵与条件概率抽样的光伏出力时间序列模拟方法 |
CN109659972A (zh) * | 2018-11-30 | 2019-04-19 | 国网浙江省电力有限公司经济技术研究院 | 基于Pair-copula理论的多点光伏出力概率预测方法和系统 |
Non-Patent Citations (2)
Title |
---|
唐巍等: "考虑用户满意度的户用型微电网日前优化调度", 《高电压技术》 * |
祝明乐: "蓄冰空调与分布式光伏协调运行关键技术研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111313403A (zh) * | 2020-02-26 | 2020-06-19 | 上海电力大学 | 一种基于马尔可夫随机场的低压配电系统网络拓扑识别方法 |
CN111313403B (zh) * | 2020-02-26 | 2023-11-03 | 上海电力大学 | 一种基于马尔可夫随机场的低压配电系统网络拓扑识别方法 |
CN112290536A (zh) * | 2020-09-23 | 2021-01-29 | 电子科技大学 | 基于近端策略优化的电-热综合能源系统在线调度方法 |
CN112290536B (zh) * | 2020-09-23 | 2022-12-23 | 电子科技大学 | 基于近端策略优化的电-热综合能源系统在线调度方法 |
CN113408924A (zh) * | 2021-06-30 | 2021-09-17 | 中国农业大学 | 基于统计机器学习的园区综合能源系统的规划方法 |
CN113408924B (zh) * | 2021-06-30 | 2024-02-02 | 中国农业大学 | 基于统计机器学习的园区综合能源系统的规划方法 |
CN115657545A (zh) * | 2022-10-24 | 2023-01-31 | 海宁云多科技有限公司 | 一种基于空调与光伏联合运行的用电调控系统和方法 |
CN115657545B (zh) * | 2022-10-24 | 2023-05-09 | 海宁云多科技有限公司 | 一种基于空调与光伏联合运行的用电调控系统和方法 |
CN116454924A (zh) * | 2023-06-16 | 2023-07-18 | 华东交通大学 | 一种基于相关系数的柔性负荷功率调控方法 |
CN116454924B (zh) * | 2023-06-16 | 2023-09-29 | 华东交通大学 | 一种基于相关系数的柔性负荷功率调控方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110323740B (zh) | 2020-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110323740A (zh) | 计及光伏出力和空调负荷动态相关性的经济调峰方法 | |
Hakimi et al. | Optimal planning of a smart microgrid including demand response and intermittent renewable energy resources | |
CN103353979B (zh) | 一种分布式电源的优化选址与定容方法 | |
CN108596525B (zh) | 冷热电气多能互补的微能源网鲁棒优化调度方法 | |
CN113964885B (zh) | 一种基于态势感知的电网无功主动预测及控制方法 | |
CN109546659A (zh) | 基于随机矩阵和智能场景匹配的配电网无功优化方法 | |
CN110991773A (zh) | 一种面向风电消纳的两阶段源荷储优化调度方法 | |
CN112419087A (zh) | 一种聚合综合能源楼宇的虚拟电厂日前优化调度方法 | |
CN110380444A (zh) | 一种基于变结构Copula的多场景下分散式风电有序接入电网的容量规划方法 | |
CN110474370A (zh) | 一种空调可控负荷、光伏储能系统的协同控制系统及方法 | |
CN106570729A (zh) | 基于空调可靠性影响因素的区域聚类方法 | |
CN110417045A (zh) | 一种针对交直流混合微电网容量配置的优化方法 | |
Munankarmi et al. | Quantification of load flexibility in residential buildings using home energy management systems | |
Han et al. | Physical-data fusion modeling method for energy consumption analysis of smart building | |
Alden et al. | Artificial intelligence method for the forecast and separation of total and hvac loads with application to energy management of smart and nze homes | |
Cui et al. | Improved normal-boundary intersection algorithm: A method for energy optimization strategy in smart buildings | |
CN114462696A (zh) | 基于td3的综合能源系统源荷协同运行优化方法 | |
Luo et al. | Optimal scheduling of micro-energy grid based on pareto frontier under uncertainty and pollutant emissions | |
CN117096856B (zh) | 考虑配网三相电价和相间电压不平衡的虚拟电厂调度方法 | |
Changliang et al. | Energy management strategy research for residential microgrid considering virtual energy storage system at demand side | |
CN109586309B (zh) | 基于大数据自由熵理论和场景匹配的配电网无功优化方法 | |
Zhang et al. | Multi–objective cluster partition method for distribution network considering uncertainties of distributed generations and loads | |
CN114819662A (zh) | 一种考虑用户热需求的蓄热电锅炉群电热调度方法及系统 | |
Maochun et al. | Typical scenario construction for fast assessment of renewable energy absorption capacity | |
Abou Houran et al. | Energy management of microgrid in smart building considering air temperature impact |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |