一种算子展开的时域物理光学（TD-PO）方法用于电大色散表面的快速近场超宽带（UWB）散射

《Applied Sciences》：An Operator-Expansion TD-PO Method for Fast Near-Field UWB Scattering from Electrically Large, Dispersive Surfaces

　　

为了评估近场地面散射对超宽带（UWB）引信性能的影响，研究人员提出了一种高效的算子展开时域物理光学（OE-TD-PO）框架。该方法将传统的远场时域物理光学（TD-PO）扩展到近场激励下的电大色散粗糙表面。通过利用局部平面波近似（LPA）和时域基尔霍夫近似（KA），将复杂的散射过程分解为独立的面元级响应，从而减少了几何与波传播之间的耦合。每个面元的散射物理特性使用闭式材料和几何算子表示。材料算子负责频率相关的色散和极化反射，而几何算子则在时域内建模面元内的延迟扩展。引入瞬态偶极子辐射公式的激励阶数展开，以将源波形与空间面元循环解耦，产生对应于激励信号的导数、比例和积分项的辐射、感应和静态分量。这种解耦降低了计算复杂度，同时保持了物理保真度。经过分析和数值基准验证，该方法有效量化了地形引起的测距偏差和起爆可靠性，为跨多种环境的超宽带引信自适应误差补偿和增益控制提供了严格基础。

**论文解读：一种算子展开的时域物理光学方法用于电大色散表面的近场超宽带散射**

**研究背景与动机**

超宽带（UWB）技术凭借其高测距精度和强抗干扰能力，已广泛应用于各类地面引信系统。UWB引信通过分析发射瞬态脉冲的回波实现精确测距与地形感知，其物理保真度取决于波束与地形的电磁交互。在低俯角条件下，宽波束照射下的地面多径路径差会引入显著的时延扩展，限制测距精度。系统工作在几米至几十米的低空，对应菲涅尔（Fresnel）区域，源照明与目标尺寸相当，天线波束不可避免地覆盖数米至数十米的地面足迹。该区域内，地面呈现显著粗糙度，足迹内的路径差导致几何回波时延扩展。因此，传统的远场平面波近似不再有效。准确高效地建模大型色散粗糙表面的瞬态散射，成为提升引信可靠性与精确高度感知的关键技术挑战。

目前，电大随机粗糙表面的散射分析在微波遥感中意义重大，但面临计算效率与数值精度之间的根本权衡：宏观上，电大尺度特征对计算资源要求极高；微观上，表面随机性与几何复杂性需准确表征统计与局部特征以保持物理保真度。全波数值方法（如矩量法（MoM）、时域有限差分法（FDTD））虽可通过加速策略扩展，但随电尺寸或频率增加，内存与计算时间仍急剧增长，难以满足UWB瞬态近场散射对效率与灵活性的要求。解析与半解析模型（如先进积分方程模型（AIEM）、二阶小斜率近似（SSA-II））虽有计算优势，但多基于频域，需在多个频点重复求解，计算成本高昂。为高效刻画瞬态散射，研究者发展了系列高频时域近似方法，如时域物理光学（TD-PO）与时域等效边缘电流（TDEEC）等，并已逐渐扩展至随机粗糙表面与复杂目标的宽带瞬态分析。然而，现有闭式TD-PO方法主要用于远场激励且采用理想导体（PEC）或简化阻抗模型，难以统一建模几何传播、距离相关衰减与材料色散的多物理耦合效应，具体体现在：几何与近场建模方面，传统公式基于远场平面波假设，忽略了球形波前曲率与距离相关衰减项；材料色散与解耦方面，简化阻抗模型不足以表征超宽带宽内的频率相关菲涅尔（Fresnel）响应，且缺乏闭式时域模型将材料色散与几何效应解耦。

为此，本文提出了一种高效的算子展开时域物理光学（OE-TD-PO）框架，用于电大色散粗糙表面的瞬态近场散射。通过将散射体表征为激励无关的系统响应，实现单卷积高保真回波合成。核心贡献在于：将闭式TD-PO从远场PEC扩展到近场色散介质，融合局部平面波近似（LPA）与时域基尔霍夫近似（KA）推导面元级等效电流；建立半解析框架解耦几何与材料响应；引入基于算子分解的激励阶数展开，实现高效的两阶段计算。论文发表在《Applied Sciences》。

**主要关键技术方法**

研究人员采用了以下关键技术方法：①局部平面波近似（LPA）与时域基尔霍夫近似（KA），将复杂散射过程分解为独立的面元级响应，减少几何与波传播耦合；②闭式几何算子与材料算子，几何算子通过Gordon积分公式将面元面积分解析转化为线积分，材料算子利用时域菲涅尔（Fresnel）反射系数描述频率相关色散与极化反射；③激励阶数展开，将瞬态偶极子辐射公式按激励信号的导数、比例与积分项分解为辐射、感应与静态分量，从而将源波形与空间面元循环解耦；④解析分数延迟补偿（AFDC）技术，在面元响应构建中消除离散采样子采样时延误差，确保数值鲁棒性。

