CN106650060B - Photovoltaic cell internal resistance attenuation coefficient prediction method - Google Patents

Photovoltaic cell internal resistance attenuation coefficient prediction method Download PDF

Info

Publication number
CN106650060B
CN106650060B CN201611124348.7A CN201611124348A CN106650060B CN 106650060 B CN106650060 B CN 106650060B CN 201611124348 A CN201611124348 A CN 201611124348A CN 106650060 B CN106650060 B CN 106650060B
Authority
CN
China
Prior art keywords
internal resistance
attenuation coefficient
photovoltaic cell
time sequence
photovoltaic
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.)
Active
Application number
CN201611124348.7A
Other languages
Chinese (zh)
Other versions
CN106650060A (en
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.)
Shenyang University of Technology
State Grid Qinghai Electric Power Co Ltd
Electric Power Research Institute of State Grid Qinghai Electric Power Co Ltd
Original Assignee
Shenyang University of Technology
State Grid Qinghai Electric Power Co Ltd
Electric Power Research Institute of State Grid Qinghai Electric Power Co Ltd
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 Shenyang University of Technology, State Grid Qinghai Electric Power Co Ltd, Electric Power Research Institute of State Grid Qinghai Electric Power Co Ltd filed Critical Shenyang University of Technology
Priority to CN201611124348.7A priority Critical patent/CN106650060B/en
Publication of CN106650060A publication Critical patent/CN106650060A/en
Application granted granted Critical
Publication of CN106650060B publication Critical patent/CN106650060B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • 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
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02SGENERATION OF ELECTRIC POWER BY CONVERSION OF INFRARED RADIATION, VISIBLE LIGHT OR ULTRAVIOLET LIGHT, e.g. USING PHOTOVOLTAIC [PV] MODULES
    • H02S50/00Monitoring or testing of PV systems, e.g. load balancing or fault identification
    • H02S50/10Testing of PV devices, e.g. of PV modules or single PV cells
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/04Ageing analysis or optimisation against ageing
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Economics (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Geometry (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Photovoltaic Devices (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a method for predicting an attenuation coefficient of internal resistance of a photovoltaic cell, which comprises the following steps: according to the parameters obtained by real-time monitoring, establishing a time sequence of a photovoltaic cell internal resistance attenuation coefficient evolution system, and establishing a photovoltaic cell internal resistance attenuation coefficient equation; performing phase space reconstruction processing on the time sequence of the evolution system; processing a battery internal resistance attenuation coefficient equation according to a wavelet network method; and carrying out prediction calculation on the internal resistance attenuation coefficient of the photovoltaic cell on the time sequence after the phase space reconstruction. The technical problems that power grid and distributed photovoltaic power generation operation data resources cannot be effectively utilized, and the evaluation accuracy and the photovoltaic utilization efficiency are low are solved. The method has the advantages of improving the reliability of the model, improving the photovoltaic utilization rate, improving the accuracy of evaluation and improving the reliability and economy of the power distribution network power system after the photovoltaic system is connected. The method is applied to detecting the photovoltaic power generation system.

Description

Photovoltaic cell internal resistance attenuation coefficient prediction method
Technical Field
The invention belongs to the technical field of photovoltaic power generation, and particularly relates to a method for predicting an attenuation coefficient of internal resistance of a photovoltaic cell.
Background
Distributed photovoltaic power generation equipment and a power distribution network in a power system form a complex system, the output efficiency of a photovoltaic cell is influenced by the internal resistance of the photovoltaic cell, and the internal resistance of the photovoltaic cell is influenced by a plurality of influence factors. The photovoltaic cell internal resistance calculation method in the prior art has the technical problems that the influences of important factors such as the investment age and the number of the cell series-parallel connection assemblies are neglected, the power grid and the distributed photovoltaic power generation operation data resources cannot be effectively utilized, and the evaluation accuracy and the photovoltaic utilization efficiency are not high. The method considers multiple influence factors, monitors the power distribution network and the operating parameters and meteorological environment parameters of the photovoltaic system in the power distribution network in real time, carries out prediction calculation on the internal resistance attenuation coefficient of the distributed photovoltaic cell according to the monitoring parameters, controls the photovoltaic power generation system and the power distribution network in real time according to the calculation result, can effectively improve the reliability of the model, and thus greatly improves the photovoltaic utilization efficiency.
Disclosure of Invention
The invention provides a method for predicting the attenuation coefficient of the internal resistance of a photovoltaic cell, which solves the technical problems that the influence of important factors such as the investment age and the number of series-parallel components of the cell is neglected, the power grid and the distributed photovoltaic power generation operation data resources cannot be effectively utilized, and the evaluation accuracy and the photovoltaic utilization efficiency are not high.
The invention is realized by the following technical scheme, and the method comprises the following steps: step 1: according to the parameters obtained by real-time monitoring, establishing a time sequence of a photovoltaic cell internal resistance attenuation coefficient evolution system, and establishing a photovoltaic cell internal resistance attenuation coefficient equation; step 2: performing phase space reconstruction processing on the time sequence of the evolution system in the step 1; and step 3: processing a battery internal resistance attenuation coefficient equation according to a wavelet network method; and 4, step 4: and (3) carrying out prediction calculation on the internal resistance attenuation coefficient of the photovoltaic cell on the time sequence after the phase space reconstruction in the step (2).
Further, in order to better implement the method, the parameters obtained by real-time monitoring are the operation parameters of the photovoltaic system and the meteorological environment parameters in the power distribution network and the power distribution network.
Further, the evolution time sequence in step 1 is an evolution time sequence established at fixed time intervals.
Further, the evolution time sequence comprises the equivalent impedance measurement value of the photovoltaic power station access point, the access point voltage, the access point active value, the ambient temperature and the ambient illumination intensity.
Further, the evolving system time series in step 1 is at a series of time instants ts1,ts2,ts3,...tsnComprises the following steps:
Figure GDA0002268647330000021
the photovoltaic cell internal resistance attenuation coefficient equation is as follows:
Figure GDA0002268647330000022
wherein n is a natural number, n is 1,2, Tr is the measured battery input time, Tr is the external temperature, Sr is the external illumination, Cr is the number of series components, and Br is the number of parallel components.
Further, the step 2 comprises the following steps:
A. establishing an optimization objective model
Figure GDA0002268647330000023
Wherein i is 1,25n
B. Constructing the time sequence { rx of the evolution system in the step 1iThe m-dimensional phase space rx ofi+1=ψ (rxi,rxi- τ,...,rxi-(m-1) τ) Wherein i ═ 1,2.. k5nτ is the delay time and m is the embedding dimension.
Further, the step 3 comprises the following steps:
A. calculating an output layer of a wavelet network model of the attenuation coefficient of the internal resistance of the battery;
B. and (4) online correction of the wavelet network model of the attenuation coefficient of the internal resistance of the battery.
Further, step a in step 3 includes the following steps:
let the time sequence of the input signal be { rcxiTherein, by
Figure GDA0002268647330000025
Calculating the output value of the hidden layer, and calculating the output value of the output layer according to the output value of the hidden layer
Figure GDA0002268647330000026
Wherein i 1,25n,j=1,2...l,ga(xzi)≥0,ga(xzi) Phi (j) is the output of the jth node of the hidden layer in the wavelet network, fjAs wavelet basis functions, αjIs fjOf the scaling factor, λjIs fjA translation factor of wijIs the size of the interconnection between the input layer and the hidden layer, l is the number of nodes in the hidden layer, wjRepresenting the connection weight between the hidden layer and the output layer.
Further, step B in step 3 includes the following steps:
according to er=yr-yccWavelet network model for on-line correction of attenuation coefficient of internal resistance of battery in factrTo predict the output, yccIs the actual measurement.
Further, the step 4 comprises the following steps:
according to the corrected wavelet neural network in the step 3, the time sequence after phase space reconstruction in the step 2 is subjected to prediction calculation of internal resistance attenuation coefficient of the photovoltaic cell, and target model correction is introducedCondition, optimization of the objective function to ya=minfar(rxi)+gar(rxi);
Wherein i is 1,25n,gar(rxi)≥0,gar(rxi) Is a constraint term of the object model, yaThe attenuation coefficient of the internal resistance of the photovoltaic cell.
Further, the method is used for detecting the photovoltaic power generation system.
Drawings
Fig. 1 is a prediction flow chart.
Detailed Description
The present invention will be described in further detail with reference to examples, but the embodiments of the present invention are not limited thereto.
Example 1:
the method for predicting the attenuation coefficient of the internal resistance of the photovoltaic cell is adopted, the prediction process is shown in figure 1, and the method comprises the following steps:
step 1:
according to the parameters, the input time, the number of the battery series-parallel connection assemblies and the standard parameters obtained by real-time monitoring, at a series of moments ts1,ts2,ts3,...tsnThe time sequence for establishing the photovoltaic cell internal resistance attenuation coefficient evolution system is as follows:
Figure GDA0002268647330000031
establishing an attenuation coefficient equation of the internal resistance of the photovoltaic cell:
Figure GDA0002268647330000041
wherein n is a natural number, n is 1,2, Tr is the measured battery input time, Tr is the external temperature, Sr is the external illumination, Cr is the number of series components, and Br is the number of parallel components.
Step 2: and (3) performing phase space reconstruction processing on the time sequence of the evolution system in the step 1:
step 2.1: establishing optimizationsObject model
Figure GDA0002268647330000042
Wherein i is 1,25n
Step 2.2: constructing the time sequence { rx of the evolution system in the step 1iThe m-dimensional phase space rx ofi+1=ψ(rxi,rxi-τ,...,rxi-(m-1)τ) Wherein i ═ 1,2.. k5n,τ=0.0152,m=5。
And step 3: processing a battery internal resistance attenuation coefficient equation according to a wavelet network method:
step 3.1: calculating an output layer of the wavelet network model of the attenuation coefficient of the internal resistance of the battery:
let the time sequence of the input signal be { rcxiTherein, by
Figure GDA0002268647330000043
Calculating the output value of the hidden layer, and calculating the output value of the output layer according to the output value of the hidden layer
Figure GDA0002268647330000044
Wherein i 1,25n,j=1,2...l,ga(xzi)≥0,ga(xzi) Phi (j) is the output of the jth node of the hidden layer in the wavelet network, fjAs wavelet basis functions, αjIs fjOf the scaling factor, λjIs fjA translation factor of wijIs the size of the interconnection between the input layer and the hidden layer, l is the number of nodes in the hidden layer, wjRepresenting the connection weight between the hidden layer and the output layer.
Step 3.2: on-line correction of a wavelet network model of the attenuation coefficient of the internal resistance of the battery:
according to er=yr-yccWavelet network model for on-line correction of attenuation coefficient of internal resistance of battery in factrTo predict the output, yccFor actual measured values, to ensure model accuracy, correction values erThe best effect is achieved when the concentration is less than or equal to 0.0001.
And 4, step 4: and (3) carrying out prediction calculation on the internal resistance attenuation coefficient of the photovoltaic cell on the time sequence after the phase space reconstruction in the step (2):
according to the wavelet neural network corrected in the step 3, the time sequence after phase space reconstruction in the step 2 is subjected to prediction calculation of the internal resistance attenuation coefficient of the photovoltaic cell, target model correction conditions are introduced, and a target function is optimized to be ya=minfar(rxi)+gar(rxi)。
Wherein i is 1,25n,gar(rxi)≥0,gar(rxi) For the constraint term of the object model,
Figure GDA0002268647330000051
yathe attenuation coefficient of the internal resistance of the photovoltaic cell.
The method is used for detecting the photovoltaic power generation system.
Compared with the prior art, the invention can obtain the following beneficial technical effects: (1) the reliability of the model is improved, (2) the photovoltaic utilization rate is improved, (3) the evaluation accuracy is improved, (4) the utilization rates of the power grid and distributed photovoltaic power generation operation data resources are improved, and (5) the reliability and the economy of the power system of the power distribution network after the photovoltaic system is connected are improved.
The above description is only a preferred embodiment of the present invention, and is not intended to limit the present invention in any way, and all simple modifications and equivalent variations of the above embodiments according to the technical spirit of the present invention are included in the scope of the present invention.

Claims (1)

1. A method for predicting an attenuation coefficient of internal resistance of a photovoltaic cell is characterized by comprising the following steps: the method comprises the following steps:
(1) according to parameters, investment time, the number of the battery series-parallel connection assemblies and standard parameters obtained through real-time monitoring, a time sequence of a photovoltaic cell internal resistance attenuation coefficient evolution system is established, and a photovoltaic cell internal resistance attenuation coefficient equation is established;
(2) performing phase space reconstruction processing on the time sequence of the evolution system in the step (1);
(3) processing a battery internal resistance attenuation coefficient equation according to a wavelet network method;
(4) performing prediction calculation on the internal resistance attenuation coefficient of the photovoltaic cell on the time sequence after the phase space reconstruction in the step (2);
the parameters obtained by real-time monitoring are power distribution network and photovoltaic system operation parameters and meteorological environment parameters in the power distribution network;
the evolution time sequence in the step (1) is an evolution time sequence established at a fixed time interval; the evolution time sequence comprises an equivalent impedance measured value of the photovoltaic power station access point, access point voltage, an access point active value, ambient temperature and ambient illumination intensity;
the evolving system time sequence in step (1) is at a series of times ts1,ts2,ts3,...tsnComprises the following steps:
Figure FDA0002268647320000011
the photovoltaic cell internal resistance attenuation coefficient equation is as follows:
Figure FDA0002268647320000012
wherein n is a natural number, n is 1,2, Tr is the measured battery input time, Tr is the external temperature, Sr is the external illumination, Cr is the number of series components, and Br is the number of parallel components;
the step (2) comprises the following steps:
(A) establishing an optimized target model minfar(rx1,rx2...rxi..rxh5n) K, where i is 1,25n
(B) Constructing the time sequence { rx of the evolution system in the step (1)iThe m-dimensional phase space rx ofi+1=ψ(rxi,rxi-τ,...,rxi-(m-1)τ) Wherein i ═ 1,2.. k5nτ is the delay time, m is the embedding dimension;
the step (3) comprises the following steps:
(A) calculating an output layer of a wavelet network model of the attenuation coefficient of the internal resistance of the battery;
(B) on-line correction of a wavelet network model of the attenuation coefficient of the internal resistance of the battery;
the step (A) includes the steps of:
let the time sequence of the input signal be { rcxiTherein, by
Figure FDA0002268647320000021
Calculating the output value of the hidden layer, and calculating the output value of the output layer according to the output value of the hidden layer
Figure FDA0002268647320000022
Wherein i 1,25n,j=1,2...l,ga(xzi)≥0,ga(xzi) Phi (j) is the output of the jth node of the hidden layer in the wavelet network, fjAs wavelet basis functions, αjIs fjOf the scaling factor, λjIs fjA translation factor of wijIs the size of the interconnection between the input layer and the hidden layer, l is the number of nodes in the hidden layer, wjRepresenting the connection weight between the hidden layer and the output layer;
wherein, yrIs a prediction output;
the step (B) comprises the steps of:
according to er=yr-yccWavelet network model for on-line correction of attenuation coefficient of internal resistance of battery in factrTo predict the output, yccIs an actual measured value;
the step (4) comprises the following steps:
according to the wavelet neural network corrected in the step (3), the time sequence after the phase space reconstruction in the step (2) is subjected to photovoltaic cell internal resistance attenuation coefficient prediction calculation, target model correction conditions are introduced, and a target function is optimized to be ya=minfar(rxi)+gar(rxi);
Wherein i is 1,25n,gar(rxi)≥0,gar(rxi) Is a constraint term of the object model, yaThe attenuation coefficient of the internal resistance of the photovoltaic cell.
CN201611124348.7A 2016-12-08 2016-12-08 Photovoltaic cell internal resistance attenuation coefficient prediction method Active CN106650060B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611124348.7A CN106650060B (en) 2016-12-08 2016-12-08 Photovoltaic cell internal resistance attenuation coefficient prediction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611124348.7A CN106650060B (en) 2016-12-08 2016-12-08 Photovoltaic cell internal resistance attenuation coefficient prediction method

Publications (2)

Publication Number Publication Date
CN106650060A CN106650060A (en) 2017-05-10
CN106650060B true CN106650060B (en) 2020-05-15

Family

ID=58819411

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611124348.7A Active CN106650060B (en) 2016-12-08 2016-12-08 Photovoltaic cell internal resistance attenuation coefficient prediction method

Country Status (1)

Country Link
CN (1) CN106650060B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107301266B (en) * 2017-05-15 2020-11-10 中山职业技术学院 LOC estimation method and system for lithium iron phosphate battery
CN107359860B (en) * 2017-06-28 2019-04-16 河海大学常州校区 A kind of perovskite solar cell electron lifetime test method based on EIS analysis
CN110458343B (en) * 2019-07-26 2023-04-07 国网山东省电力公司泰安供电公司 Method for predicting regional photovoltaic power generation capacity attenuation

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101572277A (en) * 2009-02-16 2009-11-04 武湖月 Opto-electrical energy material base and high-efficiency conversion
CN104836256A (en) * 2015-05-29 2015-08-12 国家电网公司 Calculation method and system of photovoltaic consumption capability of power distribution network
KR20150144425A (en) * 2014-06-16 2015-12-28 파워칩스주식회사 Electrical Energy Storing Device
CN105375874A (en) * 2015-10-30 2016-03-02 国家电网公司 Photovoltaic power generation maximum power tracing performance index prediction method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101572277A (en) * 2009-02-16 2009-11-04 武湖月 Opto-electrical energy material base and high-efficiency conversion
KR20150144425A (en) * 2014-06-16 2015-12-28 파워칩스주식회사 Electrical Energy Storing Device
CN104836256A (en) * 2015-05-29 2015-08-12 国家电网公司 Calculation method and system of photovoltaic consumption capability of power distribution network
CN105375874A (en) * 2015-10-30 2016-03-02 国家电网公司 Photovoltaic power generation maximum power tracing performance index prediction method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于小波分解的某些非平稳时间序列预测方法;徐科,等;《电子学报》;20010430;第566-568页 *

Also Published As

Publication number Publication date
CN106650060A (en) 2017-05-10

Similar Documents

Publication Publication Date Title
CN110070226B (en) Photovoltaic power prediction method and system based on convolutional neural network and meta-learning
Niu et al. Uncertainty modeling for chaotic time series based on optimal multi-input multi-output architecture: Application to offshore wind speed
CN106991504B (en) Building energy consumption prediction method and system based on subentry measurement time sequence and building
CN106650060B (en) Photovoltaic cell internal resistance attenuation coefficient prediction method
CN106548410B (en) Method for evaluating voltage unbalance probability of power distribution network containing distributed power supply
CN112668611B (en) Kmeans and CEEMD-PE-LSTM-based short-term photovoltaic power generation power prediction method
Kaplan et al. A novel method based on Weibull distribution for short-term wind speed prediction
Chauhan et al. Approximating parameters of photovoltaic models using an amended reptile search algorithm
CN111797132B (en) Multi-renewable energy power station power scene generation method considering space-time correlation
CN112508279A (en) Regional distributed photovoltaic prediction method and system based on spatial correlation
CN115453193A (en) Power distribution network harmonic state estimation method based on cooperation of PQM, TTU and SM measurement data
Zhang et al. Distribution system state estimation via data-driven and physics-aware deep neural networks
CN117473384A (en) Power grid line safety constraint identification method, device, equipment and storage medium
Zhu et al. Wind Speed Short-Term Prediction Based on Empirical Wavelet Transform, Recurrent Neural Network and Error Correction
Zhang et al. Short‐Term Power Prediction of Wind Power Generation System Based on Logistic Chaos Atom Search Optimization BP Neural Network
CN116400266A (en) Transformer fault detection method, device and medium based on digital twin model
CN116579386A (en) Wind turbine generator system state monitoring method, electronic equipment and storage medium
CN116258057A (en) Method for carrying out photovoltaic power generation predictive analysis based on Transformer
CN115952921A (en) Photovoltaic energy power prediction method and device, electronic equipment and storage medium
CN115907198A (en) Long-distance heat supply load intelligent prediction system
CN114298300A (en) Uncertainty prediction method, uncertainty prediction device, electronic device, and storage medium
CN107528312B (en) Power system state estimation method
Zhang et al. Wind speed prediction for wind farm based on clayton copula function and deep learning fusion
CN104037756A (en) Electric power system stability evaluation method including complex electric-power device model
Pan et al. Short-term wind speed prediction model of VMD-FOAGRNN based on Lorenz disturbance

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
GR01 Patent grant
GR01 Patent grant