A method and system for locating a stray current leakage position in rail transit

Abstract

A method and system for locating the position of stray current leakage in rail transit, the method comprising: building a basic simulation model of rail transit based on CDEGS software; according to the basic simulation model, simulating and calculating the stray current of the substation when the train is at different positions according to a first preset rule; analyzing the stray current of the substation according to a second preset rule, obtaining the weak insulation position of rail transit, and locating the stray current leakage position. The present invention analyzes the changes in the stray current of the substation and obtains the weak insulation position of rail transit by integrating multi-source information such as rail transit and power grid substations, thereby solving the problem of difficulty in locating the stray current leakage position and providing data basis for targeted management of stray current in rail transit.

Description

Technical Field

[0001] The present invention relates to the technical field of rail transit stray current, and in particular to a method and system for locating a rail transit stray current leakage position. Background Art

[0002] The single-pole operation of ultra-high and extra-high voltage DC projects, urban rail transit operation, photovoltaic, wind power, DC distribution networks, etc. will all generate DC components. The DC components flow in urban pipe corridors, energy pipe corridors and power grids, and flow into transformers, causing DC magnetic bias in transformers and transformer vibration, which not only shortens the life of transformers and endangers the safety of large power grid operations, but also causes noise disturbances, leading to complaints from nearby residents and affecting the social image of power grid companies.

[0003] The simulation model of stray current distribution of urban rail transit based on CDEGS software has been recognized and widely used by experts and scholars because it can well equate the metal structure of the line and the line operation environment. The prior art discloses a device and method for locating insulation damage of rail transit running rails, including a first detection device and a second detection device separated by a preset distance, and there is a roughly located insulation damage point between the first detection device and the second detection device, and the insulation damage point is located on the first running rail or the second running rail of the rail transit. The method includes a method for locating insulation damage of the running rail when there is no rail weld on the rails corresponding to the first detection device and the second detection device, and a method for locating insulation damage of the running rail when there is a rail weld on the rail corresponding to the first detection device or the second detection device. However, this method needs to rely on detection devices and infrared thermal imagers for detection, and can only be carried out during the subway suspension period, which affects the operation of the subway. The prior art also proposes a method for identifying high-risk areas of stray current leakage in subways based on cloud computing. The whole-line monitoring data is obtained through the stray current corrosion monitoring system, and an intelligent identification model for high-risk areas of leakage at multiple monitoring locations is established using the whole-line monitoring data. Based on the parallel cloud computing results of the intelligent identification model, the high-risk areas of stray current leakage on the whole line can be intelligently and accurately located. However, this method requires that the subway track is installed with a corrosion monitoring system, but not all rail transit is installed with this monitoring system, and this method has measurement limitations. The prior art also discloses a method for locating underground stray current sources. The propagation speed of electromagnetic waves in the earth is measured based on the relative dielectric constant of the earth. An underground stray current positioning model of five probes is established to locate the underground impact source. The array composed of five probes is a cross array. xy is defined as a plane. The central observation point A1 (0, 0, 0) of the cross array is used as the origin of the coordinates. The coordinates of the impact source S of the stray current impact are (x, y, z); the four observation points are symmetrically distributed according to the central observation point A1, and the distance from the four observation points to the central observation point A1 is a as the array element distance a. The central observation point and the four observation points are respectively provided with probes to measure the potential time domain waveform; the time delay of the impact source to each observation point is determined based on the potential time domain waveform, and the cross array method is used to locate the coordinates of the impact source S. However, this method does not involve the fusion analysis of rail transit and power grid data, and there are certain difficulties in operation, so it is not used for promotion. Therefore, there is an urgent need for a method and system for locating the position of stray current leakage in rail transit, which can locate the weak insulation position of rail transit by integrating multi-source data such as rail transit and power grid and analyzing the changes in stray current in substations. Summary of the invention

[0004] In view of the above problems, the present invention is proposed to provide a method to overcome the above problems or at least partially solve the above problems.

[0005] In order to solve the above technical problems, the embodiments of the present application disclose the following technical solutions:

[0006] In a first aspect, an embodiment of the present invention discloses a method for locating a rail transit stray current leakage position, comprising:

[0007] S100. Build a basic rail transit simulation model based on CDEGS software;

[0008] S200. According to the basic simulation model, the stray current of the substation when the train is at different positions is simulated and calculated according to the first preset rule;

