E-mail Alert Rss
 

物探与化探, 2022, 46(3): 618-627 doi: 10.11720/wtyht.2022.1189

地质调查·资源勘查

三维AMT正反演技术对贵州马坪含金刚石岩体探测的精细解释

何帅,1,2, 杨炳南,1,2,3, 阮帅4, 李永刚5, 韩姚飞1, 朱大伟1

1.贵州省地质矿产勘查开发局 103地质大队,贵州 铜仁 554300

2.自然资源部 基岩区矿产资源勘查工程技术创新中心,贵州 贵阳 550001

3.中国地质大学(武汉) 地球物理与空间信息学院,湖北 武汉 430074

4.中国地质科学院 深部探测中心,北京 100037

5.贵州省地质矿产勘查开发局101地质大队,贵州 凯里 556000

Fine Interpretation of the exploration results of diamond-bearing rock masses in Maping area, Guizhou using the 3D AMT forward modeling and inversion technologies

HE Shuai,1,2, YANG Bing-Nan,1,2,3, RUAN Shuai4, LI Yong-Gang5, HAN Yao-Fei1, ZHU Da-Wei1

1. No. 103 Geologic Team, Bureau of Geology and Mineral Exploration and Development of Guizhou, Tongren 554300, China

2. Engineering Technology Innovation Center of Resources Explorations in Basement Area of China, Ministry of Natural Resources, Guiyang 550001,China

3. Institute of Geophysics and Geomatics, China University of Geosciences, Wuhan 430074, China

4. Sinoprobe Center-China Deep Exploration Center, Chinese Academy of Geological Sciences, Beijing 100037, China

5. No. 101 Geologic Team, Bureau of Geology and Mineral Exploration and Development of Guizhou, Kaili 556000, China

通讯作者: 杨炳南(1985-),男,在读博士,高级工程师,主要从事地球物理勘查与应用工作。Email:648417001@qq.com

责任编辑: 沈效群

收稿日期: 2021-04-8   修回日期: 2021-11-7  

基金资助: 贵州省科技计划项目(黔科合支撑[2019]2868号)
贵州省地质矿产勘查开发局地质科学研究项目(黔地矿科合[2016]07号)
中国地质调查局项目“华北和扬子地区金刚石矿产调查”(DD20160059)

Received: 2021-04-8   Revised: 2021-11-7  

作者简介 About authors

何帅(1987-),男,高级工程师,主要从事地球物理勘查与应用工作。Email: 307050903@qq.com

摘要

贵州镇远马坪“东方一号”岩体是我国首次发现的含金刚石原生矿母岩,马坪地区发现的岩体属于金伯利岩浆体系的浅部相,其深部可能存在规模较大的隐伏岩管或岩筒。为揭示马坪地区深部含金刚石隐伏岩管或岩筒的空间展布特征,在区内开展了80 m×40 m高密集网度的音频大地电磁勘探工作;利用三维正演技术模拟研究区纯地形响应并在实测数据中去除,得到的定性解释结果在一定程度上恢复了被静态效应扭曲的AMT阻抗相位不变量分布形态;使用AR-QN拟牛顿反演方法对数据进行三维反演,根据研究区岩性统计结果设定地下单元的电阻率变化区间,获得了可靠的三维电性结构;最后依据地表发现的岩筒、钻孔揭露的多条岩脉等地质资料对该地电模型进行精细解释,勾画出了隐伏岩管(岩筒)的形态,为区内下一步金刚石原生矿找矿方向及预测提供了地球物理依据。

关键词： 金刚石隐伏岩筒; 音频大地电磁法; 三维正反演; 阻抗相位不变量; 贵州马坪

Abstract

The "Dongfang No.1" rock mass in the Maping area, Zhenyuan County, Guizhou is the parent rock of a primary diamond-bearing deposit discovered in China for the first time. Studies have shown that the rock mass found in the Maping area is of the shallow facies of the kimberlite magmatic system, and large-scale concealed rock pipes or buckets may exist in the deep part. To reveal the spatial distribution of deep diamond-bearing concealed rock pipes or buckets in the Maping area, this study carried out the audio-magnetotelluric (AMT) data acquisition in the area using a high grid density of 80 m × 40 m. Then it simulated the pure terrain response in the study area using the 3D forward modeling and deducted the pure terrain response from the measured data. The obtained qualitative interpretation results restored the distribution morphology of AMT impedance phase invariants to some extent that was distorted by static effects. Afterward, this study performed the 3D inversion of the data using the AR-QN quasi-Newtonian inversion method. Based on the lithologic statistical results of the study area, the resistivity variation intervals of the underground units were set during the inversion, obtaining a reliable 3D electrical structure. Finally, this study carried out a fine interpretation of the geoelectric model of this area based on geological data such as rock tubes found on the surface and multiple dikes revealed by boreholes, outlining the morphology of concealed rock pipes or buckets. This study will provide a geophysical basis for the future prospecting and prediction of primary diamond deposits in this area.

Keywords： concealed diamond tube; audio-magnetotelluric (AMT); 3D forward modeling and inversion; impedance phase invariant; Maping area of Guizhou

PDF (7036KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

何帅, 杨炳南, 阮帅, 李永刚, 韩姚飞, 朱大伟. 三维AMT正反演技术对贵州马坪含金刚石岩体探测的精细解释[J]. 物探与化探, 2022, 46(3): 618-627 doi:10.11720/wtyht.2022.1189

HE Shuai, YANG Bing-Nan, RUAN Shuai, LI Yong-Gang, HAN Yao-Fei, ZHU Da-Wei. Fine Interpretation of the exploration results of diamond-bearing rock masses in Maping area, Guizhou using the 3D AMT forward modeling and inversion technologies[J]. Geophysical and Geochemical Exploration, 2022, 46(3): 618-627 doi:10.11720/wtyht.2022.1189

0 引言

1965年,在贵州镇远马坪地区首次发现了国内含金刚石原生矿的母岩——“东方一号”岩体(D1),经过50余年的勘查,不仅在该区及毗邻区域发现了大量规模较小的钾镁煌斑岩体(群)[1],还在系统测试和综合研究的基础上,形成了一系列独特的认识:镇远地区同时具有金伯利岩和钾镁煌斑岩,镇远马坪地区的含金刚石母岩为角砾凝灰质金伯利岩[2];在镇远马坪一带发现的300余个岩体中,至少有52个岩体含金刚石矿,肯定了该岩体源区具有形成金刚石的条件[3];近几年,随着大比例尺岩体调查不断推进,不仅在马坪地区菖蒲塘地表发现了钾镁煌斑岩(D33)岩管[4],还提出马坪地区岩体呈环带状分布,推测在其中心可能存在规模较大的隐伏岩管或岩筒。而以往该地区地表发现和钻孔揭露的岩体仅仅是隐伏岩体(岩管或岩墙)的浅部出露部分,属于金伯利岩浆体系的浅部相,而非根部相[5-6]。因此,贵州镇远马坪地区金刚石原生矿找矿具有较大的潜力,在区内识别金刚石隐伏岩管或岩筒对找矿工作十分必要。

地球物理勘探是研究地球深部电性结构的一种重要手段,其中以设备轻便、勘探深度大、工作效率高的音频大地电磁法(AMT)应用最为广泛[7-8]。AMT是以天然交变电磁场为场源,在地表同步观测其电、磁场分量,并根据不同分量之间的线性关系,在频率域获取诸如视电阻率、阻抗相位等多种响应函数,最后结合已有地质资料选取合适的反演方法得到电阻率模型,进而达到勘探的目的[9]。该方法已经在沉积层覆盖的简单和复杂地质条件下均发现了不同类型的金伯利岩岩脉,在金刚石原生矿勘探中发挥了重要作用[10-11]。

近年来,随着大地电磁法三维正演技术快速发展[12],适合类似贵州山区这种复杂地电模型的三维模拟已不再是难题[13-14],能有效避免一、二维反演静位移假象的三维反演方法也越来越成熟[15-16]。故为高效地探测贵州镇远马坪地区金刚石隐伏岩管或岩筒,开展区内金刚石原生矿找矿预测与评价工作,大比例尺密集网度的音频大地电磁测深工作显得尤为关键。

1 研究区地质背景

研究区位于扬子陆块南东边缘部,EW向镇远—贵阳深大断裂与NW向松桃—三都深大断裂交汇处之南西侧。区域上地表岩体产出和分布多受EW、NE向次级断裂构造和层间虚脱空间控制[17],马坪地区所发现的各类岩体群沿EW向深冲断层(F7)及其次级断层密集分布,大部分岩体走向与断层走向一致,部分与断层呈锐角斜交;岩体呈岩墙式岩脉、岩床式岩脉、网状式岩脉侵入寒武系上统石冷水组(∈3s)及芙蓉统娄山关组(∈4ls)的白云岩中[4](图1),其形态和规模与侵入围岩的单层厚度有一定关系[2]。区内产出的金刚石多为Ⅱa型,国际上普遍认为该类型金刚石形成于地幔超深处或熔融的碳酸盐岩与上地幔的相互反应中[18⇓-20],揭示了该地区原生金刚石形成源区深度较深。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   研究区地质及AMT剖面布设

Fig.1   The geology background and location of AMT profile


2 研究区地球物理特征

研究区内偏碱性的基性-超基性煌斑岩平均密度为2.76×103 kg/m3,与围岩有0.12×103 kg/m3的密度差,磁化率为(35~126)×10-6SI,地面磁场强度最高达110 nT,相比于围岩,岩体表现为高密度、正磁性特征[21]。区域上,研究区位于1个大型圈闭似环状航磁异常内,与重力正异常有较好的对应及套合关系。该处重磁异常主要为与原生金刚石赋存有关的岩体所引起,亦或是岩筒异常引起[22];且该重磁异常显示镇远马坪—施秉地区出露的岩体为浅部岩枝异常[23],也为研究区处于“浅部相”的地质环境提供了有利的支撑。

运用直流电法仪在野外进行标本实测,包括岩石露头187处、坑道内围岩和岩体20处、钻孔岩心66件,总结了区内岩(矿)石电阻率特征(表1)。

表1   测区岩(矿)石电阻率特征

Table 1  Resistivity characteristics of rock (ore) in the study area

岩性地层代号标本
数/件		电阻率/(Ω·m)
变化范围平均值
黏土Q3073~178108.3
白云岩∈4ls、∈3s、∈3g、Z66721~3 0262 170.2
灰岩∈2q38877~1 8731 383.5
粉砂质页岩∈2p31226~607407.6
炭质页岩∈2jm、Nh2d30194~298262.8
含砾砂岩Nh3n30706~1 9561 486.2
强风化岩体/3883~313223.3
岩体/10527~826692.5

新窗口打开| 下载CSV


根据表1的实测数据,结合地层岩性综合分析,划分了研究区的电性层:寒武系娄山关组(∈4ls)、石冷水组(∈3s)、高台组(∈3g)、清虚洞组(∈2q)岩性以白云岩、灰岩为主,平均电阻率大于1 300 Ω·m,表现为中—高电阻率特征,该岩性组合层为区内第一电性层;寒武系中统杷榔组(∈2p)和九门冲组(∈2jm)岩性以粉砂质页岩为主,夹少量炭质页岩,平均电阻率为262~407 Ω·m,表现为低—中电阻率特征,为区内第二电性层;震旦系的白云岩以及南华系上统南坨组(Nh3n)的含砾砂岩,平均电阻率大于1 300 Ω·m,表现为中—高电阻率特征,为区内第三电性层;虽然南华系中统大塘坡组(Nh2d)以粉砂质页岩为主,但其埋藏较深,且地层厚度较薄(<180 m),在剖面图上很难区分开,故将该岩性层划分到区内第三电性层(图2)。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   研究区内电性结构特征

Fig.2   The resistivity structure map of study region


从地表调查和钻孔揭露资料分析,除研究区内“东风一号”岩体(D1)外,其余岩体在区内主要以强风化的形态呈现,实测电阻率较低,综合表现为低阻特征,与围岩差异较明显。

3 AMT数据采集与信噪分析

AMT测量使用加拿大凤凰公司的V8多功能电法站,剖面按SN向布设,大致垂直于深冲断层(F7),野外数据采集按80 m×40 m网度进行,采用 “十”字型观测装置,观测频率范围0.35~10 400 Hz,低通滤波(LP Filter)为1档,磁道和电道增益均为1。为增加视电阻率及阻抗相位原始数据的叠加次数,保证相应曲线圆滑连续,根据观测时间长度实验,最终确定测深点的最少观测时间为75 min。

通过2个磁场分量和2个电场分量获得了阻抗张量,图3所示为观测参考道和处理后的张量预测值之比(信噪比)。AMT测点随频率变化的电场与磁场分量之间的相干度较高,接近于1(图3b),表明工区噪声水平较低,数据质量可靠。但也存在极个别相干度稍差的测点,如L10测线的22号点(图3a),可以看到5 000~1 000 Hz频段因为“死频带”、400~130 Hz频段因为人文噪声,数据相干度稍小但总体超过0.8。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   AMT测点观测参考道和预测值间的相干度

Fig.3   Coherence map between observed reference path and predicted value of AMT measuring point


数据采集和处理过程中,采取了以下措施进行噪声压制:①野外数据采集前,开展仪器和磁棒标定工作,尽量消除仪器自身漂移;②数据采集使用极差配对小于2 mV的不极化电极,改善DC漂移;③数据处理中,采用时窗傅里叶变换互相重叠,采用基于全局相干度自动robust选谱过程,得到初步处理结果,然后基于视电阻率方差指标再进行细致选谱,最终人工剔除每个频点偏差较大的互功率谱,尽量使视电阻率曲线和相位曲线平滑。

4 基于三维正演的AMT阻抗相位校正技术

虽然近年来MT、AMT技术的应用成果表明阻抗相位数据不容易受静态位移影响[24],但本次在复杂地形下开展AMT二维、三维带地形正演研究的结果表明:在定性解释中,阻抗相位是非常重要且客观的参数,无论初始模型深部有无低阻含金刚石岩管或岩筒, TE模式(电场沿构造走向极化)和TM模式(电场垂直构造走向极化)的阻抗相位和其旋转不变量均受地形影响。鉴于此,利用区内地形数据和地质资料建立了2个三维模型(一个为背景模型,另一个为低阻岩筒异常模型),在三维正演的基础上,提出了AMT阻抗相位校正技术[25],实测资料经校正后能有效圈出低阻岩筒的水平边界。

如图4所示,研究区地形起伏非常大,在1 km2范围内最大高差达200 m,构造亦十分复杂。从典型测点曲线图和剖面拟断面图可知(图5),实测的AMT数据静态位移较严重,不利于一维、二维反演,易让解释人员得出错误结论。若不消除起伏地形的影响,无法根据实测数据的阻抗相位不变量来定性解释深部是否存在低阻岩管或岩筒。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   测区等高线及AMT剖面布设

Fig.4   The elevation map and location of AMT profile


图5

新窗口打开| 下载原图ZIP| 生成PPT

图5   实测数据静态位移

Fig.5   Static displacement graph of measured data


阻抗相位不变量校正需使用实际地形和背景电阻率数据进行三维正演,得到纯地形响应,然后在实测数据中去除。正演模型为实际测点80 m×40 m网格的纯地形模型,根据栅格化的地形数据建立 1 500 m以浅的三维网格,地表以下的网格内电阻率填充为岩性露头和钻孔岩心实测的地层电阻率,地表以上的网格内均填充为空气。对该三维模型除顶部边界外的其余5个边界,均使用2倍放大因子等比扩展9个网格,采用三维交错采样网格有限差分法[26-27],开展与实测频率一致的60个频点(0.35~10 400 Hz)的正演计算,得到TE和TM两个极化模式的三维阻抗张量,并依此来计算阻抗相位不变量(P2,如图6所示);在实测资料相同频率的阻抗相位不变量(P1)的基础上,通过简单的校正公式来获取校正后的阻抗相位不变量(P3)[25]:

P3=P1-P2+45° 。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6   纯地形背景的三维阻抗相位不变量平面

Fig.6   Three-dimensional impedance phase invariant plan view of pure terrain background


图7、图8为地形背景模型相位不变量校正前、后的成果,为提高解释分辨率,对大于60°的相位值进行了剔除。图中黑色框(T1~T5)显示该区域受地形影响引起的高相位异常被有效去除,特别是去除了T1和T4区域的高相位圈闭假异常,L07-L09线中部、L10-L11线中上部的高相位圈闭异常面积也相应缩小,南西角和 L15线剖面中部的高相位圈闭异常更加明显,其余单条剖面的零星圈闭异常也得到了很好的压制,有效避免了静位移假象。因此,从校正后的成果图(图7b、图8b)可大致判定出北西角和南东角的地层趋向于稳定,而疑似管状体的异常出现在NE向的对角线上。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7   实测数据900 Hz阻抗相位不变量校正前后水平切片

Fig.7   Measured data 900 Hz impedance phase invariant correction before and after horizontal slices


图8

新窗口打开| 下载原图ZIP| 生成PPT

图8   实测数据530 Hz阻抗相位不变量校正前后水平切片

Fig.8   Measured data 530 Hz impedance phase invariant correction before and after horizontal slices


5 基于三维反演的电性结构精细分析

采用中国地质科学院阮帅博士开发的AR-QN算法及相关代码[24,28-29],对实测数据进行三维反演。对参与反演对输入数据做如下选择:①由于大量数据在5 000~1 000 Hz的“死频带”及低于1 Hz的低频带上质量较差,故选择1 000~2 Hz的30个频点(图9); ②4个阻抗张量元素均参与反演,并按视电阻率相对误差5%设定误差门槛。反演使用与前文(第4节)一样的网格剖分和地形数据,初始模型为1 000 Ω·m均匀空间。由于AR-QN算法对其反演目标函数采用了LBFGS公式单独逼近数据拟合项,保留了正则光滑项的精确形式,大量模型合成数据和实测数据的反演测试表明该方法初始模型依赖性不强,有极好的噪声稳健性[30]。

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9   测深点反演电阻率和相位曲线拟合

Fig.9   Inversion resistivity and phase curve fitting diagram of sounding point


在本研究区的数据反演中,使用比值法自适应更新正则因子[31],最终经45次迭代得到了平均粗糙度1.5、加权残差均方根2.1的模型。必须指出的是,研究中同样尝试了初始模型为100 Ω·m和 5 000 Ω·m的模型,在自适应正则因子的AR-QN算法下,最终结果和目标函数收敛水平基本类似。将正则因子固定为很小的数(如0.01)后,虽然拟合差会小一些,但在海拔600 m以上会形成过多冗余构造。从“简单模型”的原则上考虑,最终选择了自动算法的结果。从典型的测深点曲线拟合效果(图9)分析,绝大部分频点数据拟合效果很好,仅少数高频点、低频点数据拟合效果欠佳,说明反演模型对于原始数据进行的拟合是有效的,反演成果是合理可靠的。

三维反演根据研究区岩性统计结果设定地下单元的电阻率变化区间,获得了可靠的三维电性结构(图10),整体上呈2层电性结构特征(考虑AMT测量深度,成果图未下延至物性成果的第三电性层),海拔500 m以浅的高阻层应是以碳酸盐岩为主的寒武系清虚洞组至娄山关组,为区内电性结构中的第一电性层;海拔100~500 m的低阻层对应的是以碎屑岩为主的寒武系九门冲组至杷榔组,为区内电性结构中的第二电性层;图9所示的电阻率高阻与低阻界面倾斜方向,揭示地层倾向为SWW向。

图10

新窗口打开| 下载原图ZIP| 生成PPT

图10   研究区三维AMT反演电阻率三维电性结构

Fig.10   3D AMT inversion electrical structure diagram of the study area


图10中存在一个低阻筒状体异常,在寒武系白云岩、灰岩等高阻体地层中边界差异明显,在碎屑岩中,由于电阻率较为接近,仅依据电阻率参10难分辨该低阻筒状体异常的边界,结合实测物性及深部钻孔资料分析,该低阻筒状体异常是金刚石隐伏岩管或岩筒的可能性较大。该低阻体中心顶部发育于研究区北东部,即L11~L15线之间的小号测点(7~11号点)区域,深部发育渐往南东部(9~14号点)偏移(图11所示),为研究区处于“浅部相”的地质环境提供了有利的地球物理信息。

图11

新窗口打开| 下载原图ZIP| 生成PPT

图11   实测数据三维反演综合解释成果

Fig.11   3D inversion comprehensive interpretation slice map of measured data


从图10分析可知,该低阻异常体顶部被厚度不一的高阻体覆盖,说明该异常体未发育至地表,从区内历年勘查资料及区内深部钻孔(图11)揭示,高、低阻的起伏界面即为碳酸盐岩地层与岩筒的不整合接触面;而地表零星的低阻异常区块,是区内偏碱性基性—超基性煌斑岩风化后和第四系沉积物的综合反映。

图11是电阻率三维反演后水平切片与校正后阻抗相位等值线叠合成果,为提高解释分辨率,阻抗相位等值线只绘制与电阻率同深度的高值部分。分析可知,深冲断裂F7体现在电阻率界面不连续处,下切深度大,产状陡(>80°);而发育于研究区南部、EW走向的断裂和区域断裂F6更多地体现在阻抗相位梯度变化界面。为避免因电阻率单一参数划分筒状体异常的边界,在电阻率不同深度切片、同一平面投影位置的高阻抗相位区域,圈定了一个低电阻率、高阻抗相位异常区Y1,位于F7断层的北东部,推测该异常区是金刚石隐伏岩管或岩筒的发育位置,具有浅部较小、深部较大的发育特征;区内深部钻孔均在Y1异常南西部,主要揭露了多条岩脉,未发现岩筒存在,综合分析认为受后期次级小断裂和碳酸盐岩地层“刚性”强度的影响,低阻筒状体异常区Y1附近的其他高阻抗相位异常区是岩体呈岩墙式岩脉、岩床式岩脉、网状式岩脉侵入碳酸盐岩地层的可能性较大;这些高阻抗相位异常区大多沿深冲断裂F7走向分布,且深部低电阻率异常中心伴随断裂F7发育,这说明深冲断裂F7在发育过程中切穿了隐伏岩筒,成为岩筒内岩体运移至地表的重要通道,进而形成了地表的偏碱性基性超基性煌斑岩(图1所示)沿深冲断裂F7走向分布。

6 结论

1)分析区内岩(矿)石电阻率数据,建立了区内三层电性组合模型;在复杂地形条件下,利用三维正演技术模拟研究区纯地形响应并在实测数据中去除,得到的定性解释结果在一定程度上恢复了被静态效应扭曲的AMT阻抗相位不变量分布形态,增强了定性解释的可靠性。

