Tracing infiltration and recharge of thick silt by using D and 18O isotopes of soil moisture in Xiaogan, Hubei and its ecological efffects
-
摘要:
黏性土的渗透系数极低,水分及溶质在黏性土中运移速率慢、耗时长,本次研究通过分析大别山区-江汉平原三水转换野外科学试验场(下文简称"试验场")ZK1、ZK2钻孔剖面土壤水、大气降雨D、18O同位素测试数据与孝感站(站号57482)多年年降雨量数据,确定了厚层黏性土土壤水入渗补给年份与深度的对应关系。结果表明:试验场区黏性土垂向岩性差异较小,无明显分层现象,土壤水分以"活塞流"的方式向下运移,夏、秋季的大气降雨为土壤水的主要补给来源;ZK1(取样间隔0.5~2.7 m,深度15 m)的土壤水δD、δ18O值随着埋深的增大出现周期性的波动,ZK2(取样间隔0.1 m,深度6.2 m)的土壤水δD、δ18O值随着埋深的增大出现分层波动现象;确定了黏性土层0~6.2 m深度对应的降雨入渗补给年份,并通过18O的峰值位移法计算得出降雨入渗补给在黏性土层的垂向运移速度为10.8~15.0 cm/a,年均入渗补给量为43.1~58.1 mm,占多年年均降雨量的4.01%,推算出降雨入渗补给需要近130年的时间才能穿透试验场厚层黏性土补给至地下水含水层,表明该厚层黏性土的防污性能良好。本研究所揭示大气-土壤界面下黏性土土壤水分入渗迁移历史演化特征及补给年际对应关系,对江汉平原区地下水环境保护、生态环境改善、旱涝灾害防治等具有重要意义。
Abstract:As well known, the permeability coefficient of cohesive soil is extremely low, and the transport rate of water and solute in the cohesive soil is slow and time-consuming. Stable isotopes δD and δ18O data of soil moisture and rainfall of two boreholes (ZK1, ZK2) in Dabie Mountain Area-Jianghan Plain rainfall-soil moisture-groundwater transformation scientific field test site were studied to trace the characteristics of historical infiltration and recharge of cohesive soil. The results show that there's no obvious layer boundary but just a little difference among the 15 m deep cohesive soil, and soil water moves in the way of "piston flow" vertically. The closer relation of rainfalls in summer and autumn to soil moisture in δD and δ18O value indicates that the main recharge of soil moisture comes from rainfall infiltration in summer and autumn. The δD and δ18O values of soil water in ZK1 (sampling interval 0.5-2.7 m, depth 15 m) fluctuate periodically with the increase of buried depth, while the δD and δ18O values of soil water in ZK2 (sampling interval 0.1m, depth 6.2 m) fluctuate stratified with the increase of buried depth. On the basis of determining the rainfall infiltration recharge years corresponding to the 0-6.2 m depth in the cohesive soil layer, based on 18O isotopes peak displacement method, it is calculated that the vertical migration velocity of rainfall infiltration recharge in the cohesive soil layer is 10.8-15.0 cm/a, and the annual infiltration recharge is 43.1-58.1 mm, accounting for 4.01% of the annual average rainfall. It takes about 130 years for rainfall infiltration to penetrate the thick cohesive soil for groundwater recharge, which means the thick cohesive soil has good antifouling properties. The historical evolution characteristics of soil water infiltration and migration at the air-soil interface and the inter-annual correspondence of soil water recharge are of great significance for groundwater environmental protection, ecological environment improvement, drought and flood disaster prevention and control in Jianghan Plain.
-
1. 引言
遥感地质找矿是遥感信息获取、含矿信息提取以及含矿信息成矿分析与应用的过程, 近年来, 被广泛应用于地质、矿产资源及相关环境的调查中(刘燕君等, 1993; 朱振海, 1994; 郭德方, 1995; Sabins, 1999; Everett et al., 2002; 廖崇高等, 2002; 连长云等, 2005; 刘颖璠等, 2012; 钱建平等, 2012; 谭克龙等, 2012; 杜小弟等, 2015)。航空高光谱是当前遥感的前沿技术, 通过高光谱成像所获取的地球表面的图像包含光谱维信息融合为一体, 即“图谱合一” (束炯等, 2006), 具有波段数多, 波带窄, 对地物能定量分析等优点, 是20世纪末期以来遥感领域最大的技术进展(Campbell et al., 2011)。航空高光谱遥感技术相对于星载高光谱遥感技术, 可以获取高空间分辨率的高光谱遥感数据, 因而对微小地物具有更强的识别能力(刘德长等, 2015)。通过高光谱矿物填图, 可以大面积、快速提取蚀变矿物(Clack et al., 1990)。王润生等(2010)系统总结了高光谱矿物填图技术流程、工作方法和技术体系。中国国土资源航空物探遥感中心在新疆土屋东—三岔口地区, 应用HyMap航空成像光谱数据填绘出了白云母、绿泥石、绿帘石、绿泥石和绿帘石组合、高岭石、蒙脱石、透辉石、透闪石、蛇纹石、褐铁矿、方解石等矿物或矿物组合种类(王润生等, 2011)。
本文以新疆卡拉塔格地区获取的HyMap航空高光谱数据、地面准同步定标数据、地面高光谱数据为例, 提取并筛选了基于HyMap航空高光谱影像的成矿有利蚀变信息, 在此基础上对红山铜金矿床航空-地面高光谱提取的蚀变信息进行了综合剖析, 并结合矿区地质背景, 建立了红山铜金矿床高光谱遥感地空综合找矿模型, 圈定了2处找矿有利地段。
2. 研究区地质概况
研究区大地构造位于东天山吐哈盆地南缘、大南湖—头苏泉晚古生代岛弧带(图 1)。
图 1 新疆卡拉塔格地区地质矿产图1—全新统:砂石及砾石沉积; 2—下侏罗统三工河组:粉砂质泥岩、泥质砂岩夹钙质砂岩; 3—上二叠统库莱组:砾岩、中粗粒岩屑砂岩; 4—中二叠统阿尔巴萨依组:凝灰岩、橄榄玄武岩; 5—上石炭统脐山组:玄武岩、安山岩; 6—下泥盆统大南湖组:安山岩、英安岩、熔结凝灰岩、火山角砾岩; 7—中—上志留统红柳峡组:凝灰质砾岩、晶屑岩屑凝灰岩; 8—中奥陶统荒草坡群大柳沟组:玄武岩、安山岩、英安岩; 9—志留纪二长花岗岩; 10—志留纪花岗闪长岩; 11—志留纪英云闪长岩; 12—花岗岩体; 13—地质界线; 14—角度不整合界线; 15—岩相界线; 16—逆断层; 17—性质不明断层; 18—推测断层; 19—铜金矿; 20—铜锌矿; 21—铜锌金矿; 22—铜铁矿Figure 1. Geological map of Kalatag area in Xinjiang1-Holocene, sand and gravel deposits; 2-Lower Jurassic Sangonghe Formation: river silty mudstone, argillaceous sandstone intercalated with calcareous sandstone; 3-Permian Kulai Formation: breccia, coarse grained lithic sandstone; 4-Middle Permian Aebasayi Formation: tuff, olivine basalt; 5-Carboniferous Qishan Formation: basalt, andesite; 6-Lower Devonian Dananhu Formation: andesite, dacite, welded tuff and volcanic breccia; 7-Middle-upper Silurian Hongliuxia Formation: tuffaceous conglomerate, wafer crumbs debris tuff; 8-Daliugou Formation of Ordovician Huangcaopo Group: basalt, andesite, dacite; 9-Silurian monzonitic granite; 10-Silurian granite diorite; 11-Silurian tonalite; 12-Granite rock mass; 13-Geological boundary; 14-Angular unconformity; 15-Petrofacies boundary; 16-Thrust; 17-Unknown fault; 18-Inferred fault; 19-Copper and gold deposit; 20-Copper and zinc deposit; 21-Copper zinc deposit; 22-copper deposit区内地层整体呈北西向产出, 由中奥陶统荒草坡群大柳沟组(O2Hd)、中—上志留统红柳峡组(S2-3h)、下泥盆统大南湖组(D1d)、上石炭统脐山组(C2qs)、中二叠统阿尔巴萨依组(P2a)、上二叠统库莱组(P3k)、下侏罗统三工河组(J1s)、全新统(Qh)组成。
该区经历了多期次的构造变动, 区域构造形成类型多样、特征复杂, 褶皱、断裂和火山机构发育, 区域性的断裂构造控制着该区火山岩和侵入岩的展布, 次级断裂控制着矿化蚀变带的分布, 更次一级构造裂隙带控制着矿化体的产出。主要发育北西、北北西和北东东向三组断裂构造, 其中北西向断裂为该区主要控矿构造, 同火山断裂发育, 含矿火山热液沿断裂活动强烈。
区内火山活动强烈, 火山岩分布广泛。岩石类型以中基性火山岩为主, 中酸性火山岩次之, 从奥陶系至二叠系均有出露, 其中奥陶系最为发育, 其次为泥盆系、石炭系, 二叠系局部出现。侵入岩较为发育, 主要为一套古生代侵入体, 包括中酸性闪长岩、花岗闪长岩、花岗岩等侵入岩, 中部为与铜金矿化直接相关的中生代中酸性火山机构, 主要由石英斑岩、钠长斑岩、英安岩、英安斑岩、安山玢岩、闪长玢岩和火山角砾集块岩等岩性构成。其中志留纪侵入岩由老到新依次为英云闪长岩、花岗闪长岩、二长花岗岩、花岗斑岩。该区经受了区域埋深变质作用, 变质程度较弱, 主要变质矿物为绿泥石、钠长石、葡萄石、石英、阳起石、方解石、绿帘石。
区内已知矿床自北西向南东依次产出有:红山铜金矿床→梅岭铜锌金矿床→红石铜金矿床→红海铜锌金矿床。
3. 高光谱遥感地空综合预测方法用于寻找铜金矿床的原理
3.1 高光谱遥感地空综合预测方法的依据
不同矿床类型由于岩石岩性和蚀变矿物不同, 它们的光谱特征有明显的差异(姚佛军, 2006)。铜金矿床属于内生金属矿床, 这类矿床与地表蚀变密切相关, 可以利用航空高光谱遥感影像进行蚀变矿物提取, 结合地表岩石光谱测试分析, 应用HyMap航空成像光谱仪和美国ASD公司的FieldSpec Pro FR光谱仪对不同类型矿床的光谱特征进行测量, 分析其不同类型矿床的光谱特征, 进行蚀变矿物遥感信息提取, 开展高光谱遥感地空综合预测研究, 进而指导找矿。
3.2 蚀变矿物诊断性光谱特征分析
岩矿的光谱包含有一系列特征吸收谱带。每一个特征谱带或谱带组合与岩矿内部微粒的物质属性存在一定的对应关系(Povarennykh, 1978)。在短波红外光谱区域, 褐铁矿、黄钾铁矾、白云母、绿泥石等具有可以识别的诊断性光谱特征, 特征吸收波段的深度与岩石中这些矿物的含量密切相关(刘圣伟, 2006)。矿物诊断性光谱吸收包括金属阳离子在可见光区域的电子过程以及阴离子基团在近红外区域的振动过程。蚀变矿物的可见光−近红外范围内几种官能团主要包括Fe2+、Fe3+、Al−OH、Mg−OH、OH-、CO32-。
一般阴离子诊断谱带位于2000~2500 nm光谱区域。而Fe2+、Fe3+和Mn2+诊断谱带一般位于400~1200 nm光谱区域(甘甫平等, 2003)。所以根据400~1100 nm光谱区间的Fe吸收谱带和2000~2500 nm羟基和碳酸根谱带将矿物分为含铁矿物、含羟基矿物和含碳酸根矿物以及其他矿物(甘甫平, 2003)。褐铁矿的光谱吸收谷是矿物光谱分析、遥感铁染异常提取常用的波段。其中Fe3+离子的特征吸收谷位主要在600~800 nm, 中心位置在700 nm左右, Fe2+离子的特征吸收谷位主要在900~1500 nm, 中心位置在1200 nm左右; 白云母的诊断性光谱特征是Al−OH键在2180~2230 nm之间的尖锐而深的吸收特征, 以及在2340 nm和2440 nm附近的较弱Al−OH特征, Clark et al.(1990)在高光谱分辨率下, 检测到白云母矿物在2200 nm附近的吸收波段随Al含量的增加向长波方向的移动; 绿泥石的光谱具有Fe−OH和Mg−OH的诊断性吸收特征, 波长位置分别为2235~2255 nm和2320~2360 nm, 这些吸收特征的波长随着绿泥石中铁离子含量的增加而增大。
4. HyMap航空高光谱数据获取及处理
综合考虑仪器参数、成图精度的要求、研究区自然地理、地形起伏、野外工作条件等因素, 进行飞行方案的最优设计(表 1)。将HyMap成像光谱仪搭载在Y−12飞机上, 利用国外引进的HyMap航空高光谱遥感测量系统、地面ASD光谱测量系统。采集高光谱分辨率的HyMap航空高光谱遥感数据、地面准同步的地物光谱定标数据。具体获取过程分为3个步骤:首先对HyMap成像光谱仪进行实验室定标, 保证其性能指标达到项目要求, 以获取优质的高光谱影像数据; 其次进行HyMap仪器的安装和测试, 按照仪器安装规范方法, 将HyMap安装在Y−12飞机上, 保证机下的摄影窗口满足数据获取要求, 按照航空成像光谱数据获取规范进行地面测试; 最后为HyMap航空高光谱数据获取和初步质量检查。
表 1 HyMap航空成像光谱仪主要技术参数Table 1. Main technical parameters of HyMap aviation imaging spectroradiometer在此基础上, 开展HyMap航空高光谱遥感数据辐射定标、大气校正光谱反演、几何校正和地理编码等数据处理工作。通过此项研究, 为高光谱遥感地空综合预测铜金矿床的技术方法研究和应用提供可靠的数据基础。
5. 地面高光谱数据获取及处理
地面光谱测量采用美国ASD公司生产的FieldSpec Pro FR地面光谱测量仪(表 2)。利用该仪器开展了地面光谱准同步定标。在HyMap高光谱数据获取过程中, 准同步布设黑白布定标场和明暗地物定标场, 利用FieldSpec Pro FR地面光谱测量仪进行光谱定标数据获取(图 2,图 3), 并基于经验线性模型法, 对高光谱影像进行大气校正和光谱比对分析。
表 2 FieldSpec Pro FR光谱仪主要技术参数Table 2. Main technical parameters of FieldSpec Pro FR spectrometer同时开展了典型铜金矿床地表蚀变围岩和其他相关的岩石地物的光谱测量, 完成了4条贯穿地层走向的成矿有利地质单元光谱测量和4条典型矿床光谱测量。共完成光谱测量点287个, 光谱测量曲线2870条。测量完成后, 开展了各岩石类型的岩矿光谱数据处理、光谱库建立、光谱特征分析等研究工作, 研究及分析结果为航空高光谱蚀变矿物信息提取、地面光谱验证、铜金矿床围岩蚀变信息筛选等提供重要的地面光谱信息支持。
地面高光谱数据处理采用ViewSpec Pro软件, 利用ViewSpec Pro软件挑选符合要求的地面岩矿光谱曲线, 并将这些曲线转换为txt格式的文件。由于受到太阳位置、角度条件、大气条件、地形影响及传感器本身性能的影响, 传感器所记录的数据与目标的光谱反射率或光谱辐射亮度值并不一致。需先将传感器记录的原始辐射值(DN值)转化为地物反射率(李志忠等, 2015)。
6. 蚀变矿物信息提取处理方法及分布特征
目前基于高光谱遥感影像数据提取矿物的处理方法可以分为2类: (1)通过像元光谱与矿物标准光谱曲线的匹配实现识别和量化目标矿物; (2)提取矿物光谱特征吸收波段参数作为识别和量化指标(卢燕, 2014)。前者在处理过程中光谱数据信息可能由于重定标而部分丢失。且基于单幅图像统计参数的处理方法难以做到无缝衔接, 该方法也涉及大量的运算。本次研究区蚀变矿物提取主要基于第2种方法。该方法利用ENVI软件, 基于可编程实现的特征提取处理通道, 该通道中脚本设计是基于提取光谱参数, 首先对影像光谱数据进行均值标准化、去连续统、卷积平滑等数据预处理, 在进行特征深度、特征面积、特征宽度、极小值波长、比值、算术和逻辑运算等特征参数提取。同时利用多个诊断性特征来标记或约束特定的矿物, 提取矿物的分布和相对含量, 并增强矿物识别的精度。
根据矿物光谱特征吸收波段参数的蚀变矿物提取方法, 在研究区提取了褐铁矿−黄钾铁矾、白云母−蒙脱石、绿泥石−绿帘石−碳酸盐、叶腊石等蚀变矿物或组合(图 4), 其中褐铁矿−黄钾铁矾、白云母−蒙脱石、绿泥石−绿帘石−碳酸盐等蚀变矿物组合在全区范围内分布, 且呈团块状、条带状及星散状。叶腊石则仅在研究区北东部呈星点状分布, 蚀变不明显, 对于指导找矿意义不大(黄色圈内为叶腊石分布位置)。
图 4 新疆卡拉塔格地区HyMap航空高光谱蚀变矿物丰度图a—褐铁矿-黄钾铁矾; b—白云母-蒙脱石; c—绿泥石、绿帘石、碳酸盐; d—叶腊石, 黄色圈内为叶腊石分布位置Figure 4. HyMap aviation hyperspectral alteration mineral abundance map of Kalatag area in Xinjianga-Limonite-jarosite; b-Muscovite-montmorillonite; c-Chlorite, epidote, carbonate; d-Pyrophyllite, yellow circle indicates pyrophyllite distribution position7. 红山铜金矿床高光谱遥感地空综合找矿模型的建立与找矿有利地段优选
7.1 矿区地质特征
红山铜金矿床是高硫型−浅成低温热液型矿床(董承维, 2013)。矿床产于陆相蚀变酸性火山穹窿及其通道中, 该次火山机构受近东西与北西向两组断裂控制, 且主要控岩控矿构造为北西向断裂。铜矿体受北西向断裂及其次级构造破碎带控制, 为浸染状辉铜矿化和黄铜矿化。金矿化主要与酸性次火山岩中的细粒黄铁矿有关(王瑞军等, 2014) (图 5)。
图 5 红山铜金矿区地质图(据王京彬, 2006)1—第四系残坡积物、风积物; 2—英安岩、英安斑岩、安山岩; 3—蚀变英安岩; 4—流纹斑岩、石英斑岩、霏细斑岩; 5—花岗岩; 6—花岗闪长岩脉; 7—闪长玢岩脉; 8—氧化带矿体; 9—推测断层; 10—片理化带; 11—金矿体; 12—铜矿体Figure 5. Geological map of the Hongshan Cu-Au deposit (after Wang, 2006)1-quaternary debris slope sediments, aeolian; 2-dacite, British Ann porphyry, andesite; 3-alteration dacite; 4-rhyolitic porphyry and quartz porphyry, faye fine porphyry; 5-granite; 6-granite diorite vein; 7-diorite Bin dike; 8-oxidation zone ore body; 9-speculation fault; 10-piece of physics and chemistry; 11-gold ore body; 12-copper ore body矿区地表发育黄钾铁钒化、绢云母化、硅化、褐铁矿化、黄铁矿化、绿泥石化、绿帘石化、青磐岩化。近矿蚀变以火山岩酸淋滤带蚀变体的黄钾铁矾化−明矾石化−高岭土化−石膏和中性斑岩体的绢云母化−硅化为显著特征, 远矿蚀变为青磐岩化(许英霞等, 2008)。矿区中心为泥化、硫酸盐带, 向北为硅化−伊利石化−绢云母化−褐铁矿化带。具有找矿标志的矿化蚀变为:硅化−绢云母化、伊利石化、褐铁矿化、硅化、高岭土化−明矾石化等。
红山铜金矿床矿化以铜、金矿为主, 赋矿岩性主要为流纹岩、英安岩、安山玢岩、闪长玢岩、火山角砾岩。矿石矿物主要为黄铜矿、蓝铜矿、闪锌矿、磁铁矿、赤铁矿、黄钾铁钒等。脉石矿物为石英、斜长石、绢云母、绿泥石等。矿石呈稠密浸染状、细脉状、角砾状构造(许英霞, 2010)。
7.2 红山铜金矿床高光谱遥感地空综合找矿模型的建立
充分研究区域成矿有利的构造、地层、岩浆岩条件、变质条件的基础上, 结合区域控矿因素、矿床、矿化蚀变特征, 进行航空高光谱蚀变矿物分布特征分析及评价。本文截取了研究区内红山铜金矿床及外围基于HyMap航空高光谱遥感影像的蚀变矿物分布图(图 6), 图中显示红山铜金矿床中心位置处于开采状态, 呈黑白相间且颜色突变的色调, 影纹粗糙且不均匀, 矿床外围呈灰黑色、灰绿色色调, 影纹粗糙。与矿床有关的航空高光谱蚀变矿物组合主要为褐铁矿+黄钾铁钒、绢云母+蒙脱石等, 呈大面积的团块状分布。在矿床外围分布着呈星散状、星点状的绿泥石+绿帘石蚀变矿物组合, 该组蚀变为铜金矿床外围主要发育的蚀变组合, 对于寻找该类型矿床有一定的指示作用。
在红山铜金矿区收集了具体的剖线采样和光谱数据, 图 7为红山铜金矿区的地质光谱综合剖面, 出露的地层为中奥陶统荒草坡群大柳沟组, 岩性主要为英安岩、片理化带、流纹斑岩、花岗岩、花岗闪长岩、构造蚀变岩以及铜金矿化体等。蚀变矿物主要有褐铁矿、黄钾铁钒、绢云母、石膏、高岭石、明矾石及方解石等。通过样品采集和光谱测试结果, 研究与该区铜金矿关系密切的几种赋矿岩性的光谱曲线, 分别为构造蚀变岩、花岗闪长岩、花岗岩、流纹斑岩、蚀变英安岩、及铜金矿化体。
图 7 红山铜金矿区剖线地面高光谱蚀变矿物分布图1—英安岩; 2—片理化带; 3—流纹斑岩; 4—花岗岩; 5—花岗闪长岩; 6—断裂构造; 7—地面光谱样采样点; 8—地面光谱测试点编号; 9—褐铁矿; 10—黄钾铁矾; 11—绢云母; 12—石膏; 13—高岭石; 14—明矾石; 15—方解石Figure 7. Surface hyperspectral alteration mineral distribution along geological section in the Hongshan copper gold orefield1-Dacite; 2-Schistosity zone; 3-Rhyolitic porphyry; 4-Granite; 5-Granite diorite; 6-Fracture structure; 7-Surface spectrum sampling site; 8-Serial number of spectrum test point; 9-Limonite; 10-Jarosite; 11-Sericite; 12-Gypsum; 13-Kaolinite; 14-Alunite; 15-Calcite经研究分析发现:红山铜金矿床及外围不同岩性地面光谱特征各不相同, 各类岩性具体曲线(图 8)特征如下:
图 8 红山铜金矿床及外围各岩性光谱曲线对比图a—构造蚀变岩; b—花岗闪长岩; c—花岗岩; d—流纹斑岩; e—蚀变英安岩; f—铜金矿化体Figure 8. Comparison diagram of various lithology spectral curves from the Hongshan copper gold deposit and peripheral areasa-Tectonic altered rock; b-Granodiorite; c-Granite; d-Rhyolitic porphyry; e-Altered dacite; f-Copper gold mineralization body构造蚀变岩(图 8-a)地面光谱曲线特征显示: 880 nm处的Fe3+特征吸收可能为褐铁矿化所致; 2210 nm处的弱吸收峰为发育弱绢云母化所致, 2260 nm处的较强吸收峰为岩石发育黄钾铁矾化所致。故构造蚀变岩蚀变矿物组合为:褐铁矿+黄钾铁矾+绢云母。
花岗闪长岩(图 8-b)地面光谱曲线特征显示: 880 nm处较强的对称吸收峰为岩石发育褐铁矿化所致; 2210 nm处较强的吸收峰, 及右肩2250 nm处较弱次级吸收峰, 是Al−OH的特征吸收, 为明矾石化所致; 2350 nm处强的主吸收峰, 左肩2250 nm处较强的次级吸收峰, 是Mg−OH的特征吸收, 为绿泥石化所致。故其蚀变矿物组合为:褐铁矿+明矾石+绿泥石。
花岗岩(图 8-c)地面光谱曲线特征显示: 910 nm处较弱的Fe3+特征吸收峰为褐铁矿化所致; 2210 nm处的弱Al−OH特征吸收峰为绢云母化引起, 其右肩2250 nm处的较弱吸收峰可能为明矾石化所致; 2265 nm处的较弱吸收峰可能为黄钾铁矾化所致; 2330 nm处可见较强吸收峰, 是CO32-的特征吸收峰, 为岩石发育碳酸盐化所致, 故花岗岩蚀变矿物组合为:褐铁矿+黄钾铁矾+绢云母+明矾石+方解石。
流纹斑岩(图 8-d)地面光谱曲线特征显示: 880 nm处的较强吸收峰为褐铁矿化所致; 2210 nm处的强Al−OH特征吸收峰为岩石发育绢云母化所致; 2265 nm处的弱Fe−OH特征吸收峰为岩石发育黄钾铁矾化所致, 故流纹斑岩蚀变矿物组合为:褐铁矿+黄钾铁钒+绢云母。
蚀变英安岩(图 8-e)地面光谱曲线特征显示: 880 nm处的弱吸收峰为弱褐铁矿化所致; 2210 nm处的弱Al−OH特征吸收峰为绢云母化引起; 2350 nm处的较强主吸收峰及其左肩2250 nm处的较强次级吸收峰, 是Mg−OH的特征吸收峰, 为绿泥石化所致; 2330 nm处的较强主吸收峰以及1880 nm、2060 nm、2170 nm处较弱次级吸收峰, 是CO32-的特征吸收峰, 为岩石发育碳酸盐化所致; 1836 nm、1935 nm、2213 nm、2270 nm处可见多处较强吸收峰, 为岩石发育石膏所致。故蚀变英安岩蚀变矿物组合为:褐铁矿+绢云母+绿泥石+方解石+石膏。
铜金矿化体(图 8-f)地面光谱曲线特征显示: 880 nm处的较强Fe3+特征吸收为岩石发育褐铁矿化所致; 2210 nm处的弱Al−OH特征吸收峰为岩石发育绢云母化; 2265 nm处的弱Fe−OH特征吸收峰为岩石发育黄钾铁矾化; 1540 nm、1835 nm、1920 nm、2000 nm、2165 nm、2208 nm、2250 nm、2329 nm、2374 nm、2385 nm处可见多处较强吸收峰, 可能为岩石发育黝帘石化所致。故铜金矿化体蚀变矿物组合为:褐铁矿+绢云母+黄钾铁钒+黝帘石。
根据航空高光谱和地面高光谱蚀变矿物或组合的类型和分布特征, 结合红山铜金矿床的地质特征, 进行典型矿床航空−地面高光谱综合剖析, 最终建立了红山铜金矿床高光谱遥感地空综合找矿模型(表 3)。
表 3 红山铜金矿床高光谱遥感地空综合找矿模型Table 3. Hyperspectral remote sensing air comprehensive prospecting model for the Hongshan copper gold deposit7.3 铜金矿床找矿有利地段优选
在建立红山铜金矿床高光谱遥感地空综合找矿模型的基础上, 结合地质、野外调查验证等多源地学信息, 运用高光谱遥感地空综合预测方法, 优选了2处与已知矿床航空−地面高光谱向类似的找矿有利地段, 分别为I−1和I−2(图 9)。
图 9 新疆卡拉塔格地区找矿有利地段优选图1—全新统:砂石及砾石沉积; 2—下侏罗统三工河组:粉砂质泥岩、泥质砂岩夹钙质砂岩; 3—上二叠统库莱组:砾岩、中粗粒岩屑砂岩; 4—中二叠统阿尔巴萨依组:凝灰岩、橄榄玄武岩; 5—上石炭统脐山组:玄武岩、安山岩; 6—下泥盆统大南湖组:安山岩、英安岩、熔结凝灰岩、火山角砾岩; 7—中-上志留统红柳峡组:凝灰质砾岩、晶屑岩屑凝灰岩; 8—中奥陶统荒草坡群大柳沟组:玄武岩、安山岩、英安岩; 9—志留纪二长花岗岩; 10—志留纪花岗闪长岩; 11—志留纪英云闪长岩; 12—花岗岩体; 13—地质界线; 14—角度不整合界线; 15—岩相界线; 16—逆断层; 17—性质不明断层; 18—推测断层; 19—铜金矿; 20—铜锌矿; 21—铜锌金矿; 22—铜铁矿; 23—高光谱找矿有利地段Figure 9. Favorable ore-prospecting optimization segments in Kalatag area of Xinjiang1-Holocene, sand and gravel deposits; 2-Lower Jurassic Sangonghe Formation: river silty mudstone, argillaceous sandstone intercalated with calcareous sandstone; 3-Permian Kulai Formation: breccia, coarse grained lithic sandstone; 4-Middle Permian Aebasayi Formation: tuff, olivine basalt; 5-Carboniferous Qishan Formation: basalt, andesite; 6-Lower Devonian Dananhu Formation: andesite, dacite, welded tuff and volcanic breccia; 7-Middle-upper Silurian Hongliuxia Formation: tuffaceous conglomerate, wafer crumbs debris tuff; 8-Daliugou Formation of Ordovician Huangcaopo Group: basalt, andesite, dacite; 9-Silurian monzonitic granite; 10-Silurian granite diorite; 11-Silurian tonalite; 12-Granite rock mass; 13-Geological boundary; 14-Angular unconformity; 15-Petrofacies boundary; 16-Thrust; 17-Unknown fault; 18-Inferred fault; 19-Copper and gold deposit; 20-Copper and zinc deposit; 21-Copper zinc deposit; 22-copper deposit; 23-Hyperspectral favorable ore-prospecting segmentsI−1找矿有利地段分布在玉带铜金矿点北西侧一带, 呈北北西向展布, 面积约5.2 km2。区内分布的地层为中−上志留统红柳峡组(S2-3h)火山岩, 近东西向断裂构造发育, 区内多发育褐铁矿化、硅化、绿帘石化等蚀变。微量元素分析23个, 化学分析显示钛含量最高达0.23×10-2, 金红石含量最高达0.38×10-2。为铜、金、钛等元素多金属找矿有利地段(图 10)。
I−2找矿有利地段分布在东二区南东侧, 呈近东西向展布, 面积约6.8 km2。区内分布的地层为中奥陶统荒草坡群大柳沟组(O2Hd)火山岩, 北东东向断裂构造发育, 花岗闪长岩体呈近东西向展布, 区内多发育硅化、褐铁矿化、赤铁矿化、绿帘石化、孔雀石化等蚀变。微量元素分析15个, 化学分析显示钛元素含量最高达0.21×10-2, 金红石含量最高达0.35×10-2, 且区内已有红山铜金矿床分布, 为铜、金、钛等元素多金属找矿有利地段(图 11)。
8. 结论
(1) 新疆东天山卡拉塔格铜金多金属矿床集中, 地表具有指示意义的蚀变矿物分布较为典型, 适合利用高光谱遥感地空综合预测方法在该区寻找铜金及多金属矿床, 故该方法为寻找铜金矿床的有效技术方法。
(2) 在运用HyMap航空成像光谱仪和FieldSpec Pro FR光谱仪进行高光谱数据获取和处理的基础上, 结合HyMap航空高光谱蚀变矿物提取效果分析, 最终开展典型矿床航空−地面高光谱综合剖析研究, 建立了红山铜金矿床高光谱遥感地空综合找矿模型, 并综合地质、野外调查验证等多源地学信息, 圈定了I−1和I−2等两处找矿有利地段。
(3) 与常规的地质、区域化探等地质矿产勘查方法相比, 高光谱遥感地空综合预测方法具有低成本、快速高效的优势。但由于各类矿床的找矿标志和高光谱异常信息多种多样, 从而造成蚀变异常信息的多解性, 加上利用高光谱地空综合预测方法开展多金属矿产资源调查研究还比较少, 需要今后继续完善成熟。
-
图 2 1957-2006年孝感市年降水量变化图(张涛等,2010)
Figure 2. Xiaogan's mean precipitation from 1957 to 2006 (Zhang et al., 2010)
表 1 试验场区地质结构信息表
Table 1 Geological structure information of test site
表 2 试验场区土壤岩性命名及垂向渗透系数
Table 2 Soil texture and vertical hydraulic conductivity of the test site
表 3 试验场区土壤体积含水率
Table 3 Soil volumetric moisture content (VWC) of test site
表 4 ZK1土壤水及地下水的δD、δ18O值表(2019/4/12)
Table 4 δD、δ18O of soil water in ZK1 and groundwater
表 5 ZK2土壤水及地下水的δD、δ18O值(2019/7/25)
Table 5 δD、δ18O of soil water in ZK2 and groundwater
表 6 大气降雨、土壤水及地下水的氢氧稳定同位素组成特征表
Table 6 Hydrogen and oxygen stable isotopes of soil moisture, precipitation and groundwater
表 7 孝感市1957-2018年丰水年、枯水年一览
Table 7 Low flow year and high flow year in Xiaogan from 1957 to 2018
表 8 年均入渗补给量
Table 8 Yearly average annual recharge
-
Cao Jianwen, Xia Riyuan, Fang Shangwu, Zhao Liangjie, Wang Zhe, Wang Ruofan, Yi Rui. 2019. Model and mechanism of "water exploration by cross layer" for high sulfate area in slope region of Yunnan-Guizhou Plateau[J]. Geology in China, 46(2): 235-243(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-DIZI201902003.htm
Craig H. 1961. Isotopic Variations in Meteoric Waters[J]. Science, 133(3465): 1702-1703. doi: 10.1126/science.133.3465.1702
Dai Junjie, Zhang Xinping, Luo Zidong, Wang Yue, Liu Fuji, He Xinguang. 2019. Characteristics of stable isotope in soil water and indication to soil water movement in a grove dominated by Cinnamomun camphora in Changsha[J]. Research of Environmental Sciences, 32(6): 974-983(in Chinese with English abstract). http://www.zhangqiaokeyan.com/academic-journal-cn_research-environmental-sciences_thesis/0201272816983.html
Deng Zhimin, Zhang Xiang, Pan Guoyan. 2016. Variations of hydrogen and oxygen isotopes in Meteoric precipitation in Wuhan, China[J]. Journal of Yangtze River Scientific Research Institute, 33(7): 12-17, 22(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-CJKB201607003.htm
Dong Xiaofang, Deng Huangyue, Zheng Xiangmin, Zhou Limin. 2017. Analysis of stable isotope characteristics and water vapor origins in atmospheric precipitation in the Yangtze River Basin[J]. Environmental Science & Technology, 40(4): 78-84(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-FJKS201704013.htm
Ministry of Water Resources of the People's Repulic of China. 2008. Code for Engineering Geological Investigation of Water Resources and Hydropower (GB 50487-2008)[S]. Beijing: China Planning Press(in Chinese).
He Jun, Tan Ting, Qi Zhichong, Hao Yiguo, Ma Teng. 2015. Study on contamination sources of the nitrate in groundwater by nitrogen and oxygen isotopes in Ramotswa, Botswana[J]. Geology and Mineral Resources of South China, 31(4): 421-427(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-HNKC201504012.htm
Hsieh J C, Chadwick O A, Kelly E F, Savin S M. 1998. Oxygen isotopic composition of soil water: Quantifying evaporation and transpiration[J]. Geoderma, 82(1/2/3): 269-293. http://www.researchgate.net/profile/Oliver_Chadwick/publication/229406618_Oxygen_isotope_composition_of_soil_water_Quantifying_evaporation_and_transpiration/links/00b49538df37c86d51000000.pdf
Hu Haiying, Bao Weimin, Wang Tao. 2008. Simulation and experiment on variations of hydrogen and oxygen isotopes in soil water[J]. Water Resources and Power, (4): 149-152(in Chinese with English abstract). http://www.onacademic.com/detail/journal_1000037231458410_f4de.html
Lei Zhidong, Yang Shixiu, Xie Senchuan. 1988. Soil Water Dynamics[M]. Beijing: Tsinghua University Press, 5-7(in Chinese).
Ji Wangjia, Huang Yanan, Li Bingbing, Li Zhi. 2019. Oxygen and hydrogen stable isotope compositions of soil water in deep loess profile under different land use types of northern Shaanxi, China[J]. Chinese Journal of Applied Ecology, 30(12): 4143-4149(in Chinese with English abstract). http://www.ncbi.nlm.nih.gov/pubmed/31840459
Liu Beilin, Phillips F, Hoines S, Andrew. R, Pankaj S. 1995. Water movement in desert soil traced by hydrogen and oxygen isotopes, chloride, and chloride-36, southern Arizona[J]. Journal of Hydrology, 168: 91-110. doi: 10.1016/0022-1694(94)02646-S
Liu Jun, Nie Zhenlong, Duan Baoqian, Liu Fuliang, Zhang Lin. 2016. Characteristics of stable isotope (δ2H、δ18O) in soil water in Hohhot area[J]. Journal of Arid Land Resources and Environment, 30(10): 145-150(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-GHZH201610025.htm
Liu Jun, Wei Wen, Zhang Lin, Wang Ying, Duan Baoqian, Liu Fuliang. 2012. Application on isotopes D and 18O of soil water in water movement of unsaturated zone[J]. Site Investigation Science and Technology, (5): 38-43(in Chinese with English abstract). http://www.en.cnki.com.cn/Article_en/CJFDTOTAL-KCKX201205011.htm
Ma Bin, Liang Xing, Lin Dan, Liu Shaohua, Xu Min. 2014. Tracing infiltration and recharge of the unsaturated zone using 2H, 18O Isotopes in Shijiazhuang, North China Plain[J]. Geological Science and Technology Information, 33(3): 163-168, 174(in Chinese with English abstract). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dzkjqb201403021
Ma Yingbin, Zhang Biebei, Xu Qing, Gao Deqiang, Wang Ting, Sui Mingzhen, Huang Yaru. 2018. Response of deuterium isotope in soil water to rainfall in freshwater wetland forests of Shaoxing, Zhejiang Province[J]. Scientia Silvae Sinicae, 54(10): 11-19(in Chinese with English abstract).
Prudhomme Edmunds W M, Ma Jinzhu, Li Ding, Zhang Jiawu. 2003. Groundwater recharge and climatic change during the last 1000 years from unsaturated zone of SE Badain Jaran Desert[J]. Chinese Science Bulletin, 48(14): 1469-1474. doi: 10.1360/02wd0262
Song Xianfang, Wang Shiqin, Xiao Guoqiang, Wang Zhimin, Liu Xin, Wang Peng. 2009. A study of soil water movement combining soil water potential with stable isotopes at two sites of shallow groundwater areas in the North China Plain[J]. Hydrological Process, 23(9): 1376-1388. doi: 10.1002/hyp.7267
Sun Fangqiang, Yin Lihe, Wang Xiaoyong, Ma Hongyun, Zhang Jun, Dong Jiaqiu, He Shuaijun. 2017. Determination of vertical infiltration recharge of groundwater in the thick unsaturated zone of Sangong River Basin, Xinjiang[J]. Geology in China, 44(5): 913-923(in Chinese with English abstract).
Wang Fugang, Liao Zisheng. 2007. Study of the precipitation infiltration recharge with the D、18O isotopes peak displacement Method[J]. Journal of Jilin University (Earth Science Edition), (2): 284-287, 334(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-CCDZ200702013.htm
Wu Qianqian, Ren Jiaguo, Xu Mo. 2008. Isotope features and supply sources of groundwater in the Yarkant River drainage area, Xinjiang[J]. Geology in China, (2): 331-336(in Chinese with English abstract).
Xia Pingying, Zhang Feng. 2014. Study on flood control and disaster reduction in Xiaogan City[J]. Science and Technology Innovation Herald, 11(15): 21-22(in Chinese).
Xu Yingde, Wang Jingkuan, Gao Xiaodan, Zhang Yulong. 2018. Application of hydrogen and oxygen stable isotope techniques on soil water research: A review[J]. Journal of Soil and Water Conservation, 32(3): 1-9, 15(in Chinese with English abstract). http://www.researchgate.net/publication/325965897_Application_of_Hydrogen_and_Oxygen_Stable_Isotope_Techniques_on_Soil_Water_Reaearch_A_Review
Yu Hongyang, 2016. Xiaogan 2016: Super historical flood attacks[J]. Insight China, (19): 18-23(in Chinese).
Zhang Tao, Duan Chunfang, Fang Fang, Wang Li. 2010. Precipitation variation characteristics of Xiaogan city for last 50 years[J]. Torrential Rain and Disasters, 29(1): 81-84(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-HBQX201001015.htm
Zhang Xiang, Deng Zhimin, Pan Guoyan, Wu Shaofei, Zhu Yang, Zhu Cairong. 2015. Variation in stable composition in soil water in Poyang Lake Wetland[J]. Acta Ecologica Sinica, 35(22): 7580-7588(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTotal-STXB201522032.htm
Zhang Xiaojuan, Song Weifeng, Wu Jinkui, Wang Zhuojuan. 2015. Characteristics of Hydrogen and Oxygen Isotopes of Soil Water in the Water Source Area of Yuanyang Terrace[J]. Environmental Science, 36(6): 2102-2108(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-HJKZ201506029.htm
Zimmermann U, Munnich K O, Roether W. 1967. Dowanward movement of soil moisture traced by means of hydrogen isotopes[J]. Isotope Techniques in the Hydrologic Cycle, 11: 28-36.
曹建文, 夏日元, 方尚武, 赵良杰, 王喆, 王若凡, 易瑞. 2019. 云贵高原斜坡地带典型地下水富硫酸盐地区"越层找水"模式及其机理研究[J]. 中国地质, 46(2): 235-243. http://geochina.cgs.gov.cn/geochina/ch/reader/view_abstract.aspx?file_no=20190202&flag=1 戴军杰, 章新平, 罗紫东, 王锐, 刘福基, 贺新光. 2019. 长沙地区樟树林土壤水稳定同位素特征及其对土壤水分运动的指示[J]. 环境科学研究, 32(6): 974-983. https://www.cnki.com.cn/Article/CJFDTOTAL-HJKX201906007.htm 邓志民, 张翔, 潘国艳. 2016. 武汉市大气降雨的氢氧同位素变化特征[J]. 长江科学院院报, 33(7): 12-17, 22. https://www.cnki.com.cn/Article/CJFDTOTAL-CJKB201607003.htm 董小芳, 邓黄月, 郑祥民, 周立旻. 2017. 长江流域降水中氢氧同位素特征及水汽来源[J]. 环境科学与技术, 40(4): 78-84. https://www.cnki.com.cn/Article/CJFDTOTAL-FJKS201704013.htm 中华人民共和国水利部. 2009. 水利水电工程地质勘察规范(GB 50487-2008)[S]. 北京: 中国计划出版社. 何军, 谭婷, 祁志冲, 郝义国, 马腾. 2015. 利用15N和18O识别地下水硝酸盐污染源-以博茨瓦纳东南部拉莫茨瓦为例[J]. 华南地质与矿产, 31(4): 421-427. doi: 10.3969/j.issn.1007-3701.2015.04.011 胡海英, 包为民, 王涛, 瞿思敏. 2008. 土壤水中氢氧同位素变化模拟及实验[J]. 水电能源科学, (4): 149-152. doi: 10.3969/j.issn.1000-7709.2008.04.046 姬王佳, 黄亚楠, 李冰冰, 李志. 2019. 陕北黄土区深剖面不同土地利用方式下土壤水氢氧稳定同位素特征[J]. 应用生态学报, 30(12): 4143-4149. https://www.cnki.com.cn/Article/CJFDTOTAL-YYSB201912017.htm 雷志栋, 杨诗秀, 谢森传. 1988. 土壤水动力学[M]. 北京: 清华大学出版社, 5-7. 刘君, 聂振龙, 段宝谦, 田言亮, 刘福亮, 张琳. 2016. 氢氧稳定同位素指示的呼和浩特地区土壤水的补给特征[J]. 干旱区资源与环境, 30(10): 145-150. https://www.cnki.com.cn/Article/CJFDTOTAL-GHZH201610025.htm 刘君, 卫文, 张琳, 王莹, 段宝谦, 刘福亮. 2012. 土壤水D和18O同位素在揭示包气带水分运移中的应用[J]. 勘察科学技术, (5): 38-43. doi: 10.3969/j.issn.1001-3946.2012.05.010 马斌, 梁杏, 林丹, 刘绍华, 徐敏. 2014. 应用2H、18O同位素示踪华北平原石家庄包气带土壤水入渗补给及年补给量确定[J]. 地质科技情报, 33(3): 163-168, 174. https://www.cnki.com.cn/Article/CJFDTOTAL-DZKQ201403023.htm 马迎宾, 张蓓蓓, 徐庆, 高德强, 王婷, 隋明浈, 黄雅茹. 2018. 绍兴淡水湿地森林土壤水氘同位素对降水的响应[J]. 林业科学, 54(10): 11-19. doi: 10.11707/j.1001-7488.20181002 孙芳强, 尹立河, 王晓勇, 马洪云, 张俊, 董佳秋, 贺帅军. 2017. 新疆三工河流域厚层包气带区地下水垂向补给量的厘定[J]. 中国地质, 44(5): 913-923. http://geochina.cgs.gov.cn/geochina/ch/reader/view_abstract.aspx?file_no=20170506&flag=1 王福刚, 廖资生. 2007. 应用D、18O同位素峰值位移法求解大气降雨入渗补给量[J]. 吉林大学学报(地球科学版), (2): 284-287+334. https://www.cnki.com.cn/Article/CJFDTOTAL-CCDZ200702013.htm 武倩倩, 任加国, 许模. 2008. 新疆叶尔羌河流域地下水同位素特征及其补给来源分析[J]. 中国地质, (2): 331-336 doi: 10.3969/j.issn.1000-3657.2008.02.018 夏苹应, 张峰. 2014. 孝感城市防洪与减灾对策研究[J]. 科技创新导报, 11(15): 21-22. doi: 10.3969/j.issn.1674-098X.2014.15.023 徐英德, 汪景宽, 高晓丹, 张玉龙. 2018. 氢氧稳定同位素技术在土壤水研究上的应用进展[J]. 水土保持学报, 32(3): 1-9, 15. https://www.cnki.com.cn/Article/CJFDTOTAL-TRQS201803001.htm 余弘阳. 2016. 孝感2016: 超历史洪水来袭[J]. 小康, (19): 18-23. https://www.cnki.com.cn/Article/CJFDTOTAL-CHXK201619007.htm 张涛, 段春锋, 方芳, 王俐. 2010. 近50年孝感市降水变化特征分析[J]. 暴雨灾害, 29(1): 81-84. doi: 10.3969/j.issn.1004-9045.2010.01.014 张翔, 邓志民, 潘国艳, 吴绍飞, 肖洋, 朱才荣. 2015. 鄱阳湖湿地土壤水稳定同位素变化特征[J]. 生态学报, 35(22): 7580-7588. https://www.cnki.com.cn/Article/CJFDTOTAL-STXB201522032.htm 张小娟, 宋维峰, 吴锦奎, 王卓娟. 2015. 元阳梯田水源区土壤水氢氧同位素特征[J]. 环境科学, 36(6): 2102-2108. https://www.cnki.com.cn/Article/CJFDTOTAL-HJKZ201506029.htm -
期刊类型引用(13)
1. 马庆,王元伟,高永宝,刘明. 青海省玛多县夺儿贡玛地区地球化学特征与找矿方向. 西北地质. 2025(01): 178-185 . 百度学术
2. 杨海涛,刘新伟,汪超,牛亮,胡西顺,门文辉,杨文刚. 蠎西寺沟斑岩-矽卡岩型钨钼矿物化探异常特征及找矿模型. 地质与勘探. 2022(05): 929-939 . 百度学术
3. 白德胜,李水平,纵瑞,程华,齐勇攀,张爱玲,孙进,赵华奇. 豫西董家埝构造蚀变岩型银矿物化探异常特征及找矿模型. 地质与勘探. 2021(02): 241-253 . 百度学术
4. 刘亮,王超,孙占营,张杰,马运超,赵相国,杨宇东. 青海多尔娘地区地层和岩体含矿性研究. 河南理工大学学报(自然科学版). 2021(06): 64-76 . 百度学术
5. 王占彬,宋贺民,马庆,谢志远,李保飞,王凯. 河北省怀安县朱家洼矿区地球化学特征与找矿方向. 地质与勘探. 2020(01): 102-112 . 百度学术
6. 刘亮,王超,孙占营,余长荣,赵相国. 综合物化探方法在青海珊旗根玛金锑矿点的运用. 科学技术与工程. 2020(04): 1337-1343 . 百度学术
7. 李菲. 蟒西地区中酸性岩体地质特征及找矿潜力. 科学技术与工程. 2020(18): 7176-7182 . 百度学术
8. 丁吉顺,陈伟,周恒,郭奇奇,孙渺,张祎. 西藏雄梅地区1:5万水系沉积物地球化学特征及找矿远景. 地质与勘探. 2019(01): 48-63 . 百度学术
9. 张翔,戴霜,黄万堂,赵振斌,李鸿睿,王玉玺,刘博,吴茂先. 甘肃省玛曲县大水金矿原生金矿石的发现及意义. 地质与勘探. 2019(02): 484-495 . 百度学术
10. 段吉学,刘江. 综合物化探在内蒙萤石多金属矿普查中的应用研究. 西北地质. 2019(03): 265-274 . 百度学术
11. 王伟. 新疆托克逊县阿热塔格山Ⅱ区的物化探异常特征. 黑龙江科技大学学报. 2019(05): 546-551 . 百度学术
12. 缪宇,宋文婷,何茂源,徐乐. 云南弥渡县云景地区多元地学信息集成及找矿模型. 地质与勘探. 2019(06): 1367-1378 . 百度学术
13. Shi-hong Zhang,Ke-yan Xiao,Jian-ping Chen,Jie Xiang,Ning Cui,Xiao-nan Wang. Development and future prospects of quantitative mineral assessment in China. China Geology. 2019(02): 198-210 . 必应学术
其他类型引用(0)