0 引 言

以风电、光电为代表的可再生能源具有无污染物排放、不消耗能源的优势，越来越多地渗透到电力系统中。但是风电和光电出力主要受自然条件的影响，具有较强的间歇性和随机性，电力系统的功率平衡受到了干扰。传统的火电机组在发电时会产生以二氧化碳为代表的污染物，严重影响了环境。将传统的火电机组加以改造，引入碳捕集系统，形成碳捕集电厂(Carbon Capture Plant),对产生的二氧化碳进行捕捉和封存，是减少污染物排放的有效有段。碳捕集机组的收益直接受碳排放权价格和电价的影响，在以往的研究中将碳排放权价格和电价视为固定的常数，而在实际电力市场中，碳排放权价格和电价都具有一定的波动性。电动汽车清洁性高，V2G技术为车主和电网建立了双向互动的桥梁，电动汽车具有了充电负荷和放电电源的双重角色，可以促进可再生能源的并网。采用静态分时电价虽然可以引导电动汽车改善负荷曲线，但是有可能产生新的高峰。因此，如何制定合理的电价机制，研究多种随机变量的联合作用对电力系统生产和调度的影响，是非常值得进一步研究的问题。

对于分时电价下的电动汽车和随机因素对电力系统的影响问题，相关学者已经取得了一定的研究成果。文献[1]用k均值聚类算法对时段进行划分，该分类方法属于硬分类，对于无法准确确定归属的数据难以处理。文献[2-3]用基于模糊关系的传递包方法对时段进行划分。文献[4]提出了机会约束方法，其应对随机变量时具有良好的处理能力，但是其求解过程较为复杂。文献[5]提出了条件风险理论来处理不确定因素对经济性的影响，但是无法验证该理论是否适用于电力系统领域。文献[6]考虑了负荷、发电机和网络结构不确定性对结果的影响，但是侧重于考虑随机因素本身的影响。文献[7]将超分位数方法应用到碳捕集机组场内优化运行当中，但是该文献侧重于对单个随机变量的描述，而现有的电力系统问题中，随机变量往往并不是单一的，在应对多种随机变量的联合作用时可能无法处理。文献[8]提出了考虑风、光、荷联合作用的优化调度模型，但是没有考虑电价的波动性和碳排放权交易的背景。

针对上述研究的不足，提出了基于动态分时电价并考虑多种随机因素的两阶段联合优化调度模型，分别建立了风电、光电、碳排放权价格和电价的概率分布模型，用模糊聚类分析解决峰谷时段划分，用推广的超分位数方法处理多重随机因素的联合作用，并对相关的参数进行敏感度分析，改变参数的取值以分析其对结果的影响。

黄春贵先生说：“大约用典之佳者，贵能推陈出新，无异于出自一己之创作，譬如水中煮盐，运化无迹，不使人觉。……故原本古事成辞，用典时却须重加铸造，别出心裁。……夫善纫者无隙缝，工绘者无渍痕，用典若斯，紧著题意，融化而不涩，用事而不为事使，则面目精神，方能一新。”[3]沈约虽未就用典作进一步发挥，然在其后期诗作中出现的一些作品，却在某种程度上将用典艺术推向这一高度。

1 模糊C均值聚类

聚类分析通过一定的标准将一组对象分为若干类，每个类别内部元素有一定的相似性。最初的聚类分析是硬分类方法，对事物的处理非此即彼。然而有的对象无法准确判断其归属，将模糊集理论引入到聚类分析当中，可以更好地处理数据的归属问题。模糊C均值聚类是一种基于目标函数优化的模糊聚类分析方法，用隶属度表征对象的相对归属性。因其对初始聚类中心要求不高，处理复杂分类问题时更为合理而得到了广泛应用。首先建立规划模型：

采用营养钵育苗，营养钵直径8厘米、高10厘米，每钵播1粒种子，播后盖上厚0.5～1厘米的基质，每亩大田需备营养钵约2500个。当出苗率达50%时，春季及时揭除营养钵上塑料薄膜，夏秋季通过开关通风口、开合遮阳网，合理调节棚室温度，保持23℃左右。及时做好调水、通风、换气工作。

(1)

(2)

(3)

用拉格朗日乘数法处理目标函数得：

(4)

(5)

(6)

(7)

式中n为数据的个数；c为分类数；uik为第k个数据隶属于i类的程度；dik为数据k与聚类中心i的距离；A为对称正定矩阵；m为模糊加权指数。