[0009] S300. Analyze the stray current of the substation according to the second preset rule, obtain the weak position of rail transit insulation, and locate the stray current leakage position.

[0010] Furthermore, in S100, a basic simulation model of rail transit is built based on the MALZ module of CDEGS software, and the basic simulation model at least includes a conductor model, a substation grounding grid model, a transformer model, a contact network model, a rail model, a drainage network model, a traction station model and a train model.

[0011] Furthermore, the basic simulation model sets a conductor at a preset distance, and the conductor is used to represent a train. When the train moves to the position of the conductor, the conductor at the position is set to be excited, and the conductors at other positions are not set to be excited. In the traction station model and the train model, the traction station shares a busbar, each train has a busbar, and the excitation of the traction station is equal to the sum of the excitations of all trains, but in opposite directions.

[0012] Further, in S200, the stray current of the substation when the train is at different positions is simulated and calculated according to the first preset rule, and the first preset rule at least includes: obtaining and simulating the stray current corresponding to the substation at each position of the train when the rail transit insulation is normal and the insulation is weak;

[0013] When the simulation calculation of rail transit insulation is normal, the train excitation is always set to 1000A. When the train moves along the track, the stray current corresponding to the substation at each position of the train is obtained and recorded as the reference sequence {S1, S2, …, S n},S i is the stray current of the substation when the train is at position i, 1≤i≤n;

[0014] When the simulation calculates that the insulation of the rail transit track is weak, the train excitation is always 1000A. When the train moves along the track, the stray current corresponding to the substation at each position of the train is recorded as the actual sequence {B1, B2, …, B n}, B iis the stray current of the substation when the train is at position i, 1≤i≤n.

[0015] Further, in S300, the stray current of the substation is analyzed according to a second preset rule to obtain the weak position of rail transit insulation and locate the stray current leakage position, and the second preset rule at least includes:

[0016] S301. Calculate the reference sequence {S1, S2, ..., S n} and the actual sequence {B1,B2,…,B n}, we get {CS2,CS3,…,CS n}, {CB2, CB3, …, CB n}, where CS j =S j -S j-1 , C.B. j =B j -B j-1 , 2≤j≤n;

[0017] S302. Compare the stray current magnitudes of the same positions of the reference sequence and the actual sequence. If CS j >CB j , then let the Boolean value Z j =1, otherwise let the Boolean value Z j =0,2≤j≤n;

[0018] S303. Determine the calculation interval {Q1, Q2, ..., Qj} of the correlation coefficient according to the Boolean value Zj. m};

[0019] S304. Calculate the correlation coefficient between the reference sequence and the actual sequence in each interval segment one by one;

[0020] S305. Determine whether the correlation coefficient between the reference sequence and the actual sequence corresponding to the interval segment is greater than a preset threshold. If so, return to S304 to calculate the correlation coefficient of the next interval; if less than the preset threshold, proceed to S306;

[0021] S306. Calculate the absolute value of the interval difference, and update the insulation weak position according to the interval corresponding to the absolute value of the interval difference.

[0022] Further, in S303, the calculation interval {Q1, Q2, ..., Q m The specific method of} includes: for any j, if Z j =1 and Z j+1 =0, or satisfies Z j = 0 and Z j+1=1, then the interval [j-3,j+4] is used as an interval for correlation coefficient calculation. There are 8 stray current data in this calculation interval. By traversing all j, the calculation interval of correlation coefficient {Q1,Q2,…,Q m}, where Q k represents the kth calculation interval, Q k =[q k,1 ,q k,2 ,…,q k,8 ], here q k,1 is the train position corresponding to j-3, q k,8 is the train position corresponding to j+4, 1≤k≤m.

[0023] Further, in S304, the correlation coefficients of the reference sequence and the actual sequence in the interval segment are calculated one by one. The specific method includes: for the calculation interval Q k , using the Pearson correlation coefficient to calculate the interval Q k The corresponding reference sequence and the actual sequence The correlation coefficient value is recorded as P k .

[0024] Further, in S306, the absolute value of the interval difference is calculated, and the insulation weak position is updated according to the interval corresponding to the absolute value of the interval difference; the specific method includes:

