CN105083322A - Subway train collision early warning method - Google Patents
Subway train collision early warning method Download PDFInfo
- Publication number
- CN105083322A CN105083322A CN201510150113.4A CN201510150113A CN105083322A CN 105083322 A CN105083322 A CN 105083322A CN 201510150113 A CN201510150113 A CN 201510150113A CN 105083322 A CN105083322 A CN 105083322A
- Authority
- CN
- China
- Prior art keywords
- train
- conflict
- track
- time
- subway
- 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
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000005070 sampling Methods 0.000 claims abstract description 10
- 238000013507 mapping Methods 0.000 claims abstract description 7
- 230000035945 sensitivity Effects 0.000 claims abstract description 4
- 230000008569 process Effects 0.000 claims description 32
- 238000001514 detection method Methods 0.000 claims description 21
- 239000011159 matrix material Substances 0.000 claims description 21
- 230000004907 flux Effects 0.000 claims description 18
- 238000012544 monitoring process Methods 0.000 claims description 15
- 230000008859 change Effects 0.000 claims description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 230000007246 mechanism Effects 0.000 claims description 5
- 206010027476 Metastases Diseases 0.000 claims description 3
- 238000005267 amalgamation Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000010924 continuous production Methods 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 claims description 3
- 230000009365 direct transmission Effects 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 238000005315 distribution function Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000009401 metastasis Effects 0.000 claims description 3
- 238000004451 qualitative analysis Methods 0.000 claims description 3
- 238000004445 quantitative analysis Methods 0.000 claims description 3
- 238000012549 training Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000003064 k means clustering Methods 0.000 claims 1
- 238000012545 processing Methods 0.000 claims 1
- 238000005096 rolling process Methods 0.000 abstract description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 206010033307 Overweight Diseases 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 239000002355 dual-layer Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005206 flow analysis Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 235000020825 overweight Nutrition 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L27/00—Central railway traffic control systems; Trackside control; Communication systems specially adapted therefor
- B61L27/10—Operations, e.g. scheduling or time tables
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L23/00—Control, warning or like safety means along the route or between vehicles or trains
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L27/00—Central railway traffic control systems; Trackside control; Communication systems specially adapted therefor
- B61L27/20—Trackside control of safe travel of vehicle or train, e.g. braking curve calculation
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L27/00—Central railway traffic control systems; Trackside control; Communication systems specially adapted therefor
- B61L27/40—Handling position reports or trackside vehicle data
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L27/00—Central railway traffic control systems; Trackside control; Communication systems specially adapted therefor
- B61L27/60—Testing or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Train Traffic Observation, Control, And Security (AREA)
Abstract
The invention relates to a subway train collision early warning method, which comprises the following steps: firstly, generating a topological structure chart of a rail transit network according to planned operation parameters of each train; then, based on the topological structure chart, the controllability and the sensitivity of the train flow are analyzed; generating a conflict-free running track of the multiple trains according to the planned running parameters of each train; predicting the advancing position of the train at a certain future moment at each sampling moment based on the current running state and the historical position observation sequence of the train, establishing an observer from the continuous dynamic state of the train to the discrete conflict logic, and mapping the continuous dynamic state into a conflict state expressed by a discrete observation value; when the system possibly violates the traffic control rule, the hybrid dynamic behavior of the subway traffic hybrid system is monitored, and warning information is provided for the control center. The invention predicts the subway train track in real time by rolling, effectively warns train conflict and improves the safety of subway traffic.
Description
Technical field
The present invention relates to a kind of subway train conflict method for early warning, particularly relate to a kind of subway train based on Robust Strategies conflict method for early warning.
Background technology
Along with the expanding day of China's big and medium-sized cities scale, Traffic Systems is faced with the increasing pressure, and greatly developing Rail Transit System becomes the important means solving urban traffic congestion.Country's Eleventh Five-Year Plan outline is pointed out, big city with good conditionsi and group of cities area will using track traffic as Priority settings.China is just experiencing a unprecedented track traffic development peak load conditions, and some cities have turned to the construction of net by the construction of line, and urban mass transit network is progressively formed.At Rail traffic network and the intensive complex region of train flow, still the train interval dispensing mode adopting train operation plan to combine based on subjective experience demonstrates its lag gradually, be in particular in: the formulation of (1) train operation plan timetable also reckons without the impact of various enchancement factor, easily cause flow of traffic tactics to manage crowded, reduce the safety that traffic system is run; (2) train scheduling active side overweights the personal distance kept between single train, not yet rises to the macroscopic aspect of train flow being carried out to strategy management; (3) train allocation process depends on the subjective experience of a line dispatcher more, and the selection randomness of allocating opportunity is comparatively large, lacks scientific theory and supports; (4) dispatcher's less impact considering external interference factor of allotment means of using, robustness and the availability of train programs are poor.For ensureing the safe operation of subway transportation, implement the emphasis that actv. conflict early warning just becomes subway transportation control work.Implement the emphasis that actv. subway conflict early warning just becomes subway transportation control work.
The discussion object spininess of existing documents and materials to long-distance railway transportation, and still lacks system for the Scientific Regulation scheme of the city underground traffic system under large discharge, high density and closely-spaced condition of service.Train Coordinated Control Scheme under complicated road network condition of service needs calculate the running state of single vehicles in traffic network in region and optimize on strategic level, and implements collaborative planning to the flow of traffic be made up of multiple train; Pre-tactical level solves congestion problems by actv. monitoring mechanism adjustment traffic network top subregional critical operational parameters, but at present to the prediction of subway train track and see train conflict early warning all less than scheme comparatively accurately.
Summary of the invention
The technical problem to be solved in the present invention is to provide a kind of robustness and availability good subway train conflict method for early warning, and the precision of prediction of the method to subway train is higher, the accuracy of subway train conflict early warning and ageing all better.
The technical scheme realizing the object of the invention is to provide a kind of subway train conflict method for early warning, comprises the steps:
Steps A, plan operational factor according to each train, the topology diagram of generator orbital traffic network;
Step B, topology diagram based on the Rail traffic network constructed by steps A, analyze controllability and sensivity two class feature of train flow;
Step C, plan operational factor according to each train, on the basis building Modeling Method for Train Dynamics, set up train operation conflict according to train operation conflict Coupling point and allocate model in advance, generate multiple row car Lothrus apterus running orbit;
Step D, at each sampling instant t, the running state current based on train and historical position observation sequence, predict the advanced positions in train certain moment following; Its detailed process is as follows:
Step D1, train track data pretreatment, with train at the stop position of originating station for the origin of coordinates, in each sampling instant, according to the train original discrete two-dimensional position sequence x=[x obtained
1, x
2..., x
n] and y=[y
1, y
2..., y
n], adopt first order difference method to carry out processing new train discrete location sequence Δ x=[the Δ x of acquisition to it
1, Δ x
2..., Δ x
n-1] and Δ y=[Δ y
1, Δ y
2..., Δ y
n-1], wherein Δ x
i=x
i+1-x
i, Δ y
i=y
i+1-y
i(i=1,2 ..., n-1);
Step D2, to train track data cluster, to new train discrete two-dimensional position sequence Δ x and Δ y after process, by setting cluster number M ', K-means clustering algorithm is adopted to carry out cluster to it respectively;
Step D3, HMM is utilized to carry out parameter training to the train track data after cluster, by the train operation track data Δ x after process and Δ y being considered as the aobvious observed value of hidden Markov models, upgrade period τ ' by setting hidden state number N ' and parameter, according to nearest T ' individual position detection value and adopt B-W algorithm to roll to obtain up-to-date HMM parameter lambda '; Specifically: because obtained train track sets data length is dynamic change, in order to the state variation of real-time tracking train track, be necessary at initial track HMM parameter lambda '=(π, A, B) basis is readjusted it, to infer the position of train in certain moment following more accurately; Every period τ ', the individual observed value (o of the T ' according to up-to-date acquisition
1, o
2..., o
t ') to track HMM parameter lambda '=(π, A, B) reappraise;
Step D4, foundation HMM parameter, adopt the hidden state q corresponding to Viterbi algorithm acquisition current time observed value;
Step D5, every the period
, according to the HMM parameter lambda of up-to-date acquisition '=(π, A, B) and nearest H conception of history measured value (o
1, o
2..., o
h), based on the hidden state q of train current time, at moment t, by setting prediction time domain h ', obtain the position prediction value O of future time period train;
Step e, setting up the dynamic continuously observer to discrete conflict logic from train, is the conflict situation that discrete observation value is expressed by the continuous dynamic mapping of subway transportation system; When system likely violates traffic control rule, to the Hybrid dynamics behavior implementing monitoring of subway transportation hybrid system, for subway transportation control center provides warning information timely.
Further, the detailed process of steps A is as follows:
Steps A 1, extract the site information of stopping each train travelling process from the data bank of subway transportation control center;
Steps A 2, according to positive and negative two service directions, the site information that each train is stopped to be classified, and the same site on same service direction is merged;
Steps A 3, according to website amalgamation result, according to website space layout's form with straight line connect before and after multiple website.
Further, the detailed process of step B is as follows:
Step B1, the Traffic flux detection model built in single subsegment; Its detailed process is as follows:
Step B1.1, introducing state variable Ψ, input variable u and output variable Ω, wherein Ψ represents the train quantity of certain moment existence on connected section between website, it comprises single channel section and Multiple Sections two type, u represents the Operation Measures that track traffic dispatcher implements for certain section, as adjust train speed or change train in the station time etc., Ω represents the train quantity that certain period section is left;
Step B1.2, by by time discretization, set up shape as Ψ (t+ Ψ t)=A
1Ψ (t)+B
1u (t) and Ω (t)=C
1Ψ (t)+D
1discrete time Traffic flux detection model in the single subsegment of u (t), wherein Δ t represents the sampling interval, and Ψ (t) represents the state vector of t, A
1, B
1, C
1and D
1represent the state-transition matrix of t, input matrix, output calculation matrix and direct transmission matrix respectively;
Step B2, the Traffic flux detection model built in many subsegments; Its detailed process is as follows:
Step B2.1, according to circuit space layout form and train flow historical statistical data, obtain the flow proportional parameter beta in cross link each subsegment;
Step B2.2, according to the discrete time Traffic flux detection model in flow proportional parameter and single subsegment, build shape as Ψ (t+ Δ t)=A
1Ψ (t)+B
1u (t) and Ω (t)=C
1Ψ (t)+D
1discrete time Traffic flux detection model in many subsegments of u (t);
Step B3, controllable factor matrix [B according to Controlling model
1, A
1b
1..., A
1 n-1b
1] order and the relation of numerical value n, its controllability of qualitative analysis, according to the sensitivity coefficient matrix [C of Controlling model
1(zI-A
1)
-1b
1+ D
1], its input and output sensivity of quantitative analysis, wherein n represents the dimension of state vector, I representation unit matrix, and z represents the element factor changed original discrete time Traffic flux detection model.
Further, the detailed process of step C is as follows:
Step C1, train status transfer modeling, train shows as the switching at runtime process between website along the process that track traffic road network runs, arrange according to the website in train operation plan, set up the Petri network model of single train switched and transferred between different website: E=(g, G, Pre, Post, m) be train section metastasis model, wherein g represents each sub-section between website, G represents the change-over point of train running speed state parameter, Pre and Post represents that front and back between each sub-section and website are to annexation respectively
represent the operation section residing for train, wherein m represents model identification, Z
+represent Positive Integer Set:
The modeling of step C2, train full operation profile hybrid system, the operation of train between website is considered as continuous process, from the stressed situation of train, according to the kinetics equation of energy model derivation train in the different operation phase, in conjunction with external interference factor, set up about train in a certain operation phase speed v
gmapping function v
g=λ (T
1, T
2, H, R, α), wherein T
1, T
2, H, R and α represent tractive force of train, braking force of train, train resistance, train gravity and train status random fluctuation parameter respectively;
Step C3, adopt the mode mixing emulation to infer to solve train track, by by time subdivision, the distance of utilization state continually varying characteristic Recursive Solution any time train in a certain operation phase apart from initial rest position point,
wherein J
0for initial time train is apart from the voyage of initial rest position point, Δ τ is the numerical value of time window, J(τ) for τ moment train is apart from the distance of initial rest position point, can infers thus and obtain single vehicles track;
Step C4, train, in station time probability distribution function modeling, for specific run circuit, by transferring the dwell time data of train at each station, obtain the dwell time probability distribution of different circuit different website condition Train;
The Lothrus apterus robust track allotment of step C5, multiple row car coupling, the time of conflict point is reached in advance according to each train, pass through Time segments division, at each sampling instant t, under the prerequisite incorporating random factor, according to scheduling rule, robust secondary planning is implemented to the train track not meeting personal distance requirement near conflict point.
Further, in step D, the value of cluster number M ' is 4, and the value of hidden state number N ' is 3, and it is 30 seconds that parameter upgrades period τ ', and T ' is 10,
be 30 seconds, H is 10, and prediction time domain h ' is 300 seconds.
Further, the specific implementation process of step e is as follows:
Step e 1, construct conflict hypersurface collection of functions based on regulation rule: set up hypersurface collection of functions in order to reflect the contention situation of system, wherein, continuous function h relevant to single train in conflict hypersurface
ibe I type hypersurface, the continuous function h relevant to two trains
iIit is II type hypersurface;
Step e 2, set up by train continuous state to the observer of discrete conflict situation, build the safety rule collection d that need meet when train runs in traffic network
ii(t)>=d
min, wherein d
ijt () represents the actual interval of train i and train j in t, d
minrepresent the minimum safety interval between train;
Step e 3, based on person machine system theoretical and complication system hierarchical control principle, according to train operation pattern, build people at the real-time monitoring mechanism of the train of loop, the operation of guarantee system is in safe reachable set, design the discrete monitor from conflict to conflict Resolution means, when the discrete observation vector of observer shows that safety rule rally is breached, send corresponding warning information to subway transportation control center at once.
The present invention has positive effect: (1) subway train conflict of the present invention method for early warning is under the prerequisite meeting track traffic control personal distance, and based on the real-time position information of train, maintenance data excavates means and dynamically infers train track; According to track traffic regulation rule, alarm is implemented to the conflict that may occur, better to the early warning effect of conflict, can the track of train be predicted effectively, accurately and real-time and predict train conflict, effectively improve the safety of subway transportation.
(2) the present invention is based on controllability and the susceptivity analysis result of Rail traffic network topological structure, can be the early warning of subway transportation stream and scientific basis is provided, overcome the randomness of conventional early warning scheme selection.
(3) the present invention is based on the scene monitoring mechanism of constructed " people is at loop ", effecting reaction can be made in time alternately to train inside continuous variable and the frequent of external discrete event, overcome the shortcoming of conventional open loop monitored off-line scheme.
(4) the present invention is based on constructed train operation track rolling forecast scheme, all kinds of disturbing factors in train real time execution can be incorporated in time, improve the accuracy of train trajectory predictions, overcome the shortcoming that Conventional Off-line prediction scheme accuracy rate is not high.
Accompanying drawing explanation
Fig. 1 is train flow analysis on Operating figure;
Fig. 2 is Lothrus apterus 3D robust track supposition figure;
Fig. 3 is that train operation state mixes monitoring figure.
Detailed description of the invention
(embodiment 1)
The flow-optimized control system of a kind of subway transportation, comprise wire topologies generation module, data transmission module, car-mounted terminal module, control terminal module and track monitoring module, track monitoring module is collected the status information of train and is supplied to control terminal module.
Described control terminal module comprises following submodule:
Lothrus apterus Track Pick-up module before train operation: according to the Train operation plan table time of running, first set up Modeling Method for Train Dynamics, then sets up train operation conflict according to train operation conflict Coupling point and allocates model in advance, finally generate Lothrus apterus train operation track.
Train operation Track Pick-up a middle or short term module: the train real time status information provided according to track monitoring module, utilizes data mining model, infers the running orbit of train in future time period.
Train operation situation monitoring module: at each sampling instant t, based on the track estimation result of train, when likely occurring violating the situation of safety rule when between train, provides warning information to its dynamic behaviour implementing monitoring and for control terminal.
Train collision avoidance track optimizing module: when train operation situation monitoring module sends warning information, meet train physical property, region hold stream constraint and track traffic scheduling rule prerequisite under, by setting optimizing index function, adopt Adaptive Control Theory method to carry out robust dual layer resist by control terminal module to train operation track, and by data transmission module, program results is transferred to the execution of car-mounted terminal module.Train collision avoidance track optimizing module comprises internal layer planning and outer planning two class planning process.
Apply the subway train conflict method for early warning of the flow-optimized control system of above-mentioned subway transportation, comprise the following steps:
Steps A, plan operational factor according to each train, the topology diagram of generator orbital traffic network; Its detailed process is as follows:
Steps A 1, extract the site information of stopping each train travelling process from the data bank of subway transportation control center;
Steps A 2, according to positive and negative two service directions, the site information that each train is stopped to be classified, and the same site on same service direction is merged;
Steps A 3, according to website amalgamation result, according to website space layout's form with straight line connect before and after multiple website.
Step B, topology diagram based on the Rail traffic network constructed by steps A, analyze controllability and sensivity two class feature of train flow; Its detailed process is as follows:
Step B1, see Fig. 1, build the Traffic flux detection model in single subsegment; Its detailed process is as follows:
Step B1.1, introducing state variable Ψ, input variable u and output variable Ω, wherein Ψ represents the train quantity of certain moment existence on connected section between website, it comprises single channel section and Multiple Sections two type, u represents the Operation Measures that track traffic dispatcher implements for certain section, as adjust train speed or change train in the station time etc., Ω represents the train quantity that certain period section is left;
Step B1.2, by by time discretization, set up shape as Ψ (t+ Δ t)=A
1Ψ (t)+B
1u (t) and Ω (t)=C
1Ψ (t)+D
1discrete time Traffic flux detection model in the single subsegment of u (t), wherein Δ t represents the sampling interval, and Ψ (t) represents the state vector of t, A
1, B
1, C
1and D
1represent the state-transition matrix of t, input matrix, output calculation matrix and direct transmission matrix respectively;
Step B2, the Traffic flux detection model built in many subsegments; Its detailed process is as follows:
Step B2.1, according to circuit space layout form and train flow historical statistical data, obtain the flow proportional parameter beta in cross link each subsegment;
Step B2.2, according to the discrete time Traffic flux detection model in flow proportional parameter and single subsegment, build shape as Ψ (t+ Δ t)=A
1Ψ (t)+B
1u (t) and Ω (t)=C
1Ψ (t)+D
1discrete time Traffic flux detection model in many subsegments of u (t);
Step B3, controllable factor matrix [B according to Controlling model
1, A
1b
1..., A
1 n-1b
1] order and the relation of numerical value n, its controllability of qualitative analysis, according to the sensitivity coefficient matrix [C of Controlling model
1(zI-A
1)
-1b
1+ D
1], its input and output sensivity of quantitative analysis, wherein n represents the dimension of state vector, I representation unit matrix, and z represents the element factor changed original discrete time Traffic flux detection model;
Step C, see Fig. 2, according to the plan operational factor of each train, on the basis building Modeling Method for Train Dynamics, set up train operation conflict according to train operation conflict Coupling point and allocate model in advance, generate multiple row car Lothrus apterus running orbit; Its detailed process is as follows:
Step C1, train status transfer modeling, train shows as the switching at runtime process between website along the process that track traffic road network runs, arrange according to the website in train operation plan, set up the Petri network model of single train switched and transferred between different website: E=(g, G, Pre, Post, m) be train section metastasis model, wherein g represents each sub-section between website, G represents the change-over point of train running speed state parameter, Pre and Post represents that front and back between each sub-section and website are to annexation respectively
represent the operation section residing for train, wherein m represents model identification, Z
+represent Positive Integer Set;
The modeling of step C2, train full operation profile hybrid system, the operation of train between website is considered as continuous process, from the stressed situation of train, according to the kinetics equation of energy model derivation train in the different operation phase, in conjunction with external interference factor, set up about train in a certain operation phase speed v
gmapping function v
g=λ (T
1, T
2, H, R, α), wherein T
1, T
2, H, R and α represent tractive force of train, braking force of train, train resistance, train gravity and train status random fluctuation parameter respectively;
Step C3, adopt the mode mixing emulation to infer to solve train track, by by time subdivision, the distance of utilization state continually varying characteristic Recursive Solution any time train in a certain operation phase apart from initial rest position point,
wherein J
0for initial time train is apart from the voyage of initial rest position point, Δ τ is the numerical value of time window, and J (τ), for τ moment train is apart from the distance of initial rest position point, can infers thus and obtain single vehicles track;
Step C4, train, in station time probability distribution function modeling, for specific run circuit, by transferring the dwell time data of train at each station, obtain the dwell time probability distribution of different circuit different website condition Train;
The Lothrus apterus robust track allotment of step C5, multiple row car coupling, the time of conflict point is reached in advance according to each train, pass through Time segments division, at each sampling instant t, under the prerequisite incorporating random factor, according to scheduling rule, robust secondary planning is implemented to the train track not meeting personal distance requirement near conflict point.
Step D, at each sampling instant t, the running state current based on train and historical position observation sequence, predict the advanced positions in train certain moment following; Its detailed process is as follows:
Step D1, train track data pretreatment, with train at the stop position of originating station for the origin of coordinates, in each sampling instant, according to the train original discrete two-dimensional position sequence x=[x obtained
1, x
2..., x
n] and y=[y
1, y
2..., y
n], adopt first order difference method to carry out processing new train discrete location sequence Δ x=[the Δ x of acquisition to it
1, Δ x
2..., Δ x
n-1] and Δ y=[Δ y
1, Δ y
2..., Δ y
n-1], wherein Δ x
i=x
i+1-x
i, Δ y
i=y
i+1-y
i(i=1,2 ..., n-1);
Step D2, to train track data cluster, to new train discrete two-dimensional position sequence Δ x and Δ y after process, by setting cluster number M ', K-means clustering algorithm is adopted to carry out cluster to it respectively;
Step D3, HMM is utilized to carry out parameter training to the train track data after cluster, by the train operation track data Δ x after process and Δ y being considered as the aobvious observed value of hidden Markov models, upgrade period τ ' by setting hidden state number N ' and parameter, according to nearest T ' individual position detection value and adopt B-W algorithm to roll to obtain up-to-date HMM parameter lambda '; Specifically: because obtained train track sets data length is dynamic change, in order to the state variation of real-time tracking train track, be necessary at initial track HMM parameter lambda '=(π, A, B) basis is readjusted it, to infer the position of train in certain moment following more accurately; Every period τ ', the individual observed value (o of the T ' according to up-to-date acquisition
1, o
2..., o
t ') to track HMM parameter lambda '=(π, A, B) reappraise;
Step D4, foundation HMM parameter, adopt the hidden state q corresponding to Viterbi algorithm acquisition current time observed value;
Step D5, every the period
according to the HMM parameter lambda of up-to-date acquisition '=(π, A, B) and nearest H conception of history measured value (o
1, o
2..., o
h), based on the hidden state q of train current time, at moment t, by setting prediction time domain h ', obtain the position prediction value O of future time period train;
The value of above-mentioned cluster number M ' is 4, and the value of hidden state number N ' is 3, and it is 30 seconds that parameter upgrades period τ ', and T ' is 10,
be 30 seconds, H is 10, and prediction time domain h ' is 300 seconds.
Step e, seeing Fig. 3, set up the dynamic continuously observer to discrete conflict logic from train, is the conflict situation that discrete observation value is expressed by the continuous dynamic mapping of subway transportation system; When system likely violates traffic control rule, to the Hybrid dynamics behavior implementing monitoring of subway transportation hybrid system, for subway transportation control center provides warning information timely;
The specific implementation process of described step e is as follows:
Step F 1, construct conflict hypersurface collection of functions based on regulation rule: set up hypersurface collection of functions in order to reflect the contention situation of system, wherein, continuous function h relevant to single train in conflict hypersurface
ibe I type hypersurface, the continuous function h relevant to two trains
iIit is II type hypersurface;
Step e 2, set up by train continuous state to the observer of discrete conflict situation, build the safety rule collection d that need meet when train runs in traffic network
ij(t)>=d
min, wherein d
ijt () represents the actual interval of train i and train j in t, d
minrepresent the minimum safety interval between train;
Step e 3, based on person machine system theoretical and complication system hierarchical control principle, according to train operation pattern, build people at the real-time monitoring mechanism of the train of loop, the operation of guarantee system is in safe reachable set, design the discrete monitor from conflict to conflict Resolution means, when the discrete observation vector of observer shows that safety rule rally is breached, send corresponding warning information to subway transportation control center at once.
Obviously, above-described embodiment is only for example of the present invention is clearly described, and is not the restriction to embodiments of the present invention.For those of ordinary skill in the field, can also make other changes in different forms on the basis of the above description.Here exhaustive without the need to also giving all embodiments.And these belong to spirit institute's apparent change of extending out of the present invention or change and are still among protection scope of the present invention.
Claims (6)
1. a subway train conflict method for early warning, is characterized in that comprising the steps:
Steps A, plan operational factor according to each train, the topology diagram of generator orbital traffic network;
Step B, topology diagram based on the Rail traffic network constructed by steps A, analyze controllability and sensivity two class feature of train flow;
Step C, plan operational factor according to each train, on the basis building Modeling Method for Train Dynamics, set up train operation conflict according to train operation conflict Coupling point and allocate model in advance, generate multiple row car Lothrus apterus running orbit;
Step D, at each sampling instant t, the running state current based on train and historical position observation sequence, predict the advanced positions in train certain moment following; Its detailed process is as follows:
Step D1, train track data pretreatment, with train at the stop position of originating station for the origin of coordinates, in each sampling instant, according to the train original discrete two-dimensional position sequence x=[x obtained
1, x
2..., x
n] and y=[y
1, y
2..., y
n], adopt first order difference method to carry out processing new train discrete location sequence Δ x=[the Δ x of acquisition to it
1, Δ x
2..., Δ x
n-1] and Δ y=[Δ y
1, Δ y
2..., Δ y
n-1], wherein Δ x
i=x
i+1-x
i, Δ y
i=y
i+1-y
i(i=1,2 ..., n-1);
Step D2, to train track data cluster, to new train discrete two-dimensional position sequence Δ x and Δ y after process, by setting cluster number M ', K-means clustering algorithm is adopted to carry out cluster to it respectively;
Step D3, HMM is utilized to carry out parameter training to the train track data after cluster, by the train operation track data Δ x after process and Δ y being considered as the aobvious observed value of hidden Markov models, upgrade period τ ' by setting hidden state number N ' and parameter, according to nearest T ' individual position detection value and adopt B-W algorithm to roll to obtain up-to-date HMM parameter lambda '; Specifically: because obtained train track sets data length is dynamic change, in order to the state variation of real-time tracking train track, be necessary at initial track HMM parameter lambda '=(π, A, B) basis is readjusted it, to infer the position of train in certain moment following more accurately; Every period τ ', the individual observed value (o of the T ' according to up-to-date acquisition
1, o
2..., o
t ') to track HMM parameter lambda '=(π, A, B) reappraise;
Step D4, foundation HMM parameter, adopt the hidden state q corresponding to Viterbi algorithm acquisition current time observed value;
Step D5, every the period
according to the HMM parameter lambda of up-to-date acquisition '=(π, A, B) and nearest H conception of history measured value (o
1, o
2..., o
h), based on the hidden state q of train current time, at moment t, by setting prediction time domain h ', obtain the position prediction value O of future time period train;
Step e, setting up the dynamic continuously observer to discrete conflict logic from train, is the conflict situation that discrete observation value is expressed by the continuous dynamic mapping of subway transportation system; When system likely violates traffic control rule, to the Hybrid dynamics behavior implementing monitoring of subway transportation hybrid system, for subway transportation control center provides warning information timely.
2. a kind of subway train conflict method for early warning according to claim 1, is characterized in that: the detailed process of steps A is as follows:
Steps A 1, extract the site information of stopping each train travelling process from the data bank of subway transportation control center;
Steps A 2, according to positive and negative two service directions, the site information that each train is stopped to be classified, and the same site on same service direction is merged;
Steps A 3, according to website amalgamation result, according to website space layout's form with straight line connect before and after multiple website.
3. a kind of subway train conflict method for early warning according to claim 1, is characterized in that: the detailed process of step B is as follows:
Step R1, the Traffic flux detection model built in single subsegment; Its detailed process is as follows:
Step R1.1, introducing state variable Ψ, input variable u and output variable Ω, wherein Ψ represents the train quantity of certain moment existence on connected section between website, it comprises single channel section and Multiple Sections two type, u represents the Operation Measures that track traffic dispatcher implements for certain section, as adjust train speed or change train in the station time etc., Ω represents the train quantity that certain period section is left;
Step R1.2, by by time discretization, set up shape as Ψ (t+ Δ t)=A
1Ψ (t)+B
1u (t) and Ω (t)=C
1Ψ (t)+D
1discrete time Traffic flux detection model in the single subsegment of u (t), wherein Δ t represents the sampling interval, and Ψ (t) represents the state vector of t, A
1, B
1, C
1and D
1represent the state-transition matrix of t, input matrix, output calculation matrix and direct transmission matrix respectively;
Step B2, the Traffic flux detection model built in many subsegments; Its detailed process is as follows:
Step B2.1, according to circuit space layout form and train flow historical statistical data, obtain the flow proportional parameter beta in cross link each subsegment;
Step B2.2, according to the discrete time Traffic flux detection model in flow proportional parameter and single subsegment, build shape as Ψ (t+ Δ t)=A
1Ψ (t)+B
1u (t) and Ω (t)=C
1Ψ (t)+D
1discrete time Traffic flux detection model in many subsegments of u (t);
Step B3, controllable factor matrix [B according to Controlling model
1, A
1b
1..., A
1 n-1b
1] order and the relation of numerical value n, its controllability of qualitative analysis, according to the sensitivity coefficient matrix [C of Controlling model
1(zI-A
1)
-1b
1+ D
1], its input and output sensivity of quantitative analysis, wherein n represents the dimension of state vector, I representation unit matrix, and z represents the element factor changed original discrete time Traffic flux detection model.
4. a kind of subway train conflict method for early warning according to claim 1, is characterized in that: the detailed process of step C is as follows:
Step C1, train status transfer modeling, train shows as the switching at runtime process between website along the process that track traffic road network runs, arrange according to the website in train operation plan, set up the Petri network model of single train switched and transferred between different website: E=(g, G, Pre, Post, m) be train section metastasis model, wherein g represents each sub-section between website, G represents the change-over point of train running speed state parameter, Pre and Post represents that front and back between each sub-section and website are to annexation respectively
represent the operation section residing for train, wherein m represents model identification, and Z+ represents Positive Integer Set;
The modeling of step C2, train full operation profile hybrid system, the operation of train between website is considered as continuous process, from the stressed situation of train, according to the kinetics equation of energy model derivation train in the different operation phase, in conjunction with external interference factor, set up about train in a certain operation phase speed v
gmapping function v
g=λ (T
1, T
2, H, R, α), wherein T
1, T
2, H, R and α represent tractive force of train, braking force of train, train resistance, train gravity and train status random fluctuation parameter respectively;
Step C3, adopt the mode mixing emulation to infer to solve train track, by by time subdivision, the distance of utilization state continually varying characteristic Recursive Solution any time train in a certain operation phase apart from initial rest position point,
wherein J
0for initial time train is apart from the voyage of initial rest position point, Δ τ is the numerical value of time window, and J (τ), for τ moment train is apart from the distance of initial rest position point, can infers thus and obtain single vehicles track;
Step C4, train, in station time probability distribution function modeling, for specific run circuit, by transferring the dwell time data of train at each station, obtain the dwell time probability distribution of different circuit different website condition Train;
The Lothrus apterus robust track allotment of step C5, multiple row car coupling, the time of conflict point is reached in advance according to each train, pass through Time segments division, at each sampling instant t, under the prerequisite incorporating random factor, according to scheduling rule, robust secondary planning is implemented to the train track not meeting personal distance requirement near conflict point.
5. a kind of subway train conflict method for early warning according to claim 1, it is characterized in that: in step D, the value of cluster number M ' is 4, and the value of hidden state number N ' is 3, and it is 30 seconds that parameter upgrades period τ ', and T ' is 10,
be 30 seconds, H is 10, and prediction time domain h ' is 300 seconds.
6. a kind of subway train conflict method for early warning according to claim 1, is characterized in that: the specific implementation process of step e is as follows:
Step F 1, construct conflict hypersurface collection of functions based on regulation rule: set up hypersurface collection of functions in order to reflect the contention situation of system, wherein, continuous function h relevant to single train in conflict hypersurface
ibe I type hypersurface, the continuous function h relevant to two trains
iIit is II type hypersurface;
Step e 2, set up by train continuous state to the observer of discrete conflict situation, build the safety rule collection d that need meet when train runs in traffic network
ij(t)>=d
min, wherein d
ijt () represents the actual interval of train i and train j in t, d
minrepresent the minimum safety interval between train;
Step e 3, based on person machine system theoretical and complication system hierarchical control principle, according to train operation pattern, build people at the real-time monitoring mechanism of the train of loop, the operation of guarantee system is in safe reachable set, design the discrete monitor from conflict to conflict Resolution means, when the discrete observation vector of observer shows that safety rule rally is breached, send corresponding warning information to subway transportation control center at once.
Priority Applications (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710078598.XA CN106864482A (en) | 2015-03-31 | 2015-03-31 | Subway conflict early warning method |
| CN201710078933.6A CN106926873A (en) | 2015-03-31 | 2015-03-31 | Subway train collision early warning method based on robust strategy |
| CN201510150113.4A CN105083322B (en) | 2015-03-31 | 2015-03-31 | Subway train collision early warning method |
| CN201710078599.4A CN106864483A (en) | 2015-03-31 | 2015-03-31 | Subway train collision early warning method based on robust strategy |
| CN201710078606.0A CN106741020A (en) | 2015-03-31 | 2015-03-31 | Subway conflict early warning method |
| CN201710078411.6A CN106741006A (en) | 2015-03-31 | 2015-03-31 | Subway conflict early warning method based on robust strategy |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510150113.4A CN105083322B (en) | 2015-03-31 | 2015-03-31 | Subway train collision early warning method |
Related Child Applications (5)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710078599.4A Division CN106864483A (en) | 2015-03-31 | 2015-03-31 | Subway train collision early warning method based on robust strategy |
| CN201710078411.6A Division CN106741006A (en) | 2015-03-31 | 2015-03-31 | Subway conflict early warning method based on robust strategy |
| CN201710078606.0A Division CN106741020A (en) | 2015-03-31 | 2015-03-31 | Subway conflict early warning method |
| CN201710078598.XA Division CN106864482A (en) | 2015-03-31 | 2015-03-31 | Subway conflict early warning method |
| CN201710078933.6A Division CN106926873A (en) | 2015-03-31 | 2015-03-31 | Subway train collision early warning method based on robust strategy |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105083322A true CN105083322A (en) | 2015-11-25 |
| CN105083322B CN105083322B (en) | 2017-03-08 |
Family
ID=54564882
Family Applications (6)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710078599.4A Pending CN106864483A (en) | 2015-03-31 | 2015-03-31 | Subway train collision early warning method based on robust strategy |
| CN201710078606.0A Pending CN106741020A (en) | 2015-03-31 | 2015-03-31 | Subway conflict early warning method |
| CN201710078411.6A Pending CN106741006A (en) | 2015-03-31 | 2015-03-31 | Subway conflict early warning method based on robust strategy |
| CN201710078598.XA Pending CN106864482A (en) | 2015-03-31 | 2015-03-31 | Subway conflict early warning method |
| CN201710078933.6A Pending CN106926873A (en) | 2015-03-31 | 2015-03-31 | Subway train collision early warning method based on robust strategy |
| CN201510150113.4A Active CN105083322B (en) | 2015-03-31 | 2015-03-31 | Subway train collision early warning method |
Family Applications Before (5)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710078599.4A Pending CN106864483A (en) | 2015-03-31 | 2015-03-31 | Subway train collision early warning method based on robust strategy |
| CN201710078606.0A Pending CN106741020A (en) | 2015-03-31 | 2015-03-31 | Subway conflict early warning method |
| CN201710078411.6A Pending CN106741006A (en) | 2015-03-31 | 2015-03-31 | Subway conflict early warning method based on robust strategy |
| CN201710078598.XA Pending CN106864482A (en) | 2015-03-31 | 2015-03-31 | Subway conflict early warning method |
| CN201710078933.6A Pending CN106926873A (en) | 2015-03-31 | 2015-03-31 | Subway train collision early warning method based on robust strategy |
Country Status (1)
| Country | Link |
|---|---|
| CN (6) | CN106864483A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110162912A (en) * | 2019-05-30 | 2019-08-23 | 北京航空航天大学 | Consider the network resource optimization method of resource contention coupling |
| CN111152818A (en) * | 2020-01-03 | 2020-05-15 | 中国铁道科学研究院集团有限公司通信信号研究所 | Operation diagram stage plan logic checking method |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108116456A (en) * | 2017-12-18 | 2018-06-05 | 湖南中车时代通信信号有限公司 | For the storage system of remote backup train operation monitoring device data |
| CN109552364A (en) * | 2018-11-19 | 2019-04-02 | 东莞理工学院 | A kind of rail traffic scheduling early warning system based on big data |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN203211332U (en) * | 2013-04-18 | 2013-09-25 | 南京理工大学 | Anti-collision pre-warning device for subway trains |
| DE102013210081A1 (en) * | 2013-05-29 | 2014-12-04 | Siemens Aktiengesellschaft | Method for detecting a hazard to a rail vehicle |
| CN104462856A (en) * | 2014-12-30 | 2015-03-25 | 江苏理工学院 | Ship conflict early warning method |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8812227B2 (en) * | 2011-05-19 | 2014-08-19 | Metrom Rail, Llc | Collision avoidance system for rail line vehicles |
| CN102336204A (en) * | 2011-07-31 | 2012-02-01 | 宁波市镇海西门专利技术开发有限公司 | Anti-collision system for trains |
| CN102700570A (en) * | 2012-05-22 | 2012-10-03 | 西南交通大学 | Anti-collision and early-warning system of railway vehicles |
| CN103310118B (en) * | 2013-07-04 | 2016-02-10 | 中国铁路总公司 | High speed rail train operation conflict prediction method |
-
2015
- 2015-03-31 CN CN201710078599.4A patent/CN106864483A/en active Pending
- 2015-03-31 CN CN201710078606.0A patent/CN106741020A/en active Pending
- 2015-03-31 CN CN201710078411.6A patent/CN106741006A/en active Pending
- 2015-03-31 CN CN201710078598.XA patent/CN106864482A/en active Pending
- 2015-03-31 CN CN201710078933.6A patent/CN106926873A/en active Pending
- 2015-03-31 CN CN201510150113.4A patent/CN105083322B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN203211332U (en) * | 2013-04-18 | 2013-09-25 | 南京理工大学 | Anti-collision pre-warning device for subway trains |
| DE102013210081A1 (en) * | 2013-05-29 | 2014-12-04 | Siemens Aktiengesellschaft | Method for detecting a hazard to a rail vehicle |
| CN104462856A (en) * | 2014-12-30 | 2015-03-25 | 江苏理工学院 | Ship conflict early warning method |
Non-Patent Citations (1)
| Title |
|---|
| 严波: ""城市轨道交通行车间隔时间优化模型研究"", 《城市轨道交通研究》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110162912A (en) * | 2019-05-30 | 2019-08-23 | 北京航空航天大学 | Consider the network resource optimization method of resource contention coupling |
| CN110162912B (en) * | 2019-05-30 | 2020-09-01 | 北京航空航天大学 | Network resource optimization method considering resource competition coupling |
| CN111152818A (en) * | 2020-01-03 | 2020-05-15 | 中国铁道科学研究院集团有限公司通信信号研究所 | Operation diagram stage plan logic checking method |
| CN111152818B (en) * | 2020-01-03 | 2021-12-28 | 中国铁道科学研究院集团有限公司通信信号研究所 | Operation diagram stage plan logic checking method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106926873A (en) | 2017-07-07 |
| CN106741006A (en) | 2017-05-31 |
| CN106864482A (en) | 2017-06-20 |
| CN106741020A (en) | 2017-05-31 |
| CN106864483A (en) | 2017-06-20 |
| CN105083322B (en) | 2017-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105083333B (en) | Subway traffic flow optimization control method | |
| Dakic et al. | On the design of an optimal flexible bus dispatching system with modular bus units: Using the three-dimensional macroscopic fundamental diagram | |
| Yu et al. | Real-time partway deadheading strategy based on transit service reliability assessment | |
| CN105095984A (en) | Real-time prediction method for subway train track | |
| CN105083334A (en) | Subway train collision early warning method | |
| CN105083322A (en) | Subway train collision early warning method | |
| CN105083335B (en) | Subway traffic flow optimization control method | |
| CN103106801A (en) | Self-organizing traffic signal coordination control method | |
| CN105095983A (en) | Real-time prediction method for subway train track | |
| Xiong et al. | Parallel bus rapid transit (BRT) operation management system based on ACP approach | |
| Oskarbski et al. | Multi-level transport systems model for traffic management activities | |
| CN105093929A (en) | Planning method for conflict resolution of subway train | |
| Xin et al. | Modeling bus bunching and anti-bunching control accounting for signal control and passenger swapping behavior | |
| Melki et al. | Advanced approach for the public transportation regulation system based on cybercars | |
| Hajbabaie et al. | Dynamic Speed Harmonization in Connected Urban Street Networks: Improving Mobility | |
| Li et al. | Elevating adaptive traffic signal control in semi‐autonomous traffic dynamics by using connected and automated vehicles as probes | |
| Adamski et al. | ITS: Intelligent Transportation Systems-multi-criteria control problems |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CP02 | Change in the address of a patent holder |
Address after: No. 1801 Zhong Wu Avenue, Changzhou, Jiangsu Province, Jiangsu Patentee after: Jiangsu University of Technology Address before: 213001 1801 Zhong Wu Avenue, Zhong Lou District, Changzhou, Jiangsu Patentee before: Jiangsu University of Technology |
|
| CP02 | Change in the address of a patent holder |