CN105092082A - Fiber grating temperature measurement system and method for high overload distribution transformer - Google Patents

Fiber grating temperature measurement system and method for high overload distribution transformer Download PDF

Info

Publication number
CN105092082A
CN105092082A CN201510502708.1A CN201510502708A CN105092082A CN 105092082 A CN105092082 A CN 105092082A CN 201510502708 A CN201510502708 A CN 201510502708A CN 105092082 A CN105092082 A CN 105092082A
Authority
CN
China
Prior art keywords
transformer
high overload
fiber grating
time
grating temperature
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
CN201510502708.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.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd
State Grid Shandong Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd
State Grid Shandong 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 State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd, State Grid Shandong Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201510502708.1A priority Critical patent/CN105092082A/en
Publication of CN105092082A publication Critical patent/CN105092082A/en
Pending legal-status Critical Current

Links

Abstract

The invention discloses a fiber grating temperature measurement system and method for a high overload distribution transformer. The system comprises: a broadband light source; a controllable sawtooth wave generating circuit; an F-P cavity tunable filter, used for demodulating narrow band spectra generated by the broadband light source into astigmatic spectral signals according to the drive of a sawtooth wave voltage output by the sawtooth wave generating circuit; an optical splitter, used for splitting the astigmatic spectral signals output by the F-P cavity tunable filter to a required number of optical paths and communicating with a fiber grating temperature sensor; a calibration device, used for determining the relation between the voltage of the tunable light source and the central wavelength in real time; and a photoelectric detection module, used for converting reflected wave reflected by each sensing grating and uniquely corresponding to the sensing grating from an optical signal to an electrical signal and transmitting the electrical signal to a data acquisition card. By utilizing a quasi-distributed fiber grating temperature measurement technology, high-precision measurement of hot-spot temperature of the high overload distribution transformer is realized, and aging evaluation and remaining life cycle prediction of the high overload distribution transformer are realized.

Description

The optical fiber grating temperature-measuring system of high overload substation transformer and method
Technical field
The present invention relates to equipment for power transmission and distribution on-line monitoring technique, more particularly utilize quasi-distributed fiber grating temperature sensor, measure substation transformer hot(test)-spot temperature, and carry out the prediction of substation transformer life cycle.
Background technology
Rural power grids have load, the feature that seasonal strong, average load rate is low, peak-valley difference is large, from ruuning situation analysis, rural power grids are swift and violent in particular time (during as the Spring Festival, busy farming season) load growth, each department agriculture network distribution transformer because of short-time overload heating problem very outstanding.Inside transformer hot spot temperature of winding and its insulating property and aging life-span have direct relation, the permissible value exceeding temperature can reduce the serviceable life of transformer, especially the high temperature produced under overload ruuning situation can accelerate the insulation ag(e)ing of transformer, shortens its serviceable life.Carrying out that temperature monitoring and life prediction calculate to transformer is a kind of effective way improved transformer aging performance, increase economic efficiency.
Be applied to the method that transformer temperature directly measures at present to mainly contain: 1) electric signal sensor mensuration, as thermopair, thermal resistance etc., such sensor is subject to electromagnetic interference (EMI), and serviceable life is limited, and measurement effect is undesirable; 2) IR thermometry, the method is non-cpntact measurement, is convenient to manual operation, but cannot realize On-line sampling system, and is subject to ground unrest and electromagnetic environment interference; 3) based on the distributed optical fiber temperature measuring method of RAMAN scattering or Brillouin scattering principle, its temperature measurement accuracy, mutually to restrict between spatial positioning accuracy and response time, still can not meet the requirement of inside transformer accurate temperature measurement at present; 4) the point type thermometric of fluorescence optical fiber or Semiconductor GaAs temperature probe optical fiber is adopted, measurement point is few, such Fibre Optical Sensor is existing practical application in transformer thermometric, but experience shows, at normal winding top, diverse location place may have the temperature difference being greater than 10K, only puts into 1 to 3 sensors and is difficult to detect real focus, therefore, its layout in inside transformer often needs experience to estimate to there is certain limitation in actual applications.
Optical fiber grating temperature-measuring is by the fiber grating composition optical fibre optical grating sensing array of multiple different centre wavelength of connecting on simple optical fiber, pass through wavelength-division multiplex technique, realize the application of multiple sensor on simple optical fiber, the monitoring of the careful distribution of inside transformer multiple spot can be realized, the sensing measurement of fiber grating is completed by spectrum analysis simultaneously, and the impact by light intensity is minimum; By the position of reasonable Arrangement sensor, sensor is directly installed near measuring point, can realize the high-precision temperature monitoring that inside transformer is quasi-distributed completely.
In view of the advantage of optical fiber grating temperature-measuring and the development of technology, the Case comparison completing optical fiber grating temperature-measuring at present in high-power transformer inside is common, but the case completing optical fiber grating temperature-measuring in 10kV substation transformer inside is seldom seen in report.Consider the load characteristic of high overload distribution transforming particular time, the thermometry adopted should be able to reflect high overload distribution transforming coiling hotspot temperature rise continually varying process preferably in thermometric resolution, sensitivity and precision.Therefore, high overload substation transformer temperature monitoring system technical scheme should adopt optical fiber grating temperature-measuring technology.
In existing patent, application number is 2013102675399, disclose a kind of method of carrying out transformer life prediction based on optical fiber grating temperature-measuring system, carry out temperature measurement on-line and life prediction for Large Oil immersion transformer, because its load curve is comparatively level and smooth, variational regularity is comparatively strong, do not consider sensitivity and the precision of temp measuring system in load changing situation.Substation transformer is directly connected with user, its load curve undulatory property is larger, especially rural power grids distribution transformer load has erratic time varying characteristic, and existing patent " is carried out the method for transformer life prediction " in sensitivity of thermometry and precision and can not be reflected high overload distribution transforming coiling hotspot temperature rise continually varying process in real time, exactly based on optical fiber grating temperature-measuring system.
Summary of the invention
For solving the deficiency that prior art exists, the invention discloses optical fiber grating temperature-measuring system and the method for high overload substation transformer, the present invention utilizes quasi-distributed fiber grating thermometry, realize the high-acruracy survey to high overload substation transformer hot(test)-spot temperature, the differential equation of life loss is transformed to difference equation and reduces computational complexity, significantly improve counting yield, realize the high overload substation transformer residual life cycle quick, efficiently predict.
For achieving the above object, concrete scheme of the present invention is as follows:
The optical fiber grating temperature-measuring system of high overload substation transformer, comprising:
Wideband light source, for generation of the narrow-band spectrum of continuous sweep;
Controllable sawtooth circuit for generating, for producing sawtooth wave under the control of industrial computer;
F-P cavity adjustable filter, according to controllable sawtooth circuit for generating export sawtooth voltage drive to wideband light source produce narrow-band spectrum carry out demodulation, be demodulated into astigmatism spectrum signal;
Optical branching device, is divided into the light path of requirement for the astigmatism spectrum signal exported by F-P cavity adjustable filter and communicates with fiber-optical grating temperature sensor;
Calibrating installation, for determining the relation of voltage with F-P cavity outgoing wave centre wavelength of the tunable light source that saw-tooth wave generating circuit provides in real time;
Photoelectric detection module, the unique with it corresponding reflection wave be reflected back by each sensing grating converts electric signal to from light signal and is sent to data collecting card and carries out A/D conversion, and the data after A/D conversion are sent to industrial computer and process.
Based on the high overload substation transformer life cycle Forecasting Methodology of the optical fiber grating temperature-measuring system of above-mentioned high overload substation transformer, comprising:
The relative ageing rate of setting transformer, obtains transformer life loss within a certain period of time, obtains its differential equation form to transformer life loss distortion within a certain period of time;
Discrete hot(test)-spot temperature data are monitored according to the optical fiber grating temperature-measuring system of high overload substation transformer;
Relative ageing rate V in each time interval is asked for according to ageing rate computing formula relative in " GB1094.7-2008 power transformer the 7th part: oil-immersed power transformer load directive/guide " by monitoring discrete hot(test)-spot temperature data (n);
According to the relative ageing rate V in each time interval (n)and the interval time between each time period obtains at interval of the life loss difference equation form in the time period;
According to life loss total when superposing reach section tn zone time at interval of the life loss difference equation in the time period.
L = ∫ t 1 t 2 V d t
Further, the relative ageing rate V of transformer represents, then transformer life loss L is within a certain period of time expressed as with the relation differential equation form of relative ageing rate V:
d L d t = V
Discrete hot(test)-spot temperature data are monitored according to the optical fiber grating temperature-measuring system of high overload substation transformer:
Represent the interval time between each time period with Dt, hot(test)-spot temperature is obtained by optical fiber grating temperature-measuring system measurement, can be asked for the relative ageing rate V in each time interval by hot(test)-spot temperature θ h according to ageing rate computing formula relative in load directive/guide (n), be then expressed as at interval of the life loss difference equation form in the time period:
DL (n)=V (n)×Dt
When reaching section tn zone time, total life loss is:
L (n)=L (n-1)+DL (n)
In order to obtain an accurate result, time interval Dt should be little as far as possible, can not be greater than the half of the winding time constant (being determined by Transformer Plant Test data) in transformer thermal characteristic parameter.
Beneficial effect of the present invention:
The present invention considers sensitivity and the precision of temp measuring system in distribution transforming load changing situation, temp measuring system increases industrial computer to the feedback element of saw-tooth wave generating circuit, increase calibrating installation simultaneously, can reflect that when transformer load suddenlys change the temperature rise of high overload distribution transforming coiling hotspot changes in real time, exactly.Improvement has been done to the computing method of life loss simultaneously.
The present invention utilizes quasi-distributed fiber grating thermometry, realize the high-acruracy survey to high overload substation transformer hot(test)-spot temperature, obtain substation transformer hot(test)-spot temperature Monitoring Data, Aging Assessment and residual life period forecasting thereof are carried out to high overload substation transformer; During life appraisal, the differential equation of life loss is transformed to difference equation obviously to reduce computation complexity, improve counting yield, realize the high overload substation transformer residual life cycle quick, efficiently predict, scientific and effective substation transformer O&M strategy is formulated to O&M department there is great importance.
Accompanying drawing explanation
Fig. 1 one-piece construction schematic diagram of the present invention;
In figure, 1, wideband light source, 2, optical branching device, 3, controllable sawtooth circuit for generating, 4, F-P cavity fluid filter, 5, calibrating installation, 6, grating temperature sensor, 7, photoelectric detection module, 8, data collecting card, 9, passage, 10, shunting coupler, 11,1 × n coupling mechanism, 12, industrial computer.
Embodiment:
Below in conjunction with accompanying drawing, the present invention is described in detail:
As shown in Figure 1, the optical fiber grating temperature-measuring system of high overload substation transformer is primarily of compositions such as wideband light source 1, optical branching device 2, controllable sawtooth circuit for generating 3, F-P cavity fluid filter 4, calibrating installation 5, grating temperature sensor 6, photoelectric detection module 7 and data collecting cards 8.Wherein wideband light source 1 provides the narrow-band spectrum of continuous sweep for system; Optical branching device 2 effect is the light path that the continuous spectrum sent by light source is divided into requirement; F-P cavity fluid filter 4 after controllable sawtooth voltage driven for the signal receiving to sensing grating; Calibrating installation 5 is for determining the relation of voltage with centre wavelength of tunable light source in real time; Photoelectric detection module 7 effect converts electric signal to by from light signal; Data acquisition process unit 8 mainly completes A/D conversion, carries out signal transacting.Specific works principle is as follows:
Wideband light source 1 exports discrete spectrum signal as shown in the figure after F-P cavity fluid filter 4 is tuning, the pulsed light that this wavelength is connected enters the FBG array of sensor fibre grating sensor 6, and (this n light path is parallel transmission through coupling mechanism 11 points of n roads of 1 × n, and the centre wavelength of the grating sensor in each path can overlap, this improves the multiplexing capacity of sensor-based system to a certain extent).The pulsed light row that wavelength is connected are through the Bragg grating of each Tandem, the simple venation Wavelength matched with this raster center can be reflected wash off, when pulsed light mates completely with the centre wavelength of sensing grating reflection wave, the luminous power of its reflection wave reaches maximum value (in crest shape).Like this, each sensing grating can be reflected back unique corresponding reflection wave with it.These reflected light deliver to photoelectric detection module 7 again, photoelectric detection module 7 carries out opto-electronic conversion, amplify, after the process such as filtering, faint light signal is converted into the electric signal with certain amplitude, then send into data collecting card 8 and carry out A/D conversion, last industrial computer 12 processes the wavelength and temperature value that obtain each sensor to these data.
The relative ageing rate V of transformer represents, then transformer life loss L is within a certain period of time expressed as with the relation differential equation form of relative ageing rate V:
d L d t = V
The hot(test)-spot temperature Monitoring Data of substation transformer is as shown in the table:
Represent interval time time period with Dt, the relative ageing rate V of each time period can be asked for by hot(test)-spot temperature θ h according to relative ageing rate computing formula (n), then at interval of the life loss in the time period be:
DL (n)=V (n)×Dt
Total life loss is:
L (n)=L (n-1)+DL (n)
In order to obtain an accurate result, time interval Dt should be little as far as possible, can not be greater than the half of the winding time constant (being determined by Transformer Plant Test data) in transformer thermal characteristic parameter.
The hot spot temperature of winding data obtained by this system monitoring are more accurate, the differential equation of life loss is transformed to difference equation and reduces computational complexity, realize quick, the efficient calculation of high overload substation transformer life loss, thus scientificlly and effectively assess its residual life cycle, scientific and reasonable O&M strategy is formulated to O&M department there is great importance.
By reference to the accompanying drawings the specific embodiment of the present invention is described although above-mentioned; but not limiting the scope of the invention; one of ordinary skill in the art should be understood that; on the basis of technical scheme of the present invention, those skilled in the art do not need to pay various amendment or distortion that creative work can make still within protection scope of the present invention.

Claims (5)

1. the optical fiber grating temperature-measuring system of high overload substation transformer, is characterized in that, comprising:
Wideband light source, for generation of the narrow-band spectrum of continuous sweep;
Controllable sawtooth circuit for generating, for producing sawtooth wave under the control of industrial computer;
F-P cavity adjustable filter, according to controllable sawtooth circuit for generating export sawtooth voltage drive to wideband light source produce narrow-band spectrum carry out demodulation, be demodulated into astigmatism spectrum signal;
Optical branching device, is divided into the light path of requirement for the astigmatism spectrum signal exported by F-P cavity adjustable filter and communicates with fiber-optical grating temperature sensor;
Calibrating installation, for determining the relation of voltage with F-P cavity outgoing wave centre wavelength of the tunable light source that saw-tooth wave generating circuit provides in real time;
Photoelectric detection module, the unique with it corresponding reflection wave be reflected back by each sensing grating converts electric signal to from light signal and is sent to data collecting card and carries out A/D conversion, and the data after A/D conversion are sent to industrial computer and process.
2., based on the high overload substation transformer life cycle Forecasting Methodology of the optical fiber grating temperature-measuring system of the high overload substation transformer described in the claims 1, it is characterized in that, comprising:
The relative ageing rate of setting transformer, obtains transformer life loss within a certain period of time, obtains its differential equation form to transformer life loss distortion within a certain period of time;
Discrete hot(test)-spot temperature data are monitored according to the optical fiber grating temperature-measuring system of high overload substation transformer;
Relative ageing rate V in each time interval is asked for according to relative ageing rate computing formula by monitoring discrete hot(test)-spot temperature data (n);
According to the relative ageing rate V in each time interval (n)and the interval time between each time period obtains at interval of the life loss difference equation form in the time period;
According to life loss total when superposing reach section tn zone time at interval of the life loss difference equation in the time period.
3. the high overload substation transformer life cycle Forecasting Methodology of the optical fiber grating temperature-measuring system of high overload substation transformer as claimed in claim 2, it is characterized in that, the relative ageing rate V of transformer represents, then transformer life loss L is within a certain period of time expressed as with the relation differential equation form of relative ageing rate V:
d L d t = V .
4. the high overload substation transformer life cycle Forecasting Methodology of the optical fiber grating temperature-measuring system of high overload substation transformer as claimed in claim 2, it is characterized in that, monitor discrete hot(test)-spot temperature data according to the optical fiber grating temperature-measuring system of high overload substation transformer:
Represent the interval time between each time period with Dt, hot(test)-spot temperature is obtained by optical fiber grating temperature-measuring system measurement, can be asked for the relative ageing rate V in each time interval by hot(test)-spot temperature θ h according to ageing rate computing formula relative in load directive/guide (n), be then expressed as at interval of the life loss difference equation form in the time period:
DL (n)=V (n)×Dt。
5. the high overload substation transformer life cycle Forecasting Methodology of the optical fiber grating temperature-measuring system of high overload substation transformer as claimed in claim 4, is characterized in that, when reaching section tn zone time, total life loss is:
L (n)=L (n-1)+DL (n)
Time interval Dt can not be greater than the half of the winding time constant in transformer thermal characteristic parameter.
CN201510502708.1A 2015-08-14 2015-08-14 Fiber grating temperature measurement system and method for high overload distribution transformer Pending CN105092082A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510502708.1A CN105092082A (en) 2015-08-14 2015-08-14 Fiber grating temperature measurement system and method for high overload distribution transformer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510502708.1A CN105092082A (en) 2015-08-14 2015-08-14 Fiber grating temperature measurement system and method for high overload distribution transformer

Publications (1)

Publication Number Publication Date
CN105092082A true CN105092082A (en) 2015-11-25

Family

ID=54573018

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510502708.1A Pending CN105092082A (en) 2015-08-14 2015-08-14 Fiber grating temperature measurement system and method for high overload distribution transformer

Country Status (1)

Country Link
CN (1) CN105092082A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105865363A (en) * 2016-04-08 2016-08-17 南方电网科学研究院有限责任公司 Transformer winding deformation online test apparatus and method
CN108151910A (en) * 2017-12-26 2018-06-12 安徽盛美金属科技有限公司 A kind of high intensity temperature-measuring optical fiber grating sensor of transformer winding
CN113551792A (en) * 2021-07-07 2021-10-26 贵州乌江水电开发有限责任公司 Anti-interference capability improving method for PLC temperature module

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101949744A (en) * 2010-09-06 2011-01-19 国网电力科学研究院武汉南瑞有限责任公司 Fiber grating-based transformer internal temperature detection system
CN103364658A (en) * 2013-06-28 2013-10-23 国网电力科学研究院武汉南瑞有限责任公司 Method for predicting service life of transformer based on fiber grating temperature measurement system
CN103512510A (en) * 2013-10-23 2014-01-15 中国电子科技集团公司第三十四研究所 Fiber bragg grating sensing system based on narrow-band scanning light source and operation method
CN103837178A (en) * 2013-11-29 2014-06-04 湖北工业大学 System and method for demodulating fiber bragg gratings based on liquid crystal F-P cavity adjustable filter technology
CN105004445A (en) * 2015-07-30 2015-10-28 国网山东省电力公司 Mobile distribution transformer temperature monitoring device based on infrared communication and method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101949744A (en) * 2010-09-06 2011-01-19 国网电力科学研究院武汉南瑞有限责任公司 Fiber grating-based transformer internal temperature detection system
CN103364658A (en) * 2013-06-28 2013-10-23 国网电力科学研究院武汉南瑞有限责任公司 Method for predicting service life of transformer based on fiber grating temperature measurement system
CN103512510A (en) * 2013-10-23 2014-01-15 中国电子科技集团公司第三十四研究所 Fiber bragg grating sensing system based on narrow-band scanning light source and operation method
CN103837178A (en) * 2013-11-29 2014-06-04 湖北工业大学 System and method for demodulating fiber bragg gratings based on liquid crystal F-P cavity adjustable filter technology
CN105004445A (en) * 2015-07-30 2015-10-28 国网山东省电力公司 Mobile distribution transformer temperature monitoring device based on infrared communication and method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
杨英松: "变压器绕组热电温度与油色谱在线监测应用研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105865363A (en) * 2016-04-08 2016-08-17 南方电网科学研究院有限责任公司 Transformer winding deformation online test apparatus and method
CN108151910A (en) * 2017-12-26 2018-06-12 安徽盛美金属科技有限公司 A kind of high intensity temperature-measuring optical fiber grating sensor of transformer winding
CN108151910B (en) * 2017-12-26 2019-12-31 嘉善品智联科技有限公司 High-strength temperature measurement fiber grating sensor for transformer winding
CN113551792A (en) * 2021-07-07 2021-10-26 贵州乌江水电开发有限责任公司 Anti-interference capability improving method for PLC temperature module

