999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

2,3-二羥基萘-6-磺酸鈉接枝吸附劑去除水中微量硼的機理

2015-12-13 00:50:34余肖峰胡晶晶溫曉華邵超英
上海國土資源 2015年1期
關鍵詞:化學

余肖峰,顧 華,胡晶晶,溫曉華,張 琢,邵超英

(1. 上海市地質調查研究院,上海 200072;2. 國土資源部上海資源環境監督檢測中心,上海 200072;3. 東華大學化學化工與生物工程學院,上海 201620)

2,3-二羥基萘-6-磺酸鈉接枝吸附劑去除水中微量硼的機理

余肖峰1,2,顧 華1,2,胡晶晶3,溫曉華1,2,張 琢1,2,邵超英3

(1. 上海市地質調查研究院,上海 200072;2. 國土資源部上海資源環境監督檢測中心,上海 200072;3. 東華大學化學化工與生物工程學院,上海 201620)

以11B NMR法研究了硼酸與2,3-二羥基萘-6-磺酸鈉的配位反應。結果表明,2,3-二羥基萘-6-磺酸鈉與水溶液中的硼酸根離子反應生成五元環的1:1和1:2(化學計量比)的配合物,配合物的形成常數分別為logβ1=5.22、logβ2=6.10。利用離子交換法制備了2,3-二羥基萘-6-磺酸鈉接枝硼選擇性吸附劑,吸附劑上2,3-二羥基萘-6-磺酸鈉擔載量為1.20mmol/g。以bath法研究了2,3-二羥基萘-6-磺酸鈉接枝吸附劑對水中硼的去除作用,考察了溶液pH值、初始硼濃度、溫度、吸附時間等因素對硼去除率的影響。在該吸附劑吸附硼的熱力學和動力學研究中發現,硼的吸附符合準二級速率方程,是化學吸附和物理吸附共同作用的結果,B(OH)4-以1:1配合物的形式和靜電吸引的方式結合在吸附劑上;吸附速率主要為化學吸附控制。

水環境;硼污染;2,3-二羥基萘-6-磺酸鈉;吸附劑;硼酸/硼酸鹽;去除機理

微量元素硼是生物體必需的營養元素,但過量會導致生物體產生毒理反應。隨著金屬冶煉、玻璃、日化和制藥等工業的發展,水環境中硼的污染日趨嚴重。因此,硼去除技術的研究越來越受到關注。

天然水中,硼主要以硼酸/硼酸鹽的形式存在[1]:

使用硼特效吸附劑是去除水中硼最有效的方法之一,去除機理可歸結為吸附劑對硼的吸附和其結構中多羥基官能團與硼酸/硼酸鹽絡合的共同作用。該絡合反應同時受熱力學和動力學控制[2]。由于該法有硼去除效率高、被吸附硼解吸回收便利和吸附劑可再生循環利用等優點,已被成功應用到從海水中提取硼以及從廢水中去除硼的工作中[3~8]。

本課題選擇了典型的多羥基化合物:2,3-二羥基萘-6-磺酸鈉作為吸附劑接枝化合物,以研究和開發吸附容量高、選擇性好的硼特效選擇性吸附劑。通過11B NMR分析探討了吸附機理,用bath法研究了溶液pH值及硼酸初始濃度等對吸附效果的影響,并對吸附過程的熱力學和動力學進行了分析。

1 實驗部分

1.1 儀器和試劑

儀器:JEOL JNM-GSX 500型核磁共振光譜儀(日本電子株式會社);UV 1102型紫外可見分光光度計(上海天美科學儀器有限公司);Orion 4 Star 型pH測定儀(美國熱電公司);EL 204型電子天平(梅特勒—托利多儀器(上海)有限公司);SHZ-B型水浴恒溫振蕩器(上海躍進醫療器械廠)。

試劑:DOWEX 1-X8型強堿性陰離子交換樹脂(和光純藥工業株式會社);2,3-二羥基萘-6-磺酸鈉(和光純藥工業株式會社);硼酸(上海凌峰化學試劑有限公司)。所用試劑均為市售分析純試劑,實驗用水為二次去離子水。

1.211B NMR測定

制備了含已知量硼酸(0.0185mol/L)和2,3-二羥基萘-6-磺酸鈉(0.0185mol/L)的混合溶液(離子強度I=0.2mol/ L, NaCl),溶液的pH值用CH3CO2H-CH3CO2Na、NaH2PO4-Na2HPO4、HEPES、NH3-NH4Cl緩沖溶液或少量HCl和NaOH調節。達平衡后,用于11B NMR的測定。11B NMR測試在共振頻率為160.0 MHz、10mm多核探頭、24±1℃的條件下進行。標準NMR參數如下:翻轉角度約為90?(36μs),脈沖重復時間為1s,譜寬為31kHz。化學位移以0.1mol/L的硼酸溶液作為外部參比。重疊的NMR信號用Lorentzian曲線擬合法解析。

1.3 硼選擇性吸附劑的制備

(1)接枝吸附劑制備條件試驗

實驗選擇DOWEX 1-X8樹脂(交聯度8%、總交換容量為1.2meq/mL)作為接枝樹脂。樹脂經0.1mol/L NaCl再生處

理后,濾紙壓干,備用。

稱取0.2g上述樹脂于聚乙烯瓶中,分別加入20mL pH值4~10、濃度為0.001mol/L的2,3-二羥基萘-6-磺酸鈉溶液,分別在離子強度(NaCl)為0.02~0.50mol/L、25℃條件下振蕩1~24h進行接枝樹脂的制備試驗。經濾膜過濾后,根據上清液在237nm處的紫外吸收測定其中2,3-二羥基萘-6-磺酸鈉的濃度。樹脂上2,3-二羥基萘-6-磺酸鈉的分配比kd依下式計算:

式中:C0和Ce分別為水相中2,3-二羥基萘-6-磺酸鈉初始和平衡濃度,mol/L;m為樹脂質量,g;V為溶液體積,mL。

(2)2,3-二羥基萘-6-磺酸鈉接枝吸附劑的制備

稱取6.0g2,3-二羥基萘-6-磺酸鈉溶解于250mL去離子水中,用0.5mol/L NaOH溶液調節pH值為7.4(I=0.04mol/L),加入樹脂12g,恒溫(25℃)振蕩24h。達平衡后,過濾,經濾紙壓干,稱量得接枝樹脂16.9g,即吸附劑中2,3-二羥基萘-6-磺酸鈉的擔載量為1.2mmol/g。

1.4 接枝吸附劑對硼的吸附條件試驗

稱取0.5g上述接枝吸附劑于聚乙烯瓶中,分別加入20mL濃度為0.005mol/L、0.01mol/L、0.02mol/L、0.03mol/ L、0.05mol/L、0.07mol/L的硼酸溶液(I=0.1mol/L,NaCl),調節溶液pH值至4~13、吸附時間0.5h、3h、6h、12h、24h、48h、反應溫度25℃、30℃、35℃、40℃、45℃,進行硼的吸附去除試驗。達平衡后,過濾上清液,用甲亞胺-H分光光度法測定溶液中硼的濃度。按式(3)計算接枝吸附劑上硼的吸附量:

式中:qe為硼在接枝吸附劑上的平衡吸附量,mmol/ g;C0和Ce分別為水相中硼酸的初始和平衡濃度,mol/L;m為樹脂質量,g;V為溶液體積,mL。

2 結果與討論

2.111B NMR研究

11B NMR化學位移能直觀反映溶液中硼的存在形態,因此,該技術在硼化學研究中獲得了很好的應用。由硼酸-2,3-二羥基萘-6-磺酸鈉混合溶液的11B NMR譜圖(圖1)可見,除B(OH)3/B(OH)4-的NMR信號(0ppm,-10.9ppm)外,在-5.8ppm和-11.3ppm處出現了兩個新峰。據報道[9~11],硼酸與多羥基化合物反應生成的五元環1:1和1:2配合物的化學位移分別約為-11.0ppm和-5.5ppm,因此,-11.3ppm處的NMR信號被歸結為B(OH)4-與2,3-二羥基萘-6-磺酸鈉反應形成的1:1配合物;而-5.8ppm處的NMR信號為B(OH)4-與2,3-二羥基萘-6-磺酸鈉的1:2配合物。

圖1 不同pH值條件下硼酸-2,3-二羥基萘-6-磺酸鈉混合溶液的11B NMR 圖譜a: B(OH)3/B(OH)4-;b: 1:2配合物;c: 1:1配合物Fig.111B NMR spectra for the solution of boric acid/borate and 2,3-dihydroxynaphthalene-6-sodium sulfonate at different pH value

硼酸-2,3-二羥基萘-6-磺酸鈉混合溶液中存在式(1)及(4)~(7)平衡:

式中,B(OH)2L2-表示B(OH)4-與2,3-二羥基萘-6-磺酸鈉反應生成的1:1配合物;BL23-為其1:2配合物。根據各形態硼的NMR相對信號強度及式(8)、(9),可計算出該反應生成的1:1和1:2配合物的形成常數分別為:logβ1=5.22、logβ2=6.10。

根據硼酸的Ka和β1、β2,可進一步計算一定濃度硼酸-2,3-二羥基萘-6-磺酸鈉混合溶液中各形態硼的濃度隨溶液pH值變化曲線(圖2)。可知,1:1配合物的濃度隨著pH值的上升逐漸升高,在pH為9左右達到最大值;而1:2配合物在整個pH范圍內濃度均很低。這與11B NMR測定結果基本一致。由此可得出結論,1:1配合物主要存在于堿性溶液中,而1:2配合物則難以形成。該結果為后續2,3-二羥基萘-6-磺酸鈉接枝硼選擇性吸附劑去除硼實驗條件的優化提供了有用信息,硼在該樹脂上的絡合吸附主要利用了其1:1配合物生成反應的機理。

(2) 厚層砂巖、軟弱基座和差異風化形成的凹巖腔,構造形成的剪切裂隙及組合共同構成了滑坡活動的控制性內因,裂隙水及其補給為控制性外因,施工棄土對滑坡活動具有積極推動作用。

圖2 不同pH值條件下硼酸-2,3-二羥基萘-6-磺酸鈉混合溶液中硼各存在形態的濃度Fig.2 pH dependent equilibrium concentrations of boron species for the solution of boric acid/borate and 2,3-dihydroxynaphthalene-6-sodium sulfonate

2.2 硼選擇性吸附劑制備條件的優化

樹脂上2,3-二羥基萘-6-磺酸鈉的擔載量對絡合除硼效果將產生直接影響。因此,分別考察了吸附時間、pH值和離子強度對Dowex 1-X8樹脂上2,3-羥基萘-6-磺酸鈉擔載量的影響。

吸附時間對2,3-羥基萘-6-磺酸鈉擔載量的影響如圖3所示。2,3-二羥基萘-6-磺酸鈉在樹脂上的分配比kd隨吸附時間的增長上升,至10h左右達到平衡。

圖3 吸附時間對樹脂上2,3-二羥基萘-6-磺酸鈉擔載量的影響T=25℃, C0(L) =0.001mol/L, pH=5.4, I =0.1mol/LFig.3 Effect of contact time on the amount of 2,3-dihydroxynaphtha lene-6-sodium sulfonate adsorbed on Dowex 1-X8 resin

溶液pH值對2,3-二羥基萘-6-磺酸鈉在Dowex 1-X8樹脂上擔載量的影響見圖4。在所研究的pH值范圍內,kd隨pH值的升高呈上升趨勢。這是由于2,3-二羥基萘-6-磺酸鈉的酚羥基隨pH值的升高開始解離,有利于在陰離子交換樹脂上結合。

圖4 溶液pH值對樹脂上2,3-二羥基萘-6-磺酸鈉擔載量的影響T=25℃, C0(L)=0.001mol/L, I=0.1mol/L, t=3hFig.4 Effect of pH value on the amount of 2,3-dihydroxynaphtha lene-6-sodium sulfonate adsorbed on the Dowex 1-X8 resin

