Method and device for determining subway stray current leakage points

Abstract

The present application provides a method and device for determining subway stray current leakage points, wherein the method comprises: collecting train monitoring information, station monitoring information, station mileage information, and line information under a preset scenario; constructing a subway resistance network model based on the train monitoring information, station monitoring information, station mileage information, and line information; constructing a sensitivity increment matrix equation group for changes in rail-to-ground transition resistance and rail potential using the subway resistance network model; determining the rail transition resistance distribution through an iterative calculation method based on the sensitivity increment matrix equation group, and determining the stray current leakage area based on the calculation results. This method enables the determination of subway stray current leakage areas through online monitoring data without the need for manual on-site measurements.

Description

Technical Field

[0001] The present application relates to the technical field of locating areas of insulation deterioration in subway rails, and more particularly to a method for determining subway stray current leakage points. 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 population growth in large cities, urban traffic congestion and air pollution are becoming increasingly prominent. To address this challenge, urban rail transit systems, with their environmentally friendly, low-energy consumption, high speeds, and large passenger capacity, have become the preferred choice and are experiencing rapid development. However, a significant issue is that due to on-site construction constraints, the rails cannot be completely insulated from the ground. Over the long term, rail transit systems experience varying degrees of insulation degradation due to the humid operating environment, aging insulation fasteners, and the accumulation of iron filings and dust on the rails. Consequently, the negative impact of stray currents is becoming increasingly apparent. When localized degradation of the rail insulation occurs, stray current leakage significantly increases. This current can penetrate these weak spots and accelerate corrosion of metal pipes and structural reinforcement surrounding the subway system. Therefore, while we enjoy the convenience of subways, we must address and take proactive measures to prevent electrochemical corrosion caused by stray currents and protect the safety of underground facilities.

[0003] Any localized fault in the insulation layer between the subway track and the ground accelerates the leakage of stray current. Once the leakage exceeds a threshold, drainage devices are required, but this can in turn increase the rail potential, creating an unfavorable cycle. Therefore, rather than relying solely on passive protection strategies to mitigate the hazards posed by stray currents, a proactive approach—precisely identifying and rapidly repairing damaged rail insulation—is a more effective strategy. This approach not only directly reduces stray current leakage but also mitigates the corrosion damage and DC bias caused by stray currents at their source, thereby more effectively ensuring the safety and long-term stability of the subway system.

[0004] There are two existing methods for locating insulation damage between track and ground. The first involves analyzing lines with localized insulation damage and extracting relevant characteristic values to locate the damaged point or area. However, these methods treat the transition resistance as a constant value, making it difficult to analyze and locate lines in other states. These methods primarily focus on lines with a single damaged area or point, and are unable to assess the location of multiple damage points. The second method involves measuring the transition resistance of subway lines in accordance with subway regulations. However, these methods are mostly offline, requiring workers to measure the relevant data on the line during train downtime. This is costly, inefficient, and fails to reflect the real-time status of track insulation. Additionally, some online measurement methods exist, but these require the installation of additional electrical equipment in the subway system, which is costly and difficult to calculate and assess the transition resistance status of the entire line, further hindering the ability to locate localized insulation damage between subway tracks and the ground. Summary of the Invention

[0005] In view of this, the present invention provides a method for determining subway stray current leakage points to address the technical deficiencies in the prior art. The present invention also provides an apparatus for determining subway stray current leakage points, 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, collects train monitoring information, station monitoring information, station mileage information and line information under preset scenarios;

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

[0009] S3, constructing a sensitivity increment matrix equation group regarding the change value of the rail-to-ground transition resistance and the change value of the 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, train monitoring information including train position information and train traction current information of the train at any time is collected;

[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, the line information including 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 partitions according to a preset partitioning 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 per unit length of the rail and the division accuracy;

[0019] S23, constructing a parallel resistance structure according to the rail-to-ground transition resistance 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: 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.

[0023] Optionally, the S3 includes:

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

[0025] S32, based on the station location information, subtracting the rail potential information of the ideal station and traction substation from the rail potential information of the corresponding 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 reaches 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 based on 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: Construct 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, calculating the rail transition resistance increment of the T-th iteration based on the sensitivity matrix equations of the T-th iteration by the least squares method, and summing the calculated 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 increments by 1, and executing S43 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 rail-to-ground transition resistance value per unit length and the division accuracy;

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

[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 building module is configured to build a subway resistance network model based on 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 through 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. When the instructions are executed by a processor, the steps of the method for determining subway stray current leakage points are implemented.

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

[0054] The method for determining subway stray current leakage points provided in this 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; and constructs a sensitivity increment matrix equation system for changes in rail-to-ground transition resistance and rail potential using the subway resistance network model. Based on the sensitivity increment matrix equation system, the rail transition resistance distribution is determined through an iterative calculation method, and the stray current leakage area is determined based on the calculation results. This method enables the determination of subway stray current leakage areas through online monitoring data without the need for manual on-site measurements. BRIEF DESCRIPTION OF THE DRAWINGS

[0055] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following briefly introduces the drawings required for use in the description of the embodiments. 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 creative work.

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

[0057] Figure 2 This 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 This is a schematic structural diagram of a device for determining subway stray current leakage points provided by an embodiment of the present application;

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

[0060] The following description sets forth many specific details to facilitate a thorough 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 scope of the present application. Therefore, the present application is not limited to the specific implementations disclosed below.

[0061] The terms used in one or more embodiments of the present application are for the purpose of describing specific embodiments only and are not intended to limit one or more embodiments of the present application. The singular forms "a", "the" 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 otherwise. 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," and the like 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 merely used to distinguish information of the same type from one another. For example, "first" may also be referred to as "second," and similarly, "second" may also be referred to as "first," without departing from the scope of one or more embodiments of the present application.

[0063] This application provides a method for determining subway stray current leakage points. This 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 in the following embodiments.

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

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

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

[0067] S3, constructing a sensitivity increment matrix equation group regarding the change value of the rail-to-ground transition resistance and the change value of the 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 is information including the train position and train electrical data during the train operation; station monitoring information is information including electrical data related to the station and traction substation; station mileage information is information including the location of the station on the train line; line information is information including information related to rail parameters on the train line.

[0070] Based on this, a bilateral rail-to-earth structural resistance network model (i.e., the subway resistance network model) was constructed using train monitoring information, station monitoring information, station mileage information, and line information. This model represents the equivalent circuit relationship formed by the current supplied by the train and traction substation, passing through the rails to the ground. After obtaining the subway resistance network model, the sensitivity increment matrix equations were determined by analyzing the correlation between the changes in the rail-to-ground transition resistance and the changes in the rail potential. Based on this sensitivity increment matrix equations, the rail transition resistance distribution was calculated using an iterative method. The calculated rail transition resistance distribution was analyzed to determine the areas corresponding to abnormal data and, therefore, the areas of stray current leakage.

[0071] Furthermore, 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 position 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 them with reference to the "CJJT 49-2020 Technical Standard for Corrosion Protection of Metro Stray Current". The specific setting results are not limited in this embodiment.

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

[0075] S21, divide the line between two adjacent traction substations into discrete and equidistant divisions according to a preset division accuracy to obtain multiple 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 the actual use scenario, if the total length between two adjacent traction substations is X, the number of equally spaced nodes N = n + 1, then the division accuracy L = X / n, the total number of resistance branches B = 2n + 1, the distance between the train and the traction substations at both ends is L1 and L2 respectively, after constructing the subway resistance network model, record the current source vector of the train at different node positions 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 = [D1D2…D k ], station rail potential information Where k is the total number of station nodes and traction substation nodes. At the same time, relevant line 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 marked with n+1 nodes, including traction substation nodes, station nodes and ordinary nodes; the rail is equivalent to a series resistance structure based on the longitudinal resistance of the rail combined with the accuracy, and the number of longitudinal resistors is n; then the rail-to-ground resistance is equivalent to a parallel resistance structure based on the rail transition resistance combined with the accuracy, and the number of longitudinal resistors is n+1; the traction substation power supply current 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 example. Figure 2A subway system model diagram of a method for determining subway stray current leakage points is shown in the following figure.

[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 potential of all nodes when the train is at different node positions

[0080] Furthermore, in step S3, a sensitivity increment matrix equation group for the rail-to-ground transition resistance change and the rail potential change is constructed using the subway resistance network model. In this embodiment, the specific implementation is as follows:

[0081] S31, when the train travels to 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, 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, repeating S31-S33 until the train travels to multiple different nodes, corresponding to multiple rail potential difference column vectors and multiple sensitivity matrices; S35, constructing the sensitivity increment matrix equation group according to multiple rail potential difference column vectors and multiple sensitivity matrices.

[0082] Furthermore, in step S35, the rail transition resistance distribution is determined by an iterative calculation method based on the sensitivity increment matrix equations. In this embodiment, the specific implementation method is as follows:

[0083] S351, according to the order in which the trains travel to the nodes, arrange the multiple rail potential difference column vectors to obtain the augmented rail potential difference column vector; S352, according to the order in which the trains travel to the nodes, arrange the multiple sensitivity matrices to obtain the augmented sensitivity matrix; S353, set the rail-to-ground transition resistance increment column vector; S354, according to the augmented rail potential difference column vector, the augmented sensitivity matrix and the rail-to-ground transition resistance increment column vector, construct the sensitivity increment matrix equation group.

[0084] Continuing with the above example, we use the parameter identification method 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. We make the total differential of the rail potential. Since the train traction current has nothing to do with the resistance value of each branch in the resistance 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 resistance of the rail transition resistance branch, so the sum of the change of resistance of any rail transition resistance branch multiplied by its sensitivity is U N The change in in 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 VU 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 station node i when the train is running at node v, where R is the transition resistance to be determined, and C i(c-n) is the (cn)th element in row i 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 admittance of the branch at the corresponding position in the circuit becomes larger and the 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 operating 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 node number of the train is 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, the process of determining the rail transition resistance distribution by an iterative calculation method based on the sensitivity increment matrix equation group 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 rail transition resistance increment initial value 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 squares 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, thereby obtaining 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 using a linear method. In order to solve this nonlinear problem, an iterative approach 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 and 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. T If ||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 aforementioned 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 based on the rail-to-ground transition resistance value 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: the rail-to-ground transition resistance of the operating line is R o =R s / L(R s =3Ω·km) is the threshold value. In addition, further limitation can be made based on the threshold value. For example, if it is less than 60% of the threshold value, it is identified as a transition resistance abnormal 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 shows a schematic diagram of 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 includes:

[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 based on 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 on the line running between two adjacent traction substations, the train monitoring information of the train at any time, including the train position information and the train traction current information, 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 multiple 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 ideal station and traction substation rail potential information is determined based on the station location information and the subway resistance network model; based on 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; 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; and the sensitivity increment matrix equation group is constructed based on 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, the plurality of 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, the plurality of 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] An initial rail transition resistance value is set based on the rail-to-ground transition resistance value per unit length and the division accuracy. An initial rail potential initial value vector is calculated using the node voltage method based on the rail potential information of the station and traction substation and the initial rail transition resistance value. The initial rail transition resistance increment is set to 0, and T is set to 1, where T is a positive integer. Based on the sensitivity increment matrix equation group, the rail potential difference column vector and the sensitivity matrix for the Tth iteration are calculated, and the sensitivity matrix equation group for the Tth iteration is constructed based on the calculation results. Based on the sensitivity matrix equation group for the Tth iteration, the rail transition resistance increment for the Tth iteration is calculated using the least squares method, and the calculation result is summed with the rail transition resistance value for the T-1th iteration to obtain the rail transition resistance value for 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, thereby obtaining 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 rail-to-ground transition resistance value per unit length 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 in this 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 this train monitoring information, station monitoring information, station mileage information, and line information; uses this subway resistance network model to construct a sensitivity increment matrix equation system for changes in rail-to-ground transition resistance and rail potential; and based on this sensitivity increment matrix equation system, determines the rail transition resistance distribution through an iterative calculation method, and determines the stray current leakage area based on the calculation results. This device enables the determination of subway stray current leakage areas through online monitoring data without the need for manual on-site measurements.

[0116] The above is a schematic scheme of a subway stray current leakage point determination device of this embodiment. It should be noted that the technical solution of the subway stray current leakage point determination device and the technical solution of the above-mentioned subway stray current leakage point determination method belong to the same concept. For details not described in detail in the technical solution of the subway stray current leakage point determination device, please refer to the description of the technical solution of the above-mentioned subway stray current leakage point determination method. In addition, the various components in the device embodiment should be understood as functional modules 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 4 shows a block diagram of a computing device 400 according to an embodiment of the present application. 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 a 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)) whether wired or wireless, such as an IEEE 802.11 wireless local area network (WLAN) wireless interface, a Worldwide 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. Figure 4 The computing device structure block diagram shown is for illustrative purposes only 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 can 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 computer, etc.), a mobile phone (e.g., a smartphone), a wearable computing device (e.g., a smartwatch, smart glasses, etc.), or other types of mobile devices, or a stationary computing device such as a desktop computer or PC. Computing device 400 can also be a mobile or stationary server.

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

[0122] The above is a schematic diagram of a computing device according to this embodiment. It should be noted that the technical solution of this computing device is based on the same concept as the technical solution of the aforementioned method for determining subway stray current leakage points. For details not described in detail in the technical solution of the computing device, please refer to the description of the technical solution of the aforementioned 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 subway stray current leakage points.

[0124] The above is a schematic diagram of a computer-readable storage medium according to this embodiment. It should be noted that the technical solution of this storage medium is based on the same concept as the technical solution of the aforementioned method for determining subway stray current leakage points. For details not described in detail in the technical solution of the storage medium, please refer to the description of the technical solution of the aforementioned 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 executed by the chip.

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

[0127] The computer instructions include computer program code, which may be in source code form, object code form, executable file, or some intermediate form. 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 aforementioned method embodiments, for ease of description, they are all expressed as a series of action combinations, but those skilled in the art should be aware that this application is not limited by the order of the actions described, because according to this application, certain steps can be performed in other orders or simultaneously. Secondly, those skilled in the art should also be aware that the embodiments described in this specification are all preferred embodiments, and the actions and modules involved are not necessarily required by this application.

[0129] In the above embodiments, the description of each embodiment has its own focus. 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 intended only to help illustrate the present application. The optional embodiments do not describe all details in detail, nor do they limit the invention to the specific embodiments described. Obviously, many modifications and variations can be made based on the content of this application. This application selects and describes these embodiments in detail in order to better explain the principles and practical applications of this application, so that those skilled in the art can better understand and utilize this application. This application is limited only 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, collects train monitoring information, station monitoring information, station mileage information and line information under preset scenarios; S2, constructing a subway resistance network model based on 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 rail potential change value through the subway resistance network model. Specifically, S31, when the train travels to any node of the subway resistance network model, determining the ideal station and traction substation rail potential information according to the station location information and the subway resistance network model; S32, according to the station location information, subtracting the ideal station and traction substation rail potential information from the corresponding station and traction substation rail potential information to construct a rail potential difference column vector; S33, calculating a sensitivity matrix associated with the change value of rail-to-ground transition resistance and the change value of rail potential, wherein , the rail potential change value is associated with the station node and the traction substation node; S34, repeating S31-S33 to obtain multiple rail potential difference column vectors and multiple sensitivity matrices corresponding to when the train travels to multiple different nodes; S35, according to the order in which the train travels to the nodes, arranging the multiple rail potential difference column vectors to obtain an augmented rail potential difference column vector, arranging the multiple sensitivity matrices to obtain an augmented sensitivity matrix, setting the rail-to-ground transition resistance increment column vector, and 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; 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, train monitoring information including train position information and train traction current information of the train at any time is collected; 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, the line information including 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 Said S2 comprises: S21, dividing the line between two adjacent traction substations into discrete equidistant partitions according to a preset partitioning 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 per unit length of the rail and the division accuracy; S23, constructing a parallel resistance structure according to the rail-to-ground transition resistance 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: 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.

