CN102495431A - 瞬变电磁数据静校正的数据处理方法 - Google Patents

瞬变电磁数据静校正的数据处理方法 Download PDF

Info

Publication number
CN102495431A
CN102495431A CN2011103624041A CN201110362404A CN102495431A CN 102495431 A CN102495431 A CN 102495431A CN 2011103624041 A CN2011103624041 A CN 2011103624041A CN 201110362404 A CN201110362404 A CN 201110362404A CN 102495431 A CN102495431 A CN 102495431A
Authority
CN
China
Prior art keywords
depth
degree
measuring point
static correction
landform
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
CN2011103624041A
Other languages
English (en)
Other versions
CN102495431B (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.)
Xian Research Institute Co Ltd of CCTEG
Original Assignee
Xian Research Institute Co Ltd of CCTEG
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 Xian Research Institute Co Ltd of CCTEG filed Critical Xian Research Institute Co Ltd of CCTEG
Priority to CN201110362404.1A priority Critical patent/CN102495431B/zh
Publication of CN102495431A publication Critical patent/CN102495431A/zh
Application granted granted Critical
Publication of CN102495431B publication Critical patent/CN102495431B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种对瞬变电磁数据进行静校正的数据处理方法。本发明首先确保数据消除了野外电磁干扰及边框效应的影响，根据晚期视电阻率公式计算晚期视电阻率，再根据视深度经验公式计算视深度，根据测点高程数据选择基准点和校正方式，通过各测点视深度与基准点视深度比值确定静校正因子
Figure 2011103624041100004DEST_PATH_IMAGE002
，利用静校正因子
Figure 745388DEST_PATH_IMAGE002
与
Figure 2011103624041100004DEST_PATH_IMAGE004
、的立方及平方关系校正这两个值，最后将得到的校正后
Figure DEST_PATH_IMAGE008
、带入晚期视电阻率公式计算新的
Figure DEST_PATH_IMAGE010
，此时将代入视深度经验公式就可以得到校正后的视深度，将该深度加上基准点高程，即可完成静校正过程。本发明可以快速有效消除定性解释中由于地形起伏产生的虚假异常，大大提高瞬变电磁法资料处理速度及解释精度，从而推动该方法技术的大幅度发展及应用。

Description

