999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于GIS的婁底市河流廊道寬度確定方法研究

2017-09-03 10:02:39程燦輝許乙青
福建質量管理 2017年10期
關鍵詞:功能方法模型

程燦輝 許乙青

(湖南大學 湖南 長沙 410000)

基于GIS的婁底市河流廊道寬度確定方法研究

程燦輝 許乙青

(湖南大學 湖南 長沙 410000)

根據河流廊道的特征和結構,研究河流廊道寬度確定方法;以湖南省婁底市為例,根據100年一遇的洪水淹沒范圍、行洪斷面和現有堤防范圍,綜合確定河道和河漫灘寬度,硬化的河道部分寬度被嚴重縮小,因此需大致確定其自然河道范圍,采用GIS分析硬化河道緩沖帶和整條河流自然河道緩沖帶,對兩者進行特征比較,選出與硬化河段特征大致相同的自然河道緩沖帶,通過剖面大致確定硬化河道自然邊界;同時,根據河流廊道不同的特征和問題以及土地利用情況,將河流廊道分成3段,以500m一個單位,將各段河流分成若干段,通過CAD建立每點的剖面圖和高層數據,以GIS緩沖分析和寬度與坡度關系式的結合確定各個單位河岸帶點的位置,連接各點得出河岸帶的位置,以緩沖區為界,超出的點可省略。

河流廊道;寬度;GIS定量分析;土地利用

隨著社會的發展,城市的生態環境受到嚴重破壞,城市內水體污染嚴重,河流大面積被填埋,整體河流廊道受到人為干預和破壞,尤其在城市中心區段。河流廊道的建立,可以保護洪乏區以及濱河濕地,提高河流的洪水調蓄能力;對水生和濱水棲息地生物有恢復和保護作用;對水質的保護,以及可控制泥沙和農業區的面源污染;提高環境的舒適性和休閑娛樂價值。但這些功能的實現都需要給予河流廊道一定的寬度,才能更好的發揮。

一、河流廊道的概念

河流廊道的概念出于景觀生態學理論。Forman主要從空間結構來定義,認為河流廊道是指沿河流分布而不同于基質的植被帶;從河流廊道的功能和價值方面,Gardiner和Cole認為河流廊道是河流及因河流的存在而具有現有的或者潛在價值的土地,價值包括休閑娛樂、自然資源的保護和環境的友好等;從河流廊道的連續性方面,Petersen認為河流與周圍之間的地帶,稱為河岸帶或緩沖帶,流域有一定長度時,才得稱為河流廊道。河流廊道的定義較多。本文綜合各個對河流廊道的定義,針對研究主題,將河流廊道定義為:由河岸帶、河道和河漫灘組成,具有一定連續性和寬度的線性廊道。

二、河流廊道的結構組成和功能

(一)河流廊道的結構組成

從景觀生態學的角度研究,河流廊道的結構主要由五部分組成:河道、河漫灘、邊坡、河堤,部分高地。

河流廊道有著各種各樣的形態,但是大部分的河流廊道由三部分組成:河道、河漫灘和河岸帶。河流廊道的寬度主要由這三部分決定,河道和河漫灘的寬度決定著其對洪水的緩沖能力,河岸帶(也叫河岸植被帶)的寬度決定著是否能形成足夠的內部環境,為陸生動物提供遷徙通道,這也決定著物種的豐富度。并且在一定程度上能截留污染物。

(二)河流廊道的功能

城市河流廊道總體來說具備兩大塊的功能,一個社會廊道功能,另一個自然廊道功能。從社會廊道功能出發分析,它包括經濟、歷史文化遺產、景觀等功能,景觀功能可表現為城市中的濱水公共空間可以給城市民眾提供休閑娛樂場地,歷史文化遺產功能表現為城市的歷史往往也反應河流歷史,城市河流兩岸往往也坐落著歷史性的建筑和古跡,另外對于經濟功能來說,濱水住宅、休閑商業場所、娛樂文化等大大的促進了城市的經濟發展。

從自然廊道功能分析, 首先,河岸帶對城市熱島效應有一定的抑制作用,植被通過蒸騰作用調節周圍小氣候,這稱之為氣候調節功能。河岸帶對水凈化、水土流失、降低噪聲污染等有著明顯的效益,這可稱之為其防護功能。地下水的補給、保持河流流量、對污染物的截留和凈化,以及養分循環,這叫做河流廊道的物質功能。對河流的防洪是每一個濱水城市必做的一項措施,而河流廊道中的河岸帶對城市的防洪起著重要的作用,這稱為河流廊道的防洪功能。

