CN104408916B - Based on section speed, the road traffic running status appraisal procedure of data on flows - Google Patents

Abstract

The invention discloses a kind of based on section speed, the road traffic running status appraisal procedure of data on flows, the vehicle checker speed flowrate data that section to be assessed respectively contrasts period are obtained first, using speed flowrate relation curve, it is determined that the bound of integration；And will integrate total interval division is that multiple minizone sections are processed, the accounting of total quantity as the weight of each interval integral value is accounted for using each segment data volume simultaneously, obtain respective traffic circulation state characterization value, relative superior or inferior is judged finally by comparative analysis characterization value, and obtains quantitative analysis results.Qualitative to traffic circulation state and quantitative comparative analysis is realized, from the operational efficiency in the ensemble average effect of traffic circulation state come concentrated expression section period on year-on-year basis.Using the overlying relation and its respective variation tendency of each bar speed-discharge relation curve for participating in comparative analysis come the relative superior or inferior of the intuitively qualitative respective traffic circulation state of reflection；It is actual workable, it is simple and easy to apply, it is readily appreciated that.

Description

Based on section speed, the road traffic running status appraisal procedure of data on flows

Technical field

The present invention relates to traffic circulation state analysis field, more particularly to a kind of road traffic running status appraisal procedure.

Background technology

As economic develops rapidly, the operation scale of urban transportation and highway is increasing, the operation shape of traffic Condition also becomes increasingly conspicuous for the influence degree of city overall operation efficiency, but Current traffic has a lot, such as traffic congestion Situation is increasingly severe, and traffic environment goes from bad to worse, and traffic route utilizes unreasonable etc., and this is needed in actively development traffic circulation While condition monitoring system, traffic noise prediction assay method is set up, additionally, the analysis of road section traffic volume operation conditions is commented Valency is effectively management and operation road, gives full play to the key of its potentiality, and scientific and effective evaluation traffic noise prediction can be fortune Administrative department of battalion provides analysis, management, the method and foundation of decision-making, is favorably improved section operation management level and service water It is flat.

The current evaluation to road traffic operation conditions, is all based on reflecting that many aspects of road traffic operational efficiency set Count evaluation index to characterize operation conditions, conventional evaluation index is such as：The volume of traffic, saturation degree, instantaneous velocity, average traveling speed Degree, average travel speed, traffic density, time headway and space headway, occupation rate, average travel time, average stroke delay, Link traffic congested coefficient or crowding etc..Also, majority research is by single desired value or based on many of integrated evaluating method Index comprehensive evaluation, carrys out the traffic noise prediction of comparative analysis section different times.But it is so relatively cumbersome, and show slightly inadequate Intuitively.

Therefore need it is a kind of can the concentrated expression section operational efficiency in period on year-on-year basis state evaluating method.

The content of the invention

In view of this, the technical problems to be solved by the invention are to provide a kind of road based on section speed, data on flows Road traffic circulation state appraisal procedure.

The object of the present invention is achieved like this：

It is based on section speed, the road traffic running status appraisal procedure of data on flows including following that the present invention is provided Step：

S1：Obtain section to be assessed and respectively contrast the vehicle checker data in period, and data standard is stored in tables of data In, the data item in the tables of data includes speed, the flow in each track；

S2：Vehicle checker to obtaining carries out data prediction, and each track data are processed respectively, while by after treatment It is stored in tables of data data standard；The data prediction includes two large divisions：Vehicle checker vehicle flowrate is converted into mark Quasi- car vehicle flowrate, and rejecting abnormalities data；

S3：Using by the actual vehicle checker speed in the section after data prediction and data on flows, implement section to be assessed The fitting of different times speed-discharge relation curve；Be fitted in the form of following conic section the actual speed in section- Discharge relation：

V=AQ²+BQ+C

In formula：V represents speed；Q represents flow；A, B, C are regression coefficient；

S4：It is determined that to the analysis mode of road traffic running status assessment, if carrying out qualitative analysis, can be according to curve Position relationship and its respective variation tendency provide the qualitative analysis, and terminate road traffic running status evaluation process； If carrying out quantitative analysis, go to step S5 and implement quantitative analysis；

S5：The quantitative analysis utilizes the speed-discharge relation curve for obtaining, according to effective data area, it is determined that product The bound divided；And it is that multiple minizone sections are processed that whole flow is integrated into total interval division, is obtained using integral operation Obtain quantitative analysis results.

Further, the quantitative analysis in the step S5 is comprised the following steps that：

S51：The detection performance of design speed per hour and vehicle checker according to section determines the total area of integration in whole flow rate zone Between；

S52：Away from total interval division will be integrated it is multiple minizones section by certain group, obtains multiple integrating ranges sections, and Using each interval terminal as respective range of integration value；

S53：Data sample to each curve carries out statistical analysis, it is determined that each data volume accounts for sample in the section of each by stages The ratio p of sum_i, and with ratio p_iAs the weight of each segment integrated value；

S54：For the curve of section different times, the multiple segments to above-mentioned division implement integration fortune along curve respectively Calculate, the traffic circulation state characterization value or mutual traffic circulation state in each segment are obtained with respect to characterization value, by following Formula defines two indices to characterize：

(1) traffic circulation state difference characterization value Δ S=S₂-S₁；

(2) traffic circulation state relative different ratio

Wherein, S₁、S₂Each curve corresponding traffic circulation state characterization value in certain flow is interval is represented respectively；Using The quantitative measuring and calculating of implementation is integrated to flow along curve in certain flow is interval, its quantitative measuring and calculating uses below equation：

In formula, Q_b1Represent upper limit of integral value of the contrast period 1 in certain flow interval；Q_a1Represent contrast period 1 one Lower limit of integral value in constant flow interval；A₁Represent the speed-secondary term coefficient of discharge relation formula in contrast period 1；B₁It is right to represent Than the speed-discharge relation formula Monomial coefficient in period 1；C₁Represent the speed-discharge relation formula constant term in contrast period 1； Q_b2Represent upper limit of integral value of the contrast period 2 in certain flow interval；Q_a2Represent contrast period 2 in certain flow is interval Lower limit of integral value；A₂Represent the speed-secondary term coefficient of discharge relation formula in contrast period 2；B₂Represent the speed in contrast period 2 Degree-discharge relation formula Monomial coefficient；C₂Represent the speed-discharge relation formula constant term in contrast period 2；

S55：To the curve of section different times, using respective ratio p_iTraffic circulation state to each segment is related Quantitatively characterizing value weighted calculation, obtains section different times to be assessed each speed-corresponding traffic circulation of discharge relation curve State correlation table value indicative or mutual traffic circulation state are with respect to characterization value (S1, S2, Δ S, Δ_O)；Be given in the following manner Quantitative analysis results：

(1) according to integration formula and ratio p_iIt is calculated the corresponding traffic circulation state correlation table value indicative of each curve, institute State traffic circulation state characterization value S of the traffic circulation state correlation table value indicative including different times₁、S₂, traffic circulation state difference Characterization value Δ S and traffic circulation state relative different ratio of Δ_O；

(2) according to Δ S or Δ_OWith the relation of " 0 ", the quality relatively of section different times traffic circulation state is differentiated, and Quantitative analysis result；

(3) if having selected multiple track comparative analyses, after completing the quantitative analysis in each track by above-mentioned steps, ask for each Track quantitative analysis results Δ S or Δ_OAverage.

Further, the vehicle checker vehicle flowrate in the step S2 is converted to standard vehicle vehicle flowrate is carried out as follows：

In formula, Q is the standard vehicle vehicle flowrate after conversion；q_iIt is i-th kind of flow value of vehicle；w_iIt is i-th kind of conversion of vehicle Coefficient；M is vehicle quantity.

Further, the rejecting abnormalities data in the step S2 are carried out according to the following steps：

S21：Reject the data of apparent error in the following manner according to speed, the actual effectiveness of flow, it is specific as follows：

Data to following a few class situations are rejected：

Reject speed and be more than 0, but flow is equal to 0 data；

Reject speed and be equal to 0, but data of the flow more than 0；

S22：The rejecting treatment of abnormal data is done using the average and standard deviation of vehicle checker speed data, specific rules are such as Under：

1) if there isIt is then valid data, not rejects；

2) if there isOrThen judge v_iThe corresponding data is abnormal data, should Give rejecting；Wherein, 3 σ are screening criteria, represent that the speed normal distribution three times in the section vehicle checker data for obtaining actual measurement are square Difference；It is the vehicle checker speed data average in the period to be analyzed, standard deviation is σ_v, actual vehicle checker velocity amplitude is v_i。

Further, in the step S51 the total interval of integration determines according to the following steps：

S511：It is determined that the total interval upper limit of integration, to carrying out statistical analysis respectively for being fitted the sample data of each curve, Flow maximum is obtained, as the higher limit of whole flow rate zone；

S512：It is determined that the total interval lower limit of integration, specific as follows：

Note section Maximum speed limit is V0, and speed is the longitudinal axis, at the Maximum speed limit V0 from the longitudinal axis, makees straight line V=V0；Utilize Maximum speed limit come aid in determine lower limit of integral；

When straight line V=V0 intersects with two curves to be analyzed in low flow volume region, then total interval lower limit is integrated It is Maximum speed limit V0；

When tops of the straight line V=V0 in two curves to be analyzed, then lower limit of integral value sets in the following manner：Obtain If wagon flow value is Q0 in vehicle checker detection performance requirement unit detection time, make straight line Q=Q0；Respectively with two curve intersections Point, as the corresponding lower limit of integral of each curve；

When straight line V=V0 intersects with two curves, if the joining on curve is at low discharge, respective friendship is taken The corresponding flow value of point is used as lower limit of integral value；If set lower limit owes reasonable, acquisition sets vehicle checker detection performance It is required that wagon flow value is Q0 in unit detection time, make straight line Q=Q0；Respectively with two curve intersection points, as each curve pair The lower limit of integral answered；

