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

    青藏高原雪災(zāi)高風(fēng)險(xiǎn)區(qū)飼草料儲備庫選址分析

    2020-08-12 14:23:46趙富昌劉峰貴陳永萍支澤民
    關(guān)鍵詞:雪災(zāi)儲備庫飼草料

    趙富昌,周 強(qiáng),2,劉峰貴,2,陳 瓊,2,陳永萍,支澤民

    青藏高原雪災(zāi)高風(fēng)險(xiǎn)區(qū)飼草料儲備庫選址分析

    趙富昌1,周 強(qiáng)1,2※,劉峰貴1,2,陳 瓊1,2,陳永萍1,支澤民1

    (1. 青海師范大學(xué)地理科學(xué)學(xué)院,西寧 810008;2. 高原科學(xué)與可持續(xù)發(fā)展研究院,西寧 810008)

    為了更好的發(fā)揮飼草料儲備庫在災(zāi)害風(fēng)險(xiǎn)管理中的應(yīng)用,提高高原牧區(qū)災(zāi)害韌性。該研究將雪災(zāi)高風(fēng)險(xiǎn)區(qū)與實(shí)際多年雪災(zāi)發(fā)生區(qū)相結(jié)合,采用-中值模型提出了飼草料儲備庫的選址分析方法,并利用貪婪取走啟發(fā)式算法對模型求解。以青海省雪災(zāi)高風(fēng)險(xiǎn)區(qū)為例,選址分析得出治多、瑪沁、都蘭3座飼草料儲備庫和麻多、扎朵、秋智等13座飼草料儲備點(diǎn)的分布位置與儲備規(guī)模。其中治多儲備規(guī)模為21.0萬t;瑪沁儲備規(guī)模為13.2萬 t;都蘭儲備規(guī)模為13.5萬t。布局上形成了以儲備庫為中心,儲備點(diǎn)為外圍的二元結(jié)構(gòu),驗(yàn)證了分析方法的合理性和可行性。從而為高原牧區(qū)飼草料儲備庫的選址提供參考性理論依據(jù)和可操作性方法。

    風(fēng)險(xiǎn)評估;貯備飼料;災(zāi)害;儲備庫;選址;-中值模型;雪災(zāi)高風(fēng)險(xiǎn)區(qū)

    0 引 言

    全球氣候變化背景下極端自然災(zāi)害發(fā)生的頻率、強(qiáng)度和持續(xù)時(shí)間呈現(xiàn)增加態(tài)勢,人類遭受自然災(zāi)害風(fēng)險(xiǎn)的機(jī)率在進(jìn)一步擴(kuò)大[1-2]。防災(zāi)減災(zāi)已經(jīng)成為人與自然和諧共生的永恒課題[3]。聯(lián)合國減少災(zāi)害風(fēng)險(xiǎn)辦公室(United Nations Office for Disaster Risk Reduction)在第五屆全球減災(zāi)平臺大會上表明減輕災(zāi)害風(fēng)險(xiǎn)已成為國際社會面臨的重大共同挑戰(zhàn),提出了理解災(zāi)害風(fēng)險(xiǎn),加強(qiáng)災(zāi)害風(fēng)險(xiǎn)管理等優(yōu)先行動領(lǐng)域[4]。近年來,國內(nèi)外學(xué)者通過對單災(zāi)種、災(zāi)害鏈及災(zāi)害遭遇的災(zāi)害風(fēng)險(xiǎn)研究[5-7],促使自然災(zāi)害風(fēng)險(xiǎn)評估得到了極大的發(fā)展,但將災(zāi)害風(fēng)險(xiǎn)的研究成果上升到災(zāi)害風(fēng)險(xiǎn)管理高度,在基礎(chǔ)設(shè)施規(guī)劃與防災(zāi)備災(zāi)發(fā)展方面有效融入風(fēng)險(xiǎn)評價(jià)的思想,促進(jìn)災(zāi)害風(fēng)險(xiǎn)評價(jià)成果的有效轉(zhuǎn)化應(yīng)用還亟待加強(qiáng)[8]。

    雪災(zāi)是全球十大自然災(zāi)害之一,根據(jù)雪災(zāi)的形成條件、分布范圍和表現(xiàn)形式,雪災(zāi)類型普遍被分為雪崩、風(fēng)吹雪(風(fēng)雪流)和牧區(qū)雪災(zāi)[9]。其中牧區(qū)雪災(zāi)是由于積雪過厚和維持時(shí)間太久,以至于積雪長時(shí)間掩埋牧草,使牲畜無法正常采食,導(dǎo)致大量牲畜掉膘和死亡的自然災(zāi)害,也是中國發(fā)生頻繁、影響最為嚴(yán)重的一類雪災(zāi)[10-11]。青藏高原是全國三大雪災(zāi)高發(fā)區(qū)之一[12],發(fā)生時(shí)期一般是從當(dāng)年10月開始至次年4月結(jié)束。較嚴(yán)重的雪災(zāi)已呈周期性出現(xiàn),對于以游牧或半游牧為主,靠天養(yǎng)畜,抗御自然災(zāi)害能力較差的牧業(yè)生產(chǎn),牧區(qū)雪災(zāi)已經(jīng)嚴(yán)重制約牧業(yè)生產(chǎn)發(fā)展和民族地區(qū)的經(jīng)濟(jì)振興[13-15]。2018年12月到2019年3月,位于青藏高原的玉樹和果洛藏族自治州因降雪持續(xù)時(shí)間太長,積雪過厚,導(dǎo)致5.8萬頭牲畜因凍傷死亡。持續(xù)低溫和飼草料儲備不足引起的牲畜饑荒以及道路阻斷等狀況造成的直接經(jīng)濟(jì)損失達(dá)1.92億元人民幣。青藏高原雪災(zāi)損失程度還會隨著全球變暖、極端氣候事件的加劇愈發(fā)嚴(yán)重[16-17],牧區(qū)雪災(zāi)風(fēng)險(xiǎn)在未來將進(jìn)一步加大。在理解災(zāi)害風(fēng)險(xiǎn)的基礎(chǔ)上進(jìn)行災(zāi)害風(fēng)險(xiǎn)管理是防災(zāi)減災(zāi)的關(guān)鍵性所在[18]。飼草料儲備庫的科學(xué)選址、合理儲備是提高災(zāi)害風(fēng)險(xiǎn)管理水平的重要措施之一,關(guān)乎牧區(qū)經(jīng)濟(jì)可持續(xù)發(fā)展和社會穩(wěn)定。

    飼草料合理儲備和有效調(diào)度可以提高雪災(zāi)承災(zāi)體系統(tǒng)的剛性和韌性,在有效抵抗牧區(qū)雪災(zāi)中發(fā)揮著至關(guān)重要的作用。2019年第二次青藏高原綜合科學(xué)考察中,有關(guān)專家提出了加快構(gòu)建牲畜飼草料儲備與應(yīng)急調(diào)度體系的建議,促進(jìn)飼草種植,就地儲備的備災(zāi)體系建設(shè)建立是應(yīng)對牧區(qū)雪災(zāi)的有效措施[19-20]。因而飼草料儲備庫的合理規(guī)劃與科學(xué)選址成為提高雪災(zāi)應(yīng)急能力、降低受災(zāi)損失和保障應(yīng)急救援順利實(shí)施的關(guān)鍵[21-22]。目前,學(xué)者主要從加權(quán)服務(wù)質(zhì)量期望值最大[23]、服務(wù)成本最小[24]、覆蓋范圍最大為目標(biāo)構(gòu)建-中值、-中心、最大覆蓋等選址模型[25-26],并采用啟發(fā)式、粒子群及NSGA-Ⅱ等算法對模型進(jìn)行求解[27],從而對救災(zāi)物資儲備庫進(jìn)行選址優(yōu)化分析,但這些選址方法對災(zāi)害風(fēng)險(xiǎn)等級的考慮較少。雪災(zāi)風(fēng)險(xiǎn)管理應(yīng)該建立在雪災(zāi)風(fēng)險(xiǎn)評估的基礎(chǔ)上,才能做到備災(zāi)精準(zhǔn)和救災(zāi)及時(shí),但是目前相關(guān)方面的研究卻很少見。因此該研究嘗試?yán)醚?zāi)風(fēng)險(xiǎn)評估的研究結(jié)果,采用-中值模型,對青海省雪災(zāi)高風(fēng)險(xiǎn)區(qū)飼草料儲備庫的選址進(jìn)行分析,并利用貪婪取走啟發(fā)式算法對模型求解。應(yīng)用雪災(zāi)風(fēng)險(xiǎn)評估成果進(jìn)行飼草料儲備庫選址,從而達(dá)到為高原牧區(qū)雪災(zāi)風(fēng)險(xiǎn)管理服務(wù)的目的。

    1 飼草料儲備庫選址分析方法

    1.1 分析思路

    20世紀(jì)60年代以來,在選址領(lǐng)域中經(jīng)典的選址模型有-中心模型、最大覆蓋模型、-中值模型,這些模型的中心理論是設(shè)施點(diǎn)與需求點(diǎn)之間的距離[28]。-中心模型主要強(qiáng)調(diào)公平性,一旦選定應(yīng)急物資儲備庫,系統(tǒng)內(nèi)不論哪個(gè)需求點(diǎn)出現(xiàn)災(zāi)情,都能保證在最短時(shí)間內(nèi)得到物資需求滿足,但是在實(shí)際情況中,儲備庫為體現(xiàn)出物資供應(yīng)的公平性,可能為了照顧到邊緣區(qū)域稀疏需求點(diǎn)的零星物資需求,從而造成資源浪費(fèi)[29]。以覆蓋需求最大化為目標(biāo)的最大覆蓋模型,在應(yīng)急物資儲備庫的覆蓋范圍和需求點(diǎn)數(shù)量都已知的情況下,實(shí)現(xiàn)了儲備點(diǎn)對所有需求點(diǎn)的全覆蓋,但是沒有考慮各需求點(diǎn)的需求規(guī)模,可能出現(xiàn)需求點(diǎn)“低需求、高儲備”的問題,從而造成資源浪費(fèi)[30]。-中值模型主要強(qiáng)調(diào)的是時(shí)效性,通過加權(quán)距離最短實(shí)現(xiàn)物資的運(yùn)輸速度最快,在不考慮其他因素的情況下通過加權(quán)距離的方法直接反應(yīng)了需求量與距離的關(guān)系[31]。綜上所述,基于總救助時(shí)間最短的-中值模型更適于青藏高原飼草料儲備庫的選址。利用飼草料需求量對距離加權(quán),兼顧了應(yīng)急救助效率和覆蓋度,能較好地滿足牧區(qū)雪災(zāi)應(yīng)急救助的工作特點(diǎn)。

    -中值模型常用的求解算法有免疫算法、蟻群算法和貪婪取走啟發(fā)式算法等。免疫算法利用免疫系統(tǒng)的多樣性產(chǎn)生和維持機(jī)制來保持群體的多樣性,從而克服一般尋優(yōu)過程尤其是多峰值函數(shù)中最難應(yīng)付的“早熟”問題,最終求得全局最優(yōu)解。但是存在群體數(shù)量控制和收斂速度慢等問題[32];蟻群算法在求解性能上,具有很強(qiáng)的魯棒性和搜索較好解的能力。但是參數(shù)設(shè)置不當(dāng)容易盲目的隨機(jī)搜索導(dǎo)致搜索時(shí)間較長[33]。貪婪取走啟發(fā)式算法具有步驟簡單、計(jì)算量小、節(jié)約時(shí)間和易于進(jìn)行定性分析和定量研究的優(yōu)點(diǎn)。但是由于每次都找局部最優(yōu),因此算法的效果受初始位置影響很大[34]。-中值問題屬于NP-Hard問題,要尋求一種精確高效的算法求解最優(yōu)解的可能性不大,但是貪婪取走啟發(fā)式算法可以使其解同最優(yōu)解盡可能的接近,同時(shí)又保證較好的穩(wěn)定性。因此該研究通過結(jié)合實(shí)際情況進(jìn)行多次驗(yàn)證的方式弱化該算法的弊端,從而實(shí)現(xiàn)-中值模型的快速、準(zhǔn)確求解。

    該研究基于雪災(zāi)風(fēng)險(xiǎn)評估成果,對雪災(zāi)高風(fēng)險(xiǎn)區(qū)進(jìn)行儲備庫選址分析,分析思路如圖1所示。1)根據(jù)雪災(zāi)風(fēng)險(xiǎn)圖確定雪災(zāi)高風(fēng)險(xiǎn)區(qū)并篩選出高風(fēng)險(xiǎn)區(qū)內(nèi)的縣、鄉(xiāng)、村級政府駐地;2)選擇-中值模型;3)根據(jù)-中值模型要求,進(jìn)行最佳路徑分析和飼草料需求量計(jì)算;4)運(yùn)用LINGO軟件用貪婪取走啟發(fā)式算法對-中值模型進(jìn)行計(jì)算,求解出理論值內(nèi)的縣級飼草料儲備庫和鄉(xiāng)鎮(zhèn)級飼草料儲備點(diǎn);5)根據(jù)輻射面積最大和調(diào)運(yùn)時(shí)間最短為最佳選址的原則,結(jié)合牲畜密度、交通節(jié)點(diǎn)通達(dá)性、地形條件、草地生產(chǎn)力及飼草料基地實(shí)際建設(shè)情況等選址影響因素,確定出最佳飼草料儲備庫和儲備點(diǎn)的分布位置及儲備規(guī)模(飼草料需求量)。

    圖1 飼草料儲備庫選址分析思路

    1.2 模型選擇與求解

    -中值模型是指在已知需求點(diǎn)集合位置、需求數(shù)量和設(shè)施備選點(diǎn)集合位置的情況下,分別為個(gè)設(shè)施找到合適的位置并指派每個(gè)需求點(diǎn)到一個(gè)特定的設(shè)施,使之達(dá)到連接需求點(diǎn)和設(shè)施點(diǎn)的總加權(quán)出行成本(距離、時(shí)間、費(fèi)用等)最低,適用于空間效率導(dǎo)向[35]。-中值模型的數(shù)學(xué)表達(dá)式如下[36]:

    式中表示目標(biāo)函數(shù);為個(gè)需求點(diǎn)的集合,={1,2,…,};為個(gè)設(shè)施備選點(diǎn)的集合,{1,2,…,};為需求點(diǎn)的編號,代表需求點(diǎn)的特定位置;為設(shè)施備選點(diǎn)的編號,代表可能設(shè)施點(diǎn)的位置;為設(shè)施配置的數(shù)量;W為需求點(diǎn)的需求量;D為需求點(diǎn)與設(shè)施之間的出行成本(距離、時(shí)間、費(fèi)用等);X為0/1決策變量,取1表示位置有設(shè)施設(shè)置,取0表示位置無設(shè)施設(shè)置;Y為0/1決策變量,取1表示設(shè)施為需求點(diǎn)提供服務(wù),取0表示設(shè)施不為需求點(diǎn)提供服務(wù)。

    式(1)的目標(biāo)在于使需求點(diǎn)與最近設(shè)施之間總的加權(quán)出行成本最小化;式(2)表示研究區(qū)域內(nèi)需要配置個(gè)設(shè)施;式(3)表示每個(gè)需求點(diǎn)只允許被分配到一個(gè)設(shè)施;式(4)約束需求點(diǎn)只允許被分配到有設(shè)施存在的位置;式(5)表示在位置設(shè)置或不設(shè)設(shè)施;式(6)表示需求點(diǎn)與設(shè)施可能位置點(diǎn)之間只有被服務(wù)(取值為1)與未被服務(wù)(取值為0)兩種關(guān)系。

    1.2.1 飼草料需求量分析方法

    飼草料儲備庫的選址主要服務(wù)于以草地畜牧業(yè)為主要生產(chǎn)方式的牧業(yè)區(qū)域,需求量為牲畜在發(fā)生雪災(zāi)時(shí)需要供給的飼草料。一只羊一天的飼草料需求量為5 kg,按照雪災(zāi)風(fēng)險(xiǎn)圖的等級和雪災(zāi)等級指標(biāo)可以確定雪災(zāi)持續(xù)時(shí)間,小雪災(zāi)、中雪災(zāi)、大雪災(zāi)和特大雪災(zāi)的持續(xù)時(shí)間分別為≥5、≥7、≥10和>30 d[37]。牲畜飼草料需求量公式(7)如下:

    FR=5(7)

    式中FR為飼草料需求量,kg;為通過羊單位轉(zhuǎn)換的實(shí)際牲畜數(shù)量(一匹馬折算為6只羊單位,一頭牛折算為5只羊單位,一頭驢或騾子或鹿折算為4只羊單位,一只綿羊等于1只羊單位);為雪災(zāi)持續(xù)時(shí)間,d[16]。

    1.2.2 最佳路徑分析方法

    交通路網(wǎng)是飼草料儲備庫選址的重要因素,雪災(zāi)發(fā)生時(shí),快速及時(shí)的送達(dá)飼草料是牧戶最關(guān)注的。以道路數(shù)據(jù)為基礎(chǔ)進(jìn)行道路等級賦權(quán)。以需求點(diǎn)(鄉(xiāng)鎮(zhèn)級政府駐地/村級政府駐地)為出發(fā)點(diǎn),設(shè)施點(diǎn)(縣級政府駐地/鄉(xiāng)鎮(zhèn)級政府駐地)為目的點(diǎn),運(yùn)用ArcGIS軟件中的Network Analyst求解出基于道路的最佳路徑。

    1.2.3 應(yīng)用LINGO軟件進(jìn)行算法求解

    LINGO是由美國LINDO系統(tǒng)公司研制的一款主要解決優(yōu)化問題的軟件,語言通俗易懂、計(jì)算準(zhǔn)確、容易調(diào)整,利于模型的建立和求解應(yīng)用。

    貪婪取走啟發(fā)式算法常用于解決-中值問題,該算法是一種不追求最優(yōu)解,只希望得到較為滿意解的方法[34]。貪婪取走啟發(fā)式算法步驟如下:

    1)初始化,令當(dāng)前選中設(shè)施點(diǎn)=,即將所有備選設(shè)施點(diǎn)都選中;2)將每個(gè)需求點(diǎn)指派給與運(yùn)輸成本最小的一個(gè)候選設(shè)施點(diǎn);3)選擇并取走一個(gè)設(shè)施點(diǎn),需要滿足以下條件:假如將它取走并將他的需求點(diǎn)重新分配后,總平均費(fèi)用增加量最?。?)從備選點(diǎn)中刪去取走點(diǎn),令=-1,然后轉(zhuǎn)第3)步;5)重復(fù)步驟3)、4),直到=(其中為所需選擇的配送中心數(shù)量)。

    2 實(shí)例研究

    2.1 研究區(qū)選擇

    案例中雪災(zāi)高風(fēng)險(xiǎn)區(qū)評估成果選擇本課題組前期青藏高原雪災(zāi)風(fēng)險(xiǎn)評估研究成果[38],研究基于歷史數(shù)據(jù),從雪災(zāi)致災(zāi)危險(xiǎn)性、承險(xiǎn)體物理暴露、敏感性和應(yīng)災(zāi)能力5個(gè)方面18項(xiàng)具體指標(biāo)對青藏高原雪災(zāi)風(fēng)險(xiǎn)進(jìn)行分析,結(jié)果表明青南高原部分地區(qū)的雪災(zāi)風(fēng)險(xiǎn)水平較高。該研究對前期雪災(zāi)風(fēng)險(xiǎn)圖的空值進(jìn)行修正,嘗試選取青海省雪災(zāi)高風(fēng)險(xiǎn)區(qū)進(jìn)行案例分析,以驗(yàn)證分析方法的可行性。

    2.2 研究區(qū)概況

    青海省雪災(zāi)高風(fēng)險(xiǎn)區(qū)大部分位于青海南部高原,是青海重要的牧業(yè)區(qū),其海拔3 500 m以上保持著完整的高原面,為青藏高原的腹地[39](圖2)。介于31°30′N~38°10′N和94°8′E~102°15′E,行政區(qū)域涉及玉樹藏族自治州稱多縣、囊謙縣等12縣,總面積約32萬km2,約占青海省總面積的45%。植被類型主要以高寒草甸、草原植被為主,分布著約14.62×104km2的草地,約占青海省草地面積的一半[40],畜牧業(yè)產(chǎn)值占總產(chǎn)值的60%以上[41],在全省農(nóng)業(yè)生產(chǎn)體系中具有極其重要的戰(zhàn)略地位。因其特殊的地理環(huán)境和氣候條件,冬春兩季雪災(zāi)頻發(fā),對當(dāng)?shù)氐男竽翗I(yè)生產(chǎn)產(chǎn)生了極其嚴(yán)重的影響。1961—2017年,青海省50個(gè)氣象臺(站)雪災(zāi)共發(fā)生756站次,其中輕災(zāi)486站次、中災(zāi)169站次、重災(zāi)77站次、特重度雪災(zāi)24站次,絕大多數(shù)集中分布于案例區(qū)內(nèi)。從各等級雪災(zāi)發(fā)生次數(shù)變化來看,近57年來,輕災(zāi)和重災(zāi)呈現(xiàn)微弱增多現(xiàn)象,中災(zāi)、特重度雪災(zāi)呈現(xiàn)減少現(xiàn)象[42]。

    圖2 研究區(qū)概況圖

    2019年研究區(qū)建設(shè)飼草料儲備庫總面積為2.48×104m2,可儲備飼草料7.90×104t。其中縣級飼草料儲備庫9座1.44×104m2,可儲備飼草料4.70×104t,分別建設(shè)在瑪多縣、天峻縣、稱多縣、囊謙縣、治多縣、瑪沁縣、河南縣、都蘭縣和剛察縣縣人民政府駐地;鄉(xiāng)鎮(zhèn)級飼草料儲備庫16座0.6×104m2,可儲備飼草料3.20×104t,分別建設(shè)在黃河鄉(xiāng)、扎陵湖鄉(xiāng)、扎朵鎮(zhèn)、著曉鄉(xiāng)、窩賽鄉(xiāng)、多松鄉(xiāng)等。但是其中大多數(shù)飼草料儲備庫在選址時(shí)按縣鄉(xiāng)行政區(qū)域進(jìn)行密集布點(diǎn),可能存在資源浪費(fèi)、分布不均和服務(wù)空白區(qū)等問題。該研究針對以上存在的問題對研究區(qū)飼草料儲備庫的選址做進(jìn)一步優(yōu)化。

    2.3 飼草料儲備庫選址分析過程

    2.3.1 牲畜飼草料需求分析過程

    青海省雪災(zāi)高風(fēng)險(xiǎn)區(qū)內(nèi)包含12個(gè)縣、67個(gè)鄉(xiāng)鎮(zhèn)和331個(gè)行政村。牲畜飼草料需求主要基于村級單位,通過飼草料需求量公式可計(jì)算出每個(gè)村落的牲畜飼草料需求量(表1)。該研究高風(fēng)險(xiǎn)區(qū)只包含大雪災(zāi)等級,故雪災(zāi)持續(xù)時(shí)間按10 d計(jì)算。因青藏高原牧區(qū)鄉(xiāng)鎮(zhèn)級數(shù)據(jù)不完備和村級數(shù)據(jù)較難獲取,故鄉(xiāng)鎮(zhèn)級及村級牲畜數(shù)量以縣級牲畜數(shù)量為基準(zhǔn)乘以人口占比和牧民人口占比得到。又因?yàn)轱暡萘闲枨罅繒驗(yàn)樯髷?shù)量的變化而變化,采用多年平均牲畜數(shù)量。在測算的基礎(chǔ)上,應(yīng)用ArcGIS軟件核密度工具分析求解出高風(fēng)險(xiǎn)區(qū)內(nèi)飼草料需求量的密度分布(圖3),結(jié)果顯示高密度區(qū)主要分布在剛察、都蘭、囊謙、雜多、治多、曲麻萊、稱多、達(dá)日、瑪沁、河南縣。

    表1 村級10 d飼草料需求量(部分)

    圖3 牲畜飼草料需求量密度

    2.3.2 最佳路徑分析過程

    首先在ArcGIS軟件中加載青海省雪災(zāi)高風(fēng)險(xiǎn)區(qū)內(nèi)的12個(gè)縣級、67個(gè)鄉(xiāng)鎮(zhèn)級、331個(gè)村級政府駐地以及1:100萬國家基礎(chǔ)道路矢量數(shù)據(jù)(數(shù)據(jù)來源:國家基礎(chǔ)地理信息中心)并對基礎(chǔ)道路進(jìn)行等級賦權(quán),其中道路等級權(quán)重按照各級道路行車速度(2004年《公路工程技術(shù)標(biāo)準(zhǔn)》(JTGB01-2003))確定(表2)。然后采用矢量數(shù)據(jù)結(jié)構(gòu)下交通網(wǎng)絡(luò)最短路徑算法,利用ArcGIS軟件中的Network Analyst分析模塊,分別計(jì)算出基于路網(wǎng)的縣級到鄉(xiāng)鎮(zhèn)級和鄉(xiāng)鎮(zhèn)級到村級的各需求點(diǎn)至各設(shè)施點(diǎn)(候選點(diǎn))的最佳路徑(表3)。

    表2 道路等級權(quán)重

    表3 縣級政府駐地到鄉(xiāng)級政府駐地的最佳路徑(部分)

    2.3.3 候選飼草料儲備庫設(shè)施點(diǎn)數(shù)量、容量及合理性

    因?yàn)檠芯繀^(qū)特殊的地形地貌、道路網(wǎng)的稀疏程度和防災(zāi)應(yīng)災(zāi)能力的差異,故飼草量儲備庫設(shè)施點(diǎn)建在縣、鄉(xiāng)人民政府駐地[43]。依據(jù)飼草料儲備庫選址原則,綜合考慮青藏高原地形、交通和草地生產(chǎn)力等限制因素,經(jīng)多次試驗(yàn)驗(yàn)證,飼草料儲備庫數(shù)量為3和儲備點(diǎn)數(shù)量為13時(shí),儲備庫/點(diǎn)在空間上分布均勻、達(dá)到了輻射面積最大和調(diào)運(yùn)時(shí)間最短的最佳選址要求。

    2.3.4 結(jié)果分析

    通過選址分析可得飼草料儲備庫為治多、瑪沁和都蘭飼草料儲備庫。飼草料儲備點(diǎn)為麻多、扎朵、秋智、多彩、結(jié)多、吉尼賽、白扎、拉加、香日德、龍藏、寧木特、建設(shè)和吉爾孟飼草料儲備點(diǎn)(圖4)。從案例區(qū)植被指數(shù)來看,儲備庫/點(diǎn)均分布于植被指數(shù)較高區(qū)域(圖5a),從而為飼草料的就地儲備提供了基礎(chǔ)條件。從交通便利程度來看,儲備庫/點(diǎn)均位于道路運(yùn)輸能力中高區(qū)(圖5b),從而為飼草料的及時(shí)調(diào)度提供了先決條件。

    圖4 飼草料儲備庫和儲備點(diǎn)的選址點(diǎn)、規(guī)模及服務(wù)范圍

    a. NDVI等級圖a. NDVI classification mapb. 道路運(yùn)輸能力等級圖b. Road transportation capability grade map

    治多飼草料儲備庫服務(wù)于7個(gè)飼草料儲備點(diǎn),分別為麻多、秋智、多彩、扎朵、結(jié)多、吉尼賽和白扎飼草料儲備點(diǎn),平均調(diào)運(yùn)時(shí)間和最長時(shí)間分別為3.3、5.1 h。服務(wù)人口約為17.5萬人。服務(wù)區(qū)草地面積約占青海省草地面積的33.59%。服務(wù)大牲畜約141.5萬頭(只、匹)。治多飼草料儲備庫的飼草料需求量為21.0萬 t,占儲備庫飼草料總需求量的43.94%。

    瑪沁飼草料儲備庫服務(wù)于4個(gè)飼草料儲備點(diǎn),分別為龍藏、拉加、建設(shè)和寧木特飼草料儲備點(diǎn),平均調(diào)運(yùn)時(shí)間和最長時(shí)間分別為2、2.7 h。服務(wù)人口約為11萬人。服務(wù)區(qū)草地面積占青海省草地面積的5.51%。服務(wù)大牲畜約55.7萬頭(只、匹)。瑪沁飼草料儲備庫的飼草料需求量為13.2萬t,占儲備庫飼草料總需求量的28.33%。

    都蘭飼草料儲備庫服務(wù)于2個(gè)飼草料儲備點(diǎn),分別是香日德和吉爾孟飼草料儲備點(diǎn),平均調(diào)運(yùn)時(shí)間和最長時(shí)間分別為2.3和3.5 h。服務(wù)人口約為8.5萬人。服務(wù)區(qū)草地面積占青海省草地面積的16.69%。服務(wù)大牲畜約37.3萬頭(只、匹)。都蘭飼草料儲備庫的飼草料需求量為13.5萬t,占儲備庫飼草料總需求量的27.73%。

    飼草料儲備點(diǎn)服務(wù)區(qū)域打破行政界線,實(shí)現(xiàn)統(tǒng)籌救災(zāi)的原則,以不同服務(wù)村點(diǎn)間的河流或山脈為界劃分,實(shí)現(xiàn)服務(wù)區(qū)域無空白,使得飼草量儲備點(diǎn)服務(wù)于服務(wù)區(qū)域上的所有村落。

    3 討 論

    3.1 P-中值模型在青藏高原選址中有著可行性

    從青藏高原牧區(qū)的生產(chǎn)方式來看,牧業(yè)生產(chǎn)使得牧區(qū)居民點(diǎn)相比于農(nóng)區(qū)居民點(diǎn)較為離散[44]。零散的聚落布局導(dǎo)致雪災(zāi)發(fā)生時(shí)救援的時(shí)效性顯得尤為重要,故而-中值模型計(jì)量的平均權(quán)重對于青藏高原牧區(qū)聚落的離散分布具有很大的意義[45]。從飼草料需求量來看,青藏高原牧區(qū)各需求點(diǎn)飼草料需求量各異,雪災(zāi)發(fā)生時(shí)飼草料發(fā)放優(yōu)先程度也應(yīng)不同。以需求點(diǎn)的飼草料需求量對距離進(jìn)行加權(quán),從而為需求量大的需求點(diǎn)優(yōu)先指派臨近儲備庫點(diǎn),進(jìn)而整體優(yōu)化因承災(zāi)體脆弱性不同而導(dǎo)致的損失。從飼草料分配方式來看,該模型的分配規(guī)則為個(gè)設(shè)施找到合適的位置并指派每個(gè)需求點(diǎn)到一個(gè)特定的設(shè)施。該規(guī)則可以使得青藏高原高風(fēng)險(xiǎn)區(qū)內(nèi)的每個(gè)需求點(diǎn)都能有專一的飼草料儲備庫點(diǎn)進(jìn)行服務(wù)。綜上所述,-中值模型在青藏高原飼草料儲備庫的選址中有著一定的可行性。但是該模型是單目標(biāo)模型[46],該研究主要考慮了飼草料需求量的影響,卻未考慮其他選址目標(biāo)的影響,如管理成本、服務(wù)滿意度等,故在選址綜合性考慮上有所欠缺。

    3.2 基于選址影響因素的案例區(qū)選址結(jié)果具有合理性

    案例分析結(jié)果顯示,治多、瑪沁和都蘭飼草料儲備庫,分別位于玉樹州、果洛州和海西州,在研究區(qū)內(nèi)分布均勻,覆蓋范圍輻射整個(gè)研究區(qū)。從交通便利程度看,3座飼草料儲備庫均位于交通節(jié)點(diǎn),有著交通便利的優(yōu)勢。但是研究區(qū)內(nèi)整體道路網(wǎng)密度較小,在一定程度上制約了飼草料的有效調(diào)度。故雪災(zāi)發(fā)生時(shí),保障道路通達(dá)性顯得尤為重要[47]。依據(jù)研究區(qū)高寒牧區(qū)的特點(diǎn)[48],飼草料儲備庫點(diǎn)建設(shè)在縣、鄉(xiāng)/鎮(zhèn)人民政府駐地是合理的。既可以解決地形、地質(zhì)等天然條件的影響,又可以方便政府的應(yīng)急指揮。打破行政界限進(jìn)行飼草料儲備庫/點(diǎn)的管理與服務(wù),在降低牧民遭遇雪災(zāi)的損失和提高政府應(yīng)急救災(zāi)防災(zāi)的能力方面起到顯著作用。飼草料儲備庫規(guī)模的準(zhǔn)確計(jì)算,是政府按需備料,科學(xué)防控的參考依據(jù)。治多、瑪沁和都蘭飼草料儲備庫周邊主要是純牧業(yè)鄉(xiāng)鎮(zhèn),飼草料需求量較大。對于飼草料需求量大的村落臨近分配可以縮短運(yùn)輸距離和提高應(yīng)急能力。從可利用草場面積來看,治多縣2.6×104km2、瑪沁縣1.1× 104km2、都蘭1.5×104km2在研究區(qū)草場面積中占據(jù)了很大的比例。豐富的草地資源可以就地為飼草料儲備庫/點(diǎn)提供儲備物資。結(jié)合飼草料儲備庫的選址影響因素來看,案例分析區(qū)的選址結(jié)果是合理的。該研究因牧區(qū)鄉(xiāng)鎮(zhèn)級數(shù)據(jù)不完備和村級牲畜數(shù)量較難獲取,故通過計(jì)算獲得的村級數(shù)據(jù)有所不足,未能考慮到養(yǎng)殖大戶、合作社大規(guī)模養(yǎng)殖等情況,所以在后續(xù)研究中會對此做進(jìn)一步研究。

    4 結(jié) 論

    該研究以青海省雪災(zāi)高風(fēng)險(xiǎn)區(qū)為案例,分析可得在案例區(qū)建設(shè)3座縣級飼草料儲備庫,13座鄉(xiāng)鎮(zhèn)級飼草料儲備點(diǎn),在空間上分布均勻,達(dá)到服務(wù)范圍輻射研究區(qū)的最佳選址要求,為研究區(qū)儲備庫/點(diǎn)的選址優(yōu)化提供了參考依據(jù)。從分布結(jié)構(gòu)上看,形成了以縣級儲備庫為中心,鄉(xiāng)鎮(zhèn)級儲備點(diǎn)為外圍的二元結(jié)構(gòu)。本案例驗(yàn)證了-中值模型在青藏高原進(jìn)行飼草料儲備庫選址的優(yōu)越性和該研究分析方法的可行性。如果本案例可以成功的應(yīng)用到現(xiàn)實(shí)儲備庫的選址中,將會更好的解決布局不合理和資源浪費(fèi)等問題,在防災(zāi)減災(zāi)中發(fā)揮至關(guān)重要的作用。

    [1] 史培軍,王愛慧,孫福寶,等. 全球變化人口與經(jīng)濟(jì)系統(tǒng)風(fēng)險(xiǎn)形成機(jī)制及評估研究[J]. 地球科學(xué)進(jìn)展,2016,31(8):775-781.

    Shi Peijun, Wang Aihui, Sun Fubao, et al. A study of global change population and economic system risk forming mechanism and assessment[J]. Advances in Earth Science, 2016, 31(8):775-781. (in Chinese with English abstract)

    [2] National Academies of Sciences, Engineering, and Medicine. Attribution of Extreme Weather Events in the Context of Climate Change[M]. Washington DC: The National Academies Press, 2016.

    [3] 郭躍. 災(zāi)害范式及其歷史演進(jìn)[J]. 地理科學(xué),2016,36(6):935-942.

    Guo Yue. Hazard paradigms and their historical evolution[J]. Scientia Geographica Sinica, 2016, 36(6): 935-942. (in Chinese with English abstract)

    [4] UN/ISDR (United Nations International Strategy for Disaster Reduction). Sendai framework for disaster risk reduction 2015-2030[EB/OL]. [2017-06-10].http:∥www.unisdr.org/ files/43291_sendaiframeworkfordrren.pdf.

    [5] 吳瑋,秦其明,范一大,等. 中國雪災(zāi)評估研究綜述[J]. 災(zāi)害學(xué),2013,28(4):152-158.

    Wu Wei, Qin Qiming, Fan Yida, et al. Review on snow disaster assessment in China[J]. Journal of Catastrophology, 2013, 28(4): 152-158. (in Chinese with English abstract)

    [6] 榮莉莉,劉玙婷. 基于承災(zāi)體的區(qū)域?yàn)?zāi)害鏈風(fēng)險(xiǎn)評估模型[J]. 系統(tǒng)工程學(xué)報(bào),2019,34(1):130-144.

    Rong Lili, Liu Yuting. Modeling the evaluation for regional disaster chain risk based on elements at risk[J]. Journal of Systems Engineering, 2019, 34(1): 130-144.

    [7] Dabbeek J, Silva V. Modeling the residential building stock in the Middle East for multi-hazard risk assessment[J]. Natural Hazards: Journal of the International Society for the Prevention and Mitigation of Natural Hazards, 2020, 100.

    [8] 周洪建. 當(dāng)前全球減輕災(zāi)害風(fēng)險(xiǎn)平臺的前沿話題與展望:基于2017年全球減災(zāi)平臺大會的綜述與思考[J]. 地球科學(xué)進(jìn)展,2017,32(7):688-695.

    Zhou Hongjian. Hot-topics and prospects of global platform for disaster risk reduction: Based on 2017 global platform for disaster risk reduction in Cancun, Mexico[J]. Advances in Earth Science, 2017, 32(7): 688-695. (in Chinese with English abstract)

    [9] 高慶華,馬宗晉,張業(yè)成. 自然災(zāi)害評估[M]. 北京:氣象出版社,2007.

    [10] 全國氣象防災(zāi)減災(zāi)標(biāo)準(zhǔn)化技術(shù)委員會. 牧區(qū)雪災(zāi)等級:GB/T20482-2006 [S]. 北京:中國標(biāo)準(zhǔn)出版社,2006.

    [11] 莊曉翠,周鴻奎,王磊,等. 新疆北部牧區(qū)雪災(zāi)評估指標(biāo)及其成因分析[J]. 干旱區(qū)研究,2015,32(5):1000-1006.

    Zhuang Xiaocui, Zhou Hongkui, Wanglei, et al. Evaluation and cause study on the snow disasters in pastoral areas of northern Xinjiang[J]. Arid Zone Research, 2015, 32(5): 1000-1006. (in Chinese with English abstract)

    [12] 郝璐,王靜愛,滿蘇爾,等. 中國雪災(zāi)時(shí)空變化及畜牧業(yè)脆弱性分析[J]. 自然災(zāi)害學(xué)報(bào),2002,11(4):42-48.

    Hao Lu, Wang Jingai, Man Suer, et al. Spatio-temporal change of snow disaster and analysis of vulnerability of animal husbandry in Chin[J]. Journal of Natural Disasters, 2002, 11(4): 42-48. (in Chinese with English abstract)

    [13] Dong Shikui, Sherman Ruth. Enhancing the resilience of coupled human and natural systems of alpine rangelands on the Qinghai-Tibetan Plateau[J]. The Rangeland Journal, 2015, 37(1). doi:10.1071/RJ14117.

    [14] Li Yijia, Ye Tao, Liu Weihang, et al. Linking livestock snow disaster mortality and environmental stressors in the Qinghai-Tibetan Plateau: Quantification based on generalized additive models[J]. Science of the Total Environment, 2018, 625: 87-95.

    [15] 王莘. 中國氣象災(zāi)害大典(青海卷)[M]. 北京:氣象出版社,2007.

    [16] Wang Shijin, Zhou Lanyue, Wei Yanqiang. Integrated risk assessment of snow disaster over the Qinghai-Tibet Plateau[J]. Geomatics Natural Hazards and Risk, 2019, 10(1): 740-757.

    [17] IPCC. Climate change 2014: The fourth assessment report of the inter-governmental panel on climate change[R]. UK: Cambridge University Press, 2014.

    [18] 張繼權(quán),岡田憲夫,多多納裕一. 綜合自然災(zāi)害風(fēng)險(xiǎn)管理:全面整合的模式與中國的戰(zhàn)略選擇[J]. 自然災(zāi)害學(xué)報(bào),2006(1):29-37.

    Zhang Jiquan, Okada Norio, Tatano Hirokazu. Integrated natural disaster risk management: Comprehensive and integrated model and Chinese strategy choice[J]. Journal of Natural Disasters, 2006(1): 29-37. (in Chinese with English abstract)

    [19] Wang Shijin, Chen Shengyun, Wei Yanqiang. Risk prevention and control strategies for the severely affected areas of snow disaster in the Three Rivers Source Region(TRSR)[J]. China. Sciences in Cold and Arid Regions, 2019, 11(3): 248-252.

    [20] 張強(qiáng),劉峰貴. 創(chuàng)新智慧牧區(qū)發(fā)展提升雪災(zāi)防抗救能力[N]. 中國應(yīng)急管理報(bào),2019-05-28(007).

    [21] 蔡冬雪,朱建明,王國慶. 基于情景分析的應(yīng)急裝備多層級協(xié)同布局問題研究[J]. 中國管理科學(xué),2017,25(10):72-79.

    Cai Dongxue, Zhu Jianming, Wang Guoqing. Hierarchical Collaborative Location and allocation of emergency equipment based on Scenario analysis[J]. Chinese Journal of Management Science, 2017, 25(10): 72-79. (in Chinese with English abstract)

    [22] 孫清臣,曲林遲. 考慮多類應(yīng)急物資需求不確定的選址模型[J]. 廣西大學(xué)學(xué)報(bào):自然科學(xué)版,2019,44(2):448-454.

    Sun Qingchen, Qu Linchi. A location model and with consideration of uncertain demand for multiple relief supplies[J]. Journal of Guangxi University: Natural Science Edition, 2019, 44(2): 448-454. (in Chinese with English abstract)

    [23] 俞武揚(yáng). 服務(wù)能力受損情景下的應(yīng)急設(shè)施選址模型[J]. 控制與決策,2016,31(11):1979-1984.

    Yu Wuyang. Emergency facility location model under service ability damage scenarios[J]. Control and Decision, 2016, 31(11): 1979-1984.

    [24] 付德強(qiáng),陳子豪,蹇潔,等. 應(yīng)急聯(lián)動區(qū)域下選址分配協(xié)同優(yōu)化模型研究[J]. 數(shù)學(xué)的實(shí)踐與認(rèn)識,2019,49(6):30-41.

    Fu Deqiang, Chen Zihao, Jian Jie, et al. Research of resource location allocation model in regional emergency joint action[J]. Mathematics in Practice and Theory, 2019, 49(6): 30-41. (in Chinese with English abstract)

    [25] 郗蒙浩,張靜,趙秋紅,等. 基于P-center問題的國家級應(yīng)急物資儲備設(shè)施選址優(yōu)化布局研究[J]. 自然災(zāi)害學(xué)報(bào),2019,28(3):123-129.

    Xi Menghao, Zhang Jing, Zhao Qiuhong, et al. The research of location problem of national emergency reserve facilities for emergency relief based on P-center problem[J]. Journal of Natural Disasters, 2019, 28(3): 123-129. (in Chinese with English abstract)

    [26] Cheng Z H, Wang Q, Du M, et al. Research on location model of emergency material storage warehouse based on generalized maximum coverage[J]. IOP Conference Series Materials Science and Engineering, 2020, 792: 012013.

    [27] 叢雯婧,俞武揚(yáng). 考慮臺風(fēng)情景的區(qū)域應(yīng)急物資儲備庫選址模型與算法[J]. 工業(yè)工程與管理. doi:1738.T.20191113.1758.008.html

    Cong Wenjing, Yu Wuyang. Model and algorithm for location of regional emergency material reserve site in typhoon scenarios[J]. Industrial Engineering and Management. doi:1738.T.20191113.1758.008.html

    [28] 高淑春. 需求不確定條件下的應(yīng)急物資儲備庫選址模型構(gòu)建研究[D]. 安徽:合肥工業(yè)大學(xué),2017.

    Gao Shuchun. Research on the Model of the Location of Emergency Supplies Repository under Demand Uncertainty[D]. Anhui: Hefei University of Technology, 2017. (in Chinese with English abstract)

    [29] Hakimi S L. Optimum distribution of switching centers in a communication network and some related graph theoretic problems[J]. Operations Research, 1965, 13(3): 462-475.

    [30] Church R, ReVelle C. The maximal covering location problem[J]. Papers in Regional Science, 1974, 32(1): 101-118.

    [31] Hakimi S L. Optimum locations of switching centers and the absolute centers and medians of a graph[J]. Operations Research, 1964, 12(3): 450-59.

    [32] 王磊,潘進(jìn),焦李成. 免疫算法[J]. 電子學(xué)報(bào),2000(7):74-78.

    Wang Lei, Pan Jin, Jiao Licheng. The immune algorithm[J]. Acta Electronica Sinica, 2000(7):74-78. (in Chinese with English abstract)

    [33] 李倩,張惠珍,Cesar Beltran-Royo. 求解無容量設(shè)施選址問題的混合蟻群算法[J]. 上海理工大學(xué)學(xué)報(bào),2016,38(4):367-372.

    Li Qian, Zhang Huizhen, Cesar Beltran-Royo. Hybrid ant colony algorithm for the uncapacitated facility location problem[J]. Journal of University of Shanghai for Science and Technology, 2016, 38(4): 367-372. (in Chinese with English abstract)

    [34] 關(guān)懷慶,張畢西,歐江艷. 貪婪取走啟發(fā)式算法在離散網(wǎng)絡(luò)選址中的研究[J]. 系統(tǒng)科學(xué)學(xué)報(bào),2010,18(3):49-52,56.

    Guan Huaiqing, Zhang Bixi, Ou Jiangyan. The study of greedy dropping heuristic algorithm in discrete network location[J]. Chinese Journal of Systems Science, 2010, 18(3): 49-52, 56. (in Chinese with English abstract)

    [35] 蔡臨寧. 物流系統(tǒng)規(guī)劃:建模及實(shí)例分析[M]. 北京:機(jī)械工業(yè)出版社,2003.

    [36] 宋正娜,顏庭干,劉婷,等. 新重力P中值模型及其在城市綜合醫(yī)院區(qū)位決策中的實(shí)證檢驗(yàn):以無錫市為例[J]. 地理科學(xué)進(jìn)展,2016,35(4):420-430.

    Song Zhengna, Yan Tinggan, Liu Ting, et al. A new gravity P-median model and empirical test in urban comprehensive hospital location decision making: Take Wuxi as an example[J]. Progress in Geography, 2016, 35(4): 420-430. (in Chinese with English abstract)

    [37] 周陸生,汪青春,李海紅,等. 青藏高原東部牧區(qū)大-暴雪過程雪災(zāi)災(zāi)情實(shí)時(shí)預(yù)評估方法的研究[J]. 自然災(zāi)害學(xué)報(bào),2001(2):58-65.

    Zhou Lusheng, Wang Qingchun, Li Haihong, et al. Study on real-time predictive assessment of snowstorm disaster in eastern pastoral area of Qinghai-Tibet Plateau[J]. Journal of Natural Disasters, 2001(2): 58-65. (in Chinese with English abstract)

    [38] Liu Fenggui, Mao Xufeng, Zhang Yili, et al. Risk analysis of snow disaster in the pastoral areas of the Qinghai-Tibet Plateau[J]. Journal of Geographical Sciences, 2014, 24(3): 411-426.

    [39] 張忠孝. 青海地理[M]. 西寧:青海人民出版社,2004.

    [40] 竇永紅. 基于AHP的青南高原高寒草地脆弱性評估[J]. 青海環(huán)境,2018,28(3):145-150.

    Dou Yonghong. AHP-based evaluation on alpine grassland weakness on Qingnan Plateau[J]. Journal of Qinghai Environment, 2018, 28(3): 145-150. (in Chinese with English abstract)

    [41] 吳小輝,侯光良,王記明,等. 1950-2000年青海省雪災(zāi)災(zāi)情的時(shí)空分析[J]. 青海環(huán)境,2015,25(2):59-64.

    Wu Xiaohui, Hou Guangliang, Wang Jiming, et al. Spatio-temporal analysis of snow disaster in Qinghai province from 1950 to 2000[J]. Journal of Qinghai Environment, 2015, 25(2): 59-64. (in Chinese with English abstract)

    [42] 邵全琴,劉國波,李曉東,等. 三江源區(qū)2019年春季雪災(zāi)及草地畜牧業(yè)雪災(zāi)防御能力評估[J]. 草地學(xué)報(bào),2019,27(5):1317-1327.

    Shao Quanqin, Liu Guobo, Li Xiaodong, et al. Assessment of the defense Capability of spring snow disaster and grassland animal husbandary in 2019 in the Three-river Source region[J]. Acta Agrestia Sinica, 2019, 27(5):1317-1327. (in Chinese with English abstract)

    [43] 高興川,曹小曙,李濤,等. 1976—2016年青藏高原地區(qū)通達(dá)性空間格局演變[J]. 地理學(xué)報(bào),2019,74(6):1190-1204.

    Gao Xingchuan, Cao Xiaoshu, Li Tao, et al. Evolution of accessibility spatial pattern of the Qinghai-Tibet Plateau in 1976-2016[J]. Acta Geographica Sinica, 2019, 74(6): 1190-1204. (in Chinese with English abstract)

    [44] 張海朋,何仁偉,劉運(yùn)偉,等. 青藏高原高寒牧區(qū)聚落用地適宜性評價(jià)及其重構(gòu):以藏北那曲縣為例[J]. 自然資源學(xué)報(bào),2020,35(3):698-712.

    Zhang Haipeng, He Renwei, Liu Yunwei, et al. Land suitability evaluation and reconstruction of settlements in the pastoral area of Tibetan Plateau: Aase study of nagqu county in Northern tibet [J]. Journal of Natural Resources, 2020, 35(3): 698-712. (in Chinese with English abstract)

    [45] 王雪,常志州,王效華,等. 泗洪縣車門鄉(xiāng)稻麥秸稈收儲設(shè)施選址[J]. 農(nóng)業(yè)工程學(xué)報(bào),2015,31(22):250-255.

    Wang Xue, Chang Zhizhou, Wang Xiaohua, et al. Crop straw collecting and storage facility location for Chemen township of Sihong county[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(22): 250-255. (in Chinese with English abstract)

    [46] 盧媛媛,郭志堅(jiān),王寶玲,等. 趨近理想灰關(guān)聯(lián)多約束的交通救援中心中值選址模型[J]. 公路交通科技,2018,35(7):114-119.

    Lu Yuanyuan, Guo Zhijian, Wang Baoling, et al. A traffic rescue centre median location model approaching ideal grey relational constraints[J]. Journal of Highway and Transportation Research and Development, 2018, 35(7): 114-119. (in Chinese with English abstract)

    [47] 鄒銘,袁藝,廖永豐,等.綜合風(fēng)險(xiǎn)防范:中國綜合自然災(zāi)害救助保障體系[M]. 北京:科學(xué)出版社,2011.

    [48] 劉柱. 民族地區(qū)常態(tài)災(zāi)害防災(zāi)減災(zāi)體系構(gòu)建[D]. 長沙:中南大學(xué),2010.

    Liu Zhu. Construction of Disaster Prevention and Reduction System for Normalcy Disasters in Ethnic Minority Areas[D]. Chang Sha: Central South University, 2010. (in Chinese with English abstract)

    Site selection of forage storehouse in high risk areas of snow disaster in Qinghai-Tibet Plateau

    Zhao Fuchang1, Zhou Qiang1,2※, Liu Fenggui1,2, Chen Qiong1,2, Chen Yongping1, Zhi Zemin1

    (1.,,810008,; 2.810008,)

    Snow disaster become a gradually increasing trend, due to the global climate change, particularly on the frequency, intensity, and duration in high risk area of extreme weather. In addition, the prevention and reduction of natural disaster have posed a major challenge on the international community. Therefore, it is urgent to accurately predict the disaster risk, in order to reduce the losses from the snow disaster. Being one of the most important bases for stockbreeding, the Qinghai-Tibetan Plateau is a typical region with high incidence of snow disasters in western China. In the high risk area of snow disaster,forage reserves play a vital role in the emergency response to snow disasters in pastoral areas. The objectives of this study were, 1) to evaluate the high risk area of snow disaster, and to verify the accuracy using the actual occurrence areas; 2) to propose a new method for the site selection of forage reserves by modifying P-median model; 3) to explore the optimum sites for forage reserves over the Qinghai-Tibetan Plateau. The high-risk areas of snow disaster and the occurrence areas in actual years were combined, in order to better evaluate the risk of snow disaster in disaster risk management, and thereby to improve the disaster resilience of plateau pastoral areas. The number of livestock in county-level administrative units and vector data of national basic roads were collected from government units and the National Basic Geographic Information Center. The P-median model was selected to evaluate the sites selection, compared with other models, such as P-center model, maximum coverage model, P-median model. In addition, a modified P-median model was proposed to evaluate the optimum sites of forage reserves in this study area, according to the calculation demand for livestock forage, and the optimal path from the facility to the demand of models. The high-risk area of snow disaster for the Qinghai Province can serve as a case to verify the models. The results showed that 1) Zhiduo, Maqin and Dulan counties were the optimum forage reserves in this study area, which the number of their forage sites was different, such as 7 sites in Zhiduo county (including Maduo, Qiuzhi, Duocai, Zhaduo, Jieduo, Jinisai, and Baizha towns), 4 sites in Maqin (including Longzang, Lajia, Jianshe, and Ningmute towns), and 2 sites in Dulan county (including Xiangride and Jiermeng towns); 2) There was an obvious difference among sizes, where the capacities of Zhiduo, Maqin, and Dulan forage reserves were 210 000, 132 000 and 135 000 ton of forage grass, respectively; 3) In space, a dual structural pattern was formed in “networked layout, hierarchical reserve”, where Zhidu, Maqin and Dulan reserves served as the center, Maduo, Zado, Qiuzhi and other forage reserves as the periphery, covering the entire study area. The finding demonstrated that the P-median model can be expected to perform excellent sites selection of forage reserves in the Qinghai-Tibetan Plateau. The finding can provide an available theoretical method for the selection of forage stock reserve in plateau pastoral area.

    risk assessment; silage forage; disasters; forage reserves; location selection; the-median model; the high risk area of snow disaster

    趙富昌,周強(qiáng),劉峰貴,等. 青藏高原雪災(zāi)高風(fēng)險(xiǎn)區(qū)飼草料儲備庫選址分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(13):272-279.doi:10.11975/j.issn.1002-6819.2020.13.032 http://www.tcsae.org

    Zhao Fuchang, Zhou Qiang, Liu Fenggui, et al. Site selection of forage storehouse in high risk areas of snow disaster in Qinghai-Tibet Plateau[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(13): 272-279. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2020.13.032 http://www.tcsae.org

    2020-03-19

    2020-05-12

    第二次青藏高原綜合科學(xué)考察研究項(xiàng)目(2019QZKK0906);國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2019YFA0606902)

    趙富昌,主要從事土地資源與環(huán)境變化研究。Email:15297113505@163.com

    周強(qiáng),教授,主要從事土地資源與環(huán)境變化研究。Email:zhouqiang729@163.com

    10.11975/j.issn.1002-6819.2020.13.032

    S427

    A

    1002-6819(2020)-13-0272-08

    猜你喜歡
    雪災(zāi)儲備庫飼草料
    飼草料供給對小型奶牛養(yǎng)殖戶經(jīng)濟(jì)效益的影響
    考慮飼草料質(zhì)量的家庭牧場水土草畜平衡計(jì)算方法與應(yīng)用
    青海海西地區(qū)近41年雪災(zāi)時(shí)空分布特征
    淺談發(fā)展飼草料存在的問題及建議
    浙江省糧食局直屬糧油儲備庫:人才殷倉廩 創(chuàng)新促發(fā)展
    陜西省靖邊糧食儲備庫:構(gòu)建糧食產(chǎn)后服務(wù)體系 提升服務(wù)“三農(nóng)”水平
    歷史時(shí)期以來山西省雪災(zāi)特征與氣候變化
    俄批準(zhǔn)建立金磚國家外匯儲備庫
    近51 a青藏高原雪災(zāi)時(shí)空分布特征
    2000-2010年新疆雪災(zāi)時(shí)空自相關(guān)分析
    国产精品久久久久成人av| 午夜福利在线观看吧| 男女无遮挡免费网站观看| 国产精品免费视频内射| a级毛片黄视频| 欧美乱码精品一区二区三区| 99久久精品国产亚洲精品| 啪啪无遮挡十八禁网站| 天天影视国产精品| 超碰97精品在线观看| 老司机影院成人| 欧美少妇被猛烈插入视频| 一二三四社区在线视频社区8| 日韩欧美国产一区二区入口| 亚洲中文日韩欧美视频| 首页视频小说图片口味搜索| 伊人亚洲综合成人网| 一级a爱视频在线免费观看| 成人国产av品久久久| 黄色a级毛片大全视频| 亚洲伊人久久精品综合| 两性夫妻黄色片| 国产免费现黄频在线看| 中文字幕制服av| 91老司机精品| 少妇人妻久久综合中文| 精品福利永久在线观看| 99精品欧美一区二区三区四区| 午夜免费成人在线视频| 五月天丁香电影| 日日夜夜操网爽| 国产区一区二久久| 叶爱在线成人免费视频播放| 欧美变态另类bdsm刘玥| 黄色 视频免费看| 中文字幕av电影在线播放| 男男h啪啪无遮挡| 国产在线视频一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 国产无遮挡羞羞视频在线观看| 久久久久国产精品人妻一区二区| 亚洲人成电影免费在线| 亚洲av电影在线进入| 啦啦啦免费观看视频1| 动漫黄色视频在线观看| 国产成人av激情在线播放| 人妻久久中文字幕网| 热99re8久久精品国产| 国产亚洲精品久久久久5区| 午夜福利乱码中文字幕| 国产成人欧美| 婷婷色av中文字幕| 久久人人爽人人片av| 国产成人免费无遮挡视频| www.熟女人妻精品国产| 日韩中文字幕视频在线看片| 99热国产这里只有精品6| 亚洲精品粉嫩美女一区| 欧美变态另类bdsm刘玥| 在线观看一区二区三区激情| 极品少妇高潮喷水抽搐| 午夜福利视频在线观看免费| 男人舔女人的私密视频| 亚洲精品在线美女| 欧美日韩精品网址| 黑人操中国人逼视频| 黄色a级毛片大全视频| 亚洲欧洲精品一区二区精品久久久| 精品国产一区二区久久| 如日韩欧美国产精品一区二区三区| 亚洲黑人精品在线| 欧美日韩中文字幕国产精品一区二区三区 | 成人手机av| 黄色视频,在线免费观看| 一二三四社区在线视频社区8| 亚洲久久久国产精品| 日韩 欧美 亚洲 中文字幕| 不卡一级毛片| 丝袜脚勾引网站| 丝袜美足系列| 50天的宝宝边吃奶边哭怎么回事| 久久久久精品国产欧美久久久 | 欧美亚洲日本最大视频资源| 亚洲国产看品久久| 国产一区二区激情短视频 | 99久久国产精品久久久| 在线看a的网站| 少妇裸体淫交视频免费看高清 | 大型av网站在线播放| 99re6热这里在线精品视频| av欧美777| 巨乳人妻的诱惑在线观看| 91成人精品电影| 国产精品麻豆人妻色哟哟久久| 久久人人爽人人片av| 色视频在线一区二区三区| 国产人伦9x9x在线观看| 一本色道久久久久久精品综合| 欧美日韩一级在线毛片| 免费看十八禁软件| 国产精品一区二区在线观看99| 午夜福利免费观看在线| 国产三级黄色录像| 久久中文看片网| 女人高潮潮喷娇喘18禁视频| 在线亚洲精品国产二区图片欧美| 国产精品国产三级国产专区5o| 亚洲欧美精品综合一区二区三区| 国产精品一区二区在线不卡| 日韩制服骚丝袜av| 少妇裸体淫交视频免费看高清 | 午夜福利在线免费观看网站| tube8黄色片| 欧美 日韩 精品 国产| 久久久久久久大尺度免费视频| 女性被躁到高潮视频| 精品第一国产精品| 青春草亚洲视频在线观看| 日本撒尿小便嘘嘘汇集6| 国产色视频综合| 国产成人免费观看mmmm| 亚洲av男天堂| 久久久久久久精品精品| 国产亚洲欧美精品永久| 精品久久久久久电影网| 国产国语露脸激情在线看| 一进一出抽搐动态| 在线观看免费日韩欧美大片| 免费看十八禁软件| 1024视频免费在线观看| 久久久久国内视频| 亚洲av国产av综合av卡| 亚洲黑人精品在线| 国产成人精品久久二区二区91| 少妇裸体淫交视频免费看高清 | 午夜福利在线观看吧| 高清黄色对白视频在线免费看| 欧美精品人与动牲交sv欧美| 狠狠狠狠99中文字幕| 高清黄色对白视频在线免费看| 久久人妻福利社区极品人妻图片| 久久久国产精品麻豆| 大香蕉久久网| 国产成人精品在线电影| 亚洲性夜色夜夜综合| 亚洲性夜色夜夜综合| 午夜福利在线免费观看网站| 一个人免费在线观看的高清视频 | 国产福利在线免费观看视频| 美女脱内裤让男人舔精品视频| a级片在线免费高清观看视频| 人妻 亚洲 视频| 天天影视国产精品| 久久精品熟女亚洲av麻豆精品| 亚洲第一av免费看| 国产成人系列免费观看| 91成人精品电影| 午夜影院在线不卡| 国产在线观看jvid| av在线老鸭窝| 青春草视频在线免费观看| av欧美777| 美女高潮喷水抽搐中文字幕| av有码第一页| 国产精品自产拍在线观看55亚洲 | 亚洲精品av麻豆狂野| 我要看黄色一级片免费的| 国产精品久久久久久精品电影小说| 老司机在亚洲福利影院| 国产精品久久久久久精品古装| 黄色 视频免费看| 乱人伦中国视频| 日韩视频在线欧美| 午夜福利乱码中文字幕| 1024香蕉在线观看| 女人久久www免费人成看片| 色94色欧美一区二区| 亚洲欧美精品自产自拍| 国产一区二区三区av在线| 国产97色在线日韩免费| 老司机午夜十八禁免费视频| 亚洲avbb在线观看| 亚洲精品av麻豆狂野| 欧美日韩av久久| 老汉色∧v一级毛片| 免费少妇av软件| 成年动漫av网址| 宅男免费午夜| 欧美激情极品国产一区二区三区| 亚洲精品国产精品久久久不卡| 纵有疾风起免费观看全集完整版| 后天国语完整版免费观看| bbb黄色大片| 少妇的丰满在线观看| 精品一区在线观看国产| 久久久久久久精品精品| 精品卡一卡二卡四卡免费| 不卡av一区二区三区| 精品久久久精品久久久| 亚洲欧美精品综合一区二区三区| 黄色 视频免费看| 在线亚洲精品国产二区图片欧美| 日韩 欧美 亚洲 中文字幕| 精品乱码久久久久久99久播| 久久久久精品人妻al黑| 天天添夜夜摸| 高清黄色对白视频在线免费看| 国产精品成人在线| 久久久国产成人免费| 日本av免费视频播放| 国产麻豆69| 久久99热这里只频精品6学生| 人人妻人人添人人爽欧美一区卜| 美女国产高潮福利片在线看| 妹子高潮喷水视频| 欧美精品高潮呻吟av久久| 日本欧美视频一区| 午夜福利在线免费观看网站| 丝袜人妻中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕最新亚洲高清| 在线十欧美十亚洲十日本专区| 国产无遮挡羞羞视频在线观看| 久久精品aⅴ一区二区三区四区| www.av在线官网国产| av网站免费在线观看视频| 精品一区二区三区四区五区乱码| 国产亚洲欧美精品永久| 久久国产亚洲av麻豆专区| 侵犯人妻中文字幕一二三四区| 丝袜在线中文字幕| 欧美精品高潮呻吟av久久| 日本黄色日本黄色录像| 天堂中文最新版在线下载| 国产欧美日韩一区二区精品| 久久久精品国产亚洲av高清涩受| svipshipincom国产片| 天堂8中文在线网| 国产精品熟女久久久久浪| 国产精品久久久久久精品古装| 精品人妻1区二区| 午夜福利,免费看| 精品国产乱码久久久久久男人| 免费观看a级毛片全部| 亚洲欧洲日产国产| 精品卡一卡二卡四卡免费| 99国产精品一区二区三区| 黄色视频在线播放观看不卡| 亚洲精品一卡2卡三卡4卡5卡 | 色老头精品视频在线观看| 曰老女人黄片| 侵犯人妻中文字幕一二三四区| 天天影视国产精品| 在线精品无人区一区二区三| 91av网站免费观看| 免费高清在线观看日韩| 在线观看免费日韩欧美大片| 亚洲五月婷婷丁香| 桃红色精品国产亚洲av| 国产麻豆69| 人人澡人人妻人| 在线 av 中文字幕| 国产在线一区二区三区精| kizo精华| 啦啦啦 在线观看视频| www日本在线高清视频| 老汉色∧v一级毛片| 亚洲国产毛片av蜜桃av| 人妻一区二区av| 一个人免费在线观看的高清视频 | 国产熟女午夜一区二区三区| 99久久精品国产亚洲精品| 国产人伦9x9x在线观看| 一本大道久久a久久精品| 久久人人爽人人片av| 午夜免费成人在线视频| 深夜精品福利| 纵有疾风起免费观看全集完整版| 99精品久久久久人妻精品| 大片电影免费在线观看免费| 国产成人av激情在线播放| 亚洲一区中文字幕在线| 新久久久久国产一级毛片| 1024香蕉在线观看| 午夜福利在线观看吧| 国产国语露脸激情在线看| 色播在线永久视频| 亚洲成人免费电影在线观看| 青春草视频在线免费观看| 国产欧美亚洲国产| 欧美亚洲 丝袜 人妻 在线| av片东京热男人的天堂| 91精品伊人久久大香线蕉| 国产精品免费视频内射| 丝袜喷水一区| 人人妻人人澡人人爽人人夜夜| bbb黄色大片| 青春草亚洲视频在线观看| 亚洲精品一卡2卡三卡4卡5卡 | 国产野战对白在线观看| a级片在线免费高清观看视频| 午夜福利免费观看在线| 久久久国产成人免费| 日本av手机在线免费观看| 国产精品国产三级国产专区5o| 精品国产一区二区三区四区第35| 亚洲国产日韩一区二区| 国产一区二区在线观看av| 国产极品粉嫩免费观看在线| 一本色道久久久久久精品综合| 亚洲欧美日韩高清在线视频 | 搡老乐熟女国产| 99热全是精品| 一个人免费看片子| 欧美日韩中文字幕国产精品一区二区三区 | 精品一区在线观看国产| 9191精品国产免费久久| 蜜桃国产av成人99| 考比视频在线观看| 久久精品国产亚洲av高清一级| 精品久久蜜臀av无| 国产精品久久久久成人av| 精品亚洲成a人片在线观看| 亚洲精品美女久久久久99蜜臀| 亚洲人成电影观看| av天堂久久9| 少妇 在线观看| 国产主播在线观看一区二区| 免费看十八禁软件| 午夜久久久在线观看| 色老头精品视频在线观看| 天天影视国产精品| 99精国产麻豆久久婷婷| 中文字幕人妻熟女乱码| 久久人人爽av亚洲精品天堂| 久久人人爽人人片av| 日韩电影二区| 91成年电影在线观看| 亚洲色图 男人天堂 中文字幕| 法律面前人人平等表现在哪些方面 | 黄色视频在线播放观看不卡| 欧美在线黄色| 亚洲九九香蕉| av片东京热男人的天堂| 一级毛片女人18水好多| 一边摸一边做爽爽视频免费| 日韩欧美一区二区三区在线观看 | 18禁黄网站禁片午夜丰满| 丰满迷人的少妇在线观看| 亚洲五月色婷婷综合| 久久精品aⅴ一区二区三区四区| 免费女性裸体啪啪无遮挡网站| 99精品久久久久人妻精品| 脱女人内裤的视频| 超碰97精品在线观看| 少妇被粗大的猛进出69影院| 国产精品久久久久久人妻精品电影 | 欧美变态另类bdsm刘玥| 伊人亚洲综合成人网| 国产成人啪精品午夜网站| 国产精品99久久99久久久不卡| 日韩精品免费视频一区二区三区| 满18在线观看网站| 免费久久久久久久精品成人欧美视频| 亚洲精品一卡2卡三卡4卡5卡 | 性色av乱码一区二区三区2| 亚洲国产av影院在线观看| 一级毛片精品| 午夜免费成人在线视频| 91av网站免费观看| 两人在一起打扑克的视频| 热re99久久国产66热| 国产精品影院久久| 国产亚洲精品久久久久5区| 国产成人欧美在线观看 | 久久人妻熟女aⅴ| 一级片'在线观看视频| 国产亚洲一区二区精品| 十八禁网站网址无遮挡| 欧美日韩福利视频一区二区| 久久天堂一区二区三区四区| 欧美激情久久久久久爽电影 | 久久免费观看电影| 午夜老司机福利片| 国精品久久久久久国模美| 老熟女久久久| 午夜福利在线免费观看网站| 日韩制服骚丝袜av| 老司机深夜福利视频在线观看 | 中文字幕人妻熟女乱码| 午夜日韩欧美国产| 欧美精品高潮呻吟av久久| cao死你这个sao货| 免费av中文字幕在线| 日韩电影二区| 精品久久久久久久毛片微露脸 | 午夜免费鲁丝| 老熟女久久久| 国产一区二区三区在线臀色熟女 | 免费久久久久久久精品成人欧美视频| 亚洲三区欧美一区| 在线精品无人区一区二区三| 女警被强在线播放| 女性被躁到高潮视频| 久久天堂一区二区三区四区| 亚洲激情五月婷婷啪啪| 一区二区三区四区激情视频| 久久人人97超碰香蕉20202| 最近中文字幕2019免费版| 18禁裸乳无遮挡动漫免费视频| 狂野欧美激情性bbbbbb| 日本精品一区二区三区蜜桃| 国产国语露脸激情在线看| 高清欧美精品videossex| 欧美变态另类bdsm刘玥| av福利片在线| 天天躁日日躁夜夜躁夜夜| 日本欧美视频一区| 最新的欧美精品一区二区| 少妇猛男粗大的猛烈进出视频| 丝袜喷水一区| 欧美亚洲 丝袜 人妻 在线| 国产成人精品久久二区二区免费| 亚洲国产看品久久| 超色免费av| 精品卡一卡二卡四卡免费| 精品视频人人做人人爽| 波多野结衣一区麻豆| av电影中文网址| 黄片小视频在线播放| 99热国产这里只有精品6| 我要看黄色一级片免费的| 国产一卡二卡三卡精品| 欧美日韩一级在线毛片| 国产老妇伦熟女老妇高清| 国产免费av片在线观看野外av| 99精品久久久久人妻精品| 国产在视频线精品| 久久天堂一区二区三区四区| 视频区图区小说| 手机成人av网站| 一二三四社区在线视频社区8| 少妇的丰满在线观看| 国产精品久久久久久精品古装| 国产免费现黄频在线看| √禁漫天堂资源中文www| 国产深夜福利视频在线观看| 香蕉丝袜av| 首页视频小说图片口味搜索| 国产精品一二三区在线看| 亚洲性夜色夜夜综合| 国产人伦9x9x在线观看| 99久久人妻综合| 嫁个100分男人电影在线观看| 99久久99久久久精品蜜桃| 国产免费av片在线观看野外av| 日本wwww免费看| 亚洲中文字幕日韩| 精品少妇黑人巨大在线播放| 最新在线观看一区二区三区| 国产日韩欧美视频二区| 黑丝袜美女国产一区| 亚洲情色 制服丝袜| 国产在线观看jvid| 一区在线观看完整版| 欧美国产精品va在线观看不卡| 日本精品一区二区三区蜜桃| 精品欧美一区二区三区在线| 亚洲国产欧美日韩在线播放| 久久精品国产亚洲av高清一级| 午夜免费鲁丝| 男女无遮挡免费网站观看| 国产国语露脸激情在线看| 久热爱精品视频在线9| 欧美日韩亚洲综合一区二区三区_| 欧美少妇被猛烈插入视频| 新久久久久国产一级毛片| 国产黄频视频在线观看| 高清欧美精品videossex| 国产欧美亚洲国产| svipshipincom国产片| 精品人妻一区二区三区麻豆| 欧美日韩成人在线一区二区| 老司机深夜福利视频在线观看 | 午夜久久久在线观看| 日韩制服丝袜自拍偷拍| 免费看十八禁软件| 欧美乱码精品一区二区三区| 亚洲久久久国产精品| 丰满人妻熟妇乱又伦精品不卡| 亚洲少妇的诱惑av| 麻豆国产av国片精品| 动漫黄色视频在线观看| 久久久久国内视频| 啦啦啦视频在线资源免费观看| 1024香蕉在线观看| 国产男女内射视频| 人人妻,人人澡人人爽秒播| tocl精华| 麻豆乱淫一区二区| 亚洲国产欧美一区二区综合| netflix在线观看网站| 日本黄色日本黄色录像| 国产精品一区二区在线不卡| 青春草视频在线免费观看| 精品一区二区三区四区五区乱码| 韩国高清视频一区二区三区| 老鸭窝网址在线观看| 亚洲va日本ⅴa欧美va伊人久久 | 在线观看免费日韩欧美大片| 精品福利观看| 性色av乱码一区二区三区2| 超碰成人久久| 一个人免费在线观看的高清视频 | 韩国高清视频一区二区三区| 啪啪无遮挡十八禁网站| 精品久久久精品久久久| 五月开心婷婷网| 国产成人av激情在线播放| 亚洲精品久久成人aⅴ小说| 777米奇影视久久| 91成人精品电影| 制服诱惑二区| 久久国产精品影院| 亚洲av日韩在线播放| 不卡av一区二区三区| 中文字幕色久视频| 老熟女久久久| 欧美日韩一级在线毛片| 啦啦啦啦在线视频资源| 丝袜人妻中文字幕| 国产精品1区2区在线观看. | 国产97色在线日韩免费| 丁香六月天网| 国产成人av教育| 女人被躁到高潮嗷嗷叫费观| 免费观看a级毛片全部| 一级a爱视频在线免费观看| 日本a在线网址| 咕卡用的链子| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲色图综合在线观看| 欧美午夜高清在线| 法律面前人人平等表现在哪些方面 | 在线天堂中文资源库| 青青草视频在线视频观看| 久久久精品免费免费高清| 国产男人的电影天堂91| 涩涩av久久男人的天堂| 性少妇av在线| 99久久人妻综合| 精品久久久久久电影网| 亚洲伊人色综图| 久久久精品国产亚洲av高清涩受| 久久久久久久大尺度免费视频| 欧美日韩视频精品一区| 国产视频一区二区在线看| 人成视频在线观看免费观看| 在线十欧美十亚洲十日本专区| 久久人妻福利社区极品人妻图片| 国产免费视频播放在线视频| av有码第一页| 国产av国产精品国产| 岛国毛片在线播放| 窝窝影院91人妻| 成年人免费黄色播放视频| 久久久久久久精品精品| 99re6热这里在线精品视频| 青青草视频在线视频观看| 亚洲av美国av| 两个人看的免费小视频| 1024香蕉在线观看| av天堂久久9| 黄片小视频在线播放| 精品一区在线观看国产| av视频免费观看在线观看| 欧美另类亚洲清纯唯美| 久久狼人影院| av网站在线播放免费| 巨乳人妻的诱惑在线观看| 夜夜夜夜夜久久久久| 精品人妻一区二区三区麻豆| 免费不卡黄色视频| av线在线观看网站| 三级毛片av免费| 国产在线一区二区三区精| 亚洲欧美成人综合另类久久久| 成人av一区二区三区在线看 | 男人爽女人下面视频在线观看| 男女之事视频高清在线观看| 岛国在线观看网站| 大型av网站在线播放| 国产精品一区二区精品视频观看| 超色免费av| 每晚都被弄得嗷嗷叫到高潮| 18禁黄网站禁片午夜丰满| av在线app专区| 狂野欧美激情性xxxx| 午夜影院在线不卡| 大片电影免费在线观看免费| 久久亚洲精品不卡| 久久午夜综合久久蜜桃| 日韩,欧美,国产一区二区三区| 国精品久久久久久国模美| 亚洲人成77777在线视频| 亚洲欧美清纯卡通| 亚洲av欧美aⅴ国产| 老熟妇乱子伦视频在线观看 | 亚洲成国产人片在线观看| 美女中出高潮动态图| 高清av免费在线| h视频一区二区三区| 制服人妻中文乱码| 欧美日本中文国产一区发布| 中文字幕高清在线视频|