• <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ī)城”馬鞍山 魅力黃梅戲
    算距離
    成人午夜高清在线视频| 波野结衣二区三区在线| 亚洲丝袜综合中文字幕| 色5月婷婷丁香| 天堂网av新在线| 可以在线观看毛片的网站| 赤兔流量卡办理| 国产av麻豆久久久久久久| 人妻久久中文字幕网| 高清毛片免费看| 一级av片app| 99热6这里只有精品| 麻豆av噜噜一区二区三区| 成人一区二区视频在线观看| 我的老师免费观看完整版| avwww免费| 男女之事视频高清在线观看| 成人鲁丝片一二三区免费| 中文字幕av在线有码专区| 插逼视频在线观看| 免费不卡的大黄色大毛片视频在线观看 | 亚洲人成网站在线播| 欧美3d第一页| 91精品国产九色| 国产久久久一区二区三区| 免费看光身美女| 午夜福利在线在线| 日韩欧美免费精品| 特大巨黑吊av在线直播| 伦理电影大哥的女人| 美女大奶头视频| 俄罗斯特黄特色一大片| 特大巨黑吊av在线直播| 亚洲专区国产一区二区| 亚洲综合色惰| 男人的好看免费观看在线视频| 99热这里只有精品一区| 小说图片视频综合网站| 99热只有精品国产| 国内精品宾馆在线| 成人一区二区视频在线观看| 久久久午夜欧美精品| 色哟哟·www| 在线免费十八禁| 亚洲三级黄色毛片| 又粗又爽又猛毛片免费看| 亚洲熟妇熟女久久| 亚洲欧美中文字幕日韩二区| 免费av观看视频| 中文字幕av成人在线电影| 有码 亚洲区| 亚洲国产色片| 在线播放无遮挡| 久久久久久九九精品二区国产| 亚洲精品日韩av片在线观看| 国产成人a区在线观看| 国产又黄又爽又无遮挡在线| 中文在线观看免费www的网站| 精品国内亚洲2022精品成人| 日韩精品中文字幕看吧| ponron亚洲| 国产蜜桃级精品一区二区三区| 国产精品国产三级国产av玫瑰| 亚洲内射少妇av| 性插视频无遮挡在线免费观看| 天天躁日日操中文字幕| 日日摸夜夜添夜夜爱| 亚洲aⅴ乱码一区二区在线播放| 少妇猛男粗大的猛烈进出视频 | 日本 av在线| 97超视频在线观看视频| 97在线视频观看| 日本爱情动作片www.在线观看 | 日本成人三级电影网站| 91午夜精品亚洲一区二区三区| 国产不卡一卡二| 日本精品一区二区三区蜜桃| 亚洲精品粉嫩美女一区| 校园人妻丝袜中文字幕| 简卡轻食公司| 欧美高清性xxxxhd video| 亚洲成av人片在线播放无| 午夜福利在线观看免费完整高清在 | 中文字幕熟女人妻在线| 久久天躁狠狠躁夜夜2o2o| 成年女人毛片免费观看观看9| 美女被艹到高潮喷水动态| 熟女电影av网| 亚洲性久久影院| 又粗又爽又猛毛片免费看| 日韩欧美三级三区| 成人午夜高清在线视频| 99热6这里只有精品| 久久久久久久亚洲中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 波多野结衣巨乳人妻| 小说图片视频综合网站| 2021天堂中文幕一二区在线观| 欧美最新免费一区二区三区| 日日摸夜夜添夜夜添av毛片| 久久人妻av系列| 99久久九九国产精品国产免费| 夜夜看夜夜爽夜夜摸| 亚洲精品成人久久久久久| 女生性感内裤真人,穿戴方法视频| 国产av麻豆久久久久久久| 久久久久久久久久久丰满| 日本三级黄在线观看| 97超碰精品成人国产| 啦啦啦韩国在线观看视频| 一级av片app| 国产精品伦人一区二区| 免费黄网站久久成人精品| 国产av麻豆久久久久久久| 国产午夜精品论理片| 久久精品91蜜桃| 美女免费视频网站| a级毛片免费高清观看在线播放| 一级毛片我不卡| 在线看三级毛片| 99riav亚洲国产免费| 哪里可以看免费的av片| 国产真实乱freesex| 天堂网av新在线| .国产精品久久| 久久久午夜欧美精品| 淫秽高清视频在线观看| 国产黄色视频一区二区在线观看 | 五月伊人婷婷丁香| 九九热线精品视视频播放| 免费观看的影片在线观看| 啦啦啦啦在线视频资源| 亚洲图色成人| 啦啦啦韩国在线观看视频| 午夜福利高清视频| 一进一出好大好爽视频| 男女之事视频高清在线观看| 一区二区三区免费毛片| 桃色一区二区三区在线观看| 丝袜喷水一区| 免费一级毛片在线播放高清视频| 精品熟女少妇av免费看| 欧美日韩精品成人综合77777| 国产私拍福利视频在线观看| 成人午夜高清在线视频| 中文亚洲av片在线观看爽| 婷婷精品国产亚洲av| 男女做爰动态图高潮gif福利片| 搡老妇女老女人老熟妇| 久久久国产成人精品二区| 少妇熟女欧美另类| 在现免费观看毛片| 一级毛片我不卡| 一本精品99久久精品77| 日韩,欧美,国产一区二区三区 | 99热全是精品| 国产精品久久久久久久久免| 老师上课跳d突然被开到最大视频| 久久久久久久久中文| 人人妻人人看人人澡| 免费观看在线日韩| 日韩三级伦理在线观看| 日韩欧美免费精品| 夜夜夜夜夜久久久久| 91狼人影院| 国产精品福利在线免费观看| 国产一区二区激情短视频| 亚洲欧美日韩高清专用| 国产大屁股一区二区在线视频| 欧美不卡视频在线免费观看| 国产午夜精品久久久久久一区二区三区 | 99久久中文字幕三级久久日本| 国产真实伦视频高清在线观看| 国产一区二区在线av高清观看| 在线国产一区二区在线| 精品人妻偷拍中文字幕| 黄片wwwwww| 欧美在线一区亚洲| 精品久久久久久成人av| 俄罗斯特黄特色一大片| 国产片特级美女逼逼视频| 亚洲在线观看片| 久久精品91蜜桃| 亚洲精品色激情综合| 欧美在线一区亚洲| 免费看av在线观看网站| 午夜日韩欧美国产| 日韩人妻高清精品专区| 变态另类丝袜制服| 六月丁香七月| 校园春色视频在线观看| 欧美激情国产日韩精品一区| 神马国产精品三级电影在线观看| 亚洲精品粉嫩美女一区| 精品欧美国产一区二区三| 亚洲性久久影院| 91麻豆精品激情在线观看国产| 少妇人妻精品综合一区二区 | 久久久久久久久久成人| 色吧在线观看| 又爽又黄无遮挡网站| 欧美日韩一区二区视频在线观看视频在线 | 少妇被粗大猛烈的视频| av女优亚洲男人天堂| 免费看光身美女| 国产男人的电影天堂91| 内地一区二区视频在线| 神马国产精品三级电影在线观看| 国模一区二区三区四区视频| 久久久久久久久中文| 日本 av在线| 亚洲成人久久爱视频| 少妇高潮的动态图| 舔av片在线| 国产综合懂色| 亚洲精品成人久久久久久| 国产欧美日韩一区二区精品| 亚洲经典国产精华液单| 偷拍熟女少妇极品色| 观看免费一级毛片| 日韩中字成人| 最新在线观看一区二区三区| 欧美三级亚洲精品| av在线播放精品| 欧美潮喷喷水| 国产精品一及| 国产69精品久久久久777片| 夜夜看夜夜爽夜夜摸| 午夜亚洲福利在线播放| 国产精品人妻久久久久久| www.色视频.com| 久久久色成人| 成年版毛片免费区| 听说在线观看完整版免费高清| 色av中文字幕| 国产精品1区2区在线观看.| 精品国产三级普通话版| 日本黄色视频三级网站网址| 亚洲一级一片aⅴ在线观看| 久久久久国内视频| 亚洲人成网站在线播| av天堂中文字幕网| 3wmmmm亚洲av在线观看| 亚洲七黄色美女视频| 亚洲av电影不卡..在线观看| 中国美白少妇内射xxxbb| 国产一区二区亚洲精品在线观看| 国产美女午夜福利| 人妻制服诱惑在线中文字幕| 草草在线视频免费看| а√天堂www在线а√下载| 国产精品久久电影中文字幕| 亚洲欧美成人精品一区二区| 在线观看av片永久免费下载| 十八禁网站免费在线| 春色校园在线视频观看| 国产在视频线在精品| 波野结衣二区三区在线| 亚洲成人精品中文字幕电影| 我要搜黄色片| 99热全是精品| 成人高潮视频无遮挡免费网站| 国产蜜桃级精品一区二区三区| 日日干狠狠操夜夜爽| 国产精品一二三区在线看| 国产亚洲精品综合一区在线观看| 免费av毛片视频| 国产美女午夜福利| 亚洲在线自拍视频| 国产av在哪里看| 99九九线精品视频在线观看视频| 少妇熟女欧美另类| 18+在线观看网站| 12—13女人毛片做爰片一| 男女边吃奶边做爰视频| 少妇熟女欧美另类| 国产亚洲精品av在线| 亚洲欧美日韩无卡精品| 欧美激情在线99| 天堂影院成人在线观看| 午夜影院日韩av| 能在线免费观看的黄片| 日本黄大片高清| 一个人观看的视频www高清免费观看| 亚洲欧美成人精品一区二区| 91久久精品国产一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 搡老岳熟女国产| 国产亚洲av嫩草精品影院| 99精品在免费线老司机午夜| 丰满人妻一区二区三区视频av| 麻豆乱淫一区二区| 少妇高潮的动态图| 免费观看人在逋| 赤兔流量卡办理| 黑人高潮一二区| 亚洲精品乱码久久久v下载方式| 久久人人精品亚洲av| 亚洲成人久久爱视频| 免费观看在线日韩| 禁无遮挡网站| 成年女人看的毛片在线观看| 最近手机中文字幕大全| 在线观看66精品国产| av在线播放精品| 尾随美女入室| 99九九线精品视频在线观看视频| 成人精品一区二区免费| 国产欧美日韩精品亚洲av| 中文在线观看免费www的网站| 久久精品久久久久久噜噜老黄 | 日本爱情动作片www.在线观看 | 美女内射精品一级片tv| 看片在线看免费视频| 欧美日本亚洲视频在线播放| 精品一区二区免费观看| 12—13女人毛片做爰片一| 国国产精品蜜臀av免费| 亚洲一区高清亚洲精品| 精品福利观看| 天堂√8在线中文| 日韩 亚洲 欧美在线| 伊人久久精品亚洲午夜| 日韩欧美一区二区三区在线观看| 熟女人妻精品中文字幕| 不卡视频在线观看欧美| 日日撸夜夜添| 女的被弄到高潮叫床怎么办| 91av网一区二区| 国产亚洲欧美98| 国产女主播在线喷水免费视频网站 | 高清午夜精品一区二区三区 | 国产国拍精品亚洲av在线观看| 欧美色欧美亚洲另类二区| av在线亚洲专区| 亚洲美女视频黄频| 麻豆成人午夜福利视频| 日韩国内少妇激情av| 国产不卡一卡二| 一级毛片我不卡| 亚洲精华国产精华液的使用体验 | 亚洲自偷自拍三级| 狂野欧美白嫩少妇大欣赏| 少妇高潮的动态图| 亚洲国产欧美人成| 国产精品久久电影中文字幕| 国产av一区在线观看免费| 国产亚洲精品久久久久久毛片| 一个人免费在线观看电影| 99精品在免费线老司机午夜| 色综合亚洲欧美另类图片| 性欧美人与动物交配| av在线蜜桃| 亚洲七黄色美女视频| av在线老鸭窝| 大型黄色视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲三级黄色毛片| 亚洲国产精品sss在线观看| 久久久国产成人免费| 成人高潮视频无遮挡免费网站| 九九热线精品视视频播放| 天堂av国产一区二区熟女人妻| 波多野结衣高清作品| 99热这里只有是精品50| 床上黄色一级片| 99热精品在线国产| a级毛片免费高清观看在线播放| 国产精品野战在线观看| 中文字幕熟女人妻在线| 搞女人的毛片| 搡老妇女老女人老熟妇| 免费无遮挡裸体视频| 男女视频在线观看网站免费| 波野结衣二区三区在线| 国产精品久久视频播放| 精品一区二区三区视频在线观看免费| 亚洲四区av| 国内精品美女久久久久久| 在现免费观看毛片| 亚洲激情五月婷婷啪啪| a级毛色黄片| 久久人人精品亚洲av| 欧美3d第一页| 日韩精品有码人妻一区| 乱码一卡2卡4卡精品| 国产精品国产三级国产av玫瑰| 十八禁国产超污无遮挡网站| 深夜精品福利| 国产在线男女| 欧美性猛交╳xxx乱大交人| 91狼人影院| 国产黄色小视频在线观看| 亚洲无线观看免费| av视频在线观看入口| 免费无遮挡裸体视频| 久久天躁狠狠躁夜夜2o2o| 国产精品亚洲一级av第二区| 久久九九热精品免费| 国产在线精品亚洲第一网站| 国产欧美日韩精品亚洲av| 啦啦啦韩国在线观看视频| 女生性感内裤真人,穿戴方法视频| 啦啦啦韩国在线观看视频| 国产私拍福利视频在线观看| 大又大粗又爽又黄少妇毛片口| 变态另类丝袜制服| 天堂av国产一区二区熟女人妻| 国产精品,欧美在线| 欧美另类亚洲清纯唯美| 五月伊人婷婷丁香| 亚洲一区高清亚洲精品| 久久久久久九九精品二区国产| 国产成人a∨麻豆精品| 亚洲无线观看免费| 一进一出好大好爽视频| 国产亚洲欧美98| 悠悠久久av| 97热精品久久久久久| 欧美高清成人免费视频www| 久久久久久伊人网av| 成年免费大片在线观看| 亚洲成人精品中文字幕电影| 高清毛片免费看| 色吧在线观看| 看片在线看免费视频| 国产一区二区三区在线臀色熟女| 看免费成人av毛片| 亚洲va在线va天堂va国产| 亚洲欧美精品自产自拍| 美女 人体艺术 gogo| 亚洲自拍偷在线| 精品久久久久久成人av| 亚洲av免费在线观看| 日韩一区二区视频免费看| 黄色视频,在线免费观看| 日韩国内少妇激情av| 亚洲五月天丁香| 国产午夜福利久久久久久| 可以在线观看的亚洲视频| 亚洲精品在线观看二区| 欧美极品一区二区三区四区| 99热网站在线观看| 国产av麻豆久久久久久久| 18禁在线播放成人免费| 一个人看的www免费观看视频| 亚洲va在线va天堂va国产| 日韩欧美 国产精品| 一级黄片播放器| www日本黄色视频网| 欧美zozozo另类| 成人高潮视频无遮挡免费网站| 一区二区三区免费毛片| 精品一区二区三区人妻视频| 麻豆精品久久久久久蜜桃| 婷婷色综合大香蕉| 最近最新中文字幕大全电影3| 有码 亚洲区| 久久久久免费精品人妻一区二区| 毛片一级片免费看久久久久| 婷婷精品国产亚洲av| 国产精品嫩草影院av在线观看| 男人舔奶头视频| 精品日产1卡2卡| 亚洲综合色惰| 亚洲三级黄色毛片| 我要搜黄色片| 熟女电影av网| 少妇熟女欧美另类| 国产精品av视频在线免费观看| 国产高清三级在线| 最近视频中文字幕2019在线8| 不卡视频在线观看欧美| 最近手机中文字幕大全| 人人妻,人人澡人人爽秒播| 成熟少妇高潮喷水视频| 日韩欧美三级三区| 中文资源天堂在线| 久久天躁狠狠躁夜夜2o2o| 国产女主播在线喷水免费视频网站 | 久久久久久久久久久丰满| 午夜福利18| 欧美丝袜亚洲另类| 亚洲国产精品合色在线| 91在线精品国自产拍蜜月| 亚洲av五月六月丁香网| 一个人观看的视频www高清免费观看| 看免费成人av毛片| 日本爱情动作片www.在线观看 | 亚洲国产欧美人成| 观看免费一级毛片| 国产不卡一卡二| 女的被弄到高潮叫床怎么办| 美女免费视频网站| 成年女人毛片免费观看观看9| 色在线成人网| 又爽又黄无遮挡网站| 午夜福利18| 人妻久久中文字幕网| 久久精品国产亚洲av香蕉五月| 国产精品永久免费网站| 亚洲欧美成人综合另类久久久 | 老司机福利观看| 久久中文看片网| 91在线观看av| 午夜老司机福利剧场| 欧美+亚洲+日韩+国产| 精品99又大又爽又粗少妇毛片| 十八禁网站免费在线| 欧美激情国产日韩精品一区| 国产精华一区二区三区| 国产午夜精品久久久久久一区二区三区 | 色综合色国产| 午夜爱爱视频在线播放| 一a级毛片在线观看| 国产色婷婷99| 91在线精品国自产拍蜜月| 国产一级毛片七仙女欲春2| 99热这里只有是精品在线观看| 国产伦精品一区二区三区视频9| 亚洲三级黄色毛片| 两个人视频免费观看高清| 国内精品美女久久久久久| 欧美中文日本在线观看视频| 国产精品免费一区二区三区在线| av在线亚洲专区| 久久久久九九精品影院| 色在线成人网| 久久久精品大字幕| 日本爱情动作片www.在线观看 | 国产视频内射| 丰满人妻一区二区三区视频av| 中文资源天堂在线| 国产精华一区二区三区| 国产真实伦视频高清在线观看| 丝袜喷水一区| 国产黄色小视频在线观看| 不卡一级毛片| 成熟少妇高潮喷水视频| 国产精品三级大全| 婷婷精品国产亚洲av| 国产中年淑女户外野战色| www.色视频.com| 小蜜桃在线观看免费完整版高清| 18禁在线无遮挡免费观看视频 | 亚洲精品乱码久久久v下载方式| 国产一区二区亚洲精品在线观看| 黄色日韩在线| 久久久久性生活片| 久久久久久久久久黄片| 伊人久久精品亚洲午夜| 97热精品久久久久久| 国产一级毛片七仙女欲春2| 亚洲成人久久爱视频| 亚洲成人精品中文字幕电影| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美三级三区| 99久久精品一区二区三区| 欧美一区二区亚洲| 亚洲经典国产精华液单| 99久久精品国产国产毛片| 22中文网久久字幕| 国产日本99.免费观看| 日韩中字成人| 亚洲天堂国产精品一区在线| 日韩欧美三级三区| 人妻丰满熟妇av一区二区三区| 久久精品综合一区二区三区| 91在线观看av| 欧美激情在线99| 欧美区成人在线视频| 1024手机看黄色片| 此物有八面人人有两片| 国产v大片淫在线免费观看| 久久这里只有精品中国| 国产在线精品亚洲第一网站| 国产成人福利小说| 国产精品电影一区二区三区| 成人特级黄色片久久久久久久| 麻豆国产av国片精品| 日韩人妻高清精品专区| 晚上一个人看的免费电影| 国产精品无大码| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品成人久久久久久| 级片在线观看| 国产欧美日韩精品亚洲av| 国产亚洲精品av在线| 男女边吃奶边做爰视频| 两性午夜刺激爽爽歪歪视频在线观看| 午夜老司机福利剧场| 大型黄色视频在线免费观看| 精品福利观看| 最近在线观看免费完整版| 国产精品福利在线免费观看| 亚洲国产精品久久男人天堂| 天堂影院成人在线观看| 色综合站精品国产| 婷婷亚洲欧美| 国产淫片久久久久久久久| 在线观看av片永久免费下载| 久久综合国产亚洲精品| 熟女电影av网| 可以在线观看的亚洲视频| 久久精品人妻少妇| 中文字幕免费在线视频6| 97超视频在线观看视频| 中文资源天堂在线| 亚洲高清免费不卡视频| 国产精品久久久久久久电影| 老司机影院成人| 美女大奶头视频| 淫妇啪啪啪对白视频| 最近视频中文字幕2019在线8| 日韩精品中文字幕看吧|