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

    基于InVEST模型的大連市產(chǎn)水量時空變化分析

    2019-10-11 03:50:10呂樂婷任甜甜李賽賽韓月馳
    水土保持通報 2019年4期
    關(guān)鍵詞:瓦房店市普蘭店莊河市

    呂樂婷, 任甜甜, 李賽賽, 韓月馳

    (遼寧師范大學(xué) 城市與環(huán)境學(xué)院, 遼寧 大連 116029)

    水資源是工農(nóng)業(yè)生產(chǎn)、經(jīng)濟(jì)發(fā)展和環(huán)境改善不可替代的自然資源。近年來,由于氣候變化帶來的水源供給的不確定性,嚴(yán)重威脅著生態(tài)系統(tǒng)的安全與穩(wěn)定,影響著自然景觀變化、區(qū)域人口與社會經(jīng)濟(jì)發(fā)展的布局[1]。因此,定量化地評估區(qū)域生態(tài)系統(tǒng)水源供給功能是社會經(jīng)濟(jì)發(fā)展的需要,也是當(dāng)前生態(tài)學(xué)和水文學(xué)交叉領(lǐng)域研究的主要議題[2]。隨著遙感GIS技術(shù)與水文模型的發(fā)展,越來越多的學(xué)者通過水文模型分析和評估區(qū)域生態(tài)系統(tǒng)水源供給功能[3-4],如SWAT模型、InVEST模型等。InVEST模型產(chǎn)水量模塊基于水量平衡原理,可充分考慮不同土地利用類型下土壤滲透性的空間差異性及地形、地表粗糙程度等因素對徑流的影響,以柵格為單元定量估算不同景觀類型的水源供給能力[4]。該模型基于GIS平臺,將量化的生態(tài)系統(tǒng)服務(wù)功能以地圖的形式表現(xiàn)出來,不再用繁瑣的計算公式和過多的文字贅述來表達(dá),且優(yōu)化了很多復(fù)雜的問題[5]。目前已被學(xué)者應(yīng)用于北京山區(qū)[6]、黃土高原[7]、三江源[8]、白龍江流域[9]等生態(tài)系統(tǒng)服務(wù)的評估與研究中。隨著社會經(jīng)濟(jì)發(fā)展,土地利用方式逐漸改變,流域水源供給功能等在內(nèi)的生態(tài)系統(tǒng)服務(wù)正在發(fā)生變化[10]。本文應(yīng)用InVEST模型產(chǎn)水量模塊對不同土地利用下大連市產(chǎn)水量的時空分布特征和影響因子進(jìn)行分析,為大連市水資源的開發(fā)與保護(hù)提供可視化的評估結(jié)果和依據(jù),增強(qiáng)全社會的生態(tài)保護(hù)紅線意識。

    1 研究區(qū)概況

    大連市(120°58′—123°31′E,38°43′—40°10′N)位于遼東半島南端,東臨黃海,西臨渤海,北依東北三省及內(nèi)蒙古廣闊腹地,南與山東半島隔海相望[11]。下屬區(qū)縣包括莊河市、普蘭店區(qū)、瓦房店市、金州區(qū)、長??h、市內(nèi)4區(qū)(中山區(qū)、西崗區(qū)、沙河口區(qū)、甘井子區(qū))及旅順口區(qū),總面積1.30×104km2。該區(qū)位于暖溫帶大陸性季風(fēng)氣候區(qū),冬無嚴(yán)寒,夏無酷暑,四季分明;多年平均氣溫10.5 ℃,年均降水量在550~950 mm之間,降水空間分布由東北向西南逐漸減少。流域內(nèi)耕地、林地資源豐富,耕地資源是大連市面積最大,分布最廣的資源,面積比例為46.43%。以林地、耕地為主的生態(tài)系統(tǒng)對水源的供給功能起著關(guān)鍵性的作用[12]。

    2 研究方法與數(shù)據(jù)處理

    2.1 InVEST模型產(chǎn)水量模型原理

    InVEST模型“產(chǎn)水量”子模塊基于水量平衡原理估算產(chǎn)水量。該模塊將每個柵格單元的產(chǎn)水量定義為該柵格范圍內(nèi)的降水量扣除蒸散量(包括植物蒸騰和地表蒸發(fā))之后,剩余的那一部分水量。然后,假設(shè)每一個柵格單元所有的產(chǎn)水量都通過地表徑流或者地下徑流的方式到達(dá)了流域出口;最后在子流域尺度上對每一個柵格的產(chǎn)水量進(jìn)行加和或求取其平均值。模型主要算法如下:

    (1)

    式中:Yxj——年產(chǎn)水量;Px——年均降雨量; AETxj——土地利用類型j上柵格單元x的年平均實際蒸發(fā)量,由公式(2)計算:

    (2)

    式中:Rxj——Bydyko干燥指數(shù),無量綱,表示潛在蒸發(fā)量與降雨量的比值,由公式(3)計算:

    (3)

    式中:kxj——柵格單元x中土地覆被類型j的植被蒸散系數(shù),由植被葉面積指數(shù)LAI計算獲得; ET0——潛在蒸發(fā)量,由Penman-Monteith公式計算,計算公式為:

    (4)

    式中:ET0——潛在蒸發(fā)量(mm/d);Rn——地表凈輻射〔MJ/(m·d)〕;G——土壤熱通量〔MJ/(m2·d)〕;r——干濕表常數(shù)(kPa/℃);Tmean——日平均溫度(℃);u2——2 m高處風(fēng)速(m/s);es——飽和水氣壓(kPa);ea——實際水氣壓(kPa)。

    (5)

    式中:Z——經(jīng)驗常數(shù),代表區(qū)域降水分布及其他水文地質(zhì)特征。根據(jù)相關(guān)研究,對于冬季降水為主(12月至翌年4月)的地區(qū),Z值接近10,而對于降水均勻分布的濕潤地區(qū)和夏季降水為主的地區(qū),Z值接近1[12]。Z值越高,或干旱指數(shù)較大或較小的地區(qū),模型結(jié)果受季節(jié)常數(shù)Z影響越小[13]; AWCx——柵格單元x的土壤有效含水量,由土壤深度和理化性質(zhì)決定。

    2.2 數(shù)據(jù)來源及處理

    (1) 土地利用數(shù)據(jù)。1980,1990,2000和2015年4期的土地利用利用數(shù)據(jù)(空間分辨率為1 km)均來源于中國科學(xué)院資源環(huán)境數(shù)據(jù)中心(http:∥www.resdc.cn),數(shù)據(jù)生產(chǎn)制作是以各期LandsatTM/ETM遙感影像為主要數(shù)據(jù)源,通過人工目視解譯生成。本研究區(qū)的土地利用類型根據(jù)實際情況歸并調(diào)整為6類,即耕地、林地、草地、水域、城鄉(xiāng)工礦居民用地、未利用地(表1)。

    (2) 氣象數(shù)據(jù)。本研究選取大連市1985—2015年降水、氣溫、風(fēng)速、日照時數(shù)、氣壓數(shù)據(jù),均來源于中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http:∥cdc.cma.gov.cn)。為了提高數(shù)據(jù)的代表性,將數(shù)據(jù)分為4個時段(1980—1985年,1986—1996年,1997—2007年,2008—2015年),1980,1990,2000和2015年的土地利用狀況分別代表4個時段的土地利用狀況。降水要素先進(jìn)行多年平均,然后通過反距離加權(quán)插值得到多年平均降水量柵格數(shù)據(jù);其它氣候要素均選用日值,采用FAO修正的Penman-Monteith公式計算[14]潛在蒸發(fā)量,通過反距離加權(quán)插值法得到多年平均潛在蒸發(fā)量柵格數(shù)據(jù)。

    (3) 土壤厚度和土壤有效含水量。土壤厚度數(shù)據(jù)來源于聯(lián)合國糧農(nóng)組織和維也納國際應(yīng)用系統(tǒng)研究所構(gòu)建的世界土壤數(shù)據(jù)庫,中國境內(nèi)數(shù)據(jù)源為第二次全國土地調(diào)查南京土壤所提供的1∶100萬土壤數(shù)據(jù);土壤有效含水量在分析土壤質(zhì)地(砂土—壤土—黏土)百分比含量基礎(chǔ)上,利用土壤有效含水量經(jīng)驗公式在SPAW軟件中計算獲得。

    (4) 生物物理參數(shù)表。生物物理參數(shù)表反映了研究區(qū)土地利用和土地覆蓋類型的屬性,包括土地利用編碼、植被最大深度、蒸散系數(shù)等。這些數(shù)據(jù)根據(jù)前人的研究成果[15]和InVEST模型數(shù)據(jù)庫資料獲得。

    表1 1980-2015年大連市土地利用類型%

    編號土地利用類型 1980年比例1990年比例2000年比例2015年比例1耕地47.0047.0346.4946.432林地39.7639.3438.8937.953草地1.521.471.451.384水域4.053.984.164.185城鄉(xiāng)工礦居民用地7.317.998.969.976未利用地0.350.190.060.1

    3 結(jié)果與分析

    3.1 產(chǎn)水量的時間變化特征

    大連市1980—2015年降水、氣溫、潛在蒸發(fā)、實際蒸散、干燥系數(shù)(實際蒸散量占降水量的比例)、產(chǎn)水量的年際變化趨勢如圖1所示。研究區(qū)多年平均降水量677.50 mm,最高降水量1 007.01 mm(1985年),最低降水量403.39 mm(2015年),年均降水量呈不顯著下降趨勢,多年平均氣溫9.64 ℃,呈顯著上升趨勢(p=0.025),上升速率為0.179 ℃/10 a;年均潛在蒸散量947.74 mm,呈不顯著下降趨勢,降低速率為12.67 mm/10 a;年均實際蒸散量445.78.mm,呈顯著下降趨勢(p=0.001),降低速率為19.22 mm/10 a;年均干燥系數(shù)0.72,呈顯著上升趨勢,上升速率為0.057/10 a,這意味著大部分降水通過實際蒸散返回到了大氣當(dāng)中。在氣象要素變化的影響下,研究區(qū)多年平均產(chǎn)水量226.36 mm,占降水量33.41%,呈不顯著下降趨勢,降低速率為14.7 mm/10 a。

    圖1 大連市降水、氣溫、潛在蒸散、實際蒸散、干燥系數(shù)及產(chǎn)水量時間變化

    經(jīng)過相關(guān)分析得出(表2),1980—2015年,研究區(qū)產(chǎn)水量與降水量呈顯著正相關(guān),相關(guān)系數(shù)為0.793(p<0.01);與氣溫、潛在蒸散、干燥系數(shù)呈顯著負(fù)相關(guān),相關(guān)系數(shù)分別為-0.517(p<0.01),-0.624(p<0.01),-0.486(p<0.01)。潛在蒸發(fā)反映了水分充分供給條件下最大的蒸散量,受氣溫的顯著影響(R=0.344,p<0.05)。實際蒸散代表了綜合考慮氣候要素及下墊面狀況下的蒸發(fā)能力,受降水和潛在蒸發(fā)影響,與降水呈顯著正相關(guān)(R=0.365,p<0.05),與潛在蒸發(fā)呈不顯著負(fù)相關(guān),相關(guān)系數(shù)為-0.086,也與地表覆蓋類型有關(guān)。綜合分析,本研究認(rèn)為降水量是影響大連市產(chǎn)水量的主要氣候因素,同時受氣溫、潛在蒸發(fā)與干燥系數(shù)的影響。

    表2 1980-2015年大連市產(chǎn)水量與氣候要素的相關(guān)性

    注:**表示在0.01水平,相關(guān)性顯著。

    3.2 產(chǎn)水量的空間變化特征

    大連市氣象要素、土地利用類型、產(chǎn)水量空間分布特征如圖2所示。由圖2可知,1980—2015年,大連市年均降水量自東北向西南597.57~803.29 mm。其中,莊河市年均降水量最豐富,達(dá)733.65 mm,其次為普蘭店區(qū),達(dá)672.09 mm,瓦房店市和南部7區(qū)縣的降水量次之,分別為644.73,643.94.mm。潛在蒸發(fā)空間特征與降水量相反,自東北向西南遞增858.33~1 036.39 mm。其中,南部7區(qū)縣潛在蒸發(fā)量最大,為998.94 mm;瓦房店市、普蘭店區(qū)的潛在蒸發(fā)量次之,分別為954.06,938.30 mm;莊河市的潛在蒸發(fā)量最小,為881.92 mm。實際蒸散空間分布與土地利用類型、降水基本一致,不同土地利用類型單位面積實際蒸散量不同,具體表現(xiàn)為:水域(616.74 mm)>林地(508.87 mm)>耕地(423.55 mm)>草地(422.86 mm)>居民地(306.83 mm)>未利用地(409.57 mm)。莊河市降水量居大連市之首,用地類型以水域、耕地、林地為主,蒸散發(fā)能力強(qiáng),因此年均實際蒸散量最高,達(dá)457.22 mm,分別是南部七區(qū)縣、普蘭店區(qū)、瓦房店市的1.12,1.00和1.02倍(圖2,表3)。

    受降水、蒸散發(fā)、土地利用(圖2)的綜合影響,大連市產(chǎn)水量總體變化趨勢為由東北向西南波動變化。高值區(qū)主要集中在莊河市東北部、瓦房店市南部,多年平均產(chǎn)水量在295.46~611.34 mm之間。莊河市東北部降水豐富,產(chǎn)水量大;瓦房店南部、南部7區(qū)縣中部以城市建設(shè)用地為主,蒸發(fā)能力弱,產(chǎn)水量亦最大;莊河市、瓦房店市北部及南部7區(qū)縣東北部、普蘭店區(qū)產(chǎn)水量次之,集中在168.62~295.46 mm之間;南部7區(qū)縣西南部、瓦房店市東北部以耕地和林地為主,蒸散系數(shù)高,對降水的截留能力強(qiáng),產(chǎn)水量相對較小。從產(chǎn)水總量來看,莊河市因其降水量最大,土地面積廣闊,年均產(chǎn)水量最大,達(dá)276.78 mm,占降水量的37%,產(chǎn)水能力最高,分別是南部7區(qū)縣、普蘭店區(qū)、瓦房店市的1.16,1.29,1.39倍(表3)。

    表3 1980-2015年大連市年均氣象要素及產(chǎn)水量

    4 氣候變化和土地利用變化對產(chǎn)水量的影響

    4.1 土地利用結(jié)構(gòu)及轉(zhuǎn)移方向分析

    由表4可知,大連市耕地和林地面積所占比例較大,分別為46%,38%,其次為城市建設(shè)用地,達(dá)10%,而草地、水體及未利用地所占的比例較小,分別為1.4%,4.2%和0.1%。1980—2015年,耕地、林地、未利用地面積均有不同程度的減小。其中,耕地面積減小87 km2,減幅1%,主要59%轉(zhuǎn)向建設(shè)用地及29%轉(zhuǎn)向林地;林地面積減小216 km2,減幅2%,主要95%轉(zhuǎn)向耕地。城鎮(zhèn)建設(shè)用地、草地、水域的面積均有不同程度的增加。草地、水域面積增幅較小,分別為0.1%和0.2%。城鎮(zhèn)建設(shè)用地面積增加339 km2,增幅2.7%,主要由60%以上的耕地轉(zhuǎn)移而來。從總體特征來看,大連市域是一個典型的快速城市化區(qū)域,建設(shè)用地36 a間快速增長,這主要是因為在快速城鎮(zhèn)化過程中,耕地轉(zhuǎn)為建設(shè)用地,而由于自然地理位置的原因,其它土地利用類型轉(zhuǎn)換為建設(shè)用地的面積很小[16]。

    表4 1980-2015年大連市土地利用類型轉(zhuǎn)移矩陣 km2

    圖2 大連市氣象要素、土地利用類型(2015年)、產(chǎn)水量空間分布特征

    4.2 情景設(shè)置

    為了區(qū)分1980—2015年研究區(qū)氣候變化和土地利用變化對產(chǎn)水量的影響,本研究選擇兩個階段,第一階段是1980—1990年,另一階段是2005—2015年。1980年和2015年的土地利用圖用來表示這兩個時期的土地利用模式。本研究設(shè)置了4種情景。情景1:1980—1990年的氣候和1980年的土地利用;情景2:1980—1990年的氣候和2015年的土地利用;情景3:2005—2015年的氣候和1980年的土地利用;情景4:2005—2015年的氣候和2015年的土地利用。情景1和情景4分別表示兩個時期的實際情況;情景1—2揭示土地利用變化對產(chǎn)水量的影響;情景1—3描述氣候變化對產(chǎn)水量的影響;情景4—1描述土地利用和氣候變化對產(chǎn)水量的綜合影響。

    4.3 不同情景下產(chǎn)水量的時間變化

    表5顯示了1980—2015年不同情景下研究區(qū)的產(chǎn)水量以及變化量。情景2和情景1的產(chǎn)水量差異表明,土地利用變化的影響增加了6.97 mm,占情景1條件下產(chǎn)水量的5%。情景3和情景1的產(chǎn)水量差異代表了氣候變化的效應(yīng),減少了16.92 mm,占情景1條件下產(chǎn)水量的12%。與情景1相比,情景4下土地利用和氣候變化的綜合效應(yīng)下產(chǎn)水量減少了10.09 mm,占情景1產(chǎn)水量的7%。結(jié)果表明,土地利用和氣候變化都影響了研究區(qū)的產(chǎn)水量,氣候變化的影響大于土地利用變化的影響。當(dāng)氣候變化效應(yīng)發(fā)揮主導(dǎo)作用時,土地利用變化對產(chǎn)水量的影響并不明顯。

    4.4 不同情景下產(chǎn)水量的空間變化

    圖3顯示了不同情境下大連市產(chǎn)水量的空間分布圖,表6顯示了不同情景下不同地區(qū)的產(chǎn)水量??傮w上,在不同情景下,產(chǎn)水量的空間分布趨勢基本一致。在情景1條件下,平均產(chǎn)水量142.75 mm。其中,莊河市產(chǎn)水量最大,為194.92 mm;長??h、普蘭店區(qū)的產(chǎn)水量次之,分別為180.78,131.85 mm;旅順口區(qū)的產(chǎn)水量最小,為75.54 mm。在情景2條件下,土地利用為主要驅(qū)動因素,平均產(chǎn)水量為149.72 mm,相比情景1,增加了6.97 mm,絕大部分地區(qū)的產(chǎn)水量呈增加趨勢。

    表5 不同情景下大連市的產(chǎn)水量及變化值

    圖3 不同情景下產(chǎn)水量的空間分布

    其中,市內(nèi)4區(qū)耕地面積減少了46 km2,林地面積減少24 km2,城市建設(shè)用地增加了73 km2,相應(yīng)平均產(chǎn)水量增加最大,為30 mm;其次為金州區(qū),耕地、林地面積分別減少了55和15 km2,城市建設(shè)用地增加了78 km2,相應(yīng)平均產(chǎn)水量增加了16.84 mm;旅順口區(qū)、瓦房店市、普蘭店區(qū),耕地、林地、城市建設(shè)用地面積變化較小,產(chǎn)水量亦變化較小,平均產(chǎn)水量分別增加了6.78,6.5,4.08 mm;莊河市耕地、林地面積分別減少8和48 km2,城市建設(shè)用地增加64 km2,平均產(chǎn)水量增加最小,分別為3.17 mm。綜合分析,林地和耕地面積的下降,城市建設(shè)用地的增加,對大連市產(chǎn)水量有些微的影響。在情景3條件下,氣候變化為主要驅(qū)動因素,平均產(chǎn)水量為125.83 mm,相比情景1,減少了16.92 mm,絕大部分地區(qū)的產(chǎn)水量呈減少趨勢,尤其莊河市的產(chǎn)水量減少最大,為42.82 mm,主要原因在于1980—2015年期間,氣溫增加明顯,降水量減少最大,達(dá)65 mm以上。在情景4條件下,氣候變化和土地利用為主要驅(qū)動因素,平均產(chǎn)水量為132.66 mm,相比情景1,減少了10.09 mm。莊河市、長??h、普蘭店產(chǎn)水量呈減少趨勢,其中,莊河市減少量最大,為39.69 mm。金州區(qū)、旅順口區(qū)、市內(nèi)4區(qū)、瓦房店產(chǎn)水量呈增加趨勢,其中,市內(nèi)4區(qū)產(chǎn)水量增加最大,為40.55 mm。結(jié)果表明,不同情景下,產(chǎn)水量的空間分布趨勢基本一致。研究區(qū)土地利用的變化對產(chǎn)水量的增加有些微影響,氣候變化是影響產(chǎn)水量變化的主要驅(qū)動力。

    表6 不同情景下不同地區(qū)的產(chǎn)水量mm

    地 區(qū)情景1情景2情景3情景4金州區(qū) 116.78133.62118.82135.34旅順口區(qū)75.5482.3282.1488.93普蘭店區(qū)131.85135.93118.02121.94瓦房店市115.06121.56109.63115.95長海縣 180.78180.78156.30156.30莊河市 194.92198.09152.12155.23市內(nèi)4區(qū)130.89161.89140.45171.44

    5 結(jié) 論

    (1) 1980—2015年,大連市年均降水量677.50 mm,年均潛在蒸發(fā)947.74 mm,二者不顯著下降,下降速率分別為22.22和12.67 mm/10 a。平均氣溫9.64 ℃,年均干燥系數(shù)為72%,二者顯著增加,增加速率分別為0.179和0.057/10 a;年均實際蒸散量445.78.mm,顯著下降,降低速率19.22 mm/10 a。研究區(qū)多年平均產(chǎn)水量226.36 mm,占降水量33.41%。在研究時段內(nèi),平均產(chǎn)水量不顯著下降,降低速率14.7 mm/10 a。在降水波動的影響下,不同年代,各地區(qū)潛在蒸散、實際蒸散量波動變化;產(chǎn)水量變化趨勢與降水變化趨勢一致,都呈先降低后增加的趨勢。降水量是影響大連市產(chǎn)水量的主要氣候因素,同時受潛在蒸發(fā)量與干燥系數(shù)的影響。

    (2) 1980—2015年大連市產(chǎn)水量總體上由東北向西南波動變化。高值區(qū)主要集中在莊河市東北部、瓦房店市南部、南部7區(qū)縣中部,多年平均產(chǎn)水量在295.46~611.34 mm之間;莊河市、瓦房店市北部及南部7區(qū)縣東北部、普蘭店區(qū)產(chǎn)水量次之,集中在168.62~295.46 mm之間;南部7區(qū)縣西南部、瓦房店市東北部產(chǎn)水量相對較小,集中在1.11~168.62 mm之間。從產(chǎn)水總量來看,莊河市年均產(chǎn)水總量最大,達(dá)1.10×109m3。

    (3) 1980—2015年,研究期間大連市耕地、林地面積有所減少。城市建設(shè)用地面積有所增加,主要由耕地轉(zhuǎn)移而來。其它用地類型面積變化較小。

    (4) 氣候變化和土地利用變化是引起大連市產(chǎn)水量時空變化的兩個主要因素,氣候變化的影響大于土地利用變化的影響,當(dāng)氣候變化效應(yīng)發(fā)揮主導(dǎo)作用時,土地利用變化對產(chǎn)水量的影響并不明顯。氣候因子主要通過改變降水和潛在蒸散來影響產(chǎn)水量,降水的增加可以提高產(chǎn)水量的源頭水量,潛在蒸散表征區(qū)域生態(tài)系統(tǒng)的水分消耗能力。土地利用的變化可以影響實際蒸散,進(jìn)而影響產(chǎn)水量。

    (5) InVEST模型的引入,為大中尺度流域的水源供給服務(wù)的估算與空間分布提供了定量和可視化的參考,但由于模型結(jié)構(gòu)的簡化以及缺乏長期的野外觀測數(shù)據(jù),致使研究結(jié)果存在一定的不準(zhǔn)確性。因此,建議在今后的研究工作中,加強(qiáng)野外數(shù)據(jù)的觀測、調(diào)整模型參數(shù)使其本土化,確保評估結(jié)果的可靠性。

    猜你喜歡
    瓦房店市普蘭店莊河市
    普蘭店鼓樂生存現(xiàn)狀與發(fā)展研究
    輕音樂(2022年3期)2022-12-29 06:19:48
    普蘭店區(qū)唐房革命烈士陵園
    瓦房店市農(nóng)村飲水安全存在問題與改進(jìn)措施
    基于小城鎮(zhèn)生活圈的城鄉(xiāng)公共服務(wù)設(shè)施配置優(yōu)化
    The Study on the Ability Training of Learner Autonomy in Middle School English Study
    大連普蘭店破獲假冒商標(biāo)案
    莊河市農(nóng)村留守兒童福利問題研究
    遼寧省莊河市河長制管理制度問題探析
    完善農(nóng)機(jī)監(jiān)理工作?護(hù)航現(xiàn)代農(nóng)業(yè)發(fā)展
    遼寧省瓦房店市泡崖小學(xué)第三屆“鴻潤杯”2016春聯(lián)大賽優(yōu)秀作品選
    小布老虎(2016年5期)2016-02-28 18:14:01
    午夜福利在线免费观看网站| 最新的欧美精品一区二区| 亚洲九九香蕉| 精品久久久久久电影网| 99精品在免费线老司机午夜| 国产亚洲av高清不卡| 午夜影院日韩av| 在线免费观看的www视频| 国产精品自产拍在线观看55亚洲 | 亚洲精品在线美女| 国产成人一区二区三区免费视频网站| 91精品三级在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲精品自拍成人| 久久精品国产亚洲av香蕉五月 | 国产一区在线观看成人免费| 在线观看一区二区三区激情| 巨乳人妻的诱惑在线观看| 两个人免费观看高清视频| 精品视频人人做人人爽| 一区二区日韩欧美中文字幕| 日韩三级视频一区二区三区| 如日韩欧美国产精品一区二区三区| 中文欧美无线码| 高清在线国产一区| 亚洲专区国产一区二区| 午夜免费观看网址| 国产精品偷伦视频观看了| 欧美中文综合在线视频| 女性被躁到高潮视频| 免费黄频网站在线观看国产| 色94色欧美一区二区| 一区在线观看完整版| 黄片播放在线免费| 久久精品亚洲精品国产色婷小说| av线在线观看网站| 久久中文看片网| 色精品久久人妻99蜜桃| 淫妇啪啪啪对白视频| 亚洲五月色婷婷综合| 午夜福利欧美成人| 欧美精品人与动牲交sv欧美| 法律面前人人平等表现在哪些方面| 亚洲午夜理论影院| 国产亚洲精品第一综合不卡| 午夜亚洲福利在线播放| 精品午夜福利视频在线观看一区| 欧美性长视频在线观看| 国产成人系列免费观看| 国产成人精品无人区| 成人国产一区最新在线观看| 久久久国产一区二区| 久久性视频一级片| 国产成人影院久久av| 亚洲五月天丁香| 国产区一区二久久| 中文字幕色久视频| 两人在一起打扑克的视频| 亚洲免费av在线视频| 国产亚洲欧美精品永久| 国产日韩一区二区三区精品不卡| 国产精品电影一区二区三区 | 亚洲一区二区三区欧美精品| 大码成人一级视频| 亚洲一区中文字幕在线| 欧美黑人精品巨大| 欧美亚洲 丝袜 人妻 在线| 国内久久婷婷六月综合欲色啪| 女性被躁到高潮视频| 国产又爽黄色视频| 99精国产麻豆久久婷婷| 成年女人毛片免费观看观看9 | 19禁男女啪啪无遮挡网站| 如日韩欧美国产精品一区二区三区| 黑人欧美特级aaaaaa片| 丰满饥渴人妻一区二区三| 国产免费av片在线观看野外av| 男女下面插进去视频免费观看| 女人被狂操c到高潮| 欧美 亚洲 国产 日韩一| 日韩成人在线观看一区二区三区| 男女免费视频国产| 黄色丝袜av网址大全| 亚洲中文字幕日韩| 男女床上黄色一级片免费看| 国产欧美日韩综合在线一区二区| 欧美乱码精品一区二区三区| 日韩成人在线观看一区二区三区| 18在线观看网站| 在线播放国产精品三级| 黄色a级毛片大全视频| 成人国语在线视频| 成人手机av| 精品一区二区三区四区五区乱码| 人妻久久中文字幕网| 美女福利国产在线| 老司机在亚洲福利影院| 午夜福利欧美成人| 国产成人影院久久av| 久久ye,这里只有精品| 国产精品99久久99久久久不卡| 国产真人三级小视频在线观看| 少妇裸体淫交视频免费看高清 | 欧美人与性动交α欧美精品济南到| 最新美女视频免费是黄的| 中文亚洲av片在线观看爽 | 黄色毛片三级朝国网站| 色在线成人网| 亚洲五月天丁香| 久久久久久久久免费视频了| 亚洲专区中文字幕在线| 免费高清在线观看日韩| 视频在线观看一区二区三区| 国产免费现黄频在线看| 麻豆成人av在线观看| 黄网站色视频无遮挡免费观看| 少妇 在线观看| 久久久久国产精品人妻aⅴ院 | 日韩欧美三级三区| 精品国产乱码久久久久久男人| 国产极品粉嫩免费观看在线| 欧美一级毛片孕妇| 窝窝影院91人妻| 国产精品99久久99久久久不卡| 欧美av亚洲av综合av国产av| 欧美在线黄色| www.999成人在线观看| 免费观看人在逋| 一级作爱视频免费观看| 日韩一卡2卡3卡4卡2021年| 热99re8久久精品国产| 免费在线观看亚洲国产| 午夜久久久在线观看| 色综合婷婷激情| 中文亚洲av片在线观看爽 | 日本黄色日本黄色录像| 飞空精品影院首页| 亚洲专区字幕在线| 黄色视频,在线免费观看| 国产精品久久久久久人妻精品电影| av视频免费观看在线观看| 亚洲精品粉嫩美女一区| 国产亚洲精品久久久久5区| 国产麻豆69| 国精品久久久久久国模美| 久久久精品区二区三区| 在线免费观看的www视频| 老司机影院毛片| 制服诱惑二区| 丰满人妻熟妇乱又伦精品不卡| 午夜91福利影院| 少妇粗大呻吟视频| 欧美老熟妇乱子伦牲交| 99精国产麻豆久久婷婷| 国产主播在线观看一区二区| av免费在线观看网站| 99国产精品一区二区三区| 国产亚洲精品久久久久久毛片 | 波多野结衣一区麻豆| bbb黄色大片| 国产麻豆69| 欧美日韩亚洲国产一区二区在线观看 | 高清在线国产一区| 窝窝影院91人妻| 国产aⅴ精品一区二区三区波| 悠悠久久av| 亚洲国产欧美网| 免费观看精品视频网站| 熟女少妇亚洲综合色aaa.| 丝袜美足系列| 精品少妇一区二区三区视频日本电影| 丝瓜视频免费看黄片| 国产在线观看jvid| 成人亚洲精品一区在线观看| 欧美日韩成人在线一区二区| 亚洲国产中文字幕在线视频| a级毛片黄视频| 久久天堂一区二区三区四区| 两人在一起打扑克的视频| 亚洲精品一卡2卡三卡4卡5卡| 99re6热这里在线精品视频| aaaaa片日本免费| 下体分泌物呈黄色| 久久久久久久午夜电影 | 少妇猛男粗大的猛烈进出视频| videos熟女内射| 午夜福利在线观看吧| 久久久久久久久免费视频了| 免费少妇av软件| 中国美女看黄片| 国产精品电影一区二区三区 | 久久影院123| 女人被躁到高潮嗷嗷叫费观| 中国美女看黄片| 亚洲色图综合在线观看| 日韩 欧美 亚洲 中文字幕| 色精品久久人妻99蜜桃| 欧美乱码精品一区二区三区| 99精品欧美一区二区三区四区| 久久这里只有精品19| 亚洲午夜精品一区,二区,三区| 欧美精品亚洲一区二区| 精品第一国产精品| 日本五十路高清| 亚洲av日韩在线播放| 久久热在线av| 久久国产乱子伦精品免费另类| 亚洲人成电影观看| 在线十欧美十亚洲十日本专区| 亚洲成av片中文字幕在线观看| 两人在一起打扑克的视频| av免费在线观看网站| 亚洲熟妇熟女久久| 欧美中文综合在线视频| 新久久久久国产一级毛片| 精品福利永久在线观看| 脱女人内裤的视频| av欧美777| 一区在线观看完整版| 女性被躁到高潮视频| 精品乱码久久久久久99久播| 校园春色视频在线观看| 黄色 视频免费看| 免费女性裸体啪啪无遮挡网站| 正在播放国产对白刺激| 国产不卡一卡二| 精品熟女少妇八av免费久了| 黑人巨大精品欧美一区二区mp4| 国产成人精品在线电影| 久久青草综合色| 国产不卡一卡二| 黄色片一级片一级黄色片| 国产精品免费视频内射| 亚洲一码二码三码区别大吗| 王馨瑶露胸无遮挡在线观看| 亚洲精品美女久久久久99蜜臀| 久久久久国产一级毛片高清牌| 超碰97精品在线观看| 日本a在线网址| 久久青草综合色| 久久精品aⅴ一区二区三区四区| 日本wwww免费看| 国产精品久久久久成人av| 国产1区2区3区精品| 女警被强在线播放| 怎么达到女性高潮| 亚洲精品中文字幕在线视频| 另类亚洲欧美激情| 大陆偷拍与自拍| 欧美日韩中文字幕国产精品一区二区三区 | 国产xxxxx性猛交| 51午夜福利影视在线观看| 精品午夜福利视频在线观看一区| 曰老女人黄片| 日韩欧美三级三区| 中文字幕人妻熟女乱码| 日韩大码丰满熟妇| 中文字幕最新亚洲高清| 国产精品98久久久久久宅男小说| 久热这里只有精品99| 久久精品亚洲av国产电影网| 久久国产精品影院| 国产伦人伦偷精品视频| 亚洲国产欧美网| 深夜精品福利| www.熟女人妻精品国产| 91成人精品电影| 如日韩欧美国产精品一区二区三区| 亚洲午夜精品一区,二区,三区| 国产男女超爽视频在线观看| 日日夜夜操网爽| 午夜免费观看网址| 欧美日韩亚洲国产一区二区在线观看 | 成年人免费黄色播放视频| 多毛熟女@视频| 波多野结衣av一区二区av| 操出白浆在线播放| 欧美亚洲日本最大视频资源| 中文欧美无线码| 777久久人妻少妇嫩草av网站| 成在线人永久免费视频| 色综合婷婷激情| 国产精品二区激情视频| xxx96com| 一本一本久久a久久精品综合妖精| 国产精品免费大片| 美女午夜性视频免费| 色婷婷av一区二区三区视频| 嫩草影视91久久| 桃红色精品国产亚洲av| 少妇 在线观看| 亚洲av熟女| 精品高清国产在线一区| www.999成人在线观看| 国产精品永久免费网站| 国产av一区二区精品久久| 91麻豆精品激情在线观看国产 | 欧美日韩视频精品一区| 成人特级黄色片久久久久久久| 大型av网站在线播放| 亚洲av片天天在线观看| 亚洲av日韩精品久久久久久密| 久久国产精品人妻蜜桃| 国产蜜桃级精品一区二区三区 | 日本精品一区二区三区蜜桃| 国产精品 欧美亚洲| 91成人精品电影| 午夜老司机福利片| 每晚都被弄得嗷嗷叫到高潮| 久久久国产成人免费| 日本黄色视频三级网站网址 | 欧美 日韩 精品 国产| 国产黄色免费在线视频| av免费在线观看网站| 亚洲国产精品一区二区三区在线| 91在线观看av| 久久中文字幕一级| 一本大道久久a久久精品| 国产午夜精品久久久久久| 91老司机精品| 亚洲久久久国产精品| 91大片在线观看| 国产精品香港三级国产av潘金莲| 婷婷丁香在线五月| 波多野结衣av一区二区av| aaaaa片日本免费| 亚洲 国产 在线| 国产欧美日韩一区二区精品| 一进一出抽搐gif免费好疼 | 王馨瑶露胸无遮挡在线观看| 成在线人永久免费视频| 亚洲国产欧美网| 国产一区二区激情短视频| 久久久久久久久免费视频了| 亚洲久久久国产精品| 久久国产精品人妻蜜桃| 日本精品一区二区三区蜜桃| √禁漫天堂资源中文www| 久久人人97超碰香蕉20202| 午夜影院日韩av| 中文字幕制服av| 少妇裸体淫交视频免费看高清 | 波多野结衣av一区二区av| 亚洲五月婷婷丁香| 老熟妇乱子伦视频在线观看| 亚洲av日韩精品久久久久久密| 国产精品 欧美亚洲| 老汉色∧v一级毛片| 老司机影院毛片| 国产激情欧美一区二区| 国产欧美亚洲国产| 亚洲欧洲精品一区二区精品久久久| 日韩欧美国产一区二区入口| 国产麻豆69| 好男人电影高清在线观看| 在线观看免费午夜福利视频| 韩国精品一区二区三区| x7x7x7水蜜桃| 十八禁高潮呻吟视频| 一进一出抽搐动态| 女人高潮潮喷娇喘18禁视频| 久久精品国产99精品国产亚洲性色 | 久久香蕉激情| 成人国产一区最新在线观看| 精品无人区乱码1区二区| 老熟女久久久| 夜夜夜夜夜久久久久| 午夜影院日韩av| 欧美黄色淫秽网站| 在线视频色国产色| 人妻久久中文字幕网| 亚洲国产看品久久| 欧美中文综合在线视频| 欧美午夜高清在线| 高潮久久久久久久久久久不卡| 99热只有精品国产| 亚洲精品久久午夜乱码| 岛国毛片在线播放| 热re99久久国产66热| 老司机影院毛片| 精品久久久精品久久久| 国产伦人伦偷精品视频| 老司机在亚洲福利影院| 高清毛片免费观看视频网站 | 丝瓜视频免费看黄片| 亚洲人成伊人成综合网2020| videos熟女内射| 亚洲成国产人片在线观看| 精品久久久久久,| 大码成人一级视频| 在线永久观看黄色视频| 久久久久久久久免费视频了| 亚洲中文av在线| 亚洲国产欧美网| 丝袜美腿诱惑在线| 亚洲午夜精品一区,二区,三区| 电影成人av| 国产精品98久久久久久宅男小说| 国产精华一区二区三区| 欧美 日韩 精品 国产| bbb黄色大片| av中文乱码字幕在线| 亚洲精品乱久久久久久| 久久久精品区二区三区| 欧美老熟妇乱子伦牲交| 精品久久久久久久久久免费视频 | 国产三级黄色录像| 中文字幕色久视频| 无遮挡黄片免费观看| 女人精品久久久久毛片| av线在线观看网站| 美女 人体艺术 gogo| 少妇粗大呻吟视频| 亚洲一卡2卡3卡4卡5卡精品中文| 999精品在线视频| 777米奇影视久久| 国产高清视频在线播放一区| 国产精品久久久久久人妻精品电影| 中文字幕av电影在线播放| 日韩成人在线观看一区二区三区| 十八禁人妻一区二区| 一级毛片女人18水好多| 精品视频人人做人人爽| 久久久国产成人免费| 女同久久另类99精品国产91| 亚洲七黄色美女视频| 丁香六月欧美| 一进一出好大好爽视频| 自拍欧美九色日韩亚洲蝌蚪91| 99国产极品粉嫩在线观看| 久久午夜亚洲精品久久| 亚洲欧美精品综合一区二区三区| 99精国产麻豆久久婷婷| 窝窝影院91人妻| 欧美黑人欧美精品刺激| 国内久久婷婷六月综合欲色啪| 亚洲美女黄片视频| 国产精品美女特级片免费视频播放器 | 不卡一级毛片| av福利片在线| 国产野战对白在线观看| 黑人猛操日本美女一级片| 国产精品免费一区二区三区在线 | 国产欧美日韩一区二区三| 人人澡人人妻人| 久久精品人人爽人人爽视色| 亚洲国产精品合色在线| 成人手机av| 欧美+亚洲+日韩+国产| 亚洲精品乱久久久久久| 午夜免费观看网址| 一区二区三区国产精品乱码| 国产精华一区二区三区| 亚洲综合色网址| 成年动漫av网址| 久久国产精品人妻蜜桃| 人人妻人人澡人人爽人人夜夜| 中文字幕精品免费在线观看视频| 免费看十八禁软件| 亚洲专区中文字幕在线| 欧美成人免费av一区二区三区 | 99久久精品国产亚洲精品| 美女午夜性视频免费| 欧美大码av| 国产精品.久久久| 国产在视频线精品| 国产黄色免费在线视频| 欧美+亚洲+日韩+国产| 亚洲第一欧美日韩一区二区三区| 亚洲成人国产一区在线观看| 亚洲成人手机| 亚洲aⅴ乱码一区二区在线播放 | 国产免费现黄频在线看| 又黄又爽又免费观看的视频| 香蕉国产在线看| 日韩制服丝袜自拍偷拍| 亚洲国产中文字幕在线视频| 国产一区二区三区综合在线观看| 亚洲自偷自拍图片 自拍| 这个男人来自地球电影免费观看| 伊人久久大香线蕉亚洲五| 首页视频小说图片口味搜索| 18禁国产床啪视频网站| 美女福利国产在线| 啪啪无遮挡十八禁网站| 在线观看www视频免费| x7x7x7水蜜桃| 午夜福利一区二区在线看| 国产精品电影一区二区三区 | 黄网站色视频无遮挡免费观看| 热99re8久久精品国产| 国产又色又爽无遮挡免费看| 看免费av毛片| 中文字幕人妻丝袜制服| 亚洲精品成人av观看孕妇| 视频区图区小说| 搡老岳熟女国产| 黑人欧美特级aaaaaa片| 国内毛片毛片毛片毛片毛片| 女同久久另类99精品国产91| 亚洲人成伊人成综合网2020| 国产蜜桃级精品一区二区三区 | 高清黄色对白视频在线免费看| 伦理电影免费视频| 国产一区有黄有色的免费视频| 欧美精品av麻豆av| 久久久久久久午夜电影 | 日韩精品免费视频一区二区三区| 人人澡人人妻人| 丝袜美腿诱惑在线| 天堂俺去俺来也www色官网| 视频区图区小说| 在线观看免费视频网站a站| 高潮久久久久久久久久久不卡| 国产精品av久久久久免费| 女同久久另类99精品国产91| 免费在线观看亚洲国产| 久久精品国产清高在天天线| 亚洲欧美一区二区三区久久| 亚洲一区二区三区欧美精品| 亚洲精品在线观看二区| 中文字幕最新亚洲高清| 极品教师在线免费播放| 一区二区三区激情视频| 久久久国产成人精品二区 | 男人操女人黄网站| 一级a爱片免费观看的视频| 国产日韩一区二区三区精品不卡| 久久中文字幕一级| 动漫黄色视频在线观看| 最新在线观看一区二区三区| 一本一本久久a久久精品综合妖精| 日韩免费高清中文字幕av| 黄片大片在线免费观看| 新久久久久国产一级毛片| 制服人妻中文乱码| 亚洲人成77777在线视频| 老司机在亚洲福利影院| 国产成人一区二区三区免费视频网站| 国产蜜桃级精品一区二区三区 | 久久精品91无色码中文字幕| 久久精品国产a三级三级三级| 久久久久精品国产欧美久久久| 久久久久久久午夜电影 | 91麻豆精品激情在线观看国产 | 在线观看一区二区三区激情| 免费观看精品视频网站| 亚洲黑人精品在线| 老熟女久久久| 国产aⅴ精品一区二区三区波| 丝袜美腿诱惑在线| 不卡一级毛片| 亚洲avbb在线观看| 亚洲av日韩精品久久久久久密| 成年版毛片免费区| x7x7x7水蜜桃| 亚洲,欧美精品.| 老司机深夜福利视频在线观看| 男人操女人黄网站| 国产淫语在线视频| 丝袜美足系列| 在线观看免费日韩欧美大片| 免费黄频网站在线观看国产| 免费观看精品视频网站| 老鸭窝网址在线观看| 亚洲色图av天堂| 亚洲av第一区精品v没综合| 成人av一区二区三区在线看| 久久久国产一区二区| 国产精品.久久久| 亚洲欧美日韩另类电影网站| 亚洲一区高清亚洲精品| 欧美日本中文国产一区发布| 9热在线视频观看99| 18禁黄网站禁片午夜丰满| 两个人免费观看高清视频| 9191精品国产免费久久| a在线观看视频网站| 国产亚洲精品第一综合不卡| 精品少妇久久久久久888优播| 亚洲av熟女| 日韩欧美国产一区二区入口| xxx96com| 欧美在线一区亚洲| 欧美+亚洲+日韩+国产| 成人18禁高潮啪啪吃奶动态图| 18禁美女被吸乳视频| 欧美日韩亚洲高清精品| 午夜成年电影在线免费观看| 99精国产麻豆久久婷婷| 黄色 视频免费看| 成人18禁在线播放| 中文字幕另类日韩欧美亚洲嫩草| 国产一区二区三区视频了| www.精华液| 久久中文字幕人妻熟女| 女人被躁到高潮嗷嗷叫费观| 伦理电影免费视频| 亚洲精品中文字幕一二三四区| 成人国产一区最新在线观看| xxx96com| 啪啪无遮挡十八禁网站| 午夜影院日韩av| 日日爽夜夜爽网站| 1024视频免费在线观看| 亚洲欧美色中文字幕在线| 99国产精品99久久久久| 国产精品成人在线| 国产高清激情床上av| 日韩免费高清中文字幕av| 操出白浆在线播放| 免费在线观看完整版高清| 成年版毛片免费区| 亚洲一区高清亚洲精品| 国产精品.久久久|