CN107103142A - Comprehensive traffic network operation situation towards highway and the railway network deduces emulation technology - Google Patents

Comprehensive traffic network operation situation towards highway and the railway network deduces emulation technology Download PDF

Info

Publication number
CN107103142A
CN107103142A CN201710300803.2A CN201710300803A CN107103142A CN 107103142 A CN107103142 A CN 107103142A CN 201710300803 A CN201710300803 A CN 201710300803A CN 107103142 A CN107103142 A CN 107103142A
Authority
CN
China
Prior art keywords
network
traffic
node
prediction
road
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
Application number
CN201710300803.2A
Other languages
Chinese (zh)
Other versions
CN107103142B (en
Inventor
张凡
张纪升
牛树云
周学松
李宏海
孙晓亮
刘见平
张利
朱丽丽
赵丽
崔玮
张金金
王�华
王体彬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Research Institute of Highway Ministry of Transport
Original Assignee
Research Institute of Highway Ministry of Transport
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Research Institute of Highway Ministry of Transport filed Critical Research Institute of Highway Ministry of Transport
Priority to CN201710300803.2A priority Critical patent/CN107103142B/en
Publication of CN107103142A publication Critical patent/CN107103142A/en
Application granted granted Critical
Publication of CN107103142B publication Critical patent/CN107103142B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/40Business processes related to the transportation industry

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Geometry (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Development Economics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Traffic Control Systems (AREA)

Abstract

The present invention discloses a kind of comprehensive traffic network operation situation towards highway and the railway network and deduces emulation technology, it is characterised in that including:1)Regional complex transportation network is met an urgent need simulation framework, 2)The needs estimate and forecast model combined based on trend prediction with neutral net, 3)Multi-mode comprehensive transport capability simulation model based on intelligent body;This is deduced emulation technology towards the comprehensive traffic network operation situation of highway and the railway network and utilizes existing network of highways, railway network traffic data, emphasis constructs the general frame of the emergent emulation platform of regional complex transportation network, with time-space network modeling technique, the multi-mode comprehensive traffic simulation model based on intelligent body is constructed.

Description

