CN111523800A

CN111523800A - Rapid calculation method for node conductance matrix in subway load flow calculation

Info

Publication number: CN111523800A
Application number: CN202010321745.3A
Authority: CN
Inventors: 张超; 李红波; 黄子昊; 罗文广; 刘华东; 梅文庆
Original assignee: CRRC Zhuzhou Institute Co Ltd
Current assignee: CRRC Zhuzhou Institute Co Ltd
Priority date: 2020-04-22
Filing date: 2020-04-22
Publication date: 2020-08-11
Anticipated expiration: 2040-04-22
Also published as: CN111523800B

Abstract

The invention provides a method for quickly calculating a node conductance matrix in subway load flow calculation, which comprises the following steps: sequentially sequencing and continuously numbering a traction substation, an ascending train and a descending train according to the distance relative to the starting point of the ascending train, and dividing power supply intervals; construction includes Y_SS、Y_SU、Ys_D、Y_US、Y_UU、Y_DSAnd Y_DDAnd setting each sub-matrix to zero; processing the mark vector according to the relation between the number of times of train cycle counting and the number of train vehicles and the power supply section to which the train belongs, and carrying out Y-point matching_SSPerforming a correction of Y_SSCalculating Y_SUAnd Y_UUOr calculating Y_UUOr correcting Y_SSCalculating Y_SDAnd Y_DDOr calculating Y_DD. By single layer circulationAnd elements in the node sub-conductance matrix are calculated, repeated calculation of the elements is avoided, and the calculation efficiency of the algorithm is improved.

Description

Rapid calculation method for node conductance matrix in subway load flow calculation

Technical Field

The invention relates to the technical field of subway load flow simulation calculation, in particular to a method for quickly calculating a node conductance matrix in subway load flow calculation.

Background

With the rapid development of economy in China, subways appear in most cities in China with the advantages of large passenger capacity, no congestion, no delay and the like. The method is important for subway operators to guarantee the safety of subway power supply during operation, and the load flow calculation is an important means for guaranteeing analysis.

In the current subway load flow calculation simulation, for the calculation of a node conductance matrix, the main idea is to respectively sort and calculate the node conductance matrix for an uplink contact network or a downlink contact network, and finally, draw substation nodes of the uplink contact network and the downlink contact network; in addition, multiple double-layer loop nesting calculation is adopted in the calculation of the independent uplink overhead contact system and the downlink overhead contact system, so the calculation amount of the method is very large undoubtedly.

Therefore, a method capable of quickly and accurately calculating a node matrix in the subway tide is needed.

Disclosure of Invention

In view of the above, the invention aims to provide a method for quickly calculating a node conductance matrix in subway load flow calculation, so as to solve the problems of very large subway load flow calculation amount and slow speed.

Based on the purpose, the method for quickly calculating the node conductance matrix in the subway load flow calculation comprises the following steps:

s100, providing the distance between each element and the starting point of an ascending train, the number of each element and resistance parameters, wherein each element comprises a traction substation, the ascending train and a descending train, and the resistance parameters comprise a contact network unit resistance, a resistive steel rail unit resistance and a leakage unit resistance of a rail to the ground;

s200, sequentially and respectively sequencing and numbering traction substations according to the distance relative to the starting point of the uplink train, dividing the interval between the starting point and the end point of the train line into a plurality of power supply intervals, and calculating the number and the length of the power supply intervals;

s300, constructing a node conductance matrix which comprises a node sub conductance matrix Y between the traction substation and the traction substation_SSNode conductance matrix Y between traction substation and uplink train_SUNode conductance matrix Y between traction substation and down train_SDNode conductance matrix Y between ascending train and traction substation_USAnd a node sub conductance matrix Y between the ascending trains_UUNode conductance matrix Y between down train and traction substation_DSAnd a node sub-conductance matrix Y between the downstream trains_DDAnd is combined with Y_SS、Y_SU、Ys_D、Y_US、Y_UU、Y_DSAnd Y_DDSetting zero;

s400, calculating Y according to the non-train operation_SS；

S500, sequencing and numbering the trains in sequence according to the distance relative to the starting point of the uplink train, wherein the number of the uplink train is continuous with the number of the traction substation, and the number of the downlink train is continuous with the number of the uplink train;

s600, clearing an uplink train number vector, an uplink power supply interval processing flag vector, an uplink train left end element index number and an uplink train right end element index number in a power supply interval, and setting an initial value of the uplink train cycle count number; circularly calculating the index number of the element at the left end of the ascending train and the index number of the element at the right end of the ascending train according to the relation between the number of times of circularly counting the ascending trains and the number of the ascending trains; and to Y_SSMaking a correction or a correction of Y_SSCalculating Y_SUAnd Y_UUOr calculating Y_UUOr calculating Y_US。

S700, clearing the downlink train number vector, the downlink power supply section processing flag vector, the downlink train left end element index number and the downlink train right end element index number in the power supply section, and setting the initial value of the downlink train cycle count number(ii) a Circularly calculating the index number of the element at the left end of the descending train and the index number of the element at the right end of the descending train according to the relation between the number of times of circularly counting the descending trains and the number of the descending trains; and to Y_SSMake a correction of Y_SDAnd Y_DDOr calculating Y_DDOr calculating Y_DS。

In one embodiment, step S400 includes,

by passing

Calculating Y_SSWherein i is₁And j₁Numbering nodes, j₁∈i₁Represents directly with i₁Connected nodes, r₁Is the unit resistance of the contact line, L_tIs the length of the t-th power supply interval, t is the number of the power supply interval, n is the total number of the traction substation, r₂Is the unit resistance of the running rail, r₃The unit resistance of the leakage of the steel rail to the ground.

In one embodiment, for Y_SSThe correcting comprises the following steps: when the relation between the number of times of the cycle counting of the upper train and the number of the upper train satisfies i ═ m +1, the processing flag vector dealup (j) of the power supply section of the upper train is equal to 0, and the index number of the power supply section of the upper train satisfies j ≠ 1 and j ≠ n +1, or when the relation between the number of times of the cycle counting of the train and the number of the lower train satisfies i ═ k +1, the processing flag vector dealdown (j) of the power supply section of the lower train is equal to 0, and the index number of the power supply section satisfies j ≠ 1 and j ≠ n +1,

by passing

For Y_SSCorrecting, wherein i is the number of times of cycle counting of the ascending trains or the descending trains, m is the number of the ascending trains, k is the number of the ascending trains, n is the total number of the traction substation, n +1 is the total number of the power supply interval, and L_jThe length of the jth power supply interval to which the ascending train or the descending train belongs.

In one embodiment, in step S600, Y is modified_SSCalculating Y_SUAnd Y_UUWhen the relation between the number of times of cycle counting of the train and the number of the ascending trains meets i ═ m +1, i ≠ 1, and the right-end element index number Turight (i) of the ascending train meets 0<When Turight (i-1) is less than or equal to n; or when the relation between the number of times of train cycle counting and the number of the ascending trains meets i ≠ m +1, and the index number Tupleft (i) of the element at the left end of the ascending train meets 0<When the content of Tupleft (i) is less than or equal to n,

by passing

By passing

By passing

Wherein, a, b, c and d are index numbers, a ═ i, b ═ Tupright (i-1), c ═ b, d ═ a, dis ═ i-1 distance between the i-th ascending train and the ascending train starting point — Tupright (i-1) th distance between the traction substation and the ascending train starting point, or a ═ Tupleft (i), b ═ i, c ═ a, d ═ b, dis ═ i distance between the i-th ascending train and the ascending train starting point — Tupleft (i ═ Tupleft (i), b ═ i, c ═ a, d ═ b, dis ═ i distance between the i-th ascending train and the ascending train₁) And the distance between each traction substation and the starting point of the ascending train.

In one embodiment, Y is calculated in step S600_UUIncludes that when the relation between the number of times of cycle counting of the train and the number of the ascending trains meets i ═ m +1, i ≠ 1, and the right-end element index number Tupright (i) of the ascending train meets Tupright (i-1)>When n is greater than n; or when the relation between the number of times of train cycle counting and the number of the ascending trains meets i ≠ m +1, and the index number Tupleft (i) of the ascending train left end element meets Tupleft (i)>When the number n is greater than the predetermined value,

passing through type

Calculating Y_UUWhere a ═ i, b ═ Tupright (i-1) -n, d ═ a, dis ═ i uplink vehicle is relative to the ith uplink vehicleThe distance between the starting point of the ascending train and the starting point of the ascending train is (i-1) th ascending train relative to the distance between the starting point of the ascending train or (a) tupleft (i) -n, b-i, d-b, dis) is the distance between the ith ascending train and the starting point of the ascending train and is (i-1) th ascending train relative to the distance between the starting point of the ascending train.

In one embodiment, step S600 calculates Y_USIncluding that the relation between the number of times of circularly counting the train on line and the number of the train on line satisfies i>m +1, calculating Y_SUIs transposed to obtain Y_US。

In one embodiment, step S700 calculates Y_SS、Y_SDAnd Y_DDIncludes that the relation between the number of times of train cycle counting and the number of downstream trains satisfies i ═ k +1, i ≠ 1, and the downstream train right-end element index number Tdown right (i) satisfies 0<When Tdown (i-1) is less than or equal to n; or when the relation between the number of train cycle counting and the number of the downlink train is i ≠ i +1, and the index number Tdownleft (i) of the element at the left end of the downlink train satisfies 0<When Tdownleft (i) is less than or equal to n,

by passing

By passing

By passing

Wherein, a, b, c and d are index numbers, a is i, b is tdown (i-1), c is b, d is a, dis is the distance between the i-1 th downstream train and the upstream train starting point-the tdown (i-1) th traction substation, or a is tdown (i), b is i, c is a, d is b, dis is the distance between the ith downstream train and the upstream train starting point-tdown (i)₁) And the distance between each traction substation and the starting point of the ascending train.

In one embodiment, Y is calculated in step S700_DDIncluding the number of times and the number of times of current train-train cycle countingThe relationship of the number of train vehicles satisfies i ═ k +1, i ≠ 1, and the downline right end train element index number tdown right (i) satisfies tdown right (i-1)>When n is greater than n; or the relation between the number of times of current train/queue cycle counting and the number of the downstream train is i ≠ k +1, and the index number Tdownleft (i) of the downstream train left end element satisfies Tdownleft (i)>When the number n is greater than the predetermined value,

by passing

Calculating Y_DDWherein, a ═ i, b ═ tdown right (i-1) -n, d ═ a, dis ═ ith descending vehicle distance from the ascending train starting point — i-1 th descending vehicle distance from the ascending train starting point or a ═ tdownleft (i) -n, b ═ i, d ═ b, dis ═ i descending vehicle distance from the ascending train starting point — i-1 th descending vehicle distance from the ascending train starting point.

In one embodiment, step S700 further includes that the relationship between the number of next train circulation counts and the number of next trains satisfies i>k +1, calculating Y_SDIs transposed to obtain Y_DS。

In one embodiment, the calculating the ascending train left end element index number and the ascending train right end element index number or the calculating the descending train left end element index number and the descending train right end element index number includes:

when the relation between the number of times of the upper train-train cycle counting and the number of the upper trains meets i < m +1 or the relation between the number of times of the lower train-train cycle counting and the number of the lower trains meets i < k +1, judging the power supply area of the upper train according to the distance between the upper train and the starting point of the upper train or judging the power supply area of the lower train according to the distance between the lower train and the starting point of the lower train, and according to the relation

Setting an uplink train left end element index and an uplink train right end element index or setting a downlink train left end element index and a downlink train right end element index, wherein i is the number of times of cycle counting of an uplink train or a downlink train, m is the number of the uplink trains, k is the number of the uplink trains, n is the total number of traction substations, n +1 is the total number of power supply sections, j is the index number of the power supply section to which the uplink train or the downlink train belongs, x is the j-th power supply section left index number, y is the j-th power supply section right index number, x is 0 to represent that no traction substation exists at the left end, and y is 0 to represent that no traction substation exists at the right end.

As can be seen from the above, according to the method for rapidly calculating the node conductance matrix in the subway load flow calculation provided by the invention, the traction substation, the ascending train and the descending train are sequentially numbered according to the distance from the starting point of the ascending train, and the power supply intervals are divided according to the traction substation, so that the node sub conductance matrix Y including the distance between the traction substation and the traction substation is constructed_SSNode conductance matrix Y between traction substation and uplink train_SUNode conductance matrix Y between traction substation and down train_SDNode conductance matrix Y between ascending train and traction substation_USAnd a node sub conductance matrix Y between the ascending trains_UUNode conductance matrix Y between down train and traction substation_DSAnd a node sub-conductance matrix Y between the downstream trains_DDAnd calculating each node sub-conductance matrix according to the relationship between the number of train cycle counts and the number of train vehicles, the train left end element index number and the train right end element index number. The elements in the node sub-conductance matrix are calculated through single-layer circulation, the elements in the node sub-conductance matrix of the traction substation and the node sub-conductance matrix of the traction substation are calculated only once, and meanwhile, the repeated calculation of the elements in the node conductance matrix is avoided and the calculation efficiency of the algorithm is improved by utilizing the node conductance matrix between the traction substation and the traction substation in the last time of continuous operation of the train in the same power supply interval.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.

Fig. 1 is a schematic diagram of a dc power supply network for a subway according to an embodiment of the present invention;

fig. 2 is a schematic flow chart of a method for rapidly calculating a node conductance matrix in subway power flow calculation according to an embodiment of the present invention.

Detailed Description

In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to specific embodiments and the accompanying drawings.

It is to be noted that technical terms or scientific terms used in the embodiments of the present invention should have the ordinary meanings as understood by those having ordinary skill in the art to which the present disclosure belongs, unless otherwise defined. The use of "first," "second," and similar terms in this disclosure is not intended to indicate any order, quantity, or importance, but rather is used to distinguish one element from another. The word "comprising" or "comprises", and the like, means that the element or item listed before the word covers the element or item listed after the word and its equivalents, but does not exclude other elements or items. The terms "connected" or "coupled" and the like are not restricted to physical or mechanical connections, but may include electrical connections, whether direct or indirect. "upper", "lower", "left", "right", and the like are used merely to indicate relative positional relationships, and when the absolute position of the object being described is changed, the relative positional relationships may also be changed accordingly.

As shown in fig. 1, a schematic diagram of four-line models of an uplink overhead line system, an uplink steel rail, a downlink overhead line system and a downlink steel rail, which are adopted in most current methods in a subway direct-current power supply network, is shown. Wherein, Veh represents a train including an up train and/or a down train, and TTS is a traction substation.

