CN111624623B - 基于激光雷达非均匀扫描的风场反演方法 - Google Patents
基于激光雷达非均匀扫描的风场反演方法 Download PDFInfo
- Publication number
- CN111624623B CN111624623B CN202010487703.7A CN202010487703A CN111624623B CN 111624623 B CN111624623 B CN 111624623B CN 202010487703 A CN202010487703 A CN 202010487703A CN 111624623 B CN111624623 B CN 111624623B
- Authority
- CN
- China
- Prior art keywords
- wind
- scanning
- measurement
- vector
- wind direction
- 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.)
- Active
Links
Images
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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
- G01S17/95—Lidar systems specially adapted for specific applications for meteorological use
-
- 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
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
基于激光雷达非均匀扫描的风场反演方法,包括1)激光雷达以非均匀的扫描方式对矢量风场进行测量;2)对扫描数据进行拟合;3)将拟合曲线的差值控制在预设范围;4)将剩余波束控制得数量大于有效波束数量的阈值,并反演整体测量区域的风速风向;5)取前一步所有探测距离的有效风向进行平均;6)将风向均值反馈至步骤1)并重复1)—4);7)当再一次反演出整体测量区域的风速风向之后,利用VPP方法和加权统计得到精确的目标区域矢量;8)若收到指令则结束测量,将上一步测量结果输出;否则进行风向均值计算重复步骤1)‑4)、7)。本发明在提高单台激光雷达矢量风场反演的数据更新率的同时,保证了目标测量点的测量精度。
Description
技术领域
本发明涉及一种改进的激光雷达风场反演方法,具体涉及一种激光雷达基于非均匀扫描的风场反演方法。
背景技术
大气风场是气象探测中的重要观测参数,天气预报中提供风场的基础信息,有助于天气变化的预报;航空中的阵风、风切变等情况危害很大,风场的观测可以帮助提前做好应对准备;风力发电中也需要前方风场的测定来使风机有效运行等。扫描式激光雷达可实时提供三维空间内的风场信息,时空分辨率高,体积小,方便携带,在不方便进行塔杆搭建的地方,可实现远距离探测,广泛用于气象、环境、交通、航空、海洋、风电等诸多领域。
激光雷达将激光发射到空中,与大气分子和气溶胶颗粒等大气粒子相互作用后产生回波信号,回波信号中的包含的气溶胶的径向速度信息,通过接收器接收后向推算出风场的径向速度。
扫描式激光雷达可以进行三维风场扫描,作为一种业务化的产品,需要实时高效地进行风场测量。
多部雷达联合观测,水平风场可由其几何关系导出。
两部甚至多部雷达风场反演的空间范围有限,且只有雷达同时同步扫描的空间的数据才可用于反演。
另外受到地形、经济条件等方面的限制,实际上只有单台雷达可进行风场测量。
单激光雷达反演风场的方法是以径向速度作为风场分量,反演整体风场后再对目标区域数据进行与径向分量非共线分量的计算,多个非共线分量可精确反演出包括风速、风向的风场信息。进行远距离整体风场扫描测量时常采用均匀间隔的扫描方式,连续紧密波束扫描扇形区,波束多,耗时长,造成时间和资源浪费,数据更新率不能保证;若进行大间隔扫描,减少波束会降低测量区域的风场反演精度。
发明内容
针对上述技术问题,本发明的目的是提供一种基于激光雷达非均匀扫描的风场反演方法,以弥补单激光雷达风场测量技术精度低效率低的缺点,实现高精度,高效率的矢量风场测量。
基于激光雷达非均匀扫描的风场反演方法,其特征是包括以下步骤,
(1)激光雷达以非均匀的扫描方式对矢量风场进行测量:
选一个角度θ作为中心方位角,开始非均匀扫描,与中心角度越近的波束扫描角度密集,越远的波束扫描角度稀疏,以加大中心方位角波束权重,波束扫描角度具体满足一种非线性曲线函数特征,
k为起始扫描波束到结束扫描波束数的编号,k=1,2,3......N......2N+1,N为自然数;扫描间隔f(m)为相邻两个波束扫描角度之差,满足公式(1),其中,m为从起始波束到结束波束的扫描间隔数,a为1~3范围内任意实数,m=1,2,3......N......2N,N为自然数;
f(m)=exp(a(m/N-1)2) (1)
(2)对于第k个波束,k=1,2,3......N......2N+1,根据非线性风场VAD(VelocityAzimuth Display速度方位显示)反演方法,测量其风分量数据为vk,且vk与构成一种正弦曲线函数可根据正弦函数对扫描数据进行拟合;
拟合过程是利用最小二乘法计算拟合曲线和测量数据之间的差值,如公式(3),假设每次参与拟合的各个波束数编号按升序表示为kb~ke,则差值计算方法为:
而首次拟合计算中,2N+1个波束均参与计算,即此时kb~ke表示1~2N+1之间的所有波束;
(3)判断拟合曲线和测量数据之间的差值e是否达到预设的误差阈值Te(e<Te),若大于误差阈值,说明差值过大,会降低风速反演精度,需要对差值过大的数据进行剔除后,剩余波束重复步骤(2),直到将差值控制在预设的范围内;
(4)将剩余的波束数量跟预设的有效波束数量阈值T进行比较,若数量大于T,则进行下一步;若数量小于T,则仍取步骤(1)的中心方位角θ,重复步骤(1)-(3),直至数量大于T,并反演整体测量区域的风速风向,也即反演所有探测距离的风速风向;
(5)由于与风向方向具有相同方位角的波束测量到的径向速度等于风速在径向仰角上的投影,将风向设置为中心方位角,有利于提高风速的反演精度,在地转偏向力的影响下,每个波束径向上不同探测距离的风向不同,
取步骤(4)的所有探测距离的有效风向进行平均(其目的在在于将所得的风向均值作为第二次测量的中心方位角,以自动调整激光雷达伺服扫描指向角度);
(6)返回步骤(1),并将步骤(5)所得的风向均值作为步骤(1)中θ的取值,利用步骤(1)的方式进行计算,然后进行步骤(2)、(3)、(4)的操作,以重新反演整体测量区域的风速风向;
(7)当再一次反演出整体测量区域的风速风向之后,投影到具体的目标区域,利用VPP方法,结合整场风矢量方向的径向分量和目标区域实测径向分量,求解得到与实测径向分量非共线的非共线分量结合实测分量和非共线分量反演目标区域内不同测量点的矢量风场其中:
(8)判断是否收到测量结束指令,若收到指令则结束测量,将上一步所得的精确目标区域矢量作为测量结果;
若没有收到则测量结束指令,则以与步骤(5)同样的方式对本次反演的所有探测距离的有效风向进行平均值计算,所得的风向均值作为下一次扫描的中心方位角度θ,继续进行扫描,重复步骤(1)、(2)、(3)(4)、(7),直至收到测量结束指令。
本发明的非均匀扫描的单激光雷达矢量风场反演方法在提高单台激光雷达矢量风场反演的数据更新率的同时,保证了目标测量点的测量精度。
附图说明
图1是本发明采用非均匀扫描的示意图。
图2是本发明实施的技术流程图。
其中,1.激光雷达,2.扫描波束,3.扫描角度,4.测量数据,5.目标区域。
具体实施方式
基于激光雷达非均匀扫描的风场反演方法,如图2所示,包括以下步骤,
(1)如图1所示,激光雷达以非均匀的扫描方式对矢量风场进行测量:
选一个角度θ作为中心方位角,开始非均匀扫描,与中心角度越近的波束扫描角度密集,越远的波束扫描角度稀疏,以加大中心方位角波束权重,波束扫描角度具体满足一种非线性曲线函数特征,
k为起始扫描波束到结束扫描波束数的编号,k=1,2,3......N......2N+1,N为自然数;扫描间隔f(m)为相邻两个波束扫描角度之差,满足公式(1),其中,m为从起始波束到结束波束的扫描间隔数,a为1~3范围内任意实数,m=1,2,3......N......2N,N为自然数;
f(m)=exp(a(m/N-1)2) (1)
(2)对于第k个波束,k=1,2,3......N......2N+1,根据非线性风场VAD(VelocityAzimuth Display速度方位显示)反演方法,测量其风分量数据为vk,且vk与构成一种正弦曲线函数可根据正弦函数对扫描数据进行拟合;
拟合过程是利用最小二乘法计算拟合曲线和测量数据之间的差值,如公式(3),假设每次参与拟合的各个波束数编号按升序表示为kb~ke,则差值计算方法为:
而首次拟合计算中,2N+1个波束均参与计算,即此时kb~ke表示1~2N+1之间的所有波束;
(3)判断拟合曲线和测量数据之间的差值e是否达到预设的误差阈值Te(e<Te),若大于误差阈值,说明差值过大,会降低风速反演精度,需要对差值过大的数据进行剔除后,剩余波束重复步骤(2),直到将差值控制在预设的范围内;
(4)将剩余的波束数量跟预设的有效波束数量阈值T进行比较,若数量大于T,则进行下一步;若数量小于T,则仍取步骤(1)的中心方位角θ,重复步骤(1)-(3),直至数量大于T,并反演整体测量区域的风速风向,也即反演所有探测距离的风速风向;
(5)由于与风向方向具有相同方位角的波束测量到的径向速度等于风速在径向仰角上的投影,将风向设置为中心方位角,有利于提高风速的反演精度,在地转偏向力的影响下,每个波束径向上不同探测距离的风向不同,
取步骤(4)的所有探测距离的有效风向进行平均(其目的在在于将所得的风向均值作为第二次测量的中心方位角,以自动调整激光雷达伺服扫描指向角度);
(6)返回步骤(1),并将步骤(5)所得的风向均值作为步骤(1)中θ的取值,利用步骤(1)的方式进行计算,然后进行步骤(2)、(3)、(4)的操作,以重新反演整体测量区域的风速风向;
(7)当再一次反演出整体测量区域的风速风向之后,投影到具体的目标区域,利用VPP方法,结合整场风矢量方向的径向分量和目标区域实测径向分量,求解得到与实测径向分量非共线的非共线分量结合实测分量和非共线分量反演目标区域内不同测量点的矢量风场其中:
(8)判断是否收到测量结束指令,若收到指令则结束测量,将上一步所得的精确目标区域矢量作为测量结果;
若没有收到则测量结束指令,则以与步骤(5)同样的方式对本次反演的所有探测距离的有效风向进行平均值计算,所得的风向均值作为下一次扫描的中心方位角度θ,继续进行扫描,重复步骤(1)、(2)、(3)(4)、(7),直至收到测量结束指令。
随着风向在变化,除了第一次是预设的初始方位角,以后每次扫描都是以上一轮测量中步骤(5)得到的风向作为中心角度,进行步骤(1)、(2)、(3)、(4)、(7)的操作,并调整伺服指向;直至达到设定的测量次数或收到测量结束信号。
其中加权统计如:按照目标区域内的测量点相对于目标区域中心点远近进行权重分配,测量点数量等于目标区域内覆盖的径向波束数量和每个径向上的距离点数的乘积,设目标区域测量点总数为j,具体通过对目标区域内测量点按照距离远近进行赋值,测量点风场数据为x=x1,x2,x3……xj,对应的赋值为i=i1,i2,i3……ij,距离中心点越近i值越大,距离中心点越远i值越小,测量点风场数据x1对应的权重系数为w1=i1/(i1+i2+i3+......+ij),测量点风场数据x2对应的权重系数为w2=i2/(i1+i2+i3+......+ij),以此类推,可获取每个测量点的权重,最后将测量点数据按照权重进行加权平均,得到最终目标区域风场结果为:
上述步骤(1)中,所述的中心方位角,在计算出风向之前是设置的初始方位角,在步骤(3)计算出风向后,初始方位角即被计算的风向替代,之后每次扫描都是以上一次计算的风向为中心方位角,计算波束角度并调整伺服指向。
上述步骤(2)中,第一次拟合求差时,kb=1,ke=2N+1,剔除超出范围波束数据之后再拟合求差时,剩余波束数编号按照升序排列,kb为剩余波束数的最小编号,ke为剩余波束数的最大编号。
Claims (1)
1.基于激光雷达非均匀扫描的风场反演方法,其特征是包括以下步骤,
(1)激光雷达以非均匀的扫描方式对矢量风场进行测量:
选一个角度θ作为中心方位角,开始非均匀扫描,与中心角度越近的波束扫描角度密集,越远的波束扫描角度稀疏,以加大中心方位角波束权重,波束扫描角度具体满足一种非线性曲线函数特征,
k为起始扫描波束到结束扫描波束数的编号,k=1,2,3......N......2N+1,N为自然数;扫描间隔f(m)为相邻两个波束扫描角度之差,满足公式(1),其中,m为从起始波束到结束波束的扫描间隔数,a为1~3范围内任意实数,m=1,2,3......N......2N,N为自然数;
拟合过程是利用最小二乘法计算拟合曲线和测量数据之间的差值,如公式(3),假设每次参与拟合的各个波束数编号按升序表示为kb~ke,则差值计算方法为:
而首次拟合计算中,2N+1个波束均参与计算,即此时kb~ke表示1~2N+1之间的所有波束;
(3)判断拟合曲线和测量数据之间的差值e是否达到预设的误差阈值Te,若大于误差阈值,说明差值过大,会降低风速反演精度,需要对差值过大的数据进行剔除后,剩余波束重复步骤(2),直到将差值控制在预设的范围内;
(4)将剩余的波束数量跟预设的有效波束数量阈值T进行比较,若数量大于T,则进行下一步;若数量小于T,则仍取步骤(1)的中心方位角θ,重复步骤(1)-(3),直至数量大于T,并反演整体测量区域的风速风向,也即反演所有探测距离的风速风向;
(5)取步骤(4)的所有探测距离的有效风向进行平均;
(6)返回步骤(1),并将步骤(5)所得的风向均值作为步骤(1)中θ的取值,利用步骤(1)的方式进行计算,然后进行步骤(2)、(3)、(4)的操作,以重新反演整体测量区域的风速风向;
(7)当再一次反演出整体测量区域的风速风向之后,投影到具体的目标区域,利用VPP方法,结合整场风矢量方向的径向分量和目标区域实测径向分量,求解得到与实测径向分量非共线的非共线分量结合实测分量和非共线分量反演目标区域内不同测量点的矢量风场其中:
(8)判断是否收到测量结束指令,若收到指令则结束测量,将上一步所得的精确目标区域矢量作为测量结果;
若没有收到则测量结束指令,则以与步骤(5)同样的方式对本次反演的所有探测距离的有效风向进行平均值计算,所得的风向均值作为下一次扫描的中心方位角度θ,继续进行扫描,重复步骤(1)、(2)、(3)、(4)、(7),直至收到测量结束指令。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010487703.7A CN111624623B (zh) | 2020-06-02 | 2020-06-02 | 基于激光雷达非均匀扫描的风场反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010487703.7A CN111624623B (zh) | 2020-06-02 | 2020-06-02 | 基于激光雷达非均匀扫描的风场反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111624623A CN111624623A (zh) | 2020-09-04 |
CN111624623B true CN111624623B (zh) | 2023-03-10 |
Family
ID=72272855
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010487703.7A Active CN111624623B (zh) | 2020-06-02 | 2020-06-02 | 基于激光雷达非均匀扫描的风场反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111624623B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112965084B (zh) * | 2021-01-28 | 2021-10-22 | 中国人民解放军国防科技大学 | 基于激光雷达的机场风场特征探测方法、装置和设备 |
CN115358073A (zh) * | 2022-08-19 | 2022-11-18 | 成都信息工程大学 | 一种ads-b风矢量反演方法 |
CN117420569B (zh) * | 2023-12-19 | 2024-03-12 | 南京牧镭激光科技股份有限公司 | 一种基于多普勒激光测风雷达的非均匀风场的反演方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002065155A1 (en) * | 2001-02-09 | 2002-08-22 | Commonwealth Scientific And Industrial Research Organisation | Lidar system and method |
JP2009276150A (ja) * | 2008-05-13 | 2009-11-26 | Ihi Corp | レーザレーダ及びレーザレーダの据付方向調整方法 |
CN101853335A (zh) * | 2010-06-01 | 2010-10-06 | 国家卫星海洋应用中心 | 一种海面风场反演的点方式反演方法 |
WO2011138895A1 (ja) * | 2010-05-07 | 2011-11-10 | 三菱電機株式会社 | レーザレーダ装置 |
CN103994719A (zh) * | 2014-05-30 | 2014-08-20 | 中国科学院国家天文台南京天文光学技术研究所 | 基于盖革apd阵列的高精度三维成像装置及其使用方法 |
CN104297117A (zh) * | 2014-10-23 | 2015-01-21 | 浙江省环境保护科学设计研究院 | 基于遥感技术的风景名胜区道路交通污染预警装置及方法 |
CN105467393A (zh) * | 2016-01-06 | 2016-04-06 | 北京无线电测量研究所 | 一种毫米波云雷达三维风场反演方法及系统 |
CN110456382A (zh) * | 2019-07-12 | 2019-11-15 | 中国海洋大学 | 基于单多普勒激光雷达的非均匀风场风矢量的测量方法 |
-
2020
- 2020-06-02 CN CN202010487703.7A patent/CN111624623B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002065155A1 (en) * | 2001-02-09 | 2002-08-22 | Commonwealth Scientific And Industrial Research Organisation | Lidar system and method |
JP2009276150A (ja) * | 2008-05-13 | 2009-11-26 | Ihi Corp | レーザレーダ及びレーザレーダの据付方向調整方法 |
WO2011138895A1 (ja) * | 2010-05-07 | 2011-11-10 | 三菱電機株式会社 | レーザレーダ装置 |
CN101853335A (zh) * | 2010-06-01 | 2010-10-06 | 国家卫星海洋应用中心 | 一种海面风场反演的点方式反演方法 |
CN103994719A (zh) * | 2014-05-30 | 2014-08-20 | 中国科学院国家天文台南京天文光学技术研究所 | 基于盖革apd阵列的高精度三维成像装置及其使用方法 |
CN104297117A (zh) * | 2014-10-23 | 2015-01-21 | 浙江省环境保护科学设计研究院 | 基于遥感技术的风景名胜区道路交通污染预警装置及方法 |
CN105467393A (zh) * | 2016-01-06 | 2016-04-06 | 北京无线电测量研究所 | 一种毫米波云雷达三维风场反演方法及系统 |
CN110456382A (zh) * | 2019-07-12 | 2019-11-15 | 中国海洋大学 | 基于单多普勒激光雷达的非均匀风场风矢量的测量方法 |
Non-Patent Citations (2)
Title |
---|
《All-fibermultifunction continuous-wave coherent laser radar at 1.55um for range,speed,vibration,and wind measurements》;Chirster J K et.al;《Applied Optics》;20201231;第39卷(第21期);第16-18页 * |
1.55μm全光纤多普勒测风雷达系统与试验;冯力天等;《红外与激光工程》;20110525(第05期);第844-847页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111624623A (zh) | 2020-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111624623B (zh) | 基于激光雷达非均匀扫描的风场反演方法 | |
Smalikho et al. | Lidar investigation of atmosphere effect on a wind turbine wake | |
CN106990401B (zh) | 基于全波形机载激光雷达数据二类高程误差修正方法 | |
CN110456382B (zh) | 基于单多普勒激光雷达的非均匀风场风矢量的测量方法 | |
CN104597430B (zh) | 机载气象雷达低空风切变三维回波仿真方法 | |
CN109001845B (zh) | 一种降雨量的预报方法 | |
CN102176003B (zh) | 一种机载激光雷达航测参数的优化设计方法 | |
CN109100723B (zh) | 基于多普勒天气雷达数据的高空风反演方法 | |
CN107765226A (zh) | 一种sar卫星雷达回波模拟方法、系统和介质 | |
CN109633601B (zh) | 基于地表模型的星载激光雷达脚点精确定位方法 | |
CN115236615B (zh) | 基于t矩阵法的机载极化气象雷达降水粒子回波仿真方法 | |
CN107843895A (zh) | 一种双多普勒雷达三维风场反演方法 | |
CN111638521A (zh) | 一种成像高度计遥感数据的海面风速反演方法 | |
CN109703770A (zh) | 基于测风激光雷达和cfd数据库的舰载机助降方法 | |
CN109085556B (zh) | 一种基于一二阶峰比值的高频地波雷达浪场形成方法 | |
CN111736150A (zh) | 一种面向远距离低功耗探鸟雷达的探测方法 | |
Antoniou et al. | Offshore wind profile measurements from remote sensing instruments | |
CN102426352A (zh) | 风廓线雷达风计算方法 | |
CN116520312B (zh) | 基于毫米波雷达、激光雷达数据融合的测风方法 | |
CN117452382A (zh) | 基于ppi误差寻优的三维激光测风雷达标定算法 | |
CN115508580B (zh) | 基于激光遥感技术的机场跑道虚拟风杆构建方法 | |
CN114740496B (zh) | 一种基于高阶泰勒展开的三维风场反演方法 | |
CN111190178A (zh) | 一种相频扫三坐标雷达的多波束顺序比幅测高方法 | |
CN111929702B (zh) | 空中目标大气扰动变分辨率探测方法、存储介质和系统 | |
CN107144826B (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 |