WO2023134093A1 - Virtual subway train operation big data generation method based on interval speed limit and expert experience - Google Patents
Virtual subway train operation big data generation method based on interval speed limit and expert experience Download PDFInfo
- Publication number
- WO2023134093A1 WO2023134093A1 PCT/CN2022/092115 CN2022092115W WO2023134093A1 WO 2023134093 A1 WO2023134093 A1 WO 2023134093A1 CN 2022092115 W CN2022092115 W CN 2022092115W WO 2023134093 A1 WO2023134093 A1 WO 2023134093A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- speed
- subway train
- acceleration
- interval
- stage
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000012216 screening Methods 0.000 claims abstract description 5
- 230000001133 acceleration Effects 0.000 claims description 115
- 238000004364 calculation method Methods 0.000 claims description 20
- 101100129500 Caenorhabditis elegans max-2 gene Proteins 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 8
- 238000011160 research Methods 0.000 claims description 6
- 238000011017 operating method Methods 0.000 claims description 3
- 238000013079 data visualisation Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 13
- 101100083446 Danio rerio plekhh1 gene Proteins 0.000 description 10
- 238000004422 calculation algorithm Methods 0.000 description 7
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 241000766023 Coregonus oxyrinchus Species 0.000 description 1
- 101150032862 LEF-1 gene Proteins 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012552 review Methods 0.000 description 1
Images
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/40—Handling position reports or trackside vehicle data
Definitions
- the invention relates to the technical field of data processing, in particular to a method for generating big data of virtual subway train operation based on section speed limit and expert experience.
- the operation data acquisition methods of subway trains can be divided into two types: single operation curve generated by tracking simulation and manual driving data sampling.
- the simulation to generate a single running curve has the following research: through the four steps of reading data, data preprocessing, setting the calculation interval, calculating and outputting the results, the reference running curve is simulated and generated; the genetic algorithm is used to find out the position and location of the key points of the subway train. Acceleration of subway trains, and then use the knowledge of subway train kinematics to obtain the speed curve; use non-parametric iterative learning method to optimize the speed curve of a single subway train; use PSO-CS algorithm to provide the optimal control target curve for subway trains.
- the generation process of this method of obtaining subway train operation data is relatively complicated, the optimization time is long, and the generated curve can only adapt to one type of operation time.
- the purpose of the present invention is to provide a virtual subway train operation big data generation method based on interval speed limit and expert experience, generate a large amount of virtual data based on reality, and screen according to the needs of researchers to obtain curves that meet the requirements, and perform Isochronous sampling is used to obtain the tracking data set of the curve, which is convenient for the study of train operation.
- a method for generating big data of virtual subway train operation based on section speed limit and expert experience which includes the following steps:
- Step 1 set the parameter interval of the virtual subway train according to expert experience and data, the parameter interval includes acceleration interval, speed interval, station distance interval and speed change distance interval;
- Step 2 using the five-stage operation method to segment the operation curve of the virtual subway train, and set the operation parameters for each stage respectively;
- Step 3 using the operating parameters of each stage of the set five-stage operating method to calculate the operating data of several groups of virtual subway trains, and using drawing tools to generate corresponding operating curves;
- Step 4 set screening conditions according to research requirements to obtain corresponding operating curves
- Step 5 Sampling curve data at equal time intervals for each selected curve, and outputting the state of the train at each time point.
- step 1 specifically includes the following steps:
- Step 1-1 calculate and obtain the actual subway train acceleration interval based on the actual subway train data, configure the acceleration interval of the virtual subway train to include the actual subway train acceleration interval; that is, the actual subway train acceleration interval is a subset of the acceleration interval of the virtual subway train;
- Step 1-2 calculate the speed interval of the actual subway train operation based on the actual subway train data, configure the speed interval of the virtual subway train to include the speed interval of the actual subway train operation; that is, the speed interval of the actual subway train operation is the speed interval of the virtual subway train operation a subset of velocity intervals;
- Step 1-3 based on the interval between two adjacent stations of the actual subway train line, configure the interval between two adjacent stations of the virtual subway train line to include the interval between two adjacent stations of the actual subway train line; that is, the actual subway train
- the interval between two adjacent stations of the line is a subset of the interval between two adjacent stations of the virtual subway train line;
- steps 1-4 the running process of the virtual subway train is divided into four stages of speed change and one stage of constant speed.
- the speed change distance is set to 100-600m, and the step length is 50m to form a speed change distance interval.
- step 2 divides the operation curve of the virtual subway train into four stages of variable speed process and one stage of constant speed process, and the operating parameters of each stage are as follows:
- the acceleration time t1 of the first shifting phase is calculated from the first acceleration limit speed v max1 and the first shifting acceleration a max1 , as shown in formula (1):
- the constant-speed travel time t1 of the first shift stage is calculated from the remaining preset shift distance Lef 1 of the first shift and the first acceleration speed limit v max1 , as shown in formula (2):
- the remaining preset shift distance Lef 1 of the first shift is calculated from the preset shift distance dis and the first acceleration distance s 1 , as shown in formula (3):
- the first acceleration distance s 1 is calculated from the acceleration time t 1 of the first shift stage and the first acceleration speed limit v max1 , as shown in formula (4):
- the second shift stage is the second shift stage
- the calculation of the running time of the second shifting stage is the same as the first stage, which can be calculated by formula (6) ⁇ formula (10).
- the calculation variables involved are the acceleration time t 2 of the first shifting stage, the first acceleration speed limit v max1 , the second acceleration speed limit v max2 , the second shifting acceleration a max2 , the second shifting speed
- the constant speed travel time t 2 the remaining preset shift distance Lef 2 for the second shift, the preset shift distance dis, and the second acceleration distance s 2 .
- the running time runtime of the subway train at the stage of constant speed running is calculated by the second acceleration speed limit vmax2 of the subway train and the runtance of the constant speed running distance, as shown in formula (11):
- the third shift stage is the third shift stage.
- the running time of the third shifting stage can be calculated by formula (12) ⁇ formula (16).
- the calculation variables involved are the deceleration time t 3 of the third shifting stage, the third deceleration speed limit v max3 , the second acceleration speed limit v max2 , the third shifting acceleration a max3 , the third shifting speed limit The constant speed travel time th 3 , the remaining preset shift distance Lef 3 for the third shift, the preset shift distance dis, and the third acceleration distance s 3 .
- the running time of the fourth shifting stage can be calculated by formula (17) ⁇ formula (21).
- the calculation variables involved are the deceleration time t 4 of the fourth shift stage, the third deceleration speed limit v max3 , the fourth shift acceleration a max4 , the constant speed travel time th 4 of the fourth shift stage, and the fourth The remaining preset shifting distance Lef 4 , the preset shifting distance dis , and the fourth acceleration distance s 4 .
- the running time alltime in the fastest running state is calculated from the running time T 1 , T 2 , T 3 , T 4 of the four shifting stages and the running time runtime at a constant speed, as shown in formula (22).
- step 3 the operating curve is drawn and formed with the curve corresponding to the fastest operating state as the upper bound and the curve corresponding to the slowest operating state as the lower bound.
- the data sampling interval in step 5 is 0.1 second.
- the state of the train at each time point in step 5 includes current position acceleration m/s 2 , current position elapsed running time s, current position running speed m/s, current position running distance m, current position speed limit m/s s, the time to maintain the current driving state s, the remaining shifting distance at the current position m, the remaining shifting range at the current position m/s, the remaining time at the current position at the station s, and the remaining distance at the current position at the station m.
- the present invention adopts the above technical scheme, based on human-machine hybrid intelligence, learns the advanced idea of AlphaZero, abandons the traditional method of collecting subway train operation data, generates a large amount of data according to expert experience and curve generation algorithm, and then performs screening and drawing. That is to say, the subway train driving is regarded as playing Go in one-dimensional space, and the required parameters are proposed according to the subway train speed curve generation algorithm, and the parameter range is set in combination with expert experience, so as to greatly reduce the order of magnitude of the subway train driving data and meet various running times. and interval requirements.
- the various running times covered by the virtual data can be observed by drawing the frequency distribution map of the running time of subway trains, which can meet the data requirements of different situations, and is more conducive to the research of intelligent driving algorithms for subway trains than traditional data.
- Fig. 1 is the schematic flow sheet of the present invention's virtual subway train operation big data generation method based on section speed limit and expert experience;
- Fig. 2 is five section schematic diagrams of virtual subway train operation
- Fig. 3 is the fastest running state and the slowest running state schematic diagram of virtual subway train operation
- Fig. 4 is the relation schematic diagram of subway train speed-distance
- Fig. 5 is a schematic diagram of frequency distribution of subway train running time
- Figure 6 is a schematic diagram of the operating curve under the condition of a speed limit of 80km/h;
- Figure 7 is a schematic diagram of the running time under the condition of a speed limit of 80km/h;
- Figure 8 is a schematic diagram of the operating curve under the condition of a speed limit of 100km/h;
- Figure 9 is a schematic diagram of the fastest running time under the condition of a speed limit of 100km/h;
- Figure 10 is a schematic diagram of the operation curve under the condition of a station distance of 2800m;
- Figure 11 is a schematic diagram of the running time under the condition of a station distance of 2800m and a speed limit of 100km/h;
- Figure 12 Schematic diagram of the operation curve under the condition of a station distance of 700m
- Figure 13 Schematic diagram of the running time with a station distance of 700m and a speed limit of 100km/h;
- Figure 14 Schematic diagram of positioning a single curve.
- the most critical idea of the present invention is to set the parameter range according to expert experience and data, then generate virtual data through the train operation method set in the present invention, and then filter the data set according to the researcher's requirements to obtain a small number of curves. Perform isochronous tracking on a small number of curves obtained, and export the running state data set of the curves.
- the present invention discloses a virtual subway train operation big data generation method based on section speed limit and expert experience, which includes the following steps:
- Step 1 set the parameter interval of the virtual subway train according to expert experience and data, the parameter interval includes the acceleration interval, the speed interval, the station distance interval and the speed change distance interval; further, the setting of the parameter interval in step 1 specifically includes the following steps:
- Step 1-1 calculate and obtain the actual subway train acceleration interval based on the actual subway train data, configure the acceleration interval of the virtual subway train to include the actual subway train acceleration interval; that is, the actual subway train acceleration interval is a subset of the acceleration interval of the virtual subway train;
- the present invention has inquired about the operation data of Fuzhou subway and obtained the following data: when the running speed of Fuzhou subway is 0-40km/h, the average acceleration is greater than 0.83m/h s2, when the running speed is 0-80km/h, the average acceleration is greater than 0.5m/s2. Therefore, the acceleration interval of the subway train speed curve data set generation algorithm should be greater than the actual subway train acceleration interval, that is, greater than the maximum acceleration and less than the minimum acceleration.
- the specific data is set to 0.3-1.3m/s2, and the step size is set to 0.1m/s2.
- Step 1-2 calculate the speed interval of the actual subway train operation based on the actual subway train data, configure the speed interval of the virtual subway train to include the speed interval of the actual subway train operation; that is, the speed interval of the actual subway train operation is the speed interval of the virtual subway train operation a subset of velocity intervals;
- Fuzhou Metro Line 1 As of December 2020, the total length of Fuzhou Metro Line 1 is 29.582km, all of which are underground; a total of 25 stations are set up, all of which are underground; 6-section B-type The maximum operating speed of subway trains is 80km/h[13].
- Fuzhou Metro Line 2 is initially equipped with 31 subway vehicles from Quanzhou CRRC Tangche Company. The subway trains are composed of 4 trains, 2 trains and 6 trains, with a maximum passenger capacity of 1 880 people , with a maximum operating speed of 80km/h[14].
- the speed range of the virtual data is set to be greater than the actual running speed range, that is, the virtual fastest running speed should be greater than the actual fastest running speed of 100km/h.
- the slowest running speed can not be too low, which can be set to 20km/h as an embodiment of the present invention.
- the operating speed range is set to 20-100km/h
- the step size is set to 5km/h.
- Step 1-3 based on the interval between two adjacent stations of the actual subway train line, configure the interval between two adjacent stations of the virtual subway train line to include the interval between two adjacent stations of the actual subway train line; that is, the actual subway train
- the interval between two adjacent stations of the line is a subset of the interval between two adjacent stations of the virtual subway train line;
- the starting point of Fuzhou Metro Line 1 Project (Phase I) is Xiangfeng Station, and the terminal station is Fuzhou South railway Station.
- the spacing is 1.202km; as of April 2019, Fuzhou Metro Line 2 has a total of 22 stations, all of which are underground stations.
- the maximum station spacing is 2.827km (from Houting Station to Juyuanzhou Station), and the minimum station spacing is 0.745km (Juyuanzhou Station to Hongwan Station), the average station distance is 1.392km.
- the distance between stations is set to 700-3 000m, which is smaller than the minimum station distance and greater than the maximum station distance to include the possibility of all station distances, and the step length is set to 100m; because the subway train speed generation algorithm includes four stages of variable speed process and one stage of constant speed process , so set the speed change distance as 100-600m, and set the step length as 50m.
- steps 1-4 the running process of the virtual subway train is divided into four stages of speed change and one stage of constant speed.
- the speed change distance is set to 100-600m, and the step length is 50m to form a speed change distance interval.
- Step 2 using the five-stage operation method to segment the operation curve of the virtual subway train, and set the operation parameters for each stage respectively;
- the subway train in the ideal operating state, is regarded as a rigid body, assuming that it has no friction, no air resistance, and no force on the carriage, and a five-stage subway train operation method is proposed.
- the specific implementation method is shown in Figure 2.
- the bold black dotted line in the figure is the speed limit curve of the subway train; the vertical dotted line is the division line of the operation phase; the gray point is the turning point of the operation phase.
- the curve is divided into 5 segments (divided by vertical dotted lines in Figure 1).
- the subway train accelerates to the preset first acceleration limit speed vmax1 with the preset first shifting acceleration amax1, and then travels at a constant speed to complete the first shifting remaining preset acceleration distance Lef1;
- the second shifting stage the subway train accelerates with the preset second acceleration amax2 to the preset second acceleration limit speed vmax2 of the subway train, and then moves at a constant speed, and the remaining preset acceleration distance Lef2 is left after the second speed change;
- the preset second acceleration limit speed vmax2 of the preset subway train runs at a constant speed and finishes the preset running distance runtance at a constant speed; the last two sections are decelerated running, and in the third shifting stage, the subway train first runs through the third shifting and the remaining preset shifting distance Lef3, Then decelerate with the preset third speed change acceleration amax3, so that the speed of the subway train is reduced to the preset first decel
- the five-stage operation method in step 2 divides the operation curve of the virtual subway train into the first shifting stage, the second shifting stage, the uniform speed running stage, the third shifting stage, and the fourth shifting stage, each The stage operation parameters are as follows:
- the acceleration time t1 of the first shifting phase is calculated from the first acceleration limit speed v max1 and the first shifting acceleration a max1 , as shown in formula (1):
- the constant-speed travel time th1 of the first shift stage is calculated from the remaining preset shift distance Lef 1 of the first shift and the first acceleration speed limit v max1 , as shown in formula (2):
- the remaining preset shift distance Lef 1 of the first shift is calculated from the preset shift distance dis and the first acceleration distance s 1 , as shown in formula (3):
- the first acceleration distance s 1 is calculated from the acceleration time t 1 of the first shift stage and the first acceleration speed limit v max1 , as shown in formula (4):
- T 1 t 1 +th 1 (5)
- the second shift stage is the second shift stage
- the calculation of the running time of the second shifting stage is the same as the first stage, which can be calculated by formula (6) ⁇ formula (10).
- the calculation variables involved are the acceleration time t 2 of the first shifting stage, the first acceleration speed limit v max1 , the second acceleration speed limit v max2 , the second shifting acceleration a max2 , the second shifting speed
- the constant speed travel time th 2 the remaining preset shift distance Lef 2 for the second shift, the preset shift distance dis, and the second acceleration distance s 2 .
- the running time runtime of the subway train at the stage of constant speed running is calculated by the second acceleration speed limit vmax2 of the subway train and the runtance of the constant speed running distance, as shown in formula (11):
- the third shift stage is the third shift stage.
- the running time of the third shifting stage can be calculated by formula (12) ⁇ formula (16).
- the calculation variables involved are the deceleration time t 3 of the third shifting stage, the third deceleration speed limit v max3 , the second acceleration speed limit v max2 , the third shifting acceleration a max3 , the third shifting speed limit The constant speed travel time th 3 , the remaining preset shift distance Lef 3 for the third shift, the preset shift distance dis, and the third acceleration distance s 3 .
- the running time of the fourth shifting stage can be calculated by formula (17) ⁇ formula (21).
- the calculation variables involved are the deceleration time t 4 of the fourth shift stage, the third deceleration speed limit v max3 , the fourth shift acceleration a max4 , the constant speed travel time th 4 of the fourth shift stage, and the fourth The remaining preset shifting distance Lef 4 , the preset shifting distance dis , and the fourth acceleration distance s 4 .
- the running time alltime in the fastest running state is calculated from the running time T 1 , T 2 , T 3 , T 4 of the four shifting stages and the running time runtime at a constant speed, as shown in formula (22).
- Step 3 using the operating parameters of each stage of the set five-stage operating method to calculate the operating data of several groups of virtual subway trains, and using drawing tools to generate corresponding operating curves;
- the first running state is the fastest running close to the speed limit curve of the subway train state, in this state, the first acceleration limit speed vmax1 and the first deceleration limit speed vmax3 of the subway train are set as the minimum speed limit Emin of the subway train, and the second acceleration limit speed vmax2 of the subway train is set as the minimum speed limit Emax of the subway train; the speed change stage The accelerations amax1, amax2, amax3, and amax4 are all the fastest accelerations; the second running state is to set the first acceleration limit speed vmax1 of the preset subway train and the first deceleration limit speed vmax3 to half of the preset minimum speed limit Emin of the subway train.
- the second deceleration limit speed vmax2 of the preset subway train is set to the minimum speed limit Emin of the preset subway train; the accelerations amax1, amax2, amax3, and amax4 of the speed change stage are all the slowest accelerations, and the slowest running state is obtained.
- These two operating states are the upper boundary and lower boundary of the generated curve, and other generated curves are in between, as shown in Figure 4.
- the data is classified, and the subway train operation data set can be obtained, with a total of 312 800 0 data, which can be realized by Python programming and exported as a CSV file.
- the column variables include the minimum speed limit, maximum speed limit, acceleration, preset acceleration Distance, two-station distance, four-segment acceleration distance, running time and other parameters.
- 209 005 0 pieces of data were obtained after eliminating the curves that did not conform to the actual situation (such as running time greater than 600s, not conforming to the preset running model).
- the frequency distribution of running time of these data is shown in Figure 5.
- Step 4 set screening conditions according to research requirements to obtain corresponding operating curves
- the present invention set the fastest running speed of Fuzhou Metro to 80km/h, the fastest running speed of the data set to 100km/h, the maximum station distance of Fuzhou Metro Line 2 to 2800m, and the shortest distance between the two stations of the data set. Taking the spacing of 700m as an example to screen, the following results can be obtained.
- Step 5 Sampling curve data at equal time intervals for each selected curve, and outputting the state of the train at each time point.
- the state of the train at each time point in step 5 includes current position acceleration m/s 2 , current position elapsed running time s, current position running speed m/s, current position running distance m, current position speed limit m/s s, the time to maintain the current driving state s, the remaining shifting distance at the current position m, the remaining shifting range at the current position m/s, the remaining time at the current position at the station s, and the remaining distance at the current position at the station m.
- the invention makes a breakthrough in replacing the traditional method of obtaining subway train operation data by artificially generating a large number of subway operation curves. Specifically, by investigating the actual data of subway operation and consulting relevant experts to obtain the limited parameters of the virtual data such as the traction interval, speed interval, station distance interval, and speed change interval, so that the generated virtual data has a small range and high reliability; and then according to the dynamics knowledge Assume the five-stage running mode of the subway train, so as to generate the virtual speed-distance image of the train; then use python to carry out programming and data visualization to obtain a large amount of virtual subway train running data required by the present invention. For the generated large amount of virtual data, the present invention screens them according to the conditions, and obtains the partial dotted lines required by the researcher.
- the present invention tracks it in real time, and obtains the running state of the single curve subway train every 0.1s, specifically including: current position acceleration, current position elapsed running time, current position running speed, current The running distance at the current location, the speed limit at the current location, the time to maintain the current driving state, the remaining shifting distance at the current location, the remaining shifting range at the current location, the remaining time at the current location to the station, and the remaining distance at the current location to the station.
- the present invention can understand the The real-time running status of the train is convenient for researchers to further study the subway train.
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Train Traffic Observation, Control, And Security (AREA)
Abstract
Description
本发明涉及数据处理技术领域,尤其涉及基于区间限速和专家经验的虚拟地铁列车运行大数据生成方法。The invention relates to the technical field of data processing, in particular to a method for generating big data of virtual subway train operation based on section speed limit and expert experience.
近年来,地铁列车的运行数据获取方式可分为跟踪仿真生成的单条运行曲线和人工驾驶数据的采样两种。其中,仿真生成单条运行曲线有以下研究:通过读入数据、数据预处理、设定计算区间、计算并输出结果4个步骤仿真生成参考运行曲线;采用遗传算法找出地铁列车关键点的位置和地铁列车加速度,然后利用地铁列车运动学知识得到速度曲线;采用非参数化的迭代学习方法进行优化单条地铁列车速度曲线;采用PSO-CS算法为地铁列车提供最优控制目标曲线。但这种获取地铁列车运行数据的方法生成过程比较复杂,优化时间长,并且生成的曲线只能适应一种运行时间。In recent years, the operation data acquisition methods of subway trains can be divided into two types: single operation curve generated by tracking simulation and manual driving data sampling. Among them, the simulation to generate a single running curve has the following research: through the four steps of reading data, data preprocessing, setting the calculation interval, calculating and outputting the results, the reference running curve is simulated and generated; the genetic algorithm is used to find out the position and location of the key points of the subway train. Acceleration of subway trains, and then use the knowledge of subway train kinematics to obtain the speed curve; use non-parametric iterative learning method to optimize the speed curve of a single subway train; use PSO-CS algorithm to provide the optimal control target curve for subway trains. However, the generation process of this method of obtaining subway train operation data is relatively complicated, the optimization time is long, and the generated curve can only adapt to one type of operation time.
另外,有研究从人工驾驶专家的运行曲线中获取数据进行挖掘,每0.5s对地铁列车的状态进行采样,得到地铁列车的速度-位移曲线,这种方法费时费力且效率不高,得到的数据也只能适用于该段运行区间以及所设定的运行时间,通用性不强。In addition, some studies have obtained data mining from the operating curves of human driving experts, and sampled the state of subway trains every 0.5s to obtain the speed-displacement curves of subway trains. This method is time-consuming, laborious and inefficient, and the obtained data It can only be applied to this section of the running interval and the set running time, and the versatility is not strong.
发明内容Contents of the invention
本发明的目的在于提供基于区间限速和专家经验的虚拟地铁列车运行大数据生成方法,生成大量基于现实的虚拟数据,并根据研究者的需求进行筛选,得到符合要求的曲线,并对其进行等时采样,得到曲线的跟踪数据集,便于列车运行研究。The purpose of the present invention is to provide a virtual subway train operation big data generation method based on interval speed limit and expert experience, generate a large amount of virtual data based on reality, and screen according to the needs of researchers to obtain curves that meet the requirements, and perform Isochronous sampling is used to obtain the tracking data set of the curve, which is convenient for the study of train operation.
本发明采用的技术方案是:The technical scheme adopted in the present invention is:
基于区间限速和专家经验的虚拟地铁列车运行大数据生成方法,其包括以下步骤:A method for generating big data of virtual subway train operation based on section speed limit and expert experience, which includes the following steps:
步骤1,根据专家经验和资料设置虚拟地铁列车的参数区间,参数区间包括加速度区间、速度区间、车站间距区间和变速距离区间;
步骤2,利用五段式运行法将虚拟地铁列车的运行曲线进行分段,并分别设定每个阶段运行参数;
步骤3,利用设定的五段式运行法的每个阶段运行参数计算得到若干组虚拟地铁列车的运行数据,并分别利用绘图工具生成对应运行曲线;
步骤4,根据研究需求设定筛选条件得到对应的运行曲线;
步骤5,对每条筛选出的曲线进行等时间隔的曲线数据采样,并输出每个时间点列车所处状态。Step 5: Sampling curve data at equal time intervals for each selected curve, and outputting the state of the train at each time point.
进一步地,步骤1的参数区间的设置具体包括以下步骤:Further, the setting of the parameter interval in
步骤1-1,基于实际地铁列车资料计算得到实际地铁列车加速度区间,配置虚拟地铁列车的加速度区间包含实际地铁列车加速度区间;即实际地铁列车加速度区间为虚拟地铁列车的加速度区间的子集;Step 1-1, calculate and obtain the actual subway train acceleration interval based on the actual subway train data, configure the acceleration interval of the virtual subway train to include the actual subway train acceleration interval; that is, the actual subway train acceleration interval is a subset of the acceleration interval of the virtual subway train;
步骤1-2,基于实际地铁列车资料计算得到实际地铁列车运行的速度区间,配置虚拟地铁列车的速度区间包含实际地铁列车运行的速度区间;即实际地铁列车运行的速度区间为虚拟地铁列车运行的速度区间的子集;Step 1-2, calculate the speed interval of the actual subway train operation based on the actual subway train data, configure the speed interval of the virtual subway train to include the speed interval of the actual subway train operation; that is, the speed interval of the actual subway train operation is the speed interval of the virtual subway train operation a subset of velocity intervals;
步骤1-3,基于实际地铁列车线路的相邻两个车站间距区间,配置虚拟地铁列车线路的相邻两个车站间距区间包含实际地铁列车线路的相邻两个车站间距区间;即实际地铁列车线路的相邻两个车站间距区间为虚拟地铁列车线路的相邻两个车站间距区间的子集;Step 1-3, based on the interval between two adjacent stations of the actual subway train line, configure the interval between two adjacent stations of the virtual subway train line to include the interval between two adjacent stations of the actual subway train line; that is, the actual subway train The interval between two adjacent stations of the line is a subset of the interval between two adjacent stations of the virtual subway train line;
步骤1-4,将虚拟地铁列车的运行过程分为四段变速过程和一段匀速过程,设定变速距离为100-600m,步长为50m,以形成变速距离区间。In steps 1-4, the running process of the virtual subway train is divided into four stages of speed change and one stage of constant speed. The speed change distance is set to 100-600m, and the step length is 50m to form a speed change distance interval.
进一步地,步骤2中五段式运行法将虚拟地铁列车的运行曲线分为四段变速过程和一段匀速过程,每个阶段运行参数如下:Further, the five-stage operation method in
第一次变速阶段:First shift stage:
首先,第一次变速阶段的加速时间t 1由第一次加速限速v max1和第一次变速加速度a max1计算得到,如式(1)所示: First, the acceleration time t1 of the first shifting phase is calculated from the first acceleration limit speed v max1 and the first shifting acceleration a max1 , as shown in formula (1):
第一次变速阶段的匀速行驶时间t 1由第一次变速剩余预设变速距离Lef 1和第一次加速限速v max1计算得到,如式(2)所示: The constant-speed travel time t1 of the first shift stage is calculated from the remaining preset shift distance Lef 1 of the first shift and the first acceleration speed limit v max1 , as shown in formula (2):
其中第一次变速剩余预设变速距离Lef 1由预设变速距离dis和第一次加速距离s 1计算得到,如式(3)所示: The remaining preset shift distance Lef 1 of the first shift is calculated from the preset shift distance dis and the first acceleration distance s 1 , as shown in formula (3):
Lef 1=dis-s 1 (3) Lef 1 = dis-s 1 (3)
第一次加速距离s 1由第一次变速阶段的加速时间t 1和第一次加速限速v max1计算得到,如式(4)所示: The first acceleration distance s 1 is calculated from the acceleration time t 1 of the first shift stage and the first acceleration speed limit v max1 , as shown in formula (4):
则第一阶段的总运行时间T 1由加速时间t 1和匀速行驶时间t 1计算得到,如式(5)所示: Then the total running time T1 of the first stage is calculated from the acceleration time t1 and the constant speed running time t1 , as shown in formula (5):
T 1=t 1+t 1 (5) T 1 =t 1 +t 1 (5)
第二次变速阶段:The second shift stage:
第二次变速阶段运行时间的计算与第一阶段原理相同,即可通过式(6)~式(10)计算得出。其中涉及的计算变量分别为第一次变速阶段的加速时间t 2、第一次加速限速v max1、第二次加速限速v max2、第二次变速加速度a max2、第二次变速阶段的匀速行驶时间t 2、第二次变速剩余预设变速距离Lef 2、预设变速距离dis、第二次加速距离s 2。 The calculation of the running time of the second shifting stage is the same as the first stage, which can be calculated by formula (6) ~ formula (10). The calculation variables involved are the acceleration time t 2 of the first shifting stage, the first acceleration speed limit v max1 , the second acceleration speed limit v max2 , the second shifting acceleration a max2 , the second shifting speed The constant speed travel time t 2 , the remaining preset shift distance Lef 2 for the second shift, the preset shift distance dis, and the second acceleration distance s 2 .
Lef 2=dis-s 2 (8) Lef 2 = dis-s 2 (8)
T 2=t 2+t 2 (10) T 2 =t 2 +t 2 (10)
匀速行驶阶段:Uniform speed driving stage:
匀速行驶阶段地铁列车的运行时间runtime由地铁列车第二次加速限速v max2和匀速行驶距离runtance计算得到,如式(11)所示: The running time runtime of the subway train at the stage of constant speed running is calculated by the second acceleration speed limit vmax2 of the subway train and the runtance of the constant speed running distance, as shown in formula (11):
第三次变速阶段:The third shift stage:
第三次变速阶段运行时间可通过式(12)~式(16)计算得出。其中涉及的计算变量分别为第三次变速阶段的减速时间t 3、第三次减速限速v max3、第二次加速限速v max2、第三次变速加速度a max3、第三次变速阶段的匀速行驶时间th 3、第三次变速剩余预设变速距离Lef 3、预设变速距离dis、第三次加速距离s 3。 The running time of the third shifting stage can be calculated by formula (12) ~ formula (16). The calculation variables involved are the deceleration time t 3 of the third shifting stage, the third deceleration speed limit v max3 , the second acceleration speed limit v max2 , the third shifting acceleration a max3 , the third shifting speed limit The constant speed travel time th 3 , the remaining preset shift distance Lef 3 for the third shift, the preset shift distance dis, and the third acceleration distance s 3 .
Lef 3=dis-s 3 (14) Lef 3 = dis-s 3 (14)
T 3=t 3+t 3 (16) T 3 =t 3 +t 3 (16)
第四次变速阶段fourth shift stage
第四次变速阶段运行时间可通过式(17)~式(21)计算得出。其中涉及的计算变量分别为第四次变速阶段的减速时间t 4、第三次减速限速v max3、第四次变速加速度a max4、第四次变速阶段的匀速行驶时间th 4、第四次变速剩余预设变速距离Lef 4、预设变速距离dis、第四次加 速距离s 4。 The running time of the fourth shifting stage can be calculated by formula (17) ~ formula (21). The calculation variables involved are the deceleration time t 4 of the fourth shift stage, the third deceleration speed limit v max3 , the fourth shift acceleration a max4 , the constant speed travel time th 4 of the fourth shift stage, and the fourth The remaining preset shifting distance Lef 4 , the preset shifting distance dis , and the fourth acceleration distance s 4 .
Lef 4=dis-s 4 (19) Lef 4 = dis-s 4 (19)
T 4=t 4+t 4 (21) T 4 =t 4 +t 4 (21)
最快运行状态运行时间计算:Calculation of running time in the fastest running state:
最快运行状态下的运行时间alltime由四次变速阶段运行时间T 1、T 2、T 3、T 4和匀速运行时间runtime计算得到,如式(22)。 The running time alltime in the fastest running state is calculated from the running time T 1 , T 2 , T 3 , T 4 of the four shifting stages and the running time runtime at a constant speed, as shown in formula (22).
alltime=T 1+T 2+T 3+T 4+runtime (22). alltime=T 1 +T 2 +T 3 +T 4 +runtime (22).
进一步地,步骤3中绘制形成运行曲线以最快运行状态对应曲线为上界并以最慢运行状态对应曲线作为下界。Further, in
进一步地,步骤5中数据采样的间隔为0.1秒。Further, the data sampling interval in step 5 is 0.1 second.
进一步地,步骤5中每个时间点列车所处状态包括当前位置加速度m/s 2,当前位置已运行时间s、当前位置运行速度m/s、当前位置运行距离m、当前位置限速m/s、保持当前行驶状态时间s、当前位置剩余变速距离m、当前位置剩余变速范围m/s、当前位置到站剩余时间s和当前位置到站剩余距离m。 Further, the state of the train at each time point in step 5 includes current position acceleration m/s 2 , current position elapsed running time s, current position running speed m/s, current position running distance m, current position speed limit m/s s, the time to maintain the current driving state s, the remaining shifting distance at the current position m, the remaining shifting range at the current position m/s, the remaining time at the current position at the station s, and the remaining distance at the current position at the station m.
本发明采用以上技术方案,基于人机混合智能,学习AlphaZero的先进思路,摒弃传统地铁列车运行数据采集方式,根据专家经验和曲线生成算法生成大量数据,再进行筛选、绘图。即把地铁列车驾驶看作在一维空间下围棋,根据地铁列车速度曲线生成算法提出所需参数,并结合专家经验设置参数范围,大幅度缩小地铁列车驾驶数据的数量级,并且满足各种运行时间和区间情况的需求。通过绘制地铁列车运行时间频率分布图观测到该虚拟数据覆盖的各种运行时间,可满足不同情况的数据需求,比传统数据更有利于地铁列车智能驾驶算法的研究。The present invention adopts the above technical scheme, based on human-machine hybrid intelligence, learns the advanced idea of AlphaZero, abandons the traditional method of collecting subway train operation data, generates a large amount of data according to expert experience and curve generation algorithm, and then performs screening and drawing. That is to say, the subway train driving is regarded as playing Go in one-dimensional space, and the required parameters are proposed according to the subway train speed curve generation algorithm, and the parameter range is set in combination with expert experience, so as to greatly reduce the order of magnitude of the subway train driving data and meet various running times. and interval requirements. The various running times covered by the virtual data can be observed by drawing the frequency distribution map of the running time of subway trains, which can meet the data requirements of different situations, and is more conducive to the research of intelligent driving algorithms for subway trains than traditional data.
以下结合附图和具体实施方式对本发明做进一步详细说明;The present invention will be described in further detail below in conjunction with accompanying drawing and specific embodiment;
图1为本发明基于区间限速和专家经验的虚拟地铁列车运行大数据生成方法的流程示意图;Fig. 1 is the schematic flow sheet of the present invention's virtual subway train operation big data generation method based on section speed limit and expert experience;
图2为虚拟地铁列车运行的五段示意图;Fig. 2 is five section schematic diagrams of virtual subway train operation;
图3为虚拟地铁列车运行的最快运行状态和最慢运行状态示意图;Fig. 3 is the fastest running state and the slowest running state schematic diagram of virtual subway train operation;
图4为地铁列车速度-距离的关系示意图;Fig. 4 is the relation schematic diagram of subway train speed-distance;
图5为地铁列车运行时间频率分布示意图;Fig. 5 is a schematic diagram of frequency distribution of subway train running time;
图6为限速80km/h条件下的运行曲线示意图;Figure 6 is a schematic diagram of the operating curve under the condition of a speed limit of 80km/h;
图7为限速80km/h条件下的运行时间示意图;Figure 7 is a schematic diagram of the running time under the condition of a speed limit of 80km/h;
图8为限速100km/h条件下的运行曲线示意图;Figure 8 is a schematic diagram of the operating curve under the condition of a speed limit of 100km/h;
图9为限速100km/h条件下的最快运行时间示意图;Figure 9 is a schematic diagram of the fastest running time under the condition of a speed limit of 100km/h;
图10为站距2800m条件下的运行曲线示意图;Figure 10 is a schematic diagram of the operation curve under the condition of a station distance of 2800m;
图11为站距2800m且限速100km/h条件下的运行时间示意图;Figure 11 is a schematic diagram of the running time under the condition of a station distance of 2800m and a speed limit of 100km/h;
图12站距700m条件下的运行曲线示意图;Figure 12 Schematic diagram of the operation curve under the condition of a station distance of 700m;
图13站距700m且限速100km/h的运行时间示意图;Figure 13 Schematic diagram of the running time with a station distance of 700m and a speed limit of 100km/h;
图14单条曲线定位示意图。Figure 14 Schematic diagram of positioning a single curve.
为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图对本申请实施例中的技术方案进行清楚、完整地描述。In order to make the purpose, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application.
本发明最关键的构思在于:根据专家经验和资料设置参数区间,然后通过本发明设置的列车运行方法生成虚拟数据,再根据研究者要求对数据集进行筛选,得到少量曲线。对得到的少量曲线进行等时跟踪,导出该曲线的运行状态数据集。The most critical idea of the present invention is to set the parameter range according to expert experience and data, then generate virtual data through the train operation method set in the present invention, and then filter the data set according to the researcher's requirements to obtain a small number of curves. Perform isochronous tracking on a small number of curves obtained, and export the running state data set of the curves.
如图1至图14之一所示,本发明公开了基于区间限速和专家经验的虚拟地铁列车运行大数据生成方法,其包括以下步骤:As shown in one of Figures 1 to 14, the present invention discloses a virtual subway train operation big data generation method based on section speed limit and expert experience, which includes the following steps:
步骤1,根据专家经验和资料设置虚拟地铁列车的参数区间,参数区间包括加速度区间、速度区间、车站间距区间和变速距离区间;进一步地,步骤1的参数区间的设置具体包括以下步骤:
步骤1-1,基于实际地铁列车资料计算得到实际地铁列车加速度区间,配置虚拟地铁列车的加速度区间包含实际地铁列车加速度区间;即实际地铁列车加速度区间为虚拟地铁列车的加速度区间的子集;Step 1-1, calculate and obtain the actual subway train acceleration interval based on the actual subway train data, configure the acceleration interval of the virtual subway train to include the actual subway train acceleration interval; that is, the actual subway train acceleration interval is a subset of the acceleration interval of the virtual subway train;
具体地,为了使研究产生的数据集包括各种牵引制动情况,本发明查询了福州地铁运行资料,得出以下数据:福州地铁运行速度在0-40km/h时,平均加速度大于0.83m/s2,运行速度在0-80km/h时,平均加速度大于0.5m/s2。因此地铁列车运行速度曲线数据集生成算法的加速度区间应大于实际地铁列车加速度区间,即大于最大加速度,小于最小加速度,具体数据设置为0.3-1.3m/s2,步长设置为0.1m/s2。Specifically, in order to make the data set produced by the research include various traction and braking conditions, the present invention has inquired about the operation data of Fuzhou subway and obtained the following data: when the running speed of Fuzhou subway is 0-40km/h, the average acceleration is greater than 0.83m/h s2, when the running speed is 0-80km/h, the average acceleration is greater than 0.5m/s2. Therefore, the acceleration interval of the subway train speed curve data set generation algorithm should be greater than the actual subway train acceleration interval, that is, greater than the maximum acceleration and less than the minimum acceleration. The specific data is set to 0.3-1.3m/s2, and the step size is set to 0.1m/s2.
步骤1-2,基于实际地铁列车资料计算得到实际地铁列车运行的速度区间,配置虚拟地 铁列车的速度区间包含实际地铁列车运行的速度区间;即实际地铁列车运行的速度区间为虚拟地铁列车运行的速度区间的子集;Step 1-2, calculate the speed interval of the actual subway train operation based on the actual subway train data, configure the speed interval of the virtual subway train to include the speed interval of the actual subway train operation; that is, the speed interval of the actual subway train operation is the speed interval of the virtual subway train operation a subset of velocity intervals;
具体地,经过资料查阅,得到以下信息:截至2020年12月,福州地铁1号线线路全长29.582km,全部为地下线;共设置25个车站,全部为地下车站;采用6节编组B型地铁列车,最高运行速度为80km/h[13],福州地铁2号线初期配置泉州中车唐车公司31列地铁车辆,地铁列车为4动2拖6辆编组,最大载客量1 880人,最高运行速度80km/h[14]。为使虚拟数据集概括所有的运行情况,设置虚拟数据的速度区间大于实际运行的速度区间,即虚拟最快运行速度应大于实际最快运行速度为100km/h。为考虑乘客的乘坐体验,最慢运行速度不可过低,在此作为本发明的一种实施方式可设置为20km/h。从而得到运行速度区间设置为20-100km/h,步长设置为5km/h。Specifically, after data review, the following information was obtained: As of December 2020, the total length of
步骤1-3,基于实际地铁列车线路的相邻两个车站间距区间,配置虚拟地铁列车线路的相邻两个车站间距区间包含实际地铁列车线路的相邻两个车站间距区间;即实际地铁列车线路的相邻两个车站间距区间为虚拟地铁列车线路的相邻两个车站间距区间的子集;Step 1-3, based on the interval between two adjacent stations of the actual subway train line, configure the interval between two adjacent stations of the virtual subway train line to include the interval between two adjacent stations of the actual subway train line; that is, the actual subway train The interval between two adjacent stations of the line is a subset of the interval between two adjacent stations of the virtual subway train line;
具体地,根据资料显示,调整后福州地铁1号线工程(一期)起点站为象峰站,终点站为福州火车南站,正线线路长约24.618km,共设21个车站,平均站间距为1.202km;截至2019年4月,福州地铁2号线共设置22个车站,全部为地下车站,最大站间距为2.827km(厚庭站至桔园洲站),最小站间距为0.745km(桔园洲站至洪湾站),平均站间距为1.392km。因此,设置车站间距为700-3 000m,即小于最小站间距大于最大站间距以包括所有车站间距的可能性,步长设置为100m;因为地铁列车速度生成算法包括四段变速过程和一段匀速过程,所以将变速距离设置为100-600m,步长设置为50m。Specifically, according to the data, after the adjustment, the starting point of
步骤1-4,将虚拟地铁列车的运行过程分为四段变速过程和一段匀速过程,设定变速距离为100-600m,步长为50m,以形成变速距离区间。In steps 1-4, the running process of the virtual subway train is divided into four stages of speed change and one stage of constant speed. The speed change distance is set to 100-600m, and the step length is 50m to form a speed change distance interval.
步骤2,利用五段式运行法将虚拟地铁列车的运行曲线进行分段,并分别设定每个阶段运行参数;
具体地,在理想运行状态下,把地铁列车看成一个刚体,假设其无摩擦力,无空气阻力,无车厢作用力,提出五段式地铁列车运行法,具体实现方式如图2所示,图中粗体的黑色虚线为地铁列车限速曲线;竖直虚线为运行阶段划分线;灰色点为运行阶段转折点。Specifically, in the ideal operating state, the subway train is regarded as a rigid body, assuming that it has no friction, no air resistance, and no force on the carriage, and a five-stage subway train operation method is proposed. The specific implementation method is shown in Figure 2. The bold black dotted line in the figure is the speed limit curve of the subway train; the vertical dotted line is the division line of the operation phase; the gray point is the turning point of the operation phase.
首先把曲线分为5段(图1中以竖直的虚线进行划分)。在第一次变速阶段,地铁列车以预设第一次变速加速度amax1加速到预设地铁列车首次加速限制速度vmax1,然后匀速行驶完第一次变速剩余预设加速距离Lef1;第二次变速阶段,地铁列车以预设第二次加速度amax2加速到预设地铁列车第二次加速限制速度vmax2,然后进行匀速运动,行驶完第二次变速剩 余预设加速距离Lef2;匀速行驶阶段,地铁列车保持预设地铁列车第二次加速限制速度vmax2匀速行驶完预设匀速行驶距离runtance;后两段为减速行驶,在第三次变速阶段地铁列车先行驶完第三次变速剩余预设变速距离Lef3,再以预设第三次变速加速度amax3进行减速运动,使地铁列车速度减到预设地铁列车首次减速限制速度vmax3;第四次变速阶段与第三次变速行驶相同,先匀速行驶完第四次变速剩余预设变速距离Lef4再以预设第四次变速加速度amax4减速行驶,使地铁列车速度减小到0,即地铁列车到站。First, the curve is divided into 5 segments (divided by vertical dotted lines in Figure 1). In the first shifting stage, the subway train accelerates to the preset first acceleration limit speed vmax1 with the preset first shifting acceleration amax1, and then travels at a constant speed to complete the first shifting remaining preset acceleration distance Lef1; the second shifting stage , the subway train accelerates with the preset second acceleration amax2 to the preset second acceleration limit speed vmax2 of the subway train, and then moves at a constant speed, and the remaining preset acceleration distance Lef2 is left after the second speed change; The preset second acceleration limit speed vmax2 of the preset subway train runs at a constant speed and finishes the preset running distance runtance at a constant speed; the last two sections are decelerated running, and in the third shifting stage, the subway train first runs through the third shifting and the remaining preset shifting distance Lef3, Then decelerate with the preset third speed change acceleration amax3, so that the speed of the subway train is reduced to the preset first deceleration limit speed vmax3 of the subway train; The remaining preset shifting distance Lef4 of shifting is then decelerated at the preset fourth shifting acceleration amax4, so that the speed of the subway train is reduced to 0, that is, the subway train arrives at the station.
进一步地,步骤2中五段式运行法将虚拟地铁列车的运行曲线分为第一次变速阶段、第二次变速阶段、匀速行驶阶段、第三次变速阶段、第四次变速阶段,每个阶段运行参数如下:Further, the five-stage operation method in step 2 divides the operation curve of the virtual subway train into the first shifting stage, the second shifting stage, the uniform speed running stage, the third shifting stage, and the fourth shifting stage, each The stage operation parameters are as follows:
第一次变速阶段:First shift stage:
首先,第一次变速阶段的加速时间t 1由第一次加速限速v max1和第一次变速加速度a max1计算得到,如式(1)所示: First, the acceleration time t1 of the first shifting phase is calculated from the first acceleration limit speed v max1 and the first shifting acceleration a max1 , as shown in formula (1):
第一次变速阶段的匀速行驶时间th 1由第一次变速剩余预设变速距离Lef 1和第一次加速限速v max1计算得到,如式(2)所示: The constant-speed travel time th1 of the first shift stage is calculated from the remaining preset shift distance Lef 1 of the first shift and the first acceleration speed limit v max1 , as shown in formula (2):
其中第一次变速剩余预设变速距离Lef 1由预设变速距离dis和第一次加速距离s 1计算得到,如式(3)所示: The remaining preset shift distance Lef 1 of the first shift is calculated from the preset shift distance dis and the first acceleration distance s 1 , as shown in formula (3):
Lef 1=dis-s 1 (3) Lef 1 = dis-s 1 (3)
第一次加速距离s 1由第一次变速阶段的加速时间t 1和第一次加速限速v max1计算得到,如式(4)所示: The first acceleration distance s 1 is calculated from the acceleration time t 1 of the first shift stage and the first acceleration speed limit v max1 , as shown in formula (4):
则第一阶段的总运行时间T 1由加速时间t 1和匀速行驶时间th 1计算得到,如式(5)所示: Then the total running time T 1 of the first stage is calculated from the acceleration time t 1 and the constant speed travel time th 1 , as shown in formula (5):
T 1=t 1+th 1 (5) T 1 =t 1 +th 1 (5)
第二次变速阶段:The second shift stage:
第二次变速阶段运行时间的计算与第一阶段原理相同,即可通过式(6)~式(10)计算得出。其中涉及的计算变量分别为第一次变速阶段的加速时间t 2、第一次加速限速v max1、第二次加速限速v max2、第二次变速加速度a max2、第二次变速阶段的匀速行驶时间th 2、第二次变速剩余预设变速距离Lef 2、预设变速距离dis、第二次加速距离s 2。 The calculation of the running time of the second shifting stage is the same as the first stage, which can be calculated by formula (6) ~ formula (10). The calculation variables involved are the acceleration time t 2 of the first shifting stage, the first acceleration speed limit v max1 , the second acceleration speed limit v max2 , the second shifting acceleration a max2 , the second shifting speed The constant speed travel time th 2 , the remaining preset shift distance Lef 2 for the second shift, the preset shift distance dis, and the second acceleration distance s 2 .
Lef 2=dis-s 2 (8) Lef 2 = dis-s 2 (8)
T 2=t 2+th 2 (10) T 2 =t 2 +th 2 (10)
匀速行驶阶段:Uniform speed driving stage:
匀速行驶阶段地铁列车的运行时间runtime由地铁列车第二次加速限速v max2和匀速行驶距离runtance计算得到,如式(11)所示: The running time runtime of the subway train at the stage of constant speed running is calculated by the second acceleration speed limit vmax2 of the subway train and the runtance of the constant speed running distance, as shown in formula (11):
第三次变速阶段:The third shift stage:
第三次变速阶段运行时间可通过式(12)~式(16)计算得出。其中涉及的计算变量分别为第三次变速阶段的减速时间t 3、第三次减速限速v max3、第二次加速限速v max2、第三次变速加速度a max3、第三次变速阶段的匀速行驶时间th 3、第三次变速剩余预设变速距离Lef 3、预设变速距离dis、第三次加速距离s 3。 The running time of the third shifting stage can be calculated by formula (12) ~ formula (16). The calculation variables involved are the deceleration time t 3 of the third shifting stage, the third deceleration speed limit v max3 , the second acceleration speed limit v max2 , the third shifting acceleration a max3 , the third shifting speed limit The constant speed travel time th 3 , the remaining preset shift distance Lef 3 for the third shift, the preset shift distance dis, and the third acceleration distance s 3 .
Lef 3=dis-s 3 (14) Lef 3 = dis-s 3 (14)
T 3=t 3+th 3 (16) T 3 =t 3 +th 3 (16)
第四次变速阶段fourth shift stage
第四次变速阶段运行时间可通过式(17)~式(21)计算得出。其中涉及的计算变量分别为第四次变速阶段的减速时间t 4、第三次减速限速v max3、第四次变速加速度a max4、第四次变速阶段的匀速行驶时间th 4、第四次变速剩余预设变速距离Lef 4、预设变速距离dis、第四次加速距离s 4。 The running time of the fourth shifting stage can be calculated by formula (17) ~ formula (21). The calculation variables involved are the deceleration time t 4 of the fourth shift stage, the third deceleration speed limit v max3 , the fourth shift acceleration a max4 , the constant speed travel time th 4 of the fourth shift stage, and the fourth The remaining preset shifting distance Lef 4 , the preset shifting distance dis , and the fourth acceleration distance s 4 .
Lef 4=dis-s 4 (19) Lef 4 = dis-s 4 (19)
T 4=t 4+th 4 (21) T 4 =t 4 +th 4 (21)
最快运行状态运行时间计算:Calculation of running time in the fastest running state:
最快运行状态下的运行时间alltime由四次变速阶段运行时间T 1、T 2、T 3、T 4和匀速运行时间runtime计算得到,如式(22)。 The running time alltime in the fastest running state is calculated from the running time T 1 , T 2 , T 3 , T 4 of the four shifting stages and the running time runtime at a constant speed, as shown in formula (22).
alltime=T 1+T 2+T 3+T 4+runtime (22)。 alltime=T 1 +T 2 +T 3 +T 4 +runtime (22).
步骤3,利用设定的五段式运行法的每个阶段运行参数计算得到若干组虚拟地铁列车的运行数据,并分别利用绘图工具生成对应运行曲线;
具体地,以地铁列车的最快运行状态和最慢运行状态为例,分析五段式运行方式,如图3所示:第一种运行状态是紧贴地铁列车限速曲线运行的最快运行状态,在此状态下,把地铁列车首次加速限制速度vmax1、首次减速限制速度vmax3设置为地铁列车最低限速Emin,地铁列车第二次加速限制速度vmax2设置为地铁列车最低限速Emax;变速阶段的加速度amax1、amax2、amax3、amax4均为最快加速度;第二种运行状态是把预设地铁列车首次加速限制速度vmax1、首次减速限制速度vmax3设置为预设地铁列车最低限速Emin的一半,预设地铁列车第二次减速限制速度vmax2设置为预设地铁列车最低限速Emin;变速阶段的加速度amax1、amax2、amax3、amax4均为最慢加速度,得到最慢运行状态。这两种运行状态为生成曲线的上边界和下边界,其他的生成曲线都在两者之间,如图4所示。Specifically, taking the fastest running state and the slowest running state of the subway train as an example, the five-stage running mode is analyzed, as shown in Figure 3: the first running state is the fastest running close to the speed limit curve of the subway train state, in this state, the first acceleration limit speed vmax1 and the first deceleration limit speed vmax3 of the subway train are set as the minimum speed limit Emin of the subway train, and the second acceleration limit speed vmax2 of the subway train is set as the minimum speed limit Emax of the subway train; the speed change stage The accelerations amax1, amax2, amax3, and amax4 are all the fastest accelerations; the second running state is to set the first acceleration limit speed vmax1 of the preset subway train and the first deceleration limit speed vmax3 to half of the preset minimum speed limit Emin of the subway train. The second deceleration limit speed vmax2 of the preset subway train is set to the minimum speed limit Emin of the preset subway train; the accelerations amax1, amax2, amax3, and amax4 of the speed change stage are all the slowest accelerations, and the slowest running state is obtained. These two operating states are the upper boundary and lower boundary of the generated curve, and other generated curves are in between, as shown in Figure 4.
根据上述假设进行数据分级,可得到地铁列车运行数据集,共312 800 0条数据,可通过Python编程实现并导出为CSV文件,其中列变量包括最低限速、最高限速、加速度、预设加速距离、两站距离、四段加速距离、运行时间等参数。在生成的数据集中,剔除不符合实际情况(如运行时间大于600s、不符合预设运行模型)的曲线后,得到209 005 0条数据,这些数据的运行时间频率分布情况如图5所示。According to the above assumptions, the data is classified, and the subway train operation data set can be obtained, with a total of 312 800 0 data, which can be realized by Python programming and exported as a CSV file. The column variables include the minimum speed limit, maximum speed limit, acceleration, preset acceleration Distance, two-station distance, four-segment acceleration distance, running time and other parameters. In the generated data set, 209 005 0 pieces of data were obtained after eliminating the curves that did not conform to the actual situation (such as running time greater than 600s, not conforming to the preset running model). The frequency distribution of running time of these data is shown in Figure 5.
步骤4,根据研究需求设定筛选条件得到对应的运行曲线;
具体地,为验证数据集的可用性,本发明分别以福州地铁最快运行速度80km/h、数据集设置最快运行速度100km/h、福州地铁2号线最大站间距2800m、数据集两站最短间距700m为例进行筛选,可得到以下结果。Specifically, in order to verify the usability of the data set, the present invention set the fastest running speed of Fuzhou Metro to 80km/h, the fastest running speed of the data set to 100km/h, the maximum station distance of
(4.1)以运行速度作为筛选条件:(4.1) Use the running speed as the filter condition:
设置最高限速为福州地铁实际运行速度的最高限速80km/h,其余参数设置为最低限速30km/h、两站间距1600m、加速距离300m,得到2450条数据,经过matplotlib绘图可得到如6所示的运行曲线,并得到其如图7所示的运行时间。Set the maximum speed limit to the actual operating speed of Fuzhou Metro at 80km/h, and set the other parameters to the minimum speed limit of 30km/h, the distance between two stations to 1600m, and the acceleration distance to 300m, and obtain 2450 pieces of data, which can be obtained through matplotlib drawing such as 6 As shown in the running curve, and get its running time as shown in Figure 7.
设置最高限速100km/h,最低限速50km/h、两站间距1600m、加速距离300m,得到5510条数据,经过matplotlib绘图可得到如8所示的运行曲线,并得到其如图9所示的运行时间。Set the maximum speed limit to 100km/h, the minimum speed limit to 50km/h, the distance between two stations to 1600m, and the acceleration distance to 300m, and obtain 5510 pieces of data. After matplotlib drawing, the running curve shown in Figure 8 can be obtained, and it can be obtained as shown in Figure 9 running time.
(4.2)以两站间距作为筛选条件(4.2) Take the distance between two stations as the filter condition
设置两站间距为福州地铁2号线最大站间距2800m,其余参数设置为最低限速30km/h、最高限速50km/h、加速距离500m,得到350条数据,经过matplotlib绘图可得到如10所示的运行曲线,并得到其如图11所示的运行时间。Set the distance between two stations as the maximum station distance of
设置两站间距为700m,其余参数设置为最低限速30km/h、最高限速50km/h、加速距离100m,得到350条数据,经过matplotlib绘图可得到如图12所示的运行曲线,并得到其如图13所示的运行时间。Set the distance between the two stations to 700m, and set the other parameters to the minimum speed limit of 30km/h, the maximum speed limit of 50km/h, and the acceleration distance of 100m, and obtain 350 pieces of data. After drawing with matplotlib, the running curve shown in Figure 12 can be obtained, and Its running time is shown in Figure 13.
步骤5,对每条筛选出的曲线进行等时间隔的曲线数据采样,并输出每个时间点列车所处状态。Step 5: Sampling curve data at equal time intervals for each selected curve, and outputting the state of the train at each time point.
具体地,如图14所示,对于筛选出的曲线,对每条曲线进行0.1s为间隔的数据采样,输出每个时间点列车所处状态,Specifically, as shown in Figure 14, for the screened curves, data sampling is performed at an interval of 0.1s for each curve, and the state of the train at each time point is output,
进一步地,步骤5中每个时间点列车所处状态包括当前位置加速度m/s 2,当前位置已运行时间s、当前位置运行速度m/s、当前位置运行距离m、当前位置限速m/s、保持当前行驶状态时间s、当前位置剩余变速距离m、当前位置剩余变速范围m/s、当前位置到站剩余时间s和当前位置到站剩余距离m。 Further, the state of the train at each time point in step 5 includes current position acceleration m/s 2 , current position elapsed running time s, current position running speed m/s, current position running distance m, current position speed limit m/s s, the time to maintain the current driving state s, the remaining shifting distance at the current position m, the remaining shifting range at the current position m/s, the remaining time at the current position at the station s, and the remaining distance at the current position at the station m.
本发明突破性的用人工生成大量地铁运行曲线的方式代替传统获得地铁列车运行数据的方法。具体通过调查地铁运行的实际数据和咨询有关专家来获得虚拟数据的牵引区间、速度区间、站距区间、变速区间等限定参数,使得生成的虚拟数据范围小,信度高;然后根据动力学知识假设地铁列车的五段式运行方式,以便于生成列车的虚拟速度-距离图像;再而采用python进行编程实现和数据可视化,得到本发明所需要的大量虚拟地铁列车运行数据。对于生成的大量虚拟数据,本发明对其进行条件筛选,得到研究者所需要的部分虚线。另外,对于选出的单条曲线,本发明对其进行了实时跟踪,得到了单条曲线地铁列车每0.1s的运行状态,具体包括:当前位置加速度、当前位置已运行时间、当前位置运行速度、当前位置运行距离、当前位置限速、保持当前行驶状态时间、当前位置剩余变速距离、当前位置剩余变速范围、当前位置到站剩余时间、当前位置到站剩余距离,根据这些数据本发明可以了解到地铁列车的实时运行状态,以便于研究者对地铁列车的进一步研究。The invention makes a breakthrough in replacing the traditional method of obtaining subway train operation data by artificially generating a large number of subway operation curves. Specifically, by investigating the actual data of subway operation and consulting relevant experts to obtain the limited parameters of the virtual data such as the traction interval, speed interval, station distance interval, and speed change interval, so that the generated virtual data has a small range and high reliability; and then according to the dynamics knowledge Assume the five-stage running mode of the subway train, so as to generate the virtual speed-distance image of the train; then use python to carry out programming and data visualization to obtain a large amount of virtual subway train running data required by the present invention. For the generated large amount of virtual data, the present invention screens them according to the conditions, and obtains the partial dotted lines required by the researcher. In addition, for the selected single curve, the present invention tracks it in real time, and obtains the running state of the single curve subway train every 0.1s, specifically including: current position acceleration, current position elapsed running time, current position running speed, current The running distance at the current location, the speed limit at the current location, the time to maintain the current driving state, the remaining shifting distance at the current location, the remaining shifting range at the current location, the remaining time at the current location to the station, and the remaining distance at the current location to the station. According to these data, the present invention can understand the The real-time running status of the train is convenient for researchers to further study the subway train.
显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。通常在此处附图中描述和示出的本申请实施例的组件可以以各种不同的配置来布置和设计。因此,本申请的实施例的详细描述并非旨在限制要求保护的本申请的范围,而是仅仅表示本申请的选定实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。Apparently, the described embodiments are some of the embodiments of the present application, but not all of them. In the case of no conflict, the embodiments in the present application and the features in the embodiments can be combined with each other. The components of the embodiments of the application generally described and illustrated in the figures herein may be arranged and designed in a variety of different configurations. Therefore, the detailed description of the embodiments of the present application is not intended to limit the scope of the claimed application, but merely represents selected embodiments of the application. Based on the embodiments in this application, all other embodiments obtained by persons of ordinary skill in the art without creative efforts fall within the protection scope of this application.
Claims (7)
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210044225.1 | 2022-01-14 | ||
CN202210044225.1A CN114348070A (en) | 2022-01-14 | 2022-01-14 | Virtual subway train operation big data generation method based on interval speed limit and expert experience |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2023134093A1 true WO2023134093A1 (en) | 2023-07-20 |
Family
ID=81091314
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2022/092115 WO2023134093A1 (en) | 2022-01-14 | 2022-05-11 | Virtual subway train operation big data generation method based on interval speed limit and expert experience |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN114348070A (en) |
WO (1) | WO2023134093A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN119247745A (en) * | 2024-12-03 | 2025-01-03 | 中国地质调查局水文地质环境地质调查中心 | Automatic control system of variable frequency fixed depth winch for carbon dioxide leakage sampling |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114348070A (en) * | 2022-01-14 | 2022-04-15 | 福建工程学院 | Virtual subway train operation big data generation method based on interval speed limit and expert experience |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012192825A (en) * | 2011-03-16 | 2012-10-11 | Railway Technical Research Institute | Train run curve editing method and train run curve editing system |
CN104134378A (en) * | 2014-06-23 | 2014-11-05 | 北京交通大学 | Urban rail train intelligent control method based on driving experience and online study |
CN107878510A (en) * | 2016-12-29 | 2018-04-06 | 比亚迪股份有限公司 | Automatic train control method and device, vehicle-mounted ATO |
CN114348070A (en) * | 2022-01-14 | 2022-04-15 | 福建工程学院 | Virtual subway train operation big data generation method based on interval speed limit and expert experience |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111797473B (en) * | 2020-06-28 | 2024-11-15 | 通号城市轨道交通技术有限公司 | A subway train mainline operation simulation calculation method and device |
CN113283065A (en) * | 2021-05-10 | 2021-08-20 | 中铁第一勘察设计院集团有限公司 | Subway section air shaft setting method based on train operation simulation |
-
2022
- 2022-01-14 CN CN202210044225.1A patent/CN114348070A/en active Pending
- 2022-05-11 WO PCT/CN2022/092115 patent/WO2023134093A1/en active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012192825A (en) * | 2011-03-16 | 2012-10-11 | Railway Technical Research Institute | Train run curve editing method and train run curve editing system |
CN104134378A (en) * | 2014-06-23 | 2014-11-05 | 北京交通大学 | Urban rail train intelligent control method based on driving experience and online study |
CN107878510A (en) * | 2016-12-29 | 2018-04-06 | 比亚迪股份有限公司 | Automatic train control method and device, vehicle-mounted ATO |
CN114348070A (en) * | 2022-01-14 | 2022-04-15 | 福建工程学院 | Virtual subway train operation big data generation method based on interval speed limit and expert experience |
Non-Patent Citations (1)
Title |
---|
LU YUQI U, DEWANG CHEN, ZHAOLIN ZHAO: "Algorithm for automatically generating a large number of speed curves of subway trains based on AlphaZero ", CHINESE JOURNAL OF INTELLIGENT SCIENCE AND TECHNOLOGY, vol. 3, no. 2, 1 January 2021 (2021-01-01), pages 179 - 184, XP093079026, DOI: 10.11959/j.issn.2096-6652.202118 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN119247745A (en) * | 2024-12-03 | 2025-01-03 | 中国地质调查局水文地质环境地质调查中心 | Automatic control system of variable frequency fixed depth winch for carbon dioxide leakage sampling |
Also Published As
Publication number | Publication date |
---|---|
CN114348070A (en) | 2022-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107117170B (en) | A real-time predictive cruise control system based on economical driving | |
CN105243430B (en) | The optimization method of the target velocity curve of energy-saving train operation | |
CN111216713B (en) | A kind of automatic driving vehicle speed preview control method | |
CN102981408B (en) | Running process modeling and adaptive control method for motor train unit | |
US8660723B2 (en) | Method for determining run-curves for vehicles in real-time subject to dynamic travel time and speed limit constraint | |
CN111267830B (en) | A hybrid electric bus energy management method, device and storage medium | |
US9002547B2 (en) | System and method for determining dynamically changing distributions of vehicles in a vehicle system | |
WO2023134093A1 (en) | Virtual subway train operation big data generation method based on interval speed limit and expert experience | |
US9669811B2 (en) | System and method for asynchronously controlling brakes of vehicles in a vehicle system | |
US9453735B2 (en) | System and method for determining operational group assignments of vehicles in a vehicle system | |
CN107264534A (en) | Intelligent driving control system and method, vehicle based on driver experience's model | |
CN111091721A (en) | Ramp confluence control method and system for intelligent train traffic system | |
CN108883784A (en) | The method and train driver consulting system of drive advice are provided to train driver | |
CN103955135B (en) | A Calculation Method of Locomotive Optimal Maneuvering Sequence Based on Double-layer Model Curve | |
Zhang et al. | Data-driven train operation models based on data mining and driving experience for the diesel-electric locomotive | |
CN110555476A (en) | intelligent vehicle track change track prediction method suitable for man-machine hybrid driving environment | |
CN112820126B (en) | A right-of-way priority operation control and simulation method for a non-intrusive guided transport vehicle | |
CN108829110A (en) | A kind of pilot model modeling method of cross/longitudinal movement Unified frame | |
Bin et al. | Energy-efficient speed profile optimization for high-speed railway considering neutral sections | |
CN110703757A (en) | Energy consumption optimization-oriented high-speed train speed planning method | |
Ying et al. | A sliding mode control approach to longitudinal control of vehicles in a platoon | |
CN106647269B (en) | A kind of locomotive smart steering optimized calculation method | |
Hamid et al. | Investigation into train positioning systems for saving energy with optimised train trajectories | |
CN115320616A (en) | Control method, device, equipment and medium for automatically driving vehicle speed | |
CN106585675B (en) | Train operation optimized handling method and apparatus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 22919733 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 22919733 Country of ref document: EP Kind code of ref document: A1 |