胶西北金成矿区He、Ar同位素组成及成矿流体来源研究

王金辉

(山东省地质调查院，山东 济南 250014)

胶西北金成矿区位于胶东半岛的西北部，截至目前累计查明金资源储量超过 4 000 t，成为继南非兰德和乌兹别克斯坦穆龙套之后的世界第三大岩金成矿区(宋明春，2015)，金矿床主要沿三山岛断裂带、焦家断裂带和招平断裂带3条断裂带分布，成因类型主要有焦家式破碎带蚀变岩型和玲珑式石英脉型(李士先等，2007)。金矿主成矿期成矿时代为120 Ma左右(杨进辉等，2000； 李厚民等，2003； 宋明春等，2011)，与早白垩世岩石圈拉张减薄的峰期时代一致(翟明国等，2004； 宋明春等，2011)。

前人对胶西北金成矿区成矿热液来源进行了大量的深入研究，但在矿床成因方面存在着巨大的分歧，主要有以下几种观点：变质热液(沈保丰等，1994)、岩浆期后热液(李士先等，2007)、大气降水(卢焕章等，1999)、岩浆热液为主后期大气降水加入(林文蔚等，1998)、与中基性脉岩有关的幔源流体(范宏瑞等，2005)、地幔流体参与了金矿的成矿作用(杨进辉等，2000； 张连昌等，2002； 邓军等，2004； 毛景文等，2005)等。

稀有气体尤其是He、Ar具有化学性质不活泼及不同源区的同位素组成差异明显等特点,被广泛应用于金矿成矿流体来源及成矿机制研究(胡瑞忠，1997； 毛景文等，2000； 张连昌等，2002； 张运强等，2012； 杨猛等，2012； 薛建玲等，2013； 焦学尧等，2016； 何阳阳等，2017)。前人对焦家金矿床及招平金成矿带北段的玲珑金矿均进行了He、Ar同位素研究(张连昌等，2002； 毛景文等，2005)，但是按照3条断裂成矿带来系统对比研究的基本上没有，不同断裂带控制下的矿床的成矿流体是否有差异？其对成矿作用又有什么指示意义？这对于深入了解研究大区的成矿作用及流体来源有重要意义。本文通过对招平金成矿带中南段的姜家窑金矿床、曹家洼金矿床和近年来新发现的三山岛成矿带的三山岛北部海域金矿床、焦家深部的纱岭金矿床矿石样品的载金黄铁矿中流体包裹体进行He、Ar同位素测试，结合前人在该地区其它矿床取得的He、Ar同位素数据，在分析后期地质作用对成矿流体初始He、Ar同位素造成的影响上，对胶西北金成矿区3条主要成矿带内金矿成矿流体进行示踪，探讨成矿流体的来源差异。

1 区域地质矿产特征

胶东地区位于华北克拉通东南部，郯庐断裂之东侧。区内分布着众多超大、中型金矿床，组成了举世闻名的胶东金矿集中区。胶西北金成矿区位于胶东金矿集中区西北部，是其主要组成部分。区内金矿类型主要有焦家式破碎带蚀变岩型和玲珑式石英脉型，均根据胶西北金成矿区典型矿床命名。焦家式破碎带蚀变岩型金矿床主要分布在三山岛金成矿带、焦家金成矿带和招平金成矿带中南段，矿床受规模较大的区域性断裂三山岛断裂、焦家断裂和招平断裂控制。玲珑式石英脉型金矿主要分布于招平金成矿带北段，矿床受招平断裂伴生、派生断裂控制，倾角陡，矿体呈高角度倾斜。

胶西北地层岩体主要有前寒武纪变质岩系和中-新生代地质体组成。前寒武纪变质岩系包括新太古代TTG质片麻岩，少量变辉长岩(斜长角闪岩)、中太古代唐家庄岩群及古元古代荆山群，新太古代胶东岩群呈残留包体出现。中生代岩体有侏罗纪玲珑花岗岩、白垩纪郭家岭型花岗岩、伟德山型花岗岩(艾山岩体和南宿岩体)、崂山型花岗岩(大泽山岩体)和酸性、中-基性脉岩，其中以侏罗纪玲珑花岗岩为主。新生代地质体为第四纪松散沉积层，沿海边和河流分布(图1)。区域内NE向断裂发育，自西向东分别为三山岛断裂(F1)、焦家断裂(F2)、招平断裂(F3)，具有多期、多阶段活动性，控制着蚀变带及金矿床的空间定位，金矿床严格受三条断裂及其次级断裂控制，沿断裂呈串珠状分布。

