CN104199002B - 抑制风力发电旋叶虚警目标的雷达信息处理方法 - Google Patents

抑制风力发电旋叶虚警目标的雷达信息处理方法 Download PDF

Info

Publication number
CN104199002B
CN104199002B CN201410466974.9A CN201410466974A CN104199002B CN 104199002 B CN104199002 B CN 104199002B CN 201410466974 A CN201410466974 A CN 201410466974A CN 104199002 B CN104199002 B CN 104199002B
Authority
CN
China
Prior art keywords
vane
wind
power
false
main clutter
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
Application number
CN201410466974.9A
Other languages
English (en)
Other versions
CN104199002A (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.)
Leihua Electronic Technology Research Institute Aviation Industry Corp of China
Original Assignee
Leihua Electronic Technology Research Institute Aviation Industry Corp of China
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 Leihua Electronic Technology Research Institute Aviation Industry Corp of China filed Critical Leihua Electronic Technology Research Institute Aviation Industry Corp of China
Priority to CN201410466974.9A priority Critical patent/CN104199002B/zh
Publication of CN104199002A publication Critical patent/CN104199002A/zh
Application granted granted Critical
Publication of CN104199002B publication Critical patent/CN104199002B/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
    • 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
    • G01S7/414—Discriminating targets with respect to background clutter

Landscapes

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

Abstract

本发明属于雷达信息处理技术领域，尤其涉及雷达空空下视探测过程中对地面风力发电旋叶虚警目标的雷达信息处理方法。包括以下步骤：(1)获取雷达的脉冲发射频率，积累脉冲数，共轭乘积矢量计算主杂波的多普勒频率；(2)根据主杂波的多普勒频率确定雷达多普勒谱上主杂波位置；(3)获取风力发电设备的旋叶旋转角速度，最大旋叶长度，计算风力发电设备的旋叶末端所对应的频谱滤波器宽度；(4)获取检测点所在的距离门号，检测点的功率，频谱平均噪声功率和主杂波展宽半径，以及根据判断的门限系数，得到单个距离门上以主杂波为中心，旋叶频谱展宽为半径范围内的信号点功率的二值向量分布；(5)判别与剔除虚警。

Description

抑制风力发电旋叶虚警目标的雷达信息处理方法
技术领域
本发明属于雷达信息处理技术领域，尤其涉及雷达空空下视探测过程中对地面风力发电旋叶虚警目标的雷达信息处理方法。
背景技术
现代机载脉冲多普勒(PD)雷达对空下视探测目标时，主要使用中重频(MPRF)波形工作模式。该模式采用全相参PD处理，即可实现对迎头/尾随目标的检测、跟踪，还可完成上视/下视全空域的目标监视，因而被广泛使用。
该模式目标检测在二维“距离-频率”平面上进行，主要利用目标与主杂波之间产生的多普勒频差实现频率分离，从而可使用恒虚警技术在清晰区完成目标检测，以最大化提高检测概率。
传统恒虚警检测方法使用二维滑窗计算检测点周围背景功率，如通过单元平均恒虚警(CFAR)检测处理方法，比较检测点与背景功率的大小关系来判定目标。在无地面干扰目标的情况下，对下视/上视探测状态，该方法可以获得较好的虚警率控制。
当地面存在大量的风力发电设备(尤其在北方、沿海地区)，风力发电设备叶片在风力的作用下产生旋转，从而产生与主杂波存在一定频率差的多普勒频率，该多普勒频率与叶片的长短、旋转速度、照射角度等有关。在角速度一定的情况下，旋转叶片在不同半径位置会产生不同的线速度，在雷达目标回波频谱上形成连续条带式、起伏调制的回波信号。一个地面风力发电场，一般都拥有几十、上百个独立的风力发电塔，而每一个独立风力发电塔的叶片，在风力作用下旋转，都有可能产生有别于主杂波的运动目标回波；下视状态，雷达天线波束照射到地面风力发电场，如果使用传统的恒虚警技术，风力发电旋转叶片回波将被检测为大量的空中虚警目标，从而影响雷达的检测性能与正常使用。
发明内容
本发明要解决的技术问题是提供一种抑制风力发电旋叶虚警目标的雷达信息处理方法，降低雷达下视探测引起的目标虚警概率，并保证不损失真实空中目标的检测概率。
为了解决上述技术问题，本发明包括以下步骤：
(1)获取雷达的脉冲发射频率，积累脉冲数和共轭乘积矢量计算主杂波的多普勒频率；
(2)根据主杂波的多普勒频率确定雷达多普勒谱上主杂波位置；
(3)获取风力发电设备的旋叶旋转角速度和最大旋叶长度，计算风力发电设备的旋叶末端所对应的频谱滤波器宽度；
(4)获取检测点所在的距离门号，检测点的功率，频谱平均噪声功率和主杂波展宽半径，以及根据判断的门限系数，得到单个距离门上以主杂波为中心，旋叶频谱展宽为半径范围内的信号点功率的二值向量分布；
(5)判别与剔除虚警。
本发明对真实空中目标的检测概率基本无影响。在雷达空空模式的下视探测过程中，能充分利用地面风力发电设备旋转叶片回波的频谱特征进行精确的检测定位，显著降低了雷达对空目标检测的虚警率。
附图说明
图1为风力发电旋叶频谱；
图2为抑制之前的检测图；
图3为抑制后的检测结果；
图4为采用本发明方法前的对空下视显示效果图；
图5为采用本发明方法后的对空下视显示效果图；
具体实施方式
本发明所列举的实施例，只是用于帮助理解本发明，不应理解为对本发明保护范围的限定，对于本技术领域的普通技术人员来说，在不脱离本发明思想的前提下，还可以对本发明进行改进和修饰，这些改进和修饰也落入本发明权利要求保护的范围内。
本发明包括以下步骤：
(1)获取雷达的脉冲发射频率，积累脉冲数，共轭乘积矢量，根据下列公式计算主杂波的多普勒频率：
f c = P R F 2 · π ∠ [ Σ n = 0 N a - 1 S * ( n ) S ( n + 1 ) ]
其中，fc为主杂波的多普勒频率S*(n)S(n+1)为雷达相邻方位样本数据的共轭乘积矢量，∠含义为相位运算，PRF为脉冲发射频率；Na为积累脉冲数；
(2)根据主杂波的多普勒频率确定雷达多普勒谱上主杂波位置，其计算公式如下：
N f c = mod [ f l o o r ( f c · N a P R F + N a ) , N a ]
其中，Nfc为主杂波在多普勒频谱上的位置，floor为运算符，含义为向下取整运算；
(3)获取风力发电设备的旋叶旋转角速度，最大旋叶长度，根据下列公式计算风力发电设备的旋叶末端所对应的频谱滤波器宽度：
其中，rc为风力发电设备的旋叶末端所对应的频谱滤波器宽度，为风力发电设备的旋叶旋转角速度，Lmax为最大旋叶长度，λ为雷达波长，在实际应用过程中rc和两个参数取能满足系统指标要求的统计值即可；
(4)获取检测点所在的距离门号，检测点的功率，频谱平均噪声功率和主杂波展宽半径，以及根据判断的门限系数，根据下列公式计算单个距离门上以主杂波为中心，旋叶频谱展宽为半径范围内的信号点功率的二值向量分布；
m i = 1 p N r Nf i > ( M e a n P o w V a l + K ) 0 p N r Nf i ≤ ( M e a n P o w V a l + K )
其中，Nr为检测点所在的距离门号，p(Nr,Nf)为检测点的功率，MeanPowVal为频谱平均噪声功率，K为判断的门限系数，i的取值范围为其中Δr为主杂波展宽半径。
(5)判别与剔除虚警，依据下式算判定在主杂波的距离门上是否存在地面风力发电目标，并予以剔除，
T N r , N f = 1 Σ i m i ≥ p e r * ( r c - Δ r ) 0 e l s e
其中，per为系统设定的风力发电目标判定百分比系数。
下面对本发明进行试验，验证本发明的技术效果：
参与试验的雷达系统参数为PRF＝10KHz，x波段Na＝256，vmax＝20m·s-1，Δr＝3，MeanPowVal＝70dB，K＝10dB。测得风力发电旋叶频谱如图1所示，它的横坐标是距离门的样本点序号，纵坐标是频率门的序列，图中位于70号频率门贯穿全部距离门的为主杂波功率信息。可以看出以主杂波为中心，在多个距离门上延伸出了多条频率方向的条带(图1中的白色小框指示的位置)，即地面风力发电设备的旋叶自转引起的频率展宽。通过传统方法检测可得到如图2所示的检测图。
根据系统参数和原始回波数据按照公式[1]、[2]计算得到Nfc＝62，依据公式[3]计算得到最大的旋叶频谱展宽半径rc＝16。对主杂波上的距离门数据按照公式[4]、[5]实现虚警抑制，可得到如图3所示的检测图，通过对比图2和图3检测结果之间的差别，可以看出使用本发明方法，对风力发电设备旋叶产生的虚警目标具有较好的抑制效果。
图4和图5分别为采用本发明进行雷达信息处理前、后，雷达对空下视探测的显示效果图，可见，使用传统检测方法对空下视探测时，天线波束扫描到某一个风力发电场后，会检测出非常多的虚假空中目标航迹(如图4)，采用本发明方法后，风力发电旋叶的回波虚警得到了很好的抑制，只显示出合作目标航迹(如图5)。
空中试飞结果证明，采用本发明的技术方法，在有地面风力发电旋转叶片回波干扰的环境下，平均目标航迹虚警间隔时间大于1h，远大于1min的虚警间隔时间要求。

