CN108332947A - A kind of fiber grating refractive index sensing demodulation method based on cutoff wavelength - Google Patents
A kind of fiber grating refractive index sensing demodulation method based on cutoff wavelength Download PDFInfo
- Publication number
- CN108332947A CN108332947A CN201810041633.5A CN201810041633A CN108332947A CN 108332947 A CN108332947 A CN 108332947A CN 201810041633 A CN201810041633 A CN 201810041633A CN 108332947 A CN108332947 A CN 108332947A
- Authority
- CN
- China
- Prior art keywords
- refractive index
- wavelength
- spectrum
- cut
- fiber grating
- 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.)
- Granted
Links
- 239000000835 fiber Substances 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000001228 spectrum Methods 0.000 claims abstract description 76
- 238000005253 cladding Methods 0.000 claims abstract description 64
- 238000000411 transmission spectrum Methods 0.000 claims abstract description 45
- 230000005540 biological transmission Effects 0.000 claims abstract description 14
- 238000005259 measurement Methods 0.000 claims abstract description 14
- 230000007423 decrease Effects 0.000 claims description 7
- 230000003595 spectral effect Effects 0.000 claims description 6
- 239000013307 optical fiber Substances 0.000 abstract description 3
- 238000004590 computer program Methods 0.000 abstract description 2
- 238000002474 experimental method Methods 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 4
- 239000002994 raw material Substances 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
- G01M11/0228—Testing optical properties by measuring refractive power
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
- G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
- G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
- G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
- G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
- G01D5/35338—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using other arrangements than interferometer arrangements
- G01D5/35354—Sensor working in reflection
- G01D5/35367—Sensor working in reflection using reflected light other than backscattered to detect the measured quantity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
- G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
- G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
- G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
- G01D5/36—Forming the light into pulses
- G01D5/38—Forming the light into pulses by diffraction gratings
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于截止波长的光纤光栅折射率传感解调方法,该方法以倾斜光纤光栅所处环境折射率大于或者等于包层折射率时的透射光谱作为参考光谱,将不同折射率下倾斜光纤光栅的透射光谱即传感测量光谱减去参考光谱得到解调光谱,以解调光谱中相对透射强度等于某一设定值时的最短波长作为测量标准,根据该波长与折射率函数关系解调出折射率。该方法相对于通过包层模对不同折射率的偏移量进行折射率传感解调方法,以及通过包层模的截止模式对不同折射率的偏移量进行折射率传感解调方法,能提供更高精度的解调;且易于用计算机程序和算法自动进行的解调方法。本发明方法可应用于光纤光栅折射率传感领域。
The invention discloses a fiber grating refractive index sensing and demodulation method based on the cut-off wavelength. The method uses the transmission spectrum when the refractive index of the environment where the inclined fiber grating is located is greater than or equal to the cladding refractive index as a reference spectrum, and uses different refractive index The transmission spectrum of the down-tilted fiber grating is the sensing measurement spectrum minus the reference spectrum to obtain the demodulation spectrum. The shortest wavelength when the relative transmission intensity in the demodulation spectrum is equal to a certain set value is used as the measurement standard. According to the wavelength and the refractive index function The relation demodulates the refractive index. Compared with the method of performing refractive index sensing and demodulating the offsets of different refractive indices through the cladding mode and the method of performing refractive index sensing and demodulating the offsets of different refractive indices through the cut-off mode of the cladding mode, the method, It can provide higher precision demodulation; and it is easy to automatically perform demodulation method with computer programs and algorithms. The method of the invention can be applied to the field of optical fiber grating refractive index sensing.
Description
技术领域technical field
本发明属于光纤光栅传感技术领域,涉及一种适用于倾斜光纤光栅的折射率传感的解调方法。The invention belongs to the technical field of optical fiber grating sensing, and relates to a demodulation method suitable for sensing the refractive index of an inclined optical fiber grating.
背景技术Background technique
在光纤光栅折射率传感中,目前主要采用通过测量倾斜光纤光栅透射谱中的包层模波长在不同折射率的偏移和测量包层模的截止模式波长在不同折射率的偏移这两种方法进行折射率传感解调。通过包层模对不同折射率的偏移量不同的解调方法主要是,先测定倾斜光纤光栅透射谱中的包层模波长在不同折射率环境下的偏移量,进一步根据偏移量与折射率之间的函数关系解调出待测介质的折射率。通过包层模的截止模式对不同折射率的偏移量不同的解调方法主要是,先测定倾斜光纤光栅透射谱中的包层模的截止模式波长在不同折射率环境下的偏移量,进一步根据偏移量与折射率之间的函数关系解调出待测介质的折射率。In fiber Bragg grating refractive index sensing, at present, it is mainly used to measure the shift of the cladding mode wavelength at different refractive indices in the transmission spectrum of the tilted fiber grating and measure the shift of the cutoff mode wavelength of the cladding mode at different refractive indices. A method for refractive index sensing demodulation. The method of demodulating different offsets of different refractive indices through the cladding mode is mainly to first measure the offset of the cladding mode wavelength in the transmission spectrum of the tilted fiber grating under different refractive index environments, and further according to the offset and The functional relationship between the refractive indices is demodulated to obtain the refractive index of the medium to be measured. The method of demodulating different offsets of different refractive indices by the cut-off mode of the cladding mode is mainly to first measure the offset of the cut-off mode wavelength of the cladding mode in the transmission spectrum of the inclined fiber grating under different refractive index environments, The refractive index of the medium to be measured is further demodulated according to the functional relationship between the offset and the refractive index.
通过包层模对不同折射率的偏移量进行折射率传感解调实现简单,但灵敏度低(小于30nm/RIU),分辨率低(在光谱仪分辨率为0.02nm条件下,分辨率为10-4RIU量级)。通过包层模的截止模式的解调方法折射率灵敏度高,可达到574nm/RIU。但是由于包层模式是离散的,在折射率引起的包层截止模式波长变化处于两个包层模式中间时很难测量,因此所能测量的折射率变化分辨率受限制。此外对不同折射率的偏移量进行折射率传感时,没有给出判定截止模式的具体标准,在判定截止模式的时候存在误差。Refractive index sensing and demodulation of different refractive index offsets through the cladding mode is simple to implement, but the sensitivity is low (less than 30nm/RIU) and the resolution is low (under the condition of a spectrometer resolution of 0.02nm, the resolution is 10 -4 RIU magnitude). The demodulation method through the cut-off mode of the cladding mode has high refractive index sensitivity, which can reach 574nm/RIU. However, since the cladding modes are discrete, it is difficult to measure when the wavelength change of the cladding cut-off mode caused by the refractive index is in the middle of the two cladding modes, so the resolution of the refractive index change that can be measured is limited. In addition, when the refractive index sensing is performed on the offsets of different refractive indices, no specific criteria for judging the cut-off mode are given, and there are errors in judging the cut-off mode.
发明内容Contents of the invention
本发明的目的在于克服上述现有技术的缺点,提供一种基于截止波长的光纤光栅折射率传感解调方法,其采用折射率传感时倾斜光纤光栅的透射谱减去包层模被压制的光滑透射谱得到差谱,作为解调光谱,在解调光谱中设定阈值用于判定截止波长,再根据截止波长与介质折射率的函数关系解调出折射率。该方法不仅解调精度高,而且截止波长可以很方便的通过计算机或者其他算法得到。The purpose of the present invention is to overcome the shortcoming of above-mentioned prior art, provide a kind of fiber grating refractive index sensing demodulation method based on cut-off wavelength, when it adopts refractive index sensing, the transmission spectrum of tilted fiber grating minus the cladding mode is suppressed The difference spectrum is obtained from the smooth transmission spectrum of the smooth transmission spectrum, which is used as the demodulation spectrum, and the threshold value is set in the demodulation spectrum to determine the cut-off wavelength, and then the refractive index is demodulated according to the functional relationship between the cut-off wavelength and the refractive index of the medium. This method not only has high demodulation accuracy, but also the cut-off wavelength can be easily obtained by computer or other algorithms.
本发明的目的是通过以下技术方案来解决的:The purpose of the present invention is solved by the following technical solutions:
本发明提供的一种基于截止波长的光纤光栅折射率传感解调方法,包括以下步骤:A kind of fiber grating refractive index sensing demodulation method based on the cut-off wavelength provided by the present invention comprises the following steps:
(1)测量外界环境折射率大于或者等于倾斜光纤光栅的包层的折射率时的倾斜光纤光栅透射光谱,将其作为参考光谱;(1) measure the transmission spectrum of the tilted fiber grating when the refractive index of the external environment is greater than or equal to the refractive index of the cladding of the tilted fiber grating, and use it as a reference spectrum;
(2)在需要进行传感解调的折射率变化的环境中,测量倾斜光纤光栅的透射光谱,将其作为传感测量光谱;(2) Measure the transmission spectrum of the tilted fiber grating in the environment where the refractive index changes that require sensing and demodulation, and use it as the sensing measurement spectrum;
(3)用步骤(2)中得到的传感测量光谱减去步骤(1)中得到的参考光谱,得到解调光谱;(3) subtract the reference spectrum obtained in step (1) from the sensing measurement spectrum obtained in step (2), to obtain the demodulated spectrum;
(4)根据解调光谱中的包层模式透射光谱强度,设定一个阈值,从短波长至长波长开始计算,第一个达到设定阈值的波长记作截止波长;(4) According to the cladding mode transmission spectrum intensity in the demodulation spectrum, set a threshold value, start to calculate from short wavelength to long wavelength, and the first wavelength that reaches the set threshold value is recorded as the cut-off wavelength;
(5)测量出截止波长与折射率之间的函数曲线,对倾斜光纤光栅截止波长在不同折射率下响应系数进行标定;(5) Measure the function curve between the cut-off wavelength and the refractive index, and calibrate the response coefficient of the cut-off wavelength of the inclined fiber Bragg grating under different refractive indices;
(6)根据步骤(5)中测量得到的截止波长与折射率之间的函数曲线关系以及待测环境下的透射光谱差谱,对需要进行传感解调的折射率变化的环境的折射率进行解调。(6) According to the function curve relationship between the cut-off wavelength and the refractive index measured in step (5) and the transmission spectrum difference spectrum under the environment to be measured, the refractive index of the environment where the refractive index change needs to be sensed and demodulated to demodulate.
上述技术方案中,还有进一步限定的方案。Among the above technical solutions, there are further limited solutions.
进一步的,步骤(1)中,所述外界环境折射率大于或等于包层的折射率,指的是环境中的折射率大于等于倾斜光纤光栅的包层的折射率,使得倾斜光纤光栅中的包层模式响应产生的光谱调制消失,透射谱中只有包层模式引起缓变透射损耗。Further, in step (1), the refractive index of the external environment is greater than or equal to the refractive index of the cladding, which means that the refractive index in the environment is greater than or equal to the refractive index of the cladding of the tilted fiber grating, so that in the tilted fiber grating The spectral modulation produced by the cladding mode response disappears, and only the cladding mode causes the slowly-varying transmission loss in the transmission spectrum.
进一步的,步骤(4)中,所述的包层模式透射光谱调制强度随环境折射率的增大自短波长开始而逐渐减小,调制强度开始发生减小的波长位置随环境折射率的增大而向长波长方向偏移。Further, in step (4), the modulation intensity of the cladding mode transmission spectrum gradually decreases from short wavelengths with the increase of the environmental refractive index, and the wavelength position at which the modulation intensity begins to decrease decreases with the increase of the environmental refractive index. Large and shifted to the long wavelength direction.
进一步的,步骤(4)中,所述的阈值k可以选取为n<k<m,其中m为解调光谱中包层模式强度的最大值,n为解调光谱中包层模式强度的最小值。Further, in step (4), the threshold k can be selected as n<k<m, where m is the maximum value of the cladding mode intensity in the demodulated spectrum, and n is the minimum cladding mode intensity in the demodulated spectrum value.
进一步的,步骤(4)中,所述的截止波长采用计算机确定:首先设定误差范围,再对解调光谱进行插值,使得解调光谱的离散间隔小于误差值,从短波长至长波长方向,使得解调光谱中包层模式强度减去阈值k的差值小于误差值,该强度对应的波长即为截止波长。Further, in step (4), the cut-off wavelength is determined by computer: first set the error range, and then interpolate the demodulated spectrum, so that the discrete interval of the demodulated spectrum is smaller than the error value, from the short wavelength to the long wavelength direction , so that the difference between the intensity of the cladding mode in the demodulated spectrum minus the threshold k is less than the error value, and the wavelength corresponding to the intensity is the cut-off wavelength.
进一步的,所述倾斜光纤光栅是基于飞秒激光制备的Further, the tilted fiber grating is prepared based on femtosecond laser
与现有技术相比,本发明具有以下有益效果:Compared with the prior art, the present invention has the following beneficial effects:
本发明提供的基于截止波长的光纤光栅折射率传感解调方法,将倾斜光纤光栅所处环境折射率大于或者等于包层折射率时的透射光谱作为参考光谱,解调时用传感测量光谱减去倾斜光纤光栅的参考光谱得到差谱,作为解调光谱,在解调光谱中设定阈值用于判定截止波长,解调光谱中相对透射强度等于阈值时的最短波长记作截止波长,再根据截止波长与介质折射率的函数关系解调出折射率。该方法与现有技术相比,截止波长的识别简单,准确,易于实现计算机自动识别,折射率传感的解调精度高。The cut-off wavelength-based fiber grating refractive index sensing and demodulation method provided by the present invention uses the transmission spectrum when the refractive index of the environment where the tilted fiber grating is located is greater than or equal to the cladding refractive index as the reference spectrum, and uses the sensor to measure the spectrum during demodulation Subtract the reference spectrum of the tilted fiber grating to obtain the difference spectrum, as the demodulation spectrum, set the threshold in the demodulation spectrum to determine the cut-off wavelength, the shortest wavelength when the relative transmission intensity in the demodulation spectrum is equal to the threshold is recorded as the cut-off wavelength, and then The refractive index is demodulated according to the functional relationship between the cut-off wavelength and the refractive index of the medium. Compared with the prior art, the method has the advantages of simple and accurate recognition of the cut-off wavelength, easy realization of computer automatic recognition, and high demodulation precision of the refractive index sensor.
进一步的,本发明解决了通过包层模对不同折射率的偏移量进行折射率传感解调灵敏度低(小于30nm/RIU),分辨率低(在光谱仪分辨率为0.02nm条件下,分辨率为10-4RIU量级)的问题;本发明采用的基于截止波长的光纤光栅折射率传感解调方法能够实现精度高达949.9nm/RIU。Further, the present invention solves the problem of low sensitivity (less than 30nm/RIU) and low resolution (under the condition that the resolution of the spectrometer is 0.02nm, the resolution of the refractive index sensing and demodulation of different refractive index offsets is low by using the cladding mode). rate of 10 -4 RIU level); the cut-off wavelength-based fiber grating refractive index sensing and demodulation method used in the present invention can achieve a precision as high as 949.9nm/RIU.
进一步的,本发明相对于现有技术通过包层模的截止模式的解调方法灵敏度达574nm/RIU,包层模式是离散的,在折射率引起的包层截止模式波长变化处于两个包层模式中间时很难利用该方法测量,因此所能测量的折射率变化分辨率受限制的问题;此外由于包层模式不仅调制度发生变化,包层模波长位置也发生偏移,截止模式的判别算法比较复杂。本发明的解调方法,在对不同折射率的偏移量进行折射率传感时,通过设定阈值的方法能够通过计算机程序或算法对截止波长进行自动识别,判定方法简单,容易进行定量判定。Further, compared with the prior art, the sensitivity of the demodulation method of the cut-off mode of the cladding mode is up to 574nm/RIU. The cladding mode is discrete, and the wavelength change of the cladding cut-off mode caused by the refractive index is in two cladding layers. It is difficult to use this method to measure in the middle of the mode, so the resolution of the refractive index change that can be measured is limited; in addition, because the cladding mode not only changes the modulation degree, but also shifts the wavelength position of the cladding mode, the discrimination of the cut-off mode The algorithm is more complicated. The demodulation method of the present invention can automatically identify the cut-off wavelength through a computer program or algorithm by setting a threshold when performing refractive index sensing on the offsets of different refractive indices, and the determination method is simple and easy to perform quantitative determination .
附图说明Description of drawings
图1是环境折射率分别为1.4578时,倾斜光纤光栅的透射光谱;Figure 1 is the transmission spectrum of the tilted fiber grating when the ambient refractive index is 1.4578;
图2是环境折射率为1.4401时的倾斜光纤光栅透射光谱图;Fig. 2 is the transmission spectrum diagram of the tilted fiber grating when the refractive index of the environment is 1.4401;
图3是选取阈值为0.1dB时的解调光谱图;Fig. 3 is the demodulation spectrogram when selecting threshold to be 0.1dB;
图4是选取阈值为0.1dB时,环境折射率和截止波长的对应关系;Figure 4 is the corresponding relationship between the ambient refractive index and the cut-off wavelength when the threshold is selected as 0.1dB;
图5是选取阈值为0.3dB时的解调光谱图;Fig. 5 is the demodulation spectrogram when selecting threshold to be 0.3dB;
图6是选取阈值为0.3dB时,环境折射率和截止波长的对应关系;Figure 6 is the corresponding relationship between the ambient refractive index and the cut-off wavelength when the threshold is selected as 0.3dB;
图7是选取阈值为0.9dB时的解调光谱图;Fig. 7 is the demodulation spectrogram when selecting threshold to be 0.9dB;
图8是选取阈值为0.9dB时,环境折射率和截止波长的对应关系。Figure 8 shows the corresponding relationship between the ambient refractive index and the cut-off wavelength when the threshold is selected as 0.9dB.
具体实施方式Detailed ways
下面结合附图和实施例对本发明作进一步描述。The present invention will be further described below in conjunction with the accompanying drawings and embodiments.
为了更加清晰说明本发明的目的、技术方案及优点,以下结合附图及实施例对本发明作进一步详细说明。应当理解,此处所述的具体实施例仅仅用以解释本发明,并不用于限定本发明。In order to more clearly illustrate the purpose, technical solutions and advantages of the present invention, the present invention will be further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described here are only used to explain the present invention, not to limit the present invention.
本发明的基于截止波长的光纤光栅折射率传感解调方法,具体按照以下步骤进行:The fiber grating refractive index sensing demodulation method based on the cut-off wavelength of the present invention is specifically carried out according to the following steps:
(1)将进行折射率传感的倾斜光纤光栅折射率传感器置于折射率可变的环境中,并同时测量倾斜光纤光栅的透射谱,通过改变环境折射率大于或等于光纤包层折射率,使得倾斜光纤光栅中的包层模式响应产生的光谱调制消失,透射谱中只有包层模式引起缓变透射损耗,透射谱显示为一条光滑的曲线,保存此时的透射光谱作为参考光谱。(1) Place the tilted fiber Bragg grating refractive index sensor for refractive index sensing in an environment with variable refractive index, and measure the transmission spectrum of the tilted fiber Bragg grating at the same time, by changing the refractive index of the environment to be greater than or equal to the refractive index of the fiber cladding, The spectral modulation caused by the response of the cladding mode in the tilted fiber grating disappears, and only the cladding mode in the transmission spectrum causes the slow-varying transmission loss. The transmission spectrum is displayed as a smooth curve, and the transmission spectrum at this time is saved as a reference spectrum.
上述倾斜光纤光栅是基于飞秒激光制备的。The above-mentioned tilted fiber grating is prepared based on femtosecond laser.
(2)在需要进行传感解调的待测介质中,测量倾斜光纤光栅的透射光谱,将其作为传感测量光谱。(2) In the medium to be measured that needs to be sensed and demodulated, measure the transmission spectrum of the tilted fiber grating, and use it as the sensing measurement spectrum.
(3)解调时,用步骤(2)中得到传感测量光谱减去步骤(1)中保存的参考光谱,得到解调光谱。(3) During demodulation, subtract the reference spectrum saved in step (1) from the sensing measurement spectrum obtained in step (2) to obtain the demodulation spectrum.
(4)根据解调光谱中的包层模的幅值,设定一个阈值k,阈值k可以选取为n<k<m,其中m为解调光谱中包层模式强度的最大值,n为解调光谱中包层模式强度的最小值,例如可以是0.1dB,0.3dB,0.9dB等值。包层模式透射光谱调制强度随环境折射率的增大自短波长开始而逐渐减小,调制强度开始发生减小的波长位置随环境折射率的增大而向长波长方向偏移。(4) According to the amplitude of the cladding mode in the demodulated spectrum, a threshold k is set, and the threshold k can be selected as n<k<m, where m is the maximum value of the cladding mode intensity in the demodulated spectrum, and n is The minimum value of the cladding mode intensity in the demodulated spectrum may be, for example, 0.1dB, 0.3dB, 0.9dB or the like. The modulation intensity of the cladding mode transmission spectrum decreases gradually with the increase of the ambient refractive index from the short wavelength, and the wavelength position where the modulation intensity begins to decrease shifts to the long wavelength direction with the increase of the ambient refractive index.
以阈值的大小为标准画一条虚线,从短波长至长波长方向,第一个与虚线相交的包层模对应的波长记作截止波长,不同的截止波长对应着环境中不同的折射率,可以进行折射率解调。Draw a dotted line based on the threshold value, from the short wavelength to the long wavelength direction, the first wavelength corresponding to the cladding mode that intersects with the dotted line is recorded as the cut-off wavelength, and different cut-off wavelengths correspond to different refractive indices in the environment, which can be Perform refractive index demodulation.
截止波长采用计算机或手动方式确定:首先设定误差范围,再对解调光谱进行插值,使得解调光谱的离散间隔小于误差值,从短波长至长波长方向,使得解调光谱中包层模式强度减去阈值k的差值小于误差值,该强度对应的波长即为截止波长。The cut-off wavelength is determined by computer or manually: first set the error range, and then interpolate the demodulated spectrum, so that the discrete interval of the demodulated spectrum is smaller than the error value, from the short wavelength to the long wavelength direction, so that the cladding mode in the demodulated spectrum If the difference between the intensity minus the threshold k is less than the error value, the wavelength corresponding to the intensity is the cut-off wavelength.
(5)测量出截止波长与折射率之间的函数曲线,对倾斜光纤光栅截止波长在不同折射率下响应系数进行标定。(5) Measure the function curve between the cut-off wavelength and the refractive index, and calibrate the response coefficient of the cut-off wavelength of the inclined fiber Bragg grating under different refractive indices.
(6)根据步骤(5)中测量得到的截止波长与折射率之间的函数曲线关系以及待测环境下的透射光谱差谱,对需要进行传感解调的折射率变化的环境的折射率进行解调。(6) According to the function curve relationship between the cut-off wavelength and the refractive index measured in step (5) and the transmission spectrum difference spectrum under the environment to be measured, the refractive index of the environment where the refractive index change needs to be sensed and demodulated to demodulate.
下面通过几个具体实施例来对本发明做进一步详细描述:The present invention will be described in further detail below by several specific embodiments:
本发明适用于各种倾斜光纤光栅折射率传感解调,本实施例选用飞秒激光加工的倾斜角度为3.8°倾斜光纤光栅进行实验,在解调光谱中以不同阈值进行截止波长的选取。The present invention is applicable to sensing demodulation of various inclined fiber grating refractive index. In this embodiment, the inclined angle of femtosecond laser processing is selected as 3.8° inclined fiber grating for experiments, and the cut-off wavelength is selected with different thresholds in the demodulation spectrum.
实施例1Example 1
本实施例以倾斜角度为3.8°倾斜光纤光栅为例进行实验,在解调光谱中以0.1dB为阈值进行截止波长的选取,具体如下:In this embodiment, an experiment is carried out with an inclined fiber grating whose inclination angle is 3.8° as an example, and the cut-off wavelength is selected with 0.1 dB as the threshold in the demodulation spectrum, as follows:
原始材料:飞秒激光制备的倾斜角度为3.8°倾斜光纤光栅。Raw materials: Femtosecond laser fabricated tilted fiber gratings with a tilt angle of 3.8°.
折射率传感解调步骤详细实施方式阐述如下:The detailed implementation of the refractive index sensing demodulation step is described as follows:
(1)将进行折射率传感的倾斜光纤光栅折射率传感器置于折射率可变的环境中,并同时测量倾斜光纤光栅的透射谱,通过改变环境折射率到达1.4578,使得倾斜光纤光栅中的包层模式响应产生的光谱调制消失,透射谱中只有包层模式引起缓变透射损耗,透射谱显示为一条光滑的曲线,见图1,保存该透射光谱作为参考光谱;(1) Place the tilted fiber Bragg grating refractive index sensor for refractive index sensing in an environment with variable refractive index, and measure the transmission spectrum of the tilted fiber Bragg grating at the same time. By changing the refractive index of the environment to 1.4578, the The spectral modulation produced by the cladding mode response disappears, and only the cladding mode in the transmission spectrum causes the slow-varying transmission loss, and the transmission spectrum is displayed as a smooth curve, as shown in Figure 1, and the transmission spectrum is saved as a reference spectrum;
(2)当改变环境折射率为1.4401时,测量倾斜光纤光栅折射率传感器的透射光谱,将其作为传感测量光谱,见图2;(2) When the refractive index of the environment is changed to 1.4401, the transmission spectrum of the tilted fiber grating refractive index sensor is measured, and it is used as the sensing measurement spectrum, as shown in Figure 2;
(3)用步骤(2)中得到传感测量光谱(图2)减去步骤(1)中保存的参考光谱(图1)得到的差谱,记作解调光谱,见图3;(3) subtract the difference spectrum (Fig. 1) that obtains from the reference spectrum (Fig. 1) preserved in the step (1) from the sensing measurement spectrum (Fig. 2) obtained in the step (2), and record it as the demodulation spectrum, see Fig. 3;
(4)根据解调光谱(图3)的包层模相对透射强度,设定一个阈值k为0.1dB,见图3中的1,阈值的大小在包层模最大幅值和最小幅值之间,以阈值的大小为标准画一条虚线,从短波长至长波长方向,第一个达到设定阈值的透射强度对应波长记作该折射率下的截止波长,见图3中的2,该截止波长即对应折射率为1.4401的环境;(4) According to the relative transmission intensity of the cladding mode of the demodulation spectrum (Fig. 3), setting a threshold k is 0.1dB, see 1 in Fig. 3, the size of the threshold is between the cladding mode maximum amplitude and the minimum amplitude , draw a dotted line with the threshold value as the standard, from the short wavelength to the long wavelength direction, the wavelength corresponding to the first transmission intensity that reaches the set threshold value is recorded as the cut-off wavelength under the refractive index, see 2 in Figure 3, the The cut-off wavelength corresponds to the environment with a refractive index of 1.4401;
(5)分别改变介质折射率为1.4125RIU,1.416RIU,1.4205RIU,1.4251RIU,1.4324RIU,1.4346RIU,1.4392RIU,1.4435RIU,1.4487RIU,1.4532RIU,1.4552RIU,并测量出不同折射率下的截止波长进行标定。阈值为0.1dB下的折射率和截止波长对应关系,见图4。(5) Change the refractive index of the medium to 1.4125RIU, 1.416RIU, 1.4205RIU, 1.4251RIU, 1.4324RIU, 1.4346RIU, 1.4392RIU, 1.4435RIU, 1.4487RIU, 1.4532RIU, 1.4552RIU, and measure the The cutoff wavelength is calibrated. See Figure 4 for the relationship between the refractive index and the cut-off wavelength at a threshold of 0.1dB.
(6)根据步骤(5)中测量得到的截止波长与折射率之间的函数曲线关系为:Y=423.675-281.185X(外界环境折射率范围为1.4125~1.4435),Y=1389.557-949.878X(外界环境折射率范围为1.4435~1.4552),式中Y为截至波长的偏移量,X为待测未知环境折射率。解调精度可达到949.878nm/RIU。(6) According to the function curve relation between the cut-off wavelength measured in step (5) and the refractive index is: Y=423.675-281.185X (the external environment refractive index range is 1.4125~1.4435), Y=1389.557-949.878X ( The refractive index of the external environment ranges from 1.4435 to 1.4552), where Y is the offset of the cut-off wavelength, and X is the refractive index of the unknown environment to be measured. The demodulation accuracy can reach 949.878nm/RIU.
实施例2Example 2
本实施例以倾斜角度为3.8°倾斜光纤光栅为例进行实验,在解调光谱中以0.3dB为阈值进行截止波长的选取,具体如下:In this embodiment, an experiment is carried out with an inclined fiber grating whose inclination angle is 3.8° as an example. In the demodulation spectrum, 0.3 dB is used as the threshold to select the cut-off wavelength, as follows:
原始材料:飞秒激光制备的倾斜角度为3.8°倾斜光纤光栅。Raw materials: Femtosecond laser fabricated tilted fiber gratings with a tilt angle of 3.8°.
折射率传感解调步骤详细实施方式阐述如下:The detailed implementation of the refractive index sensing demodulation step is described as follows:
(1)将进行折射率传感的倾斜光纤光栅折射率传感器置于折射率可变的环境中,并同时测量倾斜光纤光栅的透射谱,通过改变环境折射率到达1.4578,使得倾斜光纤光栅中的包层模式响应产生的光谱调制消失,透射谱中只有包层模式引起缓变透射损耗,透射谱显示为一条光滑的曲线,见图1,保存该透射光谱作为参考光谱;(1) Place the tilted fiber Bragg grating refractive index sensor for refractive index sensing in an environment with variable refractive index, and measure the transmission spectrum of the tilted fiber Bragg grating at the same time. By changing the refractive index of the environment to 1.4578, the The spectral modulation produced by the cladding mode response disappears, and only the cladding mode in the transmission spectrum causes the slow-varying transmission loss, and the transmission spectrum is displayed as a smooth curve, as shown in Figure 1, and the transmission spectrum is saved as a reference spectrum;
(2)当改变环境折射率为1.4401时,测量倾斜光纤光栅折射率传感器的透射光谱,将其作为传感测量光谱,见图2;(2) When the refractive index of the environment is changed to 1.4401, the transmission spectrum of the tilted fiber grating refractive index sensor is measured, and it is used as the sensing measurement spectrum, as shown in Figure 2;
(3)用步骤(2)中得到传感测量光谱(图2)减去步骤(1)中保存的参考光谱(图1)得到的差谱,记作解调光谱,见图5;(3) subtract the difference spectrum (Fig. 1) that obtains from the reference spectrum (Fig. 1) preserved in the step (1) from the sensing measurement spectrum (Fig. 2) obtained in the step (2), and record it as the demodulation spectrum, see Fig. 5;
(4)根据解调光谱(图5)的包层模相对透射强度,设定一个阈值为0.3dB,见图5中的3,阈值的大小在包层模最大幅值和最小幅值之间,以阈值的大小为标准画一条虚线,从短波长至长波长方向,第一个达到设定阈值的透射强度对应波长记作该折射率下的截止波长,见图5中的4,该截止波长即对应折射率为1.4401的环境;(4) According to the relative transmission intensity of the cladding mode of the demodulation spectrum (Fig. 5), a threshold is set to be 0.3dB, see 3 in Fig. 5, the size of the threshold is between the cladding mode maximum amplitude and the minimum amplitude , draw a dotted line with the threshold value as the standard, from the short wavelength to the long wavelength direction, the wavelength corresponding to the first transmission intensity reaching the set threshold value is recorded as the cut-off wavelength under the refractive index, see 4 in Figure 5, the cut-off wavelength The wavelength corresponds to the environment with a refractive index of 1.4401;
(5)分别改变介质折射率为1.4125RIU,1.416RIU,1.4205RIU,1.4251RIU,1.4324RIU,1.4346RIU,1.4392RIU,1.4435RIU,1.4487RIU,1.4532RIU,1.4552RIU,并测量出不同折射率下的截止波长进行标定。阈值为0.3dB下的折射率和截止波长对应关系,见图6。(5) Change the refractive index of the medium to 1.4125RIU, 1.416RIU, 1.4205RIU, 1.4251RIU, 1.4324RIU, 1.4346RIU, 1.4392RIU, 1.4435RIU, 1.4487RIU, 1.4532RIU, 1.4552RIU, and measure the The cutoff wavelength is calibrated. See Figure 6 for the corresponding relationship between the refractive index and the cut-off wavelength at a threshold of 0.3dB.
(6)根据步骤(5)中测量得到的截止波长与折射率之间的函数曲线关系为:Y=479.764-321.845X(外界环境折射率范围为1.416~1.4435),Y=1374.859-941.715X(外界环境折射率范围为1.4435~1.4552),式中Y为截至波长的偏移量,X为待测未知环境折射率。解调精度可达到941.715nm/RIU。(6) According to the function curve relation between the cut-off wavelength measured in step (5) and the refractive index is: Y=479.764-321.845X (the external environment refractive index range is 1.416~1.4435), Y=1374.859-941.715X ( The refractive index of the external environment ranges from 1.4435 to 1.4552), where Y is the offset of the cut-off wavelength, and X is the refractive index of the unknown environment to be measured. The demodulation accuracy can reach 941.715nm/RIU.
实施例3Example 3
本实施例以倾斜角度为3.8°倾斜光纤光栅为例进行实验,在解调光谱中以0.9dB为阈值进行截止波长的选取,具体如下:In this embodiment, an experiment is carried out with an inclined fiber grating whose inclination angle is 3.8° as an example, and the cut-off wavelength is selected with 0.9 dB as the threshold in the demodulation spectrum, as follows:
原始材料:飞秒激光制备的倾斜角度为3.8°倾斜光纤光栅。Raw materials: Femtosecond laser fabricated tilted fiber gratings with a tilt angle of 3.8°.
折射率传感解调步骤详细实施方式阐述如下:The detailed implementation of the refractive index sensing demodulation step is described as follows:
(1)将进行折射率传感的倾斜光纤光栅折射率传感器置于折射率可变的环境中,并同时测量倾斜光纤光栅的透射谱,通过改变环境折射率到达1.4578,使得倾斜光纤光栅中的包层模式响应产生的光谱调制消失,透射谱中只有包层模式引起缓变透射损耗,透射谱显示为一条光滑的曲线,见图1,保存该透射光谱作为参考光谱;(1) Place the tilted fiber Bragg grating refractive index sensor for refractive index sensing in an environment with variable refractive index, and measure the transmission spectrum of the tilted fiber Bragg grating at the same time. By changing the refractive index of the environment to 1.4578, the The spectral modulation produced by the cladding mode response disappears, and only the cladding mode in the transmission spectrum causes the slow-varying transmission loss, and the transmission spectrum is displayed as a smooth curve, as shown in Figure 1, and the transmission spectrum is saved as a reference spectrum;
(2)当改变环境折射率为1.4401时,测量倾斜光纤光栅折射率传感器的透射光谱,将其作为传感测量光谱,见图2;(2) When the refractive index of the environment is changed to 1.4401, the transmission spectrum of the tilted fiber grating refractive index sensor is measured, and it is used as the sensing measurement spectrum, as shown in Figure 2;
(3)用步骤(2)中得到传感测量光谱(图2)减去步骤(1)中保存的参考光谱(图1)得到的差谱,记作解调光谱,见图7;(3) subtract the difference spectrum (Fig. 1) that obtains from the reference spectrum (Fig. 1) preserved in the step (1) with the sensing measurement spectrum (Fig. 2) obtained in the step (2), and record it as the demodulation spectrum, see Fig. 7;
(4)根据解调光谱(图7)的包层模相对透射强度,设定一个阈值为0.9dB,见图7中的5,阈值的大小在包层模最大幅值和最小幅值之间,以阈值的大小为标准画一条虚线,从短波长至长波长方向,第一个达到设定阈值的透射强度对应波长记作该折射率下的截止波长,见图7中的6,该截止波长即对应折射率为1.4401的环境;(4) According to the relative transmission intensity of the cladding mode of the demodulation spectrum (Fig. 7), a threshold is set to be 0.9dB, see 5 in Fig. 7, the size of the threshold is between the cladding mode maximum amplitude and the minimum amplitude , draw a dotted line with the threshold value as the standard, from the short wavelength to the long wavelength direction, the wavelength corresponding to the first transmission intensity reaching the set threshold value is recorded as the cut-off wavelength under the refractive index, see 6 in Figure 7, the cut-off wavelength The wavelength corresponds to the environment with a refractive index of 1.4401;
(5)分别改变介质折射率为1.4125RIU,1.416RIU,1.4205RIU,1.4251RIU,1.4324RIU,1.4346RIU,1.4392RIU,1.4435RIU,1.4487RIU,并测量出不同折射率下的截止波长进行标定。阈值为0.9dB下的折射率和截止波长对应关系,见图8。(5) Change the refractive index of the medium to 1.4125RIU, 1.416RIU, 1.4205RIU, 1.4251RIU, 1.4324RIU, 1.4346RIU, 1.4392RIU, 1.4435RIU, 1.4487RIU, and measure the cut-off wavelengths under different refractive indices for calibration. The corresponding relationship between the refractive index and the cut-off wavelength at the threshold of 0.9dB is shown in Figure 8.
(6)根据步骤(5)中测量得到的截止波长与折射率之间的函数曲线关系为:Y=704.062-479.631X(外界环境折射率范围为1.4251~1.4487),式中Y为截至波长的偏移量,X为待测未知环境折射率。解调精度可达到479.631nm/RIU。(6) According to the function curve relation between the cut-off wavelength measured in step (5) and the refractive index is: Y=704.062-479.631X (the external environment refractive index range is 1.4251~1.4487), Y is the cut-off wavelength in the formula Offset, X is the refractive index of the unknown environment to be measured. The demodulation accuracy can reach 479.631nm/RIU.
从上述实施例可以看出,本发明采用的基于截止波长的光纤光栅折射率传感解调方法解调精度最高可达949.9nm/RIU。It can be seen from the above embodiments that the demodulation accuracy of the fiber grating refractive index sensing demodulation method based on the cut-off wavelength used in the present invention can reach up to 949.9nm/RIU.
Claims (6)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810041633.5A CN108332947B (en) | 2018-01-16 | 2018-01-16 | A Fiber Bragg Grating Refractive Index Sensing Demodulation Method Based on Cutoff Wavelength |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810041633.5A CN108332947B (en) | 2018-01-16 | 2018-01-16 | A Fiber Bragg Grating Refractive Index Sensing Demodulation Method Based on Cutoff Wavelength |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108332947A true CN108332947A (en) | 2018-07-27 |
CN108332947B CN108332947B (en) | 2019-07-19 |
Family
ID=62925181
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810041633.5A Active CN108332947B (en) | 2018-01-16 | 2018-01-16 | A Fiber Bragg Grating Refractive Index Sensing Demodulation Method Based on Cutoff Wavelength |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108332947B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5633748A (en) * | 1996-03-05 | 1997-05-27 | The United States Of America As Represented By The Secretary Of The Navy | Fiber optic Bragg grating demodulator and sensor incorporating same |
CN1305581A (en) * | 1998-04-14 | 2001-07-25 | 智能光学系统公司 | Method and system for determining wavelength of light transmitted through optical fiber |
EP1266248A1 (en) * | 2000-03-23 | 2002-12-18 | Defence Science and Technology Organisation | Method and apparatus for estimating chromatic dispersion in fibre bragg gratings |
CN201017069Y (en) * | 2007-04-09 | 2008-02-06 | 钱祥忠 | Optical fiber grating wavelength demodulating equipment |
CN101532850A (en) * | 2009-04-20 | 2009-09-16 | 浙江大学 | Method and device for sensing and demodulating Bragg fiber grating |
CN102721665A (en) * | 2012-05-16 | 2012-10-10 | 暨南大学 | Micro nano optical fiber grating refractive index sensor |
-
2018
- 2018-01-16 CN CN201810041633.5A patent/CN108332947B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5633748A (en) * | 1996-03-05 | 1997-05-27 | The United States Of America As Represented By The Secretary Of The Navy | Fiber optic Bragg grating demodulator and sensor incorporating same |
CN1305581A (en) * | 1998-04-14 | 2001-07-25 | 智能光学系统公司 | Method and system for determining wavelength of light transmitted through optical fiber |
EP1266248A1 (en) * | 2000-03-23 | 2002-12-18 | Defence Science and Technology Organisation | Method and apparatus for estimating chromatic dispersion in fibre bragg gratings |
CN201017069Y (en) * | 2007-04-09 | 2008-02-06 | 钱祥忠 | Optical fiber grating wavelength demodulating equipment |
CN101532850A (en) * | 2009-04-20 | 2009-09-16 | 浙江大学 | Method and device for sensing and demodulating Bragg fiber grating |
CN102721665A (en) * | 2012-05-16 | 2012-10-10 | 暨南大学 | Micro nano optical fiber grating refractive index sensor |
Non-Patent Citations (4)
Title |
---|
C.CAUCHETEUR等: "Demodulation Technique for Weakly Tilted Fiber Bragg Grating Refractometer", 《PHOTONICS TECHNOLOGY LETTERS》 * |
YU-CHUN LU等: "Polarization Effects in Tilted Fiber Bragg Grating Refractometers", 《JOURNAL OF LIGHTWAVE TECHNOLOGY》 * |
张卫华: "倾斜光纤光栅传感解调技术与网络的基础研究", 《万方数据知识服务平台》 * |
张卫华等: "基于倾斜光纤光栅的传感解调技术", 《纳米技术与精密工程》 * |
Also Published As
Publication number | Publication date |
---|---|
CN108332947B (en) | 2019-07-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5190330B2 (en) | Chromatic point sensor operating method and chromatic point sensor | |
US7742892B2 (en) | Calibrating an optical FMCW backscattering measurement system | |
US8928874B2 (en) | Method for identifying abnormal spectral profiles measured by a chromatic confocal range sensor | |
CN101726361A (en) | Apparatus and method for measuring optical characteristics suitable for spectrum measurement | |
CN104713660A (en) | Optical fiber sensor, manufacturing method and simultaneous measuring method for liquid temperature and refractive index | |
CN108844919A (en) | The reflection type inclined fiber grating index sensor of covering and production, measurement method | |
CN102680097B (en) | Dichroism measuring method and dichroism measurement mechanism | |
CN107850495A (en) | Temperature-measurinsubstrate substrate and temperature measurement system | |
CN108332947A (en) | A kind of fiber grating refractive index sensing demodulation method based on cutoff wavelength | |
CN118243157A (en) | Correction method and correction system for optical fiber MEMS interference spectrum | |
CN110207849B (en) | Temperature sensor based on reflective long-period fiber grating and measuring method | |
CN110530466B (en) | Intensity modulation type liquid level sensing detection method based on double-coreless optical fiber cascade connection | |
CN114894746A (en) | Refractive index amplification method, device and medium based on single FFPI vernier effect | |
CN110940370A (en) | A constant temperature and constant strain fiber grating standard device | |
US20030202176A1 (en) | Optical power meter | |
US20160146729A1 (en) | Measurement method based on an optical waveguide sensor system | |
JP5483459B2 (en) | Evaluation method for differential transimpedance amplifier, optical receiver circuit evaluation method, program and apparatus thereof | |
CN110926518A (en) | A wavelength feedback fiber grating standard device | |
JP6484125B2 (en) | Temperature measuring apparatus and temperature measuring method | |
JP7121222B1 (en) | Film thickness measuring device and film thickness measuring method | |
CN110793937B (en) | Membrane type determination method | |
WO2022244309A1 (en) | Film thickness measurement device and film thickness measurement method | |
CN111637906B (en) | Fiber grating demodulation device and method based on self-differencing coherence | |
JP3287371B2 (en) | Adjusting the amount of received light and refractometer | |
JPS6020024Y2 (en) | Optical pressure measuring device |
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 |