Comprehensive traffic network operation situation towards highway and the railway network deduces emulation technology
Technical field
The invention belongs to intelligent transportation field, specifically a kind of comprehensive traffic network towards highway and the railway network is transported Row deducing manoeuver emulation technology.
Background technology
In regional response process of emergency system, how to formulate the emergent Traffic Organization of comprehensive traffic road network ensures emergent friendship Lead to safe efficient, orderly function, for improving region emergency managerial ability, casualties and property caused by reduction disaster Loss, has important practical significance.For transportation network, particularly composite communications transport network, it is difficult to pass through mathematics Its running status changing trend under event condition of mode quick obtaining such as calculate and to different management and control measures or emergent ring Answer the visual effect of measure.Therefore, set up regional complex transportation network operation situation and deduce emulation, realize correspondence network of highways The quick loading on composite communications transport network of anxious management and control measures is emulated with deduction, is to evaluate and the emergent traffic of optimization road network The important means of organization scheme.
Single travel pattern being intended for existing traffic simulation platform, less research is more in composite communications transport network more Plant under conditions of travel pattern exists jointly and how to build emulation platform.And existing analogue system, mostly single-threaded system, Simulation velocity is partially slow, it is impossible to which adaptation is docked with the real-time of operation system.
The content of the invention
It is an object of the invention to provide one kind using existing network of highways, railway network traffic data, emphasis constructs region The general frame of the emergent emulation platform of composite communications transport network, with time-space network modeling technique, is constructed based on intelligent body Multi-mode comprehensive traffic simulation model towards highway and the railway network comprehensive traffic network operation situation deduce emulation technology.
In order to solve the above-mentioned technical problem, the technical scheme is that:
A kind of comprehensive traffic network operation situation towards highway and the railway network deduces emulation technology, including:
1) the emergent simulation framework of regional complex transportation network,
Including comprehensive travel needs estimate and prediction module, road-net node/section dynamic function Attribute tuning module, multimode Formula dynamic traffic flow, traffic flow simulation, dynamic shortest path are calculated, operation is assessed, exported and interface display module, emulation system The key data input of system includes OD demand datas, road net traffic state service data real-time with history, incident data, Network foundation data and road management and control and information distribution scheme;
2) needs estimate and forecast model combined based on trend prediction with neutral net,
Comprehensive travel needs estimate is with prediction module mainly by real-time and history OD demand datas, road grid traffic shape State service data etc. carries out pattern analysis and OD trend estimations, and Mobile state tune is entered in influence of the combination accident to transport need It is whole, generate road network Dynamic OD requirement matrix;
3) the multi-mode comprehensive transport capability simulation model based on intelligent body,
Multi-mode dynamic traffic mainly solves the Traffic Assignment Problem between a variety of travel patterns with flow module, by institute Generation, the calculating for the selecting module for going out row mode for having path set, by vehicles distribution corresponding with people, realize intelligent body Travel schedule, to each path of the intelligent body stroke one comprising timestamp sequence and sequence node.
Further, road-net node/section dynamic function attribute module is the basic data by network, with reference to dynamic Incident data and road management and control data and information distribution scheme, generation emulation road-net node/section dynamic function attribute.
Further, multi-mode dynamic traffic passes through to all generations for having path set, the choosing for going out row mode with flow module The calculating of module is selected, by vehicles distribution corresponding with people, the travel schedule of intelligent body is realized, to each intelligent body stroke one Path comprising timestamp sequence and sequence node.
Further, traffic flow simulation and dynamic shortest path module constitute the nucleus module of middle sight traffic simulation, hand over Through-flow analog module is using section transmission, node-node transmission model, on the premise of relevant constraint is met, and realization exists to vehicle The simulation shifted between section internal operation and section, section mode is related to the magnitude of traffic flow-density relations, and theory is current Ability is main by road segment classification and track quantity;Node metastasis model is related to left-hand rotation and straight trip Capacity Constraints, dynamic point Hourage of the traffic flow simulated with platform on section will be fed back in the dynamic shortest path model of next iteration, and then Decision-making foundation is provided for intelligent body selection specially designated path, meanwhile, new pathway can be fed back to next iteration.
Further, operation evaluation module is estimated to operation result from section, passage, OD to, path, region, is wrapped In terms of including efficiency, speed, reliability, energy consumption.
Further, predicted using polynomial trend and combine vehicle flowrate or trip of the kalman filter method to Ingress node The volume of the flow of passengers estimated and predicted respectively, and method is as follows:
J=road-net node sequence numbers, j=1,2 ..., N;
τ=departure time interval sequence number, τ=1,2 ...;
The stage of k=rolling forecasts, k=1,2,3 ...;
L=each rolls the departure time intervening sequence number included in the stage;
D(j,τ)The demand influxs of=node j in time interval τ;
Prior estimates of=node the j based on historical data in time interval τ demand influx;
μ(j,τ)=node j is in time interval τ D(j,τ)With prior estimateDemand disruption;
μ′(j,τ),μ″(j,τ),μ″′(j,τ)=demand disruption μ(j,τ)The single order of variable, second order, third order difference amount;
=D(j,τ)Estimate average;
(p) (j,τ)Estimation average;
Zk=stage k state variable vector;
Yk=stage k measurement vector;
Hk=measurement vector YkWith state vector ZkMatching relationship matrix;
wk=stage k processing noise;
Rk=stage k measurement error;
=the Z based on the observation prediction by the end of stage k-1kThat is E (Zk|Y1,Y2,…,Yk-1);
=the Z based on the observation estimation by the end of the k stageskThat is E (Zk|Y1,Y2,…,Yk);
Pk,k-1The Z of=stage k-1 predictionskState covariance battle array, i.e.,
Pk,k=stage k, the Z of estimationkState covariance battle array, i.e.,
Polynomial trend prediction is introduced, the prediction of τ+ζ difference vector is carried out, the function is built such as using m order polynomials Under:
μ(j,τ+ζ)=b0+b1ζ+b2ζ2+…bpζp+…+bmζm
According to Taylor expansion anticipation function μ(j,τ+ζ)It is extended to μ(j,τ)Multistage difference combination:
The coefficient of the formula of original polynomial distribution can be expressed as:
Three rank multinomial models are taken to carry out μ(j,τ+ζ)Prediction be shown below:
We are defined on forecast period k, τ=kl, and such as each stage forecast Period Length is that l includes 3 time intervals, that In the k=3 stages, τ sequence number is 9;
State-transition matrix is defined as:
The predictor formula in k+1 stages is as follows:
Zk+1=AkZk
Introduce Kalman filtering algorithm to eliminate the influence that charge station's influx measured value error band comes, because this part is straight It is identical vector value that measured value, which is connect, with estimate, therefore H=I;
Carried out first with the posterior estimate of previous stage, the prediction of current generation and the transmission of covariance:
Introduce afterwards after stage k measured value, carry out the calculating of gain:
The transmission of formation stages k posterior estimate and covariance matrix:
Pk,k=(I-KkHk)Pk,k-1
Call above-mentioned formula forecast period k+1 μ(j,τ+ζ), the node-flow for deduction can be generated based on the predicted value Enter value D(j,τ), and the multistage difference matrix formed needed for subsequent prediction.
Further, by BP neural network model can respectively in freeway net vehicle OD distribution be predicted and The OD for going out administrative staff in the railway network is distributed and is predicted, illustrates to calculate by taking the OD forecast of distribution of vehicle in freeway net as an example Method flow, algorithm flow is as follows:
Step 1:Data decimation and standardization, for n charge tiny node in road network, choose before predicted time section The OD data of m days same periods are simultaneously standardized according to formula 4.1;
Step 2:BP neural network is built, i=1 is made, for i-th of charge tiny node, the OD distributions for choosing first day are compared Example chooses the data of second day as desired output vector as input vector, the like constitute m-1 group training samples;
Step 3:Neutral net is trained, by repeatedly training, network is reached the training objective of setting;
Step 4:OD distribution proportions are predicted, for the neutral net trained, the number of input prediction objective time interval the previous day According to obtaining prediction distribution ratio, if going to step 5, otherwise go to step 2, adjust network parameter, re -training network;
Step 5:Output prediction OD distribution proportion matrixes, the output matrix if i=n, otherwise i+1 go to step 2.
Further, predict obtain road grid traffic respectively with BP neural network algorithm and polynomial trend prediction algorithm After distribution proportion and each entrance magnitude of traffic flow, two kinds are predicted the outcome and is merged, obtained under the conditions of no event, predicted on network The passenger of period or the OD matrixes of vehicle, by vehicle or passenger by node i to node j flow qijWebsite i entrances can be passed through Flow OiWith with the distribution proportion r for going to j websitesijIt is determined that, under the conditions of without event, road network predicts OD matrixes by equation below Try to achieve:
By the OD matrixes OD of above-mentioned formula respectively to vehiclevWith the OD matrixes OD of passengerpCalculated, so as to obtain pre- Survey the OD matrixes without the vehicle under the conditions of event and passenger on the period.
Further, Travel Time Budget TTB (p), starting point o are givenp, intelligent body p departure time τp, considering road Under the constraints of appearance of a street amount and Passenger Flow service ability, maximize the space-time accessibility in system scope, based on intelligence The object function of the optimization problem of the composite communications transport network dynamic distribution model of body, is shown below:
That is, the optimization aim of model is to minimize the unreachable generalized cost of user.If intelligent body p can be pre- in the given time Arrived in calculating, then inaccessible generalized cost is 0, otherwise, C (p) is by overall travel time of the intelligent body using path.
The technology of the present invention effect major embodiment is in the following areas:In regional response process of emergency system, how synthesis is formulated The emergent Traffic Organization of traffic network ensures emergent traffic safety, efficient, orderly function, for improving region contingency management energy Power, reduces the casualties and property loss caused by disaster, has important practical significance, for transportation network, special It is not composite communications transport network, becomes it is difficult to mode quick obtaining its running status under event condition such as calculate by mathematics Change situation and to different management and control measures or the visual effect of emergency response measure, therefore, set up regional complex transportation network Network emulation platform, realizes that the quick loading on composite communications transport network for management and control measures of being met an urgent need to network of highways is emulated, is to comment Valency and the important means of the emergent Traffic Organization of optimization road network.
Brief description of the drawings
Fig. 1 is the emergent simulation framework figure of regional complex transportation network of the invention;
Fig. 2 is comprehensive travel needs estimate of the invention and prediction Organization Chart;
Fig. 3 is the trip distribution modeling algorithm flow chart based on BP neural network of the invention;
Fig. 4 is the algorithm flow chart under the finite path collection based on column-generation of the present invention;
Fig. 5 is input and output design drawing of the invention.
Embodiment
A kind of comprehensive traffic network operation situation towards highway and the railway network deduces emulation technology, including
1) the emergent simulation framework of regional complex transportation network
Regional complex transportation network meets an urgent need emulation platform mainly including comprehensive travel needs estimate and prediction module, road Net node/section dynamic function Attribute tuning module, multi-mode dynamic traffic flow, traffic flow simulation, dynamic shortest path meter Calculate, run several modules such as assessment, and it is pre- with the transportation network security risk based on G I S by way of interface interchange Alert platform and Decision System of Emergency connect relatively, architecture design such as Fig. 1 institutes of the emergent emulation platform of regional complex transportation network Show.
The key data input of analogue system includes OD demand datas in real time with history, road net traffic state operation number According to incident data, network foundation data and road management and control data and information distribution scheme.
Comprehensive travel needs estimate is with prediction module mainly by real-time and history OD demand datas, road grid traffic shape State service data etc. carries out pattern analysis and OD trend estimations, and Mobile state tune is entered in influence of the combination accident to transport need It is whole, generate road network Dynamic OD requirement matrix.
Road-net node/section dynamic function attribute module is mainly the basic data by network, with reference to dynamic burst Event data and road management and control data and information distribution scheme, generation emulation road-net node/section dynamic function attribute.
Multi-mode dynamic traffic mainly solves the Traffic Assignment Problem between a variety of travel patterns with flow module, by institute Generation, the calculating for the selecting module for going out row mode for having path set, by vehicles distribution corresponding with people, realize intelligent body Travel schedule, to each path of the intelligent body stroke one comprising timestamp sequence and sequence node.
Traffic flow simulation and dynamic shortest path module constitute the nucleus module of middle sight traffic simulation, traffic flow simulation mould Block is using section transmission, node-node transmission model, on the premise of relevant constraint is met, and realizes to vehicle inside section The simulation shifted between operation and section.Section mode is related to the magnitude of traffic flow-density relations, and basic capacity is main By road segment classification and track quantity;Node metastasis model is related to left-hand rotation and straight trip Capacity Constraints.Dynamically distributes platform institute Hourage of the traffic flow of simulation on section will be fed back in the dynamic shortest path model of next iteration, and then is intelligent body Specially designated path is selected to provide decision-making foundation.Meanwhile, new pathway can be fed back to next iteration.
Operation evaluation module can be commented multiple dimensions such as, path, regions operation result from section, passage, OD Estimate, including many aspects such as efficiency, speed, reliability, energy consumption.
The input and output of the emergent emulation platform of regional complex transportation network enter shown in Fig. 5.
2) needs estimate and forecast model combined based on trend prediction with neutral net;
Needs estimate and prediction Organization Chart such as Fig. 2.
1. entrance needs estimate
Predicted using polynomial trend and combine vehicle flowrate or passenger traffic volume point of the kalman filter method to Ingress node Do not estimated and predicted.
J=road-net node sequence numbers, j=1,2 ..., N;
τ=departure time interval sequence number, τ=1,2 ...;
The stage of k=rolling forecasts, k=1,2,3 ...;
L=each rolls the departure time intervening sequence number included in the stage;
D(j,τ)The demand influxs of=node j in time interval τ;
Prior estimates of=node the j based on historical data in time interval τ demand influx;
μ(j,τ)=node j is in time interval τ D(j,τ)With prior estimateDemand disruption;
μ′(j,τ),μ″(j,τ),μ″′(j,τ)=demand disruption μ(j,τ)The single order of variable, second order, third order difference amount;
=D(j,τ)Estimate average;
(p) (j,τ)Estimation average;
Zk=stage k state variable vector;
Yk=stage k measurement vector;
Hk=measurement vector YkWith state vector ZkMatching relationship matrix;
Wk=stages k processing noise;
Rk=stage k measurement error;
=the Z based on the observation prediction by the end of stage k-1kThat is E (Zk|Y1,Y2,…,Yk-1);
=the Z based on the observation estimation by the end of the k stageskThat is E (Zk|Y1,Y2,…,Yk);
Pk,k-1The Z of=stage k-1 predictionskState covariance battle array, i.e.,
Pk,k=stage k, the Z of estimationkState covariance battle array, i.e.,
Polynomial trend prediction is introduced, the prediction of τ+ζ difference vector is carried out, the function is built such as using m order polynomials Under:
μ(j,τ+ζ)=b0+b1ζ+b2ζ2+…bpζp+…+bmζm
According to Taylor expansion anticipation function μ(j,τ+ζ)It is extended to μ(j,τ)Multistage difference combination:
The coefficient of the formula of original polynomial distribution can be expressed as:
Three rank multinomial models are taken to carry out μ(j,τ+ζ)Prediction be shown below:
We are defined on forecast period k, τ=kl, and such as each stage forecast Period Length is that l includes 3 time intervals, that In the k=3 stages, τ sequence number is 9;
State-transition matrix is defined as:
The predictor formula in k+1 stages is as follows:
Zk+1=AkZk
Introduce Kalman filtering algorithm to eliminate the influence that charge station's influx measured value error band comes, because this part is straight It is identical vector value that measured value, which is connect, with estimate, therefore H=I;
Carried out first with the posterior estimate of previous stage, the prediction of current generation and the transmission of covariance:
Introduce afterwards after stage k measured value, carry out the calculating of gain:
The transmission of formation stages k posterior estimate and covariance matrix:
Pk,k=(I-KkHk)Pk,k-1
Call above-mentioned formula forecast period k+1 μ(j,τ+ζ), the node-flow for deduction can be generated based on the predicted value Enter value D(j,τ), and the multistage difference matrix formed needed for subsequent prediction.
2. network OD distribution proportions are predicted
OD distributions that can be respectively to vehicle in freeway net by BP neural network model are predicted and to the railway network In go out the OD distributions of administrative staff and be predicted, pre- flow gauge is similar, as shown in figure 3, therefore, only with vehicle in freeway net Illustrate algorithm flow exemplified by OD forecast of distribution, algorithm flow is as follows.
Step 1:Data decimation and standardization, for n charge tiny node in road network, choose before predicted time section The OD data of m days same periods are simultaneously standardized according to formula 4.1;
Step 2:BP neural network is built, i=1 is made, for i-th of charge tiny node, the OD distributions for choosing first day are compared Example chooses the data of second day as desired output vector as input vector, the like constitute m-1 group training samples;
Step 3:Neutral net is trained, by repeatedly training, network is reached the training objective of setting;
Step 4:OD distribution proportions are predicted, for the neutral net trained, the number of input prediction objective time interval the previous day According to obtaining prediction distribution ratio rij(j=1,2,3 ..., n), ifAnd rij>=0, step 5 is gone to, step is otherwise gone to Rapid 2, adjust network parameter, re -training network;
Step 5:Output prediction OD distribution proportion matrixes, the output matrix if i=n, otherwise i+1 go to step 2.
3. comprehensive travel demand OD Matrix predictions
Predict obtain road grid traffic distribution proportion respectively with BP neural network algorithm and polynomial trend prediction algorithm Merged with after each entrance magnitude of traffic flow, it is necessary to be predicted the outcome to two kinds, under the conditions of obtaining no event, prediction period on network Passenger or vehicle OD matrixes.By vehicle or passenger by node i to node j flow qijCan be by the stream of website i entrances Measure OiWith with the distribution proportion r for going to j websitesijIt is determined that, then under the conditions of without event, road network prediction OD matrixes can be by as follows Formula is tried to achieve:
OD matrixes OD that can respectively to vehicle by above-mentioned formulavWith the OD matrixes OD of passengerpCalculated, so as to obtain The OD matrixes of vehicle and passenger under the conditions of predicted time Duan Shangwu events.
3) the multi-mode comprehensive transport capability simulation model based on intelligent body
In composite communications transport network, traveler is arrived at by different trip modes, such as railway, self-driving and group Conjunction mode (such as highway, railway transfer).The network of this multi-mode, which can be modeled as one, can carry out the multilayer of Used in Dynamic Traffic Assignment Secondary network, wherein, constraints is the crucial ginseng in the service of different aspects or the constraint of road passage capability, modeling process Number, variable symbol and implication are as shown in the table.
Given Travel Time Budget TTB (p), starting point op, intelligent body p departure time τp, considering some road capacities Under the constraints of Passenger Flow service ability, maximize the space-time accessibility in system scope, i.e., inaccessible degree Reach minimum.The object function of the optimization problem of composite communications transport network dynamic distribution model based on intelligent body, it is as follows Shown in formula.
That is, the optimization aim of model is to minimize the unreachable generalized cost of user.If intelligent body p can be pre- in the given time Arrived in calculating, then inaccessible generalized cost is 0, otherwise, C (p) is by overall travel time of the intelligent body using path. That is, object function is intended to make the quantity of destination unreachable in all travelers minimum.
There are multiple constraintss, i.e. space-time traffic flow Constraints of Equilibrium, traveler to go out row constraint, traffic flow and space in model Capacity-constrained and the constraint of railway transport of passengers service ability etc..
1. space-time traffic flow equilibrium constraint
For all p ∈ P
Each intelligent body p only one of which starting point and a terminal are represented, on spatio-temporal state summit in addition to the start and the end points only In accordance with flow equilibrium constraint, i.e., for intelligent body p, each spatio-temporal state summit only once reaches and once left.
2. intelligent body trip constraints
All p ∈ P
Intelligent body p, which is arrived at, either to be gone on a journey by spatio-temporal activity arc or is gone on a journey by virtual spacetime arc of motion.
3. traffic flow and path space capacity constraints
Traffic flow constraints, i.e., overall influx r constraints are described by seeing traffic flow model in spatial queue For:
To arbitrary (i, j) ∈ Ep, t ∈ T wherein, Ai,j(t) represent on section (i, j) Accumulative amount of reach, Di,j(t) the accumulative amount of setting out is represented, its calculation formula is respectively:
For section (i, j),The free flow travel time is represented,Represent road section length,Represent that section is blocked Density under state.
4) Passenger Flow service ability is constrained
To owning (i, j, t, s) ∈ AT
Passenger Flow service constraints:Train can only be travelled according to the service chart of formulation on specific space-time arc, in this feelings Under condition, its intelligent body quantity serviced is no more than vehicle bearing capacity.
5) rail-road transfer (arc) stand service ability constrain
To owning (i, j, t, s) ∈ AT
The service ability constraint of rail-road transfer stop:The transfer of different grades of junction is limited in one's ability, by hinge parking facility, The facilities services such as facility, safety devices ability of waiting is determined, between this part is also included from road forking node to junction Section relevant capacity consistency.
The derivation algorithm of model is the algorithm based on column-generation, and the method based on column-generation is as shown in Figure 4.
The problem of outer circulation and circulation solution internally are under finite path collection is generated based on col-generating arithmetic required The minimum user cost path of time correlation.In each outer circulation iteration m, the cost minimization routing algorithm of time correlation Applied to the path for marginal user cost's cost minimization that time correlation is found for every one O-D pairs and the time interval each left, New route (if any) is added in the subset Pm in present feasible path in m iteration, and one based on two-dimensional projection Gradient descent direction method be used to solving to define on Pm under finite path collection the problem of, if do not find new path or Default convergence standard is reached, algorithm terminates and exports the path flow of the time correlation obtained in current iteration.
The gradient descent direction method based on two-dimensional projection obtained by iteration, is formed based on being followed in col-generating arithmetic framework Ring, in each interior circulation iteration n, with the path flow rn of DNLE (d) model evaluations renewal, corresponding system user totle drilling cost TE (rn) and user cost, if the difference in two continuous iteration (i.e. TE (rn)-TE (rn-1)) between desired value Value is less than predetermined threshold value or reaches default convergence (such as n=Nmax), and interior circulation is terminated, and algorithm returns to outer circulation.
The technology of the present invention effect major embodiment is in the following areas:In regional response process of emergency system, how synthesis is formulated The emergent Traffic Organization of traffic network ensures emergent traffic safety, efficient, orderly function, for improving region contingency management energy Power, reduces the casualties and property loss caused by disaster, has important practical significance, for transportation network, special It is not composite communications transport network, becomes it is difficult to mode quick obtaining its running status under event condition such as calculate by mathematics Change situation and to different management and control measures or the visual effect of emergency response measure, therefore, set up regional complex transportation network Network emulation platform, realizes that the quick loading on composite communications transport network for management and control measures of being met an urgent need to network of highways is emulated, is to comment Valency and the important means of the emergent Traffic Organization of optimization road network.
The foregoing is only a specific embodiment of the invention, but protection scope of the present invention is not limited thereto, any The change or replacement expected without creative work, should all be included within the scope of the present invention.