Similar Documents

Publication Publication Date Title
CN102047092B (en) Beam path monitoring device, and beam path monitoring system
CN103017804B (en) High time synchronized multi-path optical fiber grating sensing system
CN102914321B (en) Ultra-low fiber bragg grating sensing system and query method thereof
CN103364658A (en) Method for predicting service life of transformer based on fiber grating temperature measurement system
CN100507455C (en) Intensity modulation type optical fiber sensor multiplexing method
CN103512510A (en) Fiber bragg grating sensing system based on narrow-band scanning light source and operation method
CN101925795A (en) Method for sensing strain in component in wind turbine, optical strain sensing system and uses thereof
US8734011B2 (en) Distributed optical fiber temperature sensor based on optical fiber delay
CN101000267A (en) Parallel distribution optical fibre raster temp. sensing method and its system
CN102818531B (en) Dynamic strain measurement instrument based on multiple overlapped gratings
CN103616091B (en) A kind of distributed fiber optic temperature and stress sensing device
CN203605976U (en) Distributed type optical fiber temperature and stress sensing device
CN103398801A (en) Optical fiber grating temperature measuring device and measuring method
CN103245369A (en) Novel fiber bragg grating demodulation method and system thereof based on multi-longitudinal mode F-P laser device
CN104457803A (en) Wave length demodulating system and method based on F-P etalon and reference gratings
CN106352905A (en) Fiber bragg grating demodulating system and demodulator
CN101876573B (en) Array waveguide grating-based temperature sensing method and temperature sensor
CN105092082A (en) Fiber grating temperature measurement system and method for high overload distribution transformer
CN103162742A (en) Novel distributed optical fiber temperature and stress sensor
CN103399191A (en) Sideband demodulation-based FBG-GMM (Fiber Bragg Grating-Giant Magnetostrictive Material) current sensor
CN203534538U (en) A fiber grating sensing system based on a narrowband scan light source
CN108196107A (en) A kind of distributed full photovoltage measuring system and its measuring method
CN206930377U (en) Based on narrow band laser demodulating fiber bragg grating temp measuring system
CN103389120B (en) Optical fiber based on wavelength-division multiplex declines and swings chamber sensing array apparatus
CN109612601B (en) Power equipment temperature and partial discharge integrated detection system and method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20151125

RJ01 Rejection of invention patent application after publication