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

    內(nèi)生性回收率與信用風(fēng)險度量研究

    2016-12-19 08:47:45吳建華王新軍
    中國管理科學(xué) 2016年1期
    關(guān)鍵詞:信用風(fēng)險度量回收率

    吳建華,王新軍,張 穎

    (1.濟南大學(xué)數(shù)學(xué)科學(xué)院, 山東 濟南 250022; 2.山東大學(xué)經(jīng)濟學(xué)院, 山東 濟南 250100)

    ?

    內(nèi)生性回收率與信用風(fēng)險度量研究

    吳建華1,王新軍2,張 穎1

    (1.濟南大學(xué)數(shù)學(xué)科學(xué)院, 山東 濟南 250022; 2.山東大學(xué)經(jīng)濟學(xué)院, 山東 濟南 250100)

    在信用風(fēng)險模型中,外生性回收率的設(shè)定會忽略回收率對損失分布尾部的影響,而且會導(dǎo)致潛在的模型風(fēng)險。本文將因子擴散過程引入結(jié)構(gòu)信用風(fēng)險模型,獲得了回收率和違約概率之間的內(nèi)在關(guān)系,利用Monte Carlo模擬方法數(shù)值分析了預(yù)期回收率對違約概率和資產(chǎn)價值波動率的依賴性,結(jié)果表明預(yù)期回收率與違約率之間具有很強的負相關(guān)關(guān)系,而且這種相關(guān)關(guān)系會受到債務(wù)人資產(chǎn)價值波動率的正向影響。在內(nèi)生性回收率下,推導(dǎo)了信用損失的概率分布,計算了信用風(fēng)險的Credit-VaR和ETF指標(biāo)。最后利用市場數(shù)據(jù)檢驗了內(nèi)生回收率信用風(fēng)險模型的有效性,結(jié)果表明該模型可以很好的描述歷史違約率和回收率的變化過程。

    內(nèi)生性回收率;因子擴散過程;信用風(fēng)險度量;數(shù)值模擬;實證檢驗

    1 引言

    對金融機構(gòu)、金融監(jiān)管和債權(quán)人來說,信用風(fēng)險度量一直是最為核心的內(nèi)容。尤其是對于結(jié)構(gòu)化信用衍生產(chǎn)品的定價和信用評級,信用風(fēng)險的度量更是一個基本的前提條件。然而,目前大多數(shù)現(xiàn)代信用風(fēng)險度量模型主要是圍繞違約概率展開,對于違約回收率的研究相對較少,從技術(shù)層面來看,違約概率的建模比較容易處理,而回收率的建模需要考慮更多的影響因素,比如借款人的信用質(zhì)量、債項的特征(債務(wù)工具優(yōu)先級)、抵押品質(zhì)量、行業(yè)分類、宏觀經(jīng)濟因素等的影響。在信用風(fēng)險管理中,違約回收率用來描述違約發(fā)生時債權(quán)人資產(chǎn)回收的程度。在商業(yè)銀行內(nèi)部評級高級法中,預(yù)期損失與非預(yù)期損失的計算以及對信貸業(yè)務(wù)的風(fēng)險調(diào)整度量都需要估計回收率風(fēng)險。同樣,在信用衍生產(chǎn)品和結(jié)構(gòu)化信用產(chǎn)品的定價計算中,也必須對標(biāo)的資產(chǎn)回收率進行估計。

    縱觀信用風(fēng)險度量模型的研究,主要以回收率的外生假設(shè)為主。從結(jié)構(gòu)化模型來看,雖然Merton[1]在第一次提出結(jié)構(gòu)化信用風(fēng)險模型時,就指出違約概率和回收率均由到期時企業(yè)的市場價值決定,但是Merton對于市場價值如何影響回收率的機制并沒有做進一步的說明。而且,在Black和Cox[2]將Merton模型推廣到首達時模型時,回收率不再設(shè)定為由企業(yè)市場價值決定,而是被獨立的模型化的,它被定義為未償付債務(wù)價值的固定比率,而且與違約概率相互獨立。之后的研究基本上是以外生性回收率作為前提,只是具體假設(shè)有所差異,比如Jokivuolle和Peura[3]的研究中假設(shè)外生的抵押品的價值是決定回收率的唯一隨機因素。Giesecke[4]的研究中回收率則直接被假設(shè)為常數(shù),Asvanunt 和Staal[5-6]的研究中回收率是由一個簡約化的方法決定的??傊^大多數(shù)結(jié)構(gòu)化模型中,無論回收率被直接假設(shè)為常數(shù),還是決定于某些因素,都是以回收率的外生性為特征的。

    從信用風(fēng)險度量的另一個路徑——簡約化模型來看,回收率則徹底被外生化。Jarrow和Turnbull[7]首次建立了一個離散形式的簡約化模型,隨后在Jarrow等[8]、Duffie和Singleton[9]的研究基礎(chǔ)上發(fā)展形成了簡約化信用風(fēng)險模型,在所有的簡約化模型中,只有在為了計算信用損失分布和進行資產(chǎn)定價的時候,才考慮回收率的問題。簡約化模型對回收率的處理方式可以概括為三種:等價回收率、市值比例回收率和面值比例回收率。Jarrow和Turnbull[1]假設(shè)當(dāng)違約發(fā)生時,債權(quán)人可以獲得固定比例的等價面值無違約風(fēng)險的零息債券。Duffie和Singleton[9]假設(shè)回收部分的價值與債券違約前的市場價值保持一個恒定的比率。Houweling和Vorst[10]將回收率看做違約債券面值的一部分。違約發(fā)生時,債權(quán)人可獲得基于面值的恒定比率的現(xiàn)金償還。這三種方法各具特點:等價回收率的形式更適合適用于公司債券;市值比例回收率的假定符合互換合約的結(jié)構(gòu)且易于拓展。面值比例回收率可以簡化模型在計算方面的復(fù)雜性。

    另外,從信用風(fēng)險的應(yīng)用模型來看,四大現(xiàn)代信用風(fēng)險度量應(yīng)用模型Credit Metrics、KMV、Credit Portfolio View和Credit Risk Plus中的違約回收率通常被看作是外生的常數(shù)或隨機變量,而且回收率和違約概率被看作兩個獨立變量。關(guān)于外生性回收率分布的假設(shè),F(xiàn)rye[11]建立了回收率的正態(tài)分布模型,Pykhtin[12]提出對數(shù)正態(tài)分布模型,Andersen和Sidenius[13]討論了正態(tài)模型、log正態(tài)模型、probit正態(tài)模型和logit-正態(tài)模型進行了對比研究。黃大海[14]對四種經(jīng)典的應(yīng)用模型的回收率情況進行了總結(jié),說明了Credit Metrics、KMV和Credit Portfolio View中的回收率通常假定服從單峰Beta分布,與違約概率互不相關(guān),Credit Risk+模型中的回收率則被假設(shè)為常數(shù)。王國棟和詹原瑞[15],汪飛星和姚磊[16]則進一步建立了回收率的雙峰Beta分布密度模型。陳暮紫等[17]利用廣義Beta回歸對影響不良貸款回收率的因素進行了實證分析,給回收率的實證研究提供了一個很好的計量經(jīng)濟學(xué)分析框架。

    關(guān)于回收率和違約概率之間的相關(guān)性關(guān)系實證研究方面存在一定的爭議。對于債券來說,Carty等[18]、Frye[19],Bakshi等[20],Hu和Perraudin[21],Cantor等[22],Carey和Gordy[23]以及Altman等[24]的研究發(fā)現(xiàn)對于企業(yè)債券來說,回收率與違約概率之間存在負相關(guān)關(guān)系。對于貸款來說,由于抵押品的存在,只要抵押品的價值相對穩(wěn)定,即使由于外在經(jīng)濟環(huán)境的影響導(dǎo)致違約概率上升,貸款的回收率也可能不會下降。但是Hamilton[25],Hu Yen-ting 等[26]通過實際的數(shù)據(jù)研究顯示,當(dāng)總的貸款違約概率較高時,回收率與違約概率之間還是存在一定的負相關(guān)性。 對于信用衍生產(chǎn)品來說,Hull和White[27]研究表明,衍生品CDO和CDS的隱含回收率與PD之間都存在顯著的負相關(guān)關(guān)系。Das和Hanouna[28]利用 CDS價差數(shù)據(jù)進行了實證分析,發(fā)現(xiàn)隱含回收率和PD呈負相關(guān)關(guān)系,尤其在PD很高時,這種負相關(guān)性更加明顯。

    綜上所述,在信用風(fēng)險度量模型中,有關(guān)回收率的研究主要集中于回收率的影響因素、回收率的分布模型、回收率與違約率的相關(guān)性等方面。目前的研究主要包括兩個方面的不足,一是回收率的外生性假設(shè)。絕大多數(shù)信用風(fēng)險模型普遍假設(shè)回收率與違約概率相互獨立,實證研究表明這會嚴重低估實際的巨額損失。而且這種主觀的假設(shè)會導(dǎo)致潛在的模型設(shè)定風(fēng)險;二是回收率分布的假設(shè),多數(shù)模型假設(shè)回收率服從貝塔分布。但是實證研究顯示,基于貝塔分布的信用VaR估計的誤差有增大的趨勢。

    然而,無論是從經(jīng)濟學(xué)直覺還是從金融實踐來看,違約資產(chǎn)的回收率都不應(yīng)該是外生的,事實上,如果某項資產(chǎn)在違約之后具有較低的回收率,那么這就意味著該項資產(chǎn)必然具有較高的違約率,反之,較高的回收率必然意味著債務(wù)人的以及債項的信用質(zhì)量較高,在正常的經(jīng)濟環(huán)境下,不會輕易的發(fā)生違約。實際上,從大量關(guān)于回收率的實證研究來看,違約概率和回收率之間不是相互獨立的。

    從上面這些文獻可知,無論是債券、貸款還是信用衍生品,回收率和違約概率之間都存在明顯的負相關(guān)關(guān)系。但是以上這些文獻只是說明了回收率和違約概率之間的負相關(guān)關(guān)系,并沒有具體說明這種負相關(guān)關(guān)系之間的內(nèi)在機理。為此,需要考慮如何將回收率進行內(nèi)生化分析。目前關(guān)于這方面的研究還較少,僅有的文獻有,Guo Xin等[29]提出在傳統(tǒng)簡約模型的框架中加入一個反映回收率的隨機過程。模型中公司違約過程、回收率過程和風(fēng)險債務(wù)定價的定量化都是基于公司的資產(chǎn)價值。Chava等[30]提出,在簡約模型中假設(shè)違約概率和回收率同時依賴于單一的公共因子,從而保留了違約和回收率之間的負相關(guān),但是該模型的缺點是待估參數(shù)過多,而且有些參數(shù)從經(jīng)濟學(xué)的角度難以解釋。不過,以上這兩篇文章對于回收率內(nèi)生化的研究,主要是在簡約模型框架內(nèi)進行的,對于結(jié)構(gòu)模型框架內(nèi)的回收率內(nèi)生化機理的研究還是個空白。

    Merton曾提出違約概率和回收率均由到期時企業(yè)的市場價值決定。但是Merton對于市場價值如何影響回收率的機制并沒有做進一步的說明。本文在結(jié)構(gòu)模型框架內(nèi),研究了多資產(chǎn)內(nèi)生性回收率的機理,在假設(shè)信用資產(chǎn)服從因子擴散過程的假設(shè)下,推導(dǎo)了回收率和違約概率之間的內(nèi)在函數(shù)關(guān)系,該函數(shù)關(guān)系僅僅依賴于一個參數(shù),我們稱之為內(nèi)生回收率。內(nèi)生回收率在保持違約率同回收率之間的負相關(guān)關(guān)系的前提下,避免了外生性回收率假設(shè)所存在的潛在的模型風(fēng)險。在內(nèi)生回收率框架內(nèi),推出了信用損失的概率分布,計算了兩個重要的信用風(fēng)險度量指標(biāo)C-VaR和ETF。利用Monte Carlo模擬數(shù)值分析,檢驗了內(nèi)生回收率的有效性。最后,利用Moody公司的部分數(shù)據(jù),檢驗了本文模型的有效性。

    2 回收率的內(nèi)生性分析

    考慮具有K個企業(yè)的債務(wù)組合,企業(yè)k(k=1,…,K)的市場價值為Vk。出于分析的簡化,假設(shè)企業(yè)資產(chǎn)結(jié)構(gòu)由所有權(quán)益和負債組成,企業(yè)k發(fā)行一種特定的零息債券,該債券的面值為Fk=F,持有期為T。假設(shè)企業(yè)資產(chǎn)市值的初始值為Vk(0)=V(0)。當(dāng)企業(yè)資產(chǎn)價值在T時刻低于面值,即Vk

    (1)

    令fVk(v)表示在到期日T,企業(yè)k的市場價值Vk的概率密度函數(shù)。那么違約概率為:

    (2)

    預(yù)期回收率為:

    (3)

    由(2)(3)式可知,違約率和回收率都是由T時刻企業(yè)的資產(chǎn)價值Vk決定的。從(3)可以看出,預(yù)期回收率E(Rk)受到違約率PD,k的直接影響,同現(xiàn)有的大多數(shù)信用風(fēng)險模型中兩者是獨立的假設(shè)相比,違約率和回收率之間的這種函數(shù)關(guān)系更貼近實際情況,我們稱之為內(nèi)生性回收率。值得一提的是,本文突出了回收率的“內(nèi)生”特征,這一點從式子(1)和(3)中可以清楚的看到,無論是回收率Rk還是預(yù)期回收率E(Rk)都是資產(chǎn)價值Vk的函數(shù),這是同已有的外生性回收率的設(shè)定本質(zhì)不同的地方,而三種常見的外生性回收率:等價回收率、市值比例回收率和面值比例回收率的設(shè)定都是在債權(quán)債務(wù)契約簽訂之前就外生設(shè)定的,而與債務(wù)契約簽訂之后的資產(chǎn)的價值是沒有任何關(guān)系的。

    3 基于因子擴散模型的內(nèi)生回收率

    3.1 企業(yè)價值的因子擴散模型

    在經(jīng)典的結(jié)構(gòu)模型中,通常假設(shè)企業(yè)的市場價值V(t)服從下面的擴散過程:

    dV(t)=V(t)μdt+V(t)σdZ

    (4)

    其中μ為漂移項,σ為波動項,dZk服從標(biāo)準(zhǔn)布朗運動。但是,這種經(jīng)典設(shè)定沒有說明造成市場價值波動的原因。從經(jīng)濟學(xué)常識我們可知,影響市場價值波動的原因無非來源于兩個方面,一是系統(tǒng)性的宏觀經(jīng)濟環(huán)境的改變;二是非系統(tǒng)性的個別企業(yè)內(nèi)部經(jīng)營管理的微觀層面的變動。因此,在擴散過程中引入公共的宏觀經(jīng)濟因子和企業(yè)自身的特殊因子,可以更好的刻畫企業(yè)價值的變化。另外,在信用資產(chǎn)組合中,不同企業(yè)之間的相關(guān)性也會影響到整個資產(chǎn)組合價值的變動,表現(xiàn)為企業(yè)之間的違約相關(guān)性,通過因子模型可以方便的描述不同企業(yè)之間的交互作用,這種相關(guān)性的描述在經(jīng)典的Merton模型中沒有給出。

    假設(shè)企業(yè)k(k=1,…,K)的市場價值Vk是由公共的宏觀經(jīng)濟因子M和個別企業(yè)的特殊因子Zk所決定的,即:

    (5)

    其中M和Zk均服從標(biāo)準(zhǔn)正態(tài)分布,Zk,k=1,…,K之間相互獨立,每一個Zk與M也相互獨立。c∈[-1,1]為常數(shù),這意味著,任何兩個企業(yè)之間的具有相同的相關(guān)結(jié)構(gòu)。這樣,我們就可以假設(shè)企業(yè)k的市場價值Vk(t)在t時刻的變動過程服從下面的隨機微分方程:

    (6)

    上式描述了一個具有線性相關(guān)特征的因子擴散過程,μ為漂移項,σ為波動項和c為相關(guān)系數(shù)。dM和dZk分別表示市場的波動性和企業(yè)k的波動性,都服從標(biāo)準(zhǔn)維納過程。顯然,由(6)式可知,不同的宏觀經(jīng)濟變量M和企業(yè)的特殊狀況Zk的實際值,會形成不同的市場價值Vk的概率密度函數(shù)fVk(v),從而根據(jù)(2)(3)我們可以計算得到不同的違約概率PD,k和預(yù)期回收率E(Rk)。

    3.2 內(nèi)生回收率的理論分析

    把債券的期限T等分為N段,Δt=T/N為離散時間的增量,這樣我們就可以得到(6)式的離散形式,對任意的t∈[0,T]:

    (7)

    其中η和εk均為標(biāo)準(zhǔn)正態(tài)分布,分別表示來自市場方面的隨機沖擊和企業(yè)自身經(jīng)營的不確定性帶來的沖擊。對(7)利用迭代運算可得在到期日T,企業(yè)k的市場價值服從如下的規(guī)律:

    (8)

    令Xk=(Vk(T)-V(0))/V(0)為企業(yè)k的市場收益率,所有K個企業(yè)在[0,T]內(nèi)的市場平均收益率記為Xm,則:

    (9)

    (10)

    (11)

    (12)

    對于企業(yè)k,利用(8)(9)(11)我們有:

    (13)

    (14)

    (15)

    由此可以得到單個企業(yè)的違約概率如下:

    (16)

    (17)

    企業(yè)的平均違約損失L=1-V/F的期望值可以如下計算:

    (18)

    從而有:

    (19)

    根據(jù)(17)我們可以將A(Xm)表達成PD的函數(shù):

    (20)

    這樣,我們就可以得到預(yù)期回收率和違約概率的函數(shù)依賴關(guān)系:

    (21)

    這就是內(nèi)生性回收率如何受到違約概率影響的內(nèi)在機理。此外,從方程(21)還可以看到預(yù)期回收率與違約率的函數(shù)關(guān)系受外生變量B的結(jié)構(gòu)性影響。為了更加直觀的理解這一關(guān)系,在圖1中給出不同的波動率下B,違約概率對預(yù)期回收率影響的變化情況。

    圖1 不同波動率下,違約概率對預(yù)期回收率的影響,B=0.1,0.2,0.3,0.4,0.5

    由圖1可知,隨著違約概率PD的增加,預(yù)期回收率對違約概率變得更加敏感,當(dāng)違約概率足夠大,趨向于1時,預(yù)期回收率會迅速趨向于0,表明預(yù)期回收率對違約概率具有高度的負的依賴性。另外,債務(wù)人資產(chǎn)的波動性B越高,預(yù)期回收率E(R(PD))就整體越高,表明資產(chǎn)價值的波動性不但影響違約概率,對于回收率的整體變動具有系統(tǒng)性的影響,并且表現(xiàn)出較強的正相關(guān)性。無論是違約概率還是波動性,它們都對預(yù)期回收率具有內(nèi)在的影響。這充分說明,回收率不應(yīng)該被人為地設(shè)定為外生的,回收率的外生性設(shè)定會帶來嚴重的模型風(fēng)險。

    3.3 對內(nèi)生回收率的數(shù)值分析

    在MC模擬中,我們考慮了方程(8)中的離散時間隨機過程。由于我們只考慮在T時刻的債務(wù)的市場價值,因此在進行模擬時,我們直接取了一大步,即Δt=T,故N=1。因此,市場價值的離散形式為:

    (22)

    每次給定一個市場波動的實現(xiàn)值η,我們模擬K=5000次得到5000個不同的εk實現(xiàn)值。對于每一次模擬的運行,我們計算市場平均收益率Xm,違約概率PD(Xm)和預(yù)期回收率E(R(Xm))。市場平均收益率Xm被定義為到期日T時刻的平均收益:

    (23)

    對于足夠大的K來說,異質(zhì)項εk會達到平均水平,而市場平均收益率Xm由η的實現(xiàn)值單獨的定義。這就是為什么我們利用市場平均收益率作為其他觀測值的參數(shù)的原因。對于違約概率做如下估計:

    PD(Xm)≈ND(Xm)/K

    (24)

    其中ND(Xm)是違約次數(shù),取值為事件Vk(T)

    (25)

    進而得到預(yù)期損失率:

    (26)

    從而可以估計預(yù)期回收率為:

    E(R(Xm))=1-E(L(Xm))

    (27)

    這里,我們假設(shè)違約的次數(shù)嚴格的非零,這對于足夠大的組合規(guī)模K來說是合理的。對市場平均收益率項的實現(xiàn)值模擬106次,這樣我們就可以得到市場平均收益率Xm,違約概率PD(Xm)和預(yù)期回收率E(R(Xm))的不同的值。

    不失一般性,假設(shè)市場價值的初始值為V0=100,零息債券的面值為F=75,到期日為T=1(比如1年)。相關(guān)系數(shù)假設(shè)為c=0.5,這同股票價格的常見的相關(guān)性是一致的。擴散過程的參數(shù)設(shè)為μ=0.05,σ=0.15,這意味著每年的平均增長率5%和波動率15%。對應(yīng)的杠桿率F/V0=75%。帶入(22)可得:

    模擬結(jié)果畫在圖2中,同時給出了預(yù)期回收率和預(yù)期損失率對違約概率的依賴性。

    圖2 預(yù)期回收率ER和預(yù)期損失率EL對違約概率PD的依賴性

    在圖2中,左圖顯示了預(yù)期回收率E(R(Xm))對違約概率PD(Xm)的依賴性,右圖顯示了預(yù)期損失率E(L(Xm))對違約概率PD(Xm)的依賴性,其中紅色的散點表示MC模擬結(jié)果,藍色的線表示分析結(jié)果。分析結(jié)果中的波動性B是由模擬結(jié)果獲得的參數(shù)。在這兩個例子中,我們分別觀察到了在MC模擬和方程(21)中的分析結(jié)果之間的高度的一致性。而且對于高違約概率,對于均值的偏離是較小的。這一點同實際的回收率是相符的,在后面第五部分的實證分析中可以得到驗證。由此可知,本文給出的基于因子擴散模型的回收率可以很好的刻畫回收率函數(shù)的內(nèi)生性特征。

    4 基于內(nèi)生回收率的信用風(fēng)險度量

    信用風(fēng)險的度量需要計算三個基本核心的參數(shù):違約概率(PD)、違約損失率 (LGD)和違約風(fēng)險敞口 (EAD)。利用這三個基本參數(shù),就可以得到任何組合債務(wù)潛在的信用損失CL,從而進一步的根據(jù)信用損失的概率分布獲得信用在險價值(Credit Value at Risk,C-VaR)和預(yù)期尾部損失(Expected Tail Loss,ETL)。

    首先討論信用損失CL的概率分布。對于具有K份合約的債務(wù)組合說,其潛在的信用損失為:

    (28)

    其中PDk、EADk和LGDk分別為債務(wù)k的違約概率、違約風(fēng)險敞口和違約損失率。對于齊次債務(wù)組合,有CL=K·PD·EAD·LGD。出于簡化,記EAD=F,LGD=L,故K份合約的齊次債務(wù)組合的潛在的信用損失為:

    CL(Xm)=K·F·PD(Xm)·L(Xm)

    (29)

    下面給出信用損失CL的概率分布。由CL的表達式可知,K,F(xiàn)都是常數(shù),而PD(Xm)和L(Xm)都是市場平均收益率Xm的函數(shù),因此,CL的概率分布fCL(l)dl同市場平均收益率的概率分布fXm(xm)dxm等價的,即

    fCL(l)dl?fXm(x)dx

    (30)

    故信用損失CL的概率分布函數(shù)為:

    (31)

    根據(jù)式子(12)可得:

    (32)

    給定組合信用損失和市場平均收益率之間的函數(shù)關(guān)系,CL(Xm)=K·F·PD(Xm)·L(Xm),我們可以將市場平均收益率Xm的密度函數(shù)轉(zhuǎn)換為信用損失CL的密度函數(shù):

    (33)

    圖3給出了信用損失CL的理論損失分布的結(jié)果和MC結(jié)果的對比。

    圖3 信用損失的理論分布和MC模擬

    從圖3中,我們觀察到信用損失的理論分布和MC模擬的高度的一致性,即使對于極端大的組合損失也是如此。下面利用上面推導(dǎo)的信用風(fēng)險損失分布來計算兩個重要的信用風(fēng)險度量指標(biāo)。

    在金融風(fēng)險管理中,VaR是目前各大金融機構(gòu)主流的風(fēng)險度量方法,而且VaR的應(yīng)用功能已經(jīng)拓展到了彌補潛在損失所需要的經(jīng)濟資本、銀行的資本充足率等資本要求的計算中,但是,作為風(fēng)險管理的主流工具的VaR并不滿足風(fēng)險度量的一致性要求[32],事實上,雖然VaR滿足單調(diào)性、正齊次性和平移不變性,但是它不滿足次可加性,因而無法充分體現(xiàn)資產(chǎn)組合的風(fēng)險分散效應(yīng),這也是本文引入預(yù)期尾部損失ETL作為對VaR信用風(fēng)險度量指標(biāo)補充的原因。

    C-VaR=WCL-ECL

    (34)

    (35)

    預(yù)期尾部損失(Expected Tail Loss,ETL)是Acebi和Tache[33]提出的風(fēng)險度量方法,它彌補了VaR方法不滿足次可加性的缺陷,是一個具有一致性的風(fēng)險度量指標(biāo)。如果X的分布函數(shù)FX(x)是連續(xù)的,在給定置信水平α下,ETL定義為:

    (36)

    對于信用損失CL來說,預(yù)期尾部損失C-ETL為:

    (37)

    下面利用MC技術(shù)進行數(shù)值計算??紤]一個由500個信用資產(chǎn)構(gòu)成的同質(zhì)的債務(wù)組合,假設(shè)市場價值的初始值為V0=100,零息債券的面值為F=75,到期日為T=1(比如1年)。假設(shè)任何兩個資產(chǎn)之間的相關(guān)系數(shù)為c=0.5,這同股票價格的常見的相關(guān)性是一致的。擴散過程的參數(shù)設(shè)為μ=0.05,σ=0.15,這意味著每年的平均增長率5%和波動率15%。對應(yīng)的杠桿率F/V0=75%。表1給出基于0.99置信水平的預(yù)期信用損失ECL,C-VaR和C-ETL的數(shù)值計算的結(jié)果。

    表1 ECL、C-VaR和ETL的數(shù)值模擬和解析結(jié)果

    注:置信水平α=0.99

    從表1可以看出,對于預(yù)期信用損失(ECL),信用在險價值(C-VaR)和預(yù)期尾部損失(ETL),解析結(jié)果同模擬結(jié)果高度一致,這說明本文所給出的基于內(nèi)生性回收率的信用風(fēng)險的度量模型,能夠很好的對信用風(fēng)險進行度量。

    5 內(nèi)生性回收率的實證檢驗

    由于國內(nèi)還沒有長期完整的債券信用違約數(shù)據(jù)庫,本文使用了山東大學(xué)金融實驗室提供的部分國外的債券信用評級數(shù)據(jù)。我們對于2000年1月1日到2011年12月31 的債券信用評級數(shù)據(jù),根據(jù)每一個債務(wù)人的違約風(fēng)險進行排序,從而構(gòu)建出同質(zhì)的信用組合。在我們的分析中,我們考慮了組合內(nèi)不同到期日和相應(yīng)的回收率的違約概率。

    我們測算了2000年1月1日到2010年12月31日的1年期債券的違約概率和回收率。第一段時期從2000年1月1日到2000年12月31日。第二段為2000年2月1日到2001年2月31日,以此類推,最后一段時期為2011年1月1日到2011年12月31日。

    由于違約事件是稀有事件,為了盡可能的獲得更多的違約數(shù)據(jù),我們考慮了初始評級較低的投機級信用資產(chǎn)組合,這些評級表明相關(guān)的信用資產(chǎn)是高風(fēng)險的。我們統(tǒng)計估算了評級為Caa1,Caa2和Caa3的優(yōu)先擔(dān)保債券的違約概率和回收率,如圖4。

    圖4 評級為Caa1,Caa2和Caa3的優(yōu)先擔(dān)保債券的違約概率和回收率

    從圖4可以看出,違約概率PD和違約回收率RR二者之間表現(xiàn)出了高度的負相關(guān)。實際的數(shù)據(jù)表明,回收率同違約概率之間是具有內(nèi)在的聯(lián)系的,而不是像大多數(shù)文獻所給出的獨立性假設(shè),故回收率是具有內(nèi)生性的特制的。

    下面我們選擇了特定評級的債券組成我們的組合,假設(shè)它們在構(gòu)成組合的相同時刻發(fā)行的。理論的預(yù)期回收率是根據(jù)式子(21)計算的,這實際上是對潛在的隨機過程實現(xiàn)值的一個平均。因此,我們對實際回收率的數(shù)據(jù)進行平均來描述預(yù)期回收率,然后觀察兩者的接近程度。為此,我們將實際的違約率數(shù)據(jù)等分成30組,然后計算每組的回收率,最后進行平均得到整個資產(chǎn)組合的平均回收率。計算結(jié)果畫在圖5中。

    圖5 兩年期和四年期違約概率與平均損失及回收率

    我們用這些數(shù)據(jù)來檢驗本文的內(nèi)生性回收率模型,從圖5 可以發(fā)現(xiàn)模型的理論結(jié)果同實際的數(shù)據(jù)高度的一致。內(nèi)生性回收率模型能夠更好的描述經(jīng)驗數(shù)據(jù)。該模型可以作為現(xiàn)存的違約概率模型的一個很好的補充。實際的債務(wù)人的資本結(jié)構(gòu)比Merton模型框架內(nèi)的更加復(fù)雜,而且違約通常在到期日之前就發(fā)生了。另外,如果我們已經(jīng)知道違約概率的分布,也可以利用方程(37)來描述組合損失分布。

    6 結(jié)語

    違約和回收率之間的交互性對大額的信用組合損失具有關(guān)鍵的影響。然而在現(xiàn)有的信用風(fēng)險模型中,違約概率和回收率經(jīng)常被獨立的模型化。本文基于經(jīng)典的擴散過程,推導(dǎo)了違約和回收率之間的函數(shù)關(guān)系。該函數(shù)關(guān)系是由一個單參數(shù)決定的。MC模擬實驗直觀的揭示了違約以及參數(shù)如何影響該內(nèi)生回收率的變化。在內(nèi)生回收率框架下,本文計算了度量信用風(fēng)險的兩個重要的指標(biāo)Credit-VaR和ETF指標(biāo)。最后利用實際數(shù)據(jù)檢驗了本文所提出的內(nèi)生回收率信用風(fēng)險模型的有效性,結(jié)果表明內(nèi)生回收率信用風(fēng)險模型可以更好的描述歷史違約率的變化過程。

    本文的研究僅僅考慮了信用風(fēng)險度量中經(jīng)典的結(jié)構(gòu)模型,對于更多的結(jié)構(gòu)模型也可以利用本文的思路進行拓展。顯然本文所提出的內(nèi)生性回收率模型可以用在任何其他的違約概率模型中,它可以改進現(xiàn)有的一系列信用風(fēng)險模型。

    此外本文沒有考慮信用風(fēng)險的定價問題,相信本文的內(nèi)生性回收率也可以用于信用風(fēng)險產(chǎn)品(比如,傳統(tǒng)的風(fēng)險債券和信用衍生品CDS,以及結(jié)構(gòu)化產(chǎn)品CDO分券等)的定價中,這也是作者下一步研究的思路和方向。

    [1] Merton R C. On the pricing of corporate debt: The risk structure of interest rates[J]. The Journal of Finance, 1974, 29(2): 449-470.

    [2] Black F, Cox J C. Valuing corporate securities: Some effects of bond indenture provisions[J]. The Journal of Finance, 1976, 31(2): 351-367.

    [3] Jokivuolle E, Peura S. A model for estimating recovery rate and collateral haircuts for bank loans [J]. European Financial Management, 2003, 9(3):611-627

    [4] Giesecke K. Credit risk modeling and valuation: An introduction[J]. Ssrn Electronic Journal, 2004:1-67.

    [5] Asvanunt A, Staal A. The corporate default probability model in Barclays Capital POINT Platform (POINT CDP)[M]. Portfolio Modeling, Barclays Capital, 2009.

    [6] Asvanunt A, Staal A. The POINT Conditional Recovery Rate (CRR) Model[M]. Portfolio Modeling, Barclays Capital, 2009.

    [7] Jarrow R A, Turnbull S M. Pricing derivatives on financial securities subject to credit risk[J]. The Journal of Finance, 1995, 50(1): 53-85.

    [8] Jarrow R A, Lando D, Turnbull S M. A Markov model for the term structure of credit risk spreads[J]. Review of Financial Studies, 1997, 10(2): 481-523.

    [9] Duffie D, Singleton K J. Modeling term structures of defaultable bonds[J]. Review of Financial Studies, 1999, 12(4): 687-720.

    [10] Houweling P, Vorst A C F. An empirical comparison of default swap pricing models[J]. Journal of labor Economics,2015,33(2):269-296.

    [11] Frye J. Depressing recoveries[J]. Risk-London-Risk Magazine Limited, 2000, 13(11): 108-111.

    [12] Pykhtin M. Recovery rates: Unexpected recovery risk[J]. Risk-London-Risk Magazine Limited, 2003, 16(8): 74-79.

    [13] Andersen L, Sidenius J. Extensions to the Gaussian copula: Random recovery and random factor loadings[J]. Journal of Credit Risk Volume, 2004, 1(1): 29-70.

    [14] 黃大海. 違約貸款回收率: 基于國外實證研究的分析[J]. 上海金融, 2006, (10): 55-58.

    [15] 王國棟, 詹原瑞. 信用風(fēng)險中回收率分布的雙 Beta 模型[J]. 中國管理科學(xué), 2011, 19(6): 10-14.

    [16] 汪飛星, 姚磊. 聚合信用風(fēng)險模型的改進和研究[J]. 價值工程, 2013, 32(5): 168-169.

    [17] 陳暮紫, 陳浩, 馬宇超, 等. 基于廣義 Beta 回歸的不良貸款回收率模型[J]. 數(shù)理統(tǒng)計與管理, 2011, 30(5): 810-823.

    [18] Carty L V, Lieberman D. Defaulted bank loan recoveries[J]. Working Paper,Moody’s Investors Service, 1996.

    [19] Frye J. Collateral damage detected[J]. Emerging Issues, 2000 (Sep).

    [20] Bakshi G, Madan D B, Zhang F X. Understanding the role of recovery in default risk models: Empirical comparisons and implied recovery rates[R]. Financial and Economics Board of Governors of the Federal Reserve System (US), 2001.

    [21] Hu Y T, Perraudin W. The dependence of recovery rates and defaults[R]. Discussion Series,Birbeck College and Bank of Working Paper, England, 2002.

    [22] Cantor R, Hamilton D T, Ou S. Default and recovery rates of corporate bond issuers[R].Working Paper, Moodys Investors Services, 2002.

    [23] Carey M, Gordy M. Systematic risk in recoveries on defaulted debt[R].Working Paper Federal Reserve Board,2003.

    [24] Altman E I,Brady B, Restiand A, et al. The link between default and recovery rates: Theory, empirical evidence and implications [J]. Journal of Business,2005,(78):2203-2228.

    [25] Hamilton D T.Default and recovery rates of corporate bond issuers: 2000 [R].Working Paper,Moody.s Investors Service,2001.

    [26] Hu Yen-Ting, Perraudin W.The dependence of recovery rate and defaults[R]. Working Paper,Birbeck College and Bank of England,2002.

    [27] Hull J C, White A D.Valuation of a CDO and an nth to default CDS without monte carlo simulation[J].The Journal of Derivative,2004,12(2): 8-23.

    [28] Das S R, Hanouna P. Implied recovery[J]. Journal of Economic Dynamics and Control, 2009, 33(11): 1837-1857.

    [29] Guo Xin, Jarrow R A, Zeng Yan. Modeling the recovery rate in a reduced form model[J]. Mathematical Finance, 2009, 19(1): 73-97.

    [30] Chava S, Stefanescu C, Turnbull S. Modeling the loss distribution[J]. Management Science, 2011, 57(7): 1267-1287.

    [31] 陸懋祖. 高等時間序列經(jīng)濟計量學(xué)[M]. 上海:上海人民出版社, 1999.

    [32] 孫立娟,風(fēng)險定量分析[M].北京:北京大學(xué)出版社,2011.

    [33] Acerbi C, Tasche D. On the coherence of expected shortfall[J]. Journal of Banking & Finance, 2002, 26(7): 1487-1503.

    Endogenous Recovery Rate and Credit Risk Measurement

    WU Jian-hua1, WANG Xin-jun2, ZHANG Ying1

    (1.School of Mathematical Sciences, University of Jinan, Jinan 250022,China;2.School of Economics, Shandong University, Jinan 250100,China)

    In credit risk models, exogenous recovery rate may neglect the impact on the tail of the loss distribution, and the exogenous specify of the recovery rate leads to the possible model risk. This paper incorporates the factor diffusion process into the structure model of default, derives the inherent relation between the recovery rate and the default probability and analyzes the dependence of expected recovery rate on the expected default probability by using the MC technology. The result shows there are strong negative correlation between expected recovery rates and default probability. Furthermore, the volatility of the asset value has positive compact on the correlation. In the framework of the endogenous recovery rate, the probability distribution of the credit loss is derved, and two index, Credit VaR and ETF, which is the measurement of the credit risk are computed. Finally, the performance of the endogenous recovery rate is tested-based on credit risk model using the market data,which shows that the model can well-character the evolution of the history default probability and recovery rates.

    endogenous recovery rate; factor diffusion process; credit risk measurement; numerical simulation; empirical analysis

    1003-207(2016)01-0001-10

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

    2014-02-19;

    2014-07-17

    教育部人文社科規(guī)劃基金資助項目(13YJAZH091); 國家社會科學(xué)基金資助項目(12BTJ015); 濟南大學(xué)社科基金資助項目(15Y1329);濟南大學(xué)優(yōu)秀人才科研基金資助項目(1008359,1008645)

    簡介:吳建華(1975-),男(漢族),山東博興人,濟南大學(xué)數(shù)學(xué)科學(xué)學(xué)院,講師,校聘A4崗,研究方向:精算與風(fēng)險控制、金融風(fēng)險量化與管理,E-mail:wu88172968@163.com.

    F832

    A

    猜你喜歡
    信用風(fēng)險度量回收率
    有趣的度量
    模糊度量空間的強嵌入
    不同形態(tài)氮肥對棉花15N回收率和產(chǎn)量的影響
    全國農(nóng)膜回收率年底前達到80%以上
    迷向表示分為6個不可約直和的旗流形上不變愛因斯坦度量
    淺析我國商業(yè)銀行信用風(fēng)險管理
    京東商城電子商務(wù)信用風(fēng)險防范策略
    地質(zhì)異常的奇異性度量與隱伏源致礦異常識別
    個人信用風(fēng)險評分的指標(biāo)選擇研究
    信用風(fēng)險的誘發(fā)成因及對策思考
    亚洲一区二区三区色噜噜| 日本与韩国留学比较| 搡女人真爽免费视频火全软件 | 国产单亲对白刺激| 搞女人的毛片| 亚洲国产精品sss在线观看| 国产黄a三级三级三级人| 日本 欧美在线| 久久草成人影院| 精品久久久久久久人妻蜜臀av| 女同久久另类99精品国产91| 草草在线视频免费看| 91麻豆av在线| 一本综合久久免费| 此物有八面人人有两片| 国产国拍精品亚洲av在线观看 | 国产伦一二天堂av在线观看| 最好的美女福利视频网| 超碰av人人做人人爽久久 | 精品国产亚洲在线| 欧美黑人欧美精品刺激| 欧美色视频一区免费| 九色国产91popny在线| 波野结衣二区三区在线 | 99久久精品国产亚洲精品| 美女被艹到高潮喷水动态| 国产极品精品免费视频能看的| 男女下面进入的视频免费午夜| 亚洲av一区综合| 免费在线观看成人毛片| 亚洲中文日韩欧美视频| 日韩人妻高清精品专区| 嫩草影院精品99| 国产麻豆成人av免费视频| 脱女人内裤的视频| 国产高潮美女av| 男女床上黄色一级片免费看| 久久久国产成人免费| 黄色丝袜av网址大全| 手机成人av网站| 久久精品国产自在天天线| 国产黄色小视频在线观看| 两人在一起打扑克的视频| 成人精品一区二区免费| 人人妻人人澡欧美一区二区| 国产黄色小视频在线观看| 亚洲一区高清亚洲精品| www日本在线高清视频| 无遮挡黄片免费观看| 俄罗斯特黄特色一大片| av片东京热男人的天堂| 日本熟妇午夜| 99久久精品热视频| 日本一二三区视频观看| 成人一区二区视频在线观看| 嫩草影视91久久| 成年女人永久免费观看视频| 国产精品美女特级片免费视频播放器| 国产激情欧美一区二区| 日韩精品青青久久久久久| 色综合站精品国产| 热99re8久久精品国产| 99久久精品一区二区三区| 国产高潮美女av| 国产一区二区在线av高清观看| 欧美色视频一区免费| 欧美最黄视频在线播放免费| 日韩精品青青久久久久久| or卡值多少钱| 美女高潮的动态| 欧美乱色亚洲激情| av国产免费在线观看| 狠狠狠狠99中文字幕| 老司机福利观看| www日本在线高清视频| 少妇丰满av| 女人十人毛片免费观看3o分钟| 欧美成人免费av一区二区三区| 国产熟女xx| 亚洲乱码一区二区免费版| 国产高清激情床上av| 国产精品精品国产色婷婷| 亚洲成av人片免费观看| 欧美成人免费av一区二区三区| 婷婷亚洲欧美| 精品国产美女av久久久久小说| 免费看十八禁软件| 99热精品在线国产| 色在线成人网| 久久亚洲真实| 日韩欧美国产在线观看| 一进一出抽搐gif免费好疼| 国产黄a三级三级三级人| 一进一出抽搐动态| 无限看片的www在线观看| 精华霜和精华液先用哪个| 男人舔女人下体高潮全视频| 桃红色精品国产亚洲av| 午夜精品一区二区三区免费看| 国产高清视频在线播放一区| 俄罗斯特黄特色一大片| 久久草成人影院| 手机成人av网站| 国产精品久久久人人做人人爽| 一边摸一边抽搐一进一小说| 美女被艹到高潮喷水动态| 久久欧美精品欧美久久欧美| 久久久久久久精品吃奶| 精品久久久久久久久久免费视频| 91久久精品电影网| 亚洲美女视频黄频| 少妇裸体淫交视频免费看高清| 国产成人欧美在线观看| 99视频精品全部免费 在线| 波多野结衣高清作品| 在线观看免费午夜福利视频| 免费看光身美女| 午夜a级毛片| 18禁裸乳无遮挡免费网站照片| 黑人欧美特级aaaaaa片| 亚洲精品在线观看二区| 免费搜索国产男女视频| 亚洲av成人精品一区久久| 一区二区三区激情视频| 丰满乱子伦码专区| 亚洲avbb在线观看| 国内精品美女久久久久久| 少妇的逼好多水| 国产成人av教育| 国产精品亚洲美女久久久| 中文字幕人成人乱码亚洲影| 禁无遮挡网站| 亚洲一区二区三区不卡视频| 精品国产美女av久久久久小说| 精品欧美国产一区二区三| 天堂av国产一区二区熟女人妻| 国产精品影院久久| 欧美成人一区二区免费高清观看| 精品久久久久久久人妻蜜臀av| 亚洲乱码一区二区免费版| 啪啪无遮挡十八禁网站| 啦啦啦韩国在线观看视频| 神马国产精品三级电影在线观看| 嫩草影院入口| 日日摸夜夜添夜夜添小说| 午夜影院日韩av| 男女视频在线观看网站免费| 啦啦啦韩国在线观看视频| 99国产极品粉嫩在线观看| 女人十人毛片免费观看3o分钟| 国产综合懂色| 精品人妻1区二区| 岛国在线免费视频观看| 淫妇啪啪啪对白视频| 精品人妻偷拍中文字幕| 日韩av在线大香蕉| 免费人成在线观看视频色| 国产精品影院久久| 十八禁人妻一区二区| 可以在线观看毛片的网站| 国产野战对白在线观看| 日本一二三区视频观看| 我要搜黄色片| 国产精品久久电影中文字幕| 中文字幕av成人在线电影| 欧美日韩乱码在线| 淫妇啪啪啪对白视频| 女人高潮潮喷娇喘18禁视频| 桃红色精品国产亚洲av| 一个人看视频在线观看www免费 | 亚洲成av人片在线播放无| 色噜噜av男人的天堂激情| 成年人黄色毛片网站| 欧美日本视频| 一进一出好大好爽视频| 成年免费大片在线观看| 日本黄色视频三级网站网址| 亚洲av第一区精品v没综合| 制服人妻中文乱码| 久久人人精品亚洲av| 嫩草影院入口| 搡老熟女国产l中国老女人| 久久亚洲精品不卡| 亚洲人与动物交配视频| 在线观看免费视频日本深夜| 成人特级av手机在线观看| 成人亚洲精品av一区二区| 999久久久精品免费观看国产| 俄罗斯特黄特色一大片| 国产色婷婷99| 热99在线观看视频| 在线看三级毛片| 欧美日韩精品网址| 嫩草影院精品99| 国产成人啪精品午夜网站| 九九在线视频观看精品| 午夜福利免费观看在线| 亚洲国产精品久久男人天堂| 国产乱人视频| 亚洲欧美日韩卡通动漫| 国产高清三级在线| 国产三级在线视频| 亚洲五月婷婷丁香| 内射极品少妇av片p| 国产精品久久久人人做人人爽| 国产欧美日韩一区二区精品| 婷婷亚洲欧美| 久久精品综合一区二区三区| 国产精品永久免费网站| 12—13女人毛片做爰片一| 偷拍熟女少妇极品色| 欧美极品一区二区三区四区| 91字幕亚洲| 搡老妇女老女人老熟妇| 男女那种视频在线观看| 村上凉子中文字幕在线| 久久精品91无色码中文字幕| 成人18禁在线播放| 亚洲av二区三区四区| 亚洲 欧美 日韩 在线 免费| 日韩av在线大香蕉| 久久久久性生活片| 亚洲在线自拍视频| 香蕉久久夜色| 成人av一区二区三区在线看| 男人舔女人下体高潮全视频| 日本免费a在线| 成人性生交大片免费视频hd| 中出人妻视频一区二区| 美女被艹到高潮喷水动态| 成年版毛片免费区| 亚洲成人中文字幕在线播放| 欧美日韩乱码在线| 国产精品亚洲av一区麻豆| 夜夜躁狠狠躁天天躁| 久久久久久久久中文| 1000部很黄的大片| 日韩欧美一区二区三区在线观看| 午夜福利欧美成人| 淫秽高清视频在线观看| 18禁美女被吸乳视频| 少妇人妻精品综合一区二区 | 每晚都被弄得嗷嗷叫到高潮| 日韩 欧美 亚洲 中文字幕| 黄片小视频在线播放| 91在线精品国自产拍蜜月 | 国产午夜精品久久久久久一区二区三区 | 天堂av国产一区二区熟女人妻| 全区人妻精品视频| 日韩欧美三级三区| 在线观看日韩欧美| 午夜视频国产福利| 亚洲真实伦在线观看| 超碰av人人做人人爽久久 | 变态另类成人亚洲欧美熟女| 亚洲美女视频黄频| 黄片大片在线免费观看| 脱女人内裤的视频| 草草在线视频免费看| 国产视频内射| 亚洲 国产 在线| 高清毛片免费观看视频网站| 婷婷精品国产亚洲av| 黄色视频,在线免费观看| 国内精品久久久久精免费| 国产主播在线观看一区二区| 老司机午夜福利在线观看视频| 91久久精品电影网| tocl精华| 精品国产亚洲在线| 十八禁人妻一区二区| 少妇熟女aⅴ在线视频| 女生性感内裤真人,穿戴方法视频| 黄色丝袜av网址大全| 美女大奶头视频| 可以在线观看的亚洲视频| 又紧又爽又黄一区二区| 色播亚洲综合网| 国产熟女xx| 变态另类丝袜制服| 中文字幕久久专区| 五月玫瑰六月丁香| 久久人妻av系列| 午夜福利欧美成人| 国产v大片淫在线免费观看| 国产乱人视频| 操出白浆在线播放| 日韩欧美 国产精品| 午夜福利在线观看吧| 国产色爽女视频免费观看| 波多野结衣高清无吗| 国产精品爽爽va在线观看网站| 色av中文字幕| 午夜福利18| 久久久国产成人免费| 欧美区成人在线视频| 女人十人毛片免费观看3o分钟| netflix在线观看网站| 亚洲精品美女久久久久99蜜臀| 国产精品久久电影中文字幕| 亚洲欧美日韩无卡精品| 在线看三级毛片| 国产精品综合久久久久久久免费| 成人午夜高清在线视频| 在线播放无遮挡| 变态另类丝袜制服| 日日夜夜操网爽| 日韩欧美三级三区| 日本一本二区三区精品| 亚洲av免费在线观看| 欧美中文综合在线视频| 欧美bdsm另类| 人妻久久中文字幕网| 久久国产乱子伦精品免费另类| 51国产日韩欧美| 国内精品美女久久久久久| 成人无遮挡网站| 国产毛片a区久久久久| av在线天堂中文字幕| 久久久久久久久中文| 床上黄色一级片| 精品一区二区三区av网在线观看| 亚洲,欧美精品.| 国内久久婷婷六月综合欲色啪| 亚洲精品成人久久久久久| 亚洲欧美激情综合另类| 最好的美女福利视频网| e午夜精品久久久久久久| 国内少妇人妻偷人精品xxx网站| 午夜精品在线福利| 国内精品美女久久久久久| 桃红色精品国产亚洲av| 又紧又爽又黄一区二区| 99riav亚洲国产免费| 黄色日韩在线| 久久天躁狠狠躁夜夜2o2o| 在线十欧美十亚洲十日本专区| 香蕉丝袜av| 在线国产一区二区在线| 黄色视频,在线免费观看| 午夜免费男女啪啪视频观看 | 搞女人的毛片| 狠狠狠狠99中文字幕| 内射极品少妇av片p| 俄罗斯特黄特色一大片| 熟妇人妻久久中文字幕3abv| 搡老妇女老女人老熟妇| 色老头精品视频在线观看| 少妇裸体淫交视频免费看高清| 91在线观看av| 日韩欧美免费精品| 亚洲欧美日韩东京热| 51午夜福利影视在线观看| 欧美绝顶高潮抽搐喷水| 日本a在线网址| 国产免费男女视频| 亚洲欧美日韩卡通动漫| 老司机在亚洲福利影院| 免费看十八禁软件| 亚洲电影在线观看av| 日本一二三区视频观看| 日韩中文字幕欧美一区二区| 日本免费a在线| 啦啦啦观看免费观看视频高清| 人人妻人人澡欧美一区二区| 国产成人影院久久av| 国产探花在线观看一区二区| 又黄又粗又硬又大视频| 亚洲成av人片在线播放无| eeuss影院久久| 亚洲精华国产精华精| 中文字幕久久专区| 精品国产三级普通话版| 最近在线观看免费完整版| 久久精品国产亚洲av涩爱 | 欧美最黄视频在线播放免费| 18禁黄网站禁片免费观看直播| 午夜福利视频1000在线观看| a级一级毛片免费在线观看| 久久精品国产清高在天天线| 国产日本99.免费观看| 国产亚洲精品久久久久久毛片| a在线观看视频网站| 国产欧美日韩精品一区二区| 在线视频色国产色| 男人和女人高潮做爰伦理| 成人国产综合亚洲| 久久精品综合一区二区三区| 高清在线国产一区| 无遮挡黄片免费观看| tocl精华| av天堂中文字幕网| 国产主播在线观看一区二区| 美女cb高潮喷水在线观看| 精品福利观看| 99国产综合亚洲精品| 欧美一级a爱片免费观看看| 3wmmmm亚洲av在线观看| 女警被强在线播放| 麻豆成人av在线观看| 国产精品女同一区二区软件 | 日韩欧美免费精品| 欧美最黄视频在线播放免费| 国产蜜桃级精品一区二区三区| 国产69精品久久久久777片| av在线蜜桃| 欧美黑人巨大hd| 国产伦一二天堂av在线观看| 老汉色av国产亚洲站长工具| 国产精品永久免费网站| 久久国产乱子伦精品免费另类| 国产麻豆成人av免费视频| 欧美又色又爽又黄视频| 免费看日本二区| 男女下面进入的视频免费午夜| 日本三级黄在线观看| 午夜a级毛片| 成人av在线播放网站| 日韩国内少妇激情av| 欧美性猛交黑人性爽| 免费av观看视频| 超碰av人人做人人爽久久 | 国产精品爽爽va在线观看网站| 好看av亚洲va欧美ⅴa在| 69人妻影院| 国产精品98久久久久久宅男小说| 欧美乱妇无乱码| 日韩欧美 国产精品| 非洲黑人性xxxx精品又粗又长| 国产精品影院久久| 色av中文字幕| 在线免费观看的www视频| 欧美最黄视频在线播放免费| 叶爱在线成人免费视频播放| 网址你懂的国产日韩在线| 久久久久九九精品影院| 日韩人妻高清精品专区| 成年女人毛片免费观看观看9| 日本 av在线| 九色国产91popny在线| 国产成人福利小说| 欧美+日韩+精品| 日韩av在线大香蕉| 欧美三级亚洲精品| 亚洲av中文字字幕乱码综合| 一进一出抽搐gif免费好疼| 夜夜爽天天搞| 特级一级黄色大片| a级毛片a级免费在线| 国产综合懂色| 在线观看一区二区三区| 在线看三级毛片| 中文字幕av成人在线电影| 久久午夜亚洲精品久久| 小说图片视频综合网站| 99久久无色码亚洲精品果冻| 特级一级黄色大片| 免费看美女性在线毛片视频| 夜夜爽天天搞| 又紧又爽又黄一区二区| 国产精品女同一区二区软件 | 国产精品久久久久久人妻精品电影| 内射极品少妇av片p| 51午夜福利影视在线观看| 手机成人av网站| 精品久久久久久久末码| 两性午夜刺激爽爽歪歪视频在线观看| 日本免费a在线| 国产精品乱码一区二三区的特点| 亚洲av熟女| 免费无遮挡裸体视频| 天天躁日日操中文字幕| 婷婷精品国产亚洲av| 免费大片18禁| 亚洲va日本ⅴa欧美va伊人久久| 制服丝袜大香蕉在线| 欧美日韩国产亚洲二区| 舔av片在线| 久久精品亚洲精品国产色婷小说| 香蕉丝袜av| 欧美乱色亚洲激情| 性色av乱码一区二区三区2| 久久精品国产综合久久久| 琪琪午夜伦伦电影理论片6080| a级一级毛片免费在线观看| 国产高清视频在线观看网站| 午夜精品一区二区三区免费看| 国产成年人精品一区二区| 免费在线观看日本一区| 床上黄色一级片| 很黄的视频免费| av中文乱码字幕在线| 97超级碰碰碰精品色视频在线观看| 97碰自拍视频| 99久久精品国产亚洲精品| 国产三级中文精品| 国产美女午夜福利| 国产乱人视频| 99热6这里只有精品| 亚洲精品亚洲一区二区| 日韩有码中文字幕| 国产高清视频在线播放一区| 久99久视频精品免费| 亚洲成人久久性| 给我免费播放毛片高清在线观看| 小说图片视频综合网站| 最近最新免费中文字幕在线| 少妇高潮的动态图| 69av精品久久久久久| 欧美激情久久久久久爽电影| 欧美日韩乱码在线| 欧美乱色亚洲激情| 超碰av人人做人人爽久久 | 2021天堂中文幕一二区在线观| 精品无人区乱码1区二区| 日韩精品青青久久久久久| 国产中年淑女户外野战色| 亚洲av成人不卡在线观看播放网| 久99久视频精品免费| 日韩免费av在线播放| bbb黄色大片| 好男人电影高清在线观看| 在线观看免费午夜福利视频| 日韩免费av在线播放| 国产v大片淫在线免费观看| 欧美黑人欧美精品刺激| 又黄又粗又硬又大视频| 欧美成狂野欧美在线观看| 亚洲自拍偷在线| 国产高清三级在线| 午夜亚洲福利在线播放| 国内精品美女久久久久久| 亚洲精品456在线播放app | 国产精品香港三级国产av潘金莲| 日韩大尺度精品在线看网址| 成人午夜高清在线视频| 国产高潮美女av| 欧美日韩亚洲国产一区二区在线观看| 床上黄色一级片| 少妇高潮的动态图| 亚洲人与动物交配视频| 国产伦精品一区二区三区视频9 | 一二三四社区在线视频社区8| 国产伦一二天堂av在线观看| 成年人黄色毛片网站| 久久久久亚洲av毛片大全| 国产精品久久久久久久久免 | 观看美女的网站| 999久久久精品免费观看国产| 老司机午夜福利在线观看视频| 国产高清视频在线观看网站| 757午夜福利合集在线观看| 欧美色视频一区免费| 日韩 欧美 亚洲 中文字幕| 啦啦啦韩国在线观看视频| 69av精品久久久久久| 亚洲专区中文字幕在线| 日本在线视频免费播放| 听说在线观看完整版免费高清| 全区人妻精品视频| 成年女人毛片免费观看观看9| 亚洲精品美女久久久久99蜜臀| 国产精品乱码一区二三区的特点| 国产真实乱freesex| 国产高清激情床上av| 午夜日韩欧美国产| 嫁个100分男人电影在线观看| 国产午夜精品久久久久久一区二区三区 | 看片在线看免费视频| 国产精品久久久久久人妻精品电影| 好男人电影高清在线观看| 精品久久久久久久久久久久久| 18禁在线播放成人免费| 久久精品91蜜桃| 特大巨黑吊av在线直播| 亚洲国产欧美人成| 欧美日韩综合久久久久久 | 日日干狠狠操夜夜爽| 两个人的视频大全免费| 欧美成人免费av一区二区三区| 亚洲国产欧美人成| 欧美日韩综合久久久久久 | 久久6这里有精品| 精品人妻偷拍中文字幕| 天堂影院成人在线观看| 国产伦人伦偷精品视频| 国产高清视频在线播放一区| 国产精品永久免费网站| 我要搜黄色片| 国产高潮美女av| 丰满人妻熟妇乱又伦精品不卡| 听说在线观看完整版免费高清| 不卡一级毛片| 亚洲不卡免费看| 叶爱在线成人免费视频播放| 网址你懂的国产日韩在线| 国产亚洲精品久久久com| 大型黄色视频在线免费观看| 国产日本99.免费观看| 国产一区在线观看成人免费| 内地一区二区视频在线| 国产精品,欧美在线| 中文字幕高清在线视频| 最近最新中文字幕大全免费视频| 全区人妻精品视频| 黄色片一级片一级黄色片| 成人三级黄色视频| 中文在线观看免费www的网站| 变态另类丝袜制服| 搞女人的毛片| 高清毛片免费观看视频网站| 大型黄色视频在线免费观看| 精品日产1卡2卡| 老汉色av国产亚洲站长工具| 成人高潮视频无遮挡免费网站|