大气探测学chpt8风的观测

1、第八章第八章 风的观测风的观测8.1 8.1 概述概述 8.2 8.2 风向的测量原理风向的测量原理 8.3 8.3 旋转风速仪原理旋转风速仪原理 8.4 8.4 散热式风速仪原理散热式风速仪原理8.5 8.5 声学风速仪原理声学风速仪原理 8.6 8.6 激光风速仪原理激光风速仪原理 8.7 8.7 风速仪的检定风速仪的检定 8.8 气象铁塔探测系统气象铁塔探测系统8.9 8.9 系留气球探测系留气球探测 8.10 8.10 风廓线雷达风廓线雷达8.118.11复习思考题复习思考题风向风速仪风速表和风向标微气象观测系统微气象观测系统风杯风速计和风向标风杯风速计和风向标风杯风速计螺旋浆式风速风

2、向计螺旋浆式风速风向计三分量风速仪三分量风速仪大气探测观测站16M风风温温湿湿梯梯度度塔塔全自动资料采集系统全自动资料采集系统太阳辐射计太阳辐射计超声湍流观测系统温湿温湿度、度、地温地温观测观测系统系统近地层大气多层风、温、湿与通量观测超声风速仪8.1 8.1 概述概述一、基本知识一、基本知识 1.1.风：大气相对地面运动的水平分量叫作风。风：大气相对地面运动的水平分量叫作风。 2.2.风速：单位时间气流运动距离叫做风速。风速：单位时间气流运动距离叫做风速。 3.3.风向：风的来向，共计风向：风的来向，共计1616个方位。个方位。二、为什么测风？二、为什么测风？三、测风仪器三、测风仪器 1.1

3、.台站：台站：ELEL电接风向风速计、手持风向风速电接风向风速计、手持风向风速计、雷达测风、气球测风；计、雷达测风、气球测风； 2.2.研究：超声风速计、三轴风速计、激光雷研究：超声风速计、三轴风速计、激光雷达、热线风速仪、气球测风；达、热线风速仪、气球测风；1、风向的测量、风向的测量1616个方位个方位风向用风向用16方位表示与方位角度的对应关系方位表示与方位角度的对应关系风速与风力等级风的测量中存在的问题 测风仪器安装：开阔、平坦测风仪器安装：开阔、平坦 距地面距地面10米高度米高度 仪器磨损：定期检定仪器磨损：定期检定 风的脉动：时间概念，峰值和统计量风的脉动：时间概念，峰值和统计量 风

4、的平均：过高效应；风的平均：过高效应； 0和和360之间的切换之间的切换 风向和风速测量的匹配风向和风速测量的匹配自动站风测量规范自动站风测量规范 1 分钟和分钟和 2 分钟平均风分钟平均风 （ 3 秒钟为步长的滑动平均）秒钟为步长的滑动平均） 10 分钟平均风（分钟平均风（1 分钟为步长的滑动平均）分钟为步长的滑动平均） 1 小时内最大风和出现时间（正点时从每小时内最大风和出现时间（正点时从每分钟的分钟的10分钟平均风速中选取最早出现的分钟平均风速中选取最早出现的最大风速和相应的风向）最大风速和相应的风向） 1 小时内极大风和出现时间（正点时从小时内极大风和出现时间（正点时从 3 秒钟平均风

5、速中选取最早出现的极大风速秒钟平均风速中选取最早出现的极大风速和相应的风向）和相应的风向）8.2 8.2 风向的测量原理风向的测量原理8.2.1.8.2.1.风向标及风向测量风向标及风向测量风向标由指向杆、平衡重锤、旋转轴和尾翼组成。风向标由指向杆、平衡重锤、旋转轴和尾翼组成。增加灵敏度：减少轴摩擦、加大尾翼、加长后杆。增加灵敏度：减少轴摩擦、加大尾翼、加长后杆。 对风向标的结构和造型主要考虑：对风向标的结构和造型主要考虑： 1.1.灵敏性：有良好的起动性能。灵敏性：有良好的起动性能。 2.2.稳定性：有良好的动态性能，迅速跟踪风向变化。稳定性：有良好的动态性能，迅速跟踪风向变化。8.2.2

6、8.2.2 风向标的动态特征风向标的动态特征当风向变动时，风向标必须迅速作出反应。假设风向标偏离当风向变动时，风向标必须迅速作出反应。假设风向标偏离风向的角度为风向的角度为 ，风尾板上受到一个力为，风尾板上受到一个力为 ，力作用中心距旋，力作用中心距旋转轴的力臂为转轴的力臂为 ，则，则单位角度风向所产生的扭力矩为：单位角度风向所产生的扭力矩为：在风力作用下，风向标转动角速度为：在风力作用下，风向标转动角速度为：风标在线速度方向相对于空气运动的速度为：风标在线速度方向相对于空气运动的速度为：/NF./ddtsin/uddt 为平均风速在线速度方向上的投影分量；为平均风速在线速度方向上的投影分量；

7、 产生的线速度（产生的线速度（ ）。作用于风标上风矢量的迎角为）。作用于风标上风矢量的迎角为 ，而，而 理论上应为理论上应为（8.2）应用了）应用了当当很小时，很小时，sinsin=tg=tg= =，coscos=1=1的假设。的假设。.11sintantancosduudtuuusinudvdt.ddt（8.1）（8.2）vFvr为在角速度为在角速度vrvvv得风向标的运动方程：得风向标的运动方程：22vvNdddJFNNNDdtudtdt上式右边第一项是气流对风标施加的扭力矩，第二项是空气对上式右边第一项是气流对风标施加的扭力矩，第二项是空气对运动风标的阻尼力矩。如果运动风标的阻尼力矩。如

