CN117154808A - 考虑电力交易的分布式新能源和多元储能聚合建模方法 - Google Patents
考虑电力交易的分布式新能源和多元储能聚合建模方法 Download PDFInfo
- Publication number
- CN117154808A CN117154808A CN202310614453.2A CN202310614453A CN117154808A CN 117154808 A CN117154808 A CN 117154808A CN 202310614453 A CN202310614453 A CN 202310614453A CN 117154808 A CN117154808 A CN 117154808A
- Authority
- CN
- China
- Prior art keywords
- energy storage
- energy
- power
- data
- wind
- 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
- 238000004146 energy storage Methods 0.000 title claims abstract description 140
- 238000000034 method Methods 0.000 title claims abstract description 99
- 230000002776 aggregation Effects 0.000 title claims abstract description 35
- 238000004220 aggregation Methods 0.000 title claims abstract description 35
- 238000009826 distribution Methods 0.000 claims abstract description 39
- 238000010248 power generation Methods 0.000 claims abstract description 38
- 238000006243 chemical reaction Methods 0.000 claims abstract description 29
- 239000001257 hydrogen Substances 0.000 claims abstract description 29
- 229910052739 hydrogen Inorganic materials 0.000 claims abstract description 29
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims abstract description 28
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 claims abstract description 19
- 239000002253 acid Substances 0.000 claims abstract description 19
- 229910052744 lithium Inorganic materials 0.000 claims abstract description 19
- 239000006185 dispersion Substances 0.000 claims abstract description 14
- 238000005516 engineering process Methods 0.000 claims abstract description 7
- 239000013598 vector Substances 0.000 claims description 22
- 230000004044 response Effects 0.000 claims description 18
- 230000008901 benefit Effects 0.000 claims description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 11
- 238000005286 illumination Methods 0.000 claims description 11
- 230000009467 reduction Effects 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000013507 mapping Methods 0.000 claims description 6
- 238000013528 artificial neural network Methods 0.000 claims description 5
- 230000000087 stabilizing effect Effects 0.000 claims description 5
- 210000002569 neuron Anatomy 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 238000006116 polymerization reaction Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 230000002159 abnormal effect Effects 0.000 claims description 2
- 238000003062 neural network model Methods 0.000 claims description 2
- 238000012847 principal component analysis method Methods 0.000 claims description 2
- 238000007670 refining Methods 0.000 claims description 2
- 238000013459 approach Methods 0.000 claims 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 47
- 239000007789 gas Substances 0.000 description 36
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 26
- 230000005611 electricity Effects 0.000 description 21
- 239000000446 fuel Substances 0.000 description 15
- 230000008569 process Effects 0.000 description 15
- 230000006870 function Effects 0.000 description 13
- 210000004027 cell Anatomy 0.000 description 12
- 238000010586 diagram Methods 0.000 description 12
- 238000010438 heat treatment Methods 0.000 description 12
- 230000005855 radiation Effects 0.000 description 12
- 238000004458 analytical method Methods 0.000 description 10
- 230000008878 coupling Effects 0.000 description 10
- 238000010168 coupling process Methods 0.000 description 10
- 238000005859 coupling reaction Methods 0.000 description 10
- 230000005540 biological transmission Effects 0.000 description 9
- 238000007599 discharging Methods 0.000 description 9
- 238000005457 optimization Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 8
- 238000003860 storage Methods 0.000 description 8
- 239000002351 wastewater Substances 0.000 description 8
- 230000000694 effects Effects 0.000 description 7
- 239000003345 natural gas Substances 0.000 description 7
- 238000012876 topography Methods 0.000 description 7
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 6
- 230000008859 change Effects 0.000 description 6
- 238000011084 recovery Methods 0.000 description 6
- 230000009897 systematic effect Effects 0.000 description 6
- 238000005868 electrolysis reaction Methods 0.000 description 5
- 238000012546 transfer Methods 0.000 description 5
- 229910002092 carbon dioxide Inorganic materials 0.000 description 4
- 239000001569 carbon dioxide Substances 0.000 description 4
- 230000001419 dependent effect Effects 0.000 description 4
- 238000005485 electric heating Methods 0.000 description 4
- 238000013178 mathematical model Methods 0.000 description 4
- 239000003595 mist Substances 0.000 description 4
- 239000002689 soil Substances 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 238000001704 evaporation Methods 0.000 description 3
- 230000008020 evaporation Effects 0.000 description 3
- 238000007726 management method Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 238000007621 cluster analysis Methods 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 2
- 238000001816 cooling Methods 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 230000005674 electromagnetic induction Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000003631 expected effect Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 238000003064 k means clustering Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 239000010755 BS 2869 Class G Substances 0.000 description 1
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000000418 atomic force spectrum Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- -1 electricity Substances 0.000 description 1
- 238000012983 electrochemical energy storage Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000020169 heat generation Effects 0.000 description 1
- 238000005338 heat storage Methods 0.000 description 1
- 150000002431 hydrogen Chemical class 0.000 description 1
- 238000005984 hydrogenation reaction Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000009347 mechanical transmission Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000013486 operation strategy Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000005406 washing Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/24—Querying
- G06F16/245—Query processing
- G06F16/2455—Query execution
- G06F16/24553—Query execution of query operations
- G06F16/24554—Unary operations; Data partitioning operations
- G06F16/24556—Aggregation; Duplicate elimination
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/27—Replication, distribution or synchronisation of data between databases or within a distributed database system; Distributed database system architectures therefor
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
- G06Q40/00—Finance; Insurance; Tax strategies; Processing of corporate or income taxes
- G06Q40/04—Trading; Exchange, e.g. stocks, commodities, derivatives or currency exchange
-
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/008—Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/28—Arrangements for balancing of the load in a network by storage of energy
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/22—The renewable source being solar energy
- H02J2300/24—The renewable source being solar energy of photovoltaic origin
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/28—The renewable source being wind energy
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/40—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation wherein a plurality of decentralised, dispersed or local energy generation technologies are operated simultaneously
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Economics (AREA)
- Evolutionary Computation (AREA)
- Strategic Management (AREA)
- Power Engineering (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Databases & Information Systems (AREA)
- Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Artificial Intelligence (AREA)
- Tourism & Hospitality (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computing Systems (AREA)
- Finance (AREA)
- Accounting & Taxation (AREA)
- Development Economics (AREA)
- Software Systems (AREA)
- Computational Linguistics (AREA)
- Quality & Reliability (AREA)
- Mathematical Physics (AREA)
- Entrepreneurship & Innovation (AREA)
- Computer Hardware Design (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Geometry (AREA)
- Technology Law (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
Abstract
本发明涉及一种考虑电力交易的分布式新能源和多元储能聚合建模方法,包括以下步骤:步骤1、整合并分析现有配电网、用户侧系统中获得的数据,研究风电、光伏发电的聚合预测,构建基于数据驱动的用户侧风/光新能源的聚合预测模型;步骤2、考虑多类型能源转换设备的非线性特性,对锂电、铅酸、氢储能等技术中常用储能的特性进行分析,建立锂电、铅酸、氢储能等多元储能数据驱动聚合模型;步骤3、分析本地区微气象分布特征和光伏分散特征,构建基于数据驱动的分布式能源的多元储能混合运行模型。本发明构建基于数据驱动的分布式新能源和多元储能聚合模型,为用户侧分布式新能源和多元储能的网源协调奠定基础。
Description
技术领域
本发明属于配电网用户侧网源协调领域,尤其是一种考虑电力交易的分布式新能源和多元储能聚合建模方法。
背景技术
分布式新能源的接入对电网稳定性产生诸多不确定性影响,尤其分布式新能源接入系统的影响在地区级电网中影响已经有所表现。然而,分布式新能源电源接入系统影响的研究还停留在表面定性研究层面,以及脱离实际测试数据的空洞简单模型上,导致相关研究与仿真结果脱离实际。
发明内容
本发明的目的在于考虑分布式新能源和多元储能的不同运行特性,构建基于数据驱动的分布式新能源和多元储能聚合模型,为用户侧分布式新能源和多元储能的网源协调奠定基础。
本发明的技术方案:
一种考虑电力交易的分布式新能源和多元储能聚合建模方法,包括以下步骤:
步骤1、整合并分析现有配电网、用户侧系统中获得的数据,研究风电、光伏发电的聚合预测,研究分布式光伏发电的电力交易,构建基于数据驱动的用户侧风/光新能源的电力交易聚合预测模型;
步骤2、考虑多类型能源转换设备的非线性特性,对锂电、铅酸、氢储能储能的特性进行分析,建立锂电、铅酸、氢储能多元储能数据驱动聚合模型;
步骤3、分析本地区微气象分布特征和光伏分散特征,构建基于数据驱动的分布式能源电力交易和多元储能混合运行模型。
优选的,所述步骤1的具体步骤包括:
步骤1.1:通过对移动互联网、智慧电网、用户侧系统获得的数据进行整合和提炼,获得电网拓扑结构、潮流计算、负荷预测数据;
步骤1.2:分析风电、光伏发电特性,研究风电、光伏的聚合预测;
步骤1.3:建立分布式光伏发电的电力交易模型,分析分布式电源参与电力交易的优势;
步骤1.4:运用数据驱动方法构建基于人工神经网络的用户侧风/光新能源的聚合预测模型。
优选的,聚合预测方法包括以下步骤:
首先提取历史数值天气中的气象因素记录,并将记录中的各个因素数据进行标准化处理;再以预测时间点的气象为核心对标准化历史气象记录进行相似聚合,筛选原始样本;然后对相似聚合中记录的气象因素进行降维,得到独立因素;最后建立独立因素和风电/光伏功率的映射关系,实现聚合预测。
优选的,气象因素包括:包括风速、风向、光照强度、气压、温度和湿度气象因素。
优选的,所述步骤1.2包括以下步骤:
步骤1.2.1:提取历史数值天气中的气象因素记录,并将记录中的各个因素数据进行标准化处理,剔除异常点并将数据线性映射到[0,1]区间内;
步骤1.2.2:以预测时间点的气象为核心对标准化历史气象记录进行相似聚合;
步骤1.2.3:对相似聚合中记录的气象因素运用主成分分析法进行降维,得到独立因素;
步骤1.2.4:建立独立因素和风电/光伏功率的映射关系,实现聚合预测。
优选的,所述步骤2的具体步骤包括:
步骤2.1:考虑多类型能源转换设备的非线性特性,对锂电、铅酸、氢储能等储能技术的启停响应速度、功率密度、能量转化效率以及成本特性进行分析总结;
步骤2.2:利用量测装置获得锂电、铅酸、氢储能等储能设备的历史能量流动数据,采用神经网络模型拟合上述储能设备的历史数据集,得到初步的储能设备非线性黑盒模型;
步骤2.3:根据具体选用的锂电池、铅酸电池、氢储能设备,通过对步骤2.2得到的初步的储能设备非线性黑盒模型进行最大荷电量额定值、极限功率等参数的调整,建立相应的储能模型。
优选的,神经网络算法是由输入层、隐含层以及输出层组成,输入层使用n个神经元,输入信号使用向量X=[x1,x2,...,xn]Τ;隐含层包含n个神经元节点,隐含层输出向量为G=[g1,g2,...,gn]Τ;输出层中包含的向量为D=[d1,d2,...,dn]Τ。
优选的,所述步骤3的具体步骤包括:
步骤3.1:掌握本地区的微气象分布特征和该区域内分散光伏特征;
步骤3.2:对数据进行建模,并利用K-中心点方法对微气象分布特征和光伏分散特征进行聚类;
步骤3.3:基于需求响应策略,以分布式能源出力为约束条件,研究计及需求响应策略的分布式能源优化运行模型;
步骤3.4:基于多元储能设备平抑风光电源出力的能量调度原则,将构建的分布式能源与多元储能的数据驱动模型进行串联,得到当地风光资源分布特性下,风光出力预测与多元储能运行特性的动态模型,构建基于数据驱动的分布式能源和多元储能混合运行模型。
优选的,掌握本地区的微气象分布特征和该区域内分散光伏特征的途径包括文献检索以及通过气象局或电网公司获得。
优选的,所述K-中心点方法为:
在n个数据对象中随机选择k个中心,集合中其他数据对象被划分到与中心最近的数据集合中去,形成了k个初始的数据集合;对于每个数据集合,计算新的中心,然后根据新的中心重新分配其他数据集合,这样不断迭代,直到所有的数据中心不再变化;该算法的具体步骤:
(1)随机选择k个初始聚类中心C={c1,c2…ck};
(2)如果s与聚类中心ci之间的距离小于该对象到cj的距离,则将s划分到聚类簇Ci中;
(3)重新计算聚类簇Ci的中心:/>重复步骤(2)和(3)直到聚类中心C不再变化。
本发明的有益效果:
1、本发明构建了基于数据驱动的分布式新能源和多元储能聚合模型,充分考虑风/光分布式新能源、水/电/气/热等多类型能源转换设备的非线性特性以及锂电/铅酸/氢储能等多元储能的运行特点,基于数据驱动构建了更全面的聚合预测模型,符合用户侧网源接入系统的实际运行情况。
2、本发明考虑到分布式能源和多元储能单元单独运行以及混合运行的不同特性,并且利用聚类技术分析运行特点,建立了区域范围内分布式能源和多元储能混合运行特性的模型,为分析用户侧网源接入对电网稳定性的不确定影响和确定网源协调运行策略提供模型基础。
附图说明
图1聚合预测方法流程图;图2神经网络算法模型图;图3K-中心点方法算法流程图;图4为本发明实施例预测精度与真实值图;图5为本发明实施例多元储能系统的特性模型结果图;图6为本发明实施多元储能系统的能效特性模型结果图;图7为本发明实施多元储能系统的运行精度图;图8本发明实施例风力发电机组模型图;图9为本发明实施例光伏电池原理图;图10为本发明实施例光伏发电系统结构图;图11为本发明实施例区域水电气热综合能源系统能量流通示意图;图12为本发明实施例电转气设备运行过程图;图13为本发明实施例热电联产机组工作原理图;图14为本发明实施例风光储联合发电系统结构图;
图8中:v—自然风速;T—风轮的转矩;ω—风轮的角速度;β—桨距角,控制有功出力的大小;Pset—为满足电网运行要求制定的当前风光系统的设定功率参考值;Pseti—第i台风电机组的设定有功参考值;Pref—每台风电机组预期发的有功。
具体实施方式
如图1-3所示,本发明的具体实施例如下:
1.1风电场聚合预测模型
风力发电的原理是:自然风首先吹向机组,使得叶轮转动,然后带动轮毂转动,风力机组的叶轮将空气流动所产生的动能转化为机械能,其次由传动结构将转化后的机械能传递到转子上,使得转子旋转发出电能,最后通过变电站将电能注入电网。其能量转换过程是:风能→机械能→电能。
风电机组的结构模型主要包括:模拟风速输入的风速模型、模拟风能转换为机械能的风力机模型、模拟传动装置的机械传动模型、模拟机械能转换为电能的风力发电机模型和风力发电机组功率控制模型,机组各模型之间的传递关系如图8所示:
首先,风光储能量分配系统根据电网调度制定的调节和控制模式得到Pset,其次,根据当前区域的风速预测结果和预先制定的能量分配策略,得到Pseti,然后对风力机的桨距角和转矩进行控制,实现风电机组的有功控制,最终输出Pref。
风能转化为机械能的风力发电机模型为:
式中:P为风轮的机械输出功率;ρ为空气密度;η为风能利用系数;S为叶片扫过的面积。
风力发电机输出的有功功率与风速的大小密切相关,所以在确定其输出功率前,对风速变化以及分布情况的研究至关重要。在实际应用中,一般有两种方法可以得到风速。一种是直接法,即直接参考某实际电厂或电站某一段时间的风速,这样所得到的数据较为准确,但可能达不到预期的效果;第二种是间接计算法,由于风速具有高度的随机性,所以我们可通过风速的概率分布来描述其统计特征,这类方法在实际应用中的研究较为广泛。
模拟风速分布的模型很多,主要有皮尔逊分布、瑞利分布和双参数威布尔(Weibull)分布等,但大量的实验数据表明精度较高的模型是Wuibull分布,其概率密度函数表达式为:
式中:v—实际风速;k—威布尔形状参数,描述函数的形状;c—威布尔尺度系数,
体现某时段的平均风速。
本文基于Weibull分布,提出使用Monte Carlo方法来模拟实际风速。
对式(1.2)进行积分,得到v的累积分布函数F(v)如式(1.3)所示:
对式(1.3)进行反演计算,则使用计算机模拟实际风速v的函数表达式为:
由于F(v)~(0,1),因此1-F(v)~(0,1),故随机抽样风速v:
当实时风速低于切入风速时,风力机组的有功出力为零;如果实时风速大于切入风速小于额定风速,风力机组的输出功率与风速成非线性关系;如果实时风速大于额定风速小于切出风速,风力机组输出额定功率,且保持不变;如果实时风速高于切出风速,风力机组的有功出力为零。
根据上述原理,可知风力发电机组的功率模型为输出功率与风速之间的分段函数:
式中:Pwt—风机额定出力;Pr—风机额定功率;vr—额定风速;vin—切入风速;vout—切出风速;
参数α、β、γ为风力发电机有功出力曲线的多项式拟合系数,其计算表达式如下所示:
1.2光伏聚合预测模型
光伏发电的原理是:我们利用电池中半导体的光生伏打效应,将辐射到电池板的太阳能转化成直流电能,再经过逆变器转化为交流电能后输出,最终注入电网中。其能量转换过程是:太阳能→电能。
其中,太阳能转换为电能的具体过程为:自然光首先照射到光伏电池上并被吸收,然后Р区和N区中的价电子在受到冲击后被释放,形成电子-空穴对。最后在内部E作用下分离,出现不同电量Q的积累,即产生光生U。如果此时接入负载,则会有Ⅰ流过,因此产生电能,这一过程可由图9表示。
光伏发电系统主要由太阳能光伏阵列、DC/DC变换器、控制器、逆变控制系统和滤波及升压电路组成,其组成结构及传递关系如图10所示。
光伏电池板通过串并联的方式组成光伏阵列,通过DC/DC变换器实现MPPT控制,然后将输出的直流电通过逆变器转换为交流电,接着经过滤波及升压电路,最终接入到三相交流电网中。
光伏发电系统的输出功率大小与光照辐照强度E,太阳能电池阵列总面积s及光电转换效率n有着直接的关系,其出力模型为:
Pst=ESη(1.10)
在电力系统的实际应用中,电池阵列总面积S及转换效率η是可以人为假定的,但是由于E是随机变化的,所以其输出功率也是随机变化的。因此,若要计算光伏发电的功率,首先要得到E的数值大小。
在实际应用中,一般有两种方法可以求得太阳光照强度。一种是直接法,即直接参考某实际电厂或电站某一段时间的太阳光照强度,这样所得到的数据较为准确,但可能达不到预期的效果;第二种是计算机模拟法,
由于光照强度具有高度的随机性,所以我们可通过计算机来模拟太阳光辐照度E,这类方法在工程实践中应用的较为广泛。
模拟太阳光辐照度E的分布很多,但贝塔分布使用最为广泛,其概率密度函数f(E)为:
式中:Em为太阳光辐照度E的最大值,λ和μ表示贝塔分布的两个形状参数。对式(1.11)积分太阳光辐照度E的概率累积分布函数F(E)如式(1.12)示:
对式(1.12)进行反演计算,则使用计算机模拟实际太阳光辐照度E的强度函数表达式为:
E=Em*F-1(E)(1.13)
由式(1.13)就可以得到太阳光辐照度E,从而求得光伏发电的输出功率。
1.3分布式电源参与电力交易模型
分布式电源的机组容量通常较小,并且发电曲线很难预测,这导致分布式电源直接参与电力市场交易在费用结算中会产生偏差考核费用。因此,分布式电源可探索由分布式电源运营商代理,聚合成1个主体参与电力市场竞争。假设聚合后的分布式电源电能量偏差结算至分布式电源运营商,其1天可以分为T个时间段,分布式电源运营商的偏差考核费用为:
式中:Rt、和ac分别为分布式绿电运营商在t时间段内的实际出力、合约电量和合约偏差系数;A、B分别为发生负偏差和正偏差时的惩罚系数;λθ为市场中平均成交价格。
分布式电源运营商(distributed generation operator,DGO)主要通过分布式光伏机组和分布式风机进行发电,在电力市场中进行电量交易获得相应的收益。在合作运行模式下,DGO不仅可以在电力市场中进行电量交易,也可以将电量卖给云储能运营商(energy storage operator,ESO)和柔性负荷聚合商(flexible load aggregator,FLA)来获得更高收益。因此1天内第个DGO主体的效益的最大化运行模型可表示为:
maxFi,DGO=Ui,g+Ui,CESO+Ui,FLA-Ci,DGO (1.16)
式中:Ui,g、Ui,ESO、Ui,FLA分别为第个DGO主体参与电力市场的收益、向ESO售电的收益、向FLA售电的收益;Ci,DGO为分布式光伏机组以及分布式风电机组的运行维护成本;为时间段内市场电能量出清价格;/>为绿色电力环境溢价;/>为分布式绿电参与电力市场的交易电量;/>分别为时段内DGO与ESO、FLA的交易电价,/>分别为时段内DGO与各个ESO、FLA的总交易电量;KW、Kpv分别为分布式风机和分布式光伏的单位发电量成本维护系数;/>分别为时段内DGO中分布式风机和分布式光伏的总发电功率;分别为t时段内DGO中分布式风机与光伏的发电量。当/>和当/>时,上述模型优化的结果为DGO单独参与市场交易时的收益。
1.4基于数据驱动的用户侧风/光新能源的聚合预测模型
风电和光伏受全年风速和光照的影响,其出力具有随机性、波动性和间歇性的特点,但又有一定规律性。风电和光伏的出力具有明显的季节性,如果直接将每个季节作为一个场景,场景划分的精确度不足,为此本文采用聚类方法为每个季度生成多个场景,提高场景精确度。首先提取表示风电、光伏出力特征的日特征值,然后把所有区域日特征值所构成的日特征向量作为聚类对象,进而根据改进的K-means算法进行场景缩减。风光出力曲线的日特征值包括日平均值、日最大值,日峰值时刻和日峰谷差。设共有n块聚合区域,d代表日期,则日特征向量可表示为:
式中,为第i块区域风电的日特征向量,/>为第i块区域光伏的日特征向量,ai,1,ai,2,ai,3和ai,4分别代表第i块区域风电和光伏的日平均值、日最大值,日峰值时刻和日峰谷差。
设某个季节共有D天,则每季节共有D个日特征向量即为D个初始场景,日特征向量是1×(n*4*2)维行向量。本文采用改进K-means算法对日特征向量
进行聚类分析,达到场景缩减的目的。
K-means是一种迭代聚类算法,通过不断更新聚类中心和分配给每类的对象来实现。为了确定最优聚类个数K,引入轮廓系数以评估聚类结果的好坏。轮廓系数的取值范围为[-1,1],越接近于1代表内聚度和分离度都相对较优。用a(i)代表第i个聚类向量的簇内不相似度,b(i)代表第i个聚类向量的簇间不相似度,则第i个聚类向量的轮廓系数为:
则通过改进K-means算法生成多场景的步骤如下:
(1)初始化日特征向量和聚类个数,日特征向量共有D个;
(2)从D个聚类对象中随机选取K个作为初始聚类中心,记为C1,0,C2,0,…Ck,0;
(3)对任意一个聚类对象,计算其到各个聚类中心的欧式距离,并将该对象归到距离最短的聚类中心所在类;
(4)利用均值法更新每类聚类对象的聚类中心,新的聚类中心记为C1,C2,…Ck若max||ci,0-ci||≤ε,则进行下一步,否则,更新聚类中心,返回(3)。其中ε是判断迭代是否终止的阈值;
(5)计算D个聚类对象所对应的轮廓系数,取D个轮廓系数的平均值作为聚类结果的轮廓系数,记为Sk;
(6)K=K+1,若K>9,进行下一步,否则返回(2);
(7)求取聚类结果轮廓系数的最小值Sk.min,则该K值对应最终聚类结果,并根据每类的聚类对象个数确定场景概率,记为Pk.min。
本文通过降低聚类向量的维度来降低计算复杂度。日出力曲线构成的聚类向量的维度为n*24*2,日特征值构成的聚类向量的维度为n*4*2,其中n代表聚合区域的个数,24是一天的小时数,2是资源的种类数。后者聚类向量的维度是前者的1/6,所以采用日特征向量进行聚类分析能降低计算复杂度。
2多元储能数据驱动的聚合建模
研究考虑水/电/气/热等多类型能源转换设备的非线性特性;研究锂电/铅酸/氢储能等快速响应机理数据驱动分析方法;研究锂电、铅酸、氢储能等多元储能数据驱动聚合建模技术。
2.1水/电/气/热等多类型能源转换设备的非线性特性
区域电-气-热-水综合能源系统是由电网、天然气网、热网和水网通过耦合元件耦合而成的联合供能系统。在综合能源系统中,耦合元件的作用是将某种形式的能源转换成另一种形式的能源,通俗一点来说就是需要以耦合设备为媒介,将电与冷、热、气、水等多能源建立更直观的联系,这样可将其他能源转为电力再去计算其成本,更直观也更有效。本文耦合设备主要包括电驱动水泵、废水利用装置、电转气设备(Power to Gas,PtG)、电锅炉、燃气锅炉、热电联产机组、水电联产机组等。以区域电-气-热综合能源系统能量流通为基础,结合水网模型及耦合设备可构成区域电-气-热-水综合能源系统能量流通示意图如图11所示:
(1)热电联产机组模型
热电联产机组(Combined heat and power,CHP)在工作时可以逐级利用各阶段能量,实现对电负荷、热负荷的联合供给,提高了能源利用率,减少碳排放量。其工作原理如图13所示。
其中热电联产机组燃烧天然气后,向区域综合能源系统所转移的电功率的数学模型如下:
式中,VMT,t—t时刻天然气的消耗量;RLHV—t时刻天然气的热值;ηMT—热电联产机组的发电效率。
其中热电联产机组产生的余热量HMT,t的表达式如下:
式中,ηL—机组的散热损失率。
热量传输效率模型:
HMT,h,t=ηMT×HMT,in,t(2.3)
式中,ηMT—交换器的传输效率;HMT,h,t—t时刻热交换器输入热能;HMT,in,t-t时刻热交换器输出热能。
热电联产机组的运行成本可以用下式表示:
式中,CMT,t—机组在1时刻的供能成本;cNG—天然气气价。
(2)供热锅炉
电制热锅炉(Electric Boiler,EB)主要以电流流通产生的损耗发热或电磁感应原理产生热能,进一步将可利用的热能通过热交换部分输出,进而供负荷侧用户使用。电制热锅炉与常规燃煤锅炉相比具有体积小、无污染等优点,同时在效率方面,电制热锅炉的电转热效率达98%以上。其数学模型可用下式表示:
HGB,t=Qgas,t×ηGB×RLHV (2.5)
Heb,t=Peb,t×ηeb (2.6)
Cgas,t=cNG×Qgas,t×Δt (2.7)
式中,HGB,t—t时段燃气锅炉t时刻所产生的热量;Heb,t—t时段电制热锅炉t时刻所产生的热量;ηGB—燃气锅炉制热效率;ηeb——制热锅炉制热效率;Cgas,t—燃气锅炉耗能成本;Qgas,t—t时段天然气的消耗量。
(3)电转气设备
电转气设备可以生成氢气和氧气,氢气也可以传输给加氢站,提供给燃料电池汽车使用,氢气也可以进一步合成甲烷,引入天然气管道,这样可以更大程度的消纳清洁能源,其运行过程如图12所示。该设备一方面可以利用能源供应富足时刻清洁能源产生的多余电能将水电解生成大量的氧气和氢气,还可以电转甲烷设备吸收二氧化碳,生成甲烷燃料和水。
在电解设备中,电解水生产的氢气可表示为:
式中,—电转气设备ptg的输入功率;/>—电转气设备ptg中电制氢气效率;—电转气设备ptg制取氢气量;ξe,g—电能转换为等值热能系数;HHVH2—氢气的高热值。
电转气设备中的电解水设备的功率通常被限制在一定的范围内:
式中,—电转气设备ptg中电解设备的最小和最大输入功率;/>—电转气设备ptg中电解设备的运行状态。
在甲烷化设备中,消耗的二氧化碳与合成的甲烷可分别表示为:
式中,—电转气设备ptg中消耗的二氧化碳;/>—电转气设备ptg中反应合成甲烷量;/>—氢气与二氧化碳之间的反应系数;/>—氢气与甲烷之间的反应系数。
(4)水电联产机组模型
由于水网系统的引入,能源系统耦合模型中新增了水泵、水电联产机组以及废水利用装置设备,其中水泵为水网中的独立元件,仅产生水资源供水,水电联产机组与废水利用装置属于水-电耦合装置。水电联产机组能够同时产生水功率与电功率,二者之间具有较强的约束,废水利用装置能够将多余的废水进行发电,产生电功率,能够实现水资源到电能的转换。以下介绍水电联产机组与废水利用装置两类耦合元件模型。
与热电联产机组类似,水电联产能够通过消耗一定的燃料,同时产生水能与电能。具体原理为:水电联产机组主要由废水利用装置与电力联产装置组成,水电联产机组通过燃烧燃料产生电能,同时将多余产生的水蒸气通过废水利用装置进行处理,进而产生水资源,从而实现水能与电能的协同产出。
通常水电联产机组的燃料成本与输出的水功率与电功率相关,其燃料成本表示如下:
式中,αc、βc、yc、ζc、ζc、ξc——分别为水电联产机组c的燃料消耗特性系数;Pc,t—表示为第t时段水电联产机组c的发电功率;Wc,t—表示为第t时段水电联产机组c的产水率;Uc,t—表示为二进制变量,当水电联产机组c开启时为1,否则为0。
2.2研究锂电/铅酸/氢储能等快速响应机理数据驱动分析方法
以锂电、铅酸、氢储能等为代表的化学储能系统,其充放电过程就是储能系统与电网完成能量交互的过程,而且它们之间的能量是双向流动的。当储能系统充电时,风电系统和光伏系统发出的多余电能被储能系统从电网中吸收;当储能系统放电时,将多发的电能向电网释放。
对于化学储能系统,其充放电能力与当前剩余电量有关,而当前剩余电量不仅与上一时刻的剩余电量有关,还与储能系统的充放电功率、充放电效率、额定容量及自放电率有关。
充电过程:
放电过程:
式中:S(t)、S(t-1)表示第t时刻、t-1时刻蓄电池的剩余电量;σ表示蓄电池的自放电率;ηc、ηd表示蓄电池的充电效率和放电效率;Pc、Pd表示蓄电池的充电功率和放电功率;E0表示储能系统的额定容量;Δt为时段t与时段t+1间的时间间隔。
除此之外,还可得到储能系统的能量公式:
E(t)=S(t-1)·E0(2.15)
根据储能系统对风电、光伏发电系统的控制作用的不同,储能系统主要分为以下四种控制模式:
(1)平滑功率输出模式
风、光的波动性和间歇性会导致其出力波动较大,给系统带来消极的作用。在该模式下,储能元件利用其秒级的充放电特性,能够平滑联合系统的总有功出力曲线,从而使得其并入电网时满足安全、可靠、稳定的运行要求。
(2)跟踪计划出力模式
在该模式下,监控系统会依照人们提前预测的出力曲线,严格调节和控制储能元件的充放电过程,使得系统发出的总有功出力曲线与计划出力曲线相似。计划出力曲线不仅可以依照系统所处区域的气候以及自然环境等预测得到,而且还可以通过负荷预测情况获得。
(3)负荷曲线削峰填谷模式
在该控制模式下,依照负荷峰谷特性的变化情况,储能元件也随之变换充放电方式。夜晚,系统负荷较低,储能元件处于充电状态,从电网中吸收并储存多余的风能和光能;白天,系统负荷较高,储能元件处于放电状态,向电网释放一定的能量,缓解负荷压力,使系统出力更加平稳。
(4)系统调频模式
在实际应用中,风力和光伏发电的有功功率随风速和太阳光照影响较大,其出力难以有效预测和控制,所以一般选用火电或水电机组来进行系统的调频。但是随着风光等新能源规模的不断扩大,用于电力系统调频的机组容量也不断降低,导致我们需要开发和研究新的调频方式。然而,储能系统具有秒级快速调节作用,可以快速地调节系统频率。
2.3研究锂电、铅酸、氢储能等多元储能数据驱动聚合建模技术
储能系统实现了同种能量在时间上的错位分布,能够调节供能在时间上的供需不匹配性,因此可以作为功能区综合能源系统的二次调节手段,弥补由于多能互补造成的能源供需不匹配性。
在区域综合能源系统运行时,以锂电、铅酸、氢储能为代表的电化学储能主要发挥平抑可再生能源出力不确定性、提高系统运行经济性的作用,可以在电网负荷低的时候储存电能,在电网负荷高的时候释放电能,用于削峰填谷,减轻电网波动,目前在可再生能源并网、分布式发电及微网领域已实现兆瓦级的示范应用。蓄热式电锅炉中的电锅炉是能量形式转换的耦合单元,一般利用光伏、风电等新能源作为电源,通过电磁感应或电阻方式发热,输出热水或高温蒸汽,加装储热罐可打破常规的“以热定电”运行模式,使得电锅炉的电负荷即具有可控性,在负荷低谷时段可增加电负荷,将多余的电量转化成热能储存起来;在负荷高峰期则降低电负荷需求,利用储存的热能进行供热。这样不仅使低谷时段的风电得到更有效的利用,而且起到削峰填谷的作用,极端情况下甚至可以将CHP机组停机,完全使用电锅炉供热,实现供热零排放。在各种冷储能方式中,冰蓄冷以将水制成冰的方式,进行冷量的储存。与水蓄冷相比,储存同样多的冷量,冰蓄冷所需的体积将比水蓄冷所需的体积小得多。
本课题假设采用化学储能、储热罐和蓄冷罐作为存储电热冷的实体设备,编号x分别为1~3,三种实体储能模型具有相似的运行特征,如充放过程存在损耗等,时段t功能区r实体储能x的存储能量Et,r,x与充/放功率间的关系为:
式中:kx为实体储能x的自损系数;Δt为调度时间间隔;和/>为实体储能x的充、放能效率。
3基于数据驱动的分布式能源电力交易和多元储能混合运行特性建模
研究基于数据驱动的微气象分布特征分析方法;研究基于数据驱动的区域内分散光伏特征聚类技术;研究计及需求响应策略的分布式能源优化运行方法;研究区域范围内分布式能源和多元储能混合运行特性建模分析方法。
3.1研究基于数据驱动的微气象分布特征分析方法
(1)山地微地形对气象参数的影响:
山地微地形对温度的影响
如果地面植被相同,在辐射型天气下,不同地形上的温度差异主要由辐射的差异引起,接受太阳辐射多的地形温度高,接受太阳辐射少的地形温度低。因此,温度分布随地形而变化的特点和规律基本与地形的辐射分布相类似的。同时,由于不同地形上的风速和湿度不同,也会在不同程度上加强或减弱这种由辐射分布特点所造成的温度差异。由于风速增大有使近地面层中温度降低的作用,冬季当有北方冷空气侵袭时,在迎风地形风速大,降温快。
在周围有山地环绕的山谷中,由于风速小和湍流交换弱,对白天增热和夜间冷却都造成有利条件,而且夜间周围山坡上的冷空气向谷底注泻,加强了谷中的冷却。在高地的顶部和坡地上,因为风速较大,湍流交换较强,且夜间地面附近的冷空气可以沿坡下拽,换来自由大气中较暖的空气,所以一般情况下,在夜间谷地的温度比高地顶部和坡地要低。
在自由大气中,热量的供应主要来源于地面的长波辐射,而大气直接吸收太阳的短波辐射相对较少,因此随着离地面高度增加,自由大气的温度一般是递减的。在中纬度地区,自由大气的年平均气温大约是每上升100m降低0.6℃。但在山区,由于地形、山体高度、干湿状况和盛行风向等影响,气温随高度的变化远比在自由大气中复杂得多。在冬季,华北地区平均气温的直减率为1.66℃。
山地微地形对湿度的影响
空气湿度的大小主要是决定于地面保水情况以及消耗水分蒸发的强弱,在微地形中,坡向对空气湿度的影响与天气和气候条件有关。一般是向风坡因降水和云雾多,空气和土壤湿度大,背风坡因降水和云雾少,空气和土壤湿度小。在比较湿润的气候条件下,因为土壤湿润,坡地上的蒸发主要决定于热力条件,所以接受太阳辐射多的坡地,随着被蒸发到空气中的水汽多,空气湿度大。接受太阳辐射少的坡地,因蒸发到空气中的水汽少,空气湿度低。在干燥的气候条件下,由于接受太阳辐射多的坡地土壤干燥,其蒸发能力虽强,但实际蒸发到空气中去的水汽少,空气湿度就反比接受太阳辐射少而土壤湿度较大的坡地低。
此外,微地形区域相对湿度的大小也决定于空气温度和水汽压的变化,并在很大程度上与云雾分布的特点有关。一般说来,凡是云雾多的地方,其相对湿度大,同时,在不同坡地上和不同高度上,由于风和温度分布的差异,也会引起空气湿度的不同。云雾随风场的运动在山谷中流动,并沿山体坡面抬升,气温按递减率下降,直到水汽达到饱和。
山地微地形对风速的影响
在一些微地形区域,中小地形能改变气流的速度和方向,同时因为各种因素所引起的增热差异或温度差异,会造成局部气压的差异,产生各种各样的局地环流,所以微地形区域的风场状况比温度和湿度的分布复杂得多。在两端开口的平直谷地内,风速的变化主要决定于自然风向与谷地走向的夹角,当风近乎平行谷地吹时,由于产生狭管效应,谷内风速比平地加强。
(2)微气象预测模型建立:
如果因变量y与自变量x之间没有确定性的函数关系,而根据知识、经验或观察,它们有一定的关联性,又有众多的、不可预测的随机因素影响着它们之间的关系,那么可以建立回归模型。由前面分析可知,影响微地形气象因子的因素有很多,而且气象因子温度、湿度、风速之间也相互影响。
因此可先将3个气象因子当中的1个因子作为因变量,其他2个量作为自变量,由此建立3者之间的回归模型为
Yi=β1+β2X2i+…+βjXki+ui(3.1)
式中,βj为偏回归系数。
多元回归分析的核心方法是利用最小二乘法求出偏回归系数。
样本回归函数为
式中,e为随机误差项。
寻求剩余平方和最小的参数估计式,即
令偏导数为0,则:
将式(3.2)两边同乘X′:
根据式(3.6)得出X′e=0,则偏回归系数的计算式为
度量回归方程整体拟合优度的统计量是可决系数R2,可以表达因变量与所有自变量之间的总体关系,即
3.2研究基于数据驱动的区域内分散光伏特征聚类技术
本课题设计了一种适用于负荷数据深度挖掘的二次聚类法。在聚类过程中,使用系统聚类法对区域内分散光伏特征进行一次聚类,并根据一次聚类得到的初始聚类中心,使用k-means聚类法对其进行二次聚类。通过设计的二次聚类法,不仅可以防止聚类效果受到初始聚类中心选择好坏的影响,而且能够使得聚类过程变得更加快速,提高分类效率。
此外,由于受到信号干扰、软件故障、设备性能等情况的影响,数据经常会出现未全面采集或者失真的现象。因此,为了避免其对后面的数据挖掘工作产生影响,需要先对获得的数据进行预处理。本课题研究主要采用不良数据识别和数据归一化的方法来对数据进行预处理。
系统聚类法:
在聚类分析中,我们假定某一类G中含有n个元素,可以通过列向量xi(i=1,2...,h)来表示,并用dij表示元素xi与xj之间的距离,用DKL表示类GK与类GK之间的距离。由于类与类之间的距离可以通过不同的方法来进行定义,因此就出现了不一样的系统聚类方法。其中常见的系统聚类法一共有六种,分别是最短距离法、最长距离法、中间距离法、重心法、类平均法和离差平方和法。在这六种系统聚类法中,离差平方和法的聚类效果最好。
离差平方和法又称为Ward法,它是把方差分析的思想用于分类上,同一个类内的离差平方和小,而类间的离差平方和应当大。类中各元素到类重心(即类均值)的平方欧氏距离之和称为类内离差平方和。设某一步GK和GL聚成一个新类GM,则GK、GL和GM的类内离差平方和分别如式(3.9)、(3.10)和(3.11)所示:
类内离差平方和反映了类内元素的分散程度。将GK与GL合并成新类GM时,类内离差平方和会有所增加,即WM-(WK+WL)>0。若GK与GL距离比较近,则增加的离差平方和应较小,于是定义GK与GL平方距离如式(3.12)所示:
类间平方距离的递推公式如式(3.13)所示:
系统聚类法的优点是原理简便,不需要进行初始参数设定。不过,在处理大数据样本时,计算量比较大,导致其计算速度比较慢。
k-means聚类:
k-means聚类给定D维空间内的n个数据点以及整数k,聚类问题为确定相同空间内(称为中心)k点的集合,使得其他数据点距离该中心的距离最小,步骤如下:
a.筛选数据点、聚类数和最大迭代数目;
b.选择k数据点作为初始聚类中心;
c.计算数据点与中心的距离,将数据按距离进行归类。利用欧氏距离进行判断:
d.计算每类的平均距离:
其中,ajm为m维数据在j类的均值;nj为j类的数据点数目;x为聚类中的数据点。
e.聚类的均值作为新的聚类中心,重复步骤3,直到达到最大迭代次数后类保持不变。
3.3研究计及需求响应策略的分布式能源优化运行方法
目前主流的需求响应策略分为基于价格的需求响应和基于激励的需求响应。基于价格的需求响应决定性因素在于需求弹性,需求响应的核心思想以实时电力供需关系为依据来调整电力价格,同时让用户以自身用电需求为前提,根据供电方所提供的实时电价来制定最优生产成本的用电模式,通常以激励电价或惩罚电价来表明供电方给出的电力价格是相对下降还是上涨。以现有的需求响应策略应用实践结果来看,当用户接受这样的用电策略时,他们的固有用电习惯将发生改变。本方法建立的计及需求响应策略的分布式能源优化运行模型如下:
能源服务商以净利润最大化为优化目标,优化运行的目标函数可以描述为:
式中T为优化周期;为t时刻能源服务商经济收益(元);/>为t时刻能源服务商运行成本(元);/>为t时刻能源服务商需求响应成本(元);/>为t时刻能源服务商环境成本(元)。
在日中,能源服务商需要下发激励来引导居民进行需求响应,其需求响应成本由下式所示:
式中NU为总用户数量(户);Ii,t为t时刻i下发用户的激励价格(元/kWh)。 为i用户在t时刻基础负荷、电热水器负荷、空调负荷、电冰箱负荷、洗衣机负荷、新能源汽车负荷的削减量。
约束条件:
式中为分布式光伏设备出力上限;/>为分布式风电设备出力上限;
3.4研究区域范围内分布式能源和多元储能混合运行特性建模分析方法
风光储系统主要包括能量管理系统、风力发电系统、光伏发电系统、储能系统、变流器、并网逆变器、直流母线和负荷。如图14所示为典型的风光储联合发电系统的结构示意图。组成部分的功能如下:
能量管理系统:对整个系统起控制作用,它依照风速情况和光照强度情况对风能和光伏发电的输出功率进行预测,与预测的负荷数据相结合,并调节储能系统的充放电,从而确保电力系统稳定运行。
风力发电系统:在风速较大的情况下,作为整个系统的主要能量来源,为电网提供一定的能源支持。
光伏发电系统:在光照较充足的情况下,作为整个发电系统的主要能量来源,除了提供给负荷外,多余的电能储存在储能元件里。
储能系统:具有平抑功率波动的作用,可以改善系统的总有功输出特性。此外,当风光充足时,储存多余的能量;当风光不足时,向负荷提供电能。
变流器:将风力、光伏和储能装置的出力转换为直流电压,并汇集到直流母线上。
并网逆变器:将直流电压转换为交流电压,并且作为风力、光伏和储能系统与电网相连的主要枢纽通道。
风光储联合发电系统的工作原理为:首先由能量管理系统根据当前风速、光照强度和负荷情况对风能和光伏发电的出力进行预测,对储能系统的充放电进行控制。然后风力、光伏和储能发电单元通过变流器将能量汇集到直流母线上,再经过并网逆变器将直流电压转换成与电网匹配的交流电压,最后经过升压变压器注入到电力网络中。
风光储联合发电系统的数学模型是以风光功率预测为基础,以风、光、储总有功功率输出与系统中发电机的有功功率输出之和等于负荷所需功率为约束,以稳定有功出力为目标,对储能单元按照一定的控制策略进行充放电。因此,建立风光储系统数学模型的关键就是控制好储能单元的充放电。其中,储能单元的充放电策略为:当风力、光伏与系统中发电机输出的有功功率之和小于负荷所需的功率时,有功功率缺额由储能系统放电补充;反之,当风力、光伏与系统中发电机输出的有功功率之和大于负荷所需的功率时,根据节约能量的原则,储能系统此时处于充电状态,将剩余的能量储存起来。
风光储系统并入电网后对系统的影响:
1)对系统有功功率波动的影响
风光储联合发电利用风光的互补特性以及储能平波的功能,当风、光、储容量配置合理的情况下,可以缓解系统的功率波动。但是在穿透率较高的情况下,虽然水电机组或燃气轮机等快速响应机组可以缓解可再生能源造成的有功功率波动,但是由于水电受资源分布等地域限制较大,而燃气轮机仍会使用化石能源,这将使得系统的经济性和环境效益受到了一定的影响。
2)对系统可靠性的影响
根据风能和太阳能间隙性和随机性的特点,加上由于目前功率预测水平还不是十分成熟,使得大量的新能源接入电网时将会对系统的可靠性造成严重的影响,而这一影响也是多方面的。
首先,由于风、光出力波动性和间隙性的特点,为了尽可能地减小负荷波动和潜在的停电事故,我们可以增大系统备用容量。其次,风光储联合发电系统中使用了大量的逆变器和变流器,这使得电网的电能质量也受到了影响。再次,输电线路的传输效率也受到消极的影响,其传输的电能还达不到需要的一半。
3)对电能质量的影响
风、光并网发电受气候因素的影响较为严重,随机波动性较大,导致系统的潮流分布发生变化,系统各点的电压也难以控制,对电能质量的影响较为严重,其影响主要包括以下几个方面:
a.电网频率
系统并网时对电网频率影响主要是由于风能、太阳能间隙性的特点。当电网容量比较大时,风光储系统并网产生的功率波动比较小,此时电网可以自身调节频率,因此也不会对电网频率造成较大的影响;而当风光储系统大量接入电网时,其所占的发电比重越来越大,产生的功率波动也较大。当供电电源与系统间的负荷平衡受到影响时,电网频率就会改变并达到一个新的稳定状态。
b.谐波污染
系统并网时会应用大量的电力电子装置来进行能量传递等功能,然而在进行操控时会造成-一系列的电网电压和电流畸变问题,这将导致谐波污染。除此之外,当这些电力电子装置的开关操作过于频繁时,开关频率附近将会产生谐波分量,这也会导致谐波污染。对于储能系统,在其并网充放电时,双向变流器的死区特性以及管压降等非线性特性也会造成谐波污染,从而影响电网的电能质量。
本发明所提方法应用于河北省张家口某新能源示范区,当地全年风力、光伏资源充足,主要分布与夏季。根据量测数据得到用户侧风/光新能源的聚合预测模型,预测精度与真实值如图4所示。
将分布式新能源与多元储能混合运行,并进行动态的能量输入/输出调节,得到多元储能系统的特性模型结果如图5(a)所示。
为提高系统能量利用率,本发明测试了多元储能系统数据模型的能效特性。搜索算法采用穷举法,以600W微型燃料电池为基础,设置分辨率为1W,极限容量为3000W,搜索目标为使燃料电池-尾气回收燃气轮机的极限功率最大。搜索结果如图6所示。从图6可看出,在燃料电池与尾气回收燃气轮机功率配比为600:1562时,多元储能的能量释放获得最大出力。由于燃料电池为了保证尾气回收燃气轮机的输出功率而在一定程度上提高电流密度增加产热量,且燃气轮机负荷率对其效率的影响,因此在最大功率点时燃料电池与尾气回收燃气轮机都并非满负荷运行,总能量利用率约为80%。考虑到燃料电池的温度限定,燃料电池、尾气回收燃气轮机与燃料电池-尾气回收燃气轮机在可行域内的运行情况分别如图5(b)、(c)、(d)所示。
为了验证所提方法并网运行的合理性,以IEEE-33节点系统为算例进行仿真分析。IEEE-33系统运行状态如图7所示。可以看到系统内各线路区段线损在理论值附近存在一定的波动幅度,验证了本文所提方法对于网络运行的安全性和经济性有可靠的保证。可以看出实际线损最大、最小值分别可比理论值高73.2%,低78.6%。其中,13号支路的精度欠佳;然而线损功率最高的2号和5号支路,受模型精度的影响并不大。这主要是因为13号支路后续节点负荷功率较小,线路传输电流不大,因此微小的数据驱动模型变化造成电流高比例波动,从而导致线损波动幅度较大;而2号、5号支路位于线路首段附近,传输功率大,数据驱动模型精准程度对其电流影响幅度较小,因此,网络运行真实状态与计算结果差异极小。
Claims (10)
1.一种考虑电力交易的分布式新能源和多元储能聚合建模方法,其特征在于:包括以下步骤:
步骤1、整合并分析现有配电网、用户侧系统中获得的数据,研究风电、光伏发电的聚合预测,构建基于数据驱动的用户侧风/光新能源电力交易的聚合预测模型;
步骤2、考虑多类型能源转换设备的非线性特性,对锂电、铅酸、氢储能储能的特性进行分析,建立锂电、铅酸、氢储能多元储能数据驱动聚合模型;
步骤3、分析本地区微气象分布特征和光伏分散特征,构建基于数据驱动的分布式能源的多元储能混合运行模型。
2.根据权利要求1所述的考虑电力交易的分布式新能源和多元储能聚合建模方法,其特征在于:所述步骤1的具体步骤包括:
步骤1.1:通过对移动互联网、智慧电网、用户侧系统获得的数据进行整合和提炼,获得电网拓扑结构、潮流计算、负荷预测数据;
步骤1.2:分析风电、光伏发电特性,研究风电、光伏的聚合预测;
步骤1.3:建立分布式光伏发电的电力交易模型,分析分布式电源参与电力交易的优势;
步骤1.4:运用数据驱动方法构建基于人工神经网络的用户侧风/光新能源的聚合预测模型。
3.根据权利要求2所述的考虑电力交易的分布式新能源和多元储能聚合建模方法,其特征在于:聚合预测方法包括以下步骤:
首先提取历史数值天气中的气象因素记录,并将记录中的各个因素数据进行标准化处理;再以预测时间点的气象为核心对标准化历史气象记录进行相似聚合,筛选原始样本;然后对相似聚合中记录的气象因素进行降维,得到独立因素;最后建立独立因素和风电/光伏功率的映射关系,实现聚合预测。
4.根据权利要求3所述的考虑电力交易的分布式新能源和多元储能聚合建模方法,其特征在于:气象因素包括:包括风速、风向、光照强度、气压、温度和湿度气象因素。
5.根据权利要求2所述的考虑电力交易的分布式新能源和多元储能聚合建模方法,其特征在于:所述步骤1.2包括以下步骤:
步骤1.2.1:提取历史数值天气中的气象因素记录,并将记录中的各个因素数据进行标准化处理,剔除异常点并将数据线性映射到[0,1]区间内;
步骤1.2.2:以预测时间点的气象为核心对标准化历史气象记录进行相似聚合;
步骤1.2.3:对相似聚合中记录的气象因素运用主成分分析法进行降维,得到独立因素;
步骤1.2.4:建立独立因素和风电/光伏功率的映射关系,实现聚合预测。
6.根据权利要求1所述的考虑电力交易的分布式新能源和多元储能聚合建模方法,其特征在于:所述步骤2的具体步骤包括:
步骤2.1:考虑多类型能源转换设备的非线性特性,对锂电、铅酸、氢储能等储能技术的启停响应速度、功率密度、能量转化效率以及成本特性进行分析总结;
步骤2.2:利用量测装置获得锂电、铅酸、氢储能等储能设备的历史能量流动数据,采用神经网络模型拟合上述储能设备的历史数据集,得到初步的储能设备非线性黑盒模型;
步骤2.3:根据具体选用的锂电池、铅酸电池、氢储能设备,通过对步骤2.2得到的初步的储能设备非线性黑盒模型进行最大荷电量额定值、极限功率等参数的调整,建立相应的储能模型。
7.根据权利要求6所述的一种考虑电力交易的分布式新能源和多元储能聚合建模方法,其特征在于:神经网络算法是由输入层、隐含层以及输出层组成,输入层使用n个神经元,输入信号使用向量X=[x1,x2,...,xn]T;隐含层包含n个神经元节点,隐含层输出向量为G=[g1,g2,...,gn]T;输出层中包含的向量为D=[d1,d2,...,dn]T。
8.根据权利要求1所述的考虑电力交易的分布式新能源和多元储能聚合建模方法,其特征在于:所述步骤3的具体步骤包括:
步骤3.1:掌握本地区的微气象分布特征和该区域内分散光伏特征;
步骤3.2:对数据进行建模,并利用K-中心点方法对微气象分布特征和光伏分散特征进行聚类;
步骤3.3:基于多元储能设备平抑风光电源出力的能量调度原则,将构建的分布式能源与多元储能的数据驱动模型进行串联,得到当地风光资源分布特性下,风光出力预测与多元储能运行特性的动态模型,构建基于数据驱动的分布式能源和多元储能混合运行模型。
9.根据权利要求8所述的考虑电力交易的分布式新能源和多元储能聚合建模方法,其特征在于:掌握本地区的微气象分布特征和该区域内分散光伏特征的途径包括文献检索以及通过气象局或电网公司获得。
10.根据权利要求8所述的考虑电力交易的分布式新能源和多元储能聚合建模方法,其特征在于:所述K-中心点方法为:
在n个数据对象中随机选择k个中心,集合中其他数据对象被划分到与中心最近的数据集合中去,形成了k个初始的数据集合;对于每个数据集合,计算新的中心,然后根据新的中心重新分配其他数据集合,这样不断迭代,直到所有的数据中心不再变化;该算法的具体步骤:
(1)随机选择k个初始聚类中心C={c1,c2…ck};
(2)如果s与聚类中心ci之间的距离小于该对象到cj的距离,则将s划分到聚类簇Ci中;
(3)重新计算聚类簇Ci的中心:/>
重复步骤(2)和(3)直到聚类中心C不再变化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310614453.2A CN117154808A (zh) | 2023-08-31 | 2023-08-31 | 考虑电力交易的分布式新能源和多元储能聚合建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310614453.2A CN117154808A (zh) | 2023-08-31 | 2023-08-31 | 考虑电力交易的分布式新能源和多元储能聚合建模方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117154808A true CN117154808A (zh) | 2023-12-01 |
Family
ID=88910724
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310614453.2A Pending CN117154808A (zh) | 2023-08-31 | 2023-08-31 | 考虑电力交易的分布式新能源和多元储能聚合建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117154808A (zh) |
-
2023
- 2023-08-31 CN CN202310614453.2A patent/CN117154808A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108599206B (zh) | 高比例不确定电源场景下的配电网混合储能配置方法 | |
CN111445067B (zh) | 一种适用于高铁站综合能源系统的多目标规划方法 | |
CN112994115B (zh) | 一种基于wgan场景模拟和时序生产模拟的新能源容量配置方法 | |
CN111313480B (zh) | 一种基于多目标寻优的多能互补系统的设计方法 | |
CN106886833A (zh) | 一种适用于复杂约束条件的风力发电机选址选型优化方法 | |
CN111144655A (zh) | 一种分布式电源选址定容与配电网网架联合优化方法 | |
CN110994606A (zh) | 一种基于复杂适应系统理论的多能源电源容量配置方法 | |
CN113937825A (zh) | 一种基于E-C-Kmeans聚类和SOP优化的DG双层优化配置方法 | |
CN112072643A (zh) | 一种基于深度确定性梯度策略的光-蓄系统在线调度方法 | |
Kang | Generation cost assessment of an isolated power system with a fuzzy wind power generation model | |
CN112801343A (zh) | 一种计及多气象场景适应成本的储能系统容量规划方法 | |
CN115940284B (zh) | 一种考虑分时电价的新能源制氢系统的运行控制策略 | |
CN110119850B (zh) | 基于光热发电调节的储热量双阶段优化调度方法 | |
CN114243766B (zh) | 一种区域多能源系统优化配置方法及系统 | |
CN117154808A (zh) | 考虑电力交易的分布式新能源和多元储能聚合建模方法 | |
Boqtob et al. | Energy scheduling of isolated microgrid with battery degradation cost using hybrid particle swarm optimization with sine cosine acceleration coefficients | |
Wan et al. | Multi-energy Complementarity Evaluation and Its Interaction with Wind and Photovoltaic Capacity optimization | |
Jinglin et al. | Multi-objective Collaborative Planning Method for Micro-energy Systems Considering Thermoelectric Coupling Clusters | |
Yang et al. | Locational marginal pricing-oriented transmission network expansion planning with consideration of line capacity sensitivity | |
Rahman et al. | Performance Evaluation of Probabilistic Clustering Techniques for Aggregating Wind Generators in Power System Dynamic Studies | |
CN112103989B (zh) | 一种考虑多类型负荷调控能力的新能源消纳贡献量化方法 | |
Xing et al. | An Optimization Capacity Design Method of Wind/Photovoltaic/Hydrogen Storage Power System Based on PSO-NSGA-II | |
Yin et al. | Optimization Method of power grid data fusion for High Proportion New Energy System | |
Li et al. | Low-Carbon Flexible Planning of Power Grid Using Alternative Energy Storage | |
Tang et al. | Design Hydro-Solar-Wind Multi-energy Complementary System via Multi-Objective Optimization |
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 |