CN112526603B - 一种有线采集和节点站采集的地震数据合成方法 - Google Patents
一种有线采集和节点站采集的地震数据合成方法 Download PDFInfo
- Publication number
- CN112526603B CN112526603B CN202011319005.2A CN202011319005A CN112526603B CN 112526603 B CN112526603 B CN 112526603B CN 202011319005 A CN202011319005 A CN 202011319005A CN 112526603 B CN112526603 B CN 112526603B
- Authority
- CN
- China
- Prior art keywords
- data
- file
- wired
- acquisition
- time
- 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.)
- Active
Links
- 238000001308 synthesis method Methods 0.000 title description 3
- 238000005070 sampling Methods 0.000 claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 24
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 9
- 238000012937 correction Methods 0.000 claims abstract description 9
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 9
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 9
- 238000006243 chemical reaction Methods 0.000 claims abstract description 6
- 230000005284 excitation Effects 0.000 claims description 16
- 238000005422 blasting Methods 0.000 claims description 6
- 230000004927 fusion Effects 0.000 abstract description 7
- 238000004519 manufacturing process Methods 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 12
- 230000003287 optical effect Effects 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/24—Recording seismic data
- G01V1/247—Digital recording of seismic data, e.g. in acquisition units or nodes
-
- 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/288—Event detection in seismic signals, e.g. microseismics
-
- 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/32—Transforming one recording into another or one representation into another
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Geophysics And Detection Of Objects (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种有线采集和节点站采集的地震数据合成方法,其包括如下步骤:第一步,获得有线采集设备和无线节点站中采集时钟的时钟差Δt:第二步,根据时钟差Δt和采样间隔T,采用三次样条插值的方法,对无线节点站的采集数据进行矫正;第三步,通过数模转换校正,将有线采集设备和无线节点站的采集数据合成,形成统一的炮集数据文件。本发明将不同系统采集的数据进行时钟同步,实现了有线采集设备和无线采集设备采集的数据进行无缝融合,通过残差分析,证明数据合成几乎不存在相差,满足生产应用的需求。
Description
技术领域
本发明属于地震数据采集技术领域,特别是涉及一种有线采集和节点站采集的地震数据合成方法。
背景技术
低功耗的无线节点地震数据采集系统是信息技术与地球物理勘探深入融合的产物,是实现低成本油气资源勘探的重要技术手段。目前有线地震数据采集系统主要是从国外进口,且成本高,受制于国外。但由于现有油气勘探企业还保存一定数量的有线采集设备,且在较长的时期生产上还会大量使用。由于有线采集设备需要铺设大量数据传输光缆,有一些复杂区域铺设光缆的难度较大,比如陡岩、城镇、江河湖泊等,导致采集数据缺失比较严重。故在未来一个时期,利用无线节点站的便捷性将其作为有线采集系统的补充,在无法布设有线采集设备的区域用无线节点站进行采集,从而形成有线和无线采集设备混合采集的地震数据采集系统。
目前,勘探设备的有线设备与节点设备多种采集设备联合采集项目日渐增多。由于各种采集设备使用的检波器类型、仪器采集幅值及时钟一致性处理均不一致,导致有线采集和节点站采集的地震数据合成效果较差,没有统一标准。
如图1所示,由于部分区域无法布设有线设备,导致采集的炮集数据存在缺失。利用无线节点站布设的便捷性,可以将无线节点站布设到缺失区域,从而保证了炮集数据的完整性。数据合成的目标就是将无线节点站采集的数据和有线设备采集数据融合,形成统一的炮集数据。
在有线设备和无线设备混合采集系统中,其中一个关键问题是要实现有线采集设备和无线采集设备分别采集的数据进行融合。但无线节点站和有线设备是相互独立的控制系统,导致整体采集系统无法实现时钟同步,进而无法保证地震数据的采集质量。所以,数据融合的难点在于如何实现不同系统采集的数据进行时钟同步。目前,现有的产品和方法研究还未出现可靠的解决方案。
发明内容
本发明的目的就是提供一种有线采集和节点站采集的地震数据合成方法,实现了有线采集设备和无线采集设备采集的数据进行无缝融合,解决了混合采集系统的关键技术问题,能完全解决上述现有技术的不足之处。
本发明的目的通过下述技术方案来实现:一种有线采集和节点站采集的地震数据合成方法,其包括如下步骤:
第一步,获得有线采集设备和无线节点站中采集时钟的时钟差Δt:
第二步,根据时钟差Δt和采样间隔T,采用三次样条插值的方法,对无线节点站的采集数据进行矫正;
第三步,通过数模转换校正,将有线采集设备和无线节点站的采集数据合成,形成统一的炮集数据文件。
作为优选方式之一,三次样条插值的方法为:利用节点站采集数据形成一个时间序列,构造三次样条函数:v=f(t),其中t为时间,v为采样值;计算出任意一个时刻的函数值,获得矫正后的采样数据。
作为优选方式之一,采样数据包含时间点信息和电压振幅值。
作为优选方式之一,通过系统启动数据记录,即放炮系统将该时刻对应的GPS时间形成记录文件,系统自动记录包含GPS时间信息即t0时间的采集日志文件。
作为优选方式之一,第二步中的Δt取值为-0.5*T~0.5*T之间,其中T为采样周期。
作为优选方式之一,将有线采集设备和无线节点站的采集数据合成的方法为:
第一步:准备好相关参数和数据,包括Seg-D炮记录文件,激发日志shotlog文件和无线节点站采集数据文件;其中,Seg-d文件提供有线设备的数据,shotlog文件提供放炮开始时间T0,无线节点采集数据文件提供打上时间戳的无线节点采集数据;
第二步:读取shotlog文件获得激发脉冲的GPS时间,并将此作为有线设备记录数据的零时刻T0,并计算出Δt;
第三步:读取Seg-D文件的头部信息,获得记录长度L和采样率T;
第四步:遍历Seg-D记录文件的每一道道头,如果当前道为有线设备记录数据则读取下一道道头;如果当前道为空道,则表明该道有线系统缺失,转到第五步;
第五步:利用T0、T和L,从节点站数据中读取该时段的记录,数据形成时间序列Q;
第六步:利用Q构建三次样条函数f(t);
第七步:利用T0、T、L和f(t),重新计算出
T0,T0+T,T0+2T,……,T0+L*T时刻的函数值;
第八步:如果当前道为最后一道,则合成结束,否则转到第四步。
作为优选方式之一,第五步中,形成时间序列Q的方法如下:无线节点仪器设备从开机后就开始记录数据,并形成数据文件,根据shotlog记录的激发时间T0,并根据采样率计算T0时间的数据记录在文件中的偏移量Sstart,根据记录长度和采样率计算数据截取时间对应数据文件的偏移量Send;根据Sstart和Send,可以从无线节点仪器数据记录文件中读取对应的数据,形成序列Q。
与现有技术相比,本发明的有益效果在于:本发明有线采集和节点站采集的地震数据合成方法,将不同系统采集的数据进行时钟同步,实现了有线采集设备和无线采集设备采集的数据进行无缝融合,通过残差分析,证明数据合成几乎不存在相差,满足生产应用的需求。
附图说明
图1是有线设备和无线设备混合采集观测系统的示意图。
图2是无线节点站和有线设备采样时钟不统一的示意图。
图3无线节点站采集信号的示意图。
图4根据有线设备时钟差矫正的信号示意图。
图5激发脉冲捕获后的示意图。
图6有线采集设备和无线节点站两者的时钟差Δt的计算表达示意图。
图7节点站采集到的信号矫正时的示意图。
图8有限采集和节点站采集的地址数据融合流程图。
图9为有线设备形成的Seg-D炮集数据文件示意图。
图10为节点站形成的炮集数据文件示意图。
图11为残差示意图。
图12为无线节点仪器数据记录文件。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1:
实现有线设备采集数据和无线节点站采集数据融合,需要以下步骤:
第一步,获得有线采集设备和无线节点站中采集时钟的差Δt:
有线采集设备系统采用的是内部时钟,实现了整个系统的时钟同步,而无线设备是利用GPS授时实现了所有无线节点站的时钟同步。两套系统的时钟不统一,从而形成了相同的地震信号,有线采集设备的内部时钟和无线节点站的GPS时钟参考时间0时刻不同。而数据合成需要将无线节点站采集设备的信号从新矫正到有线设备的内部时钟,因此,其矫正方法过程如下:
如图2所示,节点站采集信号中的中部曲线为地震信号,上方为设定的有线采集设备的采样时钟信号,下方为设定的无线节点站的采样时钟信号,两者存在一定的时钟差。由于有线采集设备和无线节点站采用了不同时钟,因此,设定两者的时钟差为Δt。为此,需要根据Δt,将图2所示的节点站采集信号矫正到图3所示的矫正信号。在同一个采样区域内,图3代表是无线仪器自己的时钟采样点的示意图,图4是校正后采样点跟有线时间重合的示意图。
如图5和图6所示,图5中激发脉冲的捕获示意图,激发信号脉冲和GPS时间点被同时记录。因为无线节点仪采用的是GPS的整毫秒作为0时刻点,而有线设备采用的是内部时钟,因此,为了获得Δt,需要捕获激发脉冲信号,并记录该时刻的GPS时间;由于激发脉冲是有线设备数据记录的参考时钟零点,因此,通过系统启动数据记录,即放炮系统将该时刻对应的GPS时间形成记录文件,系统自动记录包含GPS时间信息即t0时间的采集日志文件。
由于无线节点仪站在开机后就以GPS整毫秒作为零点启动数据记录,而Δt的GPS时间精度是微秒。故Δt取值为-0.5*T~0.5*T之间,其中T为采样周期。
Δt的计算方式:由有线设备采集时的t0时间减去无线节点的整毫秒采样的t1时间差获得,即Δt=|t0-t1|。
比如:有线设备所记录的时间为15h30m20s532ms136μs;无线设备所采集样点的位置在15h30m20s532ms,Δt=136μs。
第二步,根据时钟差Δt和采样间隔T,将无线节点站的采集数据进行矫正;无线节点站是从GPS时间的整毫秒开始记录数据,利用捕获激发脉冲获得信号记录的零时刻T0,其中GPS时间,精度为微秒,可以计算出Δt。
具体的矫正可以采用三次样条插值的方法:即,首先第一步是先把离散信号通过数模转换为连续信号;第二步再通过模数转换,将连续信号转换为时钟点对齐的离散信号;然后进行三次样条进行比对,进行矫正。如图7所示,利用节点站采集数据形成一个时间序列,构造三次样条函数:v=f(t),其中t为时间,v为采样值。从而可以计算出任意一个时刻的函数值,获得矫正后的采样数据;采样数据包含时间点信息和电压振幅值。具体的数学原理如下:
已知区间序列a<x1<x2<...<xn<b,对应的样本值为则三次样条函数s(x):
其中hi=xi+1-xi(i=0,1,...,n-1),Mi=sn(xi)为待定参数,M0,M1,...,Mn满足线性方程组:
该方程是含有n+1个未知数Mi(i=0,1,...,n)的由n-1个方程组成的线性方程组,不能定解。因此,需要补充两个条件:在插值区间两个端点各补充一个条件,通常称为端点条件。最常用的端点条件有三种类型:
s′(x0)=f′0,s′(xn)=f′n
s′(x0)=f′0,s″(xn)=f″n
s(x0)=s(xn)(j=0,1,2)
用Mi表示,这三种条件依次为:
M0=f″0,Mn=f″n
M0=Mn,λnM1+μnMn-1+2Mn=dn
其中:μn=1-λn,/>
第三步,通过数模转换校正后,将有线采集设备和无线节点站的采集数据合成,形成统一的炮集数据文件。合成流程和详细步骤如下:
有线设备系统产生炮集记录文件seg-D和激发时刻的GPS时间T0,记录长度L和采样间隔T由施工参数确定。
根据如图8所示的流程图,本方法的详细步骤如下:
第一步:准备好相关参数和数据,包括Seg-D炮记录文件,激发日志shotlog文件和无线节点站采集数据文件;其中,Seg-d文件提供有线设备的数据,shotlog文件提供放炮开始时间T0,无线节点采集数据文件提供打上时间戳的无线节点采集数据。
第二步:读取shotlog文件获得激发脉冲的GPS时间,并将此作为有线设备记录数据的零时刻T0,并计算出Δt;
第三步:读取Seg-D文件的头部信息,获得记录长度L和采样率T;
第四步:遍历Seg-D记录文件的每一道道头,如果当前道为有线设备记录数据则读取下一道道头;如果当前道为空道,则表明该道有线系统缺失,转到第五步;
第五步:利用T0、T和L,从节点站数据中读取该时段的记录,数据形成时间序列Q;具体过程如下:首先,无线节点仪器设备从开机后就开始记录数据(GPS时间整毫秒开始记录),并形成数据文件,根据shotlog记录的激发时间T0(GPS时间,但不是整毫秒)。并根据采样率计算T0时间的数据记录在文件中的偏移量Sstart,同理,根据记录长度和采样率计算数据截取时间对应数据文件的偏移量Send。根据Sstart和Send,可以从无线节点仪器数据记录文件中读取对应的数据,形成图12所示无线节点仪器数据记录文件中的序列Q。
第六步:利用Q构建三次样条函数f(t);
第七步:利用T0、T、L和f(t),重新计算出T0,T0+T,T0+2T,……,T0+L*T时刻的函数值;
第八步:如果当前道为最后一道,则合成结束,否则转到第四步。
试验通过节点站和有线设备在相同位置进行地震信号采集,并通过本方法进行数据切割,形成的炮集数据文件。结果如图9-11对比所示,图9为有线设备形成的Seg-D炮集数据文件,图10为节点站形成的炮集数据文件,图11残差。通过残差分析,证明数据合成相差小于20微秒,几乎不存在相差,满足生产应用的需求和系统设计目标。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种有线采集和节点站采集的地震数据合成方法,其特征在于,包括如下步骤:
第一步:获得有线采集设备和无线节点站中采集时钟的时钟差Δt,准备好相关参数和数据,包括Seg-D炮记录文件,激发日志shotlog文件和无线节点站采集数据文件;其中,Seg-D文件提供有线设备的数据,shotlog文件提供放炮开始时间T0,无线节点站采集数据文件提供打上时间戳的无线节点采集数据;读取shotlog文件获得激发脉冲的GPS时间,并将此作为有线设备记录数据的零时刻T0,并计算出Δt;
第二步:根据时钟差Δt和采样间隔T,采用三次样条插值的方法,对无线节点站的采集数据进行矫正:
(2.1),读取Seg-D文件的头部信息,获得记录长度L和采样率f;
(2.2),遍历Seg-D记录文件的每一道道头,如果当前道为有线设备记录数据则读取下一道道头;如果当前道为空道,则表明该道有线系统缺失,转到(2.3);
(2.3),利用T0、T和L,从节点站数据中读取该时段的记录,数据形成时间序列Q;形成时间序列Q的方法如下:无线节点仪器设备从开机后就开始记录数据,并形成数据文件,根据shotlog记录的激发时间T0,并根据采样率计算T0时间的数据记录在文件中的偏移量Sstart,根据记录长度和采样率计算数据截取时间对应数据文件的偏移量Send;根据Sstart和Send,可以从无线节点仪器数据记录文件中读取对应的数据,形成序列Q;
(2.4),利用Q构建三次样条函数v=f(t),其中,t为时间,v为采样值;计算出任意一个时刻的函数值,获得矫正后的采样数据;采样数据包含时间点信息和电压振幅值;
(2.5),利用T0、T、L和f(t),重新计算出T0,T0+T,T0+2T,……,T0+L时刻的函数值;
第三步:通过数模转换校正,将有线采集设备和无线节点站的采集数据合成,形成统一的炮集数据文件;如果当前道为最后一道,则合成结束,否则转到(2.2)。
2.根据权利要求1所述的有线采集和节点站采集的地震数据合成方法,其特征在于:通过系统启动数据记录,即放炮系统将该时刻对应的GPS时间形成记录文件,系统自动记录包含GPS时间信息即T0时间的采集日志文件。
3.根据权利要求1所述的有线采集和节点站采集的地震数据合成方法,其特征在于:第二步中的Δt取值为-0.5*T~0.5*T之间,其中T为采样间隔。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011319005.2A CN112526603B (zh) | 2020-11-23 | 2020-11-23 | 一种有线采集和节点站采集的地震数据合成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011319005.2A CN112526603B (zh) | 2020-11-23 | 2020-11-23 | 一种有线采集和节点站采集的地震数据合成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112526603A CN112526603A (zh) | 2021-03-19 |
CN112526603B true CN112526603B (zh) | 2023-12-26 |
Family
ID=74992615
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011319005.2A Active CN112526603B (zh) | 2020-11-23 | 2020-11-23 | 一种有线采集和节点站采集的地震数据合成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112526603B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108051851A (zh) * | 2017-11-27 | 2018-05-18 | 福州华虹智能科技股份有限公司 | 一种地震勘探方法及系统 |
CN109884688A (zh) * | 2019-01-18 | 2019-06-14 | 潜能恒信能源技术股份有限公司 | 一激多收地震数据采集系统 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080170469A1 (en) * | 2007-01-16 | 2008-07-17 | Aram Systems, Ltd. | Stabilizing remote clocks in a network |
US9075545B2 (en) * | 2012-08-01 | 2015-07-07 | Hewlett-Packard Development Company, L.P. | Synchronizing sensor data using timestamps and signal interpolation |
-
2020
- 2020-11-23 CN CN202011319005.2A patent/CN112526603B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108051851A (zh) * | 2017-11-27 | 2018-05-18 | 福州华虹智能科技股份有限公司 | 一种地震勘探方法及系统 |
CN109884688A (zh) * | 2019-01-18 | 2019-06-14 | 潜能恒信能源技术股份有限公司 | 一激多收地震数据采集系统 |
Also Published As
Publication number | Publication date |
---|---|
CN112526603A (zh) | 2021-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104008212B (zh) | 一种地理位置信息相关联的物联网时间序列数据存取方法 | |
CN110887460B (zh) | 一种基于激光扫描仪和里程计的隧道变形监测方法及系统 | |
CN109443188B (zh) | 一种双层多维滑坡监测方法 | |
CN110927510B (zh) | 一种输电线路双端行波故障测距的频域方法 | |
CN111596093B (zh) | 一种基于adcp的海水流速数据处理方法 | |
CN107831516B (zh) | 融合gnss和地面监测网获取大坝实时高精度位移的方法 | |
CN112526603B (zh) | 一种有线采集和节点站采集的地震数据合成方法 | |
CN104933235B (zh) | 一种融合近海多颗卫星海表高度异常数据的方法 | |
CN105589102A (zh) | 一种海上时移地震记录的潮汐校正方法 | |
CN115203625A (zh) | 一种旱涝指数数据缺失值插补方法及其插补装置 | |
CN103412338A (zh) | 伪随机码电法仪 | |
CN109143342B (zh) | 一种地震勘探无线采集数据融合方法 | |
CN103998714A (zh) | 井下可视化方法 | |
CN108375789B (zh) | 联合采集地震数据的同步匹配方法 | |
CN116989973A (zh) | 一种暂冲式高速风洞连续变迎角测压试验数据处理方法 | |
Leadabrand et al. | Preliminary results on the wavelength dependence and aspect sensitivity of radar auroral echoes between 50 and 3000 MHz | |
CN113048869B (zh) | 一种电子换挡器的标定方法、系统、设备及存储介质 | |
CN212031583U (zh) | 直流杂散电流干扰检测装置及系统 | |
CN108051851A (zh) | 一种地震勘探方法及系统 | |
CN201156089Y (zh) | Hrt波地震预报监测系统 | |
CN219452083U (zh) | 一种基于乒乓操作的微柱形电阻率同步采集系统 | |
CN117572432B (zh) | 一种海底基准站声呐阵列及定位方法 | |
CN109684317B (zh) | 一种测月雷达探测数据的数据编辑方法 | |
CN110109185B (zh) | 一种基于选频法的多日变点联测的被动场源二维电场勘探方法 | |
JP2595791B2 (ja) | 雑音データ処理方法とその装置 |
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 |