光学学报, 2020, 40 (6): 0601002, 网络出版: 2020-03-06   

拉萨光学湍流探空测量与模式分析 下载: 1036次

Optical Turbulence Measurement over Lhasa Using a Meteorological Radiosonde and Model Analysis
作者单位
1 中国科学院安徽光学精密机械研究所中国科学院大气光学重点实验室, 安徽 合肥 230031
2 中国科学技术大学研究生院科学岛分院, 安徽 合肥 230031
摘要
利用自行研制的湍流气象探空仪,对西藏拉萨的温度、风速、风向等常规气象参数廓线和 Cn2廓线进行探测。分析了湍流强度和常规气象参数随高度的变化趋势。比较了夜晚和早晨的湍流特征,发现8~15 km处两者均出现强湍流层,且8 km以上早晨的湍流强度大于夜晚的湍流强度。同时基于Hufnagel-Vally 5/7模式,利用探空数据拟合得到符合拉萨湍流特征的拉萨 Cn2经验模式。对此经验模式进行统计分析,结果进一步证明拉萨 Cn2经验模式能有效估算拉萨的湍流强度。最后,将拉萨与高美谷的探空数据进行对比分析,发现拉萨的风速较小对天文观测有利,但湍流强度较强对天文观测有一定影响。该研究为后续拉萨湍流廓线和天文台选址的研究提供了参考,并为光电工程的应用提供了技术支持。
Abstract
This paper describes the use of a self-made meteorological radiosonde for turbulence analysis to measure profiles of the turbulence parameter Cn2 and meteorological parameters such as temperature, wind speed, and wind direction in Lhasa, Tibet. The trends of variation of turbulence intensity and meteorological parameters with height are analyzed. A comparison of the turbulence characteristics in the morning and evening reveals that a strong turbulence layer appears at 8-15 km at both these durations, but the turbulence intensity above 8 km in the morning is greater than that in the evening. Furthermore, based on the Hufnagel-Vally 5/7 model, the Cn2 Lhasa model is fitted using the statistical average of all radiosonde data. The statistical analysis results show that the Cn2 Lhasa model can effectively estimate the turbulence intensity above Lhasa. Finally, a comparison of the radiosonde data of Lhasa and Gaomeigu reveals that the wind speed in Lhasa is lower and favorable for astronomical observations, but the turbulence in Lhasa is stronger and reduces the quality of astronomical observations. These results provide a foundation for the study of turbulence profiles, observatory site selection, and support for the application of electro-optical engineering in Lhasa.

1 引言

由温度改变引起的折射率的起伏,被称为大气光学湍流现象。表征大气光学湍流强度的参数为大气折射率结构常数( Cn2),其随高度的分布称为大气光学湍流廓线。大气光学湍流不仅会影响天文望远镜的观测水平,还会干扰激光传输,降低光学遥感成像的质量,影响光通信信号的稳定性等[1]。获得大气光学湍流廓线,对减弱大气光学湍流对各种仪器设备影响的选取和光学系统性能的改善至关重要。

为了避免光学湍流对天文台选址的影响,全球的天文学家都在寻找适合天文观测的地址。中国西部得天独厚的地理优势和气候优势有利于天文观测,吸引了众多天文学家的关注。目前,对中国西部湍流特性的研究已经在阿里[2-4]、那区[5]、昌都[6]等地区开展,但对拉萨湍流参数的研究仍不够深入。拉萨是西藏的政治、经济、文化和科教中心,位于青藏高原南部(91.13°E, 29.67°N),海拔约3650 m,属于高原山地气候。该地区日照时间充足,晴夜多,降水少,空气干燥。早在1989年,Wang等[7]就曾指出拉萨可作为天文观测的候选地址进行进一步观测。王红帅[8]曾利用WRF(Weather Research and Forecasting)模型计算中国地区的光学湍流参数,结果表明青藏高原南部的视宁度、大气相干时间和等晕角较好,对天文观测有利。张永婧等[9]曾根据地面风速、云量、水汽等气象参数,将中国分为10个天文候选区,结果表明青藏高原南部天文气候条件较好。

目前用于 Cn2廓线测量的设备众多,例如探空气球[10]、MASS (Multi Aperture Scintillation Sensor)[11]、SCIDAR (Scintillation Detection and Ranging)[12]、SNODAR (Surface layer Non-Doppler Acoustic Radar)、Scintillometer[13]等。其中,湍流气象探空仪作为直接测量湍流的典型手段,可以获得 Cn2廓线和常规气象参数廓线,具有高时空分辨率的特点。然而,到目前为止还没有在拉萨开展过湍流探空实验。因此,有必要对拉萨进行探空测量,获取该地区的大气湍流廓线和相关湍流参数,进一步探究该地区的湍流强度垂直分布特点。

中国科学院合肥物质科学研究院于2018年8月利用自制的湍流气象探空仪首次对西藏拉萨的常规气象参数(温度、气压、湿度、风速、风向)廓线和 Cn2廓线进行测量研究。本研究利用Hufnagel-Vally 5/7 (简写为HV 5/7)模式和AFGL AMOS 模式获得 Cn2廓线,并与实测廓线进行对比;分析了常规气象参数、 Cn2随高度的分布特征及可能的原因。对比早晨和夜晚湍流特征,利用实测的 Cn2廓线以HV 5/7模式为基础拟合得到拉萨 Cn2经验模式;最后与高美谷的 Cn2廓线、温度廓线进行了对比分析。

2 实验概况

2.1 实验介绍

2018年8月3日至8月18日,中国科学院合肥物质科学研究院利用自行研制的湍流气象探空仪在青藏高原南部的西藏拉萨气象局开展了探空实验。该探空仪主要是在GPS 电子气象探空仪的基础上附带微温脉动仪[14]。微温脉动仪的统计平均时间为5 s,频率响应范围是0.1~30 Hz,温差标准差不大于0.002 ℃,等效 Cn2噪声为 2×10-18 m-2/3。此次共施放15个热力湍流探空气球。每次施放探空气球的时间通常在早上7:30和晚上7:30。探空气球的上升速度约为5 m/s,垂直分辨率约30 m。由于信号传输和天气干扰等其他原因,本研究选取了7份探空数据,其中2份早晨数据,5份夜晚数据。每次探空结束,即可获得 Cn2、风速、气压、温度、相对湿度等参数的廓线分布。具体探空记录如表1所示。

表 1. 拉萨探空记录

Table 1. Record of balloon soundings in Lhasa used in the paper

Balloon numberLaunch dateLaunch time (local time)Termination time (local time)Termination altitude (above sea level) /m
1#2018.08.1407:3808:5931023.3
2#2018.08.1507:2408:5130217.7
3#2018.08.1419:4021:0530657.1
4#2018.08.1422:5223:5828594.2
5#2018.08.1619:1920:4331206.4
6#2018.08.0819:1620:3826900.7
7#2018.08.1319:2320:2923221.1

查看所有表

2.2 Cn2测量原理

对于充分发展的湍流,假设其满足Kolmogorov 局地均匀各向同性的湍流理论,则温度结构常数 CT2与空间两点(距离为r)的温差满足

CT2=<[T(x)-T(x+r)]2>r-2/3,(1)

式中:l0rL0;T为大气温度,单位为K;T(x)和T(x+r)分别代表空间距离为r的两点温度;<·>代表系综平均;L0l0 分别为湍流外尺度、内尺度。

在可见光和近红外波段,忽略湿度的影响,温度变化是影响折射率变化的主要因素。因此近似认为折射率结构常数 Cn2可以直接由某一高度h处的气温T、气压PCT2计算获得[15-16],即

Cn2(h)=79×10-6P(h)T2h2CT2(h),(2)

式中:h的单位为m;P的单位为hPa; Cn2的单位为m-2/3。其中,温差的变化由温度脉动仪相距为r的一对微温探头获得,两者空间距离为1 m。本实验中,微温探头为铂丝,电阻约10 Ω,直径约20 μm。通过不同位置处的微温探头,测量空间两点的温差ΔT,并将其转换为电阻的变化。再通过由一对微温探头组成的惠斯通电桥,将电阻的变化与电压的变化ΔV联系起来,得到ΔV与ΔT的正比关系,即

ΔV=a·ΔT,(3)

式中:a 为标定系数。因此,由(3)式获得两点处的温差,由(1)式惯性区距离为r的两点温差的平方平均得到 CT2,最后由(2)式得到 Cn2

3 结果和分析

3.1 实测常规气象参数廓线及 Cn2廓线分析

利用上述所有探空数据,计算对应高度处的温度梯度和风速梯度。为了定量描述,本研究每隔100 m对温度、气压、相对湿度、风速、风向、温度梯度、风速梯度、 CT2Cn2进行统计平均,得到对应廓线,如图1所示。为了便于描述,本研究中若没有特殊说明,则所有高度均指相对地面高度(above ground level)。由图可知,对流层顶出现在14 km左右。近地面层内相对湿度随高度升高,3 km处达到最大值,约为80%;随后逐渐降低。风速随高度增大,24 km左右到达最大值,近似为26 m/s。风向在近地面到15 km范围内波动较大,15 km以上趋于稳定。在5 km左右和对流层15 km左右温度梯度均出现剧烈变化,均对应图1(i)中的强湍流层。风速梯度整体趋势变化不大。 CT2在近地面较大,近地面至0.9 km随高度迅速减小,此后随高度逐渐变大。 Cn2在近地面最大,近地面至0.9 km随高度迅速递减;在对流层顶14 km附近存在强湍流层。

3.2 Cn2廓线实测值与模式估算值对比分析

表 2. 拉萨探空测量的lg(C

Table 2. Statistical analysis for lg(
ModelMRERMSEBIASRxy
Cn2 Lhasa model0.60%0.13-0.010.97

查看所有表

国内外学者对不同实验地点进行长期的测量和统计分析,提出了不同的 Cn2廓线模式,包括AFGL AMOS(Air Fore Geophysics Laboratory)模式和Hufnagel 模式。本研究利用这两种模式,对拉萨的湍流强度进行了研究。

AFGL AMOS模式是由Good等[17]和Beland等[18]基于十多年的探空气球数据获得的 Cn2廓线,而提出的 Cn2的经验公式。通过计算上述数据的算术平均,并进行多项式分段拟合,最终获得分段表达式。该探空气球施放于Mt. Haleakala 上的 AMOS观测站,具有高时空分辨率。AFGL AMOS 白天和夜间模式的具体表达式分别为

lg(Cn2)=0.0482-2.3416h-0.7211h2,3.052km<h3.54km-17.4778+0.0320h+0.0078h2+1.5066exp-0.5h-4.46030.09682,3.54km<h5.1km-16.5589-0.1424h+0.0030h2+1.5694exp-0.5h-16.63004.87572,h>5.1km,(4)lg(Cn2)=-12.412-0.4713h-0.0906h2,3.052km<h5.2km-17.1273-0.0301h-0.0010h2+0.5061exp-0.5h-15.08663.29772,h>5.2km,(5)

式中:h为海拔高度,单位km。

图 1. 实测探空数据平均廓线。(a)温度;(b)气压;(c)相对湿度;(d)风速;(e)风向;(f)温度梯度;(g)风速梯度;(h)温度结构常数;(i)折射率结构常数

Fig. 1. Average profiles of measured data by sounding balloons. (a) Temperature; (b) pressure; (c) relative humidity; (d) wind speed; (e) wind direction; (f) temperature gradient; (g) wind speed gradient; (h) CT2

下载图片 查看所有图片

Hufnagel模式最早是Hufnagel 和 Stanley[19]根据中纬度区域的气球测量和星光闪烁数据获得的,忽略了3 km 以下湍流的强度。Urich[20]根据Vally提出的建议,新增了3 km以下的湍流强度,故能对 Cn2的整层廓线进行描述,即

Cn2(h)=8.2×10-16W2(h/10)10exp(-h)+2.7×10-16exp(-h/1.5)+Aexp(-h/0.1),(6)W2=115520V2(h)dh,(7)

式中:W表5 ~20 km 高度处的方均根风速,单位是 m/s;V(h)为某一高度处的风速,单位是m/s;A是描述大气近地面光学湍流强度的参量。可调参数WA可以通过固定等晕角和相干长度来选取。当等晕角为7 μrad,相干长度为5 cm(λ=0.5 μm)时,对应的A=1.7×10-14 m-2/3,W =21 m/s,即为常用的HV 5/7模式。

图2展示了拉萨早晨和夜晚的平均湍流廓线,以及对应的AFGL AMOS和HV 5/7模式估算值。从图中可以看出,实测 Cn2廓线、AFGL AMOS夜间模式和HV 5/7模式估算值有相同的特点:近地面湍流强度通常最强, Cn2在近地面层随高度急剧减小; Cn2在平流层随高度逐渐减小;对流层均包含强湍流层,但强湍流层的位置有所不同。值得指出的是,AFGL AMOS白天模式在近地面层出现强湍流薄层,此时湍流强度最大。

图 2. 实测早晨和夜晚Cn2廓线平均值与模式值。(a)早晨;(b)夜晚

Fig. 2. Measured mean Cn2 profiles in the morning and evening and model value. (a) In the morning; (b) in the evening

下载图片 查看所有图片

3.3 Cn2经验模式的建立及统计分析

3.3.1 Cn2经验模式的建立

目前用于测量大气光学湍流廓线的方法众多,但大多数测量方法都需要花费巨大的人力、物力、财力。因此有必要发展大气光学湍流廓线模式来估算大气光学湍流强度。基于HV 5/7模式,对拉萨所有 Cn2廓线进行统计平均,并拟合得到符合该地区湍流强度变化特征的统计平均模式,在此称为拉萨 Cn2经验模式( Cn2Lhasa model),具体表达式为

Cn2=1.85×10-19h5.15exp(-h/1.82)+3.55×10-17exp(-h/7.47)+6.6×10-16exp(-h/0.0018)(8)

图3为早晨、夜晚、早晚平均的 Cn2廓线以及拟合得到的 Cn2经验模式廓线。由图中可以看出,早晚的 Cn2廓线从近地面到8 km处较为一致;8~15 km范围均出现强湍流层,但8 km以上夜晚的湍流强度小于早晨的湍流强度。与早晨相比,夜晚湍流强度波动较小,大气较为稳定。拟合得到的拉萨 Cn2经验模式,不管是从整体变化趋势还是从量级上都能基本完成对拉萨实测 Cn2廓线的估算,进而反映该地区湍流强度的变化规律。

图 3. 实测早晨、夜晚和早晚平均Cn2廓线及

Fig. 3. Cn2 profiles measured in the morning/evening, or averaged and estimated using 下载图片 查看所有图片

3.3.2 统计结果

为了定量地描述拉萨 Cn2经验模式的估算效果,本研究使用4个统计量,分别为:平均相对误差 (MRE)、方均根误差 (RMSE)、平均偏差 (BIAS)以及相关系数 (Rxy)对lg( Cn2)进行统计分析。统计量的具体表达式为

EMR=i=1NXi-YiNi=1NXiN×100%,(9)ERMS=i=1N(Xi-Yi)2N,(10)EBIAS=i=1NXi-YiN=X--Y-,(11)Rxy=i=1N(Xi-X-)(Yi-Y-)i=1N(Xi-X-)2i=1N(Yi-Y-)2,(12)

式中:XiYi分别代表第i个高度上探空气球实测和 Cn2经验模式估算的lg( Cn2);N代表对应的(Xi,Yi)数据个数; X-Y-分别代表实测值和估算值的平均值。

统计结果如表2所示: Cn2经验模式的平均相对误差(EMR)在1%以下;方均根误差(ERMS)不超过0.20;平均偏差(EBIAS)为-0.01;相关系数(Rxy)为0.97。因此,误差分析证实了拉萨 Cn2经验模式的有效性。

4 丽江高美谷和西藏拉萨探空数据对比分析

高美谷观测站坐落于云南省丽江市(100.0°E,26.7°N),海拔约3227 m,具有高原山地和亚热带季风气候。本研究对丽江高美谷与西藏拉萨的探空数据分别进行统计平均,便于对比分析。表3为高美谷探空记录。

图4为拉萨和高美谷温度、风速和 Cn2的平均廓线对比。从图4(a)可以看出,拉萨和高美谷对流层顶分别出现在14 km和15 km左右。在距地面2 km左右,高美谷地区出现逆温现象,对应于图4(c)中薄的强湍流层。由图4(c)可知,两者在15 km附近均出现强湍流层。

图 4. 拉萨和高美谷探空实测的平均廓线。(a)温度廓线;(b)风速廓线;(c) Cn2廓线

Fig. 4. Average profiles measured by balloon soundings in Gaomeigu and Lhasa. (a) Temperature; (b) wind speed; (c) Cn2 profile

下载图片 查看所有图片

表 3. 丽江高美谷探空记录

Table 3. Record of balloon soundings in Gaomeigu used in the paper

Balloon numberLaunch dateLaunch time (local time)Termination time (local time)Termination altitude(above sea level) /m
1#2013.04.0611:1112:3630148
2#2013.04.0701:2804:4431421
3#2013.04.0801:3203:1628714
4#2013.04.1323:0600:5531678
5#2013.04.1405:1707:1432467
6#2013.04.2101:0002:5230473
7#2013.04.2321:0722:5431383

查看所有表

5 结论

利用自行研制的湍流气象探空仪,获得了拉萨的 Cn2廓线、 CT2廓线和常规气象参数廓线。实验发现,在10-20~10-16m-2/3范围内,该地区大气光学湍流强度在地表附近最强,在近地面层随高度急剧减弱,在对流层顶14 km附近出现强湍流层。强湍流层的出现似乎与温度梯度的大小以及对流层顶高度相对应,再次印证了光学湍流的产生与温度梯度有关的结论[21-22]。在大气稳定的状态下,对流层顶附近会出现较强的湍流层,这可能是冷暖空气相遇、强对流天气发生、重力波等造成的湍流强度变化[23]。对比早晨和夜晚的湍流廓线发现,高于8 km夜晚湍流强度逐渐小于早晨湍流强度,但两者强湍流层均出现在8~15 km。利用本次实验有限的探空数据,基于HV 5/7模式拟合得到拉萨 Cn2经验模式,该模式能在一定程度上重现实验结果,但不能有效地反映大气光学湍流复杂的细节特征,存在一定的局限性,需要大量的探空数据结果进行不断地修正。

通过拉萨和高美谷温度、风速、 Cn2廓线的对比可以看出,两者风速和 Cn2变化趋势有所不同。拉萨无明显的高空急流,高美谷在9 km左右存在高空急流。高美谷湍流强度小于拉萨。两地 Cn2廓线和风速的差异,这可能是由于地理位置和气候的不同造成的。与高美谷相比,拉萨较小的风速适合于天文观测,但较强的光学湍流却对天文观测不利。这也体现出天文台址的选取是一个长期而复杂的过程,应该综合考虑各方面的因素。

参考文献

[1] Hutt D L. Modeling and measurements of atmospheric optical turbulence over land[J]. Optical Engineering, 1999, 38(8): 1288-1295.

[2] Liu L Y, Yao Y Q, Vernin J, et al. Atmospheric turbulence measurements at Ali observatory, Tibet[J]. Proceedings of SPIE, 2012, 8444: 844464.

[3] Ye Q Z, Su M, Li H, et al. Tibet's Ali:asia's Atacama?[J]. Monthly Notices of the Royal Astronomical Society: Letters, 2016, 457(1): L1-L4.

[4] Qian X, Yao Y Q, Wang H S, et al. Thecharacteristics at the Ali observatory based on radiosonde observations[J]. Publications of the Astronomical Society of the Pacific, 2018, 130(994): 125002.

[5] Li M S, Ma Y M, Ma W Q, et al. Analysis of turbulence characteristics over the northern Tibetan Plateau area[J]. Advances in Atmospheric Sciences, 2006, 23(4): 579-585.

[6] Bian L G, Xu X D, Lu L H, et al. Analyses of turbulence parameters in the near-surface layer at Qamdo of the Southeastern Tibetan Plateau[J]. Advances in Atmospheric Sciences, 2003, 20(3): 369-378.

[7] Wang S G, Lan S Z, Jiang S Y, et al. Preliminary investigation on potential astronomical sites in Tibet[J]. Publications of Shanxi Observatory, 1989, 12(1/2): 1-17.

[8] 王红帅. 大气光学湍流模型研究与实践[D]. 北京: 中国科学院国家天文台, 2012: 69- 74.

    Wang HS. Modeling research of practice on atmospheric optical turbulence[D]. Beijing: National Astronomical Observatories,Chinese Academy of Sciences, 2012: 69- 74.

[9] 张永婧, 王盘兴, 姚永强, 等. 中国天文气象条件的地域分布[J]. 中国科学: 物理学力学天文学, 2010, 40(10): 1302-1314.

    Zhang Y J, Wang P X, Yao Y Q, et al. Climate regionalization in China for astronomical observations[J]. Science in China Series G-Physics, Mechanics & Astronomy (in Chinese), 2010, 40(10): 1302-1314.

[10] Marks R D, Vernin J, Azouit M, et al. Antarctic site testing-microthermal measurements of surface-layer seeing at the South Pole[J]. Astronomy and Astrophysics Supplement Series, 1996, 118(2): 385-390.

[11] Kornilov V, Tokovinin A, Shatsky N, et al. Combined MASS-DIMM instruments for atmospheric turbulence studies[J]. Monthly Notices of the Royal Astronomical Society, 2007, 382(3): 1268-1278.

[12] Avila R, Vernin J, Masciadri E. Whole atmospheric-turbulence profiling with generalized scidar[J]. Applied Optics, 1997, 36(30): 7898-7905.

[13] Vernin J, Roddier F. Experimental determination of two-dimensional spatiotemporal power spectra of stellar light scintillation evidence for a multilayer structure of the air turbulence in the upper troposphere[J]. Journal of the Optical Society of America, 1973, 63(3): 270-273.

[14] 蔡俊, 李学彬, 詹国伟, 等. 一个新的海边光学湍流外尺度和 Cn2的廓线模式[J]. 物理学报, 2018, 67(1): 014206.

    Cai J, Li X B, Zhan G W, et al. A new model for the profiles of optical turbulence outer scale and Cn2 on the coast[J]. Acta Physica Sinica, 2018, 67(1): 014206.

[15] Tatarskii VI. The effects of the turbulent atmosphere on wave propagation[M]. Jerusalem: Israel Program for Scientific Translations, 1971: 66- 80.

[16] Lawrence R S, Ochs G R, Clifford S F. Measurements of atmospheric turbulence relevant to optical propagation[J]. Journal of the Optical Society of America, 1970, 60(6): 826-830.

[17] Good R E, Beland R R, Murphy E A, et al. Atmospheric models of optical turbulence[J]. Proceedings of SPIE, 1988, 0928: 165-186.

[18] Beland RR, Brown JH, Good RE, et al., 1985: AFGL-TR-88-0153.

[19] Hufnagel R E, Stanley N R. Modulation transfer function associated with image transmission through turbulent media[J]. Journal of the Optical Society of America, 1964, 54(1): 52-61.

[20] Ulrich P B. Hufnagel-Valley profiles for specified values of the coherence length and isoplanatic angle[R]. WJ Schafer Associates, 1988: MA-TN-88-013.

[21] 青春, 吴晓庆, 李学彬, 等. WRF模式估算丽江高美古大气光学湍流廓线[J]. 中国激光, 2015, 42(9): 0913001.

    Qing C, Wu X Q, Li X B, et al. Estimation ofatmospheric optical turbulence profile by WRF model at Gaomeigu[J]. Chinese Journal of Lasers, 2015, 42(9): 0913001.

[22] 青春, 吴晓庆, 李学彬, 等. 典型地区高空大气光学湍流模拟研究[J]. 光学学报, 2016, 36(5): 0501001.

    Qing C, Wu X Q, Li X B, et al. Research on simulating atmospheric optical turbulence in typical area[J]. Acta Optica Sinica, 2016, 36(5): 0501001.

[23] 吴晓庆, 钱仙妹, 黄宏华, 等. 丽江高美古视宁度、等晕角及相干时间的探空测量[J]. 天文学报, 2014, 55(2): 144-153.

    Wu X Q, Qian X M, Huang H H, et al. Measurements of seeing, isoplanatic angle, and coherence time by using balloon-borne microthermal probes at Gaomeigu[J]. Acta Astronomica Sinica, 2014, 55(2): 144-153.

韩亚娟, 吴晓庆, 罗涛, 苏昶东, 青春, 吴骕, 杨期科, 张坤. 拉萨光学湍流探空测量与模式分析[J]. 光学学报, 2020, 40(6): 0601002. Yajuan Han, Xiaoqing Wu, Tao Luo, Changdong Su, Chun Qing, Su Wu, Qike Yang, Kun Zhang. Optical Turbulence Measurement over Lhasa Using a Meteorological Radiosonde and Model Analysis[J]. Acta Optica Sinica, 2020, 40(6): 0601002.

本文已被 4 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!