Method and device for positioning electric leakage fault of branch line in front of meter, storage medium and product

Abstract

The invention discloses a pre-meter branch line electric leakage fault positioning method and device, a storage medium and a product. The positioning method comprises the steps of obtaining a transformer area total meter voltage, a user electric meter voltage, a user electric meter active power and a user electric meter reactive power in a measurement period; the virtual impedance of the user is solved; extracting global correlation features and adjacent correlation features of the user virtual impedance; calculating the symmetric relative entropy of the global correlation characteristics and the adjacent correlation characteristics according to the extracted global correlation characteristics and the extracted adjacent correlation characteristics; calculating a comprehensive abnormal score of the user virtual impedance according to the user virtual impedance and the symmetric relative entropy; and carrying out electric leakage fault positioning according to the comprehensive abnormal score of the user virtual impedance. The electric leakage fault positioning efficiency is improved, and the operation and maintenance work efficiency of the low-voltage power distribution network is effectively improved.

Description

Technical Field

[0001] The present invention belongs to the technical field of abnormal positioning of low-voltage distribution networks, and in particular relates to a method, device, storage medium and product for positioning leakage faults in branch lines before user meters in a substation. Background Art

[0002] Due to factors such as insulation moisture or aging, partial insulation damage in overhead lines or underground pipelines in the power supply line of the substation (referring to the power supply range or power supply area of ​​a transformer) leads to leakage, and the leakage location is highly concealed.

[0003] At present, leakage problems caused by insulation damage mainly rely on manual experience and traditional instrument measurement methods. That is, operation and maintenance personnel usually climb the pole from the center point of the main line with the help of tools such as multimeters and clamp ammeters, and narrow the scope of leakage fault investigation by measuring the residual current section by section. This method has limitations in real-time, accuracy and intermittent leakage detection and positioning capabilities. Summary of the invention

[0004] The purpose of the present invention is to provide a method, device, storage medium and product for locating leakage faults of branch lines before a meter, so as to solve the problems of poor real-time performance, poor accuracy and low efficiency of traditional leakage locating methods.

[0005] The present invention solves the above technical problems through the following technical solutions: a method for locating leakage faults of branch lines before a meter, comprising:

[0006] Obtain the total meter voltage of the area, the user meter voltage, the user meter active power and the user meter reactive power within the measurement period;

[0007] Based on the total meter voltage, user meter voltage, user meter active power and user meter reactive power, solve the user virtual impedance;

[0008] Extracting global correlation features and neighboring correlation features from the user virtual impedance;

[0009] The symmetric relative entropy of the global correlation feature and the neighboring correlation feature is calculated according to the extracted global correlation feature and the neighboring correlation feature;

[0010] Calculate a comprehensive abnormality score of the user virtual impedance according to the user virtual impedance and the symmetric relative entropy;

[0011] The leakage fault is located according to the comprehensive abnormal score of the user virtual impedance.

[0012] Furthermore, based on the total meter voltage, the user meter voltage, the user meter active power and the user meter reactive power, the user virtual impedance is solved, including:

[0013] Based on the user meter voltage, user meter active power and user meter reactive power, the total meter voltage estimation function is constructed. The specific expression is:

[0014] ;

[0015] in, It represents the estimated value of the total meter voltage corresponding to the nth user at the tth time, represents the active power of the nth user's meter at the tth moment, represents the reactive power of the nth user's meter at the tth moment, represents the nth user virtual resistance at the tth moment, represents the virtual reactance of the nth user at the tth moment, Indicates the voltage of the nth user's electric meter at the tth moment;

[0016] According to the total meter voltage estimation function and the total meter voltage, the user virtual impedance is solved; wherein the solution function is:

[0017] ;

[0018] ;

[0019] Where T represents the number of sampling moments in the measurement period, N represents the number of user meters in the area, represents the total meter voltage at the tth moment, It represents the virtual impedance of the nth user at the tth moment.

[0020] Further, a pre-trained neural network model is used to extract global correlation features and neighboring correlation features of the user virtual impedance;

[0021] The neural network model includes a feature extraction module, and the feature extraction module is used to extract global correlation features and neighboring correlation features of the user virtual impedance. The specific process includes:

[0022] According to the user virtual impedance, the query vector, key vector, value vector and Gaussian kernel scale parameter are calculated. The specific calculation formula is:

