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

    GMS軟件技術(shù)在原油泄漏對地下水污染狀況模擬評價中的應(yīng)用

    2017-12-20 10:49:06,,
    地下水 2017年5期
    關(guān)鍵詞:運(yùn)移水文地質(zhì)含水層

    ,,

    (海南省地質(zhì)調(diào)查院,海南 海口 570206)

    GMS軟件技術(shù)在原油泄漏對地下水污染狀況模擬評價中的應(yīng)用

    陜寧,楊永鵬,王曉林

    (海南省地質(zhì)調(diào)查院,海南 海口 570206)

    當(dāng)前中國對能源的需求日益增加,油品儲罐也向著大型化方向發(fā)展,大型的油罐裝置如果發(fā)生原油泄漏會對地下水造成嚴(yán)重污染。本文以某石化項(xiàng)目煉油廠油罐爆炸事故為例,結(jié)合煉油廠周邊的水文地質(zhì)條件,運(yùn)用GMS軟件建立地下水?dāng)?shù)值模型,研究事故工況下原油泄漏在地下水環(huán)境中的濃度變化情況,分析并預(yù)測地下水在一定時間內(nèi)受污染的情況,對原油泄漏對項(xiàng)目區(qū)地下水的污染進(jìn)行定量評價。模擬結(jié)果表明,在風(fēng)險事故工況下,以20 a為最長預(yù)測期可以發(fā)現(xiàn)原油污染的影響范圍為312.2×104m2,運(yùn)移距離為1 736.19 m,平均擴(kuò)散速率為15.61×104m2/a,平均運(yùn)移速率為86.81m/a,事故發(fā)生后,污染物進(jìn)入地下水含水層中的濃度逐漸升高,污染影響范圍隨著時間推移逐步擴(kuò)大,污染物長期平均擴(kuò)散與運(yùn)移速率會小于短期平均速率,原油在預(yù)測期內(nèi)均有超標(biāo)現(xiàn)象,污染事故具有長期而緩慢的危害性。

    環(huán)境影響評價;地下水;數(shù)值模擬;GMS軟件技術(shù);原油泄漏

    隨著中國經(jīng)濟(jì)的快速發(fā)展,對能源的需求日益增加,油品儲罐也向著大型化方向發(fā)展。大型的油罐裝置如果發(fā)生故障,可能帶來災(zāi)難性的后果。所以,評價大型油罐區(qū)的污染情況十分重要。本文以具體煉油廠油罐爆炸事故為例,采用GMS軟件,依據(jù)《環(huán)境影響評價技術(shù)導(dǎo)則—地下水環(huán)境》(HJ610-2016)[1],結(jié)合煉油廠周邊的水文地質(zhì)條件,對原油泄漏對地下水的污染進(jìn)行定量預(yù)測評價,可為類似工程提供參考。

    1 背景概況

    1.1 基本情況

    該項(xiàng)目區(qū)屬熱帶季風(fēng)海洋性氣候,全年較暖。受海洋影響,冬季少雨,夏季多雨。熱帶氣旋影響比較頻繁。年平均氣溫為24.4℃。5~10月為豐水期,總降水量為1 893.2 mm,全年降雨時間分布不均勻,7月、8月降雨偏多。

    1.2 水文地質(zhì)概況

    模擬區(qū)地形起伏較小,地貌由火山巖臺地、濱海平原組成,另外由于人類填海造地工程活動形成一種新的用地類型,暫稱為人工填海區(qū)地貌。

    區(qū)域地下水根據(jù)其賦存條件、含水巖組的巖性、水力特征,劃分為第四系松散巖類孔隙潛水、玄武巖裂隙-孔洞潛水、松散-半固結(jié)巖類孔隙承壓水和基巖裂隙水四大類。第四系松散巖類孔隙潛水含水層在整個區(qū)域均有分布,巖性為含礫粗砂、中細(xì)砂、細(xì)砂、粘土質(zhì)砂等,階地前緣較厚,后緣較薄。玄武巖裂隙水含水層分布于新英灣以北,巖性為灰、深灰-灰黑色玄武巖,巖石裂隙、氣孔發(fā)育,除德義嶺、三都一帶地表有紅土覆蓋外,其他地區(qū)玄武巖裸露。松散-半固結(jié)巖類孔隙承壓水在區(qū)域均有分布,可分為??诮M承壓含水層(組)和燈樓角組承壓含水層(組)兩個含水層(組),??诮M承壓含水層(組)巖性主要為貝殼碎屑砂礫巖,含礫中粗砂、中細(xì)砂、粘土質(zhì)砂等,砂與粘土質(zhì)砂呈松散狀,含水層孔隙度大,透水性好;燈樓角組承壓含水層巖性主要為灰綠色中粗砂、中細(xì)砂、細(xì)砂及粘土質(zhì)砂,新英灣以南,含水層頂板向北西西傾斜,新英灣以北,含水層頂板向北傾斜。

    2 地下水污染的預(yù)測與評價

    2.1 水文地質(zhì)概念模型

    2.1.1 含水層結(jié)構(gòu)概化

    由于本項(xiàng)目無重質(zhì)非水相污染物,如果污染物泄漏后主要會對潛水含水層造成污染;而對于上部存在良好隔水層的承壓含水層,則很難對其造成污染。因此,本次模型預(yù)測主要評價火山巖裂隙孔洞潛水、松散巖類孔隙潛水及與上部潛水含水層直接接觸的??诮M上部微承壓水,它們之間具有較好的水力聯(lián)系,在本次評價中統(tǒng)一概化為潛水含水組。

    2.1.2 邊界條件概化

    東北邊界:參考相關(guān)資料,模擬區(qū)東北部邊界除中段接受地下水側(cè)向徑流補(bǔ)給,概化為二類邊界的定流量邊界,其它段概化為二類邊界的零流量邊界(圖1)。

    環(huán)海邊界:整個模擬區(qū)內(nèi)地下水大致由東部大海徑流排泄,環(huán)海邊界概化為定水頭邊界(圖1)。

    圖1 模擬區(qū)邊界條件概化圖

    2.1.3 源匯項(xiàng)的確定

    模擬區(qū)的補(bǔ)給來源主要有大氣降水補(bǔ)給,地下水側(cè)向徑流補(bǔ)給;排泄方式主要為側(cè)向徑流排泄和蒸發(fā)。

    根據(jù)模擬區(qū)所在位置,地下水系統(tǒng)的內(nèi)部結(jié)構(gòu)、外部環(huán)境、邊界條件、水文地質(zhì)參數(shù)等進(jìn)行分析研究,可將模型區(qū)地下水系統(tǒng)概化為非均質(zhì)、水平各向同性、二維結(jié)構(gòu)、穩(wěn)定地下水流系統(tǒng),外部環(huán)境通過邊界進(jìn)行物質(zhì)與能量交換的二維流水文地質(zhì)概念模型。

    2.2 地下水流數(shù)值模型

    非均質(zhì)、水平各向同性、二維結(jié)構(gòu)、穩(wěn)定地下水流系統(tǒng),外部環(huán)境通過邊界進(jìn)行物質(zhì)與能量交換的二維流水文地質(zhì)概念模型,用下列偏微分方程的定解問題來描述。

    式中:K為含水層介質(zhì)的水平滲透系數(shù)(m/d);Kn為邊界面法向方向的滲透系數(shù)(m/d);h為地下水系統(tǒng)的水位標(biāo)高(m);b為潛水含水層底板標(biāo)高(m);p為潛水面的蒸發(fā)和降雨入滲強(qiáng)度等(m/d);Γ0為滲流區(qū)域的上邊界,即地下水的自由表面;h1為已知水位邊界值(m);Γ1為一類邊界,已知水位邊界;Γ2為二類邊界的定流量邊界;q為含水層二類邊界單位面積上的流量,流入為正,流出為負(fù),隔水邊界為0(m/d);Γ3為滲流區(qū)域的下邊界,即含水層底部的隔水邊界;n為邊界面的法線方向。

    2.3 模擬流場

    本次模擬以2012年3月的地下水水位作為模擬的流場。

    2.4 模擬軟件及模擬區(qū)網(wǎng)格剖分

    本次模擬采用美國環(huán)境保護(hù)局(USA EPA)開發(fā)的GMS6.5[2]。

    在模擬過程中,廠區(qū)及周邊水源地為模擬重點(diǎn)區(qū),因此,在這些地方每個單元格的邊長為50 m,其他外圍單元格邊長為200 m,總共剖分有效單元格14 195個,模型在垂向上概化為一層(圖2)。模擬區(qū)中有效觀測孔9個,可以通過觀測孔的觀測水位對模擬區(qū)的水位狀況進(jìn)行控制。

    圖2 模擬區(qū)網(wǎng)格剖分圖(剖面垂直放大25倍)

    2.5 地下水水流模型識別

    2.5.1 水文地質(zhì)參數(shù)識別

    通過對評價區(qū)地質(zhì)、水文地質(zhì)條件、地形地貌、地下水流場特征及野外抽水試驗(yàn)的計算結(jié)果的綜合分析,滲透系數(shù)K為本次模擬的主要水文地質(zhì)參數(shù)。通過模型的模擬調(diào)試,獲得校準(zhǔn)后的水文地質(zhì)參數(shù)見表1和圖3。

    2.5.2 地下水水位擬合

    圖4為模擬后地下水流場擬合情況,從圖中可以看出模型計算水位長與實(shí)際水位擬合較好。9個觀測孔計算水位與觀測值的擬合差均小于0.7 m(表2),其中除G7觀測孔外的擬合差小于0.4 m。

    表1 計算模擬參數(shù)表 m/d

    圖3 模擬區(qū)滲透系數(shù)K分區(qū)圖

    圖4 模擬區(qū)潛水含水層流場擬合圖

    2.5.3 地下水系統(tǒng)均衡分析

    通過模型識別得到地下水流數(shù)值模型的地下水水量均衡結(jié)果(表3),模擬區(qū)的地下水水量的補(bǔ)給和排泄基本平衡。

    表2 觀測水位擬合統(tǒng)計表

    表3 模擬區(qū)一個水文年潛水含水層均衡表 104·m3/a

    2.6 地下水溶質(zhì)運(yùn)移模型

    地下水溶質(zhì)運(yùn)移可通過以下方程進(jìn)行描述。

    2.6.1 控制方程

    2.6.2 初始條件

    C(x,y,z,t) = C0(x,y,z) (x,y,z)∈Ω,t = 0

    式中:C0(x,y,z)已知濃度分布;Ω為模型模擬區(qū)。

    2.6.3 邊界條件

    (1)第一類邊界-給定濃度邊界

    C(x,y,z,t)│Γ1= c (x,y,z,t) (x,y,z)∈Γ1, t ≥ 0

    式中:Γ1為給定濃度邊界;c (x,y,z,t) 為一定濃度邊界上的濃度分布。

    (2)第二類邊界-給定彌散通量邊界

    式中,Γ2為通量邊界;fi(x,y,z,t) 為Γ2邊界上已知的彌散通量函數(shù)。

    (3)第三類邊界-給定溶質(zhì)通量邊界

    式中,Γ3為混合邊界;gi(x,y,z,t) 為Γ3上已知的對流-彌散總的通量函數(shù)。

    應(yīng)用GMS 6.5中的MT3D模塊可以對以上數(shù)學(xué)模型進(jìn)行數(shù)值模擬。

    圖5 原油污染影響范圍預(yù)測圖

    污染物污染物標(biāo)準(zhǔn)/mg/L污染物入海時間/a模擬時間/a污染物最大擴(kuò)散距離/m污染物污染范圍/104m2污染物超標(biāo)范圍/104m2原油0.3161285.35113.392.145693.22168.44129.84101215.26223.52153.36201736.19312.2220.52

    2.7 預(yù)測與評價

    假定5×104m3的柴油罐的儲存量為30×104t,爆炸破壞地表防滲面積為1 000 m2,伴生二次污染事故為物料泄漏,其中95%燃燒,則地面柴油量為1.5×104t。依據(jù)廠區(qū)的水文地質(zhì)條件,柴油儲罐所在區(qū)域的包氣帶厚度約為5 m、孔隙度約為0.6,則所泄露物料的最大滲漏量為防滲層以下包氣帶的孔隙體積,即可能滲入地下水的污染物總量不超過1.5×104t。

    在調(diào)試好的地下水流場模型基礎(chǔ)上,通過GMS軟件的MT3D模塊進(jìn)行污染物遷移轉(zhuǎn)化的數(shù)值模擬計算,得到事故風(fēng)險情景下污染物運(yùn)移的預(yù)測結(jié)果。預(yù)測結(jié)果總結(jié)見表4。污染物的遷移分布見圖5,給出了風(fēng)險事故下石油類泄漏后1 a、5 a、10 a、20 a后污染物在水平方向上的運(yùn)移范圍。在風(fēng)險事故工況下,本次最長預(yù)測期20 a內(nèi),原油污染的影響范圍為312.2×104m2,運(yùn)移距離為1 736.19 m,平均擴(kuò)散速率為15.61×104m2/a,平均運(yùn)移速率為86.81 m/a。通過模擬可知,事故發(fā)生后,污染物進(jìn)入地下水,含水層中污染物濃度逐漸升高,污染影響范圍隨著時間推移逐步擴(kuò)大。原油在預(yù)測期內(nèi)均有超標(biāo)現(xiàn)象。污染物長期平均擴(kuò)散與運(yùn)移速率小于短期平均速率,因此污染的防治工作,不僅要充分考慮污染事故危害的迅速擴(kuò)展,而且還要考慮污染事故長期而緩慢的危害性。

    3 結(jié)語

    運(yùn)用GMS軟件建立某石化項(xiàng)目地下水二維數(shù)值模型,通過對含水層的概化和水文地質(zhì)參數(shù)的反復(fù)調(diào)整[3],經(jīng)過模型的識別和驗(yàn)證,模擬結(jié)果較符合實(shí)際情況,模型模擬效果較好,可以直觀說明保護(hù)目標(biāo)受到的影響程度,為類似工程的地下水環(huán)境影響評價提供參考。

    [1]環(huán)境保護(hù)部環(huán)境工程評估中心.HJ610-2016環(huán)境影響評價技術(shù)導(dǎo)則地下水環(huán)境[S].北京:中國環(huán)境科學(xué)出版社.2016.

    [2]鄭春苗, Gordon D. Benneff.孫晉玉,盧國平, 譯.地下水污染物遷移模擬[M]. 北京: 高等教育出版社.2009(2 ):186-198.

    [3]譚文清,孫春,胡婧敏,等.GMS在地下水污染質(zhì)運(yùn)移數(shù)值模擬預(yù)測中的應(yīng)用[J].東北水利水電.2008.26(5):54-55,59.

    ApplicationofGMSsoftwareinsimulationevaluationofgroundwaterpollutioncausedbycrudeoilleakage

    SHANNing,YANGYong-peng,WANGXiao-lin

    (Hainan Geological Survey Institute, Haikou, Hainan 570206)

    At present, the demand for energy in China is increasing day by day, and the oil storage tanks are developing toward large scale. Large oil tank installations can cause serious pollution to the ground water if oil spills occur. Based on a petrochemical projects refinery oil tank explosion accident as an example, combined with oil refineries around the hydrogeological conditions, using GMS software to establish a numerical model of groundwater, the oil spill accident conditions of concentration in groundwater environment changes, the analysis and prediction of groundwater polluted situation, within a certain amount of time for the oil spill pollution to the groundwater in the project area for quantitative evaluation. Simulation results show that the risk of accident conditions, with 20 years as the longest forecast period can be found that the influence of oil pollution of 312.2×104m2, migration distance of 1 736.19 m, average diffusion rate of 15.61×104m2/a, the average migration rate was 86.81 m/a, after the accident, the concentration of pollutants entering groundwater aquifer rise gradually, over time, gradually expand the pollution impact range, long-term pollutant diffusion and migration rate will be less than the short-term average rate, crude oil in the forecast period has the violation, the dangers of pollution accident has a long and slow.

    environmental impact assessment;groundwater;numerical simulation;GMS software technology;crude oil leakage

    X523

    A

    1004-1184(2017)05-0074-04

    2017-05-16

    陜寧(1985-),女,山西臨汾人,工程師,主要從事水工環(huán)地質(zhì)調(diào)查方面工作。

    猜你喜歡
    運(yùn)移水文地質(zhì)含水層
    全球多個含水層里的水正快速流失
    基于抽水試驗(yàn)計算水文地質(zhì)參數(shù)的分析與確定
    曲流河復(fù)合點(diǎn)壩砂體構(gòu)型表征及流體運(yùn)移機(jī)理
    基于GPRS實(shí)現(xiàn)自動化水文地質(zhì)抽水試驗(yàn)
    基于轉(zhuǎn)移概率的三維水文地質(zhì)結(jié)構(gòu)劃分
    東營凹陷北帶中淺層油氣運(yùn)移通道組合類型及成藏作用
    水文地質(zhì)在工程地質(zhì)勘察中的有效應(yīng)用
    美國西部奧加拉拉含水層水位下降原因初探
    開采過程中上覆急傾斜巖層運(yùn)移規(guī)律模擬與研究
    川西坳陷孝泉-新場地區(qū)陸相天然氣地球化學(xué)及運(yùn)移特征
    俄罗斯特黄特色一大片| 亚洲久久久国产精品| 国产欧美日韩一区二区三| 国产欧美亚洲国产| 又大又爽又粗| 黄色片一级片一级黄色片| 视频区图区小说| 国产麻豆69| 国产成人一区二区三区免费视频网站| av片东京热男人的天堂| 一本—道久久a久久精品蜜桃钙片| 成人国语在线视频| 午夜久久久在线观看| 午夜成年电影在线免费观看| 欧美另类亚洲清纯唯美| 成人三级做爰电影| 国产人伦9x9x在线观看| 男男h啪啪无遮挡| 欧美激情久久久久久爽电影 | 啦啦啦在线免费观看视频4| 国产人伦9x9x在线观看| 久久av网站| 国产男女超爽视频在线观看| 视频区欧美日本亚洲| 深夜精品福利| 国产av精品麻豆| 超色免费av| 在线观看免费日韩欧美大片| 国产97色在线日韩免费| 免费av中文字幕在线| 啦啦啦中文免费视频观看日本| 国产老妇伦熟女老妇高清| 麻豆成人av在线观看| 精品少妇一区二区三区视频日本电影| 精品少妇一区二区三区视频日本电影| 亚洲国产欧美日韩在线播放| av线在线观看网站| 国产又爽黄色视频| a级片在线免费高清观看视频| 亚洲精品国产精品久久久不卡| 叶爱在线成人免费视频播放| 交换朋友夫妻互换小说| 久久久久国内视频| 久久久精品国产亚洲av高清涩受| 久久久国产一区二区| 久久青草综合色| 超碰成人久久| 妹子高潮喷水视频| 亚洲国产欧美日韩在线播放| 亚洲美女黄片视频| 国产精品免费一区二区三区在线 | 日韩熟女老妇一区二区性免费视频| 制服诱惑二区| 午夜福利视频精品| 久久精品国产a三级三级三级| 亚洲精品国产精品久久久不卡| 午夜精品久久久久久毛片777| 精品国产乱子伦一区二区三区| 丰满饥渴人妻一区二区三| kizo精华| 久9热在线精品视频| 午夜免费鲁丝| av天堂在线播放| 女人久久www免费人成看片| 欧美人与性动交α欧美精品济南到| 9色porny在线观看| 久久九九热精品免费| 欧美大码av| 日韩中文字幕欧美一区二区| 女人久久www免费人成看片| 三级毛片av免费| av国产精品久久久久影院| 国产精品久久久久久精品古装| 少妇被粗大的猛进出69影院| 在线观看免费日韩欧美大片| 亚洲中文日韩欧美视频| 国产在线观看jvid| 一边摸一边抽搐一进一小说 | 亚洲精品一卡2卡三卡4卡5卡| 巨乳人妻的诱惑在线观看| 国产免费福利视频在线观看| 伦理电影免费视频| 国产视频一区二区在线看| 宅男免费午夜| 三级毛片av免费| 91九色精品人成在线观看| 新久久久久国产一级毛片| 欧美日韩亚洲综合一区二区三区_| 女人精品久久久久毛片| 久久国产精品影院| 国产三级黄色录像| 中文字幕色久视频| 在线观看免费视频日本深夜| 美女高潮到喷水免费观看| 国产成人一区二区三区免费视频网站| 亚洲一区中文字幕在线| 午夜福利一区二区在线看| 麻豆av在线久日| 久久影院123| 动漫黄色视频在线观看| 老司机影院毛片| 成人国产av品久久久| 精品少妇黑人巨大在线播放| 国产成人av教育| 日日爽夜夜爽网站| 国产真人三级小视频在线观看| 母亲3免费完整高清在线观看| 色尼玛亚洲综合影院| 黑人巨大精品欧美一区二区mp4| 人妻一区二区av| 欧美日韩视频精品一区| 女人爽到高潮嗷嗷叫在线视频| 国产精品麻豆人妻色哟哟久久| 精品福利永久在线观看| 久久影院123| 免费在线观看影片大全网站| 国产亚洲精品久久久久5区| 亚洲第一欧美日韩一区二区三区 | 国产精品av久久久久免费| 在线av久久热| 精品国产超薄肉色丝袜足j| 亚洲熟妇熟女久久| 大片电影免费在线观看免费| 亚洲精品国产区一区二| 老司机影院毛片| 老司机福利观看| 国产精品1区2区在线观看. | 18禁美女被吸乳视频| 久久久国产一区二区| 国产精品亚洲一级av第二区| 国产一区二区三区在线臀色熟女 | 18禁黄网站禁片午夜丰满| 在线观看舔阴道视频| 午夜福利欧美成人| 麻豆国产av国片精品| 咕卡用的链子| 久久中文看片网| 69精品国产乱码久久久| 欧美精品啪啪一区二区三区| 天堂动漫精品| 久久久久国产一级毛片高清牌| 久久国产精品影院| 亚洲性夜色夜夜综合| 狂野欧美激情性xxxx| 水蜜桃什么品种好| 人妻一区二区av| 视频区图区小说| 不卡av一区二区三区| 他把我摸到了高潮在线观看 | 久久久国产欧美日韩av| 精品福利永久在线观看| av天堂久久9| 亚洲国产成人一精品久久久| 69av精品久久久久久 | 99国产极品粉嫩在线观看| 天堂中文最新版在线下载| 麻豆国产av国片精品| 50天的宝宝边吃奶边哭怎么回事| 国产一区有黄有色的免费视频| 日韩人妻精品一区2区三区| 国产在线精品亚洲第一网站| 波多野结衣av一区二区av| 视频在线观看一区二区三区| 午夜91福利影院| avwww免费| 久久性视频一级片| 热99国产精品久久久久久7| kizo精华| 中文字幕色久视频| 久久久久久人人人人人| 天堂中文最新版在线下载| 精品第一国产精品| 亚洲欧美色中文字幕在线| 久久精品成人免费网站| 91九色精品人成在线观看| 两个人免费观看高清视频| 波多野结衣一区麻豆| 亚洲欧美一区二区三区久久| 欧美午夜高清在线| 十分钟在线观看高清视频www| 麻豆成人av在线观看| 成人亚洲精品一区在线观看| 最近最新中文字幕大全电影3 | 午夜福利一区二区在线看| 亚洲一卡2卡3卡4卡5卡精品中文| 757午夜福利合集在线观看| 欧美老熟妇乱子伦牲交| 精品午夜福利视频在线观看一区 | 男女下面插进去视频免费观看| 我的亚洲天堂| 极品少妇高潮喷水抽搐| 国产av又大| 咕卡用的链子| 老汉色av国产亚洲站长工具| 午夜久久久在线观看| 欧美日韩亚洲高清精品| 久久久久国内视频| 欧美国产精品va在线观看不卡| 宅男免费午夜| 精品午夜福利视频在线观看一区 | 久久久久久久久久久久大奶| 91九色精品人成在线观看| 午夜福利视频在线观看免费| 男人操女人黄网站| 老司机在亚洲福利影院| 90打野战视频偷拍视频| 久久影院123| 涩涩av久久男人的天堂| 69精品国产乱码久久久| 在线观看免费视频日本深夜| 免费观看av网站的网址| 日日摸夜夜添夜夜添小说| 美女午夜性视频免费| 久久午夜亚洲精品久久| 久热爱精品视频在线9| av有码第一页| 成人18禁高潮啪啪吃奶动态图| 亚洲伊人色综图| 视频在线观看一区二区三区| 色婷婷久久久亚洲欧美| 色94色欧美一区二区| 午夜久久久在线观看| 亚洲 国产 在线| 视频在线观看一区二区三区| 在线十欧美十亚洲十日本专区| 99精品欧美一区二区三区四区| 最新在线观看一区二区三区| 亚洲 国产 在线| 视频在线观看一区二区三区| 757午夜福利合集在线观看| 国产在线视频一区二区| 国产精品一区二区在线不卡| 午夜福利欧美成人| 成人国语在线视频| 日韩欧美国产一区二区入口| 国产欧美日韩一区二区精品| 欧美成人午夜精品| 国产成人精品无人区| 久久 成人 亚洲| 两个人免费观看高清视频| 亚洲国产av新网站| 成人av一区二区三区在线看| 国产成人影院久久av| 搡老熟女国产l中国老女人| 黄片小视频在线播放| 在线av久久热| 亚洲精品在线观看二区| 最黄视频免费看| 国产精品1区2区在线观看. | 欧美精品一区二区免费开放| 在线 av 中文字幕| 免费观看人在逋| 50天的宝宝边吃奶边哭怎么回事| 蜜桃在线观看..| 久久久久久亚洲精品国产蜜桃av| 亚洲 国产 在线| a在线观看视频网站| 国产精品久久久久成人av| 一二三四社区在线视频社区8| 欧美日韩亚洲国产一区二区在线观看 | bbb黄色大片| 亚洲欧美日韩另类电影网站| 国产成人精品无人区| 成人精品一区二区免费| 国产一区二区三区综合在线观看| 亚洲情色 制服丝袜| 丁香欧美五月| 国产97色在线日韩免费| 久久精品国产a三级三级三级| 大香蕉久久成人网| 一级黄色大片毛片| 国产成人精品无人区| 精品卡一卡二卡四卡免费| 另类亚洲欧美激情| 亚洲avbb在线观看| 国产成人系列免费观看| 日本撒尿小便嘘嘘汇集6| 激情视频va一区二区三区| 国产免费福利视频在线观看| 久久精品亚洲精品国产色婷小说| 久久精品国产亚洲av香蕉五月 | 高清在线国产一区| 久久精品aⅴ一区二区三区四区| 日本vs欧美在线观看视频| 亚洲精品在线美女| 国产精品偷伦视频观看了| 亚洲一区中文字幕在线| 亚洲av第一区精品v没综合| 国产1区2区3区精品| 亚洲成人国产一区在线观看| 正在播放国产对白刺激| 人妻 亚洲 视频| 国产区一区二久久| 精品亚洲乱码少妇综合久久| 亚洲一卡2卡3卡4卡5卡精品中文| 建设人人有责人人尽责人人享有的| 欧美精品一区二区大全| 久久精品亚洲熟妇少妇任你| 国内毛片毛片毛片毛片毛片| 天天躁日日躁夜夜躁夜夜| 日韩欧美三级三区| 搡老乐熟女国产| 国产淫语在线视频| 大码成人一级视频| 国产又爽黄色视频| 久久精品国产99精品国产亚洲性色 | 一本一本久久a久久精品综合妖精| 国产激情久久老熟女| 新久久久久国产一级毛片| 欧美日韩福利视频一区二区| 亚洲国产欧美日韩在线播放| 免费观看a级毛片全部| 一级a爱视频在线免费观看| 男女之事视频高清在线观看| 一区二区三区乱码不卡18| 五月开心婷婷网| 黄色怎么调成土黄色| 欧美午夜高清在线| 搡老乐熟女国产| 午夜福利乱码中文字幕| 国产亚洲午夜精品一区二区久久| 国产淫语在线视频| 亚洲国产中文字幕在线视频| 男女无遮挡免费网站观看| 女人精品久久久久毛片| 亚洲成人免费av在线播放| 日本一区二区免费在线视频| 少妇粗大呻吟视频| 亚洲天堂av无毛| 国产亚洲av高清不卡| 交换朋友夫妻互换小说| 777久久人妻少妇嫩草av网站| 亚洲精品中文字幕一二三四区 | 最新的欧美精品一区二区| 9热在线视频观看99| 久久精品国产亚洲av高清一级| 日韩免费高清中文字幕av| 黄色视频不卡| 在线观看免费高清a一片| 精品少妇黑人巨大在线播放| 丰满迷人的少妇在线观看| av视频免费观看在线观看| 2018国产大陆天天弄谢| 男女无遮挡免费网站观看| 中文字幕色久视频| 精品一区二区三卡| 国产精品久久电影中文字幕 | 亚洲成av片中文字幕在线观看| 搡老乐熟女国产| av视频免费观看在线观看| 欧美中文综合在线视频| 日韩视频在线欧美| av网站在线播放免费| 国产一区有黄有色的免费视频| 国产区一区二久久| 在线观看一区二区三区激情| 久久久久视频综合| 巨乳人妻的诱惑在线观看| 老司机午夜福利在线观看视频 | 国产极品粉嫩免费观看在线| 999久久久国产精品视频| 一区在线观看完整版| 丰满饥渴人妻一区二区三| 桃红色精品国产亚洲av| 另类精品久久| 亚洲天堂av无毛| 国产成人欧美| 免费在线观看黄色视频的| 纯流量卡能插随身wifi吗| 亚洲一卡2卡3卡4卡5卡精品中文| 好男人电影高清在线观看| 一进一出抽搐动态| 日韩免费av在线播放| 欧美人与性动交α欧美精品济南到| 成年女人毛片免费观看观看9 | 国产不卡av网站在线观看| 夫妻午夜视频| 精品乱码久久久久久99久播| 精品少妇久久久久久888优播| 十八禁网站网址无遮挡| 国产高清国产精品国产三级| 国产一区有黄有色的免费视频| 国产欧美日韩一区二区精品| 亚洲欧美精品综合一区二区三区| 欧美激情 高清一区二区三区| 在线观看免费午夜福利视频| 国产精品久久久久久精品电影小说| 捣出白浆h1v1| 精品欧美一区二区三区在线| 黄色 视频免费看| 久久久久久久精品吃奶| 少妇猛男粗大的猛烈进出视频| 国产主播在线观看一区二区| 首页视频小说图片口味搜索| 亚洲av成人不卡在线观看播放网| 免费人妻精品一区二区三区视频| 亚洲自偷自拍图片 自拍| 99在线人妻在线中文字幕 | 黄色a级毛片大全视频| 一个人免费看片子| 男人操女人黄网站| 97人妻天天添夜夜摸| 露出奶头的视频| 淫妇啪啪啪对白视频| 国产又爽黄色视频| 精品国产一区二区久久| 亚洲成人免费电影在线观看| 国产伦人伦偷精品视频| 国产成人免费无遮挡视频| 人妻 亚洲 视频| 国产精品久久电影中文字幕 | 19禁男女啪啪无遮挡网站| 国产1区2区3区精品| 90打野战视频偷拍视频| 国产成+人综合+亚洲专区| 男女之事视频高清在线观看| 亚洲欧美精品综合一区二区三区| 两性夫妻黄色片| 久久久久久人人人人人| 女人精品久久久久毛片| h视频一区二区三区| 高清在线国产一区| 老汉色av国产亚洲站长工具| 王馨瑶露胸无遮挡在线观看| 亚洲av日韩在线播放| a级毛片黄视频| 久久精品成人免费网站| 少妇粗大呻吟视频| 亚洲免费av在线视频| 国产成人系列免费观看| 高清在线国产一区| 啦啦啦 在线观看视频| 国产成人免费无遮挡视频| 人人妻人人澡人人爽人人夜夜| 咕卡用的链子| 夜夜爽天天搞| 国产一区二区三区视频了| 在线观看免费午夜福利视频| 91成人精品电影| 久久久久国产一级毛片高清牌| 国产主播在线观看一区二区| 精品一区二区三区av网在线观看 | 美女扒开内裤让男人捅视频| 热99久久久久精品小说推荐| 黄色怎么调成土黄色| 欧美日韩中文字幕国产精品一区二区三区 | 久久人人97超碰香蕉20202| 国产成人免费观看mmmm| 一区二区日韩欧美中文字幕| 精品国产乱子伦一区二区三区| 90打野战视频偷拍视频| 久久久久网色| 91国产中文字幕| 999久久久国产精品视频| 国产97色在线日韩免费| 考比视频在线观看| 窝窝影院91人妻| 国产精品av久久久久免费| 亚洲午夜精品一区,二区,三区| 超色免费av| 悠悠久久av| 国产区一区二久久| 国产精品国产av在线观看| 久久99热这里只频精品6学生| 国产男女超爽视频在线观看| 手机成人av网站| 亚洲七黄色美女视频| 国产亚洲精品一区二区www | 午夜久久久在线观看| 婷婷成人精品国产| av超薄肉色丝袜交足视频| 一本久久精品| 99在线人妻在线中文字幕 | 在线观看免费高清a一片| 99国产精品免费福利视频| 国产日韩一区二区三区精品不卡| 黄片小视频在线播放| 久久ye,这里只有精品| avwww免费| 亚洲综合色网址| 在线观看一区二区三区激情| 桃花免费在线播放| 亚洲专区中文字幕在线| 如日韩欧美国产精品一区二区三区| 亚洲熟女毛片儿| 久久免费观看电影| 久久久久网色| 亚洲av日韩在线播放| 老司机亚洲免费影院| 天堂中文最新版在线下载| 成人永久免费在线观看视频 | 咕卡用的链子| 国产精品自产拍在线观看55亚洲 | 欧美 日韩 精品 国产| 极品人妻少妇av视频| 日韩欧美一区视频在线观看| 日本av手机在线免费观看| 国产有黄有色有爽视频| 久热爱精品视频在线9| 日韩有码中文字幕| 精品国产一区二区三区四区第35| 两个人看的免费小视频| 欧美日韩亚洲综合一区二区三区_| 国产精品免费视频内射| 午夜老司机福利片| 午夜福利乱码中文字幕| 啦啦啦视频在线资源免费观看| 男女免费视频国产| 飞空精品影院首页| av网站在线播放免费| 亚洲久久久国产精品| 国产成+人综合+亚洲专区| 99riav亚洲国产免费| 久久精品国产a三级三级三级| 国产亚洲午夜精品一区二区久久| 岛国毛片在线播放| www日本在线高清视频| 在线观看免费日韩欧美大片| 91成人精品电影| 久久影院123| 大陆偷拍与自拍| 女同久久另类99精品国产91| 亚洲国产av新网站| 欧美午夜高清在线| 在线观看舔阴道视频| 丝袜人妻中文字幕| 涩涩av久久男人的天堂| 国产一区二区 视频在线| 成在线人永久免费视频| 亚洲熟女毛片儿| 桃花免费在线播放| 黄片小视频在线播放| 美国免费a级毛片| 一边摸一边抽搐一进一小说 | 91av网站免费观看| 热re99久久精品国产66热6| 日韩大片免费观看网站| 男女床上黄色一级片免费看| 老司机深夜福利视频在线观看| 成年人黄色毛片网站| 老熟女久久久| av一本久久久久| 国产精品自产拍在线观看55亚洲 | 一本一本久久a久久精品综合妖精| 香蕉久久夜色| 一区二区三区精品91| 亚洲成人免费av在线播放| 久久久久久久久免费视频了| 一夜夜www| 亚洲国产欧美日韩在线播放| 美女扒开内裤让男人捅视频| 亚洲av日韩精品久久久久久密| 欧美精品亚洲一区二区| 亚洲国产av影院在线观看| 国产不卡av网站在线观看| 欧美乱码精品一区二区三区| 成人手机av| 国产又爽黄色视频| 欧美变态另类bdsm刘玥| 午夜日韩欧美国产| 超色免费av| 无遮挡黄片免费观看| 91老司机精品| 大型av网站在线播放| 男人舔女人的私密视频| 亚洲精品国产色婷婷电影| 99精国产麻豆久久婷婷| 精品国产超薄肉色丝袜足j| 极品人妻少妇av视频| 国产一卡二卡三卡精品| 久久精品国产综合久久久| 久久性视频一级片| 精品一区二区三卡| 成人精品一区二区免费| 69av精品久久久久久 | 国产男女内射视频| av有码第一页| 欧美成人午夜精品| 中文字幕色久视频| 亚洲精品国产精品久久久不卡| 黄色成人免费大全| 十分钟在线观看高清视频www| 99国产精品免费福利视频| 国产精品久久久久久精品电影小说| 三级毛片av免费| 欧美激情 高清一区二区三区| 不卡av一区二区三区| 操出白浆在线播放| 国产一卡二卡三卡精品| 两性午夜刺激爽爽歪歪视频在线观看 | 黄片小视频在线播放| 99国产精品一区二区三区| 99热国产这里只有精品6| 欧美午夜高清在线| 不卡av一区二区三区| 女警被强在线播放| a级片在线免费高清观看视频| 欧美性长视频在线观看| 99久久国产精品久久久| 亚洲,欧美精品.| 色婷婷久久久亚洲欧美| 国产精品秋霞免费鲁丝片| 深夜精品福利| 欧美日韩精品网址| 久久久久久亚洲精品国产蜜桃av| 亚洲熟妇熟女久久| 精品少妇一区二区三区视频日本电影| 久久精品国产99精品国产亚洲性色 | 国产欧美日韩一区二区精品| 精品人妻在线不人妻| 国产在视频线精品| 麻豆国产av国片精品| 在线永久观看黄色视频| 无限看片的www在线观看| 日本一区二区免费在线视频|