• <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| 最好的美女福利视频网| 亚洲国产欧美人成| 国产精品久久久久久精品电影小说 | 爱豆传媒免费全集在线观看| 插阴视频在线观看视频| 日韩欧美国产在线观看| 夜夜看夜夜爽夜夜摸| 久久精品久久久久久噜噜老黄 | 亚洲成a人片在线一区二区| 性欧美人与动物交配| 老师上课跳d突然被开到最大视频| 日韩成人伦理影院| 天堂中文最新版在线下载 | 欧美丝袜亚洲另类| 国产精品不卡视频一区二区| 日日啪夜夜撸| 激情 狠狠 欧美| 国产亚洲av嫩草精品影院| 色视频www国产| 精品久久久久久久人妻蜜臀av| 国产黄片美女视频| 久久99热这里只有精品18| а√天堂www在线а√下载| 免费在线观看成人毛片| 亚洲精华国产精华液的使用体验 | 国产69精品久久久久777片| 天堂av国产一区二区熟女人妻| 亚洲国产高清在线一区二区三| 日韩中字成人| 色播亚洲综合网| 国产精品伦人一区二区| 欧美激情久久久久久爽电影| 男女下面进入的视频免费午夜| 午夜亚洲福利在线播放| 特大巨黑吊av在线直播| 一个人看视频在线观看www免费| 亚洲精品乱码久久久久久按摩| 亚洲av成人av| 精品人妻熟女av久视频| 成人午夜精彩视频在线观看| 此物有八面人人有两片| 乱人视频在线观看| 啦啦啦啦在线视频资源| 国产精品一区二区三区四区免费观看| 搡女人真爽免费视频火全软件| 亚洲av熟女| 99久国产av精品| 亚洲,欧美,日韩| 久久99精品国语久久久| 国产黄片视频在线免费观看| 九九爱精品视频在线观看| www.av在线官网国产| 国产一区二区三区av在线 | 成人亚洲欧美一区二区av| 午夜a级毛片| 欧美激情在线99| 国产精品久久久久久久久免| 人人妻人人澡人人爽人人夜夜 | 黄色日韩在线| 91麻豆精品激情在线观看国产| 亚洲精品久久国产高清桃花| 99久久成人亚洲精品观看| 日本熟妇午夜| 26uuu在线亚洲综合色| 69人妻影院| 国产一级毛片在线| h日本视频在线播放| 男人舔女人下体高潮全视频| 一区二区三区免费毛片| 精品久久久噜噜| 天堂中文最新版在线下载 | 国产久久久一区二区三区| 亚洲色图av天堂| 久久精品久久久久久噜噜老黄 | 久久精品国产鲁丝片午夜精品| 亚洲色图av天堂| 男女做爰动态图高潮gif福利片| 嫩草影院新地址| 国产精品伦人一区二区| 免费无遮挡裸体视频| 日韩精品有码人妻一区| 99久久精品国产国产毛片| 又黄又爽又刺激的免费视频.| 欧美日韩一区二区视频在线观看视频在线 | 搡老妇女老女人老熟妇| 国产v大片淫在线免费观看| 不卡视频在线观看欧美| 久久久久久久久大av| 亚洲精品影视一区二区三区av| 别揉我奶头 嗯啊视频| 国产精品1区2区在线观看.| 日韩制服骚丝袜av| 国产淫片久久久久久久久| 精品国内亚洲2022精品成人| 国产三级在线视频| 最近中文字幕高清免费大全6| 如何舔出高潮| 亚洲国产高清在线一区二区三| 国产精品人妻久久久影院| 日本欧美国产在线视频| 91精品一卡2卡3卡4卡| 亚洲真实伦在线观看| 超碰av人人做人人爽久久| 免费看av在线观看网站| 欧美激情久久久久久爽电影| 99久久精品热视频| 春色校园在线视频观看| 国产在线男女| 精品国产三级普通话版| 久久九九热精品免费| 亚洲欧美日韩高清在线视频| 国产精品久久久久久亚洲av鲁大| 女人十人毛片免费观看3o分钟| 亚洲美女搞黄在线观看| 国产精品一区二区性色av| 久久久精品94久久精品| 亚洲精品456在线播放app| 免费无遮挡裸体视频| 久久久久久伊人网av| 久久国内精品自在自线图片| 久久久久免费精品人妻一区二区| 国产蜜桃级精品一区二区三区| 一本久久中文字幕| 亚洲欧美成人精品一区二区| avwww免费| 男人的好看免费观看在线视频| 国产色婷婷99| 欧美又色又爽又黄视频| 毛片一级片免费看久久久久| 成人欧美大片| 97人妻精品一区二区三区麻豆| 99久久无色码亚洲精品果冻| a级一级毛片免费在线观看| 此物有八面人人有两片| 亚洲精品久久国产高清桃花| 91精品一卡2卡3卡4卡| 激情 狠狠 欧美| 男人狂女人下面高潮的视频| 亚洲一区二区三区色噜噜| 亚州av有码| 少妇人妻精品综合一区二区 | 国产一区二区在线av高清观看| 99热全是精品| 国产亚洲精品久久久久久毛片| 欧美精品国产亚洲| 一级av片app| 婷婷亚洲欧美| 国产亚洲91精品色在线| 一级av片app| 三级男女做爰猛烈吃奶摸视频| 国产精品人妻久久久影院| 少妇人妻精品综合一区二区 | 国产精品爽爽va在线观看网站| 91精品国产九色| 九九在线视频观看精品| 国产精品女同一区二区软件| 中文字幕人妻熟人妻熟丝袜美| 1000部很黄的大片| 丰满的人妻完整版| 午夜免费激情av| 免费观看精品视频网站| 日韩中字成人| 亚洲人与动物交配视频| 在线国产一区二区在线| 亚洲精品乱码久久久久久按摩| 搡老妇女老女人老熟妇| 婷婷色av中文字幕| 精品少妇黑人巨大在线播放 | 最后的刺客免费高清国语| 十八禁国产超污无遮挡网站| av在线观看视频网站免费| 99国产精品一区二区蜜桃av| 国产精品野战在线观看| 久久这里只有精品中国| 在线免费观看不下载黄p国产| 日韩,欧美,国产一区二区三区 | 一本久久中文字幕| 精品少妇黑人巨大在线播放 | 美女cb高潮喷水在线观看| 亚洲av中文字字幕乱码综合| 免费av毛片视频| 午夜爱爱视频在线播放| 欧美激情久久久久久爽电影| 亚洲无线在线观看| 国产乱人偷精品视频| 国产精品综合久久久久久久免费| 六月丁香七月| 亚洲丝袜综合中文字幕| 26uuu在线亚洲综合色| 日本欧美国产在线视频| 99久久九九国产精品国产免费| 99热只有精品国产| 一级av片app| av.在线天堂| 激情 狠狠 欧美| 亚洲国产精品国产精品| 人人妻人人澡欧美一区二区| 一本一本综合久久| 免费无遮挡裸体视频| 精品久久久久久久久av| 别揉我奶头 嗯啊视频| 欧美精品一区二区大全| 欧美日韩乱码在线| 亚洲天堂国产精品一区在线| 97人妻精品一区二区三区麻豆| 美女cb高潮喷水在线观看| 日韩成人伦理影院| 天堂√8在线中文| 免费看日本二区| 成年女人看的毛片在线观看| 两个人视频免费观看高清| 久久久久久久亚洲中文字幕| 国产不卡一卡二| 亚洲欧美日韩高清专用| 少妇熟女aⅴ在线视频| 国产成年人精品一区二区| 国产视频内射| 亚洲成人久久性| 狠狠狠狠99中文字幕| 联通29元200g的流量卡| 国产高清激情床上av| 国产黄色视频一区二区在线观看 | 搡女人真爽免费视频火全软件| 久99久视频精品免费| 大香蕉久久网| 免费大片18禁| 赤兔流量卡办理| 久久精品国产亚洲av香蕉五月| 亚洲人与动物交配视频| 午夜福利在线观看免费完整高清在 | 99热这里只有是精品50| 国内揄拍国产精品人妻在线| 干丝袜人妻中文字幕| 亚洲国产高清在线一区二区三| 免费观看人在逋| 久久人妻av系列| 亚洲自拍偷在线| 一本一本综合久久| 在线a可以看的网站| 淫秽高清视频在线观看| 精品国内亚洲2022精品成人| 亚洲av男天堂| 女的被弄到高潮叫床怎么办| 欧美日本亚洲视频在线播放| 毛片女人毛片| 亚洲av二区三区四区| 日韩强制内射视频| 国产探花在线观看一区二区| 国产精品,欧美在线| 青春草国产在线视频 | 两个人视频免费观看高清| 亚洲丝袜综合中文字幕| 深夜精品福利| 内地一区二区视频在线| av福利片在线观看| 欧美日韩在线观看h| 欧美性猛交黑人性爽| 久久久久网色| 国产精品野战在线观看| 成年女人看的毛片在线观看| 3wmmmm亚洲av在线观看| or卡值多少钱| 国产视频内射| 亚洲性久久影院| 99国产极品粉嫩在线观看| 中文字幕熟女人妻在线| 久久99热6这里只有精品| 搞女人的毛片| eeuss影院久久| 在线观看美女被高潮喷水网站| 女人被狂操c到高潮| 一个人看视频在线观看www免费| 黄片wwwwww| 国产成人a∨麻豆精品| 国产av一区在线观看免费| 人体艺术视频欧美日本| 爱豆传媒免费全集在线观看| 久久久久网色| 在线播放无遮挡| 亚洲国产精品sss在线观看| 美女xxoo啪啪120秒动态图| 亚洲av成人精品一区久久| 国产精品人妻久久久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 波多野结衣巨乳人妻| 99久久中文字幕三级久久日本| 中出人妻视频一区二区| 欧美日韩在线观看h| 能在线免费观看的黄片| 永久网站在线| 久久国内精品自在自线图片| 插逼视频在线观看| 久久午夜福利片| 精品少妇黑人巨大在线播放 | 国产精品嫩草影院av在线观看| 大香蕉久久网| av视频在线观看入口| 亚洲国产日韩欧美精品在线观看| 综合色av麻豆| 精品久久久久久久人妻蜜臀av| 爱豆传媒免费全集在线观看| 国产亚洲av片在线观看秒播厂 | 日本熟妇午夜| av在线天堂中文字幕| 一级毛片aaaaaa免费看小| 国产私拍福利视频在线观看| 级片在线观看| 99热网站在线观看| 又爽又黄a免费视频| 3wmmmm亚洲av在线观看| 晚上一个人看的免费电影| 神马国产精品三级电影在线观看| 国产人妻一区二区三区在| av.在线天堂| 好男人视频免费观看在线| 国产精品人妻久久久久久| 国产国拍精品亚洲av在线观看| 99久久成人亚洲精品观看| 少妇被粗大猛烈的视频| 亚洲av电影不卡..在线观看| eeuss影院久久| 成人鲁丝片一二三区免费| 国产爱豆传媒在线观看| 成人av在线播放网站| 成年女人看的毛片在线观看| 久久九九热精品免费| 三级男女做爰猛烈吃奶摸视频| av视频在线观看入口| 亚洲性久久影院| 国产视频首页在线观看| kizo精华| av在线天堂中文字幕| 午夜视频国产福利| 久久久精品94久久精品| 日韩一区二区视频免费看| 亚洲自偷自拍三级| 国内精品久久久久精免费| 色综合亚洲欧美另类图片| 一卡2卡三卡四卡精品乱码亚洲| 欧美日本亚洲视频在线播放| 干丝袜人妻中文字幕| 一进一出抽搐gif免费好疼| 午夜亚洲福利在线播放| 人妻久久中文字幕网| 久久久久久久久久久免费av| 亚洲美女视频黄频| 成人无遮挡网站| 久久久久久久亚洲中文字幕| av天堂在线播放| 99久久无色码亚洲精品果冻| 国产探花在线观看一区二区| 黄片无遮挡物在线观看| 日韩高清综合在线| 欧美变态另类bdsm刘玥| 日本与韩国留学比较| 性插视频无遮挡在线免费观看| 日韩大尺度精品在线看网址| 精品人妻视频免费看| 美女cb高潮喷水在线观看| 日韩成人伦理影院| 久久久久久久久久黄片| 秋霞在线观看毛片| 免费无遮挡裸体视频| 丰满乱子伦码专区| 国产精品久久电影中文字幕| 亚洲精品亚洲一区二区| 国模一区二区三区四区视频| 一级二级三级毛片免费看| 久久午夜亚洲精品久久| 亚洲国产日韩欧美精品在线观看| 国产一区二区在线观看日韩| 日日摸夜夜添夜夜爱| 真实男女啪啪啪动态图| 亚洲精品日韩在线中文字幕 | 久久久a久久爽久久v久久| 91午夜精品亚洲一区二区三区| 天天躁日日操中文字幕| 免费人成在线观看视频色| 在线观看一区二区三区| 51国产日韩欧美| 九九在线视频观看精品| 在线观看美女被高潮喷水网站| 成人美女网站在线观看视频| 日韩欧美精品v在线| 99久久无色码亚洲精品果冻| 中国国产av一级| 人妻系列 视频| 国产精品野战在线观看| 国产精品一区二区性色av| 精品不卡国产一区二区三区| 三级经典国产精品| 亚洲真实伦在线观看| 国产成人精品婷婷| 我要搜黄色片| 久久精品国产99精品国产亚洲性色| 亚洲国产日韩欧美精品在线观看| 国产精品福利在线免费观看| 黄色视频,在线免费观看| 免费在线观看成人毛片| 国产伦在线观看视频一区| 国产 一区精品| 久久6这里有精品| 亚洲第一电影网av| 女同久久另类99精品国产91| 在线免费观看的www视频| 国产亚洲精品久久久com| 亚洲在线自拍视频| 久久久久久九九精品二区国产| 久久久久久大精品| 国产视频内射| 少妇高潮的动态图| 国产精品一及| 在现免费观看毛片| 看十八女毛片水多多多| 国产一区二区三区在线臀色熟女| 国产日韩欧美在线精品| 男女做爰动态图高潮gif福利片| 久久精品国产99精品国产亚洲性色| 日韩一本色道免费dvd| 给我免费播放毛片高清在线观看| 岛国毛片在线播放| 成年女人永久免费观看视频| 国产又黄又爽又无遮挡在线| 成人国产麻豆网| 国产69精品久久久久777片| 22中文网久久字幕| 久久久色成人| 高清午夜精品一区二区三区 | 赤兔流量卡办理| 亚洲18禁久久av| av黄色大香蕉| 精品国产三级普通话版| 又爽又黄无遮挡网站| 欧美一区二区精品小视频在线| 男人的好看免费观看在线视频| 伦理电影大哥的女人| 一级二级三级毛片免费看| 欧美区成人在线视频| av在线天堂中文字幕| 美女被艹到高潮喷水动态| 免费不卡的大黄色大毛片视频在线观看 | .国产精品久久| 黄色视频,在线免费观看| 国产探花在线观看一区二区| 99久久精品热视频| 蜜桃亚洲精品一区二区三区| 欧美+亚洲+日韩+国产| 日本欧美国产在线视频| 老熟妇乱子伦视频在线观看| 国产一区二区亚洲精品在线观看| 尤物成人国产欧美一区二区三区| 国产精品av视频在线免费观看| 能在线免费观看的黄片| 欧美日韩一区二区视频在线观看视频在线 | 亚洲成人中文字幕在线播放| 午夜激情福利司机影院| 国产老妇女一区| 国产黄色小视频在线观看| 欧美最黄视频在线播放免费| 国产成年人精品一区二区| 国产综合懂色| 97人妻精品一区二区三区麻豆| 亚洲成a人片在线一区二区| 亚洲激情五月婷婷啪啪| 欧美xxxx黑人xx丫x性爽| 国产精品国产高清国产av| 国产黄a三级三级三级人| 久久久久久久亚洲中文字幕| 日本-黄色视频高清免费观看| 噜噜噜噜噜久久久久久91| 最近最新中文字幕大全电影3| 五月伊人婷婷丁香| 99热这里只有是精品在线观看| 亚洲一区二区三区色噜噜| 18禁在线播放成人免费| 国产成人a区在线观看| 最新中文字幕久久久久| 久久久久九九精品影院| 国产色婷婷99| 国产精品一及| 老司机福利观看| av免费观看日本| 性欧美人与动物交配| 欧美一区二区精品小视频在线| 亚洲欧美清纯卡通| 成年版毛片免费区| 精品免费久久久久久久清纯| 中国美女看黄片| 国产白丝娇喘喷水9色精品| 我要搜黄色片| 亚洲无线观看免费| 国产视频首页在线观看| 免费一级毛片在线播放高清视频| 99热只有精品国产| 观看免费一级毛片| 爱豆传媒免费全集在线观看| 乱码一卡2卡4卡精品| 国产精品久久久久久精品电影| 校园春色视频在线观看| 日本色播在线视频| 国产日本99.免费观看| 日本五十路高清| 九九爱精品视频在线观看| 天堂√8在线中文| 亚洲在久久综合| 亚洲欧美清纯卡通| 欧美日本亚洲视频在线播放| 午夜久久久久精精品| 九九热线精品视视频播放| 欧美不卡视频在线免费观看| 午夜免费激情av| 日产精品乱码卡一卡2卡三| 亚洲精品乱码久久久v下载方式| 亚洲国产欧洲综合997久久,| 国产探花极品一区二区| 久久午夜亚洲精品久久| 久久国产乱子免费精品| 天堂网av新在线| av在线天堂中文字幕| 欧美日韩综合久久久久久| 在线免费十八禁| 亚洲精品成人久久久久久| 欧美一级a爱片免费观看看| 国产精品女同一区二区软件| 高清午夜精品一区二区三区 | av免费观看日本| 国产日韩欧美在线精品| 一夜夜www| 国产成人精品久久久久久| 女人被狂操c到高潮| 国产成人91sexporn| 亚洲人成网站在线观看播放| 18禁在线无遮挡免费观看视频| 日韩av不卡免费在线播放| 蜜桃久久精品国产亚洲av| 天堂√8在线中文| 亚洲国产精品合色在线| 精品欧美国产一区二区三| 亚洲自拍偷在线| 国产精品麻豆人妻色哟哟久久 | 国产一区二区在线av高清观看| 欧美激情久久久久久爽电影| 午夜爱爱视频在线播放| 久久亚洲精品不卡| 看片在线看免费视频| 欧美成人a在线观看| 亚洲精品亚洲一区二区| 级片在线观看| 偷拍熟女少妇极品色| 精品熟女少妇av免费看| 亚洲一区高清亚洲精品| 波多野结衣高清作品| 人妻制服诱惑在线中文字幕| 1024手机看黄色片| 3wmmmm亚洲av在线观看| 天堂av国产一区二区熟女人妻| 欧美成人a在线观看| 神马国产精品三级电影在线观看| 日韩欧美国产在线观看| 青青草视频在线视频观看| 国产午夜福利久久久久久| 蜜桃亚洲精品一区二区三区| 亚洲内射少妇av| 国产白丝娇喘喷水9色精品| 麻豆成人午夜福利视频| 热99在线观看视频| 18+在线观看网站| 午夜免费激情av| 亚洲国产精品sss在线观看| 中文字幕久久专区| 一个人观看的视频www高清免费观看| 中文字幕av成人在线电影| 尤物成人国产欧美一区二区三区| 一级二级三级毛片免费看| 天堂中文最新版在线下载 | 亚洲av中文av极速乱| 老师上课跳d突然被开到最大视频| 久久综合国产亚洲精品| 日本黄大片高清| 精品久久久久久成人av| 九九热线精品视视频播放| 国产熟女欧美一区二区| 国产一级毛片在线| 色哟哟·www| 亚洲成人av在线免费| 一级黄片播放器| 亚洲av.av天堂| 内射极品少妇av片p| 简卡轻食公司| 亚洲综合色惰| 特级一级黄色大片| av女优亚洲男人天堂| 99久久精品热视频| 成人永久免费在线观看视频| 麻豆成人午夜福利视频| 久久婷婷人人爽人人干人人爱| av在线亚洲专区| 国产片特级美女逼逼视频| 男女下面进入的视频免费午夜| 欧美高清成人免费视频www| 国产成人影院久久av| 97人妻精品一区二区三区麻豆| 99久久精品一区二区三区| 97热精品久久久久久| 久久99蜜桃精品久久| 精品久久久噜噜| 亚洲五月天丁香| 日韩欧美精品v在线| 国产精品野战在线观看| 一本久久中文字幕| 久久精品国产自在天天线| 中文字幕av在线有码专区| 亚洲最大成人手机在线|