三山岛断裂位于陆地上北起三山岛，南至仓上-潘家屋子一带，向南、北两端均延伸入渤海，陆地长度为12 km，宽20～400 m不等，走向40°～50°，倾向SE，倾角30°～40°，局部达80°。断裂在平面上呈“S”形延伸，形态不规则。自南向北控制着仓上、新立、三山岛、西岭、北部海域等金矿床，皆达大型或特大型规模，累计查明金资源储量大于1 200 t(于学峰等，2016)。

焦家断裂南端起始于莱州市紫罗姬家，北到龙口市姚家，控制长度约60 km，宽度50～500 m，倾向NW，倾角30°～50°，局部达78°。断裂在平面上呈“S”形展布，形态不规则，局部膨缩，断裂下盘发育与走向平行或呈“入”字形相交的次级断裂。焦家、新城、寺庄、望儿山、河西等金矿床均受焦家断裂或其下盘次级断裂控制。近年在深部发现纱岭、腾家、朱郭李家、招贤等大型-特大型矿床。目前已探获金资源储量1 400 t(于学峰等，2016)。

图 1 胶西北地区地质简图(据宋明春，2015)Fig. 1 Simplified geological map of northwestern Jiaodong area(after Song Mingchun, 2015)1—第四系； 2—郭家岭型花岗岩； 3—玲珑型花岗岩； 4—新太古代变质岩系； 5—断裂； 6—金矿位置； F1—三山岛断裂； F2—焦家断裂；F3—招平断裂1—Quaternary； 2—Guojialing type granite； 3—Linglong type granite； 4—Neoarchean metamorphic rock complex； 5—fault； 6—gold position； F1—Sanshandao fault； F2—Jiaojia fault； F3—Zhaoping fault

招平断裂是区内规模最大的控矿断裂，南起平度宋格庄，向北经夏甸、大尹格庄、招远城区，抵玲珑矿田，由玲珑向北转为NE向至龙口市颜家沟，长约120 km，总体呈NNE向“S”形展布。按展布方向及矿化特征分为北、中、南三段，北段颜家沟-丁家庄子段(又称破头青断裂)，走向50°～70°，带宽40～300 m，最宽800 m；中段招远城-夏甸段，走向10°～20°，宽30～60 m，最宽80 m；南段为山后-山旺段(刘述敏等，2016)。招平断裂带北段及其下盘的次级断裂控制着玲珑金矿田，有玲珑、台上、东风、水旺庄、阜山、滦家河等大型-特大型金矿床；中段分布有大尹格庄-夏甸金矿田，有大尹格庄、夏甸、姜家窑、曹家洼、山后等特大型-中型金矿床；南段主断裂上矿化较弱，仅有南墅、下庄、山旺等小型金矿，下盘的次级断裂控制着旧店金矿田。成矿带累计查明金资源储量大于1 100 t(于学峰等，2016)。

2 样品及分析方法

本次研究用于He-Ar同位素测试的黄铁矿样品分3个断裂带采集。三山岛成矿带样品来自于北部海域金矿床、三山岛金矿床、仓上金矿床；焦家成矿带样品来自于纱岭金矿床、焦家金矿床；招平成矿带样品来自于姜家窑金矿床、曹家洼金矿床、玲珑金矿床、金翅岭金矿床。其中北部海域金矿床、沙岭金矿床、姜家窑金矿床及曹家洼金矿床数据为本次分析，其余为搜集前人测试数据。

研究样品在河北区域地质调查所实验测试中心经破碎、筛选后，在双目镜下挑选40～60目的新鲜、晶形完好黄铁矿颗粒，纯度达99%以上。

He-Ar同位素测试由中国地质科学院矿产资源研究所完成。采用Helix SFT稀有气体质谱仪测试，系统由压碎、纯化和质谱系统组成。测试在高真空下完成，压碎和纯化系统真空n×10-7Pa，质谱系统真空在n×10-8Pa。质谱离子源采用Nier，灵敏度对He在800 μA阱电流时好于2×10-4amps/Torr，对Ar在200 μA阱电流时好于1×10-3amps/Torr。40Ar静态上升率小于1×10-12cm3STP/min，36Ar本底小于5×10-14cm3STP。法拉第杯分辨率>400，离子计数器分辨率>700，可将3He与4He、HD+H3+与3He峰完全分开。详细的测试流程及说明参见叶先仁等(2001)。

H、O同位素分析测试由中科院地质与地球物理研究所稳定同位素地球化学实验室测试完成。硅酸盐的氧同位素分析分析方法为：称取200目干燥的硅酸盐或氧化物或全岩样品10～30 mg，在真空条件下，于550～700℃与纯净的BrF5恒温反应释放出氧气。用组合冷阱分离生成的SiF4、BrF3等杂质组分。令纯化后的氧气在700℃铂催化作用下，与碳棒逐级反应生成二氧化碳(CO2)。将纯净的二氧化碳导入Finnigan Delta S型质谱仪中测定其同位素比值。分析测试误差好于0.2‰。包体中氢同位素分析方法为：在高真空约10-3Pa条件下，通过加热(温度为600℃)将包体充分爆裂,将水之外的杂气去除。在真空中用干冰将包体中水冻住并引入金属铬反应器中。在温度800℃的条件下，包体水进入铬反应炉转变为氢气。将氢气导入Finnigan MAT 252型质谱仪中进行氢同位素的测量。分析测试误差好于3‰。

3 测试结果与讨论

3.1 分析结果

将本次研究测试的10件和前人的31件金矿床黄铁矿流体包裹体He、Ar同位素组成及相关计算参数一并列入表1。

从表1中可以发现4He含量变化较大，为0.12×10-7～10.68×10-7cm3STP/g，3He/4He值为0.09～7.14 Ra，平均为1.54 Ra。三山岛带4He含量从极少到8.06×10-7cm3STP/g，3He/4He值为1.04～7.14 Ra，平均为1.95 Ra；焦家带4He含量为2.33×10-7～10.68×10-7cm3STP/g，3He/4He值为0.76～2.36 Ra，平均为1.74 Ra；招平带4He含量为1.62×10-7～10.65×10-7cm3STP/g，3He/4He值为0.09～1.98 Ra，平均为1.17 Ra。40Ar的变化范围为0.51×10-7～25.34×10-7cm3STP/g，其中三山岛带变化范围0.51×10-7～23.14×10-7cm3STP/g；焦家带变化范围2.94×10-7～25.34×10-7cm3STP/g；招平带变化范围1.35×10-7～16.99×10-7cm3STP/g。40Ar/36Ar值变化范围为284～4 042.6，平均为848.6。其中三山岛带变化范围284～2 295，平均为614.63；焦家带变化范围350～2 955.3，平均为1 117.35；招平带变化范围321～4 042.6，平均值840.43。40ArE/4He值为0.029 6～1.89，其中三山岛带变化范围0.36～1.89，平均为0.88；焦家带变化范围0.09～1.15，平均为0.59；招平带变化范围0.029 6～0.960 2，平均为0.56。

3.2 影响因素分析

胶西北金成矿区中，黄铁矿为金矿主要载金矿物，黄铁矿相比较其他矿物不易发生漏气，因此，其是流体包裹体He、Ar同位素研究最理想的测试对象，并能反映成矿期成矿流体的原始信息(Stuartetal.，1995； 胡瑞忠等，1997b，1999； Huetal.，1998)。一些后期作用可能对流体包裹体初始He、Ar同位素组成有影响，包括： ① 成矿后热液改造，② 矿物晶格及流体包裹体内He、Ar扩散或新生放射性成因He、Ar加入，③ 样品分析过程中同位素分馏及大气污染(胡瑞忠，1997)。

