高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

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

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

基于ZEMAX二氧化硫荧光采集光路的设计仿真

刘杰辉 李鑫

引用本文:
Citation:

基于ZEMAX二氧化硫荧光采集光路的设计仿真

    作者简介: 刘杰辉(1968-), 男, 教授, 研究生导师, 主要研究方向为机械制造及自动化。E-mail:ljsgz163@163.com.
  • 中图分类号: TN202

Design and simulation of fluorescence acquisition optical path of sulfur dioxide based on ZEMAX

  • CLC number: TN202

  • 摘要: 为了实现二氧化硫分析仪精确连续地检测空气中低含量SO2, 针对传统二氧化硫检测仪在荧光汇聚时存在光线汇聚效果弱、像差大且由于透镜的焦距长增大了荧光采集光路长度等问题, 在传统仪器的基础上对荧光采集光路重新进行了设计优化。采用分解光焦度的方法改善系统的球差问题, 进行了理论分析和实验验证。改进后对其进行ZEMAX仿真并运用正交法对新光路进行修正。通过仿真两种光路的非序列结构阴影模型图、探测查看器的光照图以及点列图, 取得了光强以及光斑直径的实验数据。结果表明, 优化后的光路系统峰辐强度达到了219.41W/cm2, 弥散斑直径减少了17%。该方案能够有效解决传统光路中的不足之处, 成像质量效果更佳。
  • Figure 1.  Effect of spherical aberration

    Figure 2.  Imaging of infinite object point due to single refractive surface

    Figure 3.  Bispherical refraction system

    Figure 4.  Fluorescence acquisition system with dual convex lens

    Figure 5.  Acquisition system with symmetrical plano-convex lens

    Figure 6.  Orthogonal optimization diagram

    图 7  NSC shadow model diagram of double convex lens

    Figure 8.  NSC shadow model diagram of symmetrical plano-convex lens

    Figure 9.  View of double convex lens detector

    Figure 10.  View of symmetrical plano-convex lens detector

    Figure 11.  Point series diagram of double convex lens detector

    Figure 12.  Point series diagram of symmetrical plano-convex lens detector

    Table 1.  Comparison diagram of probe view data

    peak intensity of radial/(W·cm-2) total power/W hit the number
    biconvex lens 73.112 0.29570 2957
    symmetrical plano-convex lens 219.41 0.35410 3541
    下载: 导出CSV
  • [1]

    YANG Ch P, HE Zh J, YANG G L, et al.Study on optical properties ofatmospheric sulfur dioxide detection reaction chamber[J]. Optical Technique, 2003, 29(1): 39-41(in Chinese).
    [2]

    LI M M. Research on key technologies of optical detection of sulfur dioxide concentration[D]. Qinhuangdao: Yanshan University, 2014: 9-17(in Chinese).
    [3]

    LI X X, ZHANG Zh H, ZHANG X Y, et al. Non-contact thickness measurement of optical element based on astigmatism[J]. Laser Technology, 2019, 43(6):741-746(in Chinese).
    [4]

    LI G H. Day and night CCTV wide-angle zoom lens design[D]. Changchun: Changchun University of Science and Technology, 2013: 4-9(in Chinese).
    [5]

    WANGA L, ZHANG Y G, ZHOU X, et al. Optical sulfur dioxide sensor based on broadband absorption spectroscopy in the wavelength range of 198-222nm[J]. Sensors & Actuators, 2017, B241:146-150.
    [6]

    GHOSH A, KAREDLA N, THIELE J C, et al.Fluorescence lifetime correlation spectroscopy: Basics and applications[J]. Methods, 2018, 140/141:32-39. doi: 10.1016/j.ymeth.2018.02.009
    [7]

    AARóN C, ASTICIO V, FABIáN A, et al. Geometrical-phase lens based optical system for the spin-splitting of vector beams[J]. Optics and Lasers in Engineering, 2018, 110:401-409. doi: 10.1016/j.optlaseng.2018.06.017
    [8]

    GAO Y L. Research on optimization design of aspherical aspheric lens based on Zemax[D]. Lanzhou: Lanzhou University of Technology, 2017: 19-39(in Chinese).
    [9]

    WU X H, CHEN Y, QIU J H. Simulation analysis of deformed newton ring device based on ZEMAX[J]. University Physics Experiment, 2017, 30(1): 111-114(in Chinese).
    [10]

    CHANG W, ONO Y, KUMEMURA M, et al. On-line determination of trace sulfur dioxide in air by integrated microchip coupled with fluorescence detection[J]. Talanta, 2005, 67(3):646-650. doi: 10.1016/j.talanta.2005.03.017
    [11]

    MICHEL A P M, KAPIT J. Deep ultraviolet light emitting diode (LED)-based sensing of sulfur dioxide[J]. Applied Spectroscopy, 2017, 71(5):996-1003. doi: 10.1177/0003702816665126
    [12]

    PAREDES B L, ARAúJO H M, FROBORG F, et al. Response of photomultiplier tubes to xenon scintillation light[J]. Astroparticle Physics, 2018, 102:56-66. doi: 10.1016/j.astropartphys.2018.04.006
    [13]

    HANULIA T, INAMI W, ONO A, et al. Fluorescence lifetime measurement excited with ultraviolet surface plasmon resonance[J]. Optics Communications, 2018, 427:266-270. doi: 10.1016/j.optcom.2018.06.069
    [14]

    WU F, LI S S, ZHU H B, et al. Analysis on the optical axis error of the spherical shell in the electro-optical system[J]. Optik, 2018, 168:458-461. doi: 10.1016/j.ijleo.2018.04.093
    [15]

    CHEN X, YOU L B, YIN G Y.Pulse width measurement with ultraviolet femtosecond laser[J]. Laser Technology, 2019, 43(2):205-211(in Chinese).
  • [1] 苏杭王桂梅张振兴郭炜 . 大气SO2检测模块的光路研究. 激光技术, 2020, 44(2): 266-270. doi: 10.7510/jgjs.issn.1001-3806.2020.02.023
    [2] 赵琦樊红英李轶国蒋泽伟胡绍云赖庚辛黄燕琳耿旭 . 大能量中空光束大气传输的仿真与实验比对研究. 激光技术, 2014, 38(4): 542-545. doi: 10.7510/jgjs.issn.1001-3806.2014.04.021
    [3] 曹向广周建忠袁国定杜建钧倪敏雄 . 基于有限元模拟的激光喷丸成形正交优化设计. 激光技术, 2007, 31(6): 656-658,662.
    [4] 刘兵王巨胜杨泽后李晓锋樊冬任鹏李斌罗雄冯力天 . 基于非视线红外激光大气散射通信技术研究. 激光技术, 2014, 38(6): 854-858. doi: 10.7510/jgjs.issn.1001-3806.2014.06.028
    [5] 杨军王慧章曦 . 大气湍流对部分相干激光瑞利区间影响的研究. 激光技术, 2016, 40(3): 456-460. doi: 10.7510/jgjs.issn.1001-3806.2016.03.033
    [6] 赵琦郝红宇樊红英李建欣蒋泽伟肖星 . 部分相干cosh-Gaussian光束通过大气湍流后的聚焦特性. 激光技术, 2016, 40(5): 750-755. doi: 10.7510/jgjs.issn.1001-3806.2016.05.028
    [7] 刘玉丽 . 探测边界层大气温度的转动喇曼激光雷达. 激光技术, 2018, 42(4): 541-544. doi: 10.7510/jgjs.issn.1001-3806.2018.04.020
    [8] 刘琳琳杨杰黄见苑克娥尹凯欣胡顺星 . 激光雷达观测淮南大气SO2和NO2浓度廓线实例分析. 激光技术, 2019, 43(3): 353-358. doi: 10.7510/jgjs.issn.1001-3806.2019.03.012
    [9] 石玉立杨丰恺曹念文 . MPL反演南京北郊气溶胶光学厚度准确度的研究. 激光技术, 2015, 39(3): 372-376. doi: 10.7510/jgjs.issn.1001-3806.2015.03.020
    [10] 赵延仲孙华燕宋丰华赵琳峰 . 离焦量对“猫眼效应”反射特性的影响. 激光技术, 2008, 32(1): 71-74.
    [11] 刘秋武陈亚峰王杰王晓宾曹开法黄见胡顺星 . 差分吸收NO2激光雷达光源的设计与实现. 激光技术, 2018, 42(4): 433-439. doi: 10.7510/jgjs.issn.1001-3806.2018.04.001
    [12] 王玉诏 . 一种用于激光雷达信号的自适应分段平滑算法. 激光技术, 2014, 38(5): 599-602. doi: 10.7510/jgjs.issn.1001-3806.2014.05.005
    [13] 汤明玥李宾中 . 双曲余弦高斯光束在非Kolmogorov湍流中的湍流距离. 激光技术, 2015, 39(4): 581-584. doi: 10.7510/jgjs.issn.1001-3806.2015.04.034
    [14] 陈莉英姜洪波沈利沣赵志刚 . 脉冲腔衰荡技术探测气溶胶消光系数的研究. 激光技术, 2016, 40(2): 241-244. doi: 10.7510/jgjs.issn.1001-3806.2016.02.019
    [15] 周正兰袁扬胜束杰徐翔屈军 . 部分相干月牙形光束在非Kolmogorov谱中的漂移. 激光技术, 2019, 43(4): 579-584. doi: 10.7510/jgjs.issn.1001-3806.2019.04.025
    [16] 杨建清罗积军徐军孙琦云王婷 . 基于雨滴谱分布数据库的近红外辐射衰减研究. 激光技术, 2018, 42(2): 161-165. doi: 10.7510/jgjs.issn.1001-3806.2018.02.004
    [17] 闫顺生李海花 . 利用喇曼-米方法确定激光雷达几何因子. 激光技术, 2013, 37(4): 511-514. doi: 10.7510/jgjs.issn.1001-3806.2013.04.021
    [18] 郑卫峰李步洪谢树森林永钟周川钊 . 荧光法成像系统信噪比的实验研究. 激光技术, 2006, 30(2): 123-125.
    [19] 李蕾臧景峰 . 双狭缝扫描法测量激光光束质量. 激光技术, 2015, 39(6): 845-849. doi: 10.7510/jgjs.issn.1001-3806.2015.06.024
    [20] 张恩涛季小玲吕百达 . 内光路中大气吸收对远场光束质量的影响. 激光技术, 2006, 30(1): 96-98.
  • 加载中