瞬变电磁数据静校正的数据处理方法
技术领域   
本发明属于地球物理勘探技术领域，涉及一种对瞬变电磁数据进行静校正的数据处理方法。 
背景技术
静校正的概念较多应用于地球物理勘探领域中的地震勘探技术，随着技术发展，部分电磁法勘探中也逐渐引入了这一概念。例如在频率电磁测深法中，当地表或近地表存在局部电性不均匀二、三维地质体时，不均匀体的表面会形成一种电荷积累效应，而使电场发生畸变，且畸变值与发射频率无关，这就是静态效应，需要进行静校正处理，否则会严重影响定性处理解释结果，带来虚假异常。
瞬变电磁法简称TEM，是近年来广泛应用于资源勘查、工程物探、水文勘探等领域的一种电磁类勘探方法。由于其采集的是二次感应场的信号，而地层介质的二次感应场的变化与几何测深不同，所以之前普遍认为它适应地形复杂变化的能力较强，不需要进行静校正处理，但随着该方法的大量使用，发现在地形条件复杂地区，当实际的观测面起伏不平时，在起伏地表上观测的响应与水平面上观测的响应还是存在较大差异，TEM数据同样会受到地形造成的影响。
事实上，瞬变电磁法虽然属于时间域电磁法，但它与频率域电磁法仍有着千丝万缕的联系，在许多基本公式的推导中，都采用了先在频率域中讨论，最后通过傅氏变换转换为时间域的手段，因此TEM数据的处理解释不可避免的会像频率域电磁法一样受到地形造成的静态效应影响，需要研究有效的静校正处理技术。 
发明内容
本发明所解决的技术问题是提供一种瞬变电磁数据进行静校正的数据处理方法，解决现有瞬变电磁数据处理技术中存在的地形造成的静态效应问题。     为解决上述的技术问题，本发明采取的技术方案：
一种瞬变电磁数据进行静校正的数据处理方法，其特征在于：通过以下步骤实现： 
（1）、首先确保数据消除了野外电磁干扰及边框效应的影响；
（2）、根据晚期视电阻率公式
Figure 2011103624041100002DEST_PATH_IMAGE001
计算晚期视电阻率，式中
Figure 728314DEST_PATH_IMAGE002
亨/米，S T 为发送回线面积，S R 为接收线圈面积，t为测道时间，V(t)／I为归一化感应电动势；
（3）、根据视深度经验公式
Figure 2011103624041100002DEST_PATH_IMAGE003
计算视深度，式中，
Figure 160301DEST_PATH_IMAGE004
为视深度，
Figure 2011103624041100002DEST_PATH_IMAGE005
为经验参数，
Figure 861409DEST_PATH_IMAGE006
为视电阻率，
Figure 318936DEST_PATH_IMAGE008
为测道时间；
（4）、根据测点高程数据选择基准点和校正方式：基准点意味着校正完成后各测点数据均相当于在该点对应高程处进行了测量，校正方式可根据地形影响程度的不同选择按测区校正、按测线方向校正和按测点方向校正；
（5）、通过各测点视深度与基准点视深度比值确定静校正因子；
（6）、利用静校正因子
Figure 38368DEST_PATH_IMAGE009
与、的立方及平方关系校正这两个值，得到
Figure 783656DEST_PATH_IMAGE012
、
Figure 2011103624041100002DEST_PATH_IMAGE013
，
Figure 576031DEST_PATH_IMAGE014
为归一化感应电动势，
Figure 2011103624041100002DEST_PATH_IMAGE015
为测道时间，为地形较低测点经校正后得到的测道时间， 
Figure 848061DEST_PATH_IMAGE012
为地形较低测点经校正后得到的归一化电动势；
（7）、将得到的校正后
Figure 133549DEST_PATH_IMAGE012
、
Figure 650986DEST_PATH_IMAGE016
带入晚期视电阻率公式
计算新的
Figure 236689DEST_PATH_IMAGE018
，为地形较低测点经校正后得到的晚期视电阻率；
（8）、此时将
Figure 479582DEST_PATH_IMAGE018
代入视深度经验公式
Figure 517945DEST_PATH_IMAGE020
就可以得到校正后的视深度，
Figure 2011103624041100002DEST_PATH_IMAGE021
为地形较低测点经校正后得到的视深度，将该深度加上基准点高程，即可完成静校正过程。
与现有技术相比，本发明方法简便，易于操作，不论对于理论数据还是实测数据，均有良好的处理效果，可以快速有效消除定性解释中由于地形起伏产生的虚假异常，大大提高瞬变电磁法资料处理速度及解释精度，从而推动该方法技术的大幅度发展及应用。
附图说明
图1为本发明的总流程图；
图2为本发明需校正单点与基准点关系示意图；
图3为瞬变电磁数据静校正技术实施例的理论模型静校正前效果图；
图4为瞬变电磁数据静校正技术实施例的理论模型静校正后效果图；
图5为瞬变电磁数据静校正技术实施例的实测数据静校正前效果图；
图6为瞬变电磁数据静校正技术实施例的实测数据静校正前效果图。  
具体实施方式
下面结合具体实施方式和说明书附图对本发明进行详细说明。
参见图1，本发明通过以下步骤实现：
（1）、首先确保数据消除了野外电磁干扰及边框效应的影响；
（2）、根据晚期视电阻率公式
Figure 104653DEST_PATH_IMAGE001
计算晚期视电阻率，式中
Figure 59709DEST_PATH_IMAGE002
亨/米，S T 为发送回线面积，S R 为接收线圈面积，t为测道时间，V(t)／I为归一化感应电动势；
（3）、根据视深度经验公式
Figure 207924DEST_PATH_IMAGE003
计算视深度，式中，
Figure 733584DEST_PATH_IMAGE004
为视深度，
Figure 609136DEST_PATH_IMAGE005
为经验参数，
Figure 153119DEST_PATH_IMAGE006
为视电阻率，
Figure 721503DEST_PATH_IMAGE008
为测道时间；
（4）、根据测点高程数据选择基准点和校正方式：基准点意味着校正完成后各测点数据均相当于在该点对应高程处进行了测量，校正方式可根据地形影响程度的不同选择按测区校正、按测线方向校正和按测点方向校正；
（5）、通过各测点视深度与基准点视深度比值确定静校正因子；
（6）、参见图2，以地形较高测点为基准点，利用静校正因子
Figure 633276DEST_PATH_IMAGE009
与
Figure 516918DEST_PATH_IMAGE010
、
Figure 990625DEST_PATH_IMAGE011
的立方及平方关系校正这两个值，得到
Figure 474564DEST_PATH_IMAGE012
、
Figure 426340DEST_PATH_IMAGE013
，
Figure 430068DEST_PATH_IMAGE014
为归一化感应电动势，
Figure 74676DEST_PATH_IMAGE015
为测道时间，
Figure 281797DEST_PATH_IMAGE013
为地形较低测点经校正后得到的测道时间， 
Figure 302843DEST_PATH_IMAGE012
为地形较低测点经校正后得到的归一化电动势；
（7）、将得到的校正后
Figure 895498DEST_PATH_IMAGE012
、
Figure 711008DEST_PATH_IMAGE016
带入晚期视电阻率公式
Figure 739898DEST_PATH_IMAGE022
计算新的，为地形较低测点经校正后得到的晚期视电阻率；
（8）、此时将代入视深度经验公式
Figure 179921DEST_PATH_IMAGE020
就可以得到校正后的视深度，
Figure 2011103624041100002DEST_PATH_IMAGE023
为地形较低测点经校正后得到的视深度，将该深度加上基准点高程，即可完成静校正过程。
下面分别以理论模型和实测数据为例说明本发明效果：
给出一个瞬变电磁中心回线装置下的两层带地形的理论地电模型如图2所示，，，
Figure 687180DEST_PATH_IMAGE026
。从图3中可以看出，理论模型静校正前，在地下110m处，地形变化的位置视电阻率等值线出现一明显“台阶”，形态与地形相似，静校正后得到图4，可看到“台阶”基本消失，视电阻率分界线与模型参数吻合很好，说明了该静校正方法的有效性（图3、图4中黑色虚线为地层）。
图5为某煤矿测区内一条TEM测线静校正前的视电阻率剖面图，剖面上能明显看出，等值线形态变化与地形变化非常相似，与煤层形态很难看出关联，造成同一层位上视电阻率出现“逢沟必低”现象，即逢地形出现沟的位置必定有一个低阻异常产生；经静校正后得到图6，视电阻率等值线与地形的相关程度大大减弱，整体形态与煤层形态几乎完全吻合，低阻圈也基本消失，此时进行分析解释出的异常才是较为可信的（图5、图6中黑色点状线为地形，黑色实线为煤层）。

