CN104237162B - 一种大气折射率虚部结构常数的测量系统及测量方法 - Google Patents
一种大气折射率虚部结构常数的测量系统及测量方法 Download PDFInfo
- Publication number
- CN104237162B CN104237162B CN201410396075.6A CN201410396075A CN104237162B CN 104237162 B CN104237162 B CN 104237162B CN 201410396075 A CN201410396075 A CN 201410396075A CN 104237162 B CN104237162 B CN 104237162B
- Authority
- CN
- China
- Prior art keywords
- lens
- lni
- light
- low frequency
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000012545 processing Methods 0.000 claims abstract description 14
- 230000006641 stabilisation Effects 0.000 claims abstract description 7
- 238000011105 stabilization Methods 0.000 claims abstract description 7
- 238000001228 spectrum Methods 0.000 claims description 19
- 230000007704 transition Effects 0.000 claims description 13
- 230000003287 optical effect Effects 0.000 claims description 12
- 238000010521 absorption reaction Methods 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 238000010183 spectrum analysis Methods 0.000 claims description 5
- 230000002269 spontaneous effect Effects 0.000 claims description 4
- 239000006185 dispersion Substances 0.000 claims description 3
- 230000003321 amplification Effects 0.000 claims description 2
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 2
- 230000006870 function Effects 0.000 description 26
- 238000002474 experimental method Methods 0.000 description 4
- 238000005314 correlation function Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 239000000443 aerosol Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 239000003595 mist Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 239000011358 absorbing material Substances 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 239000005427 atmospheric aerosol Substances 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 230000031700 light absorption Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 230000002459 sustained effect Effects 0.000 description 1
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Optical Communication System (AREA)
Abstract
一种大气折射率虚部结构常数的测量系统及测量方法,包括:发射端,接收端和数据采集存储及处理单元;发射端包括:调制器,驱动电源,发光二极管和发射透镜;调制器对驱动电源进行调制;接收端包括:接收透镜、光电二极管、温度稳定衬垫、前置放大器、检波器和低频放大器;数据采集存储及处理单元包括:数据采集器、存储设备和数据处理模块。本发明能够准确客观地测量大气折射率虚部的结构常数。
Description
技术领域
本发明涉及大气探测设备,具体的说,是一种大气折射率虚部结构常数的测量系统及测量方法。
背景技术
大气折射率虚部结构常数主要用于了解大气折射率虚部的空间变化,进一步用于大气成分(气体和气溶胶)的输送通量测量。
大气有效折射率的数值与光源的波长和大气中的成分有关。大气有效折射率包括两部分:实部和虚部。大气中除气体成分外,还含有云雾和其他气溶胶粒子。光波在大气中传播的时候,会受到气体分子、云雾和气溶胶粒子的散射和吸收。散射与大气折射率实部有关,而吸收与大气折射率的虚部有关。大气折射率在空间的分布是不均匀的。大气折射率结构常数反映了大气折射率的空间变化。这个变化反映了大气湍流状况和吸收物质(如果存在吸收的话)的空间分布,因此,折射率结构常数(实部)的测量首先应用到了和光传输有关的领域,例如天文上视宁度的测量。后来又被广泛用于测量大气近地面热通量。但虚部折射率结构常数的测量和应用较少被人研究,还没有一个实用的方法。
大气折射率虚部与大气中的成分对光波的吸收有关。由于越来越严重的环境污染和人类对气候变化的不断关注,大气折射率虚部的测量显得越来越重要。关于大气虚折射率结构常数,已经开展过一些研究。早在1983年,Filho利用微波(波长为5.4毫米)在伦敦城市中心区域传输4.1千米的距离,根据闪烁频谱分析得到了衰减对闪烁的贡献,与理论估计有很好的一致性(Filho,F.C.M.,Jayasuriya,D.A.R.,Cole,R.S.,andHelmis,C.G.:Spectraldensityofmillimeterwaveamplitudescintillationsinanabsorptionregion,IEEETransactionsonAntennasandPropagation,31,672-676,1983.)。该实验所采用的波长在O2吸收线的波瓣处,而且该波长远大于大气气溶胶的直径,因此该贡献是O2吸收造成的。该实验没有给出由于O2吸收引起的折射率虚部结构常数,但是该研究为折射率虚部结构常数的后续工作打下了基础。Nieveen根据波长为940微米的近红外波传输248米得到的光闪烁数据,得到测量点处的虚折射率结构常数约为4.1×10-24米-2/3(Nieveen,J.P.,Green,A.E.,andKohsiek,W.:Usingalarge-aperturescintillometertomeasureabsorptionandrefractiveindexfluctuations,Boundary-LayerMeteorology,87,101-116,1998.)。由于该实验所采用的波长位于水汽吸收线,从而被认为是由于水汽吸收造成的。该文谈到的折射率虚部结构常数是通过分析闪烁谱密度曲线得到的。根据谱曲线确定一个过渡频率(称之为判断过渡频率法),根据过渡频率与折射率虚部结构常数的关系得到折射率虚部结构常数。由于取样时间长度的限制,谱密度在低频部分具有较大的起伏,无法客观地确定该过渡频率的位置,因而也就无法客观地给出折射率虚部结构常数的数值。
发明内容
本发明技术解决问题:克服判断过渡频率法误差大的缺点,提供一种大气折射率虚部结构常数的测量系统及测量方法,能够准确客观地测量大气折射率虚部的结构常数。
本发明技术方案,一种大气折射率虚部结构常数的测量系统,包括:发射端,接收端和数据采集存储及处理单元;
发射端包括:调制器、驱动电源、发光二极管和发射透镜;调制器对驱动电源进行调制,驱动电源向发光二极管供电,发光二极管发出的光束透过发射透镜发射至大气中;
接收端包括:接收透镜、光电二极管、温度稳定衬垫、前置放大器、检波器和低频放大器;接收端接收发射端发出的光束,经接收透镜将光束会聚到光电二极管上,光电二极管将光信号变为电信号,再依次经前置放大器放大、检波器检波、低频放大器放大到足够的幅度,温度稳定衬垫用于减少光电二极管和前置放大器的温度漂移;
数据采集存储及处理单元包括:数据采集器、存储设备和数据处理模块;数据采集器采集接收端输出的电信号,转换为光强信号存储至存储设备中;数据处理模块处理存储设备中的数据,处理过程包括:
(1)将存储设备中的光强数据取对数,即得到对数光强lnI,通过延迟相关方法或谱分析方法将对数光强起伏方差分解为高频起伏方差和低频起伏方差得到低频起伏方差
(2)利用下面的公式计算垂直于光束的横向风速v,
式中,Dt为发射透镜的直径,Dr为接收透镜的直径;为对数光强起伏方差的高频分量,WPlnI,Re为对数光强起伏的功率谱曲线中平坦过渡区的谱密度;
(3)计算低频起伏的时间延迟为τ的结构函数DlnI,Im(vτ),确定结构函数的系数S;当时间延迟距离vτ大于发射透镜直径Dt和接收透镜直径Dr时,低频起伏的时间延迟为τ的结构函数DlnI,Im(vτ)由公式计算,
式中,为对数光强的高频起伏方差,DlnI(vτ)为对数光强的时间延迟结构函数,计算公式为,
当低频起伏的时间延迟距离vτ小于且接近于湍流外尺度L0时,低频起伏的结构函数DlnI,Im(vτ)表示为,
DlnI,Im(vτ)=S·(vτ)7/6
式中,S为结构函数的系数,利用实验数据确定S的值;
(4)根据前面的计算结果,计算折射率虚部结构常数
低频起伏方差与存在着关系,
式中L0为湍流的外尺度;
低频起伏结构函数系数S与存在着关系,
联立方程,得到折射率虚部结构常数
式中,k为光源发出的光波的波数,k=2π/λ。
一种大气折射率虚部结构常数的测量方法,实现步骤如下:
(1)发射端的调制器对驱动电源进行调制,该驱动电源向发光二极管供电;发光二极管发出的光波的中心波长为λ,发光二极管发出的光束经发射透镜会聚,产生小于0.5°发散角度的光束,采用多组发光二极管和发射透镜发射,增加发射透镜的口径;一个发光二极管和一个发射透镜为一组,采用4组平行并列的方式发射直径达0.6米的光束;
(2)发射端发出的光束穿过待测量区域的大气,光束在大气中传输,会受到大气湍流的散射,气体分子及气溶胶的吸收,引起光强幅度的变化,接收端与发射端的距离为L;L为500米至10公里;
(3)接收端的接收透镜接收来自发射端的光束,将光信号会聚在光电二极管上,采用多个透镜接收,采用4个透镜并列放置,使得所有透镜接收到的光束均会聚至光电二极管,同样是为了增加接收透镜的口径;光电二极管将光信号变为电信号,首先经过前置放大器放大,再经过检波器检波,然后再经过低频放大器放大到足够的幅度,得到输出信号;
(4)接收端输出的电信号经数据采集器采集,转换为光强,然后将光强数据存储至存储设备;
(5)将存储设备中的光强数据取对数,即得到对数光强lnI,通过延迟相关方法或谱分析方法将对数光强起伏方差分解为高频起伏方差和低频起伏方差得到低频起伏方差
(6)利用下面的公式计算垂直于光束的横向风速v,
式中,Dt为发射透镜的直径,Dr为接收透镜的直径;为对数光强起伏方差的高频分量,WPlnI,Re为对数光强起伏的功率谱曲线中平坦过渡区的谱密度;
(7)计算低频起伏的时间延迟为τ的结构函数DlnI,Im(vτ),确定结构函数的系数S;当时间延迟距离vτ大于发射透镜直径Dt和接收透镜直径Dr时,低频起伏的时间延迟为τ的结构函数DlnI,Im(vτ)由公式计算,
式中,为对数光强的高频起伏方差,DlnI(vτ)为对数光强的时间延迟结构函数,计算公式为,
当低频起伏的时间延迟距离vτ小于且接近于湍流外尺度L0时,低频起伏的结构函数DlnI,Im(vτ)表示为,
DlnI,Im(vτ)=S·(vτ)7/6
式中,S为结构函数的系数,利用实验数据确定S的值;
(8)根据前面的计算结果,计算折射率虚部结构常数
低频起伏方差与存在着关系,
式中L0为湍流的外尺度;
低频起伏结构函数系数S与存在着关系,
联立方程,得到折射率虚部结构常数
式中,k为光源发出的光波的波数,k=2π/λ。
本发明与现有技术相比的优点在于:
发射端选用特定波长的光源,可以得到和气溶胶吸收相关的折射率虚部结构常数;同时为通过改变光源波长测量水汽相关的折射率虚部结构常数打下基础;
采用多组发光二极管和发射透镜发射,增加了发射透镜的口径,采用多个接收透镜接收,增加了接收透镜的口径,增加了接收信号的能量,同时也减小了对数光强的高频起伏方差,从而提高了低频闪烁信号的测量精度;
将对数光强起伏分解为高频起伏和低频起伏,得到相应的方差和结构函数,利用联立方程求解的方法得到折射率虚部结构常数,克服了现有的判断过渡频率法误差大的缺点,能够具有较高精度、客观地获得折射率虚部结构常数。
附图说明
图1为本发明测量系统的组成框图;测量系统包括发射端,接收端和数据采集存储单元。图中显示了含有空气分子,气溶胶粒子组成的待测大气;
图2为本发明测量系统中的发射端组成原理图;
图3为本发明测量系统中的接收端和数据采集组成原理图;
图4为本发明方法的流程图;
图5为利用时间延迟相关法将对数光强起伏方差分解为高频和低频分量的原理图;
图6为利用功率谱方法得到过渡平坦部分的谱密度的原理图;
图7为计算低频起伏结构函数的原理图;
图8为一次测量的实例,折射率虚部结构常数随时间的变化曲线。
具体实施方式
如图1所示,本发明实施的大气折射率虚部结构常数测量系统包括发射端,接收端和数据采集存储处理单元。
发射端包括调制器,驱动电源,发光二极管和发射透镜。调制器对驱动电源进行调制。驱动电源向发光二极管供电,发光二极管发出的光束透过发射透镜发射至大气中,如图2所示。本发明发射端采用一组LED灯作为光源,LED灯供电部分经过调制,使得发射光束避免了背景光的干扰。调制频率可以为2~120KHz。LED灯发出的光经发射透镜会聚,以便传输较远距离。
接收端包括接收透镜,光电二极管,温度稳定衬垫,前置放大器,检波器,低频放大器,如图3所示。接收端接收发射端发出的光束,经接收透镜将光束会聚到光电二极管上。光信号变为电信号,经前置放大器放大,再经检波器检波,然后再经低频放大器放大到足够的幅度。温度稳定衬垫用于减少光电二极管和前置放大器的温度漂移。
数据采集存储单元由数据采集器,存储设备组成。该单元采集接收端输出的电信号,根据光电转换系数转变为光强数据存储至存储设备。数据处理软件处理存储的光强起伏数据,得到大气折射率虚部结构常数。
本发明提供了一种大气折射率虚部结构常数测量方法,选择不同波长的光源,反映的是不同大气成分(吸收所使用的光波)的折射率虚部的空间变化。
如图4所示,本发明中提出的大气折射率虚部结构常数的测量包括以下步骤:
(1)发射端的调制器对驱动电源进行调制,该驱动电源向发光二极管供电;发光二极管发出的光波的中心波长为λ,发光二极管发出的光束经发射透镜会聚,产生小于0.5°发散角度的光束,采用多组发光二极管和发射透镜发射,增加发射透镜的口径;一个发光二极管和一个发射透镜为一组,采用4组平行并列的方式发射直径达0.6米的光束;
(2)发射端发出的光束穿过待测量区域的大气,光束在大气中传输,会受到大气湍流的散射,气体分子及气溶胶的吸收,引起光强幅度的变化,将接收端与发射端的距离为L,L可以为500米至10公里;
(3)接收端的接收透镜接收来自发射端的光束,将光信号会聚在光电二极管上,采用多个透镜接收(采用4个透镜并列放置,使得接收到的光束均会聚至光电二极管),同样是为了增加接收透镜的口径;光电二极管将光信号变为电信号,首先经过前置放大器放大,再经过检波器检波,然后再经过低频放大器放大到足够的幅度,得到输出信号;
(4)接收端输出的电信号经数据采集器采集,采样频率为500Hz。将点信号转换为光强数据,然后存储至存储设备;每20~30分钟保存一个文件;
(5)将光强数据取对数,即得到对数光强(lnI)。通过延迟相关方法将对数光强起伏方差分解为高频起伏方差(用表示,因为高频起伏方差由折射率实部引起)和低频起伏方差(用表示,因为低频起伏方差由折射率虚部引起)。计算时间延迟为τ的相关函数ClnI(τ),
式中,t为测量的时刻,上划线表示对时间t的平均。如图5所示,空心圆圈组成了相关函数曲线。根据图5的标注,利用时间延迟相关函数曲线可以计算得到高频起伏方差和低频起伏方差
(6)计算得到对数光强起伏的功率谱,如图6所示。图6的空心圆圈组成了对数光强起伏的功率谱曲线。根据对数光强起伏的功率谱曲线得到平坦过渡部分的谱密度WPlnI,Re。图6中水平虚线表示谱密度WPlnI,Re;
(7)利用下面的公式计算垂直于光束的横向风速v,
式中,Dt为发射透镜的直径,Dr为接收透镜的直径;WPlnI,Re为对数光强起伏的功率谱曲线中平坦过渡区的谱密度;
(8)计算低频起伏的时间延迟为τ的结构函数DlnI,Im(vτ),确定结构函数的系数S;
当时间延迟距离vτ大于发射透镜直径Dt和接收透镜直径Dr时,低频起伏的时间延迟为τ的结构函数DlnI,Im(vτ)由公式计算,
式中,为对数光强的高频起伏方差。DlnI(vτ)为对数光强的时间延迟结构函数,计算公式为,
如图7所示,图中的空心圆圈组成了对数光强的时间延迟结构函数曲线。水平的虚线表示根据(3)式计算得到的低频起伏的时间延迟为τ的结构函数为图中的实线。
当低频起伏的时间延迟距离vτ小于且接近于湍流外尺度L0时,低频起伏的结构函数DlnI(vτ)可以表示为,
DlnI,Im(vτ)=S·(vτ)7/6(5)
上式中的S为低频起伏结构函数的系数,利用实验数据可以确定S的数值;S即为如图7中实线所示的低频起伏结构函数DlnI,Im曲线的斜率(时间延迟距离为2~20米范围内)。
(9)根据前面的计算结果,利用下面的公式计算折射率虚部结构常数
步骤(5)中的低频起伏方差与存在着关系,
式中L0为湍流的外尺度,L为光束自发射端到接收端的传输距离。
步骤(8)中的低频起伏结构函数系数S与存在着关系,
公式(6)和(7)联立方程,得到折射率虚部结构常数
式中,k为光源发出的光波的波数,k=2π/λ,L为发射端至接收端的传输距离。
实施方式举例
以在中国科学技术大学校园内的测量为例,校园在安徽省合肥市城内,校园周边有若干交通要道,来往车辆较多。将发射端和接收端置于离地面10米高度以上,接收端离发射端几百米至数公里的距离。如果在城市开展测量,可以选择两栋楼房。发射端和接收端直线路径上不能有任何遮蔽和阻挡。发射端和接收端可大致处于同一高度。这样处理使得测量结果的应用较为简单,当测量高度较低,或者两端高差较大时,使大气折射率虚部结构常数的应用较为复杂,但并不影响结果的使用。
下面介绍测量结果并给出简单结论。
20140115~16日在中国科学技术大学校园开展了折射率虚部结构常数的测量。该实验的发射透镜和接收透镜的口径均为0.18米,光源的波长为0.620微米。发射透镜与接收透镜的距离为960米。采样频率为500Hz,每20分钟保存一个数据文件。根据前面的步骤得到20140115~16日折射率虚部结构常数随时间变化,见图8。由图8可以看到折射率虚部结构常数并不像折射率实部结构常数存在明显的日变化特征(Stull,1988),但是可以看到1月15日上午9:00左右折射率虚部结构常数最大,1月16日中午12点最大,其他时间起伏也较大。进一步分析可以说明折射率虚部结构常数不仅与湍流的变化有关,还与气溶胶的空间分布有关。
提供以上实施例仅仅是为了描述本发明的目的,而并非要限制本发明的范围。本发明的范围由所附权利要求限定。不脱离本发明的精神和原理而做出的各种等同替换和修改,均应涵盖在本发明的范围之内。
Claims (2)
1.一种大气折射率虚部结构常数的测量系统,其特征在于包括:发射端,接收端和数据采集存储及处理单元;
发射端包括:调制器、驱动电源、发光二极管和发射透镜;调制器对驱动电源进行调制,驱动电源向发光二极管供电,发光二极管发出的光束透过发射透镜发射至大气中;
接收端包括:接收透镜、光电二极管、温度稳定衬垫、前置放大器、检波器和低频放大器;接收端接收发射端发出的光束,经接收透镜将光束会聚到光电二极管上,光电二极管将光信号变为电信号,再依次经前置放大器放大、检波器检波、低频放大器放大到足够的幅度,温度稳定衬垫用于减少光电二极管和前置放大器的温度漂移;
数据采集存储及处理单元包括:数据采集器、存储设备和数据处理模块;数据采集器采集接收端输出的电信号,转换为光强信号存储至存储设备中;数据处理模块处理存储设备中的数据,处理过程包括:
(1)将存储设备中的光强数据取对数,即得到对数光强lnI,通过延迟相关方法或谱分析方法将对数光强起伏方差分解为高频起伏方差和低频起伏方差得到低频起伏方差
(2)利用下面的公式计算垂直于光束的横向风速v,
式中,Dt为发射透镜的直径,Dr为接收透镜的直径;为对数光强起伏方差的高频分量,WPlnI,Re为对数光强起伏的功率谱曲线中平坦过渡区的谱密度;
(3)计算低频起伏的时间延迟为τ的结构函数DlnI,Im(vτ),确定结构函数的系数S;当时间延迟距离vτ大于发射透镜直径Dt和接收透镜直径Dr时,低频起伏的时间延迟为τ的结构函数DlnI,Im(vτ)由公式计算,
式中,为对数光强起伏方差的高频分量,DlnI(vτ)为对数光强的时间延迟结构函数,计算公式为,
当低频起伏的时间延迟距离vτ小于且接近于湍流外尺度L0时,低频起伏的时间延迟为τ的结构函数DlnI,Im(vτ)表示为,
DlnI,Im(vτ)=S·(vτ)7/6
式中,S为结构函数的系数,利用实验数据确定S的值;t为测量的时刻,上划线表示对时间t的平均;
(4)根据前面的计算结果,计算折射率虚部结构常数
低频起伏方差与存在着关系,
式中L0为湍流的外尺度;L为光束自发射端到接收端的传输距离;
结构函数的系数S与存在着关系,
联立方程,得到折射率虚部结构常数
式中,k为光源发出的光波的波数,k=2π/λ,λ为波长,L为光束自发射端到接收端的传输距离。
2.一种大气折射率虚部结构常数的测量方法,其特征在于实现步骤如下:
(1)发射端的调制器对驱动电源进行调制,该驱动电源向发光二极管供电;发光二极管发出的光波的中心波长为λ,发光二极管发出的光束经发射透镜会聚,产生小于0.5°发散角度的光束,采用四组发光二极管和发射透镜发射,增加发射透镜的口径;一个发光二极管和一个发射透镜为一组,采用四组平行并列的方式发射直径达0.6米的光束;
(2)发射端发出的光束穿过待测量区域的大气,光束在大气中传输,会受到大气湍流的散射,气体分子及气溶胶的吸收,引起光强幅度的变化,接收端与发射端的距离为L;L为500米至10公里;
(3)接收端的接收透镜接收来自发射端的光束,将光信号会聚在光电二极管上,采用四个透镜接收,采用四个透镜并列放置,使得所有透镜接收到的光束均会聚至光电二极管,同样是为了增加接收透镜的口径;光电二极管将光信号变为电信号,首先经过前置放大器放大,再经过检波器检波,然后再经过低频放大器放大到足够的幅度,得到输出信号;
(4)接收端输出的电信号经数据采集器采集,转换为光强,然后将光强数据存储至存储设备;
(5)将存储设备中的光强数据取对数,即得到对数光强lnI,通过延迟相关方法或谱分析方法将对数光强起伏方差分解为高频起伏方差和低频起伏方差得到低频起伏方差
(6)利用下面的公式计算垂直于光束的横向风速v,
式中,Dt为发射透镜的直径,Dr为接收透镜的直径;为对数光强起伏方差的高频分量,WPlnI,Re为对数光强起伏的功率谱曲线中平坦过渡区的谱密度;
(7)计算低频起伏的时间延迟为τ的结构函数DlnI,Im(vτ),确定结构函数的系数S;当时间延迟距离vτ大于发射透镜直径Dt和接收透镜直径Dr时,低频起伏的时间延迟为τ的结构函数DlnI,Im(vτ)由公式计算,
式中,为对数光强起伏方差的高频分量,DlnI(vτ)为对数光强的时间延迟结构函数,计算公式为,
当低频起伏的时间延迟距离vτ小于且接近于湍流外尺度L0时,低频起伏的时间延迟为τ的结构函数DlnI,Im(vτ)表示为,
DlnI,Im(vτ)=S·(vτ)7/6
式中,S为结构函数的系数,利用实验数据确定S的值;t为测量的时刻,上划线表示对时间t的平均;
(8)根据前面的计算结果,计算折射率虚部结构常数
低频起伏方差与存在着关系,
式中L0为湍流的外尺度;
结构函数的系数S与存在着关系,
联立方程,得到折射率虚部结构常数
式中,k为光源发出的光波的波数,k=2π/λ,λ为波长,L为光束自发射端到接收端的传输距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410396075.6A CN104237162B (zh) | 2014-08-12 | 2014-08-12 | 一种大气折射率虚部结构常数的测量系统及测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410396075.6A CN104237162B (zh) | 2014-08-12 | 2014-08-12 | 一种大气折射率虚部结构常数的测量系统及测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104237162A CN104237162A (zh) | 2014-12-24 |
CN104237162B true CN104237162B (zh) | 2016-07-27 |
Family
ID=52225669
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410396075.6A Expired - Fee Related CN104237162B (zh) | 2014-08-12 | 2014-08-12 | 一种大气折射率虚部结构常数的测量系统及测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104237162B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105158210B (zh) * | 2015-08-31 | 2018-01-30 | 中国科学技术大学 | 一种大气气溶胶质量垂直输送通量的测量方法及系统 |
CN108760226B (zh) * | 2018-05-04 | 2020-01-21 | 西华大学 | 一种探测大气湍流结构常数的方法及装置 |
CN109187357A (zh) * | 2018-09-04 | 2019-01-11 | 南京大学 | 一种大气能见度与湍流同步测量系统及测量方法 |
CN112525826A (zh) * | 2020-11-13 | 2021-03-19 | 中国科学院合肥物质科学研究院 | 一种大气光学湍流强度测量、评估和修正方法及系统 |
CN118095964B (zh) * | 2024-04-23 | 2024-08-09 | 成都飞机工业(集团)有限责任公司 | 一种航空管路系统装配一致性的评估方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101109702A (zh) * | 2007-07-28 | 2008-01-23 | 中国科学院安徽光学精密机械研究所 | 测量大气折射率结构常数的大口径激光闪烁仪及方法 |
JP2010210384A (ja) * | 2009-03-10 | 2010-09-24 | Shimadzu Corp | 屈折率計 |
CN101957317A (zh) * | 2010-09-02 | 2011-01-26 | 长春理工大学 | 大气湍流折射率结构常数高度分布模式测量装置 |
CN101980056A (zh) * | 2010-11-06 | 2011-02-23 | 四川大学 | 基于m2因子和光闪烁指数确定大气湍流参数方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102879358B (zh) * | 2012-10-31 | 2016-12-21 | 电子科技大学 | 一种基于光延迟技术的介质折射率测量装置及其测量方法 |
-
2014
- 2014-08-12 CN CN201410396075.6A patent/CN104237162B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101109702A (zh) * | 2007-07-28 | 2008-01-23 | 中国科学院安徽光学精密机械研究所 | 测量大气折射率结构常数的大口径激光闪烁仪及方法 |
JP2010210384A (ja) * | 2009-03-10 | 2010-09-24 | Shimadzu Corp | 屈折率計 |
CN101957317A (zh) * | 2010-09-02 | 2011-01-26 | 长春理工大学 | 大气湍流折射率结构常数高度分布模式测量装置 |
CN101980056A (zh) * | 2010-11-06 | 2011-02-23 | 四川大学 | 基于m2因子和光闪烁指数确定大气湍流参数方法及装置 |
Non-Patent Citations (4)
Title |
---|
Effects of Non-Uniform Crosswind Fields on Scintillometry Measurements;H. C. Ward et al.;《Boundary-Layer Meteorol》;20110628;第141卷;第143–163页 * |
Using a large-aperture scintillometer to measure absorption and refractive index fluctuations;J. P. NIEVEEN et al.;《Boundary-Layer Meteorology》;19981231;第101–116页 * |
不同方法测量折射率结构常数的比较;袁仁民 等;《光学学报》;20000630;第20卷(第6期);第755-761页 * |
孔径接收下大气闪烁频谱的理论和实验研究;潘锋 等;《强激光与粒子束》;20060930;第18卷(第9期);第1457-1459页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104237162A (zh) | 2014-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104237162B (zh) | 一种大气折射率虚部结构常数的测量系统及测量方法 | |
CN105424631B (zh) | 一种基于紫外可见波段吸收光谱的超高灵敏度氮氧化物测量系统 | |
CN100545633C (zh) | 激光气体遥测的方法和装置 | |
CN101109702B (zh) | 测量大气折射率结构常数的大口径激光闪烁仪及方法 | |
CN102879359B (zh) | 一种大气能见度测量系统 | |
CN103983544A (zh) | 多通道气溶胶散射吸收测量仪 | |
CN105352914B (zh) | 一种基于双波长光纤环形腔的气体浓度检测系统及其方法 | |
CN103344614B (zh) | 一种高精度大气透过率测量装置及测量方法 | |
CN106596437A (zh) | 大气no3自由基浓度在线测量系统和在线测量方法 | |
CN104596645A (zh) | 复杂环境全偏振多角度散射模拟及测试系统 | |
CN103048284A (zh) | 利用镀膜玻璃管测量大气气溶胶消光系数的新方法 | |
CN105823755A (zh) | 一种基于可调谐半导体激光的自混合气体吸收传感系统 | |
CN105067563A (zh) | 开放空间气体平均浓度测量装置及测量方法 | |
CN106226783A (zh) | 基于激光雷达的大气颗粒物光学参数测量系统 | |
CN103528991B (zh) | 土壤有机质含量的测量系统及测量方法 | |
Ding et al. | Modeling and characterization of ultraviolet scattering communication channels | |
CN109507074A (zh) | 一种超低排放烟尘浓度监测装置和监测方法 | |
CN109883901A (zh) | 一种无人机机载的紫外光探测雾霾粒子系统及其探测方法 | |
CN104568818A (zh) | 一种基于光纤传导的主动式太赫兹光谱检测内窥探头 | |
CN102269699B (zh) | 禽舍硫化氢气体浓度检测方法 | |
Siegel et al. | Short-range communication with ultraviolet LEDs | |
CN106483531B (zh) | 大气拉曼-瑞利散射测温激光雷达及反演方法 | |
Li et al. | SNR and transmission error rate for remote laser communication system in real atmosphere channel | |
CN105158210B (zh) | 一种大气气溶胶质量垂直输送通量的测量方法及系统 | |
Hawe et al. | CO2 monitoring and detection using an integrating sphere as a multipass absorption cell |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160727 |
|
CF01 | Termination of patent right due to non-payment of annual fee |