CN117370917B - Urban intelligent street lamp service life prediction method and system - Google Patents

Urban intelligent street lamp service life prediction method and system Download PDF

Info

Publication number
CN117370917B
CN117370917B CN202311668068.2A CN202311668068A CN117370917B CN 117370917 B CN117370917 B CN 117370917B CN 202311668068 A CN202311668068 A CN 202311668068A CN 117370917 B CN117370917 B CN 117370917B
Authority
CN
China
Prior art keywords
street lamp
illumination intensity
intelligent street
degree
moment
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
CN202311668068.2A
Other languages
Chinese (zh)
Other versions
CN117370917A (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.)
City Light Hunan Energy Saving And Environmental Service Co ltd
Original Assignee
City Light Hunan Energy Saving And Environmental Service 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 City Light Hunan Energy Saving And Environmental Service Co ltd filed Critical City Light Hunan Energy Saving And Environmental Service Co ltd
Priority to CN202311668068.2A priority Critical patent/CN117370917B/en
Publication of CN117370917A publication Critical patent/CN117370917A/en
Application granted granted Critical
Publication of CN117370917B publication Critical patent/CN117370917B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/243Classification techniques relating to the number of classes
    • G06F18/2433Single-class perspective, e.g. one-against-all classification; Novelty detection; Outlier detection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J1/00Photometry, e.g. photographic exposure meter
    • G01J1/10Photometry, e.g. photographic exposure meter by comparison with reference light or electric value provisionally void
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M11/00Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
    • G01M11/02Testing optical properties
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/17Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/10Pre-processing; Data cleansing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2123/00Data types
    • G06F2123/02Data types in the time domain, e.g. time-series data
    • 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
    • Y02BCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B20/00Energy efficient lighting technologies, e.g. halogen lamps or gas discharge lamps
    • Y02B20/40Control techniques providing energy savings, e.g. smart controller or presence detection

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Artificial Intelligence (AREA)
  • Algebra (AREA)
  • Evolutionary Computation (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Software Systems (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Circuit Arrangement For Electric Light Sources In General (AREA)

Abstract

The invention relates to the technical field of data prediction, in particular to a life prediction method and system of an urban intelligent street lamp, comprising the following steps: acquiring actual illumination intensity time sequence data and rated illumination intensity time sequence data; acquiring the abnormal degree of the illumination intensity of the intelligent street lamp at each moment, and acquiring the first abnormal degree and the second abnormal degree of the illumination intensity of the intelligent street lamp at each moment; acquiring a plurality of intervals to obtain the fluctuation period of data points of the intelligent street lamp at each moment, and obtaining the abnormal fluctuation degree of the intelligent street lamp at each moment and the third abnormal degree of the illumination intensity of the intelligent street lamp at each moment; and correcting the actual illumination intensity time sequence data of the intelligent street lamp according to the third abnormal degree of the illumination intensity of the intelligent street lamp at each moment to obtain corrected actual illumination intensity time sequence data of the intelligent street lamp, and finally predicting the service life of the intelligent street lamp. The invention improves the accuracy of data anomaly detection on the actual illumination intensity time sequence data.

Description

Urban intelligent street lamp service life prediction method and system
Technical Field
The invention relates to the technical field of data prediction, in particular to a life prediction method and system for an urban intelligent street lamp.
Background
With the development of cities, intelligent street lamps play a key role in providing road illumination and improving urban safety as an important component of modern urban illumination systems. In order to achieve efficient street lamp management, it becomes important to accurately predict the life of intelligent street lamps. In the life prediction method of the intelligent street lamp, the light intensity attenuation of the street lamp is important data. When the light intensity attenuation of the intelligent street lamp is calculated, the light intensity data of the street lamp are required to be collected, but due to a series of reasons such as sensor errors, the collected light intensity data have a series of data noise, the accuracy of data prediction is affected, and data cleaning is required.
In order to eliminate the influence of abnormal data on intelligent street lamp life prediction, a random initial threshold value is selected through an iterative threshold algorithm in conventional implementation, and a final optimal threshold value is obtained through continuous iteration to obtain an optimal group of abnormal data.
Disclosure of Invention
The invention provides a life prediction method and system for an urban intelligent street lamp, which aims to solve the existing problems.
The invention relates to a life prediction method and a life prediction system for an urban intelligent street lamp, which adopts the following technical scheme:
the embodiment of the invention provides a life prediction method of an urban intelligent street lamp, which comprises the following steps:
acquiring actual illumination intensity time sequence data of the intelligent street lamp and corresponding rated illumination intensity time sequence data;
obtaining the abnormal degree of the illumination intensity of the intelligent street lamp at each moment according to the actual illumination intensity time sequence data and the rated illumination intensity time sequence data of the intelligent street lamp, obtaining the first abnormal degree of the illumination intensity of the intelligent street lamp at each moment according to the abnormal degree of the illumination intensity of the intelligent street lamp at each moment and the used time of the intelligent street lamp at each moment, and obtaining the second abnormal degree of the illumination intensity of the intelligent street lamp at each moment according to the first abnormal degree of the illumination intensity of the intelligent street lamp at each moment;
obtaining an actual illumination intensity curve according to actual illumination intensity time sequence data of the intelligent street lamp, obtaining a plurality of intervals according to the actual illumination intensity curve, obtaining a fluctuation period of data points of the intelligent street lamp at each moment according to the plurality of intervals, obtaining abnormal fluctuation degree of the intelligent street lamp at each moment according to the fluctuation period of the data points of the intelligent street lamp at each moment, and obtaining a third abnormal degree of the illumination intensity of the intelligent street lamp at each moment according to the abnormal fluctuation degree of the intelligent street lamp at each moment and the second abnormal degree of the illumination intensity of the intelligent street lamp at each moment;
and correcting the actual illumination intensity time sequence data of the intelligent street lamp according to the third abnormal degree of the illumination intensity of the intelligent street lamp at each moment to obtain corrected actual illumination intensity time sequence data of the intelligent street lamp, and predicting the service life of the intelligent street lamp according to the corrected actual illumination intensity time sequence data of the intelligent street lamp.
Further, the specific acquisition steps of the abnormal degree of the illumination intensity of the intelligent street lamp at each moment are as follows:
the degree of abnormality of the illumination intensity of the intelligent street lamp at each moment is equal to the absolute value of the difference value between the actual illumination intensity data and the rated illumination intensity data at each moment.
Further, the calculation formula of the first abnormal degree of the illumination intensity of the intelligent street lamp at each moment is as follows:
in the method, in the process of the invention,indicating the abnormal degree of the illumination intensity of the intelligent street lamp at the kth moment, < >>Indicating the time that the intelligent street lamp has been used up to the kth moment,/>Indicating rated life of intelligent street lamp->For a preset threshold value, ++>And the first abnormal degree of the illumination intensity of the intelligent street lamp at the kth moment is represented.
Further, the calculation formula of the second abnormal degree of the illumination intensity of the intelligent street lamp at each moment is as follows:
in the method, in the process of the invention,indicating a first abnormality degree of the illumination intensity of the intelligent street lamp at the kth time,/for the intelligent street lamp>Mean value of first abnormality degree of illumination intensity of intelligent street lamp at all time points in data point neighborhood of kth time point, +.>Variance of first abnormality degree indicating illumination intensity of intelligent street lamp at all times in data point neighborhood of kth time, +.>Indicating a second degree of abnormality of the illumination intensity of the intelligent street lamp at the kth time,/for the intelligent street lamp>Representing absolute value symbols.
Further, the step of obtaining an actual illumination intensity curve according to the actual illumination intensity time sequence data of the intelligent street lamp and obtaining a plurality of intervals according to the actual illumination intensity curve comprises the following specific steps:
performing curve fitting on the actual illumination intensity time sequence data of the intelligent street lamp by using a least square method to obtain an actual illumination intensity curve; and dividing the intervals according to extreme points in the actual illumination intensity curve to obtain a plurality of intervals.
Further, the step of obtaining the fluctuation period of the data points of the intelligent street lamp at each moment according to a plurality of intervals comprises the following specific steps:
the time length of each interval is 2 times as long as the fluctuation period of each time data point in the interval.
Further, the calculation formula of the abnormal fluctuation degree of the intelligent street lamp at each moment is as follows:
in the method, in the process of the invention,indicates the fluctuation period of the data point at the kth moment, < >>Mean value of fluctuation period of all time data points in the data point neighborhood of kth time of intelligent street lamp,/>Representing absolute value symbols, ++>And (5) representing the abnormal fluctuation degree of the intelligent street lamp at the kth moment.
Further, the calculation formula of the third abnormal degree of the illumination intensity of the intelligent street lamp at each moment is as follows:
in the method, in the process of the invention,indicating a first abnormality degree of the illumination intensity of the intelligent street lamp at the kth time,/for the intelligent street lamp>Indicating a second degree of abnormality of the illumination intensity of the intelligent street lamp at the kth time,/for the intelligent street lamp>Indicating the degree of abnormal fluctuation of the intelligent street lamp at the kth time,/->Indicating the third abnormality degree of the illumination intensity of the intelligent street lamp at the kth time,/day>Representing a linear normalization function.
Further, the correcting the actual illumination intensity time sequence data of the intelligent street lamp according to the third abnormal degree of the illumination intensity of the intelligent street lamp at each moment to obtain the corrected actual illumination intensity time sequence data of the intelligent street lamp comprises the following specific steps:
(1) Selecting the average actual illumination intensity of the actual illumination intensity corresponding to the maximum value of the third abnormal degree and the actual illumination intensity corresponding to the minimum value of the third abnormal degree in the third abnormal degree of the intelligent street lamp at all times as initial thresholds in an iterative threshold algorithm, dividing the actual illumination intensity time sequence data into two groups according to the initial thresholds, marking the two groups as an array A and an array B, obtaining the average value of the third abnormal degree of the actual illumination intensity at all times in the array A, marking the average value asObtaining the average value of the third abnormal degree of the actual illumination intensity at all times in the array B, and marking the average value as +.>Obtain->And->Is recorded as the average value ofSelecting->The actual illumination intensity which is the smallest in difference from the third abnormality degree of the actual illumination intensity is taken as a new threshold +.>
(2) When the difference between the new threshold and the initial threshold is greater than or equal to the preset threshold G, then according to the new thresholdDividing the time sequence data of the actual illumination intensity into two groups, marking the two groups as an array C and an array D, obtaining the average value of the third abnormal degree of the actual illumination intensity at all times in the array C, marking the average value as +.>Obtaining the average value of the third abnormal degree of the actual illumination intensity at all times in the array D, and marking the average value as +.>Obtain->And->Is recorded as->Selecting->The actual illumination intensity which is the smallest in difference from the third abnormality degree of the actual illumination intensity is taken as a new threshold +.>
(3) And similarly, stopping iteration until the difference of the thresholds of two adjacent iterations is smaller than a preset threshold G, and taking the threshold of the last iteration as an optimal threshold;
dividing actual illumination intensity time sequence data of the intelligent street lamp into two groups according to an optimal threshold value in the actual illumination intensity, calculating the average value of all the actual illumination intensities in each group, marking one group with the largest average value of the actual illumination intensities as abnormal data, marking the other group as normal data, deleting all the abnormal data in the actual illumination intensity time sequence data of the intelligent street lamp, and marking the deleted position of the abnormal data as a point to be interpolated; the data of the points to be interpolated are obtained through linear interpolation by two normal data at two ends, so that the corrected real illumination intensity time sequence data of the intelligent street lamp is obtained.
The invention also provides a life prediction system of the urban intelligent street lamp, which comprises a memory, a processor and a computer program stored in the memory and capable of running on the processor, wherein the processor realizes any one of the steps when executing the computer program.
The technical scheme of the invention has the beneficial effects that: according to the invention, the first abnormal degree and the second abnormal degree of the illumination intensity of the intelligent street lamp at each moment are obtained by analyzing the actual illumination intensity time sequence data and the rated illumination intensity time sequence data of the intelligent street lamp, so that the accuracy of detecting the abnormal data in the actual illumination intensity time sequence data is improved; and acquiring an optimal threshold according to the third abnormal degree of the illumination intensity of the intelligent street lamp at each moment, correcting the actual illumination intensity time sequence data of the intelligent street lamp according to the optimal threshold, and predicting the service life of the intelligent street lamp according to the corrected actual illumination intensity time sequence data of the intelligent street lamp, so that the accuracy of data prediction is improved. In the process of acquiring the optimal threshold value through the iterative threshold algorithm, acquiring the initial threshold value through the third abnormal degree, determining the initial threshold value in the iterative threshold algorithm, reducing the iterative times, further reducing the calculated data quantity and shortening the time consumed by a computer in processing the data.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart showing the steps of a life prediction method for an urban intelligent street lamp.
Detailed Description
In order to further describe the technical means and effects adopted by the invention to achieve the preset aim, the following detailed description is given below of a life prediction method and system for an urban intelligent street lamp according to the invention, which are specific embodiments, structures, features and effects thereof, with reference to the accompanying drawings and preferred embodiments. In the following description, different "one embodiment" or "another embodiment" means that the embodiments are not necessarily the same. Furthermore, the particular features, structures, or characteristics of one or more embodiments may be combined in any suitable manner.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
The invention provides a service life prediction method and a service life prediction system for urban intelligent street lamps, which are specifically described below with reference to the accompanying drawings.
Referring to fig. 1, a flowchart of steps of a method for predicting life of an urban intelligent street lamp according to an embodiment of the invention is shown, the method includes the following steps:
step S001: and acquiring actual illumination intensity time sequence data of the intelligent street lamp and corresponding rated illumination intensity time sequence data.
It should be noted that, because the illumination intensity of the intelligent street lamp can change according to the change of the traffic flow of the road, and the rated illumination intensity of the street lamp can be obtained according to the data of the traffic flow of the road and the power of the street lamp, the illumination intensity of the street lamp is adjusted, and then the illumination intensity time sequence data of each intelligent street lamp needs to be acquired for analysis.
Specifically, the illumination intensity of each intelligent street lamp is acquired by using an illumination measuring instrument at a time interval of one minute, and all data in two days are acquired to obtain actual illumination intensity time sequence data of each intelligent street lamp; and obtaining the rated illumination intensity of each intelligent street lamp at each moment according to the road traffic flow at each acquisition moment, the power of the street lamp and other factors, obtaining the rated illumination intensity of each intelligent street lamp at all moments within two days, and recording the rated illumination intensity time sequence data of each intelligent street lamp. The obtaining of the rated illumination intensity of each intelligent street lamp is a known technology, and detailed description is omitted herein. In this embodiment, an intelligent street lamp is used for analysis.
So far, the actual illumination intensity time sequence data and the rated illumination intensity time sequence data of the intelligent street lamp are obtained.
Step S002: obtaining abnormal degrees of illumination intensity of the intelligent street lamp at each moment according to actual illumination intensity time sequence data and rated illumination intensity time sequence data of the intelligent street lamp, obtaining first abnormal degrees of illumination intensity of the intelligent street lamp at each moment according to the abnormal degrees of illumination intensity of the intelligent street lamp at each moment and the time that each intelligent street lamp is used at each moment, and obtaining second abnormal degrees of illumination intensity of the intelligent street lamp at each moment according to the first abnormal degrees of illumination intensity of the intelligent street lamp at each moment.
When the iterative threshold algorithm is used for identifying and detecting abnormal data, the difference between the actual illumination intensity and the corresponding rated illumination intensity is generally used for representing the abnormal degree of each data point, and the abnormal data is identified through threshold iteration of the abnormal degree.
(1) And obtaining the abnormal degree of the illumination intensity of the intelligent street lamp at each moment according to the actual illumination intensity time sequence data and the rated illumination intensity time sequence data of the intelligent street lamp.
It should be further noted that, in the actual situation, the rated illumination intensity of the intelligent street lamp is generally determined by some known parameters around the street lamp, but because the intelligent street lamp has loss in each element in the intelligent street lamp during the use process, the loss of each element in the intelligent street lamp cannot be specifically measured, so that a difference exists between the actual illumination intensity and the rated illumination intensity of the intelligent street lamp, therefore, the analysis is performed according to the difference between the actual illumination intensity and the rated illumination intensity in the intelligent street lamp, and the abnormal data point in the abnormal data point is obtained.
Specifically, the abnormal degree of the illumination intensity at each moment is obtained according to the actual illumination intensity time sequence data and the rated illumination intensity time sequence data of each intelligent street lamp, and is expressed as follows:
in the method, in the process of the invention,representing the actual illumination intensity at the kth time in the actual illumination intensity time sequence data,/for>Rated light intensity of kth moment in the time sequence data of the rated light intensity is represented by +.>Indicating the abnormal degree of the illumination intensity of the intelligent street lamp at the kth moment, < >>Representing absolute value symbols.
Wherein, when the difference between the actual illumination intensity and the rated illumination intensity of the intelligent street lamp at each moment is larger, the data point at the moment is abnormal. The difference represents the absolute value of the difference.
(2) The first abnormal degree of the illumination intensity of the intelligent street lamp at each moment is obtained according to the abnormal degree of the illumination intensity of the intelligent street lamp at each moment and the used time of each intelligent street lamp at each moment.
It should be noted that, in the use process of the intelligent street lamp, the service life of the intelligent street lamp is gradually attenuated, and along with the attenuation of the service life of the intelligent street lamp, the actual illumination intensity of the intelligent street lamp is gradually attenuated, and the difference between the intelligent street lamp and the rated illumination intensity is gradually increased, so that the degree of abnormality is larger, and therefore, the degree of abnormality of the street lamp needs to be corrected according to the service time of the intelligent street lamp and the rated service life of the intelligent street lamp.
Specifically, a threshold value is presetWherein the present embodiment is +.>The embodiment is not particularly limited, and is described by taking 0.8 as an example, wherein +.>Depending on the particular implementation. According to the abnormal degree of the illumination intensity of the intelligent street lamp at each moment, the service time of the intelligent street lamp and the rated service life of the intelligent street lamp, the first abnormal degree of the illumination intensity of the intelligent street lamp at each moment is obtained, and the first abnormal degree is expressed as follows:
in the method, in the process of the invention,indicating the abnormal degree of the illumination intensity of the intelligent street lamp at the kth moment, < >>Indicating the time that the intelligent street lamp has been used up to the kth moment,/>Indicating rated life of intelligent street lamp->Is a preset threshold value, namely a correction coefficient,and the first abnormal degree of the illumination intensity of the intelligent street lamp at the kth moment is represented.
When the service time of the intelligent street lamp is closer to the rated service life, the illumination intensity attenuation degree of the intelligent street lamp is larger, and corresponding data points are more biased to abnormal data points.
(3) Obtaining a second abnormal degree of the illumination intensity of the intelligent street lamp at each moment according to the first abnormal degree of the illumination intensity of the intelligent street lamp at each moment.
It should be noted that, during the use of the intelligent street lamp, the illumination intensity of the intelligent street lamp gradually decays along with the use time, but the intelligent street lamp may locally violate the overall performance due to accidental factors and environmental factors. The method is characterized in that the greater the fluctuation degree of data around a data point is, the greater the abnormality degree of the data is; the greater the fluctuation around the data, the greater the degree of abnormality of the data.
It should be further noted that specific data in terms of data fluctuation are expressed as: the greater the fluctuation of the degree of abnormality near the data point, the greater the variance of the degree of abnormality around the data point, which means that the greater the fluctuation of the degree of abnormality around the data point, the more abnormal the actual illumination intensity data here; meanwhile, the greater the difference between the average value of the abnormality degree of the data point and the surrounding abnormality degree, the greater the abnormality degree of the actual illumination intensity of the data point.
Specifically, a threshold N is preset, where the embodiment is described by taking n=20 as an example, and the embodiment is not specifically limited, where N may be determined according to the specific implementation situation. Taking each data point in the actual illumination intensity time sequence data as a center point, selecting 2N+1 data points, recording the data points as data in the neighborhood of each data point, analyzing according to the data in the neighborhood of each data point to obtain a second abnormal degree of the illumination intensity of the intelligent street lamp at each moment, and expressing the second abnormal degree as follows by a formula:
in the method, in the process of the invention,indicating a first abnormality degree of the illumination intensity of the intelligent street lamp at the kth time,/for the intelligent street lamp>Mean value of first abnormality degree of illumination intensity of intelligent street lamp at all time points in data point neighborhood of kth time point, +.>Variance of first abnormality degree indicating illumination intensity of intelligent street lamp at all times in data point neighborhood of kth time, +.>Indicating a second degree of abnormality of the illumination intensity of the intelligent street lamp at the kth time,/for the intelligent street lamp>Representing absolute value symbols.
Wherein,the difference between the first abnormal degree of the illumination intensity of the intelligent street lamp at the kth moment and the average value of the first abnormal degree of the illumination intensity at all moments in the corresponding neighborhood is represented, and the larger the difference is, the more abnormal the data point is represented; />The larger the value is for the variance of the first degree of abnormality of the illumination intensity at all times in the neighborhood, the more abnormal the corresponding data point is.
So far, the second abnormal degree of the illumination intensity of the intelligent street lamp at each moment is obtained.
Step S003: obtaining an actual illumination intensity curve according to actual illumination intensity time sequence data of the intelligent street lamp, obtaining a plurality of intervals according to the actual illumination intensity curve, obtaining fluctuation periods of data points of the intelligent street lamp at each moment according to the intervals, obtaining abnormal fluctuation degrees of the intelligent street lamp at each moment according to the fluctuation periods of the data points of the intelligent street lamp at each moment, and obtaining third abnormal degrees of the illumination intensity of the intelligent street lamp at each moment according to the abnormal fluctuation degrees of the intelligent street lamp at each moment and the second abnormal degrees of the illumination intensity of the intelligent street lamp at each moment.
(1) And obtaining an actual illumination intensity curve according to the actual illumination intensity time sequence data of the intelligent street lamp, and obtaining a plurality of intervals according to the actual illumination intensity curve.
It should be noted that, when the abnormal data points appear, they generally show irregular or show strong regularity due to some specific factors, that is, the actual illumination intensity time sequence data can be curve fitted to obtain a curve, the fluctuation period is determined according to the adjacent extreme points in the curve, and the abnormal data points are analyzed according to the fluctuation period.
Specifically, performing curve fitting on actual illumination intensity time sequence data of the intelligent street lamp through a quintic polynomial by using a least square method to obtain an actual illumination intensity curve; dividing the intervals according to extreme points in the actual illumination intensity curve, namely forming one interval by adjacent extreme points, and obtaining a plurality of intervals, wherein each interval is left-open and right-closed, and the first interval is left-closed and right-closed. Wherein the extreme points include a maximum point and a minimum point. However, the fitting is not particularly limited in this embodiment, and the practitioner may be according to the specific case.
(2) And obtaining the fluctuation period of the data point of the intelligent street lamp at each moment according to the plurality of intervals, and obtaining the abnormal fluctuation degree of the intelligent street lamp at each moment according to the fluctuation period of the data point of the intelligent street lamp at each moment.
It should be noted that specific data in terms of fluctuation frequency are expressed as: the greater the frequency of fluctuation of the degree of abnormality near the data point, the faster the fluctuation of the degree of abnormality around the data point, the less stable the degree of abnormality at this time, namely the faster the variation of the difference between the actual illumination intensity and the rated illumination intensity, the more abnormal the actual illumination intensity data acquired here; meanwhile, the greater the difference between the data period where the data point is located and the surrounding fluctuation period, the greater the degree of abnormality of the illumination intensity of the data point.
Specifically, a threshold value M is preset, where the embodiment is described by taking m=10 as an example, and the embodiment is not specifically limited, where M may be determined according to the specific implementation situation. The time length of each interval is 2 times as long as the fluctuation period of each time data point in the interval. And taking each data point in the actual illumination intensity time sequence data as a center point, selecting 2M+1 data points, and recording the data points as data in the neighborhood of each data point. And obtaining the abnormal fluctuation degree of the intelligent street lamp at each moment according to the fluctuation period of the data point at each moment. Expressed by the formula:
in the method, in the process of the invention,indicates the fluctuation period of the data point at the kth moment, < >>Mean value of fluctuation period of all time data points in the data point neighborhood of kth time of intelligent street lamp,/>Representing absolute value symbols, ++>And (5) representing the abnormal fluctuation degree of the intelligent street lamp at the kth moment.
Wherein,representing the degree of difference between the fluctuation period of the kth time data point and the average value of the fluctuation periods of all the time data points in the corresponding neighborhood, wherein the larger the value is, the more abnormal the period of the data point is; />The frequency of the data point at the kth time is represented, and the larger the value is, the faster the data change around the data point is, the faster the data fluctuation is, and the higher the degree of abnormality is.
(3) And obtaining a third abnormal degree of the illumination intensity of the intelligent street lamp at each moment according to the abnormal fluctuation degree of the intelligent street lamp at each moment and the second abnormal degree of the illumination intensity of the intelligent street lamp at each moment.
Obtaining a third abnormal degree of the illumination intensity of the intelligent street lamp at each moment according to the abnormal fluctuation degree of the intelligent street lamp at each moment, the second abnormal degree and the first abnormal degree of the illumination intensity of the intelligent street lamp at each moment, and expressing the third abnormal degree as follows by a formula:
in the method, in the process of the invention,indicating a first abnormality degree of the illumination intensity of the intelligent street lamp at the kth time,/for the intelligent street lamp>Indicating a second degree of abnormality of the illumination intensity of the intelligent street lamp at the kth time,/for the intelligent street lamp>Indicating the degree of abnormal fluctuation of the intelligent street lamp at the kth time,/->Indicating the third abnormality degree of the illumination intensity of the intelligent street lamp at the kth time,/day>Representation lineAnd (5) a sex normalization function.
So far, the third abnormal degree of the illumination intensity of the intelligent street lamp at each moment is obtained.
Step S004: and correcting the actual illumination intensity time sequence data of the intelligent street lamp according to the third abnormal degree of the actual illumination intensity of the intelligent street lamp at each moment to obtain corrected actual illumination intensity time sequence data of the intelligent street lamp, and predicting the service life of the intelligent street lamp according to the corrected actual illumination intensity time sequence data of the intelligent street lamp.
A threshold value G is preset, where the embodiment is described by taking g=300 as an example, and the embodiment is not specifically limited, where G may be determined according to the specific implementation situation.
(1) Selecting the average actual illumination intensity of the actual illumination intensity corresponding to the maximum value of the third abnormal degree and the actual illumination intensity corresponding to the minimum value of the third abnormal degree in the third abnormal degree of the intelligent street lamp at all times as initial thresholds in an iterative threshold algorithm, dividing the actual illumination intensity time sequence data into two groups according to the initial thresholds, marking the two groups as an array A and an array B, obtaining the average value of the third abnormal degree of the actual illumination intensity at all times in the array A, marking the average value asObtaining the average value of the third abnormal degree of the actual illumination intensity at all times in the array B, and marking the average value as +.>Obtain->And->Is recorded as the average value ofSelecting->Actual illumination with minimum difference from third abnormality degree of actual illumination intensityIntensity as new threshold ∈ ->
(2) When the difference between the new threshold and the initial threshold is greater than or equal to the preset threshold G, then according to the new thresholdDividing the time sequence data of the actual illumination intensity into two groups, marking the two groups as an array C and an array D, obtaining the average value of the third abnormal degree of the actual illumination intensity at all times in the array C, marking the average value as +.>Obtaining the average value of the third abnormal degree of the actual illumination intensity at all times in the array D, and marking the average value as +.>Obtain->And->Is recorded as->Selecting->The actual illumination intensity which is the smallest in difference from the third abnormality degree of the actual illumination intensity is taken as a new threshold +.>
(3) And similarly, stopping iteration until the difference of the thresholds of two adjacent iterations is smaller than a preset threshold G, and taking the threshold of the last iteration as an optimal threshold.
Wherein the difference represents an absolute value of the difference; the iterative threshold algorithm is a well-known technique and will not be described in detail here.
Dividing actual illumination intensity time sequence data of the intelligent street lamp into two groups according to an optimal threshold value in the actual illumination intensity, calculating average values of all the actual illumination intensities in each group, marking one group with the largest average value of the actual illumination intensity as abnormal data, marking the other group as normal data, deleting all the abnormal data in the actual illumination intensity time sequence data of the intelligent street lamp, and marking the deleted position of the abnormal data as a point to be interpolated.
Obtaining all points to be interpolated in the actual illumination intensity time sequence data of the intelligent street lamp, wherein the points to be interpolated exist in two types in the actual illumination intensity time sequence data of the intelligent street lamp, the first type is that one point to be interpolated is contained in the middle of two normal data points, the second type is that a plurality of continuous points to be interpolated are contained in the middle of the two normal data points, namely, the data of the middle point to be interpolated obtains the data value of each point to be interpolated through linear interpolation of the two normal data at two ends, and thus, the corrected actual illumination intensity time sequence data of the intelligent street lamp is obtained.
And predicting the service life of the intelligent street lamp by using an exponential decay model according to the corrected actual illumination intensity time sequence data of the intelligent street lamp.
The embodiment provides a life prediction system of an urban intelligent street lamp, which comprises a memory, a processor and a computer program stored in the memory and capable of running on the processor, wherein the processor realizes steps S001 to S004 when executing the computer program.
The above description is only of the preferred embodiments of the present invention and is not intended to limit the invention, but any modifications, equivalent substitutions, improvements, etc. within the principles of the present invention should be included in the scope of the present invention.

Claims (8)

1. A life prediction method of an urban intelligent street lamp is characterized by comprising the following steps:
acquiring actual illumination intensity time sequence data of the intelligent street lamp and corresponding rated illumination intensity time sequence data;
obtaining the abnormal degree of the illumination intensity of the intelligent street lamp at each moment according to the actual illumination intensity time sequence data and the rated illumination intensity time sequence data of the intelligent street lamp, obtaining the first abnormal degree of the illumination intensity of the intelligent street lamp at each moment according to the abnormal degree of the illumination intensity of the intelligent street lamp at each moment and the used time of the intelligent street lamp at each moment, and obtaining the second abnormal degree of the illumination intensity of the intelligent street lamp at each moment according to the first abnormal degree of the illumination intensity of the intelligent street lamp at each moment;
obtaining an actual illumination intensity curve according to actual illumination intensity time sequence data of the intelligent street lamp, obtaining a plurality of intervals according to the actual illumination intensity curve, obtaining a fluctuation period of data points of the intelligent street lamp at each moment according to the plurality of intervals, obtaining abnormal fluctuation degree of the intelligent street lamp at each moment according to the fluctuation period of the data points of the intelligent street lamp at each moment, and obtaining a third abnormal degree of the illumination intensity of the intelligent street lamp at each moment according to the abnormal fluctuation degree of the intelligent street lamp at each moment and the second abnormal degree of the illumination intensity of the intelligent street lamp at each moment;
correcting the actual illumination intensity time sequence data of the intelligent street lamp according to the third abnormal degree of the illumination intensity of the intelligent street lamp at each moment to obtain corrected actual illumination intensity time sequence data of the intelligent street lamp, and predicting the service life of the intelligent street lamp according to the corrected actual illumination intensity time sequence data of the intelligent street lamp;
the calculation formula of the first abnormal degree of the illumination intensity of the intelligent street lamp at each moment is as follows:
in the method, in the process of the invention,indicating the abnormal degree of the illumination intensity of the intelligent street lamp at the kth moment, < >>Indicating the time that the intelligent street lamp has been used up to the kth moment,/>Indicating rated life of intelligent street lamp->For a preset threshold value, ++>The first abnormal degree of the illumination intensity of the intelligent street lamp at the kth moment is represented;
the method comprises the following specific steps of:
(1) Selecting the average actual illumination intensity of the actual illumination intensity corresponding to the maximum value of the third abnormal degree and the actual illumination intensity corresponding to the minimum value of the third abnormal degree in the third abnormal degree of the intelligent street lamp at all times as initial thresholds in an iterative threshold algorithm, dividing the actual illumination intensity time sequence data into two groups according to the initial thresholds, marking the two groups as an array A and an array B, obtaining the average value of the third abnormal degree of the actual illumination intensity at all times in the array A, marking the average value asObtaining the average value of the third abnormal degree of the actual illumination intensity at all times in the array B, and marking the average value as +.>Obtain->And->Is recorded as->Selecting->The actual illumination intensity which is the smallest in difference from the third abnormality degree of the actual illumination intensity is taken as a new threshold +.>
(2) When the difference between the new threshold and the initial threshold is greater than or equal to the preset threshold G, then according to the new thresholdDividing the time sequence data of the actual illumination intensity into two groups, marking the two groups as an array C and an array D, obtaining the average value of the third abnormal degree of the actual illumination intensity at all times in the array C, marking the average value as +.>Obtaining the average value of the third abnormal degree of the actual illumination intensity at all times in the array D, and marking the average value as +.>Obtain->And->Is recorded as->Selecting->The actual illumination intensity which is the smallest in difference from the third abnormality degree of the actual illumination intensity is taken as a new threshold +.>
(3) And similarly, stopping iteration until the difference of the thresholds of two adjacent iterations is smaller than a preset threshold G, and taking the threshold of the last iteration as an optimal threshold;
dividing actual illumination intensity time sequence data of the intelligent street lamp into two groups according to an optimal threshold value in the actual illumination intensity, calculating the average value of all the actual illumination intensities in each group, marking one group with the largest average value of the actual illumination intensities as abnormal data, marking the other group as normal data, deleting all the abnormal data in the actual illumination intensity time sequence data of the intelligent street lamp, and marking the deleted position of the abnormal data as a point to be interpolated; the data of the points to be interpolated are obtained through linear interpolation by two normal data at two ends, so that the corrected real illumination intensity time sequence data of the intelligent street lamp is obtained.
2. The life prediction method of urban intelligent street lamps according to claim 1, wherein the specific acquisition steps of the degree of abnormality of the illumination intensity of the intelligent street lamps at each moment are as follows:
the degree of abnormality of the illumination intensity of the intelligent street lamp at each moment is equal to the absolute value of the difference value between the actual illumination intensity data and the rated illumination intensity data at each moment.
3. The urban intelligent street lamp life prediction method according to claim 1, wherein the calculation formula of the second degree of abnormality of the illumination intensity of the intelligent street lamp at each moment is:
in the method, in the process of the invention,indicating a first abnormality degree of the illumination intensity of the intelligent street lamp at the kth time,/for the intelligent street lamp>Mean value of first abnormality degree of illumination intensity of intelligent street lamp at all time points in data point neighborhood of kth time point, +.>Variance of first abnormality degree indicating illumination intensity of intelligent street lamp at all times in data point neighborhood of kth time, +.>Indicating a second degree of abnormality of the illumination intensity of the intelligent street lamp at the kth time,/for the intelligent street lamp>Representing absolute value symbols.
4. The method for predicting the life of an intelligent street lamp according to claim 1, wherein the steps of obtaining an actual illumination intensity curve according to the actual illumination intensity time sequence data of the intelligent street lamp and obtaining a plurality of intervals according to the actual illumination intensity curve include the following steps:
performing curve fitting on the actual illumination intensity time sequence data of the intelligent street lamp by using a least square method to obtain an actual illumination intensity curve; and dividing the intervals according to extreme points in the actual illumination intensity curve to obtain a plurality of intervals.
5. The urban intelligent street lamp life prediction method according to claim 1, wherein the step of obtaining the fluctuation period of the intelligent street lamp data point at each moment according to a plurality of intervals comprises the following specific steps:
the time length of each interval is 2 times as long as the fluctuation period of each time data point in the interval.
6. The urban intelligent street lamp life prediction method according to claim 1, wherein the calculation formula of the abnormal fluctuation degree of the intelligent street lamp at each moment is as follows:
in the method, in the process of the invention,indicates the fluctuation period of the data point at the kth moment, < >>Mean value of fluctuation period of all time data points in the data point neighborhood of kth time of intelligent street lamp,/>Representing absolute value symbols, ++>And (5) representing the abnormal fluctuation degree of the intelligent street lamp at the kth moment.
7. The urban intelligent street lamp life prediction method according to claim 1, wherein the calculation formula of the third degree of abnormality of the illumination intensity of the intelligent street lamp at each moment is:
in the method, in the process of the invention,indicating a first abnormality degree of the illumination intensity of the intelligent street lamp at the kth time,/for the intelligent street lamp>Indicating a second degree of abnormality of the illumination intensity of the intelligent street lamp at the kth time,/for the intelligent street lamp>The abnormal fluctuation degree of the intelligent street lamp at the kth moment is represented,indicating the third abnormality degree of the illumination intensity of the intelligent street lamp at the kth time,/day>Representing a linear normalization function.
8. A city intelligent street lamp life prediction system comprising a memory, a processor and a computer program stored in the memory and running on the processor, wherein the processor, when executing the computer program, implements the steps of a city intelligent street lamp life prediction method as claimed in any one of claims 1-7.
CN202311668068.2A 2023-12-07 2023-12-07 Urban intelligent street lamp service life prediction method and system Active CN117370917B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311668068.2A CN117370917B (en) 2023-12-07 2023-12-07 Urban intelligent street lamp service life prediction method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311668068.2A CN117370917B (en) 2023-12-07 2023-12-07 Urban intelligent street lamp service life prediction method and system

Publications (2)

Publication Number Publication Date
CN117370917A CN117370917A (en) 2024-01-09
CN117370917B true CN117370917B (en) 2024-02-23

Family

ID=89393253

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311668068.2A Active CN117370917B (en) 2023-12-07 2023-12-07 Urban intelligent street lamp service life prediction method and system

Country Status (1)

Country Link
CN (1) CN117370917B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117979508B (en) * 2024-03-05 2024-07-12 六班电气有限公司 Intelligent lamplight regulation and control method based on underground environment sensing
CN118297560B (en) * 2024-06-06 2024-08-30 山东征途信息科技股份有限公司 Smart city digitized data fusion method and system

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010283274A (en) * 2009-06-08 2010-12-16 Sharp Corp Failure cause analyzer, failure cause analysis method, failure cause analysis program, and program recording medium
CN109242519A (en) * 2018-09-25 2019-01-18 阿里巴巴集团控股有限公司 A kind of abnormal behaviour recognition methods, device and equipment
CN110377447A (en) * 2019-07-17 2019-10-25 腾讯科技(深圳)有限公司 A kind of abnormal deviation data examination method, device and server
JP2020047078A (en) * 2018-09-20 2020-03-26 株式会社Screenホールディングス Data processing method, data processing device, and data processing program
WO2021056724A1 (en) * 2019-09-23 2021-04-01 平安科技(深圳)有限公司 Anomaly detection method and apparatus, electronic device and storage medium
CN114925753A (en) * 2022-04-28 2022-08-19 南通东升灯饰有限公司 Use abnormity alarm system of LED floor lamp
CN115081673A (en) * 2022-05-05 2022-09-20 国家石油天然气管网集团有限公司 Oil and gas pipeline abnormity prediction method and device, electronic equipment and medium
CN115510302A (en) * 2022-11-16 2022-12-23 西北工业大学 Intelligent factory data classification method based on big data statistics
CN116070163A (en) * 2023-03-07 2023-05-05 深圳市特安电子有限公司 Indoor harmful gas concentration anomaly monitoring data processing method
CN116400639A (en) * 2023-06-08 2023-07-07 佛山市星际云数字科技有限公司 PLC (programmable logic controller) collected data intelligent cleaning method and system
CN116595470A (en) * 2023-07-17 2023-08-15 北京云摩科技股份有限公司 Method for monitoring chain load of scraper
CN116702081A (en) * 2023-08-07 2023-09-05 西安格蒂电力有限公司 Intelligent inspection method for power distribution equipment based on artificial intelligence

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11181894B2 (en) * 2018-10-15 2021-11-23 Uptake Technologies, Inc. Computer system and method of defining a set of anomaly thresholds for an anomaly detection model

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010283274A (en) * 2009-06-08 2010-12-16 Sharp Corp Failure cause analyzer, failure cause analysis method, failure cause analysis program, and program recording medium
JP2020047078A (en) * 2018-09-20 2020-03-26 株式会社Screenホールディングス Data processing method, data processing device, and data processing program
CN109242519A (en) * 2018-09-25 2019-01-18 阿里巴巴集团控股有限公司 A kind of abnormal behaviour recognition methods, device and equipment
CN110377447A (en) * 2019-07-17 2019-10-25 腾讯科技(深圳)有限公司 A kind of abnormal deviation data examination method, device and server
WO2021056724A1 (en) * 2019-09-23 2021-04-01 平安科技(深圳)有限公司 Anomaly detection method and apparatus, electronic device and storage medium
CN114925753A (en) * 2022-04-28 2022-08-19 南通东升灯饰有限公司 Use abnormity alarm system of LED floor lamp
CN115081673A (en) * 2022-05-05 2022-09-20 国家石油天然气管网集团有限公司 Oil and gas pipeline abnormity prediction method and device, electronic equipment and medium
CN115510302A (en) * 2022-11-16 2022-12-23 西北工业大学 Intelligent factory data classification method based on big data statistics
CN116070163A (en) * 2023-03-07 2023-05-05 深圳市特安电子有限公司 Indoor harmful gas concentration anomaly monitoring data processing method
CN116400639A (en) * 2023-06-08 2023-07-07 佛山市星际云数字科技有限公司 PLC (programmable logic controller) collected data intelligent cleaning method and system
CN116595470A (en) * 2023-07-17 2023-08-15 北京云摩科技股份有限公司 Method for monitoring chain load of scraper
CN116702081A (en) * 2023-08-07 2023-09-05 西安格蒂电力有限公司 Intelligent inspection method for power distribution equipment based on artificial intelligence

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种基于深度学习的异常数据清洗算法;匡俊搴;电子与信息学报;507-513 *
数据流上异常数据的在线检测与修正;王永利;徐宏炳;董逸生;钱江波;刘学军;;应用科学学报(第03期);40-45 *

Also Published As

Publication number Publication date
CN117370917A (en) 2024-01-09

Similar Documents

Publication Publication Date Title
CN117370917B (en) Urban intelligent street lamp service life prediction method and system
CN107092582B (en) Online abnormal value detection and confidence evaluation method based on residual posterior
CN111277332A (en) Performance state detection method and device of optical module in engineering application and electronic equipment
CN111709465A (en) Intelligent identification method for rough difference of dam safety monitoring data
CN116740946B (en) Highway traffic flow big data analysis processing method
CN116308963B (en) Government affair data analysis method and system
CN110837933A (en) Leakage identification method, device, equipment and storage medium based on neural network
CN110991376B (en) Feature extraction method for partial discharge type identification
CN117195007A (en) Heat exchanger performance prediction method and system
CN117251749A (en) Data processing method of Internet of things based on incremental analysis
CN118169580A (en) Smart phone accessory production test method and system
CN108399415B (en) Self-adaptive data acquisition method based on life cycle stage of equipment
CN116576553B (en) Data optimization acquisition method and system for air conditioner
CN117313014A (en) Real-time monitoring and early warning method for abnormal energy consumption data of kiln operation
CN112013910B (en) Drainage pipe network flow detection method and device, server and storage medium
CN117828443A (en) Water quality multi-index multi-step prediction method based on sequence decomposition and autocorrelation
CN115833843B (en) Vehicle operation monitoring data storage optimization method and management platform
CN112507285A (en) Drought detection method
CN116910655A (en) Intelligent ammeter fault prediction method based on device measurement data
CN117113247A (en) Drainage system abnormality monitoring method, equipment and storage medium based on two-classification and clustering algorithm
CN110455438A (en) Fiber optic temperature distribution detector fiber lengths are calibrated automatically and automatic test approach
CN114169590A (en) Reservoir warehousing runoff forecasting and correcting method and system based on multi-scenario division
CN111241693B (en) Life prediction system based on fatigue accumulation hypothesis
CN108303196B (en) Digital demodulation and regional temperature monitoring method and system for identical fiber bragg grating
CN114186734A (en) River flood prediction method capable of explaining research strategy based on deep learning model

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