CN102607807B - Method for precisely measuring transit time of optical fiber Sagnac interferometer - Google Patents
Method for precisely measuring transit time of optical fiber Sagnac interferometer Download PDFInfo
- Publication number
- CN102607807B CN102607807B CN201210047158.5A CN201210047158A CN102607807B CN 102607807 B CN102607807 B CN 102607807B CN 201210047158 A CN201210047158 A CN 201210047158A CN 102607807 B CN102607807 B CN 102607807B
- Authority
- CN
- China
- Prior art keywords
- optical fiber
- sagnac interferometer
- transit time
- fiber sagnac
- groups
- 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.)
- Expired - Fee Related
Links
- 239000013307 optical fiber Substances 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 24
- 230000003287 optical effect Effects 0.000 claims abstract description 15
- 238000005070 sampling Methods 0.000 claims abstract description 10
- 238000006243 chemical reaction Methods 0.000 claims description 10
- 230000008878 coupling Effects 0.000 claims description 8
- 238000010168 coupling process Methods 0.000 claims description 8
- 238000005859 coupling reaction Methods 0.000 claims description 8
- 238000000691 measurement method Methods 0.000 claims 4
- 239000000835 fiber Substances 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000003111 delayed effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000003466 welding Methods 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Landscapes
- Gyroscopes (AREA)
Abstract
The invention discloses a method for precisely measuring transit time of an optical fiber Sagnac interferometer. The method comprises the steps that: a second photoelectric detector is connected at a fourth tail fiber of a coupler of the optical fiber Sagnac interferometer; the modulated phase difference introduced by an integrated phase modulator is fixed and the optical power of a light source is changed; two photoelectric detectors convert the light intensity returned by the optical fiber and the light intensity output directly by the light source respectively to electric signals, and two analogue-to-digital converters based on one identical clock are used to sample the output signals of the two photoelectric detectors respectively; and two groups of sampling sequences obtained are processed relatively to obtain the time delay between the two sequences and calculate the transverse time of the optical fiber Sagnac interferometer with the combination of the sampling rate. According to the method, the transit time of the optical fiber Sagnac interferometer can be effectively measured, and a foundation is laid for the system to accurately select the period of the modulation signal, therefore, excursion and error caused by deviation of the transit time are prevented, and the long time stability and the precision of the optical fiber Sagnac interferometer are improved.
Description
Technical field
The present invention relates to the disposal route of optical fiber Sagnac interferometer signal, especially relate to a kind of method of optical fiber Sagnac interferometer transit time precision measurement.
Background technology
1913 Frenchman's Sagnac (G.Sagnac) invented a kind of rotatable ring interferometer.The light beam that same light source is sent is decomposed into two bundles, allows them in same annular light path, along contrary direction, propagate after one week and to join, and then on screen, produces interference, Here it is Sagnac effect.Interferometer based on Sagnac effect is called optics Sagnac interferometer, and this is wherein the most representative is exactly optical fiber Sagnac interferometer instrument.
Transit time is that light is propagated one week needed time in having the light path of reciprocity.Middle high-precision optical fiber Sagnac interferometer mainly adopts the scheme of digital closed loop, first by integrated phase modulator, optical system is modulated, then by analog to digital conversion converter, the signal of optical fiber Sagnac interferometer is converted to digital quantity, by controller, carries out modulation /demodulation processing.When semiperiod of integrated phase modulators modulate signal and optical fiber Sagnac interferometer transit time are when unequal, system can produce the skew relevant with this error; This skew correspondingly changes when environmental change, thereby at optical fiber Sagnac interferometer, goes out the error component of the drift that superposes on signal, has reduced static accuracy and the resolution of system, and has weakened stability.Only have and accurately measure the actual optical fiber Sagnac interferometer transit time, then according to measurement result, be used for setting the modulation signal cycle, could improve to greatest extent the drift performance of optical fiber Sagnac interferometer.
The transit time of optical fiber Sagnac interferometer is different numerical value to different optical systems, its concrete size is decided by optical path length and the refractive index in reciprocity light path, in the production run of winding optical fiber Sagnac interferometer fiber optic loop, can obtain the estimated value of fiber lengths precision 10 meters of left and right, this estimated value precision is about 50 nanosecond orders.In high-precision optical fiber Sagnac interferometer, need more accurate numerical value, and the reason due to production technology, in the production run of optical fiber Sagnac interferometer, keep away the welding of unavoidable optical fiber, also can cause the numerical value change of transit time, so the transit time of precisely measuring optical fiber Sagnac interferometer is very necessary.
Summary of the invention
For the demand to the precisely measuring optical fiber Sagnac interferometer transit time in current high-precision optical fiber Sagnac interferometer research, and lack the present situation of simple effective ways, the object of the present invention is to provide the method for precisely measuring optical fiber Sagnac interferometer transit time a kind of.
Inventive principle:
Optical fiber Sagnac interferometer output signal under the modulation of integrated phase modulator is:
I
1(t)=I
0{1+cos[Δφ
m(t)]} (3)
In this measuring method, the phase modulation of integrated phase modulator is poor is set as zero, has Δ φ
m(t)=0, can be obtained by (3) formula:
I
1(t)=I
0{1+cos[Δφ
m(t)]}=I
0[1+cos(0)]=2I
0 (4)
This is the light intensity that described first photodetector detects.
The luminous power of second photodetector direct detection light source output, supposes that coupling mechanism splitting ratio is 1: 1, and its light intensity detecting is:
I
2(t)=2I
0=I
1(t) (5)
This is described second light intensity that photodetector detects.
Below by obtaining and processing with the method for precisely measuring optical fiber Sagnac interferometer transit time and do a concise and to the point discussion two detector electric signal.
Suppose the light intensity signal of described light source exporting change, so first photodetector and second light intensity signal I that photodetector detects respectively
1and I (t)
2(t) also change, the detector output electrical signals of answering is in contrast S thereupon
1and S (t)
2(t).Theoretically, S
1and S (t)
2(t) be two signals that have a certain constant time lag, this time delay is the transit time of optical fiber Sagnac interferometer.Electric signal to detector output carries out analog to digital conversion, and sampling rate is f
s, obtain two groups of sample sequence S
1and S (n)
2(n).Known in signal correction treatment theory, by asking for the correlation integral of two signals, can calculate the time delay between two signals.
From the computing formula (1) of relevant treatment, different Δ n values is to having a specific relevant accumulated value R (Δ n), [0, N in transit time environs
d] relevant accumulated value R (Δ n) corresponding to all delayed data point Δ n form an A (R), has a maximum value R in this sequence
m, the Δ n of its correspondence
mbe two groups of the strongest time points of sample sequence correlativity, Δ n
m/ f
snamely time delay, the i.e. transit time of optical fiber Sagnac interferometer.In order to improve the measuring accuracy of transit time, can improve sample rate f
s, work as sample rate f
swhile being set as 1GHz, transit time measuring accuracy corresponding to the present invention can reach 1ns.
The step of the technical solution used in the present invention is as follows:
Comprise optical fiber Sagnac interferometer, it is characterized in that the step of the method is as follows:
1) at the 4th tail optical fiber place of optical fiber Sagnac interferometer coupling mechanism, connect another photodetector, the phase modulation that integrated phase modulator is introduced is poor fixing, changes source device output optical power, and two photodetectors convert light intensity signal to electric signal;
2) to step 1) in the output signal of two photodetectors use the analog to digital converter based on same clock to sample, obtain two groups of sample sequences;
3) to step 2) in two groups of sample sequences carry out relevant treatment, obtain two groups of sample sequence correlativitys corresponding delay point when maximum, according to sampling rate, calculate the transit time of optical fiber Sagnac interferometer.
The phase differential that described integrated phase modulator is introduced is fixed; The implementation of described change source device output optical power is that light source powers up moment to stablize the process of Output optical power or the light source output process with additional modulation signal variation luminous power.
Described based on same clock, sampling rate is f
stwo analog to digital converters respectively photodetector output electrical signals is separately carried out to analog to digital conversion, obtain two groups of sample sequence S
1and S (n)
2(n);
To two groups of described sample sequence S
1and S (n)
2(n) carry out relevant treatment, formula is as follows:
Wherein, N
sfor sample sequence length, Δ n is the delay point of two groups of sample sequences, and span is [0, N
d], (N
d/ f
s> τ
estimate, τ wherein
estimateestimated value for the optical fiber Sagnac interferometer transit time), obtain the sequence that corresponding accumulated value R (Δ n) forms
, the maximal value R of this sequence
mrepresent, the delay point that its corresponding Δ n is two groups of sample sequence correlativity maximums, is designated as Δ n
m, now the transit time τ of optical fiber Sagnac interferometer calculates by following formula:
τ=Δn
m/f
s (2)
From above formula, the transit time τ measuring accuracy of optical fiber Sagnac interferometer depends on sample rate f
s, improve sample rate f
scan promote the measuring accuracy of transit time.
The beneficial effect that the present invention has is:
The present invention is the transit time of precisely measuring optical fiber Sagnac interferometer effectively, for system, select exactly the modulation signal cycle foundation is provided, skew and error that inhibition causes because of transit time deviation, long-time stability and the precision of raising optical fiber Sagnac interferometer.
Accompanying drawing explanation
Fig. 1 is the system principle diagram of described precisely measuring optical fiber Sagnac interferometer transit time method.
Fig. 2 carries out to two groups of sample sequences the algorithm flow chart that relevant treatment obtains the transit time.
Fig. 3 is the result curve schematic diagram after two groups of sample sequences and relevant treatment.
Embodiment
Below in conjunction with drawings and Examples, the invention will be further described.
As shown in Figure 1, comprise coupling mechanism C2, integrated phase modulator C3, fiber optic loop C4, the first detector C 5, the second detector C 6, the first analog to digital conversion channel C 7, the second analog to digital conversion channel C 8, clock C9, relevant treatment algorithm C10 and the transit time C11 that light source C1, splitting ratio are 1: 1; Splitting ratio is first being connected with light source C1 and the first detector C 5 respectively with the 3rd tail optical fiber of coupling mechanism C2 of 1: 1, and third and fourth tail optical fiber of the coupling mechanism C2 that splitting ratio is 1: 1 is connected with integrated phase modulator C3 and the second detector C 6 respectively; The output terminal of integrated phase modulator C3 is connected with fiber optic loop C4; The first detector C 5 obtains sample sequence S after first number converter C7
1(n); The second detector C 6 obtains sample sequence S after the second analog to digital converter C8
2(n), computer based relevant treatment algorithm C10 is to S
1and S (n)
2(n) process, obtain two groups of sample sequence correlativitys corresponding time delay point when maximum, then according to sampling rate, calculate the transit time C11 of optical fiber Sagnac interferometer.
The step of the method is as follows:
1) at the 4th tail optical fiber place of optical fiber Sagnac interferometer coupling mechanism, connect another photodetector, the phase modulation that integrated phase modulator is introduced is poor fixing, changes source device output optical power, and two photodetectors convert light intensity signal to electric signal;
2) to step 1) in the output signal of two photodetectors use the analog to digital converter based on same clock to sample, obtain two groups of sample sequences;
3) to step 2) in two groups of sample sequences carry out relevant treatment, obtain two groups of sample sequence correlativitys corresponding delay point when maximum, according to sampling rate, calculate the transit time of optical fiber Sagnac interferometer.
The phase differential that described integrated phase modulator is introduced is fixed; The implementation of described change source device output optical power is that light source powers up moment to stablize the process of Output optical power or the light source output process with additional modulation signal variation luminous power.
Described based on same clock, sampling rate is f
stwo analog to digital converters respectively photodetector output electrical signals is separately carried out to analog to digital conversion, obtain two groups of sample sequence S
1and S (n)
2(n);
To two groups of described sample sequence S
1and S (n)
2(n) carry out relevant treatment, formula is as follows:
Wherein, N
sfor sample sequence length, Δ n is the delay point of two groups of sample sequences, and span is [0, N
d], (N
d/ f
s> τ
estimate, τ wherein
estimateestimated value for the optical fiber Sagnac interferometer transit time), obtain the sequence that corresponding accumulated value R (Δ n) forms
, the maximal value R of this sequence
mrepresent, the delay point that its corresponding Δ n is two groups of sample sequence correlativity maximums, is designated as Δ n
m, now the transit time τ of optical fiber Sagnac interferometer calculates by following formula:
τ=Δn
m/f
s (2)
From above formula, the transit time τ measuring accuracy of optical fiber Sagnac interferometer depends on sample rate f
s, improve sample rate f
scan promote the measuring accuracy of transit time.
As shown in Figure 2, be that two groups of sample sequences are carried out relevant treatment and calculate the algorithm flow chart of transit time, in P1, Δ n is the delay point of two sequences, its span is 0-N
dinteger, N wherein
dmust meet N
d/ f
s> τ
estimate; In P2, n is the sequence number of two groups of sample sequences, and its span is 0-N
s, N wherein
s=tf
s, t is the analog to digital conversion time; S in P3
2(n) S
1(n+ Δ n) is that two groups of sample sequences exist under the condition of Δ n time delay at n product constantly; P4 is accumulated value R, cumulative n=0, and 1 ..., N
sn altogether
s+ 1 S
2(n) S
1(n+ Δ n) product value; P5 is for to judge whether n equals N
sif finish current circulation, if not continue circulation until n=N
s; P6 is by the process of accumulated value R insetion sequence, the corresponding R value of different Δ n; Whether P7 equals N for inquiring about Δ n
dif finish current circulation, if not continue circulation until Δ n=N
d; P8 is an ergodic algorithm, finds out the maximal value R in A (R)
m, according to sequence number, can know the time point Δ n that it is corresponding
m; Δ n in P9
m/ f
sthe value calculating is the transit time of optical fiber Sagnac interferometer.
As shown in Figure 3, be the result curve figure after two groups of sample sequences and relevant treatment.As we know from the figure, when two groups of sample sequences carry out after relevant treatment, the normalized correlation coefficient curve map obtaining has maximal value at the 50th time point, and this time point is just in time the delay point of sequence 2 and sequence 1, divided by sample rate f
sbe the time delay of two groups of signals, be exactly the transit time of optical fiber Sagnac interferometer in test macro.
Press the test system building of schematic diagram shown in Fig. 1, open the power supply of other system module before light source powers up, two detectors are normally worked, and analog to digital conversion is normally worked; Then light source powers up, operating characteristic due to light source self, in the process need regular hour powering up in a flash to stable output luminous power, be generally hundreds of millisecond between the several seconds, within this time period, there is certain funtcional relationship in its Output optical power and time, this has just met the requirement that changes source device output optical power, and this is the simplest method; The light intensity signal that first photodetector detects has been propagated one week around fiber optic loop, so the light intensity signal that it detects with second photodetector is compared, has a fixing time delay, and this time delay is the transit time of optical fiber Sagnac interferometer; Light source powers up front analog to digital converter and starts working, can collect light source powers up to whole signals of detector output in the whole process of stable output, two groups of sample sequences that analog to digital conversion is obtained, by Computer Design relevant treatment algorithm, calculate the transit time of optical fiber Sagnac interferometer.
In conjunction with formula (1), different Δ n values is to having a specific relevant accumulated value R (Δ n), [0, N in transit time environs
d] relevant accumulated value R (Δ n) corresponding to all delayed data point Δ n form an A (R), has a maximum value R in this sequence
m, the Δ n of its correspondence
mbe two groups of the strongest time points of sample sequence correlativity, Δ n
m/ f
snamely time delay, the i.e. transit time of optical fiber Sagnac interferometer.In order to improve the measuring accuracy of transit time, can improve sample rate f
s, work as sample rate f
swhile being set as 1GHz, transit time measuring accuracy corresponding to the present invention can reach 1ns.
Claims (4)
1. an optical fiber Sagnac interferometer transit time precision measurement method, comprise optical fiber Sagnac interferometer, the coupling mechanism of optical fiber Sagnac interferometer has four tail optical fibers, article one, tail optical fiber is connected with light source, second tail optical fiber is connected with integrated phase modulator, and the 3rd tail optical fiber is connected with a photodetector; It is characterized in that, the step of the method is as follows:
1) the 4th tail optical fiber at optical fiber Sagnac interferometer coupling mechanism connects another photodetector, and the phase modulation that integrated phase modulator is introduced is poor fixing, changes source device output optical power, and two photodetectors convert light intensity signal to electric signal;
2) to the output signal of two photodetectors in step 1), use the analog to digital converter based on same clock to sample, obtain two groups of sample sequences;
3) to step 2) in two groups of sample sequences carry out relevant treatment, obtain two groups of sample sequence correlativitys corresponding delay point when maximum, according to sampling rate, calculate the transit time of optical fiber Sagnac interferometer.
2. a kind of optical fiber Sagnac interferometer transit time precision measurement method according to claim 1, it is characterized in that: the implementation of described change source device output optical power is that light source powers up moment to stablize the process of Output optical power or the light source output process with additional modulation signal variation luminous power.
3. a kind of optical fiber Sagnac interferometer transit time precision measurement method according to claim 1, is characterized in that: described based on same clock, sampling rate is f
stwo analog to digital converters respectively photodetector output electrical signals is separately carried out to analog to digital conversion, obtain two groups of sample sequence S
1and S (n)
2(n).
4. a kind of optical fiber Sagnac interferometer transit time precision measurement method according to claim 3, is characterized in that: to two groups of described sample sequence S
1and S (n)
2(n) carry out relevant treatment, formula is as follows:
Wherein, N
sfor sample sequence length, Δ n is the delay point of two groups of sample sequences, and span is [0, N
d], N
d/ f
s> τ
estimate, τ wherein
estimatefor the estimated value of optical fiber Sagnac interferometer transit time, obtain the sequence that corresponding accumulated value R (Δ n) forms
the maximal value R of this sequence
mrepresent, the delay point that its corresponding Δ n is two groups of sample sequence correlativity maximums, is designated as Δ n
m, now the transit time τ of optical fiber Sagnac interferometer calculates by following formula:
τ=Δn
m/f
s (2)
From above formula, the transit time τ measuring accuracy of optical fiber Sagnac interferometer depends on sample rate f
s, improve sample rate f
scan promote the measuring accuracy of transit time.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210047158.5A CN102607807B (en) | 2012-02-28 | 2012-02-28 | Method for precisely measuring transit time of optical fiber Sagnac interferometer |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210047158.5A CN102607807B (en) | 2012-02-28 | 2012-02-28 | Method for precisely measuring transit time of optical fiber Sagnac interferometer |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102607807A CN102607807A (en) | 2012-07-25 |
CN102607807B true CN102607807B (en) | 2014-09-03 |
Family
ID=46525375
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210047158.5A Expired - Fee Related CN102607807B (en) | 2012-02-28 | 2012-02-28 | Method for precisely measuring transit time of optical fiber Sagnac interferometer |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102607807B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5567933A (en) * | 1995-02-14 | 1996-10-22 | Mason & Hanger National, Inc. | Optical fiber detection system with disturbance and positive cut-loop detection capabilities |
JP2005241431A (en) * | 2004-02-26 | 2005-09-08 | Tokyo Gas Co Ltd | Optical fiber interference type sensor |
CN101813238A (en) * | 2010-04-20 | 2010-08-25 | 上海大学 | Sagnac/Mach-Zehnder interferometer profile fiber sensing system and time domain positioning method thereof |
CN102132129A (en) * | 2008-08-28 | 2011-07-20 | 诺思罗普·格鲁曼·利特夫有限责任公司 | Fiber optic interferometer and method for determining physical state parameters in the interior of a fiber coil of a fiber optic interferometer |
-
2012
- 2012-02-28 CN CN201210047158.5A patent/CN102607807B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5567933A (en) * | 1995-02-14 | 1996-10-22 | Mason & Hanger National, Inc. | Optical fiber detection system with disturbance and positive cut-loop detection capabilities |
JP2005241431A (en) * | 2004-02-26 | 2005-09-08 | Tokyo Gas Co Ltd | Optical fiber interference type sensor |
CN102132129A (en) * | 2008-08-28 | 2011-07-20 | 诺思罗普·格鲁曼·利特夫有限责任公司 | Fiber optic interferometer and method for determining physical state parameters in the interior of a fiber coil of a fiber optic interferometer |
CN101813238A (en) * | 2010-04-20 | 2010-08-25 | 上海大学 | Sagnac/Mach-Zehnder interferometer profile fiber sensing system and time domain positioning method thereof |
Also Published As
Publication number | Publication date |
---|---|
CN102607807A (en) | 2012-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101484826B (en) | Method and apparatus for optoelectronic contactless range finding using the transit time principle | |
CN103926457A (en) | Method for improving closed loop feedback coefficient stability of optical fiber current transformer | |
CN109547098B (en) | Microwave photon time delay measurement calibrating device | |
CN104677396A (en) | Dynamic distributed Brillouin optical fiber sensing device and method | |
CN102914423B (en) | Measuring method for sag frequency of dispersion optical fiber | |
CN106656491B (en) | The compensation method of half-wave voltage of phase modulator in a kind of quantum key dispatching system | |
CN102798750A (en) | Method and system for measuring half-wave voltage of electro-optical modulator | |
CN102768303B (en) | Method for measuring initial phase difference of all-fiber current transformer | |
CN104914444A (en) | Long-distance laser heterodyne interference range-finding structure | |
CN101881669A (en) | Real-time closed-loop measuring and tracking method of half-wave voltage of integrated electro-optical phase modulator | |
CN102853920A (en) | High-sensitivity Mach-Zahnder interference system | |
CN101968507A (en) | Optical fiber voltage sensor and adjustment method thereof | |
Li et al. | Signal detection for optical AC and DC voltage sensors based on Pockels effect | |
CN102353452A (en) | System for measuring free spectral range of F-P (Fabry-Perot) cavity | |
CN113091617B (en) | Multimode optical fiber optical path change measuring system | |
CN103698571A (en) | Current transformer equipment with effects of self energy supply and low consumption and bus current detecting method | |
CN109831249B (en) | Device for measuring differential time delay of polarization maintaining optical fiber spindle | |
CN102901515A (en) | Rapid on-line measuring method for fiber-optic gyroscope transition time | |
CN102928647B (en) | Optical profile type voltage sensor system and corresponding iterative demodulation method | |
CN101738167A (en) | Resonant cavity frequency stabilization-based absolute distance measurement system and implementing method thereof | |
CN102607807B (en) | Method for precisely measuring transit time of optical fiber Sagnac interferometer | |
CN103344414B (en) | PZT modulation factor proving installation and method of testing | |
CN101592526A (en) | A kind of measuring method of average light wavelength and device | |
CN209979850U (en) | Large-range and high-precision absolute distance measuring instrument based on OEO (optical output interface) quick switching | |
CN103940593A (en) | Quick measuring method for fiber-optic gyroscope transition time based on digital sawtooth wave |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140903 Termination date: 20190228 |