4. The method according to claim 1, wherein 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 increment matrix equation group of the T-th iteration based on the calculation results; S45, calculating the rail transition resistance increment of the T-th iteration based on the sensitivity increment matrix equation group of the T-th iteration by the least squares method, and summing 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 increments by 1, and executing S43 until the rail transition resistance increment is less than a preset convergence accuracy, thereby obtaining the rail transition resistance distribution.

5. The method according to claim 4, 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 rail-to-ground transition resistance value per unit length and the division accuracy; S48 , comparing the rail transition resistance distribution with the transition resistance threshold, determining an abnormal transition resistance region, and using the abnormal transition resistance region as the stray current leakage region.

6. 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 building module is configured to build a subway resistance network model based on 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. Specifically, when the train travels to any node of the subway resistance network model, the ideal station and traction substation rail potential information is determined according to the station location information and the subway resistance network model; 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; and the sensitivity matrix associated with the change value of the rail-to-ground transition resistance and the change value of the rail potential is calculated. In the embodiment, the rail potential change value is associated with the station node and the traction substation node; the second construction module is repeatedly called to obtain a plurality of rail potential difference column vectors and a plurality of sensitivity matrices corresponding to the train when the train travels to a plurality of different nodes; the plurality of rail potential difference column vectors are arranged according to the order in which the train travels to the nodes to obtain an augmented rail potential difference column vector, the plurality of sensitivity matrices are arranged to obtain an augmented sensitivity matrix, a rail-to-ground transition resistance increment column vector is set, and the sensitivity increment matrix equation group is constructed according to the augmented rail potential difference column vector, the augmented sensitivity matrix, and the rail-to-ground transition resistance increment column vector; The determination module is configured to determine the rail transition resistance distribution through an iterative calculation method based on the sensitivity increment matrix equation group, and determine the stray current leakage area according to the calculation result.

7. 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 5.

8. 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 subway stray current leakage points described in any one of claims 1 to 5 are implemented.

Images