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

    基于諧振子遺傳算法的高效地下水優(yōu)化管理模型

    2015-03-07 03:36:45吳劍鋒施小清吳吉春
    關(guān)鍵詞:諧振子算例含水層

    吳 鳴,吳劍鋒,施小清,劉 杰,陳 干,吳吉春

    1.南京大學(xué)地球科學(xué)與工程學(xué)院水科學(xué)系,南京 210023 2.北京大學(xué)工學(xué)院水資源研究中心,北京 100871

    ?

    基于諧振子遺傳算法的高效地下水優(yōu)化管理模型

    吳 鳴1,吳劍鋒1,施小清1,劉 杰2,陳 干1,吳吉春1

    1.南京大學(xué)地球科學(xué)與工程學(xué)院水科學(xué)系,南京 210023 2.北京大學(xué)工學(xué)院水資源研究中心,北京 100871

    在傳統(tǒng)遺傳算法和模擬諧振子算法的基礎(chǔ)上,結(jié)合兩者的優(yōu)點(diǎn),提出了一種新型快速高效的諧振子遺傳算法。通過一個(gè)理想的水資源管理模型的算例和一個(gè)華北平原典型區(qū)地下水資源優(yōu)化的實(shí)際算例,從尋優(yōu)結(jié)果和尋優(yōu)效率兩個(gè)方面對(duì)諧振子遺傳算法、傳統(tǒng)遺傳算法和模擬諧振子算法進(jìn)行了對(duì)比分析。在兩個(gè)地下水管理模型中,與傳統(tǒng)的遺傳算法和模擬諧振子算法相比,新型的諧振子遺傳算法搜索效率達(dá)到模擬諧振子算法搜索效率的2倍以上,得到的最優(yōu)解比遺傳算法所得到的最優(yōu)解分別增加供水量1.1×103m3/d和0.47×108m3/a,說明諧振子遺傳算法具有更強(qiáng)的全局搜索能力和更好的尋優(yōu)效率。

    諧振子遺傳算法;遺傳算法;模擬諧振子算法;地下水管理模型;全局搜索能力;尋優(yōu)效率;華北平原

    0 引言

    地下水管理模型的求解是地下水資源科學(xué)管理的重要問題。地下水管理模型是尋求最少的社會(huì)經(jīng)濟(jì)環(huán)境成本來滿足既定的管理目標(biāo)[1-2]。長(zhǎng)期以來用于地下水管理模型求解的方法如單純形法、最速下降法、高斯-牛頓法等都是針對(duì)連續(xù)的、可導(dǎo)的目標(biāo)函數(shù)的,處理的問題比較簡(jiǎn)單,優(yōu)化結(jié)果的質(zhì)量強(qiáng)烈依賴初始解。由于遺傳算法(genetic algorithm, GA)[3-6]、禁忌搜索算法[7]、模擬退火算法[8]、粒子群優(yōu)化算法[2,9]、蟻群優(yōu)化算法[10]、模擬諧振子算法(simulated harmonic oscillator algorithm, SHOA)[11-13]等人工智能算法不受優(yōu)化問題滿足連續(xù)性和可導(dǎo)性的限制,且具有一定的全局尋優(yōu)能力,因而逐漸被應(yīng)用到地下水優(yōu)化管理模型求解的領(lǐng)域。然而地下水管理模型常常表現(xiàn)出多維、非線性等復(fù)雜特征,因此開發(fā)針對(duì)地下水優(yōu)化管理問題的高效算法十分必要。

    筆者針對(duì)遺傳算法和模擬諧振子算法的不足,結(jié)合傳統(tǒng)的模擬諧振子算法和遺傳算法的優(yōu)點(diǎn),首次提出了一種新型的諧振子遺傳算法(harmonic oscillator genetic algorithm, HOGA),將其與地下水管理模型相結(jié)合;并將HOGA分別應(yīng)用于理想算例和華北平原典型區(qū)地下水資源管理模型的算例中進(jìn)行對(duì)比,以期了解其在地下水管理模型方面的應(yīng)用前景。

    1 諧振子遺傳算法

    1.1 諧振子遺傳算法的基本原理

    遺傳算法屬于進(jìn)化算法的一種,它通過模仿自然界的選擇與遺傳的機(jī)理來尋找最優(yōu)解,McKiney和Lin[3]首先將其應(yīng)用到地下水管理模型的求解中,取得了很好的效果。遺傳算法具有良好的全局搜索能力,可以快速地搜索解空間,從而不會(huì)陷入局部最優(yōu)解的快速下降陷阱。但是遺傳算法的局部搜索能力較差,導(dǎo)致單純的遺傳算法比較費(fèi)時(shí),在進(jìn)化后期搜索效率較低。在實(shí)際應(yīng)用中,遺傳算法容易產(chǎn)生早熟收斂的問題。采用何種選擇方法既要使優(yōu)良個(gè)體得以保留,又要維持群體的多樣性,一直是遺傳算法中較難解決的問題。

    模擬諧振子算法(SHOA)是由模擬自然諧振子的物理規(guī)律而得到的一種通用的隨機(jī)搜索算法[11]。對(duì)于圖1所示的諧振子,當(dāng)彈簧質(zhì)點(diǎn)運(yùn)動(dòng)到左右端點(diǎn)位置時(shí),彈性勢(shì)能達(dá)到最大,等于系統(tǒng)的總能量;當(dāng)彈簧質(zhì)點(diǎn)處于振源位置時(shí),彈性勢(shì)能最小(等于0)。將彈簧質(zhì)點(diǎn)的勢(shì)能狀態(tài)與優(yōu)化問題中的目標(biāo)函數(shù)相對(duì)應(yīng),當(dāng)彈簧質(zhì)點(diǎn)到達(dá)振源位置時(shí),勢(shì)能處于最小值的基態(tài),對(duì)應(yīng)的目標(biāo)函數(shù)值即為找到的最優(yōu)解。SHOA既具有良好的局部搜索能力,又具有較好的全局搜索能力。它在整個(gè)解空間范圍內(nèi)廣泛地隨機(jī)搜索近似最優(yōu)解,而且以一定的規(guī)則接受差解,使領(lǐng)域解在局部山峰間平行跳躍,從而有效地控制了算法在全局范圍內(nèi)搜索的隨意性[11-13],但因此也給算法的尋優(yōu)效率帶來了問題。針對(duì)遺傳算法和模擬諧振子算法的不足,筆者綜合模擬諧振子算法和遺傳算法的優(yōu)點(diǎn),提出了一種新的算法----諧振子遺傳算法(HOGA)。HOGA利用模擬諧振子算法中的能級(jí)接受準(zhǔn)則,對(duì)遺傳算法中種群的選擇策略進(jìn)行了改進(jìn),同時(shí)利用模擬諧振子的領(lǐng)域解在局部山峰間的平行跳躍能力,有效地控制了遺傳算法的早熟問題。

    圖1 簡(jiǎn)諧振動(dòng)諧振子示意圖Fig.1 Oscillator of a simple harmonic oscillation

    1.2 諧振子遺傳算法求解地下水管理模型的步驟

    設(shè)諧振子遺傳算法求解的是最大值問題,則求解地下水管理模型的步驟為:

    1)隨機(jī)產(chǎn)生種群規(guī)模為p的初始群體P(m),m為進(jìn)化代數(shù),同時(shí)確定諧振子的振幅A,初始步長(zhǎng)LO、基態(tài)步長(zhǎng)LS。

    2)調(diào)用水流模型/運(yùn)移模型計(jì)算適應(yīng)度值。

    3)根據(jù)群體適應(yīng)度的最優(yōu)值和最差值建立諧振子模型,確定諧振子端點(diǎn)End和振源Init。振子端點(diǎn)End處的勢(shì)能ULO對(duì)應(yīng)適應(yīng)度的最差值f(End);基態(tài)步長(zhǎng)處的勢(shì)能為ULS,對(duì)應(yīng)為適應(yīng)度的最優(yōu)值f(LS)。通過最優(yōu)值和最差值之差的絕對(duì)值來確定諧振子的彈性系數(shù)k:

    (1)

    (2)

    求出k之后,就可以確定端點(diǎn)處勢(shì)能和基態(tài)步長(zhǎng)處的勢(shì)能:

    (3)

    (4)

    4)計(jì)算種群內(nèi)所有染色體對(duì)應(yīng)的能級(jí)Lx,對(duì)其中任意一條染色體,其適應(yīng)度值為f(Lx),則

    (5)

    (6)

    則其在諧振子中對(duì)應(yīng)的位置,也就是能級(jí)為

    (7)

    6)終止條件判斷,若不滿足終止條件,則m=m+1,形成下一代種群,轉(zhuǎn)移到第2)步,繼續(xù)進(jìn)行諧振子遺傳算法相關(guān)的能級(jí)判斷、選擇、交叉、變異等操作過程;若滿足終止條件,則輸出當(dāng)前最優(yōu)個(gè)體,算法結(jié)束。

    2 算例應(yīng)用

    2.1 理想算例

    在某潛水含水層上擬建一水源地,含水層中地下水為二維穩(wěn)定流[3],如圖2。含水層底板水平,其南部和北部邊界分別為水沼地和河流,其水位均為20 m,可作為第一類邊界;東部和西部邊界都為山體,作為隔水邊界。設(shè)含水層介質(zhì)均質(zhì)各向同性,滲透系數(shù)K=50 m/d,單位面積上的補(bǔ)給量R=0.001 m/d,共有10眼擬開采井,假設(shè)每口井的允許開采量為7 000 m3/d,求滿足該水源地含水層不致疏干情況下的最大可供水量(不考慮水井本身的釋水能力)。 管理模型可表示為:

    (8)

    約束條件為:

    (9)

    (10)

    其中:Qi為單井的抽水量(m3/d);hi為水位(m)。

    a.平面圖; b.剖面圖。根據(jù)文獻(xiàn)[3]修改。圖2 含水層結(jié)構(gòu)圖Fig.2 Structure of aquifer view

    Table 1 Total pumping rates of groundwater optimization based on different algorithms

    對(duì)比項(xiàng)目?jī)?yōu)化方法LPGASHOAHOGA最優(yōu)解/(103m3/d)59.358.058.959.1達(dá)到最優(yōu)解的計(jì)算代數(shù)7010146計(jì)算時(shí)間/min116.3170.279.7

    LP為線性規(guī)劃算法。LP最優(yōu)解引自文獻(xiàn)[3]。

    圖3 不同算法在最大化抽水量算例中的求解進(jìn)化過程Fig.3 Evolution of the total pumping rates over the generations of the different algorithms

    2.2 華北平原典型區(qū)地下水資源優(yōu)化算例

    2.2.1 研究區(qū)背景及水文地質(zhì)條件概述

    華北平原主要是指黃河以北、西倚太行山、東臨渤海的廣大平原地區(qū);行政區(qū)劃分屬北京市、天津市和河北省,并涉及山東、河南兩省的部分地域,面積約13.9 萬(wàn)km2,是我國(guó)北方重要的糧棉產(chǎn)地。該區(qū)大部分地區(qū)特別是大中城市處于地下水超采或嚴(yán)重超采狀態(tài),水資源匱乏和水環(huán)境惡化已經(jīng)成為制約社會(huì)經(jīng)濟(jì)發(fā)展的主要因素。為此,Liu等[14]以石家莊市(包括城區(qū)和7個(gè)郊縣)作為典型地區(qū)建立了該區(qū)地下水的水流模型,該水流模型不但刻畫了研究區(qū)開采前的天然狀態(tài)(補(bǔ)排均衡階段),而且再現(xiàn)了研究區(qū)不同開采時(shí)期(包括漏斗形成和發(fā)展階段[14])的水位變化過程,具有比較高的可信度,并從水量的角度建立了該區(qū)基于遺傳算法的地下水資源可持續(xù)利用的管理模型。筆者采用HOGA對(duì)Liu等[14]建立的研究區(qū)地下水管理模型進(jìn)行求解,并與GA和SHOA的求解效率進(jìn)行對(duì)比,進(jìn)一步驗(yàn)證HOGA在地下水管理模型中的可行性。

    X、Y分別表示東西方向和南北方向的距離。圖4 研究區(qū)位置及水文地質(zhì)概念模型Fig.4 Location of the study area and the hydrogeological conceptual model

    據(jù)文獻(xiàn)[14]。圖5 含水層地質(zhì)剖面示意圖Fig.5 Hydrogeologic cross section of the study area

    研究區(qū)包括石家莊市市區(qū)及正定、灤城、無極、藁城、趙縣、元氏和晉州7縣的部分區(qū)域(圖4),位于太行山南麓,屬于滹沱河沖洪積扇的中上部地帶,地勢(shì)西高東低,由西向東緩慢傾斜,面積約為4 000 km2。研究區(qū)處于滹沱河以黃壁莊為頂點(diǎn)的沖積扇的中上部,地下水主要賦存于第四系松散巖層孔隙中,自上而下可分為:全新統(tǒng)(Q4)、上更新統(tǒng)(Q3)、中更新統(tǒng)(Q2)和下更新統(tǒng)(Q1)共4個(gè)含水層組。上部的兩個(gè)含水層組之間無連續(xù)的隔水層,水力聯(lián)系密切,統(tǒng)稱為淺部含水層;下部的兩個(gè)含水層組巖性主要為含黏土卵礫石、礫石夾土層,分選性和富水性都很差,統(tǒng)稱深部含水層,如圖5所示。從總體上看:研究區(qū)含水層在水平方向上,由西向東含水層總厚度由薄變厚,粒度由粗變細(xì),層次由少增多,單層厚度變薄,分選性由差變好,富水性由強(qiáng)變?nèi)?;在垂直方向上,含水層上下部粒度較細(xì),厚度較小,中部粒度較粗,厚度也較大。隨著開采規(guī)模的不斷擴(kuò)大,實(shí)際上以上各含水層之間都存在著密切的水力聯(lián)系。

    地下水的補(bǔ)給主要有大氣降水、地表水的入滲、側(cè)向補(bǔ)給及農(nóng)田灌溉回歸水等。1965年前,石家莊市區(qū)內(nèi)地下水處于天然狀態(tài)下,由西或西北方向往東或東南方向流動(dòng),且自西向東,含水層介質(zhì)顆粒由粗變細(xì),含水單層由厚變薄,含水層數(shù)由少變多,水量和水力坡降由大變小,補(bǔ)給徑流條件和水質(zhì)由好變差,具典型的平原水文地質(zhì)特點(diǎn)。之后隨著地下水開采量的增大,地下水的排泄方式逐漸轉(zhuǎn)變?yōu)橐匀斯ら_采為主,其次為向下游的側(cè)向徑流排泄。隨著市區(qū)地下水降落漏斗范圍的擴(kuò)大和地下水位埋藏深度的增加,人工開采幾乎已成為目前唯一的地下水排泄方式[14-19]。

    2.2.2 研究區(qū)地下水模擬-優(yōu)化管理模型

    諸多研究[13-17]表明研究區(qū)淺層地下水和深層地下水之間具有統(tǒng)一的水力聯(lián)系,因此整個(gè)含水系統(tǒng)可以概化為單層的非均質(zhì)各向同性的二維水流運(yùn)動(dòng)。在邊界上,含水系統(tǒng)接受外界側(cè)向補(bǔ)給或往外排泄;在裸露區(qū)接受大氣降水入滲或農(nóng)業(yè)灌溉回滲,同時(shí)在區(qū)域內(nèi)最主要的排泄方式為人工開采。1960年以前,研究區(qū)地下水流處于天然狀態(tài),人工開采可忽略不計(jì),地下水流可近似認(rèn)為做穩(wěn)定流動(dòng);1960年以后,隨著地下水開采規(guī)模的加大,研究區(qū)地下水流不再處于天然狀況,而是對(duì)應(yīng)于不同邊界條件和源匯項(xiàng)的非穩(wěn)定運(yùn)動(dòng)。相應(yīng)地,描述研究區(qū)地下水流運(yùn)動(dòng)的微分方程可表示為

    (11)

    其中:h為水頭(m);t為時(shí)間(d);B為含水層的底板標(biāo)高(m);qS為源匯項(xiàng)(d-1);μ為給水度(無量綱)。

    再加上相應(yīng)的初始條件,這樣,微分方程與邊界條件和初始條件一起就構(gòu)成了描述研究區(qū)地下水流運(yùn)動(dòng)的數(shù)學(xué)模型[13]。模擬時(shí)間從1959年到2004年,總共63個(gè)應(yīng)力期;模型在空間上剖分成130行和135列,每個(gè)網(wǎng)格的邊長(zhǎng)均為500 m。用MODFLOW[20]對(duì)模型進(jìn)行求解,詳細(xì)的水流模型參見文獻(xiàn)[14]。

    基于以上地下水模擬模型的最大化開采量的管理模型可表述為

    (12)

    約束條件:

    (13)

    (14)

    式中:Qmax表示抽水井i的最大出水能力(m3/d);hj表示約束點(diǎn)j處的水位(m);hmin表示在點(diǎn)j處的水位最低約束值(m);NT為管理期的數(shù)目;NH,NW分別表示水位約束點(diǎn)和抽水井的數(shù)目;Δtk表示每個(gè)管理期的時(shí)長(zhǎng)(d)。

    當(dāng)然,管理模型還要耦合上述地下水模擬模型(微分方程(11)及其定解條件),可簡(jiǎn)單地表示為

    (15)

    式中:hk表示研究區(qū)地下水在k時(shí)段的水位分布;Transh代表水頭由k時(shí)段到k+1時(shí)段的狀態(tài)轉(zhuǎn)移函數(shù);T為時(shí)段總數(shù)。式(15)表示系統(tǒng)在(k+1)時(shí)刻的狀態(tài)變量h是其k時(shí)刻狀態(tài)變量h及其決策變量Q兩者的函數(shù)。

    2.2.3 地下水優(yōu)化管理模型的求解

    求解管理模型時(shí),GA和HOGA最大代數(shù)均設(shè)為100,種群數(shù)設(shè)為1 000,為便于比較,SHOA目標(biāo)函數(shù)每計(jì)算1 000個(gè)解后統(tǒng)計(jì)一次最優(yōu)結(jié)果,搜索新解的上限為10萬(wàn)個(gè),采用GA、SHOA和HOGA分別求解的優(yōu)化過程及結(jié)果如圖6,最大抽水量、地下水超采量、收斂代數(shù)和計(jì)算時(shí)間的結(jié)果如表2。

    圖6 不同算法在華北平原地下水管理模型優(yōu)化算例中的求解進(jìn)化過程Fig.6 Evolution of the total pumping rates over the generations of the different algorithms in North China Plain

    GA和HOGA都在約80代達(dá)到最優(yōu)解,而SHOA則在第89代達(dá)到最優(yōu)解,表明HOGA的搜索速度和計(jì)算效率繼承了遺傳算法的優(yōu)點(diǎn),比原來的SHOA有了很大的提高。GA最優(yōu)解為10.22×108m3/a,和2004年實(shí)際開采量[13]的14.97×108m3/a對(duì)比,地下水超采量達(dá)到4.75×108m3/a;SHOA最優(yōu)解為10.29×108m3/a,地下水超采量達(dá)4.67×108m3/a;HOGA的最優(yōu)解為10.69×108m3/a,地下水超采量為4.27×108m3/a,相比GA增加供水量0.47×108m3/a,相比SHOA增加供水量0.40×108m3/a,地下水超采量得到較大降低。由此表明HOGA的全局搜索能力有了極大的提高,均強(qiáng)于GA和SHOA的搜索能力。

    表2 華北平原地下水管理模型優(yōu)化結(jié)果

    Table 2 Result of North China Plain optimization model by different algorithms

    對(duì)比項(xiàng)目?jī)?yōu)化方法GASHOAHOGA最優(yōu)解/(108m3/a)10.2210.2910.69地下水超采量/(108m3/a)4.754.674.27達(dá)到最優(yōu)解的計(jì)算代數(shù)808980計(jì)算時(shí)間/h43.3588.9743.57

    3 結(jié)論

    針對(duì)模擬諧振子算法和遺傳算法的不足,提出一種新型的諧振子遺傳算法。諧振子系統(tǒng)勢(shì)能收斂到基態(tài)所特有的物理性質(zhì)能保證諧振子遺傳算法最后的收斂性,并通過諧振子對(duì)新解的接受準(zhǔn)則來對(duì)使種群朝著最優(yōu)的方向進(jìn)化,使得諧振子遺傳算法優(yōu)化尋優(yōu)結(jié)果的質(zhì)量和計(jì)算效率均有很大提高。

    諧振子遺傳算法應(yīng)用于地下水管理模型的求解結(jié)果表明:對(duì)于理想算例,諧振子遺傳算法表現(xiàn)出了很強(qiáng)的全局尋優(yōu)能力和搜索效率;在華北平原典型區(qū)地下水資源的優(yōu)化算例中,諧振子遺傳算法全局尋優(yōu)能力在很大程度上超過了遺傳算法和模擬諧振子算法,相對(duì)于模擬諧振子算法,搜索效率也有很大的提高。表明諧振子遺傳算法不僅對(duì)理想算例的尋優(yōu)具有很好的效果,而且對(duì)區(qū)域性復(fù)雜算例的地下水資源優(yōu)化問題也有很強(qiáng)的優(yōu)勢(shì)。

    諧振子遺傳算法融合了遺傳算法和模擬諧振子算法的優(yōu)點(diǎn),對(duì)初始解沒有太大的依賴性,通用性較強(qiáng),可推廣應(yīng)用于水文地質(zhì)參數(shù)識(shí)別、地下水污染源識(shí)別、地下水污染修復(fù)等優(yōu)化管理領(lǐng)域。

    [1] 楊蘊(yùn),吳劍鋒,吳吉春. 兩種智能算法在求解地下水管理模型中的對(duì)比[J].吉林大學(xué)學(xué)報(bào):地球科學(xué)版, 2009, 39(3):474-502. Yang Yun, Wu Jianfeng, Wu Jichun. A Comparative Study of Two Intelligent Optimization Techniques for Groundwater Management Modeling[J]. Journal of Jilin University:Earth Science Edition, 2009, 39(3): 474-502.

    [2] 駱乾坤,王佩,朱國(guó)榮. 水文地質(zhì)參數(shù)識(shí)別的快速和諧搜索算法[J].水文地質(zhì)工程地質(zhì), 2011, 38(4):14-19. Luo Qiankun,Wang Pei, Zhu Guorong. Fast Harmony Search Algorithm and Its Application to Hydrogeological Parameters Identification[J]. Hydrogeology & Engineering Geology, 2011, 38(4):14-19.

    [3] McKinney D C, Lin M D. Genetic Algorithm Solution of Groundwater Management Models[J]. Water Resources Research, 1994, 30(6): 1897-1906.

    [4] 吳鳴,吳劍鋒,林錦,等.地下油罐泄漏區(qū)污染源的自動(dòng)識(shí)別[J].環(huán)境科學(xué)學(xué)報(bào), 2013, 33(12): 3251-3259. Wu Ming, Wu Jianfeng, Lin Jin, et al. Automated Identification of the Unknown Contaminant Source in Groundwater at a Leaking Underground Storage Tank Site[J]. Acta Scientiae Circumstantiae, 2013, 33(12): 3251-3259.

    [5] 吳劍鋒,朱學(xué)愚,錢家忠,等. GASAPF方法在徐州市裂隙巖溶水資源管理模型中的應(yīng)用[J]. 水利學(xué)報(bào), 2000, 31(12):7-13. Wu Jianfeng, Zhu Xueyu, Qian Jiazhong, et al. Application of GASAPF to Optimization Model for Fracture Karst Water Resources Management[J]. Journal of Hydraulic Engineering, 2000,31(12): 7-13.

    [6] 林錦,鄭春苗,吳劍鋒,等. 基于遺傳算法的變密度條件下地下水模擬優(yōu)化模型[J]. 水利學(xué)報(bào), 2007,38(10):1236-1244. Lin Jin, Zheng Chunmiao, Wu Jianfeng, et al. Ground Water Simulation Optimization Model Based on Genetic Algorithm Under Variable Density Conditions[J]. Journal of Hydraulic Engineering, 2007,38(10): 1236-1244.

    [7] 楊蘊(yùn),吳劍鋒,于軍,等. 基于參數(shù)不確定性的地下水污染治理多目標(biāo)管理模型[J]. 環(huán)境科學(xué)學(xué)報(bào), 2013, 33(7):2059-2067. Yang Yun, Wu Jianfeng, Yu Jun, et al. A Multi-Objective Simulation-Optimization Model for Optimal Design of Groundwater Remediation Systems Under Uncertainty[J]. Acta Scientiae Circumstantiae, 33(7): 2059-2067.

    [8] 高利平,何敬宇. 模擬退火算法在邊坡穩(wěn)定分析中的應(yīng)用研究[J].內(nèi)蒙古農(nóng)業(yè)大學(xué)學(xué)報(bào), 2010, 31(4):218-221. Gao Liping, He Jingyu. The Application and Investigation of Simulated Annealing Method in the Analysis of Slope Stability[J]. Journal of Inner Mongolia Agricultural University, 2010, 31(4): 218-221.

    [9] 劉賀,張弘強(qiáng),劉斌. 基于粒子群優(yōu)化神經(jīng)網(wǎng)絡(luò)算法的深基坑變形預(yù)測(cè)方法[J].吉林大學(xué)學(xué)報(bào):地球科學(xué)版, 2014, 44(5):1609-1614. Liu He, Zhang Hongqiang, Liu Bin. A Prediction Method or the Deformation of Deep Foundation Pit Based on the Partical Swarm Optimization Neural Network[J]. Journal of Jilin University: Earth Science Edition, 2014,44(5): 1609-1614.

    [10] 秦永波,王鵬,肖黎彬,等. 量子諧振子蟻群算法[J]. 計(jì)算機(jī)應(yīng)用, 2011, 31(2):54-69. Qin Yongbo, Wang Peng, Xiao Libin, et al. Ant Colony Optimization of Quantum Harmonic Oscillators[J]. Journal of Computer Applications, 2011, 31(2): 54-69.

    [11] 倪霖,段超,鐘輝. 基于模擬諧振子算法的多項(xiàng)目調(diào)度[J]. 計(jì)算機(jī)應(yīng)用, 2011,31(9):2559-2562. Ni Lin, Duan Chao, Zhong Hui. Multi-Project Scheduling Based on Simulated Harmonic Oscillator Algorithm[J]. Journal of Computer Applications, 2011, 31(9): 2559-2562.

    [12] 王培崇,錢旭. 模擬諧振子算法及其全局收斂性分析[J]. 計(jì)算機(jī)工程, 2013,39(3):209-212. Wang Peichong, Qian Xu. Simulated Harmonic Oscillator Algorithm and Its Global Convergence Analysis[J]. Computer Engineering, 2013,39(3): 209-212.

    [13] 程勖,李文輝,劉裕斌. 基于模擬諧振子算法的服務(wù)調(diào)度技術(shù)[J].大連海事大學(xué)學(xué)報(bào). 2013, 39(2):78-81. Cheng Xu, Li Wenhui, Liu Yubin. Service Scheduling Technique Based on Simulated Harmonic Oscillator Algorithm[J]. Journal of Dalian Maritime University, 2013,39(2): 78-81.

    [14] Liu J, Zheng C M, Li Z, et al. Ground Water Sustainability: Methodology and Application to the North China Plain[J]. Ground Water, 2008, 46(6):897-909.

    [15] Cao G L, Zheng C M, Scanlon B R, et al. Use of Flow Modeling to Assess Sustainability of Groundwater Resources in the North China Plain[J]. Water Resources Research, 2013, 49: 1-17.

    [16] Wang B G, Jin M G, Nimmo J R, et al. Estimating Groundwater Recharge in Hebei Plain, China Under Varying Land Use Practices Using Tritium and Bromide Tracers[J]. Journal of Hydrology, 2008, 356:209-222.

    [17] 張光輝,費(fèi)宇紅,張行南,等. 滹沱河流域平原區(qū)地下水流場(chǎng)異常變化與原因[J].水利學(xué)報(bào). 2008, 39(6):747-752. Zhang Guanghui, Fei Yuhong, Zhang Xingnan, et al. Abnormal Variation of Groundwater Flow Field in Plain Area of Hutuo River Basin and Analysis on Its Cause[J]. Journal of Hydraulic Engineering, 2008, 39(6):747-752.

    [18] 馮慧敏,張光輝,王電龍,等.近50年來石家莊地區(qū)地下水流場(chǎng)演變驅(qū)動(dòng)力分析[J].水利學(xué)報(bào). 2014, 45(2):180-186. Feng Huimin, Zhang Guanghui, Wang Dianlong, et al. Analysis on Driving Force for Groundwater Flow Field Evolution in Shijiazhuang Area in Recent 50 Years[J]. Journal of Hydraulic Engineering, 2014, 45(2):180-186.

    [19] 費(fèi)宇紅,張光輝,劉克巖,等.華北平原河道地下調(diào)蓄與利用潛力[J].吉林大學(xué)學(xué)報(bào):地球科學(xué)版, 2011, 41(增刊1):265-271. Fei Yuhong, Zhang Guanghui, Liu Keyan, et al. The Adjustment of Groundwater Storage and Utilization Potential of River Channel in North China Plain[J]. Journal of Jilin University: Earth Science Edition, 2011, 41(Sup.1):265-271.

    [20] Harbaugh A W. MODFLOW-2005, The US Geological Survey Modular Ground-Water Model:The Ground-Water Flow Process[R]. Reston: U S Geological Survey, 2005.

    New Harmonic Oscillator Genetic Algorithm for Efficient Groundwater Optimization and Management

    Wu Ming1, Wu Jianfeng1, Shi Xiaoqing1, Liu Jie2, Chen Gan1, Wu Jichun1

    1.DepartmentofHydrosciences,SchoolofEarthSciencesandEngineering,NanjingUniversity,Nanjing210023,China2.CenterforWaterResearch,CollegeofEngineering,PekingUniversity,Beijing100871,China

    Combining the advantages of conventional genetic algorithm (GA) with simulated harmonic oscillator algorithm (SHOA), we put forward a new efficient harmonic oscillator genetic algorithm (HOGA). Based on the application to an ideal test and an optimization test for the groundwater management model of a typical region in North China Plain (NCP), the comparisons among GA, SHOA and HOGA are performed in respect both of the quality of optimal solutions and the efficiency of algorithmic optimization. Compared with GA and SHOA, HOGA achieves an efficiency over two times higher than SHOA. The best pumping schemes achieved by HOGA can increase water supply by 1.1×103m3/d and 0.47×108m3/a compared with the pumping schemes achieved by GA in the two groundwater management models. This reveals that HOGA has the best global searching ability with the highest optimization efficiency. The results suggest that the HOGA presented in this study be of extensive perspectives in application to multi-dimensional and nonlinear groundwater management and optimization.

    harmonic oscillator genetic algorithm; genetic algorithm; simulated harmonic oscillator algorithm; groundwater management model; global searching ability; optimization efficiency; North China Plain (NCP)

    10.13278/j.cnki.jjuese.201505203.

    2014-12-12

    國(guó)家自然科學(xué)基金資助項(xiàng)目(41372235, 41072175);國(guó)家“973”計(jì)劃項(xiàng)目(2010CB428803)

    吳鳴(1989--),男,博士研究生,主要從事地下水模擬優(yōu)化管理研究,E-mail:mysticwilliamnanjing@gmail.com

    吳劍鋒(1971--),男,教授,博士,主要從事水資源模擬優(yōu)化管理研究,E-mail:jfwu@nju.edu.cn。

    10.13278/j.cnki.jjuese.201505203

    P641; X523

    A

    吳鳴,吳劍鋒,施小清,等. 基于諧振子遺傳算法的高效地下水優(yōu)化管理模型.吉林大學(xué)學(xué)報(bào):地球科學(xué)版,2015,45(5):1485-1492.

    Wu Ming, Wu Jianfeng, Shi Xiaoqing, et al.New Harmonic Oscillator Genetic Algorithm for Efficient Groundwater Optimization and Management.Journal of Jilin University:Earth Science Edition,2015,45(5):1485-1492.doi:10.13278/j.cnki.jjuese.201505203.

    猜你喜歡
    諧振子算例含水層
    全球多個(gè)含水層里的水正快速流失
    諧振子支柱偏心誤差對(duì)諧振子振動(dòng)特性影響分析(英文)
    美國(guó)西部奧加拉拉含水層水位下降原因初探
    基于振蕩能量的低頻振蕩分析與振蕩源定位(二)振蕩源定位方法與算例
    互補(bǔ)問題算例分析
    全球地下含水層下降驚人:要被抽干了
    基于CYMDIST的配電網(wǎng)運(yùn)行優(yōu)化技術(shù)及算例分析
    巖溶含水層水流模型研究進(jìn)展
    燃煤PM10湍流聚并GDE方程算法及算例分析
    三維各向異性耦合諧振子體系的非形式性嚴(yán)格波函數(shù)
    在线精品无人区一区二区三| 大陆偷拍与自拍| 久久国产亚洲av麻豆专区| 我的亚洲天堂| av网站在线播放免费| 一本色道久久久久久精品综合| 亚洲一卡2卡3卡4卡5卡精品中文| 69av精品久久久久久 | 欧美黄色片欧美黄色片| 91麻豆精品激情在线观看国产 | 黄色视频不卡| 黄片小视频在线播放| 久热这里只有精品99| 男人添女人高潮全过程视频| 亚洲av成人不卡在线观看播放网 | 欧美在线一区亚洲| 午夜福利免费观看在线| 日韩,欧美,国产一区二区三区| 欧美激情 高清一区二区三区| 亚洲精品一二三| 午夜福利一区二区在线看| 一区二区日韩欧美中文字幕| 18禁国产床啪视频网站| 亚洲男人天堂网一区| 亚洲男人天堂网一区| 在线天堂中文资源库| 性高湖久久久久久久久免费观看| 国产欧美亚洲国产| 成人亚洲精品一区在线观看| 免费女性裸体啪啪无遮挡网站| 电影成人av| 桃红色精品国产亚洲av| 黄片大片在线免费观看| 亚洲精品久久久久久婷婷小说| 欧美日韩成人在线一区二区| 后天国语完整版免费观看| 国产国语露脸激情在线看| 十分钟在线观看高清视频www| 自拍欧美九色日韩亚洲蝌蚪91| 少妇的丰满在线观看| 久久99一区二区三区| 黄色 视频免费看| 蜜桃国产av成人99| 伊人久久大香线蕉亚洲五| 欧美 日韩 精品 国产| 国产精品一二三区在线看| 纯流量卡能插随身wifi吗| 乱人伦中国视频| 亚洲少妇的诱惑av| 亚洲专区中文字幕在线| 日韩欧美国产一区二区入口| 精品人妻一区二区三区麻豆| 啦啦啦免费观看视频1| 精品人妻在线不人妻| 亚洲第一青青草原| 中文字幕人妻丝袜一区二区| 亚洲精品粉嫩美女一区| 国产精品香港三级国产av潘金莲| 色视频在线一区二区三区| 爱豆传媒免费全集在线观看| 黄色片一级片一级黄色片| 欧美精品一区二区大全| 搡老乐熟女国产| 伦理电影免费视频| 欧美大码av| 欧美精品av麻豆av| 淫妇啪啪啪对白视频 | 亚洲国产精品成人久久小说| 在线十欧美十亚洲十日本专区| 亚洲国产精品成人久久小说| 日韩中文字幕欧美一区二区| 一边摸一边做爽爽视频免费| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产成人一精品久久久| 波多野结衣av一区二区av| 国产亚洲一区二区精品| 亚洲少妇的诱惑av| 国产一区有黄有色的免费视频| 人妻 亚洲 视频| 男女床上黄色一级片免费看| 亚洲av日韩精品久久久久久密| 国产伦人伦偷精品视频| 欧美精品人与动牲交sv欧美| 中文欧美无线码| 国产国语露脸激情在线看| 婷婷丁香在线五月| 中文字幕最新亚洲高清| 国产av国产精品国产| 欧美日韩av久久| 秋霞在线观看毛片| videos熟女内射| 免费在线观看完整版高清| 久久国产精品男人的天堂亚洲| 国产精品99久久99久久久不卡| 精品少妇久久久久久888优播| 99热国产这里只有精品6| 精品福利永久在线观看| 宅男免费午夜| 精品人妻一区二区三区麻豆| av免费在线观看网站| 中文字幕av电影在线播放| 日本av免费视频播放| 欧美在线一区亚洲| 欧美日韩中文字幕国产精品一区二区三区 | 日日摸夜夜添夜夜添小说| 国产成人欧美在线观看 | www.av在线官网国产| 午夜福利在线观看吧| 色婷婷av一区二区三区视频| 国产精品一区二区免费欧美 | 国产精品国产av在线观看| 大陆偷拍与自拍| 免费观看a级毛片全部| 欧美另类亚洲清纯唯美| 亚洲精品久久午夜乱码| 久久久久网色| 久久久久久久国产电影| 美女午夜性视频免费| 亚洲视频免费观看视频| 久久久国产精品麻豆| 纵有疾风起免费观看全集完整版| 亚洲色图 男人天堂 中文字幕| 欧美一级毛片孕妇| 亚洲国产欧美网| 国产视频一区二区在线看| 国产亚洲av片在线观看秒播厂| 亚洲精品一二三| 亚洲欧美一区二区三区黑人| 亚洲av国产av综合av卡| 成人黄色视频免费在线看| 制服人妻中文乱码| 人人妻人人澡人人看| 国产精品免费大片| 亚洲国产欧美网| 美女视频免费永久观看网站| 精品国产一区二区三区四区第35| 成年人免费黄色播放视频| 亚洲国产欧美一区二区综合| 亚洲伊人久久精品综合| 黄色怎么调成土黄色| 人妻人人澡人人爽人人| 免费一级毛片在线播放高清视频 | 亚洲综合色网址| 后天国语完整版免费观看| 欧美日韩亚洲综合一区二区三区_| 精品一品国产午夜福利视频| 侵犯人妻中文字幕一二三四区| 久久人人爽av亚洲精品天堂| 午夜精品久久久久久毛片777| 免费观看a级毛片全部| 亚洲一区二区三区欧美精品| 亚洲av欧美aⅴ国产| 狠狠婷婷综合久久久久久88av| 人人妻人人爽人人添夜夜欢视频| 另类亚洲欧美激情| 男人舔女人的私密视频| 中文字幕精品免费在线观看视频| 国产成人a∨麻豆精品| 亚洲精品国产色婷婷电影| 欧美一级毛片孕妇| 美女主播在线视频| 欧美日韩亚洲综合一区二区三区_| 久久国产精品影院| 久久久久国产精品人妻一区二区| 好男人电影高清在线观看| 日韩精品免费视频一区二区三区| av在线app专区| 亚洲国产欧美网| 久久久久久人人人人人| 国产又爽黄色视频| 老司机影院毛片| 欧美精品一区二区免费开放| www.精华液| 国产在视频线精品| 自线自在国产av| 宅男免费午夜| 99国产精品一区二区蜜桃av | 在线观看舔阴道视频| 老汉色∧v一级毛片| 久久久精品免费免费高清| www.av在线官网国产| a在线观看视频网站| 午夜久久久在线观看| 一本综合久久免费| 一边摸一边做爽爽视频免费| 丁香六月天网| 99国产精品免费福利视频| 亚洲一区二区三区欧美精品| 中文精品一卡2卡3卡4更新| 色播在线永久视频| 成人av一区二区三区在线看 | 欧美少妇被猛烈插入视频| 女人爽到高潮嗷嗷叫在线视频| 麻豆乱淫一区二区| 免费少妇av软件| 国产亚洲精品久久久久5区| 99热国产这里只有精品6| 国产成人av激情在线播放| 午夜两性在线视频| 一区二区日韩欧美中文字幕| 亚洲国产精品一区二区三区在线| 免费不卡黄色视频| 日日摸夜夜添夜夜添小说| 丰满饥渴人妻一区二区三| 亚洲欧美色中文字幕在线| 久久毛片免费看一区二区三区| 狠狠狠狠99中文字幕| 三级毛片av免费| 精品第一国产精品| 欧美日韩av久久| 新久久久久国产一级毛片| a级片在线免费高清观看视频| 叶爱在线成人免费视频播放| 男人操女人黄网站| 久久综合国产亚洲精品| 黑人巨大精品欧美一区二区mp4| 精品亚洲成a人片在线观看| 日韩熟女老妇一区二区性免费视频| tocl精华| 久久影院123| 又黄又粗又硬又大视频| 免费在线观看影片大全网站| 国产主播在线观看一区二区| 亚洲av电影在线进入| 视频区图区小说| 精品国内亚洲2022精品成人 | 日韩制服丝袜自拍偷拍| 女人被躁到高潮嗷嗷叫费观| 亚洲欧美日韩另类电影网站| 午夜福利在线观看吧| 国产免费现黄频在线看| 亚洲免费av在线视频| 视频区欧美日本亚洲| 国产成人系列免费观看| 视频在线观看一区二区三区| 高清黄色对白视频在线免费看| 久久中文字幕一级| 母亲3免费完整高清在线观看| 日韩制服骚丝袜av| 国产成人a∨麻豆精品| 久久久久久久大尺度免费视频| 亚洲精品一二三| 亚洲国产欧美一区二区综合| 国产精品欧美亚洲77777| 美女国产高潮福利片在线看| 欧美+亚洲+日韩+国产| 国产淫语在线视频| 满18在线观看网站| 亚洲精品乱久久久久久| 久久精品国产综合久久久| 不卡一级毛片| 亚洲成人国产一区在线观看| av网站免费在线观看视频| 久久久久国产一级毛片高清牌| av电影中文网址| 精品亚洲成国产av| 国产三级黄色录像| 久久天躁狠狠躁夜夜2o2o| 一级片'在线观看视频| xxxhd国产人妻xxx| 亚洲精品国产区一区二| 欧美黄色片欧美黄色片| 最新的欧美精品一区二区| 欧美国产精品一级二级三级| 老司机影院成人| 国产一区有黄有色的免费视频| 最近中文字幕2019免费版| 日本撒尿小便嘘嘘汇集6| 搡老熟女国产l中国老女人| 精品人妻1区二区| 久久国产精品大桥未久av| 桃花免费在线播放| 色播在线永久视频| 国产日韩欧美亚洲二区| 日韩熟女老妇一区二区性免费视频| 人人妻人人澡人人看| 午夜91福利影院| 亚洲全国av大片| 黑丝袜美女国产一区| 18禁裸乳无遮挡动漫免费视频| 中文字幕精品免费在线观看视频| 亚洲少妇的诱惑av| 国产日韩欧美亚洲二区| 亚洲午夜精品一区,二区,三区| 在线观看免费视频网站a站| 亚洲欧美精品自产自拍| 午夜91福利影院| 女警被强在线播放| 久久九九热精品免费| 欧美少妇被猛烈插入视频| 性色av一级| 一二三四社区在线视频社区8| tocl精华| 国产av国产精品国产| 亚洲国产成人一精品久久久| 岛国毛片在线播放| 一进一出抽搐动态| 一级毛片女人18水好多| 丝瓜视频免费看黄片| 国产精品av久久久久免费| 国产激情久久老熟女| 亚洲欧美精品综合一区二区三区| 国产又爽黄色视频| 97在线人人人人妻| 老鸭窝网址在线观看| 黑人欧美特级aaaaaa片| 亚洲国产欧美日韩在线播放| 亚洲九九香蕉| 欧美少妇被猛烈插入视频| 国产亚洲一区二区精品| 性高湖久久久久久久久免费观看| 一进一出抽搐动态| 亚洲精品国产av成人精品| av线在线观看网站| 中文字幕另类日韩欧美亚洲嫩草| 免费久久久久久久精品成人欧美视频| 美女视频免费永久观看网站| 欧美精品人与动牲交sv欧美| 91精品国产国语对白视频| 久久久久久久大尺度免费视频| 国产福利在线免费观看视频| 青青草视频在线视频观看| 天堂8中文在线网| 亚洲精品国产色婷婷电影| 国产精品免费视频内射| 在线永久观看黄色视频| 免费在线观看黄色视频的| 日日爽夜夜爽网站| 国产又爽黄色视频| 久久精品亚洲av国产电影网| 黄色a级毛片大全视频| 午夜精品国产一区二区电影| videosex国产| 欧美一级毛片孕妇| 久热爱精品视频在线9| 国产精品一区二区在线观看99| 精品一区在线观看国产| 韩国精品一区二区三区| a级毛片在线看网站| 国产成人精品久久二区二区免费| 亚洲精品一二三| 一本—道久久a久久精品蜜桃钙片| 亚洲成人免费电影在线观看| 丝袜美腿诱惑在线| 黄片播放在线免费| 少妇精品久久久久久久| 宅男免费午夜| av福利片在线| 精品高清国产在线一区| av在线app专区| 天天操日日干夜夜撸| 亚洲欧美日韩高清在线视频 | 国产精品99久久99久久久不卡| 久久精品国产a三级三级三级| 欧美变态另类bdsm刘玥| 一个人免费看片子| 欧美日韩亚洲综合一区二区三区_| tocl精华| 国产在视频线精品| 亚洲久久久国产精品| 国产男女内射视频| 另类精品久久| 久久久国产精品麻豆| 欧美中文综合在线视频| 国产亚洲一区二区精品| 免费高清在线观看日韩| 高清av免费在线| 中文字幕制服av| 交换朋友夫妻互换小说| 淫妇啪啪啪对白视频 | 国产精品.久久久| 日韩中文字幕视频在线看片| 婷婷成人精品国产| 满18在线观看网站| 国产成人精品久久二区二区免费| 在线永久观看黄色视频| 热99re8久久精品国产| 一级毛片女人18水好多| 亚洲伊人色综图| 午夜福利在线观看吧| 亚洲 欧美一区二区三区| 欧美97在线视频| av福利片在线| 色婷婷久久久亚洲欧美| 日本黄色日本黄色录像| 国产男女内射视频| 久久人人97超碰香蕉20202| 欧美精品人与动牲交sv欧美| 啦啦啦中文免费视频观看日本| 精品一品国产午夜福利视频| 操出白浆在线播放| 久久久久国产精品人妻一区二区| 久久中文字幕一级| 日本av手机在线免费观看| 精品人妻一区二区三区麻豆| 亚洲国产欧美一区二区综合| 色播在线永久视频| 国产伦理片在线播放av一区| 女人精品久久久久毛片| 中文字幕高清在线视频| 国产精品久久久人人做人人爽| 一进一出抽搐动态| a级毛片黄视频| 国产在线免费精品| 亚洲avbb在线观看| 日韩欧美国产一区二区入口| 国产福利在线免费观看视频| 国产免费一区二区三区四区乱码| 国产精品久久久人人做人人爽| 中文字幕人妻熟女乱码| 久久国产精品影院| 国产一区二区激情短视频 | 老熟女久久久| 亚洲国产欧美一区二区综合| 欧美久久黑人一区二区| 国产亚洲一区二区精品| 天堂中文最新版在线下载| 久久国产亚洲av麻豆专区| 国产免费av片在线观看野外av| 日韩一卡2卡3卡4卡2021年| 韩国高清视频一区二区三区| 狠狠婷婷综合久久久久久88av| 日本av手机在线免费观看| 久久亚洲精品不卡| 国产成人av教育| 97精品久久久久久久久久精品| 精品国产一区二区三区久久久樱花| 亚洲国产精品一区三区| 国产xxxxx性猛交| 麻豆国产av国片精品| 欧美日韩精品网址| 久久女婷五月综合色啪小说| 人妻人人澡人人爽人人| 日韩免费高清中文字幕av| 男女床上黄色一级片免费看| 法律面前人人平等表现在哪些方面 | 手机成人av网站| 老司机影院毛片| 看免费av毛片| 国产精品国产av在线观看| 亚洲国产av新网站| 97人妻天天添夜夜摸| 午夜福利免费观看在线| 伦理电影免费视频| 成人国产一区最新在线观看| 国产精品偷伦视频观看了| 精品福利观看| www.999成人在线观看| avwww免费| 精品久久久久久电影网| 欧美日韩国产mv在线观看视频| 国产欧美日韩一区二区三 | 精品一区在线观看国产| 午夜免费成人在线视频| 国产欧美日韩一区二区精品| 亚洲熟女毛片儿| 亚洲av欧美aⅴ国产| 精品少妇内射三级| 在线天堂中文资源库| 久久99热这里只频精品6学生| 两个人免费观看高清视频| 性少妇av在线| 妹子高潮喷水视频| 人成视频在线观看免费观看| 窝窝影院91人妻| 中文字幕人妻丝袜一区二区| 搡老乐熟女国产| 美女中出高潮动态图| 欧美日韩一级在线毛片| 一边摸一边抽搐一进一出视频| 性高湖久久久久久久久免费观看| 免费看十八禁软件| 亚洲精品乱久久久久久| 亚洲全国av大片| 国产一区二区三区av在线| 欧美日韩精品网址| 国产在线免费精品| 正在播放国产对白刺激| 免费黄频网站在线观看国产| 超色免费av| 99国产精品一区二区三区| 久久国产精品大桥未久av| 免费在线观看日本一区| 日本欧美视频一区| 少妇人妻久久综合中文| 女人被躁到高潮嗷嗷叫费观| 亚洲人成电影观看| 久久久久国内视频| 中文字幕av电影在线播放| 国产欧美日韩一区二区三 | 大片电影免费在线观看免费| 欧美+亚洲+日韩+国产| 久久精品亚洲av国产电影网| 一二三四在线观看免费中文在| 国产免费现黄频在线看| 久久人人爽av亚洲精品天堂| 久久中文看片网| 一级片'在线观看视频| 搡老岳熟女国产| 18禁裸乳无遮挡动漫免费视频| 国产精品国产三级国产专区5o| 色精品久久人妻99蜜桃| 老熟妇乱子伦视频在线观看 | 成人三级做爰电影| 午夜影院在线不卡| 国产极品粉嫩免费观看在线| 亚洲三区欧美一区| 亚洲av国产av综合av卡| 两性午夜刺激爽爽歪歪视频在线观看 | 狂野欧美激情性xxxx| 亚洲全国av大片| 免费在线观看黄色视频的| 性色av一级| 久久久久久久久久久久大奶| av电影中文网址| videosex国产| 99国产精品99久久久久| 19禁男女啪啪无遮挡网站| 久久久久久久精品精品| 亚洲男人天堂网一区| 国产成人啪精品午夜网站| 精品一品国产午夜福利视频| 成年女人毛片免费观看观看9 | 久久免费观看电影| 热99国产精品久久久久久7| 狠狠婷婷综合久久久久久88av| 亚洲成人手机| 国产三级黄色录像| 一级毛片精品| 91精品三级在线观看| 人人妻人人添人人爽欧美一区卜| 1024香蕉在线观看| 欧美国产精品一级二级三级| 午夜福利乱码中文字幕| 法律面前人人平等表现在哪些方面 | 黑丝袜美女国产一区| 日日摸夜夜添夜夜添小说| 国产欧美亚洲国产| av在线播放精品| 日韩三级视频一区二区三区| xxxhd国产人妻xxx| www.999成人在线观看| 精品国内亚洲2022精品成人 | 操美女的视频在线观看| av在线app专区| 久久久久久亚洲精品国产蜜桃av| av在线app专区| 国产成人欧美在线观看 | 亚洲欧洲精品一区二区精品久久久| 日韩欧美一区二区三区在线观看 | 在线 av 中文字幕| 人成视频在线观看免费观看| 国产精品99久久99久久久不卡| 欧美精品一区二区大全| 亚洲一区中文字幕在线| 国产男女超爽视频在线观看| 午夜日韩欧美国产| 日韩欧美免费精品| 不卡av一区二区三区| 电影成人av| 亚洲中文字幕日韩| 久9热在线精品视频| 亚洲国产成人一精品久久久| 欧美人与性动交α欧美精品济南到| 妹子高潮喷水视频| 国产极品粉嫩免费观看在线| 99久久人妻综合| 亚洲免费av在线视频| 国产亚洲精品第一综合不卡| 久久久久国产一级毛片高清牌| 精品一品国产午夜福利视频| 国产高清videossex| 中文字幕人妻丝袜一区二区| 啦啦啦中文免费视频观看日本| 色94色欧美一区二区| 欧美 亚洲 国产 日韩一| 亚洲性夜色夜夜综合| 亚洲欧美清纯卡通| 久久人妻熟女aⅴ| 美女国产高潮福利片在线看| 亚洲精品自拍成人| 热re99久久精品国产66热6| 美女主播在线视频| 午夜免费成人在线视频| 精品久久蜜臀av无| 亚洲伊人色综图| 精品欧美一区二区三区在线| 久久亚洲精品不卡| 91老司机精品| 黄色视频不卡| 亚洲欧美一区二区三区久久| 国产成人精品久久二区二区91| 亚洲专区国产一区二区| 亚洲美女黄色视频免费看| 久久久久久久国产电影| 十分钟在线观看高清视频www| 精品国产一区二区三区四区第35| 高清在线国产一区| 91精品三级在线观看| 99久久99久久久精品蜜桃| 国产免费福利视频在线观看| 国产精品一区二区精品视频观看| 亚洲美女黄色视频免费看| 丝瓜视频免费看黄片| 人妻久久中文字幕网| 国产精品二区激情视频| 国产成人a∨麻豆精品| 久久99一区二区三区| 亚洲三区欧美一区| 久久精品久久久久久噜噜老黄| 久久久国产欧美日韩av| 精品欧美一区二区三区在线| 国产精品久久久av美女十八| www.熟女人妻精品国产| 日本精品一区二区三区蜜桃|