8、果N,D为常数，得风向标运动的二阶为常数，得风向标运动的二阶微分方程：微分方程：.0JDN式中式中2.22;ddtNmu为风向标转动角加速度;J为转动惯量J=m为风标质量;D=为空气动力阻尼.（8.3）（8.4）00exp22exp()exp( 2)2ddDttiJtDttiJt其中：其中：0 0为为t t0 0时，风向标的偏离角；时，风向标的偏离角； 为风向标的阻尼谐振周期：为风向标的阻尼谐振周期：21 / 22()() 2dtNDJJ(8.4)(8.4)式的解为：式的解为：振幅衰减项振幅衰减项 谐振运动项谐振运动项dt（8.6）（8.5）(8.6)(8.6)式是一个典型的周期衰减振动。如下

9、图所示。式是一个典型的周期衰减振动。如下图所示。单衰减运动谐振衰减运动有关风向标动态特征的几个重要参数：有关风向标动态特征的几个重要参数： 120001.0,0,21,2,DDJttNJtd无阻尼;Dexp不衰减,2Jt为无阻尼谐振周期.202.,exp21,22dDDtiJJtNJ000当时风标为单纯的衰减运动,此时风标的阻尼D=D ,D 为临界阻尼,D（8.7）3.阻尼比：风标阻尼阻尼比：风标阻尼D与临界阻尼与临界阻尼D0的比值的比值000vvNDDu Dut推导：22001()()422NNNNNDDJNJNJ1/ 201()2NNDJ01/ 22()tNJ00NDt=1，D=D0，临界

10、阻尼，迅速达到平衡；1，DD0，过阻尼，缓慢达到平衡；1，DD0，欠阻尼，在平衡处振动。由由:（8.9）（8.8）P.102(5.8)在实际工作中，应用谐振周期在实际工作中，应用谐振周期t td d，无阻尼谐振周期，无阻尼谐振周期t t0 0不方便，不方便，定义一个时间尺度定义一个时间尺度t tL L, ,认为风向标经过认为风向标经过t tL L时间后时间后, ,风向标偏离风风向标偏离风向的角度由原先向的角度由原先t=0t=0时的时的0010,tLtLtL当： 时0exp()exp( 2)2dDttiJt01exp()2tLLDtLJ1ln1ln,2/2LLDLttLJDJ vNDu1ln4.

11、61lg2LvvJuLtLNNJu 0tL衰减为由由不考虑振动项不考虑振动项则则:（8.10）P.102（5.9）根据实验结果，扭力矩根据实验结果，扭力矩N的表达式为：的表达式为：212vvNu sa式中式中: 为空气密度为空气密度 ; S为风尾板有效面积为风尾板有效面积; 为扭力矩系数为扭力矩系数( ), 为空气动力系数为空气动力系数; 为为流体的动压流体的动压; ;1 1、风程、风程t t0 0为无阻尼风标的固有周期为无阻尼风标的固有周期;Ut;Ut0 0为单位时间为单位时间t t，风标振动一次，风标振动一次，气流所走距离气流所走距离( (风程风程).).由由(8.7)(8.7)式和式和(

12、8.11)(8.11)式得式得,t=0,t=0时的风程时的风程: :（8.11）aCC212u1/21/201/21/221/2222(/)(/2)vvuuuutJJN JNu sr a1/27.95()vvJr a s（8.12）P.102（5.11）31.25/kg m2 2、阻尼比、阻尼比1 / 21 / 20()22/()vvvrrrNNu tuJuJ231/ 21/ 21/ 2/ 211()()()222vvvvvru sr a sa r suJJ31/20.395()v va r sJ（8.13）P.102（5.12）3 3、Ltt时的风程 22224.61lg4.61lglg7.

13、3112LL JuL JuL JutNa Su Sa（8.14）P.102（5.13）8.2.3 8.2.3 风向标平衡重锤对风向标动态特性的影响风向标平衡重锤对风向标动态特性的影响由于平衡锤的运动方向与其空气动力矩相反，则有由于平衡锤的运动方向与其空气动力矩相反，则有wwvwr N dNNudt由于指向杆和平衡锤的受风面积很小，相对于风尾板所受的扭力由于指向杆和平衡锤的受风面积很小，相对于风尾板所受的扭力矩矩NW的作用可忽略不计。但是它的转动惯量必须考虑。风标系的作用可忽略不计。但是它的转动惯量必须考虑。风标系统对主轴的支点重力矩是平衡的，即统对主轴的支点重力矩是平衡的，即从描述风向标动态特

14、性的三个重要参数的（从描述风向标动态特性的三个重要参数的（8.12)、(8.13)和和(8.14)式中可看出都有转动惯量式中可看出都有转动惯量J。因而，风向标的平衡锤必然影响风向。因而，风向标的平衡锤必然影响风向标的动态特性。假设平衡锤在气流中也要受到一个扭力矩标的动态特性。假设平衡锤在气流中也要受到一个扭力矩NW，整，整个风标系统受到的总扭力矩为个风标系统受到的总扭力矩为NT=N-NW （因平衡锤纽力矩方向与（因平衡锤纽力矩方向与风标扭力矩方向相反）风标扭力矩方向相反）wwmm那么，整个平衡系统的转动惯量则可写作那么，整个平衡系统的转动惯量则可写作 222(1)wvvwwvvvrJm rm

15、rm rr21()11/vv vvwvm a r ssarr21(1)1/v vvvwva r sarr设：设： 为风尾单位面积的重量。并令下述因子为风为风尾单位面积的重量。并令下述因子为风向标质量因子：向标质量因子：11/vvvwvaKrr2vvva r sJK（8.15）P.103（5.14）mS因而，（因而，（8.15）式可写为）式可写为（8.16）（8.17）P.103（5.15）P.103（5.16） 将将(8.17)代入代入(8.12)、(8.13)和和(8.14)式，风向标的几个主要动式，风向标的几个主要动态特性参数改写为：态特性参数改写为：21/21/21/207.95()7.

