CN112836364A - Unified calculation method for urban rail transit stray current - Google Patents
Unified calculation method for urban rail transit stray current Download PDFInfo
- Publication number
- CN112836364A CN112836364A CN202110107753.2A CN202110107753A CN112836364A CN 112836364 A CN112836364 A CN 112836364A CN 202110107753 A CN202110107753 A CN 202110107753A CN 112836364 A CN112836364 A CN 112836364A
- Authority
- CN
- China
- Prior art keywords
- current
- rail transit
- urban rail
- voltage
- calculation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
Abstract
The invention discloses a unified calculation method for urban rail transit stray current, which comprises the following steps: step 1: monitoring and acquiring a train traction working condition, and acquiring related parameters of an urban rail transit reflux system; step 2: constructing a unified calculation model of the stray current of the urban rail transit system; and step 3: and constructing a function expression of the voltage and the current based on the model to solve, and calculating the stray current distribution along the line. The method comprehensively considers the through characteristic of the long line of the urban rail transit, the through characteristic of the uplink and the downlink and the difference of different construction modes in the long line to calculate the stray current, so that the stray current distribution along the line of the urban rail transit can be accurately obtained, and a basis is provided for the design, construction and maintenance of the stray current protection along the line of the urban rail transit system.
Description
Technical Field
The invention belongs to the field of urban rail transit traction power supply systems, and particularly relates to a unified calculation method for urban rail transit stray currents.
Background
The urban rail transit power supply system in China generally adopts direct current power supply, and a train obtains electric energy from a contact net and returns to a negative electrode of a traction substation through a steel rail. In practical engineering, the rail has a longitudinal resistance to generate a ground potential, and the rail cannot be completely insulated from the ground, resulting in leakage of a part of the current from the rail, i.e., stray current. While the stray current corrodes the self structural facilities of the urban rail transit system, the stray current can also corrode nearby buried metal pipelines, thereby causing potential safety hazards.
In construction it is common to lay drainage nets under the rails to collect a portion of the stray current leaking from the rails. And structural steel bars are built in the tunnels through open-cut backfill construction and underground excavation construction, and the structural steel bars can be used as an auxiliary drainage net to further collect partial stray current. And the viaduct and shield construction line has no structural steel bars. The actual urban rail transit system is long in line and large in span, and different construction modes are adopted along the line due to factors such as environment and the like. And the urban rail transit reflux systems under different construction modes are different, so that the distribution rules of stray current along the line are different. In addition, the urban rail transit reflux system along the line is of a through structure, and the ascending and descending steel rails form the through structure through a flow equalizing line and a traction reflux line. And a stray current calculation method comprehensively considering the actual condition of the urban rail transit system is lacked at present. Therefore, it is necessary to analyze the stray current distribution by taking different construction methods and through structures along the urban rail transit into consideration in combination with the actual situation of the urban rail transit system and adopting a unified calculation method.
Disclosure of Invention
In order to solve the problems, the invention provides a unified calculation method for urban rail transit stray current.
The invention discloses a unified calculation method for urban rail transit stray current, which comprises the following steps of:
step 1: and monitoring and obtaining the train traction working condition, and obtaining relevant parameters of the urban rail transit reflux system.
Step 2: and constructing a unified calculation model of the stray current of the urban rail transit system.
And step 3: and constructing a function expression of the voltage and the current based on the model, and calculating the distribution of the stray current.
The step 1 specifically comprises the following steps:
acquiring the position WCU of an upstream train, the position WCD of the downstream train and the current I of the upstream train by monitoringcuDownstream train current IcdAnd obtaining the reflux current ascending the traction station, the reflux current descending the traction station and the current of the current equalizing line through load flow calculation.
And acquiring the position WT of the urban rail transit system along the line traction, the position WB of the flow equalizing line and the position interval WQ of the overhead or shield construction tunnel.
And measuring and obtaining relevant parameters of the urban rail transit reflux system, including unit resistance of an uplink steel rail, unit resistance of a downlink steel rail and unit resistance of an uplink drainage network.
The step 2 specifically comprises the following steps:
s21, regarding the urban rail transit reflux system as a two-dimensional plane, and equating the two-dimensional plane to be a resistance network structure; equivalently obtaining the resistance network structures under different construction modes according to different construction modes; and the current equalizing line, the train and the traction station are equivalent to a current source.
S22 obtaining the position W of the injection current source based on WCU, WCD, WB, WTmM is 1, …, M, and the comparison results in the upward injection current JS at the position corresponding to WmThe down injection current JDm(ii) a M represents the total number of injection current sources.
S23, any two adjacent current sources are regarded as a calculation interval, N represents the serial number of the calculation interval, and N is the total number of the calculation intervals.
S24 defines a variable Q(n)If the nth calculation interval is within the WQ interval, let Q be(n)On the contrary, let Q(n)＝1。
S25 calculates the current for each structure in the model: the current of the ascending steel rail, the current of the ascending current grid, the current of the structural steel bar, the current of the underground, the current of the descending current grid and the current of the descending steel rail are respectively and sequentially represented by Iru、Isu、Ie、Ig、Isd、IrdRepresenting, calculating the voltage difference of the adjacent structures in the model: voltage between up-running steel rail and up-running drainage net, up-running drainage net andthe voltage between the structural steel bars, the voltage between the structural steel bars and the equivalent ground, the voltage between the equivalent ground and the downstream drainage network, the voltage between the downstream drainage network and the downstream steel rail and the voltage between the upstream drainage network and the equivalent ground are respectively and sequentially controlled by Uru-su、Usu-e、Ue-g、Ug-sd、Usd-rd、Usu-gAnd (4) showing.
The step 3 specifically comprises the following steps:
s31 defines a current variable I1，I2，I3，I4，I5，I6Sequentially representing the current and voltage variables U of each structure in the model12，U23，U34，U45，U56Which in turn represents the voltage difference between adjacent structures. When Q is(n)When 2, I6＝0，U56＝0。
S32 solving U at any position x in the nth calculation interval based on kirchhoff' S law12，U23，U34，U45，U56，I1，I2，I3，I4，I5，I6The differential equation expression of (a) is as follows:
s33 matrix XnAnd (4) representing a square matrix formed by the steel rail, the drainage net, the structural steel bar, the earth equivalent resistance and the transition conductance in the nth calculation interval.
S34, solving the differential equation set to obtain U at any position x12，U23，U34，U45，U56And I1，I2，I3，I4，I5，I6The function of (2) is expressed by the formula F (x)(n)：
S35 finding the characteristic direction of XValue to CnAnd solving the eigenvalue assignment of X to L(n)(L1 (n),…,L11 (n))。
S36 setting a(n)(a1 (n),…,a11 (n)) For the coefficient to be solved of the nth calculation interval, an equation set is constructed for solving the coefficient a(n)：
When m is 1: let I1 (1)(W1)＝JS1,I2 (1)(W1)＝0,I3 (1)(W1)＝0,I4 (1)(W1)＝0,I5 (1)(W1)＝0,I6 (1)(W1)＝JD1
When 1< M:
if Q is(m-1)＝Q (m)1, order
If Q is(m-1)＝Q (m)2, order
If Q is(m-1)-Q (m)1, order
If Q is(m-1)-Q (m)1, order
When M ═ M: order to
I1 (N)(WM)＝JSM,I2 (N)(WM)＝0,I3 (N)(WM)＝0,I4 (N)(WM)＝0,I5 (N)(WM)＝0,I6 (N)(WM1)＝JDM
S37 obtaining each calculation interval n by joint solution based on all the equation sets of the step S36Coefficient a(n)。
S38 solving the obtained coefficient a(n)Substitution equation set F (x)(n)。
S39, assigning values to x, and solving to obtain the output of the equation set F (x), namely U, of each calculation interval12、U23、U34、U45、U56、I1、I2、I3、I4、I5、I6。
S310 stray current distribution of arbitrary calculation interval n
The stray current along the line is distributed as Is＝[Is (1),Is (2)…,Is (n)]。
The beneficial technical effects of the invention are as follows:
the method comprehensively considers the through characteristic of the long line of the urban rail transit, the through characteristic of the uplink and the downlink and the difference of different construction modes in the long line to calculate the stray current, so that the stray current distribution along the line of the urban rail transit can be accurately obtained, and a basis is provided for the design, construction and maintenance of the stray current protection along the line of the urban rail transit system.
Drawings
Fig. 1 is a schematic flow chart of a unified calculation method for stray current of urban rail transit.
FIG. 2 is a schematic diagram of an equivalent resistance network structure of a reflux system for open cut and underground cut construction according to the present invention.
Fig. 3 is a schematic diagram of an equivalent resistance network structure of a reflux system for overhead and shield construction.
Fig. 4 is a schematic diagram of the distribution of stray current of urban rail transit in the embodiment of the present invention.
Detailed Description
The invention is described in further detail below with reference to the figures and specific embodiments.
The invention discloses a unified calculation method of urban rail transit stray current, which is shown in figure 1 and specifically comprises the following steps:
step 1: and monitoring and obtaining the train traction working condition, and obtaining relevant parameters of the urban rail transit reflux system.
Acquiring the position WCU of an upstream train, the position WCD of the downstream train and the current I of the upstream train by monitoringcuDownstream train current IcdObtaining the reflux current I ascending from the traction station through load flow calculationstTraction station down reflux current IdtAnd current I of the current equalizerb。
And acquiring the position WT of the urban rail transit system along the line traction, the position WB of the flow equalizing line and the position interval WQ of the overhead or shield construction tunnel.
And measuring and obtaining relevant parameters of the urban rail transit reflux system. These parameters include: unit resistance value R of ascending steel railruUnit resistance value R of down railrdUnit resistance R of up-flow drainage networksuUnit resistance value R of down-flow drainage networksdStructural steel bar unit resistance ReEquivalent resistance R of earthdThe conductance G of the downstream drainage network of the upstream steel rail pairaudThe electric conductance G of the upstream drainage network of the downstream steel rail pairaduConductance G of ascending steel rail pair ascending drainage networkauThe down rail pair down drainage network conductance GadElectric conductance G of structural steel bar of ascending steel rail pairduElectric conductance G of down steel rail pair structural steel barddStructural steel bar conductance G of up-flow drainage net pairbuDown current drainage net pair structure steel bar electric conductance GbdGround electric conductor G of ascending steel railfuAnd the down going rail is earthed and electrically conductive GfdAnd the up-flow grid rail to earth electric conductor GeuAnd the down-flow grid rail to ground electric conductor GedElectric conductance G of ascending rail to descending railgConductance G of upstream drainage network to downstream drainage networkhStructural steel bar ground conductance Gc。
Step 2: and constructing a unified calculation model of the stray current of the urban rail transit system.
1. The urban rail transit backflow system is regarded as a two-dimensional plane and is equivalent to a resistance network structure. The reflux system of the overhead or shield construction is equivalent to a three-layer structure of a steel rail, a drainage net and the ground, and the reflux system of the open excavation method or the underground excavation method construction is equivalent to a four-layer structure of the steel rail, the drainage net, a structural steel bar and the ground. And the current equalizing line, the train and the traction station are equivalent to a current source injected into the reflux system, and as shown in fig. 2, the structure schematic diagram of the reflux system equivalent resistance network for open excavation and underground excavation construction is shown. Fig. 3 is a schematic diagram of an equivalent resistance network structure of a backflow system in overhead and shield construction. .
2. A virtual train is added to the downstream position corresponding to the upstream train position, and the train current is set to 0, and similarly, a virtual train is added to the upstream position corresponding to the downstream train position, and the train current is set to 0. The up and down train positions in the calculation model are represented by WC. I iscu-xRepresenting the up-train current in the calculation model, Icd-xRepresenting the downstream train current in the computational model.
3. Two-dimensional matrices CL _ CU (2, NC), CL _ CD (2, NC), two-dimensional matrix CL _ T (2, NT), two-dimensional matrix CL _ B (2, NB) are defined. NC, NT and NB are respectively the number of the upstream and downstream trains, the traction station and the flow equalizing line in the calculation model. Let CL _ CU (1, I) wc (I), CL _ CU (2, I) Icu-x(i) And in the same way, other two-dimensional matrixes are assigned with values.
4. The CL _ CU and the CL _ T, CL _ B are combined in groups to form a new matrix CU, and the CL _ CD and the CL _ T, CL _ B are combined in groups to form a new matrix CD. The CL is sorted in ascending order of rows based on the data of the first row of CL as a sorting criterion. Let Wm＝CU(1,m)，JSm＝CU(2,m)，JDmCD (2, M) M is 1,2, …, M. M represents the total number of injected current sources, WmTo inject the current source at the location, JSmFor upward injection of current, JDmThe current is injected downstream.
5. Any two adjacent current sources form a calculation interval, and N (1, …, N) represents a calculation interval number. N represents the total number of calculation intervals.
6. Defining variable Q(n)If the nth calculation interval is within the WQ interval, let Q be(n)On the contrary, let Q(n)＝1。
7. Upward steel rail current, upward current-discharging net current and knotThe current of the steel structure bar, the current in the ground, the current of the down current drainage network and the current of the down steel rail are respectively Iru、Isu、Ie、Ig、Isd、IrdThe voltage between the ascending steel rail and the ascending drainage network, the voltage between the ascending drainage network and the structural steel bar, the voltage between the structural steel bar and the equivalent ground, the voltage between the equivalent ground and the descending drainage network, the voltage between the descending drainage network and the descending steel rail and the voltage between the ascending drainage network and the equivalent ground are respectively represented by Uru-su、Usu-e、Ue-g、Ug-sd、Usd-rd、Usu-gAnd (4) showing.
And step 3: and constructing a function expression of the voltage and the current based on the model, and calculating the distribution of the stray current.
1. Defining a current variable I1，I2，I3，I4，I5，I6U, voltage variable12，U23，U34，U45，U56。
When Q is(n)When 1, I1＝Ig，I2＝Isu，I3＝Ie，I4＝Ig，I5＝Isd，I6＝Ird，U12＝Uru-su，U23＝Usu-e，U34＝Ue-g，U45＝Ug-sd，U56＝Usd-rd；
When Q is(n)When 2, I1＝Ig，I2＝Isu，I3＝Ig，I4＝Isd，I5＝Ird，I6＝0，U12＝Uru-su，U23＝Usu-g，U34＝Ug-sd，U45＝Usd-rd，U56＝0。
2. Based on kirchhoff's law, solving U at any position x in nth calculation interval12，U23，U34，U45，U56，I1，I2，I3，I4，I5，I6The differential equation expression of (a) is as follows:
3. matrix XnThe kirchhoff voltage relation and the kirchhoff current relation in the reaction resistance network are square matrixes formed by steel rails, drainage nets, structural steel bars, earth equivalent resistance and transition conductance in the nth calculation interval and are represented by Q(n)When 1 is taken as an example, XnThe assignment process is as follows:
wherein; x1＝-(Gau+Gdu+Gfu+Gaud+Gg)，X2＝-(Gdu+Gfu+Gaud+Gg)，X3＝-(Gfu+Gaud+Gg)，X4＝-(Gbu+Geu+Gh+Gadu)，X5＝-(Geu+Gh+Gadu)，X6＝-(Gc+Gbd+Gdd)，X7＝Gc+Geu+Gfu，X8＝Gbd+Gh+Gaud，X9＝Gdd+Gadu+Gg，X10＝Gad+Gfd+Gdd+Gadu+Gg。
4. Solving the differential equation set to obtain U at any position x12，U23，U34，U45，U56And I1，I2，I3，I4，I5，I6The function of (2) is expressed by the formula F (x)(n)：
5. Evaluating the eigenvector of X to CnAnd solving the eigenvalue assignment of X to L(n)(L1 (n),…,L11 (n))。
6、a(n)(a1 (n),…,a11 (n)) And calculating the coefficient to be solved of the interval for the nth. Constructing a system of equations for solving the coefficients a(n)：
When m is 1: let I1 (1)(W1)＝JS1,I2 (1)(W1)＝0,I3 (1)(W1)＝0,I4 (1)(W1)＝0,I5 (1)(W1)＝0,I6 (1)(W1)＝JD1
When 1< M:
if Q is(m-1)＝Q (m)1, order
If Q is(m-1)＝Q (m)2, order
If Q is(m-1)-Q (m)1, order
If Q is(m-1)-Q (m)1, order
When M ═ M: order to
I1 (N)(WM)＝JSM,I2 (N)(WM)＝0,I3 (N)(WM)＝0,I4 (N)(WM)＝0,I5 (N)(WM)＝0,I6 (N)(WM1)＝JDM
7. Based on the equation set in the step 6, the coefficient a of any calculation interval n is obtained through combined solution(n)。
8. The coefficient a obtained by solving(n)Substitution equation set F (x)(n)。
9. Assigning values to x, and solving to obtain the output of an equation set F (x), namely U, of each calculation interval12、U23、U34、U45、U56、I1、I2、I3、I4、I5、I6。
10. Stray current distribution for arbitrary calculation interval n
The stray current along the line is distributed as Is＝[Is (1),Is (2)…,Is (n)]。
Examples
By adopting the method, the stray current distribution of a certain urban rail transit system is calculated.
The parameters of the urban rail transit reflux system are shown in Table 1
TABLE 1 calculation of urban rail traffic reflux system parameter selection in an embodiment
The train traction conditions in the calculation example are shown in table 2:
TABLE 2 train traction condition selection
The up and down reflux currents and the current of the current equalizing line of the traction station obtained by load flow calculation are shown in table 3:
TABLE 3 calculation of the traction station up and down reflux currents and current of the current equalizer
A unified calculation model of stray current of an urban rail transit line is constructed according to parameters in the table 1, train traction working condition parameters are further obtained as shown in the table 2, and an upper return cable, a lower return cable and a current equalizing cable of a traction station are obtained through load flow calculation as shown in the table 3. The stray current distribution along the line calculated according to the method of the invention is shown in fig. 4.
Claims (4)
1. A unified calculation method for urban rail transit stray current is characterized by comprising the following steps:
step 1: monitoring and acquiring a train traction working condition, and acquiring related parameters of an urban rail transit reflux system;
step 2: constructing a unified calculation model of the stray current of the urban rail transit system;
and step 3: and constructing a function expression of the voltage and the current based on the model, and calculating the distribution of the stray current.
2. The unified calculation method for the stray current of the urban rail transit according to claim 1, wherein the step 1 specifically comprises:
acquiring the position WCU of an upstream train, the position WCD of the downstream train and the current I of the upstream train by monitoringcuDownstream train current IcdObtaining the reflux current ascending the traction station, the reflux current descending the traction station and the current of the current equalizing line through load flow calculation;
acquiring a position WT of a traction position along a line of the urban rail transit system, a position WB of a flow equalizing line and a position interval WQ of an overhead or shield construction tunnel;
and measuring and obtaining relevant parameters of the urban rail transit reflux system, including unit resistance of an uplink steel rail, unit resistance of a downlink steel rail and unit resistance of an uplink drainage network.
3. The unified calculation method for urban rail transit stray currents according to claim 2, wherein the step 2 specifically comprises:
s21, regarding the urban rail transit reflux system as a two-dimensional plane, and equating the two-dimensional plane to be a resistance network structure; equivalently obtaining the resistance network structures under different construction modes according to different construction modes; and the current equalizing line, the train and the traction station are equivalent to a current source;
s22 obtaining the position W of the injection current source based on WCU, WCD, WB, WTmM is 1, …, M, and the comparison results in the upward injection current JS at the position corresponding to WmThe down injection current JDm(ii) a M represents the total number of injection current sources;
s23, taking any two adjacent current sources as a calculation interval, wherein N represents the serial number of the calculation interval, and N is the total number of the calculation intervals;
s24 defines a variable Q(n)If the nth calculation interval is within the WQ interval, let Q be(n)On the contrary, let Q(n)＝1；
S25 calculates the current for each structure in the model: the current of the ascending steel rail, the current of the ascending current grid, the current of the structural steel bar, the current of the underground, the current of the descending current grid and the current of the descending steel rail are respectively and sequentially represented by Iru、Isu、Ie、Ig、Isd、IrdRepresenting, calculating the voltage difference of the adjacent structures in the model: the voltage between the ascending steel rail and the ascending drainage network, the voltage between the ascending drainage network and the structural steel bar, the voltage between the structural steel bar and the equivalent ground, the voltage between the equivalent ground and the descending drainage network, the voltage between the descending drainage network and the descending steel rail and the voltage between the ascending drainage network and the equivalent ground are respectively and sequentially measured by a Uru-su、Usu-e、Ue-g、Ug-sd、Usd-rd、Usu-gAnd (4) showing.
4. The unified calculation method for urban rail transit stray currents according to claim 3, wherein the step 3 specifically comprises:
s31 defines a current variable I1，I2，I3，I4，I5，I6Sequentially representing the current and voltage variables U of each structure in the model12，U23，U34，U45，U56Sequentially representing the voltage difference between adjacent structures; when Q is(n)When 2, I6＝0，U56＝0；
S32 solving U at any position x in the nth calculation interval based on kirchhoff' S law12，U23，U34，U45，U56，I1，I2，I3，I4，I5，I6The differential equation expression of (a) is as follows:
s33 matrix XnRepresenting a square matrix formed by the steel rail, the drainage net, the structural steel bar, the earth equivalent resistance and the transition conductance in the nth calculation interval;
s34, solving the differential equation set to obtain U at any position x12，U23，U34，U45，U56And I1，I2，I3，I4，I5，I6The function of (2) is expressed by the formula F (x)(n)：
S35 finding the characteristic vector of X and assigning the characteristic vector to CnAnd solving the eigenvalue assignment of X to L(n)(L1 (n),…,L11 (n))；
S36 setting a(n)(a1 (n),…,a11 (n)) For the coefficient to be solved of the nth calculation interval, an equation set is constructed for solving the coefficient a(n)：
When m is 1: let I1 (1)(W1)＝JS1,I2 (1)(W1)＝0,I3 (1)(W1)＝0,I4 (1)(W1)＝0,I5 (1)(W1)＝0,I6 (1)(W1)＝JD1
When 1< M:
if Q is(m-1)＝Q(m)1, order
If Q is(m-1)＝Q(m)2, order
If Q is(m-1)-Q(m)1, order
If Q is(m-1)-Q(m)1, order
When M ═ M: order to
I1 (N)(WM)＝JSM,I2 (N)(WM)＝0,I3 (N)(WM)＝0,I4 (N)(WM)＝0,I5 (N)(WM)＝0,I6 (N)(WM1)＝JDM
S37, based on all the equation sets in the step S36, the coefficient a of each calculation interval n is obtained through combined solution(n)；
S38 solving the obtained coefficient a(n)Substitution equation set F (x)(n)；
S39, assigning values to x, and solving to obtain the output of the equation set F (x), namely U, of each calculation interval12、U23、U34、U45、U56、I1、I2、I3、I4、I5、I6；
S310 stray current distribution of arbitrary calculation interval n
The stray current along the line is distributed as Is＝[Is (1),Is (2)…,Is (n)]。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110107753.2A CN112836364B (en) | 2021-01-27 | 2021-01-27 | Unified calculation method for stray current of urban rail transit |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110107753.2A CN112836364B (en) | 2021-01-27 | 2021-01-27 | Unified calculation method for stray current of urban rail transit |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112836364A true CN112836364A (en) | 2021-05-25 |
| CN112836364B CN112836364B (en) | 2022-08-19 |
Family
ID=75931783
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110107753.2A Active CN112836364B (en) | 2021-01-27 | 2021-01-27 | Unified calculation method for stray current of urban rail transit |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112836364B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113452007A (en) * | 2021-06-29 | 2021-09-28 | 西南交通大学 | Neutral point direct current calculation method for transformer of urban power grid |
| CN115600458A (en) * | 2022-10-11 | 2023-01-13 | 苏州大学（Cn） | Stray current dynamic distribution calculation method considering three-dimensional ground system |
| CN115600458B (en) * | 2022-10-11 | 2024-05-10 | 苏州大学 | Stray current dynamic distribution calculation method considering three-dimensional ground system |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102692885A (en) * | 2012-05-15 | 2012-09-26 | 西安理工大学 | Monitoring drainage control method of stray current |
| CN105303021A (en) * | 2014-07-24 | 2016-02-03 | 南车株洲电力机车研究所有限公司 | Adaptive real-time dynamic mathematic modeling method for direct-current traction power supply system in urban railway traffic |
| CN109256771A (en) * | 2018-10-16 | 2019-01-22 | 国网湖南省电力有限公司 | Metro stray current and its caused transformer neutral point DC current calculation method |
| CN110244810A (en) * | 2019-05-30 | 2019-09-17 | 深圳供电局有限公司 | The acquisition methods of Stray current values in tractive power supply system |
| CN110427703A (en) * | 2019-08-03 | 2019-11-08 | 国网山西省电力公司电力科学研究院 | A kind of stray electrical current emulation modelling method based under the conditions of multi-train movement |
| CN110763898A (en) * | 2019-10-16 | 2020-02-07 | 西南交通大学 | Method for acquiring stray current distribution of metro vehicle section |
| CN111274754A (en) * | 2020-03-05 | 2020-06-12 | 苏州大学 | Method and system for calculating power supply parameters of rail transit under non-uniform rail insulation condition |
| CN111324966A (en) * | 2020-03-04 | 2020-06-23 | 西南交通大学 | CDEGS-based urban rail transit stray current simulation calculation method |
| WO2020158527A1 (en) * | 2019-01-30 | 2020-08-06 | 株式会社日立製作所 | Stray current estimation system, stray current estimation method, and stray current estimation program |
| CN112230040A (en) * | 2020-10-13 | 2021-01-15 | 广州市扬新技术研究有限责任公司 | Stray current evaluation method for calculating direct current traction power supply system |
-
2021
- 2021-01-27 CN CN202110107753.2A patent/CN112836364B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102692885A (en) * | 2012-05-15 | 2012-09-26 | 西安理工大学 | Monitoring drainage control method of stray current |
| CN105303021A (en) * | 2014-07-24 | 2016-02-03 | 南车株洲电力机车研究所有限公司 | Adaptive real-time dynamic mathematic modeling method for direct-current traction power supply system in urban railway traffic |
| CN109256771A (en) * | 2018-10-16 | 2019-01-22 | 国网湖南省电力有限公司 | Metro stray current and its caused transformer neutral point DC current calculation method |
| WO2020158527A1 (en) * | 2019-01-30 | 2020-08-06 | 株式会社日立製作所 | Stray current estimation system, stray current estimation method, and stray current estimation program |
| CN110244810A (en) * | 2019-05-30 | 2019-09-17 | 深圳供电局有限公司 | The acquisition methods of Stray current values in tractive power supply system |
| CN110427703A (en) * | 2019-08-03 | 2019-11-08 | 国网山西省电力公司电力科学研究院 | A kind of stray electrical current emulation modelling method based under the conditions of multi-train movement |
| CN110763898A (en) * | 2019-10-16 | 2020-02-07 | 西南交通大学 | Method for acquiring stray current distribution of metro vehicle section |
| CN111324966A (en) * | 2020-03-04 | 2020-06-23 | 西南交通大学 | CDEGS-based urban rail transit stray current simulation calculation method |
| CN111274754A (en) * | 2020-03-05 | 2020-06-12 | 苏州大学 | Method and system for calculating power supply parameters of rail transit under non-uniform rail insulation condition |
| CN112230040A (en) * | 2020-10-13 | 2021-01-15 | 广州市扬新技术研究有限责任公司 | Stray current evaluation method for calculating direct current traction power supply system |
Non-Patent Citations (4)
| Title |
|---|
| ROMAN ŠPRLÁK 等: "Analysis of Stray Currents for Traction Vehicles", 《2014 14TH INTERNATIONAL CONFERENCE ON ENVIRONMENT AND ELECTRICAL ENGINEERING》 * |
| SHENG LIN 等: "Infinitesimal Method Based Calculation of Metro Stray Current in Multiple Power Supply Sections", 《IEEE ACCESS》 * |
| 楼锦君 等: "城轨钢轨电位与杂散电流动态分布", 《北京交通大学学报》 * |
| 薛皓: "基于可变电阻模块的地铁轨道电位及杂散电流模拟系统研究", 《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113452007A (en) * | 2021-06-29 | 2021-09-28 | 西南交通大学 | Neutral point direct current calculation method for transformer of urban power grid |
| CN113452007B (en) * | 2021-06-29 | 2022-04-29 | 西南交通大学 | Neutral point direct current calculation method for transformer of urban power grid |
| CN115600458A (en) * | 2022-10-11 | 2023-01-13 | 苏州大学（Cn） | Stray current dynamic distribution calculation method considering three-dimensional ground system |
| CN115600458B (en) * | 2022-10-11 | 2024-05-10 | 苏州大学 | Stray current dynamic distribution calculation method considering three-dimensional ground system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112836364B (en) | 2022-08-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Charalambous et al. | Dynamic stray current evaluations on cut-and-cover sections of DC metro systems | |
| WO2020158527A1 (en) | Stray current estimation system, stray current estimation method, and stray current estimation program | |
| Charalambous et al. | A holistic stray current assessment of bored tunnel sections of DC transit systems | |
| CN107609208B (en) | Traction network modeling method considering tunnel section comprehensive grounding system structure | |
| Charalambous et al. | Modeling for preliminary stray current design assessments: The effect of crosstrack regeneration supply | |
| CN110763898B (en) | Method for acquiring stray current distribution of metro vehicle section | |
| CN106019077B (en) | A kind of current mode traveling wave ranging device Optimizing method | |
| CN104901030A (en) | Construction method of resistance reduction and grounding device of transformer station | |
| CN112836365B (en) | Unified calculation method for flow distribution of urban rail transit traction power supply system | |
| CN112836364B (en) | Unified calculation method for stray current of urban rail transit | |
| Zhou et al. | A uniform model for stray current of long-line DC metro systems | |
| Cerman et al. | Resistive-type network model of stray current distribution in railway DC traction system | |
| CN109149563B (en) | Method for determining stray current value in alternating current power grid | |
| CN109975596A (en) | Earth current is distributed research method under a kind of monopole the earth method of operation | |
| CN112348359A (en) | Method for determining safety design index of steel pipe tower grounding device of urban power transmission line | |
| CN113452007B (en) | Neutral point direct current calculation method for transformer of urban power grid | |
| Aslam et al. | Design analysis and optimization of ground grid mesh of extra high voltage substation using an intelligent software | |
| CN115422691A (en) | Rail transit stray current model obtaining method considering train traction characteristics | |
| Hoger et al. | Influence of high voltage power lines on the propagation of stray currents from DC traction | |
| Höger et al. | Investigating the impacts of high voltage powerlines on ground current fields | |
| CN117113733B (en) | Method and device for acquiring bias current of direct current near zone of power grid | |
| CN114709782B (en) | Method for reducing potential difference in transformer substation grounding grid caused by impact current | |
| Nahman et al. | Safety analysis at overhead line towers in close proximity to the substation | |
| CN114707351B (en) | Subway transfer station grounding grid equivalent method suitable for CDEGS simulation | |
| Pradeep et al. | Comprehensive Design of an Earthing System for Broadlands Hydropower Station Switchyard |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |