Subway stray current leakage point determination method and device

Abstract

The invention provides a metro stray current leakage point determination method and device. The metro stray current leakage point determination method comprises the steps that train monitoring information, station monitoring information, station mileage information and line information in a preset scene are collected; constructing a subway resistance network model according to the train monitoring information, the station monitoring information, the station mileage information and the line information; constructing a sensitivity increment matrix equation set about a steel rail to-ground transition resistance change value and a steel rail potential change value through the subway resistance network model; and based on the sensitivity increment matrix equation set, rail transition resistance distribution is determined through an iterative calculation method, and a stray current leakage area is determined according to a calculation result. Manual field measurement is not needed, and the subway stray current leakage area is determined through on-line monitoring data.

Description

Technical Field

[0001] The present application relates to the technical field of subway rail insulation deterioration area positioning, and in particular to a subway stray current leakage point determination method. The present application also relates to a subway stray current leakage point determination device, a computing device, and a computer-readable storage medium. Background Art

[0002] With the rapid growth of population in large cities, urban traffic congestion and air pollution are becoming increasingly prominent. To meet this challenge, urban rail transit systems have become the first choice and have developed rapidly due to their environmental friendliness, low energy consumption, high speed and large passenger capacity. However, one problem that cannot be ignored is that due to the limitations of on-site construction, the rails cannot be completely insulated from the ground. With the long-term operation of the rail transit system, the rail-to-ground insulation performance has been deteriorated to varying degrees due to the humid operating environment, aging of insulation fasteners, and accumulation of iron filings and dust on the rails. Therefore, the negative impact of the stray current problem has gradually emerged. When the rail insulation layer deteriorates locally, the stray current leakage will increase significantly, and the stray current will accelerate the corrosion of metal pipes and structural steel bars around the subway through these weaknesses. Therefore, while we enjoy the convenience brought by the subway, we must face up to and take active measures to prevent the electrochemical corrosion caused by stray currents and protect the safety of underground facilities.

[0003] Any local fault in the insulation layer between the subway track and the ground will accelerate the leakage of stray current. Once the leakage exceeds the threshold, it is necessary to intervene with the help of a drainage device, but this may cause the rail potential to rise, leading to an unfavorable cycle effect. Therefore, compared with relying solely on passive protection strategies to control the hazards caused by stray currents, taking an active strategy, that is, accurately identifying and quickly repairing the damaged locations of track insulation, should be a more effective strategy. This method can not only directly reduce the leakage of stray currents, but also reduce the corrosion damage and DC bias magnetism caused by stray currents from the source, thereby more efficiently ensuring the safety and long-term stability of the subway system.

[0004] In the prior art, there are two methods for locating insulation damage points between the track and the ground. First, by analyzing the line with local insulation damage, extracting relevant characteristic values, and realizing the positioning target of the damage point or area, however, this method regards the transition resistance as a fixed value for analysis, and it is difficult to analyze and locate the line in other states. This method basically takes the line with a single damaged area or damaged point as the research object, and cannot realize the location and evaluation of multiple damages on the line. Second, the subway transition resistance measurement method stipulated in the subway regulations, however, most of these methods are offline measurements, requiring workers to measure relevant data on the line during the train stop period, which is costly but inefficient and cannot reflect the real-time state of track insulation. In addition, there are some online measurement methods, however, this method requires the installation of other electrical facilities in the subway system, the economic cost is high, and it is difficult to calculate and evaluate the transition resistance state of the entire line, and it is further difficult to realize the location of the local insulation damage point of the subway track to the ground. Summary of the invention

[0005] In view of this, the embodiment of the present application provides a method for determining a subway stray current leakage point to solve the technical defects in the prior art. The embodiment of the present application also provides a subway stray current leakage point determination device, a computing device, and a computer-readable storage medium.

[0006] According to a first aspect of an embodiment of the present application, a method for determining a subway stray current leakage point is provided, comprising:

[0007] S1, collect train monitoring information, station monitoring information, station mileage information and line information under preset scenarios;

[0008] S2, constructing a subway resistance network model according to the train monitoring information, the station monitoring information, the station mileage information and the line information;

[0009] S3, constructing a group of sensitivity increment matrix equations about the change value of rail-to-ground transition resistance and the change value of rail potential through the subway resistance network model;

[0010] S4, based on the sensitivity increment matrix equation group, determine the rail transition resistance distribution through an iterative calculation method, and determine the stray current leakage area according to the calculation result.

[0011] Optionally, the S1 includes:

[0012] When there is only one train running between two adjacent traction substations on the line, S11, collecting train monitoring information of the train at any time, including train position information and train traction current information;

[0013] S12, collecting the station monitoring information including the power supply current information of the traction substation and the rail potential information of the station and the traction substation;

[0014] S13, collecting the station mileage information including the station location information;

[0015] S14, including the line information of the longitudinal resistance value of the rail per unit length and the transition resistance value of the rail per unit length to the ground.

[0016] Optionally, the S2 includes:

[0017] S21, dividing the line between two adjacent traction substations into discrete equidistant divisions according to a preset division accuracy to obtain a plurality of nodes, wherein the nodes include traction substation nodes, station nodes and common nodes;

[0018] S22, constructing an equivalent series resistance structure according to the longitudinal resistance value of the rail per unit length and the division accuracy;

