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 PDF

Info

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
Application number
PCT/CN2022/092115
Other languages
French (fr)
Chinese (zh)
Inventor
陈德旺
芦玉琦
Original Assignee
福建工程学院
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 福建工程学院 filed Critical 福建工程学院
Publication of WO2023134093A1 publication Critical patent/WO2023134093A1/en

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B61RAILWAYS
    • B61LGUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
    • B61L27/00Central railway traffic control systems; Trackside control; Communication systems specially adapted therefor
    • B61L27/40Handling 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

Disclosed in the present invention is a virtual subway train operation big data generation method based on an interval speed limit and expert experience. Limited parameters of virtual data, such as a traction interval, a speed interval, a station spacing interval and a speed change interval, are obtained by investigating actual data of subway operation and consulting related experts, such that the generated virtual data is in a small range and has high reliability. Then, a five-section operation mode of a subway train is assumed according to dynamic knowledge, so as to generate a virtual speed-distance image of the train. Afterwards, programming implementation and data visualization are performed by using Python, so as to obtain a large amount of required virtual subway train operation data. Conditional screening is performed on a large amount of generated virtual data, so as to obtain some dotted lines required by a researcher. In addition, a selected single curve is tracked in real time, so as to obtain an operation state of the subway train on the single curve every 0.1 second, such that the real-time operation state of the subway train is learned, thereby facilitating the researcher in further researching the subway train.

Description

基于区间限速和专家经验的虚拟地铁列车运行大数据生成方法Big data generation method for virtual subway train operation based on section speed limit and expert experience 技术领域technical field

本发明涉及数据处理技术领域,尤其涉及基于区间限速和专家经验的虚拟地铁列车运行大数据生成方法。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.

背景技术Background technique

近年来,地铁列车的运行数据获取方式可分为跟踪仿真生成的单条运行曲线和人工驾驶数据的采样两种。其中,仿真生成单条运行曲线有以下研究:通过读入数据、数据预处理、设定计算区间、计算并输出结果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,根据专家经验和资料设置虚拟地铁列车的参数区间,参数区间包括加速度区间、速度区间、车站间距区间和变速距离区间;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;

步骤2,利用五段式运行法将虚拟地铁列车的运行曲线进行分段,并分别设定每个阶段运行参数;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;

步骤3,利用设定的五段式运行法的每个阶段运行参数计算得到若干组虚拟地铁列车的运行数据,并分别利用绘图工具生成对应运行曲线;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;

步骤4,根据研究需求设定筛选条件得到对应的运行曲线;Step 4, set screening conditions according to research requirements to obtain corresponding operating curves;

步骤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 step 1 specifically includes the following steps:

步骤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 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:

第一次变速阶段: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):

Figure PCTCN2022092115-appb-000001
Figure PCTCN2022092115-appb-000001

第一次变速阶段的匀速行驶时间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):

Figure PCTCN2022092115-appb-000002
Figure PCTCN2022092115-appb-000002

其中第一次变速剩余预设变速距离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):

Figure PCTCN2022092115-appb-000003
Figure PCTCN2022092115-appb-000003

则第一阶段的总运行时间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 2The 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 .

Figure PCTCN2022092115-appb-000004
Figure PCTCN2022092115-appb-000004

Figure PCTCN2022092115-appb-000005
Figure PCTCN2022092115-appb-000005

Lef 2=dis-s 2        (8) Lef 2 = dis-s 2 (8)

Figure PCTCN2022092115-appb-000006
Figure PCTCN2022092115-appb-000006

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):

Figure PCTCN2022092115-appb-000007
Figure PCTCN2022092115-appb-000007

第三次变速阶段:The third shift stage:

第三次变速阶段运行时间可通过式(12)~式(16)计算得出。其中涉及的计算变量分别为第三次变速阶段的减速时间t 3、第三次减速限速v max3、第二次加速限速v max2、第三次变速加速度a max3、第三次变速阶段的匀速行驶时间th 3、第三次变速剩余预设变速距离Lef 3、预设变速距离dis、第三次加速距离s 3The 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 .

Figure PCTCN2022092115-appb-000008
Figure PCTCN2022092115-appb-000008

Figure PCTCN2022092115-appb-000009
Figure PCTCN2022092115-appb-000009

Lef 3=dis-s 3         (14) Lef 3 = dis-s 3 (14)

Figure PCTCN2022092115-appb-000010
Figure PCTCN2022092115-appb-000010

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 4The 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 .

Figure PCTCN2022092115-appb-000011
Figure PCTCN2022092115-appb-000011

Figure PCTCN2022092115-appb-000012
Figure PCTCN2022092115-appb-000012