对于第①种情况，本次所测试的黄铁矿晶形完好，未变形，无明显的后期改造，说明黄铁矿流体包裹体中的He在流体包裹体被圈闭后无明显的丢失(Baptiste and Fouquet，1996)。对于第②种情况，胶西北成矿区内黄铁矿为不含钾矿物，U、Th含量也极低，故放射成因40Ar的含量可以忽略不计(Qiu and Groves，1999)；另外放射成因的40Ar也难以进入流体包裹体中，晶格中由放射性元素衰变产生的4He很难进入包裹体(胡瑞忠等，1998； 薛春纪等，2003)；胶西北成矿区岩石中缺乏含锂矿物，由其产生的放射性成因3He的后生叠加影响可排除(胡瑞忠等，1998)。黄铁矿晶格及流体包裹体内He、Ar的扩散系数均很低(Burnardetal.，1999)，以黄铁矿为载体的流体包裹体中He、Ar的扩散丢失对He、Ar同位素组成影响也很小，甚至可忽略不计(胡瑞忠，1997)。对于第③种情况，由于宇宙成因的3He只产生在近地表1.5 m的范围内(胡瑞忠，1997； 丁德建等，2014)，本文的样品均采自地下坑道和钻孔中，且暴露时间很短，因此可以排除样品中存在宇宙成因3He的可能(Stuartetal.，1995)。氦在大气中的含量极低，不足以对地壳流体中氦的含量和同位素组成产生明显的影响。另样品中若含有大气氦组分，则F4He=1(Kendricketal.，2001； 毛景文等，2005)。由表1中可以发现，胶西北金矿床的F4He范围为53～72 659，远远大于1，所以大气He的混染作用可以排除。本次测得的黄铁矿包裹体的He、Ar同位素数据，可以代表原生流体包裹体的初始值。成矿流体中的氦只有地幔和地壳两个来源。另外由于稀有气体的特殊性，无论在成矿时捕获还是真空压碎法取得流体包裹体样品，均不会产生明显的同位素分馏(胡瑞忠等，1999)。

表 1 胶西北地区金矿黄铁矿中流体包裹体的He、Ar同位素组成
Table 1 He, Ar isotopic composition of fluid inclusions in pyrite from the gold deposits in northwestern Jiaodong area

注： Ra为空气的3He/4He比值(取1.40×10-6)； R/Ra表示样品的3He/4He与空气的3He/4He(取1.40×10-6)比值； He(%)为地幔成因的He同位素，He(%)=[(R-Rc)/(Rm-Rc)]×100，Rm、 Rc、 R分别代表地慢流体、 地壳流体以及样品的氦同位素组成； 40Ar*(%)表示成矿流体中放射成因40Ar的相对含量，40Ar*(%)=[1-295.5/(40Ar/36Ar)样品]×100； 40ArE表示扣除空气40Ar后的过剩氩，40ArE=(40Ar)样品-295.5×36Ar样品； F4He=(4He/36Ar)样品/(4He/36Ar)大气，大气中的4He/36Ar值为0.165 5。

3.3 He和Ar同位素分析结果讨论

热液流体中He、Ar同位素有4种来源，它们分别是大气来源、饱和大气水来源、地幔来源和地壳来源(Burnardetal.，1999； 李晓峰等，2003)。其中，饱和大气水中3He/4He值为1.39×10-6，40Ar/36Ar值为295.5；幔源流体3He/4He值为6～9 Ra，40ArE/4He值为0.33～0.56，40Ar/36Ar值大于4 000；地壳流体的3He/4He值为0.01～0.05 Ra，40ArE/4He值为0.156；大气Ar的同位素40Ar/36Ar值为295.5(Dunai and Touret, 1995； Stuartetal.，1995)。由于He在大气中的含量极低，对地壳流体中的He丰度和同位素组成不会产生明显的影响(Martyetal.，1989)。

本次研究测试和收集的41件样品中，除三山岛金矿1件样品由于He含量低于检出限，1件样品出现极高值外，其余39件黄铁矿样品的3He/4He值变化范围为0.09～2.36 Ra，平均为1.40 Ra，高于地壳3He/4He值(0.01～0.05 Ra)，低于地幔3He/4He值(6～9 Ra)(Stuartetal.，1995)，说明成矿流体的He同位素是地幔和地壳两种来源的混合。三山岛成矿带11件样品3He/4He均值为1.95 Ra，焦家成矿带11件样品3He/4He均值为1.74 Ra，招平成矿带18件样品3He/4He均值为1.17 Ra，从西往东逐渐降低，说明地幔流体参与成矿自西往东逐渐降低。在氦同位素组成分布图(图2)上，数据点均落在地幔与地壳之间的过渡区。

图 2 胶西北地区金矿黄铁矿流体的He同位素组成(据Mamyrin and Tolstikhin，1984)Fig. 2 He isotope composition of fluid inclusions in pyrite from gold deposits in northwestern Jiaodong area (after Mamyrin and Tolstikhin, 1984)