Claims (1)

1.一种瞬变电磁数据进行静校正的数据处理方法，其特征在于：通过以下步骤实现： 
（1）、首先确保数据消除了野外电磁干扰及边框效应的影响；
（2）、根据晚期视电阻率公式
计算晚期视电阻率，式中
Figure 225704DEST_PATH_IMAGE002
亨/米，S T 为发送回线面积，S R 为接收线圈面积，t为测道时间，V(t)／I为归一化感应电动势；
（3）、根据视深度经验公式
Figure 2011103624041100001DEST_PATH_IMAGE003
计算视深度，式中，为视深度，
Figure 2011103624041100001DEST_PATH_IMAGE005
为经验参数，为视电阻率，
Figure 504742DEST_PATH_IMAGE008
为测道时间；
（4）、根据测点高程数据选择基准点和校正方式：基准点意味着校正完成后各测点数据均相当于在该点对应高程处进行了测量，校正方式可根据地形影响程度的不同选择按测区校正、按测线方向校正和按测点方向校正；
（5）、通过各测点视深度与基准点视深度比值确定静校正因子
Figure 2011103624041100001DEST_PATH_IMAGE009
；
（6）、利用静校正因子
Figure 381431DEST_PATH_IMAGE009
与、
Figure 549555DEST_PATH_IMAGE008
的立方及平方关系校正这两个值，得到
Figure 545193DEST_PATH_IMAGE012
、
Figure 2011103624041100001DEST_PATH_IMAGE013
，为归一化感应电动势，为测道时间，
Figure 17555DEST_PATH_IMAGE013
为地形较低测点经校正后得到的测道时间， 
Figure 303043DEST_PATH_IMAGE016
为地形较低测点经校正后得到的归一化电动势；
（7）、将得到的校正后
Figure 649842DEST_PATH_IMAGE016
、
Figure 2011103624041100001DEST_PATH_IMAGE017
带入晚期视电阻率公式
计算新的
Figure 2011103624041100001DEST_PATH_IMAGE019
，
Figure 508131DEST_PATH_IMAGE020
为地形较低测点经校正后得到的晚期视电阻率；
（8）、此时将
Figure 280915DEST_PATH_IMAGE019
代入视深度经验公式
Figure 2011103624041100001DEST_PATH_IMAGE021
就可以得到校正后的视深度，
Figure 634667DEST_PATH_IMAGE022
为地形较低测点经校正后得到的视深度，将该深度加上基准点高程，即可完成静校正过程。
CN201110362404.1A 2011-11-16 2011-11-16 瞬变电磁数据静校正的数据处理方法 Active CN102495431B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110362404.1A CN102495431B (zh) 2011-11-16 2011-11-16 瞬变电磁数据静校正的数据处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110362404.1A CN102495431B (zh) 2011-11-16 2011-11-16 瞬变电磁数据静校正的数据处理方法