河流廊道的功能還有很多,如維持生物多樣性,農田灌溉等,從核心來講河流廊道功能主要是社會廊道功能和自然廊道功能。這些功能的發揮都與河流廊道的寬度息息相關,如果沒有適當的寬度,這些功能也不能得到很好的發揮。因此河流廊道寬度的合理確定也應成為一個重要的研究內容。

三、河流廊道寬度確定現狀研究

就目前來說,關于定量確定河流廊道寬度的方法很少,因為涉及到眾多的因素,比如鄰近河流的土地利用狀況、水質狀況、周邊植被生長狀況、土壤地質狀況等,因此要將其量化需要花費大量的精力。

目前要確定河流廊道各部分的寬度,主要根據其所承載的各種功能綜合確定,由于河道和河漫灘內部密切的水文和生物聯系,是一個不可分割的部分,所以將其看作一個整體,共同確定其寬度。另一部分則是河岸帶寬度的確定。這一部分寬度的確定較為靈活和復雜,大部分對其范圍的確定只是通過經驗,沒有一個統一的計算標準,并且,針對不同的河流廊道和周邊環境等因素有不同的確定方法。

(一)河道和河漫灘的寬度確定方法

河道和河漫灘的主要功能是,水生和濕地生物的保護以及對洪水的控制。寬度可根據行洪要求和天然河道的游蕩范圍較明確的確定,通常通過洪水分析或者采用現有的標準,主要有:通過100年一遇的洪水淹沒范圍、當地防洪規劃確定的行洪斷面、現有堤防范圍和恢復天然游蕩河道所需要的范圍等。

本文的案例中中心區段的河道完全被硬化,現有河道邊界遠小于自然邊界,因此需要借助GIS分析,首先分析出硬化段河流的緩沖區,同時也建立河流自然段的緩沖區,通過對比分析,大致確定硬化段河流的邊界。

(二)河岸帶寬度確定方法

第一種方法是基于數學模型的河岸緩沖帶寬度確定方法,主要有河岸帶生態管理模型(簡稱REMM);農田管理系統化合物、徑流及侵蝕模型(簡稱CREAMS);植被過濾帶模型(簡稱VFSMOD)。

河岸帶生態管理模型需要四類數據才能運行,分別是河岸帶土壤數據、氣候數據、田間數據以及植被數據,通過這些數據加上若干指標可以模擬河岸緩沖帶的水文特征,不同緩沖帶的寬度、坡度、土壤和植被類型對水質的影響;所以,當植被土壤類型和坡度大小固定時,可通過REMM模型的模擬找到最優的緩沖帶寬度。

農田管理系統化合物、徑流及侵蝕模型是在20世紀70年代由美國農業部農業研究所建立的,是基于田塊尺度的模型,它主要包括三個模塊:水文模塊、侵蝕模塊和化合物模塊,水文模塊根據降雨量計算徑流量,進而通過水文模型輸出的結果和降雨量的大小、土壤可蝕性計算侵蝕模塊的侵蝕量,最后結合前兩者輸出的結果及化合物的特性估算化合物模塊的化合物流失量,三者是相互聯系相互作用的過程。該模型根據植被、土壤和降雨的特征,結合它的三個模型來模擬田塊的徑流、侵蝕和農用化合物流失過程,評價不同耕作措施對面源污染負荷的影響,以此來確定河岸緩沖帶寬度和污染物去除率之間的關系。

植被過濾帶模型是計算河岸緩沖帶對坡面地表徑流中泥沙凈化效果的數學模型,該模型能夠計算緩沖帶的出流水量、下滲水量和泥沙截留效果。不同寬度的緩沖帶泥沙過濾效果不一樣。

上述方法都是基于物理過程的模擬模型,模擬整個污染物沉積的過程,從而推出對于預期污染物去除效率所需要的河岸緩沖帶寬度。

第二種方法基于回歸方程的應用,定量確定河岸帶寬度,主要有Phillips水文模型、Phillips時間模型、Nieswand模型、Mander模型、SWAT模型。Phillips水文模型只考慮河岸帶攔截地標徑流中的沉積物和污染物的物理過程,除污效率低;Phillips時間模型綜合考慮了地表徑流和地下徑流的水文過程,既考慮了攔截吸附等物理過程,也考慮土壤微生物的化學效率,除污效果較好。Nieswand模型主要考慮寬度與坡面坡度關系,當坡度大于15%時,緩沖帶寬度需要加寬。Mander模型的計算,很大成分偏向經驗,比如模型中的坡長因子,就是一個純粹的經驗值。SWAT模型需經過現場實測,并用實測數據校正不同植被、寬度和坡度下回歸方程中的變化常數、寬度變化指數和坡度變化指數。

另外還有夏繼紅等總結的統計回歸模型計算法、動態基質的數學模型計算法以及綜合優化模型計算法,通過建立河岸帶寬度和坡度、氮素、坡面流速和時間、泥沙沉降速度等因素的數學關系式,以定量的方式確定河岸帶寬度。

第三種方法是基于GIS的分析方法,通過建立流域內的水文、地形、土壤特征、土地利用和土地覆蓋信息數據庫等一系列數據,以上述回歸方程為模型,把模型需要的數據疊加到同一個圖層上,用GIS得出緩沖帶寬度,然后把矢量轉化成柵格形式得到河岸帶范圍。這種方法目前的運用較少,GIS定量計算河岸帶寬度的方法逐漸成為重要的工具,并且針對不同的環境和條件,GIS的運用也不同。

本文在前人研究的基礎上,運用GIS定量分析和定量公式結合,研究婁底市河流廊道河岸帶的寬度問題。同時針對婁底市河流的問題和現狀,確定河岸帶寬度的方法。主要的創新點在于先用GIS確定河岸帶的緩沖區,再結合Mayer等針對河岸帶的坡度和寬度的定量公式,算出所選點的河岸帶寬度,以緩沖區作為最大控制范圍,連接各個點得出河岸帶邊界。

四、婁底市河流廊道寬度確定研究

(一)研究區概況

婁底市位于湖南省中部,其地理位置東經110°45′40″~112°31′05″,北緯27°12′28″~ 28°14′27″。轄區總面積為8117 km2,婁底市區位于市域東部。

婁底市境內河流密布,水系較多,主要河流有:東部漣水河,為湘江中游一大支流,源于新邵觀音山,自西向東,流經漣源市、婁星區、雙峰縣,經湘鄉至湘潭縣河口入湘江,境內全長85.85 km,沿途納孫水河、湄江河、測水河等1~4級支流89條,控制流域面積3906 km2。西部資水河,由南向北,流經冷水江、新化,經安化柘溪,過益陽注入洞庭湖,貫穿境內西半部,境內流程112 km,有1~4級支流100條,控制流域面積3985 km2。

根據婁底市水系概況可知,婁底境內的主干河流(或稱為一級河流)為漣水河,由西向東橫跨婁底市市區,它是婁底市最具有代表性的河流廊道,因此,以漣水河為研究對象,確定其河流廊道寬度計算的方法。

根據3中所述,將整個河流廊道寬度分成兩部分來確定,以河道和河漫灘作為一部分,河岸帶作為另外一部分來研究。

(二)確定河道和河漫灘寬度的方法

河道也稱為河床,是指河流流經的路線,它的寬度就是自然河流兩邊邊界點的長度,但是婁底市河流的某些河段被填過,導致有些河段并不是自然河流的寬度,尤其在土地價值較大的地方,如市中心區地段,大部分河段河渠被硬化,河道的自然邊界線嚴重被縮小。因此要采用適當的方法確定河道自然邊界線,具體操作是:在地形圖中找到這一段被硬化的河道,以500m作為一個單位等分,將這一段導入GIS建立緩沖區。然后可得到每一段等分后的河流廊道的平面橫斷面剖面線,利用CAD中GeoSection工具得出每段剖面線的剖面圖及剖面圖上的高層點數據。將這些數據導入GIS中。

將未被硬化的天然河段地形圖導入GIS,按同樣的距離建立緩沖區。在其中找尋一段與被硬化河道緩沖區相似特征的,對比分析。并在CAD中標注這一段天然河道,也同樣以500m為一個等分單位,得出剖面圖和高層點數據。大致確定每個點的自然邊界線位置。然后進一步對比分析,進而大致得出硬化河道的自然邊界線。

其他的河段按照現有寬度計算,確定完河道寬度后,進一步綜合100年一遇的洪水淹沒范圍、當地防洪規劃確定的行洪斷面、現有堤防范圍。

圖一 河段對比示意圖

(三)確定河岸帶寬度的方法

河岸帶寬度確定是最為復雜和最為關鍵的,河岸帶具有許多至關重要的功能,它的確定要考慮許多因素,在前面已經詳細明述其各個設計方法,但不同的方法適用于不同的河岸帶環境,并不是任何一種方法都會適用于所有的條件,針對不同的環境,確定河岸帶寬度的方法也不一樣。

婁底市漣水河貫穿湖南省婁底市和湘鄉市境內,婁星區內漣水河位于城區中部,為漣鋼和中心城區的分界線,兩岸用地以居住和商業為主,現狀只有婁底一大橋至二大橋北岸建設了防洪堤,南岸部分地段建設防洪墻,按50年一遇洪水位標準建設。

經實地調查,漣水河每一段呈現的問題不同,比如中心區段主要問題是水體的污染和建筑沿岸而建,而在周邊的郊區地段,主要的問題在于河流廊道寬度太小,甚至有些河流直接被填埋,水體污染較小,在郊野地段,雖然水質較好,但河岸帶寬度較小,兩邊的植被單一雜亂。因此,根據問題分析以及漣水河不同河段的問題特點,可將其分成三段,針對不同的問題確定其河流廊道寬度計算方法(見下表)。

區段具體地段問題中心區段漣鋼大橋—婁底大橋—梅子灣大橋水質污染嚴重,氮磷含量超標郊區段漣鋼大橋—漣鋼北大橋—婁漣公路大橋;梅子灣大橋—吉星大橋—江龍灘橋大量河床被填埋,河流河岸帶被占用甚至荒廢,致使河流廊道斷裂郊野段婁漣公路大橋—石井鎮的黃泥塘街道江龍灘橋—大埠橋河岸帶植被單一,滿足不了生物多樣性的要求,提供良好的陸生生物棲息地

1.中心區段河岸帶寬度確定方法研究

中心區段河岸帶以降低周邊氮磷為主要目標,達到改善水質的目的。同時提高河岸帶環境質量,提供休閑娛樂場地。并且,中心城區用地較為緊張,很多建筑沿岸而建,拆遷難度較大,因此,河岸帶的寬度受到一定限制。對于限定在建筑與河漫灘之間的河岸帶范圍,要達到降低氮磷的效果,需要一些靈活的方式。

首先,可通過定量公式計算中心區段河流廊道河岸帶的理想寬度,將其與實際被限定的河岸帶范圍進行比較,一般來說,理論寬度會大于實際河岸帶范圍,因此,需要改變其他因素來彌補所缺少的河岸帶寬度,如通過坡度的改變、植被的配置等。

氮素去除率和河岸帶寬度的回歸關系式:WRZ=(15.24 S1/2),RZN=39.5WRZ0.1644;S:河岸平均坡度;WRZ:河岸帶寬度;RRN:氮素去除率。從公式可知,河岸帶寬度和坡度以及與氮素去除率成正比,即當寬度越大時,也就是水流通過河岸帶的時間越長,氮素去除率越大。就婁底中心區段而言,只能通過改變植物種植密度或植被覆蓋度來降低水流的速度。

因此,河岸帶寬度就是建筑后退之后的邊界到河漫灘的距離,寬度雖達不到要求,但可通過在有限的范圍里達到降低氮磷的效果,即通過增加植被密度。

圖二 區段示意圖

2.郊野段河岸帶寬度確定方法研究

郊野區段以滿足河岸帶生物棲息地為目標,河岸帶擴展范圍較為自由,周邊基本以農田和山體為主,因此可根據GIS定量確定。

首先對郊野段的河流廊道做緩沖區分析,緩沖區距離設為3km(因河岸帶寬度基本在這個范圍以內),得到緩沖區的目的是為了給河岸帶寬度一個上限值。得到緩沖區距離后,以500m為單位將郊野段河流廊道分為若干段,將等分后的圖導入CAD中與地形圖中的等高線疊加,這樣就可得到每一段等分后河流廊道的平面橫斷面剖面線,進而利用CAD中GeoSection工具得出每段剖面線的剖面圖及剖面圖上的高層點數據。

得到剖面圖和高層點數據之后,計算每個剖面的坡度,在一定程度上為了更好削減氮磷的含量,同時截留和吸收一些污染物。因此可根據公式WRZ=(15.24 S1/2)計算出寬度,將每個計算的寬度值標在橫斷面剖面上與高層點值相對應,最后將對應的點數據導入GIS中,建立河岸帶邊界點圖層,與所得出的緩沖區疊加。超過緩沖區范圍的點以緩沖區邊界的點替代。將在緩沖區內的點連接成線,即為河岸帶的范圍。

