離子交換樹脂對D-甘油酸的吸附熱力學和動力學

王寶貝，蒲洋，林麗芹，李清彪,，盧英華

?

離子交換樹脂對D-甘油酸的吸附熱力學和動力學

王寶貝1，蒲洋2，林麗芹1，李清彪1,2，盧英華2

（1泉州師范學院海洋與食品學院，福建泉州 362000；2廈門大學化學與化工學院，福建廈門 361005）

D-甘油酸是一種重要的甘油衍生物，具有解酒護肝的功能。關于D-甘油酸生產(chǎn)方法的研究已有不少，但對其分離方法及分離機理的研究卻鮮見報道。通過靜態(tài)吸附實驗，研究了D-甘油酸在201×7陰離子交換樹脂上的等溫熱力學和動力學特性。結果表明，D-甘油酸在201×7陰離子交換樹脂上的最大平衡吸附容量隨pH的增加而降低，其吸附等溫線符合Freundlich模型。在293～308 K下，吸附焓變?yōu)?4.77 kJ·mol?1，表明該吸附過程為吸熱過程。升高溫度有利于提高吸附速率，但對最大平衡吸附容量影響不大。同時，采用動邊界模型描述D-甘油酸在該樹脂上的交換行為，分別考察了料液濃度、樹脂粒徑和溫度對交換過程的影響。交換過程的吸附速率隨D-甘油酸濃度和溫度的增加而增大，但隨樹脂粒徑的增大而減小。研究表明該離子交換過程的速率控制步驟為顆粒擴散過程，交換過程的反應速率常數(shù)0為1.22×10?3，反應級數(shù)為0.631，表觀活化能a為14.90 kJ·mol?1，并得到了動力學總方程。

D-甘油酸；離子交換；吸附；動力學；熱力學

引 言

甘油酸（glyceric acid，化學名2, 3-dihydroxypropionic acid）是一種重要的甘油衍生物。它在化工合成、精細化工、醫(yī)藥等行業(yè)中具有廣泛的應用[1-2]。其右旋產(chǎn)品，D-甘油酸，具有解酒護肝的功效[3]，是一種重要的醫(yī)藥原料。雖然化學合成的甘油酸產(chǎn)量較高，但產(chǎn)品中D型和L型各占50%，而其中只有D-甘油酸才具有生物學活性。當前，藥用D-甘油酸主要采用生物發(fā)酵法生產(chǎn)[1, 4]，再利用電滲析法進行分離[5]。雖然該法能有效從發(fā)酵液中分離得到純度較高的D-甘油酸，但其生產(chǎn)效率有限，且使用過程中發(fā)酵液中的細胞常黏附于膜上，大大縮短了膜的使用壽命。此外，電滲析設備成本高、能耗大，不利于規(guī)?；a(chǎn)使用[6]。相比之下，離子交換吸附分離方法設備投資少、工藝簡單、污染小、能耗小，并且易于中試放大，因此工業(yè)生產(chǎn)中常用該法分離純化各種羧酸[7-8]、弱酸[9-10]、金屬離子等產(chǎn)品[11-12]。D-甘油酸具有弱酸性，在堿性條件下可以完全電離，電離出的甘油酸根離子可與陰離子交換樹脂上的活性基團發(fā)生離子交換，因此具備了采用離子交換分離的條件。有研究成功采用陰離子交換樹脂（IRA400和IRN78）從裂解液中分離出D-甘油酸和乳酸[13]。然而，當前對采用離子交換技術分離D-甘油酸的研究較少，對離子交換吸附D-甘油酸的動力學和熱力學特性的研究鮮有報道。為此，研究了D-甘油酸在201×7陰離子交換樹脂（與IRA400性質(zhì)相近）上的吸附動力學和熱力學，為今后采用離子交換法分離D-甘油酸的生產(chǎn)實踐提供理論基礎。

1 實驗材料和方法

1.1 試劑

D-甘油酸標準品購自Sigma公司；201×7型陰離子樹脂購自西隴化工股份有限公司；其他試劑均為分析純試劑。

1.2 樹脂預處理

樹脂先用2 mol·L-1的NaOH和HCl溶液反復浸泡3次后，用2 mol·L-1的NaOH轉成OH-型，最終用去離子水洗至中性，然后于60℃干燥2 h備用。

1.3 吸附等溫線的測定

準確稱取0.5000 g預處理樹脂，分別在pH為8、9、10、11和12條件下，加入15 ml不同濃度的D-甘油酸溶液，于25℃下100 r·min-1振蕩24 h，確保達到吸附平衡。檢測上清液中D-甘油酸濃度e，按式（1）計算平衡吸附量e，作e-e圖，得吸附等溫線。

1.4 吸附動力學實驗

本實驗采用間歇法，即準確稱取0.5000 g預處理的樹脂于錐形瓶中，加入已知濃度的D-甘油酸溶液，于恒溫恒速搖床中進行反應，定時測定溶液中的D-甘油酸濃度。分別在不同濃度、不同粒徑和不同溫度條件下測定溶液濃度隨時間的變化關系，得到相關的動力學數(shù)據(jù)。第次取樣時D-甘油酸的交換度為

式中，0為初始D-甘油酸濃度，g·L?1；C為第次取樣時溶液中D-甘油酸濃度，g·L?1；為料液體積，L；m為樹脂飽和吸附量，g·kg?1；為樹脂質(zhì)量，kg。

1.5 檢測方法

D-甘油酸采用高效液相色譜（Agilent 1200）分析檢測，以0.5 mmol·L-1H2SO4為流動相，流速為0.8 ml·min?1，分離柱為Shodex SH1011[4]。

2 實驗結果與討論

2.1 吸附等溫線

2.1.1 不同pH下的吸附等溫線 在離子交換分離過程中，料液溶質(zhì)濃度越高，溶質(zhì)與樹脂接觸的機會越多，就越容易進入樹脂的網(wǎng)孔內(nèi)進行離子交換。這可在一定程度上對樹脂的交換容量和交換速度產(chǎn)生影響。但料液溶質(zhì)濃度太高也可能引起樹脂表面及內(nèi)部交聯(lián)網(wǎng)孔收縮，進而影響離子的進入，從而降低樹脂的交換容量。謝友坪等[14]、程亮等[12]通過對色氨酸、鎘離子的離子交換分離研究均發(fā)現(xiàn)料液溶質(zhì)濃度是影響樹脂吸附容量、吸附速率的一個重要因素。此外，Caetano等[9]、Xie等[15]在對苯酚、色氨酸的離子交換分離機理的研究中發(fā)現(xiàn)料液pH會通過改變?nèi)苜|(zhì)的解離平衡和溶質(zhì)在樹脂上的離子交換吸附平衡，影響吸附選擇性及樹脂交換容量。鑒于料液濃度和pH在離子交換過程中的重要作用，本研究分別考察了不同料液濃度和不同料液pH對D-甘油酸的離子交換吸附平衡的影響。由圖1可見，在該樹脂上，D-甘油酸的平衡吸附量隨著料液濃度的增加而逐漸增大。這主要是由于提高料液中D-甘油酸的濃度，增加了吸附過程的推動力，從而促進了吸附過程。此外，隨著料液pH的降低，樹脂的最大平衡吸附量逐漸增加。當pH為8時，樹脂的最大平衡吸附量達到最大。D-甘油酸屬于弱酸，堿性條件有利于D-甘油酸的解離，進而促進離子交換吸附。然而，本實驗所用樹脂為OH-型，OH-是離子交換吸附過程的產(chǎn)物之一。因此，太高的料液pH反而會使離子交換吸附朝反方向進行，進而降低最大平衡吸附量。本研究分別采用Langmuir[式(3)]、Freundlich[式(4)]，對以上靜態(tài)吸附數(shù)據(jù)進行擬合

(4)

兩方程的擬合結果如圖1所示，回歸結果及其相關系數(shù)見表1。由表可知，F(xiàn)reundlich方程回歸曲線的相關系數(shù)2均比Langmuir方程高。在所考察的pH中，F(xiàn)reundlich模型對D-甘油酸的靜態(tài)吸附都有很好的擬合效果（2>0.98）。不同pH條件下，F(xiàn)reundlich模型中的1/均介于0與1之間，表明樹脂201×7對D-甘油酸具有良好的吸附效果。相反，Langmuir模型并不適用于所有pH條件：在弱堿條件下，D-甘油酸的靜態(tài)吸附并不滿足Langmuir模型（2<0.95），但是隨著pH的升高，Langmuir模型的擬合度越來越高。在pH=12時，其擬合度達到2=0.99。

表1 Langmuir和Freundlich方程回歸結果

