基于包络约束的幅相联合净空测量方法.pdf

80第40 卷第6 期2023年6 月真机计仿算文章编号：10 0 6-9 348（2 0 2 3)0 6-0 0 8 0-0 5基于包络约束的幅相联合净空测量方法邓佳欣，马丁峰？，徐云山？（1.中国飞行试验研究院，陕西西安7 10 0 8 9；2.中国人民解放军9 59 6 6 部队，黑龙江哈尔滨150 0 6 0）摘要：针对多径和塔筒震动引起的杂波区叶片检测困难问题，基于叶片的先验边界约束，应用距离多普勒二维幅相联合掩膜检测方法，能够剔除非叶片的杂波信息，同时利用叶片的近似线性变化特性，结合最小二乘拟合和单元平均进一步提升目标的相位估计精度，保证叶片的净空估计连续性和准确性。所提方法能够显著改善叶片在强噪声和杂波环境下的干涉相位估计精度，提升系统对于叶片净空的测量稳定性。利用风场采集的地面雷达实测数据对所提方法进行了验证分析，上述算法能够有效剔除虚假点，检测叶片目标，并在各种工况下准确提取叶片目标的干涉相位，同时所提方法计算复杂度较低，对于低成本化设计具有很高的工程价值。关键词：民用雷达应用；多径和虚警点抑制；干涉相位估计；叶片净空监测中图分类号：TP301.6文献标识码：BA Distance Measurement Method Based on EnvelopeConstraint with Amplitude and Phase InformationDENG Jia-xin,MA Ding-feng,XU Yun-shan?(1.Chinese Flight Test Establishment,Xian Shaanxi 710089,China;2.95966 Troops of the PLA,Harbin Heilongjiang 150060,China)ABSTRACT:Aiming at the difficulty of blade detection in clutter area caused by multipath and tower vibration,based on the prior boundary constraint of the blade,this method applies the range-doppler two-dimensional amplitudeand phase joint mask detection method to eliminate the false alarm points of non-blade targets.At the same time,u-sing the approximate linear variation characteristics of the blade,combined with the least square fitting and unit aver-age,the phase estimation accuracy of the target is further improved to ensure the continuity and accuracy of the bladeclearance estimation.The proposed method can significantly improve the accuracy of interference phase estimation ofblades in strong noise and clutter environment,and the stability of blade clearance measurement.The method is veri-fied and analyzed by using the measured data of ground radar collected by the wind field.The algorithm can effective-ly eliminate false points,detect blade targets,and accurately extract the interference phase of blade targets under vari-ous working conditions.Meanwhile,the proposed method has high engineering value for its low-cost and low computa-tional complexity design.KEYWORDS:Civil radar applications;False alarm point suppression;Interferometric phase estimation;Bladeclearance monitoring1引言由于国家十四五期间大力倡导低碳能源政策，因此大量用户会进一步加大新能源（风能和光伏发电）的投入力度，并结合微电网综合管理系统在用电的高峰和低谷期尽可能与收稿日期：2 0 2 1-0 9-0 3修回日期：2 0 2 1-0 9-2 3传统能源行业协同工作，为低碳环保做出贡献。其中风力发电机由于其峰值发电功率的不断提升，其风机的规模会越来越大，塔筒和叶片均达到了近百米量级，因此相应的也存在了运行过程当中的潜在风险。根据用户对大型风力发电机机组运行过程中的安全监控需求，需要高精度监测叶片在高速旋转过程中与风机塔筒之间的水平距离，评估叶片扫塔概率，在发生紧急状况之前给后台提供预警信息，以防安全事81故的发生。为了得到叶片在风机坐标系下的实时三维坐标，需要采用雷达监控设备获取目标的三坐标信息，从而精确测量目标的径向距离、方位角、俯仰角信息。文献1 介绍了三坐标雷达系统的发展历程，其相对于两坐标雷达，具有更强的目标观测能力。文献2,3 分别介绍了俄罗斯和荷兰比较成熟的三坐标雷达系统，并对其雷达性能进行了全面的分析。文献4-10分别从三坐标雷达系统设计、扫描体制、多维跟踪以及目标连续观测的高度闪烁等层面提供了不同的解决方案，为三坐标雷达的检测、跟踪及测高等提供了技术途径，为后续的雷达系统提供了大量借鉴思路。对于风力发电机的叶片净空测量系统，由于叶片下半段与塔筒之间的净空测量受方位维测量的影响较小，其测量误差对于7 0-10 0 米叶片的精度损失在0.2 米量级，因此可以借助分时和波束技术实现方位维测量的粗估计11,12 。然而俯仰方向的角度测量精度对叶片净空的估计精度影响非常大，对于7 0-10 0 米叶片的净空估计精度，1度对应的估计偏差达到了1米量级，因此必须对俯仰向进行高精度测角。文献13-16 从多个角度提供了干涉测角的测量方法。文献13从雷达干涉仪的测量误差、基线影响、测角模糊以及测量精度等多个角度对干涉仪雷达系统进行了比较全面的分析。文献14 提出了一种可以应用于可全方位同时探测的米波雷达天线及相应的测角算法，天线尺寸小、结构简单。文献15 提出了一种新的干涉测角数据处理算法,该算法对相位干涉仪测角中因相位解模糊出错产生的错误角度值进行野值剔除，并引人2 阶修正卡尔曼滤波器对剔除野值之后的角度信息进行平滑修正，保证了稳定跟踪的连续性。文献16提出了一种利用空间谱估计对干涉测角系统进行解算的方法，该方法在多重信号分类（MUSIC)算法的基础上，对协方差矩阵的估计、特征值分解、谱峰搜索等计算量较大的步骤进行了改进，提高了运算速度以及解模糊的正确率，但上述方法对系统的通道自由度和一致性要求较高。文献17针对多基线相位干涉仪的相位模糊问题,提出了一种基于参差基线解模糊算法原理的改进解模糊算法，该算法具有计算量小、实时计算以及正确概率高的优点。本文在充分借鉴上述文献思想的基础上，从风机叶片净空估计的实际应用背景出发，充分考虑工程可实现性、空间约束、稳定性以及运算复杂度几个方面，在双接收通道雷达观测的条件下提出了基于边缘约束的幅相联合约束叶片净空测量方法。本文的结构如下：第二节简单介绍了连续波雷达的信息体制以及风机叶片观测的物理模型；第三节则主要介绍了所提方法的主要思想；第四节描述了算法实现的具体步骤；第五节通过连续波雷达实际采集的叶片观测数据验证所提方法在有多径扰动和震动引起的杂波环境扰动下的叶片识别及净空估计效果；最后一节总结了本文所提方法的优势及其工程价值。其中2信号模型由于观测场景覆盖范围较小，且对于雷达系统的成本控制严格，因此叶片净空测量系统采用了线性调频连续波雷达。利用线性调频连续波的解线频调技术能够大幅降低接收系统的带宽，以较低的采样率录取回波数据，同时由于采用了连续波体制，雷达峰值功率可以大幅降低。雷达的发射信号波形如下所示s(t,tm)=a(t)e(2mf+im2)(1)其中t和tm分别代表雷达系统的全时间和慢时间,a(t)为发射信号的包络，一般在发射信号带宽内，该包络引起的幅度变化不会超过3dB，可以忽略不计。f。为雷达系统的中心频率，为线性调频信号的调频率。X0YOT2(XT,YT,ZT)111T1(XT,yT,ZT)1X111Y1z 1图1风机叶片净空测量场景示意图对于风机的叶片净空测量，采用如图1所示的地面雷达测量方式。其中以叶片机舱平面与塔筒轴线分别为YO轴和ZO轴，并根据右手原则确定初始坐标系的XO轴。而地面雷达的坐标系对应的坐标系分别为X1轴、Y1轴和Z1轴，两个坐标系之间的关系由现场的实际安装环境决定。T1和T2分别为雷达波束范围内照射到的叶片最近端和最远端，其相对于雷达的距离分别为R1和R2。假定雷达发射的信号带宽为B，考虑到距离向加窗的影响，其距离向分辨率p,为CP=(2)2B其中c为电磁波的传播速度，为加窗扩展系数。则叶片沿着雷达径向方向按照距离分辨率可以划分成如下离散的散射点集R,=R,+(i-1)p,i 1,2,.,L(3)(R,-R,)L=ceil(4)P,82ceil(）为向上取整函数。对叶片的回波进行解线频调处理之后，其回波形式如下所示a(t-T,)e(2m/g+2mfatm)(5)scenei=1其中TiB(6)T。20元ifa=(7)入T,=2R/c为叶片散射点的回波延时，,为叶片散射点相对于雷达的径向速度。由于场景中除了叶片自身的回波，还包括从天线旁瓣进人的地杂波、塔筒回波以及叶片和塔筒多径等回波。其中塔筒和地杂波回波由于位于零多普勒通道，因此对检测影响较小，但是其距离多普勒两维积分旁瓣分量和多径回波对目标的检测造成了比较强的干扰，因此需要对其进行检测剔除。雷达的总回波形式如下Stotal(t,tm）)=Ssenee(t,tm)+Selu(t,tm)+Snoise(t,tm)(8)3所提方法的实现方法根据大量实测数据的分析以及系统脉内/脉冲间的相位起伏，杂波和噪声的分布具有随机分布特征，由于噪声是高斯独立分布的，无法通过滤波器进行滤除，因此对于目标检测而言，系统必须保证信噪比。其次，旁瓣积分杂波和多径效应在距离、方位以及多通道间的相位也会具有比较大的随机性。另外，叶片的回波能量是由叶片散射点集的回波组成，其整体结构呈现近似线性的变化趋势，在一个相干积累时间之内，每一个距离分辨单位内的散射点多普勒扩展随着与雷达之间的距离呈现逐渐变大的趋势，同时，结合目标叶片散射点集与两个通道之间的固定相位差，可以利用目标在距离、速度和通道干涉相位图的包络约束特性对叶片进行检测，同时剔除噪声和杂波虚警点。假设连续波雷达系统采用双通道干涉相位体制，在相干积累时间内共发射了M组脉冲信号，感兴趣距离门总数为L。对解线频调之后的接收信号进行两维FFT处理得到距离和多普勒两维积累之后的回波。第1个距离单元、第个脉冲、第n个通道信号两维积累的处理如下Xn,m,l=FFT.(FFT(9)totaa.2其中FFTrange和FFTai分别是对目标沿着距离和方位分别进行FFT处理，range和a是距离和方位的窗函数，再得到上述距离和方位积累回波之后，可以计算两个通道之间的干涉相位，其处理方式如下=angle(x(10)结合目标在距离、速度平面以及干涉相位平面的包络约束特征，通过图像计算函数分别对距离多普勒图和干涉相位图进行包络检测，并计算其交集部分，处理过程如下clu,=cluster(x1.m,1)(11)cluz=cluster(m,t)clu=clu,ncluz(12)得到目标的置信区间clu之后，对每一个距离门内检测到的多个数据进行单元平均，进一步通过非相干积累提升信杂噪比，其距离门和干涉相位信息如下m1nn,clun.m.l(13)m,clum=1mt,clu1m，l(14)ml,clum=1由于叶片的近似线性分布,并结合其与雷达之间的几何分布模型，可以对多个连续距离的干涉相位信息进行二阶拟合处理，其操作如下=polyfit(t)(15)=polyval(R;,),i=1,2,.,L(16)在得到拟合干涉相位之后，可以根据两个通道之间与叶片仰角之间的相位关系反算得到目标叶片不同散射点集相对于雷达的仰角，其换算公式如下0,=90-asin(17)360d其中d为两个通道之间的阵元中心间距，入为发射信号的波长，asin（）为反余弦函数。结合仰角、距离以及雷达与塔筒之间的几何关系，可以进一步得到叶片与塔筒间的净空值，算法的整个流程图如图2 所示。1通道距离向脉冲基于先验的样本剔方位向FFT压缩除（保留7 0-110 m）雷达录取2通道距离向脉冲基于先验的样本剔方位向FFT原始压缩除（保留7 0-110 m）数据3通道距离向脉冲基于先验的样本剔方位向FFT压缩除（保留7 0-110 m）基于边缘目标三多通道分时能量输出叶片净约束的联维坐标俯仰干集方位测空值合CFAR虚转换涉测角角警剔除图2基于边缘约束的幅相联合叶片净空测量流程图4仿真研究为了充分验证所提方法的实际应用效果，利用某研究所开发的毫米波净空监测雷达在风力发电厂的地面端进行了实际测量实验，录取不同风况、不同朝向的数据，其现场测试环境如图3、图4所示。录取的实际测量数据涵盖了从2 r/min到12 r/min的数据样本,并且因风速的实时变化，风机在测量过程中出现了不同的朝向，其场景会出现一定程度的变化，刚好为验证本文所提方法的稳定性提供了数据多样性。因此，这组实测数83图3风电现场实际测量环境图4测量过程中风机的实时转速据可以较好的验证所提方法的性能。通道1的距离多普勒图50060-10-2070-3080-4090-50100-60110-70120-8002004006008001000多普勒频率/Hz(a)叶片距离多普勒幅度分布图两通道的干涉相位图5015060100705088-50100-100110-15012002004006008001000多普勒频率/Hz(b）叶片两通道干涉相位图图5叶片距离多普勒幅度分布图图5给出了风机处于10 r/min情况下录取的实测数据距离多普勒两维分布情况。其中图5（a）是回波的能量分布图，图5(b)是回波的两通道干涉相位图。如果对上图的结果直接进行检测处理，除了能检测到叶片的目标，还能够检测到很多扩展的杂波以及强噪声点。如果将上述所有检测点共同应用于叶片与塔筒之间的距离评估，会导致系统的测量误差增加，其目标序列对应的净空估计结果如表1中所示，出现了净空的跳变现象，不符合叶片净空的连续性现象。然而，从上图中可以看出叶片回波的幅度和干涉相位均存在比较明显的集中现象，因此可以基于上述先验信息，利用基于图像能量阈值的边缘约束方法提取叶片的边缘轮廓。另外，由于塔筒与叶片之间存在镜面/非镜面的多径效应，因此在目标的能量和干涉分布集中区域出现了干涉相位跳变的现象，这是由于多径目标的空间观测仰角发生了变化，使得干涉相位发生了变化，因此可以在边缘约束的基础上，利用幅相二维联合检测的方法，对边缘轮廓内的目标进行二次剔除，从而提升叶片目标的识别精度，为提取叶片的净空信息提供保障。通道1的距离多普勒图50-15-1660-1770-1880-19-2090-21100-22-23110-24120-2502004006008001000多普勒频率/Hz(a)幅度分布图两通道的干涉相位图5015060100705080090-50100-100110-15012002004006008001000多普勒频率/Hz(b)干涉相位图图6基于边缘约束的距离多普勒幅度分布与两通道干涉相位图图6 给出了利用能量阈值控制和边缘检测的结果，从图6(b)的干涉相位图掩膜可以明显发现,能够有效提取叶片观测仰角一直的区域，且随着距离的变化呈现连续性变化趋势。图7 在上述检测结果的基础上应用最小二乘拟合和检测点单元非相干叠加技术，在提升叶片目标干涉相位估计精8493.59392.59291.59190.59089.58988.57580859095100105110相对于雷达的径向距离(a）目标相对于雷达的仰角1514.814.614.414.21413.880859095100105相对于雷达的径向距离(b）目标相对于塔筒边缘的距离图7基于边缘约束的叶片净空测量值度的同时，得到了叶片相对于塔筒边缘的净空值。由图7（a)可以发现,由于在边缘约束范围内,无法剔除所有的非叶片目标杂波点，因此需要利用最小二程拟合技术对上述结果进行平滑处理，从而保证系统的稳定度。图7(b)给出了考虑塔筒模型和雷达间几何关系的转换后净空测量值，与激光测量仪的实际结果相差0.5m，测量结果完全符合工业应用要求。详细的净空点测量信息如表1所示，通过对比发现，采用本文所提的方法能够显著改善因杂波和强噪声环境引起的净空估计跳变问题，并且提升净空的估计精度。因此，本文所提算法相比常规的检测和干涉相位测角技术,具有更有的稳定性和测量精度，提升系统针对风机复杂工作场景下的干扰点剔除能力。另外，本文所提方法的计算复杂度较低，运算量增加可以忽略不计，因此其带来的系统增效在大多数应用场景下具备实际的工程价值。表1不同目标点相对于塔筒的距离测量值检测目标序号1234567809常规检测方法15.315.714.314.114.914.314.513.613.6本文所提方法15.315.114.914.714.614.414.214.013.9结论本文所提方法提供了一种应用于地面安装雷达测量风机叶片相对于塔筒净空的测量方法，能够有效克服叶片与塔筒多次弹射和塔筒晃动引起的杂波扩散及虚警检测问题。该方法基于叶片的整体先验信息以及图像能量阈值边界约束方法，联合目标的距离多普勒幅相二维分布图，应用掩膜检测的方法能够尽可能剔除非叶片的杂波信息，同时利用叶片的近似线性的变化特性，能够进一步剔除虚警点，最终利用最小二乘拟合和单元平均提升目标的相位估计精度，从而保证了叶片的净空估计连续性和准确性。所提算法能够在叶片高速旋转过程中塔筒震动导致的杂波扩散分布图中有效的识别出叶片分量，同时借助幅相二维检测图，能够改善在强杂波和噪声环境中的检测性能，适用于自适应叶片目标的选择与参数估计；另外所提方法应用的图像阈值约束和二维掩膜测量运算量极低，能够显著降低系统的运算复杂度，对于叶片实时监测系统的低成本化具有非常大的工程意义。参考文献：1吴道庆.地面对空情报三坐标雷达技术发展研究J.现代雷达,2 0 0 8,30(11):1-4.2周正喜，李钢.荷兰SMART系列三坐标雷达J.电子工程信息,2 0 15,(4):40-46.3Kostenko A A,Nosich A I,Tishchenko I A.Development of the firstSoviet three-coordinate L-band pulsed radar in Kharkov beforeWWIIJ.IEEE Antennas and Propagation Magazine,2001,43(3):29-48.4林晓斌，张承志，王宝，等.三坐标雷达目标高度突跳问题研究J.现代电子技术,2 0 13,11.5尹园威.三坐标DBF雷达SystemVue仿真建模方法研究J.计算机测量与控制,2 0 19,2 7(11)：197-2 0 1.6盛琥，赵温波，张远.基于转换量测高斯化的三坐标雷达跟踪算法J.电子学报,2 0 17,46（3）：56 2-56 8.7王冬杰.小型三坐标搜索警戒系统信号处理技术研究与实现D.西安电子科技大学,2 0 19.8刘伟超，成思渺，齐放，等。一种电扫描三坐标雷达测高系统设计J.现代雷达,2 0 17-1:9-10,2 1.9Liu W,Cheng S,Qi F,et al.A Design of Three-dimension RadarSystemJ.Modern Radar,2017:1.10Zhang S,Li J,Wu L,et al.Separated Maneuvering TargetTracking Algorithm Along with Three C artesian Coordinates forThree-Dimension RadarJ.Asian Journal of Control,2013,15(6):1809-1820.11王晓楠.基于相控阵雷达波束扫描的目标测角误差分析J.舰船电子对抗，2 0 18-3：7 4-7 8.12朱晓丹，朱伟强，陈卓，等。基于方向图匹配的单天线波束扫描测角方法J.航天电子对抗，2 0 17,4.13 Hanssen R F.Radar interferometry:data interpretation and erroranalysis M.Springer Science&Business Media,2001.14赵怀坤，张容权，王盛鳌.基于干涉仪的全向探测体制测角算法J.现代雷达，2 0 15-10：43-45,55.(下转第2 0 6 页）206上接第8 4页）numeration in High-Resolution Array Processing J.IEEE Trans-actions on Signal Processing,2007,55:5658-5667.9B Nadler.Nonparametric Detection of Signals by Information Theo-retic Criteria:Performance Analysis and an Improved EstimatorJ.IEEE Transactions on Signal Processing,2010,58:2746-2756.10 1L Huang,H C So.Source Enumeration Via MDL Criterion Basedon Linear Shrinkage Estimation of Noise Subspace Covariance Ma-trixJ.IEEE Transactions on Signal Processing,2013,61:4806-4821.11 M Wax,T Kailath.Detection of the number of coherent signals bythe MDL principle J.IEEE Trans.Acoust.Speech SignalProcess,1989,37:1190-1196.12M Wax.Detection and localization of multiple sources via the sto-chastic signals modelJ.IEEE Trans.Signal Process,1991,39:2450-2456.13 A Eguizabal,C Lameiro,D Ramirez,P J Schreier.Source Enumer-ation in the Presence of Colored Noise J.IEEE SignalProcessing Letters,2015,26:475-479.14谢纪岭，司锡才.基于协方差矩阵对角加载的信源数估计方法J.系统工程与电子技术，2 0 0 8：46-49.15B D Calson.Coveriance matrix estimation errors and diagonalloading in adaptive arrays J.IEEE Trans.on Aerospace and E-lectronics Systems,1988,24:397-401.16Ning Ma,Joo Thiam Goh.Efficient method to determine diagonalloading value J.Acoustics Speech and Signal Processing,2003,15:341-344.17L Huang,T Long,E Mao,H C So.MMSE-based MDL method forrobust estimation of number of sources without eigendecompositionJ.IEEE Trans.Signal Process,2009,57:4135-4142.18许佳奇，王川川，曾勇虎，汪连栋。盖尔圆定理和最小描述长度准则相结合的信源数目估计方法研究J.信号处理，2017,33(S1),19褚鼎立，陈红，蔡晓霞.基于盖尔圆准则的信源数目估计改进算法J.探测与控制学报，2 0 1 8,40（4).15 朱新国，崔鬼一种新的干涉测角数据处理算法J北京理工大學學，2 0 0 8,2 8（6）：541-544.16贾学振，孙武空间谱估计在干涉测角中的应用J.现代防御技术,2 0 1 1-4.17居易，张学成，邵文建.一种改进的多基线相位干涉仪解模糊算法J.舰船电子对抗，2 0 1 8-4.20H T Wu,J F Yang,F K Chen.Source number estimators usingtransformed Gerschgorin radiiJ.IEEE Transaction on SignalProcessing,1995,43:1325-1333.21R R Nadakuditi,J W Silverstein.Fundamental Limit of SampleGeneralized Eigenvalue Based Detection of Signals in Noise UsingRelatively Few Signal-Bearing and Noise-Only Samples J.IEEE Journal of Selected Topics in Signal Processing,2010,4:468-480.22J A Diaz-santos,K E Wage.An improved whiteningtransformation for snapshot-deficient scenarios C.OCEANS2015-MTS/IEEE Washington,Washington,DC,USA,2015:1-6.23 JJ A Diaz-santos,K E Wage.Improving whitening filter design u-sing broadband snapshots C.OCEANS 2016 MTS/IEEEMonterey,Monterey,CA,USA,2016:1-7.24 J A Diaz-santos,K E Wage.Statistical characterization of the lar-gest DMR-whitened eigenvalue for source enumeration C.OCEANS 2018 MTS/IEEE Charleston,Charleston,SC,USA,2018:1-8.25L M Johnstone.Multivariate analysis and Jacobi ensembles:Largest eigenvalue,Tracy-Widom limits and rates of convergenceJ.Ann.Statist,2008,36(6):2638-2716.26金芳晓，邱天爽，王鹏，夏楠，李景春。基于1 _1 稀疏正则化的信源个数估计新算法J.通信学报，2 0 1 6,37（1 0）.作者简介黄泽钦（1 996-），男（汉族），广东恩平人，硕士研究生，主要研究领域为信号处理、模式识别、深度学习。潘晴（1 97 5-），男（汉族），江苏宜兴人，副教授，硕士研究生导师，主要研究领域为信号处理、模式识别、深度学习（通信作者）。田妮莉（1 98 2-），女（汉族），湖北武汉人，讲师，博士，主要研究领域为图像及视频处理、时频分析。作者简介邓佳欣（1 990-），女（汉族），陕西省西安市人，工程师，主要从事雷达信号处理、目标散射特性研究和自卫对抗试飞工作。马丁峰（1 97 8-），男（汉族），山西省运城市人，高级工程师，主要研究方向为试飞技术。徐云山（1 98 8-），男（满族），辽宁省辽阳人，试飞员，主要研究方向为试飞技术。