When straight line V=V0 is in low flow volume region and a certain bar curve intersection, then take the flow value and make with the average of Q0 It is lower limit of integral value；

Vehicle checker inspection is set when straight line V=V0 is obtained when low flow volume region is located at below two curves or is not intersected, then It is Q0 to survey wagon flow value in performance requirement unit detection time, makees straight line Q=Q0；Respectively with two curve intersection points, as each The corresponding lower limit of integral of curve；

There is no actual data point at the lower limit of integral value for determining as stated above, then judge whether the lower limit of integral value is big In Q0, if judged result is " being more than ", the flow minimum in pretreated data sample of learning from else's experience is used as lower limit of integral value；

It is Q0 otherwise to obtain and set wagon flow value in vehicle checker detection performance requirement unit detection time, makees straight line Q=Q0；Point Not with two curve intersection points, as the corresponding lower limit of integral of each curve；If lower limit of integral value does not have real data when being Q0 yet Point, then the flow minimum in pretreated data sample of learning from else's experience is used as lower limit of integral value.

It is further, described that each total interval use of the corresponding integration of speed-discharge relation curve will be whole in the period of different The mode of individual flow rate zone segment processing is calculated, specific as follows：

If 1) need not ensure that each by stages Duan Jun is contrasted under equal flow, each contrast period curve product That divides is total interval different, and the integration that the method according to step S51 asks for each curve is always interval.

[Q_a(i),Q_b(i)]；

Wherein, i=1,2 ..., represent that the curve integration of different times is total interval, that is, contrast the curve integrating range in period i.

2) if desired ensure that each by stages Duan Jun is contrasted under equal flow, i.e., each contrast period curve integration Take same integration total interval.The corresponding interval of each curve is obtained, each smaller lower limit in interval is selected as under the total interval of integration Limit, selects larger higher limit as total interval upper limit is integrated, and the total interval of the integration of each curve is taken as into the total interval of same integration, Total interval [the Q of same integration_a,Q_b] establish rules really, it is：

If A. Q_a1≥Q_a2, and have Q_a1≤Q_a2, then taking the total interval of integration is：[Q_a2,Q_b2]；

If B. Q_a1＞ Q_a2, and have Q_b1＞ Q_b2, then taking the total interval of integration is：[Q_a2,Q_b1]；

If C. Q_a1＜ Q_a2, and have Q_b1＜ Q_b2, then taking the total interval of integration is：[Q_a1,Q_b2]；

If D. Q_a1≤Q_a2, and have Q_b1≥Q_b2, then taking the total interval of integration is：[Q_a1,Q_b1]。

Further, the total interval [Q of the integration_a(i),Q_b(i)] (or [Q_a,Q_b]) by stages is determined according to the following steps Bound：

Each curve is integrated into total interval [Q away from m by certain group_a(i),Q_b(i)] (or the total interval [Q of same integration_a,Q_b]) draw Be divided into the small by stages sections of n, each small by stages section using each interval terminal as respective range of integration value, specifically It is as follows：

1) if each total interval of curve integration is not to take same interval, point integrating range for each dividing is as follows：

[Q_a(i),Q_a(i)+m]；

[Q_a(i)+m,Q_a(i)+2m]；

[Q_a(i)+2m,Q_a(i)+3m]；

……；

[Q_a(i)+(n-2)m,Q_a(i)+(n-1)m]；

[Q_a(i)+(n-1)m,Q_b(i)]。

Wherein, i=1,2 ..., represent the curve i of above-mentioned institute by stages correspondence different times, i.e., above-mentioned institute by stages is right Than each point of integrating range section that period i curve is integrated.

If 2) each curve integrates total interval and takes same interval [Q_a,Q_b], then point integrating range for dividing is as follows：

[Q_a,Q_a+m]；

[Q_a+m,Q_a+2m]；

[Q_a+2m,Q_a+3m]；

……；

[Q_a+(n-2)m,Q_a+(n-1)m]；

[Q_a+(n-1)m,Q_b]。

The beneficial effects of the present invention are：Present invention employs a kind of road traffic based on section speed, data on flows Running status appraisal procedure, it is bent using the speed obtained based on section vehicle checker data (speed, flow)-flow actual relationship Line carrys out comprehensive analysis traffic circulation state, comparative analysis qualitative to traffic circulation state and quantitative is realized, from traffic circulation shape The operational efficiency in the ensemble average effect of state comes concentrated expression section period on year-on-year basis.Speed based on section vehicle checker data acquisition Degree-flow actual relationship curve, realizes comparative analysis qualitative to traffic circulation state and quantitative.Using participation to score It is each that the overlying relation and its respective variation tendency of each bar speed-discharge relation curve of analysis carry out intuitively qualitative reflection From the relative superior or inferior of traffic circulation state；For quantitative analysis, press certain in speed-flow coordinate plane according to each curve The certain area of regular (effective ranges of data) interception carrys out comparative analysis, first to curve in certain Limit of J-validity (flow rate zone) is divided into multiple segments, and segmentation is integrated along curve to flow, while being accounted for each segment data volume The accounting of total quantity then obtains respective traffic circulation as the weight of each interval integral value to each subsection integral value weighted calculation State table value indicative, relative superior or inferior is judged finally by comparative analysis characterization value, and obtain quantitative analysis results.Analysis result can To aid in illustrating the validity of other assessment mode (such as Utilization assessment metrics evaluation) assessment results, it is also possible to relief road traffic Fortune pipe personnel are to the running status of road, history management and control measures are evaluated, and are that daily management, maintenance plan and planning construction are carried For reference, and then guidance management decision-making, formulate or optimization management and control measures, promote the raising of operation management level and service level.

Using the speed-discharge relation curve obtained based on section vehicle checker speed, data on flows come analysis and assessment section The traffic circulation state of different times, more intuitively can on the whole reflect traffic circulation state, and without analyzing multinomial commenting Valency index, and the validity of other assessment mode (such as Utilization assessment metrics evaluation) assessment results can be aided in illustrating, additionally, actual It is workable, it is simple and easy to apply, it is readily appreciated that.

Brief description of the drawings

In order that the object, technical solutions and advantages of the present invention are clearer, below in conjunction with accompanying drawing the present invention is made into The detailed description of one step, wherein：

Fig. 1 is the road traffic running status assessment side based on section speed, data on flows provided in an embodiment of the present invention Method implementing procedure figure；

Fig. 2 is provided in an embodiment of the present invention based on section speed, the speed of data on flows-flow actual relationship curve (descending from Chongqing highway Yu Wu roads)；

Fig. 3 is section speed provided in an embodiment of the present invention-discharge relation tracing analysis schematic diagram.

Specific embodiment

Hereinafter with reference to accompanying drawing, the preferred embodiments of the present invention are described in detail.It should be appreciated that preferred embodiment Only for the explanation present invention, rather than in order to limit the scope of the invention.

Embodiment 1

As shown in Figure 1-Figure 3, what the present invention was provided is assessed based on section speed, the road traffic running status of data on flows Method, comprises the following steps：

S1：Obtain section to be assessed and respectively contrast the vehicle checker data in period, and data standard is stored in tables of data In, the data item in the tables of data includes speed, the flow in each track；

S2：Vehicle checker to obtaining carries out data prediction, and each track data are processed respectively, while by after treatment It is stored in tables of data data standard；The data prediction includes two large divisions：Vehicle checker vehicle flowrate is converted into mark Quasi- car vehicle flowrate, and rejecting abnormalities data；

S3：Using by the actual vehicle checker speed in the section after data prediction and data on flows, implement section to be assessed The fitting of different times speed-discharge relation curve；Be fitted in the form of following conic section the actual speed in section- Discharge relation：

V=AQ²+BQ+C

In formula：V represents speed；Q represents flow；A, B, C are regression coefficient；

S4：It is determined that to the analysis mode of road traffic running status assessment, if carrying out qualitative analysis, can be according to curve Position relationship and its respective variation tendency provide the qualitative analysis, and terminate road traffic running status evaluation process； If carrying out quantitative analysis, go to step S5 and implement quantitative analysis；

S5：The quantitative analysis utilizes the speed-discharge relation curve for obtaining, according to effective data area, it is determined that product The bound divided；And it is that multiple minizone sections are processed that whole flow is integrated into total interval division, is obtained using integral operation Obtain quantitative analysis results.

Quantitative analysis in the step S5 is comprised the following steps that：

S51：The detection performance of design speed per hour and vehicle checker according to section determines the total area of integration in whole flow rate zone Between；

S52：Away from total interval division will be integrated it is multiple minizones section by certain group, obtains multiple integrating ranges sections, and Using each interval terminal as respective range of integration value；

S53：Data sample to each curve carries out statistical analysis, it is determined that each data volume accounts for sample in the section of each by stages The ratio p of sum_i, and with ratio p_iAs the weight of each segment integrated value；

S54：For the curve of section different times, the multiple segments to above-mentioned division implement integration fortune along curve respectively Calculate, the traffic circulation state characterization value or mutual traffic circulation state in each segment are obtained with respect to characterization value, by following Formula defines two indices to characterize：

(1) traffic circulation state difference characterization value Δ S=S₂-S₁；

(2) traffic circulation state relative different ratio

Wherein, S₁、S₂Each curve corresponding traffic circulation state characterization value in certain flow is interval is represented respectively；Using The quantitative measuring and calculating of implementation is integrated to flow along curve in certain flow is interval, its quantitative measuring and calculating uses below equation：

In formula, Q_b1Represent upper limit of integral value of the contrast period 1 in certain flow interval；Q_a1Represent contrast period 1 one Lower limit of integral value in constant flow interval；A₁Represent the speed-secondary term coefficient of discharge relation formula in contrast period 1；B₁It is right to represent Than the speed-discharge relation formula Monomial coefficient in period 1；C₁Represent the speed-discharge relation formula constant term in contrast period 1； Q_b2Represent upper limit of integral value of the contrast period 2 in certain flow interval；Q_a2Represent contrast period 2 in certain flow is interval Lower limit of integral value；A₂Represent the speed-secondary term coefficient of discharge relation formula in contrast period 2；B₂Represent the speed in contrast period 2 Degree-discharge relation formula Monomial coefficient；C₂Represent the speed-discharge relation formula constant term in contrast period 2；