2)在考虑AMT有效探测深度的基础上, 对实测数据的阻抗张量全参数进行了三维AR-QN反演并得到了较可靠的研究区电性结构。三维电性结构揭示了一个带低阻筒状体、具有2层结构特征的三维地电模型。

3)在电阻率三维反演水平切片成果配合校正后阻抗相位等值线解释的基础上,利用低电阻率、高阻抗相位结合的模式,对含金刚石隐伏岩管或岩筒发育位置进行了初步判断,为区内进一步的矿产勘查和科学研究提供了地球物理信息。

致谢

中国地质调查局南京地质调查中心杨献忠教授级高级工程师、中国地质大学(武汉)高凌峰博士提出了宝贵的修改意见,在此表示衷心感谢。

参考文献

杨光忠.

贵州镇远地区钾镁煌斑岩产出控制因素浅析

[J]. 地质与勘探, 2013, 49(4):696-702.

[本文引用: 1]

Yang G Z.

Controlling factors of lamprophyre’s occurrence in the Zhenyuan area of Guizhou Province

[J]. Geology and Exploration, 2013, 49(4): 696-702.

[本文引用: 1]

李永刚, 向璐, 黄远成, 等.

贵州镇远地区含金刚石母岩再认识

[J]. 地质通报, 2019, 38(1):103-109.

[本文引用: 2]

Li Y G, Xiang L, Huang Y C, et al.

Re-understanding of diamond-bearing parent rocks in Zhenyuan area, Guizhou Province

[J]. Geological Bulletin of China, 2019, 38(1): 103-109.

[本文引用: 2]

李永刚. 贵州镇远地区含金刚石母岩岩石学及含矿性研究[D]. 武汉: 中国地质大学(武汉), 2019.

[本文引用: 1]

Li Y G. Petrological and ore-bearingproperties of the diamond-bearing pluton in Zhenyuan Area of Guizhou Province[D]. Wuhan: China University of Geosciences(Wuhan), 2019.

[本文引用: 1]

黄远成, 石睿, 林泽渊, 等.

贵州镇远苍蒲塘钾镁煌斑岩管发现及找矿意义

[J]. 贵州地质, 2015, 32(1):32-36.

[本文引用: 2]

Huang Y C, Shi R, Lin Z Y, et al.

Discovery of lamproite tube and Its prospecting significance in Cangputang of Zhenyuan, Guizhou Province

[J]. Guizhou Geology, 2015, 32(1): 32-36.

[本文引用: 2]

张锡贵, 石睿, 吴寿宁, 等.

贵州施秉翁哨地区钾镁煌斑岩的新发现及其金刚石找矿意义

[J]. 贵州地质, 2015, 32(1):37-40.

[本文引用: 1]

Zhang X G, Shi R, Wu S N, et al.

New discovery of lamproite and its significance for diamond exploration in Wengshao Area of Shibing,Guizhou Province

[J]. Guizhou Geology, 2015, 32(1): 37-40.

[本文引用: 1]

黄远成, 李志翔, 丘志力, 等.

贵州镇远钾镁煌斑岩原生及砂矿金刚石矿物学特征及其找矿意义

[J]. 中山大学学报:自然科学版, 2016, 55(5):108-118.

[本文引用: 1]

Huang Y C, Li Z X, Qiu Z L, et al.

Mineralogical characteristics oflamproite-hosted and placer diamonds from Zhenyuan, Guizhou and their significance for primary deposit prospecting

[J]. Acta Scientiarum Naturalium Universitatis Sunyatseni, 2016, 55(5): 108-118.

[本文引用: 1]

田占峰, 毛星, 罗旭, 等.

音频大地电磁测深法在电性结构研究中的应用——以郯庐断裂带宿迁段为例

[J]. 物探与化探, 2016, 40(4):732-736.

[本文引用: 1]

Tian Z F, Mao X, Luo X, et al.

The application of AMT method to the electrical structure study: A case study of Suqian sector of Tan-Lu Fault Zone