Claims (6)

1.抑制风力发电旋叶虚警目标的雷达信息处理方法，其特征在于，包括以下步骤：
(1)获取雷达的脉冲发射频率，积累脉冲数和共轭乘积矢量计算主杂波的多普勒频率；
(2)根据主杂波的多普勒频率确定雷达多普勒谱上主杂波位置；
(3)获取风力发电设备的旋叶旋转角速度和最大旋叶长度，计算风力发电设备的旋叶末端所对应的频谱滤波器宽度；
(4)获取检测点所在的距离门号，检测点的功率，频谱平均噪声功率和主杂波展宽半径，以及根据判断的门限系数，得到单个距离门上以主杂波为中心，旋叶频谱展宽为半径范围内的信号点功率的二值向量分布；
(5)判别与剔除虚警。
2.根据权利要求1所述的抑制风力发电旋叶虚警目标的雷达信息处理方法，其特征在于，雷达的多普勒频率由下式算出：
其中，fc为主杂波的多普勒频率S*(n)S(n+1)为雷达相邻方位样本数据的共轭乘积矢量，∠含义为相位运算，PRF为脉冲发射频率；Na为积累脉冲数。
3.根据权利要求1所述的抑制风力发电旋叶虚警目标的雷达信息处理方法，其特征在于，主杂波位置由下式算出：
其中，Nfc为主杂波在多普勒频谱上的位置，floor为运算符，含义为向下取整运算。
4.根据权利要求1所述的抑制风力发电旋叶虚警目标的雷达信息处理方法，其特征在于，风力发电设备的旋叶末端所对应的频谱滤波器宽度由下式算出：
其中，rc为风力发电设备的旋叶末端所对应的频谱滤波器宽度，为风力发电设备的旋叶旋转角速度，Lmax为最大旋叶长度，λ为雷达波长，在实际应用过程中rc和两个参数取能满足系统指标要求的统计值即可。
5.根据权利要求1所述的抑制风力发电旋叶虚警目标的雷达信息处理方法，其特征在于，由下式得到单个距离门上以主杂波为中心，旋叶频谱展宽为半径范围内的信号点功率的二值向量分布：
其中，Nr为检测点所在的距离门号，p(Nr,Nf)为检测点的功率，MeanPowVal为频谱平均噪声功率，K为判断的门限系数，i的取值范围为其中Δr为主杂波展宽半径。
6.根据权利要求1所述的抑制风力发电旋叶虚警目标的雷达信息处理方法，其特征在于，依据下式算判定在主杂波的距离门上是否存在地面风力发电目标，并予以剔除，
其中，per为系统设定的风力发电目标判定百分比系数。
CN201410466974.9A 2014-09-12 2014-09-12 抑制风力发电旋叶虚警目标的雷达信息处理方法 Active CN104199002B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410466974.9A CN104199002B (zh) 2014-09-12 2014-09-12 抑制风力发电旋叶虚警目标的雷达信息处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410466974.9A CN104199002B (zh) 2014-09-12 2014-09-12 抑制风力发电旋叶虚警目标的雷达信息处理方法

Publications (2)

