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

    基于遺傳算法的可燃毒物設(shè)計(jì)優(yōu)化方法研究

    2021-08-02 03:02:54劉仕倡李滿倉(cāng)周冰燕王連杰陳義學(xué)
    原子能科學(xué)技術(shù) 2021年8期
    關(guān)鍵詞:燃耗毒物堆芯

    肖 鵬,王 健,劉仕倡,*,李滿倉(cāng),周冰燕,王連杰,陳義學(xué)

    (1.中國(guó)核動(dòng)力研究設(shè)計(jì)院 核反應(yīng)堆系統(tǒng)設(shè)計(jì)技術(shù)重點(diǎn)實(shí)驗(yàn)室,四川 成都 610213; 2.華北電力大學(xué) 核科學(xué)與工程學(xué)院,北京 102206)

    核反應(yīng)堆燃料組件設(shè)計(jì)是提高反應(yīng)堆燃料利用率、堆芯功率展平及反應(yīng)性控制的關(guān)鍵。其中,可燃毒物用于降低初始堆芯剩余反應(yīng)性和冷卻劑的硼濃度,保證反應(yīng)堆具有負(fù)的慢化劑溫度系數(shù)。可燃毒物組件與控制棒組件不同,為堆芯不動(dòng)部件??扇级疚锝M件是按物理設(shè)計(jì)分散布置于堆芯適當(dāng)位置的燃料組件中,有利于展平堆芯中子注量率的分布,可提高功率密度[1-2]。通過(guò)可燃毒物材料、含量、毒物棒數(shù)量等參數(shù),在滿足燃料和冷卻劑等反應(yīng)性系數(shù)在壽期內(nèi)為負(fù)、功率峰因子低于限值等約束條件下,實(shí)現(xiàn)降低初始反應(yīng)性、反應(yīng)性變化較平緩、堆芯壽期盡量長(zhǎng)、毒物在壽期內(nèi)消耗速率平緩、壽期末毒物殘留少、反應(yīng)性懲罰較小等優(yōu)化目標(biāo)[3]。對(duì)于核反應(yīng)堆燃料優(yōu)化設(shè)計(jì),由于是燃料棒、燃料組件、堆芯的多尺度問(wèn)題,同時(shí)包含燃料棒、毒物棒、控制棒等多個(gè)對(duì)象,而且涉及到不同溫度、密度等各種工況,甚至還是隨堆芯燃耗深度不斷變化的動(dòng)態(tài)問(wèn)題,因此核反應(yīng)堆燃料優(yōu)化設(shè)計(jì)具有很大的挑戰(zhàn)性。

    目前,常規(guī)的反應(yīng)堆燃料設(shè)計(jì)過(guò)程通常是設(shè)計(jì)者根據(jù)經(jīng)驗(yàn)給出一個(gè)初步方案,再針對(duì)此方案進(jìn)行堆芯物理計(jì)算分析以驗(yàn)證其設(shè)計(jì)的性能,如不滿足設(shè)計(jì)目標(biāo),則調(diào)整設(shè)計(jì)參數(shù)給出新方案并再次進(jìn)行計(jì)算,如此反復(fù)[4-5]。該方法能否得到優(yōu)化方案,很大程度上取決于設(shè)計(jì)者的經(jīng)驗(yàn),對(duì)于設(shè)計(jì)經(jīng)驗(yàn)缺乏人員往往耗時(shí)巨大且很難得到滿意解,即使是經(jīng)驗(yàn)豐富的設(shè)計(jì)者,進(jìn)行優(yōu)化迭代設(shè)計(jì)也非常耗時(shí)耗力。

    核反應(yīng)堆燃料優(yōu)化設(shè)計(jì)是一多輸入多目標(biāo)問(wèn)題。遺傳算法作為一種有效的隨機(jī)搜索算法,具有智能性、并行性、通用性等特點(diǎn)[6]。它在許多領(lǐng)域得到了應(yīng)用,取得了良好的效果。本文將遺傳算法和蒙特卡羅粒子輸運(yùn)方法相結(jié)合,研究遺傳算法在壓水堆燃料組件可燃毒物多目標(biāo)優(yōu)化設(shè)計(jì)中的理論模型和實(shí)現(xiàn)方法。

    1 壓水堆燃料組件中的可燃毒物設(shè)計(jì)

    本文研究的可燃毒物類型為可燃毒物及燃料混合在一起的整體型可燃毒物。優(yōu)化對(duì)象為13×13的壓水堆組件,其中包含3種棒:燃料棒、毒物棒和中心導(dǎo)管。組件結(jié)構(gòu)如圖1所示。

    圖1 壓水堆組件結(jié)構(gòu)模型 Fig.1 PWR assembly model

    1.1 可燃毒物優(yōu)化變量

    優(yōu)化變量包括:1) 毒物材料類型,包括6種氧化物,即Gd2O3、Er2O3、Sm2O3、Eu2O3、Dy2O3、HfO2(表1);2) 毒物含量,采用的范圍是毒物質(zhì)量占含可燃毒物燃料棒整體質(zhì)量(毒物+UO2)的1%~20%共20種質(zhì)量百分比情況;3) 排布類型(圖2);4) 軸向分層,在三維燃耗計(jì)算中將軸向分為5層,排列方式1表示整根棒全為可燃毒物,2表示兩端各有1層Zr層,中間3層為可燃毒物,3表示兩端4層為Zr層,只有中間1層為可燃毒物(圖3)。

    圖2 含可燃毒物燃料組件布置Fig.2 Fuel assembly layouts with burnable poison

    圖3 軸向分區(qū)布置Fig.3 Axial division

    表1 候選可燃毒物的密度和宏觀吸收截面Table 1 Theoretical density and macroscopic neutron absorption cross section of candidate burnable poison

    1.2 可燃毒物設(shè)計(jì)目標(biāo)

    從中子學(xué)角度,可燃毒物設(shè)計(jì)應(yīng)遵循3個(gè)原則:1) 壽期初,可燃毒物引入足夠大的負(fù)反應(yīng)性;2) 壽期內(nèi),可燃毒物緩慢釋放反應(yīng)性,控制功率分布,降低功率峰因子;3) 壽期末,可燃毒物的反應(yīng)性懲罰盡可能小。

    2 多目標(biāo)優(yōu)化遺傳算法

    遺傳算法是一種基于自然選擇的搜索算法[7],它通過(guò)繁殖將最適合的個(gè)體組合起來(lái),從而找到一個(gè)多元問(wèn)題的近似最優(yōu)解[8],并將潛在的新思想和變異性引入突變?nèi)后w中。該算法具有很強(qiáng)的并行性,因此本文將遺傳算法應(yīng)用于毒物選型優(yōu)化這種搜索空間非常復(fù)雜的多目標(biāo)優(yōu)化問(wèn)題[9]。

    針對(duì)多目標(biāo)優(yōu)化問(wèn)題[10],目前采用最多的是帕累托最優(yōu)解方法。帕累托前沿是以設(shè)計(jì)目標(biāo)為軸的優(yōu)化圖,由在圖上創(chuàng)建的線或平面顯示有沖突設(shè)計(jì)目標(biāo)之間的優(yōu)化權(quán)衡。在帕累托前沿,對(duì)于1個(gè)或多個(gè)目標(biāo),每個(gè)單獨(dú)的最優(yōu)值可能較任何其他最優(yōu)值都好,但不是所有的目標(biāo)都是最優(yōu)的。帕累托最優(yōu)的成員也被稱為“非占優(yōu)”[11]。

    本文的多目標(biāo)優(yōu)化采用NSGA-Ⅲ算法。NSGA-Ⅲ隨機(jī)生成含有N個(gè)個(gè)體的初始種群,算法進(jìn)行迭代直至終止條件滿足。在當(dāng)前種群的基礎(chǔ)上,通過(guò)隨機(jī)選擇產(chǎn)生子代種群。為從種群中選擇最好的N個(gè)解進(jìn)入下一代,首先利用基于帕累托前沿的非支配排序?qū)⑵浞譃槿舾刹煌姆侵鋵?。然后算法?gòu)建1個(gè)新的種群,構(gòu)建方法是逐次將各非支配層的解加入到新種群中,直至新種群大小等于N,或首次大于N,丟棄第N層之后的解。對(duì)于最后一層非支配層解的選取,NSGA-Ⅱ用的是基于擁擠距離的方法,而NSGA-Ⅲ用的是基于參考點(diǎn)的方法。擁擠距離度量并不適合求解3個(gè)及更多目標(biāo)的多目標(biāo)優(yōu)化問(wèn)題。因此NSGA-Ⅲ采用了新的選擇機(jī)制,該機(jī)制會(huì)通過(guò)所提供的參考點(diǎn),對(duì)新種群中的個(gè)體進(jìn)行更加系統(tǒng)的分析,以選擇最后的非支配層中部分解進(jìn)入新種群。因此本文選擇NSGA-Ⅲ作為多目標(biāo)優(yōu)化算法[12]。

    同時(shí),為提高計(jì)算效率,本文采用并行遺傳算法,將種群內(nèi)的個(gè)體分配到不同的計(jì)算核心。常規(guī)遺傳算法和并行遺傳算法的基本流程示于圖4。圖4中,COREn表示第n個(gè)計(jì)算核心。不同的計(jì)算核心處理不同的單個(gè)子種群,種群間互相獨(dú)立進(jìn)行進(jìn)化,種群間進(jìn)行個(gè)體遷移和種群競(jìng)爭(zhēng)。因此并行遺傳算法可顯著提高計(jì)算效率。

    圖4 并行遺傳算法基本流程圖Fig.4 Flow chart of parallel genetic algorithm

    3 可燃毒物優(yōu)化數(shù)學(xué)模型

    3.1 決策變量與目標(biāo)函數(shù)

    二維組件燃耗計(jì)算的輸入?yún)?shù)共3個(gè),故定義決策變量維數(shù)為3。第1個(gè)輸入變量為可燃毒物組件的排列方式(共7種),表示含可燃毒物燃料組件布置,具體排布見(jiàn)1.1節(jié);第2個(gè)輸入變量為毒物類型,包括6種不同材料;第3個(gè)輸入變量是可燃毒物的含量,為1%~20%共20種質(zhì)量占比。三維燃耗計(jì)算增加了1個(gè)決策變量,第4個(gè)決策變量是毒物棒的軸向分層,具體裝載情況見(jiàn)1.1節(jié)。4個(gè)變量以隨機(jī)數(shù)的形式拼接在1條染色體基因鏈上。由于各自的范圍不同,4個(gè)變量在基因鏈上的長(zhǎng)度亦不同。因?yàn)榛蜴準(zhǔn)怯啥M(jìn)制編碼的,所以4個(gè)輸入變量在染色體基因鏈上的長(zhǎng)度分別為3、3、7、2字節(jié),基因鏈總長(zhǎng)為15字節(jié)。

    燃耗計(jì)算目標(biāo)函數(shù)共5個(gè),分別表示為F1、F2、F3、F4、F5。F1為初始有效增殖因數(shù)(keff),F(xiàn)2為各燃耗步最大功率峰因子(ppf),F(xiàn)3為可燃毒物在燃耗中期的含量,F(xiàn)4為可燃毒物在燃耗末期的含量,F(xiàn)5為各燃耗步中有效增殖因數(shù)的最小值(表2)。其中,毒物含量為毒物氧化物中吸收截面最大的同位素剩余百分比。總?cè)己臅r(shí)長(zhǎng)為1 000 d,燃耗中期含量設(shè)為第500 d的毒物含量與初始毒物含量的比值;末期含量為第1 000 d毒物含量與初始毒物含量的比值。為降低初始反應(yīng)性,展平徑向功率分布,減少壽期末毒物殘留,分別設(shè)置優(yōu)化目標(biāo)為F1、F2和F4最小化。為使可燃毒物緩慢釋放反應(yīng)性,達(dá)到最大的循環(huán)長(zhǎng)度,分別設(shè)置優(yōu)化目標(biāo)為F3和F5最大化。

    表2 燃耗計(jì)算決策變量和目標(biāo)函數(shù)Table 2 Decision variable and objective function of burnup calculation

    3.2 約束條件

    組件燃耗計(jì)算中選用約束條件的方法對(duì)目標(biāo)結(jié)果加以限制。在燃耗計(jì)算中,設(shè)置了4個(gè)約束條件。首先使初始有效增殖因數(shù)大于1,即F1>1。二維燃耗計(jì)算采用最大功率峰因子F2<1.4的約束條件。而在三維燃耗計(jì)算中,因三維燃耗的加入使得功率不均勻程度增大,因此設(shè)置約束條件使各燃耗步最大功率峰因子F2<3.0。此外,燃耗中期要保證一定量的可燃毒物剩余,燃耗末期剩余的可燃毒物不能過(guò)多,故增加2個(gè)約束使燃耗中期剩余毒物的含量大于1%和燃耗末期剩余毒物的含量小于10%,即F3>0.01和F4<0.1。

    4 計(jì)算結(jié)果

    本文的遺傳算法采用Python的Geatpy遺傳算法庫(kù)[13]。Geatpy是一個(gè)高性能實(shí)用型進(jìn)化算法工具箱,提供許多已實(shí)現(xiàn)的進(jìn)化算法中各項(xiàng)重要操作的庫(kù)函數(shù),并提供高度模塊化、耦合度低的面向?qū)ο蟮倪M(jìn)化算法框架,利用“定義問(wèn)題類+調(diào)用算法模板”的模式進(jìn)行進(jìn)化優(yōu)化,可用于求解單目標(biāo)優(yōu)化、多目標(biāo)優(yōu)化、復(fù)雜約束優(yōu)化、組合優(yōu)化、混合編碼進(jìn)化優(yōu)化等。組件臨界計(jì)算及燃耗計(jì)算采用清華大學(xué)REAL團(tuán)隊(duì)研發(fā)的堆用蒙特卡羅程序RMC[14],RMC具備內(nèi)耦合的燃耗計(jì)算功能[15-17]。

    4.1 二維組件燃耗計(jì)算

    遺傳算法的種群數(shù)為50,迭代次數(shù)為50。RMC計(jì)算的源迭代采用的粒子數(shù)為2 000/代,非活躍代為5代,總代數(shù)為15,keff的統(tǒng)計(jì)標(biāo)準(zhǔn)差約為0.5%。比功率為30 W/gHM,燃料棒中235U的富集度為3%。不同棒數(shù)具體排列方式見(jiàn)1.1節(jié)。燃耗計(jì)算中采用的燃耗步數(shù)為6,步長(zhǎng)分別為100、100、150、150、250、250 d,共1 000 d。最大平均燃耗深度為30 MW·d/kgHM。不考慮帶軸向分層,毒物棒為整根全為燃料和毒物混合的單一材料。優(yōu)化方案及結(jié)果列于表3,最終的優(yōu)化方案中所用的材料均為Er2O3,因此Er2O3較其他材料更優(yōu)。該結(jié)果與通過(guò)人工搜索優(yōu)化得到的結(jié)論一致。二維優(yōu)化方案的keff隨燃耗時(shí)間的變化如圖5所示。

    圖5 二維組件各燃耗步keff隨燃耗時(shí)間的變化Fig.5 Change of keff with burnup time of two-dimensional assembly

    表3 二維組件燃耗計(jì)算優(yōu)化方案及結(jié)果Table 3 Optimization scheme and result of two-dimensional assembly burnup calculation

    4.2 三維組件燃耗計(jì)算

    考慮加入軸向分層后的燃耗計(jì)算,遺傳算法的種群數(shù)和迭代次數(shù)、蒙特卡羅計(jì)算的粒子數(shù)和代數(shù)均與二維算例相同。keff的統(tǒng)計(jì)標(biāo)準(zhǔn)差約為0.4%。

    優(yōu)化后得到40個(gè)非支配解,其中Gd2O3有13個(gè)、Er2O3有21個(gè)、Sm2O3有2個(gè)、Eu2O3有4個(gè)(表4)。與二維燃耗的只有1種材料Er2O3相比,通過(guò)軸向分層使得候選材料種類增加。各材料的keff隨燃耗時(shí)間的變化如圖6所示。

    圖6 三維組件各燃耗步keff隨時(shí)間的變化Fig.6 Change of keff with burnup time of three-dimensional assembly

    表4 三維組件燃耗計(jì)算優(yōu)化方案及結(jié)果Table 4 Optimization scheme and result of three-dimensional assembly burnup calculation

    三維燃耗ppf隨不同軸向分層的分布如圖7所示??煽闯觯謱臃绞?(即兩端各有1層Zr層、中間3層為可燃毒物)的燃耗ppf較分層方式1(即整根棒全為可燃毒物)的小??梢?jiàn)改變分層方式,以部分毒物代替整體全為可燃毒物的布置方式可減小ppf值。

    圖7 三維燃耗ppf隨軸向分層分布的變化Fig.7 Distribution of ppf in three-dimensional burnup with axial division

    5 結(jié)論

    本文將多目標(biāo)并行遺傳算法應(yīng)用于壓水堆組件毒物選型優(yōu)化,以反應(yīng)性控制、功率分布和不同時(shí)期燃耗剩余等為目標(biāo),對(duì)可燃毒物材料類型、含可燃毒物燃料棒排列方式、毒物含量、軸向分層等決策變量進(jìn)行了優(yōu)化。針對(duì)二維和三維燃耗計(jì)算,分別篩選了13和40種優(yōu)化方案。計(jì)算結(jié)果表明:Er2O3用作毒物的綜合效果最好;Gd2O3、Eu2O3和Sm2O3的應(yīng)用需結(jié)合堆芯方案開展進(jìn)一步研究;HfO2和Dy2O3不適合用作可燃毒物。該結(jié)果與通過(guò)人工搜索優(yōu)化得到的結(jié)論基本一致。同時(shí),三維軸向分層可為優(yōu)化提供更多可選的材料種類方案,以部分毒物的分層布置方式可減小ppf值。

    通過(guò)本文計(jì)算,初步驗(yàn)證了智能優(yōu)化算法在進(jìn)行毒物選型優(yōu)化時(shí),具有一定的協(xié)助工程師進(jìn)行優(yōu)化設(shè)計(jì)方案篩選的潛力。

    猜你喜歡
    燃耗毒物堆芯
    快跑!有毒物 絕密毒藥報(bào)告
    毒物——水銀
    應(yīng)用CDAG方法進(jìn)行EPR機(jī)組的嚴(yán)重事故堆芯損傷研究
    基于Hoogenboom基準(zhǔn)模型的SuperMC全堆芯計(jì)算能力校驗(yàn)
    核技術(shù)(2016年4期)2016-08-22 09:05:32
    基于切比雪夫有理逼近方法的蒙特卡羅燃耗計(jì)算研究與驗(yàn)證
    核技術(shù)(2016年4期)2016-08-22 09:05:28
    AMDIS在法醫(yī)毒物分析中的應(yīng)用
    IFBA/WABA 可燃毒物元件的燃耗特性分析
    壓水堆堆芯中應(yīng)用可燃毒物的兩個(gè)重要實(shí)驗(yàn)
    低價(jià)值控制棒中子吸收體材料燃耗相關(guān)數(shù)據(jù)的制作及驗(yàn)證研究
    基于SOP規(guī)程的大亞灣堆芯冷卻監(jiān)測(cè)系統(tǒng)改造
    99国产综合亚洲精品| 在线观看三级黄色| av在线播放精品| 不卡av一区二区三区| 自线自在国产av| 久久亚洲国产成人精品v| 久久久国产精品麻豆| 只有这里有精品99| 在线观看国产h片| 亚洲一区中文字幕在线| 青春草视频在线免费观看| 免费高清在线观看日韩| 亚洲av男天堂| 久久精品久久精品一区二区三区| 悠悠久久av| 久久亚洲国产成人精品v| 天天躁夜夜躁狠狠久久av| av片东京热男人的天堂| 午夜福利影视在线免费观看| 亚洲美女搞黄在线观看| 国产精品.久久久| 久久性视频一级片| xxxhd国产人妻xxx| 国产欧美日韩一区二区三区在线| 国产精品.久久久| 波多野结衣一区麻豆| 制服人妻中文乱码| 国产精品久久久久久精品古装| av视频免费观看在线观看| 久久人妻熟女aⅴ| 亚洲国产av新网站| 精品国产国语对白av| 亚洲精品国产av蜜桃| 欧美黄色片欧美黄色片| 亚洲美女黄色视频免费看| 久久ye,这里只有精品| a 毛片基地| 在线观看一区二区三区激情| 成人毛片60女人毛片免费| 日韩大片免费观看网站| av女优亚洲男人天堂| tube8黄色片| 岛国毛片在线播放| 日本一区二区免费在线视频| 亚洲国产欧美网| 中文字幕亚洲精品专区| 久久女婷五月综合色啪小说| 午夜老司机福利片| 日韩免费高清中文字幕av| 伊人久久国产一区二区| 天天添夜夜摸| 777米奇影视久久| 亚洲欧美清纯卡通| 亚洲国产看品久久| 男女边摸边吃奶| 国产精品三级大全| 亚洲精品在线美女| 宅男免费午夜| 99精品久久久久人妻精品| 天美传媒精品一区二区| 亚洲美女搞黄在线观看| 男女之事视频高清在线观看 | 欧美精品人与动牲交sv欧美| 黑人猛操日本美女一级片| a级毛片黄视频| 久久国产精品男人的天堂亚洲| 免费观看av网站的网址| 夫妻午夜视频| 久久精品久久久久久久性| 亚洲第一青青草原| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美乱码精品一区二区三区| 搡老岳熟女国产| 天堂中文最新版在线下载| 18禁裸乳无遮挡动漫免费视频| 精品国产国语对白av| 精品亚洲成a人片在线观看| 叶爱在线成人免费视频播放| 2018国产大陆天天弄谢| 国产精品国产三级国产专区5o| 国产精品国产av在线观看| 亚洲av国产av综合av卡| 日本爱情动作片www.在线观看| 婷婷色av中文字幕| 国产成人精品在线电影| 可以免费在线观看a视频的电影网站 | 少妇 在线观看| 午夜福利一区二区在线看| 一级,二级,三级黄色视频| 又粗又硬又长又爽又黄的视频| 成年人午夜在线观看视频| 欧美97在线视频| 天天操日日干夜夜撸| 夫妻午夜视频| 女人被躁到高潮嗷嗷叫费观| 人体艺术视频欧美日本| 欧美中文综合在线视频| 热re99久久精品国产66热6| 午夜福利一区二区在线看| 亚洲av电影在线进入| 国产在线免费精品| 在线免费观看不下载黄p国产| 狂野欧美激情性bbbbbb| 老熟女久久久| 免费看av在线观看网站| 亚洲av在线观看美女高潮| 18禁观看日本| 国产欧美亚洲国产| 1024视频免费在线观看| 亚洲av成人精品一二三区| 各种免费的搞黄视频| 美女扒开内裤让男人捅视频| 男女边吃奶边做爰视频| 最近中文字幕高清免费大全6| 久久久久网色| 亚洲av综合色区一区| 久久国产精品男人的天堂亚洲| 美女午夜性视频免费| 欧美激情极品国产一区二区三区| 男女边摸边吃奶| 午夜激情久久久久久久| 亚洲精品国产av成人精品| 下体分泌物呈黄色| 91成人精品电影| 老司机影院毛片| 爱豆传媒免费全集在线观看| 黄网站色视频无遮挡免费观看| 亚洲,欧美,日韩| 观看av在线不卡| 亚洲,一卡二卡三卡| 国产深夜福利视频在线观看| a级毛片在线看网站| 成人影院久久| 中文字幕制服av| 捣出白浆h1v1| 免费不卡黄色视频| 国产精品一二三区在线看| 日韩,欧美,国产一区二区三区| 性少妇av在线| 麻豆精品久久久久久蜜桃| 成年人免费黄色播放视频| 亚洲欧美成人综合另类久久久| av国产久精品久网站免费入址| 波多野结衣av一区二区av| 久久久久久久久久久久大奶| 亚洲色图综合在线观看| 欧美97在线视频| 色精品久久人妻99蜜桃| 高清不卡的av网站| 新久久久久国产一级毛片| 亚洲av欧美aⅴ国产| 亚洲av电影在线观看一区二区三区| 欧美成人精品欧美一级黄| 亚洲成人av在线免费| 大码成人一级视频| 欧美成人精品欧美一级黄| 精品国产乱码久久久久久男人| √禁漫天堂资源中文www| 国产精品久久久久久精品古装| 国产男女超爽视频在线观看| 国产视频首页在线观看| 91成人精品电影| 男女边摸边吃奶| 一区二区三区四区激情视频| 国产淫语在线视频| 久久精品久久久久久噜噜老黄| 人妻人人澡人人爽人人| 久久毛片免费看一区二区三区| 欧美日韩一级在线毛片| 下体分泌物呈黄色| 夫妻性生交免费视频一级片| 久久久久人妻精品一区果冻| 免费黄频网站在线观看国产| 欧美亚洲 丝袜 人妻 在线| 久久 成人 亚洲| 在线观看一区二区三区激情| 看非洲黑人一级黄片| 我要看黄色一级片免费的| 蜜桃在线观看..| 国产欧美日韩综合在线一区二区| 免费在线观看视频国产中文字幕亚洲 | 国产激情久久老熟女| 亚洲精品国产区一区二| 亚洲一码二码三码区别大吗| 毛片一级片免费看久久久久| a级片在线免费高清观看视频| 日韩精品有码人妻一区| 久久人人爽人人片av| 成人免费观看视频高清| 一边亲一边摸免费视频| av线在线观看网站| 亚洲精品,欧美精品| 天天躁日日躁夜夜躁夜夜| 久久久精品免费免费高清| 最黄视频免费看| 亚洲欧美一区二区三区久久| 亚洲精品一二三| 亚洲国产av影院在线观看| 水蜜桃什么品种好| av.在线天堂| 青春草亚洲视频在线观看| 亚洲伊人色综图| 天天躁夜夜躁狠狠久久av| 91精品国产国语对白视频| avwww免费| 日本色播在线视频| 一二三四在线观看免费中文在| 美女高潮到喷水免费观看| 欧美日韩亚洲国产一区二区在线观看 | 热99国产精品久久久久久7| 亚洲国产精品999| 精品久久久久久电影网| 亚洲av综合色区一区| 亚洲第一区二区三区不卡| 国产精品女同一区二区软件| 男女无遮挡免费网站观看| 欧美日韩一级在线毛片| 青春草国产在线视频| 丝袜美腿诱惑在线| 日日摸夜夜添夜夜爱| 精品一品国产午夜福利视频| 国产欧美日韩综合在线一区二区| 一区二区三区激情视频| 欧美在线一区亚洲| 一区二区三区乱码不卡18| 亚洲精品一二三| 日韩大片免费观看网站| 九九爱精品视频在线观看| av网站免费在线观看视频| 卡戴珊不雅视频在线播放| 一区二区日韩欧美中文字幕| 十分钟在线观看高清视频www| 一区二区三区激情视频| 亚洲av日韩在线播放| 看免费av毛片| 久久天堂一区二区三区四区| 成人三级做爰电影| 国产日韩欧美亚洲二区| 成人国产麻豆网| 侵犯人妻中文字幕一二三四区| 久久国产精品大桥未久av| 亚洲国产精品一区三区| 两个人免费观看高清视频| 日本wwww免费看| 国产乱来视频区| 十八禁网站网址无遮挡| 日本vs欧美在线观看视频| 欧美人与性动交α欧美软件| 熟女av电影| 精品一区二区免费观看| 国产精品嫩草影院av在线观看| 黄片小视频在线播放| 69精品国产乱码久久久| 热re99久久国产66热| 亚洲第一av免费看| 黄频高清免费视频| 中文字幕高清在线视频| 国产在线视频一区二区| 久久热在线av| 免费看不卡的av| 18禁裸乳无遮挡动漫免费视频| 国产精品久久久久久久久免| 亚洲av国产av综合av卡| 九色亚洲精品在线播放| 色网站视频免费| 久久久久久久精品精品| 国产毛片在线视频| 欧美中文综合在线视频| 91精品国产国语对白视频| 亚洲欧美成人精品一区二区| svipshipincom国产片| 亚洲成人免费av在线播放| 男人操女人黄网站| 国产亚洲av高清不卡| 国产高清不卡午夜福利| 少妇猛男粗大的猛烈进出视频| 免费看不卡的av| 精品午夜福利在线看| 国产深夜福利视频在线观看| 操出白浆在线播放| 欧美成人午夜精品| 亚洲精品日韩在线中文字幕| www日本在线高清视频| 久久鲁丝午夜福利片| 男女国产视频网站| 日韩精品免费视频一区二区三区| 亚洲综合色网址| 亚洲精品一区蜜桃| 人人妻人人澡人人看| 麻豆乱淫一区二区| 伊人亚洲综合成人网| 99久久综合免费| 国产探花极品一区二区| 黄片播放在线免费| 久久av网站| 多毛熟女@视频| 波多野结衣av一区二区av| 亚洲欧美成人综合另类久久久| 激情视频va一区二区三区| 亚洲国产欧美一区二区综合| 午夜激情av网站| 免费女性裸体啪啪无遮挡网站| 91精品伊人久久大香线蕉| 99国产精品免费福利视频| 亚洲欧洲精品一区二区精品久久久 | 性高湖久久久久久久久免费观看| 日本爱情动作片www.在线观看| 少妇被粗大的猛进出69影院| 在线观看三级黄色| 欧美在线黄色| 亚洲av福利一区| 日日啪夜夜爽| bbb黄色大片| 国产成人欧美| 69精品国产乱码久久久| 久久久久久久大尺度免费视频| 超色免费av| 一二三四在线观看免费中文在| 亚洲一码二码三码区别大吗| 国产片特级美女逼逼视频| 美女脱内裤让男人舔精品视频| 一二三四中文在线观看免费高清| www.av在线官网国产| 欧美久久黑人一区二区| 国产成人一区二区在线| 9色porny在线观看| 在线精品无人区一区二区三| 观看美女的网站| 亚洲视频免费观看视频| 成人18禁高潮啪啪吃奶动态图| 国产精品香港三级国产av潘金莲 | 最近手机中文字幕大全| 国产有黄有色有爽视频| 我要看黄色一级片免费的| 国产av精品麻豆| 毛片一级片免费看久久久久| 亚洲人成77777在线视频| 高清不卡的av网站| 色94色欧美一区二区| 国产99久久九九免费精品| 欧美黄色片欧美黄色片| 人人澡人人妻人| 熟女av电影| 制服人妻中文乱码| 一本一本久久a久久精品综合妖精| 男女高潮啪啪啪动态图| 宅男免费午夜| 综合色丁香网| 久久久久精品人妻al黑| 亚洲 欧美一区二区三区| 深夜精品福利| 国产99久久九九免费精品| 深夜精品福利| 亚洲精品久久午夜乱码| 99国产精品免费福利视频| 欧美日韩视频高清一区二区三区二| 久热这里只有精品99| 日本猛色少妇xxxxx猛交久久| 亚洲欧美色中文字幕在线| 一级毛片 在线播放| 男女午夜视频在线观看| 精品一区二区免费观看| 国产亚洲精品第一综合不卡| 久久久久视频综合| 国产色婷婷99| √禁漫天堂资源中文www| 久久精品aⅴ一区二区三区四区| 一级片'在线观看视频| 美女午夜性视频免费| 色综合欧美亚洲国产小说| 天天影视国产精品| 国产亚洲欧美精品永久| 国产av一区二区精品久久| 男女边吃奶边做爰视频| 久久亚洲国产成人精品v| 丰满迷人的少妇在线观看| 免费在线观看完整版高清| 亚洲情色 制服丝袜| 深夜精品福利| 看免费av毛片| 亚洲成色77777| 午夜91福利影院| 最近中文字幕2019免费版| 欧美日韩国产mv在线观看视频| 最近中文字幕2019免费版| 一区二区三区乱码不卡18| 中文字幕人妻丝袜制服| 国产老妇伦熟女老妇高清| 永久免费av网站大全| 97人妻天天添夜夜摸| 亚洲精品视频女| 麻豆精品久久久久久蜜桃| 国产精品无大码| 欧美 日韩 精品 国产| 日韩伦理黄色片| 成年av动漫网址| 久久这里只有精品19| videos熟女内射| 亚洲成国产人片在线观看| 日本欧美视频一区| 国产精品欧美亚洲77777| 国产精品国产av在线观看| 国产一区二区 视频在线| 亚洲精品美女久久久久99蜜臀 | 飞空精品影院首页| 精品久久久久久电影网| 日韩视频在线欧美| 欧美少妇被猛烈插入视频| 久久毛片免费看一区二区三区| 丰满迷人的少妇在线观看| 捣出白浆h1v1| 亚洲国产最新在线播放| 午夜福利网站1000一区二区三区| 一级黄片播放器| 欧美精品一区二区免费开放| 韩国精品一区二区三区| xxx大片免费视频| 久久精品熟女亚洲av麻豆精品| 国产精品三级大全| 深夜精品福利| 大片免费播放器 马上看| 国产熟女欧美一区二区| 午夜激情av网站| 一区二区三区精品91| 国产成人欧美| 七月丁香在线播放| 男人添女人高潮全过程视频| 亚洲精品国产色婷婷电影| 午夜福利,免费看| 亚洲,一卡二卡三卡| 亚洲人成77777在线视频| 欧美人与性动交α欧美软件| 色精品久久人妻99蜜桃| 悠悠久久av| 日韩制服骚丝袜av| 欧美在线黄色| 国产一区二区 视频在线| 黄色视频不卡| 婷婷色综合www| 交换朋友夫妻互换小说| 久久久精品区二区三区| 日本wwww免费看| 巨乳人妻的诱惑在线观看| 国产片内射在线| 国产成人一区二区在线| 亚洲婷婷狠狠爱综合网| 欧美在线一区亚洲| 婷婷色av中文字幕| av片东京热男人的天堂| 久热这里只有精品99| 91老司机精品| 国产男女超爽视频在线观看| 亚洲欧美中文字幕日韩二区| 在线观看免费日韩欧美大片| av线在线观看网站| 免费久久久久久久精品成人欧美视频| 高清视频免费观看一区二区| 久久精品国产综合久久久| 亚洲欧美一区二区三区黑人| 麻豆精品久久久久久蜜桃| 少妇被粗大猛烈的视频| 一级爰片在线观看| 在线观看国产h片| 性高湖久久久久久久久免费观看| 亚洲一区二区三区欧美精品| 在线免费观看不下载黄p国产| 国产熟女午夜一区二区三区| 日本91视频免费播放| 国产激情久久老熟女| 两个人免费观看高清视频| 人人妻人人添人人爽欧美一区卜| 最近中文字幕2019免费版| 黄色视频在线播放观看不卡| 欧美精品一区二区免费开放| 男女午夜视频在线观看| 亚洲精品国产色婷婷电影| 久久鲁丝午夜福利片| 国语对白做爰xxxⅹ性视频网站| 国产女主播在线喷水免费视频网站| 国产精品一二三区在线看| 久久久精品国产亚洲av高清涩受| avwww免费| 人人妻人人添人人爽欧美一区卜| 国产av码专区亚洲av| 国产精品一区二区精品视频观看| 国产 一区精品| 国产野战对白在线观看| 欧美精品av麻豆av| videosex国产| 国产精品99久久99久久久不卡 | 涩涩av久久男人的天堂| 亚洲精品国产色婷婷电影| 五月开心婷婷网| www日本在线高清视频| 青草久久国产| 久久久国产欧美日韩av| 搡老乐熟女国产| 人人妻人人澡人人爽人人夜夜| 免费日韩欧美在线观看| 两性夫妻黄色片| 国产亚洲av片在线观看秒播厂| 一边亲一边摸免费视频| 啦啦啦 在线观看视频| 国产精品久久久久成人av| 最新的欧美精品一区二区| 日韩制服丝袜自拍偷拍| 国产精品 欧美亚洲| 99re6热这里在线精品视频| av电影中文网址| 男女边摸边吃奶| 自线自在国产av| 日本猛色少妇xxxxx猛交久久| 日韩,欧美,国产一区二区三区| av不卡在线播放| 国产黄频视频在线观看| 国产欧美日韩一区二区三区在线| 人人妻人人澡人人看| 久久人人97超碰香蕉20202| www日本在线高清视频| 五月开心婷婷网| 高清av免费在线| 超碰97精品在线观看| 婷婷色麻豆天堂久久| 国产精品国产三级国产专区5o| 亚洲人成电影观看| 水蜜桃什么品种好| 婷婷色麻豆天堂久久| 精品国产超薄肉色丝袜足j| 国产 一区精品| 亚洲精品国产av蜜桃| 美女主播在线视频| 国产精品蜜桃在线观看| 国产精品国产三级专区第一集| 亚洲欧美激情在线| 欧美亚洲 丝袜 人妻 在线| 性色av一级| 日本av免费视频播放| 18禁动态无遮挡网站| 成人国产麻豆网| av在线app专区| 亚洲精品乱久久久久久| 青春草视频在线免费观看| 国产精品 欧美亚洲| 亚洲欧美成人精品一区二区| 国产精品女同一区二区软件| 成年人免费黄色播放视频| 美女福利国产在线| 欧美日韩综合久久久久久| 不卡av一区二区三区| 亚洲av欧美aⅴ国产| 国产亚洲av高清不卡| 国产精品无大码| 午夜精品国产一区二区电影| 日本wwww免费看| 国产成人系列免费观看| 高清av免费在线| 亚洲精品一二三| 久久精品久久精品一区二区三区| 亚洲成色77777| 我要看黄色一级片免费的| 久久免费观看电影| 另类精品久久| 日本wwww免费看| 亚洲视频免费观看视频| 日韩免费高清中文字幕av| 国产 一区精品| 成年动漫av网址| 卡戴珊不雅视频在线播放| 麻豆av在线久日| 人人妻,人人澡人人爽秒播 | 欧美中文综合在线视频| 午夜福利在线免费观看网站| 天堂俺去俺来也www色官网| 成年女人毛片免费观看观看9 | 成人毛片60女人毛片免费| 999久久久国产精品视频| 男女国产视频网站| 你懂的网址亚洲精品在线观看| 亚洲精品中文字幕在线视频| 国产男人的电影天堂91| 亚洲av日韩精品久久久久久密 | 亚洲美女搞黄在线观看| 久久午夜综合久久蜜桃| 制服丝袜香蕉在线| 久久久久久人妻| 中文字幕最新亚洲高清| 成年美女黄网站色视频大全免费| 纯流量卡能插随身wifi吗| 黑人猛操日本美女一级片| 精品国产露脸久久av麻豆| 婷婷成人精品国产| 夫妻午夜视频| 精品国产乱码久久久久久男人| 亚洲国产欧美日韩在线播放| 欧美国产精品一级二级三级| 国产男女超爽视频在线观看| 超碰97精品在线观看| 男女边吃奶边做爰视频| 欧美日韩亚洲综合一区二区三区_| 亚洲国产欧美日韩在线播放| 久久久精品94久久精品| 叶爱在线成人免费视频播放| 如何舔出高潮| 亚洲欧洲精品一区二区精品久久久 | 国产麻豆69| 欧美人与善性xxx| 欧美日韩av久久| 一级毛片 在线播放| 欧美日本中文国产一区发布| 亚洲欧美日韩另类电影网站| 国产一区亚洲一区在线观看| 久久精品人人爽人人爽视色| 啦啦啦啦在线视频资源| 国产精品欧美亚洲77777| 欧美日韩一区二区视频在线观看视频在线|