一种提高风功率预测准确率的储能系统容量配置方法
一种提高风功率预测准确率的储能系统容量配置方法
【技术领域】
[0001] 本发明涉及实现提高风功率预测准确率的储能系统的容量配置问题领域，特别涉 及一种提高风功率预测准确率的储能系统容量配置方法。
【背景技术】
[0002] 风功率预测准确率是电网对风电场功率预测能力考核的重要指标之一。电网根据 风电场的功率预测准确率确定风电场的计划发电量，直接影响风电场的经济效益。按照目 前允许风电场并网发电的考核要求，风电场功率预测准确率提高的程度将决定风电场并网 发电的排名，可以不同程度的增加风电场的发电计划。所以确定不同风功率预测准确率提 高指标所需的储能容量对风电场有着重要的经济意义。
[0003] 根据我国《风电场功率预测预报管理暂行办法》规定，日预报要求并网风电场每日 在规定时间前按规定要求向电网调度机构提交次日〇时到24时每15分钟共96个时间节 点风电有功功率预测数据和装机容量。
[0004] 风电场发电日预测预报准确率为公式（1)。
[0006] 式中，ri为预测计划曲线准确率，N为日考核总时段数（取96点减去免考核点数），Pmk为k时段实际平均功率，Ppk为k时段的预测平均功率，Cap为风电场装机容量。
[0007] 利用储能提高风功率预测准确率已经有了一些研宄。考虑傅里叶变换法有在一个 分析周期内除直流分量，每一个频率的正弦分量都是完整的周波量的特点。若储能按照傅 里叶变换法得到的正弦分量出力，不考虑储能充放电过程中的能量损失，一个分析周期里 储能的充放电电量相当，储能净充（放）电电量为零。一个分析周期结束，储能能量状态可 以回归初始时刻。本专利将重点放在傅里叶变换法的使用，解决了储能装置一段时间后要 矫正能量状态回归初始状态的问题。通过合理组合傅里叶变换得到的幅频分量作为储能出 力，计算不同风功率预测准确率提高指标下所需的储能容量。
【发明内容】
[0008] 本发明的目的是提供一种提高风功率预测准确率的储能系统容量配置方法，利用 傅里叶变换法处理风电场历史数据，计算得到风电场不同功率预测准确率提高指标所需的 储能容量配置。
[0009] 本发明解决其技术问题所采用的技术方案是：
[0010] 一种提高风功率预测准确率的储能系统容量配置方法，利用傅里叶变换法分析预 测误差功率，在傅里叶变换的结果上用试差法计算得到不同风功率预测准确率提高指标所 需的储能容量，具体包括以下几个步骤：
[0011] 步骤1:傅里叶变换法分析预测误差功率：
[0012] (1)确定时间分析尺度，根据我国《风电场功率预测预报管理暂行办法》的相关规 定及傅里叶变换法的特点，选取合理的时间分析尺度对风电场历史数据进行分析；
[0013] (2)确定储能出力分量，根据傅里叶变换法的特点及对储能出力分量的要求，确定 可作为储能出力分量的傅里叶变换结果；
[0014] 步骤2 :试差法确定储能出力，用试差法从傅里叶变换法得到的若干幅频分量中 选取合理的满足风功率预测准确率提高指标的分量组合方式，在时域叠加这些分量可以得 到分析周期内的储能输出功率曲线；
[0015] 步骤3 :计算储能容量：
[0016] (1)根据储能的输出功率，计算得到储能的能量波动状态；
[0017] (2)根据储能的能量波动状态计算所需的储能容量，这里不考虑储能的充放电效 率和储能电池的工作范围；
[0018] 步骤4:算法编写，通过编写程序，计算多个风功率预测准确率提高指标所需的储 能容量配置，得到两者的关系曲线，通过曲线可以得到任一风功率预测准确率提高指标所 需的储能容量配置。
[0019] 进一步地，选取合理的时间分析尺度为一天。
[0020] 进一步地，幅频分量从高频分量开始选择，逐渐向低频段延伸。
[0021] 本发明和现有技术相比，具有以下优点和效果：本发明通过合理组合傅里叶变换 得到的幅频分量作为储能出力，可以得到风电场不同风功率预测准确率提高指标所需的储 能容量配置，并且可以不考虑储能在长时间运行过程中的能量状态矫正问题。电网根据风 电场的功率预测准确率确定风电场的计划发电量，大大提高了风电场的经济效益。
【附图说明】
[0022] 图1是储能容量计算流程图；
[0023] 图2是风电场样本数据及预测误差功率曲线；
[0024] 图3是风电场预测误差频谱分析结果；
[0025] 图4是储能输出功率曲线；
[0026] 图5是储能能量状态曲线；
[0027] 图6是储能容量计算结果；
[0028] 图7是储能最大输出功率计算结果。
【具体实施方式】
[0029] 下面结合附图并通过实施例对本发明作进一步的详细说明，以下实施例是对本发 明的解释而本发明并不局限于以下实施例。
[0030] 一种提高风功率预测准确率的储能系统容量配置方法，利用傅里叶变换法分析预 测误差功率，在傅里叶变换的结果上用试差法计算得到不同风功率预测准确率提高指标所 需的储能容量。采用风电场实际历史数据进行计算，得到风功率预测准确率提高量与储能 容量的关系曲线，最终得到提高不同风功率预测准确率提高指标所需储能容量的大小。
[0031] 1、傅里叶变换法分析误差功率
[0032] (1)离散傅里叶变换可以将一个采样得到的数字信号变换到频域。对功率信号P =[P[1]，. . .，P[n]，. . .，P[N]]进行离散傅里叶变换可以获得N个对应不同频率f的幅值M结果，如公式（2)所示。
[0034] 式中，P[n]为待处理的第n个数据；N为数据总数;M[n]是傅里叶变换结果中第n 个频率f[n]对应的幅值。f[n]的计算公式如公式（3)。
[0036] 式中，Ts是数据分析周期。
[0037] 根据《风电场功率预测预报管理暂行办法》对风功率预测准确率计算方法的规定， 考虑傅里叶变换所得幅频分量个数与分析数据个数相关的特点，选取本发明的时间分析尺 度为一天。
[0038] (2)《风电场功率预测预报管理暂行办法》规定一日风电场需提交的风功率预测数 据为96个。由傅里叶变换原理可知，对一日风电场预测误差功率的历史数据进行傅里叶变 换可以得到96个幅频分量。由傅里叶变换的对称性 可知，只需考虑0~1/2TS频率范围的 48个幅频分量即可。本发明要求储能出力分量在分析周期内是一个完整的周波分量，故直 流分量不能作为储能出力分量。综上考虑，共有47个幅频分量可作为储能出力分量。
[0039] 2、试差法确定储能出力
[0040] 试差法是在两个未知量中，先设一个量为已知，代入方程求出另一未知量，然后 进行验算。若所得结果与所设不符，重新修正假设再行计算，直到验算正确为止。本情境 可以先设确定的储能补偿范围，代入求得储能补偿后的功率预测准确率，若与预测目标准 确率不符，重新修正储能补偿范围再进行计算，直到达到风功率预测目标准确率。
[0041] 本发明对傅里叶变换结果从高频分量开始选择，逐渐向低频段延伸。对选中的幅 频分量在时域求和作为储能出力，直到储能的补偿功率可以使风功率预测准确率提高目标 指标。通过试差法可以得到提高目标风功率预测准确率指标时储能各个时刻的输出功率。 储能输出功率P0的表达式如公式（4)。
[0042] P0 = [P01，P02，P03,…P095，P096] (4)
[0043] 3、计算储能容量
[0044] (1)基于确定的储能输出功率，对各个时刻的储能充放电电量进行累计，可获得不 同采样时刻储能相对于初始状态的能量波动。亦即该时刻之前全部采样周期内储能累计充 放电能量之和。k时刻储能能量状态E0k的计算方法如公式（5)。
[0046] 式中，At为风功率预测的数据采样时间间隔，《办法》规定为15分钟。
[0047] (2)根据储能在整个分析周期内的能量状态曲线，储能在分析周期内能量波动的 最大值、最小值之差即为不考虑充放电效率和储能电池工作范围时的储能容量大小。储能 容量S的计算方法如公式（6)。
[0048] S=max(E0) -min(E0) (6)
[0049] 4、算法编写
[0050] 本专利具体的操作流程，如图1所示。
[0051] 5、储能配置结果
[0052] 采用一个装机容量为49. 5WM的风电场的一日历史数据进行计算。该天风功率预 测准确率为86. 2%。
[0053] 首先需要将风电场一日历史数据导入到算法中，计算出风电场一日功率预测误差 曲线。用傅里叶变换法分析误差功率，得到可作为储能出力的47个不同频率的正弦分量。 用试差法计算满足风功率预测准确率提高指标的储能输出功率。从而可以最终计算得到提 高相应风功率预测准确率提高指标所需的储能容量。
[0054] 风电场一日历史数据及风功率预测误差功率如图2。对风功率预测误差进行傅里 叶变换，〇~1/2TS频率范围内的幅频特性曲线如图3。
[0055] 通过仿真计算得到风功率预测准确率提高0. 91 %，需要储能容量为5. 284MW.h， 即0. 1067pu;需要储能最大输出功率为6. 0255WM，即0. 1217pu。储能输出功率曲线如图4， 储能能量波动状态如图5。
[0056] 计算多个风功率预测准确率提高量所需的储能容量和储能最大功率，可以得到储 能容量与准确率提高量关系如图6,储能最大输出功率与准确率提高量关系如图7。由储能 容量，储能最大输出功率与准确率提高量关系曲线可知，不同风功率预测准确率提高指标 所需的储能容量配置情况如表1所示。
[0057] 表1.不同风功率预测准确率提高指标所需储能配置结果
[0059] 本说明书中所描述的以上内容仅仅是对本发明所作的举例说明。本发明所属技术 领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式 替代，只要不偏离本发明说明书的内容或者超越本权利要求书所定义的范围，均应属于本 发明的保护范围。
【主权项】
1. 一种提高风功率预测准确率的储能系统容量配置方法，其特征在于：利用傅里叶变 换法分析预测误差功率，在傅里叶变换的结果上用试差法计算得到不同风功率预测准确率 提高指标所需的储能容量，具体包括以下步骤： 步骤1:傅里叶变换法分析预测误差功率： (1) 确定时间分析尺度，根据我国《风电场功率预测预报管理暂行办法》的相关规定及 傅里叶变换法的特点，选取合理的时间分析尺度对风电场历史数据进行分析； (2) 确定储能出力分量，根据傅里叶变换法的特点及对储能出力分量的要求，确定可作 为储能出力分量的傅里叶变换结果； 步骤2 :试差法确定储能出力，用试差法从傅里叶变换法得到的若干幅频分量中选取 合理的满足风功率预测准确率提高指标的分量组合方式，在时域叠加这些分量可以得到分 析周期内的储能输出功率曲线； 步骤3 :计算储能容量： (1) 根据储能的输出功率，计算得到储能的能量波动状态； (2) 根据储能的能量波动状态计算所需的储能容量，这里不考虑储能的充放电效率和 储能电池的工作范围； 步骤4 :算法编写，通过编写程序，计算多个风功率预测准确率提高指标所需的储能容 量配置，得到两者的关系曲线，通过曲线可以得到任一风功率预测准确率提高指标所需的 储能容量配置。2. 根据权利要求1所述的一种提高风功率预测准确率的储能系统容量配置方法，其特 征在于：选取合理的时间分析尺度为一天。3. 根据权利要求1所述的一种提高风功率预测准确率的储能系统容量配置方法，其特 征在于：幅频分量从高频分量开始选择，逐渐向低频段延伸。
【专利摘要】本发明提供一种提高风功率预测准确率的储能系统容量配置方法，属于实现提高风功率预测准确率的储能系统的容量配置问题领域。利用傅里叶变换法分析预测误差功率，在傅里叶变换的结果上用试差法计算得到不同风功率预测准确率提高指标所需的储能容量，具体包括以下几个步骤：步骤1：傅里叶变换法分析预测误差功率；步骤2：试差法确定储能出力；步骤3：计算储能容量；步骤4：算法编写。本发明通过合理组合傅里叶变换得到的幅频分量作为储能出力，可以得到风电场不同风功率预测准确率提高指标所需的储能容量配置，并且可以不考虑储能在长时间运行过程中的能量状态矫正问题。
【IPC分类】G06F19/00
【公开号】CN104899460
【申请号】CN201510333672
【发明人】黄梅, 苏粟, 谢桦, 龚敏明
【申请人】北京亿利智慧能源科技有限公司
【公开日】2015年9月9日
【申请日】2015年6月16日
【技术领域】
[0001] 本发明涉及实现提高风功率预测准确率的储能系统的容量配置问题领域，特别涉 及一种提高风功率预测准确率的储能系统容量配置方法。
【背景技术】
[0002] 风功率预测准确率是电网对风电场功率预测能力考核的重要指标之一。电网根据 风电场的功率预测准确率确定风电场的计划发电量，直接影响风电场的经济效益。按照目 前允许风电场并网发电的考核要求，风电场功率预测准确率提高的程度将决定风电场并网 发电的排名，可以不同程度的增加风电场的发电计划。所以确定不同风功率预测准确率提 高指标所需的储能容量对风电场有着重要的经济意义。
[0003] 根据我国《风电场功率预测预报管理暂行办法》规定，日预报要求并网风电场每日 在规定时间前按规定要求向电网调度机构提交次日〇时到24时每15分钟共96个时间节 点风电有功功率预测数据和装机容量。
[0004] 风电场发电日预测预报准确率为公式（1)。
[0006] 式中，ri为预测计划曲线准确率，N为日考核总时段数（取96点减去免考核点数），Pmk为k时段实际平均功率，Ppk为k时段的预测平均功率，Cap为风电场装机容量。
[0007] 利用储能提高风功率预测准确率已经有了一些研宄。考虑傅里叶变换法有在一个 分析周期内除直流分量，每一个频率的正弦分量都是完整的周波量的特点。若储能按照傅 里叶变换法得到的正弦分量出力，不考虑储能充放电过程中的能量损失，一个分析周期里 储能的充放电电量相当，储能净充（放）电电量为零。一个分析周期结束，储能能量状态可 以回归初始时刻。本专利将重点放在傅里叶变换法的使用，解决了储能装置一段时间后要 矫正能量状态回归初始状态的问题。通过合理组合傅里叶变换得到的幅频分量作为储能出 力，计算不同风功率预测准确率提高指标下所需的储能容量。
【发明内容】
[0008] 本发明的目的是提供一种提高风功率预测准确率的储能系统容量配置方法，利用 傅里叶变换法处理风电场历史数据，计算得到风电场不同功率预测准确率提高指标所需的 储能容量配置。
[0009] 本发明解决其技术问题所采用的技术方案是：
[0010] 一种提高风功率预测准确率的储能系统容量配置方法，利用傅里叶变换法分析预 测误差功率，在傅里叶变换的结果上用试差法计算得到不同风功率预测准确率提高指标所 需的储能容量，具体包括以下几个步骤：
[0011] 步骤1:傅里叶变换法分析预测误差功率：
[0012] (1)确定时间分析尺度，根据我国《风电场功率预测预报管理暂行办法》的相关规 定及傅里叶变换法的特点，选取合理的时间分析尺度对风电场历史数据进行分析；
[0013] (2)确定储能出力分量，根据傅里叶变换法的特点及对储能出力分量的要求，确定 可作为储能出力分量的傅里叶变换结果；
[0014] 步骤2 :试差法确定储能出力，用试差法从傅里叶变换法得到的若干幅频分量中 选取合理的满足风功率预测准确率提高指标的分量组合方式，在时域叠加这些分量可以得 到分析周期内的储能输出功率曲线；
[0015] 步骤3 :计算储能容量：
[0016] (1)根据储能的输出功率，计算得到储能的能量波动状态；
[0017] (2)根据储能的能量波动状态计算所需的储能容量，这里不考虑储能的充放电效 率和储能电池的工作范围；
[0018] 步骤4:算法编写，通过编写程序，计算多个风功率预测准确率提高指标所需的储 能容量配置，得到两者的关系曲线，通过曲线可以得到任一风功率预测准确率提高指标所 需的储能容量配置。
[0019] 进一步地，选取合理的时间分析尺度为一天。
[0020] 进一步地，幅频分量从高频分量开始选择，逐渐向低频段延伸。
[0021] 本发明和现有技术相比，具有以下优点和效果：本发明通过合理组合傅里叶变换 得到的幅频分量作为储能出力，可以得到风电场不同风功率预测准确率提高指标所需的储 能容量配置，并且可以不考虑储能在长时间运行过程中的能量状态矫正问题。电网根据风 电场的功率预测准确率确定风电场的计划发电量，大大提高了风电场的经济效益。
【附图说明】
[0022] 图1是储能容量计算流程图；
[0023] 图2是风电场样本数据及预测误差功率曲线；
[0024] 图3是风电场预测误差频谱分析结果；
[0025] 图4是储能输出功率曲线；
[0026] 图5是储能能量状态曲线；
[0027] 图6是储能容量计算结果；
[0028] 图7是储能最大输出功率计算结果。
【具体实施方式】
[0029] 下面结合附图并通过实施例对本发明作进一步的详细说明，以下实施例是对本发 明的解释而本发明并不局限于以下实施例。
[0030] 一种提高风功率预测准确率的储能系统容量配置方法，利用傅里叶变换法分析预 测误差功率，在傅里叶变换的结果上用试差法计算得到不同风功率预测准确率提高指标所 需的储能容量。采用风电场实际历史数据进行计算，得到风功率预测准确率提高量与储能 容量的关系曲线，最终得到提高不同风功率预测准确率提高指标所需储能容量的大小。
[0031] 1、傅里叶变换法分析误差功率
[0032] (1)离散傅里叶变换可以将一个采样得到的数字信号变换到频域。对功率信号P =[P[1]，. . .，P[n]，. . .，P[N]]进行离散傅里叶变换可以获得N个对应不同频率f的幅值M结果，如公式（2)所示。
[0034] 式中，P[n]为待处理的第n个数据；N为数据总数;M[n]是傅里叶变换结果中第n 个频率f[n]对应的幅值。f[n]的计算公式如公式（3)。
[0036] 式中，Ts是数据分析周期。
[0037] 根据《风电场功率预测预报管理暂行办法》对风功率预测准确率计算方法的规定， 考虑傅里叶变换所得幅频分量个数与分析数据个数相关的特点，选取本发明的时间分析尺 度为一天。
[0038] (2)《风电场功率预测预报管理暂行办法》规定一日风电场需提交的风功率预测数 据为96个。由傅里叶变换原理可知，对一日风电场预测误差功率的历史数据进行傅里叶变 换可以得到96个幅频分量。由傅里叶变换的对称性 可知，只需考虑0~1/2TS频率范围的 48个幅频分量即可。本发明要求储能出力分量在分析周期内是一个完整的周波分量，故直 流分量不能作为储能出力分量。综上考虑，共有47个幅频分量可作为储能出力分量。
[0039] 2、试差法确定储能出力
[0040] 试差法是在两个未知量中，先设一个量为已知，代入方程求出另一未知量，然后 进行验算。若所得结果与所设不符，重新修正假设再行计算，直到验算正确为止。本情境 可以先设确定的储能补偿范围，代入求得储能补偿后的功率预测准确率，若与预测目标准 确率不符，重新修正储能补偿范围再进行计算，直到达到风功率预测目标准确率。
[0041] 本发明对傅里叶变换结果从高频分量开始选择，逐渐向低频段延伸。对选中的幅 频分量在时域求和作为储能出力，直到储能的补偿功率可以使风功率预测准确率提高目标 指标。通过试差法可以得到提高目标风功率预测准确率指标时储能各个时刻的输出功率。 储能输出功率P0的表达式如公式（4)。
[0042] P0 = [P01，P02，P03,…P095，P096] (4)
[0043] 3、计算储能容量
[0044] (1)基于确定的储能输出功率，对各个时刻的储能充放电电量进行累计，可获得不 同采样时刻储能相对于初始状态的能量波动。亦即该时刻之前全部采样周期内储能累计充 放电能量之和。k时刻储能能量状态E0k的计算方法如公式（5)。
[0046] 式中，At为风功率预测的数据采样时间间隔，《办法》规定为15分钟。
[0047] (2)根据储能在整个分析周期内的能量状态曲线，储能在分析周期内能量波动的 最大值、最小值之差即为不考虑充放电效率和储能电池工作范围时的储能容量大小。储能 容量S的计算方法如公式（6)。
[0048] S=max(E0) -min(E0) (6)
[0049] 4、算法编写
[0050] 本专利具体的操作流程，如图1所示。
[0051] 5、储能配置结果
[0052] 采用一个装机容量为49. 5WM的风电场的一日历史数据进行计算。该天风功率预 测准确率为86. 2%。
[0053] 首先需要将风电场一日历史数据导入到算法中，计算出风电场一日功率预测误差 曲线。用傅里叶变换法分析误差功率，得到可作为储能出力的47个不同频率的正弦分量。 用试差法计算满足风功率预测准确率提高指标的储能输出功率。从而可以最终计算得到提 高相应风功率预测准确率提高指标所需的储能容量。
[0054] 风电场一日历史数据及风功率预测误差功率如图2。对风功率预测误差进行傅里 叶变换，〇~1/2TS频率范围内的幅频特性曲线如图3。
[0055] 通过仿真计算得到风功率预测准确率提高0. 91 %，需要储能容量为5. 284MW.h， 即0. 1067pu;需要储能最大输出功率为6. 0255WM，即0. 1217pu。储能输出功率曲线如图4， 储能能量波动状态如图5。
[0056] 计算多个风功率预测准确率提高量所需的储能容量和储能最大功率，可以得到储 能容量与准确率提高量关系如图6,储能最大输出功率与准确率提高量关系如图7。由储能 容量，储能最大输出功率与准确率提高量关系曲线可知，不同风功率预测准确率提高指标 所需的储能容量配置情况如表1所示。
[0057] 表1.不同风功率预测准确率提高指标所需储能配置结果
[0059] 本说明书中所描述的以上内容仅仅是对本发明所作的举例说明。本发明所属技术 领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式 替代，只要不偏离本发明说明书的内容或者超越本权利要求书所定义的范围，均应属于本 发明的保护范围。
【主权项】
1. 一种提高风功率预测准确率的储能系统容量配置方法，其特征在于：利用傅里叶变 换法分析预测误差功率，在傅里叶变换的结果上用试差法计算得到不同风功率预测准确率 提高指标所需的储能容量，具体包括以下步骤： 步骤1:傅里叶变换法分析预测误差功率： (1) 确定时间分析尺度，根据我国《风电场功率预测预报管理暂行办法》的相关规定及 傅里叶变换法的特点，选取合理的时间分析尺度对风电场历史数据进行分析； (2) 确定储能出力分量，根据傅里叶变换法的特点及对储能出力分量的要求，确定可作 为储能出力分量的傅里叶变换结果； 步骤2 :试差法确定储能出力，用试差法从傅里叶变换法得到的若干幅频分量中选取 合理的满足风功率预测准确率提高指标的分量组合方式，在时域叠加这些分量可以得到分 析周期内的储能输出功率曲线； 步骤3 :计算储能容量： (1) 根据储能的输出功率，计算得到储能的能量波动状态； (2) 根据储能的能量波动状态计算所需的储能容量，这里不考虑储能的充放电效率和 储能电池的工作范围； 步骤4 :算法编写，通过编写程序，计算多个风功率预测准确率提高指标所需的储能容 量配置，得到两者的关系曲线，通过曲线可以得到任一风功率预测准确率提高指标所需的 储能容量配置。2. 根据权利要求1所述的一种提高风功率预测准确率的储能系统容量配置方法，其特 征在于：选取合理的时间分析尺度为一天。3. 根据权利要求1所述的一种提高风功率预测准确率的储能系统容量配置方法，其特 征在于：幅频分量从高频分量开始选择，逐渐向低频段延伸。
【专利摘要】本发明提供一种提高风功率预测准确率的储能系统容量配置方法，属于实现提高风功率预测准确率的储能系统的容量配置问题领域。利用傅里叶变换法分析预测误差功率，在傅里叶变换的结果上用试差法计算得到不同风功率预测准确率提高指标所需的储能容量，具体包括以下几个步骤：步骤1：傅里叶变换法分析预测误差功率；步骤2：试差法确定储能出力；步骤3：计算储能容量；步骤4：算法编写。本发明通过合理组合傅里叶变换得到的幅频分量作为储能出力，可以得到风电场不同风功率预测准确率提高指标所需的储能容量配置，并且可以不考虑储能在长时间运行过程中的能量状态矫正问题。
【IPC分类】G06F19/00
【公开号】CN104899460
【申请号】CN201510333672
【发明人】黄梅, 苏粟, 谢桦, 龚敏明
【申请人】北京亿利智慧能源科技有限公司
【公开日】2015年9月9日
【申请日】2015年6月16日
最新回复(0)