Publication Number Publication Date
CN104199002A CN104199002A (zh) 2014-12-10
CN104199002B true CN104199002B (zh) 2016-09-14

Family

ID=52084314

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410466974.9A Active CN104199002B (zh) 2014-09-12 2014-09-12 抑制风力发电旋叶虚警目标的雷达信息处理方法

Country Status (1)

Country Link
CN (1) CN104199002B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109557514B (zh) * 2019-01-14 2023-05-02 三峡大学 一种风电机叶片回波准确求解方法
CN111060880B (zh) * 2019-12-30 2021-09-07 无锡市雷华科技有限公司 一种基于恒虚警检测原理的气象杂波抑制方法
CN111123236B (zh) * 2019-12-30 2021-10-29 无锡市雷华科技有限公司 地面固定动目标自适应抑制方法
CN111198359B (zh) * 2020-01-11 2022-07-05 桂林电子科技大学 一种基于目标多普勒宽度的旋翼无人机检测方法
CN115234447B (zh) * 2022-07-08 2023-07-25 南京牧镭激光科技股份有限公司 一种基于海上净空雷达的叶片检测方法及系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5111778B2 (ja) * 2006-04-28 2013-01-09 株式会社東芝 レーダ装置とその信号処理方法
CN101251596B (zh) * 2008-03-21 2010-09-22 北京理工大学 一种风廓线雷达非相参无线电干扰抑制方法
FR2954516B1 (fr) * 2009-12-18 2012-09-21 Thales Sa Procede de detection radar, notamment pour des radars aeroportes mettant en oeuvre une fonction de detection et d'evitement d'obstacles
CN103197301B (zh) * 2013-03-19 2014-11-19 中国人民解放军海军航空工程学院 海面微动目标Radon-线性正则变换长时间相参积累检测方法
CN103344949B (zh) * 2013-06-18 2015-03-18 中国人民解放军海军航空工程学院 基于Radon-线性正则模糊函数的雷达微弱动目标检测方法

Also Published As

Publication number Publication date
CN104199002A (zh) 2014-12-10

Similar Documents

Publication Publication Date Title
CN104199002B (zh) 抑制风力发电旋叶虚警目标的雷达信息处理方法
CN103454624B (zh) 基于降维稀疏重构空时谱的直接数据域动目标检测方法
CN103837863B (zh) 基于梯度投影的距离‑速度同步拖引欺骗干扰识别
CN104215951B (zh) 一种在海杂波背景下慢速小目标的检测系统及其检测方法
CN102508219B (zh) 风廓线雷达湍流目标检测方法
CN103439697B (zh) 基于动态规划的目标检测方法
CN102879767B (zh) 一种用于空时自适应处理的干扰目标检测方法
CN105093215A (zh) 基于多普勒信息的雷达对低空慢速小目标的跟踪方法
CN104931949B (zh) 雷达扫描模式下风轮机杂波背景下飞机目标检测方法
US9268008B1 (en) Detection of low observable objects in clutter using non-coherent radars
CN105487060B (zh) 一种双通道四斜率调制的多目标提取方法
CN104714225B (zh) 一种基于广义似然比的动态规划检测前跟踪方法
CN104166135A (zh) 一种宽带雷达目标的原始点迹凝聚处理方法
CN103558597B (zh) 基于谱峭度的海杂波中弱目标检测方法
CN103176187B (zh) 一种机载预警雷达地面高速公路目标过滤方法
CN102721954A (zh) 基于动态杂波图的风电场杂波识别与抑制方法
CN102298142A (zh) 雷达测距模糊条件下多微弱目标ht-stc-tbd检测方法
CN106199548B (zh) 基于四极化通道融合的海面漂浮微弱雷达目标的检测方法
CN105116387A (zh) 基于位置和多普勒速度信息的pd雷达抗速度拖引方法
Zheng et al. Gust front statistical characteristics and automatic identification algorithm for CINRAD
CN112198487A (zh) 一种风电场杂波背景下的目标检测方法
CN102621536A (zh) 基于relax的空中多机动目标检测与参数估计方法
CN106468772A (zh) 一种基于距离‑多普勒量测的多站雷达人体跟踪方法
CN202433521U (zh) 风廓线雷达湍流目标检测处理板
CN106443660A (zh) 一种线性调频连续波雷达多运动目标匹配方法

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