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

    大壩變形監(jiān)控模型識(shí)別的R-OC準(zhǔn)則

    2023-05-26 09:12:18吳震宇陳建康
    工程科學(xué)與技術(shù) 2023年3期
    關(guān)鍵詞:監(jiān)測(cè)數(shù)據(jù)時(shí)效過(guò)度

    張 博,劉 健,吳震宇*,陳建康,尹 川

    (1.四川大學(xué) 水力學(xué)與山區(qū)河流開(kāi)發(fā)保護(hù)國(guó)家重點(diǎn)實(shí)驗(yàn)室 水利水電學(xué)院,四川 成都 610065;2.雅礱江流域水電開(kāi)發(fā)有限公司,四川 成都 610051)

    中國(guó)目前已修建水庫(kù)大壩9.8萬(wàn)余座,帶來(lái)巨大社會(huì)經(jīng)濟(jì)效益的同時(shí),其安全風(fēng)險(xiǎn)也不容忽視[1]。變形監(jiān)測(cè)是大壩長(zhǎng)期運(yùn)行安全的重要手段[2-3],常通過(guò)分析變形監(jiān)測(cè)數(shù)據(jù)來(lái)診斷壩體結(jié)構(gòu)性態(tài),預(yù)測(cè)大壩發(fā)展趨勢(shì)。而構(gòu)建監(jiān)控模型是分析大壩變形監(jiān)測(cè)數(shù)據(jù)的常用方法,通常根據(jù)監(jiān)控模型預(yù)測(cè)的預(yù)警界限對(duì)大壩變形進(jìn)行異常預(yù)警,因此,構(gòu)建大壩變形監(jiān)控模型對(duì)于大壩安全管理具有重要意義。

    大壩變形監(jiān)測(cè)數(shù)據(jù)序列一般由兩部分組成,即由水庫(kù)水位和環(huán)境溫度變化引起的可逆變形和大壩隨時(shí)間演變的不可逆變形?;诒O(jiān)測(cè)資料,應(yīng)用回歸方法進(jìn)行分析,可將影響大壩變形的因素歸結(jié)為水壓、溫度和時(shí)效分量[4]。在3種分量中,水壓分量一般采用上游相對(duì)水深3次冪多項(xiàng)式的形式[5]。溫度分量可由多個(gè)周期諧波的三角函數(shù)[6]、變形測(cè)量前多段平均氣溫的線性組合[7]構(gòu)成,或?qū)误w溫度計(jì)監(jiān)測(cè)數(shù)據(jù)進(jìn)行主成分分析,將主成分的線性組合作為溫度分量[8]。時(shí)效分量一般通過(guò)預(yù)設(shè)不同函數(shù)進(jìn)行模擬,多采用線性函數(shù)和對(duì)數(shù)函數(shù)組合[9-11]、線性函數(shù)和指數(shù)函數(shù)組合[12-14]、指數(shù)和多項(xiàng)式組合[15-16]等形式構(gòu)建,也可采用EMD方法分離出大壩變形監(jiān)測(cè)序列時(shí)效分量,使非平穩(wěn)的大壩變形時(shí)間序列平穩(wěn)化,能有效提高預(yù)測(cè)精度[17-18]。

    對(duì)于同一個(gè)變形測(cè)點(diǎn),通過(guò)不同形式的水壓、溫度和時(shí)效分量的組合,可以建立多種大壩變形監(jiān)控模型,因此,需要從所有待選模型中識(shí)別出最優(yōu)模型。大壩變形監(jiān)控模型的優(yōu)劣主要體現(xiàn)為模型泛化能力的高低。模型泛化能力即模型對(duì)訓(xùn)練集外樣本的預(yù)測(cè)能力,欠擬合與過(guò)度擬合均會(huì)導(dǎo)致監(jiān)控模型的泛化能力不足[19]。欠擬合是指建立的監(jiān)控模型不能從歷史監(jiān)測(cè)數(shù)據(jù)充分學(xué)習(xí)和準(zhǔn)確描述大壩變形監(jiān)測(cè)量的統(tǒng)計(jì)變化規(guī)律,具體表現(xiàn)為監(jiān)控模型對(duì)歷史監(jiān)測(cè)數(shù)據(jù)的回歸顯著性較低,擬合誤差較大,常通過(guò)復(fù)相關(guān)系數(shù)R、均方根誤差RMSE、平均絕對(duì)誤差MAE、平均絕對(duì)百分比誤差MAPE等指標(biāo)進(jìn)行評(píng)判[20-21]。R越大,表明監(jiān)控模型對(duì)歷史監(jiān)測(cè)數(shù)據(jù)的回歸顯著性越高,RMSE、MAE、MAPE的值越小,表明模型擬合誤差越小。而過(guò)度擬合表現(xiàn)為模型對(duì)監(jiān)測(cè)序列的過(guò)度學(xué)習(xí),甚至?xí)蹲降接?xùn)練集中單個(gè)樣本的特征擬合抽樣誤差,從而導(dǎo)致模型泛化能力降低,出現(xiàn)對(duì)監(jiān)測(cè)數(shù)據(jù)異常的誤判。AIC、BIC準(zhǔn)則被廣泛應(yīng)用于監(jiān)控模型識(shí)別[22-25],均采用似然函數(shù)度量模型擬合誤差,同時(shí)引入模型參數(shù)數(shù)量的懲罰項(xiàng)。經(jīng)實(shí)際應(yīng)用檢驗(yàn),AIC準(zhǔn)則適用于樣本量較小的情況,BIC準(zhǔn)則適用于樣本量較大的情況。但AIC、BIC準(zhǔn)則不能定量比較和評(píng)價(jià)模型的過(guò)度擬合程度,在某些情況下,識(shí)別出的最優(yōu)模型可能存在嚴(yán)重的過(guò)度擬合問(wèn)題。因此,在進(jìn)行模型選擇時(shí),需要能夠量化模型過(guò)度擬合程度的識(shí)別準(zhǔn)則。

    本文定義了能夠量化模型過(guò)度擬合程度的指標(biāo)—過(guò)度擬合系數(shù)OC,同時(shí)采用復(fù)相關(guān)系數(shù)R刻畫(huà)模型的擬合精度,提出了大壩變形監(jiān)控模型識(shí)別的R-OC準(zhǔn)則。將R-OC準(zhǔn)則應(yīng)用于某大壩垂直位移測(cè)點(diǎn)的變形監(jiān)控模型優(yōu)選,并與AIC和BIC準(zhǔn)則進(jìn)行對(duì)比,驗(yàn)證了R-OC準(zhǔn)則的有效性。

    1 大壩變形監(jiān)控模型

    大壩變形監(jiān)控模型一般由水壓分量、溫度分量和時(shí)效分量構(gòu)成,表達(dá)式如下:

    式中,Y(t) 為 壩體變形監(jiān)測(cè)值在時(shí)間t的統(tǒng)計(jì)估計(jì)值,F(xiàn)1[H(t)]為 壩體變形的水壓分量,F(xiàn)2[T(t)]為壩體變形的溫度分量,F(xiàn)3[θ(t)]為 壩體變形的時(shí)效分量,C為待定常數(shù)項(xiàng)。

    1.1 水壓分量

    大壩變形的水壓分量與上、下游水位的乘冪相關(guān),對(duì)水位進(jìn)行歸一化處理[26],構(gòu)造形式如下:

    式中,aui、adi分別為上、下游水位因子的待定回歸系數(shù),Ht、ht分 別為t時(shí) 刻的上、下游水位,Hmin、Hmax分別為大壩運(yùn)行時(shí)期最低、最高上游庫(kù)水位,hmin、hmax分別為大壩運(yùn)行時(shí)期最低、最高下游庫(kù)水位。

    當(dāng)水位對(duì)大壩變形影響存在滯后時(shí),水壓分量與觀測(cè)日前j天的平均水位呈線性關(guān)系:

    1.2 溫度分量

    1.2.1 傳統(tǒng)溫度分量

    1)多個(gè)周期諧波三角函數(shù)

    引起壩體位移的溫度因素主要是邊界溫度及壩體混凝土的水化熱。對(duì)于運(yùn)行時(shí)間較長(zhǎng)的大壩,水化熱可忽略不計(jì)。壩體溫度基本呈年周期性變化。因此,可采用多個(gè)周期諧波的三角函數(shù)作為統(tǒng)計(jì)模型的溫度分量,具體構(gòu)造如下:

    式(4)~(5)中,bi(i=1,2,3,4)為 待定回歸系數(shù),t0為設(shè)定的1月1日的初始時(shí)間。

    2)多段平均氣溫線性組合

    壩體邊界溫度的變化都是由于氣溫的季節(jié)性變化引起。當(dāng)外部環(huán)境氣溫監(jiān)測(cè)序列完整時(shí),溫度分量可采用多段平均氣溫的線性組合,函數(shù)構(gòu)造如下:

    式中,bi(i=1,2,3,4)為 待定回歸系數(shù),Ti代表觀測(cè)日當(dāng)天及前5、15、60 d的平均大氣氣溫。

    1.2.2 基于主成分分析的溫度分量

    當(dāng)壩體埋設(shè)的溫度計(jì)較多且監(jiān)測(cè)數(shù)據(jù)連續(xù)時(shí),可采用壩體溫度監(jiān)測(cè)序列的主成分構(gòu)造大壩變形的溫度分量,能有效減少模型的變量和參數(shù),降低模型的過(guò)度擬合程度。

    假設(shè)壩體埋設(shè)了p支溫度計(jì),每支溫度計(jì)有n(n>p)個(gè)監(jiān)測(cè)數(shù)據(jù),表達(dá)為矩陣形式:

    式中,xi=(x1i,x2i,···,xni)T,i=1,2,···,p。

    并計(jì)算相關(guān)系數(shù)矩陣:

    為消除量綱不同帶來(lái)的影響,對(duì)數(shù)據(jù)X進(jìn)行標(biāo)準(zhǔn)化處理,使變量的平均值為0,方差為1,標(biāo)準(zhǔn)化公式如下:

    式中:

    計(jì)算相關(guān)系數(shù)矩陣的特征值和對(duì)應(yīng)的特征向量,依次排列大小,相應(yīng)的單位特征向量如下:

    則第i個(gè)主成分PCi(i=1,2,···,p)的表達(dá)式為:

    式中,PCi(i=1,2,···,p)為對(duì)篩選出的壩體溫度計(jì)進(jìn)行主成分分析后得到的主成分因子。

    通過(guò)式(14)計(jì)算各主成分貢獻(xiàn)率,并由大到小累加,根據(jù)式(15)計(jì)算累積貢獻(xiàn)率。

    式(14)~(15)中,λi為 特征值,ηj為各主成分貢獻(xiàn)率,Mm為主成分累積貢獻(xiàn)率。

    選取合適的主成分PCi作為因子,進(jìn)行溫度分量的構(gòu)造,構(gòu)造形式如下:

    式中:bi(i=1,2,3,4)為 待定回歸系數(shù);m為所選取的主成分因子個(gè)數(shù),一般選取累積貢獻(xiàn)率大于85%或95%的主成分個(gè)數(shù)。

    1.3 時(shí)效分量

    1.3.1 傳統(tǒng)時(shí)效分量

    時(shí)效分量反映了大壩的不可逆變形。常采用如下兩種形式來(lái)表示:

    式(17)~(18)中,θ=(t-t0)/100,t0為監(jiān)測(cè)點(diǎn)初始監(jiān)測(cè)日的時(shí)間,c1、c2為待定回歸系數(shù)。

    1.3.2 基于經(jīng)驗(yàn)?zāi)B(tài)分解的時(shí)效分量

    通過(guò)經(jīng)驗(yàn)?zāi)B(tài)分解法(EMD)將大壩變形監(jiān)測(cè)效應(yīng)量分解并重構(gòu)為周期變化項(xiàng)和時(shí)效趨勢(shì)項(xiàng),選擇合適的函數(shù)形式對(duì)時(shí)效趨勢(shì)項(xiàng)進(jìn)行擬合,構(gòu)建時(shí)效分量模型。該方法能克服傳統(tǒng)方法在未知時(shí)效分量實(shí)際變化規(guī)律和趨勢(shì)情況下預(yù)設(shè)模型形式的缺點(diǎn)。

    EMD的步驟如下:

    1)找出x(t)上所有極值點(diǎn)xmax(t)、xmin(t)。

    2)通過(guò)3次樣條差值函數(shù)將極值點(diǎn)連接,繪制上、下包絡(luò)線。上、下包絡(luò)線的均值m1(t)可以表示為:

    3)計(jì)算原始數(shù)據(jù)和m1(t)之 間的差值h1(t):

    4)用h1(t)替 換x(t),重復(fù)上述3個(gè)步驟,直到滿足相應(yīng)的終止條件:①數(shù)據(jù)序列的零點(diǎn)數(shù)m和極值點(diǎn)數(shù)n滿足不等式|m-n|≤1;②上下包絡(luò)線的均值m1(t)趨于0。

    第1個(gè)imf1(t)分量可表示為:

    5)計(jì)算原始數(shù)據(jù)和imf1(t)之 間的差值r1(t):

    6)將步驟1)中的x(t)替 換為r1(t),重復(fù)步驟1)~5),直到rn(t)殘差單調(diào)或小于預(yù)設(shè)誤差,EMD分解終止。

    變形監(jiān)測(cè)數(shù)據(jù)時(shí)間序列EMD分解的數(shù)學(xué)表達(dá)式如下:

    式(19)~(23)中,x(t) 為變形監(jiān)測(cè)測(cè)點(diǎn)在時(shí)間t的監(jiān)測(cè)效應(yīng)量,imfi(t) 為 通過(guò)EMD分解得到的n階IMF分量,rn(t) 為 監(jiān)測(cè)數(shù)據(jù)序列x(t)的時(shí)效趨勢(shì)項(xiàng)。

    則時(shí)效分量F3(θ(t))可表示為:

    對(duì)于時(shí)效分量F3(θ(t)),本文采用如下復(fù)合指數(shù)函數(shù)進(jìn)行擬合:

    式中,a、b、c、d為待定擬合參數(shù),θ=(t-t0)/100,t0為監(jiān)測(cè)點(diǎn)初始監(jiān)測(cè)日的時(shí)間。

    2 模型識(shí)別準(zhǔn)則

    2.1 信息準(zhǔn)則

    2.1.1 AIC準(zhǔn)則

    AIC準(zhǔn)則是衡量統(tǒng)計(jì)模型擬合優(yōu)良性的一種標(biāo)準(zhǔn),由日本統(tǒng)計(jì)學(xué)家Akaike在1974年提出,公式如下:

    式中,k為模型參數(shù)個(gè)數(shù),L為似然函數(shù)。

    從待選模型中選擇最佳模型時(shí),通常選擇AIC值最小的模型。一般而言,當(dāng)模型復(fù)雜度提高時(shí),似然函數(shù)L增大,從而AIC值減小;但當(dāng)k過(guò)大時(shí),似然函數(shù)增速減緩,導(dǎo)致AIC增大。因此模型過(guò)于復(fù)雜容易造成過(guò)度擬合現(xiàn)象,而 2k可以作為對(duì)參數(shù)的懲罰。AIC準(zhǔn)則的實(shí)質(zhì)就是由候選模型取最大似然后加上一個(gè)懲罰項(xiàng)得到的,可以衡量模型參數(shù)復(fù)雜度與擬合優(yōu)度之間的平衡關(guān)系,并選擇模型參數(shù)較少但擬合度較好的模型作為最優(yōu)模型。

    2.1.2 BIC準(zhǔn)則

    BIC準(zhǔn)則是Schwartz于1978年根據(jù)Bayes理論提出的判別準(zhǔn)則,引入了后驗(yàn)概率的思想,其計(jì)算公式如下:

    式中,k為模型參數(shù)個(gè)數(shù),n為 樣本數(shù)量,L為似然函數(shù)。

    BIC準(zhǔn)則與AIC準(zhǔn)則的差別主要體現(xiàn)在懲罰項(xiàng),考慮到樣本個(gè)數(shù)會(huì)對(duì)模型選擇起到一定的影響作用,所以在BIC準(zhǔn)則中加入樣本量n作為模型懲罰項(xiàng)的一部分;當(dāng)樣本數(shù)量過(guò)多時(shí),可有效防止模型精度過(guò)高造成的模型復(fù)雜度過(guò)高。最優(yōu)模型通常為 B IC值最小的模型。

    2.2 R-OC準(zhǔn)則

    過(guò)度擬合是導(dǎo)致模型泛化能力降低的主要原因之一,通常表現(xiàn)為模型對(duì)歷史數(shù)據(jù)的擬合精度很高,但預(yù)測(cè)精度明顯低于擬合精度。因此,可以通過(guò)比較模型預(yù)測(cè)誤差和擬合誤差來(lái)確定模型的過(guò)度擬合程度。本文提出模型過(guò)度擬合系數(shù)OC,其表達(dá)式如下:

    式中,RMSE為均方根誤差,MAE為平均絕對(duì)誤差,MAPE為平均絕對(duì)百分比誤差,下標(biāo)V和下標(biāo)F分別代表驗(yàn)證時(shí)段和擬合時(shí)段。相應(yīng)的RMSE、MAE及MAPE的值定量反映了監(jiān)控模型預(yù)測(cè)誤差和擬合誤差的大小。

    RMSE、MAE和MAPE是常用的誤差評(píng)價(jià)指標(biāo)[21],但各有優(yōu)缺點(diǎn)。RMSE能夠很好地反映擬合值與監(jiān)測(cè)數(shù)據(jù)的偏差,但對(duì)于遠(yuǎn)離監(jiān)測(cè)數(shù)據(jù)均值的異常值很敏感,其表達(dá)式如下:

    MAE克服了RMSE的缺點(diǎn),能更好地反映擬合值誤差的實(shí)際情況,但不能合理反映模型擬合的相對(duì)誤差,其表達(dá)式如下:

    MAPE適合評(píng)價(jià)相對(duì)誤差,缺點(diǎn)是當(dāng)實(shí)測(cè)值為0時(shí),將無(wú)法計(jì)算,其表達(dá)式如下:

    當(dāng)模型過(guò)度擬合時(shí),驗(yàn)證時(shí)段的模型預(yù)測(cè)誤差將大于擬合時(shí)段的擬合誤差,即:當(dāng)OC>1,且OC的值越大,表明模型過(guò)度擬合程度越高;反之,當(dāng)OC<1時(shí),且OC的值越小,表明模型預(yù)測(cè)精度高于擬合精度,不存在過(guò)度擬合,預(yù)測(cè)精度相對(duì)越好;當(dāng)OC=1時(shí),表示預(yù)測(cè)精度和擬合精度相同。因此,OC能定量反映監(jiān)控模型的過(guò)度擬合程度。

    而進(jìn)行模型選擇時(shí),既要保證模型驗(yàn)證時(shí)段數(shù)據(jù)的預(yù)測(cè)精度,不能存在過(guò)度擬合;也要保證模型擬合時(shí)段數(shù)據(jù)的擬合精度,不能存在欠擬合。因此本文選用復(fù)相關(guān)系數(shù)R定量反映模型擬合時(shí)段數(shù)據(jù)的擬合精度。

    復(fù)相關(guān)系數(shù)R反映一個(gè)因變量與多個(gè)自變量之間的相關(guān)程度,其表達(dá)式如下:

    式中,R2為可決系數(shù),SSR為回歸平方和,SST為總偏差平方和,Yi為觀測(cè)數(shù)據(jù),Y為觀測(cè)數(shù)據(jù)均值,Y?i為擬合數(shù)據(jù)。

    回歸平方和SSR反映回歸方程中全部自變量的“方差貢獻(xiàn)”。可決系數(shù)R2表示這種貢獻(xiàn)在總偏差平方和SST中所占的比重,R2越大,說(shuō)明自變量對(duì)因變量的解釋程度越高,引起的變動(dòng)占總變動(dòng)的比重越大,模型擬合程度就越高;反之,說(shuō)明模型對(duì)樣本觀測(cè)值的擬合效果越差。因此,復(fù)相關(guān)系數(shù)R越大,監(jiān)控模型的擬合時(shí)段數(shù)據(jù)的擬合精度就越高。

    本文基于反映模型擬合精度的指標(biāo)R和反映過(guò)度擬合程度的指標(biāo)OC,提出了一種新的大壩變形監(jiān)控模型識(shí)別的R-OC準(zhǔn)則。在進(jìn)行模型選擇時(shí),可將指標(biāo)R和OC繪制為2維散點(diǎn)圖,從而可以直觀地挑選出擬合精度高、過(guò)度擬合程度低的模型。

    3 案例分析

    3.1 工程概況

    國(guó)內(nèi)某水電站工程,屬二等大(2)型工程,由左岸重力式擋水壩段、河床式電站廠房壩段、河床4孔泄洪閘壩段、右岸導(dǎo)流明渠內(nèi)3孔泄洪閘壩段、右岸重力式擋水壩段等建筑物組成。水庫(kù)正常蓄水位為1 015 m,死水位為1 012 m,壩頂高程1 020 m,最大壩高69.5 m,水庫(kù)具有日調(diào)節(jié)性能。

    選取10#壩段真空激光準(zhǔn)直系統(tǒng)LA14測(cè)點(diǎn)垂直位移構(gòu)建位移監(jiān)控模型,測(cè)點(diǎn)布置情況如圖1所示。10#壩段長(zhǎng)30 m,順?biāo)鞣较蜷L(zhǎng)60 m,最大閘壩高60 m,堰頂高程994 m。為監(jiān)測(cè)壩體運(yùn)行期溫度場(chǎng)的變化情況,在10#壩段壩體內(nèi)共埋設(shè)30支溫度計(jì),封存停測(cè)后納入運(yùn)行期統(tǒng)計(jì)共計(jì)14支,具體儀器布置情況如圖2所示。

    圖1 LA14測(cè)點(diǎn)布置及工程鳥(niǎo)瞰圖Fig.1 LA14 measurement point layout and engineering aeroview

    圖2 10#壩段溫度計(jì)布置圖Fig.2 10# dam section thermometer layout

    3.2 位移監(jiān)控模型

    在建模時(shí)段內(nèi)大壩正常運(yùn)行,工作性態(tài)良好,巡視檢查未發(fā)現(xiàn)異常的變形跡象。收集LA14測(cè)點(diǎn)垂直位移自2017年7月6日至2021年3月30日的完整監(jiān)測(cè)數(shù)據(jù),測(cè)量頻率為1 次/d。剔除異常值后,獲得1 238組實(shí)測(cè)序列,用以構(gòu)建模型的樣本觀測(cè)值。為檢驗(yàn)監(jiān)控模型的預(yù)測(cè)效果,將數(shù)據(jù)序列分成擬合時(shí)段和驗(yàn)證時(shí)段數(shù)據(jù)。擬合時(shí)段數(shù)據(jù)用于擬合構(gòu)建監(jiān)控模型,通過(guò)監(jiān)控模型對(duì)LA14測(cè)點(diǎn)位移進(jìn)行預(yù)測(cè)得到預(yù)測(cè)數(shù)據(jù),并同驗(yàn)證時(shí)段內(nèi)的LA14測(cè)點(diǎn)位移的實(shí)測(cè)值進(jìn)行對(duì)比分析。其中,擬合時(shí)段起止時(shí)間為2017年7月6日至2020年12月31日,驗(yàn)證時(shí)段起止時(shí)間為2021年1月1日至2021年3月30日。LA14測(cè)點(diǎn)全過(guò)程時(shí)間序列曲線如圖3所示。

    圖3 LA14測(cè)點(diǎn)全過(guò)程時(shí)間序列曲線Fig.3 Whole process time series curve of LA14 measuring point

    3.2.1 水壓分量

    水壓分量根據(jù)第1.1節(jié)中的形式進(jìn)行構(gòu)建,統(tǒng)一對(duì)上下游水位進(jìn)行歸一化處理,且需考慮下游水位的滯后性。

    3.2.2 溫度分量

    1)傳統(tǒng)溫度分量:由于該水電站運(yùn)行多年,水化熱已釋放完畢,壩體內(nèi)部溫度場(chǎng)已趨向穩(wěn)定,故在回歸計(jì)算時(shí),水化熱可不予考慮;且外部環(huán)境氣溫監(jiān)測(cè)序列完整。因此,溫度分量可分別采用第1.2節(jié)中的形式進(jìn)行構(gòu)建。

    2)基于PCA的溫度分量:通過(guò)數(shù)據(jù)篩選,選擇T01-XHZ、T04-XHZ、T05-XHZ、T10-XHZ、T14-XHZ、T15-XHZ、T16-XHZ、T17-XHZ、T21-XHZ、T23-XHZ、T28-XHZ、T29-XHZ共12支溫度計(jì),形成溫度主成分PC1、PC2、···、PC12,各主成分貢獻(xiàn)率及累積貢獻(xiàn)率如圖4所示。

    圖4 溫度監(jiān)測(cè)數(shù)據(jù)主成分的貢獻(xiàn)率及累積貢獻(xiàn)率Fig.4 Contribution rates and cumulative contribution rates of principal components in temperature monitoring data

    由圖4可得,當(dāng)選取前4支主成分PC1、PC2、PC3、PC4時(shí),累積貢獻(xiàn)率達(dá)到99.60%,可近似表示上述12支溫度計(jì)在所選時(shí)段的所有包含信息,達(dá)到降維處理的目的。因此,將主成分PC1、PC2、PC4作為溫度分量因子。

    3.2.3 時(shí)效分量

    1)傳統(tǒng)時(shí)效分量:時(shí)效分量可采用第1.3.1節(jié)中的形式進(jìn)行構(gòu)建。

    2)基于EMD的時(shí)效分量:對(duì)LA14測(cè)點(diǎn)監(jiān)測(cè)數(shù)據(jù)序列進(jìn)行EMD分解,根據(jù)EMD提取的時(shí)效位移時(shí)間序列,選擇復(fù)合指數(shù)函數(shù)構(gòu)建時(shí)效分量,采用遺傳算法確定復(fù)合指數(shù)函數(shù)的系數(shù)。LA14測(cè)點(diǎn)EMD時(shí)效分量及復(fù)合函數(shù)擬合效果如圖5所示。

    圖5 LA14測(cè)點(diǎn)EMD分解重構(gòu)及擬合Fig.5 EMD decomposition reconstruction and fitting of LA14 measuring point

    擬合結(jié)果顯示,復(fù)相關(guān)系數(shù)R>0.99,均方根誤差RMSE<0.01,即復(fù)合指數(shù)函數(shù)對(duì)分解出的壩頂時(shí)效位移擬合效果良好。

    3.3 待選監(jiān)控模型

    大壩變形監(jiān)控模型一般包含水壓、溫度和時(shí)效分量3部分,通過(guò)選取上述構(gòu)建的水壓、溫度和時(shí)效分量不同的函數(shù)形式進(jìn)行組合,并確定是否考慮下游水位滯后,即可構(gòu)建不同的監(jiān)控模型,本文總計(jì)建立18種大壩變形待選監(jiān)控模型。各監(jiān)控模型表達(dá)式見(jiàn)表1。

    表1 LA14測(cè)點(diǎn)垂直位移監(jiān)控模型匯總Tab.1 Summary of LA14 vertical displacement monitoring model

    3.4 模型識(shí)別結(jié)果及討論

    采用全回歸方法對(duì)LA14測(cè)點(diǎn)進(jìn)行擬合。在評(píng)判預(yù)測(cè)效果時(shí),使用核密度估計(jì)法進(jìn)行分布擬合,其表達(dá)式如下:

    式中,f?h(x)為 殘差樣本x的 分布概率,K、h分別為核函數(shù)、平滑參數(shù),xi為參與統(tǒng)計(jì)的數(shù)據(jù)信息。

    根據(jù)監(jiān)控模型對(duì)LA14測(cè)點(diǎn)位移的估計(jì)值和估計(jì)誤差的概率分布確定異常監(jiān)測(cè)數(shù)據(jù)的預(yù)警界限。其中,預(yù)警上限為U CL=E+Δ97.5%,預(yù)警下限為L(zhǎng)CL=E+Δ2.5%,E為L(zhǎng)A14測(cè)點(diǎn)位移的模型估計(jì)值,Δ2.5%和Δ97.5%分別為模型估計(jì)誤差的2.5%和97.5%分位值。當(dāng)新獲得的監(jiān)測(cè)數(shù)據(jù)超出預(yù)警界限時(shí),則判斷為異常數(shù)據(jù),并發(fā)出預(yù)警。

    由于案例分析中的重力壩在驗(yàn)證時(shí)段內(nèi)處于正常運(yùn)行狀態(tài),巡視檢查未見(jiàn)任何異常變形跡象,驗(yàn)證時(shí)段內(nèi)的監(jiān)測(cè)數(shù)據(jù)不應(yīng)超出預(yù)警界限。因此,本文將監(jiān)控模型的誤警率定義為驗(yàn)證時(shí)段內(nèi)被誤判為異常數(shù)據(jù)的數(shù)量與該時(shí)段監(jiān)測(cè)數(shù)據(jù)總數(shù)比值。

    模型評(píng)價(jià)指標(biāo)及誤警率FAR計(jì)算結(jié)果見(jiàn)表2。

    表2 模型評(píng)價(jià)指標(biāo)及誤警率統(tǒng)計(jì)結(jié)果Tab.2 Statistical results of model evaluation indenes and false alarm rates

    1)將表2中過(guò)度擬合系數(shù)OC、AIC值、BIC值和誤警率FAR分別繪制成散點(diǎn)圖,以O(shè)C=1.0和OC=1.2的直線將模型分為無(wú)過(guò)度擬合、低過(guò)度擬合、高過(guò)度擬合3種類型;以FAR=0和FAR=0.25的直線將模型分為無(wú)誤警、低誤警、高誤警3種類型,結(jié)果如圖6和7所示。

    圖6 監(jiān)控模型的過(guò)度擬合系數(shù)與誤警率散點(diǎn)圖Fig.6 Scatter plot of over-fitting coefficients and false alarm rates of monitoring model

    由圖6可知:過(guò)度擬合系數(shù)OC和誤警率FAR的散點(diǎn)圖呈較明顯的線性關(guān)系;當(dāng)無(wú)過(guò)度擬合現(xiàn)象時(shí)(OC≤0),模型誤警率 FAR均為0;當(dāng)過(guò)度擬合程度較低時(shí)(1.00<OC≤1.20),模型為低誤警(0<FAR≤0.25);當(dāng)過(guò)度擬合系數(shù)較高時(shí)(OC>1.20),模型為高誤警(FAR>0.25)。因此,模型誤警率FAR與過(guò)度擬合系數(shù)OC大致呈正相關(guān)關(guān)系。

    由圖7可知,AIC、BIC準(zhǔn)則同誤警率FAR無(wú)明顯關(guān)系。1#、2#、5#、6#等模型的AIC、BIC值相對(duì)較低,但誤警率FAR均高于50%,屬于高誤警類型;3#、7#、11#、12#等模型AIC、BIC值雖相對(duì)較高,但誤警率FAR均為0,屬于無(wú)誤警類型。因此,AIC、BIC準(zhǔn)則不能良好地反映模型誤警率FAR,與AIC、BIC值越小,模型效果越好的規(guī)則相矛盾。

    圖7 監(jiān)控模型的AIC和BIC計(jì)算值與誤警率散點(diǎn)圖Fig.7 Scatter plot of AIC and BIC values and false alarm rates of monitoring model

    因此,從模型誤警率FAR的角度出發(fā),R-OC準(zhǔn)則中的過(guò)度擬合系數(shù)OC比AIC、BIC準(zhǔn)則更能反映大壩變形監(jiān)控預(yù)警的真實(shí)情況。

    2)當(dāng)采用AIC、BIC準(zhǔn)則選取模型時(shí),根據(jù)AIC、BIC值最小原則,最優(yōu)模型為16#模型,其擬合段精度最高為0.974,但驗(yàn)證時(shí)段數(shù)據(jù)的預(yù)測(cè)效果并非最好,過(guò)度擬合系數(shù)為0.58。

    采用R-OC準(zhǔn)則選取模型時(shí),將復(fù)相關(guān)系數(shù)R和過(guò)度擬合系數(shù)OC繪制成散點(diǎn)圖,選出的較優(yōu)模型用紅線圈出,如圖8所示。

    圖8 監(jiān)控模型的R和OC指標(biāo)的散點(diǎn)圖Fig.8 Scatter plot of R and OC indexes of the monitoring model

    圖8中左上角區(qū)域(OC>1.0,R<0.9)中的模型既有欠擬合又存在過(guò)度擬合;右上角區(qū)域(OC>1.0,R≥0.9)中的模型無(wú)欠擬合但存在過(guò)度擬合;左下角區(qū)域(OC≤1.0,R<0.9)中的模型有欠擬合但無(wú)過(guò)度擬合;右下角區(qū)域(OC≤1.0,R≥0.9)中的模型既無(wú)欠擬合又無(wú)過(guò)度擬合。因此,滿足OC≤1.0、R≥0.9的模型為較優(yōu)模型。當(dāng)所有模型均不滿足此條件時(shí),可以從兩個(gè)方面進(jìn)行優(yōu)化:①改進(jìn)監(jiān)控模型的水壓、溫度和時(shí)效分量的函數(shù)形式;②采用深度學(xué)習(xí)算法[27]對(duì)歷史監(jiān)測(cè)數(shù)據(jù)進(jìn)行學(xué)習(xí),建立監(jiān)控模型。在LA14測(cè)點(diǎn)較優(yōu)模型9#、10#、13#、14#、15#、16#模型中,各模型的復(fù)相關(guān)系數(shù)R均在0.9以上,滿足模型精度的要求;對(duì)比各模型的過(guò)度擬合系數(shù)OC,并從小到大進(jìn)行排序,次序?yàn)?#、10#、13#、15#、16#、14#,最優(yōu)模型為9#模型。

    3)根據(jù)AIC、BIC準(zhǔn)則選擇模型時(shí),選擇最優(yōu)模型的準(zhǔn)確度受模型數(shù)量限制。由表2可知:①當(dāng)只構(gòu)建前12種模型時(shí),AIC、BIC準(zhǔn)則會(huì)將5#模型判定為最優(yōu)模型,雖然其復(fù)相關(guān)系數(shù)R為0.973,擬合精度最高,但其過(guò)度擬合系數(shù)為1.75,屬于高過(guò)度擬合類型,誤警率FAR高達(dá)70.79%,5#模型預(yù)警界限如圖9所示。但R-OC準(zhǔn)則會(huì)選擇9#模型,其復(fù)相關(guān)系數(shù)R為0.928,擬合精度滿足要求,同時(shí)過(guò)度擬合系數(shù)為0.24,無(wú)過(guò)度擬合現(xiàn)象,誤警率FAR為0。9#模型預(yù)警界限如圖10所示。②當(dāng)對(duì)18種模型進(jìn)行選擇時(shí),根據(jù)第2)點(diǎn)分析,AIC、BIC準(zhǔn)則將16#模型判定為最優(yōu)模型;此時(shí),R-OC準(zhǔn)則仍會(huì)選擇9#模型。

    圖9 基于AIC、BIC準(zhǔn)則識(shí)別的最優(yōu)模型(5#模型)的預(yù)警界限Fig.9 Early warning limits of optimal model (5# model)based on AIC and BIC criteria

    圖10 基于R-OC準(zhǔn)則識(shí)別的最優(yōu)模型(9#模型)的預(yù)警界限Fig.10 Early warning limits of an optimal model (9# model) based on R-OC criterion

    因此,對(duì)于不同數(shù)量的待選模型,R-OC準(zhǔn)則均能識(shí)別出擬合和預(yù)測(cè)精度都較高的模型,而AIC和BIC準(zhǔn)則識(shí)別出的最優(yōu)模型可能會(huì)存在較嚴(yán)重的過(guò)度擬合現(xiàn)象。

    4 結(jié) 論

    本文定義了能夠量化過(guò)度擬合程度的指標(biāo),即過(guò)度擬合系數(shù)OC。同時(shí),為滿足既不能欠擬合又不能過(guò)度擬合的模型識(shí)別原則,結(jié)合復(fù)相關(guān)系數(shù)R,提出了大壩變形監(jiān)控模型識(shí)別方法R-OC準(zhǔn)則。

    案例分析表明,過(guò)度擬合系數(shù)OC與誤警率FAR存在密切關(guān)系,當(dāng)過(guò)度擬合系數(shù)OC≤1時(shí),監(jiān)控模型誤警率FAR為0;當(dāng)過(guò)度擬合系數(shù)OC>1時(shí),監(jiān)控模型誤警率FAR與過(guò)度擬合系數(shù)OC呈正相關(guān)關(guān)系。R-OC模型識(shí)別準(zhǔn)則彌補(bǔ)了傳統(tǒng)模型識(shí)別方法AIC、BIC準(zhǔn)則不能定量比較和評(píng)價(jià)模型過(guò)度擬合程度的缺陷;通過(guò)R和OC兩個(gè)指標(biāo),既保證了模型的擬合精度,又限制了模型的過(guò)度擬合程度,合理選擇最優(yōu)模型,提高了模型識(shí)別的準(zhǔn)確度。同時(shí),對(duì)于不同數(shù)量的待選模型,R-OC準(zhǔn)則均能識(shí)別出擬合和預(yù)測(cè)精度都較高的模型。

    猜你喜歡
    監(jiān)測(cè)數(shù)據(jù)時(shí)效過(guò)度
    中藥煎煮前不宜過(guò)度泡洗
    過(guò)度減肥導(dǎo)致閉經(jīng)?
    GSM-R接口監(jiān)測(cè)數(shù)據(jù)精確地理化方法及應(yīng)用
    希望你沒(méi)在這里:對(duì)過(guò)度旅游的強(qiáng)烈抵制
    J75鋼的時(shí)效處理工藝
    一種新型耐熱合金GY200的長(zhǎng)期時(shí)效組織與性能
    上海金屬(2016年3期)2016-11-23 05:19:47
    環(huán)保執(zhí)法如何把握對(duì)違法建設(shè)項(xiàng)目的追責(zé)時(shí)效?
    過(guò)度加班,咋就停不下來(lái)?
    GPS異常監(jiān)測(cè)數(shù)據(jù)的關(guān)聯(lián)負(fù)選擇分步識(shí)別算法
    基于小波函數(shù)對(duì)GNSS監(jiān)測(cè)數(shù)據(jù)降噪的應(yīng)用研究
    黑人高潮一二区| 日韩 亚洲 欧美在线| 亚洲国产成人一精品久久久| 亚洲精华国产精华液的使用体验| 精品人妻偷拍中文字幕| 2018国产大陆天天弄谢| 亚洲精品日本国产第一区| av线在线观看网站| 少妇高潮的动态图| 精品国产乱码久久久久久小说| 91久久精品国产一区二区三区| 老司机亚洲免费影院| 十八禁高潮呻吟视频 | 我要看日韩黄色一级片| 日韩亚洲欧美综合| 亚洲丝袜综合中文字幕| 精品熟女少妇av免费看| 国产男人的电影天堂91| 国产精品久久久久久精品电影小说| 乱码一卡2卡4卡精品| 岛国毛片在线播放| 日韩精品有码人妻一区| 18禁在线无遮挡免费观看视频| 欧美成人精品欧美一级黄| 尾随美女入室| 观看av在线不卡| 一级毛片久久久久久久久女| 国产一级毛片在线| 女性生殖器流出的白浆| 极品人妻少妇av视频| 日日摸夜夜添夜夜爱| 一级av片app| av视频免费观看在线观看| 日本免费在线观看一区| 2022亚洲国产成人精品| 精品久久久久久电影网| 一级毛片 在线播放| 天美传媒精品一区二区| 老司机亚洲免费影院| 伊人久久国产一区二区| 内射极品少妇av片p| 在线天堂最新版资源| 国产有黄有色有爽视频| 久久精品国产a三级三级三级| 久久精品国产a三级三级三级| 久久久精品免费免费高清| 亚洲成色77777| 国产精品偷伦视频观看了| 啦啦啦在线观看免费高清www| 久久97久久精品| 亚洲欧美一区二区三区国产| 蜜桃久久精品国产亚洲av| 国产在线男女| 国产熟女欧美一区二区| 国产一区二区在线观看日韩| 国产精品欧美亚洲77777| 亚洲婷婷狠狠爱综合网| 少妇 在线观看| 狠狠精品人妻久久久久久综合| 另类亚洲欧美激情| 99视频精品全部免费 在线| 最黄视频免费看| 久久人人爽av亚洲精品天堂| 精品一区二区三区视频在线| 狠狠精品人妻久久久久久综合| 日韩精品免费视频一区二区三区 | 国产一区二区三区av在线| 精品一区在线观看国产| 亚洲婷婷狠狠爱综合网| 激情五月婷婷亚洲| 免费av不卡在线播放| 免费少妇av软件| 欧美97在线视频| av天堂中文字幕网| 少妇的逼好多水| 国产片特级美女逼逼视频| 男女无遮挡免费网站观看| 久久久久久久久久久久大奶| 国产精品女同一区二区软件| 亚洲美女搞黄在线观看| 亚洲第一区二区三区不卡| 精品午夜福利在线看| 女性被躁到高潮视频| 国产色婷婷99| 国产黄频视频在线观看| 久久6这里有精品| 黑人巨大精品欧美一区二区蜜桃 | 观看免费一级毛片| 亚洲精品亚洲一区二区| av国产精品久久久久影院| 特大巨黑吊av在线直播| 国产熟女欧美一区二区| 黄色一级大片看看| 国产av精品麻豆| 十分钟在线观看高清视频www | 性高湖久久久久久久久免费观看| 亚洲精品乱久久久久久| 国产日韩欧美视频二区| 伦理电影大哥的女人| 国产一区二区在线观看av| 亚洲欧美清纯卡通| 午夜视频国产福利| 国产一区二区在线观看日韩| 日韩亚洲欧美综合| 亚洲在久久综合| 国内揄拍国产精品人妻在线| 老司机影院成人| 成人国产av品久久久| 久久精品国产亚洲网站| 亚洲国产日韩一区二区| 自拍欧美九色日韩亚洲蝌蚪91 | 看非洲黑人一级黄片| 91精品国产国语对白视频| 亚洲精品乱码久久久v下载方式| 人体艺术视频欧美日本| 婷婷色麻豆天堂久久| 免费人成在线观看视频色| 欧美成人午夜免费资源| 两个人免费观看高清视频 | 免费大片18禁| 中文字幕精品免费在线观看视频 | 日韩大片免费观看网站| 国产美女午夜福利| 久久久国产精品麻豆| 亚洲av.av天堂| 男人爽女人下面视频在线观看| 80岁老熟妇乱子伦牲交| 国产午夜精品久久久久久一区二区三区| 日本av手机在线免费观看| 久久久久久人妻| 精品国产露脸久久av麻豆| 亚洲av福利一区| 97精品久久久久久久久久精品| 人妻一区二区av| 久久久久精品久久久久真实原创| 水蜜桃什么品种好| 色94色欧美一区二区| 精品午夜福利在线看| 少妇人妻精品综合一区二区| 老司机影院毛片| 国产国拍精品亚洲av在线观看| 秋霞在线观看毛片| 大陆偷拍与自拍| 亚洲av福利一区| 麻豆乱淫一区二区| 中文字幕人妻丝袜制服| 丰满迷人的少妇在线观看| 中文字幕免费在线视频6| 免费观看在线日韩| 欧美日韩视频高清一区二区三区二| 狠狠精品人妻久久久久久综合| 免费观看在线日韩| 免费观看性生交大片5| 日日啪夜夜爽| av在线观看视频网站免费| 一二三四中文在线观看免费高清| a级一级毛片免费在线观看| 国产精品无大码| 亚洲成人手机| 午夜福利视频精品| 国产乱人偷精品视频| 国产免费一级a男人的天堂| 偷拍熟女少妇极品色| 丝瓜视频免费看黄片| 亚洲欧美一区二区三区黑人 | 国产精品.久久久| 成人国产av品久久久| 欧美精品高潮呻吟av久久| 这个男人来自地球电影免费观看 | 老熟女久久久| a级毛片在线看网站| 亚洲精品久久午夜乱码| 精品视频人人做人人爽| 少妇人妻久久综合中文| 99热全是精品| av福利片在线观看| 人人澡人人妻人| 99久久人妻综合| 在线观看av片永久免费下载| 性色avwww在线观看| 99久久精品一区二区三区| 国产熟女午夜一区二区三区 | 大又大粗又爽又黄少妇毛片口| 亚洲真实伦在线观看| 18禁在线无遮挡免费观看视频| 一级毛片aaaaaa免费看小| 免费观看a级毛片全部| 国产精品久久久久成人av| 国产精品国产av在线观看| 人人妻人人澡人人看| 国产片特级美女逼逼视频| 国产一区二区三区综合在线观看 | 搡女人真爽免费视频火全软件| 日韩亚洲欧美综合| 十八禁网站网址无遮挡 | 亚洲精华国产精华液的使用体验| 国产精品国产三级国产av玫瑰| 亚洲精品成人av观看孕妇| 国产白丝娇喘喷水9色精品| 亚洲欧美成人精品一区二区| 街头女战士在线观看网站| 热re99久久国产66热| 99久久精品一区二区三区| 日韩精品有码人妻一区| 三级国产精品片| 人体艺术视频欧美日本| 日韩成人av中文字幕在线观看| 国产一区有黄有色的免费视频| 欧美精品一区二区免费开放| 黄色视频在线播放观看不卡| 欧美日韩综合久久久久久| 国国产精品蜜臀av免费| 国产日韩欧美亚洲二区| 美女主播在线视频| 最近的中文字幕免费完整| 成人亚洲精品一区在线观看| 视频中文字幕在线观看| 亚洲人成网站在线观看播放| 最近最新中文字幕免费大全7| 少妇熟女欧美另类| 国产成人精品久久久久久| 免费大片18禁| 另类亚洲欧美激情| 五月玫瑰六月丁香| 亚洲精品一区蜜桃| av免费观看日本| 久久久久精品性色| 欧美三级亚洲精品| 精品国产国语对白av| 另类亚洲欧美激情| 大片电影免费在线观看免费| 国产乱人偷精品视频| 自拍偷自拍亚洲精品老妇| 婷婷色综合www| 国产探花极品一区二区| 夫妻午夜视频| 欧美精品高潮呻吟av久久| 香蕉精品网在线| 亚洲av二区三区四区| 黑人猛操日本美女一级片| 伊人久久国产一区二区| 久久精品国产自在天天线| 国产黄色视频一区二区在线观看| 天堂中文最新版在线下载| 观看美女的网站| 免费少妇av软件| 26uuu在线亚洲综合色| 狂野欧美激情性bbbbbb| 寂寞人妻少妇视频99o| 日日爽夜夜爽网站| 高清视频免费观看一区二区| 一级毛片电影观看| 国产日韩欧美在线精品| 七月丁香在线播放| av线在线观看网站| 最新中文字幕久久久久| 亚洲国产欧美日韩在线播放 | 国产一级毛片在线| 日韩av不卡免费在线播放| 成人无遮挡网站| 国产毛片在线视频| 欧美另类一区| 成年女人在线观看亚洲视频| 国产成人免费观看mmmm| 男人狂女人下面高潮的视频| 日产精品乱码卡一卡2卡三| 在现免费观看毛片| 熟女电影av网| 国产免费一区二区三区四区乱码| 久久鲁丝午夜福利片| 激情五月婷婷亚洲| 秋霞伦理黄片| 国产精品国产三级国产av玫瑰| 女的被弄到高潮叫床怎么办| 日本与韩国留学比较| 成年人午夜在线观看视频| 欧美最新免费一区二区三区| 国产成人免费无遮挡视频| 人人妻人人看人人澡| 亚洲一级一片aⅴ在线观看| 美女视频免费永久观看网站| 亚洲精品,欧美精品| 国产淫片久久久久久久久| 亚洲成色77777| 日韩成人av中文字幕在线观看| 久久这里有精品视频免费| 中国三级夫妇交换| 乱码一卡2卡4卡精品| 欧美bdsm另类| 国产精品国产三级国产专区5o| 高清黄色对白视频在线免费看 | 22中文网久久字幕| 亚洲精品日本国产第一区| 国产男女超爽视频在线观看| av卡一久久| 国产毛片在线视频| 99九九线精品视频在线观看视频| 国产午夜精品一二区理论片| 免费av中文字幕在线| 国产日韩欧美在线精品| 18禁裸乳无遮挡动漫免费视频| 插阴视频在线观看视频| 成人亚洲精品一区在线观看| 日日爽夜夜爽网站| 欧美激情极品国产一区二区三区 | 一边亲一边摸免费视频| 中文字幕久久专区| 大码成人一级视频| 国产日韩欧美在线精品| 各种免费的搞黄视频| 欧美成人午夜免费资源| 精品久久久久久久久亚洲| 国产69精品久久久久777片| 男女边摸边吃奶| 成年人午夜在线观看视频| 精品国产露脸久久av麻豆| 卡戴珊不雅视频在线播放| 国产成人一区二区在线| 偷拍熟女少妇极品色| 久久毛片免费看一区二区三区| 日韩三级伦理在线观看| 在线观看一区二区三区激情| 亚洲精品乱久久久久久| 午夜影院在线不卡| 国产欧美另类精品又又久久亚洲欧美| 午夜影院在线不卡| 国产成人精品婷婷| 国产日韩欧美在线精品| 久久久国产一区二区| 精品亚洲成国产av| av免费在线看不卡| 国产高清三级在线| 亚洲精品乱久久久久久| 成人美女网站在线观看视频| 国产亚洲91精品色在线| 观看免费一级毛片| 精品午夜福利在线看| 少妇精品久久久久久久| 亚州av有码| 国产亚洲av片在线观看秒播厂| 视频区图区小说| 亚洲国产成人一精品久久久| 视频区图区小说| 亚洲国产成人一精品久久久| 亚洲欧洲国产日韩| 乱人伦中国视频| 国产欧美日韩综合在线一区二区 | 蜜臀久久99精品久久宅男| 亚洲欧美一区二区三区国产| 这个男人来自地球电影免费观看 | 欧美+日韩+精品| 亚洲国产精品一区三区| 精品少妇黑人巨大在线播放| 久久6这里有精品| 欧美日本中文国产一区发布| 一区二区三区四区激情视频| 国产黄片美女视频| 99久久综合免费| 伦理电影大哥的女人| 中文字幕av电影在线播放| 成年人免费黄色播放视频 | 另类亚洲欧美激情| 亚洲av电影在线观看一区二区三区| kizo精华| 日本av手机在线免费观看| 99国产精品免费福利视频| 偷拍熟女少妇极品色| 国产毛片在线视频| 亚洲欧美日韩卡通动漫| 欧美最新免费一区二区三区| 我要看黄色一级片免费的| 国产有黄有色有爽视频| 日韩亚洲欧美综合| 日韩av不卡免费在线播放| 亚洲欧美一区二区三区国产| 久久99精品国语久久久| 美女视频免费永久观看网站| 亚洲欧洲国产日韩| 久久久国产欧美日韩av| 少妇的逼好多水| 国产伦精品一区二区三区四那| 有码 亚洲区| 亚洲丝袜综合中文字幕| av天堂久久9| 亚洲无线观看免费| 久久6这里有精品| 五月开心婷婷网| 国产成人免费无遮挡视频| 国产白丝娇喘喷水9色精品| 毛片一级片免费看久久久久| 大陆偷拍与自拍| 春色校园在线视频观看| 我的老师免费观看完整版| 简卡轻食公司| 自拍偷自拍亚洲精品老妇| 亚洲av综合色区一区| 精品人妻偷拍中文字幕| 大香蕉97超碰在线| 免费观看在线日韩| 欧美精品亚洲一区二区| 国产精品一区二区性色av| 久久 成人 亚洲| 免费大片黄手机在线观看| 日日摸夜夜添夜夜爱| 人人妻人人澡人人爽人人夜夜| 王馨瑶露胸无遮挡在线观看| 亚洲经典国产精华液单| 18禁在线播放成人免费| 午夜视频国产福利| 色5月婷婷丁香| 少妇丰满av| 另类精品久久| 久久免费观看电影| 国内精品宾馆在线| 免费av不卡在线播放| 久久 成人 亚洲| 熟女av电影| 国产成人91sexporn| 超碰97精品在线观看| 观看av在线不卡| 亚洲精品乱码久久久v下载方式| 久久久久久久久大av| 久久久久久久亚洲中文字幕| 色哟哟·www| 免费看不卡的av| 偷拍熟女少妇极品色| 国产欧美亚洲国产| 最近手机中文字幕大全| 国产永久视频网站| 亚洲精品456在线播放app| 亚洲美女视频黄频| 国产免费又黄又爽又色| 亚洲欧洲精品一区二区精品久久久 | 国产亚洲精品久久久com| 国产成人精品久久久久久| 亚洲欧美一区二区三区国产| 欧美日韩视频精品一区| 乱人伦中国视频| av又黄又爽大尺度在线免费看| 久久婷婷青草| 色视频在线一区二区三区| 欧美日本中文国产一区发布| 一区二区三区乱码不卡18| 纵有疾风起免费观看全集完整版| 国产亚洲一区二区精品| 国产精品秋霞免费鲁丝片| 乱码一卡2卡4卡精品| 2018国产大陆天天弄谢| 在线精品无人区一区二区三| videossex国产| 亚洲欧美一区二区三区黑人 | 日本wwww免费看| 丝袜喷水一区| 日韩强制内射视频| 五月开心婷婷网| 欧美成人午夜免费资源| 亚洲欧美成人综合另类久久久| 丝袜喷水一区| 搡老乐熟女国产| 大香蕉久久网| 日韩亚洲欧美综合| 亚洲欧美精品专区久久| 高清在线视频一区二区三区| 五月玫瑰六月丁香| 亚洲性久久影院| 91精品伊人久久大香线蕉| av福利片在线观看| 水蜜桃什么品种好| av福利片在线| 亚洲国产毛片av蜜桃av| 狂野欧美激情性xxxx在线观看| 九草在线视频观看| 99久久精品国产国产毛片| 亚洲精品视频女| 91精品国产九色| av免费在线看不卡| 人人妻人人看人人澡| 高清不卡的av网站| 大话2 男鬼变身卡| 热re99久久精品国产66热6| 免费播放大片免费观看视频在线观看| 黄色欧美视频在线观看| 观看美女的网站| 美女脱内裤让男人舔精品视频| 一区二区三区免费毛片| 午夜福利在线观看免费完整高清在| 99热6这里只有精品| 欧美日韩精品成人综合77777| 亚洲美女搞黄在线观看| 男女边吃奶边做爰视频| 亚洲欧美日韩东京热| 中文字幕精品免费在线观看视频 | 我的老师免费观看完整版| 久久精品久久久久久久性| 极品人妻少妇av视频| 在线观看免费视频网站a站| 亚洲精品久久午夜乱码| 永久网站在线| 性色av一级| 极品人妻少妇av视频| 人妻 亚洲 视频| 最近中文字幕高清免费大全6| 中文字幕人妻丝袜制服| 建设人人有责人人尽责人人享有的| 免费播放大片免费观看视频在线观看| 丝袜喷水一区| 9色porny在线观看| 国产精品久久久久久久电影| 久久99热6这里只有精品| 国产成人91sexporn| 少妇被粗大的猛进出69影院 | 色婷婷av一区二区三区视频| 国产精品国产av在线观看| 少妇猛男粗大的猛烈进出视频| 午夜福利网站1000一区二区三区| 欧美成人精品欧美一级黄| 亚洲精品色激情综合| 中文资源天堂在线| 成人无遮挡网站| 婷婷色麻豆天堂久久| av线在线观看网站| 久久ye,这里只有精品| 色94色欧美一区二区| 亚洲欧美一区二区三区国产| 一区二区av电影网| 人妻 亚洲 视频| av卡一久久| 日韩一本色道免费dvd| 内地一区二区视频在线| 精品久久国产蜜桃| 十八禁网站网址无遮挡 | 在线观看免费日韩欧美大片 | av福利片在线观看| 99热全是精品| 一级毛片 在线播放| 中文资源天堂在线| 中文字幕精品免费在线观看视频 | 亚洲精品亚洲一区二区| 亚洲精品一区蜜桃| 青春草亚洲视频在线观看| 久久久久久久久久人人人人人人| 午夜激情福利司机影院| 九九久久精品国产亚洲av麻豆| 人妻 亚洲 视频| 六月丁香七月| 汤姆久久久久久久影院中文字幕| 伊人久久国产一区二区| 一级二级三级毛片免费看| 插逼视频在线观看| 婷婷色av中文字幕| 国产亚洲欧美精品永久| 成人免费观看视频高清| 一级黄片播放器| 99久久中文字幕三级久久日本| 国产毛片在线视频| 色婷婷久久久亚洲欧美| 国产欧美日韩综合在线一区二区 | 五月玫瑰六月丁香| 人妻制服诱惑在线中文字幕| 乱系列少妇在线播放| 日本wwww免费看| 国产精品不卡视频一区二区| 插阴视频在线观看视频| 国产 一区精品| 久久精品熟女亚洲av麻豆精品| 久久 成人 亚洲| 亚洲精品自拍成人| 在线观看人妻少妇| 日本爱情动作片www.在线观看| 成人毛片a级毛片在线播放| 精品99又大又爽又粗少妇毛片| 色哟哟·www| 成人亚洲欧美一区二区av| 国产精品国产三级国产专区5o| 亚洲精品乱久久久久久| 九草在线视频观看| 99热6这里只有精品| 日韩大片免费观看网站| 精品人妻熟女毛片av久久网站| 五月天丁香电影| 少妇人妻 视频| 中文资源天堂在线| 亚洲精品日本国产第一区| 亚洲欧洲国产日韩| 国产伦理片在线播放av一区| 老司机影院成人| 一区二区av电影网| 色5月婷婷丁香| 亚洲国产精品成人久久小说| 九九久久精品国产亚洲av麻豆| 国产毛片在线视频| 国产精品久久久久久精品电影小说| 日韩制服骚丝袜av| 亚洲无线观看免费| 国产淫片久久久久久久久| 国产一区二区在线观看av| 麻豆成人av视频| 欧美日韩av久久| 十分钟在线观看高清视频www | 久久影院123| 久久99一区二区三区| 久久精品国产a三级三级三级| 久久人人爽av亚洲精品天堂| 边亲边吃奶的免费视频| 十八禁高潮呻吟视频 | 天堂中文最新版在线下载| 中文字幕免费在线视频6| 国产日韩一区二区三区精品不卡 | 婷婷色综合大香蕉| 欧美高清成人免费视频www| 99久久精品国产国产毛片| 最新中文字幕久久久久| 国产毛片在线视频| 熟女人妻精品中文字幕| 久久久久久久久大av|