圖三 河岸帶邊界選取示意圖

3.郊區段河岸帶寬度確定方法研究

郊區段以提高對地表和地下水的保護功能以及提供休憩空間為目標和基礎,盡量拓寬河岸帶寬度。這一段不存在上一段土地利用緊張的問題,但是仍然有一些建筑緊鄰河岸而建,大多數河段的土地利用較為充裕。對于建筑緊鄰河岸的河段,采用與中心區段同樣的方法,即通過改變坡度和增加植被密度。

對于土地利用較為充裕的地段采用與郊野段同樣的方法,即通過GIS和CAD確定范圍,主要的區別在于,郊區段河流廊道緩沖區的范圍不同,因為周邊有一些待開發用地,盡量減少對待開發土地的肆意利用,因此緩沖區的范圍設為800m,滿足基本寬度要求。

五、結論

確定河流廊道的寬度要考慮的因素眾多,并且不同的河流廊道特征所涉及的因素會不同,如本文分析的婁底市,周邊土地的利用對其寬度有決定性的影響,同時還有植被狀況,以及土壤地質特征等。

本文利用河流廊道范圍確定的,相關研究結合GIS技術,研究婁底市河流廊道寬度的確定方法。河道和河漫灘寬度的確定通過百年一遇的洪水淹沒范圍和行洪斷面以及現有堤防綜合確定,被硬化的河道部分,利用GIS分析硬化河道緩沖帶和整條河流自然河道緩沖帶,對比分析得出結果,這種方法存在一定的缺陷,依據GIS分析之后的緩沖區特點,找尋與硬化河道相同特征的自然河段并不容易,并且有些偏向主觀性;河岸帶寬度根據所要建立的河岸帶功能確定,功能不同,河岸帶的寬度也不同。在確定河岸帶功能后,采用定量公式或者GIS的方法綜合確定。但是針對不同環境和不同特征的河流廊道,采用的定量公式也會有所差異。

[1]孫金偉,許文盛.河岸植被緩沖帶生態功能及其過濾機理的研究進展[J].長江科學院院報,2017,(03):40-44.

[2]陳星.城市濱水緩沖帶景觀設計研究[D].安徽農業大學,2016.

[3]蒙倩彬.基于生物多樣性保護的城市生態廊道研究[D].北京林業大學,2016.

[4]尚琴.景觀生態規劃中GIS技術的應用[D].西安建筑科技大學,2014.

[5]李莉,姜允芳.國內外城市河流綠道的理論與實踐研究進展[J].中國人口.資源與環境,2014,(S1):309-312.

[6]侯利萍.河湖濱岸緩沖帶保護紅線劃定技術研究[D].河北師范大學,2013.

[7]韓路,王海珍,于軍.河岸帶生態學研究進展與展望[J].生態環境學報,2013,(05):879-886

[8]趙警衛.河岸帶景觀結構、功能及其關系研究[D].華東師范大學,2013.

[9]侯利萍,何萍,錢金平,賈俊艷,徐杰.河岸緩沖帶寬度確定方法研究綜述[J].濕地科學,2012,(04):500-506.

[10]肖紅霞.河岸緩沖帶生態護岸模式劃分及景觀設計[D].北京林業大學,2012.

[11]王家生,孔麗娜,林木松,閔鳳陽.河岸帶特征和功能研究綜述[J].長江科學院院報,2011,(11):28-35.

[12]李靜羽.基于3S技術的城市綠色河流廊道景觀格局分析——以廣州市番禺區為例[J].吉林農業,2011,(02):132-133.

[13]張東旭,郭晉平.我國河岸帶生態學的研究進展[J].山西林業科技,2010,(04):26-30.

[14]閻麗鳳.河岸緩沖帶對氮、磷污染物的去除效果研究[D].遼寧工程技術大學,2011.

[15]孫汝斌.城市化過程中河流景觀的生態優化途徑[D].河南農業大學,2008.

[16]呂海燕,李政海,李建東,宋國寶.廊道研究進展與主要研究方法[J].安徽農業科學,2007,(15):4480-4482+4484.

[17]吳亞偉.基于生態理念的城市綠化空間系統建構初探[D].西安建筑科技大學,2007.

[18]蒙倩彬.基于生物多樣性保護的城市生態廊道研究[D].北京林業大學,2016.

