蔣志國，李 斌，*，王燕華，張海德（.海南大學食品學院，海南?？?70228；2.海南大學信息學院，海南?？?70228）

菠蘿蜜果皮多酚超聲微波協(xié)同提取工藝優(yōu)化及抗氧化活性研究

蔣志國1，李 斌1，*，王燕華2，*，張海德1
（1.海南大學食品學院，海南?？?70228；2.海南大學信息學院，海南?？?70228）

采用超聲-微波協(xié)同提取技術（UMAE）對菠蘿蜜果皮中多酚的提取工藝進行優(yōu)化，并對抗氧化活性進行了評價。以單因素實驗為基礎，根據(jù)Box-Behnken中心組合設計原理，選取乙醇體積分數(shù)、料液比、微波功率和微波時間4因素3水平進行響應曲面分析，建立多酚得率的二次多項數(shù)學模型，分析各因素的顯著性和交互作用，得到多酚提取工藝的最佳條件為：乙醇體積分數(shù)70%、料液比1∶40、微波功率75 W、微波時間12 min，多酚得率為7.19 mg/g。在該條件下，超聲-微波協(xié)同提取方法提取效率優(yōu)于傳統(tǒng)水浴回流法（1.04 mg/g）、微波輔助法（5.23 mg/g）和超聲輔助法（5.89 mg/g）?？寡趸钚匝芯勘砻?，菠蘿蜜果皮多酚提取物對DPPH自由基和ABTS自由基均有較強的清除能力，呈量效關系，其EC50值分別為101.39 μg/mL和106.60 μg/mL，表明多酚是菠蘿蜜果皮抗氧化活性的物質基礎。

菠蘿蜜果皮，多酚，超聲-微波協(xié)同萃取，響應面法，抗氧化

表1 超聲波-微波協(xié)同提取波菠蘿蜜多酚物質單因素實驗Table1 Univariate experimrnts of UMAE for extracting polyphenols from jackfruit peel

菠蘿蜜（Artocarpus heterophyllus Lam.）屬于桑科木菠蘿屬，是世界著名的熱帶水果之一，大量種植于東南亞、印度大陸、中國南部、中非和拉丁美洲［1］，國內以海南種植最多、最集中。

菠蘿蜜果肉爽脆、蜜甜、芳香濃郁，營養(yǎng)成分十分豐富，因此常被加工成果干、果汁、果酒等多種產(chǎn)品遠銷各地。然而，研究者對占果實總重約42%的菠蘿蜜下腳料果皮的研究卻甚少，主要集中在果膠的提取和利用上［2-3］，但對果皮中其他活性成分的研究鮮有報道，特別是多酚類化合物。由于近年來發(fā)現(xiàn)多酚類化合物能夠延緩衰老、抑制人類腫瘤細胞的發(fā)生和發(fā)展［4］、降血糖和降血脂［5］、抗病毒和抗氧化等作用［6］，因而激發(fā)了人們研究和開發(fā)利用多酚類化合物的極大興趣。

目前，對多酚化合物的提取常采用水浴振蕩法［7］、超聲波輔助法和微波輔助法。超聲-微波協(xié)同萃取法是固體萃取中超聲波法及微波法發(fā)展而來的新技術。此法充分利用超聲波振動的空化作用及微波的高能作用，克服了常規(guī)超聲波和微波的萃取不足，即超聲波產(chǎn)生的振蕩和攪拌作用能有效彌補微波傳熱、傳質不均等缺陷，而微波極佳的熱效應能有效地彌補超聲波產(chǎn)熱不足的問題，實現(xiàn)了優(yōu)勢互補［8-10］。國內外很多學者采用超聲-微波協(xié)同萃取各種原料中的功能因子，與傳統(tǒng)提取方法相比，均發(fā)現(xiàn)萃取時間更短、得率更高、效率更好［11-12］。

目前，關于超聲-微波協(xié)同萃取法提取菠蘿蜜果皮多酚還未見報道。因此，本文以菠蘿蜜果皮為原料，利用響應曲面法對菠蘿蜜果皮中多酚的提取工藝進行優(yōu)化，并對抗氧化活性進行了初步評價，為實現(xiàn)菠蘿蜜的綜合利用，提高產(chǎn)品附加值，充分發(fā)揮海南省菠蘿蜜的資源優(yōu)勢，并為尋找開發(fā)菠蘿蜜的新途徑提供一定的理論基礎，為其工業(yè)化生產(chǎn)提供理論指導。

