朱普周,段宇洲,楊森森
(1.河南工業大學 土木工程學院,河南 鄭州 450001;2.長春工程學院 管理學院,吉林 長春 130021)
隨著我國經濟水平的不斷提高,城市交通擁堵現象日益增多,其中非機動車道內自行車及電動自行車造成的擁堵現象也日益嚴重。為此,國內外學者進行了大量的研究,其中賈順平[1-2]等通過實地調查,分析了電動自行車交通流三參數之間的關系,對非機動車進行了定量分析;孫明正等[3]在大量實測數據的基礎上,對信號控制交叉口作為間斷流交通運行特征的自行車流展開了較為深入的研究,分析了自行車的膨脹特性;周晨靜等[4]分析了信號交叉口自行車對臨近車道機動車通行能力的影響;Schleinitz等[5]發現不同類型的電動自行車與人力自行車的平均行駛速度有顯著差異;徐程[6-7]等針對傳統自行車和電動自行車混行路段,構建了基于高斯混合模型的速度分布函數,并通過引入機動車速度-密度關系模型,提出了非機動車交通流速度-密度關系模型。
本文在前人研究的基礎上,通過信號交叉口非機動車交通數據的調查與分析,研究信號交叉口非機動車速度與流量的關系,為非機動車交通流的研究提供思路,以減少交通沖突與擁堵。
利用Adobe After Effects CS4軟件,設定數據提取精度為25幀/s,時間精度為0.04 s,于2019年4月22日—26日共5個工作日的下午5:30—6:30的交通高峰時期,對鄭州市區黃河路與花園路、緯五路與花園路的兩個交叉路口進行數據提取,獲得了兩個信號交叉口非機動車流量、速度等數據。所調查的兩個信號交叉口南北方向禁止機動車、非機動車左轉,其中花園路與黃河路交叉口示意圖如圖1所示。

圖1 花園路與黃河路交叉口示意圖Fig.1 Schematic diagram of the intersection of Huayuan Road and Huanghe Road
在交通流參數統計時,經常使用車輛換算系數這一概念,車輛換算系數是指交通流中各類車型化為標準車的當量值。因此對于非機動車,其換算系數即為相同道路條件下單位電動自行車等效于單位自行車的數量。關于非機動車換算系數,葉曉飛等[8]基于有效面積和沖突事件模型提出兩種計算模型,最終得到電動自行車的換算系數為1.229 5;徐興等[9]結合非機動車有效行駛面積計算模型與超車率法,最終得出電動自行車對自行車的換算系數為1.346 2。從上述研究可以看到,電動自行車對自行車的換算系數都大于1.2且小于1.4,因此本文參考上述文獻最終選定電動自行車對自行車的換算系數為1.3。
由于不同寬度的非機動車道,非機動車的通行能力是不同的[10]。參照機動車流量定義,非機動車流量定義為單位時間內通過單位寬度非機動車道的非機動車數量,即[11]
(1)
式中:q為非機動車流量,輛/(s·m);Q為特定區域單位時間內通過非機動車道的車輛,輛;W為非機動車道的寬度,m;t為調查時間,s。
非機動車群通過信號交叉口時,受到信號燈的影響,集群性比較明顯。當綠燈剛剛啟亮,非機動車以飽和流率釋放,不規則地行駛在非機動車道上。由于不同寬度的非機動車道,飽和通過的非機動車流量是不同的,由此可將非機動車的飽和流率定義為單位時間內單位寬度的非機動車道飽和釋放的非機動車數量。
綠燈剛啟亮時,信號交叉口停車線處的非機動車以飽和狀態完全釋放,一般按照10 s的時間來計算飽和通過的非機動車數量。根據信號交叉口非機動車道的寬度,結合文獻[12]交通流研究與機動車交通流調查方法,非機動車飽和流率的計算公式為[11]:

(2)
式中:Qb為非機動車飽和流率,輛/(s·m);Nb為非機動車以飽和狀態通過信號交叉口的數量,輛;Δt為Nb輛非機動車以飽和狀態通過信號交叉口所用的時間,s。
在實際計算非機動車飽和流率時,由于非機動車具有集群性,在綠燈前期的短時間內就會釋放完成,因此通常選取高峰時段綠燈啟亮5 s的時間內通過的非機動車輛作為調查對象;另外,由于非機動車騎行者在綠燈啟亮前通常就會越過停車線占據到信號交叉口內,因此非機動車道寬度可以根據信號交叉口前排車輛實際占用的平均寬度確定。
非機動車的速度主要指自行車與電動自行車混合下的各自速度,通常考察指標有區間速度和平均速度。
2.3.1 區間速度
區間速度指的是非機動車通過某一區間的速度,又稱為行程車速,公式如下:
(3)
式中:vs為區間速度,m/s;t1、t2分別為非機動車駛入、駛出區間的時刻,s;L為該區間的距離,m。
2.3.2 平均速度
平均速度指的是非機動車行駛一定距離H與平均行駛時間的比值,如下式所示:
(4)

