Power supply reliability analysis method and system based on operation and maintenance characteristics of distribution network

Abstract

The present invention discloses a power supply reliability analysis method and system based on the operation and maintenance characteristics of a distribution network, specifically related to the technical field of distribution network power supply analysis. It collects harmonic data of each node in the distribution network system in real time, including harmonic current, harmonic voltage, and frequency components. Combining with the distribution network topology structure, and by optimizing and solving a sparse linear equation set through the L1 norm, the main sources of harmonics are identified. By calculating the contribution degree of each nonlinear load to the harmonic distortion of the distribution network system, a list of key harmonic sources is screened out using Pareto analysis. Based on the identified key harmonic sources, a filter configuration optimization model is constructed, with the harmonic distortion level and filter installation cost as constraint functions, and the installation position and parameter configuration of the filter as decision variables. An optimization algorithm is used to solve it, and the recommended installation position and parameter configuration of the filter are output. The harmonic suppression effect under different load conditions is evaluated, effectively solving the problem that harmonic distortion reduces power supply reliability.

Description

Technical Field

[0001] The present invention relates to the technical field of power distribution network power supply analysis, and more specifically, to a power supply reliability analysis method and system based on the operation and maintenance characteristics of a power distribution network. Background Art

[0002] In modern power distribution network systems, the widespread use of non-linear loads such as industrial frequency converters, electric vehicle chargers, and household electronic devices has significantly increased the probability of harmonic distortion; harmonic distortion can cause equipment overheating, insulation aging, power system oscillations, and even failures, seriously threatening power supply reliability. For example, in a 10 kV power distribution network, the connection of an electric vehicle charger may cause the total harmonic distortion of a local node to rise from 3% to 15%, directly affecting the power supply quality of adjacent residential areas and highlighting the urgent need for harmonic governance.

[0003] Existing harmonic governance technologies have various deficiencies and are difficult to meet the requirements of dynamic operation scenarios of power distribution networks; firstly, traditional methods mostly adopt passive monitoring and post-fault repair strategies, such as locating problems through power quality monitors after a fault occurs, lacking the ability to actively identify and prevent harmonic sources; existing technologies mainly rely on static models in filter configuration optimization and fail to fully consider the dynamics of power distribution network operation characteristics. For example, filter parameters are not adjusted according to real-time load fluctuations, which easily leads to resonance risks or over-compensation phenomena.

[0004] Therefore, how to accurately identify the main sources of harmonics in a power distribution network through data analysis technology and dynamically optimize the layout position of filters to reduce the impact of harmonic distortion on power supply reliability has become an urgent technical problem to be solved. Summary of the Invention

[0005] In order to overcome the above-mentioned defects of the prior art, the present invention provides a power supply reliability analysis method and system based on the operation and maintenance characteristics of a power distribution network. By constructing a harmonic contribution matrix and adopting an L1 norm optimization method, key harmonic sources can be effectively identified, and filter configuration can be optimized in combination with a dynamic adjustment strategy, significantly improving power supply reliability to solve the problems proposed in the above background art.

[0006] To achieve the above object, the present invention provides the following technical solution: A power supply reliability analysis method based on the operation and maintenance characteristics of a power distribution network, including the following steps:

[0007] Step S11: Use smart meters and power quality monitors to collect harmonic data of each node in the power distribution network system in real time, including harmonic current, harmonic voltage, and frequency components;

[0008] Further, the collected data needs to be cleaned and standardized, and outliers and noise are removed to ensure data quality;

[0009] Step S12. Identify the main sources of harmonics: Combine the topology of the distribution network and solve the sparse linear equations by L1 norm optimization to identify the main sources of harmonics; Calculate the contribution degree of each nonlinear load to the harmonic distortion of the distribution network system, and use Pareto analysis to screen out the list of key harmonic sources, including the node numbers and the corresponding contribution degrees.

[0010] Step S13. Optimize the filter configuration: Based on the identified key harmonic sources, construct an optimization model for filter configuration, with the harmonic distortion level and the filter installation cost as constraint functions, and the installation location and parameter configuration of the filter as decision variables, and use an optimization algorithm to solve it, and output the recommended installation location and parameter configuration of the filter.

[0011] Step S14. Verify and dynamically adjust the optimization scheme: Verify the filter configuration and evaluate the harmonic suppression effect under different load conditions; If the verification result shows that the harmonic distortion level does not meet the standard or there is over-compensation phenomenon, then adjust the parameters of the filter configuration optimization model and solve it again; Combine the real-time monitoring data and update the filter configuration optimization model regularly to achieve dynamic optimization.

[0012] Preferably, based on the located harmonic sources, construct a harmonic contribution matrix to quantify the impact of each harmonic source on the harmonic distortion of each node in the distribution network; Use i and j to represent the sequential numbers of the nodes in the distribution network system, and calculate the contribution degree of harmonic distortion through the following formula:

[0013] ; where represents the contribution degree of node j to the harmonic voltage distortion of node i, is the harmonic impedance between node i and node j, is the harmonic current of node j, is the harmonic voltage of node i;

[0014] Based on the contribution matrix, calculate the cumulative contribution degree of each harmonic source to the total harmonic distortion of the distribution network system; Use Pareto analysis (80 / 20 principle) to screen out the harmonic sources whose cumulative contribution degree to the total harmonic distortion accounts for the top 80%, and generate a list of key harmonic sources.

[0015] Preferably, in step S11, introduce an adaptive data cleaning method based on wavelet transform to perform wavelet decomposition on the harmonic data and extract high-frequency and low-frequency components;

[0016] Use information entropy analysis to distinguish normal mutations and noise;

[0017] Filter out the noise components, retain the mutation characteristics, and reconstruct the cleaned harmonic data set;

