CN114002746B - 基于可控源音频大地电磁测深法建立成矿模型的方法、装置和系统 - Google Patents
基于可控源音频大地电磁测深法建立成矿模型的方法、装置和系统 Download PDFInfo
- Publication number
- CN114002746B CN114002746B CN202111288926.1A CN202111288926A CN114002746B CN 114002746 B CN114002746 B CN 114002746B CN 202111288926 A CN202111288926 A CN 202111288926A CN 114002746 B CN114002746 B CN 114002746B
- Authority
- CN
- China
- Prior art keywords
- ore
- csamt
- inversion
- resistivity
- contour map
- 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
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
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (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
本申请提出一种基于可控源音频大地电磁测深法建立成矿模型的方法、装置和系统，所述方法包括：依据地形高程设置初始纯地形模型，对CSAMT观测数据进行反演，得到所述纯地形模型的电阻率值；对反演得到的所述电阻率值进行处理，形成CSAMT反演电阻率等值线图；将所述CSAMT反演电阻率等值线图和已知矿地质剖面进行对比分析，确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征；根据所述CSAMT反演电阻率等值线形态特征，判断生成热液型矿体的地质环境和接触带型矿与层间矿的结果，以此建立成矿模型。本发明正确判断可能生成热液型矿体的地质环境和接触带型矿与层间矿的方法，由此建立成矿模型。
Description
技术领域
本发明涉及地球物理勘探领域，具体涉及一种基于可控源音频大地电磁测深法建立成矿模型的方法、装置和系统。
背景技术
随着我国经济的快速发展，矿产资源需求量也随之增长，为满足日益增长的矿产需求，发现新矿产和在现有老矿山寻找盲矿和深部隐伏新矿及增加矿产资源储量是目前矿产资源勘查的工作重点。研究现有矿区的成矿模式对于找矿勘查具有理论指导意义，建立符合实际的成矿模式对于深部隐伏矿的发现有着现实意义。而以往通过地质工程推测成矿地质模型，往往需要大量的地质钻孔和其他的地质信息；也有通过地球物理技术找矿，由于受地球物理技术的限制不能从成矿机理分析，只能提供电性资料进行定性解释推断是否有矿的可能性，随着地球物理方法技术的迅速发展也极大地推动了成矿模式的研究，为深部矿的找寻提供科学证据。
可控源音频大地电磁法(CSAMT)是一种人工源频率域电磁测深方法，其特点是以人工场源来增加电磁信号强度，弥补天然场源信号微弱、不易观测等缺点，已被广泛地应用到矿产普查、油气勘探、水文环境等领域。通过近些年对仪器和数据处理的抗干扰技术研究，使得数据质量有了较大提升，在反演技术上研究了带地形带源的二维反演技术，以此克服了静态效应、地形影响和场源影响，使得反演电阻率断面逐渐接近实际地质结构。随着CSAMT观测的数据质量和反演技术的提升，提高了反演结果的可靠性，为研究成矿模型奠定了基础。
发明内容
本发明所要解决的技术问题是提供一种基于可控源音频大地电磁测深法建立成矿模型的方法、装置和系统，以便建立成矿模型，为找矿提供依据。
为了实现上述发明目的，本发明采取的技术方案如下：
本发明提供一种基于可控源音频大地电磁测深法建立成矿模型的方法，包括：
依据地形高程设置初始纯地形模型，对CSAMT观测数据进行反演，得到所述纯地形模型的电阻率值；
对反演得到的所述电阻率值进行处理，形成CSAMT反演电阻率等值线图；
在有已知矿地质资料的剖面上将所述CSAMT反演电阻率等值线图和已知矿地质剖面进行对比分析，确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征；
根据所述CSAMT反演电阻率等值线形态特征，判断生成热液型矿体的地质环境和接触带型矿与层间矿的结果；
在未知剖面上利用磁法判断深部岩体的存在后并用CSAMT反演电阻率等值线形态特征进行生成热液型矿体的地质环境和接触带型矿与层间矿的判断，并建立成矿模型。
可选地，依据地形高程设置初始纯地形模型，对CSAMT观测数据进行反演，得到所述纯地形模型的电阻率值包括：
获取带地形剖分的初始模型，设定所述初始模型的初始电阻率值，经过反演和迭代，修正所述初始模型的电阻率值，直到达到预设条件后停止反演，得到纯地形模型每个区域的电阻率值，并得到每个点不同深度的反演电阻率。
可选地，纯地形模型为带地形的地电模型。
可选地，对反演得到的所述电阻率值进行处理包括以下至少之一：
对反演得到的所述电阻率值进行网格化、加地形或加色标。
可选地，确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征包括：
针对CSAMT反演电阻率等值线图，确定与热液型矿体形成对应的深部断裂和形成接触带型矿的花岗岩形态，以及层间矿在电阻率等值线图上的表现特征。
可选地，根据所述CSAMT反演电阻率等值线形态特征，判断生成热液型矿体的地质环境和接触带型矿与层间矿的结果包括：根据所述CSAMT反演电阻率等值线图判断深部的断裂与花岗岩岩体连接情况，确定深部存在热液上升通道的断裂带，所述CSAMT反演电阻率等值线图对断裂带反应为低阻带，对花岗岩岩体反应为高阻，对接触带型矿的花岗岩表现为在高阻隆起形态的局部凹陷或隆起旁侧低处存在凸起的特征形态，这种高阻隆起形态的局部凹陷或隆起旁侧低处存在凸起有利于含矿热液顺着断裂上升的途径中在此形态部位富集成矿；对层间矿的反应为在高阻隆起旁侧的呈现水平形态特征，这种形态特征有利于含矿热液沿微裂隙充填，形成层间氧化矿。
本发明还提供一种基于可控源音频大地电磁测深法建立成矿模型的装置，包括：
反演模块，依据地形高程设置初始纯地形模型，对CSAMT观测数据进行反演，得到所述纯地形模型的电阻率值；
成像模块，设置为对反演得到的所述电阻率值进行处理，形成CSAMT反演电阻率等值线图；
比对模块，设置为在有已知矿地质资料的剖面上将所述CSAMT反演电阻率等值线图和已知矿地质剖面进行对比分析，确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征；
生成模块，设置为根据所述CSAMT反演电阻率等值线形态特征，判断生成热液型矿体的地质环境和接触带型矿与层间矿的结果；
建模模块，设置为在未知剖面上利用磁法判断深部岩体的存在后并用CSAMT反演电阻率等值线形态特征进行生成热液型矿体的地质环境和接触带型矿与层间矿的判断，并建立成矿模型。
可选地，所述比对模块确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征包括：
针对CSAMT反演电阻率等值线图，确定与热液型矿体形成对应的深部断裂和形成接触带型矿的花岗岩形态，以及层间矿在电阻率等值线图上的表现特征。
本发明还提供一种基于可控源音频大地电磁测深法建立成矿模型的系统，包括：可控源音频大地电磁测深器、CSAMT反演器和中央处理器；
所述可控源音频大地电磁测深器包括：发射系统和接收系统，其中，所述发射系统通过一对电偶电极A极、B极向地下发射发送不同频率的交变电流，所述接收系统通过两个不极化电极M极、N极接收反馈的水平电场信号Ex，以及与接收电场相交的水平磁场信号Hy；由观测的电场信号和磁场信号计算出视电阻率和相位；
所述CSAMT反演器设置为：依据地形高程设置初始纯地形模型，对CSAMT观测数据进行反演，得到所述纯地形模型的电阻率值；
所述中央处理器设置为：对反演得到的所述电阻率值进行处理，形成CSAMT反演电阻率等值线图；在有已知矿地质资料的剖面上将所述CSAMT反演电阻率等值线图和已知矿地质剖面进行对比分析，确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征；根据所述CSAMT反演电阻率等值线形态特征，判断生成热液型矿体的地质环境和接触带型矿与层间矿的结果；在未知剖面上利用磁法判断深部岩体的存在后并用CSAMT反演电阻率等值线形态特征进行生成热液型矿体的地质环境和接触带型矿与层间矿的判断，并建立成矿模型。
可选地，所述中央处理器确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征包括：
针对CSAMT反演电阻率等值线图，确定与热液型矿体形成对应的深部断裂和形成接触带型矿的花岗岩形态，以及层间矿在电阻率等值线图上的表现特征。
本发明和现有技术相比，具有如下有益效果：
本发明利用CSAMT反演的电阻率形态建立成矿模型，其成本较以往地质工程要经济很多。
本发明依据地形高程设置初始纯地形模型，对CSAMT观测数据进行反演，得到模型的电阻率值，利用绘制等值线图软件对反演的电阻率值进行网格化、加地形和加色标等手段，形成反演电阻率等值线图，将所述CSAMT反演电阻率等值线图和已知矿剖面进行对比分析，确定不同地质成矿模型在CSAMT反演电阻率等值线形态特征：接触带型矿在反演电阻率等值线形态表现为高阻异常顶部有局部凹陷形态或高阻隆起旁侧的低处存在凸起形态；层间矿在深部隆花岗岩岩体旁侧的相对低阻等值线为水平形态。
本发明通过反演电阻率形态与已知地质剖面对比，给出了以CSAMT反演电阻率等值线形态判断是否有矿体存在，生成接触带型矿和层间矿的电阻率表现特征，以便正确判断可能生成热液型矿体的地质环境和接触带型矿与层间矿的方法，由此建立成矿模型；在未知剖面上可用磁法判断深部岩体的存在后并直接用CSAMT反演电阻率等值线形态特征进行生成热液型矿体的地质环境和接触带型矿与层间矿的判断，并建立成矿模型。
本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
图1为本发明实施例的基于CSAMT测深建立成矿模型的方法的流程图；
图2为本发明实施例的基于CSAMT测深建立成矿模型的装置的结构示意图；
图3为本发明实施例的可控源音频大地电磁测深器的示意图；
图4是本申请实施例的CSAMT反演电阻率断面图；
图5是本申请实施例的已知地质资料和地质剖面；
图6是本申请实施例的地质成矿模型图；
图7是本申请实施例建立的成矿模型图。
具体实施方式
为使本发明的发明目的、技术方案和有益效果更加清楚明了，下面结合附图对本发明的实施例进行说明，需要说明的是，在不冲突的情况下，本申请中的实施例和实施例中的特征可以相互任意组合。
在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行。并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
如图1所示，本发明实施例提供一种基于可控源音频大地电磁测深建立成矿模型的方法，可以包括如下步骤S101至步骤S104：
S101、依据地形高程设置初始纯地形模型，对CSAMT观测数据进行反演，得到所述纯地形模型的电阻率值；
S102、对反演得到的所述电阻率值进行处理，形成CSAMT反演电阻率等值线图；
S103、在有已知矿地质资料的剖面上将所述CSAMT反演电阻率等值线图和已知矿地质剖面进行对比分析，确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征；
S104、根据所述CSAMT反演电阻率等值线形态特征，判断生成热液型矿体的地质环境和接触带型矿与层间矿的结果；
S105、在未知剖面上利用磁法判断深部岩体的存在后并用CSAMT反演电阻率等值线形态特征进行生成热液型矿体的地质环境和接触带型矿与层间矿的判断，并建立成矿模型。
本发明实施例，在未知剖面上可用磁法判断深部岩体的存在后并用CSAMT反演电阻率等值线形态特征判断形成矿体的地质环境和接触带型矿与层间矿；根据所述CSAMT反演电阻率等值线形态特征，判断生成热液型矿体的地质环境和接触带型矿与层间矿的结果，以此建立成矿模型；在未知剖面上可用磁法判断深部岩体的存在后并用CSAMT反演电阻率等值线形态特征进行生成热液型矿体的地质环境和接触带型矿与层间矿的判断，并建立成矿模型。
本发明实施例中，步骤S101依据地形高程设置初始纯地形模型，对CSAMT观测数据进行反演，得到所述纯地形模型的电阻率值包括：
获取带地形剖分的初始模型，设定所述初始模型的初始电阻率值，经过反演和迭代，修正所述初始模型的电阻率值，直到达到预设条件后停止反演，得到纯地形模型每个区域的电阻率值，并得到每个点不同深度的反演电阻率。本发明实施例中的“剖分”，将地形模型中大的区域划分为多个矩形或三角形，可以是均等的也可以是非均等大小的。把模型剖分成很多网格，每个网格可以设定统一的电阻率值，例如给每个网格初始值为100Ω·m。
本发明实施例中纯地形模型可以为带地形的地电模型。
本发明实施例中和已知矿地质剖面进行对比的目的是为了说明CSAMT反演电阻率等值线形态和成矿空间的关系是正确的，由此可以从CSAMT反演电阻率等值线形态在没有地质剖面上也可以进行成矿地质环境和成矿的判断和成矿模型的建立。
本发明实施例中，步骤S102对反演得到的所述电阻率值进行处理包括以下至少之一：
对反演得到的所述电阻率值进行网格化、加地形或加色标。
本发明实施例中，步骤S103中确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征包括：
针对CSAMT反演电阻率等值线图，确定与热液型矿体形成对应的深部断裂和形成接触带型矿的花岗岩形态，以及层间矿在电阻率等值线图上的表现特征。
本发明实施例中，步骤S104中根据所述CSAMT反演电阻率等值线形态特征，判断生成热液型矿体的地质环境和接触带型矿与层间矿的结果包括：
根据所述CSAMT反演电阻率等值线图判断深部的断裂与花岗岩岩体连接情况，确定深部存在热液上升通道的断裂带，所述CSAMT反演电阻率等值线图对断裂带反应为低阻带，对花岗岩岩体反应为高阻，对接触带型矿的花岗岩表现为在高阻隆起形态的局部凹陷或隆起旁侧低处存在凸起的特征形态，这种高阻隆起形态的局部凹陷或隆起旁侧低处存在凸起有利于含矿热液顺着断裂上升的途径中在此形态部位富集成矿；对层间矿的反应为在高阻隆起旁侧的呈现水平形态特征，这种形态特征有利于含矿热液沿微裂隙充填，形成层间氧化矿。
如图2所示，本发明实施例基于可控源音频大地电磁测深法建立成矿模型的装置，包括：
反演模块，依据地形高程设置初始纯地形模型，对CSAMT观测数据进行反演，得到所述纯地形模型的电阻率值；
成像模块，设置为对反演得到的所述电阻率值进行处理，形成CSAMT反演电阻率等值线图；
比对模块，设置为在有已知矿地质资料的剖面上将所述CSAMT反演电阻率等值线图和已知矿地质剖面进行对比分析，确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征；
生成模块，设置为根据所述CSAMT反演电阻率等值线形态特征，判断生成热液型矿体的地质环境和接触带型矿与层间矿的结果；
建模模块，设置为在未知剖面上利用磁法判断深部岩体的存在后并用CSAMT反演电阻率等值线形态特征进行生成热液型矿体的地质环境和接触带型矿与层间矿的判断，并建立成矿模型。
本发明实施例中，如果在未知剖面进行建模时，可以省略比对模块。
本发明实施例中，所述比对模块确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征包括：
针对CSAMT反演电阻率等值线图，确定与热液型矿体形成对应的深部断裂和形成接触带型矿的花岗岩形态，以及层间矿在电阻率等值线图上的表现特征。
本发明实施例的基于可控源音频大地电磁测深法建立成矿模型的系统，包括：可控源音频大地电磁测深器、CSAMT反演器和中央处理器；
如图3所示，所述可控源音频大地电磁测深器包括：发射系统和接收系统，其中，所述发射系统通过一对电偶电极A极、B极向地下发射发送不同频率的交变电流，所述接收系统通过两个不极化电极M极、N极接收反馈的水平电场信号Ex，以及与接收电场相交的水平磁场信号Hy；由观测的电场信号和磁场信号计算出视电阻率和相位；
所述CSAMT反演器设置为：依据地形高程设置初始纯地形模型，对CSAMT观测数据进行反演，得到所述纯地形模型的电阻率值；
所述中央处理器设置为：对反演得到的所述电阻率值进行处理，形成CSAMT反演电阻率等值线图；在有已知矿地质资料的剖面上将所述CSAMT反演电阻率等值线图和已知矿地质剖面进行对比分析，确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征；根据所述CSAMT反演电阻率等值线形态特征，判断生成热液型矿体的地质环境和接触带型矿与层间矿的结果；在未知剖面上利用磁法判断深部岩体的存在后并用CSAMT反演电阻率等值线形态特征进行生成热液型矿体的地质环境和接触带型矿与层间矿的判断，并建立成矿模型。
本发明实施例中，所述中央处理器确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征包括：
针对CSAMT反演电阻率等值线图，确定与热液型矿体形成对应的深部断裂和形成接触带型矿的花岗岩形态，以及层间矿在电阻率等值线图上的表现特征。
实施例一
由于CSAMT法在野外施工时，受发射机功率的限制，无法将发射源放在距接收点很远的地方，一般而言收发距是探测深度的4-6倍，目前实际工作中对于30kw的发射机的最大收发距为8-15km，这样采集到的数据为全区数据(包含远区、过渡区和近区)。利用带地形带源反演软件对CSAMT观测数据进行反演，最终得到模型的电阻率值；利用等值线图软件对反演的电阻率值绘制反演电阻率等值线图；将所述CSAMT反演电阻率等值线图和已知地质剖面对比；确定不同地质成矿模型在CSAMT反演电阻率等值线形态特征；再以CSAMT反演电阻率等值线形态特征，结合成矿规律建立成矿模型，确定成矿类型。
图4为云南某锡、铜多金属矿区D9线的CSAMT数据经过带起伏地形带源二维反演电阻率断面图。矿区内出露的地层以中生界三叠系碳酸盐为主，其中三叠系中统个旧组(T2g)为本区的主要熔矿层位。由灰岩、白云质灰岩、灰质白云岩和白云岩互层，局部夹薄层泥质灰岩组成。总厚度超过1000m，尤其是个旧组下段(T2g1 1-6)是整个个旧矿区熔矿最为丰富的层位。矿体与围岩具有明显的物性差异，无金属矿化的大理岩、白云岩、花岗岩电阻率较高(274-5841Ω·m)、极化率较低，而含金属矿化的各类岩石电阻率(299-1318Ω·m)、极化率最高，硫化矿为低阻(0.7-7.5Ω·m)高极化率，氧化矿则以高阻(639-1444Ω·m)、低极化率为特征。
图4反演电阻率断面反映出低阻-高阻-低阻-高阻的四个电性层，是断面内三叠系中统个旧组下段(T2g1 1-4)地层的清晰反映。三叠系地层整体上倾向西，倾角较缓。浅部电性层较为凌乱，有多处低阻与高阻相间异常(750点～890点、1340点～1460点、1790点～1890点)，和已知F1、F2、F3断裂相对应，其位置及产状与地质资料相符(图5)。高程1400m左右的高阻层形态与花岗岩隆起形态吻合，此层之上的低阻层为个旧组下段T2g1 1地层，是矿区熔矿最为丰富的层位。在1090点下高程1450m处和1640点下高程1300m处存在明显的局部低阻异常DY1和DY2，为矿体的反映，与矿体位置及形态相符，局部低阻异常DY1和DY2显示出较好的矿致异常信息，而左侧出现的局部低阻异常DY3无矿体显示。
结合地质断面和地质资料，在认识电阻率等值线形态特征的基础上，建立成矿模型：锡矿的成矿热液包括成矿物质主要来源于花岗岩浆演化热液，花岗岩浆沿着断裂裂隙上升，含有矿物元素的热液在花岗岩顶部的凹陷部位或受花岗岩凸起的阻挡汇集在一起，从而形成接触带硫化物型锡矿，而在隆起花岗岩岩体右侧的含矿热液沿白云岩中脆弱的微裂隙充填的形成层间氧化矿型，而在左侧的低阻异常形如斜坡，无任何褶皱或凸起阻挡，从而不具备矿物元素聚集的环境，无法形成矿体。
建立起该剖面的成矿模型如图7所示，在深部的成矿热液沿着F1、F2断裂裂隙上升，在花岗岩顶部的凹陷部位形成接触带硫化物型锡矿，在隆起花岗岩岩体右侧因含矿热液沿白云岩中脆弱的微裂隙充填的形成层间氧化矿，地表浅部在断裂附近形成的有砂矿。以上的成矿模型推断与汪志芬总结的个旧矿区成矿模式图(图6)基本一致，后期探矿过程中依据该成矿模式在岩体右侧1600m至1900m地段见3层层间氧化锡矿。图6中，(1)砂锡；
(2)细脉带型矿床；(3)含锂云母萤石脉；(4)电气石脉；(5)绿柱石黑电气石脉；(6)绿柱石长石脉；(7)矽卡岩脉；(8)正接触带矽卡岩；(9)接触带矽卡岩硫化物型矿床；(10)锡石-硫化物型矿床；(11)东西向锡铅矿带；(12)含锡白云岩型矿床；(13)变辉绿岩型铜矿床；(14)顶部岩浆熔离带；(15)亲氧富挥发分熔浆；(16)重金属硫化物熔浆；(17)晚期富集硫化物熔浆；(18)燕山中晚期黑云母花岗岩。
由CSAMT反演电阻率等值线形态建立成矿模型的重点要素有：深部的断裂与花岗岩岩体连接情况，深部花岗岩岩体形态在其隆起部是否局部凹陷或旁侧存在凸起，深部花岗岩岩体旁侧形态是否是水平形态。
本发明实施例与通常通过大量的钻探钻孔资料进行的地质成矿建模会付出高额的经济成本，通常利用地球物理的电磁法依靠电阻率低值异常判断矿体的可能性，由此脱离成矿环境形态分析可导致在模型建立时形成错误的推断结果，利用本方法可以正确建立成矿模型。可以依据本文中所公开的实施方法较容易的进行成矿模型建立，后期再结合地质资料形成丰富的地质成矿模型，更好地为地质找矿服务。
虽然本发明所揭示的实施方式如上，但其内容只是为了便于理解本发明的技术方案而采用的实施方式，并非用于限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所揭示的核心技术方案的前提下，可以在实施的形式和细节上做任何修改与变化，但本发明所限定的保护范围，仍须以所附的权利要求书限定的范围为准。
Claims (10)
1.基于可控源音频大地电磁测深法建立成矿模型的方法，其特征在于，包括：
依据地形高程设置初始纯地形模型，对CSAMT观测数据进行反演，得到所述纯地形模型的电阻率值；
对反演得到的所述电阻率值进行处理，形成CSAMT反演电阻率等值线图；
在有已知矿地质资料的剖面上将所述CSAMT反演电阻率等值线图和已知矿地质剖面进行对比分析，确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征；
根据所述CSAMT反演电阻率等值线图中的形态特征，判断生成热液型矿体的地质环境和接触带型矿与层间矿的结果；
在未知剖面上利用磁法判断深部岩体的存在后并用CSAMT反演电阻率等值线图中的形态特征进行生成热液型矿体的地质环境和接触带型矿与层间矿的判断，并建立成矿模型。
2.如权利要求1所述的方法，其特征在于：依据地形高程设置初始纯地形模型，对CSAMT观测数据进行反演，得到所述纯地形模型的电阻率值包括：
获取带地形剖分的初始模型，设定所述初始模型的初始电阻率值，经过反演和迭代，修正所述初始模型的电阻率值，直到达到预设条件后停止反演，得到纯地形模型每个区域的电阻率值，并得到每个点不同深度的反演电阻率。
3.如权利要求1所述的方法，其特征在于：纯地形模型为带地形的地电模型。
4.如权利要求1所述的方法，其特征在于：对反演得到的所述电阻率值进行处理包括以下至少之一：
对反演得到的所述电阻率值进行网格化、加地形或加色标。
5.如权利要求1所述的方法，其特征在于：确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征包括：
针对CSAMT反演电阻率等值线图，确定与热液型矿体形成对应的深部断裂和形成接触带型矿的花岗岩形态，以及层间矿在电阻率等值线图上的表现特征。
6.如权利要求5所述的方法，其特征在于，根据所述CSAMT反演电阻率等值线图中的形态特征，判断生成热液型矿体的地质环境和接触带型矿与层间矿的结果包括：
根据所述CSAMT反演电阻率等值线图判断深部的断裂与花岗岩岩体连接情况，确定深部存在热液上升通道的断裂带，所述CSAMT反演电阻率等值线图对断裂带反应为低阻带，对花岗岩岩体反应为高阻，对接触带型矿的花岗岩表现为在高阻隆起形态的局部凹陷或隆起旁侧低处存在凸起的特征形态，这种高阻隆起形态的局部凹陷或隆起旁侧低处存在凸起有利于含矿热液顺着断裂上升的途径中在此形态部位富集成矿；对层间矿的反应为在高阻隆起旁侧的呈现水平形态特征，这种形态特征有利于含矿热液沿微裂隙充填，形成层间氧化矿。
7.基于可控源音频大地电磁测深法建立成矿模型的装置，其特征在于，包括：
反演模块，依据地形高程设置初始纯地形模型，对CSAMT观测数据进行反演，得到所述纯地形模型的电阻率值；
成像模块，设置为对反演得到的所述电阻率值进行处理，形成CSAMT反演电阻率等值线图；
比对模块，设置为在有已知矿地质资料的剖面上将所述CSAMT反演电阻率等值线图和已知矿地质剖面进行对比分析，确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征；
生成模块，设置为根据所述CSAMT反演电阻率等值线图中的形态特征，判断生成热液型矿体的地质环境和接触带型矿与层间矿的结果；
建模模块，设置为在未知剖面上利用磁法判断深部岩体的存在后并用CSAMT反演电阻率等值线图中的形态特征进行生成热液型矿体的地质环境和接触带型矿与层间矿的判断，并建立成矿模型。
8.如权利要求7所述的装置，其特征在于，所述比对模块确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征包括：
针对CSAMT反演电阻率等值线图，确定与热液型矿体形成对应的深部断裂和形成接触带型矿的花岗岩形态，以及层间矿在电阻率等值线图上的表现特征。
9.基于可控源音频大地电磁测深法建立成矿模型的系统，其特征在于，包括：可控源音频大地电磁测深器、CSAMT反演器和中央处理器；
所述可控源音频大地电磁测深器包括：发射系统和接收系统，其中，所述发射系统通过一对电偶电极A极、B极向地下发射发送不同频率的交变电流，所述接收系统通过两个不极化电极M极、N极接收反馈的水平电场信号Ex，以及与接收电场相交的水平磁场信号Hy；由观测的电场信号和磁场信号计算出视电阻率和相位；
所述CSAMT反演器设置为：依据地形高程设置初始纯地形模型，对CSAMT观测数据进行反演，得到所述纯地形模型的电阻率值；
所述中央处理器设置为：对反演得到的所述电阻率值进行处理，形成CSAMT反演电阻率等值线图；在有已知矿地质资料的剖面上将所述CSAMT反演电阻率等值线图和已知矿地质剖面进行对比分析，确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征；根据所述CSAMT反演电阻率等值线图中的形态特征，判断生成热液型矿体的地质环境和接触带型矿与层间矿的结果；在未知剖面上利用磁法判断深部岩体的存在后并用CSAMT反演电阻率等值线图中的形态特征进行生成热液型矿体的地质环境和接触带型矿与层间矿的判断，并建立成矿模型。
10.根据权利要求9所述的系统，其特征在于，所述中央处理器确定接触带型矿和层间矿在所述CSAMT反演电阻率等值线图中的形态特征包括：
针对CSAMT反演电阻率等值线图，确定与热液型矿体形成对应的深部断裂和形成接触带型矿的花岗岩形态，以及层间矿在电阻率等值线图上的表现特征。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111288926.1A CN114002746B (zh) | 2021-11-02 | 2021-11-02 | 基于可控源音频大地电磁测深法建立成矿模型的方法、装置和系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111288926.1A CN114002746B (zh) | 2021-11-02 | 2021-11-02 | 基于可控源音频大地电磁测深法建立成矿模型的方法、装置和系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114002746A CN114002746A (zh) | 2022-02-01 |
| CN114002746B true CN114002746B (zh) | 2022-06-17 |
Family
ID=79926475
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111288926.1A Active CN114002746B (zh) | 2021-11-02 | 2021-11-02 | 基于可控源音频大地电磁测深法建立成矿模型的方法、装置和系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114002746B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104597519A (zh) * | 2015-02-04 | 2015-05-06 | 山东省物化探勘查院 | 一种深部金矿阶梯式找矿方法 |
| CN107037492A (zh) * | 2017-05-26 | 2017-08-11 | 贵州省地质矿产勘查开发局0三地质大队 | 一种地质数据分析建模方法 |
| CN108020865A (zh) * | 2016-11-02 | 2018-05-11 | 核工业北京地质研究院 | 一种花岗岩型铀矿深部有利成矿空间识别及定位方法 |
| CN113325486A (zh) * | 2021-06-03 | 2021-08-31 | 青海省第一地质勘查院 | 覆盖区下的构造蚀变岩型矿物勘查方法、系统及装置 |
| CN113391358A (zh) * | 2021-06-03 | 2021-09-14 | 中国科学院地质与地球物理研究所 | 一种宽频大地电磁数据的采集方法和采集装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8014988B2 (en) * | 2006-04-06 | 2011-09-06 | Exxonmobil Upstream Research Co. | Method for obtaining resistivity from controlled source electromagnetic data |
-
2021
- 2021-11-02 CN CN202111288926.1A patent/CN114002746B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104597519A (zh) * | 2015-02-04 | 2015-05-06 | 山东省物化探勘查院 | 一种深部金矿阶梯式找矿方法 |
| CN108020865A (zh) * | 2016-11-02 | 2018-05-11 | 核工业北京地质研究院 | 一种花岗岩型铀矿深部有利成矿空间识别及定位方法 |
| CN107037492A (zh) * | 2017-05-26 | 2017-08-11 | 贵州省地质矿产勘查开发局0三地质大队 | 一种地质数据分析建模方法 |
| CN113325486A (zh) * | 2021-06-03 | 2021-08-31 | 青海省第一地质勘查院 | 覆盖区下的构造蚀变岩型矿物勘查方法、系统及装置 |
| CN113391358A (zh) * | 2021-06-03 | 2021-09-14 | 中国科学院地质与地球物理研究所 | 一种宽频大地电磁数据的采集方法和采集装置 |
Non-Patent Citations (3)
| Title |
|---|
| CSAMT在强干扰矿区的找矿效果;雷达;《首届全国矿产勘查大会论文集》;20211012;第253-257页 * |
| 坦桑尼亚环维多利亚湖地区金矿勘查地球物理方法应用研究;白德胜;《地球物理学进展》;20201109;第36卷(第3期);第1029-1045页 * |
| 藏南扎西康矿集区地下结构与成矿模式探讨;梁斯琪;《中国优秀硕士学位论文全文数据库 工程科技I辑》;20180215(第02期);第3-4章 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114002746A (zh) | 2022-02-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Mao et al. | Three-dimensional prospectivity modeling of the Jiaojia-type gold deposit, Jiaodong Peninsula, Eastern China: A case study of the Dayingezhuang deposit | |
| CN106569287B (zh) | 一种基于复合砂体构型模式的扇三角洲前缘储层预测方法 | |
| CN109738947B (zh) | 一种圈定砂岩型铀矿找矿远景区的物化探组合方法 | |
| CN113109889B (zh) | 一种基于“双阶段双模式”成矿模型的砂岩型铀矿找矿方法 | |
| Malinowska et al. | The impact of deep underground coal mining on earth fissure occurrence | |
| CN106503834A (zh) | 一种湖相超低孔渗储层致密油甜点区的预测方法 | |
| CN108802851A (zh) | 一种沉积变质型铁-石墨矿勘探方法 | |
| CN111967631A (zh) | 一种砂岩型铀矿勘探早期成矿远景区预测的方法 | |
| CN105988143A (zh) | 一种适用于矿田深部铀资源潜力定位预测的方法 | |
| CN114386464B (zh) | 一种瞬变电磁激电信息的深度学习提取方法 | |
| Perring et al. | A new fluid-flow model for the genesis of banded iron formation-hosted martite-goethite mineralization, with special reference to the North and South Flank deposits of the Hamersley Province, Western Australia | |
| CN110489809A (zh) | 一种盆地油气资源整体评价方法及装置 | |
| CN114002746B (zh) | 基于可控源音频大地电磁测深法建立成矿模型的方法、装置和系统 | |
| CN107479100B (zh) | 一种陡倾斜深部隐伏铁矿找矿方法 | |
| CN111983719A (zh) | 一种适用于煤田勘查区的砂岩型铀矿快速评价方法 | |
| CN110118992B (zh) | 一种全隐蔽深部煤田煤炭资源勘查方法 | |
| Niiranen et al. | Central Lapland Greenstone Belt 3D modeling project final report | |
| Spencer | Magnetite deposits of the Cornwall type in Pennsylvania | |
| CN115861551A (zh) | 一种地浸采铀的数字建井方法 | |
| CN115238966A (zh) | 一种预测油-铀叠置区砂岩型铀矿有利区带的方法 | |
| Ashton et al. | Discovery of the Tara Deep Zn-Pb Mineralization at the Boliden Tara Mine, Navan, Ireland: Success with Modern Seismic Surveys. | |
| Leney | Field studies in iron ore geophysics | |
| CN109188518B (zh) | 基于地震分频技术的煤系地层砂岩识别方法和系统 | |
| Yatini et al. | Application of polarization method (IP) to delineate the gold mineralization zones in Cihonje Areas, Banyumas Regency, Province of Central Java | |
| Minglei et al. | The Enrichment Conditions and Resource Potential of Marine-Continental Transitional Coal-Measure Shale Gas: A Case Study of the Permo-Carboniferous Systerm in the Huanghebei Coalfield of North China |
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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Lei Da Inventor after: Di Qingyun Inventor after: Fu Changmin Inventor after: Wang Ruo Inventor after: An Zhiguo Inventor after: Zhen Qihui Inventor after: Ren Hao Inventor before: Lei Da Inventor before: Di Qingyun Inventor before: Fu Changmin Inventor before: Wang Ruo Inventor before: An Zhiguo Inventor before: Zhen Qihui Inventor before: Ren Hao |