16、95()7.95()vvvvvvvvvasJuta sa sKK13321/21/220.395()0.3950.395()vvvvasa SKKJa St=t0时的风程：时的风程： 222lglg( )lg( )7.377.377.37LvvvLL a SJLutr a sa SKK阻尼比：阻尼比：t=tL时的风程：时的风程：（8.18）P.103（5.17）P.103（5.18）P.103（5.19）（8.19）（8.20）希望风标系统动态跟踪风向变动的性能要好，即希望风标系统动态跟踪风向变动的性能要好，即 的值较小，关键的值较小，关键在于有数值较大的质量因子在于有数值较大的质量因子 ，从

17、式（，从式（8.16）可分析得到：）可分析得到：1、 值较大，即气流动压值较大，即气流动压 能有效作用到风尾板上；能有效作用到风尾板上；2、 值较小，单位面积风尾板的质量小。意味着在制造风尾板时，值较小，单位面积风尾板的质量小。意味着在制造风尾板时， 必须选用强度高、比重小的材料；必须选用强度高、比重小的材料；3、 1，临界阻尼=1，欠阻尼11时：时：21/221/2000021/2exp()cos(1)sin(1)(1)ttt 当=1时：000exp()(1)tt当1时：2 1/20002 1/21exp()cos(1)(1)tt （8.23）（8.24）（8.25）21/221/20:ar

18、ctan(1)(1)其中为观测的角频率，实际频率的为观测的角频率，实际频率的2 2倍，代入倍，代入 (5.25)(5.25)式式, ,则有：则有：021/221/21expcos()(1)(1)tt121/ 22(1)2exp(1)nnnn2 1/22exp()(1)（5.27） （8.28）因为方程因为方程 ） 描述的是一个以描述的是一个以 为周期的振动运动，为周期的振动运动，所以，相邻两个最大值之比可表示为：所以，相邻两个最大值之比可表示为：（8.28）2（8.29）（5.26） 它的解可以表示为通解与特解之和它的解可以表示为通解与特解之和,即：即：2220012sinddAdtdt21/

19、221/2100200222 1/200exp(1)exp(1)sin()4()1 () ct itct itAt 其中：02202arctan()当二阶微分方程有一个外力函数当二阶微分方程有一个外力函数 作用其上时作用其上时,即：即： sinf tA（8.30）（8.31）（8.32）当1时，风向标指示振幅大于实际振幅时，风向标指示振幅大于实际振幅; G1时，风向标指示振幅小于实际振幅时，风向标指示振幅小于实际振幅.（8.33）（8.34）当=0，阻尼比很小时，易于发生共振.G G与风标特性参数的关系见与风标特性参数的关系见P105P105页图页图5.65.6所示所示, ,由图表明由图表明

20、1.在/0 =1， 1时，风向标发生共振; 2.在/0 1，1，G 1时，风向标振幅小于实际风向振幅.8.2.5 8.2.5 风向标动态参数的选择风向标动态参数的选择lg7.37LvLutKlglg 2.718287.377.371.242.58vLLKut一个性能良好的风向标应能准确地反应连续变化的风向一个性能良好的风向标应能准确地反应连续变化的风向.另外另外,为了使为了使风向资料具有可比较性风向资料具有可比较性,世界气象组织原则上作了如下规定世界气象组织原则上作了如下规定:1. 在风速为2.58米/秒时,风向标应在1s内使风向偏差衰减到1/e,即001,.102,1,tLLtLtLPteL