The inventor finds that the calculation speed of the load flow algorithm is particularly critical aiming at the characteristic that the structure of a subway direct-current power supply network changes constantly in long-term subway load flow calculation work, and the speed of the algorithm mainly depends on the formation of a network node conductance matrix and the iterative cycle of voltage. Aiming at the formation of the node conductance matrix, the research method has larger difference and more direct influence on the algorithm.

The inventor also finds that for the structural design and operation of the node conductance matrix, 4-node modeling of a train and a traction substation is considered in the method, namely, upper and lower contact network nodes, upper and lower steel rail nodes and upper and lower nodes of the traction substation are connected through small resistors, so that the modeling matrix in the method is very large in scale, and the potentials of a large number of virtual nodes are calculated; in the other method, 2-node modeling of a train and a traction substation is considered, an uplink contact network node and a steel rail node or a downlink contact network node and the steel rail node are considered, and finally the uplink traction substation node and the downlink traction substation node are combined. Such patterns may result in repeated operation of the traction substation nodes. Meanwhile, for the element calculation of the node conductance matrix, the traditional method mainly adopts a traction substation and a train to carry out uniform sequencing and numbering, and finally adopts double-layer circulation to calculate the elements in the node conductance matrix one by one, so that a large amount of zero elements can be calculated.

The inventor provides a method for quickly calculating a node conductance matrix in subway load flow calculation, wherein the node conductance matrix is preferentially numbered according to the characteristic of fixed position of a traction substation, and the numbering sequence of network elements is as follows: traction substation, ascending train and descending train. And the node matrix is calculated in blocks, which is specifically divided into: the node conductance matrix between the traction substation and the traction substation, the node conductance matrix between the traction substation and the ascending train, the node conductance matrix between the traction substation and the descending train, the node conductance matrix between the ascending train and the ascending train, and the node conductance matrix between the descending train and the descending train. And (4) pre-calculating a node conductance matrix between the traction substation and the traction substation, and dynamically modifying in a train insertion mode. And establishing a train left index vector and a train right index vector by utilizing the characteristic that the train is connected only by the nodes on the left and the right. The calculation of the node conductance matrix can be optimized, the speed of load flow calculation can be improved, and the formation and calculation of the node conductance matrix in subway load flow calculation can be greatly optimized.

Referring to fig. 2, a method for rapidly calculating a node conductance matrix in a subway power flow calculation according to an embodiment of the present invention includes:

s400, calculating Y according to the non-train operation_SS；

s600, clearing an uplink train number vector, an uplink power supply interval processing flag vector, an uplink train left end element index number and an uplink train right end element index number in a power supply interval, and setting an initial value of the uplink train cycle count number; circularly calculating the index number of the element at the left end of the ascending train and the index number of the element at the right end of the ascending train according to the relation between the number of times of circularly counting the ascending trains and the number of the ascending trains; and to Y_SSMaking a correction or a correction of Y_SSCalculating Y_SUAnd Y_UUOr calculating Y_UUOr calculating Y_US；

S700, clearing a downlink train number vector, a downlink power supply interval processing flag vector, a downlink train left end element index number and a downlink train right end element index number in a power supply interval, and setting an initial value of the downlink train cycle count number; circularly calculating the index number of the element at the left end of the descending train and the index number of the element at the right end of the descending train according to the relation between the number of times of circularly counting the descending trains and the number of the descending trains; and to Y_SSMaking a correction or a correction of Y_SSCalculating Y_SDAnd Y_DDOr calculating Y_DDOr calculating Y_DS。

The method for rapidly calculating the node conductance matrix in the subway load flow calculation sequentially numbers a traction substation, an ascending train and a descending train according to the distance relative to the starting point of the ascending train, divides power supply intervals according to the traction substation, and constructs the node conductance matrix Y between the traction substation and the traction substation_SSNode conductance matrix Y between traction substation and uplink train_SUNode conductance matrix Y between traction substation and down train_SDNode conductance matrix Y between ascending train and traction substation_USAnd a node sub conductance matrix Y between the ascending trains_UUNode conductance matrix Y between down train and traction substation_DSAnd a node sub-conductance matrix Y between the downstream trains_DDThe node sub-conductance matrix of (2) is based on the relation between the number of train cycle counts and the number of train vehicles, and the train left end elementAnd calculating the sub conductance matrix of each node in the index number and the index number of the train right-end element. The elements in the node sub-conductance matrix are calculated through single-layer circulation, the elements in the node sub-conductance matrix of the traction substation and the node sub-conductance matrix of the traction substation are calculated only once, and meanwhile, the repeated calculation of the elements in the node conductance matrix is avoided and the calculation efficiency of the algorithm is improved by utilizing the node conductance matrix between the traction substation and the traction substation in the last time of continuous operation of the train in the same power supply interval.

The index number of the defined element is i, namely i is the index number of a traction substation, a power supply area or an uplink train or a downlink train. The distance between the element and the start point of the ascending train is defined as ds, that is, ds is the distance between the traction substation, the power supply section, the ascending train or the descending train and the start point of the ascending train. Each element has two nodes, specifically, two ends of each element correspond to one node, that is, each traction substation has an upper node and a lower node, and each ascending train or descending train has an upper node and a lower node. In the formula of the loop calculation, for example, in the variable + constant/parameter, the variable on the left side of the formula and the variable on the right side of the formula are different values of the same variable, that is, the variable on the left side of the formula is the value of the variable calculated in the current loop, the variable on the right side of the formula is the value of the variable calculated in the previous loop, and when the variable on the left side of the formula is the value calculated in the first loop, the variable on the right side of the formula is the initial value.

In step S200, the traction substations may be numbered sequentially according to the distance from the start point of the ascending train from small to large or the distance from the start point of the ascending train from large to small. For example, when the total number of traction substations is n, the traction substations are numbered 1 to n in sequence, and the section between the start point and the end point of the train line is divided into n +1 power supply sections.

Specifically, when the traction substation is numbered according to the distance from the start point of the ascending train from small to large, the division of the power supply section, that is, the distance from the power supply section to the start point of the ascending train is calculated by the following formula (1). Wherein n is the total number of the traction substation, n +1 is the total number of the power supply interval, sr_iRepresents the firstDistance, ds, of i power supply sections from the start of the ascending train_iRepresents the distance, ds, of the ith traction substation from the start of the up-train_i+1Represents the distance, ds, of the i +1 th traction substation from the start of the upstream train_i-1And represents the distance between the i-1 st traction substation and the starting point of the ascending train.

Further, the length of the 2 nd to nth power feeding sections in the power feeding sections is calculated by the formula (2), wherein L_iFor the length of the supply interval, L_i＝ds_i-ds_i-1,1<i is less than or equal to n (2). The length of the 1 st power supply section is the distance from the 1 st traction substation to the starting point of the ascending train, for example, the starting point of the ascending train, and the length of the nth power supply section is the distance from the nth traction substation to the end point of the ascending train.

In step S300, the form of the node conductance matrix is shown in formula (3). And Y is a node conductance matrix of the whole subway direct-current power supply network. Y is_SSThe node sub conductance matrix between the traction substation and the traction substation. Y is_SUThe node sub-conductance matrix between the traction substation and the uplink train is formed. Y is_SDThe node sub-conductance matrix between the traction substation and the downlink train. Y is_USIs a node sub-conductance matrix, Y, between an up train and a traction substation_UUIs a node conductance matrix, Y, between an ascending train and an ascending train_DSFor the node conductance matrix, Y, between the down train and the traction substation_DDThe node sub-conductance matrix is between the downstream trains.

