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

    核設(shè)施退役過程中的輻射場重構(gòu)與拆除路徑優(yōu)化

    2017-10-13 09:21:58宋英明梁燁葉凱萱張震宇張秋楠王寧顏佳偉朱志超鄒樹梁
    核技術(shù) 2017年5期
    關(guān)鍵詞:核設(shè)施器件遺傳算法

    宋英明 梁燁 葉凱萱 張震宇 張秋楠 王寧 顏佳偉 朱志超 鄒樹梁

    ?

    核設(shè)施退役過程中的輻射場重構(gòu)與拆除路徑優(yōu)化

    宋英明1梁燁1葉凱萱1張震宇1張秋楠1王寧2顏佳偉1朱志超1鄒樹梁1

    1(南華大學(xué)核科學(xué)技術(shù)學(xué)院 衡陽 421001) 2(中國核動(dòng)力研究設(shè)計(jì)院信息中心 成都610213)

    輻射場分布和路徑選擇是核設(shè)施退役過程中影響工作人員吸收劑量的主要因素。針對(duì)核設(shè)施退役過程中多源項(xiàng)的三維輻射場,采用點(diǎn)核積分的方法進(jìn)行重構(gòu),重構(gòu)結(jié)果與蒙特卡羅程序計(jì)算結(jié)果吻合很好,驗(yàn)證了點(diǎn)核積分法的可行性。將核退役拆除路徑問題抽象為一種類旅行商問題數(shù)學(xué)模型,構(gòu)造不同拆除路徑下所受外照射劑量對(duì)應(yīng)的劑量矩陣,根據(jù)輻射防護(hù)ALARA (As Low As Reasonably Achievable)原則,利用遺傳算法進(jìn)行尋優(yōu)。對(duì)于多源項(xiàng)退役拆除實(shí)例,計(jì)算給出了最優(yōu)化拆除路徑和三維可視化顯示,并對(duì)優(yōu)化效果進(jìn)行了討論。設(shè)計(jì)開發(fā)了面向核設(shè)施退役過程的輻射場重構(gòu)與拆除路徑優(yōu)化功能軟件。

    核設(shè)施退役,輻射場重構(gòu),拆除路徑優(yōu)化,功能軟件

    核設(shè)施退役是一項(xiàng)周期長、涉及面廣、投資高的系統(tǒng)工程,基于安全性、經(jīng)濟(jì)性、合理性的考慮,在實(shí)施各項(xiàng)核設(shè)施退役工作前,必須借助計(jì)算機(jī)仿真手段對(duì)退役技術(shù)方案和關(guān)鍵步驟進(jìn)行驗(yàn)證、評(píng)估和演練,優(yōu)化拆除工藝,以降低對(duì)人員和環(huán)境的 危害。

    國際上對(duì)核設(shè)施退役仿真系統(tǒng)進(jìn)行了廣泛和深入的研究。法國原子能委員會(huì)(French Atomic Energy Commission, CEA)[1]開發(fā)了基于虛擬現(xiàn)實(shí)技術(shù)框架與點(diǎn)核積分方法的NARVEOS系統(tǒng)。該軟件支持情景設(shè)計(jì)、快速劑量計(jì)算,是適用于核輻射評(píng)估的ALARA (As Low As Reasonably Achievable)分析工具。經(jīng)濟(jì)合作與發(fā)展組織核能局(Nuclear Energy Agency, Organization for Economic Co-operation and Development, OECD/NEA)、日本原子力科學(xué)研究機(jī)構(gòu)(Japan Atomic Energy Agency, JAEA)和哈爾登(Halden)反應(yīng)堆工程中心[2]共同研發(fā)了一套人員劑量評(píng)估系統(tǒng)VRDose,可實(shí)現(xiàn)評(píng)估劑量、優(yōu)化退役計(jì)劃等功能。韓國原子能研究所和忠南大學(xué)[3]聯(lián)合對(duì)韓國KRR-1&2反應(yīng)堆退役開展工人輻射劑量計(jì)算機(jī)模擬,能實(shí)現(xiàn)基于MCNP-4C[4]輻射輸運(yùn)程序下的人的輻射劑量計(jì)算。隨著我國核能事業(yè)的高速發(fā)展,對(duì)核設(shè)施退役技術(shù)研究與工程準(zhǔn)備工作的要求更加體系化、專業(yè)化、標(biāo)準(zhǔn)化。目前,國內(nèi)已有部分單位在輻射場可視化軟件開發(fā)和核設(shè)施虛擬退役等方面開展了相關(guān)研究,取得了較好的成果[5?9],但對(duì)于核設(shè)施退役中拆除路徑最優(yōu)化方法的研究相對(duì)較少且缺乏針對(duì)性。本文驗(yàn)證并采用點(diǎn)核積分法重構(gòu)核設(shè)施退役過程中多源項(xiàng)的三維輻射場,然后利用類旅行商問題(Traveling Salesman Problem, TSP)數(shù)學(xué)模型求解核退役拆除路徑問題,構(gòu)造不同拆除路徑下所受外照射劑量矩陣,通過遺傳算法尋找最優(yōu)路徑。設(shè)計(jì)開發(fā)了軟件實(shí)現(xiàn)相應(yīng)功能并進(jìn)行了實(shí)例分析。

    1 核退役設(shè)施三維輻射場重構(gòu)

    定量和可視化描述核退役過程中工作人員所受輻射照射,需要構(gòu)建核設(shè)施的三維輻射場。由于核退役設(shè)施一般規(guī)模較大、源項(xiàng)較多,利用MCNP軟件求解核設(shè)施內(nèi)部輻射場會(huì)導(dǎo)致運(yùn)算時(shí)間很長,且難以滿足誤差要求,因此發(fā)展基于點(diǎn)核積分法的核退役設(shè)施三維輻射場求解方法。

    1.1 點(diǎn)核積分重構(gòu)三維輻射場

    點(diǎn)核積分計(jì)算核設(shè)施退役場所內(nèi)輻射場分布的方法如下[10]:

    1.2 蒙特卡羅程序校驗(yàn)點(diǎn)核積分結(jié)果

    利用MCNP5驗(yàn)證簡單源項(xiàng)情況下點(diǎn)核積分的結(jié)果,為簡化計(jì)算,將結(jié)果歸一化處理,如圖1所示。其中,圖1(a)是退役設(shè)施內(nèi)為一點(diǎn)源的結(jié)果,圖1(b)是核設(shè)施內(nèi)為一立方體源(10cm×10cm× 10cm)的結(jié)果。從圖1可以發(fā)現(xiàn),點(diǎn)源和體源情況下,點(diǎn)核積分結(jié)果都與MCNP計(jì)算結(jié)果吻合很好。這可在相同計(jì)算精度下,大幅節(jié)省最優(yōu)退役拆除路徑的搜尋時(shí)間。

    2 基于TSP的核退役拆除路徑優(yōu)化模型

    2.1 核退役拆除路徑問題的簡述

    TSP是著名的數(shù)學(xué)圖論問題,應(yīng)用背景廣泛,問題描述如下:旅行商從某城市出發(fā),遍歷各城市后返回到出發(fā)城市,尋找旅行商所走總距離最小的路線。

    在核設(shè)施退役模型中,設(shè)施內(nèi)有個(gè)器件,則拆除所有器件的路徑有!條,路徑數(shù)目隨增加而急劇上升,最優(yōu)路徑的求解難度增加。為得到最優(yōu)退役拆除路徑,將核設(shè)施退役拆除路徑問題假設(shè)成類似TSP問題,特點(diǎn)在于拆除路徑起點(diǎn)可不同、可不為閉合回路,即:

    1) 退役過程中需拆除的器件作為TSP問題中的城市;

    2) 退役過程中的優(yōu)化目標(biāo)為工作人員所受總劑量,作為TSP問題中的總路程。由于在設(shè)施中行走所需的時(shí)間遠(yuǎn)小于拆除過程所需的時(shí)間,忽略工作人員行走中所受劑量。

    2.2 不同路徑下總劑量的求解

    建立時(shí)間矩陣:

    劑量率矩陣和時(shí)間矩陣相乘,得劑量矩陣,矩陣各元素求和即為拆除人員受到的總劑量()。

    2.3 器件單元化的推廣

    實(shí)際操作中,一個(gè)器件可能分多次拆除。該情況下,設(shè)定器件單元數(shù),即一個(gè)器件平均分為個(gè)單元,各單元的活度和產(chǎn)生的劑量率為器件的1/,拆除各單元所需的時(shí)間為拆除器件總時(shí)間的1/設(shè)同一器件中各單元位置相同,則同一器件不同單元輻射效果相同。單元總數(shù)為′個(gè),拆除人員需遍歷所有單元。

    (8)

    3 遺傳算法尋優(yōu)計(jì)算與分析

    旅行商問題有多種解法,目前較為常用的算法有神經(jīng)網(wǎng)絡(luò)法和遺傳算法。遺傳算法[11]是最近幾年發(fā)展起來的全局優(yōu)化算法,借用生物遺傳學(xué)觀點(diǎn),通過選擇、遺傳、變異和免疫等作用機(jī)制,使每個(gè)個(gè)體的適應(yīng)性提高。由于全局搜索的特性,遺傳算法解決TSP問題有特定優(yōu)勢(shì)。故采用遺傳算法求解上述問題,并通過圖形用戶界面(Graphical User Interface, GUI)設(shè)計(jì)了求解本問題的可視化軟件。

    3.1 算法流程

    利用遺傳算法求解核退役拆除路徑優(yōu)化問題的步驟如下:

    1) 種群初始化

    由于核設(shè)施問題假設(shè)成多目標(biāo)旅行商問題,個(gè)體編碼方法采用實(shí)數(shù)編碼。在核設(shè)施退役中,將器件組順序?qū)崝?shù)編碼為1?(設(shè)施內(nèi)的總器件數(shù))或1?(設(shè)施內(nèi)的總單元數(shù))的隨機(jī)排列,初始化的參數(shù)有種群個(gè)數(shù)、染色體基因個(gè)數(shù)、迭代次數(shù)c、交叉概率c以及變異概率m。

    2) 建立適應(yīng)度函數(shù)

    拆除某一器件時(shí)所受劑量率由點(diǎn)核積分已知,器件組隨機(jī)全排列對(duì)應(yīng)的總劑量可求,總劑量越小,種群的適應(yīng)度函數(shù)越好。將總劑量的倒數(shù)作為適應(yīng)度函數(shù),由求min()變?yōu)榍髆ax[1/()]。

    3) 選擇操作

    用輪盤賭法和兩兩競爭的方法進(jìn)行選擇操作,采用基于適應(yīng)度比例的選擇策略,即適應(yīng)度越好的個(gè)體被選擇的概率越大,并保存所選適應(yīng)度最高的個(gè)體。

    4) 交叉操作

    利用部分映射交叉(Partial Mapped Crossover, PMX)算法進(jìn)行交叉操作,先對(duì)父代進(jìn)行常規(guī)的兩點(diǎn)交叉,再根據(jù)交叉區(qū)域內(nèi)各基因值之間的映射關(guān)系來修改交叉區(qū)域之外的各個(gè)基因值,該算法可以較多地產(chǎn)生父代遍歷路徑所沒有的部分新路徑,收斂速度較低,但其全局搜索能力較強(qiáng)。

    5) 變異操作

    對(duì)于變異操作,隨機(jī)選取個(gè)體,同時(shí)隨機(jī)選取個(gè)體的兩個(gè)基因進(jìn)行交換以實(shí)現(xiàn)變異操作。

    算法整體流程如圖2所示。

    圖2 求解核設(shè)施退役拆除路徑優(yōu)化問題的算法整體流程

    3.2 求解實(shí)例

    當(dāng)核設(shè)施內(nèi)源項(xiàng)個(gè)數(shù)為1時(shí),不同的器件單元數(shù)必然對(duì)應(yīng)同一優(yōu)化路徑。經(jīng)程序驗(yàn)證發(fā)現(xiàn),源項(xiàng)較少時(shí),不同的器件單元數(shù)趨于對(duì)應(yīng)同一優(yōu)化路徑,即單元數(shù)為1對(duì)應(yīng)的路徑。故源項(xiàng)較少情況下,器件可分多次拆除無明顯優(yōu)勢(shì);源項(xiàng)較多時(shí),分多次拆除同一器件的路徑較為復(fù)雜,為便于觀看,建立一個(gè)10cm×10cm×10cm的廠房,隨機(jī)放置30個(gè)正方體形狀的放射性器件,器件的尺寸在(1?8)×10?3m3內(nèi)隨機(jī)分布,器件的活度在(1.85?5.5)×1010內(nèi)隨機(jī)分布,假設(shè)拆除時(shí)間與拆除器件的體積成正比:t=/(其中:為體積,m3;為比例系數(shù),以下取=3為例)。

    各器件單元數(shù)=1時(shí),利用MATLAB實(shí)現(xiàn)遺傳算法求解基于旅行商問題下的最優(yōu)退役路線如圖3(a)所示,其中初始種群為100個(gè),迭代次數(shù)為1000次。由圖3(b)可知,該例中拆除人員的總劑量在迭代200次左右達(dá)到收斂,最優(yōu)路線對(duì)應(yīng)的受照總劑量為。

    單元數(shù)設(shè)為5時(shí),拆除人員可分多次對(duì)某一器件進(jìn)行拆除,最優(yōu)退役路徑如圖4(a)所示。由圖4(b)可知,該例中拆除人員的總劑量在迭代900次左右達(dá)到收斂,最優(yōu)路徑對(duì)應(yīng)受照總劑量為。

    由計(jì)算結(jié)果可知,最優(yōu)路徑可能存在于分多次拆除同一器件情況下。另外,針對(duì)核設(shè)施拆除過程中可能造成的二次污染,可通過劑量估算后,再次計(jì)算輻射場,利用功能軟件進(jìn)行二次尋優(yōu)。

    4 功能軟件開發(fā)

    利用MATLAB實(shí)現(xiàn)了基于遺傳算法求解旅行商問題下的最優(yōu)退役路線,并通過GUI設(shè)計(jì)了面向核設(shè)施退役過程的輻射場重構(gòu)與拆除路徑優(yōu)化功能軟件。

    圖5為新建廠房數(shù)據(jù)界面,在該界面可進(jìn)行三維情景設(shè)計(jì),包括放射性器件數(shù)據(jù)如器件位置、形狀、尺寸、活度、拆除所需時(shí)間及廠房相關(guān)數(shù)據(jù)并保存,廠房及任意房間的布局可視。

    圖6為最優(yōu)路徑求解界面,可在該界面更換保存的廠房數(shù)據(jù),調(diào)整器件單元數(shù)、迭代次數(shù)、種群數(shù)量等參數(shù),并顯示遺傳算法求解最優(yōu)路徑的進(jìn)程,所得結(jié)果可保存為文本文件或Excel文件。

    圖5 新建廠房數(shù)據(jù)界面

    圖6 求解最優(yōu)路徑界面

    圖7為最終的結(jié)果,給出了迭代次數(shù)與拆除人員受到總劑量的曲線,可以得到收斂后所受的最小劑量值。圖7還給出了最優(yōu)拆除路徑,最優(yōu)路徑窗口中按最優(yōu)路徑依次顯示了每次要拆除的器件以及所要拆除的百分比。

    圖7 最終結(jié)果

    圖8為最優(yōu)拆除路徑的動(dòng)態(tài)顯示,顯示速度可調(diào)。正被拆除的器件的顏色將不停變化,尺寸將不斷變小。該界面以一條移動(dòng)的短線顯示拆除人員在器件間的運(yùn)動(dòng)過程。

    圖8 最優(yōu)路徑的動(dòng)態(tài)顯示圖

    5 結(jié)語

    針對(duì)核設(shè)施退役過程中多源項(xiàng)的三維輻射場,發(fā)展了基于點(diǎn)核積分的輻射場快速重構(gòu)方法,并利用蒙特卡羅程序驗(yàn)證。核退役拆除路徑問題抽象為類TSP問題數(shù)學(xué)模型,構(gòu)造不同拆除路徑下所受外照射的劑量矩陣,根據(jù)輻射防護(hù)ALARA原則,利用遺傳算法尋優(yōu)計(jì)算。計(jì)算給出了多源項(xiàng)拆除實(shí)例最優(yōu)化拆除路徑和三維顯示,并對(duì)優(yōu)化效果進(jìn)行了討論?;谝陨戏椒ㄔO(shè)計(jì)開發(fā)了面向核設(shè)施退役過程的輻射場重構(gòu)與拆除路徑優(yōu)化功能軟件。據(jù)此可進(jìn)一步研究多旅行商問題下的核退役拆除路徑優(yōu)化問題,并考慮核設(shè)施退役過程中的個(gè)人受照劑量、集體受照劑量、經(jīng)濟(jì)代價(jià)、環(huán)境影響和公眾心理等眾多影響因素,實(shí)現(xiàn)核設(shè)施退役過程中的輻射防護(hù)最優(yōu)化,為核設(shè)施退役提供更合理、有效的決策 方案。

    1 Leservot A, Chodorge L. Chavir: a virtual site simulation environment[C]. European Nacigation Conference, Versailles, French, 2005.

    2 Rodenas J, Zarza I, Burgos M C,. Developing a virtual reality application for training nuclear power plant operators: setting up a database containing dose rates in the refuelling plant[J]. Radiation Protection Dosimetry, 2004, 111(2): 173?180. DOI: 10.1093/rpd/nch043.

    3 Kim Y H, Park W M. Use of simulation technology for prediction of radiation dose in nuclear power plant[M]. Springer Berlin Heidelberg: Computational and Information Science, 2004: 413?418.

    4 Park H S, Kim G H, Lee K W,. Development of animation and simulation module for evaluation of Worker’s dose[C]. Waste Management Symposium, Tucson, United States, 2007.

    5 戴波, 張永領(lǐng), 周斌, 等. 反應(yīng)堆退役仿真技術(shù)研究方案[J]. 核動(dòng)力工程, 2013, 34(3): 168?171. DAI Bo, ZHANG Yongling, ZHOU Bin,. Study scheme of reactor decommissioning simulation technology[J]. Nuclear Power Engineering, 2013, 34(3): 168?171.

    6 楊子輝. 大規(guī)模輻射場景三維實(shí)時(shí)仿真關(guān)鍵技術(shù)研 究[D]. 合肥: 中國科學(xué)技術(shù)大學(xué), 2016. YANG Zihui. Study on key technologies of three-dimensional real-time simulation for large-scale radiation scenario[D]. Hefei: University of Science and Technology of China, 2016.

    7 晁楠, 劉永闊, 李夢(mèng)堃, 等. 核設(shè)施退役場景的仿真實(shí)現(xiàn)技術(shù)研究[J]. 應(yīng)用科技, 2015, 42(6): 62?66. CHAO Nan, LIU Yongkuo, LI Mengkun,. A scene simulation system for decommissioning of nuclear facilities[J]. Applied Science and Technology, 2015, 42(6): 62?66.

    8 謝小龍, 胡國輝. 虛擬現(xiàn)實(shí)技術(shù)在核設(shè)施退役中的應(yīng)用[J]. 科技創(chuàng)新導(dǎo)報(bào), 2013, 260(8): 93?94. XIE Xiaolong, HU Guohui. Application of virtual reality technology in the decommissioning nuclear facilities[J]. Science and Technology Innovation Herald, 2013, 260(8): 93?94.

    9 吳宜燦, 何桃, 胡麗琴, 等. 核與輻射安全仿真系統(tǒng)SuperMC/RVIS2.3研發(fā)與應(yīng)用[J]. 原子能科學(xué)技術(shù), 2015, 49(增1): 7?15. DOI: 10.7538/yzk.2015.49.S0.0007. WU Yican, HE Tao, HU Liqin,. Development of virtual reality-based simulation system for nuclear and radiation safety[J]. Atomic Energy Science and Technology, 2015, 49(Suppl 1): 7?15. DOI: 10.7538/yzk. 2015.49.S0.0007.

    10 席福坤, 方邦城. 計(jì)算g射線點(diǎn)核積分的蒙特卡羅方 法[J]. 核動(dòng)力工程, 1987, 8(2): 66?74. XI Fukun, FANG Bangcheng. A Monte Carlo method of the point kernel integration for gamma rays calculation[J]. Nuclear Power Engineering, 1987, 8(2): 66?74.

    11 Holland J H. Adaptation in natural and artificial systems[M]. Cambridge, Massachusetts: MIT Press, 2015.

    Radiation field reconstruction and dismantling route optimization during nuclear decommissioning

    SONG Yingming1LIANG Ye1YE Kaixuan1ZHANG Zhenyu1ZHANG Qiunan1WANG Ning2YAN Jiawei1ZHU Zhichao1ZOU Shuliang1

    1(School of Nuclear Science and Technology, University of South China, Hengyang 421001, China) 2(Center of Information, Nuclear Power Institute of China, Chengdu 610213, China)

    Background: Decommissioning nuclear facilities with less cost is a key part of nuclear reactor development, which mainly depends on the dismantling route based on the distribution of radiation field. Purpose: Thisstudy attempts to reconstruct the radiation field of decommissioning nuclear facilities and find out the optimal dismantling route. Methods: The method of point-kernel integration is adopted to reconstruct the three-dimensional radiation field of multi-source terms during the decommissioning of nuclear facilities and the results are compared with the Monte Carlo program. A functional software is designed and developed with functions of radiation field reconstruction and removal route optimization during the nuclear facility decommissioning process. The problem of optimizing dismantling route is abstracted to a mathematical model of the Traveling Salesman Problem (TSP), and the Dose Matrix under different routes are constructed. According to ALARA (As Low As Reasonably Achievable) principle of the radiation protection, the methodology is optimized by genetic algorithm. Results: The reconstructed radiation field is in good agreement with the result calculated by the Monte Carlo program. Upon the instance of decommissioning multi-source terms’ removal, the optimal route and three-dimensional visualization are calculated and obtained, and the optimal results are discussed. Conclusion: The method of point-kernel integrationand the TSP mathematical model can be used during nuclear decommissioning.

    Nuclear decommissioning, Radiation field reconstruction, Dismantling route optimization, Functional software

    TL72

    10.11889/j.0253-3219.2017.hjs.40.050502

    宋英明,男,1980年出生,2009年于清華大學(xué)獲博士學(xué)位,研究領(lǐng)域?yàn)楹嗽O(shè)施退役與輻射安全

    2016-09-20,

    2016-11-07

    Supported by Special Research Fund for the Doctoral Program of Higher Education (No.20134324120003), Key Project Fund of Hunan Provincial Department of Education (No.14A120)

    高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金(No.20134324120003)、湖南省教育廳重點(diǎn)項(xiàng)目基金(No.14A120)資助

    SONG Yingming, male, born in 1980, graduated from Tsinghua University with a doctoral degree in 2009, focusing on the field of nuclear facilities and radiation safety

    2016-09-20, accepted date: 2016-11-07

    猜你喜歡
    核設(shè)施器件遺傳算法
    核設(shè)施用爆破片ASME與GB/T 567的對(duì)比分析
    核設(shè)施退役產(chǎn)業(yè)全球經(jīng)驗(yàn)探析
    能源(2018年7期)2018-09-21 07:56:36
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    核設(shè)施退役去污技術(shù)的現(xiàn)狀及發(fā)展
    中國核電(2017年2期)2017-08-11 08:01:17
    一種基于遺傳算法的聚類分析方法在DNA序列比較中的應(yīng)用
    基于遺傳算法和LS-SVM的財(cái)務(wù)危機(jī)預(yù)測
    旋涂-蒸鍍工藝制備紅光量子點(diǎn)器件
    基于改進(jìn)的遺傳算法的模糊聚類算法
    面向高速應(yīng)用的GaN基HEMT器件
    一種加載集總器件的可調(diào)三維周期結(jié)構(gòu)
    国产又色又爽无遮挡免| 真人做人爱边吃奶动态| 肉色欧美久久久久久久蜜桃| 国产精品久久久久久精品古装| 国产视频首页在线观看| 国产成人精品在线电影| 久热这里只有精品99| 精品视频人人做人人爽| 精品福利观看| 黄色怎么调成土黄色| 日韩人妻精品一区2区三区| av有码第一页| avwww免费| 悠悠久久av| 91精品伊人久久大香线蕉| 婷婷丁香在线五月| 美女主播在线视频| 国产福利在线免费观看视频| 欧美黄色片欧美黄色片| 韩国精品一区二区三区| 在线观看免费高清a一片| e午夜精品久久久久久久| 欧美日韩成人在线一区二区| 国产欧美亚洲国产| 99久久人妻综合| 大片免费播放器 马上看| 国产一区二区在线观看av| 国产精品一区二区在线观看99| 新久久久久国产一级毛片| 一本大道久久a久久精品| 2018国产大陆天天弄谢| 国产精品久久久久久人妻精品电影 | 久久久久久久精品精品| 欧美日韩视频精品一区| 国产伦人伦偷精品视频| 国产午夜精品一二区理论片| 97人妻天天添夜夜摸| 青春草亚洲视频在线观看| 电影成人av| 五月开心婷婷网| av天堂在线播放| 麻豆乱淫一区二区| 在现免费观看毛片| 男女高潮啪啪啪动态图| 久久精品亚洲av国产电影网| 国产成人一区二区在线| 久久久久久人人人人人| 高清欧美精品videossex| 欧美国产精品一级二级三级| 国产精品偷伦视频观看了| 99国产精品免费福利视频| 乱人伦中国视频| 侵犯人妻中文字幕一二三四区| 黄色a级毛片大全视频| 精品熟女少妇八av免费久了| 国产男女超爽视频在线观看| videos熟女内射| 久久久久久免费高清国产稀缺| 精品国产国语对白av| 午夜激情久久久久久久| av福利片在线| 久久久亚洲精品成人影院| 精品久久久久久久毛片微露脸 | 别揉我奶头~嗯~啊~动态视频 | 人人妻人人澡人人爽人人夜夜| 免费看不卡的av| 国产亚洲av片在线观看秒播厂| 亚洲精品成人av观看孕妇| 精品久久久精品久久久| 亚洲av电影在线观看一区二区三区| 日韩制服骚丝袜av| 亚洲五月婷婷丁香| 国产精品久久久久久人妻精品电影 | 国产精品 国内视频| 大香蕉久久网| 亚洲,欧美精品.| 天天影视国产精品| 人人妻,人人澡人人爽秒播 | 丁香六月天网| 国产无遮挡羞羞视频在线观看| av电影中文网址| 久久99热这里只频精品6学生| 各种免费的搞黄视频| 亚洲天堂av无毛| 大香蕉久久成人网| av网站免费在线观看视频| 亚洲av电影在线进入| av视频免费观看在线观看| 日本欧美视频一区| 中文字幕亚洲精品专区| 老司机在亚洲福利影院| 精品国产国语对白av| www日本在线高清视频| 亚洲自偷自拍图片 自拍| av一本久久久久| 色综合欧美亚洲国产小说| 国产欧美日韩综合在线一区二区| 免费高清在线观看视频在线观看| 男女之事视频高清在线观看 | cao死你这个sao货| 精品国产一区二区三区久久久樱花| 亚洲成av片中文字幕在线观看| 亚洲精品久久成人aⅴ小说| 国产日韩欧美亚洲二区| 亚洲精品国产区一区二| 日本欧美国产在线视频| 男男h啪啪无遮挡| 亚洲 欧美一区二区三区| 中国国产av一级| 国产av国产精品国产| 精品少妇久久久久久888优播| 国产高清videossex| 成人三级做爰电影| 婷婷丁香在线五月| 国产成人欧美在线观看 | 欧美日韩福利视频一区二区| 欧美日韩视频精品一区| 成人手机av| 青草久久国产| 久久性视频一级片| 超色免费av| 久久久亚洲精品成人影院| 亚洲黑人精品在线| 人成视频在线观看免费观看| 99精国产麻豆久久婷婷| videosex国产| 国产91精品成人一区二区三区 | 日韩中文字幕视频在线看片| 自拍欧美九色日韩亚洲蝌蚪91| 色婷婷av一区二区三区视频| 别揉我奶头~嗯~啊~动态视频 | 天天躁夜夜躁狠狠久久av| 男女免费视频国产| 免费女性裸体啪啪无遮挡网站| 亚洲成色77777| 视频区欧美日本亚洲| 欧美激情极品国产一区二区三区| 欧美激情 高清一区二区三区| 搡老岳熟女国产| 免费少妇av软件| 国产一区有黄有色的免费视频| 精品亚洲乱码少妇综合久久| 亚洲av成人不卡在线观看播放网 | 久久精品久久精品一区二区三区| 久久热在线av| 亚洲av电影在线观看一区二区三区| 青草久久国产| 亚洲五月婷婷丁香| 日韩人妻精品一区2区三区| 日本a在线网址| 欧美日韩亚洲综合一区二区三区_| 大陆偷拍与自拍| 99国产综合亚洲精品| 国产男女超爽视频在线观看| 人体艺术视频欧美日本| 成人免费观看视频高清| 在线精品无人区一区二区三| 欧美黄色片欧美黄色片| 国产精品二区激情视频| 五月开心婷婷网| 性高湖久久久久久久久免费观看| 亚洲国产日韩一区二区| 免费在线观看完整版高清| 亚洲国产精品999| 久久精品人人爽人人爽视色| 午夜福利在线免费观看网站| 久久国产精品大桥未久av| 欧美日韩亚洲综合一区二区三区_| 久久女婷五月综合色啪小说| 国产免费现黄频在线看| 亚洲成色77777| 久久精品国产亚洲av高清一级| 免费女性裸体啪啪无遮挡网站| 国产高清不卡午夜福利| 亚洲,欧美,日韩| 一二三四社区在线视频社区8| 欧美变态另类bdsm刘玥| 熟女av电影| 深夜精品福利| 黄色视频不卡| 成人亚洲欧美一区二区av| 国产有黄有色有爽视频| 亚洲人成电影观看| 久久精品亚洲av国产电影网| 操美女的视频在线观看| e午夜精品久久久久久久| 十分钟在线观看高清视频www| 精品久久久久久电影网| 黄色片一级片一级黄色片| 18禁裸乳无遮挡动漫免费视频| 中文精品一卡2卡3卡4更新| 欧美日韩亚洲国产一区二区在线观看 | 男人爽女人下面视频在线观看| 一边摸一边抽搐一进一出视频| 国产视频一区二区在线看| 超碰97精品在线观看| 亚洲国产成人一精品久久久| 一个人免费看片子| 欧美变态另类bdsm刘玥| 亚洲精品久久成人aⅴ小说| 1024香蕉在线观看| 一二三四社区在线视频社区8| 人人妻人人澡人人爽人人夜夜| 纵有疾风起免费观看全集完整版| 最近最新中文字幕大全免费视频 | 天天躁日日躁夜夜躁夜夜| 丰满少妇做爰视频| 国产片特级美女逼逼视频| 国产日韩欧美视频二区| 飞空精品影院首页| 欧美变态另类bdsm刘玥| 黄片小视频在线播放| 性高湖久久久久久久久免费观看| 亚洲一码二码三码区别大吗| 后天国语完整版免费观看| 国产一区二区在线观看av| 亚洲欧洲日产国产| 久久综合国产亚洲精品| 精品一区在线观看国产| 自拍欧美九色日韩亚洲蝌蚪91| 免费av中文字幕在线| 日本av手机在线免费观看| 少妇的丰满在线观看| 国产欧美日韩精品亚洲av| 欧美黄色片欧美黄色片| 老司机影院毛片| 国产精品国产三级专区第一集| 精品一区二区三区四区五区乱码 | 久久久精品国产亚洲av高清涩受| 久久影院123| 国产成人啪精品午夜网站| 国产1区2区3区精品| 手机成人av网站| 亚洲欧美一区二区三区国产| 另类精品久久| 欧美成人午夜精品| 涩涩av久久男人的天堂| 久久精品国产亚洲av涩爱| 亚洲一码二码三码区别大吗| 99久久精品国产亚洲精品| 亚洲 欧美一区二区三区| 欧美亚洲日本最大视频资源| 高清视频免费观看一区二区| 国产女主播在线喷水免费视频网站| 亚洲伊人久久精品综合| 脱女人内裤的视频| 三上悠亚av全集在线观看| a级片在线免费高清观看视频| 久久精品成人免费网站| 国产男女内射视频| 亚洲欧美中文字幕日韩二区| 无遮挡黄片免费观看| 亚洲av电影在线进入| 欧美成狂野欧美在线观看| 丰满人妻熟妇乱又伦精品不卡| 99国产精品99久久久久| 久久久久久免费高清国产稀缺| 国产精品久久久人人做人人爽| 少妇裸体淫交视频免费看高清 | 在线观看免费视频网站a站| 亚洲av日韩精品久久久久久密 | 国产男人的电影天堂91| 一区二区三区乱码不卡18| 欧美日韩亚洲高清精品| 考比视频在线观看| 少妇裸体淫交视频免费看高清 | 亚洲欧美激情在线| 在线观看一区二区三区激情| 巨乳人妻的诱惑在线观看| 久久久久久久国产电影| 国产精品人妻久久久影院| 亚洲国产看品久久| 亚洲精品国产av成人精品| 久久精品成人免费网站| 亚洲久久久国产精品| av在线播放精品| 精品一区二区三卡| 日本vs欧美在线观看视频| 亚洲一码二码三码区别大吗| 一本—道久久a久久精品蜜桃钙片| 欧美日韩视频高清一区二区三区二| 国产精品熟女久久久久浪| 中文字幕制服av| kizo精华| 久久久久久久久久久久大奶| 狂野欧美激情性xxxx| 久久 成人 亚洲| 麻豆乱淫一区二区| 国产真人三级小视频在线观看| 久久女婷五月综合色啪小说| 如日韩欧美国产精品一区二区三区| 亚洲精品国产一区二区精华液| 一级毛片我不卡| 欧美亚洲日本最大视频资源| 亚洲精品久久成人aⅴ小说| 久久国产精品大桥未久av| 亚洲色图 男人天堂 中文字幕| 2018国产大陆天天弄谢| 女人爽到高潮嗷嗷叫在线视频| 国产无遮挡羞羞视频在线观看| 精品亚洲乱码少妇综合久久| 亚洲第一青青草原| 亚洲av日韩在线播放| 国产国语露脸激情在线看| 国产激情久久老熟女| 最近手机中文字幕大全| 欧美成人午夜精品| 精品视频人人做人人爽| 一级毛片黄色毛片免费观看视频| 菩萨蛮人人尽说江南好唐韦庄| 欧美精品av麻豆av| 久久久久网色| 搡老岳熟女国产| 国产精品一区二区在线观看99| 亚洲五月婷婷丁香| 欧美日韩亚洲国产一区二区在线观看 | 国产99久久九九免费精品| 国产片特级美女逼逼视频| 精品亚洲乱码少妇综合久久| 亚洲欧美色中文字幕在线| 亚洲精品美女久久av网站| 亚洲精品自拍成人| 精品欧美一区二区三区在线| 亚洲av片天天在线观看| 韩国高清视频一区二区三区| 国产午夜精品一二区理论片| 欧美xxⅹ黑人| 女人精品久久久久毛片| av国产精品久久久久影院| 婷婷成人精品国产| 国产精品麻豆人妻色哟哟久久| 性色av乱码一区二区三区2| 久久久久久亚洲精品国产蜜桃av| 脱女人内裤的视频| 亚洲精品中文字幕在线视频| 久久国产精品男人的天堂亚洲| 交换朋友夫妻互换小说| 亚洲国产欧美在线一区| 欧美性长视频在线观看| 日韩 欧美 亚洲 中文字幕| 欧美精品人与动牲交sv欧美| 色综合欧美亚洲国产小说| netflix在线观看网站| www.自偷自拍.com| 美女高潮到喷水免费观看| 欧美日韩视频精品一区| 黄色a级毛片大全视频| 麻豆乱淫一区二区| 欧美激情 高清一区二区三区| 精品免费久久久久久久清纯 | av网站免费在线观看视频| 国产97色在线日韩免费| 男男h啪啪无遮挡| 亚洲精品国产av蜜桃| 亚洲国产日韩一区二区| 亚洲av国产av综合av卡| 男女国产视频网站| 狂野欧美激情性bbbbbb| 男女无遮挡免费网站观看| 免费高清在线观看日韩| 制服诱惑二区| 一级毛片女人18水好多 | 久久久欧美国产精品| 大香蕉久久网| 欧美日韩视频精品一区| av一本久久久久| 水蜜桃什么品种好| 亚洲精品久久成人aⅴ小说| 电影成人av| 男女高潮啪啪啪动态图| 手机成人av网站| 国产成人系列免费观看| 亚洲七黄色美女视频| 日本vs欧美在线观看视频| 亚洲av片天天在线观看| 乱人伦中国视频| 亚洲国产日韩一区二区| 欧美人与性动交α欧美精品济南到| av国产精品久久久久影院| 999精品在线视频| 人妻一区二区av| 亚洲精品日韩在线中文字幕| 国产人伦9x9x在线观看| 又紧又爽又黄一区二区| 宅男免费午夜| 一区二区三区精品91| 9色porny在线观看| 国产一级毛片在线| 色94色欧美一区二区| 精品欧美一区二区三区在线| av不卡在线播放| 国产一区有黄有色的免费视频| 国产欧美亚洲国产| 你懂的网址亚洲精品在线观看| 欧美黄色淫秽网站| 国产成人精品久久二区二区免费| 婷婷色综合www| 晚上一个人看的免费电影| 黄色视频不卡| 亚洲天堂av无毛| www.熟女人妻精品国产| 最近手机中文字幕大全| 国产在线观看jvid| 欧美黄色片欧美黄色片| 汤姆久久久久久久影院中文字幕| √禁漫天堂资源中文www| 女性被躁到高潮视频| 日韩精品免费视频一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩亚洲国产一区二区在线观看 | 激情视频va一区二区三区| 国产91精品成人一区二区三区 | 亚洲午夜精品一区,二区,三区| 在线观看一区二区三区激情| 老司机影院毛片| 自拍欧美九色日韩亚洲蝌蚪91| 天天躁日日躁夜夜躁夜夜| 日本五十路高清| 亚洲七黄色美女视频| 狂野欧美激情性xxxx| 777久久人妻少妇嫩草av网站| 亚洲成人国产一区在线观看 | 男女国产视频网站| 久久精品久久久久久久性| 人妻一区二区av| 久久久久久久久久久久大奶| 又大又爽又粗| 亚洲av男天堂| 亚洲精品久久午夜乱码| 啦啦啦啦在线视频资源| 一本一本久久a久久精品综合妖精| 色精品久久人妻99蜜桃| 国产精品一国产av| 99九九在线精品视频| 一区福利在线观看| 亚洲一区二区三区欧美精品| 午夜免费鲁丝| 国产爽快片一区二区三区| 交换朋友夫妻互换小说| 成人18禁高潮啪啪吃奶动态图| 一级a爱视频在线免费观看| 考比视频在线观看| 久久九九热精品免费| 99国产精品99久久久久| 无遮挡黄片免费观看| 久久久亚洲精品成人影院| 黄片播放在线免费| 丁香六月天网| 大陆偷拍与自拍| 亚洲国产av影院在线观看| 国产精品久久久久久精品电影小说| 80岁老熟妇乱子伦牲交| 精品熟女少妇八av免费久了| 亚洲欧美一区二区三区久久| 亚洲精品在线美女| 人人妻人人澡人人爽人人夜夜| 欧美成狂野欧美在线观看| av国产精品久久久久影院| 精品亚洲乱码少妇综合久久| 免费高清在线观看日韩| 色综合欧美亚洲国产小说| 国产日韩欧美视频二区| 19禁男女啪啪无遮挡网站| 狂野欧美激情性xxxx| 激情视频va一区二区三区| 久久女婷五月综合色啪小说| 美女扒开内裤让男人捅视频| 天堂俺去俺来也www色官网| 女人久久www免费人成看片| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久久久精品电影小说| 国产在线视频一区二区| 狂野欧美激情性xxxx| 一区二区三区激情视频| 国产精品久久久av美女十八| 亚洲激情五月婷婷啪啪| 亚洲精品国产色婷婷电影| 中文字幕制服av| 日韩一本色道免费dvd| 亚洲综合色网址| av不卡在线播放| 国产在线视频一区二区| 国产精品人妻久久久影院| 国产精品国产三级国产专区5o| 亚洲中文av在线| 国产精品秋霞免费鲁丝片| 一级毛片我不卡| 午夜福利一区二区在线看| 亚洲欧美成人综合另类久久久| 9热在线视频观看99| 巨乳人妻的诱惑在线观看| 国产精品一区二区在线观看99| 18禁黄网站禁片午夜丰满| 9191精品国产免费久久| 欧美在线黄色| 免费人妻精品一区二区三区视频| 热re99久久国产66热| 久久狼人影院| 麻豆国产av国片精品| 亚洲欧美一区二区三区久久| 赤兔流量卡办理| 久久精品久久久久久噜噜老黄| 777久久人妻少妇嫩草av网站| 日韩,欧美,国产一区二区三区| 免费人妻精品一区二区三区视频| 日日摸夜夜添夜夜爱| 青春草视频在线免费观看| 国产伦人伦偷精品视频| 大香蕉久久成人网| 国产高清国产精品国产三级| 男女午夜视频在线观看| 欧美精品人与动牲交sv欧美| 亚洲欧美日韩高清在线视频 | 国产99久久九九免费精品| 少妇精品久久久久久久| 日韩人妻精品一区2区三区| www.熟女人妻精品国产| 亚洲伊人色综图| 亚洲激情五月婷婷啪啪| 日本猛色少妇xxxxx猛交久久| 飞空精品影院首页| 日韩熟女老妇一区二区性免费视频| 最新在线观看一区二区三区 | www.自偷自拍.com| 女人久久www免费人成看片| 日韩中文字幕视频在线看片| 在线观看www视频免费| 欧美日韩亚洲综合一区二区三区_| 青青草视频在线视频观看| 国产亚洲av片在线观看秒播厂| 男女免费视频国产| 久久国产精品人妻蜜桃| 女人被躁到高潮嗷嗷叫费观| 日韩熟女老妇一区二区性免费视频| 成人国语在线视频| 婷婷色av中文字幕| 久久99一区二区三区| 亚洲国产欧美在线一区| 少妇粗大呻吟视频| 亚洲黑人精品在线| 国产1区2区3区精品| 人妻 亚洲 视频| 亚洲色图 男人天堂 中文字幕| www日本在线高清视频| 美女高潮到喷水免费观看| 男女高潮啪啪啪动态图| 欧美日韩亚洲高清精品| 久久精品国产亚洲av高清一级| 日日夜夜操网爽| 91麻豆精品激情在线观看国产 | 秋霞在线观看毛片| 大陆偷拍与自拍| 欧美另类一区| 无遮挡黄片免费观看| 国产熟女欧美一区二区| 丰满迷人的少妇在线观看| 久热这里只有精品99| 一区二区av电影网| 在线观看www视频免费| 国产熟女欧美一区二区| 国产精品一区二区在线不卡| 天堂俺去俺来也www色官网| 色网站视频免费| 亚洲九九香蕉| 啦啦啦在线观看免费高清www| 超碰97精品在线观看| 亚洲第一青青草原| 国产精品人妻久久久影院| 午夜视频精品福利| 国产黄色免费在线视频| 亚洲成人手机| 久久av网站| 婷婷丁香在线五月| 伊人久久大香线蕉亚洲五| 18在线观看网站| 色婷婷久久久亚洲欧美| 视频区欧美日本亚洲| 美女扒开内裤让男人捅视频| 无遮挡黄片免费观看| 国产一区二区在线观看av| av不卡在线播放| 五月开心婷婷网| 国产亚洲欧美精品永久| 91字幕亚洲| 亚洲国产欧美一区二区综合| 亚洲自偷自拍图片 自拍| 成人免费观看视频高清| 50天的宝宝边吃奶边哭怎么回事| 蜜桃在线观看..| 亚洲激情五月婷婷啪啪| 欧美国产精品va在线观看不卡| 国产午夜精品一二区理论片| 两人在一起打扑克的视频| 精品久久久精品久久久| 亚洲人成77777在线视频| 深夜精品福利| 日韩中文字幕视频在线看片| 欧美日韩黄片免| 国产免费现黄频在线看| 自拍欧美九色日韩亚洲蝌蚪91| videos熟女内射| a级毛片在线看网站| 一区二区三区四区激情视频| 久久人人97超碰香蕉20202| 久久精品成人免费网站| 日韩免费高清中文字幕av| 宅男免费午夜| 自拍欧美九色日韩亚洲蝌蚪91| 欧美+亚洲+日韩+国产| 日本av免费视频播放| 免费日韩欧美在线观看| 国产精品av久久久久免费| 9色porny在线观看| 男人添女人高潮全过程视频| 久久精品久久久久久噜噜老黄|