CN111951572A - Time interval division optimization method for multi-time interval signal control scheme of urban road intersection - Google Patents

Time interval division optimization method for multi-time interval signal control scheme of urban road intersection Download PDF

Info

Publication number
CN111951572A
CN111951572A CN202010645121.7A CN202010645121A CN111951572A CN 111951572 A CN111951572 A CN 111951572A CN 202010645121 A CN202010645121 A CN 202010645121A CN 111951572 A CN111951572 A CN 111951572A
Authority
CN
China
Prior art keywords
time
time interval
phase
rate
period
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010645121.7A
Other languages
Chinese (zh)
Other versions
CN111951572B (en
Inventor
毛秀星
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Traffic Police Brigade Of Yongjia County Public Security Bureau
Original Assignee
Traffic Police Brigade Of Yongjia County Public Security Bureau
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 Traffic Police Brigade Of Yongjia County Public Security Bureau filed Critical Traffic Police Brigade Of Yongjia County Public Security Bureau
Priority to CN202010645121.7A priority Critical patent/CN111951572B/en
Publication of CN111951572A publication Critical patent/CN111951572A/en
Application granted granted Critical
Publication of CN111951572B publication Critical patent/CN111951572B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/07Controlling traffic signals
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0125Traffic data processing
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0125Traffic data processing
    • G08G1/0129Traffic data processing for creating historical data or processing based on historical data
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0137Measuring and analyzing of parameters relative to traffic conditions for specific applications

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Traffic Control Systems (AREA)

Abstract

A time interval division optimization method for a multi-time interval signal control scheme of an urban road intersection starts with queuing and evacuating time from an inlet to a hierarchy, a phase corresponding to minimum time granularity, a flow direction vacancy rate and an unbalance rate are calibrated, then the utilization condition of adjacent time intervals in green time of different flow directions and the requirement unbalance condition of different flow directions in the same phase are analyzed, and the feasibility of combining the adjacent time intervals is researched and judged on the basis. By multiple analysis and combination, the original time interval division scheme is optimized, and the signal optimization target under the unsaturated and saturated supersaturated states of the intersection is met, so that the signal control time interval division with higher reliability and higher stability is realized.

Description

Time interval division optimization method for multi-time interval signal control scheme of urban road intersection
Technical Field
The invention relates to the field of road traffic control, in particular to a time interval division optimization method for a multi-time interval signal control scheme of an urban road intersection.
Background
In order to meet traffic demands in different time periods throughout the day, urban road traffic signal control generally adopts a multi-period signal control scheme, namely, a plurality of time periods are divided, and the same signal is adopted to control time-sharing parameters such as a period, a phase sequence, a green signal ratio and the like in the same time period; for signal single-point optimization, time interval division optimization is a primary link and is a precondition and basis for optimization of subsequent signal timing parameters.
The traditional time interval division method takes traffic police control experience as a leading factor, divides the time intervals of high peak, flat peak and low peak according to the traffic load capacity or queuing characteristics of intersections, and further configures a corresponding fixed signal control scheme for each time interval. Under a data-oriented control mode, an automatic time interval division method based on traffic flow operation detection data mostly extracts traffic operation similarity characteristics of an intersection entrance lane level or a lane level according to basic traffic flow parameters such as historical traffic flow, queuing length and the like and data mining methods such as a threshold value method, a statistical method, a clustering method and the like, and divides signal control time intervals according to the similarity.
The traffic flow and the queuing length are directly reactions to the traffic demand of an entrance way or a flow direction level, but the maximum utilization rate of single-point signal optimization is the optimization target when most of the single-point signal optimization is green in a non-saturation period, and the maximum release rate of the single-point signal optimization is the optimization target when most of the single-point signal optimization is in a saturated supersaturation state; therefore, in performing the period optimization, only the traffic demand is considered, and the release of the signal control phase for each flow direction is ignored, which is not comprehensive for the control period division.
Disclosure of Invention
Aiming at the problems in the background technology, the invention provides a time interval division optimization method for a multi-time interval signal control scheme of an urban road intersection, which aims at queuing and evacuating time from an inlet to a hierarchy, calibrates a phase corresponding to minimum time granularity, an idle discharge rate of a flow direction and an unbalance rate, further analyzes the utilization condition of adjacent time intervals in green time with different flow directions and the requirement unbalance condition of different flow directions in the same phase, and researches and judges the feasibility of combining the adjacent time intervals on the basis. By multiple analysis and combination, the original time interval division scheme is optimized, and the signal optimization target under the unsaturated and saturated supersaturated states of the intersection is met, so that the signal control time interval division with higher reliability and higher stability is realized.
A time interval division optimization method for a multi-time interval signal control scheme of an urban road intersection comprises the following steps:
step 1, butting intersection channelized information, entrance traffic flow detection data and a multi-period traffic signal control scheme;
step 2, calibrating the phase vacancy rate, the phase unbalance rate and the flow direction vacancy rate of the short time interval based on actually measured inlet road traffic flow detection data of the short time interval and a multi-time-interval traffic signal control scheme, and performing data preprocessing to obtain a phase vacancy rate and phase unbalance rate time sequence in a control time interval to be divided;
step 3, extracting a short time interval sequence F to be optimized from the preprocessed phase vacancy rate and phase imbalance rate time sequence in the to-be-divided control time period obtained in the step 2 based on the difference analysis of the green time utilization condition and the different flow direction vehicle passing imbalance conditions of the same phase in the same control time periodk
Step 4, for a short time interval sequence F to be optimizedkIf any two time continuous intervals exist in the interval, performing similarity analysis on the phase free space rate, the phase unbalance rate and the flow direction free space rate in the interval, judging the merging feasibility, and obtaining a secondary merging interval sequence F' to be optimizedk
Step 5, the secondary merging interval sequence F' to be optimized formed after the processing of the step 4kIntervals in the sequence being t', the start moments of the intervals within the sequence being marked as slot cut points p(k,t″)(ii) a Arranging all time-interval dividing points p in time sequence(k,t)、p(k,t″)And forming a plurality of sub-periods in the original period by the dividing points, if the length of the sub-period is less than the time period duration threshold T1Then further detecting the lengths of the front and back adjacent time periods, and canceling the division point between the sub time period and the adjacent time period with smaller length; thereby forming an optimized period PkInner time segmentation point sequence m(k,t)Accordingly, the original period PkCarrying out optimized division to form an optimized scheme Pm(ii) a Wherein the time interval duration threshold T1Usually not less than 15 min;
step 6, for the optimization scheme PmAny two adjacent periods P inm、Pm+1Wherein m denotes an optimization scheme PmIn the middle time period, if the original time periods of the adjacent time periods are different, the running characteristic similarity of the adjacent time periods is analyzed based on the flow direction free space rate, and if the merging condition is met, the adjacent time periods P are subjected to the combinationm、Pm+1Merging, otherwise, keeping the original division point; obtaining the formation optimization scheme Pv
Step 7, extracting an optimization scheme PvMedium length less than time length threshold T2And combining successive ones of the time periods of (a) and (b), a time threshold of duration T2The value of (A) is not less than 30 min; for a period in which there is a dispersion, detecting whether f exists in its adjacent periodtShort time intervals of 1, if present, are compared with the presence of ftCombining the time intervals with short time intervals and short total time length; otherwise, merging the time interval into the time interval with shorter total time length in the adjacent time interval; the combined time period is the final optimized division scheme Pf
Further, in step 1, the channelized information includes a corresponding relationship between the flow direction j of each inlet lane and the lane; the data for detecting the traffic flow of the entrance road comprise the time interval h of each turning saturated locomotivejAnd flow rate q in short time interval(j,t)Wherein t refers to a short time interval, and takes 5min as the short time interval duration; the multi-period signal control scheme comprises a traffic signal control period P of the intersection all daykAnd the start and end time of the period and the signal control cycle C corresponding to the periodkPhase sequence, phase green time duration g(k,i)And calculating the duration G of each flow direction green light according to the flow direction composition condition in the phase(k,j)Where k denotes a control period number, i denotes a phase number, and j denotes a flow direction number in the phase i.
Further, the specific steps of step 2 include:
step 2-1, based on the idle discharge time delta of each flow direction of the intersection(j,t)Calculating the free space rate A of each phase in a short time interval(i,t)Each phase imbalance ratio U(i,t)Flow direction empty rate R(i,t)
Wherein the flow direction is idle for a long time
Figure BDA0002572829390000041
Wherein the short time interval t is included in the control period PkWithin a time interval of (c);
phase free space ratio
Figure BDA0002572829390000042
In the formula
Figure BDA0002572829390000043
The shortest idle time in the flow direction contained in the phase i is pointed;
phase imbalance ratio
Figure BDA0002572829390000044
In the formula
Figure BDA0002572829390000045
The average value of the idle discharge time length of all the flow directions of the phase in the time period is obtained, and N is the flow direction number contained in the phase i;
flow direction empty rate
Figure BDA0002572829390000046
The parameters in the formula are as before;
step 2-2, calculating the obtained empty discharge rate A of each phase according to the step 2-1(i,t)Each phase imbalance ratio U(i,t)Time series of the original traffic signal control period PkJudging the integrity of the data in the time interval, and if the missing proportion of the phase vacancy rate data and the unbalance rate data exceeds a threshold value D, marking the control time interval as not to be segmented; otherwise, the control time interval is marked to be divided, and phase vacancy rate and phase imbalance rate parameter time sequences of all short time intervals in the control time interval are extracted; the value of the threshold D is not more than 50 percent;
step 2-3, performing exponential smoothing processing on the extracted time series of the short-time interval phase vacancy rate and the phase imbalance rate parameters, and determining an exponential smoothing coefficient according to the range of the time series of the phase imbalance rate, wherein the range is different
Figure BDA0002572829390000051
I.e. the difference between the maximum and minimum values of the phase imbalance ratio.
Further, in step 3, a sequence of short time intervals F to be optimizedkComprises the following steps:
step 3-1, if the phase position empty rate A(i,t)>α1And phase imbalance ratio U(i,t)≤μ1Then the optimized characteristic parameter f of the short time interval i is settLabeled 0; otherwise, ftLabeled 1; wherein, the optimization characteristic parameters are defined as:
Figure BDA0002572829390000052
α1mu is the splitting threshold of the phase space ratio and the unbalance ratio respectively, and the value is alpha1=0.1,μ1=0.5;
Step 3-2, extracting all marks ftShort time interval of 1, and further judging
Figure BDA0002572829390000053
Or
Figure BDA0002572829390000054
If yes, extracting the short time interval t, and combining the short time intervals belonging to the same control time interval into an interval sequence F to be optimizedk(ii) a If not, optimizing the characteristic parameter ftThe value is assigned to 0; wherein, FkCorresponding to the original traffic signal control period PkBeta is a change rate threshold value, and adjustment is carried out according to the actual measurement condition of the traffic operation of the intersection;
step 3-3, extracting all ftShort time interval of 0, merging the continuous intervals existing in the short time interval, and marking the start time point of the merged interval as a cutting point p(k,t)
Further, the step 4 comprises the following steps:
step 4-1, for any two consecutive short time intervals t, t +1, if for any phase i, A within the interval(i,t)>α2If both are true, combining the short time intervals t and t +1 to form a primary combined interval sequence F to be optimizedk'; wherein alpha is2A primary merging threshold;
step 4-2, for FkPerforming convolution smoothing treatment on the original phase imbalance rate in each time period;
step 4-3, for Fk' inInterval, if there are any two time-consecutive intervals t ', t' +1, the imbalance ratio of the flow direction satisfies R(j,t')·R(j,t'+1)Greater than 0, and the unbalanced flow direction corresponds to the phase imbalance exponent | U(i,t')-U(i,t'+1)|≤μ2Then the intervals t ', t' +1 are merged to form a secondary merged interval sequence F ″ to be optimizedk(ii) a Wherein t' represents the sequence F after the merging treatment in the step 4-1kTime interval of2The imbalance index similarity threshold is quadratic combined.
Further, in step 6, the merging conditions include:
condition 1, arbitrary short time interval t in two time periodsmFree flow direction free space ratio of
Figure BDA0002572829390000061
Absolute threshold gamma of flow direction free space ratio in formula1The value is usually not less than 0.8;
condition 2, all flow directions in adjacent intervals satisfy
Figure BDA0002572829390000062
Condition 3, all flow directions in adjacent intervals satisfy
Figure BDA0002572829390000063
Relative threshold gamma of flow direction free space ratio in the formula2Values are typically not greater than 0.3;
and merging when the three conditions are simultaneously met.
The invention achieves the following beneficial effects: the traffic operation characteristic similarity is analyzed based on the idle discharge condition and the flow direction unbalance condition of the phase, and the original signal control time interval division scheme is optimized by studying, judging and combining a plurality of conditions with the minimum granularity time interval, so that the common signal optimization targets in the unsaturated and saturated supersaturated states of the intersection can be met, and the applicability is strong; because the optimization is based on the original time interval division scheme, the reliability and stability of the optimized scheme can be guaranteed.
Drawings
Fig. 1 is a flowchart illustrating steps of a time interval division optimizing method according to an embodiment of the present invention.
FIG. 2 is a flowchart illustrating steps in a data collection phase according to an embodiment of the present invention.
FIG. 3 is a flowchart illustrating steps of the data preprocessing stage according to an embodiment of the present invention.
Fig. 4 is a flowchart illustrating steps of an interval extraction phase to be optimized according to an embodiment of the present invention.
Fig. 5 is a flowchart of steps of an interval merging phase to be optimized according to an embodiment of the present invention.
Fig. 6 is a flowchart of steps in the optimization scheme generation phase according to an embodiment of the present invention.
Detailed Description
The technical scheme of the invention is further explained in detail by combining the drawings in the specification.
A time interval division optimization method for a multi-time interval signal control scheme of an urban road intersection comprises the following steps:
s1, connecting intersection channelized information, entrance road traffic flow detection data and a multi-period traffic signal control scheme.
Wherein, the channelized information comprises the corresponding relation between the flow direction j of each inlet channel and the lane; the data for detecting the traffic flow of the entrance road comprise the time interval h of each turning saturated locomotivejAnd flow rate q in short time interval(j,t)Wherein t is a short time interval, usually 5min is used as the short time interval; the multi-period signal control scheme comprises a traffic signal control period P of the intersection all daykStarting and ending time and signal control cycle C corresponding to the time periodkPhase sequence, phase green time duration g(k,i)And extracting the duration G of each flow direction green light according to the time(k,j)Where k denotes a control period number, i denotes a phase number, and j denotes a flow direction number in the phase i.
In the embodiment, taking a common intersection as an example, the value of the flow direction j is [1,12], and the numbering sequence refers to { north straight, north left, north right, south straight, south left, south right, west straight, west left, west right, east straight, east left, east right }.
S2, calibrating the phase vacancy rate, the phase unbalance rate and the flow direction vacancy rate of the short time interval based on the actually measured short time interval entrance road traffic flow detection data and the multi-period traffic signal control scheme, and performing data preprocessing; the method comprises the following specific steps:
s21, empty time length delta based on each flow direction of intersection(j,t)Calculating the free space rate A of each phase in a short time interval(i,t)Each phase imbalance ratio U(i,t)Flow direction empty rate R(i,t)
Wherein the flow direction is idle for a long time
Figure BDA0002572829390000081
Wherein the short time interval t is included in the control period PkWithin a time interval of (a).
Phase free space ratio
Figure BDA0002572829390000082
In the formula
Figure BDA0002572829390000083
The shortest free time in the flow direction included in phase i.
Phase imbalance ratio
Figure BDA0002572829390000084
In the formula
Figure BDA0002572829390000085
And N is the average value of the time length of all the flow directions of the phase in the period, and is the flow direction number contained in the phase i.
Flow direction empty rate
Figure BDA0002572829390000086
The parameters in the formula are as before.
S22, calculating the obtained empty discharge rate A of each phase according to S21(i,t)Each phase imbalance ratio U(i,t)Time series of the original traffic signal control period PkThe data integrity is judged, if the missing proportion of the phase space ratio and the unbalance ratio data exceeds a threshold value D, the phase space ratio and the unbalance ratio data are judged to be out of balanceThe control period is marked as not sliced; otherwise, the control time interval is not marked, and phase vacancy rate and phase imbalance rate parameter time sequences of all short time intervals in the control time interval are extracted; typically, the threshold D takes on a value of no more than 50%.
S23, performing exponential smoothing processing on the extracted time series of the short-time interval phase free space rate and the phase imbalance rate parameter, and determining an exponential smoothing coefficient according to the time series range of the phase imbalance rate, wherein the range is different
Figure BDA0002572829390000091
I.e. the difference between the maximum and minimum values of the phase imbalance ratio.
In the embodiment, the relationship between the range of the phase imbalance ratio and the selection of the exponential smoothing mode and the smoothing coefficient is shown in the following table:
RUi exponential smoothing mode Coefficient of smoothing
≤0.04 First order exponential smoothing 0.1
(0.04,0.15] First order exponential smoothing 0.4
>0.15 Quadratic exponential smoothing 0.9
And S3, extracting a short time interval to be optimized based on the difference analysis of the utilization condition of the same phase green time and the different flow direction vehicle passing unbalance conditions of the same phase green time in the same control period from the preprocessed phase empty rate and phase unbalance rate time sequence in the unmarked control period in the S2. The extraction conditions were:
s31, phase free space ratio A(i,t)>α1And phase imbalance ratio U(i,t)≤μ1Then the optimized characteristic parameter f of the short time interval i is settLabeled 0; the optimization characteristic parameter is defined as:
Figure BDA0002572829390000092
α1mu is the splitting threshold of the phase vacancy rate and the unbalance rate respectively, and is usually taken as alpha1=0.1,μ10.5; otherwise, ftLabeled 1.
S32, extracting all further judgment ftShort time interval of 1, and further judging
Figure BDA0002572829390000093
Or
Figure BDA0002572829390000094
If yes, extracting the short time interval t, and combining the short time intervals belonging to the same control time interval into an interval sequence F to be optimizedk(ii) a In the formula, FkCorresponding to the original traffic signal control period PkBeta is a change rate threshold value, the value in the embodiment is 0.15, and the adjustment can be carried out according to the actual measurement condition of the traffic operation of the intersection; if not, optimizing the characteristic parameter ftThe value is assigned to 0.
S33, extracting all ftShort time interval of 0, merging the continuous intervals existing in the short time interval, and marking the start time point of the merged interval as a cutting point p(k,t)
S4, for interval sequence F to be optimizedkIf there are any two time-consecutive intervalsIf the interval is not the same, carrying out similarity analysis on the phase free space rate, the phase unbalance rate and the flow direction free space rate in the interval, and judging the feasibility of combination; the method comprises the following specific steps:
s41, for any two consecutive short time intervals t, t +1, if for any phase i, A in the interval(i,t)>α2If both are true, combining the short time intervals t and t +1 to form a primary combined interval sequence F to be optimizedk'; in the formula, alpha2A primary merging threshold; in an embodiment, the threshold α2The value is 0.8.
S42, pair Fk' the original phase imbalance ratio of each time interval is processed by convolution smoothing.
S43 for FkThe interval in (1), if there are any two time-consecutive intervals t ', t' +1, the flow imbalance ratio satisfies R(j,t')·R(j,t'+1)Greater than 0, and the unbalanced flow direction corresponds to the phase imbalance exponent | U(i,t')-U(i,t'+1)|≤μ2Then the intervals t ', t' +1 are merged to form a secondary merged interval sequence F ″ to be optimizedk(ii) a Wherein t' represents the sequence F after S41 combinationkTime interval of2The imbalance index similarity threshold is combined for the second time; in the examples, μ2The value is 0.2.
S5 time period F to be optimizedkSequence F ″, formed after treatment with S4kIntervals in the sequence being t', the start moments of the intervals within the sequence being marked as slot cut points p(k,t″)(ii) a Arranging all time-interval dividing points p in time sequence(k,t)、p(k,t″)And forming a plurality of sub-periods in the original period by the dividing points, if the length of the sub-period is less than the time period duration threshold T1Then further detecting the lengths of the front and back adjacent time periods, and canceling the division point between the sub time period and the adjacent time period with smaller length; thereby forming an optimized period PkInner time segmentation point sequence m(k,t)Accordingly, the original period PkCarrying out optimized division to form an optimized scheme Pm(ii) a In the formula, the time interval duration threshold T1Usually not less than 15 min.
S6 for renTwo adjacent periods Pm、Pm+1Wherein m denotes an optimization scheme PmAnd if the original time periods of the time periods in the time period are different, performing characteristic similarity analysis on the operation of the time periods based on the flow direction free space rate.
T within any short time interval in two time periodsmFree flow direction free space ratio of
Figure BDA0002572829390000111
Absolute threshold gamma of flow direction free space ratio in formula1The value is usually not less than 0.8; in an embodiment the threshold value is 0.8.
All flow directions in adjacent intervals satisfy
Figure BDA0002572829390000112
All flow directions in adjacent intervals all satisfy
Figure BDA0002572829390000113
Relative threshold gamma of flow direction free space ratio in the formula2Values are typically not greater than 0.3; in an embodiment the threshold value is 0.2.
If the above conditions are satisfied simultaneously, the adjacent time interval P is definedm、Pm+1Merging; otherwise, keeping the original division point; form an optimization plan Pv
S7, extracting an optimization scheme PvMedium length less than time length threshold T2And combining successive ones thereof, typically by a time threshold T2The value of (A) is not less than 30 min; for a period in which there is a dispersion, detecting whether f exists in its adjacent periodtShort time intervals of 1, if present, are compared with the presence of ftCombining the time intervals with short time intervals and short total time length; otherwise, merging the time interval into the time interval with shorter total time length in the adjacent time interval; the combined time period is the final optimized division scheme Pf
The above description is only a preferred embodiment of the present invention, and the scope of the present invention is not limited to the above embodiment, but equivalent modifications or changes made by those skilled in the art according to the present disclosure should be included in the scope of the present invention as set forth in the appended claims.

Claims (6)

1. A time interval division optimization method for a multi-time interval signal control scheme of an urban road intersection is characterized by comprising the following steps: the method comprises the following steps:
step 1, butting intersection channelized information, entrance traffic flow detection data and a multi-period traffic signal control scheme;
step 2, calibrating the phase vacancy rate, the phase unbalance rate and the flow direction vacancy rate of the short time interval based on actually measured inlet road traffic flow detection data of the short time interval and a multi-time-interval traffic signal control scheme, and performing data preprocessing to obtain a phase vacancy rate and phase unbalance rate time sequence in a control time interval to be divided;
step 3, extracting a short time interval sequence F to be optimized from the preprocessed phase vacancy rate and phase imbalance rate time sequence in the to-be-divided control time period obtained in the step 2 based on the difference analysis of the green time utilization condition and the different flow direction vehicle passing imbalance conditions of the same phase in the same control time periodk
Step 4, for a short time interval sequence F to be optimizedkIf any two time continuous intervals exist in the interval, performing similarity analysis on the phase free space rate, the phase unbalance rate and the flow direction free space rate in the interval, judging the merging feasibility, and obtaining a secondary merging interval sequence F' to be optimizedk
Step 5, forming a secondary merging interval sequence F to be optimized after the processing of the step 4k", the intervals in the sequence are t", the starting point instant of each interval in the sequence is marked as a period segmentation point p(k,t”)(ii) a Arranging all time-interval dividing points p in time sequence(k,t)、p(k,t”)And forming a plurality of sub-periods in the original period by the dividing points, if the length of the sub-period is less than the time period duration threshold T1Then further detecting the lengths of the front and back adjacent time periods, and canceling the division point between the sub time period and the adjacent time period with smaller length; thereby forming an optimized period PkInner time segmentation point sequence m(k,t)Accordingly, the original period PkCarrying out optimized division to form an optimized scheme Pm(ii) a Wherein the time interval duration threshold T1Usually not less than 15 min;
step 6, for the optimization scheme PmAny two adjacent periods P inm、Pm+1Wherein m denotes an optimization scheme PmIn the middle time period, if the original time periods of the adjacent time periods are different, the running characteristic similarity of the adjacent time periods is analyzed based on the flow direction free space rate, and if the merging condition is met, the adjacent time periods P are subjected to the combinationm、Pm+1Merging, otherwise, keeping the original division point; obtaining the formation optimization scheme Pv
Step 7, extracting an optimization scheme PvMedium length less than time length threshold T2And combining successive ones of the time periods of (a) and (b), a time threshold of duration T2The value of (A) is not less than 30 min; for a period in which there is a dispersion, detecting whether f exists in its adjacent periodtShort time intervals of 1, if present, are compared with the presence of ftCombining the time intervals with short time intervals and short total time length; otherwise, merging the time interval into the time interval with shorter total time length in the adjacent time interval; the combined time period is the final optimized division scheme Pf
2. The method for optimizing time interval division of the urban road intersection multi-time interval signal control scheme according to claim 1, wherein: in the step 1, channelizing information comprises the corresponding relation between the flow direction j of each inlet channel and a lane; the data for detecting the traffic flow of the entrance road comprise the time interval h of each turning saturated locomotivejAnd flow rate q in short time interval(j,t)Wherein t refers to a short time interval, and takes 5min as the short time interval duration; the multi-period signal control scheme comprises a traffic signal control period P of the intersection all daykAnd the start and end time of the period and the signal control cycle C corresponding to the periodkPhase sequence, phase green time duration g(k,i)And calculating the duration G of each flow direction green light according to the flow direction composition condition in the phase(k,j)Where k denotes a control period number, i denotes a phaseThe number j refers to the flow number in phase i.
3. The method for optimizing time interval division of the urban road intersection multi-time interval signal control scheme according to claim 1, wherein: the specific steps of step 2 include:
step 2-1, based on the idle discharge time delta of each flow direction of the intersection(j,t)Calculating the free space rate A of each phase in a short time interval(i,t)Each phase imbalance ratio U(i,t)Flow direction empty rate R(i,t)
Wherein the flow direction is idle for a long time
Figure FDA0002572829380000031
Wherein the short time interval t is included in the control period PkWithin a time interval of (c);
phase free space ratio
Figure FDA0002572829380000032
In the formula
Figure FDA0002572829380000033
The shortest idle time in the flow direction contained in the phase i is pointed;
phase imbalance ratio
Figure FDA0002572829380000034
In the formula
Figure FDA0002572829380000035
The average value of the idle discharge time length of all the flow directions of the phase in the time period is obtained, and N is the flow direction number contained in the phase i;
flow direction empty rate
Figure FDA0002572829380000036
The parameters in the formula are as before;
step 2-2, calculating the obtained empty discharge rate A of each phase according to the step 2-1(i,t)Each phase imbalance ratio U(i,t)Time series of the original traffic signalPeriod PkJudging the integrity of the data in the time interval, and if the missing proportion of the phase vacancy rate data and the unbalance rate data exceeds a threshold value D, marking the control time interval as not to be segmented; otherwise, the control time interval is marked to be divided, and phase vacancy rate and phase imbalance rate parameter time sequences of all short time intervals in the control time interval are extracted; the value of the threshold D is not more than 50 percent;
step 2-3, performing exponential smoothing processing on the extracted time series of the short-time interval phase vacancy rate and the phase imbalance rate parameters, and determining an exponential smoothing coefficient according to the range of the time series of the phase imbalance rate, wherein the range is different
Figure FDA0002572829380000037
I.e. the difference between the maximum and minimum values of the phase imbalance ratio.
4. The method for optimizing time interval division of the urban road intersection multi-time interval signal control scheme according to claim 1, wherein: in step 3, a short time interval sequence F to be optimizedkComprises the following steps:
step 3-1, if the phase position empty rate A(i,t)>α1And phase imbalance ratio U(i,t)≤μ1Then the optimized characteristic parameter f of the short time interval i is settLabeled 0; otherwise, ftLabeled 1; wherein, the optimization characteristic parameters are defined as:
Figure FDA0002572829380000041
α1mu is the splitting threshold of the phase space ratio and the unbalance ratio respectively, and the value is alpha1=0.1,μ1=0.5;
Step 3-2, extracting all marks ftShort time interval of 1, and further judging
Figure FDA0002572829380000042
Or
Figure FDA0002572829380000043
Whether or not to becomeAnd if so, extracting the short time interval t, and combining the short time intervals belonging to the same control time period into an interval sequence F to be optimizedk(ii) a If not, optimizing the characteristic parameter ftThe value is assigned to 0; wherein, FkCorresponding to the original traffic signal control period PkBeta is a change rate threshold value, and adjustment is carried out according to the actual measurement condition of the traffic operation of the intersection;
step 3-3, extracting all ftShort time interval of 0, merging the continuous intervals existing in the short time interval, and marking the start time point of the merged interval as a cutting point p(k,t)
5. The method for optimizing time interval division of the urban road intersection multi-time interval signal control scheme according to claim 1, wherein: the step 4 comprises the following steps:
step 4-1, for any two consecutive short time intervals t, t +1, if for any phase i, A within the interval(i,t)>α2If all the time intervals are satisfied, merging the short time intervals t and t +1 to form a primary merging interval sequence F 'to be optimized'k(ii) a Wherein alpha is2A primary merging threshold;
step 4-2, p of F'kCarrying out convolution smoothing treatment on the original phase imbalance rate in each time interval;
step 4-3, for F'kIf there are any two time-consecutive intervals t ', t' +1, the imbalance ratio of the flow direction satisfies R(j,t')·R(j,t'+1)Greater than 0, and the unbalanced flow direction corresponds to the phase imbalance exponent | U(i,t')-U(i,t'+1)|≤μ2Then the intervals t ', t' +1 are merged to form a secondary merged interval sequence F ″ to be optimizedk(ii) a Wherein t ' represents a sequence F ' after the merging treatment of the step 4-1 'kTime interval of (1), mu2The imbalance index similarity threshold is quadratic combined.
6. The method for optimizing time interval division of the urban road intersection multi-time interval signal control scheme according to claim 1, wherein: in step 6, the merging conditions include:
condition 1, arbitrary short time interval t in two time periodsmFree flow direction free space ratio of
Figure FDA0002572829380000051
Absolute threshold gamma of flow direction free space ratio in formula1The value is usually not less than 0.8;
condition 2, all flow directions in adjacent intervals satisfy
Figure FDA0002572829380000052
Condition 3, all flow directions in adjacent intervals satisfy
Figure FDA0002572829380000053
Relative threshold gamma of flow direction free space ratio in the formula2Values are typically not greater than 0.3;
and merging when the three conditions are simultaneously met.
CN202010645121.7A 2020-07-07 2020-07-07 Time interval division optimization method for multi-time interval signal control scheme of urban road intersection Active CN111951572B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010645121.7A CN111951572B (en) 2020-07-07 2020-07-07 Time interval division optimization method for multi-time interval signal control scheme of urban road intersection

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010645121.7A CN111951572B (en) 2020-07-07 2020-07-07 Time interval division optimization method for multi-time interval signal control scheme of urban road intersection

Publications (2)

Publication Number Publication Date
CN111951572A true CN111951572A (en) 2020-11-17
CN111951572B CN111951572B (en) 2022-11-08

Family

ID=73341064

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010645121.7A Active CN111951572B (en) 2020-07-07 2020-07-07 Time interval division optimization method for multi-time interval signal control scheme of urban road intersection

Country Status (1)

Country Link
CN (1) CN111951572B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112991729A (en) * 2021-02-25 2021-06-18 杭州海康威视数字技术股份有限公司 Time interval dividing method and device and computer storage medium
CN114999184A (en) * 2022-04-21 2022-09-02 上海商汤智能科技有限公司 Method, device, equipment, storage medium and program product for updating signal time interval

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016091207A (en) * 2014-10-31 2016-05-23 住友電気工業株式会社 Traffic signal controller, computer program, and traffic signal control method
CN105608510A (en) * 2015-12-31 2016-05-25 连云港杰瑞电子有限公司 Traffic period automatic division method based on Fisher algorithm
CN105788236A (en) * 2014-12-26 2016-07-20 浙江大华技术股份有限公司 Traffic control method and traffic control device
CN105976062A (en) * 2016-05-13 2016-09-28 腾讯科技(深圳)有限公司 Signal lamp duration data mining method, and travel service implementation method and device
CN109087519A (en) * 2018-08-10 2018-12-25 江苏智通交通科技有限公司 Joint control crossing traffic signaling plan divided stages and it is green when configuration method
CN109410574A (en) * 2018-10-29 2019-03-01 东南大学 A kind of timing parameter optimization method towards stage-phase signal control program
CN110136444A (en) * 2019-06-06 2019-08-16 南京慧尔视智能科技有限公司 A kind of intersection green light sky puts time calculation method and device
CN110189529A (en) * 2019-03-21 2019-08-30 江苏智通交通科技有限公司 Traffic signal control scheme switching adjustment management system and method
CN110956804A (en) * 2019-11-19 2020-04-03 中电海康集团有限公司 Intersection signal timing diagnosis method
CN111161550A (en) * 2020-01-02 2020-05-15 合肥学院 Intersection control time period and phase timing design and closed-loop verification method
CN111192465A (en) * 2020-01-07 2020-05-22 上海宝康电子控制工程有限公司 Method for realizing signal timing scheme group division processing based on flow data

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016091207A (en) * 2014-10-31 2016-05-23 住友電気工業株式会社 Traffic signal controller, computer program, and traffic signal control method
CN105788236A (en) * 2014-12-26 2016-07-20 浙江大华技术股份有限公司 Traffic control method and traffic control device
CN105608510A (en) * 2015-12-31 2016-05-25 连云港杰瑞电子有限公司 Traffic period automatic division method based on Fisher algorithm
CN105976062A (en) * 2016-05-13 2016-09-28 腾讯科技(深圳)有限公司 Signal lamp duration data mining method, and travel service implementation method and device
CN109087519A (en) * 2018-08-10 2018-12-25 江苏智通交通科技有限公司 Joint control crossing traffic signaling plan divided stages and it is green when configuration method
CN109410574A (en) * 2018-10-29 2019-03-01 东南大学 A kind of timing parameter optimization method towards stage-phase signal control program
CN110189529A (en) * 2019-03-21 2019-08-30 江苏智通交通科技有限公司 Traffic signal control scheme switching adjustment management system and method
CN110136444A (en) * 2019-06-06 2019-08-16 南京慧尔视智能科技有限公司 A kind of intersection green light sky puts time calculation method and device
CN110956804A (en) * 2019-11-19 2020-04-03 中电海康集团有限公司 Intersection signal timing diagnosis method
CN111161550A (en) * 2020-01-02 2020-05-15 合肥学院 Intersection control time period and phase timing design and closed-loop verification method
CN111192465A (en) * 2020-01-07 2020-05-22 上海宝康电子控制工程有限公司 Method for realizing signal timing scheme group division processing based on flow data

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
熊睿成等: "单点交叉口信号控制方案时段划分方法", 《第十四届中国智能交通年会论文集》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112991729A (en) * 2021-02-25 2021-06-18 杭州海康威视数字技术股份有限公司 Time interval dividing method and device and computer storage medium
CN112991729B (en) * 2021-02-25 2022-05-20 杭州海康威视数字技术股份有限公司 Time interval dividing method and device and computer storage medium
CN114999184A (en) * 2022-04-21 2022-09-02 上海商汤智能科技有限公司 Method, device, equipment, storage medium and program product for updating signal time interval

Also Published As

Publication number Publication date
CN111951572B (en) 2022-11-08

Similar Documents

Publication Publication Date Title
CN108596202B (en) Method for calculating personal commuting time based on mobile terminal GPS positioning data
CN111951572B (en) Time interval division optimization method for multi-time interval signal control scheme of urban road intersection
CN108109382B (en) Congestion point, congestion line and congestion area discovery method based on composite network
CN104464304A (en) Urban road vehicle running speed forecasting method based on road network characteristics
CN106408939A (en) Traffic flow sequence classification method based on density peak value clustering
CN109377753B (en) Trunk line coordination optimization method for repeatedly releasing in coordination direction
CN110751311B (en) Data extraction and real-time prediction method for sporadic traffic jam duration
CN113222382B (en) Method for determining passing capacity of heterogeneous traffic flow lane change influence road sections in Internet of vehicles environment
CN112508237B (en) Rain type region division method based on data analysis and real-time rain type prediction method
CN109754598A (en) A kind of congestion is formed a team recognition methods and system
CN106683441B (en) Intersection signal timing scheme evaluation method
CN111325392B (en) Precipitation forecasting system in tropical cyclone process
CN106683442A (en) Multi-index based intersection signal timing plan evaluation method
CN109697849A (en) Intelligent transportation Time segments division method based on sliding average algorithm
CN108492560A (en) A kind of Road Detection device missing data complementing method and device
CN110827537B (en) Method, device and equipment for setting tidal lane
CN106846808B (en) A kind of vehicle parking based on license plate data time number calculating method
CN112652164B (en) Traffic time interval dividing method, device and equipment
CN104715160A (en) Soft measurement modeling data outlier detecting method based on KMDB
CN116704750B (en) Traffic state identification method based on clustering algorithm, electronic equipment and medium
CN117173889A (en) Urban road multi-intersection and associated road section comprehensive evaluation method
CN115100847B (en) Queuing service time estimation method for low-permeability network-connected track data
CN113034940A (en) Fisher ordered clustering-based single-point signalized intersection optimization timing method
CN110119884A (en) A kind of high-speed railway passenger flow Time segments division method based on neighbour's propagation clustering
CN115691170A (en) Urban trunk road subregion division method based on number plate data

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant