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

    離子液體在甲醇/硫酸介質(zhì)中對Q235鋼表面的緩蝕性能

    2020-06-06 01:27:46鄭天宇王璐劉金彥王佳
    化工學(xué)報 2020年5期
    關(guān)鍵詞:緩蝕劑硫酸原子

    鄭天宇,王璐,劉金彥,王佳

    (內(nèi)蒙古科技大學(xué)化學(xué)與化工學(xué)院,內(nèi)蒙古包頭014010)

    引 言

    甲醇作為石油的替代燃料引起了人們的廣泛關(guān)注,甲醇環(huán)境負(fù)荷低,不像其他化石燃料,會造成空氣污染和溫室效應(yīng)。在汽車行業(yè),燃料電池電動汽車?yán)眉状甲鳛闅湓碵1]。存放甲醇的罐體材料出于安全考慮,應(yīng)該對甲醇具有較高的緩蝕性。由于甲醇是偶極溶劑,可以作為很多物質(zhì)的良溶劑。靜態(tài)下,溶解的化學(xué)物質(zhì)很容易電離并產(chǎn)生電解液。因此,含有雜質(zhì)的甲醇很容易導(dǎo)致金屬容器中的腐蝕。如果是酸作為污染物存在于甲醇中,其腐蝕性會增加,因此針對甲醇的罐體材料的防腐研究很重要[2]。

    前些年出現(xiàn)了金屬在甲醇/水的酸性介質(zhì)中腐蝕行為的研究[3-5],但最近此方面的研究并不多見。在某些盛有甲醇/水酸性物質(zhì)的鐵容器中加入少量的緩蝕劑可以起到抗腐蝕作用,由于緩蝕劑在金屬表面的物理化學(xué)吸附作用,會在金屬表面形成保護(hù)層,所以對甲醇/水介質(zhì)環(huán)境不產(chǎn)生影響。

    目前離子液體作為環(huán)境友好的綠色溶劑出現(xiàn)在人們的視野中,由于其具有“可設(shè)計性”,低毒、無揮發(fā)性、熱穩(wěn)定性好、使用方便等一系列優(yōu)點[6],近年來陸續(xù)出現(xiàn)了其作為金屬緩蝕劑的研究[7]。文獻(xiàn)[8]報道通過1-甲基-3-(2-氨基乙烷)咪唑溴化銨離子液體改性檸檬酸基碳點,結(jié)果能夠有效抑制碳鋼在鹽酸以及氯化鈉溶液中的腐蝕;段騰龍等[9]研究表明1-甲基-3-(4-硫酸基丁基)咪唑硫酸鹽離子液體在鹽酸溶液中對銅有良好的緩蝕效果;張光華等[10]研究表明辛烷基二甲基芐基季銨鹽離子液體在鹽酸溶液中對碳鋼的緩蝕效率可達(dá)95%左右,但作為有機(jī)溶液金屬表面緩蝕劑的研究卻不多見。

    Q235 鋼作為低碳鋼,廣泛應(yīng)用于工業(yè)生產(chǎn)中,研究其在甲醇/硫酸水溶液中的腐蝕行為以及離子液體緩蝕劑對其腐蝕過程影響具有非常重要的意義。因此本文通過靜態(tài)失重法、電化學(xué)測試、掃描電子顯微鏡(SEM)研究甲醇中硫酸含量增加對Q235鋼的腐蝕影響,并探究1-丁基-3-甲基-咪唑氯鹽作為緩蝕劑時,Q235 鋼在甲醇/硫酸水溶液中的耐腐蝕性能,并通過量子化學(xué)計算和分子動力學(xué)模擬分析其緩蝕機(jī)理。

    1 實驗材料和方法

    1.1 實驗材料

    實驗所用金屬材料為Q235 鋼,其化學(xué)成分(質(zhì)量分?jǐn)?shù)):C 0.18%,Si 0.02%,Mn 0.45%,S 0.02%,P 0.01%,F(xiàn)e 余量[11]。所有試樣均需通過120#、240#、400#、600#、800#金相砂紙打磨,再用1200#金相砂紙在流動水下打磨至光滑鏡面。

    實驗所用試劑:1-丁基-3-甲基-咪唑氯鹽([Bmim]Cl,C8H15N2Cl,分子量174.67,99%純,中國科學(xué)院蘭州化學(xué)物理研究所),結(jié)構(gòu)式如圖1 所示;其余藥品:硫酸(98%),無水甲醇,無水乙醇,丙酮,乙二胺,環(huán)氧樹脂(E-44),均為AR,阿拉丁試劑公司(上海)。

    圖1 [Bmim]Cl的結(jié)構(gòu)式Fig.1 Chemical structure of[Bmim]Cl

    主要儀器:AutoLab302N 型電化學(xué)工作站,中國瑞士萬通公司;環(huán)境掃描電子顯微鏡Quanta400,荷蘭菲利浦公司。

    1.2 靜態(tài)失重法

    靜態(tài)失重法所使用的Q235 鋼試樣尺寸為20 mm×30 mm×5 mm,將打磨好的試樣用無水乙醇、丙酮依次沖洗,吹干,稱重,將稱重后的試樣分別浸入加與未加甲醇、緩蝕劑的溶液中,測試溫度為30℃,浸泡72 h 后,依次用無水乙醇、去離子水擦洗,干燥后稱重。腐蝕速率v和緩蝕效率η1分別按式(1)、式(2)計算[12]:

    式中,W為試樣實驗前后質(zhì)量的差值,g;A為試樣表面積,m2;t為實驗時間,h;v0、v′分別為空白試樣與加入緩蝕劑試樣的腐蝕速率,g·(m2·h)-1。

    1.3 電化學(xué)法

    電化學(xué)法分別為電化學(xué)阻抗譜測試和極化曲線法。在電化學(xué)工作站上使用三電極體系:Q235鋼(工作電極)、鉑片電極(輔助電極)、飽和甘汞電極(參比電極)。將尺寸為10 mm×10 mm×5 mm Q235鋼電極與銅導(dǎo)線進(jìn)行焊接,使用環(huán)氧樹脂與乙二胺,按照16∶1 的比例密封于PVC 管中,有效工作面積為1 cm2。實驗電極在測試溶液中浸泡1 d后進(jìn)行電化學(xué)測試。實驗結(jié)果通過Nova 和Origin 進(jìn)行擬合處理。

    動電位極化測試的參數(shù)設(shè)置:掃描速率0.5 mV·s-1,電位變化范圍-250~350 mV,緩蝕效率η2通過式(3)計算得出[13]:

    電化學(xué)阻抗譜參數(shù):掃描頻率范圍為0.1~105Hz,采用10 mV 開路電壓擾動信號,緩蝕效率η3通過式(4)計算得出[14]:

    1.4 掃描電子顯微鏡

    采用掃描電子顯微鏡(Quanta 400,20 kV 條件)分別觀察拋光后的Q235 鋼以及未添加與添加[Bmim]Cl 緩蝕劑的甲醇/硫酸混合水溶液中保持30℃浸泡3 d后的樣品表面的微觀結(jié)構(gòu)。

    1.5 量子化學(xué)計算

    通過Material Studio 7.0軟件中的DMol3模塊,采用GGA/PW91 算法的DND 基組對緩蝕劑的全電子進(jìn)行計算。收斂能量為2.0×105Ha;最大迭代次數(shù)為200;SCF收斂標(biāo)準(zhǔn)為1×105;最大SCF周期為200。對緩蝕劑分子進(jìn)行幾何結(jié)構(gòu)優(yōu)化,并計算分子的最高已占軌道能EHOMO、分子最低未占軌道能ELUMO、能量差ΔE、Fukui 指數(shù)來分析緩蝕劑分子結(jié)構(gòu)與緩蝕機(jī)理的關(guān)系[15-16]。

    1.6 分子動力學(xué)模擬

    分子動力學(xué)模擬(MD)通過Materials Studio 7.0中的Forcite 模塊完成。模擬過程如下:首先分別優(yōu)化并切割Fe(110)和Fe2O3(110)表面,然后建立一個Fe(110)為11×11×8、Fe2O3(110)為5×5×8 的超晶胞,然后使用A-morphous cell 模塊構(gòu)建包含1 個緩蝕劑分子,500個水分子,100個甲醇分子,2個硫酸分子,3個硫酸根離子,2個氫離子,4個水合氫根離子的溶液層。真空層厚度為20 ?(1 ?=0.1 nm),先對整個體系結(jié)構(gòu)進(jìn)行優(yōu)化,然后采用正則系統(tǒng)(NVT),模擬時間步長為1 fs,時間總長為300 ps,使用Andersen控溫器控制溫度在303 K,采用COMPASS 力場,固定所有Fe原子進(jìn)行動力學(xué)模擬。

    2 實驗結(jié)果與討論

    2.1 甲醇濃度的確定

    通過靜態(tài)失重法研究在純甲醇中添加不同體積的0.05 mol·L-1H2SO4溶液(保持腐蝕溶液總量為100 ml)后,對Q235 鋼腐蝕速率的影響,結(jié)果如表1所示。Vsulphuricacid表示100 ml 腐蝕溶液中0.05 mol·L-1H2SO4的添加量,Cmethanol表示添加硫酸后甲醇的濃度。

    表1 加入不同體積的0.05 mol·L-1 H2SO4后對Q235鋼在甲醇中的腐蝕速率Table 1 Corrosion rate of Q235 steel in methanol after adding 0.05 mol·L-1 H2SO4

    由表1可以看出,在純的甲醇中,幾乎對鋼片不腐蝕,但是隨著加入濃度為0.05 mol·L-1硫酸溶液的體積增大,腐蝕速率(v)加快,當(dāng)甲醇濃度達(dá)到10 mol·L-1時,v為1.30 g·(m2·h)-1,甲醇濃度為5 mol·L-1時,v為1.59 g·(m2·h)-1。由于質(zhì)子在甲醇中以H3O+的形式存在[17],當(dāng)硫酸的含量增加,導(dǎo)致質(zhì)子導(dǎo)電率和擴(kuò)散系數(shù)增加,其腐蝕速率增大。為了研究方便,在下面的研究中確定H2SO4的濃度為0.05 mol·L-1,甲醇的濃度為10 mol·L-1。

    2.2 [Bmim]Cl對Q235的緩蝕作用

    2.2.1 靜態(tài)失重法 表2 為添加了不同濃度的[Bmim]Cl 后對Q235 鋼在含有59.51 ml 0.05 mol·L-1H2SO4的甲醇溶液中(即10 mol·L-1甲醇硫酸水溶液)中的靜態(tài)失重數(shù)據(jù)。

    表2 不同濃度[Bmim]Cl對Q235鋼的靜態(tài)失重數(shù)據(jù)Table 2 Static mass-loss parameters of Q235 steel with different concentrations of[Bmim]Cl

    由表2 可以看出,[Bmim]Cl 能夠有效地減緩Q235鋼在甲醇/硫酸水溶液中的腐蝕,隨著[Bmim]Cl濃度的升高,緩蝕效率η1越來越好。當(dāng)濃度達(dá)到0.6 mol·L-1時,失重量(W)最小,η1達(dá)到最大值,為90.63%,說明此時Q235 鋼表面吸附量接近飽和,隨著濃度進(jìn)一步上升,η1開始下降,這可能是由于吸附達(dá)到過飽和而脫落所導(dǎo)致的。

    2.2.2 動電位極化曲線法 采用動電位極化曲線法測定不同濃度的[Bmim]Cl 對Q235 鋼在10 mol·L-1甲醇硫酸水溶液中的緩蝕性能,其極化曲線如圖2所示,相應(yīng)的擬合參數(shù)如表3所示。

    表3 中Ecorr為腐蝕電壓,Icorr為腐蝕電流密度,βc、βa分別為陰、陽極Tafel 斜率。由圖2 和表3 可知,空白溶液的Icorr為121.33 μA·cm-2,與未加緩蝕劑相比,從濃度為0.05 mol·L-1開始,隨著[Bmim]Cl 濃度的增加,Icorr逐漸降低,緩蝕效率η2越來越好;當(dāng)濃度達(dá)到0.6 mol·L-1時,Icorr達(dá)到最小值,為18.13 μA·cm-2,η2達(dá)到最高,為85.06%,說明此時Q235 鋼表面吸附的緩蝕劑接近飽和狀態(tài);隨著[Bmim]Cl 濃度的進(jìn)一步增大,Icorr開始升高,η2開始下降,說明此時吸附達(dá)到過飽和而脫落。βa和βc分別減少、增大,說明[Bmim]Cl對Q235 鋼在甲醇/硫酸水溶液中起到抑制作用,同時,Ecorr朝著正方向移動,且加入緩蝕劑與不加緩蝕劑的Ecorr差值小于80 mV,說明該緩蝕劑為陽極抑制為主的混合型抑制[18]。與上述靜態(tài)失重法測定結(jié)果相一致。

    圖2 不同濃度[Bmim]Cl對Q235鋼的極化曲線Fig.2 Polarization curves of Q235 steel with different concentrations of[Bmim]Cl

    表3 Q235鋼在不同[Bmim]Cl濃度下的極化曲線擬合參數(shù)Table 3 Polarization parameters for Q235 steel with different concentrations of[Bmim]Cl

    2.2.3 電化學(xué)阻抗法 采用電化學(xué)阻抗法測定在10 mol·L-1甲醇硫酸水溶液中不同濃度的[Bmim]Cl對Q235 鋼的Nyquist 譜圖如圖3 所示,相應(yīng)的擬合參數(shù)見表4。

    擬合的電化學(xué)等效電路圖如圖4 所示,圖中Rs為溶液電阻,Rp為電荷轉(zhuǎn)移電阻。χ2為擬合誤差,CPE 為電極和電解質(zhì)溶液之間的雙電層電容[19],CPE的阻抗值(ZCPE)由式(5)計算得出[20]:

    式中,Y0表示CPE 常數(shù),j表示虛根,ω表示角頻率,n表示彌散指數(shù)。當(dāng)n=1 時,常相位角元件表示純電容,此時電極表面達(dá)到最均勻的理想狀態(tài);當(dāng)n=0 時,常相位角元件表示純電阻。表4 中的n介于0.71~0.80 之間,說明[Bmim]Cl 緩蝕劑在Q235 鋼表面形成的吸附膜較均勻。

    圖3 不同濃度的[Bmim]Cl對Q235的Nyquist圖(a)、Bode圖(b)和相位角圖(c)Fig.3 Nyquist plots(a),Bode plots(b),and phase angle plots(c)for Q235 steel with different concentrations of[Bmim]Cl

    圖4 電化學(xué)等效電路圖Fig.4 Electrochemical equivalent circuit

    表4 不同濃度的[Bmim]Cl對Q235鋼的阻抗譜擬合參數(shù)Table 4 EIS parameters for Q235 steel with different concentrations of[Bmim]Cl

    從圖3(a)和表4 可得,空白組的Rp為157.41 Ω·cm-2,與不加緩蝕劑的空白組相比,從0.05 mol·L-1開始,隨著[Bmim]Cl 濃度的增大,容抗弧逐漸增大,相應(yīng)的Rp值也增大,緩蝕效率η3越來越好;當(dāng)濃度達(dá)到0.6 mol·L-1時,容抗弧達(dá)到最大,Rp值為1052.90 Ω·cm-2,η3達(dá)到85.05%,說明此時[Bmim]Cl 對Q235鋼的η3接近最佳值,隨著濃度進(jìn)一步的增大,η3開始下降。Rs隨著[Bmim]Cl 濃度的增加逐漸降低,說明溶液中[Bmim]Cl 的存在能夠有效地提高離子電導(dǎo)率。CPE 隨著濃度的增大先減小后增大,說明隨著[Bmim]Cl 含量增加,[Bmim]Cl 在金屬表面取代水分子的比例增大,吸附產(chǎn)生的雙電容層越致密,使Q235 鋼在溶液中的腐蝕得到了抑制。交流阻抗法與上述靜態(tài)失重法、極化曲線法的結(jié)果相吻合。在圖3(b)Bode 圖中,低頻處(0.1~102Hz)的阻抗絕對值(|Z|)也隨著[Bmim]Cl 濃度的增大而增加,濃度為0.6 mol·L-1時,|Z|達(dá)到最大值,說明此時緩蝕效果接近最佳狀態(tài);當(dāng)f>102Hz時,|Z|均趨于穩(wěn)定。圖3(c)相位角變化圖中,隨著[Bmim]Cl濃度的逐漸增加,相位角變得更負(fù)、更寬,說明[Bmim]Cl的加入有助于提高Q235 鋼的耐腐蝕性[21],Bode 圖和相位角變化與Nyquist圖結(jié)果一致。

    2.3 SEM分析

    圖5 Q235鋼的SEM形貌圖Fig.5 SEM micrographs of Q235 steel

    為了觀察不同狀態(tài)的Q235鋼表面腐蝕情況,將不同溶液浸泡后的碳鋼鍍金處理,再用環(huán)境掃描電鏡觀測表面。圖5 分別為不同狀態(tài)的Q235 鋼的SEM 形貌圖。圖5(a)為打磨光滑的Q235 鋼,表面基本平整;圖5(b)為Q235 鋼在10 mol·L-1甲醇硫酸水溶液中腐蝕后的形貌圖,可以看出表面腐蝕非常嚴(yán)重,表面極不平整;圖5(c)為添加了最佳濃度0.6 mol·L-1的[Bmim]Cl 作為緩蝕劑的形貌圖,可見低碳鋼表面腐蝕坑和裂痕明顯減少,表明[Bmim]Cl 具有良好的緩蝕效果;圖5(d)為添加了濃度為0.8 mol·L-1的[Bmim]Cl作為緩蝕劑的形貌圖,圖5(d)中位置1的表面形貌與圖5(c)表面類似,但圖5(d)中位置2、3 出現(xiàn)了與圖5(b)類似的腐蝕產(chǎn)物,說明在該濃度下緩蝕效率較低,這與之前的緩蝕性能測定相一致。另外,圖5(d)中位置4 出現(xiàn)了該緩蝕劑脫附的現(xiàn)象,說明可能是由于緩蝕劑濃度過高,導(dǎo)致其在金屬表面先逐漸達(dá)到吸附飽和狀態(tài),然后因吸附過飽和而脫落。

    2.4 量子化學(xué)計算

    為了進(jìn)一步研究[Bmim]Cl 的分子結(jié)構(gòu)對Q235鋼的緩蝕機(jī)理的影響,對[Bmim]Cl 單體分子以及[Bmim]+進(jìn)行幾何結(jié)構(gòu)優(yōu)化,圖6 分別為[Bmim]Cl 分子和[Bmim]+的幾何優(yōu)化結(jié)果,但考慮到離子液體在溶液中以離子狀態(tài)存在,故對[Bmim]+進(jìn)行HOMO、LUMO 電荷密度計算,結(jié)果如圖7 所示,相應(yīng)的量子化學(xué)參數(shù)如表5 所示。從圖7 可以看出,[Bmim]+的HOMO 主要分布在C C、C N 以及C—N 所形成的共軛π 鍵上,而LUMO 都主要分布在C C、C 原子、N 原子上,電子云密度都較大,且HOMO 和LUMO 電荷密度分布呈“不均勻”狀態(tài),說明緩蝕劑分子的咪唑環(huán)易形成強(qiáng)的活性吸附中心,而鐵原子的電子軌道分布為[Ar]4s23d6,最外層3d 軌道是未充滿的,4s軌道是充滿電子的,這種分布可使咪唑環(huán)與Fe原子上3d 空軌道形成配位鍵,同時又能夠接受Fe 原子上4s 軌道的電子形成反饋鍵,促使緩蝕劑分子緊密吸附于試樣表面,起到保護(hù)的作用。

    圖6 幾何優(yōu)化結(jié)構(gòu)圖Fig.6 Optimized structure

    圖7 [Bmim]+的前線分子軌道分布Fig.7 Frontier molecular orbital distribution of[Bmim]+

    表5 [Bmim]+的量子化學(xué)計算參數(shù)Table 5 Quantum chemical parameters for[Bmim]+

    根據(jù)量子化學(xué)的前線軌道理論[22],EHOMO與分子供電子能力有關(guān),其值越大,表示越可能提供電子給低能級或有空軌道的電子受體;ELUMO表示分子接收電子的能力,其值越小,分子越容易接受電子;二者軌道能之差?E=ELUMO-EHOMO是非常重要的穩(wěn)定性指標(biāo),其差值的絕對值越小,穩(wěn)定性越差,活性越強(qiáng),越容易發(fā)生交互作用[23]。其中Fe前線軌道能量來自文獻(xiàn)[24]。由表5 可知,由[Bmim]+的LUMO 躍遷到Fe 的HOMO所需的?E遠(yuǎn)遠(yuǎn)小于由[Bmim]+的HOMO 躍遷到Fe的LUMO 所需的?E,即[ELUMO([Bmim]+)-EHOMO(Fe)]<[ELUMO(Fe)-EHOMO([Bmim]+)],因此[Bmim]+相對于提供電子與Fe 原子作用而言,更易接受電子與Fe 原子作用,這也與[Bmim]+為正離子所一致[25]。

    為進(jìn)一步確定該緩蝕劑分子的反應(yīng)活性位點,表6 給出了[Bmim]Cl 分子的主要原子的Fukui 指數(shù)參數(shù)。f(r)-和f(r)+分別為親電攻擊指數(shù)和親核攻擊指數(shù),該值越大,表示分子中的原子給/得電子能力越強(qiáng)。從表6中可以看出,f(r)-較大的是N(1)、C(2)、N(5),說明在這些原子上形成了一個多中心的供電子區(qū)域[26],可向Fe原子的3d空軌道提供電子形成配位鍵;f(r)+較大的是C(3)、C(4)原子,說明這些原子形成了一個多中心的得電子區(qū)域,可接受Fe原子4s軌道上的電子形成反饋鍵。由此表明發(fā)生反應(yīng)的活性位點主要發(fā)生在咪唑環(huán)上,這與前線軌道分布結(jié)果相吻合。

    表6 由Hirshfeld電荷分布計算得到[Bmim]+的Fukui指數(shù)Table 6 Fukui index of[Bmim]+calculated by Hirshfeld charge distribution

    2.5 分子動力學(xué)模擬

    對[Bmim]Cl 進(jìn)行了分子動力學(xué)模擬,模擬是在Fe(110)和Fe2O3(110)表面上進(jìn)行的,因為它是最穩(wěn)定的[27],目的是為了更加接近實驗時的真實情況。

    圖8 [Bmim]Cl在Fe(110)和Fe2O3(110)上動力學(xué)模擬吸附的最低能量構(gòu)型Fig.8 Lowest energy configuration for[Bmim]Cl adsorption on Fe(110)and Fe2O3(110)

    表7 [Bmim]Cl在甲醇/硫酸水溶液中不同表面上的吸附能Table 7 Adsorption energy of[Bmim]Cl on different surfaces in methanol/sulfuric acid aqueous solutions

    [Bmim]Cl 分子在甲醇/硫酸水分子存在情況下,在Fe(110)和Fe2O3(110)表面吸附行為的模擬結(jié)果如圖8所示,由圖中可以看出,所有硫酸分子以及硫酸根離子都吸附在Fe(110)表面,說明硫酸會與Fe發(fā)生反應(yīng),造成Fe 表面的腐蝕現(xiàn)象;有部分甲醇分子吸附在金屬表面,這與前面討論2.1 節(jié)中的結(jié)果一致,甲醇分子的含量越多,金屬的腐蝕速率越??;緩蝕劑分子的陽離子[Bmim]+最終平行吸附在Fe(110)表面,這樣有助于增加表面覆蓋度,提高緩蝕效率;而陰離子Cl-則不吸附于Fe(110)表面,擴(kuò)散在溶液中,說明緩蝕劑分子中的Cl-未起到緩蝕的作用,原因是鐵表面在硫酸存在的情況下帶負(fù)電,不利于帶負(fù)電荷的Cl-吸附,因此脫附擴(kuò)散在溶液中。[Bmim]Cl 分子在甲醇/硫酸水分子存在情況下在Fe2O3(110)表面吸附行為與Fe(110)表面的吸附行為相同,全部硫酸分子和部分甲醇分子吸附在Fe2O3(110)表面上,而陽離子[Bmim]+的咪唑環(huán)平行吸附于金屬表面,有利于具有疏水性的可彎曲的烷基鏈在低碳鋼表面形成保護(hù)膜,防止溶液中的腐蝕介質(zhì)(水、硫酸等)擴(kuò)散到試樣表面,增加離子液體緩蝕劑膜的致密性,從而提高低碳鋼的耐腐蝕性[28]。陰離子Cl-同樣擴(kuò)散到溶液中。

    緩蝕劑與金屬表面之間的吸附能是衡量緩蝕劑分子是否具有良好吸附能力的指標(biāo),由式(6)計算得出[29],該值越負(fù),說明體系越穩(wěn)定,分子的緩蝕性能也就越好。

    根據(jù)分子動力學(xué)模擬所獲得的吸附能列于表7。由表中數(shù)據(jù)可知,緩蝕劑分子在Fe(110)表面的吸附能為-109.222 kcal·mol-1,在Fe2O3(110)表面的吸附能為-79.828 kcal·mol-1,數(shù)值均為負(fù)值,表明[Bmim]Cl吸附屬于自發(fā)過程[30]。

    3 結(jié) 論

    本文通過實驗和理論計算研究了離子液體[Bmim]Cl 在甲醇/硫酸介質(zhì)中對Q235 型鋼表面的緩蝕性能,得出如下結(jié)論。

    (1)碳鋼的腐蝕速率隨甲醇中硫酸含量的增加而增加,在10 mol·L-1甲醇硫酸水溶液中,[Bmim]Cl對Q235 鋼有很好的緩蝕效果,濃度為0.6 mol·L-1時,緩蝕效率達(dá)90%以上,靜態(tài)失重法和電化學(xué)法測試結(jié)果一致。[Bmim]Cl是以陽極抑制為主的混合型緩蝕劑。

    (2)量子化學(xué)計算表明[Bmim]Cl 具有很強(qiáng)的活性吸附中心,其咪唑環(huán)中含孤對電子的N 原子和共軛π 鍵為Q235 鋼試樣表面Fe 原子的3d 軌道提供電子形成配位鍵,同時能接受Fe 原子4s 軌道的電子形成反饋鍵,另外咪唑環(huán)的優(yōu)先吸附能使具有疏水性的烷基鏈在低碳鋼表面上形成保護(hù)膜,從而提高低碳鋼的耐腐蝕性。分子動力學(xué)模擬結(jié)果表明,緩蝕劑分子通過陽離子[Bmim]+平行吸附在金屬表面、Cl-擴(kuò)散到溶液中的方式,起到緩蝕的效果。

    符 號 說 明

    A——試樣表面積,m2

    C——腐蝕溶液中緩蝕劑的濃度,mol·L-1

    Cmethanol——添加硫酸后甲醇的濃度,mol·L-1

    CPE——雙電層電容,μF·cm-2

    Eadsorption——緩蝕劑的吸附能,kcal·mol-1

    Ecorr——腐蝕電壓,mV

    EHOMO——緩蝕劑分子的最高占有軌道能,eV

    Einhibitor/solution——緩蝕劑與腐蝕溶液的能量,kcal·mol-1

    ELUMO——緩蝕劑分子的最低未占軌道能,eV

    Esolution——腐蝕溶液的能量,kcal·mol-1

    Esurface/solution——鐵表面與腐蝕溶液的能量,kcal·mol-1

    Etotal——包含1 個緩蝕劑分子和金屬表面、腐蝕溶液的總能量,kcal·mol-1

    ?E——EHOMO和ELUMO的差值,eV

    f(r)-,f(r)+——分別為親電攻擊指數(shù)和親核攻擊指數(shù)

    ,I′corr——分別為不加與加緩蝕劑的腐蝕電流密度,μA·cm-2

    j——虛根

    n——彌散指數(shù)

    ,R′p——分別為不加與加緩蝕劑的電荷轉(zhuǎn)移電阻,Ω·cm-2

    Rs——溶液電阻,Ω·cm-2

    t——實驗時間,h

    Vsulphuricacid——100 ml 腐 蝕 溶 液 中0.05 mol·L-1H2SO4的 添加量,ml

    v,v′——分別為空白試樣與加入緩蝕劑試樣的腐蝕速率,g·(m2·h)-1

    W——試樣實驗前后質(zhì)量的差值,g

    ZCPE——CPE的阻抗值

    |Z|——阻抗絕對值,Ω·cm-2

    βa,βc——分別為陽極極化曲線斜率和陰極極化曲線斜率,mV·dec-1

    η——緩蝕效率,%

    χ2——擬合誤差

    ω——角頻率

    猜你喜歡
    緩蝕劑硫酸原子
    原子究竟有多???
    原子可以結(jié)合嗎?
    帶你認(rèn)識原子
    JRHS-2二氧化碳緩蝕劑的研究與應(yīng)用
    云南化工(2020年11期)2021-01-14 00:50:46
    硫酸渣直接還原熔分制備珠鐵
    2019年硫酸行業(yè)運行情況
    2018磷復(fù)肥硫酸大數(shù)據(jù)發(fā)布
    新型緩蝕劑的性能研究
    硫酸很重要特點要知道
    一種適用于鐵基脫硫液的緩蝕劑
    山西化工(2016年6期)2016-04-09 07:17:41
    69人妻影院| 两个人视频免费观看高清| 欧美日韩国产亚洲二区| 亚洲精品456在线播放app | 国产成人影院久久av| 色吧在线观看| av福利片在线观看| 乱码一卡2卡4卡精品| 亚洲欧美精品综合久久99| 成人av在线播放网站| 九九热线精品视视频播放| 日本 av在线| 国产精品98久久久久久宅男小说| 在线免费十八禁| 97热精品久久久久久| 国产高清有码在线观看视频| 波多野结衣高清作品| 美女大奶头视频| 免费观看在线日韩| 黄色一级大片看看| 欧美日韩国产亚洲二区| 亚洲第一电影网av| 成熟少妇高潮喷水视频| 国产真实乱freesex| 中国美女看黄片| 午夜影院日韩av| 国产毛片a区久久久久| 99热这里只有精品一区| 啦啦啦啦在线视频资源| 18禁在线播放成人免费| 久久人人精品亚洲av| 国产亚洲精品综合一区在线观看| 无人区码免费观看不卡| 禁无遮挡网站| 欧美在线一区亚洲| 一级毛片久久久久久久久女| 日韩欧美三级三区| 1000部很黄的大片| 欧美精品啪啪一区二区三区| 免费观看精品视频网站| 亚洲第一电影网av| 精品午夜福利视频在线观看一区| 成年人黄色毛片网站| 国产黄色小视频在线观看| 欧美性猛交黑人性爽| 在线免费观看的www视频| 亚洲精品成人久久久久久| 日本与韩国留学比较| 精品人妻1区二区| 天堂影院成人在线观看| 亚洲国产色片| 久久精品影院6| 亚洲色图av天堂| 国产69精品久久久久777片| 性色avwww在线观看| 国产av在哪里看| 欧美潮喷喷水| 国产在视频线在精品| 成人性生交大片免费视频hd| 男女啪啪激烈高潮av片| 99久久成人亚洲精品观看| 午夜福利成人在线免费观看| 亚洲一区高清亚洲精品| 国产精品美女特级片免费视频播放器| 国产精品美女特级片免费视频播放器| 久9热在线精品视频| 日韩亚洲欧美综合| 久久九九热精品免费| 亚洲精品久久国产高清桃花| 亚洲av日韩精品久久久久久密| 无遮挡黄片免费观看| 不卡一级毛片| 动漫黄色视频在线观看| 欧美最新免费一区二区三区| 一级黄色大片毛片| 国产白丝娇喘喷水9色精品| 国产免费一级a男人的天堂| 色在线成人网| 五月伊人婷婷丁香| 91av网一区二区| 亚洲av中文av极速乱 | 婷婷色综合大香蕉| 国产三级中文精品| 精品福利观看| 亚洲人成网站高清观看| 午夜爱爱视频在线播放| 少妇裸体淫交视频免费看高清| 亚洲av免费在线观看| 亚洲乱码一区二区免费版| 欧美区成人在线视频| 国产毛片a区久久久久| 国产激情偷乱视频一区二区| 久久久国产成人精品二区| 久久久色成人| 嫩草影院入口| 国产av不卡久久| 午夜福利18| 精品一区二区三区视频在线观看免费| 国产高清不卡午夜福利| 欧美精品啪啪一区二区三区| 老司机深夜福利视频在线观看| 国产探花在线观看一区二区| 中文资源天堂在线| 婷婷精品国产亚洲av| 国内精品美女久久久久久| 国产高清视频在线播放一区| 两个人的视频大全免费| 深夜a级毛片| 精品久久久久久,| 此物有八面人人有两片| 亚洲乱码一区二区免费版| 少妇高潮的动态图| 国产在线男女| netflix在线观看网站| 九九在线视频观看精品| 波多野结衣高清无吗| 此物有八面人人有两片| 婷婷色综合大香蕉| 成年女人永久免费观看视频| 精品久久国产蜜桃| 高清日韩中文字幕在线| 国产高清视频在线观看网站| 淫秽高清视频在线观看| 一个人免费在线观看电影| 在线观看美女被高潮喷水网站| 免费大片18禁| 美女xxoo啪啪120秒动态图| 他把我摸到了高潮在线观看| 国产 一区精品| 我的女老师完整版在线观看| 久久久久久久午夜电影| 岛国在线免费视频观看| 亚洲欧美日韩高清专用| 国产v大片淫在线免费观看| 亚洲人成网站高清观看| 精品久久久噜噜| 久久久久久久亚洲中文字幕| 男女那种视频在线观看| 最近最新中文字幕大全电影3| 男女那种视频在线观看| 亚洲在线观看片| 欧美一级a爱片免费观看看| 小蜜桃在线观看免费完整版高清| 麻豆一二三区av精品| 中文字幕av在线有码专区| 久久久久久久久中文| 一个人看视频在线观看www免费| 成人鲁丝片一二三区免费| 一本一本综合久久| 国产亚洲91精品色在线| 国产欧美日韩精品亚洲av| 一进一出抽搐gif免费好疼| 老司机深夜福利视频在线观看| 尾随美女入室| 亚洲av免费高清在线观看| 一边摸一边抽搐一进一小说| 国产亚洲精品久久久com| 99久国产av精品| 亚洲成a人片在线一区二区| 69av精品久久久久久| 狂野欧美激情性xxxx在线观看| 欧美性猛交╳xxx乱大交人| 日本免费a在线| 男女之事视频高清在线观看| 男人和女人高潮做爰伦理| 高清毛片免费观看视频网站| 毛片女人毛片| 午夜激情欧美在线| 国产精品伦人一区二区| 国产色爽女视频免费观看| 大又大粗又爽又黄少妇毛片口| 国产在线男女| 国产白丝娇喘喷水9色精品| 不卡一级毛片| 欧美一区二区精品小视频在线| 日韩中字成人| 日本精品一区二区三区蜜桃| 午夜免费激情av| 亚洲第一电影网av| 精品福利观看| 3wmmmm亚洲av在线观看| 久久精品久久久久久噜噜老黄 | 蜜桃亚洲精品一区二区三区| 好男人在线观看高清免费视频| 69人妻影院| 观看美女的网站| 国产精品日韩av在线免费观看| 精品乱码久久久久久99久播| 在线免费观看的www视频| 免费不卡的大黄色大毛片视频在线观看 | 成人性生交大片免费视频hd| 嫁个100分男人电影在线观看| 成人无遮挡网站| 亚洲黑人精品在线| 精品久久久久久久末码| 狠狠狠狠99中文字幕| 国模一区二区三区四区视频| 欧美色视频一区免费| 热99re8久久精品国产| 亚洲欧美精品综合久久99| 国产在线精品亚洲第一网站| 中文字幕av在线有码专区| 国产熟女欧美一区二区| 一区二区三区激情视频| 亚洲av五月六月丁香网| 久久精品国产亚洲av涩爱 | 3wmmmm亚洲av在线观看| 久久久久久伊人网av| 欧美bdsm另类| 久久精品国产清高在天天线| 欧美又色又爽又黄视频| www日本黄色视频网| 国产黄a三级三级三级人| 国产精品女同一区二区软件 | 赤兔流量卡办理| 国产毛片a区久久久久| 国产高清视频在线播放一区| 美女大奶头视频| 亚洲一级一片aⅴ在线观看| 欧美性猛交黑人性爽| 国产精品久久视频播放| 一级黄片播放器| 91久久精品国产一区二区成人| 夜夜爽天天搞| 乱人视频在线观看| 最后的刺客免费高清国语| 狠狠狠狠99中文字幕| 国模一区二区三区四区视频| 国产精品伦人一区二区| 国产人妻一区二区三区在| 成人欧美大片| 久久久久久久久久久丰满 | 三级国产精品欧美在线观看| 搡女人真爽免费视频火全软件 | 国产精品综合久久久久久久免费| 久久精品国产自在天天线| 国产精品无大码| 男女做爰动态图高潮gif福利片| 亚洲一级一片aⅴ在线观看| 亚洲国产欧洲综合997久久,| 一级毛片久久久久久久久女| 亚洲成人久久爱视频| 日本在线视频免费播放| 少妇高潮的动态图| 91av网一区二区| 一级毛片久久久久久久久女| 真人一进一出gif抽搐免费| 午夜精品一区二区三区免费看| 三级国产精品欧美在线观看| 最新在线观看一区二区三区| 麻豆一二三区av精品| 成人av一区二区三区在线看| 久久久久久伊人网av| 色吧在线观看| av在线天堂中文字幕| 久久久色成人| 亚洲av不卡在线观看| 国产精品一区二区免费欧美| 欧美激情久久久久久爽电影| 人人妻,人人澡人人爽秒播| 亚洲国产欧洲综合997久久,| 俺也久久电影网| 如何舔出高潮| 国产免费男女视频| 国产女主播在线喷水免费视频网站 | 亚洲精品一区av在线观看| 91麻豆精品激情在线观看国产| 黄色女人牲交| 麻豆国产97在线/欧美| 欧美日韩中文字幕国产精品一区二区三区| 听说在线观看完整版免费高清| 我要搜黄色片| 在线观看午夜福利视频| 久久久国产成人免费| 成年免费大片在线观看| 成人亚洲精品av一区二区| 亚洲av成人精品一区久久| 中文资源天堂在线| 日韩中文字幕欧美一区二区| 亚洲真实伦在线观看| 国产精品,欧美在线| 久久精品综合一区二区三区| 一夜夜www| 男人的好看免费观看在线视频| 中文字幕精品亚洲无线码一区| 美女高潮的动态| av在线老鸭窝| 国产伦精品一区二区三区四那| 看黄色毛片网站| 精品一区二区三区视频在线观看免费| 欧美一区二区精品小视频在线| 日韩欧美一区二区三区在线观看| 深爱激情五月婷婷| 亚洲成人久久爱视频| 国产高清三级在线| 国产精品一区二区三区四区免费观看 | 免费电影在线观看免费观看| 国产私拍福利视频在线观看| 有码 亚洲区| 亚洲黑人精品在线| aaaaa片日本免费| 亚洲av成人精品一区久久| 国产毛片a区久久久久| 一进一出抽搐gif免费好疼| 亚洲精品成人久久久久久| 淫秽高清视频在线观看| 1024手机看黄色片| 欧洲精品卡2卡3卡4卡5卡区| 搡老熟女国产l中国老女人| 一区二区三区高清视频在线| 久久精品国产自在天天线| 99九九线精品视频在线观看视频| 亚洲熟妇熟女久久| 天美传媒精品一区二区| 我要搜黄色片| 国产精品一区www在线观看 | 成人国产综合亚洲| 国产精品综合久久久久久久免费| 神马国产精品三级电影在线观看| eeuss影院久久| 国产精品一区二区三区四区久久| АⅤ资源中文在线天堂| 最新在线观看一区二区三区| 中出人妻视频一区二区| 人妻制服诱惑在线中文字幕| av.在线天堂| 国产爱豆传媒在线观看| 亚洲成人精品中文字幕电影| or卡值多少钱| 日本一二三区视频观看| 啦啦啦啦在线视频资源| 搡老熟女国产l中国老女人| 久久精品人妻少妇| 哪里可以看免费的av片| 春色校园在线视频观看| 久久精品国产鲁丝片午夜精品 | 久久久久久久久久黄片| videossex国产| 欧洲精品卡2卡3卡4卡5卡区| 亚洲va日本ⅴa欧美va伊人久久| 日本免费一区二区三区高清不卡| 国产私拍福利视频在线观看| 麻豆成人午夜福利视频| 国产av一区在线观看免费| 91在线观看av| 又粗又爽又猛毛片免费看| 91久久精品国产一区二区三区| 精品久久久久久久末码| 久久久国产成人精品二区| 国产亚洲精品av在线| 12—13女人毛片做爰片一| 久久99热这里只有精品18| 最近中文字幕高清免费大全6 | 国产久久久一区二区三区| 网址你懂的国产日韩在线| 午夜福利18| 一夜夜www| 国内精品宾馆在线| 精品久久久噜噜| 欧美日本亚洲视频在线播放| 97热精品久久久久久| 午夜视频国产福利| 日日啪夜夜撸| 色播亚洲综合网| 国产精品野战在线观看| 日韩中字成人| 免费观看的影片在线观看| 成人精品一区二区免费| 亚洲精品一区av在线观看| 国内毛片毛片毛片毛片毛片| 午夜激情福利司机影院| 国产午夜精品论理片| 亚洲第一区二区三区不卡| 一a级毛片在线观看| 日本免费a在线| 97超级碰碰碰精品色视频在线观看| 国产亚洲精品av在线| 午夜免费成人在线视频| 女生性感内裤真人,穿戴方法视频| 51国产日韩欧美| 看黄色毛片网站| 人人妻人人澡欧美一区二区| 18禁在线播放成人免费| 美女高潮的动态| 伊人久久精品亚洲午夜| 日日撸夜夜添| 国产一区二区激情短视频| 日本 欧美在线| 国产精品伦人一区二区| 伦精品一区二区三区| 欧美日韩精品成人综合77777| 少妇裸体淫交视频免费看高清| 免费看日本二区| 又黄又爽又免费观看的视频| 91av网一区二区| 亚洲av一区综合| 最新在线观看一区二区三区| 熟妇人妻久久中文字幕3abv| 国产 一区精品| 一区二区三区四区激情视频 | 成人精品一区二区免费| 神马国产精品三级电影在线观看| 国产一区二区三区在线臀色熟女| aaaaa片日本免费| 欧美日韩国产亚洲二区| 最新中文字幕久久久久| 色精品久久人妻99蜜桃| 欧美丝袜亚洲另类 | 女生性感内裤真人,穿戴方法视频| 免费看av在线观看网站| 变态另类成人亚洲欧美熟女| 亚洲精品国产成人久久av| 欧美国产日韩亚洲一区| 国产精品乱码一区二三区的特点| 在线免费观看的www视频| 成年免费大片在线观看| 国产一级毛片七仙女欲春2| 在线观看美女被高潮喷水网站| 久久久久久久亚洲中文字幕| 国产一区二区三区av在线 | 久久久午夜欧美精品| 男插女下体视频免费在线播放| 国产探花极品一区二区| 亚洲av美国av| 成人美女网站在线观看视频| 国产精品三级大全| 亚洲国产日韩欧美精品在线观看| 校园人妻丝袜中文字幕| 老熟妇乱子伦视频在线观看| 日韩欧美精品免费久久| 色哟哟哟哟哟哟| 成人性生交大片免费视频hd| 99国产极品粉嫩在线观看| 97超级碰碰碰精品色视频在线观看| 2021天堂中文幕一二区在线观| 亚洲欧美激情综合另类| 91av网一区二区| 成人亚洲精品av一区二区| 亚洲美女视频黄频| 国产大屁股一区二区在线视频| 香蕉av资源在线| 国产精品一区二区免费欧美| 亚洲熟妇中文字幕五十中出| 中文字幕av在线有码专区| 欧美性猛交╳xxx乱大交人| 男女那种视频在线观看| 一进一出抽搐gif免费好疼| 欧美xxxx性猛交bbbb| 日日啪夜夜撸| 人妻久久中文字幕网| 久久99热这里只有精品18| 久久人妻av系列| 久久人人爽人人爽人人片va| 国产亚洲精品av在线| 国产黄色小视频在线观看| 老熟妇乱子伦视频在线观看| 美女cb高潮喷水在线观看| 三级国产精品欧美在线观看| 亚洲人与动物交配视频| 亚洲美女黄片视频| 亚洲精品乱码久久久v下载方式| 他把我摸到了高潮在线观看| 免费在线观看日本一区| 欧美性猛交╳xxx乱大交人| 少妇熟女aⅴ在线视频| 中亚洲国语对白在线视频| 夜夜看夜夜爽夜夜摸| 别揉我奶头~嗯~啊~动态视频| 国产精品久久久久久亚洲av鲁大| 如何舔出高潮| 亚洲一级一片aⅴ在线观看| 啦啦啦韩国在线观看视频| 综合色av麻豆| 少妇熟女aⅴ在线视频| 欧美在线一区亚洲| 夜夜看夜夜爽夜夜摸| 成人性生交大片免费视频hd| 亚洲自拍偷在线| 内地一区二区视频在线| 久久久久久伊人网av| 免费在线观看日本一区| 精品久久国产蜜桃| 亚洲三级黄色毛片| 久久精品国产自在天天线| 99久国产av精品| 亚洲自拍偷在线| 欧美丝袜亚洲另类 | 成人av在线播放网站| 一级黄色大片毛片| 99国产精品一区二区蜜桃av| 欧美日本亚洲视频在线播放| 亚洲人与动物交配视频| 亚洲专区中文字幕在线| 禁无遮挡网站| 在线观看66精品国产| 两个人视频免费观看高清| 亚洲性久久影院| 国产精品自产拍在线观看55亚洲| av专区在线播放| 国产黄a三级三级三级人| 国产色婷婷99| 亚洲18禁久久av| 中文字幕人妻熟人妻熟丝袜美| 一区福利在线观看| 午夜福利成人在线免费观看| 18禁黄网站禁片免费观看直播| 日本黄色片子视频| av天堂在线播放| 无人区码免费观看不卡| 亚洲专区中文字幕在线| 国产精品一区www在线观看 | 俺也久久电影网| 在线观看免费视频日本深夜| 久久久久国内视频| 九九在线视频观看精品| 欧美xxxx黑人xx丫x性爽| 白带黄色成豆腐渣| 观看美女的网站| 日韩在线高清观看一区二区三区 | 欧美色视频一区免费| 久久精品91蜜桃| 国模一区二区三区四区视频| 精品不卡国产一区二区三区| 校园人妻丝袜中文字幕| 搡老熟女国产l中国老女人| 欧美日韩中文字幕国产精品一区二区三区| 成人综合一区亚洲| 十八禁网站免费在线| 黄色欧美视频在线观看| 欧美黑人欧美精品刺激| 天堂影院成人在线观看| 精品一区二区免费观看| 久久精品国产鲁丝片午夜精品 | 中出人妻视频一区二区| 日本色播在线视频| 免费看av在线观看网站| 亚洲精品日韩av片在线观看| 国产v大片淫在线免费观看| 美女免费视频网站| 成人亚洲精品av一区二区| eeuss影院久久| 免费黄网站久久成人精品| 精品国产三级普通话版| 久久久久性生活片| 亚洲真实伦在线观看| 最新在线观看一区二区三区| 亚洲美女黄片视频| 极品教师在线视频| 久久精品国产自在天天线| 欧洲精品卡2卡3卡4卡5卡区| 久久精品国产亚洲av香蕉五月| 国产精品一区二区性色av| 村上凉子中文字幕在线| 午夜激情欧美在线| 国产毛片a区久久久久| 99久久无色码亚洲精品果冻| 日本熟妇午夜| 99久久无色码亚洲精品果冻| h日本视频在线播放| 丝袜美腿在线中文| 国产精品国产高清国产av| 亚洲黑人精品在线| 婷婷精品国产亚洲av| 成人无遮挡网站| 成年女人看的毛片在线观看| 真人一进一出gif抽搐免费| 免费在线观看影片大全网站| 简卡轻食公司| 久久欧美精品欧美久久欧美| 色视频www国产| 久久欧美精品欧美久久欧美| 免费观看的影片在线观看| 国产白丝娇喘喷水9色精品| 国产精品日韩av在线免费观看| 观看免费一级毛片| 床上黄色一级片| 中文字幕av在线有码专区| or卡值多少钱| 日韩,欧美,国产一区二区三区 | 一边摸一边抽搐一进一小说| 不卡一级毛片| 中文字幕精品亚洲无线码一区| 日韩人妻高清精品专区| 亚洲最大成人手机在线| 成人国产综合亚洲| 欧美性感艳星| 精品国内亚洲2022精品成人| 男人舔女人下体高潮全视频| 热99re8久久精品国产| 久久精品91蜜桃| 少妇猛男粗大的猛烈进出视频 | 最近最新免费中文字幕在线| 国产精品不卡视频一区二区| 俄罗斯特黄特色一大片| 免费看光身美女| 中亚洲国语对白在线视频| 夜夜看夜夜爽夜夜摸| 成人性生交大片免费视频hd| 免费看a级黄色片| 国产精品久久电影中文字幕| 国产精品国产三级国产av玫瑰| 免费av观看视频| 国产亚洲av嫩草精品影院| 久久人人精品亚洲av| 国产精品一区二区免费欧美| 成人午夜高清在线视频| 日本免费一区二区三区高清不卡| 久久精品国产鲁丝片午夜精品 | 麻豆成人av在线观看| 欧美丝袜亚洲另类 | 国产午夜精品论理片| 好男人在线观看高清免费视频| 久久精品夜夜夜夜夜久久蜜豆| 国产亚洲精品综合一区在线观看| 久久国内精品自在自线图片|