CN113156519B - 一种音频大地电磁阵列高效施工勘探方法 - Google Patents

一种音频大地电磁阵列高效施工勘探方法 Download PDF

Info

Publication number
CN113156519B
CN113156519B CN202110448574.5A CN202110448574A CN113156519B CN 113156519 B CN113156519 B CN 113156519B CN 202110448574 A CN202110448574 A CN 202110448574A CN 113156519 B CN113156519 B CN 113156519B
Authority
CN
China
Prior art keywords
field component
station
electric field
measuring point
measuring
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
CN202110448574.5A
Other languages
English (en)
Other versions
CN113156519A (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 Huadong 814 Geophysical Exploration Co ltd Team 814 East China Geological Survey Bureau Of Nonferrous Metals Jiangsu Province
Original Assignee
Jiangsu Huadong 814 Geophysical Exploration Co ltd Team 814 East China Geological Survey Bureau Of Nonferrous Metals Jiangsu Province
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 Huadong 814 Geophysical Exploration Co ltd Team 814 East China Geological Survey Bureau Of Nonferrous Metals Jiangsu Province filed Critical Jiangsu Huadong 814 Geophysical Exploration Co ltd Team 814 East China Geological Survey Bureau Of Nonferrous Metals Jiangsu Province
Priority to CN202110448574.5A priority Critical patent/CN113156519B/zh
Publication of CN113156519A publication Critical patent/CN113156519A/zh
Application granted granted Critical
Publication of CN113156519B publication Critical patent/CN113156519B/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
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • 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)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Electromagnetism (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

本发明公开了一种音频大地电磁阵列高效施工勘探方法，包括如下步骤：S1.垂直构造走向部署作为测线方向；S2.在测区内建立观测主站，主站部署矢量方式，主站两侧建立辅站，辅站沿测线方向部署标量方式；S3.主站测点沿X方向和Y方向分别部署磁道并通过磁道进行磁场观测，测点还进行电道观测，辅站共用主站磁道；S4.利用测点处观测的磁场值计算磁场分量，利用磁场分量和电场分量计算测点处的卡尼亚视电阻率；S5.借用主站测点的电场分量，对辅站标量的原始时间序列进行处理，每个辅站分解计算生成两个正常测点，两个正常测点累加计算生成加密点。本发明可提高仪器的使用率，并提高音频大地电磁阵列的施工效率和精度。

Description

一种音频大地电磁阵列高效施工勘探方法
技术领域
本发明涉及地球物理勘探技术领域，具体涉及一种音频大地电磁阵列高效施工勘探方法。
背景技术
频率域电磁测深在深部地质构造的研究方面具有很好的勘探效果，其中频率为音频范围的频率域电磁测深由于具有勘探深度适中、勘探时间短的特点在矿产及工程方面应用较广，适中的勘探深度为1-3km左右，勘探较短时间一般为1-3小时，根据场源的不同又可以分为可控源音频大地电磁和音频大地电磁。
音频大地电磁仪较可控源音频大地电磁仪轻便，同时不存在可控源音频大地电磁的“近场”影响以及收发极之间的“阴影”效应，因此理论上应该更加适用。但在实际生产中，可控源音频大地电磁应用更广。其原因一方面在于可控源音频大地电磁抗干扰能力更强，另外一方面在于可控源音频大地电磁采用阵列式采集，具有施工效率高的特点。
现有技术中，由于实际生产中投入的音频大地电磁仪器少，无法像大地电磁阵列一样进行阵列式布极，而且目前的音频大地电磁计算软件均不支持标量计算，所以音频大地电磁无法像可控源音频大地电磁一样只测量一个方向的电场分量，这在一定程度造成了勘探精度的降低和施工效率不高的特点。而音频大地电磁如果想提高效率，必须进行标量采集。在早期大地电磁阵列方法中，软件支持标量采集，但后期仪器升级更新不再支持标量采集，在大地电磁阵列勘探中也存在不少对原始时间序列数据进行重新读取计算的方法，但其目标侧重于去除噪声或进行静态效应校正上。
发明内容
为解决上述问题，本发明提出了一种音频大地电磁阵列高效施工勘探方法，提高仪器的使用率，并提高音频大地电磁阵列的施工效率和精度。
本发明采用的技术方案是：
本申请提出一种音频大地电磁阵列高效施工勘探方法，包括如下步骤：
S1.垂直构造走向部署作为测线方向，以测线方向为X方向，以垂直测线方向为Y方向；
S2.在测区内建立观测主站，主站呈十字形部署矢量方式并在主站建立测点，主站两侧建立辅站，辅站沿测线方向部署标量方式；
S3.主站测点沿X方向和Y方向分别部署磁道并通过磁道进行磁场观测，测点还进行电道观测，辅站共用主站磁道；
S4.利用测点处观测的磁场值计算磁场分量，利用磁场分量和电场分量计算测点处的卡尼亚视电阻率；
S5.借用主站测点的电场分量，对辅站标量的原始时间序列进行处理，每个辅站分解计算生成两个正常测点，两个正常测点累加计算生成加密点；
步骤S5具体包括如下步骤：
S51.辅站对标量原始时间序列进行重新读取计算时需借用主站的Y向电场分量，即垂直方向辅站的卡尼亚视电阻率ρyx与主站测点垂直方向的卡尼亚视电阻率ρyx相同；
S52.每个辅站可分解计算出两个正常测点，将主站原始信号的X向电场分量更换为相同时间段辅站的其中一个正常测点的X向电场分量，其中一个正常测点的卡尼亚视电阻率为：
其中，Ex1为辅站其中一个正常测点处的X向电场分量，ρxy12极距为该正常测点至辅站的距离a；
S53.将主站原始信号的X向电场分量更换为相同时间段辅站的另一个正常测点的X向电场分量，另一个正常测点的卡尼亚视电阻率为：
其中，Ex2为辅站该正常测点处的X向电场分量，ρxy10极距为该正常测点至辅站的距离a；
S54.将主站原始信号的X向电场分量更换为相同时间段辅站的两个正常测点的X向电场分量的累加和，生成加密点，加密点的卡尼亚视电阻率为：
其中，ρxy11极距为2a。
优先地，步骤S2中，主站呈十字形进行矢量正常布极，分别部署N、S、E、W四个接线柱，N向接线柱接X正方向，S向接线柱接X负方向，E向接线柱接Y正方向，W向接线柱接Y负方向。
优先地，步骤S2中，辅站呈一字型进行标量布极，分别部署N1、S1、E1、W1四个接线柱，S向接线柱和E向接线柱接地，N向接线柱接X正方向，W向接线柱接X负方向。
优先地，步骤S4中，主站测点处的卡尼亚视电阻率为：
其中，f为发射电磁波的频率，Ex和Ey分别为测点处的X向电场分量和Y向电场分量，Hx和Hy分别为测点处的X向磁场分量和Y向磁场分量。
优先地，辅站数量至少为一个，用于提高工作效率。
优先地，沿测线方向的接线柱首尾相连或断开连接。
本发明的有益效果是：
1.辅站通过标量部署及对原始时间序列进行重新读取计算，使辅站由一台仪器只能作一个点扩大到分解为两个正常测点和一个加密点，提高了仪器的使用率，极大提高了音频大地电磁阵列的施工效率和精度；
2.一个主站可配置多个辅站，提高工作效率。
附图说明
附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例一起用于解释本发明，并不构成对本发明的限制。在附图中：
图1是本发明的音频大地电磁阵列的数据采集示意图；
图2是本发明的时间序列数据处理流程图；
图3是本发明的实际应用音频大地电磁阵列数据处理断面图。
图中附图标记为：1.AMT仪器，2.接线柱，3.磁棒。
具体实施方式
本申请提出一种音频大地电磁阵列高效施工勘探方法，包括如下步骤：
如图1所示，S1.垂直构造走向部署作为测线方向，以测线方向为X方向，以垂直测线方向为Y方向。
如图1所示，S2.在测区内建立观测主站，主站呈十字形部署矢量方式并在主站建立测点SITE14，主站两侧建立辅站，辅站沿测线方向部署标量方式。主站和辅站均设立AMT仪器1，主站呈十字形进行矢量正常布极，分别部署N、S、E、W四个接线柱2，N向接线柱2接X正方向，S向接线柱2接X负方向，E向接线柱2接Y正方向，W向接线柱2接Y负方向。辅站呈一字型进行标量布极，分别部署N、S、E、W四个接线柱2，S向接线柱2和E向接线柱2接地，N向接线柱2接X正方向，W向接线柱2接X负方向。
如图1所示，S3.主站测点SITE14沿X方向和Y方向分别部署磁道并通过磁道进行磁场观测，并通过磁棒3检测磁场分量，测点SITE14还进行电道观测，辅站共用主站磁道。
如图1所示，S4.利用测点SITE14处观测的磁场值计算磁场分量，利用磁场分量和电场分量计算测点SITE14处的卡尼亚视电阻率。主站测点SITE14处的卡尼亚视电阻率为：
其中，f为发射电磁波的频率，Ex和Ey分别为测点SITE14处的X向电场分量和Y向电场分量，Hx和Hy分别为测点SITE14处的X向磁场分量和Y向磁场分量。
如图1-2所示，S5.借用主站测点SITE14的电场分量，对辅站标量的原始时间序列进行处理，每个辅站分解计算生成两个正常测点，两个正常测点累加计算生成加密点。
具体包括如下步骤：
S51.辅站对标量原始时间序列进行重新读取计算时需借用主站的Y向电场分量，即垂直方向辅站的卡尼亚视电阻率ρyx与主站测点垂直方向的卡尼亚视电阻率ρyx相同；
S52.每个辅站可分解计算出两个正常测点，将主站原始信号的X向电场分量更换为相同时间段辅站的其中一个正常测点SITE12的X向电场分量，其中一个正常测点SITE12的卡尼亚视电阻率为：
其中，Ex1为辅站其中一个正常测点SITE12处的X向电场分量，ρxy12极距为该正常测点SITE12 X向电场分量Ex1的距离a；
S53.将主站原始信号的X向电场分量更换为相同时间段辅站的另一个正常测点SITE10的X向电场分量，另一个正常测点SITE10的卡尼亚视电阻率为：
其中，Ex2为辅站该正常测点SITE10处的X向电场分量，ρxy10极距为该正常测点SITE10 X向电场分量Ex2的距离a；
S54.将主站原始信号的X向电场分量更换为相同时间段辅站的两个正常测点的X向电场分量的累加和，生成加密点SITE11，加密点SITE11的卡尼亚视电阻率为：
其中，ρxy11极距为两个正常测点SITE12和SITE10的电场分量的累加和Ex1+Ex2的距离2a。
如图1所示，辅站数量至少为一个，用于提高工作效率，一般一个主站配备1-2个辅站的工作效率较高。
如图1所示，沿测线方向的接线柱2首尾相连或断开连接，较为灵活。断开连接状态更便于施工。
辅站的标量采集除了提高了施工的效率，还对困难地形地区施工特别适宜。本申请施工实例如下：2011年在秦岭地区国家某引水隧道工程进行了音频大地电磁法(AMT)勘探工作，10条线，长27km，点距50m。采取一个主站和一个辅站方式，主站仪器布置矢量，辅站仪器布置标量，辅站计算时借用主站垂直分量并利用程序使之分解为两个点，该次未使用加密点。如图3所示为某线处理断面图，高阻为花岗闪长岩区，低阻为含水区或断裂破碎带，与地表地质及实际验证情况较为吻合。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (6)

1.一种音频大地电磁阵列高效施工勘探方法，其特征在于：包括如下步骤：
S1.垂直构造走向部署作为测线方向，以测线方向为X方向，以垂直测线方向为Y方向；
S2.在测区内建立观测主站，主站呈十字形部署矢量方式并在主站建立测点，主站两侧建立辅站，辅站沿测线方向部署标量方式；
S3.主站测点沿X方向和Y方向分别部署磁道并通过磁道进行磁场观测，测点还进行电道观测，辅站共用主站磁道；
S4.利用测点处观测的磁场值计算磁场分量，利用磁场分量和电场分量计算测点处的卡尼亚视电阻率；
S5.借用主站测点的电场分量，对辅站标量的原始时间序列进行处理，每个辅站分解计算生成两个正常测点，两个正常测点累加计算生成加密点；
其中，步骤S5的具体步骤为：
S51.辅站对标量原始时间序列进行重新读取计算时需借用主站的Y向电场分量，即垂直方向辅站的卡尼亚视电阻率与主站测点垂直方向的卡尼亚视电阻率/>相同；
S52.每个辅站可分解计算出两个正常测点，将主站原始信号的X向电场分量更换为相同时间段辅站的其中一个正常测点的X向电场分量，其中一个正常测点的卡尼亚视电阻率为：
；
；
其中，为其中一个正常测点测线方向的视电阻率，/>为其中一个正常测点垂直测线方向的视电阻率，/>为发射电磁波的频率，/>为测点处的Y向电场分量，/>和/>分别为测点处的X向磁场分量和Y向磁场分量，/>为辅站其中一个正常测点处的X向电场分量，/>极距为该正常测点X向电场分量的距离a；
S53.将主站原始信号的X向电场分量更换为相同时间段辅站的另一个正常测点的X向电场分量，另一个正常测点的卡尼亚视电阻率为：
；
；
其中，为另一个正常测点测线方向的视电阻率，/>为另一个正常测点垂直测线方向的视电阻率，/>为发射电磁波的频率，/>为测点处的Y向电场分量，/>和/>分别为测点处的X向磁场分量和Y向磁场分量，/>为辅站该正常测点处的X向电场分量，/>极距为该正常测点X向电场分量的距离a；
S54.将主站原始信号的X向电场分量更换为相同时间段辅站的两个正常测点的X向电场分量的累加和，生成加密点，加密点的卡尼亚视电阻率为：
；
；
其中，为加密点测线方向的视电阻率，/>为加密点垂直测线方向的视电阻率，极距为两个正常测点的电场分量的累加和的距离2a。
2.根据权利要求1所述的音频大地电磁阵列高效施工勘探方法，其特征在于：步骤S2中，主站呈十字形进行矢量正常布极，分别部署N、S、E、W四个接线柱，N向接线柱接X正方向，S向接线柱接X负方向，E向接线柱接Y正方向，W向接线柱接Y负方向。
3.根据权利要求2所述的音频大地电磁阵列高效施工勘探方法，其特征在于：步骤S2中，辅站呈一字型进行标量布极，分别部署N1、S1、E1、W1四个接线柱，S向接线柱和E向接线柱接地，N向接线柱接X正方向，W向接线柱接X负方向。
4.根据权利要求1所述的音频大地电磁阵列高效施工勘探方法，其特征在于：步骤S4中，主站测点处的卡尼亚视电阻率为：
；
；
其中，为发射电磁波的频率，/>和/>分别为测点处的X向电场分量和Y向电场分量，和/>分别为测点处的X向磁场分量和Y向磁场分量，/>为测线方向的视电阻率，/>为垂直测线方向的视电阻率。
5.根据权利要求1所述的音频大地电磁阵列高效施工勘探方法，其特征在于：辅站数量至少为一个，用于提高工作效率。
6.根据权利要求3所述的音频大地电磁阵列高效施工勘探方法，其特征在于：沿测线方向的接线柱首尾相连或断开连接。
CN202110448574.5A 2021-04-25 2021-04-25 一种音频大地电磁阵列高效施工勘探方法 Active CN113156519B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110448574.5A CN113156519B (zh) 2021-04-25 2021-04-25 一种音频大地电磁阵列高效施工勘探方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110448574.5A CN113156519B (zh) 2021-04-25 2021-04-25 一种音频大地电磁阵列高效施工勘探方法

Publications (2)

Publication Number Publication Date
CN113156519A CN113156519A (zh) 2021-07-23
CN113156519B true CN113156519B (zh) 2023-08-11

Family

ID=76870321

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110448574.5A Active CN113156519B (zh) 2021-04-25 2021-04-25 一种音频大地电磁阵列高效施工勘探方法

Country Status (1)

Country Link
CN (1) CN113156519B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113759426B (zh) * 2021-08-02 2022-09-02 山东大学 一种基于参考道的人工源电磁勘探方法及系统

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5043667A (en) * 1989-04-21 1991-08-27 Amoco Corporation Method of magnetotelluric exploration using areal arrays
AU2002328690B2 (en) * 2001-10-11 2007-10-25 Cgg Data Services Ag Airborne geophysical measurements
CN102183792A (zh) * 2011-03-01 2011-09-14 吉林大学 人工源频率域电磁测深装置及测量方法
CN103809204A (zh) * 2012-11-06 2014-05-21 核工业北京地质研究院 一种野外音频大地电磁的数据采集方法
JP2014159991A (ja) * 2013-02-19 2014-09-04 National Institute Of Information & Communication Technology 複数の電磁波源からの電磁波の比吸収率の測定方法
CN104656156A (zh) * 2013-11-22 2015-05-27 中国石油天然气集团公司 音频大地电磁测深三维采集资料的磁参考处理方法
CN205384379U (zh) * 2015-12-22 2016-07-13 张雷 一种同步阵列电法仪系统
CN108897050A (zh) * 2018-02-28 2018-11-27 长江大学 一种可控源音频大地电磁测深法最佳收发距确定方法及系统
CN110488365A (zh) * 2019-09-11 2019-11-22 湖南五维地质科技有限公司 一种多极化大地电磁测深方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110133727A (zh) * 2019-06-10 2019-08-16 周丹 超深油气储层高精度电磁频谱探测方法
CN110531422B (zh) * 2019-07-25 2021-04-02 中国科学院地质与地球物理研究所 一种张量人工源电磁信号数据采集处理方法及装置

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5043667A (en) * 1989-04-21 1991-08-27 Amoco Corporation Method of magnetotelluric exploration using areal arrays
AU2002328690B2 (en) * 2001-10-11 2007-10-25 Cgg Data Services Ag Airborne geophysical measurements
CN102183792A (zh) * 2011-03-01 2011-09-14 吉林大学 人工源频率域电磁测深装置及测量方法
CN103809204A (zh) * 2012-11-06 2014-05-21 核工业北京地质研究院 一种野外音频大地电磁的数据采集方法
JP2014159991A (ja) * 2013-02-19 2014-09-04 National Institute Of Information & Communication Technology 複数の電磁波源からの電磁波の比吸収率の測定方法
CN104656156A (zh) * 2013-11-22 2015-05-27 中国石油天然气集团公司 音频大地电磁测深三维采集资料的磁参考处理方法
CN205384379U (zh) * 2015-12-22 2016-07-13 张雷 一种同步阵列电法仪系统
CN108897050A (zh) * 2018-02-28 2018-11-27 长江大学 一种可控源音频大地电磁测深法最佳收发距确定方法及系统
CN110488365A (zh) * 2019-09-11 2019-11-22 湖南五维地质科技有限公司 一种多极化大地电磁测深方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
可控源音频大地电磁测深法在永定河隧道勘察中的应用;陈煊;胡树林;郭志;李冬泉;刘金彪;;工程勘察(第01期);第81-86页 *

Also Published As

Publication number Publication date
CN113156519A (zh) 2021-07-23

Similar Documents

Publication Publication Date Title
CN100480734C (zh) 一种高分辨率去静态频率域大地电磁法
CN102788991B (zh) 基于z分量的透射槽波高密度快速探测方法
Liu et al. Three-dimensional seismic ahead-prospecting method and application in TBM tunneling
AU2020100984A4 (en) A Method and Apparatus for Ground-tunnel Wide Field Electromagnetic Surveying
CN104656157B (zh) 一种识别页岩气甜点区的方法及装置
CN101216566B (zh) 一种适合于准远区的电偶源频率域电测深方法及装置
CN113156519B (zh) 一种音频大地电磁阵列高效施工勘探方法
CN102176064A (zh) 煤矿井下瞬变电磁仪去除电磁干扰信号的方法
CN110879418B (zh) 一种频率域电性源电磁场视电导率计算方法
CN104656156A (zh) 音频大地电磁测深三维采集资料的磁参考处理方法
CN105301663A (zh) 时频电磁勘探数据空中采集装置及系统
CN104977619A (zh) 一种四方位全覆盖三维可控源电磁数据采集方法
CN103424769A (zh) 一种采空区联合多波地震勘探方法
CN102720484A (zh) 一种随钻声波测井仪器及测井方法
CN102998698A (zh) 一种束状海底电缆宽方位角地震数据采集方法
CN102608665A (zh) 物探时间域瞬变电磁系统测量数据资料精细处理技术
CN102998700B (zh) 一种垂直地震剖面资料的地表一致性能量补偿方法
CN104459789A (zh) 一种地震波初至拾取的方法
CN205139386U (zh) 时频电磁勘探数据空中采集装置及系统
CN105093270A (zh) 黄土塬区沿沟激发不规则煤田三维地震勘探技术
CN104502986B (zh) 物探激电测深数据层析法处理方法
CN111323825B (zh) 一种地-坑广域电磁探测方法及装置
CN105204074A (zh) 一种矿用本安型聚焦双频激电法激电信号接收系统
CN112882124A (zh) 一种三维铀成矿有利部位的圈定方法
Zhao et al. Preliminary structure framework and concealed intrusive rocks in the Qingchengzi ore field in northeast China disclosed by large-scale two-dimensional audio-magnetotelluric sounding

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