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

    基于SWAM模型的二元水循環(huán)與區(qū)域水資源配置匹配度研究

    2023-08-28 11:43:41楊木易
    水利水電快報(bào) 2023年7期
    關(guān)鍵詞:松花江流域

    楊木易

    摘要:為了實(shí)現(xiàn)水資源的高效利用和可持續(xù)發(fā)展,通過(guò)SWAM模型,獲得合理準(zhǔn)確的區(qū)域水資源配置和二元水循環(huán)分布二者匹配度綜合評(píng)價(jià)結(jié)果。在此基礎(chǔ)上,對(duì)松花江流域的水資源合理配置與水循環(huán)分布的合理性做出評(píng)判。運(yùn)用SWAM模型完成研究區(qū)域的水循環(huán)模擬,并且在各類約束下選擇最優(yōu)的多目標(biāo)平衡穩(wěn)定的水資源配置。結(jié)果表明:使用區(qū)域二元水資源循環(huán)系統(tǒng)匹配度最大的區(qū)域水資源分配方法,在P=50%時(shí),其所對(duì)應(yīng)的天然流域水資源總需求量為145.5億m3,總?cè)彼繛?.9億m3,總?cè)彼蕿?.3%;在P=90%時(shí),其所對(duì)應(yīng)的天然流域水資源總需求量為165.7億m3,總?cè)彼繛?6.8億m3,總?cè)彼蕿?1.3%。SWAM模型所模擬的二元水循環(huán)與水資源配置方案的匹配度可以全方位揭示經(jīng)濟(jì)用水過(guò)程和天然水文循環(huán)過(guò)程間相互作用,在保證二元水循環(huán)系統(tǒng)正常運(yùn)作的基礎(chǔ)上,為區(qū)域水資源配置提供參考。

    關(guān)鍵詞:水資源配置; SWAM模型; 二元水循環(huán); 松花江流域

    中圖法分類號(hào):TV213.4

    文獻(xiàn)標(biāo)志碼:A

    DOI:10.15974/j.cnki.slsdkb.2023.07.003

    文章編號(hào):1006-0081(2023)07-0021-06

    0 引 言

    流域水資源空間分配不均、水利基礎(chǔ)設(shè)施不足、工程管護(hù)管理不當(dāng),以及隨著社會(huì)經(jīng)濟(jì)高速發(fā)展而產(chǎn)生的需水量迅速增長(zhǎng)等因素,導(dǎo)致全域及局部地區(qū)產(chǎn)生了嚴(yán)重缺水和水生態(tài)污染等各種水資源問(wèn)題[1]。需運(yùn)用相應(yīng)的方法和策略,以保證人類未來(lái)經(jīng)濟(jì)生活的總需水量與水資源的供應(yīng)在時(shí)間和空間上的分配平衡,改變生態(tài)環(huán)境以適應(yīng)生產(chǎn)、生活、生態(tài)環(huán)境等3個(gè)領(lǐng)域的水資源需求,實(shí)現(xiàn)協(xié)調(diào)平衡的水資源配置。水資源配置是一個(gè)多目標(biāo)、階段區(qū)別明顯、主體多樣化的復(fù)雜問(wèn)題。水資源配置體系包括自然水資源循環(huán)系統(tǒng)、經(jīng)濟(jì)發(fā)展、人口、產(chǎn)業(yè)結(jié)構(gòu)和分布、生產(chǎn)形式、科學(xué)技術(shù)、自然資源狀況等要素,是復(fù)雜的大體系,需要運(yùn)用模型進(jìn)行研究[2]。水資源配置建模通??煞殖?種,分別是仿真模型、混合模型和優(yōu)化模型[3]。仿真模型通常按照節(jié)點(diǎn)的來(lái)水、供應(yīng)、調(diào)蓄、回歸,按照由上行至下行的次序完成逐節(jié)點(diǎn)水平衡運(yùn)算[4]。雖然此種仿真方法較為簡(jiǎn)單,但由于每次都只是仿真一種水資源配置方案,而同時(shí)模擬多個(gè)水資源配置方案時(shí)需要計(jì)算機(jī)進(jìn)行相關(guān)操作,還需要耗費(fèi)人力進(jìn)行調(diào)整、對(duì)比和選擇各方案,往往無(wú)法迅速得到最適宜的方案?;旌夏P鸵苑抡婺P偷挠?jì)算為前提,自動(dòng)進(jìn)行人工對(duì)比選擇方案過(guò)程,不足之處是計(jì)算繁瑣且工程量較大。而優(yōu)化模型則按照所選擇的對(duì)象數(shù)量,分成了單總體目標(biāo)綜合優(yōu)化建模[5]和多目標(biāo)綜合優(yōu)化建模[6]。兩種設(shè)計(jì)模式都使用優(yōu)化后的設(shè)計(jì)方法得到最優(yōu)化配置方法,優(yōu)點(diǎn)是能夠自行得到優(yōu)化方法,不足之處是運(yùn)算周期很長(zhǎng),運(yùn)算效率低,并且在解決某些效益問(wèn)題時(shí)程序的計(jì)算復(fù)雜程度會(huì)大大提高,但由于現(xiàn)在計(jì)算機(jī)技術(shù)更新速度快,上述的缺點(diǎn)隨著算法的優(yōu)化能夠得到大幅改善。

    水循環(huán)是水資源產(chǎn)生和發(fā)展的前提,是水資源合理配置的基礎(chǔ)。有研究者對(duì)基于水循環(huán)網(wǎng)絡(luò)的水資源合理配置模式進(jìn)行探索,如廣義自然資源配置模式,包括水資源合理配置、水循環(huán)仿真和水資源環(huán)境模擬共3個(gè)模式,分布式水資源合理配置仿真DTVGM-WEAR[7],概念性半分配式水資源整合仿真與調(diào)度模式WAS[8];基于傳統(tǒng)農(nóng)業(yè)水文模式的灌溉水優(yōu)化配置模式[9];采用改進(jìn)的SWAT模式,并嵌合傳統(tǒng)水資源分配模型的分布式水資源調(diào)度模式。上述研究成果中主要使用單目標(biāo)配置模型,僅少數(shù)研究成果中使用了多目標(biāo)配置模型,但由于在多目標(biāo)配置模型求解過(guò)程中常常采用人為的給定加權(quán),將其結(jié)果轉(zhuǎn)換為單目標(biāo)方程來(lái)解決,其結(jié)果易受到主觀影響而產(chǎn)生偏差[10]。同時(shí)如果僅采用一個(gè)指數(shù)對(duì)多目標(biāo)模型實(shí)現(xiàn)單目標(biāo)轉(zhuǎn)換,既無(wú)法反映眾多管理者的選擇需求,也無(wú)法衡量選擇方法對(duì)多種維度的綜合影響。因此,有些研究者把各類優(yōu)選方法都融入了多目標(biāo)優(yōu)選模式的計(jì)算中,如推出了一種通過(guò)目標(biāo)排序計(jì)算適配點(diǎn)的多目標(biāo)遺傳算法,并經(jīng)過(guò)多次迭代實(shí)現(xiàn)了水資源的多目標(biāo)優(yōu)選配置。本文通過(guò)SWAM模型,獲得合理準(zhǔn)確的區(qū)域水資源配置和二元水循環(huán)分布二者匹配度的綜合評(píng)價(jià)結(jié)果,在此基礎(chǔ)上對(duì)松花江流域的水資源合理配置與水循環(huán)分布二者的合理性做出評(píng)判,為當(dāng)?shù)厮Y源優(yōu)化調(diào)度提供參考依據(jù)。

    1 二元水循環(huán)與水資源配置聯(lián)系

    水資源配置和水循環(huán)相互依賴,水資源配置方案的制定要遵循水循環(huán)的原則。水的自然水文循環(huán)和社會(huì)循環(huán)共同構(gòu)成了整個(gè)水循環(huán)系統(tǒng)。將水資源體系的所有子系統(tǒng)緊密結(jié)合在一起,二元社會(huì)自然水循環(huán)影響了所有子系統(tǒng)相互之間的關(guān)系。

    自然水文循環(huán)是二元水循環(huán)的主體,在其中居于主導(dǎo)地位,它是指在太陽(yáng)輻射和地球引力的作用下,地球上的水不斷進(jìn)行相態(tài)轉(zhuǎn)化的循環(huán)過(guò)程。正因?yàn)樽匀凰难h(huán)的自發(fā)進(jìn)行,使自然界中的總體水量保持相對(duì)平衡,再加上循環(huán)中各個(gè)環(huán)節(jié)的相互作用,使得水體不斷的變化和流動(dòng),這為人類的生存和發(fā)展提供了最基本的物質(zhì)來(lái)源。社會(huì)水文循環(huán)是指社會(huì)生活中人的衣食住行對(duì)水的需要,通過(guò)取、用、耗、排水等方式,在自然水文循環(huán)上疊加的過(guò)程。取水循環(huán)是指利用人工干預(yù)的方法把自然界水文循環(huán)系統(tǒng)中的水體帶入社會(huì)水文循環(huán)系統(tǒng)的開(kāi)端階段,其目的就是讓天然水體創(chuàng)造出經(jīng)濟(jì)價(jià)值。用水環(huán)節(jié)和耗水環(huán)節(jié)正是實(shí)現(xiàn)這種目的最有效途徑,也是整個(gè)社會(huì)水文循環(huán)的核心。

    “二元水循環(huán)”的“二元”代表水循環(huán)在自然、社會(huì)兩個(gè)不同體系同時(shí)發(fā)生,隨著社會(huì)發(fā)展,兩者之間的聯(lián)系愈發(fā)緊密,整個(gè)循環(huán)環(huán)環(huán)相扣。由于經(jīng)濟(jì)發(fā)展,人口增多等社會(huì)因素,農(nóng)業(yè)、工業(yè)、日常用水等對(duì)水資源的需求不斷提高,這意味著愈來(lái)愈多以“自然”為單一水循環(huán)系統(tǒng)的水資源過(guò)渡為社會(huì)和自然綜合的水循環(huán)系統(tǒng),而在短期內(nèi)水資源的循環(huán)系統(tǒng)發(fā)生變化,隨之而來(lái)的問(wèn)題也愈發(fā)嚴(yán)重。由于水資源的過(guò)量開(kāi)采、水資源分配不當(dāng)?shù)纫幌盗幸蛩貙?dǎo)致自然環(huán)境破壞,地方缺水情況加重。因此,合理配置水資源,保證自然水資源循環(huán)所依附的生態(tài)環(huán)境不因水資源缺乏而被破壞,合理分配各用水單位的水資源供應(yīng)量,保證社會(huì)經(jīng)濟(jì)能夠正常發(fā)展,進(jìn)而能夠保證水資源在自然、社會(huì)雙體系下正常穩(wěn)定地循環(huán)。在社會(huì)水文循環(huán)中,農(nóng)業(yè)生產(chǎn)和工業(yè)生產(chǎn)是最主要的需水對(duì)象,因此,清楚水資源循環(huán)和區(qū)域水資源配置之間的匹配度,在此基礎(chǔ)上選擇匹配度優(yōu)先級(jí)最高的水資源配置方案,既能保證農(nóng)業(yè)、工業(yè)生產(chǎn)正常用水,也能維持自然、社會(huì)水循環(huán)系統(tǒng)的穩(wěn)定。

    2 SWAM模型

    水循環(huán)和水資源配置模型以二元水循環(huán)為基礎(chǔ),對(duì)水資源的配置進(jìn)行多目標(biāo)均衡最優(yōu)化處理,在優(yōu)化過(guò)程中考慮各類限制條件和系統(tǒng)規(guī)則,因此最終的水資源配置結(jié)果不僅能夠最優(yōu)化配置水資源,同時(shí)也能滿足現(xiàn)實(shí)的客觀實(shí)際需求。該模型由多個(gè)模型模塊共同組成,其中有需水、多水源分質(zhì)供水、水循環(huán)仿真、多目標(biāo)分配等模型。水循環(huán)仿真模型能夠計(jì)算出合理的水資源供應(yīng)量,這是多水源分質(zhì)供水模型所需要的重要參數(shù)。當(dāng)多目標(biāo)配置模型計(jì)算了水資源分配方法后,篩選出推薦方法,再展開(kāi)水循環(huán)計(jì)算,直到結(jié)果合理。

    2.1 水循環(huán)模塊

    自然水循環(huán)部分是以網(wǎng)格中新安江模型為底層原理加以改進(jìn)。模型的蒸擴(kuò)散采取三層蒸發(fā)模式,充分討論了土質(zhì)含水百分率沿垂直方向分布的特異性。產(chǎn)流部分的模擬計(jì)算以上述模型的蓄滿產(chǎn)流為主要方法。對(duì)于匯流部分,流域單元研究區(qū)域的土壤水流量和地下水流量選取單位線法計(jì)算單元水資源出流量,依次計(jì)算各個(gè)研究單元的匯流順序,分別堆疊,直到流域產(chǎn)出。

    2.2 多水源分質(zhì)供水模塊

    可供用水總量和各水源的取水控制總量的約束條件主要來(lái)自于流域的水資源管理用水的控制紅線和水資源循環(huán)系統(tǒng)的有效控制。實(shí)際取水量主要取決于所計(jì)算單位的需水總量,以及飲用水來(lái)源(如河道、水庫(kù)、地下含水層等)的可供給量。

    地表水資源(含外調(diào)水)供應(yīng):地表水(含外調(diào)水)的可供應(yīng)量取決于水資源的利用水量、取水工程建設(shè)(如引水渠道、給排水管線、機(jī)井等)的水資源供應(yīng)效率和水源地(包含境內(nèi)地表水與外調(diào)水)的飲用水管理紅線等制約要求。

    水庫(kù)水資源供應(yīng):水庫(kù)庫(kù)容的計(jì)算根據(jù)上一時(shí)段末庫(kù)容總量和本時(shí)期水庫(kù)上游來(lái)水、庫(kù)面蒸發(fā)、水庫(kù)滲漏量和水庫(kù)供水情況共同確定。

    再生水資源供應(yīng):由于水資源匱乏情況日益突出,因此將非常規(guī)水資源引入供給系統(tǒng)十分必要。目前,中國(guó)對(duì)非常規(guī)水資源的使用尚處在探索時(shí)期,非常規(guī)水資源常用作供應(yīng)給生態(tài)、農(nóng)業(yè)、工業(yè)等,為實(shí)行水資源分質(zhì)供應(yīng)奠定了基礎(chǔ)。

    2.3 需水模塊

    2.3.1 社會(huì)經(jīng)濟(jì)發(fā)展數(shù)據(jù)

    本模型還需提供現(xiàn)狀年的社區(qū)數(shù)據(jù)和發(fā)展規(guī)劃水平年增長(zhǎng)資料。社區(qū)數(shù)據(jù)包含:常駐居民(鄉(xiāng)村、城鎮(zhèn)居民)、大家畜、小家畜,工業(yè)、建設(shè)、服務(wù)業(yè)生產(chǎn)總值、糧食作物(麥子、包谷、水稻、經(jīng)濟(jì)作物)栽培面積、設(shè)計(jì)灌溉面積、森林草地面積、畜牧業(yè)規(guī)模和增長(zhǎng)率等。關(guān)于社區(qū)經(jīng)濟(jì)規(guī)模層次的發(fā)展資料包含:人均增長(zhǎng)率、GDP增長(zhǎng)速度、農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)的變化量、糧食作物播種面積應(yīng)變?cè)鲩L(zhǎng)率等。

    2.3.2 水資源需求量

    水資源需求量計(jì)算包括畜牧、工業(yè)、服務(wù)業(yè)與建筑業(yè)、生活、農(nóng)業(yè)種植業(yè)。農(nóng)業(yè)的水資源需求量為林業(yè)、草業(yè)、漁業(yè)、灌溉水資源需求量的總和。

    2.4 多目標(biāo)配置模塊

    本模型共設(shè)定了3個(gè)主要目標(biāo):經(jīng)濟(jì)總體目標(biāo)、社會(huì)總體目標(biāo)與水量總體目標(biāo)。經(jīng)濟(jì)總體目標(biāo)以區(qū)域效益最高來(lái)表現(xiàn),經(jīng)濟(jì)社會(huì)總體目標(biāo)以因水資源不足而導(dǎo)致的糧食減產(chǎn)量最少來(lái)表現(xiàn),水量總體目標(biāo)以水資源缺乏量最少來(lái)表現(xiàn)。

    2.5 信息聚合

    信息聚合也可理解為數(shù)據(jù)整合,是指將水資源供應(yīng)和用水單位的水資源消耗的數(shù)據(jù)連接運(yùn)算。水源供水核算是水資源配置模型和水循環(huán)模塊之間的關(guān)鍵階段。對(duì)各個(gè)計(jì)算模塊內(nèi)的各種水資源進(jìn)行統(tǒng)計(jì),再同步對(duì)各個(gè)用水單位的水資源消耗量調(diào)整控制。將所有用水單位的水資源供應(yīng)量求和,求和所得的水資源供應(yīng)總量數(shù)值是水循環(huán)模型對(duì)水循環(huán)進(jìn)行仿真的重要參數(shù)。

    3 研究區(qū)概況

    松花江流域(圖1)地處中國(guó)大陸的東北區(qū)域,位于北緯41°39′~51°40′,東經(jīng)119°49′~132°29′。松花江流域東北以小興安嶺山地為界,西達(dá)大興安嶺,東側(cè)地區(qū)以完達(dá)山余脈、老爺嶺、張大才嶺、長(zhǎng)白山等為界,西北部的山麓地區(qū)則為松花江與遼河2個(gè)主要流域面積的分水嶺。

    松花江是中國(guó)七大河流之一,總長(zhǎng)1 931 km,東西橫跨930 km,南北寬約1 080 km,橫跨中國(guó)內(nèi)蒙古、吉林省、黑龍江三省,面積約56.79萬(wàn)km2,占黑龍江省總面積(約185.5萬(wàn)km2)的30.61%,年均徑流量773億m3。

    本次研究選取的城市分布在松花江流域干、支流水源區(qū)域,這些城市分別在人口、耕地面積、工業(yè)程度上具有典型代表性,如吉林市區(qū)、松原市區(qū)、扶余縣、哈爾濱、佳木斯等地區(qū),對(duì)上述地區(qū)進(jìn)行供水預(yù)測(cè)能夠在一定程度上體現(xiàn)松花江流域在規(guī)劃水平年的供水情況。

    3.1 水資源總量

    1960~2021年松花江流域多年平均值水資源量為941.7億m3,同時(shí)地表水資源量平均值為792.5億m3,其他自然資源量為348.9億m3,地下水資源和地表水資源的不重復(fù)總量均為169.3億m3。

    3.2 水資源可利用量

    1960~2021年松花江流域水資源需求量為511.8億m3,實(shí)際可利用率為55.3%;地表水的可開(kāi)發(fā)利用總量為378.2億m3,可開(kāi)發(fā)利用度為48.1%;平原區(qū)淺層地下水可開(kāi)采量181.0億m3,占總補(bǔ)給量的80.9%。

    4 供水預(yù)測(cè)

    本文基礎(chǔ)年份定為2021年,而規(guī)劃年份則設(shè)定在2025年??紤]來(lái)水頻率分別為P=50%,P=90%時(shí)的情況。

    根據(jù)1960~2021年松花江流域60 a降水量排頻結(jié)果,來(lái)水率為50%的對(duì)應(yīng)的年份為2014年、來(lái)水率為90%的對(duì)應(yīng)年份為2006年。將2014年和2006年的降水量數(shù)據(jù)導(dǎo)入模型,根據(jù)水循環(huán)模擬數(shù)據(jù),依次計(jì)算該流域各個(gè)用水單位的可供水量,見(jiàn)表1。從表1中得出結(jié)果:當(dāng)來(lái)水頻率為50%時(shí),總水資源可供應(yīng)量為146.9億m3,包括地表水資源總供應(yīng)量95.0億m3,地下水水資源總供應(yīng)量51.4億m3,其他水資源總供應(yīng)量為0.5億m3;當(dāng)來(lái)水頻率為90%時(shí),總水資源可供應(yīng)量為144.3億m3,包括地表水資源總供應(yīng)量90.7億m3,地下水水資源總供應(yīng)量53.3億m3,其他水資源總供應(yīng)量為0.3億m3。

    5 結(jié)果與討論

    5.1 水循環(huán)轉(zhuǎn)化關(guān)系

    對(duì)松花江流域的水循環(huán)系統(tǒng)進(jìn)行仿真,得出在來(lái)水頻率為50%時(shí),全流域的總降雨量為1 436.0億m3,其中山區(qū)降雨量為812.0億m3,平原地帶降雨量為535.0億m3;全流域的總蒸散發(fā)總量為855.0億m3,其中山區(qū)的總蒸發(fā)能力為544.0億m3,平原地帶的總蒸發(fā)能力為364.0億m3;全流域的下滲總量為152.3億m3,其中山區(qū)的總下滲量為54.9億m3,平原地帶的總下滲量為88.5億m3;全流域面積年末總土地水蓄變量為2.4億m3,進(jìn)境供水為146.0億m3,過(guò)境供水為318.0億m3;全流域面積經(jīng)濟(jì)社會(huì)總體取水量為142.6億m3,基本匯集于平原區(qū),其中地表水取水量為89.1億m3,淺層地下水總體取水量為43.6億m3。在來(lái)水頻率為90%時(shí),全流域面積的總降雨量為1 092.0億m3,其中山地降雨量為659.7億m3,平原地帶降雨量為443.5億m3;全流域面積的總蒸散發(fā)量為954.5億m3,其中山地的總蒸發(fā)能力為586.7億m3,平原地帶的總蒸發(fā)能力為388.8億m3;山地的下滲量為16.3億m3,平川區(qū)的下滲量為25.6億m3??梢?jiàn),松花江流域主要的水資源供給來(lái)源為自然大氣降水和地表徑流,地下水的水資源供應(yīng)源則為自然雨水、河流湖庫(kù)的滲入等;流域面積內(nèi)最主要的水資源消耗項(xiàng)為蒸散發(fā)和經(jīng)濟(jì)取水,所占比依次為87.6%,14.6%。在旱災(zāi)時(shí)期,由于降水少,導(dǎo)致土地水蓄變量下降,流域徑流量大幅下降。旱災(zāi)時(shí)期因社會(huì)取水增多,造成地下水開(kāi)采量上升。

    5.2 水資源配置結(jié)果

    復(fù)雜水資源分配體系的影響因子很多,包括天然水循環(huán)(自然來(lái)水的隨機(jī)性與不確定性)、社會(huì)經(jīng)濟(jì)文化因素、城市供水能力、低碳技術(shù)、城市化、生態(tài)環(huán)境因素等,而大量的半結(jié)構(gòu)性、無(wú)結(jié)構(gòu)性因素也要求水資源配置方案的決策者做出重新評(píng)估與選擇,而決策者的偏好也將改變水資源配置體系。

    本部分將介紹SWAM模型的應(yīng)用結(jié)果,包括來(lái)水頻率P為50%和90%時(shí)的水資源配置方案結(jié)果。在SWAM模型的基礎(chǔ)上,本研究采用了需水量、供水量、缺水量和缺水率等指標(biāo),對(duì)水資源配置方案進(jìn)行評(píng)價(jià)。需水量是指區(qū)域內(nèi)各個(gè)用水部門(mén)的水需求量,供水量是指水資源供給的總量,缺水量則是指需水量與供水量之差,缺水率則是指缺水量與需水量的比值,具體數(shù)值見(jiàn)表2和表3。

    規(guī)劃水平年至2025年,松花江流域的水循環(huán)變化情況如下:當(dāng)來(lái)水頻率為50%時(shí),該流域面積的自然環(huán)境降水補(bǔ)給量為1 436.0億m3,蒸散發(fā)量為855.0億m3,下滲水量為152.3億m3,土壤蓄水能力變化為2.4億m3,入境供水為146.0億m3,出境供水為318.0億m3;在沖積平原區(qū),經(jīng)濟(jì)社會(huì)總?cè)」┧疄?42.0億m3,同時(shí),地面環(huán)境取供水為89.1億m3,淺層地下水平均取供水為53.6億m3。在來(lái)水頻率為90%時(shí),該流域面積的環(huán)境降水補(bǔ)給量為1 092.0億m3,蒸散發(fā)量為954.5億m3,流域年末土壤蓄水能力變化為-5.3億m3;平原地區(qū)水資源總?cè)∮昧繛?49.9億m3,地下水取使用量為84.0億m3,淺層地下水水資源總?cè)∮昧繛?7.2億m3。

    由于水資源配置結(jié)果由多目標(biāo)的結(jié)果共同體現(xiàn),針對(duì)單一目標(biāo)進(jìn)行優(yōu)化只能得出在該目標(biāo)上的最佳配置方案,因此為了體現(xiàn)水資源配置方案的綜合效益,對(duì)各現(xiàn)有方案進(jìn)行多目標(biāo)的配置決策評(píng)價(jià),給每個(gè)水資源配置方案最終優(yōu)化目標(biāo)設(shè)定權(quán)重,而pareto解集計(jì)算所得最優(yōu)解能夠避免各個(gè)目標(biāo)權(quán)重在設(shè)定時(shí)存在主觀性,通過(guò)最終的多目標(biāo)權(quán)重的博弈選擇出更加適配地方特性的水資源配置方案[11-12]。通過(guò)在經(jīng)濟(jì)效益目標(biāo)、社會(huì)目標(biāo)、用水量目標(biāo)之間的博弈,如表4所示,共得出了14組非劣解,以經(jīng)濟(jì)效益最好、糧食減產(chǎn)量最少為基礎(chǔ),篩選出來(lái)水頻率為50%,90%時(shí)與二元水循環(huán)系統(tǒng)匹配度最高的水資源配置方案為最優(yōu)方法,其相應(yīng)的流域內(nèi)水資源總供應(yīng)量為145.5億m3,總?cè)彼繛?.9億m3,總?cè)彼蕿?.3%。在該方法下生活供水量(不含生態(tài)需水量)、工業(yè)生產(chǎn)用水量和農(nóng)業(yè)需水量分別為10.8億、17.5億、99.9億m3。在各種供應(yīng)水源地中,地表水供水量最大,為99.5億m3,其次為地下水,為53.6億m3,其余水源中再生水最少,為0.6億m3,在該方法下的經(jīng)濟(jì)損失為2.6億元,糧食減產(chǎn)量為62.0萬(wàn)t。當(dāng)來(lái)水頻率為90%,松花江流域的水資源總需求量為165.7億m3,水資源總供應(yīng)量為149.0億m3,總?cè)彼繛?6.8億m3,總用水虧缺率為11.3%。

    6 結(jié) 語(yǔ)

    本文建立了水循環(huán)和水資源配置模型,系統(tǒng)地介紹了建??蚣?、建模流程、求解方式。模塊中主要的功能模塊如下:水循環(huán)模塊、多源分質(zhì)供水功能、需水功能、多目標(biāo)配置功能。同時(shí),本文以松花江流域?yàn)槔?,基于SWAM模型,以二元水循環(huán)系統(tǒng)和水資源配置方案的匹配度為依據(jù),對(duì)研究區(qū)域的水資源配置方案進(jìn)行優(yōu)化。通過(guò)SWAM模型而進(jìn)行的二元水循環(huán)與水資源使用的匹配性分析可以充分反映經(jīng)濟(jì)用水活動(dòng)和天然水文循環(huán)活動(dòng)間的相互作用,了解水資源循環(huán)和區(qū)域水資源配置之間的匹配度,在此基礎(chǔ)上選擇匹配度優(yōu)先級(jí)最高的水資源配置方案。這樣既能保證農(nóng)業(yè)、工業(yè)生產(chǎn)正常用水,也能維持自然、社會(huì)水循環(huán)系統(tǒng)的穩(wěn)定。通過(guò)這種方式,在保證二元水循環(huán)系統(tǒng)正常運(yùn)作的基礎(chǔ)上,為區(qū)域水資源配置提供參考。

    參考文獻(xiàn):

    [1] 趙晶,王濤,畢彥杰,等.水循環(huán)與水資源合理配置模型SWAM(Ⅰ):模型原理與構(gòu)建[J/OL].華北水利水電大學(xué)學(xué)報(bào)(自然科學(xué)版):1-10[2023-01-07].http:∥kns.cnki.net/kcms/detail/41.1432.TV.20221012.1528.003.html.

    [2] 劉建軍.南疆地區(qū)水資源保障問(wèn)題和對(duì)策分析[J].水資源開(kāi)發(fā)與管理,2022,8(8):8-11,5.

    [3] 王煜,彭少明,周翔南,等.基于增量動(dòng)態(tài)均衡配置的黃河分水方案優(yōu)化方法研究[J].水資源保護(hù),2022,38(4):48-55.

    [4] 李致家,張珂,王棟,等.現(xiàn)代水文模擬與預(yù)報(bào)技術(shù)(第2版)[M].南京:河海大學(xué)出版社,2021.

    [5] 楊明智,許繼軍,桑連海,等.基于水循環(huán)的分布式水資源調(diào)配模型開(kāi)發(fā)與應(yīng)用[J].水利學(xué)報(bào),2022,53(4):456-470.

    [6] 劉移勝,閆少鋒.基于拓?fù)潢P(guān)系的水資源配置模型研究及應(yīng)用[J].中國(guó)農(nóng)村水利水電,2022(11):26-31.

    [7] 曹永強(qiáng),李玲慧.基于水資源綜合模擬與調(diào)配一體化模型的北京市水資源模擬分析[J].水利水電科技進(jìn)展,2021,41(6):25-31.

    [8] 張真,張光明,張盼月.基于多目標(biāo)不確定性機(jī)會(huì)約束規(guī)劃的礦區(qū)水資源配置研究[J].水資源與水工程學(xué)報(bào),2022,33(5):53-61.

    [9] 樊紅梅,劉曉民,劉廷璽,等.基于空間均衡的水資源合理配置研究[J].水資源與水工程學(xué)報(bào),2022,33(2):61-67.

    [10] 尋昊.三亞市水資源承載力分析及水資源配置優(yōu)化研究[J].中國(guó)給水排水,2022,38(13):59-63.

    [11] 卞錦宇,夏玉林,毛學(xué)謙,等.太湖流域水資源承載力與經(jīng)濟(jì)社會(huì)發(fā)展適應(yīng)性評(píng)價(jià)[J].水利水電快報(bào),2022,43(4):31-37.

    [12] 童坤,孫偉,王小軍.基于水環(huán)境承載能力與壓力的產(chǎn)業(yè)布局引導(dǎo)研究——以長(zhǎng)江干流39個(gè)城市為例[J].人民長(zhǎng)江,2021,52(7):54-59.

    (編輯:江 文)

    猜你喜歡
    松花江流域
    松花江流域吉林段土壤侵蝕變化分析
    清代松花江流域洪澇災(zāi)害及政府救助(1736—1911年)
    松花江百里生態(tài)長(zhǎng)廊生境網(wǎng)絡(luò)構(gòu)建研究
    自然災(zāi)害、農(nóng)戶傳統(tǒng)風(fēng)險(xiǎn)管理手段與天氣指數(shù)保險(xiǎn)需求
    吉林省松花江流域產(chǎn)業(yè)生態(tài)系統(tǒng)適應(yīng)性分析
    基于一維水動(dòng)力模型的洪水頂托影響分析
    農(nóng)業(yè)面源氨氮污染減排體系的構(gòu)建及減排控制對(duì)策
    尼爾基 豐滿水庫(kù)攔蓄作用分析
    松花江流域沿江重點(diǎn)監(jiān)控企業(yè)水環(huán)境潛在污染風(fēng)險(xiǎn)分析
    松花江流域不同空間尺度典型流域泥沙輸移比及其影響因素
    午夜免费观看网址| 视频区图区小说| 99久久国产精品久久久| 男男h啪啪无遮挡| 黑人巨大精品欧美一区二区mp4| 精品久久久久久,| 日韩熟女老妇一区二区性免费视频| 欧美日韩成人在线一区二区| 欧美大码av| 亚洲欧美一区二区三区久久| 精品人妻1区二区| 99国产精品99久久久久| 麻豆乱淫一区二区| 久久精品aⅴ一区二区三区四区| 俄罗斯特黄特色一大片| 性色av乱码一区二区三区2| aaaaa片日本免费| 美女福利国产在线| 亚洲少妇的诱惑av| 免费女性裸体啪啪无遮挡网站| 黄色片一级片一级黄色片| 精品少妇久久久久久888优播| 日本黄色视频三级网站网址 | 欧美成人免费av一区二区三区 | 18禁国产床啪视频网站| 国产成人精品无人区| 久久精品成人免费网站| 色综合欧美亚洲国产小说| 久久精品成人免费网站| 欧美乱码精品一区二区三区| 国产1区2区3区精品| 亚洲欧洲精品一区二区精品久久久| 黄色丝袜av网址大全| av国产精品久久久久影院| www日本在线高清视频| 不卡av一区二区三区| 自线自在国产av| 女人久久www免费人成看片| 亚洲熟女精品中文字幕| 在线观看免费高清a一片| 亚洲专区中文字幕在线| 久久人妻福利社区极品人妻图片| 久久久国产一区二区| 又黄又爽又免费观看的视频| 久久狼人影院| 国产欧美日韩一区二区三| 亚洲一码二码三码区别大吗| 夜夜夜夜夜久久久久| 色播在线永久视频| 啦啦啦免费观看视频1| 亚洲精品中文字幕一二三四区| 可以免费在线观看a视频的电影网站| 亚洲专区字幕在线| 80岁老熟妇乱子伦牲交| 大香蕉久久成人网| 99精品在免费线老司机午夜| av天堂在线播放| 成在线人永久免费视频| e午夜精品久久久久久久| 搡老熟女国产l中国老女人| 国产野战对白在线观看| 天堂中文最新版在线下载| 少妇裸体淫交视频免费看高清 | 日韩视频一区二区在线观看| 久久这里只有精品19| 在线观看免费高清a一片| 午夜精品国产一区二区电影| 国产av精品麻豆| 最新在线观看一区二区三区| 精品高清国产在线一区| 国产无遮挡羞羞视频在线观看| 国产欧美日韩一区二区三| 19禁男女啪啪无遮挡网站| 成人黄色视频免费在线看| av国产精品久久久久影院| 午夜精品国产一区二区电影| 丝袜美腿诱惑在线| 国内久久婷婷六月综合欲色啪| 在线观看66精品国产| 18禁国产床啪视频网站| 国产精品美女特级片免费视频播放器 | 日韩熟女老妇一区二区性免费视频| 国产精品成人在线| 亚洲欧美一区二区三区久久| 午夜福利乱码中文字幕| a在线观看视频网站| 久久久久久久久免费视频了| 成人av一区二区三区在线看| 90打野战视频偷拍视频| 久久久国产成人精品二区 | 久久国产乱子伦精品免费另类| 在线天堂中文资源库| 国产精品美女特级片免费视频播放器 | 操出白浆在线播放| 午夜日韩欧美国产| 日韩欧美一区视频在线观看| 亚洲色图 男人天堂 中文字幕| 亚洲色图综合在线观看| tube8黄色片| 交换朋友夫妻互换小说| av国产精品久久久久影院| 日韩制服丝袜自拍偷拍| 亚洲一码二码三码区别大吗| 搡老乐熟女国产| 999久久久精品免费观看国产| 精品人妻1区二区| 国产一区二区三区综合在线观看| 国产成人免费观看mmmm| 麻豆成人av在线观看| 日本精品一区二区三区蜜桃| 国产欧美日韩综合在线一区二区| 国产免费男女视频| 一本一本久久a久久精品综合妖精| 精品熟女少妇八av免费久了| a级毛片黄视频| 国产一卡二卡三卡精品| 性少妇av在线| 久久香蕉精品热| 99在线人妻在线中文字幕 | 日本五十路高清| 男女之事视频高清在线观看| 久久人妻av系列| 纯流量卡能插随身wifi吗| 大陆偷拍与自拍| 久久久久国产一级毛片高清牌| 亚洲专区国产一区二区| 91成年电影在线观看| 可以免费在线观看a视频的电影网站| 亚洲色图av天堂| 一本一本久久a久久精品综合妖精| 身体一侧抽搐| 亚洲午夜精品一区,二区,三区| 成年动漫av网址| 久久精品aⅴ一区二区三区四区| 日韩欧美一区二区三区在线观看 | 高潮久久久久久久久久久不卡| 12—13女人毛片做爰片一| 精品久久久久久久毛片微露脸| 欧美黄色淫秽网站| av免费在线观看网站| 天堂动漫精品| 一级片'在线观看视频| 久久国产精品影院| 成年人免费黄色播放视频| 国产亚洲欧美98| 精品福利观看| 日韩三级视频一区二区三区| 欧美中文综合在线视频| 天天躁夜夜躁狠狠躁躁| 老鸭窝网址在线观看| 精品少妇久久久久久888优播| 国产视频一区二区在线看| 一区二区三区激情视频| 欧美+亚洲+日韩+国产| 亚洲av成人不卡在线观看播放网| 亚洲精品自拍成人| 天天影视国产精品| 一进一出抽搐gif免费好疼 | 成年版毛片免费区| 免费看a级黄色片| 女性生殖器流出的白浆| 欧美色视频一区免费| 黑人欧美特级aaaaaa片| 黑人巨大精品欧美一区二区mp4| 国产精品香港三级国产av潘金莲| 中文字幕人妻丝袜制服| bbb黄色大片| 中文字幕av电影在线播放| 国产精品久久久久久精品古装| 国产欧美日韩一区二区三| 老鸭窝网址在线观看| 一级片'在线观看视频| 免费日韩欧美在线观看| 亚洲中文字幕日韩| 日本欧美视频一区| 亚洲av电影在线进入| 中文字幕人妻丝袜制服| 夜夜躁狠狠躁天天躁| 正在播放国产对白刺激| 国产成+人综合+亚洲专区| 国产亚洲av高清不卡| 国产真人三级小视频在线观看| 色94色欧美一区二区| 久久人妻福利社区极品人妻图片| 精品久久久久久电影网| 又黄又粗又硬又大视频| 欧美日韩精品网址| 国产成人欧美| 亚洲成av片中文字幕在线观看| 后天国语完整版免费观看| 精品一区二区三区av网在线观看| 亚洲国产中文字幕在线视频| 女人爽到高潮嗷嗷叫在线视频| 妹子高潮喷水视频| 另类亚洲欧美激情| 亚洲五月天丁香| 色综合欧美亚洲国产小说| 中文字幕精品免费在线观看视频| 天天躁狠狠躁夜夜躁狠狠躁| 嫁个100分男人电影在线观看| 久久久久国产一级毛片高清牌| 丁香六月欧美| 国产视频一区二区在线看| 亚洲精品国产色婷婷电影| 满18在线观看网站| 欧美中文综合在线视频| 999久久久精品免费观看国产| 在线永久观看黄色视频| 人人妻人人澡人人爽人人夜夜| 操出白浆在线播放| 国产单亲对白刺激| 黑人操中国人逼视频| 99久久人妻综合| 女人久久www免费人成看片| 亚洲人成电影免费在线| 最新美女视频免费是黄的| 国产一区二区三区视频了| 99在线人妻在线中文字幕 | 欧美乱色亚洲激情| 99国产极品粉嫩在线观看| 国产精品自产拍在线观看55亚洲 | 国产日韩欧美亚洲二区| 狠狠婷婷综合久久久久久88av| 人妻久久中文字幕网| 在线观看免费视频日本深夜| 18禁国产床啪视频网站| 后天国语完整版免费观看| 日日夜夜操网爽| 亚洲精品在线观看二区| 99久久99久久久精品蜜桃| 欧美日韩亚洲高清精品| 女人久久www免费人成看片| 欧美日韩一级在线毛片| xxx96com| 王馨瑶露胸无遮挡在线观看| 一区二区三区激情视频| a级片在线免费高清观看视频| a在线观看视频网站| 最新美女视频免费是黄的| 免费在线观看影片大全网站| 五月开心婷婷网| 男人的好看免费观看在线视频 | а√天堂www在线а√下载 | 建设人人有责人人尽责人人享有的| 亚洲美女黄片视频| 免费高清在线观看日韩| 黑人猛操日本美女一级片| 18禁观看日本| 日韩三级视频一区二区三区| 一级,二级,三级黄色视频| 大香蕉久久网| 国产精品美女特级片免费视频播放器 | 亚洲 欧美一区二区三区| 老司机靠b影院| 久久久久精品人妻al黑| 中文字幕人妻熟女乱码| 日日摸夜夜添夜夜添小说| 亚洲 国产 在线| 精品国产一区二区三区久久久樱花| 99久久综合精品五月天人人| 国产亚洲精品一区二区www | 69精品国产乱码久久久| 国产99久久九九免费精品| av天堂在线播放| 欧美日韩亚洲综合一区二区三区_| 女人被狂操c到高潮| 成年版毛片免费区| 欧美激情 高清一区二区三区| 中文字幕人妻熟女乱码| 欧美国产精品va在线观看不卡| tube8黄色片| 黄片大片在线免费观看| 窝窝影院91人妻| 久久人人97超碰香蕉20202| 一夜夜www| 高清在线国产一区| 身体一侧抽搐| 久久国产精品大桥未久av| 成人三级做爰电影| 精品国产亚洲在线| 国产精品自产拍在线观看55亚洲 | 国产亚洲精品第一综合不卡| 满18在线观看网站| svipshipincom国产片| 国产精品99久久99久久久不卡| 黄片小视频在线播放| 别揉我奶头~嗯~啊~动态视频| 在线国产一区二区在线| 国产不卡一卡二| 午夜精品在线福利| 韩国av一区二区三区四区| 国产精品香港三级国产av潘金莲| 香蕉国产在线看| 欧美激情久久久久久爽电影 | 天堂俺去俺来也www色官网| 天堂动漫精品| 老司机午夜十八禁免费视频| 国产精品一区二区精品视频观看| 91麻豆精品激情在线观看国产 | 在线十欧美十亚洲十日本专区| 国产免费现黄频在线看| 搡老乐熟女国产| 91成人精品电影| 国产黄色免费在线视频| 久久 成人 亚洲| 欧美午夜高清在线| 久久ye,这里只有精品| avwww免费| 免费在线观看黄色视频的| 99国产精品一区二区三区| 久久中文看片网| 18禁裸乳无遮挡动漫免费视频| 一进一出抽搐gif免费好疼 | 国产成人精品无人区| 亚洲欧洲精品一区二区精品久久久| 日韩制服丝袜自拍偷拍| 亚洲aⅴ乱码一区二区在线播放 | 精品免费久久久久久久清纯 | 他把我摸到了高潮在线观看| 在线观看免费日韩欧美大片| 99精品在免费线老司机午夜| 国产在线观看jvid| 久久精品亚洲av国产电影网| 亚洲av片天天在线观看| 久久久国产欧美日韩av| 亚洲精品美女久久久久99蜜臀| 人人妻人人澡人人看| www日本在线高清视频| 久久中文看片网| av福利片在线| 亚洲精品一卡2卡三卡4卡5卡| 男女床上黄色一级片免费看| 日本一区二区免费在线视频| 亚洲一区中文字幕在线| 深夜精品福利| 成人国语在线视频| 国产精品亚洲一级av第二区| 老司机影院毛片| √禁漫天堂资源中文www| 飞空精品影院首页| 男女下面插进去视频免费观看| 免费在线观看日本一区| 91成人精品电影| 欧美日韩福利视频一区二区| 麻豆乱淫一区二区| 欧美亚洲 丝袜 人妻 在线| 麻豆成人av在线观看| 久久久久久久精品吃奶| 国产精品99久久99久久久不卡| 狠狠狠狠99中文字幕| 亚洲自偷自拍图片 自拍| 国产精品自产拍在线观看55亚洲 | av网站免费在线观看视频| 日韩免费高清中文字幕av| 亚洲男人天堂网一区| 久久草成人影院| 亚洲av成人不卡在线观看播放网| 很黄的视频免费| 91字幕亚洲| 一边摸一边抽搐一进一小说 | 亚洲性夜色夜夜综合| 国产成人啪精品午夜网站| 欧美在线黄色| 熟女少妇亚洲综合色aaa.| 国产深夜福利视频在线观看| 欧美日韩成人在线一区二区| 中文字幕精品免费在线观看视频| 黄色视频不卡| 窝窝影院91人妻| 亚洲精品美女久久av网站| 久久青草综合色| 99国产精品免费福利视频| 人人澡人人妻人| 在线观看一区二区三区激情| 两人在一起打扑克的视频| 国产欧美日韩一区二区三区在线| 日本wwww免费看| 国产单亲对白刺激| 大型黄色视频在线免费观看| 美女 人体艺术 gogo| 999久久久精品免费观看国产| 亚洲人成电影观看| 香蕉久久夜色| 99精国产麻豆久久婷婷| 又黄又粗又硬又大视频| 亚洲av第一区精品v没综合| 亚洲中文日韩欧美视频| 亚洲av成人av| 久久国产亚洲av麻豆专区| 午夜福利免费观看在线| 欧美人与性动交α欧美软件| 久久久久久久精品吃奶| 国产精品1区2区在线观看. | 久久午夜综合久久蜜桃| 久久久国产欧美日韩av| 99精品久久久久人妻精品| 女同久久另类99精品国产91| 国产成人精品久久二区二区91| 韩国av一区二区三区四区| av有码第一页| 搡老熟女国产l中国老女人| 亚洲av美国av| 人妻久久中文字幕网| av天堂在线播放| 欧美丝袜亚洲另类 | 首页视频小说图片口味搜索| 新久久久久国产一级毛片| 国产亚洲欧美98| 免费少妇av软件| 久9热在线精品视频| 欧美成人午夜精品| 亚洲五月天丁香| 国产精品98久久久久久宅男小说| 久久性视频一级片| 中文字幕制服av| 很黄的视频免费| 亚洲 国产 在线| 色精品久久人妻99蜜桃| 两性午夜刺激爽爽歪歪视频在线观看 | 1024香蕉在线观看| 国产高清videossex| 成人国产一区最新在线观看| 男人的好看免费观看在线视频 | 欧美日韩亚洲国产一区二区在线观看 | 日本a在线网址| 看免费av毛片| a级片在线免费高清观看视频| 久久久久视频综合| 久久久久久人人人人人| 日韩人妻精品一区2区三区| 一区二区三区国产精品乱码| 日韩欧美三级三区| 大码成人一级视频| 12—13女人毛片做爰片一| 交换朋友夫妻互换小说| 色播在线永久视频| 深夜精品福利| 男女高潮啪啪啪动态图| 亚洲少妇的诱惑av| 18禁黄网站禁片午夜丰满| 亚洲三区欧美一区| 最近最新免费中文字幕在线| 国产蜜桃级精品一区二区三区 | 一进一出抽搐gif免费好疼 | 在线播放国产精品三级| 日韩 欧美 亚洲 中文字幕| 欧美一级毛片孕妇| 久久精品国产a三级三级三级| 免费在线观看黄色视频的| 午夜两性在线视频| 在线观看午夜福利视频| 水蜜桃什么品种好| 精品免费久久久久久久清纯 | 高清黄色对白视频在线免费看| 好看av亚洲va欧美ⅴa在| 亚洲av日韩精品久久久久久密| 亚洲精品久久午夜乱码| 成人免费观看视频高清| 女性被躁到高潮视频| 在线观看午夜福利视频| 男人舔女人的私密视频| 欧美一级毛片孕妇| 99re在线观看精品视频| 777米奇影视久久| 久热这里只有精品99| 黑人欧美特级aaaaaa片| 天天添夜夜摸| 99精品在免费线老司机午夜| 99精品欧美一区二区三区四区| 国产精品 欧美亚洲| 国产在线观看jvid| 久久久国产成人精品二区 | 另类亚洲欧美激情| 多毛熟女@视频| 国产真人三级小视频在线观看| 精品久久蜜臀av无| videos熟女内射| 最新在线观看一区二区三区| 乱人伦中国视频| 欧美成人免费av一区二区三区 | 一边摸一边抽搐一进一小说 | 亚洲av日韩在线播放| 又黄又粗又硬又大视频| 国产精品久久久久久人妻精品电影| 午夜福利影视在线免费观看| 亚洲av电影在线进入| 9191精品国产免费久久| 精品人妻熟女毛片av久久网站| 精品少妇久久久久久888优播| 美女高潮喷水抽搐中文字幕| 两个人看的免费小视频| 日本黄色日本黄色录像| 国产成+人综合+亚洲专区| 一边摸一边抽搐一进一小说 | 亚洲人成电影免费在线| 新久久久久国产一级毛片| 免费久久久久久久精品成人欧美视频| 国产一区二区三区综合在线观看| 国产成人免费无遮挡视频| 19禁男女啪啪无遮挡网站| 一级毛片精品| 久久国产精品大桥未久av| 久久久久国产精品人妻aⅴ院 | 久久久国产欧美日韩av| 成人影院久久| 极品教师在线免费播放| 国产精品亚洲一级av第二区| 美女 人体艺术 gogo| 久久亚洲精品不卡| 日日爽夜夜爽网站| 交换朋友夫妻互换小说| 三上悠亚av全集在线观看| 欧美黄色淫秽网站| 亚洲av美国av| 亚洲一区二区三区不卡视频| 嫁个100分男人电影在线观看| 婷婷精品国产亚洲av在线 | 国产午夜精品久久久久久| 欧美日韩瑟瑟在线播放| 欧美成狂野欧美在线观看| 精品久久久久久,| av天堂久久9| 黄色视频,在线免费观看| 亚洲精品成人av观看孕妇| 操出白浆在线播放| 精品人妻在线不人妻| 一级毛片精品| av超薄肉色丝袜交足视频| 午夜精品国产一区二区电影| 夫妻午夜视频| 亚洲av电影在线进入| 国产精品一区二区精品视频观看| 桃红色精品国产亚洲av| 老汉色∧v一级毛片| 国产极品粉嫩免费观看在线| 国产在线一区二区三区精| 亚洲成国产人片在线观看| 亚洲国产精品sss在线观看 | 俄罗斯特黄特色一大片| 欧美日韩精品网址| av福利片在线| 国产主播在线观看一区二区| 99国产精品99久久久久| 欧美在线一区亚洲| 999久久久国产精品视频| 亚洲精品国产一区二区精华液| 一二三四社区在线视频社区8| 老司机深夜福利视频在线观看| 最新的欧美精品一区二区| 好男人电影高清在线观看| 法律面前人人平等表现在哪些方面| 久久久水蜜桃国产精品网| 亚洲国产精品合色在线| 免费观看a级毛片全部| 色在线成人网| 亚洲欧美日韩高清在线视频| 久久人妻熟女aⅴ| 久久婷婷成人综合色麻豆| 婷婷丁香在线五月| 亚洲av成人av| 9热在线视频观看99| 欧美日韩福利视频一区二区| 欧美黑人精品巨大| 99国产精品99久久久久| 久久精品成人免费网站| 新久久久久国产一级毛片| 欧美大码av| 欧美日韩乱码在线| 亚洲中文字幕日韩| 欧美日韩一级在线毛片| av天堂在线播放| 亚洲情色 制服丝袜| 免费女性裸体啪啪无遮挡网站| 婷婷成人精品国产| 美女 人体艺术 gogo| 亚洲欧美色中文字幕在线| 制服人妻中文乱码| 欧美激情高清一区二区三区| 热re99久久精品国产66热6| 女警被强在线播放| 久久国产精品男人的天堂亚洲| 国产日韩一区二区三区精品不卡| 亚洲欧美一区二区三区黑人| 大型黄色视频在线免费观看| 高清在线国产一区| 狠狠婷婷综合久久久久久88av| 国产高清视频在线播放一区| 多毛熟女@视频| 欧美日韩亚洲综合一区二区三区_| 亚洲精品国产区一区二| 高清av免费在线| 国产精品自产拍在线观看55亚洲 | 操美女的视频在线观看| 十八禁网站免费在线| 美女国产高潮福利片在线看| 99久久国产精品久久久| 午夜福利,免费看| www日本在线高清视频| 美女 人体艺术 gogo| 亚洲男人天堂网一区| 欧美激情久久久久久爽电影 | 中文字幕人妻熟女乱码| 日韩欧美一区视频在线观看| 久久精品亚洲精品国产色婷小说| 99精品久久久久人妻精品| 国产av一区二区精品久久| 一本大道久久a久久精品| 国产精品美女特级片免费视频播放器 | 久久久久久亚洲精品国产蜜桃av| 日韩欧美一区视频在线观看| 精品国产乱子伦一区二区三区| 免费观看人在逋| 九色亚洲精品在线播放| 精品熟女少妇八av免费久了| tube8黄色片|