• 中国出版政府奖提名奖

    中国百强科技报刊

    湖北出版政府奖

    中国高校百佳科技期刊

    中国最美期刊

    留言板

    尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

    姓名
    邮箱
    手机号码
    标题
    留言内容
    验证码

    以SiO2含量为辅助因子的ASTER热红外遥感硅化信息提取

    刘道飞 陈圣波 陈磊 马明

    刘道飞, 陈圣波, 陈磊, 马明, 2015. 以SiO2含量为辅助因子的ASTER热红外遥感硅化信息提取. 地球科学, 40(8): 1396-1402. doi: 10.3799/dqkx.2015.124
    引用本文: 刘道飞, 陈圣波, 陈磊, 马明, 2015. 以SiO2含量为辅助因子的ASTER热红外遥感硅化信息提取. 地球科学, 40(8): 1396-1402. doi: 10.3799/dqkx.2015.124
    Liu Daofei, Chen Shengbo, Chen Lei, Ma Ming, 2015. Silicification Information Extraction Based on the Content of SiO2 from ASTER TIR Data. Earth Science, 40(8): 1396-1402. doi: 10.3799/dqkx.2015.124
    Citation: Liu Daofei, Chen Shengbo, Chen Lei, Ma Ming, 2015. Silicification Information Extraction Based on the Content of SiO2 from ASTER TIR Data. Earth Science, 40(8): 1396-1402. doi: 10.3799/dqkx.2015.124

    以SiO2含量为辅助因子的ASTER热红外遥感硅化信息提取

    doi: 10.3799/dqkx.2015.124
    基金项目: 

    国家高技术研究发展计划(863计划)项目 2012AA12A308

    中国地质调查局项目 1212011220469

    国家自然科学基金项目 41402293

    详细信息
      作者简介:

      刘道飞(1989-), 男, 硕士研究生, 主要从事热红外遥感研究.E-mail: liudaofei1989@126.com

      通讯作者:

      马明, E-mail: maming521521@163.com

    • 中图分类号: P627

    Silicification Information Extraction Based on the Content of SiO2 from ASTER TIR Data

    • 摘要: 硅化蚀变是岩石矿物蚀变中一种重要的矿化蚀变类型, 与很多金矿的形成有着密切的关系, 且硅化信息常作为野外重要的找矿标志.通过分析硅化蚀变矿物在先进星载热发射和反射辐射仪(advanced spaceborne thermal emission and reflection radiometer, ASTER)热红外波段的光谱特征, 依据硅化作用与SiO2含量间的关系, 选取了SiO2含量作为提取硅化信息的辅助因子, 提出了ASTER热红外遥感硅化信息提取方法.以内蒙古二连浩特市北部地区为例, 完成了该地区硅化信息提取工作.通过野外实地勘察验证, 发现在39个野外实地硅化蚀变采样点中33个采样点在蚀变图像中得到验证, 精度达到86.14%.

       

    • 硅化蚀变作用是指岩石在热液作用下,产生含有石英、玉髓、蛋白石、似碧玉等蚀变矿物的过程,是岩石矿物蚀变中一种重要的矿化蚀变类型,且石英脉和硅化带等硅化信息也常作为野外重要的找矿标志(于正军等,2010).很多金矿的形成与硅化蚀变密切相关,一般硅化越强,矿化越显著(肖晋等,1989).

      岩石在发生硅化作用后形成的蚀变矿物,其主要成分为SiO2,由此岩石在发生硅化蚀变作用前后SiO2含量会发生巨大的变化,而且蚀变程度越强,岩石所含SiO2含量越高.这就使得发生硅化作用的岩石在SiO2含量方面成为一种异常,偏离相对均一的地质体背景,所以依据SiO2含量可以用来辅助对硅化信息的提取.

      自从先进星载热发射和反射辐射仪(advanced spaceborne thermal emission and reflection radiometer,ASTER)传感器问世以来,国内外诸多学者针对造岩矿物在ASTER热红外波段的发射光谱特征与SiO2含量的关系做了大量研究,Ninomiya(2002)实现了从定性到定量的突破;Lyon(1965)测得25个火成岩样品(酸性-超基性)的发射率光谱及其SiO2含量,结果表明随着SiO2含量的增加,发射率吸收位置会向短波方向偏移.在国内,杨长保和朱群(2009)从USGS、JPL和JHU 3个岩石矿物光谱库中选择62个火成岩样本,并利用SPSS软件进行多元逐步回归分析从而得出了定量反演SiO2含量的回归方程.在硅化信息提取研究方面,利用遥感可见光-近红外数据对其进行提取的技术方法已日趋成熟,杨波等(2005)利用Landsat/TM数据建立了矿化信息定量提取模型,对鹰嘴山金矿区的硅化信息进行了提取.

      前人对于SiO2含量反演及硅化信息提取是分开独立研究,并没有考虑依据二者间的关系.利用SiO2含量和热红外遥感数据来对硅化信息的提取的研究还尚未见报道.因此,本文选取ASTER热红外数据,对蚀变矿物在ASTER热红外波段的光谱特征进行分析,以SiO2含量为辅助因子,寻求硅化信息的提取方法.

      ASTER热红外数据包含5个波段(波段10~14),波段范围为8.125~11.650 μm,空间分辨率为90 m,辐射度分辨率≤0.3 k,量化级别为12 bits.

      ASTER热红外数据的预处理包括辐射定标、大气校正、几何校正以及发射率反演.笔者基于TIR AtmosCorrect模型对ASTER热红外数据进行大气校正,消除了大气对地物发射率信息的影响;同时在ASTER可见光-近红外波段影像上选取与基准图像相对应的控制点,通过控制点坐标信息找到对应于热红外波段影像上的点,完成了ASTER热红外数据的几何校正.利用发射率归一化法对ASTER热红外波段数据进行了发射率反演(徐州等,2006),得到了ASTER的10~14波段发射率影像图.

      在研究ASTER热红外数据硅化信息提取方法的过程中,对USGS标准波谱数据库中硅化蚀变矿物石英及蛋白石发射率曲线重采样到ASTER热红外波段,综合分析后,笔者发现石英和蛋白石在9.0 μm处有明显的吸收谷(对应ASTER数据的B12),在10.65 μm处有明显的发射峰特征(对应ASTER数据的B13)(图 1).

      图  1  硅化蚀变矿物ASTER热红外发射率曲线
      Fig.  1.  The ASTER TIR emissivity spectral of silicification minerals

      依据蚀变矿物在ASTER热红外波段的发射率曲线特征,结合Rowan and Mars(2003)针对ASTER热红外数据使用的比值法,笔者发现可以突出硅化信息,对硅化信息按照蚀变强度分为3级进行初步提取.

      在对硅化信息初步提取的基础上,依据发生硅化作用的岩石在SiO2含量大小上明显高于其他地质背景,同时硅化蚀变越强SiO2含量越高.依据此,可以将SiO2含量定义为提取硅化信息的辅助因子,对初步提取的硅化信息做进一步纠正,以完善ASTER热红外遥感硅化信息提取方法,方法具体流程如图 2所示.

      图  2  ASTER热红外遥感硅化信息提取方法流程
      Fig.  2.  The procedure of method for ASTER TIR silicification information extraction

      为验证提出的ASTER热红外硅化信息提取方法是否具有可行性,笔者选取两处成像于2005年9月15日的ASTER热红外影像,在内蒙古二连浩特市北部区域做实地验证分析.

      研究区位于内蒙古二连浩特市北部,地理坐标范围111°45′00″~112°15′00″E,44°10′00″~44°45′00″N.由于研究区地处戈壁荒漠区,植被覆盖度较低,岩石裸露条件较好,所以被选作为低覆盖草原蚀变信息提取方法及地质填图方法研究区.

      研究区出露地层以石炭系-二叠系为主,零星发育中、晚奥陶世、早石炭世、晚石炭世和晚侏罗世的地层,多为第四系覆盖.区内褶皱构造和断裂构造同等发育,所显示的构造线方向以NE方向为主.区内岩浆岩分布相对广泛,以华力西期和燕山期的花岗岩类侵入岩石为主.华力西期岩体侵入石炭系宝力高庙组或中、下奥陶统乌宾敖包组中,以二长花岗岩、正长花岗岩、砂砾岩以及黑云母花岗闪长岩为主,华力西晚期的岩体被侏罗纪地层不整合覆盖,或被燕山期花岗岩侵入(图 3).

      图  3  研究区地质图
      1.湖积层:现代湖积淤泥沉积;2.冲洪积层: 由砂岩及砾石层组成;3.冲洪积层: 由粗砂和砾石层组成;4.阿巴嘎组: 灰紫色、紫褐色气孔状玄武岩及伊丁玄武岩、安山玄武岩;5.宝格达乌拉组:砖红色泥岩夹含砾粗砂岩;6.伊尔丁曼哈组:红色粘土和黄色砂砾岩;7.大磨拐河组:页岩、泥岩、砂岩、砂砾岩及砾岩组成夹褐煤;8.白音高老组:流纹质岩屑晶屑凝灰岩、流纹岩及流纹质溶结凝灰岩等酸性火山岩;9.玛尼吐组:安山岩、粗安岩、石英粗安岩、安山玢岩安山质角砾凝灰岩、英安玢岩及灰黑色及灰紫色玄武岩;10.宝力高庙组二段:灰-灰褐色安山岩、溶解凝灰岩,黄褐色灰绿色流纹质含角砾晶屑凝灰岩、凝灰质砂岩英安岩等中酸性火山岩及火山碎屑岩,含植物化石;11.宝力高庙组二段:(变质)长石砂岩、板岩、砾岩、硬砂岩为主夹中酸性岩屑晶屑凝灰岩、安山玢岩及灰岩透镜体;12.泥鳅河组:为浅海相碎屑岩夹灰岩组合,岩性为灰色、灰绿色、褐灰色(变质)粉砂质泥、砂岩夹灰岩;13.乌宾敖包组:灰褐色、灰绿色板岩、绢云母板岩、分砂质板岩、长石砂岩、变泥岩夹安山玢岩及灰岩透镜体;14.肉红色中细粒正长花岗岩;15.肉红色中细粒云母二长花岗岩;16.肉红色斑状中细粒黑云母二长花岗岩;17.肉红色中细粒-细中粒碱长花岗岩;18.灰绿色角闪花岗闪长岩;19.灰绿-暗灰绿色中细粒闪长岩、石英闪长岩;20.灰、灰白、灰黄色中细粒正长花岗岩、碱长花岗岩局部有闪石碱性花岗岩;21.灰白-灰色斑状中细粒黑云母正长花岗岩;22.肉红色中细粒-细粒二长花岗岩;23.灰白-灰色中细粒黑云母二长花岗岩;24.灰白-灰色斑状中细粒黑云母二长花岗岩及少量花岗闪长岩;25.灰白-灰色花岗闪长岩及少量石英闪长岩、石英二长岩;26.晚侏罗世次粗面斑岩;27.石英脉;28.花岗岩脉;29.碱性花岗斑岩脉;30.闪长岩脉;31.地质界线;32.实测地层不整合界线、火山喷发不整合地质界线;33.实测正断层;34.实测逆断层;35.实测平移断层;36.实测性质不明断层;37.角岩化;38.构造破碎带
      Fig.  3.  Geological sketch of the research area

      陈江和王安建(2007)将ASU波谱库的矿物波谱重采样至ASTER热红外波段,对矿物的波谱进行波段比值处理,选择波段比值与SiO2含量最大相关系数,进行对数模拟,从而确定了发射率波谱与SiO2含量的数值关系.其反演公式如下:

      $$ \begin{array}{l} \;\;\;\;\;\;\;{\rm{Si}}{{\rm{O}}_2} = 28.76 \times {\rm{ln(}}6.56 \times {E_{13}} \times {E_{14}}/({E_{10}} \times \\ {\mathit{E}_{12}}){\rm{), }} \end{array} $$ (1)

      式(1)中,E10E12E13E14对应ASTER热红外波段中心波长处的发射率值.笔者依据公式(1),在对ASTER热红外数据做完预处理的基础上,利用波段计算工具得出了研究区的SiO2含量结果(图 4).图 4中含量高值对应于地质图上的红色粘土和黄色砂砾岩,含量低值对应于地质图上的安山岩及粗安岩等.

      图  4  研究区SiO2含量
      Fig.  4.  The sketch showing the content of SiO2

      利用比值法(B13/B12)对研究区发射率影像进行处理,得到突出硅化信息的图像.对其进行滤波处理后,统计滤波后图像的均值和标准差,对初步提取的硅化信息按蚀变强度分为3级.硅化信息初步结果如图 5所示.

      图  5  硅化蚀变信息初步结果
      Fig.  5.  The original silicification information results of study area

      对反演得到的SiO2含量结果与初步提取的硅化信息结果进行叠加对比分析,将图像中SiO2含量明显高于周围背景而没有在蚀变结果中得到体现的区域,认定为存在硅化蚀变.由此得到研究区的最终硅化信息结果(图 6).对比分析图 5图 6,可以明显地发现在区域A~D中,通过引入SiO2含量这一辅助因子使得在最终得到的结果图中硅化信息得以补充完善.尤其是在区域AC中,完善的硅化信息在野外实地勘察中得到了实地验证.

      图  6  硅化蚀变信息最终结果
      Fig.  6.  The ultimate silicification information results of study area

      通过实验区野外考察,笔者发现在研究区裸露岩石处硅化蚀变现象十分明显,硅化程度由弱到强,在39个野外硅化采样点中,33个采样点在结果图中得到验证,精度达到86.14%.然而对于初始硅化信息提取结果,其没有考虑SiO2含量与硅化信息间的关系,这就使得有些野外实际存在的岩石硅化蚀变在结果图中没有被提取出来,提取结果精度为76.92%.总的来说,通过以SiO2含量为辅助因子对硅化信息纠正,使得硅化信息提取结果精度提高了9.22%.

      通过分析硅化蚀变矿物在ASTER热红外波段的波谱特征,利用比值法(B13/B12)得到初步硅化信息结果,结合硅化蚀变与SiO2含量间的关系对硅化信息进行纠正,完善了提取方法.通过在内蒙古二连浩特市北部区域对方法进行实地应用,笔者发现硅化信息提取结果与野外实地勘察吻合度较高,精度达到86.14%,与未利用SiO2含量来对硅化信息进行纠正的结果相比,精度提高了9.22%.结果表明,以SiO2含量为辅助因子的硅化信息提取方法能够实现对野外硅化信息进行较精确的提取,丰富了利用遥感技术对硅化信息提取的手段,对遥感地质找矿工作带来了帮助.

    • 图  1  硅化蚀变矿物ASTER热红外发射率曲线

      Fig.  1.  The ASTER TIR emissivity spectral of silicification minerals

      图  2  ASTER热红外遥感硅化信息提取方法流程

      Fig.  2.  The procedure of method for ASTER TIR silicification information extraction

      图  3  研究区地质图

      1.湖积层:现代湖积淤泥沉积;2.冲洪积层: 由砂岩及砾石层组成;3.冲洪积层: 由粗砂和砾石层组成;4.阿巴嘎组: 灰紫色、紫褐色气孔状玄武岩及伊丁玄武岩、安山玄武岩;5.宝格达乌拉组:砖红色泥岩夹含砾粗砂岩;6.伊尔丁曼哈组:红色粘土和黄色砂砾岩;7.大磨拐河组:页岩、泥岩、砂岩、砂砾岩及砾岩组成夹褐煤;8.白音高老组:流纹质岩屑晶屑凝灰岩、流纹岩及流纹质溶结凝灰岩等酸性火山岩;9.玛尼吐组:安山岩、粗安岩、石英粗安岩、安山玢岩安山质角砾凝灰岩、英安玢岩及灰黑色及灰紫色玄武岩;10.宝力高庙组二段:灰-灰褐色安山岩、溶解凝灰岩,黄褐色灰绿色流纹质含角砾晶屑凝灰岩、凝灰质砂岩英安岩等中酸性火山岩及火山碎屑岩,含植物化石;11.宝力高庙组二段:(变质)长石砂岩、板岩、砾岩、硬砂岩为主夹中酸性岩屑晶屑凝灰岩、安山玢岩及灰岩透镜体;12.泥鳅河组:为浅海相碎屑岩夹灰岩组合,岩性为灰色、灰绿色、褐灰色(变质)粉砂质泥、砂岩夹灰岩;13.乌宾敖包组:灰褐色、灰绿色板岩、绢云母板岩、分砂质板岩、长石砂岩、变泥岩夹安山玢岩及灰岩透镜体;14.肉红色中细粒正长花岗岩;15.肉红色中细粒云母二长花岗岩;16.肉红色斑状中细粒黑云母二长花岗岩;17.肉红色中细粒-细中粒碱长花岗岩;18.灰绿色角闪花岗闪长岩;19.灰绿-暗灰绿色中细粒闪长岩、石英闪长岩;20.灰、灰白、灰黄色中细粒正长花岗岩、碱长花岗岩局部有闪石碱性花岗岩;21.灰白-灰色斑状中细粒黑云母正长花岗岩;22.肉红色中细粒-细粒二长花岗岩;23.灰白-灰色中细粒黑云母二长花岗岩;24.灰白-灰色斑状中细粒黑云母二长花岗岩及少量花岗闪长岩;25.灰白-灰色花岗闪长岩及少量石英闪长岩、石英二长岩;26.晚侏罗世次粗面斑岩;27.石英脉;28.花岗岩脉;29.碱性花岗斑岩脉;30.闪长岩脉;31.地质界线;32.实测地层不整合界线、火山喷发不整合地质界线;33.实测正断层;34.实测逆断层;35.实测平移断层;36.实测性质不明断层;37.角岩化;38.构造破碎带

      Fig.  3.  Geological sketch of the research area

      图  4  研究区SiO2含量

      Fig.  4.  The sketch showing the content of SiO2

      图  5  硅化蚀变信息初步结果

      Fig.  5.  The original silicification information results of study area

      图  6  硅化蚀变信息最终结果

      Fig.  6.  The ultimate silicification information results of study area

    • [1] Chen, J., Wang, A.J., 2007. The Pilot Study on Petrochemistry Components Mapping with ASTER Thermal Infrared Remote Sensing Data. Journal of Remote Sensing, 11(4): 1-3 (in Chinese with English abstract). http://www.researchgate.net/publication/284923747_The_pilot_study_on_petrochemistry_components_mapping_with_ASTER_thermal_infrared_remote_sensing_data
      [2] Lyon, R.J.P., 1965. Analysis of Rocks by Spectral Infrared Emission (8 to 25 Microns). Economic Geology, 60(4): 715-736. doi: 10.2113/gsecongeo.60.4.715
      [3] Ninomiya, Y., 2002. Mapping Quartz, Carbonate Minerals, and Mafic-Ultramafic Rocks Using Remotely Sensed Multispectral Thermal Infrared ASTER Data. Proc. SPIE 4710, Thermosense XXIV, 191. doi: 10.1117/12.459566
      [4] Rowan, L.C., Mars, J.C., 2003. Lithologic Mapping in the Mountain Pass, California Area Using Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Data. Remote Sensing of Environment, 84(3): 350-366. doi: 10.1016/S0034-4257(02)00127-X
      [5] Xiao, J., Han, L.F., Tang, D.L., 1989. The Genesis of Silicified Zones in Southern Hunan and Their Ore-Prospecting Significance. Journal of Ore Deposit Geology, 8(4): 9 (in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-KCDZ198904009.htm
      [6] Xu, Z., Zhao, H.J., Li, N., 2006. Study of Inversing Land Surface Emissivity from ASTER Data. Infrared and Laser Engineering, 35(Suppl. ): 517-522(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-HWYJ2006S4097.htm
      [7] Yang, B., Wu, D.W., Lai, J.Q., 2005. The Set Up of Quantization Remote Sensing Models for Mineralization and a Literation—Take Silicification in Yingzuishan Gold Mine as an Example. Journal of Remote Sensing, 9(6): 717-724(in Chinese with English abstract). http://www.oalib.com/paper/1470905
      [8] Yang, C.B., Zhu, Q., 2009. Quantitative Inversion of SiO2 Content in the Surface Rock from ASTER Thermal Infrared Data. Geology and Exploration, 6(11): 692-694(in Chinese with English abstract).
      [9] Yu, Z.J., Xu, S.M., Wang, J.D., et al., 2010. A New Type of Reservoir of Paleozoic Buried Hill in Zhuanghai Area: Discovery and Characteristics of Silicified Rock. Earth Science—China University of Geosciences, 35(1): 87-96 (in Chinese with English abstract). doi: 10.3799/dqkx.2010.009
      [10] 陈江, 王安建, 2007. 利用ASTER热红外遥感数据开展岩石化学成分填图的初步研究. 遥感学报, 11(4): 1-3. https://www.cnki.com.cn/Article/CJFDTOTAL-YGXB200704024.htm
      [11] 肖晋, 韩兰生, 汤大立, 1989. 论湘南地区"硅化带"的成因及其找矿意义. 矿床地质, 8(4): 9. https://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ198904009.htm
      [12] 徐州, 赵慧洁, 李娜, 2006. 基于ASTER数据的地物光谱比辐射率的反演研究. 红外与激光工程, 35(增刊): 517-522. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ2006S4097.htm
      [13] 杨波, 吴德文, 赖健清, 等, 2005. 矿化信息提取定量遥感模型的建立——以鹰嘴山硅化蚀变为例. 遥感学报, 9(6): 717-724. https://www.cnki.com.cn/Article/CJFDTOTAL-YGXB200506012.htm
      [14] 杨长保, 朱群, 2009. ASTER热红外遥感地表岩石的二氧化硅含量定量反演. 地质与勘探, 6(11): 692-694. https://www.cnki.com.cn/Article/CJFDTOTAL-DZKT200906010.htm
      [15] 于正军, 许淑梅, 王金铎, 等, 2010. 桩海地区下古生界潜山内新一类储层: 硅化岩储层的发现及特征. 地球科学——中国地质大学学报, 35(1): 87-96. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX201001012.htm
    • 期刊类型引用(9)

      1. 汪大明,栗旭升,魏佳林,曹思琦,汪翡翠,仝云霄,闫国强. 虎头山地区多平台同步热红外遥感稀有金属岩脉识别. 地球科学. 2024(06): 2242-2252 . 本站查看
      2. 裴秋明,沈家乐,王世明,房大任,高永璋,李典,马少兵. 多源遥感卫星数据在脉状萤石矿床中的找矿预测应用:以内蒙古水头萤石矿床为例. 西北地质. 2024(04): 121-134 . 百度学术
      3. 曹海玲,王正海,秦昊洋,孙袁超,尹国盼. 基于加权融合TIR数据的LST反演与硅化信息提取——以广东仁差盆地铀多金属矿区为例. 遥感学报. 2023(07): 1691-1701 . 百度学术
      4. 茹菲娜·阿力木江,陈川,高玲玲,李顺达. 基于ASTER数据的遥感蚀变信息提取——以西天山穆龙套地区为例. 地球学报. 2022(02): 235-245 . 百度学术
      5. 张弘,高鹏鑫,高卿楠. 热红外反射光谱技术在石英含量评价中的应用. 岩矿测试. 2021(05): 710-719 . 百度学术
      6. 代晶晶,赵龙贤,姜琪,王海宇,刘婷玥. 热红外高光谱技术在地质找矿中的应用综述. 地质学报. 2020(08): 2520-2533 . 百度学术
      7. 张元涛,余长发,潘蔚,田青林. WorldView-3影像与ASTER热红外影像在内蒙古卫境地区铀矿勘查中的应用. 铀矿地质. 2020(05): 408-417 . 百度学术
      8. 王玉山,刘玉红,周勇,鲁立辉,陈鑫,许荣科,郑有业. 基于ASTER热红外数据的SiO_2及岩性类别信息提取——以青海省绿梁山—双口山地区为例. 地质找矿论丛. 2019(04): 596-600 . 百度学术
      9. 郭帮杰,张杰林,武鼎. 热红外高光谱遥感回归分析定量反演石英含量. 科学技术与工程. 2018(17): 125-130 . 百度学术

      其他类型引用(5)

    • 加载中
    图(6)
    计量
    • 文章访问数:  2814
    • HTML全文浏览量:  174
    • PDF下载量:  336
    • 被引次数: 14
    出版历程
    • 收稿日期:  2015-03-14
    • 刊出日期:  2015-08-01

    目录

    /

    返回文章
    返回