CN111126772A - 一种模拟电力现货市场平衡的方法及*** - Google Patents
一种模拟电力现货市场平衡的方法及*** Download PDFInfo
- Publication number
- CN111126772A CN111126772A CN201911174612.1A CN201911174612A CN111126772A CN 111126772 A CN111126772 A CN 111126772A CN 201911174612 A CN201911174612 A CN 201911174612A CN 111126772 A CN111126772 A CN 111126772A
- Authority
- CN
- China
- Prior art keywords
- power
- output
- quotation
- probability density
- renewable energy
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000010248 power generation Methods 0.000 claims abstract description 46
- 230000011218 segmentation Effects 0.000 claims abstract description 24
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 15
- 238000000605 extraction Methods 0.000 claims description 19
- 238000004146 energy storage Methods 0.000 claims description 18
- 238000012360 testing method Methods 0.000 claims description 17
- 230000005540 biological transmission Effects 0.000 claims description 12
- 230000005611 electricity Effects 0.000 claims description 10
- 238000007599 discharging Methods 0.000 claims description 6
- 238000012886 linear function Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000002093 peripheral effect Effects 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 5
- 238000000638 solvent extraction Methods 0.000 claims 1
- 238000004088 simulation Methods 0.000 abstract description 5
- 238000005286 illumination Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 5
- 238000013507 mapping Methods 0.000 description 5
- 238000005315 distribution function Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/067—Enterprise or organisation modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q30/00—Commerce
- G06Q30/02—Marketing; Price estimation or determination; Fundraising
- G06Q30/0201—Market modelling; Market analysis; Collecting market data
- G06Q30/0206—Price or cost determination based on market factors
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S50/00—Market activities related to the operation of systems integrating technologies related to power network operation or related to communication or information technologies
- Y04S50/14—Marketing, i.e. market research and analysis, surveying, promotions, advertising, buyer profiling, customer management or rewards
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Physics & Mathematics (AREA)
- Finance (AREA)
- Marketing (AREA)
- Accounting & Taxation (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Educational Administration (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本申请公开了一种模拟电力现货市场平衡的方法及***,该方法包括:首先获取可再生能源的出力模型,然后根据出力模型确定可再生能源的概率密度函数,利用多元线性蒙特卡洛模拟方法对该概率密度函数进行线性分割,使可再生能源出力与概率密度函数相匹配,再根据线性分割结果和节点报价平衡条件计算得出不同节点的发电出力和报价,最后,根据不同节点的发电出力和报价,利用公式计算得出电力***报价匹配平衡结果。该***包括:可再生能源出力模型获取模块、概率密度函数确定模块、线性分割模块、节点发电出力和报价模块以及电力***报价模块。通过本申请,能够大大提高模拟结果的准确性。
Description
技术领域
本申请涉及电力市场运行管理技术领域,特别是涉及一种模拟电力现货市场平衡的方法及***。
背景技术
随着电力市场的发展,电力现货市场逐渐成为我国电力市场的关键组成部分。如何模拟电力现货市场的平衡,从而对市场运行的稳定性提供有效指导,是个重要的技术问题。
目前模拟电力现货市场平衡的方法,主要采用蒙特卡洛模拟仿真方法,该方法是一种数据随机抽取的方法。具体地,针对10KW、20KW、30KW等不同发电***中,光伏、风机出力以及用户用电需求有各自的取值范围,每一个取值有与其相对应的取值概率。在供应侧发电出力确定的情况下,对用户侧用电需求进行随机抽取;在用户侧用户需求确定的情况下,对供应侧的风速、光照强度等进行抽取。在供应侧发电出力和用户侧用户需求都不确定的情况下,对供应侧和用户侧随机抽取。
然而,目前模拟电力现货市场平衡的方法中,当供应侧或用户侧单侧需求确定时,抽取结果相对较准确。当供应侧和用户侧两侧都进行抽取时,由于供应侧和用户侧随机抽取,极有可能出现采用风速的概率密度抽取光照强度,或者采用光照强度的概率密度抽取风速,这种概率不匹配的情况,会导致模拟结果不够准确。
发明内容
本申请提供了一种模拟电力现货市场平衡的方法及***,以解决现有技术中的模拟电力现货市场平衡的方法导致模拟结果准确性较低的问题。
为了解决上述技术问题,本申请实施例公开了如下技术方案:
一种模拟电力现货市场平衡的方法,所述方法包括:
获取可再生能源的出力模型,所述可再生能源包括风能和光伏,所述出力模型包括:风机出力模型和光伏发电出力模型;
根据所述出力模型,确定所述可再生能源的概率密度函数;
利用多元线性蒙特卡洛模拟方法,对所述可再生能源的概率密度函数进行线性分割,使所述可再生能源的出力与所述可再生能源的概率密度函数相匹配;
根据线性分割结果和节点报价平衡条件计算得出不同节点的发电出力和报价,所述节点为电力***的节点,输电容量约束为Vxy(ρx,t-ρy,t)≤TRxy,qxh,t为火电机组在t时间段和在节点x的成交出力之和,qxj,bt、qxu,bt为风机j、光伏机组u在x节点、t时间段的成交出力之和,nz为x周边节点数量,ρx,t和ρy,t为节点x和y的电压相位角,Dx,t为节点x处的负荷,该等式在任意的t时间段内成立;
根据不同节点的发电出力和报价,利用公式计算得出电力***报价匹配平衡结果,其中,qxh,bt=q1h,bt+q2h,bt+…+qxh,bt,qxj,bt=q1j,bt+q2j,bt+…+qxj,bt,qxu,bt=q1u,bt+p2u,bt+…+qxu,bt,报价约束为:qxy,bt≤qx max,0≤qjy,bt≤qj max,0≤quy,bt≤qu max,qi,t为火电机组出力,qj,t、qu,t分别为风电、光伏在t时间段内的发电出力,qLOSS为输电功率损耗,Dtotal为***总负荷,DCUT为***切负荷,qo,t为储能设备充放电功率,当储能设备放电时,取值为正,储能设备充电时取值为负。
可选地,所述利用多元线性蒙特卡洛模拟方法,对所述可再生能源的概率密度函数进行线性分割,使所述可再生能源的出力与所述可再生能源的概率密度函数相匹配的方法,包括:
确定输入变量和输出变量的概率密度函数,所述输入变量包括:用户用电需求报价、火电机组报价、风电机组报价以及光伏发电机组报价,所述输出变量包括:用户用电负荷;
利用公式Yi=Y0-i+Ji -1[Xi-μ(Xi)],对负荷以及所述可再生能源的概率密度函数的取值概率区间进行区域划分,其中,i=1,2,3……m,用户用电负荷的概率密度函数被分割为m部分,Yi是第i部分的***输出指标,Y0-i是电力电量平衡公式中用户用电负荷的抽取概率空间,且gs(Y0-i)=μ(Xi),Ji是Y0-i的雅克比矩阵,Xi为第i部分的随机输入变量,μ(Xi)是第i部分的输入变量的平均值;
利用公式gs(Y0-i)=μ(Xi)计算得出用户用电负荷的抽取概率空间进行线性化,其中,gs为线性函数。
可选地,当i取值为1时,利用公式Y=Y0+J0 -1[X-μ(X)],对负荷以及所述可再生能源的概率密度函数的取值概率区间进行区域划分。
一种模拟电力现货市场平衡的***,所述***包括:
可再生能源出力模型获取模块,用于获取可再生能源的出力模型,所述可再生能源包括风能和光伏,所述出力模型包括:风机出力模型和光伏发电出力模型;
概率密度函数确定模块,用于根据所述出力模型,确定所述可再生能源的概率密度函数;
线性分割模块,用于利用多元线性蒙特卡洛模拟方法,对所述可再生能源的概率密度函数进行线性分割,使所述可再生能源的出力与所述可再生能源的概率密度函数相匹配;
节点发电出力和报价模块,用于根据线性分割结果和节点报价平衡条件计算得出不同节点的发电出力和报价,所述节点为电力***的节点,输电容量约束为Vxy(ρx,t-ρy,t)≤TRxy,qxh,t为火电机组在t时间段和在节点x的成交出力之和,qxj,bt、qxu,bt分别为风机j、光伏机组u在x节点、t时间段的成交出力之和,nz为x周边节点数量,ρx,t和ρy,t为节点x和y的电压相位角,Dx,t为节点x处的负荷,该等式在任意的t时间段内成立;
电力***报价模块,用于根据不同节点的发电出力和报价,利用公式计算得出电力***报价匹配平衡结果,其中,qxh,bt=q1h,bt+q2h,bt+…+qxh,bt,qxj,bt=q1j,bt+q2j,bt+…+qxj,bt,qxu,bt=q1u,bt+p2u,bt+…+qxu,bt,报价约束为:qxy,bt≤qx max,0≤qjy,bt≤qj max,0≤quy,bt≤qu max,qi,t为火电机组出力,qj,t、qu,t为风电、光伏在t时间段内的发电出力,qLOSS为输电功率损耗,Dtotal为***总负荷,DCUT为***切负荷,qo,t为储能设备充放电功率,当储能设备放电时,取值为正,充电时取值为负。
可选地,所述线性分割模块包括:
取值概率区间获取单元,用于获取输入输出变量概率密度所确定的取值概率区间,所述输入变量包括:用户用电需求报价、火电机组报价、风电机组报价以及光伏发电机组报价,所述输出变量包括:用户用电负荷;
区域划分单元,用于利用公式Yi=Y0-i+Ji -1[Xi-μ(Xi)],对负荷以及所述可再生能源的概率密度函数的取值概率区间进行区域划分,其中,i=1,2,3……m,用户用电负荷的概率密度函数被分割为m部分,Yi是第i部分的***输出指标,Y0-i是电力电量平衡公式中用户用电负荷的抽取概率空间,且gs(Y0-i)=μ(Xi),Ji是Y0-i的雅克比矩阵,Xi为第i部分的随机输入变量,μ(Xi)是第i部分的输入变量的平均值;
线性化单元,用于利用公式gs(Y0-i)=μ(Xi)计算得出用户用电负荷的抽取概率空间进行线性化,其中,gs为线性函数。
本申请的实施例提供的技术方案可以包括以下有益效果:
本申请提供一种模拟电力现货市场平衡的方法,该方法首先获取可再生能源的出力模型,然后根据出力模型确定可再生能源的概率密度函数,利用多元线性蒙特卡洛模拟方法对该概率密度函数进行线性分割,使可再生能源出力与概率密度函数相匹配,再根据线性分割结果和节点报价平衡条件计算得出不同节点的发电出力和报价,最后,根据不同节点的发电出力和报价,利用公式计算得出电力***报价匹配平衡结果。本实施例通过对现有的蒙特卡洛模型进行改进,建立可再生能源出力的概率区间与用户用电负荷概率区间之间的映射关系,且利用该映射关系指导用户进行报价,有利于提高蒙特卡洛模型的计算精度,通过线性切割的方法将用户用电需求报价、风电、光伏发电出力报价的概率区间切割成相互对应的若干区间,在对应的区间内按照正态分布概率及一定的规律进行随机抽取,抽取结果更加科学合理,更加与电力市场运行的现实情况更加相符,因此,这种方法有利于提高抽取结果和报价结果的准确性,从而提高对电力现货市场平衡的模拟结果的准确性。
本申请还提供一种模拟电力现货市场平衡的***,该***主要包括:可再生能源出力模型获取模块、概率密度函数确定模块、线性分割模块、节点发电出力和报价模块以及电力***报价模块。通过线性分割模块,对可再生能源的概率密度函数进行线性分割,从而使可再生能源的出力的概率区间与用户用电负荷概率区间之间建立映射关系,提高两种概率区间的匹配度,有利于提高模拟结果的准确性。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本申请。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本申请的实施例,并与说明书一起用于解释本申请的原理。
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,对于本领域普通技术人员而言,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例所提供的一种模拟电力现货市场平衡的方法的流程示意图;
图2为本申请实施例中概率密度区间线性分割示意图;
图3为本申请实施例中多元线性MC方法计算流程图;
图4为本申请实施例所提供的一种模拟电力现货市场平衡的***的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
为了更好地理解本申请,下面结合附图来详细解释本申请的实施方式。
实施例一
参见图1,图1为本申请实施例所提供的一种模拟电力现货市场平衡的方法的流程示意图。由图1可知,本实施例中模拟电力现货市场平衡的方法主要包括如下步骤:
S1:获取可再生能源的出力模型。
本实施例中的可再生能源包括风能和光伏,可再生能源的出力模型包括:风机出力模型和光伏发电出力模型。
本实施例中风机出力模型采用:其中,QW(v)为风机出力功率,Qr为风机的额定功率,v为风速,vc为切入风速,vr为额定风速,vf为切出风速,风速v是随机性较强的因素,本实施例中采用韦伯分布来确定v的概率分布:
本实施例中光伏发电出力模型采用:QPV=LMθ,其中,QPV为光伏发电出力功率,L为光照强度,M为受光面积,θ为发电效率。其中θ主要受到环境温度影响,其计算公式为:θ=θTEST[1-ψ(T-TTEST)],该公式中,QTEST为标准测试条件下的转化效率,ψ为功率温度系数,T为实际温度,TTEST为测试标准情况下的温度。光照强度L是随机性较强的因素,本实施例中采用贝塔分布函数确定L的概率分布。贝塔分布函数如下: 其中,fL(L)为光照强度的概率密度函数,Lmax为太阳光照强度的最大偏差值,ξ和η分别为平均偏差值和标准偏差值。
继续参见图1可知,获取到可再生能源的出力模型之后,执行步骤S2:根据出力模型,确定可再生能源的概率密度函数。
S3:利用多元线性蒙特卡洛模拟方法,对可再生能源的概率密度函数进行线性分割,使可再生能源的出力与可再生能源的概率密度函数相匹配。
本实施例中的蒙特卡洛模拟方法即MC(Monte Carlo,蒙特卡洛)方法,具体地,步骤S3又包括如下过程:
S31:确定输入变量和输出变量的概率密度函数。
其中,输入变量包括:用户用电需求报价、火电机组报价、风电机组报价以及光伏发电机组报价。输出变量包括:用户用电负荷。
本实施例中,用户用电负荷即用户用电需求,用户用电需求报价是指用户根据自身的用电负荷而上报的用电电价。
S32:利用公式Yi=Y0-i+Ji -1[Xi-μ(Xi)],对负荷以及可再生能源的概率密度函数的取值概率区间进行区域划分。
其中,i=1,2,3……m,用户用电负荷的概率密度函数被分割为m部分,Yi是第i部分的***输出指标,Y0-i是电力电量平衡公式中用户用电负荷的抽取概率空间,且gs(Y0-i)=μ(Xi),Ji是Y0-i的雅克比矩阵,Xi为第i部分的随机输入变量,μ(Xi)是第i部分的输入变量的平均值。概率密度函数的取值概率区间即pdf参数。
进一步地,本实施例中当i取值为1时,利用公式Y=Y0+J0 -1[X-μ(X)],对负荷以及所述可再生能源的概率密度函数的取值概率区间进行区域划分。电力电量平衡公式是以Y0为中心进行线性化,对于单一不确定变量存在的情况,可以采用公式Y=Y0+J0 -1[X-μ(X)]进行线性化。当i≥2时,具备多元线性化特征,利用公式Yi=Y0-i+Ji -1[Xi-μ(Xi)]进行线性化。
本实施例中概率密度区间线性分割示意图如图2所示。由图2可知,在多元线性在多元线性MC方法中,用户用电需求Dx,t的概率分布可以被分割为m个部分,每个部分以均值为中心对称分布,每个区域为R0-i,与之相对应的是Y0-i(i=1,2,3,…,m),当i=1时,将按照公式Y=Y0+J0 -1[X-μ(X)]的线性MC方法进行求解。当i≥2时,具备多元线性特征。当Dx,t的概率分布函数确定后,风电、光伏两类报价不确定机组的出力概率分布函数也将划分为m个部分,并且与Dx,t的概率分布一一对应。其中m个部分中,每个部分采用Y=Y0+J0 -1[X-μ(X)]的方法进行线性化:
M1:Y1=Y0-1+J1 -1[X1-μ(X1)],
M2:Y2=Y0-2+J2 -1[X2-μ(X2)],
……,
Mm:Ym=Y0-m+Jm -1[Xm-μ(Xm)。
S33:利用公式gs(Y0-i)=μ(Xi)计算得出用户用电负荷的抽取概率空间进行线性化。其中,gs为线性函数。
通过以上步骤S31-S33可知,在电力现货市场中,市场交易频率增大,不确定变量增多,本实施例采用多元线性蒙特卡洛模拟方法,对可再生能源的概率密度函数进行线性分割,有利于提高抽取匹配度。本实施例利用多元线性化蒙特卡洛模拟方法,能够建立可再生能源发电出力的概率区间与用户用电负荷概率区间之间的映射关系,利用该映射关系指导用户进行报价。由于对可再生能源的概率密度函数进行线性分割,不同的可再生能源发电出力概率区间匹配其相应的用户用电概率密度区间,抽查时避免随机性,有利于大大提高计算结果的准确性。
本实施例中节点为电力***的节点,每个用于接入的地方可以是一个节点,每个发电机组接入的地方也可以是一个节点,每个变电站也可以是一个节点。
输电容量约束为Vxy(ρx,t-ρy,t)≤TRxy,qxh,t为火电机组在t时间段和在节点x的成交出力之和,qxj,bt、qxu,bt为风机j、光伏机组u在x节点、t时间段的成交出力之和,nz为x周边节点数量,ρx,t和ρy,t为节点x和y的电压相位角,Dx,t为节点x处的负荷,该等式在任意的t时间段内成立。
其中,qxh,bt=q1h,bt+q2h,bt+…+qxh,bt,qxj,bt=q1j,bt+q2j,bt+…+qxj,bt,qxu,bt=q1u,bt+p2u,bt+…+qxu,bt,报价约束为:qxy,bt≤qx max,0≤qjy,bt≤qj max,0≤quy,bt≤qu max,qi,t为火电机组出力,qj,t、qu,t分别为风电、光伏在t时间段内的发电出力,t以每小时为一个时间段,qLOSS为输电功率损耗,Dtotal为***总负荷,DCUT为***切负荷,qo,t为储能设备充放电功率,当储能设备放电时,取值为正,储能设备充电时取值为负。
本实施例中采用多元线性蒙特卡洛方法模拟电力现货市场平衡时,多元线性MC方法的计算流程可以参见图3,图3中K为试验次数,Kmax为设定的试验次数最大值。
实施例二
在图1-图3所示实施例的基础之上参见图4,图4为本申请实施例所提供的一种模拟电力现货市场平衡的***的结构示意图。由图4可知,本实施例中模拟电力现货市场平衡的***,主要包括:可再生能源出力模型获取模块、概率密度函数确定模块、线性分割模块、节点发电出力和报价模块以及电力***报价模块。
其中,可再生能源出力模型获取模块,用于获取可再生能源的出力模型,可再生能源包括风能和光伏,出力模型包括:风机出力模型和光伏发电出力模型。概率密度函数确定模块用于根据出力模型,确定可再生能源的概率密度函数。线性分割模块,用于利用多元线性蒙特卡洛模拟方法,对可再生能源的概率密度函数进行线性分割,使可再生能源的出力与可再生能源的概率密度函数相匹配。节点发电出力和报价模块,用于根据线性分割结果和节点报价平衡条件计算得出不同节点的发电出力和报价,节点为电力***的节点,输电容量约束为Vxy(ρx,t-ρy,t)≤TRxy,qxh,t为火电机组在t时间段和在节点x的成交出力之和,qxj,bt、qxu,bt分别为风机j、光伏机组u在x节点、t时间段的成交出力之和,nz为x周边节点数量,ρx,t和ρy,t为节点x和y的电压相位角,Dx,t为节点x处的负荷,该等式在任意的t时间段内成立。电力***报价模块,用于根据不同节点的发电出力和报价,利用公式计算得出电力***报价匹配平衡结果,其中,qxh,bt=q1h,bt+q2h,bt+…+qxh,bt,qxj,bt=q1j,bt+q2j,bt+…+qxj,bt,qxu,bt=q1u,bt+p2u,bt+…+qxu,bt,报价约束为:qxy,bt≤qx max,0≤qjy,bt≤qj max,0≤quy,bt≤qu max,qi,t为火电机组出力,qj,t、qu,t为风电、光伏在t时间段内的发电出力,qLOSS为输电功率损耗,Dtotal为***总负荷,DCUT为***切负荷,qo,t为储能设备充放电功率,当储能设备放电时,取值为正,储能设备充电时取值为负。
进一步地,线性分割模块又包括:取值概率区间获取单元、区域划分单元和线性化单元。其中,取值概率区间获取单元,用于获取输入输出变量概率密度所确定的取值概率区间,输入变量包括:用户用电需求报价、火电机组报价、风电机组报价以及光伏发电机组报价,输出变量包括:用户用电负荷。区域划分单元,用于利用公式Yi=Y0-i+Ji -1[Xi-μ(Xi)],对负荷以及可再生能源的概率密度函数的取值概率区间进行区域划分,其中,i=1,2,3……m,用户用电负荷的概率密度函数被分割为m部分,Yi是第i部分的***输出指标,Y0-i是电力电量平衡公式中用户用电负荷的抽取概率空间,且gs(Y0-i)=μ(Xi),Ji是Y0-i的雅克比矩阵,Xi为第i部分的随机输入变量,μ(Xi)是第i部分的输入变量的平均值。线性化单元,用于利用公式gs(Y0-i)=μ(Xi)计算得出用户用电负荷的抽取概率空间进行线性化,其中,gs为线性函数。
该实施例中模拟电力现货市场平衡的***的工作原理和工作方法,在图1-图3所示的实施例中已经详细阐述,两个实施例之间可以互相参照,在此不再赘述。
以上所述仅是本申请的具体实施方式,使本领域技术人员能够理解或实现本申请。对这些实施例的多种修改对本领域的技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下,在其它实施例中实现。因此,本申请将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (7)
1.一种模拟电力现货市场平衡的方法,其特征在于,所述方法包括:
获取可再生能源的出力模型,所述可再生能源包括风能和光伏,所述出力模型包括:风机出力模型和光伏发电出力模型;
根据所述出力模型,确定所述可再生能源的概率密度函数;
利用多元线性蒙特卡洛模拟方法,对所述可再生能源的概率密度函数进行线性分割,使所述可再生能源的出力与所述可再生能源的概率密度函数相匹配;
根据线性分割结果和节点报价平衡条件计算得出不同节点的发电出力和报价,所述节点为电力***的节点,输电容量约束为Vxy(ρx,t-ρy,t)≤TRxy,qxh,t为火电机组在t时间段和在节点x的成交出力之和,qxj,bt、qxu,bt为风机j、光伏机组u在x节点、t时间段的成交出力之和,nz为x周边节点数量,ρx,t和ρy,t为节点x和y的电压相位角,Dx,t为节点x处的负荷,该等式在任意的t时间段内成立;
根据不同节点的发电出力和报价,利用公式计算得出电力***报价匹配平衡结果,其中,qxh,bt=q1h,bt+q2h,bt+…+qxh,bt,qxj,bt=q1j,bt+q2j,bt+…+qxj,bt,qxu,bt=q1u,bt+p2u,bt+…+qxu,bt,报价约束为:qxy,bt≤qx max,0≤qjy,bt≤qj max,0≤quy,bt≤qu max,qi,t为火电机组出力,qj,t、qu,t分别为风电、光伏在t时间段内的发电出力,qLOSS为输电功率损耗,Dtotal为***总负荷,DCUT为***切负荷,qo,t为储能设备充放电功率,当储能设备放电时,取值为正,储能设备充电时取值为负。
3.根据权利要求1所述的一种模拟电力现货市场平衡的方法,其特征在于,所述利用多元线性蒙特卡洛模拟方法,对所述可再生能源的概率密度函数进行线性分割,使所述可再生能源的出力与所述可再生能源的概率密度函数相匹配的方法,包括:
确定输入变量和输出变量的概率密度函数,所述输入变量包括:用户用电需求报价、火电机组报价、风电机组报价以及光伏发电机组报价,所述输出变量包括:用户用电负荷;
利用公式Yi=Y0-i+Ji -1[Xi-μ(Xi)],对负荷以及所述可再生能源的概率密度函数的取值概率区间进行区域划分,其中,i=1,2,3……m,用户用电负荷的概率密度函数被分割为m部分,Yi是第i部分的***输出指标,Y0-i是电力电量平衡公式中用户用电负荷的抽取概率空间,且gs(Y0-i)=μ(Xi),Ji是Y0-i的雅克比矩阵,Xi为第i部分的随机输入变量,μ(Xi)是第i部分的输入变量的平均值;
利用公式gs(Y0-i)=μ(Xi)计算得出用户用电负荷的抽取概率空间进行线性化,其中,gs为线性函数。
4.根据权利要求3所述的一种模拟电力现货市场平衡的方法,其特征在于,当i取值为1时,利用公式Y=Y0+J0 -1[X-μ(X)],对负荷以及所述可再生能源的概率密度函数的取值概率区间进行区域划分。
5.一种模拟电力现货市场平衡的***,其特征在于,所述***包括:
可再生能源出力模型获取模块,用于获取可再生能源的出力模型,所述可再生能源包括风能和光伏,所述出力模型包括:风机出力模型和光伏发电出力模型;
概率密度函数确定模块,用于根据所述出力模型,确定所述可再生能源的概率密度函数;
线性分割模块,用于利用多元线性蒙特卡洛模拟方法,对所述可再生能源的概率密度函数进行线性分割,使所述可再生能源的出力与所述可再生能源的概率密度函数相匹配;
节点发电出力和报价模块,用于根据线性分割结果和节点报价平衡条件计算得出不同节点的发电出力和报价,所述节点为电力***的节点,输电容量约束为Vxy(ρx,t-ρy,t)≤TRxy,qxh,t为火电机组在t时间段和在节点x的成交出力之和,qxj,bt、qxu,bt分别为风机j、光伏机组u在x节点、t时间段的成交出力之和,nz为x周边节点数量,ρx,t和ρy,t为节点x和y的电压相位角,Dx,t为节点x处的负荷,该等式在任意的t时间段内成立;
电力***报价模块,用于根据不同节点的发电出力和报价,利用公式计算得出电力***报价匹配平衡结果,其中,qxh,bt=q1h,bt+q2h,bt+…+qxh,bt,qxj,bt=q1j,bt+q2j,bt+…+qxj,bt,qxu,bt=q1u,bt+p2u,bt+…+qxu,bt,报价约束为:qxy,bt≤qx max,0≤qjy,bt≤qj max,0≤quy,bt≤qu max,qi,t为火电机组出力,qj,t、qu,t为风电、光伏在t时间段内的发电出力,qLOSS为输电功率损耗,Dtotal为***总负荷,DCUT为***切负荷,qo,t为储能设备充放电功率,当储能设备放电时,取值为正,储能设备充电时取值为负。
6.根据权利要求5所述的一种模拟电力现货市场平衡的***,其特征在于,所述线性分割模块包括:
取值概率区间获取单元,用于获取输入输出变量概率密度所确定的取值概率区间,所述输入变量包括:用户用电需求报价、火电机组报价、风电机组报价以及光伏发电机组报价,所述输出变量包括:用户用电负荷
区域划分单元,用于利用公式Yi=Y0-i+Ji -1[Xi-μ(Xi)],对负荷以及所述可再生能源的概率密度函数的取值概率区间进行区域划分,其中,i=1,2,3……m,用户用电负荷的概率密度函数被分割为m部分,Yi是第i部分的***输出指标,Y0-i是电力电量平衡公式中用户用电负荷的抽取概率空间,且gs(Y0-i)=μ(Xi),Ji是Y0-i的雅克比矩阵,Xi为第i部分的随机输入变量,μ(Xi)是第i部分的输入变量的平均值;
线性化单元,用于利用公式gs(Y0-i)=μ(Xi)计算得出用户用电负荷的抽取概率空间进行线性化,其中,gs为线性函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911174612.1A CN111126772A (zh) | 2019-11-26 | 2019-11-26 | 一种模拟电力现货市场平衡的方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911174612.1A CN111126772A (zh) | 2019-11-26 | 2019-11-26 | 一种模拟电力现货市场平衡的方法及*** |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111126772A true CN111126772A (zh) | 2020-05-08 |
Family
ID=70496770
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911174612.1A Pending CN111126772A (zh) | 2019-11-26 | 2019-11-26 | 一种模拟电力现货市场平衡的方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111126772A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111784203A (zh) * | 2020-07-28 | 2020-10-16 | 南方电网能源发展研究院有限责任公司 | 适用于发电机组参与的电力现货市场风险模拟分析方法 |
-
2019
- 2019-11-26 CN CN201911174612.1A patent/CN111126772A/zh active Pending
Non-Patent Citations (1)
Title |
---|
曾鸣等: "基于蒙特卡罗方法的电力***电源侧协调规划模拟仿真研究", 《华北电力大学学报(自然科学版)》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111784203A (zh) * | 2020-07-28 | 2020-10-16 | 南方电网能源发展研究院有限责任公司 | 适用于发电机组参与的电力现货市场风险模拟分析方法 |
CN111784203B (zh) * | 2020-07-28 | 2021-03-16 | 南方电网能源发展研究院有限责任公司 | 适用于发电机组参与的电力现货市场风险模拟分析方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Optimal location planning of renewable distributed generation units in distribution networks: an analytical approach | |
CN104156892B (zh) | 一种有源配电网电压跌落仿真与评估方法 | |
Mosadeghy et al. | A time-dependent approach to evaluate capacity value of wind and solar PV generation | |
CN105337315B (zh) | 一种风光蓄互补独立微电网高维多目标优化配置方法 | |
CN105656031A (zh) | 基于高斯混合分布特征的含风电电力***安全风险评估方法 | |
CN108364117B (zh) | 一种考虑光伏电站元件可靠性的电网风险评估方法 | |
CN110112790B (zh) | 有源配电网运行效率的评价方法、装置、设备及存储介质 | |
CN104778352A (zh) | 基于stft算子的七参数光伏电池输出特性建模方法 | |
CN103530473A (zh) | 一种含大规模光伏电站的电力***随机生产模拟方法 | |
Wu et al. | Reliability evaluation of the solar power system based on the Markov chain method | |
CN107562992A (zh) | 一种基于svm和粒子群算法的光伏列阵最大功率跟踪方法 | |
CN111523204B (zh) | 一种并网型综合能源网电-气储能***优化配置求解方法 | |
CN110795691A (zh) | 基于改进拟蒙特卡洛法的配电网风险评估方法 | |
CN111985805A (zh) | 综合能源***动态需求响应的方法及*** | |
CN107834540B (zh) | 一种基于概率约束的分布式光伏接入容量的确定方法 | |
CN105743093B (zh) | 一种考虑风光出力不确定性的配电网潮流计算方法 | |
CN105262146A (zh) | 含风电的电力***备用容量计算方法和*** | |
CN113609699A (zh) | 辐射状配电网交流潮流模型的计算方法和*** | |
CN111126772A (zh) | 一种模拟电力现货市场平衡的方法及*** | |
Bai et al. | Probabilistic production simulation of a wind/photovoltaic/energy storage hybrid power system based on sequence operation theory | |
CN115796721B (zh) | 高比例新能源接入的配电网运行状态智能感知方法及*** | |
Singh et al. | Reliability assessment of a remote hybrid renewable energy system using Monte Carlo simulation | |
CN113609686B (zh) | 一种新能源置信容量分析方法及*** | |
CN113591224A (zh) | 一种城市电网连锁故障风险评估方法及*** | |
Chenxu et al. | An improved cumulant method for probabilistic load flow calculation |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200508 |
|
RJ01 | Rejection of invention patent application after publication |