[go: up one dir, main page]

CN106679795B - 电磁探测噪声测量系统及降噪方法 - Google Patents

电磁探测噪声测量系统及降噪方法 Download PDF

Info

Publication number
CN106679795B
CN106679795B CN201710039321.6A CN201710039321A CN106679795B CN 106679795 B CN106679795 B CN 106679795B CN 201710039321 A CN201710039321 A CN 201710039321A CN 106679795 B CN106679795 B CN 106679795B
Authority
CN
China
Prior art keywords
noise
signal
electrode
receiving
electromagnetic detection
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
CN201710039321.6A
Other languages
English (en)
Other versions
CN106679795A (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.)
Beijing University of Technology
Original Assignee
Beijing University 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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201710039321.6A priority Critical patent/CN106679795B/zh
Publication of CN106679795A publication Critical patent/CN106679795A/zh
Application granted granted Critical
Publication of CN106679795B publication Critical patent/CN106679795B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H11/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties
    • G01H11/06—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties by electric means
    • 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

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Noise Elimination (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明提供一种电磁探测噪声测量系统及降噪方法，该系统包括一个与发射机连接的发射电极和与发射电极距离预设测量间距处的接收电极组，每个接收电极均与信号采集装置连接；还包括设置在相邻接收电极之间的噪声测量电极，噪声测量电极与信号采集装置连接；发射电极与接收电极处于同一测线方向上；噪声测量电极的布设方向与测线方向正交。本发明的电磁探测噪声测量系统及降噪方法，通过在相邻接收电极之间布设噪声测量电极，使发射电极与接收电极处于同一测线方向上，噪声测量电极的布设方向与测线方向正交，达到既有效采集接收端接收信号，又能采集接收端的背景噪声信号，后续根据预设降噪模型对各种信号进行降噪处理，得到所需的真实场信号。

Description

电磁探测噪声测量系统及降噪方法
技术领域
本发明涉及地球物理勘探技术领域，尤其涉及一种电磁探测噪声测量系统及降噪方法。
背景技术
电磁探测具有成本低、方便快捷、分辨率高的特点，尤其适用于金属矿藏探测，目前已成为地球物理勘探领域中的研究热点。
然而一方面在电磁探测中有效信号的强度随着收发距离的增加迅速衰减，在接收端电场强度一般只能达到mV级别；另一方面，随着经济社会的发展，各种背景噪声越来越严重，在部分人文活动频繁的工区，尖峰干扰与电力干扰强度极大，有效信号往往被这些背景噪声所淹没。
目前常见的降噪手段是增大发射功率、多次发射进行数据叠加、使用工频陷波等数字信号处理手段进行噪声滤除。但对于大收发距情况下，增大发射功率带来的效果微乎其微。多次发射进行数据叠加只能对随机噪声进行压制，无法滤除人文噪声，而当噪声含量复杂时，一般的数字滤波技术也很难取得较好的降噪效果。
发明内容
本发明提供一种电磁探测噪声测量系统及降噪方法，用于解决现有技术中降噪效果不好的问题。
第一方面，本发明提供一种电磁探测噪声测量系统，包括：一个与发射机连接的发射电极和与所述发射电极距离预设测量间距处的接收电极组，所述接收电极组包括多个接收电极，每个接收电极均与信号采集装置连接；还包括设置在相邻接收电极之间的噪声测量电极，所述噪声测量电极与信号采集装置连接；
其中，所述发射电极与所述接收电极处于同一测线方向上；所述噪声测量电极的布设方向与所述测线方向正交。
可选地，所述信号采集装置包括示波器。
可选地，所述预设测量间距为4-8km。
可选地，所述发射电极的布设距离为300-600m。
可选地，所述信号采集装置内设存储单元。
第二方面，本发明提供一种电磁探测信号降噪方法，包括：
获取电磁探测区域接收端处的混合信号和噪声信号，所述混合信号包括偶极子场信号和噪声信号，所述混合信号和所述噪声信号均为采用上述系统采集到的信号；
根据混合信号、噪声信号和预设的降噪模型获得对应于所述电磁探测接收端处的偶极子场信号。
可选地，所述降噪模型包括：
噪声信号的计算公式为：mn(t)＝hn(t)*n(t)；
混合信号的计算公式为：
mr(t)＝hr(t)*(n(t)+s(t))＝hr(t)*n(t)+hr(t)*s(t)；
降噪公式为：mr(t)-mn(t)＝hr(t)*n(t)-hn(t)*n(t)+hr(t)*s(t)，其中，hn(t)＝hr(t)，则mr(t)-mn(t)＝hr(t)*s(t)＝hs(t)；
反褶积公式为：f(t)*hs(t)＝f(t)*hr(t)*s(t)＝δ(t)*s(t)＝s(t)；
其中，mn(t)为噪声信号，hn(t)为噪声测量电极的冲激响应，mr(t)为混合信号，hr(t)是所用接收电极的冲激响应，n(t)为真实的噪声信号，mr(t)为混合信号，s(t)为真实的偶极子场信号，*为卷积符号，f(t)为反算子，δ(t)为单位冲激信号。
由上述技术方案可知，本发明的电磁探测噪声测量系统及降噪方法，通过在相邻接收电极之间布设噪声测量电极，使发射电极与接收电极处于同一测线方向上，噪声测量电极的布设方向与测线方向正交，达到既有效采集接收端接收信号，又能采集接收端的背景噪声信号，后续根据预设降噪模型对各种信号进行降噪处理，得到所需的真实场信号。
附图说明
图1为本发明实施例1提供的电磁探测噪声测量系统的结构分布示意图；
图2为本发明实施例提供的偶极子场的分布示意图；
图3为接收电极采集的信号示意图；
图4为噪声测量电极采集的噪声信号示意图；
图5为降噪后接收信号的示意图；
图6为降噪前后辨识出的大地脉冲响应曲线示意图；
图7为本发明实施例2提供的电磁探测信号降噪方法的流程示意图。
具体实施方式
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
本发明提供一种电磁探测噪声测量系统，包括：一个与发射机连接的发射电极和与所述发射电极距离预设测量间距处的接收电极组，所述接收电极组包括多个接收电极，每个接收电极均与信号采集装置连接；还包括设置在相邻接收电极之间的噪声测量电极，所述噪声测量电极与信号采集装置连接；其中，所述发射电极与所述接收电极处于同一测线方向上；所述噪声测量电极的布设方向与所述测线方向正交。
如图1所示实施例1提供一种电磁探测噪声测量系统，包括一个与发射机连接的发射电极和与该发射电极距离预设测量间距处的接收电极组，所述接收电极组包括两个接收电极，每个接收电极均与信号采集装置连接。
本系统还包括设置在该两个接收电极之间的噪声测量电极，所述噪声测量电极与信号采集装置连接。
其中，所述发射电极与所述接收电极处于同一测线方向上；所述噪声测量电极的布设方向与所述测线方向正交。
在本发明实施例中，由于电磁探测中的收发距离一般长达数公里，故在本实施例中预设测量间距可为4-8km。所述发射电极的布设距离可为300-600m。该布设距离相对于收发距离来说，所述的发射电极便可等效于一对电偶极子。如图2所示，偶极子场沿轴线呈对称分布，即在轴线方向上存在电势差。这也是噪声测量电极所测量的电压信号，而在轴线两侧(正交方向)的对称位置处为等电势点，电势差为零。本发明所述噪声测量电极的测量点为图2中的等电势点，接收电极的测量点为图2中的电势差点。
在理想状态下，由于等电势点间无电势差，噪声测量电极观测不到任何信号，但实际上，由于存在背景噪声，且其分布不具备对称性，因此噪声测量电极所记录的信号即为全部的背景噪声。
在本发明实施例中，所述信号采集装置可为示波器，用于采集电压信号。
在本发明实施例中，所述信号采集装置可内设存储单元，具有存储采集到的数据的功能。
由于在电磁探测过程中，接收电极采集电压信号并存储，噪声测量电极采集电压信号并存储。噪声测量电极采集的是噪声信号，接收电极采集的是混合信号，即包括偶极子场信号和噪声信号。因此，要得到最终的电磁探测信号需要对获取到的噪声信号及混合信号进行处理。
因此，在处理过程中，各信号采集装置需要将采集到的电压信号发送给用于进行信号处理的处理器中。
该处理器获得信号后，会按照预设的处理步骤对信号进行处理，最终得到所需的电测探测信号。
下面对信号处理步骤进行解释说明：
首先由于偶极子场的分布特性，噪声测量电极布设方向必须正交于测线方向，这样才能保证所测数据全部为背景噪声。
在接收端包含两种信号：沿测线方向的偶极子场信号，即场信号s(t)以及沿正交方向的背景噪声信号n(t)。其中接收电极采集的是二者的混合信号s(t)+n(t)，噪声测量电极采集的是背景噪声信号n(t)。
由于任何测量装置都存在带宽、衰减以及波形畸变，因此测量装置记录到的信号波形往往存在失真，即“真实信号”与“记录信号”之间存在偏差。根据信号与系统理论，“真实信号”与测量装置的“冲激响应”的卷积就是“记录信号”。
为便于区分，用以下符号表示真实信号与记录信号：
s(t)-真实的场信号
n(t)-真实的噪声信号
mn(t)-记录的噪声信号
mr(t)-记录的混合信号
hn(t)-噪声测量电极的冲激响应
hr(t)-接收电极的冲激响应
*-卷积符号
根据信号与系统理论可知：
记录的噪声信号等于噪声测量电极冲激响应与真实噪声信号的卷积，即：
mn(t)＝hn(t)*n(t) (1)
同理，记录的混合信号为：
mr(t)＝hr(t)*(n(t)+s(t))＝hr(t)*n(t)+hr(t)*s(t) (2)
显然(2)式中hr(t)*n(t)为噪声分量，从mr(t)中减去mn(t)以削弱该噪声分量，即：
mr(t)-mn(t)＝hr(t)*n(t)-hn(t)*n(t)+hr(t)*s(t) (3)
可见，只有hn(t)＝hr(t)时，即噪声测量电极与接收电极的冲激响应相同时，噪声分量才能有效剔除，此时(3)式可写作：
mr(t)-mn(t)＝hr(t)*s(t) (4)
如(4)式中，其结果已不包含噪声分量，说明本发明所述噪声测量电极所记录的噪声信号通过信号处理后实现了噪声剔除。
然而(4)式中的hr(t)*s(t)并不是真实的场信号s(t)，直接使用会降低数据质量，因此还需通过反褶积从hr(t)*s(t)中恢复出s(t)。
记hs(t)＝hr(t)*s(t)。
反褶积的目的是从hs(t)中恢复出s(t)。
由于hr(t)是所用接收电极的冲激响应，是已知的。因此可求解一反算子f(t)，使得hr(t)*f(t)＝δ(t)，其中δ(t)为单位冲激信号，δ(t)与任何信号的卷积都等于该信号本身。
对f(t)与hs(t)作卷积运算，有：
f(t)*hs(t)＝f(t)*hr(t)*s(t)＝δ(t)*s(t)＝s(t) (5)
由此可求得真实的场信号s(t)，该信号可用于后期的数据处理且具备较高的信噪比。
具体实施例：为验证本发明所述系统的可行性，在保津高速和廊沧高速相交的西南处，省道S334北东侧进行了一次验证试验。测线南端距离文安县城约4km，测线南北端经纬度分别为38°52'26.15"北；116°24'50.86"东和38°54'13.02"北；116°21'22.65"东。本次试验发射电流30A，发射极布设距离为300m，接收极布设距离为50m，收发距离为4km。
由于测线紧邻省道，因此尖峰干扰严重，另外接收区有输电线及电力变压器，存在严重的电力干扰。图3为接收数据，有用信号几乎完全被噪声淹没。图4为本发明所述噪声测量电极测得的背景噪声数据。对图4数据进行降噪，见图5。
从图5可以看出，进行降噪能够有效压制噪声，降噪后的接收信号能够清楚分辨出发射源信号的包络。分别使用原始接收数据以及降噪后的接收数据对大地脉冲响应曲线进行辨识，辨识结果对比见图6。
可以看出经过本发明所述系统采集数据及对数据降噪后，大地脉冲响应曲线的精度与平滑度得到了提高，经计算，本实施例中脉冲响应曲线的辨识信噪比提高了40dB。
图7示出了本发明实施例2提供一种电磁探测信号降噪方法，包括：
S21、获取电磁探测区域接收端处的混合信号和噪声信号，所述混合信号包括偶极子场信号和噪声信号，所述混合信号和所述噪声信号均为采用上述实施例1所述系统采集到的信号；
S22、根据混合信号、噪声信号和预设的降噪模型获得对应于所述电磁探测接收端处的偶极子场信号。
所述降噪模型包括：
噪声信号的计算公式为：mn(t)＝hn(t)*n(t)；
混合信号的计算公式为：
mr(t)＝hr(t)*(n(t)+s(t))＝hr(t)*n(t)+hr(t)*s(t)；
降噪公式为：mr(t)-mn(t)＝hr(t)*n(t)-hn(t)*n(t)+hr(t)*s(t)，其中，hn(t)＝hr(t)，则mr(t)-mn(t)＝hr(t)*s(t)＝hs(t)；
反褶积公式为：f(t)*hs(t)＝f(t)*hr(t)*s(t)＝δ(t)*s(t)＝s(t)；
其中，mn(t)为噪声信号，hn(t)为噪声测量电极的冲激响应，mr(t)为混合信号，hr(t)是所用接收电极的冲激响应，n(t)为真实的噪声信号，mr(t)为混合信号，s(t)为真实的偶极子场信号，*为卷积符号，f(t)为反算子，δ(t)为单位冲激信号。
针对本实施例中的降噪模型在上述实施例1中有详细解释说明，故在此不对上述各种信息的降噪处理进行再说明。
本发明的电磁探测信号降噪方法，通过在相邻接收电极之间的噪声测量电极，使发射电极与接收电极处于同一测线方向上，噪声测量电极的布设方向与测线方向正交，达到既有效采集接收端接收信号，又能接收端的背景噪声信号，后续根据预设降噪模型对各种信号进行降噪处理，得到所需的真实场信号。
此外，本领域的技术人员能够理解，尽管在此所述的一些实施例包括其它实施例中所包括的某些特征而不是其它特征，但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如，在下面的权利要求书中，所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
应该注意的是上述实施例对本发明进行说明而不是对本发明进行限制，并且本领域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施例。在权利要求中，不应将位于括号之间的任何参考符号构造成对权利要求的限制。单词“包含”不排除存在未列在权利要求中的元件或步骤。位于元件之前的单词“一”或“一个”不排除存在多个这样的元件。本发明可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的单元权利要求中，这些装置中的若干个可以是通过同一个硬件项来具体体现。单词第一、第二、以及第三等的使用不表示任何顺序。可将这些单词解释为名称。
本领域普通技术人员可以理解：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明权利要求所限定的范围。

Claims (1)

1.一种电磁探测信号降噪方法，其特征在于，包括：
获取电磁探测区域接收端处的混合信号和噪声信号，所述混合信号包括偶极子场信号和噪声信号，所述混合信号和所述噪声信号均为采用电磁探测噪声测量系统采集到的信号；
根据混合信号、噪声信号和预设的降噪模型获得对应于所述电磁探测接收端处的偶极子场信号；
所述降噪模型包括：
噪声信号的计算公式为：mn(t)＝hn(t)*n(t)；
混合信号的计算公式为：
mr(t)＝hr(t)*(n(t)+s(t))＝hr(t)*n(t)+hr(t)*s(t)；
降噪公式为：mr(t)-mn(t)＝hr(t)*n(t)-hn(t)*n(t)+hr(t)*s(t)，其中，hn(t)＝hr(t)，则mr(t)-mn(t)＝hr(t)*s(t)＝hs(t)；
反褶积公式为：f(t)*hs(t)＝f(t)*hr(t)*s(t)＝δ(t)*s(t)＝s(t)；
其中，mn(t)为噪声信号，hn(t)为噪声测量电极的冲激响应，mr(t)为混合信号，hr(t)是所用接收电极的冲激响应，n(t)为真实的噪声信号，s(t)为真实的偶极子场信号，*为卷积符号，f(t)为反算子，δ(t)为单位冲激信号；
所述电磁探测噪声测量系统包括：一个与发射机连接的发射电极和与所述发射电极距离预设测量间距处的接收电极组，所述接收电极组包括多个接收电极，每个接收电极均与信号采集装置连接；还包括设置在相邻接收电极之间的噪声测量电极，所述噪声测量电极与信号采集装置连接；
其中，所述发射电极与所述接收电极处于同一测线方向上；所述噪声测量电极的布设方向与所述测线方向正交。
CN201710039321.6A 2017-01-18 2017-01-18 电磁探测噪声测量系统及降噪方法 Active CN106679795B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710039321.6A CN106679795B (zh) 2017-01-18 2017-01-18 电磁探测噪声测量系统及降噪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710039321.6A CN106679795B (zh) 2017-01-18 2017-01-18 电磁探测噪声测量系统及降噪方法

Publications (2)

Publication Number Publication Date
CN106679795A CN106679795A (zh) 2017-05-17
CN106679795B true CN106679795B (zh) 2023-10-03

Family

ID=58860743

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710039321.6A Active CN106679795B (zh) 2017-01-18 2017-01-18 电磁探测噪声测量系统及降噪方法

Country Status (1)

Country Link
CN (1) CN106679795B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4890264A (en) * 1988-03-21 1989-12-26 Atlantic Richfield Company Seismic exploration method and apparatus for cancelling non-uniformly distributed noise
CN104020496A (zh) * 2014-06-27 2014-09-03 吉林大学 一种轴向同线方式的地面可控源电磁勘探方法
WO2016105478A1 (en) * 2014-12-23 2016-06-30 Axell Wireless Ltd. Harmonizing noise aggregation and noise management in distributed antenna system
CN206563613U (zh) * 2017-01-18 2017-10-17 北京工业大学 电磁探测噪声测量系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0616784D0 (en) * 2006-08-24 2006-10-04 Mtem Ltd Reduction of noise in electrical field measurements
RU2425399C2 (ru) * 2007-05-08 2011-07-27 Общество с ограниченной ответственностью "Сибирская геофизическая научно-производственная компания" Способ морской электроразведки и устройство для морской электроразведки в движении судна
US10012746B2 (en) * 2014-02-24 2018-07-03 Cgg Services Sas Systems and methods for pre-correlation noise attentuation using random sweeps in seismic data processing

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4890264A (en) * 1988-03-21 1989-12-26 Atlantic Richfield Company Seismic exploration method and apparatus for cancelling non-uniformly distributed noise
CN104020496A (zh) * 2014-06-27 2014-09-03 吉林大学 一种轴向同线方式的地面可控源电磁勘探方法
WO2016105478A1 (en) * 2014-12-23 2016-06-30 Axell Wireless Ltd. Harmonizing noise aggregation and noise management in distributed antenna system
CN206563613U (zh) * 2017-01-18 2017-10-17 北京工业大学 电磁探测噪声测量系统

Also Published As

Publication number Publication date
CN106679795A (zh) 2017-05-17

Similar Documents

Publication Publication Date Title
CN103389514B (zh) 一种大地电磁去噪装置及大地电磁去噪方法
Ji et al. Noise reduction of time domain electromagnetic data: Application of a combined wavelet denoising method
CN105866852B (zh) 一种基于相关性检测的远参考大地电磁阻抗计算方法
Chun et al. Magnitude estimation and source discrimination: A close look at the 2006 and 2009 North Korean underground nuclear explosions
CN103033849B (zh) 多通道核磁共振地下水探测仪及其野外工作方法
CN106814402B (zh) 瞬变电磁信号叠前去噪方法
CN203759262U (zh) 瞬变电磁信号接收装置
CN106291719B (zh) 一种阵列人工源磁场频率测深方法
CN106199733B (zh) 用于评价地热富水地质异常体的方法及装置
Cao et al. Speeded‐up robust features based single‐ended travelling wave fault location: a practical case study in Yunnan power grid of China
CN105204065A (zh) 一种拾取初至波的方法和装置
CN102830170B (zh) 一种基于超声测试获取煤样横波信号的控制方法及装置
Lei et al. Anti-interference test for the new SEP instrument: CSAMT study at Dongguashan Copper Mine, China
WO2024255521A1 (zh) 一种钻孔瞬变电磁法三维展开处理方法及系统
CN103698673A (zh) 基于定向天线阵列的套管局部放电带电检测装置
CN104749643B (zh) 一种用于电磁干扰地区的csamt测量方法和装置
CN206563613U (zh) 电磁探测噪声测量系统
CN106679795B (zh) 电磁探测噪声测量系统及降噪方法
Tian et al. Permutation fuzzy entropy based ICEEMDAN de-noising for inertial sensors
Cousins Direction finding on whistlers and related VLF signals
Li et al. Time-domain spectral induced polarization based on pseudo-random sequence
US12117583B1 (en) Method and device for mobile rapid exploration of mineral resources
Kumar et al. Determination and identification of focal mechanism solutions for Himalayan earthquakes from waveform inversion employing ISOLA software
CN110058319A (zh) 一种大地电磁数据采集方法、装置及终端设备
JP5562237B2 (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