**研究结果**

**4.1 Numerical Accuracy and LPA Validity（数值精度与LPA有效性）**
通过对单一面元在不同面元尺寸与源-面元距离下的相位误差分布进行二维扫描，验证了LPA的有效边界。结果表明，当离散化准则设定为d_max/λ_2dB ≤ 1/10时，绝对相位误差严格位于安全区域内，为框架正确性提供理论支撑。

**4.2 Scattering Fidelity for Canonical Targets（标准目标散射保真度）**
针对1.0 m × 1.0 m的PEC平板在3 m、20 m和100 m距离下的宽带雷达散射截面（RCS）比较：在深近场，框架准确捕捉球形波前相位曲率引起的振荡干涉模式；随距离增加，结果平滑收敛至理论远场特性。对色散介质（采用一阶德拜（Debye）模型）的宽带反射系数验证表明，模拟结果与理论菲涅尔反射系数在整个0.5–6 GHz带宽内高度一致，精确再现了不同湿度土壤的反射强度变化与频率依赖衰减。

**4.3 Analysis of Numerical Convergence and Efficiency（数值收敛性与效率分析）**
对比四种延迟补偿方案（最近邻插值NNI、样条插值CSI、频域相移FDPS与提出的AFDC），AFDC方法实现了高保真脉冲重建，相对误差控制在1%以内，而NNI的峰值误差超过5.5%。在计算复杂度与内存消耗方面，测量运行时间随面元数N呈严格线性增长，验证了算子分离带来的O(N)复杂度。与传统直接TD-PO及近场频域PO（NF-FD-PO）相比，本方法在宽带问题中显著降低计算时间，打破了大规模建模的效率瓶颈。

**4.4 Sensing Sensitivity Analysis for Random Rough Surfaces（随机粗糙表面感知灵敏度分析）**
在55°波束宽度下分析表面均方根（RMS）高度δ_h对回波峰值幅度的影响。结果显示：当δ_h较小时（0.1–0.3 λ），低空峰值幅度稳定，足迹扩展的几何增益补偿距离衰减；随着δ_h增大，散射机制从相干反射转为非相干漫散射，峰值幅度快速衰减并趋于理论自由空间1/R²规律。在二维特征空间中，回波总能量E_total与均方根脉宽τ_rms构成近似正交映射：介电常数主要调制总能量，表面粗糙度主要控制脉冲展宽，支持近场UWB散射的多参数反演。

**5.1 Ranging Accuracy and Uncertainty（测距精度与不确定性）**
采用能量重心（EC）准则估计回波到达时间。结果表明，测距偏差ΔR与统计不确定度σ_R均随表面RMS高度δ_h单调递增：当δ_h=0.6 λ时，系统偏差达14 cm。这为地形自适应延迟补偿提供了定量依据。

**5.2 Triggering Reliability Assessment（触发可靠性评估）**
通过触发概率P_trig评估UWB引信在不同地面类型（水面、土壤、沙地）下的起爆可靠性。触发概率受介电常数ε_r与表面粗糙度δ_h共同调制：沙地表面（ε_r=4）的P_trig从0.25（δ_h=0.2 λ）降至近零（δ_h=1.4 λ），识别出临界感知失效区域，为自适应触发逻辑提供了量化基础。

**讨论与结论**

论文提出了并全面验证了一种高效的OE-TD-PO仿真方法，突破了同时实现高精度与计算效率在建模近场UWB信号与电大色散粗糙表面交互中的固有折衷。基于LPA与时域KA，该方法实现了对近场球形波前曲率与材料色散效应的高保真建模。主要贡献与发现总结如下：

1. 通过将时域辐射积分重写为算子分离形式，长持续时间激励信号对应的卷积操作完全与面元级循环解耦，将计算瓶颈转移，显著降低计算复杂度。
2. 在面元响应构建中融入解析分数延迟补偿（AFDC）技术，有效消除离散采样引起的子采样时延误差，使归一化均方根误差（NRMSE）始终保持在1%以下。
3. 大量数值实验证明了框架对关键物理现象的准确表征，包括近-远场过渡行为与介质色散。对于涉及数百万面元的大尺度场景，该算法在仅CPU计算条件下实现了实测加速比，凸显其在工程应用中的潜力。
4. 框架通过预测地形引起的测距偏差与统计起爆可靠性，实现了传感器性能的定量评估，为开发UWB引信的自适应误差补偿与增益控制策略提供了坚实的理论基础。

由于其模块化的算子架构，该框架可自然扩展融入时域迭代物理光学（TD-IPO）以处理多次散射效应。除具体UWB引信应用外，算子展开策略为其他领域中几何介导波现象与宽带瞬态的建模提供了概念上的参考。未来工作将探索GPU加速与分布式内存架构的扩展应用。