CN108764528A - A kind of Daily treatment cost interval prediction method of allowed for influencing factors analysis - Google Patents

A kind of Daily treatment cost interval prediction method of allowed for influencing factors analysis Download PDF

Info

Publication number
CN108764528A
CN108764528A CN201810368777.1A CN201810368777A CN108764528A CN 108764528 A CN108764528 A CN 108764528A CN 201810368777 A CN201810368777 A CN 201810368777A CN 108764528 A CN108764528 A CN 108764528A
Authority
CN
China
Prior art keywords
prediction
load
interval
trend
item
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
Application number
CN201810368777.1A
Other languages
Chinese (zh)
Inventor
陈崇明
吴瀛
吴焕
谭炳源
姚栋方
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Maintenance and Test Center of Extra High Voltage Power Transmission Co
Original Assignee
Maintenance and Test Center of Extra High Voltage Power Transmission Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Maintenance and Test Center of Extra High Voltage Power Transmission Co filed Critical Maintenance and Test Center of Extra High Voltage Power Transmission Co
Priority to CN201810368777.1A priority Critical patent/CN108764528A/en
Publication of CN108764528A publication Critical patent/CN108764528A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Development Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a kind of Daily treatment cost interval prediction methods of allowed for influencing factors analysis, first, the trend curve of load data is extracted using median filtering method, then principal component analysis is carried out to impact factor and trend term load, several principal components are integrated out as input, trend term load realizes the trend term prediction of load as output in support vector machines;Then trend term is subtracted from load curve, obtains random entry, and interval prediction is carried out using Markov model to random entry;Then trend term predicted value is added with random entry forecast interval to get to the forecast interval value of load;The upper and lower envelope of historical load curve is extracted using cubic spline interpolation, point in sampling originally carries out Markov interval prediction, the upper and lower limit of forecast interval and upper and lower envelope are compared, calculate error, finally the error is added in original forecast interval to get to final forecast interval.Prediction result of the present invention is stablized, and precision of prediction is high.

Description

Daily maximum load interval prediction method based on influence factor analysis
Technical Field
The invention relates to the technical field of energy prediction, in particular to a daily maximum load interval prediction method based on influence factor analysis.
Background
Because various uncertain factors are contained in the power system, decision-making work necessarily faces risks to a certain degree, and uncertainty of power demand must be considered during decision-making. The results of traditional deterministic prediction methods do not reflect the uncertainty of the demand, and interval prediction can meet this objective requirement. The result of the interval prediction is not a simple deterministic value, but an interval, and the interval corresponds to a certain level of probability confidence level, which can describe the possible range of future prediction results. According to the interval prediction result, a power system decision maker can better recognize uncertainty and risk factors which may exist in future loads when performing production planning, system safety analysis and other work, so that a more reasonable decision can be made in time. Therefore, the method has important theoretical significance and practical value for analyzing the change rule of the load of the power system, researching the prediction method of the power load interval and realizing the uncertainty prediction of the power load.
The current interval prediction method has the following defects:
1. the calculation is complex;
2. the hypothesis is strong;
3. the calculation time is long.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide a daily maximum load interval prediction method based on influence factor analysis, so that the load of a power distribution area can be accurately predicted in the middle period.
In order to achieve the purpose, the technical scheme of the invention is as follows:
the daily maximum load interval prediction method based on influence factor analysis is characterized by comprising the following steps
Downloading historical load data to obtain a historical load curve;
finding out a trend curve of the historical load curve by using a median filtering method to obtain a trend item load;
carrying out principal component analysis on the influence factors and the trend item loads in the power system, integrating a plurality of principal components, taking the plurality of principal components as input and the trend item loads as output, and realizing the trend item prediction of the loads on a support vector machine;
subtracting the trend item load from the historical load curve to obtain a random item; averagely dividing the curve of the random item into a plurality of regions according to the maximum and minimum amplitude as the upper and lower limits, and calculating a transition probability matrix, wherein the region with the maximum transition probability is the prediction region of the random item;
carrying out interval prediction on the random item by using a Markov model to obtain a trend item prediction value of the load;
and adding the predicted value of the trend item of the load and the prediction interval of the random item to obtain an original prediction interval.
The method further comprises the following steps:
extracting upper and lower envelope lines of the historical load curve by a cubic spline interpolation method;
comparing the upper limit and the lower limit of the original prediction interval with the upper envelope line and the lower envelope line, and calculating an error;
and adding the error into the original prediction interval to obtain a final prediction interval.
The impact factors include the area covered by the load, and the industrial/commercial/residential power consumption ratio.
The process of analyzing the main components of the influence factors and the trend item loads in the power system and integrating a plurality of main components is as follows:
the first step is as follows: assuming that n samples and p variables are provided, the observation data matrix is:
wherein X is an observation data matrix, XijIs the value of the jth variable for the ith sample.
The second step is that: calculating a matrix of correlation coefficients for the samples
Assuming that the original data is still represented by X after normalization, the correlation coefficient of the normalized data is:
wherein r isijIs the correlation coefficient of the ith variable and the jth variable.
The third step: calculating characteristic value (lambda) of correlation coefficient matrix R by using Jacobian method12,…,λp) And the corresponding featuresEigenvector ai=[ai1,ai2,…,aip],i=1,2,…,p;
The fourth step: selecting important principal components, and writing a principal component expression;
the first k principal components are selected according to the magnitude of the accumulated contribution rate of each principal component, wherein the contribution rate refers to the proportion of the variance of a certain principal component to the total variance, namely the proportion of a certain characteristic value to the total characteristic value, namely
Wherein λ isiIs the characteristic value of the ith variable.
The larger the contribution rate is, the stronger the information of the original variable contained in the principal component is; the selection of the number k of the principal components is mainly determined according to the accumulated contribution rate of the principal components, and the accumulated contribution rate is required to reach more than 85 percent;
the principal component expression is as follows:
Zi=ai1×x1+ai2×x2+…+aip×xp
wherein Z isiIs the ith main component, xjIs the jth variable, aijIs the feature vector of the jth variable in the ith principal component.
Compared with the prior art, the invention has the beneficial effects that:
1. the prediction result is stable, and the prediction precision is high;
2. the prediction method is a data-driven and self-adaptive method, and the prediction result does not depend on the prior knowledge of a user;
3. the invention has simple, visual and easy operation;
4. the method has strong guiding significance on the actual situations of power grid dispatching, planning and the like.
Drawings
FIG. 1 is a flowchart of a daily maximum load interval prediction method based on impact factor analysis according to the present invention;
fig. 2 is a graph of ideal forecast results and actual data.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and detailed description.
Example (b):
as shown in fig. 1, the method for predicting the daily maximum load interval based on the impact factor analysis provided in this embodiment specifically includes the following steps:
1) and finding out a trend curve of the historical load curve by using a median filtering method to obtain a trend item load, and subtracting the trend item load from the historical load curve to obtain a random item. Assuming that the history curve f (t), t is 1,2, …, N and the trend curve is M (t), then
Where 2d +1 is the window length.
2) Carrying out principal component analysis on the influence factors (including the area covered by the load, the industrial/commercial/residential electricity utilization ratio and the like) and the trend item load to integrate a plurality of principal components;
the first step is as follows: assuming that n samples and p variables are provided, the observation data matrix is:
wherein X is an observation data matrix, XijIs the value of the jth variable for the ith sample.
The second step is that: calculating a matrix of correlation coefficients for the samples
Assuming that the original data is still represented by X after normalization, the correlation coefficient of the normalized data is:
wherein r isijIs the correlation coefficient of the ith variable and the jth variable.
The third step: calculating characteristic value (lambda) of correlation coefficient matrix R by using Jacobian method12,…,λp) And corresponding feature vector ai=[ai1,ai2,…,aip],i=1,2,…,p;
The fourth step: selecting important principal components, and writing a principal component expression;
however, because the variance of each principal component is decreased and the amount of information contained is also decreased, in the actual analysis, the first k principal components are selected according to the magnitude of the cumulative contribution rate of each principal component, instead of selecting p principal components, where the contribution rate refers to the proportion of the variance of a principal component to the total variance, and actually, the proportion of a feature value to the total feature value, that is, the proportion of the feature value to the total feature value
Wherein λ isiIs the characteristic value of the ith variable.
The larger the contribution rate is, the stronger the information of the original variable contained in the principal component is; the selection of the number k of the principal components is mainly determined according to the accumulated contribution rate of the principal components, and the accumulated contribution rate is generally required to reach more than 85 percent, so that the comprehensive variables can be ensured to include most of information of the original variables;
the principal component expression is as follows:
Zi=ai1×x1+ai2×x2+…+aip×xp
wherein Z isiIs the ith main component, xjIs the jth variable, aijIs the feature vector of the jth variable in the ith principal component.
3) Taking the principal component and the historical load data at the time t-1 as input, taking the historical load data at the time t as output, and training a model; based on the trained model, taking principal components and load data immediately before the prediction point as input, and predicting by the model to obtain a trend item prediction result of the load;
4) and averagely dividing the random item curve into 3 intervals according to the maximum and minimum amplitude as the upper limit and the lower limit, and calculating a transition probability matrix, wherein the interval with the maximum transition probability is the prediction interval of the random item.
Wherein omega1,Ω2,Ω3Represents interval 1, interval 2, interval 3.
The transition probability matrix can be expressed as
Wherein A isij(i-1, 2, 3; j-1, 2,3) indicates the probability that the value S (t-1) at the previous time is located in the i interval and the value S (t) at the next time is located in the j interval. Namely:
wherein, BijAnd (t) judging whether the value S (t-1) at the previous moment is in the interval i and the value S (t) at the next moment is in the variable of the interval j, if so, the value is 1, and if not, the value is 0.
If the load value of the day before the forecast is in the ith interval, forecasting the load random term interval S of the dayΩ(t) is Ai1,Ai2,Ai3The interval in which the maximum value of (1) is located.
5) The sum of the predicted value of the trend item of the load and the prediction interval of the random item is an original prediction interval.
Wherein,for trend term prediction, SΩ(t) a load random term interval of the prediction day,is the original prediction interval.
6) Finding out the upper and lower envelope lines y (t) ═ y (y) of the load curve by a Cubic Spline Interpolation (Spline Interpolation for short), which is a process of obtaining a curve function set by solving a three-bending moment equation set mathematically through a smooth curve of a series of shape value points)l(t)yh(t)];
7) And calculating the error between the upper envelope line and the lower envelope line and the upper limit and the lower limit of the original prediction interval.
8) The sum of the original prediction interval and the error is the final prediction interval.
In the prediction model, the prediction accuracy of the prediction method is evaluated in terms of both accuracy and precision.
1. Accuracy of
If the actual value of the predicted point falls between the upper limit and the lower limit of the prediction result, the prediction is accurate.
Accuracy c is defined as follows:
where m is the total number of predictions.
2. Accuracy of measurement
Adjusting the distance proportion:
the distance proportion is adjusted downwards:
taking the average absolute value, standard deviation, maximum value and minimum value of the up/down distance proportion as evaluation standards, and respectively marking the evaluation standards asstd(D1),D1max,D1minAndstd(D2),D2max,D2min
among the eight indexes, the smaller the value is, the smaller the threshold value between the predicted value and the true value is, the more stable the prediction is, and the accuracy of the prediction is embodied.
From the above analysis, it can be seen that there is a mutual constraint relationship between accuracy and precision, and a model with high accuracy will inevitably result in low precision, while a model with high precision will be limited in accuracy. How to make a trade-off selection between the two requires analysis of the specific situation.
Fig. 2 is a graph of ideal forecast results and actual data. As can be seen from the figure, in most cases, the true value is located between the upper limit and the lower limit of the prediction interval, and the accuracy is extremely high. The upper limit and the lower limit of the prediction interval are well fit with the fluctuation condition of a real curve, the error is small, and the accuracy is high.
The above embodiments are only for illustrating the technical concept and features of the present invention, and the purpose thereof is to enable those skilled in the art to understand the contents of the present invention and implement the present invention accordingly, and not to limit the protection scope of the present invention accordingly. All equivalent changes or modifications made in accordance with the spirit of the present disclosure are intended to be covered by the scope of the present disclosure.

Claims (4)

1. The daily maximum load interval prediction method based on influence factor analysis is characterized by comprising the following steps
Downloading historical load data to obtain a historical load curve;
finding out a trend curve of the historical load curve by using a median filtering method to obtain a trend item load;
carrying out principal component analysis on the influence factors and the trend item loads in the power system, integrating a plurality of principal components, taking the plurality of principal components as input and the trend item loads as output, and realizing the trend item prediction of the loads on a support vector machine;
subtracting the trend item load from the historical load curve to obtain a random item; averagely dividing the curve of the random item into a plurality of regions according to the maximum and minimum amplitude as the upper and lower limits, and calculating a transition probability matrix, wherein the region with the maximum transition probability is the prediction region of the random item;
carrying out interval prediction on the random item by using a Markov model to obtain a trend item prediction value of the load;
and adding the predicted value of the trend item of the load and the prediction interval of the random item to obtain an original prediction interval.
2. The impact factor analysis-based daily maximum load interval prediction method of claim 1, further comprising:
extracting upper and lower envelope lines of the historical load curve by a cubic spline interpolation method;
comparing the upper limit and the lower limit of the original prediction interval with the upper envelope line and the lower envelope line, and calculating an error;
and adding the error into the original prediction interval to obtain a final prediction interval.
3. The method for predicting the daily maximum load interval based on the impact factor analysis as claimed in claim 1 or 2, wherein the impact factors include the area covered by the load and the industrial/commercial/residential power consumption ratio.
4. The daily maximum load interval prediction method based on impact factor analysis according to claim 1 or 2, wherein the main component analysis is performed on the impact factors and the trend item loads in the power system, and the process of integrating a plurality of main components is as follows:
the first step is as follows: assuming that n samples and p variables are provided, the observation data matrix is:
wherein X is an observation data matrix, XijIs the value of the jth variable for the ith sample.
The second step is that: calculating a matrix of correlation coefficients for the samples
Assuming that the original data is still represented by X after normalization, the correlation coefficient of the normalized data is:
wherein r isijIs the correlation coefficient of the ith variable and the jth variable.
The third step: calculating characteristic value (lambda) of correlation coefficient matrix R by using Jacobian method12,…,λp) And corresponding feature vector ai=[ai1,ai2,…,aip],i=1,2,…,p;
The fourth step: selecting important principal components, and writing a principal component expression;
the first k principal components are selected according to the magnitude of the accumulated contribution rate of each principal component, wherein the contribution rate refers to the proportion of the variance of a certain principal component to the total variance, namely the proportion of a certain characteristic value to the total characteristic value, namely
Wherein λ isiIs the characteristic value of the ith variable.
The larger the contribution rate is, the stronger the information of the original variable contained in the principal component is; the selection of the number k of the principal components is mainly determined according to the accumulated contribution rate of the principal components, and the accumulated contribution rate is required to reach more than 85 percent;
the principal component expression is as follows:
Zi=ai1×x1+ai2×x2+…+aip×xp
wherein Z isiIs the ith main component, xjIs the jth variable, aijIs the feature vector of the jth variable in the ith principal component.
CN201810368777.1A 2018-04-23 2018-04-23 A kind of Daily treatment cost interval prediction method of allowed for influencing factors analysis Pending CN108764528A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810368777.1A CN108764528A (en) 2018-04-23 2018-04-23 A kind of Daily treatment cost interval prediction method of allowed for influencing factors analysis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810368777.1A CN108764528A (en) 2018-04-23 2018-04-23 A kind of Daily treatment cost interval prediction method of allowed for influencing factors analysis

Publications (1)

Publication Number Publication Date
CN108764528A true CN108764528A (en) 2018-11-06

Family

ID=64011598

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810368777.1A Pending CN108764528A (en) 2018-04-23 2018-04-23 A kind of Daily treatment cost interval prediction method of allowed for influencing factors analysis

Country Status (1)

Country Link
CN (1) CN108764528A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109801097A (en) * 2018-12-14 2019-05-24 深圳壹账通智能科技有限公司 Analysis method, device, storage medium and the analytical equipment of operation data
CN112288187A (en) * 2020-11-18 2021-01-29 国网湖南省电力有限公司 Big data-based electricity sales amount prediction method
CN112651545A (en) * 2020-12-07 2021-04-13 华电郑州机械设计研究院有限公司 Industrial heat load prediction method based on characteristic decomposition
CN115310212A (en) * 2022-10-12 2022-11-08 中汽研(天津)汽车工程研究院有限公司 Method for sampling characteristic data of automobile shock absorber

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1786740A (en) * 2005-11-15 2006-06-14 郑睿敏 Apparatus for controlling thickness of digital optical film
CN103860201A (en) * 2014-03-06 2014-06-18 西安交通大学 Method for extracting perfusion time intensity curve based on spread beam contrast imaging
CN105303468A (en) * 2015-11-20 2016-02-03 国网天津市电力公司 Comprehensive evaluation method of smart power grid construction based on principal component cluster analysis
CN106529197A (en) * 2016-12-07 2017-03-22 中国地质大学(武汉) Landslide stability time-varying law analysis method
CN106600038A (en) * 2016-11-18 2017-04-26 华南理工大学 Load interval prediction method based on Markov model
CN106600037A (en) * 2016-11-18 2017-04-26 华南理工大学 Multi-parameter auxiliary load forecasting method based on principal component analysis

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1786740A (en) * 2005-11-15 2006-06-14 郑睿敏 Apparatus for controlling thickness of digital optical film
CN103860201A (en) * 2014-03-06 2014-06-18 西安交通大学 Method for extracting perfusion time intensity curve based on spread beam contrast imaging
CN105303468A (en) * 2015-11-20 2016-02-03 国网天津市电力公司 Comprehensive evaluation method of smart power grid construction based on principal component cluster analysis
CN106600038A (en) * 2016-11-18 2017-04-26 华南理工大学 Load interval prediction method based on Markov model
CN106600037A (en) * 2016-11-18 2017-04-26 华南理工大学 Multi-parameter auxiliary load forecasting method based on principal component analysis
CN106529197A (en) * 2016-12-07 2017-03-22 中国地质大学(武汉) Landslide stability time-varying law analysis method

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109801097A (en) * 2018-12-14 2019-05-24 深圳壹账通智能科技有限公司 Analysis method, device, storage medium and the analytical equipment of operation data
CN112288187A (en) * 2020-11-18 2021-01-29 国网湖南省电力有限公司 Big data-based electricity sales amount prediction method
CN112651545A (en) * 2020-12-07 2021-04-13 华电郑州机械设计研究院有限公司 Industrial heat load prediction method based on characteristic decomposition
CN115310212A (en) * 2022-10-12 2022-11-08 中汽研(天津)汽车工程研究院有限公司 Method for sampling characteristic data of automobile shock absorber
CN115310212B (en) * 2022-10-12 2023-03-24 中汽研(天津)汽车工程研究院有限公司 Method for sampling characteristic data of automobile shock absorber