1 材料與方法

1.1 材料與儀器

菠蘿蜜 購于?？谑腥f福新村農貿市場；2，2-二苯基-1-苦基肼（DPPH自由基） Sigma公司；2，2′-聯(lián)氮雙（3-乙基苯并噻唑啉-6-磺酸）二銨鹽（ABTS自由基） Sigma公司；沒食子酸 Sigma公司；過硫酸鉀 國藥集團；福林酚試劑 國藥集團化學試劑有限公司；其他試劑均為分析純。

CW-2000微波超聲協(xié)同萃取儀 上海新拓分析儀器科技有限公司；FW177中草藥粉碎機 成都一科儀器設備有限公司；AL204電子天平 梅特勒-托利多儀器（上海）有限公司；722型可見分光光度計 上海奧譜勒儀器有限公司；SHZ-DⅢ循環(huán)水式真空泵 鞏義市予華儀器有限責任公司。

1.2 實驗方法

1.2.1 原料的預處理 新鮮的菠蘿蜜果肉和果皮分離后，取菠蘿蜜果皮清洗干凈，并于60℃烘干、粉碎、過40目篩備用。

1.2.2 總酚含量測定 以沒食子酸（GA）為對照品，參考Fang等采用的Folin-酚法［13］對菠蘿蜜果皮多酚提取物及純化物中的總酚含量進行測定，以GA當量表示樣品的總酚含量。GA標準曲線：Y=0.0729X+0.0008，R2=0.9985，其中X為GA質量濃度（μg/mL），Y為吸光度值。

稱取1.0 g菠蘿蜜果皮粉末，超聲-微波協(xié)同萃取得多酚提取液，經(jīng)抽濾后用蒸餾水定容至100 mL，即為樣品溶液。取1.5 mL樣品溶液于25 mL比色管中，分別加入1 mL福林-酚試劑和2 mL質量分數(shù)為10% 的Na2CO3溶液，然后加入蒸餾水定容至25 mL，室溫避光靜置40 min，在760 nm處測定吸光值。將吸光值代入標準曲線方程和式（1）計算樣品中多酚得率。

式中，Y為多酚得率（mg/g）；C為沒食子酸含量（mg/mL）；V為提取液體積（mL）；N為稀釋倍數(shù)；M為樣品質量（g）。

1.2.3 超聲-微波協(xié)同萃取菠蘿蜜果皮多酚單因素實驗 稱取1.0 g菠蘿蜜果皮粉末，進行單因素實驗（表1），研究各因素對菠蘿蜜果皮多酚得率的影響。

1.2.4 響應面實驗設計 在單因素實驗基礎上，根據(jù)Box-Benhnken中心組合實驗設計原理，安排4因素3水平實驗，確定超聲-微波協(xié)同萃取菠蘿蜜果皮多酚的最佳工藝條件，實驗因子及水平如表2。

表2 Box-Benhnken實驗設計因素水平Table2 Factors and levels table of Box-Benhnken experimental design

1.2.5 抗氧化活性測定

1.2.5.1 DPPH自由基清除率測定 參照Y Liu等［14］的方法，略有改動。取3 mL樣品溶液、0.1 mmol/L DPPH乙醇溶液1 mL加入具塞試管中，混合振蕩搖勻，室溫避光放置30 min，以無水乙醇為空白在517 nm測定其吸光度Ai；取樣品液3 mL、無水乙醇1 mL同上方法測吸光度Aj；取無水乙醇3 mL、DPPH液1 mL同上方法測吸光度A0。以BHT為陽性對照，每組實驗均平行3次。樣品液對DPPH的清除率（%）=［1-（Ai-Aj）/A0］×100。

1.2.5.2 ABTS自由基清除率測定 采用Gan CY等［15］的方法，略有改動。7.0 mmol/L的ABTS水溶液和2.45 mmol/L的過硫酸鉀溶液1∶1（v∶v）混合后，避光反應16 h，用95%乙醇稀釋使其在734 nm處吸光值為0.700±0.005（95%乙醇為對照），作為工作液備用。取樣品液0.4 mL、ABTS工作液1.6 mL于具塞試管，混合振蕩搖勻，在黑暗中靜置6 min，以95%乙醇為空白在734 nm處測定吸光度Ai；取樣品液0.4 mL、95%乙醇1.6 mL吸光度Aj；取95%乙醇0.4 mL、ABTS液1.6 mL吸光度A0。以BHT為陽性對照，每組實驗均平行3次。樣品液對ABTS的清除率（%）=［1-（Ai-Aj）/A0］×100。

