1. 引言
雪峰弧形构造隆起带位于扬子板块与华夏板块之间的过渡地段，大致沿通道–安化–桃江一线呈北西凸出的弧形展布，其在湖南境内延长约600 km，宽可达数十千米(图1)。雪峰弧形构造隆起带是华南重要的金成矿带，湖南省已有的金矿床(点)约80%分布于该区冷家溪群和板溪群浅变质岩中 [1] 。桃源县沅水以南地区(以下称研究区)位于雪峰弧形构造隆起带的转折部位，在研究区周边先后发现沅陵县沃溪金锑钨矿床、安化县符竹溪金锑矿床、桃江县板溪锑矿床等。上世纪60~70年代研究区先后进行了多次地质勘查工作，圈出金、锑、钨重砂异常数十处，发现了西安钨金矿 [2] [3] 、冷家溪 [4] [5] 、小桃源 [6] 、马家坡 [7] 等小型金矿，另外还有数十个金矿(化)点。
上世纪80年代以来，许多专家和学者对研究区及周边矿床的成矿物质来源、成矿时代、矿区构造等方面的问题，提出了许多不同的见解和观点。成矿物质来源目前主要有三种观点：鲍振襄 [8] 认为赋矿围岩即矿源层；彭建堂 [9] 认为成矿物质主要来源于下伏更成熟的陆壳基底；毛景文 [10] 则认为成矿物质一部分来源于中晚元古宙地层，另一部分来自深部或燕山期花岗岩体。成矿时代认识争议较大，譬如关于沃溪金矿成矿时代，顾雪祥 [11] 认为形成于新元古代雪峰期，彭建堂 [9] 认为形成于加里东期，彭渤 [12] 认为形成于中生代，李彬 [13] 认为形成于燕山期的两阶段幕式成矿作用。关于矿区构造，目前的主要认识是该区处于雪峰弧形构造带北东段，是一个多期叠加的变形带，矿床受韧性剪切带、节理裂隙带、劈理化带和褶皱构造带等控制明显 [14] [15] 。
前人研究侧重于矿点、矿区等大比例尺的工作，中比例尺成矿预测方面的工作较为薄弱。研究区为侵蚀剥蚀低山—丘陵地貌，植被茂盛，基岩出露较差，地表较难发现矿化信息，土壤地球化学测量对于植被覆盖区找矿或寻找隐伏矿床有很好效果 [16] [17] [18] 。为了明确研究区资源潜力，湖南省地质调查所开展了1:50000土壤地球化学测量(重点矿区加密至1:10000)。本次研究基于土壤地球化学测量扫面成果，采用含量—面积分形法确定土壤地球化学异常，结合地质背景对异常进行评价，圈定找矿靶区，以推进该区地质找矿工作取得新突破。
2. 区域地质背景
2.1. 地层与岩浆岩
Figure 1. Geological and mineral distribution diagram of the Xuefeng arc-shaped structural uplift zone
图1. 雪峰弧形构造隆起带地质及矿点分布简图
研究区出露的地层有新元古界冷家溪群、板溪群、南华系，下古生界寒武系，上古生界中泥盆统、中上石炭统，中生界白垩系和新生界第三系、第四系(图2)。冷家溪群主要出露于西南部木旺溪、牯牛山一带，东部向阳村、西北部竹园村一带零星出露，为浅海相变质砂泥质碎屑岩建造。板溪群在研究区南部广泛分布，为浅海相变质砂泥质、凝灰质碎屑岩建造。南华系在铁山溪村、沙堤村一带呈带状分布，为含砾砂岩、冰碛砾岩建造。震旦系、寒武系毗邻南华系分布，为硅质岩、碳质页岩建造。上古生界泥盆系、石炭系在研究区中部羯羊铺村零星出露，范围很小。中生界白垩系发育完整，在研究区中部、北部广泛分布，为河湖相沉积。新生界分布在沅水两侧，为河流沉积。
岩浆岩出露于茶庵铺镇北西方向的走马岗和梅子溪两处，为辉绿岩脉群。走马岗有辉绿岩小岩体12个，呈雁行状排列；梅子溪出露辉绿岩体2个。辉绿岩脉侵入于新元古界地层中，一般长1000~1500 m，宽30~50 m，呈岩床和岩墙产出，常见分支复合现象。
2.2. 构造
研究区处于雪峰隆起的东北端及中–新生代洞庭凹陷的西南角，影响较大的构造运动有武陵运动、雪峰运动、加里东运动和燕山运动。按构造线展布方向可分为东西向构造、北东向构造、北北东向构造等，研究区以东西向构造、北东向构造为主，构造线总体方向为60˚~90˚。东西向构造由一系列走向大致平行的北东–东西向压扭性断裂和褶皱组成，该构造带控制了区内金矿床的形成，是省内重要的金矿成矿远景区。北东向构造由一系列大致互相平行的压性、压扭性断裂和韧性剪切带组成。
2.3. 矿产
Figure 2. Geological and mineral resources map
图2. 研究区地质矿产图
研究区矿产资源丰富，主要有石煤、铁、钒、铜、铅、锌、白钨、锑、金、黄铁矿、重晶石、金刚石等十多种。金矿是区内的优势矿产，岩金以冷家溪、小桃源、西安、沙坪、寺坪等金矿点组成东西向成矿带。砂金一般沿沅江河床、第四系砂砾层分布。另外，沿东西向西安溪压性断裂带有一条白钨矿带分布，西起磨子坪，向东经西安溪直到桃安，长约15 km。
3. 样品采集与测试方法
调查区总面积1477.33 km2，共采集测试了2193个面积性土壤样品(图3)。其中，农用地面积202.55 km2，实际采样1148个，采样密度5.67个/km2。林地面积1152.47 km2，实际采样1031个，采样密度0.89个/km2。其他用地(采矿用地、空闲地、裸土地等) 3.26 km2，实际采样14个，采样密度4.29个/km2。
样品采集按照多目标区域地球化学调查规范执行，在布设采样点周围20米范围内采用梅花形布点法(平缓区)或蛇形布点法(低山丘陵区)等量采集5个子样组成一个样品。农用地采集地表至20 cm深处的土柱，林地、园地采集深度0~60 cm，采集时去除杂草、草根、砾石、砖块、肥料团块等杂物。
Figure 3. Land use status and sampling point map
图3. 研究区土地利用现状及采样点位图
样品加工按照《土地质量地球化学评价规范(DZ/T 0295-2016)》执行，样品检测由湖南省地质调查所测试中心(具备中国地质调查局多目标地球化学调查样品测试能力(54项)资质)完成，检测过程依据《区域地球化学样品分析方法(DZ/T 0279-2016)》执行。该项目测定了Au、Ag、Sb、W、As、Pb、Zn、Cu、Mn、Mo、Hg、Cr、Ni等元素。其中Cu、Mn、Mo、Pb、Zn、Ni、Cr采用电感耦合等离子体发射光谱法测定，Au、Sb、W采用电感耦合等离子体质谱法测定，Ag采用垂直电极—发射光谱法测定，As采用氢化物发生—原子荧光光谱法测定，Hg采用蒸气发生—冷原子荧光光谱法测定。
4. 统计结果
4.1. 成矿元素的统计特征
运用SPSS软件的描述性统计功能对数据进行统计，统计参数包括最大值、最小值、极差、平均值、中位数、标准差、方差、变异系数、偏度、峰度等。统计结果见表1。
Table 1. Statistical results of soil metallogenic elements (n = 2193)
表1. 土壤成矿元素统计结果(n = 2193)
注：Au单位ng/g，其它元素为μg/g，下同。
变异系数(Cv)主要描述元素分异特征，变异系数大，表明元素在某一地区或某种地质体中的分布分配不均匀，离散程度大，易于活动迁移形成矿(化)体或强异常；变异系数小，表明元素分布分配较均匀，离散程度小，多形成弱异常或低背景，一般仅具有某种地质意义，元素本身找矿意义不大。本次研究将变异程度分为4种类型：CV < 30%为均匀分布，30% ≤ CV < 60%为弱分异，60% ≤ CV < 100%为较强分异，CV ≥ 100%为强分异。针对原始数据的变异系数，参照这一分类标准显示，Au、W、Hg、Mo、As、Ag、Sb、Pb为强分异，Mn为较强分异，Cu、Zn、Ni、Cr为弱分异。
原始数据的变异系数为CV1，采用3倍标准差法(平均值 ± 3倍标准差)处理异常值后的变异系数为CV2，CV1/CV2的比值可以反映数据处理时对特高值特低值的削平程度 [19] 。运用CV1和CV1/CV2制作元素变异系数散点图，可以对元素的成矿性进行评价 [20] 。由图4可以看出，Au、W高值多，分异强，成矿的可能性较大；Hg、Mo、As、Ag、Sb、Pb高值较多，分异较强，具有一定的成矿潜力；Mn、Cu、Zn、Ni、Cr高值少，分异弱，成矿的可能性小。
Figure 4. Scatter plot of element variation coefficient
图4. 元素变异系数散点图
4.2. 成矿元素组合特征
4.2.1. R型聚类分析
对于由n个变量组成的总体，变量之间有一定的内在联系，因此可以用相似性系数作为R型聚类分析的一种统计量来定量描述各个变量之间的相似程度，相似性系数越接近1，变量之间的相似程度越高 [21] [22] 。通过R型聚类分析可以较好地了解元素间的共生组合关系。
Figure 5. R-type cluster analysis genealogy
图5. R型聚类分析谱系图
从图5可以看出，在r = 0.3的相似水平上，13个元素可明显聚为3簇：第1簇为Au-Sb-W-Hg，主体由低温元素组成，表明Au、Sb、Hg元素在迁移、富集过程中存在协同关系。在r = 0.5的相似水平上，该簇又细分为Au、Sb-W、Hg三小簇，表明Sb-W的化学行为相关性更强，同时伴有高温元素W，暗示元素富集与岩浆活动有关，研究区Au成矿温度较高。第2簇为As-Pb-Zn-Cu-Cr-Ni-Mn，在r = 0.7的相似水平上，As、Pb、Mn、Zn-Cu-Cr-Ni形成四簇；在r = 0.8的相似水平上，Zn-Cu-Cr-Ni又细分为Zn-Cu、Cr-Ni两小簇，Zn、Cu为中低温亲硫元素，Cr、Ni为高温元素，表明研究区可能存在两次或两次以上的构造—岩浆活动。第3簇为Ag-Mo，在r = 0.2的相似水平上，Ag-Mo与其他元素分为两簇，表明Ag-Mo在迁移、富集过程中具有很强的独立性。
4.2.2. 因子分析
因子分析是作为一种从众多变量数据中提取若干共性因子的多元统计，能很好地揭示元素的共生组合特点，为矿床成因和异常解释提供依据 [23] [24] 。
KMO (Kaiser-Meyer-Olkin)值为检验变量之间相关关系的指标，KMO指标 > 0.6则表明变量间相关性强，适合因子分析；Bartlett’s球度检验值表示各变量间的独立性，其值 < 0.05时适合做因子分析。检验结果显示，KMO值为0.826，大于0.6；Bartlett球形检验给出的相伴概率Sig.为0.000***，小于显著性水平0.05，表明各元素之间具有相关关系，适合进行因子分析且具有较高的准确性。
以旋转载荷平方和大于1、累计方差贡献率68.529%、因子载荷绝对值大于0.66为阈值，可以取出3个主因子(表2)，这3个因子(F1、F2、F3)基本反映了区内13种元素的地球化学信息。
Table 2. Component matrix after rotation
表2. 旋转后的成分矩阵
因子F1的方差贡献率为36.541%，Cu、Ni、Cr、Zn、Mn、Pb、As有较高的正载荷。Cu、Zn、Mn、Pb、As高值区与冷家溪群、板溪群、南华系等老地层高度重合，低值区与白垩系、第三系、第四系等中新生代地层高度重合。Ni、Cr高值区近东西向，与区内断裂走向重合。说明F1为与地层、构造、区域变质等地质背景相关的因子。
因子F2的方差贡献率为17.845%，Au、Hg、Sb、W有较高的正载荷。Au的高值中心分布于金矿产地或金矿化点附近，W的高值中心分布于西安钨金矿附近，Hg、Sb高值区分布于冷家溪金矿附近。说明F2为与金矿、钨金矿等矿脉有关的因子。
因子F3的方差贡献率为14.144%，Ag、Mo有较高的正载荷。Ag、Mo除了在震旦纪—寒武纪的黑色岩系中含量较高，在农业生产区含量值也较高，可能与上世纪七十年代开展的磷矿粉肥效试验有关，因此F3为与农业生产活动相关的因子。
以上分析结果表明，研究区Au、Hg、Sb、W的成矿潜力最大。
5. 异常识别与找矿靶区确定
5.1. 异常下限的确定
传统的地球化学异常下限的确定主要是使用统计学方法，以样品数据呈正态分布为假设前提，通过计算数据的统计学参数对异常下限进行确定，一般是以平均值与K倍的标准离差之和作为地球化学异常下限值。成秋明 [25] [26] [27] 等认识到这种方法存在局限性，从分形的观点(地球化学背景和异常的形成是两个独立的过程)认为地球化学背景值和异常值具有各自独立的幂律关系，由此导致了一种多重分形分布。该方法较传统方法具有一定优势，异常较强时，多重分形方法所确定的异常下限较传统方法确定的异常下限偏高，可以弱化高背景异常，缩小异常查证范围；异常较弱时，多重分形方法所确定的异常下限较传统方法确定的异常下限偏低，可以强化低背景异常，不会遗漏弱异常。
本次研究采用含量—面积模型处理Au、Hg、Sb、W数据，处理过程如下：首先将原始数据网格化，绘制元素含量等值线图。然后以r为等值线含量值，N(r)为r等值线所围成的平面面积，这样就得到N(r1)，N(r2)，…，N(rn)和(r1, r2, …, rn)的数据对，将这些数据代入分形模型公式中，应用最小二乘法分段拟合成多段直线。最后选择各个区间的拟合直线和原数据间的残差平方和最小的值作为分界值。最终拟合结果见表3。
Table 3. Fitting results of content-area model
表3. 含量–面积模型拟合结果
应用上述方法，将分界点值进行反对数换算，还原为真数，即为分形方法确定的异常下限值。研究区Au、Hg、Sb、W的异常下限计算值分别为44.03 ng/g、0.61 μg/g、8.53 μg/g、22.93 μg/g。结合研究区实际情况对计算值取整数作为元素异常下限，Au、Hg、Sb、W异常下限最终取值分别为40 ng/g、1 μg/g、10 μg/g、20 μg/g。
5.2. 单元素异常特征
本次研究按元素异常下限的1、2、4倍圈定单元素浓度分带，外带为1倍异常下限值，中带为2倍异常下限值，内带为4倍异常下限值。Au、Hg、Sb、W异常特征分述如下：
研究区Au异常强烈，最高值达2418.00 ng/g，是异常下限的60倍。Au总异常面积46.34 km2，浓度分带明显(图6(a))。异常区外带含量介于40~80 ng/g，面积23.83 km2，中带含量80~160 ng/g，面积14.21 km2，异常内带含量大于160 ng/g，面积8.30 km2。面积最大的Au异常区分布在杨溪桥镇中南部煌山村、冷家溪村、竹山村，近东西向展布，其他异常分布规律不明显。该异常区地层主体为新元古界冷家溪群小木坪组，异常展布与断裂密切相关，两条北东东向断裂贯穿整个异常区，异常浓集中心正处在两条断裂之间，异常浓集中心与冷家溪金矿点高度吻合。
研究区Hg异常强度中等，最高值达21.86 μg/g，是异常下限的20倍。Hg总异常面积8.16 km2，浓度分带明显(图6(b))。异常区外带含量介于1~2 μg/g，面积5.50 km2，中带含量2~4 μg/g，面积1.39 km2，异常内带含量大于4 μg/g，面积1.27 km2。Hg异常集中分布在杨溪桥镇朝阳庵村和落马洞村交界处，Hg异常形态与Au异常形态相近，但两者浓集中心有一定偏离。
研究区Sb异常中等，最高值达130.00 μg/g，是异常下限的13倍。Sb总异常面积61.00 km2，浓度分带明显(图6(c))。异常区外带含量介于10~20 μg/g，面积44.69 km2，中带含量20~40 μg/g，面积8.75 km2，异常内带含量大于40 μg/g，面积7.56 km2。Sb异常分布形态和Hg高度重合，浓集中心也十分相近。Hg、Sb是典型的低温成矿元素，迁移能力强，一般位于原生晕的前缘晕，常作为寻找金矿的远程指示元素，同时也是金矿的伴生元素。
研究区W异常强烈，最高值达761.47 μg/g，是异常下限的38倍。W总异常面积45.82 km2，浓度分带明显(图6(d))。异常区外带含量介于20~40 μg/g，面积24.41 km2，中带含量40~80 μg/g，面积14.86 km2，异常内带含量大于80 μg/g，面积6.55 km2。面积最大的W异常区分布在西安镇的西安村和东安村，另外在茶庵铺镇的木旺溪村异常面积也较大，其他异常均较小。W异常区地层主体为新元古界冷家溪群黄浒洞组，异常展布与断裂密切相关，异常浓集中心与西安钨金矿点、木旺溪金矿化点高度重合。
(a) 
(b) 
(c) 
(d)
Figure 6. Single element anomaly map of Au, Hg, Sb, and W
图6. Au、Hg、Sb、W单元素异常图
5.3. 典型矿床地质特征及找矿靶区
5.3.1. 西安钨金矿床地质特征
矿体(脉)赋存于新元古界板溪群马底驿组下段白云质灰岩(大理岩)中，矿脉分布在矿区北部、东南部和西南部，呈带状展布。矿区有上、下两层矿体。上矿层位于马底驿组下段灰紫色板岩与灰绿色板岩交界处的白云质灰岩中，含矿层厚0~11 m；下矿层位于板溪群马底驿组与下伏冷家溪群不整合面上，或在不整合面上相距数米到十余米的白云质灰岩中，含矿层厚0.5~10 m。
矿石类型为硫化物–石英–方解石网脉带型白钨矿，细分网脉带型白钨矿、石英脉型白钨矿两种类型，以网脉带型白钨矿为主。
围岩蚀变与钨矿化关系密切，围岩蚀变以硅化、褪色化为主，其次为绿泥石化、黄铁矿化、绢云母化。在硅化强烈的含矿层上其盖层(围岩)褪色化宽度大的部位，常可找到钨矿体。
结合上述矿床地质特征和图7可以看出，西安型钨金矿矿化富集规律如下：① 矿床赋存于新元古界板溪群马底驿组底部与下伏冷家溪群不整合面上的白云质灰岩及马底驿组下段的白云质灰岩中。地层岩性控矿十分明显。② 矿床严格受东西向构造带中东安溪逆冲断裂带控制，为该矿床的导矿构造。次级的东西向、北东向断裂破碎带为矿体的容矿构造。③ 以W异常为主，异常面积大，强度高，Au异常面积小，异常强度中等。Hg、Sb异常背景值稍高，但异常微弱。④ 矿体矿石类型较单一，以网脉带型白钨矿为主，其次为石英脉型白钨矿。⑤ 硅化强烈和褪色化宽度大的部位，往往指示有富矿体存在。硅化和褪色化蚀变可以作为寻找钨矿的找矿标志。
5.3.2. 冷家溪金矿床地质特征
矿体(脉)主要赋存于新元古界板溪群马底驿组中上部紫红色板岩、砂质板岩中，在板溪群五强溪组中亦有产出。该矿区共发现大小含金石英脉75条，矿脉群(带)东西长4 km，宽1.0 km左右。矿体(脉)长10~800 m，矿脉间距数米到60 m不等，呈斜列雁行状。
矿石类型：主要有石英–白钨矿–黄铁矿型、自然金–黄铁矿–石英型和石英–碳酸盐型三种类型，其中自然金–黄铁矿–石英型为最主要的矿石类型。
根据矿物组构和矿体空间分布特征，该矿区成矿大体可分三个阶段：第一阶段石英白钨矿期，围岩蚀变主要有褪色化、硅化、黄铁矿化等；第二阶段为自然金硫化物期，蚀变以黄铁矿化为主；第三阶段为石英碳酸盐期，蚀变以碳酸盐化为主，与金矿化关系不甚明显。
结合上述矿床地质特征和图7可以看出，冷家溪型金矿成矿富集规律如下：① 矿体主要赋存于板溪群马底驿组中上部紫红色绢云母板岩、砂质板岩中，少数矿体产于板溪群五强溪组凝灰质板岩、凝灰岩中，受地层岩性控制明显。② 矿床受东西向冷家溪深大断裂带控制，为控矿构造，次级东西向断裂破碎带、节理密集带为储矿构造。③ Hg、Sb异常组合形态相似、位置相近，异常面积小，异常强度中等，是前缘晕的指示元素。W的背景值较高，但异常微弱。④ 矿体矿石以自然金–黄铁矿–石英型为主。⑤ 矿化作用与围岩褪色化、硅化、黄铁矿化关系密切，各种蚀变作用叠加时形成矿体。
5.3.3. 找矿靶区
Figure 7. Comprehensive anomalies of Au, Sb, Hg, and W and prospecting targets
图7. 研究区Au、Sb、Hg、W综合异常及找矿靶区
西安钨金矿区、冷家溪金矿区均已有探矿权，不再作为找矿靶区。综合上述异常特征，结合地层、岩性、构造等有利控矿因素评判，在其他研究区共圈出三个Au-Hg-Sb-W找矿靶区。
① 寺坪Au找矿靶区
该找矿靶区综合异常为以Au为主的综合异常，位于研究区北部，1:20万化探异常显示有较好的金元素异常，在寺坪附近浓集中心明显，与已知矿点吻合度好。在江里溪、寺坪、齐家冲等地见有多条金矿化带，金矿体主要产于板溪群马底驿组和五强溪组断裂破碎带中，另在冷家溪小木坪组中也见有金矿化带，岩性主要为板岩、粉砂质板岩、砂质板岩，呈脉状及透镜状产出。
② 雷打岭Au找矿靶区
该找矿靶区综合异常为以Au为主的综合异常，位于研究区东部，1:20万砂金重砂异常与水系沉积物异常重叠性好，异常强度高，浓集中心明显，与已知矿点吻合度好。在雷打岭等地见有多条金矿化带，金矿体产于板溪群马底驿组和五强溪组板岩、粉砂质板岩、砂质板岩的层间裂隙及断裂破碎带中，多呈脉状产出。
③ 木旺溪W-Au找矿靶区
该找矿靶区综合异常为以W为主的综合异常，位于茶庵铺镇，W异常规模大、强度高。在木旺溪等地见有多条白钨矿化带，产于冷家溪群中，岩性主要为板岩、粉砂质板岩、砂质板岩。
5.4. 异常实证
罗大富 [28] [29] 、王操 [30] 等在雷打岭金矿普查项目中调查发现，金矿体产于近东西向及北东向蚀变矿化带中，矿体主要由含金石英脉、含金构造角砾岩及含金蚀变破碎粉砂质板岩与粉砂岩组成，矿体的产出严格受断裂构造控制。矿石品位最大值2.91 g/t，最小值0.96 g/t，平均值1.81 g/t。矿石主要具碎裂结构、镶嵌结构；角砾状、块状、网脉状及蜂窝状构造。围岩蚀变主要见硅化、黄铁矿化、绿泥石化、绢云母化、碳酸盐化等，地表矿脉带中具较强的褐铁矿化，其中硅化、黄铁矿化与金矿化关系最为密切。蚀变强度由金矿体向外整体逐渐变弱，蚀变分带性不明显。调查结果认为雷打岭异常区勘探深度低，矿体具有半隐伏特征，初步预测矿区及周边深部找矿潜力巨大。
6. 结论
1) 研究区土壤地球化学数据聚类分析、因子分析以及元素成矿潜力研究结果表明，Au、Hg、Sb、W的成矿潜力最大；
2) 应用含量—面积分析法计算Au、Hg、Sb、W的异常下限值分别为40 ng/g、1 μg/g、10μg/g、20 μg/g；
3) 单元素异常和组合异常图分析显示，研究区金矿有两种类型，一种是冷家溪型金矿，主要特征是以Au异常为主，Hg、Sb异常吻合度高，为前缘晕指示元素，W异常低，矿产类型以金矿为主。另一种是西安型金矿，主要特征是以W异常为主，Au异常中等，Hg、Sb异常低，矿产类型以钨矿为主，金矿为辅；
4) 通过分析西安钨金矿区、冷家溪金矿区矿床地质特征，总结成矿富集规律，圈定了3个找矿靶区，普查结果发现雷打岭附近致矿异常较好，找矿潜力巨大。
致谢
本文在编写过程中得到了吉林大学地球探测科学与技术学院郝立波教授、赵新运博士和湖南省自然资源事务中心贺安生教授的指导与帮助，在此表示衷心的感谢!
基金项目
湖南省自然资源厅项目“湖南省桃源县地球化学与土地质量调查(20200805)”资助。
NOTES
*通讯作者。