与3He/4He比值相比，黄铁矿流体包裹体中的40Ar/36Ar比值为284～4 042.6，分布范围较广，平均848.6，除仓上Csh5样品分析结果外都大于大气Ar同位素均值。40Ar/36Ar值高于大气值说明流体中存在来源于地壳的放射成因Ar(40Ar*)，经计算获得放射性成因Ar所占比例范围为1.17%～92.69%。其中三山岛成矿带1.17%～87.12%，平均33.04%，焦家成矿带15.57%～90.00%，平均55.75%，招平成矿带7.94%～92.69%，平均44.23%，各成矿带内及成矿区不同样品放射性成因Ar所占比例变化较大。

胶西北成矿区具有较高的3He/4He和40Ar/36Ar比值，说明胶西北成矿区地幔流体参与了部分成矿。同时40Ar/36Ar比值变化范围较大，说明了深部来源的幔源流体在成矿作用过程中，少量加入了浅部低温海水或者大气降水。

流体幔源He所占比例范围在0.97%～29.87%，平均17.57%，显示出幔源流体对成矿有贡献。其中三山岛成矿带平均24.65%，焦家成矿带平均22.01%，招平成矿带平均14.62%，说明胶西北金矿床成矿物质来源主要来源于地壳，且有少量幔源物质的加入，自西往东幔源物质加入逐渐减少。

将He和Ar同位素组成投入R/Ra-40Ar/36Ar图解(图3)，投点集中在地壳流体与地幔流体之间，并且靠近大气饱和水，说明成矿流体以地壳流体与地幔流体参与成矿作用外,大量的大气降水也参与其中。另外从图3可以发现，除去个别样品外，基本上所有数据点也沿平行x坐标轴的方向分布，表明成矿流体中Ar同位素的变化基本上不受氦的影响，说明胶西北成矿区内成矿流体中He和Ar同位素的组成变化相互独立，指示了区内金矿床具有相似的流体来源及混合过程，只是在后期的演化中由于成矿构造环境的差异而导致不同。

图 3 胶西北地区金矿成矿流体R/Ra-40Ar/36Ar图解(据宋世明等，2007)Fig. 3 R/Ra versus 40Ar/36Ar diagram of ore-forming fluids in northwestern Jiaodong area (after Song Shiming et al., 2007)

40ArE/4He比值除2件较低(0.029 6，0.09)外，其余样品比值范围0.232 2～1.89，平均0.65，其中三山岛带平均0.88，焦家带平均0.59，招平带平均0.56，普遍大于地壳值(0.156)，与地幔值(0.33～0.56)(Stuartetal., 1995；胡瑞忠等，1997a)相近并略高，成矿流体中40ArE/4He值高于正常40Ar/4He值，说明大气降水在流经地壳岩石过程中捕获了放射性成因40Ar。由于黄铁矿中氦的封闭温度低(200℃)，而氩的封闭温度高(>200℃)(Ballentineetal.，2002)，高的40ArE/4He值表明成矿流体温度较高(>200℃)，这与该区主成矿阶段温度为170～335℃，一般高于200℃的事实一致(范宏瑞等，2005)。同时也反映出胶西北成矿区成矿流体主要来源于地壳，并与深部幔源流体发生不同程度的壳幔混合，大气降水也参与成矿作用。将He和Ar同位素组成投入R/Ra-40ArE/4He图解(图4)，投影点落在地幔流体和地壳流体之间斑岩铜矿床区域左上方，显示成矿流体中有地幔流体不同比例的混入。

3.4 与成矿流体δ18O和δD组成的比较

为了进一步对成矿流体来源的合理解释，分析了三山岛北部海域、姜家窑及曹家洼金矿床成矿流体的H、O同位素组成(表2)。搜集了前人部分资料，综合地质和成矿流体He-Ar同位素系统进行综合分析。

测试获得胶西北金成矿区δ18O范围为7.9‰～13.0‰，三山岛带均值12.08‰，焦家带均值10.37‰，招平带均值9.68‰。δ18OH2O范围1.57‰～8.55‰，三山岛带均值7.26‰，焦家带均值5.43‰，招平带均值3.63‰，自三山岛带往东至招平带逐渐降低,说明幔源物质参与逐渐降低。氢同位素组成δD值变化于48‰～-86‰，三山岛带均值-60.10‰，焦家带均值-76.67‰，招平带均值-66.84‰。将上述数据投影到氢氧同位素关系图中(图5)，可以看出，投点主要落于变质水下方，岩浆水附近，自西部三山岛成矿带向东部招平成矿带逐渐靠近大气降水。δD分布范围相对较窄，表明成矿流体源于初始岩浆水，晚期有少量大气降水的参与，这与前人的认识一致(毛景文等，2005)。由此可见，胶西北成矿区金矿床成矿流体H-O和He-Ar同位素系统所示踪的流体来源具有一致性。

图 4 胶西北地区金矿成矿流体R/Ra-40ArE/4He图解(底图Ballentine et al.，2002)Fig. 4 R/Ra versus 40ArE/4He diagram of ore-forming fluids in northwestern Jiaodong area (after Ballentineet al., 2002)

3.5 成矿作用过程

依据前人研究成果和本文得出的结论，推断胶西北成矿区成矿流体的形成、演化经历了如下过程：在侏罗纪，古太平洋板块自SE向NW向欧亚大陆下方俯冲，胶东地壳增厚并发生地壳重熔，在160～150 Ma，玲珑型花岗岩侵位，古老陆壳中的金初步富集。在白垩纪早期，中国东部构造体制发生转折，岩石圈强烈减薄，地幔物质上涌，构造、岩浆活动活跃，壳-幔发生相互作用。130～126 Ma，壳幔混源的郭家岭型花岗岩侵位。同时，受早白垩世壳幔混合型岩浆作用影响，玲珑花岗岩岩基隆升，其上前寒武纪变质岩发生拆离滑脱，形成发育在玲珑岩基顶部的三山岛断裂、焦家断裂和招平断裂等拆离断层。在120 Ma左右，岩石圈已经大幅度减薄，大量煌斑岩、辉绿玢岩、闪长玢岩等中基性脉岩沿NNE断裂侵位，地幔物质(包括流体)不断向地壳涌入，地幔流体在向上运移的过程中，与地壳流体发生混合，并有大气降水的加入。碳、氧、氢、硫同位素和流体包裹体研究反映它们由来自同一的流体库，通过水岩反应和流体混合作用，金最终卸载成矿。

表 2 胶西北金矿石英流体包裹体氢氧同位素组成
Table 2 H and O isotope composition of fluid inclusions in quartz from gold deposits, northwestern Jiaodong

注： t/℃依据石英流体包裹体测温结果；通过石英-水之间的氧同位素分馏方程： 1 000 lnαQ-W=3.38×106/t2-2.90(Clayton et al.，1972)，估算成矿热液中水的氧同位素组成δ18OH2O，即δ18OH2O=δ18Ov-smow-1 000 lnαQ-W，其中αQ-W为石英与水之间的分馏系数。

图 5 胶西北地区金矿床成矿流体δD-δ18OH2O分布图(底图据Hugh and Taylor, 1974)Fig. 5 δD-δ18OH2O of ore-forming fluid in northwestern Jiaodong area (after Hugh and Taylor, 1974)

4 结论

(1) 胶西北成矿区3条金成矿带具有相似的流体来源及混合过程，成矿流体主要来源于地壳，并与深部幔源流体发生不同程度的壳幔混合，大气降水也参与成矿作用。后期受不同构造环境的影响和制约，从而在不同的构造部位形成不同的金矿床。

(2) 黄铁矿流体包裹体揭示三山岛成矿带3He/4He均值为1.95 Ra，焦家成矿带3He/4He均值为1.74 Ra，招平成矿带3He/4He均值为1.54 Ra，从西往东逐渐降低，说明地幔流体参与成矿自西往东逐渐降低，与H、O稳定同位素研究得出的结论一致。

(3) 黄铁矿流体包裹体中放射性成因Ar所占比例范围三山岛成矿带为1.17%～87.12%，平均33.04%，焦家成矿带为15.57%～90.00%，平均55.75%，招平成矿带为7.94%～92.69%，平均44.23%，说明在成矿过程中，地壳成矿流体参与成矿作用焦家成矿带最高，招平成矿带次之，三山岛成矿带最低。

致谢本文在成文过程中得到了审稿专家的学术指导和中肯建议，在此表示诚挚的谢意。