CN113296154A - 一种基伸缩调频同步挤压地震储层预测方法 - Google Patents

一种基伸缩调频同步挤压地震储层预测方法 Download PDF

Info

Publication number
CN113296154A
CN113296154A CN202110557394.0A CN202110557394A CN113296154A CN 113296154 A CN113296154 A CN 113296154A CN 202110557394 A CN202110557394 A CN 202110557394A CN 113296154 A CN113296154 A CN 113296154A
Authority
CN
China
Prior art keywords
frequency
time
seismic
signal
synchronous extrusion
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
Application number
CN202110557394.0A
Other languages
English (en)
Inventor
胡英
陈旭平
朱冰雪
陈辉
李蕊
周怀来
方玉霞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN202110557394.0A priority Critical patent/CN113296154A/zh
Publication of CN113296154A publication Critical patent/CN113296154A/zh
Pending legal-status Critical Current

Links

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. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures
    • 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. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • 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. analysis, for interpretation, for correction
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/40—Transforming data representation
    • G01V2210/45—F-x or F-xy domain
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/624—Reservoir parameters

Abstract

本发明公开了基伸缩调频同步挤压地震储层预测方法，包括以下步骤：S1、输入待分析的原始二维地震剖面信号s(x；t)；S2、利用基伸缩自适应调频变换对原始二维地震剖面信号s(x；t)进行分解，获取时频变换结果SBCT(x；f，t)；S3、根据上述时频变换结果SBCT(x；t，f)估计第x道地震信号在各个时频位置处的瞬时频率；S4、根据同步挤压原理，在该时频域上构建一种以信号瞬时频率为曲线中心的同步挤压算子，用于将原时频谱能量挤压至该曲线上获得新的时频系数；S5、对上述时频系数取模，得到基伸缩调频同步挤压变换后的时频谱。S6、采用傅氏变换计算原始地震信号的主频范围，确定高频值和低频值；再从S5得到的时频谱中分别抽取高频值和低频值对应的频谱，得到两个共频率剖面；最后通过对比两个共频率剖面中地震信号的衰减情况来确定该地震剖面中是否存在储层。本发明能够显著提高时频能量聚焦性，提高地震储层的识别精度。

Description