模糊C均值聚类算法的具体流程为：(1)初始化。确定模糊加权指数、分类数、迭代停止阈值和迭代次数，随机生成初始隶属度矩阵；(2)按照式(6)修正聚类中心；(3)按照式(5)更新隶属度矩阵。按照上述步骤反复迭代直到满足收敛条件；(4)收敛条件判定(b为迭代次数，ε为迭代停止阈值)。

2 推广的α超分位数方法

超分位数方法是Rockafellar R.T.和Royset J.O.在2009年提出的一种可以有效处理随机问题的方法[9]。它满足一致性公理、具有次可加性和凸性，解决了α分位数的缺陷。它代表的是一种平均水平，可以精确衡量潜在的风险，能将含随机变量的复杂规划问题进行简化求解。α超分位数定义为：状态函数不低于分位数条件下的期望值。超分位数方法刚提出时用于解决单个随机变量问题，现在将该方法推广以处理多个随机变量问题。

根据《布鲁塞尔条例Ⅱa》第15条不方便法院原则的特殊立法规定和欧洲法院在Child and Family Agency案对不方便法院原则的特殊适用，欧盟家事诉讼中不方便法院原则的进步性与合理性在于：

设决策变量x∈X⊂Rn，随机变量w∈Rm，v∈Rl，u∈Rs，y∈Rk，设状态函数为g(x,w,u,v,y)，w、v、u、y的联合概率密度函数为ρ(w,u,v,y)，则状态函数不超过阈值β的分布函数为：

碳捕集电厂在原有火电厂的基础上引入碳捕集系统，碳捕集系统能够对产生的CO2进行捕捉和封存，等效于减少了系统的CO2排放量[12]。碳捕集电厂的数学模型为：

生19：老师，其实不一定要是直角，只要∠C1AB=∠C2BA=∠EDF就有△AED∽△BDF，只要当∠C1AB=∠C2BA=∠EDF ，D是AB中点就有△AED∽△DEF∽△BDF．

ATB混合料碾压应遵循“紧跟慢压、高频低幅”原则，压路机应紧跟摊铺机后进行碾压，避免混合料热量散失过多，初压使用13T双钢轮振动压路机静压2遍；复压采用13T振动压路机振动碾压3～4遍，碾压温度控制为100～110℃；终压使用11T钢轮压路机静压2～3遍以消除轮印，碾压温度不低于80℃。

g(x,w,v,u,y)≤β

(8)

设置信水平为α，则分布函数的α分位数为：

qα(x)=min{β∈R,ψ(x,β)≥α}

(9)

α超分位数定义为超过α分位数条件下状态函数的期望值：

=∬∬g(x,w,v,u,y)ρ(w,v,u,y)dwdvdudy

(10)

直接根据解析式求解比较困难，通过引入中间变量简化问题求解：

z0}ρ(w,v,u,y)dwdvdudy

(11)

式中的积分项难以计算，对式(11)进行转化：

1936年，美国参议院颁发了第264号文件。文件中包括这样一段警告：“一个惊人的事实是，在数百万英亩的土地上再也没有足够的矿物质，而从这些土地上出产的食物(水果、蔬菜和粮食)正在让我们变得饥饿——无论我们吃多少东西都难以消除这种饥饿。今天的人们永远不可能摄取足够量的蔬菜和水果，为他的身体系统提供塑造完美体格所需要的矿物质，因为他的胃不可能装得下。”

(12)

通过样本均值对上式进行估计，采用蒙特卡洛随机抽样，设各随机变量的样本分别为(w1,w2,...,wm)，(v1,v2,...,vl)，(u1,u2,...,us)，(y1,y2,...,yk)，则有：

（4）调控微生物群落的原理。高比例的产酸菌极大地促进了废弃矿业地的酸化，并且容易导致返酸现象的发生。通过添加石灰急剧改变废弃矿业地的pH，由于调节后的pH超出这些产酸菌的生长pH范围，因而产酸菌比例可以得到大幅降低。通过添加的微生物菌剂，引入与培育有益微生物，进一步降低产酸菌比例，发挥有益微生物的作用持续改良土壤。

(13)

由于式(13)中含有非光滑函数，仍然难以直接求解，再次引入中间变量：

z=max{0,g(x,wj,vr,uo,yq)-z0}

(14)

对式(14)中的等式约束进行松弛处理，可以得到最终的优化模型：

z≥g(Pevt,Pwtj,Pstr,Cto,Stq)-z0

2007—2017年世界纸和纸板出口额排名前5位的国家包括加拿大、芬兰、中国、德国、加蓬、英国、斐济和瑞典，2007年依次为加拿大、芬兰、德国、加蓬和英国，2017年为中国、芬兰、德国、英国和加蓬。各国世界占比除2007年最高者加拿大为10%、最低者英国为5%以及2007年最高者加拿大为9%、最低者斐济为5%外，其余年份均差别不大，前5位国家均在5%～8%区间。

o=1,2,...,s; q=1,2,...,k

(15)

3 两阶段联合优化调度模型

3.1 风速分布函数

风电出力大小受风速直接影响，相关研究表明，风速服从Weibull分布、正态分布或Rayleigh分布。由于后者使用的更多，文中假设风速服从Rayleigh分布[10]：

(16)

式中v为风速；σw为分布参数。

z≥0,j=1,2,...,m; r=1,2,...,l;

风电出力与风速的关系可以由分段函数来表示：

(17)

(18)

电动汽车在时间段t不仅要与电网互动，自身的行驶还需要消耗一部分能量，其在相邻时间段的电池容量关系可以用下式来描述：

3.2 光照强度分布函数

光电出力大小直接受光照强度影响，相关研究表明，某时刻光照强度服从Beta分布[11]。光照强度概率密度函数为：

(19)

式中 为时段t最大光照强度；α、 β为分布参数。

光电出力与光照调度的关系可以表示为：

Pst=rtηsA

(20)

式中ηs为光电转换效率；A为光电电池方阵的总面积。

由于光电出力与光照强度呈线性关系，因此光电在某一时段出力也服从Beta分布：

(21)

3.3 碳排放权价格分布函数

相关的碳政策规定当系统的CO2排放量超标时，会给予对应的经济惩罚，而排放量较低时，可以将多余的额度拿到电力市场出售。每个发电系统有自己的碳排放额度，含碳捕集电厂系统的碳排放权收益可以用以下模型描述：

(22)

(23)

Ei(PGit)=δiPGit

(24)

(25)

(26)

(27)

式中T为调度周期的时段数目；ηc为单位发电量碳排放限额；PDt为单位时段内系统总的发电量；N为常规机组的数目；M为碳捕集电厂机组的数目；Ei(PGit)为单位时段内常规机组i的碳排放量；为单位时段碳捕集机组j的净出力；为单位时段内碳捕集机组j的碳排放量；PGit为单位时段常规机组i的发电量；PCjt为单位时段碳捕集电厂机组j的实际发电量；Ct为当前的碳排放权价格；Ep为调度周期内总的CO2排放量；δi和δj分别为常规机组i和碳捕集机组j单位发电量的碳排放强度；F1为碳排放惩罚成本。

ψ(x,β)=∬∬ρ(w,v,u,y)dwdvdudy,

E′=E-Ec

(1)太行山区矿山企业数量和从业人员数多。2017年河北省持证矿山企业数量太行山区1 030个、冀北729个、冀东584个，分别占矿山企业总数的33.59%、23.78%及19.05%；按从业人员数量排序，位居第一的是太行山区7.82万人，其次是冀东5.75万人， 分别占矿山企业从业人员总数的41.59%和30.6%。统计数据显示，太行山区具有矿山企业数量和从业人员数多的特点。

(28)

g′=g-gN

(29)

Ec=γE

(30)

E=δg

(31)

gN=αEc

式中FR和FC分别为常规机组和碳捕集机组的发电成本；ai、bi、ci和aj、bj、cj分别为常规机组和碳捕集机组的发电成本系数；Fw和Fs分别为风电和光电的惩罚成本；Cw和Cs分别为风电和光电的惩罚成本系数； 和分别为风电和光电在时段t的出力预测值；Fev为电动汽车充放电成本；St为时段t的电价。

(32)

(33)

式中E′为机组的等效CO2排放量；E为机组发电总的CO2产生量；Ec为机组的CO2捕集量；g′为机组的净出力；g为机组的实际出力；gN为用于捕捉CO2而损耗的出力； γ为机组的CO2捕集率，其取值一般在85%～95%之间[12]；δ为机组的单位电量碳排放强度；α为捕集单位CO2所损耗的能量。

碳排放权价格受许多因素的影响，波动性较大，相关研究用AR-GARCH模型来描述其分布规律[13]：

(34)

式中Ct为第t日的碳排放权价格；Rt为收益率；ut是服从正态分布的随机变量，均值为0,方差为有限值；δt、σt为中间变量，其它参数均为常数。

若已知t-1,t-2,t-3日的碳排放权价格，可以根据上式求出δt-1和σt-1，进一步可以求出t日的碳排放权价格：

(35)

由于式中存在随机变量，这就决定了碳排放权价格也是按照某种规律分布，具有一定的波动性，设Ct和ut满足如下的函数关系：

Ct=h(ut)

(36)

3.4 电价分布函数

电力市场受随机因素的影响，电价难以准确确定。原有的研究都假设电价服从正态分布，但是近年来越来越多的研究发现电价的分布呈现“尖峰厚尾”的特点，这表明电价并不是服从简单的正态分布，但其本质上还是一种正态型曲线，此处用加权双高斯模型来更加精确地刻画电价分布特性[14]：

过了好久，竹韵的泪水才涌出眼眶，她不知道龙斌为什么要自杀，八年多过来了，夫妻俩过得好好的，他会有什么想不通的呢？她感到很累很无助，听着海力充满温情的劝慰，抬头看了一眼他焦急而又关切的面容，突然感到他是那么亲切，那么令她感动，如果她不是个有夫之妇，真想扑进他怀里痛哭一场，然后靠在他厚实宽大的肩膀上睡上一觉，跟他在一起总让人感到安全踏实……

(37)

式中w1、w2为权系数；μ1、μ2分别为两个正态分布的期望；σ1、σ2分别为两个正态分布的标准差。

3.5 两阶段优化调度模型

风电和光电有较强的波动性，电动汽车入网可以协调这些能源的并网，并且可以与电网双向互动。电动汽车在一天中有90%的时间处于停驶状态[15]，可以作为电源被电网调度。在调度周期内的t时段可供调度的电动汽车数量Navt可以表示为：

Navt=Nevpparkt

(38)

式中Nev为电动汽车的总数量；pparkt为时段t电动汽车停驶的概率。

式中vt 、vi、 vr 、vo分别为t时段的风速、切入风速、额定风速和切出风速；PN为风机额定功率。

(39)

式中ηev为电动汽车充放电效率；Pevt为电动汽车的调度计划出力，若为正值表示电动汽车作为电源向电网放电，若为负值说明电动汽车作为负荷充电；Δt为单位时间长度；Qkm为电动汽车单位距离所消耗的电能；为电动汽车平均行驶速率。

通过价格机制引导电动汽车可以促进削峰填谷，但是现有电价机制难以反映实际供需状况，静态分时电价可能会产生新的负荷高峰。在这里采取动态分时电价机制，动态分时电价是指高峰、低谷电价固定，而峰谷时段变化的电价机制[16]。第一阶段根据负荷曲线和风电、光电出力的预测值，按照第1节提出的方法对峰谷时段进行有效划分，其中分类数为3类，聚类中心分别取负荷最高峰、最低谷和居中的负荷值。

在此基础上，第二阶段对多重随机因素的调度问题进行优化。碳捕集机组的成本和常规机组的成本类似。风电和光电出力具有不确定性，当实际出力大于计划出力时，造成弃光、弃风，实际出力小于计划出力，需要增加额外的备用容量，对该偏差引入惩罚系数来衡量不确定因素造成的损失。电动汽车的充放电成本值为正时，表示电动汽车向电网放电，电网需要向车主支付相应的费用。电网总成本的目标函数可以表示为：

(40)

(41)

(42)

(43)

(44)

Fev=StPevt

(45)

In the RF measurement system,following the usual signalflows in a digital architecture of the Intermediate Frequency(IF)receiver,the incoming signal composed of direct and reflected signals in the presence of the multipath effect can be expressed as:15,18

通过多年来的掘进实践调查，复杂地质对煤巷机械化掘进效能的影响因素是多方面的，且各因素间相对独立又相互影响。不同的地区地质差异很大，须分析具体地质状况来采取预防措施。合理的技术措施可保证煤矿的安全掘进，并提高采掘效率。因此，如何做好煤矿采掘过程的支护工作，确保巷道的畅通运行，保证围岩的稳定是我们面临的重要问题。

由上述分析可以得到风电功率与风速、光电出力与光照强度、碳排放权价格与ut以及电价本身之间的关系：

(46)

其逆变换的雅克比行列式为：

若风速、光照强度、ut和电价的联合概率密度为f(vt,rt,ut,St),则风电出力、光电出力、碳排放权价格和电价的联合概率密度为：

(47)

现用推广的α超分位数方法对上述模型进行变换，状态函数为：

g(Pevt,Pwt,Pst,Ct,St)=Fw+Fs+Fev+F1

小径薄壁管材内表面检测，传统方式是采用目视和聚光灯配合，有异样状况时，则会采用工业视频内窥镜来检查其缺陷状况。我单位最常使用的为IPLEX LX系列，其焦距为 4.0～6.0m，镜头直径为 4—8.5mm，探头软管长度为3.5m，采用USB接口，CD成像技术。光纤视频内窥镜价格昂贵，通过探头拍摄到管子内部情况在液晶显示器上成像，由此可初步判定管材内表面异常情况，为技术人员和质量检验人员体统分析依据。如图1、图2所示。

(48)

状态函数不超过阈值β的分布函数为：

ψ(Pevt,β)=∬∬ρ(Pwt,Pst,Ct,St)dPwtdPstdCtdSt,

g(Pevt,Pwt,Pst,Ct,St)≤β

(49)

设α为给定的置信水平，则分位数为：

qα(Pevt ) =min{β∈R,ψ(Pevt,β)≥α}

(50)

其超分位数刻画为：

∬∬ g(Pevt,Pwt,Pst,Ct,St)ρ(Pwt,Pst ,Ct,St)dPwt dPstdCt dSt

(51)

引入中间变量：

(52)

为了对上述式子进行估计，再引入中间变量，其中Pwtj、Pstr、Cto、Stq为各随机变量的样本。

z=max{0,g(x,Pwtj,Pstr,Cto,Stq)-z0}

(53)

松弛上述等式约束，目标函数变为：

(54)

约束条件为：

(1)系统约束

(55)

(2)常规机组和碳捕集机组约束

(56)

(3)电动汽车

(57)

(4)超分位数约束

z≥g(Pevt,Pwtj,Pstr,Cto,Stq)-z0

z≥0,j=1,2,...,m;r=1,2,...,l;

o=1,2,...,s;q=1,2,...,k

(58)

式中PLt为t时段总的负荷量；PGimin、PGimax分别为常规机组i的出力下限和上限；rGiup、rGidown分别为机组i出力上升速率和下降速率上限；PCjmin、PCjmax分别为碳捕集机组j的出力下限和上限；rCjup、rCjdown分别为碳捕集机组j出力上升速率和下降速率上限；Ru、Rd分别为系统的正、负旋转备用率；αev和βev分别为电动汽车电池容量下限和上限系数；qevmax为电动汽车电池最大容量；pc、pd分别为电动汽车充放电功率。

4 算例分析

4.1 算例系统参数

采用6台常规机组[17]、10 000辆电动汽车、10单位风机和光伏发电系统参与调度的算例来进行分析，调度周期为一天24 h，将机组6改造为碳捕集机组。电动汽车参数、停驶概率和负荷参数如表1、图1、表2所示。风电机组参数、光伏发电系统参数、分布参数、风速和光照强度参数参考文献[18]。α取0.31(MW·h)/t，η取0.798， γ取85%，正负旋转备用率取5%，风电和光电的惩罚系数取200元/(MWh)，待优化日前三天的碳排放权价格依次为140元/t、105元/t、119元/t,基础置信水平取0.90。

表1 电动汽车参数 Tab.1 Parameters of electric vehicle

参数数值电池最大储存容量(kWh)30电池放电下限系数0.2电池充电上限系数0.8充放电效率0.95充放电功率(kW)5平均行驶速率(km/h)40行驶单位距离耗能(kWh)0.2

图1 电动汽车停驶概率 Fig.1 EV parking probability

表2 系统基础负荷 Tab.2 Basic loads of system

时段负荷时段负荷时段负荷11759325172502187.510350182753212.511362.5193004237.5123752035052501335021325627514325222757287.51530023225830016262.524200

4.2 不同模式的结果

根据第一节的方法对峰谷时段进行划分，其中峰时段为8:00～14:00、19：00～21:00，谷时段为22:00到次日5:00、15:00～18:00，其余为平时段。现在定义3种模式：模式1，同时考虑所有随机因素的波动；模式2，不考虑电价的波动性；模式3，不考虑碳排放权价格波动性。模式1和模式2电动汽车在调度周期的总出力(MWh)如表3所示，模式1和模式3常规机组和碳捕集机组在调度周期内的总出力(MWh)如表4所示。

表3 电动汽车出力 Tab.3 Output of electric vehicles

模式12电动汽车出力-81.3-71.8

表4 机组出力 Tab.4 Output of units

模式13G11 425.71 362.4G21 212.51 103.7G31 004.6943.7G4784.3707.1G5706.4674.3G6868.4789.3

表3中电动汽车总出力为负表示电动汽车充电电量有一部分用于自身消耗，一部分用于放电调度。从表中可以看出在考虑电价的波动性后，电动汽车的充电电量增加，电网用于支付电动汽车的成本降低，这表明为了应对电价的波动性，电力市场及时响应电价的变化，调整电动汽车的出力，这样使得目标函数趋于最优。

系统中的常规机组和碳捕集机组会向环境中排放二氧化碳，这些决策变量在目标函数中也与随机变量碳排放权价格直接相关。从表4可以看出考虑碳排放权价格波动性时，常规机组和碳捕集机组的出力都有所降低，这表明在碳交易市场背景下，为了响应碳排放权价格的波动性，调度决策会及时调整与其相关的变量来保证总的成本最低，在这里通过降低常规机组出力来降低碳交易成本。考虑电价和碳排放权波动性得到的结果与电力市场的实际情况更为接近，说明了模型的合理性。

4.3 置信水平的影响

置信水平决定了调度结果的可靠度，不同置信水平对结果的影响不同，图2给出了在不同置信水平下总成本值的变化情况，基础置信水平取0.90，其它置信水平分别取0.92、0.94、0.96、0.98、1。

从图2中可以看出，随着置信水平的增大，系统的总成本逐渐增加。置信水平决定了调度中随机因素的波动情况，置信水平越高，系统供电可靠性越高，能够在供需紧张时第一时间提供所需的电能，这时的电能会比平时具有更高的价值，电价会升高，也就是说在提高系统可靠性的同时需要付出经济性成本。

图2 置信水平对成本影响 Fig.2 Effect of confidence level on cost

4.4 电价波动水平影响

峰谷电价对电动汽车的行为有一定的引导作用，为了研究这一作用的程度，现定义3种情形：峰谷电价相对于平时段的电价浮动分别为20%、30%、40%，各个时段电价本身也具有波动性，如图3所示。

图3 电价对电动汽车出力影响 Fig.3 Effect of power price on EV outputs

电价浮动越大时，谷时段电价越低，峰时段电价越高。谷时段电价低，充电成本低，电动汽车用户越希望在谷时段尽可能让电动汽车充电；峰时段成本高，放电利润越大，这种趋势可以让电动汽车放电获取更多的收益。从图中看出，情形2和3相比基本情形来说，谷时段充电功率和峰时段的放电功率都增加，这说明电价的改变可以有效地刺激电汽车的充放电行为。

5 结束语

将模糊C均值聚类用于峰谷时段划分，结果更为合理；提出了推广的超分位数方法用于处理含多种随机变量的问题，将复杂的非线性规划问题进行简化，使得问题易于求解。分析了随机因素波动性、置信水平、电价浮动程度对于调度结果的影响，结果表明，考虑变量随机性的模型更符合实际情况，置信水平和电价对调度结果均有一定的影响，不仅说明了模型的合理性，也为电力系统的调度决策提供了参考。

参 考 文 献

[1] 侯建朝,胡群丰,谭忠富.计及需求响应的风电-电动汽车协同调度多目标优化模型[J].电力自动化设备,2016,(7): 22-27.

Hou Jianchao,Hun Qunfeng,Tan Zhongfu.Multi-objective optimization model of collaborative WP-EV dispatch considering demand response[J].Electric Power Automation Equipment,2016,36(7): 22-27.

[2] 程瑜,翟娜娜.基于用户响应的分时电价时段划分[J].电力系统自动化,2012,36(9): 42-46.

Cheng Yu,Di Nana.Electricity price peak and valley periods division based on customer response[J].Automation of Electric Power Systems,2012,36(9): 42-46.

[3] 赵洪山,王莹莹,陈松.需求响应对配电网供电可靠性的影响[J].电力系统自动化,2015,39(17): 49-55.

Zhao Hongshan,Wang Yingying,Chen Song.Impact of demand response on distribution system reliability[J].Automation of Electric Power Systems,2015,39(17): 49-55.

[4] 雷亚洲,王伟胜,印永华,等.含风电场电力系统的有功优化潮流[J].电网技术,2002,26(6): 18-21.

Lei Yazhou,Wang Weisheng,Yin Yonghua,et al.Optimal real power flow in wind power integrated system[J].Power System Technology,2002,26(6): 18-21.

[5] 周任军,姚龙华,童小娇,等.采用条件风险方法的含风电系统安全经济调度[J].中国电机工程学报,2011,31(1): 56-63.

Zhou Renjun,Yao Longhua,Tong Xiaojiao,et al.Security economic dispatch in wind power integrated system using conditional risk method[J].Proceedings of the CSEE,2011,31(1): 56-63.

[6] 董雷,杨以涵,张传成,等.综合考虑网络结构不确定性的概率潮流计算方法[J].电工技术学报,2012,27(1): 210-216.

Dong Lei,Yang Yihan,Zhang Chuancheng,et al.Probabilistic load flow considering network configuration uncertainties[J].Transactions of China Electrotechnical Society,2012,27(1): 210-216.

[7] 刘阳升,周任军,李星朗,等.碳排放权交易下碳捕集机组的厂内优化运行[J].电网技术,2013,37(2): 295-300.

Liu Yangsheng,Zhou Renjun,Li Xinglang,et al.Inside-plant optimal operation of carbon capture unit under carbon emission right trade[J].Power System Technology,2013,37(2): 295-300.

[8] 冉晓洪,苗世洪,刘阳升,等.考虑风光荷联合作用下的电力系统经济调度建模[J].中国电机工程学报,2014,34(16): 2552-2560.

Ran Xiaohong,Miao Shihong,Liu Yangsheng,et al.Modeling of economic dispatch of power system considering joint effect of wind power,solar energy and load[J].Proceedings of the CSEE,2014,34(16): 2552-2560.

[9] Rockafellar R T,Royset J O.On buffered failure probability in design and optimization of structures[J].Reliability Engineering and System Safety,2010,95(5): 499-510.

[10]HERTZER J,YU D C,BHATTARAI K.An economic dispatch model incorporating wind power[J].IEEE Trans on Energy Conversion,2008,23(2): 603-611.

[11]SALAMEH Z M,BOROWY B S,AMIN A R A.Photovoltaic module-site matching based on the capacity factors[J].IEEE Trans on Energy Conversion,1995,10(2): 326-332.

[12]IPCC.Special report on carbon dioxide capture and storage[D].Cambridge,UK: Cambridge University Press,2005.

[13]Benz E,Trück S.Modeling the price dynamics of C02 emission allowance[J].Energy Economics,2009,31(1): 4-15.

[14]刘西陲,沈炯,李益国.系统边际电价概率分布检验及模型研究[J].中国电机工程学报,2009,35(4): 72-77.

Liu Xichui,Shen Jiong,Li Yiguo.Study on probability distribution and model of system marginal price[J].Proceedings of the CSEE,2009,35(4): 72-77.

[15]KEMPTON W,LETENDRE S.Electric vehicles as a new power source for electric utilities[J].Transportation research: Part D,1997,2(3): 157-175.

[16]杨晓东,张晓兵,翁国庆,等.基于虚拟电价的电动汽车充放电优化调度及其实现机制研究[J].电工技术学报,2016,31(17): 52-62.

Yang Xiaodong,Zhang Xiaobing,Weng Guoqing,et al.Virtual Time-of-use tariffs based optimal scheduling and implementation mechanism of electric vehicles charging and discharging[J].Transactions of China Electrotechnical Society,2016,31(17): 52-62.

[17]张晓辉,闫柯柯,卢志刚,等.基于碳交易的含风电系统低碳经济调度[J].电网技术,2013,37(10): 2697-2704.

Zhang Xiaohui,Yan Keke,Lu Zhigang,et al.Carbon trading based low-carbon economic dispatching for power grid integrated with wind power system[J].Power System Technology,2013,37(10): 2697-2704.

[18]王贵斌,赵俊华,文福拴,等.配电系统中电动汽车与可再生能源的随机协同调度[J].电力系统自动化,2012,36(19): 22-29.

Wang Guibin,Zhao Junhua,Wen Fushuan,et al.Stochastic optimization dispatching of plug-in hybrid electric vehicles in coordination with renewable generation in distribution syetems[J].Automation of Electric Power Systems,2012,36(19): 22-29.