CN115079284B - 一种激电中梯异常提取方法 - Google Patents
一种激电中梯异常提取方法 Download PDFInfo
- Publication number
- CN115079284B CN115079284B CN202210661903.9A CN202210661903A CN115079284B CN 115079284 B CN115079284 B CN 115079284B CN 202210661903 A CN202210661903 A CN 202210661903A CN 115079284 B CN115079284 B CN 115079284B
- Authority
- CN
- China
- Prior art keywords
- induced
- measuring point
- polarization
- apparent
- ladder
- 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
Links
- 230000010287 polarization Effects 0.000 title claims abstract description 44
- 238000000605 extraction Methods 0.000 title claims description 8
- 230000002159 abnormal effect Effects 0.000 claims abstract description 29
- 238000005192 partition Methods 0.000 claims abstract description 20
- 230000000007 visual effect Effects 0.000 claims abstract description 20
- 230000005856 abnormality Effects 0.000 claims abstract description 13
- 238000000034 method Methods 0.000 claims abstract description 12
- 238000000638 solvent extraction Methods 0.000 claims abstract description 4
- 238000005259 measurement Methods 0.000 abstract description 5
- 230000005284 excitation Effects 0.000 abstract description 3
- 230000015572 biosynthetic process Effects 0.000 abstract description 2
- 230000002349 favourable effect Effects 0.000 abstract description 2
- 238000005065 mining Methods 0.000 description 7
- 239000011532 electronic conductor Substances 0.000 description 6
- JQJCSZOEVBFDKO-UHFFFAOYSA-N lead zinc Chemical compound [Zn].[Pb] JQJCSZOEVBFDKO-UHFFFAOYSA-N 0.000 description 5
- 239000002184 metal Substances 0.000 description 5
- 229910052751 metal Inorganic materials 0.000 description 5
- UCKMPCXJQFINFW-UHFFFAOYSA-N Sulphide Chemical compound [S-2] UCKMPCXJQFINFW-UHFFFAOYSA-N 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- QAOWNCQODCNURD-UHFFFAOYSA-L Sulfate Chemical compound [O-]S([O-])(=O)=O QAOWNCQODCNURD-UHFFFAOYSA-L 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
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/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- 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
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种激电中梯异常提取方法。该方法包括以下步骤：①确定勘测区域的测点，通过激电测量得到每个测点的视极化率和视电阻率；②对勘测区域进行分区，计算出每个分区内测点的视极化率平均值和视电阻率平均值；③结合每个测点的实测视电阻率，计算每个分区内的每个测点视极化率期望值；④结合每个测点的实测视极化率值，计算每个分区不同测点的衍生指数；⑤根据上步骤计算得到的每个分区不同测点的衍生指数绘制等值线平面图，在图中圈定出衍生指数异常区域，衍生指数异常区域即为激电中梯异常区域。本发明可以进一步缩小激电中梯异常范围，更加准确的圈定出最终异常地段，从而精确预测成矿的有利地段。
Description
技术领域
本发明涉及一种提取方法，具体地说是一种激电中梯异常提取方法。
背景技术
铅锌矿在有色金属工业中占有重要地位，是我国重要的战略性矿产资源。为保证我国铅锌工业可持续发展，长期以来不断加大对铅锌矿的找矿与勘探工作。随着地质勘查工作程度的逐步提高，浅部资源的发现越加困难，找矿难度日益增大，因此开展深部铅锌矿资源空间的勘查十分迫切。
由于铅锌矿主要以硫化物和硫酸盐类形式存在，导致矿石与围岩之间的电性存在差异，所以人们通过提取激电中梯异常来对地层或矿源进行评价，从而初步圈定金属矿区特别是含硫化物金属矿区，为后期找矿工作提供重要依据。
目前，激电中梯异常提取方法，多以视极化率和视电阻率为基础，提取视极化率异常，然后再结合视电阻率特征对地层、岩性、断裂进行评价。但是，该方法不能突出与电子导体相关的低阻极化体异常，且不能压制与电子导体不相关的激电异常，导致金属矿区圈定范围不准或者过大，使得后期找矿工作的不能顺利开展。
发明内容
本发明的目的就是提供一种激电中梯异常提取方法，以解决采用现有的激电中梯异常提出到的激电中梯不能突出于电子导体相关的低阻化体异常，且不能压制与电子导体不相关的激电异常，导致金属矿区圈定范围不准或者过大，使得后期找矿工作的不能顺利开展的问题。
本发明是这样实现的：一种激电中梯异常提取方法，包括以下步骤：
a、根据勘测区域地质图上的地层、岩性确定勘测区域的测点，通过激电测量得到每个测点的视极化率和视电阻率；
b、根据地层、岩性对勘测区域进行分区，并计算出每个分区内测点的视极化率平均值和视电阻率平均值/>
c、根据步骤b计算所得各地层、岩性分区的计算视极化率平均值和视电阻率平均值，结合每个测点的实测视电阻率，计算每个分区内的每个测点视极化率期望值ηq，ηq的计算公式为：
其中，ρi代表不同分区第i个测点实测视电阻率；
d、根据每个分区不同测点的视极化率期望值，结合每个测点的实测视极化率值，每个分区不同测点的衍生指数μ，μ的计算公式为：
其中，ηi代表不同分区第i个测点实测视极化率；
e、根据上步骤计算得到的每个分区不同测点的衍生指数绘制等值线平面图，在图中圈定出衍生指数异常区域，衍生指数异常区域即为激电中梯异常区域，即完成激电中梯异常提取工作。
步骤e中，等值线平面图通过绘图软件绘制得到。
本发明提供了一种激电中梯异常提取的新方法，该方法通过激电测量得到勘测区域测点的视极化率和视电阻率值，然后根据地层和岩性将勘测区域分成若干个分区，并计算出每个分区内测点的视极化率平均值和视电阻率平均值；再次计算每个分区内的每个测点视极化率期望值，并将每个测点的视极化率期望值与实测视极化率值结合，得到每个分区不同测点的衍生指数，最后根据得到的衍生指数绘制等值线平面图，等值线平面图中的衍生指数异常区域即为激电中梯异常区域。本发明可以突出与电子导体相关的低阻极化体异常，且能压制与电子导体不相关的激电异常，可以进一步缩小激电中梯异常范围，更加准确的圈定出最终异常地段，从而精确预测成矿的有利地段。
附图说明
图1是本发明的流程图。
图2是本发明实施例一中衍生指数等值线平面图。
图3是本发明实施例二中衍生指数等值线平面图。
具体实施方式
如图1所示，本发明提供了一种勘测区域激电中梯异常的提取方法，该方法包括以下步骤：
S1、首先根据勘测区域地质图上的地层、岩性确定勘测区域的测点，采用激电测量得到每个测点的视极化率和视电阻率；
S2、根据地层、岩性对勘测区域进行分区，并计算每个分区内测点的视极化率平均值和视电阻率平均值/>和/>的具体计算公式如下：
其中，ρi代表第i个测点实测视电阻率，ηi代表第i个测点实测视极化率，n代表同一分区测点个数；
S3、根据步骤S2计算所得各地层、岩性分区的计算视极化率平均值和视电阻率平均值，结合每个测点的实测视电阻率，计算每个分区内的每个测点视极化率期望值ηq，ηq的计算公式为：
S4、根据每个分区不同测点的视极化率期望值，结合每个测点的实测视极化率值，计算每个分区不同测点的衍生指数μ，μ的计算公式为：
其中，ηi代表不同分区第i个测点实测极化率；
S5、根据上步骤计算得到的每个分区不同测点的衍生指数绘制等值线平面图，在图中圈定出衍生指数异常区域，衍生指数异常区域即为激电中梯异常区域，即完成激电中梯异常提取工作。本步骤中的等值线平面图通过绘图软件绘制得到。
工作人员采用本方法在两个勘测区域进行了勘测实验，图1和图2分别是两个勘测区域内通过本方法得到的衍生指数等值线平面图，从图中可以看出见矿孔均在衍生指数异常圈定范围内，即采用本方法可以准确提取到的激电中梯异常。因此，通过提取勘测区衍生指数异常，结合激电中梯异常综合分析研究，能够为勘测区找矿提供依据和参考。
Claims (2)
1.一种激电中梯异常提取方法，其特征是，包括以下步骤：
a、根据勘测区域地质图上的地层、岩性确定勘测区域的测点，通过激电中梯测量得到每个测点的视极化率和视电阻率；
b、根据地层、岩性对勘测区域进行分区，并计算出每个分区内测点的视极化率平均值和视电阻率平均值/>
c、根据步骤b计算所得各地层、岩性分区的计算视极化率平均值和视电阻率平均值，结合每个测点的实测视电阻率，计算每个分区内的每个测点视极化率期望值ηq，ηq的计算公式为：
其中，ρi代表不同分区第i个测点实测视电阻率；
d、根据每个分区不同测点的视极化率期望值，结合每个测点的实测视极化率值，计算每个分区不同测点的衍生指数μ，μ的计算公式为：
其中，ηi代表不同分区第i个测点实测视极化率；
e、根据上步骤计算得到的每个分区不同测点的衍生指数绘制等值线平面图，在图中圈定出衍生指数异常区域，衍生指数异常区域即为激电中梯异常区域，即完成激电中梯异常提取工作。
2.根据权利要求1所述的激电中梯异常提取方法，其特征是，步骤e中，所述等值线平面图通过绘图软件绘制得到。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210661903.9A CN115079284B (zh) | 2022-06-13 | 2022-06-13 | 一种激电中梯异常提取方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210661903.9A CN115079284B (zh) | 2022-06-13 | 2022-06-13 | 一种激电中梯异常提取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115079284A CN115079284A (zh) | 2022-09-20 |
| CN115079284B true CN115079284B (zh) | 2024-06-21 |
Family
ID=83250786
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210661903.9A Active CN115079284B (zh) | 2022-06-13 | 2022-06-13 | 一种激电中梯异常提取方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115079284B (zh) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110118995A (zh) * | 2019-04-29 | 2019-08-13 | 山东省地质矿产勘查开发局第六地质大队 | 一种隐伏多金属矿找矿方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110133737B (zh) * | 2019-06-26 | 2020-03-06 | 中国科学院地质与地球物理研究所 | 一种隐伏矿的电磁学预测方法 |
| CN112965141B (zh) * | 2021-02-06 | 2024-03-08 | 核工业北京地质研究院 | 一种铀多金属矿的成矿有利地段的圈定方法 |
| CN113447992B (zh) * | 2021-06-22 | 2023-11-03 | 陕西地矿第二综合物探大队有限公司 | 一种运用时间域激发极化法进行矿体勘探的方法及系统 |
| CN114384593A (zh) * | 2022-01-18 | 2022-04-22 | 中国地质科学院 | 一种分布式三维激电数据采集处理装置及方法 |
-
2022
- 2022-06-13 CN CN202210661903.9A patent/CN115079284B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110118995A (zh) * | 2019-04-29 | 2019-08-13 | 山东省地质矿产勘查开发局第六地质大队 | 一种隐伏多金属矿找矿方法 |
Non-Patent Citations (1)
| Title |
|---|
| 激电中梯测量在云南岩脚铅锌矿区找矿中的应用;崔中良;洪托;刘洋;;煤矿机械;20160915(第09期);80-82 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115079284A (zh) | 2022-09-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108868805B (zh) | 基于统计分析与XGboost结合的盾构纠偏方法 | |
| CN104329079B (zh) | 一种识别气测录井油气层的方法以及系统 | |
| CN103233741B (zh) | 定量化评价超大型地下洞室群施工期围岩稳定性的方法 | |
| CN104636980B (zh) | 针对河道砂油藏类型油气汇集条件的地球物理表征方法 | |
| CN111458767B (zh) | 一种基于交汇图法识别岩性的方法及系统 | |
| CN107977480B (zh) | 一种页岩气储层产气性快速评价方法 | |
| CN111090709A (zh) | 一种砂岩型铀矿成矿预测的大数据地质分析方法 | |
| CN101236257A (zh) | 油井位置确定技术方法 | |
| CN107728227B (zh) | 一种井网成熟区快速判别地下河道砂体方法 | |
| CN106033127B (zh) | 基于横波速度变化率的地应力方位地震预测方法 | |
| CN115185015B (zh) | 一种深部锂铍矿勘查方法 | |
| CN107798189A (zh) | 一种基于透明空间地球物理的精准动态防突方法 | |
| CN111986321B (zh) | 三维地质建模方法、装置、设备及存储介质 | |
| CN112859183A (zh) | 一种城市地下管线的深度测量系统及方法 | |
| CN109324171A (zh) | 一种基于岩性统计的沉积相定量识别方法 | |
| CN103590828B (zh) | 一种录井dck指数法评价地层压力的方法 | |
| Ścigała | The identification of parameters of theories used for prognoses of post mining deformations by means of present software | |
| CN110968840A (zh) | 一种基于大地电磁测深电阻率判定隧围岩等级的方法 | |
| CN115079284B (zh) | 一种激电中梯异常提取方法 | |
| CN115542406A (zh) | 一种激电测量异常提取方法 | |
| CN104749618B (zh) | 泥页岩缓倾角裂缝叠后概率定量表征方法 | |
| CN105204066A (zh) | 一种基于谱分解的煤层火成岩侵入位置直接指示方法 | |
| CN115205479A (zh) | 一种多数据结合预测成矿远景区的方法 | |
| CN113804853A (zh) | 煤矿瓦斯含量检测方法及系统 | |
| CN111815769B (zh) | 逆冲推覆构造带构造的建模方法、计算设备及存储介质 |
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 |