對不同類型非機動車分別進行統計,使用區間速度公式求取不同類型非機動車的速度;根據所求速度由平均速度公式得到不同類型非機動車間的速度差,以驗證不同比例非機動車流對交通流的影響。
各信號交叉口非機動車流量如表1所示。由表1可以看出,信號交叉口電動自行車占非機動車總數的比例均在70%以上,其中花園路—黃河路南進口道的電動自行車所占比例為87.69%。

表1 信號交叉口不同類型非機動車流量構成Table 1 Composition of traffic flow of different types of non-motor vehicles at signalized intersections
高峰時段綠燈啟亮5 s內通過信號交叉口4個進口道的非機動車數量如表2所示。由表2可以看出,非機動車飽和流率主要集中在0.76~0.97輛/(s·m),平均飽和流率為0.90輛/(s·m)。

表2 信號交叉口飽和流率統計表Table 2 Statistics of saturated flow rate at signalized intersections
非機動車速度數據如表3所示,不同路段非機動車的速度分布如圖2、圖3所示。由表3可以看出:黃河路—花園路的自行車平均速度為12.81 km/h、85%位車速為16.52 km/h;電動自行車平均速度為19.39 km/h、85%位車速為24.47 km/h。緯五路—花園路的自行車平均速度為13.73 km/h、85%位車速為18.15 km/h;電動自行車平均速度21.17 km/h、85%位車速為25.94 m/h。顯然,非機動車通過交叉口的速度差異較大,主要表現為:(1)電動自行車和自行車的平均速度差異較大,最大相差達7.44 km/h;(2)在綠燈中端或末端時,非機動車沒有加速啟動時間與非機動車群限制,速度較大,最大速度超過30 km/h。雖然如此,但是所調查信號交叉口電動自行車的85%位車速在25 km/h左右,基本滿足最新政策規定的電動自行車限速到25 km/h的要求。

表3 非機動車速度特性 Table 3 Speed characteristics of non-motor vehicles

圖2 花園路—黃河路不同類型非機動車速度分布情況Fig.2 The speed distribution of different types of non-motor vehicles on Huayuan Road-Huanghe Road

圖3 花園路—緯五路不同類型非機動車速度分布情況Fig.3 The speed distribution of different types of non-motor vehicles on Huayuan Road-Weiwu Road
由圖2、圖3可以看出:約90%的自行車速度小于18 km/h,約90%的電動自行車速度小于26 km/h;相較于花園路—黃河路,花園路—緯五路上的非機動車速度較大,最大速度超過34 km/h。顯然,這與兩者的非機動車道的寬度和流量有著較大聯系。
以花園路—黃河路為例進行速度-流量的回歸分析。在回歸分析之前要對電動自行車與自行車的流量進行統一化,即根據非機動車的不同換算系數,把兩者進行統一化;然后再根據統一化的流量數據,得到花園路—黃河路信號交叉口非機動車的速度-流量散點圖,如圖4所示。分析圖4,發現其速度-流量之間具有明顯的冪函數關系,因此建立數學模型如下:
vu=a1Sa2
(5)
式中:a1、a2為回歸系數;vu表示非機動車統一化后的速度,km/h;S表示非機動車統一化后的流量,由實際數據擬合得到,輛/(h·m)。

圖4 速度-流量關系擬合圖Fig.4 Fitting diagram of velocity-flow relationship
利用最小二乘法計算回歸模型中的具體參數值,并對回歸方程進行統計檢驗。模型摘要如表4所示,回歸方程的顯著性檢驗如表5所示,回歸系數的顯著性檢驗如表6所示,圖5為標準化殘差。
表5回歸方程檢驗、表6回歸系數檢驗(B值表示回歸系數與截距)的顯著性均為0.000,小于0.05,說明速度與流量之間存在相應關系且具有顯著意義,圖5回歸標準殘差曲線服從正態分布,能夠較好地反映回歸方程相應變量的變化規律;表4在未去除全部奇異值時,相關系數為0.868,判定系數R2為0.753>0.7,表明擬合效果較好。

表4 模型摘要Table 4 Model summary

表5 回歸方程的顯著性檢驗Table 5 Significance test of regression coefficients

表6 模型回歸系數分析Table 6 Model regression coefficient analysis
根據調查地點的速度-流量數據,得到回歸方程為v=199.34S-0.352。

圖5 回歸方程標準化殘差Fig.5 Standardized residuals of regression equations
通過對非機動車交通流速度、流量的調查分析,發現其速度-流量之間具有明顯的冪函數關系,回歸方程為v=199.34S-0.352,R2為0.753﹥0.7,顯著性<0.05,擬合效果較優。