The node conductance matrix in the form is constructed by only considering two nodes in a network of the traction substation, the uplink train and the downlink train and enabling the node numbers in the uplink network and the downlink network of the same traction substation to be the same, so that the memory space of the node conductance matrix can be reduced, the formation of the network node conductance matrix is accelerated, and the calculation speed of the algorithm is further improved.

Calculating the node conductance matrix Y of the whole subway direct-current power supply network, firstly calculating Y when no train runs_SSAnd sub-conductance matrix Y of other nodes_SU、Ys_D、Y_US、Y_UU、Y_DSAnd Y_DDSet to zero (i.e., step S400). Then, numbering the ascending trains and the descending trains in sequence according to the condition that the trains run (namely, the form of inserting trains) (namely, step S500); circularly calculating the index number of the left end element and the index number of the right end element of the ascending train, and circularly correcting Y_SSCalculating Y_SS、Y_SU、Y_UUAnd Y_US(i.e., step S600); circularly calculating the index number of the left end element and the index number of the right end element of the downlink train, and circularly correcting Y_SSCalculate Ys_D、Y_DDAnd Y_DS(i.e., step S700).

In step S400, specifically, Y is calculated by the equation (4)_SS. Wherein i₁And j₁Numbering nodes, j₁∈i₁Represents directly with i₁Connected nodes, r₁Is the unit resistance of the contact line, L_tIs the length of the t-th power supply interval, t is the number of the power supply interval, n is the total number of the traction substation, r₂Is the unit resistance of the running rail, r₃The unit resistance of the leakage of the steel rail to the ground.

In step S500, the trains may be numbered according to the same rule as the sequence numbers of the traction substations. For example, when the traction substation numbers the distance from the start point of the ascending train from small to large, the trains are also numbered from small to large. When the trains are numbered, the serial numbers of the ascending trains and the serial numbers of the descending trains are sequentially numbered, so that the serial numbers of the ascending trains are continuous with the serial numbers of the traction substation, and the serial numbers of the descending trains are continuous with the serial numbers of the ascending trains.

In step S600, before the uplink train left and right element index numbers are calculated in a loop, the number of upper train vehicles in the power supply section, the uplink power supply section processing flag vector, and the uplink train element left end index vector and the uplink train element right end index vector are all set to zero, and the initial value of the number of times of the uplink train loop count is set to 1. Specifically, the setting may include a clear uplink train number vector nreptrain in the power supply interval, the dimension of which is (n +1) × 1, and a clear uplink power supply interval processing flag vector dealup, the dimension of which is (n +1) × 1, where n is the number of traction substations; clearing a left index vector Tupleft of the uplink train element, wherein the dimension is mx 1; clearing the right index vector Turight of the uplink train element, wherein the dimension is mx 1; the initial value i of the number of times of the up train loop count is defined to be 1. Wherein m is the number of vehicles in the upper row.

And after the setting is finished, calculating the index number of the element at the left end of the train, the index number of the element at the right end of the train and the node sub conductance matrix according to the relation between the number of the real-time uplink train cycle counting times and the number of the uplink trains.

The specific calculation includes:

circularly judging the relation between the number i of the real-time uplink train cycle counting and the number m of the uplink trains;

s610, firstly, for i>m + 1: if yes, the process proceeds to step S670, where Y is calculated_US(ii) a If the judgment result is no, that is, i is less than or equal to m +1, the process proceeds to step S620, and continues to judge that i is m + 1;

s620, if the judgment is negative, i<And m +1, judging the jth power supply section to which the ith (namely the index number of the real-time uplink train represented by the number i of real-time train cycle counting) train belongs according to the formula (1), and setting the left index number and the right index number of the power supply section to which the uplink train belongs, and the left end element index number and the right end element index number of the uplink train according to the formula (5). Wherein x is the left index number of the jth power supply interval, y is the right index number of the jth power supply interval, and the numerical value is equal to 0 to represent the left or rightNo traction substation; and n is the total number of the traction substations.

When the initial upstream train number vector nrupprain (j) is equal to 0, the actually counted upstream train number vector nrupprain (j) satisfies nrupprain (j) ═ nrupprain (j) + 1; tupleft (i) x, tupright (i) y; nreptrain (j) 1 when nreptrain (j) is not equal to 0; tupleft (i) n + i-1, Tupright (i) y, Tupright (i-1) n + i. Wherein n is the total number of the traction substation, i is the number of cycle counting, j is the number of the power supply section to which the uplink train belongs, x is the index number of the element at the left end of the uplink train, and y is the index number at the right end of the uplink train.

If yes, i is m +1, the process proceeds to step S630, a processing flag vector dealup (j) of the power supply section to which the up train belongs and an index number j of the power supply section to which the up train belongs are determined, and if dealup (j) is equal to 0, j ≠ 1 and j ≠ n +1 are satisfied, Y obtained in step S400 is determined by equation (6)_SSMaking a correction, after which dealup (j) is equal to 1; not correcting Y when dealup (j) is equal to 0, but j ≠ 1 and j ≠ n +1 are not satisfied_SS。

In the formula (6), i is the number of times of the cycle counting of the uplink trains, m is the number of the uplink trains, k is the number of the uplink trains, n is the total number of the traction substation, n +1 is the total number of the power supply section, and L_jThe length of the jth power supply interval to which the uplink train belongs. Y is_SS(2j-3,2j-1)、Y_SS(2j-1,2j-3)、Y_SS(2j-2,2j)、Y_SS(2j,2j-2)、Y_SS(2j-3,2j-3)、Y_SS(2j-1,2j-1)、Y_SS(2j-2 ) and Y_SS(2j,2j) are each Y_SSOf the different elements.

S640, determining that i is m +1, and when the determination is no, i is m +1, performing further determination and calculation in step S640; if the determination is yes, i is m +1, the process proceeds to step S650, and the determination is continued with i being 1, and if the determination is no, i is m +1 and i is not equal to 1, the further determination and calculation in step S650 are performed. If yes, i is 1, the process proceeds to step S660.

The further judgment and calculation in step S640 includes judgment of the index number tupleft (i) of the train at the left end of the ascending train:

when tupleft (i) is judged to be 0, no processing is performed, namely no calculation is performed;

when the judgment is 0<When Tupleft (i) is less than or equal to n, calculating Y by formula (7)_SSCalculating Y by equation (8)_SU(ii) a Calculating Y by equation (10)_UU(ii) a The distance between a, b, c and d is an index number, a is tupleft (i), b is i, c is a, d is b, dis is the distance between the ith ascending train and the ascending train starting point-tupleft (i) distance between the traction substation and the ascending train starting point;

when the judgment is Tupleft (i)>When n is in the above range, Y is calculated from the above formulas (9) and (10)_UU(ii) a Wherein a, b, c and d are index numbers, a is tupleft (i) -n, b is i, d is b, dis is the distance between the ith ascending vehicle and the ascending train starting point-the distance between the (i-1) th ascending vehicle and the ascending train starting point;

the further determination and calculation in step S650 includes determining the up train right end element index number turight (i):

when the judgment result is that Tupright (i-1) is 0, no processing is performed, that is, no calculation is performed;

when the judgment is 0<When Turight (i-1) is not more than n, calculating Y by formula (7)_SSCalculating Y by equation (8)_SU(ii) a Calculating Y by equation (10)_UU(ii) a The distance between the i-1 th uplink train and the start point of the uplink train-the distance between the i-1 th traction substation and the start point of the uplink train is-Tupright (i-1);

when it is judged to be Turight (i-1)>When n is in the above range, Y is calculated from the above formulas (9) and (10)_UU(ii) a Wherein, a, b, c and d are index numbers, a is i, b is Tupright (i-1) -n, d is a, dis is the distance from the ith ascending vehicle to the ascending train starting point-the distance from the (i-1) th ascending vehicle to the ascending train starting point.

S660: the number of times i of the up train loop count is set to i +1, and the process returns to step S610 to perform loop determination and calculation.

S670: calculating Y_US: by calculating Y_SUIs transposed to obtain Y_US。

That is, in step S600, when the relationship between the number of train cycles counted in the upper train and the number of ascending trains satisfies i ═ m +1, the processing flag vector dealup (j) of the ascending train power supply section is equal to 0, and the index number of the ascending train power supply section satisfies j ≠ 1 and j ≠ n +1, Y is calculated by equation (6)_SSCorrecting;

when the relation between the number of times of the cycle counting of the upper train and the number of the ascending trains meets i ═ m +1, and i ≠ 1, and meanwhile, the ascending train right-end element index number Turight (i) meets 0<When Turight (i-1) is less than or equal to n; or when the relation between the number of times of train cycle counting and the number of the ascending trains meets i ≠ m +1, and the index number Tupleft (i) of the element at the left end of the ascending train meets 0<When Tupleft (i) is less than or equal to n, calculating Y by formula (7)_SSCalculating Y by equation (8)_SU(ii) a Calculating Y by equation (10)_UU。

When the relation between the number of times of cycle counting of the train on the upper row and the train on the lower row meets i ═ m +1, and i ≠ 1, and the train on the upper right end element index number Tupright (i) meets Tupright (i-1)>When n is greater than n; or when the relation between the number of times of train cycle counting and the number of the ascending trains meets the condition that i is not equal to m +1, and the left end element index of the ascending trainNumber Tupleft (i) satisfies Tupleft (i)>When n is in the above range, Y is calculated from the above formulas (9) and (10)_UU。

In step S700, before the calculation of the index numbers of the left and right train elements in the down train cycle is started, the number of vehicles in the down train in the power supply section, the down power supply section processing flag vector, and the left and right train element index vectors are all set to zero, and the initial value of the number of times of the down train cycle count is set to 1.

Specifically, the setting may include: clearing a downlink train vector nRdowntrain in the power supply interval, wherein the dimensionality is (n +1) multiplied by 1, clearing a processing flag vector dealdown of the downlink power supply interval, the dimensionality is (n +1) multiplied by 1, and n is the number of traction substations; clearing a left index vector Tdownleft of the downlink train element, wherein the dimensionality is k multiplied by 1; and clearing a right index vector Tdown of the element of the downstream train, wherein the dimension is k multiplied by 1, and an initial value i of the number of times of the cycle count of the downstream train is defined to be 1, wherein k is the number of vehicles in the downstream row.

And after the setting is finished, calculating the index number of the element at the left end of the descending train, the index number of the element at the right end of the descending train and the node sub conductance matrix according to the relation between the real-time number of the circulation counting of the descending train and the number of the descending trains. The specific calculation includes:

circularly judging the relation between the real-time number i of the down train circular counting and the number k of the down train vehicles;

s710, firstly, for i>k + 1: if yes, the process proceeds to step S770, where Y is calculated_US(ii) a If the determination is no, that is, i is not greater than k +1, the process proceeds to step S720, and continues to determine that i is k + 1;

and S720, when the judgment result is negative, namely i is less than k +1, judging the jth power supply section to which the ith (namely the index number of the real-time descending train represented by the number i of the real-time train cycle counting) train belongs according to the formula (1), and setting the left index number and the right index number of the power supply section to which the descending train belongs, and the left end element index number and the right end element index number of the descending train according to the formula (5). Wherein x is the left index number of the jth power supply interval, y is the right index number of the jth power supply interval, and the numerical value is equal to 0, which represents that no traction substation exists on the left or right; and n is the total number of the traction substations.

When the initial downlinktrain number vector nrdownlinktrain (j) is equal to 0, the downlinktrain number vector nrdownlinktrain (j) at the time of actual counting satisfies nrdownlinktrain (j) of nrdownlinktrain (j) + 1; tdowleft (i) x, tdowright (i) y; when the initial nrdowntain (j) is not equal to 0, nrdowntain (j) ═ nrdowntain (j) + 1; tdowleft (i) n + i-1, tdowright (i) y, tdowright (i-1) n + i. Wherein n is the total number of the traction substation, i is the number of cycle counting, j is the number of the power supply section to which the downstream train belongs, x is the index number of the element at the left end of the downstream train, and y is the index number at the right end of the downstream train.

If yes, i.e., if i is k +1, the process proceeds to step S730, a processing flag vector dealup (j) of the power supply section to which the downstream train belongs and an index number j of the power supply section to which the downstream train belongs are determined, and if dealup (j) is equal to 0, and j is not equal to 1 and j is not equal to n +1, Y obtained in step S600 is determined by equation (6)_SSMaking a correction, after which dealup (j) is equal to 1; not correcting Y when dealup (j) is equal to 0, but j ≠ 1 and j ≠ n +1 are not satisfied_SS。

In formula (6) in this step, L_jFor the length of the jth power supply interval to which the downlink train belongs, other parameters are as described in step S600, and are not described herein again.

S740, determining that i is k +1, and if the determination is negative, that is, i is less than k +1, performing further determination and calculation in step S740; if the determination is yes, i is k +1, the process proceeds to step S750, and the determination is continued with i being 1, and if the determination is no, i < k +1 and i ≠ 1, the further determination and calculation in step S750 are performed. If yes, i is 1, the process proceeds to step S760.

The further judgment and calculation in step S740 includes judgment of the index number tdown (i) of the downstream train left end element:

when the judgment result is Tdownleft (i)0, no processing is carried out, namely no calculation is carried out;

when the judgment is 0<When Tdownleft (i) is less than or equal to n, calculating Y by the formula (7)_SSY is calculated by equation (11)_SD(ii) a Calculating Y by equation (13)_DD(ii) a Wherein a, b, c and d are index numbers, a is tdownleft (i), b is i, c is a, d is b, dis is the distance from the starting point of the ith downlink train to the starting point of the uplink train, and-tdownleft (i) the distance from the starting point of the traction substation to the starting point of the uplink train;

when the judgment is Tdownleft (i)>When n is in the above range, Y is calculated from the above formulas (12) and (13)_DD(ii) a Wherein a, b, c and d are index numbers, a is tdownleft (i) -n, b is i, d is b, dis is the distance from the starting point of the ith descending vehicle to the starting point of the ascending train-the distance from the starting point of the (i-1) th descending vehicle to the starting point of the ascending train;

the further determination and calculation in step S750 includes determining the downstream train right end element index number tdown (i):

when the Tdown (i-1) is judged to be 0, no processing is carried out, namely no calculation is carried out;

when the judgment is 0<When Tdown right (i-1) is not more than n, Y is calculated by the formula (7)_SSY is calculated by equation (11)_SD(ii) a Calculating Y by equation (13)_DD(ii) a Wherein, a is i, b is tdown (i-1), c is b, d is a, dis is the distance from the i-1 th downstream train to the upstream train starting point — the distance from the tdown (i-1) th traction substation to the upstream train starting point;

when the judgment is Tdown right (i-1)>When n is in the above range, Y is calculated from the above formulas (12) and (13)_DD(ii) a Wherein, a, b, c and d are index numbers, a is equal to i, b is equal to tdowright (i-1) -n, d is equal to a, dis is equal to the distance between the ith descending vehicle and the ascending train starting point-the distance between the (i-1) th descending vehicle and the ascending train starting point.

S760: the number of times i of the up train loop count is set to i +1, and the process returns to step S710 to perform loop determination and calculation.

S770: calculating Y_DSBy calculating Y_SDTransposing to obtain Y_DS。

That is, in step S700, when the relationship between the number of train cycles counted in the next train and the number of downstream trains satisfies i ═ k +1, the processing flag vector dealup (j) of the downstream train power supply section is equal to 0, and the index number of the downstream train power supply section satisfies j ≠ 1 and j ≠ n +1, Y is calculated by equation (6)_SSCorrecting;

when the relationship between the number of times of train cycle counting and the number of downstream trains satisfies i ═ k +1, i ≠ 1, and the downstream train right-end element index number tdown right (i) satisfies 0<When Tdown (i-1) is less than or equal to n; or the relation between the number of times of current train/queue cycle counting and the number of the downstream train is satisfied with i ≠ k +1, and the index number Tdownleft (i) of the downstream train left end element satisfies 0<When Tdownleft (i) is less than or equal to n, calculating Y by the formula (7)_SSY is calculated by equation (11)_SD(ii) a Calculating Y by equation (13)_DD。

When the relationship between the number of times of train cycle counting and the number of downstream trains satisfies i ═ k +1, i ≠ 1, and the downstream train right-end element index number Tdown (i) satisfies Tdown (i-1)>When n is greater than n; or the relation between the number of times of current train/queue cycle counting and the number of the downstream train is i ≠ k +1, and the index number Tdownleft (i) of the downstream train left end element satisfies Tdownleft (i)>When n is in the above range, Y is calculated from the above formulas (12) and (13)_DD。

By defining the left and right index vectors of the train, the circular calculation can be changed into single-layer calculation, and the calculation efficiency is greatly improved.

By solving the node conductance matrix of the whole subway direct-current power supply network in blocks, the symmetry of the node conductance can be utilized, and the obtaining form is as follows: sub-node conductance matrixes of the traction substation, the ascending train and the lower row workshop; the sub-node conductance matrixes of the ascending train and the ascending train, the sub-node conductance matrixes of the descending train and the node conductance matrixes of the ascending train, the descending train and the traction substation can avoid repeated calculation of the elements of the node conductance matrixes.

Examples

A method for rapidly calculating a node conductance matrix in subway load flow calculation is characterized by comprising the following steps:

1) inputting position and number information and network parameters:

the element position and number information comprises the distance between a traction substation and the starting point of an uplink train, the total number of the traction substations, the distance between the uplink train and the starting point of the uplink train, the number of uplink trains, the distance between a downlink train and the starting point of the uplink train and the number of downlink trains;

the network parameters comprise a contact network unit resistance, a walking steel rail unit resistance and a leakage unit resistance of the rail to the ground.

2) Numbering traction substations and dividing power supply intervals:

2-1) numbering traction substations from small to large according to the distance relative to the starting point of the ascending train, wherein n is the total number of the traction substations;

2-2) n traction substations generate n +1 power supply intervals, and the range of the ith power supply interval is as follows:

in the formula, sr represents the distance between the power supply section and the start point of the uplink train, ds represents the distance between the traction substation and the start point of the uplink train, n is the total number of the traction substations, and i is the index number of the traction substation or the power supply section.

And calculating the length of the 2 nd to the nth power supply intervals, wherein the calculation formula is L_i＝ds_i-ds_i-1,1<i is less than or equal to n (2). Wherein i is the index number of the traction substation, ds_iRepresenting the distance of the ith traction substation from the start of the upstream train.

3) Pre-calculating a node conductance matrix, wherein the matrix form is as follows:

in the formula, S represents a traction substation, U represents an up train, D represents a down train, and each subarray is assigned a value of 0.

Calculating Yss between the traction substation and the traction substation according to the network no-operation train, wherein the calculation formula is as follows:

in the formula i₁And j₁Numbering nodes, j₁∈i₁Represents directly with i₁Connected nodes, r₁Is the unit resistance of the contact line, L_tIs the length of the t-th power supply interval, t is the number of the power supply interval, n is the total number of the traction substation, r₂Is the unit resistance of the running rail, r₃The unit resistance of the leakage of the steel rail to the ground.

4) Numbering the trains:

4-1) sequencing the ascending trains from small to large according to the distance relative to the starting point of the ascending trains, and sequentially numbering the ascending trains as n + 1-n + m, wherein n is the total number of traction substations, and m is the number of the ascending trains;

4-2) sequencing the descending trains from small to large according to the distance relative to the starting point of the ascending train, and sequentially numbering the descending trains as n + m + 1-n + m + k, wherein n is the total number of traction substations, m is the number of the ascending trains, and k is the number of the descending trains;

5) calculating the index number of the left and right elements of the ascending train and the node sub conductance matrix:

5-1) clearing an uplink train number vector nRupTrain in the power supply interval, wherein the dimensionality is (n +1) multiplied by 1, and clearing an uplink power supply interval processing flag vector dealup, the dimensionality is (n +1) multiplied by 1, wherein n is the number of traction substations; clearing a left index vector Tupleft of the uplink train element, wherein the dimension is mx 1; clearing the right index vector Turight of the uplink train element, wherein the dimension is mx 1, m is the number of vehicles in the upper row and the lower row, and the cycle flag i is made to be 1;

5-2) judging that i is greater than m +1, wherein m is the number of the ascending trains, if the condition is met, entering 5-9), and if not, entering 5-3);

5-3) judging that i is equal to m +1, and entering 5-4 if the condition is met); otherwise, judging the jth power supply interval to which the ith ascending vehicle belongs according to the formula (1); and sets left and right indexes according to the formula (5)

In the formula, x is the left index number of the jth power supply interval, y is the right index number of the jth power supply interval, and if the number is equal to 0, the left side or the right side of the power supply interval is not provided with a traction substation; and n is the total number of the traction substations.

If nRupTrain (j) is equal to 0; nreptrain (j) ═ nreptrain (j) + 1; tupleft (i) x, tupright (i) y; otherwise, nreptrain (j) ═ nreptrain (j) + 1; n is the total number of traction substations, i is a cycle mark, j is a power supply interval number, x is an uplink train left end element index number, and y is an uplink train right end index number.

5-4) if dealup (j) is equal to 0, and j is not equal to 1 and j is not equal to n +1, n is the total number of traction substations, correcting Yss, dealup (j) is equal to 1; otherwise, not correcting the Yss, wherein the correction formula is as follows:

r₁contact net unit resistance r₂Is the unit resistance of the running rail, r₃For rail leakage to ground, unit resistance, L_jThe length of the jth power supply interval.

5-5) if i is equal to m +1, turning to 5-6), otherwise, judging as follows:

if tupleft (i) is 0, do not process;

if 0< tupleft (i) <n, dis ═ i distance between the i-th ascending train and the ascending train starting point-the i-th traction substation and the ascending train starting point distance, a ═ tupleft (i), b ═ i, c ═ a, d ═ b; calculating (7), (8) and (10);

if tupleft (i) > n, dis is the distance between the ith ascending vehicle and the starting point of the ascending train-the distance between the (i-1) th ascending vehicle and the starting point of the ascending train; a ═ tupleft (i) -n, b ═ i, d ═ b; calculating (9) and (10);