1.2.6 超聲-微波協(xié)同萃取高效作用分析 為了評價提取效率，本文將超聲-微波協(xié)同萃取與單獨超聲波輔助、單獨微波輔助和傳統(tǒng)水浴回流提取法［7］進行比較，具體對比參數(shù)設置如表3所示。

1.3 數(shù)據(jù)處理

使用Excel和Design-Expert 8.0軟件對數(shù)據(jù)進行處理和分析，每組實驗均重復3次，實驗數(shù)據(jù)表示為平均值。

2 結果與分析

2.1 料液比對菠蘿蜜果皮多酚得率的影響

圖1 料液比對菠蘿蜜果皮多酚得率的影響Fig.1 The effect of solid-liquid ratio on polyphenols yield of jackfruit peel

在固液提取中，一般提取的活性物質會隨著提取溶劑的增多而增多，但當擴散達到平衡后，活性物質的量不再隨著溶劑的量增加而增加，而且，提取溶劑的增加還會對后續(xù)處理造成困難，因此需要選擇合適的料液比來提取目標物。由圖1可看出，開始時菠蘿蜜果皮中的多酚類化合物得率隨溶劑量的增大而增大，當料液比超過1∶30后，出現(xiàn)了隨著料液比的增加多酚提取量反而減小的現(xiàn)象。因此選擇的最佳料液比為1∶30。

2.2 乙醇體積分數(shù)對菠蘿蜜果皮多酚得率的影響

圖2 乙醇體積分數(shù)對菠蘿蜜果皮多酚得率的影響Fig.2 The effect of ethanol volume on polyphenols yield of jackfruit peel

溶劑提取主要是利用被提取物與提取溶劑的相似相溶原理，當提取溶劑的極性與菠蘿蜜果皮中多酚類化合物的極性相似時，能夠最有效的將多酚類化合物提取出來［16］。由圖2可知，隨著乙醇體積分數(shù)的增大，多酚得率呈先增大后減小的趨勢，當乙醇體積分數(shù)為60%時，多酚得率達到最大值。超過60%后，多酚得率開始逐漸降低，這是由于溶液極性增強，溶劑與多酚極性的相似程度降低，導致多酚類物質的溶解度下降［17-18］。所以60%為最優(yōu)乙醇體積分數(shù)。

2.3 微波功率對菠蘿蜜果皮多酚得率的影響

圖3 微波功率對菠蘿蜜果皮多酚得率的影響Fig.3 The effect of mcrowave power on polyphenols yield of jackfruit peel

表3 超聲-微波協(xié)同萃取與其他工藝提取效率的比較Table3 Extraction efficency comparison among different processes

由圖3可知，菠蘿蜜果皮多酚得率剛開始隨微波功率的增大而增加，這是由于微波的破壁效應、熱效應，超聲波的機械振動作用和空化作用，以及超聲波和微波的協(xié)同作用，使得菠蘿蜜果皮細胞內的多酚可以更好地釋放、進入溶劑并充分混合［19］。但是微波功率超過50 W，得率反而下降，可能是由于超聲-微波的熱效應產(chǎn)生的高溫使提取物的熱不穩(wěn)定性成分發(fā)生了分解的緣故。因此選用微波功率50 W比較適宜。

2.4 提取時間對菠蘿蜜果皮多酚得率的影響

圖4 超聲微波時間對菠蘿蜜果皮多酚得率的影響Fig.4 The effect of extraction time on polyphenols yield of jackfruit peel

由圖4可知，當微波超聲時間達到9 min時，多酚得率達到最大值；之后隨時間的延長，多酚得率反而減少。說明微波超聲9 min時多酚已提取完全，而隨著多酚在空氣中暴露時間的增加，有部分分解，使其得率有所降低。因此，最佳微波超聲時間為9 min。

2.5 響應面法優(yōu)化菠蘿蜜果皮多酚提取工藝