[0025] Calculate the absolute value of the interval difference D k , The absolute value of the interval difference D k Compared with the difference absolute value initialization value D0, the difference absolute value initialization value D0 is 0; if D k >D0, then D k Assign the value of D0 to D0, update the value of D0, and k The corresponding interval Q k =[q k,1 ,q k,2 ,…,q k,8 ] to update the insulation weakness position A = q k,4 ; Determine whether the calculation of all intervals is completed. If not, return to S304. If the calculation of all intervals is completed, output the insulation weak position.

[0026] In a second aspect, an embodiment of the present invention discloses a system for locating a rail transit stray current leakage position, comprising: a rail transit basic simulation model building unit, a substation stray current calculation unit and a stray current leakage position locating unit; wherein:

[0027] Rail transit basic simulation model building unit, used to build rail transit basic simulation model based on CDEGS software;

[0028] A substation stray current calculation unit, used to simulate and calculate the stray current of the substation when the train is at different positions according to the basic simulation model and a first preset rule;

[0029] The stray current leakage position positioning unit is used to analyze the stray current of the substation according to the second preset rule, obtain the weak position of the rail transit insulation, and locate the stray current leakage position.

[0030] In a third aspect, an embodiment of the present invention discloses an electronic device, including:

[0031] one or more processors;

[0032] A memory for storing one or more programs;

[0033] When the one or more programs are executed by the one or more processors, the one or more processors implement the method for locating the leakage position of molten iron stray current.

[0034] The beneficial effects of the above technical solution provided by the embodiment of the present invention include at least:

[0035] The present invention discloses a method and system for locating the position of stray current leakage in rail transit, the method comprising: building a basic simulation model of rail transit based on CDEGS software; according to the basic simulation model, simulating and calculating the stray current of the substation when the train is at different positions according to a first preset rule; analyzing the stray current of the substation according to a second preset rule, obtaining the weak insulation position of rail transit, and locating the stray current leakage position. The present invention can solve the problem of difficulty in locating the stray current leakage position by integrating multi-source information such as rail transit and power grid substations, analyzing the changes in the stray current of the substation, and obtaining the weak insulation position of rail transit, thereby providing a data basis for the targeted management of stray current in rail transit.

[0036] The technical solution of the present invention is further described in detail below through the accompanying drawings and embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] The accompanying drawings are used to provide a further understanding of the present invention and constitute a part of the specification. Together with the embodiments of the present invention, they are used to explain the present invention and do not constitute a limitation of the present invention. In the accompanying drawings:

[0038] Figure 1 This is a flow chart of a method for locating a stray current leakage position in rail transit in Embodiment 1 of the present invention;

[0039] Figure 2 This is a flow chart of calculating the insulation weak position in Embodiment 1 of the present invention;

[0040] Figure 3 This is a schematic diagram of building a rail transit basic simulation model in Example 1 of the present invention;

[0041] Figure 4 This is a schematic diagram of the structure of an electronic device in Example 3 of the present invention. DETAILED DESCRIPTION

[0042] The exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although the exemplary embodiments of the present disclosure are shown in the accompanying drawings, it should be understood that the present disclosure can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided to enable a more thorough understanding of the present disclosure and to fully convey the scope of the present disclosure to those skilled in the art.

[0043] In order to solve the problems existing in the prior art, the embodiments of the present invention provide.

[0044] Example 1

[0045] The present invention discloses a method for locating a stray current leakage position in rail transit. Figure 1 ,include:

[0046] S100. Build a basic simulation model of rail transit based on CDEGS software; in S100 of this embodiment, build a basic simulation model of rail transit based on the MALZ module of CDEGS software, and the basic simulation model at least includes a conductor model, a substation grounding network model, a transformer model, a contact network model, a rail model, a drainage network model, a traction station model and a train model. The basic simulation model sets a conductor at a preset distance, and the conductor is used to represent the train. When the train moves to the position of the conductor, the conductor at the position is set to excitation, and the conductors at other positions are not set to excitation; in the traction station model and the train model, the traction station shares a busbar, and each train has a busbar, and the excitation of the traction station is equal to the sum of the excitations of all trains, but in the opposite direction.

[0047] Specifically, CDEGS is the abbreviation of the first letters of the English initials of Current Distribution, Electromagnetic Fields, Grounding and Soil Structure Analysis. It is a powerful tool software for solving engineering problems such as power system grounding, electromagnetic fields and electromagnetic interference, and can solve problems such as cathodic protection. The software contains multiple engineering application modules, among which the MALZ module is a tool specifically used to analyze the current distribution, electromagnetic fields and grounding systems of underground and above-ground power systems. The MALZ module of CDEGS software is used to build a basic simulation model of rail transit.

[0048] S200. According to the basic simulation model, the stray current of the substation when the train is at different positions is simulated and calculated according to the first preset rule; in S200 of this embodiment, the stray current of the substation when the train is at different positions is simulated and calculated according to the first preset rule, and the first preset rule at least includes: respectively obtaining and simulating the stray current corresponding to the substation at each position of the train when the insulation of the rail transit is normal and the insulation is weak;

[0049] When the simulation calculation of rail transit insulation is normal, the train excitation is always set to 1000A. When the train moves along the track, the stray current corresponding to the substation at each position of the train is obtained and recorded as the reference sequence {S1, S2, …, S n},S i is the stray current of the substation when the train is at position i, 1≤i≤n;

[0050] When the simulation calculates that the insulation of the rail transit track is weak, the train excitation is always 1000A. When the train moves along the track, the stray current corresponding to the substation at each position of the train is recorded as the actual sequence {B1, B2, …, B n}, B i is the stray current of the substation when the train is at position i, 1≤i≤n.

[0051] S300. Analyze the stray current of the substation according to the second preset rule, obtain the weak position of rail transit insulation, and locate the stray current leakage position.

[0052] In S300 of this embodiment, the stray current of the substation is analyzed according to the second preset rule to obtain the weak position of rail transit insulation and locate the stray current leakage position, such as Figure 2 , the second preset rule at least includes:

[0053] S301. Calculate the reference sequence {S1, S2, ..., S n} and the actual sequence {B1,B2,…,Bn}, we get {CS2,CS3,…,CS n}, {CB2, CB3, …, CB n}, where CS j =S j -S j-1 , C.B. j =B j -B j-1 , 2≤j≤n;

[0054] S302. Compare the stray current magnitudes of the same positions of the reference sequence and the actual sequence. If CS j >CB j , then let the Boolean value Z j =1, otherwise let the Boolean value Z j =0, 2≤j≤n;

[0055] S303. According to the Boolean value Z j Determine the calculation interval of the correlation coefficient {Q1,Q2,…,Q m};

[0056] S304. Calculate the correlation coefficient between the reference sequence and the actual sequence in each interval segment one by one;

[0057] S305. Determine whether the correlation coefficient between the reference sequence and the actual sequence corresponding to the interval segment is greater than a preset threshold. If so, return to S304 to calculate the correlation coefficient of the next interval; if less than the preset threshold, proceed to S306; in this embodiment, the preset threshold is preferably 0.98.

[0058] S306. Calculate the absolute value of the interval difference, and update the insulation weak position according to the interval corresponding to the absolute value of the interval difference.

[0059] In S303 of some preferred embodiments, the calculation interval {Q1, Q2, ..., Q m The specific method of} includes: for any j, if Z j =1 and Z j+1 =0, or satisfies Z j = 0 and Z j+1 =1, then the interval [j-3,j+4] is used as an interval for correlation coefficient calculation. There are 8 stray current data in this calculation interval. By traversing all j, the calculation interval of correlation coefficient {Q1,Q2,…,Q m}, where Q k represents the kth calculation interval, Q k =[q k,1 ,q k,2 ,…,q k,8], here q k,1 is the train position corresponding to j-3, q k,8 is the train position corresponding to j+4, 1≤k≤m.

[0060] In S304 of some preferred embodiments, the correlation coefficients of the reference sequence and the actual sequence in the interval segments are calculated one by one. The specific method includes: for calculating the interval Q k , using the Pearson correlation coefficient to calculate the interval Q k The corresponding reference sequence and the actual sequence The correlation coefficient value is recorded as P k .

[0061] 9. In some preferred embodiments, in S306, the absolute value of the interval difference is calculated, and the insulation weak position is updated according to the interval corresponding to the absolute value of the interval difference; the specific method includes:

[0062] Calculate the absolute value of the interval difference D k , The absolute value of the interval difference D k Compared with the difference absolute value initialization value D0, the difference absolute value initialization value D0 is 0; if D k >D0, then D k Assign the value of D0 to D0, update the value of D0, and k The corresponding interval Q k =[q k,1 ,q k,2 ,…,q k,8 ] to update the insulation weakness position A = q k,4 ; Determine whether the calculation of all intervals is completed. If not, return to S304. If the calculation of all intervals is completed, output the insulation weak position.

[0063] In order to better understand the present embodiment, the technical solution and effects of the present invention are described in detail with a specific embodiment below:

[0064] Build a basic simulation model of multi-source DC components, such as Figure 2 As shown, the thick line segment represents rail transit, and the thin line segment represents the power grid topology;

[0065] According to the basic simulation model, the stray current of the substation when the train is at different positions is simulated and calculated, as shown in Table 1 (only some examples)

[0066] Table 1

[0067]

[0068] According to the substation stray current obtained in the first step, the substation stray current is analyzed, the weak position of rail transit insulation is calculated, the first-order difference of the reference sequence and the actual sequence is calculated, and the size relationship is compared, as shown in Table 2:

[0069] Table 2

[0070]

[0071] The interval correlation coefficient and interval difference insulation value are calculated, as shown in Table 3. When the interval correlation coefficient is less than 0.98, the position with the largest interval difference insulation value is 162. Therefore, it can be judged that the position of stray current leakage is 162, which is consistent with the actual situation.

[0072] Table 3

[0073]

[0074]

[0075] The present embodiment discloses a method for locating the position of stray current leakage in rail transit, including: building a basic simulation model of rail transit based on CDEGS software; according to the basic simulation model, simulating and calculating the stray current of the substation when the train is at different positions according to a first preset rule; analyzing the stray current of the substation according to a second preset rule, obtaining the weak insulation position of rail transit, and locating the stray current leakage position. The present invention analyzes the changes in the stray current of the substation and obtains the weak insulation position of rail transit by integrating multi-source information such as rail transit and power grid substations, thereby solving the problem of difficulty in locating the stray current leakage position and providing data basis for targeted management of stray current in rail transit.

[0076] Example 2

[0077] Based on the same inventive concept, this embodiment discloses a system for locating the position of stray current leakage in rail transit, including: a rail transit basic simulation model building unit, a substation stray current calculation unit and a stray current leakage position locating unit; wherein:

[0078] Rail transit basic simulation model building unit, used to build rail transit basic simulation model based on CDEGS software;

[0079] A substation stray current calculation unit, used to simulate and calculate the stray current of the substation when the train is at different positions according to the basic simulation model and a first preset rule;

[0080] The stray current leakage position positioning unit is used to analyze the stray current of the substation according to the second preset rule, obtain the weak position of rail transit insulation, and locate the stray current leakage position.

[0081] Among them, the specific working methods of the rail transit basic simulation model building unit, the substation stray current calculation unit and the stray current leakage position positioning unit have been described in detail in Example 1, and will not be repeated in this embodiment.

[0082] Example 3

[0083] Based on the same inventive concept, an embodiment of the present disclosure also provides an electronic device. Figure 4 FIG. 1 is a schematic diagram of the structure of an electronic device according to an embodiment of the present disclosure. Figure 4 As shown, an embodiment of the present disclosure provides an electronic device including: one or more processors 101, a memory 102, and one or more I / O interfaces 103. The memory 102 stores one or more programs, and when the one or more programs are executed by the one or more processors, the one or more processors implement any optimization method in the above embodiments; the one or more I / O interfaces 103 are connected between the processor and the memory, and are configured to implement information interaction between the processor and the memory.

[0084] Among them, the processor 101 is a device with data processing capabilities, including but not limited to a central processing unit (CPU), etc.; the memory 102 is a device with data storage capabilities, including but not limited to random access memory (RAM, more specifically SDRAM, DDR, etc.), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory (FLASH); the I / O interface (read-write interface) 103 is connected between the processor 101 and the memory 102, and can realize information interaction between the processor 101 and the memory 102, including but not limited to a data bus (Bus), etc.

[0085] In some embodiments, the processor 101 , the memory 102 , and the I / O interface 103 are connected to each other via a bus 104 , and further connected to other components of the computing device.

[0086] In some embodiments, the one or more processors 101 include a field programmable gate array.

[0087] According to an embodiment of the present disclosure, a computer-readable medium is further provided, wherein a computer program is stored on the computer-readable medium, wherein when the program is executed by a processor, the steps in any optimization method in the above-mentioned embodiment are implemented.

