CN114813633A - 一种激光多气体混叠光谱解调和非线性补偿的浓度检测方法 - Google Patents
一种激光多气体混叠光谱解调和非线性补偿的浓度检测方法 Download PDFInfo
- Publication number
- CN114813633A CN114813633A CN202210380189.6A CN202210380189A CN114813633A CN 114813633 A CN114813633 A CN 114813633A CN 202210380189 A CN202210380189 A CN 202210380189A CN 114813633 A CN114813633 A CN 114813633A
- Authority
- CN
- China
- Prior art keywords
- gas
- concentration
- harmonic
- concentrations
- signal
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/39—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using tunable lasers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Mathematical Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Pathology (AREA)
- Immunology (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Biology (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Operations Research (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Chemical & Material Sciences (AREA)
- Algebra (AREA)
- Health & Medical Sciences (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Optics & Photonics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种激光多气体混叠光谱解调和非线性补偿的浓度检测方法,包括以下步骤:首先采集待测气体的二次谐波信号,基于非负最小二乘方法对其严重混叠光谱信号进行解调,可以得到多种气体测量时的浓度结果,但是由于非线性效应的产生,利用三次多项式对浓度结果进行了补偿,获得了更加精确的多种气体浓度测量结果,以煤矿中CO和CH4为例进行了说明。通过上述方式,本发明能够实现多种气体二次谐波混叠光谱的解调,并基于解调得到的结果对浓度进行非线性补偿,得到精确的气体浓度信息,满足特殊工况环境中多种气体同时在线监测的需求。
Description
技术领域
本发明涉及气体检测技术领域,特别是涉及一种激光多气体混叠光谱解调和非线性补偿的浓度检测方法。
背景技术
采用激光光谱技术进行多种气体同时测量时,谱线之间的交叉干扰问题一直以来都是困扰的难题,比如:CO和CH4气体是煤矿生产作业环境中的两种常见气体。CO气体和人体的血红蛋白具有非常强的亲和力,这会导致血红蛋白输氧的能力大大下降,从而导致中毒症状的产生。CH4气体又称瓦斯,是煤矿开采中产生的主要气体,是一种易燃易爆气体,严重威胁着煤矿的生产安全。因此对煤矿中的CO和CH4气体进行实时在线高精度检测是非常有必要的。一般情况下,矿井中允许的CO浓度不超过24ppm,允许的CH4浓度一般不超过5×104ppm。可见二者的浓度相差了3个甚至4个数量级。
电化学传感器是检测煤矿中CO和CH4的常用方法,要针对两种气体选择独立的电化学传感器,但是在使用过程中,电化学传感器不具备良好的选择性,并且其易中毒老化,使用寿命短,一般需要半年或者一年更换一次,其测量精度也不高。另一种常用方法是气相色谱技术,其操作复杂,设备体积庞大,分析周期长,难以实现实时原位在线检测。所以传统方法对于煤矿中微量CO和高浓度CH4无法实现实时在线、高精度同时检测。
可调谐二极管激光吸收光谱(Tunable Diode Laser Absorption Spectroscopy,TDLAS)技术采用窄线宽可调谐的激光器作为光源,通过控制温度和电流来获得合适的激光输出波长,当激光穿过对该波长有吸收的气体介质时,能量被吸收。根据Beer-Lambert吸收定律,可通过测量透射光谱信号,得到光束路径上的气体吸收强度,进而反演得到气体浓度参数。实际应用中,为了能够测量得到浓度更低的微量气体,往往使用波长调制技术。波长调制技术具有抗干扰能力强、灵敏度高等优点,但是波长调制技术一般要求其测量的气体吸光度要小于0.05,否则会产生明显的非线性。而在煤矿中微量CO的吸光度不会超过0.05,但是高浓度CH4的吸光度往往会远远超过0.05。并且在测量过程中,微量CO的二次谐波信号将被高浓度CH4混合气体的二次谐波信号严重干扰,导致无法精确测量CO浓度。中科院等相关研究机构已将波长调制技术应用到了CO、CO2等气体的浓度检测上,也有对煤矿等工况条件下的多种气体测量相关文献,但是都未讨论两种或者多种气体浓度差别巨大条件下的光谱解调和浓度反演。还未见有对煤矿等特殊工况中的微量CO和高浓度CH4气体(浓度差别3-4个数量级)混叠光谱信号进行解调,并同时测量两种气体浓度的文献报道。
在特殊工况中,多种气体同时存在,并且浓度差别在3-4个数量级。例如:在煤矿中,微量CO和高浓度CH4是同时存在的,并且高浓度CH4会严重干扰微量CO的信号,这就导致使用同一套系统测量对两者进行测量成为了一个难题。
因此亟需提供一种新型的多种气体同时检测的方法来解决上述问题。
发明内容
本发明所要解决的技术问题是提供一种激光多气体混叠光谱解调和非线性补偿的浓度检测方法,能够实现多种气体二次谐波混叠光谱的解调,并基于解调得到的结果对浓度进行非线性补偿,得到精确的气体浓度信息。
为解决上述技术问题,本发明采用的一个技术方案是:提供一种激光多气体混叠光谱解调和非线性补偿的浓度检测方法,所述待测气体为多气体混叠光谱信号,包括气体A、气体B,且气体A为单一谱线吸收气体、吸光度小于0.05,气体B的吸光度大于0.05,所述检测方法包括以下步骤:
S1:采集气体A、气体B标准气体的二次谐波作为参考信号,此时的气体A和气体B的浓度分别记为CR气体A、CR气体B,二次谐波信号分别记为向量SR气体A、SR气体B;
S2:采集待测气体的实时二次谐波信号,记为S;
S3:根据步骤S2中的实时二次谐波信号S和步骤S1中的气体A、气体B参考信号SR气体A、SR气体B,基于非负最小二乘法求出气体A、气体B的浓度,分别为CM气体A、CM气体B;
S4:采集四组不同标准浓度C1、C2、C3、C4的气体B的二次谐波信号,利用最小二乘法计算这四组气体B的浓度,记为M1、M2、M3、M4;
S5:对步骤S4中的四组标准浓度C1、C2、C3、C4和四组计算浓度M1、M2、M3、M4进行三次多项式拟合,得到气体B浓度的非线性修正模型;
S6:根据步骤S5中的非线性修正模型对步骤S3中求出的气体B浓度CM气体B进行修正,得到最终的气体B浓度C气体B;
S7:根据修正前后的气体B浓度之间的比例关系,对混合气体中的气体A浓度进行修正,得到最终的气体A浓度C气体A。
在本发明一个较佳实施例中,在步骤S1、S2和S4中,采集气体的二次谐波信号均进行数据预处理,包括背景光强消除、光谱校正、多次平均、S-G滤波。
在本发明一个较佳实施例中,在步骤S3中,计算实时二次谐波信号中气体A、气体B的浓度CM气体A、CM气体B的具体步骤包括:
首先实时二次谐波信号S记为步骤S1中标准气体A和气体B气体的二次谐波的比例叠加,表示为以下形式:
S=k1×SR气体A+k2×SR气体B=[k1,k2]T[SR气体A,SR气体B]
其中,k1、k2为比例系数;
然后基于非负最小二乘法求解比例系数k1、k2,公式如下:
minβ||[k1,k2]T[SR气体A,SR气体B]-S||2 2s.t.k1≥0,k2≥0;
最后根据下列两式求出实时二次谐波信号中气体A、气体B的浓度CM气体A、CM气体B:
CM气体A=k1×CR气体A
CM气体B=k2×CR气体B。
在本发明一个较佳实施例中,步骤S4的具体步骤包括:
将四组不同标准浓度分别为C1、C2、C3、C4的气体B气体的二次谐波信号与步骤S1中的气体B参考信号进行线性最小二乘拟合,得到系数为gi(其中i=1,2,3,4),进而根据下式求出对应的结果分别为M1、M2、M3、M4,
Mi=gi×CR气体B,i=1,2,3,4。
在本发明一个较佳实施例中,在步骤S5中,得到气体B浓度的非线性修正模型的具体步骤包括:
根据泰勒展开公式,利用三次多项式进行拟合,最终的非线性修正模型表示为:
M=f0+f1×C+f2×C2+f3×C3
其中,M表示计算得到的浓度结果,C表示标准气体的浓度,f0,f1,f2,f3分别为三次多项式的0次、1次、2次、3次系数。
在本发明一个较佳实施例中,在步骤S7中,对混合气体中的气体A浓度进行修正的公式如下:
本发明完成混叠光谱解调的原理是基于非负最小二乘算法,完成非线性补偿是基于比尔朗伯定律和波长调制技术的原理及其泰勒展开公式。
本发明的有益效果是:
(1)本发明能够实现多种气体二次谐波混叠光谱的解调,并基于解调得到的结果对浓度进行非线性补偿,得到精确的气体浓度信息;该方法能够弥补现有方法的不足,能够满足特殊工况环境中多种气体同时在线监测的需求;
(2)本发明利用一套波长调制系统完成了煤矿中CO和CH4的同时在线监测,该方法的创新型在于:1.利用非负最小二乘法完成了两种气体二次谐波混叠光谱的解调;2.对波长调制技术进行探究,发现其在测量气体浓度时,如果气体吸光度超过0.05,根据泰勒展开公式,就必须引入后面二次项、三次项甚至更高次项来满足精确测量的需求,基于此,本发明采用三次多项式对浓度进行非线性补偿。
附图说明
图1是本发明一种激光多气体混叠光谱解调和非线性补偿的浓度检测方法的流程图;
图2是单一CO气体、单一CH4气体、两者混合气体的二次谐波信号示意图;
图3是CH4浓度的非线性修正模型示意图。
具体实施方式
下面结合附图对本发明的较佳实施例进行详细阐述,以使本发明的优点和特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚明确的界定。
请参阅图1和图2,本发明实施例包括:
一种激光多气体混叠光谱解调和非线性补偿的浓度检测方法,包括以下步骤:
本示例以煤矿等特殊工况中的微量CO和高浓度CH4气体(浓度差别为3—4个数量级)的混叠光谱信号为例,同样本发明所述浓度检测方法还适用于其它多气体混叠光谱信号,其中具有一种气体吸光度小于0.05的微量(体积分数)气体(A)、一种气体吸光度远超0.05的高浓度气体(B),微量气体收到高浓度气体的干扰,具体如下表所示:
S1:采集CO、CH4标准气体的二次谐波作为参考信号,此时的CO和CH4浓度分别记为CRCO、CRCH4,二次谐波信号分别记为向量SRCO、SRCH4;
其中,采集已知浓度的CO和CH4标准气体的二次谐波信号进行了数据预处理,包括背景光强消除、光谱校正、多次平均、S-G滤波。以下步骤S2、S4中采集的二次谐波信号均进行同样的数据预处理。
S2:采集待测气体的实时二次谐波信号,记为S;
在一定条件下,可以理解为此时得到的二次谐波信号S是步骤S1中标准CO和CH4气体的二次谐波的比例叠加,可以表示为以下形式:
S=k1×SRCO+k2×SRCH4=[k1,k2]T[SRCO,SRCH4]
其中,k1、k2为比例系数。
S3:根据步骤S2中的实时二次谐波信号S和步骤S1中的CO、CH4参考信号SRCO、SRCH4,基于非负最小二乘法求出CO、CH4的浓度,分别为CMCO、CMCH4;
非负最小二乘方法是一种求解步骤S2中k1、k2两个比例系数的带有约束的线性回归方法,其可以表示为:
利用该方法求解得到两个比例系数的值之后,根据下列两式求出CO、CH4的浓度。
CMCO=k1×CRCO
CMCH4=k2×CRCH4
S4:采集四组不同标准浓度C1、C2、C3、C4的CH4气体的二次谐波信号,利用最小二乘法计算这四组CH4气体的浓度,记为M1、M2、M3、M4;具体步骤包括:
采集四组不同标准浓度CH4气体(浓度分别为C1,C2,C3,C4)的二次谐波信号,并进行和步骤S1中相同的数据预处理;将这四组信号分别与步骤1中的CH4参考信号,利用最小二乘法求出斜率系数分别为g1,g2,g3,g4,进而根据下式求出对应的结果分别为M1,M2,M3,M4。
Mi=gi×CRCH4 i=1,2,3,4
S5:对步骤S4中的四组标准浓度C1、C2、C3、C4和四组计算浓度M1、M2、M3、M4进行三次多项式拟合,得到CH4浓度的非线性修正模型,参阅图3;具体步骤包括:
M=f0+f1×C+f2×C2+f3×C3
其中,M表示计算得到的浓度结果,C表示标准气体的浓度,f0,f1,f2,f3分别为三次多项式的0次、1次、2次、3次系数。
S6:利用步骤S5中得到的非线性修正模型,将步骤S3中计算得到的CH4浓度CMCH4带入到以上公式中,即可计算出对应更加接近标准结果的更精确的CH4浓度,即为待测的CH4浓度,记录为CCH4
S7:根据修正前的CH4浓度CMCH4和修正后的CH4浓度CCH4之间的比例关系,对混合气体中的CO浓度进行修正,公式如下,得到最终的CO浓度CCO。
由于此时计算得到的CO浓度误差主要来源于高浓度的CH4所产生的非线性,因此利用修正前后的CH4浓度之间的比例关系对CO浓度进行修正,能够得到更加精确的结果。
本发明利用非负最小二乘方法对微量CO和CH4混合气体的二次谐波混叠光谱进行处理,得到两者各自的浓度,之后基于二次谐波的非线性,采用三次多项式对得到的浓度进行补偿,最终得精确的气体浓度数据。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (6)
1.一种激光多气体混叠光谱解调和非线性补偿的浓度检测方法,其特征在于,所述待测气体为多气体混叠光谱信号,包括气体A、气体B,且气体A为单一谱线吸收气体、吸光度小于0.05,气体B的吸光度大于0.05,所述检测方法包括以下步骤:
S1:采集气体A、气体B标准气体的二次谐波作为参考信号,此时的气体A和气体B的浓度分别记为CR气体A、CR气体B,二次谐波信号分别记为向量SR气体A、SR气体B;
S2:采集待测气体的实时二次谐波信号,记为S;
S3:根据步骤S2中的实时二次谐波信号S和步骤S1中的气体A、气体B参考信号SR气体A、SR气体B,基于非负最小二乘法求出气体A、气体B的浓度,分别为CM气体A、CM气体B;
S4:采集四组不同标准浓度C1、C2、C3、C4的气体B的二次谐波信号,利用最小二乘法计算这四组气体B的浓度,记为M1、M2、M3、M4;
S5:对步骤S4中的四组标准浓度C1、C2、C3、C4和四组计算浓度M1、M2、M3、M4进行三次多项式拟合,得到气体B浓度的非线性修正模型;
S6:根据步骤S5中的非线性修正模型对步骤S3中求出的气体B浓度CM气体B进行修正,得到最终的气体B浓度C气体B;
S7:根据修正前后的气体B浓度之间的比例关系,对混合气体中的气体A浓度进行修正,得到最终的气体A浓度C气体A。
2.根据权利要求1所述的激光多气体混叠光谱解调和非线性补偿的浓度检测方法,其特征在于,在步骤S1、S2和S4中,采集气体的二次谐波信号均进行数据预处理,包括背景光强消除、光谱校正、多次平均、S-G滤波。
3.根据权利要求1所述的激光多气体混叠光谱解调和非线性补偿的浓度检测方法,其特征在于,在步骤S3中,计算实时二次谐波信号中气体A、气体B的浓度CM气体A、CM气体B的具体步骤包括:
首先实时二次谐波信号S记为步骤S1中标准气体A和气体B气体的二次谐波的比例叠加,表示为以下形式:
S=k1×SR气体A+k2×SR气体B=[k1,k2]T[SR气体A,SR气体B]
其中,k1、k2为比例系数;
然后基于非负最小二乘法求解比例系数k1、k2,公式如下:
minβ||[k1,k2]T[SR气体A,SR气体B]-S||2 2s.t.k1≥0,k2≥0;
最后根据下列两式求出实时二次谐波信号中气体A、气体B的浓度CM气体A、CM气体B:
CM气体A=k1×CR气体A
CM气体B=k2×CR气体B。
4.根据权利要求1所述的激光多气体混叠光谱解调和非线性补偿的浓度检测方法,其特征在于,步骤S4的具体步骤包括:
将四组不同标准浓度分别为C1、C2、C3、C4的气体B气体的二次谐波信号与步骤S1中的气体B参考信号进行线性最小二乘拟合,得到系数为gi(其中i=1,2,3,4),进而根据下式求出对应的结果分别为M1、M2、M3、M4,
Mi=gi×CR气体B,i=1,2,3,4。
5.根据权利要求1所述的激光多气体混叠光谱解调和非线性补偿的浓度检测方法,其特征在于,在步骤S5中,得到气体B浓度的非线性修正模型的具体步骤包括:
根据泰勒展开公式,利用三次多项式进行拟合,最终的非线性修正模型表示为:
M=f0+f1×C+f2×C2+f3×C3
其中,M表示计算得到的浓度结果,C表示标准气体的浓度,f0,f1,f2,f3分别为三次多项式的0次、1次、2次、3次系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210380189.6A CN114813633A (zh) | 2022-04-12 | 2022-04-12 | 一种激光多气体混叠光谱解调和非线性补偿的浓度检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210380189.6A CN114813633A (zh) | 2022-04-12 | 2022-04-12 | 一种激光多气体混叠光谱解调和非线性补偿的浓度检测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114813633A true CN114813633A (zh) | 2022-07-29 |
Family
ID=82534190
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210380189.6A Pending CN114813633A (zh) | 2022-04-12 | 2022-04-12 | 一种激光多气体混叠光谱解调和非线性补偿的浓度检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114813633A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115753678A (zh) * | 2023-01-10 | 2023-03-07 | 杭州泽天春来科技有限公司 | 气体分析方法 |
CN116735528A (zh) * | 2023-08-16 | 2023-09-12 | 中国科学院合肥物质科学研究院 | 一种烷烃类气体混叠光谱解析方法 |
CN117686459A (zh) * | 2023-12-08 | 2024-03-12 | 数量级(无锡)信息技术有限公司 | 一种多气体同步测量方法及装置 |
-
2022
- 2022-04-12 CN CN202210380189.6A patent/CN114813633A/zh active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115753678A (zh) * | 2023-01-10 | 2023-03-07 | 杭州泽天春来科技有限公司 | 气体分析方法 |
CN116735528A (zh) * | 2023-08-16 | 2023-09-12 | 中国科学院合肥物质科学研究院 | 一种烷烃类气体混叠光谱解析方法 |
CN116735528B (zh) * | 2023-08-16 | 2024-01-05 | 中国科学院合肥物质科学研究院 | 一种烷烃类气体混叠光谱解析方法 |
CN117686459A (zh) * | 2023-12-08 | 2024-03-12 | 数量级(无锡)信息技术有限公司 | 一种多气体同步测量方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114813633A (zh) | 一种激光多气体混叠光谱解调和非线性补偿的浓度检测方法 | |
US4271124A (en) | Non-dispersive infrared gas analyzer for testing gases containing water-vapor | |
CN103323115A (zh) | 基于波长调制的气体吸收谱线线宽和线型系数的测量方法 | |
US3924442A (en) | Pollutant variation correcting system | |
CN109557045B (zh) | 一种用于提取tdlas技术中直接吸收谱线的拟合方法 | |
CN104990886B (zh) | 一种傅里叶变换红外光谱仪波数漂移的校正方法 | |
CN102680020A (zh) | 一种基于波长调制光谱技术的气体参数在线测量方法 | |
CN104730021A (zh) | 一种非分光红外气体传感器的标定方法 | |
CN102183468A (zh) | 多组分气体分析的干扰修正与浓度反演方法 | |
CN105203460B (zh) | 红外激光光谱痕量水汽检测系统及其检测方法 | |
CN107091818A (zh) | 一种多气室复杂组分气体分析系统及方法 | |
CN103543124A (zh) | 一种基于软件锁相的可调激光吸收光谱气体检测方法 | |
CN107478593A (zh) | 一种低浓度no和so2混合气体的浓度检测方法 | |
CN105067564A (zh) | 一种具有温度补偿能力的光纤气体浓度检测方法 | |
CN116735528B (zh) | 一种烷烃类气体混叠光谱解析方法 | |
CN114942235B (zh) | 一种用于复杂背景环境提取有效吸收信息的方法 | |
CN103558182B (zh) | 一种用于激光气体在线分析仪确定气体浓度的方法 | |
CN109632680A (zh) | 一种基于排列熵的水体中磷的检测方法 | |
Wang et al. | Multi-component gas measurement aliasing spectral demodulation method for interference separation in laser absorption spectroscopy | |
Zhao et al. | Method for demodulating the overlapping absorption spectra of CO and CH 4 | |
CN212364040U (zh) | 一种用于高精度测量气体浓度的滤光片 | |
CN114323507B (zh) | 一种西林瓶密封完整性测量装置和方法 | |
CN111220545A (zh) | 一种用于高精度测量气体浓度的滤光片 | |
KR20090051244A (ko) | 산화 에틸렌을 제조하기 위한 공정의 개선된 제어 및 최적화 방법 | |
CN114970168A (zh) | 一种变压器油在线监测装置评估方法 |
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 |