[0018] The cleaned harmonic data set is directly input into step S12 for harmonic source location.

[0019] Preferably, it includes an anomaly detection step based on the time series characteristics of harmonic data, including:

[0020] Obtain the time series data of the harmonic of the distribution network to be detected;

[0021] Identify the abnormal features in the harmonic time series data and generate a time series anomaly profile corresponding to the abnormal features;

[0022] According to the time series anomaly profile, determine whether the harmonic data contains time series anomalies;

[0023] If the harmonic data contains time series anomalies, an anomaly detection result is generated based on the long short-term memory network model, and the input data of the filter configuration optimization scheme is adjusted to improve the system's response ability to dynamic anomalies.

[0024] Preferably, based on the dynamic changes of the distribution network topology, a graph neural network is used to update the harmonic contribution matrix. Based on the change of the edge weights of the dynamic collaborative network graph, the harmonic contribution matrix is updated in real time to generate a list of key harmonic sources adapted to the topological changes, enhancing the adaptability of harmonic source identification.

[0025] Preferably, when the load of the distribution network suddenly changes or the topology changes, a dynamic harmonic source identification model is started to identify the dynamic influence of the harmonic source and the collaborative effect between the harmonic source and the load;

[0026] Input the real-time harmonic data of each node in the distribution network system, including harmonic voltage, harmonic current and the distribution network topology;

[0027] Output an enhanced list of key harmonic sources, including node numbers and phase-time delay comprehensive contribution degrees, for subsequent filter configuration optimization.

[0028] Preferably, the operation process of the dynamic harmonic source identification model includes the following steps:

[0029] Step S21: Collect real-time harmonic data, extract voltage phase, current phase and time delay features, and construct a spatio-temporal feature matrix to ensure the data accuracy of dynamic collaborative analysis; by constructing a spatio-temporal feature matrix, ensure that subsequent analysis can accurately capture the dynamic propagation characteristics of harmonic sources; use and to represent the harmonic current sequences of nodes i and j. By comparing and the similarity at different time offsets, determine the time delay feature of the harmonic propagating from node i to node j;

[0030] Step S22: Analyze the spatio-temporal feature matrix, and generate a dynamic collaborative network graph containing node collaborative relationships and propagation paths through phase-time delay joint clustering, revealing the spatio-temporal propagation path of harmonic sources and providing a basis for contribution degree calculation;

[0031] The distance between node i and node j is calculated by the following distance function :

[0032]

[0033] where and are the weights of phase and time delay respectively, is the maximum value of the time delay feature after normalization; represents the voltage phase difference;

[0034] An improved K-means algorithm is adopted, the spatio-temporal feature matrix and the distance between nodes are input, the number of clusters K is set, and after iterative grouping, the nodes are divided into several collaborative groups;

[0035] According to the clustering results, the nodes are divided into several collaborative groups, which reflect the dynamic interaction mode between harmonic sources and loads. Based on the collaborative groups and the distribution network topology structure, a dynamic collaborative network diagram is constructed. The nodes in the diagram are the distribution network nodes, the edges in the diagram are the node collaborative relationships, and weights are assigned to each edge based on the distance;

[0036] Step S23: Process the dynamic collaborative network diagram, calculate the phase-time delay comprehensive contribution degree, and output an enhanced list of key harmonic sources to improve the pertinence of filter configuration;

[0037] Construct a dynamic collaborative network diagram, and simultaneously collect the topology structure data of each node and the corresponding line in the distribution network. Calculate the phase-time delay comprehensive contribution degree of each node through the following formula :

[0038] where is the set of adjacent nodes of node i, is the harmonic impedance between nodes; is the time delay attenuation factor, is the attenuation coefficient; represent the voltage phases of node i and node j respectively;

[0039] Step S24: Use Pareto analysis to screen and obtain an enhanced list of key harmonic sources.

[0040] Preferably, the acquisition method of the time delay feature is as follows:

[0041] Use the fast Fourier transform to perform frequency-domain decomposition on the harmonic voltage and current sequences of each node, extract the phase angles of the main harmonics, denoted as the voltage phase and the current phase, where i and j are node numbers;

[0042] Adopt the cross-correlation function Calculate the time-delay characteristic T of harmonic propagation between (topologically connected) adjacent nodes ij , defined as:

[0043]

[0044] where represents the continuous time series of the harmonic current value of node i at time s (such as the 5th harmonic current), represents the time-lag variable; is the function that takes the value to make reach the maximum value.

[0045] Preferably, it includes the steps of optimizing the filter configuration based on the harmonic propagation between nodes, including:

[0046] Obtain the harmonic data of the distribution network nodes to be analyzed;

[0047] Identify the harmonic propagation relationship between each node in the distribution network and generate the non-linear characteristic profile of the harmonic propagation between nodes;

[0048] According to the non-linear characteristic profile, determine whether the harmonic propagation between nodes has a significant non-linear effect;

[0049] If the harmonic propagation between nodes has a significant non-linear effect, generate a non-linear effect analysis result based on the support vector machine model, optimize the filter configuration scheme, and improve the suppression ability of the distribution network to complex harmonic propagation.

[0050] Preferably, the filter configuration optimization in step S13 includes: through the harmonic impedance frequency response analysis, adding the resonance risk constraint conditions to ensure that the filter parameters do not cause system resonance, outputting the recommended configuration that meets the constraints, and improving the safety of the filter configuration.

[0051] Preferably, use the long short-term memory network to predict the future load and harmonic trends, adjust the filter configuration based on the prediction results, generate a dynamic scheduling scheme and verify it through simulation, including:

[0052] Use the historical load and harmonic data to train the long short-term memory network to predict the load and harmonic trends in the future time (such as within 24 hours);

[0053] According to the prediction results, adjust and optimize the model parameters to generate the filter configuration adapted to the future load;