2.5.1 響應面實驗設計方案及結果 根據(jù)單因素實驗結果，采用時間（X1）、功率（X2）、乙醇濃度（X3）、料液比（X4）4因素3水平設計Box-Behnken實驗中心組合實驗，實驗方案及結果見表4。

2.5.2 數(shù)學模型的建立和方差分析 利用Design-Expert 8.0統(tǒng)計軟件對表4數(shù)據(jù)進行多元回歸擬合，得到多酚得率（Y）的二次多項回歸模型方程：

此模型的方差分析結果見表5。由表5可知，模型p值小于0.0001，說明二次回歸模型是極顯著的。失擬項p值為0.077（＞0.05），差異不顯著，未知因素對實驗結果干擾小，說明殘差均由隨機誤差引起。相關系數(shù)R2為0.9944，說明二次回歸模型能夠較好的反應響應值的變化，模型與實際實驗擬合程度好。

從各個因素的顯著性水平差異可知，對多酚得率的影響次序為：料液比＞功率＞乙醇濃度＞提取時間。其中，微波功率、乙醇濃度、液料比的一次項，時間與微波功率的交互項、時間與乙醇濃度的交互項、時間與料液比的交互項、微波功率與料液比的交互項對多酚得率的影響都達到了極顯著水平（p＜0.01）；微波功率與乙醇濃度的交互項對多酚得率的影響達到了顯著水平（0.01＜p＜0.05）。其因素間交互作用的響應面如圖5所示，與方差分析結果一致。

表4 Box-Behnken設計方案及實驗結果Table4 Experimental design and result of Box-Behnken

2.5.3 超聲-微波協(xié)同工藝的優(yōu)化組合及其驗證 根據(jù)以上的實驗結果分析以及模型擬合，利用Desig-Expert 8.0軟件分析并修正后得出菠蘿蜜果皮多酚提取的最佳工藝為：料液比1∶40，乙醇濃度70%，微波功率75 W，微波超聲時間12 min，菠蘿蜜果皮多酚得率的預期值約為7.23 mg/g。根據(jù)最佳工藝條件，進行三次平行實驗所得的菠蘿蜜果皮多酚得率平均值為7.19 mg/g，與預期值相接近，說明該響應面法得到的回歸模型具有一定的可靠性。

2.6 抗氧化活性測定

菠蘿蜜提取物的抗氧化活性測定結果如表6所示，從表6可以看出，利用超聲-微波協(xié)同萃取的菠蘿蜜果皮多酚提取物對DPPH自由基和ABTS自由基均有較強的清除能力，且隨濃度的增加而增強，呈量效關系，其EC50值分別為101.39 μg/mL和106.60 μg/mL，抑制效果略低于陽性對照物BHT（89.29 μg/mL）。并且，多酚提取物濃度與DPPH自由基清除率和ABTS自由基的清除率之間具有較好的線性關系，R2分別為0.9657和0.9913，表明多酚是菠蘿蜜果皮中具有抗氧化活性的物質基礎。

表5 回歸模型的方差分析Table5 Analysis of variance of the regression model

表6 菠蘿蜜果皮多酚提取物抗氧化活性實驗結果Table6 The antioxidant activity results of the polyphenol extracts from jackfruit peel

2.7 超聲-微波協(xié)同萃取高效作用分析

為了評價超聲波-微波協(xié)同提取效率，本文將其與單獨超聲波輔助、單獨微波輔助和傳統(tǒng)水浴回流提取法進行比較，測定結果如圖6所示。由圖6a可知，協(xié)同最佳工藝所得多酚得率顯著高于其他三種提取工藝（p＜0.05），超聲、微波、水浴回流的多酚得率分別為協(xié)同提取物的86.87%、83.44%和14.38%，其中傳統(tǒng)水浴回流提取得率最低，為1.04 mg/g。同時由圖6 b可見，協(xié)同提取最佳工藝所得提取物DPPH、ATBS值也顯著高于其他三種提取工藝所得提取物（p＜0.05）。對比結果說明，超聲-微波協(xié)同萃取方法提取效率最高，展現(xiàn)出更高的多酚得率和抗氧化活性，原因是超聲-微波協(xié)同萃取法結合了超聲振動的空化作用與微波的高能作用，達到優(yōu)勢互補，因而對植物的有效成分能快速高效的提取，且不破壞所萃取有機分子的的結構［20］。陳衛(wèi)云等［21］從荔枝果肉中提取多糖時，通過比較超聲-微波協(xié)同法與傳統(tǒng)熱水法、超聲法、微波法，也得出了相似的結論。