Here with two curves in contrast period only for explanation this method.

S55：To the curve of section different times, using respective ratio p_iTraffic circulation state to each segment is related Quantitatively characterizing value weighted calculation, obtains section different times to be assessed each speed-corresponding traffic circulation of discharge relation curve State correlation table value indicative or mutual traffic circulation state are with respect to characterization value (S1, S2, Δ S, Δ_O)；Be given in the following manner Quantitative analysis results：

(1) according to integration formula and ratio p_iIt is calculated the corresponding traffic circulation state correlation table value indicative of each curve, institute State traffic circulation state characterization value S of the traffic circulation state correlation table value indicative including different times₁、S₂, traffic circulation state difference Characterization value Δ S and traffic circulation state relative different ratio of Δ_O；Traffic circulation state correlation table value indicative S₁、S₂, Δ S and Δ_O It is calculated as follows：

If 1) integration of each curve is total interval different, that is, there is each by stages section to differ, then not directly ask for Δ S_iWith Δ_Oi(the traffic circulation state difference characterization value in i.e. each by stages, relative different ratio) is, it is necessary to first ask for each subregion of each curve Between section S_1(i)And S_2(i), recycle the S of each by stages section of each curve_1(i)And S_2(i)Ask for Δ S and Δ_O；

In formula：Q_b1(i)、Q_a1(i)、p_1(i)The upper limit of integral value, integration of 1 curve of contrast period, i-th by stages are represented respectively Lower limit, sample data volume accounting；Q_b2(i)、Q_a2(i)、p_2(i)Represent respectively on the integration of 2 curve of contrast period, i-th by stages Limit value, lower limit of integral value, sample data volume accounting；S_1(i)And S_2(i)I-th by stages section of each contrast period curve is represented respectively Traffic circulation state characterization value；n₁Represent that the curve of contrast period 1 integrates the by stages segment number that total interval includes；n₂Represent contrast The curve of period 2 integrates the by stages segment number that total interval includes；

Explanation：Here with two curves in contrast period only for illustration method.

If 2) integration of each curve is total interval identical, that is, there is each by stages section identical, then have：

In formula, Δ S_iRepresent the traffic circulation state difference characterization value in each by stages；Δ_OiRepresent the friendship in each by stages Logical running status relative different ratio；N represents the by stages segment number that the total interval of integration of each curve includes, now each curve By stages hop count is identical；Other specification implication is similar to above-mentioned.

(2) according to Δ S or Δ_OWith the relation of " 0 ", the quality relatively of section different times traffic circulation state is differentiated, and Quantitative analysis result (Δ S or Δ_O)；Method of discrimination is as follows：

If A. Δ S ＞ 0 or Δ_O＞ 0, then it represents that traffic circulation state of the section in period 2 is better than period 1, and lifting values It is Δ S or relative lifting Δs_O；

If B. Δ S ＜ 0 or Δ_O＜ 0, then it represents that traffic circulation state of the section in period 2 is worse than period 1, and negative lifting It is Δ S or relative lifting Δs to be worth_O；

If C. Δ S=0 or Δ_O=0, then it represents that section is substantially suitable with period 1 in the traffic circulation state in period 2.This When, it is necessary to further do qualitative discrimination, if curve substantially completely overlaps, illustrate two traffic circulation states in period of contrast Substantially it is suitable；If curve not exclusively overlaps, and exist intersect and cannot distinguish between any bar curve it is more above when, then pass through The gradual degree of curve judges that curve is gentler, represents that speed reduces more slow, the traffic circulation in the period when flow increases More steady, i.e., traffic circulation state is relatively slightly good.

(3) if having selected multiple track comparative analyses, after completing the quantitative analysis in each track by above-mentioned steps, ask for each Track quantitative analysis results Δ S or Δ_OAverage.

Vehicle checker vehicle flowrate in the step S2 is converted to standard vehicle vehicle flowrate and carries out as follows：

In formula, Q is the standard vehicle vehicle flowrate after conversion；q_iIt is i-th kind of flow value of vehicle；w_iIt is i-th kind of conversion of vehicle Coefficient；M is vehicle quantity.

Rejecting abnormalities data in the step S2 are carried out as follows：

S21：Reject the data of apparent error in the following manner according to speed, the actual effectiveness of flow, it is specific as follows：

Data to following a few class situations are rejected：

Reject speed and be more than 0, but flow is equal to 0 data；

Reject speed and be equal to 0, but data of the flow more than 0；

S22：The rejecting treatment of abnormal data is done using the average and standard deviation of vehicle checker speed data, specific rules are such as Under：

1) if there isIt is then valid data, not rejects；

2) if there isOrThen judge v_iThe corresponding data is abnormal data, should Give rejecting；Wherein, 3 σ are screening criteria, represent that the speed normal distribution three times in the section vehicle checker data for obtaining actual measurement are square Difference；It is the vehicle checker speed data average in the period to be analyzed, standard deviation is σ_v, actual vehicle checker velocity amplitude is v_i。

The total interval of integration in the step S51 determines according to the following steps：

S511：It is determined that the total interval upper limit of integration, to carrying out statistical analysis respectively for being fitted the sample data of each curve, Flow maximum is obtained, as the higher limit of whole flow rate zone；

S512：It is determined that the total interval lower limit of integration, specific as follows：

Note section Maximum speed limit is V0, and speed is the longitudinal axis, at the Maximum speed limit V0 from the longitudinal axis, makees straight line V=V0；Utilize Maximum speed limit come aid in determine lower limit of integral；

When straight line V=V0 intersects with two curves to be analyzed in low flow volume region, then total interval lower limit is integrated It is Maximum speed limit V0；

When tops of the straight line V=V0 in two curves to be analyzed, then lower limit of integral value sets in the following manner：Obtain If wagon flow value is Q0 in vehicle checker detection performance requirement unit detection time, make straight line Q=Q0；Respectively with two curve intersections Point, as the corresponding lower limit of integral of each curve；

When straight line V=V0 intersects with two curves, if the joining on curve is at low discharge, respective friendship is taken The corresponding flow value of point is used as lower limit of integral value；If set lower limit owes reasonable, acquisition sets vehicle checker detection performance It is required that wagon flow value is Q0 in unit detection time, make straight line Q=Q0；Respectively with two curve intersection points, as each curve pair The lower limit of integral answered；

When straight line V=V0 is in low flow volume region and a certain bar curve intersection, then take the flow value and make with the average of Q0 It is lower limit of integral value；

Vehicle checker inspection is set when straight line V=V0 is obtained when low flow volume region is located at below two curves or is not intersected, then It is Q0 to survey wagon flow value in performance requirement unit detection time, makees straight line Q=Q0；Respectively with two curve intersection points, as each The corresponding lower limit of integral of curve；

There is no actual data point at the lower limit of integral value for determining as stated above, then judge whether the lower limit of integral value is big In Q0, if judged result is " being more than ", the flow minimum in pretreated data sample of learning from else's experience is used as lower limit of integral value；

It is Q0 otherwise to obtain and set wagon flow value in vehicle checker detection performance requirement unit detection time, makees straight line Q=Q0；Point Not with two curve intersection points, as the corresponding lower limit of integral of each curve；If lower limit of integral value does not have real data when being Q0 yet Point, then the flow minimum in pretreated data sample of learning from else's experience is used as lower limit of integral value.

It is described that each the total interval of the corresponding integration of speed-discharge relation curve is used whole flow in the period of different The mode of Concourse Division treatment is calculated, specific as follows：

If 1) need not ensure that each by stages Duan Jun is contrasted under equal flow, each contrast period curve product That divides is total interval different, and the integration that the method according to step S51 asks for each curve is always interval.

[Q_a(i),Q_b(i)]；

Wherein, i=1,2 ..., represent that the curve integration of different times is total interval, that is, contrast the curve integrating range in period i.

2) if desired ensure that each by stages Duan Jun is contrasted under equal flow, i.e., each contrast period curve integration Take same integration total interval.The corresponding interval of each curve is obtained, each smaller lower limit in interval is selected as under the total interval of integration Limit, selects larger higher limit as total interval upper limit is integrated, and the total interval of the integration of each curve is taken as into the total interval of same integration, Total interval [the Q of same integration_a,Q_b] establish rules really, it is：

If A. Q_a1≥Q_a2, and have Q_a1≤Q_a2, then taking the total interval of integration is：[Q_a2,Q_b2]；

If B. Q_a1＞ Q_a2, and have Q_b1＞ Q_b2, then taking the total interval of integration is：[Q_a2,Q_b1]；

If C. Q_a1＜ Q_a2, and have Q_b1＜ Q_b2, then taking the total interval of integration is：[Q_a1,Q_b2]；

If D. Q_a1≤Q_a2, and have Q_b1≥Q_b2, then taking the total interval of integration is：[Q_a1,Q_b1]。

Total interval [the Q of integration_a(i),Q_b(i)] (or [Q_a,Q_b]) bound of by stages is determined according to the following steps：

By certain group away from each total interval [Q of curve integration of m_a(i),Q_b(i)] (or [Q_a,Q_b]) it is divided into n small by stages Section, each small by stages section is specific as follows using each interval terminal as respective range of integration value：

1) if each total interval of curve integration is not to take same interval, point integrating range for each dividing is as follows：

[Q_a(i),Q_a(i)+m]；

[Q_a(i)+m,Q_a(i)+2m]；

[Q_a(i)+2m,Q_a(i)+3m]；

……；

[Q_a(i)+(n-2)m,Q_a(i)+(n-1)m]；

[Q_a(i)+(n-1)m,Q_b(i)]。

Wherein, i=1,2 ..., represent the curve i of above-mentioned institute by stages correspondence different times, i.e., above-mentioned institute by stages is right Than each point of integrating range section that period i curve is integrated.

If 2) each curve integrates total interval and takes same interval [Q_a,Q_b], then point integrating range for dividing is as follows：

[Q_a,Q_a+m]；

[Q_a+m,Q_a+2m]；

[Q_a+2m,Q_a+3m]；

……；

[Q_a+(n-2)m,Q_a+(n-1)m]；

[Q_a+(n-1)m,Q_b]。

Embodiment 2

The present embodiment is differed only in embodiment 1：

Present embodiments provide a kind of based on section speed, the road traffic running status appraisal procedure of data on flows, profit With the speed-flow actual relationship curve obtained based on section vehicle checker data (speed, flow) come comprehensive analysis traffic circulation State, from the operational efficiency in the ensemble average effect of traffic circulation state come concentrated expression section period on year-on-year basis.Based on section car Examine the speed-flow actual relationship curve of device data acquisition, it is possible to achieve qualitative to traffic circulation state and quantitative to score Analysis.Analysis result can aid in illustrating the validity of other assessment mode (such as Utilization assessment metrics evaluation) assessment results, also may be used Pipe personnel are transported with relief road traffic to evaluate the running status of road, history management and control measures, be daily management, maintenance meter Draw and planning construction provides reference, and then guidance management decision-making, formulate or optimization management and control measures, promote the horizontal kimonos of operation management The raising of business level.

The present invention solves above-mentioned technical problem by following technological means：

Based on section speed, the road traffic running status appraisal procedure of data on flows, comprise the following steps：

Step 1：Obtain section to be assessed and respectively contrast the vehicle checker data in period, and data standard is stored in data In table.Vehicle checker data select 5 minutes (or 3 minutes) data, and data item must include speed, the flow in each track, and select Data come from one or more larger vehicle checkers of flow on section.

In at home and abroad studying speed flowrate relational model, it is generally recognized that be within 15 minutes and 5 minutes analysis traffic flow system The statistical interval of analysis model is counted, and《HCM HCM2000》Also so recommend.According to correlation statistical analysis Research, using 3 minutes, the statistical interval of 5 minutes and analytical effect of the analytical effect with 15 minutes that obtains is suitable even more excellent, Additionally, the statistical interval for using is shorter, the effective range for obtaining matched curve is bigger, therefore is more prone to use shorter system Meter interval.The actual traffic running status in section can be more characterized using the vehicle checker data with big flow.

Step 2：Vehicle checker data to obtaining are pre-processed, and each track data are processed respectively, while will treatment It is stored in tables of data data standard afterwards.Data prediction mainly includes two large divisions：By the conversion of vehicle checker vehicle flowrate It is standard vehicle vehicle flowrate, and rejecting abnormalities data.

(1) vehicle checker data on flows is converted into standard vehicle wagon flow value

Different types of vehicle is to the influence degree difference of traffic circulation state, it is necessary to it to be all converted to standard wagon flow Amount, to reflect the influence of different automobile types.Different automobile types flow value is generally converted into standard wagon flow using vehicle conversion coefficient Amount.

According to correlation technique handbook such as《HCM》、《Highway capacity handbook》Deng, and combine car test Device equipment detects the division of vehicle, with every 5 minutes for a data cell, to the flow value vehicle conversion coefficient of different automobile types Weighted sum, to obtain the standard vehicle flow value of every 5 minutes.Processing method sees below formula.

In formula, Q is the standard vehicle vehicle flowrate after conversion；q_iIt is i-th kind of flow value of vehicle；w_iIt is i-th kind of conversion of vehicle Coefficient；M is vehicle quantity.

(2) rejecting abnormalities data

In detection process, due to the influence of vehicle checker equipment reason and other factors in itself, may there is part Detection data is abnormal data, and abnormal data can not reflect actual traffic circulation state, and to speed-discharge relation The fitting effect influence of curve is larger, it is necessary to which it is rejected.Herein, the rejecting of abnormal data is completed in two steps：

1) data of apparent error first, are rejected according to speed, the actual effectiveness of flow；

2) rejecting of the average and standard deviation of vehicle checker speed data to do abnormal data is recycled to process.

Step 3：Using by the actual vehicle checker speed in the section after data prediction and data on flows, implement road to be assessed The fitting of section different times speed-discharge relation curve.The data of analysis be typically chosen data on the big track of flow or Treatment Analysis are distinguished to each track data, and according to correlative study, the reality in section is characterized in the form of conic section Speed-discharge relation.

For some sections, its actual speed-discharge relation curve curve form or the equal energy of conic section form Preferably reflect actual traffic circulation state, but according to the speed to being obtained based on section vehicle checker data (speed, flow) The Fitting Analysis of degree-flow actual relationship curve, a curve form is better than using the fitting effect of conic section form, The actual traffic running status of road can preferably be reflected.Additionally,《Highway capacity handbook》Under recommending each calculating speed Traffic flow character curve in the form of conic section (specific discuss refer to the handbook).Therefore, using following Mathematical Modeling To describe speed-discharge relation.

V=AQ²+BQ+C

In formula：V represents speed；Q represents flow；A, B, C are regression coefficient.

Fitting to curve can be realized by MATLAB or other analysis softwares.

Step 4：Determine qualitative analysis or quantitative analysis mode.If only needing to qualitative analysis, can be according to the position of curve Put relation and its respective variation tendency provides analysis result, and analysis terminates.Otherwise, quantitative analysis can directly be implemented, additionally, According to qualitative analysis can not judge well section respectively contrast period traffic circulation state it is fine or not when, then implement quantitative point Analysis, goes to step 5.When the multiple track comparative analyses of selection, track need to be distinguished, the speed in each track of period li section will be contrasted Degree-discharge relation curve is analyzed respectively.

Step 5：Based on the speed-discharge relation curve for obtaining, according to effective data area, it is determined that integration is upper and lower Limit.In quantitative analysis, quantitative analysis results are obtained using integral operation, and whole flow rate zone (integration is total interval) is drawn It is divided into multiple minizone sections to be processed, therefore, primary is to determine integrating range.The step of determining integrating range mainly includes Following components：

(1) first, design speed per hour and the detection performance of vehicle checker according to section determines that a reasonable interval (is entirely flowing Integration in amount interval is always interval)；

(2) secondly, as certain group away from (as with 10pcu/ln/5min as group away from) by previous step obtained by interval division be Multiple minizones section, obtains multiple integrating ranges sections, and using each interval terminal as respective range of integration value.

Because there is upper lower limit value in the actual flow in section in the regular period, the integration fortune in quantitative analysis is implemented During calculation, integrating range can not select [0 ,+∞).Additionally, in whole flow rate zone, it is impossible to which all of data are broadly all used In integration, it is necessary to be selected to concentrate the data point of overall reflection actual traffic running status, also, different flow is whole to section The influence degree of body traffic circulation state is different, and in different flow rate zones, the data volume for characterizing traffic circulation state is different, Also difference, it is necessary to be handled differently, otherwise can produce influence to its extent of reaction to traffic circulation state on result.Therefore, Whole flow rate zone is divided into multiple minizones sections and is handled differently by needs, and then determines rational integrating range, in terms of improving The reliability of calculation.

Step 6：Data sample to each curve carries out statistical analysis, it is determined that each each by stages section in data volume or its Account for the ratio p of total sample number_i(frequency), and with p_iAs the weight of each segment integrated value.

Different flow is different to the influence degree of section entirety traffic circulation state, need to reflect respective otherness.Though Data point so in big flow segment is little, but still reflects actual traffic operation shape of the road in some times State, and the overall quality of road traffic running status by big flow when influenceed also larger, it is necessary to reasonable consideration, but it is unsuitable directly Connect and whole big flow by stages section is integrated.Therefore, accounting for the ratio of total sample number as each area by the use of each interval censored data amount Between section integrated value weight, to reflect the otherness of traffic circulation state under different flow.

Step 7：For the curve of section different times, the multiple segments to above-mentioned division implement integration along curve respectively Computing, obtains the traffic circulation state characterization value or mutual traffic circulation state in each segment with respect to characterization value.And it is fixed Adopted two indices are characterized：

(1) traffic circulation state difference characterization value Δ S=S₂-S₁；

(2) traffic circulation state relative different ratio

Wherein, S1, S2 represent each curve (two, section different times) corresponding traffic in certain flow is interval respectively Running status characterization value.The quantitative measuring and calculating of implementation, its quantitative survey are integrated to flow along curve using in certain flow is interval Calculation method is：

When have selected multiple track comparative analyses, track need to be distinguished, the speed in each track of period li section will be contrasted Degree-discharge relation curve carries out subsection integral treatment respectively.

Step 8：For the curve of section different times, using respective p_iTraffic circulation state to each segment is related Quantitatively characterizing value weighted calculation, obtains section different times to be assessed each speed-corresponding traffic circulation of discharge relation curve State correlation table value indicative or mutual traffic circulation state are with respect to characterization value (S1, S2, Δ S, Δ_O).Then quantitative point is implemented Analysis, and quantitative analysis result.When have selected multiple track comparative analyses, track need to be distinguished, period li will be contrasted Speed-discharge relation the curve in each track of section is analyzed respectively, finally can further provide each track quantitative result Average come characterize road section traffic volume operation integral status.Implementation steps mainly include following components：

(1) according to integration formula and p_iIt is calculated the corresponding traffic circulation state correlation table value indicative of each curve：When different The traffic circulation state characterization value of phase is (such as：S₁、S₂), traffic circulation state difference characterization value Δ S and traffic circulation state it is relative Difference ratio of Δ_O。

(2) according to Δ S or Δ_OWith the relation of " 0 ", the quality relatively of section different times traffic circulation state is differentiated, and Quantitative analysis result.

(3) if having selected multiple track comparative analyses, after completing the quantitative analysis in each track by above-mentioned steps, one can be entered Step asks for each track quantitative analysis results Δ S or Δ_OAverage.