[0054] Verify the effectiveness of the predicted configuration through simulation, update the scheduling scheme regularly, and optimize the dynamic performance of the filter.

[0055] Preferably, under the condition of periodic fluctuations of the distribution network load, the following steps are executed to optimize harmonic suppression:

[0056] Step S31: Decompose the harmonic data in Step S11, extract the harmonic components on short-time, medium-time, and long-time scales using the empirical mode decomposition method to form a multi-time scale harmonic data set, enhancing the understanding of the harmonic propagation dynamics;

[0057] Step S32: Analyze the multi-time scale harmonic data set, construct a state transition model based on Markov chain to predict the harmonic propagation trend, generate a multi-scale harmonic propagation prediction map, and reveal potential distortion risks;

[0058] Step S33: Process the multi-scale harmonic propagation prediction map, adjust the constraint conditions of the optimization model by combining the key harmonic source list in Step S12 to output a filter configuration strategy on multi-time scales, improving the adaptability of the suppression effect;

[0059] Step S34: Verify the filter configuration strategy on multi-time scales, evaluate the harmonic suppression effect on each time scale using the simulation software in Step S14, generate a dynamically optimized filter scheduling scheme, and ensure the stability of the distribution network system.

[0060] Preferably, the method includes the following steps:

[0061] Step S41: Real-time collect the harmonic data of each node in the distribution network system, including harmonic current, harmonic voltage, frequency components, and load power fluctuation data;

[0062] Step S42: Construct a dynamic harmonic source identification model, including:

[0063] Space-time feature extraction: Extract the phase difference of harmonic voltage / current through the fast Fourier transform (FFT), and calculate the harmonic propagation time-delay feature between adjacent nodes using the cross-correlation function;

[0064] Dynamic collaborative network modeling: Based on the phase difference and time-delay features between nodes, construct a space-time feature matrix, and dynamically update the harmonic propagation path weights in combination with the graph neural network (GNN);

[0065] Phase-time delay comprehensive contribution calculation: Quantify the contribution of node i to the harmonic distortion of the system through a formula;

[0066] Step S43: Based on the output of the dynamic harmonic source identification model, screen the key harmonic source list, and optimize the filter configuration using the adaptive reinforcement learning algorithm to minimize the harmonic distortion rate and installation cost;

[0067] Step S44: Verify the harmonic suppression effect through the digital twin platform simulation, and dynamically update the model parameters in combination with the real-time monitoring data to form a closed loop.

[0068] To achieve the above object, the present invention provides the following technical solution: A power supply reliability analysis system based on the operation and maintenance characteristics of the distribution network, including:

[0069] A data acquisition module uses smart meters and power quality monitors to collect harmonic data of each node in the distribution network system in real time, including harmonic current, harmonic voltage, and frequency components.

[0070] A harmonic source identification module receives the real-time harmonic data set output by the data acquisition module. Combining with the distribution network topology structure, it optimally solves the sparse linear equations through the L1 norm, models the harmonic current injection as a sparse vector, and identifies the main sources of harmonics; calculates the contribution degree of each nonlinear load to the harmonic distortion of the distribution network system, and uses Pareto analysis to screen out the list of key harmonic sources, including node numbers and corresponding contribution degrees, and outputs the list of key harmonic sources.

[0071] A filter configuration optimization module receives the list of key harmonic sources output by the harmonic source identification module, constructs a filter configuration optimization model based on this, uses the harmonic distortion level and filter installation cost as constraint functions, uses the installation location and parameter configuration of the filter as decision variables, solves using an optimization algorithm, and outputs the recommended installation location and parameter configuration of the filter, and outputs data: filter configuration scheme.

[0072] A verification and adjustment module receives the filter configuration scheme output by the filter configuration optimization module, verifies it in combination with the distribution network simulation software, and evaluates the harmonic suppression effect under different load conditions; if the verification result shows that the harmonic distortion level does not meet the standard or there is over-compensation, the parameters of the filter configuration optimization model are adjusted, re-solved, and combined with real-time monitoring data to set a regular update frequency to achieve dynamic optimization.

[0073] The technical effects and advantages of the present invention:

[0074] (1) The power supply reliability analysis method based on the operation and maintenance characteristics of the distribution network provided by the present invention combines the distribution network topology structure and real-time harmonic data, uses the L1 norm optimization and graph neural network to update the harmonic contribution matrix, generates a key harmonic source list adapted to dynamic changes, and identifies the synergistic effect between the harmonic source and the load through phase-time delay joint clustering, effectively solving the problem of insufficient identification accuracy of traditional methods in scenarios of topological changes and load mutations.

[0075] (2) The power supply reliability analysis method based on the operation and maintenance characteristics of the distribution network provided by the present invention introduces a multi-objective optimization algorithm (such as non-dominated sorting genetic algorithm) and resonance risk constraint in the filter configuration optimization, and optimizes the configuration scheme through nonlinear influence analysis and time-series anomaly detection, significantly improving the harmonic suppression effect and system stability, and solving the problems of low filter configuration efficiency and potential safety hazards in the prior art.

[0076] (3)The power supply reliability analysis method based on the operation and maintenance characteristics of the distribution network provided by the present invention can adapt to the changes in the operation characteristics of the distribution network and ensure long-term power supply reliability by predicting future load trends through a long short-term memory network, analyzing and optimizing harmonic suppression at multiple time scales, and a dynamic adjustment mechanism updated regularly (such as step S14), thoroughly overcoming the limitations of traditional static analysis methods that cannot cope with complex dynamic environments. BRIEF DESCRIPTION OF THE DRAWINGS

[0077] Figure 1 It is a flowchart of the power supply reliability analysis method of the present invention.

