湘中秋旺冲金矿区脉岩岩石学与地球化学特征及其地质意义
2020-08-10彭建堂鲁玉龙徐军伟
周 溪,彭建堂,鲁玉龙,徐军伟,谢 帅,谢 青
(1.湖南省有色地质勘查局 二总队,湖南 湘潭 411102;2.中南大学 地球科学与信息物理学院 有色金属成矿预测与地质环境监测教育部重点实验室,湖南 长沙 410083;3.中国科学院 地球化学研究所 矿床地球化学国家重点实验室,贵州 贵阳 550002)
湘中地区以盛产金、锑矿闻名,且在多个矿区发现有脉岩与金、锑矿化相伴生的现象(刘继顺,1996;鲍振襄等,2002;赵军红等,2005;鲁玉龙等,2017a),如龙山金锑矿、符竹溪金矿、板溪锑矿、廖家坪金锑钨矿、包金山金矿等。前人对湘中地区部分矿区中脉岩与金、锑成矿的关系进行过初步探讨,并认为这些脉岩中长英质脉岩与金、锑矿化存在紧密的时空和成因联系;长英质脉岩主要形成于印支晚期,其侵位时间大多与成矿时间相近,它们为本区金、锑成矿有的提供重要热源和成矿流体,有的直接作为成矿母岩,其侵入作用本身就是成矿作用(刘继顺,1996;鲍振襄等,2002;鲁玉龙等,2017a)。但总体上对该区脉岩研究不多,缺乏系统性的研究,这些脉岩的岩石成因、形成的构造环境等方面的资料仍相当少(赵军红等,2005;鲁玉龙等,2017a)。
秋旺冲金矿区位于湘中EW向白马山-紫云山成矿带东端、紫云山花岗岩体的北东侧。前人在该矿区及其周边地区,对紫云山岩体及该矿区西侧的包金山金矿中的脉岩进行了少量年代学和地球化学研究(刘凯等,2014;鲁玉龙等,2017a,2017b),结果显示,紫云山花岗岩体与包金山金矿中的花岗闪长斑岩脉侵位时间都在印支晚期,且在包金山金矿区,花岗闪长斑岩脉与金成矿时间相近,空间上相互穿插,为成矿期脉岩(鲁玉龙等,2017a)。2015~2016年,笔者项目组在秋旺冲金矿区开展地质填图、钻探等地质勘查工作时,新发现了较多的长英质脉岩(周溪等,2017)(1)周 溪,黄建伟,曾晓霞,等.2017.湖南省湘乡市秋旺冲矿区金矿预查报告.湖南省有色地质勘查局二总队.,且这些脉岩可能与金成矿有关,但因发现较晚,且矿区地质工作程度偏低,目前尚未引起学者的关注,其岩石类型和地球化学特征如何,是否与金成矿关系密切,均不清楚。因此,本文拟对该矿区脉岩进行岩石学和全岩地球化学研究,并结合区域上前人的研究成果,初步探讨这些脉岩的成因、形成的构造环境及其与金矿化在空间上的对应关系,以期为进一步深入认识整个湘中地区脉岩的成因及其与金、锑成矿的关系提供资料。
1 脉岩的分布与岩石学特征
湘中地区大地构造上处于扬子陆块与华夏陆块汇聚带的北侧(图1),区域构造较复杂,整体表现为穹盆构造,岩浆活动也较频繁,花岗岩分布非常广泛(鲁玉龙等,2017b)。受构造-岩浆活动的影响,该地区形成了丰富的金、锑、钨、铅锌等矿产,是湖南重要的贵金属、有色金属成矿区(刘继顺,1996;赵军红等,2005;鲁玉龙等,2017b)。秋旺冲金矿区位于湘中盆地东缘,矿区内出露的地层主要为新元古界板溪群马底驿组(Pt3bnm),为一套浅变质岩,主要为砂质板岩、含钙质粉砂质板岩和变质砂岩。区内断裂较发育,以EW-NWW向为主,NE向和SN向次之,其中EW-NWW向断裂构造是本区金矿的容矿构造。矿区西南部为湘中地区规模较大的紫云山岩体(图1),出露面积达280 km2,岩性主要为似斑状花岗岩、二长花岗岩、黑云母花岗岩(刘凯等,2014;鲁玉龙等,2017b)。矿区脉岩主要出露在南部靠近紫云山花岗岩体界线约500 m范围内,目前发现有5条脉岩出露地表,另据矿区钻孔资料反映深部还有多条隐伏脉岩,规模都不大,以花岗岩脉(Mγ)为主。
该矿区脉岩较发育,均为浅色的长英质脉岩,出露于金矿化带附近的脉岩,大体可分为似斑状花岗岩、含电气石花岗岩和文象花岗岩3类(图2)。这些脉岩的岩石学特征分述如下:
似斑状花岗岩,岩石呈灰白色,似斑状结构,块状构造(图2a),斑晶含量一般为10%~20%。主要矿物为石英(35%~40%)、钾长石(30%~40%)和斜长石(An=2~8,钠长石,20%~25%),副矿物有磷灰石、磁铁矿等,斑晶主要由斜长石和石英组成。石英为他形粒状,与长石互嵌、互为伸入,并交代长石。钾长石以半自形-他形为主,与石英互嵌、互为伸入。斜长石呈半自形-自形板条状,弱绢云母化,边缘常被熔蚀,发育聚片双晶(图2b、2c)。
图1 湘中秋旺冲金矿区地质简图Fig.1 Geological sketch map of the Qiuwangchong gold deposit,central Hunan
含电气石花岗岩,岩石呈灰白色,中细粒花岗结构,块状构造(图2d)。主要矿物为石英(30%~35%)、钾长石(35%~40%)和斜长石(An=5~11,钠长石-更长石,15%~20%),次要矿物为电气石(5%~10%)、白云母(1%~3%)和少量石榴子石,副矿物有磷灰石、锆石、磁铁矿等。石英为他形粒状,与钾长石互嵌、互为伸入,并交代斜长石。钾长石以他形为主,少量具格子双晶,可见卡氏双晶,表面泥化显混浊,与石英互嵌、互为伸入(图2e、2f)。斜长石呈半自形-自形板条状,弱绢云母化,边缘常被熔蚀,发育聚片双晶(图2e)。电气石多呈短柱-长柱状(图2e),偶见楔状,横切面呈不规则状(图2f),横向裂理发育,具强的吸收性、多色性,见彩色环带,在正交偏光下其干涉色常被掩盖,总体呈较均匀分散分布在长石、石英粒间。白云母呈片状,具鲜艳的干涉色,零散分布,部分碳酸盐化。石榴子石颗粒极小,呈圆形颗粒分布于石英中,单偏光镜下无色,突起较高,正交偏光下全消光。
文象花岗岩,岩石呈灰白-肉红色,石英在钾长石中呈有规则象形文字产出,构成显微文象结构(图2h、2i),块状构造。主要矿物为钾长石(60%~65%)、石英(25%~35%),次要矿物有斜长石(An=4~15,钠长石-更长石,3%~5%),少量黑云母和白云母,副矿物有磷灰石、锆石、磁铁矿等。钾长石颗粒粗大,具条纹结构,在其中石英呈象形文字状嵌布,表面泥化混浊状。石英为它形,呈象形文字状嵌布在钾长石中,在正交偏光镜下同时消光。斜长石呈半自形板条状,见聚片双晶发育,边缘常被熔蚀呈不规则状,见绢云母化。黑云母和白云母分布不均匀,局部位置可见少量;前者呈自形、半自形板状、片状,后者呈片状,具鲜艳的干涉色。
2 样品及测试
从秋旺冲金矿区V14金矿化带附近发现的3类脉岩中采集了3组共10个脉岩样品(采样位置见图1),来进行岩石地球化学分析,样品均采自钻孔岩芯,新鲜,无风化现象。
主量元素分析采用X射线荧光光谱法(XRF),在广州澳实矿物实验室完成。总氧化物分析误差为1%~3%。实验流程为:称取适量样品→加入适量硼酸锂-硝酸锂高温熔融成玻璃片→在XRF上用外标法测定氧化物含量。
图2 秋旺冲金矿区脉岩手标本和显微照片(矿物缩写符号据沈其韩,2009)Fig.2 Hand specimens and microscopic photographs of felsic dykes in the Qiuwangchong gold deposit (mineral abbreviation after Shen Qihan,2009)a—似斑状花岗岩手标本;b—斜长石聚片双晶(+);c—斜长石聚片双晶、它形粒状石英(+);d—含电气石花岗岩手标本;e—斜长石聚片双晶、短柱状电气石(+);f—电气石彩色环带(+);g—文象花岗岩手标本;h—石英、钾长石文象结构(+);i—石英、钾长石文象结构 (+);Pl—斜长石;Kfs—钾长石;Qtz—石英;Tur—电气石a—hand specimens of porphyritic granite;b—polysynthetic twin of plagioclase (+);c—polysynthetic twin of plagioclase and xenomorphic granular quartz (+);d—hand specimens of tourmaline-bearing granite;e—polysynthetic twin of plagioclase and prismatic tourmaline (+);f—colored band of tourmaline (+);g—hand specimens of hebraic granite;h—graphic texture formed by quartz and potassium feldspar (+);i—graphic texture formed by quartz and potassium feldspar (+);Pl—plagioclase;Kfs—potassium feldspar;Qtz—quartz;Tur—tourmaline
微量元素分析采用电感耦合等离子质谱法(ICP-MS),在贵州同微测试科技有限公司完成。测试所使用仪器型号为德国热电(Thermo Fisher)ICP-MS X2,稀释后的上机溶液含4个内标,原始数据用昆士兰大学俸月星博士编写的程序计算含量,每个样品每个元素均测试4位,RSD值为0.12%~8.66%。
3 地球化学特征
3.1 主量元素特征
主量元素测试分析结果(表1)显示,3种长英质脉岩的SiO2含量介于72.45%~84.52%之间(平均为76.88%),属酸性-超酸性岩;在全碱-硅(TAS)分类图解(图3a)中所有样品都落到花岗岩区域。脉岩的Al2O3含量介于8.37%~14.92%之间,均值为12.58%;A/CNK值介于0.94~1.18之间,除一个样品略大于1.1外,其它均小于1.1;A/NK值为1.05~1.26,在A/NK-A/CNK图解中(图3b),所有样品均落入准铝质-弱过铝质区域,因此,秋旺冲金矿区的脉岩应为准铝质-弱过铝质系列。脉岩的全碱(K2O+Na2O)介于5.18%~11.11%之间,均值为8.25%;K2O/Na2O值变化较大,为0.32~3.41。样品里特曼组合指数 值为0.65~4.14,大多小于3.3,均值为2.18,总体属于钙碱性系列;该区脉岩在MALI-SiO2(碱钙指数)图解(图4)中也大多落入钙碱性区域,与里特曼组合指数所反映的特征一致。该区脉岩的分异指数(DI)普遍高,为93.66~98.12,显示岩浆分异演化程度高。
表1 秋旺冲金矿区脉岩主量元素分析结果 wB/%Table 1 Major element compositions of felsic dykes in the Qiuwangchong gold deposit
图3 秋旺冲金矿区脉岩TAS图解(a,据Middlemost,1994)和A/NK-A/CNK图解(b,据Maniar and Piccoli,1989)Fig.3 TAS diagram (a,after Middlemost,1994) and A/NK-A/CNK diagram (b,after Maniar and Piccoli,1989) of felsic dykes in the Qiuwangchong gold deposit
图4 秋旺冲金矿区脉岩MALI-SiO2图解(据Frost et al.,2001)Fig.4 MALI-SiO2 diagram of felsic dykes in the Qiuwangchong gold deposit (after Frost et al.,2001)
在主量元素Harker图解中(图5),该区脉岩SiO2与Al2O3呈明显负相关性,与K2O也存在一定程度的负相关,与其他主量元素相关性不明显,表明在岩浆结晶过程中有富铝、富钾矿物的结晶分离(刘凯等,2014),且在岩浆演化过程中可能存在斜长石的分馏结晶作用(赵军红等,2005)。脉岩样品中TiO2、MnO、P2O5含量都很低,部分样品含量甚至低于0.01%,说明可能在脉岩结晶过程的早中期大部分钛矿物、磷灰石等矿物就已发生结晶分离。脉岩中SiO2含量高而CaO、MgO、Fe2O3含量较低,整体表现出酸度高、相对贫钙、贫镁、贫铁的特点。3组样品主量元素之间这些相同的特点,也反映了该矿区这3种类型的脉岩应为同一母岩分异演化形成的。
3.2 微量元素特征
微量元素测试分析结果见表2。在原始地幔标准化蛛网图中(图6a),10件样品微量元素配分型式基本一致,有2个样品略有差异(富钾,与主量元素特征一致),均属于右倾型。3种脉岩样品都具有较高的Rb含量及Rb/Sr值,强不相容元素Rb高度富集,表明岩浆分异作用进行得比较充分(王永磊等,2007;尹光候等,2010),花岗岩在形成过程中经历了较高程度的演化。岩石中强烈亏损Ba、Sr、P、Ti元素,呈明显的V型谷,其中P、Ti元素亏损说明在岩浆分异过程中可能存在磷灰石和钛铁矿的分离结晶作用,Ba、Sr元素亏损可能是与斜长石的分离结晶作用有关(赵军红等,2005),Sr的负异常可能还与岩浆分异或岩浆中钙含量偏低有关。岩石中富集Rb、K等大离子亲石元素(LILE),Th、U、Ce、Nd、Zr、Hf这些高场强元素(HFSE)也都表现出一定程度的富集。该区脉岩的岩石微量元素正负异常(峰、谷)都比较突出,反映了其较强的分异演化特点(张利军,2015)。
该区脉岩的K/Rb、Zr/Hf、Nb/Ta值均较小,暗示其岩浆分异演化过程中有流体的参与(Shaw,1968),并发生了壳、幔之间流体和熔体的相互作用,同时说明本区长英质脉岩的岩浆分异演化程度相对较高(Irber,1999;刘朕语,2018)。其微量元素特征与紫云山岩体相比,相对富集或亏损的元素基本相同,曲线分布型式大体相似,但K/Rb、Nb/Ta值都明显降低(刘凯等,2014;鲁玉龙等,2017b),具有更高的演化程度。
3.3 稀土元素特征
稀土元素测试分析结果(表2)显示,该区脉岩的稀土元素总量相当低,介于10.4×10-6~94.8×10-6之间,均值为51.1×10-6。在稀土元素球粒陨石标准化配分型式图(图6b)中,所有样品的稀土元素配分型式基本相同,类似海鸥型,说明秋旺冲矿区3种脉岩的岩浆应是由同一母岩浆演化而来(张舒等,2009);轻稀土元素较为平坦,重稀土元素趋势总体平坦略微左倾,且轻、重稀土元素分馏都不明显。轻稀土元素总量为4.2×10-6~29.9×10-6,均值11.5×10-6,重稀土元素总量为6.2×10-6~83.0×10-6,均值39.6×10-6;轻、重稀土元素比值为0.14~0.85(平均为0.35),(La/Yb)N值介于0.17~1.04间,大多小于1.0,表明重稀土元素相对较富集,而轻稀土元素相对亏损。
脉岩δEu值为0.01~0.15,均值为0.07,强烈亏损Eu,稀土元素配分型式图上呈明显V型谷,说明岩浆演化过程中有大量的斜长石分离结晶,暗示该区脉岩的岩浆经历了强烈结晶分异作用,分异演化程度很高,这与其主量、微量元素所反映的地球化学特征相吻合;Ce异常不明显,δ Ce值大都在1.0附近。与邻近的紫云山岩体相比,该区脉岩的LREE及LREE/HREE明显降低(刘凯等,2014;张利军,2015;鲁玉龙等,2017b),应是岩浆分异作用进行得更加充分。
图5 秋旺冲金矿区脉岩Harker图解Fig.5 Harker diagrams of felsic dykes in the Qiuwangchong gold deposit
4 讨论
4.1 脉岩成因类型的探讨
花岗岩分类方案有多种,其中I、S、M、A分类方案较为常用。前述已对秋旺冲金矿区脉岩做了岩石学与地球化学特征分析,其主量元素具有富Si、Na、K而贫Al、Ca、Mg的特点;微量元素具有富集Rb、Th、U、Ce、Nd、Zr、Hf而亏损Ba、Sr、P、Ti的特点;全碱(K2O+Na2O)含量较高,里特曼指数大多小于3.3,总体属于钙碱性花岗岩;矿区脉岩的A/CNK值大多小于1.1,绝大部分样品落入准铝质-弱过铝质区域(图3b),总体上具有I型花岗岩的特征(Chappell,1999),与紫云山岩体基本一致(刘凯等,2014;鲁玉龙等,2017b)。在Y-Ga/Al及Zr-Ga/Al判别图解(图7a、7b)中,该区脉岩的样品主要落入I &S型花岗区域,只有个别样品落入A型花岗岩范围;在Ce-SiO2判别图解(图7c)中,所有样品点均落入I型花岗岩范围;在(K2O+Na2O)/CaO- (Zr+Nb+Ce+Y)判别图解(图7d)中,样品点则主要落入高分异I型花岗区域,表明本区脉岩可能兼具I型和A型花岗岩的部分特征。但矿区脉岩中并未见有碱性暗色矿物,且10 000 Ga/Al值介于1.44~2.73之间,除一个样品值略大于2.6外,其它均小于2.6,故该区脉岩不应是A型花岗岩,其部分与A型花岗岩相似的特征,很可能是脉岩岩浆高度分异(图8)所致因此,根据秋旺冲金矿区脉岩的地质和地球化学特征以及成因类型判别研究可知,该区花岗岩脉应属于高分异的Ⅰ型花岗岩,应是其母岩紫云山花岗岩体进一步高度分异演化的产物。
表2 秋旺冲金矿区脉岩微量及稀土元素分析结果 wB/10-6Table 2 Trace element and REE compositions of felsic dykes in the Qiuwangchong gold deposit
图6 秋旺冲金矿区脉岩微量元素原始地幔标准化蛛网图(a)和稀土元素球粒陨石标准化配分型式图(b)(原始地幔和球粒陨石数据据Sun and McDonough,1989)Fig.6 Spidergrams of the primitive mantle-normalized trace elements (a) and chondrite-normalized REE patterns (b) of felsic dykes in the Qiuwangchong gold deposit (primitive mantle-normalized values and chondrite-normalized values after Sun and McDon- ough,1989)
图7 秋旺冲金矿区脉岩Y-Ga/Al、Zr-Ga/Al、Ce-SiO2和(K2O+Na2O)/CaO-(Zr+Nb+Ce+Y)成因类型判别图解(a、b 据Whalen et al.,1987;c 据Collins et al.,1982;d 据Eby,1992)Fig.7 Y-Ga/Al,Zr-Ga/Al,Ce-SiO2 and (K2O+Na2O)/CaO-(Zr+Nb+Ce+Y) genetic discrimination diagrams of felsic dykes in the Qiuwangchong gold deposit (a,b after Whalen et al.,1987;c after Collins et al.,1982;d after Eby,1992)
图8 秋旺冲金矿区脉岩Rb-Ba-Sr三元图解(据谢磊等,2013)Fig.8 Rb-Ba-Sr diagrams of the felsic dykes in the Qiuwangchong gold deposit (after Xie Lei et al.,2013)
4.2 脉岩形成的构造环境探讨
秋旺冲金矿区脉岩大都出露在距紫云山花岗岩体500 m的范围内,因而对其形成的构造环境的探讨应与紫云山岩体相结合。以往的研究显示,紫云山岩体属于华南印支晚期(锆石U-Pb SIMS年龄为227~222 Ma)花岗岩,为碰撞晚期或碰撞后构造背景的产物(刘凯等,2014;鲁玉龙等,2017b)。
在微量元素Nb-Y、Rb-(Yb+Ta)构造环境判别图解(图9a、9b)中,所有样品点均落入同碰撞花岗岩(Syn-COLG)或火山弧花岗岩(VAG)区域内;在Rb-(Y+Nb)判别图解(图9c)中,样品点主要落入碰撞后花岗岩(Post-collision)区域内。在TFeO/(TFeO+MgO)-SiO2判别图解(图9d)中,矿区脉岩样品点主要落入造山后花岗岩(POG)区域。上述构造环境判别图解表明,秋旺冲矿区脉岩的总体上属于碰撞后花岗岩,同时也与造山作用有关。
华南陆块位于秦岭-大别山和Song Ma两条印支期缝合带之间,而本区在大地构造上则位于华南陆块内部的扬子陆块与华夏陆块的拼合的边缘薄弱地带(图1)。在印支运动中,印支陆块与华南陆块南缘发生碰撞拼合,所产生的挤压应力由南向北传递(刘凯等,2014),华南陆块内部的薄弱地带随之发生地壳加厚而产生造山作用(刘凯等,2014;张利军,2015),并由此伴随着强烈的构造-岩浆作用,而该区脉岩很可能就形成于此时的造山后伸展环境中。
4.3 脉岩与金矿化的关系
在秋旺冲金矿区,已发现的金矿(化)体主要赋存于EW-NWW向的断裂构造中(图1),在该组断裂构造中,附近有脉岩密集发育并在空间上与之相互穿插的地段也是金矿化较强之处(V14中段)。矿区地质勘查资料显示,本文所研究的3类脉岩其钻孔原生晕数据中Au含量为4.5×10-9~14.7×10-9(石墨炉原子吸收光谱法测定)(周溪等,2017)(2)周 溪,黄建伟,曾晓霞,等.2017.湖南省湘乡市秋旺冲矿区金矿预查报告.湖南省有色地质勘查局二总队.,是中国陆壳Au丰度(2.25×10-9)(黎彤,1994)的2.0~6.5倍;钻孔一般岩矿分析数据中脉岩的Au含量为0.03×10-6~0.20×10-6(泡沫塑料富集火焰法测定)(周溪等,2017)(3)周 溪,黄建伟,曾晓霞,等.2017.湖南省湘乡市秋旺冲矿区金矿预查报告.湖南省有色地质勘查局二总队.,且这两种方法测定的数据都显示,越靠近金矿(化)体的脉岩其金含量越高,在空间上与金矿(化)体相互穿插的脉岩本身就是弱金矿化体。总之,该区脉岩在空间上与金矿化存在对应关系,对研究金矿的成因及下一步的找矿勘查,均具有一定指示意义。
秋旺冲金矿区邻近的包金山金矿区,长英质的脉岩也相当发育,且矿体与脉岩存在相互切穿关系,脉岩与金成矿关系密切(鲁玉龙等,2017a)。最近的定年结果亦显示,包金山金矿形成于印支期,其成矿年龄与长英质脉岩及紫云山岩体的侵位年龄相当吻合(鲁玉龙等,2017a,2017b)。因此,可以推断,包括秋旺冲、包金山等矿区在内的湘中地区的长英质脉岩,在成因上,与该区的金成矿也应存在密切联系。
5 结论
(1) 秋旺冲金矿区脉岩岩性主要有似斑状花岗岩、含电气石花岗岩和文象花岗岩,都为浅色长英质脉岩。
(2) 该区脉岩为准铝-弱过铝质、钙碱性花岗岩;富硅、富碱,富集Rb、Th、U、K 、Ce、Nd、Zr和Hf元素,强亏损Ba、Sr、P和Ti元素;稀土元素总量低,轻重稀土元素分馏均不明显,重稀土元素相对较富集,Eu负异常极为明显;分异指数高,显示岩浆分异演化程度高。
(3) 该区脉岩成因类型属于高分异的Ⅰ型花岗岩,应是其母岩紫云山花岗岩体进一步高度分异结晶的产物。
图9 秋旺冲金矿区脉岩的构造环境判别图解(a、b 据Pearce et al.,1984;c 据Pearce,1996;d 据Maniar and Piccoli,1989)Fig.9 Tectonic discrimination diagram of felsic dykes in the Qiuwangchong gold deposit (a,b after Pearce et al.,1984;c after Pearce,1996;d after Maniar and Piccoli,1989)ORG—洋中脊花岗岩;WPG—板内花岗岩;VAG—火山弧花岗岩;Syn-COLG—同碰撞花岗岩;Post-collision—碰撞后花岗岩;IAG—岛弧花岗岩;CAG—大陆弧花岗岩;CCG—大陆碰撞带花岗岩;POG—造山后花岗岩;RRG—裂谷花岗岩;CEUG—大陆造陆隆起花岗岩ORG—ocean ridge granites;WPG—within plate granites;VAG—volcanic arc granites;Syn-COLG—syn-collision granites;Post-collision—post-collision granites;IAG—island arc granites;CAG—continental arc granites;CCG—continental collision granites;POG—post orogeny granites;RRG—related granites;CEUG—continental epeirogenic uplift granites
(4) 该区脉岩是在碰撞造山后伸展构造环境中形成的,属于碰撞后花岗岩。
(5) 该区脉岩与金矿化存在空间对应关系,在成因上也可能关系密切,这些脉岩对金矿床的成因及找矿勘查均有一定的指示意义。
致谢匿名审稿人和编辑在文章修改过程中提供了许多建设性的意见,在此致以诚挚的谢意!