• <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ā)成因及對策思考
    午夜福利18| 最近最新中文字幕大全免费视频| 一级作爱视频免费观看| 欧美久久黑人一区二区| 亚洲 欧美 日韩 在线 免费| 国产免费av片在线观看野外av| 99国产综合亚洲精品| 日韩有码中文字幕| 可以在线观看毛片的网站| 亚洲一区中文字幕在线| 久久亚洲精品不卡| 色播亚洲综合网| 久久久国产成人精品二区| 一级毛片高清免费大全| 久久欧美精品欧美久久欧美| 亚洲精品美女久久久久99蜜臀| 丝袜人妻中文字幕| 国产又黄又爽又无遮挡在线| 精品国产亚洲在线| 亚洲精品色激情综合| 中亚洲国语对白在线视频| 女人被狂操c到高潮| 色综合婷婷激情| 悠悠久久av| 精品乱码久久久久久99久播| 禁无遮挡网站| 精品国产乱子伦一区二区三区| 精品不卡国产一区二区三区| 深夜精品福利| 麻豆成人午夜福利视频| 久久 成人 亚洲| 老汉色∧v一级毛片| 午夜免费观看网址| 香蕉国产在线看| 制服人妻中文乱码| 欧美亚洲日本最大视频资源| 天堂√8在线中文| 成人手机av| 国产又黄又爽又无遮挡在线| 丝袜在线中文字幕| 女生性感内裤真人,穿戴方法视频| 欧美成狂野欧美在线观看| 中文在线观看免费www的网站 | 亚洲av成人不卡在线观看播放网| 免费高清在线观看日韩| 国内精品久久久久久久电影| 日日摸夜夜添夜夜添小说| 国产一级毛片七仙女欲春2 | 在线国产一区二区在线| 国产亚洲精品第一综合不卡| 免费女性裸体啪啪无遮挡网站| av片东京热男人的天堂| 日韩国内少妇激情av| 村上凉子中文字幕在线| 欧美日韩亚洲综合一区二区三区_| 很黄的视频免费| 久久久久久久久中文| 天天躁夜夜躁狠狠躁躁| 亚洲av日韩精品久久久久久密| 欧美乱妇无乱码| 成人亚洲精品av一区二区| 中文字幕av电影在线播放| 美国免费a级毛片| 国产成人欧美| av欧美777| 国产精品日韩av在线免费观看| 免费女性裸体啪啪无遮挡网站| 好男人在线观看高清免费视频 | 午夜福利欧美成人| 波多野结衣高清作品| 欧美性猛交╳xxx乱大交人| 桃红色精品国产亚洲av| 99热6这里只有精品| 中文字幕最新亚洲高清| 人人妻人人澡欧美一区二区| 国产一区在线观看成人免费| 日韩欧美三级三区| 51午夜福利影视在线观看| 亚洲精品久久国产高清桃花| 热99re8久久精品国产| 亚洲av中文字字幕乱码综合 | 狂野欧美激情性xxxx| 黄色a级毛片大全视频| 国产亚洲精品av在线| 日本成人三级电影网站| 亚洲人成网站在线播放欧美日韩| 嫩草影视91久久| 亚洲真实伦在线观看| 男人舔女人下体高潮全视频| 亚洲精品美女久久久久99蜜臀| 欧美激情极品国产一区二区三区| 满18在线观看网站| 99在线视频只有这里精品首页| 久久草成人影院| 国产精品一区二区免费欧美| 一本大道久久a久久精品| 日本成人三级电影网站| 村上凉子中文字幕在线| 99热只有精品国产| 极品教师在线免费播放| 给我免费播放毛片高清在线观看| 中文字幕久久专区| 大型黄色视频在线免费观看| 一夜夜www| 精品日产1卡2卡| 欧美日韩精品网址| 99热只有精品国产| 国产成人av教育| 99久久99久久久精品蜜桃| 亚洲美女黄片视频| 嫁个100分男人电影在线观看| 高潮久久久久久久久久久不卡| 丝袜在线中文字幕| 搡老妇女老女人老熟妇| 在线av久久热| 精品高清国产在线一区| 国产又色又爽无遮挡免费看| 亚洲性夜色夜夜综合| 一进一出抽搐动态| 成年版毛片免费区| 国产欧美日韩精品亚洲av| 在线观看免费午夜福利视频| 亚洲一区中文字幕在线| 亚洲av美国av| 国产97色在线日韩免费| 欧美日韩精品网址| 亚洲精品美女久久av网站| 中文字幕精品亚洲无线码一区 | 亚洲熟妇中文字幕五十中出| 女人被狂操c到高潮| 国产精品98久久久久久宅男小说| 国产精品九九99| 国产日本99.免费观看| 亚洲va日本ⅴa欧美va伊人久久| 久久亚洲精品不卡| 无人区码免费观看不卡| 变态另类成人亚洲欧美熟女| 久久久久久国产a免费观看| 精品无人区乱码1区二区| 99国产精品一区二区三区| 黄色女人牲交| 精品无人区乱码1区二区| 少妇熟女aⅴ在线视频| 久久亚洲精品不卡| 欧美人与性动交α欧美精品济南到| 亚洲av片天天在线观看| 香蕉av资源在线| 久久久国产成人免费| 美国免费a级毛片| 久久精品aⅴ一区二区三区四区| 在线观看免费日韩欧美大片| 丁香六月欧美| 国产成人系列免费观看| 国产av又大| 亚洲国产日韩欧美精品在线观看 | 人妻丰满熟妇av一区二区三区| 97碰自拍视频| 大型av网站在线播放| 男女那种视频在线观看| 婷婷精品国产亚洲av在线| 最近最新中文字幕大全免费视频| 丁香六月欧美| 他把我摸到了高潮在线观看| 国产日本99.免费观看| 老司机在亚洲福利影院| 天天添夜夜摸| 我的亚洲天堂| 亚洲精品久久成人aⅴ小说| 久久九九热精品免费| 国产亚洲精品综合一区在线观看 | a级毛片在线看网站| xxx96com| 一级a爱片免费观看的视频| 俺也久久电影网| 精品高清国产在线一区| 精品乱码久久久久久99久播| 特大巨黑吊av在线直播 | 国产视频内射| 在线观看www视频免费| 久久久久久大精品| 日韩欧美国产在线观看| 久久久国产欧美日韩av| 免费看十八禁软件| 久久人妻av系列| 亚洲电影在线观看av| 久久久久久国产a免费观看| 99久久综合精品五月天人人| 亚洲精品国产区一区二| 一区二区日韩欧美中文字幕| 日韩大码丰满熟妇| 一级毛片精品| 午夜免费成人在线视频| 中国美女看黄片| 中文字幕av电影在线播放| 99re在线观看精品视频| 女警被强在线播放| 亚洲电影在线观看av| av在线天堂中文字幕| 国产精品永久免费网站| 亚洲欧美日韩无卡精品| 女警被强在线播放| 久久精品91蜜桃| 嫩草影院精品99| 亚洲国产精品久久男人天堂| 亚洲真实伦在线观看| www.www免费av| 色综合婷婷激情| 久久99热这里只有精品18| 欧美日韩福利视频一区二区| 国产精品国产高清国产av| 91成人精品电影| 最近最新中文字幕大全免费视频| 国产欧美日韩一区二区精品| 久久99热这里只有精品18| 色综合站精品国产| 久久精品成人免费网站| 在线观看舔阴道视频| 午夜福利一区二区在线看| 亚洲成人久久爱视频| 午夜激情av网站| 色综合站精品国产| 夜夜躁狠狠躁天天躁| 亚洲熟女毛片儿| 不卡一级毛片| 婷婷丁香在线五月| 我的亚洲天堂| 国产精品亚洲一级av第二区| 久久国产精品男人的天堂亚洲| 19禁男女啪啪无遮挡网站| 91国产中文字幕| 99久久综合精品五月天人人| 亚洲午夜精品一区,二区,三区| 久久人人精品亚洲av| 久久久久久大精品| 女警被强在线播放| 色播在线永久视频| 欧美一区二区精品小视频在线| 岛国视频午夜一区免费看| 午夜免费激情av| 亚洲欧美日韩高清在线视频| 男人舔奶头视频| 最近最新免费中文字幕在线| 色婷婷久久久亚洲欧美| 亚洲性夜色夜夜综合| 亚洲人成网站高清观看| 一区二区日韩欧美中文字幕| 欧美色欧美亚洲另类二区| 精品国产亚洲在线| 这个男人来自地球电影免费观看| 曰老女人黄片| 91av网站免费观看| 精品高清国产在线一区| 哪里可以看免费的av片| 观看免费一级毛片| 男女做爰动态图高潮gif福利片| 免费在线观看成人毛片| 亚洲 国产 在线| 国产成人影院久久av| 香蕉久久夜色| 一个人免费在线观看的高清视频| 啦啦啦观看免费观看视频高清| 日韩 欧美 亚洲 中文字幕| bbb黄色大片| 精品不卡国产一区二区三区| 精品免费久久久久久久清纯| 女人爽到高潮嗷嗷叫在线视频| 88av欧美| 观看免费一级毛片| 色综合亚洲欧美另类图片| 最新在线观看一区二区三区| 好男人电影高清在线观看| 日韩大码丰满熟妇| 亚洲精品一卡2卡三卡4卡5卡| 国产精品av久久久久免费| 久久午夜综合久久蜜桃| 夜夜看夜夜爽夜夜摸| 亚洲国产毛片av蜜桃av| 波多野结衣av一区二区av| 叶爱在线成人免费视频播放| 99久久精品国产亚洲精品| 高潮久久久久久久久久久不卡| 亚洲国产高清在线一区二区三 | 亚洲五月婷婷丁香| 国产精品1区2区在线观看.| 久久中文字幕一级| 亚洲午夜理论影院| 岛国视频午夜一区免费看| 亚洲男人的天堂狠狠| 午夜免费成人在线视频| x7x7x7水蜜桃| 亚洲专区中文字幕在线| 久9热在线精品视频| 一级a爱视频在线免费观看| 18禁美女被吸乳视频| 波多野结衣巨乳人妻| 日本精品一区二区三区蜜桃| 91大片在线观看| 午夜a级毛片| 久久久国产成人精品二区| 嫩草影视91久久| 真人一进一出gif抽搐免费| 777久久人妻少妇嫩草av网站| 中文字幕另类日韩欧美亚洲嫩草| 亚洲午夜理论影院| 琪琪午夜伦伦电影理论片6080| 久久午夜综合久久蜜桃| 国产91精品成人一区二区三区| 老司机午夜福利在线观看视频| 欧美又色又爽又黄视频| 黄色片一级片一级黄色片| 精品国产超薄肉色丝袜足j| 精品免费久久久久久久清纯| 中文字幕精品免费在线观看视频| 久久中文字幕人妻熟女| 国产亚洲精品一区二区www| 在线天堂中文资源库| 久久午夜综合久久蜜桃| 久久中文看片网| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产精品999在线| 免费在线观看亚洲国产| 国产男靠女视频免费网站| 在线观看日韩欧美| 每晚都被弄得嗷嗷叫到高潮| 精品欧美国产一区二区三| 亚洲熟女毛片儿| 国产精品 国内视频| 757午夜福利合集在线观看| 免费高清在线观看日韩| 韩国av一区二区三区四区| 女生性感内裤真人,穿戴方法视频| 亚洲一区高清亚洲精品| 身体一侧抽搐| 国产欧美日韩一区二区精品| 亚洲成人久久性| 一本一本综合久久| 每晚都被弄得嗷嗷叫到高潮| 操出白浆在线播放| 天天躁夜夜躁狠狠躁躁| 99热6这里只有精品| 91字幕亚洲| 久9热在线精品视频| 999久久久精品免费观看国产| 在线观看午夜福利视频| 亚洲 欧美 日韩 在线 免费| 好男人在线观看高清免费视频 | xxxwww97欧美| 人人妻人人澡欧美一区二区| 国产成人欧美在线观看| 国产真人三级小视频在线观看| 成年女人毛片免费观看观看9| 亚洲美女黄片视频| 欧美一级毛片孕妇| 国产高清激情床上av| 人妻久久中文字幕网| 97人妻精品一区二区三区麻豆 | 伊人久久大香线蕉亚洲五| 中文字幕人妻熟女乱码| 国产精品综合久久久久久久免费| 一进一出抽搐动态| av中文乱码字幕在线| aaaaa片日本免费| 人人妻,人人澡人人爽秒播| 精品一区二区三区四区五区乱码| av中文乱码字幕在线| 桃色一区二区三区在线观看| 美女免费视频网站| 国产久久久一区二区三区| 一a级毛片在线观看| 91麻豆av在线| 午夜a级毛片| 欧洲精品卡2卡3卡4卡5卡区| 免费电影在线观看免费观看| 日韩大尺度精品在线看网址| 色综合站精品国产| 国产伦在线观看视频一区| 91av网站免费观看| 亚洲成人久久爱视频| 91字幕亚洲| 中文字幕精品亚洲无线码一区 | 欧美成狂野欧美在线观看| 亚洲三区欧美一区| 一本一本综合久久| 伊人久久大香线蕉亚洲五| 久久久精品国产亚洲av高清涩受| 国产精品影院久久| 男人的好看免费观看在线视频 | 一本大道久久a久久精品| 亚洲熟妇熟女久久| 亚洲精品国产一区二区精华液| 97碰自拍视频| 亚洲五月天丁香| 丁香六月欧美| av天堂在线播放| cao死你这个sao货| 亚洲中文日韩欧美视频| 久久99热这里只有精品18| 久久久国产成人免费| 成人永久免费在线观看视频| 宅男免费午夜| 少妇的丰满在线观看| 一进一出抽搐gif免费好疼| 久久久久久久精品吃奶| 一级毛片高清免费大全| 久久精品国产亚洲av香蕉五月| 中文字幕精品亚洲无线码一区 | 91麻豆精品激情在线观看国产| 女人被狂操c到高潮| 国产精品99久久99久久久不卡| 欧美一级a爱片免费观看看 | 女同久久另类99精品国产91| 亚洲无线在线观看| 美女扒开内裤让男人捅视频| 动漫黄色视频在线观看| 国产成人精品无人区| 一级毛片女人18水好多| 一a级毛片在线观看| 欧美成狂野欧美在线观看| 国产精品98久久久久久宅男小说| 日韩欧美免费精品| 麻豆一二三区av精品| 亚洲av成人av| 国内久久婷婷六月综合欲色啪| 欧美+亚洲+日韩+国产| 亚洲五月婷婷丁香| 亚洲午夜精品一区,二区,三区| 中文字幕人妻丝袜一区二区| 日本撒尿小便嘘嘘汇集6| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品粉嫩美女一区| 国内久久婷婷六月综合欲色啪| ponron亚洲| 国产成人影院久久av| 757午夜福利合集在线观看| 亚洲av成人av| 黄色丝袜av网址大全| 97人妻精品一区二区三区麻豆 | 亚洲一区二区三区不卡视频| 久久中文看片网| 国产日本99.免费观看| 精品一区二区三区视频在线观看免费| 亚洲一码二码三码区别大吗| 看片在线看免费视频| 国产亚洲欧美在线一区二区| 香蕉久久夜色| 中国美女看黄片| 亚洲国产欧美日韩在线播放| 国产日本99.免费观看| 国产一区在线观看成人免费| 黑人欧美特级aaaaaa片| 成人亚洲精品一区在线观看| 日韩成人在线观看一区二区三区| 精品欧美国产一区二区三| 夜夜看夜夜爽夜夜摸| www.www免费av| 侵犯人妻中文字幕一二三四区| 免费人成视频x8x8入口观看| 日本熟妇午夜| 少妇裸体淫交视频免费看高清 | 国产精品久久电影中文字幕| 亚洲国产毛片av蜜桃av| 两人在一起打扑克的视频| 久久久久久九九精品二区国产 | 两人在一起打扑克的视频| 自线自在国产av| 亚洲成人久久爱视频| 一进一出抽搐gif免费好疼| 香蕉国产在线看| 精品一区二区三区视频在线观看免费| 久久精品成人免费网站| 国语自产精品视频在线第100页| av视频在线观看入口| 欧美绝顶高潮抽搐喷水| 听说在线观看完整版免费高清| 久久久久久久久中文| 久久久久久久午夜电影| 亚洲欧美日韩高清在线视频| 国产成年人精品一区二区| 欧美一级a爱片免费观看看 | 狂野欧美激情性xxxx| 久久天堂一区二区三区四区| 99国产精品99久久久久| 免费人成视频x8x8入口观看| 免费在线观看影片大全网站| 可以免费在线观看a视频的电影网站| avwww免费| 一级毛片女人18水好多| 精品免费久久久久久久清纯| 欧美三级亚洲精品| 国产97色在线日韩免费| 18禁国产床啪视频网站| 亚洲中文日韩欧美视频| 母亲3免费完整高清在线观看| 身体一侧抽搐| 久久草成人影院| 99热只有精品国产| 久久久久免费精品人妻一区二区 | 欧美性长视频在线观看| 国产精华一区二区三区| 午夜福利18| 91大片在线观看| 国产一区二区三区在线臀色熟女| 亚洲色图av天堂| 别揉我奶头~嗯~啊~动态视频| 国产精品久久久av美女十八| 成人18禁高潮啪啪吃奶动态图| 国产一区二区三区在线臀色熟女| 亚洲专区国产一区二区| 久久久精品国产亚洲av高清涩受| 男女做爰动态图高潮gif福利片| 成人一区二区视频在线观看| 亚洲精品中文字幕一二三四区| 在线观看66精品国产| bbb黄色大片| 国内揄拍国产精品人妻在线 | 黄片小视频在线播放| 黑人操中国人逼视频| 满18在线观看网站| 每晚都被弄得嗷嗷叫到高潮| 成人三级做爰电影| 黄色视频不卡| 亚洲色图av天堂| 日韩有码中文字幕| 国产高清激情床上av| 国产1区2区3区精品| 成人午夜高清在线视频 | 久久久国产成人精品二区| 日韩欧美三级三区| 欧美绝顶高潮抽搐喷水| 精品欧美国产一区二区三| 美女大奶头视频| 久久久久久久久免费视频了| 搡老妇女老女人老熟妇| 老汉色∧v一级毛片| 欧美一级毛片孕妇| 日韩三级视频一区二区三区| 成年版毛片免费区| 午夜老司机福利片| 久久精品国产清高在天天线| 国产亚洲精品一区二区www| 男人的好看免费观看在线视频 | 亚洲欧美一区二区三区黑人| 悠悠久久av| 亚洲精品国产一区二区精华液| 女生性感内裤真人,穿戴方法视频| 白带黄色成豆腐渣| 日本一区二区免费在线视频| 国内揄拍国产精品人妻在线 | 亚洲久久久国产精品| 99国产精品一区二区三区| 99在线视频只有这里精品首页| 国产又色又爽无遮挡免费看| 中文在线观看免费www的网站 | 免费在线观看黄色视频的| 久久亚洲精品不卡| 精品卡一卡二卡四卡免费| 18禁裸乳无遮挡免费网站照片 | 人人妻人人澡人人看| 国产高清激情床上av| 亚洲国产欧美网| 桃红色精品国产亚洲av| 亚洲国产精品999在线| 脱女人内裤的视频| 日本 欧美在线| 国产男靠女视频免费网站| 黄片播放在线免费| 51午夜福利影视在线观看| 一进一出抽搐gif免费好疼| 午夜成年电影在线免费观看| 国产一区二区三区在线臀色熟女| 亚洲中文字幕日韩| 91麻豆av在线| 天堂动漫精品| 亚洲色图av天堂| 欧美色欧美亚洲另类二区| 亚洲精品色激情综合| 岛国在线观看网站| 非洲黑人性xxxx精品又粗又长| 亚洲成人国产一区在线观看| 欧美乱色亚洲激情| 在线观看免费午夜福利视频| 一进一出好大好爽视频| 老司机深夜福利视频在线观看| 在线观看日韩欧美| 国产欧美日韩一区二区三| 中文字幕精品免费在线观看视频| 一级作爱视频免费观看| 久久99热这里只有精品18| 亚洲一区二区三区不卡视频| 久久香蕉激情| 一区二区日韩欧美中文字幕| 女人爽到高潮嗷嗷叫在线视频| 好男人在线观看高清免费视频 | 国产激情偷乱视频一区二区| 午夜福利在线观看吧| 久久国产乱子伦精品免费另类| 99精品在免费线老司机午夜| 夜夜爽天天搞| 人人妻人人澡欧美一区二区| 免费观看精品视频网站| 国产精品久久视频播放| 欧美av亚洲av综合av国产av| 免费一级毛片在线播放高清视频| 视频在线观看一区二区三区| 午夜福利18| 黄色毛片三级朝国网站| 一本久久中文字幕| 亚洲国产看品久久| 欧美日本视频| 别揉我奶头~嗯~啊~动态视频| 一进一出好大好爽视频| 在线十欧美十亚洲十日本专区| 国产成人av教育| 男女床上黄色一级片免费看| 欧美不卡视频在线免费观看 |