CN113552616B - 一种检波器数据的补偿方法及系统 - Google Patents

一种检波器数据的补偿方法及系统 Download PDF

Info

Publication number
CN113552616B
CN113552616B CN202010332346.7A CN202010332346A CN113552616B CN 113552616 B CN113552616 B CN 113552616B CN 202010332346 A CN202010332346 A CN 202010332346A CN 113552616 B CN113552616 B CN 113552616B
Authority
CN
China
Prior art keywords
detector
signal
sensitivity
signals
compensation factor
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
Application number
CN202010332346.7A
Other languages
English (en)
Other versions
CN113552616A (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.)
JIANGSU EARTHQUAKE ADMINISTRATION
Original Assignee
JIANGSU EARTHQUAKE ADMINISTRATION
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 JIANGSU EARTHQUAKE ADMINISTRATION filed Critical JIANGSU EARTHQUAKE ADMINISTRATION
Priority to CN202010332346.7A priority Critical patent/CN113552616B/zh
Publication of CN113552616A publication Critical patent/CN113552616A/zh
Application granted granted Critical
Publication of CN113552616B publication Critical patent/CN113552616B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种检波器数据的补偿方法及系统，属于地震勘探领域中的地震数据的信号处理技术领域。所述方法利用布置在同一个检波点处的高灵敏度检波器和常规检波器同时采集信号，然后利用高灵敏度检波器采集到的信号重建常规检波器的信号，从而恢复和补偿部分高频信号。本发明通过高灵敏度检波器信号和常规检波器信号之间的关系，找出其高频补偿因子，对低频保存较好的常规检波器信号补偿部分高频信号，使其信号的频道更宽，高频更丰富，特别是在高分辨率的勘探中，它们可以提供更丰富的反演结果。

Description

一种检波器数据的补偿方法及系统
技术领域
本发明属于地震勘探领域中的地震数据的信号处理技术领域，具体涉及一种检波器数据的补偿方法及系统。
背景技术
目前地震勘探领域中，对地震勘探检波器信号一般都是做低频恢复和补偿。因为反演过程总是从低频到高频，如果低频成份缺失，则很难得到满意的反演结果。实际上，高频信号也非常重要，特别是在高分辨率的勘探中，它们可以提供更丰富的反演结果。
目前，高频补偿恢复的方法有：利用VSP数据求取不同深度地层的反滤波器补偿高频成分；利用地震微测井进行地震信号高频补偿方法，这些方法存在以下不足：在求取不同深度地层的反滤波器时，会受初至波信噪比的影响。如果初至波信噪比低，则地震信号高频补偿效果就会不好。
发明内容
本发明的目的在于解决上述现有技术中存在的难题，提供一种检波器数据的补偿方法及系统，恢复和补偿检波器数据中的部分高频信号，使其频带更宽、更丰富，提高信号的高频质量。
本发明是通过以下技术方案实现的：
一种检波器数据的补偿方法，所述方法利用布置在同一个检波点处的高灵敏度检波器和常规检波器同时采集信号，然后利用高灵敏度检波器采集到的信号重建常规检波器的信号，从而恢复和补偿部分高频信号。
所述方法包括：
(1)采集常规检波器和高灵敏度检波器的信号：在一个检波点处同时设置常规检波器和高灵敏度检波器，并且高灵敏度检波器的埋深和常规检波器的埋深相同，同时采集常规检波器的信号和高灵敏度检波器的信号；
(2)对高灵敏检波器的信号和常规检波器的信号建立差平方的关系式；
(3)求取所述差平方的关系式的极小值，得到关于补偿因子的Wiener-Hopf方程；
(4)求解所述Wiener-Hopf方程得到补偿因子；
(5)利用所述补偿因子对常规检波器的信号进行重建，得到匹配补偿后的信号。
所述步骤(1)中的同时采集常规检波器的信号和高灵敏度检波器的信号的操作包括：
在采集排列滚动时，对位于同一排或者同一列中的所有检波器同步发送命令，并同步接收常规检波器采集到的地震数据和高灵敏度检波器采集到的地震数据，即获得了常规检波器的信号和高灵敏度检波器的信号。
所述步骤(2)的操作包括：
假设常规检波器的信号为xj(t)，高灵敏度检波器的信号为yj(t)，构造补偿因子pj(t)，建立以下差平方的关系式：
E1(t)＝||pj(t)*yj(t)-xj(t)||2 (1.1)
所述步骤(3)的操作包括：
应用最小二乘原理对所述差平方的关系式取极小值：
求解得到关于补偿因子pj(t)的Wiener-Hopf方程为：
Rxxpj＝Rzx(1.3)
其矩阵形式为：
其中，Rxx是托普利兹(Toeplitz)矩阵，其值为xj(t)的自相关函数；pj(t)是补偿因子；Rzx是yj(t)和xj(t)的互相关函数，N表示补偿因子个数。
所述步骤(4)的操作包括：
先将Rxx和Rzx的每个值利用傅里叶算法转成频域得到频域的振幅谱；
然后对所述振幅谱进行统计得到振幅谱的高斯分布；
最后取振幅谱高的一段范围作为分频匹配的范围，获得该分频匹配的范围内的不同振幅值的个数，该个数即为N；
然后利用Wiener-Hopf方程计算得到补偿因子。
所述取振幅谱高的一段范围作为分频匹配的范围的操作包括：
将距离高斯分布的均值0的左、右的概率为0.8的范围作为分频匹配的范围。
所述步骤(5)的操作包括：
利用下式对常规检波器的信号进行重建，得到匹配补偿后的信号
本发明还提供一种检波器数据的补偿系统，所述系统包括：
采集单元：在一个检波点处同时设置常规检波器和高灵敏度检波器，并且高灵敏度检波器的埋深和常规检波器的埋深相同，同时采集常规检波器的信号和高灵敏度检波器的信号；
差平方的关系式建立单元：与所述采集单元连接，对高灵敏检波器的信号和常规检波器的信号建立差平方的关系式；
方程建立单元：与所述差平方的关系式建立单元连接，求取所述差平方的关系式的极小值，得到关于补偿因子的Wiener-Hopf方程；
补偿因子获取单元：与所述方程建立单元连接，求解所述Wiener-Hopf方程得到补偿因子；
匹配补偿单元：与所述补偿因子获取单元连接，利用所述补偿因子对常规检波器的信号进行重建，得到匹配补偿后的信号。
本发明还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机可执行的至少一个程序，所述至少一个程序被所述计算机执行时使所述计算机执行本发明的检波器数据的补偿方法中的步骤。
与现有技术相比，本发明的有益效果是：本发明通过高灵敏检波器信号和常规检波器信号之间的关系，找出其高频补偿因子，对低频保存较好的常规检波器信号补偿部分高频信号。使其信号的频道更宽，高频更丰富，特别是在高分辨率的勘探中，它们可以提供更丰富的反演结果。
附图说明
图1本发明方法的步骤框图；
图2本发明系统的组成结构图。
具体实施方式
下面结合附图对本发明作进一步详细描述：
为了提高信号的高频质量，本发明首先同时采集高灵敏检波器信号和常规检波器信号，再对高灵敏检波器信号和常规检波器信号建立误差关系式；然后求其误差最小值，得到其补偿因子；最后使用补偿因子，重建常规检波器信号。从而恢复和补偿部分高频信号，使其频带更宽、更丰富。
如图1所示，本发明方法的步骤包括：
(1)采集常规检波器和高灵敏度检波器的信号：在一个检波点同时设置常规检波器和高灵敏度检波器，并且高灵敏度检波器的埋深和常规检波器的埋深相同，在采集排列滚动时，对位于同一排或者同一列中的所有检波器(包括常规检波器和高灵敏度检波器)同步发送命令，并同步接收各个检波器的地震数据，最终获得常规检波器采集到的地震数据和高灵敏度检波器采集到的地震数据，即获得了常规检波器的信号和高灵敏度检波器的信号；
(2)对高灵敏检波器信号和常规检波器信号建立差平方的关系式，具体如下：
假设常规检波器的信号为xj(t)，高灵敏度检波器的信号为yj(t)，构造补偿因子pj(t)，使得目标泛函：
E1(t)＝||pj(t)*yj(t)-xj(t)||2 (1.1)
公式(1.1)即为差平方的关系式。
(3)使用最小二乘原理，求上述差平方的关系式的极小值，得到关于补偿因子的Wiener-Hopf方程，具体如下：
应用最小二乘原理对式(1.1)取极小值：
求解得到关于补偿因子pj(t)的Wiener-Hopf方程为：
Rxxpj＝Rzx (1.3)
其矩阵形式为：
其中，Rxx是托普利兹(Toeplitz)矩阵，其值为xj(t)的自相关函数；pj(t)向量是补偿因子；Rzx是yj(t)和xj(t)的互相关函数，N表示补偿因子个数。
(4)使用Levinson-Durbin算法求解方程(1.4)，即可得到补偿因子，具体如下：
方程(1.4)中的N没有确定，本发明通过将方程(1.3)转换到频域来分析常规数据和高灵敏度数据在频域的差别，确定分频匹配的范围，并求得该分频匹配的范围内的补偿因子。
具体的，先将Rxx和Rzx的每个值利用傅里叶算法转成频域即可得到频域的振幅谱，然后对振幅谱进行统计得到振幅谱的高斯分布，最后取振幅谱高的一段范围(即把距离高斯分布的均值0的左、右的概率为0.8的范围(即从概率为-0.8到概率为0.8的范围))作为分频匹配的范围，获得该分频匹配的范围内的不同振幅值的个数，该个数即为N。然后使用上述Wiener-Hopf方程计算出补偿因子。
(5)把补偿因子和常规检波器信号进行卷积，即可得到高频段更宽和更丰富的信号，即对目标数据进行全区记录补偿，具体如下：
高灵敏度补偿后数据如下：
其中，为匹配补偿后数据，xj(t)为匹配前数据，即常规检波器j所获得的信号为，pj(t)为分频匹配滤波器，此处将步骤(4)求得的补偿因子作为分频匹配滤波器。
如图2所示，本发明还提供一种检波器数据的补偿系统，所述系统包括：
采集单元10：在一个检波点处同时设置常规检波器和高灵敏度检波器，并且高灵敏度检波器的埋深和常规检波器的埋深相同，同时采集常规检波器的信号和高灵敏度检波器的信号；
差平方的关系式建立单元20：与所述采集单元10连接，对高灵敏检波器信号和常规检波器信号建立差平方的关系式；
方程建立单元30：与所述差平方的关系式建立单元20连接，求取所述差平方的关系式的极小值，得到关于补偿因子的Wiener-Hopf方程；
补偿因子获取单元40：与所述方程建立单元30连接，求解所述Wiener-Hopf方程得到补偿因子；
匹配补偿单元50：与所述补偿因子获取单元40连接，利用所述补偿因子对常规检波器的信号进行重建，得到匹配补偿后的信号。
利用本发明方法使得信号高频段更宽和更丰富，特别是在高分辨率的勘探中，它们可以提供更丰富的反演结果。
上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本发明上述具体实施方式所描述的方法，因此前面描述的方式只是优选的，而并不具有限制性的意义。