图(12) / 表(1)
计量
  • 文章访问数:  428
  • HTML全文浏览量:  316
  • PDF下载量:  2
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-04-22
  • 录用日期:  2019-05-29
  • 刊出日期:  2020-03-25

基于ZEMAX二氧化硫荧光采集光路的设计仿真

    作者简介: 刘杰辉(1968-), 男, 教授, 研究生导师, 主要研究方向为机械制造及自动化。E-mail:ljsgz163@163.com
  • 河北工程大学 机械与装备工程学院, 邯郸 056038

摘要: 为了实现二氧化硫分析仪精确连续地检测空气中低含量SO2, 针对传统二氧化硫检测仪在荧光汇聚时存在光线汇聚效果弱、像差大且由于透镜的焦距长增大了荧光采集光路长度等问题, 在传统仪器的基础上对荧光采集光路重新进行了设计优化。采用分解光焦度的方法改善系统的球差问题, 进行了理论分析和实验验证。改进后对其进行ZEMAX仿真并运用正交法对新光路进行修正。通过仿真两种光路的非序列结构阴影模型图、探测查看器的光照图以及点列图, 取得了光强以及光斑直径的实验数据。结果表明, 优化后的光路系统峰辐强度达到了219.41W/cm2, 弥散斑直径减少了17%。该方案能够有效解决传统光路中的不足之处, 成像质量效果更佳。

