第29卷第5期 2008年l0月 暨南大学学报(自然科学版) V01.29 No.5 Oct. 2oo8 Journal of Jinan University(Natural Science) 电力市场一类分段线性供应函数纳 什均衡解的存在性 周四清 (暨南大学经济学院,广东广州510632) [摘要] 供应函数均衡模型是研究电力市场发电商竟标行为的一种重要工具.不同于Coumot模型描述的产量 竞争和Bertrand模型描述的价格竞争,供应函数均衡模型更适合描述电力市场中竞标函数竞争行为.当电力市场 需求函数表示为市场价格与有界随机需求冲击的非线性函数时,每个发电商在其发电容量范围内向电力市场提交 一个分段线性供应函数,以获取最大化利润.在这种电力市场环境下构造了一类分段线性供应函数,证明了分段线 性供应函数纳什均衡解的存在性,从理论上获得了相关新结论,为电力市场供应函数均衡近似数值求解算法设计 提供理论基础. [关键词] 电力市场; 分段线性供应函数; 纳什均衡 [中图分类号】F416.61;F224.32 [文献标识码]A [文章编号】1000—9965(2008)05—0451—06 Research on existence of piecewise linear supply function equilibria in electricity markets ZHOU Si—qing (CoHege of Economics,Jinan University,Guangzhou,510632,China) [Abstract]The supply function equilibrium model has been widely used to study generators’bidding behavior in electricity wholesale markets over the past few years.In this paper there is a market where demand is uncertain,but can be represented as a nonlinear function of price together wiht a bounded rail— dom shock.Each generator offers a piecewise linear supply function,within his capacity of generation, into the market and wishes to maximize his expected profit.By constructing a different piecewise linear supply function,the existence of piecewise linear supply function equilibria in electricity market is proved.From the theoretical point of view,new results regarding existence of htese equilibria are presen— ted. [Key words]electricity market;piecewise linear supply function;Nash equilibrium [收稿日期]2008—05—23 [基金项目】 国家社科基金资助项目(06BJL066);暨南大学引进优秀人才科研启动基金资助项目(51205068) [作者简介]周四清(1964一),男,副教授,博士(后),主要研究方向:能源金融,智能信息处理 452 暨南大学学报(自然科学版) 第29卷 Klemperer&Meyer¨ 在分析英格兰和威尔士电 所有定理和结论的证明. 力现货市场交易竞争行为时提出了供应函数均衡概 念,首次求解了需求不确定情况下寡头电力市场供 应函数均衡解,并证明了只有当需求无限时才存在 唯一纳什均衡解.Green&Newbery 研究了发电机 1 问题描述 本文假设所有发电商成本函数是其发电量的二 次函数,发电商i的成本函数表示为c (q)=C g+ 组容量约束对电力市场供应函数均衡的影响,其容 量约束会大大降低纳什均衡解的个数;Green[3 将供 0.5D g2,其中ci,D 是给定的常数,q≥o表示发电 量,发电商i的发电容量为 ,i=1,…,n.在投 应函数模拟为线性供应函数,从此线性供应函数均 衡模型就成为了电力市场中一个研究热点;Baldick 等-4 研究了容量约束和限价约束对线性供应函数 的影响,证明了只有线性供应函数均衡点是稳定的, 特别对于具有非零截距线性边际成本的非对称发电 商,只有当市场价格大于其边际成本截距时才愿意 发电,由于不同发电商一般具有不同边际成本截距, 因此不同价格段内愿意发电的发电商组成与数目可 能很不相同,也就是不同价格段内可能存在不同线 性供应函数均衡.Green&Newbery 在应用供应函 数均衡模型研究英格兰和威尔士电力市场时,通过 分析每天提交的调度方案,发现需求的不确定性等 价于13时变需求.自此以后,供应函数均衡模型广泛 应用于电力现货市场发电商的竞标行为研 究 。 一 . Baldick&Hogan_5 研究了存在价格上限约束时 供应函数均衡点的存在性.Anderson&Hu_1。。从数 值近似求解角度提出了一种供应函数均衡解逼近方 法,通过构造一类分段线性供应函数来研究均衡解 的存在性,但均衡问题求解规模太大,同时在求解问 题的均衡解时,为了保证逼近解收敛于真实的供应 函数,需要验证全局最优. 本文研究电力市场供应函数均衡近似数值求解 的基础理论问题,当电力市场需求函数表示为市场 价格与有界随机需求冲击的非线性函数时,每个发 电商在其发电容量范围内向电力市场提交一个分段 线性供应函数,以获取最大化利润.在这一电力市场 环境下,本文构造了不同于文献[1O]的一类分段线 性供应函数,研究了分段线性供应函数纳什均衡解 的存在性,为电力市场供应函数均衡近似数值求解 算法设计提供理论基础. 第1节给出了问题的描述及相关记号,第2节 研究分段线性供应函数均衡解的存在性,附录给出 标决策过程中,发电商将成本函数表示为分段线性 函数,因而得出的边际成本也是阶梯函数(分段常 数函数). 假设需求函数形如D(P,8)=D(P)+ =A一 印+占,其中A,曰是给定常数,B表示发电商剩余需 求价格弹性系数,P是电力价格, 是需求随机需求 冲击变量且有界占≤ ≤否,满足已知累积分布F( ), 将随机需求冲击区间[ ,否]离散化: : <占:<… ,,. t、 <占 < K=否,其中K>1, =F l、 一上,l (k=l, …, ). 当价格上限为 时,发电商i必须提交所有可 用发电容量口i,供应函数s (q)在0<P≤卢时有定 义,特别当P= 时,8 (g)= ,(i=l,…,几),可进一 步假设 >max{c +D : 1,…,n},价格上限 严格大于所有发电商在其发电容量时边际成本.如 果存在某发电商 满足 <C + 及市场价格为 ,那么发电商为了自身的最大利益不会全力供电, 这样的话,发电商可能就会反对这种市场规则.事实 上,在实际运营的电力市场中,价格上限大于任何发 电商全力发电时的边际成本,如澳大利亚国家电力 市场,价格上限是10000澳元MWh;同时假设已有 容量满足价格响应需求,即D( , )≤∑ : . 给定整数K>1,发电商i的 对供应价格一电 量(pm,q )(k=1,…,K),其中P <P +t,0≤q*≤ q +,≤ ,后=1,…, ,构造分段线性供应函数如下: U’P≤pn且qil U q a-qg 1+il(P--Pil),p≤p 1且g l>0 p一pf1s (P)= q + (p-P ̄),’ (1) P ≤p≤p +1且1≤.i}≤ 一1 qil(+ +—P—p-Pile),,p≥p 且q p≥p 且q <qi一 iK-Pig-1 (p 第5期 周四清:电力市场一类分段线性供应函数纳什均衡解的存在性 453 图1表示了一个满足q 。>0,q <亏 的函数.由 (1)定义的函数正好连接了给定的价格一电量对, 拓广了价格定义区间[P¨P ],并且保证发电量不 为负或大于发电容量.当g >0及g*< 时,由(1) 定义的函数s (P)满足:存在6>0,当P∈[p 一6,p + ]时,s (P)>0及s (P)<口 .当其它发电商响应 发电商的供应函数时,发电商i供应函数仍然趋于 P 和P 由(1)式定义的函数s (p)在 处不光滑, k:1,…。 图1 连接最优响应点的逼近方法 图2最优响应点光滑的逼近方法 图2是文献[10]中函数s (P)的一种构造方 法,分段线性函数s (p)在p 处是光滑的,其中 是 相邻两分段的交点. 本文将得出与图1和图2所定义两种逼近方法 完全不同的供应函数均衡解. 给定价格P 和随机需求冲击占 ,发电商i的利 润为兀;(Pm,F ):p q —C (q ),其中q =O(p , s )一∑ sj(p )是竞争对手留给发电商i的剩余 供应电量.如文献[11]所建立的供应函数均衡框 架,假设发电商i已知竞争对手的供应函数s (p), ≠ ,选择价格P ≤ 使得它的利润兀i(P )最大, 并在此价格下满足剩余需求(i=1,…,n),也就是对 于发电商 ,求解下面利润最大化问题 max (P*,占 ) pit s.t.P ≤ 且 0≤D(p 占 )一∑ (p )≤ (2) 其中 是随机需求冲击,s (p)由(1)式定义.注意 到(2)是一维最大化问题,并且分段光滑.此外,在 (2)中所有函数都有左导数和右导数.所以,P (不 是区间的端点)是(2)的局部解当且仅当 兀 +(p )≤O和兀 i一(p )≤O. 下一节将在上述假设条件下研究分段线性供应 函数纳什均衡的存在性. 2分段线性供应函数纳什均衡解的存 在性 假设随机需求冲击区间[ , ]的一个离散化样 本为 =( 一,占 ),其中 曼= 1<占2<…< JI(一1<l?K=否. 分段线性供应函数3i (P)集:由(1)式定义并 满足(p主,g主)(i=1,…,n,k=1,…, )是 的一个 纳什均衡,如果对每个随机需求冲击 ,当i=l, …,n时,所有价格p萎都是相同的,记为p ,对于i= l,…,,l和给定的s (p)(J≠ ),Pk‘是问题(2)的最 优解.如果没有特别说明,对于 =1,…,n和k=l, …, ,价格一电量对(p ,g )在p主<P :+ ≤ 和 O≤q主≤q + ≤口下是有意义的. 在讨论由最大化问题(2)定义一阶优化条件之 前,先导出一个条件意义下优化模型: ≤元,如果 <元,则Y=0;如果 = ,则Y=1.等价于下列优化问 题的解: maxY S.t.0≤y≤1,Y(元一 )=0 其中 ≤元,一阶优化条件变成 l+ (元一 )+ 1一 2=0,0≤),≤1,t,1≥O, 21>0, 1Y=0, 2(1一Y)=0 其中 , , :分别是对应约束的拉格朗日乘子. 考虑下面3种情形, 第1种情形:pK≤ ,如果P = ,则q =q (i= 1,…,n). 第2种情形:q ≥0,如果q =0,则s (P )=0. 暨南大学学报(自然科学版) 第29卷 . 第3种情形:类似于第2种情形,但k=K当时, q =茸 . 明这种优化曲线是存在的.由文献[7]中定理l4有 以下结论成立. 又因为兀 +(P , ):q +[P 一C (‰)]・ 性质2 假设n个发电商有相同的凸成本函 【D ( )一∑ s +(p )],兀 一(p , )有类似的 数c(・),则供应函数均衡存在,每个参与者提交供 表达式,所以最大化问题(2)在P 处优化条件可以 写成(3): q 一(p —C —D 【B+∑ ( + ~qjk)/ (P +i—P )J+ 一 m≤0 一q +(p 一C —Diq )【B+∑ ( 一 )/ (P 一P 一1)l— + ≤0 ∑ :。 =A+占 一8p (q —q )/(p 一P )≥0 ≤P≤ ,0≤q*≤ ≥0, ≥0 (亏;一q )=0, q =0 i=1,…, ,k=1,…,K (3) 其中q =0,q = .由(3)中第1个不等式得出方 程组的解P ≥C ,因为q ≥0,B+∑ ( +。一 q )/(p +1一P )>0,此夕 ,Pl<…<P < . 方程组(2)有 +, 未知数:p ,q i=1,…,n, k=1,…, 对于(3)的解,由占…> 有P +1>P 及q + ≥‰事实上,如果存在某发电商i和随机需 求冲击 使得q肌 <q ,则P…<P^以及对任意 ≠ , +1≤ ,因为(qm+1一q )/(p +1一Pk)>t0.但 由于D( + )=A+ 一印 =∑ + < ∑ =A+s 一印 =O(p , ),不可能推出 >s 及P +1<P . 定理1假设(3)的解为(p ,q*)(i=1,…,n,k =l,…, ),则 (1)每个发电商都有一个来自竞争对手的凸聚 集供应函数,它的利润是凹的. (2)其解是供应函数竞标博弈的纳什均衡. 证明见附录. 当只有两个竞标者或者发电商对称时,定理1 的结论表明任何均衡供应函数都是凸的,下面将证 应函数s(・)= 一 (・),其中r(q)=q I口+(凡 L ’ 1 一I)f[cq +( )/x ] l,J nq 是最大需求,任意选 择的常数0满足口≥c +(q )(q )卜 . 于是得到下面性质3 性质3假设n个发电商的成本函数都为cg+ 0.5Dq ,s(p)和口分别是性质2中的均衡供应函数 和常数,则下面结论成立: (1)当 =2时,供应函数 (・)在其定义域内 是凸的. (2)当n>2且口≤c (g )(q ) 一 +D (q ) /(n一2)时,供应函数s(・)在其定义域内 是凸的. (3)当n>2且0>c (q )(q )卜 +D (q ) /(n一2)时,供应函数 (・)在其定义域内 是凹的. 证明见附录. 如果存在光滑供应函数均衡使得每个发电商有 一个凸聚集供应函数,则当在端点 和 上合理地 选择供应函数的值时,供应曲线上的点就构成竞标 博弈{(2)} i=1的纳什均衡.即有下面定理4成立. 定理4设{s (P) 是区间[ ,卢]上的光滑 供应函数均衡,n个发电商的光滑凸成本函数为c (q),使得∑ sj (p)是p的凸函数, =1,…,凡.对 于任意给定有穷个点P 满足笆≤_p <P <…<pK≤ 卢,令需求函数D(p)满足D (P)<0,对于随机需求 冲击区问内的点 。,…, 满足∑ : s (p )= O(p )+s ,则{P1,s (P ),…,P ,5 (PK)} :l是满 足(3)韵解,特别当P =卫或P =卢时,可以忽略在巳 或 时右导数条件或左导数条件. 证明见附录 从定理1和性质3,(3)并没有包含当发电商对 称时供应函数均衡是凹的情形.下面定义另一方程 组处理此种情况,其中并不直接处理随机需求冲击, 而是用任意两个相邻随机需求冲击的中点来定义 第5期 周四清:电力市场一类分段线性供应函数纳什均衡解的存在性 455 Q 一( —C —Diq )f B+∑ ( + 一 )/ 1 s ’(P )} : 是满足(4)的解. 证明见附录 (P 一P )f>0 J —3小结 当利用供应函数均衡模型来研究发电商竞标策 q*+(p 一C —Diq )f B+∑ ( 一 一 )/ 1 (P 一P )I>0 J ∑ =A+(占 + )/2一却 (9 —g 一1)/(p 一P 一1)≥0 ≤P ≤ ,0≤g*+( +1一s )/2≤ f i=1,…,n,k=l,…,K (4) 其中 Q =A+% 一却 一∑ =q +( , 一 )/2 i=1,…,凡,k=1,…, 由(4)中第一个严格不等式表明,发电商的利 润函数随市场价格P严格递增,其中占 是随机需 求冲击,聚集供应函数是 ∑ [ +(p—p )( + 一 )/(p 一Pk)], P∈[ ,P ] 而(4)中第二个严格不等式表明,发电商的利 润函数随价格P严格递减, 是随机需求冲击,聚 集供应函数是 ∑ [ 一 +(p—Pk-1)( 一 一。)/( — P¨)],P∈[ ,P ] 所以当随机需求冲击为 ㈨时,发电商的最优 价格P∈(P ,P ),因为发电商i的利润函数在区间 [P ,P川]内是凹的.当 一 , 一 时, 是随机 需求冲击,所有发电商有相同的最优价格. 类似于定理4,有如下定理成立. 定理5设{s,(P)} : 是区间[ , ]上的光滑 供应函数均衡,n个发电商的光滑凸成本函数为c (g),且满足∑ (p)是p的凹函数,i=1,…,n. 对于任意给定有穷多个点P。满 ≤P。<P:<… <P ≤ ,令需求函数D(p)满足D (P)<0,随机需 求冲击区间内的点 。。, 满足∑ : s (P )= D(p )+( +占 +1)/2,贝U{P1,s (P1),…,P , 略行为时,对于需求不确定的电力市场,如果市场需 求可以表示为市场价格和随机需求冲击的非线性函 数,本文通过构造不同于文献[10]的一类分段线性 供应函数,从理论上证明了分段线性供应函数均衡 解的存在性,下一步证明这一类分段线性供应函数 均衡解的收敛性,这些理论结果将为设计电力市场 分段线性供应函数数值求解方法提供理论基础. 4 附录 定理1的证明:对于i=1,…,n和k=1,…, , 由(3)中前两个不等式得出 ∑ ( + 一 )/(p 一p )>∑ ( 一 g —1)/(p 一P 一1) 即每个发电商在(3)的解将有一个聚集凸上升 分段线性供应曲线,发电商有全局拟凸利润函数,从 而解就是纳什均衡.并由上述表达式发现,两个端点 的选择影响解的存在性.证毕. 性质3的证明:由C(q)=Cq+0.5Dq 直接计 算得, 当n=2时,T(q)=C+g(n—C/q +Dlnq )一 Dqlnq 当凡>2时, T(q)=[(n一2)C+(n一1)Dg]/(n一2)+9 一 [n一[(n一2)C+(n一1)Dg ](g ) 一 /(n一2)] 当0<q<q 时,由性质2得 (q)>0.计算 (・)的二阶导数得, 当n=2时, (g)=一D/q 当凡>2时, (q)=(n一1)(,£一2)q 一 [a—c (g )(q ) 一“——D(q ) 一“/(,z——2)] 又s(p)=T一(P),得5”(p)= (s(P))/[T (s (P))] ,所以当口≥c (q )(g )卜 和0<q<q 时, 得出 ”(q)= ,>0, 口≥c (q )(q ) ,凡=2 {≥o,口≤c (q )(9 ) 一 +D(g ) 一“/(n一2),n>2 【<0,n>Ct(g )(q’) 一 +D(q ) 一 /(n-2),凡>2 456 暨南大学学报(自然科学版) 第29卷 由上式知,定理结论成立.证毕. 定理4的证明:首先,对于k=2,…,K,有 ∑ ( )一 ( )]/(p )]≥0 又Q“ =q +(8 一 )/2>q ,(4)的第1个严 格不等式成立.类似可以证明(4)式的第2个严格 不等式.证毕. [参考文献] ∑ ( )≥∑ ( ) ̄-ID(p )<D(p ), < .其次,对每个i,P是 从而有 1<占2<…< 下列最大化凸问题的最优解 m a x (D(p)+ 一∑, (p))p—c (D(p)+ s一∑ (p)) 对任意 和最优解条件,s。 (P)=一[P— c (s (p))](D 一∑ sj (p)),从而推出p≥ c (s,(p)),如果s (p)>o,由D 一∑ (p) <0知严格不等式成立.又∑ 々 (p)是凸的,对 任意P‘<P…有 ∑ (p )≤ ∑ [ (p ) (p )]/( 一p )≤∑ si ( ) (3)式在p 关于q = (p )的右导数条件 s (p )一[p 一c,。(q )]_Dip+ ∑, [ (p + )一s/(p )]/(p“。一p )]= [p —c’ (g )]I∑ (p )一 ∑ [ ( +t)一 ( )t/(p 一 )l≤o 类似地可证明左导数条件成立.证毕. 定理5的证明:首先,类似于定理4的证明,有 81<s2 <… < 一1 < 和s (p) =一[p— (s『(p))](D'F一∑ (p)),p≥ ( (p))(如 果si’(p)>o,则严格不等式成立).因为 ∑ (p)是凹的,对任意 <p㈧,有 ∑ (p )≥∑』 [ (p )一s,.(p )]/ (p 一pk)≥∑ ( ) 从而 si (p )一[p 一c (g )][~D + ∑ [ )一 ( )]/( 一Pk)J= [p (q )]【∑ (p )一 [1] KLEMPERER P D,MEYER M A.Supply function equi- libria in oligopoly under uncertainty[J].Econometirca, 1989,57(6):1243—1277. [2]GREEN R J,NEWBERY D M.Competition in the briitsh electricity spot market[J].Journal of Political Economy, 1992,100(5):929—953. [3] GREEN R J.Increasing competition in the british elec・ tricity spot market[J].Journal of Industiral Economics, 1996,44(2):205—216. [4]BALDICK R,GRANT R,KAHN E.Linear supply fune- tion equilibrimn:generalizations,application,and limita— tions[EB/OL].http://www.ucei.berkeley.edu/ucei/ pwrpubs/pwp078.html,2000. [5]HOGAN W.Capacity constrianed supply funciton equi- librium models of electricity markets:stability,nonde— creasing constraints,and fnuction space iterations[R]. Paper PWP-089,University of California Energy Institu・ te,2001. [6] BALDICK R.Electriciyt market equiilbrium models:The efect of parametirzation[J].IEEE Transactions on Pow・ er System,2002,17(4):1170—1176. [7] RUDKEVICH A.Supply function equiilbrium in Poolco— ytpe power markets:learning all the way[R].Technical Report TCA 0699—1701,Tabors,Caramanis,and Asso- ciates,1999. [8] WOLAK F.An empiircal analysis of impact of hedge con— tracts on bidding behavior in a competitive electirciyt mal'- ket[R].NBER Working Paper No.8212,April 2001. [9] BORENSTEIN S.The trouble with electricity markets: nuderstanding California’s restructuring disaster[J]. Journal of Economic Perspectives,2002,16:191—212. [10]ANDERSON E,HU X.Computing supply function equi・ libria in electriciyt markets[R].Austrlaina Graduate School of Management,University of New South Wales, Sydney,2005. [11]HOLMBERG P.Asymmetirc supply function equililibri— um with constant marginal costs[R].Department of Eco- nomics,Uppsala Universitet,Sweden,2005. [责任编辑:王景周]