[0078] Figure 2 It is a flowchart of the power supply reliability analysis method based on the dynamic harmonic source identification model of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0079] Hereinafter, exemplary embodiments of the present disclosure will be described in more detail with reference to the drawings. Although the exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided so that the present disclosure can be more thoroughly understood and the scope of the present disclosure can be fully conveyed to those skilled in the art.

[0080] At the same time, it should be understood that, for the sake of description, the sizes of the various parts shown in the drawings are not drawn in actual proportional relationships.

[0081] The following description of at least one exemplary embodiment is merely illustrative in nature and in no way limits the application or use of the present application.

[0082] Known technologies, methods, and devices for those of ordinary skill in the relevant art may not be discussed in detail, but where appropriate, the technologies, methods, and devices should be regarded as part of the specification.

[0083] It is explained that the distribution network topology structure described in the embodiments of the present invention refers to the connection relationship and layout mode among various electrical components (such as transformers, lines, switches, loads, etc.) in the distribution network, which is represented as a set of nodes and edges in graph theory in mathematics; each node in the distribution network refers to a specific point with electrical significance in the distribution network, which is the connection point of the line, the installation point of the device, or the access point of the load, and the node can be a power supply point, a load point, or an intermediate connection point; the distribution network system refers to an overall composed of the distribution network topology structure, nodes, electrical equipment (such as transformers, switches, filters), and a control system, which is used to distribute electric energy from the power supply to users.

[0084] Example 1, refer to Figure 1Flowchart of the power supply reliability analysis method. An embodiment of the present invention provides a power supply reliability analysis method based on the operation and maintenance characteristics of the distribution network, including the following steps:

[0085] Step S11: Use smart meters and power quality monitors to collect harmonic data of each node in the distribution network system in real time, including harmonic current, harmonic voltage, and frequency components; the collected data needs to be cleaned and standardized to remove outliers and noise to ensure data quality;

[0086] Step S12: Identify the main sources of harmonics: Combine the distribution network topology structure and solve the sparse linear equations through L1 norm optimization to identify the main sources of harmonics (model the injection of harmonic current as a sparse vector and use L1 norm regularization to solve); by calculating the contribution of each nonlinear load to the harmonic distortion of the distribution network system, use Pareto analysis to screen out a list of key harmonic sources, including node numbers and corresponding contribution degrees;

[0087] Explanation: Based on the located harmonic sources, construct a harmonic contribution matrix to quantify the impact of each harmonic source on the harmonic distortion of each node in the distribution network; use i and j to represent the sequential numbers of nodes in the distribution network system, and calculate the contribution degree of harmonic distortion through the following formula:

[0088] ; where represents the contribution degree of node j to the harmonic voltage distortion of node i, is the harmonic impedance between node i and node j, is the harmonic current of node j, is the harmonic voltage of node i;

[0089] Based on the contribution matrix, calculate the cumulative contribution degree of each harmonic source to the total harmonic distortion of the distribution network system; use Pareto analysis (80 / 20 principle) to screen out the harmonic sources whose cumulative contribution degree to the total harmonic distortion accounts for the top 80%, and generate a list of key harmonic sources;

[0090] Step S13: Optimize the filter configuration: Based on the key harmonic sources identified in step S12, construct an optimization model for filter configuration, use the harmonic distortion level and filter installation cost as constraint functions, and use the installation location and parameter configuration of the filter (such as capacitance, inductance value) as decision variables, and use an optimization algorithm (such as genetic algorithm, particle swarm optimization) to solve, and output the recommended installation location and parameter configuration of the filter;

[0091] Use the non-dominated sorting genetic algorithm for multi-objective optimization to generate a set of non-dominated solution sets that simultaneously optimize the total harmonic distortion level and installation cost;

[0092] Step S14: Verification and dynamic adjustment of the optimization scheme: Verify the filter configuration obtained in step S13 and evaluate the harmonic suppression effect under different load conditions; if the verification result shows that the harmonic distortion level does not meet the standard or there is overcompensation, adjust the parameters of the filter configuration optimization model and solve it again; combine real-time monitoring data to update the filter configuration optimization model regularly to achieve dynamic optimization;

[0093] For example, use a distribution network simulation software for simulation verification. For example, use the MATLAB / Simulink distribution network simulation software to verify the filter configuration obtained in step S13. If the total harmonic distortion (THD) under any load condition exceeds 5% or overcompensation is detected, trigger the subsequent adjustment mechanism; combine real-time monitoring data to set the frequency of automatically updating the filter configuration optimization model to adapt to the changes in the operating characteristics of the distribution network and achieve dynamic optimization.

[0094] Foreground summary: Harmonic data in the distribution network may generate instantaneous outliers due to load mutations (such as rapid connection or disconnection of electric vehicle charging piles) or external interferences (such as lightning strikes). Traditional cleaning methods (such as mean filtering) are difficult to distinguish normal mutations from noise, resulting in data distortion or omission of key features; existing technologies mostly focus on the elimination of outliers in the long-term trend and ignore the impact of instantaneous harmonic mutations on data quality, and such mutations may be important signals for the dynamic changes of harmonic sources. Based on this:

[0095] In a possible embodiment, in step S11, an adaptive data cleaning method based on wavelet transform is introduced to perform wavelet decomposition on the harmonic data and extract high-frequency and low-frequency components;

[0096] Use information entropy analysis to distinguish normal mutations (high information entropy) from noise (low information entropy);

[0097] Filter out the noise components, retain the mutation features, and reconstruct the cleaned harmonic data set;

[0098] The cleaned harmonic data set is directly input into step S12 for harmonic source location.

[0099] In the embodiments of the present invention, it needs to be further explained that the method includes an outlier detection step based on the time series characteristics of harmonic data, including:

[0100] Obtain the time series data of the distribution network harmonics to be detected;

