亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        證據(jù)權重方法在企業(yè)信用風險評估應用

        2020-02-22 06:52:25危明鑄沈鳳山袁峰麥偉杰
        科技創(chuàng)新導報 2020年29期
        關鍵詞:企業(yè)信用信用風險

        危明鑄 沈鳳山 袁峰 麥偉杰

        摘? 要:本文以廣東省各個行政部門的企業(yè)數(shù)據(jù)為基礎,研究了證據(jù)權重在企業(yè)信用風險等級評估中的應用,并且根據(jù)國家“企業(yè)信用風險分類評價標準”系統(tǒng)建立了企業(yè)信用評價指標,成功地將證據(jù)權重邏輯回歸算法應用到真實的企業(yè)信用數(shù)據(jù)集,建立企業(yè)信用風險評估模型,使得監(jiān)管部門能夠準確地掌握企業(yè)的信用情況。通過與經(jīng)典的邏輯算法比較,驗證了該方法的有效性。

        關鍵詞:證據(jù)權重? 邏輯回歸? 信用風險? 企業(yè)信用

        中圖分類號:O212? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻標識碼:A? ? ? ? ? ? ? ? ? ? 文章編號:1674-098X(2020)10(b)-0137-07

        Abstract: Based on the enterprise data of various administrative departments in Guangdong Province, this paper studies the application of weight of evidence in the assessment of enterprise credit risk levels, and establishes enterprise credit evaluation indicators based on the national "Enterprise Credit Risk Classification and Evaluation Standards". The weight-of-evidence logistic regression algorithm is applied to real enterprise credit data sets to establish enterprise credit risk assessment model, enabling regulators to accurately grasp the enterprise credit situation. Compared with the classical logic algorithm, the validity of the method is verified.

        Key Words: Weight of evidence; Logistic regression; Credit risk; Credit of enterprise

        企業(yè)信用是市場經(jīng)濟的產(chǎn)物,是對各類市場參與主體履行相應經(jīng)濟契約的能力及其企業(yè)整體的可信程度所進行的一種綜合分析和測定,是企業(yè)的一項重要無形資產(chǎn)。隨著市場主體“寬進嚴管”改革的不斷深入,國務院相繼頒發(fā)了《國務院辦公廳關于推廣隨機抽查規(guī)范事中事后監(jiān)管的通知》(國辦發(fā)[2015]58號)[1]、《國務院關于“先照后證”改革后加強事中事后監(jiān)管的意見》(國發(fā)[2015]62號)[2]和《國務院關于印發(fā)2016年推進簡政放權放管結合優(yōu)化服務改革工作要點的通知》(國發(fā)[2016]30號)[3]等相關文件,企業(yè)信用在社會經(jīng)營活動中變成一個有效的“身份證”,塑造了一個企業(yè)在社會的面目和形象。

        針對目前市場主體規(guī)模龐大、難以把握重點監(jiān)管對象、雙隨機抽查的靶向性不強等不足,有些學者已經(jīng)發(fā)表了自己的研究成果。Odom等人[4]早在1990年把人工神經(jīng)網(wǎng)絡應用到企業(yè)信用風險評估上,并將其與經(jīng)典的多元回歸分析比較,實驗表明人工神經(jīng)網(wǎng)絡具有更優(yōu)越的性能;Prinzie等人把邏輯回歸(Logistic)引入隨機森林算法并對其進行優(yōu)化以及改進,然后將改進的算法應用于預測公司的信用風險[5];Lin F等人研究隨機森林與KMV模型結合,提出將違約距離作為隨機森林的輸入,實驗表明對企業(yè)信用風險預測性能更有效[6];Traskin等人利用隨機森林具有篩選重要變量的特征,提出將其應用在保險公司償付判別中[7];吳麗麗運用Logistic回歸模型深入探討了我國商業(yè)銀行信用風險監(jiān)管的問題[8];郭玉華根據(jù)微型企業(yè)的特征,運用Logit模型進行實證分析,銀行可以借助該模型對微企的信用風險進行評估[9];方匡南,范新妍等人指出傳統(tǒng)的Logistic回歸建立企業(yè)運行風險預警模型效果不夠好,提出了基于網(wǎng)絡結構關系的Logistic模型[10];楊俊等人使用Gradient Boosting算法對中國建設銀行上海分行的企業(yè)貸款客戶數(shù)據(jù)建立模型,并和邏輯回歸以及專家規(guī)則進行橫向比較,結果表明Gradient Boosting算法的模型要優(yōu)于另外兩種模型[11];熊正得等人利用因子分析法對深滬A股上市的制造企業(yè)財務數(shù)據(jù)構建風險評價體系,并在違約測度階段應用Logistic回歸對不同組樣本進行測度[12];劉丹等人使用證據(jù)權重、逐步回歸對信用評價指標進行篩選,構建一套具有區(qū)分違約能力的信用風險模型[13];劉麗君、韓靜磊等人運用WOE法評估了生活垃圾焚燒廠固化飛灰中重金屬的非致癌健康風險,并將其與傳統(tǒng)的非致癌健康風險評價方法進行比較[14];趙雅迪等人通過信息值(information value,Ⅳ)及證據(jù)權重轉化(weight of evidence,WOE)結合邏輯回歸算法構建用電客戶電費風險預測模型[15];陳超等人采用卡方分箱法和WOE編碼判別確定影響轉爐理想終點目標的關鍵工藝參數(shù),并且運用邏輯回歸算法對編碼后數(shù)據(jù)構建的轉爐操作工藝評價模型[16];危明鑄等人在企業(yè)信用風險預測上綜合運用了各種機器學習算法比較各種方法的優(yōu)缺點[17]。

        本文對過去一定時期(如1年)出現(xiàn)過信用風險事件(如偷稅、行政處罰等)的企業(yè)數(shù)據(jù)集,根據(jù)國家“企業(yè)信用風險分類評價標準”系統(tǒng)建立企業(yè)信用評價指標體系。應用機器學習方法預處理評價指標數(shù)據(jù)集,如特征選擇(information value,IV)、缺失值處理、異常值處理,并將評價指標以證據(jù)權重(Weight of Evidence,WOE)編碼方式離散化形成woe數(shù)據(jù)集,之后結合邏輯回歸算法學習企業(yè)過去發(fā)生信用風險事件與否的規(guī)律,建模企業(yè)信用風險評分模型。

        1? 相關技術

        1.1 邏輯回歸

        邏輯回歸(Logistic regression)是一種比較流行的二分類的機器學習算法。例如,探討企業(yè)信用風險與那些變量有關,可將企業(yè)數(shù)據(jù)集標記為兩組,假設存在信用風險的一組企業(yè)標記為1,那么沒有信用風險的一組為0,并用Y表示響應變量,X表示自變量,邏輯回歸用(1)式表達:

        上面(8)式即為邏輯回歸的最后表達式,且有:,即與變量X為線性關系,而(9)式可以通過最大似然估計及梯度上升法[18]求解。

        1.2 信息價值和證據(jù)權重

        信息價值(Information Value,IV)是訓練模型之前所提取的有效的信息量,亦即提取有用的變量信息作為建模。當經(jīng)典的邏輯回歸算法用于自變量多的企業(yè)數(shù)據(jù)建立分類模型時,經(jīng)常存在諸多不足。譬如變量共線比較敏感、很難擬合真實的數(shù)據(jù)分布,因此需要對自變量進行特征篩選,因此引入信息價值。

        IV可以根據(jù)量化指標的大小來衡量自變量的預測能力,它基于信息熵作為測量單位。信息熵用來描述“一個系統(tǒng)的混亂程度”,通常是度量樣本集合純度的一種指標。所謂“純”,就是讓分類器的一個分支擁有相同的屬性。當熵為1時,表示企業(yè)數(shù)據(jù)中某個變量趨于一個分段,此時信息價值IV較小;當熵為0時,表示企業(yè)數(shù)據(jù)中某個變量存在多個分段,此時信息價值IV較大。

        設X是x個企業(yè)數(shù)據(jù)樣本的集合,具有個不同的連續(xù)變量,其概率密度函數(shù)為,它的信息熵表示如下[19]:

        (10)式中E表示數(shù)學期望。

        由于信息熵能夠精確地描述企業(yè)信用有風險與無風險樣本分布,本文將其引入企業(yè)信用風險評估中。對于給定的變量x,按照企業(yè)數(shù)據(jù)集將其分成兩部分,對應有信用風險的數(shù)據(jù)密度概率為,對應無信用風險的數(shù)據(jù)密度概率為,即IV值由有信用風險的數(shù)據(jù)分布與無信用風險的數(shù)據(jù)分布之比的信息熵加上無信用風險的數(shù)據(jù)分布與有信用風險的數(shù)據(jù)分布之比的信息熵[20]:

        (13)式表明WOE為企業(yè)無信用風險概率與企業(yè)有信用風險概率之比取自然對數(shù)。到此,可以將(12)、(13)式寫在一起,得到IV與WOE的關系式:

        在企業(yè)信用風險評估模型中,我們需要選擇具有好的預測能力變量,即要求信息價值值大。值越大,說明該變量對無信用風險的企業(yè)與有信用風險企業(yè)的區(qū)分度越大。從(14)關系式中可知,為了使得較大,必須保證該變量的WOE與分布函數(shù)具有一致性,亦即WOE與企業(yè)的數(shù)據(jù)分布單調一致,否則需要根據(jù)實際業(yè)務中的情況進行具體的分析。

        2? 建模與過程

        本文的企業(yè)信用風險評估使用評分卡的方式建立模型。評分卡是當前運用比較廣泛的信用風險評價方法,其原理是將評價指標以證據(jù)權重(WOE)編碼方式離散化之后,再運用邏輯回歸進行模型訓練。該方法的特點為可解釋性強,模型結果穩(wěn)定。

        整個建模過程包括數(shù)據(jù)獲取和整合、目標確定和訓練窗口選擇、變量分段并計算信息價值、模型訓練及模型評價。

        2.1 數(shù)據(jù)獲取和整合

        根據(jù)企業(yè)該領域的特殊性,文本在建模中抽取以下三類變量:

        (1)原始變量。這些變量從數(shù)據(jù)庫中提取,概括了基本情況、各類原始明細等。通常來說,這些變量易于理解,但往往不是對模型最有效的;

        (2)衍生業(yè)務指標。這些字段來源于原始變量,有明確業(yè)務含義(如企業(yè)過去3年被投訴舉報數(shù)量),這些變量通常比原始變量有更好地預測能力,也是評分卡最常用指標。然而,這些數(shù)據(jù)通常不自然存在于數(shù)據(jù)庫中,需要在實施階段通過計算得到;

        (3)衍生模型指標。這些變量由分析人員對原始數(shù)據(jù)進行多層轉換和計算,通常由其他模型計算得出(如輿情風險指數(shù)),這些變量的預測能力一般是很好的,但和衍生業(yè)務指標相比,最大的區(qū)別在于缺少簡單清晰的業(yè)務含義,由模型高度抽象后計算得出。

        2.2 目標確定和訓練窗口選擇

        針對過去出現(xiàn)過信用風險事件的企業(yè),其目標變量標記為1,過去未出現(xiàn)過信用風險事件的企業(yè),其目標變量標記為0。

        評分卡模型是用特定時間段的數(shù)據(jù)訓練的,這個特定時間段被稱為建模窗口。選用窗口過短,企業(yè)數(shù)據(jù)在短期內(nèi)波動頻繁,容易引起模型預測結果的偏差,窗口過長則模型對企業(yè)及周圍環(huán)境的反應容易產(chǎn)生滯后現(xiàn)象。本文以1年為上限,即模型將預測企業(yè)在未來一年內(nèi)是否發(fā)生信用風險事件的概率。

        2.3 變量分段并計算信息價值

        實際的企業(yè)信用風險數(shù)據(jù)集中變量有兩種情況,分別是定性變量和定量變量。

        對于類別變量,已分好段,直接計算IV值;對于連續(xù)變量和離散變量,假設將變量X2分成K段,我們需要尋找K-1分點使得WOE與該變量保持單調性。必須注意的是變量X2可能有多種劃分保證這種單調性,這時我們選取IV的最大值作為最優(yōu)分箱結果。

        分段后,設PBK、PGK分別表示變量X2第K段對應目標變量為1(有信用風險)和0(無信用風險)的比例, 由(13)、(14)式得到:

        其中,

        有了每個分段的WOEK后,變量X2的信息價值的定義如下:

        計算完IV后,根據(jù)實際企業(yè)業(yè)務情況,選擇IV≥0.01的變量入模,因為IV在0.01以下的變量幾乎毫無預測能力,可舍棄。

        2.4 模型訓練

        模型的理論基礎為邏輯回歸算法,計算模型事件發(fā)生(本文y=1,即出現(xiàn)過信用風險事件)的概率,有(8)式展開得:

        這里,表示截距,為邏輯回歸中的系數(shù),由最大似然法(ML)求解,為原始變量經(jīng)轉換后的WOE值。

        模型開始訓練時,通常會選擇用逐步回歸對變量持續(xù)進行篩選,每一步都移入對模型預測能力有幫助的變量,同時移除對模型無增益的變量。

        除了模型本身對變量的選擇外,還考慮其它因素對變量進行選擇,如下幾項:

        (1)變量獨立:即共線問題?;貧w中的多重共線性是一個當模型中一些預測變量與其他預測變量高度相關時發(fā)生的條件。嚴重的多重共線性可能會產(chǎn)生問題,因為它可以增大回歸系數(shù)的方差,使它們變得不穩(wěn)定。

        (2)變量一致性:即變量訓練出來的參數(shù)正負性,和變量與目標變量的相關系數(shù)正負性必須一致,否則說明變量有偏差,需剔除;

        (3)變量可解釋性:即變量及其變化趨勢是可以被業(yè)務理解和使用的,而不是完全黑盒不可解釋,或者變量趨勢無業(yè)務含義。

        2.5 模型評價

        模型性能的好壞決定其在實際中的應用效果,良好的模型會真實地反映出企業(yè)的發(fā)展趨勢。本文采用“柯爾莫哥洛夫-斯米爾諾夫曲線”和“受試者工作特征”對企業(yè)信用風險模型進行性能評價。

        (1)柯爾莫哥洛夫-斯米爾諾夫曲線(Kolmogorov-Smirnov,K-S):將總體樣本進行n等分(通常選用較大數(shù)值,如1000),并按模型計算出的違約概率降序排序,計算每一等分中信用風險事件發(fā)生的累計百分比(Bad)和信用風險事件未發(fā)生的累計百分比(Good),繪制出兩者的差異,驗證兩者的一致性。

        (2)受試者工作特征(Receiver Operating Characteristic,ROC)曲線:通過在0到1之間改變用于創(chuàng)建混淆矩陣(Confusion Matrix)的臨界值,繪制Sensitivity(靈敏度)與1-Specificity(1-特殊度)曲線。ROC曲線下的面積為AUC統(tǒng)計量,一般認為AUC統(tǒng)計量≥0.75時,建立的模型具有較好的預測能力。

        3? 實例分析

        3.1 數(shù)據(jù)說明與指標體系

        本節(jié)先根據(jù)“國家‘互聯(lián)網(wǎng)+監(jiān)管系統(tǒng)企業(yè)信用風險分類評價標準”建立企業(yè)信用風險指標體系,再從廣東省各個行政部門獲取過去一年內(nèi)的企業(yè)有關數(shù)據(jù)集共1000條記錄。我們隨機選擇700條記錄作為模型訓練集;剩余300條記錄作為測試集,用來驗證模型的性能及有效性。其數(shù)據(jù)格式與指標體系分別如表1、表2所示。

        3.2 前期工作

        對1000條企業(yè)信用記錄,鑒于每條記錄擁有51個變量,可以考慮在建模前對了進行一些預處理工作。如下是實例分析前對數(shù)據(jù)集進行的預處理說明。

        (1)變量的相同率:指的是某個數(shù)據(jù)集中某個變量的值有多少個是相同的。高度一致的變量值包含的信息量少,無法區(qū)分各條記錄的差異,因此會導致該變量對模型的建立不起作用。本文定義變量的相同率閾值identical_limit0.9,高于0.9的變量直接刪除;

        (2)缺失值:這是數(shù)據(jù)挖掘建模中不可避免的步驟,造成數(shù)據(jù)缺失的原因是多方面的,文本根據(jù)企業(yè)業(yè)務情況,某變量的缺失率閾值missing_limit0.8,高于0.8的變量直接刪除,其余用到如圖1的方式處理。

        3.3 結果分析

        實例分別運用經(jīng)典邏輯回歸及證據(jù)權重邏輯回歸對企業(yè)信用風險數(shù)據(jù)集建模,并在訓練集、測試集采用(Kolmogorov-Smirnov,K-S)、(Receiver Operating Characteristic,ROC)曲線檢驗模型的有效性。

        圖2、圖3為證據(jù)權重邏輯回歸法在訓練集、測試集的K-S、ROC曲線圖。其中,圖中的虛線為訓練集、測試集的Good數(shù)據(jù)(企業(yè)無信用風險)與Bad(企業(yè)有信用風險)的累積概率分布,其意義表示兩個分布函數(shù)是否有差異,用p值衡量,并設置顯著水平為λ=0.05。這里,模型在訓練集、測試集均有,說明兩者來自同一分布。此外,K-S值越大模型性能越好,而模型在訓練集、測試集中分別有K-S=0.5541、K-S=0.5404,查看表3的參考值可知,模型具有非常好的分辨力能。再者,模型在訓練集、測試集的ROC曲線都有不錯的表現(xiàn),由其與橫坐標圍成的圖形面積值分別為AUC=0.8279、AUC=0.8097,均大于0.75,說明模型的二分類性能良好。最后,我們根據(jù)二分類的混合矩陣分析模型在預測陰性(文本為‘無信用風險企業(yè))和陽性(文本為‘有信用風險企業(yè))的單側能力,模型的混合矩陣如表4、表5所示。從表中看到模型在訓練集中識別“有信用風險”的記錄能力為88.1%,201條‘有風險的記錄只有24條記錄被錯誤預測為‘無風險;在測試集中識別“有信用風險”的記錄能力為83.0%,87條‘有風險的記錄只有14條記錄被錯誤預測為‘無風險;同時,模型在上述數(shù)據(jù)集的整體準確率分別為85.6%、82.5%。

        圖4、圖5為經(jīng)典邏輯回歸算法在訓練集、測試集的K-S、ROC曲線圖,其在訓練集和測試集K-S的值分別為0.5225、0.4306;ROC中的AUC值分別為0.8387、0.7667,很明顯經(jīng)典邏輯算法在訓練集上有過擬合現(xiàn)象。

        不論從K-S值或ROC中的AUC值比較,實例證明帶證據(jù)權重邏輯回歸法在企業(yè)信用風險模型評估中顯然由于經(jīng)典的邏輯回歸算法。

        4? 結語

        本文基于廣東省有關行政部門的真實企業(yè)數(shù)據(jù),應用證據(jù)權重邏輯回歸對其進行信用風險建模。通過對模型的K-S值和ROC中AUC值作了深入分析,并與經(jīng)典的邏輯回歸做比較,實例驗證了帶證據(jù)權重邏輯回歸法在預測企業(yè)信用方面的有效性。展望未來,我們會應用決策樹算法的信息增益結合證據(jù)權重再次探討它們在企業(yè)信用風險評估中的效果,然后做一個全面的分析、歸納。

        參考文獻

        [1] 國務院辦公廳關于推廣隨機抽查規(guī)范事中事后監(jiān)管的通知[EB/OL].http://www.gov.cn/zhengce/content/2015-08/05/content_10051.htm.

        [2] 國務院關于“先照后證”改革后加強事中事后監(jiān)管的意見[EB/OL].http://www.gov.cn/zhengce/content/2015-11/03/content_10263.htm.

        [3] 國務院關于印發(fā)2016年推進簡政放權放管結合優(yōu)化服務改革工作要點的通知[EB/OL].http://www.gov.cn/zhengce/content/2016-05/24/content_5076241.htm.

        [4] Odom M D,Sharda R.A neural network model for bankruptcy prediction[C]// IJCNN International Joint Conference on Neural Networks.IEEE, 1990:163-168.

        [5] Prinzie A,Poel D V D. Random forest for multiclass classification: random multinomial logit[J]. Working Papers of Faculty of Economics & Business Administration Ghent University Belgium, 2008,34(3):1721-1732.

        [6] Yeh C C, Lin F, Hsu C Y. A hybrid KMV model, random forests and rough set theory approach for credit rating[J]. Knowledge-Based Systems, 2012, 33(3):166-172.

        [7] Kartasheva A V, Traskin M. Insurers insolvency prediction using random forest classification[J]. Social Science Electronic Publishing, 2013, 10(3): 16-62.

        [8] 吳麗麗. 基于Logistic回歸模型的商業(yè)銀行信用風險管理研究[D].哈爾濱:哈爾濱工業(yè)大學,2007.

        [9] 郭玉華. 微型企業(yè)信用風險評估——基于Logit模型的分析[J]. 經(jīng)濟論壇,2011(11):213-216.

        [10] 方匡南,范新妍,馬雙鴿.基于網(wǎng)絡結構Logistic模型的企業(yè)信用風險預警[J].統(tǒng)計研究,2016,33(4):50-55.

        [11] 楊俊,夏晨琦.基于Gradient Boosting算法的小企業(yè)信用風險評估[J].浙江金融,2017(9):44-50.

        [12] 熊正德,張帆,熊一鵬.引入WFCM算法能提高信用違約測度模型準確率嗎?——以滬深A股制造業(yè)上市公司為樣本的實證研究[J]. 財經(jīng)理論與實踐,2018(1):147-153.

        [13] 基于WOE-Probit逐步回歸的信用指標組合篩選模型及應用[J]. 管理科學,2018,48(2):76-87.

        [14] 劉麗君,韓靜磊,錢益斌,等.利用靶器官毒性劑量法(TTD)和證據(jù)權重分析法(WOE) 評估固化飛灰中重金屬非致癌健康風險[J]. 環(huán)境化學,2019,38(5):1014-1020.

        [15]趙雅迪,吳釗,李慶兵,等.電費回收風險預測的大數(shù)據(jù)方法應用[J]. 電信科學,2019,35(2):125-133.

        [16] 陳超,王楠,于海洋,等.基于卡方分箱法和邏輯回歸算法的轉爐操作工藝評價模型[J]. 材料與冶金學報,2019,18(2):87-91.

        [17] 危明鑄,麥偉杰,袁峰,等.基于機器學習企業(yè)運行風險研究[J]. 軟件,2019,40(8):29-37.

        [18] H?rdleW, Sperlich S, et al. Nonparametric and semiparametric models[M]. Berlin: Springer, 2004:145-165.

        [19] Siddiqi N. Credit risk scorecards[M]. John Wiley & Sons, 2006: 70-125.

        [20] Good I J. Weight of evidence: a brief survey[J]. Bayesian statistics. North Holland, Amsterdam 1983, 2: 249-269.

        猜你喜歡
        企業(yè)信用信用風險
        揚州市穩(wěn)步推進安全生產(chǎn)領域企業(yè)信用修復
        泰州市推行企業(yè)信用修復全鏈條服務模式
        “涉軍”企業(yè)信用評價擴大試點工作即將啟動
        中國軍轉民(2017年6期)2018-01-31 02:22:13
        影響企業(yè)信用政策制定的因素分析
        消費導刊(2017年20期)2018-01-03 06:27:45
        “涉軍”企業(yè)信用評論擴大試點工作即將啟動
        中國軍轉民(2017年8期)2017-12-13 08:37:01
        淺析我國商業(yè)銀行信用風險管理
        京東商城電子商務信用風險防范策略
        PPP項目發(fā)行中期票據(jù)的可能性與信用風險分析
        個人信用風險評分的指標選擇研究
        上市公司信用風險測度的不確定性DE-KMV模型
        日韩女同在线免费观看| 久久精品国产91久久性色tv| 国产亚洲欧美另类久久久| 白色月光免费观看完整版| 欧美精品色婷婷五月综合| 久久精品麻豆日日躁夜夜躁| 国产农村妇女毛片精品久久久| 色中文字幕视频在线观看| 人成在线免费视频网站| 99精品国产一区二区三区不卡| 亚洲国产精品久久久久久久| 亚欧免费无码AⅤ在线观看| 精品乱色一区二区中文字幕| 九九久久自然熟的香蕉图片| 香蕉视频www.5.在线观看| 亚洲男人在线无码视频| 风流熟女一区二区三区| 无码视频在线观看| 久久AⅤ无码精品为人妻系列| 亚洲区一区二区中文字幕| 日本午夜理论片在线观看| 久激情内射婷内射蜜桃| 久久久亚洲欧洲日产国产成人无码| 亚洲成生人免费av毛片| 亚洲youwu永久无码精品| 日韩精品无码视频一区二区蜜桃 | 婷婷亚洲综合五月天小说| 国产一区二区三区高清视频| 国产免费二区三区视频| 中文字幕+乱码+中文字幕一区| 白浆出来无码视频在线| 蜜桃视频永久免费在线观看| 精品av熟女一区二区偷窥海滩| 无码人妻一区二区三区免费| 亲少妇摸少妇和少妇啪啪| 青青草骚视频在线观看 | 秋霞午夜无码鲁丝片午夜精品| 中文字幕日韩人妻在线| 中文字幕亚洲乱码成熟女1区| 日韩精品无码视频一区二区蜜桃| 亚洲欧美变态另类综合|