Claims (8)

1.一种检波器数据的补偿方法，其特征在于：所述方法利用布置在同一个检波点处的高灵敏度检波器和常规检波器同时采集信号，然后利用高灵敏度检波器采集到的信号重建常规检波器的信号，从而恢复和补偿部分高频信号；
所述方法包括：
(1)采集常规检波器和高灵敏度检波器的信号：在一个检波点处同时设置常规检波器和高灵敏度检波器，并且高灵敏度检波器的埋深和常规检波器的埋深相同，同时采集常规检波器的信号和高灵敏度检波器的信号；
(2)对高灵敏检波器的信号和常规检波器的信号建立差平方的关系式；
(3)求取所述差平方的关系式的极小值，得到关于补偿因子的Wiener-Hopf方程；
(4)求解所述Wiener-Hopf方程得到补偿因子；
(5)利用所述补偿因子对常规检波器的信号进行重建，得到匹配补偿后的信号；
所述步骤(4)的操作包括：
先将Rxx和Rzx的每个值利用傅里叶算法转成频域得到频域的振幅谱；
然后对所述振幅谱进行统计得到振幅谱的高斯分布；
最后取振幅谱高的一段范围作为分频匹配的范围，获得该分频匹配的范围内的不同振幅值的个数，该个数即为N，N表示补偿因子个数；
然后利用Wiener-Hopf方程计算得到补偿因子；
其中，Rxx是托普利兹矩阵，其值为xj(t)的自相关函数；Rzx是yj(t)和xj(t)的互相关函数；xj(t)为常规检波器的信号，yj(t)为高灵敏度检波器的信号。
2.根据权利要求1所述的检波器数据的补偿方法，其特征在于：所述步骤(1)中的同时采集常规检波器的信号和高灵敏度检波器的信号的操作包括：
在采集排列滚动时，对位于同一排或者同一列中的所有检波器同步发送命令，并同步接收常规检波器采集到的地震数据和高灵敏度检波器采集到的地震数据，即获得了常规检波器的信号和高灵敏度检波器的信号。
3.根据权利要求2所述的检波器数据的补偿方法，其特征在于：所述步骤(2)的操作包括：
假设常规检波器的信号为xj(t)，高灵敏度检波器的信号为yj(t)，构造补偿因子pj(t)，建立以下差平方的关系式：
E1(t)＝||pj(t)*yj(t)-xj(t)||2(1.1)。
4.根据权利要求3所述的检波器数据的补偿方法，其特征在于：所述步骤(3)的操作包括：
应用最小二乘原理对所述差平方的关系式取极小值：
求解得到关于补偿因子pj(t)的Wiener-Hopf方程为：
Rxxpj＝Rzx(1.3)
其矩阵形式为：
其中，Rxx是托普利兹矩阵，其值为xj(t)的自相关函数；pj(t)是补偿因子；Rzx是yj(t)和xj(t)的互相关函数，N表示补偿因子个数。
5.根据权利要求4所述的检波器数据的补偿方法，其特征在于：所述取振幅谱高的一段范围作为分频匹配的范围的操作包括：
将距离高斯分布的均值0的左、右的概率为0.8的范围作为分频匹配的范围。
6.根据权利要求5所述的检波器数据的补偿方法，其特征在于：所述步骤(5)的操作包括：
利用下式对常规检波器的信号进行重建，得到匹配补偿后的信号
7.一种检波器数据的补偿系统，其特征在于：所述系统包括：
采集单元：在一个检波点处同时设置常规检波器和高灵敏度检波器，并且高灵敏度检波器的埋深和常规检波器的埋深相同，同时采集常规检波器的信号和高灵敏度检波器的信号；
差平方的关系式建立单元：与所述采集单元连接，对高灵敏检波器的信号和常规检波器的信号建立差平方的关系式；
方程建立单元：与所述差平方的关系式建立单元连接，求取所述差平方的关系式的极小值，得到关于补偿因子的Wiener-Hopf方程；
补偿因子获取单元：与所述方程建立单元连接，求解所述Wiener-Hopf方程得到补偿因子；
匹配补偿单元：与所述补偿因子获取单元连接，利用所述补偿因子对常规检波器的信号进行重建，得到匹配补偿后的信号；
补偿因子获取单元进行以下操作：
先将Rxx和Rzx的每个值利用傅里叶算法转成频域得到频域的振幅谱；
然后对所述振幅谱进行统计得到振幅谱的高斯分布；
最后取振幅谱高的一段范围作为分频匹配的范围，获得该分频匹配的范围内的不同振幅值的个数，该个数即为N，N表示补偿因子个数；
然后利用Wiener-Hopf方程计算得到补偿因子；
其中，Rxx是托普利兹矩阵，其值为xj(t)的自相关函数；Rzx是yj(t)和xj(t)的互相关函数；xj(t)为常规检波器的信号，yj(t)为高灵敏度检波器的信号。
8.一种计算机可读存储介质，其特征在于：所述计算机可读存储介质存储有计算机可执行的至少一个程序，所述至少一个程序被所述计算机执行时使所述计算机执行权利要求1-6任一项所述的检波器数据的补偿方法中的步骤。
CN202010332346.7A 2020-04-24 2020-04-24 一种检波器数据的补偿方法及系统 Active CN113552616B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010332346.7A CN113552616B (zh) 2020-04-24 2020-04-24 一种检波器数据的补偿方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010332346.7A CN113552616B (zh) 2020-04-24 2020-04-24 一种检波器数据的补偿方法及系统

Publications (2)

Publication Number Publication Date
CN113552616A CN113552616A (zh) 2021-10-26
CN113552616B true CN113552616B (zh) 2024-03-15

Family

ID=78129633

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010332346.7A Active CN113552616B (zh) 2020-04-24 2020-04-24 一种检波器数据的补偿方法及系统

Country Status (1)

Country Link
CN (1) CN113552616B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5774415A (en) * 1996-12-06 1998-06-30 Geo Space Corporation Geophone normalization process
CN1673775A (zh) * 2005-03-11 2005-09-28 浙江大学 利用地震微测井进行地震信号高频补偿方法
CN201732173U (zh) * 2010-07-26 2011-02-02 中国石油天然气股份有限公司 一种非零炮检距地震信号能量校准装置及系统
CN102200587A (zh) * 2010-11-03 2011-09-28 北京三迭系石油勘探技术开发有限责任公司 一种地震信号接收及其处理方法
CN102323616A (zh) * 2011-06-08 2012-01-18 浙江大学 提高灰岩出露区地震数据分辨率的分频匹配方法
WO2016141598A1 (zh) * 2015-03-12 2016-09-15 中国矿业大学(北京) 共姿态道集地表一致性振幅补偿方法及装置
CN106814394A (zh) * 2015-11-27 2017-06-09 中国石油化工股份有限公司 台站式检波器与常规检波器联合低频恢复方法
CN107144879A (zh) * 2017-05-12 2017-09-08 重庆大学 一种基于自适应滤波与小波变换结合的地震波降噪方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8077541B2 (en) * 2007-10-19 2011-12-13 Westerngeco L.L.C. Testing a sensor to produce a filter for noise attenuation

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5774415A (en) * 1996-12-06 1998-06-30 Geo Space Corporation Geophone normalization process
CN1673775A (zh) * 2005-03-11 2005-09-28 浙江大学 利用地震微测井进行地震信号高频补偿方法
CN201732173U (zh) * 2010-07-26 2011-02-02 中国石油天然气股份有限公司 一种非零炮检距地震信号能量校准装置及系统
CN102200587A (zh) * 2010-11-03 2011-09-28 北京三迭系石油勘探技术开发有限责任公司 一种地震信号接收及其处理方法
CN102323616A (zh) * 2011-06-08 2012-01-18 浙江大学 提高灰岩出露区地震数据分辨率的分频匹配方法
WO2016141598A1 (zh) * 2015-03-12 2016-09-15 中国矿业大学(北京) 共姿态道集地表一致性振幅补偿方法及装置
CN106814394A (zh) * 2015-11-27 2017-06-09 中国石油化工股份有限公司 台站式检波器与常规检波器联合低频恢复方法
CN107144879A (zh) * 2017-05-12 2017-09-08 重庆大学 一种基于自适应滤波与小波变换结合的地震波降噪方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
A regularized Wiener–Hopf filter for inverting models with magnetic susceptibility;Shib Sankar Ganguli等;《Geophysical Prospecting》;第64卷(第2期);第456-468页 *
Seismic absorption compensation: A least squares inverse scheme;Changjun Zhang, Tadeusz J. Ulrych;《GEOPHYSICS》;第72卷(第6期);第109-114页 *
利用微测井资料补偿地震数据的高频成分;田钢等;《石油地球物理勘探》;第40卷(第05期);第546-549页 *
利用数字耦合滤波器消除沙漠层对地震信号的影响;石战结等;《中国地球物理.2003——中国地球物理学会第十九届年会论文集》;第99页 *
地震信号匹配处理方法研究;龙云;《CNKI优秀硕士学位论文全文库》;第1-71页 *
基于Gabor变换的地震资料高分辨率处理方法研究;徐倩茹;孙成禹;乔志浩;唐杰;《断块油气田》;第23卷(第04期);第56-60页 *
沙漠地区地震检波器耦合的高频信号匹配滤波技术;石战结等;《石油物探》;第44卷(第03期);第261-263页 *

Also Published As

Publication number Publication date
CN113552616A (zh) 2021-10-26

Similar Documents

Publication Publication Date Title
CN105116443B (zh) 一种低频信号的能量补偿方法及装置
CN108845357B (zh) 一种基于同步挤压小波变换估计地层等效品质因子的方法
CN106019376B (zh) 一种频率驱动空变q值模型构建的地震波补偿方法
CN113805234B (zh) 在被动源地震数据中增强面波的处理方法
US11372122B2 (en) High-resolution processing method for seismic data based on inverse multi-resolution singular value decomposition
CN104635264B (zh) 叠前地震数据的处理方法及设备
CN113552616B (zh) 一种检波器数据的补偿方法及系统
CN112558159B (zh) 一种声波测井初至拾取方法
CN109765608A (zh) 一种基于联合字典的煤层巷道锚杆振动噪声压制方法
CN104778342B (zh) 一种基于小波奇异熵的心音特征提取方法
CN109143345B (zh) 基于模拟退火的品质因子q非线性反演方法及系统
Mukherjee et al. S-transform based processing of noisy surface wave record for recovering high-resolution spectrum
CN110749923A (zh) 一种基于范数方程提高分辨率的反褶积方法
CN111856559B (zh) 基于稀疏贝叶斯学习理论的多道地震谱反演方法及系统
CN109991657B (zh) 基于逆二分递推奇异值分解的地震资料高分辨率处理方法
CN114371505A (zh) 一种基于地震分频技术的多子波反演方法及系统
CN114428348A (zh) 一种提高地震资料分辨率的方法及系统
CN116413805A (zh) 一种子波一致性时变校正方法、装置、设备和相关系统
Feng et al. Transform learning in the synchrosqueezing frequency domain—A novel denoising strategy for optical fiber seismic records
CN113219532B (zh) 基于低频数据补偿的全波形反演速度建模方法
CN112526611A (zh) 表层地震波品质因子的提取方法及装置
CN117111152A (zh) 一种基于分裂Bregman迭代与宽带低通滤波器的地震记录低频补偿方法
CN113589384B (zh) 基于信号随偏移距变化特征的叠前道集保幅去噪方法
CN113138419B (zh) 提取下行子波和衰减参数的方法、装置
CN111830565B (zh) 一种基于kl-dsw的tsp多波场分离及噪声压制方法

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Kang Qingqing

Inventor after: Gu Qinping

Inventor after: Zhao Qiguang

Inventor after: Meng Ke

Inventor after: Yang Hao

Inventor before: Gu Qinping

Inventor before: Kang Qingqing

Inventor before: Zhao Qiguang

Inventor before: Meng Ke

Inventor before: Yang Hao

GR01 Patent grant
GR01 Patent grant