The present invention is simple and easy to apply, it is readily appreciated that, and in road traffic service data, vehicle checker data are relatively rich Richness, the method is easily realized under actual conditions.Based on speed-flow reality that section vehicle checker data (speed, flow) obtains Relation curve carrys out the traffic circulation state of comparative analysis same section different times, can apply to commenting for road running status Estimate in analysis, analysis result can aid in fortune pipe portion door accurately to hold road traffic running status on the whole, also can be auxiliary Help the validity for illustrating other assessment mode (such as Utilization assessment metrics evaluation) assessment results.

Embodiment 3

The present embodiment is differed only in embodiment 1：

Referring to Fig. 1, based on section speed, the road traffic running status appraisal procedure of data on flows, comprise the following steps：

Step 1：Obtain section to be assessed and respectively contrast the vehicle checker data in period, and data standard is stored in data In table.Vehicle checker data select 5 minutes (or 3 minutes) data, and data item must include speed, the flow in each track, and select Data come from one or more larger vehicle checkers of flow on section.

In at home and abroad studying speed flowrate relational model, it is generally recognized that be within 15 and 5 minutes analysis traffic flow statistics point The statistical interval of model is analysed, and《HCM HCM2000》Also so recommend.Ground according to correlation statistical analysis Study carefully, using 3 minutes, the statistical interval of 5 minutes and analytical effect of the analytical effect with 15 minutes that obtains is suitable even more excellent, this Outward, the statistical interval of use is shorter, and the effective range for obtaining matched curve is bigger, therefore is more prone to use shorter statistics Interval.The actual traffic running status in section can be more characterized using the vehicle checker data with big flow.

Step 2：Vehicle checker data to obtaining are pre-processed, and each track data are processed respectively, while will treatment It is stored in tables of data data standard afterwards.Data prediction mainly includes two large divisions：By the conversion of vehicle checker vehicle flowrate It is standard vehicle vehicle flowrate, and rejecting abnormalities data.

(1) vehicle checker data on flows is converted into standard vehicle wagon flow value

Different types of vehicle is to the influence degree difference of traffic circulation state, it is necessary to it to be all converted to standard wagon flow Amount, to reflect the influence of different automobile types.Different automobile types flow value is generally converted into standard wagon flow using vehicle conversion coefficient Amount.

According to correlation technique handbook such as《HCM》、《Highway capacity handbook》Deng, and combine car test Device equipment detects the division of vehicle, with every 5 minutes for a data cell, to the flow value vehicle conversion coefficient of different automobile types Weighted sum, to obtain the standard vehicle flow value of every 5 minutes.Processing method sees below formula.

In formula, Q is the standard vehicle vehicle flowrate after conversion；q_iIt is i-th kind of flow value of vehicle；w_iIt is i-th kind of conversion of vehicle Coefficient；M is vehicle quantity.

(2) rejecting abnormalities data

In detection process, due to the influence of vehicle checker equipment reason and other factors in itself, may there is part Detection data is abnormal data, and abnormal data can not reflect actual traffic circulation state, and to speed-discharge relation The fitting effect influence of curve is larger, it is necessary to which it is rejected.Herein, the rejecting of abnormal data is completed in two steps, first, The data of apparent error are rejected according to speed, the actual effectiveness of flow, the average and standard of vehicle checker speed data is recycled Rejecting of the difference to do abnormal data is processed.

First, the data to following a few class situations are rejected：

1) speed is more than 0, but flow is equal to 0；

2) speed is equal to 0, but flow is more than 0.

Secondly, further data processing is done.Studied according to correlation statistical analysis, the speed in the section vehicle checker data of actual measurement Angle value Normal Distribution, and by the knowledge of probability theory, the sample of normal population there is deviation more than twice mean square deviation (2 σ) and The possibility of three times mean square deviation (3 σ) is respectively smaller than 5% and 0.3%, is a small probability event.If the deviation of outlier is big In 2 σ or 3 σ, then have reason it to be judged into exceptional value and is rejected.Therefore, selecting 3 σ as screening criteria, and set the period to be analyzed Interior vehicle checker speed data average isStandard deviation is σ_v, actual vehicle checker velocity amplitude is v_i, then abnormal data rejecting rule It is then：

1) if there isThen think that the data is valid data, not reject；

2) if there isOrThen judge v_iThe corresponding data is abnormal data, should Rejected.

Step 3：Using by the actual vehicle checker speed in the section after data prediction and data on flows, implement road to be assessed The fitting of section different times speed-discharge relation curve.The data of analysis be typically chosen data on the big track of flow or Treatment Analysis are distinguished to each track data, and according to correlative study, the reality in section is characterized in the form of conic section Speed-discharge relation.

For some sections, its actual speed-discharge relation curve curve form or the equal energy of conic section form Preferably reflect actual traffic circulation state, but according to the Fitting Analysis to section actual speed-discharge relation curve, Curve form is better than using the fitting effect of conic section form, can preferably reflect the actual traffic operation of road State.Additionally,《Highway capacity handbook》Recommend the traffic flow character curve under each calculating speed using the shape of conic section Formula (specific discussion refers to the handbook).Therefore, describing speed-discharge relation with following Mathematical Modeling.

V=AQ²+BQ+C

In formula：V represents speed；Q represents flow；A, B, C are regression coefficient.

Fitting to curve can be realized by MATLAB or other analysis softwares.Based on section vehicle checker data (speed, Flow) the section actual speed-discharge relation curve of fitting illustrates referring to Fig. 2-Chongqing highway Yu Wulu (descending) two The speed in individual contrast period-discharge relation curve.

Step 4：Determine qualitative analysis or quantitative analysis mode.If only needing to qualitative analysis, can be according to the position of curve Put relation and its respective variation tendency provides analysis result, and analysis terminates.Otherwise, quantitative analysis can directly be implemented.It is elected When selecting multiple track comparative analyses, track need to be distinguished, the speed-discharge relation curve in each track of period li section will be contrasted It is analyzed respectively.

If qualitative analysis is implemented in current selection, analysis method is as follows：

Overlying relation according to each curve judges to treat the excellent of comparative evaluation section contrast period traffic circulation state Bad, judgment rule is as follows：

If 1) speed-discharge relation curve for participating in contrast differs greatly, hence it is evident that there is curve major part region position In another top of curve, then it is considered that the traffic circulation state that curve above is characterized is better than on the whole The traffic circulation state that lower curve is characterized.Referring to Fig. 2, the curve in period 2 is wholly or largely on the curve in period 1 Side, can qualitatively think that running status of the section in period 2 is better than period 1 on the whole；

If 2) two curves substantially completely overlap, the basic phase of traffic circulation state in contrast period can be qualitatively thought When；

If 3) curve not exclusively overlaps, and exist intersect and cannot distinguish between any bar curve it is more above when, if certain Curve is gentler, represents that speed reduction is more slow when flow increases, and the traffic circulation in the period is more steady more preferable；

4) can not well judge that section respectively contrasts the fine or not or current of period traffic circulation state according to qualitative analysis Without only doing qualitative analysis, then implement quantitative analysis, go to step 5.

Step 5：Based on the speed-discharge relation curve for obtaining, according to effective data area, it is determined that integration is upper and lower Limit.In quantitative analysis, analysis result is obtained using integral operation, and whole flow rate zone (integration is total interval) is divided into Multiple minizone sections are processed, therefore, primary is to determine integrating range.First, the design speed per hour according to section and car test The detection performance of device determines a reasonable interval (integration in whole flow rate zone is always interval)；Secondly, by certain group away from (as with 10pcu/ln/5min as group away from) by previous step obtained by interval division be multiple minizones section, obtain multiple integrated areas Between section, and using each interval terminal as respective range of integration value.

Because there is upper lower limit value in the actual flow in section in the regular period, the integration fortune in quantitative analysis is implemented During calculation, integrating range can not select [0 ,+∞).Additionally, in whole flow rate zone, it is impossible to which all of data are broadly all used In integration, it is necessary to be selected to concentrate the data point of overall reflection actual traffic running status, also, different flow is whole to section The influence degree of body traffic circulation state is different, and in different flow rate zones, the data volume for characterizing traffic circulation state is different, Also difference, it is necessary to be handled differently, otherwise can produce influence to its extent of reaction to traffic circulation state on result.Therefore, Whole flow rate zone is divided into multiple minizones sections and is handled differently by needs, and then determines rational integrating range, in terms of improving The reliability of calculation.

Referring to Fig. 3, if two, certain section contrasts period corresponding curve as illustrated, wherein, dash area is table in figure Show the difference between the contrast period traffic circulation state of two, certain section.Meanwhile, speed is represented with V, represent flow with Q.

If section is A in 1 speed of period-corresponding regression coefficient of discharge relation curve₁、B₁、C₁, the corresponding mathematics of curve Expression formula is：

V=A₁Q²+B₁Q+C₁

Section is A in 2 speed of period-corresponding regression coefficient of discharge relation curve₂、B₂、C₂, the corresponding mathematical table of curve It is up to formula：

V=A₂Q²+B₂Q+C₂

The step of determining range of integration is specific as follows：

(1) design speed per hour and the detection performance of vehicle checker according to section determines rational data area, to determine integration Whole reasonable interval.

A. the determination of the total upper limit of integrating range.To carrying out statistical analysis respectively for being fitted the sample data of each curve, obtain Flow maximum is taken, in this, as the higher limit of whole flow rate zone, such as Qb in figure, that is, makees straight line Q=Qb, respectively with two Bar curve intersection, with the corresponding reasonable interval higher limit of each curve integration of this determination.

B. the determination of the total lower limit of integrating range.It is V0 that speed per hour (or Maximum speed limit) is designed in note section, from the longitudinal axis (speed axle) On V0 make straight line V=V0.Aided in using speed per hour (or Maximum speed limit) is designed determine lower limit of integral, when straight line V=V0 with treat Two curves of analysis can be caused in integrating range comprising most valid data when low flow volume region intersects.