[0019] S23, constructing a parallel resistance structure according to the rail-to-ground transition resistance value per unit length and the division accuracy;

[0020] S24, constructing a first equivalent current source according to the power supply current information of the traction substation;

[0021] S25, constructing a second equivalent current source according to the train position information and the train traction current information;

[0022] S26, constructing the subway resistance network model based on the equivalent series resistance structure, the parallel resistance structure, the first equivalent current source and the second equivalent current source.

[0023] Optionally, the S3 includes:

[0024] S31, when the train travels to any of the nodes, determining the rail potential information of the ideal station and the traction substation according to the station location information and the subway resistance network model;

[0025] S32, according to the station location information, subtract the rail potential information of the ideal station and traction substation from the corresponding rail potential information of the station and traction substation to construct a rail potential difference column vector;

[0026] S33, calculating a sensitivity matrix associating the rail transition resistance change value and the rail potential change value, wherein the rail potential change value is associated with the station node and the traction substation node;

[0027] S34, repeatedly executing S31-S33 until the train travels to a plurality of different nodes, corresponding to a plurality of rail potential difference column vectors and a plurality of sensitivity matrices;

[0028] S35, constructing the sensitivity increment matrix equation group according to the multiple rail potential difference column vectors and the multiple sensitivity matrices.

[0029] Optionally, the S35 includes:

[0030] S351, arranging the plurality of rail potential difference column vectors according to the order in which the trains travel to the nodes to obtain an augmented rail potential difference column vector;

[0031] S352, arranging the plurality of sensitivity matrices according to the order in which the trains travel to the nodes to obtain an augmented sensitivity matrix;

[0032] S353, setting the rail-to-ground transition resistance increment column vector;

[0033] S354, constructing the sensitivity increment matrix equation group according to the augmented rail potential difference column vector, the augmented sensitivity matrix and the rail-to-ground transition resistance increment column vector.

[0034] Optionally, in S4, determining the rail transition resistance distribution by an iterative calculation method based on the sensitivity increment matrix equation group includes:

[0035] S41, setting an initial rail transition resistance value according to the rail-to-ground transition resistance value per unit length and the division accuracy;

[0036] S42, calculating an initial rail potential initial value vector by a node voltage method based on the rail potential information of the station and the traction substation and the initial rail transition resistance value;

[0037] S43, setting the initial value of the rail transition resistance increment to 0, and setting T=1, where T is a positive integer;

[0038] S44, calculating the rail potential difference column vector and the sensitivity matrix of the T-th iteration according to the sensitivity increment matrix equation group, and constructing the sensitivity matrix equation group of the T-th iteration based on the calculation results;

[0039] S45, based on the sensitivity matrix equation group of the T-th iteration, calculate the rail transition resistance increment of the T-th iteration by the least square method, and sum the calculation result with the rail transition resistance value of the T-1-th iteration to obtain the rail transition resistance value of the T-th iteration;

[0040] S46, updating the rail potential difference column vector and the sensitivity matrix, T is incremented by 1, and S43 is executed until the rail transition resistance increment is less than a preset convergence accuracy, thereby obtaining the rail transition resistance distribution.

[0041] Optionally, in S4, determining the stray current leakage area according to the calculation result includes:

[0042] S47, determining a transition resistance threshold value according to the transition resistance value of the rail to the ground per unit length and the division accuracy;

[0043] S48, comparing the rail transition resistance distribution with the transition resistance threshold, determining a transition resistance abnormal area, and using the transition resistance abnormal area as the stray current leakage area.

[0044] According to a second aspect of an embodiment of the present application, a device for determining a subway stray current leakage point is provided, comprising:

[0045] A collection module is configured to collect train monitoring information, station monitoring information, station mileage information and line information under preset scenarios;

[0046] A first construction module is configured to construct a subway resistance network model according to the train monitoring information, the station monitoring information, the station mileage information and the line information;

[0047] The second construction module is configured to construct a sensitivity increment matrix equation group about the change value of the rail-to-ground transition resistance and the change value of the rail potential through the subway resistance network model;

[0048] The determination module is configured to determine the rail transition resistance distribution by an iterative calculation method based on the sensitivity increment matrix equation group, and determine the stray current leakage area according to the calculation result.

[0049] According to a third aspect of an embodiment of the present application, a computing device is provided, including:

[0050] Memory and processor;

[0051] The memory is used to store computer executable instructions, and the processor implements the steps of the method for determining subway stray current leakage points when executing the computer executable instructions.

[0052] According to a fourth aspect of an embodiment of the present application, a computer-readable storage medium is provided, which stores computer-executable instructions, and when the instructions are executed by a processor, the steps of the method for determining a subway stray current leakage point are implemented.

[0053] According to a fifth aspect of an embodiment of the present application, a chip is provided, which stores a computer program, and when the computer program is executed by the chip, the steps of the method for determining a subway stray current leakage point are implemented.

[0054] The method for determining subway stray current leakage points provided by the present application collects train monitoring information, station monitoring information, station mileage information and line information under preset scenarios; constructs a subway resistance network model based on the train monitoring information, station monitoring information, station mileage information and line information; constructs a sensitivity increment matrix equation group about the change value of rail-to-ground transition resistance and rail potential change value through the subway resistance network model; based on the sensitivity increment matrix equation group, determines the rail transition resistance distribution through an iterative calculation method, and determines the stray current leakage area based on the calculation results. It realizes the determination of subway stray current leakage areas through online monitoring data without manual on-site measurement. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0056] Figure 1 It is a flow chart of a method for determining subway stray current leakage points provided by an embodiment of the present application;

[0057] Figure 2 It is a subway system model diagram of a method for determining subway stray current leakage points provided by an embodiment of the present application;

[0058] Figure 3 It is a structural schematic diagram of a subway stray current leakage point determination device provided by an embodiment of the present application;

[0059] Figure 4 It is a structural block diagram of a computing device provided in one embodiment of the present application. DETAILED DESCRIPTION

[0060] Many specific details are described in the following description to facilitate a full understanding of the present application. However, the present application can be implemented in many other ways than those described herein, and those skilled in the art can make similar generalizations without violating the connotation of the present application, so the present application is not limited by the specific implementation disclosed below.

[0061] The terms used in one or more embodiments of the present application are only for the purpose of describing specific embodiments, and are not intended to limit one or more embodiments of the present application. The singular forms of "a", "said" and "the" used in one or more embodiments of the present application and the appended claims are also intended to include plural forms, unless the context clearly indicates other meanings. It should also be understood that the term "and / or" used in one or more embodiments of the present application refers to and includes any or all possible combinations of one or more associated listed items.

[0062] It should be understood that, although the terms first, second, etc. may be used to describe various information in one or more embodiments of the present application, such information should not be limited to these terms. These terms are only used to distinguish the same type of information from each other. For example, without departing from the scope of one or more embodiments of the present application, the first may also be referred to as the second, and similarly, the second may also be referred to as the first.

[0063] In the present application, a method for determining a subway stray current leakage point is provided. The present application also relates to a subway stray current leakage point determination device, a computing device, and a computer-readable storage medium, which are described in detail one by one in the following embodiments.

[0064] Figure 1 A flowchart of a method for determining a subway stray current leakage point provided in accordance with an embodiment of the present application is shown, which specifically includes the following steps:

[0065] S1, collect train monitoring information, station monitoring information, station mileage information and line information under preset scenarios;

[0066] S2, constructing a subway resistance network model according to the train monitoring information, the station monitoring information, the station mileage information and the line information;

[0067] S3, constructing a group of sensitivity increment matrix equations about the change value of rail-to-ground transition resistance and the change value of rail potential through the subway resistance network model;

[0068] S4, based on the sensitivity increment matrix equation group, determine the rail transition resistance distribution through an iterative calculation method, and determine the stray current leakage area according to the calculation result.

[0069] Among them, train monitoring information refers to information including train position and train electrical data during train operation; station monitoring information refers to information including electrical data related to the station and traction substation; station mileage information refers to information including the location of the station on the train line; and line information refers to information including parameters related to rails on the train line.

[0070] Based on this, a bilateral power supply rail-earth structure resistance network model, namely the subway resistance network model, is constructed through train monitoring information, station monitoring information, station mileage information and line information to characterize the equivalent circuit relationship formed by the current provided by the train and traction substation through the rail to the earth. After obtaining the subway resistance network model, the sensitivity increment matrix equation group is determined by analyzing the correlation between the change value of the rail-to-ground transition resistance and the change value of the rail potential, and the rail transition resistance distribution is calculated by iterative calculation method based on the sensitivity increment matrix equation group. The calculation results of the rail transition resistance distribution are analyzed to determine the area corresponding to the abnormal data, and then the stray current leakage area is determined.

[0071] Further, regarding the process of collecting train monitoring information, station monitoring information, station mileage information and line information under a preset scenario in step S1, in this embodiment, the specific implementation method is as follows:

[0072] When there is only one train running between two adjacent traction substations on the line, S11, collects the train monitoring information of the train at any time, including the train position information and the train traction current information; S12, collects the station monitoring information including the traction substation power supply current information, the station and traction substation rail potential information; S13, collects the station mileage information including the station location information; S14, collects the line information including the longitudinal resistance value of the unit length of the rail and the transition resistance value of the unit length of the rail to the ground.

[0073] Among them, the station location information represents the topological structure of the train line, and the longitudinal resistance value per unit length of the rail and the transition resistance value per unit length of the rail to the ground are determined by the actual usage scenario, such as setting with reference to the "CJJT 49-2020 Technical Standard for Corrosion Protection of Metro Stray Current", etc. The specific setting results are not limited in this embodiment.

[0074] Furthermore, in step S2, the process of constructing a subway resistance network model according to the train monitoring information, the station monitoring information, the station mileage information and the line information is specifically implemented as follows in this embodiment:

[0075] S21, divide the line between two adjacent traction substations into discrete and equidistant parts according to a preset division accuracy to obtain a plurality of nodes, wherein the nodes include traction substation nodes, station nodes and ordinary nodes; S22, construct an equivalent series resistance structure according to the longitudinal resistance value of the unit length rail and the division accuracy; S23, construct a parallel resistance structure according to the transition resistance value of the unit length rail to the ground and the division accuracy; S24, construct a first equivalent current source according to the power supply current information of the traction substation; S25, construct a second equivalent current source according to the train position information and the train traction current information; S26, construct the subway resistance network model based on the equivalent series resistance structure, the parallel resistance structure, the first equivalent current source and the second equivalent current source.

