CN115437006A - 一种煤岩界面折射波的多阶折射时间求取方法 - Google Patents
一种煤岩界面折射波的多阶折射时间求取方法 Download PDFInfo
- Publication number
- CN115437006A CN115437006A CN202210907743.1A CN202210907743A CN115437006A CN 115437006 A CN115437006 A CN 115437006A CN 202210907743 A CN202210907743 A CN 202210907743A CN 115437006 A CN115437006 A CN 115437006A
- Authority
- CN
- China
- Prior art keywords
- refraction
- refracted
- order
- wave
- coal
- 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
- 239000003245 coal Substances 0.000 title claims abstract description 57
- 239000011435 rock Substances 0.000 title claims abstract description 24
- 238000004364 calculation method Methods 0.000 title claims description 5
- 238000000034 method Methods 0.000 claims abstract description 33
- 239000011159 matrix material Substances 0.000 claims abstract description 26
- 239000013598 vector Substances 0.000 claims description 23
- 230000005284 excitation Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 4
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 5
- 230000000737 periodic effect Effects 0.000 description 7
- 230000002159 abnormal effect Effects 0.000 description 6
- 238000005065 mining Methods 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- 238000001514 detection method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
- G01V1/305—Travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B17/00—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations
- G01B17/02—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations for measuring thickness
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种煤岩界面折射波的多阶折射时间求取方法。当震源波长大于煤厚时,相邻阶的折射波相互混叠,常规反褶积方法难以对其稀疏化。本发明利用各阶折射波相似性高、阶数已知的特点,借鉴正交匹配追踪算法,能够从严重混叠的折射波信号中准确提取各阶折射时间。本发明考虑到了各阶折射波之间的相位差异,通过构造包含相位旋转角的字典矩阵,使各阶折射波能准确匹配对应子波,获得准确的脉冲信号,提高所得时间点的准确性。
Description
技术领域
本发明涉及一种多阶折射时间求取方法,属于煤层地质探测技术领域,具体是涉及一种煤岩界面折射波的多阶折射时间求取方法。
背景技术
煤层厚度是采煤过程的重要参数。在回采前准确探测出工作面内煤层厚度的分布情况,是采煤机动态调整位置和采高,从而实现煤矿工作面智能开采的前提和基础,对提高回采效率、减少资源浪费以及保障煤矿安全生产均具有重要意义。
震源在煤层内激发后,产生的地震波将在煤层与顶板和底板的两个分界面之间产生多次反射,其中仅有以临界角入射的地震波能够产生沿界面传播的折射波;因此当煤层的厚度稳定时,接收到的地震信号的折射波部分将具有明显的周期性,即在围岩速度对应的时间点上出现初至折射波之后,每隔固定的时间周期出现续至折射波;当工作面内存在薄煤带等煤厚异常区域时,续至折射波出现的时间点会发生移动,偏离正常的周期,因此求出各阶折射波出现的时间,就可以由此计算出异常煤层的厚度。
专利“一种基于折射波周期振幅衰减的薄煤带探测方法及系统”(公开号:202111434047.5)利用在煤层与围岩界面上传播的折射波的周期特性,探测采煤工作面内薄煤带等煤厚异常区的位置和边界。此专利方法仅能圈定出薄煤带的范围,不能准确估计出薄煤带内具体的煤厚及其分布情况。这是由于此专利中对折射波信号进行反褶积的方法较为简单,反褶积后的折射波信号不是稀疏的,从中无法准确求出多阶折射发生的具体时间,因此不能确定各个折射点对应的煤层厚度。为此,需要针对折射波特性采用更有效的算法对折射波信号稀疏化,并从中求取多阶折射时间,为准确估计煤厚分布创造条件。
发明内容
以下给出一个或多个方面的简要概述以提供对这些方面的基本理解。此概述不是所有构想到的方面的详尽综览,并且既非旨在指认出所有方面的关键性或决定性要素亦非试图界定任何或所有方面的范围。其唯一的目的是要以简化形式给出一个或多个方面的一些概念以为稍后给出的更加详细的描述之序。
本发明主要的目的是解决现有技术中所存在的技术问题,提供了一种煤岩界面折射波的多阶折射时间求取方法,该方法简单易用,适应于现有槽波探测工程的观测系统与施工条件,能够将基于折射波的煤厚探测方法的适用性条件由原来的巨厚煤层推广到厚煤层与中厚煤层中。
为解决上述问题,本发明的方案是:
一种煤岩界面折射波的多阶折射时间求取方法,包括:
计算折射波的折射阶数K,并将子波wn(t)合成为字典矩阵W;
按以下迭代步骤计算得到稀疏信号r(t):
迭代步骤1,初始化,令k=1,A0为空集,u1=s;
迭代步骤3,基于下式求取迭代后的Ak:
Ak=Ak-1∪aik;
迭代步骤5,更新残差向量uk+1=s-Akrk;
迭代步骤6,令k=k+1,如果k≤K则返回至迭代步骤2,否则将当前向量rk作为稀疏信号r(t)输出;
迭代步骤7,拾取稀疏信号r(t)的K个脉冲峰值对应的时间,得到折射波的多阶折射时间;
其中,各迭代步骤中,k为当前迭代次数,rk表示第k次迭代的稀疏向量,uk表示第k次迭代的残差向量,ik表示第k次迭代找到的索引,所述索引为W矩阵中对应的列序号,aik表示矩阵W中以ik为序号的列向量,Ak为一个由所有挑选出的列向量aik所组成的集合矩阵。
优选的,上述的一种煤岩界面折射波的多阶折射时间求取方法,基于下式计算折射阶数K:
由已知的激发点位置(xs,ys)和接收点位置(xr,yr),计算折射波的折射阶数K:
式中,地震波速度vr、煤层的地震波速度vc,T为折射波周期。
优选的,上述的一种煤岩界面折射波的多阶折射时间求取方法,基于下式计算折射波周期:
其中,围岩的地震波速度vr、煤层的地震波速度vc。
优选的,上述的一种煤岩界面折射波的多阶折射时间求取方法,基于下式产生子波:
其中t∈[0,NΔt],周期相位旋转角Φ,折射波周期T,t表示时间,t∈[0,NΔt],N为采集的折射波信号的长度,Δt是采样间隔。
优选的,上述的一种煤岩界面折射波的多阶折射时间求取方法,基于下式产生字典矩阵W:
其中符号∪表示集合并运算。
其中,
主频为fc的雷克子波w0(τ),τ表示时间,τ∈[-T,T],T为折射波周期;
H(·)为希尔伯特变换,τ∈[-T,T]。
优选的,上述的一种煤岩界面折射波的多阶折射时间求取方法,基于下式构建主频为fc的雷克子波w0(τ):
其中,τ表示时间,τ∈[-T,T],T为折射波周期。
因此,相对于现有技术,本发明的优点是:
(1)本发明可实现震源波长大于煤厚条件下的折射波处理。当震源波长大于煤厚时,相邻阶的折射波相互混叠,常规反褶积方法难以对其稀疏化。本发明利用各阶折射波相似性高、阶数已知的特点,借鉴正交匹配追踪算法,能够从严重混叠的折射波信号中准确提取各阶折射时间。
(2)本发明考虑到了各阶折射波之间的相位差异,通过构造包含相位旋转角的字典矩阵,使各阶折射波能准确匹配对应子波,获得准确的脉冲信号,提高所得时间点的准确性。
附图说明
并入本文并形成说明书的一部分的附图例示了本发明的实施例,并且附图与说明书一起进一步用于解释本发明的原理以及使得所属领域技术人员能够制作和使用本公开。
图1为煤岩界面折射波的多阶折射时间求取方法的总流程图。
图2为工作面模型的示意图,其中,(a)为模型及观测系统平面图,(b)为模型剖面图;
图3为工作面模型不同主频震源的三维弹性波正演模拟结果,其中,(a)为200Hz主频震源产生的折射波信号,(b)为400Hz主频震源产生的折射波信号,(c)为R15接收到的两种主频的折射波信号的对比图;
图4为震源附近折射波信号及其频谱,其中,(a)为折射波信号,(b)为频谱;
图5为以200Hz为主频构建雷克子波;
图7为R15道折射波信号s(t)与用Φ对子波w0进行周期相位旋转得到的合成信号的对比图;
图8为字典矩阵W对应的图像;
图9为s(t)与得到的稀疏信号r(t)的示意图;
图10为对所有数据逐道处理,得到的稀疏信号;
将参照附图描述本发明的实施例。
具体实施方式
实施例
当震源波长大于煤厚时,相邻阶的折射波相互混叠,常规反褶积方法难以对其稀疏化,因而本实施例的方法利用各阶折射波相似性高、阶数已知的特点,借鉴正交匹配追踪算法,并且考虑到了各阶折射波之间的相位差异,通过构造包含相位旋转角的字典矩阵,使各阶折射波能准确匹配对应子波,获得准确的脉冲信号,提高所得时间点的准确性。
下面结合附图对本实实施例进行详细说明。如图1所示,本实施例的煤岩界面折射波的多阶折射时间求取方法包括:
步骤1,折射波周期计算步骤。由已知的围岩的地震波速度vr、煤层的地震波速度vc,煤层平均厚度h,计算出折射波周期T:
步骤2,由已知的激发点位置(xs,ys)和接收点位置(xr,yr),计算折射波的折射阶数K:
步骤3,估计震源主频fc。
选取激发点附近的折射波信号s0(t),其中t表示时间,t∈[0,NΔt],N为采集的折射波信号的长度,Δt是采样间隔。用傅里叶变换求出其振幅谱,寻找到最大振幅对应的频率,就是震源主频fc。
步骤4,构建主频为fc的雷克子波w0(τ)。
其中τ表示时间,τ∈[-T,T],T为折射波周期。
步骤5,估计出周期相位旋转角Φ,即折射波信号每经过一个周期,其相位发生旋转的角度,具体包括:
步骤5.1,求出折射波的初至时间t0:
其中:
H(·)为希尔伯特变换,τ∈[-T,T]。
步骤6,产生一个N×N的字典矩阵W,N为信号长度,包括以下步骤:
步骤6.1令n=0,w为空集;
步骤6.2,产生子波wn(t):
其中t∈[0,NΔt]。
步骤6.3,将子波wn(t)记为N维向量wn;
步骤6.4,令n=n+1,如果n<N则返回至第6.2步;
步骤6.5,将所有wn合并成一个矩阵W:
其中符号∪表示集合并运算。
步骤7,将折射阶数K作为输入参数,采用正交匹配追踪算法用字典矩阵W对折射信号s(t)进行稀疏反褶积,得到稀疏信号r(t)。具体实施为:
将长度为N的折射波信号s(t)记为N维向量s,设k为当前迭代次数,rk表示第k次迭代的稀疏向量,uk表示第k次迭代的残差向量,ik表示第k次迭代找到的索引(W矩阵中对应的列序号),aik表示矩阵W中以ik为序号的列向量,Ak为一个由所有挑选出的列向量aik所组成的集合矩阵。符号∪表示集合并运算,符号<·,·>表示向量内积。按下述步骤迭代:
步骤7.1,初始化,令k=1,A0为空集,u1=s;
步骤7.3,Ak=Ak-1∪aik;
步骤7.5,更新残差向量uk+1=s-Akrk;
步骤7.6,令k=k+1,如果k≤K则返回至第(2)步,否则将当前向量rk作为稀疏信号r(t)输出;
步骤8,稀疏信号r(t)具有K个脉冲,拾取这K个脉冲峰值对应的时间,得到折射波的多阶折射时间。
下面以正演模型数据为例说明本发明效果:
模型由顶板、煤层、底板三层组成,顶板与底板岩性相同。围岩纵波速度4000m/s,围岩横波速度2300m/s,密度2.56g/cm3;煤层纵波速度2000m/s,煤层横波速度1050m/s,密度1.4g/cm3。煤层厚度为10m,煤层中包含两条平行巷道,两巷间距200m。在煤层中存在一处薄煤区,大小为100m×50m,煤厚为5m。图2(a)为模型及观测系统平面图,图2(b)为模型剖面图。
激发点和接收点分别布置在两条巷道内,位置如图2(a)中所示。采用不同主频的雷克子波作为震源进行三维弹性波数值模拟,200Hz和400Hz主频震源产生的折射波信号分别如图3(a)和3(b)所示。图3(c)为R15接收到的两种频率的折射波信号的对比图,可以看出400Hz信号中各阶折射波彼此分隔,能直观的分辨出各阶折射波并从中找到最大振幅对应的各阶折射时间;而200Hz折射波则相互混叠,不能直观分辨出各阶折射波波形,因此无法找到各阶折射时间。
现依据本发明步骤对R15接收的折射波信号s(t)进行处理,求出各阶折射时间。
执行步骤1,由模型参数vr=4000m/s、vc=2000m/s,h=10m,求出折射波周期T=86.6ms。
执行步骤2,由激发点位置(xs,ys)=(400,55)和接收点R15位置(xr,yr)=(140,250),求得K=5。
执行步骤3,震源附近折射波信号s0(t),及其频谱分别如图4(a)和4(b)所示,得到主频fc为200Hz。
执行步骤4,以200Hz为主频构建雷克子波w0,如图5所示。
执行步骤5(1),计算出折射波初至时间t0=71ms。
执行步骤5(3),由图5估计出周期相位旋转角Φ=1.24。图7为R15道折射波信号s(t)与用Φ对子波w0进行周期相位旋转得到的合成信号的对比图。
执行步骤6,合成字典矩阵W。由于截取的折射波信号长度为1600个采样点,因此产生的W矩阵尺寸为1600×1600。将W中各元素作为图像中各点,则W矩阵可输出为一个1600×1600的图像,如图8所示,横坐标为矩阵的行,纵坐标为矩阵的列,图像中的值为矩阵中元素对应的颜色(例如灰度值)。
执行步骤7,对折射波信号s(t)进行如下处理:
执行步骤7(1),完成初始化,令k=1,将A0设为空集,将信号s(t)中数据设为u1;
执行步骤7(2),将W中各列与u1求内积,找到使内积最大的列为第1108列,则ik=1108;
执行步骤7(3),将A0与W中的第1108列合并成A1;
执行步骤7(4),计算出r1;
执行步骤7(5),计算出u1;
执行步骤7(7),由于k≤K,则令k=k+1,返回步骤7(2),直到k>K;s(t)与得到的稀疏信号r(t)如图9所示。
执行步骤8,从图9中的稀疏信号中,提取出各个最大值点对应的时间,分别为(76,849,938,102.4,1108,1174)ms。
可以看出虽然各阶折射波之间发生了混叠,但采用本发明获得的各阶折射时间相差基本为86ms,即一个折射波周期。因此,本发明能够准确获得的各阶折射时间。对所有数据逐道按上述步骤处理,得到的稀疏信号如图10所示,对比图3a,可以很明显看出右侧圈出区域内各阶折射时间存在异常。这种异常是由模型中的薄煤区引起的。
总体而言,本发明利用煤岩界面上传播的折射波的周期性对其进行稀疏反褶积,并从稀疏信号中获取多阶折射时间。准确的多阶折射时间是圈定煤厚异常区,并估计异常区内煤层厚度的先决条件,对保障煤矿工作面的智能化开采,提高回采效率和采出率,减少资源浪费具有重要意义。
本实施例中,尽管为使解释简单化将上述方法图示并描述为一系列动作,但是应理解并领会,这些方法不受动作的次序所限,因为根据一个或多个实施例,一些动作可按不同次序发生和/或与来自本文中图示和描述或本文中未图示和描述但本领域技术人员可以理解的其他动作并发地发生。
注意到,说明书中对“一个实施例”、“实施例”、“示例实施例”、“一些实施例”等的引用指示所描述的实施例可以包括特定特征、结构或特性,但是每个实施例可以不必包括所述特定特征、结构或特性。而且,这样的短语不必指代同一实施例。此外,当结合实施例描述特定特征、结构或特性时,无论是否明确描述,结合其他实施例来实现这样的特征、结构或特性将在所属领域的技术人员的知识范围内。
提供对本公开的先前描述是为使得本领域任何技术人员皆能够制作或使用本公开。对本公开的各种修改对本领域技术人员来说都将是显而易见的,且本文中所定义的普适原理可被应用到其他变体而不会脱离本公开的精神或范围。由此,本公开并非旨在被限定于本文中所描述的示例和设计,而是应被授予与本文中所公开的原理和新颖性特征相一致的最广范围。
Claims (7)
1.一种煤岩界面折射波的多阶折射时间求取方法,其特征在于,包括:
计算折射波的折射阶数K,并将子波wn(t)合成为字典矩阵W;
按以下迭代步骤计算得到稀疏信号r(t):
迭代步骤1,初始化,令k=1,A0为空集,u1=s;
迭代步骤3,基于下式求取迭代后的Ak:
Ak=Ak-1∪aik;
迭代步骤5,更新残差向量uk+1=s-Akrk;
迭代步骤6,令k=k+1,如果k≤K则返回至迭代步骤2,否则将当前向量rk作为稀疏信号r(t)输出;
迭代步骤7,拾取稀疏信号r(t)的K个脉冲峰值对应的时间,得到折射波的多阶折射时间;
其中,各迭代步骤中,k为当前迭代次数,rk表示第k次迭代的稀疏向量,uk表示第k次迭代的残差向量,ik表示第k次迭代找到的索引,所述索引为W矩阵中对应的列序号,aik表示矩阵W中以ik为序号的列向量,Ak为一个由所有挑选出的列向量aik所组成的集合矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210907743.1A CN115437006A (zh) | 2022-07-29 | 2022-07-29 | 一种煤岩界面折射波的多阶折射时间求取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210907743.1A CN115437006A (zh) | 2022-07-29 | 2022-07-29 | 一种煤岩界面折射波的多阶折射时间求取方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115437006A true CN115437006A (zh) | 2022-12-06 |
Family
ID=84242953
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210907743.1A Pending CN115437006A (zh) | 2022-07-29 | 2022-07-29 | 一种煤岩界面折射波的多阶折射时间求取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115437006A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117092700A (zh) * | 2023-08-25 | 2023-11-21 | 徐州恒佳机械科技有限公司 | 一种煤岩界面识别方法 |
CN118091753A (zh) * | 2023-11-21 | 2024-05-28 | 中南大学 | 基于稀疏反褶积与时频分析的强干扰环境震源初至拾取方法、装置、设备及介质 |
-
2022
- 2022-07-29 CN CN202210907743.1A patent/CN115437006A/zh active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117092700A (zh) * | 2023-08-25 | 2023-11-21 | 徐州恒佳机械科技有限公司 | 一种煤岩界面识别方法 |
CN117092700B (zh) * | 2023-08-25 | 2024-05-17 | 徐州恒佳机械科技有限公司 | 一种煤岩界面识别方法 |
CN118091753A (zh) * | 2023-11-21 | 2024-05-28 | 中南大学 | 基于稀疏反褶积与时频分析的强干扰环境震源初至拾取方法、装置、设备及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115437006A (zh) | 一种煤岩界面折射波的多阶折射时间求取方法 | |
US20190113641A1 (en) | Seismic travel time tomographic inversion method based on two point ray tracing | |
US7120541B2 (en) | Sonic well logging methods and apparatus utilizing parametric inversion dispersive wave processing | |
AU2008239658B2 (en) | Inverse-vector method for smoothing dips and azimuths | |
US9720130B2 (en) | Guided bayesian experimental design | |
NO334218B1 (no) | Behandling av målinger på lydbølgeformer for å bestemme langsomheten | |
CN107290722B (zh) | 微震源的定位方法和装置 | |
EP2713185B1 (en) | Method and apparatus to detect and analyze seismic signals | |
CN103424777A (zh) | 一种提高地震成像分辨率的方法 | |
US20170023687A1 (en) | Fracture Surface Extraction from Image Volumes Computed from Passive Seismic Traces | |
WO2014084751A1 (en) | A method for processing acoustic waveforms | |
US10459102B2 (en) | Signal detection in semblance methods | |
CN112305591B (zh) | 隧道超前地质预报方法、计算机可读存储介质 | |
CN103048690B (zh) | 基于最优地震子波提取的快速匹配投影分解的地层反射拾取技术 | |
Zhang et al. | Microseismic hydraulic fracture imaging in the Marcellus Shale using head waves | |
Bruno et al. | Groundwater characterization in arid regions using seismic and gravity attributes: Al Jaww Plain, UAE | |
CN115903024A (zh) | 一种基于梯度下降法的横波分裂分析方法 | |
CN114791633A (zh) | 监测页岩气压裂的方法、系统及介质 | |
WO2019224354A1 (en) | Method for establishing position and timing of seismic recording devices that are deployed on a seafloor | |
Eaton et al. | Enhanced downhole microseismic processing using matched filtering analysis (MFA) | |
WO2014084752A1 (en) | A method for processing acoustic waveforms | |
US11927711B2 (en) | Enhanced-resolution sonic data processing for formation body wave slowness with full offset waveform data | |
CN110858005B (zh) | 一种基于基追踪横向多道约束的各向异性参数反演方法 | |
Cheng et al. | 3D trajectory inversion of an adjacent well using scattered P-wave | |
Wang et al. | Broadband dispersion extraction of borehole acoustic modes via sparse Bayesian learning |
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 |