If two tops of curve, now, lower limit of integral value can be set to " 0 " to the straight line, according to actual conditions point It is to owe rational to analyse this, it is thus impossible to directly delimit lower limit of integral with this.If vehicle checker detects performance requirement unit detection time Interior (such as 5 minutes), when vehicle flowrate is less than certain value, are designated as Q0, and the degree of accuracy of detection is relatively low, now, are delimited using the value and integrated Qa=Q0 in lower limit, such as figure.Make straight line Q=Qa, respectively with two curve intersections, the corresponding product of each curve is determined with straight line Qa Divide lower limit；

When straight line V=V0 intersects with two curves, if the joining on curve is at low discharge, i.e., in reasonable Scope (the corresponding flow value of intersection point is located near Q0), can take the flow value corresponding to respective intersection point as lower limit of integral Value, it is also possible to take both 1/2nd location points.If set lower limit owes reasonable, weighed in a manner described in conjunction with Q0 New settings lower limit of integral.

When straight line V=V0 is in low flow volume region and a certain bar curve intersection, then take the flow value and make with the average of Q0 It is lower limit of integral value, to ensure to have actual data point at lower limit of integral as far as possible.

When straight line V=V0 is when low flow volume region is located at below two curves or is not intersected, can be drawn directly in conjunction with Q0 Lower limit of the definite integral value, is taken as Qa=Q0.

If additionally, there is no actual data point at the lower limit of integral value for determining as stated above, judging the lower limit of integral Whether value is more than Q0, if judged result is " being more than ", the flow minimum in the pretreated data sample that can directly learn from else's experience is made It is lower limit of integral value, otherwise takes Q0.If but lower limit of integral value when being Q0 also without actual data point, also learn from else's experience pretreated Flow minimum in data sample is used as lower limit of integral value.

It is derived from the corresponding integration of each curve total interval.If two periods are contrasted, respective speed-discharge relation The corresponding total interval of integration of curve is respectively：[Q_a1,Q_b1] and [Q_a2,Q_b2].If being contrasted to multiple periods, each period curve product Divide total interval for [Q_a(i),Q_b(i)], wherein, i=1,2 ... represent that the curve integration of different times is total interval, that is, contrast period i Curve integrating range.

Due in quantitative calculating, being calculated by the way of by whole flow rate zone segment processing, if desired ensure Each by stages Duan Jun is contrasted under equal flow, can be after the corresponding reasonable interval of each curve is obtained, and selection is each Interval smaller lower limit selects larger higher limit as the total interval upper limit of integration as the total interval lower limit of integration, will be each The total interval of integration of curve is taken as same, so may be such that the segment that each curve is divided is identical.But the part for so obtaining Or even the total interval of integration of each curve can be interval slightly larger than actual, there is certain error, because the interval being extended is distributed in The two ends (being predominantly located at big flow area) of curve, data volume is relatively little, the error for existing also very little, and it is also reasonable to be processed as 's.Therefore, the same total interval [Q of integration_a,Q_b] establish rules really, it is：

If A. Q_a1≥Q_a2, and have Q_a1≤Q_a2, then taking the total interval of integration is：[Q_a2,Q_b2]；

If B. Q_a1＞ Q_a2, and have Q_b1＞ Q_b2, then taking the total interval of integration is：[Q_a2,Q_b1]；

If C. Q_a1＜ Q_a2, and have Q_b1＜ Q_b2, then taking the total interval of integration is：[Q_a1,Q_b2]；

If D. Q_a1≤Q_a2, and have Q_b1≥Q_b2, then taking the total interval of integration is：[Q_a1,Q_b1]。

(2) by certain group away from m (such as with m=10pcu/ln/5min as group away from) by the interval [Q obtained by previous step_a(i), Q_b(i)] (or [Q_a,Q_b]) be divided into multiple minizones section n, obtain multiple integrating ranges sections, and using each interval terminal as Respective range of integration value.The range of integration of thus obtained each by stages is respectively：

1) if each total interval of curve integration is not to take same interval, point integrating range for each dividing is as follows：

[Q_a(i),Q_a(i)+m]；

[Q_a(i)+m,Q_a(i)+2m]；

[Q_a(i)+2m,Q_a(i)+3m]；

……；

[Q_a(i)+(n-2)m,Q_a(i)+(n-1)m]；

[Q_a(i)+(n-1)m,Q_b(i)]。

Wherein, i=1,2 ..., represent the curve i of above-mentioned institute by stages correspondence different times, i.e., above-mentioned institute by stages is right Than each point of integrating range section that period i curve is integrated.

If 2) each curve integrates total interval and takes same interval [Q_a,Q_b], will not distinguish, then point integrating range for dividing It is as follows：

[Q_a,Q_a+m]；

[Q_a+m,Q_a+2m]；

[Q_a+2m,Q_a+3m]；

……；

[Q_a+(n-2)m,Q_a+(n-1)m]；

[Q_a+(n-1)m,Q_b]。

Step 6：Data sample to each curve carries out statistical analysis, it is determined that each each by stages section in data volume or its Account for the ratio p of total sample number_i(frequency), and with p_iAs the weight of each segment integrated value.

Different flow is different to the influence degree of section entirety traffic circulation state, need to reflect respective otherness.Though Data point so in big flow segment is little, but still reflects actual traffic operation shape of the road in some times State, and the overall quality of road traffic running status by big flow when influenceed also larger, it is necessary to reasonable consideration, but it is unsuitable directly Connect and whole big flow by stages section is integrated.Therefore, accounting for the ratio of total sample number as each area by the use of each interval censored data amount Between section integrated value weight, to reflect the otherness of traffic circulation state under different flow.

Step 7：For the curve of section different times, the multiple segments to above-mentioned division implement integration along curve respectively Computing, obtains the traffic circulation state characterization value or mutual traffic circulation state in each segment with respect to characterization value.And it is fixed Adopted two indices are characterized：Traffic circulation state difference characterization value Δ S, traffic circulation state relative different ratio of Δ_O；Work as selection During multiple track comparative analyses, track need to be distinguished, the speed-discharge relation curve in each track of period li section will be contrasted Subsection integral treatment is carried out respectively.

(1) quantitative target, and quantitative measuring method are defined

A. traffic circulation state difference characterization value Δ S.The traffic circulation state difference between different comparative analysis objects is represented, Quantitatively characterizing is come by the difference of respective characterization value, i.e., quantitatively measure formula is：Δ S=S₂-S₁；

B. traffic circulation state relative different ratio of Δ_O.Represent the phase of traffic circulation state between different comparative analysis objects To difference, quantitatively characterizing is come by the difference and the percentage of contrast period characterization value of respective characterization value, i.e., quantitatively measure formula is：

Each curve corresponding traffic circulation state in certain flow interval is calculated according to such as lower integral formula to characterize Value S1, S2 and Δ S and Δ_O。

If the range of integration corresponding to each curve is identical, that is, sets each integrating range and be [Qa, Qb], then have：

Above formula is to carry out averagely, characterizing section in whole flow rate zone with respect to lifting values by speed under equal flow The average relative mistake opposite sex of traffic circulation state between different times.

(2) quantitative measuring and calculating is implemented to each segment

If 1) integration of each curve is total interval different, that is, there is each by stages section to differ, then not directly ask for Δ S_iWith Δ_Oi, it is necessary to first ask for the S of each by stages section of each curve_1(i)And S_2(i)。

If 2) integration of each curve is total interval identical, that is, there is each by stages section identical, then have：

Step 8：For the curve of section different times, using respective p_iTraffic circulation state to each segment is related Quantitatively characterizing value weighted calculation, obtains section different times to be assessed each speed-corresponding traffic circulation of discharge relation curve State correlation table value indicative or mutual traffic circulation state are with respect to characterization value (S1, S2, Δ S, Δ_O).Then quantitative point is implemented Analysis, and quantitative analysis result.When have selected multiple track comparative analyses, track need to be distinguished, period li will be contrasted Speed-discharge relation the curve in each track of section is analyzed respectively, finally can further provide each track quantitative result Average come characterize road section traffic volume operation integral status.Implementation steps mainly include following components：

(1) p is utilized_iTo the traffic circulation state characterization value weighted calculation of each segment, to obtain traffic circulation state pair Than the correlation table value indicative analyzed.

If 1) integration of each curve is total interval different, that is, there is each by stages section to differ, due to not asking for Δ S directly_iWith Δ_Oi, it is necessary to using the S of each by stages section of each curve_1(i)And S_2(i)Ask for.

If 2) integration of each curve is total interval identical, that is, there is each by stages section identical, then have：

(2) according to Δ S or Δ_ODifferentiate and quantitative analysis result.

If A. Δ S ＞ 0 or Δ_O＞ 0, then it represents that traffic circulation state of the section in period 2 is better than period 1, and lifting values It is Δ S or relative lifting Δs_O；

If B. Δ S ＜ 0 or Δ_O＜ 0, then it represents that traffic circulation state of the section in period 2 is worse than period 1, and negative lifting It is Δ S or relative lifting Δs to be worth_O；

If C. Δ S=0 or Δ_O=0, then it represents that section is substantially suitable with period 1 in the traffic circulation state in period 2.This When, it is necessary to further do qualitative discrimination, if curve substantially completely overlaps, illustrate two traffic circulation states in period of contrast Substantially it is suitable；If curve not exclusively overlaps, and exist intersect and cannot distinguish between any bar curve it is more above when, then pass through The gradual degree of curve judges that curve is gentler, represents that speed reduces more slow, the traffic circulation in the period when flow increases More steady, i.e., traffic circulation state is relatively slightly good.

(3) if having selected multiple track comparative analyses, after completing the quantitative analysis in each track by above-mentioned steps, one can be entered Step asks for each track quantitative analysis results Δ S or Δ_OAverage.

Finally illustrate, the above embodiments are merely illustrative of the technical solutions of the present invention and it is unrestricted, although by ginseng According to the preferred embodiments of the present invention, invention has been described, it should be appreciated by those of ordinary skill in the art that can To make various changes to it in the form and details, without departing from the spirit and scope that the present invention is limited.

Claims (6)

1. section speed, the road traffic running status appraisal procedure of data on flows are based on, it is characterised in that：Including following step Suddenly：