Claims (9)

1. a kind of comprehensive traffic network operation situation towards highway and the railway network deduces emulation technology, it is characterised in that including:
1) the emergent simulation framework of regional complex transportation network,
Moved including comprehensive travel needs estimate and prediction module, road-net node/section dynamic function Attribute tuning module, multi-mode State traffic flow distribution, traffic flow simulation, dynamic shortest path are calculated, operation is assessed, exported and interface display module, analogue system Key data input includes OD demand datas, road net traffic state service data in real time with history, incident data, network Basic data and road management and control and information distribution scheme;
2) needs estimate and forecast model combined based on trend prediction with neutral net,
Comprehensive travel needs estimate is with prediction module mainly by being transported to real-time and history OD demand datas, road net traffic state Row data etc. carry out pattern analysis and OD trend estimations, and Mobile state adjustment is entered in influence of the combination accident to transport need, Generate road network Dynamic OD requirement matrix;
3) the multi-mode comprehensive transport capability simulation model based on intelligent body,
Multi-mode dynamic traffic mainly solves the Traffic Assignment Problem between a variety of travel patterns with flow module, by having to all The generation of path set, go out row mode selecting module calculating, by the vehicles it is corresponding with people distribution, realize the stroke of intelligent body Scheme, to each path of the intelligent body stroke one comprising timestamp sequence and sequence node.
2. the comprehensive traffic network operation situation towards highway and the railway network deduces emulation technology according to claim 1, its It is characterised by:Road-net node/section dynamic function attribute module is the basic data by network, with reference to dynamic accident Data and road management and control data and information distribution scheme, generation emulation road-net node/section dynamic function attribute.
3. the comprehensive traffic network operation situation towards highway and the railway network deduces emulation technology according to claim 2, its It is characterised by:Multi-mode dynamic traffic passes through to all selecting modules for having the generation of path set, going out row mode with flow module Calculate, by vehicles distribution corresponding with people, realize the travel schedule of intelligent body, the time is included to each intelligent body stroke one Stab the path of sequence and sequence node.
4. the comprehensive traffic network operation situation towards highway and the railway network deduces emulation technology according to claim 3, its It is characterised by:Traffic flow simulation and dynamic shortest path module constitute the nucleus module of middle sight traffic simulation, traffic flow simulation Module is using section transmission, node-node transmission model, on the premise of relevant constraint is met, and realizes to vehicle inside section The simulation shifted between operation and section, section mode is related to the magnitude of traffic flow-density relations, and basic capacity is main By road segment classification and track quantity;Node metastasis model is related to left-hand rotation and straight trip Capacity Constraints, dynamically distributes platform institute Hourage of the traffic flow of simulation on section will be fed back in the dynamic shortest path model of next iteration, and then is intelligent body Specially designated path is selected to provide decision-making foundation, meanwhile, new pathway can be fed back to next iteration.
5. the comprehensive traffic network operation situation towards highway and the railway network deduces emulation technology according to claim 4, its It is characterised by:Operation evaluation module is estimated to operation result from section, passage, OD to, path, region, including efficiency, speed In terms of degree, reliability, energy consumption.
6. the comprehensive traffic network operation situation towards highway and the railway network deduces emulation technology according to claim 1, its It is characterised by, is predicted using polynomial trend and combine vehicle flowrate or passenger traffic volume point of the kalman filter method to Ingress node Do not estimated and predicted, method is as follows:
J=road-net node sequence numbers, j=1,2 ..., N;
τ=departure time interval sequence number, τ=1,2 ...;
The stage of k=rolling forecasts, k=1,2,3 ...;
L=each rolls the departure time intervening sequence number included in the stage;
D(j,τ)The demand influxs of=node j in time interval τ;
μ′(j,τ),μ″(j,τ),μ″′(j,τ)=demand disruption μ(j,τ)The single order of variable, second order, third order difference amount;
Estimate average;
Estimation average;
Zk=stage k state variable vector;
Yk=stage k measurement vector;
Hk=measurement vector YkWith state vector ZkMatching relationship matrix;
wk=stage k processing noise;
Rk=stage k measurement error;
Pk,k-1The Z of=stage k-1 predictionskState covariance battle array, i.e.,
Pk,k=stage k, the Z of estimationkState covariance battle array, i.e.,
Polynomial trend prediction is introduced, the prediction of τ+ζ difference vector is carried out, builds the function using m order polynomials as follows:
μ(j,τ+ζ)=b0+b1ζ+b2ζ2+…bpζp+…+bmζm
According to Taylor expansion anticipation function μ(j,τ+ζ)It is extended to μ(j,τ)Multistage difference combination:
The coefficient of the formula of original polynomial distribution can be expressed as:
Three rank multinomial models are taken to carry out μ(j,τ+ζ)Prediction be shown below:
We are defined on forecast period k, τ=kl, and such as each stage forecast Period Length is that l includes 3 time intervals, then In the k=3 stages, τ sequence number is 9;
State-transition matrix is defined as:
The predictor formula in k+1 stages is as follows:
Zk+1=AkZk
Introduce Kalman filtering algorithm to eliminate the influence that charge station's influx measured value error band comes, because this part is directly surveyed Value is identical vector value with estimate, therefore H=I;
Carried out first with the posterior estimate of previous stage, the prediction of current generation and the transmission of covariance:
Introduce afterwards after stage k measured value, carry out the calculating of gain:
The transmission of formation stages k posterior estimate and covariance matrix:
Pk,k=(I-KkHk)Pk,k-1
Call above-mentioned formula forecast period k+1 μ(j,τ+ζ), can be generated based on the predicted value and flow into value for the node of deduction D(j,τ), and the multistage difference matrix formed needed for subsequent prediction.
7. the comprehensive traffic network operation situation towards highway and the railway network deduces emulation technology according to claim 6, its It is characterised by, OD distributions that can be respectively to vehicle in freeway net by BP neural network model are predicted and to the railway network In go out administrative staff OD distribution be predicted, illustrate algorithm flow by taking the OD forecast of distribution of vehicle in freeway net as an example, calculate Method flow is as follows:
Step 1:Data decimation and standardization, for n charge tiny node in road network, choose before predicted time section m days The OD data of same period are simultaneously standardized according to formula 4.1;
Step 2:BP neural network is built, i=1 is made, for i-th of charge tiny node, the OD distribution proportions for choosing first day are made For input vector, choose the data of second day as desired output vector, the like constitute m-1 group training samples;
Step 3:Neutral net is trained, by repeatedly training, network is reached the training objective of setting;
Step 4:OD distribution proportions predict that, for the neutral net trained, the data of input prediction objective time interval the previous day are obtained To prediction distribution ratio, if going to step 5, step 2 is otherwise gone to, network parameter, re -training network is adjusted;
Step 5:Output prediction OD distribution proportion matrixes, the output matrix if i=n, otherwise i+1 go to step 2.
8. the comprehensive traffic network operation situation towards highway and the railway network deduces emulation technology according to claim 7, its It is characterised by, predicts obtain road grid traffic distribution proportion respectively with BP neural network algorithm and polynomial trend prediction algorithm After each entrance magnitude of traffic flow, two kinds are predicted the outcome and is merged, under the conditions of obtaining no event, the trip of prediction period on network Visitor or the OD matrixes of vehicle, by vehicle or passenger by node i to node j flow qijCan be by the flow O of website i entrancesi With with the distribution proportion r for going to j websitesijIt is determined that, under the conditions of without event, road network prediction OD matrixes are tried to achieve by equation below:
By the OD matrixes OD of above-mentioned formula respectively to vehiclevWith the OD matrixes OD of passengerpCalculated, so as to obtain in prediction Between vehicle under the conditions of Duan Shangwu events and the OD matrixes of passenger.
9. the comprehensive traffic network operation situation towards highway and the railway network deduces emulation technology according to claim 1, its It is characterised by, gives Travel Time Budget TTB (p), starting point Op, intelligent body p departure time τp, consider road capacity and Under the constraints of Passenger Flow service ability, maximize the space-time accessibility in system scope, the synthesis based on intelligent body The object function of the optimization problem of transportation network dynamically distributes model, is shown below:
That is, the optimization aim of model is to minimize the unreachable generalized cost of user.If intelligent body p can be in given time budget Arrive at, then inaccessible generalized cost is 0, otherwise, C (p) is by overall travel time of the intelligent body using path.
CN201710300803.2A 2017-07-11 2017-07-11 Emulation technology is deduced towards the comprehensive traffic network operation situation of highway and the railway network Active CN107103142B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710300803.2A CN107103142B (en) 2017-07-11 2017-07-11 Emulation technology is deduced towards the comprehensive traffic network operation situation of highway and the railway network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710300803.2A CN107103142B (en) 2017-07-11 2017-07-11 Emulation technology is deduced towards the comprehensive traffic network operation situation of highway and the railway network

Publications (2)

Publication Number Publication Date
CN107103142A true CN107103142A (en) 2017-08-29
CN107103142B CN107103142B (en) 2019-12-03

Family

ID=59657979

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710300803.2A Active CN107103142B (en) 2017-07-11 2017-07-11 Emulation technology is deduced towards the comprehensive traffic network operation situation of highway and the railway network

Country Status (1)

Country Link
CN (1) CN107103142B (en)

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108133291A (en) * 2017-12-22 2018-06-08 北京交通大学 A kind of correlation analysis between social event and passenger traffic demand
CN108446759A (en) * 2018-02-11 2018-08-24 陕西师范大学 Tourist flows amount prediction technique based on seasonal shot and long term memory network model
CN108491979A (en) * 2018-04-03 2018-09-04 中南大学 High-speed railway travelling time-dependent demand prediction technique based on equal strength distribution
CN108520353A (en) * 2018-04-09 2018-09-11 交通运输部公路科学研究所 A kind of composite communications transport network security risk point discrimination method
CN109118412A (en) * 2018-08-15 2019-01-01 北京交通大学 Urban rail transit network passenger flow on-line control system
CN109785625A (en) * 2019-02-18 2019-05-21 北京航空航天大学 A kind of " red green " region recognition and evaluation method of city road network operating status
CN110162931A (en) * 2019-06-21 2019-08-23 北京交通大学 Urban track traffic large-scale road network high-speed simulation system
CN110490365A (en) * 2019-07-12 2019-11-22 四川大学 A method of based on the pre- survey grid of multisource data fusion about vehicle order volume
CN110555223A (en) * 2018-06-04 2019-12-10 北京交通大学 High-speed train system reliability assessment method based on seepage theory
CN110599760A (en) * 2019-10-17 2019-12-20 东南大学 Travel behavior simulation method under multi-mode traffic network
CN111951553A (en) * 2020-08-17 2020-11-17 上海电科智能系统股份有限公司 Prediction method based on traffic big data platform and mesoscopic simulation model
CN112101610A (en) * 2020-07-29 2020-12-18 北京交通大学 Multi-mode network design method for maximizing reachable flow based on reference dependence theory
CN112201033A (en) * 2020-09-08 2021-01-08 同济大学 Road traffic operation active control strategy selection method
CN112598177A (en) * 2020-12-22 2021-04-02 北京市地铁运营有限公司 Online passenger flow prediction and simulation system for urban rail transit emergency
CN112633740A (en) * 2020-12-30 2021-04-09 北京城市轨道交通咨询有限公司 Urban rail transit network safety evaluation method
CN112753041A (en) * 2018-09-26 2021-05-04 科斯莫科技 Method for regulating a multi-mode transport network
CN113065215A (en) * 2021-03-22 2021-07-02 广州弘毅软件科技有限公司 Method and device for predicting operation situation of expressway
CN113129582A (en) * 2019-12-31 2021-07-16 阿里巴巴集团控股有限公司 Traffic state prediction method and device
CN113191528A (en) * 2021-04-02 2021-07-30 东南大学 Traffic state deduction method integrating macroscopic traffic simulation and short-term traffic prediction
CN113393044A (en) * 2021-06-22 2021-09-14 中远海运科技(北京)有限公司 Logistics resource emergency optimization system based on big data integration
CN116341882A (en) * 2023-05-29 2023-06-27 中铁第四勘察设计院集团有限公司 High-speed rail transfer passenger traffic analysis method and system based on generalized benefit model
CN112753041B (en) * 2018-09-26 2024-06-28 科斯莫科技 Method for controlling a multimode transport network

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103593535A (en) * 2013-11-22 2014-02-19 南京洛普股份有限公司 Urban traffic complex self-adaptive network parallel simulation system and method based on multi-scale integration
CN103793760A (en) * 2014-01-24 2014-05-14 同济大学 Multi-mode optimization method for comprehensive transportation junction interior passenger flow transfer distribution ratio
CN104217250A (en) * 2014-08-07 2014-12-17 北京市交通信息中心 Rail transit new line opening passenger flow prediction method based on historical data
CN105787858A (en) * 2016-03-29 2016-07-20 交通运输部公路科学研究所 Situation deduction method for expressway network

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103593535A (en) * 2013-11-22 2014-02-19 南京洛普股份有限公司 Urban traffic complex self-adaptive network parallel simulation system and method based on multi-scale integration
CN103793760A (en) * 2014-01-24 2014-05-14 同济大学 Multi-mode optimization method for comprehensive transportation junction interior passenger flow transfer distribution ratio
CN104217250A (en) * 2014-08-07 2014-12-17 北京市交通信息中心 Rail transit new line opening passenger flow prediction method based on historical data
CN105787858A (en) * 2016-03-29 2016-07-20 交通运输部公路科学研究所 Situation deduction method for expressway network

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
孟梦: ""组合出行模式下城市交通流分配模型与算法"", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 *
赵微: ""考虑路径约束的城市多模式交通流分配模型与算法研究"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108133291B (en) * 2017-12-22 2021-05-11 北京交通大学 Correlation analysis method between social events and passenger traffic demands
CN108133291A (en) * 2017-12-22 2018-06-08 北京交通大学 A kind of correlation analysis between social event and passenger traffic demand
CN108446759A (en) * 2018-02-11 2018-08-24 陕西师范大学 Tourist flows amount prediction technique based on seasonal shot and long term memory network model
CN108491979A (en) * 2018-04-03 2018-09-04 中南大学 High-speed railway travelling time-dependent demand prediction technique based on equal strength distribution
CN108520353A (en) * 2018-04-09 2018-09-11 交通运输部公路科学研究所 A kind of composite communications transport network security risk point discrimination method
CN110555223A (en) * 2018-06-04 2019-12-10 北京交通大学 High-speed train system reliability assessment method based on seepage theory
CN109118412A (en) * 2018-08-15 2019-01-01 北京交通大学 Urban rail transit network passenger flow on-line control system
CN109118412B (en) * 2018-08-15 2021-08-03 北京交通大学 Urban rail transit network passenger flow online control system
CN112753041A (en) * 2018-09-26 2021-05-04 科斯莫科技 Method for regulating a multi-mode transport network
CN112753041B (en) * 2018-09-26 2024-06-28 科斯莫科技 Method for controlling a multimode transport network
CN109785625B (en) * 2019-02-18 2020-09-08 北京航空航天大学 Method for identifying and evaluating red and green areas of urban road network operation state
CN109785625A (en) * 2019-02-18 2019-05-21 北京航空航天大学 A kind of " red green " region recognition and evaluation method of city road network operating status
CN110162931A (en) * 2019-06-21 2019-08-23 北京交通大学 Urban track traffic large-scale road network high-speed simulation system
CN110162931B (en) * 2019-06-21 2021-06-22 北京交通大学 Large-scale road network rapid simulation system for urban rail transit
CN110490365A (en) * 2019-07-12 2019-11-22 四川大学 A method of based on the pre- survey grid of multisource data fusion about vehicle order volume
CN110599760A (en) * 2019-10-17 2019-12-20 东南大学 Travel behavior simulation method under multi-mode traffic network
CN113129582A (en) * 2019-12-31 2021-07-16 阿里巴巴集团控股有限公司 Traffic state prediction method and device
CN113129582B (en) * 2019-12-31 2023-05-09 阿里巴巴集团控股有限公司 Traffic state prediction method and device
CN112101610B (en) * 2020-07-29 2023-10-17 北京交通大学 Multimode network design method for maximizing reachable flows based on reference dependency theory
CN112101610A (en) * 2020-07-29 2020-12-18 北京交通大学 Multi-mode network design method for maximizing reachable flow based on reference dependence theory
CN111951553A (en) * 2020-08-17 2020-11-17 上海电科智能系统股份有限公司 Prediction method based on traffic big data platform and mesoscopic simulation model
CN112201033A (en) * 2020-09-08 2021-01-08 同济大学 Road traffic operation active control strategy selection method
CN112201033B (en) * 2020-09-08 2022-06-07 同济大学 Road traffic operation active control strategy selection method
CN112598177A (en) * 2020-12-22 2021-04-02 北京市地铁运营有限公司 Online passenger flow prediction and simulation system for urban rail transit emergency
CN112598177B (en) * 2020-12-22 2023-10-24 北京市地铁运营有限公司 Urban rail transit emergency online passenger flow prediction and simulation system
CN112633740A (en) * 2020-12-30 2021-04-09 北京城市轨道交通咨询有限公司 Urban rail transit network safety evaluation method
CN112633740B (en) * 2020-12-30 2024-04-26 北京城市轨道交通咨询有限公司 Urban rail transit network safety evaluation method
CN113065215A (en) * 2021-03-22 2021-07-02 广州弘毅软件科技有限公司 Method and device for predicting operation situation of expressway
CN113191528B (en) * 2021-04-02 2022-02-22 东南大学 Traffic state deduction method integrating macroscopic traffic simulation and short-term traffic prediction
CN113191528A (en) * 2021-04-02 2021-07-30 东南大学 Traffic state deduction method integrating macroscopic traffic simulation and short-term traffic prediction
CN113393044A (en) * 2021-06-22 2021-09-14 中远海运科技(北京)有限公司 Logistics resource emergency optimization system based on big data integration
CN116341882A (en) * 2023-05-29 2023-06-27 中铁第四勘察设计院集团有限公司 High-speed rail transfer passenger traffic analysis method and system based on generalized benefit model

Also Published As

Publication number Publication date
CN107103142B (en) 2019-12-03

Similar Documents

Publication Publication Date Title
CN107103142B (en) Emulation technology is deduced towards the comprehensive traffic network operation situation of highway and the railway network
CN109308543B (en) Subway short-term passenger flow prediction method based on LS-SVM and real-time big data
CN111369181B (en) Train autonomous scheduling deep reinforcement learning method and device
Wang et al. Passenger-demands-oriented train scheduling for an urban rail transit network
CN103632212B (en) System and method for predicating time-varying user dynamic equilibrium network-evolved passenger flow
Chen et al. Uncertainty in urban mobility: Predicting waiting times for shared bicycles and parking lots
CN106485359A (en) A kind of urban track traffic section passenger flow estimation method based on train schedule
CN105787858B (en) Expressway network situation deduction method
CN102346964A (en) Real-time jam prediction and intelligent management system for road traffic network area
CN115527369B (en) Large passenger flow early warning and evacuation method under large-area delay condition of airport hub
Chandra et al. Cross-correlation analysis and multivariate prediction of spatial time series of freeway traffic speeds
CN107392389A (en) Taxi dispatching processing method based on ARIMA models
JP5376465B2 (en) Congestion situation prediction program, computer-readable recording medium recording congestion situation prediction program and congestion situation prediction apparatus, navigation program, computer-readable recording medium recording navigation program, and navigation apparatus
CN115204478A (en) Public traffic flow prediction method combining urban interest points and space-time causal relationship
Chen et al. A Deep Learning Model with Conv‐LSTM Networks for Subway Passenger Congestion Delay Prediction
Yedavalli et al. Microsimulation analysis for network traffic assignment (MANTA) at metropolitan-scale for agile transportation planning
CN110246332A (en) The real-time passenger flow method for monitoring of rail traffic and system based on multisource data fusion
Utku et al. Multi-layer perceptron based transfer passenger flow prediction in Istanbul transportation system
Kiani et al. A novel intelligent traffic recovery model for emergency vehicles based on context-aware reinforcement learning
CN114580751B (en) Method, system, storage medium and terminal for predicting evacuation time of passengers arriving at airport
Adhikari et al. An integrated solution approach for the time minimization evacuation planning problem.
Si et al. Urban transit assignment model based on augmented network with in-vehicle congestion and transfer congestion
Mutlu et al. COVID-19 transmission risk minimization at public transportation stops using Differential Evolution algorithm
Lv et al. An origin–destination passenger flow prediction system based on convolutional neural network and passenger source-based attention mechanism
CN113641749A (en) Population mobility network estimation method and system based on bipartite graph

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Zhang Fan

Inventor after: Zhao Li

Inventor after: Cui Wei

Inventor after: Zhang Jinjin

Inventor after: Wang Hua

Inventor after: Wang Tibin

Inventor after: Niu Shuyun

Inventor after: Zhang Jisheng

Inventor after: Zhou Xuesong

Inventor after: Li Honghai

Inventor after: Sun Xiaoliang

Inventor after: Liu Jianping

Inventor after: Zhang Li

Inventor after: Zhu Lili

Inventor before: Zhang Fan

Inventor before: Zhao Li

Inventor before: Cui Wei

Inventor before: Zhang Jinjin

Inventor before: Wang Hua

Inventor before: Wang Tibin

Inventor before: Zhang Jisheng

Inventor before: Niu Shuyun

Inventor before: Zhou Xuesong

Inventor before: Li Honghai

Inventor before: Sun Xiaoliang

Inventor before: Liu Jianping

Inventor before: Zhang Li

Inventor before: Zhu Lili

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant