• 中文核心期刊
  • 中国科技核心期刊
  • 中国科学引文数据库核心期刊

大兴安岭北段古生代扎兰屯岛弧增生演化的锆石U-Pb年龄证据

汪岩, 秦涛, 张立东, 臧延庆, 钱程, 赵春强

汪岩, 秦涛, 张立东, 臧延庆, 钱程, 赵春强. 2021: 大兴安岭北段古生代扎兰屯岛弧增生演化的锆石U-Pb年龄证据. 地质通报, 40(12): 2186-2188. DOI: 10.12097/gbc.dztb-40-12-2186
引用本文: 汪岩, 秦涛, 张立东, 臧延庆, 钱程, 赵春强. 2021: 大兴安岭北段古生代扎兰屯岛弧增生演化的锆石U-Pb年龄证据. 地质通报, 40(12): 2186-2188. DOI: 10.12097/gbc.dztb-40-12-2186
WANG Yan, QIN Tao, ZHANG Lidong, ZANG Yanqing, QIAN Cheng, ZHAO Chunqiang. 2021: Zircon U-Pb age evidence for the accretionary evolution of the Paleozoic Zalantun island arc in the northern part of the Great Xing'an Range. Geological Bulletin of China, 40(12): 2186-2188. DOI: 10.12097/gbc.dztb-40-12-2186
Citation: WANG Yan, QIN Tao, ZHANG Lidong, ZANG Yanqing, QIAN Cheng, ZHAO Chunqiang. 2021: Zircon U-Pb age evidence for the accretionary evolution of the Paleozoic Zalantun island arc in the northern part of the Great Xing'an Range. Geological Bulletin of China, 40(12): 2186-2188. DOI: 10.12097/gbc.dztb-40-12-2186

大兴安岭北段古生代扎兰屯岛弧增生演化的锆石U-Pb年龄证据

基金项目: 

中国地质调查局项目《松辽盆地西缘扎赉特—巴林左旗区域地质调查》 DD20190039

详细信息
    作者简介:

    汪岩(1968-), 男, 教授级高级工程师, 从事区域地质与火成岩研究。E-mail: wy68413@163.com

    通讯作者:

    秦涛(1985-), 男, 硕士, 高级工程师, 从事区域地质与火成岩研究。E-mail: 465468905@qq.com

Zircon U-Pb age evidence for the accretionary evolution of the Paleozoic Zalantun island arc in the northern part of the Great Xing'an Range

  • 柴达木盆地北缘(柴北缘)构造带处于南祁连地块与柴达木地块的拼合部位[1-2],研究表明其为构造复杂、物质组成多样、时间跨度大的多单元复合构造带[3],由南向北主要由鱼卡河-沙柳河高压-超高压变质带、滩间山群火山-沉积岩系、欧龙布鲁克地块等组成[2-4]。随着研究程度的提高,绿梁山地区发现早古生代弧后盆地型蛇绿混杂岩[2, 5]及原岩为大洋蛇绿岩的榴辉岩[6-9],说明早古生代蛇绿岩也是柴北缘构造带的重要组成部分。郝国杰等[10]指出,在古生代之前的南华纪—震旦纪,受Rodinia超级大陆裂解事件的影响,柴北缘发生了拉张与裂解,在乌兰以北形成具洋壳特点的蛇绿岩,代表了初始洋盆的形成;沿柴北缘构造带分布的一系列形成于消减带环境下的岩石[2, 11-13]表明,柴北缘古大洋在晚寒武世(约515Ma)已经开始俯冲消减,并至少持续至中奥陶世晚期(约460Ma);吴才来等[12, 14-16]通过对团鱼山、柴达木山等地同造山陆-陆碰撞成因的S型花岗岩的研究,认为祁连洋于晚奥陶世闭合,祁连地块向南逆冲到柴达木地块之上并形成陆陆碰撞带;对区内大量高压/超高压变质岩的研究表明,柴北缘在志留纪(440~420Ma)主要处于陆-陆碰撞和深俯冲阶段[17-19];Song等[20]认为,晚志留世—早泥盆世大陆深俯冲引起先期大洋地壳折返并发生减压熔融,在柴北缘形成TTG岩浆;而具典型伸展型磨拉石建造的牦牛山组的形成是柴北缘早古生代主造山作用结束的标志[21]

    基性岩脉是由源于地幔的玄武质岩浆及其分异的和受地壳混染作用影响的岩浆充填张性空间而成,是岩石圈伸展和地壳拉张的产物[22-23],是一种特殊的构造岩浆类型,在大陆地壳演化中具有特殊而重要的研究意义[24-25]。近年来笔者等在柴达木盆地西北缘开展1:5万区域地质调查时,在柴北缘构造带西端欧龙布鲁克地块西北缘达肯大坂岩群中识别出大量基性岩脉,然而这些基性岩脉的形成时代、源区性质、构造环境及对柴北缘构造演化的意义尚不清楚。鉴于此,本次选择欧龙布鲁克地块西北缘达肯大坂岩群中的辉长岩脉为研究对象,通过岩石学、地球化学、年代学及Lu-Hf同位素特征研究,确定辉长岩脉的形成时代、源区性质、构造环境,并结合区域地质背景,为柴北缘构造演化的研究提供约束。

    研究区位于阿尔金断裂以南,柴北缘构造带西端欧龙布鲁克地块的西北缘,其北为红柳沟-拉配泉蛇绿构造混杂岩,南部被第四系冲洪积物覆盖(图 1-a)。区内出露欧龙布鲁克陆块基底岩系古元古代达肯大坂岩群,该岩群主要由混合片麻岩段(Pt1DK1.)、条带状片麻岩段(Pt1DK2.)、片麻岩夹片岩段(Pt1DK3.)、片岩段(Pt1DK4.)、大理岩段(Pt1DK5.)等组成[27]。达肯大坂岩群发育大量基性和中酸性岩脉,基性岩脉以(变)辉长岩脉为主,宽窄不一(1~ 5m)、延伸不远(5~50m),走向多为北东东向,局部呈北西—南东向延伸;中酸性岩脉主要为闪长岩、二长花岗岩、浅色花岗岩、花岗伟晶岩等(图 1-b),其中以花岗伟晶岩脉分布最广,但其分布无规律,一般长20~200m,宽2~5cm。本文研究的辉长岩脉多呈脉状、透镜状沿片理、片麻理侵入达肯大坂岩群片麻岩中,呈明显的侵入接触关系(图 2-a)。

    图  1  柴北缘地质简图(a,底图据参考文献[26]修改)及研究区地质图(b)
    1—达肯大坂岩群;2—大理岩透镜体;3—晚奥陶世花岗闪长岩;4—早二叠世英云闪长岩;5—中二叠世花岗闪长岩;6—中更新统冰碛堆积;7—上更新统冲洪积堆积;8—全更新统洪冲积堆积;9—基性-超基性岩脉;10—伟晶岩脉;11—辉长岩脉;12—花岗岩脉;13—韧性剪切带;14—一般断裂;15—隐伏断裂;16—样品位置
    Figure  1.  Sketch geological map of the northern margin of Qaidam Basin (a) and geological map of the study area(b)
    图  2  辉长岩岩相学特征
    a—辉长岩脉野外露头;b—角闪石、辉石等暗色矿物定向排列(正交偏光);c—残留辉石(单偏光);d—变余辉长结构(正交偏光);Hbl—角闪石;Px—辉石;Pl-斜长石
    Figure  2.  Petrologic features of metagabbro

    本次用于同位素及地球化学研究的辉长岩样品(编号为15-DK)采自冷湖镇北西82km处,地理坐标北纬39°07′42"、东经92°32′16"。样品风化色呈灰褐色,新鲜面绿黑色,具块状、弱片麻状构造,粒状变晶结构。岩石主要由角闪石(约50%)、斜长石(约38%)、辉石(约10%)、石英(< 1%)、榍石(< 1%),以及少量绢云母(< 1%)等组成。角闪石矿物晶体多呈粒状、粒柱状,颗粒长轴大致平行分布,多数由辉石变质形成(图 2-b),粒径大小一般为0.5~2.5mm。镜下角闪石呈绿色,具多色性,C∧Ng’≈24°,属普通角闪石。斜长石多为粒状或不规则粒状(图 2-b),粒径0.2~1.5mm,部分发育绢云母化、黝帘石化。辉石多数已发生角闪石化,残留辉石呈港湾状(图 2-c),局部可见变余辉长结构(图 2-d),单偏光镜下高突起明显,两组节理发育。

    样品的主量、微量及稀土元素测试分析在中国地质调查局西安地质调查中心实验测试中心完成,其中主量元素采用SX45型X荧光光谱仪(XRF)分析,分析误差小于1%;微量和稀土元素利用SX50型电感耦合等离子体光谱仪(ICP-MS)测定,分析误差小于5%~10%。锆石挑选在河北廊坊诚信地质服务有限公司完成,锆石制靶及反射光、阴极发光照相在自然资源部岩浆作用成矿与找矿重点实验室完成,测试点的选取首先根据锆石反射光和透射光照片进行初选,再与阴极发光图像反复对比,力求避开内部裂隙和包裹体,以获得较准确的年龄信息。LA-ICP-MS锆石微区U-Pb测年在自然资源部岩浆作用成矿与找矿重点实验室完成,采用193nmArF准分子(excimer)激光器的Geo Las200M剥蚀系统,ICP-MS为Agilent 7700,激光束斑直径24μm,以GJ-1为同位素监控标样,91500为年龄标定标样,NIST610为元素含量标样进行校正,普通铅校正依据实测204Pb进行校正。

    采用Glitter(ver4.0,Macquarie University)程序对锆石的同位素比值及元素含量进行计算,并按照Anderson的方法[28],用LAMICPMS Common Lead Correction(ver3.15)对其进行了普通铅校正,年龄计算及谐和图采用Isoplot(ver3.0)完成[29]

    锆石原位Lu- Hf同位素分析在配备了Geolas2500激光剥蚀系统的Nu Plasma HR多接收电感耦合等离子体质谱仪(MC-ICP-MS)上完成,激光剥蚀脉冲频率为10Hz,激光束斑直径为44μm,剥蚀时间约50s。用176Lu/175Lu = 0.02669和176Yb/172Yb=0.5886进行同量异位干扰校正计算,测定样品的176Lu/177Hf和176Hf/177Hf值[30]ɛHft)值计算采用176Lu衰变常数为1. 867×10-11a[31],球粒陨石现今的176Hf/177Hf=0.282785、176Lu/177Hf=0.0336[32];Hf亏损地幔模式年龄(tDM1)计算采用现今的亏损地幔176Hf/177Hf=0.28325和176Lu/177Hf=0.0384值[33]

    达肯大坂岩群中辉长岩脉的主量元素分析结果见表 1。岩石SiO2含量为50.86%~52.03%,平均51.22%;Al2O3、CaO含量较高,平均含量分别为16.00%、10.92%;低P2O5、TiO2,二者含量分别为0.07%~0.12%、0.81%~0.95%。在Nb/Y-Zr/TiO2 × 0.0001图解中,样品点均落入玄武岩区域(图 3)。全碱(K2O+Na2O)含量为3.86%~4.26%,Na2O/K2O值介于2.91~4.00之间,TFeO含量为8.66%~9.33%,MgO含量为5.10%~6.48%,在AFM图解中,样品点均落入钙碱性玄武岩范围(图 4)。Mg#[Mg#=100Mg/(Mg+Fe2+)]=55.0~59.9,小于原始岩浆玄武岩Mg#值(68~78),表明岩浆在演化过程中经历了一定程度的橄榄石、辉石等结晶分异作用[34]

    表  1  辉长岩主量、微量及稀土元素含量分析结果
    Table  1.  Major, trace elements and REE compositions of metagabbro
    下载: 导出CSV 
    | 显示表格
    图  3  辉长岩Nb/Y-Zr/TiO2×0.0001图
    Figure  3.  Nb/Y-Zr/TiO2×0.0001 diagram of metagabbro
    图  4  辉长岩AFM图解
    Figure  4.  AFM diagram of metagabbro

    辉长岩稀土和微量元素分析结果见表 1。样品稀土元素总量(ΣREE)较低,为43.17 × 10-6~ 57.53×10-6,为球粒陨石的13~17倍;轻稀土元素LREE=28.14 × 10-6~40.74 × 10-6,重稀土元素HREE=13.77×10-6~20.21×10-6,LREE/HREE=1.82~2.77,平均2.29,(La/Yb)N=1.07~2.05,表明轻稀土元素富集,暗示辉长岩原始岩浆演化分异程度较低[35];(La/Sm)N=1.16~2.07,(Gd/Yb)N=0.86~1.06,表明样品轻、重稀土元素内部分馏较弱。稀土元素分布模式曲线(图 5-a)整体表现为右倾型。δEu=0.80~ 0.83,显示弱的负异常,表明岩石演化过程中发生了弱的斜长石分离结晶作用[35],此外,δCe=0.98~ 1.02,δCe的波动范围较小,暗示样品具有一致的源区和相似的演化过程[36]。样品微量元素原始地幔标准化蛛网图(图 5-b)显示,所有样品均具有相似的配分型式,相对富集Rb、Th、Nd等元素,亏损Ba、Nb、Zr、Sr、P、Ti等元素。

    图  5  辉长岩稀土元素球粒陨石标准化图解(a)和微量元素原始地幔标准化蛛网图(b)[37]
    Figure  5.  Chondrite-normalized REE patterns (a) and primitive-mantle normalized spidergrams (b) for metagabbro

    对达肯大坂岩群辉长岩进行锆石阴极发光(CL)图像观察和LA-ICP-MS锆石U-Pb定年。锆石阴极发光图像见图 6,U-Pb同位素比值和表面年龄测试数据列于表 2。大多数锆石镜下呈无色透明,长柱状、短柱状晶形明显(长50~250μm),长宽比为1:1~3:1,CL图像显示该类锆石发育宽缓或条带状振荡环带结构(1、4、10、13、20号等),部分锆石具有后期蚀变或变质产生的窄浅色增生边(2、8、35号),显示岩浆结晶锆石的特征。本次选择具有代表性的39粒锆石进行LA-ICP-MS测年,分析结果显示数据点均沿谐和线分布,部分数据点位置偏下(图 7-a)。206Pb/238U表面年龄介于348±4~368±4Ma之间。39个测点谐和度高且年龄值较集中,206Pb/238U年龄加权平均值为357±4Ma,代表了辉长岩的形成年龄(图 7-b)。

    图  6  辉长岩代表性锆石阴极发光图像
    (圆圈表示测点位置,圆圈中的数字表示测点号)
    Figure  6.  Representative cathodoluminescence images of the zircons from metagabbro
    表  2  辉长岩LA-ICP-MS锆石U-Th-Pb同位素测年结果
    Table  2.  LA-ICP-MS zircon U-Th-Pb isotopic analyses of metagabbro
    下载: 导出CSV 
    | 显示表格
    图  7  辉长岩锆石U-Pb谐和图(a)和206Pb/238U年龄加权平均值(b)
    Figure  7.  U-Pb concordia diagram (a) and weighted mean 206Pb/238U ages (b) of zircons from metagabbro

    在U-Pb测年的基础上, 对达肯大坂岩群辉长岩的39颗锆石进行原位Hf同位素分析,其中有效测点34个。所有锆石的176Lu/177Hf值均小于0.002(表 3),表明因Lu衰变形成的放射性成因Hf同位素的积累量极少,Hf同位素比值可代表锆石形成时的初始Hf同位素比值。176Hf/177Hf值为0.282813~ 0.282979,各单颗粒锆石U-Pb年龄计算的εHft)值介于9.29~15.16之间,平均为12.25,在年龄-εHft)图解(图 8)上,所有数据点均位于亏损地幔演化线附近或之下。测试数据显示,达肯大坂岩群辉长岩Hf模式年龄(TDM)分布在378~610Ma之间,平均为494Ma。

    表  3  变辉长岩锆石Hf同位素组成
    Table  3.  Zircon Hf isotope compositions of metagabbro
    下载: 导出CSV 
    | 显示表格
    图  8  辉长岩锆石年龄-εHft)图解(底图据参考文献[45])
    Figure  8.  t-εHf(t) diagram for zircon of metagabbro

    蚀变作用和变质作用通常导致大离子亲石元素如K、Rb、Sr、Ba、Cs、Pb2+、Eu2+,具有明显的活动性,而稀土元素及部分高场强元素(Nb、Ta、Zr、Hf、Th、REE、Ce、U、Ti)甚至在高级变质作用中亦能相对稳定[38-40]。所以,本文主要利用不活动元素进行相关讨论。

    众所周知,基性岩浆岩源自地幔,而在Nb/YZr/TiO2×0.0001图解(图 3)中,辉长岩样品均落入玄武岩区域,表明辉长岩岩浆为玄武质岩浆。Ormerod等[41]指出,Zr/Ba值大于0.2的玄武岩源区来自软流圈地幔,Zr/Ba值小于0.2则来自岩石圈地幔,达肯大坂岩群辉长岩的Zr/Ba值为0.25~0.61,暗示其岩浆源区来自更深的软流圈地幔。此外,非常高的原始地幔标准化Th/Nb值(≫1)[42]和低Nb/La (< 1) [24]是地壳混染作用的2个可靠的微量元素指标。岩石地球化学分析数据显示,柴达木盆地西北缘达肯大坂岩群辉长岩样品具有高的ThN/NbN值(3.75~8.03)和低的Nb/La值(0.44~0.84)(图 9);另外,Ti元素在后期地质过程中不易受蚀变影响[43],且具有负Ti异常(相对于Eu),通常被认为是陆壳特征之一[44],因此负Ti异常也常被作为判断幔源岩浆是否遭受壳源混染的指标之一。辉长岩样品中TiO2含量为0.81%~0.95%,在微量元素原始地幔标准化蛛网图中呈现明显的负异常。上述主量、微量元素特征表明,辉长岩在形成过程中受到较强的壳源组分混染。辉长岩样品稀土元素总量较低,富集轻稀土元素,轻、重稀土元素内部分馏程度较弱,稀土元素配分模式表现为右倾型(图 5),也暗示辉长岩形成过程中遭受一定程度的壳源组分混染。

    图  9  辉长岩ThN/NbN-Nb/La图解(底图据参考文献[48])
    Figure  9.  ThN/NbN-Nb/La diagram for metagabbro

    锆石具有很高的Hf同位素体系封闭温度,可以记录岩浆源区不同性质的源岩特征,成为探讨岩浆起源、岩石成因甚至地壳演化及壳幔相互作用过程的重要工具[45]。研究表明,不同地球化学储库的176Lu/177Hf值不同,亏损地幔和球粒陨石的176Hf/177Hf值较高(≥0.282722),其对应的εHft)值为正值或零;地壳和富集地幔具有较低的176Hf/177Hf值,地壳的εHft)值通常为负值,而富集地幔的εHft)值可能为较小的正值或负值[45-46]。柴北缘达肯大坂岩群辉长岩的εHft)值为9.25~15.11,176Hf/177Hf值介于0.282815~0.283186之间,均大于0.282722,表明锆石母岩浆具有亏损地幔属性;然而,年龄-εHft)图解(图 8)显示,εHft)值均低于锆石形成时的亏损地幔εHft)值,暗示原始母岩浆来自于受富集组分影响的地幔[47]。锆石Hf模式年龄反映寄主岩石从亏损地幔中分离的时间,而锆石U-Pb年龄则代表锆石结晶年龄,如果锆石母岩浆直接来源于未受任何影响的亏损地幔,那么锆石结晶年龄应近似等于锆石Hf模式年龄。柴达木盆地西北缘达肯大坂岩群辉长岩锆石的Hf模式年龄(TDM)为378~610Ma,平均494Ma,明显大于其锆石结晶年龄357Ma,同样指示岩浆源区可能受到富集组分影响。综合上述岩石地球化学及锆石Hf同位素分析,认为柴达木盆地西北缘达肯大坂岩群辉长岩的岩浆源区为软流圈地幔并可能受到壳源富集组分的混染。

    研究表明,早古生代以来柴北缘经历了洋壳俯冲[2, 11-13] -大洋闭合及陆陆碰撞[49-50]碰撞后板块折返[12, 18]-后造山陆内伸展[21, 51-53]的完整的造山旋回。前已述及,沉积时代上限为晚志留世—早泥盆世的具典型伸展型磨拉石建造的牦牛山组被认为是柴北缘早古生代主造山作用结束的标志[21, 54]。李睿华等[53]通过柴北缘西段牛鼻子梁地区闪长岩的年代学、地球化学及Hf同位素研究,认为牛鼻子梁闪长岩的形成时代为361±3Ma,且形成于南祁连向柴达木板块俯冲碰撞后伸展的环境。吴才来等[12, 51]在柴北缘中段锡铁山—大柴旦—嗷唠山一带厘定出一系列晚泥盆世(约372Ma)具I型、S型成因的中酸性岩浆岩,并认为是受后造山伸展作用影响,造山带不同块体之间的均衡调整而产生的滑塌、拉伸,引起地壳的部分熔融形成的。Wang等[52]认为,柴北缘东段都兰地区侵位结晶年龄为386~360Ma的花岗质岩浆岩具后碰撞特征,其岩浆是造山带去根和软流圈上涌,在形成具有亏损地幔地球化学特征的镁铁质岩浆后,与中上地壳相互作用并引发中上地壳熔融的产物。上述分析说明,中—晚泥盆世(386~ 360Ma)整个柴北缘均处于后造山伸展的构造演化阶段。

    一般认为,基性岩脉的形成与岩石圈伸展和地壳拉张具有直接联系,代表了伸展拉张的构造环境[55]。本文辉长岩侵位于柴北缘构造带西端欧龙布鲁克地块西北缘达肯大坂岩群中,LA-ICP-MS锆石U-Pb定年结果表明,其形成时代为357±4Ma,即早石炭世。结合前文关于柴北缘构造环境的阐述及辉长岩代表的伸展环境,认为早石炭世柴北缘仍处于后造山伸展的构造演化阶段。此外,岩石地球化学及锆石Hf同位素相关分析表明,辉长岩岩浆源区为软流圈地幔,可能是早石炭世后造山伸展阶段柴北缘造山带去根、软流圈地幔上涌作用的产物。

    (1)欧龙布鲁克地块西北缘辉长岩中锆石206Pb/238U年龄加权平均值为357±4Ma,代表辉长岩的形成年龄。

    (2)早石炭世(约357Ma)柴北缘仍处于后造山伸展的构造演化阶段。

    (3)研究区辉长岩脉是后造山伸展阶段柴北缘造山带去根、软流圈地幔上涌作用形成的,其岩浆源区为软流圈地幔并可能受到壳源富集组分的混染。

    致谢: 野外地质调查及成文期间得到中国地质调查局沈阳地质调查中心付俊彧、杨晓平教授级高工的细心指导,在此表示衷心感谢。
  • 图  1   扎兰屯地区大地构造位置图及地质简图

    A—中亚造山带及邻区构造简图;B—大兴安岭及邻区构造单元划分简图;C—扎兰屯地区早古生代地质简图及岛弧结构图

    Figure  1.   Geotectonic location map and geological sketch map of Zhalantun area

  • Liu Y, Li W, Feng Z, et al. A review of the paleozoic tectonics in the eastern part of central asian orogenic belt[J]. Gondwana Research, 2017, 43: 123-148. doi: 10.1016/j.gr.2016.03.013

    杨现力. 扎兰屯浅变质岩系地质特征及碎屑锆石年代学研究[D]. 吉林大学博士学位论文, 2007.
    苗来成, 刘敦一, 张福勤, 等. 大兴安岭韩家园子和新林地区兴华渡口群和扎兰屯群锆石SHRIMP U-Pb年龄[J]. 科学通报, 2007, 52(5): 591-591. doi: 10.3321/j.issn:0023-074X.2007.05.016
    周建波, 王斌, 曾维顺, 等. 大兴安岭地区扎兰屯变质杂岩的碎屑锆石U-Pb年龄及其大地构造意义[J]. 岩石学报, 2014, 30(7): 1879-1888. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201407004.htm

    Qin T, Guo R R, Qian C, et al. Recognition of Early Paleozoic Magmatisms in the Supposed Proterozoic Basements of Zhalantun, Great Xing'an Range, NE China[J]. Acta Geologica Sinica, 2019, 93(3): 433-439. http://qikan.cqvip.com/Qikan/Article/Detail?id=7100185519

    Zhao C, Qin K Z, Song G X, et al. Early palaeozoic high-mg basalt-andesite suite in the duobaoshan porphyry cu deposit, NE China: constraints on petrogenesis, mineralization, and tectonic setting[J]. Gondwana Research, 2019, 71: 91-116. doi: 10.1016/j.gr.2019.01.015

图(1)
计量
  • 文章访问数:  2318
  • HTML全文浏览量:  472
  • PDF下载量:  1599
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-04-20
  • 修回日期:  2021-09-09
  • 网络出版日期:  2023-08-15
  • 刊出日期:  2021-12-14

目录

/

返回文章
返回