CN107870332A - 一种风廓线雷达大气湍流能量探测系统及方法 - Google Patents

一种风廓线雷达大气湍流能量探测系统及方法 Download PDF

Info

Publication number
CN107870332A
CN107870332A CN201610851730.1A CN201610851730A CN107870332A CN 107870332 A CN107870332 A CN 107870332A CN 201610851730 A CN201610851730 A CN 201610851730A CN 107870332 A CN107870332 A CN 107870332A
Authority
CN
China
Prior art keywords
data
atmospheric turbulance
component
energy
component data
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.)
Granted
Application number
CN201610851730.1A
Other languages
English (en)
Other versions
CN107870332B (zh
Inventor
任迎新
沈江林
王志锐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Radio Measurement
Original Assignee
Beijing Institute of Radio Measurement
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Radio Measurement filed Critical Beijing Institute of Radio Measurement
Priority to CN201610851730.1A priority Critical patent/CN107870332B/zh
Publication of CN107870332A publication Critical patent/CN107870332A/zh
Application granted granted Critical
Publication of CN107870332B publication Critical patent/CN107870332B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及大气湍流能量探测领域，尤其涉及一种风廓线雷达大气湍流能量探测系统及方法。该系统包括数字中频接收机，用于采集大气湍流雷达回波信号，将雷达回波信号进行数字化处理得到时域数据并发送到第一存储模块中；第一存储模块，用于存储数字中频接收机发送的时域数据；信号处理模块，用于对第一存储模块中的时域数据进行信号处理，生成大气湍流回波功率谱数据并发送到第二存储模块中。本发明提供了一种用于实现对大气湍流能量和摩擦速度进行观察、探测和处理的风廓线雷达大气湍流能量探测系统和方法，能实时获得数据，数据准确，精确度高，真实还原大气湍流，可应用于科研、大气观测和农业应用等多个领域。

Description

一种风廓线雷达大气湍流能量探测系统及方法
技术领域
本发明涉及大气湍流能量探测领域，尤其涉及一种风廓线雷达大气湍流能量探测系统及方法。
背景技术
大气湍流是大气中一种不规则的随机运动，湍流每一点上的压强、速度、温度等物理特性等随机涨落。大气湍流的存在使大气中的动量、热量、水气和污染物的垂直和水平交换作用明显增强，远大于分子运动的交换强度，同时对光波、声波和电磁波在大气中的传播产生一定的干扰作用。大气湍流最常发生的3个区域是：大气底层的边界层内，对流云的云体内部，大气对流层上部的西风急流区内。
大气湍流对飞行器的飞行性能、结构载荷、飞行安全的影响很大；大气边界层中，是地气相互作用的重要过程,同时也对农业、森林和大气污染等实际应用问题起至关重要的作用；另外,近地面层湍流通量的确定还是大气环流和数值模拟等相关科研工作的重要问题。对于大气湍流的观测和研究一直是有关专业领域的重要内容，但是，当前湍流理论尚不完备，湍流观测仪器和实验手段存在局限性。
在气象观测中，用风廓线雷达探测大气湍流速度已经得到普遍应用，能够探测大气边界层、对流层等区域，得到大气湍流风场的水平风速、水平风向和垂直气流速度。风廓线雷达具有时空分辨力高、精度高等优点，是全天候无人值守探测方式。随着风廓线雷达的推广应用，获得了大量观测数据，风廓线雷达的回波功率谱中含有丰富的大气湍流信息，在此基础上，能够更加深入和精细化地研究大气湍流；其中，大气湍流能量TKE(TurbulenceKinetic Energy)和摩擦速度FV(Friction Velocity)是大气湍流研究中的重要参量，表征了湍流尺度、能量传输、随机性运动规律等大气湍流的特征，湍流能量很大时，随机性的湍流可能形成风切变，研究湍流与风切变的关系对预警风切变有一定的意义。
受实际探测技术的限制，以前对于大气湍流的这两个参量主要是在理论、实验和数值模拟研究方面有一些成果，而原有风廓线雷达的数据处理和数据产品中，还没有对大气湍流能量(TKE)和摩擦速度(FV)进行探测和处理的先例。
发明内容
本发明所要解决的技术问题是针对现有技术的不足，提供一种风廓线雷达大气湍流能量探测系统及方法。
本发明解决上述技术问题的技术方案如下：
一种风廓线雷达大气湍流能量探测系统，包括：
数字中频接收机，用于采集大气湍流雷达回波信号，将雷达回波信号进行数字化处理得到时域数据并发送到第一存储模块中；
第一存储模块，用于存储数字中频接收机发送的时域数据；
信号处理模块，用于对第一存储模块中的时域数据进行信号处理，生成大气湍流回波功率谱数据并发送到第二存储模块中；
第二存储模块，用于存储信号处理模块发送的大气湍流回波功率谱数据；
数据处理模块，用于对第二存储模块中的大气湍流回波功率谱数据进行谱矩处理，得到多个大气风场的分量数据，将多个分量数据按获取的时间顺序进行有重叠分组，按分组的顺序处理分量数据，得到按时间顺序排列的大气湍流能量和摩擦速度并储存。
本发明的有益效果：数字中频接收机采集信号并处理得到时域数据，第一存储模块存储时域数据，方便查询原始记录，信号处理模块处理时域数据得到功率谱数据，第二存储模块存储功率谱数据，数据处理模块对功率谱数据进行谱矩处理，得到按时间顺序排列的大气湍流能量和摩擦速度，提供了一种用于实现对大气湍流能量(TKE)和摩擦速度(FV)进行探测和处理的风廓线雷达大气湍流能量探测系统，该系统实时获得数据，数据准确，精确度高，真实还原大气湍流，可应用于科研、大气观测和农业应用等多个领域。
进一步的，所述数据处理模块将多个分量数据按获取的时间顺序进行有重叠分组，具体为，假设组内分量数据的个数为N，则分量数据重叠的个数为N-1，其中重叠部分为连续获取的分量数据，每得到一个分量数据，与所述分量数据之前的N-1个分量数据分为一组。
采用上述进一步方案的有益效果是：数据处理模块通过将多个分量数据按获取的时间顺序进行有重叠分组，充分利用数字中频接收机采集的信号，多方面多角度探测、处理得到按时间顺序排列的大气湍流能量和摩擦速度，风廓线雷达大气湍流能量探测系统实时获得数据，数据准确，精确度高，真实还原大气湍流。
进一步，所述组内分量数据的个数为5，得到大气风场的第一组分量数据，处理得到一个大气湍流能量和摩擦速度的数值，按获得分量数据的时间顺序再获得一个分量数据，所述分量数据与所述分量数据之前的4个分量数据分为另一组，再处理得到一个大气湍流能量和摩擦速度的数值，依照分组的顺序处理分量数据得到按时间顺序排列的大气湍流能量和摩擦速度并储存。
采用上述进一步方案的有益效果是：数据处理模块通过5个分量数据按获取的时间顺序处理得到一个大气湍流能量和摩擦速度的数值，按获得分量数据的时间顺序再获得一个分量数据，与所述分量数据之前的4个分量数据再处理得到一个大气湍流能量和摩擦速度的数值，依照分组的顺序处理分量数据，得到按时间顺序排列的大气湍流能量和摩擦速度，多方面多角度探测、处理得到按时间顺序排列的大气湍流能量和摩擦速度，风廓线雷达大气湍流能量探测系统实时获得数据，数据准确，精确度高，真实还原大气湍流。
进一步，所述分量数据包括u分量、v分量、w分量，所述u分量、v分量、w分量是三个相互的垂直方向上的径向速度分量。
进一步，还包括系统控制模块，所述系统控制模块用于在开始探测前对所述信号处理模块进行加载参数处理，对所述数字中频接收机进行初始化和加载参数处理。
采用上述进一步方案的有益效果是：系统控制模块分别对信号处理模块和数字中频接收机进行相应处理，提供了一种用于实现对大气湍流能量(TKE)和摩擦速度(FV)进行观察、探测和处理的风廓线雷达大气湍流能量探测系统。
进一步，还包括显示模块，所述显示模块用于显示数据处理模块发送的大气湍流能量和摩擦速度。
采用上述进一步方案的有益效果是：显示模块显示大气湍流能量和摩擦速度，提供了一种用于实现对大气湍流能量(TKE)和摩擦速度(FV)进行观察的风廓线雷达大气湍流能量探测系统。
进一步，所述数字中频接收机、第一存储模块、信号处理模块和系统控制模块设置在工业控制计算机内；第二存储模块、数据处理模块和显示模块设置在数据终端计算机内；所述工业控制计算机内各个部分通过外设部件互连总线进行通讯，数据终端计算机内各个部分通过外设部件互连总线进行通讯，工业控制计算机和数据终端计算机通过以太网卡和网线进行通讯。
采用上述进一步方案的有益效果是：数字中频接收机、第一存储模块、信号处理模块和系统控制模块设置在工业控制计算机内；第二存储模块、数据处理模块和显示模块设置在数据终端计算机内。各个部分实现通信，提供了一种用于实现对大气湍流能量(TKE)和摩擦速度(FV)进行观察、探测和处理的风廓线雷达大气湍流能量探测系统。
本发明还提供一种风廓线雷达大气湍流能量探测方法，包括如下步骤：
S1,采集大气湍流雷达回波信号，将雷达回波信号进行数字化处理得到时域数据；
S2,存储时域数据；
S3,对时域数据进行信号处理，生成大气湍流回波功率谱数据；
S4,存储大气湍流回波功率谱数据；
S5,对大气湍流回波功率谱数据进行谱矩处理，得到多个大气风场的分量数据，将多个分量数据按获取的时间顺序进行有重叠分组，按分组的顺序处理分量数据，得到按时间顺序排列的大气湍流能量和摩擦速度并储存。
本发明的有益效果：采集信号并处理得到时域数据，存储时域数据，方便查询原始记录，处理时域数据得到功率谱数据，存储功率谱数据，对功率谱数据进行谱矩处理，得到按时间顺序排列的大气湍流能量和摩擦速度，提供了一种用于实现对大气湍流能量(TKE)和摩擦速度(FV)进行探测和处理的风廓线雷达大气湍流能量探测方法，该方法实时获得数据，数据准确，精确度高，真实还原大气湍流，可应用于科研、大气观测和农业应用等多个领域。
进一步，所述S5中将多个分量数据按获取的时间顺序进行有重叠分组，具体为，假设组内分量数据的个数为N，则分量数据重叠的个数为N-1，其中重叠部分为连续获取的分量数据，每得到一个分量数据，与所述分量数据之前的N-1个分量数据分为一组。
采用上述进一步方案的有益效果是：通过将多个分量数据按获取的时间顺序进行有重叠分组，充分利用采集的信号，多方面多角度探测、处理得到按时间顺序排列的大气湍流能量和摩擦速度，风廓线雷达大气湍流能量探测方法实时获得数据，数据准确，精确度高，真实还原大气湍流。
进一步，所述组内分量数据的个数为5，得到大气风场的第一组分量数据，处理得到一个大气湍流能量和摩擦速度的数值，按获得分量数据的时间顺序再获得一个分量数据，所述分量数据与所述分量数据之前的4个分量数据分为另一组，再处理得到一个大气湍流能量和摩擦速度的数值，依照分组的顺序处理分量数据得到按时间顺序排列的大气湍流能量和摩擦速度并储存。
采用上述进一步方案的有益效果是：通过5个分量数据按获取的时间顺序处理得到一个大气湍流能量和摩擦速度的数值，按获得分量数据的时间顺序再获得一个分量数据，与所述分量数据之前的4个分量数据再处理得到一个大气湍流能量和摩擦速度的数值，依照分组的顺序处理分量数据，得到按时间顺序排列的大气湍流能量和摩擦速度，多方面多角度探测、处理得到按时间顺序排列的大气湍流能量和摩擦速度，风廓线雷达大气湍流能量探测方法实时获得数据，数据准确，精确度高，真实还原大气湍流。
附图说明
图1为本发明实施例所述一种风廓线雷达大气湍流能量探测系统的结构示意图；
图2为本发明实施例所述一种风廓线雷达大气湍流能量探测方法的流程示意图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。
本发明实施例提供一种风廓线雷达大气湍流能量探测系统。
如图1所示，一种风廓线雷达大气湍流能量探测系统，包括：
数字中频接收机，用于采集大气湍流雷达回波信号，将雷达回波信号进行数字化处理得到时域数据并发送到第一存储模块中；
具体地，所述数字中频接收机按照加载参数对雷达回波的中频信号进行模数转换和采集，将得到的时域数据存入本地先进先出存储器并发送到第一存储模块中。
第一存储模块，用于存储数字中频接收机发送的时域数据；
信号处理模块，用于对第一存储模块中的时域数据进行信号处理，生成大气湍流回波功率谱数据并发送到第二存储模块中；
第二存储模块，用于存储信号处理模块发送的大气湍流回波功率谱数据；
数据处理模块，用于对第二存储模块中的大气湍流回波功率谱数据进行谱矩处理，得到多个大气风场的分量数据，将多个分量数据按获取的时间顺序进行有重叠分组，按分组的顺序处理分量数据，得到按时间顺序排列的大气湍流能量和摩擦速度并储存。
数字中频接收机采集信号并处理得到时域数据，第一存储模块存储时域数据，方便查询原始记录，信号处理模块处理时域数据得到功率谱数据，第二存储模块存储功率谱数据，数据处理模块对功率谱数据进行谱矩处理，得到按时间顺序排列的大气湍流能量和摩擦速度，提供了一种用于实现对大气湍流能量(TKE)和摩擦速度(FV)进行探测和处理的风廓线雷达大气湍流能量探测系统，该系统实时获得数据，数据准确，精确度高，真实还原大气湍流，可应用于科研、大气观测和农业应用等多个领域。
所述数据处理模块将多个分量数据按获取的时间顺序进行有重叠分组，具体为，假设组内分量数据的个数为N，则分量数据重叠的个数为N-1，其中重叠部分为连续获取的分量数据，每得到一个分量数据，与所述分量数据之前的N-1个分量数据分为一组。
数据处理模块通过将多个分量数据按获取的时间顺序进行有重叠分组，充分利用数字中频接收机采集的信号，多方面多角度探测、处理得到按时间顺序排列的大气湍流能量和摩擦速度，风廓线雷达大气湍流能量探测系统实时获得数据，数据准确，精确度高，真实还原大气湍流。
所述组内分量数据的个数为5，得到大气风场的第一组分量数据，处理得到一个大气湍流能量和摩擦速度的数值，按获得分量数据的时间顺序再获得一个分量数据，所述分量数据与所述分量数据之前的4个分量数据分为另一组，再处理得到一个大气湍流能量和摩擦速度的数值，依照分组的顺序处理分量数据得到按时间顺序排列的大气湍流能量和摩擦速度并储存。
具体地，所述数据处理模块对大气湍流回波功率谱数据进行谱矩处理，处理出探测波束所指方向上各高度层的径向风速度，再由径向风速度和探测波束的方向角、俯仰角计算出水平风速在两个正交水平方向上的分量，其中东方向上为u分量，北方向上为v分量，在垂直天顶方向上的径向速度是直接测量得到的，为w分量。5个分量数据为一组，5个分量数据为不同高度的分量数据，高度分辨率最小为50米，一个分量数据包括u分量、v分量、w分量。按获得分量数据的时间顺序依次得到5个分量数据，时间分辨率最小为3分钟，每个分量数据分别包括u分量、v分量、w分量，为大气风场的第一组分量数据，处理得到一个大气湍流能量和摩擦速度的数值。一组分量数据，即5个分量数据得到一个大气湍流能量和摩擦速度的数值。按获得分量数据的时间顺序再获得一个分量数据R，所述分量数据R与所述分量数据R之前相邻的4个分量数据为第二组分量数据，第二组分量数据处理得到一个大气湍流能量和摩擦速度的数值，按获得分量数据的时间顺序又获得一个分量数据Y，分量数据Y与分量数据R、分量数据R之前相邻的三个分量数据为第三组分量数据，第三组分量数据处理得到一个大气湍流能量和摩擦速度的数值，依照分组的顺序处理分量数据得到按时间顺序排列的大气湍流能量和摩擦速度。
u、v、w的计算公式，单位米/秒(m/s)
w(h)＝VRZ(h) (3)
其中，θ为倾斜探测波束的天顶角，即倾斜波束与天顶垂直方向的夹角。
VRZ(h)、VRE(h)、VRN(h)为风廓线雷达在天顶方向、偏东方向、偏北方向的带有高度值的测得的径向速度，u(h)和v(h)分别为水平风在东和北方向的分量，w(h)为大气垂直运动速度。
风廓线雷达测得的径向速度均以朝向雷达方向为正速度。
大气湍流能量TKE的计算公式如下，单位平方米/平方秒(m2/s2)
TKE＝(u'2+v'2+w'2)/2 (4)
其中，u'、v'、w'是分别由u、v、w计算得到的，计算公式如下，
其中：N＝5，ui、vi、wi分别表示每组中的每个u、v、w，分别表示每组u、v、w的平均值。每个u、v、w是同一个高度上的径向风速度分量，每个组中的每个分量数据可以是不同高度的分量数据，但是一定是按照获取分量数据的时间顺序获得的。
摩擦速度FV的计算公式，单位米/秒(m/s)
FV2＝[(duw)2+(dvw)2]1/2 (13)
FV＝[FV2]1/2 (14)
其中：N＝5，ui、vi、wi分别表示每组中的每个u、v、w，分别表示每组u、v、w的平均值，计算公式同(8)、(9)、(10)。
风廓线雷达(Wind Profiling Radar,WPR)是一种多普勒脉冲雷达，是通过探测大气湍流对电磁波的反射和散射来连续探测大气风场等物理量的一种全天候大气遥感探测设备，具有很高的探测实时性和精度，是中小尺度大气探测的关键设备之一，在我国得到了快速发展和应用。
本发明根据大气湍流能量(TKE)和摩擦速度(FV)参量的计算公式，对风廓线雷达探测获取到的丰富的大气湍流风场数据，作进一步的处理和计算，实时得到高时空分辨力和高精度的两个参量数据产品，用于航空气象保障和气候变化等领域的研究。
数据处理模块通过5个分量数据按获取的时间顺序处理得到一个大气湍流能量和摩擦速度的数值，按获得分量数据的时间顺序再获得一个分量数据，与所述分量数据之前的4个分量数据再处理得到一个大气湍流能量和摩擦速度的数值，依照分组的顺序处理分量数据得到按时间顺序排列的大气湍流能量和摩擦速度，多方面多角度探测、处理得到按时间顺序排列的大气湍流能量和摩擦速度，风廓线雷达大气湍流能量探测系统实时获得数据，数据准确，精确度高，真实还原大气湍流。
所述分量数据包括u分量、v分量、w分量，所述u分量、v分量、w分量是三个相互的垂直方向上的径向速度分量。
一种风廓线雷达大气湍流能量探测系统还包括系统控制模块，所述系统控制模块用于在开始探测前对所述信号处理模块进行加载参数处理，对所述数字中频接收机进行初始化和加载参数处理。
系统控制模块分别对信号处理模块和数字中频接收机进行相应处理，提供了一种用于实现对大气湍流能量(TKE)和摩擦速度(FV)进行观察、探测和处理的风廓线雷达大气湍流能量探测系统。
一种风廓线雷达大气湍流能量探测系统还包括显示模块，所述显示模块用于显示数据处理模块发送的大气湍流能量和摩擦速度。
显示模块显示大气湍流能量和摩擦速度，提供了一种用于实现对大气湍流能量(TKE)和摩擦速度(FV)进行观察的风廓线雷达大气湍流能量探测系统。
所述数字中频接收机、第一存储模块、信号处理模块和系统控制模块设置在工业控制计算机内；第二存储模块、数据处理模块和显示模块设置在数据终端计算机内；所述工业控制计算机内各个部分通过外设部件互连总线进行通讯，数据终端计算机内各个部分通过外设部件互连总线进行通讯，工业控制计算机和数据终端计算机通过以太网卡和网线进行通讯。
数字中频接收机、第一存储模块、信号处理模块和系统控制模块设置在工业控制计算机内；第二存储模块、数据处理模块和显示模块设置在数据终端计算机内。各个部分实现通信，提供了一种用于实现对大气湍流能量(TKE)和摩擦速度(FV)进行观察、探测和处理的风廓线雷达大气湍流能量探测系统。
本发明实施例还提供一种风廓线雷达大气湍流能量探测方法。
如图2所示，一种风廓线雷达大气湍流能量探测方法，包括如下步骤：
S1,采集大气湍流雷达回波信号，将雷达回波信号进行数字化处理得到时域数据；
S2,存储时域数据；
S3,对时域数据进行信号处理，生成大气湍流回波功率谱数据；
S4,存储大气湍流回波功率谱数据；
S5,对大气湍流回波功率谱数据进行谱矩处理，得到多个大气风场的分量数据，将多个分量数据按获取的时间顺序进行有重叠分组，按分组的顺序处理分量数据，得到按时间顺序排列的大气湍流能量和摩擦速度并储存。
采集信号并处理得到时域数据，存储时域数据，方便查询原始记录，处理时域数据得到功率谱数据，存储功率谱数据，对功率谱数据进行谱矩处理，得到按时间顺序排列的大气湍流能量和摩擦速度，提供了一种用于实现对大气湍流能量(TKE)和摩擦速度(FV)进行探测和处理的风廓线雷达大气湍流能量探测方法，该方法实时获得数据，数据准确，精确度高，真实还原大气湍流，可应用于科研、大气观测和农业应用等多个领域。
所述S5中将多个分量数据按获取的时间顺序进行有重叠分组，具体为，假设组内分量数据的个数为N，则分量数据重叠的个数为N-1，其中重叠部分为连续获取的分量数据，每得到一个分量数据，与所述分量数据之前的N-1个分量数据分为一组。
通过将多个分量数据按获取的时间顺序进行有重叠分组，充分利用采集的信号，多方面多角度探测、处理得到按时间顺序排列的大气湍流能量和摩擦速度，风廓线雷达大气湍流能量探测方法实时获得数据，数据准确，精确度高，真实还原大气湍流。
所述组内分量数据的个数为5，得到大气风场的第一组分量数据，处理得到一个大气湍流能量和摩擦速度的数值，按获得分量数据的时间顺序再获得一个分量数据，所述分量数据与所述分量数据之前的4个分量数据分为另一组，再处理得到一个大气湍流能量和摩擦速度的数值，依照分组的顺序处理分量数据得到按时间顺序排列的大气湍流能量和摩擦速度并储存。
通过5个分量数据按获取的时间顺序处理得到一个大气湍流能量和摩擦速度的数值，按获得分量数据的时间顺序再获得一个分量数据，与所述分量数据之前的4个分量数据再处理得到一个大气湍流能量和摩擦速度的数值，依照分组的顺序处理分量数据，得到按时间顺序排列的大气湍流能量和摩擦速度，多方面多角度探测、处理得到按时间顺序排列的大气湍流能量和摩擦速度，风廓线雷达大气湍流能量探测方法实时获得数据，数据准确，精确度高，真实还原大气湍流。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种风廓线雷达大气湍流能量探测系统，其特征在于，包括：
数字中频接收机，用于采集大气湍流雷达回波信号，将雷达回波信号进行数字化处理得到时域数据并发送到第一存储模块中；
第一存储模块，用于存储数字中频接收机发送的时域数据；
信号处理模块，用于对第一存储模块中的时域数据进行信号处理，生成大气湍流回波功率谱数据并发送到第二存储模块中；
第二存储模块，用于存储信号处理模块发送的大气湍流回波功率谱数据；
数据处理模块，用于对第二存储模块中的大气湍流回波功率谱数据进行谱矩处理，得到多个大气风场的分量数据，将多个分量数据按获取的时间顺序进行有重叠分组，按分组的顺序处理分量数据，得到按时间顺序排列的大气湍流能量和摩擦速度并储存。
2.根据权利要求1所述的风廓线雷达大气湍流能量探测系统，其特征在于，所述数据处理模块将多个分量数据按获取的时间顺序进行有重叠分组，具体为，假设组内分量数据的个数为N，则分量数据重叠的个数为N-1，其中重叠部分为连续获取的分量数据，每得到一个分量数据，与所述分量数据之前的N-1个分量数据分为一组。
3.根据权利要求2所述的风廓线雷达大气湍流能量探测系统，其特征在于，所述组内分量数据的个数为5，得到大气风场的第一组分量数据，处理得到一个大气湍流能量和摩擦速度的数值，按获得分量数据的时间顺序再获得一个分量数据，所述分量数据与所述分量数据之前的4个分量数据分为另一组，再处理得到一个大气湍流能量和摩擦速度的数值，依照分组的顺序处理分量数据得到按时间顺序排列的大气湍流能量和摩擦速度并储存。
4.根据权利要求1-3任意一项所述的风廓线雷达大气湍流能量探测系统，其特征在于，所述分量数据包括u分量、v分量、w分量，所述u分量、v分量、w分量是三个相互的垂直方向上的径向速度分量。
5.根据权利要求1-3任意一项所述的风廓线雷达大气湍流能量探测系统，其特征在于，还包括系统控制模块，所述系统控制模块用于在开始探测前对所述信号处理模块进行加载参数处理，对所述数字中频接收机进行初始化和加载参数处理。
6.根据权利要求5所述风廓线雷达大气湍流能量探测系统，其特征在于，还包括显示模块，所述显示模块用于显示数据处理模块发送的大气湍流能量和摩擦速度。
7.根据权利要求6所述风廓线雷达大气湍流能量探测系统，其特征在于，所述数字中频接收机、第一存储模块、信号处理模块和系统控制模块设置在工业控制计算机内；第二存储模块、数据处理模块和显示模块设置在数据终端计算机内；所述工业控制计算机内各个部分通过外设部件互连总线进行通讯，数据终端计算机内各个部分通过外设部件互连总线进行通讯，工业控制计算机和数据终端计算机通过以太网卡和网线进行通讯。
8.一种风廓线雷达大气湍流能量探测方法，其特征在于，包括如下步骤：
S1,采集大气湍流雷达回波信号，将雷达回波信号进行数字化处理得到时域数据；
S2,存储时域数据；
S3,对时域数据进行信号处理，生成大气湍流回波功率谱数据；
S4,存储大气湍流回波功率谱数据；
S5,对大气湍流回波功率谱数据进行谱矩处理，得到多个大气风场的分量数据，将多个分量数据按获取的时间顺序进行有重叠分组，按分组的顺序处理分量数据，得到按时间顺序排列的大气湍流能量和摩擦速度并储存。
9.根据权利要求8所述的风廓线雷达大气湍流能量探测方法,其特征在于,所述S5中将多个分量数据按获取的时间顺序进行有重叠分组，具体为，假设组内分量数据的个数为N，则分量数据重叠的个数为N-1，其中重叠部分为连续获取的分量数据，每得到一个分量数据，与所述分量数据之前的N-1个分量数据分为一组。
10.根据权利要求9所述的风廓线雷达大气湍流能量探测方法,其特征在于,所述组内分量数据的个数为5，得到大气风场的第一组分量数据，处理得到一个大气湍流能量和摩擦速度的数值，按获得分量数据的时间顺序再获得一个分量数据，所述分量数据与所述分量数据之前的4个分量数据分为另一组，再处理得到一个大气湍流能量和摩擦速度的数值，依照分组的顺序处理分量数据得到按时间顺序排列的大气湍流能量和摩擦速度并储存。
CN201610851730.1A 2016-09-26 2016-09-26 一种风廓线雷达大气湍流能量探测系统及方法 Active CN107870332B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610851730.1A CN107870332B (zh) 2016-09-26 2016-09-26 一种风廓线雷达大气湍流能量探测系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610851730.1A CN107870332B (zh) 2016-09-26 2016-09-26 一种风廓线雷达大气湍流能量探测系统及方法

Publications (2)

Publication Number Publication Date
CN107870332A true CN107870332A (zh) 2018-04-03
CN107870332B CN107870332B (zh) 2019-11-22

Family

ID=61751950

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610851730.1A Active CN107870332B (zh) 2016-09-26 2016-09-26 一种风廓线雷达大气湍流能量探测系统及方法

Country Status (1)

Country Link
CN (1) CN107870332B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109188566A (zh) * 2018-08-15 2019-01-11 河海大学 利用沃克环流指数判断沃克环流强弱的方法
CN109254273A (zh) * 2018-11-01 2019-01-22 中国气象科学研究院 风廓线雷达回波信号的处理方法和装置
CN109581384A (zh) * 2019-01-28 2019-04-05 中国气象局气象探测中心 基于多普勒天气雷达的晴空垂直风廓线探测方法及系统
CN109828276A (zh) * 2019-03-06 2019-05-31 南京微麦科斯电子科技有限责任公司 一种基于干涉算法的风廓线雷达
CN112611921A (zh) * 2020-12-09 2021-04-06 上海无线电设备研究所 一种大气声场模拟装置及其电磁散射特性测试方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030151591A1 (en) * 2002-02-11 2003-08-14 Harpen James M. Weather data delivery system
US20080001808A1 (en) * 2004-12-30 2008-01-03 Passarelli Richard E Jr System and method for processing data in weather radar
US20090243918A1 (en) * 2005-08-17 2009-10-01 Kelly Terence F System and Method for Presenting Wind Speed Information
CN102323574A (zh) * 2011-06-25 2012-01-18 中国航天科工集团第二研究院二十三所 一种风廓线雷达信号处理方法
CN102508219A (zh) * 2011-10-17 2012-06-20 中国人民解放军理工大学气象学院 风廓线雷达湍流目标检测方法
CN202735511U (zh) * 2012-07-12 2013-02-13 南京信息工程大学 一种风廓线雷达信号处理系统
CN104133216A (zh) * 2014-07-17 2014-11-05 北京无线电测量研究所 一种获取低空风廓线的雷达探测方法及装置
CN104237892A (zh) * 2014-09-17 2014-12-24 北京无线电测量研究所 一种风廓线雷达大气温度探测系统及方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030151591A1 (en) * 2002-02-11 2003-08-14 Harpen James M. Weather data delivery system
US20080001808A1 (en) * 2004-12-30 2008-01-03 Passarelli Richard E Jr System and method for processing data in weather radar
US20090243918A1 (en) * 2005-08-17 2009-10-01 Kelly Terence F System and Method for Presenting Wind Speed Information
CN102323574A (zh) * 2011-06-25 2012-01-18 中国航天科工集团第二研究院二十三所 一种风廓线雷达信号处理方法
CN102323574B (zh) * 2011-06-25 2013-03-06 中国航天科工集团第二研究院二十三所 一种风廓线雷达信号处理方法
CN102508219A (zh) * 2011-10-17 2012-06-20 中国人民解放军理工大学气象学院 风廓线雷达湍流目标检测方法
CN202735511U (zh) * 2012-07-12 2013-02-13 南京信息工程大学 一种风廓线雷达信号处理系统
CN104133216A (zh) * 2014-07-17 2014-11-05 北京无线电测量研究所 一种获取低空风廓线的雷达探测方法及装置
CN104237892A (zh) * 2014-09-17 2014-12-24 北京无线电测量研究所 一种风廓线雷达大气温度探测系统及方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109188566A (zh) * 2018-08-15 2019-01-11 河海大学 利用沃克环流指数判断沃克环流强弱的方法
CN109188566B (zh) * 2018-08-15 2020-07-17 河海大学 利用沃克环流指数判断沃克环流强弱的方法
CN109254273A (zh) * 2018-11-01 2019-01-22 中国气象科学研究院 风廓线雷达回波信号的处理方法和装置
CN109581384A (zh) * 2019-01-28 2019-04-05 中国气象局气象探测中心 基于多普勒天气雷达的晴空垂直风廓线探测方法及系统
CN109581384B (zh) * 2019-01-28 2021-09-17 中国气象局气象探测中心 基于多普勒天气雷达的晴空垂直风廓线探测方法及系统
CN109828276A (zh) * 2019-03-06 2019-05-31 南京微麦科斯电子科技有限责任公司 一种基于干涉算法的风廓线雷达
CN112611921A (zh) * 2020-12-09 2021-04-06 上海无线电设备研究所 一种大气声场模拟装置及其电磁散射特性测试方法

Also Published As

Publication number Publication date
CN107870332B (zh) 2019-11-22

Similar Documents

Publication Publication Date Title
CN107870332B (zh) 一种风廓线雷达大气湍流能量探测系统及方法
Sathe et al. A six-beam method to measure turbulence statistics using ground-based wind lidars
Wiggs et al. The role of streamline curvature in sand dune dynamics: evidence from field and wind tunnel measurements
Kiemle et al. Estimation of boundary layer humidity fluxes and statistics from airborne differential absorption lidar (DIAL)
Bye et al. Estimating forest canopy parameters from satellite waveform LiDAR by inversion of the FLIGHT three-dimensional radiative transfer model
Fukao et al. Extensive studies of large‐amplitude Kelvin–Helmholtz billows in the lower atmosphere with VHF middle and upper atmosphere radar
CN107102309B (zh) 风场空间分布转化为风场时间方向的方法及气溶胶时空分布与风场时空分布的叠加方法
Sathe et al. Estimating turbulence statistics and parameters from ground-and nacelle-based lidar measurements: IEA Wind expert report
CN108959705B (zh) 一种预测亚热带森林生物量的方法
CN102998369B (zh) 一种二维损伤定量化监测方法
CN102262220A (zh) 一种基于非线性拟合的微震源或声发射源的定位方法
CN102902977A (zh) 基于风场扰动特征的空中目标分类方法
Doyle et al. Three-dimensional characteristics of stratospheric mountain waves during T-REX
Takahashi et al. Land–ocean differences in the warm‐rain formation process in satellite and ground‐based observations and model simulations
CN104792672A (zh) Pm2.5质量浓度的高度分布探测装置及其探测方法
CN107390229A (zh) 一种测风激光雷达信号的处理方法及其处理装置
Mayor et al. Two-component horizontal aerosol motion vectors in the atmospheric surface layer from a cross-correlation algorithm applied to scanning elastic backscatter lidar data
Tonttila et al. Cloud base vertical velocity statistics: a comparison between an atmospheric mesoscale model and remote sensing observations
CN104819963A (zh) 一种大气垂直能见度的测定方法及监控系统
Dallman et al. Flow and turbulence in an industrial/suburban roughness canopy
Anguelova et al. Using energy dissipation rate to obtain active whitecap fraction
Dodson et al. Droplet inhomogeneity in shallow cumuli: The effects of in-cloud location and aerosol number concentration
Egerer et al. Estimating turbulent energy flux vertical profiles from uncrewed aircraft system measurements: exemplary results for the MOSAiC campaign
CN103759915A (zh) 局部水收集系数的试验测量方法
CN105891078B (zh) 基于风廓线雷达的沙尘暴沙尘质量浓度定量反演估算方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant