• 全国中文核心期刊
  • 中国科学院引文数据库核心期刊(CSCD)
  • 中国科技核心期刊
  • F5000优秀论文来源期刊
  • 荷兰《文摘与引文数据库》(Scopus)收录期刊
  • 美国《化学文摘》收录期刊
  • 俄罗斯《文摘杂志》收录期刊
高级检索

雄安新区地热水的化学场特征及影响因素分析

邢一飞, 王慧群, 李捷, 滕彦国, 张保健, 李燕燕, 王贵玲

邢一飞, 王慧群, 李捷, 滕彦国, 张保健, 李燕燕, 王贵玲. 雄安新区地热水的化学场特征及影响因素分析[J]. 中国地质, 2022, 49(6): 1711-1722. DOI: 10.12029/gc20220601
引用本文: 邢一飞, 王慧群, 李捷, 滕彦国, 张保健, 李燕燕, 王贵玲. 雄安新区地热水的化学场特征及影响因素分析[J]. 中国地质, 2022, 49(6): 1711-1722. DOI: 10.12029/gc20220601
XING Yifei, WANG Huiqun, LI Jie, TENG Yanguo, ZHANG Baojian, LI Yanyan, WANG Guiling. Chemical field of geothermal water in Xiong'an New Area and analysis of influencing factors[J]. GEOLOGY IN CHINA, 2022, 49(6): 1711-1722. DOI: 10.12029/gc20220601
Citation: XING Yifei, WANG Huiqun, LI Jie, TENG Yanguo, ZHANG Baojian, LI Yanyan, WANG Guiling. Chemical field of geothermal water in Xiong'an New Area and analysis of influencing factors[J]. GEOLOGY IN CHINA, 2022, 49(6): 1711-1722. DOI: 10.12029/gc20220601

雄安新区地热水的化学场特征及影响因素分析

基金项目: 

国家重点研发项目 2018YFC0604302

中国地质调查局项目 DD20189114

详细信息
    作者简介:

    邢一飞, 男, 1984年生, 博士, 高级工程师, 主要从事地热地质、水文地质等相关研究; E-mail: xingyifei@cags.ac.cn

    通讯作者:

    李捷, 女, 1988年生, 博士, 讲师, 主要从事水文地质相关研究工作; E-mail: lijie_lm@163.com

  • 中图分类号: P314

Chemical field of geothermal water in Xiong'an New Area and analysis of influencing factors

Funds: 

the projects of National Key Research and Development 2018YFC0604302

China Geological Survey DD20189114

