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

    基于BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型的城市安全距離研究
    ——以蕪湖和馬鞍山為例

    2016-10-13 19:44:26寧,董欣,寧雄,劉
    中國(guó)環(huán)境科學(xué) 2016年6期
    關(guān)鍵詞:馬鞍山建成區(qū)間距

    賈 寧,董 欣,寧 雄,劉 毅

    ?

    基于BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型的城市安全距離研究
    ——以蕪湖和馬鞍山為例

    賈 寧,董 欣,寧 雄,劉 毅*

    (清華大學(xué)環(huán)境學(xué)院,北京 100084)

    快速的城鎮(zhèn)化使城市間距不斷減小,上下游城市間水質(zhì)的相互影響愈發(fā)顯著.本研究提出基于水質(zhì)的城市安全距離概念,建立基于BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型的城市安全距離量化方法,并選擇長(zhǎng)江沿岸相鄰的蕪湖、馬鞍山兩市為案例,評(píng)估未來(lái)兩市建成區(qū)擴(kuò)張后城市間距的安全性,計(jì)算兩市建成區(qū)的最小安全距離.結(jié)果顯示,2020年蕪湖與馬鞍山4.6km的間距屬于安全距離,能夠保證下游城市馬鞍山上游控制斷面地表II類(lèi)的水功能要求.但與2010年相比,控制斷面水質(zhì)變差,COD與氨氮濃度分別提高了29.2%與23.2%.為了保證控制斷面的水功能要求,蕪湖與馬鞍山兩市的最小安全距離為3.2km.

    城市安全距離;城市空間增長(zhǎng)邊界;BP神經(jīng)網(wǎng)絡(luò);水質(zhì)模型

    中國(guó)正經(jīng)歷著快速的城鎮(zhèn)化[1],城市規(guī)模逐漸擴(kuò)大,建設(shè)用地不斷擴(kuò)張,城市連片分布,隨之而來(lái)的,是日益凸顯的資源環(huán)境壓力[2-3].研究表明,土地由非建設(shè)用地轉(zhuǎn)變?yōu)榻ㄔO(shè)用地,將加大水資源需求[4],加劇暴雨徑流[5],增大污水排放[6],成為入湖河流水質(zhì)污染的主要來(lái)源[7].與此同時(shí),城市的擴(kuò)張使得相鄰城市建成區(qū)的間距不斷縮小,上游城市的污染排放對(duì)下游城市飲用水安全的威脅不能通過(guò)河流自?xún)舳?上下游城市間的水安全影響越發(fā)顯著.由此可見(jiàn),在設(shè)定城市空間增長(zhǎng)邊界(UGB)并將其作為一種控制城市無(wú)序擴(kuò)張、引導(dǎo)城市合理發(fā)展的手段[8]時(shí),不僅需要考慮城市擴(kuò)張對(duì)本地水資源與水環(huán)境的壓力,還需要考慮對(duì)下游城市水安全的影響,特別是在城市密集的區(qū)域.

    現(xiàn)有針對(duì)城市空間增長(zhǎng)邊界的研究大多是對(duì)單個(gè)城市空間邊界的量化[9-11],模擬過(guò)程中以本地的生態(tài)安全格局[12-14]和資源環(huán)境承載能力[15]等為約束,很少涉及城市間的相互影響.賈濱洋等[16]提出在確定城市空間增長(zhǎng)邊界時(shí),應(yīng)當(dāng)考慮城市間的最小生態(tài)安全距離.

    在上述城市空間增長(zhǎng)邊界的研究中,對(duì)城市間相互影響的考慮基本集中在概念探討與定性分析層面,缺少定量化的研究.

    基于此,本研究在城市空間增長(zhǎng)邊界設(shè)定中引入安全距離的概念,綜合考慮城市擴(kuò)張對(duì)其自身污染排放的影響及上下游城市間的水質(zhì)響應(yīng)關(guān)系,以下游城市水環(huán)境功能定位為約束目標(biāo),利用基于BP神經(jīng)網(wǎng)絡(luò)的水質(zhì)模型[17-19],確定上下游城市間的合理距離,進(jìn)而確定城市空間增長(zhǎng)的合理邊界.

    1 基于BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型的城市安全距離量化方法

    安全距離的概念常見(jiàn)于石油化工行業(yè)建設(shè)項(xiàng)目的環(huán)境影響評(píng)價(jià)[20-21].在本研究中,安全距離是基于水環(huán)境質(zhì)量的安全距離,指上游城市的正常排放不會(huì)影響下游城市水環(huán)境功能的城市間距.對(duì)于具體城市而言,本研究所定義的安全距離不具有唯一性,但具有最小值:當(dāng)上下游城市間距小于最小安全距離時(shí),上游城市正常的污染排放就會(huì)威脅到下游城市的水環(huán)境質(zhì)量.本研究構(gòu)建的BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型的城市安全距離量化方法如圖1所示,包括3部分:城市污染排放特征概化、BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型構(gòu)建和城市間距安全性判斷.

    1.1 城市污染排放特征概化

    城市污染排放特征概化是基于城市土地利用變化模擬結(jié)果對(duì)研究區(qū)域點(diǎn)源與面源排放強(qiáng)度及空間分布進(jìn)行的量化.點(diǎn)源包括生活源與工業(yè)源兩類(lèi),其排放量通過(guò)研究區(qū)域的居民生活強(qiáng)度與工業(yè)行業(yè)特征確定.非點(diǎn)源指雨水沖刷建成區(qū)和城市之間非建成區(qū)產(chǎn)生的地表徑流帶入水體的污染物,其排放量通過(guò)場(chǎng)次降雨平均沖刷模型進(jìn)行計(jì)算.在本研究中,假設(shè)點(diǎn)源均通過(guò)污水管網(wǎng)收集至污水處理廠,處理后就近排入城市水環(huán)境;非點(diǎn)源則以線源的形式沿程排入水體.

    1.2 BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型構(gòu)建

    河流某斷面的污染物濃度與斷面所在位置以及上游污染物排放位置和排放量P之間存在著某種函數(shù)關(guān)系,如式(1)所示.BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型就是通過(guò)監(jiān)測(cè)數(shù)據(jù)的樣本訓(xùn)練確定函數(shù)關(guān)系,進(jìn)而實(shí)現(xiàn)斷面污染物濃度的預(yù)測(cè).

    =(,P,) (1)

    在本研究中,BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型的輸入即其輸入層包括:河流上游來(lái)水的污染物濃度與流量,點(diǎn)源排放量及排放位置,非點(diǎn)源線型排放強(qiáng)度及排放起始點(diǎn)位置,以及模擬斷面的位置;輸出即其輸出層為模擬斷面的污染物濃度;隱含層神經(jīng)元的個(gè)數(shù)依據(jù)輸入輸出神經(jīng)元的個(gè)數(shù)具體確定.本研究采用串行訓(xùn)練的方式對(duì)BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型的所有訓(xùn)練樣本同時(shí)進(jìn)行訓(xùn)練,如果總體誤差值滿(mǎn)足收斂要求,則訓(xùn)練終止;如果不滿(mǎn)足要求,則進(jìn)行反向的權(quán)值修正過(guò)程,具體過(guò)程見(jiàn)圖2.

    1.3 城市間距安全性判斷

    利用 BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型,對(duì)下游城市控制斷面的污染物濃度進(jìn)行預(yù)測(cè).如果控制斷面污染物濃度滿(mǎn)足下游城市水環(huán)境功能的要求,則上下游城市的間距屬于安全距離,也就是說(shuō),城市間河段的環(huán)境容量與自?xún)裟芰梢跃徑庥缮舷掠纬鞘锌臻g擴(kuò)張、間距縮小而引發(fā)的水環(huán)境壓力.如果控制斷面污染物濃度不能滿(mǎn)足下游城市水環(huán)境功能的要求,則需要重新選擇上下游城市空間增長(zhǎng)模式,重新設(shè)定城市空間增長(zhǎng)邊界,直至下游城市控制斷面水質(zhì)滿(mǎn)足要求,此時(shí)得到的上下游城市間距為安全距離.

    2 案例應(yīng)用與結(jié)果

    2.1 案例區(qū)域基本情況

    選取蕪湖和馬鞍山兩市作為研究對(duì)象.兩市位于安徽省東南部,是安徽省重要的經(jīng)濟(jì)中心,也是長(zhǎng)江下游重要的港口城市.2012年,兩市人口分別為356.6萬(wàn)和218.8萬(wàn),城鎮(zhèn)化率分別達(dá)到58.0%和61.2%.2000年以來(lái),兩市建設(shè)用地增長(zhǎng)迅速.蕪湖市建設(shè)用地總面積由2000年的68km2增長(zhǎng)到2010年的135km2,馬鞍山市的增長(zhǎng)更為迅速,從36km2增長(zhǎng)至92km2.

    快速的城鎮(zhèn)化使得兩市建成區(qū)呈現(xiàn)沿江相向擴(kuò)張的趨勢(shì)、間距不斷縮小.2010年國(guó)家發(fā)改委批復(fù)的《皖江城市帶承接產(chǎn)業(yè)轉(zhuǎn)移示范區(qū)規(guī)劃》[22]中提出,“加快蕪(湖)馬(鞍山)同城化進(jìn)程”,“建設(shè)現(xiàn)代化濱江組團(tuán)式大城市”.可以預(yù)見(jiàn),未來(lái)兩地建成區(qū)將進(jìn)一步相向擴(kuò)張,城市間距將進(jìn)一步縮小.而兩地主要的飲用水源和污水去向均為長(zhǎng)江干流,因此這種城市間距的不斷減小會(huì)使上游蕪湖的污染排放給下游馬鞍山的水安全帶來(lái)潛在威脅.

    本研究綜合考慮蕪湖市與馬鞍山市建成區(qū)的分布,提取長(zhǎng)江東岸1063km2的城市用地密集區(qū)作為研究區(qū)域.該區(qū)域2020年的用地情況如圖3所示,建設(shè)用地預(yù)計(jì)達(dá)到490km2.

    2.2 案例區(qū)域污染排放特征概化

    研究區(qū)域內(nèi)包括8個(gè)水質(zhì)監(jiān)測(cè)斷面與5個(gè)污水處理廠服務(wù)區(qū).根據(jù)蕪湖市與馬鞍山市建成區(qū)的范圍、污水處理廠排放口等信息,可將研究區(qū)域污染排放的空間特征概化,如圖4所示.研究選取下游城市馬鞍山市的上游斷面為控制斷面,以此斷面水質(zhì)表征城市間距變化對(duì)水環(huán)境質(zhì)量的影響.本研究以最上游監(jiān)測(cè)斷面為起點(diǎn),江岸任意點(diǎn)至起始點(diǎn)的沿江距離為該點(diǎn)的坐標(biāo).由于長(zhǎng)江江面較寬,因此,本研究假設(shè)岸邊的水質(zhì)僅受同側(cè)岸邊污染排放的影響,忽略江對(duì)岸排放的影響.

    研究區(qū)域內(nèi)的點(diǎn)源為5座污水處理廠的尾水.污水廠的污水來(lái)源于生活排放和工業(yè)排放.

    生活排放的計(jì)算基于人口,如式(2):

    1365(2)

    式中:1為生活污染物排放總量;為區(qū)域人口數(shù)量;為人均日用水量,L/d;為產(chǎn)污系數(shù),即污水排放量與用水量的比值;為排放的污水中污染物的濃度,mg/L.

    工業(yè)排放基于工業(yè)用地面積計(jì)算,如式(3):

    2(3)

    式中:2為工業(yè)污染物排放總量;為單位工業(yè)用地面積的污染物排放量,t/km2,考慮技術(shù)進(jìn)步因素,排放強(qiáng)度呈逐年下降趨勢(shì);為工業(yè)用地面積,基于實(shí)際2005、2010年以及預(yù)測(cè)得到的2020年數(shù)據(jù),采用線性插值法得到逐年工業(yè)的用地面積.

    假設(shè)生活污水和工業(yè)廢水均經(jīng)處理后達(dá)一級(jí)B標(biāo)準(zhǔn)排放[23].自2006年以來(lái),研究區(qū)域內(nèi)的5座污水處理廠均是在原有基礎(chǔ)上不斷擴(kuò)建,并未改變位置,因此,假設(shè)到2020年,5座污水處理廠仍在原位置擴(kuò)建以滿(mǎn)足城市擴(kuò)張的需求,即點(diǎn)源的排放量發(fā)生變化,但位置不變.基于此,在后續(xù)構(gòu)建的BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型中,對(duì)于點(diǎn)源,僅將5個(gè)排放量作為輸入變量,數(shù)據(jù)見(jiàn)表1.

    區(qū)域內(nèi)的非點(diǎn)源排放以行政區(qū)為單元進(jìn)行量化.本文采取場(chǎng)次降雨徑流平均濃度計(jì)算非點(diǎn)源污染物的排放量,見(jiàn)式(4):

    ×EMC (4)

    式中:指污染物排放總量;指徑流系數(shù),指建設(shè)用地面積;指降雨量,數(shù)據(jù)來(lái)自地方部門(mén);EMC為場(chǎng)次降雨平均濃度.

    假設(shè)2個(gè)區(qū)域的非點(diǎn)源分別以均勻線型源的形式沿城市邊界排入長(zhǎng)江,則非點(diǎn)源沿江排放強(qiáng)度的計(jì)算方法見(jiàn)式(5):

    /(5)

    式中:為沿江排放強(qiáng)度;為建成區(qū)非點(diǎn)源污染物排放總量;為建成區(qū)沿江跨度.

    此外, 還要考慮非點(diǎn)源污染物排放的位置.研究區(qū)域的建成區(qū)可分為蕪湖和馬鞍山兩段.運(yùn)用GIS平臺(tái)將建成區(qū)轉(zhuǎn)換成矩形,則建成區(qū)的起點(diǎn)、終點(diǎn)可用坐標(biāo)表示.基于實(shí)際2005、2010年以及預(yù)測(cè)2020年的坐標(biāo)值,采用線性插值的方法獲得其他年份的坐標(biāo)信息.

    因此,在后續(xù)構(gòu)建的BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型中,對(duì)于非點(diǎn)源,需要將排放強(qiáng)度及建成區(qū)的起、終點(diǎn)信息作為輸入變量,數(shù)據(jù)見(jiàn)表2.

    2.3 “蕪湖—馬鞍山”段BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型構(gòu)建

    表1 研究區(qū)域內(nèi)不同年份點(diǎn)源排放強(qiáng)度(t/a)Table 1 Discharge intensities of regional point sources in different years (t/a)

    根據(jù)上述污染特征概化結(jié)果,“蕪湖—馬鞍山”段水質(zhì)模型的輸入變量共有14個(gè):作為點(diǎn)源的5座污水處理廠的污染排放量(1~5),非點(diǎn)源污染線型排放強(qiáng)度(1,2)及排放區(qū)域起點(diǎn)終點(diǎn)坐標(biāo)(1~4),上游來(lái)水流量(0)及污染物濃度(0),斷面位置.模型訓(xùn)練和驗(yàn)證時(shí),斷面位置輸入圖4所示8個(gè)監(jiān)測(cè)斷面之一的坐標(biāo)(1,2, …,或8).用模型預(yù)測(cè)時(shí),斷面位置輸入控制斷面的坐標(biāo)(’).水質(zhì)模型的輸出變量為斷面污染物濃度.模型訓(xùn)練和驗(yàn)證時(shí),輸出為某個(gè)監(jiān)測(cè)斷面的污染物濃度(1,1,…,或8);預(yù)測(cè)時(shí),輸出為控制斷面濃度(’).模型對(duì)COD和氨氮兩種污染物進(jìn)行模擬.所構(gòu)建的BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型如圖5所示.

    整理2006~2013年間長(zhǎng)江干流蕪湖馬鞍山段8個(gè)監(jiān)測(cè)斷面的COD和氨氮年均濃度數(shù)據(jù),剔除部分未檢出的情況,共得到37組“年份-斷面位置-濃度”數(shù)據(jù).每輪隨機(jī)選出30組樣本用于模型訓(xùn)練,剩余7組樣本用于模型驗(yàn)證.重復(fù)進(jìn)行500輪訓(xùn)練,每輪產(chǎn)生的30組訓(xùn)練樣本不完全相同.設(shè)置訓(xùn)練的收斂條件為:總誤差小于0.01,或迭代次數(shù)達(dá)2000次,滿(mǎn)足其中之一即可.

    模型訓(xùn)練過(guò)程中模擬誤差的收斂軌跡如圖6(a)所示.500輪訓(xùn)練的結(jié)果表明,樣本的選擇對(duì)模型模擬誤差的影響不大,COD的模擬誤差集中在0.01~0.03之間,氨氮的模擬誤差集中在0.01~ 0.04之間,見(jiàn)圖6(b).可見(jiàn)該模型具有較高的穩(wěn)定性和可靠性.

    (a)模型訓(xùn)練誤差收斂過(guò)程

    (b) 500輪訓(xùn)練的收斂誤差

    圖6 “蕪湖—馬鞍山”段BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型訓(xùn)練誤差
    Fig.6 Training errors of the BP neural network water quality model in Wuhu-Ma'anshan River Reach

    從圖7可知,COD與氨氮濃度的模擬誤差均可以控制在10%以?xún)?nèi),模型性能較好.

    2.4 “蕪湖—馬鞍山”間距安全性判斷

    以蕪湖市建成區(qū)末端與馬鞍山市建成區(qū)起始點(diǎn)之間的距離表征2個(gè)城市的間距.2010年,兩城市間距為9.85km;2020年,兩城市間距則縮短為4.64km.利用上述BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型可以對(duì)控制斷面的水質(zhì)進(jìn)行預(yù)測(cè),結(jié)果表明,2020年控制斷面的COD年均值為12.80mg/L,氨氮濃度年均值為0.31mg/L.馬鞍山市對(duì)控制斷面的水功能定位為地表II類(lèi),從COD和氨氮2個(gè)指標(biāo)的年均值來(lái)看,2020年控制斷面水質(zhì)能夠滿(mǎn)足水功能要求,也就是說(shuō)未來(lái)蕪湖市與馬鞍山市4.6km的間距屬于安全距離,現(xiàn)階段兩市城鎮(zhèn)化過(guò)程中建成區(qū)相向擴(kuò)張的模式較為合理.然而與2010年相比,2020年控制斷面COD提高了29.2%,氨氮的濃度提高了23.2%.由此可見(jiàn),雖然2020年控制斷面的水質(zhì)能夠滿(mǎn)足當(dāng)?shù)氐乃δ芤?但與現(xiàn)狀相比,惡化較快,如果按照現(xiàn)有的城市擴(kuò)張模式發(fā)展,未來(lái)控制斷面的水質(zhì)將不能滿(mǎn)足水功能要求.

    上述模型結(jié)果表明,控制斷面的COD增長(zhǎng)較快,因此,本研究以COD為控制目標(biāo),進(jìn)一步對(duì)蕪湖與馬鞍山兩市間的最小安全距離進(jìn)行了分析.假設(shè)2020年2個(gè)城市的社會(huì)經(jīng)濟(jì)發(fā)展情景不變,可以認(rèn)為2個(gè)城市的點(diǎn)源與非點(diǎn)源排放總量不變,僅調(diào)整上游蕪湖市建成區(qū)的終點(diǎn)位置,即調(diào)整兩個(gè)城市建成區(qū)的間距,模擬分析不同城市間距下控制斷面的COD,結(jié)果如圖8所示.

    從圖8可以看出,控制斷面的COD濃度隨著兩市建成區(qū)間距的減小而升高.當(dāng)兩市建成區(qū)的間距小于2.5km時(shí),下游馬鞍山市的控制斷面將不能滿(mǎn)足地表II類(lèi)水的功能要求.考慮到突發(fā)性污染事故、城鎮(zhèn)化速度過(guò)快等因素對(duì)控制斷面水質(zhì)的影響,選取水功能要求濃度限值的90%對(duì)應(yīng)的建成區(qū)間距為2個(gè)城市的最小安全距離,可以得到蕪湖與馬鞍山兩市建成區(qū)的最小安全距離為3.2km.也就是說(shuō),為了保證控制斷面的水功能要求,兩城市在未來(lái)城市邊界增長(zhǎng)過(guò)程中應(yīng)保持大于3.2km的間距.已有研究對(duì)蕪湖與馬鞍山兩市2020年城市擴(kuò)張情景進(jìn)行了預(yù)測(cè),預(yù)計(jì)到2020年兩市建成區(qū)間距為4.6km,可見(jiàn)按照現(xiàn)有的城市建設(shè)模式發(fā)展,到2020年馬鞍山控制斷面的水質(zhì)能夠滿(mǎn)足當(dāng)?shù)氐乃δ苄枨?仍有空間可供兩個(gè)城市相向擴(kuò)張,見(jiàn)圖9.

    在城市邊界增長(zhǎng)的過(guò)程中考慮水環(huán)境質(zhì)量與上下游城市間距的響應(yīng)關(guān)系,量化基于水環(huán)境質(zhì)量的城市安全距離,既可以判斷城市間距離是否合理、是否預(yù)留足夠的緩沖地帶,也能夠?yàn)樽R(shí)別城市發(fā)展空間格局的環(huán)境影響、實(shí)現(xiàn)可持續(xù)的城鎮(zhèn)化提供技術(shù)支持與決策依據(jù).

    3 結(jié)語(yǔ)

    本研究在提出基于水環(huán)境質(zhì)量的城市安全距離的基礎(chǔ)上,建立了包括污染排放特征概化、BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型構(gòu)建及城市間距安全性判斷三步的基于BP神經(jīng)網(wǎng)絡(luò)水質(zhì)模型的城市安全距離量化方法,并以長(zhǎng)江流域“蕪湖—馬鞍山”區(qū)域?yàn)榘咐?進(jìn)行了兩市建成區(qū)間距安全性的分析,量化了兩市建成區(qū)的最小安全距離.結(jié)果顯示,2020年蕪湖與馬鞍山4.6km的間距屬于安全距離,能夠保證下游城市馬鞍山上游控制斷面地表II類(lèi)的水功能要求.但與2010年相比,控制斷面水質(zhì)變差,COD與氨氮濃度分別提高了29.2%與23.2%.為了保證控制斷面的水功能要求,蕪湖與馬鞍山兩市的最小安全距離為3.2km.研究為蕪湖、馬鞍山2個(gè)城市的空間邊界增長(zhǎng)與空間格局開(kāi)發(fā)提出了約束.

    Gregory B. Aspects of Urbanization in China: Shanghai, Hong Kong, Guangzhou [M]. Amsterdam: Amsterdam University Press, 2012:13-24.

    Tang Z, Engel B A, Pijanowski B C, et al. Forecasting land use change and its environmental impact at a watershed scale [J]. Journal of environmental management, 2005,76(1):35-45.

    Duh J D, Shandas V, Chang H, et al. Rates of urbanisation and the resiliency of air and water quality [J]. Science of the Total Environment, 2008,400(1):238-256.

    Liu G, Yang Z, Chen B, et al. Emergy-based dynamic mechanisms of urban development, resource consumption and environmental impacts [J]. Ecological Modelling, 2014,271:90- 102.

    Fletcher T, Burns M. Urban stormwater runoff: a new class of environmental flow problem [J]. PLoS ONE [P]. 2012,7(9):1-10.

    Dong Y, Liu Y, Chen J. Will urban expansion lead to an increase in future water pollution loads?—a preliminary investigation of the Haihe River Basin in northeastern China [J]. Environmental Science and Pollution Research, 2014,21(11):7024-7034.

    孫金華,曹曉峰,黃 藝.滇池流域土地利用對(duì)入湖河流水質(zhì)的影響 [J]. 中國(guó)環(huán)境科學(xué), 2011,31(12):2052-2057.

    Benfield F K, Terris J, Vorsanger N. Solving sprawl: Models of smart growth in communities across America [M]. Washington, DC: Island Press, 2003:137-138.

    徐 康,吳紹華,陳東湘,等.基于水文效應(yīng)的城市增長(zhǎng)邊界的確定——以鎮(zhèn)江新民洲為例 [J]. 地理科學(xué), 2013,8:979-985.

    Tayyebi A, Pijanowski B C, Tayyebi A H. An urban growth boundary model using neural networks, GIS and radial parameterization: An application to Tehran, Iran [J]. Landscape and Urban Planning, 2011,100(1):35-44.

    Bhatta B. Modelling of urban growth boundary using geoinformatics [J]. International Journal of Digital Earth, 2009, 2(4):359-381.

    李詠華.生態(tài)視角下的城市增長(zhǎng)邊界劃定方法——以杭州市為例 [J]. 城市規(guī)劃, 2011,12:83-90.

    周 銳,王新軍,蘇海龍,等.基于生態(tài)安全格局的城市增長(zhǎng)邊界劃定——以平頂山新區(qū)為例 [J]. 城市規(guī)劃學(xué)刊, 2014,4:57- 63.

    張曉瑞,賀巖丹,方創(chuàng)琳,等.城市生態(tài)環(huán)境脆弱性的測(cè)度分區(qū)與調(diào)控 [J]. 中國(guó)環(huán)境科學(xué), 2015,35(7):2200-2208.

    王振波,張 薔,張曉瑞,等.基于資源環(huán)境承載力的合肥市增長(zhǎng)邊界劃定 [J]. 地理研究, 2013,12:2302-2311.

    賈濱洋,曾九利,李 玫,等.“多規(guī)融合”下的城市開(kāi)發(fā)邊界與最小生態(tài)安全距離 [J]. 環(huán)境保護(hù), 2015,Z1:23-26.

    宋國(guó)浩.人工神經(jīng)網(wǎng)絡(luò)在水質(zhì)模擬與水質(zhì)評(píng)價(jià)中的應(yīng)用研究 [D]. 重慶:重慶大學(xué), 2009.

    Zhao Y, Nan J, Cui F, et al. Water quality forecast through application of BP neural network at Yuqiao reservoir [J]. Journal of Zhejiang University SCIENCE A, 2007,8(9):1482-1487.

    王 祎,李靜文,邵 雪,等.基于計(jì)算智能的流域污染排放優(yōu)化模式研究 [J]. 中國(guó)環(huán)境科學(xué), 2012,32(1):173-180.

    李 娜,張 震,冷 飛.油田生產(chǎn)設(shè)施環(huán)境安全距離理論與案例分析 [J]. 油氣田環(huán)境保護(hù), 2012,22(6):76-80.

    劉艷菊,張 震,李 巍,等.油田生產(chǎn)設(shè)施環(huán)境安全距離理論與估算——以勝利油田東辛采油廠為例 [J]. 安全與環(huán)境工作, 2006,13(1):77-84.

    國(guó)家發(fā)展改革委.關(guān)于印發(fā)皖江城市帶承接產(chǎn)業(yè)轉(zhuǎn)移示范區(qū)規(guī)劃的通知[EB/OL]. http://www.gov.cn/zwgk/2010-03/24/ content_1563507.htm, 2010-3-24/2015-12-1.

    GB18918-2002 城鎮(zhèn)污水處理廠污染物排放標(biāo)準(zhǔn) [S].

    *責(zé)任作者, 教授, yi.liu@tsinghua.edu.cn

    Safe distance between cities based on the BP neural network water quality model: a study on Wuhu and Ma'anshan

    JIA Ning, DONG Xin, NING Xiong, LIU Yi*

    (School of Environment, Tsinghua University, Beijing 10084, China)., 2016,36(6):1905~1912

    As a result of decreasing distance between cities during the rapid urbanization, the water quality of cities located in upper and lower reaches has closer interations. This paper proposes the concept of safe distance between cities to ensure the water quality in downstream cities, which is quantified based on the BP neural network model for water quality. Two adjacent cities along the Yangtze River, Wuhu and Ma'anshan, are chosen as the representative case to evaluate the safety of water quality and quantify the minimum safe distance after city expansion. The results reveal a safe distance of 4.6km between the two cities in 2020, which could ensure the water quality of the control section in upper reaches of Ma'anshan (the downstream city) to meet the class II standard of surface water. However, compared with the year 2010, the water quality of the control section will decline, where the COD concentration is projected to increase by 29.2% and NH3-N by 23.2%. In order to ensure the water function of the control section, the minimum safe distance between the two cities needs to be 3.2km.

    safe distance between cities;urban growth boundary;BP neural network;water quality model

    X321

    A

    1000-6923(2016)06-1905-08

    賈 寧(1992-),女,山西太原人,清華大學(xué)博士研究生,主要研究方向?yàn)閼?zhàn)略環(huán)境影響評(píng)價(jià).發(fā)表論文1篇.

    2015-12-10

    國(guó)家自然科學(xué)基金資助項(xiàng)目(71473148)

    猜你喜歡
    馬鞍山建成區(qū)間距
    宋代君臣殿上間距考論
    基于VBANUI夜間燈光調(diào)節(jié)指數(shù)的河南省城市擴(kuò)張研究
    大眾科技(2022年2期)2022-04-28 08:19:50
    馬鞍山鄭蒲港新區(qū)
    江淮法治(2022年2期)2022-03-17 08:54:16
    2013—2018年金普新區(qū)建成區(qū)空間形態(tài)擴(kuò)張規(guī)律
    高速公路指揮中心小間距LED應(yīng)用探討
    基于Landsat的南通市建成區(qū)擴(kuò)展研究
    成自瀘高速馬鞍山隧道機(jī)電工程維護(hù)淺析
    “詩(shī)城”馬鞍山 魅力黃梅戲
    算距離
    日韩中字成人| 国产欧美另类精品又又久久亚洲欧美| 啦啦啦中文免费视频观看日本| 水蜜桃什么品种好| av专区在线播放| 亚洲第一区二区三区不卡| 少妇被粗大猛烈的视频| 亚洲aⅴ乱码一区二区在线播放| 成年人午夜在线观看视频 | 大陆偷拍与自拍| 91aial.com中文字幕在线观看| 久久久精品欧美日韩精品| 午夜福利在线观看吧| 国产一级毛片在线| 免费看光身美女| 只有这里有精品99| 观看美女的网站| 午夜激情福利司机影院| 国产欧美日韩精品一区二区| 亚洲精品亚洲一区二区| 国产精品日韩av在线免费观看| 99久国产av精品国产电影| 中文在线观看免费www的网站| 网址你懂的国产日韩在线| 自拍偷自拍亚洲精品老妇| 久久久久久国产a免费观看| 一本久久精品| 亚洲av成人精品一区久久| 成年版毛片免费区| 国产成人91sexporn| 午夜福利视频精品| 久久久久久久国产电影| 九草在线视频观看| 亚洲欧洲国产日韩| 精品久久久精品久久久| 十八禁国产超污无遮挡网站| 永久免费av网站大全| 一级毛片 在线播放| 亚洲精品影视一区二区三区av| 女的被弄到高潮叫床怎么办| 国产精品.久久久| 精华霜和精华液先用哪个| 亚洲av一区综合| 国产亚洲一区二区精品| kizo精华| 免费高清在线观看视频在线观看| 日本免费a在线| 成年免费大片在线观看| 九草在线视频观看| 亚洲av成人精品一二三区| 日日干狠狠操夜夜爽| 超碰97精品在线观看| 国产午夜精品论理片| 只有这里有精品99| 久久精品久久精品一区二区三区| 女的被弄到高潮叫床怎么办| 国产有黄有色有爽视频| 欧美精品国产亚洲| 97精品久久久久久久久久精品| 啦啦啦韩国在线观看视频| 成年女人看的毛片在线观看| 日本猛色少妇xxxxx猛交久久| av女优亚洲男人天堂| 精品国内亚洲2022精品成人| 亚洲精品乱码久久久久久按摩| 精品久久久久久久久av| 国产亚洲一区二区精品| 亚洲av二区三区四区| 久久精品国产自在天天线| 日韩成人av中文字幕在线观看| 韩国高清视频一区二区三区| 欧美激情久久久久久爽电影| 我的女老师完整版在线观看| 在现免费观看毛片| 天堂网av新在线| 麻豆久久精品国产亚洲av| 国产男人的电影天堂91| 久久久久性生活片| 国产毛片a区久久久久| 日韩中字成人| 国内精品宾馆在线| 日韩在线高清观看一区二区三区| 国产极品天堂在线| 超碰av人人做人人爽久久| 亚洲精品中文字幕在线视频 | 久久精品综合一区二区三区| 成人二区视频| 全区人妻精品视频| 欧美精品国产亚洲| 干丝袜人妻中文字幕| 99re6热这里在线精品视频| 久久国产乱子免费精品| 97热精品久久久久久| 久久久久久久久久久丰满| 久久久久久久久久久免费av| 国产 亚洲一区二区三区 | 久久午夜福利片| 乱码一卡2卡4卡精品| 国产av码专区亚洲av| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 大香蕉久久网| 哪个播放器可以免费观看大片| 欧美zozozo另类| 欧美成人午夜免费资源| 高清午夜精品一区二区三区| 五月玫瑰六月丁香| 韩国av在线不卡| 日日摸夜夜添夜夜添av毛片| 国产午夜精品久久久久久一区二区三区| 日韩成人av中文字幕在线观看| 午夜精品国产一区二区电影 | 国产免费一级a男人的天堂| 有码 亚洲区| 在线免费观看不下载黄p国产| 91精品一卡2卡3卡4卡| 久久99蜜桃精品久久| 日本爱情动作片www.在线观看| 九九爱精品视频在线观看| 久久久久久久久久人人人人人人| xxx大片免费视频| 简卡轻食公司| 国产成人精品福利久久| 中文精品一卡2卡3卡4更新| 久99久视频精品免费| 在线播放无遮挡| 亚洲熟女精品中文字幕| 久久国内精品自在自线图片| 啦啦啦韩国在线观看视频| 国产av在哪里看| 最近中文字幕高清免费大全6| 亚洲精品乱码久久久v下载方式| 欧美日韩综合久久久久久| 五月玫瑰六月丁香| 成人二区视频| 三级国产精品欧美在线观看| 午夜亚洲福利在线播放| 又黄又爽又刺激的免费视频.| 精品一区二区三区人妻视频| 天美传媒精品一区二区| 老师上课跳d突然被开到最大视频| 亚洲精品aⅴ在线观看| 日韩强制内射视频| 久久精品国产鲁丝片午夜精品| 久久这里有精品视频免费| 91久久精品电影网| 一个人免费在线观看电影| 日本三级黄在线观看| 亚洲成色77777| 国产精品人妻久久久影院| 水蜜桃什么品种好| 亚洲最大成人中文| 国产成人91sexporn| 精品一区二区三区视频在线| 免费看光身美女| 色吧在线观看| 日本一二三区视频观看| 亚洲在久久综合| 国产伦精品一区二区三区视频9| 大香蕉97超碰在线| 国产亚洲精品av在线| 精品国产三级普通话版| 最近的中文字幕免费完整| 99re6热这里在线精品视频| 99热网站在线观看| 亚洲一区高清亚洲精品| 国产午夜福利久久久久久| 综合色av麻豆| 国产色爽女视频免费观看| 久久久亚洲精品成人影院| 嫩草影院新地址| 高清在线视频一区二区三区| 男的添女的下面高潮视频| 国产高清国产精品国产三级 | 免费在线观看成人毛片| 亚洲经典国产精华液单| 91久久精品国产一区二区三区| 中文字幕免费在线视频6| 九九在线视频观看精品| 丰满乱子伦码专区| 老师上课跳d突然被开到最大视频| 日韩av在线免费看完整版不卡| 亚洲aⅴ乱码一区二区在线播放| 天堂网av新在线| 免费看a级黄色片| 在线观看免费高清a一片| 成人综合一区亚洲| 久久精品夜夜夜夜夜久久蜜豆| 一本一本综合久久| 日韩,欧美,国产一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 国产亚洲一区二区精品| 精品一区二区三区人妻视频| 亚洲婷婷狠狠爱综合网| 欧美+日韩+精品| 日韩欧美三级三区| 国产亚洲一区二区精品| 成人毛片a级毛片在线播放| 久久综合国产亚洲精品| 观看免费一级毛片| 日韩中字成人| 国产视频内射| 纵有疾风起免费观看全集完整版 | 亚洲精华国产精华液的使用体验| 国产欧美另类精品又又久久亚洲欧美| 欧美3d第一页| 男的添女的下面高潮视频| 一区二区三区四区激情视频| 国产午夜精品一二区理论片| 五月玫瑰六月丁香| 亚洲国产精品国产精品| 久久国产乱子免费精品| 少妇人妻精品综合一区二区| 免费观看av网站的网址| 日日摸夜夜添夜夜爱| 日韩一区二区视频免费看| 婷婷色综合大香蕉| 非洲黑人性xxxx精品又粗又长| 网址你懂的国产日韩在线| 色综合站精品国产| 久久久a久久爽久久v久久| 中文欧美无线码| 亚洲av免费高清在线观看| 十八禁网站网址无遮挡 | av在线亚洲专区| 亚洲综合色惰| 精品99又大又爽又粗少妇毛片| 在线观看av片永久免费下载| 国产在视频线精品| 精品酒店卫生间| 欧美另类一区| 伊人久久国产一区二区| 91av网一区二区| 久久久久精品性色| 欧美日韩综合久久久久久| 美女xxoo啪啪120秒动态图| 真实男女啪啪啪动态图| 国产 一区精品| 在现免费观看毛片| h日本视频在线播放| 男女国产视频网站| 精品久久久噜噜| av免费观看日本| 亚洲国产精品成人久久小说| av.在线天堂| 天美传媒精品一区二区| 国产精品不卡视频一区二区| 国产精品三级大全| videos熟女内射| 国产午夜精品久久久久久一区二区三区| 一本一本综合久久| 欧美日韩一区二区视频在线观看视频在线 | 91aial.com中文字幕在线观看| 在现免费观看毛片| 你懂的网址亚洲精品在线观看| 国产精品精品国产色婷婷| 中文在线观看免费www的网站| 欧美日韩在线观看h| 国产老妇伦熟女老妇高清| 亚洲经典国产精华液单| 亚洲欧美日韩卡通动漫| 国产探花极品一区二区| 美女xxoo啪啪120秒动态图| 国产色婷婷99| 黑人高潮一二区| 18禁动态无遮挡网站| 久久久久久久久中文| 午夜老司机福利剧场| 国产成人精品一,二区| 亚洲欧美日韩无卡精品| 99热6这里只有精品| 亚洲成人av在线免费| 嫩草影院精品99| 精品人妻一区二区三区麻豆| 哪个播放器可以免费观看大片| 成人性生交大片免费视频hd| 全区人妻精品视频| 男女视频在线观看网站免费| 一个人看视频在线观看www免费| 精品久久久久久久久av| 精品亚洲乱码少妇综合久久| 久久精品久久久久久噜噜老黄| 国产精品一区www在线观看| av福利片在线观看| 啦啦啦中文免费视频观看日本| 日韩一区二区视频免费看| 亚洲伊人久久精品综合| 成人亚洲精品av一区二区| 日韩一本色道免费dvd| 国产在视频线精品| 精品不卡国产一区二区三区| 嫩草影院入口| 网址你懂的国产日韩在线| 18禁在线播放成人免费| 亚洲精华国产精华液的使用体验| 亚洲欧洲日产国产| 啦啦啦韩国在线观看视频| www.色视频.com| 嘟嘟电影网在线观看| 22中文网久久字幕| 好男人在线观看高清免费视频| 国产白丝娇喘喷水9色精品| 听说在线观看完整版免费高清| 大又大粗又爽又黄少妇毛片口| 亚洲国产精品专区欧美| 国产91av在线免费观看| 精品久久久久久久久亚洲| 一级黄片播放器| 大话2 男鬼变身卡| 一本一本综合久久| 亚洲欧美一区二区三区黑人 | 国产伦在线观看视频一区| 一边亲一边摸免费视频| 99热这里只有是精品50| 国产精品国产三级国产av玫瑰| 久久久久久久国产电影| 男人舔女人下体高潮全视频| 亚洲av不卡在线观看| 日本与韩国留学比较| 国产黄色视频一区二区在线观看| 午夜福利网站1000一区二区三区| 大香蕉97超碰在线| 午夜福利在线观看免费完整高清在| 国产免费视频播放在线视频 | 午夜福利视频精品| 免费无遮挡裸体视频| 久久综合国产亚洲精品| 伦精品一区二区三区| 色综合站精品国产| 亚洲人成网站在线观看播放| a级毛片免费高清观看在线播放| 午夜精品一区二区三区免费看| 国产成人免费观看mmmm| 国产在线一区二区三区精| 精品久久久久久久久亚洲| 日本黄色片子视频| 亚洲精品,欧美精品| 精品人妻熟女av久视频| 国内少妇人妻偷人精品xxx网站| 美女cb高潮喷水在线观看| 三级经典国产精品| 欧美成人a在线观看| 日韩国内少妇激情av| 激情 狠狠 欧美| 美女被艹到高潮喷水动态| 亚洲精品乱久久久久久| 国产成人精品一,二区| 精品一区二区免费观看| 亚洲精品一区蜜桃| 亚洲av二区三区四区| 婷婷色综合大香蕉| 亚洲最大成人av| 久久久久九九精品影院| 99热这里只有是精品在线观看| 免费看光身美女| 深爱激情五月婷婷| 国产成人一区二区在线| 久久综合国产亚洲精品| 久久精品国产亚洲av涩爱| 偷拍熟女少妇极品色| 国产精品不卡视频一区二区| 久久综合国产亚洲精品| 天天一区二区日本电影三级| 国精品久久久久久国模美| 久久久久久久久久成人| 国产av国产精品国产| 久久久久久久久久成人| 国产高清有码在线观看视频| www.av在线官网国产| 69人妻影院| 人妻一区二区av| 一级毛片电影观看| 久久久久免费精品人妻一区二区| 一级二级三级毛片免费看| 欧美日本视频| 2022亚洲国产成人精品| 国产精品国产三级国产专区5o| 日韩av在线免费看完整版不卡| 51国产日韩欧美| 男女边摸边吃奶| 乱系列少妇在线播放| 欧美精品一区二区大全| 99热6这里只有精品| freevideosex欧美| 女人被狂操c到高潮| 亚洲自拍偷在线| 亚洲欧美一区二区三区国产| 亚洲精品乱码久久久久久按摩| 国产人妻一区二区三区在| 亚洲欧美精品自产自拍| 国产老妇伦熟女老妇高清| 美女脱内裤让男人舔精品视频| 性色avwww在线观看| 九九在线视频观看精品| 男人和女人高潮做爰伦理| 亚洲第一区二区三区不卡| 亚洲四区av| 身体一侧抽搐| 在线免费观看不下载黄p国产| 中文天堂在线官网| 少妇人妻精品综合一区二区| 国产午夜精品久久久久久一区二区三区| 成人鲁丝片一二三区免费| 黄片无遮挡物在线观看| 女的被弄到高潮叫床怎么办| www.色视频.com| 国产大屁股一区二区在线视频| 菩萨蛮人人尽说江南好唐韦庄| 欧美日韩亚洲高清精品| 国产不卡一卡二| 精品99又大又爽又粗少妇毛片| 免费观看a级毛片全部| a级毛色黄片| 午夜福利在线在线| 国产单亲对白刺激| 男人舔女人下体高潮全视频| 中文在线观看免费www的网站| 黑人高潮一二区| 青春草视频在线免费观看| 九色成人免费人妻av| 少妇的逼水好多| 美女xxoo啪啪120秒动态图| 国产探花极品一区二区| 亚洲婷婷狠狠爱综合网| 日本免费在线观看一区| av女优亚洲男人天堂| 精品人妻偷拍中文字幕| 女人被狂操c到高潮| 欧美bdsm另类| 丝袜喷水一区| 国产精品三级大全| 国产精品99久久久久久久久| 大片免费播放器 马上看| 精品国产一区二区三区久久久樱花 | 三级经典国产精品| 国内精品宾馆在线| 在线免费观看的www视频| 亚洲天堂国产精品一区在线| 国产成年人精品一区二区| 啦啦啦韩国在线观看视频| 淫秽高清视频在线观看| 尾随美女入室| 天天躁日日操中文字幕| 肉色欧美久久久久久久蜜桃 | 一级毛片aaaaaa免费看小| 亚洲欧洲日产国产| 熟妇人妻久久中文字幕3abv| 91精品伊人久久大香线蕉| 精品国产露脸久久av麻豆 | 精品国产露脸久久av麻豆 | 人体艺术视频欧美日本| 人妻少妇偷人精品九色| 国产真实伦视频高清在线观看| 亚洲精品乱码久久久v下载方式| 免费不卡的大黄色大毛片视频在线观看 | 日韩精品有码人妻一区| av黄色大香蕉| 精品久久久久久久久av| 十八禁网站网址无遮挡 | 男的添女的下面高潮视频| 久久99热6这里只有精品| 男女下面进入的视频免费午夜| 最近最新中文字幕大全电影3| 亚洲国产欧美人成| 黄色配什么色好看| 亚洲av免费高清在线观看| 国产黄色免费在线视频| 久久久久久久大尺度免费视频| 七月丁香在线播放| 国产午夜精品一二区理论片| 久久综合国产亚洲精品| 黄色配什么色好看| 亚洲av免费高清在线观看| 亚洲精品久久久久久婷婷小说| 中文字幕制服av| 亚州av有码| 国产精品人妻久久久久久| 最近中文字幕高清免费大全6| 插逼视频在线观看| 秋霞在线观看毛片| 免费播放大片免费观看视频在线观看| 国产色爽女视频免费观看| 精品一区在线观看国产| 国产毛片a区久久久久| 国产亚洲精品久久久com| 婷婷色麻豆天堂久久| 欧美 日韩 精品 国产| 国产成人免费观看mmmm| 免费看不卡的av| 久久韩国三级中文字幕| 男的添女的下面高潮视频| 免费看日本二区| 我的女老师完整版在线观看| 亚洲欧美成人精品一区二区| 国产精品久久视频播放| 亚洲激情五月婷婷啪啪| 熟妇人妻不卡中文字幕| 麻豆国产97在线/欧美| 欧美人与善性xxx| 久久久久九九精品影院| av在线观看视频网站免费| 别揉我奶头 嗯啊视频| 久久人人爽人人爽人人片va| 天天一区二区日本电影三级| 日本欧美国产在线视频| 国产午夜精品久久久久久一区二区三区| 精品久久久久久久久久久久久| 一边亲一边摸免费视频| 岛国毛片在线播放| 亚洲欧美中文字幕日韩二区| 国产一区二区在线观看日韩| 亚洲美女搞黄在线观看| 国产亚洲午夜精品一区二区久久 | 黄片wwwwww| 国产精品麻豆人妻色哟哟久久 | 国产精品1区2区在线观看.| 嫩草影院入口| 天天躁夜夜躁狠狠久久av| 欧美 日韩 精品 国产| 日日撸夜夜添| 联通29元200g的流量卡| av在线播放精品| 亚洲av日韩在线播放| 国产单亲对白刺激| 亚洲av福利一区| 国产伦在线观看视频一区| 亚洲精品国产av成人精品| av在线观看视频网站免费| av黄色大香蕉| 欧美性感艳星| 亚洲国产高清在线一区二区三| 韩国av在线不卡| 亚洲av男天堂| 国产午夜精品论理片| 婷婷色综合www| 黄片无遮挡物在线观看| 亚洲在线自拍视频| 国产精品女同一区二区软件| 熟女电影av网| 婷婷色av中文字幕| 亚洲av福利一区| 久久鲁丝午夜福利片| 亚洲四区av| 97人妻精品一区二区三区麻豆| 黄色欧美视频在线观看| 免费观看无遮挡的男女| 国产乱人视频| 亚洲精品一区蜜桃| 欧美性猛交╳xxx乱大交人| 高清午夜精品一区二区三区| 女人被狂操c到高潮| 精品午夜福利在线看| 国产一区二区三区综合在线观看 | 丰满少妇做爰视频| 亚洲av国产av综合av卡| 亚洲图色成人| 国产视频首页在线观看| 国产精品久久久久久久电影| 亚洲一区高清亚洲精品| 午夜激情久久久久久久| 久久久久久久久久久免费av| 少妇熟女欧美另类| 亚洲av男天堂| 内射极品少妇av片p| 中文乱码字字幕精品一区二区三区 | 狂野欧美白嫩少妇大欣赏| 亚洲伊人久久精品综合| 国产成人一区二区在线| 五月天丁香电影| 精品人妻偷拍中文字幕| 久久99热6这里只有精品| av免费观看日本| 99久久人妻综合| 在线免费观看的www视频| 日日摸夜夜添夜夜添av毛片| 国内精品宾馆在线| 久久精品久久久久久噜噜老黄| 欧美变态另类bdsm刘玥| 丰满少妇做爰视频| 少妇人妻一区二区三区视频| 成人欧美大片| 午夜老司机福利剧场| 亚洲欧美日韩无卡精品| 可以在线观看毛片的网站| 午夜福利在线观看吧| 国产午夜精品久久久久久一区二区三区| 亚洲在线观看片| 女人十人毛片免费观看3o分钟| 男女边摸边吃奶| 久久久a久久爽久久v久久| 两个人视频免费观看高清| 亚洲激情五月婷婷啪啪| 国产色爽女视频免费观看| 久久97久久精品| 天天躁夜夜躁狠狠久久av| 亚洲天堂国产精品一区在线| 嫩草影院新地址| 男女啪啪激烈高潮av片| 中文字幕人妻熟人妻熟丝袜美| 青青草视频在线视频观看| 一区二区三区高清视频在线| 亚洲国产精品国产精品| 日本免费a在线| 午夜福利高清视频| 日日啪夜夜爽| 一夜夜www| 亚洲精品视频女| 丝袜喷水一区| 亚洲国产精品成人综合色| 日日啪夜夜撸| 欧美97在线视频| 国模一区二区三区四区视频| 美女大奶头视频| 午夜免费男女啪啪视频观看| 亚洲四区av| 日韩制服骚丝袜av|