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

    基于隨機(jī)森林算法的CPTu土類識(shí)別模型研究及其在不同區(qū)域的應(yīng)用

    2023-12-29 01:53:28伍圣超王睿張建民
    關(guān)鍵詞:粉土砂土土層

    伍圣超 ,王睿 ,張建民

    (1. 清華大學(xué) 水沙科學(xué)與水利水電工程國家重點(diǎn)試驗(yàn)室,北京,100084;2. 城市軌道交通綠色與安全建造技術(shù)國家工程試驗(yàn)室,北京,100084;3. 清華大學(xué) 土木水利學(xué)院,北京,100084)

    CPTu 具有測試快捷、成本低、可靠性高以及測試數(shù)據(jù)沿深度連續(xù)的特點(diǎn),在巖土工程勘察中得到了廣泛應(yīng)用[1-3]。CPTu 無法獲取土樣,因此,不能采用傳統(tǒng)基于顆粒級(jí)配、塑性指數(shù)、液限的規(guī)范方法判別土類,但CPTu 在工程中應(yīng)用廣泛。已有研究提出一些利用CPTu數(shù)據(jù)來識(shí)別土類的方法,如:利用CPTu測試物理量組成二維圖表,以圖表不同區(qū)域代表不同土類。工程中常用的有ROBERTSON[2,4-5]提出的土類指數(shù)分類圖、SBT分類圖、SBTn 分類圖以及ESLAMI 等[6]提出的雙對數(shù)模型分類圖等。這類圖表法的缺點(diǎn)是其均基于土的力學(xué)行為分類,與基于土的物理特性的分類標(biāo)準(zhǔn)(如美國統(tǒng)一土質(zhì)分類標(biāo)準(zhǔn)和中國國家標(biāo)準(zhǔn)《巖土工程勘察規(guī)范》分類標(biāo)準(zhǔn))缺乏一致性,無法直接對應(yīng)土的工程分類。同時(shí),這些土分類圖表法常用的歸一化參數(shù)(歸一化錐尖阻力(Qt)、歸一化側(cè)壁摩阻力(Fr)、歸一化孔隙水壓力(Bq))對于淺層土體常出現(xiàn)偏大的結(jié)果,造成淺層土體土類識(shí)別出現(xiàn)誤差。一些學(xué)者提出了基于概率統(tǒng)計(jì)理論的土類識(shí)別方法[7-10],但這些方法多仍以SBT/SBTn圖表法的框架為基礎(chǔ)。蔡國軍等[11]利用聚類分析方法實(shí)現(xiàn)了基于CPTu 的土層界面劃分。劉松玉等[12]采用SBTn 圖表法建立了基于CPTu 的中國實(shí)用土分類方法,其分類結(jié)果符合中國標(biāo)準(zhǔn)的土的工程分類結(jié)果。實(shí)際上,土的工程分類和CPTu測試物理量之間的映射關(guān)系非常復(fù)雜,二維圖表無法充分反映這種映射關(guān)系,且無法充分利用CPTu數(shù)據(jù)的所有信息。機(jī)器學(xué)習(xí)方法適合解決這類復(fù)雜映射關(guān)系問題,且適用于多維輸入,能夠充分挖掘CPTu數(shù)據(jù)中關(guān)于土類的信息。應(yīng)用機(jī)器學(xué)習(xí)的方法解決巖土工程中的復(fù)雜映射關(guān)系問題受到了普遍的關(guān)注,如基于CPTu 預(yù)測黏性土應(yīng)力歷史[13]、液化判別[14-15]、建立黏性土不排水抗剪強(qiáng)度預(yù)測模型[16-17]等。但目前對于基于CPTu 的土分類問題,缺乏基于全球多地區(qū)數(shù)據(jù)庫的有效機(jī)器學(xué)習(xí)模型。

    針對已有研究中的不足,本文作者基于3個(gè)地區(qū)(新西蘭、奧地利、德國)的“CPTu+鉆孔”數(shù)據(jù)庫,建立一個(gè)具有良好泛化性能的土類識(shí)別機(jī)器學(xué)習(xí)模型,探討多地區(qū)廣泛適用的土類識(shí)別模型的可行性。該模型采用CPTu數(shù)據(jù)的統(tǒng)計(jì)特征作為輸入,使得模型擁有良好的泛化性能;以基于USCS 標(biāo)準(zhǔn)對標(biāo)土的工程分類為模型輸出。同時(shí),為了選擇最優(yōu)算法,研究對比隨機(jī)森林(RF)、支持向量機(jī)(SVM)、BP 神經(jīng)網(wǎng)絡(luò)(BPANN)和最近鄰(KNN) 4種機(jī)器學(xué)習(xí)算法對本問題的適用性;將建立的模型與工程中廣泛采用的SBTn圖表法性能進(jìn)行對比,以闡明新方法的優(yōu)越性,并通過4個(gè)工程實(shí)例論證該模型應(yīng)用于工程實(shí)踐的有效性。

    1 “CPTu+鉆孔”數(shù)據(jù)庫

    1.1 新西蘭巖土工程數(shù)據(jù)庫(NZGD)

    通過數(shù)據(jù)預(yù)處理[19],將CPTu 數(shù)據(jù)拆分成若干0.3 m 長、無數(shù)據(jù)突變、土類相同的數(shù)據(jù)段,單個(gè)數(shù)據(jù)段作為土類識(shí)別模型的一個(gè)樣本。數(shù)據(jù)預(yù)處理的作用如下:1) 刪除土層界面附近的CPTu 數(shù)據(jù),消除界面附近的噪聲數(shù)據(jù)對模型訓(xùn)練的影響;2) 消除數(shù)據(jù)突變(如CPTu 測試過程中碰到石塊、貝殼等)對模型訓(xùn)練的影響;3) 單個(gè)CPTu 數(shù)據(jù)段可包含16~31個(gè)數(shù)據(jù),因此,可采用這些數(shù)據(jù)的統(tǒng)計(jì)值作為土類識(shí)別模型的輸入,以提取更多特征訓(xùn)練模型,增強(qiáng)模型性能。訓(xùn)練集、驗(yàn)證集(基督城)、驗(yàn)證集(北島)中樣本數(shù)量如表1所示,共包含礫石(GW、GP)、砂土(SW、SP)、粉質(zhì)砂土(SM)、粉土(ML、MH) 4種土類。

    表1 NZGD數(shù)據(jù)庫中訓(xùn)練集和驗(yàn)證集樣本數(shù)量Table 1 Number of samples for training and validation sets in NZGD

    1.2 “ 奧地利+ 德國”數(shù)據(jù)庫(Premstaller Geotechnik Database)

    Premstaller Geotechnik 數(shù)據(jù)庫由OBERHOLLENZER 等[20]建立,包含由Premstaller Geotechnik ZT GmbH 公司在奧地利和德國開展的靜力觸探測試結(jié)果(CPT、CPTu、SCPT、SCPTu)和鉆孔測試結(jié)果。采用EN ISO 14688-1[21]對鉆孔中土進(jìn)行分類。該數(shù)據(jù)庫中,將測點(diǎn)距離在50 m 以內(nèi)的“CPTu+鉆孔”作為1對數(shù)據(jù),并認(rèn)為其具有相同的土層剖面,共有46對“CPTu+鉆孔”測點(diǎn)。采用與NZGD 數(shù)據(jù)庫相同的數(shù)據(jù)預(yù)處理方法提取各類土的樣本,樣本數(shù)量如表2所示。該數(shù)據(jù)庫包含礫石(Gr, sa, si→Gr, co)、砂土(Sa, si→Sa, gr, si)、粉土(Si, sa, cl→Si, sa, gr)、黏土(Cl, Si, fsa→Cl, gr)4 種土類。由于NZGD 數(shù)據(jù)庫中缺乏黏土樣本,Premstaller Geotechnik 數(shù)據(jù)庫能夠成為NZGD 數(shù)據(jù)庫的有效補(bǔ)充。

    表2 Premstaller Geotechnik數(shù)據(jù)庫中訓(xùn)練集和驗(yàn)證集樣本數(shù)量Table 2 Number of samples for training and validation sets in Premstaller Geotechnik database

    NZGD 數(shù)據(jù)庫和Premstaller Geotechnik 數(shù)據(jù)庫分別采用USCS和EN ISO 14688-1土分類標(biāo)準(zhǔn),這2種標(biāo)準(zhǔn)存在一定的差別。例如:對于砂土顆粒和粉土顆粒的界限,EN ISO 14688-1 標(biāo)準(zhǔn)為63 μm,USCS 標(biāo)準(zhǔn)為75 μm;對于礫石顆粒和砂土顆粒的界限,EN ISO 14688-1 標(biāo)準(zhǔn)為2 mm,USCS 標(biāo)準(zhǔn)為4.75 mm。EN ISO 14688-1 標(biāo)準(zhǔn)中,fSa(fine sand,細(xì)砂)的粒徑范圍是63~200 μm,但此部分fSa 在USCS 標(biāo)準(zhǔn)中會(huì)被歸類為粉土;EN ISO 14688-1 標(biāo)準(zhǔn)中fGr(fine gravel,細(xì)礫)的粒徑范圍是2~6.3 mm,但此部分fGr在USCS標(biāo)準(zhǔn)中會(huì)被歸類為砂土。為了使得Premstaller Geotechnik數(shù)據(jù)庫和NZGD數(shù)據(jù)庫能夠合并使用,需要對Premstaller Geotechnik數(shù)據(jù)庫中樣本進(jìn)行處理,因此,本研究刪除了Premstaller Geotechnik數(shù)據(jù)庫中以fSa和fGr為主要成分的土類樣本,以實(shí)現(xiàn)兩個(gè)數(shù)據(jù)庫的統(tǒng)一。

    2 機(jī)器學(xué)習(xí)土類識(shí)別模型

    2.1 模型建立

    建立機(jī)器學(xué)習(xí)模型需要考慮輸入選擇、輸出選擇、學(xué)習(xí)器選擇、參數(shù)優(yōu)化等問題。本研究將一段CPTu 數(shù)據(jù)(包含16~31 個(gè)數(shù)據(jù)點(diǎn))的統(tǒng)計(jì)值作為模型輸入,同時(shí)在模型輸入中增加了相應(yīng)深度處的豎向應(yīng)力和靜水壓力。因此,模型共有8個(gè)輸入特征,分別為修正后錐尖阻力qt的平均值、摩阻比Rf的平均值、孔隙水壓力σv0的平均值、初始豎向應(yīng)力u0的平均值、初始孔隙水壓力的平均值、歸一化錐尖阻力Qtn的鄰近偏差、歸一化側(cè)壁摩阻力Fr的鄰近偏差、歸一化孔隙水壓力Bq的鄰近偏差s*。其中,鄰近偏差的定義為:

    其中:xi為CPTu 數(shù)據(jù)段中第i個(gè)數(shù)據(jù)點(diǎn)(Qtn、Fr或Bq);n為CPTu 數(shù)據(jù)段中數(shù)據(jù)點(diǎn)總數(shù)。鄰近偏差表示CPTu數(shù)據(jù)段中數(shù)據(jù)的離散程度,是與土類相關(guān)的重要數(shù)據(jù)特征參數(shù)。

    該土類識(shí)別模型共有4 個(gè)輸出:礫石、砂土、粉土、黏土。SM根據(jù)細(xì)粒含量可劃分為砂土或粉土,該細(xì)粒含量閾值為37%[19]。采用總體精確度(ra)、單類土精確度(rs)和Kappa 系數(shù)(κ) 3 個(gè)指標(biāo)來衡量模型性能。其中,κ用來衡量真實(shí)值和模型預(yù)測值之間的一致性程度,其取值范圍為-1~1;ra、rs和κ越大表示模型性能越好,同時(shí),κ是評估模型對不平衡分類問題性能的指標(biāo)。ra、rs和κ的定義見式(2)~(4)。

    其中:ni為某一土類(礫石、砂土或粉土)分類正確的數(shù)量;Ni為相應(yīng)土類的樣本總數(shù);Po和Pe分別為觀測結(jié)果和預(yù)測結(jié)果的一致性概率。

    達(dá)沙替尼和伊馬替尼的不良反應(yīng)信號(hào)檢測研究…………………………………………………… 吳邦華等(20):2840

    2.2 不同機(jī)器學(xué)習(xí)方法性能對比

    隨機(jī)森林(RF)模型是若干決策樹學(xué)習(xí)器的集合,該集合中的每個(gè)決策樹學(xué)習(xí)器均可估計(jì)一個(gè)分類,而隨機(jī)森林綜合每個(gè)決策樹的分類決定最終的分類結(jié)果,隨機(jī)森林通過樣本選擇和特征選擇提高模型預(yù)測性能。支持向量機(jī)(SVM)基于訓(xùn)練集在樣本空間中找到一個(gè)劃分超平面,將不同類別的樣本區(qū)分開,當(dāng)超平面距離最近數(shù)據(jù)樣本的間隔越大時(shí),分類的魯棒性越好。BP 神經(jīng)網(wǎng)絡(luò)(BPANN)關(guān)鍵在于學(xué)習(xí)大量數(shù)據(jù)樣本中隱藏的輸入-輸出模式映射關(guān)系,學(xué)習(xí)規(guī)則采用最速下降法,通過反向傳播來不斷調(diào)整網(wǎng)絡(luò)的權(quán)值和閾值。最近鄰(KNN)原理是根據(jù)預(yù)測樣本與訓(xùn)練樣本的距離排序,并選取k個(gè)距離最近的訓(xùn)練樣本,根據(jù)這些樣本分類判斷最終分類結(jié)果。

    本文通過對比RF、SVM、BPANN 和KNN 這4種機(jī)器學(xué)習(xí)算法,以選擇最適用于本研究問題的算法。模型訓(xùn)練采用K折交叉驗(yàn)證(K=5),采用貝葉斯算法優(yōu)化模型參數(shù),4 種算法分別訓(xùn)練的4 個(gè)土類識(shí)別模型優(yōu)化后的參數(shù)見表3,各模型在訓(xùn)練集、驗(yàn)證集(基督城)、驗(yàn)證集(北島)、驗(yàn)證集(Premstaller Geotechnik)上的性能分別如表4~7 所示??梢姡焊鶕?jù)ra、rs和κ這3 個(gè)指標(biāo),RF 方法具有最優(yōu)的性能,其次是BPANN、KNN、SVM。對于訓(xùn)練集、驗(yàn)證集(基督城)、驗(yàn)證集(北島)、驗(yàn)證集(Premstaller Geotechnik),RF 的κ分別為0.98(接近KNN 的1.00)、0.84、0.67、0.72,在4 個(gè)模型中分別排名第1、第1、第3、第1。雖然對于驗(yàn)證集(北島),RF 的κ(0.67)比SVM(κ=0.81)和BPANN(κ=0.78)的低,但觀察rs指標(biāo)可以發(fā)現(xiàn),3個(gè)模型針對每類土的預(yù)測精度差別較小。

    表3 4個(gè)土類識(shí)別機(jī)器學(xué)習(xí)模型優(yōu)化后的參數(shù)Table 3 Optimized parameters of four kinds of machine learning soil classification models

    表4 4種機(jī)器學(xué)習(xí)方法性能對比(訓(xùn)練集)Table 4 Performance comparison of four kinds of machine learning methods(training set)

    表5 4種機(jī)器學(xué)習(xí)方法性能對比(驗(yàn)證集(基督城))Table 5 Performance comparison of four kinds of machine learning methods(validation set(Christchurch))

    表6 4種機(jī)器學(xué)習(xí)方法性能對比(驗(yàn)證集(北島))Table 6 Performance comparison of four kinds of machine learning methods(validation set(North Island))

    表7 4種機(jī)器學(xué)習(xí)方法性能對比(驗(yàn)證集(Premstaller Geotechnik))Table 7 Performance comparison of four kinds of machine learning methods(validation set(Premstaller Geotechnik))

    本研究中數(shù)據(jù)集屬于典型的不平衡數(shù)據(jù)集,砂土樣本數(shù)量是黏土的5~6倍、粉土的2倍、礫石的5~6倍。訓(xùn)練數(shù)據(jù)的不均衡會(huì)導(dǎo)致對于各類土的預(yù)測精度不均衡,而選擇合適的算法以適應(yīng)該不平衡數(shù)據(jù)集至關(guān)重要。對比上述4 種算法可以發(fā)現(xiàn):RF 對于4 種土類的識(shí)別效果較為均衡,其次是KNN和BPANN,SVM的識(shí)別效果非常不均衡。因此,RF 適用于本研究中的不平衡分類問題,而SVM則不適用。本研究最終采用RF作為模型訓(xùn)練的學(xué)習(xí)器。RF 通過訓(xùn)練樣本隨機(jī)采樣規(guī)則提高了算法針對不平衡數(shù)據(jù)集的適應(yīng)性。

    2.3 模型驗(yàn)證

    采用RF訓(xùn)練土類識(shí)別模型,并采用驗(yàn)證集(基督城)、 驗(yàn)證集( 北島)、 驗(yàn)證集(Premstaller Geotechnik)對其性能進(jìn)行驗(yàn)證,模型訓(xùn)練和驗(yàn)證的結(jié)果以混淆矩陣的形式展示,見表8~11。

    表8 土類識(shí)別模型在訓(xùn)練集上的混淆矩陣Table 8 Confusion matrix of soil classification model on training set

    從表8可見:土體識(shí)別模型對于4種土的識(shí)別率均達(dá)到了96%以上,κ為0.98。在驗(yàn)證集(基督城)上,對礫石和砂土的識(shí)別率達(dá)到了95%左右,對粉土的識(shí)別率為85%,κ為0.84(表9)。在驗(yàn)證集(北島)上,對砂土的識(shí)別率達(dá)到了89%左右,對粉土的識(shí)別率為80%,κ為0.67(表10)。在驗(yàn)證集(Premstaller Geotechnik)上,對礫石的識(shí)別率為78%,對砂土的識(shí)別率為75%,對粉土和黏土的識(shí)別率達(dá)到了85%左右,κ為0.72(表11)。說明模型在3個(gè)驗(yàn)證集上均取得滿意的效果。

    表9 土類識(shí)別模型在驗(yàn)證集(基督城)上的混淆矩陣Table 9 Confusion matrix of soil classification model on validation set(Christchurch)

    表10 土類識(shí)別模型在驗(yàn)證集(北島)上的混淆矩陣Table 10 Confusion matrix of soil classification model on validation set(North Island)

    表11 土類識(shí)別模型在驗(yàn)證集(Premstaller Geotechnik)上的混淆矩陣Table 11 Confusion matrix of soil classification model on validation set(Premstaller Geotechnik)

    ROC曲線和AUC被用來進(jìn)一步評估模型的性能。ROC曲線是在不同的閾值(某一類土的識(shí)別率)下,以縱坐標(biāo)為真正率(如粉土識(shí)別為粉土的概率),以橫坐標(biāo)為對應(yīng)的假正率(如非粉土識(shí)別為粉土的概率),對于單獨(dú)某一類土繪制而成。綜合ROC 曲線整體考慮了各類土的識(shí)別效果,是一個(gè)針對模型整體性能的評價(jià)指標(biāo)。理論上,最優(yōu)模型的ROC曲線是頂點(diǎn)在坐標(biāo)(0,1)處的兩段直線,而預(yù)測結(jié)果隨機(jī)分布的模型,其ROC 曲線則為通過原點(diǎn)的1∶1 直線。AUC 是指ROC 曲線和坐標(biāo)軸所包圍的面積,是一個(gè)評價(jià)土類識(shí)別模型性能的量化指標(biāo),AUC越高,表明模型性能越好。

    不同區(qū)域驗(yàn)證集數(shù)據(jù)的ROC 曲線和AUC 如圖1所示。對于驗(yàn)證集(基督城),其各類土的ROC曲線和綜合ROC 曲線接近最優(yōu)模型的ROC 曲線,AUC 達(dá)到0.98 左右,如圖1(a)所示,說明土類識(shí)別模型在驗(yàn)證集(基督城)上的性能表現(xiàn)優(yōu)異。對于驗(yàn)證集(北島),其各類土的ROC 曲線和綜合ROC曲線仍與最優(yōu)模型的ROC曲線較接近,AUC達(dá)到0.90左右,如圖1(b)所示,土類識(shí)別模型在驗(yàn)證集(北島)上的性能表現(xiàn)良好。驗(yàn)證集(Premstaller Geotechnik)數(shù)據(jù)來源于奧地利和德國,相比于驗(yàn)證集(北島),各類土的ROC 曲線和綜合ROC 曲線更接近最優(yōu)模型的ROC曲線,黏土和礫石的AUC在0.98 左右,粉土和砂土的AUC 為0.95 左右,綜合ROC 曲線AUC 為0.96,如圖1(c)所示,土類識(shí)別模型在驗(yàn)證集(Premstaller Geotechnik)上的性能表現(xiàn)同樣優(yōu)異。

    3 同SBTn圖表法對比

    為了證明本文模型的優(yōu)越性,將本文建立的土類識(shí)別模型與工程中應(yīng)用最為廣泛的SBTn圖表法[5]進(jìn)行對比。將SBTn 圖表法和本文模型在3 個(gè)不同區(qū)域(基督城(新西蘭)、北島(新西蘭)、奧地利+德國)的土類預(yù)測效果同時(shí)展示在Qtn-Fr圖中,見圖2~4。鉆孔土分類結(jié)果在Qtn-Fr圖中用不同的顏色表示,本文模型土類預(yù)測結(jié)果用不同的符號(hào)表示,因此,驗(yàn)證集中的每一個(gè)樣本在Qtn-Fr圖中表示為一個(gè)帶有顏色的符號(hào)。若某一樣本的所屬顏色和符號(hào)指向同一土類,則說明本文模型預(yù)測正確,若符號(hào)指向其他土類,則模型錯(cuò)誤預(yù)測為其他土類,例如,黃色上三角的樣本代表礫石正確預(yù)測為礫石,黃色圓形樣本代表礫石錯(cuò)誤預(yù)測為砂土。SBTn 圖表法的預(yù)測結(jié)果在Qtn-Fr圖中用數(shù)字1~9代表的不同區(qū)域表示,樣本點(diǎn)落在某區(qū)域則說明SBTn圖表法將該樣本預(yù)測為該區(qū)域數(shù)字代表的土類。

    圖2 機(jī)器學(xué)習(xí)土類識(shí)別模型和SBTn圖表法對比(基督城)Fig. 2 Comparation between machine learning soil classification model and SBTn method(Christchurch)

    從圖2~4 可以發(fā)現(xiàn),在驗(yàn)證數(shù)據(jù)所在地區(qū)(新西蘭、奧地利、德國),本文建立的土類識(shí)別模型相對于SBTn圖表法具有明顯的優(yōu)勢。對于新西蘭基督城地區(qū)(圖2),本文模型對于礫石和粉土的識(shí)別率高于SBTn圖表法,在SBTn圖表法中有1/3的礫石樣本落在了區(qū)域6(粉質(zhì)砂土-砂土),而本文模型將其正確識(shí)別為礫石;有1/4的粉土樣本落在了區(qū)域3、6、8、9(黏土-粉質(zhì)黏土、粉質(zhì)砂土-砂土、黏質(zhì)砂土-極硬砂土、極硬細(xì)砂),而本文模型將這些被SBTn 圖表法錯(cuò)分類的粉土樣本正確識(shí)別。對于新西蘭北島地區(qū)(圖3),本文模型對粉土的識(shí)別率相對于SBTn 圖表法有明顯提升,約1/3粉土樣本在SBTn 圖表中落在了區(qū)域3(黏土-粉質(zhì)黏土),本文模型均將這些粉土樣本正確識(shí)別;另有一些粉土樣本落在了區(qū)域9(極硬細(xì)砂),本文模型亦成功識(shí)別。對于奧地利和德國(圖4),本文模型同樣較大地提升了礫石和粉土的識(shí)別率,另外,該地區(qū)砂土樣本在SBTn圖表中約一半落在了區(qū)域5(砂質(zhì)粉土-粉質(zhì)砂土),而本文模型將其直接識(shí)別為砂土,減少了土類識(shí)別的不確定性。

    圖3 機(jī)器學(xué)習(xí)土類識(shí)別模型和SBTn圖表法對比(北島)Fig. 3 Comparation between machine learning soil classification model and SBTn method(North Island)

    圖4 機(jī)器學(xué)習(xí)土類識(shí)別模型和SBTn圖表法對比(奧地利+德國)Fig. 4 Comparation between machine learning soil classification model and SBTn method(Austria & Germany)

    4 工程應(yīng)用

    本文的土類識(shí)別模型是基于CPTu數(shù)據(jù)段的統(tǒng)計(jì)特征進(jìn)行預(yù)測,因此,在實(shí)際應(yīng)用于完整的CPTu 測試曲線時(shí),需要提前確定土層界面位置以將CPTu 測試曲線劃分成若干段,再基于各段CPTu 數(shù)據(jù)的統(tǒng)計(jì)特征預(yù)測土類。采用小波變換模極大值方法(WTMM)[22]確定土層界面。首先,通過WTMM方法確定土層數(shù)量和每層土的厚度,然后采用土類識(shí)別模型預(yù)測每層土的類別,最終重構(gòu)CPTu測點(diǎn)處的土層分布。本文分別以新西蘭和“奧地利+德國”區(qū)域?qū)嶋H工程中的CPTu 數(shù)據(jù)為例,闡述基于CPTu重構(gòu)土層分布的應(yīng)用流程,并以鄰近的鉆孔數(shù)據(jù)驗(yàn)證其效果。

    4.1 新西蘭

    對新西蘭地區(qū)2 個(gè)CPTu 測點(diǎn)(CPT_19118,CPT_6527,均為NZGD中測點(diǎn)代號(hào))進(jìn)行土層界面劃分和土類識(shí)別,重構(gòu)測點(diǎn)處的土層分布,并對比重構(gòu)的土層分布和鄰近鉆孔中土層分布,以說明方法的實(shí)際效果。2 個(gè)CPTu 測點(diǎn)處的土層重構(gòu)結(jié)果分別如圖5和圖6所示,可見重構(gòu)土層分布與鉆孔土層分布之間具有較好的一致性。

    圖5 CPT_19118的土層重構(gòu)結(jié)果Fig. 5 Soil stratification results of CPTu log CPT_19118

    圖6 CPT_6527的土層重構(gòu)結(jié)果Fig. 6 Soil stratification result of CPTu log CPT_6527

    以鉆孔中土層分布作為真實(shí)值,以重構(gòu)土層中與鉆孔土層一致的土段長度占總長度的比例作為成功重構(gòu)率,2 個(gè)CPTu 測點(diǎn)的成功重構(gòu)率分別為97%和95%。對于測點(diǎn)CPT_19118(圖5),土層界面識(shí)別誤差小于0.3 m,上層粉土、下層礫石的土層結(jié)構(gòu)被成功識(shí)別;對于測點(diǎn)CPT_6527(圖6),雖然鉆孔所得土層分布較為復(fù)雜,但本方法重構(gòu)的土層分布與鉆孔土層分布基本一致,粉質(zhì)砂土(SM)被識(shí)別為粉土或砂土,因其屬于過渡土類。

    4.2 奧地利+德國

    對“奧地利+德國”地區(qū)2個(gè)CPTu測點(diǎn)(CPTu_1172, CPTu_1143,均為Premstaller Geotechnik 數(shù)據(jù)庫中測點(diǎn)代號(hào))進(jìn)行土層界面劃分和土類識(shí)別,重構(gòu)測點(diǎn)處的土層分布,并對比重構(gòu)的土層分布和鄰近鉆孔中的土層分布。2 個(gè)CPTu 測點(diǎn)的土層重構(gòu)結(jié)果分別如圖7和圖8所示,可見重構(gòu)土層分布與鉆孔土層分布之間良好的一致性。

    圖7 CPTu_1172的土層重構(gòu)結(jié)果Fig. 7 Soil stratification results of CPTu log CPTu_1172

    圖8 CPTu_1143的土層重構(gòu)結(jié)果Fig. 8 Soil stratification results of CPTu log CPTu_1143

    2 個(gè)CPTu 測點(diǎn)的成功重構(gòu)率分別為93%和96%。對于測點(diǎn)CPTu_1172(圖7),砂土、粉土、黏土互層的土層結(jié)構(gòu)被成功識(shí)別,在9.0 m 深度處,觀察CPTu數(shù)據(jù)可以發(fā)現(xiàn),qt、fs曲線在該深度處發(fā)生了突變,說明在該處可能存在薄夾層。對于測點(diǎn)CPTu_1143(圖8),礫石-砂土-粉土的土層結(jié)構(gòu)被成功識(shí)別,界面位置的預(yù)測誤差小于0.3 m,在3.3 m深度處,一個(gè)薄粉土層同樣被識(shí)別。

    5 結(jié)論

    1) 在僅區(qū)分礫石、砂土、粉土、黏土四大類的情況下,利用隨機(jī)森林算法的強(qiáng)非線性映射能力,可基于跨地區(qū)“CPTu+鉆孔”數(shù)據(jù)庫建立一個(gè)多地區(qū)適用的土類識(shí)別模型。

    2) 基于隨機(jī)森林算法的CPTu土類識(shí)別模型在新西蘭、奧地利、德國具有良好的泛化性能。結(jié)合相應(yīng)的土層界面確定方法,能夠成功重構(gòu)測點(diǎn)處的土層分布,重構(gòu)土層分布與鄰近鉆孔中土層分布之間具有很好的一致性。

    3) 基于隨機(jī)森林算法的CPTu土類識(shí)別模型相比SBTn圖表法具備顯著優(yōu)越的性能,尤其在礫石和粉土的預(yù)測方面優(yōu)勢明顯。

    4) 對于本文中的土類識(shí)別問題,相比于SVM、BPANN、KNN 算法,RF 算法具有最優(yōu)良的性能。RF 對不平衡分類效果良好,對于各類土的識(shí)別效果較好,而SVM在不平衡分類問題上效果較差。

    猜你喜歡
    粉土砂土土層
    乳化瀝青穩(wěn)定粉土性能及其農(nóng)村道路鋪裝實(shí)踐
    石油瀝青(2023年5期)2023-12-08 08:35:16
    水泥-瀝青粉復(fù)合穩(wěn)定粉土的強(qiáng)度及滲透特性
    土釘噴錨在不同土層的支護(hù)應(yīng)用及效果分析
    飽和粉土三軸等p應(yīng)力路徑試驗(yàn)特性研究
    水道港口(2021年3期)2021-08-24 07:08:40
    飽和砂土地層輸水管道施工降水方案設(shè)計(jì)
    土層 村與人 下
    土層——伊當(dāng)灣志
    土層 沙與土 上
    龍之中華 龍之砂土——《蟠龍壺》創(chuàng)作談
    城市淺埋隧道穿越飽和砂土復(fù)合地層時(shí)適宜的施工工法
    欧美另类一区| 男人爽女人下面视频在线观看| 婷婷色综合大香蕉| 亚洲激情五月婷婷啪啪| 国产高清不卡午夜福利| 黄片小视频在线播放| 国产一区有黄有色的免费视频| www日本在线高清视频| 日本猛色少妇xxxxx猛交久久| 日韩制服骚丝袜av| 日韩人妻精品一区2区三区| 国产精品一区二区在线观看99| 无遮挡黄片免费观看| 在线观看www视频免费| 王馨瑶露胸无遮挡在线观看| 久久久精品94久久精品| 国产精品一国产av| 免费高清在线观看日韩| 国产午夜精品一二区理论片| 日韩制服丝袜自拍偷拍| 下体分泌物呈黄色| 美女扒开内裤让男人捅视频| xxxhd国产人妻xxx| 大型av网站在线播放| 日本五十路高清| 人妻 亚洲 视频| 天天躁日日躁夜夜躁夜夜| 日韩av免费高清视频| 99香蕉大伊视频| 嫁个100分男人电影在线观看 | 精品福利永久在线观看| videosex国产| a级片在线免费高清观看视频| 日本91视频免费播放| 亚洲精品中文字幕在线视频| 精品少妇一区二区三区视频日本电影| 啦啦啦在线免费观看视频4| av又黄又爽大尺度在线免费看| 精品一区二区三区四区五区乱码 | 国产成人欧美在线观看 | 国产成人精品无人区| 国产亚洲欧美在线一区二区| 1024香蕉在线观看| 欧美av亚洲av综合av国产av| 欧美精品高潮呻吟av久久| 婷婷成人精品国产| 韩国精品一区二区三区| 国产精品国产三级国产专区5o| 久久这里只有精品19| 性少妇av在线| 欧美黄色淫秽网站| 搡老乐熟女国产| 国产99久久九九免费精品| av国产久精品久网站免费入址| 高清欧美精品videossex| 真人做人爱边吃奶动态| 亚洲欧美一区二区三区久久| 一二三四在线观看免费中文在| 韩国精品一区二区三区| 赤兔流量卡办理| 看十八女毛片水多多多| 国产福利在线免费观看视频| 制服诱惑二区| 两性夫妻黄色片| 精品国产乱码久久久久久小说| 国产成人一区二区三区免费视频网站 | 老汉色∧v一级毛片| 女人爽到高潮嗷嗷叫在线视频| 我的亚洲天堂| 制服人妻中文乱码| 91老司机精品| 亚洲欧洲精品一区二区精品久久久| 久久99热这里只频精品6学生| 女人爽到高潮嗷嗷叫在线视频| 亚洲专区国产一区二区| 女人被躁到高潮嗷嗷叫费观| 国产日韩欧美亚洲二区| 日韩 欧美 亚洲 中文字幕| www.999成人在线观看| 免费看不卡的av| 人人妻人人澡人人爽人人夜夜| 国产伦人伦偷精品视频| 看十八女毛片水多多多| 搡老岳熟女国产| 日韩大片免费观看网站| 国产在线免费精品| a级片在线免费高清观看视频| 蜜桃在线观看..| 免费看不卡的av| 久久人人97超碰香蕉20202| 亚洲av综合色区一区| 亚洲精品成人av观看孕妇| 一级黄色大片毛片| 91九色精品人成在线观看| 亚洲国产最新在线播放| 一区二区三区四区激情视频| 国产人伦9x9x在线观看| 亚洲成国产人片在线观看| 免费不卡黄色视频| 大码成人一级视频| 成人国语在线视频| 操美女的视频在线观看| 观看av在线不卡| 青草久久国产| 亚洲激情五月婷婷啪啪| 色精品久久人妻99蜜桃| 欧美精品人与动牲交sv欧美| 久久99一区二区三区| 日韩av不卡免费在线播放| 亚洲精品自拍成人| 超色免费av| 日本av免费视频播放| 欧美人与性动交α欧美精品济南到| 99精品久久久久人妻精品| 亚洲av成人精品一二三区| 日韩一卡2卡3卡4卡2021年| 性少妇av在线| 少妇的丰满在线观看| 精品久久蜜臀av无| 欧美中文综合在线视频| 国产高清videossex| 国产av国产精品国产| 久久这里只有精品19| 成人黄色视频免费在线看| 国产色视频综合| 51午夜福利影视在线观看| 日韩欧美一区视频在线观看| 亚洲精品在线美女| 久久人人爽人人片av| 黄网站色视频无遮挡免费观看| 晚上一个人看的免费电影| av天堂在线播放| 久久久国产欧美日韩av| 一个人免费看片子| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲第一av免费看| 午夜福利乱码中文字幕| 亚洲国产精品999| 国产不卡av网站在线观看| 欧美激情极品国产一区二区三区| 夫妻性生交免费视频一级片| 免费观看人在逋| a级毛片在线看网站| 两性夫妻黄色片| 国产欧美日韩综合在线一区二区| 99热国产这里只有精品6| 十分钟在线观看高清视频www| 欧美大码av| 精品欧美一区二区三区在线| www.精华液| 中文字幕最新亚洲高清| 91老司机精品| av在线播放精品| 国语对白做爰xxxⅹ性视频网站| 国产成人a∨麻豆精品| 国产精品久久久久成人av| 久久久久网色| 中文乱码字字幕精品一区二区三区| 久久国产精品影院| 亚洲精品国产区一区二| 亚洲专区中文字幕在线| 老熟女久久久| 国产又爽黄色视频| 另类亚洲欧美激情| 国产有黄有色有爽视频| 精品免费久久久久久久清纯 | 亚洲精品国产av蜜桃| 黑人巨大精品欧美一区二区蜜桃| 亚洲国产毛片av蜜桃av| 国产精品久久久人人做人人爽| 久久精品国产亚洲av涩爱| 国产有黄有色有爽视频| 欧美精品一区二区大全| 嫁个100分男人电影在线观看 | 久久久久久免费高清国产稀缺| 嫩草影视91久久| 在线观看免费视频网站a站| 国产精品二区激情视频| 丝袜美腿诱惑在线| av在线app专区| 精品人妻在线不人妻| 亚洲成av片中文字幕在线观看| 叶爱在线成人免费视频播放| 欧美在线一区亚洲| 69精品国产乱码久久久| 一级毛片 在线播放| 中文字幕最新亚洲高清| 男男h啪啪无遮挡| 精品福利观看| 亚洲国产欧美日韩在线播放| videos熟女内射| 两个人免费观看高清视频| 大香蕉久久网| 大型av网站在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 精品一区二区三区四区五区乱码 | 国产1区2区3区精品| 女性被躁到高潮视频| 久久精品久久久久久久性| 极品少妇高潮喷水抽搐| 亚洲国产看品久久| 欧美老熟妇乱子伦牲交| 久久久国产精品麻豆| 大型av网站在线播放| 两性夫妻黄色片| 一本—道久久a久久精品蜜桃钙片| 日本猛色少妇xxxxx猛交久久| 看免费av毛片| 夜夜骑夜夜射夜夜干| 9191精品国产免费久久| 高清欧美精品videossex| 亚洲av电影在线进入| 黄色 视频免费看| 免费少妇av软件| 女人久久www免费人成看片| 亚洲美女黄色视频免费看| av一本久久久久| 亚洲免费av在线视频| 少妇人妻 视频| 亚洲精品国产av成人精品| 无遮挡黄片免费观看| 欧美日韩综合久久久久久| 国产成人91sexporn| 在线观看免费高清a一片| 久久久久久久大尺度免费视频| 三上悠亚av全集在线观看| 欧美国产精品一级二级三级| 亚洲国产精品一区三区| 丝袜美足系列| 一级黄片播放器| 国产xxxxx性猛交| 女人久久www免费人成看片| 后天国语完整版免费观看| 激情五月婷婷亚洲| 高清欧美精品videossex| 国产亚洲一区二区精品| 亚洲欧美色中文字幕在线| 五月开心婷婷网| 日韩大码丰满熟妇| av网站在线播放免费| 一级毛片女人18水好多 | 天天影视国产精品| 久久久久精品人妻al黑| 亚洲精品久久久久久婷婷小说| 日韩中文字幕视频在线看片| 脱女人内裤的视频| 丰满饥渴人妻一区二区三| 欧美在线黄色| 国产欧美亚洲国产| 99九九在线精品视频| 亚洲天堂av无毛| 午夜久久久在线观看| 99九九在线精品视频| 999久久久国产精品视频| 日本五十路高清| 国产在线视频一区二区| 亚洲一码二码三码区别大吗| 99精品久久久久人妻精品| av天堂久久9| 成年人午夜在线观看视频| 亚洲国产精品国产精品| 国产在线视频一区二区| 一区二区三区精品91| av网站在线播放免费| 丝袜美足系列| 蜜桃国产av成人99| 久久亚洲国产成人精品v| 久久久久久久久久久久大奶| 可以免费在线观看a视频的电影网站| 真人做人爱边吃奶动态| 精品少妇一区二区三区视频日本电影| 欧美在线黄色| 成人亚洲欧美一区二区av| 亚洲一卡2卡3卡4卡5卡精品中文| 久久ye,这里只有精品| 一边摸一边抽搐一进一出视频| 无遮挡黄片免费观看| 深夜精品福利| 男女边摸边吃奶| 国产男女内射视频| 免费在线观看视频国产中文字幕亚洲 | 高清黄色对白视频在线免费看| 欧美黄色淫秽网站| 啦啦啦视频在线资源免费观看| 亚洲精品国产区一区二| 99re6热这里在线精品视频| 男女边摸边吃奶| www.熟女人妻精品国产| 五月天丁香电影| 亚洲欧美色中文字幕在线| 国产男女超爽视频在线观看| 亚洲天堂av无毛| 久久国产精品人妻蜜桃| 又紧又爽又黄一区二区| 国产精品久久久av美女十八| 国产成人精品在线电影| 18禁黄网站禁片午夜丰满| av国产精品久久久久影院| 免费看av在线观看网站| 国产成人一区二区在线| 亚洲人成电影观看| 国产精品 欧美亚洲| 精品福利观看| 欧美黑人精品巨大| 久久精品国产亚洲av高清一级| 人人妻人人澡人人看| 亚洲色图综合在线观看| 精品国产超薄肉色丝袜足j| 久久天堂一区二区三区四区| 久久精品成人免费网站| 多毛熟女@视频| 欧美精品亚洲一区二区| 在线观看免费视频网站a站| 免费观看人在逋| 桃花免费在线播放| 午夜福利影视在线免费观看| 午夜福利视频在线观看免费| 一级毛片黄色毛片免费观看视频| 十分钟在线观看高清视频www| 亚洲精品国产av成人精品| 巨乳人妻的诱惑在线观看| 男女免费视频国产| 久久久久视频综合| 99国产精品一区二区三区| 免费少妇av软件| 在线看a的网站| 日韩av免费高清视频| 下体分泌物呈黄色| h视频一区二区三区| 99国产精品一区二区蜜桃av | 高清av免费在线| 丝袜脚勾引网站| 国产不卡av网站在线观看| 国产黄色免费在线视频| 国产日韩一区二区三区精品不卡| 日韩一区二区三区影片| 亚洲成人国产一区在线观看 | 又紧又爽又黄一区二区| 国产午夜精品一二区理论片| 国产黄频视频在线观看| 亚洲欧美日韩高清在线视频 | 美女主播在线视频| 中文字幕高清在线视频| 好男人电影高清在线观看| 成人亚洲欧美一区二区av| 亚洲欧洲日产国产| 亚洲精品第二区| 青春草视频在线免费观看| 最近最新中文字幕大全免费视频 | 丁香六月天网| 天堂中文最新版在线下载| 熟女少妇亚洲综合色aaa.| 在现免费观看毛片| 成年动漫av网址| 一级片免费观看大全| 成人18禁高潮啪啪吃奶动态图| 欧美亚洲日本最大视频资源| 九草在线视频观看| 久久久精品国产亚洲av高清涩受| 国产1区2区3区精品| 国产野战对白在线观看| 香蕉丝袜av| 中文精品一卡2卡3卡4更新| 欧美日韩精品网址| 亚洲av美国av| 欧美日韩一级在线毛片| 亚洲成人手机| 久久亚洲国产成人精品v| 久久99精品国语久久久| 亚洲情色 制服丝袜| 在线亚洲精品国产二区图片欧美| 波野结衣二区三区在线| 首页视频小说图片口味搜索 | 啦啦啦 在线观看视频| 免费在线观看黄色视频的| 色播在线永久视频| 精品一区二区三卡| 国产精品偷伦视频观看了| 国产99久久九九免费精品| 嫩草影视91久久| 亚洲一码二码三码区别大吗| 亚洲欧洲精品一区二区精品久久久| 麻豆乱淫一区二区| 深夜精品福利| 亚洲伊人久久精品综合| 少妇猛男粗大的猛烈进出视频| 自线自在国产av| 男女之事视频高清在线观看 | 一本综合久久免费| 午夜影院在线不卡| 永久免费av网站大全| 青青草视频在线视频观看| 男人舔女人的私密视频| 精品国产一区二区三区久久久樱花| av又黄又爽大尺度在线免费看| 建设人人有责人人尽责人人享有的| av有码第一页| 又粗又硬又长又爽又黄的视频| 韩国高清视频一区二区三区| 亚洲成色77777| 一级片'在线观看视频| 另类精品久久| 免费观看人在逋| 在线观看国产h片| 日韩制服丝袜自拍偷拍| 纵有疾风起免费观看全集完整版| 国产日韩一区二区三区精品不卡| 久久精品国产a三级三级三级| 免费女性裸体啪啪无遮挡网站| 青春草亚洲视频在线观看| 男人舔女人的私密视频| 老鸭窝网址在线观看| 亚洲 欧美一区二区三区| 丝袜喷水一区| 十八禁人妻一区二区| 亚洲欧美一区二区三区久久| av视频免费观看在线观看| 国产欧美日韩综合在线一区二区| 亚洲自偷自拍图片 自拍| 少妇被粗大的猛进出69影院| 黄色片一级片一级黄色片| 免费日韩欧美在线观看| 久久久欧美国产精品| 久久精品久久久久久噜噜老黄| 老司机亚洲免费影院| 亚洲精品一二三| e午夜精品久久久久久久| 国产亚洲av高清不卡| 夫妻性生交免费视频一级片| 欧美少妇被猛烈插入视频| www.自偷自拍.com| 交换朋友夫妻互换小说| 黑人欧美特级aaaaaa片| 国产野战对白在线观看| 亚洲精品日韩在线中文字幕| 精品高清国产在线一区| 国产av一区二区精品久久| 黄色片一级片一级黄色片| 黄色视频在线播放观看不卡| 在线天堂中文资源库| 2018国产大陆天天弄谢| 少妇粗大呻吟视频| 亚洲七黄色美女视频| 中国国产av一级| 国产精品麻豆人妻色哟哟久久| 超色免费av| 国产在视频线精品| 久久精品国产亚洲av涩爱| 国产亚洲欧美精品永久| 久久久久久人人人人人| 国产精品 国内视频| 亚洲中文日韩欧美视频| 久久毛片免费看一区二区三区| 久久久精品国产亚洲av高清涩受| 亚洲精品一区蜜桃| 一区二区三区乱码不卡18| 成人午夜精彩视频在线观看| 99热国产这里只有精品6| 手机成人av网站| 90打野战视频偷拍视频| √禁漫天堂资源中文www| 久久久久久久国产电影| 欧美精品人与动牲交sv欧美| 黄色毛片三级朝国网站| 国产淫语在线视频| 老司机午夜十八禁免费视频| 久久精品久久久久久久性| 国产女主播在线喷水免费视频网站| 久久国产亚洲av麻豆专区| 一本一本久久a久久精品综合妖精| 啦啦啦在线观看免费高清www| 热re99久久精品国产66热6| 国产有黄有色有爽视频| 日本av手机在线免费观看| 国产在线一区二区三区精| 国产黄色免费在线视频| 亚洲欧美清纯卡通| 国产精品二区激情视频| 大片电影免费在线观看免费| 在线亚洲精品国产二区图片欧美| 日韩人妻精品一区2区三区| 久久精品国产综合久久久| 尾随美女入室| 2018国产大陆天天弄谢| 国产日韩欧美亚洲二区| 十八禁人妻一区二区| 亚洲精品日本国产第一区| 女性生殖器流出的白浆| www.精华液| 日韩制服丝袜自拍偷拍| 十八禁人妻一区二区| 色婷婷久久久亚洲欧美| 欧美人与性动交α欧美精品济南到| 亚洲欧洲国产日韩| 亚洲专区国产一区二区| 久久国产精品大桥未久av| 国产一区二区 视频在线| 精品一品国产午夜福利视频| 国产欧美日韩综合在线一区二区| 丝袜脚勾引网站| 精品国产一区二区三区久久久樱花| 人人妻人人添人人爽欧美一区卜| 久久毛片免费看一区二区三区| 在线av久久热| 天天躁夜夜躁狠狠躁躁| 大香蕉久久成人网| 97人妻天天添夜夜摸| 免费观看人在逋| 精品少妇内射三级| 一二三四社区在线视频社区8| 美女扒开内裤让男人捅视频| av在线app专区| 国产成人av教育| 啦啦啦啦在线视频资源| 免费看十八禁软件| 日本欧美国产在线视频| 夫妻性生交免费视频一级片| 精品少妇久久久久久888优播| 久久久欧美国产精品| 国产精品免费大片| 搡老岳熟女国产| 国产精品秋霞免费鲁丝片| 老汉色∧v一级毛片| 母亲3免费完整高清在线观看| 亚洲国产最新在线播放| 国产一区二区三区综合在线观看| 9色porny在线观看| 欧美黑人精品巨大| 色综合欧美亚洲国产小说| 无遮挡黄片免费观看| 精品国产一区二区久久| 国产高清视频在线播放一区 | 少妇裸体淫交视频免费看高清 | 国产高清videossex| 国产免费又黄又爽又色| 最新在线观看一区二区三区 | 国产一区二区三区综合在线观看| 欧美另类一区| 精品一区二区三卡| 成年女人毛片免费观看观看9 | 下体分泌物呈黄色| 精品国产乱码久久久久久男人| 高清不卡的av网站| 亚洲久久久国产精品| 国产一区二区 视频在线| av视频免费观看在线观看| 免费久久久久久久精品成人欧美视频| 天天躁狠狠躁夜夜躁狠狠躁| 七月丁香在线播放| 亚洲av在线观看美女高潮| 99热国产这里只有精品6| 晚上一个人看的免费电影| 欧美日韩视频精品一区| 韩国高清视频一区二区三区| 亚洲精品久久久久久婷婷小说| 欧美日韩av久久| 大话2 男鬼变身卡| 亚洲自偷自拍图片 自拍| 一级毛片黄色毛片免费观看视频| 国产免费视频播放在线视频| 啦啦啦 在线观看视频| 一区二区三区四区激情视频| 亚洲第一青青草原| 亚洲欧洲精品一区二区精品久久久| 午夜福利一区二区在线看| 菩萨蛮人人尽说江南好唐韦庄| 国产黄色视频一区二区在线观看| 亚洲欧美精品自产自拍| 丁香六月欧美| 欧美黄色片欧美黄色片| 满18在线观看网站| 精品久久蜜臀av无| 亚洲欧美一区二区三区久久| 纯流量卡能插随身wifi吗| 国产又色又爽无遮挡免| videos熟女内射| 亚洲久久久国产精品| 99国产综合亚洲精品| 少妇的丰满在线观看| 国产野战对白在线观看| 亚洲国产av新网站| 亚洲精品第二区| 欧美亚洲 丝袜 人妻 在线| 一区二区三区激情视频| 大香蕉久久成人网| 国产精品久久久久久人妻精品电影 | 国产精品 欧美亚洲| 久热这里只有精品99| 男女高潮啪啪啪动态图| 国产91精品成人一区二区三区 | 欧美黄色片欧美黄色片| 50天的宝宝边吃奶边哭怎么回事| 精品少妇黑人巨大在线播放| 啦啦啦中文免费视频观看日本| 19禁男女啪啪无遮挡网站| 丁香六月欧美| 成人亚洲欧美一区二区av| 夜夜骑夜夜射夜夜干| 男女免费视频国产| 水蜜桃什么品种好| 免费观看av网站的网址| 男人爽女人下面视频在线观看| 永久免费av网站大全| 侵犯人妻中文字幕一二三四区| 久久性视频一级片| 欧美+亚洲+日韩+国产| 操美女的视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产欧美在线一区| 波多野结衣av一区二区av| 狂野欧美激情性xxxx| 欧美在线黄色| 伦理电影免费视频| 亚洲成av片中文字幕在线观看| 丁香六月欧美|