Measurement uncertainty of temperature pulsation meter based on instrument amplifier
-
摘要:
表征大气湍流强度的关键参数一般是折射率结构常数，其由温度脉动仪测量得到。为了研究该常数的测量不确定度，基于温度脉动仪的工作原理，采用直径为10 μm、电阻值为90 Ω、温度变化系数为4.5×10-3/K的钨丝作为测量传感器，获取温度差系综平均，并将同步监测到的压强、温度、距离等参数进行综合分析计算，得到折射率结构常数的测量不确定度。结果表明，当近地面层的大气折射率结构常数的真实值范围在10-13~10-16时，计算得到的折射率结构常数测量不确定度可控制在1%之内，具有较高的可信度。这一结果对于提高温度脉动仪的测量误差具有很大的作用，同时也为控制其测量不确定度提供了新的解决途径。
Abstract:The refractive index structure constant, as one of the key indicators characterizing atmospheric turbulence intensity, is measured by the temperature pulsation meter. Research on its measurement uncertainty is of significant importance. Based on the principle of the temperature pulsation meter, a tungsten filament with a diameter of 10 μm, a resistance value of 90 Ω and a temperature coefficient of variation of 4.5×10-3/K were used as the measuring sensor to obtain the average temperature different coefficient. The measurement uncertainty of refractive index structure constant was derived from an integrated analysis of parameters including pressure, temperature and distance. Experimental data indicates that within the true value range of 10-13 to 10-16 near the ground, the calculated measurement uncertainty remains within 1%, demonstrating its reliability. The result is very effective for improving the measurement resolution of temperature pulsation meter, and it is also a new idea for controlling the measurement uncertainty of temperature pulsation meter.
-
0. 引言
光波在大气中传输时，因大气湍流的存在，会对光波产生光斑的扩展、光强闪烁和光束的传播方向发生偏折等一系列现象[1-3]，从而限制了光电设备的性能。光波受湍流的影响主要是温度起伏引起的[4]，温度起伏引起折射率的起伏，造成光束波前的畸变。通过对温度起伏的直接测量，可以计算出大气折射率结构常数Cn2[5-6]。测量温度起伏的仪器为温度脉动仪[7-8]，通常利用10 μm左右的金属丝来测量单点的大气湍流强度，具有响应快速、线性度好、分辨率高等特点。将温度脉动仪放置在近地面的光电设备附近，可以测量近地面的Cn2[9-12]; 通过搭载温度脉动仪的测量塔，可以测量表面层的湍流强度随高度分布情况；将温度脉动仪搭载到探空气球上，可以测量光学站址的大气Cn2廓线[13-14]。
折射率结构常数Cn2作为一个随机量，其测量值的不确定度对于准确测量外场湍流强度、评估光电设备的性能具有重要意义。本文作者从温度脉动仪的测量原理出发，结合测量不确定度的定义，分析了影响温度脉动仪测量结果的各个因素，并定量计算了温度脉动仪的测量不确定度。
1. 温度脉动仪测量原理
1.1 测量原理
大气湍流强度的表征参数Cn2由科尔莫戈罗夫提出[15]，对局部均匀和各向同性的涡流，用折射率结构函数Dn(r)来表征折射率的变化，它与测量距离r存在正比关系，与r的2/3次方成正比：
Dn(r)=⟨[n(r1+r)−n(r1)]2⟩=C2n(h)×r2/3,(l≪r≪L) (1) 式中：l为湍流内尺寸; L为湍流外尺寸; n(r1)为r1处的折射率; n(r1+r)为r1+r处的折射率; Dn(r)是根据2个观测点间折射率增量，先计算均方值再统计平均值得到; Cn2(h)为h处的比例常数，是大气湍流强度的表示方式，表示h处的平均湍流强度，量纲为m－2/3，称为折射率结构常数。
折射率的变化是由于温度起伏而产生的，在湍流的惯性区间内，2个观测点的温度结构函数DT(r)只与2个观测点间的距离r2/3有关，与2个观测点的相对方向和位置没有关系[3]。如果大气中2个观测点之间的距离在湍流惯性区的范围内，则会满足以下表达式:
DT(r)=C2T(h)×r2/3,(l≪r≪L) (2) 式中：CT2(h)是温度结构常数，单位为K2·m－2/3; r的单位是m。
空间两点的温度结构函数为:
DT(r)=⟨[T(x)−T(x+r)]2⟩ (3) 式中: T(x)和T(x+r)分别为空间中两个不同位置点的温度; 〈[T(x)－T(x+r)]2〉表示空间中两个点温度差的系综平均。
由上面的两个式子可以得出大气的温度结构函数为:
C2T(r)=⟨[T(x)−T(x+r)]2⟩×r−2/3 (4) 在光学波段，Cn2和CT2的关系[10]可以描述为：
C2n=(79.2×10−6×pT2)2×C2T (5) 式中: p是空间测量点的气压(单位为Pa); T是绝对温度(单位为K)。温度脉动仪通过测量距离r、温度T、气压p和温差的系综平均值等数值，可以算出折射率结构常数。测量温度差系综平均值[16-17]的是金属丝传感器[18]，如铂丝、钨丝等。金属丝传感器的电阻值随温度变化呈线性变化。通过测量金属丝电阻值的变化，就可以推导出温度的变化。金属丝传感器的热惯性系数小，能够响应大气温度的微小和快速变化。
1.2 测量电路设计
温度脉动仪的电路示意图[19]如图 1a所示。利用惠斯通电桥，经过仪表放大器、滤波、采样和计算，测量距离为r的空间两点的金属丝传感器温度差的系综平均，同时测量压强和气温，就可由前面的公式计算大气折射率结构常数。图 1b为电路实物图。
![图 1 a—温度脉动仪的电路示意图b—电路实物图]()
图 1 a—温度脉动仪的电路示意图b—电路实物图Figure 1. a—circuit diagram of temperature pulsation meter b—circuit diagram2. 温度脉动仪测量不确定度分析
测量不确定度是依照所得到的数据信息，用来表征赋予被测量的量值分散性的非负参量。测量不确定度以定量的方式表示测量结果的可信度[20-22]。测量不确定度小，说明对应的特定测量获得的可能量值分布范围小，因而可信程度高，使用价值高。标准不确定度[23]通过用标准偏差的测量结果的不确定度来定义。就温度脉动仪而言，Cn2作为被测量的量，是一个随机量，因此无法通过独立重复测量的方法来估计测量不确定度，只能通过对影响Cn2的各输入量进行估计，得到合成标准不确定度。
根据式(4)，影响Cn2的输入量有压强p、距离r、温度T和两点间温度差的系综平均, 则:
U(Cn2)=√U2(p)+U2(T)+U2(r)+U2(⟨[T(x)−T(x+r)]2⟩) (6) 式中: U(p)表示压强引入的不确定度，是由压强传感器的测量不确定度给出; U(T)表示温度引入的不确定度，由温度传感器的测量不确定度给出; U(CT2)表示温度结构常数的不确定度，由温度脉动的测量不确定度和距离r的测量不确定度给出。
对于压强、温度和距离r，选用合适的传感器和测量工具，各输入变量的不确定度分别为:
U(p)=±1 Pa (7) U(T)=±0.3 K (8) U(r)=±0.01 m (9) 在近地面一般测量条件下，压强p=105 Pa，距离r=1 m，温度T=20 ℃(293 K)，不同的温度起伏对应的折射率结构常数见表 1。
表 1 近地面条件下不同温度起伏对应的折射率结构常数Table 1. Refractive index structure constants corresponding to different temperature fluctuations under near ground conditionsCn2/m－2/3 p/Pa T/K r/m 〈[T(x)－T(x+r)]2〉 8.5×10-19 1×105 293 1 10-6 8.5×10-18 1×105 293 1 10-5 8.5×10-17 1×105 293 1 10-4 8.5×10-16 1×105 293 1 10-3 8.5×10-15 1×105 293 1 10-2 8.5×10-14 1×105 293 1 10-1 8.5×10-13 1×105 293 1 100 8.5×10-12 1×105 293 1 101 8.5×10-11 1×105 293 1 102 两点温度差的系综平均值由被测量大气温度起伏和噪声共同决定。其中噪声因素包括金属丝传感器及其支架的热惯性、滤波电路的影响、采样电路中的噪声、金属丝中电流的加热效应、太阳直射对金属丝传感器的加热等[24]。金属丝及其支架的惯性对高频湍流有响应衰减，金属丝中电流的加热效应和太阳直射的加热是一个稳定量，可以通过电路的高通滤波过滤去除。滤波电路的影响对电路频率的影响可通过校正系数去除。采样电路中的噪声为选用芯片的噪声决定。
采用钨丝传感器，直径为10 μm，钨丝的温度变化系数为4.5×10-3/K, 电阻值为90 Ω。表 2中给出了两点之间不同的温度起伏，对应的电阻值变化和压强变化。
表 2 温度起伏对应的电阻值和电压变化Table 2. Changes of resistance and voltage corresponding to temperature fluctuations[T(x)-T(x+r)]/K resistance change value/Ω current value/mA voltage change value/mV 10-3 0.000405 1 0.000405 10-2 0.00405 1 0.00405 10-1 0.0405 1 0.0405 100 0.405 1 0.405 101 4.05 1 4.05 将钨丝传感器放置在密封环境中，作为没有温度变化情况，测量得到电路噪声，如图 2所示。经多次重复测量，电路电压起伏可以控制在0.00025 mV以下，对应的温度起伏为0.0003 K。图 3为一次典型的噪声测量结果，因此0.0003 K可以作为温度起伏的测量不确定度。
U(T(x)−T(x+r))=±0.0003 K (10) ![图 2 放置于封闭的纸箱内的温度脉动仪钨丝传感器]()
图 2 放置于封闭的纸箱内的温度脉动仪钨丝传感器Figure 2. Tungsten wire sensor of temperature pulsation meter in the closed cardboard box利式(5)对各参量求偏微分，可以得到各参量的不确定度引起的折射率结构常数的不确定度，如表 3所示。此时, p=105 Pa，T=293 K，r=1 m，温度起伏分别为10-3 K~100 K。可以看出，折射率结构常数的测量不确定度主要由电路的噪声决定。当温度起伏在10-3 K时，电路噪声几乎和温度起伏引起的电压起伏相当，此时相对不确定度可以达到60%。随着湍流强度的增强，相对不确定度也在逐步降低。
表 3 影响温度脉动仪测量不确定度的各参量对折射率结构常数不确定度的影响Table 3. Uncertainty of refractive index structure constant caused by various parameters affecting the measurement uncertainty of temperature pulsation meter[T(x)-T(x+r)]/K U(p) U(T) U(r) U(T(x)－T(x+r)) Cn2/m-2/3 U(Cn2) relative uncertainty/% 10-3 1.7×10-23 3.5×10-21 5.7×10-21 5.1×10-19 8.5×10-19 5.1×10-19 60.0 10-2 1.7×10-19 3.5×10-19 5.7×10-19 5.1×10-18 8.5×10-17 5.1×10-18 6.0 10-1 1.7×10-17 3.5×10-17 5.7×10-17 5.1×10-17 8.5×10-15 8.4×10-17 1.0 100 1.7×10-15 3.5×10-15 5.7×10-15 5.1×10-16 8.5×10-13 6.7×10-15 0.8 101 1.7×10-13 3.5×10-13 5.7×10-13 5.1×10-15 8.5×10-11 6.7×10-14 0.8 3. 结论
从温度脉动仪的测量原理出发，结合测量不确定度的定义，分析了影响温度脉动仪测量结果的各个因素，并定量计算了温度脉动仪的测量不确定度。研究表明，近地面层的大气折射率结构常数的范围通常在10－13 m-2/3~10－16 m-2/3之间[25]，此时温度脉动仪的测量不确定小于1%，有比较高的可信度。该研究对于提高温度脉动仪的测量分辨率有一定效果，同时也为控制其测量不确定度提供了新的解决途径。
-
![]()
图 1 a—温度脉动仪的电路示意图b—电路实物图
Figure 1. a—circuit diagram of temperature pulsation meter b—circuit diagram
![]()
图 2 放置于封闭的纸箱内的温度脉动仪钨丝传感器
Figure 2. Tungsten wire sensor of temperature pulsation meter in the closed cardboard box
表 1 近地面条件下不同温度起伏对应的折射率结构常数
Table 1 Refractive index structure constants corresponding to different temperature fluctuations under near ground conditions
Cn2/m－2/3 p/Pa T/K r/m 〈[T(x)－T(x+r)]2〉 8.5×10-19 1×105 293 1 10-6 8.5×10-18 1×105 293 1 10-5 8.5×10-17 1×105 293 1 10-4 8.5×10-16 1×105 293 1 10-3 8.5×10-15 1×105 293 1 10-2 8.5×10-14 1×105 293 1 10-1 8.5×10-13 1×105 293 1 100 8.5×10-12 1×105 293 1 101 8.5×10-11 1×105 293 1 102 
下载: 导出CSV
表 2 温度起伏对应的电阻值和电压变化
Table 2 Changes of resistance and voltage corresponding to temperature fluctuations
[T(x)-T(x+r)]/K resistance change value/Ω current value/mA voltage change value/mV 10-3 0.000405 1 0.000405 10-2 0.00405 1 0.00405 10-1 0.0405 1 0.0405 100 0.405 1 0.405 101 4.05 1 4.05
下载: 导出CSV
表 3 影响温度脉动仪测量不确定度的各参量对折射率结构常数不确定度的影响
Table 3 Uncertainty of refractive index structure constant caused by various parameters affecting the measurement uncertainty of temperature pulsation meter
[T(x)-T(x+r)]/K U(p) U(T) U(r) U(T(x)－T(x+r)) Cn2/m-2/3 U(Cn2) relative uncertainty/% 10-3 1.7×10-23 3.5×10-21 5.7×10-21 5.1×10-19 8.5×10-19 5.1×10-19 60.0 10-2 1.7×10-19 3.5×10-19 5.7×10-19 5.1×10-18 8.5×10-17 5.1×10-18 6.0 10-1 1.7×10-17 3.5×10-17 5.7×10-17 5.1×10-17 8.5×10-15 8.4×10-17 1.0 100 1.7×10-15 3.5×10-15 5.7×10-15 5.1×10-16 8.5×10-13 6.7×10-15 0.8 101 1.7×10-13 3.5×10-13 5.7×10-13 5.1×10-15 8.5×10-11 6.7×10-14 0.8
下载: 导出CSV
-
[1] 黄永平, 赵光普. 激光通信在大气湍流中的传输[J]. 大众科技, 2009(7): 39-41. HUANG Y P, ZHAO G P. The transmission of laser communication in atmospheric turbulence[J]. Volkswagen Technology, 2009 (7): 39-41(in Chinese).
[2] 孙刚, 翁宁泉, 肖黎明. 合肥地区大气折射率结构常数高度分布模式[J]. 强激光与粒子束, 2008, 20(2): 183-188. SUN G, WENG N Q, XIAO L M. Height distribution model of atmospheric refractive index structure constant in Hefei region[J]. High Power Laser and Particle Beams, 2008, 20(2): 183-188(in Chinese).
[3] 张志刚, 段哲民, 强希文, 等. 基于真均方根转换器的温度脉动仪设计[J]. 计算机测量与控制, 2012, 20(9): 2577-2580. ZHANG Zh G, DUAN Zh M, QIANG X W, et al. Design of temperature pulsation meter based on true root mean square converter[J]. Computer Measurement and Control, 2012, 20(9): 2577-2580(in Chinese).
[4] 姜楠, 李晓英, 牛春晖, 等. 大气湍流对激光空间传输特性影响的实验研究[J]. 激光技术, 2022, 46(5): 708-712. DOI: 10.7510/jgjs.issn.1001-3806.2022.05.022 JIANG N, LI X Y, NIU Ch H, et al. Experimental study on the influence of atmospheric turbulence on laser spatial transmission characteristics[J]. Laser Technology, 2022, 46(5): 708-712(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2022.05.022
[5] 翁宁泉, 曾宗泳, 肖黎明, 等. 大气光学湍流测量中平均时间和原始数据的筛选[J]. 强激光与粒子束, 2004, 16(9): 1101-1105. WENG N Q, ZENG Z Y, XIAO L M, et al. Screening of average time and raw data in atmospheric optical turbulence measurement[J]. High Power Laser and Particle Beams, 2004, 16(9): 1101-1105(in Chinese).
[6] 强希文, 宗飞, 翟胜伟, 等. 大气湍流的实验模拟与测量[J]. 量子电子学报, 2020, 37(4): 506-512. QIANG X W, ZONG F, ZHAI Sh W, et al. Experimental simulation and measurement of atmospheric turbulence[J]. Journal of Quantum Electronics, 2020, 37(4): 506-512(in Chinese).
[7] 曾宗泳, 袁仁民, 谭锟, 等. 复杂地形近地面温度谱[J]. 量子电子学报, 1998, 14(2): 134-139. ZENG Z Y, YUAN R M, TAN K, et al. Near surface temperature spectra of complex terrains[J]. Journal of Quantum Electronics, 1998, 14(2): 134-139(in Chinese).
[8] 许满满, 邵士勇, 刘庆, 等. 复杂地形下海陆风对大气湍流的影响[J]. 光学学报, 2020, 40(12): 1201002. XU M M, SHAO Sh Y, LIU Q, et al. The impact of sea and land breeze on atmospheric turbulence in complex terrain[J]. Acta Optica Sinica, 2020, 40(12): 1201002(in Chinese).
[9] 侯金良. 天文选址的主要参数及测量方法[J]. 天文学进展, 1994, 12(2): 126-132. HOU J L. The main parameters and measurement methods for astronomical site selection[J]. Progress in Astronomy, 1994, 12 (2): 126-132(in Chinese).
[10] 吴晓庆, 王英俭, 曾宗泳, 等. 大气温度起伏谱的测量[J]. 强激光与粒子束, 2002, 14(2): 233-237. WU X Q, WANG Y J, ZENG Z Y, et al. Measurement of atmospheric temperature fluctuation spectrum[J]. High Power Laser and Particle Beams, 2002, 14(2): 233-237(in Chinese).
[11] 王红帅, 李志, 刘威卫. 用于测量圆顶视宁度的微温脉动仪的研制[J]. 天文研究与技术, 2010, 7(1): 60-64. WANG H Sh, LI Zh, LIU W W. Development of a microthermal pulsation meter for measuring the visual acuity of a dome[J]. Astronomical Research and Technology, 2010, 7(1): 60-64(in Chinese).
[12] 杜金林, 张宏升, 张霭琛. 钨铼丝温度脉动仪的制作与应用[J]. 气象水文海洋仪器, 2000(2): 1-6. DU J L, ZHANG H Sh, ZHANG A Ch. The production and application of tungsten rhenium wire temperature pulsation instrument[J]. Meteorological, Hydrological and Oceanic Instruments, 2000(2): 1-6(in Chinese).
[13] AZOUIT M, VERNIN J. Optical turbulence profiling with balloons relevant to astronomy and atmospheric physics[J]. Publications of the Astronomical Society of the Pacific, 2005, 117(831): 536-543. DOI: 10.1086/429785
[14] BARLETTI R, LEMMETTI P, PATERNO L. Balloon-borne radiosonde for the measurement of the atmospheric thermal structure coefficient[J]. Review of Scientific Instruments, 1974, 45(12): 1563-1565. DOI: 10.1063/1.1686562
[15] 宋正方. 应用大气光学基础[M]. 北京: 气象出版社, 1990: 48-51. SONG Zh F. Fundamentals of applied atmospheric optics[M]. Beijing: Meteorological Publishing House, 1990: 48-51(in Chinese).
[16] 任鹏飞, 崔朝龙, 刘庆, 等. 基于宽频带温度脉动仪的大气温度起伏谱测量[J]. 强激光与粒子束, 2019, 31(8): 081001. REN P F, CUI Ch L, LIU Q, et al. Measurement of atmospheric temperature fluctuation spectrum based on broadband temperature pulsation meter[J]. High Power Laser and Particle Beams, 2019, 31(8): 081001(in Chinese).
[17] BEASON M, ANDREWS L, TOSELLI I. Calculating structure Cn2 function constant from measured in non-Kolmogorov and including anisotropic turbulence inner scale effects[J]. Applied Optics, 2019, 58(25): 6813-6819. DOI: 10.1364/AO.58.006813
[18] 李红亮, 许文忠, 王雷, 等. 一种应用于大气边界层大气湍流强度测量的系统[J]. 大气与环境光学学报, 2008, 3(5): 349-354. LI H L, XU W Zh, WANG L, et al. A system for measuring atmospheric turbulence intensity in the atmospheric boundary layer[J]. Journal of Atmospheric and Environmental Optics, 2008, 3(5): 349-354(in Chinese).
[19] 张宏升, 陈家宜, 张霭琛. 自制铂丝温度脉动仪与超声风温仪温度脉动观测的对比研究[J]. 高原气象, 1995, 14(2): 243-249. ZHANG H Sh, CHEN J Y, ZHANG A Ch. Comparative study on temperature fluctuation observation between self-made platinum wire temperature pulsation meter and ultrasonic wind thermometer[J]. Plateau Meteorology, 1995, 14 (2): 243-249(in Chinese).
[20] ISO, OIML, BIPM, et al. Guide to the expression of uncertainty in measurement[M]. Madrid, Spain: Aenor, 1993: 78-79.
[21] 张海鹏, 黄克磊, 马跃. 测量不确定度的评定与表示[J]. 中国医疗器械信息, 2021, 27(5): 19-21. ZHANG H P, HUANG K L, MA Y. Evaluation and representation of measurement uncertainty[J]. China Medical Device Information, 2021, 27 (5): 19-21(in Chinese).
[22] 刘黎, 徐清溪, 肖常国, 等. 食品中菌落总数检测的不确定度评定[J]. 食品安全导刊, 2023(5): 91-93. LIU L, XU Q X, XIAO Ch G, et al. Evaluation of uncertainty in the detection of total bacterial count in food[J]. Journal of Food Safety, 2023(5): 91-93(in Chinese).
[23] 倪育才. 实用测量不确定度评定[M]. 北京: 中国计量出版社, 2005: 31-52. NI Y C. Evaluation of practical measurement uncertainty[M]. Beijing: China Metrology Press, 2005: 31-52(in Chinese).
[24] 姚远成. 影响温度脉动仪测量的因素分析[D]. 合肥: 中国科学院合肥物质科学研究院, 2013: 56-58. YAO Y Ch. Analysis of factors affecting temperature pulsation meter measurement[D]. Hefei: Hefei Institute of Physical Sciences, Chinese Academy of Sciences, 2013: 56-58(in Chinese).
[25] 封双连, 张志刚, 强希文, 等. 戈壁地区近地面大气折射率结构常数的统计分析[J]. 强激光与粒子束, 2012, 24(1): 39-42. FENG Sh L, ZHANG Zh G, QIANG X W, et al. Statistical analysis of the refractive index structure constant of the near surface atmosphere in the Gobi region[J]. High Power Laser and Particle Beams, 2012, 24 (1): 39-42(in Chinese).
计量
- 文章访问数: 2
- HTML全文浏览量: 0
- PDF下载量: 2