[0101] Identify the outlier features in the harmonic time series data and generate a time series outlier profile corresponding to the outlier features;

[0102] Determine whether the harmonic data contains time series outliers according to the time series outlier profile;

[0103] If the harmonic data contains temporal anomalies, an anomaly detection result is generated based on the long short-term memory network model, the input data of the filter configuration optimization scheme is adjusted, and the system's response ability to dynamic anomalies is enhanced;

[0104] The anomaly detection result includes the time point of anomaly occurrence, the duration, the degree of anomaly (such as the exceeded value of the total harmonic distortion level THD), and the possible cause (for example, the abnormal injection of a specific harmonic source).

[0105] Foreground summary: In step S12 (harmonic source identification), the distribution network topology may change dynamically due to switch operations or fault isolation, resulting in distorted analysis results of the contribution degree of harmonic sources; traditional methods assume a fixed topology, do not consider the impact of dynamic changes on harmonic source location, and ignore the common dynamic adjustments (such as load transfer) in the operation and maintenance of the distribution network system, which are common in actual distribution networks; based on this, in step S12, a graph neural network is introduced to track harmonic sources under dynamic topologies.

[0106] In the embodiments of the present invention, it needs to be further explained that based on the dynamic changes of the distribution network topology, the harmonic contribution matrix is updated using a graph neural network. Based on the change of the edge weights of the dynamic collaborative network graph, the harmonic contribution matrix is updated in real time to generate a list of key harmonic sources adapted to topological changes, enhancing the adaptability of harmonic source identification.

[0107] In a possible embodiment, the rapid access of nonlinear loads (such as electric vehicle chargers, industrial frequency converters) in the distribution network system and the dynamic adjustment of the distribution network topology (such as isolating a faulty line by switch operation) will cause complex changes in the spatio-temporal distribution of harmonic sources; traditional methods mostly identify harmonic sources based on static features (such as amplitude) or a single dimension (such as spatial topology), and it is difficult to adapt to the dual impacts brought by load mutations and topological dynamics; by integrating phase features and time-delay features, a spatio-temporal collaborative analysis framework is constructed to accurately identify key harmonic sources and optimize filter configuration, improving the power supply reliability of the distribution network;

[0108] When there are load mutations (such as a sharp increase in power demand during the peak period of an electric vehicle charging station) or topological changes (such as a switch cutting off a certain line after a fault) in the distribution network, a dynamic harmonic source identification model is started to identify the dynamic impact of harmonic sources and the collaborative effect between harmonic sources and loads;

[0109] Input the real-time harmonic data of each node in the distribution network system (collected by the smart meters and power quality monitors in step S11), including harmonic voltage, harmonic current, and the distribution network topology;

[0110] Output an enhanced list of key harmonic sources, including node numbers and the comprehensive contribution degree of phase-time delay, for subsequent filter configuration optimization.

[0111] In the embodiments of the present invention, it should be further explained that the operation process of the dynamic harmonic source identification model includes the following steps:

[0112] Step S21: Collect real-time harmonic data, extract voltage phase, current phase and time-delay characteristics, construct a spatio-temporal feature matrix to ensure the data accuracy of dynamic collaborative analysis; by constructing the spatio-temporal feature matrix, ensure that the subsequent analysis can accurately capture the dynamic propagation characteristics of the harmonic source; use and to represent the harmonic current sequences of nodes i and j, and determine the time-delay characteristics of the harmonic propagating from node i to node j by comparing the and similarity at different time offsets;

[0113] Step S22: Analyze the spatio-temporal feature matrix, generate a dynamic collaborative network diagram containing node collaborative relationships and propagation paths through phase-time-delay joint clustering, reveal the spatio-temporal propagation path of the harmonic source, and provide a basis for contribution degree calculation;

[0114] Calculate the distance between node i and node j through the following distance function :

[0115]

[0116] where , are the weights of phase and time-delay respectively (the weights are adjusted according to the load mutation degree), is the maximum value of the normalized time-delay characteristics; represents the voltage phase difference;

[0117] Adopt an improved K-means algorithm, input the spatio-temporal feature matrix and the distance between nodes, set the number of clusters K, and after iterative grouping, divide the nodes into several collaborative groups;

[0118] According to the clustering results, divide the nodes into several collaborative groups, which reflect the dynamic interaction mode between the harmonic source and the load. Based on the collaborative groups and the distribution network topology, construct a dynamic collaborative network diagram. The nodes in the diagram are the distribution network nodes, and the edges in the diagram are the node collaborative relationships. Assign weight values to each edge based on the distance;

[0119] Step S23: Process the dynamic collaborative network diagram, calculate the phase-time-delay comprehensive contribution degree, and output an enhanced list of key harmonic sources to improve the pertinence of filter configuration;

[0120] Explanation: Processing the dynamic collaborative network diagram includes calculating the comprehensive contribution degree between each node based on the latest phase difference and time-delay characteristics, and determining the weight of each connection in the network accordingly; pruning the connections with weights lower than the preset threshold to remove weak associations and obtain a sparsified network structure; normalizing the weights of the remaining connections to eliminate the interference of magnitude differences on the analysis results; applying graph filtering or spectral clustering methods to the normalized network diagram to evaluate the collaborative contribution degree of each node, so as to identify the most influential nodes in the network; outputting the processed dynamic collaborative network diagram and the results of the collaborative contribution degree of each node, providing a basis for key harmonic source location and subsequent optimization;

[0121] Construct a dynamic collaborative network diagram, and simultaneously collect the topological structure data of each node and the corresponding lines in the distribution network, and calculate the phase-time-delay comprehensive contribution degree of each node through the following formula : ;