Publications (2)

Publication Number Publication Date
CN102495431A true CN102495431A (zh) 2012-06-13
CN102495431B CN102495431B (zh) 2014-05-21

Family

ID=46187271

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110362404.1A Active CN102495431B (zh) 2011-11-16 2011-11-16 瞬变电磁数据静校正的数据处理方法

Country Status (1)

Country Link
CN (1) CN102495431B (zh)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103064124A (zh) * 2013-01-31 2013-04-24 江苏大学 一种校正电磁勘探中地形影响的比值方法
CN104422970A (zh) * 2013-09-04 2015-03-18 中国石油天然气集团公司 一种电磁数据的保幅静校正方法
CN104698502A (zh) * 2015-03-22 2015-06-10 山西煤炭进出口集团科技发展有限公司 由井出发的瞬变电磁勘探地形校正处理方法
CN104749647A (zh) * 2015-04-09 2015-07-01 中国科学院地质与地球物理研究所 一种地质分界面构造的探测方法、装置和系统
CN104903750A (zh) * 2014-07-23 2015-09-09 王雅苹 一种电磁数据的保幅静校正方法
CN105572748A (zh) * 2015-12-08 2016-05-11 中国科学院地质与地球物理研究所 一种瞬变电磁探测时间窗口的确定方法和装置
CN105676300A (zh) * 2016-03-11 2016-06-15 陕西省煤田物探测绘有限公司 一种消除瞬变电磁观测数据地形影响的方法和装置
CN106646666A (zh) * 2017-01-16 2017-05-10 中南大学 基于平面波电磁测深的静态效应校正方法
CN107024722A (zh) * 2017-04-19 2017-08-08 吉林大学 一种基于异常环的低温超导磁源瞬变电磁地形校正方法
CN108535783A (zh) * 2018-04-16 2018-09-14 中煤科工集团西安研究院有限公司 一种对电阻率断面进行层位校正的方法及装置
CN109541695A (zh) * 2018-11-16 2019-03-29 吉林大学 人工场源频率域电场梯度远区视电阻率快速成像方法
US20190391293A1 (en) * 2018-06-20 2019-12-26 Institute of Geology and Geophysics, Chinese Acade my of Sciences Correction method for transient electromagnetic method-based prediction of water-bearing structure ahead of tunnel wall
CN110717138A (zh) * 2019-09-12 2020-01-21 吉林大学 一种拖曳式瞬变电磁数据预处理方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6903553B2 (en) * 2002-09-06 2005-06-07 Baker Hughes Incorporated Method and apparatus for a quadrupole transmitter for directionally sensitive induction tool
CN101706586A (zh) * 2009-12-03 2010-05-12 成都理工大学 一种时间域视电阻率的计算方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6903553B2 (en) * 2002-09-06 2005-06-07 Baker Hughes Incorporated Method and apparatus for a quadrupole transmitter for directionally sensitive induction tool
CN101706586A (zh) * 2009-12-03 2010-05-12 成都理工大学 一种时间域视电阻率的计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《工程地球物理学报》 20050831 曹国栋等 TEM小线框技术在近地表结构中的应用--以徐州市淮塔地基结构研究为例 270-275 1 第2卷, 第4期 *
曹国栋等: "TEM小线框技术在近地表结构中的应用——以徐州市淮塔地基结构研究为例", 《工程地球物理学报》 *
范涛: "煤田电法勘探中的TEM地形校正方法", 《物探与化探》 *

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103064124A (zh) * 2013-01-31 2013-04-24 江苏大学 一种校正电磁勘探中地形影响的比值方法
CN103064124B (zh) * 2013-01-31 2015-09-16 江苏大学 一种校正电磁勘探中地形影响的比值方法
CN104422970A (zh) * 2013-09-04 2015-03-18 中国石油天然气集团公司 一种电磁数据的保幅静校正方法
CN104422970B (zh) * 2013-09-04 2017-03-15 中国石油天然气集团公司 一种电磁数据的保幅静校正方法
CN104903750A (zh) * 2014-07-23 2015-09-09 王雅苹 一种电磁数据的保幅静校正方法
WO2016011627A1 (zh) * 2014-07-23 2016-01-28 王雅苹 一种电磁数据的保幅静校正方法
CN104698502A (zh) * 2015-03-22 2015-06-10 山西煤炭进出口集团科技发展有限公司 由井出发的瞬变电磁勘探地形校正处理方法
CN104749647A (zh) * 2015-04-09 2015-07-01 中国科学院地质与地球物理研究所 一种地质分界面构造的探测方法、装置和系统
CN105572748A (zh) * 2015-12-08 2016-05-11 中国科学院地质与地球物理研究所 一种瞬变电磁探测时间窗口的确定方法和装置
CN105676300A (zh) * 2016-03-11 2016-06-15 陕西省煤田物探测绘有限公司 一种消除瞬变电磁观测数据地形影响的方法和装置
CN106646666A (zh) * 2017-01-16 2017-05-10 中南大学 基于平面波电磁测深的静态效应校正方法
CN106646666B (zh) * 2017-01-16 2018-08-14 中南大学 基于平面波电磁测深的静态效应校正方法
CN107024722A (zh) * 2017-04-19 2017-08-08 吉林大学 一种基于异常环的低温超导磁源瞬变电磁地形校正方法
CN107024722B (zh) * 2017-04-19 2018-09-28 吉林大学 一种基于异常环的低温超导磁源瞬变电磁地形校正方法
CN108535783A (zh) * 2018-04-16 2018-09-14 中煤科工集团西安研究院有限公司 一种对电阻率断面进行层位校正的方法及装置
CN108535783B (zh) * 2018-04-16 2019-07-05 中煤科工集团西安研究院有限公司 一种对电阻率断面进行层位校正的方法及装置
US20190391293A1 (en) * 2018-06-20 2019-12-26 Institute of Geology and Geophysics, Chinese Acade my of Sciences Correction method for transient electromagnetic method-based prediction of water-bearing structure ahead of tunnel wall
US10823872B2 (en) * 2018-06-20 2020-11-03 Institute Of Geology And Geophysics, Chinese Academy Of Sciences Correction method for transient electromagnetic method-based prediction of water-bearing structure ahead of tunnel wall
CN109541695A (zh) * 2018-11-16 2019-03-29 吉林大学 人工场源频率域电场梯度远区视电阻率快速成像方法
CN110717138A (zh) * 2019-09-12 2020-01-21 吉林大学 一种拖曳式瞬变电磁数据预处理方法
CN110717138B (zh) * 2019-09-12 2023-06-30 吉林大学 一种拖曳式瞬变电磁数据预处理方法

Also Published As

Publication number Publication date
CN102495431B (zh) 2014-05-21

Similar Documents

Publication Publication Date Title
CN102495431B (zh) 瞬变电磁数据静校正的数据处理方法
Stopin et al. Multiparameter waveform inversion of a large wide-azimuth low-frequency land data set in Oman
US8564296B2 (en) Systems and methods for remote electromagnetic exploration for mineral and energy resources
Pedersen et al. Groundwater exploration using combined controlled-source and radiomagnetotelluric techniques
Niederleithinger et al. Evaluation of geophysical techniques for dike inspection
Bastani et al. CSRMT measurements in the frequency range of 1–250 kHz to map a normal fault in the Volvi basin, Greece
Xia et al. Resolution of high-frequency Rayleigh-wave data
Jiu‐Long et al. Synthetic aperture imaging in advanced detection of roadway using mine transient electromagnetic method
US9395463B2 (en) EM processing using field ratios
Liu et al. Mapping water-abundant zones using transient electromagnetic and seismic methods when tunneling through fractured granite in the Qinling Mountains, China
Bastani et al. Delineating hydrothermal stockwork copper deposits using controlled-source and radio-magnetotelluric methods: A case study from northeast Iran
Li et al. Application of grounded electrical source airborne transient electromagnetic (GREATEM) system in goaf water detection
Dennis et al. Transient electromagnetic surveys for the measurement of near-surface electrical anisotropy
CN109917466A (zh) 一种基于电磁场垂直分量的电阻率测量方法
An et al. Application of the CSAMT method for exploring deep coal mines in Fujian Province, Southeastern China
Yu et al. Detection of water-enriched goaf via the semi-airborne transient electromagnetic method
Liu et al. Metallic mineral exploration by using ambient noise tomography in Ashele copper mine, Xinjiang, China
WO2021096877A1 (en) Removing the effect of near-surface inhomogeneities in surface-to-borehole measurements
Davydycheva et al. Focused-source EM survey versus time-domain and frequency-domain CSEM
Oskooi et al. Resolution of airborne VLF data
Malik et al. Ground‐Penetrating Radar Investigations along Hajipur Fault: Himalayan Frontal Thrust—Attempt to Identify Near Subsurface Displacement, NW Himalaya, India
Ezersky et al. Integrated study of the sinkhole development site on the Western shores of the Dead Sea using geophysical methods
Huang et al. A novel apparent resistivity for land-based controlled-source electromagnetic method
Wang et al. 3D pre-stack reverse time migration of ground penetrating radar for subsurface imaging
Yao et al. Detecting fracture zones in granite with seismic and audio-magnetotelluric methods

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent of invention or patent application
CB02 Change of applicant information

Address after: 710075 Xi'an high tech Industrial Development Zone, Shaanxi, Jin Jin Road, No. 82

Applicant after: Xi'an Research Institute, CCTEG

Address before: 710075 Xi'an high tech Industrial Development Zone, Shaanxi, Jin Jin Road, No. 82

Applicant before: China Coal Technology & Engineering Group Xi'an Research Institute

COR Change of bibliographic data

Free format text: CORRECT: APPLICANT; FROM: CHINA COAL TECHNOLOGY + ENGINEERING GROUP XI AN RESEARCH INSTITUTE TO: XI AN RESEARCH INSTITUTE OF CHINA COAL TECHNOLOGY + ENGINEERING GROUP CO., LTD.

C14 Grant of patent or utility model
GR01 Patent grant