21、L0根据页,定义一个风向标的时间尺度t 认为风向标经过时间之后风向标偏离风向的角度由原先t=0时刻的衰减为由8.20式质量因子取1.25K2. 阻尼比取阻尼比取 : 0.3 1.75. 1.75.8.3 旋转风速计原理旋转风速计原理2004/11/228.3.1 8.3.1 测风速仪器原理概述测风速仪器原理概述1.1.蜗杆风速表蜗杆风速表( (齿轮转速齿轮转速风速的关系风速的关系););2.2.电感式电感式( (电感、电压电感、电压风速的关系）风速的关系）; ;3.3.光电式（光电频率光电式（光电频率风速的关系）风速的关系）; ;4.4.悬浮式（光电频率悬浮式（光电频率风速的关系）风速的关系）

22、; ;5.5.超声风速仪（接收声波时间超声风速仪（接收声波时间与风速的关系）与风速的关系）8.3 旋转风速计原理旋转风速计原理8.3.2 8.3.2 旋转风杯风速计的感应原理旋转风杯风速计的感应原理 旋转风杯风速计是风速测量的最常用的仪器旋转风杯风速计是风速测量的最常用的仪器. .它的感它的感应部分是一个固定在旋转轴上的感应部件应部分是一个固定在旋转轴上的感应部件. .一般由三个半一般由三个半球型的空心杯壳组成球型的空心杯壳组成, ,杯壳固定在互成杯壳固定在互成120120度的支架上度的支架上, ,杯杯的凹面顺着一个方向排列的凹面顺着一个方向排列, ,支架固定在垂直的旋转轴上。支架固定在垂直的

23、旋转轴上。 在稳定的风力作用下在稳定的风力作用下, ,风杯受到扭力矩而开始旋转风杯受到扭力矩而开始旋转, ,它的转速与风速成一定关系。它的转速与风速成一定关系。设：风速为设：风速为u，第，第i个风杯与空气相对运动速度个风杯与空气相对运动速度uri2cosriiuunRn：风杯每秒转数：风杯每秒转数; R：旋转半径：旋转半径 i：气流与杯内法线夹角气流与杯内法线夹角. .单位时间气流对风杯作用的有效质量：单位时间气流对风杯作用的有效质量： A：风杯横切面积：风杯横切面积 Cn(i)：风杯为：风杯为i的压力系数的压力系数()niACu作用于互成作用于互成120三个风杯上的风压力为：三个风杯上的风压

24、力为：(2)uiiPA u auRbn000000( )(120 )(240 )( )cos(120 )cos(120 )(240 )cos(240 )inininiiniiniiniiaCCCbCCC（8.35）其中其中 风杯的压力系数风杯的压力系数; 为风杯的阻力系数为风杯的阻力系数; 为空气密度为空气密度.在风杯旋转时在风杯旋转时,它所受到的风压随风杯所处的它所受到的风压随风杯所处的 角度而变化角度而变化,但每转但每转过过120则恢复到则恢复到0时的状态时的状态.如果取如果取0120范围内风压的平范围内风压的平均值均值,(8.35)式就可以简化为式就可以简化为;iaibi取am,bm为

25、和 在 变化120时均值。2(2)ummPAa uRb un在风压下，组件受到的扭力矩：在风压下，组件受到的扭力矩：22(2)2mmMA R a uRb unNuDun（8.36）iaibi（8.37）当风速恒定时当风速恒定时,风杯组件的转速应为某个固定数值。此时风杯组件的转速应为某个固定数值。此时,组件所受组件所受合力为零合力为零,即扭力矩即扭力矩M正好与它的机械系统的动摩擦力矩正好与它的机械系统的动摩擦力矩B1n以及以及静摩擦静摩擦B0相抵消。相抵消。2102BnBNuDnu2012NuBnBDu（8.38）（8.39）当摩擦阻力矩很小可忽略不计时,上式可简化为：2NunD当(8.39)式

26、中n=0时，(非u=0)，则：2min020NuBumin为风杯风速计的启动风速（为风杯风速计的启动风速（风杯风速仪的检定曲线风杯风速仪的检定曲线）。从(8.41)式可看出,起动风速取决于静摩擦力矩B0和扭力矩N。所以,高灵敏度低启动风速的风杯风速计,一般采用低静摩擦的悬浮式。虽然增大扭力矩也能减小 umin 但增大N,势必增大风杯的横臂长度和质量。0min2BuN（8.40）对某一风杯风速计,扭力矩N和空气动力阻尼D一定,风杯的转速与风速成正比关系。（8.41）8.3.3 8.3.3 风杯风速计的惯性风杯风速计的惯性响应过程中风杯风速计的转动方程为：响应过程中风杯风速计的转动方程为：2011

27、22ndnJB nBNuD u nd2101111112( )( )22NuBBDuBDudnndJBDuJ 11nnT其中：111122JJLTBDuDuu2 JLDn1：实际转速：实际转速 ; n(t)：环境实际风速计转速。：环境实际风速计转速。T为风速计的时间为风速计的时间常数常数;L为尺度常数为尺度常数;时间常数越小越好时间常数越小越好,即风速越能及时反应实际即风速越能及时反应实际环境风速环境风速.(8.43)式与温度计的散热方程式与温度计的散热方程 很相似很相似。转动扭力矩转动扭力矩 动摩擦力矩动摩擦力矩 静摩擦力矩静摩擦力矩 扭力矩扭力矩 空气阻尼力矩空气阻尼力矩整理后的:（8.4

28、2）（8.43）1dTdT 8.4 8.4 散热式风速计原理散热式风速计原理一、旁热式热线风速计一、旁热式热线风速计（结构图示）（结构图示）210.24tdQi R20.241()i Rt其中：其中：Rt为加热后金属丝阻值；为加热后金属丝阻值；t为金属丝温度；为金属丝温度； R为为金属丝在环境温度下阻值；：环境温度。 被加热的金属丝在气流中散热，它的散热率与风速的被加热的金属丝在气流中散热，它的散热率与风速的大小有密切的关系，利用这一特性可制成测量风速的仪器。大小有密切的关系，利用这一特性可制成测量风速的仪器。 假若以一定电功率的电流连续加热金属丝，热丝的温度与气温存假若以一定电功率的电流连续

29、加热金属丝，热丝的温度与气温存在一定的差值，差值的大小随风速的大小而异，利用热丝的温差电动在一定的差值，差值的大小随风速的大小而异，利用热丝的温差电动势可测定风速，利用这种原理制成的风速仪称作热线风速仪。热线风势可测定风速，利用这种原理制成的风速仪称作热线风速仪。热线风速仪有速仪有旁热式和直热式之分。旁热式和直热式之分。 焦尔热使具有电阻值的通电流金属丝温度升高，与此同时，金属焦尔热使具有电阻值的通电流金属丝温度升高，与此同时，金属丝将通过它的表面向四周空气散热。电流提供的热功率为：丝将通过它的表面向四周空气散热。电流提供的热功率为：（8.44）在风速在风速v的气流中，一根垂直于气流的金属丝的

30、散热率为：的气流中，一根垂直于气流的金属丝的散热率为：2()()dQAB v t其中：其中：A：分子散热作用，：分子散热作用，Bv1/2：气流作用，：气流作用，t-：热线与气温的差热线与气温的差。当当v 较大时，可忽略分子散热作用项较大时，可忽略分子散热作用项A，且热交换达到平衡时，则：，且热交换达到平衡时，则：2120.24()dQdQi RB v t由空气动力学原理，热线在气流中散热率为：由空气动力学原理，热线在气流中散热率为：（8.45）（8.45）() ()nvDclkt() ()nnnc lkDvt() ()nKvt其中：()nekcRDl tD2()()dQhF th Dl t,n

31、ennkhCRkDvDKC lkDe为空气分子导热系数;R其中：h为散热系数；D为热线直径；l为热线长度；C和n为实验确定的系数（P.115，表5.2);为空气动力学粘性系数。（8.46）当dQ1=dQ2时：20.24() ()ni RKvt由于测温热电偶的电动势正比于冷热端温差，则，上式写为：由于测温热电偶的电动势正比于冷热端温差，则，上式写为：20.24()nti RKv对上式两边取对数：20.24lglg()lg,ti RcnvKRxcnyxycn热 丝 电 阻加 热 电 流 i,分 子导 热 系 数 K及 n为 常 数 ,则 上 式 简 化 为 ;（8.47）（8.48）（8.49）（

32、8.50）在双对数坐标纸上，纵轴为风速取对数在双对数坐标纸上，纵轴为风速取对数lgv，横轴为电动势取对数，横轴为电动势取对数lgv即可得到热线风速仪的风速与输出电动势的线性关系。即可得到热线风速仪的风速与输出电动势的线性关系。在在1-2米米/秒的风速处出现秒的风速处出现折线，折线，是是因为因为雷诺数雷诺数Re在在40左右不连续。左右不连续。 旁热式热线风速计检定曲线图示旁热式热线风速计检定曲线图示旁热式热线风速计检定应用时应注意两点：旁热式热线风速计检定应用时应注意两点： 1.检定时和测量时的空气密度校正（检定时和测量时的空气密度校正（使用前进行检定使用前进行检定） v=v=0 0v vr r

33、 v=v v=vr r0 0/ / 2. 2.方向校正，热线感应的风速方向校正，热线感应的风速u u和环境实际风速和环境实际风速v v的关系为：的关系为：22222(sincos)uvA2222,1sincosuvAA热线越长热线越长A值越小。值越小。（8.50）（8.51）二二.直热式热线微风仪：结构示意图直热式热线微风仪：结构示意图 直热式热线风速仪的感应部分是一根直直热式热线风速仪的感应部分是一根直径约为径约为5-10 m的铂金属丝，紧绷在支架的铂金属丝，紧绷在支架上长度约几个至上长度约几个至20mm。由于较大的电。由于较大的电流流经铂丝，它的温度要比环境空气温流流经铂丝，它的温度要比环

34、境空气温度高度高200-500。 直热式热线风速仪的铂金属丝，一丝两直热式热线风速仪的铂金属丝，一丝两用，它既用来感应风速，又以它的电阻用，它既用来感应风速，又以它的电阻值确定热线的温度，它的感应方程为：值确定热线的温度，它的感应方程为：22110 .2 410 .2 4ntttntttiRKvtRRtiRKvRRKKR（8.52）（8.53）（8.54）（8.55）超声风速仪超声风速仪8.58.5 声学风速仪原理声学风速仪原理超声风速仪：超声风速仪：sonic anemometer声波在大气中传播速度与空气的温度、风速有关。静止空气中的声波在大气中传播速度与空气的温度、风速有关。静止空气中的

35、声速为声速为20.067svcTTsv：声学绝对虚温：声学绝对虚温(10.32)sveTTp式中式中e为水汽压；为水汽压；P为大气压力；气温为大气压力；气温20时，干空气中声速时，干空气中声速343.5m/s. 设设d为超声风速仪的路经长度为超声风速仪的路经长度;v为风速为风速;t1为第一个发声为第一个发声源到达接收端所需的时间源到达接收端所需的时间;t2为第二个发声源到达接收端的为第二个发声源到达接收端的时间时间.只要求出只要求出t1和和t2即可求出风速。即可求出风速。 （8.56）（8.57）smcCTPeTqTTRTCCcttdcttdvvctdvctdvvvvp/5 .34320)32

