分布式星载SAR系统时间同步和波束同步误差分析

李航舰 王 宇 邓云凯 王 伟 张 衡

①(中国科学院电子学研究所 北京 100190)

②(中国科学院大学 北京 100039)

1 引言

合成孔径雷达(Synthetic Aperture Radar,SAR)可以实现全天时、全天候的高分辨率成像[1]，星载SAR系统为全球地理遥感、地形测绘、形变监测等应用提供了丰富的高分辨数据。分布式星载SAR系统是一种特殊模式的合成孔径雷达系统，雷达发射机和接收机位于不同平台，这种收发分置的特点带来了很多优势，比如灵活的基线配置，多样的接收模式，可以提供更精细的地物分类信息等，具有单基情况无可比拟的优势。但是收发分置的特点也带来了很多新的挑战，尤其是时间同步、相位同步和波束同步等三大同步问题[2,3]。由于收发雷达不共频率源，导致了时间同步和相位同步问题[4]，时间同步问题是指，发射雷达和接收雷达使用不同的晶振产生脉冲信号，导致了触发脉冲不一致，进而导致信号发射时刻与接收窗起始时刻不一致，使距离向的录取回波减少，时间同步误差会在回波中引入相位误差，并且晶振偏差会随时间累积，对回波数据产生较大影响，严重降低雷达图像质量[5,6]。为达到时间同步的要求，Tandem-X SAR系统采用频率一致性指标为10–8量级的晶振，并采用跳跃PRI的方法，确保600 s内，时间同步误差小于5 μs[7]。相位同步问题是指，由于晶振频率的差异，导致了发射信号调制频率和接收信号解调频率不一致，造成解调回波中存在沿方位时间变化的相位误差，该误差在成像后传递到图像域，导致方位向位置偏移、主瓣展宽、积分旁瓣比恶化等问题，并且相位误差最终会传递到干涉相位，造成干涉测高精度损失。TanDEM-X SAR系统采用脉冲交替对传的同步方案解决相位同步的问题。波束同步问题是指，由于收发天线的控制精度不同，使得收发天线无法在同一时刻完全覆盖地面同一区域，导致双基合成天线方向图的增益降低，同时接收天线无法获得全部的目标散射信号，造成有效观测幅宽损失。为保证系统系能，TanDEM-X SAR系统的波束指向精度控制在0.01°，最差情况下导致的天线增益损失小于0.1 dB。

为了消除相位同步带来的影响，可以建立专门的同步链路，文献[8]提出了卫星间连续波对传的同步方法，文献[9]和文献[10]进一步提出了脉冲交替对传的方法，文献[11]提出一种基于信号交换的调频连续波合成孔径雷达相位同步方法。无论采用哪种同步链路方式，目的都是获得相位误差信息，最后在进行数据处理时提取并补偿相位误差，这也是解决分布式卫星SAR相位同步的优选方法[12]。文献[13]建立了一种时间同步误差模型，分析了时间同步误差对成像和干涉测高性能的影响，但是没有分析时间同步和相位同步的联系。时间同步和相位同步来自同一误差源，都是由于收发平台采用了独立的晶振系统来校准系统时钟[14]，所以两者之间存在一定的关系。波束同步的研究主要集中在波束同步控制方法，以及波束同步误差对雷达成像性能的影响，文献[15]分析了几种可行的飞行控制策略，同时分析了波束同步对收发平台和姿态控制精度的要求，文献[16]分析了系统波束控制精度需要满足的要求。但是，在不同的卫星姿态控制精度和波束指向精度情况下，波束同步误差对分布式星载SAR系统的综合影响还需要进一步分析。

本文结构安排如下。第2节介绍了时间同步误差模型和波束同步误差模型。第3节基于误差模型分析了时间同步和相位同步导致的相位误差之间的关系，并分析晶振稳定度需要满足的条件。第4节分析了波束同步误差对合成天线方向图、双基图像信噪比、多普勒去相干以及波束覆盖偏差的具体影响。第5节，进行仿真验证，为分布式SAR系统设计提供参考。

2 误差模型

2.1 时间同步误差模型

收发雷达之间的时间同步误差模型可以表示为常数项、线性项和随机项3部分[13]，随机项由频率源的短期稳定性决定，由于随机的特点，该项影响结果不具备积累性，一般只有在频率源稳定度较低或者所需相干积累时间较长时有较大影响[17–20]，因此该模型可以表达为：

其中，Δt为常数项，表示初始状态主辅星间雷达累积的定时信号之差；βt为线性项，由主辅星晶振频率偏差导致，β可以表示为：