[J]. Geophysical and Geochemical Exploration, 2016, 40(4): 732-736.

[本文引用: 1]

何帅, 杨炳南, 李核良, 等.

音频大地电磁法对渝东南Ⅳ级地堑构造的识别及意义

[J]. 地质科技情报, 2019, 38(1):270-276.

[本文引用: 1]

He S, Yang B N, Li H L, et al.

Identification of Ⅳ graben tectionics of southeast Chongqing by AMT method and its significanc

[J]. Geological Science and Technology Information, 2019, 38(1): 270-276.

[本文引用: 1]

Simpson F, Bahr K. Practical magnetotellurics[M]. Cambridge: Cambridge University Press, 2005.

[本文引用: 1]

Saraev A K, Larionov K A.

音频大地电磁测深在金伯利岩勘探中的应用

[J]. 石油地球物理勘探, 2004, 39(S1):144-145.

[本文引用: 1]

Saraev A K, Larionov K A.

Application of audio geomagnetic bathymetry in Kimberley rock exploration

[J]. Petroleum Geophysical Exploration, 2004, 39(S1): 144-145.

[本文引用: 1]

Saraev A K, Antaschuk K M, Nikiforov A B, 等.

AMT测深法在金刚石矿勘探中的应用

[J]地球物理学报, 2010, 53(3): 657-676.