Similar Documents

Publication Publication Date Title
CN108764528A (en) A kind of Daily treatment cost interval prediction method of allowed for influencing factors analysis
CN110443002A (en) A kind of Deformation of Steep Slopes prediction technique and system
CN112365029B (en) Missing value processing method for air conditioner load prediction and air conditioner load prediction system
CN108877223A (en) A kind of Short-time Traffic Flow Forecasting Methods based on temporal correlation
CN106600037B (en) Multi-parameter auxiliary load prediction method based on principal component analysis
CN107765347A (en) A kind of Gaussian process returns and the short-term wind speed forecasting method of particle filter
CN107358318A (en) Based on GM(1,1)The urban power consumption Forecasting Methodology of model and Grey Markov chain predicting model
CN107944612B (en) Bus net load prediction method based on ARIMA and phase space reconstruction SVR
CN110751318A (en) IPSO-LSTM-based ultra-short-term power load prediction method
CN110232483A (en) Deep learning load forecasting method, device and terminal device
CN105321345A (en) Road traffic flow prediction method based on ARIMA model and kalman filtering
CN113130014B (en) Rare earth extraction simulation method and system based on multi-branch neural network
CN110428270A (en) The potential preference client recognition methods of the channel of logic-based regression algorithm
CN116307215A (en) Load prediction method, device, equipment and storage medium of power system
Zhao et al. On-line least squares support vector machine algorithm in gas prediction
CN113449919B (en) Power consumption prediction method and system based on feature and trend perception
CN109637128A (en) A kind of grey Verhulst Short-time Traffic Flow Forecasting Methods and system based on Markov
CN112149896A (en) Attention mechanism-based mechanical equipment multi-working-condition fault prediction method
CN115545333A (en) Method for predicting load curve of multi-load daily-type power distribution network
CN114548494A (en) Visual cost data prediction intelligent analysis system
CN106600038A (en) Load interval prediction method based on Markov model
Ballı et al. An application of artificial neural networks for prediction and comparison with statistical methods
CN112700050A (en) Method and system for predicting ultra-short-term 1 st point power of photovoltaic power station
CN108665090B (en) Urban power grid saturation load prediction method based on principal component analysis and Verhulst model
CN116484675A (en) Crack propagation life prediction method and system for ship engine blade

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20181106