CN117092453B - 三芯电缆的故障定位方法、装置、设备及存储介质 - Google Patents
三芯电缆的故障定位方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN117092453B CN117092453B CN202311352631.5A CN202311352631A CN117092453B CN 117092453 B CN117092453 B CN 117092453B CN 202311352631 A CN202311352631 A CN 202311352631A CN 117092453 B CN117092453 B CN 117092453B
- Authority
- CN
- China
- Prior art keywords
- fault
- core cable
- reflection coefficient
- free
- conductor
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 239000004020 conductor Substances 0.000 claims abstract description 253
- 239000011159 matrix material Substances 0.000 claims abstract description 61
- 230000008878 coupling Effects 0.000 claims abstract description 29
- 238000010168 coupling process Methods 0.000 claims abstract description 29
- 238000005859 coupling reaction Methods 0.000 claims abstract description 29
- 238000012545 processing Methods 0.000 claims abstract description 19
- 230000006870 function Effects 0.000 claims description 119
- 239000002184 metal Substances 0.000 claims description 28
- 238000004590 computer program Methods 0.000 claims description 16
- 230000009466 transformation Effects 0.000 claims description 8
- 230000004888 barrier function Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000004807 localization Effects 0.000 claims description 7
- 230000002159 abnormal effect Effects 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 11
- 238000004891 communication Methods 0.000 description 8
- 230000007547 defect Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 238000010408 sweeping Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000009413 insulation Methods 0.000 description 4
- 238000011282 treatment Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 229910052755 nonmetal Inorganic materials 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000001953 sensory effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/083—Locating faults in cables, transmission lines, or networks according to type of conductors in cables, e.g. underground
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/11—Locating faults in cables, transmission lines, or networks using pulse reflection methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Locating Faults (AREA)
Abstract
本发明公开了一种三芯电缆的故障定位方法、装置、设备及存储介质。该方法包括：基于多导体电磁耦合理论建立无故障三芯电缆的波动方程，并获取波动方程的第一通解；根据第一通解和无故障三芯电缆在末端阻抗不匹配点处的反射系数矩阵，确定波动方程的第二通解；获取待检测三芯电缆中各导体的首端故障反射系数函数；根据第二通解和扫频信号确定各导体的首端无故障反射系数函数，并根据首端故障反射系数函数和首端无故障反射系数函数确定故障信息函数；对各故障信息函数进行处理获得待检测三芯电缆中所有导体的故障定位波形；根据故障定位波形确定待检测三芯电缆的故障位置，实现了三芯电缆的故障定位。
Description
技术领域
本发明涉及故障定位技术领域，尤其涉及一种三芯电缆的故障定位方法、装置、设备及存储介质。
背景技术
由于海上地形复杂架空线路难以敷设，为了保证电能安全稳定传输，通常使用海底电缆进行输电。常用的海底电缆包括三相单芯电缆与三相三芯电缆，与三相单芯高压电缆相比，三相三芯高压电缆电能损耗更小、传输容量更大、占地面积更小，且运输与安装的成本更低，因此得到广泛应用。
三芯海底电缆受到诸多不利因素与金属护套环流的影响后，电气性能会逐渐变差，一些破坏、老化等缺陷会慢慢变成故障，大大缩减了三芯电缆的寿命。目前已有的针对电缆的故障定位方法只限于线芯导体，但是电缆并非只有线芯一根导体，高压三芯海底电缆中含有多根导体，比如管道内三根同轴电缆的线芯与金属套，以及管道中的铠装等。这些导体之间存在复杂的电磁耦合关系，仅仅基于线芯导体建立的微元电路模型无法正确定位电缆的故障。
发明内容
本发明提供了一种三芯电缆的故障定位方法、装置、设备及存储介质，以解决三芯电缆中存在的多个导体之间存在着复杂的电磁耦合关系，无法实现准确确定三芯电缆的故障位置的问题。
根据本发明的一方面，提供了一种三芯电缆的故障定位方法，包括：
基于多导体电磁耦合理论建立无故障三芯电缆的波动方程，并获取所述波动方程的第一通解；所述第一通解包括：所述无故障三芯电缆中各导体的第一电压模量和第一电流模量；
根据所述第一通解和所述无故障三芯电缆在末端阻抗不匹配点处的反射系数矩阵，确定所述波动方程在末端阻抗不匹配情况下的第二通解；所述第二通解包括：所述无故障三芯电缆中各导体的第二电压模量和第二电流模量；
向待检测三芯电缆中的各导体发送不同频率的扫频信号，并获取所述待检测三芯电缆中各导体的首端故障反射系数函数；
根据所述第二通解和所述扫频信号确定所述无故障三芯电缆中各导体的首端无故障反射系数函数，并根据所述首端故障反射系数函数和所述首端无故障反射系数函数确定故障信息函数；
对各所述故障信息函数进行处理获得所述待检测三芯电缆中所有导体的故障定位波形；根据所述故障定位波形确定所述待检测三芯电缆的故障位置。
根据本发明的另一方面，提供了一种三芯电缆的故障定位装置，包括：
第一通解确定模块，用于基于多导体电磁耦合理论建立无故障三芯电缆的波动方程，并获取所述波动方程的第一通解；所述第一通解包括：所述无故障三芯电缆中各导体的第一电压模量和第一电流模量；
第二通解确定模块，用于根据所述第一通解和所述无故障三芯电缆在末端阻抗不匹配点处的反射系数矩阵，确定所述波动方程在末端阻抗不匹配情况下的第二通解；所述第二通解包括：所述无故障三芯电缆中各导体的第二电压模量和第二电流模量；
扫频模块，用于向待检测三芯电缆中的各导体发送不同频率的扫频信号，并获取所述待检测三芯电缆中各导体的首端故障反射系数函数；
信息函数确定模块，用于根据所述第二通解和所述扫频信号确定所述无故障三芯电缆中各导体的首端无故障反射系数函数，并根据所述首端故障反射系数函数和所述首端无故障反射系数函数确定故障信息函数；
故障定位模块，用于对各所述故障信息函数进行处理获得所述待检测三芯电缆中所有导体的故障定位波形；根据所述故障定位波形确定所述待检测三芯电缆的故障位置。
根据本发明的另一方面，提供了一种电子设备，所述电子设备包括：
至少一个处理器；以及
与所述至少一个处理器通信连接的存储器；其中，
所述存储器存储有可被所述至少一个处理器执行的计算机程序，所述计算机程序被所述至少一个处理器执行，以使所述至少一个处理器能够执行本发明任一实施例所述的三芯电缆的故障定位方法。
根据本发明的另一方面，提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使处理器执行时实现本发明任一实施例所述的三芯电缆的故障定位方法。
本发明实施例的技术方案，通过基于多导体电磁耦合理论建立无故障三芯电缆的波动方程，并获取波动方程的第一通解；第一通解包括：无故障三芯电缆中各导体的第一电压模量和第一电流模量； 根据第一通解和无故障三芯电缆在末端阻抗不匹配点处的反射系数矩阵，确定波动方程在末端阻抗不匹配情况下的第二通解；第二通解包括：无故障三芯电缆中各导体的第二电压模量和第二电流模量；向待检测三芯电缆中的各导体发送不同频率的扫频信号，并获取待检测三芯电缆中各导体的首端故障反射系数函数；根据第二通解和扫频信号确定无故障三芯电缆中各导体的首端无故障反射系数函数，并根据首端故障反射系数函数和首端无故障反射系数函数确定故障信息函数；对各故障信息函数进行处理获得待检测三芯电缆中所有导体的故障定位波形；根据故障定位波形确定待检测三芯电缆的故障位置，采用频域反射技术和多导体电磁耦合理论确定故障定位波形，从而实现了准确定位三芯电缆的故障位置，解决了由于三芯电缆中存在的多个导体之间存在着复杂的电磁耦合关系，无法实现准确定位三芯电缆的故障位置的问题。
应当理解，本部分所描述的内容并非旨在标识本发明的实施例的关键或重要特征，也不用于限制本发明的范围。本发明的其它特征将通过以下的说明书而变得容易理解。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例一提供的一种三芯电缆的故障定位方法的流程图；
图2是三芯电缆的层状结构示意图；
图3是三芯电缆的接地方式的示意图；
图4是本发明实施例二提供的一种三芯电缆的故障定位方法的流程图；
图5是无故障三芯电缆的微元阻抗电路模型的示意图；
图6是无故障三芯电缆的微元导纳电路模型的示意图；
图7是本发明实施例三提供的一种三芯电缆的故障定位装置的结构示意图；
图8是实现本发明实施例的三芯电缆的故障定位方法的电子设备的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、***、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
频域反射技术（Frequency Domain Reflectometer，FDR）是依据电缆首端反射系数实现故障定位的有效方法，但是以往关于电缆FDR故障定位的方法一般是针对同轴电缆建立计算模型，而三芯电缆自身包含多根导体，多根导体之间存在复杂的电磁耦合关系。采用同轴电缆的单导体计算模型不能正确反映三芯电缆在FDR方法下的定位结果。同时，同一个的故障会引起不同导体定位波形发生不同程度的变化，而现有故障定位方法仅测量线芯导体的定位波形，这样不仅难以检测到远离线芯导体的故障，还难以定位三芯电缆的故障位置。针对上述问题，本发明提出的基于多导体电磁耦合理论的三芯电缆FDR故障定位方法，可以更加准确的给出三芯电缆多根导体在FDR方法下的定位波形，根据定位波形确定电缆的缺陷或故障发生的位置以及所在电缆的结构。
实施例一
图1是本发明实施例一提供了一种三芯电缆的故障定位方法的流程图，本实施例可适用于针对存在多导体的三芯电缆进行故障定位的情况，该方法可以由三芯电缆的故障定位装置来执行，该三芯电缆的故障定位装置可以采用硬件和/或软件的形式实现，该三芯电缆的故障定位装置可配置于电子设备中。如图1所示，该方法包括：
S110、基于多导体电磁耦合理论建立无故障三芯电缆的波动方程，并获取波动方程的第一通解；第一通解包括：无故障三芯电缆中各导体的第一电压模量和第一电流模量。
其中，三芯电缆中的ABC三相同轴电缆在管道中彼此紧密接触，因此除了同相中的线芯与金属护套之间存在电磁耦合关系，不同相之间的金属护套、金属护套与铠装之间同样存在电磁耦合关系。这些导体共同构成了三芯电缆的多导体***，它们之间存在着复杂的电磁耦合关系。示例性的，图2是三芯电缆的层状结构示意图。如图2所示，三芯电缆由外层向内层依次包括：外皮层12、铠装11、内垫层10、包带9、填充8和ABC三相同轴电缆；每一相的同轴电缆由外层向内层依次包括：非金属护套7、金属护套6、半导电阻水带5、外半导电屏蔽层4、绝缘层3、内半导电屏蔽层2和线芯1。图3是三芯电缆的接地方式的示意图，如图3所示，三芯电缆的首端开路，末端线芯正常接负载，金属护套与铠装正常接地。三芯电缆一般包含7根导体，分别是内部三根同轴电缆的线芯1与金属护套6以及最外层的铠装11。
其中，无故障三芯电缆的波动方程是针对无故障的三芯电缆的多导体微元电路模型（如各导体的阻抗、导纳、电压和电流之间的关系）所建立的能够表示无故障的三芯电缆中的各导体的电压和电流的波动情况的方程。
在本实施例中，将多导体电磁耦合理论应用于电缆首端反射系数的理论分析，依据电缆的相邻导体之间构成的回路和电磁耦合关系构建三芯电缆在无故障状态下的波动方程。并根据单位阻抗矩阵与单位导纳矩阵进一步求解无故障三芯电缆的波动方程的第一通解。第一通解包括：无故障三芯电缆中各导体在任意处的电压模量和电流模量。
S120、根据第一通解和无故障三芯电缆在末端阻抗不匹配点处的反射系数矩阵，确定波动方程在末端阻抗不匹配情况下的第二通解；第二通解包括：无故障三芯电缆中各导体的第二电压模量和第二电流模量。
其中，三芯电缆在阻抗不匹配的情况下会引起行波折反射，从而产生对应的反射系数。因此，还需要依据行波折反射原理所产生的反射系数矩阵进一步确定波动方程在末端阻抗不匹配情况下的第二通解。
在本实施例中，依据行波折反射原理能够计算无故障三芯电缆在末端阻抗不匹配点处的反射系数矩阵；根据反射系数矩阵对波动方程的第一通解进行改写，可以得到在考虑末端阻抗不匹配情况下波动方程的第二通解。
S130、向待检测三芯电缆中的各导体发送不同频率的扫频信号，并获取待检测三芯电缆中各导体的首端故障反射系数函数。
其中，待检测三芯电缆是指需要进行故障定位的电缆。首端无故障反射系数为无故障的三芯电缆在末端阻抗不匹配的情况下首端产生的行波折反射对应的反射系数。在本实施例中，通过扫频仪每次分别向无故障三芯电缆中的一个导体发送一个预设频率的扫频信号，并测量得到无故障三芯电缆中该导体的首端故障反射系数；根据不同频率的扫频信号得到的首端故障反射系数确定该导体的首端故障反射系数函数。
S140、根据第二通解和扫频信号确定待检测三芯电缆中各导体的首端无故障反射系数函数，并根据首端故障反射系数函数和首端无故障反射系数函数确定故障信息函数。
其中，首端故障反射系数是指待检测三芯电缆在末端阻抗不匹配的情况下首端产生的行波折反射对应的反射系数。故障信息函数能够反映待检测三芯电缆在末端阻抗不匹配的情况下与无故障的三芯电缆，在行波折反射上呈现的波形差距。
在本实施例中，对于待检测三芯电缆中每个导体，根据第二通解和扫频信号确定导体的首端无故障反射系数函数，并根据首端无故障反射系数函数和首端故障反射系数确定故障信息函数。
S150、对各故障信息函数进行处理获得待检测三芯电缆中所有导体的故障定位波形；根据故障定位波形确定待检测三芯电缆的故障位置。
其中，待检测三芯电缆的故障位置可以包括：定位三芯电缆在横截面上的故障（例如可以包括三芯电缆中的导体故障，也可以包括导体之间的结构层故障）或者三芯电缆在长度方向上的故障位置。
在本实施例中，通过对故障信息函数进行去噪和时域定位等处理可以获得待检测三芯电缆中各导体的故障定位波形；根据预设故障定位策略对各导体的故障定位波形进行分析，确定待检测三芯电缆的故障位置。
在一个可选的实施例中，待检测三芯电缆在长度方向上的故障定位的方法可以是：根据待检测三芯电缆中任意一个导体的故障定位波形确定。若导体的故障定位波形出现异常的尖峰，则将尖峰所在位置确定为待检测三芯电缆在长度方向上的故障位置。
在另一个可选的实施中。待检测三芯电缆在横截面上的故障定位的方法可以是：获取故障导体的数目；故障导体可以是在对应的故障定位波形中出现异常尖峰的导体；根据故障导体的数目和各故障导体的位置关系定位待检测三芯电缆的故障层。
本发明实施例的技术方案，通过基于多导体电磁耦合理论建立无故障三芯电缆的波动方程，并获取波动方程的第一通解；第一通解包括：无故障三芯电缆中各导体的第一电压模量和第一电流模量；根据第一通解和无故障三芯电缆在末端阻抗不匹配点处的反射系数矩阵，确定波动方程在末端阻抗不匹配情况下的第二通解；第二通解包括：无故障三芯电缆中各导体的第二电压模量和第二电流模量；向待检测三芯电缆中的各导体发送不同频率的扫频信号，并获取待检测三芯电缆中各导体的首端故障反射系数函数；根据第二通解和扫频信号确定无故障三芯电缆中各导体的首端无故障反射系数函数，并根据首端故障反射系数函数和首端无故障反射系数函数确定故障信息函数；对各故障信息函数进行处理获得待检测三芯电缆中所有导体的故障定位波形；根据故障定位波形确定待检测三芯电缆的故障位置。采用频域反射技术和多导体电磁耦合理论确定故障定位波形，从而实现了准确定位三芯电缆的故障位置，解决了由于三芯电缆中存在的多个导体之间存在着复杂的电磁耦合关系，无法实现准确定位三芯电缆的故障位置的问题。
实施例二
图4是本发明实施例二提供的一种三芯电缆的故障定位方法的流程图，本实施例对上述实施例的各导体的故障信息函数的确定方式进一步限定。如图4所示，该方法包括：
S210、基于多导体电磁耦合理论建立无故障三芯电缆的波动方程，并获取波动方程的第一通解；第一通解包括：无故障三芯电缆中各导体的第一电压模量和第一电流模量。
在本实施例的一个可选实施例中，基于多导体电磁耦合理论建立无故障三芯电缆的波动方程，并获取所述波动方程的第一通解，包括：
S211、根据所述无故障三芯电缆的微元电路模型，建立所述无故障三芯电缆的单位阻抗矩阵、各导体的电流与电压之间的第一微分方程，以及所述无故障三芯电缆的单位导纳矩阵、各导体的电流与电压之间的第二微分方程；
S212、基于多导体电磁耦合理论，根据所述第一微分方程和所述第二微分方程，确定所述无故障三芯电缆的波动方程；
S213、使用相模变换计算方法，计算所述波动方程的第一通解。
图5是无故障三芯电缆的微元阻抗电路模型的示意图；图6是无故障三芯电缆的微元导纳电路模型的示意图。在本实施例中，S211的具体步骤可以为：（1）将无故障三芯电缆进行微元划分，每段长度无限小，记为；（2）根据如图5所示的微元阻抗电路模型建立无故障三芯电缆的单位阻抗矩阵与各导体的电流和电压之间的第一微分方程，如
公式1：；
公式2：
公式3：；
公式4：；
公式5：
公式6：；/>；
公式7：。
其中，表示无故障三芯电缆的三相线芯的电压，c表示线芯，/>表示无故障三芯电缆的A、B和C三相；/>表示无故障三芯电缆的三相金属护套的电压，s表示金属护套；/>是三芯电缆最外层的铠装的电压；/>表示无故障三芯电缆的三相线芯的电流；/>表示无故障三芯电缆的三相金属护套的电流；/>表示无故障三芯电缆最外层的铠装的电流。单位阻抗矩阵/>为对称矩阵；/>表示无故障三芯电缆的三相线芯与金属护套的自阻抗，/>；表示障三芯电缆的铠装的自阻抗；/>表示无故障三芯电缆的同相线芯与金属护套之间的互阻抗，/>，/>且/>；/>表示无故障三芯电缆的铠装与三相线芯以及三相金属护套之间的互阻抗。
、/>、/>、/>是为了方便描述三芯电缆单位阻抗矩阵而定义的中间变量，具体计算公式如式2至公式7所示。其中，/>表示三芯电缆三相中线芯及金属护套导体内表面单位阻抗；/>表示三芯电缆的三相线芯及金属护套的外表面单位阻抗；/>表示三芯电缆三相的金属护套的内外表面互阻抗；/>、/>、/>分别表示三芯电缆最外层铠装的内、外表面单位阻抗及单位互阻抗；/>表示三芯电缆中三相线芯、金属护套和铠装间的绝缘构成的单位绝缘阻抗；/>表示三芯电缆最外层铠装与大地之间绝缘构成的单位绝缘阻抗；表示第/>相金属护套与第/>相金属护套之间的单位阻抗。/>是大地单位阻抗。
（3）根据如图6所示的微元导纳电路模型建立无故障三芯电缆的单位导纳矩阵与各导体的电流和电压之间的第二微分方程，如
公式8：；
公式9：；
公式10：；
公式11：；
公式12：
公式13：；
公式14：；
其中，单位导纳矩阵Y为称矩阵。表示无故障三芯电缆ABC三相线芯与金属护套的自导纳；/>是无故障三芯电缆的铠装的自导纳。/>表示无故障三芯电缆不同相金属护套之间的互导纳。/>表示无故障三芯电缆的铠装与ABC三相金属护套之间的互导纳。
、/>、/>、/>是为了方便描述三芯电缆的单位导纳矩阵而定义的中间变量，具体计算公式如公式9至公式14所示。其中，/>表示三芯电缆三相中线芯与金属护套之间的导纳；/>表示三芯电缆三相中金属护套与铠装之间的导纳；/>表示三芯电缆铠装与大地之间的导纳；/>表示三芯电缆中相邻相金属护套间的导纳。
S212的具体步骤可以为：将上述第一微分方程和第二微分方程中的导体电压变量矩阵记为U，导体电流矩阵记为I，单位阻抗矩阵记为Z，单位导纳矩阵记为Y；考虑了多导体间电磁耦合作用的无故障三芯电缆的微分方程可以为简记为如下形式：
公式15：
公式16：
于是，根据公式1至公式2以及公式15至公式16，可以将无故障三芯电缆的波动方程变化为：
公式17：
公式18：
S213的具体步骤可以为：相模变换需足如下关系式：
公式19：
公式20：；/>。
其中，表示导体电压相量变换矩阵；/>表示导体电流相量变换矩阵；/>表示对角矩阵，/>表示导体电压模量；/>表示导体电流模量。
根据公式17至公式20，可以得到相模变换后的无故障三芯电缆的波动方程如：
公式21：
公式22：
对相模变换后的无故障三芯电缆的波动方程进行求解，得到第一通解为：
公式23：
公式24：
其中，表示电缆的模量传输系数，/>表示模量特征阻抗，/>和分别表示无故障三芯电缆首端的前行电压波与反行电压波，它们的数值由边界条件给定。
S220、根据第一通解和无故障三芯电缆在末端阻抗不匹配点处的反射系数矩阵，确定波动方程在末端阻抗不匹配情况下的第二通解；第二通解包括：无故障三芯电缆中各导体的第二电压模量和第二电流模量。
在本发明实施例的一个可选实施例中，根据所述第一通解和所述无故障三芯电缆在末端阻抗不匹配点处的反射系数矩阵，确定所述波动方程在末端阻抗不匹配情况下的第二通解，包括：
S221、建立所述无故障三芯电缆在末端不匹配阻抗情况下的电路模型；
S222、基于行波折反射原理，根据所述电路模型确定所述无故障三芯电缆的末端反射系数矩阵和首端反射系数矩阵；
S223、将所述末端反射系数矩阵和所述首端反射系数矩阵代入所述第一通解，获得所述波动方程在末端阻抗不匹配情况下的第二通解。
在本实施例中，S221的具体步骤可以为：建立一段均匀无故障的三芯电缆与末端不匹配阻抗的电路模型，设定电缆长度为l，末端阻抗为，且满足三芯电缆的波动方程的第一通解。从而，根据行波折反射原理，可以得到在末端x=l处电缆末端多导体电压/>和电流/>，如：
公式25：；
公式26：
其中， 和/>分别表示无故障三芯电缆的末端的前行波电压与反行波电压。
S222的具体步骤可以为：
由欧姆定律易知，无故障三芯电缆的末端的电压和电流/>满足：
公式27：；其中，/>表示无故障三芯电缆的末端连接的负载的阻抗矩阵。
末端的电压中的前行波电压/>和反行波电压/>之间满足行波折反射原理，具体表现为：
28：；其中，F表示无故障三芯电缆的末端反射系数矩阵。
根据公式25至公式28，可得无故障三芯电缆的末端反射系数矩阵为：
公式29：
在首端x=0处，无故障三芯电缆的首端电压和首端电流/>满足：
公式30：
公式31：
首端电压中的前行波电压/>和反行波电压/>间满足行波折反射原理可以表示为：
公式32：其中，K是表示无故障三芯电缆的首端反射系数矩阵。
根据公式28、公式30和公式32，可以得到首端反射系数矩阵K为：
公式33：
根据公式32，可以进一步将公式30改写为：
公式34：其中，E是n阶单位矩阵，n与无故障三芯电缆中的导体个数对应。
S223的具体步骤可以为：根据无故障三芯电缆的波动方程的第一通解确定波动方程在末端阻抗不匹配情况下的第二电压模量和第二电流模量，如：
公式35：
公式36：
S230、向待检测三芯电缆中的各导体发送不同频率的扫频信号，并获取待检测三芯电缆中各导体的首端故障反射系数函数。
在本实施例中，通过扫频仪向待检测三芯电缆中的各导体发送不同频率的扫频信号，并根据获取到不同频率的扫频信号在待检测三芯电缆中的各导体形成的首端故障反射系数确定首端故障反射系数函数。
S240、计算各扫频信号在无故障三芯电缆中各导体上产生的反射信号。
在本实施例中，记扫频信号源（如扫频仪）发出的扫频信号为e，同时扫频信号源存在内阻抗。由于扫频信号源对每个导体分别发出不同频率的扫频信号，且每次仅对无故障三芯电缆中的一根导体进行扫频，因此仅在扫频信号射向的导体存在信号与内阻抗，其余参数均为0。示例性的，扫频信号源发出信号的幅值为10V，起始频率为1MHz，截止频率为100MHz，共10001个采样点。
在本实施例中，根据公式31和公式32，进一步确定无故障三芯电缆的首端电流与前行波电压/>之间的关系式，可以表示为：
公式37：；
利用首端电压与首端电流/>确定无故障三芯电缆在末端不匹配阻抗情况下的首端输入阻抗矩阵/>，可以表示为：
公式38：。
记扫频信号源发出的信号为e，同时扫频信号源存在内阻抗。由于每次仅对三芯电缆其中的一根导体打入扫频信号，因此仅在打入信号的那根导体处存在信号与内阻抗，其余参数均为0。采用扫频信号源对无故障三芯电缆进行扫频的过程中，扫频信号源内的阻抗与电缆为串联关系，满足欧姆定律，结合公式30可以推导出扫频信号e与无故障三芯电缆的电压/>之间的关系，表示为：
公式39：；
考虑了扫频信号源信号e的电缆各导体的电压和电流/>，可以表示为：
公式40：
公式41：
根据无故障三芯电缆的首端电压、首端电流/>、信号源内阻抗/>算出扫频信号源接收到的反射信号/>，可以表示为：
公式42：。
S250、根据各反射信号和第二通解确定无故障三芯电缆中各导体的首端无故障反射系数函数；首端无故障反射系数函数是首端无故障反射系数随频率变化的函数。
在本实施例中，扫频信号源接收到的发射信可分成前行波电压/>、反行波电压/>，可以表示为：
公式43：。
根据扫频信号源的内阻抗与三芯电缆的首端输入阻抗/>算出反射系数/>，可以表示为：
公式44：；
扫频信号源接收到的发射信号是前行波电压与反行波电压/>的叠加，前行波电压与反行波电压之间的关系满足行波折反射原理，具体可以表示为：
公式45：；
将公式43至公式45代入公式42，计算得到前行波与反行波/>，具体可以表示为：
公式46：；
公式47：；
由于与/>均是模量参数，需要根据公式20对其进行转换，获得实际前行波与反行波/>，具体可以表示为：
公式48：；
公式49：；
由于扫频信号每次仅扫频无故障三芯电缆中的一个导体，将被扫频的导体作为目标导体。目标导体前行波是矩阵/>中的第/>个目标导体对应的元素，目标导体反行波是矩阵/>是中的第/>个目标导体对应的元素。因此，无故障三芯电缆中的目标导体的首端反射系数/>，可以表示为：
公式50：。
选定扫频信号源的扫频信号e、起始频率与截止频率/>以及频率测试点个数N，循环执行通过扫频信号源向无故障三芯电缆中的每个导体分别发出不同频率的扫频信号，并根据扫频信号源的扫频信号计算扫频信号源接收到的反射信号；根据反射信号和第二通解确定无故障三芯电缆中的每个导体的首端反射系数，最终获得无故障三芯电缆中的每根导体在不同预设频率下的首端反射系数/>随频率变化的函数/>。
S260、对于待检测三芯电缆中的各导体，将首端故障反射系数函数和首端无故障反射系数函数的差，确定为导体的故障信息函数。
在本实施例中，将首端故障反射系数函数和首端无故障反射系数函数的差确定为故障信息函数/>，可以表示为：/>
公式51：。
S270、对于每个导体，对导体的故障信息函数进行加窗处理；对加窗后的故障信息函数进行傅里叶逆变换处理，获得故障信息时域定位函数。
在本实施例中，可以选择切比雪夫窗对原始故障信息函数实部进行处理，处理后的故障信息函数实部函数，可以表示为：
公式52：；
其中，real表示取出括号内数据的实部，chebwin表示切比雪夫窗函数，表示同规模矩阵对应元素相乘。
采用离散傅里叶逆变换处理加窗后的故障信息实部函数，故障信息时域定位函数/>，可以表示为：
公式53：
其中，N表示故障信息实部函数数组的长度，/>表示故障信息实部函数数组中数据的变量，/>是与信号起始频率/>与截止频率/>有关的时间数组，且满足，/>，/>；/>为时间分辨率，/>为扫频信号往返一周用的总时间，/>为待检测三芯电缆的总长度，/>为扫频信号在待检测三芯电缆中的波速。
S280、将故障信息时域定位函数转换到空间域，获得导体的故障定位波形。
在本实施例中，根据扫频信号的波速与时间、距离的关系将故障信息时域定位函数转换到空间域，获得待检测三芯电缆中的所有导体的定位波形。
扫频信号走过的距离、波速/>、时间/>之间的关系，可以表示为：
公式54：。
由于扫频信号源接收到的是在电缆里往返一周的信号，因此信号实际走过距离为两倍的电缆总长度。据此，计算出时间数组对应的空间数组/>，如下所示，空间数组的作用是将每个时间对应的时域定位信息与具体的位置对应起来，空间数组/>可以表示为：
公式55：；
用空间数组替换时间数组/>，最终可以获得待检测三芯电缆所有导体的故障信息空间域函数/>，该函数便是定位波形。
S290、获取故障导体的数目；根据故障导体的数目和各故障导体的位置关系定位待检测三芯电缆的故障层。
其中，故障导体是在对应的故障定位波形中出现异常尖峰的导体。
在本实施中，从分析所有导体的定位波形尖峰开始，再研究线芯导体、金属护套导体、铠装导体的定位波形之间的关系，获得待检测三芯电缆中故障的位置与所在结构，具体步骤包括：
（1）确定定位波形出现异常尖峰的故障导体的数目。待检测三芯电缆存在故障会导致导体的定位波形出现异于常态的尖峰。检测线芯导体、金属护套导体、铠装导体的定位波形尖峰值，只要某一位置的尖峰数量大于等于1，则认为该位置的待检测三芯电缆存在故障。
（2）根据故障导体之间的位置关系确定故障层。
若仅存在一个故障导体，那么该故障导体即为故障层。
若存在两个故障导体，当两个故障导体不相邻时，故障单独发生在两个故障导体上；当两根导体相邻时，故障发生在两个故障导体以及两个故障导体之间的结构层上。
若存在三个故障导体，当三个故障导体均不相邻时，故障分别发生在ABC三相的线芯上；当两个故障导体之间存在相邻关系时，故障发生在第一、二根相邻导体及其中间所夹的结构层上，且第三个故障导体上单独存在故障；当三个故障导体之间存在相邻关系时，故障发生在贯穿三个故障导体及其中间的结构层上。
若存在四个故障导体，当四个故障导体均不相邻时，故障分别发生在ABC三相的线芯以及铠装上；当两个故障导体之间存在相邻关系时，故障发生在第一、二个相邻故障导体及其中间所夹的结构层上，且第三、四个故障导体上单独存在缺陷或故障；当三个故障导体之间存在连续的相邻关系时，故障发生在第一、二、三个连续的相邻故障导体及其中间所夹的结构层上，且第四个故障导体上单独存在缺陷或故障；当四个故障导体之间存在连续的相邻关系时，故障发生在贯穿四个故障导体及其中间的结构层上；
若存在五个故障导体，由于电缆结构限制，至少会有三个故障导体之间存在连续的相邻关系，当三根导体之间存在连续的相邻关系时，故障发生在第一、二、三个连续的相邻故障导体及其中间所夹的结构上，且第四、五故障导体上单独存在缺陷或故障；当四个故障导体之间存在相邻关系时，故障发生在第一、二、三、四根连续的相邻故障导体及其中间所夹的结构层上，且第五个故障导体上单独存在缺陷或故障；当五个故障导体之间存在连续的相邻关系时，故障发生在贯穿五个故障导体及其中间的结构层上；
若存在六个故障导体，由于电缆结构限制，至少会有五个故障导体之间存在连续的相邻关系，且故障发生在第一到五根连续的相邻故障导体及其中间所夹的结构上，且第六个故障导体上单独存在故障；当六个故障导体之间存在连续的相邻关系时，故障发生在贯穿六个故障导体及其中间的结构层上；
若存在七个故障导体，由于电缆结构限制，七个导体均之间必然存在连续的相邻关系，故障发生在贯穿七个故障导体及其中间的结构层上。
实施例三
图7是本发明实施例三提供的一种三芯电缆的故障定位装置的结构示意图。如图7所示，该装置包括：第一通解确定模块310、第二通解确定模块320、扫频模块330、信息函数确定模块340和故障定位模块350；其中，
第一通解确定模块310，用于基于多导体电磁耦合理论建立无故障三芯电缆的波动方程，并获取所述波动方程的第一通解；所述第一通解包括：所述无故障三芯电缆中各导体的第一电压模量和第一电流模量；
第二通解确定模块320，用于根据所述第一通解和所述无故障三芯电缆在末端阻抗不匹配点处的反射系数矩阵，确定所述波动方程在末端阻抗不匹配情况下的第二通解；所述第二通解包括：所述无故障三芯电缆中各导体的第二电压模量和第二电流模量；
扫频模块330，用于向待检测三芯电缆中的各导体发送不同频率的扫频信号，并获取所述待检测三芯电缆中各导体的首端故障反射系数函数；
信息函数确定模块340，用于根据所述第二通解和所述扫频信号确定所述无故障三芯电缆中各导体的首端无故障反射系数函数，并根据所述首端故障反射系数函数和所述首端无故障反射系数函数确定故障信息函数；
故障定位模块350，用于对各所述故障信息函数进行处理获得所述待检测三芯电缆中所有导体的故障定位波形；根据所述故障定位波形确定所述待检测三芯电缆的故障位置。
可选的，所述第一通解确定模块，具体用于：
根据所述无故障三芯电缆的微元电路模型，建立所述无故障三芯电缆的单位阻抗矩阵、各导体的电流与电压之间的第一微分方程，以及所述无故障三芯电缆的单位导纳矩阵、各导体的电流与电压之间的第二微分方程；
基于多导体电磁耦合理论，根据所述第一微分方程和所述第二微分方程，确定所述无故障三芯电缆的波动方程；
使用相模变换计算方法，计算所述波动方程的第一通解。
可选的，所述第二通解确定模块，具体用于：
建立所述无故障三芯电缆在末端不匹配阻抗情况下的电路模型；
基于行波折反射原理，根据所述电路模型确定所述无故障三芯电缆的末端反射系数矩阵和首端反射系数矩阵；
将所述末端反射系数矩阵和所述首端反射系数矩阵代入所述第一通解，获得所述波动方程在末端阻抗不匹配情况下的第二通解。
可选的，所述信息函数确定模块，具体用于：
计算各所述扫频信号在所述无故障三芯电缆中各导体上产生的反射信号；
根据各所述反射信号和所述第二通解确定所述无故障三芯电缆中各导体的首端无故障反射系数函数；所述首端无故障反射系数函数是所述首端无故障反射系数随频率变化的函数；
对于所述待检测三芯电缆中的各导体，将所述首端故障反射系数函数和所述首端无故障反射系数函数的差确定为所述导体的故障信息函数。
可选的，所述故障定位模块，具体用于：
对于每个导体，对所述导体的故障信息函数进行加窗处理；
对加窗后的故障信息函数进行傅里叶逆变换处理，获得故障信息时域定位函数；
将所述故障信息时域定位函数转换到空间域，获得所述导体的故障定位波形。
可选的，所述故障定位模块，还用于：
获取故障导体的数目；所述故障导体是在对应的故障定位波形中出现异常尖峰的导体；
根据所述故障导体的数目和各所述故障导体的位置关系定位所述待检测三芯电缆的故障层。
本发明实施例所提供的三芯电缆的故障定位装置可执行本发明任意实施例所提供的三芯电缆的故障定位方法，具备执行方法相应的功能模块和有益效果。
实施例四
图8示出了可以用来实施本发明的实施例的电子设备10的结构示意图。电子设备旨在表示各种形式的数字计算机，诸如，膝上型计算机、台式计算机、工作台、个人数字助理、服务器、刀片式服务器、大型计算机、和其它适合的计算机。电子设备还可以表示各种形式的移动装置，诸如，个人数字处理、蜂窝电话、智能电话、可穿戴设备（如头盔、眼镜、手表等）和其它类似的计算装置。本文所示的部件、它们的连接和关系、以及它们的功能仅仅作为示例，并且不意在限制本文中描述的和/或者要求的本发明的实现。
如图8所示，电子设备10包括至少一个处理器11，以及与至少一个处理器11通信连接的存储器，如只读存储器（ROM）12、随机访问存储器（RAM）13等，其中，存储器存储有可被至少一个处理器执行的计算机程序，处理器11可以根据存储在只读存储器（ROM）12中的计算机程序或者从存储单元18加载到随机访问存储器（RAM）13中的计算机程序，来执行各种适当的动作和处理。在RAM 13中，还可存储电子设备10操作所需的各种程序和数据。处理器11、ROM 12以及RAM 13通过总线14彼此相连。输入/输出（I/O）接口15也连接至总线14。
电子设备10中的多个部件连接至I/O接口15，包括：输入单元16，例如键盘、鼠标等；输出单元17，例如各种类型的显示器、扬声器等；存储单元18，例如磁盘、光盘等；以及通信单元19，例如网卡、调制解调器、无线通信收发机等。通信单元19允许电子设备10通过诸如因特网的计算机网络和/或各种电信网络与其他设备交换信息/数据。
处理器11可以是各种具有处理和计算能力的通用和/或专用处理组件。处理器11的一些示例包括但不限于中央处理单元（CPU）、图形处理单元（GPU）、各种专用的人工智能（AI）计算芯片、各种运行机器学习模型算法的处理器、数字信号处理器（DSP）、以及任何适当的处理器、控制器、微控制器等。处理器11执行上文所描述的各个方法和处理，例如三芯电缆的故障定位方法。
在一些实施例中，三芯电缆的故障定位方法可被实现为计算机程序，其被有形地包含于计算机可读存储介质，例如存储单元18。在一些实施例中，计算机程序的部分或者全部可以经由ROM 12和/或通信单元19而被载入和/或安装到电子设备10上。当计算机程序加载到RAM 13并由处理器11执行时，可以执行上文描述的三芯电缆的故障定位方法的一个或多个步骤。备选地，在其他实施例中，处理器11可以通过其他任何适当的方式（例如，借助于固件）而被配置为执行三芯电缆的故障定位方法。
本文中以上描述的***和技术的各种实施方式可以在数字电子电路***、集成电路***、场可编程门阵列（FPGA）、专用集成电路（ASIC）、专用标准产品（ASSP）、芯片上***的***（SOC）、负载可编程逻辑设备（CPLD）、计算机硬件、固件、软件、和/或它们的组合中实现。这些各种实施方式可以包括：实施在一个或者多个计算机程序中，该一个或者多个计算机程序可在包括至少一个可编程处理器的可编程***上执行和/或解释，该可编程处理器可以是专用或者通用可编程处理器，可以从存储***、至少一个输入装置、和至少一个输出装置接收数据和指令，并且将数据和指令传输至该存储***、该至少一个输入装置、和该至少一个输出装置。
用于实施本发明的方法的计算机程序可以采用一个或多个编程语言的任何组合来编写。这些计算机程序可以提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器，使得计算机程序当由处理器执行时使流程图和/或框图中所规定的功能/操作被实施。计算机程序可以完全在机器上执行、部分地在机器上执行，作为独立软件包部分地在机器上执行且部分地在远程机器上执行或完全在远程机器或服务器上执行。
在本发明的上下文中，计算机可读存储介质可以是有形的介质，其可以包含或存储以供指令执行***、装置或设备使用或与指令执行***、装置或设备结合地使用的计算机程序。计算机可读存储介质可以包括但不限于电子的、磁性的、光学的、电磁的、红外的、或半导体***、装置或设备，或者上述内容的任何合适组合。备选地，计算机可读存储介质可以是机器可读信号介质。机器可读存储介质的更具体示例会包括基于一个或多个线的电气连接、便携式计算机盘、硬盘、随机存取存储器（RAM）、只读存储器（ROM）、可擦除可编程只读存储器（EPROM或快闪存储器）、光纤、便捷式紧凑盘只读存储器（CD-ROM）、光学储存设备、磁储存设备、或上述内容的任何合适组合。
为了提供与用户的交互，可以在电子设备上实施此处描述的***和技术，该电子设备具有：用于向用户显示信息的显示装置（例如，CRT（阴极射线管）或者LCD（液晶显示器）监视器）；以及键盘和指向装置（例如，鼠标或者轨迹球），用户可以通过该键盘和该指向装置来将输入提供给电子设备。其它种类的装置还可以用于提供与用户的交互；例如，提供给用户的反馈可以是任何形式的传感反馈（例如，视觉反馈、听觉反馈、或者触觉反馈）；并且可以用任何形式（包括声输入、语音输入或者、触觉输入）来接收来自用户的输入。
可以将此处描述的***和技术实施在包括后台部件的计算***（例如，作为数据服务器）、或者包括中间件部件的计算***（例如，应用服务器）、或者包括前端部件的计算***（例如，具有图形用户界面或者网络浏览器的用户计算机，用户可以通过该图形用户界面或者该网络浏览器来与此处描述的***和技术的实施方式交互）、或者包括这种后台部件、中间件部件、或者前端部件的任何组合的计算***中。可以通过任何形式或者介质的数字数据通信（例如，通信网络）来将***的部件相互连接。通信网络的示例包括：局域网（LAN）、广域网（WAN）、区块链网络和互联网。
计算***可以包括客户端和服务器。客户端和服务器一般远离彼此并且通常通过通信网络进行交互。通过在相应的计算机上运行并且彼此具有客户端-服务器关系的计算机程序来产生客户端和服务器的关系。服务器可以是云服务器，又称为云计算服务器或云主机，是云计算服务体系中的一项主机产品，以解决了传统物理主机与VPS服务中，存在的管理难度大，业务扩展性弱的缺陷。
应该理解，可以使用上面所示的各种形式的流程，重新排序、增加或删除步骤。例如，本发明中记载的各步骤可以并行地执行也可以顺序地执行也可以不同的次序执行，只要能够实现本发明的技术方案所期望的结果，本文在此不进行限制。
上述具体实施方式，并不构成对本发明保护范围的限制。本领域技术人员应该明白的是，根据设计要求和其他因素，可以进行各种修改、组合、子组合和替代。任何在本发明的精神和原则之内所作的修改、等同替换和改进等，均应包含在本发明保护范围之内。
Claims (9)
1.一种三芯电缆的故障定位方法，其特征在于，包括：
基于多导体电磁耦合理论建立无故障三芯电缆的波动方程，并获取所述波动方程的第一通解；所述第一通解包括：所述无故障三芯电缆中各导体的第一电压模量和第一电流模量；其中，三芯电缆的多导体包括：三根同轴电缆的线芯、金属护套以及铠装；
根据所述第一通解和所述无故障三芯电缆在末端阻抗不匹配点处的反射系数矩阵，确定所述波动方程在末端阻抗不匹配情况下的第二通解；所述第二通解包括：所述无故障三芯电缆中各导体的第二电压模量和第二电流模量；
向待检测三芯电缆中的各导体发送不同频率的扫频信号，并获取所述待检测三芯电缆中各导体的首端故障反射系数函数；
根据所述第二通解和所述扫频信号确定所述无故障三芯电缆中各导体的首端无故障反射系数函数，并根据所述首端故障反射系数函数和所述首端无故障反射系数函数确定故障信息函数；首端故障反射系数是指待检测三芯电缆在末端阻抗不匹配的情况下首端产生的行波折反射对应的反射系数；所述首端故障反射系数函数是所述首端故障反射系数随频率变化的函数；首端无故障反射系数为无故障的三芯电缆在末端阻抗不匹配的情况下首端产生的行波折反射对应的反射系数；所述首端无故障反射系数函数是所述首端无故障反射系数随频率变化的函数；故障信息函数能够反映待检测三芯电缆在末端阻抗不匹配的情况下与无故障的三芯电缆，在行波折反射上呈现的波形差距；
对各所述故障信息函数进行处理获得所述待检测三芯电缆中所有导体的故障定位波形；根据所述故障定位波形确定所述待检测三芯电缆的故障位置；
其中，所述根据所述第一通解和所述无故障三芯电缆在末端阻抗不匹配点处的反射系数矩阵，确定所述波动方程在末端阻抗不匹配情况下的第二通解，包括：
建立所述无故障三芯电缆在末端不匹配阻抗情况下的电路模型；
基于行波折反射原理，根据所述电路模型确定所述无故障三芯电缆的末端反射系数矩阵和首端反射系数矩阵；
将所述末端反射系数矩阵和所述首端反射系数矩阵代入所述第一通解，获得所述波动方程在末端阻抗不匹配情况下的第二通解。
2.根据权利要求1所述的方法，其特征在于，所述基于多导体电磁耦合理论建立无故障三芯电缆的波动方程，并获取所述波动方程的第一通解，包括：
根据所述无故障三芯电缆的微元电路模型，建立所述无故障三芯电缆的单位阻抗矩阵、各导体的电流与电压之间的第一微分方程，以及所述无故障三芯电缆的单位导纳矩阵、各导体的电流与电压之间的第二微分方程；
基于多导体电磁耦合理论，根据所述第一微分方程和所述第二微分方程，确定所述无故障三芯电缆的波动方程；
使用相模变换计算方法，计算所述波动方程的第一通解。
3.根据权利要求1所述的方法，其特征在于，所述根据所述第二通解和所述扫频信号，确定所述无故障三芯电缆中各导体的首端无故障反射系数函数，包括：
计算各所述扫频信号在所述无故障三芯电缆中各导体上产生的反射信号；
根据各所述反射信号和所述第二通解确定所述无故障三芯电缆中各导体的首端无故障反射系数函数。
4.根据权利要求3所述的方法，其特征在于，所述根据所述首端故障反射系数和所述首端无故障反射系数确定故障信息函数，包括：
对于所述待检测三芯电缆中的各导体，将所述首端故障反射系数函数和所述首端无故障反射系数函数的差确定为所述导体的故障信息函数。
5.根据权利要求1所述的方法，其特征在于，所述对各所述故障信息函数进行处理获得所述待检测三芯电缆中所有导体的故障定位波形，包括：
对于每个导体，对所述导体的故障信息函数进行加窗处理；
对加窗后的故障信息函数进行傅里叶逆变换处理，获得故障信息时域定位函数；
将所述故障信息时域定位函数转换到空间域，获得所述导体的故障定位波形。
6.根据权利要求1或5所述的方法，其特征在于，所述根据所述故障定位波形确定所述待检测三芯电缆的故障位置，包括：
获取故障导体的数目；所述故障导体是在对应的故障定位波形中出现异常尖峰的导体；
根据所述故障导体的数目和各所述故障导体的位置关系定位所述待检测三芯电缆的故障层。
7.一种三芯电缆的故障定位装置，其特征在于，包括：
第一通解确定模块，用于基于多导体电磁耦合理论建立无故障三芯电缆的波动方程，并获取所述波动方程的第一通解；所述第一通解包括：所述无故障三芯电缆中各导体的第一电压模量和第一电流模量；其中，三芯电缆的多导体包括：三根同轴电缆的线芯、金属护套以及铠装；
第二通解确定模块，用于根据所述第一通解和所述无故障三芯电缆在末端阻抗不匹配点处的反射系数矩阵，确定所述波动方程在末端阻抗不匹配情况下的第二通解；所述第二通解包括：所述无故障三芯电缆中各导体的第二电压模量和第二电流模量；
扫频模块，用于向待检测三芯电缆中的各导体发送不同频率的扫频信号，并获取所述待检测三芯电缆中各导体的首端故障反射系数函数；
信息函数确定模块，用于根据所述第二通解和所述扫频信号确定所述无故障三芯电缆中各导体的首端无故障反射系数函数，并根据所述首端故障反射系数函数和所述首端无故障反射系数函数确定故障信息函数；首端故障反射系数是指待检测三芯电缆在末端阻抗不匹配的情况下首端产生的行波折反射对应的反射系数；所述首端故障反射系数函数是所述首端故障反射系数随频率变化的函数；首端无故障反射系数为无故障的三芯电缆在末端阻抗不匹配的情况下首端产生的行波折反射对应的反射系数；所述首端无故障反射系数函数是所述首端无故障反射系数随频率变化的函数；故障信息函数能够反映待检测三芯电缆在末端阻抗不匹配的情况下与无故障的三芯电缆，在行波折反射上呈现的波形差距；
故障定位模块，用于对各所述故障信息函数进行处理获得所述待检测三芯电缆中所有导体的故障定位波形；根据所述故障定位波形确定所述待检测三芯电缆的故障位置；
其中，所述第二通解确定模块，具体用于：
建立所述无故障三芯电缆在末端不匹配阻抗情况下的电路模型；
基于行波折反射原理，根据所述电路模型确定所述无故障三芯电缆的末端反射系数矩阵和首端反射系数矩阵；
将所述末端反射系数矩阵和所述首端反射系数矩阵代入所述第一通解，获得所述波动方程在末端阻抗不匹配情况下的第二通解。
8.一种电子设备，其特征在于，所述电子设备包括：
至少一个处理器；以及
与所述至少一个处理器通信连接的存储器；其中，
所述存储器存储有可被所述至少一个处理器执行的计算机程序，所述计算机程序被所述至少一个处理器执行，以使所述至少一个处理器能够执行权利要求1-6中任一项所述的三芯电缆的故障定位方法。
9.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使处理器执行时实现权利要求1-6中任一项所述的三芯电缆的故障定位方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311352631.5A CN117092453B (zh) | 2023-10-19 | 2023-10-19 | 三芯电缆的故障定位方法、装置、设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311352631.5A CN117092453B (zh) | 2023-10-19 | 2023-10-19 | 三芯电缆的故障定位方法、装置、设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117092453A CN117092453A (zh) | 2023-11-21 |
| CN117092453B true CN117092453B (zh) | 2023-12-19 |
Family
ID=88777216
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311352631.5A Active CN117092453B (zh) | 2023-10-19 | 2023-10-19 | 三芯电缆的故障定位方法、装置、设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117092453B (zh) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108445350A (zh) * | 2018-02-23 | 2018-08-24 | 北京交通大学 | 基于输入导纳的直流输电线路故障定位方法 |
| CN111413586A (zh) * | 2020-03-26 | 2020-07-14 | 成都飞机工业（集团）有限责任公司 | 电缆检测设备 |
| CN113075501A (zh) * | 2021-03-26 | 2021-07-06 | 华中科技大学 | 一种基于阻抗谱周期特征的电缆故障定位方法及*** |
| CN113466612A (zh) * | 2021-06-02 | 2021-10-01 | 湖南省计量检测研究院 | 一种电缆故障识别方法 |
| CN113655342A (zh) * | 2021-09-23 | 2021-11-16 | 国网江苏省电力有限公司连云港供电分公司 | 一种基于多导体耦合模型的三芯电缆单相故障定位方法 |
| CN114966320A (zh) * | 2022-05-27 | 2022-08-30 | 华东交通大学 | 一种基于深度学习的单端故障行波测距新算法 |
| CN115047287A (zh) * | 2022-05-06 | 2022-09-13 | 广东核电合营有限公司 | 一种核电仪控电缆的故障缆芯定位方法和*** |
| CN115184718A (zh) * | 2022-04-20 | 2022-10-14 | 中国矿业大学 | 一种基于时域分析的三芯电缆故障定位方法 |
| CN116776626A (zh) * | 2023-07-04 | 2023-09-19 | 西安理工大学 | 地-各向异性电离层波导中甚低频电波场强计算方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7512503B2 (en) * | 2003-05-12 | 2009-03-31 | Simmonds Precision Products, Inc. | Wire fault detection |
| NO341197B1 (no) * | 2012-10-24 | 2017-09-11 | Wirescan As | Fremgangsmåte og system for overvåkning av tilstanden til elektriske kabler |
| FR3032804B1 (fr) * | 2015-02-16 | 2017-02-10 | Commissariat Energie Atomique | Procede de caracterisation d'un defaut non franc dans un cable |
| FR3050036B1 (fr) * | 2016-04-08 | 2020-06-19 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Procede de detection de defauts non francs dans un cable par fusion de donnees |
-
2023
- 2023-10-19 CN CN202311352631.5A patent/CN117092453B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108445350A (zh) * | 2018-02-23 | 2018-08-24 | 北京交通大学 | 基于输入导纳的直流输电线路故障定位方法 |
| CN111413586A (zh) * | 2020-03-26 | 2020-07-14 | 成都飞机工业（集团）有限责任公司 | 电缆检测设备 |
| CN113075501A (zh) * | 2021-03-26 | 2021-07-06 | 华中科技大学 | 一种基于阻抗谱周期特征的电缆故障定位方法及*** |
| CN113466612A (zh) * | 2021-06-02 | 2021-10-01 | 湖南省计量检测研究院 | 一种电缆故障识别方法 |
| CN113655342A (zh) * | 2021-09-23 | 2021-11-16 | 国网江苏省电力有限公司连云港供电分公司 | 一种基于多导体耦合模型的三芯电缆单相故障定位方法 |
| CN115184718A (zh) * | 2022-04-20 | 2022-10-14 | 中国矿业大学 | 一种基于时域分析的三芯电缆故障定位方法 |
| CN115047287A (zh) * | 2022-05-06 | 2022-09-13 | 广东核电合营有限公司 | 一种核电仪控电缆的故障缆芯定位方法和*** |
| CN114966320A (zh) * | 2022-05-27 | 2022-08-30 | 华东交通大学 | 一种基于深度学习的单端故障行波测距新算法 |
| CN116776626A (zh) * | 2023-07-04 | 2023-09-19 | 西安理工大学 | 地-各向异性电离层波导中甚低频电波场强计算方法 |
Non-Patent Citations (3)
| Title |
|---|
| Online localization method for local defects in power cables based on electromagnetic coupling;H.jiaxing et al.;2022 18th ACDC;第1-6页 * |
| 基于频域反射法的电缆故障定位测量***研究;张巍;中国优秀硕士学位论文全文数据库工程科技Ⅱ辑(第10期);第23-55页 * |
| 基于频域反射系数谱的电缆故障定位与故障类型识别方法研究;梁钟颖 等;电工电能新技术;第41卷(第8期);第79-88页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117092453A (zh) | 2023-11-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Smaïl et al. | Detection and location of defects in wiring networks using time-domain reflectometry and neural networks | |
| Bang et al. | Classification of faults in multicore cable via time–frequency domain reflectometry | |
| CN111679159B (zh) | 一种频域反射法中阻抗变化类型判断的方法 | |
| Zhang et al. | Analysis of time-domain reflectometry combined with wavelet transform for fault detection in aircraft shielded cables | |
| NO341197B1 (no) | Fremgangsmåte og system for overvåkning av tilstanden til elektriske kabler | |
| CN113253046B (zh) | 一种基于阻抗谱技术的电缆水树枝故障定位方法 | |
| CN110146784B (zh) | 一种基于阻抗相位变化比率的电缆局部缺陷定位方法 | |
| US11061060B2 (en) | Method and system for locating defects on an electric cable | |
| CN113702754A (zh) | 采用加窗傅里叶变换的配电电缆缺陷的定位算法 | |
| CN110261739B (zh) | 一种电缆软故障定位装置及定位方法 | |
| CN114217166A (zh) | 一种基于fdr频域波形的变电站低压电缆局部缺陷定位方法 | |
| CN115389877A (zh) | 电缆绝缘故障的定位方法、装置、终端及存储介质 | |
| Tang et al. | A frequency sweep location method for soft faults of power cables based on MUSIC-pseudospectrum | |
| CN117092453B (zh) | 三芯电缆的故障定位方法、装置、设备及存储介质 | |
| CN116859182A (zh) | 考虑电缆接头的频域反射谱高压电缆缺陷定位方法及*** | |
| CN115128403B (zh) | 一种用于电力电缆的在线信号测量方法 | |
| CN116540016A (zh) | 基于反射系数离散序列的电缆缺陷诊断方法 | |
| CN114675128A (zh) | 基于护层电流和电压的海底电缆绝缘故障在线定位方法 | |
| Shi et al. | Automated wire fault location using impedance spectroscopy and differential evolution | |
| CN110426598B (zh) | 一种通信电缆屏蔽层故障定位的方法及*** | |
| Minullin et al. | Model-Experimental Detection of Faults on Overhead Power Transmission Lines by the Location Method | |
| Hassen et al. | Shielding damage characterization in twisted pair cables using OMTDR-based reflectometry and inverse problems | |
| Kingston et al. | Spread Spectrum Time Domain Reflectometry (SSTDR) digital twin simulation of photovoltaic systems for fault detection and location | |
| Zhang et al. | Research on the Influence Mechanism of Thermal Aging Degree of Cable on Frequency Domain Reflectometry and Defect Location Accuracy | |
| Shirkoohi | Modelling and simulation of fault detection in shielded twisted pair cables |
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 |