Lef 4=dis-s 4        (19) Lef 4 = dis-s 4 (19)

Figure PCTCN2022092115-appb-000013
Figure PCTCN2022092115-appb-000013

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 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.

进一步地,步骤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.

附图说明Description of drawings

以下结合附图和具体实施方式对本发明做进一步详细说明;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.

具体实施方式Detailed ways

为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图对本申请实施例中的技术方案进行清楚、完整地描述。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的参数区间的设置具体包括以下步骤: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:

步骤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 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]. In order to make the virtual data set summarize all the running conditions, 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. In order to consider the riding experience of the passengers, the slowest running speed can not be too low, which can be set to 20km/h as an embodiment of the present invention. Thus, the operating speed range is set to 20-100km/h, and the step size is set to 5km/h.

步骤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 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. Therefore, 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.

步骤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,利用五段式运行法将虚拟地铁列车的运行曲线进行分段,并分别设定每个阶段运行参数;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;

具体地,在理想运行状态下,把地铁列车看成一个刚体,假设其无摩擦力,无空气阻力,无车厢作用力,提出五段式地铁列车运行法,具体实现方式如图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):

Figure PCTCN2022092115-appb-000014
Figure PCTCN2022092115-appb-000014

第一次变速阶段的匀速行驶时间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):

Figure PCTCN2022092115-appb-000015
Figure PCTCN2022092115-appb-000015

其中第一次变速剩余预设变速距离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):

Figure PCTCN2022092115-appb-000016
Figure PCTCN2022092115-appb-000016

则第一阶段的总运行时间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 2The 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 .

Figure PCTCN2022092115-appb-000017
Figure PCTCN2022092115-appb-000017

Figure PCTCN2022092115-appb-000018
Figure PCTCN2022092115-appb-000018

Lef 2=dis-s 2       (8) Lef 2 = dis-s 2 (8)

Figure PCTCN2022092115-appb-000019
Figure PCTCN2022092115-appb-000019

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):

Figure PCTCN2022092115-appb-000020
Figure PCTCN2022092115-appb-000020

第三次变速阶段:The third shift stage:

第三次变速阶段运行时间可通过式(12)~式(16)计算得出。其中涉及的计算变量分别为第三次变速阶段的减速时间t 3、第三次减速限速v max3、第二次加速限速v max2、第三次变速加速度a max3、第三次变速阶段的匀速行驶时间th 3、第三次变速剩余预设变速距离Lef 3、预设变速距离dis、第三次加速距离s 3The 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 .

Figure PCTCN2022092115-appb-000021
Figure PCTCN2022092115-appb-000021

Figure PCTCN2022092115-appb-000022
Figure PCTCN2022092115-appb-000022

Lef 3=dis-s 3          (14) Lef 3 = dis-s 3 (14)

Figure PCTCN2022092115-appb-000023
Figure PCTCN2022092115-appb-000023

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 4The 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 .

Figure PCTCN2022092115-appb-000024
Figure PCTCN2022092115-appb-000024

Figure PCTCN2022092115-appb-000025
Figure PCTCN2022092115-appb-000025

Lef 4=dis-s 4          (19) Lef 4 = dis-s 4 (19)

Figure PCTCN2022092115-appb-000026
Figure PCTCN2022092115-appb-000026

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,利用设定的五段式运行法的每个阶段运行参数计算得到若干组虚拟地铁列车的运行数据,并分别利用绘图工具生成对应运行曲线;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;

具体地,以地铁列车的最快运行状态和最慢运行状态为例,分析五段式运行方式,如图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,根据研究需求设定筛选条件得到对应的运行曲线;Step 4, set screening conditions according to research requirements to obtain corresponding operating curves;

具体地,为验证数据集的可用性,本发明分别以福州地铁最快运行速度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 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.

(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 Fuzhou Metro Line 2 to 2800m, and set the other parameters as the minimum speed limit of 30km/h, the maximum speed limit of 50km/h, and the acceleration distance of 500m to obtain 350 pieces of data, which can be obtained through matplotlib drawing such as 10 The running curve shown, and get its running time as shown in Figure 11.

设置两站间距为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)