圖5 因素交互作用對總酚含量影響的響應面Fig.5 Response surface for interactive effects four process parameters on polyphenols yield

圖6 超聲-微波協(xié)同萃取與其他工藝提取效率的比較Fig.6 Extraction efficency comparison among different processes

3 結論

采用超聲-微波協(xié)同技術（UMAE），通過響應面實驗設計對菠蘿蜜果皮中多酚的提取工藝進行了優(yōu)化，其最佳工藝條件為料液比1∶40，乙醇體積分數(shù)70%，微波功率75 W，提取時間12 min。在該條件下，多酚得率為7.19 mg/g，其提取效率優(yōu)于傳統(tǒng)水浴回流（1.04 mg/g）、微波輔助法（5.23 mg/g）和超聲輔助法（5.89 mg/g）?？寡趸钚匝芯勘砻?，菠蘿蜜果皮多酚提取物對DPPH自由基和ABTS自由基均有較強的清除能力，呈量效關系，其EC50值分別為101.39 μg/mL 和106.60 μg/mL，表明多酚是菠蘿蜜果皮抗氧化活性的物質基礎。因此，該方法的應用對進一步開發(fā)菠蘿蜜下腳料的應用價值具有重大的借鑒意義。

［1］Kategunya R，Sanguansri C.Physico-chemical properties of jackfruit seed starch（Artocarpus heterophyllus）and its application as a thickener and stabilizer in chilli sauce［J］.LWT-Food Science and Technology，2011（4）：1309-1313.

［2］姚定，董明.菠蘿蜜果皮中果膠的提取工藝優(yōu)化［J］.食品科學，2009，30（20）：165-168.

［3］臺建祥，范鴻雁，薛慧，等.菠蘿蜜果皮果膠提取工藝的優(yōu)化［J］.食品科技，2014，39（1）：236-239.

［4］Chen AY，Chen YC.A review of the dietary flavonoid，kaempferol on human health and cancer chemoprevention［J］.Food Chemistry，2013，138（4）：2099-2107.

［5］陳亮，孫鵬，王婷，等.天竺桂多酚提取物降血糖活性研究［J］.中成藥，2014，36（2）：229-234.

［6］Agourram A，Ghirardello D，Rantsiou K，et al.Phenolic content，antioxidant potential，and antimicrobial activities of fruit and vegetable by-product extracts［J］.International Journal of Food Properties，2013，16（5）：1092-1104.

［7］王佳珺，文俊，何鋼.菠蘿蜜果皮中多酚的測定及其應用初探［J］.廣州化工，2014，42（11）：124-125.

［8］Cristina L，Timothy JM.Microwave and ultrasonic processing：now a realistic option for industry［J］.Chemical Engineering and Processing：Process Iintensification，2010，49（9）：885-900.

［9］Chen Y，Gu X，Huang SQ，et al.Optimization of ultrasonic/ microwave assisted extraction（UMAE）of polysaccharides from inonotus obliquus and evaluation of its anti-tumor activities［J］.International Journal of Biological Macromolecules，2010，46（4）：429-435.

［10］Zhang LF，Liu ZL.Optimization and comparison of ultrasound/ microwave assisted extraction（UMAE）and ultrasonic assisted extraction（UAE） of lycopene from tomatoes［J］.Ultrasonics Sonochemistry，2008，15（5）：731-737.

［11］Azmira J，Zaidul ISM，Rahmana MM，et al.Techniques for extraction of bioactive compounds from plant materials：A review ［J］.Journal of Food Engineering，2013，117（4）：426-436.

［12］Liu RL，Song SH，Wu M，et al.Rapid analysis of fatty acid profiles in raw nuts and seeds by microwave-ultrasonic synergistic in situ extraction-derivatisation and gas chromatography-mass spectronmetry［J］.Food Chemistry，2013，141（4）：4269-4277.

［13］Fang XS，Wang JH，Wang YZ，et al.Optimization of ultrasonicassisted extraction of wedelolactone and antioxidant polyphenols from Eclipta prostrate L.using response surface methodology［J］.Separation and Purification Technology，2014，138：55-64.

［14］Liu Y，Wei SL，Liao MC.Optimization of ultrasonic extraction of phenolic compounds from Euryale ferox seed shells using response surface methodology［J］.Industrial Crops and Products，2013，49：837-843.

［15］Gan CY，Latiff AA.Optimisation of the solvent extraction of bioactive compounds from Parkia speciosa pod using response surface methodology［J］.Food chemistry，2011，124（3）：1277-1283.

［16］柯增光，許偉，邵榮，等.蓖麻堿的超聲波輔助提取及其抗氧化性能［J］.南京工業(yè)大學學報：自然科學版，2015，37（1）：14-19.

［17］趙曉云，趙博，邢媛媛，等.荷葉多酚的微波輔助提取工藝優(yōu)化及其抗氧化活性研究［J］.中國釀造，2010，216（3）：79-83.

［18］陳義勇，竇祥龍，黃友如，等.響應面法優(yōu)化超聲-微波協(xié)同輔助提取茶多糖工藝［J］.食品科學，2012，33（4）：100-103.

［19］江慎華，趙錦，劉常金，等.采用超聲-微波協(xié)同技術提取丁香總多酚工藝優(yōu)化及其高效原因分析［J］.食品工業(yè)科技，2014，35（24）：243-247.

［20］蘇堯堯，童群義.超聲微波協(xié)同萃取垂柳葉中α-葡萄糖苷酶抑制劑的研究［J］.食品工業(yè)科技，2014，35（2）：235-241.

［21］陳衛(wèi)云，張名位，廖森泰，等.荔枝多糖超聲微波酶解協(xié)同提取工藝優(yōu)化［J］.中國食品學報，2013，13（5）：77-84.

Optimization of ultrasonic-microwave assisted extraction of polyphenols from jackfruit peel and evaluation of its antioxidant activities

JIANG Zhi-guo1，LI Bin1，*，WANG Yan-hua2，*，ZHANG Hai-de1
（1.College of Food Science，Hainan University，Haikou 570228，China；2.College of Information，Hainan University，Haikou 570228，China）

In order to maximize the yield of polyphenols from jackfruit peel，response surface methodology（RSM）was employed to optimize the ultrasonic-microwave assisted extraction（UMAE）conditions，and the antioxidant activities of the polyphenol extracts were evaluated in vitro.On the basis of single-factor experiments，the Box-Behnken center-composite experimental design principle was used to design the experiments and the response surface analysis with 4 factors and 3 levels was adopted.A second order quadratic equation for extraction of polyphenols was created.After analyzed the significance of various factors and their interactions，the optimal extraction conditions for polyphenols from jackfruit peel were determined as follows:70%aqueous ethanol as the extraction solvent，1∶40 as solid-liquid ratio，extraction time 12 min and microwave power 75 W.The extraction yield of polyphenols was determined to be 7.19 mg/g under the optimized extraction conditions，which were above that of traditional hot water extraction，microwave assisted extraction and ultrasonic assisted extraction（1.04 mg/g，5.23mg/g and 5.89 mg/g，respectively）.These results confirmed that UMAE of polyphenols had great potential and efficiency compared with traditonal extraction technologies.Furthermore，the experiments of antioxidant activity showed that polyphenol extracts of jackfruit peel exhibited strong antioxidant activities by the DPPH and ABTS radical scavenging assays，and the EC50were 101.39 μg/mL and 106.60 μg/mL，respectively.The present study showed that the jackfruit peel were strong radical scavengers and could be considered as a good source of natural antioxidants for medicinal and commercial uses.

jackfruit peel；polyphenols；microwave-ultrasonic assisted extraction；response surface methodology；antioxidant activity

TS209

B

1002-0306（2016）02-0270-06

10.13386/j.issn1002-0306.2016.02.046

2015-07-02

蔣志國（1977-），男，博士，副教授，研究方向：食品化學與營養(yǎng)，E-mail：jzg6666@163.com。

*通訊作者：李斌（1966-），男，大學本科，實驗師，研究方向：食品化學與營養(yǎng)，E-mail：lzc66@21cn.com。

王燕華（1980-），女，大學本科，講師，研究方向：統(tǒng)計分析，E-mail：wangyanhua0007@163.com。

海南省自然科學基金（20153157）；海南省高等學校科學研究重點項目（Hnky2015ZD-3）。