浅覆盖区隐伏萤石矿找矿预测:以内蒙古赤峰俄力木台为例
唐利1,张寿庭1,*,王亮1,裴秋明2,方乙3,曹华文4,邹灏5,尹少波6
1中国地质大学(北京)地球科学与资源学院
2西南交通大学地球科学与环境工程学院
3西南石油大学地球科学与技术学院
4中国地质调查局成都地质调查中心
5成都理工大学地球科学学院
6淮矿西部煤矿投资管理有限责任公司
作者简介:唐利,博士,副教授,主要从事固体矿产资源勘查评价方面的教学和研究工作。
通信作者简介:张寿庭,博士,教授,博士生导师,主要从事固体矿产资源勘查评价方面的教学和研究工作。
导读:
内蒙古赤峰俄力木台地区第四系覆盖严重,是开展浅覆盖区隐伏矿体找矿预测方法实践的理想地区。
萤石矿属于我国战略性矿产,需加强勘查找矿。
萤石矿为断裂控制的热液充填型脉状矿,规模小,且密度、磁性、电性等物性特征与围岩差异不大,寻找覆盖层下隐伏矿难度较大。
作者利用地质特征-地球物理-地球化学-遥感解译(地-物-化-遥)综合找矿方法进行浅覆盖区萤石矿体的定位预测研究,钻探验证找矿效果较好,值得读者阅读!
关键词:地-物-化-遥找矿方法;成矿预测;隐伏萤石矿。
内容提纲:
0 前言
1 区域地质背景
2 俄力木台萤石矿化地质特征
3 遥感-物探-化探找矿方法实践
4 讨论
5 结论
0 前言
萤石作为重要的矿物原料,广泛应用于冶金、制冷、建材、光学工业、新能源、新材料等现代工业领域。萤石是美国、中国、欧盟和日本等国的关键矿产,具有战略资源供应安全等重要意义。中国是萤石资源大国,根据美国地调局最新统计数据,我国萤石基础储量仅次于南非,居世界第2位。我国的萤石矿具有品质较高、开采条件较好等优势,在全球萤石行业中占有重要地位。然而,我国萤石产量长期居世界第1位,约占世界总产量的1/2,过度开采导致我国的萤石资源优势逐渐削弱,现面临严峻的资源形势。随着地表和浅部萤石资源逐渐枯竭,覆盖区隐伏萤石矿体及深部找矿将成为未来萤石接替资源的首要目标。
矿体深部定位预测是矿产勘查领域的主要难题和研究热点之一。覆盖区深部找矿主要面临深度大、干扰强、精度低、深部成矿规律难以总结等挑战。在大多数情况下,获取深部地质地球化学信息较为困难,尤其是热液脉型萤石矿体。针对热液脉型萤石矿床的找矿方法及深部成矿预测,勘查地球化学方法主要涉及原生晕圈定和水系沉积物测量。原生晕和次生晕中的深部矿化信息和元素分带异常容易受到围岩或覆盖沉积物中元素特征、后期风化和外来运积物的影响。地球物理方法主要包括甚低频电磁法和双频激发极化法,但容易受到各种因素干扰,其结果具有多解性。由于地质特征的复杂性,依靠单一物探或化探方法揭示矿致异常难以在矿产资源预测中得到可靠的结论。随着萤石找矿难度的不断加大,以往单一的找矿技术已不能满足隐伏萤石找矿的要求。针对热液脉型萤石矿床,利用地-物-化-遥综合找矿方法进行矿体定位预测研究,对寻找隐伏矿床具有重要意义。
内蒙古赤峰地区位处环太平洋成矿域和古亚洲洋成矿域的叠加复合部位,因多期构造-岩浆-热液活动叠加,成矿地质条件优越,形成了一系列受断裂带控制的热液脉型萤石矿床。俄力木台地区是林西水头萤石成矿带的南段延伸,成矿和找矿潜力巨大。俄力木台地区第四系覆盖严重,是评价不同找矿技术方法并寻找浅覆盖区隐伏萤石矿体的理想地区。本文利用野外地质调查、甚低频电磁法、便携式X射线荧光分析、偏提取地球化学和WorldView-2高分辨率多光谱遥感等地-物-化-遥综合找矿方法,基于萤石矿体的垂向分带规律、遥感异常信息解译和物化探矿致异常信息提取,对俄力木台地区开展隐伏萤石矿体定位预测,为浅覆盖区受断裂控制的脉状萤石矿体的找矿勘查提供借鉴。
1 区域地质背景
大兴安岭南段位于二连—赫根山断裂以南,西拉木伦断裂以北,大地构造位置处于中亚造山带东段(图1a)。中亚造山带位于华北板块和西伯利亚板块之间,经历晚古生代至早中生代多期次板块拼合和构造运动,并伴有强烈的岩浆活动,为区内多种类型矿床的形成提供了优越的成矿条件。研究区同时位于环太平洋成矿带与古亚洲洋成矿带的交合部位,受古亚洲洋成矿域和环太平洋成矿域的叠加影响,区内发育一系列与构造-岩浆-热液活动相关的斑岩型、夕卡岩型、隐爆角砾岩型、热液脉型有色金属-多金属矿床和非金属矿床。
区域出露地层主要有前寒武系变质结晶基底,下古生界西别河组,二叠系寿山沟组、大石寨组、哲斯组和林西组,侏罗系新民组、满克头鄂博组和玛尼吐组,白垩系白音高老组和梅勒图组等。二叠系—侏罗系沉积-火山碎屑岩是区域萤石矿床的主要赋矿围岩,该地层主体呈北北东向展布,普遍经受浅变质作用。二叠系地层主要以黑色板岩、粉砂岩和砂岩为主,侏罗系地层以长英质火山碎屑岩、凝灰岩和流纹岩为主区内断裂十分发育,受古亚洲洋构造域、蒙古—鄂霍茨克构造域和环太平洋构造域的复合影响,形成了区域复杂的构造展布,其中EW向断裂带控制本区的主要构造格局,如西拉木伦断裂和索伦一林西断裂。NE-NNE向断裂带(如赫根山断裂和嫩江断裂)为区内最重要的区域性深大断裂,其次级断裂为重要的控岩-控矿断裂带,控制着区内大多数地质体及矿床的空间展布。区内中生代岩浆活动频繁,包括早侏罗世—早白垩世大面积的火山喷发和大规模的岩浆侵入,形成一系列的中酸性侵入岩和碱性岩,与区内的多金属-萤石矿床具有密切的时空联系。
2 俄力木台萤石矿化地质特征
水头萤石成矿带位于赤峰市林西县西北部。区内出露地层主要为二叠系寿山沟组、大石寨组、哲斯组和林西组,主要由黑色板岩、粉砂岩、砂岩、海相中酸性熔岩及凝灰岩组成,为一套海相火山-沉积碎屑岩系。水头萤石成矿带内目前已发现近20条萤石矿体,矿体严格受SN向、NNE向断裂控制,呈舒缓波状断续分布达10km。自北向南,水头萤石成矿带可划分为北段、中段和南段(图1b)。北段以赛波罗沟萤石矿床为代表,萤石矿体断续延伸2.5km;中段位于水头矿区,萤石矿体长约1.5km;南段为旱泡子一俄力木台一带,地表萤石矿化带断续延伸约4km。
区内萤石矿体主要呈透镜状和脉状赋存于断裂中,矿体产状与控矿断裂一致,总体走向NW345°至NNE30°,倾角65°〜85°。单个矿体厚度多为0.5〜8m,已控矿体延深可达250m以上。萤石呈紫色、浅绿色、白色,透明-半透明。萤石结构以粗晶自形-半自形粒状结构、细晶半自形-它形粒状结构为主,矿石构造主要有块状构造、条带状构造、角砾状构造和网脉状构造。根据野外地质特征和矿物共生次序,该区萤石成矿作用可划分为成矿早阶段(萤石+石英阶段)、成矿主阶段(石英+萤石阶段)和成矿晚阶段(碳酸盐阶段)3个阶段。矿石矿物为萤石,脉石矿物以石英为主,其次为蛋白石、方解石、高岭石、黄铁矿、重晶石、绿泥石和绢云母等。蚀变类型主要有硅化、绢云母化、黄铁矿化、碳酸盐化、高岭石化、绿泥石化和重晶石化等。
俄力木台地区属于水头萤石成矿带的南段延伸(图1b)。萤石矿化和热液蚀变在地表断续分布,表现为强硅化角砾岩、次生石英岩和石英网脉,萤石矿化少量发育。该区第四系覆盖严重,前期勘探程度较低,矿体展布情况不明(图2)。基于详细的野外地质调查,区内发现多处萤石矿化和蚀变露头,由北往南5处典型露头的矿化地质特征如下。
图2 俄力木台地质概况及勘探线布置图
(1)1号矿化露头。该点位于山脊鞍部,地表产出强硅化角砾岩、次生石英岩和石英网脉,主矿化蚀变带宽8〜10m(图3a)。矿化带受断裂控制,倾向为266°〜270°,倾角67°〜75°。矿物组合以石英为主,含少量蛋白石和萤石,见少量褐铁矿(由黄铁矿次生氧化)。萤石矿化主要呈团块状,少量发育细脉状(图3b)、次生淋滤孔洞状构造。
a,b—1号矿化露头;c,d—2号矿化露头;e,f一3号矿化露头;g,h—4号矿化露头;i—5号矿化露头。
图3 俄力木台矿化地质特征
(2)2号矿化露头。位于1号露头南侧小山梁,南沿山坡有一探槽揭露。地表特征与1号矿化露头相似,以产出强硅化角砾岩、次生石英岩和石英网脉为特征,主矿化蚀变带宽大于5m(图3c)。矿物组合以石英为主,见少量蛋白石、萤石。萤石矿化主要呈团块状,少量发育细脉状(图3d)。
(3)3号矿化露头。位于2号露头南侧山谷谷底,有一小型坑探工程揭露。与1、2号矿化露头相比,该点位于相对低洼处,高差30〜50m。该露头除发育强硅化和石英网脉外,有萤石粗脉体产出,萤石矿脉发育于矿化带近底板一侧(图3e)。矿化蚀变带宽达3〜5m,主脉带宽约1.5〜2m,萤石矿脉宽度20〜30cm。矿物组合以石英为主、萤石次之(图3f)。萤石呈脉状、小扁豆状产出,矿石构造以块状、正条带状和角砾状为主。
(4)4号矿化露头。位于3号露头南侧山坡,地表产出强硅化角砾岩、次生石英岩和石英网脉,矿化蚀变主脉带宽大于3m。矿物组成以石英为主,局部见少量萤石,萤石矿化弱(图3g,h)。
(5)5号矿化露头。位于4号露头南侧山坡,地表以强硅化角砾岩带和石英脉带产出(图3i),矿化蚀变带和石英脉宽度均比前者小,推断其属于矿化蚀变带的走向狭缩地段。
3 遥感-物探-化探找矿方法实践
3.1 WorldView-2高分辨率多光谱遥感
WorldView-2高分辨率多光谱遥感技术可以提供0.5m高空间分辨率的多光谱影像和0.46m分辨率的全色影像,能够精确展现地表微观地质信息(如地层、构造、矿体等)。在一定地质条件下,WorldView-2遥感影像能够将含矿地质体以其线、环、色、带、块等要素特征,更为清晰地反映于地表。在萤石矿区,控矿断裂及矿化蚀变带在WorldView-2遥感影像上一般表现出明显的线性影像特征。地表出露的硅质顶盖或者石英脉的抗风化能力较强,因此在浅覆盖区易形成断续分布的脊状地貌特征,具体在遥感影像图上表现为坑洼状和颗粒状。俄力木台及其周边区域的WorldView-2遥感影像清楚地揭示了该区域的已知矿带(脉)的发育部位与分布特征,并对该区的基本构造格局以及含矿断裂带的空间展布与变化特征也有较明显的反映(图4)。
图4 俄力木台地区WorldView-2遥感影像图
3.2甚低频电磁测量
甚低频电磁法利用不同地质体(岩体、矿体、构造等)之间的电性差异来定位地质异常,可查明隐伏-半隐伏含矿断裂带的矿化强度、规模和连续性,揭示低阻异常地质体的展布特征。林西地区不同类型萤石矿矿石和围岩的磁性、电性、密度等物性参数见表1。针对受断裂控制的热液脉型萤石矿床,利用甚低频电磁法能够清晰地展示萤石矿区低阻带的浅层分布特征,并进一步评价异常带在地表的延伸情况。本次研究利用重庆地质仪器厂生产的DDS-I型甚低频电磁仪,针对测区第四系覆盖严重,以及开展隐伏矿体的探测与空间定位的研究目的,采用磁倾角(D)法。测线方向与矿区含矿构造带方向近于垂直,剖面测量点距10m,系统扫面线距100m。
表1 林西地区岩矿石物性特征
甚低频电磁测量的低阻异常平面分布图显示,本区构造格局主要发育NNE-SN向、NW-NNW向两组构造,前者属本区主要的控矿构造,后者应属成矿后断裂并对前者的连续性有一定的破坏作用(图5)。甚低频电磁测量对测区北部的已知萤石矿带具有明显的反映,在测区东部发育有一条SN向贯通整个测区的较强的低阻异常带,其东侧还有一条次级异常带(图5)。测区断续岀露的矿化蚀变带与上述SN向的低阻异常带相吻合,强异常带较连续地分布于中段,延长达1km以上。
图5 俄力木台地区甚低频低阻异常平面分布图
3.3 便携式X射线荧光分析
便携式X射线荧光仪能够快速分析土壤或岩石中的元素含量,在矿产勘查领域应用广泛。在萤石找矿勘查中,Ca元素高值异常可以有效地指示萤石矿体的存在。本次研究对粒径>20目的原生土壤进行测量,采集腐殖层以下数十厘米深度(10〜50cm)的土壤,样品点距5〜10m。本次研究采用成都理工大学自主研发的便携式X射线荧光仪,型号为IED2000T。原生土壤样品压平或用压样器压实后进行测量,测量时间100s。
便携式X射线荧光分析元素含量异常剖面图中,10线Ca元素高值异常主要分布于矿化露头带并有向山坡下方迁移的趋势;20线Ca元素含量整体偏高,受地形影响,矿化露头带对应Ca元素高值异常向西偏移(图6,7)。上述剖面中的元素含量变化特征与萤石矿体位置有较好的对应关系,Ca元素高值异常能够指示萤石矿体的赋存部位。
图6 俄力木台10号勘探线物化探综合异常剖面图
图7 俄力木台20号勘探线物化探综合异常剖面图
3.4偏提取地球化学分析
土壤偏提取地球化学分析是针对覆盖层沉积物中残留或吸附态离子的提取和测量方法。偏提取地球化学法的野外工作测线方向、线距和点距与其他方法相同。样品为40cm左右深度的土壤,样品重量100g,过100目筛并粉碎至200目。准确称取2.5g样品,加入25mL盐酸羟胺-盐酸混合溶液。充分搅拌,静置24h,采用抽滤方法用5.0Mm的微孔滤出清液即得提取液,然后对提取液采用等离子质谱法分析F、Ca、Zn、Hg、Ag、Pb、Rb、Y、Ba、Sr、Be等元素,氟离子选择性电极测定F。
俄力木台萤石矿10线、20线偏提取地球化学结果显示(图6,7),F元素含量曲线整体平缓,在矿化体正上方以及NW侧(底板)有明显的异常,这些异常应属断裂与矿化异常的综合反应。
4 讨论
4.1 矿化分带特征与预测评价
研究矿化分带规律对隐伏-半隐伏矿体的找矿预测具有重要意义。萤石矿体空间分带(垂向、横向和纵向分带)是含矿热液沿特定负压构造破碎空间上升运移过程中,由于物理化学条件的递变,导致成矿组分在三维空间内按序晶出与沉淀富集,导致萤石矿体不同部位的矿石标型、元素组合及围岩蚀变呈现差异性变化和分带特征。徐旃章和张寿庭通过对浙江省萤石矿的综合研究,把受构造控制的萤石矿体在垂向上的分带特征总结为自上而下的硅质顶盖、头部矿体、中部矿体和尾部矿体四个部分(图8)。
内蒙古赤峰地区的典型萤石矿床具有类似的垂向分带特征。该区的硅质顶盖(1、2和4号矿化露头)以发育强硅化、次生石英岩和石英、蛋白石脉为特征,属萤石矿成矿早阶段的产物。硅质顶盖以硅质组分高,抗风化能力强,呈现脊骨状隆起、带状分布的正地形地貌为特征(图3)。在浅地表相对开放的断裂破碎带扩容空间,由于温度、压力等物化条件变化,以及成矿热液中SiO2、CaF2等不同组分的密度、结晶温度的差异,导致成矿热液运移前锋以SiO2为主,结晶分异的早期产物以石英或蛋白石为主。硅质顶盖属于主成矿期含矿热液沿控矿断裂带上升运移、充填的前锋产物。俄力木台地区的3号矿化露头为头部矿体,以发育石英网脉和萤石脉为特征。与硅质顶盖相比,头部矿体中萤石渐趋发育萤石多呈紫色。头部矿体多位于主干断裂带上部低级别、低序次构造的发育地段,由含矿热液充填形成的单脉体构成宽度小、密度大的脉体群。俄力木台北侧水头萤石矿床的开采深度达地下80〜120m开采对象主要为中部矿体。中部矿体为萤石主矿体的膨大部位,矿石品位较高,萤石含量可达50%以上。在控矿断裂的膨大部位,较大容量的成矿热液聚集导致温度降低缓慢,萤石结晶程度高。林西地区萤石矿的勘探开发程度总体较低,且尚未揭露到矿体尾部。总体而言,尾部矿体表现为萤石矿体的狭缩尖灭,成矿晚期碳酸盐阶段的产物如方解石脉以及相伴的矿化蚀变类型渐趋发育,尤以强碳酸盐化为特征。综合以上垂向矿化分带特征,俄力木台地区的硅质顶盖露头和头部矿体露头高差30〜50m由浅及深萤石矿化渐趋发育,深部可能发育萤石主矿体。
4.2 物-化-遥综合异常与预测评价
基于林西地区萤石矿找矿勘查结果,World-View-2遥感影像的萤石矿解译标志包括:(1)萤石矿化带和石英脉在遥感影像上呈现明显的线性分布特征;(2)地表出露的硅质顶盖或者石英脉的抗风化能力较强,因此在浅覆盖区易形成断续分布的脊状地貌特征,在遥感影像图上表现为坑洼状和颗粒状;(3)萤石矿带在遥感影像上颗粒感较强,石英脉带在局部表现为坑洼线性特征。俄力木台地区属于水头萤石成矿带的南延组成部分,由北向南,矿化带的走向方位变化明显,中-北段呈近SN向发育,南部至克什克腾旗境内主要呈NNE-NE向(图4)。本区成矿后NW-NNW向断裂发育,并表现为左行位移的错动特征,本区东部NNE-SN向山脉也有明显的节位错特征,由此推断它们对本区的萤石矿带也会有位错。上述构造在遥感影像上均呈现明显的线性分布特征,其中俄力木台北部的已知萤石矿带具有较强的颗粒感,中南部预测矿带的坑洼线性特征指示地表石英脉带的分布(图4)。
在萤石找矿勘查中,甚低频电磁法主要用于探测萤石矿体的控矿断裂或矿化蚀变带的空间展布,既可以进行重点地段单剖面测量,也可以开展大比例尺扫面,获取剖面沿线浅层电性结构以及低阻带的平-剖面空间展布。甚低频电磁测量在萤石矿控矿断裂带通常表现为较明显的低阻异常带。在俄力木台的勘探剖面沿线与已知矿化露头相对应的部位,甚低频电磁测量的低阻异常特征最为明显(图6,7)。甚低频电磁测量的低阻异常平面分布图中,俄力木台北部的已知萤石矿带具有较强的甚低频低阻异常特征(图5)。在测区东部发育有一条SN向的较强的甚低频低阻异常带,其东侧还有一条次级异常带,指示俄力木台中南段甚低频低阻异常带可能发育萤石矿化带。
便携式X射线荧光分析广泛应用于萤石矿找矿勘查,主要基于含矿地质体或残积覆盖物中较高的钙和其他伴生元素的含量,通过现场直接测定Ca等元素,可以快速追索含矿地质体。在萤石找矿勘查中,土壤吸附态F离子含量对萤石矿化带具有良好的指示意义,在含矿带上有明显的F元素正异常,可用于隐伏萤石矿体的定位预测。由于X射线荧光分析技术只能探测与萤石密切相关的Ca等元素的含量,不能用于F元素的分析,土壤偏提取地球化学测量可作为补充手段。俄力木台10线Ca元素含量较高,Ca元素高值异常主要分布于矿化露头带并向山坡下方迁移,20线Ca元素含量整体偏高且异常不明显。俄力木台地区出露的第四系沉积物主要为原生岩石的风化堆积物,排除外来沉积物的影响。因此,俄力木台整体偏高的Ca元素含量在一定程度上受林西组钙质沉积岩的影响,在矿化露头带及附近呈现较弱的Ca元素高值异常,可指示萤石矿体的存在。本次研究的土壤偏提取地球化学测量结果在测线东段与甚低频低阻异常带对应部位均显示较明显的F元素高值异常。另外,在10线东段甚低频低阻异常带、便携式X射线荧光分析土壤Ca元素高值异常带和土壤偏提取地球化学测量F元素高值异常带高度吻合(图9),该异常组合可作为本区含矿断裂和隐伏萤石矿的重要预测标志。
图9 俄力木台预测区物化探综合异常图
4.3 综合成矿预测与工程验证
俄力木台地区第四系覆盖严重,地表零星发育硅化强烈的硅质顶盖。结合地表地质填图、地质-甚低频电磁测量、便携式X射线荧光分析、土壤偏提取地球化学测量等综合剖面和甚低频电磁扫面,研究结果认为俄力木台深部可能发育隐伏-半隐伏萤石矿体。地质和物化探综合剖面表明,便携式X荧光分析Ca元素高值异常在10线对应地段较明显,主要分布于矿化露头带及其山坡下方;土壤偏提取地球化学测量,在测线东段与甚低频低阻异常带对应部位均显示有较明显的F元素高值异常。甚低频电磁测量的低阻异常平面分布图指示本区构造基本格局主要由NNE-SN向和NW-NNW向两组断裂组成,前者属主要控矿断裂,后者属该区成矿后断裂并对前者的连续性具有一定的破坏作用。在剖面沿线与已知矿化露头相对应的部位,甚低频电磁测量的低阻异常最为明显,可以作为本区隐伏构造矿化带定位预测的重要依据。甚低频电磁测量对测区北部的已知萤石矿带具有明显的反映,在测区东部发育有一条SN向贯通整个测区的较强的低阻异常带,其东侧还有一条次级异常带。该SN向低阻异常带、F和Ca剖面异常与断续出露的矿化蚀变带相吻合,预示了本区构造矿化异常带的成矿远景,因此在俄力木台地区可能发育两条萤石矿带(图9)。
基于上述地-物-化-遥综合找矿预测结果,本次研究在重点预测地带开展了深部探矿工程,共布置钻孔14个,总计孔深3143m,其中5个钻孔见到工业矿体(CaF2加权平均品位36.89%),6个钻孔见到萤石矿化体。工程验证结果表明,俄力木台地区不同标高均可见萤石矿化(图10)。萤石主要以脉状产出,呈无色、紫色、浅绿色,围岩及夹石中可见黄铁矿化。较高品位的萤石矿体尚未揭露,其萤石主矿体可能埋藏较深,目前揭露的部位仍为硅质顶盖、头部矿体或者构造狭缩带,但其深部具有较好的资源潜力,找矿前景良好。
图10 俄力木台重点预测地带钻孔揭露的萤石矿化特征
5 结论
(1)内蒙古赤峰地区萤石矿体的垂向分带特征与浙江地区较为一致,受断裂控制的萤石矿体在垂向上自上而下包括硅质顶盖、头部矿体、中部矿体和尾部矿体四个部分。
(2)WorldView-2遥感影像揭示俄力木台地区含矿断裂带的空间展布与线性变化特征,矿化带中-北段呈近SN向,南部主要呈NNE-NE向。
地质和物化探综合剖面中甚低频低阻异常带、便携式X射线荧光分析土壤Ca元素高值异常带和土壤偏提取地球化学测量F元素高值异常带高度吻合,该异常组合可作为本区含矿断裂和隐伏萤石矿的重要预测标志。
地-物-化-遥综合研究结果指示俄力木台可能发育隐伏-半隐伏萤石矿体。重点预测地带的钻孔验证见矿效果良好,目前揭露的萤石矿化和蚀变主要为硅质顶盖、头部矿体或者构造狭缩带,俄力木台深部具有较好的资源潜力和找矿前景。
感谢赵鹏大院士在矿产勘查领域的重要贡献,谨以此文祝贺赵鹏大院士九十华诞。林西县自然资源局潘俊峰、郭建龙以及中国地质大学(北京)林西萤石项目组成员在研究过程中提供了帮助,在此一并表示衷心感谢。
来源:唐利,张寿庭,王亮等/地学前缘(EarthScienceFrontiers)2011,28(3)
联系客服