36、2. 01 ()518. 01 (067.20)()11(2)11(202/1212121时2222 2()()()xyzxv tyv tzv tc t由左图，当逆风发射时，由左图，当逆风发射时，T1发射，发射，R1接收：接收：222221111()(0)(0)xyzdv tv tv tc t2222211()()xyzdv tcvvt2211xdv tcvt假设气流速度为假设气流速度为v的三个分量为的三个分量为vx、vy、vz，声波从坐标原点到达某一，声波从坐标原点到达某一位相面（位相面（x,y,z)所需的时间为所需的时间为t时，则有：时，则有：设设y和和z等于零，等位相面到达点等于零，等位

37、相面到达点(d1,0,0)和和点点(d2,0,0)的时间为的时间为t1和和t2.2212222()xxdcvvdtcvvcv当顺风发射时，当顺风发射时，T2发射，发射，R2接收：接收：222222222()(0)(0)xyzdv tv tv tc t2222222()()xyzdv tcvvt2222xdv tcvt2222222()xxdcvvdtcvcvv（8.58）（8.59）12/ttd c由（由（8.58和和(8.59)式可见，当风速为零时，声波到达两个方向相反，式可见，当风速为零时，声波到达两个方向相反，距离相等的接收器所需时间相等。即距离相等的接收器所需时间相等。即 而有风时，而

38、有风时，v 0,vd 0,t1 t2.于是就利用这一原理，测量相对方向上声波脉冲的传于是就利用这一原理，测量相对方向上声波脉冲的传播时间之差，以测得风速，即：播时间之差，以测得风速，即： （8.60）2122222222121xxxxxdvdvtttcVA cctvdVAVc 22式中 ,当c时,A=12122222222121yyyyydvdvtttcVA cctvdVAVc 22式中 ,当c时,A=1同理同理,y方向的风速分量方向的风速分量vy的计算式为的计算式为:（8.61）yytv计算示意图2122222222121zzzzzdvdvtttcVA cctvdVAVc 22式中 ,当c时