一种基伸缩调频同步挤压地震储层预测方法
技术领域
本发明属于信号处理领域，具体涉及一种基伸缩调频同步挤压地震储层预测方法。
背景技术
时频分析是地震数据处理的有力工具，它为我们提供了一种从二维角度解释地震信号时频特性的方法，即一维时间信号可以通过时频处理表示成时域和频域的联合分布函数此外，作为信号处理的有效手段，时频分析清楚地描述了信号频率和时间之间的紧密关系。时频分析方法作为地震储层预测的一种重要手段，能有效表征地震信号的非平稳特征，揭示信号频率随时间的变化关系。
但大多数时频表征方法在没有先验条件的情况下，都无法处理复杂的多分量信号。而基伸缩自适应调频变换(SBCT)方法通过构造基函数从而自适应匹配多分量信号。通过在一定的窗口长度范围内对时间中心及其周围的基函数进行缩放来扩展chirp率，从而生成随时间和频率变化的线性调频，使得整个窗口长度范围内的多分量信号都能匹配每个瞬时频率轨迹的斜率。与现有方法相比，SBCT方法所获得的时频表征可以实现更高的时频能量聚集性。
基伸缩调频同步挤压地震储层预测方法的基本思想是先通过SBCT获取信号的时频表征结果，再在该时频域上构建一种估计信号真实瞬时频率的同步挤压算子用于对原时频谱能量进行挤压，将时频能量重新排列至信号真实的瞬时频率脊线上以显著提高信号时频表征能量分布的聚焦性。
发明内容
针对现有技术中的上述不足，本发明提供的一种基伸缩调频同步挤压地震储层预测方法。本发明能够显著提高时频能量聚焦性，提高地震储层的识别精度。
为了达到上述发明目的，本发明采用的技术方案为：一种基伸缩调频同步挤压地震储层预测方法，包括以下步骤：
S1、输入待分析的原始二维地震剖面信号s(x；t)；
S2、利用基伸缩自适应调频变换对原始二维地震剖面信号s(x；t)进行分解，获取时频变换结果SBCT(x；f，t)；
S3、根据上述得到的时频变换结果SBCT(x；f，t)估计第x道地震信号在各个时频位置处的瞬时频率；
S4、根据同步挤压原理，在该时频域上构建一种以信号瞬时频率为曲线中心的同步挤压算子用于将原时频谱能量挤压至该曲线上获得新的时频系数；
S5、对上述时频系数取模，得到基伸缩调频同步挤压变换后的时频谱；
S6、采用傅氏变换计算原始地震信号的主频范围，将主频范围中的频率最大值的85％设置为高频值，将主频范围中的频率最小值设置为低频值，然后从S5得到的时频谱中抽取高频值对应的频谱，得到高频值的共频率剖面；从S5得到的时频谱中抽取低频值对应的频谱，得到低频值的共频率剖面。通过对比两个共频率剖面中地震信号的衰减情况来确定该地震剖面中是否存在储层。
作为优选，所述步骤S2中原始二维地震剖面信号s(x；t)的基伸缩自适应调频变换的时频变换结果为：
Figure BDA0003077674470000021
其中，t表示时间，f表示频率中心，h(x；τ)表示高斯函数。
Figure BDA0003077674470000022
表示相位函数，τ为时间变量，所述相位函数
Figure BDA0003077674470000031
为：
Figure BDA0003077674470000032
其中，k＝1，2，...，n，n为相位函数
Figure BDA0003077674470000033
的阶数，(a1，a2，...，an)是相位函数的自适应参数，可根据峰度理论确定，即：
Figure BDA0003077674470000034
其中，argmax(·)表示函数取到最大值时的参数集。
作为优选，根据得到的时频变换结果SBCT(x；t，f)估计信号在各个时频位置(t，f)处的瞬时频率：
Figure BDA0003077674470000035
其中，
Figure BDA0003077674470000036
表示窗函数h(x；t)对t的偏导，
Figure BDA0003077674470000037
表示在窗函数
Figure BDA0003077674470000038
下的SBCT结果；
Figure BDA0003077674470000039
良示在函数tkh(x；t)下的SBCT结果。
作为优选，所述步骤S4中获取新的时频系数的具体方法为：根据同步挤压原理，在该时频域上构建一种以信号瞬时频率为曲线中心的同步挤压算子STO(x；f，t)用于将原时频谱能量挤压至该曲线上获得新的时频系数T(x；f，t)，同步挤压算子STO(x；f，t)满足：
Figure BDA00030776744700000310
其中，δ(·)为单位脉冲函数。新的时频系数T(x；f，t)为：
Figure BDA00030776744700000311
作为优选，所述步骤S5可采用下式对(6)中的T(x；f，t)进行逆变换，重构地震信号s(x；t).
Figure BDA00030776744700000312
其中
Figure BDA0003077674470000041
本发明的思路为：
首先，输入待分析的原始二维地震剖面信号s(x；t)；
第二，对输入信号s(x；t)进行基伸缩自适应调频变换，获取时频变换结果SBCT(x；t，f)
Figure BDA0003077674470000042
其中，t表示时间，f表示频率中心，h(x；τ)表示高斯函数，
Figure BDA0003077674470000043
表示相位函数，τ为时间变量；
第三，根据得到的时频变换结果SBCT(x；t，f)估计第x道地震信号在各个时频位置(t，f)处的瞬时频率
Figure BDA0003077674470000044
第四，根据同步挤压原理，在该时频域上构建一种以信号瞬时频率曲线为中心的同步挤压算子STO(x；f，t)用于将原时频谱能量挤压至该曲线上获得新的时频系数；
第五，对第四步获得的时频系数取模，得到基伸缩调频同步挤压变换后的时频谱。
第六，采用傅氏变换计算原始地震信号的主频范围，将主频范围中的频率最大值的85％设置为高频值，将主频范围中的频率最小值设置为低频值，然后从第五步得到的时频谱中抽取高频值对应的频谱，得到高频值的共频率剖面；从第五步得到的时频谱中抽取低频值对应的频谱，得到低频值的共频率剖面。通过对比两个共频率剖面中地震信号的衰减情况来确定该地震剖面中是否存在储层。
本发明的工作原理为：采集待分析的原始二维地震剖面信号s(x；t)；利用基伸缩自适应调频变换对原始二维地震剖面信号s(x；t)进行分解，获取时频变换结果SBCT(x；f，t)；根据上述时频变换结果SBCT(x；f，t)估计第x道信号在各个时频位置处的瞬时频率；根据同步挤压原理，在该时频域上构建一种以信号瞬时频率为曲线中心的同步挤压算子用于将原时频谱能量挤压至该曲线上获得新的时频系数；对上述时频系数取模，得到基伸缩调频同步挤压变换后的时频谱。本发明能够显著提高时频能量聚焦性，提高地震储层的识别精度。
本发明针对同步挤压变换和基伸缩自适应调频变换在多分量非平稳非线性复杂信号中出现的识别不够清晰、聚集性不好等问题，提出了基伸缩调频同步挤压地震储层预测方法，首先对信号进行基伸缩自适应调频变换处理得到其表征结果，然后根据时频谱的相位信息估计信号在各个时频位置处的瞬时频率。再根据同步挤压原理，在该时频域上构建一种估计信号真实瞬时频率的同步挤压算子用于对原时频谱能量进行挤压，获得新的时频系数；最后对上述时频系数取模，得到基伸缩调频同步挤压变换后的时频谱。本发明在重构能力以及计算效率上都表现出良好的效果，能够显著提高时频能量聚焦性，提高地震储层的识别精度。
附图说明
图1为本发明流程图；
图2为中南西部地区四川盆地中江气田剖面图
图3为经本发明方法处理的中江气田的低频值共频率剖面
图4为经本发明方法处理的中江气田的高频值共频率剖面
具体实施方式
下面将结合附图对本发明作进一步说明。
实施例1：参见图1，一种基伸缩调频同步挤压地震储层预测方法，包括以下步骤：
S1、输入待分析的原始二维地震剖面信号s(x；t)；
S2、利用基伸缩自适应调频变换对原始二维地震剖面信号s(x；t)进行分解，获取原始二维地震剖面信号s(x；t)的基伸缩自适应调频变换的时频变换结果；
S3、根据上述时频变换结果SBCT(x；t，f)估计信号在各个时频位置处的瞬时频率；
S4、根据同步挤压原理，在该时频域上构建一种以信号瞬时频率为曲线中心的同步挤压算子STO(x；f，t)用于将原时频谱能量挤压至该曲线上获得新的时频系数；
S5、对上述时频系数取模，得到基伸缩调频同步挤压变换后的时频谱。
S6、采用傅氏变换计算原始地震信号的主频范围，将主频范围中的频率最大值的85％设置为高频值，将主频范围中的频率最小值设置为低频值，然后从S5得到的时频谱中抽取高频值对应的频谱，得到高频值的共频率剖面；从S5得到的时频谱中抽取低频值对应的频谱，得到低频值的共频率剖面。通过对比两个共频率剖面中地震信号的衰减情况来确定该地震剖面中是否存在储层。
作为优选，所述步骤S2中二维地震剖面信号s(x；t)的基伸缩自适应调频变换的时频谱为：
Figure BDA0003077674470000061
其中，t表示时间，f表示频率中心，h(x；τ)表示高斯函数，
Figure BDA0003077674470000062
表示相位函数，τ为时间变量；
所述相位函数
Figure BDA0003077674470000071
为：
Figure BDA0003077674470000072
其中，k＝1，2，...，n，n为相位函数
Figure BDA0003077674470000073
的阶数，(a1，a2，...，an)是相位函数的自适应参数，可根据峰度理论确定，即：
Figure BDA0003077674470000074
其中，argmax(·)表示函数取到最大值时的参数集。
作为优选，所述步骤S3中根据时频谱相位信息估计信号在各个时频位置(t，f)处的瞬时频率为：
Figure BDA0003077674470000075
其中，
Figure BDA0003077674470000076
表示窗函数h(x；t)对t的偏导，
Figure BDA0003077674470000077
表示在窗函数
Figure BDA0003077674470000078
下的SBCT结果；
Figure BDA0003077674470000079
表示在函数tkh(x；t)下的SBCT结果。
作为优选，所述步骤S4中获取新的时频系数的具体方法为：根据同步挤压原理，在该时频域上构建一种以信号瞬时频率为曲线中心的同步挤压算子STO(x；f，t)用于将原时频谱能量挤压至该曲线上获得新的时频系数T(x；f，t)，同步挤压算子STO(x；f，t)满足：
Figure BDA00030776744700000710
其中，δ(·)为单位脉冲函数。新的时频系数T(x；f，t)为：
Figure BDA00030776744700000711
作为优选，所述步骤S5可采用下式对(6)中的T(x；t，η)进行逆变换，重构地震信号s(x；t)；
Figure BDA0003077674470000081
其中
Figure BDA0003077674470000082
本发明适用于重构地震信号。
参见图1到图4，我们具体以一个地震剖面为例，该二维地震剖面图如图2所示。图3-图4分别是从地震剖面进行基伸缩自适应调频变换得到的时频谱中抽取的低频值共频率剖面和高频值共频率剖面。图中，横坐标表示地震道数，纵坐标表示时间，右侧色标表示能量值。实施例证明，经本发明方法处理后得到的结果图时频分辨率更高，能量更加集中，提高了地震储层的识别精度。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (6)

1.一种基伸缩调频同步挤压地震储层预测方法，其特征在于，包括以下步骤：
S1、输入待分析的原始二维地震剖面信号s(x；t)；其中x表示地震道，t表示时间；
S2、对对S1输入的二维地震剖面信号s(x；t)进行基伸缩自适应调频变换，获取时频变换结果SBCT(x；f,t)，计算方法如下：
Figure FDA0003077674460000011
其中，t表示时间中心，f表示频率中心，h(x；τ)表示高斯函数，
Figure FDA0003077674460000012
表示相位函数，τ为时间变量；
S3、根据S2得到的时频SBCT(x；t,f)估计第x道地震信号在各个时频位置(t,f)处的瞬时频率
Figure FDA0003077674460000013
S4、根据同步挤压原理，在该时频域上构建一种以信号瞬时频率曲线为中心的同步挤压算子STO(x；t,f)用于将原时频谱能量挤压至该曲线上获得新的时频系数T(x；f,t)；
S5、对S4获得的时频系数取模，得到基伸缩调频同步挤压变换后的时频谱；
S6、采用傅氏变换计算原始地震信号的主频范围，将主频范围中的频率最大值的85％设置为高频值，将主频范围中的频率最小值设置为低频值，然后从S5得到的时频谱中抽取高频值对应的频谱，得到高频值的共频率剖面；从S5得到的时频谱中抽取低频值对应的频谱，得到低频值的共频率剖面。通过对比两个共频率剖面中地震信号的衰减情况来确定该地震剖面中是否存在储层。
2.根据权利要求1所述的基伸缩调频同步挤压地震储层预测方法，其特征在于，所述步骤S2中二维地震剖面信号s(x；t)的基伸缩自适应调频变换的时频谱为：
Figure FDA0003077674460000021
其中，t表示时间中心，f表示频率中心，h(x；τ)表示高斯函数，
Figure FDA0003077674460000022
表示相位函数，τ为时间变量，所述相位函数
Figure FDA0003077674460000023
为：
Figure FDA0003077674460000024
其中，k＝1,2,…,n，n为相位函数
Figure FDA0003077674460000025
的阶数，(a1,a2,…,an)是相位函数
Figure FDA0003077674460000026
的自适应参数，可根据峰度理论确定，即：
Figure FDA0003077674460000027
其中，argmax(·)表示函数取到最大值时的参数集。
3.根据权利要求1所述的基伸缩调频同步挤压地震储层预测方法，其特征在于，所述步骤S3中根据S2得到的时频变换结果SBCT(x；f,t)估计信号在各个时频位置处(t,f)的瞬时频率
Figure FDA0003077674460000028
Figure FDA0003077674460000029
其中，
Figure FDA00030776744600000210
表示窗函数h(x；t)对时间t的偏导，
Figure FDA00030776744600000211
表示在窗函数
Figure FDA00030776744600000212
下的SBCT结果；
Figure FDA00030776744600000213
表示在函数tkh(x；t)下的SBCT结果。
4.根据权利要求1所述的基伸缩调频同步挤压地震储层预测方法，其特征在于，所述步骤S4中获取新的时频系数T(x；f,t)的具体方法为：根据同步挤压原理，在该时频域上构建一种以信号瞬时频率曲线为中心的同步挤压算子STO(x；t,f)用于对原时频谱能量挤压至该曲线上获得新的时频系数T(x；f,t)；其中，同步挤压算子STO(x；t,f)满足：
Figure FDA0003077674460000031
其中δ(·)为单位脉冲函数，新的时频系数T(x；f,t)：
Figure FDA0003077674460000032
5.根据权利要求1所述的基伸缩调频同步挤压地震储层预测方法，其特征在于，对所述步骤S4中获得的时频系数取模。
6.根据权利要求1所述基伸缩调频同步挤压地震储层预测方法，其特征在于，所述步骤S5可采用下式对(7)中的T(x；t,η)进行逆变换，重构地震信号s(x；t)；
Figure FDA0003077674460000033
其中
Figure FDA0003077674460000034
本发明适用于重构地震信号。
CN202110557394.0A 2021-05-21 2021-05-21 一种基伸缩调频同步挤压地震储层预测方法 Pending CN113296154A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110557394.0A CN113296154A (zh) 2021-05-21 2021-05-21 一种基伸缩调频同步挤压地震储层预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110557394.0A CN113296154A (zh) 2021-05-21 2021-05-21 一种基伸缩调频同步挤压地震储层预测方法

Publications (1)

Publication Number Publication Date
CN113296154A true CN113296154A (zh) 2021-08-24

Family

ID=77323600

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110557394.0A Pending CN113296154A (zh) 2021-05-21 2021-05-21 一种基伸缩调频同步挤压地震储层预测方法

Country Status (1)

Country Link
CN (1) CN113296154A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114252915A (zh) * 2021-11-03 2022-03-29 成都理工大学 一种基于二阶水平多重同步挤压变换的油气储层识别方法
CN114280671A (zh) * 2021-09-17 2022-04-05 成都理工大学 一种广义w变换地震含气储层表征方法
CN114355442A (zh) * 2022-01-12 2022-04-15 成都理工大学 一种三参数w变换的地震储层识别时频分析方法
CN114563824A (zh) * 2022-02-25 2022-05-31 成都理工大学 一种二阶多重同步挤压多项式chirplet变换薄储层识别方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114280671A (zh) * 2021-09-17 2022-04-05 成都理工大学 一种广义w变换地震含气储层表征方法
CN114252915A (zh) * 2021-11-03 2022-03-29 成都理工大学 一种基于二阶水平多重同步挤压变换的油气储层识别方法
CN114355442A (zh) * 2022-01-12 2022-04-15 成都理工大学 一种三参数w变换的地震储层识别时频分析方法
CN114563824A (zh) * 2022-02-25 2022-05-31 成都理工大学 一种二阶多重同步挤压多项式chirplet变换薄储层识别方法
CN114563824B (zh) * 2022-02-25 2024-01-30 成都理工大学 一种二阶多重同步挤压多项式chirplet变换薄储层识别方法

Similar Documents

Publication Publication Date Title
CN113296154A (zh) 一种基伸缩调频同步挤压地震储层预测方法
CN106597532B (zh) 一种结合井资料与层位资料的叠前地震数据频带拓展方法
Chen et al. Comparison of methods for non-stationary hydrologic frequency analysis: Case study using annual maximum daily precipitation in Taiwan
CN110361778B (zh) 一种基于生成对抗网络的地震数据重建方法
CN109343020B (zh) 一种基于改进窗函数的s变换时频分析方法
CN111856562B (zh) 一种广义高阶同步挤压地震信号时频分解与重构方法
CN107083957B (zh) 钻井液随钻信号的泵冲干扰消除方法及系统
CN113296155A (zh) 一种基伸缩调频同步提取地震储层预测方法
CN110503060B (zh) 一种光谱信号去噪方法及其系统
CN103345920B (zh) 基于Mel-KSVD稀疏表示的自适应内插加权谱模型的语音转换及重构方法
CN106707334A (zh) 一种提高地震资料分辨率的方法
CN109100687B (zh) 一种雷达设备lfm脉冲信号pslr参数确定方法
Li Speech denoising based on improved discrete wavelet packet decomposition
CN112101089B (zh) 信号降噪方法、装置、电子设备及存储介质
Han et al. Noise reduction method for chaotic signals based on dual-wavelet and spatial correlation
CN103915102A (zh) 一种lfm水声多途信号的噪声抑制方法
CN111474582A (zh) 生成高精度时频谱的精准s变换方法
CN111198357A (zh) 一种基于可调窗函数的s变换时频分析方法
CN114638267A (zh) 基于爆破振动多元特征信息比对的电子雷管盲炮识别方法
CN109270573B (zh) 一种快速保频保幅s变换方法
CN109597042B (zh) 基于奇异谱分析的目标进动频率估计方法
Dadu et al. Multistep lead time forecasting of hydrologic time series using daubechies wavelet–neural network hybrid model
Huang Wave parameters and functions in wavelet analysis with filtering
CN112505640A (zh) 基于参数自适应的扩展b分布脉冲信号时频分析方法
CN103854658B (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