-
在吸收光谱技术中, 使用光学腔增长激光与气体介质的作用路径, 可提升探测灵敏度. 然而, 高反射率腔镜会存在双折射效应, 导致光学腔产生两个本征偏振态, 入射光在两个偏振方向相移的不同会导致腔模的分裂, 会引起腔增强光谱信号以及腔衰荡光谱信号的扭曲. 本文观测到了双折射效应下腔增强信号的频率分裂现象, 并给出了函数模型, 拟合结果表明其可以准确得到透射腔模中不同偏振光的比例. 根据上述比例, 可给出考虑不同耦合效率、双折射效应下的腔衰荡信号模型, 实验结果表明相较于传统e指数模型, 本文模型可更精确描述腔衰荡信号, 得到拟合残差的标准偏差最大抑制了9倍. 该分析有利于改善腔衰荡信号信噪比和不确定性, 提升其浓度反演准确度.
-
关键词:
- Fabry-Pérot腔 /
- 双折射 /
- 腔衰荡光谱
In laser absorption spectroscopy, in order to improve gas detection sensitivity, optical cavity with high finesse is used to prolong the interaction path between the laser and the absorber. However, the birefringence of high reflectivity cavity mirrors generates two polarization eigenstates, and owing to the different phase shifts along the two directions, the cavity mode will be split. In this work, we first measure the cavity enhanced signal under birefringence and observe the mode split. And a model to mimic cavity enhanced spectroscopy under birefringent effect is presented, which can accurately fit the different polarization ratios at transmission. Finally, we propose a cavity ring-down signal model considering different coupling efficiencies of the two polarization directions of the cavity. Comparing with the conventional exponential model, the standard deviation of residual maximum suppression is as high as 9 times. And this analysis is helpful in improving the signal-to-noise ratio and uncertainty of cavity ring-down signal and increasing the accuracy of concentration inversion.-
Keywords:
- Fabry–Pérot cavity /
- birefringence /
- cavity ring-down spectroscopy
1. 引 言
Fabry-Pérot腔(FP腔)是光学谐振腔的一种, 由两个平行的高反镜组成, 具有体积小、结构稳定、便于调节等优点, 在物理学众多领域都有广泛的应用, 包括在激光器中用于提高光束的方向性和单色性[1]; 在光学钟中提供窄线宽的频率参考[2]; 在引力波探测器中用于增强信号[3].
在激光吸收光谱技术中, 光学腔可以增长激光与腔内气体介质的作用路径, 从而大大增强吸收信号, 提升探测灵敏度[4-7]. 其中路径增长因子等于2F/π, F是光学腔的精细度, 可以用腔镜的反射率r表示为
π√r/(1−r) . 随着近十年高反镜镀膜工艺的发展, 基于离子束溅射法(IBS)和单晶镀膜法(CMC)等技术, 近红外波段的反射率可以大于99.999%, 对应增长因子可以达到100万以上. 基于光学腔发展了腔增强直接吸收光谱技术(CEAS)[4,8]、积分腔输出光谱技术(ICOS)[9,10]、腔衰荡光谱技术(CRDS)[11-14]、噪声免疫腔增强光外差分子光谱技术(NICE-OHMS)[15]等, 对痕量气体的探测灵敏度最高可以达到10–14 cm–1, 最小可探测浓度达到百万亿分之一量级[16]. 其中CRDS通过测量腔内光场的衰减时间来反演气体浓度, 不受激光强度噪声的影响, 因此在系统复杂程度没有大幅增加的前提下, 可以获得很高的探测灵敏度, 同时, CRDS还具有浓度免定标的优势, 从而得到了广泛的应用.然而, 在FP腔中, 由于腔镜折射率的各向异性会引起双折射效应, 表现为腔镜存在两个本征偏振方向[17], 光沿两个偏振方向传播的相位延迟不同, 从而导致光学腔在这两个方向的自由光谱区以及纵模频率不再相同. 造成高反镜双折射效应的主要原因是残余应力中的本征应力, 涉及基底材料的微观形成机制和镀膜过程的化学反应, 其大小主要和镀膜沉淀工艺有关. 该现象在基于IBS和CMC的高反镜中都有被观察到, 并且随着反射率的增大, 双折射效应会越发明显[18,19].
在CEAS中, 双折射效应使透射腔模在两个偏振方向出现分裂的现象, CEAS信号由于干涉发生形变. 而对于CRDS, 腔镜的双折射效应以及腔后光学器件(包括透镜、探测器窗口[20]等)轻微的偏振效应, 会使得腔衰荡信号不再按e指数形式衰减, 引入衰荡时间的误差, 导致光谱系统准确性和探测灵敏度降低. 这限制了CEAS和CRDS在痕量气体检测中的应用. 当腔镜还同时存在二向色性时, 对于不同偏振态反射率不同, 将进一步恶化系统的性能. 因此分析腔镜的双折射效应及其影响对于提升CEAS和CRDS的可靠性和灵敏度都有重要意义.
Huang等[21]将FP腔的两个腔镜看成两个单独的双折射器件, 利用琼斯矩阵以及光腔内往返光场公式, 建立了完备的FP腔双折射模型, 并且分析了其对衰荡信号的影响; 2016年, Fleisher等[22]使用更加精简的模型, 将光学腔看成一个整体, 将腔内模式分为两个正交的本征偏振方向进行分析, 把两个本征频率差作为拟合参量, 最后针对腔衰荡光谱建立了非e指数模型(NEM), 测量结果显示该模型可以较好地反映双折射影响下的腔衰荡信号. 然而, 他们没有考虑两个偏振态耦合效率不同对腔衰荡信号的影响, 并且目前还没有相关工作对双折射效应对CEAS的影响进行建模和评估.
本文基于一个窄线宽光纤激光器以及基于IBS腔镜搭建的FP腔, 首先测量了双折射效应下的CEAS信号, 即透射腔模信号, 并建立模型对其进行了拟合; 通过透射腔模分裂现象并结合调制边带获得了两个偏振态的本征频率差; 同时测量并且分析了本征频率差与腔内气体压强的关系; 之后基于NEM模型, 考虑了耦合效率对于衰荡信号的影响; 最后对衰荡信号进行准确拟合, 研究双折射效应对衰荡时间的影响.
2. 实验装置
研究FP腔双折射效应及其对CEAS和CRDS影响的实验装置如图1所示, 激光器(Koheras Adjustik E15)输出波长为1531 nm的单纵模线偏光, 激光线宽为1 kHz. 输出光先后经过一个光纤电光调制器(EOM)和光纤声光调制器(AOM), 前者用于给腔的两个本征偏振方向的频率差进行定标, 后者作为光学开关用来快速斩断激光获得腔衰荡信号. 之后激光又经过一个匹配透镜(ML)、一个二分之波片(λ/2)和一个偏振分束器(PBS), 其中λ/2和PBS配合用于提高输出线偏光的偏振度以及调节腔前光强, 实测PBS的消光比达1∶20000. 出射的线偏光通过一个四分之波片(λ/4)变为椭圆偏振光后射入一个FP腔, 通过旋转λ/4可以改变光投影在FP腔两个本征偏振方向的比例.
FP腔由两个高反镜组成, 高反镜采用IBS镀膜技术, 腔的精细度F大约为3000, 腔长为39.4 cm, 对应自由光谱区为380 MHz, 腔纵模线宽大约为120 kHz. FP腔的一个本征偏振方向与PBS的透射偏振方向一致. 腔体采用具有低温膨胀系数的殷钢材料. 腔的透射信号射入探测器(PD1), PD1输出电压分为两路: 一路通过阈值电路, 用于判断腔模幅度继而输出脉冲信号控制AOM关断, 以触发腔衰荡事件; 另一路输出通过数据采集卡进入电脑. 当进行CRDS结果分析时, 为了更清楚说明双折射对CRDS的影响, 并且保持与前人工作的一致性, 会在透射探测器前放一个偏振片(PA).
3. 双折射效应对CEAS的影响
3.1 双折射频率分裂间隔
通过旋转λ/4波片使得其偏振方向与PBS出射方向成45°夹角, 此时腔入射光为圆偏光, 其在腔的两个本征偏振方向的分量相等. 由于腔镜双折射效应, 激光在两个偏振方向的相移不再相同, 表现为腔模信号在两个本征偏振方向的分裂. 实验中用一个频率为10 Hz的三角波扫描腔长, 观测到的透射腔模如图2中黑线所示, 曲线中心处可以清楚观测到分裂的腔模信号, 两个峰分别对应于两个本征偏振方向的模式, 并且两个模式的幅度接近. 另外, 为了对腔模频率分裂间隔进行频率定标, 我们使用一个频率为1 MHz的射频信号对EOM进行调制, 调制系数为0.3, 此时在透射腔模上可以观测到两个调制边带, 同时也可以清楚看到边带模式的分裂. 利用边带到中心腔模的距离和中心腔模两个峰值距离的比值可以定标出频率分裂间隔为98.9 kHz.
3.2 双折射频率间隔与腔内压力的关系
在腔衰荡光谱和腔增强光谱技术中, 为了减少腔内介质的干扰, 一般会进行抽真空处理. 而镀膜表面压力变化会导致腔镜发生形变, 诱导镀膜内部应力发生变化, 引起折射率分布的各向异性, 从而引起双折射效应的变化.
为了研究腔内气压与双折射频率分裂的关系, 向腔内冲入不同气压的氮气, 利用上述频率定标方法测量了不同气压下的频率间隔. 气压范围从0.0014—8.83 Torr (1 Torr = 133.322 Pa), 跨越了3个数量级. 频率间隔随气压的变化结果如图3所示, 在低气压时, 由于腔镜存在较大形变, 双折射引起的频率间隔较大, 超过104 kHz, 而随着气压升高达到低应力条件, 频率间隔会逐渐降低. 采用幂函数
f=a⋅xb 拟合, 结果如图3红线, 拟合参量a=100.97 ,b=−5.56×10−3 .3.3 透射腔模拟合
基于光学腔传递函数构建双折射透射腔模模型. FP腔的透射函数
T(ν) 可以用洛伦兹函数表示为ItI0=T(ν)=11+(2Fπ )2sin2(π ΔνFSR), (1) 式中
Δν 为激光频率失谐量, FSR表示腔的自由光谱区. 当存在双折射效应的情况下, 腔模存在分裂, 并且频率分裂间隔为δν , 并加入两个偏振方向幅度分量P1和P2, 透射函数可以表示为两个正交分量的叠加:ItI0=T(ν)=P11+(2Fπ )2sin2(π ΔνFSR)+P21+(2Fπ )2sin2(π (Δν−δν)FSR). (2) 用(2)式作为模型拟合透射腔模可以简单直接获得更精确的
Px(x=1,2) 值, 即各偏振分量在透射端的比例.为了验证理论模型的准确性, 测量了不同偏振下的CEAS信号, 并对其进行拟合. 通过调节腔前λ/4波片, 改变其与PBS偏振方向的夹角
θ , 可以改变腔前入射光的两个偏振态的幅度比例(即椭圆度). 选择了4个不同的值θ , 分别为0°, 30°, 45°和60°. 测量得到的透射腔模如图4黑色点所示, 采用(2)式进行拟合, 得到结果为红线所示. 拟合残差用图中下面板的黑线所示, 图中同时标明了拟合得到的P1和P2值. 从图4可以看出, 当θ 为0°时, 腔入射光接近线偏光, 并且方向与腔的其中一个本征方向重合, 因此腔透射光仍然为线偏光; 当θ 为30°时, 腔入射光为椭圆偏振光, 在FP腔的两个本征偏振方向都存在分量, 因此透射腔模出现分裂的现象, 并且两个峰值幅度不等, 得到P1 = 0.79, P2 = 0.21; 当θ 为45°时, 腔入射光在腔的两个本征偏振方向分量近似相等, 因此在透射端可以观测到等高的腔模信号; 当θ 为60°时, 腔入射光再次成为椭圆光, 透射腔模两个峰的幅度也再次不相等, P1 = 0.38, P2 = 0.62. 同时从拟合结果可以看出, 本工作给出的CEAS模型可以很好地匹配实验结果, 得到的拟合信噪比大于60.4. 双折射效应对CRDS信号的影响
4.1 双折射效应下的腔衰荡模型
基于NEM模型[22], 将腔透射光电场写成两偏振分量的叠加:
E(t)=E1(t)exp[−i(ω1t−ϕ1)]+E2(t)exp[−i(ω2t−ϕ2)], (3) 其中,
Ex(t) (x = 1, 2)为光在腔的两个本征偏振方向分量的电场幅度, 本征角频率ωx=2π υx ,ϕx 为相位. E1和E2的引入不仅考虑了光强在两个偏振方向的不同, 同时也包含了在两个方向激光到腔耦合效率的不同. 这是因为激光存在频率噪声, 并且光学腔由于压电陶瓷的存在, 在温度漂移和振动噪声等的影响下, 造成腔模频率不稳定, 从而导致激光到腔的耦合效率在两个偏振方向不同, 并且随着时间变化. 当电场经过线性PA, 且PA偏振方向与腔的其中一个本征偏振方向夹角为γ 时, 两偏振分量在PA轴上的投影可以写成:E(γ,t)=E1(t)exp[−i(ω1t−ϕ1)]cos(γ)+E2(t)exp[−i(ω2t−ϕ2)]sin(γ). (4) 因此光电探测器位置处的光强可以表示为
I(γ,t)=E(γ,t)E∗(γ,t)=I1(t)cos2(γ)+I2(t)sin2(γ)+√I1(t)I2(t)sin(2γ)cos(2\pi δνt+δϕ), (5) 其中,
2\pi δν 是双折射引起的频率分裂的角频率,δϕ 是相位差. 引入两偏振分量的比值p=I2/I1 来反映两个偏振方向光强的不同. 结合光强e指数衰减模型, 在考虑两个偏振方向光强不同的情况下, 腔衰荡光强信号可以写成:I(γ,t)=B+I1exp(−t/τ)[cos2(γ)+psin2(γ)+√psin(2γ)cos(2\pi δνt+δϕ)], (6) 其中, B表示偏置, 包括光电探测器暗电流、杂散光等影响.
4.2 实际拟合腔衰荡信号
实验测量的腔衰荡信号如图5中黑色点所示. 同时使用传统CRDS的e指数模型和(6)式所示的NEM模型对测量信号进行拟合, 得到的结果分别如图5(a), (b)中红线所示. 图5(c), (d)分别为拟合残差. 其中在NEM模型拟合中, 使用了通过拟合CEAS信号获取的
δν 值(98.9 kHz)以及p值(0.59),γ 值通过PA与λ/4的夹角获得为60°, 从而进一步简化了拟合, 并且提升了拟合的准确度. 从图5(a)可以清楚地看到测量的衰荡信号不再按照单e指数衰减, 在图5(c)中可以看到残差存在很大的结构性噪声. 而从图5(b)可以看出, 使用 (6) 式可以很好地拟合测量结果, 拟合残差中的结构性噪声被大大抑制, 拟合残差的标准偏差抑制了5倍, 从2.4×10–3减小到4.6×10–4. 根据两个模型得到的衰荡时间τ 分别为1.39 µs和1.37 µs, 偏差大于1.5%, 证实了双折射效应对CRDS进行痕量气体检测产生了很大的影响.为了验证(6)式对不同偏振光衰荡信号的适用性, 并进一步分析双折射效应对腔衰荡的影响, 以20°为步长旋转PA的角度, 范围为360°, 测量得到一系列腔衰荡信号, 分别使用e指数模型和NEM两种模型拟合, 得到衰荡时间
τ 如图6(a)所示. 图中黑色点线是单e指数拟合的结果, 可以看出衰荡时间起伏较大, 尤其在PA为80°和280°附近有突起, 衰荡时间最大值达到2.63 µs, 是平均值的1.9倍, 说明此处衰荡信号受双折射效应影响导致的衰荡线型扭曲最为严重. 红色点线为使用NEM模型拟合得到的结果, 衰荡时间基本保持不变, 极差小于0.16 µs. 图6(b)中显示了两种模型拟合衰荡信号获得的残差的标准偏差. 黑色表示的单e指数模型拟合残差的标准偏差与 图6(a)中趋势一致, 且整体始终大于5.8×10−4 , 在上述两突起处拟合残差较大, 最大可以达到3.7×10−3 . 红色表示的NEM模型拟合残差的标准偏差整体平稳, 均小于4.3×10−4 . 两种模型的拟合效果在PA为80°时差距最为明显, 拟合残差的标准偏差最大抑制了9倍. 因此可得出结论, 使用NEM模型可以更准确地描述腔衰荡信号.5. 总 结
本文研究了FP腔腔镜的双折射效应及其对腔增强吸收光谱和腔衰荡光谱技术的影响. 直接观测到了双折射效应下腔模透射信号的分裂现象, 通过使用EOM增加调制边带的方式对腔模分裂频率间隔进行了定标, 并且分析了腔内气体压强对频率分裂间隔的影响. 给出了双折射下腔增强信号模型, 并对实际信号进行了拟合, 验证其可以准确拟合出透射光中不同偏振成分的比例. 之后建立了用于分析腔衰荡信号的NEM模型, 考虑了耦合效率不同下的腔衰荡信号模型, 并对实际测量信号进行了准确拟合, 相较于传统e指数模型, 拟合残差的标准偏差最大抑制了9倍. 最后使用不同偏振下腔衰荡信号对模型进了验证, 证实了新模型对双折射下衰荡信号拟合的适用性.
[1] Strekalov D V, Thompson R J, Baumgartel L M, Grudinin I S, Yu N 2011 Opt. Express 19 14495
Google Scholar
[2] Kessler T, Hagemann C, Grebing C, Legero T, Sterr U, Riehle F, Martin M J, Chen L, Ye J 2012 Nature. Photon. 6 687
Google Scholar
[3] Jordan B C, William K, Martin M F, Eric G 2001 Appl. Opt. 40 3753
Google Scholar
[4] Ishibashi C, Sasada H 1999 Jpn. J. Appl. Phys. 38 920
Google Scholar
[5] Robert C 2007 Appl. Opt. 46 5408
Google Scholar
[6] Herriott D, Kogelnik H, Kompfner R 1964 Appl. Opt. 3 523
Google Scholar
[7] Liu J, Zhou Y, Guo S, Hou J, Zhao G, Ma W, Wu Y, Dong L, Zhang L, Yin W, Xiao L, Axner O, Jia S 2019 Opt Express. 27 1249
Google Scholar
[8] Livio G, Richard F W, Leo H 1999 J. Opt. Soc. Am. B 16 2247
Google Scholar
[9] 赵卫雄, 高晓明, 张为俊, 黄 腾 2006 光学学报 26 1260
Google Scholar
Zhao W X, Gao X M, Zhang W J, Huang T 2006 Acta Opt. Sin. 26 1260
Google Scholar
[10] 董美丽, 赵卫雄, 程跃, 胡长进, 顾学军, 张为俊 2012 物理学报 61 060702
Google Scholar
Dong M L, Zhao W X, Cheng Y, Hu C J, Gu X J, Zhang W J 2012 Acta. Phys. Sin. 61 060702
Google Scholar
[11] Zalicki P, Zare R N 1995 J. Chem. Phys. 102 2708
Google Scholar
[12] Zhao G, Bailey D M, Fleisher A J, Hodges J T, Lehmann K K 2020 Phys. Rev. A 101 062509
Google Scholar
[13] 胡纯栋, 焉镜洋, 王 艳, 梁立振 2018 光谱学与光谱分析 38 346
Hu C D, Y J Y, W Y, Liang L Z 2018 Spectrosc. Spect. Anal. 38 346
[14] Li Z X, Ma W G, Fu X F, Tan W, Zhao G, Dong L, Zhang L, Yin W B, Jia S T 2013 Appl. Phys. Express 6 072402
Google Scholar
[15] Ye J, Ma L S, Hall J 1996 Opt. Lett. 21 1000
Google Scholar
[16] Zhao G, Hausmaninger T, Ma W, Axner O 2018 Opt. Lett. 43 715
Google Scholar
[17] 肖石磊, 李斌成 2020 光电工程 47 190068
Xiao S L, Li B C 2020 Opto-Electronic Engineering 47 190068
[18] Winkler G, Perner L W, Truong G W, Zhao G, Bachmann D, Mayer A S, Fellinger J, Follman D, Heu P, Deutsch C, Bailey D M, Peelaers H, Puchegger S, Fleisher A J, Cole G D, Heckl O H 2021 Optica 8 686
Google Scholar
[19] Xiao S, Li B, Wang J 2020 Appl. Opt. 59 A99
Google Scholar
[20] 付小芳, 赵 刚, 马维光, 谭巍, 李志新, 董 磊, 张雷, 尹王保, 贾锁堂 2014 光谱学与光谱分析 34 1456
Google Scholar
Fu X F, Zhao G, Ma W G, Tan W, Li Z X, Dong L, Zhang L, Yin W B, Jia S T 2014 Spectrosc. Spect. Anal. 34 1456
Google Scholar
[21] Huang H F, Lehmann K K 2008 Appl. Opt. 47 3817
Google Scholar
[22] Fleisher A J, Long D A, Liu Q N, Hodges J T 2016 Phys. Rev. A 93 013833
Google Scholar
-
-
[1] Strekalov D V, Thompson R J, Baumgartel L M, Grudinin I S, Yu N 2011 Opt. Express 19 14495
Google Scholar
[2] Kessler T, Hagemann C, Grebing C, Legero T, Sterr U, Riehle F, Martin M J, Chen L, Ye J 2012 Nature. Photon. 6 687
Google Scholar
[3] Jordan B C, William K, Martin M F, Eric G 2001 Appl. Opt. 40 3753
Google Scholar
[4] Ishibashi C, Sasada H 1999 Jpn. J. Appl. Phys. 38 920
Google Scholar
[5] Robert C 2007 Appl. Opt. 46 5408
Google Scholar
[6] Herriott D, Kogelnik H, Kompfner R 1964 Appl. Opt. 3 523
Google Scholar
[7] Liu J, Zhou Y, Guo S, Hou J, Zhao G, Ma W, Wu Y, Dong L, Zhang L, Yin W, Xiao L, Axner O, Jia S 2019 Opt Express. 27 1249
Google Scholar
[8] Livio G, Richard F W, Leo H 1999 J. Opt. Soc. Am. B 16 2247
Google Scholar
[9] 赵卫雄, 高晓明, 张为俊, 黄 腾 2006 光学学报 26 1260
Google Scholar
Zhao W X, Gao X M, Zhang W J, Huang T 2006 Acta Opt. Sin. 26 1260
Google Scholar
[10] 董美丽, 赵卫雄, 程跃, 胡长进, 顾学军, 张为俊 2012 物理学报 61 060702
Google Scholar
Dong M L, Zhao W X, Cheng Y, Hu C J, Gu X J, Zhang W J 2012 Acta. Phys. Sin. 61 060702
Google Scholar
[11] Zalicki P, Zare R N 1995 J. Chem. Phys. 102 2708
Google Scholar
[12] Zhao G, Bailey D M, Fleisher A J, Hodges J T, Lehmann K K 2020 Phys. Rev. A 101 062509
Google Scholar
[13] 胡纯栋, 焉镜洋, 王 艳, 梁立振 2018 光谱学与光谱分析 38 346
Hu C D, Y J Y, W Y, Liang L Z 2018 Spectrosc. Spect. Anal. 38 346
[14] Li Z X, Ma W G, Fu X F, Tan W, Zhao G, Dong L, Zhang L, Yin W B, Jia S T 2013 Appl. Phys. Express 6 072402
Google Scholar
[15] Ye J, Ma L S, Hall J 1996 Opt. Lett. 21 1000
Google Scholar
[16] Zhao G, Hausmaninger T, Ma W, Axner O 2018 Opt. Lett. 43 715
Google Scholar
[17] 肖石磊, 李斌成 2020 光电工程 47 190068
Xiao S L, Li B C 2020 Opto-Electronic Engineering 47 190068
[18] Winkler G, Perner L W, Truong G W, Zhao G, Bachmann D, Mayer A S, Fellinger J, Follman D, Heu P, Deutsch C, Bailey D M, Peelaers H, Puchegger S, Fleisher A J, Cole G D, Heckl O H 2021 Optica 8 686
Google Scholar
[19] Xiao S, Li B, Wang J 2020 Appl. Opt. 59 A99
Google Scholar
[20] 付小芳, 赵 刚, 马维光, 谭巍, 李志新, 董 磊, 张雷, 尹王保, 贾锁堂 2014 光谱学与光谱分析 34 1456
Google Scholar
Fu X F, Zhao G, Ma W G, Tan W, Li Z X, Dong L, Zhang L, Yin W B, Jia S T 2014 Spectrosc. Spect. Anal. 34 1456
Google Scholar
[21] Huang H F, Lehmann K K 2008 Appl. Opt. 47 3817
Google Scholar
[22] Fleisher A J, Long D A, Liu Q N, Hodges J T 2016 Phys. Rev. A 93 013833
Google Scholar
计量
- 文章访问数: 6004
- PDF下载量: 119