2.1.2 不同溫度下的靜態(tài)吸附等溫線 在293、298、303和308 K時，分別研究了樹脂201×7對D-甘油酸的吸附等溫線，并采用擬合效果較好的Freundlich模型擬合，結果如圖2、表2所示。由圖2可見，隨著溫度的升高，該樹脂對D-甘油酸的吸附速率越快，但溫度對最大平衡吸附量影響不大。由表2可知，F(xiàn)reundlich模型中的1/隨溫度升高而減小，這表明溫度越高樹脂對D-甘油酸的吸附效果越好。

表2 Freundlich方程回歸結果

2.2 靜態(tài)吸附D-甘油酸的熱力學參數(shù)

樹脂201×7對D-甘油酸的吸附規(guī)律可以通過計算吸附自由能變（Δ）、吸附焓變（Δ）、吸附熵變（Δ）等熱力學函數(shù)加以解釋。根據(jù)Freundlich經(jīng)驗模型的吸附系數(shù)與溫度、吸附焓存在如下關系[16]

式中，為熱力學溫度，K；0為常數(shù)；為摩爾氣體常數(shù)，=8.314 J·mol?1·K?1；Δ為吸附焓變，kJ·mol?1。

以ln對1/作圖，進行線性擬合，通過斜率求得Δ。此外，吸附自由能變Δ可由式(6)得到[17]

Δ= ?(6)

式中，為Freundlich常數(shù)；Δ為吸附自由能變，kJ·mol?1。

吸附熵變Δ（kJ·mol?1·K?1），可根據(jù)Gibbs-Helmhotz方程計算

Δ=(ΔΔ)/(7)

由式(5)～式(7)計算出熱力學函數(shù)Δ、Δ、Δ，結果如表3所示。

表3 樹脂201×7吸附D-甘油酸的熱力學函數(shù)值

由表3可以看出，D-甘油酸在樹脂上的吸附焓變Δ>0，表明體系中D-甘油酸的吸附產(chǎn)生的放熱效應小于溶劑水的解吸等帶來的吸熱效應。因此，D-甘油酸在該樹脂上的吸附過程是吸熱過程，適當提高溫度有利于增大吸附速率。該分析結果與圖2實驗結果一致，即溫度越高樹脂對D-甘油酸的吸附效果越好。Δ>0表明樹脂吸附D-甘油酸的過程是熵增的過程，這說明達到平衡時，該吸附體系的自由度升高。Δ<0表明被吸附到樹脂上的D-甘油酸與溶液中的D-甘油酸相比，自由運動得到更大限制，且該吸附過程是自發(fā)進行的。

2.3 D-甘油酸與樹脂作用的動力學研究

為了研究樹脂201×7吸附D-甘油酸的動力學過程，并探知其速率控制步驟，采用動邊界模型（shell progressive model, SPM）擬合D-甘油酸的動力學吸附數(shù)據(jù)。該模型將離子交換吸附分為3個步驟：①離子由溶液經(jīng)液膜擴散到樹脂表面（即液膜擴散過程）；②離子由樹脂表面向樹脂內(nèi)部擴散（即顆粒擴散過程）；③離子在樹脂內(nèi)活性基團的位置發(fā)生交換反應（即化學反應過程）。因此，離子交換的動力學過程受到以上3個步驟的共同影響，其中速率最慢的一步為該離子交換過程的速率控制步驟。當液膜擴散、顆粒擴散和化學反應分別為速率控制步驟時，其模型方程分別為[18]：

膜擴散

顆粒擴散

(9)

化學反應

分別以、3?3(1?)2/3?2和1?(1?)1/3與進行線性擬合。擬合結果如圖3所示，上述各方程的線性相關系數(shù)R分別為0.926、0.985、0.946，即3?3(1?)2/3?2與的擬合線性關系最好，故在樹脂201×7吸附D-甘油酸的過程中，顆粒擴散為速率控制步驟。

圖3 吸附過程速率控制步驟的擬合

Fig.3 Estimation of controlling step of adsorption process

2.3.1 料液濃度對交換過程的影響 料液濃度對離子交換吸附過程有很大的影響，該樹脂對D-甘油酸的平衡吸附量隨著料液濃度的升高而逐漸增加（圖1）。本研究考察了3種D-甘油酸濃度對離子交換過程的影響，通過動邊界模型顆粒擴散控制方程式（9）擬合數(shù)據(jù)，結果如圖4所示，線性回歸結果見表4。從表4可見不同料液濃度下，各直線斜率，即表觀速率常數(shù)（×10?4）分別為：6.128、4.475、2.883。

表4 不同料液濃度的顆粒擴散方程

由于表觀速率常數(shù)與濃度冪函數(shù)呈正比，即

ln=ln0(11)

其中是D-甘油酸濃度的反應級數(shù)，對ln和ln0進行線性擬合，可得ln= ?9.2820.631ln0（2=0.99998），反應級數(shù)為0.631。

2.3.2 樹脂粒徑對交換過程的影響 D-甘油酸在樹脂中的擴散速率受到樹脂粒徑的直接影響，而上述研究表面顆粒擴散是該離子交換過程的速率控制步驟，因此樹脂粒徑對離子交換速率的影響至關重要。本研究考察了兩種不同粒徑的樹脂對D-甘油酸吸附過程的影響。所得動力學數(shù)據(jù)用顆粒擴散控制方程式（3）進行擬合，結果如圖5所示，線性回歸結果見表5。

表5 不同粒徑的顆粒擴散方程

由表可知，粒徑較小的樹脂對D-甘油酸的吸附速率更快。這是由于粒徑較小的樹脂，溶液中D-甘油酸根離子通過樹脂表面進入內(nèi)部未反應區(qū)的路徑也較短，從而加快了樹脂對D-甘油酸的吸附速率。同時，通過方程式（6）可以發(fā)現(xiàn)，表觀速率常數(shù)和樹脂粒徑平方的倒數(shù)呈正比。為驗證該關系，將圖6中的數(shù)據(jù)用式（12）[14]處理

結果如圖6所示。當樹脂直徑分別為0.4、0.5 mm時，擬合直線的2為0.9920。這說明與1/2呈較好的正比關系，進一步表明D-甘油酸的吸附過程受顆粒內(nèi)擴散控制。

圖6 不同樹脂粒徑下[3?3(1?)2/3?2]2與的關系

Fig.6 Relation of [3?3(1?)2/3?2]2andat different particle diameters

2.3.3 溫度對交換過程的影響 考察4種溫度下樹脂對D-甘油酸的動態(tài)吸附過程，通過顆粒內(nèi)擴散模型對數(shù)據(jù)進行擬合，3?3(1?)2/3?2與呈線性關系，如圖7所示，所得線性回歸見表6。由表6可知，溫度升高，擴散速率加快，有利于吸附過程的進行。從表中各直線斜率可得到不同溫度下D-甘油酸交換過程的表觀速率常數(shù)，將其代入Arrhenius方程

ln= ?a/+(13)

以ln對1/作圖，擬合得到方程：ln= ?1792.6/?1.407（R=0.9706）。由此，可得D-甘油酸吸附過程的表觀活化能a為14.90 kJ·mol?1。

表6 不同溫度的顆粒擴散方程

2.3.4 離子交換過程的動力學方程 通過研究D-甘油酸濃度、反應溫度以及樹脂粒徑對離子交換過程的影響，以及通過顆粒擴散方程得到一系列的表觀速率常數(shù)，可表示為D-甘油酸濃度、溫度以及樹脂粒徑的函數(shù)關系[19]

求得D-甘油酸吸附過程的速率常數(shù)0為1.22×10?3，因此該吸附過程的動力學總方程式為

(15)

3 結 論

（1）D-甘油酸在201×7型陰離子交換樹脂上的吸附等溫線符合Freundlich等溫方程，而Langmuir等溫方程只有在pH為12時擬合效果才較好。

（2）離子交換樹脂201×7吸附D-甘油酸的熱力學研究表明：靜態(tài)平衡吸附等溫線較好地符合Freundlich模型，其中Δ>0，Δ<0，Δ>0，表明樹脂對D-甘油酸的吸附是一個吸熱的過程，且是自發(fā)進行的，吸附過程是熵增加的過程，升溫有利于提高吸附速率。

（3）動力學研究表明，離子交換過程的速率控制步驟為顆粒擴散，其表觀活化能為14.90 kJ·mol?1，D-甘油酸濃度的反應級數(shù)為0.631，表觀速率常數(shù)和樹脂粒徑平方呈反比，并推導出該離子交換過程的動力學總方程為。

符 號 說 明

a——反應級數(shù) Ce——吸附平衡后溶液中D-甘油酸的濃度，g·L?1 Ci——第i次取樣時，溶液中D-甘油酸的濃度，g·L?1 C0——料液中初始D-甘油酸的濃度，g·L?1 D——顆粒有效擴散系數(shù)，m2·s?1 Ea——表觀活化能，kJ·mol?1 F——交換度 ΔG——吸附自由能變，kJ·mol?1 ΔH——吸附焓變，kJ·mol?1 K——Freundlich吸附平衡常數(shù) KL——Langmuir常數(shù)，L·g?1 k——表觀速率常數(shù) m——樹脂質(zhì)量，g n——Freundlich方程的特征常數(shù) qe——樹脂平衡吸附量，g·g?1 qm——樹脂飽和吸附量，g·g?1 R——氣體常數(shù)，J·mol?1·K?1 r——樹脂顆粒半徑，m ΔS——吸附熵變，kJ·mol?1·K?1 T——熱力學溫度，K V——料液體積，L

References

[1] HABE H, FUKUOKA T, KITAMOTO D,. Biotechnological production of D-glyceric acid and its application [J]. Applied Microbiology and Biotechnology, 2009, 84(3): 445-452.

[2] FUKUOKA T, IKEDA S, HABE H,. Synthesis and interfacial properties of monoacyl glyceric acids as a new class of green surfactants [J]. Journal of Oleo Science, 2012, 61(16): 343-348.

[3] HABE H, SATO H, FUKUOKA T,. Effect of glyceric acid calcium salt on the viability of ethanol-dosed gastric cells [J]. Journal of Oleo Science, 2011, 60(11): 585-590.

[4] HABE H, FUKUOKA T, KITAMOTO D,. Biotransformation of glycerol to D-glyceric acid by acetobacter tropicalis [J]. Applied Microbiology and Biotechnology, 2009, 81(6): 1033-1039.

[5] HABE H, FUKUOKA T, KITAMOTO, D,. Application of electrodialysis to glycerate recovery from a glycerol containing model solution and culture broth [J]. Journal of Bioscience and Bioengineering, 2009, 107(4): 425-428.

[6] MOLDES A B, ALONSO J L, PARAJO J C. Recovery of lactic acid from simultaneous saccharification and fermentation media using anion exchange resins [J]. Bioprocess Biosystems Engineering, 2003, 25(6): 357-363.

[7] LOPEZ-GARZON C S, STRAATHOF A J. Recovery of carboxylic acids produced by fermentation [J]. Biotechnology Advances, 2014, 32(5): 873-904.

[8] ATAEI S A, VASHEGHANI-FARAHANI E.separation of lactic acid from fermentation broth using ion exchange resins [J]. Journal of Industrial Microbiology and Biotechnology, 2008, 35(11): 1229-1233.

[9] CAETANO M, VALDERRAMA C, FARRAN A,. Phenol removal from aqueous solution by adsorption and ion exchange mechanisms onto polymeric resins [J]. Journal of Colloid and Interface Science, 2009, 338(2): 402-409.

[10] CHEN K, LYU H, HAO S,. Separation of phenolic compounds with modified adsorption resin from aqueous phase products of hydrothermal liquefaction of rice straw [J]. Bioresource Technology, 2015, 182(1): 160-168.

[11] 張金麗, 孫道華, 景孝廉, 等. Au（Ⅲ）離子在黑曲霉菌上的吸附熱力學和動力學特性[J]. 化工學報, 2013, 64(4): 1283-1292. ZHANG J L, SUN D H, JING X L,. Thermodynamic and kinetic characteristic of Au(Ⅲ) biosorption onto[J]. CIESC Journal, 2013, 64(4): 1283-1292.

[12] 程亮, 侯翠紅, 徐麗, 等. 納米腐殖酸動態(tài)吸附廢水中鎘離子及其洗脫特性[J]. 化工學報, 2016, 67(4): 1348-1356. CHEN L, HOU C H, XU L,. Dynamic adsorption and de-sorption characteristics of wastewater containing cadmium ion on nanoscale humic acid [J]. CIESC Journal, 2016, 67(4): 1348-1356.

[13] JACKSON J E, MILLER D J, MARINCEAN S. Process for the preparation of lactic acid and glyceric acid: US 7459584 (B) [P]. 2008-12-02.

[14] 謝友坪, 敬科舉, 姚傳義, 等. 離子交換樹脂吸附L-色氨酸的動力學研究[J]. 高校化工學報, 2012, 26(4): 710-715.XIE Y P, JING K J, YAO C Y,. Study on the kinetics of L-tryptophan adsorption by ion exchange resin [J]. Journal of Chemical Engineering of Chinese Universities, 2012, 26(4): 710-715.