S1：Obtain section to be assessed and respectively contrast the vehicle checker data in period, and data standard is stored in tables of data, institute Stating the data item in tables of data includes speed, the flow in each track；

S2：Vehicle checker data to obtaining carry out data prediction, and each track data are processed respectively, while by after treatment It is stored in tables of data data standard；The data prediction includes two large divisions：Vehicle checker vehicle flowrate is converted into mark Quasi- car vehicle flowrate, and rejecting abnormalities data；

S3：Using by the actual vehicle checker speed in the section after data prediction and data on flows, implement section to be assessed different The fitting of period speed-discharge relation curve；Actual speed-the flow in section is fitted in the form of following conic section Relation：

V=AQ²+BQ+C

In formula：V represents speed；Q represents flow；A, B, C are regression coefficient；

S4：It is determined that to the analysis mode of road traffic running status assessment, if carrying out qualitative analysis, according to the position of curve Relation and its respective variation tendency provide the qualitative analysis, and terminate road traffic running status evaluation process；If entered Row quantitative analysis, goes to step S5 and implements quantitative analysis；

S5：The quantitative analysis utilizes the speed-discharge relation curve for obtaining, according to effective data area, it is determined that integrate Bound；And it is that multiple minizone sections are processed that whole flow is integrated into total interval division, is determined using integral operation Amount analysis result；

Quantitative analysis in the step S5 is comprised the following steps that：

S51：The detection performance of design speed per hour and vehicle checker according to section determines that the integration in whole flow rate zone is always interval；

S52：Away from total interval division will be integrated it is multiple minizones section by certain group, obtains multiple integrating ranges sections, and with each Interval terminal is used as respective range of integration value；

S53：Data sample to each curve carries out statistical analysis, it is determined that each data volume accounts for total sample number in the section of each by stages Ratio p_i′, and with ratio p_i′As the weight of each segment integrated value；Ratio p_i′Calculated as follows：

$p_{i^{\′}}=\frac{N_{i^{\′}}}{N_0};$

In formula, N₀It is certain contrast period speed-flow curve total sample number of section；N_i′It is certain contrast period curve of section Sample data volume in each by stages section；

S54：For the curve of section different times, the multiple segments to above-mentioned division implement integral operation respectively along curve, Traffic circulation state characterization value or mutual traffic circulation state in each segment are obtained with respect to characterization value, as follows Two indices are defined to characterize：

(1) traffic circulation state difference characterization value Δ S=S₂-S₁；

(2) traffic circulation state relative different ratio

Wherein, S₁、S₂Each curve corresponding traffic circulation state characterization value in certain flow is interval is represented respectively；Using one The quantitative measuring and calculating of implementation is integrated to flow along curve in constant flow is interval, its quantitative measuring and calculating uses below equation：

$S_1={\∫}_{Q_{a1}}^{Q_{b1}}(A_1{Q}^2+B_{1}Q+C_1)dQ;$

$S_2={\∫}_{Q_{a2}}^{Q_{b2}}(A_2{Q}^2+B_{2}Q+C_2)dQ;$

In formula, Q_b1Represent upper limit of integral value of the contrast period 1 in certain flow interval；Q_a1Represent contrast period 1 in a constant current Lower limit of integral value in amount interval；A₁Represent the speed-secondary term coefficient of discharge relation formula in contrast period 1；B₁When representing contrast The speed of phase 1-discharge relation formula Monomial coefficient；C₁Represent the speed-discharge relation formula constant term in contrast period 1；Q_b2Table Show upper limit of integral value of the contrast period 2 in certain flow interval；Q_a2Represent integration of the contrast period 2 in certain flow interval Lower limit；A₂Represent the speed-secondary term coefficient of discharge relation formula in contrast period 2；B₂Represent the speed-flow in contrast period 2 Relational expression Monomial coefficient；C₂Represent the speed-discharge relation formula constant term in contrast period 2；

S55：To the curve of section different times, using respective ratio p_i′Traffic circulation state to each segment is related quantitative Characterization value weighted calculation, obtains section different times to be assessed each speed-corresponding traffic circulation state of discharge relation curve Correlation table value indicative or mutual traffic circulation state are with respect to characterization value S1, S2, Δ S, Δ_O；Quantitative divides in the following manner Analysis result：

(1) according to integration formula and ratio p_i′It is calculated the corresponding traffic circulation state correlation table value indicative of each curve, the friendship Logical running status correlation table value indicative includes the traffic circulation state characterization value S of different times₁、S₂, traffic circulation state difference characterize Value Δ S and traffic circulation state relative different ratio of Δ_O；Traffic circulation state correlation table value indicative S₁、S₂, Δ S and Δ_OBy such as Under type is calculated：

1) if the integration of each curve is total interval different, that is, there is each by stages to differ, then not directly ask for Δ S_iAnd Δ_Oi, i.e., Traffic circulation state difference characterization value, traffic circulation state relative different ratio in each by stages is, it is necessary to first ask for each curve The S of each by stages section_1(i)And S_2(i), recycle the S of each by stages section of each curve_1(i)And S_2(i)Ask for Δ S and Δ_O；

$S_{1\left(i\right)}={\∫}_{Q_{a1\left(i\right)}}^{Q_{b1\left(i\right)}}(A_1{Q}^2+B_{1}Q+C_1)dQ$

$S_{2\left(i\right)}={\∫}_{Q_{a2\left(i\right)}}^{Q_{b2\left(i\right)}}(A_2{Q}^2+B_{2}Q+C_2)dQ$

$\mathrm{\Δ}S=S_2-S_1=\underset{i=1}{\overset{n_2}{\Σ}}(p_{2\left(i\right)}*S_{2\left(i\right)})-\underset{i=1}{\overset{n_1}{\Σ}}(p_{1\left(i\right)}*S_{1\left(i\right)})$

${\mathrm{\Δ}}_O=\frac{S_2-S_1}{S_1}=\frac{\underset{i=1}{\overset{n_2}{\Σ}}(p_{2\left(i\right)}*S_{2\left(i\right)})-\underset{i=1}{\overset{n_1}{\Σ}}(p_{1\left(i\right)}*S_{1\left(i\right)})}{\underset{i=1}{\overset{n_1}{\Σ}}(p_{1\left(i\right)}*S_{1\left(i\right)})}$

In formula：Q_b1(i)、Q_a1(i)、p_1(i)The upper limit of integral value, lower limit of integral of 1 curve of contrast period, i-th by stages are represented respectively Value, sample data volume accounting；Q_b2(i)、Q_a2(i)、p_2(i)Represent respectively 2 curve of contrast period i-th by stages upper limit of integral value, Lower limit of integral value, sample data volume accounting；S_1(i)And S_2(i)The traffic fortune of each i-th by stages of contrast period curve is represented respectively Row state table value indicative；n₁Represent that the curve of contrast period 1 integrates the by stages segment number that total interval includes；n₂Represent that contrast period 2 is bent The by stages segment number that the total interval of line integral includes；

If 2) integration of each curve is total interval identical, that is, there is each by stages section identical, then have：

$S_{1\left(i\right)}={\∫}_{Q_{a1\left(i\right)}}^{Q_{b1\left(i\right)}}(A_1{Q}^2+B_{1}Q+C_1)dQ=p_{1\left(i\right)}{\∫}_{Q_{a\left(i^{\′}\right)}}^{Q_{b\left(i^{\′}\right)}}(A_1{Q}^2+B_{1}Q+C_1)dQ$

$S_{2\left(i\right)}={\∫}_{Q_{a2\left(i\right)}}^{Q_{b2\left(i\right)}}(A_2{Q}^2+B_{2}Q+C_2)dQ=p_{2\left(i\right)}{\∫}_{Q_{a\left(i^{\′}\right)}}^{Q_{b\left(i^{\′}\right)}}(A_2{Q}^2+B_{2}Q+C_2)dQ$

$\begin{array}{c}{\mathrm{\ΔS}}_i=S_{2\left(i\right)}-S_{1\left(i\right)}=p_{2\left(i\right)}*{\∫}_{Q_{a\left(i^{\′}\right)}}^{Q_{b\left(i^{\′}\right)}}(A_2{Q}^2+B_{2}Q+C_2)dQ-p_{1\left(i\right)}*{\∫}_{Q_{a\left(i^{\′}\right)}}^{Q_{b\left(i^{\′}\right)}}(A_1{Q}^2+B_{1}Q+C_1)dQ\\ =p_{2\left(i\right)}*{\∫}_{Q_{a\left(i^{\′}\right)}}^{Q_{b\left(i^{\′}\right)}}\[\left(A_2{Q}^2+B_{2}Q+C_2\right)-\left(\frac{p_{1\left(i\right)}}{p_{2\left(i\right)}}\right)*\left(A_1{Q}^2+B_{1}Q+C_1\right)\]dQ\end{array}$

$\begin{array}{c}{\mathrm{\Δ}}_{Oi}=\frac{S_{2\left(i\right)}-S_{1\left(i\right)}}{S_{1\left(i\right)}}\×100\%\\ =\frac{p_{2\left(i\right)}*{\∫}_{Q_{a\left(i^{\′}\right)}}^{Q_{b\left(i^{\′}\right)}}\left(A_2{Q}^2+B_{2}Q+C_2\right)dQ-p_{1\left(i\right)}*{\∫}_{Q_{a\left(i^{\′}\right)}}^{Q_{b\left(i^{\′}\right)}}\left(A_1{Q}^2+B_{1}Q+C_1\right)dQ}{p_{1\left(i\right)}*{\∫}_{Q_{a\left(i^{\′}\right)}}^{Q_{b\left(i^{\′}\right)}}\left(A_1{Q}^2+B_{1}Q+C_1\right)dQ}\×100\%\\ =\left(\frac{p_{2\left(i\right)}*{\∫}_{Q_{a\left(i^{\′}\right)}}^{Q_{b\left(i^{\′}\right)}}\left(A_2{Q}^2+B_{2}Q+C_2\right)dQ}{p_{1\left(i\right)}*{\∫}_{Q_{a\left(i^{\′}\right)}}^{Q_{b\left(i^{\′}\right)}}\left(A_1{Q}^2+B_{1}Q+C_1\right)dQ}-1\right)\×100\%\\ =\[\left(\frac{p_{2\left(i\right)}}{p_{1\left(i\right)}}\right)*\left(\frac{{\∫}_{Q_{a\left(i^{\′}\right)}}^{Q_{b\left(i^{\′}\right)}}\left(A_2{Q}^2+B_{2}Q+C_2\right)dQ}{{\∫}_{Q_{a\left(i^{\′}\right)}}^{Q_{b\left(i^{\′}\right)}}\left(A_1{Q}^2+B_{1}Q+C_1\right)dQ}\right)-1\]\×100\%\end{array}$

