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

    快堆MOX燃料氧與金屬原子比控制原理

    2014-08-08 06:37:42屈哲昊尹邦躍楊啟法梁雪元
    原子能科學技術 2014年7期
    關鍵詞:芯塊點缺陷平衡常數(shù)

    屈哲昊,尹邦躍,楊啟法,梁雪元

    (1.中國原子能科學研究院 反應堆工程研究設計所,北京 102413;

    2.核工業(yè)標準化研究所,北京 100091)

    氧與金屬原子比(O/M比,余同)指MOX燃料即(U1-yPuy)O2±x混合氧化物燃料晶格中氧原子與金屬原子個數(shù)的比值,即O/M比=2±x。O/M比是MOX燃料最重要的性能指標之一,對MOX燃料的熔點、晶格常數(shù)、理論密度、熱物性、擴散系數(shù)、氧在燃料中的化學位即氧勢以及堆內輻照性能均有很大影響。而PuO2含量、Pu分布均勻性、燒結溫度、燒結氣氛氧分壓、粘接劑種類和添加量、批次質量等因素均對O/M比有影響。Toci等[1]指出,富钚顆粒的氧勢一般要高于钚均勻分布的MOX固溶體的氧勢。中國實驗快堆MOX燃料的PuO2質量分數(shù)為25%,O/M比為1.96~1.99,如何從理論計算和工藝控制方面達到所要求的O/M比是一個重要的研究課題。

    國外從20世紀50年代開始研究核燃料的氧勢和O/M比,現(xiàn)已得到數(shù)千個數(shù)據(jù),但這些數(shù)據(jù)通常是在有限的“氧勢-溫度-成分”范圍內重復進行,由于實驗測量條件不同,數(shù)據(jù)分散性較大。各國也提出了很多氧勢理論模型,但由于采用不同的離子平衡假設和不同的熱力學數(shù)據(jù),計算結果與實驗數(shù)據(jù)存在較大差異,至今沒有一個普遍公認的氧勢模型[2-7]。根據(jù)經驗,在MOX燃料芯塊生產過程中,一般通過控制燒結氣氛的氧分壓即氣體氧勢來調整固體芯塊的氧勢和O/M比。顯然,了解MOX燃料的O/M比控制原理將有利于優(yōu)化MOX燃料芯塊的燒結工藝參數(shù),本文將進行MOX燃料的O/M比控制原理的研究。

    1 氧勢模型

    氧化物陶瓷芯塊在保溫燒結過程中,固體芯塊內的氧勢與燒結氣氛的氧勢存在平衡,根據(jù)熱力學定律,氧化學位為:

    (1)

    (2)

    MOX芯塊的O/M比是固體內的化學位即氧勢、溫度和Pu含量的函數(shù),為建立O/M比與MOX燃料芯塊燒結工藝(燒結溫度、氧分壓和PuO2含量)之間的關系,需通過熱力學參數(shù)或離子反應平衡常數(shù)精確計算氧勢,并建立O/M比與氧勢之間的關系式。國外早在20世紀50年代就已開始進行核燃料的氧勢研究,從晶體結構缺陷、熱力學行為角度出發(fā),發(fā)展了燃料的多種氧勢熱力學理論和數(shù)學模型[8]。公認的氧勢模型主要有Rand-Markin模型、Blackburn模型、相平衡模型、點缺陷模型。這些理論模型采用了一些假設和有差異的熱力學數(shù)據(jù),測量數(shù)據(jù)由于原理或技術的原因存在誤差,使氧勢模型的應用有限。

    1) Rand-Markin價態(tài)模型

    Rand-Markin模型是基于實驗結果的經驗模型,它把MOX視為UO2+x與PuO2-x形成的理想固溶體(U1-yPuy)O2±x。該模型假設:當MOX芯塊為超化學計量時,即(U1-yPuy)-O2+x的O/M比大于2.00時,氧勢為鈾平均價態(tài)的函數(shù);當其為亞化學計量時,即(U1-yPuy)-O2-x的O/M比小于2.00時,氧勢為钚平均價態(tài)的函數(shù),且不受溫度和Pu/U比的影響[9]。

    對于亞化學計量(U1-yPuy)O2-x而言,U離子的價態(tài)VU是+4,但Pu離子的平均價態(tài)VPu不是整數(shù)。

    VPu=4-2x/y

    (3)

    對于超化學計量(U1-yPuy)O2+x而言,Pu離子的價態(tài)為+4,但U離子的平均價態(tài)不是整數(shù):

    VU=4+2x/(1-y)

    (4)

    2) Blackburn離子平衡模型

    Blackburn模型[12]可認為是半經驗性模型,它認為MOX固溶體中存在確定的U、Pu離子價態(tài),各價態(tài)離子相互平衡,構成一系列平衡方程。MOX燃料通常在強還原性氣氛中進行燒結,離子平衡反應要考慮U2+-U4+、U4+-U6+、Pu3+-Pu4+、Pu2+-Pu3+等反應體系,其對應的反應平衡常數(shù)分別為:

    (5)

    (6)

    (7)

    (8)

    式中,[U2+]、[U4+]、[U6+]、[Pu2+]、[Pu3+]、[Pu4+]、[O2-]分別為U2+、U4+、U6+、Pu2+、Pu3+、Pu4+、O2-的濃度,余同。

    對于確定成分與O/M比的(U1-yPuy)O2±x燃料,x和y是已知量,有:

    [U2+]+[U4+]+[U6+]=1-y

    (9)

    [Pu2+]+[Pu3+]+[Pu4+]=y

    (10)

    2[U2+]+4[U4+]+6[U6+]+2[Pu2+]+

    3[Pu3+]+4[Pu4+]=2±x=2[O2-]

    (11)

    當反應平衡常數(shù)已知,聯(lián)立解上述7個公式即可求得pO2;或當pO2和y已知,同樣可根據(jù)式(5)~(8)求解x,從而得到混合氧化物燃料的O/M比。圖1為利用Blackburn模型計算得到的MOX燃料氧勢與O/M比的關系[13]。

    Blackburn模型計算過程較簡單,方法易掌握,計算量少,精確度受反應平衡常數(shù)的影響較大。Blackburn發(fā)布模型時,利用相圖中UO2/U4O9、UO2-x/U(l)平衡數(shù)據(jù)分別計算得到了KU24、KU46。也可利用已知的pO2、x、y等數(shù)據(jù),擬合得到平衡常數(shù)。受Stan等[11,14]的啟發(fā),若結合量子力學理論與分子動力學計算技術進行反應平衡常數(shù)的理論估算,再利用實驗數(shù)據(jù)進行修正,可得到更準確的反應平衡常數(shù),從而提高Blackburn模型的計算準確度。

    圖1 Blackburn模型計算得到的MOX燃料的氧勢與O/M比的關系[13]

    3) 相平衡模型

    相平衡模型假定MOX燃料中U4O9與UO2、Pu2O3與PuO2之間存在相互平衡,且U4O9、Pu2O3這兩種含非化學計量缺陷的化合物可對應地溶解于無缺陷的化學計量比化合物UO2、PuO2中。

    8UO2+O2(g)=2U4O9

    2Pu2O3+O2(g)=4PuO2

    形成(U1-yPuy)Ox固溶體的標準吉布斯自由能為:

    (12)

    相平衡模型是Hoch等[15]在錒系氧化物中采用的模型,Besmann等[16-17]和Krishnaiah[18]將其應用于MOX燃料。目前很少采用該方法計算MOX燃料的氧勢,僅Varamban等[19]利用此方法預測MOX燃料的氧勢。

    4) 點缺陷模型

    由于非化學計量化學物中必然存在間隙離子、空位、自由電子、電子空隙等結構缺陷,它們之間也會構成平衡,因此點缺陷模型認為可通過計算缺陷平衡反應來得到固溶體的氧勢。

    MOX燃料可能存在以下缺陷反應(使用Kroger-Vink符號表示缺陷):

    假設上述6個反應平衡常數(shù)已知,則可通過計算MOX固溶體的氧空位與氧間隙原子濃度計算O/M比:

    (13)

    利用點缺陷模型求解MOX燃料氧勢的計算結果與實驗結果符合較好,此方法被Kato等[20-21]廣泛應用于MOX燃料和嬗變燃料的研究中。但目前尚無文獻報道如何通過理論計算得到平衡常數(shù)。Kato等[22]利用實際測量的氧勢擬合計算得到U0.8Pu0.2O2±x的各平衡常數(shù)如下:

    (14)

    (15)

    (16)

    (17)

    (18)

    (19)

    Blackburn模型和點缺陷模型是目前應用最多的兩種理論模型,并已積累了大量實驗數(shù)據(jù)。利用國外現(xiàn)有離子反應平衡常數(shù)和熱力學數(shù)據(jù),采用Blackburn模型可精確研究MOX燃料的氧勢,但目前這些數(shù)據(jù)較少,今后應開展更為精確的點缺陷模型計算研究。

    2 MOX燃料O/M比控制工藝計算

    2.1 氧分壓計算

    為使MOX燃料的O/M比達到指標設計要求,要根據(jù)氧勢模型計算出氧勢的控制范圍,再根據(jù)氧勢精確計算結果計算在某燒結溫度下的氧分壓,最后根據(jù)氧分壓控制氣體露點。

    假設快堆MOX燃料芯塊的設計成分分別為(U0.7Pu0.3)O2-x、(U0.75Pu0.25)O2-x、(U0.8Pu0.2)O2-x,分別利用Blackburn模型和點缺陷模型計算得到不同O/M比對應的氧分壓,結果分別列于表1和2。采用Blackburn模型進行計算,對于(U0.7Pu0.3)O2-x燃料,當設計O/M比為1.97時,若在1 700 ℃、Ar-5%H2氣中燒結,則要求將氧分壓控制為7.70×10-6Pa,或將氧勢控制為-381.98 kJ/mol。對于(U0.75Pu0.25)O2-x燃料,當設計O/M比為1.97時,若在1 700 ℃、Ar-5%H2氣中燒結,則要求將氧分壓控制為3.05×10-6Pa,或將氧勢控制為-397.18 kJ/mol??梢姡_到相同的亞化學計量O/M比,PuO2含量高則對燒結氣氛的氧分壓控制要求可適當放松。

    采用點缺陷模型計算得到了MOX燃料的燒結氣氛氧分壓,對于(U0.7Pu0.3)O2-x燃料,當設計O/M比為1.97時,如果在1 700 ℃、Ar-5%H2氣中燒結,則要求將氧分壓控制在3.52×10-6Pa(而Blackburn模型計算值為7.70×10-6Pa),或將氧勢控制在-394.87 kJ/mol。對于(U0.75Pu0.25)O2-x燃料,當設計O/M比為1.97時,如果在1 700 ℃、Ar-5%H2氣中燒結,則要求將氧分壓控制在1.23×10-6Pa(而Blackburn模型計算值為3.05×10-6Pa),或將氧勢控制在-412.11 kJ/mol。

    表1 Blackburn模型計算得到的燒結氣氛的氧分壓和氧勢

    表2 點缺陷模型計算得到的燒結氣氛的氧分壓和氧勢

    續(xù)表2

    比較Blackburn模型與點缺陷模型的計算結果可發(fā)現(xiàn),兩種模型的計算結果變化規(guī)律基本一致(圖2),但Blackburn模型的氧勢和氧分壓隨O/M比增大而增大的速度較快。圖2顯示,當O/M比≤1.96時,Blackburn模型的氧勢和氧分壓計算結果均相應地低于點缺陷模型的;而當O/M比>1.96時,Blackburn模型的氧勢和氧分壓計算結果均相應地高于點缺陷模型的。今后應根據(jù)自有實驗數(shù)據(jù)選擇合適模型以用于理論預測我國MOX芯塊燒結氣氛氧分壓。由于目前文獻報道的點缺陷模型平衡常數(shù)有限,盡管國外研究者認為其與實驗符合較好,但在本文后續(xù)計算中仍選用Blackburn模型計算數(shù)據(jù)。

    圖2 Blackburn模型與點缺陷模型計算結果的比較

    2.2 露點計算

    如果燒結氣氛中存在水,存在2H2+O2→2H2O反應,則:

    (20)

    式中:p為各氣態(tài)物質的分壓;Keq為該反應平衡常數(shù),有:

    (21)

    ΔG0=ΔH0-TΔS0=-492 880+109.6T

    (22)

    pH2/pH2O=

    (23)

    根據(jù)Blackburn模型計算結果,以(U0.75-Pu0.25)O1.97為例進行計算分析,1 700 ℃燒結時要求將氧分壓控制在3.05×10-6Pa,于是計算得到:

    pH2/pH2O=

    若pH2=0.1 MPa,則pH2O=2 529.08 Pa。實際MOX芯塊燒結氣氛為Ar-5%H2混合氣體,若混合氣體的壓力為0.1 MPa,pH2=5 kPa,則pH2O=126.45 Pa,即混合氣體中的水分含量應控制為1 264.5 ppm。

    通過冷凝器冰浴[23]可調節(jié)氣體pH2/pH2O,對于冰面體系,飽和水蒸氣壓與溫度的關系可通過以下方程計算:

    lnpH2O=-6 024.528 2T-1+29.327 07+

    1.061 386 8×10-2T-1.319 882 5×

    10-5T2-0.493 825 77lnT

    (24)

    通過式(23)計算得到pH2O后,求解式(24)即可得到冷凝器的工作溫度即露點T為291.58 K(18.58 ℃)??梢?,將露點為18.58 ℃的0.1 MPa純H2氣或水分含量為1 264.5 ppm、露點為255.28 K(-17.72 ℃)的0.1 MPa Ar-5%H2混合氣通過冷凝器冰浴脫水后,再通入燒結爐內在1 700 ℃燒結,理論上可得到成分為(U0.75Pu0.25)O1.97的MOX芯塊。若燒結溫度提高至1 750 ℃,要使(U0.75Pu0.25)O2-x的O/M比控制在1.97,則要求將0.1 MPa Ar-5%H2氣體的水分含量控制在1 633.7 ppm、露點控制在258.02 K(-14.98 ℃)。

    采用上述同樣方法可計算得到1 750 ℃燒結具有不同成分和O/M比的(U0.75Pu0.25)O1.95、(U0.75Pu0.25)O1.96、(U0.75Pu0.25)O1.98、(U0.75Pu0.25)-O1.99、(U0.75Pu0.25)O1.995芯塊時,必須將0.1 MPa Ar-5%H2氣體的水分含量分別控制在370.4、739.8、4 403.6、17 855.4、43 064.8 ppm,露點分別控制在-30.10、-23.27、-3.77、13.83、26.16 ℃(圖3)。其他成分MOX燃料的燒結氣體水分含量和露點計算結果列于表3。

    表3 MOX芯塊燒結氣氛中水分含量及露點

    圖3 (U0.75Pu0.25)O2-x的O/M比與燒結氣體水分含量和露點的關系

    圖4 氧勢、燒結溫度、氣體水分含量和露點對(U0.8Pu0.2)O2-x燃料O/M比的影響[24]

    圖4所示為美國阿貢國家實驗室的(U0.8Pu0.2)O2-x燃料的O/M比與氧勢、溫度和氣體露點的關系[24]??梢?,若要求O/M=1.97,在1 700 ℃燒結時氣體的氧勢應控制在約-420 kJ/mol,較本文計算的氧勢偏低(Blackburn模型計算的氧勢為-397.18 kJ/mol,點缺陷模型的為-412.11 kJ/mol),露點約為-20 ℃。本文計算得到(U0.8Pu0.2)O1.97燃料在1 700 ℃燒結時,要求控制氣體水分含量為690.2 ppm、露點為-23.98 ℃;在1 750 ℃燒結時,要求控制氣體水分含量為893.4 ppm、露點為-21.35 ℃,這與文獻[24]的結論基本一致。俄羅斯專家指出快堆MOX芯塊在1 750 ℃、Ar-8%H2氣中燒結,當H2氣中H2O含量分別為0.03%、0.3%時,即氣體露點分別為-40、-16 ℃時,對應O/M比可分別達到1.95、1.98,但其MOX燃料的PuO2含量未知。本文計算得到(U0.75Pu0.25)O1.95、(U0.75Pu0.25)O1.98燃料分別在1 750 ℃、Ar-8%H2氣中燒結時,要求控制氣體露點分別為-32.34、-6.5 ℃,比俄羅斯專家提出的露點均高約8 ℃。

    3 結論

    1) Blackburn模型和點缺陷模型是兩種常用的核燃料氧勢模型,離子反應平衡常數(shù)和熱力學數(shù)據(jù)是影響氧勢模型精確度的主要因素。采用精確的平衡常數(shù)和熱力學數(shù)據(jù),以及大量實測數(shù)據(jù)對氧勢模型進行修正,可提高氧勢計算和O/M比預測的精確度。

    2) 在PuO2含量和燒結溫度一定的情況下,MOX芯塊的O/M比可通過控制氣相中氧分壓或露點來調整。將保護氣體通過冷凝器冰浴脫水后再進入燒結爐內,可定量調節(jié)氣體氧勢和精確控制芯塊的O/M比。

    3) 當要求(U0.75Pu0.25)O2-x燃料的O/M比為1.97時,若在1 750 ℃、0.1 MPa Ar-5%H2氣中燒結,采用Blackburn模型進行計算,則理論上要求將氧分壓控制在1.07×10-5Pa,或將氧勢控制在-386.15 kJ/mol。采用點缺陷模型計算,要求將氧分壓控制在0.7×10-5Pa,或將氧勢控制在-393.22 kJ/mol。

    4) (U0.75Pu0.25)O2-x燃料在1 750 ℃、0.1 MPa Ar-5%H2氣中燒結,當要求O/M比分別為1.95、1.96、1.97、1.98、1.99、1.995時,理論上應將氣體中的水分含量分別控制在370.4、739.8、16 33.7、4 403.6、17 855.4、43 064.8 ppm,或將氣體露點分別控制在-30.10、-23.27、-14.98、-3.77、13.83、26.16 ℃。

    參考文獻:

    [1] TOCI F, MANES L, EWART F T, et al. On the determination of oxygen potentials and O/M ratios in mixed U-Pu oxides by means of solid state galvanic cells[J]. Journal of Applied Electrochemistry, 1979, 9(3): 399-401.

    [2] KUTTY T R G, HEGDE P V, KAHN K B, et al. Densification behavior of UO2in six different atmospheres[J]. Journal of Nuclear Materials, 2002, 305(2-3): 159-168.

    [3] LAY K W, CARTER B E. Role of the O/U ratio on the sintering of UO2[J]. Journal of Nuclear Materials, 1969, 30(1-2): 74-87.

    [4] HARADA Y. UO2sintering in controlled oxygen atmospheres of three-stage process[J]. Journal of Nuclear Materials, 1997, 245(2-3): 217-223.

    [5] CHEVREL H, DEHAUDT P, FRANCOIS B, et al. Influence of surface phenomena during sintering of over stoichiometric uranium dioxide UO2+x[J]. Journal of Nuclear Materials, 1992, 189(2): 175-182.

    [6] ASSMANN H, DOERR W, PEEHS M. Oxide fuels with controlled microstructure[J]. Journal of the American Ceramic Society, 1984, 67(9): 631-636.

    [7] FUHRMAN N, HOWER L D, HOLDEN R B. Low-temperature sintering of uranium dioxide[J]. Journal of the American Ceramic Society, 1963, 46(3): 114-121.

    [8] 唐納德·奧蘭德. 核反應堆燃料元件基本問題[M]. 北京:原子能出版社,1983:231.

    [9] RAND M H, MARKIN T L. Some thermodynamic aspects of (U, Pu)O2solid solutions and their use as nuclear fuels[C]∥Thermodynamics of Nuclear Materials 1967 Symposium Proceedings. Vienna: [s. n.], 1968: 637-650.

    [10] FUJINO T, SATO N, YAMADA K. A refined analysis of oxygen potential of MyU1-yO2+x(M=M3+and M2+) by lattice statistics based on the grand partition function and the Flory methods[J]. Journal of Nuclear Materials, 1995, 223(1): 6-19.

    [11] STAN M, CRISTEA P, Sr WALLACE T C. Thermochemistry of defects and oxygen diffusion in PuO2-x[C]∥Proceedings of 11th Symposium on Thermodynamics of Nuclear Materials. Karlsruhe, Germany: [s. n.], 2004.

    [12] BLACKBURN P E, JOHNSON C E. Oxygen pressures over fast breeder reactor fuel (Ⅰ): A model for UO2±x[J]. Journal of Nuclear Materials, 1973, 46(3): 244-252.

    [13] 梁雪元. MOX燃料混合工藝及熱力學分析[D]. 北京:中國原子能科學研究院,2006.

    [14] STAN M. Theory-based models and simulations of materials for fusion[C]∥Fusion Power Renew Workshop. US: University of California, 2009.

    [15] HOCH M, FURMAN F J. Non-stoichiometry of UO2-CeO2: The system UO2-CeO2-CeO1.5at 900 to 1 200 ℃[C]∥Thermodynamics Proceedings of the Symposiumn on Thermodynamics With Emphasis on Nuclear Materials and Atomic Transport in Solids. Vienna: IAEA, 1966: 517-532.

    [16] BESMANN T M, LINDEMER T B. Chemical thermodynamic representation of UO2±x[C]∥ Proceedings of Conference on Thermodynamics. Hamilton, Canada: [s. n.], 1984.

    [17] BESMANN T M, LINDEMER T B. Chemical thermodynamic representation of PuO2-xand U1-zPuzOw[C]∥ Proceedings of Conference on Chemical Thermodynamics and 39 Calorimetry Conference. Hamilton, Canada: [s. n.], 1984.

    [18] KRISHINAIAH M V, SRIRAMAMURTI P. Computational model for the oxygen potentials of mixed uranium-plutonium oxide[J]. Journal of the American Ceramic Society, 1984, 67(8): 568-571.

    [19] VANA VARAMBAN S, GANESAN V, VASUDEVA RAO P R. Estimation of oxygen potential for U-Pu-O mixed oxide fuel[J]. Journal of Nuclear Materials, 2009, 389(1): 72-77.

    [20] KATO M, NAKAMICHI S, TAKEUCHI K, et al. Measurement of oxygen potential of (U0.8-Pu0.2)O2±xat 1 773 and 1 873 K, and its analysis based on point defect chemistry[J]. CALPHAD, 2011, 35(4): 623-626.

    [21] KATO M, KONASHI K, NAKAE N. Analysis of oxygen potential of (U0.7Pu0.3)O2±xand (U0.8-Pu0.2)O2±xbased on point defect chemistry[J]. Journal of Nuclear Materials, 2009, 389(1): 164-169.

    [22] KATO M. Oxygen potentials and defect chemistry in nonstoichiometric (U, Pu)O2[R]∥Stoichiometry and Materials Science: When Numbers Matter. [S. l.]: [s. n.], 2012: 203-218.

    [23] RASMUSSEN D E, BENECKE M W, JENTZEN W R, et al. Oxygen to metal ratio control during fabrication of mixed oxide fast breeder reactor fuel pellets, HEDL-SA-1763FP[R]. Atlanta: Hanford Engineering Development Laboratory, 1979.

    [24] JAVED N A. Thermodynamic behavior of (U, Pu) mixed-oxide fuels[J]. Journal of Nuclear Materials, 1973, 47(3): 336-344.

    猜你喜歡
    芯塊點缺陷平衡常數(shù)
    真空燒結U3Si2燃料芯塊的微觀組織與導熱性能
    化學平衡常數(shù)的計算
    常見壓強平衡常數(shù)Kp的幾種類型的計算
    常見壓強平衡常數(shù)K的幾種類型的計算
    金紅石型TiO2中四種點缺陷態(tài)研究
    場輔助燒結二氧化鈾基燃料芯塊研究進展
    陶瓷學報(2020年2期)2020-10-27 02:15:42
    環(huán)形燃料芯塊一維穩(wěn)態(tài)溫度場計算方法研究
    核技術(2020年6期)2020-06-15 03:01:14
    五大平衡常數(shù)的比較和應用
    Fe-Cr-Ni合金中點缺陷形成及相互作用的第一性原理研究
    GaN中質子輻照損傷的分子動力學模擬研究
    国产精品影院久久| 亚洲国产高清在线一区二区三 | 人人妻人人爽人人添夜夜欢视频| 在线永久观看黄色视频| АⅤ资源中文在线天堂| 欧美日韩中文字幕国产精品一区二区三区 | 麻豆成人av在线观看| 欧美精品啪啪一区二区三区| 九色亚洲精品在线播放| 男女午夜视频在线观看| 国产高清激情床上av| 久久久久久大精品| 18禁裸乳无遮挡免费网站照片 | 在线观看www视频免费| 人人妻人人爽人人添夜夜欢视频| 午夜影院日韩av| 999久久久国产精品视频| 欧美日韩精品网址| 国产熟女xx| 热re99久久国产66热| 在线国产一区二区在线| 亚洲人成电影免费在线| 欧美人与性动交α欧美精品济南到| 亚洲精品久久国产高清桃花| 国产精品久久久av美女十八| 韩国精品一区二区三区| 一进一出抽搐gif免费好疼| 免费在线观看影片大全网站| 淫妇啪啪啪对白视频| 亚洲第一电影网av| 黄色片一级片一级黄色片| 99国产精品一区二区蜜桃av| 村上凉子中文字幕在线| av免费在线观看网站| 亚洲一区二区三区不卡视频| 欧美黄色淫秽网站| 国产亚洲av嫩草精品影院| 久久久久久久久久久久大奶| 国产精品久久视频播放| 1024视频免费在线观看| 亚洲黑人精品在线| 狠狠狠狠99中文字幕| 免费女性裸体啪啪无遮挡网站| 1024香蕉在线观看| 午夜福利18| 国产成人系列免费观看| 美女大奶头视频| 亚洲自拍偷在线| 国产精品一区二区三区四区久久 | 热99re8久久精品国产| 国产精品 欧美亚洲| 久久精品国产综合久久久| 中文亚洲av片在线观看爽| 午夜福利18| 亚洲第一青青草原| 中亚洲国语对白在线视频| 人人妻人人澡人人看| 国产精品永久免费网站| 97超级碰碰碰精品色视频在线观看| 欧美日韩福利视频一区二区| 十八禁人妻一区二区| 国产一区二区三区综合在线观看| 老汉色av国产亚洲站长工具| 两个人免费观看高清视频| 最好的美女福利视频网| 亚洲自拍偷在线| 涩涩av久久男人的天堂| 国产野战对白在线观看| 日本五十路高清| 好男人电影高清在线观看| 精品人妻1区二区| 男女午夜视频在线观看| 一进一出好大好爽视频| 黄色视频,在线免费观看| 午夜成年电影在线免费观看| 精品福利观看| 啦啦啦免费观看视频1| 久久精品影院6| 视频在线观看一区二区三区| 欧美国产日韩亚洲一区| 国产精品久久电影中文字幕| 亚洲一区二区三区不卡视频| 妹子高潮喷水视频| 国产精品,欧美在线| 一二三四社区在线视频社区8| 露出奶头的视频| 日韩高清综合在线| 中文字幕色久视频| 可以在线观看毛片的网站| 精品高清国产在线一区| 欧美日韩精品网址| 免费少妇av软件| 91成年电影在线观看| 丝袜人妻中文字幕| 校园春色视频在线观看| 亚洲精品国产色婷婷电影| 免费在线观看视频国产中文字幕亚洲| 亚洲伊人色综图| 国产精品二区激情视频| www.www免费av| 成熟少妇高潮喷水视频| 手机成人av网站| 欧美一区二区精品小视频在线| 亚洲三区欧美一区| 熟妇人妻久久中文字幕3abv| 久久热在线av| 热99re8久久精品国产| 又大又爽又粗| 免费搜索国产男女视频| 天天添夜夜摸| 亚洲av熟女| 一个人观看的视频www高清免费观看 | 大陆偷拍与自拍| 日韩av在线大香蕉| 免费看十八禁软件| 别揉我奶头~嗯~啊~动态视频| 真人一进一出gif抽搐免费| av网站免费在线观看视频| 久久久国产成人精品二区| 久久热在线av| 亚洲欧美日韩另类电影网站| 一区福利在线观看| 免费女性裸体啪啪无遮挡网站| 国产伦人伦偷精品视频| 搡老妇女老女人老熟妇| 久久久国产成人精品二区| 窝窝影院91人妻| 亚洲人成网站在线播放欧美日韩| 久热这里只有精品99| 亚洲av第一区精品v没综合| 日韩一卡2卡3卡4卡2021年| 久久影院123| 老司机在亚洲福利影院| 国产精品久久久久久精品电影 | 精品国产乱码久久久久久男人| 亚洲一卡2卡3卡4卡5卡精品中文| 精品久久久久久成人av| 成人永久免费在线观看视频| 乱人伦中国视频| 亚洲,欧美精品.| 美女午夜性视频免费| avwww免费| 日本欧美视频一区| 国产高清激情床上av| 国产欧美日韩精品亚洲av| 免费观看人在逋| 一级片免费观看大全| 一区在线观看完整版| 久久久久久亚洲精品国产蜜桃av| 久久久久国产一级毛片高清牌| 91老司机精品| 91大片在线观看| 精品久久久久久久人妻蜜臀av | 亚洲欧洲精品一区二区精品久久久| 亚洲色图综合在线观看| 国产伦一二天堂av在线观看| 村上凉子中文字幕在线| 夜夜爽天天搞| 九色国产91popny在线| 精品不卡国产一区二区三区| 美女高潮喷水抽搐中文字幕| 国产亚洲av嫩草精品影院| 国产片内射在线| 亚洲午夜精品一区,二区,三区| 美女免费视频网站| 日本 欧美在线| 欧洲精品卡2卡3卡4卡5卡区| 亚洲专区国产一区二区| 欧美黑人精品巨大| 国产精品 国内视频| 三级毛片av免费| 99精品久久久久人妻精品| 老司机午夜福利在线观看视频| 看黄色毛片网站| 手机成人av网站| 亚洲va日本ⅴa欧美va伊人久久| 悠悠久久av| av天堂在线播放| 日本欧美视频一区| 免费搜索国产男女视频| 亚洲色图 男人天堂 中文字幕| 国产精品一区二区在线不卡| 老熟妇仑乱视频hdxx| 久久久久久人人人人人| 久久亚洲精品不卡| 亚洲aⅴ乱码一区二区在线播放 | 欧美精品啪啪一区二区三区| 精品日产1卡2卡| 午夜激情av网站| 99国产精品免费福利视频| 精品熟女少妇八av免费久了| 男人的好看免费观看在线视频 | 午夜久久久在线观看| 啦啦啦韩国在线观看视频| 国产亚洲精品综合一区在线观看 | 深夜精品福利| 大陆偷拍与自拍| 18美女黄网站色大片免费观看| 国产国语露脸激情在线看| 色尼玛亚洲综合影院| 脱女人内裤的视频| 好看av亚洲va欧美ⅴa在| 国产99白浆流出| 美女午夜性视频免费| 亚洲av第一区精品v没综合| 久久久久国产一级毛片高清牌| 亚洲色图综合在线观看| 欧美不卡视频在线免费观看 | 欧美乱码精品一区二区三区| 日韩精品免费视频一区二区三区| 久久亚洲真实| 国产成人av激情在线播放| 好男人电影高清在线观看| 国产精品av久久久久免费| 欧美中文日本在线观看视频| 91大片在线观看| 久久人妻av系列| 国产精品亚洲av一区麻豆| 窝窝影院91人妻| 9191精品国产免费久久| 欧美成人午夜精品| 久久国产精品男人的天堂亚洲| av有码第一页| 日本免费a在线| 后天国语完整版免费观看| 欧美中文综合在线视频| 亚洲在线自拍视频| 欧美不卡视频在线免费观看 | 亚洲无线在线观看| 嫁个100分男人电影在线观看| 国产区一区二久久| 欧美日韩亚洲国产一区二区在线观看| 亚洲在线自拍视频| 国产亚洲精品综合一区在线观看 | 老汉色∧v一级毛片| 国产不卡一卡二| 在线永久观看黄色视频| 久久久国产精品麻豆| 国产av精品麻豆| 老司机深夜福利视频在线观看| 国产成+人综合+亚洲专区| 亚洲一区中文字幕在线| 国产亚洲av高清不卡| 999久久久精品免费观看国产| 国产三级在线视频| 午夜福利高清视频| 亚洲午夜精品一区,二区,三区| 中文亚洲av片在线观看爽| av天堂在线播放| 老司机靠b影院| 嫩草影视91久久| 国产精品香港三级国产av潘金莲| 亚洲成av片中文字幕在线观看| 日韩av在线大香蕉| 午夜福利在线观看吧| 国产欧美日韩综合在线一区二区| 久久国产精品人妻蜜桃| 黑人欧美特级aaaaaa片| 亚洲av五月六月丁香网| 国产成人欧美| 亚洲av成人不卡在线观看播放网| 免费久久久久久久精品成人欧美视频| 中出人妻视频一区二区| 人妻丰满熟妇av一区二区三区| 男女午夜视频在线观看| 久久久水蜜桃国产精品网| 三级毛片av免费| 精品国产乱子伦一区二区三区| 男人舔女人的私密视频| 久久人妻熟女aⅴ| 法律面前人人平等表现在哪些方面| 国产99白浆流出| 国产精品日韩av在线免费观看 | 淫妇啪啪啪对白视频| 亚洲全国av大片| 一夜夜www| 久久午夜综合久久蜜桃| 亚洲精品在线观看二区| 久久九九热精品免费| 99久久99久久久精品蜜桃| 亚洲成人国产一区在线观看| 天堂√8在线中文| 巨乳人妻的诱惑在线观看| 欧美激情 高清一区二区三区| 国产麻豆成人av免费视频| 精品一区二区三区av网在线观看| 亚洲情色 制服丝袜| 国产精品亚洲av一区麻豆| 不卡av一区二区三区| 精品人妻1区二区| 亚洲中文av在线| 国产色视频综合| 午夜免费鲁丝| 日韩精品青青久久久久久| 国产精品爽爽va在线观看网站 | 国产精品九九99| 欧美丝袜亚洲另类 | 亚洲欧洲精品一区二区精品久久久| 妹子高潮喷水视频| 国产一区二区三区在线臀色熟女| 中出人妻视频一区二区| 日本黄色视频三级网站网址| 亚洲精品粉嫩美女一区| 亚洲免费av在线视频| 久热爱精品视频在线9| 国产精品自产拍在线观看55亚洲| 亚洲av片天天在线观看| 露出奶头的视频| 亚洲色图综合在线观看| 在线天堂中文资源库| 欧美久久黑人一区二区| 久热这里只有精品99| 亚洲av电影不卡..在线观看| 精品国产一区二区三区四区第35| 亚洲成国产人片在线观看| 精品电影一区二区在线| 国产国语露脸激情在线看| 免费在线观看影片大全网站| 如日韩欧美国产精品一区二区三区| 欧美精品亚洲一区二区| 色播亚洲综合网| 在线十欧美十亚洲十日本专区| 国产成人精品在线电影| 丝袜美足系列| 91成年电影在线观看| 国产精品野战在线观看| 女性生殖器流出的白浆| 麻豆一二三区av精品| 99在线人妻在线中文字幕| 亚洲在线自拍视频| 国产精品一区二区三区四区久久 | 国产成人免费无遮挡视频| 亚洲久久久国产精品| 精品久久久久久成人av| 男女做爰动态图高潮gif福利片 | 婷婷精品国产亚洲av在线| 99香蕉大伊视频| 午夜两性在线视频| 成人亚洲精品一区在线观看| 亚洲精品国产精品久久久不卡| 精品免费久久久久久久清纯| bbb黄色大片| 久久精品亚洲精品国产色婷小说| 国产又爽黄色视频| 十八禁人妻一区二区| 久久精品国产综合久久久| 国产免费av片在线观看野外av| 国产精品美女特级片免费视频播放器 | 人妻久久中文字幕网| 欧美 亚洲 国产 日韩一| 欧美激情 高清一区二区三区| 很黄的视频免费| 亚洲欧美日韩另类电影网站| 免费高清在线观看日韩| 中文字幕最新亚洲高清| 日韩欧美国产一区二区入口| 一级,二级,三级黄色视频| 真人做人爱边吃奶动态| 99re在线观看精品视频| 性少妇av在线| 亚洲成人精品中文字幕电影| 中文字幕色久视频| 一区福利在线观看| 午夜久久久在线观看| 亚洲精品久久成人aⅴ小说| 熟女少妇亚洲综合色aaa.| 午夜福利欧美成人| 国产精品 欧美亚洲| 久久婷婷人人爽人人干人人爱 | 国产精品九九99| 可以在线观看毛片的网站| 在线天堂中文资源库| 亚洲精品av麻豆狂野| 国产一区二区激情短视频| 国产99久久九九免费精品| АⅤ资源中文在线天堂| 成人精品一区二区免费| 久久中文字幕一级| 又黄又爽又免费观看的视频| 最近最新免费中文字幕在线| 亚洲 欧美 日韩 在线 免费| 一级作爱视频免费观看| 亚洲 国产 在线| 久久这里只有精品19| 国产av在哪里看| 日韩一卡2卡3卡4卡2021年| av天堂在线播放| 操美女的视频在线观看| 国产一区在线观看成人免费| 久久久久久亚洲精品国产蜜桃av| 12—13女人毛片做爰片一| 亚洲精品久久国产高清桃花| 午夜久久久久精精品| 国产精品久久电影中文字幕| 好男人电影高清在线观看| 国产成人影院久久av| 怎么达到女性高潮| 日本欧美视频一区| 国产欧美日韩综合在线一区二区| ponron亚洲| 国产成人影院久久av| 欧美 亚洲 国产 日韩一| 亚洲精品粉嫩美女一区| 亚洲情色 制服丝袜| 久久婷婷人人爽人人干人人爱 | 成人国产综合亚洲| 免费高清在线观看日韩| 啦啦啦观看免费观看视频高清 | 亚洲精品国产色婷婷电影| 国产亚洲精品第一综合不卡| 波多野结衣一区麻豆| 久久青草综合色| 午夜福利一区二区在线看| 国产亚洲欧美98| 18禁观看日本| 亚洲精品久久国产高清桃花| 久久午夜综合久久蜜桃| 亚洲国产毛片av蜜桃av| 国产麻豆成人av免费视频| 777久久人妻少妇嫩草av网站| 国产精品精品国产色婷婷| 国产精品一区二区免费欧美| 国产av一区二区精品久久| 欧美日韩福利视频一区二区| 国产精品亚洲一级av第二区| 亚洲国产精品成人综合色| 真人做人爱边吃奶动态| 长腿黑丝高跟| 午夜福利成人在线免费观看| 精品电影一区二区在线| 两个人视频免费观看高清| 日日夜夜操网爽| 午夜福利一区二区在线看| 久久久国产成人精品二区| 亚洲自偷自拍图片 自拍| ponron亚洲| 亚洲九九香蕉| 亚洲精品av麻豆狂野| 欧美在线黄色| 18禁黄网站禁片午夜丰满| 黑人巨大精品欧美一区二区mp4| 欧美老熟妇乱子伦牲交| 久久久精品国产亚洲av高清涩受| 999久久久精品免费观看国产| 制服人妻中文乱码| 99国产精品一区二区三区| 久久精品影院6| 人妻丰满熟妇av一区二区三区| 国产av一区二区精品久久| 法律面前人人平等表现在哪些方面| 亚洲国产毛片av蜜桃av| 午夜精品在线福利| 动漫黄色视频在线观看| 电影成人av| 色老头精品视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 美女高潮到喷水免费观看| 精品久久久久久成人av| 日本欧美视频一区| av片东京热男人的天堂| 久久久久久人人人人人| 久久性视频一级片| 国产亚洲av高清不卡| av网站免费在线观看视频| 女人爽到高潮嗷嗷叫在线视频| 1024香蕉在线观看| 少妇的丰满在线观看| 午夜精品国产一区二区电影| 真人一进一出gif抽搐免费| 一二三四在线观看免费中文在| 久久中文字幕一级| 亚洲专区字幕在线| 免费观看人在逋| 亚洲av成人一区二区三| 美女高潮到喷水免费观看| xxx96com| 亚洲精品在线观看二区| 这个男人来自地球电影免费观看| 麻豆一二三区av精品| 免费av毛片视频| 三级毛片av免费| 一进一出抽搐gif免费好疼| 国内毛片毛片毛片毛片毛片| 久久久久久久久中文| 欧美激情极品国产一区二区三区| 国产精品电影一区二区三区| 欧美黑人欧美精品刺激| 亚洲国产中文字幕在线视频| 亚洲国产看品久久| 精品欧美一区二区三区在线| 久久人人97超碰香蕉20202| 国产av又大| 成人特级黄色片久久久久久久| 麻豆av在线久日| 黄网站色视频无遮挡免费观看| 一本久久中文字幕| 亚洲一区二区三区色噜噜| 一本久久中文字幕| 在线观看www视频免费| 国产色视频综合| 国产精品亚洲美女久久久| 琪琪午夜伦伦电影理论片6080| 成人18禁在线播放| 身体一侧抽搐| 午夜精品在线福利| 久久久久久人人人人人| 黄色片一级片一级黄色片| 91字幕亚洲| 很黄的视频免费| 色在线成人网| 欧美日韩亚洲国产一区二区在线观看| 在线观看免费视频网站a站| 99久久综合精品五月天人人| 亚洲少妇的诱惑av| 非洲黑人性xxxx精品又粗又长| 成人国产一区最新在线观看| 国产欧美日韩综合在线一区二区| 大型黄色视频在线免费观看| 麻豆国产av国片精品| 欧美亚洲日本最大视频资源| 久久香蕉激情| 久久久久久人人人人人| 在线观看日韩欧美| 久久人人97超碰香蕉20202| 国产成人av教育| 嫩草影视91久久| 十分钟在线观看高清视频www| 精品福利观看| 啦啦啦 在线观看视频| 欧美最黄视频在线播放免费| 日韩欧美一区视频在线观看| 午夜免费激情av| 精品国产亚洲在线| 香蕉国产在线看| 999久久久精品免费观看国产| 别揉我奶头~嗯~啊~动态视频| 亚洲精品在线观看二区| 可以免费在线观看a视频的电影网站| 国产在线观看jvid| 国产亚洲欧美精品永久| 久热这里只有精品99| 精品一区二区三区av网在线观看| 久久久久久久午夜电影| 韩国av一区二区三区四区| 91av网站免费观看| www日本在线高清视频| 国产成人影院久久av| 成人av一区二区三区在线看| 国产欧美日韩综合在线一区二区| 日韩欧美国产一区二区入口| 首页视频小说图片口味搜索| 男女下面进入的视频免费午夜 | 黄片大片在线免费观看| 久久久久久免费高清国产稀缺| 国产欧美日韩一区二区精品| 大型黄色视频在线免费观看| 两人在一起打扑克的视频| 久久精品国产亚洲av高清一级| 韩国av一区二区三区四区| 日本免费a在线| 国产精品免费视频内射| 91国产中文字幕| 成人国产一区最新在线观看| 精品欧美一区二区三区在线| 男女做爰动态图高潮gif福利片 | 日韩免费av在线播放| 男人的好看免费观看在线视频 | 国内久久婷婷六月综合欲色啪| 亚洲精品在线观看二区| av片东京热男人的天堂| av在线播放免费不卡| 国产日韩一区二区三区精品不卡| 久久久久久国产a免费观看| 好看av亚洲va欧美ⅴa在| 久久久久久久精品吃奶| 欧美激情久久久久久爽电影 | 亚洲欧美精品综合一区二区三区| 久久久久国内视频| 丝袜在线中文字幕| 久久久精品国产亚洲av高清涩受| 亚洲国产精品999在线| www.精华液| 午夜福利,免费看| 波多野结衣巨乳人妻| 黄频高清免费视频| 亚洲欧美精品综合一区二区三区| tocl精华| 亚洲精品一区av在线观看| 欧美成人免费av一区二区三区| 国产精品久久久久久精品电影 | 熟妇人妻久久中文字幕3abv| 麻豆av在线久日| 欧美一区二区精品小视频在线| 国产精品秋霞免费鲁丝片| 国产又爽黄色视频| 亚洲情色 制服丝袜| 久久精品人人爽人人爽视色| 看免费av毛片| 日本免费a在线| 在线观看66精品国产| 一级作爱视频免费观看| 精品久久久久久,| 国产精品一区二区在线不卡| 国产午夜精品久久久久久| 午夜福利影视在线免费观看| 一夜夜www| 亚洲三区欧美一区| 日本撒尿小便嘘嘘汇集6| 国产主播在线观看一区二区| 欧美大码av| 国产一区二区三区在线臀色熟女| 激情在线观看视频在线高清| 欧美乱妇无乱码| av中文乱码字幕在线| 在线播放国产精品三级|