39、,A=1同理同理,z方向的风速分量方向的风速分量vz的计算式为的计算式为:（8.62）2222xyzVvvvzztv计算示意图 由由(8.60)(8.60)可得结论可得结论: :声波沿声程方向顺风传播和逆风传播的声波沿声程方向顺风传播和逆风传播的时间差与沿声程方向的风速成正比时间差与沿声程方向的风速成正比. .由此可知由此可知, ,声学风速仪的测声学风速仪的测速技术就归之于测量声波传播的时间差速技术就归之于测量声波传播的时间差. .用适当的电子线路和用适当的电子线路和装置实现此项时间差的测量装置实现此项时间差的测量, ,便可测得所需风速及其分量便可测得所需风速及其分量. .下面下面是一实用是一

40、实用原理图原理图: : 通常超声风速仪有三组测风探头通常超声风速仪有三组测风探头, ,以测量以测量x,y,zx,y,z三个风速分量。超声风速仪响应速度很快三个风速分量。超声风速仪响应速度很快, ,可每秒测量可每秒测量2020次次, ,测量范围可从测量范围可从0.01-300.01-30米米/ /秒。秒。 8.6 8.6 激光风速仪原理激光风速仪原理一、概述：一、概述： 激光风速仪（激光风速仪（Laser Doppler Anemometer LDA; LaserDoppler Velocimeter LDV) 是建立在激光技术和多普勒频移原理基础上，是建立在激光技术和多普勒频移原理基础上，通过

41、频率测量来测定风速的（通过频率测量来测定风速的（FV）。 激光通过大气层时，大气层中的气溶胶粒子对入射光有散射效应，激光通过大气层时，大气层中的气溶胶粒子对入射光有散射效应，而运行的气溶胶粒子将使散射光的频率产生多普勒频移效应。在接收而运行的气溶胶粒子将使散射光的频率产生多普勒频移效应。在接收器内比较发射的参考光和散射光的频差，就可确定运载气溶胶粒子器内比较发射的参考光和散射光的频差，就可确定运载气溶胶粒子（Aerosol particle)的气流速度。的气流速度。 激光风速仪（激光风速仪（LDA）有两大优点：）有两大优点： 1、可远距离遥测，响应速度快；、可远距离遥测，响应速度快； 2、完全

42、不干扰自然流场。、完全不干扰自然流场。 激光风速仪的机构：激光风速仪的机构： 1、激光器；、激光器；2、入射光学单元；、入射光学单元；3、收集光学单元；、收集光学单元；4、多普勒信号处、多普勒信号处理单元、理单元、5、数据处理系统。、数据处理系统。差动式频移差动式频移LDV系统示意图系统示意图1.激光器：激光器： He-Ne , = 632.8nm= 632.8nm Ar , Ar , =488nm(blue) , =488nm(blue) , =514.5nm(green)=514.5nm(green)2.入射光学单元：入射光学单元： 多普勒频移频率多普勒频移频率2sin2,2sin2DDu

43、fuf3.光学接收单元光学接收单元 接收运动流体中粒子的散射光，得到多普勒频移频率接收运动流体中粒子的散射光，得到多普勒频移频率fD4.多普勒信号处理器多普勒信号处理器 频率跟踪器，计数式处理器，光子相关器频率跟踪器，计数式处理器，光子相关器5.数据处理系统数据处理系统 模拟量、数字量转换成风速模拟量、数字量转换成风速二二.激光测速原理激光测速原理 设照射光的频率为设照射光的频率为f0，粒子，粒子P的运动速度为的运动速度为u，fs为粒子为粒子多勒效应频率多勒效应频率, e0为粒子入射光单位向量，为粒子入射光单位向量，es为粒子散射光单为粒子散射光单位向量，位向量，c为介质中的光速。根据相对论变

44、换公式，经多普勒为介质中的光速。根据相对论变换公式，经多普勒效应后，粒子接收到的光波频率为：效应后，粒子接收到的光波频率为：010201/1(/)uecffuec当当ue0 c 时，可得时，可得(8.63)式的近似式式的近似式粒子作为运动的光源，以此频率向四周发射散射光。当静止粒子作为运动的光源，以此频率向四周发射散射光。当静止的观测者（光检测器）从某一方向上观测粒子的散射光时，的观测者（光检测器）从某一方向上观测粒子的散射光时，由于它们之间存在相对运动，接收器接收到的散射光频率又由于它们之间存在相对运动，接收器接收到的散射光频率又会与粒子所接收到的不同，其频率为：会与粒子所接收到的不同，其频

45、率为：cueff0011（8.63）这就是在静止的光源和运动的粒子条件下，经过一次多普勒这就是在静止的光源和运动的粒子条件下，经过一次多普勒效应后的频率关系式。效应后的频率关系式。（8.64）cueffss11（8.65）es为粒子散射光单位向量，为粒子散射光单位向量，c为介质中的光速。为介质中的光速。 取正是因取正是因为选择为选择es向量由粒子朝向光检测器。将向量由粒子朝向光检测器。将(8.64)代入代入(8.65)式，式，当当u c 时，忽略高次项，可得经两次多普勒效应后的频率时，忽略高次项，可得经两次多普勒效应后的频率关系式为关系式为001() / ssffu eec（8.66）（8.6

46、7）cues最后可得它与光源频率之差（多普勒频移），即最后可得它与光源频率之差（多普勒频移），即00001fceeueeufffsssD 因此，可设计一种光路，能使得接收器同时接收到两种因此，可设计一种光路，能使得接收器同时接收到两种频率的光，由光电接收器检测出差频来。频率的光，由光电接收器检测出差频来。 如下图所示，由两束频率相同的入射光，照射到运动的如下图所示，由两束频率相同的入射光，照射到运动的粒子粒子P上。上。由前面多普勒频移原理：由前面多普勒频移原理：00()/Dsff u eec对第一束入射光，接收器接收到的对第一束入射光，接收器接收到的散射光的频率为：散射光的频率为：1111()

