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

    超臨界水工況下鐵素體馬氏體鋼氧化速率常數(shù)的計(jì)算方法與分析

    2018-03-06 01:52:40閻維平
    動(dòng)力工程學(xué)報(bào) 2018年2期
    關(guān)鍵詞:尖晶石磁鐵礦活度

    孫 利, 閻維平

    (華北電力大學(xué) 能源動(dòng)力與機(jī)械工程學(xué)院,河北保定 071003)

    為了提高發(fā)電機(jī)組的效率,國(guó)內(nèi)新建了許多超臨界或超超臨界發(fā)電機(jī)組.在高溫高壓水或超臨界水(SCW)[1]環(huán)境下,材料的抗氧化性是機(jī)組安全運(yùn)行的關(guān)鍵要求之一.作為過熱器管材之一,鐵素體馬氏體鋼具有良好的高溫強(qiáng)度和抗蠕變性、高導(dǎo)熱、低熱膨脹系數(shù)等特性[2],因此受到廣泛關(guān)注.

    目前,關(guān)于超臨界水環(huán)境下鐵素體馬氏體鋼抗氧化性的研究主要集中在實(shí)驗(yàn)方面,采用氧化增重法、X射線衍射(XRD)和掃描電子顯微鏡(SEM)等對(duì)氧化皮的形貌、結(jié)構(gòu)組成和元素分布情況進(jìn)行觀察和分析,推測(cè)出鐵素體馬氏體鋼的氧化機(jī)制.鐘祥玉等[3]采用腐蝕增重法和XRD法研究了400 ℃/25 MPa超臨界水環(huán)境下304不銹鋼和鐵素體馬氏體鋼P(yáng)92(以下簡(jiǎn)寫為P92鋼)的腐蝕行為,分析其腐蝕機(jī)理.尹開鋸等[4]采用XRD、SEM和能譜分析(EDX)技術(shù)研究了500 ℃/25 MPa和550 ℃/25 MPa超臨界水環(huán)境下9Cr-1Mo馬氏體鋼(以下簡(jiǎn)稱T91鋼)和P92鋼氧化膜的多孔性成因.馬云海等[5]研究了4種鍋爐管材料在600 ℃/25 MPa超臨界水環(huán)境下的氧化行為, 發(fā)現(xiàn)均形成了雙層氧化膜,可用楔入機(jī)制模型解釋P91鋼氧化膜的剝落機(jī)制.Ampornrat等[6]在超臨界水環(huán)境下進(jìn)行了鐵素體馬氏體鋼的氧化實(shí)驗(yàn),分析了SCW溫度和溶解氧對(duì)氧化速率的影響.Betova等[7]在500~700 ℃和30 MPa除氧的超臨界水環(huán)境下對(duì)P91鋼進(jìn)行了維持100 h的氧化實(shí)驗(yàn),觀察到氧化物的表面、組分和結(jié)構(gòu),發(fā)現(xiàn)與同一溫度范圍內(nèi)濕空氣氧化的結(jié)果類似.Angell等[8]研究了在500~650 ℃蒸汽環(huán)境下壓力對(duì)9Cr-1Mo鋼氧化的影響,發(fā)現(xiàn)鋼的氧化速率遵循拋物線的動(dòng)力學(xué)規(guī)律,氧化膜厚度隨壓力的升高而增加.Zhang等[9]研究了550 ℃、25 MPa超臨界水環(huán)境下不同溶解氧體積分?jǐn)?shù)對(duì)P92鋼氧化速率的影響,發(fā)現(xiàn)溶解氧體積分?jǐn)?shù)越高,氧化增重越大.Ampornrat等[6]、Bischoff等[10]和Zhong等[11]針對(duì)含9%~12%Cr的鐵素體馬氏體鋼在SCW中的腐蝕行為分別提出了不同的氧化生長(zhǎng)機(jī)制,表明超臨界水環(huán)境下鐵素體馬氏體鋼氧化膜遵循拋物型氧化動(dòng)力學(xué)規(guī)律,氧化皮由內(nèi)、外氧化層和過渡層組成.Rouillard等[12]認(rèn)為550 ℃下9Cr-1Mo鋼在CO2環(huán)境下形成了包含外部磁鐵礦層和內(nèi)部Cr-Fe尖晶石層的雙層氧化皮,空洞誘導(dǎo)氧化機(jī)制解釋了內(nèi)部尖晶石氧化物層的形成.

    目前關(guān)于鐵素體馬氏體鋼氧化的動(dòng)力學(xué)規(guī)律研究著重于實(shí)驗(yàn),模擬計(jì)算相對(duì)較少.筆者著重模擬鐵素體馬氏體鋼在超臨界水環(huán)境下氧化皮生長(zhǎng)的動(dòng)力學(xué)規(guī)律.基于材料界面反應(yīng)的熱力學(xué)分析,通過計(jì)算界面氧活度模擬磁鐵礦層和鐵鉻尖晶石層的增長(zhǎng)速率,將模擬值與文獻(xiàn)[13]中在超臨界水環(huán)境下2種鐵素體馬氏體鋼(HCM12A鋼和NF616鋼)的實(shí)驗(yàn)值進(jìn)行比較.

    1 計(jì)算模型

    不同F(xiàn)e-Cr合金鋼在高溫蒸汽、CO2和液態(tài)LBE環(huán)境下氧化,生成了雙層氧化皮,其氧化特征非常相似,“可用空間”模型可用于解釋形成雙層氧化皮的氧化機(jī)制[14].“可用空間”模型是指內(nèi)外氧化層厚度比隨時(shí)間是恒定的,內(nèi)層的生長(zhǎng)不受氧擴(kuò)散的限制,氧通過短路等微通道到達(dá)氧化物/金屬界面.該機(jī)制中限速步驟是Fe的擴(kuò)散,在氧化物/基體界面形成氧化物生長(zhǎng)的可用空間.由于外層受Fe擴(kuò)散的限制,內(nèi)、外氧化層的形成是相關(guān)的,因此氧化物生長(zhǎng)取決于Fe在內(nèi)外氧化層中的擴(kuò)散.一般來說,合金在超臨界水中氧化,F(xiàn)e向外遷移,在超臨界水/外層界面上生長(zhǎng)形成外層氧化皮,氧向內(nèi)遷移,在金屬/內(nèi)層界面上生長(zhǎng)形成內(nèi)層氧化皮.

    為模擬鐵素體馬氏體鋼在超臨界水環(huán)境下氧化生長(zhǎng)的動(dòng)力學(xué)規(guī)律,提出以下假設(shè):

    (1) 鐵素體馬氏體鋼在超臨界水環(huán)境下生成雙層氧化皮,外層為磁鐵礦層,在超臨界水/磁鐵礦界面生長(zhǎng);內(nèi)層為鐵鉻尖晶石層,在磁鐵礦層/鐵鉻尖晶石層界面生長(zhǎng),形成磁鐵礦層的鐵通量等于形成鐵鉻尖晶石氧化物的鐵通量.2個(gè)氧化層間的界面達(dá)到穩(wěn)定的狀態(tài).

    (2) 整個(gè)氧化層存在局部熱力學(xué)平衡.

    (3) 內(nèi)層和外層界面組分質(zhì)量濃度恒定.

    (4) 氧化物生長(zhǎng)遵循拋物線規(guī)律.

    (5) 界面處的氧活度保持恒定.

    (6) 鐵鉻尖晶石層組分偏離化學(xué)計(jì)量比很小.

    (7) 不考慮空洞、裂紋等的影響.

    (8) 不考慮氧化物的溶解.

    一般來說,消耗的T91鋼厚度等于生成的鐵鉻尖晶石層厚度,鐵鉻尖晶石層厚度與磁鐵礦層厚度的關(guān)系為[14]:

    (1)

    式中:hsp、hmag分別為鐵鉻尖晶石層和磁鐵礦層的厚度,mm;cFe,mag、cFe,a和cFe,sp分別為磁鐵礦、鐵素體馬氏鋼和鐵鉻尖晶石中Fe的物質(zhì)的量濃度,mol/cm3.

    對(duì)于給定的氧化時(shí)間,氧化皮厚度為:

    h2=kpt

    (2)

    式中:h為氧化皮厚度,mm;kp為氧化速率常數(shù);t為氧化時(shí)間.

    根據(jù)Wagner理論,鐵離子向外擴(kuò)散并在合金表面形成Fe3O4,氧化速率常數(shù)為[15]:

    (3)

    (4)

    Fe的擴(kuò)散系數(shù)大小與氧化物晶格中存在的結(jié)構(gòu)缺陷有關(guān),而缺陷濃度與氧活度有關(guān),因此Fe的擴(kuò)散系數(shù)取決于氧化物中的氧活度.磁鐵礦中Fe的擴(kuò)散系數(shù)DFe,mag可用氧活度和溫度的函數(shù)表示[16]:

    (5)

    A=DVKV/12

    (6)

    B=4DIKI/3

    (7)

    式中:DV、DI分別為磁鐵礦層中陽(yáng)離子的空位擴(kuò)散系數(shù)和間隙擴(kuò)散系數(shù);KV、KI分別為空位擴(kuò)散和間隙擴(kuò)散的常數(shù),均為溫度的函數(shù)[16].

    式(5)中第1項(xiàng)代表Fe的空位擴(kuò)散,第2項(xiàng)代表Fe的間隙擴(kuò)散,均為溫度和氧活度的函數(shù).對(duì)于低氧活度,F(xiàn)e優(yōu)先發(fā)生間隙擴(kuò)散;對(duì)于高氧活度,F(xiàn)e優(yōu)先發(fā)生空位擴(kuò)散.

    鐵鉻尖晶石中,Fe的擴(kuò)散系數(shù)可用氧活度和溫度的函數(shù)表示[12]:

    (8)

    RV=DTopfer,sp,V/DTopfer,mag,V

    (9)

    RI=DTopfer,sp,I/DTopfer,mag,I

    (10)

    式中:RV為根據(jù)T?pfer等[17]的實(shí)驗(yàn)數(shù)據(jù)外推得到的鐵鉻尖晶石層與磁鐵礦層中陽(yáng)離子空位擴(kuò)散系數(shù)的比值;RI為鐵鉻尖晶石層與磁鐵礦層中陽(yáng)離子間隙擴(kuò)散系數(shù)的比值;DTopfer,sp,V為鐵鉻尖晶石層中Fe的空位擴(kuò)散系數(shù);DTopfer,mag,V為磁鐵礦層中Fe的空位擴(kuò)散系數(shù);DTopfer,sp,I為鐵鉻尖晶石層中Fe的間隙擴(kuò)散系數(shù);DTopfer,mag,I為磁鐵礦層中Fe的間隙擴(kuò)散系數(shù).

    根據(jù)假設(shè)(1), 鐵鉻尖晶石層內(nèi)Fe的擴(kuò)散通量等于磁鐵礦內(nèi)Fe的擴(kuò)散通量.鐵素體馬氏體鋼/鐵鉻尖晶石界面的氧活度由不同溫度下相圖的相邊界氧分壓確定,超臨界水/磁鐵礦界面的氧活度由熱力學(xué)平衡條件確定,鐵鉻尖晶石/磁鐵礦界面的氧活度是未知的.根據(jù)式(1)~式(10),可得到磁鐵礦層和鐵鉻尖晶石層的氧化速率常數(shù)、各層厚度以及鐵鉻尖晶石/磁鐵礦界面的氧活度.

    2 參數(shù)的確定

    選取HCM12A鋼和NF616鋼為研究對(duì)象,分別模擬400 ℃/25 MPa、500 ℃/25 MPa和600 ℃/25 MPa超臨界水環(huán)境下的氧化速率常數(shù).

    2.1 材料各成分的質(zhì)量分?jǐn)?shù)

    HCM12A鋼和NF616鋼各成分的質(zhì)量分?jǐn)?shù)如表1所示[13].

    表1 HCM12A鋼和NF616鋼各成分的質(zhì)量分?jǐn)?shù)

    2.2 Fe在磁鐵礦層和鐵鉻尖晶石層的擴(kuò)散系數(shù)

    Backhaus等[16]給出了溫度范圍為900~1 400 ℃時(shí)Fe在磁鐵礦中擴(kuò)散系數(shù)的實(shí)驗(yàn)數(shù)據(jù),可推算較低溫度時(shí)的擴(kuò)散系數(shù).T?pfer等[17]給出了溫度為1 200 ℃時(shí)Fe在不同鐵鉻尖晶石層擴(kuò)散系數(shù)的實(shí)驗(yàn)數(shù)據(jù),可推算出較低溫度時(shí)的擴(kuò)散系數(shù).

    在同一溫度下,F(xiàn)e在磁鐵礦和鐵鉻尖晶石中的擴(kuò)散系數(shù)具有相關(guān)性,即Fe在磁鐵礦和鐵鉻尖晶石中的擴(kuò)散系數(shù)的比值為定值.分區(qū)域(空位機(jī)制或間隙機(jī)制)求各個(gè)實(shí)驗(yàn)點(diǎn)Fe的空位擴(kuò)散系數(shù)與間隙鐵擴(kuò)散系數(shù)的比值,進(jìn)而求算術(shù)平均值,即可得到RV和RI值.NF616鋼內(nèi)層氧化物的組分為Fe2.4Cr0.6O4,RV值為0.644,RI值為0.31;HCM12A鋼內(nèi)層氧化物的組分為Fe2CrO4,RV值為0.5,RI值為0.039.

    2.3 超臨界水/磁鐵礦界面和合金基體/鐵鉻尖晶石界面的氧活度

    根據(jù)合金組成和氧分壓可確定生成氧化物的類型,由Fe-Cr-O優(yōu)勢(shì)圖[17]可查到400 ℃、500 ℃和600 ℃ 時(shí)鐵素體馬氏體鋼/鐵鉻尖晶石界面的氧分壓,即相邊界氧分壓,分別為3×10-30Pa、10-24Pa和2×10-20Pa.

    Fe在鐵鉻尖晶石層擴(kuò)散,與H2O反應(yīng),且隨著氧化時(shí)間的延長(zhǎng),F(xiàn)e3O4層逐漸變厚.在磁鐵礦/超臨界水界面,磁鐵礦的晶界通過反應(yīng)(11)形成了Fe2O3.氧化物沿晶界和短路通道向內(nèi)傳輸,F(xiàn)e3O4層和鐵鉻尖晶石層變厚.磁鐵礦/超臨界水界面處參與反應(yīng)的氧分壓值不是溶解氧的分壓值,而是滿足反應(yīng)平衡的等效氧分壓值.

    (11)

    利用Gaskell[20]提供的數(shù)據(jù),反應(yīng)的標(biāo)準(zhǔn)吉布斯自由能可表示為:

    ΔGT=-3 440+81.7T

    (12)

    式中:ΔGT為吉布斯自由能,J;T為溫度,K.

    溫度為600 ℃時(shí),ΔGT值為67 884 J,平衡時(shí)蒸汽與氫氣分壓的比值pH2O/pH2為11 532,由水的離解方程、吉布斯自由能和水離解反應(yīng)速率可得到600 ℃時(shí)磁鐵礦/超臨界水界面的等效氧分壓.不同溫度下磁鐵礦/超臨界水界面處的等效氧分壓如表2所示.

    2.4 氧化皮增重與氧化皮內(nèi)外層厚度的關(guān)系

    氧的吸收導(dǎo)致氧化皮增重,且氧化皮增重與氧化皮厚度成線性關(guān)系,氧化皮的平均密度為定值,即隨氧化時(shí)間的增加,內(nèi)外層氧化皮厚度增加,但氧化皮平均密度保持不變.氧化皮增重與氧化皮內(nèi)外層厚度的關(guān)系可表示為[9]:

    (13)

    式中:Δw為單位面積吸氧量,mg/cm2;ρouter和ρinner分別為外層和內(nèi)層氧化皮的密度,g/cm3;houter和hinner分別為外層和內(nèi)層氧化物的厚度,cm;Mo,outer和Mo,inner分別為外層和內(nèi)層氧化物氧的摩爾質(zhì)量,g/mol;MFe3O4和MFe3O4,FeCr2O4分別為外層和內(nèi)層氧化物的摩爾質(zhì)量,g/mol.

    表2 熱力學(xué)數(shù)據(jù)

    3 模擬結(jié)果及討論

    3.1 Fe在磁鐵礦和鐵鉻尖晶石中的擴(kuò)散系數(shù)

    利用式(5)~式(10)可求出不同溫度(400 ℃、500 ℃和600 ℃)下Fe在磁鐵礦和鐵鉻尖晶石中的擴(kuò)散系數(shù).以氧活度的對(duì)數(shù)為橫坐標(biāo),F(xiàn)e擴(kuò)散系數(shù)的對(duì)數(shù)為縱坐標(biāo),氧活度與Fe擴(kuò)散系數(shù)的關(guān)系如圖1和圖2所示.圖1中標(biāo)出了600 ℃下Fe在磁鐵礦中空位擴(kuò)散和間隙擴(kuò)散區(qū)域.

    圖1 不同溫度下Fe在HAM12A鋼雙層氧化皮內(nèi)擴(kuò)散系數(shù)與氧活度的關(guān)系

    Fig.1 Iron diffusion coefficients in duplex scale of HAM12A vs. oxygen activity at different temperatures

    氧化生長(zhǎng)速率與界面處的氧活度密切相關(guān),F(xiàn)e的擴(kuò)散系數(shù)與氧化層中氧活度有較高的相關(guān)性.氧活度的值會(huì)改變Fe的擴(kuò)散機(jī)制,高氧活度遵循空位機(jī)制,低氧活度遵循間隙機(jī)制.鐵素體馬氏體鋼/氧化物界面處氧活度低,F(xiàn)e擴(kuò)散遵循間隙機(jī)制.

    圖2 500 ℃和600 ℃下Fe在NF616鋼雙層氧化皮內(nèi)擴(kuò)散系數(shù)和氧活度的關(guān)系

    Fig.2 Iron diffusion coefficients in duplex scale of NF616 vs. oxygen activity at 500 ℃ and 600 ℃

    在鐵鉻尖晶石/磁鐵礦界面,F(xiàn)e的擴(kuò)散系數(shù)不連續(xù).由于鐵鉻尖晶石/磁鐵礦界面的氧活度需經(jīng)過計(jì)算得出,所以無法預(yù)測(cè)該界面的擴(kuò)散機(jī)制.磁鐵礦/超臨界水界面處的氧活度較高,通過空位機(jī)制進(jìn)行擴(kuò)散.

    3.2 模擬和實(shí)驗(yàn)的氧化增重及其誤差分析

    根據(jù)式(1)~式(10)、式(13)可模擬400 ℃/25 MPa、500 ℃/25 MPa和600 ℃/25 MPa環(huán)境下HCM12A鋼和NF616鋼內(nèi)外層氧化皮的氧化速率及雙層氧化皮界面的氧分壓.模擬與實(shí)驗(yàn)的增重曲線如圖3和圖4所示.

    圖3 不同溫度下HCM12A鋼的模擬與實(shí)驗(yàn)增重曲線

    Fig.3 Simulated and experimental weight gain curves of HCM12A at different temperatures

    文獻(xiàn)[22]中的實(shí)驗(yàn)是在壓力為25 MPa,溫度分別為400 ℃、500 ℃和 600 ℃,溶解氧體積分?jǐn)?shù)小于2.5×10-8的環(huán)境下進(jìn)行的,氧化時(shí)間為1~3 000 h.

    對(duì)于HCM12A鋼,氧化速率的模擬值均低于實(shí)驗(yàn)值,600 ℃時(shí)模擬與實(shí)驗(yàn)的氧化增重最大差值為104 mg/dm2,折合成厚度約為7.5 μm;500 ℃時(shí)模擬與實(shí)驗(yàn)的氧化增重最大差值為181 mg/dm2,折合成厚度約為13 μm;400 ℃時(shí)模擬與實(shí)驗(yàn)的氧化增重最大差值為2 mg/dm2,折合成厚度約為1.57 μm.對(duì)于NF616鋼,600 ℃時(shí)氧化速率的模擬值高于實(shí)驗(yàn)值,模擬與實(shí)驗(yàn)的氧化增重最大差值為306 mg/dm2,折合成厚度約為22 μm,原因可能是在該溫度下外氧化層不含F(xiàn)e2O3[9],模擬過程中選取的平衡氧分壓比實(shí)際值大;500 ℃時(shí)氧化速率的模擬值低于實(shí)驗(yàn)值,模擬與實(shí)驗(yàn)的氧化增重最大差值為76 mg/dm2,折合成厚度約為5.4 μm.氧化3 000 h后,不同溫度下HCM12A鋼和NF616鋼氧化皮增重模擬值與實(shí)驗(yàn)值的相對(duì)誤差如圖5和圖6所示.

    圖4 500 ℃、600 ℃下NF616鋼氧化皮的模擬與實(shí)驗(yàn)增重曲線

    Fig.4 Simulated and experimental weight gain curves of NF616 at 500 ℃ and 600 ℃

    圖5 不同溫度下HCM12A鋼氧化皮增重模擬值與實(shí)驗(yàn)值的相對(duì)誤差

    Fig.5 Relative error between simulated and experimental weight gain of HCM12A at different temperatures

    NF616鋼(500 ℃和600 ℃)和HCM12A鋼(600 ℃)的氧化皮增重模擬值與實(shí)驗(yàn)值的相對(duì)誤差基本小于20%,而400 ℃和500 ℃時(shí)HCM12A氧化皮增重模擬值與實(shí)驗(yàn)值的相對(duì)誤差大于40%.

    磁鐵礦層和鐵鉻尖晶石層的厚度遵循拋物線增長(zhǎng)規(guī)律,根據(jù)式(1)和式(13)可將氧化皮厚度轉(zhuǎn)換為氧化增重,模擬的氧化增重與t0.5成一定比例.由圖3可知,500 ℃時(shí)HCM12A鋼的實(shí)驗(yàn)氧化增重時(shí)間指數(shù)n為0.497[22],接近0.5.因此相對(duì)誤差為常量,如圖5所示.由圖5可知,400 ℃和600 ℃時(shí)HCM12A鋼的實(shí)驗(yàn)氧化增重時(shí)間指數(shù)n分別為0.402和0.418[22],模擬相對(duì)誤差逐漸減小.

    圖6 500 ℃、600 ℃下NF616鋼氧化皮增重模擬值與實(shí)驗(yàn)值的相對(duì)誤差

    Fig.6 Relative error between simulated and experimental weight gain of NF616 at 500 ℃ and 600 ℃

    3.3 氧化速率常數(shù)的比較與分析

    表3給出了模擬得到的界面氧分壓、內(nèi)外層氧化皮的氧化速率(kmag、ksp)和總氧化速率(kox)以及美國(guó)電力研究院(EPRI)的相關(guān)數(shù)據(jù)[23].

    由表3可知,NF616鋼(500 ℃和600 ℃)和HCM12A鋼(600 ℃)的氧化速率與EPRI的氧化速率在同一數(shù)量級(jí).400 ℃時(shí)HCM12A鋼的氧化速率比EPRI的氧化速率低2個(gè)數(shù)量級(jí),500 ℃時(shí)HCM12A鋼的氧化速率比EPRI的氧化速率低1個(gè)數(shù)量級(jí),其原因可能是將高溫?cái)U(kuò)散系數(shù)的數(shù)據(jù)推算至低溫時(shí)產(chǎn)生了較大誤差.600 ℃時(shí)NF616鋼的氧化速率比500 ℃時(shí)近似高2個(gè)數(shù)量級(jí).600 ℃時(shí)HCM12A鋼的氧化速率比500 ℃時(shí)高1個(gè)數(shù)量級(jí).由于存在內(nèi)層鉻,鐵鉻尖晶石層的氧化速率低于磁鐵礦的氧化速率.相同溫度下,由于NF616鋼的Cr含量低于HCM12A鋼的Cr含量,因此NF616鋼的氧化速率高于HCM12A鋼的氧化速率.NF616鋼在600 ℃下氧化時(shí),鐵鉻尖晶石層的氧活度從2×10-25變化至1.41×10-23,其擴(kuò)散過程遵循間隙機(jī)制;磁鐵礦層的氧活度從1.41×10-23變化至2.57×10-16,其擴(kuò)散過程由間隙機(jī)制轉(zhuǎn)變?yōu)榭瘴粰C(jī)制.HCM12A鋼在600 ℃下氧化時(shí),鐵鉻尖晶石層的氧活度從2×10-25變化至9.59×10-17,其擴(kuò)散過程由間隙機(jī)制轉(zhuǎn)變?yōu)榭瘴粰C(jī)制;磁鐵礦層的氧活度從9.59×10-17變化至2.57×10-16,擴(kuò)散過程遵循空位機(jī)制.

    表3 內(nèi)外層氧化皮的氧化速率(kmag、ksp)、總氧化速率(kox)以及界面氧分壓模擬結(jié)果

    4 結(jié) 論

    (1)推算出較低溫下Fe在磁鐵礦和鐵鉻尖晶石層中的擴(kuò)散系數(shù),可用于模擬超臨界水條件下氧化的動(dòng)力學(xué)規(guī)律.

    (2)NF616 鋼(500 ℃、600 ℃)和HCM12A鋼(600 ℃)氧化皮增重的模擬值與實(shí)驗(yàn)值較接近,說明基本假設(shè)是有效的.

    (3)氧化皮的氧化速率模擬值取決于合金/內(nèi)氧化層界面和磁鐵礦/超臨界水界面氧分壓的取值,用反應(yīng)(11)的平衡氧分壓作為磁鐵礦/超臨界水界面的氧分壓,氧化速率常數(shù)的模擬值與實(shí)驗(yàn)值接近.

    (4)在400 ℃、500 ℃條件下HCM12A鋼氧化皮增重的模擬值與實(shí)驗(yàn)值誤差為40%,氧化速率比EPRI的氧化速率低1~2個(gè)數(shù)量級(jí),原因可能是磁鐵礦/超臨界水界面的氧分壓不是反應(yīng)(11)的平衡氧分壓,或是根據(jù)T?pfer實(shí)驗(yàn)數(shù)據(jù)外推到低溫的擴(kuò)散系數(shù)有較大誤差.

    [1] KRITZER P. Corrosion in high-temperature and supercritical water and aqueous solutions: a review[J].TheJournalofSupercriticalFluids, 2004, 29(1/2): 1-29.

    [2] KLUEH R L, NELSON A T. Ferritic/martensitic steels for next-generation reactors[J].JournalofNuclearMaterials, 2007, 371(1/3): 37-52.

    [3] 鐘祥玉, 吳欣強(qiáng), 韓恩厚. 核級(jí)不銹鋼和鐵素體-馬氏體耐熱鋼在400 ℃/25 MPa超臨界水中的腐蝕行為[J].金屬學(xué)報(bào), 2011, 47(7): 932-938.

    ZHONG Xiangyu, WU Xinqiang, HAN Enhou. Corrosion behaviors of nuclear-grade stainless steel and ferritic-martensitic steel in supercritical water[J].ActaMetallurgicaSinica, 2011, 47(7): 932-938.

    [4] 尹開鋸, 邱紹宇, 唐睿, 等. 鐵素體-馬氏體鋼P(yáng)91和P92在超臨界水中腐蝕后氧化膜多孔性分析[J].中國(guó)腐蝕與防護(hù)學(xué)報(bào),2010, 30(1): 1-5.

    YIN Kaiju, QIU Shaoyu, TANG Rui, et al. Characterization of the porosity of the oxide scales on ferritic-martensitic steel P91 and P92 exposed in supercritical water[J].JournalofChineseSocietyforCorrosionandProtection, 2010, 30(1): 1-5.

    [5] 馬云海, 鄭開云. 四種鍋爐管材料650 ℃/26 MPa的蒸汽氧化研究[J].動(dòng)力工程學(xué)報(bào), 2015, 35(6): 509-516.

    MA Yunhai, ZHENG Kaiyun. Study on steam oxidation of 4 boiler tube materials at 650 ℃/26 MPa[J].JournalofChineseSocietyofPowerEngineering, 2015, 35(6): 509-516.

    [6] AMPORNRAT P, WAS G S. Oxidation of ferritic-martensitic alloys T91, HCM12A and HT-9 in supercritical water[J].JournalofNuclearMaterials, 2007, 371(1/3): 1-17.

    [7] BETOVA I, BOJINOV M, KINNUNEN P, et al. Composition, structure, and properties of corrosion layers on ferritic and austenitic steels in ultra super critical water[J].JournaloftheElectrochemicalSociety, 2006, 153(11): B464-B473.

    [8] ANGELL M G, LISTER S K, RUDGE A. The effect of steam pressure on the oxidation behavior of annealed 9Cr 1Mo boiler tubing materials[C]//15thInternationalConferenceonthePropertiesofWaterandSteam(ICPWSXV).Berlin, Germany: [s.n.], 2008.

    [9] ZHANG Naiqiang, XU Hong, LI Baorong, et al. Influence of the dissolved oxygen content on corrosion of the ferritic-martensitic steel P92 in supercritical water[J].CorrosionScience, 2012, 56: 123-128.

    [10] BISCHOFF J, MOTTA A T. Oxidation behavior of ferritic-martensitic and ODS steels in supercritical water[J].JournalofNuclearMaterials, 2012, 424(1/3): 261-276.

    [11] ZHONG Xiangyu, WU Xinqiang, HAN Enhou. Effects of exposure temperature and time on corrosion behavior of a ferritic-martensitic steel P92 in aerated supercritical water[J].CorrosionScience, 2015, 90: 511-521.

    [12] ROUILLARD F, MARTINELLI L. Corrosion of 9Cr steel in CO2at intermediate temperature III: modelling and simulation of void-induced duplex oxide growth[J].OxidationofMetals, 2012, 77(1/2): 71-83.

    [13] TAN L, REN X, ALLEN T R. Corrosion behavior of 9-12% Cr ferritic-martensitic steels in supercritical water[J].CorrosionScience, 2010, 52(4): 1520-1528.

    [14] MARTINELLI L, BALBAUD-CéLéRIER F, TERLAIN A, et al. Oxidation mechanism of an Fe-9Cr-1Mo steel by liquid Pb-Bi eutectic alloy at 470 ℃ (part II)[J].CorrosionScience, 2008, 50(9): 2537-2548.

    [15] ATKINSON A, O'DWYER M L, TAYLOR R I.55Fe diffusion in magnetite crystals at 500 ℃ and its relevance to oxidation of iron[J].JournalofMaterialsScience, 1983, 18(8): 2371-2379.

    [16] BACKHAUS R M, DIECKMANN R. Defects and cation diffusion in magnetite (VII): diffusion controlled formation of magnetite during reactions in the iron-oxygen system[J].BerichteDerBunsengesellschaftFürPhysikalischeChemie, 1986, 90(8): 690-698.

    [17] T?PFER J, AGGARWAL S, DIECKMANN R. Point defects and cation tracer diffusion in (CrxFe1-x)3-δO4spinels[J].SolidStateIonics, 1995, 81(3/4): 251-266.

    [19] FURUKAWA T, MüLLER G, SCHUMACHER G, et al. Effect of oxygen concentration and temperature on compatibility of ODS steel with liquid, stagnant Pb45Bi55[J].JournalofNuclearMaterials, 2004, 335(2): 189-193.

    [20] GASKELL D R. Introduction to metallurgical thermo dynamics [M].2nd ed. London: Taylor & Francis, 1981.

    [21] TAN L, REN X, ALLEN T R. Corrosion behavior of 9-12% Cr ferritic-martensiticsteels in supercritical water[J].CorrosionScience, 2010,52(4): 1520-1528.

    [22] AMPORNRAT P. Determination of oxidation mechanisms of ferritic-martensitic alloys in supercritical water[D]. Michigan, USA: The University of Michigan, 2011.

    [23] DOOLEY R B, WRIGHT I G, TOETORELLI P F, et al. Program on technology innovation: oxide growth and exfoliation on alloys exposed to steam[R]. California, USA: Electric Power Research Institute, 2007.

    猜你喜歡
    尖晶石磁鐵礦活度
    HISMELT SRV環(huán)境下剛玉尖晶石材料抗侵蝕性能研究
    山東冶金(2022年4期)2022-09-14 08:58:10
    鎂鋁尖晶石種類對(duì)尖晶石-方鎂石復(fù)相材料燒結(jié)性能的影響
    耐火材料(2022年4期)2022-08-28 03:01:10
    尖晶石的資源與商貿(mào)現(xiàn)狀
    釩鈦磁鐵礦冶煉工藝比較分析
    山東冶金(2019年1期)2019-03-30 01:35:06
    二次鋁灰燒結(jié)制備鈣鋁黃長(zhǎng)石/鎂鋁尖晶石復(fù)相材料
    CaO-SiO2-FeO-P2O5-Al2O3脫磷渣系中組元活度的計(jì)算
    釩鈦磁鐵礦球團(tuán)氧化焙燒行為和固結(jié)特性
    核電廠惰性氣體排放活度濃度的估算
    釩鈦磁鐵礦行業(yè)標(biāo)準(zhǔn)將制定
    鉆井液處理劑溶液活度測(cè)量方法對(duì)比
    亚洲国产日韩一区二区| 国语对白做爰xxxⅹ性视频网站| 国产精品一二三区在线看| 丝袜喷水一区| 精品国产一区二区三区久久久樱花| 美女主播在线视频| 天美传媒精品一区二区| 婷婷成人精品国产| 亚洲国产日韩一区二区| 如日韩欧美国产精品一区二区三区| www.熟女人妻精品国产 | 熟女av电影| 成人毛片60女人毛片免费| 久久这里只有精品19| 国产综合精华液| a级毛片黄视频| 免费黄频网站在线观看国产| 国产高清不卡午夜福利| 国产欧美日韩综合在线一区二区| 九色成人免费人妻av| 亚洲精品国产av成人精品| 欧美成人午夜免费资源| 老司机影院毛片| 欧美bdsm另类| av国产久精品久网站免费入址| 狠狠精品人妻久久久久久综合| 91精品三级在线观看| 2022亚洲国产成人精品| 综合色丁香网| 午夜精品国产一区二区电影| 国产男女超爽视频在线观看| 亚洲成国产人片在线观看| 成人无遮挡网站| 久久人人97超碰香蕉20202| 不卡视频在线观看欧美| 久久久久精品久久久久真实原创| 精品少妇内射三级| 久久久a久久爽久久v久久| 亚洲成人av在线免费| 亚洲美女视频黄频| 大片免费播放器 马上看| 国产爽快片一区二区三区| 久久鲁丝午夜福利片| 国产精品无大码| 黑人欧美特级aaaaaa片| 妹子高潮喷水视频| 欧美激情国产日韩精品一区| 国产极品粉嫩免费观看在线| 精品久久久精品久久久| 国产淫语在线视频| av在线观看视频网站免费| www.熟女人妻精品国产 | 波多野结衣一区麻豆| 人人澡人人妻人| 久久久精品区二区三区| 五月开心婷婷网| 国产一区二区在线观看日韩| av国产久精品久网站免费入址| 国产乱人偷精品视频| 免费高清在线观看日韩| 国产成人91sexporn| 欧美日韩国产mv在线观看视频| 国产乱人偷精品视频| 熟女电影av网| 国产精品欧美亚洲77777| 97人妻天天添夜夜摸| 日韩欧美一区视频在线观看| 国产永久视频网站| 成人漫画全彩无遮挡| 久久99热6这里只有精品| 熟女av电影| 两性夫妻黄色片 | 久久精品aⅴ一区二区三区四区 | 亚洲综合色惰| 国语对白做爰xxxⅹ性视频网站| 国产亚洲精品第一综合不卡 | 欧美日韩av久久| 亚洲国产看品久久| 日韩av不卡免费在线播放| 免费av中文字幕在线| 日韩 亚洲 欧美在线| 午夜免费男女啪啪视频观看| 日本与韩国留学比较| 日本爱情动作片www.在线观看| 黄色 视频免费看| 久久鲁丝午夜福利片| 亚洲欧洲国产日韩| 国内精品宾馆在线| 亚洲精华国产精华液的使用体验| 欧美亚洲 丝袜 人妻 在线| 两个人看的免费小视频| 老熟女久久久| 高清黄色对白视频在线免费看| 国产国拍精品亚洲av在线观看| 国产 精品1| 热99久久久久精品小说推荐| 日韩伦理黄色片| 久久午夜福利片| 涩涩av久久男人的天堂| 母亲3免费完整高清在线观看 | 91精品伊人久久大香线蕉| 日韩制服丝袜自拍偷拍| 天天操日日干夜夜撸| 丰满乱子伦码专区| 三上悠亚av全集在线观看| 免费av不卡在线播放| 99热网站在线观看| 51国产日韩欧美| 日韩成人伦理影院| 日韩中字成人| av不卡在线播放| 秋霞伦理黄片| 欧美日韩精品成人综合77777| 在线观看免费高清a一片| 亚洲精品一区蜜桃| 日韩av不卡免费在线播放| 丰满迷人的少妇在线观看| 女的被弄到高潮叫床怎么办| 黑丝袜美女国产一区| 精品久久蜜臀av无| 丰满饥渴人妻一区二区三| 亚洲精品美女久久av网站| 日韩制服骚丝袜av| 亚洲人与动物交配视频| 国产毛片在线视频| 精品一区二区三区四区五区乱码 | 日韩熟女老妇一区二区性免费视频| 国产麻豆69| 亚洲av免费高清在线观看| 80岁老熟妇乱子伦牲交| 国产一区二区三区av在线| 91精品国产国语对白视频| 国产亚洲精品第一综合不卡 | 如何舔出高潮| 亚洲成色77777| 亚洲精品乱久久久久久| 少妇人妻精品综合一区二区| 夫妻性生交免费视频一级片| 男女边吃奶边做爰视频| 自线自在国产av| 日本免费在线观看一区| 国产精品一区二区在线不卡| 国产亚洲精品第一综合不卡 | 亚洲国产成人一精品久久久| 午夜福利网站1000一区二区三区| 看免费av毛片| 国产一区有黄有色的免费视频| 日韩一本色道免费dvd| 青青草视频在线视频观看| 久久人人97超碰香蕉20202| 日本-黄色视频高清免费观看| 91久久精品国产一区二区三区| 丝袜在线中文字幕| 免费观看在线日韩| 久久青草综合色| 在线观看人妻少妇| 国产精品成人在线| 亚洲av免费高清在线观看| 熟女av电影| 久久这里只有精品19| 亚洲国产色片| 性高湖久久久久久久久免费观看| 精品久久国产蜜桃| 国产精品一二三区在线看| 成年动漫av网址| 国产不卡av网站在线观看| 黑人欧美特级aaaaaa片| 卡戴珊不雅视频在线播放| 亚洲精品中文字幕在线视频| 亚洲,欧美,日韩| 日韩av不卡免费在线播放| 欧美精品国产亚洲| 综合色丁香网| 嫩草影院入口| 在线观看www视频免费| 亚洲国产毛片av蜜桃av| 中国美白少妇内射xxxbb| 日韩av不卡免费在线播放| 女人精品久久久久毛片| 欧美日韩视频高清一区二区三区二| 亚洲色图综合在线观看| 国语对白做爰xxxⅹ性视频网站| 蜜桃国产av成人99| 亚洲国产毛片av蜜桃av| 午夜福利在线观看免费完整高清在| 婷婷色综合大香蕉| 有码 亚洲区| 男男h啪啪无遮挡| 水蜜桃什么品种好| 国产乱来视频区| 丝瓜视频免费看黄片| 久久狼人影院| 波野结衣二区三区在线| 国产精品成人在线| 一级毛片电影观看| 国产综合精华液| 丝袜美足系列| 欧美日韩亚洲高清精品| 亚洲少妇的诱惑av| 亚洲中文av在线| 在线观看免费日韩欧美大片| 国产欧美日韩综合在线一区二区| 观看美女的网站| 王馨瑶露胸无遮挡在线观看| 免费日韩欧美在线观看| 在线亚洲精品国产二区图片欧美| 成年女人在线观看亚洲视频| 侵犯人妻中文字幕一二三四区| videos熟女内射| 欧美激情极品国产一区二区三区 | 观看美女的网站| 99热6这里只有精品| 亚洲国产精品国产精品| 久久av网站| 成人亚洲精品一区在线观看| 插逼视频在线观看| 七月丁香在线播放| 久久97久久精品| 三级国产精品片| 青春草视频在线免费观看| 在线天堂中文资源库| 亚洲国产精品一区二区三区在线| 精品国产一区二区久久| 中国国产av一级| www日本在线高清视频| 亚洲伊人色综图| 大香蕉久久网| 永久网站在线| 26uuu在线亚洲综合色| 91在线精品国自产拍蜜月| 妹子高潮喷水视频| 国产午夜精品一二区理论片| 26uuu在线亚洲综合色| 国产免费一区二区三区四区乱码| 丁香六月天网| 日韩人妻精品一区2区三区| 久久人人爽av亚洲精品天堂| 欧美成人午夜精品| 国产在线一区二区三区精| h视频一区二区三区| 另类亚洲欧美激情| 免费播放大片免费观看视频在线观看| 美国免费a级毛片| 亚洲人成网站在线观看播放| 90打野战视频偷拍视频| av网站免费在线观看视频| av.在线天堂| xxx大片免费视频| 日韩成人av中文字幕在线观看| 啦啦啦中文免费视频观看日本| 亚洲av综合色区一区| 欧美精品亚洲一区二区| 好男人视频免费观看在线| 国产69精品久久久久777片| 伊人久久国产一区二区| 精品国产一区二区三区四区第35| 国产精品一区二区在线观看99| 国产精品无大码| 亚洲精品国产色婷婷电影| 水蜜桃什么品种好| 赤兔流量卡办理| 男女边摸边吃奶| 丰满迷人的少妇在线观看| 国产精品国产三级国产专区5o| 久久久精品免费免费高清| 亚洲国产精品一区二区三区在线| 国产一区二区在线观看av| 啦啦啦中文免费视频观看日本| 黄网站色视频无遮挡免费观看| 色5月婷婷丁香| 91aial.com中文字幕在线观看| 久久青草综合色| 国产爽快片一区二区三区| freevideosex欧美| 丝袜脚勾引网站| 国产男女超爽视频在线观看| av卡一久久| 欧美人与善性xxx| 少妇的逼好多水| 涩涩av久久男人的天堂| www.色视频.com| 飞空精品影院首页| 成人黄色视频免费在线看| 午夜激情av网站| 免费观看av网站的网址| 日韩欧美精品免费久久| 成人亚洲精品一区在线观看| 观看美女的网站| 久久毛片免费看一区二区三区| 精品午夜福利在线看| 国产在线一区二区三区精| 日本色播在线视频| 亚洲精品中文字幕在线视频| 香蕉国产在线看| 51国产日韩欧美| 在线看a的网站| 卡戴珊不雅视频在线播放| 人人妻人人爽人人添夜夜欢视频| 中国三级夫妇交换| 天堂俺去俺来也www色官网| 日本-黄色视频高清免费观看| 日韩三级伦理在线观看| 王馨瑶露胸无遮挡在线观看| 色婷婷av一区二区三区视频| 王馨瑶露胸无遮挡在线观看| videos熟女内射| 女性生殖器流出的白浆| 日本猛色少妇xxxxx猛交久久| 欧美亚洲日本最大视频资源| 肉色欧美久久久久久久蜜桃| 有码 亚洲区| 三上悠亚av全集在线观看| 国产成人精品无人区| 丰满迷人的少妇在线观看| 美女视频免费永久观看网站| 99国产精品免费福利视频| 日本av免费视频播放| videossex国产| 久久这里有精品视频免费| 欧美日韩精品成人综合77777| 9191精品国产免费久久| 欧美精品高潮呻吟av久久| 国产女主播在线喷水免费视频网站| 免费观看无遮挡的男女| 日韩精品免费视频一区二区三区 | av网站免费在线观看视频| av一本久久久久| a 毛片基地| 久久久久久久久久成人| 国产熟女午夜一区二区三区| 99久国产av精品国产电影| 伦精品一区二区三区| 亚洲天堂av无毛| 色视频在线一区二区三区| 十八禁高潮呻吟视频| 亚洲av欧美aⅴ国产| 搡老乐熟女国产| 国产又爽黄色视频| 国产精品一国产av| 99久久综合免费| 亚洲少妇的诱惑av| 午夜av观看不卡| 国内精品宾馆在线| 午夜日本视频在线| 日本猛色少妇xxxxx猛交久久| 性色av一级| 少妇高潮的动态图| 亚洲av男天堂| 丝袜脚勾引网站| 亚洲成av片中文字幕在线观看 | 国产男女超爽视频在线观看| 一边亲一边摸免费视频| 999精品在线视频| 国产熟女欧美一区二区| 亚洲国产精品一区三区| 成年av动漫网址| 永久免费av网站大全| 久久久久人妻精品一区果冻| av线在线观看网站| 国产精品偷伦视频观看了| 国产精品久久久久成人av| 色5月婷婷丁香| 人人妻人人澡人人看| 美女国产高潮福利片在线看| 国产成人午夜福利电影在线观看| 中文字幕人妻熟女乱码| 亚洲综合精品二区| 宅男免费午夜| 成人国产麻豆网| 青春草国产在线视频| 亚洲国产精品一区二区三区在线| 精品少妇久久久久久888优播| 日日啪夜夜爽| 日本午夜av视频| 最后的刺客免费高清国语| 少妇被粗大的猛进出69影院 | 大片免费播放器 马上看| 午夜视频国产福利| 国产精品国产三级国产av玫瑰| 一边摸一边做爽爽视频免费| 久久久久国产网址| 自拍欧美九色日韩亚洲蝌蚪91| 天天躁夜夜躁狠狠躁躁| 这个男人来自地球电影免费观看 | 永久免费av网站大全| 久久精品久久精品一区二区三区| 久久精品国产a三级三级三级| 天天操日日干夜夜撸| 高清不卡的av网站| 在线观看人妻少妇| 久久99蜜桃精品久久| 一区二区三区乱码不卡18| 亚洲精品一二三| 美女国产视频在线观看| 国产精品一区www在线观看| 国产探花极品一区二区| 最近最新中文字幕免费大全7| 秋霞在线观看毛片| 欧美xxⅹ黑人| 啦啦啦视频在线资源免费观看| 免费黄频网站在线观看国产| 蜜臀久久99精品久久宅男| 最后的刺客免费高清国语| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 我要看黄色一级片免费的| 精品视频人人做人人爽| 制服诱惑二区| 成年人午夜在线观看视频| 免费av不卡在线播放| 综合色丁香网| 日韩精品免费视频一区二区三区 | 国产成人免费无遮挡视频| 热99久久久久精品小说推荐| 久久这里只有精品19| 久久精品久久久久久久性| 亚洲精品日韩在线中文字幕| 精品少妇内射三级| 天天影视国产精品| 精品一区二区免费观看| 水蜜桃什么品种好| 色网站视频免费| 国国产精品蜜臀av免费| 精品人妻一区二区三区麻豆| 日本午夜av视频| 午夜激情av网站| 久久国产精品大桥未久av| 亚洲精品国产色婷婷电影| tube8黄色片| 亚洲伊人久久精品综合| 韩国高清视频一区二区三区| 亚洲精品,欧美精品| 99热这里只有是精品在线观看| 国产av码专区亚洲av| 男女国产视频网站| 久久久a久久爽久久v久久| 极品人妻少妇av视频| 啦啦啦啦在线视频资源| 亚洲欧美精品自产自拍| 精品久久久精品久久久| 久久久国产精品麻豆| 久久ye,这里只有精品| 亚洲欧美色中文字幕在线| 两性夫妻黄色片 | 考比视频在线观看| 国产精品欧美亚洲77777| videos熟女内射| 国产精品一国产av| 国产成人aa在线观看| 交换朋友夫妻互换小说| 欧美国产精品va在线观看不卡| 国产av精品麻豆| 亚洲 欧美一区二区三区| 国产伦理片在线播放av一区| 国产精品国产三级国产av玫瑰| 在线 av 中文字幕| 黄色一级大片看看| 美女xxoo啪啪120秒动态图| 人人妻人人爽人人添夜夜欢视频| 欧美精品av麻豆av| 国产精品久久久av美女十八| 精品一区二区免费观看| 日韩制服丝袜自拍偷拍| 又大又黄又爽视频免费| 高清毛片免费看| 制服诱惑二区| 久久99热6这里只有精品| 日韩av在线免费看完整版不卡| 亚洲av免费高清在线观看| 熟女电影av网| av国产久精品久网站免费入址| 伦理电影大哥的女人| 国精品久久久久久国模美| 国产毛片在线视频| 综合色丁香网| 国产欧美亚洲国产| 亚洲色图综合在线观看| 91精品三级在线观看| 欧美国产精品一级二级三级| 精品亚洲成a人片在线观看| 成人国语在线视频| 久久久久久久亚洲中文字幕| 国产精品人妻久久久影院| 少妇精品久久久久久久| 亚洲国产av新网站| av黄色大香蕉| 免费不卡的大黄色大毛片视频在线观看| 日韩制服丝袜自拍偷拍| 又大又黄又爽视频免费| 纯流量卡能插随身wifi吗| 日本色播在线视频| 精品一品国产午夜福利视频| 十分钟在线观看高清视频www| 久久精品国产亚洲av天美| 精品一区二区三区四区五区乱码 | 在线观看www视频免费| 人成视频在线观看免费观看| 欧美精品亚洲一区二区| 精品福利永久在线观看| 老司机亚洲免费影院| 高清黄色对白视频在线免费看| 99九九在线精品视频| 国产免费现黄频在线看| 性色avwww在线观看| 人体艺术视频欧美日本| 热re99久久国产66热| 国产男女超爽视频在线观看| 性色avwww在线观看| 全区人妻精品视频| 春色校园在线视频观看| 亚洲欧美一区二区三区黑人 | 国产片特级美女逼逼视频| 狠狠婷婷综合久久久久久88av| 99久久综合免费| 成人免费观看视频高清| 久久97久久精品| 精品少妇黑人巨大在线播放| 夫妻午夜视频| 少妇的逼水好多| 天堂8中文在线网| 国产探花极品一区二区| 黄色怎么调成土黄色| 免费看光身美女| 日日爽夜夜爽网站| www.av在线官网国产| 久久99热这里只频精品6学生| 久久久久国产网址| 一级毛片我不卡| 伊人亚洲综合成人网| 性色av一级| 久久这里只有精品19| 国产色婷婷99| 亚洲欧洲日产国产| 欧美人与善性xxx| 黄色怎么调成土黄色| 精品人妻熟女毛片av久久网站| 91在线精品国自产拍蜜月| 日韩,欧美,国产一区二区三区| av国产精品久久久久影院| 欧美日本中文国产一区发布| 黄色毛片三级朝国网站| 久久久久久久久久人人人人人人| 99久久精品国产国产毛片| 一级毛片黄色毛片免费观看视频| 日韩欧美精品免费久久| 国产成人精品无人区| 三级国产精品片| 最新中文字幕久久久久| 黑人巨大精品欧美一区二区蜜桃 | 久久精品国产亚洲av天美| 99久久中文字幕三级久久日本| 精品一区二区三卡| 两个人看的免费小视频| 久久鲁丝午夜福利片| 丝袜美足系列| 你懂的网址亚洲精品在线观看| a 毛片基地| 如日韩欧美国产精品一区二区三区| 国产深夜福利视频在线观看| 最近最新中文字幕免费大全7| 精品亚洲成国产av| av免费在线看不卡| 亚洲国产毛片av蜜桃av| 一级毛片电影观看| 看免费成人av毛片| 我要看黄色一级片免费的| 天天操日日干夜夜撸| 韩国高清视频一区二区三区| 日本午夜av视频| 中文字幕免费在线视频6| 日韩不卡一区二区三区视频在线| 国产日韩欧美亚洲二区| 国产亚洲欧美精品永久| 另类亚洲欧美激情| 久久97久久精品| videosex国产| 美国免费a级毛片| 天天躁夜夜躁狠狠久久av| 国产69精品久久久久777片| 国产欧美日韩综合在线一区二区| 亚洲伊人色综图| 国产一区亚洲一区在线观看| 国产一区有黄有色的免费视频| 国产精品三级大全| 国产毛片在线视频| 大片免费播放器 马上看| 国产免费一级a男人的天堂| 国产老妇伦熟女老妇高清| 最新的欧美精品一区二区| 九九爱精品视频在线观看| 18禁裸乳无遮挡动漫免费视频| 成年动漫av网址| 日本av手机在线免费观看| 欧美最新免费一区二区三区| www日本在线高清视频| 国产精品国产三级专区第一集| 三上悠亚av全集在线观看| av免费在线看不卡| 国产色婷婷99| 亚洲情色 制服丝袜| 久久精品人人爽人人爽视色| 大香蕉97超碰在线| 亚洲精品自拍成人| 中文字幕亚洲精品专区| av片东京热男人的天堂| 久久久精品免费免费高清| 亚洲精品久久午夜乱码| 精品国产露脸久久av麻豆| 午夜激情av网站| 国产熟女午夜一区二区三区| 日本91视频免费播放| 欧美另类一区| 在线天堂中文资源库| 中文字幕精品免费在线观看视频 | 熟女人妻精品中文字幕| 日本午夜av视频|