考虑潮汐资源评估与流速方向的潮汐流能发电场规划方法

摘要

本发明公开了考虑潮汐资源评估与流速方向的潮汐流能发电场规划方法，主要步骤为：1)获取潮汐流速、发电场成本、集电系统、粒子群算法的输入参数。2)基于所述输入参数，确定不同方向、不同幅值的流速出现的概率。3)计算规划海域内潮汐资源储量Pc，并根据潮汐资源储量Pc确定潮汐流能发电场的机组数量。4)建立潮汐流能发电场优化规划模型。5)基于潮汐流能发电场优化规划模型，规划出潮汐流能发电机组和海上变电站的初始位置。6)计算发电场日发电量，并更新潮汐流能发电机组和变电站坐标。7)对发电场内的机组进行分组，避免线路的跨区域交叉。8)完成集电系统拓扑规划。9)反复迭代，得到最优机组布局和集电系统规划方案。本发明可广泛应用于考虑潮汐资源评估和潮汐流速方向的潮汐流能发电场规划，能够为与潮汐流能发电系统相关的规划问题分析提供有益的参考。

权利要求

1.考虑潮汐资源评估与流速方向的潮汐流能发电场规划方法，其特征在于，主要包括以下步骤：

1)获取集电系统的输入参数；

2)基于所述输入参数，计算不同方向、不同幅值的流速出现的概率，主要步骤如下：

2.1)按照角度，将流速方向平均分为gd组；根据流速方向，对实测潮汐流速幅值vtm进行分组；设定每一组内的潮汐流速方向无差别；

实测潮汐流速幅值vtm落入第i组内的概率pi如下所示：式中，ndi为落入第i组内实测潮汐流速幅值数据的个数；i＝1,2,…,gd；nt为实测潮汐流速数据样本数量；

2.2)对同一方向内实测潮汐流速幅值进行聚类；

采用k均值聚类算法，将属于同一组内的潮汐流速幅值分为kd类；第j类典型流速为第j类潮汐流速幅值的平均值；

第i组内第j个典型流速的概率pij如下所示：式中，nij为落入第i组中第j类的数据个数；j＝1,2,…,kd；

其中，类别总数kd如下所示：

kd＝1+lg ndi；                  (3)

3)计算规划海域内潮汐资源储量Pc，并根据潮汐资源储量Pc确定潮汐流能发电场的机组数量NTCT；

规划海域内潮汐资源储量Pc如下所示：

Pc＝kexmaxymaxPdensity；                (4)式中，ke为环保系数；xmax与ymax分别为发电场规划海域的长和宽；Pdensity为海域潮汐能量密度；

其中，海域潮汐能量密度Pdensity如下所示：式中，vt为第t个时刻的潮汐流速幅值；t∈(0,T)；T为潮汐流速观测时长；Cp为潮汐流能发电机组的获能系数；ρ为海水密度；

潮汐流能发电场机组数量NTCT如下所示：

式中，NTCT为发电场内的机组总数，Prated为第n台机组的额定功率；

4)建立潮汐流能发电场优化规划模型，主要步骤如下：

4.1)确定潮汐流能发电场优化模型的目标函数min LCOE，即：式中，LCOE为发电场的平准化发电成本；Ctotal,TCT为机组的日等效投资成本；CECS为集电系统的日等效投资成本；Cdis为设备在退役后的拆解费用；Res为设备的残值；ET为发电场的日发电量；

机组的日等效投资成本Ctotal.TCT如下所示式中，NTCT为发电场内的机组总数；Cunit为机组单价，CTrans为机组配套变压器的单价，Cinstall为机组安装费用；r和L分别为设备的折现率和经济寿命；

集电系统的日等效投资成本CECS如下所示：

CECS＝C1+C2+C3+C4+C5；        (9)式中，C1为中压电缆日等效成本，C2为中压电缆的挖掘和维护日等效成本，C3为中压电缆日网损成本，C4为海上变电站的日等效成本，C5为高压电缆的日等效成本；

发电场的日发电量ET如下所示：

ET＝PT×24；        (10)

式中，PT为发电场的总出力期望值；

4.2)确定潮汐流能发电场优化模型的约束条件，所述约束条件主要包括第m台机组的坐标约束、机组和机组之间的距离约束、机组和变电站之间的距离约束、在第i个方向分组中的第j个典型流速下，第s条线路传输的功率约束和在第i个方向分组中的第j个典型流速下，第s条线路传输的压降约束；

5)基于潮汐流能发电场优化规划模型和粒子群算法，规划出潮汐流能发电机组和海上变电站的初始位置；

6)计算发电场日发电量ET，主要步骤如下：

6.1)根据公式10调整发电机组坐标，使流速方向沿x轴正半轴；

式中，(xm’,ym’)和(xm,ym)分别是调整后与调整前的第m台发电机组的坐标；θ为流速方向与x正半轴的夹角；

6.2)计算第m台机组产生的尾流与第n台机组叶片的重合面积Amn，即：式中，ymn为第m台和第n台机组纵坐标之差的绝对值，θw,mn和θ0,mn表示角度；r0为第n台机组的叶片半径；σmn为第m台机组产生的尾流区域在第n台机组处的半径；

半径σmn如下所示：

式中，I0代表湍流强度，xmn代表第m台机组和第n台机组横坐标之差的绝对值；

湍流强度函数c(I0)如下所示：

6.3)计算第n台机组处的流速vn，即：

式中，W(n)为第n台机组的上游机组集合，Nn为集合W(n)中的元素数量，A0为第n台机组叶片扫过的面积，CT为机组的推力系数；V0为自由流速；

6.4)计算第n台机组的输出功率Pn,ij，即：式中，vn,ij为在第i个流速方向分组中的第j个典型流速下，第n台机组处的流速；vin，vout和vrated分别为潮汐流能发电机组的切入流速、切出流速和额定流速；Cp为机组的获能系数；ρ为海水密度；Prated为第n台机组的额定功率；Pn,ij的概率为pij；

6.5)计算第n台机组出力的期望值PEn，即：式中，gd为流速方向分组数，kd为第i组内流速的类数；

潮汐流能发电场的出力期望值如下所示：

式中，NTCT为潮汐流能发电场中机组的总数；

6.6)根据输出功率期望值，计算出优化模型目标函数中的发电场日发电量ET；

7)基于发电场日发电量ET，更新潮汐流能发电机组和变电站坐标，主要步骤如下：

7.1)更新各个粒子的最佳位置和种群的最佳位置；

7.2)按照公式19和公式20更新粒子位置：式中， 和 分别为第d个粒子在第t代的移动速度和位置； 为第d个粒子到第t代为止的最佳位置；globalt为种群到第t代为止取得的最佳位置；w为惯性因子；l1和l2为学习因子；r1和r2为[0,1]之间的随机数；

式中， 和 为第d个粒子在第t+1代的移动速度和位置；

8)对潮汐流能发电机组进行分组，主要步骤如下：

8.1)确定分组数；

每一组内机组数上限nmax如下所示：

式中，Imv,max为中压海底电缆的最大载流量，Irated为每台机组的额定电流；

机组的分组数Ng如下所示：

式中，NTCT为发电场内的机组总数；

8.2)计算各台发电机组的夹角；

以变电站为原点建立笛卡尔坐标系；定义机组的夹角为机组与变电站的连线与x正半轴的夹角；第m台机组的夹角记为θm；θm∈[0,360°]，m＝1,2,…,NTCT；

8.3)对发电机组进行分组；

随机生成Ng条以原点为顶点的射线作为聚类中心，第g个聚类中心的角度为θg；θg∈[0,

360°]，g＝1,2,…,k；

分别计算各台机组与各聚类中心的夹角，并将机组归入夹角最小的组中；

判断是否有机组数越限情况产生；若第g组内的机组数ng>nmax，则在第g组中，选择与第g个聚类中心夹角最大的ng-nmax台机组，将其分入夹角第二小的组内，直到所有组均不越限为止；

对分入各个组内的机组与x正半轴的夹角值求平均，得到新的聚类中心；重复分组与调整过程，直至聚类中心不再变化；

9)完成集电系统拓扑规划，主要步骤如下：

9.1)针对布局Xd1，随机产生Nt个生成树连接各台机组，用以组成初始种群，并保证初始种群的线路不交叉；

9.2)对于各个已经连接完成的拓扑，采用Prufer数列对拓扑进行编码；

9.3)采用如下更新策略对拓扑粒子进行更新：式中， 和 为第h个粒子在第t代时的位置和速度；h＝1,2,…,Nt； 和 为第h个粒子在第t+1代时的位置和速度， 为第h个粒子到第t代为止的最佳位置，globalt为种群到第t代为止取得的最佳位置，w’为惯性因子，l′1、l'2为学习因子，r′1、r′2为[0,1]之间的随机数；

基于上述更新策略，每次迭代对种群中的粒子进行更新；

9.4)每次迭代更新粒子之后，对粒子进行Prufer解码，使其还原为拓扑结构；

9.5)根据解码完成后得到的拓扑结构，从拓扑末端的机组开始，沿线路向变电站侧推算，根据相应机组的出力期望值计算出每条线路中传输的功率，根据功率确定该条线路使用的电缆类型，进而计算集电系统的各项成本；

9.6)集电系统优化的迭代；

令集电系统规划的最大迭代次数为imax，当相邻两次迭代的集电系统总成本CECS差的绝对值小于临界阈值ε1，或迭代次数大于imax时迭代停止，该布局下的最佳集电系统拓扑及电缆选型即已得到；若不满足上述两个条件之一，则继续迭代；

将得到的最佳拓扑下的集电系统总成本CECS储存在公式7中的CECS中，用于下一次机组布局的迭代；

9.7)机组布局优化的迭代

机组布局种群中的所有粒子均得到最优集电系统之后，即完成了一次机组布局迭代；

令机组布局优化的最大迭代次数为Imax，则当相邻两次迭代的LCOE之差的绝对值小于邻接阈值ε2，或迭代次数大于Imax时迭代停止；若不满足上述两个条件之一，则继续迭代。

2.根据权利要求1所述的考虑潮汐资源评估与流速方向的潮汐流能发电场规划方法，其特征在于：所述输入参数主要包括实测潮汐流速幅值vtm与流速方向vtd、潮汐流速方向的分组数gd、潮汐流能发电机组的切入流速vin、额定流速vrated、切出流速vout、第n台机组的额定功率Prated、获能系数Cp、推力系数CT、叶片直径D、海水密度ρ、中压海底电缆的最大载流量Imv,max、海底电缆的单位电阻Rci、海底电缆的单位成本cmv,ci、潮汐流能发电场长方形规划海域的长xmax、宽ymax，机组最小安全距离Dmin、中压集电系统的额定电压U、海底电缆最大允许压降ΔUmax、粒子群算法种群规模Np、惯性因子w、学习因子l1、学习因子l2、最大迭代次数Imax、潮汐流能发电机组单价Cunit、发电机组配套小型变压器单价CTrans、机组安装费用Cinstall、海上变电站总价csub、高压海底电缆单价chv、设备折现率r、设备经济寿命L、电缆单位挖掘成本ct、单位网损成本cp和维修成本系数a；

其中，tm＝1,2,…,nt，td＝1,2,…,nt，nt为潮汐流速实测数据样本数；流速方向vtd与实测潮汐流速幅值vtm一一对应；ci＝1,2,…,nc，nc为海底电缆的类别数量。

3.根据权利要求1或2所述的考虑潮汐资源评估与流速方向的潮汐流能发电场规划方法，其特征在于，中压电缆日等效成本C1如下所示：式中，nl代表中压电缆条数；kv代表典型流速的数量；cf,mv为第c类中压电缆的单价；f＝

1,2,…,nf,nf为电缆种类数；ls为第s条中压电缆的长度；s＝1,2,…,nl；

中压电缆的挖掘和维护日等效成本C2如下所示：式中，ct为单位长度电缆的运输与铺设成本；a代表电缆维护系数；

中压电缆日网损成本C3如下所示：

式中，cp为单位电量损耗成本；PEn为第n台机组的出力期望值；Np,s为功率流过第s条电缆的机组数；U为中压集电系统额定电压；rf为第f类中压电缆的单位长度电阻；

海上变电站的日等效成本C4如下所示：

式中，csub为海上升压变电站总价；

高压电缆的日等效成本C5如下所示：

式中，chv为高压海底电缆的单价；lhv为高压电缆的长度；

中压电缆条数nl如下所示：

nl＝NTCT；       (31) 。

4.根据权利要求1所述的考虑潮汐资源评估与流速方向的潮汐流能发电场规划方法，其特征在于，第m台机组的横坐标xm约束如下所示：

0＜xm＜xmax；      (32)

第m台机组的纵坐标ym约束如下所示：

0＜ym＜ymax；        (33)机组和机组之间的距离约束如下所示：

式中，(xm,ym)为第m台机组的坐标；(xn,yn)为第n台机组的坐标；m,n＝1,2,…,NTCT，且m≠n；Dmin1为机组与机组之间的最小安全距离；

机组和变电站之间的距离约束如下所示：

式中，Dmin2为机组与变电站之间的最小安全距离； 为变电站的坐标；

在第i个方向分组中的第j个典型流速下，第s条线路传输的功率Ps,ij约束如下所示：Ps,ij≤Ps,max；         (36)式中，Ps,max为第s条线路的最大允许传输功率；

在第i个方向分组中的第j个典型流速下，第s条线路传输的压降ΔUs,ij约束如下所示：ΔUs,ij≤ΔUmax；        (37)式中，ΔUmax为最大允许压降。

说明书

考虑潮汐资源评估与流速方向的潮汐流能发电场规划方法

技术领域

[0001] 本发明涉及电力系统规划领域，具体是考虑潮汐资源评估与流速方向的潮汐流能发电场规划方法。

背景技术

[0002] 面积广阔的海洋可以为人类提供多种可再生能源，潮汐流能作为其中的一种，由于其具有环境友好性、功率密度大、规律性强等优点，近年来发展势头良好且关注度逐渐提高，已有多个国际组织及国际公司开展潮汐流能发电技术的研究，并已有潮汐流能发电场投产运行或处于建设之中。

[0003] 潮汐流能发电场主要由发电机组和集电系统构成；发电机组目前主要有两种，分别是利用潮汐势能发电的坝式机组，以及利用潮汐动能发电的潮汐流能发电机组；由于潮汐流能发电机组建设工程量较小，并且不阻碍航道，因此该种机组更加适用；集电系统主要由海底电缆、海上变电站、开关等设备构成，其作用是将机组发出的电能收集并输送到陆上电网。

[0004] 发电场的收入依靠售电，为了提高发电场的发电量，需要着重考虑发电场规划区域内的潮汐流能资源、尾流效应以及潮汐流速的特征，具体为：

[0005] 1)规划海域内，潮汐流能的储量是有限的，结合资源储量确定发电场的容量是充分挖掘规划区域和机组发电潜能的有效手段；

[0006] 2)尾流效应是指机组从潮汐流中获取能量的同时，在其下游形成潮汐流速下降的尾流区；增大机组间距能够削弱尾流效应对机组发电量的影响，但增大间距的同时，也会增加海底电缆的长度和成本，以及运行成本；与海上风电相比，潮汐具有极强的方向性，基于潮汐流速的方向性合理规划发电场，是提高发电场产出，并降低投资与运行成本的另一种有效手段。

[0007] 在潮汐流能发电场规划中，若不计及潮汐资源储量与方向性，将严重影响规划方法的工程实用价值，削弱发电场的盈利能力；然而，在目前的研究中，尚未见同时考虑潮汐资源储量与潮汐流速方向性的规划方法，尚属空白。

发明内容

[0008] 本发明的目的是解决现有技术中存在的问题。

[0009] 为实现本发明目的而采用的技术方案是这样的，考虑潮汐资源评估与流速方向的潮汐流能发电场规划方法，主要包括以下步骤：

[0010] 1)获取集电系统的输入参数。

[0011] 所述输入参数主要包括实测潮汐流速幅值vtm与流速方向vtd、潮汐流速方向的分组数gd、潮汐流能发电机组的切入流速vin、额定流速vrated、切出流速vout、额定输出功率Prated、获能系数Cp、推力系数CT、叶片直径D、海水密度ρ、中压海底电缆的最大载流量Capci、海底电缆的单位电阻Rci、海底电缆的单位成本cmv,ci、潮汐流能发电场长方形规划海域的长xmax、宽ymax，机组最小安全距离Dmin、中压集电系统的额定电压U、海底电缆最大允许压降ΔUmax、粒子群算法种群规模Np、惯性权重w、学习因子l1、学习因子l2、最大迭代次数Imax、潮汐流能发电机组单价Cunit、发电机组配套小型变压器单价CTrans、机组安装费用Cinstall、海上变电站总价csub、高压海底电缆单价chv、设备折现率r、设备经济寿命L、电缆单位挖掘成本ct、单位网损成本cp和维修成本系数a。

[0012] 其中，tm＝1,2,…,nt，td＝1,2,…,nt，nt为潮汐流速实测数据样本数。流速方向vtd与实测潮汐流速幅值vtm一一对应。ci＝1,2,…,nc，nc为海底电缆的类别数量。

[0013] 2)基于所述输入参数，计算不同方向、不同幅值的流速出现的概率，主要步骤如下：

[0014] 2.1)按照角度，将流速方向平均分为gd组。根据流速方向，对实测潮汐流速幅值vtm进行分组。设定每一组内的潮汐流速数据方向无差别。

[0015] 实测潮汐流速幅值vtm落入第i组内的概率pi如下所示：

[0016]

[0017] 式中，ndi为落入第i组内实测潮汐流速幅值数据的个数。i＝1,2,…,gd。nt为实测潮汐流速数据样本数量。

[0018] 2.2)对同一方向内实测潮汐流速幅值进行聚类。

[0019] 采用k均值聚类算法，将属于同一组内的潮汐流速幅值分为kd类。第j类典型流速为第j类潮汐流速幅值的平均值。

[0020] 第i组内第j个典型流速的概率pij如下所示：

[0021]

[0022] 式中，nij为落入第i组中第j类的数据个数。j＝1,2,…,kd。

[0023] 其中，类别总数kd如下所示：

[0024] kd＝1+lgndi；   (3)

[0025] 3)计算规划海域内潮汐资源储量Pc，并根据潮汐资源储量Pc确定潮汐流能发电场的机组数量NTCT。

[0026] 海域潮汐能量密度Pdensity如下所示：

[0027]

[0028] 式中，vt为第t个时刻的潮汐流速幅值。t∈(0,T)。T为潮汐流速观测时长。

[0029] 规划海域内潮汐资源储量Pc如下所示：

[0030] Pc＝kexmaxymaxPdensity。   (5)

[0031] 式中，ke为环保系数。xmax与ymax分别为发电场规划海域的长和宽。Pdensity为海域潮汐能量密度。

[0032] 潮汐流能发电场机组数量NTCT如下所示：

[0033]

[0034] 式中，NTCT为发电场内的机组总数，Prated为单台机组的额定容量。

[0035] 4)建立潮汐流能发电场优化规划模型，主要步骤如下：

[0036] 4.1)确定潮汐流能发电场优化模型的目标函数min LCOE，即：

[0037]

[0038] 式中，LCOE为发电场的平准化发电成本。Ctotal.TCT为机组的日等效投资成本。CECS为集电系统的日等效投资成本。Cdis为设备在退役后的拆解费用。Res为设备的残值。ET为发电场的日发电量。

[0039] 机组的日等效投资成本Ctotal.TCT如下所示

[0040]

[0041] 式中，NTCT为发电场内的机组总数。Cunit为机组单价，CTrans为机组配套变压器的单价，Cinstall为机组安装费用。r和L分别为设备的折现率和经济寿命。

[0042] 集电系统的日等效投资成本CECS如下所示：

[0043] CECS＝C1+C2+C3+C4+C5。   (9)

[0044] 式中，C1为中压电缆日等效成本，C2为中压电缆的挖掘和维护日等效成本，C3为中压电缆日网损成本，C4为海上变电站的日等效成本，C5为高压电缆的日等效成本。

[0045] 中压电缆日等效成本C1如下所示：

[0046]

[0047] 式中，nl代表中压电缆条数。kv代表典型流速的数量。cf.mv为第c类中压电缆的单价。f＝1,2,…,nf,nf为电缆种类数。ls为第s条中压电缆的长度。s＝1,2,…,nl。

[0048] 中压电缆的挖掘和维护日等效成本C2如下所示：

[0049]

[0050] 式中，ct为单位长度电缆的运输与铺设成本。a代表电缆维护系数。

[0051] 中压电缆日网损成本C3如下所示：

[0052]

[0053] 式中，cp为单位电量损耗成本。PEn为第n台机组的出力期望值。Np,s为功率流过第s条电缆的机组数。U为中压集电系统额定电压。rf为第f类中压电缆的单位长度电阻。

[0054] 海上变电站的日等效成本C4如下所示：

[0055]

[0056] 式中，csub为海上升压变电站总价。

[0057] 高压电缆的日等效成本C5如下所示：

[0058]

[0059] 式中，chv为高压海底电缆的单价。lhv为高压电缆的长度。

[0060] 发电场的日发电量ET如下所示：

[0061] ET＝PT×24。   (15)

[0062] 式中，PT为发电场的总出力期望值。

[0063] 中压电缆条数nl如下所示：

[0064] nl＝NTCT。   (16)

[0065] 4.2)确定潮汐流能发电场优化模型的约束条件，所述约束条件主要包括第m台机组的坐标约束、机组和机组之间的距离约束、机组和变电站之间的距离约束、在第i个方向分组中的第j个典型流速下，第s条线路传输的功率约束和在第i个方向分组中的第j个典型流速下，第s条线路传输的压降约束。

[0066] 第m台机组的横坐标xm约束如下所示：

[0067] 0＜xm＜xmax。   (17)

[0068] 第m台机组的纵坐标ym约束如下所示：

[0069] 0＜ym＜ymax。   (18)

[0070] 机组和机组之间的距离约束如下所示：

[0071]

[0072] 式中，(xm,ym)为第m台机组的坐标。(xn,yn)为第n台机组的坐标。m,n＝1,2,…,NTCT，且m≠n。Dmin1为机组与机组之间的最小安全距离。

[0073] 机组和变电站之间的距离约束如下所示：

[0074]

[0075] 式中，Dmin2为机组与变电站之间的最小安全距离。 为变电站的坐标。

[0076] 在第i个方向分组中的第j个典型流速下，第s条线路传输的功率Ps,ij约束如下所示：

[0077] Ps,ij≤Ps,max。   (21)

[0078] 式中，Ps,max为第s条线路的最大允许传输功率。

[0079] 在第i个方向分组中的第j个典型流速下，第s条线路传输的压降ΔUs.ij约束如下所示：

[0080] ΔUs,ij≤ΔUmax。   (22)

[0081] 式中，ΔUmax为最大允许压降。

[0082] 5)基于潮汐流能发电场优化规划模型和粒子群算法，规划出潮汐流能发电机组和海上变电站的初始位置。

[0083] 生成潮汐流能发电机组与海上变电站初始位置的主要步骤如下：

[0084] 5.1)采用粒子群算法的粒子表示潮汐流能发电场中每台机组和海上变电站的坐标。

[0085] 粒子种群中共有Np个粒子，其中第d个粒子Xd如下所示：

[0086]

[0087] 式中， 分别代表第1台、第2台、…、第NTCT台机组的横纵坐标。 代表海上变电站的横纵坐标。

[0088] 5.2)基于公式23，随机生成Np个粒子的坐标，从而构成初始粒子种群。

[0089] 6)计算发电场日发电量ET，主要步骤如下：

[0090] 6.1)根据公式24调整发电机组坐标，使流速方向沿x轴正半轴。

[0091]

[0092] 式中，(xm’,ym’)和(xm,ym)分别是调整后与调整前的第m台发电机组的坐标。θ为流速方向与x正半轴的夹角。

[0093] 6.2)计算第m台机组产生的尾流与第n台机组叶片的重合面积Amn，即：

[0094]

[0095]

[0096] 式中，ymn为第m台和第n台机组纵坐标之差的绝对值，θw,mn和θ0,mn表示角度。r0为第n台机组的叶片半径。σmn为第m台机组产生的尾流区域在第n台机组处的半径。

[0097] 半径σmn如下所示：

[0098]

[0099] 式中，I0代表湍流强度，xmn代表第m台机组和第n台机组横坐标之差的绝对值。

[0100] 湍流强度函数c(I0)可以由下式计算：

[0101]

[0102] 6.3)计算第n台机组处的流速vn，即：

[0103]

[0104] 式中，W(n)为第n台机组的上游机组集合，Nn为集合W(n)中的元素数量，A0为第n台机组叶片扫过的面积，CT为机组的推力系数。V0为自由流速。

[0105] 6.4)计算第n台机组的输出功率Pn,ij，即：

[0106]

[0107] 式中，vn,ij为在第i个流速方向分组中的第j个典型流速下，第n台机组处的流速。vin，vout和vrated分别为潮汐流能发电机组的切入流速、切出流速和额定流速。Cp为机组的获能系数。ρ为海水密度。Prated为第n台机组的额定功率。Pn,ij的概率为pij。

[0108] 6.5)计算第n台机组出力的期望值PEn，即：

[0109]

[0110] 式中，gd为流速方向分组数，kd为第i组内流速的类数。

[0111] 潮汐流能发电场的出力期望值如下所示：

[0112]

[0113] 式中，NTCT为潮汐流能发电场中机组的总数。

[0114] 6.6)根据输出功率期望值，计算出优化模型目标函数中的发电场日发电量ET。

[0115] 7)基于发电场日发电量ET，更新潮汐流能发电机组和变电站坐标，主要步骤如下：

[0116] 7.1)更新各个粒子的最佳位置和种群的最佳位置。

[0117] 7.2)按照公式33和公式34更新粒子位置：

[0118]

[0119] 式中， 和 分别为第d个粒子在第t代的移动速度和位置。 为第d个粒子到第t代为止的最佳位置。globalt为种群到第t代为止取得的最佳位置。w为惯性因子。l1和l2为学习因子。r1和r2为[0,1]之间的随机数。

[0120]

[0121] 式中， 和 为第d个粒子在第t+1代的移动速度和位置。

[0122] 8)对潮汐流能发电机组进行分组，主要步骤如下：

[0123] 8.1)确定分组数。

[0124] 每一组内机组数上限nmax如下所示：

[0125]

[0126] 式中，Imv,max为中压电缆的最大载流量，Irated为每台机组的额定电流。

[0127] 机组的分组数Ng如下所示：

[0128]

[0129] 式中，NTCT为发电场内的机组总数。

[0130] 8.2)计算各台发电机组的夹角。

[0131] 以变电站为原点建立笛卡尔坐标系。定义机组的夹角为机组与变电站的连线与x正半轴的夹角。第m台机组的夹角记为θm。θm∈[0,360°]，m＝1,2,…,NTCT。

[0132] 8.3)对发电机组进行分组。

[0133] 随机生成Ng条以原点为顶点的射线作为聚类中心，第g个聚类中心的角度为θg。θg∈[0,360°]，g＝1,2,…,k。

[0134] 分别计算各台机组与各聚类中心的夹角，并将机组归入夹角最小的组中。

[0135] 判断是否有机组数越限情况产生。若第g组内的机组数ng>nmax，则在第g组中，选择与第g个聚类中心夹角最大的ng-nmax台机组，将其分入夹角第二小的组内，直到所有组均不越限为止。

[0136] 对分入各个组内的机组与x正半轴的夹角值求平均，得到新的聚类中心。重复分组与调整过程，直至聚类中心不再变化。

[0137] 9)完成集电系统拓扑规划，主要步骤如下：

[0138] 9.1)针对布局Xd1，随机产生Nt个生成树连接各台机组，用以组成初始种群，并保证初始种群的线路不交叉。

[0139] 9.2)对于各个已经连接完成的拓扑，采用Prufer数列对拓扑进行编码。

[0140] 9.3)采用如下更新策略对拓扑粒子进行更新：

[0141]

[0142]

[0143] 式中， 和 为第h个粒子在第t代时的位置和速度。h＝1,2,…,Nt。 和为第h个粒子在第t+1代时的位置和速度， 为第h个粒子到第t代为止的最佳位置，globalt为种群到第t代为止取得的最佳位置，w’为惯性因子，l′1、l′2为学习因子，r′1、r′2为[0,1]之间的随机数。

[0144] 基于上述更新策略，每次迭代对种群中的粒子进行更新。

[0145] 9.4)每次迭代更新粒子之后，对粒子进行Prufer解码，使其还原为拓扑结构。

[0146] 9.5)根据解码完成后得到的拓扑结构，从拓扑末端的机组开始，沿线路向变电站侧推算，根据相应机组的出力期望值计算出每条线路中传输的功率，根据功率确定该条线路使用的电缆类型，进而计算集电系统的各项成本。

[0147] 9.6)集电系统优化的迭代。

[0148] 令集电系统规划的最大迭代次数为imax，当相邻两次迭代的集电系统总成本CECS差的绝对值小于临界阈值ε1，或迭代次数大于imax时迭代停止，该布局下的最佳集电系统拓扑及电缆选型即已得到。若不满足上述两个条件之一，则继续迭代。

[0149] 将得到的最佳拓扑下的集电系统总成本CECS储存在公式(7)中的CECS中，用于下一次机组布局的迭代。

[0150] 9.7)机组布局优化的迭代

[0151] 机组布局种群中的所有粒子均得到最优集电系统之后，即完成了一次机组布局迭代。令机组布局优化的最大迭代次数为Imax，则当相邻两次迭代的LCOE之差的绝对值小于邻接阈值ε2，或迭代次数大于Imax时迭代停止。若不满足上述两个条件之一，则继续迭代。

[0152] 本发明的技术效果是毋庸置疑的。本发明方法充分计及规划海域潮汐资源储量对发电场规划的影响，该方法能够充分挖掘规划海域的潮汐发电潜能以及潮汐流能发电机组的利用率，从而有效提高发电场的发电量，降低发电成本。

[0153] 本发明在规划机组布局和集电系统时，充分考虑了潮汐流速的方向，不仅能够提高尾流和机组出力的计算精度，还能使机组按照潮汐流速方向布局，进一步提高发电场的发电量。

[0154] 本发明可广泛应用于考虑潮汐资源评估和潮汐流速方向的潮汐流能发电场规划，能够为与潮汐流能发电系统相关的规划问题分析提供有益的参考。

附图说明

[0155] 图1为基于实测流速数据做出的潮汐玫瑰图；

[0156] 图2为尾流区域与机组叶片重合面积示意图；

[0157] 图3为算例中用于说明编码与解码方法的拓扑示意图；

[0158] 图4(a)为机组布局方案；

[0159] 图4(b)为集电系统的拓扑；

[0160] 图5为方法流程图。

具体实施方式

[0161] 下面结合实施例对本发明作进一步说明，但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下，根据本领域普通技术知识和惯用手段，做出各种替换和变更，均应包括在本发明的保护范围内。

[0162] 实施例1：

[0163] 参见图1至图5，考虑潮汐资源评估与流速方向的潮汐流能发电场规划方法，主要包括以下步骤：

[0164] 1)获取集电系统，也即电力网络的输入参数。

[0165] 所述输入参数主要包括实测潮汐流速幅值vtm与流速方向vtd、潮汐流速方向的分组数gd、潮汐流能发电机组的切入流速vin、额定流速vrated、切出流速vout、额定输出功率Prated、获能系数Cp、推力系数CT、叶片直径D、海水密度ρ、中压海底电缆的最大载流量Capci、海底电缆的单位电阻Rci、海底电缆的单位成本cmv,ci、潮汐流能发电场长方形规划海域的长xmax、宽ymax，机组最小安全距离Dmin、中压集电系统的额定电压U、海底电缆最大允许压降ΔUmax、粒子群算法种群规模Np、惯性权重w、学习因子l1、学习因子l2、最大迭代次数Imax、潮汐流能发电机组单价Cunit、发电机组配套小型变压器单价CTrans、机组安装费用Cinstall、海上变电站总价csub、高压海底电缆单价chv、设备折现率r、设备经济寿命L、电缆单位挖掘成本ct、单位网损成本cp和维修成本系数a。

[0166] 其中，tm＝1,2,…,nt，td＝1,2,…,nt，nt为潮汐流速实测数据样本数。流速方向vtd与实测潮汐流速幅值vtm一一对应。ci＝1,2,…,nc，nc为海底电缆的类别数量。

[0167] 2)基于所述输入参数，评估规划海域流速特性，从而计算不同方向、不同幅值的流速出现的概率，主要步骤如下：

[0168] 2.1)将[0,360°]的方向范围平均分为gd组。根据流速方向，对实测潮汐流速幅值vtm进行分组。设定每一组内的潮汐流速数据方向无差别。

[0169] 实测潮汐流速幅值vtm落入第i组内的概率pi如下所示：

[0170]

[0171] 式中，ndi为落入第i组内实测潮汐流速幅值数据的个数。i＝1,2,…,gd。nt为实测潮汐流速数据样本数量。

[0172] 2.2)对同一方向内实测潮汐流速幅值进行聚类。

[0173] 采用k均值聚类算法，将属于同一组内的潮汐流速幅值分为kd类。第j类典型流速为第j类潮汐流速幅值的平均值。

[0174] 第i组内第j个典型流速的概率pij如下所示：

[0175]

[0176] 式中，nij为落入第i组中第j类的数据个数。j＝1,2,…,kd。

[0177] 其中，类别总数kd如下所示：

[0178] kd＝1+lgndi；   (3)

[0179] 2.3)基于概率pi、概率pij和第j类典型流速，建立潮汐玫瑰图。所述潮汐玫瑰图表示不同方向、不同幅值的流速出现的概率。

[0180] 3)计算规划海域内潮汐资源储量Pc，并根据潮汐资源储量Pc确定潮汐流能发电场的机组数量NTCT。

[0181] 海域潮汐能量密度Pdensity如下所示：

[0182]

[0183] 式中，vt为第t个时刻的潮汐流速幅值。t∈(0,T)。T为潮汐流速观测时长，本实施例设定为24。

[0184] 规划海域内潮汐资源储量Pc如下所示：

[0185] Pc＝kexmaxymaxPdensity。   (5)

[0186] 式中，ke为环保系数。xmax与ymax分别为发电场规划海域的长和宽。Pdensity为海域潮汐能量密度。

[0187] 潮汐流能发电场机组数量NTCT如下所示：

[0188]

[0189] 式中，NTCT为发电场内的机组总数，Prated为单台机组的额定容量。

[0190] 4)建立潮汐流能发电场优化规划模型，主要步骤如下：

[0191] 4.1)确定潮汐流能发电场优化模型的目标函数min LCOE，即：

[0192]

[0193] 式中，LCOE为发电场的平准化发电成本。Ctotal.TCT为机组的日等效投资成本。CECS为集电系统的日等效投资成本。Cdis为设备在退役后的拆解费用。Res为设备的残值。ET为发电场的日发电量。计算时认为Cdis＝Res。

[0194] 机组的日等效投资成本Ctotal.TCT如下所示

[0195]

[0196] 式中，NTCT为发电场内的机组总数。Cunit为机组单价，CTrans为机组配套变压器的单价，Cinstall为机组安装费用。r和L分别为设备的折现率和经济寿命。

[0197] 集电系统的日等效投资成本CECS如下所示：

[0198] CECS＝C1+C2+C3+C4+C5。   (9)

[0199] 式中，C1为中压电缆日等效成本，C2为中压电缆的挖掘和维护日等效成本，C3为中压电缆日网损成本，C4为海上变电站的日等效成本，C5为高压电缆的日等效成本。

[0200] 中压电缆日等效成本C1如下所示：

[0201]

[0202] 式中，nl代表中压电缆条数。kv代表典型流速的数量。cf.mv为第c类中压电缆的单价。f＝1,2,…,nf,nf为电缆种类数。ls为第s条中压电缆的长度。s＝1,2,…,nl。

[0203] 中压电缆的挖掘和维护日等效成本C2如下所示：

[0204]

[0205] 式中，ct为单位长度电缆的运输与铺设成本。a代表电缆维护系数。

[0206] 中压电缆日网损成本C3如下所示：

[0207]

[0208] 式中，cp为单位电量损耗成本。PEn为第n台机组的出力期望值。Np,s为功率流过第s条电缆的机组数。U为中压集电系统额定电压。rf为第f类中压电缆的单位长度电阻。

[0209] 海上变电站的日等效成本C4如下所示：

[0210]

[0211] 式中，csub为海上升压变电站总价。

[0212] 高压电缆的日等效成本C5如下所示：

[0213]

[0214] 式中，chv为高压海底电缆的单价。lhv为高压电缆的长度。

[0215] 发电场的日发电量ET如下所示：

[0216] ET＝PT×24。   (15)

[0217] 式中，PT为发电场的总出力期望值。

[0218] 中压电缆条数nl如下所示：

[0219] nl＝NTCT。   (16)

[0220] 4.2)确定潮汐流能发电场优化模型的约束条件，所述约束条件主要包括第m台机组的坐标约束、机组和机组之间的距离约束、机组和变电站之间的距离约束、在第i个方向分组中的第j个典型流速下，第s条线路传输的功率约束和在第i个方向分组中的第j个典型流速下，第s条线路传输的压降约束。

[0221] 第m台机组的横坐标xm约束如下所示：

[0222] 0＜xm＜xmax。   (17)

[0223] 第m台机组的纵坐标ym约束如下所示：

[0224] 0＜ym＜ymax。   (18)

[0225] 机组和机组之间的距离约束如下所示：

[0226]

[0227] 式中，(xm,ym)为第m台机组的坐标。(xn,yn)为第n台机组的坐标。m,n＝1,2,…,NTCT，且m≠n。Dmin1为机组与机组之间的最小安全距离。Dmin1＝Dmin2＝54m。

[0228] 机组和变电站之间的距离约束如下所示：

[0229]

[0230] 式中，Dmin2为机组与变电站之间的最小安全距离。 为变电站的坐标。

[0231] 在第i个方向分组中的第j个典型流速下，第s条线路传输的功率Ps,ij约束如下所示：

[0232] Ps,ij≤Ps,max。   (21)

[0233] 式中，Ps,max为第s条线路的最大允许传输功率。Ps,ij为在第i个方向分组中的第j个典型流速下，第s条线路传输的功率。

[0234] 在第i个方向分组中的第j个典型流速下，第s条线路传输的压降ΔUs.ij约束如下所示：

[0235] ΔUs,ij≤ΔUmax。   (22)

[0236] 式中，ΔUmax为最大允许压降。ΔUmax＝0.33kV。ΔUs.ij为在第i个方向分组中的第j个典型流速下，第s条线路传输的压降。

[0237] 5)基于潮汐流能发电场优化规划模型和粒子群算法，规划出潮汐流能发电机组和海上变电站的初始位置。

[0238] 生成潮汐流能发电机组与海上变电站初始位置的主要步骤如下：

[0239] 5.1)采用粒子群算法的粒子表示潮汐流能发电场中每台机组和海上变电站的坐标。

[0240] 粒子种群中共有Np个粒子，其中第d个粒子Xd如下所示：

[0241]

[0242] 式中， 分别代表第1台、第2台、…、第NTCT台机组的横纵坐标。 代表海上变电站的横纵坐标。

[0243] 5.2)基于公式23，随机生成Np个粒子的坐标，从而构成初始粒子种群。

[0244] 6)计算发电场日发电量ET，主要步骤如下：

[0245] 6.1)根据公式24调整发电机组坐标，使流速方向沿x轴正半轴。

[0246]

[0247] 式中，(xm’,ym’)和(xm,ym)分别是调整后与调整前的第m台发电机组的坐标。θ为流速方向与x正半轴的夹角。

[0248] 6.2)计算第m台机组产生的尾流与第n台机组叶片的重合面积Amn，即：

[0249]

[0250]

[0251] 式中，ymn为第m台和第n台机组纵坐标之差的绝对值，θw,mn和θ0,mn表示角度。r0为第n台机组的叶片半径。σmn为第m台机组产生的尾流区域在第n台机组处的半径。

[0252] 半径σmn如下所示：

[0253]

[0254] 式中，I0代表湍流强度，xmn代表第m台机组和第n台机组横坐标之差的绝对值。

[0255] 湍流强度函数c(I0)可以由下式计算：

[0256]

[0257] 6.3)计算第n台机组处的流速vn，即：

[0258]

[0259] 式中，W(n)为第n台机组的上游机组集合，Nn为集合W(n)中的元素数量，A0为第n台机组叶片扫过的面积，CT为机组的推力系数。V0为自由流速。

[0260] 6.4)计算第n台机组的输出功率Pn,ij，即：

[0261]

[0262] 式中，vn,ij为在第i个流速方向分组中的第j个典型流速下，第n台机组处的流速。vin，vout和vrated分别为潮汐流能发电机组的切入流速、切出流速和额定流速。Cp为机组的获能系数。ρ为海水密度。Prated为第n台机组的额定功率。Pn,ij的概率为pij。

[0263] 6.5)计算第n台机组出力的期望值PEn，即：

[0264]

[0265] 式中，gd为流速方向分组数，kd为第i组内流速的类数。

[0266] 潮汐流能发电场的出力期望值如下所示：

[0267]

[0268] 式中，NTCT为潮汐流能发电场中机组的总数。

[0269] 6.6)根据输出功率期望值，计算出优化模型目标函数中的发电场日发电量ET。

[0270] 7)基于发电场日发电量ET，更新潮汐流能发电机组和变电站坐标，主要步骤如下：

[0271] 7.1)更新各个粒子的最佳位置和种群的最佳位置。

[0272] 7.2)按照公式33和公式34更新粒子位置：

[0273]

[0274] 式中， 和 分别为第d个粒子在第t代的移动速度和位置。 为第d个粒子到第t代为止的最佳位置。globalt为种群到第t代为止取得的最佳位置。w为惯性因子。l1和l2为学习因子。r1和r2为[0,1]之间的随机数。

[0275]

[0276] 式中， 和 为第d个粒子在第t+1代的移动速度和位置。

[0277] 8)对潮汐流能发电机组进行分组，主要步骤如下：

[0278] 8.1)确定分组数。

[0279] 每一组内机组数上限nmax如下所示：

[0280]

[0281] 式中，Imv,max为中压电缆的最大载流量，Irated为每台机组的额定电流。

[0282] 机组的分组数Ng如下所示：

[0283]

[0284] 式中，NTCT为发电场内的机组总数。

[0285] 8.2)计算各台发电机组的夹角。

[0286] 以变电站为原点建立笛卡尔坐标系。定义机组的夹角为机组与变电站的连线与x正半轴的夹角。第m台机组的夹角记为θm。θm∈[0,360°]，m＝1,2,…,NTCT。

[0287] 8.3)对发电机组进行分组。

[0288] 随机生成Ng条以原点为顶点的射线作为聚类中心，第g个聚类中心的角度为θg。θg∈[0,360°]，g＝1,2,…,k。

[0289] 分别计算各台机组与各聚类中心的夹角，并将机组归入夹角最小的组中。

[0290] 判断是否有机组数越限情况产生。若第g组内的机组数ng>nmax，则在第g组中，选择与第g个聚类中心夹角最大的ng-nmax台机组，将其分入夹角第二小的组内，直到所有组均不越限为止。

[0291] 对分入各个组内的机组与x正半轴的夹角值求平均，得到新的聚类中心。重复分组与调整过程，直至聚类中心不再变化。

[0292] 9)完成集电系统拓扑规划，主要步骤如下：

[0293] 9.1)针对布局Xd1，随机产生Nt个生成树连接各台机组，用以组成初始种群，并保证初始种群的线路不交叉。

[0294] 9.2)对于各个已经连接完成的拓扑，采用Prufer数列对拓扑进行编码。

[0295] Prufer编码的步骤如下：

[0296] a)定义三个向量U、P和Q，其中，U包含本组内的所有节点。节点数量为nf，节点编号为1～nf。P用来储存Prufer编码，Q代表没有包含在P中的节点。

[0297] b)从U中找到只与一个节点相连的节点，并选择其中编号最小的节点b定义为叶子节点。

[0298] c)将节点b所连接的节点d放入向量P中，并且每次均将新加入的节点放在最右边。

[0299] d)从Q中删去叶子节点b并从拓扑中删去相应的线路(b,d)。

[0300] e)重复步骤b)-d)，直到在拓扑中仅剩一条线路为止；此时，向量P由nf-2个元素组成，将U内未在向量P中出现的所有元素放入Q中。

[0301] 9.3)采用如下更新策略对拓扑粒子进行更新：

[0302]

[0303]

[0304] 式中， 和 为第h个粒子在第t代时的位置和速度。h＝1,2,…,Nt。 和 为第h个粒子在第t+1代时的位置和速度， 为第h个粒子到第t代为止的最佳位置，globalt为种群到第t代为止取得的最佳位置，w’为惯性因子，l′1、l′2为学习因子，r′1、r′2为[0,1]之间的随机数。w′为惯性因子。

[0305] 更新策略中的概念与运算符定义如下：

[0306] a)交换序：粒子速度由一系列交换序组成，交换序(ik,jk)表示交换粒子中值为ik与jk分量的位置。

[0307] b)加法操作 粒子位置与速度之间的加法操作，其结果为粒子位置；假设粒子k在第t代的位置为 则 表示交换 中相应分量的位置，即得到第t+1代粒子的位置 粒子速度与速度之间的加法操作，其结果为粒子速度，将第二个速度的交换序连接到第一个速度交换序的末尾。

[0308] c)减法操作Θ：粒子位置与粒子位置之间的减法操作，其结果为粒子速度，即一组交换序列。

[0309] d)乘法操作 对实数c∈(0,1)，假设速度vt为一组共k个交换序，乘法操作表示截取速度列表，使得新的速度长度等于[c×k]个交换序列，即选择交换序列表中的前[c×k]个交换序；对于大于1的实数c，令c＝k+c’，c′∈(0,1)，k为整数，则

[0310] 基于上述更新策略，每次迭代对种群中的粒子进行更新。

[0311] 9.4)每次迭代更新粒子之后，对粒子进行Prufer解码，使其还原为拓扑结构。

[0312] Prufer解码的过程如下：

[0313] a)向量U、P、Q分别为节点编号的全集、Prufer数列和未出现在P中的节点编号；从Q中选择编号最小的节点u，从P中选择最左边的节点w，将线路(u,w)加入集电系统的拓扑中。

[0314] b)从Q中删除u，从P中删除w；如果节点w不再出现在P中，则将w加入Q中。

[0315] 重复步骤a)-b)，直到P为空集为止；当P为空集时，仍有两个节点u0和w0留在Q中；将线路(u0,w0)加入集电系统拓扑中；至此，Prufer解码完成，Prufer数列被还原为拓扑。

[0316] 9.5)根据解码完成后得到的拓扑结构，从拓扑末端的机组开始，沿线路向变电站侧推算，根据相应机组的出力期望值计算出每条线路中传输的功率，根据功率确定该条线路使用的电缆类型，进而计算集电系统的各项成本。

[0317] 9.6)集电系统优化的迭代。

[0318] 令集电系统规划的最大迭代次数为imax，当相邻两次迭代的集电系统总成本CECS差的绝对值小于临界阈值ε1，或迭代次数大于imax时迭代停止，该布局下的最佳集电系统拓扑及电缆选型即已得到。若不满足上述两个条件之一，则继续迭代。将得到的最佳拓扑下的集电系统总成本CECS储存在公式(7)中的CECS中，用于下一次机组布局的迭代。

[0319] 9.7)机组布局优化的迭代

[0320] 机组布局种群中的所有粒子均得到最优集电系统之后，即完成了一次机组布局迭代。令机组布局优化的最大迭代次数为Imax，则当相邻两次迭代的LCOE之差的绝对值小于邻接阈值ε2，或迭代次数大于Imax时迭代停止。若不满足上述两个条件之一，则继续迭代。

[0321] 实施例2：

[0322] 在考虑规划海域潮汐资源储量和潮汐流速的基础上，在我国东南沿海地区建设一座潮汐流能发电场，具体实施步骤如下：

[0323] 1)输入参数：输入90天每天24个时段的实测潮汐流速幅值vtm与流速方向vtd，(tm,td＝1,2,…,2160)，且流速方向与流速幅值一一对应；潮汐流速方向分为12组。潮汐流能发电机组的切入流速vin＝0.7m/s、额定流速vrated＝3.0m/s、切出流速vout＝5.0m/s、额定输出功率Prated＝1.5MW、获能系数Cp＝0.45、推力系数CT＝0.89、叶片直径D＝18m。海水密度ρ＝3 3
1.025×10kg/m。中压电缆数据如表1所示：

[0324] 表1中压电缆参数表

[0325] 电缆编号 电缆类型 最大载流量(A) 单位电阻(mΩ/m)1 1/0 150 0.5482
2 4/0 211 0.2741
3 500kcmil 332 0.1184
4 750kcmil 405 0.0813
5 1000kcmil 462 0.0633

[0326] 注：1kcmil＝0.507mm2

[0327] 潮汐流能发电场长方形规划海域的长xmax＝700m、宽ymax＝400m、机组最小安全距离Dmin＝54m；中压集电系统的额定电压U＝33kV；海底电缆最大允许压降ΔUmax＝0.33kV。粒子群算法种群规模Np＝50、惯性权重w＝0.8、学习因子l1＝2，l2＝2、最大迭代次数Imax＝100。

[0328] 潮汐流能发电机组单价Cunit＝$799712、发电机组配套小型变压器单价CTrans＝$7
35920、机组安装费用Cinstall＝$386181、海上变电站总价csub＝$1.05×10 、高压海底电缆单价chv＝755$/m、设备折现率r＝0.03、设备经济寿命L＝20年、电缆单位挖掘成本ct＝50$/m、单位网损成本cp＝0.52$/kWh。

[0329] 2)规划海域流速特性评估方法

[0330] 2.1)流速方向的分组：将[0,360°]的方向范围平均分为12组，即：[0,30°),[30,60°),…,[300,330°),[330,360°)；按照流速方向数据vtd，将对应的vtm分入相应方向范围内；每个方向范围内的概率记为Pi。