[0023] , , , ;

[0024] in, , , and They represent the query vector, key vector, value vector and Gaussian kernel scale parameter of the lth layer in the feature extraction module respectively; Represents the input amount of the lth layer in the feature extraction module. When l=1, ; , , and They represent the weight parameter matrix of the lth layer in the feature extraction module respectively;

[0025] According to the query vector , key vector Sum value vector Calculate the global correlation features. The specific calculation formula is:

[0026] ;

[0027] in, represents the global correlation feature of the lth layer in the feature extraction module, represents the normalization function, the superscript T represents the matrix transpose, The input dimension of the lth layer in the feature extraction module;

[0028] According to the Gaussian kernel scale parameter Calculate the neighboring association features. The specific calculation formula is:

[0029] ;

[0030] in, represents the neighboring correlation features of the lth layer in the feature extraction module, represents the normalization function, T represents the number of sampling moments in the measurement period, , They respectively represent the sequence numbers of the time points in the measurement period.

[0031] Furthermore, the neural network model also includes a first calculation module, which uses the first calculation module to calculate the symmetric relative entropy of the global correlation feature and the neighboring correlation feature according to the extracted global correlation feature and the neighboring correlation feature. The specific calculation formula is:

[0032] ;

[0033] in, represents the symmetric relative entropy between global correlation features and neighboring correlation features, L represents the number of layers of the feature extraction module, Represents relative entropy.

[0034] Further, calculating the comprehensive abnormal score of the user virtual impedance according to the user virtual impedance and the symmetric relative entropy specifically includes:

[0035] The reconstruction error is calculated according to the user virtual impedance and the symmetric relative entropy. The specific calculation formula is:

[0036] ;

[0037] in, represents the reconstruction error; represents the virtual impedance of the nth user at the tth moment; Indicates the virtual impedance of the nth user at time t The predicted value of represents the loss coefficient; represents the F-norm; represents the 1 norm;

[0038] The abnormal score of the user virtual impedance is calculated according to the symmetric relative entropy and the reconstruction error. The specific calculation formula is:

[0039] ;

[0040] in, represents the abnormal score of the virtual impedance of the nth user at the tth time, represents element-wise product, represents the normalization function;

[0041] The anomaly scores of the user virtual impedance at all times are summed to obtain the comprehensive anomaly score of the user virtual impedance.

[0042] Further, the leakage fault is located according to the comprehensive abnormal score of the user virtual impedance, specifically including:

[0043] The leakage fault location is determined based on the number of anomalies of a single leakage fault location and the comprehensive anomaly score of all user virtual impedances.

[0044] Furthermore, the method further comprises:

[0045] After eliminating the leakage fault according to the leakage fault location result, obtain the residual current of the substation, and judge whether the leakage fault still exists according to the residual current of the substation; if so, repeat the leakage fault location step.

[0046] Based on the same concept, the present invention also provides an electronic device, including a memory, a processor and a computer program / instruction stored in the memory, wherein the processor executes the computer program / instruction to implement the method for locating leakage fault of a branch line before a meter as described above.

[0047] Based on the same concept, the present invention also provides a computer-readable storage medium on which a computer program / instruction is stored. When the computer program / instruction is executed by a processor, the method for locating leakage fault of a branch line before a meter as described above is implemented.

[0048] Based on the same concept, the present invention also provides a computer program product, including a computer program / instruction, which, when executed by a processor, implements the method for locating leakage fault of a branch line before a meter as described above.

[0049] Compared with the prior art, the advantages of the present invention are:

[0050] The leakage fault locating method of the present invention does not require human intervention. The leakage fault of the user branch line can be located in real time only through the total meter voltage, user meter voltage, user meter active power and user meter reactive power, thereby improving the leakage fault locating efficiency and effectively improving the operation and maintenance efficiency of the low-voltage distribution network. The present invention integrates user virtual impedance modeling and global-adjacent two-dimensional correlation feature analysis to accurately capture the user virtual impedance mutation characteristics caused by leakage, thereby significantly improving the leakage detection accuracy. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] In order to more clearly illustrate the technical solution of the present invention, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only one embodiment of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0052] Figure 1 It is a flow chart of a method for locating leakage faults in branch lines before user meters in a substation area according to an embodiment of the present invention; Figure 2 is a heat map of the neighboring correlation characteristics of the virtual impedance of a normal user in an embodiment of the present invention; Figure 3 is a global correlation feature distribution heat map of normal user virtual impedance in an embodiment of the present invention; Figure 4 is a heat map of the neighboring correlation characteristics of the virtual impedance of the leakage user in the embodiment of the present invention; Figure 5 is a global correlation feature distribution heat map of the virtual impedance of the leakage user in the embodiment of the present invention; Figure 6 It is a visualization diagram of the leakage fault location result in the embodiment of the present invention. DETAILED DESCRIPTION