圖5 離子強度對樹脂上2,3-二羥基萘-6-磺酸鈉擔載量的影響T=25℃, C0(L)=0.001mol/L, pH=7.4, t=3hFig.5 Effect of ionic strength on the amount of 2,3-dihydroxynaphtha lene-6-sodium sulfonate adsorbed on the Dowex 1-X8 resin

離子強度I對2,3-二羥基萘-6-磺酸鈉在Dowex 1-X8樹脂上擔載量的影響見圖5。可見,log kd隨log I的增大減小,二者呈良好的線性關系,線性相關系數R2=0.998。說明高離子強度環境不利于2,3-二羥基萘-6-磺酸鈉接枝吸附劑的制備。因此,實驗選擇了在pH 7.4、I=0.04mol/L的條件下制備2,3-二羥基萘-6-磺酸鈉接枝硼選擇性吸附劑。

2.3 2,3-二羥基萘-6-磺酸鈉接枝吸附劑上硼的吸附行為研究

(1)pH值對除硼效果的影響

圖6是不同pH條件下2,3-二羥基萘-6-磺酸鈉接枝吸附劑上硼吸附量的測試結果。由圖可知,接枝吸附劑上硼的吸附量先升后降,在pH值8.9處呈現峰值。這與圖2顯示的1:1配合物的理論分布一致。

為了比較,測試了相同條件下硼在Dowex1-X8樹脂上的吸附量(圖6)。可見,盡管Dowex1-X8樹脂對硼酸(硼酸根離子)也有一定的吸附能力,但對于飽和吸附2,3-二羥基萘-6-磺酸鈉的接枝樹脂而言,對硼吸附量的貢獻主要源于2,3-二羥基萘-6-磺酸鈉的作用。

圖6 溶液pH值對2,3-二羥基萘-6-磺酸鈉接枝吸附劑和Dowex 1-X8樹脂上硼吸附量的影響T=25℃, C0(B)=0.05mol/L, I=0.1mol/L, t=48hFig.6 Effect of pH value on boron adsorption capacities on the adsorbent presorbed 2,3-dihydroxynaphthalene-6-sodium sulfonate and Dowex 1-X8 resin

(2)吸附等溫實驗研究

溶液中硼酸初始濃度對2,3-二羥基萘-6-磺酸鈉接枝吸附劑上硼吸附量的影響試驗結果如圖7所示。圖7表明,樹脂上的硼吸附量隨硼酸初始濃度的增大上升,在C0=0.05mol/L處最大,此后呈下降趨勢。

圖7 硼酸初始濃度對吸附劑上硼吸附量的影響T=25℃, pH=8.9, I=0.1mol/L, t=48hFig.7 Effect of concentration of boric acid on boron adsorption capacity on the adsorbent presorbed 2,3-dihydroxynaphthalene-6-sodium sulfonate

分別用Langmuir吸附等溫式(10)和Freundlich吸附等溫式(11)對圖7的實驗結果擬合,結果見表1。

式中,Ce為平衡溶液中硼酸的濃度、qmax為飽和吸附量、KL為與溫度及吸附熱有關[12]的吸附系數、KF為Freundlich吸附系數、n為Freundlich常數。

表1 Langmuir和Freundlich吸附等溫模型的常數和擬合相關系數Table 1 Constants and correlation coefficients of Langmuir and Freundich models

由對Langmuir和Freundlich吸附模型擬合的相關系數(表1)可見,該吸附過程與上述兩種吸附模型均不相符,說明B(OH)4

-在2,3-二羥基萘-6-磺酸鈉接枝吸附劑上并非是簡單的化學式吸附。

根據圖7所得的硼最佳初始濃度(0.05mol/L)和制備樹脂中接枝2,3-二羥基萘-6-磺酸鈉的濃度(約0.03mol/L),可估算出pH8.9的條件下約生成1:1配合物0.03mol/L,即吸附劑上吸附量約為1.20mmol/g;而實際測得硼的最大吸附量為0.67mmol/g(圖6),約占絡合吸附的56%。

(3)吸附熱力學研究

式中:kd為硼酸在沒食子酸接枝樹脂上的分配比,ΔS0為標準熵變,ΔH0為標準焓變,R為氣體常數,T為溫度。

以lnkd對1/T作圖(圖8)并由直線的斜率和截距得ΔH0和ΔS0。將所得的ΔH0、ΔS0以及根據ΔG0=ΔH0–TΔS0計算的ΔG0列于表2。可見,在所研究的溫度范圍內,標準吉布斯自由能ΔG0均為負值,說明B(OH)4-在2,3-二羥基萘-6-磺酸鈉接枝吸附劑上的吸附配位反應是熱力學可行的自發過程。表2也顯示了吸附為放熱過程(ΔH0<0)。一般地,當ΔH0的絕對值為5~10kJ/mol時,吸附屬物理吸附;而當ΔH0的絕對值在30~70kJ/mol范圍內時,吸附屬于化學吸附[14]。

圖8 溫度對2,3-二羥基萘-6-磺酸鈉接枝吸附劑上硼吸附量的影響pH=8.9, C0(B)=0.05mol/L, I=0.1mol/L, t=24hFig.8 Effect of temperature on boron adsorption capacity on the adsor bent presorbed 2,3-dihydroxynaphthalene-6-sodium sulfonate

本實驗中ΔH0的絕對值介于10~30kJ/mol之間,表明該吸附過程可能既有化學吸附,也存在物理吸附。化學吸附即B(OH)4-與吸附劑中的多羥基官能團的絡合作用,B(OH)4-主要以1:1配合物的形式結合在2,3-二羥基萘-6-磺酸鈉接枝吸附劑上;而物理吸附可能是樹脂中的季胺離子R4N+與B(OH)4-之間產生的靜電吸引作用的結果。

表2 硼在2,3-二羥基萘-6-磺酸鈉接枝吸附劑上的吸附熱力學參數Table 2 Thermodynamic parameters for adsorption of boron on the adsorbent presorbed 2,3-dihydroxynaphthalene-6-sodium sulfonate

(4)吸附動力學研究

動力學研究結果表明(圖9),在硼的去除過程中,初期,吸附劑上硼的吸附量迅速升高,前6h的硼累計吸附量已達飽和吸附量的80%以上;之后吸附速率逐漸減小,至24小時左右趨于平衡。

圖9 硼在2,3-二羥基萘-6-磺酸鈉接枝吸附劑上的動力學曲線T=25℃, pH=8.9, C0(B)=0.5mol/L, I=0.1mol/LFig.9 Adsorption kinetics of boron on the adsorbent presorbed 2,3-dihydroxynaphthalene-6-sodium sulfonate

將圖9的實驗結果用準二級速率方程[15]擬合:

式中,qt為t時間的硼吸附量、k2為準二級速率常數。

擬合結果(圖10)顯示,B(OH)4-在該吸附劑上的吸附過程符合準二級速率方程(R2=0.999)。根據t/qt~t直線的截距,得k2=1.461g/mmol/h。由此可以推斷,吸附速率被化學吸附機理所控制,該化學吸附涉及到吸附劑與吸附質間的電子共用或電子轉移[12]。

圖10 硼在2,3-二羥基萘-6-磺酸鈉接枝吸附劑上的準二級速率方程擬合曲線T=25℃, pH=8.9, C0(B)=0.5mol/L, I=0.1mol/LFig.10 Pseudo-second-order adsorption kinetics of boron on the adsor bent presorbed 2,3-dihydroxynaphthalene-6-sodium sulfonate

References)

[1] Yoshimura K, Miyazaki Y, Ota F, et al. Complexation of boric acid with the N-methyl-D-glucamine group in solution and in crosslinked polymer[J]. Journal of the Chemical Society, Faraday Trans, 1998,94(5): 683-689.

[2] Shao C Y, Matsuoka S, Miyazaki Y, et al. Equilibrium and kinetic studies on the complexation of boric acid with chromotropic acid[J]. Journal of the Chemical Society, Dalton Trans, 2000: 3136-3142.

[3] Yoshimura K, Miyazaki Y, Sawada S, et al.11B NMR studies on complexation of borate with linear and crosslinked polysaccharides[J]. Journal of the Chemical Society, Faraday Trans, 1996,92(4): 651-656.

[4] Kabay N, Yilmaz I, Yamac S, et al. Removal and recovery of boron from geothermal wastewater by selective ion-exchange resins—II. Field tests[J]. Desalination, 2004, 167: 427-438.

[5] Jacob C. Seawater desalination: Boron removal by ion exchange technology[J]. Desalination, 2007, 205: 47-52.

[6] Ipeki Y, Kabay N, Yuksel M. Modeling of fixed bed column studies for removal of boron from geothermal water by selective chelating ion exchange resins[J]. Desalination, 2013, 310: 151-157.

[7] Hilal N, Kim G J, Somerfield C. Boron removal from saline water: A comprehensive review[J]. Desalination, 2011, 273: 23-25.

[8] Arias M F C, Bru L V I, Rico D P, et al. Comparison of ion exchange resins used in reduction of boron in desalinated water for human consumption[J]. Desalination, 2011, 278: 244-249.

[9] Dawber J G, Green S I E. An 11B nuclear magnetic resonance study of the reaction of the tetrahydroxyborate ion with polyhydroxy compounds[J]. Journal of the Chemical Society, Faraday Trans, 1986, 82: 3407-3413.

[10] Dawber J G, Green S I E, Dawber J C, et al. A polarimetric and 11B and 13C nuclear magnetic resonance study of the reaction of the tetrahydroxyborate ion with polyols and carbohydrates[J]. Journal of the Chemical Society, Faraday Trans, 1988,84(1): 41-56.

[11] Shao C Y, Matsuoka S, Miyazaki Y, et al. Studies on the complexation of boric acid with polyhydroxyl compounds[J]. Analytic Science, 2001,17(S1): 1475-1478.

[12] 張宏,張敬華. 生物吸附的熱力學平衡模型和動力學模型綜述[J].天中學刊,2009,24(5):19-22. Zhang H, Zhang J H. Review on thermodynamic equilibrium isotherm and kinetic model of biosorption[J]. Journal of Tianzhong, 2009,24(5):19-22.

[13] Yurdako? M, Seki Y, Karahan S, et al. Kinetic and thermodynamic studies of boron removal by Siral 5, Siral 40, and Siral 80[J]. Journal of Colloid & Interface Science, 2005, 286: 440-446

[14] Chen F X, Zhou C R, Li G P, et al. Thermodynamics and kinetics of glyphosate adsorption on resin D301[J]. Arabian Journal of Chemistry, 2012, 4,22.

[15] Ho Y S, Mckay G. Sorption of dye from aqueous solution by peat[J]. Chemical Engineering Journal, 1998, 70: 115-124.

Mechanism for the removal of trace amounts of boron by adsorbent presorbed 2,3-dihydroxynaphthalene-6-sodium sulfonate

YU Xiao-Feng1,2, GU Hua1,2, HU Jing-Jing3, WEN Xiao-Hua1,2, ZHANG Zhuo1,2, SHAO Chao-Ying3
(1. Shanghai Institute of Geological Survey, Shanghai 200072, China; 2. Shanghai Center for Intendance and Testing of Resources and Environment, Ministry of Land and Resources of China, Shanghai 200072,China; 3. College of Chemistry, Chemical Engineering and Biotechnology, Donghua University, Shanghai 201620, China)

The complexation of boric acid with 2,3-dihydroxynaphthalene-6-sodium sulfonate in aqueous solution was examined using11B NMR measurements. This confirms that boric acid reacts with 2,3-dihydroxynaphthalene-6-sodium sulfonate to form two complexes with 1:1 and 1:2 stoichiometries. The formation constants for these complexes are logβ1= 5.22 and logβ2= 6.10. A novel boron selective adsorbent was prepared by adsorbing 2,3-dihydroxynaphthalene-6-sodium sulfonate onto an ion-exchange resin. The loading amount of 2,3-dihydroxynaphthalene-6-sodium sulfonate on the adsorbents was 1.20 mmol/g. Bath experiments were conducted to evaluate the performance of the adsorbent for removal of boron from water, and the effects of pH and boron solution concentration, contact time, and bath temperature on the removal efficiency were investigated. The adsorption of boron onto the adsorbent conforms to the pseudo-second-order rate equation. It is a chemical and physical adsorption process, whereby B(OH)4?bonds to the adsorbent mainly in the form of 1:1 complex and electrostatic attraction. The adsorption rate is controlled by chemical adsorption mechanisms.

water environment; boron pollution; 2,3-dihydroxynaphthalene-6-sodium sulfonate; adsorbent; boric acid/ borate; mechanism of removal

P641.3

A

2095-1329(2015)01-0100-05

2014-08-08

2014-12-02

余肖峰(1986-),男,助理工程師,主要從事水文地球化學研究.

電子郵箱: yuxiaofeng86822@126.com

聯系電話: 021-56612395

中國地質調查局地質調查項目(12120113001900);上海市科學技術委員會科技專項(14231200504)

10.3969/j.issn.2095-1329.2015.01.023

猜你喜歡
化學
化學與日常生活
奇妙的化學
奇妙的化學
奇妙的化學
奇妙的化學
奇妙的化學
化學:我有我“浪漫”
化學:舉一反三,有效學習
考試周刊(2016年63期)2016-08-15 22:51:06
化學與健康
絢麗化學綻放
主站蜘蛛池模板: 国产女人18毛片水真多1| 亚瑟天堂久久一区二区影院| 国产尤物在线播放| 波多野结衣一区二区三区88| 国产99久久亚洲综合精品西瓜tv| 手机在线看片不卡中文字幕| 国产特级毛片aaaaaaa高清| 久久www视频| 欧美精品成人| 日本亚洲欧美在线| 国产成人精品2021欧美日韩| 天天综合天天综合| 亚洲综合精品第一页| 国产精品分类视频分类一区| 欧美日本在线一区二区三区| 伦精品一区二区三区视频| 久久人搡人人玩人妻精品| 亚洲欧美人成人让影院| 中文字幕日韩丝袜一区| 极品私人尤物在线精品首页 | 18禁高潮出水呻吟娇喘蜜芽| 国产69囗曝护士吞精在线视频| 91啦中文字幕| 精品无码人妻一区二区| 亚洲综合色区在线播放2019| 国产原创自拍不卡第一页| 国产va在线观看| 亚洲AⅤ综合在线欧美一区| 久久夜色精品| 欧美成人综合视频| 亚洲天天更新| 亚洲第一页在线观看| 国产精品区网红主播在线观看| 中文字幕在线欧美| 欧美日本在线播放| 58av国产精品| 国产91精品久久| 就去色综合| 亚洲综合九九| 国产真实二区一区在线亚洲| 久久这里只有精品免费| 欧美日韩午夜| 九九香蕉视频| 久久精品丝袜高跟鞋| 精品久久综合1区2区3区激情| 成人在线综合| 亚洲视频二| 国产免费a级片| 欧美成一级| 久久国产精品影院| 国内精品九九久久久精品| 2024av在线无码中文最新| 亚洲AV无码久久精品色欲| 国产在线91在线电影| 无码人中文字幕| 伊人91在线| 成人国产精品视频频| 国产91麻豆视频| 日韩精品亚洲一区中文字幕| AⅤ色综合久久天堂AV色综合 | A级毛片无码久久精品免费| 国产成人亚洲无码淙合青草| 老司机精品99在线播放| 久久综合亚洲色一区二区三区| 99九九成人免费视频精品 | 亚洲日韩精品综合在线一区二区| 5555国产在线观看| 亚洲视频三级| 91蝌蚪视频在线观看| 青青草原国产一区二区| 真实国产精品vr专区| 国产人妖视频一区在线观看| 日韩欧美视频第一区在线观看| 广东一级毛片| 国产免费网址| 国产综合另类小说色区色噜噜| 波多野结衣亚洲一区| 98超碰在线观看| 67194成是人免费无码| 丁香婷婷激情综合激情| 伊在人亚洲香蕉精品播放| 国产成人精品三级|