基于区间限速和专家经验的虚拟地铁列车运行大数据生成方法,其特征在于:其包括以下步骤:The big data generation method of virtual subway train operation based on interval speed limit and expert experience is characterized in that: it comprises the following steps: 步骤1,根据专家经验和资料设置虚拟地铁列车的参数区间,参数区间包括加速度区间、速度区间、车站间距区间和变速距离区间;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; 步骤2,利用五段式运行法将虚拟地铁列车的运行曲线进行分段,并分别设定每个阶段运行预设参数;Step 2, use the five-stage operation method to segment the operation curve of the virtual subway train, and set the operation preset parameters for each stage respectively; 步骤3,利用设定的五段式运行法的每个阶段运行参数计算得到若干组虚拟地铁列车的运行数据,并分别利用绘图工具生成对应运行曲线;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; 步骤4,根据研究需求设定筛选条件得到对应的运行曲线;Step 4, set screening conditions according to research requirements to obtain corresponding operating curves; 步骤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所述的基于区间限速和专家经验的虚拟地铁列车运行大数据生成方法,其特征在于:步骤1的参数区间的设置具体包括以下步骤:The virtual subway train operation big data generation method based on interval speed limit and expert experience according to claim 1, is characterized in that: the setting of the parameter interval of step 1 specifically comprises the following steps: 步骤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. 根据权利要求1所述的基于区间限速和专家经验的虚拟地铁列车运行大数据生成方法,其特征在于:步骤2中五段式运行法将虚拟地铁列车的运行曲线依次分为第一次变速阶段、第二次变速阶段、匀速行驶阶段、第三次变速阶段和第四次变速阶段,四段变速过程和一段匀速过程的预设参数包括预设变速距离dis、第一次加速限速v max 1、表示第二次变速加速度a max 2和第一次变速加速度a max 1;匀速行驶距离runtance。 According to claim 1, the virtual subway train operation big data generation method based on interval speed limit and expert experience is characterized in that: in step 2, the five-stage operation method divides the operation curve of the virtual subway train into the first speed change successively stage, the second speed change stage, the constant speed driving stage, the third speed change stage and the fourth speed change stage, the preset parameters of the four-stage speed change process and the one-stage constant speed process include the preset speed change distance dis, the first acceleration speed limit v max 1 means the second speed change acceleration a max 2 and the first time speed change acceleration a max 1 ; the running distance runtance at a constant speed. 根据权利要求3所述的基于区间限速和专家经验的虚拟地铁列车运行大数据生成方法,其特征在于:步骤2中每个阶段运行参数的计算如下:According to claim 3, based on section speed limit and expert experience virtual subway train operation big data generation method, it is characterized in that: in step 2, the calculation of each stage operation parameter is as follows: 第一次变速阶段:First shift stage: 第一次变速阶段的加速时间t 1由第一次加速限速v max 1和第一次变速加速度a max 1计算得到,如式(1)所示: The acceleration time t1 of the first shift stage is calculated from the first acceleration limit speed v max 1 and the first shift acceleration a max 1 , as shown in formula (1):
Figure PCTCN2022092115-appb-100001
Figure PCTCN2022092115-appb-100001
第一次变速阶段的匀速行驶时间th 1由第一次变速剩余预设变速距离Lef 1和第一次加速限速v max 1计算得到,如式(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 max 1 , as shown in formula (2):
Figure PCTCN2022092115-appb-100002
Figure PCTCN2022092115-appb-100002
第一次变速剩余预设变速距离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 max 1计算得到,如式(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 max 1 , as shown in formula (4):
Figure PCTCN2022092115-appb-100003
Figure PCTCN2022092115-appb-100003
第一阶段的总运行时间T 1由加速时间t 1和匀速行驶时间th 1计算得到,如式(5)所示: 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: 各个变量的计算公式分别为The calculation formulas for each variable are
Figure PCTCN2022092115-appb-100004
Figure PCTCN2022092115-appb-100004
Figure PCTCN2022092115-appb-100005
Figure PCTCN2022092115-appb-100005
Lef 2=dis-s 2  (8) Lef 2 = dis-s 2 (8)
Figure PCTCN2022092115-appb-100006
Figure PCTCN2022092115-appb-100006
T 2=t 2+th 2  (10) T 2 =t 2 +th 2 (10) 其中,t 2表示第一次变速阶段的加速时间、v max 1表示第一次加速限速、v max 2表示第二次加速限速、a max 2表示第二次变速加速度、th 2表示第二次变速阶段的匀速行驶时间、Lef 2表示第二次变速剩余预设变速距离、dis表示预设变速距离、s 2表示第二次加速距离; Among them, t 2 represents the acceleration time of the first shift stage, v max 1 represents the first acceleration speed limit, v max 2 represents the second acceleration speed limit, a max 2 represents the second speed change acceleration, and th 2 represents the first speed limit The constant-speed driving time in the second shift stage, Lef 2 represents the remaining preset shift distance for the second shift, dis represents the preset shift distance, and s 2 represents the second acceleration distance; 匀速行驶阶段:Uniform speed driving stage: 匀速行驶阶段地铁列车的运行时间runtime由地铁列车第二次加速限速v max 2和匀速行驶距离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 v max 2 of the subway train and the runtance of the constant speed running distance, as shown in formula (11):
Figure PCTCN2022092115-appb-100007
Figure PCTCN2022092115-appb-100007
第三次变速阶段:The third shift stage: 各个计算变量的计算公式分别为The calculation formulas for each calculation variable are
Figure PCTCN2022092115-appb-100008
Figure PCTCN2022092115-appb-100008
Figure PCTCN2022092115-appb-100009
Figure PCTCN2022092115-appb-100009
Lef 3=dis-s 3  (14) Lef 3 = dis-s 3 (14)
Figure PCTCN2022092115-appb-100010
Figure PCTCN2022092115-appb-100010
T 3=t 3+th 3  (16) T 3 =t 3 +th 3 (16) 其中,涉及的计算变量分别为第三次变速阶段的减速时间t 3、第三次减速限速v max 3、第二次加速限速v max 2、第三次变速加速度a max 3、第三次变速阶段的匀速行驶时间th 3、第三次变速剩余预设变速距离Lef 3、预设变速距离dis、第三次加速距离s 3Among them, the calculation variables involved are the deceleration time t 3 of the third shift stage, the third deceleration speed limit v max 3 , the second acceleration speed limit v max 2 , the third speed change acceleration a max 3 , the third speed limit The constant speed travel time th 3 of the secondary shifting stage, the remaining preset shifting distance Lef 3 of the third shifting, the preset shifting distance dis, and the third acceleration distance s 3 ; 第四次变速阶段:The fourth shift stage: 各个计算变量的计算公式分别为The calculation formulas for each calculation variable are
Figure PCTCN2022092115-appb-100011
Figure PCTCN2022092115-appb-100011
Figure PCTCN2022092115-appb-100012
Figure PCTCN2022092115-appb-100012
Lef 4=dis-s 4  (19) Lef 4 = dis-s 4 (19)
Figure PCTCN2022092115-appb-100013
Figure PCTCN2022092115-appb-100013
T 4=t 4+th 4  (21) T 4 =t 4 +th 4 (21) 其中,涉及的计算变量分别为第四次变速阶段的减速时间t 4、第三次减速限速v max 3、第四次变速加速度a max 4、第四次变速阶段的匀速行驶时间th 4、第四次变速剩余预设变速距离Lef 4、预设变速距离dis、第四次加速距离s 4Among them, the calculation variables involved are the deceleration time t 4 of the fourth shift stage, the third deceleration speed limit v max 3 , the fourth shift acceleration a max 4 , the constant speed travel time th 4 of the fourth shift stage, The remaining preset shift distance Lef 4 for the fourth shift, the preset shift distance dis, and the fourth acceleration distance s 4 ; 最快运行状态运行时间计算: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).
根据权利要求1所述的基于区间限速和专家经验的虚拟地铁列车运行大数据生成方法,其特征在于:步骤3中绘制形成运行曲线以最快运行状态对应曲线为上界并以最慢运行状态对应曲线作为下界。According to claim 1, the virtual subway train operation big data generation method based on interval speed limit and expert experience, is characterized in that: in step 3, draw and form the operation curve with the fastest operation state corresponding curve as the upper bound and the slowest operation The state-corresponding curve is used as the lower bound. 根据权利要求1所述的基于区间限速和专家经验的虚拟地铁列车运行大数据生成方法,其特征在于:步骤5中数据采样的间隔为0.1秒。The method for generating big data of virtual subway train operation based on section speed limit and expert experience according to claim 1, characterized in that: the interval of data sampling in step 5 is 0.1 second. 根据权利要求1所述的基于区间限速和专家经验的虚拟地铁列车运行大数据生成方法,其特征在于:步骤5中每个时间点列车所处状态包括当前位置加速度m/s 2,当前位置已运行时间s、当前位置运行速度m/s、当前位置运行距离m、当前位置限速m/s、保持当前行驶状态时间s、当前位置剩余变速距离m、当前位置剩余变速范围m/s、当前位置到站剩余时间s和当前位置到站剩余距离m。 According to claim 1, based on section speed limit and expert experience virtual subway train running big data generation method, it is characterized in that: in step 5, the state of train at each time point 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, current driving state time s, current position remaining shifting distance m, current position remaining shifting range m/s, The remaining time to the station at the current position s and the remaining distance to the station at the current position m.
PCT/CN2022/092115 2022-01-14 2022-05-11 Virtual subway train operation big data generation method based on interval speed limit and expert experience WO2023134093A1 (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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