CN116500665A - 傅里叶尺度变换特性的微剂量探测器组织等效换算方法 - Google Patents
傅里叶尺度变换特性的微剂量探测器组织等效换算方法 Download PDFInfo
- Publication number
- CN116500665A CN116500665A CN202310561795.2A CN202310561795A CN116500665A CN 116500665 A CN116500665 A CN 116500665A CN 202310561795 A CN202310561795 A CN 202310561795A CN 116500665 A CN116500665 A CN 116500665A
- Authority
- CN
- China
- Prior art keywords
- deposition energy
- micro
- energy
- dose detector
- equivalent conversion
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/36—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T7/00—Details of radiation-measuring instruments
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明公开了基于傅里叶尺度变换特性的微剂量探测器组织等效换算方法,属于辐射剂量测量技术领域,包括:获取给定条件下入射辐射在微剂量探测器中的沉积能量分布D(x),并计算入射辐射在微剂量探测器和生物组织中的平均沉积能量比;将D(x)进行尺度变换,并将D(x)、进行尺度变换后得到的H(x)进行傅里叶变换,得到D(x)、H(x)分别在傅里叶空间中的模MD、MH;根据平均沉积能量比对模MD、MH差值进行修正,然后将修正后的差值进行傅里叶逆变换、归一化处理;最后将归一化处理的结果与H(x)进行卷积,得到生物组织中的沉积能谱T(x)。本发明能够准确可靠的将微剂量探测器测得的沉积能量谱转换为生物组织中沉积能量谱,有效改善当前组织等效换算方法在低能区换算时存在偏差的问题,实现测量目标。
Description
技术领域
本发明涉及辐射剂量测量与辐射防护技术领域,尤其涉及一种基于傅里叶尺度变换特性的微剂量探测器组织等效换算方法。
背景技术
固体微剂量探测器由于能直接将探测器灵敏体积制成生物组织大小,可以准确表征入射辐射在微体积内的相互作用,在近年来引发了广泛的关注。不过,由于探测器材料(主要为硅和金刚石)与生物组织材料之间存在差值,如物理特性、密度、电离能等,致使放射性粒子或射线在微剂量探测器中的沉积能量分布与生物组织中的沉积能量分布并不一致。因此,为获得生物组织微观体积内的沉积能量分布,需要对微剂量探测器的测量结果进行组织等效换算。
当前,关于微剂量探测器测得的沉积能谱向同等大小生物组织的换算,有学者提出了一种基于谱图分布变换的换算方法。由于在换算过程中未对低能沉积事件的特点加以考虑,致使低能区的换算结果并不理想,即无法有效将微剂量探测器测得的沉积能量谱转换为生物组织中沉积能量谱,特别在低能区换算时存在较大误差。
发明内容
本发明的目的在于克服现有技术的问题,提供了一种基于傅里叶尺度变换特性的微剂量探测器组织等效换算方法。
本发明的目的是通过以下技术方案来实现的:傅里叶尺度变换特性的微剂量探测器组织等效换算方法,该方法包括以下步骤:
获取给定条件(一定能量射线粒子辐照场)下入射辐射在微剂量探测器中的沉积能量分布D(x);
计算当前条件下入射辐射在微剂量探测器和生物组织中的平均沉积能量比λ;
根据平均沉积能量比λ对测得的沉积能量分布D(x)进行尺度变换,得到尺度变换后的沉积能量分布H(x);
将D(x)、H(x)进行傅里叶变换,得到D(x)、H(x)分别在傅里叶空间中的模MD、MH;
计算模MD和MH之间的差值,并根据平均沉积能量比对差值进行修正得到ME;
将修正后的差值ME进行傅里叶逆变换得到沉积能量分布差E(x),再对E(x)进行归一化处理得到C(x);
将归一化处理的结果C(x)与尺度变换后的沉积能量分布H(x)进行卷积,得出同等条件下与微剂量探测器大小相同的生物组织中的沉积能谱T(x)。
在一示例中,所述平均沉积能量比λ的计算式为:
式中,ε表示入射辐射在介质中的沉积能量;ε(Ein,微剂量探测器)表示入射能量为Ein时,微剂量探测器中的沉积能量;ε(Ein,组织)表示入射能量为Ein时,生物组织中的沉积能量。
在一示例中,粒子平均沉积能量ε的计算式为:
ΔE=Ei-Ei+1
ε=E0-En
式中,L表示入射辐射在介质(微剂量探测器或生物组织)中穿过的距离;ΔE表示能量间隔;E0表示入射辐射进入介质之前的能量;En表示入射辐射离开介质时的能量;i表示第i个能量间隔点;S(Ei)表示能量大小为Ei时,介质对入射粒子的线性阻止本领。
在一示例中,对测得的沉积能量分布D(x)进行尺度变换的计算式为:
H(x)=λD(λx)
式中,x表征粒子沉积能量。
在一示例中,计算模MD和MH之间的差值,并根据平均沉积能量比对差值进行修正的计算式为:
在一示例中,沉积能谱T(x)的计算式为:
T(x)=C(x)*H(x)
式中,x表征粒子沉积能量。
在一示例中,傅里叶变换、傅里叶空间中模的计算公式为:
式中,1≤k≤N,表示粒子沉积能量数据点x在傅里叶空间对应的数据点,N表示数据长度,x表征粒子沉积能量。
在一示例中,傅里叶逆变换计算公式为:
式中,1≤k≤N,表示粒子沉积能量数据点x在傅里叶空间对应的数据点,N表示数据维度,x表征粒子沉积能量。
在一示例中,剔除沉积能量分布差E(x)中的负值。
在一示例中,剔除E(x)负值后进行归一化处理的计算式为:
式中,E1(x)表示剔除E(x)分布中负值后得到的分布;C(x)表示归一化后的分布,x表征粒子沉积能量。
需要进一步说明的是,上述方法各示例对应的技术特征可以相互组合或替换构成新的技术方案。
与现有技术相比,本发明有益效果是:
本发明方法能够广泛应用于固体微剂量探测器的组织等效修正中,使探测器适用于粒子全射程内测量;结合仿真结果可知,本发明方法能够准确可靠的将微剂量探测器测得的沉积能量谱转换为生物组织中沉积能量谱,有效改善当前组织等效换算方法在低能区换算时存在偏差的问题,实现测量目标;整个过程操作简便,在一定程度上降低测量成本。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明,此处所说明的附图用来提供对本申请的进一步理解,构成本申请的一部分,在这些附图中使用相同的参考标号来表示相同或相似的部分,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。
图1为50MeV质子在骨骼体膜中的百分深度剂量曲线图;
图2为本发明公开的一种基于傅里叶尺度变换特性的微剂量探测器组织等效换算方法的流程图;
图3为百分深度剂量曲线坪区金刚石微剂量探测器的等效换算结果与骨骼中沉积能谱的对比图;
图4为百分深度剂量曲线布拉格峰前端处金刚石微剂量探测器的等效换算结果与骨骼中沉积能谱的对比图;
图5为百分深度剂量曲线布拉格峰后端处金刚石微剂量探测器的等效换算结果与骨骼中沉积能谱的对比图。
具体实施方式
下面结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要说明的是,属于“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方向或位置关系为基于附图所述的方向或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,使用序数词(例如,“第一和第二”、“第一至第四”等)是为了对物体进行区分,并不限于该顺序,而不能理解为指示或暗示相对重要性。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,属于“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
此外,下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。
本实施例以10×10×10μm3大小金刚石微剂量探测器在50MeV质子辐照场中测得的能谱向同等大小生物组织中沉积能量等效换算为例,采用本发明方法对金刚石微剂量探测器在骨骼百分深度剂量曲线坪区、布拉格峰前端、布拉格峰后端测得的能量沉积分布进行组织等效换算。其中,50MeV质子照射下骨骼的百分深度剂量曲线如图1所示,它表示射线粒子(50MeV质子)在骨骼中某一深度处的吸收剂量与参考深度处吸收剂量之比的百分数,是描述射线中心轴不同深度处相对剂量分布的物理量。
此时,如图2所示,本发明一种基于傅里叶尺度变换特性的微剂量探测器组织等效换算方法包括以下步骤:
S1:获取金刚石微剂量探测器在百分深度剂量曲线坪区、布拉格峰前端、布拉格峰后端处的沉积能量分布D1(x)、D2(x)、D3(x);
本实施例中,选取的布拉格曲线坪区、布拉格峰前端、布拉格峰后端测量点的三个测量点在体膜中的具体深度参数如表1所示:
表1选取的测量点在骨骼体膜中的深度位置
S2:计算在测量点1、测量点2、测量点3处入射能量下金刚石微剂量探测器与骨骼之间的平均沉积能量比,λ1、λ2、λ3;
S3:根据计算得到的平均沉积能量比将测量点1、测量点2、测量点3处测得的沉积能量分布D1(x)、D2(x)、D3(x)进行尺度变换,得到H1(x)、H2(x)、H3(x);
S4:将分布D1(x)、D2(x)、D3(x)和H1(x)、H2(x)、H3(x)进行傅里叶变换并求得变换后结果在傅里叶空间里的模MD1、MD2、MD3和MH1、MH2、MH3;
S5:计算模MD1和MH1、MD2和MH2、MD3和MH3之间的差异并根据平均沉积能量比对差异进行修正得到ME1、ME2、ME3;
S6:将ME1、ME2、ME3进行傅里叶逆变换得到E1(x)、E2(x)、E3(x),并对E1(x)、E2(x)、E3(x)进行归一化处理得到C1(x)、C2(x)、C3(x);
S7:将归一化处理的结果C1(x)、C2(x)、C3(x)分别与尺度变换后的沉积能量分布H1(x)、H2(x)、H3(x)进行卷积计算,得出测量点1、测量点2、测量点3处骨骼中的沉积能量分布T1(x)、T2(x)、T3(x)。
为验证本发明方法的组织等效换算效果,通过软件模拟获取各测量点处质子在相同大小骨骼中的沉积能量分布B1(x)、B2(x)、B3(x),将运用本发明提出的方法得到的沉积能量分布T1(x)、T2(x)、T3(x)与B1(x)、B2(x)、B3(x)进行对比,结果如图3-5所示,图3-5中横坐标ε表示骨骼中沉积的能量,纵坐标F(ε)表示沉积能量的概率分布。其中图3为测量点1处金刚石微剂量探测器的等效换算结果与骨骼中沉积能量分布的对比;图4为测量点2处金刚石微剂量探测器的等效换算与骨骼中沉积能量分布的对比;图5为测量点3处金刚石微剂量探测器的等效换算结果与骨骼中沉积能量分布的对比。综合上述图3-5可以看出,通过本发明方法对微剂量探测器进行组织等效换算后得到的沉积能量分布与相同条件下骨骼中的沉积能量分布吻合度较好,尤其是在测量点1、2处低能区,换算结果与实际结果之间的重合度较高,即:本发明方法能将探测器测得的沉积能量分布准确转换为相同条件下生物组织体中的沉积能量分布,且与现有换算方法相比,该方法可以有效消除低能区的换算偏差,提高换算结果与实际结果之间的吻合度。
以上具体实施方式是对本发明的详细说明,不能认定本发明的具体实施方式只局限于这些说明,对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演和替代,都应当视为属于本发明的保护范围。
Claims (10)
1.一种傅里叶尺度变换特性的微剂量探测器组织等效换算方法,其特征在于:其包括以下步骤:
获取给定条件下入射辐射在微剂量探测器中的沉积能量分布D(x);
计算当前条件下入射辐射在微剂量探测器和生物组织中的平均沉积能量比λ;
根据平均沉积能量比λ对测得的沉积能量分布D(x)进行尺度变换,得到尺度变换后的沉积能量分布H(x);
将D(x)、H(x)进行傅里叶变换,得到D(x)、H(x)分别在傅里叶空间中的模MD、MH;
计算模MD和MH之间的差值,并根据平均沉积能量比对差值进行修正得到ME;
将修正后的差值ME进行傅里叶逆变换得到沉积能量分布差E(x),再对E(x)进行归一化处理得到C(x);
将归一化处理的结果C(x)与尺度变换后的沉积能量分布H(x)进行卷积,得出同等条件下与微剂量探测器大小相同的生物组织中的沉积能谱T(x)。
2.根据权利要求1所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法,其特征在于:所述平均沉积能量比λ的计算式为:
式中,ε表示入射辐射在介质中的沉积能量;ε(Ein,微剂量探测器)表示入射能量为Ein时,微剂量探测器中的沉积能量;ε(Ein,组织)表示入射能量为Ein时,生物组织中的沉积能量。
3.根据权利要求1所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法,其特征在于:粒子平均沉积能量ε的计算式为:
ΔE=Ei-Ei+1
ε=E0-En
式中,L表示入射辐射在介质中穿过的距离;ΔE表示能量间隔;E0表示入射辐射进入介质之前的能量;En表示入射辐射离开介质时的能量;i表示第i个能量间隔点;S(Ei)表示能量大小为Ei时,介质对入射粒子的线性阻止本领。
4.根据权利要求1所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法,其特征在于:对测得的沉积能量分布D(x)进行尺度变换的计算式为:
H(x)=λD(λx)
式中,x表征粒子沉积能量。
5.根据权利要求1所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法,其特征在于:计算模MD和MH之间的差值,并根据平均沉积能量比对差值进行修正的计算式为:
6.根据权利要求1所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法,其特征在于:沉积能谱T(x)的计算式为:
T(x)=C(x)*H(x)
式中,x表征粒子沉积能量。
7.根据权利要求1所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法,其特征在于:傅里叶变换、傅里叶空间中模的计算公式为:
式中,1≤k≤N,表示粒子沉积能量数据点x在傅里叶空间对应的数据点,N表示数据长度,x表征粒子沉积能量。
8.根据权利要求1所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法,其特征在于:傅里叶逆变换计算公式为:
式中,1≤k≤N,表示粒子沉积能量数据点x在傅里叶空间对应的数据点,N表示数据维度,x表征粒子沉积能量。
9.根据权利要求1所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法,其特征在于:剔除沉积能量分布差E(x)中的负值。
10.根据权利要求9所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法,其特征在于:剔除E(x)负值后进行归一化处理的计算式为:
式中,E1(x)表示剔除E(x)分布中负值后得到的分布;C(x)表示归一化后的分布,x表征粒子沉积能量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310561795.2A CN116500665A (zh) | 2023-05-17 | 2023-05-17 | 傅里叶尺度变换特性的微剂量探测器组织等效换算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310561795.2A CN116500665A (zh) | 2023-05-17 | 2023-05-17 | 傅里叶尺度变换特性的微剂量探测器组织等效换算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116500665A true CN116500665A (zh) | 2023-07-28 |
Family
ID=87318230
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310561795.2A Pending CN116500665A (zh) | 2023-05-17 | 2023-05-17 | 傅里叶尺度变换特性的微剂量探测器组织等效换算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116500665A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005119295A1 (en) * | 2004-06-04 | 2005-12-15 | Bc Cancer Agency | Method and apparatus for verifying radiation dose distributions |
US20100090118A1 (en) * | 2006-12-19 | 2010-04-15 | University Of Wollongong | Method and apparatus for tissue equivalent solid state microdosimetry |
CN114034721A (zh) * | 2021-11-10 | 2022-02-11 | 成都大学 | 一种微剂量探测器的组织等效修正方法 |
CN114518589A (zh) * | 2022-01-28 | 2022-05-20 | 西北核技术研究所 | 基于厚放射源实现气体正比探测器能量刻度的方法 |
-
2023
- 2023-05-17 CN CN202310561795.2A patent/CN116500665A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005119295A1 (en) * | 2004-06-04 | 2005-12-15 | Bc Cancer Agency | Method and apparatus for verifying radiation dose distributions |
US20100090118A1 (en) * | 2006-12-19 | 2010-04-15 | University Of Wollongong | Method and apparatus for tissue equivalent solid state microdosimetry |
CN114034721A (zh) * | 2021-11-10 | 2022-02-11 | 成都大学 | 一种微剂量探测器的组织等效修正方法 |
CN114518589A (zh) * | 2022-01-28 | 2022-05-20 | 西北核技术研究所 | 基于厚放射源实现气体正比探测器能量刻度的方法 |
Non-Patent Citations (2)
Title |
---|
余松科 等: "固体微剂量探测器组织等效换算教学实验设计", 实验技术与管理, vol. 40, no. 2, 28 February 2023 (2023-02-28), pages 170 - 174 * |
雷鸣 等: "SOI硅微剂量探测器对中子和伽马辐射场线能谱测量的GEANT4模拟研究", 辐射防护, vol. 37, no. 3, 31 May 2017 (2017-05-31), pages 169 - 173 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Whalin et al. | Photodisintegration of Deuterium by 60-to 250-MeV X-Rays | |
WO2021069922A1 (en) | Methods and apparatus for particle beam dose profile measurement | |
Kaune et al. | Inclusive cross sections for pion and proton production by photons using collimated coherent bremsstrahlung | |
Casnati et al. | Measurements of the total atomic differential cross section of elastic scattering of 59.54-keV photons | |
CN116500665A (zh) | 傅里叶尺度变换特性的微剂量探测器组织等效换算方法 | |
CN114034721B (zh) | 一种微剂量探测器的组织等效修正方法 | |
Amann et al. | Inelastic π+ scattering from C 12 and Si at low energies | |
CN114185078B (zh) | 一种1MeV以上单能中子注量率的测量方法 | |
Yu et al. | A method for converting microdosimetric spectra in diamond to tissue in proton therapy | |
Bianchi et al. | Microdosimetric measurements for LET monitoring in proton therapy. The development of engineered mini-TEPCs for clinical applications: First results | |
Burcham et al. | LXXVI. The emission of short-range alpha particles from light elements under proton bombardment.—I. Experimental method and the reaction 10B (pα) 7Be | |
Livesey et al. | An investigation of the neutrons produced in the deuteron-deuteron reaction | |
CN110221338B (zh) | 一种重构屏蔽体周边区域辐射场的方法 | |
Lam et al. | Study of the E p= 1.724 MeV Resonance in the 27Al (p, γ) 28Si Reaction | |
CN118447968A (zh) | 非组织等效材料的组织等效修正方法 | |
CN110412639A (zh) | 基于纳剂量学获得离子束辐照方案的方法 | |
Ahmad et al. | The energy dependence of fission fragment anisotropy in fast-neutron-induced fission of uranium-235 | |
KR20070093400A (ko) | 고속 전자들의 로컬 방사율 프로파일을 결정하는 방법 | |
Ma et al. | Differential and integral cross sections for the valence-shell excitations in D 2 O studied by fast electron impact | |
Shaw et al. | The Clatterbridge high-energy neutron therapy facility: measurements of beam parameters for clinical use | |
US8898038B2 (en) | Determination of beam parameters for unflattened photon beams | |
Friedman | Studies of Neutron Thermalization in H2O by the Pulsed Source and “Non 1/ν” Absorbers Method | |
EP3821276B1 (en) | Silicon carbide ionizing radiation detector | |
CN114371498A (zh) | 一种建立低剂量率水平辐射质的方法 | |
Sheppard et al. | A theory for heterogeneous ion chambers and neutron response of Victoreen dosimeters |
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 |