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

    基于層次分析法的德州地面沉降易發(fā)區(qū)的可拓學(xué)劃分應(yīng)用

    2016-08-05 10:05:32張永偉李紅霞商婷婷賈超梁浩楊亞賓
    山東國(guó)土資源 2016年7期
    關(guān)鍵詞:評(píng)價(jià)

    張永偉,李紅霞,商婷婷,賈超,梁浩,楊亞賓

    (1.山東省地質(zhì)環(huán)境監(jiān)測(cè)總站,山東 濟(jì)南 250014;2.山東大學(xué),山東 濟(jì)南 250061;3.濟(jì)南市育賢小學(xué),山東 濟(jì)南 250024;4.山東省魯北地質(zhì)工程勘察院,山東 德州 253015)

    ?

    技術(shù)方法

    基于層次分析法的德州地面沉降易發(fā)區(qū)的可拓學(xué)劃分應(yīng)用

    張永偉1,2,李紅霞3,商婷婷1,賈超2,梁浩1,楊亞賓4

    (1.山東省地質(zhì)環(huán)境監(jiān)測(cè)總站,山東 濟(jì)南250014;2.山東大學(xué),山東 濟(jì)南250061;3.濟(jì)南市育賢小學(xué),山東 濟(jì)南250024;4.山東省魯北地質(zhì)工程勘察院,山東 德州253015)

    摘要:基于物元模型的可拓學(xué)理論分析,提出華北平原德州地面沉降易發(fā)區(qū)的可拓學(xué)劃分方法。利用地面沉降易發(fā)性分級(jí)、分類(lèi)標(biāo)準(zhǔn)和影響因子分析,構(gòu)造出經(jīng)典域物元和節(jié)域物元,應(yīng)用物元和可拓集合中的關(guān)聯(lián)函數(shù),建立了易發(fā)性等級(jí)綜合評(píng)判的可拓評(píng)價(jià)模型,通過(guò)基于層次分析法的可拓學(xué)評(píng)價(jià)分析,采用正方形等間距(2km×2km)剖分了135個(gè)網(wǎng)格,再將剖分單元轉(zhuǎn)化成面元,經(jīng)過(guò)等值差分,得到了地面沉降易發(fā)區(qū)劃結(jié)果。計(jì)算結(jié)果顯示,可拓學(xué)方法能夠?qū)崿F(xiàn)定量化、多因子評(píng)價(jià)地面沉降易發(fā)性等級(jí)問(wèn)題,從而科學(xué)合理地指導(dǎo)地面沉降綜合分區(qū)防治。

    關(guān)鍵詞:地面沉降;可拓學(xué);物元模型;易發(fā)區(qū)劃;層次分析法

    引文格式:張永偉,李紅霞,賈超,等.基于層次分析法的德州地面沉降易發(fā)區(qū)的可拓學(xué)劃分應(yīng)用[J].山東國(guó)土資源,2016,32(7):58-64.ZHANG Yongwei,LI Hongxia,SHANG Tingting,etc. Application of Extenincs Theory for Partition of Land Subsidence Susceptible Area of Dezhou City Based on Analytic Hierarchy Process[J].Shandong Land and Resources, 2016,32(7):58-64.

    地面沉降是一種致使地面高程降低的緩變性地質(zhì)災(zāi)害[1]?;趯哟畏治龅目赏貙W(xué)是較為科學(xué)的定量化評(píng)價(jià)方法,采用這種方法對(duì)華北平原地面沉降最為嚴(yán)重的地區(qū)之一德州進(jìn)行易發(fā)區(qū)劃分,對(duì)于該區(qū)的地面沉降綜合防治具有一定的指導(dǎo)作用。

    1研究區(qū)概況

    德州市位于黃河下游沖積平原區(qū),1965年深層地下水開(kāi)采以來(lái),深層地下水位降落漏斗迅速在縱、橫向上拓展,漏斗中心水位埋深已超過(guò)140m。20世紀(jì)80年代,德州市范圍內(nèi)產(chǎn)生了不同程度的地面沉降。德州是華北平原地面沉降出現(xiàn)最早、影響最大和最為嚴(yán)重的地區(qū)之一,已與河北省形成連片沉降區(qū)[2]。

    德州市地面沉降各監(jiān)測(cè)點(diǎn)以每年16.9~62.5mm(1991—2010年)的速率下沉。沉降中心位于城區(qū)的國(guó)棉一廠院內(nèi)(D62),19年累計(jì)沉降量為-1186.9mm,年均沉降量為-62.5mm。城區(qū)邊緣的木材公司院內(nèi)(D9)累計(jì)沉降量(19年)為-636.9mm,年均沉降量(19年平均值)-33.5mm。外圍(YD1)累計(jì)沉降量為-321.7mm,年均沉降量(19年平均值)-16.9mm(表1)。由此看出,監(jiān)測(cè)區(qū)各觀測(cè)點(diǎn)存在明顯的不均勻沉降。隨著觀測(cè)點(diǎn)不均勻沉降的發(fā)生,地面點(diǎn)高差會(huì)越來(lái)越大[3]。

    2可拓學(xué)模型

    2.1可拓學(xué)簡(jiǎn)介

    可拓學(xué)(extenics)是一門(mén)橫斷學(xué)科和交叉學(xué)科,由廣東工業(yè)大學(xué)蔡文為首的中國(guó)學(xué)者創(chuàng)立,是研究事物拓展的可能性和開(kāi)拓創(chuàng)新的規(guī)律與方法[4]。目前,可拓學(xué)應(yīng)用于工業(yè)[5]、醫(yī)學(xué)、農(nóng)業(yè)、土地等級(jí)評(píng)價(jià)等領(lǐng)域,在地質(zhì)學(xué)上應(yīng)用如圍巖質(zhì)量分類(lèi)[6]、泥石流危險(xiǎn)度評(píng)價(jià)、地質(zhì)災(zāi)害脆弱性評(píng)價(jià)和損毀山體影響度評(píng)價(jià)等[7]。

    假設(shè)事物的等級(jí)為N,它關(guān)于特征c的量值為v,以此三元組作為表述事物的基本元,為物元R構(gòu)成的三要素,R=(N,c,v)。

    可拓學(xué)應(yīng)用的基本思路為:將評(píng)價(jià)目標(biāo)根據(jù)數(shù)據(jù)值進(jìn)行等級(jí)劃分,構(gòu)建可拓模型,再將評(píng)價(jià)指標(biāo)依次代入各等級(jí)的集合中進(jìn)行多指標(biāo)評(píng)定,根據(jù)評(píng)定結(jié)果與各等級(jí)集合的關(guān)聯(lián)度大小進(jìn)行比較來(lái)確定評(píng)定等級(jí)。

    表1 德州市1991—2010年各監(jiān)測(cè)點(diǎn)高程及沉降量統(tǒng)計(jì)

    2.2可拓學(xué)模型建立

    2.2.1確定經(jīng)典域

    (1)

    式中:N0j表示所劃分的第j(j=1,2,3,…,m)個(gè)等級(jí),ci(i=1,2,3,…,n)表示易發(fā)性等級(jí)N0j的因素,Voji=﹤aoji,boji﹥?yōu)镹0j關(guān)于因素ci所確定的量值范圍——經(jīng)典域。

    2.2.2確定節(jié)域

    (2)

    式中:P為等級(jí)的全體Vpi=﹤api,bpi﹥,為P關(guān)于因素i所取量值的范圍,即P的節(jié)域。

    2.2.3確定待評(píng)物元

    (3)

    式中:p為待評(píng)價(jià)目標(biāo),ci為影響易發(fā)等級(jí)的因素,vi為p關(guān)于因素ci的量值。

    2.2.4確定評(píng)價(jià)指標(biāo)關(guān)于各等級(jí)的關(guān)聯(lián)度

    各單項(xiàng)評(píng)價(jià)指標(biāo)Vi關(guān)于各等級(jí)j的關(guān)聯(lián)度為

    (4)

    2.2.5待評(píng)物元關(guān)于各類(lèi)別等級(jí)的關(guān)聯(lián)度

    待評(píng)目標(biāo)p關(guān)于評(píng)價(jià)等級(jí)j的關(guān)聯(lián)度為:

    (5)

    2.2.6可拓評(píng)判等級(jí)

    (6)

    (7)

    式中:j*為P的級(jí)別變量特征值,從j*中可看出偏向另一類(lèi)的程度。

    3地面沉降易發(fā)性等級(jí)可拓評(píng)價(jià)

    3.1評(píng)價(jià)體系的建立

    3.1.1評(píng)價(jià)因子的確定

    在綜合分析地面沉降影響因素的基礎(chǔ)上,參考有關(guān)研究成果[8-10],選擇地面沉降狀況因子和地質(zhì)環(huán)境條件因子(包括地下水開(kāi)采狀況)作為評(píng)價(jià)因子。地面沉降狀況因子包括累計(jì)沉降量、沉降速率和預(yù)測(cè)沉降量[11-14]。地質(zhì)環(huán)境條件因子包括水位允許降深、新生界厚度、軟土層厚度和地面高程。

    3.1.2評(píng)價(jià)指標(biāo)體系

    (1)體系的層次構(gòu)成。評(píng)價(jià)體系由評(píng)價(jià)目標(biāo)層、影響因子層(準(zhǔn)則層)、要素屬性層和要素指標(biāo)層構(gòu)成見(jiàn)表2。

    (2)評(píng)價(jià)指標(biāo)標(biāo)準(zhǔn):易發(fā)程度高為Ⅰ級(jí);易發(fā)程度中等為Ⅱ級(jí);易發(fā)程度低為Ⅲ級(jí)。

    3.2確定評(píng)價(jià)區(qū)并制作基礎(chǔ)網(wǎng)格

    該次評(píng)價(jià)范圍確定為德州市德城區(qū),面積為539km2。該次評(píng)價(jià)采用正方形等間距網(wǎng)格進(jìn)行剖分,網(wǎng)格大小為2km×2km,共剖分為135個(gè)網(wǎng)格。然后將剖分單元轉(zhuǎn)化成面元,并對(duì)每一剖分單元進(jìn)行連續(xù)編號(hào),將單元編號(hào)作為各面元的“ID”屬性賦值。將單元的各種指標(biāo)值分別賦予各單元的屬性值,建立單元屬性庫(kù)。

    3.3單元數(shù)據(jù)的獲取

    根據(jù)野外實(shí)測(cè)值、觀測(cè)值、地形圖、地質(zhì)圖量算值等方式獲取單元數(shù)據(jù)。

    3.4權(quán)值的確定

    權(quán)值的確定方法采用層次分析法(Analytic Hierarchy Process,簡(jiǎn)稱(chēng)AHP)[15]。層次分析定權(quán)法判斷矩陣標(biāo)度分級(jí)依據(jù)表3。權(quán)值計(jì)算采取自上而下(目標(biāo)層A—因子層B—屬性層C)的方法。

    表2 評(píng)價(jià)指標(biāo)體系與指標(biāo)等級(jí)劃分

    表3 層次分析定權(quán)法的判斷矩陣標(biāo)度分級(jí)及其含義

    3.4.1構(gòu)建目標(biāo)層A的綜合判斷矩陣

    為實(shí)現(xiàn)B層對(duì)目標(biāo)層A層的描述,對(duì)B層指標(biāo)進(jìn)行相對(duì)重要性?xún)蓛杀容^,確定地面沉降狀況(B1)和地質(zhì)環(huán)境條件(B2)的組合群對(duì)總目標(biāo)的權(quán)重關(guān)系。目標(biāo)層A的判斷矩陣為:A={bij│i,j=1-2}2×2在上述工作的基礎(chǔ)上,構(gòu)建的地面沉降危害性評(píng)價(jià)權(quán)重綜合判斷矩陣見(jiàn)表4。

    表4 地面沉降危害性評(píng)價(jià)權(quán)重綜合判斷矩陣

    3.4.2構(gòu)建2個(gè)因子層B的判斷矩陣

    分別根據(jù)因子層(B1、B2)的要求構(gòu)建B層和C層的判斷矩陣(表5、表6)。

    表5 地面沉降狀況因子相對(duì)重要性判斷矩陣

    表6 地質(zhì)環(huán)境狀況因子相對(duì)重要性判斷矩陣

    3.4.3層次排序及求解權(quán)向量

    根據(jù)判斷矩陣,利用線(xiàn)性代數(shù),求出矩陣的最大特征根所對(duì)應(yīng)的特征向量。所求特征向量即為各評(píng)價(jià)因素的重要性排序,歸一化后,也就是權(quán)數(shù)分配。

    對(duì)于一個(gè)正向量W=(W1,W2,…,Wn)T,其標(biāo)準(zhǔn)化向量W為:

    (8)

    式中:Ws(Ws1,Ws2,…,Wsn)T為同一層次相應(yīng)因子對(duì)于上一層次某個(gè)指標(biāo)相對(duì)重要性的排序權(quán)值。

    目前常見(jiàn)的計(jì)算因子排序權(quán)重向量的方法主要有方根法、和積法、特征根法、最小二乘法和對(duì)數(shù)最小二乘法等。鑒于方根法計(jì)算方便,該次采用方根法近似求解,并用Excel進(jìn)行計(jì)算[16]。計(jì)算結(jié)果見(jiàn)表7。

    表7 矩陣權(quán)值計(jì)算成果

    3.4.4判斷矩陣一致性檢驗(yàn)

    單層排序權(quán)值是否合理,需進(jìn)行一致性檢驗(yàn),公式為:

    (9)

    (10)

    式中:CR為一致性比例;RI為平均一致性指標(biāo),其數(shù)據(jù)如表8所示。

    CI為一致性指標(biāo)。當(dāng)CR<0.1時(shí),達(dá)到一致性;否則,需重新調(diào)整。判斷矩陣一致性檢驗(yàn)結(jié)果如表9所示,CR值均<0.1,判斷矩陣均具有滿(mǎn)意的一致性。

    表8 層次分析法的平均隨機(jī)一致性指標(biāo)值

    表9 矩陣一致性檢驗(yàn)結(jié)果

    3.4.5確定層次總排序權(quán)值

    目標(biāo)層A包含2個(gè)因子B1和B2,其權(quán)值分別為a1和a2;B層包含7個(gè)因子C1,C2,C3,...,Ck,...,C7;它們對(duì)因子Bj的層次排序權(quán)值分別為b1j,b2j,...,bnj(當(dāng)Ck與Bj無(wú)關(guān)時(shí),bnj=0)。C層總排序權(quán)值為二層權(quán)值乘積代數(shù)和。計(jì)算公式為:

    (11)

    計(jì)算結(jié)果見(jiàn)表10。

    表10 評(píng)價(jià)體系權(quán)值計(jì)算結(jié)果

    3.5可拓評(píng)價(jià)

    根據(jù)式(6)計(jì)算出每個(gè)單元網(wǎng)格的K值,確定每個(gè)單位網(wǎng)格的可拓評(píng)判等級(jí),分高、中、低三級(jí)評(píng)價(jià)其易發(fā)程度,最后根據(jù)所有單元網(wǎng)格可拓等級(jí),并綜合考慮地質(zhì)環(huán)境條件和地面沉降狀況綜合確定地面沉降高、中、低易發(fā)區(qū)。評(píng)價(jià)過(guò)程簡(jiǎn)述如下。

    3.5.1物元的構(gòu)造

    將易發(fā)程度等級(jí)Ⅰ,Ⅱ,Ⅲ分別記為N01,N02,N03;可以構(gòu)造出易發(fā)程度各等級(jí)的經(jīng)典物元和節(jié)域物元。

    3.5.2確定待評(píng)物元

    Ri為第i個(gè)待評(píng)的單元格,vij為第i個(gè)待評(píng)單元格第j個(gè)預(yù)測(cè)指標(biāo)實(shí)際值。按公式(4)計(jì)算其評(píng)價(jià)指標(biāo)C1關(guān)于各易發(fā)程度等級(jí)的關(guān)聯(lián)度。按照公式(5)計(jì)算出待評(píng)物元關(guān)于各類(lèi)別等級(jí)的關(guān)聯(lián)度,評(píng)價(jià)單元(沉降中心D62監(jiān)測(cè)點(diǎn)所在單元格)計(jì)算結(jié)果見(jiàn)表11。

    表11 關(guān)聯(lián)度計(jì)算結(jié)果

    表中K01(P)=0.184;K02(P)=-0.579;K03(P)=-0.727,經(jīng)比較分析,D62所在單元格易發(fā)程度屬于Ⅰ級(jí),即D62所在單元為地面沉降高易發(fā)區(qū),其他單元按照同樣的方法計(jì)算確定其易發(fā)級(jí)別。根據(jù)所有單元格的評(píng)價(jià)結(jié)果,再將各剖分單元轉(zhuǎn)化成面元,最后經(jīng)過(guò)等值差分,結(jié)合地質(zhì)環(huán)境條件等綜合劃分出不同級(jí)別的地面沉降地質(zhì)災(zāi)害易發(fā)區(qū)。

    4地面沉降易發(fā)區(qū)劃

    根據(jù)地面沉降地質(zhì)災(zāi)害現(xiàn)狀、發(fā)育程度、危害程度、發(fā)展趨勢(shì)、地質(zhì)環(huán)境條件,結(jié)合地下水開(kāi)采等評(píng)價(jià)因子,按照基于層次分析的可拓學(xué)綜合評(píng)價(jià)等級(jí)大小進(jìn)行易發(fā)區(qū)劃分,地面沉降分為高、中和低3類(lèi)易發(fā)區(qū)。評(píng)價(jià)區(qū)易發(fā)程度分區(qū)結(jié)果如圖1所示。

    圖1 地面沉降易發(fā)程度分區(qū)圖

    4.1高易發(fā)區(qū)(Ⅰ)

    位于德州市中心地帶,面積約200km2。該區(qū)深層地下水位埋深多大于100m,是深層地下水的集中開(kāi)采區(qū)和地面沉降重點(diǎn)區(qū)域。地面沉降累計(jì)沉降量局部大于1000mm,1991—2010年年平均沉降速率為35~62.5mm。該區(qū)深層地下水主要為工業(yè)、第三產(chǎn)業(yè)和城市生活用水。

    該區(qū)應(yīng)嚴(yán)格執(zhí)行深層地下水壓采方案,調(diào)整產(chǎn)業(yè)結(jié)構(gòu);充分利用南水北調(diào)、黃河水源補(bǔ)充水資源;實(shí)施生態(tài)補(bǔ)償原則和用途管制,實(shí)行階梯水價(jià);對(duì)第Ⅲ,Ⅳ含水層組實(shí)施人工回灌;對(duì)地下水、地面沉降進(jìn)行長(zhǎng)期監(jiān)測(cè),加強(qiáng)沉降風(fēng)險(xiǎn)評(píng)估,建立地面沉降監(jiān)測(cè)預(yù)警信息平臺(tái)。

    4.2中易發(fā)區(qū)(Ⅱ)

    位于高易發(fā)區(qū)東南側(cè),面積139km2,是工業(yè)和經(jīng)濟(jì)較為發(fā)達(dá)的地段。該區(qū)深層地下水位埋深90~110m,2010年地面沉降累計(jì)量500~1000mm,1991—2010年年平均沉降速率為20~35mm,是深層地下水開(kāi)采較集中區(qū),是地面沉降和深層地下水降落漏斗的發(fā)展地段。該區(qū)深層地下水主要為工業(yè)、農(nóng)業(yè)、第三產(chǎn)業(yè)和城鄉(xiāng)生活使用。

    該區(qū)宜作為后備水源或應(yīng)急水源,實(shí)行控制性開(kāi)采,盡快建立分層標(biāo)等地面沉降自動(dòng)監(jiān)測(cè)體系。

    4.3低易發(fā)區(qū)(Ⅲ)

    位于研究區(qū)東南邊緣地段以及東部京滬高鐵沿線(xiàn)一帶,面積200km2,是德城區(qū)農(nóng)業(yè)和第三產(chǎn)業(yè)未來(lái)發(fā)展的重要地段。該區(qū)深層地下水位埋深50~90m,2010年地面沉降累計(jì)沉降量300~450mm,1991—2010年年平均沉降速率為10~20mm,深層地下水開(kāi)采比較分散,是地面沉降和深層地下水降落漏斗的周邊地帶。該區(qū)深層地下水主要滿(mǎn)足鄉(xiāng)鎮(zhèn)人民的生活用水,其次是部分小企業(yè)的工業(yè)用水。

    該區(qū)深層地下水實(shí)行控制開(kāi)采,開(kāi)發(fā)利用淺層地下水,適量引用地表水;提倡農(nóng)業(yè)節(jié)水灌溉;可在基巖埋深小的地段建立1組地面沉降監(jiān)測(cè)基巖標(biāo)。

    5結(jié)論

    可拓學(xué)地面沉降易發(fā)性評(píng)價(jià)具有定量化、多因子、綜合性的特點(diǎn)。它引入了描述客觀事物性質(zhì)變化的關(guān)聯(lián)函數(shù)工具,對(duì)實(shí)變函數(shù)中距離的概念加以拓展,為質(zhì)、量兼顧解決2大類(lèi)7小類(lèi)地面沉降易發(fā)性評(píng)價(jià)因子間的矛盾問(wèn)題提供了定量手段。

    該次地面沉降物元模型的建立,實(shí)現(xiàn)了在沉降易發(fā)區(qū)等級(jí)劃分時(shí),將三元組(易發(fā)區(qū)評(píng)價(jià)等級(jí)N,關(guān)于評(píng)價(jià)因子c的量值v)放在一個(gè)統(tǒng)一體中進(jìn)行綜合評(píng)價(jià)。

    綜合考慮地質(zhì)環(huán)境條件、地下水開(kāi)采規(guī)劃所預(yù)測(cè)的累計(jì)沉降量等值線(xiàn)變化趨勢(shì)與可拓學(xué)劃分的地面沉降易發(fā)性分區(qū)形態(tài)具有一定的相似性,反映出地面沉降具有區(qū)域緩變性和受地下水開(kāi)采影響明顯的特點(diǎn)。這一相似性規(guī)律需要進(jìn)一步研究。

    綜上,基于層次分析法的可拓學(xué)理論進(jìn)行地面沉降易發(fā)區(qū)劃,方法是可行的,結(jié)果基本符合當(dāng)?shù)貙?shí)際情況。

    參考文獻(xiàn):

    [1]林良俊,孫曉明,柳富田.新時(shí)期沿海地區(qū)國(guó)土開(kāi)發(fā)面臨的形勢(shì)和主要環(huán)境地質(zhì)問(wèn)題及對(duì)策研究[J].水文地質(zhì)工程地質(zhì),2010,37,(3):1-3.

    [2]孟慶峰,戴魯旗.德州市地面沉降現(xiàn)狀及防治對(duì)策[J].山東國(guó)土資源,2008,24(11):8-9.

    [3]張永偉.華北平原德州地面沉降成生機(jī)理、監(jiān)測(cè)預(yù)警與可控性研究[D].濟(jì)南:山東大學(xué),2014.

    [4]蔡文.物元模型及應(yīng)用[M].北京:科學(xué)技術(shù)文獻(xiàn)出版社,1994.

    [5]XU Y, ZHU QX. A new extension theory-based production operation method in industrial process[J]. Chin J Chem Eng., 2013, 21(1):44-54.

    [6]賈超、肖樹(shù)芳、劉寧,可拓學(xué)理論在洞室?guī)r體質(zhì)量評(píng)價(jià)中的應(yīng)用[J]. 巖石力學(xué)與工程學(xué)報(bào),2003,(5):751-756.

    [7]ZHANG YW, LI SC & MENG FQ. Application of extenics theory for evaluating effect degree of damaged mountains based on analytic hierarchy process[J]. Environ Earth Sci, 2014,71:4463-4471. DOI 10.1007/s12665-013-2840-y.

    [8]賈德旺,岳躍華,葉進(jìn)霞,等.濟(jì)寧城區(qū)地面沉降演變規(guī)律及經(jīng)濟(jì)損失估算方法[J].山東國(guó)土資源,2006,22(11):38-41.

    [9]陳杰,朱國(guó)榮,顧阿明,等.Biot固結(jié)理論在地面沉降計(jì)算中的應(yīng)用[J].水文地質(zhì)工程地質(zhì),2003,(2):28-31.

    [10]龔士良.集對(duì)分析及在城市地而沉降研究中的應(yīng)用[J].水文地質(zhì)工程地質(zhì),1998,(5):33-40.

    [11]李斌,劉峰.遺傳神經(jīng)網(wǎng)絡(luò)在地面沉降趨勢(shì)預(yù)測(cè)中的應(yīng)用.山東國(guó)土資源,2013,29(8):47-49.

    [12]郭小萌,宮輝力,朱峰,等.基于矩陣的GM(1,1)模型預(yù)測(cè)地面沉降——以北京市平原區(qū)為例[J].水文地質(zhì)工程地質(zhì),2013,40(06):101-105.

    [13]楊玉堂,劉清德.聊城市主要地質(zhì)災(zāi)害現(xiàn)狀及防治工作進(jìn)展[J].山東國(guó)土資源,2008,24(1-2):20-21.

    [14]王小剛,鄒祖光,王秀芹,等.東營(yíng)市城區(qū)地面沉降影響因素[J].山東國(guó)土資源,2006,22(5):50-53.

    [15]王哲,易發(fā)成.基于層次分析法的綿陽(yáng)市地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)[J].自然災(zāi)害學(xué)報(bào),2009,18(1):14-23.

    [16]許紹雙,Excel在層次分析法中的應(yīng)用[J]中國(guó)管理信息化,2006,11(9):17-19.

    收稿日期:2015-12-17;

    修訂日期:2016-02-20;編輯:王敏

    基金項(xiàng)目:山東省地面沉降監(jiān)測(cè)與防治研究(工作項(xiàng)目編碼:1212011220182)、華北平原德州地面沉降監(jiān)測(cè)預(yù)警系統(tǒng)建設(shè)(魯國(guó)土資發(fā)[2012]112號(hào))

    作者簡(jiǎn)介:張永偉(1972—),男,山東泗水人,研究員,主要從事水工環(huán)及地質(zhì)災(zāi)害防治等方面的工作;E-mail:zyw526@126.com

    中圖分類(lèi)號(hào):P642.26

    文獻(xiàn)標(biāo)識(shí)碼:B

    Application of Extenincs Theory for Partition of Land Subsidence Susceptible Area of Dezhou City Based on Analytic Hierarchy Process

    ZHANG Yongwei1,2,LI Hongxia3,SHANG Tingting1, JIA Chao2,LIANG Hao1,YANG Yabin3

    (1. Shandong Monitoring Center of Geo-Environment,Jinan,250014,China;2. Shandong University, Jinan,250061, China. 3. Jinan Yuxian primary school,Jinan,250024,China ;4. Lubei Geo-engineering Exploration Institute,Dezhou, 253015,China)

    Abstract:An evaluating method is proposed to the partition of susceptible area of Dezhou city in Huabei plain based on the matter element model as deduced from the theories of extenics. Using the susceptible degree, classification criteria and influencing factors of land subsidence, the elements in classical and sectorized field are defined, whereby an extension-evaluating model is developed for comprehensive evaluation of the classification of susceptible degree of land subsidence by the correlation function of matter elements and extension sets. The analytic hierarchy process methodology for calculating weight aggregation of evaluating factors is proposed. All region are divided into 135 grids using the square(2km×2km), then the split unit is transformed into surface element after equivalent difference. The evaluating results are thus obtained through such an evaluation. The calculation results show that the present method is effective and applicable to forecast the classification of susceptible degree of land subsidence on multiple-factor quantitatively, which lends the support for control planning for land subsidence.

    Key words:Land subsidence; extenics; matter element model; partition of susceptible area; analytic hierarchy process

    猜你喜歡
    評(píng)價(jià)
    SBR改性瀝青的穩(wěn)定性評(píng)價(jià)
    石油瀝青(2021年4期)2021-10-14 08:50:44
    中藥治療室性早搏系統(tǒng)評(píng)價(jià)再評(píng)價(jià)
    自制C肽質(zhì)控品及其性能評(píng)價(jià)
    寫(xiě)作交流與評(píng)價(jià):詞的欣賞
    基于Moodle的學(xué)習(xí)評(píng)價(jià)
    關(guān)于項(xiàng)目后評(píng)價(jià)中“專(zhuān)項(xiàng)”后評(píng)價(jià)的探討
    HBV-DNA提取液I的配制和應(yīng)用評(píng)價(jià)
    有效評(píng)價(jià)讓每朵花兒都綻放
    模糊數(shù)學(xué)評(píng)價(jià)法在水質(zhì)評(píng)價(jià)中的應(yīng)用
    治淮(2013年1期)2013-03-11 20:05:18
    保加利亞轉(zhuǎn)軌20年評(píng)價(jià)
    欧美国产日韩亚洲一区| 欧美又色又爽又黄视频| 国产精品亚洲美女久久久| 色综合欧美亚洲国产小说| 最近最新免费中文字幕在线| 免费看a级黄色片| 精品一区二区三区av网在线观看| 欧美成人性av电影在线观看| 人人妻人人澡欧美一区二区| 国产精品国产高清国产av| 亚洲成a人片在线一区二区| 天美传媒精品一区二区| 精品一区二区三区人妻视频| 国产亚洲精品av在线| 欧美另类亚洲清纯唯美| 亚洲精品影视一区二区三区av| 国产一区二区在线观看日韩| 白带黄色成豆腐渣| 欧美日本亚洲视频在线播放| 日本a在线网址| 免费在线观看影片大全网站| 欧美最新免费一区二区三区 | 18禁裸乳无遮挡免费网站照片| 丰满的人妻完整版| 免费看a级黄色片| 久久香蕉精品热| 免费黄网站久久成人精品 | 国语自产精品视频在线第100页| 天堂网av新在线| 极品教师在线视频| 99在线视频只有这里精品首页| 嫩草影院新地址| 色综合站精品国产| 久久这里只有精品中国| 久久热精品热| 九九在线视频观看精品| 99久久99久久久精品蜜桃| 性欧美人与动物交配| 18美女黄网站色大片免费观看| 国产私拍福利视频在线观看| 国产精品综合久久久久久久免费| 成人高潮视频无遮挡免费网站| 真人做人爱边吃奶动态| 亚洲男人的天堂狠狠| 亚洲欧美日韩高清专用| 欧美成人一区二区免费高清观看| a级毛片免费高清观看在线播放| 久久精品国产自在天天线| 一区二区三区高清视频在线| 黄色日韩在线| 婷婷精品国产亚洲av| 日日夜夜操网爽| 天堂影院成人在线观看| 狂野欧美白嫩少妇大欣赏| 99久久无色码亚洲精品果冻| 成年女人永久免费观看视频| 少妇被粗大猛烈的视频| 夜夜看夜夜爽夜夜摸| 欧美绝顶高潮抽搐喷水| 午夜精品一区二区三区免费看| 成人亚洲精品av一区二区| 午夜老司机福利剧场| av在线天堂中文字幕| 免费av毛片视频| 少妇丰满av| 又紧又爽又黄一区二区| 成人美女网站在线观看视频| 午夜激情福利司机影院| 男人的好看免费观看在线视频| 听说在线观看完整版免费高清| 国产探花在线观看一区二区| 一区二区三区激情视频| 久久久久久久亚洲中文字幕 | 赤兔流量卡办理| 精品人妻一区二区三区麻豆 | 熟女人妻精品中文字幕| 亚洲国产精品成人综合色| 日韩欧美三级三区| 波多野结衣高清无吗| 变态另类成人亚洲欧美熟女| 欧美激情国产日韩精品一区| 亚洲天堂国产精品一区在线| 美女高潮的动态| 18美女黄网站色大片免费观看| 美女xxoo啪啪120秒动态图 | 好男人在线观看高清免费视频| 超碰av人人做人人爽久久| 色吧在线观看| 日韩人妻高清精品专区| 两人在一起打扑克的视频| 性色avwww在线观看| 亚洲av电影在线进入| 国产精品三级大全| 99视频精品全部免费 在线| 波多野结衣巨乳人妻| 中文字幕精品亚洲无线码一区| 亚洲精品影视一区二区三区av| 亚洲熟妇中文字幕五十中出| 一级a爱片免费观看的视频| 亚洲成人免费电影在线观看| 在线国产一区二区在线| 国产成人啪精品午夜网站| 国产精华一区二区三区| 色播亚洲综合网| 偷拍熟女少妇极品色| 国产精品一及| 午夜福利高清视频| 熟妇人妻久久中文字幕3abv| 国产精品电影一区二区三区| 免费看日本二区| 精品一区二区三区视频在线| 18美女黄网站色大片免费观看| 嫩草影院新地址| a级毛片a级免费在线| 欧美午夜高清在线| 美女免费视频网站| 久久欧美精品欧美久久欧美| 精品一区二区三区视频在线观看免费| 日韩精品中文字幕看吧| www.www免费av| 国产高清有码在线观看视频| 亚洲最大成人av| а√天堂www在线а√下载| 超碰av人人做人人爽久久| 精品免费久久久久久久清纯| 一区二区三区四区激情视频 | 美女xxoo啪啪120秒动态图 | 真人做人爱边吃奶动态| 国产欧美日韩精品一区二区| 亚洲午夜理论影院| 国产熟女xx| 欧美精品啪啪一区二区三区| 欧美激情久久久久久爽电影| 亚洲最大成人中文| 少妇熟女aⅴ在线视频| 欧美高清成人免费视频www| 丰满人妻熟妇乱又伦精品不卡| 夜夜夜夜夜久久久久| 美女 人体艺术 gogo| 无人区码免费观看不卡| av视频在线观看入口| 五月玫瑰六月丁香| 国产不卡一卡二| www日本黄色视频网| 少妇被粗大猛烈的视频| 午夜福利在线在线| 日韩欧美国产一区二区入口| 91久久精品电影网| 精品人妻偷拍中文字幕| 国产一级毛片七仙女欲春2| 久久香蕉精品热| 我的老师免费观看完整版| 全区人妻精品视频| 国产高潮美女av| 日本免费一区二区三区高清不卡| 三级国产精品欧美在线观看| 麻豆成人av在线观看| 熟妇人妻久久中文字幕3abv| 赤兔流量卡办理| 天堂网av新在线| 国产av一区在线观看免费| 97超级碰碰碰精品色视频在线观看| 欧美成人一区二区免费高清观看| 搡老妇女老女人老熟妇| 亚洲自偷自拍三级| 国产高清视频在线观看网站| 亚洲综合色惰| 他把我摸到了高潮在线观看| 五月伊人婷婷丁香| 国产精品久久久久久亚洲av鲁大| ponron亚洲| 日本黄色视频三级网站网址| 在线观看免费视频日本深夜| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久久久久久成人| 亚洲欧美日韩东京热| 亚洲国产精品合色在线| 亚洲国产日韩欧美精品在线观看| 99热精品在线国产| 亚洲欧美日韩高清在线视频| xxxwww97欧美| 亚洲av电影在线进入| 99久久无色码亚洲精品果冻| 97超级碰碰碰精品色视频在线观看| 欧美色欧美亚洲另类二区| 看片在线看免费视频| 国产毛片a区久久久久| 亚洲专区国产一区二区| 2021天堂中文幕一二区在线观| 尤物成人国产欧美一区二区三区| 十八禁国产超污无遮挡网站| 日日夜夜操网爽| netflix在线观看网站| 特级一级黄色大片| 在线a可以看的网站| 一区二区三区高清视频在线| 久久国产乱子免费精品| 欧美黑人欧美精品刺激| 国产黄a三级三级三级人| 极品教师在线免费播放| 性色avwww在线观看| 国产中年淑女户外野战色| 蜜桃久久精品国产亚洲av| 久久久久久久亚洲中文字幕 | avwww免费| 宅男免费午夜| 欧美一区二区国产精品久久精品| 亚洲黑人精品在线| 国产真实乱freesex| 一级av片app| 亚洲美女搞黄在线观看 | АⅤ资源中文在线天堂| 我要看日韩黄色一级片| 97热精品久久久久久| 国产成人aa在线观看| 午夜福利成人在线免费观看| 天堂网av新在线| 欧美激情国产日韩精品一区| 成年女人看的毛片在线观看| 亚洲精品在线观看二区| 成人av一区二区三区在线看| 亚洲精品成人久久久久久| 免费观看人在逋| 亚洲男人的天堂狠狠| 国产精品1区2区在线观看.| 国产免费av片在线观看野外av| 午夜福利在线观看吧| 一本精品99久久精品77| 三级国产精品欧美在线观看| 中文字幕熟女人妻在线| 如何舔出高潮| 成人国产综合亚洲| 国产午夜精品久久久久久一区二区三区 | 看十八女毛片水多多多| 亚洲精品成人久久久久久| 一级毛片久久久久久久久女| 一卡2卡三卡四卡精品乱码亚洲| 韩国av一区二区三区四区| 日本免费一区二区三区高清不卡| 精品乱码久久久久久99久播| 97热精品久久久久久| 欧美另类亚洲清纯唯美| 老司机深夜福利视频在线观看| 亚洲精品成人久久久久久| 亚洲狠狠婷婷综合久久图片| 日本撒尿小便嘘嘘汇集6| 国产熟女xx| 婷婷精品国产亚洲av在线| a级毛片a级免费在线| 51国产日韩欧美| 十八禁国产超污无遮挡网站| 日韩中文字幕欧美一区二区| 九九久久精品国产亚洲av麻豆| 亚洲七黄色美女视频| 桃色一区二区三区在线观看| 级片在线观看| 成人国产综合亚洲| 最好的美女福利视频网| 色吧在线观看| 99久久无色码亚洲精品果冻| 麻豆av噜噜一区二区三区| 亚洲第一电影网av| 99久国产av精品| 国产v大片淫在线免费观看| 亚洲人成网站高清观看| 91九色精品人成在线观看| 亚洲在线自拍视频| 少妇的逼水好多| 国产综合懂色| 免费av观看视频| 三级男女做爰猛烈吃奶摸视频| 日韩精品青青久久久久久| 久久久成人免费电影| 亚洲av日韩精品久久久久久密| 毛片一级片免费看久久久久 | 18禁黄网站禁片免费观看直播| www.999成人在线观看| 一个人免费在线观看电影| 国产欧美日韩一区二区精品| 国产v大片淫在线免费观看| 国产成人a区在线观看| 欧美绝顶高潮抽搐喷水| 麻豆久久精品国产亚洲av| 亚洲成人精品中文字幕电影| 宅男免费午夜| 观看美女的网站| 俺也久久电影网| 国产中年淑女户外野战色| 99久国产av精品| 又黄又爽又刺激的免费视频.| 国产亚洲精品综合一区在线观看| av福利片在线观看| 国产三级在线视频| 欧美黄色片欧美黄色片| 岛国在线免费视频观看| 18禁黄网站禁片午夜丰满| 精品人妻一区二区三区麻豆 | 亚洲国产精品999在线| 亚洲熟妇熟女久久| 男人的好看免费观看在线视频| 99热这里只有是精品在线观看 | 久久精品国产亚洲av香蕉五月| 极品教师在线免费播放| av女优亚洲男人天堂| 在线播放国产精品三级| 欧美日韩国产亚洲二区| 免费在线观看影片大全网站| 他把我摸到了高潮在线观看| 男女之事视频高清在线观看| 一级作爱视频免费观看| 欧美国产日韩亚洲一区| 欧美激情久久久久久爽电影| 国产精品影院久久| 哪里可以看免费的av片| 嫩草影院入口| 欧美3d第一页| 又爽又黄a免费视频| 中文资源天堂在线| 我的老师免费观看完整版| 黄色女人牲交| 禁无遮挡网站| 国内精品久久久久精免费| 18+在线观看网站| 色哟哟哟哟哟哟| www.色视频.com| 51午夜福利影视在线观看| 国产精品美女特级片免费视频播放器| 黄色丝袜av网址大全| 国产精品99久久久久久久久| 午夜福利在线在线| 国产国拍精品亚洲av在线观看| www.www免费av| 国产v大片淫在线免费观看| a级一级毛片免费在线观看| 国产极品精品免费视频能看的| 久久这里只有精品中国| 最近视频中文字幕2019在线8| 国产精华一区二区三区| 真实男女啪啪啪动态图| 国产精品乱码一区二三区的特点| 1000部很黄的大片| 宅男免费午夜| 听说在线观看完整版免费高清| 亚洲欧美日韩高清专用| 老司机午夜福利在线观看视频| 日本一本二区三区精品| 亚洲最大成人中文| 成人av一区二区三区在线看| 久久这里只有精品中国| 欧美日韩中文字幕国产精品一区二区三区| 精品欧美国产一区二区三| 精品久久久久久久久久免费视频| 91久久精品电影网| 日日干狠狠操夜夜爽| 精品久久久久久久久久免费视频| 日韩欧美国产一区二区入口| 成人午夜高清在线视频| 老司机午夜福利在线观看视频| 天天一区二区日本电影三级| 村上凉子中文字幕在线| 听说在线观看完整版免费高清| 非洲黑人性xxxx精品又粗又长| 中文字幕久久专区| 久久热精品热| 波多野结衣高清作品| 久久久色成人| 亚洲av成人精品一区久久| 国产精品影院久久| 综合色av麻豆| 精品日产1卡2卡| 日本黄色视频三级网站网址| 中文字幕熟女人妻在线| 91麻豆av在线| 国产午夜福利久久久久久| 嫩草影院精品99| 欧美高清性xxxxhd video| 成人毛片a级毛片在线播放| 欧美性猛交黑人性爽| 久久草成人影院| 亚洲中文字幕日韩| 色综合站精品国产| 免费黄网站久久成人精品 | 综合色av麻豆| 高清在线国产一区| 此物有八面人人有两片| 亚洲五月天丁香| 最新在线观看一区二区三区| 女同久久另类99精品国产91| 久久精品国产清高在天天线| 久久午夜亚洲精品久久| 国产精品电影一区二区三区| 日韩欧美精品v在线| 99国产精品一区二区蜜桃av| 一级黄片播放器| 成人永久免费在线观看视频| 男女下面进入的视频免费午夜| 最后的刺客免费高清国语| 国产高清三级在线| 人妻夜夜爽99麻豆av| 久久久久久大精品| 色av中文字幕| 中文亚洲av片在线观看爽| 久久精品影院6| 精品无人区乱码1区二区| 欧美+日韩+精品| 日本 av在线| 亚洲五月婷婷丁香| 精品国内亚洲2022精品成人| 国产乱人视频| 日本免费a在线| 亚洲人成网站在线播| 一级黄色大片毛片| 亚洲精品色激情综合| 99热这里只有是精品50| 亚洲人成网站在线播放欧美日韩| 97超视频在线观看视频| 男人和女人高潮做爰伦理| 99久久精品一区二区三区| 人人妻人人澡欧美一区二区| 嫩草影院新地址| 久久久久精品国产欧美久久久| 人妻夜夜爽99麻豆av| 综合色av麻豆| 欧美成人免费av一区二区三区| 国产精品98久久久久久宅男小说| 我的女老师完整版在线观看| 成年人黄色毛片网站| 看片在线看免费视频| 深夜a级毛片| 欧美成人一区二区免费高清观看| 日韩大尺度精品在线看网址| av天堂在线播放| 亚洲av熟女| 日本一本二区三区精品| 亚洲精品亚洲一区二区| 日本a在线网址| 亚洲第一电影网av| 成人av一区二区三区在线看| 免费大片18禁| 国产在视频线在精品| 嫩草影视91久久| 岛国在线免费视频观看| 99久久成人亚洲精品观看| 免费在线观看影片大全网站| 18禁在线播放成人免费| 色尼玛亚洲综合影院| 久久久久九九精品影院| 欧美成人一区二区免费高清观看| x7x7x7水蜜桃| 国语自产精品视频在线第100页| 天堂动漫精品| 麻豆av噜噜一区二区三区| 黄色日韩在线| 有码 亚洲区| 免费黄网站久久成人精品 | 最近中文字幕高清免费大全6 | 男女之事视频高清在线观看| 天堂影院成人在线观看| 好男人电影高清在线观看| 久久6这里有精品| 亚洲人成网站高清观看| 国产男靠女视频免费网站| 性插视频无遮挡在线免费观看| 欧美一级a爱片免费观看看| 此物有八面人人有两片| 日本 欧美在线| 真人做人爱边吃奶动态| eeuss影院久久| 亚洲av免费高清在线观看| 搡老熟女国产l中国老女人| 精品一区二区免费观看| 一进一出抽搐动态| 国产精品国产高清国产av| 国产av不卡久久| 特大巨黑吊av在线直播| 精品一区二区三区人妻视频| 性色av乱码一区二区三区2| 天堂av国产一区二区熟女人妻| 国产欧美日韩精品亚洲av| 国产综合懂色| 精品一区二区免费观看| 蜜桃亚洲精品一区二区三区| 成人国产一区最新在线观看| 久久午夜亚洲精品久久| 麻豆成人av在线观看| 狠狠狠狠99中文字幕| 国产av麻豆久久久久久久| 搡女人真爽免费视频火全软件 | 久久久久久久久中文| 日本 av在线| 两人在一起打扑克的视频| 亚洲电影在线观看av| 91久久精品国产一区二区成人| 亚洲午夜理论影院| 欧美黄色片欧美黄色片| 久久草成人影院| 欧美+亚洲+日韩+国产| 欧美成人一区二区免费高清观看| 欧美日本视频| 天堂√8在线中文| 别揉我奶头 嗯啊视频| 午夜福利在线在线| 高清毛片免费观看视频网站| 国内精品一区二区在线观看| 嫩草影院精品99| 国产探花在线观看一区二区| 国产伦精品一区二区三区视频9| 久久99热6这里只有精品| 欧美三级亚洲精品| 欧美在线一区亚洲| 午夜免费男女啪啪视频观看 | 欧美日韩亚洲国产一区二区在线观看| 成人国产一区最新在线观看| 午夜老司机福利剧场| 国产精品乱码一区二三区的特点| 国产三级黄色录像| 国产大屁股一区二区在线视频| 婷婷色综合大香蕉| 国产中年淑女户外野战色| 黄色日韩在线| 亚洲激情在线av| 亚洲成a人片在线一区二区| 精品一区二区三区人妻视频| 日日干狠狠操夜夜爽| 波多野结衣巨乳人妻| 神马国产精品三级电影在线观看| 午夜亚洲福利在线播放| 免费在线观看日本一区| 精品人妻1区二区| 日韩欧美精品免费久久 | 九色成人免费人妻av| 又爽又黄a免费视频| 婷婷六月久久综合丁香| 亚洲一区高清亚洲精品| 在线观看av片永久免费下载| 18禁裸乳无遮挡免费网站照片| 黄色一级大片看看| 舔av片在线| 国产探花在线观看一区二区| 国产伦精品一区二区三区视频9| 国产精品久久久久久久久免 | 九色成人免费人妻av| 国产 一区 欧美 日韩| 欧美bdsm另类| 欧美成人a在线观看| 亚洲av成人精品一区久久| 他把我摸到了高潮在线观看| 如何舔出高潮| 午夜免费成人在线视频| 亚洲欧美日韩高清专用| 亚洲 欧美 日韩 在线 免费| a在线观看视频网站| 黄色日韩在线| 91在线精品国自产拍蜜月| 国产精品不卡视频一区二区 | 日韩欧美在线二视频| eeuss影院久久| 免费看日本二区| 自拍偷自拍亚洲精品老妇| 日韩 亚洲 欧美在线| 国产一级毛片七仙女欲春2| 在线观看午夜福利视频| 亚洲最大成人中文| 中文字幕熟女人妻在线| 99久久精品一区二区三区| 天堂动漫精品| av在线老鸭窝| 国产精品,欧美在线| 久久久久精品国产欧美久久久| 久久婷婷人人爽人人干人人爱| 亚洲精品乱码久久久v下载方式| 国产真实伦视频高清在线观看 | 亚洲成人久久爱视频| 日本免费一区二区三区高清不卡| 熟女人妻精品中文字幕| 国内少妇人妻偷人精品xxx网站| 国产人妻一区二区三区在| 欧美+亚洲+日韩+国产| 99热6这里只有精品| 色在线成人网| 国产久久久一区二区三区| 99久久精品热视频| 性欧美人与动物交配| 亚洲成人中文字幕在线播放| 色综合欧美亚洲国产小说| 久久亚洲精品不卡| 国产av一区在线观看免费| 色吧在线观看| 淫妇啪啪啪对白视频| 精品久久国产蜜桃| 中文字幕av在线有码专区| 国产精品精品国产色婷婷| 9191精品国产免费久久| 我要搜黄色片| 婷婷精品国产亚洲av| 男插女下体视频免费在线播放| 身体一侧抽搐| 好男人在线观看高清免费视频| 国产午夜福利久久久久久| 国产麻豆成人av免费视频| 悠悠久久av| 成年免费大片在线观看| 亚洲片人在线观看| av在线老鸭窝| 一个人免费在线观看的高清视频| 久久精品国产亚洲av涩爱 | 成人无遮挡网站| 变态另类成人亚洲欧美熟女| 高清毛片免费观看视频网站| 老熟妇仑乱视频hdxx| 99久久精品国产亚洲精品| 久久精品国产自在天天线| 久久人妻av系列| 成人特级av手机在线观看| 少妇熟女aⅴ在线视频|