基于二端口网络的输电杆塔雷击瞬态特性建模：理论、方法与应用

基于二端口网络的输电杆塔雷击瞬态特性建模：理论、方法与应用一、引言1.1研究背景与意义在现代社会，电力是支撑社会运转和经济发展的关键能源。输电线路作为电力系统的重要基础设施，承担着将电能从发电厂高效、可靠地传输到各个用电终端的重任，是电力系统不可或缺的“生命线”。随着经济的飞速发展和社会对电力需求的持续增长，输电线路的规模不断扩大，电压等级逐步提高，其安全稳定运行对于保障电力供应的可靠性和稳定性愈发重要。然而，输电线路通常分布广泛，大多暴露在自然环境中，不可避免地会遭受各种自然因素的威胁，其中雷击是影响输电线路安全运行的主要因素之一。雷电是一种强大的自然放电现象，其瞬间释放的巨大能量会在输电杆塔上产生极高的过电压和强大的冲击电流。当输电杆塔遭受雷击时，可能引发一系列严重后果。一方面，雷击产生的高电压可能会击穿绝缘子串，导致线路发生反击闪络，引起线路跳闸，从而造成大面积停电事故，严重影响电力系统的正常供电，给社会生产和人们的生活带来极大不便。例如，在一些工业生产中，突然的停电可能导致生产线中断，造成大量产品报废和设备损坏，带来巨大的经济损失；在医疗领域，停电可能危及患者的生命安全。另一方面，雷击产生的冲击电流会通过杆塔接地系统流入大地，若接地系统性能不佳，会使杆塔电位急剧升高，对线路设备的绝缘造成严重威胁，缩短设备的使用寿命，增加设备的维护成本。此外，雷击还可能引发火灾等次生灾害，进一步加剧对输电线路和周边环境的破坏。准确掌握输电杆塔在雷击作用下的瞬态特性，对于提高输电线路的防雷性能和保障电力系统的稳定运行具有至关重要的意义。基于二端口网络的建模方法，能够将输电杆塔抽象为具有特定端口特性的网络模型，通过对端口参数的分析和计算，可以深入研究杆塔在雷击时的电气响应特性，如电压分布、电流传输等。这有助于揭示雷击过程中杆塔内部的电磁暂态过程，为防雷措施的优化设计提供理论依据。通过建立精确的二端口网络模型，可以更加准确地评估不同防雷措施的效果，如避雷针、避雷线的布置，接地电阻的优化等，从而有针对性地改进和完善防雷方案，提高输电线路抵御雷击的能力。此外，基于二端口网络的建模研究还可以为电力系统的运行维护提供技术支持，通过对杆塔雷击瞬态特性的监测和分析，及时发现潜在的安全隐患，提前采取措施进行处理，有效降低雷击事故的发生率，保障电力系统的安全稳定运行。因此，开展基于二端口网络的输电杆塔雷击瞬态特性建模研究具有重要的理论意义和实际应用价值。1.2国内外研究现状雷击对输电线路的影响一直是电力领域的研究重点，国内外学者围绕输电杆塔雷击特性及相关建模方法开展了大量研究工作。在国外，早期研究主要聚焦于雷击参数的统计分析，如IEEE相关标准对雷电流幅值、波形等参数进行了详细的统计和规范，为后续研究提供了重要的基础数据。随着计算技术的发展，基于电磁场理论的数值计算方法逐渐应用于输电杆塔雷击特性研究。例如，有限元法（FEM）和时域有限差分法（FDTD）被用于分析杆塔在雷击时的电磁场分布和暂态响应，能够较为准确地模拟复杂的电磁环境，但计算量较大，对计算机性能要求较高。在建模方面，二端口网络模型逐渐受到关注，通过将杆塔等效为二端口网络，利用端口参数描述杆塔的电气特性，简化了分析过程。一些学者通过实验测量和理论推导相结合的方式，确定二端口网络的参数，建立了不同类型杆塔的二端口模型，并应用于雷击过电压的计算和分析。国内在输电杆塔雷击特性研究方面也取得了丰硕成果。研究人员通过对实际输电线路的运行数据进行分析，深入了解了雷击事故的发生规律和影响因素。在理论研究方面，提出了多种分析方法，如行波法，该方法能够清晰地描述雷击瞬间电流和电压在行波传播过程中的变化规律，对于理解杆塔内部的电磁暂态过程具有重要意义。在建模技术上，不断改进和完善二端口网络建模方法。运用矢量匹配法对杆塔的阻抗特性进行拟合，得到更为精确的频域有理多项式，进而建立更准确的二端口网络模型。同时，结合实际工程需求，开展了大量的仿真研究和现场试验，验证了模型的有效性和可靠性。然而，当前研究仍存在一些不足之处。一方面，现有模型在考虑复杂环境因素对杆塔雷击特性的影响时还不够全面，如地形地貌、气候条件等因素对雷击放电过程和杆塔电气参数的影响研究尚显薄弱。另一方面，模型的通用性和适应性有待提高，不同类型杆塔的模型参数确定方法缺乏统一的标准，难以快速应用于实际工程中的各种杆塔。此外，在模型验证方面，虽然进行了一些现场试验，但由于试验条件的限制，数据的完整性和准确性仍需进一步提升。本文针对当前研究的不足，旨在深入研究基于二端口网络的输电杆塔雷击瞬态特性建模方法。全面考虑复杂环境因素对杆塔雷击特性的影响，通过理论分析、数值仿真和实验研究相结合的方式，建立更为通用、准确的二端口网络模型。并通过大量的现场试验和实际案例分析，对模型进行验证和优化，为输电线路的防雷设计和运行维护提供更加可靠的理论依据和技术支持。1.3研究目标与内容本研究旨在通过深入研究，建立基于二端口网络的精确输电杆塔雷击瞬态特性模型，全面揭示输电杆塔在雷击作用下的电气响应规律，为输电线路的防雷设计与运行维护提供坚实可靠的理论依据和技术支持。围绕这一目标，具体研究内容如下：输电杆塔雷击相关理论研究：系统梳理雷电放电的物理过程和基本特性，深入分析雷电流的幅值、波形、波头时间、波尾时间等参数的变化规律及其对输电杆塔雷击特性的影响。研究雷击输电杆塔时的电磁暂态过程，包括电流和电压在行波传播过程中的变化规律，明确杆塔各部分之间的电磁耦合关系。对二端口网络理论进行深入剖析，掌握二端口网络的基本概念、参数定义和分析方法，为后续建立输电杆塔的二端口网络模型奠定理论基础。基于二端口网络的输电杆塔建模方法研究：结合输电杆塔的结构特点和电气特性，运用二端口网络理论，建立输电杆塔的等效二端口网络模型。通过对杆塔的几何尺寸、材料特性、接地条件等因素的综合考虑，确定二端口网络模型的参数，如阻抗、导纳、传输矩阵等。研究复杂环境因素对杆塔雷击特性和二端口网络模型参数的影响，考虑地形地貌（如山地、平原、丘陵等）、气候条件（如湿度、温度、风速等）对雷击放电过程和杆塔电气参数的影响，建立考虑复杂环境因素的二端口网络模型修正方法。运用矢量匹配法等先进技术，对杆塔的阻抗特性进行拟合，得到精确的频域有理多项式，提高二端口网络模型的准确性和精度。模型验证与实验研究：搭建输电杆塔雷击实验平台，采用模拟雷击装置对不同类型的输电杆塔进行雷击实验，测量杆塔在雷击作用下的电压、电流等电气参数。将实验测量数据与基于二端口网络模型的仿真计算结果进行对比分析，验证模型的准确性和可靠性。对模型中存在的误差进行分析，找出误差产生的原因，并对模型进行修正和优化。开展现场实测研究，选取实际运行的输电线路杆塔，在雷击发生时，利用现场监测设备采集杆塔的雷击数据，进一步验证模型在实际工程中的有效性和适用性。基于模型的输电杆塔雷击特性分析与应用：利用建立的二端口网络模型，深入分析输电杆塔在不同雷击条件下的瞬态特性，如电压分布、电流传输、能量损耗等，揭示雷击过程中杆塔内部的电磁暂态变化规律。根据模型分析结果，评估不同防雷措施对输电杆塔雷击特性的影响，如避雷针、避雷线的布置方式，接地电阻的大小，绝缘子的绝缘性能等，为防雷措施的优化设计提供科学依据。将模型应用于实际输电线路的防雷设计和运行维护中，通过对输电线路杆塔的雷击风险评估，提出针对性的防雷改进措施，提高输电线路的防雷性能和安全运行水平。1.4研究方法与技术路线为深入研究基于二端口网络的输电杆塔雷击瞬态特性建模，本研究将综合运用理论分析、数值仿真和实验测量三种方法，从不同角度对输电杆塔雷击特性进行全面、系统的研究，确保研究结果的科学性、准确性和可靠性。理论分析：深入研究雷电放电的物理过程、雷电流特性以及雷击输电杆塔时的电磁暂态过程，为建立基于二端口网络的输电杆塔模型提供坚实的理论基础。通过对二端口网络理论的深入剖析，掌握二端口网络的参数定义和分析方法，结合输电杆塔的结构特点和电气特性，推导二端口网络模型的参数计算方法。运用数学工具对雷击过程中的电流、电压分布进行分析，建立数学模型，揭示输电杆塔在雷击作用下的电气响应规律。数值仿真：利用专业的电磁暂态仿真软件，如ATP-EMTP、PSCAD等，建立输电杆塔的精细化数值模型。在仿真模型中，考虑杆塔的几何结构、材料特性、接地条件以及复杂环境因素（如地形地貌、气候条件等）对雷击特性的影响，通过设置不同的雷击参数（如雷电流幅值、波形、入射角等），模拟输电杆塔在各种雷击情况下的瞬态响应，获取杆塔上的电压、电流分布等数据。对仿真结果进行分析和处理，与理论分析结果相互验证，深入研究输电杆塔雷击瞬态特性的影响因素和变化规律。实验测量：搭建输电杆塔雷击实验平台，采用模拟雷击装置对不同类型的输电杆塔进行雷击实验。在实验过程中，利用高精度的测量仪器，如罗氏线圈、高压分压器等，测量杆塔在雷击作用下的电压、电流等电气参数，并记录实验数据。将实验测量数据与理论分析和数值仿真结果进行对比分析，验证理论模型和仿真模型的准确性和可靠性，找出模型中存在的误差和不足之处，为模型的修正和优化提供依据。开展现场实测研究，选取实际运行的输电线路杆塔，在雷击发生时，利用现场监测设备采集杆塔的雷击数据，进一步验证模型在实际工程中的有效性和适用性。基于上述研究方法，本研究的技术路线如图1所示：图1技术路线图首先，进行输电杆塔雷击相关理论研究，包括雷电放电特性、电磁暂态过程以及二端口网络理论等，为后续研究奠定理论基础。在理论研究的基础上，结合输电杆塔的结构和电气特性，运用二端口网络理论建立输电杆塔的等效二端口网络模型，并考虑复杂环境因素对模型参数的影响，建立模型修正方法。利用矢量匹配法等技术对杆塔的阻抗特性进行拟合，提高模型的准确性。然后，利用数值仿真软件对建立的二端口网络模型进行仿真分析，研究输电杆塔在不同雷击条件下的瞬态特性。同时，搭建实验平台进行雷击实验，对不同类型的输电杆塔进行实验测量，获取实验数据。将实验测量数据与数值仿真结果进行对比验证，对模型进行修正和优化。最后，将优化后的模型应用于实际输电线路的防雷设计和运行维护中，通过对输电线路杆塔的雷击风险评估，提出针对性的防雷改进措施，提高输电线路的防雷性能和安全运行水平。二、二端口网络理论基础2.1二端口网络的基本概念二端口网络，又称为双口网络，是指端口数n等于2的多端网络。在一个具有四个引出端的网络中，将这四个引出端划分为两对，若在任意时刻，每对引出端中流入一个端子的电流恒等于从另一个端子流出的电流，即满足i_1=i_1′、i_2=i_2′，则这样的四端子网络就构成了两对端口，被称为二端口网络。其中，接电源的端口称为输入端口，一般记为1-1′；接负载的端口称为输出端口，通常记为2-2′。端口上的电压V_1、V_2和电流i_1、i_2被统称为端口变量，它们是描述二端口网络特性的关键物理量。二端口网络具有多种类型，根据其内部是否含有独立电源，可分为无源二端口网络和有源二端口网络。无源二端口网络内部不包含独立电源，仅由电阻、电感、电容等无源元件组成，其特性主要取决于这些无源元件的参数和连接方式。而有源二端口网络内部含有独立电源，如电压源或电流源，这些电源会对网络的电气特性产生重要影响，使其分析更为复杂。按照网络元件的性质，二端口网络又可分为线性二端口网络和非线性二端口网络。线性二端口网络由线性元件构成，其端口变量之间满足线性关系，即满足叠加定理。对于线性二端口网络，当多个输入信号同时作用时，其输出响应等于各个输入信号单独作用时所产生响应的叠加。而非线性二端口网络中包含非线性元件，如二极管、晶体管等，其端口变量之间的关系是非线性的，不满足叠加定理，分析时需要采用特殊的方法。此外，根据网络特性是否随时间变化，二端口网络还可分为时不变二端口网络和时变二端口网络。时不变二端口网络的参数不随时间变化，其电气特性在不同时刻保持一致，分析相对较为简单。时变二端口网络的参数会随时间发生变化，例如某些含有变容二极管或时变电感的网络，其分析需要考虑时间因素，难度较大。在实际应用中，许多电路都可以抽象为二端口网络进行分析。以变压器为例，它是电力系统中常用的电气设备，用于实现电压的变换和电能的传输。从电路模型的角度来看，变压器可以看作是一个二端口网络，其初级绕组和次级绕组分别对应二端口网络的输入端口和输出端口。通过对变压器的二端口网络模型进行分析，可以研究其电压比、电流比、阻抗变换等特性，为电力系统的设计和运行提供重要依据。又如放大器，它是电子电路中用于放大信号的关键部件，同样可以等效为一个二端口网络。输入信号从放大器的输入端口输入，经过内部的放大电路处理后，从输出端口输出放大后的信号。通过分析放大器的二端口网络模型，可以确定其电压增益、电流增益、输入阻抗、输出阻抗等参数，从而优化放大器的性能，满足不同的应用需求。2.2二端口网络的参数与等效电路在二端口网络的研究中，为了深入分析其电气特性，常使用多种参数来描述网络的端口特性。其中，阻抗参数（Z参数）、导纳参数（Y参数）和传输参数（T参数）是最为常用的三种参数。阻抗参数（Z参数）以端口电流作为自变量，端口电压作为因变量来描述二端口网络。对于一个线性二端口网络，其端口电压U_1、U_2与电流I_1、I_2之间的关系可以表示为：\begin{cases}U_1=Z_{11}I_1+Z_{12}I_2\\U_2=Z_{21}I_1+Z_{22}I_2\end{cases}写成矩阵形式为：\begin{bmatrix}U_1\\U_2\end{bmatrix}=\begin{bmatrix}Z_{11}&Z_{12}\\Z_{21}&Z_{22}\end{bmatrix}\begin{bmatrix}I_1\\I_2\end{bmatrix}其中，Z_{11}表示端口2开路（I_2=0）时，端口1的输入阻抗，即Z_{11}=\frac{U_1}{I_1}\big|_{I_2=0}；Z_{12}表示端口1开路（I_1=0）时，端口2到端口1的转移阻抗，Z_{12}=\frac{U_1}{I_2}\big|_{I_1=0}；Z_{21}表示端口2开路时，端口1到端口2的转移阻抗，Z_{21}=\frac{U_2}{I_1}\big|_{I_2=0}；Z_{22}表示端口1开路时，端口2的输入阻抗，Z_{22}=\frac{U_2}{I_2}\big|_{I_1=0}。对于互易网络，满足Z_{12}=Z_{21}；对于对称网络，除了满足互易条件外，还满足Z_{11}=Z_{22}。导纳参数（Y参数）则是以端口电压作为自变量，端口电流作为因变量，其方程为：\begin{cases}I_1=Y_{11}U_1+Y_{12}U_2\\I_2=Y_{21}U_1+Y_{22}U_2\end{cases}矩阵形式为：\begin{bmatrix}I_1\\I_2\end{bmatrix}=\begin{bmatrix}Y_{11}&Y_{12}\\Y_{21}&Y_{22}\end{bmatrix}\begin{bmatrix}U_1\\U_2\end{bmatrix}其中，Y_{11}是端口2短路（U_2=0）时端口1的输入导纳，Y_{11}=\frac{I_1}{U_1}\big|_{U_2=0}；Y_{12}是端口1短路（U_1=0）时端口1对端口2的转移导纳，Y_{12}=\frac{I_1}{U_2}\big|_{U_1=0}；Y_{21}是端口2短路时端口2对端口1的转移导纳，Y_{21}=\frac{I_2}{U_1}\big|_{U_2=0}；Y_{22}是端口1短路时端口2的输入导纳，Y_{22}=\frac{I_2}{U_2}\big|_{U_1=0}。对于互易网络，有Y_{12}=Y_{21}；对于对称网络，Y_{11}=Y_{22}。传输参数（T参数），又称ABCD参数，用端口2的电压U_2和电流-I_2作为自变量，端口1的电压U_1和电流I_1作为因变量，其线性方程组为：\begin{cases}U_1=A_{11}U_2+A_{12}(-I_2)\\I_1=A_{21}U_2+A_{22}(-I_2)\end{cases}矩阵形式为：\begin{bmatrix}U_1\\I_1\end{bmatrix}=\begin{bmatrix}A_{11}&A_{12}\\A_{21}&A_{22}\end{bmatrix}\begin{bmatrix}U_2\\-I_2\end{bmatrix}其中，A_{11}表示T2开路时电压的转移参数，A_{11}=\frac{U_1}{U_2}\big|_{I_2=0}；A_{12}表示T2短路时转移阻抗，A_{12}=\frac{U_1}{-I_2}\big|_{U_2=0}；A_{21}表示T2开路时转移导纳，A_{21}=\frac{I_1}{U_2}\big|_{I_2=0}；A_{22}表示T2短路时电流的转移参数，A_{22}=\frac{I_1}{-I_2}\big|_{U_2=0}。对于互易网络，满足A_{11}A_{22}-A_{12}A_{21}=1；对于对称网络，A_{11}=A_{22}。这些参数之间存在着一定的转换关系。以Z参数和Y参数为例，若已知Z参数矩阵[Z]，则可通过公式[Y]=[Z]^{-1}求得Y参数矩阵[Y]，反之亦然。具体计算时，对于2\times2的矩阵[Z]=\begin{bmatrix}Z_{11}&Z_{12}\\Z_{21}&Z_{22}\end{bmatrix}，其逆矩阵[Z]^{-1}=\frac{1}{Z_{11}Z_{22}-Z_{12}Z_{21}}\begin{bmatrix}Z_{22}&-Z_{12}\\-Z_{21}&Z_{11}\end{bmatrix}，从而得到Y参数矩阵。同样，Z参数与T参数、Y参数与T参数之间也可以通过相应的公式进行转换，这些转换关系在实际分析中为根据不同的已知条件选择合适的参数进行电路分析提供了便利。根据不同的参数，可以得到二端口网络的常见等效电路形式。基于Z参数的T型等效电路，由三个阻抗元件组成，如图2所示：图2基于Z参数的T型等效电路其中，Z_{a}=Z_{11}-Z_{12}，Z_{b}=Z_{12}，Z_{c}=Z_{22}-Z_{12}。对于互易网络，由于Z_{12}=Z_{21}，该T型等效电路仅包含三个独立的阻抗元件。若二端口网络内部含有受控源，其四个参数相互独立，等效电路会更加复杂，但依然可以通过对端口电压、电流关系的分析来确定等效电路参数。基于Y参数的Π型等效电路，由三个导纳元件组成，如图3所示：图3基于Y参数的Π型等效电路其中，Y_{a}=Y_{11}+Y_{12}，Y_{b}=-Y_{12}，Y_{c}=Y_{22}+Y_{12}。对于互易网络，同样由于Y_{12}=Y_{21}，该Π型等效电路中的元件参数关系相对简单。若网络内部含有受控源，参数独立性影响等效电路的复杂度。这些等效电路在电路分析中具有重要作用。以分析输电杆塔的雷击瞬态特性为例，将输电杆塔等效为二端口网络后，利用其等效电路可以方便地计算杆塔在雷击时端口的电压、电流响应。通过对等效电路中元件参数的调整和分析，能够研究不同因素（如杆塔结构、接地电阻等）对雷击特性的影响。例如，在研究接地电阻对杆塔雷击过电压的影响时，可以通过改变等效电路中与接地相关的元件参数，模拟不同接地条件下的电路响应，从而为优化接地设计提供理论依据。2.3二端口网络在电力系统建模中的应用原理在电力系统建模领域，二端口网络理论凭借其独特的优势发挥着关键作用。其应用原理主要基于对网络端口特性的深入研究，通过特定的参数来精准描述电力系统元件或子系统的电气特性。在输电线路建模中，将其视为二端口网络时，端口1可定义为输电线路的首端，端口2为末端。端口变量中的电压V_1、V_2分别表示首端和末端的电压，电流i_1、i_2则表示首端流入和末端流出的电流。通过确定二端口网络的参数，如阻抗参数（Z参数）、导纳参数（Y参数）或传输参数（T参数），可以建立起输电线路的数学模型，从而深入分析其电气性能。以Z参数为例，对于输电线路，Z_{11}表示线路末端开路时首端的输入阻抗，它反映了首端对电源的负载特性。当研究输电线路的功率传输时，Z_{11}的值会影响电源与线路之间的功率匹配情况。若Z_{11}与电源内阻不匹配，会导致功率传输效率降低，出现功率反射等问题。Z_{12}是线路首端开路时末端到首端的转移阻抗，它体现了线路两端之间的电磁耦合关系。在分析输电线路的电压分布时，Z_{12}会影响首端电压对末端电压的传递特性。当线路遭受雷击等外部干扰时，Z_{12}会对干扰信号在两端之间的传输产生作用，进而影响线路的正常运行。同样，Z_{21}和Z_{22}也分别具有特定的物理意义，它们共同构成了描述输电线路电气特性的Z参数矩阵。基于二端口网络的输电线路建模方法具有显著优势。这种方法能够将复杂的输电线路简化为具有明确端口特性的网络模型，大大降低了分析的复杂性。相比于直接对输电线路的详细物理结构进行分析，二端口网络模型只需关注端口变量和参数，避免了对线路内部复杂电磁过程的详细求解，提高了分析效率。在计算输电线路的稳态电压和电流分布时，利用二端口网络的参数方程，可以快速准确地得到结果，为电力系统的规划和运行提供了有力支持。二端口网络模型具有良好的通用性和灵活性。不同类型、不同长度的输电线路，只需确定其相应的二端口网络参数，就可以采用统一的分析方法进行研究。在实际工程中，对于新建的输电线路或对现有线路进行改造时，都可以方便地应用二端口网络模型进行电气性能的预测和分析。与其他常见的输电线路建模方法相比，二端口网络建模方法具有明显的差异。例如，传统的电路分析法，如基尔霍夫定律法，需要对输电线路的每一个元件、每一条支路进行详细的分析和计算，建立复杂的电路方程。这种方法虽然能够精确描述线路内部的电气过程，但计算量巨大，对于大规模的输电网络，计算过程繁琐且容易出错。而二端口网络建模方法则从整体端口特性出发，简化了分析过程。在计算一条较长的输电线路时，采用基尔霍夫定律法需要对线路上众多的电阻、电感、电容元件进行逐一分析，列出大量的方程。而使用二端口网络建模，只需确定线路的二端口网络参数，通过简单的矩阵运算即可得到端口的电压、电流等信息。另一种常用的建模方法是有限元法（FEM），它通过将输电线路划分为众多微小的单元，对每个单元进行电磁分析，然后通过数值计算得到整个线路的电气特性。FEM能够精确模拟输电线路的复杂几何形状和电磁边界条件，但计算成本高，对计算机性能要求苛刻。相比之下，二端口网络建模方法不需要进行复杂的网格划分和数值迭代计算，计算速度快，更适合于工程实际应用。在分析含有复杂地形地貌的输电线路时，FEM需要精确考虑地形对电磁场的影响，进行复杂的建模和计算。而二端口网络建模方法可以通过对参数的修正来近似考虑地形等因素的影响，计算过程相对简单。三、输电杆塔雷击瞬态特性分析3.1雷击对输电杆塔的影响机制雷击是一种极具破坏力的自然现象，其发生过程涉及复杂的物理变化和电磁过程。当雷电发生时，雷云与大地之间形成一个巨大的电场，随着电场强度的不断增强，空气的绝缘性能逐渐下降，最终导致空气被击穿，形成导电通道，即雷电通道。在这个过程中，大量的电荷通过雷电通道快速释放，形成强大的雷电流。雷电流的幅值可高达数十千安甚至更高，其波形具有快速上升和缓慢下降的特点，波头时间通常在数微秒以内，波尾时间则在数十微秒到数百微秒之间。当输电杆塔遭受雷击时，雷电流会瞬间注入杆塔。由于杆塔本身具有一定的电阻和电感，雷电流在杆塔中流动时会产生电压降，导致杆塔电位急剧升高。这一过程中，杆塔的冲击阻抗起着关键作用。冲击阻抗是指在冲击电流作用下，杆塔对电流呈现的等效阻抗，它与杆塔的结构、材料以及接地条件等因素密切相关。对于高耸的铁塔结构杆塔，其电感较大，在雷击瞬间，电感会阻碍电流的快速变化，使得杆塔上的电压迅速上升。接地电阻的大小也会对杆塔的冲击阻抗产生显著影响。若接地电阻较大，雷电流在流入大地时会受到较大的阻碍，导致杆塔电位进一步升高。雷击瞬间，杆塔上会产生极高的过电压。这一过电压主要由两部分组成：一部分是雷电流通过杆塔电阻和电感产生的压降，称为电阻压降和电感压降；另一部分是由于雷电通道与输电线路之间的电磁耦合，在导线上感应出的过电压。电阻压降与雷电流幅值和杆塔电阻成正比，电感压降则与雷电流的变化率和杆塔电感成正比。当雷电流以极高的速度上升时，电感压降会非常显著。电磁耦合感应过电压的大小与雷电通道与输电线路的距离、相对位置以及雷电流的变化率等因素有关。在靠近雷电通道的一侧导线上，感应过电压可能会达到很高的数值。雷击产生的过电压和冲击电流对输电线路的安全运行构成严重威胁，可能引发多种线路故障。其中，绝缘子闪络是最为常见的故障之一。当杆塔上的过电压超过绝缘子的绝缘耐受水平时，绝缘子表面的空气会被击穿，形成导电通道，即发生闪络现象。绝缘子闪络会导致线路瞬间短路，产生强大的短路电流，可能损坏绝缘子、导线以及其他线路设备。雷击还可能导致线路跳闸。当线路发生绝缘子闪络后，如果继电保护装置检测到短路电流，会迅速动作，使线路开关跳闸，以保护电力系统的安全。频繁的线路跳闸会严重影响电力系统的供电可靠性，导致大面积停电事故，给社会生产和生活带来巨大损失。在一些极端情况下，雷击产生的强大能量还可能使导线熔断、杆塔损坏，甚至引发火灾，进一步加剧对输电线路和周边环境的破坏。3.2输电杆塔雷击瞬态特性的研究方法3.2.1理论分析方法理论分析方法是研究输电杆塔雷击瞬态特性的重要手段之一，它基于经典的电磁学理论，通过数学推导和分析来揭示雷击过程中杆塔的电气特性变化规律。在早期的研究中，人们常采用集中参数电路模型来分析输电杆塔的雷击特性。将杆塔等效为一个集中电阻和电感的串联电路，利用基尔霍夫定律来计算塔顶过电压。假设雷电流为i(t)，杆塔电阻为R，电感为L，则塔顶过电压u(t)可表示为：u(t)=Ri(t)+L\frac{di(t)}{dt}这种方法虽然简单直观，但它忽略了杆塔的分布参数特性，无法准确描述雷击过程中电流和电压在杆塔上的分布情况，计算结果与实际情况存在较大偏差。随着研究的深入，基于电磁场理论的分析方法逐渐被应用。从麦克斯韦方程出发，考虑杆塔的几何结构和材料特性，建立杆塔的电磁场模型，通过求解电磁场方程来得到杆塔的冲击阻抗等参数。对于一个简单的垂直杆塔，假设其半径为r，高度为h，周围土壤的电导率为\sigma，磁导率为\mu，根据电磁场理论，可以推导出杆塔的冲击阻抗Z的计算公式为：Z=\sqrt{\frac{j\omega\mu}{2\pi\sigma}}\ln(\frac{2h}{r})其中，\omega为角频率，j为虚数单位。该公式考虑了杆塔的几何尺寸和土壤参数对冲击阻抗的影响，相比集中参数模型，能够更准确地描述杆塔的电气特性。然而，这种基于电磁场理论的分析方法也存在一定的局限性。它通常需要对杆塔的结构进行简化假设，对于复杂的杆塔结构，如带有横担、绝缘子串等的实际杆塔，求解电磁场方程的过程非常复杂，甚至难以得到解析解。而且，在实际应用中，土壤的电导率和磁导率等参数往往难以准确测量，这也会影响计算结果的准确性。此外，该方法在考虑雷击过程中的非线性因素，如土壤电离、绝缘子闪络等方面存在困难，而这些非线性因素在实际雷击过程中对杆塔的瞬态特性有着重要影响。因此，理论分析方法虽然在揭示输电杆塔雷击瞬态特性的基本原理方面具有重要作用，但在实际应用中，往往需要结合其他方法来进行综合研究。3.2.2数值仿真计算方法数值仿真计算方法是研究输电杆塔雷击瞬态特性的重要手段，它借助计算机技术，通过建立数学模型对雷击过程进行模拟分析，能够深入研究各种因素对杆塔瞬态特性的影响。目前，常用的数值仿真软件有ATP-EMTP、PSCAD/EMTDC、COMSOLMultiphysics等。ATP-EMTP（AlternativeTransientsProgram-ElectromagneticTransientsProgram）是一款广泛应用于电力系统电磁暂态分析的软件。在输电杆塔雷击特性研究中，使用ATP-EMTP建立仿真模型的流程如下：首先，根据输电杆塔的实际结构和参数，在软件中搭建杆塔的等效电路模型，将杆塔的不同部分（如塔身、横担、绝缘子等）用相应的电路元件（电阻、电感、电容等）来表示。对于塔身，可以等效为多个集中电感和电阻的串联，电感反映了塔身的电感特性，电阻则考虑了塔身的电阻损耗。然后，设置雷击源参数，如雷电流的幅值、波形（常用双指数函数波形i(t)=I_0(1-e^{-\alphat})e^{-\betat}，其中I_0为雷电流幅值，\alpha、\beta为时间常数）、波头时间和波尾时间等。还需设置输电线路的参数，包括导线的型号、长度、弧垂，以及线路的电阻、电感、电容等参数。在设置这些参数时，需要参考相关的电力工程标准和实际线路的测量数据，以确保模型的准确性。完成模型搭建和参数设置后，进行仿真计算，得到杆塔在雷击作用下的电压、电流分布等结果。PSCAD/EMTDC（PowerSystemsComputerAidedDesign/ElectromagneticTransientsinDCSystems）也是一款功能强大的电磁暂态仿真软件。在利用PSCAD/EMTDC进行输电杆塔雷击仿真时，同样需要先建立杆塔的详细模型，该软件提供了丰富的元件库和建模工具，可以方便地构建复杂的输电杆塔模型。通过设置不同的仿真工况，如不同的雷击位置（塔顶、避雷线、导线等）、不同的接地电阻值、不同的土壤电阻率等，来研究这些因素对杆塔雷击瞬态特性的影响。在研究接地电阻对杆塔雷击过电压的影响时，可以在仿真模型中设置一系列不同的接地电阻值，从较小的数值逐渐增大，观察杆塔上不同位置的电压变化情况。通过对仿真结果的分析，可以得到接地电阻与杆塔雷击过电压之间的关系曲线，为优化接地设计提供依据。COMSOLMultiphysics是一款多物理场耦合仿真软件，它能够考虑电磁场、热场、流场等多种物理场的相互作用。在输电杆塔雷击特性研究中，COMSOLMultiphysics可以建立精确的三维电磁场模型，充分考虑杆塔的复杂几何形状和材料特性。利用该软件的有限元分析功能，将杆塔及其周围空间离散为众多微小的单元，对每个单元进行电磁场分析，然后通过数值计算得到整个杆塔在雷击时的电场和磁场分布。这种方法能够精确模拟雷击过程中电磁场的传播和变化，对于研究雷击时杆塔内部的电磁暂态过程具有重要意义。数值仿真计算方法的优点在于可以灵活地改变各种参数，模拟不同的雷击条件和工况，快速得到大量的计算结果，为研究输电杆塔雷击瞬态特性提供了丰富的数据支持。通过仿真分析，可以深入了解雷击过程中杆塔的电气响应特性，如电压分布、电流传输、能量损耗等，为防雷措施的优化设计提供理论依据。然而，数值仿真结果的准确性和可靠性依赖于模型的合理性和参数的准确性。如果模型建立不合理，忽略了一些重要的因素，或者参数设置不准确，与实际情况偏差较大，那么仿真结果就会失去参考价值。因此，在进行数值仿真计算时，需要对模型进行充分的验证和校准，将仿真结果与实际测量数据或理论分析结果进行对比，不断优化模型和参数，以提高仿真结果的准确性和可靠性。3.2.3实验测量方法实验测量方法是研究输电杆塔雷击瞬态特性的重要手段，它通过实际的实验操作，直接获取输电杆塔在雷击作用下的电气参数，为理论分析和数值仿真提供验证依据。在实验测量中，通常采用模拟雷击装置来产生与实际雷电流相似的冲击电流，对输电杆塔进行雷击实验。模拟雷击装置主要由冲击电流发生器、波形调整装置和测量系统等部分组成。冲击电流发生器是产生冲击电流的核心设备，它通过电容储能、快速放电的方式产生高幅值、短持续时间的冲击电流。波形调整装置用于调整冲击电流的波形，使其尽可能接近实际雷电流的波形，如双指数函数波形。测量系统则包括罗氏线圈、高压分压器、示波器等设备，用于测量雷击过程中杆塔上的电流和电压。罗氏线圈用于测量冲击电流，它基于电磁感应原理，将被测电流转换为感应电压，通过测量感应电压来计算电流值。高压分压器用于测量杆塔上的高电压，将高电压按一定比例降低后，输入示波器进行测量和记录。示波器则用于实时显示和记录电流和电压的波形，以便后续分析。实验测量的过程如下：首先，搭建实验平台，将输电杆塔模型安装在合适的位置，并连接好模拟雷击装置和测量系统。对输电杆塔模型进行设计和制作时，需要尽量模拟实际杆塔的结构和尺寸，采用相同的材料和连接方式，以保证实验结果的可靠性。然后，根据实验要求，设置模拟雷击装置的参数，如冲击电流的幅值、波形、波头时间和波尾时间等。在设置这些参数时，需要参考实际雷电流的统计数据和相关标准，如国际电工委员会（IEC）制定的雷电冲击标准。启动模拟雷击装置，对输电杆塔进行雷击实验，同时利用测量系统记录杆塔上的电流和电压数据。在实验过程中，需要多次重复雷击实验，以获取稳定可靠的数据，并对实验数据进行统计分析，减小测量误差。将实验测量结果与理论分析和数值仿真结果进行对比，可以发现它们之间存在一定的差异。一方面，实验测量结果可能受到测量设备精度、实验环境等因素的影响。测量设备本身存在一定的误差，如罗氏线圈的测量精度、高压分压器的分压比误差等，都会导致测量结果与实际值存在偏差。实验环境中的电磁干扰、温度变化等因素也可能对实验结果产生影响。另一方面，理论分析和数值仿真模型在建立过程中往往进行了一些简化假设，忽略了一些实际因素的影响。在理论分析中，可能忽略了杆塔的非线性特性、土壤的复杂电磁特性等；在数值仿真中，可能由于模型参数的不准确或模型结构的简化，导致仿真结果与实际情况存在差异。通过对这些差异的分析，可以进一步改进理论分析和数值仿真模型，提高对输电杆塔雷击瞬态特性的研究精度。例如，如果发现实验测量得到的杆塔冲击阻抗与理论计算值存在较大偏差，可以进一步研究杆塔的结构细节、材料特性以及土壤参数等因素对冲击阻抗的影响，对理论模型进行修正和完善。3.3现有研究方法的不足与二端口网络建模的优势在输电杆塔雷击瞬态特性研究领域，现有研究方法虽然取得了一定成果，但仍存在诸多不足。理论分析方法中，集中参数电路模型过于简化，忽略了杆塔的分布参数特性，使得计算结果与实际情况偏差较大，无法准确反映雷击时电流和电压在杆塔上的分布情况。基于电磁场理论的分析方法，虽考虑了杆塔的几何结构和材料特性，但对于复杂杆塔结构，求解电磁场方程极为复杂，往往需要进行大量简化假设，且实际中土壤参数难以精确测量，这严重影响了计算结果的准确性。此外，该方法在处理雷击过程中的非线性因素时存在较大困难，而这些非线性因素对杆塔瞬态特性有着重要影响。数值仿真计算方法中，尽管ATP-EMTP、PSCAD/EMTDC等软件能够模拟多种工况，但模型的准确性高度依赖于参数设置。实际输电杆塔的结构和电气参数复杂多变，难以精确获取，若参数设置不合理，仿真结果的可靠性将大打折扣。COMSOLMultiphysics虽能建立精确的三维电磁场模型，但计算成本高昂，对计算机性能要求苛刻，在实际工程应用中受到很大限制。实验测量方法虽能直接获取输电杆塔在雷击作用下的电气参数，但实验条件往往难以完全模拟实际情况。模拟雷击装置产生的冲击电流与实际雷电流在波形、幅值等方面可能存在差异，实验环境中的电磁干扰、温度变化等因素也会对测量结果产生影响。而且，实验测量成本较高，难以进行大规模、多工况的测试，数据的代表性有限。与现有方法相比，二端口网络建模在研究输电杆塔雷击瞬态特性方面具有显著优势。二端口网络建模能够将复杂的输电杆塔结构简化为具有明确端口特性的网络模型，仅需关注端口变量和参数，大大降低了分析的复杂性。在计算杆塔的雷击过电压时，通过确定二端口网络的参数，利用简单的矩阵运算即可得到端口的电压、电流等信息，避免了对杆塔内部复杂电磁过程的详细求解，提高了分析效率。二端口网络模型具有良好的通用性和灵活性。不同类型、不同结构的输电杆塔，只需确定其相应的二端口网络参数，就可以采用统一的分析方法进行研究。这使得该模型能够快速应用于实际工程中的各种杆塔，为工程设计和运行维护提供便利。二端口网络建模方法在考虑复杂环境因素对杆塔雷击特性的影响方面具有独特优势。通过对端口参数的修正，可以方便地考虑地形地貌、气候条件等因素对杆塔电气特性的影响。在山地地形中，由于土壤电阻率的变化和杆塔周围电场分布的改变，会对杆塔的冲击阻抗产生影响。利用二端口网络建模，可以通过调整端口参数来反映这种变化，从而更准确地研究雷击特性。四、基于二端口网络的输电杆塔雷击瞬态特性建模方法4.1矢量匹配法在杆塔阻抗特性拟合中的应用4.1.1矢量匹配法原理矢量匹配法是一种在电路分析和系统建模领域广泛应用的强大拟合技术，由B.Gustavsen提出。其核心目的是通过有理函数对给定的频域响应数据进行高精度近似，以获得简洁且准确的数学模型，从而有效简化复杂系统的分析过程。在实际应用中，许多系统的频域响应数据往往呈现出复杂的变化趋势，难以直接用简单的函数进行描述。矢量匹配法通过巧妙的数学处理，能够将这些复杂的数据拟合为有理函数形式，为后续的分析和计算提供便利。考虑一个线性时不变系统，其频域响应可以用传递函数H(s)来描述，其中s=\sigma+j\omega为复频率。矢量匹配法旨在找到一组合适的有理函数系数，使得拟合后的有理函数R(s)能够尽可能准确地逼近原传递函数H(s)。该方法的实现主要分为两个关键步骤：极点定位和留数确定。在极点定位阶段，首先需要在感兴趣的频率范围内合理分布一组初始实数极点\{p_{k0}\}。这些初始极点的选择对最终拟合结果的准确性和计算效率有着重要影响。通常，根据系统的特性和频域响应的特点，选择在关键频率点附近或响应变化剧烈的区域分布极点。引入一个随频率变化的未知变标参数\alpha(s)，构建变标函数M(s)，其表达式为：M(s)=\sum_{k=1}^{n}\frac{\alpha(s)}{s-p_{k0}}通过调整变标参数\alpha(s)，使变标函数M(s)与指定的初始极点能够精确匹配。在这个过程中，利用最小二乘法等优化算法，不断迭代求解变标参数，使得变标函数与原频域响应在给定频率点上的误差最小化。经过迭代计算，得到一组修正后的极点\{p_{k}\}，这些修正极点能够更好地反映原系统频域响应的特征。在留数确定阶段，使用修正后的极点\{p_{k}\}对未变标函数进行匹配。此时，将有理函数近似式表示为：R(s)=r+\sum_{k=1}^{n}\frac{q_{k}}{s-p_{k}}其中，r为常数项，q_{k}为对应极点p_{k}的留数。通过求解线性最小二乘方程，确定未知留数q_{k}的值。具体而言，将已知的频域响应数据代入上述有理函数近似式，构建线性方程组：\begin{bmatrix}1&\frac{1}{s_{1}-p_{1}}&\frac{1}{s_{1}-p_{2}}&\cdots&\frac{1}{s_{1}-p_{n}}\\1&\frac{1}{s_{2}-p_{1}}&\frac{1}{s_{2}-p_{2}}&\cdots&\frac{1}{s_{2}-p_{n}}\\\vdots&\vdots&\vdots&\ddots&\vdots\\1&\frac{1}{s_{m}-p_{1}}&\frac{1}{s_{m}-p_{2}}&\cdots&\frac{1}{s_{m}-p_{n}}\end{bmatrix}\begin{bmatrix}r\\q_{1}\\q_{2}\\\vdots\\q_{n}\end{bmatrix}=\begin{bmatrix}H(s_{1})\\H(s_{2})\\\vdots\\H(s_{m})\end{bmatrix}其中，s_{i}为给定的频率点，H(s_{i})为对应的频域响应值，m为频率点的数量。通过求解这个线性方程组，即可得到使有理函数R(s)最佳拟合原频域响应数据的留数q_{k}和常数项r。矢量匹配法在拟合频域函数时具有显著的优势。它能够有效地避免在有理函数逼近过程中出现的病态矩阵问题，这使得该方法在处理复杂系统的频域响应数据时更加稳定可靠。许多传统的拟合方法在面对高次有理函数拟合或频域响应数据存在噪声的情况时，容易出现矩阵奇异或数值不稳定的问题，导致拟合结果不准确甚至无法得到有效解。而矢量匹配法通过合理的极点定位和迭代求解策略，能够在保证拟合精度的同时，确保计算过程的稳定性。该方法的收敛速度较快，通常只需经过有限步骤的迭代就能得到满足要求的极点和留数，大大提高了计算效率。这使得矢量匹配法在处理大规模数据和实时性要求较高的应用场景中具有明显的优势。4.1.2杆塔阻抗特性曲线的测量与获取准确测量和获取输电杆塔的阻抗特性曲线是研究其雷击瞬态特性的关键环节，对于建立精确的二端口网络模型具有重要意义。目前，常用的测量方法主要有现场实测法和模拟实验法。现场实测法是在实际运行的输电线路杆塔上进行测量。这种方法能够直接获取真实工况下杆塔的阻抗特性，数据具有较高的真实性和可靠性。在现场实测时，通常采用冲击电流注入法。使用专门的冲击电流发生器向杆塔注入特定波形和幅值的冲击电流，同时利用罗氏线圈、高压分压器等高精度测量设备，实时测量杆塔上不同位置的电压和电流。通过测量得到的电压和电流数据，根据阻抗的定义Z=\frac{U}{I}，计算出不同频率下杆塔的阻抗值。在注入冲击电流时，需要根据实际雷电流的统计特性，选择合适的波形和幅值。常用的雷电流波形为双指数函数波形，如i(t)=I_0(1-e^{-\alphat})e^{-\betat}，其中I_0为雷电流幅值，\alpha、\beta为时间常数。通过调整冲击电流发生器的参数，使其产生的冲击电流波形和幅值尽可能接近实际雷电流。现场实测法虽然能够获取真实的杆塔阻抗特性，但也存在一些局限性。测量过程受到现场环境的影响较大，如电磁干扰、气候条件等。在实际输电线路周围，存在着各种电气设备和电磁信号，这些干扰可能会影响测量设备的准确性，导致测量数据出现误差。气候条件的变化，如温度、湿度、风速等，也会对杆塔的电气特性产生影响，进而影响测量结果的准确性。现场实测需要停电进行，这会对电力系统的正常运行造成一定的影响，增加了测量的成本和难度。模拟实验法是在实验室环境中搭建输电杆塔模型，通过模拟雷击过程来测量杆塔的阻抗特性。在搭建杆塔模型时，需要尽量模拟实际杆塔的结构和尺寸，采用相同的材料和连接方式。利用模拟雷击装置产生与实际雷电流相似的冲击电流，对杆塔模型进行雷击实验。同样使用罗氏线圈、高压分压器等测量设备，测量杆塔模型在冲击电流作用下的电压和电流，计算出阻抗值。为了保证模拟实验的准确性，还需要对实验环境进行严格控制，减少外界干扰的影响。模拟实验法的优点是可以在可控的环境下进行测量，能够方便地调整实验参数，多次重复实验，获取更丰富的数据。通过改变杆塔模型的结构、材料、接地条件等参数，研究这些因素对杆塔阻抗特性的影响。该方法不受实际输电线路运行状态的限制，不会对电力系统的正常运行造成影响。然而，模拟实验法也存在一定的不足。由于杆塔模型与实际杆塔在尺寸、材料特性等方面可能存在差异，导致测量结果与实际情况存在一定的偏差。在模拟雷击过程中，虽然尽量模拟实际雷电流的波形和幅值，但与真实的雷击过程仍存在一定的差距，这也会影响测量结果的准确性。影响杆塔阻抗特性曲线的因素众多，杆塔的结构是一个重要因素。不同类型的杆塔，如单杆、双杆、铁塔等，其几何形状、尺寸和材料分布不同，会导致杆塔的电感、电容和电阻等参数不同，从而影响阻抗特性。对于高耸的铁塔结构，其电感较大，在高频段，电感对阻抗的影响更为显著，使得杆塔的阻抗呈现出较大的感性分量。接地条件对杆塔阻抗特性也有着关键影响。接地电阻的大小直接影响雷电流的散流路径和杆塔的电位升高。若接地电阻较大，雷电流在流入大地时会受到较大阻碍，导致杆塔电位升高，从而使杆塔的阻抗增大。土壤的电阻率、电导率等特性也会影响接地效果，进而影响杆塔的阻抗特性。在土壤电阻率较高的地区，接地电阻难以降低，杆塔的阻抗特性会受到较大影响。此外，雷击电流的频率特性也会对杆塔阻抗产生影响。随着雷击电流频率的变化，杆塔的电感、电容等参数的电抗值会发生变化，导致杆塔的阻抗特性曲线发生改变。在高频段，杆塔的电容效应逐渐显现，会对阻抗特性产生重要影响。4.1.3基于矢量匹配法的有理函数拟合过程利用矢量匹配法对杆塔阻抗特性曲线进行拟合，能够得到精确的频域有理多项式，为建立准确的二端口网络模型提供关键支持。其具体拟合过程如下：首先，对测量得到的杆塔阻抗特性曲线数据进行预处理。由于实际测量过程中可能存在噪声干扰和测量误差，需要对原始数据进行滤波和去噪处理。可以采用数字滤波器，如巴特沃斯滤波器、切比雪夫滤波器等，对数据进行平滑处理，去除高频噪声。还需对数据进行归一化处理，将阻抗值和频率值映射到合适的范围内，以提高拟合算法的稳定性和精度。在极点定位阶段，根据杆塔阻抗特性曲线的特点和变化趋势，在感兴趣的频率范围内合理分布初始实数极点。对于输电杆塔的阻抗特性，通常在低频段和高频段都需要重点关注。在低频段，主要考虑杆塔的电阻和电感特性，因此可以在低频区域适当分布较多的极点。在高频段，电容效应逐渐显著，需要在高频区域合理分布极点，以准确捕捉阻抗特性的变化。引入变标参数进行迭代计算，不断调整极点位置，使变标函数与杆塔阻抗特性曲线在给定频率点上的误差最小化。在迭代过程中，可以设置合适的收敛条件，如误差小于某个阈值或迭代次数达到一定值，以确保极点定位的准确性和计算效率。完成极点定位后，进入留数确定阶段。将修正后的极点代入有理函数近似式，通过求解线性最小二乘方程来确定留数。在构建线性方程组时，需要将测量得到的杆塔阻抗特性数据代入方程。为了提高留数计算的准确性，可以采用加权最小二乘法，对不同频率点的数据赋予不同的权重。对于阻抗特性变化较为剧烈的频率点，赋予较大的权重，以突出这些点在拟合过程中的重要性。通过求解线性方程组，得到使有理函数最佳拟合杆塔阻抗特性曲线的留数和常数项。得到拟合后的有理函数后，需要对拟合结果进行评价。常用的评价指标包括均方误差（MSE）、平均绝对误差（MAE）和决定系数（R^2）等。均方误差（MSE）的计算公式为：MSE=\frac{1}{n}\sum_{i=1}^{n}(Z_{i}^{meas}-Z_{i}^{fit})^2其中，n为测量数据点的数量，Z_{i}^{meas}为第i个测量数据点的阻抗值，Z_{i}^{fit}为拟合函数在第i个数据点的阻抗值。MSE反映了拟合值与测量值之间的平均误差平方，MSE值越小，说明拟合结果越准确。平均绝对误差（MAE）的计算公式为：MAE=\frac{1}{n}\sum_{i=1}^{n}|Z_{i}^{meas}-Z_{i}^{fit}|MAE衡量了拟合值与测量值之间的平均绝对误差，它对误差的绝对值进行求和，更直观地反映了拟合结果与实际数据的偏差程度。决定系数（R^2）的计算公式为：R^2=1-\frac{\sum_{i=1}^{n}(Z_{i}^{meas}-Z_{i}^{fit})^2}{\sum_{i=1}^{n}(Z_{i}^{meas}-\overline{Z}^{meas})^2}其中，\overline{Z}^{meas}为测量数据的平均值。R^2表示拟合模型对数据的拟合优度，取值范围在0到1之间。R^2越接近1，说明拟合模型对数据的解释能力越强，拟合效果越好。在实际应用中，通过综合分析这些评价指标，可以全面评估拟合结果的准确性和可靠性。若MSE、MAE值较小，R^2值接近1，则说明拟合结果较好，有理函数能够准确地描述杆塔的阻抗特性曲线。反之，则需要调整拟合参数，重新进行拟合，以提高拟合精度。4.2基于(p+g)端口回转器的杆塔二端口网络模型构建4.2.1(p+g)端口回转器的特性与功能(p+g)端口回转器是一种特殊的多端口网络元件，在基于二端口网络的输电杆塔雷击瞬态特性建模中发挥着关键作用。它具有独特的电气特性，能够实现电压与电流之间的特殊转换关系。从端口特性来看，(p+g)端口回转器的电压和电流关系可以用特定的方程来描述。对于一个具有p个输入端口和g个输出端口的回转器，其端口电压U_{i}（i=1,2,\cdots,p）和电流I_{j}（j=1,2,\cdots,g）之间存在如下关系：U_{i}=\sum_{j=1}^{g}r_{ij}I_{j}I_{j}=\sum_{i=1}^{p}g_{ji}U_{i}其中，r_{ij}和g_{ji}分别为回转电阻和回转电导，它们决定了回转器的回转特性。这种电压与电流之间的转换关系使得回转器在电路中能够实现一些特殊的功能。在二端口网络模型中，(p+g)端口回转器主要起到阻抗变换和能量传输的作用。它可以将一个端口的阻抗特性转换为另一个端口的不同阻抗特性。在输电杆塔的二端口网络模型中，通过合理设置回转器的参数，可以将杆塔的复杂阻抗特性进行等效变换，使其更便于分析和计算。当杆塔的某一端口连接的负载阻抗与电源阻抗不匹配时，利用回转器的阻抗变换功能，可以调整端口的等效阻抗，实现更好的功率传输。回转器在能量传输过程中能够保持能量守恒，它本身并不消耗能量，只是将电能从一个端口高效地传输到另一个端口。这一特性在研究输电杆塔的雷击能量传输过程中具有重要意义，有助于深入了解雷击瞬间能量在杆塔中的分布和传输规律。例如，在雷击发生时，通过分析回转器在二端口网络中的能量传输特性，可以确定杆塔各部分吸收的能量大小，为评估杆塔的损坏程度提供依据。与传统的二端口网络元件相比，(p+g)端口回转器具有明显的优势。传统的电阻、电感、电容等元件只能实现单一的电气功能，如电阻主要用于消耗电能，电感用于储存磁场能量，电容用于储存电场能量。而回转器能够实现电压与电流的相互转换，以及阻抗的变换，功能更加多样化。在处理复杂的电路问题时，回转器可以简化电路结构，减少元件数量。在一些需要实现特定阻抗匹配的电路中，使用回转器可以避免使用多个电阻、电感和电容进行复杂的组合，从而降低电路的成本和复杂度。此外，回转器的响应速度快，能够快速适应电路中电压和电流的变化，在研究雷击这种瞬态过程时，能够更准确地反映电路的动态特性。4.2.2杆塔二端口网络模型的拓扑结构设计基于(p+g)端口回转器构建的输电杆塔二端口网络模型，其拓扑结构设计需要综合考虑输电杆塔的实际结构和电气特性。在设计过程中，首先将输电杆塔的不同部分进行合理划分，分别对应二端口网络模型中的不同元件。将杆塔的塔身等效为一个具有分布参数的电感和电阻的组合。塔身的电感反映了其在雷击瞬间对电流变化的阻碍作用，电阻则考虑了电流通过塔身时的能量损耗。根据杆塔的高度、材料以及几何形状等因素，可以确定电感和电阻的具体参数。对于一个高度为h、半径为r的圆柱形杆塔，其电感L可以通过公式L=\mu_0\frac{h}{2\pi}\ln(\frac{2h}{r})进行估算，其中\mu_0为真空磁导率。电阻R则可以根据杆塔材料的电阻率\rho、长度h和横截面积S=\pir^2，利用公式R=\rho\frac{h}{S}计算得到。将杆塔的横担等效为电阻和电容的组合。横担的电阻主要考虑其自身的电阻损耗，电容则反映了横担与周围物体之间的电容效应。横担与导线、塔身以及大地之间都存在一定的电容，这些电容在雷击过程中会对电流和电压的分布产生影响。横担与导线之间的电容C_{1}可以通过公式C_{1}=2\pi\epsilon_0\frac{l}{\ln(\frac{d}{r_w})}估算，其中\epsilon_0为真空介电常数，l为横担与导线的平行长度，d为横担与导线之间的距离，r_w为导线半径。在拓扑结构中引入(p+g)端口回转器，将塔身和横担的等效电路连接起来。回转器的端口连接方式需要根据杆塔的实际电气特性进行设计，以实现对杆塔阻抗特性的准确模拟。可以将回转器的一个端口连接到塔身等效电感和电阻的连接处，另一个端口连接到横担等效电阻和电容的连接处。通过合理设置回转器的回转电阻和回转电导参数，能够调整塔身和横担之间的电气耦合关系，使二端口网络模型更符合实际杆塔的电气特性。除了上述主要元件外，还需考虑杆塔的接地部分。将杆塔的接地装置等效为一个接地电阻R_g，它反映了雷电流流入大地时所受到的阻碍。接地电阻的大小对杆塔的雷击特性有着重要影响，通常需要根据土壤电阻率、接地极的形状和尺寸等因素进行计算。对于一个垂直打入地下的圆柱形接地极，其接地电阻R_g可以通过公式R_g=\frac{\rho}{2\pil}\ln(\frac{4l}{d})计算，其中\rho为土壤电阻率，l为接地极长度，d为接地极直径。将接地电阻连接到二端口网络模型中与塔身相连的一端，以完整地模拟杆塔的接地特性。通过以上设计，构建出的基于(p+g)端口回转器的杆塔二端口网络模型能够较为准确地反映输电杆塔的电气特性。在实际应用中，可以根据不同类型杆塔的特点，对模型中的元件参数进行调整和优化，以提高模型的准确性和适用性。对于不同高度、结构和材料的杆塔，通过合理调整塔身电感、电阻，横担电容、电阻以及回转器参数等，使模型能够更好地模拟实际杆塔在雷击时的瞬态特性。4.2.3模型参数的确定与计算确定基于(p+g)端口回转器的杆塔二端口网络模型参数是建立准确模型的关键环节，其计算方法需要综合考虑多种因素。对于塔身等效电感L和电阻R的计算，前文已给出相关公式。在实际计算时，需要准确测量杆塔的高度h、半径r以及材料的电阻率\rho。若杆塔由钢材制成，其电阻率\rho可根据钢材的具体型号从材料手册中获取。通过精确测量和计算，能够得到较为准确的塔身等效电感和电阻参数。横担等效电阻R_{t}和电容C_{t}的计算也有相应的方法。横担电阻R_{t}可根据横担材料的电阻率、长度和横截面积进行计算。横担电容C_{t}除了考虑与导线之间的电容外，还需考虑与塔身和大地之间的电容。对于横担与塔身之间的电容C_{2}，可通过类似的公式C_{2}=2\pi\epsilon_0\frac{l_{t-t}}{\ln(\frac{d_{t-t}}{r_{t}})}估算，其中l_{t-t}为横担与塔身的平行长度，d_{t-t}为横担与塔身之间的距离，r_{t}为横担的等效半径。横担与大地之间的电容C_{3}可根据横担的高度和土壤的介电常数进行估算。将这些电容进行合理的组合，得到横担的等效电容C_{t}。(p+g)端口回转器的回转电阻r和回转电导g的确定相对复杂。可以通过对实际输电杆塔进行实验测量，获取杆塔在不同频率下的阻抗特性数据。利用这些数据，结合矢量匹配法等拟合技术，对回转器的参数进行优化计算。具体而言，将测量得到的杆塔阻抗特性数据代入包含回转器参数的二端口网络模型方程中，通过最小二乘法等优化算法，不断调整回转电阻r和回转电导g的值，使得模型计算得到的阻抗特性与实际测量数据之间的误差最小化。在迭代计算过程中，可以设置合适的收敛条件，如误差小于某个阈值，以确保计算结果的准确性和计算效率。以某实际输电杆塔为例，该杆塔高度为30m，塔身采用钢材，半径为0.5m，横担长度为5m，与导线的距离为2m，与塔身的距离为1m。通过测量得到土壤电阻率为100\Omega\cdotm。根据上述公式计算得到塔身等效电感L=\mu_0\frac{30}{2\pi}\ln(\frac{2\times30}{0.5})\approx3.3\times10^{-4}H，塔身等效电阻R=\rho\frac{30}{\pi\times0.5^2}\approx38.2\Omega。横担等效电阻R_{t}根据横担材料计算得到为5\Omega。横担与导线之间的电容C_{1}=2\pi\epsilon_0\frac{5}{\ln(\frac{2}{r_w})}（假设导线半径r_w=0.01m），计算得C_{1}\approx1.1\times10^{-10}F。横担与塔身之间的电容C_{2}=2\pi\epsilon_0\frac{5}{\ln(\frac{1}{r_{t}})}（假设横担等效半径r_{t}=0.1m），计算得C_{2}\approx1.8\times10^{-10}F。通过实验测量杆塔的阻抗特性，并利用矢量匹配法进行拟合计算，得到回转器的回转电阻r=100\Omega，回转电导g=0.01S。接地电阻R_g=\frac{100}{2\pi\times3}\ln(\frac{4\times3}{0.1})\approx12.8\Omega（假设接地极长度为3m，直径为0.1m）。通过这些参数的计算和确定，建立起该输电杆塔的二端口网络模型，为后续的雷击瞬态特性分析提供了基础。4.3考虑土壤特性和接地方式的模型修正4.3.1土壤特性对杆塔雷击特性的影响分析土壤特性是影响输电杆塔雷击特性的关键因素，其中土壤电阻率和介电常数对杆塔雷击过电压和冲击阻抗有着显著的影响。土壤电阻率是表征土壤导电性能的重要参数，它与土壤的类型、含水量、温度等因素密切相关。不同类型的土壤，如砂土、黏土、壤土等，其电阻率差异较大。砂土的颗粒较大，孔隙率高，水分容易流失，导致其电阻率相对较高，一般在100-1000Ω・m之间。而黏土的颗粒细小，保水性好，电阻率相对较低，通常在10-100Ω・m范围内。土壤含水量对电阻率的影响也十分显著，随着含水量的增加，土壤中的离子浓度增大，导电性能增强，电阻率降低。在干旱地区，土壤含水量低，电阻率可高达数千Ω・m；而在湿润地区，土壤含水量高，电阻率可能降至几十Ω・m。温度的变化也会影响土壤电阻率，一般来说，温度升高，土壤中的水分蒸发，电阻率会相应增大。当土壤电阻率发生变化时，杆塔的冲击阻抗会随之改变。根据电磁场理论，杆塔的冲击阻抗与土壤电阻率的平方根成正比。当土壤电阻率增大时，杆塔的冲击阻抗也会增大，这意味着雷电流在流入大地时会受到更大的阻碍，导致杆塔电位升高。在土壤电阻率较高的山区，雷击时杆塔的冲击阻抗可能是低电阻率地区的数倍，杆塔电位升高更为明显。杆塔电位的升高会使雷击过电压增大，从而增加绝缘子闪络和线路跳闸的风险。当杆塔电位升高到一定程度，超过绝缘子的绝缘耐受水平时，绝缘子就会发生闪络，引发线路故障。土壤介电常数反映了土壤在电场作用下储存电能的能力，它与土壤的成分、结构等因素有关。不同土壤的介电常数也存在差异，一般来说，黏土的介电常数相对较高，约为10-30；砂土的介电常数较低，在3-10之间。土壤介电常数对雷击过电压的影响主要体现在电磁耦合过程中。当雷击发生时，雷电通道与输电线路之间会产生电磁耦合，在导线上感应出过电压。土壤介电常数的变化会影响电磁场的分布和传播，从而改变感应过电压的大小。在介电常数较高的土壤中，电磁场的传播速度会变慢，感应过电压的峰值可能会降低，但持续时间会延长。而在介电常数较低的土壤中，感应过电压的峰值可能会较高，但持续时间相对较短。这些变化都会对输电线路的绝缘造成不同程度的威胁。4.3.2不同接地方式下的杆塔模型调整策略常见的输电杆塔接地方式主要有水平接地极、垂直接地极和复合接地极等。不同的接地方式具有各自的特点，对杆塔的电气性能产生不同的影响，因此在建立基于二端口网络的杆塔模型时，需要根据不同的接地方式进行相应的调整。水平接地极通常采用扁钢或圆钢等材料，在地面下水平铺设。其优点是施工简单，成本较低。水平接地极的长度和截面积会影响接地电阻的大小。根据接地电阻的计算公式R=\frac{\rho}{2\piL}\ln(\frac{2L}{d})（其中\rho为土壤电阻率，L为水平接地极长度，d为接地极直径），增加水平接地极的长度可以有效降低接地电阻。但当水平接地极长度过长时，其末端的散流效果会逐渐减弱，存在一个有效长度。在实际应用中，一般根据土壤电阻率和杆塔的防雷要求来确定水平接地极的长度和截面积。在土壤电阻率较高的地区，可能需要增加水平接地极的长度或采用多根水平接地极并联的方式来降低接地电阻。在基于二端口网络的杆塔模型中，对于水平接地极接地方式，将接地极等效为一个电阻元件连接到杆塔模型的接地端。根据上述接地电阻计算公式，确定电阻元件的阻值，以准确模拟水平接地极的接地特性。垂直接地极是将接地极垂直打入地下，常见的有角钢、钢管等。垂直接地极的优点是能够深入到土壤深层，利用深层土壤的低电阻率特性，降低接地电阻。其接地电阻主要与接地极的长度、直径以及土壤电阻率有关。垂直接地极的长度越长，接地电阻越小。在土壤电阻率变化较大的地区，采用垂直接地极可以更好地适应不同土层的电阻率，提高接地效果。在杆塔模型中，对于垂直接地极接地方式，同样将接地极等效为电阻元件。根据垂直接地极的接地电阻计算公式R=\frac{\rho}{2\pil}\ln(\frac{4l}{d})（其中l为垂直接地极长度，d为接地极直径），计算电阻元件的阻值，并连接到杆塔模型的接地端。复合接地极则是将水平接地极和垂直接地极相结合，充分发挥两者的优势，以获得更好的接地效果。在一些对防雷要求较高的输电杆塔中，常采用复合接地极。复合接地极的设计需要综合考虑水平接地极和垂直接地极的参数，如长度、截面积、数量等。通过合理的设计，可以使复合接地极的接地电阻更低，且在不同土壤条件下都能保持较好的稳定性。在杆塔模型中，对于复合接地极接地方式，将水平接地极和垂直接地极分别等效为电阻元件，然后按照实际的连接方式连接到杆塔模型的接地端。根据各自的接地电阻计算公式，确定每个电阻元件的阻值，以准确模拟复合接地极的接地特性。4.3.3修正后模型的验证与分析为了验证考虑土壤特性和接地方式修正后的杆塔二端口网络模型的准确性，进行了实验验证和对比分析。在实验验证中，搭建了模拟输电杆塔实验平台。该平台包括模拟输电杆塔、模拟雷击装置、测量系统等部分。模拟输电杆塔按照实际杆塔的结构和尺寸进行设计制作，采用相同的材料和连接方式。模拟雷击装置能够产生与实际雷电流相似的冲击电流，其波形、幅值等参数可以根据实验需求进行调整。测量系统由罗氏线圈、高压分压器、示波器等设备组成，用于测量杆塔在雷击作用下的电流和电压。在实验过程中，设置了不同的土壤特性和接地方式。通过改变土壤的成分和含水量，模拟不同电阻率和介电常数的土壤。采用水平接地极、垂直接地极和复合接地极等不同的接地方式，并调整接地极的参数，如长度、截面积等。对每种实验工况进行多次雷击实验，记录杆塔上的电流和电压数据。将实验测量数据与修正后模型的仿真计算结果进行对比分析。从对比结果可以看出，修正后模型的计算结果与实验测量数据具有较