[0076] Among them, in actual use scenarios, if the total length between two adjacent traction substations is X, the number of equally spaced nodes is N = n + 1, then the division accuracy is L = X / n, the total number of resistance branches is B = 2n + 1, and the distances between the train and the traction substations at both ends are L and 1 and L 2 After constructing the subway resistance network model, the current source vector of the train at different node locations is recorded. j is the node position number of the train, and Supply current to the traction substation, I j is the train traction current, and the relationship between the traction substation power supply current and the train traction current is: and

[0077] Record the station node location and traction substation node location information on the line D = [D 1 D 2 …D k ], station rail potential information Where k is the total number of station nodes and traction substation nodes. At the same time, the line related parameters are obtained, including the longitudinal resistance of the rail per unit length R r , rail transition resistance R s .

[0078] Then, the rail texture is set to be uniform, and the line between two adjacent traction substations is discretely divided into n equidistant intervals according to a certain accuracy and n+1 nodes are marked, including traction substation nodes, station nodes and ordinary nodes; the rail is equivalent to a series resistance structure according to the longitudinal resistance of the rail combined with the accuracy, and the number of longitudinal resistances is n; then the rail-to-ground resistance is equivalent to a parallel resistance structure according to the rail transition resistance combined with the accuracy, and the number of longitudinal resistances is n+1; the power supply current of the traction substation is equivalent to a current source, and the traction current of the train at different node positions is equivalent to a current source, and a bilateral power supply rail-earth structure resistance network model is constructed, as shown in the following figure. Figure 2 A subway system model diagram of a method for determining subway stray current leakage points is shown in the provided diagram.

[0079] Specifically, read the current source vector of the train at different node positions Based on the bilateral power supply rail-earth structure resistance network model, the node admittance matrix Y is formed N , Y N Each matrix element in is related to the subway line resistance / conductance value of the corresponding node, and the subway line resistance / conductance value of any node is associated with the set positioning accuracy, that is, the division accuracy. Specifically, according to the positioning accuracy L, the longitudinal conductance Y of the rails between nodes can be obtained. g =1 / (R r *L) and rail transition conductance Y d =1 / (R s / L), substitute the longitudinal conductance of the rails between nodes and the transition conductance of the rails into the branch admittance matrix Y b The corresponding position in the branch admittance matrix Y is obtained b ; According to the topological relationship of the node network, the node network association matrix A is established; the node admittance matrix Y is obtained from the branch admittance matrix and the node network association matrix N =AY b -1 A T ; Rail potentials at all nodes when the train is at different node positions

[0080] Furthermore, in step S3, the process of constructing a group of sensitivity increment matrix equations about the change value of rail-to-ground transition resistance and rail potential change value through the subway resistance network model is specifically implemented as follows in this embodiment:

[0081] S31, when the train reaches any of the nodes, the ideal station and traction substation rail potential information is determined according to the station location information and the subway resistance network model; S32, according to the station location information, the ideal station and traction substation rail potential information is subtracted from the corresponding station and traction substation rail potential information to construct a rail potential difference column vector; S33, calculate a sensitivity matrix associating the rail transition resistance change value and the rail potential change value, wherein the rail potential change value is associated with the station node and the traction substation node; S34, repeat S31-S33 until the train reaches multiple different nodes, corresponding to multiple rail potential difference column vectors and multiple sensitivity matrices; S35, construct the sensitivity increment matrix equation group according to multiple rail potential difference column vectors and multiple sensitivity matrices.

[0082] Furthermore, in step S35, based on the sensitivity increment matrix equation group, the process of determining the rail transition resistance distribution by an iterative calculation method is specifically implemented as follows in this embodiment:

[0083] S351, according to the order in which the trains travel to the nodes, the plurality of rail potential difference column vectors are arranged to obtain an augmented rail potential difference column vector; S352, according to the order in which the trains travel to the nodes, the plurality of sensitivity matrices are arranged to obtain an augmented sensitivity matrix; S353, a rail-to-ground transition resistance increment column vector is set; S354, according to the augmented rail potential difference column vector, the augmented sensitivity matrix and the rail-to-ground transition resistance increment column vector, the sensitivity increment matrix equation group is constructed.

[0084] Using the above example, the parameter identification method is used to reflect the influence of the resistance change of each branch on the rail potential at each node when the train is running at a certain node. The total differential of the rail potential is calculated. Since the train traction current has nothing to do with the resistance value of each branch in the resistor network model, I N With each branch resistance R i It doesn't matter. The above formula can be transformed into By Y N -1 Y N =1 (1 is the unit matrix), we can get So,

[0085] Because U N The change of is related to each resistor in the rail transition resistance branch, so the change of resistance of any rail transition resistance branch multiplied by its sensitivity is U N The change in in VU is the sensitivity of rail potential change caused by branch transition resistance change. The differential equation is written in the form of an incremental matrix to obtain N =MVR, where M is the sensitivity matrix, and the elements in M It represents the sensitivity of the resistance change of the cth rail transition resistance branch to the rail potential change at the station node i when the train is running at the node v, where R is the transition resistance to be determined, let C i(c-n) is the (cn)th element in the i-th row of C.