[本文引用: 1]

Saraev A K, Antaschuk K M, Nikiforov A B, et al.

Audiomagnetotelluric soundings for the diamond exploration

[J]. Chinese J. Geophys., 2010, 53(3): 657-676.

[本文引用: 1]

Zhdanov M S, Fang S, Hursan G.

Electromagnetic inversion using quasi-linear approximation

[J]. Geophysics, 2000b, 65(5): 1501-1513.

DOI:10.1190/1.1444839      URL     [本文引用: 1]

Siripunvaraporn W, Egbert G, Lenbury Y.

Numerical accuracy of magnetotelluric modeling:a comparison of finite difference approximations

[J]. Earth Planets Space, 2002, 54(6): 721-725.

DOI:10.1186/BF03351724      URL     [本文引用: 1]

梁苗, 刘双, 胡祥云.

基于断层模型的AMT数值模拟及其应用

[J]. 地质科技情报, 2016, 35(5):231-237.

[本文引用: 1]

Liang M, Liu S, Hu X Y.

AMT numerical simulation of fault model and its applications

[J]. Geological Science and Technology Information, 2016, 35(5): 231-237.

[本文引用: 1]

Zhdanov M S, Tolstaya E.

Minimum support nonlinear parametrization in the solution of a 3D magnetotelluric inverse problem

[J]. Inverse Problems, 2004, 20(3): 937-952.

DOI:10.1088/0266-5611/20/3/017      URL     [本文引用: 1]

汤井田, 张林成, 王显莹.

庐枞矿集区矾山—将军庙地区AMT三维反演及地质结构解释

[J]. 地球物理学报, 2018, 61(4):1576-1587.

[本文引用: 1]

Tang J T, Zhang L C, Wang X Y, et al.

Subsurface electrical structure of the Fanshan-Jiangjunmiao region in the Lujiang-Zongyang Ore District derived from 3-D inversion of audio-magnetotelluric data

[J]. Chinese J. Geophys., 2018, 61(4): 1576-1587.

[本文引用: 1]

孙士军, 杨松平.

贵州金刚石成矿条件初探

[J]. 贵州地质, 1998, 15(1):1-8.

[本文引用: 1]

Sun S J, Yang S P.

A discussion on prerequisites to search for primary diamondsin Guizhou

[J]. Guizhou Geology, 1998, 15(1): 1-8.

[本文引用: 1]

Hayman P C, Kopylova M G, Kaminsky F V.

Lower mantle diamonds from Rio Soriso (Juina area, Mato Grosso,Brazil)

[J]. Contributions to Mineralogy & Petrology, 2005, 149(4): 430-445.

[本文引用: 1]

Tappert R, Stachel T, Harris J W, et al.

Diamonds from Jagersfontein (South Africa): Messengers from the sublithospheric mantle

[J]. Contributions to Mineralogy & Petrology, 2005, 150(5): 505-522.

[本文引用: 1]

Bulanova G P, Walter M J, Smith C B, et al.

Mineral inclusions in sublithospheric diamonds from Collier 4 kimberlite pipe, Juina, Brazil: Subducted protoliths, carbonated melts and primary kimberlite magmatism

[J]. Contributions to Mineralogy and Petrology, 2010, 160(4): 489-510.

