3、寄生参数提取原理：场求解器基础（2D/2.5D/3D），提取流程（几何识别、网格划分、矩阵求解）
寄生参数提取（Parasitic Extraction, PEX）是连接物理版图与电路仿真的桥梁。其核心任务是通过求解电磁场方程，将互连线的几何形状与材料属性转化为集总或分布式的电阻（R）、电容（C）和电感（L）网络。本章深入剖析场求解器的分类与提取流程，为后续校准精度优化奠定理论基础。
3.1 场求解器分类：2D / 2.5D / 3D
根据求解维度与假设条件，场求解器可分为三类。选择哪种求解器取决于互连结构的几何复杂度、频率范围以及精度要求。
| 类型 | 核心假设 | 适用场景 | 精度 | 速度 |
|---|---|---|---|---|
| 2D 场求解器 | 导体在第三维（长度方向）无限长且截面均匀 | 长直互连线、时钟树、电源网格 | 中等（忽略末端效应） | 极快 |
| 2.5D 场求解器 | 导体在垂直方向分层，水平方向可任意；基于“薄导体”假设 | 多层金属层、MIM电容、电感螺旋线圈 | 高（适用于平面工艺） | 快 |
| 3D 场求解器 | 无几何简化，全三维空间求解麦克斯韦方程 | TSV、焊球、键合线、不规则结构、高频效应 | 最高（全波精度） | 慢 |
工程经验： 在28nm及以下节点，由于互连线的宽高比和侧壁电容占比显著增加，2D求解器在耦合电容上的误差可达15%-20%。建议关键路径（如敏感模拟信号、高速时钟）强制使用2.5D或3D求解器进行校准。
3.2 提取流程：从版图到寄生网络
无论采用何种求解器，寄生参数提取的底层流程均遵循“几何识别 → 网格划分 → 矩阵求解”三大步骤。下图展示了完整的提取流水线。
3.3 步骤一：几何识别（Geometry Recognition）
几何识别是将版图数据（GDSII/OASIS）转换为求解器可理解的物理对象的过程。该阶段需要处理以下关键任务：
- 层叠映射：将版图中的层号（如M1、M2、V1）映射到工艺定义的物理厚度、介电常数（εr）和电导率（σ）。
- 图案布尔运算：对多边形进行AND/OR/NOT操作，识别出导体、介质孔和阻挡层。
- 通孔/接触孔处理：将通孔视为垂直方向的导体柱，并计算其与上下金属层的接触电阻。
- 密度均匀化：对于非关键区域，采用等效介质模型（Effective Medium）简化几何，减少计算量。
校准要点： 几何识别阶段的常见误差源包括：
• 版图与工艺的层叠厚度偏差（如ILD厚度波动±5%）。
• 通孔阵列的简化模型（如将多个小通孔合并为一个大通孔）导致电阻误差。
• 金属侧壁的粗糙度被忽略，导致高频下表面阻抗计算不准。
• 版图与工艺的层叠厚度偏差（如ILD厚度波动±5%）。
• 通孔阵列的简化模型（如将多个小通孔合并为一个大通孔）导致电阻误差。
• 金属侧壁的粗糙度被忽略，导致高频下表面阻抗计算不准。
3.4 步骤二：网格划分（Mesh Generation）
网格划分将连续的几何空间离散化为有限个单元，是求解精度与效率的平衡点。网格类型与策略直接影响矩阵的条件数。
| 网格类型 | 单元形状 | 优势 | 劣势 | 典型应用 |
|---|---|---|---|---|
| 结构化网格 | 矩形/六面体 | 生成快，内存占用低 | 难以贴合弯曲边界 | 规则互连线、电源平面 |
| 非结构化网格 | 三角形/四面体 | 可精确拟合复杂几何 | 生成慢，矩阵稀疏性差 | TSV、焊球、不规则结构 |
| 自适应网格 | 混合（h-自适应） | 在电场梯度大的区域自动加密 | 迭代次数多 | 高精度电容提取 |
网格划分的工程准则：
- 最小网格尺寸：通常设置为最小线宽的1/5 ~ 1/10。对于先进节点（≤7nm），建议网格尺寸 ≤ 0.1μm。
- 长宽比约束：三角形单元的长宽比应 ≤ 5:1，否则会导致矩阵病态。
- 边界层网格：在导体表面生成2~3层加密网格，以捕捉趋肤效应（Skin Effect）引起的电流拥挤。
3.5 步骤三：矩阵求解（Matrix Solver）
网格划分完成后，电磁场问题被转化为大型稀疏线性方程组。求解器的选择决定了提取速度与可扩展性。
3.5.1 常用数值方法
- 有限元法（FEM）：基于变分原理，将偏微分方程离散化。适用于复杂几何和多介质，但内存消耗大（O(N^1.5) ~ O(N^2)）。
- 边界元法（BEM）：仅对导体表面进行离散，将问题降维。适用于电容提取，矩阵为稠密矩阵，但未知量少。
- 矩量法（MoM）：基于积分方程，常用于3D全波电磁求解。高频下精度高，但低频时存在收敛问题。
- 快速多极子法（FMM）：加速BEM/MoM的矩阵-向量乘积，将复杂度从O(N^2)降至O(N log N)，是大规模提取的标准算法。
3.5.2 求解器收敛性检查
在实际工程中，矩阵求解的收敛性可通过以下指标监控：
| 指标 | 定义 | 可接受范围 | 不收敛时的对策 |
|---|---|---|---|
| 残差（Residual） | ||Ax - b|| / ||b|| | < 1e-6 | 增加迭代次数或改用直接求解器 |
| 矩阵条件数 | κ(A) = ||A|| · ||A⁻¹|| | < 1e4 | 重新划分网格，改善单元长宽比 |
| 能量守恒偏差 | 电容矩阵各行之和与地电容的偏差 | < 1% | 检查悬浮导体或浮空节点的处理 |
常见陷阱： 当提取结果出现负电容或负电阻时，通常不是物理现象，而是矩阵求解不收敛或网格质量差导致的数值伪影。此时应优先检查网格划分，而非直接调整SPICE模型参数。
3.6 从场解到SPICE网表：寄生网络生成
矩阵求解完成后，需要将场解结果转换为SPICE可识别的集总或分布式元件网络。转换过程遵循以下原则：
- RC网络：对于低频数字电路，通常提取π型或T型RC网络。每个互连线段被建模为一个电阻和两个对地电容。
- RLC网络：对于高频或模拟电路，需要提取自感和互感，形成RLC耦合网络。互感矩阵通常通过部分元等效电路（PEEC）方法获得。
- 分布式模型：当互连线长度超过信号上升沿有效长度的1/10时，必须使用传输线模型（W-element或RLGC模型）。
3.7 本章小结
寄生参数提取的精度取决于三个环节的协同：
- 求解器选择：根据结构复杂度在2D/2.5D/3D之间权衡，高频或3D结构必须使用全波求解器。
- 网格质量：自适应网格是平衡精度与效率的关键，建议在导体边缘和介质界面设置加密区。
- 矩阵求解收敛性：监控残差和条件数，确保数值解是物理真实的。
在后仿真校准中，提取流程的每个步骤都可能引入误差。下一阶段我们将深入探讨如何通过测试结构（Testkey）和硅片测量数据，对这些提取结果进行系统性校准。