47、/ siisffu eec对第二束光，接收器接收到的散对第二束光，接收器接收到的散射光的频率为：射光的频率为：2221()/ siisffu eec这两部分散射光在接收器上混频后，可测得差频为：这两部分散射光在接收器上混频后，可测得差频为：111112()(),isiisDssiiiu eefu eecfffcf接收器接收器（8.68）当当es1 =es2，多普勒频移，多普勒频移fD是两部分光波频率之差。由两束是两部分光波频率之差。由两束入射光的单位矢量入射光的单位矢量ei1 , ei2与风速与风速u的点乘积关系：的点乘积关系：11()2sin2isu eeu2sin2,2 sin2iDDiu

48、fuf 是两束入射光的夹角，将是两束入射光的夹角，将（8.69）代入（）代入（8.68）得：）得：（8.69）（8.70）由（由（8.70）可见，当可见，当 和和 一定时，粒子运动速度与一定时，粒子运动速度与多普勒频多普勒频移移fD成正比。就是说，可以通过测量多普勒频移成正比。就是说，可以通过测量多普勒频移fD而得到而得到风速风速u，这就是激光测量风速的原理。，这就是激光测量风速的原理。三、激光测速仪存在的问题：三、激光测速仪存在的问题： 激光风速仪测得的风速并不是真正的大气流动速度，而是激光风速仪测得的风速并不是真正的大气流动速度，而是悬浮于气流中的散射粒子的运动速度。因为粒子与流体的密度悬

49、浮于气流中的散射粒子的运动速度。因为粒子与流体的密度有显著差异，粒子不可能完全跟随流体运动，其速度总会有差异，有显著差异，粒子不可能完全跟随流体运动，其速度总会有差异，这就造成了激光测量的误差。这就造成了激光测量的误差。100100(1),100100fppffuuuuupuf若以为粒子速度,u 为流体速度,则由于粒子不完全跟随流体运动而造成的相对误差为;0,0,0,0,这就是激光测速技术中的粒子跟随问题.=1, =1,完全跟随;粒子超前于流体运动;粒子滞后流体运动. 一、旋臂机一、旋臂机 图示图示 旋臂旋转时，与空气产生相对运动（空气假定静止），旋臂旋转时，与空气产生相对运动（空气假定静止）

50、，其相对运动速度为：其相对运动速度为： 式中式中n n为旋臂转动圈数；为旋臂转动圈数；R R为旋臂的半径；为旋臂的半径；t t为旋臂转动为旋臂转动n n圈圈时的时间。时的时间。 同时记录风速仪的输出，如热线风速仪的输出电动势同时记录风速仪的输出，如热线风速仪的输出电动势做出做出V V 的关系图或回归方程。的关系图或回归方程。 8.7 8.7 风速仪的检定风速仪的检定2 RnVttt 优点：转速稳定可调；可检定低风速，例如优点：转速稳定可调；可检定低风速，例如0.2m/s;0.2m/s; 缺点：低风速时保证室内静风（门窗不能漏风，缺点：低风速时保证室内静风（门窗不能漏风， 人不能走动）；人不能走

51、动）； 高风速时要保证旋臂卷起的风要消除（旋高风速时要保证旋臂卷起的风要消除（旋 臂周围树立栅栏）以消除大转速时，室内臂周围树立栅栏）以消除大转速时，室内 气流造成的检定误差。气流造成的检定误差。 做出纵坐标为做出纵坐标为V V（m/s)m/s)与横坐标为输出电动势（与横坐标为输出电动势（mv)mv)的的 检定曲线，并作出回归方程。检定曲线，并作出回归方程。 检定曲线图示检定曲线图示二、风洞实验装置环境学院大气边界层风洞，宽环境学院大气边界层风洞，宽3m，高，高2m，长，长32m1 1、风洞的结构、风洞的结构 收缩段、工作段、马达、电力部分收缩段、工作段、马达、电力部分 风洞风速的大小，由电动

52、机的转速决定，风扇驱动风洞风速的大小，由电动机的转速决定，风扇驱动电机有两个绕组，一组绕在旋转铁心上，成为转子；另一电机有两个绕组，一组绕在旋转铁心上，成为转子；另一组绕在外围铁壳上，称为定子；当定子被一个恒定电压供组绕在外围铁壳上，称为定子；当定子被一个恒定电压供电，改变与转子间的磁场，即可改变马达的转速，达到调电，改变与转子间的磁场，即可改变马达的转速，达到调节风洞风速的目的。节风洞风速的目的。 风洞中风速的测量是借助风洞收缩段两个静压口的压风洞中风速的测量是借助风洞收缩段两个静压口的压力差来测量。力差来测量。 设收缩段静压为设收缩段静压为P P1 1，动压为，动压为 ；工作段静压为；工作

53、段静压为P P2 2，动，动压为压为 ；得：；得：212V222V 221221222222222222111222,1112221111212211,0.45PVPVVCV CVVPPVCPPVVCC 为收缩比;为工作段风速.式中.12TPP2a只要测定P1和P2,再测量当时的大气压力和温度,计算出大气密度P=,以及确定 值,就可由上式计算出工作段的风速VR的测量由微压计测量.（1） 微压计的结构如图所示：微压计的结构如图所示： 由（由（1 1）和（）和（2 2）式即可计算出风洞工作段的风速）式即可计算出风洞工作段的风速V V2 2，若，若在风洞工作段中放入风速仪，即可对风速进行标定，即：在

