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

    國際糧食貿(mào)易網(wǎng)絡(luò)魯棒性分析

    2021-08-12 17:51:42董志良楊巧然
    當(dāng)代經(jīng)濟(jì)管理 2021年6期
    關(guān)鍵詞:復(fù)雜網(wǎng)絡(luò)糧食安全魯棒性

    董志良 楊巧然

    [摘 要] 目前國際新冠疫情形式嚴(yán)峻,各國對糧食供給的安全性關(guān)注度較高,為了解國際糧食貿(mào)易的基本特征,研究國際糧食貿(mào)易是否安全穩(wěn)定,文章采用復(fù)雜網(wǎng)絡(luò)理論研究國際糧食貿(mào)易網(wǎng)絡(luò)的網(wǎng)絡(luò)特征并分析其魯棒性。首先,利用聯(lián)合國統(tǒng)計(jì)司發(fā)布的2019年糧食貿(mào)易數(shù)據(jù)構(gòu)建了全球糧食貿(mào)易網(wǎng)絡(luò),通過計(jì)算網(wǎng)絡(luò)特征參數(shù)分析國際糧食貿(mào)易的特性;然后,評價(jià)貿(mào)易網(wǎng)絡(luò)中的重要節(jié)點(diǎn),選擇魯棒性分析指標(biāo);最后,計(jì)算不同攻擊方式下的魯棒性指標(biāo),分析不同攻擊下的糧食貿(mào)易網(wǎng)絡(luò)魯棒性。結(jié)果表明:糧食貿(mào)易網(wǎng)絡(luò)的節(jié)點(diǎn)在面對隨機(jī)攻擊時(shí)擁有較好的魯棒性,在目標(biāo)攻擊下表現(xiàn)出脆弱性。通過研究,可以為國際糧食貿(mào)易參與國糧食供給安全提供參考。

    [關(guān)鍵詞]糧食貿(mào)易;復(fù)雜網(wǎng)絡(luò);糧食安全;魯棒性

    [中圖分類號(hào)] F316.11[文獻(xiàn)標(biāo)識(shí)碼] A[文章編號(hào)]1673-0461(2021)06-0073-06

    一、引 言

    在經(jīng)濟(jì)全球化背景下,國際糧食貿(mào)易量呈現(xiàn)快速增長態(tài)勢,越來越多的國家參與到糧食進(jìn)出口的國際糧食貿(mào)易中來,糧食的自由貿(mào)易使糧食資源在全球范圍內(nèi)重新分配,糧食貿(mào)易已成為連接全球糧食資源豐富地區(qū)和匱乏地區(qū)的重要紐帶口,這不僅滿足了很多國家糧食供給缺失的問題,也強(qiáng)化了各國之間的經(jīng)濟(jì)聯(lián)系和相互依賴。由此,全球糧食安全正日漸依賴于一些糧食的國際貿(mào)易,全球糧食的穩(wěn)定供應(yīng)影響著世界糧食安全,糧食供給穩(wěn)定是關(guān)乎國計(jì)民生的頭等大事,但伴隨著新冠疫情在全球范圍內(nèi)的傳播,國際貿(mào)易的環(huán)境處于極其不穩(wěn)定的狀態(tài),國際糧食貿(mào)易也受到此次疫情的影響,因此,研究國際貿(mào)易網(wǎng)絡(luò)的魯棒性對分析世界糧食安全具有重要意義。

    糧食作為國家的戰(zhàn)略物資和人民的生活必需品,很多學(xué)者都對糧食進(jìn)行了不同方面的研究。其中,有些學(xué)者對糧食生產(chǎn)方面進(jìn)行研究,研究了糧食生產(chǎn)潛力、生產(chǎn)格局以及糧食生產(chǎn)技術(shù)等一系列問題[1-3],有些學(xué)者對糧食流通方面展開研究,包括糧食流通體制和糧食流通中的各種問題[4-5],糧食貿(mào)易作為糧食流通的一個(gè)重要環(huán)節(jié),是被學(xué)者廣泛研究的一個(gè)重要方向。封志明等[6]基于糧食生產(chǎn)與貿(mào)易數(shù)據(jù),利用集中化指數(shù)與回歸分析法,對中國-東盟糧食生產(chǎn)與貿(mào)易的格局和貿(mào)易演變進(jìn)行分析,Aparicio等 [7]通過揭示小麥、玉米和大米貿(mào)易結(jié)構(gòu)特征研究其市場運(yùn)作情況,分析了貿(mào)易雙方的進(jìn)出口流動(dòng)狀況。在糧食貿(mào)易中,大家廣泛關(guān)注的是在貿(mào)易中能否保證糧食供給穩(wěn)定,從而保障世界各國的糧食安全,對此,國內(nèi)外很多學(xué)者展開了研究。其中,余志剛等[8]通過分析國際糧食貿(mào)易的現(xiàn)狀、趨勢以及中國糧食貿(mào)易和國際糧食貿(mào)易的關(guān)系,分析了國際糧食貿(mào)易市場對中國糧食供給安全的影響;張晶等[9]運(yùn)用Johansen協(xié)整檢驗(yàn)與P-T模型對3個(gè)大米主產(chǎn)國大米價(jià)格變化的共同因素進(jìn)行了分析,從而分析價(jià)格對糧食安全的作用機(jī)制;蘇珊珊[10]通過利用1990—2009年的跨國面板數(shù)據(jù),建立固定效應(yīng)模型研究糧食安全和糧食貿(mào)易的關(guān)系;Wang等[11]分析了糧食危機(jī)后糧食供給政策的變化和根源,并提出我國在糧食危機(jī)后可以實(shí)行的措施。

    復(fù)雜網(wǎng)絡(luò)是研究許多現(xiàn)實(shí)世界系統(tǒng)的一個(gè)有價(jià)值的工具,現(xiàn)實(shí)世界中的許多系統(tǒng)都可以用復(fù)雜網(wǎng)絡(luò)來描述[12],利用復(fù)雜網(wǎng)絡(luò)所具有的結(jié)構(gòu)特征和網(wǎng)絡(luò)功能來描述現(xiàn)實(shí)系統(tǒng)的特征和結(jié)構(gòu),目前被學(xué)者們廣泛應(yīng)用于鐵路安全和原油鋼鐵等大宗商品的貿(mào)易網(wǎng)絡(luò)等各個(gè)領(lǐng)域[13]。另一方面,國內(nèi)外也有學(xué)者將這種方法應(yīng)用到糧食研究中[14],學(xué)者們通過將貿(mào)易雙方作為節(jié)點(diǎn),貿(mào)易流向作為邊,構(gòu)建各種糧食作物的貿(mào)易網(wǎng)絡(luò),探究糧食網(wǎng)絡(luò)的結(jié)構(gòu)特征,以及糧食貿(mào)易動(dòng)態(tài)演化的趨勢,但是目前還沒有人利用復(fù)雜網(wǎng)絡(luò)研究國際糧食貿(mào)易網(wǎng)絡(luò)的魯棒性。

    因此,本文將全球糧食貿(mào)易體系抽象成一個(gè)復(fù)雜網(wǎng)絡(luò),以參與貿(mào)易的各經(jīng)濟(jì)體為節(jié)點(diǎn),以各經(jīng)濟(jì)體之間的貿(mào)易關(guān)系為邊,通過探究拓?fù)渚W(wǎng)絡(luò)的特性分析國際糧食貿(mào)易的特征,以及通過對網(wǎng)絡(luò)進(jìn)行模擬攻擊分析國際糧食貿(mào)易在受到目標(biāo)中斷以及隨機(jī)中斷時(shí)的安全問題。

    二、數(shù)據(jù)與方法

    (一)數(shù)據(jù)來源

    大米、小麥、玉米是世界三大糧食作物,雖然中國國家統(tǒng)計(jì)局對糧食的定義中包含豆類,但大量的國際糧食安全研究和糧農(nóng)組織、IMF等數(shù)據(jù)統(tǒng)計(jì)中,一般未將大豆統(tǒng)計(jì)在內(nèi)[15],本文研究的是國際貿(mào)易中糧食的進(jìn)出口貿(mào)易情況,故選擇小麥、玉米和大米進(jìn)行研究。數(shù)據(jù)來源為聯(lián)合國貿(mào)易統(tǒng)計(jì)數(shù)據(jù)庫,由于本文主要研究糧食網(wǎng)絡(luò)面對攻擊的魯棒性,因此只選擇2019年的糧食數(shù)據(jù)來構(gòu)建網(wǎng)絡(luò)進(jìn)行研究。

    (二)方法

    以貿(mào)易國家作為網(wǎng)絡(luò)節(jié)點(diǎn),以國家之間的貿(mào)易關(guān)系為邊,貿(mào)易流向?yàn)榫W(wǎng)絡(luò)中的箭頭方向,構(gòu)建復(fù)雜網(wǎng)絡(luò)模型。對全球糧食貿(mào)易網(wǎng)絡(luò)F進(jìn)行如下定義:F=(M,N,W,W)。其中,M代表所有節(jié)點(diǎn)(國家)的集合,N為所有邊(兩國間貿(mào)易聯(lián)系)的集合,W代表所有的節(jié)點(diǎn)屬性(國家的貿(mào)易伙伴數(shù)量)的集合,W表示所有邊屬性(兩國之間的貿(mào)易總量)的集合。分別構(gòu)建小麥、玉米以及大米的國際貿(mào)易網(wǎng)絡(luò)。

    這樣便將糧食貿(mào)易網(wǎng)絡(luò)轉(zhuǎn)換成一個(gè)有向加權(quán)連通圖,利用復(fù)雜網(wǎng)絡(luò)的特征參數(shù),從度與度分布、聚類系數(shù)、平均路徑長度等方面對糧食貿(mào)易網(wǎng)絡(luò)進(jìn)行說明。

    1.度與度分布

    節(jié)點(diǎn)度(K)是指在拓?fù)渚W(wǎng)絡(luò)中與某個(gè)特定節(jié)點(diǎn)直接相連的節(jié)點(diǎn)數(shù)量。在有向網(wǎng)絡(luò)中,度又分為出度和入度。在糧食貿(mào)易網(wǎng)絡(luò)中,出度是一個(gè)國家的所有出口國數(shù)量,入度則是一個(gè)國家對應(yīng)的所有進(jìn)口國的數(shù)量,度值大小代表這個(gè)國家在貿(mào)易網(wǎng)絡(luò)中與之發(fā)生貿(mào)易的國家數(shù)量多少,節(jié)點(diǎn)度越大,表示與該國發(fā)生貿(mào)易的國家數(shù)量越多,該國在貿(mào)易網(wǎng)絡(luò)中的地位也就相對重要。度的計(jì)算公式如下:

    Ki=∑nj=1eij(1)

    度分布是可以用函數(shù)P(K)表示,P(K)是指在網(wǎng)絡(luò)中隨機(jī)選一個(gè)節(jié)點(diǎn),它的度恰好為K的概率。

    2.網(wǎng)絡(luò)密度

    網(wǎng)絡(luò)密度(Q)反映的是一個(gè)網(wǎng)絡(luò)的緊密程度。在一個(gè)有向網(wǎng)絡(luò)中是n(n-1),網(wǎng)絡(luò)密度為所有節(jié)點(diǎn)均兩兩相連時(shí)的連邊數(shù)量和網(wǎng)絡(luò)中實(shí)際存在的連邊數(shù)量之比,在糧食貿(mào)易網(wǎng)絡(luò)中,網(wǎng)絡(luò)密度越大,代表國際糧食貿(mào)易越活躍,國與國之間的貿(mào)易往來越密切。網(wǎng)絡(luò)密度的計(jì)算公式如下:

    Q=mn(n-1)(2)

    其中,m是網(wǎng)絡(luò)實(shí)際存在的連邊數(shù)量,n(n-1)是網(wǎng)絡(luò)中可能的最大連邊數(shù)量。

    3.聚類系數(shù)

    聚類系數(shù)(G)反映的是一個(gè)大網(wǎng)絡(luò)中各集聚的小網(wǎng)絡(luò)分布和相互聯(lián)系的狀況。在一個(gè)網(wǎng)絡(luò)中,如果兩個(gè)節(jié)點(diǎn)分別與一個(gè)相同節(jié)點(diǎn)直接相連,那么這兩個(gè)節(jié)點(diǎn)很有可能直接相連,稱為復(fù)雜網(wǎng)絡(luò)的聚類特性。在糧食貿(mào)易網(wǎng)絡(luò)中,與同一個(gè)國家存在貿(mào)易關(guān)系的兩個(gè)國家也可能存在貿(mào)易關(guān)系。聚類系數(shù)越大,節(jié)點(diǎn)周圍的節(jié)點(diǎn)聚集程度越高,整個(gè)貿(mào)易更加緊密,代表國際糧食貿(mào)易越密切。聚類系數(shù)的計(jì)算公式為:

    Gi=2NiKi(Ki-1)(3)

    其中,Gi表示節(jié)點(diǎn)i的聚類系數(shù),Ki表示節(jié)點(diǎn)i的節(jié)點(diǎn)度,Ni表示節(jié)點(diǎn)i與所有鄰居節(jié)點(diǎn)間的邊的數(shù)量。

    4.平均路徑長度

    平均路徑長度(L)是貿(mào)易網(wǎng)絡(luò)評估商品貿(mào)易效率的一個(gè)重要指標(biāo)。在一個(gè)網(wǎng)絡(luò)中,一個(gè)節(jié)點(diǎn)可以通過一條或者多條不同的路徑到達(dá)另一個(gè)節(jié)點(diǎn),連接兩個(gè)節(jié)點(diǎn)的最少的邊的數(shù)量稱為節(jié)點(diǎn)間的最短路徑(dij),平均路徑長度是所有相連節(jié)點(diǎn)間最短路徑的平均值。網(wǎng)絡(luò)中任意兩個(gè)節(jié)點(diǎn)之間距離的最大值稱為網(wǎng)絡(luò)的直徑(D)。在糧食貿(mào)易網(wǎng)絡(luò)中,平均路徑長度和網(wǎng)絡(luò)直徑越小,說明全球糧食貿(mào)易交往越緊密。平均路徑長度(L)和網(wǎng)絡(luò)直徑(D)的計(jì)算公式如下:

    L=1n(n-1)∑i∑jdij(4)

    D=maxdij(5)

    5.網(wǎng)絡(luò)效率

    在評估網(wǎng)絡(luò)的魯棒性時(shí),網(wǎng)絡(luò)效率(E)經(jīng)常被用來作為評價(jià)指標(biāo),它能夠很好地反映網(wǎng)絡(luò)的拓?fù)湫阅?,尤其是?jié)點(diǎn)之間的連通性及整體效率,被學(xué)者們應(yīng)用于交通網(wǎng)絡(luò)以及能源網(wǎng)絡(luò)的魯棒性評價(jià)中。因此,本文選取網(wǎng)絡(luò)效率作為糧食貿(mào)易網(wǎng)絡(luò)的評價(jià)指標(biāo),用來衡量糧食網(wǎng)絡(luò)中貿(mào)易流通效率,網(wǎng)絡(luò)效率越高,則全球糧食貿(mào)易聯(lián)系越緊密,代表國際糧食的供應(yīng)和需求越能得以滿足;當(dāng)網(wǎng)絡(luò)效率很低時(shí),則說明網(wǎng)絡(luò)的連通性和傳輸效率都很低,這在糧食貿(mào)易網(wǎng)絡(luò)中表現(xiàn)為國與國之間的貿(mào)易關(guān)系數(shù)量很少,不能維持大部分國家的糧食需求,此時(shí)全球糧食安全受到威脅。網(wǎng)絡(luò)效率是指網(wǎng)絡(luò)中所有節(jié)點(diǎn)對之間距離倒數(shù)之和的平均值。計(jì)算公式為:

    E=1n(n-1)∑1dij(6)

    三、網(wǎng)絡(luò)整體特征及其重要節(jié)點(diǎn)分析

    (一)糧食網(wǎng)絡(luò)整體特征

    利用Gephi計(jì)算出拓?fù)渚W(wǎng)絡(luò)的特征參數(shù),結(jié)果如表1所示,全球大部分的國家或地區(qū)都參與國際間的糧食貿(mào)易,其中,共有228個(gè)國家或地區(qū)參與大米的貿(mào)易,214個(gè)國家或地區(qū)參與玉米貿(mào)易,199個(gè)國家或地區(qū)參與小麥貿(mào)易;3種糧食網(wǎng)絡(luò)的復(fù)雜程度都很高,其中,大米貿(mào)易網(wǎng)絡(luò)中貿(mào)易聯(lián)系達(dá)3 875條,平均1個(gè)國家與17個(gè)國家都進(jìn)行大米貿(mào)易,與12—13個(gè)國家進(jìn)行玉米貿(mào)易,與10—11個(gè)國家進(jìn)行小麥貿(mào)易。3種糧食貿(mào)易的網(wǎng)絡(luò)密度在0.05—0.075之間,成員貿(mào)易結(jié)構(gòu)相對松散。網(wǎng)絡(luò)的平均路徑長度在2—3之間,網(wǎng)絡(luò)的直徑均在5—6之間,表明在糧食貿(mào)易網(wǎng)絡(luò)中,有少數(shù)國家或地區(qū)在網(wǎng)絡(luò)中擔(dān)任橋梁的作用,他們在糧食貿(mào)易網(wǎng)絡(luò)中處于重要地位。

    在復(fù)雜網(wǎng)絡(luò)中,如果少部分的節(jié)點(diǎn)度值較大,大部分的節(jié)點(diǎn)度值較小,即少部分國家或地區(qū)在網(wǎng)絡(luò)中和大多數(shù)國家或地區(qū)均存在貿(mào)易關(guān)系,那么網(wǎng)絡(luò)節(jié)點(diǎn)的度分布服從冪律分布,具有冪律分布特征的網(wǎng)絡(luò)屬于無標(biāo)度網(wǎng)絡(luò),可以利用無標(biāo)度網(wǎng)絡(luò)的一些基本特征對糧食貿(mào)易網(wǎng)絡(luò)進(jìn)行解釋。為了驗(yàn)證糧食網(wǎng)絡(luò)是否服從冪律分布,對國際糧食貿(mào)易網(wǎng)絡(luò)中的節(jié)點(diǎn)度及累積度分布數(shù)據(jù)在雙對數(shù)坐標(biāo)系中進(jìn)行擬合,結(jié)果如圖1所示。

    由此可知,3種糧食貿(mào)易網(wǎng)絡(luò)的度分布均服從冪律分布,R2均在0.8上下,擬合效果較好,進(jìn)而可以知道糧食貿(mào)易網(wǎng)絡(luò)是無標(biāo)度網(wǎng)絡(luò)。在Albert的研究結(jié)論中[12]可以知道,無標(biāo)度網(wǎng)絡(luò)在受到隨機(jī)攻擊時(shí)可以表現(xiàn)出較好的魯棒性,而對目標(biāo)攻擊表現(xiàn)出脆弱性,進(jìn)一步為本文提供了更加可靠的理論依據(jù),后文將進(jìn)行網(wǎng)絡(luò)的仿真攻擊來分析糧食網(wǎng)絡(luò)的魯棒性。

    (二)重要節(jié)點(diǎn)分析

    對于國際糧食貿(mào)易網(wǎng)絡(luò)而言,一個(gè)節(jié)點(diǎn)的重要性并不完全取決于節(jié)點(diǎn)連接度的數(shù)值大小,還與邊的權(quán)重大小有關(guān),它代表著貿(mào)易雙方的貿(mào)易數(shù)量,貿(mào)易數(shù)量越大,說明它在網(wǎng)絡(luò)中的地位也越重要,同時(shí),對于無標(biāo)度網(wǎng)絡(luò)而言,中介中心性和特征向量中心性同樣是衡量節(jié)點(diǎn)重要性的關(guān)鍵指標(biāo)。中介中心性在糧食貿(mào)易網(wǎng)絡(luò)中衡量了一個(gè)國家或地區(qū)作為媒介者的能力,也就是在其它國家或地區(qū)之間擔(dān)任貿(mào)易中轉(zhuǎn)的角色,在網(wǎng)絡(luò)中,中介中心性越高,則它擔(dān)任中轉(zhuǎn)的角色越多,這樣的國家或地區(qū)在糧食貿(mào)易網(wǎng)中就更加重要。特征向量中心性是根據(jù)網(wǎng)絡(luò)中鄰居節(jié)點(diǎn)的價(jià)值來衡量該節(jié)點(diǎn)的價(jià)值,在糧食貿(mào)易網(wǎng)中,與一個(gè)國家或地區(qū)進(jìn)行貿(mào)易的所有國家或地區(qū)在網(wǎng)絡(luò)中是重要的,那么該國家或地區(qū)的特征向量中心性就會(huì)越高,這個(gè)國家也越重要。為了能綜合體現(xiàn)網(wǎng)絡(luò)中節(jié)點(diǎn)的重要程度,借助文獻(xiàn)[16]中提出的節(jié)點(diǎn)重要性評價(jià)模型,構(gòu)建了如下的評價(jià)模型:

    Yi=KiPi(a1Xi1+a2Xi2)(7)

    其中,Yi是節(jié)點(diǎn)i的重要性指標(biāo);a1和a2為權(quán)重系數(shù),且滿足a1+a2=1;Ki為節(jié)點(diǎn)i的度值;Pi為節(jié)點(diǎn)i的進(jìn)出口總額;Xi1是節(jié)點(diǎn)i的中介中心性;Xi2為節(jié)點(diǎn)i的特征向量中心性。

    本文確定權(quán)重系數(shù)a1、a2取值的方法如下:a1和a2的取值在[0,1]間變化,分別計(jì)算11種不同的a1、a2取值下,網(wǎng)絡(luò)中節(jié)點(diǎn)重要性的排名情況,根據(jù)該排名進(jìn)行節(jié)點(diǎn)的目標(biāo)攻擊試驗(yàn),分別計(jì)算每次攻擊后的網(wǎng)絡(luò)效率,為最大程度反映節(jié)點(diǎn)的重要度,與初始網(wǎng)絡(luò)效率值相比較,網(wǎng)絡(luò)效率變化最大的一組值即為a1、a2的取值。經(jīng)過試驗(yàn)得到:當(dāng)a1=1,a2=0時(shí),3個(gè)網(wǎng)絡(luò)的網(wǎng)絡(luò)效率值降到0.0001的速度都是最快的,故選用該權(quán)重組合作為評價(jià)指標(biāo)的權(quán)重系數(shù)。根據(jù)計(jì)算結(jié)果,不同糧食網(wǎng)絡(luò)下的節(jié)點(diǎn)重要性排名前10的國家或地區(qū)如表2所示。

    美國、歐盟27國在3種糧食作物貿(mào)易中均為重要的貿(mào)易國家或地區(qū),法國、德國以及巴西等國家在兩種糧食作物的進(jìn)出口中處于重要的地位,還有部分國家如印度、泰國是某一種糧食作物的生產(chǎn)大國。在全球糧食貿(mào)易網(wǎng)絡(luò)中,幾個(gè)核心的貿(mào)易國家?guī)缀跖c世界上的所有國家都存在著貿(mào)易關(guān)系,這些國家不僅擁有大量的貿(mào)易合作伙伴,并且糧食交易量總額也在所有貿(mào)易國家中位居前列,它們在國際糧食貿(mào)易中起著至關(guān)重要的作用。

    四、結(jié)果分析

    (一)魯棒性分析

    復(fù)雜網(wǎng)絡(luò)的節(jié)點(diǎn)攻擊是指當(dāng)某一節(jié)點(diǎn)受到攻擊時(shí),移除該節(jié)點(diǎn)和與該節(jié)點(diǎn)相連的所有邊,形成一個(gè)新的拓?fù)渚W(wǎng)絡(luò),節(jié)點(diǎn)攻擊方式通常有隨機(jī)攻擊和目標(biāo)攻擊兩種。隨機(jī)攻擊是指以某種概率隨機(jī)地破壞網(wǎng)絡(luò)中的節(jié)點(diǎn),與網(wǎng)絡(luò)節(jié)點(diǎn)在網(wǎng)絡(luò)中的位置和重要性無關(guān),在國際貿(mào)易網(wǎng)絡(luò)中,隨機(jī)攻擊表現(xiàn)為某些國家或地區(qū)可能受到天氣、自然災(zāi)害或國際經(jīng)貿(mào)動(dòng)蕩等原因?qū)е录Z食在某一年無法進(jìn)行出口貿(mào)易,這種情況發(fā)生伴隨著不確定性。目標(biāo)攻擊是指對網(wǎng)絡(luò)中的節(jié)點(diǎn)進(jìn)行針對性地破壞,這種攻擊是按照節(jié)點(diǎn)的重要程度進(jìn)行破壞,對網(wǎng)絡(luò)造成的危害往往更大,在糧食貿(mào)易網(wǎng)絡(luò)中可能體現(xiàn)在某一個(gè)國家或地區(qū)的進(jìn)出口政策發(fā)生巨大改變,使該國的糧食供應(yīng)或需求發(fā)生中斷。在新冠疫情影響下,部分國家或地區(qū)出臺(tái)的有關(guān)停止國際糧食貿(mào)易的政策可以看成是目標(biāo)攻擊造成的供給中斷。

    復(fù)雜網(wǎng)絡(luò)的邊攻擊是指當(dāng)某一條連邊受到攻擊時(shí),移除該邊,保持該節(jié)點(diǎn)和該節(jié)點(diǎn)與其它節(jié)點(diǎn)的連邊不變,形成一個(gè)新的拓?fù)渚W(wǎng)絡(luò)。在糧食貿(mào)易中,兩個(gè)國家或地區(qū)之間的貿(mào)易中斷都可以體現(xiàn)為網(wǎng)絡(luò)的邊攻擊。

    本文在進(jìn)行對網(wǎng)絡(luò)的節(jié)點(diǎn)仿真破壞時(shí),采取度值攻擊、按節(jié)點(diǎn)重要性進(jìn)行攻擊以及隨機(jī)攻擊3種方式。在進(jìn)行目標(biāo)攻擊時(shí),分別按照節(jié)點(diǎn)的度值和節(jié)點(diǎn)重要性依次刪除網(wǎng)絡(luò)中的節(jié)點(diǎn)并觀察整個(gè)網(wǎng)絡(luò)的網(wǎng)絡(luò)效率變化情況。隨機(jī)攻擊是利用Matlab程序生成隨機(jī)序列,隨機(jī)攻擊網(wǎng)絡(luò)中的節(jié)點(diǎn),由于隨機(jī)攻擊的不確定性較大,本文的隨機(jī)攻擊取10次隨機(jī)攻擊結(jié)果的平均值進(jìn)行研究。圖2—圖4分別為小麥、玉米以及大米網(wǎng)絡(luò)的仿真結(jié)果。在進(jìn)行邊的仿真模擬時(shí),采用根據(jù)邊權(quán)重開展目標(biāo)攻擊的方式進(jìn)行。3種糧食網(wǎng)絡(luò)的仿真結(jié)果如圖5所示。

    對節(jié)點(diǎn)進(jìn)行攻擊時(shí),目標(biāo)攻擊和隨機(jī)攻擊對節(jié)點(diǎn)的影響程度不同,與目標(biāo)攻擊方式相比,面對隨機(jī)攻擊時(shí),網(wǎng)絡(luò)效率下降較為緩慢,在網(wǎng)絡(luò)效率均下降至0.1以下時(shí),目標(biāo)攻擊所攻擊的糧食網(wǎng)絡(luò)節(jié)點(diǎn)數(shù)量僅為隨機(jī)攻擊的節(jié)點(diǎn)數(shù)量的25%,糧食貿(mào)易網(wǎng)絡(luò)在面對隨機(jī)攻擊時(shí)具有較好的魯棒性。在全球小麥貿(mào)易網(wǎng)絡(luò)中,當(dāng)攻擊的節(jié)點(diǎn)數(shù)量占網(wǎng)絡(luò)中節(jié)點(diǎn)的7%時(shí),網(wǎng)絡(luò)效率已經(jīng)下降到0.1以下;當(dāng)攻擊節(jié)點(diǎn)數(shù)占所有節(jié)點(diǎn)的20%時(shí),網(wǎng)絡(luò)效率下降至0.01以下,網(wǎng)絡(luò)間的連通基本崩潰。說明核心貿(mào)易國在貿(mào)易中的作用不可或缺,當(dāng)這些國家的糧食安全出現(xiàn)問題時(shí),將會(huì)影響全球小麥貿(mào)易的正常運(yùn)行。在玉米貿(mào)易網(wǎng)絡(luò)和大米貿(mào)易網(wǎng)絡(luò)中,當(dāng)攻擊的節(jié)點(diǎn)數(shù)量在各自節(jié)點(diǎn)數(shù)量12%時(shí),網(wǎng)絡(luò)效率下降為0.1以下。3種糧食貿(mào)易網(wǎng)絡(luò)在目標(biāo)攻擊下,根據(jù)度值進(jìn)行攻擊和根據(jù)節(jié)點(diǎn)重要性程度進(jìn)行攻擊的網(wǎng)絡(luò)效率變化情況相差不大。這是因?yàn)樵趯?jié)點(diǎn)進(jìn)行度值和節(jié)點(diǎn)重要性排序時(shí),各個(gè)國家的排序情況并未有較大的變化,同時(shí)說明,一個(gè)合作伙伴較多的國家,它的進(jìn)出口糧食數(shù)量總和也很龐大,同時(shí)它在網(wǎng)絡(luò)中所處的位置也靠近網(wǎng)絡(luò)中心。在對網(wǎng)絡(luò)中的邊進(jìn)行攻擊時(shí),網(wǎng)絡(luò)效率下降的非常緩慢,以大米網(wǎng)絡(luò)為例,在攻擊的連邊數(shù)量達(dá)到1 044條,即全世界的糧食貿(mào)易有1 000多對貿(mào)易聯(lián)系被切斷時(shí),整個(gè)網(wǎng)絡(luò)的網(wǎng)絡(luò)效率依然保持在0.25以上,可以看出糧食貿(mào)易網(wǎng)絡(luò)在面對邊攻擊時(shí)能夠具有較強(qiáng)的魯棒性。

    (二)世界糧食安全分析

    世界糧食貿(mào)易受到各種各樣不同因素的影響,貿(mào)易是否保持穩(wěn)定是不確定的。

    在新冠疫情的影響下,多個(gè)國家為了保障本國的糧食供應(yīng)穩(wěn)定,開始禁止或開始限制糧食的出口貿(mào)易,如小麥出口大國哈薩克斯坦曾頒布政策禁止出口面粉,大米的生產(chǎn)大國越南也曾于2020年3月份宣布停止大米的出口,俄羅斯作為小麥的核心貿(mào)易國之一也表示會(huì)根據(jù)未來情況來決定是否出臺(tái)出口禁令,并實(shí)行配額管理。這些糧食生產(chǎn)大國對出口的阻斷是政策上的阻斷,在糧食貿(mào)易網(wǎng)絡(luò)中表現(xiàn)為對節(jié)點(diǎn)進(jìn)行目標(biāo)攻擊,盡管這些國家目前已解除了這些政策,但是為了防范糧食出口大國的糧食中斷,對這些國家停止出口后可能給其它國家?guī)淼挠绊戇M(jìn)行分析。以越南為例,越南作為大米的出口大國之一,2019年出口的糧食超600萬噸,出口的國家或地區(qū)超100個(gè)。其中,菲律賓的進(jìn)口量高達(dá)270萬噸,加納、科特迪瓦、馬來西亞、中國、盧森堡、新加坡等的進(jìn)口量也超過100萬噸。一旦越南停止出口,中斷這些貿(mào)易聯(lián)系,將會(huì)對這些進(jìn)口國或地區(qū)產(chǎn)生很大的影響,并且對于進(jìn)口量不是很龐大但是糧食自給率低、糧食依賴于進(jìn)口的國家或地區(qū)來說,出口國的貿(mào)易中斷無疑會(huì)引發(fā)該國糧食安全問題。

    國際間影響糧食貿(mào)易的因素還有很多,自然災(zāi)害引起的糧食產(chǎn)量不足導(dǎo)致供給不足或者需求增加在網(wǎng)絡(luò)中表現(xiàn)為隨機(jī)攻擊,由國際間動(dòng)蕩的經(jīng)濟(jì)因素引起的糧食貿(mào)易受到?jīng)_擊可能表現(xiàn)為隨機(jī)攻擊或者是國與國之間的連邊攻擊,具體取決于該國面對沖擊所頒發(fā)的政策,需要根據(jù)政策的指向性對沖擊的類型進(jìn)行分析。

    (三)中國在糧食貿(mào)易中的地位

    在2019年國際糧食貿(mào)易網(wǎng)絡(luò)中,中國的小麥貿(mào)易伙伴有41個(gè),其中出口貿(mào)易伙伴占23個(gè),進(jìn)口貿(mào)易伙伴為18個(gè);玉米貿(mào)易伙伴63個(gè);大米貿(mào)易伙伴93個(gè),其中71個(gè)為出口合作關(guān)系,22個(gè)是進(jìn)口貿(mào)易關(guān)系。中國在全球糧食貿(mào)易網(wǎng)絡(luò)中扮演著十分重要的角色,每年在出口大量糧食的同時(shí)也在進(jìn)口大量糧食,在糧食出口上和世界大范圍的國家都有聯(lián)系,影響范圍較大,中國糧食保持供給和需求穩(wěn)定有利于世界糧食安全。

    五、結(jié) 論

    為了探討糧食貿(mào)易受到新冠疫情沖擊和可能面對的其它沖擊時(shí)能否保證供給安全,本文通過構(gòu)建3種糧食作物的貿(mào)易網(wǎng)絡(luò),分析國際糧食貿(mào)易網(wǎng)絡(luò)的基本特征,模擬糧食貿(mào)易網(wǎng)絡(luò)在受到攻擊時(shí)的魯棒性,從而分析國際糧食貿(mào)易在受到意外沖擊的情況下能否保持供應(yīng)安全,本文的研究結(jié)論如下:

    第一,國際糧食貿(mào)易參與國眾多,貿(mào)易聯(lián)系密切,所構(gòu)成的糧食貿(mào)易網(wǎng)絡(luò)均服從冪律分布,具有無標(biāo)度網(wǎng)絡(luò)的特征,即少部分的糧食大國與網(wǎng)絡(luò)中大部分的國家存在糧食貿(mào)易,這些國家對全球糧食供需影響顯著。

    第二,當(dāng)網(wǎng)絡(luò)中節(jié)點(diǎn)受到隨機(jī)攻擊時(shí),能保持較強(qiáng)的穩(wěn)定性,在遭受目標(biāo)攻擊時(shí)網(wǎng)絡(luò)更加脆弱,當(dāng)遭受攻擊的國家超過一定數(shù)量時(shí),可能會(huì)導(dǎo)致國際糧食貿(mào)易無法正常進(jìn)行。因此,世界上幾個(gè)國家的貿(mào)易中斷往往不會(huì)對國際糧食貿(mào)易安全產(chǎn)生大的危害,但是會(huì)對其貿(mào)易進(jìn)口合作國造成影響,當(dāng)退出國際貿(mào)易的國家達(dá)到一定數(shù)量后,就會(huì)導(dǎo)致整個(gè)糧食貿(mào)易網(wǎng)絡(luò)的崩潰,從而影響世界糧食安全。各個(gè)國家都應(yīng)關(guān)注環(huán)境和國際貿(mào)易中的動(dòng)蕩因素,識(shí)別潛在風(fēng)險(xiǎn),以便調(diào)整糧食進(jìn)出口貿(mào)易關(guān)系,保證糧食安全。

    第三,在糧食貿(mào)易網(wǎng)絡(luò)中,核心貿(mào)易國的出口量非常龐大,與它們存在貿(mào)易關(guān)系的國家數(shù)量也非常大,這些國家的糧食出口是否穩(wěn)定會(huì)直接對國家糧食貿(mào)易造成影響。因此,作為進(jìn)口國,為了防范進(jìn)口來源國的貿(mào)易惡意中斷,可以采取進(jìn)出口貿(mào)易多元化戰(zhàn)略,將糧食進(jìn)口分散化,并且,建立多種糧食作物的貿(mào)易關(guān)系。同時(shí),加強(qiáng)多國之間的貿(mào)易合作,以便保障本國糧食安全。另一方面,出口國中斷出口貿(mào)易,自給自足的方式不能從根本上解決問題,世界是一個(gè)經(jīng)濟(jì)共同體,在糧食生產(chǎn)基本能夠保障的情況下,應(yīng)該繼續(xù)保持與其它國家的貿(mào)易關(guān)系,保證整個(gè)世界糧食供應(yīng)穩(wěn)定,保障世界糧食安全。

    第四,糧食供應(yīng)安全問題是我國政府一直以來都很重視的問題,中國的糧食安全又和世界糧食安全緊密相關(guān)。中國在國際糧食網(wǎng)絡(luò)中居于重要地位;此外,中國糧食的對外依存度較低,在新冠疫情的影響下,國內(nèi)的糧食供應(yīng)保持穩(wěn)定,糧食價(jià)格也并未出現(xiàn)大幅度波動(dòng)。但作為糧食進(jìn)口大國之一,在呼吁人民節(jié)約糧食并增加國內(nèi)糧食產(chǎn)量的同時(shí),中國應(yīng)該積極推動(dòng)國際糧食貿(mào)易正常進(jìn)行,建立多元化的貿(mào)易關(guān)系,增加糧食貿(mào)易網(wǎng)絡(luò)的魯棒性,防止單一進(jìn)口國貿(mào)易中斷帶來供給危機(jī)。

    [參考文獻(xiàn)]

    [1]孫致陸,李先德.聯(lián)合國可持續(xù)發(fā)展目標(biāo)下的非洲糧食生產(chǎn)演變趨勢與發(fā)展?jié)摿Ψ治鯷J].中國農(nóng)業(yè)大學(xué)學(xué)報(bào),2020,25(2):160-170.

    [2]LU CHUNXIA,LIU AIMIN,XIAO YU,et al.Changes in Chinas grain production pattern and the effects of urbanization and dietary structure[J].Journal of resources and ecology,2020,11(4):358-365.

    [3]馬林靜,王雅鵬,吳娟.中國糧食生產(chǎn)技術(shù)效率的空間非均衡與收斂性分析[J].農(nóng)業(yè)技術(shù)經(jīng)濟(jì),2015(4):4-12.

    [4]曹寶明,劉婷,虞松波.中國糧食流通體制改革:目標(biāo)、路徑與重啟[J].農(nóng)業(yè)經(jīng)濟(jì)問題,2018(12):33-38.

    [5]孔貴寶,張趙晉.我國糧食流通現(xiàn)狀及高質(zhì)量發(fā)展路徑探析[J].商業(yè)經(jīng)濟(jì)研究,2020(5):177-180.

    [6]封志明,肖池偉,李鵬.中國-東盟自由貿(mào)易區(qū)糧食生產(chǎn)與貿(mào)易的時(shí)空格局演變[J].自然資源學(xué)報(bào),2017,32(6):937-950.

    [7]APARICIO G, PINILLA V. International trade in cereals and the collapse of the first wave of globalization, 1900-1938[J].Glob hist ,2019,14(1):44-67.

    [8]余志剛,王亞.供給側(cè)改革背景下國際糧食貿(mào)易形勢及對中國糧食安全的影響[J].世界農(nóng)業(yè),2017(8):4-9.

    [9]張晶,周海川.國際大米價(jià)格互動(dòng)性與中國糧食安全研究[J].中國人口·資源與環(huán)境,2014,24(10):163-169.

    [10]蘇珊珊.世界糧食貿(mào)易視角下的中國糧食安全問題研究[J].區(qū)域經(jīng)濟(jì)評論,2015(4):65-71.

    [11]WANG RUI,LIU CHUN-YU. International grain trade policy after the crisis: a limited return of protectionism and its roots [J]. Research of agricultural modernization,2015,6:954-959.

    [12]ALBERT R, JEONG H,BARABSI AL. Error and attack tolerance of complex networks [J].Nature,2000,406:378-382.

    [13]劉嬋娟,胡志華.基于復(fù)雜網(wǎng)絡(luò)的全球海運(yùn)網(wǎng)絡(luò)魯棒性研究[J].廣西大學(xué)學(xué)報(bào)(自然科學(xué)版),2016,41(5):1441-1448.

    [14]周墨竹,王介勇.基于復(fù)雜網(wǎng)絡(luò)的全球稻米貿(mào)易格局演化及其啟示[J].自然資源學(xué)報(bào),2020,35(5):1055-1067.

    [15]王祥,強(qiáng)文麗,牛叔文,等.全球農(nóng)產(chǎn)品貿(mào)易網(wǎng)絡(luò)及其演化分析[J].自然資源學(xué)報(bào),2018,33(6):940-953.

    [16]張振江,張玉召,王小榮.鐵路快捷貨運(yùn)網(wǎng)絡(luò)魯棒性分析[J].中國安全科學(xué)學(xué)報(bào),2020,30(3):150-156.

    Analysis on the Robustness of International Grain Trade Network

    Dong? Zhiliang,? Yang? Qiaoran

    (Hebei GEO University, Shijiazhuang 050031, China)

    Abstract: At present, the COVID-19 is still wreaking havoc on the human activities in the world. And the security of grain supply, the most important part of human life, is drawing increasing attention in many countries. In order to understand the basic characteristics of international grain trade and study whether the international grain trade is safe and stable, on the basis of complex network theory, this paper studied the network characteristics of international grain trade network and analyzed its robustness. The authors constructed the global grain trade network by using the grain trade data of 2019 released by the United Nations Statistics Department. They analyzed the characteristics of international grain trade by calculating the network characteristic parameters. Then, the important nodes in the trade network were evaluated and the robustness analysis index was selected. Finally, the robustness index under different attack modes was calculated to analyze the grain trade network under different attacks. The results show that: the nodes of grain trade network have good robustness when facing random attacks, but is vulnerable under targeted attacks. Through this study, we hope to provide a reference for the grain supply security for the countries participating in the international grain trade.

    Key words: grain trade; complex network; grain security; Robustness

    (責(zé)任編輯:李 萌)

    收稿日期:2020-12-02

    基金項(xiàng)目:國家社會(huì)科學(xué)基金項(xiàng)目《“互聯(lián)網(wǎng)+”對京津冀傳統(tǒng)產(chǎn)業(yè)價(jià)值鏈的破壞性重構(gòu)作用機(jī)理研究》(17BGL202)。

    作者簡介:董志良(1975—),男,河北元氏人,教授,研究方向?yàn)榧夹g(shù)經(jīng)濟(jì)及管理;楊巧然(1997—),女,河北秦皇島人,碩士研究生,研究方向?yàn)榧夹g(shù)經(jīng)濟(jì)及管理。

    猜你喜歡
    復(fù)雜網(wǎng)絡(luò)糧食安全魯棒性
    荒漠綠洲區(qū)潛在生態(tài)網(wǎng)絡(luò)增邊優(yōu)化魯棒性分析
    基于確定性指標(biāo)的弦支結(jié)構(gòu)魯棒性評價(jià)
    基于圖熵聚類的重疊社區(qū)發(fā)現(xiàn)算法
    廣義水資源利用效率綜合評價(jià)指數(shù)的
    城鎮(zhèn)化加速背景下河北省糧食儲(chǔ)備形式轉(zhuǎn)變分析
    強(qiáng)化農(nóng)田水利改革 確保糧食生產(chǎn)安全
    農(nóng)村科學(xué)儲(chǔ)糧與糧食安全研究
    基于復(fù)雜網(wǎng)絡(luò)理論的通用機(jī)場保障網(wǎng)絡(luò)研究
    城市群復(fù)合交通網(wǎng)絡(luò)復(fù)雜性實(shí)證研究
    科技視界(2016年20期)2016-09-29 11:19:34
    基于非支配解集的多模式裝備項(xiàng)目群調(diào)度魯棒性優(yōu)化
    亚洲 欧美一区二区三区| 亚洲电影在线观看av| 久久久久久免费高清国产稀缺| 国产精品香港三级国产av潘金莲| 国产精品久久久久久人妻精品电影| 午夜激情福利司机影院| 日日爽夜夜爽网站| 99久久久亚洲精品蜜臀av| 国内揄拍国产精品人妻在线 | 两人在一起打扑克的视频| 国产99久久九九免费精品| 成年版毛片免费区| 丰满的人妻完整版| 啪啪无遮挡十八禁网站| 亚洲avbb在线观看| 手机成人av网站| 最好的美女福利视频网| 丝袜美腿诱惑在线| 亚洲欧美精品综合久久99| 91九色精品人成在线观看| 天天躁夜夜躁狠狠躁躁| 日韩欧美三级三区| 香蕉国产在线看| 一本精品99久久精品77| 久久国产精品影院| 国产一卡二卡三卡精品| 在线观看午夜福利视频| 精品一区二区三区av网在线观看| 99久久久亚洲精品蜜臀av| www.999成人在线观看| 丝袜人妻中文字幕| 91老司机精品| 妹子高潮喷水视频| 亚洲中文av在线| 久久午夜亚洲精品久久| 窝窝影院91人妻| 亚洲国产欧美日韩在线播放| 久久久久久国产a免费观看| 好男人在线观看高清免费视频 | 露出奶头的视频| 国产精品美女特级片免费视频播放器 | 波多野结衣高清作品| www.www免费av| 亚洲国产精品久久男人天堂| 老司机在亚洲福利影院| 午夜免费观看网址| 午夜福利视频1000在线观看| 亚洲 欧美一区二区三区| 老司机靠b影院| 亚洲av成人av| 很黄的视频免费| 麻豆一二三区av精品| 99久久久亚洲精品蜜臀av| 啦啦啦 在线观看视频| 男人操女人黄网站| 热re99久久国产66热| 亚洲专区中文字幕在线| 丝袜美腿诱惑在线| 亚洲欧美精品综合久久99| 欧美精品啪啪一区二区三区| 亚洲精品久久成人aⅴ小说| 夜夜爽天天搞| 国产精品一区二区免费欧美| 男人舔奶头视频| 欧美性猛交黑人性爽| 国产麻豆成人av免费视频| 99国产精品99久久久久| 女人高潮潮喷娇喘18禁视频| 国产亚洲精品久久久久5区| 久久久久久亚洲精品国产蜜桃av| 亚洲午夜精品一区,二区,三区| 国产真人三级小视频在线观看| 在线免费观看的www视频| 男人舔女人下体高潮全视频| 欧美日本亚洲视频在线播放| 亚洲国产看品久久| 国产高清激情床上av| 天堂√8在线中文| 黄频高清免费视频| 亚洲色图 男人天堂 中文字幕| 国产v大片淫在线免费观看| 精品一区二区三区av网在线观看| 99国产精品一区二区三区| 黄色 视频免费看| avwww免费| 91成年电影在线观看| 日韩大码丰满熟妇| aaaaa片日本免费| 久久久久久久久免费视频了| 白带黄色成豆腐渣| av中文乱码字幕在线| 日韩欧美一区二区三区在线观看| 特大巨黑吊av在线直播 | 国产久久久一区二区三区| 国产成人av教育| 老司机福利观看| 99精品在免费线老司机午夜| 老熟妇乱子伦视频在线观看| 999久久久国产精品视频| 青草久久国产| 青草久久国产| 欧美另类亚洲清纯唯美| 久久国产亚洲av麻豆专区| 日本成人三级电影网站| 波多野结衣高清无吗| 此物有八面人人有两片| 免费在线观看完整版高清| 亚洲国产欧美一区二区综合| 夜夜躁狠狠躁天天躁| 每晚都被弄得嗷嗷叫到高潮| 精品免费久久久久久久清纯| 国产成人av激情在线播放| 99久久久亚洲精品蜜臀av| or卡值多少钱| 国产精品 欧美亚洲| 一边摸一边做爽爽视频免费| 亚洲av电影在线进入| 可以免费在线观看a视频的电影网站| 色老头精品视频在线观看| 久久婷婷人人爽人人干人人爱| 在线观看免费日韩欧美大片| 日日爽夜夜爽网站| 亚洲欧美日韩无卡精品| a级毛片在线看网站| 91大片在线观看| 久久久久久人人人人人| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品国产区一区二| 欧美激情高清一区二区三区| 日韩欧美在线二视频| 亚洲人成网站在线播放欧美日韩| 国产一区二区三区在线臀色熟女| 正在播放国产对白刺激| 国产男靠女视频免费网站| 亚洲欧美一区二区三区黑人| 级片在线观看| 国产男靠女视频免费网站| 欧美zozozo另类| 亚洲精品国产一区二区精华液| 亚洲熟女毛片儿| 天堂√8在线中文| 99久久精品国产亚洲精品| 少妇熟女aⅴ在线视频| 少妇被粗大的猛进出69影院| 很黄的视频免费| 欧美日韩瑟瑟在线播放| 亚洲精品国产精品久久久不卡| 国产精品亚洲av一区麻豆| 亚洲一区中文字幕在线| 黄色视频不卡| 精品久久久久久成人av| 久久婷婷成人综合色麻豆| 亚洲人成77777在线视频| 精品久久蜜臀av无| 12—13女人毛片做爰片一| 丝袜美腿诱惑在线| 中文字幕最新亚洲高清| 日韩精品青青久久久久久| 18禁美女被吸乳视频| 精品久久久久久久久久久久久 | www.www免费av| 久久 成人 亚洲| 国产又爽黄色视频| 午夜精品久久久久久毛片777| 亚洲国产日韩欧美精品在线观看 | 真人做人爱边吃奶动态| 久久狼人影院| 在线观看日韩欧美| 欧美 亚洲 国产 日韩一| 欧美午夜高清在线| 桃色一区二区三区在线观看| 一级a爱视频在线免费观看| 国产av不卡久久| 亚洲免费av在线视频| 白带黄色成豆腐渣| 日韩视频一区二区在线观看| 欧洲精品卡2卡3卡4卡5卡区| 最好的美女福利视频网| 国产欧美日韩一区二区精品| 色老头精品视频在线观看| 波多野结衣高清作品| 欧美丝袜亚洲另类 | 色老头精品视频在线观看| 国产成年人精品一区二区| x7x7x7水蜜桃| 在线观看66精品国产| 国产精品九九99| 亚洲第一欧美日韩一区二区三区| 亚洲专区中文字幕在线| 一本久久中文字幕| 久久久久国产精品人妻aⅴ院| 国产精品一区二区三区四区久久 | 免费在线观看影片大全网站| 老鸭窝网址在线观看| 一级毛片精品| 国内揄拍国产精品人妻在线 | 亚洲电影在线观看av| 欧美性猛交╳xxx乱大交人| 久久伊人香网站| 黄色 视频免费看| 51午夜福利影视在线观看| 日本一本二区三区精品| 国产伦一二天堂av在线观看| 中出人妻视频一区二区| 听说在线观看完整版免费高清| 欧美zozozo另类| 精品欧美一区二区三区在线| 好看av亚洲va欧美ⅴa在| 亚洲av成人不卡在线观看播放网| 琪琪午夜伦伦电影理论片6080| 欧美日韩福利视频一区二区| 国产精品电影一区二区三区| 国产高清视频在线播放一区| 19禁男女啪啪无遮挡网站| 后天国语完整版免费观看| 正在播放国产对白刺激| 欧美成人午夜精品| 91麻豆精品激情在线观看国产| 又大又爽又粗| 最新美女视频免费是黄的| 精品国产超薄肉色丝袜足j| 国产私拍福利视频在线观看| 18美女黄网站色大片免费观看| 一级毛片女人18水好多| 两个人免费观看高清视频| 国产极品粉嫩免费观看在线| 亚洲中文字幕日韩| 91av网站免费观看| 久久热在线av| 亚洲精品久久成人aⅴ小说| 99久久国产精品久久久| 精品国产超薄肉色丝袜足j| 日本a在线网址| 亚洲 欧美一区二区三区| 亚洲片人在线观看| 久久精品国产99精品国产亚洲性色| 999久久久精品免费观看国产| 国产高清有码在线观看视频 | 久久久水蜜桃国产精品网| 国产区一区二久久| 黄色 视频免费看| 国产精品久久久久久人妻精品电影| 波多野结衣巨乳人妻| 丰满的人妻完整版| 一卡2卡三卡四卡精品乱码亚洲| 18禁观看日本| 国产av又大| 老司机在亚洲福利影院| 亚洲精品中文字幕一二三四区| 中文字幕av电影在线播放| 久久婷婷人人爽人人干人人爱| xxxwww97欧美| 午夜福利高清视频| 午夜精品在线福利| 女性被躁到高潮视频| 欧美激情 高清一区二区三区| 免费看日本二区| 日本免费一区二区三区高清不卡| 色精品久久人妻99蜜桃| 色婷婷久久久亚洲欧美| 制服丝袜大香蕉在线| 欧美最黄视频在线播放免费| 宅男免费午夜| 在线播放国产精品三级| 天天躁夜夜躁狠狠躁躁| 日本a在线网址| 国产国语露脸激情在线看| 国产一区二区激情短视频| 夜夜爽天天搞| 国产亚洲欧美精品永久| 激情在线观看视频在线高清| 国产精品1区2区在线观看.| 99精品欧美一区二区三区四区| 亚洲成人久久爱视频| 免费一级毛片在线播放高清视频| 无遮挡黄片免费观看| 黄频高清免费视频| 天堂影院成人在线观看| 亚洲 欧美 日韩 在线 免费| 久久香蕉国产精品| 婷婷丁香在线五月| 亚洲熟女毛片儿| 亚洲精品美女久久久久99蜜臀| 欧美日韩一级在线毛片| av电影中文网址| 一a级毛片在线观看| 国产精华一区二区三区| 亚洲一区二区三区不卡视频| 国产精品久久久久久精品电影 | 国产亚洲精品一区二区www| 国产av一区二区精品久久| 久久精品成人免费网站| 日韩大尺度精品在线看网址| 老熟妇仑乱视频hdxx| 听说在线观看完整版免费高清| 一卡2卡三卡四卡精品乱码亚洲| 国产精品98久久久久久宅男小说| 亚洲精品色激情综合| 别揉我奶头~嗯~啊~动态视频| 男男h啪啪无遮挡| 美女午夜性视频免费| 一本一本综合久久| 好看av亚洲va欧美ⅴa在| 黄片小视频在线播放| 免费在线观看成人毛片| 国产黄片美女视频| 日本精品一区二区三区蜜桃| 很黄的视频免费| 国产精品99久久99久久久不卡| 亚洲精品粉嫩美女一区| 夜夜夜夜夜久久久久| 国产成人精品久久二区二区免费| 国产精品久久久人人做人人爽| 热99re8久久精品国产| 日日夜夜操网爽| 日韩精品青青久久久久久| 亚洲精品美女久久av网站| 真人一进一出gif抽搐免费| 亚洲av片天天在线观看| 妹子高潮喷水视频| 国产av不卡久久| 精品欧美一区二区三区在线| 亚洲人成网站在线播放欧美日韩| 午夜免费观看网址| 久久久久九九精品影院| 人妻丰满熟妇av一区二区三区| 18禁观看日本| 久久中文字幕人妻熟女| 精品久久久久久成人av| 亚洲精品国产区一区二| 亚洲精品中文字幕一二三四区| 日本一本二区三区精品| 韩国精品一区二区三区| 黄频高清免费视频| 91大片在线观看| 婷婷精品国产亚洲av在线| 嫁个100分男人电影在线观看| 国产精品野战在线观看| 搡老岳熟女国产| 在线观看日韩欧美| 777久久人妻少妇嫩草av网站| videosex国产| 亚洲欧洲精品一区二区精品久久久| 一二三四在线观看免费中文在| 亚洲,欧美精品.| 在线观看日韩欧美| 久久久久国产精品人妻aⅴ院| 男女床上黄色一级片免费看| 免费在线观看成人毛片| 久久久久免费精品人妻一区二区 | 日韩欧美免费精品| 亚洲国产中文字幕在线视频| 精品少妇一区二区三区视频日本电影| 色综合站精品国产| 伦理电影免费视频| 2021天堂中文幕一二区在线观 | 女性被躁到高潮视频| 亚洲成人国产一区在线观看| 中文字幕精品亚洲无线码一区 | 国产在线观看jvid| 后天国语完整版免费观看| 美国免费a级毛片| 欧美激情 高清一区二区三区| 在线看三级毛片| 亚洲人成网站在线播放欧美日韩| 琪琪午夜伦伦电影理论片6080| 国产精品久久电影中文字幕| 成人国产一区最新在线观看| 亚洲五月天丁香| 午夜日韩欧美国产| 一区二区日韩欧美中文字幕| 老司机午夜十八禁免费视频| 黄网站色视频无遮挡免费观看| 久久人妻福利社区极品人妻图片| 夜夜躁狠狠躁天天躁| 欧美av亚洲av综合av国产av| 99久久无色码亚洲精品果冻| 中文字幕久久专区| 国产不卡一卡二| bbb黄色大片| 国产精品1区2区在线观看.| 91av网站免费观看| 免费观看精品视频网站| 国产三级在线视频| 黄片大片在线免费观看| 欧美性长视频在线观看| 国产成人av教育| 久久久久国产精品人妻aⅴ院| 亚洲精品国产精品久久久不卡| 黑人欧美特级aaaaaa片| 两个人免费观看高清视频| 最新在线观看一区二区三区| 午夜两性在线视频| 亚洲性夜色夜夜综合| 伦理电影免费视频| 成人三级黄色视频| 非洲黑人性xxxx精品又粗又长| 日韩中文字幕欧美一区二区| av视频在线观看入口| 亚洲专区字幕在线| 国产免费av片在线观看野外av| 亚洲人成伊人成综合网2020| 精品少妇一区二区三区视频日本电影| 国产精品日韩av在线免费观看| 少妇的丰满在线观看| 国产成人欧美| 又大又爽又粗| 又黄又爽又免费观看的视频| 国产精品爽爽va在线观看网站 | 一区二区三区高清视频在线| 男人的好看免费观看在线视频 | 一级黄色大片毛片| 免费女性裸体啪啪无遮挡网站| 国产精品野战在线观看| 宅男免费午夜| 亚洲第一欧美日韩一区二区三区| 在线观看www视频免费| av在线天堂中文字幕| 一本精品99久久精品77| 无遮挡黄片免费观看| 欧美色视频一区免费| 亚洲中文字幕日韩| 久久亚洲真实| 99久久99久久久精品蜜桃| 99精品欧美一区二区三区四区| 成人三级做爰电影| 午夜福利在线在线| 国产国语露脸激情在线看| 国产精品亚洲一级av第二区| 99久久无色码亚洲精品果冻| 国产亚洲精品av在线| 19禁男女啪啪无遮挡网站| 欧美国产日韩亚洲一区| 91成人精品电影| 男人操女人黄网站| 一级a爱视频在线免费观看| 一个人免费在线观看的高清视频| 麻豆久久精品国产亚洲av| 亚洲熟妇熟女久久| 高清在线国产一区| 巨乳人妻的诱惑在线观看| 1024手机看黄色片| 在线观看www视频免费| 欧美黑人巨大hd| 精品国产国语对白av| 免费高清视频大片| 国产在线精品亚洲第一网站| 久热这里只有精品99| 91麻豆av在线| 人人妻,人人澡人人爽秒播| 无人区码免费观看不卡| 精品一区二区三区视频在线观看免费| 69av精品久久久久久| 可以在线观看毛片的网站| 亚洲第一青青草原| 午夜久久久在线观看| 韩国精品一区二区三区| 别揉我奶头~嗯~啊~动态视频| 亚洲电影在线观看av| 一本久久中文字幕| 日本成人三级电影网站| 日日干狠狠操夜夜爽| 看片在线看免费视频| 十八禁网站免费在线| 中文在线观看免费www的网站 | 日韩成人在线观看一区二区三区| 欧美色欧美亚洲另类二区| 午夜免费成人在线视频| 国产真实乱freesex| 久久婷婷成人综合色麻豆| 黄片大片在线免费观看| 国产乱人伦免费视频| 露出奶头的视频| 成人一区二区视频在线观看| 黄色视频不卡| 国产成人精品久久二区二区免费| 久久久久久久午夜电影| 一级毛片女人18水好多| 亚洲精品美女久久久久99蜜臀| 久久久久精品国产欧美久久久| www日本黄色视频网| 色综合欧美亚洲国产小说| 男女做爰动态图高潮gif福利片| 精品久久久久久久久久久久久 | 国产极品粉嫩免费观看在线| 久久久久久免费高清国产稀缺| 日韩高清综合在线| 日韩精品中文字幕看吧| 极品教师在线免费播放| 看黄色毛片网站| 久久久水蜜桃国产精品网| 一级毛片女人18水好多| 中文字幕人妻丝袜一区二区| xxx96com| 日韩欧美三级三区| 一本精品99久久精品77| 久久精品aⅴ一区二区三区四区| √禁漫天堂资源中文www| av欧美777| 精品久久久久久成人av| av欧美777| 国产v大片淫在线免费观看| 国产又爽黄色视频| 国产区一区二久久| 91字幕亚洲| 国产欧美日韩一区二区精品| 欧美一区二区精品小视频在线| 成人手机av| 悠悠久久av| 操出白浆在线播放| 免费看a级黄色片| www国产在线视频色| 国内揄拍国产精品人妻在线 | 婷婷亚洲欧美| 免费女性裸体啪啪无遮挡网站| 亚洲国产欧美网| 亚洲人成77777在线视频| 欧美乱码精品一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 久久久久免费精品人妻一区二区 | 非洲黑人性xxxx精品又粗又长| 欧美成人一区二区免费高清观看 | 1024手机看黄色片| 亚洲国产欧美日韩在线播放| av福利片在线| 日韩欧美国产一区二区入口| 午夜老司机福利片| 999久久久国产精品视频| 亚洲成人久久爱视频| 日本五十路高清| 精品欧美一区二区三区在线| 久久99热这里只有精品18| 欧美国产日韩亚洲一区| 午夜免费观看网址| 午夜久久久久精精品| 国产成年人精品一区二区| 怎么达到女性高潮| 亚洲av电影在线进入| 超碰成人久久| 90打野战视频偷拍视频| 久久精品国产综合久久久| 99国产综合亚洲精品| 成人永久免费在线观看视频| 精品国产超薄肉色丝袜足j| 一边摸一边抽搐一进一小说| 国产亚洲精品久久久久久毛片| 禁无遮挡网站| 黄色成人免费大全| 大香蕉久久成人网| 亚洲一区二区三区色噜噜| 国产av又大| 日韩免费av在线播放| 成年免费大片在线观看| 亚洲精品久久成人aⅴ小说| 亚洲欧洲精品一区二区精品久久久| 日韩欧美一区视频在线观看| 日韩欧美国产在线观看| www.999成人在线观看| 日日摸夜夜添夜夜添小说| 在线天堂中文资源库| 国产成人精品久久二区二区免费| 99国产综合亚洲精品| 不卡av一区二区三区| 亚洲国产精品成人综合色| 一级黄色大片毛片| 真人一进一出gif抽搐免费| 国内少妇人妻偷人精品xxx网站 | 亚洲av成人不卡在线观看播放网| 亚洲一区高清亚洲精品| 国产又色又爽无遮挡免费看| 精品久久久久久久久久久久久 | 久99久视频精品免费| 久久久久亚洲av毛片大全| www国产在线视频色| 国产高清激情床上av| 国产在线精品亚洲第一网站| 久久久久久国产a免费观看| 久久久久免费精品人妻一区二区 | 欧美国产日韩亚洲一区| 天堂动漫精品| 曰老女人黄片| 很黄的视频免费| 亚洲男人天堂网一区| 亚洲一区中文字幕在线| 亚洲人成伊人成综合网2020| 亚洲男人天堂网一区| 色老头精品视频在线观看| 久久久精品国产亚洲av高清涩受| 免费在线观看完整版高清| 亚洲国产精品合色在线| 午夜激情福利司机影院| 久久久久久久久久黄片| 757午夜福利合集在线观看| 国产成人一区二区三区免费视频网站| 黄片小视频在线播放| 又黄又爽又免费观看的视频| 三级毛片av免费| 88av欧美| 国产黄片美女视频| 十分钟在线观看高清视频www| 亚洲国产欧洲综合997久久, | 性色av乱码一区二区三区2| 波多野结衣巨乳人妻| 国产亚洲精品久久久久久毛片| a在线观看视频网站| 亚洲一区二区三区色噜噜| 国产欧美日韩一区二区精品| 91字幕亚洲| 国产麻豆成人av免费视频| 国产在线观看jvid| 国产精品一区二区免费欧美| 欧美激情极品国产一区二区三区| 国产成人精品久久二区二区免费| 狠狠狠狠99中文字幕| 欧美成人午夜精品|