[0122] where is the set of adjacent nodes of node i, is the harmonic impedance between nodes; is the time-delay attenuation factor, is the attenuation coefficient, and its value range is from 0.1 to 0.5, which is used to adjust the weight of the influence of time delay on harmonic propagation; respectively represent the voltage phases of node i and node j;

[0123] Step S24: Use Pareto analysis to screen and obtain an enhanced list of key harmonic sources.

[0124] In a possible embodiment, under the condition of a sharp increase in non-linear loads in the distribution network, the following steps are executed to identify the collaborative influence between harmonic sources and loads:

[0125] Collect the harmonic data in step S11, extract the harmonic voltage phase and current phase characteristics of each node to form a phase characteristic data set, and improve the characterization accuracy of the dynamic relationship between harmonic sources and loads;

[0126] Analyze the phase characteristic data set, apply the phase difference clustering algorithm to identify the collaborative action mode between harmonic sources and load nodes, generate a source-load collaborative relationship diagram, and reveal the hidden harmonic propagation path;

[0127] Process the source-load collaborative relationship diagram, calculate the phase collaborative contribution degree in combination with the distribution network topology structure, output an enhanced list of key harmonic sources, and optimize the pertinence of subsequent filter configuration.

[0128] In a possible embodiment, an integrated learning method based on multi-source heterogeneous data is used to solve the problem that the identification of harmonic sources in a distribution network is affected by environmental and operation and maintenance characteristics; for example, based on harmonic data, real-time meteorological data (such as temperature, humidity), and load fluctuation data of each node in the distribution network; an integrated learning model is constructed to fuse multi-source data through feature weighting, output an enhanced harmonic contribution matrix, and generate a more robust list of key harmonic sources.

[0129] In a possible embodiment, referring to Figure 2 the flowchart of the power supply reliability analysis method based on the dynamic harmonic source identification model, the method includes the following steps:

[0130] Step S41: Real-time collect harmonic data of each node in the distribution network system, including harmonic current, harmonic voltage, frequency components, and load power fluctuation data;

[0131] Step S42: Construct a dynamic harmonic source identification model, including:

[0132] Space-time feature extraction: Extract the phase difference of harmonic voltage / current through fast Fourier transform, and calculate the harmonic propagation time-delay feature between adjacent nodes using the cross-correlation function;

[0133] Dynamic collaborative network modeling: Based on the phase difference and time-delay features between nodes, construct a space-time feature matrix, and dynamically update the harmonic propagation path weights in combination with the graph neural network;

[0134] Further, the space-time feature matrix S is denoted as , where Q represents the total number of node indices, represents the harmonic phase difference feature between node i and node j, represents the harmonic propagation time-delay feature between node i and node j.

[0135] Phase-time-delay comprehensive contribution calculation: Quantify the contribution of node i to the harmonic distortion of the system through a formula;

[0136] Step S43: Based on the output of the dynamic harmonic source identification model, screen the list of key harmonic sources, and use an adaptive reinforcement learning algorithm to optimize the filter configuration to minimize the harmonic distortion rate and installation cost;

[0137] Step S44: Verify the harmonic suppression effect through the digital twin platform, and dynamically update the model parameters in combination with real-time monitoring data to form a closed loop.

[0138] Embodiment 2. The difference between the embodiment of the present invention and Embodiment 1 is that the method further includes a filter configuration optimization step based on the harmonic propagation between nodes, including:

[0139] Obtain the harmonic data of the distribution network nodes to be analyzed;

[0140] Identify the harmonic propagation relationship between nodes in the distribution network and generate a non-linear feature profile of harmonic propagation between nodes;

[0141] Based on the non-linear feature profile, determine whether the harmonic propagation between nodes has a significant non-linear impact;

[0142] If the harmonic propagation between nodes has a significant non-linear impact, generate a non-linear impact analysis result based on the support vector machine model, optimize the filter configuration scheme, and enhance the suppression ability of the distribution network against complex harmonic propagation.

[0143] In a possible embodiment, the optimization of the filter configuration in step S13 includes: through harmonic impedance frequency response analysis, adding resonance risk constraint conditions to ensure that the filter parameters do not cause system resonance, outputting a recommended configuration that meets the constraints, and improving the safety of the filter configuration.

[0144] Explanation: After completing the harmonic impedance frequency response analysis, for the frequency range that may cause system resonance, set resonance risk constraint conditions; the constraint conditions at least include: setting an upper limit on the impedance amplitude at each frequency point to avoid generating over-amplified harmonics; monitoring the impedance phase mutation points to identify possible resonance phenomena; imposing restrictions on the frequency segments exceeding the threshold to ensure that the resonance risk is controllable during system operation.

[0145] In a possible embodiment, use the long short-term memory network to predict the future load and harmonic trends, adjust the filter configuration based on the prediction results, generate a dynamic scheduling scheme and verify it through simulation, including:

[0146] Use historical load and harmonic data to train the long short-term memory network to predict the load and harmonic trends in the future time (such as within 24 hours);

[0147] According to the prediction results, adjust and optimize the model parameters to generate a filter configuration suitable for the future load;

[0148] Verify the effectiveness of the predicted configuration through simulation, update the scheduling scheme regularly, and optimize the dynamic performance of the filter.

[0149] In a possible embodiment, under the condition of periodic fluctuations in the distribution network load, perform the following steps to optimize harmonic suppression:

[0150] Step S31: Decompose the harmonic data in step S11, adopt the empirical mode decomposition method to extract the harmonic components on short-term, medium-term, and long-term scales, form a multi-time scale harmonic data set, and enhance the understanding of the harmonic propagation dynamics;

[0151] Step S32: Analyze the multi-time scale harmonic data set, construct a state transition model based on Markov chain to predict the harmonic propagation trend, generate a multi-scale harmonic propagation prediction map, and reveal potential distortion risks;

[0152] Step S33: Process the multi-scale harmonic propagation prediction map, adjust the filter configuration to optimize the constraint conditions of the model in combination with the key harmonic source list in Step S12, and output a filter configuration strategy for multiple time scales to improve the adaptability of the suppression effect;

[0153] Explanation: The filter configuration optimization model includes:

[0154] Define the installation positions and corresponding parameters (such as capacitance, inductance, etc.) of each filter as decision variables;

[0155] Construct an objective function by weighted combination of the total harmonic distortion level and the total filter installation cost, where the harmonic distortion level and the installation cost are respectively assigned different weight coefficients to balance power quality and investment expenditure;

[0156] Set two types of constraints on the objective function: one is to ensure that the total harmonic distortion level does not exceed the preset maximum allowable value; the other is to ensure that the total filter installation cost does not exceed the budget upper limit;

[0157] Set boundary constraints on the minimum and maximum values of each filter parameter to limit the range of optional device specifications.

[0158] Step S34: Verify the filter configuration strategy for multiple time scales, use the simulation software in Step S14 to evaluate the harmonic suppression effect at each time scale, and generate a dynamically optimized filter scheduling scheme to ensure the stability of the distribution network system.

[0159] The embodiment of the present invention provides a power supply reliability analysis system based on the operation and maintenance characteristics of the distribution network, including:

[0160] A data acquisition module, using smart meters and power quality monitors, to collect harmonic data of each node in the distribution network system in real time, including harmonic current, harmonic voltage and frequency components; the collected data is subjected to cleaning and standardization processing to remove outliers and noise, and output the cleaned real-time harmonic data set (including harmonic current, harmonic voltage and frequency components);

[0161] A harmonic source identification module, receiving the real-time harmonic data set output by the data acquisition module, combining the distribution network topology structure, solving the sparse linear equations through L1 norm optimization, modeling the harmonic current injection as a sparse vector and identifying the main sources of harmonics; calculating the contribution degree of each nonlinear load to the harmonic distortion of the distribution network system, and using Pareto analysis to screen out the key harmonic source list, including node numbers and corresponding contribution degrees, and output the key harmonic source list;

[0162] The filter configuration optimization module receives the list of key harmonic sources output by the harmonic source identification module, constructs a filter configuration optimization model based on this, uses the harmonic distortion level and the filter installation cost as constraint functions, uses the installation location and parameter configuration of the filter (such as capacitance and inductance values) as decision variables, and solves using an optimization algorithm (such as genetic algorithm or particle swarm optimization), and outputs the recommended installation location and parameter configuration of the filter. Output data: filter configuration scheme (recommended installation location and parameter configuration);

[0163] The verification and adjustment module receives the filter configuration scheme output by the filter configuration optimization module, verifies it in combination with a distribution network simulation software (such as MATLAB / Simulink), and evaluates its harmonic suppression effect under different load conditions; if the verification result shows that the harmonic distortion level does not meet the standard (for example, THD exceeds 5%) or there is an over-compensation phenomenon, the parameters of the filter configuration optimization model are adjusted, re-solved, and combined with real-time monitoring data, and a regular update frequency (such as every 24 hours) is set to achieve dynamic optimization.

[0164] Finally: The above are only the preferred embodiments of the present invention and are not used to limit the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A power supply reliability analysis method based on the operation and maintenance characteristics of a distribution network, characterized in that It includes the following steps: Collect harmonic data of each node in the distribution network system in real time, including harmonic current, harmonic voltage and frequency components; Identify the main sources of harmonics: combine the distribution network topology structure, and through the L1 norm optimization to solve the sparse linear equations to identify the main sources of harmonics; extract the phase difference of harmonic voltage / current through the fast Fourier transform, and calculate the harmonic propagation time-delay characteristics between adjacent nodes by using the cross-correlation function; Based on the phase difference and time-delay characteristics between nodes, construct a spatio-temporal feature matrix, and combine the graph neural network to dynamically update the harmonic propagation path weights; by calculating the contribution degree of each nonlinear load to the harmonic distortion of the distribution network system, use the Pareto analysis to screen out the list of key harmonic sources, including the node numbers and the corresponding contribution degrees; Optimize the filter configuration: based on the key harmonic sources, construct an optimization model for filter configuration, use the harmonic distortion level and the filter installation cost as constraint functions, use the installation location and parameter configuration of the filter as decision variables, and use the optimization algorithm to solve and output the recommended installation location and parameter configuration of the filter; Verify the harmonic suppression effect through the digital twin platform, and dynamically update the model parameters in combination with the real-time monitoring data to form a closed loop.

2. The power supply reliability analysis method based on the operation and maintenance characteristics of the distribution network according to claim 1, characterized in that Based on the dynamic changes of the distribution network topology structure, use the graph neural network to update the harmonic contribution matrix. Based on the change of the edge weights of the dynamic collaborative network graph, the harmonic contribution matrix is updated in real time to generate a list of key harmonic sources adapted to the topological changes, enhancing the adaptability of harmonic source identification.

3. The power supply reliability analysis method based on the operation and maintenance characteristics of the distribution network according to claim 2, characterized in that, When the load of the distribution network suddenly changes or the topology changes, start the dynamic harmonic source identification model to identify the dynamic influence of harmonic sources and the synergistic effect between harmonic sources and loads; Input the real-time harmonic data of each node in the distribution network system, including harmonic voltage, harmonic current and the distribution network topology structure; Output an enhanced list of key harmonic sources, including node numbers and phase-time-delay comprehensive contribution degrees, for subsequent filter configuration optimization.

4. The power supply reliability analysis method based on the operation and maintenance characteristics of the distribution network according to claim 3, characterized in that, The operation process of the dynamic harmonic source identification model includes the following steps: Step S21: Collect real-time harmonic data, extract voltage phase, current phase and time-delay characteristics, construct a spatio-temporal feature matrix, and determine the time-delay characteristics of harmonic propagation from node i to node j by comparing the similarity of harmonic current sequences under different time offsets; Step S22: Analyze the spatio-temporal feature matrix, and generate a dynamic collaborative network graph including node collaborative relationships and propagation paths through phase-time-delay joint clustering; Step S23: Process the dynamic collaborative network graph, calculate the phase-time-delay comprehensive contribution degree, and output an enhanced list of key harmonic sources: Construct a dynamic collaborative network diagram and simultaneously collect the topological structure data of each node and the corresponding line in the distribution network. Calculate the comprehensive contribution degree of phase-time delay for each node through the following formula : ; where is the set of adjacent nodes of node i, is the harmonic impedance between nodes; is the time-delay attenuation factor, is the attenuation coefficient; represent the voltage phases of node i and node j respectively.

5. The power supply reliability analysis method based on the operation and maintenance characteristics of the distribution network according to claim 4, wherein The acquisition method of the time-delay characteristics is: Use the fast Fourier transform to perform frequency-domain decomposition on the harmonic voltage and current sequences of each node, and extract the voltage phase and current phase of the main harmonics, where i and j are node numbers; Using the cross-correlation function Calculate the time-delay characteristic T of harmonic propagation between adjacent nodes ij , which is defined as: ; among them, represents the continuous time series of the harmonic current value of node i at time s, represents the time lag variable; is to take the function that makes reach the maximum value.

6. The power supply reliability analysis method based on the operation and maintenance characteristics of the distribution network according to claim 1, wherein It includes steps for optimizing the filter configuration based on harmonic propagation between nodes, including: Obtain the harmonic data of the distribution network nodes to be analyzed; Identify the harmonic propagation relationships between nodes in the distribution network and generate a non-linear feature profile of harmonic propagation between nodes; According to the non-linear feature profile, determine whether there is a significant non-linear influence on harmonic propagation between nodes; If the harmonic propagation between nodes has a significant nonlinear effect, generate the nonlinear effect analysis results based on the support vector machine model and optimize the filter configuration scheme.

7. The power supply reliability analysis method based on the operation and maintenance characteristics of the distribution network according to claim 6, wherein The filter configuration optimization includes: through the harmonic impedance frequency response analysis, adding resonance risk constraint conditions to ensure that the filter parameters do not cause system resonance, outputting the recommended configuration that meets the constraints, and improving the safety of the filter configuration.

8. The power supply reliability analysis method based on the operation and maintenance characteristics of the distribution network according to claim 6, wherein Use the long short-term memory network to predict the future load and harmonic trends, adjust the filter configuration based on the prediction results, generate a dynamic scheduling scheme and verify it through simulation, including: Use historical load and harmonic data to train the long short-term memory network to predict the load and harmonic trends in the future time; According to the prediction results, adjust and optimize the model parameters to generate a filter configuration suitable for the future load; Verify the effectiveness of the predicted configuration through simulation, update the scheduling scheme regularly, and optimize the dynamic performance of the filter.

9. The power supply reliability analysis method based on the operation and maintenance characteristics of the distribution network according to claim 8, wherein, Under the condition of periodic fluctuations of the distribution network load, perform the following steps to optimize harmonic suppression: Step S31: For the decomposed harmonic data, use the empirical mode decomposition method to extract the harmonic components on short, medium, and long time scales, form a multi-time scale harmonic data set, and enhance the understanding of the harmonic propagation dynamics; Step S32: Analyze the multi-time scale harmonic data set, construct a state transition model based on the Markov chain to predict the harmonic propagation trend, generate a multi-scale harmonic propagation prediction map, and reveal the potential distortion risk; Step S33: Process the multi-scale harmonic propagation prediction map, combine the list of key harmonic sources to adjust the constraint conditions of the filter configuration optimization model, and output the filter configuration strategy on multiple time scales to improve the adaptability of the suppression effect; Step S34: Verify the filter configuration strategy on multiple time scales, use simulation software to evaluate the harmonic suppression effect on each time scale, generate a dynamically optimized filter scheduling scheme, and ensure the stability of the distribution network system.

10. A power supply reliability analysis system based on the operation and maintenance characteristics of a distribution network, characterized in that, Including: Data acquisition module: Use smart meters and power quality monitors to collect the harmonic data of each node in the distribution network system in real time, including harmonic current, harmonic voltage, and frequency components; Harmonic source identification module: Receive the real-time harmonic data set output by the data acquisition module, combine the distribution network topology structure, solve the sparse linear equations through L1 norm optimization, model the harmonic current injection as a sparse vector and identify the main sources of harmonics; calculate the contribution degree of each nonlinear load to the harmonic distortion of the distribution network system, and use Pareto analysis to screen out the list of key harmonic sources, including node numbers and corresponding contribution degrees, and output the list of key harmonic sources; Filter configuration optimization module: Receive the list of key harmonic sources output by the harmonic source identification module, construct a filter configuration optimization model based on this, use the harmonic distortion level and filter installation cost as constraint functions, use the installation location and parameter configuration of the filter as decision variables, solve using an optimization algorithm, and output the recommended installation location and parameter configuration of the filter, and output data: filter configuration scheme; The verification and adjustment module receives the filter configuration scheme output by the filter configuration optimization module, verifies it in combination with the distribution network simulation software, and evaluates the harmonic suppression effect under different load conditions; if the verification result shows that the harmonic distortion level does not meet the standard or there is an over-compensation phenomenon, it adjusts the parameters of the filter configuration optimization model, re-solves, and combines the real-time monitoring data to set a regular update frequency to achieve dynamic optimization.

Images