[0088] It should be understood that the specific order or hierarchy of steps in the disclosed process is an example of an exemplary method. Based on design preferences, it should be understood that the specific order or hierarchy of steps in the process can be rearranged without departing from the scope of protection of the present disclosure. The attached method claims present the elements of the various steps in an exemplary order and are not intended to be limited to the specific order or hierarchy described.

[0089] In the above detailed description, various features are grouped together in a single embodiment to simplify the disclosure. This method of disclosure should not be interpreted as reflecting an intention that the embodiments of the claimed subject matter require more features than are clearly stated in each claim. On the contrary, as reflected in the appended claims, the invention is in a state of having less than all the features of the disclosed individual embodiments. Therefore, the appended claims are hereby expressly incorporated into the detailed description, with each claim standing on its own as a separate preferred embodiment of the invention.

[0090] Those skilled in the art will also appreciate that the various illustrative logic blocks, modules, circuits, and algorithmic steps described in conjunction with the embodiments herein can all be implemented as electronic hardware, computer software, or a combination thereof. In order to clearly illustrate the interchangeability between hardware and software, various illustrative components, blocks, modules, circuits, and steps are generally described above around their functions. Whether such functions are implemented as hardware or software depends on specific applications and the design constraints imposed on the entire system. A skilled person can implement the described functions in an alternative manner for each specific application, but such implementation decisions should not be interpreted as departing from the scope of protection of the present disclosure.

[0091] The steps of the method or algorithm described in conjunction with the embodiments herein may be directly embodied as hardware, a software module executed by a processor, or a combination thereof. The software module may be located in a RAM memory, a flash memory, a ROM memory, an EPROM memory, an EEPROM memory, a register, a hard disk, a mobile disk, a CD-ROM, or any other form of storage medium well known in the art. An exemplary storage medium is connected to the processor so that the processor can read information from the storage medium and can write information to the storage medium. Of course, the storage medium may also be an integral part of the processor. The processor and the storage medium may be located in an ASIC. The ASIC may be located in a user terminal. Of course, the processor and the storage medium may also be present in a user terminal as discrete components.

[0092] For software implementation, the techniques described in this application can be implemented with modules (e.g., procedures, functions, etc.) that perform the functions described in this application. These software codes can be stored in a memory unit and executed by a processor. The memory unit can be implemented within the processor or outside the processor. In the latter case, it is coupled to the processor in a communication manner via various means, which are well known in the art.

[0093] The above description includes examples of one or more embodiments. Of course, it is impossible to describe all possible combinations of components or methods for the purpose of describing the above embodiments, but it should be recognized by those skilled in the art that the various embodiments may be further combined and arranged. Therefore, the embodiments described herein are intended to cover all such changes, modifications and variations that fall within the scope of protection of the appended claims. In addition, with respect to the term "comprising" used in the specification or claims, the word is covered in a manner similar to the term "including", just as "including," is explained as a transitional word in the claims. In addition, any term "or" used in the specification of the claims is intended to mean "non-exclusive or".

Claims