[0053] The following is a clear and complete description of the technical solutions in the present invention in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0054] The technical solution of the present application is described in detail with specific embodiments below. The following specific embodiments can be combined with each other, and the same or similar concepts or processes may not be described in detail in some embodiments.

[0055] Embodiment 1

[0056] In order to solve the problem that the traditional leakage fault location mainly relies on manual experience and traditional instrument measurement, resulting in poor real-time performance and accuracy, and poor intermittent leakage detection and positioning capabilities, the present invention provides a leakage fault location method for branch lines before user meters in a substation area, which is applied to leakage identification and positioning of branch lines before user meters in low-voltage substation areas. Figure 1 As shown, the leakage fault location method includes the following steps:

[0057] Step 1: When leakage fault location is required, obtain the total meter voltage of the substation area, the user meter voltage, the user meter active power, and the user meter reactive power within the measurement period.

[0058] Real-time monitoring of the residual current in the substation area. When the residual current in the substation area is greater than the current threshold, it indicates that there is leakage in the substation area and leakage fault location is required. The current threshold is the tripping threshold of the leakage protection element. The current threshold in this embodiment is 300mA.

[0059] During the measurement period, according to the sampling interval, the user information collection system is used to collect the total meter voltage, user meter voltage, user meter active power and user meter reactive power of the area, that is, the total meter voltage at each sampling time t within the measurement period T can be obtained. , the voltage of each user's meter n at each sampling time t (i.e., user meter voltage), active power of each user meter n at each sampling time t and reactive power (i.e. the active power of the user's meter and the reactive power of the user's meter).

[0060] Step 2: Based on the total meter voltage , User meter voltage , User meter active power And the user's meter reactive power , solve for the user virtual impedance .

[0061] There is a shortest power supply path between each user's electricity meter and the transformer. The loop impedance on the shortest power supply path is defined as the user's virtual impedance. The specific expression is:

[0062] (1)

[0063] in, represents the virtual impedance of the nth user at the tth moment; represents the nth user virtual resistance at the tth moment; represents the virtual reactance of the nth user at the tth moment, and j represents the imaginary unit. The voltage drop on the shortest power supply path The longitudinal component Add lateral component The form is expressed as:

[0064] (2)

[0065] (3)

[0066] (4)

[0067] in, Indicates the voltage drop of the nth user at the tth moment; It represents the longitudinal component of the voltage drop of the nth user at the tth moment; Represents the transverse component of the voltage drop of the nth user at the tth moment.

[0068] According to formula (1) to formula (4), the total meter voltage estimation function can be constructed with the user virtual impedance as the independent variable and the user meter voltage, user meter active power and user meter reactive power as the dependent variables:

[0069] (5)

[0070] in, It represents the total voltage estimation value corresponding to the nth user at the tth time.

[0071] The problem of solving the user virtual impedance is transformed into an objective function optimization problem of solving the minimum difference between the total meter voltage in the substation area and the estimated value of the total meter voltage. The specific expression is:

[0072] (6)

[0073] Where T represents the number of sampling moments in the measurement period, N represents the number of user meters in the area, represents the objective function. Optimize and solve formula (6) to obtain , , and then according to formula (1) the user virtual impedance can be solved . Thus, the virtual impedance of each user at each sampling moment in the measurement period can be obtained:

[0074] (7)

[0075] Formula (7) is flattened into a one-dimensional form to facilitate subsequent feature extraction operations.

[0076] Step 3: Extract global correlation features and neighboring correlation features of user virtual impedance.

[0077] In a specific implementation of the present invention, a pre-trained neural network model is used to perform feature extraction, symmetric relative entropy calculation and anomaly score calculation. The neural network model of this embodiment is based on the original Transformer model and adds a feature extraction module, a first calculation module and a second calculation module. The newly added feature extraction module is used to extract global correlation features and neighboring correlation features of the user virtual impedance. The newly added first calculation module is used to calculate the symmetric relative entropy of the global correlation features and the neighboring correlation features according to the extracted global correlation features and neighboring correlation features. The newly added second calculation module is used to calculate the comprehensive anomaly score of the user virtual impedance according to the user virtual impedance and the symmetric relative entropy. The original Transformer model is used to predict the user virtual impedance to obtain the user virtual impedance prediction value. .

[0078] In a specific embodiment of the present invention, a newly added feature extraction module is used to extract global correlation features and neighboring correlation features of user virtual impedance, specifically including:

[0079] Step 3.1: Calculate the query vector, key vector, value vector and Gaussian kernel scale parameter according to the user virtual impedance. The specific calculation formula is:

[0080] (8)

[0081] (9)

[0082] (10)

[0083] (11)

[0084] in, , , and They represent the query vector, key vector, value vector and Gaussian kernel scale parameter of the lth layer in the feature extraction module respectively; Represents the input amount of the lth layer in the feature extraction module. When l=1, ; , , and They respectively represent the weight parameter matrix of the lth layer in the feature extraction module.

[0085] Step 3.2: Based on the query vector , key vector Sum value vector Calculate the global correlation features. The specific calculation formula is:

[0086] (12)

[0087] in, represents the global correlation feature of the lth layer in the feature extraction module, represents the normalization function, the superscript T represents the matrix transpose, It represents the input dimension of the lth layer in the feature extraction module. The global correlation feature uses the multi-head attention mechanism in the original Transformer model to quantify the inherent electrical interconnection properties between different users.

[0088] Step 3.3: According to the Gaussian kernel scale parameter Calculate the neighboring association features. The specific calculation formula is:

[0089] (13)

[0090] in, represents the neighboring correlation features of the lth layer in the feature extraction module, represents the normalization function, T represents the number of sampling moments in the measurement period, They represent the sequence numbers of the time points in the measurement cycle. The neighboring association feature uses an adaptive Gaussian kernel to capture the neighboring leakage users that are most affected by the leakage point.

[0091] Figure 2 and Figure 3 The neighboring correlation feature distribution and global correlation feature distribution of normal user virtual impedance are shown respectively. Figure 4 and Figure 5 The neighboring correlation feature distribution and the global correlation feature distribution of the virtual impedance of the leakage user are shown respectively, where the numbers represent the correlation feature values, and the closer to 1, the greater the correlation.

[0092] The pre-training process of the neural network model is similar to the leakage fault location process. The weight parameter matrix of the feature extraction module is determined through pre-training. , , and As well as the parameters in the original Transformer model. In order to obtain training samples, a real distribution network laboratory was used for simulation. There were 117 simulated users in the simulation area (i.e., N=117), the measurement period was 1 month, and the sampling interval was 15 min / time (T=2880). The total meter voltage, simulated user meter voltage, simulated user meter active power, and simulated user meter reactive power were obtained. According to the scenario classification in Table 1, each scenario was divided into training set and test set in a ratio of 7:3, and the total number of samples was 24750.

[0093] Table 1 All leakage scenarios included in the sample

[0094]

[0095] The pre-training of the neural network model includes the training of the original Transformer model and the training of the feature extraction module. The original Transformer model is trained first, and then the feature extraction module is trained. The training process of the original Transformer model is:

[0096] Calculate the simulated user virtual impedance according to the total meter voltage, the simulated user meter voltage, the simulated user meter active power and the simulated user meter reactive power (such as formula (5) and formula (6));

[0097] The original Transformer model is trained and tested with the total meter voltage, simulated user meter voltage, simulated user meter active power and simulated user meter reactive power as input quantities and the calculated simulated user virtual impedance as the true label.

[0098] The original Transformer model after training and testing can be used to predict the user's virtual impedance and obtain the user's virtual impedance prediction value. After completing the training and testing of the original Transformer model, the feature extraction module is trained. The specific training process of the feature extraction module is as follows:

[0099] Extract global correlation features and neighboring correlation features from the calculated simulated user virtual impedance; calculate the symmetric relative entropy of global correlation features and neighboring correlation features based on the extracted global correlation features and neighboring correlation features; calculate the comprehensive anomaly score of the simulated user virtual impedance based on the simulated user virtual impedance and the symmetric relative entropy; determine the predicted leakage simulated user based on the comprehensive anomaly score of the simulated user virtual impedance; reversely adjust the parameters in the feature extraction module based on the predicted leakage simulated user and the real leakage simulated user to achieve pre-training of the feature extraction module. If the number of predicted leakage simulated users is equal to the number of anomalies located in a single leakage fault, and the predicted leakage simulated users are the same as the real leakage simulated users, it indicates that the prediction is correct; otherwise, adjust the parameters in the feature extraction module. Table 2 shows the parameter settings during model training.

[0100] Table 2 Model training parameter settings

[0101]

[0102] After the training is completed, the test set is used to evaluate the neural network model. The evaluation index results are shown in Table 3.

[0103] Table 3 Results of neural network model evaluation indicators

[0104]

[0105] Step 4: Calculate the symmetric relative entropy of the global correlation features and the neighboring correlation features based on the extracted global correlation features and neighboring correlation features.

[0106] The symmetric relative entropy is used to quantify the feature differences between normal users and leakage users. The first calculation module is used to calculate the symmetric relative entropy of the global correlation feature and the neighboring correlation feature according to the extracted global correlation feature and neighboring correlation feature. The specific calculation formula is:

[0107] (14)

[0108] in, represents the symmetric relative entropy between global correlation features and neighboring correlation features, L represents the number of layers of the feature extraction module, Represents relative entropy.

[0109] Step 5: Calculate the comprehensive anomaly score of the user virtual impedance based on the user virtual impedance and the symmetric relative entropy.

[0110] In a specific embodiment of the present invention, the second calculation module is used to calculate the comprehensive abnormal score of the user virtual impedance according to the user virtual impedance and the symmetric relative entropy, specifically including:

[0111] Step 5.1: Calculate the reconstruction error based on the user virtual impedance and the symmetric relative entropy. The specific calculation formula is:

[0112] (15)

[0113] in, represents the reconstruction error; Indicates the virtual impedance of the nth user at time t The predicted value of Predicted by the original Transformer model; represents the loss coefficient; represents the F-norm, Represents the 1-norm.

[0114] The reconstruction error calculation enhances the model's ability to extract neighboring correlation features and global correlation features.

[0115] Step 5.2: Calculate the abnormal score of the user virtual impedance based on the symmetric relative entropy and the reconstruction error. The specific calculation formula is:

[0116] (16)

[0117] in, represents the abnormal score of the virtual impedance of the nth user at the tth time, represents element-wise product, Represents the normalized function. The anomaly score of the virtual impedance of each user at each moment combines the symmetric relative entropy and the reconstruction error, which improves the scoring accuracy.

[0118] Step 5.3: Sum the abnormal scores of the user virtual impedance at all times to obtain the comprehensive abnormal score of the user virtual impedance.

[0119] For each user, the anomaly scores of the virtual impedance at all times are summed up to obtain the total anomaly score of the user's virtual impedance, and the leakage location problem of the branch line before the meter is converted into a time series anomaly detection problem.

[0120] Step 6: Locate the leakage fault based on the comprehensive abnormal score of the user's virtual impedance.

[0121] In a specific implementation of the present invention, the leakage fault location is determined according to the number of anomalies of a single leakage fault location and the comprehensive anomaly score of all user virtual impedances.

[0122] The number of abnormalities in a single leakage fault location refers to the number of users with leakage faults determined during a single leakage fault location. When leakage fault location is required, there is at least one user with leakage faults, so the number of abnormalities in a single leakage fault location is greater than or equal to 1.

[0123] For example, the number of users in the substation is 117, and the number of abnormalities in a single leakage fault location is 1, which means that when the leakage fault location method of the present invention is used to locate a leakage fault, one leakage fault user needs to be located. Therefore, the maximum comprehensive abnormality score is selected from the comprehensive abnormality scores of all user virtual impedances, and the user front branch line corresponding to the maximum comprehensive abnormality score is the leakage fault location, such as Figure 6 shown. Figure 6 In the example, when the number of abnormalities in a single leakage fault location is 1, the method of the present invention can be used to calculate the comprehensive abnormality scores of the virtual impedances of 117 users, and locate the branch line before the meter of the 40th user as the leakage fault location; Figure 6 The threshold δ in is determined according to the number of anomalies located in a single leakage fault and the comprehensive anomaly score of all user virtual impedances. The threshold δ in this embodiment is 0.21, that is, when the comprehensive anomaly score of the user virtual impedance exceeds the threshold δ, the corresponding user meter front branch line is the leakage fault location.

[0124] Step 7: After eliminating the leakage fault according to the leakage fault location result of step 6, obtain the residual current of the substation, and determine whether the leakage fault still exists based on the residual current of the substation. If so, repeat steps 1 to 7.

[0125] Embodiment 2

[0126] An embodiment of the present invention also provides an electronic device, which includes: a memory, a processor and a computer program / instructions stored in the memory, and the processor executes the computer program / instructions to implement the method for locating leakage faults in branch lines before user meters in a substation in an embodiment of the present application.

[0127] Although not shown, the electronic device includes a processor that can perform various appropriate operations and processes according to the programs and / or data stored in the read-only memory (ROM) or the programs and / or data loaded from the storage part into the random access memory (RAM). The processor can be a multi-core processor or can include multiple processors. In some embodiments, the processor can include a general-purpose main processor and one or more special coprocessors, such as a central processing unit, a graphics processing unit (GPU), a neural network processor (NPU), a digital signal processor (DSP), etc. In RAM, various programs and data required for device operation are also stored. The processor, ROM, and RAM are connected to each other via a bus. The input / output (I / O) interface is also connected to the bus.

[0128] The processor and memory are used together to execute the program / instructions stored in the memory. When the program / instructions are executed by the computer, the methods, steps or functions described in the above embodiments can be implemented.

[0129] Although not shown, an embodiment of the present invention further provides a computer-readable storage medium having a computer program / instruction stored thereon, which, when executed by a processor, implements the method for locating leakage faults in branch lines before user meters in a substation in an embodiment of the present application.

[0130] Readable storage media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer-readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer-readable media does not include temporary computer-readable media (transitory media), such as modulated data signals and carrier waves.

[0131] Although not shown, an embodiment of the present invention further provides a computer program product, including: a computer program / instruction, which, when executed by a processor, implements the method for locating leakage faults in branch lines before user meters in a substation in an embodiment of the present application.

[0132] What is disclosed above is only a specific implementation mode of the present invention, but the protection scope of the present invention is not limited thereto. Any technician familiar with the technical field can easily think of changes or modifications within the technical scope disclosed in the present invention, which should be covered within the protection scope of the present invention.

Claims

1. A method for locating leakage fault of branch line before meter, characterized in that: The fault location method comprises: Obtain the total meter voltage of the area, the user meter voltage, the user meter active power and the user meter reactive power within the measurement period; Based on the total meter voltage, user meter voltage, user meter active power and user meter reactive power, solve the user virtual impedance; Extracting global correlation features and neighboring correlation features from the user virtual impedance; The symmetric relative entropy of the global correlation feature and the neighboring correlation feature is calculated according to the extracted global correlation feature and the neighboring correlation feature; Calculate a comprehensive abnormality score of the user virtual impedance according to the user virtual impedance and the symmetric relative entropy; The leakage fault is located according to the comprehensive abnormal score of the user virtual impedance.

2. The method for locating leakage fault of branch line before the meter according to claim 1 is characterized in that: Based on the total meter voltage, user meter voltage, user meter active power and user meter reactive power, the user virtual impedance is solved, including: Based on the user meter voltage, user meter active power and user meter reactive power, the total meter voltage estimation function is constructed. The specific expression is: ； in, It represents the estimated value of the total meter voltage corresponding to the nth user at the tth time, represents the active power of the nth user's meter at the tth moment, represents the reactive power of the nth user's meter at the tth moment, represents the nth user virtual resistance at the tth moment, represents the virtual reactance of the nth user at the tth moment, Indicates the voltage of the nth user's electric meter at the tth moment; According to the total meter voltage estimation function and the total meter voltage, the user virtual impedance is solved; wherein the solution function is: ； ； Where T represents the number of sampling moments in the measurement period, N represents the number of user meters in the area, represents the total meter voltage at the tth moment, It represents the virtual impedance of the nth user at the tth moment.

3. The method for locating leakage fault of branch line before the meter according to claim 1 is characterized in that: Extracting global correlation features and neighboring correlation features of the user virtual impedance using a pre-trained neural network model; The neural network model includes a feature extraction module, and the feature extraction module is used to extract global correlation features and neighboring correlation features of the user virtual impedance. The specific process includes: According to the user virtual impedance, the query vector, key vector, value vector and Gaussian kernel scale parameter are calculated. The specific calculation formula is: ， ， ， ； in, , , and They represent the query vector, key vector, value vector and Gaussian kernel scale parameter of the lth layer in the feature extraction module respectively; Represents the input amount of the lth layer in the feature extraction module. When l=1, ; , , and They represent the weight parameter matrix of the lth layer in the feature extraction module respectively; According to the query vector , key vector Sum value vector Calculate the global correlation features. The specific calculation formula is: ； in, represents the global correlation feature of the lth layer in the feature extraction module, represents the normalization function, the superscript T represents the matrix transpose, Represents the input dimension of the lth layer in the feature extraction module; According to the Gaussian kernel scale parameter Calculate the neighboring association features. The specific calculation formula is: ； in, represents the neighboring correlation features of the lth layer in the feature extraction module, represents the normalization function, T represents the number of sampling moments in the measurement period, , They respectively represent the sequence numbers of the time points in the measurement period.

4. The method for locating leakage fault of branch line before the meter according to claim 3 is characterized in that: The neural network model also includes a first calculation module, which uses the first calculation module to calculate the symmetric relative entropy of the global correlation feature and the neighboring correlation feature according to the extracted global correlation feature and the neighboring correlation feature. The specific calculation formula is: ； in, represents the symmetric relative entropy between global correlation features and neighboring correlation features, L represents the number of layers of the feature extraction module, Represents relative entropy.

5. The method for locating leakage fault of branch line before the meter according to claim 1, characterized in that: Calculating a comprehensive abnormality score of the user virtual impedance according to the user virtual impedance and the symmetric relative entropy includes: The reconstruction error is calculated according to the user virtual impedance and the symmetric relative entropy. The specific calculation formula is: ； in, represents the reconstruction error; represents the virtual impedance of the nth user at the tth moment; Indicates the virtual impedance of the nth user at time t The predicted value of represents the loss coefficient; represents the F-norm; represents the symmetric relative entropy between the global correlation feature and the neighboring correlation feature, represents the 1 norm; The abnormal score of the user virtual impedance is calculated according to the symmetric relative entropy and the reconstruction error. The specific calculation formula is: ； in, represents the abnormal score of the virtual impedance of the nth user at the tth time, represents element-wise product, represents the normalization function; The anomaly scores of the user virtual impedance at all times are summed to obtain the comprehensive anomaly score of the user virtual impedance.

6. The method for locating leakage fault of branch line before the meter according to claim 1 is characterized in that: The leakage fault is located according to the comprehensive abnormal score of the user virtual impedance, specifically including: The leakage fault location is determined based on the number of anomalies of a single leakage fault location and the comprehensive anomaly score of all user virtual impedances.

7. The method for locating leakage fault of a branch line before a meter according to any one of claims 1 to 6, characterized in that: The fault location method further includes: After eliminating the leakage fault according to the leakage fault location result, obtain the residual current of the substation, and judge whether the leakage fault still exists according to the residual current of the substation; if so, repeat the leakage fault location step.

8. An electronic device comprising a memory, a processor, and a computer program / instruction stored in the memory, characterized in that: The processor executes the computer program / instructions to implement the method for locating leakage faults of branch lines before the meter as described in any one of claims 1 to 7.

9. A computer-readable storage medium having a computer program / instruction stored thereon, characterized in that: When the computer program / instruction is executed by a processor, the method for locating leakage fault of a branch line before a meter as described in any one of claims 1 to 7 is implemented.

10. A computer program product comprising a computer program / instructions, characterized in that When the computer program / instruction is executed by a processor, the method for locating leakage fault of a branch line before a meter as described in any one of claims 1 to 7 is implemented.

Images