[0086] Set the stray current leakage area, that is, set the transition resistance at a certain node to be smaller, that is, the branch admittance matrix Y b The branch admittance at the corresponding position in the equation becomes larger and we get Y bp , get the node admittance matrix Y after leakage Np =AY bp -1 A T ;

[0087] Obtain the rail potential of the station node and the traction substation nodes at both ends when the train is running at different node positions, and obtain the voltage difference under the ideal running state Where k is the total number of station nodes and traction substations at both ends, j is the node number where the train is located, and the elements in the column vector The voltage difference when the train is running at different node positions (the nodes where the train is located are numbered from 1 to N) is arranged into a column according to the order of the nodes where the train is located, and the augmented rail potential difference column vector is obtained. ΔU N =MΔR is written as the expanded equation system

[0088]

[0089] Furthermore, in step S4, based on the sensitivity increment matrix equation group, the process of determining the rail transition resistance distribution by an iterative calculation method is specifically implemented as follows in this embodiment:

[0090] S41, according to the unit length rail to ground transition resistance value and the division accuracy, set the initial rail transition resistance value; S42, based on the rail potential information of the station and traction substation and the initial rail transition resistance value, calculate the initial rail potential initial value vector by the node voltage method; S43, set the initial value of the rail transition resistance increment to 0, and set T = 1, and T is a positive integer; S44, according to the sensitivity increment matrix equation group, calculate the rail potential difference column vector and the sensitivity matrix of the Tth iteration, and based on the calculation results Construct the sensitivity matrix equation group of the T-th iteration; S45, based on the sensitivity matrix equation group of the T-th iteration, calculate the rail transition resistance increment of the T-th iteration by the least square method, and sum the calculation result with the rail transition resistance value of the T-1-th iteration to obtain the rail transition resistance value of the T-th iteration; S46, update the rail potential difference column vector and the sensitivity matrix, T is incremented by 1, and S43 is executed until the rail transition resistance increment is less than the preset convergence accuracy to obtain the rail transition resistance distribution.

[0091] Continuing with the above example, since the above equations are incremental matrix equations; and ΔU n There is a nonlinear relationship between and ΔR, that is, the elements in M ​​are not constants, so the augmented fault diagnosis equation cannot be solved directly by linear methods. In order to solve this nonlinear problem, an iterative measure is taken. In order to simplify the description, the subsequent X T Represents the T-th iteration result of X.

[0092] Obtain the rail potential U at the station node and traction substation node at different node positions of the train N According to the initial value of rail transition resistance The node voltage method is used to calculate the initial value vector U of the rail potential N,0 The initial value of each element in the transition resistance increment ΔR is set to 0. Let the iteration counter T = 1 to calculate the rail potential difference column vector ΔU at the Tth iteration. N,T =U N -U N,T-1 and the sensitivity matrix M T , construct the incremental matrix equation ΔU at the Tth iteration N,T =M T ΔR T #. Then use the least squares method to find the Tth iteration to satisfy the formula ΔU N,T =U N -U N,T-1 The optimal solution of the objective function described by ΔR T ,in,

[0093] Minf(ΔR n+1，T , ΔR n+2，T, ..., ΔR 2n+1，T )=||M T ΔR T -ΔU N，T || 2 #.

[0094] Calculate the branch resistance vector R at the Tth iteration T =R T-1 +ΔR T If ||ΔR T || 2 <ε, then exit the iteration, R T This is the final diagnosis result, recorded as R D , where ε is a small positive number, representing the convergence accuracy. If the formula ||ΔR T || 2 <ε is not satisfied, then the R of this iteration result T Get M T+1 , the rail potential vector U is calculated using the node voltage method N,T , return to the above construction ΔU N,T+1 =M T+1 ΔR T+1 During the process, ΔR is calculated T+1 , and so on.

[0095] Furthermore, in step S4, the process of determining the stray current leakage area according to the calculation result is specifically implemented as follows in this embodiment:

[0096] S47, determining a transition resistance threshold value according to the transition resistance value of the rail to the ground per unit length and the division accuracy; S48, comparing the rail transition resistance distribution with the transition resistance threshold value, determining an abnormal transition resistance area, and using the abnormal transition resistance area as the stray current leakage area.

[0097] Using the above example, extract R T The values ​​of the elements in the figure are as follows: R o =R s / L(R s =3Ω·km) is the threshold value. In addition, it can be further limited according to the threshold value. For example, if it is less than 60% of the threshold value, it is identified as an abnormal transition resistance area, and the positioning area accuracy L=X / n.

[0098] Corresponding to the above method embodiment, the present application also provides an embodiment of a device for determining a subway stray current leakage point, Figure 3 FIG. 1 is a schematic diagram showing a structure of a device for determining a subway stray current leakage point provided by an embodiment of the present application. Figure 3 As shown, the device comprises:

[0099] The collection module 302 is configured to collect train monitoring information, station monitoring information, station mileage information and line information under a preset scenario;

[0100] A first construction module 304 is configured to construct a subway resistance network model according to the train monitoring information, the station monitoring information, the station mileage information and the line information;

[0101] The second construction module 306 is configured to construct a sensitivity increment matrix equation group about the change value of the rail-to-ground transition resistance and the change value of the rail potential through the subway resistance network model;

[0102] The determination module 308 is configured to determine the rail transition resistance distribution by an iterative calculation method based on the sensitivity increment matrix equation group, and determine the stray current leakage area according to the calculation result.

[0103] In an optional embodiment, the acquisition module 302 is further configured to:

[0104] When there is only one train running between two adjacent traction substations on the line, the train monitoring information including the train position information and the train traction current information of the train at any time is collected; the station monitoring information including the traction substation power supply current information, the station and traction substation rail potential information is collected; the station mileage information including the station position information is collected; and the line information including the longitudinal resistance value of the unit length of the rail and the transition resistance value of the unit length of the rail to the ground is collected.

[0105] In an optional embodiment, the first building module 304 is further configured to:

[0106] The line between two adjacent traction substations is discretely and equidistantly divided according to a preset division accuracy to obtain a plurality of nodes, wherein the nodes include traction substation nodes, station nodes and ordinary nodes; an equivalent series resistance structure is constructed according to the longitudinal resistance value of the unit length rail and the division accuracy; a parallel resistance structure is constructed according to the transition resistance value of the unit length rail to the ground and the division accuracy; a first equivalent current source is constructed according to the power supply current information of the traction substation; a second equivalent current source is constructed according to the train position information and the train traction current information; and the subway resistance network model is constructed based on the equivalent series resistance structure, the parallel resistance structure, the first equivalent current source and the second equivalent current source.

[0107] In an optional embodiment, the second building module 306 is further configured to:

[0108] When the train reaches any of the nodes, the rail potential information of the ideal station and traction substation is determined according to the station location information and the subway resistance network model; according to the station location information, the rail potential information of the ideal station and traction substation is subtracted from the corresponding rail potential information of the station and traction substation to construct a rail potential difference column vector; a sensitivity matrix associating the rail transition resistance change value and the rail potential change value is calculated, wherein the rail potential change value is associated with the station node and the traction substation node; the second construction module 306 is repeatedly called until the train reaches multiple different nodes, corresponding to multiple rail potential difference column vectors and multiple sensitivity matrices; the sensitivity increment matrix equation group is constructed according to the multiple rail potential difference column vectors and multiple sensitivity matrices.

[0109] In an optional embodiment, the second building module 306 is further configured to:

[0110] According to the order in which the trains travel to the nodes, a plurality of the rail potential difference column vectors are arranged to obtain an augmented rail potential difference column vector; according to the order in which the trains travel to the nodes, a plurality of the sensitivity matrices are arranged to obtain an augmented sensitivity matrix; a rail-to-ground transition resistance increment column vector is set; and according to the augmented rail potential difference column vector, the augmented sensitivity matrix and the rail-to-ground transition resistance increment column vector, the sensitivity increment matrix equation group is constructed.

[0111] In an optional embodiment, the determining module 308 is further configured to:

[0112] According to the unit length rail-to-ground transition resistance value and the division accuracy, the initial rail transition resistance value is set; based on the rail potential information of the station and traction substation and the initial rail transition resistance value, the initial rail potential initial value vector is calculated by the node voltage method; the initial value of the rail transition resistance increment is set to 0, and T=1 is set, and T is a positive integer. According to the sensitivity increment matrix equation group, the rail potential difference column vector and the sensitivity matrix of the Tth iteration are calculated, and the sensitivity matrix equation group of the Tth iteration is constructed based on the calculation results; based on the sensitivity matrix equation group of the Tth iteration, the rail transition resistance increment of the Tth iteration is calculated by the least squares method, and the calculation result is summed with the rail transition resistance value of the T-1th iteration to obtain the rail transition resistance value of the Tth iteration; the rail potential difference column vector and the sensitivity matrix are updated, T is incremented by 1, and the determination module 308 is called until the rail transition resistance increment is less than the preset convergence accuracy to obtain the rail transition resistance distribution.

[0113] In an optional embodiment, the determining module 308 is further configured to:

[0114] A transition resistance threshold is determined according to the transition resistance value of the rail per unit length to the ground and the division accuracy; the rail transition resistance distribution is compared with the transition resistance threshold to determine an abnormal transition resistance area, and the abnormal transition resistance area is used as the stray current leakage area.

[0115] The device for determining subway stray current leakage points provided by the present application collects train monitoring information, station monitoring information, station mileage information and line information under preset scenarios; constructs a subway resistance network model based on the train monitoring information, station monitoring information, station mileage information and line information; constructs a sensitivity increment matrix equation group about the change value of the rail-to-ground transition resistance and the change value of the rail potential through the subway resistance network model; based on the sensitivity increment matrix equation group, determines the rail transition resistance distribution through an iterative calculation method, and determines the stray current leakage area based on the calculation results. It realizes the determination of the subway stray current leakage area through online monitoring data without manual on-site measurement.

[0116] The above is a schematic scheme of a subway stray current leakage point determination device of the present embodiment. It should be noted that the technical scheme of the subway stray current leakage point determination device and the technical scheme of the subway stray current leakage point determination method mentioned above belong to the same concept. For the details not described in detail in the technical scheme of the subway stray current leakage point determination device, please refer to the description of the technical scheme of the subway stray current leakage point determination method mentioned above. In addition, each component in the device embodiment should be understood as a functional module that must be established to implement each step of the program flow or each step of the method, and each functional module is not an actual functional division or separation definition. The device claim defined by such a group of functional modules should be understood as a functional module architecture that mainly implements the solution through the computer program recorded in the specification, and should not be understood as a physical device that mainly implements the solution through hardware.

[0117] Figure 4 The block diagram of a computing device 400 according to an embodiment of the present application is shown. The components of the computing device 400 include but are not limited to a memory 410 and a processor 420. The processor 420 is connected to the memory 410 via a bus 430, and the database 450 is used to store data.

[0118] The computing device 400 also includes an access device 440 that enables the computing device 400 to communicate via one or more networks 460. Examples of these networks include a public switched telephone network (PSTN), a local area network (LAN), a wide area network (WAN), a personal area network (PAN), or a combination of communication networks such as the Internet. The access device 440 may include one or more of any type of network interface (e.g., a network interface card (NIC)) wired or wireless, such as an IEEE 802.11 wireless local area network (WLAN) wireless interface, a World Wide Interoperability for Microwave Access (Wi-MAX) interface, an Ethernet interface, a universal serial bus (USB) interface, a cellular network interface, a Bluetooth interface, a near field communication (NFC) interface, and the like.

[0119] In one embodiment of the present application, the above components of the computing device 400 and Figure 4 Other components not shown in the figure may also be connected to each other, for example, via a bus. It should be understood that Figure 4 The computing device structure block diagram shown is only for the purpose of illustration, and is not intended to limit the scope of the present application. Those skilled in the art may add or replace other components as needed.

[0120] Computing device 400 may be any type of stationary or mobile computing device, including a mobile computer or mobile computing device (e.g., a tablet computer, a personal digital assistant, a laptop computer, a notebook computer, a netbook, etc.), a mobile phone (e.g., a smart phone), a wearable computing device (e.g., a smart watch, smart glasses, etc.), or other types of mobile devices, or a stationary computing device such as a desktop computer or PC. Computing device 400 may also be a mobile or stationary server.

[0121] The processor 420 is used to execute computer executable instructions of each step of the method for determining subway stray current leakage points.

[0122] The above is a schematic scheme of a computing device of this embodiment. It should be noted that the technical scheme of the computing device and the technical scheme of the above-mentioned method for determining subway stray current leakage points belong to the same concept, and the details not described in detail in the technical scheme of the computing device can be referred to the description of the technical scheme of the above-mentioned method for determining subway stray current leakage points.

[0123] An embodiment of the present application further provides a computer-readable storage medium storing computer instructions, which, when executed by a processor, are used to execute the steps of the method for determining a subway stray current leakage point.

[0124] The above is a schematic scheme of a computer-readable storage medium of this embodiment. It should be noted that the technical scheme of the storage medium and the technical scheme of the above-mentioned method for determining subway stray current leakage points belong to the same concept, and the details not described in detail in the technical scheme of the storage medium can be referred to the description of the technical scheme of the above-mentioned method for determining subway stray current leakage points.

[0125] An embodiment of the present application further provides a chip storing a computer program, which implements the steps of the method for determining subway stray current leakage points when the computer program is executed by the chip.

[0126] The above describes specific embodiments of the present application. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recorded in the claims can be performed in an order different from that in the embodiments and still achieve the desired results. In addition, the processes depicted in the accompanying drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0127] The computer instructions include computer program codes, which may be in source code form, object code form, executable files or some intermediate forms, etc. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, USB flash drive, mobile hard disk, magnetic disk, optical disk, computer memory, read-only memory (ROM), random access memory (RAM), electric carrier signal, telecommunication signal and software distribution medium, etc. It should be noted that the content contained in the computer-readable medium may be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, computer-readable media do not include electric carrier signals and telecommunication signals.

[0128] It should be noted that, for the above-mentioned method embodiments, for the sake of simplicity of description, they are all expressed as a series of action combinations, but those skilled in the art should be aware that the present application is not limited by the described action sequence, because according to the present application, certain steps can be performed in other sequences or simultaneously. Secondly, those skilled in the art should also be aware that the embodiments described in the specification are all preferred embodiments, and the actions and modules involved are not necessarily required by the present application.

[0129] In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0130] The preferred embodiments of the present application disclosed above are only used to help explain the present application. The optional embodiments do not describe all the details in detail, nor do they limit the invention to the specific implementation methods described. Obviously, many modifications and changes can be made according to the content of the present application. The present application selects and specifically describes these embodiments in order to better explain the principles and practical applications of the present application, so that those skilled in the art can understand and use the present application well. The present application is only limited by the claims and their full scope and equivalents.

Claims

1. A method for determining subway stray current leakage points, characterized in that: include: S1, collecting train monitoring information, station monitoring information, station mileage information and line information under preset scenarios; S2, constructing a subway resistance network model according to the train monitoring information, the station monitoring information, the station mileage information and the line information; S3, constructing a group of sensitivity increment matrix equations about the change value of rail-to-ground transition resistance and the change value of rail potential through the subway resistance network model; S4, based on the sensitivity increment matrix equation group, determine the rail transition resistance distribution through an iterative calculation method, and determine the stray current leakage area according to the calculation result.

2. The method according to claim 1, characterized in that Said S1 comprises: When there is only one train running between two adjacent traction substations on the line, S11, collecting train monitoring information of the train at any time, including train position information and train traction current information; S12, collecting the station monitoring information including the power supply current information of the traction substation and the rail potential information of the station and the traction substation; S13, collecting the station mileage information including the station location information; S14, including the line information of the longitudinal resistance value of the rail per unit length and the transition resistance value of the rail per unit length to the ground.

3. The method according to claim 2, characterized in that The S2 comprises: S21, dividing the line between two adjacent traction substations into discrete equidistant divisions according to a preset division accuracy to obtain a plurality of nodes, wherein the nodes include traction substation nodes, station nodes and common nodes; S22, constructing an equivalent series resistance structure according to the longitudinal resistance value of the rail per unit length and the division accuracy; S23, constructing a parallel resistance structure according to the rail-to-ground transition resistance value per unit length and the division accuracy; S24, constructing a first equivalent current source according to the power supply current information of the traction substation; S25, constructing a second equivalent current source according to the train position information and the train traction current information; S26, constructing the subway resistance network model based on the equivalent series resistance structure, the parallel resistance structure, the first equivalent current source and the second equivalent current source.

4. The method according to claim 3, characterized in that The S3 includes: S31, when the train travels to any of the nodes, determining the rail potential information of the ideal station and the traction substation according to the station location information and the subway resistance network model; S32, according to the station location information, subtract the rail potential information of the ideal station and traction substation from the corresponding rail potential information of the station and traction substation to construct a rail potential difference column vector; S33, calculating a sensitivity matrix associating the rail transition resistance change value and the rail potential change value, wherein the rail potential change value is associated with the station node and the traction substation node; S34, repeatedly executing S31-S33 until the train travels to a plurality of different nodes, corresponding to a plurality of rail potential difference column vectors and a plurality of sensitivity matrices; S35, constructing the sensitivity increment matrix equation group according to the multiple rail potential difference column vectors and the multiple sensitivity matrices.

5. The method according to claim 4, characterized in that The S35 comprises: S351, arranging the plurality of rail potential difference column vectors according to the order in which the trains travel to the nodes to obtain an augmented rail potential difference column vector; S352, arranging the plurality of sensitivity matrices according to the order in which the trains travel to the nodes to obtain an augmented sensitivity matrix; S353, setting the rail-to-ground transition resistance increment column vector; S354, constructing the sensitivity increment matrix equation group according to the augmented rail potential difference column vector, the augmented sensitivity matrix and the rail-to-ground transition resistance increment column vector.

6. The method according to claim 5, characterized in that In S4, the rail transition resistance distribution is determined by an iterative calculation method based on the sensitivity increment matrix equation group, including: S41, setting an initial rail transition resistance value according to the rail-to-ground transition resistance value per unit length and the division accuracy; S42, calculating an initial rail potential initial value vector by a node voltage method based on the rail potential information of the station and the traction substation and the initial rail transition resistance value; S43, setting the initial value of the rail transition resistance increment to 0, and setting T=1, where T is a positive integer; S44, calculating the rail potential difference column vector and the sensitivity matrix of the T-th iteration according to the sensitivity increment matrix equation group, and constructing the sensitivity matrix equation group of the T-th iteration based on the calculation results; S45, based on the sensitivity matrix equation group of the T-th iteration, calculate the rail transition resistance increment of the T-th iteration by the least square method, and sum the calculation result with the rail transition resistance value of the T-1-th iteration to obtain the rail transition resistance value of the T-th iteration; S46, updating the rail potential difference column vector and the sensitivity matrix, T is incremented by 1, and S43 is executed until the rail transition resistance increment is less than a preset convergence accuracy, thereby obtaining the rail transition resistance distribution.

7. The method according to claim 6, characterized in that In S4, determining the stray current leakage area according to the calculation result includes: S47, determining a transition resistance threshold value according to the transition resistance value of the rail to the ground per unit length and the division accuracy; S48, comparing the rail transition resistance distribution with the transition resistance threshold, determining a transition resistance abnormal area, and using the transition resistance abnormal area as the stray current leakage area.

8. A device for determining subway stray current leakage points, characterized in that: include: A collection module is configured to collect train monitoring information, station monitoring information, station mileage information and line information under preset scenarios; A first construction module is configured to construct a subway resistance network model according to the train monitoring information, the station monitoring information, the station mileage information and the line information; The second construction module is configured to construct a sensitivity increment matrix equation group about the change value of the rail-to-ground transition resistance and the change value of the rail potential through the subway resistance network model; The determination module is configured to determine the rail transition resistance distribution by an iterative calculation method based on the sensitivity increment matrix equation group, and determine the stray current leakage area according to the calculation result.

9. A computing device, characterized in that include: Memory and processor; The memory is used to store computer executable instructions, and the processor is used to execute the computer executable instructions to implement the steps of the method for determining subway stray current leakage points as described in any one of claims 1 to 7.

10. A computer-readable storage medium storing computer instructions, characterized in that: When the instruction is executed by the processor, the steps of the method for determining the subway stray current leakage point described in any one of claims 1 to 7 are implemented.

Images