一种基于主从博弈的直购电双边交易的合同电价制定方法
【技术领域】
[0001] 本发明涉及合同电价制定方法，特别是涉及一种基于主从博弈的直购电双边交易 的合同电价制定方法。
【背景技术】
[0002] 随着中国电力改革的深入，电网垄断经销电量的格局正逐渐被打破。作为我国电 力市场化改革的重要组成部分，大用户与发电企业之间的直接交易是开展电力市场双边交 易的突破口。在电力市场环境下，现货市场的存在给市场参与者创造了获取较高利润的机 会，然而伴随着这种高利润的是高风险。而双边合同的持续时间一般较长，且合同价格在合 同签订的那一刻就决定了的，所以双边合同在能令合同签订的双方规避电价波动所带来的 风险的同时，也能保证用户的长期电力需求以及电力公司的长期负荷需求。协商式的双边 交易方式，作为目前国内大用户直接交易试点采用的比较通用的交易方式，具有充分体现 双方意愿，有助于买卖双方根据自身需要进行灵活的交易；买卖双方直接见面，交易简便易 行，不需要建立复杂的电子交易系统，技术条件要求较低等优点。可以预见，在未来开放程 度越来越高的电力市场中，双边交易将大量存在。
[0003]目前，国内外学者对于双边交易合同电价的研究，多是研究单一发电商与单一大 用户进行讨价还价式交易下的合同电价确定方法。鲜有文章对涉及多个发电商与多个大用 户的双边交易的合同电价进行研究。然而，在开放的电力市场环境下，大用户可以选择多个 发电商签订双边合同，而发电商也可以售电给多个大用户。
[0004] 由此可见，当前尚未有考虑多个发电商与多个大用户进行双边交易时的合同电价 制定方法。
【发明内容】
[0005] 本发明的目的是提出一种能够兼顾大用户和发电商双方利益的基于主从博弈的 直购电双边交易的合同电价制定方法。
[0006] 本发明所述的基于主从博弈的直购电双边交易的合同电价制定方法，包括如下步 骤：
[0007] Sl :基于主从博弈的直购电双边交易的合同电价决策的上层博弈模型
[0008] 在该博弈模型中，发电商有I个，大用户有J个，双边合同交易持续的时间为T ; 处于主从博弈上层的是发电商，发电商们在一个特定的计划范围内确定双边交易的合同电 价，在这个决策过程中，发电商的合同售电电量，由大用户根据发电商的合同报价来进行优 化决定；
[0009] 根据大用户优化合同购电量的结果，建立发电商i合同售电利润模型如下：
[0010] SI. 1:目标函数
[0012] 式⑴中，减号的左边表示的是发电商i通过合同售电得到的收益，减号右边代表 发电商i的发电成本，\ ,和b i分别为发电商i对大用户j合同报价曲线的截距和斜率，E u 为IX J的单位矩阵，
其中（为t时刻（t = 1，2,…，T)大用户j在 发电商i处的合同购电量，AjP B身别为发电商i发电成本的二次函数的一次项系数和二 次项系数；在该函数中，发电商通过修改自身合同报价曲线的截距来使得自己的合同售电 利润最大；
[0013] SI. 2 :发电商i合同售电利润模型达到最大化利润的约束条件为：
[0017] 其中qiniax为发电商i的发电上限；&和4分别为发电商i合同报价曲线的截距的 下限和上限；
[0018] S2 :大用户购电成本优化的下层规划模型
[0019] S2. 1 :大用户从发电商处的购电成本
[0020] 对于发电商i (i = 1，2,…，I)给大用户j的合同报价曲线
当大用 户j与发电商i签订的合同电量为
_.时，大用户j在所有发电商处的总购电成本为：
[0022] 其中，
为所有发电商对大用户j的合同报价曲线截距所组 成的矩阵；
为t时刻大用户j在所有发电商处合同购电量所组成的 矩阵；b = diagfti，b2，…，Idi]是一个对数矩阵，该矩阵的元素代表的是发电商对大用户合 同报价曲线的斜率；
[0023] S2. 2 :大用户从现货市场处的购电成本
[0024] 基于现货市场的预测电价M，当通过合同购电不能满足其用电需求时，大用户也 可以从现货市场买电，如果t时刻大用户j从现货市场购买的电量为，那么合同周期T 内，其从现货市场购电的总成本为：
[0026] S2. 3 :大用户购电成本建模
[0027] 大用户购电成本优化的下层规划模型可以表示为：
[0029] 约束条件：
[0033] 其中，式⑶是大用户的用电约束，C表示t时刻大用户j的需求电量，式（9)和 (10)为变量的非负约束；
[0034] S3 :直购电双边交易的合同电价决策的主从博弈模型
[0035] 综合发电商和大用户的最优规划模型，可得到主从博弈模型如下：
[0037] 约束条件：
[0041] 其中自变量qii j属于式（15)的函数的最优解集合：
[0043] 约束条件：
[0047] S4 :将合同电价决策的主从博弈模型转化为等价非合作博弈模型
[0048] 将含非合作博弈的合同电价决策的主从博弈模型转化成只含非合作博弈的模型， 具体步骤如下：
[0049] S4. 1 :利用大用户的KKT条件替代每个大用户的规划模型：
[0050] 设1^为大用户j在t时刻的购电成本优化的下层规划模型的拉格朗日函数，L ,的 表达式为：
[0052] 式（19)中，λ = [λ1; λ 2,…，A1], λ和u均为KKT条件的乘子，将Lj对式（19) 对所有的变量求偏导，就可以得到大用户j购电成本优化的KKT条件：
[0057] 解式（20)就可以得到购电成本优化后t时刻的大用户j的购电分布，从而得出大 用户在整个合同周期内的购电分布；
[0058] S4. 2 :发电商合同报价的非合作博弈模型
[0059] 将式（20)作为发电商合同售电利润模型的约束，得到直购电双边交易的合同电 价决策的主从博弈模型的等价非合作博弈模型如下：
[0061] 约束条件：
[0066] 与现有技术相比，本发明的有益效果在于：
[0067] 本发明将主从博弈的思想引入到大用户直购电的双边交易当中，根据主从博弈将 涉及多个发电商与多个大用户的双边交易的合同定价问题转化为发电商之间的非合作博 弈问题并进行求解，得到满足纳什均衡条件的发电商最优合同报价。大用户也可以按照本 模型，制定最优的购电分布使得自己的购电成本最小。
【附图说明】
[0068] 图1为本发明所涉及交易场景的决策关系图；
[0069] 图2为寻找本发明所建立的博弈模型的纳什均衡解的流程图；
[0070] 图3为发电商纳什均衡报价下的大用户的最优购电分布图。
[0071] 其中：图1中，氣』表不大用户j与发电商i签订的最优合同电量，f，.表不大用户 j与除发电商i以外其他所有发电商签订的最优合同电量，表示除大用户j以外其他所 有大用户与发电商i签订的最优合同电量，a;= (a aii2,…，Siil)表不发电商i对所有大 用户合同报价曲线的截距集合，a ;= (a a2,…，a; ai+1,…，aj表示除发电商i以外其他 所有发电商合同报价曲线的截距集合。a u表示除发电商i以外其他所有发电商对大用户 j的合同报价曲线的截距的集合，\ ,则表示发电商i对除了大用户j以外的其他所有大 用户的合同报价曲线的截距集合。
【具体实施方式】
[0072] 下面结合附图对本发明作更进一步的说明。
[0073] 本发明提出一种基于主从博弈的直购电双边交易的合同电价制定方法，为了兼顾 解决双边交易双方各自面临的问题，本发明采用主从博弈模型来研究发电商的最优合同报 价决策。本发明涉及的双边交易的决策关系如图1所示，首先发电商会给欲与其签订双边 合同的大用户一个合同报价曲线；然后，当大用户从所有的发电商处得到合同报价曲线之 后，在结合预测的合同周期内的现货电价以及自身的用电需求，来对自身的购电决策进行 优化使得购电成本最小，并将优化后的购电决策返回给发电商；最后，发电商根据大用户返 回来的购电决策来与大用户签订双边电力合同。
[0074] 在本发明的交易场景中，发电商的合同售电收入由与各大用户签订的合同电量和 合同电价决定，而大用户的购电成本则是由各发电商的合同报价曲线以及现货市场的电价 决定，因此，发电商与大用户之间是存在着主从博弈关系的。而发电商通过修改自己的合同 报价曲线，会影响到大用户的购电决策，进而会影响到其他发电商的合同售电利润，所以发 电商之间存在为竞争合同电量而产生的非合作博弈关系。
[0075] 本发明的方法包括：1)大用户根据各发电商的报价曲线，以及合同周期内的现货 市场的电价，并结合自身的用电需求来建立基于主从博弈的双边交易的下层优化模型；2) 发电商根据自身的合同报价曲线，发电成本函数，以及大用户的购电决策来建立基于主从 博弈的双边交易的上层博弈模型；3)综合发电商的博弈模型和大用户的优化模型，建立包 含多个发电商和多个大用户的两个利益相关方的完整的主从博弈模型；4)用大用户的购 电成本优化模型的KKT条件替换大用户购电成本的优化模型，并将该KKT条件作为发电商 售电利润优化模型的约束条件，从而将主从博弈模型转换为发电商之间的非合作博弈模 型，并对该非合作博弈模型进行求解得到各发电商纳什均衡条件下的最优合同报价。本发 明的具体步骤如下：
[0076] Sl :基于主从博弈的直购电双边交易的合同电价决策的上层博弈模型
[0077] 在该博弈模型中，发电商有I个，大用户有J个，双边合同交易持续的时间为T ; 处于主从博弈上层的是发电商，发电商们在一个特定的计划范围内确定双边交易的合同电 价，在这个决策过程中，发电商的合同售电电量，由大用户根据发电商的合同报价来进行优 化决定；
[0078] 根据大用户优化合同购电量的结果，建立发电商i合同售电利润模型如下：
[0079] SI. 1:目标函数
[0081] 式⑴中，减号的左边表示的是发电商i通过合同售电得到的收益，减号右边代表 发电商i的发电成本，\ ,和b i分别为发电商i对大用户j合同报价曲线的截距和斜率，E u 为IXJ的单位矩阵，
其中为t时刻（t = 1，2, "·，Τ)大用户j在 发电商i处的合同购电量，AjP B身别为发电商i发电成本的二次函数的一次项系数和二 次项系数；在该函数中，发电商通过修改自身合同报价曲线的截距来使得自己的合同售电 利润最大；
[0082] SI. 2 :发电商i合同售电利润模型达到最大化利润