More Information
    Author Bio:

    XING Yifei, male, born in 1984, Ph.D., senior engineer, mainly engages in geothermal geology and hydrogeology; E-mail: xingyifei@cags.ac.cn

    Corresponding author:

    LI Jie, female, born in 1988, Ph.D., lecturer, majors in hydrogeology; E-mail: lijie_lm@163.com

  • 摘要:
    研究目的 

    雄安新区地热水化学场特征及其影响因素对了解区域地热资源聚敛机制具有重要意义。

    研究方法 

    在区内5个主要构造单元(容城凸起、牛驼镇凸起、霸县凹陷、保定凹陷和高阳低凸起)开展地热水水化学采样和测试分析工作。

    研究结果 

    研究区内主要分布有砂岩热储和岩溶热储,水化学类型以Cl-Na型和Cl·HCO3-Na为主。从浅入深,地热水TDS值呈增大趋势。自西向东,地热水TDS含量表现为增大的趋势,指示着地下水流向。受断裂带导水作用影响,部分深部雾迷山组与馆陶组地热水发生混合现象。雾迷山组岩溶热储地热水的变质系数和脱硫系数低,热储封闭性较好,处于相对还原状态,而馆陶组和明化镇组砂岩热储地热水封闭性相对较差。

    结论 

    区内地下水化学场主要受水岩相互作用程度、断裂、储层封闭性等因素影响。研究成果对于认识雄安新区地热资源、科学合理开发利用地热能,推进北方地区冬季清洁取暖具有重要意义。

    创新点:采用雄安新区及其周边的近100组地下水化学样品,构建了雄安新区地热水化学场,探讨了其主要影响因素为水岩相互作用程度、储层封闭性等。

    Abstract:

    This paper is the result of geothermal geological survey engineering.

    Objective 

    The hydrochemical formation and evolution of the geothermal water in the Xiong'an New Area have great significant implications to understand the formation mechanism of geothermal resource.

    Methods 

    Geothermal water and cold groundwater in the Rongcheng uplift, Niutuozhen uplift, Baxian depression, Baoding depression, and Gaoyang low uplift were sampled for hydrochemical analyses.

    Results 

    There are two thermal reservoirs: Sandstone and karst thermal reservoirs. The main hydrochemical types were Cl·HCO3-Na and Cl-Na. The TDS values of geothermal water increased with well depths. Mixing of the geothermal water of the Wumishan Formation and the Guantao Formation was observed in the fault zone. The metamorphic coefficients and desulfurization coefficients of geothermal water in the Wumishan Formation were low, indicating a relatively closed thermal reservoir and reduced state, while the sealing of the reservoirs in the Guantao Formation and Minghuazhen Formation was relatively poor.

    Conclusions 

    The chemical field of geothermal water in the study area is controlled by the extent of water-rock interaction, faults, the sealing of the reservoirs, and so on. Our result is of great significance for understanding the geothermal resources in Xiong'an New Area and promoting clean heating in winter in northern China.

  • 西南“三江”地区构造位置属于特提斯构造域(图 1a),地质构造及成矿作用极其复杂地区,许多学者认为西南“三江”南段地区古特提斯构造演化阶段存在多个陆壳地块与洋盆相间排列的多岛洋构造格局(莫宣学等,1993刘本培等,2002俞赛赢等,2003王立全等, 2008, 2013;潘桂棠等,2013;任飞等,2017尹福光等,2017刘军平等,2017王保弟等,2018Liu et al., 20202021李生喜等,2023),而昌宁—孟连古特提斯洋的构造演化及其转换时限是目前三江地区构造演化研究中最热门的科学问题之一。

    图  1  西南三江地区大地构造位置图(a、b,据王保弟等,2018)、地质简图及采样点(c,据云南省地质调查院,2016
    1—第四系;2—古生代—中生代盖层;3—低级变质带;4—高级变质带;5—糜棱岩化花岗质岩石;6—花岗岩;7—二叠纪二长花岗岩;8—二叠纪花岗闪长岩;9—三叠纪二长花岗岩;10—白垩纪二长花岗岩;11—崇山岩群;12—断裂带;13—侵入界线;14—渐变过渡界线;15—地名;16—测年样品号及采样位置;17—研究区
    Figure  1.  Geotectonic location map(a, b, after Wang Baodi et al., 2018) and simplified geological map (c, modifed from Yunnan Institute of Geological Survey, 2016) in Sanjiang area
    1-Quaternary; 2-Paleozoic-Mesozoic caprock; 3-Low metamorphic zone; 4-High metamorphic zone; 5-Mylonitized granitic rock; 6-Granite; 7-Permian adamellite; 8-Permian granodiorite; 9-Triassic adamellite; 10-Cretaceous adamellite; 11-Chongshan group; 12-Fault zone; 13-Intrusion boundary; 14-Gradual transition boundary; 15-Toponym; 16-Dating sample number and location; 17-Research area

    目前昌宁—孟连古特提斯构造演化研究比较成熟,经历了早古生代中晚期-晚古生代特提斯洋俯冲消减、晚二叠世—早三叠世主碰撞汇聚、晚三叠世晚碰撞造山等阶段的演化过程(刘本培等,2002李文昌等,2010孔会磊等,2012;潘桂棠等,2013;任飞等,2017尹福光等,2017刘军平等,2017王保弟等,2018)。晚古生代特提斯洋俯冲消减过程沉积了以鲕粒-生物碎屑灰岩为主的石炭系—二叠系渔塘寨组(CPy),以洋盆深水硅质岩、硅质泥岩、放射虫硅质岩、大洋玄武岩为主的石炭系—二叠系光色组(CPg),及中酸性火山岩为主的石炭系—二叠系龙洞河组(CPl)总体构成了洋盆-洋岛-海山-弧后盆地古地理格局;而在其东侧的临沧—勐海地块,由于昌宁—孟连洋壳向东俯冲消减,在临沧—景洪一带出现中酸性岩浆活动,锆石U-Pb年龄主要集中在280~334 Ma,显示了汇聚型大陆边缘岩浆弧花岗岩的特点(魏君奇等,2008;Peng et al.,2008;孔会磊等,2012;徐桂香等,2015;孙载波等,2015刘军平等,2017王保弟等,2018),但昌宁—孟连古特提斯洋弧-陆俯冲向碰撞、汇聚转换的具体时限及岩石学证据目前尚未报道。

    笔者首次在风庆县小湾地区发现一套中二叠世晚期弱过铝质花岗岩岩体,目前并未有类似岩体的报道。本次运用锆石U-Pb LA-ICP-MS定年、岩石地球化学分析及岩石学对其进行详细研究,探讨了其岩石组合、侵位时代、岩浆演化、构造背景及源区,进而讨论了与昌宁—孟连古特提斯洋弧-陆俯冲向碰撞汇聚转换的关系,极大地丰富了昌宁—孟连古特提斯洋构造演化时限及岩浆岩记录。

    研究区位于西南“三江”南段,构造上位于羌塘—三江造山系(Ⅶ)和班公湖—双湖—怒江—昌宁—孟连对接带(Ⅷ),地层属于华南地层大区、羌北—昌都—思茅地层区、兰坪—思茅地层分区的澜沧地层小区及漾濞地层小区❶❷。研究区出露地层仅有古元古界崇山岩群,岩浆岩以二叠纪、三叠纪花岗岩为主,少量白垩纪花岗岩(图 1c);小湾花岗岩岩体与古元古界崇山岩群呈侵入接触关系,岩体中富含黑云母和暗色残留体包体,被后期三叠纪黑云二长花岗岩、白垩纪细粒钾长花岗岩侵入。1:25万凤庆幅把小湾花岗岩岩体时代置于中元古代,仅因其普遍糜棱岩化,发育片麻状构造等现象,但未有同素位年龄约束;其主要岩石类型以糜棱岩化片麻状似斑状中粒黑云二长花岗岩为主,少量糜棱岩化片麻状细粒花岗闪长岩及片麻状粗粒二长花岗质变晶初糜棱岩,岩石结构上由细粒向粗粒演化,变化的有序性和单向性明显,且在空间上紧密共生,形成时间、成分及结构变化上表现出清楚的亲缘和演化关系,说明它们是同一岩浆热事件的产物。野外三者为渐变过渡接触,整个岩体由中心至边部矿物颗粒由粗粒变为中细粒;本次以小湾花岗岩岩体为研究对象,采集了同位素、薄片、地球化学样品,野外特征见图 2a

    图  2  云南凤庆地区小湾花岗岩岩体野外照片和镜下照片特征
    a—野外露头;b—黑云二长花岗岩手标本;c—黑云二长花岗岩正交偏光照片;d—花岗闪长斑岩正交偏光照片;Kfs—钾长石;Qtz—石英;Pl—斜长石;Bt—黑云母
    Figure  2.  Photographs and granite features for the Xiaowan granite in Fengqing area, Yunnan
    a-Field outcrop; b-Biotite feldspar granite hand specimen; c-Biotite feldspar granite orthogonal polarization photo; d-Granodiorite orthogonal polarization photo; Kfs-K-feldspar; Qtz-Quartz; Pl-Plagioclase; Bt-Biotite

    糜棱岩化片麻状似斑状中粒黑云二长花岗斑岩:具变余似斑状花岗结构,中粒结构,片麻状构造。主要矿物成分由斜长石(27%)、钾长石(37%)、石英(31%)、黑云母(5%)组成,副矿物有磷灰石、锆石。矿物粒径多在2~5 mm。斜长石呈自形、半自形板状,内见聚片双晶、环带结构,已不同程度钠黝帘石化;钾长石呈他形、半自形板状,以具格子双晶微斜长石为主,部分为条纹长石,内见自形较好的斜长石包裹体及石英的穿孔交代结构,边缘见蠕英石包裹体,少部分粒度较粗大,构成斑晶;石英他形不规则粒状、集合状,充填于长石的间隙中。黑云母自形、半自形的片状,Ng-棕红色、Np-淡黄色,绿泥石化。手标本见图 2b及镜下特征见图 2c

    糜棱岩化片麻状细粒花岗闪长斑岩:具变余似斑状花岗结构,细粒结构,网结-网纹状定向构造。矿物成分主要由钾长石(15%~20%)、斜长石(50%~55%)、石英(20%~25%)、黑云母(3%~4%)、白云母(1%~2%)、长英质(1%~7%)组成,其中斑晶以斜长石为主,少量石英;基质主要为石英、钾长石及黑云母。副矿物主要有磷灰石、磁铁矿。镜下特征见图 2d

    片麻状粗粒二长花岗质变晶初糜棱岩:具变余初糜棱结构,粗粒结构,分异流动状、条痕状构造。矿物成分主要为斜长石(25%)、钾长石(39%)、石英(30%)、黑云母(1%)、电气石(< 5%);副矿物有磷灰石、锆石。矿物粒径多在0.1~2 mm。斜长石他形微细粒变晶状(部分不同程度地保留其半自形的板状特征,残余岩浆组构特征),部分粒度较粗大,保留原透镜状、眼球状碎斑外形,内见聚片双晶、环带结构,已不同程度钠黝帘石化。钾长石他形微细粒变晶状,以具格子双晶微斜长石为主,部分为条纹长石,部分粒度较粗大,保留原眼球状、透镜状碎斑外形,边缘普遍见蠕英石及斜长石包裹体。石英大部分呈多晶条带、矩形条带产出,重结晶。黑云母显微鳞片状,沿糜棱面理呈断续线痕状聚集产出,Ng-褐绿色、Np-淡黄色。电气石自形、半自形的柱状,沿长轴方向沿糜棱面理定向分布,为后期电气石化作用形成。

    锆石分选在南京宏创地质勘查技术服务有限公司完成,原岩样品经人工粉碎后,经人工淘洗后去除轻矿物部分,将得到的重砂部分经电磁选后得到含有少量杂质的锆石样品,最后在双目镜下挑选出晶型较好的锆石,制成锆石样品靶。对锆石进行反射光、透射光显微照相和阴极发光(CL)图像分析,最后根据反射光、透射光及锆石CL图像选择代表性的锆石颗粒和区域进行U-Pb测年。

    U-Pb同位素定年在湖北省地质实验室测试中心重点实验室利用LA-ICP-MS分析完成。测试仪器采用的是由美国Coherent Inc公司生产的GeoLasPro全自动版193 nm ArF准分子激光剥蚀系统(LA)和美国Agilent公司生产的7700X型电感耦合等离子质谱仪(ICP-MS)联用构成的激光剥蚀电感耦合等离子体质谱分析系统(LA-ICP-MS)。样品的同位素比值和元素含量采用ICPMSDataCal 9.0进行处理分析,加权平均年龄的计算及锆石年龄谐和图的绘制采用Isoplot3.0(Ludwing,2003)来完成。采用年龄为206Pb/238U年龄,其加权平均值的误差为2σ206Pb/238U(和207Pb/206Pb)平均年龄误差为95%置信度(王海然等,2013刘军平等,2018)。

    选择10件岩体样品进行主量元素、稀土和微量元素分析。样品磨碎至200目后,在武汉上谱分析科技有限公司岩石矿物研究室进行主量和微量元素分析测试。主量元素使用X-射线荧光光谱仪(XRF-1500)法测试。用0.6 g样品和6 g四硼酸锂制成的玻璃片在ShimadzuXRF-1500上测定氧化物的质量分数值,精度优于2%~3%。微量元素及稀土元素利用酸溶法制备样品,使用ICP-MS (ElementⅡ)测试,分析精度(按照GSR-1和GSR-2国家标准):当元素质量分数值大于10×10-6时,精度优于5%,当质量分数值小于10×10-6时,精度优于10%。

    本次用于锆石U-Pb LA-ICP-MS年龄测试的样品采集位置见图 1,样品分析数据见表 1。样品糜棱岩化似斑状中粒黑云二长花岗斑岩(D01-1)为浅灰色,具不等粒半自形—自形短柱状结构,片麻状构造。将样品先经手工粉碎,后按常规重力及电磁法浮选出锆石颗粒,最后在实体镜下挑选出晶型较好的锆石约150余粒。锆石多为浅紫红色,个别呈褐色,粒状、短柱状、碎粒状,金刚光泽,透明,部分具磨蚀特征,锆石长100~150 μm,少数达200 μm。

    表  1  云南凤庆地区小湾花岗岩岩体中锆石LA-ICP-MS U-Pb同位素分析结果
    Table  1.  LA-ICP-MS U-Pb data of zircons from the Xiaowan granite in Fengqing area, Yunnan
    下载: 导出CSV 
    | 显示表格

    根据样品阴极发光(CL)图像锆石自形程度较好,棱角分明,发育较宽的韵律环带,Th/U比值较大,为典型的酸性岩浆结晶锆石(Hoskin and Schaltegger, 2003吴元保等,2004刘军平等,2018);部分锆石外围发育较亮较窄的变质增生边,可能是后期发生变质作用所致(图 3)。

    图  3  云南凤庆小湾地区花岗岩岩体代表性锆石阴极发光(CL)图像
    Figure  3.  The CL images of zircons from the Xiaowan granite in Fengqing area, Yunnan

    本次对糜棱岩化片麻状似斑状中粒黑云二长花岗岩(D01-1)共完成18个测点的锆石年龄测定(表 1),获得有效数据18个,18个有效分析点均具有较好的谐和度,谐和度均大于91%,从图 4a可以看出,分析的18个点均沿谐和线分布且较集中,且振荡环带清楚,Th/U比值较大(值为0.4~0.8),均为岩浆成因锆石,18个分析点获得206Pb/238U加权平均年龄为(260.4±1.30)Ma(MSWD=0.32,n=18),代表黑云二长花岗岩侵位时代,说明其形成于中二叠世晚期,而非前人认为的中元古代。

    图  4  云南凤庆地区小湾花岗岩岩体锆石U-Pb年龄谐和图
    Figure  4.  Zircon U-Pb concordia diagrams of the Xiaowan granite in Fengqing area, Yunnan

    小湾花岗岩的主量元素分析数据列于表 2。岩石化学成分SiO2含量为67.54%~78.66%,Al2O3含量为11.34%~14.30%,K2O+Na2O含量为6.20%~8.52%,里特曼指数σ为1.31~2.51,显示岩石属过铝质钙碱性岩类;在深成岩的SiO2-(Na2O+K2O)图解中(图 5a)样品投影点多数样品落入花岗岩区内,一件样品投影点落入花岗闪长岩区,与镜下矿物命名一致,均为亚碱性岩石;岩体A/CNK-A/NK图解中大部分样品落入准铝质区-过铝质区过渡区域,表现出弱过铝质的特点,仅1件样品落入过铝质区(图 5b)。在SiO2-K2O图解中大部分落入高钾钙碱性系列区(图 6),1件样品落入钙碱性系列。样品CIPW标准矿物显示,岩石为高硅富铝的特点,与岩石地球化学含量较为一致。

    表  2  云南凤庆地区小湾花岗岩岩体的主量元素(%)和微量元素(10-6)组成及有关参数
    Table  2.  Analytical results of major elements (%) and trace elements (10-6) of the Xiaowan granite in Fengqing area, Yunnan
    下载: 导出CSV 
    | 显示表格
    图  5  云南凤庆小湾地区花岗岩岩体SiO2-(Na2O+K2O)(a) and (A/CNK)-(A/NK)(b)图解(a,据Middlemost et al., 1994;b,据Maniar and Piccoli, 1989
    Figure  5.  SiO2-(Na2O+K2O)(a) and (A/CNK)-(A/NK) (b) diagrams of the Xiaowan granite in Fengqing area, Yunnan (a, after Middlemost et al., 1994; b, after Maniar and Piccoli, 1989)
    图  6  云南凤庆地区小湾花岗岩岩体SiO2-K2O图解(据Richwood, 1989
    Figure  6.  SiO2-K2O diagrams of the Xiaowan granite in Fengqing area, Yunnan (after Richwood, 1989)

    稀土元素、微量元素分析成果及特征参数见表 1。小湾花岗岩稀土总量(ΣREE)为125.64×10-6~300.04×10-6,平均为198.54×10-6,含量较低;轻稀土元素(LREE)含量为83.70×10-6~243.04×10-6,重稀土元素(HREE)含量为31.23×10-6~55.54×10-6,LREE/HREE为2.04~4.93,显示轻稀土元素相对于重稀土元素有一定程度的富集;(La/Yb)N=3.86~15.84,平均10.99;(La/Sm)N=2.66~7.08,平均5.37;(Gd/Yb)N=1.02~1.98,平均1.42;上述特征显示轻、重稀土之间具有明显分异特征,轻稀土元素内部的分异程度较高,稀土元素分异程度高,反映源区地壳具有较高的成熟度。小湾花岗岩δEu介于0.27~0.81,平均为0.53,在球粒陨石标准化图解上(图 7a)配分曲线Eu负异常明显,具有向右缓倾的V型特征,暗示岩石可能经历了强烈的斜长石分离结晶作用,长石矿物的分离结晶作用也可以解释花岗岩中Ba、Sr的亏损。且普遍具有弱程度的负铈异常(δCe=0.89~1.32,平均为0.96),暗示了岩浆形成于较高的氧逸度环境(刘洪等,2016刘军平等, 2017, 2020a, 2020b)。

    图  7  云南凤庆地区小湾花岗岩岩体稀土元素配分形式(a)与微量元素蛛网图(b)(原始地幔标准化数据、球粒陨石标准化数据引自文献Sun and McDonough, 1989
    Figure  7.  Chondrite-normalized rare earth elements patterns (a) and primitive mantle-normalized trace elements patterns (b) of the Xiaowan granite in Fengqing area, Yunnan (Chondrite and primitive mantle normalized data after Sun and McDonough, 1989)

    小湾花岗岩原始地幔标准化微量元素蛛网图解(图 7b)显示,样品具有明显的负Eu异常,相对富集Rb、Th、Ce等大离子亲石元素(LILE),相对亏损Ta、Nb、Zr、Ti、Hf、Y等高场强元素(HFSE)的特征;样品曲线形态趋势相近,它们应该具有相似的源区(刘洪等,2016张洪亮等,2019Liu et al., 2023)。

    研究区小湾花岗岩A/CNK值均在1.0以上,最高达1.47,表现出弱过铝质特点。CIPW标准矿物中出现刚玉分子(平均值为0.94%),Eu负异常明显,δEu平均为0.53,且具有低FeOT/MgO、低P2O5(0.03%~0.16%)和高K2O+Na2O(6.20%~8.52%)等地球化学特征,与弧火山岩及同碰撞花岗岩有部分类似特点(吴福元等,2007),可能与昌宁—孟连洋壳俯冲消亡向弧-陆碰撞汇聚对接转换有关,为过渡类型。

    不同性质的源岩形成了不同构造环境中的弱过铝质花岗岩;判别弱过铝质花岗岩的形成环境则需要判别其形成的源岩性质,对于判别弱过铝质花岗岩的源岩性质,斜长石是一个很好的指示矿物,泥质源岩贫斜长石(< 5%),而杂砂岩则富斜长石(> 5%)。一般CaO/Na2O < 0.3表明原岩为泥质岩石的部分熔融,同时岩石Rb/Sr比值> 0.1和Rb/Ba比值> 0.3,反之若CaO/Na2O < 0.3表明原岩为杂砂岩的部分熔融,其Rb/Sr和Rb/Ba比值较低(Barbarin,1999邓晋福等,2004)。小湾花岗岩的CaO/Na2O比值大部分 < 0.3(3个样品大于0.3),Rb/Sr比值均 > 0.1和大部分Rb/Ba比值> 0.3为特征,结合CaO/(MgO+FeOT)-Al2O3/(MgO+FeOT)图解(图 8a),显示其应为泥质岩石夹部分杂砂岩局部熔融而形成。

    图  8  云南凤庆地区小湾花岗岩岩体CaO/(MgO+FeOT)-Al2O3/(MgO+FeOT) (a)及Yb-Sr(b)构造判别图(a,据Alter et al., 2000;b,据张旗等,2006)
    Figure  8.  (CaO/(MgO+FeOT)-Al2O3/(MgO+FeOT)) (a)-(Yb-Sr) (b) tectonic discrimination diagrams of the Xiaowan granite in Fengqing area, Yunnan (a, after Alter et al., 2000; b, after Zhang Qi et al., 2006)

    研究区的小湾花岗岩在张旗的Yb-Sr图解(图b)中样品落入极低Sr高Yb区和低Sr高Yb(图 9a),样品表现为低Sr高Yb型向高Sr低Yb型演化趋势,指示出地壳厚度增厚压力增大,可能与昌宁—孟连洋俯冲、碰撞造山有关。小湾花岗岩岩体6件样品(D01-1~D01-6),DI指数为93.94~90.26,分异程度较强,另一组4件样品(D01-7~D01-10),DI指数为71.66~88.65,分异程度弱或基本未分异,可能与本区花岗岩具弧花岗岩及同碰撞花岗岩双重过渡性质有关,应处于两者间的过渡类型。

    图  9  云南凤庆地区小湾花岗岩岩体(Y+Nb)-Rb(a)及Yb-Ta(b)图解(据Pearce et al., 1984)
    Figure  9.  (Y+Nb)-Rb(a) and Yb-Ta (b) diagram of the Xiaowan granite in Fengqing area, Yunnan(after Pearce et al., 1984)

    花岗质岩石可以形成于多种构造环境,比如活动大陆边缘、岛弧造山带、大陆碰撞带等构造部位。

    小湾花岗岩岩体富集轻稀土及大离子亲石元素,亏损高场强元素,且具弱过铝质特点,在稀土元素和微量元素图解方面显示出岛弧花岗岩及同碰撞花岗岩双重特征,在(Y+Nb)-Rb图解(图 9a)及Yb-Ta图解(图 9b)中样品落入火山弧花岗岩及同碰撞花岗岩交汇处,应具两者双重性质,为两者过渡类型。微量元素Sr的负异常和Th的正异常也显示了岛弧岩浆岩的特征,暗示着岩浆源区可能受到过俯冲板块流体或熔体的交代作用。中晚二叠世时期为昌宁—孟连洋向东俯冲未期,由于俯冲板块体或熔体的交代作用导致楔形地幔源区发生部分熔融作用,幔源岩浆上涌,从而导致上地壳物质升温发生熔融,研究区中二叠世晚期具有典型岛弧地球化学特征的小湾花岗岩即为该事件的岩石学记录,相对于大规模俯冲消减作用的弧岩浆岩记录(280~305 Ma),本期花岗岩滞后了约20 Ma,此时区域上开始转入弧-陆碰撞阶段,由于本期花岗岩侵位时代处于二次构造事件转换阶段,故同时兼具弧花岗岩及同碰撞花岗岩双重地球化学特征。结合区域资料,昌宁—孟连带在早二叠世末—早三叠世全区进入弧-陆俯冲向碰撞汇聚转换阶段(280~245 Ma);早三叠世—晚三叠世全区进入主碰撞阶段(245~210 Ma);之后全区进入伸展拉张阶段(~210 Ma),主要沉积了以早侏罗统芒汇河组,小红桥组为代表的火山岩夹碎屑岩沉积(刘军平, 2017, 2020c王保弟等,2018吕留彦等,2019)。由些,笔者认为小湾花岗岩可能是昌宁—孟连洋壳俯冲消亡向弧-陆碰撞、汇聚对接作用的产物(Collins and Richards, 2008王保弟等,2018),而非碰撞后花岗岩。

    石炭纪—晚二叠世是昌宁—孟连古特提斯洋盆扩张、俯冲消减及弧-陆碰撞、汇聚转换的重要时期;其岩浆作用在凤庆、澜沧、景洪等地区均有该期侵入岩及火山岩的报道,该期典型代表有曼秀岩体、早石炭世平掌组(C1pz)火山岩,显示岛弧岩浆岩的地球化学特征,可能代表了昌宁—孟连洋壳发生(洋内)俯冲消减开始的记录,其岩浆记录主要集中在334 Ma以及320~348 Ma锆石U-Pb年龄(孔会磊等, 2012, 孙载波等,2015王保弟等,2018);具有典型的俯冲形成的岛弧型岩浆岩年龄主要集中在280~305 Ma,代表性岩体主要有南联山岩体锆石年龄为300~305 Ma(徐桂香等,2016);雅口岩体堆晶辉长岩的锆石年龄为296 Ma;半坡岩体锆石年龄为284~288 Ma(Wang et al., 2010);吉岔岩体锆石年龄280~297 Ma(魏君奇等,2008)。

    本次发现的小湾花岗岩锆石U-Pb LA-ICP-MS年龄为260.4 Ma,为弱过铝质花岗岩,具弧火山岩及同碰撞花岗岩双重特征,应为过渡类型;260.4 Ma是目前临沧岩浆弧发现的最为年轻的具弧岩浆岩及同碰撞花岗岩双重性质的酸性岩浆记录,与云县漫湾地区龙洞河组(CPl)中二叠世弧火山岩年龄(262.5±3.1 Ma)刘军平等,2022)及景洪大勐龙怕冷具岛弧性质的辉长-辉绿岩岩体258~262 Ma的年龄相当王保弟等,2018)。因些,结合前人的研究成果,昌宁—孟连洋构造体制由弧-陆碰撞到汇聚对接作用的转换可能发生在280~258 Ma,并引发一定规模的构造热事件,本文所报道的小湾花岗岩体应是该期事件的产物,其转换时限可能为260.4 Ma。笔者推测在云南境内昌宁—孟连洋盆向东俯冲消减作用最早碰撞发生在凤庆小湾、云县漫湾及景洪大勐龙一带,后逐渐相向延伸。

    (1)小湾花岗岩岩体结晶年龄为(260.4±1.3)Ma(MSWD=0.32,n=18),表明岩体形成于中二叠世晚期,而非前人认为的中元古代,是目前临沧岩浆弧发现的最为年轻的具弧花岗岩及同碰撞花岗岩双重性质的酸性岩浆记录。

    (2)小湾花岗岩岩体是在昌宁—孟连洋壳向东俯冲向汇聚转换背景下,上地壳泥质岩石夹部分杂砂岩成分发生部分熔融作用形成的弱过铝质亚碱性花岗岩。

    (3)小湾花岗岩岩体应是昌宁—孟连洋壳俯冲消亡向弧-陆碰撞、汇聚对接作用的产物,其转换时限可能为260.4 Ma;结合区域资料,云南境内昌宁—孟连洋盆向东俯冲消减作用最早碰撞可能发生在凤庆小湾、云县漫湾及景洪大勐龙一带,后逐渐相向延伸。

    注释

    ❶云南省地质调查院. 2008.云南1∶25万凤庆幅区域地质矿产调查报告[R].

    ❷云南省地质调查院. 2016.云南1∶5万诗礼幅、大河幅、犀牛街幅、老家库幅、哨街幅区域地质矿产调查报告[R].

    ❸云南省地质调查院. 2012.云南1∶25万澜沧幅、景洪幅、勐腊幅、勐海幅区域地质调查报告[R].

  • 图  1   采样点位置分布图

    1—第四系冷水;2—明化镇组地热水;3—馆陶组地热水;4—雾迷山组地热水

    Figure  1.   Distribution map of sampling points

    1-Quaternary cold water; 2- Geothermal water of Minghuazhen Formation (Nm); 3- Geothermal water of Guantao Formation (Ng); 4-Geothermal water of Wumishan Formation (Jxw)

    图  2   地热水(a~e)和第四系地下水(f)Piper三线图

    1—雾迷山组;2—馆陶组;3—明化镇组;4—燕山;5—容城凸起;6—高阳低凸起;7—牛驼镇凸起;8—霸县凹陷

    Figure  2.   Ternary diagrams of water samplings from geothermal water (a-e) and groundwater(f)

    1-Wumishan Formation (Jxw); 2-Guantao Formation (Ng); 3-Minghuazhen Formation (Nm); 4 -Yanshan Mountains; 5-Rongcheng uplift; 6— Gaoyang low uplift; 7-Niutuozhen uplift; 8-Baxian depression

    图  3   研究区南部(a)和北部(b)不同地层地热水中TDS分布

    Figure  3.   Distribution characteristics of TDS in geothermal water of different strata in the southern area (a) and northern area (b)

    图  4   地热水离子浓度

    1—容城凸起;2—牛驼镇凸起;3—霸县凹陷;4—保定凹陷;5—高阳低凸起

    Figure  4.   Ion concentration of geothermal water

    1-Rongcheng uplift; 2-Niutuozhen uplift; 3-Baxian depression; 4-Baoding depression; 5-Gaoyang low uplift

    图  5   研究区南部(a)和北部(b)SO42-与井深的关系

    Figure  5.   Relationship between SO42- and well depths in the southern area (a) and northern area (b)

    图  6   地热水γNa+/γCl-(a)、100×γSO42-/γCl-(b)与Cl-的关系

    Figure  6.   Relationship between γNa+/γCl-(a), 100×γSO42-/γCl-(b) and Cl-

    图  7   雄安新区地热水水化学成因模式示意图(修改自Cui et al., 2019)

    Figure  7.   Conceptual model of hydrochemical mechanism of Xiong'an geothermal system (modified from Cui et al., 2019)

    表  1   各单元水化学特征(平均值)

    Table  1   Hydrochemical characteristics of each unit (average value)

    下载: 导出CSV
  • Chen Moxiang. 1988. Geothermal Energy in North China[M]. Beijing: Science Press.

    Cui Yue, Zhu Chuanqing, Qiu Nansheng, Tang Boning, Guo Sasa. 2019. Radioactive heat production and terrestrial heat flow in the Xiong'an Area, North China[J]. Energies, 12: 4608. doi: 10.3390/en12244608

    Cui Yue. 2020. Genetic Mechanism of Geothermal Resource Heat Source in Wumishan Formation of Xiong'an New Area[D]. Beijing: China University of Petroleum (Beijing).

    Feng Xin, Zhang Yazhe. 2014. Analysis and research on ion proportion coefficient of groundwater in Shenzhou area[J]. China Rural Water and Hydropower, (4): 18-20, 24(in Chinese with English abstract). doi: 10.3969/j.issn.1007-2284.2014.04.006

    Giggenbach W F. 1988. Geothermal solute equilibria. Derivation of Na-K-Mg-Ca geoindicators[J]. Geochimica et Cosmochimica Acta, 52(12): 2749-2765. doi: 10.1016/0016-7037(88)90143-3

    Guo Qi, Pang Zhonghe, Wang Yingchun, Tian Jiao. 2017. Fluid geochemistry and geothermometry applications of the Kangding high-temperature geothermal system in eastern Himalayas[J]. Applied Geochemistry, 81: 63-75. doi: 10.1016/j.apgeochem.2017.03.007

    Guo Qinghai, Wang Yanxin, Liu Wei. 2007. Major hydrogeochernical processes in the two reservoirs of the Yangbajing geothermal field, Tibet, China[J]. Journal of Volcanology and Geothermal Research, 166(3/4): 255-268.

    Guo Sasa, Zhu Chuanqing, Qiu Nansheng, Tang Boning, Cui Yue, Zhang Jiatang, Zhao Yuhang. 2019. Present geothermal characteristics and influencing factors in the Xiong'an New Area, North China[J]. Energies, 12(20): 1-22.

    Guo Sasa, Zhu Chuanqing, Qiu Nansheng, Tang Boning, Cui Yue. 2020. Formation conditions and favorable areas for the deep geothermal resources in the Xiong'an New Area[J]. Acta Geologica Sinica, 94(7): 2026-2035(in Chinese with English abstract). doi: 10.3969/j.issn.0001-5717.2020.07.011

    Han Zheng, Cui Yiqiao, Wang Shufang, Li Xiao, Sun Ying. 2015. Geothermal resource assessment based on Monte-Carlo method——A case study of geothermal field in Xiong County of Hebei Province[J]. Urban Geology, 4: 58-62(in Chinese with English abstract).

    Jiang Haiyang, Wang Shuxing, Liu Lian, Cao Yanling, Shi Meng, Li Tingting, Sun Wenguang. 2018. Hydrochemical and environmental isotope characteristics of geothermal resources in Linyi section of Yishu fault zone[J]. Land and Resources of Shanghai, 39(4): 90-94(in Chinese with English abstract).

    Kong Yanlong, Pang Zhonghe, Pang Jumei, Li Jie, Pan Sheng. 2020. Fault-affected fluid circulation revealed by hydrochemistry and isotopes in a large-scale utilized geothermal reservoir[J]. Geofulids, 2020(24): 1-13.

    Li Hu, Xia Lingxiao, Jiang Guosheng. 2015. Analysis and research on ion proportion coefficient of Ordovician geothermal fluid in Tianjin[J]. Groundwater, 37(5): 35-37(in Chinese with English abstract). doi: 10.3969/j.issn.1004-1184.2015.05.011

    Li Kewen, Bian Huiyuan, Liu Changwei, Zhang Danfeng, Yang Yanan. 2015. Comparison of geothermal with solar and wind power generation systems[J]. Renewable and Sustainable Energy Reviews, 42: 1464-1474. doi: 10.1016/j.rser.2014.10.049

    Li Yanyan, Zhang Baojian, Xing Yifei, Wang Guiling. 2021. Study on the fracture law of carbonate rocks of Gaoyuzhuang Formation under different confining pressures in Gaoyang geothermal field of Xiong'an New Area[J]. Geology in China: 1-15[2021-04-19] https://kns.cnki.net/kcms/detail/11.1167.P.20210419.1328.006.html (in Chinese with English abstract).

    Liu Mingliang, He Tong, Wu Qifan, Guo Qinghai. 2020. Hydrogeochemistry of geothermal waters from Xiong'an New Area and its indicating significance[J]. Earth Science, 45(6): 2221-2231.

    Ma Feng, Wang Guiling, Zhang Wei, Zhu Xi, Zhang Hanxiong, Yue Gaofan. 2020. Thermal storage space structure and resource potential of Rongcheng geothermal field in Xiong'an New Area[J]. Acta Geologica Sinica, 94(7): 1981-1990(in Chinese with English abstract). doi: 10.3969/j.issn.0001-5717.2020.07.007

    Pang Jumei, Pang Zhonghe, Lü Min, Tian Jiao, Kong Yanlong. 2018. Geochemical and isotopic characteristics of fluids in the Niutuozhen geothermal field, North China[J]. Environmental Earth Sciences, 77(1): 12. doi: 10.1007/s12665-017-7171-y

    Pang Zhonghe, Kong Yanlong, Pang Jumei, Hu Shengbiao, Wang Jiyang. 2017. Research on geothermal resources and development and utilization in Xiong'an new area[J]. Journal of Chinese Academy of Sciences, 32(11): 1224-1230(in Chinese with English abstract).

    Wang Guiling, Gao Jun, Zhang Baojian, Xing Yifei, Zhang Wei, Ma Feng. 2020. Research on thermal reservoir characteristics and high productivity geothermal well parameters of Wumishan Formation in Gaoyang low uplift area of Xiong'an New Area[J]. Acta Geologica Sinca, 94(7): 1970-1980(in Chinese with English abstract).

    Wang Guiling, Liu Yanguang, Zhu Xi, Zhang Wei. 2020. Current situation and development trend of geothermal resources in China[J]. Earth Science Frontiers, 27(1): 1-9(in Chinese with English abstract).

    Wang Guiling, Wang Wanli, Zhang Wei, Ma Feng, Liu Feng. 2020. The status quo and prospect of geothermal resources exploration and development in Beijing-Tianjin-Hebei region in China[J]. China Geology, 3(1): 173-181. doi: 10.31035/cg2020013

    Wang Guiling, Zhang Wei, Lin Wenjing, Liu Feng, Zhu Xi, Liu Yanguang, Li Jun. 2017. Study on the accumulation model and potential of geothermal resources in Beijing Tianjin Hebei region[J]. Geology in China, 44(6): 1074-1085(in Chinese with English abstract).

    Wang Jiyang, Gong Yulie, Lu Zhenneng, Ma Weibin. 2013. Looking at China's geothermal development and utilization from the perspective of European geothermal development[J]. Advances in New and Renewable Energy, 1(1): 1-6(in Chinese with English abstract).

    Wang Jiyang, Pang Zhonghe. 2021. Geothermal energy welcomes unprecedented development opportunities[N]. China Science Daily, 2021-01-06(1).

    Wang Jiyang, Qiu Nansheng, Hu Shengbiao, He Lijuan. 2017. Progress and development trend of geothermal research in China's oil fields[J]. Earth Science Frontiers, 24(3): 1-12(in Chinese with English abstract).

    Wang Jiyang. 2015. Geothermy and its Application[M]. Beijing: Science Press: 257-376(in Chinese).

    Wang Jiyang. 2020. Clean geothermal heating is promising[N]. China Science Daily, 2020-09-02(2).

    Wang Zhuting, Zhang Chao, Jiang Guangzheng, Hu Jie, Tang Xianchun, Hu Shengbiao. 2019. Characteristics and genetic mechanism of current geothermal field in Xiong'an New Area[J]. Chinese Journal of Geophysics, 62(11): 4313-4322(in Chinese with English abstract). doi: 10.6038/cjg2019M0326

    Xia Liangyu, Zhang Yabo. 2019. An overview of world geothermal power generation and a case study on China-The resource and market perspective[J]. Renewable and Sustainable Energy Reviews, 112: 411-423. doi: 10.1016/j.rser.2019.05.058

    Yan Xiaoxue, Gan Hao, male, Yue Gaofan. 2019. Hydrogeochemical characteristics and genesis analysis of typical geothermal fields in Huizhou Conghua, Guangdong[J]. Geological Review, 65(3): 743-754(in Chinese with English abstract).

    Zhang Meng, Lin Wenjing, Liu Zhao, Liu Zhiming, Hu Xiancai, Wang Guiling. 2014. Hydrogeochemical characteristics and genetic model of Gulu high-temperature geothermal system in Tibet[J]. Journal of Chengdu University of Technology (Science &Technology Edition), 41(3): 382-392(in Chinese with English abstract). doi: 10.3969/j.issn.1671-9727.2014.03.15

    Zhao Xingang, Wan Guan. 2014. Current situation and prospect of China's geothermal resources[J]. Renewable and Sustainable Energy Reviews, 32: 651-661. doi: 10.1016/j.rser.2014.01.057

    Zhou Cheng, Elham Doroodchi, Behdad Moghtaderi. 2013. An in-depth assessment of hybrid solar-geothermal power generation[J]. Energy Conversion and Management, 74: 88-101. doi: 10.1016/j.enconman.2013.05.014

    Zhou Zongying, Liu Shiliang, Liu Jinxia. 2015. Characteristics and development countermeasures of China's geothermal resources[J]. Journal of Natural Resources, 30(7): 1210-1221(in Chinese with English abstract).

    陈墨香. 1988. 华北地热[M]. 北京: 科学出版社.
    崔悦. 2020. 雄安新区雾迷山组地热资源热源成因机制[D]. 北京: 中国石油大学(北京).
    冯欣, 张亚哲. 2014. 深州地区地下水离子比例系数分析研究[J]. 中国农村水利水电, (4): 18-20, 24. doi: 10.3969/j.issn.1007-2284.2014.04.006
    郭飒飒, 朱传庆, 邱楠生, 唐博宁, 崔悦. 2020. 雄安新区深部地热资源形成条件与有利区预测[J]. 地质学报, 94(7): 2026-2035. doi: 10.3969/j.issn.0001-5717.2020.07.011
    韩征, 崔一娇, 王树芳, 李潇, 孙颖2015. 基于蒙特卡洛法的地热资源评价-以河北省雄县地热田为例[J]. 城市地质, 4: 58-62 doi: 10.3969/j.issn.1007-1903.2015.04.011
    江海洋, 王树星, 刘连, 曹艳玲, 史猛, 李婷婷, 孙文广. 2018. 沂沭断裂带临沂段地热资源水化学及环境同位素特征[J]. 上海国土资源, 39(4): 90-94. doi: 10.3969/j.issn.2095-1329.2018.04.019
    李虎, 夏玲晓, 江国胜. 2015. 天津地区奥陶系地热流体离子比例系数分析研究[J]. 地下水, 37(5): 35-37. doi: 10.3969/j.issn.1004-1184.2015.05.011
    李燕燕, 张保建, 邢一飞, 王贵玲. 2021. 雄安新区高阳地热田热储高于庄组碳酸盐岩不同围压下破碎规律研究[J]. 中国地质: 1-15[2021-04-19] https://kns.cnki.net/kcms/detail/11.1167.P.20210419.1328.006.html.
    马峰, 王贵玲, 张薇, 朱喜, 张汉雄, 岳高凡. 2020. 雄安新区容城地热田热储空间结构及资源潜力[J]. 地质学报, 94(7): 1981-1990. doi: 10.3969/j.issn.0001-5717.2020.07.007
    庞忠和, 孔彦龙, 庞菊梅, 胡圣标, 汪集暘. 2017. 雄安新区地热资源与开发利用研究[J]. 中国科学院院刊, 32(11): 1224-1230.
    汪集暘, 龚宇烈, 陆振能, 马伟斌. 2013. 从欧洲地热发展看我国地热开发利用问题[J]. 新能源进展, 1(1): 1-6. https://www.cnki.com.cn/Article/CJFDTOTAL-XNYJ201301002.htm
    汪集暘, 庞忠和. 2021. 地热能迎空前发展机遇[N]. 中国科学报, 2021-01-06(1).
    汪集暘, 邱楠生, 胡圣标, 何丽娟. 2017. 中国油田地热研究的进展和发展趋势[J]. 地学前缘, 24(3): 1-12. https://www.cnki.com.cn/Article/CJFDTOTAL-DXQY201703002.htm
    汪集暘. 2015. 地热学及其应用[M]. 北京: 科学出版社: 257-376.
    汪集暘. 2020. 地热清洁取暖大有可为[N]. 中国科学报, 2020-09-02(2).
    王贵玲, 高俊, 张保建, 邢一飞, 张薇, 马峰. 2020. 雄安新区高阳低凸起区雾迷山组热储特征与高产能地热井参数研究[J]. 地质学报, 94(7): 1970-1980. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXE202007006.htm
    王贵玲, 刘彦广, 朱喜, 张薇. 2020. 中国地热资源现状及发展趋势[J]. 地学前缘, 27(1): 1-9. https://www.cnki.com.cn/Article/CJFDTOTAL-DXQY202001002.htm
    王贵玲, 张薇, 蔺文静, 刘峰, 朱喜, 刘彦广, 李郡. 2017. 京津冀地区地热资源成藏模式与潜力研究[J]. 中国地质, 44(6): 1074-1085. http://geochina.cgs.gov.cn/geochina/article/abstract/20170603?st=search
    王朱亭, 张超, 姜光政, 胡杰, 唐显春, 胡圣标. 2019. 雄安新区现今地温场特征及成因机制[J]. 地球物理学报, 62(11): 4313-4322. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201911026.htm
    闫晓雪, 甘浩男, 岳高凡. 2019. 广东惠州-从化典型地热田水文地球化学特征及成因分析[J]. 地质论评, 65(3): 743-754. https://www.cnki.com.cn/Article/CJFDTOTAL-DZLP201903020.htm
    张萌, 蔺文静, 刘昭, 刘志明, 胡先才, 王贵玲. 2014. 西藏谷露高温地热系统水文地球化学特征及成因模式[J]. 成都理工大学学报(自然科学版), 41(3): 382-392. https://www.cnki.com.cn/Article/CJFDTOTAL-CDLG201403015.htm
    周总瑛, 刘世良, 刘金侠. 2015. 中国地热资源特点与发展对策[J]. 自然资源学报, 30(7): 1210-1221. https://www.cnki.com.cn/Article/CJFDTOTAL-ZRZX201507013.htm
  • 期刊类型引用(6)

    1. 王卓群,陈扣平,吴吉春,王保战,柳旭. 太湖周边地下水中高硝酸盐负荷促进反硝化功能菌群富集. 中国环境科学. 2025(01): 477-486 . 百度学术
    2. 李丽君,李旭光. 西辽河平原浅层地下水中“三氮”分布特征及健康风险评价. 地质与资源. 2024(01): 90-97 . 百度学术
    3. 杨开丽,刘淑亮,万惠芸,阎琨,陶昱君,王洪松. 基于氮氧同位素和MixSIAR模型的烟台北部地表水硝酸盐来源分析. 中国地质. 2024(06): 2066-2076 . 本站查看
    4. 蒋宗妤,马圣祥,刘莹,马兴冠. 复杂地貌地区地下水硝酸盐生态及健康风险评价. 给水排水. 2024(S1): 7-11 . 百度学术
    5. 许飞青,李潇,李凯,于喆,郭亚杉,赵志强. 随机森林回归模型在地下水水质评价的新应用. 地质与勘探. 2023(02): 408-417 . 百度学术
    6. 王国华,张妍,缑倩倩,张仲伍,孙九林. 黑河流域中游绿洲边缘地表水和地下水水化学特征分析. 地理科学. 2022(10): 1818-1828 . 百度学术

    其他类型引用(5)

图(7)  /  表(1)
计量
  • 文章访问数:  2157
  • HTML全文浏览量:  984
  • PDF下载量:  2419
  • 被引次数: 11
出版历程
  • 收稿日期:  2021-05-28
  • 修回日期:  2021-09-22
  • 网络出版日期:  2023-09-25
  • 刊出日期:  2022-12-24

目录

/

返回文章
返回
x 关闭 永久关闭