54、风洞工作段中放入风速仪，即可对风速进行标定，即： sin ,0.4,0.2,0.10.05;,Lgg003压力计具有一定的斜率,共分和四挡,视风速大小选择斜率.设风速为零时,P1-P2=0,液柱读数为L当液柱读数为 时,测得静压差为:P1-P2= L-L式中 为酒精密度, =800kg/m为重力加速度.（2）022 LLgV 做出纵坐标为做出纵坐标为V V（m/s)m/s)与横坐标为输出电与横坐标为输出电动势（动势（mv)mv)的检定曲线或者是风速仪转动圈数的检定曲线或者是风速仪转动圈数/ /秒，秒，并作出回归方程。并作出回归方程。 检定曲线图示检定曲线图示 8.8 气象铁塔探测系统气象铁塔探

55、测系统 气象铁塔的架设和仪器配置气象铁塔的架设和仪器配置 梯度测量梯度测量 温度测量温度测量 湿度测量湿度测量 风速测量风速测量 风向测量风向测量气象塔架设的基本要求气象塔架设的基本要求 建于平坦、开阔的地面上建于平坦、开阔的地面上 周围任何障碍物与塔体的距离至少为障碍物本周围任何障碍物与塔体的距离至少为障碍物本身高度的十倍以上身高度的十倍以上 塔与观测机房之间应有一电缆沟相连，以便铺塔与观测机房之间应有一电缆沟相连，以便铺设电源电缆和信号传输总线设电源电缆和信号传输总线 塔体安装后应坚固、可靠，保证十年以上的正塔体安装后应坚固、可靠，保证十年以上的正常使用常使用 塔体要能承受当地最恶劣气象条

56、件，如大风、塔体要能承受当地最恶劣气象条件，如大风、沙尘暴、暴雨等沙尘暴、暴雨等 塔体要符合避雷安全要求塔体要符合避雷安全要求气象铁塔的架设和仪器配置气象铁塔的架设和仪器配置 塔结构：三角形结构，空心，塔结构：三角形结构，空心， 主要指标：塔高和边长主要指标：塔高和边长 层次设计：对数或指数分布层次设计：对数或指数分布 仪器架设臂长：与直径成比例仪器架设臂长：与直径成比例 (3倍倍) 钢制底座和混凝土地基，钢丝拉线加固钢制底座和混凝土地基，钢丝拉线加固 便于工作人员拆卸、维修仪器设备便于工作人员拆卸、维修仪器设备气象观测塔观测系统硬件部分各主要单元结构气象观测塔观测系统硬件部分各主要单元结构

57、配套地面观测系统硬件部分各主要单元结构配套地面观测系统硬件部分各主要单元结构塔上仪器配置塔上仪器配置 每层架设风向标、风杯风速计、铂金属温每层架设风向标、风杯风速计、铂金属温度计、湿度计，进行慢响应平均场梯度测度计、湿度计，进行慢响应平均场梯度测量，每量，每20秒记录一组数据。秒记录一组数据。 在需要的高度上架设超声风速温度仪，进在需要的高度上架设超声风速温度仪，进行湍流场测量，每秒采集行湍流场测量，每秒采集10次或次或20次，得次，得到三维风速和温度脉动数据。到三维风速和温度脉动数据。 根据需要架设各种仪器设备。根据需要架设各种仪器设备。梯度测量梯度测量 在不同高度上进行：在不同高度上进行：

58、 温度测量：铂电阻温度计温度测量：铂电阻温度计 湿度测量：湿敏电容湿度计湿度测量：湿敏电容湿度计 风速测量：旋转式风速计风速测量：旋转式风速计 风向测量：风向标风向测量：风向标温度、湿度传感器和防辐射罩温度、湿度传感器和防辐射罩大气湍流测量大气湍流测量 1、三维湍流脉动风速测量、三维湍流脉动风速测量 超声风速计超声风速计 2、温度脉动测量、温度脉动测量 超声风速计测温度脉动超声风速计测温度脉动 铂金属温度脉动仪铂金属温度脉动仪 3、湿度脉动测量、湿度脉动测量 紫外湿度计紫外湿度计 红外湿度计红外湿度计 超声风速仪超声风速仪超声风速仪性能参数超声风速仪性能参数 声程：声程：15 cm（或（或20

59、）（）（X、Y、Z三轴向）三轴向） 量程：风速：量程：风速：0.01 - 30 m/s（或（或60） 风向：风向：0 - 360 温度：温度：-30 - 50 采样速率：最大每秒采样速率：最大每秒20次次 (20Hz) 传输速率：传输速率：400 kb/s 传输距离：传输距离：400 m超声风速仪的安装超声风速仪的安装 底座保持水平，防止水平风对垂直方向测底座保持水平，防止水平风对垂直方向测量的影响量的影响 阴影效应和减少塔体影响阴影效应和减少塔体影响 取决于探头轴线与风向夹角和风速取决于探头轴线与风向夹角和风速 塔体迎风时会产生较大影响塔体迎风时会产生较大影响 方向性和角度修正方向性和角度修正 磁罗盘法或北极星法定方位磁罗盘法或北极星法定方位超声风速仪的资料采集超声风速仪的资料采集 记录数据：记录数据：U、V、W、T （T(Pt） 生成数据：三维平均风速生成数据：三维平均风速 平均风向平均风向 温度温度 三维均方差三维均方差 平均动量通量、热量通量平均动量通量、热量通量2*wuuNiiuVuNu12)(1122*uuwM* TuCtwCHpp三分量风速仪三分量风速仪湿度脉动测量湿度脉动测量 紫外湿度计紫外湿度计 利用水汽拉曼利用水汽拉曼-阿尔法吸收线（波长阿尔法吸收线