其中，f0表示标称频率，Δf1和Δf2分别表示发射雷达和接收雷达的频率偏移。当存在同步误差时，在以发射信号为参考的时间坐标系中，不同脉冲的回波在接收雷达采样窗中的位置会发生偏移[21]。

2.2 波束同步误差模型

卫星姿态误差会影响卫星和地面目标的波束传播距离，同时影响回波多普勒频率，如图1所示，X轴从星下点指向卫星方向，Y轴指向方位向方向，Z轴指向距离向方向，卫星姿态误差可以分解为偏航角、俯仰角和横滚角，其中Rs表示斜距，R0表示最近斜距，P0表示姿态导引后期望的波束中心地面落点，P表示波束中心实际落点，θy表示偏航角，θp表示俯仰角，θr表示横滚角。可见，偏航角和俯仰角分别导致波束中心线绕X轴和Z轴旋转，进而偏离零多普勒面，均会导致斜视角。横滚角导致波束中心线绕Y轴旋转，所以横滚角引起波束在距离向移动，在距离向产生波束指向误差。

根据几何关系可以得到，偏航角和俯仰角引起的斜视角分别为：

其中，θl为下视角。横滚角引起的距离向波束指向误差为：

3 时间同步误差分析

主星发射的线性调频脉冲信号可以表示为：

其中，rect[·]为矩形窗函数，τ为距离向快时间，t为方位向慢时间，Tp为脉冲宽度，f0为发射信号载频，k为调频率。将时间同步误差模型代入发射信号，经辅星雷达接收解调后含有时间同步误差的回波信号为：

其中，c为光速，为波长，tc为波束中心偏离时间，Rt为雷达和卫星之间的波束传播距离，

其中，R1和R2分别表示收发最近斜距，t01和t02分别表示波束中心穿越时刻，V1和V2分别表示收发卫星的速度。为简化分析，下文中，假设两星同时穿越波束中心时刻且以该时刻为方位时间零点，即t01=t02=0。

以主星发射信号，辅星接收信号为例，采用距离多普勒聚焦算法，推导时间同步误差对SAR成像的影响，假设主星雷达晶振为标准频率源，辅星雷达晶振含有误差项，即时间同步误差全部表现在辅星获得的图像上。

3.1 固定时间误差分析

辅星雷达只存在固定时间误差的情况下，即

成像处理之后的辅星图像将存在时间同步误差项，具体可以表示为：

其中，Bd是多普勒带宽，fac是多普勒中心频率。由式(10)可以得到，固定时间同步误差在辅星图像中引入距离向主瓣偏移和固定相位误差，固定相位误差在经过干涉处理之后会产生常量干涉相位误差。主瓣偏移量和干涉相位误差分别为：

常量相位误差会影响绝对测高精度。

3.2 时间同步与相位同步关系

辅星雷达只存在线性时间误差的情况下，即

对式(7)中带有线性时间误差的相位项进一步推导可以得到

由于晶振频率存在偏移，t时刻主星的晶振频率为f1=f0+Δf1，辅星的晶振频率为f2=f0+Δf2，根据文献[9]中对脉冲对传相位同步方案的分析，辅星需要补偿的相位偏差为：

其中，fd是多普勒频率，τsys是同步脉冲的发射延迟，nφ1和nφ2分别代表晶振的相位噪声，φSNR1和φSNR2分别代表主星和辅星的接收机噪声，Δφsys(t)代表通道链路的相位噪声，Δφant(t)代表天线方向图引入的相位噪声。

其中，Ka为方位向调频率。式(16)中β(fa/Ka)项会造成距离徙动校正误差，使方位信号弥散在若干距离单元上，导致主瓣展宽。如果要求距离徙动校正误差造成的主瓣展宽小于2%，那么要求未被校正的徙动量小于半个分辨率单元[22]，即

根据时频关系可以得到

其中Tsar为合成孔径时间。

满足式(18)的情况下，经方位压缩后的SAR信号为：

此时主要导致主瓣展宽，严重情况可能导致散焦。

4 波束同步误差分析

理想情况下，SAR系统处于匀速运动状态，但是在实际工作情况下，卫星的运动受到控制系统精度的影响，会产生姿态误差，所以必须保证分布式系统的姿态控制精度[23]。波束同步误差可以分为两个方面，一方面是卫星姿态引起的误差，另一方面是波束指向误差。波束同步误差会导致天线照射范围不一致，导致信噪比损失，多普勒带宽偏移，方位向天线增益降低等问题[24]。同时，由于照射区域的偏差，还会导致距离向重叠幅宽减小。

4.1 对方位向影响

本节首先引入文献[7]中介绍的基于椭圆轨道的双星一发双收工作模式下的姿态导引方式，在此基础上，分析波束同步误差对多普勒去相干、天线增益以及双基图像信噪比损失的影响。

双星有两种姿态导引方式[7]，一种是只针对主星进行零多普勒姿态导引，此时，场景内的回波多普勒中心频率不为0。另一种是主辅双星均进行零多普勒姿态导引，此时辅星与主星照射的共同区域等效为正侧视，即场景内多普勒中心频率为0。图2和图3给出了在主辅星顺轨基线从0到3 km情况下，不同姿态导引方式下方位向合成天线方向图的情况，其它仿真参数如表1所示。

表1 仿真参数表Tab. 1 Simulation parameters

图2中，由于仅主星进行零多普勒姿态导引，导致合成后的方位天线方向图指向等效于斜视，即场景回波的多普勒中心频率不为0，产生方位多普勒去相干，导致相干系数发生变化，如图4所示。

考虑到顺轨基线和波束同步误差的综合影响，相干系数的计算方法可以表示为：

ΔfDop1表示由顺轨基线导致的多普勒频率偏移，计算公式为[7]：

其中，dalong为基线长度，r为平均斜距。

ΔfDop2表示由波束指向误差和卫星姿态误差导致的多普勒频率偏移，

图3中，由于对主辅双星均进行零多普勒姿态导引，合成后的方位天线方向图等效于正侧视，此时对多普勒去相干影响很小。考虑到主辅双星之间顺轨基线的影响，合成后的天线方向图幅度会有一定的下降，从仿真中可以看出，在顺轨基线小于3 km的情况下，导致的幅度衰减小于0.05 dB，可以忽略。在不同顺轨基线情况下，主辅星同时进行零多普勒导引的合成天线方向图对频谱加权的计算公式为[7]：

其中，l为天线长度，v为卫星速度，ΔfDop是由顺轨基线导致的多普勒频率偏移。

在此基础上，加入波束指向误差和卫星姿态误差的影响，顺轨基线和波束同步误差的综合加权公式可以表示为：

由于波束同步误差，导致成像后辅星图像信噪比降低，信噪比损失量为[25]：

其中，φ为方位向离轴角，θbw1为方位向波束宽度。

4.2 对距离向影响

针对距离向，波束控制误差主要影响主、辅星的重叠幅宽，当主、辅星距离向同时存在相反方向的指向偏差时，距离向的重叠区域(图5和图6中黄色区域)会大幅减小，影响后续的干涉处理。根据几何模型，可以推导出覆盖区域偏差的计算公式：

其中，H为卫星高度，Re为地球半径，θl为下视角，θbw2为俯仰向波束宽度，Δθr为距离向波束同步误差。

5 仿真结果

仿真参数如表1所示。

5.1 时间同步误差

仿真中加入参数为Δt=100ns的固定时间误差，点目标分析结果如图7所示，通过计算可以得到理论上的距离向位置偏移为Δr=Δt·c/2=15m，仿真得到的对应偏差为15 m，如图7(b)所示。

时间同步误差会引起主瓣展宽，如式(18)所示，仿真参数下，当β＜1.2×10-9时，引起的展宽可以忽略。当β=5×10-8时，距离向和方位向剖面图如图8所示，由仿真结果可以看到明显的主瓣展宽。

5.2 波束同步误差

本节的仿真给出了当顺轨基线为3 km时，波束指向误差和卫星姿态误差对系统性能的影响。图9和图10分别比较了方位向波束指向精度在0.01°到0.05°情况下，对合成天线方向图及相干系数的影响，仿真结果表明，与无误差情况对比，0.05°波束指向误差造成方位天线增益下降小于0.05 dB，造成相干系数减小0.02。

图11和图12分别给出了卫星姿态控制精度在0.005°到0.020°的情况下，对方位向合成天线方向图及相干系数的影响，仿真结果表明，与无误差情况对比，0.020°的姿态控制误差造成方位天线增益下降小于0.02 dB，造成相干系数减小0.01。

不同波束指向和卫星姿态误差情况下，对合成天线方向图及相干系数的综合影响结果如图13和图14所示。波束指向精度为0.05°，卫星姿态控制精度0.02°的情况下，和无误差情况对比，造成方位天线增益下降小于0.1 dB，造成相干系数减小0.03。

不同波束指向和卫星姿态误差情况下，辅星图像信噪比损失情况如表2所示，0.05°的波束指向误差和0.02°的姿态误差导致的辅星图像信噪比损失小于0.1 dB。

0°～0.06°的波束指向误差对距离向重叠幅宽的影响如图15所示。横滚角会在距离向产生波束控制误差，对距离向覆盖偏差产生影响，0.02°的姿态控制误差对重叠幅宽的影响如图16所示。横滚角从0°步进至0.02°，距离向波束指向精度从0°步进至0.06°情况下，两者对距离向重叠幅宽的综合影响如图17所示。上述结果表明，0.02°的姿态控制误差对距离向重叠幅宽的影响为368 m, 0.06°的波束指向误差对距离向重叠幅宽的影响为1100 m，由于波束指向和姿态控制导致的误差均为随机误差，所以综合影响采用几何平均数的形式，两者综合之后对距离向重叠幅宽的影响为1164 m。考虑到轨道控制精度的影响，在重轨工作模式下，两次测量中卫星的位置会在一定范围内偏移，在卫星姿态误差和波束指向误差基础上，系统设计中应综合考虑卫星位置偏移的影响。

表2 波束同步误差造成的辅星图像信噪比损失(dB)Tab. 2 Influence of beam synchronization error on SNR (dB)

由上述分析可知，在仿真精度情况下，波束同步精度对多普勒去相干、天线增益以及双基图像信噪比损失的影响较小，但对距离向重叠幅宽影响较大，因此需要在波位设计中保留一定的幅宽余量，或者根据幅宽需求，对波束同步精度提出要求。不同卫星姿态控制精度和波束指向精度对重叠幅宽的综合影响总结在表3可供系统设计参考。针对不同的卫星系统，可以根据系统参数结合式(27)求得对波束同步误差的限制。

表3 波束同步误差对重叠幅宽的影响(m)Tab. 3 Influence of beam synchronization error to overlapping area (m)

6 结论

本文分析了分布式SAR系统中时间同步误差和波束同步误差的影响，根据固定项和线性项的时间同步误差模型，推导时间同步和相位同步的关系，并分析线性项需要满足的条件。波束同步部分，本文分析了波束指向误差和卫星姿态误差对方位向合成天线方向图、多普勒去相干以及双基图像信噪比损失情况的影响，并根据几何关系推导距离向重叠幅宽的损失量。根据仿真参数，给出实验结果，以便在系统设计中充分考虑同步问题的影响，提出指标要求，为系统设计提供指导。

[1]Curlander J C and Mcdonough R N. Synthetic Aperture Radar: Systems and Signal Processing[M]. New York:Wiley, 1991, 2077: 3–7.

[2]Krieger G, Hajnsek I, Papathanassiou K P,et al..Interferometric Synthetic Aperture Radar (SAR) missions employing formation flying[J].Proceedings of the IEEE,2010, 98(5): 816–843. DOI: 10.1109/JPROC.2009.2038948.

[3]Wang Wen-qin. GPS-based time & phase synchronization processing for distributed SAR[J].IEEE Transactions on Aerospace and Electronic Systems, 2009, 45(3): 1040–1051.DOI: 10.1109/TAES.2009.5259181.

[4]Krieger G and Younis M. Impact of oscillator noise in bistatic and multistatic SAR[J].IEEE Geoscience and Remote Sensing Letters, 2006, 3(3): 424–428. DOI:10.1109/LGRS.2006.874164.

[5]Rutman J. Characterization of phase and frequency instabilities in precision frequency sources: Fifteen years of progress[J].Proceedings of the IEEE, 1978, 66(9):1048–1075. DOI: 10.1109/PROC.1978.11080.

[6]Wang W Q, Cai J Y, and Yang Y W. Extracting phase noise of microwave and millimetre-wave signals by deconvolution[J].IEE Proceedings-Science,Measurement and Technology, 2006, 153(1): 7–12. DOI: 10.1049/ipsmt:20050013.

[7]Krieger G, Moreira A, Fiedler H,et al.. TanDEM-X: A satellite formation for high-resolution SAR interferometry[J].IEEE Transactions on Geoscience and Remote Sensing, 2007, 45(11): 3317–3341. DOI:10.1109/TGRS.2007.900693.

[8]Eineder M. Ocillator clock drift compensation in bistatic interferometric SAR[C]. Proceedings of 2003 IEEE International Geoscience and Remote Sensing Symposium,Toulouse, France, 2003: 1449–1451.

[9]Younis M, Metzig R, Krieger G,et al.. Performance prediction and verification for the synchronization link of TanDEM-X[C]. Proceedings of 2007 IEEE International Geoscience and Remote Sensing Symposium, Barcelona,Spain, 2007: 5206–5209.

[10]Younis M, Metzig R, and Krieger G. Performance prediction of a phase synchronization link for bistatic SAR[J].IEEE Geoscience and Remote Sensing Letters,2006, 3(3): 429–433. DOI: 10.1109/LGRS.2006.874163.

[11]He Zhi-hua, Jin Guang-hu, Feng He,et al... Phase synchronization method for distributed spaceborne fmcw SAR system[J].Proceedings of 2016 IEEE International Geoscience and Remote Sensing Symposium,Beijing,China,2016: 1776–7379. DOI: 10.1109/IGARSS.2016.7730924.

[12]Marquez J, Mak K, Notter M,et al.. Design of passive noncooperative spaceborne SAR payloads-challenges and strategies[C]. Proceedings of the 11th European Conference on Synthetic Aperture Radar, Hamburg, Germany, 2016:1138–1143.

[13]Zhang Yong-sheng, Liang Dian-nong, and Dong Zhen.Analysis of time and frequency synchronization errors in spaceborne parasitic InSAR system[C]. Proceedings of 2006 IEEE International Conference on Geoscience and Remote Sensing Symposium, Dresden, Germany, 2006: 3047–3050.

[14]Cherniakov M. Bistatic Radars: Emerging Technology[M].New York: John Wiley & Sons, 2008.

[15]Knedlik S, Loffeld O, and Gebhardt U. On position and attitude determination requirements for future bistatic SAR experiments[C]. Proceedings of 2006 IEEE International Conference on Geoscience and Remote Sensing Symposium,Denver, Colorado, USA, 2006: 1216–1219.

[16]D’Errico M and Moccia A. Attitude and antenna pointing design of bistatic radar formations[J].IEEE Transactions on Aerospace and Electronic Systems, 2003, 39(3): 949–960.DOI: 10.1109/TAES.2003.1238748.

[17]Wang Wen-qin. Clock timing jitter analysis and compensation for bistatic synthetic aperture radar systems[J].Fluctuation and Noise Letters, 2007, 7(3):L341–L350. DOI: 10.1142/S0219477507003982.

[18]杨威, 李春升, 陈杰, 等. 时钟同步对星载多通道合成孔径雷达成像影响[J]. 电波科学学报, 2011, 26(4): 636–640. DOI: 10.13443/j.cjors.2011.04.012.Yang Wei, Li Chun-sheng, Chen Jie,et al.. Effects of time synchronization on image quality of spaceborne multichannels SAR[J].Chinese Journal of Radio Science, 2011,26(4): 636–640. DOI: 10.13443/j.cjors.2011.04.012.

[19]何志华. 分布式卫星SAR半实物仿真关键技术研究[D]. [博士论文], 国防科学技术大学, 2011.He Zhi-hua. Study on key techniques of distributed spaceborne SAR hardware-in-the-loop simulation[D]. [Ph.D.dissertation], National University of Defense Technology,2011.

[20]眭明. 星机双基地SAR同步技术研究[D]. [硕士论文], 电子科技大学, 2014.Sui Ming. Research on synchronization for spaceborne/airborne bistatic hybrid SAR[D]. [Master dissertation],University of Electronics Science and Technology of China,2014.

[21]周鹏, 皮亦鸣, 戴永寿, 等. 星机双基地SAR的时间同步误差分析[J]. 电子学报, 2011, 39(6): 1467–1470.Zhou Peng, Pi Yi-ming, Dai Yong-shou,et al.. Analysis of time synchronization errors in spaceborne/airborne hybrid bistatic SAR[J].Acta Electronica Sinica, 2011, 39(6):1467–1470.

[22]Cumming I G and Wong F H. Digital Processing of Synthetic Aperture Radar Data: Algorithms and Implementation[M].Norwood, MA, USA: Artech House, 2005.

[23]Wang Wen-qing. Inflight antenna pattern measurement for bistatic synthetic aperture radar systems[J].IEEE Antennas and Wireless Propagation Letters, 2007, 6(11): 432–435.DOI: 10.1109/LAWP.2007.905013.

[24]Geudtner D, Zink M, Gierull C,et al.. Interferometric alignment of the X-SAR antenna system on the space shuttle radar topography mission[J].IEEE Transactions on Geoscience and Remote Sensing, 2002, 40(5): 995–1006.DOI: 10.1109/TGRS.2002.1010887.

[25]张永俊, 黄海风, 张永胜, 等. 最大相干波束同步方法及其性能验证研究[J]. 电子学报, 2012, 40(8): 1564–1571. DOI: 10.3969/j.issn.0372-2112.2012.08.011.Zhang Yong-jun, Huang Hai-feng, Zhang Yong-sheng,et al..Maximum coherenc illumination synchronization method and its performance verification[J].Acta Electronica Sinica,2012, 40(8): 1564–1571. DOI: 10.3969/j.issn.0372-2112.2012.08.011.