English Abstract

    • 迄今为止,空气中二氧化硫的污染相对较为严重。人们在关注经济快速发展的同时也开始注重环境的保护。为了更加准确地了解到我国现在环境的污染状况, 研制出能够实时在线检测的二氧化硫仪器势在必行。二氧化硫检测仪的原理是紫外荧光法,这是国家标准GB3095-2012中所规定的检测方法。紫外荧光法可以突破传统化学检测方法的种种缺陷,能够更准确地检测二氧化硫[1]。其原理是:SO2气体被波长范围为220nm附近的紫外光照射时,SO2分子就从基态转化为激发态,激发态的SO2是不稳定的,在其转化为基态的时候产生荧光,对其产生的荧光进行收集利用就可以得出空气中二氧化硫的含量[2]。荧光强度和光电倍增管输出的电信号除去一些干扰以后是呈线性关系的。所以,对光电倍增管输出的电信号进行一系列的放大、滤波等处理后,就可以得出二氧化硫的具体含量。紫外荧光法与传统方法相比存在许多优势,不仅能够在线实时检测,并且可以降低误差,提高检测的准确性。利用荧光光谱法对SO2含量进行检测的过程中,最小的检测范围已经达到了10-9数量级[3]

      目前研发的二氧化硫检测仪器采用的荧光采集光路是利用双凸透镜对荧光点光源进行了汇聚。该系统存在很多不足,不仅光线的汇聚效果不理想、存在较大像差, 而且由于透镜的焦距长导致增大了荧光采集光路的长度[4]。基于上述光路的种种缺陷,本文中提出了一种新型对称式平凸结构,可以有效解决上述光路存在的问题,并通过ZEMAX仿真验证了对称式平凸结构的合理性。

    • 在实际运用中几乎不可能达到完美的成像质量。由球差效果图可知,光线并不在光轴上聚于同一个焦点,而在镜头边缘入射的光线与光轴的交点要比接近光轴入射的光线与光轴的交点离镜头的距离近,如图 1所示。这是很常见的像差之一,称之为球差[5]。球差的大小取决于光线在入瞳上的高度。消除球差可以采用非球面,非球面虽然不能完全消除像差,但可以使像差达到最小[6]。减小球差的另外一个方法是将光焦度进行分解[7]

      Figure 1.  Effect of spherical aberration

    • 在光学系统中,想要将球差减小到最小,首先应该明白球差产生的原因是什么。从理论上来说,如果一个系统中不存在球差,那么系统中任意光线的光程是永远相等的。

      在建模过程中首先需要建立一个适当的坐标系。在光学系统中,取其中一条单面子午截线的对称中心作为原点建立直角坐标系。所以,无限远处的点光源发出的平行光照经过该系统后,选取经过单面子午截线上点P(xy)的远轴光线和近轴光线作为研究对象,如图 2所示。按照上述理论可知,近轴光线和远轴光线的光程应该相等。远轴光线经过点P(xy)交x轴于点b,交y轴于点a,近轴光线是与x轴重合的光线,与远轴光线共同交于点bnn′分别代表的是空气的折射率和光学系统的折射率;xPyP也就是P点的坐标值;f′是近轴光线经过的路程。对于任意一个给定的光学系统,如果这个系统不存在球差,那么一切经过这个特定光学系统的光线都应该在光轴上交于同一个点。也就是说,这两条光线的光程是恒等的。远轴光线光程的具体公式是:

      $ n x_{P}+n^{\prime} \sqrt{\left(f^{\prime}-x_{P}\right)^{2}+y_{P}^{2}} $

      (1)

      Figure 2.  Imaging of infinite object point due to single refractive surface

      近轴光线光程则为nf′,因此对于任意子午截线上的点,光程恒等的曲面是:

      $ n^{\prime} f^{\prime}=n x_{P}+n^{\prime} \sqrt{\left(f^{\prime}-x_{P}\right)^{2}+y_{P}^{2}} $

      (2)

      化简得:

      $ \left(1-\frac{n^{2}}{n^{\prime 2}}\right) x_{P}^{2}+y_{P}^{2}-2\left(1-\frac{n}{n^{\prime}}\right) f^{\prime} x_{P}=0 $

      (3)

      由(3)式可得:当f′一定时, n′>n为椭圆;n′/n=-1为抛物面;n′/n=1为平面(yP=0)。

      除去n′/n=-1的情况后,该方程仍然是二次的,所研究的光路系统nn′分别代表空气和石英的折射率,所以不符合上述条件。即将球面透镜非球面化以后可以达到消球差的目的。

      非球面上各个点的曲率半径均不相同,其面形由高次多项式决定。从光学角度而言,非球面指折射面为球面以外的所有面。设计者可以通过改变其曲率来获得理想的成像质量。

    • 图 3为双球面折射系统图。图中, u1表示入射张角,u2表示出射张角,h为光学系统入射点到主轴的距离,n1n2n3分别表示不同介质的折射率。由图 3可知,双球面系统光焦度F的表达式为:

      $ \begin{array}{c} F=F_{1}+F_{2}-F_{2} F_{1} \frac{d}{n_{2}}=\frac{n_{2}-n_{1}}{r_{1}}+ \\ \frac{n_{3}-n_{2}}{r_{2}}-\frac{\left(n_{2}-n_{1}\right)\left(n_{3}-n_{2}\right) d}{r_{1} r_{2} n_{2}} \end{array} $

      (4)

      Figure 3.  Bispherical refraction system

      式中, O1折射面的光焦度为F1, O2折射面的光焦度为F2d为两个折射面之间的距离;r为曲率半径,光从左向右传播,以球面和主光轴的交点为准,球面的球心在该点以左,则曲率半径为负,反之,球心在在该点以右,则曲率半径为正。r1为第1个折射面的曲率半径,r2为第2个折射面的曲率半径。

      如果是双凸透镜系统,石英透镜的折射率要大于空气折射率,则n1=n3代表空气中的折射率;n2表示石英透镜的折射率;曲率半径r1=-r2,且r1>0, 那么化简可得:

      $ F=2 \times \frac{n_{2}-n_{1}}{r_{1}}-\frac{\left(n_{3}-n_{2}\right)^{2} d}{r_{1}^{2} n_{2}} $

      (5)

      同理, 如果是对称式平凸透镜,则n1=n3代表石英透镜的折射率;n2表示空气的折射率;曲率半径r1=-r2,且r2>0, 化简得:

      $ F=2 \times \frac{n_{3}-n_{2}}{r_{2}}-\frac{\left(n_{3}-n_{2}\right)^{2} d}{r_{2}^{2} n_{2}} $

      (6)

      由于两个光路都对光线进行了汇聚,所以光焦度F的值应该大于0。在本次设计过程中,对称式平凸透镜和双凸透镜的曲率半径相同;材料和波长相同,所以空气和石英透镜的折射率是相同的;双凸透镜的d值明显大于对称式平凸透镜的d值,所以F对称式平凸>F双凸。由理论可知,在光焦度大于0的情况下,F的值越大,汇聚效果越好。

      因为非球面价格昂贵、没有足够的精度以及随温度变化的不稳定性,故其在实际领域里难以普及。所以减少像差可以在光学系统中增加光学元件,光焦度可以在几个元件间进行分解,这样可以降低每个表面上的入射角度,从而达到的目的。

    • 目前的二氧化硫监测仪器采用一个双凸透镜和一个滤光片组成荧光采集光路,如图 4所示。由荧光发出的点光源经过双凸透镜进行汇聚,在通过滤光片滤光后入射到探测器上。为了减少像差而采用的对称式平凸结构,在两个透镜中设有可以限制较大入射角光线的隔镜环,在探测器前面设有滤光片,如图 5所示。将以上两种光路运用ZEMAX仿真进行对比。

      Figure 4.  Fluorescence acquisition system with dual convex lens

      Figure 5.  Acquisition system with symmetrical plano-convex lens

    • 设计或者优化任何一个光路结构,都有特定的要求,例如焦距、入瞳直径、视场、波长、材料、分辨率、渐晕和调制传递函数(modulation transfer function, MTF)等,根据不同系统的简易程度要求也各不相同[8]。由于荧光波段在330nm强度最高,所以设定的波长是330nm;在240nm~420nm的波长范围内,石英透镜对二氧化硫发射出荧光的透过率比普通玻璃透镜大很多,所以材料选用石英玻璃;其中视场为1,入瞳直径为20mm,分辨率为2048×2048。

      在确定光路的具体路径后对两个光路进行局部优化,得到最小的光斑半径。ZEMAX中局部优化算法有阻尼最小二乘法(damped least squares, DLS)和正交法(orthogonal descent, OD),在纯非序列系统中由于探测结果是在被像素化的探测器上获得,其评价函数是不连续的,使用阻尼最小二乘法的运用效果并不理想,故采用正交法更为合适。

      正交法对变量进行标准化,并对解空间进行离散化抽样来降低评价函数,并不计算评价函数的数值导数。相比较于阻尼最小二乘法,正交法对含有噪声评价系统的函数(比如非系列系统)更有优势。正交法能够很好的完成照度最大化、亮度增强以及均匀化的优化问题。图 6即为采用正交法后的优化结果。可以看到,将光斑半径明显减少达到了预期的优化效果。

      Figure 6.  Orthogonal optimization diagram

    • 以下是运用ZEMAX对两种光学系统进行仿真优化后的具体对比。

      (1) 首先进行两种系统的光路仿真,光线从光源发出经过透镜后汇聚于探测器上,如图 7图 8所示。对比两种光路的非序列结构(non-sequential compoonent, NSC)阴影模型图,观察到图 7中经过对称式双凸透镜入射到探测器上的光线较图 8中双凸透镜的汇聚效果更加集中,光线散射程度更小。对称式双凸透镜光路设计更加合理。

      图  7  NSC shadow model diagram of double convex lens

      Figure 8.  NSC shadow model diagram of symmetrical plano-convex lens

      (2) 为了更好地比较两种光路汇聚效果,设置光线条数为1000,跟踪分析光线探测器,得到探测器视图,如图 9图 10所示。对比两种光路的峰辐强度、总功率以及撞击次数(见表 1)[9],发现优化后的峰辐强度达到了219.41W/cm2, 是之前的3倍,总功率达到了0.35410W。由此可见,对称式平凸透镜的汇聚效果要比双凸透镜的汇聚效果好, 即在相同的荧光强度下,对称式平凸透镜可以更好地进行汇聚,以便于接收器接收到更多的信号[10]

      Figure 9.  View of double convex lens detector

      Figure 10.  View of symmetrical plano-convex lens detector

      Table 1.  Comparison diagram of probe view data

      peak intensity of radial/(W·cm-2) total power/W hit the number
      biconvex lens 73.112 0.29570 2957
      symmetrical plano-convex lens 219.41 0.35410 3541

      (3) 进行点列图仿真分析。在理想状态下,由一个点光源发出的所有光线在通过一个比较理想的光学系统之后,会汇聚到一个点,这就是这个物点的像点。[11]但是在实际光学系统的应用过程中,因为不可能不存在像差,所以经过光学系统用探测器查看后会形成一个光斑,称其为弥散斑。探测器上的视图称为点列图。

      点列图实际代表的是光斑分布的密集状况,越密集说明成像质量越好。衡量点列图光斑分布情况可以用几何最大半径值(geometric maximum radius, GEO)和均方根(root mean square, RMS)半径值来进行表示[12]。几何最大半径值就是以参考光线点为中心,包含所有光线的最大圆的半径,而均方根半径则是每条光线交点与参考光线点的距离的平方,除以光线条数后再开方[13]。也就是说均方根半径值反映了光能的集中程度,几何最大半径反映了像差的最大值。这两个的值越小,说明成像质量越好[14]

      图 11图 12是两种系统光学仿真后得到的点列图。由图 11可知, 对称式平凸透镜点列图的GEO为1182.71μm,RMS为671.872μm。图 12中双凸透镜点列图的GEO为6798.75μm,RMS为2501.57μm。比较两种系统数据,对称式平凸透镜的几何最大半径值和均方根半径值远小于双凸透镜点列图的值,两个系统的几何最大半径之比为:1182.71μm/6798.75μm=17%。所以对称式平凸透镜的弥散斑直径更小,成像质量更佳[15]

      Figure 11.  Point series diagram of double convex lens detector

      Figure 12.  Point series diagram of symmetrical plano-convex lens detector

      经过上述对比不难发现, 从NSC阴影模型图中光线的汇聚、探测器探测视图以及点列图中弥散斑半径大小均可得到对称式平凸透镜的成像质量要优于一个双凸透镜的结果。

    • 运用光学仿真软件ZEMAX,在目前监测仪器的荧光采集光路基础上对其进行了优化设计。对比探测视图和点列图后发现优化后的峰辐强度是之前的3倍,总功率达到了0.35410W,有效半径和均方根半径都大幅度减小,达到了优化的目的。在与光源距离相同的条件下,对称式平凸透镜相较于双凸透镜不仅结构简单、安装方便、光能损失较小而且成像质量相对较好。因此,对称式平凸透镜结构能够很好地应用于荧光的检测。

参考文献 (15)

目录

    /

    返回文章
    返回