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

    CAP1000核電廠堆外探測器響應(yīng)函數(shù)計算方法研究

    2016-12-25 08:53:48丁謙學(xué)夏春梅梅其良
    核科學(xué)與工程 2016年2期
    關(guān)鍵詞:響應(yīng)函數(shù)共軛堆芯

    丁謙學(xué),夏春梅,梅其良

    (上海核工程研究設(shè)計院,上海200233)

    CAP1000核電廠堆外探測器響應(yīng)函數(shù)計算方法研究

    丁謙學(xué),夏春梅,梅其良

    (上海核工程研究設(shè)計院,上海200233)

    堆外探測器響應(yīng)函數(shù)代表了堆芯活性區(qū)各組件對堆外探測器讀數(shù)的貢獻,反映了堆芯功率分布與探測器讀數(shù)的關(guān)系。本文利用二維離散縱標法(SN)程序DORT,研究其共軛輸運方法,建立CAP1000反應(yīng)堆模型,分析其堆外探測器徑向和軸向響應(yīng)函數(shù)及其特性,并與采用DORT程序正向輸運計算的結(jié)果進行比較。研究表明,共軛輸運方法可以極大簡化計算量,且計算結(jié)果與正向輸運方法結(jié)果符合較好。

    離散縱標方法;共軛輸運;堆外探測器響應(yīng)函數(shù);DORT

    在核電廠中,由中子探測器和測量儀表組 成的堆外核測儀表系統(tǒng)可以測量從反應(yīng)堆堆芯泄漏的中子注量率,向操縱員提供反應(yīng)堆裝換料、停堆、啟動和功率運行狀態(tài)下的信息。堆外中子探測器的響應(yīng)函數(shù)反映了堆芯功率分布與探測器讀數(shù)的關(guān)系,可以用來連續(xù)監(jiān)測反應(yīng)堆功率、功率水平的變化和功率分布等,為堆芯在線監(jiān)測系統(tǒng)提供輸入。

    計算探測器響應(yīng)函數(shù)的方法主要有正向輸運法和共軛輸運法。與正向輸運方法相比,共軛輸運法可以大大減少計算量,提高計算效率。國內(nèi)外相關(guān)單位多采用MC方法程序進行計算,如VVER反應(yīng)堆采用MCNP程序(蒙特卡洛粒子輸運程序)通過正向輸運針對每盒燃料組件計算其響應(yīng),該方法較容易估計統(tǒng)計誤差,但計算量偏大;CPR1000反應(yīng)堆采用MCNP程序的多群共軛輸運方法進行計算,大大提高計算效率。但考慮到探測器響應(yīng)計算屬于大尺度深穿透問題,采用MC方法進行計算時,結(jié)果的收斂需要重點考慮。

    鑒于此,本文基于美國橡樹嶺國家實驗室開發(fā)的二維離散縱標法(SN)程序DORT以及配套使用的BUGLE-96數(shù)據(jù)庫,采用共軛輸運方法,針對CAP1000反應(yīng)堆計算其堆外探測器響應(yīng)函數(shù),提高計算效率,也有效避免計算結(jié)果收斂困難的問題;與SN程序正向輸運法的計算結(jié)果進行比較,驗證了DORT程序共軛輸運方法在CAP1000反應(yīng)堆探測器響應(yīng)計算中的正確性。

    1 計算原理

    堆芯裂變產(chǎn)生的中子經(jīng)過反應(yīng)堆堆芯、堆內(nèi)構(gòu)件、反射層等,發(fā)生吸收、散射等反應(yīng),最終到達堆外探測器,與探測器靈敏區(qū)發(fā)生反應(yīng)而產(chǎn)生信號。探測器響應(yīng)函數(shù)代表了堆芯內(nèi)不同位置處單位強度的裂變中子源對探測器讀數(shù)的貢獻份額。

    因此,通過粒子輸運方法進行響應(yīng)函數(shù)計算的具體原理可概括如下。

    1.1 正向輸運法

    正向輸運法的計算原理是:在堆芯內(nèi)ri處的單位體積內(nèi)設(shè)一個單位強度的裂變中子源,通過輸運計算求得該位置處的裂變中子源對堆芯外r0處探測器的響應(yīng),對于堆芯內(nèi)不同的位置,需分別進行計算,并最終將結(jié)果進行歸一。

    典型的正向輸運穩(wěn)態(tài)中子輸運方程為:

    基于此可得到堆外探測器的讀數(shù)為:

    式中:?——正向輸運計算得到的中子注量率;

    Σd——探測器所在位置處靈敏區(qū)的中子響應(yīng)截面;

    Vd——探測器靈敏區(qū)體積。

    在離散縱標法輸運計算中,假設(shè)離散后的堆芯內(nèi)網(wǎng)格數(shù)為Np,能群數(shù)為G,則堆芯活性區(qū)內(nèi)第i個網(wǎng)格點歸一化的探測器響應(yīng)函數(shù)wi為:

    式中:Σd,g——探測器靈敏區(qū)的中子響應(yīng)截面;

    ?i(r0,Eg)——堆芯內(nèi)第i個網(wǎng)格點、單位強度的裂變中子源在探測器所在位置r0處產(chǎn)生的第g群中子注量率;

    ΔVi——堆芯內(nèi)第i個網(wǎng)格點體積。

    1.2 共軛輸運法

    根據(jù)2.1節(jié)所述,基于正向輸運法,若要得到堆芯內(nèi)Np個位置的探測器響應(yīng)函數(shù),則需進行Np次正向計算,而采用共軛輸運方法,則只進行一次計算即可,可大大提高計算效率,減少計算量。

    若L表示輸運算符,則公式(1)可表示為:

    其共軛形式為:

    式中:L+——共軛輸運算符;

    ?+——共軛中子注量率,與空間位置、角度及能量有關(guān);

    S+——共軛源。

    式中:Σd(E)——在r0處的探測器的能量響應(yīng)截面。

    根據(jù)共軛算符的性質(zhì),可得到堆外探測器的讀數(shù)為:

    式中:?+——堆芯活性區(qū)各位置的共軛中子注量率;

    S——堆芯活性區(qū)裂變中子源;

    V c——堆芯活性區(qū)體積。

    各向同性的堆芯裂變中子源S可表示為:

    式中:χ(E)——堆芯內(nèi)的裂變中子能譜;

    S(r)——堆芯內(nèi)r處的裂變中子源強。

    因此,堆芯內(nèi)r i處、單位強度的裂變中子源對r0處探測器讀數(shù)的貢獻可表示為:

    離散縱標法計算中,堆芯活性區(qū)內(nèi)第i個網(wǎng)格點的歸一化后的探測器響應(yīng)函數(shù)為:

    2 計算模型和分析流程

    2.1 反應(yīng)堆模型

    選取CAP1000反應(yīng)堆的功率量程堆外探測器分析其響應(yīng)函數(shù),這些探測器分別布置在堆腔內(nèi)22.5°、112.5°、202.5°和292.5°方向,分上下兩節(jié),其靈敏區(qū)中心分別位于堆芯活性區(qū)底部以上318.77 cm和107.95 cm??紤]對稱性,只需選取22.5°位置處的探測器進行分析。

    堆芯活性區(qū)高度為426.72 cm;布置157盒燃料組件,堆芯圍板厚度為2.54 cm;吊籃厚度為5.10 cm;中子屏蔽板中心位于0°、90°、180°和270°方位角處;RPV內(nèi)表面堆焊層內(nèi)半徑為201.9 cm,厚度為0.60 cm;RPV母材厚度為21.35 cm;RPV母材外為堆腔。

    根據(jù)以上模型,建立的二維(r,θ)和(r,z)模型如圖1和圖2所示。

    圖1 CAP1000反應(yīng)堆計算模型示意圖-(r,θ)模型Fig.1 Sketch of CAP1000 reactor calculation model-(r,θ)

    圖2 CAP1000反應(yīng)堆計算模型示意圖-(r,z)模型Fig.2 Sketch of CAP1000 reactor calculation model-(r,z)

    假設(shè)所分析的探測器為10B(n,α)電離室探測器,共軛源的能譜選擇10B(n,α)反應(yīng)截面。

    2.2 分析流程

    (1)堆外探測器徑向響應(yīng)函數(shù)計算方法

    通過DORT(r,θ)計算得到堆芯區(qū)域內(nèi)(r,θ)分布的共軛中子注量率,基于公式(10)通過幾何轉(zhuǎn)換、裂變譜加權(quán)以及考慮功率分布的空間積分可以得到堆芯中任意一個組件對探測器的響應(yīng)函數(shù),將所有組件的響應(yīng)函數(shù)歸一化,即可得到當前堆芯運行狀態(tài)下,在探測器處產(chǎn)生一次反應(yīng)時,每個組件貢獻份額的徑向分布。

    需要說明的是,反應(yīng)堆運行一段時間后,由于燃耗的加深,堆內(nèi)重核增多,裂變譜變硬。此時,若只考慮235U裂變譜,會導(dǎo)致一定的偏差,本文考慮了由于燃耗不同引起的裂變譜變化從而導(dǎo)致的對探測器響應(yīng)函數(shù)的偏差。

    組件混合裂變譜通過如下公式得到:

    式中:χg——混合裂變譜;

    χg,n——核素n的裂變譜,考慮235U、238U、239Pu、240Pu、241Pu、242Pu六種可裂變核素;

    f n——各核素的裂變份額;

    υn——核素n每次裂變產(chǎn)生的中子數(shù);

    g——能群號。

    表1 CAP1000核電廠堆外探測器徑向響應(yīng)函數(shù)-共軛輸運方法Table 1 Ex-core detector response functions for CAP1000 nuclear power plant-adjoint transport method

    表2 CAP1000核電廠堆外探測器徑向響應(yīng)函數(shù)-正向輸運方法Table 2 Ex-core detector response functions for CAP1000 nuclear power plant-forward transport method

    續(xù)表

    (2)堆外探測器軸向響應(yīng)函數(shù)計算方法

    通過DORT(r,z)計算得到堆芯區(qū)域內(nèi)(r,z)分布的共軛中子注量率,考慮到堆芯裂變產(chǎn)生的中子經(jīng)過活性區(qū)、反射層以及堆內(nèi)構(gòu)件等區(qū)域的屏蔽衰減,最后對探測器計數(shù)產(chǎn)生貢獻的高能中子主要來自堆芯外圍組件,因此,基于公式(10),取堆芯外圍組件的共軛中子注量率通過裂變譜加權(quán)以及空間積分可以得到堆芯活性區(qū)軸向分布的響應(yīng)函數(shù),對其進行歸一化處理,即可得到當前堆芯運行狀態(tài)下,在探測器處產(chǎn)生一次反應(yīng)時,堆芯活性區(qū)不同軸向位置的貢獻份額。

    (3)計算方法驗證

    為驗證本文在計算響應(yīng)函數(shù)中采用的離散縱標共軛輸運方法,尤其是在將共軛中子注量率轉(zhuǎn)換為各組件的徑向響應(yīng)函數(shù)時,所作的幾何轉(zhuǎn)換、燃耗相關(guān)混合裂變譜處理等的正確性,本文采用DORT(r,θ)的正向輸運方法分別計算各組件單獨在探測器處產(chǎn)生的10B(n,α)反應(yīng)率,并對各組件的貢獻進行歸一化處理,即可得到基于正向輸運計算的每個組件對探測器響應(yīng)貢獻份額的徑向分布,并將結(jié)果與上述2.2第(1)節(jié)中所述的共軛計算結(jié)果進行比較。

    3 計算結(jié)果和分析

    表1和表2分別給出了采用離散縱標共軛輸運方法和正向輸運方法計算得到的反應(yīng)堆內(nèi)各組件對位于堆腔內(nèi)的10B(n,α)電離室探測器的徑向響應(yīng)函數(shù)。

    由表1和表2可見:

    (1)本文采用兩種方法計算得到的徑向響應(yīng)函數(shù)符合較好。其中,對探測器響應(yīng)貢獻較大的外圍組件偏差在1%之內(nèi);反應(yīng)堆中心部位組件的偏差相對較大,從4%~10%不等,但這部分組件對探測器響應(yīng)的貢獻很小,最大僅為0.01%,故對最終結(jié)果的影響可以忽略。

    (2)對探測器響應(yīng)的貢獻大部分來自反應(yīng)堆外圍組件,反應(yīng)堆堆芯最外面一圈組件的貢獻之和達到95%,其中靠近探測器的反應(yīng)堆外圍近端組件的貢獻達到90%左右。

    (3)遠離探測器的反應(yīng)堆外圍組件對響應(yīng)函數(shù)的貢獻比中心組件大,這是因為這部分組件裂變產(chǎn)生的中子經(jīng)過反射層的散射等達到探測器處的概率較大,而相比之下,堆芯中部組件產(chǎn)生的中子經(jīng)過堆芯活性區(qū)的自屏蔽,基本不會對探測器處的中子注量率產(chǎn)生影響。

    (4)本文采用正向輸運方法計算時,針對堆芯活性區(qū)每盒組件進行計算,需進行157次計算;而共軛輸運方法只需進行一次計算,并通過后續(xù)對堆芯活性區(qū)內(nèi)共軛中子注量率的裂變譜加權(quán)、空間積分等即可得到結(jié)果。與正向輸運方法相比,共軛輸運方法可以極大提高計算效率。

    本文同樣采用離散縱標共軛輸運方法計算了反應(yīng)堆堆芯活性區(qū)各位置對位于堆腔內(nèi)的10B(n,α)電離室探測器的軸向響應(yīng)函數(shù),結(jié)果見圖3。由于探測器分上下兩部分,軸向響應(yīng)函數(shù)也分別按照上部探測器和下部探測器分別給出。

    圖3 堆外探測器上下兩節(jié)的軸向響應(yīng)函數(shù)Fig.3 Axial response functions of lower and upper ex-core detectors

    從圖3可見,上下兩節(jié)探測器的軸向響應(yīng)函數(shù)并不完全對稱,這是由于堆芯半環(huán)板的不對稱布置導(dǎo)致的。

    軸向響應(yīng)函數(shù)分布如圖3。

    4 結(jié)論和建議

    本文研究了使用離散縱標共軛輸運方法計算堆外探測器響應(yīng)函數(shù)的方法,針對CAP1000反應(yīng)堆模型,使用DORT程序,分別采用正向輸運和共軛輸運方法計算堆外探測器響應(yīng)函數(shù)。

    各組件的堆外探測器響應(yīng)函數(shù)與組件位置密切相關(guān),同時,還取決于反應(yīng)堆內(nèi)各部位的材料成分、結(jié)構(gòu)尺寸等參數(shù)。

    與正向輸運法相比,本文所研究的共軛輸運法與正向輸運法計算結(jié)果吻合良好,且可以通過一次計算得到堆芯內(nèi)各位置的堆外探測器響應(yīng)函數(shù),避免了正向輸運法中對每個組件均要進行計算而導(dǎo)致的計算量過大的問題。

    [1] R.E.Maerker,Accounting for Changing Source Distribution s in Light Water Reactor Surveillance Dosimetry Analysis,Nuclear Science and Engineering,94,291-308(1986).

    [2] 周旭華,高溫氣冷堆堆外探測器空間響應(yīng)函數(shù)的計算和特性分析,核動力工程,2011年12月.

    [3] NUREG/CR-6115,PWR and BWR Pressure Vessel Fluence Calculation Benchmark Problems and Solutions,Brookhaven National Laboratory,BNL-NUREG-52395,2001.9.

    [4] MCNP-4B,Los Alamos,RSICC Computer Code Collection,Monte Carlo N-Particle Transport Code System.

    [5] DORT,A Two-Dimensional Discrete Ordinates Transport Code,ORNL/TM-11778.

    [6] BUGLE-96.Coupled 47 neutron,20 Gamma-Ray Group Cross Section Library Derived from ENDF/B-VI for LWR Shielding and Pressure Vessel Dosimetry Applications[R].DLC-185,1996.

    Research of Ex-core Detector Response Function Calculation method for CAP1000 Nuclear Power Plant

    DING Qian-xue,XIA Chun-mei,MEI Qi-liang
    (Shanghai Nuclear Engineering Research and Design Institute,Shanghai 200233,China)

    Ex-core detector response function represents the contribution to the detector tally from each reactor core assemblies,reflects the relationship between reactor core power distribution and detector tally.In this paper,using the two-dimensional discrete coordinate method(SN)program DORT to make the research of adjoint transport method,establish the CAP1000 reactor model,calculate its ex-core detector response function and analyse its characteristic in radial and axial direction.The research shows that the adjoint calculation could increase the calculation efficiency,and its calculation results are nearly accordant with the adjoint calculation results.

    discrete coordinate method;adjoint transport;ex-core detector response function;DORT

    TL328

    A

    0258-0918(2016)01-0257-06

    2015-09-21

    丁謙學(xué)(1984—),男,河北衡水人,工程師,碩士,核電廠源項分析和屏蔽設(shè)計

    猜你喜歡
    響應(yīng)函數(shù)共軛堆芯
    一類具有Beddington-DeAngelis響應(yīng)函數(shù)的階段結(jié)構(gòu)捕食模型的穩(wěn)定性
    一個帶重啟步的改進PRP型譜共軛梯度法
    一個改進的WYL型三項共軛梯度法
    巧用共軛妙解題
    一種自適應(yīng)Dai-Liao共軛梯度法
    應(yīng)用CDAG方法進行EPR機組的嚴重事故堆芯損傷研究
    相機響應(yīng)函數(shù)定標的正則化方法
    克服動態(tài)問題影響的相機響應(yīng)函數(shù)標定
    秦嶺太白山地區(qū)樹輪寬度對氣候變化的響應(yīng)
    基于Hoogenboom基準模型的SuperMC全堆芯計算能力校驗
    核技術(shù)(2016年4期)2016-08-22 09:05:32
    很黄的视频免费| 十八禁网站免费在线| 久久久精品大字幕| 国产人妻一区二区三区在| 一边摸一边抽搐一进一小说| 欧美xxxx黑人xx丫x性爽| 欧美在线一区亚洲| 男女视频在线观看网站免费| 舔av片在线| 免费av毛片视频| 国产乱人视频| 色5月婷婷丁香| 亚洲中文日韩欧美视频| 韩国av一区二区三区四区| 91九色精品人成在线观看| 老司机深夜福利视频在线观看| 搞女人的毛片| 成人永久免费在线观看视频| 免费一级毛片在线播放高清视频| 丁香六月欧美| 午夜视频国产福利| 亚洲中文日韩欧美视频| 久久久久久久久久成人| 老熟妇乱子伦视频在线观看| 网址你懂的国产日韩在线| 国产视频一区二区在线看| 久久热精品热| 老熟妇仑乱视频hdxx| 亚洲国产精品sss在线观看| 日本 av在线| 欧美精品国产亚洲| 亚洲性夜色夜夜综合| 99国产精品一区二区蜜桃av| 亚洲五月婷婷丁香| 亚洲欧美日韩高清专用| 欧美色欧美亚洲另类二区| 国产精品永久免费网站| 亚洲欧美日韩高清在线视频| 午夜免费男女啪啪视频观看 | 久久久久久久久中文| 99精品在免费线老司机午夜| 午夜精品一区二区三区免费看| 日韩欧美国产在线观看| 啦啦啦观看免费观看视频高清| 天堂网av新在线| 成人精品一区二区免费| 美女被艹到高潮喷水动态| 麻豆一二三区av精品| 人人妻人人看人人澡| 一本精品99久久精品77| 高潮久久久久久久久久久不卡| 高潮久久久久久久久久久不卡| 成人精品一区二区免费| 看黄色毛片网站| 内射极品少妇av片p| 亚洲最大成人中文| 日韩国内少妇激情av| 变态另类成人亚洲欧美熟女| 欧美日韩综合久久久久久 | 亚洲av成人精品一区久久| 成人精品一区二区免费| 日韩国内少妇激情av| 在线观看66精品国产| 动漫黄色视频在线观看| 国产男靠女视频免费网站| 美女大奶头视频| 亚洲激情在线av| 欧美一区二区国产精品久久精品| 91在线精品国自产拍蜜月| 国产成人影院久久av| 女人十人毛片免费观看3o分钟| 18美女黄网站色大片免费观看| 日韩欧美精品免费久久 | 免费人成在线观看视频色| 夜夜躁狠狠躁天天躁| 色综合婷婷激情| 成人特级黄色片久久久久久久| 亚洲电影在线观看av| 午夜福利在线在线| www日本黄色视频网| 男女做爰动态图高潮gif福利片| 国产一区二区三区在线臀色熟女| 青草久久国产| 夜夜爽天天搞| 一级黄片播放器| 国产高清三级在线| 在线播放国产精品三级| 成熟少妇高潮喷水视频| 中文字幕高清在线视频| 韩国av一区二区三区四区| 九九在线视频观看精品| 午夜福利在线在线| 男女床上黄色一级片免费看| 成人午夜高清在线视频| 午夜激情欧美在线| 免费在线观看亚洲国产| 久久国产精品人妻蜜桃| 亚洲欧美日韩东京热| 99国产精品一区二区三区| 老熟妇乱子伦视频在线观看| 赤兔流量卡办理| 美女被艹到高潮喷水动态| 久久精品91蜜桃| 不卡一级毛片| 婷婷六月久久综合丁香| 国产av不卡久久| 久久午夜福利片| 日韩亚洲欧美综合| 国产av一区在线观看免费| 看片在线看免费视频| 午夜两性在线视频| 老熟妇仑乱视频hdxx| 色在线成人网| 亚洲av成人av| 亚洲欧美日韩卡通动漫| 亚洲人成电影免费在线| 老司机午夜福利在线观看视频| 尤物成人国产欧美一区二区三区| 久久久久久久久中文| 又紧又爽又黄一区二区| 亚洲精品一卡2卡三卡4卡5卡| 男人舔奶头视频| 美女免费视频网站| 香蕉av资源在线| 亚洲内射少妇av| 午夜老司机福利剧场| 欧美极品一区二区三区四区| 免费人成在线观看视频色| 99视频精品全部免费 在线| eeuss影院久久| 亚洲最大成人av| 亚洲国产色片| 别揉我奶头 嗯啊视频| 欧美激情国产日韩精品一区| 麻豆久久精品国产亚洲av| 欧美日本亚洲视频在线播放| 色综合亚洲欧美另类图片| 天美传媒精品一区二区| 国产美女午夜福利| 午夜福利视频1000在线观看| 国产麻豆成人av免费视频| 中文字幕免费在线视频6| 久久久久精品国产欧美久久久| 亚洲av中文字字幕乱码综合| 久久久久免费精品人妻一区二区| 亚洲一区二区三区色噜噜| 91av网一区二区| 精品久久国产蜜桃| 深爱激情五月婷婷| 一级作爱视频免费观看| 两人在一起打扑克的视频| 高潮久久久久久久久久久不卡| 亚洲精品在线美女| 美女大奶头视频| 国产精品一区二区三区四区久久| 长腿黑丝高跟| 亚洲av电影在线进入| 99riav亚洲国产免费| 狂野欧美白嫩少妇大欣赏| 尤物成人国产欧美一区二区三区| 欧美精品啪啪一区二区三区| 精品国产亚洲在线| 我的女老师完整版在线观看| 免费搜索国产男女视频| 色播亚洲综合网| 亚洲av电影不卡..在线观看| 免费av观看视频| 内射极品少妇av片p| 成年女人看的毛片在线观看| 精品一区二区免费观看| 欧美zozozo另类| 亚洲不卡免费看| 人人妻人人看人人澡| 日韩欧美国产一区二区入口| 国产亚洲欧美在线一区二区| 男人和女人高潮做爰伦理| 国产成人啪精品午夜网站| 国产在线男女| 我的老师免费观看完整版| 免费高清视频大片| 啪啪无遮挡十八禁网站| 蜜桃亚洲精品一区二区三区| 国产精品人妻久久久久久| 亚洲精品在线观看二区| 精品一区二区三区视频在线观看免费| 精品久久久久久久久亚洲 | 久久天躁狠狠躁夜夜2o2o| 久久精品国产99精品国产亚洲性色| 真人做人爱边吃奶动态| 成人毛片a级毛片在线播放| 久久人人爽人人爽人人片va | 中文字幕免费在线视频6| 老司机深夜福利视频在线观看| 亚洲国产精品sss在线观看| 国产亚洲欧美在线一区二区| 18禁裸乳无遮挡免费网站照片| 一卡2卡三卡四卡精品乱码亚洲| 老司机午夜福利在线观看视频| 国产精品亚洲美女久久久| 麻豆久久精品国产亚洲av| 搡老岳熟女国产| 制服丝袜大香蕉在线| 中文字幕精品亚洲无线码一区| 亚洲不卡免费看| 婷婷亚洲欧美| 亚洲成人精品中文字幕电影| 欧美性猛交黑人性爽| 一进一出抽搐动态| 久久久久久久亚洲中文字幕 | 观看美女的网站| 在线免费观看不下载黄p国产 | 国产高清有码在线观看视频| 久久久精品欧美日韩精品| 色精品久久人妻99蜜桃| 欧美成人免费av一区二区三区| 人妻制服诱惑在线中文字幕| av中文乱码字幕在线| 久久久久久久久久黄片| 天天躁日日操中文字幕| 99精品久久久久人妻精品| 国产精品电影一区二区三区| 国产精品亚洲一级av第二区| 在线免费观看的www视频| 午夜福利在线观看吧| www.熟女人妻精品国产| 激情在线观看视频在线高清| 亚洲成人中文字幕在线播放| 国产免费av片在线观看野外av| 嫩草影院新地址| 日本黄色片子视频| 免费人成在线观看视频色| 久久久久亚洲av毛片大全| 露出奶头的视频| 特大巨黑吊av在线直播| 婷婷精品国产亚洲av| 色综合亚洲欧美另类图片| 18禁裸乳无遮挡免费网站照片| 国产欧美日韩精品一区二区| 日韩中文字幕欧美一区二区| 国产91精品成人一区二区三区| 国产伦人伦偷精品视频| 欧美日韩亚洲国产一区二区在线观看| 波野结衣二区三区在线| 搡女人真爽免费视频火全软件 | 小说图片视频综合网站| 亚洲av美国av| 啪啪无遮挡十八禁网站| 男女那种视频在线观看| 小蜜桃在线观看免费完整版高清| 嫩草影视91久久| 噜噜噜噜噜久久久久久91| 欧美丝袜亚洲另类 | 午夜激情福利司机影院| 国内精品一区二区在线观看| avwww免费| 日韩高清综合在线| 国内精品久久久久精免费| 大型黄色视频在线免费观看| 一夜夜www| 免费在线观看成人毛片| 两性午夜刺激爽爽歪歪视频在线观看| 宅男免费午夜| 非洲黑人性xxxx精品又粗又长| 在线十欧美十亚洲十日本专区| 草草在线视频免费看| 夜夜爽天天搞| 国产欧美日韩精品一区二区| av中文乱码字幕在线| 国产三级黄色录像| 校园春色视频在线观看| 免费无遮挡裸体视频| 少妇的逼好多水| 蜜桃亚洲精品一区二区三区| 亚洲精品在线美女| 亚洲av第一区精品v没综合| 看免费av毛片| av中文乱码字幕在线| 一本久久中文字幕| 日韩欧美在线乱码| 男人舔女人下体高潮全视频| 欧美日韩中文字幕国产精品一区二区三区| 蜜桃久久精品国产亚洲av| 精品一区二区三区视频在线观看免费| 国产一区二区在线av高清观看| 亚洲欧美日韩高清专用| 一本综合久久免费| 国产精品一及| 国产精品av视频在线免费观看| 久久精品国产亚洲av香蕉五月| 精品国产三级普通话版| 嫩草影院入口| 欧美黄色淫秽网站| 在线观看免费视频日本深夜| 亚洲色图av天堂| 国产精品影院久久| 深爱激情五月婷婷| 国产综合懂色| 午夜影院日韩av| 一区二区三区免费毛片| 欧美日本视频| 久久精品久久久久久噜噜老黄 | 久久久久免费精品人妻一区二区| 精品久久久久久,| 亚洲aⅴ乱码一区二区在线播放| 窝窝影院91人妻| 欧美在线黄色| 老司机深夜福利视频在线观看| 久久久久国产精品人妻aⅴ院| 少妇丰满av| 亚洲av熟女| 久久精品国产99精品国产亚洲性色| 国产亚洲精品av在线| 久久国产乱子伦精品免费另类| 18禁在线播放成人免费| 久久久久九九精品影院| 99热这里只有是精品50| 日韩欧美精品v在线| 欧美成人一区二区免费高清观看| 免费人成在线观看视频色| 精品人妻1区二区| 久久精品国产清高在天天线| 熟女人妻精品中文字幕| 真人做人爱边吃奶动态| 赤兔流量卡办理| 婷婷六月久久综合丁香| 亚洲五月婷婷丁香| 亚洲三级黄色毛片| 日本黄大片高清| 久久精品国产亚洲av香蕉五月| 日韩av在线大香蕉| 日日摸夜夜添夜夜添小说| avwww免费| 一本久久中文字幕| 欧美+日韩+精品| 一进一出好大好爽视频| 三级毛片av免费| 黄色配什么色好看| 宅男免费午夜| 日韩有码中文字幕| 国内精品一区二区在线观看| 国产一区二区在线观看日韩| 国产日本99.免费观看| 欧美一区二区精品小视频在线| 一夜夜www| 天堂av国产一区二区熟女人妻| 国产真实伦视频高清在线观看 | 亚洲不卡免费看| 一进一出抽搐gif免费好疼| 亚洲最大成人av| 很黄的视频免费| 国内精品美女久久久久久| 欧美+亚洲+日韩+国产| 亚洲av日韩精品久久久久久密| 男人的好看免费观看在线视频| 搡老熟女国产l中国老女人| 久久国产乱子免费精品| 久久久久久久精品吃奶| 久久久久精品国产欧美久久久| 亚洲黑人精品在线| 日韩欧美在线乱码| 精品乱码久久久久久99久播| 国产亚洲欧美98| 亚洲电影在线观看av| a级毛片免费高清观看在线播放| 亚洲美女搞黄在线观看 | 久久九九热精品免费| 亚洲欧美日韩东京热| 99热只有精品国产| 又爽又黄无遮挡网站| 久久精品91蜜桃| 最近视频中文字幕2019在线8| 国产又黄又爽又无遮挡在线| 国产老妇女一区| 一个人免费在线观看电影| www.熟女人妻精品国产| 国产精品久久久久久亚洲av鲁大| 日日干狠狠操夜夜爽| 99久久99久久久精品蜜桃| 18禁黄网站禁片免费观看直播| 国产美女午夜福利| 国产精品一区二区三区四区免费观看 | 小说图片视频综合网站| 在线观看美女被高潮喷水网站 | 欧美日韩瑟瑟在线播放| 99在线人妻在线中文字幕| 亚洲男人的天堂狠狠| 久久久久久九九精品二区国产| 午夜福利视频1000在线观看| 国产精品乱码一区二三区的特点| 在线免费观看的www视频| 成熟少妇高潮喷水视频| 亚洲乱码一区二区免费版| 亚洲av免费高清在线观看| 老熟妇乱子伦视频在线观看| 国产精品精品国产色婷婷| 日本三级黄在线观看| 神马国产精品三级电影在线观看| 欧美成人一区二区免费高清观看| 国产精品不卡视频一区二区 | 中文字幕久久专区| 久久久久久久久大av| 性欧美人与动物交配| 又粗又爽又猛毛片免费看| 嫩草影院入口| 亚洲一区二区三区不卡视频| 亚洲自偷自拍三级| 校园春色视频在线观看| 亚洲精品在线观看二区| 日韩亚洲欧美综合| 国产精品一区二区三区四区免费观看 | 淫妇啪啪啪对白视频| 一区福利在线观看| 麻豆av噜噜一区二区三区| 亚洲专区中文字幕在线| 午夜福利视频1000在线观看| 99久久九九国产精品国产免费| 在线观看美女被高潮喷水网站 | 国内揄拍国产精品人妻在线| 日本撒尿小便嘘嘘汇集6| 国产午夜福利久久久久久| 两个人视频免费观看高清| av专区在线播放| 美女高潮喷水抽搐中文字幕| 欧美xxxx性猛交bbbb| 少妇高潮的动态图| 99热只有精品国产| 久久精品国产亚洲av天美| 久久久久国内视频| 国产免费男女视频| 波多野结衣高清无吗| 禁无遮挡网站| 免费无遮挡裸体视频| 欧美性感艳星| 99热这里只有是精品50| 草草在线视频免费看| 亚洲avbb在线观看| 国产男靠女视频免费网站| 久久久久九九精品影院| 久久久久国产精品人妻aⅴ院| 亚洲最大成人av| 国产色婷婷99| 亚洲av.av天堂| 热99re8久久精品国产| 999久久久精品免费观看国产| 69av精品久久久久久| 赤兔流量卡办理| 亚洲国产欧美人成| 最近最新中文字幕大全电影3| 3wmmmm亚洲av在线观看| 久99久视频精品免费| 又紧又爽又黄一区二区| 国产黄色小视频在线观看| 大型黄色视频在线免费观看| 国产毛片a区久久久久| 丰满乱子伦码专区| 成人美女网站在线观看视频| 国产免费av片在线观看野外av| 一卡2卡三卡四卡精品乱码亚洲| 国产野战对白在线观看| 国产成人aa在线观看| 在线观看免费视频日本深夜| 免费人成视频x8x8入口观看| 757午夜福利合集在线观看| 少妇裸体淫交视频免费看高清| 久久久久久久亚洲中文字幕 | 亚洲国产精品久久男人天堂| 亚洲国产精品sss在线观看| 亚洲欧美激情综合另类| 91字幕亚洲| 嫁个100分男人电影在线观看| 看十八女毛片水多多多| 日韩av在线大香蕉| 人妻制服诱惑在线中文字幕| 精品免费久久久久久久清纯| 搞女人的毛片| 国产精品自产拍在线观看55亚洲| 看黄色毛片网站| 热99在线观看视频| 国产精品久久电影中文字幕| 亚洲欧美日韩无卡精品| 少妇人妻精品综合一区二区 | 色在线成人网| 亚洲成a人片在线一区二区| av中文乱码字幕在线| 很黄的视频免费| 神马国产精品三级电影在线观看| 精品人妻熟女av久视频| 1000部很黄的大片| 真人一进一出gif抽搐免费| 国产精品久久久久久久电影| 久久亚洲真实| 欧美激情久久久久久爽电影| 九九热线精品视视频播放| 国产伦精品一区二区三区四那| 少妇的逼好多水| 成人一区二区视频在线观看| 国产国拍精品亚洲av在线观看| 亚洲av成人精品一区久久| 高清日韩中文字幕在线| 又爽又黄无遮挡网站| 久久人人精品亚洲av| 桃色一区二区三区在线观看| 亚洲真实伦在线观看| 精品熟女少妇八av免费久了| 免费在线观看影片大全网站| 久久久久国内视频| 女同久久另类99精品国产91| 嫩草影院入口| 三级男女做爰猛烈吃奶摸视频| 国模一区二区三区四区视频| 麻豆一二三区av精品| 久久久久久大精品| 99热这里只有是精品50| 久久精品国产清高在天天线| 午夜久久久久精精品| 日韩欧美三级三区| 久久久精品大字幕| 欧美黑人欧美精品刺激| 亚洲国产精品久久男人天堂| 草草在线视频免费看| 成人欧美大片| 亚洲乱码一区二区免费版| 欧美丝袜亚洲另类 | 亚洲真实伦在线观看| 51国产日韩欧美| 亚洲乱码一区二区免费版| 亚洲第一欧美日韩一区二区三区| 久久人人精品亚洲av| 免费看日本二区| 成人三级黄色视频| 色综合欧美亚洲国产小说| 少妇被粗大猛烈的视频| 欧美日韩亚洲国产一区二区在线观看| 日韩中字成人| 欧洲精品卡2卡3卡4卡5卡区| 中文字幕人成人乱码亚洲影| 内地一区二区视频在线| 草草在线视频免费看| 99久久成人亚洲精品观看| 搡女人真爽免费视频火全软件 | 精品人妻视频免费看| 99久久精品国产亚洲精品| 亚洲欧美日韩无卡精品| 如何舔出高潮| 热99在线观看视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 啪啪无遮挡十八禁网站| 网址你懂的国产日韩在线| 十八禁人妻一区二区| 国产老妇女一区| 亚洲最大成人手机在线| ponron亚洲| aaaaa片日本免费| 此物有八面人人有两片| 18禁黄网站禁片午夜丰满| 亚洲成人精品中文字幕电影| 国产视频内射| 亚洲美女黄片视频| 色综合欧美亚洲国产小说| 精品国产三级普通话版| 十八禁国产超污无遮挡网站| 亚洲熟妇熟女久久| 国产成人a区在线观看| 国产欧美日韩精品亚洲av| 日韩欧美在线乱码| 免费av不卡在线播放| 精品久久国产蜜桃| 欧美绝顶高潮抽搐喷水| 欧美黑人巨大hd| 免费看光身美女| 免费人成视频x8x8入口观看| 国产精品98久久久久久宅男小说| 国产成人影院久久av| 成人高潮视频无遮挡免费网站| www.色视频.com| 最新在线观看一区二区三区| 在线观看一区二区三区| 欧美在线一区亚洲| 小蜜桃在线观看免费完整版高清| 欧美日韩乱码在线| 男女视频在线观看网站免费| 91狼人影院| 精品免费久久久久久久清纯| 免费无遮挡裸体视频| 一区二区三区免费毛片| 91午夜精品亚洲一区二区三区 | 国产精品久久久久久人妻精品电影| 亚洲国产精品合色在线| 国产高清有码在线观看视频| 国产黄片美女视频| 18禁裸乳无遮挡免费网站照片| 国产亚洲精品综合一区在线观看| 欧美最新免费一区二区三区 | 热99re8久久精品国产| 大型黄色视频在线免费观看| 又爽又黄无遮挡网站| 国产免费av片在线观看野外av| 亚洲三级黄色毛片| 噜噜噜噜噜久久久久久91| 亚洲中文日韩欧美视频| 一区二区三区四区激情视频 | 亚洲无线在线观看| 成人国产综合亚洲| 亚洲五月婷婷丁香| 亚洲精品亚洲一区二区| 国内精品久久久久久久电影| 少妇人妻一区二区三区视频| 国产精品久久电影中文字幕| 亚州av有码| 成人精品一区二区免费| 嫩草影院新地址| 在线天堂最新版资源| 能在线免费观看的黄片| 亚洲欧美日韩无卡精品| 又爽又黄a免费视频|