• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    有樁自行車(chē)共享系統(tǒng)吞吐率的近似模型及算法

    2023-12-17 01:59:32王靜妍張勇
    山東科學(xué) 2023年6期
    關(guān)鍵詞:運(yùn)營(yíng)效率

    王靜妍 張勇

    摘要:通過(guò)探討用戶(hù)需求和騎行時(shí)間不確定的有樁自行車(chē)共享系統(tǒng)(docked bike-sharing system, DBSS)建立吞吐率的近似模型及算法。一個(gè)具有固定自行車(chē)數(shù)量的DBSS可視為封閉的排隊(duì)網(wǎng)絡(luò),每個(gè)站點(diǎn)都是有限的M/M/1隊(duì)列,由此建立了DBSS吞吐率的近似模型及其算法。該算法不僅能夠計(jì)算道路上期望的自行車(chē)數(shù)量、騎行時(shí)間、車(chē)站的期望庫(kù)存及停留時(shí)間,還能計(jì)算最優(yōu)自行車(chē)投放量,即吞吐率最大值對(duì)應(yīng)的最小的自行車(chē)投放量。同時(shí),給出了給定用戶(hù)需求、路由矩陣和車(chē)樁分配下站點(diǎn)自行車(chē)集聚與空缺的判斷方法。將該近似算法在真實(shí)的DBSS中進(jìn)行了應(yīng)用。結(jié)果表明,隨著自行車(chē)投放量的增加,系統(tǒng)吞吐率呈階梯形遞增但存在上限;自行車(chē)投放量一旦超過(guò)最優(yōu)數(shù)量將產(chǎn)生閑置,并且自行車(chē)集聚與空缺站點(diǎn)分布也將固定。

    關(guān)鍵詞:有樁自行車(chē)共享系統(tǒng);運(yùn)營(yíng)效率;封閉排隊(duì)網(wǎng)絡(luò);吞吐率;空滿(mǎn)樁站點(diǎn);自行車(chē)投放量

    中圖分類(lèi)號(hào):U491?? 文獻(xiàn)標(biāo)志碼:A?? 文章編號(hào):1002-4026(2023)06-0074-12

    An approximate model and algorithm for throughput rate of

    a docked bike-sharing system

    WANG Jingyan, ZHANG Yong*

    (School of Rail Transportation, Soochow University, Suzhou 215131, China)

    Abstract∶In this paper, an approximate model and algorithm for the throughput rate are established by studying a docked bike-sharing system (DBSS) using stochastic user demands, routing matrix, and cycling times. A DBSS with a fixed number of bikes can be considered a closed queuing network with a buffered M/M/1 queue at each station, thus establishing an approximate model and algorithm for the throughput rate of DBSS. This algorithm can calculate the average number of bikes on roads and at stations. Moreover, it can estimate the average cycling time on roads and bike dwell time at stations and further determine the optimal number of bikes achieving the maximum throughput rate in the DBSS. Additionally, this paper proposes a method to determine whether a station is a bike surplus station or a bike deficient station under given user demands, routing matrix, cycling time matrix, and dock allocation. Finally, the approximate algorithm is verified in a real-world DBSS. The results show that the throughput rate of the DBSS increases in a step-wise manner with the increasing bike input under an superior limit. When the number of bike inputs exceeds the optimal quantity, there will be idle bikes, and the spatial distribution of bike surplus stations and bike deficient stations will remain unchanged.

    Key words∶docked bike-sharing system; operational efficiency; closed queuing network;system throughput rate; bike surplus and deficient stations; bike input

    有樁自行車(chē)共享系統(tǒng)(docked bike-sharing system,DBSS)是基于固定站點(diǎn)提供借還自行車(chē)服務(wù)的服務(wù)系統(tǒng)[1]。近年來(lái)在全球得到了廣泛應(yīng)用,截至2022年12月,全球建設(shè)了1 900多個(gè)有樁自行車(chē)共享系統(tǒng),自行車(chē)運(yùn)營(yíng)數(shù)量超過(guò)897萬(wàn)輛[2]。公共自行車(chē)對(duì)于解決公共交通“最后一公里”難題,倡導(dǎo)綠色出行發(fā)揮了重要作用。而現(xiàn)實(shí)中自行車(chē)站點(diǎn)常常出現(xiàn)無(wú)車(chē)可借、無(wú)樁可還的情況。此時(shí),如何計(jì)算自行車(chē)系統(tǒng)的服務(wù)能力以及在規(guī)劃DBSS時(shí),如何判斷哪些站點(diǎn)將發(fā)生空滿(mǎn)樁現(xiàn)象,對(duì)這些問(wèn)題的研究對(duì)于科學(xué)開(kāi)展DBSS規(guī)劃與運(yùn)營(yíng)調(diào)度具有重要的意義。

    目前對(duì)于DBSS的研究主要聚焦主題有:自行車(chē)道網(wǎng)絡(luò)設(shè)計(jì)[3]、自行車(chē)站選址設(shè)計(jì)[4]、車(chē)隊(duì)規(guī)模設(shè)計(jì)[5]、自行車(chē)搬遷[6-7]、自行車(chē)庫(kù)存水平管理[8-9]以及與公共交通的整合[10]。這些研究試圖通過(guò)優(yōu)化算法實(shí)現(xiàn)自行車(chē)共享系統(tǒng)的最優(yōu)設(shè)計(jì)與運(yùn)營(yíng)。DBSS有兩個(gè)顯著的特性:(1)用戶(hù)到達(dá)站點(diǎn)的時(shí)間間隔是隨機(jī)的,自行車(chē)從租車(chē)點(diǎn)到還車(chē)點(diǎn)間的騎行時(shí)間因用戶(hù)騎行的異質(zhì)性也呈隨機(jī)性,故用戶(hù)還車(chē)間隔也是隨機(jī)的;(2)投入到自行車(chē)共享系統(tǒng)的自行車(chē)總數(shù)是固定的,且只能在站點(diǎn)停留或在站點(diǎn)間循環(huán)。因此,所有的自行車(chē)只有兩種狀態(tài),在站點(diǎn)停留和在道路騎行。從自行車(chē)角度來(lái)看,這兩個(gè)特征決定了DBSS本質(zhì)上屬于封閉排隊(duì)網(wǎng)絡(luò)。

    為了提高DBSS的運(yùn)營(yíng)效率,目前的研究主要考慮不確定的有樁自行車(chē)共享系統(tǒng)規(guī)劃。George等[11]基于封閉排隊(duì)網(wǎng)絡(luò)研究了自行車(chē)可用性與車(chē)隊(duì)規(guī)模問(wèn)題,基于利潤(rùn)最大化給出了最優(yōu)車(chē)隊(duì)規(guī)模。Fricker等[5]提出了同質(zhì)自行車(chē)系統(tǒng)的隨機(jī)模型,研究了用戶(hù)到達(dá)不確定性對(duì)“空滿(mǎn)樁”現(xiàn)象的影響及最優(yōu)車(chē)隊(duì)規(guī)模問(wèn)題。這些研究盡管考慮了需求的不確定性,但忽略了站點(diǎn)停車(chē)樁容量限制,從而無(wú)法準(zhǔn)確分析有樁自行車(chē)共享系統(tǒng)服務(wù)能力和水平。為此,部分研究引入了停車(chē)樁容量限制,考慮借還車(chē)失敗次數(shù)[12-13]、站點(diǎn)空滿(mǎn)樁持續(xù)時(shí)間[14-15]等約束,尋求用戶(hù)需求最大化目標(biāo)下的最優(yōu)站點(diǎn)選址與容量分配方案。但大部分研究未能考慮站點(diǎn)容量限制,無(wú)法揭示自行車(chē)系統(tǒng)的服務(wù)水平與自行車(chē)、停車(chē)樁分配之間的數(shù)學(xué)關(guān)系,也就無(wú)法探索“無(wú)車(chē)可借,無(wú)樁可還”的發(fā)生問(wèn)題。即使有些研究考慮了站點(diǎn)容量限制,但是還沒(méi)有公開(kāi)的解決方案和方法來(lái)量化DBSS的吞吐率,與現(xiàn)實(shí)仍存在一定的差距。

    本文對(duì)DBSS的服務(wù)能力進(jìn)行了定義,構(gòu)建了近似的DBSS模型,開(kāi)發(fā)出一種用于計(jì)算系統(tǒng)的服務(wù)績(jī)效近似算法;同時(shí),探索了自行車(chē)系統(tǒng)的運(yùn)營(yíng)性質(zhì),包括系統(tǒng)吞吐率與自行車(chē)投放量的關(guān)系,以及站點(diǎn)出現(xiàn)自行車(chē)集聚與空缺的判斷條件;最后通過(guò)一個(gè)實(shí)際案例驗(yàn)證提出的近似算法和結(jié)論,展示其在DBSS中的初步應(yīng)用。

    1 模型構(gòu)建

    1.1 問(wèn)題描述及假設(shè)

    考慮一個(gè)有樁自行車(chē)共享系統(tǒng),衡量其服務(wù)效率最直接指標(biāo)為單位時(shí)間內(nèi)服務(wù)的有效用戶(hù)人數(shù),即封閉排隊(duì)網(wǎng)絡(luò)中的吞吐率[16],因此本文將在封閉排隊(duì)網(wǎng)絡(luò)理論框架下求解有樁自行車(chē)共享系統(tǒng)的吞吐率,即用戶(hù)借車(chē)成功離開(kāi)站點(diǎn)、真正得到站點(diǎn)提供服務(wù)的速率。根據(jù)上述自行車(chē)借還過(guò)程,為便于建模特作如下假設(shè):

    A1 有樁自行車(chē)共享系統(tǒng)由m個(gè)站點(diǎn)和n輛自行車(chē)組成,每個(gè)站點(diǎn)有Bi(i =1,2,…,m)個(gè)停車(chē)樁,并且任意一對(duì)自行車(chē)站點(diǎn)之間至少可以雙向路徑連通;

    A2 假設(shè)用戶(hù)到達(dá)站點(diǎn)i的時(shí)間間隔服從參數(shù)λi(i=1,2,…,m)的負(fù)指數(shù)分布,各站點(diǎn)用戶(hù)相互獨(dú)立,用戶(hù)在站點(diǎn)i以先到先服務(wù)規(guī)則借還車(chē);

    A3 用戶(hù)從站點(diǎn)i騎自行車(chē)到站點(diǎn)j (i, j=1,2,…,m)的選擇概率為pij≥0,且∑mj=1pij=1,本文稱(chēng)P=pij為路由矩陣;

    A4 用戶(hù)在站點(diǎn)i、 j間騎行的旅行時(shí)間服從參數(shù)為τij的獨(dú)立負(fù)指數(shù)分布;

    A5 對(duì)于用戶(hù)借還車(chē)活動(dòng),若站點(diǎn)沒(méi)有可用的自行車(chē)時(shí),用戶(hù)會(huì)離開(kāi)系統(tǒng);對(duì)于用戶(hù)還車(chē)活動(dòng),若目的地站點(diǎn)滿(mǎn)樁時(shí),用戶(hù)會(huì)騎行至該站點(diǎn)最近站點(diǎn)還車(chē),直至成功為止。

    1.2 封閉排隊(duì)網(wǎng)絡(luò)分析

    對(duì)于有樁自行車(chē)共享系統(tǒng)構(gòu)成的封閉排隊(duì)網(wǎng)絡(luò),用戶(hù)在站點(diǎn)i完成借車(chē)后以概率pijj=1,…,m騎行至站點(diǎn)j,該概率滿(mǎn)足∑mj=1pij=1 (i=1,…,m),即路由矩陣P是一個(gè)馬爾可夫轉(zhuǎn)移矩陣,假定該矩陣不可約,以π=π1,π2,…,πm記這個(gè)馬爾科夫鏈的平穩(wěn)概率,即π是

    πj=∑mi=1πipij,? ∑mj=1πj=1(1)

    的唯一正解。

    如果記站點(diǎn)j的自行車(chē)平均到達(dá)率或平均離開(kāi)率為aj(n),其中j=1,2,…,m,那么滿(mǎn)足

    aj(n)=∑mi=1ai(n)pij。(2)

    由平穩(wěn)概率方程可以得到

    aj(n)=a(n)πj,其中a(n)=∑mj=1aj(n)。(3)

    從該方程可以看到,a(n)是整個(gè)系統(tǒng)的自行車(chē)平均離開(kāi)速率,其包括了系統(tǒng)的有效吞吐速率和還車(chē)失敗速率兩部分。

    引理1 如果X1,…,Xq是服從均值為λi的獨(dú)立的泊松隨機(jī)變量(q=1, 2, …, m),那么Y:=∑qi=1Xi服從均值為∑qi=1λi的泊松分布。

    證明 采用數(shù)學(xué)歸納法證明。設(shè)X1~Possionλ1,X2~Possionλ2,X1和X2相互獨(dú)立,那么對(duì)于任意y∈N,

    pX1+X2y=∑xpX1kpX2y-k=∑yk=0λk1e-λ1k!·λy-k2e-λ2y-k!

    =e-λ1+λ2·λy2·∑yk=0λ1λ2kk!y-k! =e-λ1+λ2y!λ1+λ2y。(4)

    因此,X1+X2~Possionλ1+λ2?,F(xiàn)在,假設(shè)X1,…,Xq-1是服從均值為λi( i =1,2,…,q-1)的獨(dú)立的泊松隨機(jī)變量,那么X′:=∑q-1i=1Xi服從均值為∑q-1i=1λi的泊松分布。若Xq是獨(dú)立于X′且服從均值為λq的泊松隨機(jī)變量,則由公式(4)可得,

    pX′+Xqy=e-∑q-1i=1λi+λqy!∑q-1i=1λi+λqy,

    p∑qi=1Xiy=e-∑qi=1λiy!∑qi=1λiy。

    因此,∑qi=1Xi~Possion∑qi=1λi,證明完畢。

    命題1 在有樁自行車(chē)共享系統(tǒng)中,如果用戶(hù)到達(dá)站點(diǎn)i服從速率為λi的獨(dú)立的泊松分布,那么到達(dá)站點(diǎn)i的自行車(chē)數(shù)量也是獨(dú)立的泊松隨機(jī)變量。

    證明 根據(jù)自行車(chē)借還過(guò)程可知,到站的自行車(chē)分為兩種類(lèi)型:第一種是用戶(hù)自愿歸還至目的地站點(diǎn)的自行車(chē);第二種是用戶(hù)因目的地站點(diǎn)滿(mǎn)樁而被迫歸還至該站點(diǎn)的自行車(chē)。為了證明自行車(chē)到達(dá)站點(diǎn)的過(guò)程服從泊松分布,需通過(guò)引理1證明自愿用戶(hù)和被迫用戶(hù)的到達(dá)過(guò)程都服從獨(dú)立的泊松分布。

    如果用戶(hù)成功將自行車(chē)歸還到意愿目的地站點(diǎn),則稱(chēng)這類(lèi)用戶(hù)為自愿用戶(hù)。本文首先證明自愿用戶(hù)到達(dá)站點(diǎn)服從獨(dú)立的泊松分布。假設(shè)用戶(hù)按速率為λi (i=1,2,…,m) 的泊松過(guò)程到達(dá)站點(diǎn)i,忽略用戶(hù)從借到車(chē)至離開(kāi)站點(diǎn)所花費(fèi)的時(shí)間。一個(gè)在時(shí)刻ss≤t到達(dá)站點(diǎn)i的用戶(hù),以概率Prent, 1is借車(chē)成功,以概率Prent, 2is借不到車(chē)(其中∑2j=1Prent,? jis=1)。利用文獻(xiàn)[16]中命題5.3可知,直到時(shí)刻t為止,在站點(diǎn)i成功借車(chē)用戶(hù)數(shù)Nrent, 1it是具有均值ENrent, 1it=λi∫t0Prent, 1isds的獨(dú)立泊松隨機(jī)變量。同樣,將在站點(diǎn)i (i=1,2,…,m)成功借車(chē)的用戶(hù)分為m 類(lèi),每類(lèi)用戶(hù)是從站點(diǎn)i自愿還車(chē)至每個(gè)站點(diǎn)的用戶(hù)。由于借車(chē)成功的用戶(hù)從每個(gè)站點(diǎn)出發(fā)是泊松過(guò)程,因此從站點(diǎn)i(i=1,2,…,m)至還車(chē)站點(diǎn)j(j=1,2,…,m)的自愿用戶(hù)的到達(dá)過(guò)程也是獨(dú)立泊松過(guò)程。由引理1可知,任意站點(diǎn)的自愿用戶(hù)的到達(dá)過(guò)程也是一個(gè)獨(dú)立的泊松過(guò)程。

    同理可得,任意站點(diǎn)的被迫用戶(hù)的到達(dá)過(guò)程也是一個(gè)獨(dú)立的泊松過(guò)程。綜上,任意站點(diǎn)的自行車(chē)到達(dá)數(shù)為自愿用戶(hù)數(shù)和被迫用戶(hù)數(shù)的總和。利用引理1可得,任意站點(diǎn)的自行車(chē)到達(dá)服從泊松分布。

    值得注意的是,由命題1可知自行車(chē)到達(dá)站點(diǎn)過(guò)程服從泊松分布,那么自行車(chē)到達(dá)站點(diǎn)的時(shí)間間隔服從負(fù)指數(shù)分布。因此有樁自行車(chē)共享系統(tǒng)在時(shí)間長(zhǎng)度上是一個(gè)無(wú)記憶的系統(tǒng)。

    基于經(jīng)典的平均值分析法[17],可通過(guò)兩個(gè)主要定理(到達(dá)定理和Little定律)推導(dǎo)有樁自行車(chē)共享系統(tǒng)的吞吐率的近似算法。根據(jù)到達(dá)定理[16]可知,在有n輛自行車(chē)的封閉網(wǎng)絡(luò)系統(tǒng)中,騎自行車(chē)到達(dá)站點(diǎn)i的到達(dá)者所看到的系統(tǒng)的分布與只有n-1輛自行車(chē)的同樣的網(wǎng)絡(luò)系統(tǒng)的平穩(wěn)分布相同。因此,令ENi(n)為網(wǎng)絡(luò)中有n輛自行車(chē)時(shí)停在站點(diǎn)i的平均自行車(chē)數(shù)量;ETi(n)為網(wǎng)絡(luò)中有n輛自行車(chē)時(shí)在站點(diǎn)i的自行車(chē)平均停留時(shí)間。對(duì)一輛到達(dá)自行車(chē)看到在站點(diǎn)i的自行車(chē)數(shù)量Ci取期望值,推出

    ETi(n)=1+E[Ci]λi=1+ENi(n-1)λi。(5)

    系統(tǒng)中的自行車(chē)包括在站點(diǎn)停放和正在被使用兩類(lèi)自行車(chē),故n=∑mi=1ENi(n)+∑mi=1∑mj=1ENij(n)。而在站點(diǎn)i的自行車(chē)停留的時(shí)間為ETi(n),其他站點(diǎn)騎行至站點(diǎn)i的平均旅行時(shí)間為∑mj=1pjiτji。因此,系統(tǒng)的自行車(chē)平均到達(dá)率為

    an=n/∑mi=1πiETi(n)+∑mi=1∑mj=1πipijτij。(6)

    因此站點(diǎn)的有效吞吐速率為

    Hi(n)=min a(n)πi,λi,

    其中

    Hn=∑mi=1Hi(n)。(7)

    利用Little公式可以得到站點(diǎn)i的平均自行車(chē)數(shù)量為

    ENi(n)=min Hi(n)ETi(n),Bi。(8)

    根據(jù)命題1可知,對(duì)于擁有Bi個(gè)停車(chē)樁站點(diǎn),自行車(chē)的借還過(guò)程可以看作是用戶(hù)與自行車(chē)的雙排隊(duì)系統(tǒng)。若站點(diǎn)無(wú)自行車(chē),則用戶(hù)到達(dá)站點(diǎn)后即刻離開(kāi)系統(tǒng),故站點(diǎn)空樁就不存在顧客排隊(duì)。因此,任意站點(diǎn)i可以看作是M/M/1/Bi排隊(duì)系統(tǒng),狀態(tài)轉(zhuǎn)移如圖1所示。

    在該排隊(duì)系統(tǒng)中,將用戶(hù)視為服務(wù)員,到達(dá)的自行車(chē)視為作業(yè),則站點(diǎn)i處可停放的自行車(chē)上限為Bi個(gè),因此有

    λn=ai,????? n=0,1,…,Bi-1 0,??????? n=Bi,? Bi+1,… ,

    μn=λi,???? n=0,1,…。

    根據(jù)流入和流出每個(gè)狀態(tài)的期望速率必須相等的條件,可以得到狀態(tài)k的平衡方程:

    aiPxi=k-1+λiPxi=k+1=ai+λiPxi=k,??? k=1,2,…,Bi-1。

    結(jié)合∑Bik=0Pxi=k=1,即可得到站點(diǎn)i無(wú)車(chē)和滿(mǎn)樁的概率分別為

    Pxi=0=1-ρi1-ρiBi+1,?????? ρi≠1 1Bi+1,??????????? ρi=1 ?,(9)

    Pxi=Bi=1-ρiρiBi1-ρiBi+1,??? ρi≠1 1Bi+1,??????????? ρi=1 。(10)

    因此,當(dāng)系統(tǒng)的總自行車(chē)投放量超過(guò)站點(diǎn)的最小容量時(shí),站點(diǎn)開(kāi)始存在用戶(hù)還車(chē)失敗的情形。根據(jù)假設(shè)A5可知用戶(hù)需要去目的地站點(diǎn)的最近站點(diǎn)繼續(xù)還車(chē),自行車(chē)網(wǎng)絡(luò)的實(shí)際轉(zhuǎn)移概率發(fā)生變化,即站點(diǎn)j的車(chē)輛平均到達(dá)率由任意站點(diǎn)的意愿轉(zhuǎn)移速率和被迫轉(zhuǎn)移速率組成。故利用公式(9)~(10)可得由站點(diǎn)i向站點(diǎn)j的實(shí)際自行車(chē)轉(zhuǎn)移概率為

    pij=λi1-Pxi=0pij+aiPxi=Biηij∑mj=1λi1-Pxi=0pij+aiPxi=Biηij=λiρi-ρiBi+1pij+ai1-ρiρiBiηij∑mj=1λiρi-ρiBi+1pij+ai1-ρiρiBiηij, (11)

    其中,若站點(diǎn)j 是站點(diǎn)i的最近站點(diǎn),則ηij=1;否則ηij=0。

    1.3 吞吐率的近似算法

    利用公式(5)~(11)可以建立計(jì)算自行車(chē)網(wǎng)絡(luò)中站點(diǎn)的平均自行車(chē)數(shù)量、平均等待時(shí)間和吞吐率的近似算法,步驟如下:

    步驟1 計(jì)算初始的馬爾科夫鏈平穩(wěn)概率,πj=∑mi=1πipij,∑mj=1πj=1;

    步驟2 Let ENi(0)=0,i=0,…,m and a(0)=0;

    步驟3 For k=1,…,n do

    ETi(k)=ENi(k-1)+1/λi;

    If k≤min (Bi)then

    ak=k/∑mi=1πiETi(k)+∑mi=1∑mj=1πipijτij;

    ai(k)=a(k)πi;

    Else

    pij=λiai(k-1)λi-ai(k-1)λiBi+1pij+ai(k-1)1-ai(k-1)λiai(k-1)λiBiηij∑mj=1λiai(k-1)λi-ai(k-1)λiBi+1pij+ai(k-1)1-ai(k-1)λiai(k-1)λiBiηij;

    πj=∑mi=1πipij,??? ∑mi=1πj=1;

    ak=k/∑mi=1πiETi(k)+∑mi=1∑mj=1πipijτij;

    ai(k)=a(k)πi;

    End

    Hi(k)=min ai(k),λi;

    Hk=∑mi=1Hik;

    ENik=min Hi(k)ETi(k),Bi;

    End

    1.4 系統(tǒng)的運(yùn)營(yíng)性質(zhì)

    本小節(jié)通過(guò)考察系統(tǒng)吞吐率與投放自行車(chē)數(shù)量、停車(chē)樁數(shù)之間的相互關(guān)系,揭示站點(diǎn)產(chǎn)生“無(wú)車(chē)可借,無(wú)樁可還”的發(fā)生機(jī)制,并給出最優(yōu)自行車(chē)數(shù)量的定義,為運(yùn)營(yíng)商提升有樁自行車(chē)共享系統(tǒng)的服務(wù)效率,促進(jìn)交通綠色可持續(xù)發(fā)展提供決策依據(jù)。

    1.4.1 吞吐率

    命題2 在一個(gè)有樁自行車(chē)共享系統(tǒng)中,如果自行車(chē)到達(dá)站點(diǎn)的時(shí)間越早,那么用戶(hù)騎自行車(chē)離開(kāi)該站點(diǎn)的時(shí)間也越早。

    證明 令A(yù)ik為第k 輛自行車(chē)到達(dá)站點(diǎn)i的時(shí)間,Dik為到達(dá)站點(diǎn)i的第k輛自行車(chē)的離開(kāi)時(shí)間,Wik為到達(dá)站點(diǎn)i的第k輛自行車(chē)的停留時(shí)間。命題2意味著,對(duì)于站點(diǎn)i,若An+1ik≤Anik,k=1,2,…,K,則Dn+1ik≤Dnik,k=1,2,…,K。命題2采用歸納法進(jìn)行證明。對(duì)于K=1, Dn+1i1=An+1i1+Wi1≤Ani1+Wi1=Dni1。假設(shè)引理針對(duì)K=l成立,即Dn+1il≤Dnil。那么對(duì)于K=l + 1, Dn+1i(l+1)=max Dn+1il,An+1i(l+1)+Wi(l+1)≤max Dnil,Ani(l+1)+Wi(l+1)=Dni(l+1),證明完畢。

    性質(zhì)1 在給定停車(chē)樁的情況下,有樁自行車(chē)共享系統(tǒng)吞吐率隨自行車(chē)數(shù)量的增加呈非遞減性,即對(duì)于所有的t≥0,有H(n+1)≥H(n)。

    證明 采用歸納法,在有n輛公共自行車(chē)和m個(gè)站點(diǎn)的系統(tǒng)中,令τn1<τn2<…為至少一輛自行車(chē)完成服務(wù)的瞬時(shí)時(shí)間,因此將tn定義為在有n輛自行車(chē)和n + 1輛自行車(chē)系統(tǒng)中至少有一輛自行車(chē)完成服務(wù)的時(shí)間序列,其中t0:=0;ts:=min minτniτni>ts-1,minτn+1iτn+1i>ts-1,s≥1。為了證明DBSS吞吐率隨系統(tǒng)自行車(chē)數(shù)量的增加呈非遞減性,即對(duì)于t≥0有H(n+1)≥H(n),由于H(n)=limt→∞∑mi=1Dni(t)/t,意味著要證明Dn+1i(t)≥Dni(t),i=1,2,…,m。

    當(dāng)tk= t0 = 0時(shí),Dn+1i(0)=Dni(0)=0,因此Dn+1i(t)≥Dni(t)成立。假設(shè)Dn+1i(t)≥Dni(t)對(duì)于t0, t1,…,tq成立。令Slj(t)為站點(diǎn)l中第j個(gè)離開(kāi)的自行車(chē)在服務(wù)完成后將去往的站點(diǎn)。若i=j,則s(i, j)=1;否則為0。Tijk為從站點(diǎn)i到站點(diǎn)j的第k輛自行車(chē)的行程時(shí)間,Xi(0) 為站點(diǎn)i的初始自行車(chē)數(shù)量。利用命題2,對(duì)于k=1,2,…,q; b = 1, 2,…, m有

    An+1i(tk)=∑ml=1∑Dn+1l(tk-Tlij)j=1σ(Slj(tk-Tlij),i)+Xi(0)+1·σ(b,i)

    ≥∑ml=1∑Dnl(tk-Tlij)j=1σ(Slj(tk-Tlij),i)+Xi(0)=Ani(tk)。

    因?yàn)锳n+1i(tk)和Ani(tk)在[tk-1,tk)是恒定的,因此對(duì)于t≤tq,An+1i(t)≥Ani(t)。所以得到An+1ij≤Anij,j=1,…,Ani(tq)。根據(jù)命題2可知Dn+1ij≤Dnij,j=1,…,Ani(tq)。由于Wij>0,Anij≤tq,對(duì)于所有的j有Dnij≤tq+1,那么也遵循Dn+1ij≥Dnij,對(duì)于所有的站點(diǎn)j 有Dnij≤tq+1。因此Dn+1i(t)≥Dni(t)對(duì)于tq+1也成立,證明完畢。

    性質(zhì)1表明有樁自行車(chē)共享系統(tǒng)的有效吞吐率隨著自行車(chē)數(shù)量的遞增呈非遞減性。然而,本文將在第2.2節(jié)中看到,有樁自行車(chē)共享系統(tǒng)的有效吞吐率隨著自行車(chē)數(shù)量的增加而增加,但存在一個(gè)上限。因此,可能存在一個(gè)自行車(chē)臨界值阻止吞吐率的增加。通過(guò)以下的定義,可以得出在DBSS中最優(yōu)的自行車(chē)數(shù)量以及如何尋找該最優(yōu)值。

    定義1 最優(yōu)自行車(chē)數(shù)量:對(duì)于一個(gè)有樁自行車(chē)共享系統(tǒng),若滿(mǎn)足以下條件,則可以確定系統(tǒng)的最優(yōu)自行車(chē)數(shù)量

    nOPT=max kHk-Hk-1>0? 且? Hk+1-Hk=0且k∈argmax H(l),l∈Z+,k∈Z+。

    值得注意的是,在本文設(shè)計(jì)的近似算法中,迭代計(jì)算了H(k)和H(k-1),因此系統(tǒng)的最優(yōu)自行車(chē)數(shù)量很容易確定。

    1.4.2 空滿(mǎn)樁的判斷

    在有樁自行車(chē)共享系統(tǒng)中,站點(diǎn)所有的車(chē)樁如果長(zhǎng)時(shí)間出現(xiàn)空樁,用戶(hù)將無(wú)法借車(chē);如果長(zhǎng)時(shí)間出現(xiàn)滿(mǎn)樁,用戶(hù)將無(wú)法還車(chē)。如果自行車(chē)共享系統(tǒng)多個(gè)站點(diǎn)長(zhǎng)時(shí)間空滿(mǎn)樁將顯著惡化用戶(hù)體驗(yàn),可能迫使公共自行車(chē)用戶(hù)轉(zhuǎn)移到其他交通方式。因此,有必要診斷有樁自行車(chē)共享系統(tǒng)站點(diǎn)是否會(huì)發(fā)生空滿(mǎn)樁。

    定義2 對(duì)于一個(gè)帶有Bi個(gè)停車(chē)樁的站點(diǎn)i,若滿(mǎn)足

    (1) 站點(diǎn)i的平均自行車(chē)數(shù)量低于φ1Bi的概率大于ω1,即Pxi<φ1Bi」>ω1,則稱(chēng)該站點(diǎn)為自行車(chē)空缺的;

    (2) 站點(diǎn)i的平均自行車(chē)數(shù)量高于φ2Bi的概率大于ω2,即Pxi>φ2Bi」>ω2,則稱(chēng)該站點(diǎn)為自行車(chē)集聚的。

    其中φ1和φ2為比例參數(shù),且0<φ1,φ2,ω1,ω2<1;符號(hào)“·」”表示向下取整數(shù)。

    命題3 自行車(chē)集聚與空缺的判斷條件:對(duì)于任意站點(diǎn)i,當(dāng)自行車(chē)到達(dá)率ai和用戶(hù)到達(dá)率λi滿(mǎn)足下列條件(1)~ (3)中的任意一條,可判斷該站點(diǎn)的自行車(chē)庫(kù)存情況:

    (1) 若滿(mǎn)足下面的條件,則站點(diǎn)i會(huì)產(chǎn)生自行車(chē)集聚情況

    aφ2Bi」+1i·λBi-φ2Bi」i-aBi+1iλBi+1i-aBi+1i>ω2,當(dāng)aiλi≠1時(shí)Bi-φ2Bi」Bi+1>ω2,???????????????????? 當(dāng)aiλi=1時(shí);(12)

    (2) 若滿(mǎn)足下面的條件,則站點(diǎn)i會(huì)產(chǎn)生自行車(chē)空缺情況

    λBi+1i-aφ1Bi」i·λBi-φ1Bi」+1iλBi+1i-aBi+1i>ω1,當(dāng)aiλi≠1時(shí)φ1Bi」Bi+1>ω1,??????????????????????????? 當(dāng)aiλi=1時(shí);(13)

    (3) 若滿(mǎn)足下面的條件,則站點(diǎn)i自行車(chē)既不集聚也不空缺

    λBi+1i-aφ1Bi」i·λBi-φ1Bi」+1iλBi+1i-aBi+1i≤ω1且aφ2Bi」+1i·λBi-φ2Bi」i-aBi+1iλBi+1i-aBi+1i≤ω2, 當(dāng)aiλi≠1時(shí)φ1Bi」Bi+1≤ω1且Bi-φ2Bi」Bi+1≤ω2,??????????????????????????????????????????????? 當(dāng)aiλi=1時(shí)。(14)

    證明 在平穩(wěn)狀態(tài)下,若站點(diǎn)i自行車(chē)產(chǎn)生集聚,需要滿(mǎn)足

    Pxi>φ2Bi」>ω2。(15)

    當(dāng)ρi≠1(其中ρi=ai/λi)時(shí),則不等式(15)可化簡(jiǎn)為Pxi>φ2Bi」=ρφ2Bi」+1i-ρBi+1i/1-ρBi+1i>ω2,即滿(mǎn)足條件ρφ2Bi」+1i-ρBi+1i/1-ρBi+1i>ω,

    則站點(diǎn)i自行車(chē)會(huì)產(chǎn)生集聚。當(dāng)ρi=1時(shí),則不等式(15)可化簡(jiǎn)為Pxi>φ2Bi」=Bi-φ2Bi」/Bi+1>ω2,即滿(mǎn)足條件Bi-φ2Bi」/Bi+1>ω2,則站點(diǎn)i自行車(chē)會(huì)產(chǎn)生集聚。站點(diǎn)自行車(chē)產(chǎn)生空缺或既不集聚也不空缺的證明與上述過(guò)程類(lèi)似,不再贅述。

    推論1 若定義自行車(chē)空缺站點(diǎn)為平均自行車(chē)數(shù)量為0的概率大于ω,即Pxi=0>ω;自行車(chē)集聚站點(diǎn)為平均自行車(chē)數(shù)量為Bi的概率大于ω,即Pxi=Bi>ω。那么對(duì)于任意站點(diǎn)i,當(dāng)自行車(chē)到達(dá)率ai和用戶(hù)到達(dá)率λi滿(mǎn)足下列條件(1)~ (3)中的任意一條,可判斷該站點(diǎn)的自行車(chē)庫(kù)存情況:

    (1) 若滿(mǎn)足下面的條件,則站點(diǎn)i會(huì)產(chǎn)生自行車(chē)集聚情況

    aBii∑Bik=0λBi-kiaki>ω,當(dāng)aiλi≠1時(shí)1Bi+1>ω,???????????? 當(dāng)aiλi=1時(shí);(16)

    (2) 若滿(mǎn)足下面的條件,則站點(diǎn)i會(huì)產(chǎn)生自行車(chē)空缺情況

    λBii∑Bik=0λBi-kiaki>ω,當(dāng)aiλi≠1時(shí)1Bi+1>ω,???????????? 當(dāng)aiλi=1時(shí);(17)

    (3) 若滿(mǎn)足下面的條件,則站點(diǎn)i自行車(chē)既不集聚也不空缺

    λBii∑Bik=0λBi-kiaki≤ω且aBii∑Bik=0λBi-kiaki≤ω,當(dāng)aiλi≠1時(shí)1Bi+1≤ω,????????????????????????????????????????????? 當(dāng)aiλi=1時(shí)。(18)

    證明過(guò)程與命題3類(lèi)似,在此不再贅述。

    值得注意的是,為了避免空滿(mǎn)樁現(xiàn)象,有樁自行車(chē)共享系統(tǒng)均須開(kāi)展自行車(chē)的再分布調(diào)度,以平衡各個(gè)站點(diǎn)的自行車(chē)數(shù)量。命題3及其推論1能夠揭示有樁自行車(chē)共享系統(tǒng)中哪些站點(diǎn)發(fā)生自行車(chē)空缺和集聚。這將為自行車(chē)共享系統(tǒng)的站點(diǎn)選址、自行車(chē)規(guī)劃提供決策支持,也能為公共自行車(chē)再分布調(diào)度中的取送起訖點(diǎn)及路徑規(guī)劃、空滿(mǎn)樁站點(diǎn)自行車(chē)取送數(shù)量等提供關(guān)鍵支撐。

    2 案例應(yīng)用

    2.1 案例介紹

    本研究以蘇州獨(dú)墅湖科教創(chuàng)新區(qū)的自行車(chē)共享系統(tǒng)作為例,驗(yàn)證本文模型。該系統(tǒng)共有37個(gè)自行車(chē)站點(diǎn),停車(chē)樁總數(shù)為1 260個(gè),平均自行車(chē)總數(shù)為442個(gè),站點(diǎn)和停車(chē)樁的具體分布如圖2所示。系統(tǒng)運(yùn)營(yíng)商提供了該區(qū)域內(nèi)2019年10月12日24 h的2 078次自行車(chē)出行數(shù)據(jù),這些數(shù)據(jù)記錄了高峰時(shí)期用戶(hù)借還自行車(chē)的時(shí)間及站點(diǎn)、車(chē)樁及用戶(hù)ID號(hào)。利用該數(shù)據(jù)統(tǒng)計(jì)得到自行車(chē)共享系統(tǒng)平均有效吞吐率為86.58人次/h。由于站點(diǎn)空樁情況下會(huì)出現(xiàn)用戶(hù)流失,故利用該數(shù)據(jù)統(tǒng)計(jì)各站點(diǎn)的空樁時(shí)段,將非空樁期間的用戶(hù)到達(dá)各站點(diǎn)次數(shù)按照時(shí)間比例填補(bǔ)到空樁時(shí)段,將其作為各個(gè)站點(diǎn)實(shí)際用戶(hù)到達(dá)率。利用原始數(shù)據(jù)還統(tǒng)計(jì)了行程時(shí)間矩陣、路由矩陣和各站點(diǎn)的最近站點(diǎn)(最近站點(diǎn)見(jiàn)表1所示)。由于篇幅原因,本文不再列出行程時(shí)間矩陣和路由矩陣。

    2.2 結(jié)果分析

    2.2.1 系統(tǒng)吞吐率

    根據(jù)近似模型計(jì)算得到的系統(tǒng)有效吞吐率為88.36人次/h,與真實(shí)吞吐率相對(duì)誤差為2.04%。圖3給出了站點(diǎn)的平均用戶(hù)到達(dá)率和有效吞吐率。可以看到,只有4個(gè)站點(diǎn)(4號(hào)、5號(hào)、7號(hào)、34號(hào))的用戶(hù)需求完全得到了滿(mǎn)足,即站點(diǎn)的有效吞吐率與用戶(hù)平均到達(dá)率相等。其余站點(diǎn)的供給與用戶(hù)需求不匹配,其中,1號(hào)、9號(hào)、22號(hào)、26號(hào)、27號(hào)、37號(hào)站點(diǎn)供不應(yīng)求的情況極為顯著,平均每小時(shí)有5名以上用戶(hù)不能借到自行車(chē)。這說(shuō)明目前研究區(qū)域內(nèi)有樁自行車(chē)共享系統(tǒng)的自行車(chē)分布極不平衡,系統(tǒng)內(nèi)自行車(chē)自循環(huán)結(jié)果在很大程度上無(wú)法滿(mǎn)足用戶(hù)的需求。

    圖4給出了系統(tǒng)有效吞吐速率與自行車(chē)投放量之間的關(guān)系。由圖可知,隨著自行車(chē)投放量的增加,系統(tǒng)有效吞吐率呈現(xiàn)階梯式的遞增趨勢(shì),當(dāng)自行車(chē)投放量大于225輛后,系統(tǒng)有效吞吐率保持不變,驗(yàn)證了性質(zhì)1。根據(jù)定義1可知,該有樁自行車(chē)共享系統(tǒng)的最優(yōu)自行車(chē)投放量為225輛,此時(shí)系統(tǒng)有效吞吐率達(dá)到最大值88.4人次/h。但是該系統(tǒng)實(shí)際的自行車(chē)投放量為442輛,故系統(tǒng)自行車(chē)數(shù)量明顯過(guò)剩。

    2.2.2 空滿(mǎn)樁t分析

    一般來(lái)說(shuō),如果自行車(chē)共享系統(tǒng)具有較好的服務(wù)水平,當(dāng)站點(diǎn)處于還車(chē)高峰時(shí),站點(diǎn)仍需保留20%的空樁;反之,當(dāng)站點(diǎn)處于借車(chē)高峰時(shí),站點(diǎn)也需要配備20%的車(chē)樁應(yīng)有的自行車(chē)數(shù)[18]。因此,結(jié)合空滿(mǎn)樁的定義,本文將定義2的比例參數(shù)設(shè)置為φ1=0.2,φ2=0.8,概率參數(shù)ω1=ω2=0.8。

    圖5給出了站點(diǎn)自行車(chē)庫(kù)存的平均值以及平衡自行車(chē)數(shù)的區(qū)間??梢钥闯?,絕大多數(shù)站點(diǎn)自行車(chē)庫(kù)存極為緊缺,少量站點(diǎn)庫(kù)存爆滿(mǎn),還有極少數(shù)的站點(diǎn)的平均庫(kù)存處于平衡區(qū)間內(nèi)。圖6為自行車(chē)投放量442輛下各站點(diǎn)自行車(chē)平均停留時(shí)間??梢钥闯?,自行車(chē)在站點(diǎn)4、5、7、34的平均停留時(shí)間超過(guò)10 h,在其他站點(diǎn)的停留時(shí)間在0.1~2.0 h范圍內(nèi)。上述進(jìn)一步驗(yàn)證了自行車(chē)在站點(diǎn)4、5、7、34的嚴(yán)重集聚狀態(tài)以及其他站點(diǎn)自行車(chē)庫(kù)存普遍緊缺情況。

    自行車(chē)投放量的增加會(huì)改變自行車(chē)集聚與空缺站點(diǎn)的空間分布。利用命題3計(jì)算得到圖4中a、b、c、d、e、f等6個(gè)點(diǎn)對(duì)應(yīng)的自行車(chē)投放量下,站點(diǎn)的集聚和空缺概率,具體的空間分布如圖7所示。可以看出,隨著自行車(chē)投放量的增加,站點(diǎn)4、5、7、34從空缺站點(diǎn)逐漸轉(zhuǎn)化為集聚站點(diǎn),站點(diǎn)29、36、39從空缺站點(diǎn)轉(zhuǎn)化為自行車(chē)既不集聚也不空缺站點(diǎn),其他站點(diǎn)始終處于自行車(chē)空缺狀態(tài)。隨著自行車(chē)投放量的增長(zhǎng)超過(guò)225輛,自行車(chē)集聚與空缺站點(diǎn)分布基本不再變化,投入的大量的自行車(chē)將迫使大量用戶(hù)在最近站點(diǎn)間騎行,試圖尋找最近站點(diǎn)還車(chē)。

    3 總結(jié)

    本文基于封閉排隊(duì)網(wǎng)絡(luò)構(gòu)造了有樁自行車(chē)共享系統(tǒng)吞吐率近似算法,考察了系統(tǒng)的運(yùn)營(yíng)性質(zhì)。理論上證明了有樁自行車(chē)共享系統(tǒng)吞吐率隨自行車(chē)投放量的增加呈非遞減性,給出了站點(diǎn)產(chǎn)生自行車(chē)集聚與空缺的判斷依據(jù)。以蘇州市為例進(jìn)行了上述算法的應(yīng)用與結(jié)論的驗(yàn)證,給出了該有樁自行車(chē)共享系統(tǒng)的服務(wù)績(jī)效。研究表明,隨著自行車(chē)投放量增長(zhǎng),吞吐率呈階梯式遞增且達(dá)到最大值后保持不變;自行車(chē)投放量一旦超過(guò)最優(yōu)數(shù)量將產(chǎn)生閑置,并且自行車(chē)集聚與空缺站點(diǎn)分布也將固定。

    依托本文建立的近似模型,可進(jìn)一步從多個(gè)方面進(jìn)行拓展,例如:本研究考慮的用戶(hù)到達(dá)自行車(chē)站點(diǎn)借車(chē)為平穩(wěn)隨機(jī)過(guò)程,未來(lái)可以考慮用戶(hù)到達(dá)存在早晚高峰的非平穩(wěn)過(guò)程;本研究中自行車(chē)最優(yōu)投放量并未考慮車(chē)站選址、車(chē)站間距、用戶(hù)騎行距離等因素,未來(lái)應(yīng)考慮這些因素對(duì)有樁自行車(chē)共享系統(tǒng)的最優(yōu)規(guī)劃設(shè)計(jì)問(wèn)題。

    參考文獻(xiàn):

    [1]ZHANG J, MENG M, WONG Y D, et al. A data-driven dynamic repositioning model in bicycle-sharing systems[J]. International Journal of Production Economics, 2021, 231: 107909. DOI: 10.1016/j.ijpe.2020.107909.

    [2]MEDDINR, DEMIAO P. The bike-sharing world map [EB/OL]. [2022-12-31].https://www.metrobike.net/.

    [3]LIN J J, YU C J. A bikeway network design model for urban areas[J]. Transportation, 2013, 40(1): 45-68. DOI: 10.1007/s11116-012-9409-6.

    [4]MIX R, HURTUBIA R, RAVEAU S. Optimal location of bike-sharing stations: a built environment and accessibility approach[J]. Transportation Research Part A: Policy and Practice, 2022, 160: 126-142. DOI: 10.1016/j.tra.2022.03.022.

    [5]FRICKER C, GAST N. Incentives and redistribution in homogeneous bike-sharing systems with stations of finite capacity[J]. EURO Journal on Transportation and Logistics, 2016, 5(3): 261-291. DOI: 10.1007/s13676-014-0053-5.

    [6]NEUMANN-SAAVEDRA B A, MATTFELD D C, HEWITT M. Assessing the operational impact of tactical planning models for bike-sharing redistribution[J]. Transportation Research Part A: Policy and Practice, 2021, 150: 216-235. DOI: 10.1016/j.tra.2021.06.003.

    [7]LV C, ZHANG C Y, LIAN K L, et al. A two-echelon fuzzyclustering based heuristic for large-scale bike sharing repositioning problem[J]. Transportation Research Part B: Methodological, 2022, 160: 54-75. DOI: 10.1016/j.trb.2022.04.003.

    [8]FU C Y, MA S F, ZHU N, et al. Bike-sharing inventory management for market expansion[J]. Transportation Research Part B: Methodological, 2022, 162: 28-54. DOI: 10.1016/j.trb.2022.05.009.

    [9]GAMMELLI D, WANG Y H, PRAK D, et al. Predictive and prescriptive performance of bike-sharing demand forecasts for inventory management[J]. Transportation Research Part C: Emerging Technologies, 2022, 138: 103571. DOI: 10.1016/j.trc.2022.103571.

    [10]LUO X L, GU W H, FAN W B. Joint design of shared-bike and transit services in corridors[J]. Transportation Research Part C: Emerging Technologies, 2021, 132: 103366. DOI: 10.1016/j.trc.2021.103366.

    [11]GEORGE D K, XIA C H. Fleet-sizing and service availability for a vehicle rental system via closed queueing networks[J]. European Journal of Operational Research, 2011, 211(1): 198-207. DOI: 10.1016/j.ejor.2010.12.015.

    [12]ELEBI D, YRSN A, IK H. Bicycle sharing system design with capacity allocations[J]. Transportation Research Part B: Methodological, 2018, 114: 86-98. DOI: 10.1016/j.trb.2018.05.018.

    [13]MALEKI VISHKAEI B, MAHDAVI I, MAHDAVI-AMIRI N, et al. Balancing public bicycle sharing system using inventory critical levels in queuing network[J]. Computers & Industrial Engineering, 2020, 141: 106277. DOI: 10.1016/j.cie.2020.106277.

    [14]KASPI M, RAVIV T, TZUR M. Bike-sharing systems: User dissatisfaction in the presence of unusable bicycles[J]. IISE Transactions, 2017, 49(2): 144-158. DOI: 10.1080/0740817x.2016.1224960.

    [15]CAGGIANI L, CAMPOREALE R, MARINELLI M, et al. User satisfaction based model for resource allocation in bike-sharing systems[J]. Transport Policy, 2019, 80: 117-126. DOI: 10.1016/j.tranpol.2018.03.003.

    [16]ROSS S M. Introduction to probability models[M]. 10th ed. Amsterdam: Elsevier Academic Press, 2010.

    [17]REISER M, LAVENBERG SS. Mean-value analysis of closed multichain queuing networks[J]. Journal of the ACM, 1980, 27(2): 313-322. DOI: 10.1145/322186.322195.

    [18]蔣琳. 城市公共自行車(chē)系統(tǒng)調(diào)度模型及算法研究[D].蘭州:蘭州交通大學(xué), 2017.

    猜你喜歡
    運(yùn)營(yíng)效率
    我國(guó)零售業(yè)上市公司運(yùn)營(yíng)效率實(shí)證分析
    基于DEA模型的農(nóng)業(yè)科技園區(qū)運(yùn)營(yíng)效率評(píng)價(jià)
    我國(guó)商業(yè)銀行運(yùn)營(yíng)效率的財(cái)務(wù)指標(biāo)分析
    提升校園眾創(chuàng)空間建設(shè)運(yùn)營(yíng)效率的對(duì)策建議
    商情(2017年5期)2017-03-30 23:19:27
    中通赴美上市探究
    時(shí)空視角下省級(jí)政府政務(wù)微博運(yùn)營(yíng)效率評(píng)價(jià)
    關(guān)于企業(yè)精益化管理的研究與思考
    提高國(guó)有資本運(yùn)營(yíng)效率對(duì)策及探究
    內(nèi)部控制與慈善組織運(yùn)營(yíng)效率
    事務(wù)所組織形式、運(yùn)營(yíng)效率與審計(jì)質(zhì)量
    91久久精品电影网| 国产精品电影一区二区三区| 波多野结衣高清无吗| 国产成人精品一,二区 | 国产精品电影一区二区三区| 高清毛片免费看| 国产一区二区三区在线臀色熟女| 国产伦精品一区二区三区视频9| ponron亚洲| 老熟妇乱子伦视频在线观看| 国产激情偷乱视频一区二区| 国产久久久一区二区三区| 好男人在线观看高清免费视频| or卡值多少钱| 一区福利在线观看| 欧美+亚洲+日韩+国产| 一个人看的www免费观看视频| 中文字幕久久专区| 精品日产1卡2卡| 午夜激情欧美在线| 18禁裸乳无遮挡免费网站照片| 色哟哟·www| 18禁黄网站禁片免费观看直播| 欧美三级亚洲精品| 人妻制服诱惑在线中文字幕| 亚洲国产精品国产精品| 国产麻豆成人av免费视频| 亚洲精品影视一区二区三区av| 亚洲国产欧洲综合997久久,| 男人舔女人下体高潮全视频| 久久精品人妻少妇| 插逼视频在线观看| 日韩,欧美,国产一区二区三区 | 精品人妻一区二区三区麻豆| 热99re8久久精品国产| 欧美色视频一区免费| 毛片女人毛片| 久久韩国三级中文字幕| 精品免费久久久久久久清纯| 欧美性猛交╳xxx乱大交人| 成人特级av手机在线观看| 国产成人精品婷婷| 高清日韩中文字幕在线| 男人舔女人下体高潮全视频| 五月伊人婷婷丁香| 国产精品福利在线免费观看| 只有这里有精品99| 村上凉子中文字幕在线| 国产精品久久久久久亚洲av鲁大| 亚洲第一区二区三区不卡| 亚洲经典国产精华液单| 不卡一级毛片| 熟女电影av网| 伊人久久精品亚洲午夜| 欧美日韩国产亚洲二区| 少妇的逼水好多| 国产综合懂色| 色尼玛亚洲综合影院| 成年av动漫网址| 99久久成人亚洲精品观看| 国产一级毛片在线| 精品99又大又爽又粗少妇毛片| 国产亚洲欧美98| 日本黄色片子视频| 色尼玛亚洲综合影院| 观看免费一级毛片| 国产精品野战在线观看| 不卡视频在线观看欧美| 成人亚洲欧美一区二区av| 国产免费一级a男人的天堂| 国产成人精品久久久久久| 悠悠久久av| 国产极品天堂在线| 日韩欧美精品v在线| 精品国产三级普通话版| 亚洲综合色惰| 国产精品一及| 亚州av有码| 久久久久久久久大av| 美女黄网站色视频| 国产精品,欧美在线| 欧美精品国产亚洲| 亚洲内射少妇av| 长腿黑丝高跟| 波多野结衣高清作品| 国产成年人精品一区二区| 色综合站精品国产| 中文精品一卡2卡3卡4更新| 国产一区二区激情短视频| 麻豆国产97在线/欧美| 亚洲av成人av| 变态另类丝袜制服| 哪里可以看免费的av片| 亚洲最大成人av| 国产不卡一卡二| 亚洲精品久久久久久婷婷小说 | 午夜激情欧美在线| 一本一本综合久久| 在线观看午夜福利视频| 欧美zozozo另类| 欧美高清成人免费视频www| 女的被弄到高潮叫床怎么办| 色视频www国产| 国产久久久一区二区三区| 亚洲第一电影网av| 国产又黄又爽又无遮挡在线| 又黄又爽又刺激的免费视频.| 99热这里只有精品一区| 极品教师在线视频| 校园春色视频在线观看| 国产精品伦人一区二区| 男人狂女人下面高潮的视频| 精品人妻一区二区三区麻豆| 91狼人影院| 亚洲经典国产精华液单| 人妻夜夜爽99麻豆av| 美女大奶头视频| 国产一区二区激情短视频| 欧美色视频一区免费| 国产久久久一区二区三区| 一级毛片aaaaaa免费看小| 国产激情偷乱视频一区二区| 欧美变态另类bdsm刘玥| 久久久久久久久久成人| 亚洲成av人片在线播放无| 中国美白少妇内射xxxbb| 国产精品国产高清国产av| 夜夜看夜夜爽夜夜摸| 欧美一区二区精品小视频在线| 男女啪啪激烈高潮av片| 此物有八面人人有两片| 亚洲av中文av极速乱| 久久久午夜欧美精品| av天堂在线播放| 亚洲国产精品sss在线观看| 成人国产麻豆网| 美女国产视频在线观看| 亚洲国产欧洲综合997久久,| 久久久久久久久久黄片| 亚洲国产欧洲综合997久久,| 少妇的逼好多水| 长腿黑丝高跟| 狂野欧美激情性xxxx在线观看| 日韩强制内射视频| 欧美zozozo另类| 99riav亚洲国产免费| 精品人妻熟女av久视频| 日本撒尿小便嘘嘘汇集6| 亚洲国产精品久久男人天堂| 人妻系列 视频| 九九爱精品视频在线观看| 精品久久久噜噜| 国产在线男女| 亚洲欧美日韩东京热| 国产黄片美女视频| 国产成人freesex在线| 国产探花极品一区二区| 久久久久久久久久成人| 亚洲成人av在线免费| 麻豆久久精品国产亚洲av| 国产精品1区2区在线观看.| 亚洲电影在线观看av| 日韩人妻高清精品专区| 天天躁夜夜躁狠狠久久av| 亚洲精品亚洲一区二区| 亚洲,欧美,日韩| 九草在线视频观看| av.在线天堂| 日韩精品青青久久久久久| 欧美日本亚洲视频在线播放| 草草在线视频免费看| 国产高清视频在线观看网站| 一级毛片久久久久久久久女| 欧美色视频一区免费| 人妻系列 视频| 亚洲国产精品久久男人天堂| 欧美潮喷喷水| 99久久九九国产精品国产免费| 国产黄片美女视频| 夜夜看夜夜爽夜夜摸| 大型黄色视频在线免费观看| 亚洲av中文av极速乱| 亚洲精品久久久久久婷婷小说 | 女的被弄到高潮叫床怎么办| 最近最新中文字幕大全电影3| 国产成人精品一,二区 | 成人国产麻豆网| 亚洲真实伦在线观看| 人妻夜夜爽99麻豆av| .国产精品久久| 麻豆成人av视频| h日本视频在线播放| 在线观看66精品国产| 91久久精品电影网| 深夜a级毛片| 久久久久久久久久黄片| 91久久精品国产一区二区三区| av在线天堂中文字幕| av视频在线观看入口| 午夜精品在线福利| 人妻系列 视频| 亚洲图色成人| 午夜精品一区二区三区免费看| 国产欧美日韩精品一区二区| 国产精品乱码一区二三区的特点| 欧美一区二区亚洲| 色播亚洲综合网| 国产精品久久久久久精品电影小说 | 六月丁香七月| 精品久久久久久久末码| 赤兔流量卡办理| 人人妻人人看人人澡| 黄色一级大片看看| 天堂中文最新版在线下载 | 亚洲欧美日韩高清在线视频| 尾随美女入室| 人人妻人人澡人人爽人人夜夜 | 国产亚洲91精品色在线| 亚洲七黄色美女视频| 一本一本综合久久| 99热这里只有精品一区| 大又大粗又爽又黄少妇毛片口| av在线亚洲专区| 一进一出抽搐gif免费好疼| 日韩欧美精品v在线| 99热只有精品国产| 偷拍熟女少妇极品色| 性插视频无遮挡在线免费观看| 国产精品久久久久久精品电影小说 | 亚洲av中文字字幕乱码综合| 国产伦一二天堂av在线观看| 成人毛片a级毛片在线播放| 91aial.com中文字幕在线观看| 波多野结衣高清无吗| 色吧在线观看| 国产高清有码在线观看视频| 干丝袜人妻中文字幕| 国产伦精品一区二区三区四那| 久久精品国产亚洲av涩爱 | 免费搜索国产男女视频| 99在线视频只有这里精品首页| 别揉我奶头 嗯啊视频| 国产午夜精品一二区理论片| 精品熟女少妇av免费看| 成人三级黄色视频| АⅤ资源中文在线天堂| 欧美激情在线99| 精品人妻熟女av久视频| 欧洲精品卡2卡3卡4卡5卡区| 天堂av国产一区二区熟女人妻| 成人亚洲欧美一区二区av| 高清日韩中文字幕在线| 最近手机中文字幕大全| 国产又黄又爽又无遮挡在线| 久久精品国产亚洲av涩爱 | 观看美女的网站| 亚洲va在线va天堂va国产| 一级av片app| 一级毛片aaaaaa免费看小| 免费av不卡在线播放| 久久99精品国语久久久| 亚洲av成人av| videossex国产| 中出人妻视频一区二区| 亚洲丝袜综合中文字幕| 国产熟女欧美一区二区| 黑人高潮一二区| 国产人妻一区二区三区在| 免费看光身美女| 久久婷婷人人爽人人干人人爱| 亚洲最大成人手机在线| 亚洲人成网站高清观看| 麻豆久久精品国产亚洲av| 97超视频在线观看视频| 国产麻豆成人av免费视频| 色综合色国产| 亚洲欧洲日产国产| 日本黄色片子视频| 精品熟女少妇av免费看| 国产黄色小视频在线观看| 亚洲在线观看片| 国产高潮美女av| 国产黄色小视频在线观看| 亚洲av成人精品一区久久| 国产av在哪里看| 18禁在线无遮挡免费观看视频| 成人二区视频| 看十八女毛片水多多多| 高清日韩中文字幕在线| 狠狠狠狠99中文字幕| 99热这里只有是精品50| 亚洲乱码一区二区免费版| 18禁在线无遮挡免费观看视频| 午夜精品一区二区三区免费看| 国产精品无大码| 18禁裸乳无遮挡免费网站照片| 日韩在线高清观看一区二区三区| 精品一区二区免费观看| 一进一出抽搐动态| 伊人久久精品亚洲午夜| 久久精品国产清高在天天线| 国产真实乱freesex| 久久久久久久久久成人| 美女国产视频在线观看| 麻豆成人午夜福利视频| 中文资源天堂在线| 18禁黄网站禁片免费观看直播| 春色校园在线视频观看| av福利片在线观看| 国产精品日韩av在线免费观看| 国产精品久久久久久久久免| 亚洲色图av天堂| 亚洲真实伦在线观看| 精品人妻偷拍中文字幕| 欧美性感艳星| 亚洲欧美日韩高清在线视频| 亚洲精品乱码久久久久久按摩| 亚洲七黄色美女视频| 中文字幕人妻熟人妻熟丝袜美| 久久欧美精品欧美久久欧美| 一进一出抽搐动态| 亚洲激情五月婷婷啪啪| 午夜激情福利司机影院| 亚洲av中文av极速乱| 精品人妻熟女av久视频| 男女啪啪激烈高潮av片| 少妇人妻精品综合一区二区 | 久久人人爽人人爽人人片va| 国产色婷婷99| 免费不卡的大黄色大毛片视频在线观看 | 18禁在线播放成人免费| 欧美高清成人免费视频www| 成年女人永久免费观看视频| 欧美变态另类bdsm刘玥| 国产精品爽爽va在线观看网站| a级一级毛片免费在线观看| 少妇的逼好多水| 久久久精品大字幕| 九色成人免费人妻av| 亚洲欧美日韩高清在线视频| 哪里可以看免费的av片| 日韩大尺度精品在线看网址| 亚洲一级一片aⅴ在线观看| 精品人妻一区二区三区麻豆| 在线免费观看的www视频| 乱码一卡2卡4卡精品| 极品人妻少妇av视频| 久久久久久久国产电影| 九九久久精品国产亚洲av麻豆| 99热这里只有精品一区| 男的添女的下面高潮视频| 久久青草综合色| 丝袜在线中文字幕| 曰老女人黄片| 少妇被粗大的猛进出69影院 | 日日啪夜夜爽| 男女啪啪激烈高潮av片| 一区二区三区四区激情视频| 中文字幕亚洲精品专区| 在线天堂最新版资源| 亚洲内射少妇av| 亚洲成人一二三区av| 亚洲人成77777在线视频| 国产男女内射视频| 久久久久久人妻| 肉色欧美久久久久久久蜜桃| 日韩免费高清中文字幕av| 高清午夜精品一区二区三区| 日韩亚洲欧美综合| 纵有疾风起免费观看全集完整版| 九色亚洲精品在线播放| 国产精品国产三级专区第一集| 九色亚洲精品在线播放| 亚洲av免费高清在线观看| 美女脱内裤让男人舔精品视频| 午夜91福利影院| 久久久久久久久久人人人人人人| 我要看黄色一级片免费的| 日本免费在线观看一区| 女性被躁到高潮视频| 99国产综合亚洲精品| 日韩欧美精品免费久久| 我的老师免费观看完整版| 亚洲欧美中文字幕日韩二区| 久久人妻熟女aⅴ| 亚洲av免费高清在线观看| 日本免费在线观看一区| 99九九线精品视频在线观看视频| 精品亚洲乱码少妇综合久久| 色网站视频免费| 亚洲av成人精品一区久久| 看免费成人av毛片| 国产亚洲最大av| 日本91视频免费播放| 男人操女人黄网站| 99热6这里只有精品| 纯流量卡能插随身wifi吗| 久久久国产精品麻豆| 日本黄大片高清| 国产综合精华液| 亚洲精华国产精华液的使用体验| 日韩不卡一区二区三区视频在线| 99热6这里只有精品| 亚洲成人av在线免费| 国产av精品麻豆| 国产精品女同一区二区软件| 精品一区二区三卡| 极品人妻少妇av视频| 啦啦啦中文免费视频观看日本| 在线 av 中文字幕| 精品一品国产午夜福利视频| 五月玫瑰六月丁香| 亚洲不卡免费看| 日本色播在线视频| 久久精品夜色国产| 亚洲欧美日韩另类电影网站| 涩涩av久久男人的天堂| 天天操日日干夜夜撸| 亚洲国产最新在线播放| 伦理电影大哥的女人| 亚洲,一卡二卡三卡| 国产在线视频一区二区| 国产综合精华液| 蜜桃在线观看..| 国国产精品蜜臀av免费| 国产淫语在线视频| 免费久久久久久久精品成人欧美视频 | 老司机影院毛片| 欧美性感艳星| 久久久久久久亚洲中文字幕| 精品一区二区免费观看| 中文字幕制服av| 黄片无遮挡物在线观看| a级毛色黄片| 久久久午夜欧美精品| a级毛片黄视频| 国产成人精品无人区| 啦啦啦中文免费视频观看日本| 国产淫语在线视频| 成人国产av品久久久| 免费黄色在线免费观看| 在线天堂最新版资源| 大片免费播放器 马上看| 女人久久www免费人成看片| 久久影院123| 亚洲欧美成人精品一区二区| 成人国语在线视频| 91精品国产九色| 街头女战士在线观看网站| 久久精品国产亚洲网站| freevideosex欧美| 黄色视频在线播放观看不卡| 男女免费视频国产| 啦啦啦视频在线资源免费观看| 人人妻人人爽人人添夜夜欢视频| 青春草亚洲视频在线观看| 另类亚洲欧美激情| 亚洲精品国产色婷婷电影| 91精品一卡2卡3卡4卡| 久久午夜综合久久蜜桃| 国产av一区二区精品久久| 免费久久久久久久精品成人欧美视频 | 99久久综合免费| 国产成人精品在线电影| 亚洲色图综合在线观看| 欧美激情 高清一区二区三区| 欧美少妇被猛烈插入视频| 美女内射精品一级片tv| 内地一区二区视频在线| 一区在线观看完整版| 免费看av在线观看网站| 黄片无遮挡物在线观看| 丝袜脚勾引网站| 免费av不卡在线播放| av国产精品久久久久影院| 日韩在线高清观看一区二区三区| 亚洲精品国产色婷婷电影| h视频一区二区三区| 亚州av有码| 午夜精品国产一区二区电影| 亚洲av中文av极速乱| 国产 精品1| 久久99一区二区三区| 色婷婷久久久亚洲欧美| 天堂中文最新版在线下载| 精品一区二区三区视频在线| 亚洲av日韩在线播放| 日本av手机在线免费观看| kizo精华| 欧美国产精品一级二级三级| 日韩电影二区| 五月玫瑰六月丁香| 国产精品国产av在线观看| 免费观看无遮挡的男女| 免费大片黄手机在线观看| av在线播放精品| 精品人妻熟女毛片av久久网站| 精品久久国产蜜桃| 国产精品久久久久久精品古装| 亚洲精品亚洲一区二区| 亚洲国产精品成人久久小说| 永久免费av网站大全| 国产精品嫩草影院av在线观看| 五月开心婷婷网| 免费观看性生交大片5| 在线精品无人区一区二区三| 亚洲婷婷狠狠爱综合网| 交换朋友夫妻互换小说| 三上悠亚av全集在线观看| 精品久久蜜臀av无| 99国产精品免费福利视频| 婷婷色综合www| 国产极品天堂在线| 新久久久久国产一级毛片| 丰满乱子伦码专区| 久久精品人人爽人人爽视色| 狂野欧美激情性bbbbbb| 成人国产麻豆网| 秋霞伦理黄片| 91精品伊人久久大香线蕉| av有码第一页| 国产精品嫩草影院av在线观看| 如日韩欧美国产精品一区二区三区 | 另类亚洲欧美激情| 精品国产国语对白av| 亚洲激情五月婷婷啪啪| 久久精品国产鲁丝片午夜精品| 丰满乱子伦码专区| 亚洲人成77777在线视频| 午夜激情福利司机影院| 日本vs欧美在线观看视频| 国产不卡av网站在线观看| 精品国产国语对白av| freevideosex欧美| videos熟女内射| 成人漫画全彩无遮挡| 国产熟女欧美一区二区| 青春草国产在线视频| 日本vs欧美在线观看视频| 亚洲精华国产精华液的使用体验| 午夜老司机福利剧场| a级毛片在线看网站| 免费黄色在线免费观看| 国产午夜精品久久久久久一区二区三区| 午夜福利视频精品| 国产精品蜜桃在线观看| 五月伊人婷婷丁香| 狠狠精品人妻久久久久久综合| 99热国产这里只有精品6| 亚洲精品一区蜜桃| 亚洲精品乱久久久久久| 国产色婷婷99| 黑人高潮一二区| 亚洲av成人精品一区久久| 校园人妻丝袜中文字幕| .国产精品久久| 国产精品熟女久久久久浪| 天天躁夜夜躁狠狠久久av| 免费黄色在线免费观看| 黄片无遮挡物在线观看| 亚洲人成网站在线播| 久久久久久伊人网av| 伊人久久国产一区二区| 午夜激情久久久久久久| 日韩一区二区视频免费看| 丰满迷人的少妇在线观看| 97超碰精品成人国产| 美女福利国产在线| 亚洲精品日韩在线中文字幕| 午夜福利视频在线观看免费| 综合色丁香网| 交换朋友夫妻互换小说| 国产深夜福利视频在线观看| 成人18禁高潮啪啪吃奶动态图 | 777米奇影视久久| 另类亚洲欧美激情| 一区二区三区四区激情视频| 亚洲成人av在线免费| 亚洲欧美中文字幕日韩二区| 国产精品国产三级专区第一集| 美女视频免费永久观看网站| 久久久久久久久大av| 亚洲国产毛片av蜜桃av| 精品一区在线观看国产| 国产一区二区三区综合在线观看 | 在现免费观看毛片| 久久久久久久久久久久大奶| 久久精品国产亚洲网站| 国产精品成人在线| 特大巨黑吊av在线直播| 亚洲怡红院男人天堂| 国产精品麻豆人妻色哟哟久久| 国产片内射在线| 久久久久国产网址| av电影中文网址| 久久热精品热| 丝袜在线中文字幕| 亚洲国产av影院在线观看| 91久久精品国产一区二区三区| 男女无遮挡免费网站观看| 日韩强制内射视频| 国产精品99久久久久久久久| 观看美女的网站| 性高湖久久久久久久久免费观看| 日韩亚洲欧美综合| 欧美日韩视频精品一区| 乱人伦中国视频| 成人国产av品久久久| 99国产综合亚洲精品| 香蕉精品网在线| 亚洲欧洲精品一区二区精品久久久 | 国产亚洲av片在线观看秒播厂| 有码 亚洲区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 人成视频在线观看免费观看|