1. A method for locating the position of stray current leakage in rail transit, characterized in that: include: S100. Build a basic simulation model of rail transit based on CDEGS software; In S100, build a basic simulation model of rail transit based on the MALZ module of CDEGS software, the basic simulation model of rail transit includes at least a conductor model, a substation grounding network model, a transformer model, a contact network model, a rail model, a drainage network model, a traction station model and a train model; the basic simulation model of rail transit sets a conductor at a preset distance, the conductor is used to represent the train, when the train moves to the position of the conductor, the conductor at the position is set to excitation, and the conductors at other positions are not set to excitation; in the traction station model and the train model, the traction station shares a busbar, each train has a busbar, and the excitation of the traction station is equal to the sum of the excitations of all trains, but in opposite directions; S200. According to the rail transit basic simulation model, the stray current of the substation when the train is at different positions is simulated and calculated according to the first preset rule; The first preset rule at least includes: obtaining and calculating the stray current corresponding to the substation at each position of the train when the rail transit insulation is normal and when the insulation is weak; When the simulation calculation of rail transit insulation is normal, the train excitation is always set to 1000A. When the train moves along the track, the stray current corresponding to the substation at each position of the train is obtained and recorded as the reference sequence {S1, S2, …, S n },S i is the stray current of the substation when the train is at position i, 1≤i≤n; When the simulation calculates that the insulation of the rail transit track is weak, the train excitation is always 1000A. When the train moves along the track, the stray current corresponding to the substation at each position of the train is recorded as the actual sequence {B1, B2, …, B n }, B i is the stray current of the substation when the train is at position i, 1≤i≤n; S300. Analyze the stray current of the substation according to the second preset rule, obtain the weak position of rail transit insulation, and locate the stray current leakage position; the second preset rule includes at least: S301. Calculate the reference sequence {S1, S2, ..., S n } and the actual sequence {B1,B2,…,B n }, we get {CS2,CS3,…,CS n }, {CB2, CB3, …, CB n }, where CS j =S j -S j-1 , C.B. j =B j -B j-1 , 2≤j≤n; S302. Compare the stray current magnitudes of the same positions of the reference sequence and the actual sequence. If CS j >CB j , then let the Boolean value Z j =1, otherwise let the Boolean value Z j =0,2≤j≤n; S303. According to the Boolean value Z j Determine the calculation interval of the correlation coefficient {Q1,Q2,…,Q m }; S304. Calculate the correlation coefficient between the reference sequence and the actual sequence in each interval segment one by one; S305. Determine whether the correlation coefficient between the reference sequence and the actual sequence corresponding to the interval segment is greater than a preset threshold. If so, return to S304 to calculate the correlation coefficient of the next interval; if less than the preset threshold, proceed to S306; S306. Calculate the absolute value of the interval difference, and update the insulation weak position according to the interval corresponding to the absolute value of the interval difference.

2. A method for locating a rail transit stray current leakage position as claimed in claim 1, characterized in that: In S303, according to the Boolean value Z j Determine the calculation interval of the correlation coefficient {Q1,Q2,…,Q m The specific method of} includes: for any j, if Z j =1 and Z j+1 =0, or satisfies Z j = 0 and Z j+1 =1, then the interval [j-3,j+4] is used as an interval for correlation coefficient calculation. There are 8 stray current data in this calculation interval. By traversing all j, the calculation interval of correlation coefficient {Q1,Q2,…,Q m }, where Q k represents the kth calculation interval, Q k =[q k,1 ,q k,2 ,…,q k,8 ], here q k,1 is the train position corresponding to j-3, q k,8 is the train position corresponding to j+4, 1≤k≤m.

3. A method for locating a rail transit stray current leakage position as claimed in claim 1, characterized in that: In S304, the correlation coefficients between the reference sequence and the actual sequence in each interval are calculated one by one. The specific method includes: for the calculation interval Q k , using the Pearson correlation coefficient to calculate the interval Q k The corresponding reference sequence and the actual sequence The correlation coefficient value is recorded as P k .

4. A method for locating a rail transit stray current leakage position as claimed in claim 1, characterized in that: In S306, the absolute value of the interval difference is calculated, and the insulation weak position is updated according to the interval corresponding to the absolute value of the interval difference. The specific method includes: Calculate the absolute value of the interval difference D k , The absolute value of the interval difference D k Compared with the difference absolute value initialization value D0, the difference absolute value initialization value D0 is 0; if D k >D0, then D k Assign the value of D0 to D0, update the value of D0, and k The corresponding interval Q k =[q k,1 ,q k,2 ,…,q k,8 ] to update the insulation weakness position A = q k,4 ; Determine whether the calculation of all intervals is completed. If not, return to S304. If the calculation of all intervals is completed, output the insulation weak position.

5. A system for locating the position of stray current leakage in rail transit, using the method for locating the position of stray current leakage in rail transit as claimed in any one of claims 1 to 4, characterized in that: include: Rail transit basic simulation model building unit, substation stray current calculation unit and stray current leakage location positioning unit; Among them: Rail transit basic simulation model building unit, used to build rail transit basic simulation model based on CDEGS software; A substation stray current calculation unit, used to simulate and calculate the stray current of the substation when the train is at different positions according to the rail transit basic simulation model and a first preset rule; The stray current leakage position positioning unit is used to analyze the stray current of the substation according to a second preset rule, obtain the weak position of rail transit insulation, and locate the stray current leakage position.

6. An electronic device comprising: one or more processors; A memory for storing one or more programs; When the one or more programs are executed by the one or more processors, the one or more processors implement the method for locating the position of stray current leakage in rail transit as described in any one of claims 1 to 4.

Images