CN103927595A - 基于自学习复合数据源光伏发电功率超短期预测方法 - Google Patents
基于自学习复合数据源光伏发电功率超短期预测方法 Download PDFInfo
- Publication number
- CN103927595A CN103927595A CN201410163054.XA CN201410163054A CN103927595A CN 103927595 A CN103927595 A CN 103927595A CN 201410163054 A CN201410163054 A CN 201410163054A CN 103927595 A CN103927595 A CN 103927595A
- Authority
- CN
- China
- Prior art keywords
- data
- model
- photovoltaic generation
- input
- generation power
- 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
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于自学***均模型参数;输入光伏发电功率预测所需输入数据到根据自回归滑动平均模型的参数确定的自回归滑动平均模型中得到预测结果,对预测结果进行后评估,即分析预测值与实测值之间的误差,如预测误差大于允许的最大误差,则从新进行模型定阶和模型参数估计。通过对光伏发电过程中的光伏发电功率进行预测,为新能源发电实时调度、新能源发电日前计划、新能源发电月度计划、新能源发电能力评估和弃光电量估计提供关键信息。通过引入复合数据源有效提高光伏发电功率超短期预测精度,从而实现在保障电网安全稳定经济运行的前提下有效提高新能源上网电量目的。
Description
技术领域
本发明涉及新能源发电过程中光伏发电功率预测技术领域,具体地,涉及一种基于自学***均模型的光伏发电功率超短期预测方法。
背景技术
我国风电进入规模化发展阶段以后所产生的大型新能源基地多数位于“三北地区”(西北、东北、华北),大型新能源基地一般远离负荷中心,其电力需要经过长距离、高电压输送到负荷中心进行消纳。由于风、光资源的间歇性、随机性和波动性,导致大规模新能源基地的风电、光伏发电出力会随之发生较大范围的波动,进一步导致输电网络充电功率的波动,给电网运行安全带来一系列问题。
截至2014年4月,光伏发电装机容量已达到435万千瓦,约占甘肃电网总装机容量的13%,同时甘肃成为我国光伏发电装机规模最大的省份。目前,甘肃电网风电、光伏发电装机超过甘肃电网总装机容量的1/3。随着新能源并网规模的不断提高,光伏发电不确定性和不可控性给电网的安全稳定经济运行带来诸多问题。准确预估可利用的发电光资源是对大规模光伏发电优化调度的基础。对光伏发电过程中的光伏发电功率进行预测,可为新能源发电实时调度、新能源发电日前计划、新能源发电月度计划、新能源发电能力评估和弃光电量估计提供关键信息。
发明内容
本发明的目的在于,针对上述问题,提出一种基于自学习复合数据源光伏发电功率超短期预测方法,以实现高精度光伏发电功率超短期预测的优点。
为实现上述目的,本发明采用的技术方案是:
一种基于自学***均模型参数;
输入光伏发电功率预测所需输入数据到根据上述自回归滑动平均模型的参数确定的自回归滑动平均模型中得到预测结果;
对预测结果进行后评估,即分析预测值与实测值之间的误差,如预测误差大于允许的最大误差,则从新进行模型定阶和模型参数估计;
所述输入数据得到自回归滑动平均模型参数包括,步骤101、输入模型训练基础数据;
步骤102、模型定阶;
步骤103、采用矩估计方法对定阶的ARMA(p,q)模型参数进行估计。
根据本发明的优选实施例,所述步骤101输入模型训练基础数据,输入数据包括,光伏电站基础信息、历史辐照度数据、历史功率数据和地理信息***数据。
根据本发明的优选实施例,所述步骤102模型定阶:
采用残差方差图法进行模型定阶,具体为设xt为需要估计的项,xt-1,xt-2,...,xt-n为已知历史功率序列,对于ARMA(p,q)模型,模型定阶即确定模型中参数p和q的值;
用系列阶数逐渐递增的模型拟合原始序列,每次都计算残差平方和然后画出阶数和的图形,当阶数由小增大时,会显著下降,达到真实阶数后的值会逐渐趋于平缓,甚至反而增大,
实际观测值个数指拟合模型时实际使用的观察值项数,对于具有N个观察值的序列,拟合AR(p)模型,则实际使用的观察值最多为N-p,模型参数个数指所建立的模型中实际包含的参数个数,对于含有均值的模型,模型参数个数为模型阶数加1,对于N个观测值的序列,ARMA模型的残差估计式为:
根据本发明的优选实施例,所述步骤103采用矩估计方法对定阶的ARMA(p,q)模型参数进行估计具体步骤为:
将光伏电站历史功率数据利用数据序列x1,x2,...,xt表示,其样本自协方差定义为
其中,k=0,1,2,...,n-1,xt和xt-k均为数据序列x1,x2,...,xt中的数值;
则
则历史功率数据样本自相关函数为:
其中,k=0,1,2,...,n-1。
AR部分的矩估计为,
令
则协方差函数为
用的估计代替γk,
可得参数
对MA(q)模型系数θ1,θ2,...,θq采用矩估计有
其中k=1,2,...,m,
以上m+1个方程非线性方程,采用迭代法进行求解即得到自回归滑动平均模型参数。
根据本发明的优选实施例,所述输入光伏发电功率预测所需输入数据到根据上述自回归滑动平均模型的参数确定的自回归滑动平均模型中得到预测结果的步骤包括,
步骤201、输入功率预测基础数据;
步骤202、对输入的基础数据进行噪声滤波及数据预处理;
步骤203、根据确定的参数建立自回归滑动平均模型,并将处理后的数据输入从而得到预测结果;
步骤204、将预测结果输出至数据库中,并通过图表及曲线展示预测结果、并展示预测与实测结果的对比。
根据本发明的优选实施例,所述输入功率预测基础数据包括资源监测***数据和运行监测***数据,所述资源监测***数据包含光资源监测数据;所述运行监测***数据包括光伏组件监测数据、升压站监测数据和数据采集与监视控制***数据。
根据本发明的优选实施例,所述噪声滤波及数据预处理具体为:噪声滤波模块对监测***实时采集得到的带有噪声的数据进行滤波处理,去除坏数据和奇异值;数据预处理模块对数据进行对齐、归一化处理和分类筛选处理。
根据本发明的优选实施例,所述自回归滑动平均模型为:
其中,和θj(1≤j≤q)是系数,αt是白噪声序列。
本发明的技术方案具有以下有益效果:
本发明的技术方案通过对光伏发电过程中的光伏发电功率进行预测,为新能源发电实时调度、新能源发电日前计划、新能源发电月度计划、新能源发电能力评估和弃光电量估计提供关键信息。通过引入复合数据源有效提高光伏发电功率超短期预测精度,从而实现在保障电网安全稳定经济运行的前提下有效提高新能源上网电量目的。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明实施例所述的基于自学习复合数据源光伏发电功率超短期预测方法的原理框图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
一种基于自学***均模型参数;
输入光伏发电功率预测所需输入数据到根据自回归滑动平均模型的参数确定的自回归滑动平均模型中得到预测结果;
对预测结果进行后评估,即分析预测值与实测值之间的误差,如预测误差大于允许的最大误差,则从新进行模型定阶和模型参数估计;
输入数据得到自回归滑动平均模型参数包括,步骤101、输入模型训练基础数据;
步骤102、模型定阶;
步骤103、采用矩估计方法对定阶的ARMA(p,q)模型参数进行估计。
含大规模光伏发电的电力***运行依赖庞大的、准确的数据集,而光伏发电功率预测若能将这些数据有效融合利用则可有效提高预测精度。与常规电力***SCADA监测不同,在各类电气、机械和热力等数据之外,光伏发电监测数据还包含大量的资源监测、运行监测及地理信息等。
如图1所示,本发明技术方案提出的光伏电功率超短期预测可分为两个阶段:模型训练阶段和功率预测阶段。
阶段1:模型训练
步骤1.1:模型训练基础数据输入
光伏发电功率预报***模型训练所需输入数据包括,光伏电站基础信息、历史辐照度数据、历史功率数据,地理信息***(GIS)数据(光伏电站坐标、测光站坐标、升压站坐标等)。将基础数据输入到预测模型中进行模型训练。
步骤1.2:模型定阶
由于事先无法确定需要使用多少已知时间序列的项来建立估计函数,所以需要对模型进行定阶判断。
设xt为需要估计的项,xt-1,xt-2,...,xt-n为已知历史功率序列,对于ARMA(p,q)模型,模型定阶就是确定模型中参数p和q的值。
采用残差方差图法进行模型定阶。假定模型是有限阶自回归模型,如果设置的阶数小于真实阶数,则是一种不足拟合,因而拟合残差平方和必定偏大,此时通过提高阶数可以显著降低残差平方和。反之,如果阶数已经达到真实值,那么再增加阶数,就是过度拟合,此时增加阶数不会令残差平方和显著减小,甚至会略有增加。
这样用系列阶数逐渐递增的模型来拟合原始序列,每次都计算残差平方和然后画出阶数和的图形。当阶数由小增大时,会显著下降,达到真实阶数后的值会逐渐趋于平缓,有时甚至反而增大。残差方差的估计式为:
“实际观测值个数”是指拟合模型时实际使用的观察值项数,对于具有N个观察值的序列,拟合AR(p)模型,则实际使用的观察值最多为N-p。
“模型参数个数”是指所建立的模型中实际包含的参数个数,对于含有均值的模型,模型参数个数为模型阶数加1。对于N个观测值的序列,相应ARMA模型的残差估计式为:
(式1)
其中,公式中,Q为拟合误差的平方和函数,和θj(1≤j≤q)是模型系数,N是观测序列长度,是模型参数中的常数项,的常识值是根据不同的和θj(1≤j≤q)进行变化的常数项,不同的和θj(1≤j≤q)对照不同的值。
步骤1.3:模型参数估计
采用矩估计方法对ARMA(p,q)的模型参数进行估计。首先,将光伏电站历史功率数据利用数据序列x1,x2,...,xt表示,其样本自协方差定义为
其中,k=0,1,2,...,n-1,xt和xt-k均为数据序列x1,x2,...,xt中的数值。
特别的,
则历史功率数据样本自相关函数为:
AR部分的矩估计为
(式5)令
(式6)
则协方差函数为
(式7)
用的估计代替γk,有
(式8)可得参数
对MA(q)模型系数θ1,θ2,...,θq采用矩估计有
……
……
……
其中k=1,2,...,m。
以上共包含m+1个方程,对其参数而言,方程为非线性,采用迭代法进行求解。
具体步骤如下,将方程变形为:
给定θ1,θ2,...,θq和的一组初始值,如
θ1=θ2=...=θq=0,(式13)
代入以上两式右边,左边所得到的值为第一步迭代值,记为
再将该值依次代入上两式的右侧,便得到第二步迭代值,
依次类推,直到相邻两次迭代结果小于给定阈值时,取所得的结果作为参数的近似解。
通过上述求解过程发现,要求解时间序列模型的阶数,就要得到时间序列的预测值;要得到时间序列的预测值,必须先建立具体的预测函数;要建立具体的预测函数,必须知道模型的阶数。
根据实践得知,时间序列模型阶数一般不超过5阶。所以在该算法具体实现时,可以首先假设模型为1阶,利用步骤1.3中的参数估计方法得到一阶模型的参数,进而建立估计函数便可以求得一阶模型时间序列模型估计得到各个项的预测值,从而求得一阶模型的残差方差;之后,假设模型为二阶,用上述方法求得二阶模型的残差;以此类推,可以得到1到5阶模型的残差,选残差最小的模型的阶数作为最终模型的阶数。确定模型阶数后,便可计算得到参数θ1,θ2,...,θq的值。
阶段2:功率预测
步骤2.1:功率预测基础数据输入
光伏发电功率预测所需输入数据包括资源监测***数据和运行监测***数据两部分,其中,资源监测***数据包含光资源监测数据;运行监测***数据包括光伏组件监测数据、升压站监测数据和数据采集与监视控制***(SCADA)数据等。
步骤2.2:噪声滤波及数据预处理
噪声滤波模块对实时监测***采集得到的带有噪声的进行滤波处理,去除坏数据和奇异值;数据预处理模块对数据进行对齐、归一化处理和分类筛选等操作,以便使得输入的数据可以为模型所用。
步骤2.3:超短期功率预测
将模型参数估计出来之后,结合已估计的模型阶数,便可得到用于光伏发电功率超短期预测的时间序列方程。根据上述步骤2和步骤3得出的p和q值,以及θ1,θ2,...,θq的值建立自回归滑动平均模型;
自回归滑动平均模型如下:
(式14)
其中,和θj(1≤j≤q)是系数,αt是白噪声序列。
步骤2.4:预测结果输出及展示
将预测结果输出至数据库中,并通过图表及曲线展示预测结果、展示预测与实测结果的对比。
步骤3:预测结果后评估及模型修正
首先对预测结果进行后评估,分析预测值与实测值之间的误差。如果预测误差大于允许的最大误差,则跳转到模型训练过程,从新进行模型定阶和模型参数估计。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种基于自学***均模型参数;
输入光伏发电功率预测所需输入数据到根据上述自回归滑动平均模型的参数确定的自回归滑动平均模型中得到预测结果;
对预测结果进行后评估,即分析预测值与实测值之间的误差,如预测误差大于允许的最大误差,则从新进行模型定阶和模型参数估计;
所述输入数据得到自回归滑动平均模型参数包括,步骤101、输入模型训练基础数据;
步骤102、模型定阶;
步骤103、采用矩估计方法对定阶的ARMA(p,q)模型参数进行估计。
2.根据权利要求1所述的基于自学习复合数据源光伏发电功率超短期预测方法,其特征在于,所述步骤101输入模型训练基础数据,输入数据包括,光伏电站基础信息、历史辐照度数据、历史功率数据和地理信息***数据。
3.根据权利要求2所述的基于自学习复合数据源光伏发电功率超短期预测方法,其特征在于,所述步骤102模型定阶:
采用残差方差图法进行模型定阶,具体为设xt为需要估计的项,xt-1,xt-2,...,xt-n为已知历史功率序列,对于ARMA(p,q)模型,模型定阶即确定模型中参数p和q的值;
用系列阶数逐渐递增的模型拟合原始序列,每次都计算残差平方和然后画出阶数和的图形,当阶数由小增大时,会显著下降,达到真实阶数后的值会逐渐趋于平缓,甚至反而增大,
实际观测值个数指拟合模型时实际使用的观察值项数,对于具有N个观察值的序列,拟合AR(p)模型,则实际使用的观察值最多为N-p,模型参数个数指所建立的模型中实际包含的参数个数,对于含有均值的模型,模型参数个数为模型阶数加1,对于N个观测值的序列,ARMA模型的残差估计式为:
其中,Q为拟合误差的平方和函数,和θj(1≤j≤q)是模型系数,N是观测序列长度,是模型参数中的常数项。
4.根据权利要求3所述的基于自学习复合数据源光伏发电功率超短期预测方法,其特征在于,所述步骤103采用矩估计方法对定阶的ARMA(p,q)模型参数进行估计具体步骤为:
将光伏电站历史功率数据利用数据序列x1,x2,...,xt表示,其样本自协方差定义为
其中,k=0,1,2,...,n-1,xt和xt-k均为数据序列x1,x2,...,xt中的数值;
则
则历史功率数据样本自相关函数为:
其中,k=0,1,2,...,n-1;
AR部分的矩估计为,
令
则协方差函数为
用的估计代替γk,
可得参数
对MA(q)模型系数θ1,θ2,...,θq采用矩估计有
其中k=1,2,...,m,
以上m+1个方程非线性方程,采用迭代法进行求解即得到自回归滑动平均模型参数。
5.根据权利要求4所述的基于自学***均模型的参数确定的自回归滑动平均模型中得到预测结果的步骤包括,
步骤201、输入功率预测基础数据;
步骤202、对输入的基础数据进行噪声滤波及数据预处理;
步骤203、根据确定的参数建立自回归滑动平均模型,并将处理后的数据输入从而得到预测结果。
6.根据权利要求5所述的基于自学习复合数据源光伏发电功率超短期预测方法,其特征在于,还包括,
步骤204、将预测结果输出至数据库中,并通过图表及曲线展示预测结果、并展示预测与实测结果的对比。
7.根据权利要求6所述的基于自学习复合数据源光伏发电功率超短期预测方法,其特征在于,所述输入功率预测基础数据包括资源监测***数据和运行监测***数据,所述资源监测***数据包含光资源监测数据;所述运行监测***数据包括光伏组件监测数据、升压站监测数据和数据采集与监视控制***数据。
8.根据权利要求6所述的基于自学习复合数据源光伏发电功率超短期预测方法,其特征在于,所述噪声滤波及数据预处理具体为:噪声滤波模块对监测***实时采集得到的带有噪声的数据进行滤波处理,去除坏数据和奇异值;数据预处理模块对数据进行对齐、归一化处理和分类筛选处理。
9.根据权利要求6所述的基于自学***均模型为:
其中,和θj(1≤j≤q)是系数,αt是白噪声序列。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410163054.XA CN103927595A (zh) | 2014-04-22 | 2014-04-22 | 基于自学习复合数据源光伏发电功率超短期预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410163054.XA CN103927595A (zh) | 2014-04-22 | 2014-04-22 | 基于自学习复合数据源光伏发电功率超短期预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103927595A true CN103927595A (zh) | 2014-07-16 |
Family
ID=51145811
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410163054.XA Pending CN103927595A (zh) | 2014-04-22 | 2014-04-22 | 基于自学习复合数据源光伏发电功率超短期预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103927595A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104850914A (zh) * | 2015-05-29 | 2015-08-19 | 东南大学 | 一种基于特征建模的新能源发电量预测方法 |
CN105184097A (zh) * | 2015-09-29 | 2015-12-23 | 国网甘肃省电力公司 | 基于网架结构分群的大型光伏发电集群弃光电量评估方法 |
CN113378459A (zh) * | 2021-06-02 | 2021-09-10 | 兰州交通大学 | 基于卫星和物联网信息的光伏电站超短期功率预测方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103473322A (zh) * | 2013-09-13 | 2013-12-25 | 国家电网公司 | 基于时间序列模型的光伏发电功率超短期预测方法 |
-
2014
- 2014-04-22 CN CN201410163054.XA patent/CN103927595A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103473322A (zh) * | 2013-09-13 | 2013-12-25 | 国家电网公司 | 基于时间序列模型的光伏发电功率超短期预测方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104850914A (zh) * | 2015-05-29 | 2015-08-19 | 东南大学 | 一种基于特征建模的新能源发电量预测方法 |
CN104850914B (zh) * | 2015-05-29 | 2018-07-17 | 东南大学 | 一种基于特征建模的新能源发电量预测方法 |
CN105184097A (zh) * | 2015-09-29 | 2015-12-23 | 国网甘肃省电力公司 | 基于网架结构分群的大型光伏发电集群弃光电量评估方法 |
CN113378459A (zh) * | 2021-06-02 | 2021-09-10 | 兰州交通大学 | 基于卫星和物联网信息的光伏电站超短期功率预测方法 |
CN113378459B (zh) * | 2021-06-02 | 2022-09-02 | 兰州交通大学 | 基于卫星和物联网信息的光伏电站超短期功率预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103927695A (zh) | 基于自学习复合数据源的风电功率超短期预测方法 | |
Zhang et al. | Modeling conditional forecast error for wind power in generation scheduling | |
CN103996073A (zh) | 测光网络实时校正自学习arma模型光伏功率预测方法 | |
WO2017054537A1 (zh) | 一种长时间尺度光伏出力时间序列建模方法及装置 | |
Tabone et al. | Modeling variability and uncertainty of photovoltaic generation: A hidden state spatial statistical approach | |
CN103559561A (zh) | 一种光伏电站辐照度的超短期预测方法 | |
CN103927597A (zh) | 基于自回归滑动平均模型的风电功率超短期预测方法 | |
CN106779223A (zh) | 一种光伏***发电功率实时预测方法及装置 | |
CN103473322A (zh) | 基于时间序列模型的光伏发电功率超短期预测方法 | |
Raptis et al. | Total power quality index for electrical networks using neural networks | |
CN104037790B (zh) | 一种基于时序蒙特卡罗模拟的新能源接纳能力评估方法 | |
CN104573876A (zh) | 基于时序长记忆模型的风电场短期风速预测方法 | |
CN103984988A (zh) | 测光网络实时校正arma模型光伏功率超短期预测方法 | |
Wu et al. | Transmission network expansion planning based on chronological evaluation considering wind power uncertainties | |
CN103258103A (zh) | 基于偏最小二乘回归的戴维南等值参数辨识方法 | |
CN103984986A (zh) | 实时校正的自学习arma模型风电功率超短期预测方法 | |
Chen et al. | Distributional analysis for model predictive deferrable load control | |
Abdullah et al. | A noniterative method to estimate load carrying capability of generating units in a renewable energy rich power grid | |
Mohan et al. | Solar energy disaggregation using whole-house consumption signals | |
CN103984987B (zh) | 一种测风网络实时校正的arma模型风电功率超短期预测方法 | |
CN103927595A (zh) | 基于自学习复合数据源光伏发电功率超短期预测方法 | |
CN103927598A (zh) | 基于自学习径向基核函数的光伏发电功率预测方法 | |
Alzubaidi et al. | Identification of Suitable Probability Density Function for Wind Speed Profiles in Power System Studies | |
CN103927594A (zh) | 基于自学习复合数据源自回归模型的风电功率预测方法 | |
Maitanova et al. | Quantifying power and energy fluctuations of photovoltaic systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140716 |
|
RJ01 | Rejection of invention patent application after publication |