[0331] 2.2)采用k均值聚类算法，对同一流速方向内流速幅值进行聚类。各个典型流速的数值及概率如表2和表3所示。其中，第一列代表分组范围，第一行为每一组内的典型流速编号；

[0332] 2.3)根据上述概率及典型流速数据，可以做出如附图1所示的潮汐玫瑰图，其中同心圆代表不同方向下流速出现的概率。

[0333] 表2典型流速数值(单位：m/s)

[0334]  1 2 3 4 5
[0,30°) 0.0766 0.1614 0.2901 — —
[30,60°) 0.0892 0.1811 0.3025 — —
[60,90°) 0.2854 0.7892 1.5562 1.9097 —
[90,120°) 0.3392 0.7796 1.2880 1.7764 2.2099
[120,150°) 0.2559 0.5864 0.9314 1.2480 —
[150,180°) 0.1066 0.2363 0.3839 — —
[180,210°) 0.1018 0.2366 0.3768 — —
[210,240°) 0.1192 0.2496 0.4128 — —
[240,270°) 0.1552 0.3484 0.5340 — —
[270,300°) 0.4619 1.1378 1.9914 2.7549 3.4329
[300,330°) 0.4648 1.2665 2.0403 2.8344 3.5177
[330,360°) 0.0829 0.1881 0.3097 — —

[0335] 表3典型流速的概率(％)

[0336]  1 2 3 4 5 合计
[0,30°) 0.2064 0.1559 0.0367 — — 0.3990
[30,60°) 0.1880 0.1972 0.0459 — — 0.4311
[60,90°) 1.4677 0.3302 0.0596 0.0459 — 1.9034
[90,120°) 3.9215 9.0492 10.1087 10.1041 2.3437 35.5272
[120,150°) 2.0181 1.5640 0.2798 0.0046 — 3.8665
[150,180°) 0.3853 0.4678 0.0963 — — 0.9494
[180,210°) 0.4036 0.3211 0.0963 — — 0.8210
[210,240°) 0.3669 0.3807 0.0642 — — 0.8118
[240,270°) 0.6559 0.6375 0.0321 — — 1.3255
[270,300°) 6.8477 5.1277 4.2838 3.5775 1.5044 21.3411
[300,330°) 5.0727 7.4806 9.9252 7.3247 2.0731 31.8763
[330,360°) 0.2798 0.3440 0.0688 — — 0.6926

[0337] 3)潮汐流能发电场的定容

[0338] 3.1)规划海域潮汐能量密度Pdensity＝1.4855kW/m2；

[0339] 3.2)规划海域潮汐资源储量计算

[0340] 潮汐资源储量Pc计算方法如下式：

[0341] Pc＝kexmaxymaxPdensity   (39)

[0342] 式中，ke＝0.15为环保系数；xmax＝400m与ymax＝700m分别为发电场规划海域的长和宽；计算可得Pc＝62.391MW；

[0343] 3.3)潮汐流能发电场机组数量的确定：

[0344]

[0345] 式中，Prated＝1.5MW为单台机组的额定容量，Pc＝62.391MW；计算可得发电场内的机组总数NTCT＝41台。

[0346] 4)建立潮汐流能发电场优化规划模型

[0347] 4.1)优化模型的目标函数

[0348] 潮汐流能发电场规划模型的目标函数如下式所示：

[0349]

[0350]

[0351] CECS＝C1+C2+C3+C4+C5   (43)

[0352] ET＝PT×24   (44)

[0353] 式中，NTCT＝41为发电场内的机组总数，Cunit＝$799712为机组单价，CTrans＝$35920为机组配套变压器的单价，Cinstall＝$386181为机组安装费用，r＝0.03和L＝20年分别为设备的折现率和经济寿命；C1、C2、C3、C5均基于电缆的具体长度求得；nl＝40代表中压电缆条数，kv＝44代表典型流速的数量；cf.mv为第f类中压电缆的单价，f＝1,2,…,5，数据可见表1，ct＝50$/m为单位长度电缆的运输与铺设成本；cp为单位电量损耗成本，csub＝$1.05×107为海上升压变电站总价，chv＝705$/m为高压海底电缆的单价；ls为第s条中压电缆的长度，s＝1,2,…,40，需根据具体拓扑确定，lhv为高压电缆的长度，需根据具体布局确定；U＝33kV为中压集电系统额定电压，rf为第f类中压电缆的单位长度电阻，具体数据见表1，pt为第t类典型流速出现的概率，具体数据见表3。

[0354] 4.2)确定优化模型的约束条件。

[0355] 5)潮汐流能发电机组与海上变电站初始位置生成。

[0356] 5.1)初始位置生成

[0357] 采用粒子群算法的粒子表示潮汐流能发电场中每台机组和海上变电站的坐标；粒子种群中共有Np＝50个粒子，其中第d个粒子可由下式表示：

[0358] Xd＝[x1,y1,x2,y2,…,x41,y41,x42,y42],d＝1,2,...,50   (45)

[0359] 式中，(x1,y1)、(x2,y2)、(x41,y41)分别代表第1台、第2台、第41台机组的横纵坐标，(x42,y42)代表海上变电站的横纵坐标；

[0360] 基于上式，随机生成50个粒子的坐标，即构成初始粒子种群；

[0361] 5.2)粒子位置约束条件的判断

[0362] 得到初始种群后，依次判断粒子种群中50个粒子是否满足约束条件(17)-(20)，若满足，则保留该粒子，若不满足，则舍弃该粒子并重新生成，直到种群中所有粒子均满足约束条件为止。

[0363] 6)计算发电场日发电量

[0364] 为了得到粒子群中每个粒子所确定的发电场的日发电量，需要计算潮汐玫瑰图中每一典型流速下各台机组的输出功率，获取每台机组的期望输出功率；可按照如下步骤计算；

[0365] 6.1)根据流速方向调整机组坐标

[0366] 计算尾流时需要判断机组的上下游位置关系，为了简化判断方法，将规划区域及机组放在笛卡尔坐标系中，计算潮汐玫瑰图中任意典型流速下的尾流分布时，首先按照下式调整机组坐标，使流速方向沿x轴正半轴：

[0367]

[0368] 机组坐标调整完成后，由于潮汐流从x负方向流向x正方向，因此可按照机组横坐标的大小判断机组的上下游位置关系；

[0369] 6.2)重合面积的计算

[0370] 第m台机组产生的尾流与第n台机组叶片的重合面积Amn可由下式计算：

[0371]

[0372]

[0373] 式中，ymn为第m台和第n台机组纵坐标之差的绝对值，θw,mn和θ0,mn所示角度见附图2；r0＝9m为第n台机组的叶片半径；σmn为第m台机组产生的尾流区域在第n台机组处的半径，可以由下式计算：

[0374]

[0375] 式中，I0＝0.1代表湍流强度，xmn代表第m台机组和第n台机组横坐标之差的绝对值，c(I0)可以由下式计算：

[0376]

[0377] 根据式(47)-(50)，可以计算出任一上游机组产生的尾流区域与下游机组的重叠面积；

[0378] 6.3)各典型流速下各台机组流速的计算

[0379] 在潮汐玫瑰图中各个典型流速下，将典型流速设为自由流速V0，第n台机组处的流速vn可以按下式计算：

[0380]

[0381] 式中，W(n)为第n台机组的上游机组集合，Nn为集合W(n)中的元素数量，A0＝254.47m2为第n台机组叶片扫过的面积，CT＝0.89为机组的推力系数；

[0382] 6.4)各典型流速下机组出力的计算

[0383] 基于上述尾流模型，可以求得不同典型流速下的各台机组处的流速，进而可根据下式求得第i个流速方向分组中的第j个典型流速下，第n台机组的输出功率Pn,ij：

[0384]

[0385] 式中，vn,ij为在第i个流速方向分组中的第j个典型流速下，第n台机组处的流速；vin＝0.7m/s，vout＝5.0m/s和vrated＝3.0m/s分别为潮汐流能发电机组的切入流速、切出流速和额定流速；Cp＝0.45为机组的获能系数；ρ＝1.025×103kg/m3为海水密度；Prated＝1.5MW为第n台机组的额定功率；Pn,ij的概率为pij，具体数据见表3；

[0386] 6.5)机组与发电场出力期望的计算

[0387] 第n台机组出力的期望值可由下式计算：

[0388]

[0389] 式中，gd＝12为流速方向分组数，kd为第i组内流速的类数；

[0390] 潮汐流能发电场的出力期望值可由下式计算：

[0391]

[0392] 式中，NTCT＝41为潮汐流能发电场中机组的总数；

[0393] 根据输出功率期望值，即可计算出优化模型目标函数中的发电场日发电量ET。

[0394] 7)机组及变电站坐标的更新

[0395] 得到各个粒子确定的发电场日发电量后，更新各个粒子的最佳位置，以及种群的最佳位置，并按照如下两个公式更新粒子位置：

[0396]

[0397]

[0398] 式中， 和 分别为第d个粒子在第t代的移动速度和位置， 和 为第d个粒子在第(t+1)代时的移动速度和位置， 为第d个粒子到第t代为止的最佳位置，globalt为种群到第t代为止取得的最佳位置，w＝0.8为惯性因子，l1＝2、l2＝2为学习因子，r1、r2为[0,1]之间的随机数。

[0399] 8)潮汐流能发电机组的分组

[0400] 8.1)确定分组数

[0401] 由中压海底电缆的最大载流量与机组的额定电流确定每一组内机组数上限nmax，计算公式如下：

[0402]

[0403] 式中，Imv,max＝486A为中压电缆的最大载流量，Irated＝26.24A为每台机组的额定电流；计算可得nmax＝17；

[0404] 机组的分组数Ng可由如下公式计算：

[0405]

[0406] 式中，NTCT＝41为发电场内的机组总数；计算可得Ng＝3；

[0407] 8.2)计算各台机组的夹角

[0408] 以更新后的粒子所确定的机组与变电站坐标为操作对象，以变电站为原点建立笛卡尔坐标系；定义机组的夹角为：机组与变电站的连线与x正半轴的夹角；计算各台机组的夹角，并令第m台机组的夹角为θm(θm∈[0,360°]，m＝1,2,…,41)；

[0409] 8.3)机组分组与调整

[0410] 随机生成Ng＝3条以原点为顶点的射线作为聚类中心，第g个聚类中心的角度为θg(θg∈[0,360°]，g＝1,2,3)；分别计算各台机组与各聚类中心的夹角，并将机组归入夹角最小的组中；

[0411] 各台机组归入相应的组后，判断是否有机组数越限情况产生；若第g组内的机组数ng>17，则在第g组中，选择与第g个聚类中心夹角最大的(ng-17)台机组，将其分入夹角第二小的组内，直到所有组均不越限为止；

[0412] 对分入各个组内的机组与x正半轴的夹角值求平均，得到新的聚类中心；重复分组与调整过程，直至聚类中心不再变化；至此，一个粒子所确定的发电场即分组完毕；不妨设粒子Xd被划分为Xd1,Xd2,Xd3共3个组。

[0413] 9)集电系统拓扑规划及电缆选型

[0414] 对于完成分组的粒子，各个组均按如下步骤规划集电系统拓扑和电缆选型；现以布局Xd1为例介绍详细步骤；

[0415] 9.1)生成初始拓扑

[0416] 针对布局Xd1，随机产生Nt＝50个生成树，连接各台机组，用以组成初始种群，并保证初始种群的线路不交叉；

[0417] 9.2)Prufer编码

[0418] 对于各个已经连接完成的拓扑，采用Prufer数列对拓扑进行编码；将拓扑编码为Prufer数列的过程称为Prufer编码；假设一随机生成的拓扑如附图3所示，其中方块代表海上变电站，圆圈代表机组，线段代表海底电缆，阿拉伯数字为机组/变电站编号，罗马数字为线路的编号；Prufer编码的步骤如下：

[0419] a)定义三个向量U、P和Q，其中，U包含本组内的所有节点，即U＝[1,2,3,4,5,6,7]，P用来储存Prufer编码，初始为空，Q代表没有包含在P中的节点，初始Q＝U；

[0420] b)首次操作时，从U中找出只与一个节点相连的节点，即节点1、3、6、7，并选择其中编号最小的节点1定义为叶子节点；

[0421] c)将节点1所连接的节点2放入向量P中，并从Q中删去叶子节点1，从拓扑中删去线路(1,2)；

[0422] d)重复步骤b)-c)，直到拓扑中仅剩一条线路为止(即线路(4,7))；此时，向量P由5个元素组成，将U内未在向量P中出现的所有元素放入Q中；

[0423] 编码后向量分别为P＝[2,2,4,5,4]，Q＝[1,3,6,7]；

[0424] 9.3)拓扑粒子更新

[0425] 由于拓扑粒子中的元素均为正整数，因此采用如下更新策略对拓扑粒子进行更新：

[0426]

[0427]

[0428] 式中， 和 为第h个粒子在第t代时的位置和速度(h＝1,2,…,50)， 和为第h个粒子在第t+1代时的位置和速度， 为第h个粒子到第t代为止的最佳位置，globalt为种群到第t代为止取得的最佳位置，w’＝0.5为惯性因子，l′1＝2、l′2＝2为学习因子，r′1、r′2为[0,1]之间的随机数；

[0429] 假设上式中 globalt＝(2,3,4,1)，r1＝0.5，r2＝0.5，更新策略中的概念与运算过程如下：

[0430] a)交换序：粒子速度 由两个交换序(1,4)、(3,1)组成；

[0431] b)减法操作Θ：

[0432] c)乘法操作 由于w’＝0.8， 中共有k＝2个交换序，因此 的长度为[w’×2]＝2个交换序，即 同理，

[0433] 9.4)Prufer解码

[0434] 每次迭代更新粒子之后，对粒子进行解码，使其还原为拓扑结构；将Prufer数列还原为拓扑的过程称为Prufer解码；同样以附图3为例，Prufer解码的过程如下：

[0435] a)向量U、P、Q分别为节点编号的全集、Prufer数列和未出现在P中的节点编号，且U＝[1,2,3,4,5,6,7]，P＝[2,2,4,5,4]，Q＝[1,3,6,7]；

[0436] b)首次操作时，从Q中选择编号最小的节点u＝1，从P中选择最左边的节点w＝2，将线路(1,2)加入集电系统的拓扑中；

[0437] c)从Q中删除u＝1，从P中删除w＝2；如果节点w不再出现在P中，则将w加入Q中；

[0438] d)重复步骤b)-c)，直到P为空集为止；当P为空集时，仍有两个节点u0＝4和w0＝7留在Q中；将线路(4,7)加入集电系统拓扑中；至此，Prufer解码完成，Prufer数列被还原为拓扑；

[0439] 9.5)集电系统成本的计算

[0440] 根据解码完成后得到的拓扑结构，从拓扑末端的机组开始，沿线路向变电站侧推算，根据相应机组的出力期望值计算出每条线路中传输的功率，根据功率确定该条线路使用的电缆类型；仍以附图3所示拓扑为例，表4列出了每条线路中流过的电流(每台机组的出力均记为额定功率)以及选取的与表1相对应的电缆型号；由此，可以计算出集电系统的中压电缆成本、挖掘成本和网损成本，至此集电系统规划的一次迭代即已完成；

[0441] 表4附图3电缆选型

[0442]线路编号 传输电流(A) 电缆型号
I 157.46 2
II 26.24 1
III 104.97 1
IV 52.48 1
V 26.24 1
VI 26.24 1

[0443] 9.6)集电系统优化的迭代

[0444] 令集电系统规划的最大迭代次数为imax＝100，则当相邻两次迭代的集电系统总成本CECS差的绝对值小于临界阈值ε1＝1×10-6，或迭代次数大于imax时迭代停止，该布局下的最佳集电系统拓扑及电缆选型即已得到；若不满足上述两个条件之一，则继续迭代；

[0445] 将得到的最佳拓扑下的集电系统总成本CECS储存在公式(7)中的CECS中，用于下一次机组布局的迭代；

[0446] 9.7)机组布局优化的迭代

[0447] 机组布局种群中的所有粒子均得到最优集电系统之后，即完成了一次机组布局迭代；令机组布局优化的最大迭代次数为Imax＝100，则当相邻两次迭代的LCOE之差的绝对值小于临界阈值ε2＝1×10-6，或迭代次数大于Imax时迭代停止；若不满足上述两个条件之一，则继续迭代；

[0448] 达到收敛条件后，可得最优规划方案如附图4所示，其中图4(a)为机组布局方案，图4(b)为集电系统的拓扑；成本结果如表5所示，其中，C1为中压电缆日等效成本，C2为中压电缆的挖掘和维护日等效成本，C3为中压电缆日网损成本，C4为海上变电站的日等效成本，C5为高压电缆的日等效成本，ET为发电场的日发电量，Ctotal.TCT为机组的日等效投资成本，Lmv为中压电缆的总长，LCOE为发电场的平准化发电成本。

[0449] 表5发电场规划成本

[0450]