[15] XIE Y P, JING K J, LU Y H. Kinetics, equilibrium and thermodynamic studies of L-tryptophan adsorption using a cation exchange resin [J] Chemical Engineering Journal, 2011, 171(3): 1227-1233.

[16] 程亮, 侯翠紅, 劉國際, 等. 納米腐殖酸對鎘離子的吸附熱力學及動力學 [J]. 高?；瘜W工程學報, 2015, 29(1): 72-77. CHENG L, HOU C H, LIU G J,. Thermodynamics and kinetics of cadmium ions adsorption by nanoscale humic acid [J]. Journal of Chemical Engineering of Chinese Universities, 2015, 29(1): 72-77.

[17] MARK S, PINHUA R, IRENE L. Effects of hardness and alkalinity on the removal of arsenic (Ⅴ) from humic acid-deficient and humic acid-rich groundwater by zero-valent iron [J]. Water Research, 2009, 43(17): 4296-4304.

[18] 楊莉麗, 康海彥, 李娜, 等. 離子交換樹脂吸附鎘的動力學研究[J]. 離子交換與吸附, 2004, 20(2): 138-143. YANG L L, KANG H Y, LI N,. Study on the kinetics of adsorption of cadmium by ion-exchange resin [J]. Ion Exchange and Adsorption, 2004, 20(2): 138-143.

[19] 王永利, 馬曉建, 常春. 離子交換樹脂吸附乙酰丙酸的動力學研究[J]. 離子交換與吸附, 2007, 23(4): 330-336. WANG Y L, MA X J, CHANG C. Study on the kinetics of adsorption of levulinic acid by ion-exchange resin [J]. Ion Exchange and Adsorption, 2007, 23(4): 330-336.

Thermodynamics and kinetics of D-glyceric acid adsorption on ion exchange resin

WANG Baobei1, PU Yang2, LIN Liqin1, LI Qingbiao1, 2, LU Yinghua2

(1College of Oceanography and Food Science, Quanzhou Normal University, Quanzhou 362000, Fujian, China;2Department of Chemical and Biochemical Engineering, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, Fujian, China)

D-glyceric acid, one of the promising glycerol derivatives, has various biological functions such as accelerating ethanol and acetaldehyde oxidation. Plenty studies on the synthesis but few studies on downstream processes to recover D-glyceric acid were available. A static equilibrium adsorption was performed to study isothermal thermodynamics and kinetics of D-glyceric acid adsorption onto 201×7 strong basic anion exchange resin. The results showed that the adsorption process followed the Freundlich isothermal model at different pH and temperature conditions, which the maximum adsorption capacity decreased with increasing initial pH. The change of adsorption enthalpy was 14.77 kJ·mol?1at 293—308 K, suggesting an endo-thermic adsorption process. Raising temperature improved adsorption rate with limited influence on equilibrium adsorption capacity. The Dynamic Boundary Model was used to investigate adsorption kinetics with a focus on effects of D-glyceric acid concentration, resin particle diameter and adsorption temperature. The results suggested the particle diffusion was rate-limiting step of the adsorption process. The adsorption rate increased with the increase of D-glyceric acid concentration and temperature, but decreased with the increase of resin particle diameter. The general kinetics equation was obtained with the rate constant of exchange (0), reaction order () and the apparent activation energy of reaction (a) as following0=1.22×10?3,=0.631,a=14.90 kJ·mol?1.

D-glyceric acid; ion exchange; adsorption; kinetics; thermodynamics

2016-05-18.

Prof. LU Yinghua, ylu@xmu.edu.cn

10.11949/j.issn.0438-1157.20160685

TQ 028.3

A

0438—1157（2016）11—4671—07

王寶貝（1982—），女，博士，講師。

福建省高校產(chǎn)學合作項目（2015N5006）；福建省教育廳中青年教師教育科研項目（JA15393）；泉州師范學院青年博士預研基金項目（2015QBKJ04）。

2016-05-18收到初稿，2016-07-12收到修改稿。

聯(lián)系人：盧英華。

supported by the Industry-University Cooperation Project of Fujian Province (2015N5006), the Young PhD Pre-research Fund of Quanzhou Normal University (2015QBKJ04), and the Education Research Fund for Young and Middle-aged Teachers in Fujian Province (JA15393).