$\mathrm{\Δ}S=\underset{i=1}{\overset{n}{\Σ}}\left({\mathrm{\ΔS}}_i\right)$

${\mathrm{\Δ}}_O=\underset{i=1}{\overset{n}{\Σ}}\left({\mathrm{\Δ}}_{Oi}\right)$

In formula, Δ S_iRepresent the traffic circulation state difference characterization value in each by stages；Δ_OiRepresent the traffic fortune in each by stages Row state relative different ratio；N represents the by stages segment number that the total interval of integration of each curve includes, now the subregion of each curve Between hop count it is identical；

[Q_a(i′),Q_b(i′)] represent that the curve integration of contrast period i ' is total interval, wherein, i '=1,2 ...；

(2) according to Δ S or Δ_OWith the relation of " 0 ", the quality relatively of section different times traffic circulation state is differentiated, and provide fixed Amount analysis result Δ S or Δ_O；Method of discrimination is as follows：

If A. Δ S ＞ 0 or Δ_O＞ 0, then it represents that traffic circulation state of the section in period 2 is better than period 1, and lifting values are Δ S Or relative lifting Δ_O；

If B. Δ S ＜ 0 or Δ_O＜ 0, then it represents that traffic circulation state of the section in period 2 is worse than period 1, and negative lifting values are Δ S or relative lifting Δs_O；

If C. Δ S=0 or Δ_O=0, then it represents that section is substantially suitable with period 1 in the traffic circulation state in period 2；Now, need Qualitative discrimination is further done, if curve substantially completely overlaps, illustrates that the traffic circulation state in two periods of contrast is basic Quite；If curve not exclusively overlaps, and exist intersect and cannot distinguish between any bar curve it is more above when, then by curve Gradual degree judge that curve is gentler, represent that speed reduces more slow when flow increases, the traffic circulation in period 2 is more flat Surely, i.e., traffic circulation state is relatively slightly good；

(3) if having selected multiple track comparative analyses, after completing the quantitative analysis in each track by above-mentioned steps, each track is asked for Quantitative analysis results Δ S or Δ_OAverage.

2. according to claim 1 based on section speed, the road traffic running status appraisal procedure of data on flows, it is special Levy and be：Vehicle checker vehicle flowrate in the step S2 is converted to standard vehicle vehicle flowrate and carries out as follows：

$Q=\underset{i^{\′\′}=1}{\overset{m}{\Σ}}(w_{i^{\′\′}}\·q_{i^{\′\′}})$

In formula, Q is the standard vehicle vehicle flowrate after conversion；q_i″" flow value of vehicle is planted for i-th；w_i″" conversion of vehicle is planted for i-th Coefficient；M is vehicle quantity.

3. according to claim 1 based on section speed, the road traffic running status appraisal procedure of data on flows, it is special Levy and be：Rejecting abnormalities data in the step S2 are carried out as follows：

S21：Reject the data of apparent error in the following manner according to speed, the actual effectiveness of flow, it is specific as follows：

Data to following a few class situations are rejected：

Reject speed and be more than 0, but flow is equal to 0 data；

Reject speed and be equal to 0, but data of the flow more than 0；

S22：The rejecting treatment of abnormal data is done using the average and standard deviation of vehicle checker speed data, specific rules are as follows：

1) if there isIt is then valid data, not rejects；

2) if there isOrThen judge v_IThe corresponding data is abnormal data, should be picked Remove；Wherein, 3 σ are screening criteria, represent the speed normal distribution three times mean square deviation in the section vehicle checker data for obtaining actual measurement； It is the vehicle checker speed data average in the period to be analyzed, standard deviation is σ_v, actual vehicle checker velocity amplitude is v_I。

4. according to claim 1 based on section speed, the road traffic running status appraisal procedure of data on flows, it is special Levy and be：The total interval of integration in the step S51 determines according to the following steps：

S511：It is determined that the total interval upper limit of integration, to carrying out statistical analysis respectively for being fitted the sample data of each curve, obtains Flow maximum, as the higher limit of whole flow rate zone；

S512：It is determined that the total interval lower limit of integration, specific as follows：

Note section Maximum speed limit is V0, and speed is the longitudinal axis, at the Maximum speed limit V0 from the longitudinal axis, makees straight line V=V0；Using highest Speed limit come aid in determine lower limit of integral；

When straight line V=V0 intersects with two curves to be analyzed in low flow volume region, then total interval lower limit is integrated for most Speed limit V0 high；

When tops of the straight line V=V0 in two curves to be analyzed, then lower limit of integral value sets in the following manner：Acquisition sets car Wagon flow value is Q0 in inspection device detection performance requirement unit detection time, makees straight line Q=Q0；Respectively with two curve intersection points, As the corresponding lower limit of integral of each curve；

When straight line V=V0 intersects with two curves, if the joining on curve is at low discharge, respective intersection point institute is taken Corresponding flow value is used as lower limit of integral value；If set lower limit owes reasonable, acquisition sets vehicle checker detection performance requirement Wagon flow value is Q0 in unit detection time, makees straight line Q=Q0；It is corresponding as each curve respectively with two curve intersection points Lower limit of integral；

When straight line V=V0 is in low flow volume region and a certain bar curve intersection, then take equal with Q0 with curve intersection point flow value Value is used as lower limit of integral value；

Vehicle checker detection property is set when straight line V=V0 is obtained when low flow volume region is located at below two curves or is not intersected, then Can require that wagon flow value is Q0 in unit detection time, makees straight line Q=Q0；Respectively with two curve intersection points, as each curve Corresponding lower limit of integral；

There is no actual data point at the lower limit of integral value for determining as stated above, then judge whether the lower limit of integral value is more than Q0, if judged result is " being more than ", the flow minimum in pretreated data sample of learning from else's experience is used as lower limit of integral value；

It is Q0 otherwise to obtain and set wagon flow value in vehicle checker detection performance requirement unit detection time, makees straight line Q=Q0；Respectively with Two curve intersection points, as the corresponding lower limit of integral of each curve；If lower limit of integral value does not have actual data point when being Q0 yet, The flow minimum learnt from else's experience in pretreated data sample is used as lower limit of integral value.

5. according to claim 4 based on section speed, the road traffic running status appraisal procedure of data on flows, it is special Levy and be：It is described that each the total interval of the corresponding integration of speed-discharge relation curve is used whole flow area in the period of different Between the mode of segment processing calculated, it is specific as follows：

If 1) need not ensure that each by stages Duan Jun is contrasted under equal flow, each contrast period curve integration It is total interval different, and it is always interval to ask for the integration of each curve；

[Q_a(i′),Q_b(i′)]；

Wherein, i '=1,2 ..., represent that the curve integration of different times is total interval, that is, contrast the curve integrating range in period i '；

2) if desired ensure that each by stages Duan Jun is contrasted under equal flow, i.e., each contrast period curve integration takes same One integration is total interval；The corresponding interval of each curve is obtained, each smaller lower limit in interval is selected as the total interval lower limit of integration, choosing Larger higher limit is selected as total interval upper limit is integrated, the total interval of the integration of each curve is taken as the total interval of same integration, it is described The same total interval [Q of integration_a,Q_b] establish rules really, it is：

If A. Q_a1≥Q_a2, and have Q_a1≤Q_a2, then taking the total interval of integration is：[Q_a2,Q_b2]；

If B. Q_a1＞ Q_a2, and have Q_b1＞ Q_b2, then taking the total interval of integration is：[Q_a2,Q_b1]；

If C. Q_a1＜ Q_a2, and have Q_b1＜ Q_b2, then taking the total interval of integration is：[Q_a1,Q_b2]；

If D. Q_a1≤Q_a2, and have Q_b1≥Q_b2, then taking the total interval of integration is：[Q_a1,Q_b1]。

6. according to claim 5 based on section speed, the road traffic running status appraisal procedure of data on flows, it is special Levy and be：Total interval [the Q of integration_a(i′),Q_b(i′)] or [Q_a,Q_b] bound of by stages is determined according to the following steps：

Each curve is integrated into total interval [Q away from m by certain group_a(i′),Q_b(i′)] or the total interval [Q of same integration_a,Q_b] it is divided into n Individual small by stages section, each small by stages section is specific as follows using each interval terminal as respective range of integration value：

1) if each total interval of curve integration is not to take same interval, point integrating range for each dividing is as follows：

[Q_a(i′),Q_a(i′)+m]；

[Q_a(i′)+m,Q_a(i′)+2m]；

[Q_a(i′)+2m,Q_a(i′)+3m]；

……；

[Q_a(i′)+(n-2)m,Q_a(i′)+(n-1)m]；

[Q_a(i′)+(n-1)m,Q_b(i′)]；

Above-mentioned institute by stages is each point of integrating range section of contrast period i ' curves integration, wherein, i '=1,2 ...；

If 2) each curve integrates total interval and takes same interval [Q_a,Q_b], then point integrating range for dividing is as follows：

[Q_a,Q_a+m]；

[Q_a+m,Q_a+2m]；

[Q_a+2m,Q_a+3m]；

……；

[Q_a+(n-2)m,Q_a+(n-1)m]；

[Q_a+(n-1)m,Q_b]。