DOI:10.1007/s00410-010-0490-6      URL     [本文引用: 1]

王亮, 陶平.

贵州东南部含金刚石钾镁煌斑岩找矿远景区预测

[J]. 地质与勘探, 2012, 48(4):775-783.

[本文引用: 1]

Wang L, Tao P.

Ore-search prospecting areas of diamond-bearing lamproite in southeastern Guizhou Province

[J]. Geology and Exploration, 2012, 48(4): 775-783.

[本文引用: 1]

王亮, 陶平.

利用区域物探重磁资料圈定黔东金刚石母岩钾镁煌斑岩的尝试

[J]. 贵州地质, 2011, 28(4):254-259.

[本文引用: 1]

Wang L, Tao P.

Determination of Diamond Parent Rock in Eastern Guizhou By Regional Geophysical Prospecting Gravity and Magnetic Information

[J]. Guizhou Geology, 2011, 28(4): 254-259.

[本文引用: 1]

樊洪富, 王亮.

黔东南火山构造特征与金刚石找矿前景分析

[J]. 工程地球物理学报, 2016, 13(3):389-398.

[本文引用: 1]

Fan H F, Wang L.

Volcano-tectonic characteristics of diamond ore prospects in southeast Guizhou

[J]. Chinese Journal of Engineering Geophysics, 2016, 13(3): 389-398.

[本文引用: 1]

阮帅. 三维大地电磁有限内存拟牛顿反演[D]. 成都: 成都理工大学, 2015.

[本文引用: 2]

Ruan S. Three-dimensional geomagnetic finite memory is intended for Newtonian inversion[D]. Chengdu: Chengdu University of Technology, 2015.

[本文引用: 2]

阮帅, 张炯, 孙远彬, 等.

基于三维正演的音频大地电磁阻抗相位不变量校正技术

[J]. 地球物理学报, 2015, 58(2):685-696.

[本文引用: 2]

Ruan S, Zhang J, Sun Y B, et al.

AMT impedance phase invariant correction based on 3D MT modeling technology

[J]. Chinese J. Geophy., 2015, 58(2): 685-696.

[本文引用: 2]

Mackie R L, Madden T R, Wannamaker P E.

Three-dimensional magnetotelluric modeling using difference equations: Theory and comparisons to integral equation solutions

[J]. Geophysics, 1993, 58(2): 215-226.

DOI:10.1190/1.1443407      URL     [本文引用: 1]

Newman G A, Alumbaugh D L.

Three-dimensional magnetotelluric inversion usingnon-linear conjugate gradients

[J]. Geophys. J. Int., 2000, 140(2): 410-424.

DOI:10.1046/j.1365-246x.2000.00007.x      URL     [本文引用: 1]

Byrd R H, Nocedal J, Schnabel R B.

Representations of quasi-Newton matrices and their use in limited memory methods

[J]. Mathematical Programming, 1994, 63(1-3): 129-156.

DOI:10.1007/BF01582063      URL     [本文引用: 1]

邓琰, 汤吉, 阮帅.

三维大地电磁自适应正则化有限内存拟牛顿反演

[J]. 地球物理学报, 2019, 62(9):3601-3614.

[本文引用: 1]

Deng Y, Tang J, Ruan S.

Adaptive regularized three-dimensional magnetotelluric inversion based on the LBFGS quasi-Newton method

[J]. Chinese J. Geophy., 2019, 62(9): 3601-3614.

[本文引用: 1]

阮帅, 汤吉, 陈小斌, 等.

三维大地电磁自适应L1范数正则化反演

[J]. 地球物理学报, 2020, 63(10):3896-3911.

[本文引用: 1]

Ruan S, Tang J, Chen X B, et al.

Three-dimensional magnetotelluric inversion based on adaptive L1-norm regularization

[J]. Chinese J Geophy, 2020, 63(10): 3896-3911.

[本文引用: 1]

陈小斌, 赵国泽, 汤吉, 等.

大地电磁自适应正则化反演算法

[J]. 地球物理学报, 2005, 48(4):937-946.

[本文引用: 1]

Chen X B, Zhao G Z, Tang J, et al.

An adaptive regularized inversion algorithm for magnetotelluric data

[J]. Chinese J. Geophy., 2005, 48(4): 937-946.

[本文引用: 1]

/

〈 〉

京ICP备05055290号-3
版权所有 © 2021《物探与化探》编辑部
通讯地址：北京市学院路29号航遥中心 邮编：100083
电话：010-62060192；62060193 E-mail：whtbjb@sina.com