5-6) if i is equal to 1, turning to 5-7), otherwise, judging as follows:

if Tupright (i-1) ═ 0, no treatment;

if 0< tuprint (i-1) < ═ n, dis ═ i-1 the distance between the starting point of the uplink relative to the i-th traction substation-the distance between the starting point of the uplink relative to the i-th traction substation; a ═ i, b ═ Tupright (i-1), c ═ b, d ═ a; calculating (7), (8) and (10);

if Turight (i-1) > n, dis is the distance between the ith ascending vehicle and the starting point of the ascending train-the distance between the (i-1) th ascending vehicle and the starting point of the ascending train; a ═ i, b ═ Tupright (i-1) -n, d ═ a; calculating (9) and (10);

(7) in the formula (1) to (10), r₁Contact net unit resistance r₂Is the unit resistance of the running rail, r₃For the unit resistance of the rail to the ground leakage, dis is the distance between adjacent elements, and a, b, c and d are index numbers.

5-7) let i ═ i +1, go to 5-2).

5-8) calculating Y_USA value of Y_SUThe transposing of (1).

6) Calculating the index number of the left and right elements of the downlink train and the node sub-conductance matrix:

6-1) clearing a downlink train vector nRdowntrain in the power supply interval, wherein the dimensionality is (n +1) multiplied by 1, clearing a processing flag vector dealdown of the downlink power supply interval, the dimensionality is (n +1) multiplied by 1, and n is the number of traction substations; clearing a left index vector Tdownleft of the downlink train element, wherein the dimensionality is k multiplied by 1; clearing a right index vector Tdown of the downstream train element, wherein the dimensionality is k multiplied by 1, and k is the number of downstream trains and the circulating mark i is 1;

6-2) judging that i is greater than k +1, wherein k is the number of the downstream trains, if the condition is met, entering 6-9), and if not, entering 6-3);

6-3) judging that i is equal to k +1, wherein k is the number of vehicles going downwards, and if the number of vehicles going downwards meets the requirement, entering 6-4); otherwise, judging the jth power supply interval to which the ith ascending vehicle belongs according to the formula (1); and setting a left index and a right index according to the formula (5):

if nRdowntrain (j) is equal to 0; nrdowntain (j) ═ nrdowntain (j) + 1; tdowleft (i) x, tdowright (i) y; otherwise, nrdowntain (j) ═ nrdowntain (j) + 1; tdowleft (i) n + i-1, tdowright (i) y, tdowright (i-1) n + i. n is the total number of the traction substation, i is a cycle mark, j is a power supply interval number, x is a downlink train left end index number, and y is a downlink train right end index number.

6-4) if dealdown (j) is equal to 0, and j is not equal to 1 and j is not equal to n +1, n is the total number of the traction substation, correcting Y_SSDealdown (j) 1; otherwise Y_SSThe correction formula is shown as (6) without correction:

6-5) if i is equal to k +1, k is the number of vehicles going downwards, turning to 5-6), otherwise, judging as follows:

if tdwonleft (i) ═ 0;

if 0< tdownleft (i) < ═ n, dis ═ i distance between the i-th descending vehicle and the starting point of the ascending train, the t-downteft (i) distance between the traction substation and the starting point of the ascending train, a ═ tdownleft (i), b ═ i, c ═ a, d ═ b; calculating (7), (11) and (13);

if Tdownleft (i) > n, dis is the distance from the starting point of the ith descending vehicle to the starting point of the ascending train-the distance from the starting point of the (i-1) th descending vehicle to the starting point of the ascending train; a ═ tdownleft (i) -n, b ═ i, d ═ b; calculating (12) and (13);

6-6) if i is equal to 1, turning to 6-7), otherwise, judging as follows:

if Tdown right (i-1) ═ 0, do not process;

if 0< Tdown (i-1) n, dis is the distance between the ith uplink vehicle and the starting point of the uplink train-the distance between the Tdown (i-1) traction substation and the starting point of the uplink train; a, b, t down (i-1), c, d, a; calculating (7), (11) and (13);

if Tdownright (i-1) > n, dis is the distance between the ith ascending vehicle and the ascending train starting point-the distance between the (i-1) th ascending vehicle and the ascending train starting point; a, b, Tdownright (i-1) -n, d, a; calculating (12) and (13);

(11) in the formulae (1) to (13), r₁Contact net unit resistance r₂Is the unit resistance of the running rail, r₃For the unit resistance of the rail to the ground leakage, dis is the distance between adjacent elements, and a, b, c and d are index numbers.

6-7) let i ═ i +1, go to 6-2).

6-8) calculating Y_DSA value of Y_SDThe transposing of (1).

7) An output node conductance matrix.

The method for rapidly calculating the node conductance matrix in the subway load flow calculation has the following effects: 1) only two nodes are considered in a network of a traction substation, an uplink train and a downlink train, wherein the nodes in the uplink network and the downlink network of the same traction substation are the same in number, and the memory space of a node conductance matrix is reduced; 2) the method for solving the subarray of the node conductance matrix in a blocking mode specifically comprises the following steps: the node conductance matrixes of the ascending train and the ascending train, the node conductance matrixes of the descending train and the descending train, the traction substation and the traction substation, the ascending train and the lower row and column workshop. By utilizing the symmetry of the node conductance, a node conductance matrix of an uplink train, a downlink train and a traction substation is obtained, and the repeated calculation of elements is avoided; 3) considering that the contact network nodes or the steel rail nodes of the train are only connected with two nodes, defining left and right index vectors of the train, and directly finishing the calculation of the electric conductance sub-matrix of the nodes, wherein the double-layer circulation of the traditional method is reduced into a single layer; 4) and (3) solving a node conductance matrix between the traction substation and the traction substation in advance, and dynamically modifying the mutual conductance of the node conductance matrix between the traction substation and the traction substation each time in a train insertion mode. In the algorithm, in the multiple power flow calculation, the node conductance matrixes of the traction substation and the traction substation are calculated only once, the characteristic that a train continuously runs for a period of time in the same power supply interval can be considered, the repeated calculation of the node conductance matrix elements is avoided by utilizing the previous node conductance matrix between the traction substation and the traction substation, and the calculation efficiency of the algorithm is improved.

It should be noted that the method of the embodiment of the present invention may be executed by a single device, such as a computer or a server. The method of the embodiment can also be applied to a distributed scene and completed by the mutual cooperation of a plurality of devices. In the case of such a distributed scenario, one of the multiple devices may only perform one or more steps of the method according to the embodiment of the present invention, and the multiple devices interact with each other to complete the method.

The foregoing description has been directed to specific embodiments of this disclosure. Other embodiments are within the scope of the following claims. In some cases, the actions or steps recited in the claims may be performed in a different order than in the embodiments and still achieve desirable results. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In some embodiments, multitasking and parallel processing may also be possible or may be advantageous.

Those of ordinary skill in the art will understand that: the discussion of any embodiment above is meant to be exemplary only, and is not intended to intimate that the scope of the disclosure, including the claims, is limited to these examples; within the idea of the invention, also features in the above embodiments or in different embodiments may be combined, steps may be implemented in any order, and there are many other variations of the different aspects of the invention as described above, which are not provided in detail for the sake of brevity.

In addition, well known power/ground connections to Integrated Circuit (IC) chips and other components may or may not be shown within the provided figures for simplicity of illustration and discussion, and so as not to obscure the invention. Furthermore, devices may be shown in block diagram form in order to avoid obscuring the invention, and also in view of the fact that specifics with respect to implementation of such block diagram devices are highly dependent upon the platform within which the present invention is to be implemented (i.e., specifics should be well within purview of one skilled in the art). Where specific details (e.g., circuits) are set forth in order to describe example embodiments of the invention, it should be apparent to one skilled in the art that the invention can be practiced without, or with variation of, these specific details. Accordingly, the description is to be regarded as illustrative instead of restrictive.

While the present invention has been described in conjunction with specific embodiments thereof, many alternatives, modifications, and variations of these embodiments will be apparent to those of ordinary skill in the art in light of the foregoing description. For example, other memory architectures (e.g., dynamic ram (dram)) may use the discussed embodiments.

The embodiments of the invention are intended to embrace all such alternatives, modifications and variances that fall within the broad scope of the appended claims. Therefore, any omissions, modifications, substitutions, improvements and the like that may be made without departing from the spirit and principles of the invention are intended to be included within the scope of the invention.

Claims

1. A method for rapidly calculating a node conductance matrix in subway load flow calculation is characterized by comprising the following steps:

s400, calculating Y according to the non-train operation_SS；

s600, clearing an uplink train number vector, an uplink power supply interval processing flag vector, an uplink train left end element index number and an uplink train right end element index number in a power supply interval, and setting an initial value of the uplink train cycle count number; and according to the number of times of the circulation counting of the uplink train andcircularly calculating the index number of the element at the left end of the ascending train and the index number of the element at the right end of the ascending train according to the relation of the number of the vehicles in the upper row and the lower row; and to Y_SSMaking a correction or a correction of Y_SSCalculating Y_SUAnd Y_UUOr calculating Y_UUOr calculating Y_US；

S700, clearing a downlink train number vector, a downlink power supply interval processing flag vector, a downlink train left end element index number and a downlink train right end element index number in a power supply interval, and setting an initial value of the downlink train cycle count number; circularly calculating the index number of the element at the left end of the descending train and the index number of the element at the right end of the descending train according to the relation between the number of times of circularly counting the descending trains and the number of the descending trains; and to Y_SSMake a correction of Y_SDAnd Y_DDOr calculating Y_DDOr calculating Y_DS。

2. The method for rapidly calculating the node conductance matrix in the subway flow calculation as claimed in claim 1, wherein the step S400 comprises,

by passing

3. The method for rapidly calculating the node conductance matrix in the subway load flow calculation as claimed in claim 2, wherein Y is_SSThe correcting comprises the following steps: when the relation between the number of times of the cycle counting of the upper train and the lower train and the number of the upper train meets i ═ m +1, the processing mark vector dealup (j) of the power supply interval of the upper train is equal to 0, and the cable of the power supply interval of the upper train is connected with the processing mark vector dealup (j) of the power supply interval of the upper trainWhen the index number is j ≠ 1 and j ≠ n +1, or when the relationship between the number of train cycle counts and the number of downstream trains is i ═ k +1, the processing flag vector dealdown (j) of the downstream train power supply section is equal to 0, and the index number of the power supply section is j ≠ 1 and j ≠ n +1,

by passing

4. The method for rapidly calculating the node conductance matrix in the subway load flow calculation as claimed in claim 3, wherein in step S600, Y is corrected_SSCalculating Y_SUAnd Y_UUWhen the relation between the number of times of cycle counting of the train and the number of the ascending trains meets i ═ m +1, i ≠ 1, and the right-end element index number Turight (i) of the ascending train meets 0<When Turight (i-1) is less than or equal to n; or when the relation between the number of times of train cycle counting and the number of the ascending trains meets i ≠ m +1, and the index number Tupleft (i) of the element at the left end of the ascending train meets 0<When the content of Tupleft (i) is less than or equal to n,

by passing

Correcting the Yss;

by passing

Calculating Y_SU；

By passing

Calculating Y_UU，

Wherein a, b, c and d are index numbers, a is i, b is Tupright (i-1) The distance between c and d is a, dis is the distance between the i-th ascending train and the starting point of the ascending train-the starting point of the first traction substation-the starting point of the ascending train, or a is Tupleft (i), b is i, c is a, d is b, dis is the distance between the i-th ascending train and the starting point of the ascending train-Tupleft (i)₁) And the distance between each traction substation and the starting point of the ascending train.

5. The method for rapidly calculating the node conductance matrix in the subway load flow calculation as claimed in claim 4, wherein in step S600, Y is calculated_UUIncludes that when the relation between the number of times of cycle counting of the train and the number of the ascending trains meets i ═ m +1, i ≠ 1, and the right-end element index number Tupright (i) of the ascending train meets Tupright (i-1)>When n is greater than n; or when the relation between the number of times of train cycle counting and the number of the ascending trains meets i ≠ m +1, and the index number Tupleft (i) of the ascending train left end element meets Tupleft (i)>When the number n is greater than the predetermined value,

passing through type

Calculating Y_UUWherein, a ═ i, b ═ Tupright (i-1) -n, d ═ a, dis ═ i ascending vehicle relative ascending train starting point distance-i-1) th ascending vehicle relative ascending train starting point distance or a ═ tupleft (i) -n, b ═ i, d ═ b, dis ═ i ascending vehicle relative ascending train starting point distance-i-1 th ascending vehicle relative ascending train starting point distance.

6. The method for rapidly calculating the node conductance matrix in the subway load flow calculation as claimed in claim 4, wherein in step S600, Y is calculated_USIncluding that the relation between the number of times of circularly counting the train on line and the number of the train on line satisfies i>m +1, calculating Y_SUIs transposed to obtain Y_US。

7. The method for rapidly calculating the node conductance matrix in the subway load flow calculation as claimed in claim 2, wherein in step S700, Y is calculated_SS、Y_SDAnd Y_DDIncludes that the relation between the number of times of train cycle counting and the number of downstream trains satisfies i ═ k +1, i ≠ 1, and the downstream train right-end element index number Tdown right (i) satisfies 0<When Tdown (i-1) is less than or equal to n; or when the relation between the number of train cycle counting and the number of the downlink train is i ≠ i +1, and the index number Tdownleft (i) of the element at the left end of the downlink train satisfies 0<When Tdownleft (i) is less than or equal to n, the medicine is passed

Calculating Yss;

by passing

Calculating Y_SD；

By passing

Calculating Y_DD，

8. The method for rapidly calculating the node conductance matrix in the subway load flow calculation as claimed in claim 7, wherein calculating Y in step S700_DDIncludes that the relation between the number of times of train cycle counting and the number of downstream trains satisfies i ═ k +1, i ≠ 1, and the downstream train right-end element index number Tdown (i) satisfies Tdown (i-1)>When n is greater than n; or the relation between the number of times of current train/queue cycle counting and the number of the downstream train is i ≠ k +1, and the index number Tdownleft (i) of the downstream train left end element satisfies Tdownleft (i)>When the number n is greater than the predetermined value,

by passing

9. The method for rapidly calculating the node conductance matrix in the subway load flow calculation as claimed in claim 7, wherein calculating Y in step S700_DSIncluding that the relation between the number of times of current train-train cycle counting and the number of downstream trains meets i>k +1, calculating Y_SDIs transposed to obtain Y_DS。

10. The method for rapidly calculating the node conductance matrix in the subway load flow calculation according to claim 1, wherein said calculating the index number of the left end element of the ascending train and the index number of the right end element of the ascending train or said calculating the index number of the left end element of the descending train and the index number of the right end element of the descending train comprises: