CN112698115A - 一种多层复合介质暂态过程中的电场分布确定方法及系统 - Google Patents

一种多层复合介质暂态过程中的电场分布确定方法及系统 Download PDF

Info

Publication number
CN112698115A
CN112698115A CN202011405249.2A CN202011405249A CN112698115A CN 112698115 A CN112698115 A CN 112698115A CN 202011405249 A CN202011405249 A CN 202011405249A CN 112698115 A CN112698115 A CN 112698115A
Authority
CN
China
Prior art keywords
equivalent
medium
branch
layer
fitting
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011405249.2A
Other languages
English (en)
Other versions
CN112698115B (zh
Inventor
黄猛
吴延宇
宋翰林
吕玉珍
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
Original Assignee
North China Electric Power University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by North China Electric Power University filed Critical North China Electric Power University
Priority to CN202011405249.2A priority Critical patent/CN112698115B/zh
Publication of CN112698115A publication Critical patent/CN112698115A/zh
Application granted granted Critical
Publication of CN112698115B publication Critical patent/CN112698115B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/12—Measuring electrostatic fields or voltage-potential
    • G01R29/14—Measuring field distribution

Abstract

本发明公开了一种多层复合介质暂态过程中的电场分布确定方法及系统，所述方法具体包括：分别测量多层复合介质的每一层介质的介电频谱和极化电流；对第j层介质的介电频谱进行拟合，获得第j层介质的n个第一等效支路；对第j层介质的极化电流进行拟合，获得第j层介质的m个第二等效支路；根据第j层介质的第一等效支路和第二等效支路，建立第j层介质的单一介质等效模型；将每一层介质的单一介质等效模型进行串联连接，获得复合介质等效模型；基于复合介质等效模型获得多层复合介质的电场分布。本发明通过建立介质的等效模型，并结合电路分析方法即可快速算出多层介质暂态过程中的电场分布，避免了繁琐的数值方法，极大降低了计算量和计算时间。

Description

一种多层复合介质暂态过程中的电场分布确定方法及系统
技术领域
本发明涉及电气设备绝缘测试技术领域，特别涉及一种多层复合介质暂态过程中的电场分布确定方法及系统。
背景技术
为了绝缘配合、发挥不同介质的特性、设备需求等原因，实际中多层绝缘介质经常用到，例如变压器的油纸绝缘、电缆接头、穿墙套管等。对于绝缘安全而言，每层介质中的电场需要合理，否则可能造成设备放电，引发安全事故。
随着新能源大规模接入、电力电子化设备大量应用，设备上的暂态电压较多，当前对于多层介质电场的计算方法，通常交流电压下电场按介电常数计算，直流电压下按电阻率计算，冲击电压同样按介电常数计算，但不适用于一般变化缓慢的电压计算。其它暂态电压也没有通用的计算方法，一般是通过有限元仿真、电荷运动微观方程计算，过程复杂繁琐、耗时较长，且不能计算复杂结构，不适合电气设备绝缘结构电场分布计算。
如何提供一种快速的电气设备绝缘结构电场分布计算方法，成为一个亟待解决的技术问题。
发明内容
本发明的目的是提供一种多层复合介质暂态过程中的电场分布确定方法及系统，以提供一种快速的电气设备绝缘结构电场分布计算方法。
为实现上述目的，本发明提供了如下方案：
一种多层复合介质暂态过程中的电场分布确定方法，所述确定方法包括如下步骤：
分别测量多层复合介质的每一层介质的介电频谱和极化电流，根据介电频谱的工频复电容实部得到每一层介质的介质等效电容，根据极化电流稳态值得到每一层介质的介质等效电阻；所述介电频谱为：在不同频率的交流电压激励下的频域复介电常数；
对第j层介质的介电频谱进行拟合，获得第j层介质的n个第一等效支路的第一支路等效电阻和第一支路等效电容；所述第一等效支路为第一支路等效电阻和第一支路等效电容串联连接形成的支路；
对第j层介质的极化电流进行拟合，获得第j层介质的m个第二等效支路的第二支路等效电阻和第二支路等效电容；所述第二等效支路为第二支路等效电阻和第二支路等效电容串联连接形成的支路；
根据第j层介质的极化电流的时间变化曲线确定从阶跃电压施加至极化电流稳定之间的时间作为第j层介质的极化电流时间；
判断第j层介质的极化电流时间是否满足公式td＜2/fl，获得第一判断结果；其中，td表示极化电流时间，fl表示介电频谱能够测量的最低频率；
若所述第一判断结果表示是，则将第j层介质的介质等效电阻和介质等效电容以及n个第一等效支路并联连接，形成第j层介质的单一介质等效模型；
若所述第一判断结果表示否，则将第j层介质的介质等效电阻和介质等效电容以及n个第一等效支路和m个第二等效支路并联连接，形成第j层介质的单一介质等效模型；
将每一层介质的单一介质等效模型进行串联连接，并将每一层介质的介质等效电阻、介质等效电容、第一支路等效电阻、第一支路等效电容、第二支路等效电阻和第二支路等效电容换算成介电常数和电导率参数，获得复合介质等效模型；
对复合介质等效模型施加暂态电压，计算得到复合介质等效模型输出的电场的分布，作为多层复合介质的电场分布。
可选的，所述分别测量多层复合介质的每一层介质的介电频谱和极化电流，根据介电频谱的工频复电容实部得到每一层介质的介质等效电容，根据极化电流稳态值得到每一层介质的介质等效电阻，具体包括：
测量每一层介质在不同频率的交流电压激励下的频域复介电常数；
测量每一层介质在阶跃直流电压激励下的极化电流。
可选的，所述对第j层介质的介电频谱进行拟合，获得第j层介质的n个第一等效支路的第一支路等效电阻和第一支路等效电容，具体包括：
根据所述介电频谱和电容随频率变化的关系函数
Figure BDA0002813879000000031
确定不同频率下的电容的测量值；
其中，C(ω)为频率ω下的复电容的测量值，ε(ω)表示频率ω下的复介电常数，S为第j层介质的横截面积，d为第j层介质的厚度；
令n的数值等于1；
根据不同频率下的电容的测量值，利用公式
Figure BDA0002813879000000032
对电容与介质的等效电容的差值，进行函数拟合，获得n个支路的拟合函数Nn(ω)；
其中，Ci为第i个第一等效支路的等效电容值，Ri表示第i个第一等效支路的等效电阻值；
选取K个频率点，利用公式
Figure BDA0002813879000000033
计算每个频率点的函数拟合误差值；
其中，Nn(ωk)表示第k个频率点的n个支路的拟合函数值，C'(ωk)表示第k个频率点下的复电容的实部，Co表示第j层介质工频下的复电容实部，ek表示第k个频率点的函数拟合误差值；
判断N个频率点的函数拟合误差值的平均值是否小于或等于1％，获得第二判断结果；
若所述第二判断结果表示否，则令n的数值增加1，返回步骤“根据不同频率下的电容的测量值，利用公式
Figure BDA0002813879000000034
对电容与介质的等效电容的差值，进行函数拟合，获得n个支路的拟合函数Nn(ω)”；
若所述第二判断结果表示是，输出n个支路的拟合函数的每个第一等效支路的等效电阻和等效电容。
可选的，所述对第j层介质的极化电流进行拟合，获得第j层介质的m个第二等效支路的第二支路等效电阻和第二支路等效电容，具体包括：
令h的数值等于1；
利用公式
Figure BDA0002813879000000041
对所述极化电流与极化电流最终达到的稳态值的差值，进行函数拟合，获得h个支路的拟合函数Mh(t)；
其中，Ar表示h个支路的拟合函数Mm(t)的第r个等效支路的拟合函数的系数，τr表示h个支路的拟合函数Mm(t)的第r等效支路的拟合函数的时间常数；
选取L个时间点，利用公式
Figure BDA0002813879000000042
计算每个时间点的函数拟合误差值；
其中，el表示第l个时间点的函数拟合误差值，Mh(tl)表示第l个时间点的h个支路的拟合函数值，id(tl)表示第l个时间点的极化电流值，IR表示极化电流最终达到的稳态值；
判断L个时间点的函数拟合误差值的平均值是否小于或等于1％，获得第三判断结果；
若所述第三判断结果表示否，则令n的数值增加1，返回步骤“利用公式
Figure BDA0002813879000000043
对所述极化电流与极化电流最终达到的稳态值的差值，进行函数拟合，获得h个支路的拟合函数Mh(t)”；
若所述第三判断结果表示是，则利用公式
Figure BDA0002813879000000044
和τr＝RrCr，计算h个支路的拟合函数的每个等效支路的支路等效电阻值和支路等效电容值；
其中，Uo表示测量极化电流时所加的阶跃电压幅值，Rr和Cr分别表示第r个等效支路的等效电阻和等效电容。
将τr从小到大排列，依次从中取出满足τr≥2/fl的m个支路构成第二等效支路，相应的Rr和Cr分别表示第二等效支路等效电阻和第二等效支路等效电容。
一种多层复合介质暂态过程中的电场分布确定系统，所述确定系统包括：
测量模块，用于分别测量多层复合介质的每一层介质的介电频谱和极化电流，根据介电频谱的工频复电容实部得到每一层介质的介质等效电容，根据极化电流稳态值得到每一层介质的介质等效电阻；所述介电频谱为：在不同频率的交流电压激励下的频域复介电常数；
介电频谱拟合模块，用于对第j层介质的介电频谱进行拟合，获得第j层介质的n个第一等效支路的第一支路等效电阻和第一支路等效电容；所述第一等效支路为第一支路等效电阻和第一支路等效电容串联连接形成的支路；
极化电流拟合模块，用于对第j层介质的极化电流进行拟合，获得第j层介质的m个第二等效支路的第二支路等效电阻和第二支路等效电容；所述第二等效支路为第二支路等效电阻和第二支路等效电容串联连接形成的支路；
极化电流时间确定模块，用于根据第j层介质的极化电流的时间变化曲线确定从阶跃电压施加至极化电流稳定之间的时间作为第j层介质的极化电流时间；
第一判断模块，用于判断第j层介质的极化电流时间是否满足公式td＜2/fl，获得第一判断结果；其中，td表示极化电流时间，fl表示介电频谱能够测量的最低频率；
第一单一介质等效模型建立模块，用于若所述第一判断结果表示是，则将第j层介质的介质等效电阻和介质等效电容以及n个第一等效支路并联连接，形成第j层介质的单一介质等效模型；
第二单一介质等效模型建立模块，用于若所述第一判断结果表示否，则将第j层介质的介质等效电阻和介质等效电容以及n个第一等效支路和m个第二等效支路并联连接，形成第j层介质的单一介质等效模型；
复合介质等效模型建立模块，用于将每一层介质的单一介质等效模型进行串联连接，并将每一层介质的介质等效电阻、介质等效电容、第一支路等效电阻、第一支路等效电容、第二支路等效电阻和第二支路等效电容换算成介电常数和电导率参数，获得复合介质等效模型；
电场分布计算模块，用于对复合介质等效模型施加暂态电压，计算得到复合介质等效模型输出的电场的分布，作为多层复合介质的电场分布。
可选的，所述测量模块，具体包括：
频域复介电常数测量子模块，用于测量每一层介质在不同频率的交流电压激励下的频域复介电常数；
极化电流测量子模块，用于测量每一层介质在阶跃直流电压激励下的极化电流。
可选的，所述介电频谱拟合模块，具体包括：
电容测量值确定子模块，用于根据所述介电频谱和电容随频率变化的关系函数
Figure BDA0002813879000000061
确定不同频率下的电容的测量值；
其中，C(ω)为频率ω下的复电容的测量值，ε(ω)表示频率ω下的复介电常数，S为第j层介质的横截面积，d为第j层介质的厚度；
第一初始化子模块，用于令n的数值等于1；
函数拟合子模块，用于根据不同频率下的电容的测量值，利用公式
Figure BDA0002813879000000062
对电容与介质的等效电容的差值，进行函数拟合，获得n个支路的拟合函数Nn(ω)；
其中，Ci为第i个第一等效支路的等效电容值，Ri表示第i个第一等效支路的等效电阻值；
函数拟合误差值计算子模块，用于选取K个频率点，利用公式
Figure BDA0002813879000000063
计算每个频率点的函数拟合误差值；
其中，Nn(ωk)表示第k个频率点的n个支路的拟合函数值，C'(ωk)表示第k个频率点下的复电容的实部，Co表示第j层介质工频下的复电容实部，ek表示第k个频率点的函数拟合误差值；
第二判断子模块，用于判断N个频率点的函数拟合误差值的平均值是否小于或等于1％，获得第二判断结果；
第一返回子模块，用于若所述第二判断结果表示否，则令n的数值增加1，返回步骤“根据不同频率下的电容的测量值，利用公式
Figure BDA0002813879000000071
对电容与介质的等效电容的差值，进行函数拟合，获得n个支路的拟合函数Nn(ω)”；
第一拟合结果输出子模块，用于若所述第二判断结果表示是，输出n个支路的拟合函数的每个第一等效支路的等效电阻和等效电容。
可选的，所述极化电流拟合模块，具体包括：
第二初始化子模块，用于令h的数值等于1；
第二函数拟合子模块，用于利用公式
Figure BDA0002813879000000072
对所述极化电流与极化电流最终达到的稳态值的差值，进行函数拟合，获得h个支路的拟合函数Mh(t)；
其中，Ar表示h个支路的拟合函数Mh(t)的第r个等效支路的拟合函数的系数，τr表示h个支路的拟合函数Mh(t)的第r个等效支路的拟合函数的时间常数；
第二函数拟合误差计算子模块，用于选取L个时间点，利用公式
Figure BDA0002813879000000073
计算每个时间点的函数拟合误差值；
其中，el表示第l个时间点的函数拟合误差值，Mh(tl)表示第l个时间点的h个支路的拟合函数值，id(tl)表示第l个时间点的极化电流值，IR表示极化电流最终达到的稳态值；
第三判断子模块，用于判断L个时间点的函数拟合误差值的平均值是否小于或等于1％，获得第三判断结果；
第二返回子模块，用于若所述第三判断结果表示否，则令h的数值增加1，返回步骤“利用公式
Figure BDA0002813879000000081
对所述极化电流与极化电流最终达到的稳态值的差值，进行函数拟合，获得h个支路的拟合函数Mh(t)”；
第二拟合结果输出子模块，用于若所述第三判断结果表示是，则利用公式
Figure BDA0002813879000000082
和τr＝RrCr，计算h个支路的拟合函数的每个等效支路的支路等效电阻值和支路等效电容值；
其中，Uo表示测量极化电流时所加的阶跃电压，Rr和Cr分别表示第r个等效支路的等效电阻和等效电容。
第二等效支路建立子模块，用于将τr从小到大排列，依次从中取出满足τr≥2/fl的m个支路构成第二等效支路，相应的Rr和Cr分别表示第二等效支路等效电阻和第二等效支路等效电容。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明公开了一种多层复合介质暂态过程中的电场分布确定方法，所述确定方法包括如下步骤：分别测量多层复合介质的每一层介质的介电频谱和极化电流，根据介电频谱的工频复电容实部得到每一层介质的介质等效电容，根据极化电流稳态值得到每一层介质的介质等效电阻；获得第j层介质的介质等效电阻和介质等效电容，对第j层介质的介电频谱进行拟合，获得第j层介质的n个第一等效支路的第一支路等效电阻和第一支路等效电容；所述第一等效支路为第一支路等效电阻和第一支路等效电容串联连接形成的支路；对第j层介质的极化电流进行拟合，获得第j层介质的m个第二等效支路的第二支路等效电阻和第二支路等效电容；所述第二等效支路为第二支路等效电阻和第二支路等效电容串联连接形成的支路；根据第j层介质的极化电流的时间变化曲线确定从阶跃电压施加至极化电流稳定之间的时间作为第j层介质的极化电流时间；判断第j层介质的极化电流时间是否满足公式td＜2/fl，获得第一判断结果；若所述第一判断结果表示是，则将第j层介质的介质等效电阻和介质等效电容以及n个第一等效支路并联连接，形成第j层介质的单一介质等效模型；若所述第一判断结果表示否，则将第j层介质的介质等效电阻和介质等效电容以及n个第一等效支路和m个第二等效支路并联连接，形成第j层介质的单一介质等效模型；将每一层介质的单一介质等效模型进行串联连接；将每一层介质的介质等效电阻、介质等效电容、第一支路等效电阻、第一支路等效电容、第二支路等效电阻和第二支路等效电容计算出等效介电常数和电导率，获得复合介质等效模型；对复合介质等效模型施加暂态电压，计算得到复合介质等效模型输出的电场的分布，作为多层复合介质的电场分布。本发明通过建立介质的等效模型，并结合电路分析方法即可快速算出多层介质暂态过程中的电场分布，避免了繁琐的数值方法，极大降低了计算量和计算时间。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明提供的一种多层复合介质暂态过程中的电场分布确定方法的流程图；
图2为本发明具体实施例提供的油纸绝缘结构示意图；
图3为本发明具体实施例提供的变压器油的介电频谱；
图4为本发明具体实施例提供的油浸绝缘纸的介电频谱；
图5为本发明具体实施例提供的变压器油的极化电流曲线；
图6为本发明具体实施例提供的油浸绝缘纸的极化电流曲线；
图7为本发明具体实施例提供的第一单一介质等效模型图；
图8为本发明具体实施例提供的第二单一介质等效模型图；
图9为本发明具体实施例提供的复合介质等效模型的仿真图。
具体实施方式
本发明的目的是提供一种多层复合介质暂态过程中的电场分布确定方法及系统，以提供一种快速的电气设备绝缘结构电场分布计算方法。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对发明作进一步详细的说明。
如图1所示，本发明提供一种多层复合介质暂态过程中的电场分布确定方法，所述确定方法包括如下步骤：
步骤101，分别测量多层复合介质的每一层介质的介电频谱和极化电流，根据介电频谱的工频复电容实部得到每一层介质的介质等效电容，根据极化电流稳态值得到每一层介质的介质等效电阻；所述介电频谱为：在不同频率的交流电压激励下的频域复介电常数。
步骤101所述分别测量多层复合介质的每一层介质的介电频谱和极化电流，根据介电频谱的工频复电容实部得到每一层介质的介质等效电容，根据极化电流稳态值得到每一层介质的介质等效电阻，具体包括：测量每一层介质在不同频率的交流电压激励下的频域复介电常数；测量每一层介质在阶跃直流电压激励下的极化电流。通过测量的结果就可以获得介质本身的介质等效电阻和介质等效电容，电容是50Hz下的介电常数算来的，电阻是用极化电流测量的阶跃电压比稳态电流。
由于介电频谱对于太低的电流频率下的测量结果不准确，频率为0时无法测量；同时，在极化电流的测量过程中，初始时刻的极化电流测量难度较大，采样率不高，导致在高频段的数据难以采集。介电频谱更能反映介质在高频电流下的信息，而极化电流在低频下的测量效果更好，因此采用介电频谱和极化电流相结合的测量手段。
介电频谱测量：
取厚度为d，横截面积为S的第j层介质，对两端施加频率为ω(包括从低频fl到高频fh)的交流电压，此时其内部会出现相应的极化弛豫现象，通过仪器对经过样品电流的检测，可计算得到不同频率下介质自身的频域复介电常数，进而得到复电容参数C*(包括实部和虚部)，实部表示为C'，虚部表示为C”。第j层介质在工频(50Hz)下测得的介电常数为εo，根据
Figure BDA0002813879000000111
可以得到该介质的介质等效电容C0，绘制电容参数随频率变化的曲线。
极化电流测量：
取厚度为d，横截面积为S的第j层介质，在其两端施加阶跃直流电压，用电流表测量在阶跃电压作用下的极化电流。测量开始时，先向被测试品施加电压为Uo的阶跃激励，直至电流持续稳定不变，记录整个过程流过介质的极化电流为id(t)，极化过程的时间为Tp，以及电流达到稳态值Ip。根据
Figure BDA0002813879000000112
可以得到第j层介质的介质等效电阻R0。
步骤102，对第j层介质的介电频谱进行拟合，获得第j层介质的n个第一等效支路的第一支路等效电阻和第一支路等效电容；所述第一等效支路为第一支路等效电阻和第一支路等效电容串联连接形成的支路；
步骤102所述对第j层介质的介电频谱进行拟合，获得第j层介质的n个第一等效支路的第一支路等效电阻和第一支路等效电容，具体包括：根据所述介电频谱和电容随频率变化的关系函数
Figure BDA0002813879000000113
确定不同频率下的电容的测量值；其中，C(ω)为频率ω下的复电容的测量值，ε(ω)表示频率ω下的复介电常数，S为第j层介质的横截面积，d为第j层介质的厚度；令n的数值等于1；根据不同频率下的电容的测量值，利用公式
Figure BDA0002813879000000114
对电容与介质的等效电容的差值，进行函数拟合，获得n个支路的拟合函数Nn(ω)；其中，Ci为第i个第一等效支路的等效电容值，Ri表示第i个第一等效支路的等效电阻值；选取K个频率点，利用公式
Figure BDA0002813879000000115
计算每个频率点的函数拟合误差值；其中，Nn(ωk)表示第k个频率点的n个支路的拟合函数值，C'(ωk)表示第k个频率点下的复电容的实部，Co表示第j层介质工频下的复电容实部，ek表示第k个频率点的函数拟合误差值；判断N个频率点的函数拟合误差值的平均值是否小于或等于1％，获得第二判断结果；若所述第二判断结果表示否，则令n的数值增加1，返回步骤“根据不同频率下的电容的测量值，利用公式
Figure BDA0002813879000000121
对电容与介质的等效电容的差值，进行函数拟合，获得n个支路的拟合函数Nn(ω)”；若所述第二判断结果表示是，输出n个支路的拟合函数的每个第一等效支路的等效电阻和等效电容。
步骤103，对第j层介质的极化电流进行拟合，获得第j层介质的m个第二等效支路的第二支路等效电阻和第二支路等效电容；所述第二等效支路为第二支路等效电阻和第二支路等效电容串联连接形成的支路。
步骤103所述对第j层介质的极化电流进行拟合，获得第j层介质的m个第二等效支路的第二支路等效电阻和第二支路等效电容，具体包括：令h的数值等于1；利用公式
Figure BDA0002813879000000122
对所述极化电流与极化电流最终达到的稳态值的差值，进行函数拟合，获得h个支路的拟合函数Mh(t)；其中，Ar表示h个支路的拟合函数Mh(t)的第r个等效支路的拟合函数的系数，τr表示h个支路的拟合函数Mh(t)的第r个等效支路的拟合函数的时间常数；选取L个时间点，利用公式
Figure BDA0002813879000000123
计算每个时间点的函数拟合误差值；其中，el表示第l个时间点的函数拟合误差值，Mh(tl)表示第l个时间点的h个支路的拟合函数值，id(tl)表示第l个时间点的极化电流值，IR表示极化电流最终达到的稳态值；判断L个时间点的函数拟合误差值的平均值是否小于或等于1％，获得第三判断结果；若所述第三判断结果表示否，则令h的数值增加1，返回步骤“利用公式
Figure BDA0002813879000000124
对所述极化电流与极化电流最终达到的稳态值的差值，进行函数拟合，获得h个支路的拟合函数Mm(t)”；若所述第三判断结果表示是，则利用公式
Figure BDA0002813879000000125
和τr＝RrCr，计算m个支路的拟合函数的每个等效支路的支路等效电阻值和支路等效电容值；其中，Uo表示测量极化电流时所加的阶跃电压，Rr和Cr分别表示第r个第二等效支路的等效电阻和等效电容；将τr从小到大排列，依次从中取出满足τr≥2/fl的m个支路构成第二等效支路，相应的Rr和Cr分别表示第二等效支路等效电阻和第二等效支路等效电容。
步骤104，根据第j层介质的极化电流的时间变化曲线确定从阶跃电压施加至极化电流稳定之间的时间作为第j层介质的极化电流时间；
步骤105，判断第j层介质的极化电流时间是否满足公式td＜2/fl，获得第一判断结果；其中，td表示极化电流时间，fl表示介电频谱能够测量的最低频率；
步骤106，若所述第一判断结果表示是，则将第j层介质的介质等效电阻和介质等效电容(几何电容)以及n个第一等效支路并联连接，形成第j层介质的单一介质等效模型；
步骤107，若所述第一判断结果表示否，则将第j层介质的介质等效电阻和介质等效电容(几何电容)以及n个第一等效支路和m个第二等效支路并联连接，形成第j层介质的单一介质等效模型；
步骤108，将每一层介质的单一介质等效模型进行串联连接，并将每一层介质的介质等效电阻、介质等效电容、第一支路等效电阻、第一支路等效电容、第二支路等效电阻和第二支路等效电容换算成介电常数和电导率参数，获得复合介质等效模型；
步骤109，对复合介质等效模型施加暂态电压，计算得到复合介质等效模型输出的电场的分布，作为多层复合介质的电场分布。
本发明还提供一种多层复合介质暂态过程中的电场分布确定系统，所述确定系统包括：
测量模块，用于分别测量多层复合介质的每一层介质的介电频谱和极化电流，根据介电频谱的工频复电容实部得到每一层介质的介质等效电容，根据极化电流稳态值得到每一层介质的介质等效电阻；所述介电频谱为：在不同频率的交流电压激励下的频域复介电常数。
所述测量模块，具体包括：频域复介电常数测量子模块，用于测量每一层介质在不同频率的交流电压激励下的频域复介电常数；极化电流测量子模块，用于测量每一层介质在阶跃直流电压激励下的极化电流。
介电频谱拟合模块，用于对第j层介质的介电频谱进行拟合，获得第j层介质的n个第一等效支路的第一支路等效电阻和第一支路等效电容；所述第一等效支路为第一支路等效电阻和第一支路等效电容串联连接形成的支路。
所述介电频谱拟合模块，具体包括：电容测量值确定子模块，用于根据所述介电频谱和电容随频率变化的关系函数
Figure BDA0002813879000000141
确定不同频率下的电容的测量值；其中，C(ω)为频率ω下的复电容的测量值，ε(ω)表示频率ω下的复介电常数，S为第j层介质的横截面积，d为第j层介质的厚度；第一初始化子模块，用于令n的数值等于1；函数拟合子模块，用于根据不同频率下的电容的测量值，利用公式
Figure BDA0002813879000000142
对电容与介质的等效电容的差值，进行函数拟合，获得n个支路的拟合函数Nn(ω)；其中，Ci为第i个第一等效支路的等效电容值，Ri表示第i个第一等效支路的等效电阻值；函数拟合误差值计算子模块，用于选取K个频率点，利用公式
Figure BDA0002813879000000143
计算每个频率点的函数拟合误差值；其中，Nn(ωk)表示第k个频率点的n个支路的拟合函数值，C'(ωk)表示第k个频率点下的复电容的实部，Co表示第j层介质工频下的复电容实部，ek表示第k个频率点的函数拟合误差值；第二判断子模块，用于判断N个频率点的函数拟合误差值的平均值是否小于或等于1％，获得第二判断结果；第一返回子模块，用于若所述第二判断结果表示否，则令n的数值增加1，返回步骤“根据不同频率下的电容的测量值，利用公式
Figure BDA0002813879000000144
对电容与介质的等效电容的差值，进行函数拟合，获得n个支路的拟合函数Nn(ω)”；第一拟合结果输出子模块，用于若所述第二判断结果表示是，输出n个支路的拟合函数的每个第一等效支路的等效电阻和等效电容。
极化电流拟合模块，用于对第j层介质的极化电流进行拟合，获得第j层介质的m个第二等效支路的第二支路等效电阻和第二支路等效电容；所述第二等效支路为第二支路等效电阻和第二支路等效电容串联连接形成的支路。
所述极化电流拟合模块，具体包括：第二初始化子模块，用于令h的数值等于1；第二函数拟合子模块，用于利用公式
Figure BDA0002813879000000151
对所述极化电流与极化电流最终达到的稳态值的差值，进行函数拟合，获得h个支路的拟合函数Mh(t)；其中，Ar表示h个支路的拟合函数Mh(t)的第r个等效支路的拟合函数的系数，τr表示h个支路的拟合函数Mh(t)的第r个等效支路的拟合函数的时间常数；第二函数拟合误差计算子模块，用于选取L个时间点，利用公式
Figure BDA0002813879000000152
计算每个时间点的函数拟合误差值；其中，el表示第l个时间点的函数拟合误差值，Mh(tl)表示第l个时间点的h个支路的拟合函数值，id(tl)表示第l个时间点的极化电流值，IR表示极化电流最终达到的稳态值；第三判断子模块，用于判断L个时间点的函数拟合误差值的平均值是否小于或等于1％，获得第三判断结果；第二返回子模块，用于若所述第三判断结果表示否，则令h的数值增加1，返回步骤“利用公式
Figure BDA0002813879000000153
对所述极化电流与极化电流最终达到的稳态值的差值，进行函数拟合，获得h个支路的拟合函数Mh(t)”；第二拟合结果输出子模块，用于若所述第三判断结果表示是，则利用公式
Figure BDA0002813879000000154
和τr＝RrCr，计算m个支路的拟合函数的每个等效支路的支路等效电阻值和支路等效电容值；其中，Uo表示测量极化电流时所加的阶跃电压，Rr和Cr分别表示第r个等效支路的等效电阻和等效电容。第二等效支路建立子模块，用于将τr从小到大排列，依次从中取出满足τr≥2/fl的m个支路构成第二等效支路，相应的Rr和Cr分别表示第二等效支路等效电阻和第二等效支路等效电容。
极化电流时间确定模块，用于根据第j层介质的极化电流的时间变化曲线确定从阶跃电压施加至极化电流稳定之间的时间作为第j层介质的极化电流时间；
第一判断模块，用于判断第j层介质的极化电流时间是否满足公式td＜2/fl，获得第一判断结果；其中，td表示极化电流时间，fl表示介电频谱能够测量的最低频率；
第一单一介质等效模型建立模块，用于若所述第一判断结果表示是，则将第j层介质的介质等效电阻和介质等效电容以及n个第一等效支路并联连接，形成第j层介质的单一介质等效模型；
第二单一介质等效模型建立模块，用于若所述第一判断结果表示否，则将第j层介质的介质等效电阻和介质等效电容以及n个第一等效支路和m个第二等效支路并联连接，形成第j层介质的单一介质等效模型；
复合介质等效模型建立模块，用于将每一层介质的单一介质等效模型进行串联连接，并将每一层介质的介质等效电阻、介质等效电容、第一支路等效电阻、第一支路等效电容、第二支路等效电阻和第二支路等效电容换算成介电常数和电导率参数，获得复合介质等效模型；
电场分布计算模块，用于对复合介质等效模型施加暂态电压，计算得到复合介质等效模型输出的电场的分布，作为多层复合介质的电场分布。
为了说明本发明的确定方法及系统的实现方式和技术效果，本发明还提供了一个具体的实施例。
如图2所示，以油纸绝缘结构，一层油一层纸的情况为例，求油和油浸纸的暂态过程中的电场分布。
具体步骤如下：
步骤1：取厚度分别为do、dp,横截面积分别为So，Sp的变压器油和油浸纸，分别测量变压器油和油浸纸从低频到高频下的介电频谱和极化电流：
由于介电频谱对于太低的电流频率下的测量结果不准确，频率为0时无法测量；同时，在极化电流的测量过程中，初始时刻的极化电流测量难度较大，采样率不高，导致在高频段的数据难以采集。介电频谱更能反映介质在高频电流下的信息，而极化电流在低频下的测量效果更好，因此采用介电频谱和极化电流相结合的测量手段。
介电频谱测量：
分别对两种介质两端施加频率为ω(包括从低频fl到高频fh)的交流电压，此时其内部会出现相应的极化弛豫现象，通过仪器对经过样品电流的检测，可计算得到不同频率下介质自身的频域复介电常数，进而得到复电容参数C*(包括实部和虚部)，实部表示为C'，虚部表示为C”。变压器油在工频(50Hz)下测得的电容参数实部(变压器油的介质等效电容)为Co，油浸纸在工频(50Hz)下测得的电容参数实部(油浸纸的介质等效电容)为Cp，绘制电容参数随频率变化的曲线，如图3和图4。
极化电流测量：
分别在厚度分别为do、dp，横截面积分别为So，Sp的变压器油和油浸纸两种介质两端施加阶跃直流电压，用电流表测量在阶跃电压作用下的极化电流。测量开始时，先向被测试品施加电压为Uo的阶跃激励，直至电流持续稳定不变，记录整个过程流过介质的极化电流为id(t)，如图5和图6，极化过程的时间为Tp，以及电流达到稳态值Ip。根据
Figure BDA0002813879000000174
可以得到变压器油和油浸纸的介质等效电阻Ro、Rp。
步骤2：对步骤1中测量的介电频谱、极化电流进行拟合，得到拟合参数：
(1)变压器油和绝缘纸介电频谱拟合
拟合公式：
Figure BDA0002813879000000171
Figure BDA0002813879000000172
C(ω)为介质在不同频率下的复电容，ε(ω)为不同频率下的复介电常数，ω为频率，S为被测试样的横截面积，d为被测试样的厚度，C’(ω)为不同频率下的复电容的实部，Co为介质等效电容。
拟合方法：
测量介电频谱的过程中得到介电常数随频率的变化关系ε(ω)，通过
Figure BDA0002813879000000173
得到电容随频率变化的关系C(ω)，对电容的实部进行函数拟合，设置最大允许误差1％，令n＝1，对C’(ω)-Co值进行函数拟合，此时记做
Figure BDA0002813879000000181
N1(ω)中选取若干个点，若
Figure BDA0002813879000000182
的平均误差小于等于1％，则结束，否则令n＝2，继续循环，直至误差小于1％，结束循环。得到表示介质松弛极化过程的n个不同时间常数下的电阻参数R1、R2、…、Rn和n个电容参数C1、C2、…、Cn。
油和纸的拟合方法相同。
(2)变压器油和绝缘纸极化电流拟合
拟合公式：
Figure BDA0002813879000000183
Figure BDA0002813879000000184
τi＝RiCi
id(t)为极化电流随时间的变化情况，IR为极化电流最终达到的稳态值，t为时间，Ai和τi为拟合过程中的未知参数，Uo为测量极化电流时所加的阶跃电压，Ri、Ci为表示松弛极化过程的不同时间常数下的电容电阻参数。
拟合过程：
测量极化电流的过程中得到极化电流随时间的变化关系id(t)，设置最大允许误差1％，令n＝1，对id(t)-IR的值进行指数函数拟合，得到A1、τ1，此时记做
Figure BDA0002813879000000185
M1(t)中选取若干个点，若
Figure BDA0002813879000000186
的平均误差小于等于1％，则结束，否则令n＝2，继续循环，直至误差小于1％，结束循环。得到A1、A2、…An和τ1、τ2、…τn，通过
Figure BDA0002813879000000187
和τi＝RiCi得到表示介质松弛极化过程的h个不同时间常数下的电阻参数R1’、R2’、…、Rh’和n个电容参数C1’、C2’、…、Ch’。
油和纸的拟合方法相同。
步骤3：根据步骤2极化电流和介电频谱拟合的电容电阻参数建立单一介质的等值参数模型。
判断极化时间td和介电频谱测量的最低频fl的关系，若
Figure BDA0002813879000000188
那么等效模型由介质本身的等效电阻和等效电容以及介电频谱拟合得到的不同时间常数下的等效电容电阻确定，如图7所示。
若
Figure BDA0002813879000000191
则模型由介质本身的等效电阻和等效电容、介电频谱拟合以及极化电流拟合得到的m个时间常数大于等于2/fl下的等效电容电阻建立，如图8所示，其中，R1、R2和Rn为第一支路等效电阻，C1、C2和Cn为第一支路等效电容；R1’、R2’和Rm’为第二支路等效电阻，C1’、C2’和Cm’为第二支路等效电容。
步骤4：由步骤3得到的单一介质模型串联连接得到复合介质等效模型。
计算电容电阻参数的电阻率和介电常数：
将步骤3得到的每一层介质等效模型串联连接，再将电容电阻参数通过公式
Figure BDA0002813879000000192
和
Figure BDA0002813879000000193
计算出等效电导率和介电常数，S表示被测试样的横截面积，d表示被测试样的厚度。
油和纸的等效介电常数和电导率参数如表1和表2所示，其中，σo0、εo0分别为油的等效电导率、等效介电常数，σo1、σo2为油的支路1、支路2的等效电导率，εo1、εo2为支路1、支路2的等效介电常数；σp0、εp0分别为油浸纸的等效电导率、等效介电常数，σp1、σp2、σp3为油的支路1、支路2、支路3的等效电导率，εp1、εp2、εp3为支路1、支路2、支路3的等效介电常数。
表1变压器油的等效电导率和等效介电常数参数表
Figure BDA0002813879000000194
表2油浸纸的等效电导率和等效介电常数参数表
Figure BDA0002813879000000195
Figure BDA0002813879000000201
因为介质内部任意一点电场和电流密度的关系，再结合电流连续性方程和外施电压，即可求出多层介质的电场分布。由于实施例中油和油浸纸采用平行结构，且在复合介质两端施加电压，根据实际介质的横截面积和厚度，通过
Figure BDA0002813879000000202
和
Figure BDA0002813879000000203
求得待求电场的介质的等效电容、等效电阻和不同支路的介质等效电容、介质等效电阻，如图9所示，构成复合介质的等值电路模型。其中，Ro、Co分别为油的等效电阻、等效电容，Ro1、Ro2为油的支路1、支路2的等效电阻，Co1、Co2为油的支路1、支路2的等效电容；Rp、Cp分别为油浸纸的等效电阻、等效电容，Rp1、Rp2、Rp3为油浸纸的支路1、支路2、支路3的等效电阻，Cp1、Cp2、Cp3为油浸纸的支路1、支路2、支路3的等效电容；V1为外施电压。
步骤4输出结果：复合介质的等效模型。应用在步骤5使用电路分析软件求暂态过程的电场分布。
步骤5：对步骤4得到的复合介质的等效模型求解。
将以如图9所示的电路在电路分析软件中建模，运用电路分析软件求解复合介质在暂态过程中的电场分布。
以mulitisim软件为例
1.新建工程文件
2.根据步骤4得到的电路模型以及电容电阻参数，在软件中加入相应的元器件以及设置相应的电容电阻值，如图9所示。
3.电路两端根据实际需要施加电压，这里以2kV的直流电压为例，给电路加上开关，并添加虚拟示波器元件，分别在表示绝缘油和油浸纸的模型两端连接上示波器。
4.点击运行选项，合上开关，根据电流连续可计算出油和纸上的电压，进而根据电场等于电压除以厚度得到电场分布。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明公开了一种多层复合介质暂态过程中的电场分布确定方法，所述确定方法包括如下步骤：分别测量多层复合介质的每一层介质的介电频谱和极化电流，根据介电频谱的工频复电容实部得到每一层介质的介质等效电容，根据极化电流稳态值得到每一层介质的介质等效电阻；获得第j层介质的介质等效电阻和介质等效电容；对第j层介质的介电频谱进行拟合，获得第j层介质的n个第一等效支路的第一支路等效电阻和第一支路等效电容；所述第一等效支路为第一支路等效电阻和第一支路等效电容串联连接形成的支路；对第j层介质的极化电流进行拟合，获得第j层介质的m个第二等效支路的第二支路等效电阻和第二支路等效电容；所述第二等效支路为第二支路等效电阻和第二支路等效电容串联连接形成的支路；根据第j层介质的极化电流的时间变化曲线确定从阶跃电压施加至极化电流稳定之间的时间作为第j层介质的极化电流时间；判断第j层介质的极化电流时间是否满足公式td＜2/fl，获得第一判断结果；若所述第一判断结果表示是，则将第j层介质的介质等效电阻和介质等效电容以及n个第一等效支路并联连接，形成第j层介质的单一介质等效模型；若所述第一判断结果表示否，则将第j层介质的介质等效电阻和介质等效电容以及n个第一等效支路和m个第二等效支路并联连接，形成第j层介质的单一介质等效模型；将每一层介质的单一介质等效模型进行串联连接，获得复合介质等效模型；对复合介质等效模型施加暂态电压，计算得到复合介质等效模型输出的电场的分布，作为多层复合介质的电场分布。本发明通过建立介质的等效模型，并结合电路分析方法即可快速算出多层介质暂态过程中的电场分布，避免了繁琐的数值方法，极大降低了计算量和计算时间。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例，基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

Claims (8)

1.一种多层复合介质暂态过程中的电场分布确定方法，其特征在于，所述确定方法包括如下步骤：
分别测量多层复合介质的每一层介质的介电频谱和极化电流，根据介电频谱的工频复电容实部得到每一层介质的介质等效电容，根据极化电流稳态值得到每一层介质的介质等效电阻；所述介电频谱为在不同频率的交流电压激励下的频域复介电常数；
对第j层介质的介电频谱进行拟合，获得第j层介质的n个第一等效支路的第一支路等效电阻和第一支路等效电容；所述第一等效支路为第一支路等效电阻和第一支路等效电容串联连接形成的支路；
对第j层介质的极化电流进行拟合，获得第j层介质的m个第二等效支路的第二支路等效电阻和第二支路等效电容；所述第二等效支路为第二支路等效电阻和第二支路等效电容串联连接形成的支路；
根据第j层介质的极化电流的时间变化曲线确定从阶跃电压施加至极化电流稳定之间的时间作为第j层介质的极化电流时间；
判断第j层介质的极化电流时间是否满足公式td＜2/fl，获得第一判断结果；其中，td表示极化电流时间，fl表示介电频谱能够测量的最低频率；
若所述第一判断结果表示是，则将第j层介质的介质等效电阻和介质等效电容以及n个第一等效支路并联连接，形成第j层介质的单一介质等效模型；
若所述第一判断结果表示否，则将第j层介质的介质等效电阻和介质等效电容以及n个第一等效支路和m个第二等效支路并联连接，形成第j层介质的单一介质等效模型；
将每一层介质的单一介质等效模型进行串联连接，并将每一层介质的介质等效电阻、介质等效电容、第一支路等效电阻、第一支路等效电容、第二支路等效电阻和第二支路等效电容换算成介电常数和电导率参数，获得复合介质等效模型；
对复合介质等效模型施加暂态电压，计算得到复合介质等效模型输出的电场的分布，作为多层复合介质的电场分布。
2.根据权利要求1所述的多层复合介质暂态过程中的电场分布确定方法，其特征在于，所述分别测量多层复合介质的每一层介质的介电频谱和极化电流，根据介电频谱的工频复电容实部得到每一层介质的介质等效电容，根据极化电流稳态值得到每一层介质的介质等效电阻，具体包括：
测量每一层介质在不同频率的交流电压激励下的频域复介电常数；
测量每一层介质在阶跃直流电压激励下的极化电流。
3.根据权利要求1所述的多层复合介质暂态过程中的电场分布确定方法，其特征在于，所述对第j层介质的介电频谱进行拟合，获得第j层介质的n个第一等效支路的第一支路等效电阻和第一支路等效电容，具体包括：
根据所述介电频谱和电容随频率变化的关系函数
Figure FDA0002813878990000021
确定不同频率下的电容的测量值；
其中，C(ω)为频率ω下的复电容的测量值，ε(ω)表示频率ω下的复介电常数，S为第j层介质的横截面积，d为第j层介质的厚度；
令n的数值等于1；
根据不同频率下的电容的测量值，利用公式
Figure FDA0002813878990000022
对电容与介质的等效电容的差值，进行函数拟合，获得n个支路的拟合函数Nn(ω)；
其中，Ci为第i个第一等效支路的等效电容值，Ri表示第i个第一等效支路的等效电阻值；
选取K个频率点，利用公式
Figure FDA0002813878990000023
计算每个频率点的函数拟合误差值；
其中，Nn(ωk)表示第k个频率点的n个支路的拟合函数值，C'(ωk)表示第k个频率点下的复电容的实部，Co表示第j层介质工频下的复电容实部，ek表示第k个频率点的函数拟合误差值；
判断N个频率点的函数拟合误差值的平均值是否小于或等于1％，获得第二判断结果；
若所述第二判断结果表示否，则令n的数值增加1，返回步骤“根据不同频率下的电容的测量值，利用公式
Figure FDA0002813878990000031
对电容与介质的等效电容的差值，进行函数拟合，获得n个支路的拟合函数Nn(ω)”；
若所述第二判断结果表示是，则输出n个支路的拟合函数的每个第一等效支路的等效电阻和等效电容。
4.根据权利要求1所述的多层复合介质暂态过程中的电场分布确定方法，其特征在于，所述对第j层介质的极化电流进行拟合，获得第j层介质的m个第二等效支路的第二支路等效电阻和第二支路等效电容，具体包括：
令h的数值等于1；
利用公式
Figure FDA0002813878990000032
对所述极化电流与极化电流最终达到的稳态值的差值，进行函数拟合，获得h个支路的拟合函数Mh(t)；
其中，Ar表示h个支路的拟合函数Mh(t)的第r个支路的拟合函数的系数，τr表示h个支路的拟合函数Mh(t)的第r个支路的拟合函数的时间常数；
选取L个时间点，利用公式
Figure FDA0002813878990000033
计算每个时间点的函数拟合误差值；
其中，el表示第l个时间点的函数拟合误差值，Mh(tl)表示第l个时间点的h个支路的拟合函数值，id(tl)表示第l个时间点的极化电流值，IR表示极化电流最终达到的稳态值；
判断L个时间点的函数拟合误差值的平均值是否小于或等于1％，获得第三判断结果；
若所述第三判断结果表示否，则令h的数值增加1，返回步骤“利用公式
Figure FDA0002813878990000034
对所述极化电流与极化电流最终达到的稳态值的差值，进行函数拟合，获得h个支路的拟合函数Mh(t)”；
若所述第三判断结果表示是，则利用公式
Figure FDA0002813878990000041
和τr＝RrCr，计算h个支路的拟合函数的每个等效支路的支路等效电阻值和支路等效电容值；
其中，Uo表示测量极化电流时所加的阶跃电压幅值，Rr和Cr分别表示第r个等效支路的等效电阻和等效电容；
将τr从小到大排列，依次从中取出满足τr≥2/fl的m个支路构成第二等效支路，相应的Rr和Cr分别表示第二支路等效电阻和第二支路等效电容。
5.一种多层复合介质暂态过程中的电场分布确定系统，其特征在于，所述确定系统包括：
测量模块，用于分别测量多层复合介质的每一层介质的介电频谱和极化电流，根据介电频谱的工频复电容实部得到每一层介质的介质等效电容，根据极化电流稳态值得到每一层介质的介质等效电阻；所述介电频谱为：在不同频率的交流电压激励下的频域复介电常数；
介电频谱拟合模块，用于对第j层介质的介电频谱进行拟合，获得第j层介质的n个第一等效支路的第一支路等效电阻和第一支路等效电容；所述第一等效支路为第一支路等效电阻和第一支路等效电容串联连接形成的支路；
极化电流拟合模块，用于对第j层介质的极化电流进行拟合，获得第j层介质的m个第二等效支路的第二支路等效电阻和第二支路等效电容；所述第二等效支路为第二支路等效电阻和第二支路等效电容串联连接形成的支路；
极化电流时间确定模块，用于根据第j层介质的极化电流的时间变化曲线确定从阶跃电压施加至极化电流稳定之间的时间作为第j层介质的极化电流时间；
第一判断模块，用于判断第j层介质的极化电流时间是否满足公式td＜2/fl，获得第一判断结果；其中，td表示极化电流时间，fl表示介电频谱能够测量的最低频率；
第一单一介质等效模型建立模块，用于若所述第一判断结果表示是，则将第j层介质的介质等效电阻和介质等效电容以及n个第一等效支路并联连接，形成第j层介质的单一介质等效模型；
第二单一介质等效模型建立模块，用于若所述第一判断结果表示否，则将第j层介质的介质等效电阻和介质等效电容以及n个第一等效支路和m个第二等效支路并联连接，形成第j层介质的单一介质等效模型；
复合介质等效模型建立模块，用于将每一层介质的单一介质等效模型进行串联连接，并将每一层介质的介质等效电阻、介质等效电容、第一支路等效电阻、第一支路等效电容、第二支路等效电阻和第二支路等效电容换算成介电常数和电导率参数，获得复合介质等效模型；
电场分布计算模块，用于对复合介质等效模型施加暂态电压，计算得到复合介质等效模型输出的电场的分布，作为多层复合介质的电场分布。
6.根据权利要求5所述的多层复合介质暂态过程中的电场分布确定系统，其特征在于，所述测量模块，具体包括：
频域复介电常数测量子模块，用于测量每一层介质在不同频率的交流电压激励下的频域复介电常数；
极化电流测量子模块，用于测量每一层介质在阶跃直流电压激励下的极化电流。
7.根据权利要求5所述的多层复合介质暂态过程中的电场分布确定系统，其特征在于，所述介电频谱拟合模块，具体包括：
电容测量值确定子模块，用于根据所述介电频谱和电容随频率变化的关系函数
Figure FDA0002813878990000051
确定不同频率下的电容的测量值；
其中，C(ω)为频率ω下的复电容的测量值，ε(ω)表示频率ω下的复介电常数，S为第j层介质的横截面积，d为第j层介质的厚度；
第一初始化子模块，用于令n的数值等于1；
函数拟合子模块，用于根据不同频率下的电容的测量值，利用公式
Figure FDA0002813878990000052
对电容与介质的等效电容的差值，进行函数拟合，获得n个支路的拟合函数Nn(ω)；
其中，Ci为第i个第一等效支路的等效电容值，Ri表示第i个第一等效支路的等效电阻值；
函数拟合误差值计算子模块，用于选取K个频率点，利用公式
Figure FDA0002813878990000061
计算每个频率点的函数拟合误差值；
其中，Nn(ωk)表示第k个频率点的n个支路的拟合函数值，C'(ωk)表示第k个频率点下的复电容的实部，Co表示第j层介质工频下的复电容实部，ek表示第k个频率点的函数拟合误差值；
第二判断子模块，用于判断N个频率点的函数拟合误差值的平均值是否小于或等于1％，获得第二判断结果；
第一返回子模块，用于若所述第二判断结果表示否，则令n的数值增加1，返回步骤“根据不同频率下的电容的测量值，利用公式
Figure FDA0002813878990000062
对电容与介质的等效电容的差值，进行函数拟合，获得n个支路的拟合函数Nn(ω)”；
第一拟合结果输出子模块，用于若所述第二判断结果表示是，输出n个支路的拟合函数的每个第一等效支路的等效电阻和等效电容。
8.根据权利要求5所述的多层复合介质暂态过程中的电场分布确定系统，其特征在于，所述极化电流拟合模块，具体包括：
第二初始化子模块，用于令h的数值等于1；
第二函数拟合子模块，用于利用公式
Figure FDA0002813878990000063
对所述极化电流与极化电流最终达到的稳态值的差值，进行函数拟合，获得h个支路的拟合函数Mh(t)；
其中，Ar表示h个支路的拟合函数Mh(t)的第r个等效支路的拟合函数的系数，τr表示h个支路的拟合函数Mh(t)的第r个等效支路的拟合函数的时间常数；
第二函数拟合误差计算子模块，用于选取L个时间点，利用公式
Figure FDA0002813878990000071
计算每个时间点的函数拟合误差值；
其中，el表示第l个时间点的函数拟合误差值，Mh(tl)表示第l个时间点的h个支路的拟合函数值，id(tl)表示第l个时间点的极化电流值，IR表示极化电流最终达到的稳态值；
第三判断子模块，用于判断L个时间点的函数拟合误差值的平均值是否小于或等于1％，获得第三判断结果；
第二返回子模块，用于若所述第三判断结果表示否，则令h的数值增加1，返回步骤“利用公式
Figure FDA0002813878990000072
对所述极化电流与极化电流最终达到的稳态值的差值，进行函数拟合，获得h个支路的拟合函数Mh(t)”；
第二拟合结果输出子模块，用于若所述第三判断结果表示是，则利用公式
Figure FDA0002813878990000073
和τr＝RrCr，计算h个支路的拟合函数的每个等效支路的支路等效电阻值和支路等效电容值；
其中，Uo表示测量极化电流时所加的阶跃电压，Rr和Cr分别表示第r个等效支路的等效电阻和等效电容；
第二等效支路建立子模块，用于将τr从小到大排列，依次从中取出满足τr≥2/fl的m个支路构成第二等效支路，相应的Rr和Cr分别表示第二等效支路等效电阻和第二等效支路等效电容。
CN202011405249.2A 2020-12-03 2020-12-03 一种多层复合介质暂态过程中的电场分布确定方法及系统 Active CN112698115B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011405249.2A CN112698115B (zh) 2020-12-03 2020-12-03 一种多层复合介质暂态过程中的电场分布确定方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011405249.2A CN112698115B (zh) 2020-12-03 2020-12-03 一种多层复合介质暂态过程中的电场分布确定方法及系统

Publications (2)

Publication Number Publication Date
CN112698115A true CN112698115A (zh) 2021-04-23
CN112698115B CN112698115B (zh) 2021-11-23

Family

ID=75506800

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011405249.2A Active CN112698115B (zh) 2020-12-03 2020-12-03 一种多层复合介质暂态过程中的电场分布确定方法及系统

Country Status (1)

Country Link
CN (1) CN112698115B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0271167A (ja) * 1988-09-06 1990-03-09 Victor Co Of Japan Ltd 表面電位分布検出装置
CN103091611A (zh) * 2013-01-14 2013-05-08 江苏省电力公司电力科学研究院 油纸绝缘电力设备绝缘老化状态检测方法
CN107860980A (zh) * 2017-10-17 2018-03-30 武汉泽电新材料有限公司 一种时频域结合快速介电响应测试方法
CN207662967U (zh) * 2017-11-03 2018-07-27 国家电网公司 频域介电响应测试的温度校正设备及测试设备
CN109408885A (zh) * 2018-09-19 2019-03-01 上海电力学院 一种高压直流下绝缘子空间电荷密度模型优化方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0271167A (ja) * 1988-09-06 1990-03-09 Victor Co Of Japan Ltd 表面電位分布検出装置
CN103091611A (zh) * 2013-01-14 2013-05-08 江苏省电力公司电力科学研究院 油纸绝缘电力设备绝缘老化状态检测方法
CN107860980A (zh) * 2017-10-17 2018-03-30 武汉泽电新材料有限公司 一种时频域结合快速介电响应测试方法
CN207662967U (zh) * 2017-11-03 2018-07-27 国家电网公司 频域介电响应测试的温度校正设备及测试设备
CN109408885A (zh) * 2018-09-19 2019-03-01 上海电力学院 一种高压直流下绝缘子空间电荷密度模型优化方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
应鹏宇: "《TiO2纳米粒子浓度对油纸复合绝缘击穿特性和》", 《中国电机工程学报》 *

Also Published As

Publication number Publication date
CN112698115B (zh) 2021-11-23

Similar Documents

Publication Publication Date Title
Gafvert et al. Dielectric spectroscopy in time and frequency domain applied to diagnostics of power transformers
Faranda et al. A new simplified model of double-layer capacitors
CN107609251B (zh) 一种基于fds的油纸电容式套管串联极化模型建模方法
Zhao et al. High frequency electric circuit modeling for transformer frequency response analysis studies
CN112051310B (zh) 一种油纸绝缘套管中x蜡检测及含量评估方法
Cui et al. Understanding moisture dynamics and its effect on the dielectric response of transformer insulation
Sonerud et al. Dielectric response measurements utilizing semi-square voltage waveforms
Smith et al. A dielectric frequency response model to evaluate the moisture content within an oil impregnated paper condenser bushing
Guerrero et al. Transformer oil diagnosis based on a capacitive sensor frequency response analysis
Arikan et al. Prediction of dielectric parameters of an aged MV cable: A comparison of curve fitting, decision tree and artificial neural network methods
CN112698115B (zh) 一种多层复合介质暂态过程中的电场分布确定方法及系统
Wang et al. Hybrid methodology based on extension theory for partial discharge fault diagnosis of power capacitors
Abu-Siada High frequency transformer modelling using state space representation for FRA studies
Kostiukov Measurement of dissipation factor of inner layers of insulation in three-core belted cables
Pradhan et al. Estimation of dielectric dissipation factor of cellulosic parts in oil-paper insulation by frequency domain spectroscopy
Linde et al. Comparison of Dielectric Loss Measuring Methods on Epoxy Samples under Harmonic Distorted Voltages
Saha et al. Effects of temperature on time-domain dielectric diagnostics of transformers
Onal A study for examining dissipation factors of various insulations and test transformers in the wide range of frequency
Saha et al. Some precautions for the field users of PDC measurement for transformer insulation condition assessment
Suo et al. Dynamic characteristics analysis on interface polarization and depolarization of nonlinear double-layered dielectrics
de Almeida Arantes et al. Relative permittivity meter using a capacitive sensor and an oscillating current source
Tao et al. Modelling and calculation for dielectric response circuit of oil-paper insulation transformers
CN111551792B (zh) 一种绝缘电介质无穷高频相对介电常数测量方法
CN111579879B (zh) 一种绝缘电介质无穷高频相对介电常数测量方法
Zhang et al. Dielectric frequency response of oil-paper insulation under uneven damp

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant