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

    基于網(wǎng)絡(luò)視角的銀行業(yè)系統(tǒng)性風(fēng)險(xiǎn)度量方法

    2016-12-28 02:04:02譚照林王宗堯
    中國管理科學(xué) 2016年5期
    關(guān)鍵詞:傳染系統(tǒng)性沖擊

    隋 聰,譚照林,王宗堯

    (1.東北財(cái)經(jīng)大學(xué)金融學(xué)院,遼寧 大連 116025;2.東北財(cái)經(jīng)大學(xué)商品市場與行為決策研究中心,遼寧 大連 116025;3.東北財(cái)經(jīng)大學(xué)薩里國際學(xué)院,遼寧 大連 116025)

    ?

    基于網(wǎng)絡(luò)視角的銀行業(yè)系統(tǒng)性風(fēng)險(xiǎn)度量方法

    隋 聰1,2,譚照林1,王宗堯3

    (1.東北財(cái)經(jīng)大學(xué)金融學(xué)院,遼寧 大連 116025;2.東北財(cái)經(jīng)大學(xué)商品市場與行為決策研究中心,遼寧 大連 116025;3.東北財(cái)經(jīng)大學(xué)薩里國際學(xué)院,遼寧 大連 116025)

    網(wǎng)絡(luò)模型已經(jīng)成為研究銀行系統(tǒng)性風(fēng)險(xiǎn)的重要方法。然而現(xiàn)有研究忽視了銀行系統(tǒng)性風(fēng)險(xiǎn)的小概率特點(diǎn),同時(shí)也缺少度量銀行系統(tǒng)性風(fēng)險(xiǎn)的統(tǒng)一標(biāo)準(zhǔn)。為此,本文提出了基于網(wǎng)絡(luò)模型的銀行系統(tǒng)性風(fēng)險(xiǎn)度量方法:銀行系統(tǒng)性風(fēng)險(xiǎn)VaR和銀行系統(tǒng)性風(fēng)險(xiǎn)ES。首先,本文采用蒙特卡洛模擬方法,模擬銀行外部沖擊造成銀行間網(wǎng)絡(luò)損失的大樣本。在銀行間網(wǎng)絡(luò)損失大樣本中,估計(jì)銀行系統(tǒng)性風(fēng)險(xiǎn)VaR和銀行系統(tǒng)性風(fēng)險(xiǎn)ES。這兩個(gè)測度能夠捕捉到銀行間網(wǎng)絡(luò)損失的尾部特征,解決了對(duì)比隨機(jī)沖擊結(jié)果無法反映銀行系統(tǒng)性風(fēng)險(xiǎn)的問題。其次,在模擬實(shí)驗(yàn)中,本文利用真實(shí)銀行間網(wǎng)絡(luò)結(jié)構(gòu)參數(shù),對(duì)模擬的三種銀行間網(wǎng)絡(luò)進(jìn)行校準(zhǔn),保證了研究結(jié)論真實(shí)性和可靠性。最后,在模擬實(shí)驗(yàn)中發(fā)現(xiàn):(1)外部沖擊會(huì)引發(fā)違約傳染的連鎖反應(yīng),并導(dǎo)致銀行間網(wǎng)絡(luò)損失分布從近似正態(tài)分布轉(zhuǎn)變成尖峰厚尾分布,最后變成雙峰分布。(2)網(wǎng)絡(luò)集中度越高發(fā)生違約傳染連鎖反應(yīng)的概率越小,但是傳染的破壞力會(huì)更大。(3)銀行間網(wǎng)絡(luò)的潛在傳染作用會(huì)極大的放大銀行系統(tǒng)的風(fēng)險(xiǎn),而且違約傳染效應(yīng)是呈指數(shù)增長的。

    系統(tǒng)性風(fēng)險(xiǎn);銀行間網(wǎng)絡(luò);蒙特卡洛模擬;VaR

    1 引言

    銀行間同業(yè)市場為銀行間調(diào)劑流動(dòng)性提供便利的同時(shí),也為危機(jī)的蔓延提供了傳染渠道。通過銀行間同業(yè)市場形成的銀行間的借貸網(wǎng)絡(luò),構(gòu)成了復(fù)雜的債權(quán)債務(wù)網(wǎng)絡(luò)——銀行間網(wǎng)絡(luò)。銀行間網(wǎng)絡(luò)在系統(tǒng)性風(fēng)險(xiǎn)的蔓延過程中扮演著重要的角色,越來越受到金融學(xué)家的重視。Iori等[1]研究了銀行間市場作為安全網(wǎng)的作用,研究發(fā)現(xiàn)銀行間市場的結(jié)構(gòu)會(huì)影響其發(fā)揮相互保護(hù)的作用。如果銀行是同質(zhì)的,銀行間市場有利于穩(wěn)定系統(tǒng);如果銀行是異質(zhì)的,銀行間市場穩(wěn)定系統(tǒng)的作用就會(huì)下降。近些年,涌現(xiàn)出大量利用銀行間網(wǎng)絡(luò)模型研究銀行系統(tǒng)性風(fēng)險(xiǎn)的文獻(xiàn)。

    Allen和Gale[2]最早提出了網(wǎng)絡(luò)視角下的銀行間風(fēng)險(xiǎn)傳染的問題。他們的研究認(rèn)為完全結(jié)構(gòu)網(wǎng)絡(luò)比不完全結(jié)構(gòu)網(wǎng)絡(luò)更穩(wěn)定。Nier等[3]研究了銀行間網(wǎng)絡(luò)的連接度對(duì)風(fēng)險(xiǎn)傳染的影響,并發(fā)現(xiàn)網(wǎng)絡(luò)連接度存在閾值。當(dāng)網(wǎng)絡(luò)連接度超過閾值時(shí),違約傳染的概率下降;網(wǎng)絡(luò)連接度低于閾值時(shí),風(fēng)險(xiǎn)傳染發(fā)生的概率隨平均度上升而上升。隨后,May和Arinaminpathy[4]利用平均場近似技術(shù)對(duì)Nier等[3]的研究結(jié)果提供了一個(gè)完整細(xì)致的解釋。Gai和Kapadia[5]利用矩母函數(shù)技術(shù),對(duì)發(fā)生違約傳染的平均度相變位置作了深入的理論分析。隋聰?shù)萚6]研究發(fā)現(xiàn)集中程度越高的網(wǎng)絡(luò)由于傳染而倒閉的銀行數(shù)量就越多,風(fēng)險(xiǎn)傳染的差異并不明顯。Wells[7]對(duì)英國的銀行間市場的傳染風(fēng)險(xiǎn)做了分析。研究發(fā)現(xiàn),一家銀行違約引起其他銀行違約的概率很小,但是會(huì)降低銀行系統(tǒng)的資本充足率。損失的嚴(yán)重程度與銀行間拆借規(guī)模和違約損失率密切相關(guān)。Lublóy[8]在研究匈牙利銀行業(yè)時(shí)發(fā)現(xiàn),銀行間的網(wǎng)絡(luò)結(jié)構(gòu)對(duì)違約傳染具有重要影響。匈牙利的銀行間債務(wù)比例低于其他國家,而且銀行間網(wǎng)絡(luò)的集中程度適中。所以,發(fā)生連鎖傳染的概率低,連鎖傳染的破壞力也低。Mistrilli[9]的研究結(jié)果進(jìn)一步證實(shí)了銀行間網(wǎng)絡(luò)結(jié)構(gòu)會(huì)影響違約傳染。并且發(fā)現(xiàn)當(dāng)違約損失率較高(大于0.8)時(shí),在最大熵方法估計(jì)的完全結(jié)構(gòu)網(wǎng)絡(luò)中,違約傳染的嚴(yán)重程度會(huì)更高。這一結(jié)論不同于Allen和Gale[2]的完全結(jié)構(gòu)網(wǎng)絡(luò)更穩(wěn)定的觀點(diǎn)。黃聰、賈彥東[10]研究了中國的銀行間網(wǎng)絡(luò),研究發(fā)現(xiàn)初始沖擊銀行的規(guī)模決定了傳染的影響力。只有四大行受到?jīng)_擊時(shí),才會(huì)對(duì)銀行間網(wǎng)絡(luò)產(chǎn)生極端影響。馬君潞等[11]同樣發(fā)現(xiàn)初始受到?jīng)_擊的銀行能夠決定傳染的力度和范圍。陳庭強(qiáng)和何建敏[12]構(gòu)建了基于復(fù)雜網(wǎng)絡(luò)的信用風(fēng)險(xiǎn)傳染模型,通過仿真實(shí)驗(yàn)發(fā)現(xiàn)網(wǎng)絡(luò)結(jié)構(gòu)異質(zhì)性越大,風(fēng)險(xiǎn)傳染的概率越小,傳染的影響范圍越小。此外,歐陽紅兵和劉曉東[13]利用最小生成樹和平面極大過濾圖方法,對(duì)銀行間網(wǎng)絡(luò)節(jié)點(diǎn)的系統(tǒng)重要性進(jìn)行了研究。

    上面這些基于網(wǎng)絡(luò)模型的研究都發(fā)現(xiàn),銀行間網(wǎng)絡(luò)結(jié)構(gòu)是影響風(fēng)險(xiǎn)傳染和銀行系統(tǒng)性風(fēng)險(xiǎn)的關(guān)鍵因素。但是,在網(wǎng)絡(luò)結(jié)構(gòu)是如何影響銀行系統(tǒng)性風(fēng)險(xiǎn)的問題上,研究的結(jié)論還存在差異。這主要有兩個(gè)方面的原因。

    第一,缺少系統(tǒng)性風(fēng)險(xiǎn)度量方法?;诰W(wǎng)絡(luò)模型的研究普遍采用違約銀行數(shù)量或違約銀行比率作為衡量系統(tǒng)性風(fēng)險(xiǎn)的指標(biāo)。在Nier等[3]的研究中,每次沖擊一家銀行,假設(shè)該銀行的外部資產(chǎn)全部損失。依次,對(duì)所有銀行單獨(dú)沖擊,并比較違約銀行數(shù)量的平均值以及觀測結(jié)果的95%范圍。針對(duì)每一個(gè)參數(shù),每組實(shí)驗(yàn)進(jìn)行100次,用違約銀行數(shù)量衡量系統(tǒng)性風(fēng)險(xiǎn)狀況。May和Arinaminpathy[4]在Nier等[3]的研究基礎(chǔ)上,增加考慮了違約銀行的頻率分布。Gai和Kapadia[5]采用了隨機(jī)的方式,每次選出一家銀行進(jìn)行沖擊,同樣假設(shè)該銀行的外部資產(chǎn)全部損失,并利用違約銀行比例衡量傳染范圍。同時(shí),Gai和Kapadia[5]將傳染超過5家定義為出現(xiàn)傳染,即發(fā)生系統(tǒng)性風(fēng)險(xiǎn),并利用出現(xiàn)傳染的頻率衡量系統(tǒng)性風(fēng)險(xiǎn)。隋聰?shù)萚6]同樣采用了隨機(jī)的沖擊方式,但與Gai和Kapadia[5]不同。首先,隨機(jī)選擇一家銀行,假設(shè)該銀行外部資產(chǎn)損失。然后,再隨機(jī)增加一家銀行進(jìn)行沖擊。同理,將沖擊不斷進(jìn)行下去。最后,利用違約銀行數(shù)量衡量系統(tǒng)性風(fēng)險(xiǎn)狀況。Lublóy[8]模擬了匈牙利的39家銀行50天的情景。每次沖擊一家銀行,假設(shè)該銀行損失其全部一級(jí)資本。利用第一輪和第二輪倒閉銀行數(shù)量度量系統(tǒng)性風(fēng)險(xiǎn)。黃聰和賈彥東[10],馬君潞等[11]模擬特定銀行發(fā)出流動(dòng)性沖擊,利用違約銀行比率、傳染輪次度量系統(tǒng)性風(fēng)險(xiǎn)。綜上,在銀行間網(wǎng)絡(luò)模型的研究中沒有統(tǒng)一的系統(tǒng)性風(fēng)險(xiǎn)測度方法。不同的研究對(duì)比的標(biāo)準(zhǔn)不同。而且,無論是隨機(jī)方式還是特定方式的沖擊都遠(yuǎn)不能涵蓋現(xiàn)實(shí)中所有可能出現(xiàn)的情景。事實(shí)上,銀行系統(tǒng)性風(fēng)險(xiǎn)往往是在極端情況下發(fā)生的,是一個(gè)典型的小概率事件。

    第二,網(wǎng)絡(luò)結(jié)構(gòu)差異大。Allen和Gale[2]利用4家銀行構(gòu)建的簡單網(wǎng)絡(luò)模型,比較了完全結(jié)構(gòu)網(wǎng)絡(luò)和不完全結(jié)構(gòu)網(wǎng)絡(luò)對(duì)風(fēng)險(xiǎn)傳染的影響。但是只有4家銀行的網(wǎng)絡(luò)太過簡單很難反映真實(shí)的銀行間網(wǎng)絡(luò)特征。Nier等[3],May和Arinaminpathy[4],Gai和Kapadia[5]以隨機(jī)網(wǎng)絡(luò)為基礎(chǔ),系統(tǒng)分析了違約傳染與影響因素的關(guān)系以及傳染過程,比較了不同連接度對(duì)系統(tǒng)性風(fēng)險(xiǎn)的影響。隋聰?shù)萚6]用無標(biāo)度網(wǎng)絡(luò)構(gòu)造銀行間網(wǎng)絡(luò),并研究了不同標(biāo)度參數(shù)的無標(biāo)度網(wǎng)絡(luò)對(duì)風(fēng)險(xiǎn)傳染以及系統(tǒng)性風(fēng)險(xiǎn)的作用。Wells[7]網(wǎng)絡(luò)模型的參數(shù)數(shù)據(jù)是基于英國銀行間市場的真實(shí)數(shù)據(jù)。Lublóy[8]采用了匈牙利的銀行間業(yè)務(wù)數(shù)據(jù)。同時(shí)發(fā)現(xiàn)匈牙利的銀行系統(tǒng)呈現(xiàn)一種多重貨幣中心結(jié)構(gòu),由10—15家大銀行扮演貨幣中心的角色。Mistrilli[9]利用意大利銀行間的雙邊債務(wù)數(shù)據(jù),對(duì)比了真實(shí)債務(wù)矩陣和利用最大熵方法估計(jì)的債務(wù)矩陣兩種銀行間網(wǎng)絡(luò)的違約傳染過程。黃聰和賈彥東[10]利用中國人民銀行支付系統(tǒng)的銀行間支付數(shù)據(jù),構(gòu)建銀行間網(wǎng)絡(luò)。馬君潞等[11]采用的是中國商業(yè)銀行的資產(chǎn)負(fù)債表數(shù)據(jù),并利用最大熵方法構(gòu)建銀行間網(wǎng)絡(luò)。需要指出,最大熵方法構(gòu)建的網(wǎng)絡(luò)是完全結(jié)構(gòu)網(wǎng)絡(luò),并不是真實(shí)的銀行間網(wǎng)絡(luò)??梢?,不同研究所采用的網(wǎng)絡(luò)結(jié)構(gòu)差異很大。

    針對(duì)這些問題,本文提出了銀行系統(tǒng)性風(fēng)險(xiǎn)VaR(Value at Risk)和銀行系統(tǒng)性風(fēng)險(xiǎn)ES(Expected Shortfall)兩種測度方法。該方法利用蒙特卡洛模擬確定銀行間網(wǎng)絡(luò)的損失分布,利用VaR和ES度量銀行系統(tǒng)性風(fēng)險(xiǎn)。通過蒙特卡洛模擬方法模擬外部沖擊可能出現(xiàn)的所有情景,利用銀行系統(tǒng)性風(fēng)險(xiǎn)VaR和銀行系統(tǒng)性風(fēng)險(xiǎn)ES度量極端情況下銀行間網(wǎng)絡(luò)的損失,即系統(tǒng)性風(fēng)險(xiǎn)。此外,Boss等[14]、Soram?ki 等[15]、Becher等[16]、Krause和Giansante[17]等文獻(xiàn)對(duì)不同國家的銀行間網(wǎng)絡(luò)結(jié)構(gòu)進(jìn)行了卓有成效的研究。因此,為了解決隨機(jī)網(wǎng)絡(luò)與規(guī)則網(wǎng)絡(luò)與真實(shí)銀行間網(wǎng)絡(luò)差異很大的問題,本文采用銀行間網(wǎng)絡(luò)結(jié)構(gòu)的相關(guān)研究結(jié)論對(duì)網(wǎng)絡(luò)進(jìn)行校準(zhǔn)。從而保證了研究結(jié)論符合真實(shí)銀行間網(wǎng)絡(luò)結(jié)構(gòu)的特點(diǎn)。

    2 基于蒙特卡洛模擬的系統(tǒng)性風(fēng)險(xiǎn)度量

    2.1 基礎(chǔ)模型

    2.1.1 債務(wù)矩陣

    為了研究銀行間的風(fēng)險(xiǎn)傳染,本文假設(shè)每家銀行的資產(chǎn)包括同業(yè)資產(chǎn)(Interbank Assets,IA)和外部資產(chǎn)(External Assets,EA),負(fù)債包括同業(yè)負(fù)債(Interbank Liability,IL)和存款(Deposits,D),以及股東權(quán)益(Equity,E)。

    銀行間的債務(wù)聯(lián)系將銀行連接成一個(gè)相互交織的銀行間網(wǎng)絡(luò)。利用債務(wù)矩陣可以概括銀行間網(wǎng)絡(luò)的基本要素。債務(wù)矩陣L表示為:

    (1)

    由于銀行間交易的私密性特點(diǎn),具體的銀行間貸款lij數(shù)據(jù)很難獲得。所以,現(xiàn)有研究確定銀行間債務(wù)矩陣有兩種方法。一種方法是全完采用模擬。另一種方法是根據(jù)真實(shí)銀行的總量數(shù)據(jù)(同業(yè)資產(chǎn)和同業(yè)資產(chǎn)),利用最大熵方法估計(jì)債務(wù)矩陣。最大熵方法最早由Sheldon和Maurer[18]引入銀行系統(tǒng)性風(fēng)險(xiǎn)研究領(lǐng)域,并受到學(xué)者們的廣泛認(rèn)可和使用。

    2.1.2 違約動(dòng)態(tài)

    (1)沖擊和傳染

    外部沖擊是指銀行系統(tǒng)外部事件對(duì)銀行系統(tǒng)造成的沖擊影響,導(dǎo)致某些銀行違約。本文利用蒙特卡洛模擬方法,模擬每家銀行外部資產(chǎn)受到不同程度的損失,來反映外部沖擊。在外部沖擊下,銀行違約可以表示為:

    (2)

    (3)

    (2)違約清算

    由于銀行間債務(wù)是一個(gè)復(fù)雜的網(wǎng)絡(luò),當(dāng)網(wǎng)絡(luò)中同時(shí)存在基礎(chǔ)違約和傳染違約銀行時(shí),確定銀行間債務(wù)的清算支付就變得十分困難。Eisenberg & Noe(2001)[19]提出了銀行間債務(wù)網(wǎng)絡(luò)的清算支付方法。銀行違約清算時(shí)按債務(wù)比例償還,因此債務(wù)矩陣中的每筆同業(yè)負(fù)債除以其總的同業(yè)負(fù)債,可以得到一個(gè)償還比例矩陣Π,其中的元素為:

    (4)

    銀行能夠收回的同業(yè)資產(chǎn)IA*取決于其他銀行能夠償還多少同業(yè)負(fù)債IL*,即:

    (5)

    因此計(jì)算能夠收回的同業(yè)資產(chǎn)IA*就轉(zhuǎn)換成計(jì)算能夠償還多少同業(yè)負(fù)債IL*。銀行i的償還支付表示為:

    (6)

    2.2 蒙特卡洛模擬

    銀行業(yè)的系統(tǒng)性風(fēng)險(xiǎn)是小概率事件。現(xiàn)有研究側(cè)重正常狀態(tài)下不同銀行間網(wǎng)絡(luò)結(jié)構(gòu)的比較分析,很難有效捕捉這一小概率事件,不足以解釋極端情況下的系統(tǒng)性風(fēng)險(xiǎn)的嚴(yán)重程度。研究這一小概率事件,需要大量的樣本,而現(xiàn)實(shí)世界不可能提供這么多的樣本。因此,本文利用蒙特卡洛模擬方法,模擬銀行間網(wǎng)絡(luò)遭受沖擊的多種情景,從而可以幫助我們研究系統(tǒng)性風(fēng)險(xiǎn)這樣極端事件。

    蒙特卡洛模擬基本思路:首先,選擇外部資產(chǎn)損失變動(dòng)的隨機(jī)過程和分布;其次,針對(duì)外部資產(chǎn)的每次變動(dòng),統(tǒng)計(jì)銀行間網(wǎng)絡(luò)的損失;最后,重復(fù)前述操作,得到足夠多的樣本數(shù)量。

    在蒙特卡洛模擬中,確定外部資產(chǎn)損失變動(dòng)的分布形式是基礎(chǔ)。由于本文旨在研究銀行間網(wǎng)絡(luò)的風(fēng)險(xiǎn)傳染,銀行外部資產(chǎn)的損失情況不是本文的研究重點(diǎn)。因此,本文將銀行外部資產(chǎn)的損失變動(dòng)看成是外生的,這并不會(huì)影響本文的研究結(jié)論。我們借鑒Iori等[1]的研究思路,假設(shè)銀行外部資產(chǎn)損失服從半正態(tài)分布。銀行外部資產(chǎn)的變動(dòng)表示為:

    EAi=EAi-|εi|EAi,εi~N(0,τi2)

    (7)

    其中,εi表示銀行i的外部資產(chǎn)的損失比例,εi服從標(biāo)準(zhǔn)差為τi的正態(tài)分布。稱τi為銀行資產(chǎn)的波動(dòng)率,它反映了銀行系統(tǒng)外部沖擊的大小。

    為簡化而不失一般性,假設(shè)每家銀行的εi是獨(dú)立且同分布的。因此式(7)等價(jià)于:

    EAi=EAi-|εi|EAi,εi~N(0,τ2)

    (8)

    銀行間網(wǎng)絡(luò)風(fēng)險(xiǎn)傳染的蒙特卡洛模擬基本步驟:

    1)確定外部沖擊大小(即標(biāo)準(zhǔn)差τ),根據(jù)式(8)對(duì)每家銀行的外部資產(chǎn)EAi進(jìn)行隨機(jī)抽樣。

    2)根據(jù)式(6)計(jì)算每家銀行的償還支付,然后分別依據(jù)式(2)和(3)判斷基礎(chǔ)違約和傳染違約。進(jìn)而統(tǒng)計(jì)總違約銀行數(shù)量、基礎(chǔ)違約銀行的數(shù)量、傳染違約銀行的數(shù)量,用以衡量銀行系統(tǒng)的損失。

    3)重復(fù)上述步驟M次,獲得M個(gè)銀行間網(wǎng)絡(luò)損失的樣本。

    2.3 銀行系統(tǒng)性風(fēng)險(xiǎn)度量

    由于銀行系統(tǒng)性風(fēng)險(xiǎn)的小概率特點(diǎn),因此研究銀行間網(wǎng)絡(luò)損失分布的尾部特征顯得格外重要?,F(xiàn)有同類研究都關(guān)注不同沖擊強(qiáng)度的差異、不同網(wǎng)絡(luò)結(jié)構(gòu)的差異等問題,卻忽略了沖擊的極端情況,即損失分布的尾部特征。由于沖擊是隨機(jī)的,沖擊造成的結(jié)果也是隨機(jī)的。所以對(duì)隨機(jī)結(jié)果的比較并不具有代表性,也不能反映系統(tǒng)性風(fēng)險(xiǎn)的小概率特點(diǎn)。

    本文提出利用風(fēng)險(xiǎn)價(jià)值VaR(Value at Risk)和期望損失ES(Expected Shortfall)兩個(gè)測度構(gòu)建銀行系統(tǒng)性風(fēng)險(xiǎn)的度量模型。同時(shí),利用違約銀行數(shù)量衡量銀行間網(wǎng)絡(luò)損失。銀行系統(tǒng)性風(fēng)險(xiǎn)VaR的含義可以概括為,在某一概率水平下,銀行間網(wǎng)絡(luò)的最大可能損失,即最壞情況時(shí)銀行間網(wǎng)絡(luò)的損失。用公式表示:

    P(Δ≤VaRα)=1-α

    (9)

    其中,Δ表示銀行間網(wǎng)絡(luò)損失的統(tǒng)計(jì)量,α表示顯著性水平,1-α則表示為置信水平。式(9)的含義是銀行間網(wǎng)絡(luò)損失小于VaRα的可能性為1-α。

    在實(shí)際應(yīng)用中,VaR不滿足次可加性,而且表現(xiàn)并不穩(wěn)定。所以本文同時(shí)構(gòu)建了另一個(gè)度量模型,銀行系統(tǒng)性風(fēng)險(xiǎn)ES。銀行系統(tǒng)性風(fēng)險(xiǎn)ES含義是,當(dāng)銀行間網(wǎng)絡(luò)損失超過VaR閥值時(shí)所遭受的平均損失程度。ES在VaR的基礎(chǔ)上進(jìn)一步考慮了出現(xiàn)極端情況時(shí)的平均損失程度,可以更完整地衡量銀行間網(wǎng)絡(luò)的極端損失風(fēng)險(xiǎn)。用公式表示:

    (10)

    銀行系統(tǒng)性風(fēng)險(xiǎn)VaR和銀行系統(tǒng)性風(fēng)險(xiǎn)ES衡量了銀行間網(wǎng)絡(luò)的極端損失風(fēng)險(xiǎn),體現(xiàn)了銀行系統(tǒng)性風(fēng)險(xiǎn)的小概率特點(diǎn)。從而能夠解決現(xiàn)有研究比較隨機(jī)沖擊結(jié)果無法反映系統(tǒng)性風(fēng)險(xiǎn)的問題。

    3 模擬實(shí)驗(yàn)

    3.1 銀行間網(wǎng)絡(luò)的參數(shù)校準(zhǔn)

    銀行間風(fēng)險(xiǎn)傳染是通過銀行間債權(quán)債務(wù)網(wǎng)絡(luò)實(shí)現(xiàn)的。所以,模擬實(shí)驗(yàn)之前需要對(duì)銀行間債務(wù)網(wǎng)絡(luò)的相關(guān)參數(shù)進(jìn)行設(shè)置和校準(zhǔn)。首先,構(gòu)造符合銀行間網(wǎng)絡(luò)特征的網(wǎng)絡(luò)結(jié)構(gòu),即確定鄰接矩陣。其次,確定每家銀行的同業(yè)資產(chǎn)和同業(yè)負(fù)債的總量。再次,根據(jù)同業(yè)資產(chǎn)和同業(yè)負(fù)債與總資產(chǎn)和股東權(quán)益的統(tǒng)計(jì)關(guān)系,確定每家銀行的總資產(chǎn)、股東權(quán)益、外部資產(chǎn)、存款等資產(chǎn)負(fù)債表數(shù)據(jù)。最后,根據(jù)鄰接矩陣、同業(yè)資產(chǎn)和同業(yè)負(fù)債,利用交叉熵估計(jì)銀行間債務(wù)矩陣。

    3.1.1 銀行間網(wǎng)絡(luò)結(jié)構(gòu)

    大量文獻(xiàn)表明,銀行間網(wǎng)絡(luò)是無標(biāo)度網(wǎng)絡(luò)[14-17]。為了保證模擬實(shí)驗(yàn)符合真實(shí)銀行間網(wǎng)絡(luò)結(jié)構(gòu)特征,本文根據(jù)文獻(xiàn)[6]提出的無標(biāo)度網(wǎng)絡(luò)構(gòu)建方法模擬了三種規(guī)模為200家銀行的無標(biāo)度網(wǎng)絡(luò),并利用鄰接矩陣分別表示三種銀行間網(wǎng)絡(luò)的結(jié)構(gòu)。根據(jù)網(wǎng)絡(luò)集中度從大到小排列,三種銀行間網(wǎng)絡(luò)的基本特征參數(shù)列入表1。

    表1 三種網(wǎng)絡(luò)的特征參數(shù)

    從表1可以看出,三種網(wǎng)絡(luò)的平均度和平均路徑長度差別不大,而網(wǎng)絡(luò)差異主要體現(xiàn)在聚集系數(shù)和集中度。聚集系數(shù)和集中度都是反映網(wǎng)絡(luò)的集中程度的指標(biāo)。本文以網(wǎng)絡(luò)集中度來區(qū)別三個(gè)銀行間網(wǎng)絡(luò)的差異。

    3.1.2 同業(yè)資產(chǎn)和同業(yè)負(fù)債

    Barrat 等[20]把與一個(gè)節(jié)點(diǎn)的所有連線的權(quán)重之和定義為節(jié)點(diǎn)強(qiáng)度。在銀行間網(wǎng)絡(luò)中,銀行i的同業(yè)資產(chǎn)IAi和同業(yè)負(fù)債ILi分別代表了銀行i的出度強(qiáng)度和入度強(qiáng)度,銀行i的債務(wù)銀行數(shù)量和債權(quán)銀行數(shù)量分別代表了銀行i的出度和入度。Soram?ki等[15]的研究表明銀行間網(wǎng)絡(luò)中節(jié)點(diǎn)度和節(jié)點(diǎn)強(qiáng)度符合下面的冪函數(shù)關(guān)系:

    s=αkβ, k∈[1,+∞)

    (11)

    由3.1.1模擬的三種銀行間網(wǎng)絡(luò)的鄰接矩陣,可以得到節(jié)點(diǎn)出度和入度。并利用式(11)確定銀行間網(wǎng)絡(luò)中的同業(yè)資產(chǎn)IAi和同業(yè)負(fù)債ILi。

    3.1.3 總資產(chǎn)和股東權(quán)益

    Soram?ki等[15]的研究還發(fā)現(xiàn)總資產(chǎn)與同業(yè)資產(chǎn)之間也存在著式(11)所表示的冪函數(shù)關(guān)系。本文利用取自Bankscope的2012年的中國的110家商業(yè)銀行的數(shù)據(jù)進(jìn)行了檢驗(yàn),得到了下面的回歸結(jié)果:

    adj_R2=0.9101

    (12)

    (13)

    其中,TAi表示總資產(chǎn),IAi、ILi分別表示同業(yè)資產(chǎn)和同業(yè)負(fù)債,Ei表示股東權(quán)益。括號(hào)內(nèi)是t統(tǒng)計(jì)量,adj_R2是調(diào)整后的R2統(tǒng)計(jì)量。

    根據(jù)式(12)和(13)確定銀行間網(wǎng)絡(luò)中每家銀行的總資產(chǎn)和股東權(quán)益。然后根據(jù)銀行資產(chǎn)負(fù)債表的平衡性確定銀行外部資產(chǎn)和存款,從而獲得銀行間網(wǎng)絡(luò)中每家銀行的完整資產(chǎn)負(fù)債表信息。

    3.1.4 債務(wù)矩陣的估計(jì)

    根據(jù)3.1.1確定的三種網(wǎng)絡(luò)的鄰接矩陣、3.1.2確定的同業(yè)資產(chǎn)與負(fù)債數(shù)據(jù),利用交叉熵方法可以估計(jì)債務(wù)矩陣,即式(1)。交叉熵方法可以表示為下面的優(yōu)化問題:

    (14)

    同時(shí),滿足約束條件:

    (15)

    其中,L為要估計(jì)的債務(wù)矩陣,lij∈L。E表示銀行間網(wǎng)絡(luò)的鄰接矩陣,eij∈E,eij等于0或1,eij=0表示銀行i沒有貸款給銀行j,eij=1表示銀行i有貸款給銀行j。

    利用信息熵中RAS算法求解式(14)和(15)可以獲得銀行間債務(wù)矩陣L。由于交叉熵方法引入了鄰接矩陣,從而保證了校準(zhǔn)后的銀行間網(wǎng)絡(luò)保持無標(biāo)度網(wǎng)絡(luò)特征,避免了最大熵方法只能獲得完全結(jié)構(gòu)網(wǎng)絡(luò)的問題。

    應(yīng)當(dāng)指出,本文中銀行間債務(wù)網(wǎng)絡(luò)相關(guān)參數(shù)的設(shè)置和校準(zhǔn)更加符合真實(shí)的銀行間債務(wù)網(wǎng)絡(luò)。第一,無標(biāo)度網(wǎng)絡(luò)設(shè)計(jì)符合真實(shí)銀行間網(wǎng)絡(luò)結(jié)構(gòu)特征。第二,同業(yè)資產(chǎn)、同業(yè)負(fù)債、總資產(chǎn)、股東權(quán)益等資產(chǎn)負(fù)債表信息是根據(jù)真實(shí)銀行業(yè)數(shù)據(jù)回歸結(jié)果確定的。第三,交叉熵方法估計(jì)的債務(wù)矩陣保持了無標(biāo)度網(wǎng)絡(luò)結(jié)構(gòu),避免了完全結(jié)構(gòu)網(wǎng)絡(luò)的問題。

    3.2 銀行間網(wǎng)絡(luò)的蒙特卡洛模擬

    3.2.1 銀行間網(wǎng)絡(luò)損失與違約傳染概率

    針對(duì)3.1建立的三種銀行間網(wǎng)絡(luò),首先根據(jù)不同的銀行外部資產(chǎn)損失的標(biāo)準(zhǔn)差τ,確定25種外部沖擊情景,其中τ∈(0,0.1]。其次,針對(duì)每一種沖擊情景,根據(jù)式(8)對(duì)每家銀行的外部資產(chǎn)EAi進(jìn)行10000次隨機(jī)抽樣。并利用清算支付向量式(6)、違約類型判別條件式(2)和(3),記錄每次沖擊的總違約銀行數(shù)量、基礎(chǔ)違約銀行數(shù)量、傳染違約銀行的量。最后,統(tǒng)計(jì)三種銀行間網(wǎng)絡(luò)在25種外部沖擊情景下的損失分布。

    根據(jù)國際清算銀行的定義,銀行系統(tǒng)性風(fēng)險(xiǎn)的一個(gè)特征是違約傳染發(fā)生連鎖反應(yīng),形成多米諾效應(yīng)。為了研究這一特征,本文將傳染違約銀行的數(shù)量達(dá)到10家(5%)以上,確定為銀行系統(tǒng)發(fā)生了違約傳染的連鎖反應(yīng)。在10000模擬中統(tǒng)計(jì)出現(xiàn)違約傳染連鎖反應(yīng)的頻率,并將其定義為違約傳染概率。圖1展示了三種銀行間網(wǎng)絡(luò)在25種沖擊情景下的違約傳染概率。其中,橫坐標(biāo)代表不同外部沖擊情景,縱坐標(biāo)代表違約傳染概率,三角代表網(wǎng)絡(luò)1,圓圈代表網(wǎng)絡(luò)2,星號(hào)代表網(wǎng)絡(luò)3。

    由圖1可以發(fā)現(xiàn),當(dāng)外部沖擊較小時(shí)(τ≤0.03),違約傳染概率很小,在三種銀行間網(wǎng)絡(luò)中幾乎不會(huì)發(fā)生違約傳染的連鎖反應(yīng)。隨著外部沖擊的加大,違約傳染概率也在逐漸增加。當(dāng)外部沖擊達(dá)到τ=0.05時(shí),違約傳染概率達(dá)到了20%。平均每5次沖擊,就會(huì)出現(xiàn)1次違約傳染的連鎖反應(yīng)。當(dāng)外部沖擊較大時(shí)(τ≥0.07),幾乎可以確定銀行間網(wǎng)絡(luò)必然會(huì)出現(xiàn)違約傳染連鎖反應(yīng)。更重要的是,我們發(fā)現(xiàn)網(wǎng)絡(luò)集中度越高發(fā)生違約傳染的概率越低。也就是說,銀行間網(wǎng)絡(luò)越集中,越不容易發(fā)生違約傳染。

    圖1 傳染違約發(fā)生的概率

    為了對(duì)比三種銀行間網(wǎng)絡(luò)的損失情況,我們統(tǒng)計(jì)的每種沖擊情景下10000次模擬中,銀行間網(wǎng)絡(luò)的平均損失程度(總違約銀行數(shù)量的平均值、基礎(chǔ)違約銀行數(shù)量的平均值、傳染違約銀行數(shù)量的平均值)。圖2展示了三種銀行間網(wǎng)絡(luò)不同沖擊情景下的銀行間網(wǎng)絡(luò)的平均損失程度,其中代表三種銀行間網(wǎng)絡(luò)的符號(hào)與圖1相同。

    由圖2可以發(fā)現(xiàn),三種銀行間網(wǎng)絡(luò)的總違約銀行數(shù)量的平均值、基礎(chǔ)違約銀行數(shù)量的平均值、傳染違約銀行數(shù)量的平均值幾乎沒有任何差別。平均來看,銀行間網(wǎng)絡(luò)的損失程度與網(wǎng)絡(luò)集中度沒有明顯關(guān)系。可見,在現(xiàn)有研究中用隨機(jī)的沖擊結(jié)果來比較銀行間網(wǎng)絡(luò)的損失程度是沒有意義的。

    圖2 平均違約銀行數(shù)量

    3.2.2 銀行間網(wǎng)絡(luò)的損失分布

    銀行間網(wǎng)絡(luò)的平均損失相似,并不代表損失分布是相同。接下來我們將考察三種銀行間網(wǎng)絡(luò)下的損失分布情況。我們選取了六種沖擊情景,τ=0.03,0.04,…,0.08。每一種沖擊情景都有10000次的模擬結(jié)果。我們利用每次模擬結(jié)果中總違約銀行數(shù)量衡量銀行間網(wǎng)絡(luò)的損失。六種沖擊情景下的不同銀行間網(wǎng)絡(luò)損失的描述性統(tǒng)計(jì)量被列入表2中。

    根據(jù)每種沖擊情景下10000次的模擬結(jié)果,畫出總違約銀行數(shù)量的頻率直方圖,見圖3。為了更加直觀的展示損失分布的特征,圖3中每個(gè)子圖還包括了銀行間網(wǎng)絡(luò)損失分布的概率密度估計(jì)曲線、正態(tài)分布密度曲線。概率密度估計(jì)曲線是根據(jù)銀行間網(wǎng)絡(luò)損失數(shù)據(jù)估計(jì)出來的密度曲線,它反映了損失數(shù)據(jù)的分布特征,在圖3中用虛線表示。正態(tài)分布密度曲線則是對(duì)應(yīng)銀行間網(wǎng)絡(luò)損失數(shù)據(jù)均值和標(biāo)準(zhǔn)差的正態(tài)分布密度曲線,它反映了正態(tài)分布特征,在圖3中用實(shí)線表示。圖3中的每一行子圖依次代表了不同的沖擊情景,τ=0.03,0.04,…,0.08。每一列子圖按照順序分別代表了三種銀行間網(wǎng)絡(luò)。

    表2 銀行間網(wǎng)絡(luò)損失的描述性統(tǒng)計(jì)量

    圖3 銀行間網(wǎng)絡(luò)損失的頻率分布直方圖

    從表2和圖3可以看出銀行間網(wǎng)絡(luò)損失分布的兩個(gè)明顯特征。第一,隨著沖擊強(qiáng)度的增加,銀行間網(wǎng)絡(luò)損失分布從近似正態(tài)分布轉(zhuǎn)變成尖峰厚尾分布,最后變成雙峰分布。當(dāng)外部沖擊較小時(shí)(τ=0.03),銀行間網(wǎng)絡(luò)損失分布近似正態(tài)分布。圖3中第一行子圖顯示,概率密度估計(jì)曲線與正態(tài)分布密度曲線相對(duì)接近。同時(shí),從表2第1行可以看出,此時(shí)損失分布的偏度接近0,峰度接近3,相對(duì)接近正態(tài)分布。從對(duì)圖1的分析可知當(dāng)τ=0.03,銀行間網(wǎng)絡(luò)不會(huì)出現(xiàn)違約傳染連鎖反應(yīng)??梢?,此時(shí)銀行間網(wǎng)絡(luò)中違約的銀行幾乎都是基礎(chǔ)違約。由于沖擊是隨機(jī)的,所以銀行間網(wǎng)絡(luò)損失會(huì)近似正態(tài)分布。圖3中第二、三、四行子圖顯示,概率密度估計(jì)曲線與正態(tài)分布密度曲線相比,峰部更尖尾部更厚。同時(shí),從表2第2、3、4行可以看出,損失分布的峰度明顯高于3,說明當(dāng)外部沖擊加大后(τ=0.04,0.05,0.06),銀行間網(wǎng)絡(luò)損失具有尖峰厚尾的分布特征。當(dāng)外部沖擊加大后,銀行間網(wǎng)絡(luò)出現(xiàn)違約傳染連鎖反應(yīng)的概率也相應(yīng)的增加。盡管基礎(chǔ)違約的銀行數(shù)量仍然近似正態(tài)分布,但是疊加了違約傳染連鎖效應(yīng)后導(dǎo)致了銀行間網(wǎng)絡(luò)損失呈現(xiàn)出尖峰厚尾分布特征。當(dāng)外部沖擊達(dá)到必然引發(fā)違約傳染連鎖效應(yīng)的程度后(τ=0.07,0.08),違約傳染徹底改變了銀行間網(wǎng)絡(luò)損失的分布形態(tài)。從圖3最后兩行子圖的概率密度估計(jì)曲線(虛線)可以明顯看出,銀行間網(wǎng)絡(luò)損失分布呈現(xiàn)雙峰分布特征。其中,左邊大峰部是基礎(chǔ)違約和傳染違約疊加后的結(jié)果,而右邊小峰部是傳染違約影響的結(jié)果。也就是說,此時(shí)銀行間的違約傳染導(dǎo)致銀行間網(wǎng)絡(luò)的極端損失(尾部分布)不再是小概率事件。可見,傳染違約是導(dǎo)致銀行系統(tǒng)性風(fēng)險(xiǎn)爆發(fā)的關(guān)鍵。

    第二,當(dāng)銀行間網(wǎng)絡(luò)損失分布呈現(xiàn)尖峰厚尾特征時(shí),網(wǎng)絡(luò)集中度越高銀行間網(wǎng)絡(luò)損失分布的厚尾現(xiàn)象越明顯。從表2第2、3、4行可以看出,當(dāng)銀行間網(wǎng)絡(luò)損失具有尖峰厚尾的分布特征時(shí)(τ=0.04,0.05,0.06),網(wǎng)絡(luò)集中度越高,銀行間網(wǎng)絡(luò)損失的峰度越大??梢姡y行間網(wǎng)絡(luò)集中度越高,損失分布的厚尾現(xiàn)象越明顯。但是,從表2第1、5、6行來看,銀行間網(wǎng)絡(luò)損失分布呈現(xiàn)近似正態(tài)分布(τ=0.03)和雙峰分布(τ=0.07,0.08)時(shí),峰度統(tǒng)計(jì)量的差異并不明顯。

    為了進(jìn)一步比較不同網(wǎng)絡(luò)結(jié)構(gòu)下銀行系統(tǒng)性風(fēng)險(xiǎn)的差異,下面將采用本文提出的兩種系統(tǒng)性風(fēng)險(xiǎn)測度(銀行系統(tǒng)性風(fēng)險(xiǎn)VaR和銀行系統(tǒng)性風(fēng)險(xiǎn)ES)進(jìn)行對(duì)比。

    3.3 銀行系統(tǒng)性風(fēng)險(xiǎn)的測度結(jié)果

    根據(jù)每種外部沖擊下的10000次隨機(jī)模擬的銀行間網(wǎng)絡(luò)損失結(jié)果,可以估計(jì)出對(duì)應(yīng)置信水平銀行系統(tǒng)性風(fēng)險(xiǎn)VaR值和銀行系統(tǒng)性風(fēng)險(xiǎn)ES值。表3列出了部分98%置信水平下的VaR估計(jì)結(jié)果和ES估計(jì)結(jié)果。表4列出了部分99%置信水平下的VaR估計(jì)結(jié)果和ES估計(jì)結(jié)果。

    表3和表4的前三列是三種網(wǎng)絡(luò)下VaR的估計(jì)結(jié)果。其含義是在對(duì)應(yīng)置信水平下銀行間網(wǎng)絡(luò)最大的損失(即總違約銀行數(shù)量)。第4到6列是對(duì)應(yīng)的傳染違約銀行數(shù)量。例如,表3第5行第1列數(shù)值是64,第5行第4列數(shù)值是24。這表示在第一種銀行間網(wǎng)絡(luò)中98%的概率下,最多會(huì)有64家銀行違約,其中傳染違約的銀行數(shù)量是24家。

    表3 銀行間系統(tǒng)性風(fēng)險(xiǎn)的測度結(jié)果(98%置信水平)

    表4 銀行間系統(tǒng)性風(fēng)險(xiǎn)的測度結(jié)果(99%置信水平)

    表3和表4的第7到9列是三種網(wǎng)絡(luò)下ES的估計(jì)結(jié)果。其含義是在對(duì)應(yīng)置信水平下銀行間網(wǎng)絡(luò)損失超過VaR閥值時(shí)的平均損失程度(即平均總違約銀行數(shù)量)。第10到12列是對(duì)應(yīng)的平均傳染違約銀行數(shù)量。例如,表3第5行第7列數(shù)值是95,第5行第10列數(shù)值是55。這表示在第一種銀行間網(wǎng)絡(luò)中98%的概率下,超過VaR閥值時(shí),平均會(huì)有95家銀行發(fā)生違約,其中傳染違約的銀行數(shù)量是55家。

    為了更加直觀,我們同時(shí)用圖形展示銀行系統(tǒng)性風(fēng)險(xiǎn)VaR和銀行系統(tǒng)性風(fēng)險(xiǎn)ES的估計(jì)結(jié)果。圖4中兩個(gè)子圖分別展示了98%置信水平下的VaR估計(jì)結(jié)果和ES估計(jì)結(jié)果。圖5對(duì)應(yīng)的是99%置信水平下估計(jì)結(jié)果。橫坐標(biāo)對(duì)應(yīng)的是不同的外部沖擊情景,縱坐標(biāo)對(duì)應(yīng)的是銀行系統(tǒng)性風(fēng)險(xiǎn)估計(jì)結(jié)果。圖中三角、圓圈、星號(hào)依次代表三種網(wǎng)絡(luò)。

    圖4 銀行間系統(tǒng)性風(fēng)險(xiǎn)VaR與ES(98%置信水平)

    圖5 銀行間系統(tǒng)性風(fēng)險(xiǎn)VaR與ES(99%置信水平)

    圖4和圖5中每個(gè)子圖都有三組線。最上面一組線代表了三種網(wǎng)絡(luò)下銀行系統(tǒng)性風(fēng)險(xiǎn)VaR(左邊子圖)和銀行系統(tǒng)性風(fēng)險(xiǎn)ES(右邊子圖),也就是銀行系統(tǒng)中總的違約銀行數(shù)量。中間的一組線表示傳染違約銀行數(shù)量,最下面一組線表示基礎(chǔ)違約銀行數(shù)量。

    從表3和表4可以看出,當(dāng)外部沖擊較小時(shí)(τ≤0.03),銀行間網(wǎng)絡(luò)的損失很小,而且在不同的網(wǎng)絡(luò)中沒有差異。τ=0.02,最多有1家銀行違約;τ=0.028最多有6到7家銀行違約。當(dāng)外部沖擊加大后(τ>0.03),銀行間網(wǎng)絡(luò)的損失增加,而且在不同的網(wǎng)絡(luò)中也體現(xiàn)出差異??偟谋憩F(xiàn)為,銀行間網(wǎng)絡(luò)越集中,損失越大。這種規(guī)律在銀行系統(tǒng)性風(fēng)險(xiǎn)ES上表現(xiàn)的更加穩(wěn)定和明顯。同時(shí)可以明顯的看出,不同網(wǎng)絡(luò)下銀行系統(tǒng)性風(fēng)險(xiǎn)的差異主要來自違約傳染的影響差異。

    圖4和圖5更直觀的展現(xiàn)出表3和表4的結(jié)論。網(wǎng)絡(luò)越集中,銀行系統(tǒng)性風(fēng)險(xiǎn)VaR和銀行系統(tǒng)性風(fēng)險(xiǎn)ES(總違約銀行數(shù)量)越大。而且在圖4和圖5中看出,基礎(chǔ)違約銀行數(shù)量差異不大(ES表現(xiàn)更明顯)。這再一次證明了銀行系統(tǒng)性風(fēng)險(xiǎn)在不同網(wǎng)絡(luò)中的差異是來自于違約傳染。銀行間網(wǎng)絡(luò)越集中,違約傳染的銀行數(shù)量越多。我們發(fā)現(xiàn)了一個(gè)有趣的結(jié)論,網(wǎng)絡(luò)集中度越高發(fā)生違約傳染連鎖反應(yīng)的概率越小(圖1),但是網(wǎng)絡(luò)集中度越高違約傳染連鎖反應(yīng)的破壞力越大(圖4、5)。

    從表3、4和圖4、5中還發(fā)現(xiàn)兩個(gè)結(jié)論。第一,外部沖擊在0.04<τ<0.06范圍內(nèi),銀行系統(tǒng)性風(fēng)險(xiǎn)有一個(gè)加速上升的過程。而且,銀行間違約傳染的連鎖反應(yīng)加速過程更加明顯。這說明外部沖擊對(duì)銀行系統(tǒng)性風(fēng)險(xiǎn)的貢獻(xiàn)不是線性的,而是近似指數(shù)型形式的加速過程。這個(gè)加速過程是由銀行間的違約傳染造成的。銀行間債務(wù)關(guān)系形成的銀行間網(wǎng)絡(luò),構(gòu)成了潛在的違約傳染渠道。在外部沖擊下,銀行間網(wǎng)絡(luò)的潛在傳染作用會(huì)極大的放大銀行系統(tǒng)的風(fēng)險(xiǎn),即爆發(fā)銀行系統(tǒng)風(fēng)險(xiǎn)。第二,外部沖擊在τ=0.06左右,銀行系統(tǒng)內(nèi)違約傳染影響力達(dá)到最大。本文模擬了很極端的外部沖擊,在τ=0.06時(shí),銀行間網(wǎng)絡(luò)會(huì)出現(xiàn)大量的銀行違約。外部沖擊再提高就會(huì)直接造成銀行基礎(chǔ)違約,所以傳染違約銀行會(huì)下降。但這并不代表違約傳染影響力會(huì)下降。

    4 結(jié)語

    本文提出了基于網(wǎng)絡(luò)模型的銀行系統(tǒng)性風(fēng)險(xiǎn)度量方法:銀行系統(tǒng)性風(fēng)險(xiǎn)VaR和銀行系統(tǒng)性風(fēng)險(xiǎn)ES。大量銀行違約引發(fā)銀行系統(tǒng)性風(fēng)險(xiǎn)是極為罕見的,可供研究的樣本非常少。為此,我們采用蒙特卡洛模擬方法,模擬銀行外部沖擊造成的銀行間網(wǎng)絡(luò)損失的大樣本。利用銀行間網(wǎng)絡(luò)損失大樣本估計(jì)銀行系統(tǒng)性風(fēng)險(xiǎn)VaR和銀行系統(tǒng)性風(fēng)險(xiǎn)ES,能夠有效的捕捉銀行系統(tǒng)性風(fēng)險(xiǎn)的尾部特征,解決對(duì)比隨機(jī)沖擊結(jié)果無法反映銀行系統(tǒng)性風(fēng)險(xiǎn)的問題。

    在模擬實(shí)驗(yàn)中,本文利用真實(shí)銀行間網(wǎng)絡(luò)結(jié)構(gòu)參數(shù),對(duì)模擬的三種銀行間網(wǎng)絡(luò)進(jìn)行校準(zhǔn),保證了研究結(jié)論真實(shí)性和可靠性。在對(duì)三種不同集中度的銀行間網(wǎng)絡(luò)的模擬實(shí)驗(yàn)中,我們發(fā)現(xiàn)了一些對(duì)監(jiān)管當(dāng)局頗有意義的結(jié)論。第一,隨著外部沖擊加大,銀行間網(wǎng)絡(luò)損失分布從近似正態(tài)分布轉(zhuǎn)變成尖峰厚尾分布,最后變成雙峰分布。這是由于外部沖擊加大引起了違約傳染的連鎖反應(yīng)。第二,網(wǎng)絡(luò)集中度越高發(fā)生違約傳染連鎖反應(yīng)的概率越小,但是網(wǎng)絡(luò)集中度越高違約傳染連鎖反應(yīng)的破壞力越大。第三,違約傳染是導(dǎo)致銀行系統(tǒng)性風(fēng)險(xiǎn)爆發(fā)的關(guān)鍵。銀行間網(wǎng)絡(luò)的潛在傳染作用會(huì)極大的放大銀行系統(tǒng)的風(fēng)險(xiǎn)。而且違約傳染效應(yīng)是呈指數(shù)增長的。

    [1] Iori G, Jafarey S, Padilla F.Systemic risk on the interbank market [J].Journal of Economic Behavior & Organization, 2006, 61(4): 525-542.

    [2] Allen F, Gale D.Financial contagion [J].Journal of Political Economy, 2000, 108(1): 1-33.

    [3] Nier E, Yang Jing, Yorulmazer T, et al.Network models and financial stability [J].Journal of Economic Dynamics and Control, 2007, 31(6): 2033-2060.

    [4] May R M, Arinaminpathy N.Systemic risk: The dynamics of model banking systems [J].Journal of the Royal Society Interface, 2010, 7(46): 823-838.

    [5] Gai P, Kapadia S.Contagion in financial networks [J].Proceedings:Macthematical,Physical and Engineering Sciences,2010,466(2120):2401-2423.

    [6] 隋聰, 遲國泰, 王宗堯.網(wǎng)絡(luò)結(jié)構(gòu)與銀行系統(tǒng)性風(fēng)險(xiǎn)[J].管理科學(xué)學(xué)報(bào), 2014,17(4): 57-70.

    [7] Wells S.UK interbank exposures: Systemic risk implications [J].Financial Stability Review, 2002, 13(12): 175-182.

    [9] Mistrulli P.E.Assessing financial contagion in the interbank market: Maximum entropy versus observed interbank lending patterns [J].Journal of Banking & Finance, 2011, 35(5): 1114-1127.

    [10] 黃聰, 賈彥東.金融網(wǎng)絡(luò)視角下的宏觀審慎管理——基于銀行間支付結(jié)算數(shù)據(jù)的實(shí)證分析[J].金融研究, 2010(4): 1-14.

    [11] 馬君潞, 范小云, 曹元濤.中國銀行間市場雙邊傳染的風(fēng)險(xiǎn)估測及其系統(tǒng)性特征分析[J].經(jīng)濟(jì)研究, 2007(1): 68-78.

    [12] 陳庭強(qiáng), 何建敏.基于復(fù)雜網(wǎng)絡(luò)的信用風(fēng)險(xiǎn)傳染模型研究[J].中國管理科學(xué), 2014, 22(11): 1-10.

    [13] 歐陽紅兵, 劉曉東.中國金融機(jī)構(gòu)的系統(tǒng)重要性及系統(tǒng)性風(fēng)險(xiǎn)傳染機(jī)制分析——基于復(fù)雜網(wǎng)絡(luò)的視角[J].中國管理科學(xué), 2015, 23(10): 30-37.

    [14] Boss M, Elsinger H, Summer M, et al.Network topology of the interbank market [J].Quantitative Finance, 2004, 4(6): 677-684.

    [15] Soram?ki K, Bech M L, Arnold J, et al.The topology of interbank payment flows [J].Physica A, 2007, 379(1): 317-333.

    [16] Becher C, Millard S, Soram?ki K.The network topology of CHAPS Sterling [R].Working Paper UK: Bank of England, 2008.

    [17] Krause A, Giansante S.Interbank lending and the spread of bank failures: A network model of systemic risk [J].Journal of Economic Behavior and Organization, 2012, 83(3): 583-608.

    [18] Sheldon G, Maurer M.Interbank lending and systemic risk: an empirical analysis for Switzerland [J].Swiss Society of Economics and Statistics, 1998, 134(4): 685-704.

    [19] Eisenberg L, Noe T H.Systemic risk in financial systems [J].Management Science, 2001, 47(2): 236-249.

    [20] Barrat A, Barthelemy M, Pastor-Satorras R,et al.The architecture of complex weighted networks [J].Proceedings of the National Academy of Sciences, 2004, 101(11): 3747-3752.

    A Network Perspective Measurement Method for Banking Systemic Risk

    SUI Cong1,2, TAN Zhao-lin1, WANG Zong-yao3

    (1.School of Finance, Dongbei University of Finance and Economics, Dalian 116025, China;2.Center for Commodity Markets and Behavioral Decision Research, Dongbei University of Finance and Economics, Dalian 116025, China; 3.Surry International Institute, Dongbei University of Finance and Economics, Dalian 116025, China)

    The interbank network is convenient for liquidity adjustment of the interbank market.But, the network configuration also increases risk contagion among the interbank market.However, the relationship between the network configuration and the systemic risk is a disputable issue in the research area.On one hand, the systemic risk is a typical small probability event which only happens in extreme circumstance.On the other hand, simulated interbank market network in the research area is different from real interbank market configuration.So building a realistic network model and researching the behavior of real network model in extreme circumstance are the key issues.Based on reality interbank network model, two evaluation parameters of the systemic risk: VaR and ES are presented in this paper.Firstly, Monte Carlo method is utilized to simulate the external impact of interbank system.Then, the systemic risk VaR and ES which can reflect the small probability characters of the systemic risk are estimated and the tail properties of interbank system loss are captured.Secondly, real bank parameters are utilized to calibrate three kinds of interbank network in simulation.Such a method ensures the reality and reliability of simulation results.Finally, three valuable conclusions are drawn: (1) External impact will trigger contagion.The interbank system loss will change from norm distribution to heavy tail distribution and then to bimodal distribution.(2) The contagion probability of high density interbank network is smaller than that of low density network, but the destruction is much higher.(3) The potential contagion will enlarge the systemic risk and default contagion effect will increase exponentially.A model which can evaluate the extent of the destruction of the systemic risk in extreme condition is presented.Furthermore, the simulation results comprehensively reveal the relationship between the network configuration and the systemic risk.

    systemic risk; interbank network; monte carlo simulation; VaR

    1003-207(2016)05-0054-11

    10.16381/j.cnki.issn1003-207x.2016.05.007

    2014-12-19;

    2015-12-07

    國家自然科學(xué)基金項(xiàng)目(71571034,61304180);教育部人文社會(huì)科學(xué)基金項(xiàng)目(12YJCZH211);遼寧省高等學(xué)校優(yōu)秀人才支持計(jì)劃資助項(xiàng)目(WJQ2015012)

    簡介:隋聰(1978-),男(漢族),遼寧沈陽人,東北財(cái)經(jīng)大學(xué)金融學(xué)院博士,副教授,碩士生導(dǎo)師,研究方向:銀行系統(tǒng)性風(fēng)險(xiǎn)、網(wǎng)絡(luò)模型、金融工程、銀行風(fēng)險(xiǎn)管理,Email: suicong2004@163.com.

    F830

    A

    猜你喜歡
    傳染系統(tǒng)性沖擊
    Our Mood Can Affect Others
    對(duì)于單身的偏見系統(tǒng)性地入侵了我們的生活?
    聽說,笑容是會(huì)“傳染”的
    傳染
    一類具有非線性傳染率的SVEIR模型的定性分析
    超聲引導(dǎo)經(jīng)直腸“10+X”點(diǎn)系統(tǒng)性穿刺前列腺的診療體會(huì)
    奧迪Q5換擋沖擊
    奧迪A8L換擋沖擊
    一汽奔騰CA7165AT4尊貴型車換擋沖擊
    巴菲特給我沖擊最大
    欧美日韩亚洲国产一区二区在线观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av五月六月丁香网| 99国产极品粉嫩在线观看| 久久天堂一区二区三区四区| 母亲3免费完整高清在线观看| 国产精品亚洲av一区麻豆| 真人做人爱边吃奶动态| 老汉色∧v一级毛片| 亚洲精品久久成人aⅴ小说| 国产精品爽爽va在线观看网站 | 亚洲成av人片免费观看| www.熟女人妻精品国产| 又紧又爽又黄一区二区| 午夜a级毛片| 日韩欧美免费精品| 成人三级做爰电影| 免费无遮挡裸体视频| 无遮挡黄片免费观看| 国产精品九九99| 午夜亚洲福利在线播放| 欧美黄色淫秽网站| 丝袜美足系列| 国产成人欧美| 女性生殖器流出的白浆| 国产国语露脸激情在线看| 人妻丰满熟妇av一区二区三区| 女同久久另类99精品国产91| 久久中文看片网| 1024视频免费在线观看| 色尼玛亚洲综合影院| 久久久久国产精品人妻aⅴ院| 黑人巨大精品欧美一区二区蜜桃| 黑丝袜美女国产一区| bbb黄色大片| 成人欧美大片| 波多野结衣高清无吗| 亚洲在线自拍视频| 亚洲成人免费电影在线观看| 午夜免费观看网址| 淫秽高清视频在线观看| 丝袜在线中文字幕| 性少妇av在线| 日本 欧美在线| 免费搜索国产男女视频| 免费高清视频大片| 欧美黑人精品巨大| 精品人妻在线不人妻| 国产精品精品国产色婷婷| 一本久久中文字幕| 国内精品久久久久精免费| 妹子高潮喷水视频| 亚洲久久久国产精品| 禁无遮挡网站| 两个人看的免费小视频| 国产精品国产高清国产av| 自线自在国产av| 黄色女人牲交| 亚洲精品久久国产高清桃花| 午夜免费鲁丝| ponron亚洲| 国产午夜精品久久久久久| 精品欧美国产一区二区三| 国产成人av激情在线播放| 日本黄色视频三级网站网址| 美女高潮到喷水免费观看| 久99久视频精品免费| 窝窝影院91人妻| 日本免费一区二区三区高清不卡 | av天堂久久9| 又紧又爽又黄一区二区| 亚洲精品国产精品久久久不卡| 操出白浆在线播放| 国产伦一二天堂av在线观看| av欧美777| a级毛片在线看网站| 丰满人妻熟妇乱又伦精品不卡| 免费无遮挡裸体视频| 欧美黑人精品巨大| 久久久久九九精品影院| 亚洲精品av麻豆狂野| 亚洲全国av大片| 嫩草影院精品99| 桃色一区二区三区在线观看| 手机成人av网站| 久久婷婷人人爽人人干人人爱 | 亚洲精品国产区一区二| 90打野战视频偷拍视频| 婷婷丁香在线五月| 一区二区三区激情视频| 色综合站精品国产| 97超级碰碰碰精品色视频在线观看| 国产伦一二天堂av在线观看| 美女免费视频网站| 桃红色精品国产亚洲av| 精品国产国语对白av| 亚洲自拍偷在线| 国产精品精品国产色婷婷| av有码第一页| 国产精品,欧美在线| 一进一出好大好爽视频| 精品乱码久久久久久99久播| 美女免费视频网站| 精品国产美女av久久久久小说| 宅男免费午夜| 亚洲中文av在线| av网站免费在线观看视频| 欧美日韩一级在线毛片| 法律面前人人平等表现在哪些方面| 成人欧美大片| 亚洲人成网站在线播放欧美日韩| 亚洲国产精品sss在线观看| 天堂动漫精品| av网站免费在线观看视频| 欧美色欧美亚洲另类二区 | 久久精品国产清高在天天线| 在线观看免费午夜福利视频| 在线十欧美十亚洲十日本专区| 国产主播在线观看一区二区| 欧美成狂野欧美在线观看| 国产1区2区3区精品| 亚洲少妇的诱惑av| 久久久久国内视频| 咕卡用的链子| 黄色女人牲交| 国内精品久久久久久久电影| av免费在线观看网站| 久久中文看片网| 亚洲人成电影免费在线| 757午夜福利合集在线观看| 久久婷婷人人爽人人干人人爱 | 亚洲色图av天堂| 99久久综合精品五月天人人| avwww免费| 亚洲 欧美 日韩 在线 免费| 免费看a级黄色片| 亚洲熟妇中文字幕五十中出| 母亲3免费完整高清在线观看| 欧美丝袜亚洲另类 | 十分钟在线观看高清视频www| 亚洲欧美日韩高清在线视频| 一二三四在线观看免费中文在| 亚洲国产毛片av蜜桃av| 亚洲人成77777在线视频| 无人区码免费观看不卡| 国产精品免费视频内射| or卡值多少钱| 日韩欧美国产在线观看| 亚洲全国av大片| 亚洲av第一区精品v没综合| 午夜久久久在线观看| 在线观看免费视频日本深夜| 国产精品日韩av在线免费观看 | 国产精品免费一区二区三区在线| av有码第一页| 一本大道久久a久久精品| 视频区欧美日本亚洲| 国产成年人精品一区二区| 免费高清在线观看日韩| 天天添夜夜摸| 国产精品精品国产色婷婷| 极品教师在线免费播放| a在线观看视频网站| 亚洲欧美精品综合久久99| 在线国产一区二区在线| 12—13女人毛片做爰片一| 国产精品,欧美在线| 午夜福利在线观看吧| 亚洲av电影不卡..在线观看| 亚洲午夜理论影院| 亚洲专区中文字幕在线| 88av欧美| 久久久久九九精品影院| 91大片在线观看| 伊人久久大香线蕉亚洲五| 亚洲欧美一区二区三区黑人| 99精品在免费线老司机午夜| 无遮挡黄片免费观看| 久久久久久久久久久久大奶| 久久国产精品人妻蜜桃| 亚洲国产毛片av蜜桃av| 国产精华一区二区三区| av电影中文网址| 午夜福利欧美成人| 一级作爱视频免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲人成伊人成综合网2020| 亚洲国产高清在线一区二区三 | 最近最新中文字幕大全免费视频| 法律面前人人平等表现在哪些方面| 搡老熟女国产l中国老女人| av在线播放免费不卡| 大陆偷拍与自拍| 十分钟在线观看高清视频www| 国产亚洲av高清不卡| 色综合站精品国产| 日韩欧美在线二视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲av熟女| 高潮久久久久久久久久久不卡| 91成人精品电影| 18禁裸乳无遮挡免费网站照片 | 精品福利观看| 成在线人永久免费视频| 9色porny在线观看| 午夜福利欧美成人| 精品无人区乱码1区二区| 亚洲欧美精品综合一区二区三区| 90打野战视频偷拍视频| 亚洲欧美日韩高清在线视频| 日韩欧美一区视频在线观看| 97超级碰碰碰精品色视频在线观看| 色婷婷久久久亚洲欧美| 精品国产国语对白av| 日韩国内少妇激情av| 老司机在亚洲福利影院| 一进一出抽搐gif免费好疼| 午夜福利欧美成人| 精品人妻在线不人妻| 久久久水蜜桃国产精品网| 色播亚洲综合网| 亚洲欧美精品综合久久99| 国产精品乱码一区二三区的特点 | e午夜精品久久久久久久| 在线国产一区二区在线| 在线播放国产精品三级| 亚洲一区高清亚洲精品| 一进一出抽搐动态| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品亚洲一级av第二区| 人人妻人人爽人人添夜夜欢视频| 免费在线观看黄色视频的| 午夜福利成人在线免费观看| 亚洲aⅴ乱码一区二区在线播放 | 亚洲,欧美精品.| 黄色视频,在线免费观看| 国产成人av教育| 精品久久久久久久人妻蜜臀av | 亚洲熟妇中文字幕五十中出| 如日韩欧美国产精品一区二区三区| 国产激情欧美一区二区| 一区二区三区高清视频在线| 人人妻人人爽人人添夜夜欢视频| 嫁个100分男人电影在线观看| svipshipincom国产片| 日本黄色视频三级网站网址| 中文字幕人妻丝袜一区二区| 精品欧美国产一区二区三| 国产人伦9x9x在线观看| 激情视频va一区二区三区| 亚洲久久久国产精品| www.www免费av| 欧美国产日韩亚洲一区| 久久久久久人人人人人| 丝袜美腿诱惑在线| 欧美中文日本在线观看视频| 97人妻天天添夜夜摸| 午夜免费观看网址| 日韩有码中文字幕| 99国产精品一区二区三区| 免费av毛片视频| 国产麻豆69| 9热在线视频观看99| 日韩欧美一区视频在线观看| 母亲3免费完整高清在线观看| 欧美国产日韩亚洲一区| 中文亚洲av片在线观看爽| 在线视频色国产色| 日韩精品中文字幕看吧| 精品久久久久久,| 日韩三级视频一区二区三区| 色哟哟哟哟哟哟| 19禁男女啪啪无遮挡网站| 日本a在线网址| 一区在线观看完整版| 国产成人啪精品午夜网站| 超碰成人久久| 午夜福利欧美成人| 亚洲国产精品成人综合色| 法律面前人人平等表现在哪些方面| 在线观看66精品国产| 一进一出好大好爽视频| 日韩大码丰满熟妇| 亚洲精品国产色婷婷电影| 久热爱精品视频在线9| 国产一卡二卡三卡精品| 两性夫妻黄色片| 99国产精品99久久久久| 午夜精品国产一区二区电影| 日韩欧美免费精品| 国产亚洲精品av在线| 国产一区二区三区综合在线观看| 亚洲在线自拍视频| 99热只有精品国产| 亚洲国产中文字幕在线视频| 国产xxxxx性猛交| 性少妇av在线| 午夜免费成人在线视频| 热99re8久久精品国产| 国产99白浆流出| 午夜福利在线观看吧| www日本在线高清视频| 国产一区二区三区在线臀色熟女| 热re99久久国产66热| 91精品三级在线观看| 国产人伦9x9x在线观看| 国产欧美日韩一区二区三| 精品日产1卡2卡| 最好的美女福利视频网| 性少妇av在线| 国产三级在线视频| 久久久久久亚洲精品国产蜜桃av| av中文乱码字幕在线| 在线观看免费午夜福利视频| 欧美一区二区精品小视频在线| 男人的好看免费观看在线视频 | 成人亚洲精品av一区二区| 人人妻,人人澡人人爽秒播| 亚洲色图 男人天堂 中文字幕| 一边摸一边做爽爽视频免费| 十八禁人妻一区二区| 国产精品99久久99久久久不卡| 精品国产乱码久久久久久男人| 男男h啪啪无遮挡| 99在线人妻在线中文字幕| 两个人看的免费小视频| 男女床上黄色一级片免费看| 亚洲在线自拍视频| 长腿黑丝高跟| 免费在线观看日本一区| 88av欧美| 免费观看精品视频网站| 19禁男女啪啪无遮挡网站| 精品一区二区三区av网在线观看| 欧美老熟妇乱子伦牲交| 人人澡人人妻人| 美女免费视频网站| 亚洲色图 男人天堂 中文字幕| 高清黄色对白视频在线免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 成人永久免费在线观看视频| 在线观看免费视频网站a站| 欧美成狂野欧美在线观看| 亚洲一区中文字幕在线| 亚洲精品美女久久av网站| 国产成人欧美在线观看| 久久精品91蜜桃| 国产一区二区三区在线臀色熟女| 国产精品 欧美亚洲| 精品国产亚洲在线| 日韩高清综合在线| 午夜激情av网站| 欧美成人性av电影在线观看| 老熟妇乱子伦视频在线观看| a级毛片在线看网站| 97人妻天天添夜夜摸| 欧美在线一区亚洲| 在线观看66精品国产| 色精品久久人妻99蜜桃| 一边摸一边做爽爽视频免费| 亚洲成av人片免费观看| 亚洲国产欧美日韩在线播放| 女警被强在线播放| 国产av在哪里看| 日本撒尿小便嘘嘘汇集6| 日韩精品青青久久久久久| 三级毛片av免费| 亚洲av成人不卡在线观看播放网| 欧美色欧美亚洲另类二区 | 国产极品粉嫩免费观看在线| 伊人久久大香线蕉亚洲五| 国产成人av教育| 日韩精品青青久久久久久| 热99re8久久精品国产| 嫩草影视91久久| 国产精品美女特级片免费视频播放器 | 天堂√8在线中文| 国产精品美女特级片免费视频播放器 | 国产xxxxx性猛交| a级毛片在线看网站| 日日爽夜夜爽网站| 美女免费视频网站| 男人操女人黄网站| 欧美乱色亚洲激情| 欧美老熟妇乱子伦牲交| 国产成人欧美在线观看| 一级a爱视频在线免费观看| 女性生殖器流出的白浆| 亚洲男人的天堂狠狠| 亚洲一区二区三区色噜噜| 国产亚洲精品久久久久5区| 国产成人av激情在线播放| 国产一区在线观看成人免费| 免费在线观看亚洲国产| 亚洲国产精品999在线| www国产在线视频色| 日本欧美视频一区| 91麻豆av在线| 色综合亚洲欧美另类图片| 最近最新中文字幕大全电影3 | 国产精品精品国产色婷婷| 美国免费a级毛片| 51午夜福利影视在线观看| 美女国产高潮福利片在线看| 99riav亚洲国产免费| 精品国产亚洲在线| 国产欧美日韩综合在线一区二区| 婷婷精品国产亚洲av在线| 欧美中文日本在线观看视频| 一区二区日韩欧美中文字幕| 级片在线观看| 国产黄a三级三级三级人| 亚洲中文日韩欧美视频| 欧美日韩乱码在线| 亚洲欧美精品综合一区二区三区| 成人手机av| 亚洲国产欧美日韩在线播放| 欧美色欧美亚洲另类二区 | 亚洲情色 制服丝袜| 欧美日本中文国产一区发布| 9191精品国产免费久久| 男女下面插进去视频免费观看| 亚洲欧美精品综合一区二区三区| 香蕉国产在线看| 18禁观看日本| 欧美中文日本在线观看视频| 黄色女人牲交| 午夜免费鲁丝| 成在线人永久免费视频| 欧美不卡视频在线免费观看 | 国产成人精品无人区| 一二三四在线观看免费中文在| 大香蕉久久成人网| 嫩草影视91久久| 性欧美人与动物交配| 99re在线观看精品视频| 脱女人内裤的视频| av在线天堂中文字幕| 色婷婷久久久亚洲欧美| 久久久久九九精品影院| 久久人人97超碰香蕉20202| 十分钟在线观看高清视频www| 美女午夜性视频免费| 国产高清videossex| 免费高清视频大片| 午夜福利影视在线免费观看| 国产精品久久视频播放| 国产高清视频在线播放一区| 欧美日韩黄片免| 日日干狠狠操夜夜爽| 九色亚洲精品在线播放| 午夜精品久久久久久毛片777| 亚洲精品在线观看二区| 亚洲男人的天堂狠狠| 男女做爰动态图高潮gif福利片 | 国产高清有码在线观看视频 | 亚洲国产精品sss在线观看| 久久久国产成人免费| 久久精品91蜜桃| 成人手机av| 欧美成人午夜精品| 一级,二级,三级黄色视频| 777久久人妻少妇嫩草av网站| 男女下面插进去视频免费观看| 无限看片的www在线观看| 变态另类丝袜制服| 国产精品久久电影中文字幕| 亚洲国产精品sss在线观看| 亚洲人成伊人成综合网2020| 每晚都被弄得嗷嗷叫到高潮| 一区二区三区精品91| 欧美精品亚洲一区二区| 国产一区二区三区综合在线观看| 色综合欧美亚洲国产小说| 乱人伦中国视频| 99久久国产精品久久久| 国产亚洲av高清不卡| 夜夜爽天天搞| 咕卡用的链子| 色综合站精品国产| 一级毛片精品| 美女高潮到喷水免费观看| 看免费av毛片| 亚洲第一青青草原| 亚洲av电影不卡..在线观看| 成人av一区二区三区在线看| 18禁裸乳无遮挡免费网站照片 | 国内精品久久久久久久电影| 18禁观看日本| 别揉我奶头~嗯~啊~动态视频| 精品一区二区三区视频在线观看免费| 久久久久久久久久久久大奶| 色老头精品视频在线观看| 最近最新中文字幕大全免费视频| 久久精品亚洲熟妇少妇任你| 日本免费a在线| 精品国产国语对白av| videosex国产| 窝窝影院91人妻| 美女高潮到喷水免费观看| 啦啦啦观看免费观看视频高清 | 精品国产乱子伦一区二区三区| 国产精品久久久久久精品电影 | 男人舔女人的私密视频| 精品午夜福利视频在线观看一区| 黄片小视频在线播放| 国产成人精品久久二区二区免费| 日日干狠狠操夜夜爽| 夜夜看夜夜爽夜夜摸| 精品一区二区三区四区五区乱码| 精品国产一区二区三区四区第35| 怎么达到女性高潮| 国产伦人伦偷精品视频| 成人特级黄色片久久久久久久| 亚洲欧美激情综合另类| 男女做爰动态图高潮gif福利片 | 黄片小视频在线播放| 看黄色毛片网站| 国产精品 国内视频| 少妇粗大呻吟视频| 欧美另类亚洲清纯唯美| 国产黄a三级三级三级人| 亚洲国产精品合色在线| 69精品国产乱码久久久| 男人舔女人下体高潮全视频| 亚洲在线自拍视频| 少妇裸体淫交视频免费看高清 | 91老司机精品| 免费看美女性在线毛片视频| 激情视频va一区二区三区| 97超级碰碰碰精品色视频在线观看| av在线播放免费不卡| 日韩大尺度精品在线看网址 | 999久久久精品免费观看国产| 欧美黄色淫秽网站| 男女床上黄色一级片免费看| 亚洲va日本ⅴa欧美va伊人久久| 久久午夜亚洲精品久久| 国产xxxxx性猛交| 日日夜夜操网爽| 人成视频在线观看免费观看| 久久人人97超碰香蕉20202| 久久伊人香网站| 女生性感内裤真人,穿戴方法视频| 99国产精品一区二区三区| 这个男人来自地球电影免费观看| 精品少妇一区二区三区视频日本电影| 国产色视频综合| 91老司机精品| 亚洲国产欧美日韩在线播放| 午夜精品久久久久久毛片777| 午夜视频精品福利| 久9热在线精品视频| 黄色 视频免费看| 操出白浆在线播放| 日韩精品青青久久久久久| 日本撒尿小便嘘嘘汇集6| 久久国产精品男人的天堂亚洲| 精品久久久久久成人av| 国产精品久久久久久精品电影 | 多毛熟女@视频| 亚洲成人国产一区在线观看| 精品国内亚洲2022精品成人| 淫秽高清视频在线观看| 免费在线观看视频国产中文字幕亚洲| 久久久久久久精品吃奶| 成年版毛片免费区| 男人的好看免费观看在线视频 | 正在播放国产对白刺激| 国产av精品麻豆| 亚洲中文字幕日韩| 精品午夜福利视频在线观看一区| 久久香蕉国产精品| 麻豆成人av在线观看| 操美女的视频在线观看| 国产精品一区二区三区四区久久 | 一区二区三区高清视频在线| 多毛熟女@视频| 极品教师在线免费播放| 欧美亚洲日本最大视频资源| 午夜福利成人在线免费观看| 91成年电影在线观看| 日韩中文字幕欧美一区二区| 可以免费在线观看a视频的电影网站| 国产一区二区三区综合在线观看| av免费在线观看网站| 国产亚洲欧美98| 91麻豆精品激情在线观看国产| 亚洲自偷自拍图片 自拍| 69精品国产乱码久久久| 久热这里只有精品99| 黄色毛片三级朝国网站| 十八禁人妻一区二区| 国产又爽黄色视频| 看黄色毛片网站| 久久九九热精品免费| 一区福利在线观看| 好看av亚洲va欧美ⅴa在| 国产熟女午夜一区二区三区| 国产免费av片在线观看野外av| 欧美乱色亚洲激情| 日韩有码中文字幕| 国产一卡二卡三卡精品| 女性被躁到高潮视频| 亚洲人成伊人成综合网2020| 亚洲成av片中文字幕在线观看| 免费无遮挡裸体视频| e午夜精品久久久久久久| 天堂√8在线中文| 又紧又爽又黄一区二区| 中出人妻视频一区二区| 久久久久久大精品| 日本黄色视频三级网站网址| 亚洲男人的天堂狠狠| 国产精品av久久久久免费|