[19]尚琴.景觀生態規劃中GIS技術的應用[D].西安建筑科技大學,2014.

[20]李莉,姜允芳.國內外城市河流綠道的理論與實踐研究進展[J].中國人口.資源與環境,2014,(S1):309-312.

[21]侯利萍.河湖濱岸緩沖帶保護紅線劃定技術研究[D].河北師范大學,2013.

[22]李靜羽.基于3S技術的城市綠色河流廊道景觀格局分析——以廣州市番禺區為例[J].吉林農業,2011,(02):132-133.

[23]袁思.甘淑,吳文強,彭文啟.運用GIS技術確定河岸帶范圍的方法研究[D].人民黃河,2013.

[24]夏繼紅,鞠蕾,林俊強,陳永明,王為林,韓玉玲,胡玲.河岸帶事宜寬度要求與確定方法[D].河海大學學報(自然科學版),2013.

程燦輝(1992.12-),男,漢族,江西景德鎮人,碩士研究生,湖南大學建筑學院,丘陵城市生態規劃;許乙青(1965.11-),男,漢族,湖南長沙人,碩士生導師,高級工程師,湖南大學建筑學院,丘陵城市規劃設計。

猜你喜歡
功能方法模型
一半模型
也談詩的“功能”
中華詩詞(2022年6期)2022-12-31 06:41:24
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
關于非首都功能疏解的幾點思考
3D打印中的模型分割與打包
用對方法才能瘦
Coco薇(2016年2期)2016-03-22 02:42:52
四大方法 教你不再“坐以待病”!
Coco薇(2015年1期)2015-08-13 02:47:34
捕魚
中西醫結合治療甲狀腺功能亢進癥31例
主站蜘蛛池模板: 亚洲黄色视频在线观看一区| 亚洲综合18p| 色哟哟色院91精品网站| 亚洲午夜国产精品无卡| 中文字幕免费播放| 日韩免费中文字幕| 亚洲精品国产精品乱码不卞| 亚洲第一成年人网站| 亚洲欧洲自拍拍偷午夜色| 成人免费午夜视频| 无码在线激情片| 日本91在线| 91毛片网| 欧美人人干| 国产综合色在线视频播放线视| 久久国产亚洲欧美日韩精品| 亚洲码一区二区三区| 97精品伊人久久大香线蕉| 青青草原偷拍视频| 亚洲精品第一页不卡| 99久久精品久久久久久婷婷| 国产高颜值露脸在线观看| 黄色福利在线| 国产欧美精品一区aⅴ影院| 国产人妖视频一区在线观看| 国产成人盗摄精品| 在线看片中文字幕| 秋霞国产在线| 久草视频中文| AV天堂资源福利在线观看| 伊人AV天堂| 国产精品刺激对白在线| 成人午夜天| 国产一级精品毛片基地| 69av免费视频| 国产SUV精品一区二区| 999福利激情视频| 夜夜爽免费视频| 国产第一页第二页| 青青草一区| 国产乱人伦精品一区二区| 亚洲无码久久久久| 欧美午夜视频| 无码国产伊人| 国产专区综合另类日韩一区| 福利片91| 色噜噜狠狠色综合网图区| 亚洲国产一成久久精品国产成人综合| 美女无遮挡拍拍拍免费视频| 国产欧美视频综合二区| 免费无码又爽又黄又刺激网站| 内射人妻无码色AV天堂| 国产美女精品人人做人人爽| 九九九精品成人免费视频7| 中国精品自拍| 国产丝袜啪啪| 青青热久免费精品视频6| 亚洲第一视频网| 亚洲欧洲日韩综合色天使| 992tv国产人成在线观看| 波多野结衣第一页| 国产91色在线| 久综合日韩| 无码福利日韩神码福利片| 99热这里只有免费国产精品| 欧美精品高清| 动漫精品啪啪一区二区三区| 亚洲欧美精品在线| 国产黄网站在线观看| 99精品一区二区免费视频| 波多野吉衣一区二区三区av| 久久96热在精品国产高清| 亚洲精品无码AⅤ片青青在线观看| 亚洲中文字幕23页在线| AV天堂资源福利在线观看| 亚洲AⅤ波多系列中文字幕| 啦啦啦网站在线观看a毛片| 亚洲精品欧美日韩在线| 女人天堂av免费| 性色一区| 国产亚洲成AⅤ人片在线观看| 亚洲AV无码久久精品色欲|