CN115586577A - 一种定源瞬变电磁非中心点观测数据全时转换方法 - Google Patents
一种定源瞬变电磁非中心点观测数据全时转换方法 Download PDFInfo
- Publication number
- CN115586577A CN115586577A CN202211149374.0A CN202211149374A CN115586577A CN 115586577 A CN115586577 A CN 115586577A CN 202211149374 A CN202211149374 A CN 202211149374A CN 115586577 A CN115586577 A CN 115586577A
- Authority
- CN
- China
- Prior art keywords
- point
- model
- observation
- conversion
- electromotive force
- 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.)
- Withdrawn
Links
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 69
- 230000001052 transient effect Effects 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000011160 research Methods 0.000 claims abstract description 7
- 238000005259 measurement Methods 0.000 claims abstract description 3
- 238000012545 processing Methods 0.000 claims abstract description 3
- 238000004364 calculation method Methods 0.000 claims description 19
- 238000003384 imaging method Methods 0.000 abstract description 4
- 238000004088 simulation Methods 0.000 abstract description 2
- 230000009466 transformation Effects 0.000 description 10
- 230000002159 abnormal effect Effects 0.000 description 8
- 238000001514 detection method Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 8
- 230000004044 response Effects 0.000 description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 7
- 230000002547 anomalous effect Effects 0.000 description 5
- 238000000418 atomic force spectrum Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 239000003245 coal Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000005672 electromagnetic field Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
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
-
- 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
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
-
- 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/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
- G01V2003/086—Processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Electromagnetism (AREA)
- Geophysics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开一种定源瞬变电磁非中心点观测数据全时转换方法，具体步骤包括：先依据研究区钻孔电阻率测井曲线构建层状地电模型，通过对层状地电模型进行数值模拟，获取源场观测空间各观测点转换系数，包括依据现场定源瞬变电磁数据采集方式设置模型线框大小、模型测点点距、模型观测和关断时间等；再提取包括中心点在内的各观测点全时感应电动势数据，按公式进行比值处理，获取转换系数并将其代入实测瞬变电磁数据中将非中心点感应电动势转换为中心点感应电动势。基于以上非中心点观测数据的转换方法，可大幅度提高定源回线装置内部观测面积，获得高覆盖比观测数据，同时改善了早时间段观测数据的质量，进而实现浅层地电结构的高精度成像。
Description
技术领域
本发明属于地球物理探测领域，具体涉及一种定源瞬变电磁非中心点观测数据全时转换方法。
背景技术
定源瞬变电磁法是一种时间域电磁法，其因效率高、抗干扰能力强、探测深度大等优点，被广泛应用于探测采空区、水文地质、资源调查等领域。近年来，该方法也逐渐应用于城市、交通、水利等浅部地层探测领域。
定源回线装置进行数据采集时，其观测数据覆盖范围一般为1/3回线面积，即在发射回线中心区域1/3范围内进行观测。然而由于该方法布线难度较大，为了一次布设采集更多数据，得到更高的工作效率，必须提高观测数据的覆盖比。此外，该方法观测数据大部分分布于非中心区域，而非中心数据无全时高精度解析解。
针对瞬变电磁全时视电阻率计算问题，目前已有许多学者进行了大量研究，提出了多项式逼近、连分式求根、逆样条插值、数值选代等方法，但是涉及定源瞬变电磁非中心数据全时视电阻率计算研究较少，以往一般采用中心方式的晚期近似处理方法，这种解释方法虽然简单，但早中期数据被舍弃，具有一定的浅部探测盲区，并且随着偏移距的增大，视电阻率计算误差也逐渐增大。2011年研究人员将非中心点各分量表示成与中心点垂直分量类似的多项式形式，采用待定系数法，完成了非中心点视电阻率定义，但是对于早期道，校正后的视电阻率值与真实值之间误差相对较大；2012年研究人员采用电偶极子积分法计算任意场点的电磁响应，并采用迭代法提取任意点全时视电阻率值，这种方法求解非中心位置视电阻率较为准确，但计算速度较慢；2015年研究人员利用贝塞尔函数展开法计算大回线源非中心点垂直分量与水平分量的电磁响应，通过分析非中心点处的多分量电磁响应特征，得出了一种基于磁感应强度的大回线源全域视电阻率计算方法：在中心点及框内偏移距较小处，直接利用插值方法算出视电阻率的唯一解，对于偏移距较大的非中心点，采用无约束点的最小曲率差分选代法求解，计算过程较为复杂。
综上可知，目前已提出的定源瞬变电磁非中心数据全时视电阻率计算方法大都存在计算过程复杂、速度慢、早期计算结果误差大等缺点。针对上述问题，现设计一种定源瞬变电磁非中心点观测数据全时转换方法。
发明内容
针对现有技术的不足，本发明的目的在于提供一种定源瞬变电磁非中心点观测数据全时转换方法，大幅提高定源回线装置观测区间大小，获得高覆盖比观测数据，提高其野外工作效率，并实现观测区浅层地电结构高精度成像。
本发明的目的可以通过以下技术方案实现：
一种定源瞬变电磁非中心点观测数据全时转换方法，具体包括如下步骤：
S1、收集研究区钻孔电阻率测井曲线；
S2、构建层状地质与地电模型；
S3、模型数值计算获得观测点感应电动势；
S4、提取非中心点转换系数；
S5、转换实测非中心点数据。
进一步的，具体转换系数获取方法如下：
①、依据定源瞬变电磁数据采集方式，设置模型线框大小；
②、依据定源瞬变电磁数据采集方式，设置模型测点点距；
③、依据定源瞬变电磁数据采集方式，设置模型观测和关断时间；
④、采用时域有限差分方法对各观测点进行正演计算；
⑤、提取包括中心点在内的各观测点全时感应电动势数据，并与中心点感应电动势数据按公式
进行比值处理，其结果作为转换系数，式中Ki(t)为转换系数，Vi(t)为模型非中心点感应电动势值，V0(t)为模型中心点感应电动势值，i为点号。
进一步的，转换系数代入实测瞬变电磁数据中进行非中心点感应电动势转换，转换公式为
式中Ki(t)为转换系数，Vi′(t)为实测非中心点感应电动势值，V′0(t)为实测非中心点转换后感应电动势值，i为点号。
本发明的有益效果：
1、本发明提出的定源瞬变电磁非中心点观测数据全时转换方法，可大幅度提高定源回线装置内部观测面积，获得高覆盖比观测数据，提高其野外工作效率；
2、本发明提出的定源瞬变电磁非中心点观测数据全时转换方法，可将非中心点感应电动势转换为中心点感应电动势，为非中心点感应电动势直接求解全时视电阻率提供了途径；
3、本发明提出的定源瞬变电磁非中心点观测数据全时转换方法，可明显改善定源瞬变电磁观测数据的质量，对早时间段观测数据质量的提升尤为突出，为实现浅层地电结构的高精度成像建立了基础。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例的整体流程示意图；
图2是本发明实施例的数值模拟正演流程示意图；
图3是本发明实施例的测点选取位置示意图；
图4是本发明实施例的均匀地电模型A通过正演计算获取的感应电动势曲线图；
图5是本发明实施例的层状地电模型B通过正演计算获取的感应电动势曲线图；
图6是本发明实施例的含异常体模型C中1号点感应电动势曲线转换图；
图7是本发明实施例的含异常体模型C中2号点感应电动势曲线转换图；
图8是本发明实施例的含异常体模型C中3号点感应电动势曲线转换图；
图9是本发明实施例的含异常体模型C中4号点感应电动势曲线转换图；
图10是本发明实施例的含异常体模型C中5号点感应电动势曲线转换图；
图11是本发明实施例的含异常体模型C中6号点感应电动势曲线转换图；
图12是本发明实施例的含异常体模型D中1号点感应电动势曲线转换图；
图13是本发明实施例的含异常体模型D中2号点感应电动势曲线转换图；
图14是本发明实施例的含异常体模型D中3号点感应电动势曲线转换图；
图15是本发明实施例的含异常体模型D中4号点感应电动势曲线转换图；
图16是本发明实施例的含异常体模型D中5号点感应电动势曲线转换图；
图17是本发明实施例的含异常体模型D中6号点感应电动势曲线转换图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
如图1所示，本发明的整体流程包括如下步骤：
S1、收集研究区钻孔电阻率测井曲线；
S2、构建层状地质与地电模型；
S3、模型数值计算获得观测点感应电动势；
S4、提取非中心点转换系数；
S5、转换实测非中心点数据。
具体实施例如下：
本发明以某矿定源瞬变电磁现场试验探测为实施例，该地区采煤地面塌陷导致生态环境问题频繁发生，需对其主采煤层上覆岩土层的含水层分布及其赋水性进行精准探测，如图2所示，具体实现方式如下：
(1)、在正演软件中依据研究区钻孔电阻率测井曲线构建层状地电模型B，各地层由上而下视电阻率分别为60Ω·m、25Ω·m、60Ω·m、180Ω·m、200Ω·m、800Ω·m、300Ω·m，厚度分别为20m、70m、60m、3m、20m、7m、90m，空间范围为600m×600m×300m；
(2)、依据定源瞬变电磁数据采集方式，设置模型线框大小为250m×250m；
(3)、依据定源瞬变电磁数据采集方式，设置模型测点点距为20m；
(4)、依据定源瞬变电磁数据采集方式，设置模型观测时间为0.00000586s～0.00930413s；
(5)、因瞬变场的对称性，模型仅选取6个非中心观测点(P1、P2、P3、P4、P5、P6)和1个中心观测点(P0)进行正演计算，如图3所示；
(6)、为对比层状地电模型与均匀地电模型的转换效果，以研究区第二层位地质与地电特征为参考，构建均匀地电模型A，设视电阻率为25Ω·m，空间范围为600m×600m×300m，并选取与层状地电模型B相同的观测点进行正演计算；
均匀地电模型A与层状地电模型B感应电动势曲线如图4、图5所示，可见，均匀和层状地层条件下非中心观测点早中期瞬变场响应存在误差，若用原始瞬变电磁数据直接进行计算反演，会导致探测结果误差较大，难以准确反映非中心点下方浅部异常体信息；
(7)、提取均匀地电模型A与层状地电模型B中6个非中心观测点(P1、P2、P3、P4、P5、P6)的转换系数，
式中Ki(t)为转换系数，Vi(t)为模型非中心点感应电动势值，V0(t)为模型中心点感应电动势值，i为点号；
(8)、针对实际现场探测中浅部地层存在富水异常情况构建两种待转换模型(含异常体模型C和含异常体模型D)，并选取与层状地电模型B相同的观测点进行正演计算，具体模型参数如下：
(8—1)、含异常体模型C在层状地电模型B的基础上于发射回线中心正下方第二层位设空间大小为60m×30m×70m的局部富水区，富水区视电阻率为0.625Ω·m；
(8—2)、含异常体模型D模拟了实际探测中偏移距较大的非中心点下方存在异常体情况，其在含异常体模型C的基础上于发射回线左下方第二层位增设一个空间大小为60m×40m×70m的局部富水区，富水区视电阻率为0.625Ω·m；
(9)、将转换系数分别代入两种待转换模型中，转换公式为
式中Ki(t)为转换系数，Vi′(t)为转换模型非中心点感应电动势值，V′0(t)为转换模型非中心点转换后感应电动势值，i为点号，对转换模型非中心点感应电动势进行转换；
(10)、通过移动回线测得两种待转换模型非中心点作中心点时实际感应电动势，并与转换模型非中心点转换后感应电动势进行对比，根据相对误差
(式中，λ为相对误差)大小来判断模型转换效果，两种转换模型转换结果对比如图6—图17所示。
分别对比分析图6—图11和图12—图17中均匀地电模型A与层状地电模型B对两种含异常体待转换模型的转换效果，综上可知：
两种转换模型对非中心点瞬变场响应曲线均有一定转换效果，在早期，由层状地电模型B转换后的V′0B(t)曲线较均匀地电模型A转换后的V′0A(t)曲线与V(t)曲线吻合度更高。从曲线来看，同一测点在早期同一时刻，层状地电模型B的转换误差λB小于均匀地电模型A的转换误差λA；两种转换模型对偏移距小的非中心点转换结果相似，对偏移距较大的非中心点转换结果差异明显；层状地电模型B转换最大误差值显著低于均匀地电模型A转换最大误差值。由结果可见，总体上层状地电模型B对早期非中心点瞬变场响应的转换效果优于均匀地电模型A。
实施案例效果表明，本发明方法有效转换了定源瞬变电磁法非中心点早期瞬变电磁场响应，此外，通过对比两种转换模型的转换效果，可以得出层状地电模型的转换效果明显优于均匀地电模型。因此，在采用定源瞬变电磁法探测相对平滑的地层时，通过建立研究区层状地电模型对实测数据进行转换，可大幅度提高定源回线装置内部观测面积，获得高覆盖比观测数据，提高其野外工作效率，同时改善了早时间段观测数据的质量，进而实现浅层地电结构的高精度成像。
在本说明书的描述中，参考术语“一个实施例”、“示例”、“具体示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。
Claims (3)
1.一种定源瞬变电磁非中心点观测数据全时转换方法，其特征在于，具体包括如下步骤：
S1、收集研究区钻孔电阻率测井曲线；
S2、构建层状地质与地电模型；
S3、模型数值计算获得观测点感应电动势；
S4、提取非中心点转换系数；
S5、转换实测非中心点数据。
2.根据权利要求1所述的一种定源瞬变电磁非中心点观测数据全时转换方法，其特征在于，具体转换系数获取方法如下：
①、依据定源瞬变电磁数据采集方式，设置模型线框大小；
②、依据定源瞬变电磁数据采集方式，设置模型测点点距；
③、依据定源瞬变电磁数据采集方式，设置模型观测和关断时间；
④、采用时域有限差分方法对各观测点进行正演计算；
⑤、提取包括中心点在内的各观测点全时感应电动势数据，并与中心点感应电动势数据按公式
进行比值处理，其结果作为转换系数，式中Ki(t)为转换系数，Vi(t)为模型非中心点感应电动势值，V0(t)为模型中心点感应电动势值，i为点号。
3.根据权利要求1所述的一种定源瞬变电磁非中心点观测数据全时转换方法，其特征在于，转换系数代入实测瞬变电磁数据中进行非中心点感应电动势转换，转换公式为
式中Ki(t)为转换系数，Vi′(t)为实测非中心点感应电动势值，V0′(t)为实测非中心点转换后感应电动势值，i为点号。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211149374.0A CN115586577A (zh) | 2022-09-21 | 2022-09-21 | 一种定源瞬变电磁非中心点观测数据全时转换方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211149374.0A CN115586577A (zh) | 2022-09-21 | 2022-09-21 | 一种定源瞬变电磁非中心点观测数据全时转换方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115586577A true CN115586577A (zh) | 2023-01-10 |
Family
ID=84778986
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211149374.0A Withdrawn CN115586577A (zh) | 2022-09-21 | 2022-09-21 | 一种定源瞬变电磁非中心点观测数据全时转换方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115586577A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117454675A (zh) * | 2023-12-26 | 2024-01-26 | 中国地质科学院地球物理地球化学勘查研究所 | 一种定源瞬变电磁响应校正方法、系统及设备 |
-
2022
- 2022-09-21 CN CN202211149374.0A patent/CN115586577A/zh not_active Withdrawn
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117454675A (zh) * | 2023-12-26 | 2024-01-26 | 中国地质科学院地球物理地球化学勘查研究所 | 一种定源瞬变电磁响应校正方法、系统及设备 |
| CN117454675B (zh) * | 2023-12-26 | 2024-04-09 | 中国地质科学院地球物理地球化学勘查研究所 | 一种定源瞬变电磁响应校正方法、系统及设备 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105785455B (zh) | 一种基于b样条插值的二维地面核磁共振反演方法 | |
| CN102508303B (zh) | 地下工程聚焦层析激发极化超前探测方法 | |
| CN104537714B (zh) | 磁共振与瞬变电磁空间约束联合反演方法 | |
| CN111927552A (zh) | 基于天空地勘察技术的复杂山区长大深埋隧道勘察方法 | |
| CN106022339B (zh) | 一种复垦土地浅埋地埋管深度的提取方法 | |
| CN112949134B (zh) | 基于非结构有限元方法的地-井瞬变电磁反演方法 | |
| CN102901985B (zh) | 一种适用于起伏地表的深度域层速度修正方法 | |
| Xue et al. | Extracting the virtual reflected wavelet from TEM data based on regularizing method | |
| Zhou et al. | Migration velocity analysis and prestack migration of common-transmitter GPR data | |
| CN113341467A (zh) | 基于多插值方法的矿井瞬变电磁三维显示方法 | |
| CN206378448U (zh) | 基于综合物探方法的污染土检测测线布置结构 | |
| CN108169802B (zh) | 一种粗糙介质模型的时域电磁数据慢扩散成像方法 | |
| Steiner et al. | Improved estimation of ice and water contents in Alpine permafrost through constrained petrophysical joint inversion: The Hoher Sonnblick case study | |
| CN115586577A (zh) | 一种定源瞬变电磁非中心点观测数据全时转换方法 | |
| CN114386464B (zh) | 一种瞬变电磁激电信息的深度学习提取方法 | |
| CN115292890A (zh) | 基于多源辅助数据开发的场地土壤污染物浓度三维空间预测方法 | |
| CN112666612B (zh) | 基于禁忌搜索的大地电磁二维反演方法 | |
| CN113406707A (zh) | 一种大地电磁多尺度、多时段探测方法 | |
| CN114896673B (zh) | 一种干热岩储层水力压裂改造体积预测方法 | |
| CN106483570A (zh) | 一种大地电磁场物探方法及装置 | |
| CN115857028A (zh) | 基于非接触观测法的水文地质参数的时空定量表达方法 | |
| CN112882124B (zh) | 一种三维铀成矿有利部位的圈定方法 | |
| CN114488327B (zh) | 基于地面基点的水平磁场与井中垂直磁场联合测量方法 | |
| CN105549101A (zh) | 一种瞬变电磁数据微分电导解释方法 | |
| Abbas et al. | Implementation of ground penetrating radar and electrical resistivity tomography for inspecting the Greco-Roman Necropolis at Kilo 6 of the Golden Mummies Valley, Bahariya Oasis, Egypt |
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 | ||
| WW01 | Invention patent application withdrawn after publication | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20230110 |