CN115826056A - 深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法 - Google Patents

深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法 Download PDF

Info

Publication number
CN115826056A
CN115826056A CN202310132204.XA CN202310132204A CN115826056A CN 115826056 A CN115826056 A CN 115826056A CN 202310132204 A CN202310132204 A CN 202310132204A CN 115826056 A CN115826056 A CN 115826056A
Authority
CN
China
Prior art keywords
depth
seismic
streamer
towed
deep
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
Application number
CN202310132204.XA
Other languages
English (en)
Other versions
CN115826056B (zh
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.)
Shandong University of Science and Technology
Original Assignee
Shandong University of Science and 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 Shandong University of Science and Technology filed Critical Shandong University of Science and Technology
Priority to CN202310132204.XA priority Critical patent/CN115826056B/zh
Publication of CN115826056A publication Critical patent/CN115826056A/zh
Application granted granted Critical
Publication of CN115826056B publication Critical patent/CN115826056B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Geophysics And Detection Of Objects (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明公开深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法，属于地球物理技术领域，用于拖缆水听器阵列的定位，包括使用深拖地震数据的海面反射波或海底反射波初至时间利用海水平均速度计算水听器阵列的初始深度，根据拖缆水听器阵列为曲线的特点，将水听器阵列初始深度与震源组合，采用加权最小二乘法多项式拟合的方法组合后的拟合优度或均方根误差，通过更新计算时的海水速度使得拟合优度最大或均方根误差最小，最后将拟合结果作为最终水听器阵列的各道深度值。本发明能够获得定位精度达到小于一个采样间隔的深拖高分辨率地震数据处理的要求，考虑了拖缆阵列的曲率平滑，为后续深拖地震数据精细处理提供准确的各道深度。

Description

深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法
技术领域
本发明公开深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法，属于地球物理技术领域。
背景技术
拖缆水听器阵列的高精度定位是深拖式高分辨率地震探测技术及后续数据处理中最重要的问题之一。深拖式高分辨率地震探测系统的震源——拖缆水听器阵列由于频繁收放缆（保持在距海底一定高度或保障系统安全，需要收放电缆）、拖曳速度的不稳定和深海环境变化等因素的影响，震源和拖缆水听器阵列在数据采集过程中不可避免的上下浮动。对室内高分辨率地震数据处理而言，拖缆水听器阵列定位精度的高低，是数据处理的结果质量的关键所在，一般要求垂向上的精度达到一个采样间隔Δt。由于深拖高分辨率地震探测系统的震源和拖缆水听器阵列作业时拖曳于海平面以下300 m，所以基于电磁波的无线电导航系统GPS无法对深拖震源和拖缆水听器阵列进行定位。目前深拖拖缆水听器阵列定位主要采用走时定位的方法，然而在计算时采用的海水水体的平均速度，导致计算结果达不到数据处理的要求，且没有考虑走时的相位畸变和拖缆水听器阵列的曲率平滑。
发明内容
本发明的目的在于提供深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法，以解决现有技术中，深拖拖缆水听器阵列定位的计算结果达不到数据处理的要求，且没有考虑走时的相位畸变和拖缆水听器阵列的曲率平滑的问题。
深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法，包括：
S1.计算获取拖缆水听器阵列的初水平位置和垂直深度，获得拖缆水听器初始坐标
Figure SMS_1
；
S2.确定海水水体速度的范围；
S3.将震源深度与S1的拖缆水听器初始坐标组成数据点集，利用加权最小二乘法多项式拟合方法，对数据点集进行拟合，并求均方根误差或拟合优度，给走时相位畸变的点赋予小的权重，对于畸变严重的点权重给0；
S4.取拟合优度最大或均方根误差最小的f(xi)为拟合后的拖缆水听器阵列修正的最终深度值。
S1包括：
计算拖缆水听器阵列的两项参数：
初水平位置
Figure SMS_2
:
Figure SMS_3
(1)；
垂直深度
Figure SMS_4
:
Figure SMS_5
(2)；
其中，
Figure SMS_6
为地震记录中拾取的第i道直达波的走时值，
Figure SMS_7
为海面反射波或海底反射波的走时值，
Figure SMS_8
为地震波在测线区域海水中传播的速度，
Figure SMS_9
为震源深度，即拖体深度，
Figure SMS_10
为地震波系统深度的海水声波速度，由温盐深测量仪测得。
S3包括：
S3.1.将震源坐标(0，
Figure SMS_11
)与（
Figure SMS_12
，
Figure SMS_13
）组成数据点集[(0，
Figure SMS_14
)，
Figure SMS_15
]；
S3.2.利用加权最小二乘法多项式拟合方法对数据点集进行拟合：
Figure SMS_16
（3）；
其中，
Figure SMS_17
为多项式f(x)的系数，令偏差
Figure SMS_18
与权重
Figure SMS_19
的乘积的平方的和最小：
Figure SMS_20
（4）；
式中，
Figure SMS_21
表示第i个
Figure SMS_22
；
S3.3.确定多项式f(x)中的系数
Figure SMS_23
，根据最小二乘原则：
Figure SMS_24
（5）；
求出系数
Figure SMS_25
后得到拟合方程，给相位走时相位畸变点计算得到的水听器坐标赋予较小的权重，
Figure SMS_26
=0；
S3.4.计算拟合优度R2和均方根误差S，评价拟合的质量：
Figure SMS_27
（6）；
Figure SMS_28
（7）；
其中
Figure SMS_29
是走时定位计算的水听器垂直深度
Figure SMS_30
，
Figure SMS_31
为走时定位计算的水听器垂直深度
Figure SMS_32
的平均值。
相对比现有技术，本发明具有以下有益效果：采用加权最小二乘法多项式拟合的方法，对走时定位结果进行了修正，解决了使用平均速度计算的拖缆水听器阵列带来的深度误差，并解决了个别道存在的相位畸变问题和拖缆曲率平滑问题；本发明最终定位精度优于0.15m，达到了深拖高分辨率阵列定位的精度要求；采用本方法计算的拖缆水听器阵列定位结果应用于深拖数据处理的浮动基准面处理后，速度普能量团聚焦的更好，提高了速度分析的精度，且叠加成像的信噪比和分辨率得到了有效的提高，能够提供更为详实的与海底沉积和构造细节。
附图说明
图1为海面或海底反射波走时定位原理图；
图2为海面反射波地震记录图；
图3为拖缆水听器阵列初始位置图；
图4为最小二乘法多项式拟合后与现有方法计算的初始位置对比图；
图5为不同海水速度计算拟合后的均方根误差变化图；
图6为不同海水速度计算拟合后的拟合优度变化图；
图7为常规方法拖缆水听器阵列深度误差图；
图8为本发明方法的拖缆水听器阵列深度误差图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法，包括：
S1.计算获取拖缆水听器阵列的初水平位置和垂直深度，获得拖缆水听器初始坐标
Figure SMS_33
；
S2.确定海水水体速度的范围；
S3.将震源深度与S1的拖缆水听器初始坐标组成数据点集，利用加权最小二乘法多项式拟合方法，对数据点集进行拟合，并求均方根误差或拟合优度，给走时相位畸变的点赋予小的权重，对于畸变严重的点权重给0；
S4.取拟合优度最大或均方根误差最小的f(xi)为拟合后的拖缆水听器阵列修正的最终深度值。
S1包括：
计算拖缆水听器阵列的两项参数：
初水平位置
Figure SMS_34
:
Figure SMS_35
(1)；
垂直深度
Figure SMS_36
:
Figure SMS_37
(2)；
其中，
Figure SMS_38
为地震记录中拾取的第i道直达波的走时值，
Figure SMS_39
为海面反射波或海底反射波的走时值，
Figure SMS_40
为地震波在测线区域海水中传播的速度，
Figure SMS_41
为震源深度，即拖体深度，
Figure SMS_42
为地震波系统深度的海水声波速度，由温盐深测量仪测得。
S3包括：
S3.1.将震源坐标(0，
Figure SMS_43
)与（
Figure SMS_44
，
Figure SMS_45
）组成数据点集[(0，
Figure SMS_46
)，
Figure SMS_47
]；
S3.2.利用加权最小二乘法多项式拟合方法对数据点集进行拟合：
Figure SMS_48
（3）；
其中，
Figure SMS_49
为多项式f(x)的系数，令偏差
Figure SMS_50
与权重
Figure SMS_51
的乘积的平方的和最小：
Figure SMS_52
（4）；
式中，
Figure SMS_53
表示第i个
Figure SMS_54
；
S3.3.确定多项式f(x)中的系数
Figure SMS_55
，根据最小二乘原则：
Figure SMS_56
（5）；
求出系数
Figure SMS_57
后得到拟合方程，给相位走时相位畸变点计算得到的水听器坐标赋予较小的权重，
Figure SMS_58
=0；
S3.4.计算拟合优度R2和均方根误差S，评价拟合的质量：
Figure SMS_59
（6）；
Figure SMS_60
（7）；
其中
Figure SMS_61
是走时定位计算的水听器垂直深度
Figure SMS_62
，
Figure SMS_63
为走时定位计算的水听器垂直深度
Figure SMS_64
的平均值。
本发明实施例中，通过更新海水水体的速度值，重新计算拖缆水听器阵列的深度值并回到S2，将拟合优度最大或均方根误差最小的拟合结果作为最终拖缆水听器阵列深度定位结果，将深拖震源和各道检波点置于一个起伏较为平缓的浮动基准面上，获取深拖地震数据浮动基准面校正量；
本实施例中，深拖式高分辨率多道地震探测系统拖缆水听器阵列海面或海底反射波走时定位原理如图1所示，其中，S为电火花震源，
Figure SMS_65
为拖缆水听器阵列的第i道，
Figure SMS_66
为拖缆水听器阵列的第i道相对海面或海底的镜像点；
Figure SMS_67
为第i道水听器的垂直深度；
Figure SMS_68
为第i道水听器的水平偏移距；L为海面或海底反射波走时距离。图2为某炮点的海面反射波地震记录，为浮动基准面校正前的速度谱，图3为利用公式（1）和（2）计算得到的拖缆水听器阵列初始位置图，为浮动基准面校正后的速度谱。经对比可以看出，经过浮动基准面校正后，速度谱能量团聚焦效果非常好，有利于精细速度分析。
实施例中
Figure SMS_69
的取值范围为1480~1490m/s，本实施例中选取的一炮利用计算的拖缆水听器阵列初始值。图4为最小二乘法多项式拟合后与初始位置（现有方法计算）对比图。拖缆水听器阵列不光滑且与震源不在同一曲线上，可见最优速度定位结果和拟合曲线均与震源在同一曲线上。
本实施例中选取完成了步骤S3的一炮数据，图5为均方根误差图，可见当海水水体速度为1485.315m/s时，均方根误差最小，约为0.1355，即拟合曲线最优。图6为拟合优度图，可见当速度为1485.315m/s时，拟合优度达到了0.9995，即拟合最优。
本实施例中选取完成了步骤S4的一炮数据，图7为常规方法拖缆水听器阵列深度误差图，图8为本发明方法的拖缆水听器阵列深度误差图，可见修真后的拖缆水听器阵列各道深度误差绝对值均小于0.15m(误差时间约为0.1ms)，小于系统数据的采样间隔0.125ms。定位修真后的速度谱能量团更聚焦，有利于提高速度分析，定位修真后的叠加剖面同相轴连续性较好，信噪比高，可见修正后的定位精度完全满足了深拖地震数据处理对拖缆水听器阵列定位精度的要求。
综上所述，本发明通过综合引入了加权最小二乘法多项式拟合的方法实现了深拖拖缆水听器阵列的高精度定位，避免定位误差、相位畸变和拖缆水听器阵列曲率平滑的问题，定位结果用于深拖地震数据后处理，能够实现海底以下地层构造形态的高分辨率、高信噪比、高保真度成像。
应当理解的是，本发明方法各个步骤的执行并有严格的顺序限制，这些步骤必须按照顺序执行。
以上实施例仅用于说明本发明的技术方案，而非对其限制，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换，而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (3)

1.深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法，其特征在于，包括：
S1.计算获取拖缆水听器阵列的初水平位置和垂直深度，获得拖缆水听器初始坐标（
Figure QLYQS_1
，
Figure QLYQS_2
）；
S2.确定海水水体速度的范围；
S3.将震源深度与S1的拖缆水听器初始坐标组成数据点集，利用加权最小二乘法多项式拟合方法，对数据点集进行拟合，并求均方根误差或拟合优度，给走时相位畸变的点赋予小的权重，对于畸变严重的点权重给0；
S4.取拟合优度最大或均方根误差最小的f(xi)为拟合后的拖缆水听器阵列修正的最终深度值。
2.根据权利要求1所述的深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法，其特征在于，S1包括：
计算拖缆水听器阵列的两项参数：
初水平位置
Figure QLYQS_3
:
Figure QLYQS_4
(1)；
垂直深度
Figure QLYQS_5
:
Figure QLYQS_6
(2)；
其中，
Figure QLYQS_7
为地震记录中拾取的第i道直达波的走时值，
Figure QLYQS_8
为海面反射波或海底反射波的走时值，
Figure QLYQS_9
为地震波在测线区域海水中传播的速度，
Figure QLYQS_10
为震源深度，即拖体深度，
Figure QLYQS_11
为地震波系统深度的海水声波速度，由温盐深测量仪测得。
3.根据权利要求2所述的深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法，其特征在于，S3包括：
S3.1.将震源坐标(0，
Figure QLYQS_13
)与（
Figure QLYQS_14
，
Figure QLYQS_15
）组成数据点集[(0，
Figure QLYQS_16
)，(
Figure QLYQS_17
，
Figure QLYQS_18
)… (
Figure QLYQS_19
，
Figure QLYQS_12
)]；
S3.2.利用加权最小二乘法多项式拟合方法对数据点集进行拟合：
Figure QLYQS_20
（3）；
其中，
Figure QLYQS_21
…
Figure QLYQS_22
为多项式f(x)的系数，令偏差
Figure QLYQS_23
与权重
Figure QLYQS_24
的乘积的平方的和最小：
Figure QLYQS_25
（4）；
式中，
Figure QLYQS_26
表示第i个
Figure QLYQS_27
；
S3.3.确定多项式f(x)中的系数
Figure QLYQS_28
，根据最小二乘原则：
Figure QLYQS_29
（5）；
求出系数
Figure QLYQS_30
后得到拟合方程，给相位走时相位畸变点计算得到的水听器坐标赋予较小的权重，
Figure QLYQS_31
=0；
S3.4.计算拟合优度R2和均方根误差S，评价拟合的质量：
Figure QLYQS_32
（6）；
Figure QLYQS_33
（7）；
其中
Figure QLYQS_34
是走时定位计算的水听器垂直深度
Figure QLYQS_35
，
Figure QLYQS_36
为走时定位计算的水听器垂直深度
Figure QLYQS_37
的平均值。
CN202310132204.XA 2023-02-20 2023-02-20 深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法 Active CN115826056B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310132204.XA CN115826056B (zh) 2023-02-20 2023-02-20 深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310132204.XA CN115826056B (zh) 2023-02-20 2023-02-20 深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法

Publications (2)

Publication Number Publication Date
CN115826056A true CN115826056A (zh) 2023-03-21
CN115826056B CN115826056B (zh) 2023-06-16

Family

ID=85521767

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310132204.XA Active CN115826056B (zh) 2023-02-20 2023-02-20 深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法

Country Status (1)

Country Link
CN (1) CN115826056B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5148406A (en) * 1989-05-31 1992-09-15 Geco A.S. Method for simultaneous collection of seismic data from shallow and deep targets
US5648937A (en) * 1995-01-18 1997-07-15 Atlantic Richfield Company Method and apparatus for correlating geological structure horizons from velocity data to well observations
CN101852867A (zh) * 2009-04-03 2010-10-06 中国石油集团东方地球物理勘探有限责任公司 一种矢量合成检波点二次定位方法
US20120147700A1 (en) * 2010-12-14 2012-06-14 Svein Arne Frivik Determining Streamer Depth and Sea Surface Profile
US20150003195A1 (en) * 2013-07-01 2015-01-01 Pgs Geophysical As Variable depth multicomponent sensor streamer
CN109765620A (zh) * 2019-02-20 2019-05-17 自然资源部第二海洋研究所 近底拖曳式随机接收电缆地震数据采集系统与方法
CN110780350A (zh) * 2019-10-25 2020-02-11 中国海洋大学 一种海底地震检波器动态定位方法、系统及检波器
CN112505750A (zh) * 2020-11-18 2021-03-16 广州海洋地质调查局 一种深拖多道地震拖缆姿态确定方法及处理终端
CN114460649A (zh) * 2022-04-14 2022-05-10 自然资源部第一海洋研究所 一种深海近底拖曳式多道地震接收阵列形态重建方法
CN114779340A (zh) * 2022-06-21 2022-07-22 自然资源部第一海洋研究所 一种深拖式高分辨率多道地震数据精细处理方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5148406A (en) * 1989-05-31 1992-09-15 Geco A.S. Method for simultaneous collection of seismic data from shallow and deep targets
US5648937A (en) * 1995-01-18 1997-07-15 Atlantic Richfield Company Method and apparatus for correlating geological structure horizons from velocity data to well observations
CN101852867A (zh) * 2009-04-03 2010-10-06 中国石油集团东方地球物理勘探有限责任公司 一种矢量合成检波点二次定位方法
US20120147700A1 (en) * 2010-12-14 2012-06-14 Svein Arne Frivik Determining Streamer Depth and Sea Surface Profile
US20150003195A1 (en) * 2013-07-01 2015-01-01 Pgs Geophysical As Variable depth multicomponent sensor streamer
CN109765620A (zh) * 2019-02-20 2019-05-17 自然资源部第二海洋研究所 近底拖曳式随机接收电缆地震数据采集系统与方法
CN110780350A (zh) * 2019-10-25 2020-02-11 中国海洋大学 一种海底地震检波器动态定位方法、系统及检波器
CN112505750A (zh) * 2020-11-18 2021-03-16 广州海洋地质调查局 一种深拖多道地震拖缆姿态确定方法及处理终端
CN114460649A (zh) * 2022-04-14 2022-05-10 自然资源部第一海洋研究所 一种深海近底拖曳式多道地震接收阵列形态重建方法
CN114779340A (zh) * 2022-06-21 2022-07-22 自然资源部第一海洋研究所 一种深拖式高分辨率多道地震数据精细处理方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
YANLIANG PEI等: "Development of a high-resolution deep-towed multi-channel seismic exploration system: Kuiyang ST2000" *
余文坤等: "基于多项式曲线拟合的海上地震勘探拖缆定位", pages 773 - 779 *
孔繁达;何涛;: "深拖曳多道地震系统阵列几何形态的高精度反演" *
陈桂秀: "用程序解加权最小二乘拟合多项式的系数", pages 48 - 49 *
魏峥嵘等: "深拖式多道高分辨率地震探测系统在南海首次应用", pages 967 - 969 *

Also Published As

Publication number Publication date
CN115826056B (zh) 2023-06-16

Similar Documents

Publication Publication Date Title
US8531912B2 (en) Method and device to acquire marine seismic data
MXPA02012547A (es) Un metodo para procesar datos sismicos.
EP3078991B1 (en) Method for swell effect and mis-tie correction in high-resolution seismic data using multi-beam echo sounder data
GB2517329A (en) Method and device to acquire seismic data
CN114460649B (zh) 一种深海近底拖曳式多道地震接收阵列形态重建方法
CN111487678B (zh) 一种确定高分辨率小多道地震最小偏移距和系统延迟的分析方法
CN110132281B (zh) 一种基于询问应答模式的水下高速目标高精度自主声学导航方法
CN114779340B (zh) 一种深拖式高分辨率多道地震数据精细处理方法
US8717845B2 (en) Quality-based steering methods and systems for 4D geophysical surveys
CN111220146B (zh) 一种基于高斯过程回归学习的水下地形匹配定位方法
WO2019043452A1 (en) ESTIMATING SOURCE-RECEIVER POSITION BY DIRECT ARRIVAL MODELING AND INVERSION
CN112327362A (zh) 速度域的海底多次波预测与追踪衰减方法
CN116879901A (zh) 一种矢量水听器垂直阵列深海宽带声源深度估计方法
CN110221278B (zh) 一种基于多传感器组合的合成孔径声呐运动补偿方法
CN115826056A (zh) 深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法
US20140249757A1 (en) Apparatus and method for determination of far-field signature from variable-depth seismic data
Dosso et al. Array element localization for towed marine seismic arrays
CN113075732B (zh) 一种消除高分辨率小多道地震地层异常起伏的方法
CN116358544A (zh) 基于声学特征匹配定位修正惯性导航误差的方法和系统
CN109490962A (zh) 一种浅地层剖面非零偏移距的消除方法
CN113009417B (zh) 利用声场干涉特性的海底声学阵列阵形估计方法
CN113189646B (zh) 一种去除拖拽式浅剖地层异常起伏的方法
CN112083466B (zh) 一种顾及时间偏差的海底应答器定位方法及系统
CN113703050B (zh) 一种深海地震垂直缆二次定位方法
CN113075734B (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
GR01 Patent grant
GR01 Patent grant