基于多項(xiàng)式回歸分析的摩爾多瓦葡萄富硒效應(yīng)

李偉 蔣代華 吳代東 白厚義 龐東 吳艷艷 黃偉華 黃永才

摘要：【目的】探討基于多項(xiàng)式回歸分析的摩爾多瓦葡萄富硒效應(yīng)，為開展廣西富硒土原位生產(chǎn)富硒葡萄提供參考依據(jù)。【方法】以摩爾多瓦葡萄為試驗(yàn)材料，采用三元二次飽和D-最優(yōu)設(shè)計(jì)，以石灰石粉（x1）、AM菌劑（x2）、營(yíng)養(yǎng)調(diào)節(jié)劑（x3）為因變量，以試驗(yàn)結(jié)果中葡萄硒含量為自變量（y），設(shè)11個(gè)處理（處理1～處理11），經(jīng)回歸分析明確試驗(yàn)因素對(duì)葡萄硒含量的影響，優(yōu)選生產(chǎn)富硒摩爾多瓦葡萄的外源物質(zhì)施用量。【結(jié)果】施用石灰石粉、AM菌劑和營(yíng)養(yǎng)調(diào)節(jié)劑各處理（除處理1外）摩爾多瓦葡萄的硒含量比對(duì)照（CK）提高1.1～1.8倍，最高增幅達(dá)87.2%。多項(xiàng)式回歸方程中單因素效應(yīng)分析結(jié)果顯示，當(dāng)石灰石粉、AM菌劑和營(yíng)養(yǎng)調(diào)節(jié)劑施用量分別為2400.00、69.00和39.98 kg/ha時(shí)，摩爾多瓦葡萄硒含量最高，為11.56 μg/kg；互作效應(yīng)分析結(jié)果顯示，兩兩互作對(duì)摩爾多瓦葡萄硒含量的影響排序?yàn)閤1x3>x1x2>x2x3。【結(jié)論】摩爾多瓦葡萄分別開溝施用石灰石粉、AM菌劑和營(yíng)養(yǎng)調(diào)節(jié)劑2400.00、69.00和39.98 kg/ha，其果實(shí)硒含量達(dá)富硒水平。因此，施用石灰石粉、AM菌劑和營(yíng)養(yǎng)調(diào)節(jié)劑等外源物質(zhì)可作為生產(chǎn)富硒葡萄的新型技術(shù)手段。

關(guān)鍵詞： 摩爾多瓦葡萄；富硒土壤；富硒葡萄；二次飽和D-最優(yōu)設(shè)計(jì)；營(yíng)養(yǎng)調(diào)節(jié)劑

中圖分類號(hào)： S663.1 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：2095-1191（2017）05-0849-05

Selenium enriched Moldova grape research based on

polynomial regression analysis

LI Wei1，2， JIANG Dai-hua1 *， WU Dai-dong 2 *， BAI Hou-yi1， PANG Dong 1， WU Yan-yan2， HUANG Wei-hua2， HUANG Yong-cai2， LIU Jie-yun2， MU Hai-fei2，

WEI Shao-long 2， HUANG Yan-fei3， LIU Yong-xian3

（1 College of Agriculture， Guangxi University， Nanning 530004， China； 2 Biotechnology Research Institute， Guangxi Academy of Agricultural Sciences， Nanning 530007， China； 3 Institute of Agricultural Resources and Environment，

Guangxi Academy of Agricultural Sciences， Nanning 530007， China）

Abstract：【Objective】Selenium enriched Moldova grape was researched based on polynomial regression analysis， so as to provide reference for planting selenium enriched grapes in Guangxi selenium enriched field without adding selenium substance. 【Method】Moldova grape was used as the material in experiment. With ternary quadratic saturation D-optimum design，limestone powder（x1），AM fungal inoculum（x2）， nutrient regulators（x3） as dependent variables， selenium content in grape as independent variable（y）， 11 treatments（treatment 1-11） were set to conduct regression analysis. Through the experiment， effects of experiment factors on selenium content in Moldova grape were studied， in order to determine application rates of allogenic materials for selenium enriched grapes. 【Result】In all treatments applying limestone powder， AM fungal inoculum and nutrient regulators（except treatment 1）， selenium contents in Moldova grape increased 1.1 to 1.8 times compared with CK， the maximum increasement reached 87.2%. Single factor effect results showed that selenium content in grapes was the highest when applying 2400.00 kg/ha limestone powder， 69.00 kg/ha AM fungal inoculum and 39.98 kg/ha nutrient regulators， and the selenium content in Moldova grape was 11.56 μg/kg. Interaction effects analysis indicated that effects rank of two factors interaction on Moldova grape were as follows： x1x3>x1x2>x2x3. 【Conclusion】Applying 2400.00 kg/ha limestone powder， 69.00 kg/ha AM fungal inoculum and 39.98 kg/ha nutrient regulators to Moldova grape plants， the fruits can reach selenium enriched level. Therefore， application of limestone powder， AM fungal inoculum and nutrient regulators can be utilized as a new technique for planting selenium enriched grapes.

Key words： Moldova grape； selenium enriched soil； selenium enriched grape； quadratic saturation D-optimum design； nutrient regulator

0 引言

【研究意義】廣西以亞熱帶氣候?yàn)橹鳎瑢俑邷囟嘤甑貐^(qū)，也是我國(guó)葡萄特殊優(yōu)勢(shì)種植區(qū)，近年來(lái)葡萄產(chǎn)量年均增長(zhǎng)10%以上（白先進(jìn)等，2010），栽培面積已超過(guò)2.67萬(wàn)ha（孫欣等，2015），以鮮食葡萄居多，品種以巨峰和夏黑為主，但這兩個(gè)品種的抗病性、成花率和坐果率在廣西地區(qū)表現(xiàn)欠佳。廣西農(nóng)業(yè)科學(xué)院通過(guò)引種赤霞珠、西拉、黑玫瑰、陽(yáng)光玫瑰、美人指和摩爾多瓦等品種，對(duì)該地區(qū)葡萄品種結(jié)構(gòu)進(jìn)行調(diào)整并取得了一系列成果（周詠梅等，2015）。摩爾多瓦葡萄原產(chǎn)摩爾多瓦國(guó)，我國(guó)于1997年從羅馬尼亞引入（張潰珍和劉煥，2015），已在廣西、云南、貴州、福建、河南和寧夏等地推廣種植，表現(xiàn)高抗葡萄霜霉病，高萌芽率，成熟晚，耐貯運(yùn)，極豐產(chǎn)，成活率較高，有良好的發(fā)展前景（彭羽等，2010；何志山，2013）。硒是人體必需元素，具有抗氧化、抗癌、拮抗重金屬等作用。廣西富硒土壤面積達(dá)212多萬(wàn)ha，是我國(guó)多目標(biāo)區(qū)域地球化學(xué)調(diào)查工作圈定的最大面積連片富硒土壤區(qū)域，已開發(fā)出富硒水稻、富硒玉米和富硒金桔果糕等（何雪梅等，2015；劉永賢等，2015），但針對(duì)葡萄富硒方面的研究較少。因此，開展摩爾多瓦葡萄富硒效應(yīng)研究，既可為廣西富硒葡萄研究提供參考依據(jù)，也可為廣西地區(qū)硒資源開發(fā)和樹立地方富硒品牌打下基礎(chǔ)。【前人研究進(jìn)展】人類對(duì)自然界中無(wú)機(jī)硒的吸收主要通過(guò)以動(dòng)物、植物和微生物為載體進(jìn)行富硒食品研究開發(fā)而獲得，植物富硒是通過(guò)以植物為載體、施用硒肥來(lái)實(shí)現(xiàn)無(wú)機(jī)硒有機(jī)化進(jìn)而生產(chǎn)富硒食品（陳彥霖和范艷麗，2008）。朱麗琴等（2007）研究表明，在不施用硒肥前提下，葡萄對(duì)硒的累積順序?yàn)楦?葉>果實(shí)，施硒條件下其累積規(guī)律為葉>根>果實(shí)。王海波等（2011）研究表明，葉面噴施硒肥后，葡萄對(duì)硒累積的順序?yàn)槿~片>果肉>葉柄，而作土施處理，其累積順序?yàn)楣?葉片>葉柄。陳劍俠（2013）研究發(fā)現(xiàn)，葉面噴硒后，硒在葡萄器官的分布表現(xiàn)為葉>果>根。酸性環(huán)境中土壤硒主要以四價(jià)存在，中性或堿性條件下以六價(jià)硒存在，由于六價(jià)態(tài)硒性質(zhì)較穩(wěn)定，通常將其視為植物主要吸收價(jià)態(tài)（張華華和康玉凡，2013）。王忠（2015）在溫克葡萄盛花期和初果期噴施瓜果型鋅硒葆葉面肥，生產(chǎn)的葡萄達(dá)富硒標(biāo)準(zhǔn)。【本研究切入點(diǎn)】廣西地區(qū)土壤呈酸性或弱酸性，通過(guò)施用石灰、AM菌劑和營(yíng)養(yǎng)調(diào)節(jié)劑調(diào)節(jié)其pH等理化性質(zhì)以實(shí)現(xiàn)硒的四、六價(jià)間價(jià)態(tài)轉(zhuǎn)換進(jìn)而提高葡萄對(duì)硒吸收效率的研究鮮見報(bào)道。【擬解決的關(guān)鍵問(wèn)題】以摩爾多瓦葡萄為試驗(yàn)材料，以石灰石粉、AM菌劑和營(yíng)養(yǎng)調(diào)節(jié)劑為試驗(yàn)因素，采用三元二次飽和D-最優(yōu)設(shè)計(jì)，檢驗(yàn)廣西富硒土原位生產(chǎn)富硒葡萄的可能性，為富硒葡萄生產(chǎn)提供新的研究手段。

1 材料與方法

1. 1 試驗(yàn)地概況

試驗(yàn)點(diǎn)位于廣西南寧市武鳴區(qū)雙橋鎮(zhèn)楊李新村合作葡萄試驗(yàn)示范基地，土壤為赤紅土，pH 4.5，有機(jī)質(zhì)含量33.13 g/kg，堿解氮含量400.0 mg/kg，速效磷（P2O5）含量13.2 mg/kg，速效鉀（K2O）含量110.5 mg/kg，全硒含量1.0 mg/kg。

1. 2 試驗(yàn)材料

供試摩爾多瓦葡萄樹為一年生移栽樹；石灰石粉（方解石粉末）購(gòu)自武鳴區(qū)石粉加工廠，其CaCO3含量在95%以上；AM（叢枝菌根真菌）菌劑由廣西農(nóng)業(yè)科學(xué)院微生物研究所提供，孢子量為150個(gè)/g；營(yíng)養(yǎng)調(diào)節(jié)劑（主要成分為胡敏酸鉀）制作方法由廣西大學(xué)白厚義教授提供。

1. 3 試驗(yàn)方法

1. 3. 1 試驗(yàn)設(shè)計(jì) 以石灰石粉（x1）、AM菌劑（x2）、營(yíng)養(yǎng)調(diào)節(jié)劑（x3）為試驗(yàn)因素，參照白厚義（2003）的方法，采用三元二次飽和D-最優(yōu)設(shè)計(jì)進(jìn)行處理方案設(shè)計(jì)。試驗(yàn)因素的上、下水平及變化間距見表1，試驗(yàn)處理方案見表2。田間試驗(yàn)采用隨機(jī)排列，每3株葡萄為一個(gè)小區(qū)，共11個(gè)處理（處理1～處理11），3次重復(fù)。于葡萄萌芽期，在葡萄樹種植穴兩邊開溝，將各處理所施物料與5.0 kg土混合均勻后施于溝中填土掩埋，按常規(guī)田間管理方法進(jìn)行管理。以按常規(guī)方法施肥（底肥）為對(duì)照（CK）。

1. 3. 2 測(cè)定方法 土壤pH采用pH計(jì)測(cè)定（水土比為1∶1），有機(jī)質(zhì)含量采用重鉻酸鉀外加熱法測(cè)定，速效磷采用NH4F-HCl比色法測(cè)定，速效鉀采用火焰光度法測(cè)定，全硒含量參考瞿建國(guó)等（1998）的方法進(jìn)行測(cè)定；葡萄全硒含量參考國(guó)家標(biāo)準(zhǔn)GB/T 5009.93-2003進(jìn)行測(cè)定。

1. 4 統(tǒng)計(jì)分析

試驗(yàn)數(shù)據(jù)采用Excel 2007進(jìn)行統(tǒng)計(jì)，采用SAS 9.2進(jìn)行多項(xiàng)式回歸分析和差異顯著性檢驗(yàn)，采用Qbasic進(jìn)行摩爾多瓦葡萄硒含量尋優(yōu)。

2 結(jié)果與分析

2. 1 多項(xiàng)式回歸方程的構(gòu)建及差異顯著性檢測(cè)結(jié)果

2. 1. 1 多項(xiàng)式回歸方程的建立 表3為不同水平石灰石粉、AM菌劑和營(yíng)養(yǎng)調(diào)節(jié)劑處理下的葡萄硒含量，以x1、x2和x3 3個(gè)因素編碼值為控制變量，葡萄硒含量為依變量（y），經(jīng)回歸分析得到三元二次多項(xiàng)式回歸方程：

y =11.55352 - 0.17479x1 + 0.76398x2 + 0.07329x3 +

0.85746x1x2-0.85134x1x3-0.04591x2x3-0.72185x12-

1.77310x22-1.97152x32 ①

2. 1. 2 多項(xiàng)式回歸方程的差異顯著性檢測(cè)結(jié)果 經(jīng)SAS 9.2運(yùn)算得R=0.770>R0.01=0.514，F(xiàn)=3.61>F0.01=3.26，方程達(dá)極顯著差異水平，說(shuō)明進(jìn)行多項(xiàng)式回歸分析得到的三元二次多項(xiàng)式回歸方程擬合度高，能反映各試驗(yàn)因素與葡萄硒含量間的效應(yīng)關(guān)系，石灰石粉、AM菌劑、營(yíng)養(yǎng)調(diào)節(jié)劑與葡萄含硒量相關(guān)極顯著（P<0.01）。

從表3還可看出，除處理1外，其余處理葡萄的硒含量分別比CK提高1.1～1.8倍。其中，處理5即石灰石粉、營(yíng)養(yǎng)調(diào)節(jié)劑和AM菌劑分別為300.00、68.66和46.72 kg/ha的硒含量最高，為11.29 μg/kg，比CK高87.2%，而處理1即石灰石粉、營(yíng)養(yǎng)調(diào)節(jié)劑和AM菌劑分別為300.00、15.00和4.95 kg/ha時(shí)，葡萄中硒含量最低，為5.51 μg/kg，比CK低8.06%，說(shuō)明石灰石粉、營(yíng)養(yǎng)調(diào)節(jié)劑、AM菌劑不同配比方案對(duì)土壤硒生物有效性影響較大。

2. 2 多項(xiàng)式回歸方程的單因素效應(yīng)分析結(jié)果

將多項(xiàng)式回歸方程進(jìn)行降維可得到以下3個(gè)單因素效應(yīng)回歸方程yi（i=1，2，3）（1代表石灰石粉，2代表AM菌劑，3代表營(yíng)養(yǎng)調(diào)節(jié)劑）：

y1=11.55352-0.17479x1-0.72185x12②

y2=11.55352+0.76398x2-1.77310x22③

y3=11.55352+0.07329x3-1.97152x32④

分別令yi為0，然后分別對(duì)x求導(dǎo)可得：x1=-0.121，x2=0.215，x3=0.018。

上述單因素效應(yīng)回歸方程中，一次項(xiàng)系數(shù)絕對(duì)值0.76398>0.17479>0.07329，即各因素作用的排序依次為x2>x1>x3，將表3中相應(yīng)編碼值代入單因素回歸方程式中計(jì)算，產(chǎn)生各因素對(duì)產(chǎn)量的單效應(yīng)圖（圖1）。由圖1可看出，x1、x2和x3在編碼水平為[-1，0]時(shí)，葡萄硒含量隨著各因素施用量的增加而增高；在編碼水平為[0，1]時(shí)，葡萄含硒量隨各因素施入量增加而逐漸降低。在石灰石粉的效應(yīng)中，編碼水平為-1、[-1，0]和1時(shí)，對(duì)應(yīng)的pH分別為6.3～6.5、6.5～7.0和7.0～7.4。說(shuō)明摩爾多瓦葡萄對(duì)硒吸收的土壤pH適宜范圍為6.5～7.0。

x1、x2和x3在編碼水平分別為-0.121、0.215和0.018時(shí)葡萄硒含量最高，分別為11.56、11.64和11.55 μg/kg，此3組硒含量均達(dá)富硒標(biāo)準(zhǔn)，且編碼值均趨于0，說(shuō)明試驗(yàn)設(shè)計(jì)梯度合理有效。

2. 3 多項(xiàng)式回歸方程的兩因素交互效應(yīng)分析結(jié)果

對(duì)三元二次多項(xiàng)式回歸方程采用降維法得到二元二次多項(xiàng)式回歸方程：

y12=11.55352-0.17479x1+0.76398x2+0.85746x1x2-

0.72185x12-1.77310x22 ⑤

y13=11.55352-0.17479x1+0.07329x3-0.85134x1x3-

0.72185x12-1.97152x32 ⑥

y23=11.55352+0.76398x2+0.07329x3-0.04591x2x3-

1.77310x22-1.97152x32 ⑦

2. 3. 1 x1與x2的交互效應(yīng)分析 將各編碼（-1、-0.2912、0、0.1925、1）代入方程式⑤得到x1與x2的交互效應(yīng)（表4）。由表4可看出，當(dāng)x1與x2處于低量（-1、-1）時(shí)葡萄硒含量最低（9.32 μg/kg），隨著x1和x2用量增加，葡萄硒含量也逐漸增加，當(dāng)x1和x2用量增加至（0.1925、0.1925）時(shí)，葡萄硒含量最高（11.61 μg/kg）；當(dāng)x1和x2用量繼續(xù)增加時(shí)，硒含量反而下降；x1和x2的組合為（0、0.1925）時(shí)，硒含量達(dá)最高（11.63 μg/kg）。

2. 3. 2 x1與x3的交互效應(yīng)分析 將各編碼（-1、-0.2912、0、0.1925、1）代入方程式⑥得到x1與x3的交互效應(yīng)（表5）。由表5可看出，當(dāng)x1和x3處于低量（-1、-1）時(shí)，葡萄硒含量最低（8.11 μg/kg），隨著x1和x3用量增加，葡萄硒含量也逐漸增加，當(dāng)增加至（0、0）時(shí)，葡萄的硒含量達(dá)最高（11.55 μg/kg）；當(dāng)x1和x3用量繼續(xù)增加時(shí)，硒含量反而下降，且下降至最低含量（9.96 μg/kg）；當(dāng)x1和x3組合為（0、0）時(shí)，摩爾多瓦葡萄的硒含量最高（11.55 μg/kg），石灰石粉與營(yíng)養(yǎng)調(diào)節(jié)劑的互作效應(yīng)最大。

2. 3. 3 x2與x3的交互效應(yīng)分析 將各編碼（-1、-0.2912、0、0.1925、1）代入方程式⑦得到x2與x3的交互效應(yīng)（表6）。由表6可看出，當(dāng)x2和x3處于低量（-1、-1）時(shí)，葡萄硒含量最低（6.93 μg/kg）；隨著x2和x3用量增加，葡萄硒含量也逐漸增加，當(dāng)增加至（0.1925、0.1925）時(shí)，葡萄硒含量較高（11.57 μg/kg）；當(dāng)x2和x3用量繼續(xù)增加時(shí)，硒含量反而下降（降至7.16 μg/kg）；而x2和x3組合為（0、0.1925）時(shí)，硒含量最高（11.63 μg/kg）。說(shuō)明AM菌劑與營(yíng)養(yǎng)調(diào)節(jié)劑分別為60.00和46.72 kg/ha時(shí)，二者的互作效應(yīng)最大。

借助Qbasic進(jìn)行尋優(yōu)，獲得摩爾多瓦葡萄最高硒含量即y為11.56 μg/kg時(shí)石灰石粉、AM菌劑和營(yíng)養(yǎng)調(diào)節(jié)劑3個(gè)試驗(yàn)因素對(duì)應(yīng)的編碼為0（2400.00 kg/ha）、0.2（69.00 kg/ha）和0（39.98 kg/ha）。

3 討論

本研究尋優(yōu)結(jié)果顯示，摩爾多瓦葡萄最高硒含量為11.56 μg/kg，達(dá)到陳劍俠等（2013）推薦的富硒葡萄硒含量（10.00～100.00 μg/kg）水平，且比陳劍俠（2013）對(duì)夏黑葡萄噴施加康硒富硒葉面肥的處理結(jié)果（硒含量為11.20 μg/kg）略高，表明通過(guò)改變富硒土壤理化性質(zhì)能提高葡萄對(duì)硒的吸收效率，克服富硒土壤補(bǔ)充外源硒導(dǎo)致硒過(guò)量的風(fēng)險(xiǎn)，也說(shuō)明直接以施用外源物質(zhì)改變富硒土壤理化性質(zhì)作為新型富硒手段具有可行性。本研究中摩爾多瓦葡萄硒含量提高了1.1～1.8倍，低于王忠（2015）對(duì)溫克葡萄葉面噴施瓜果型鋅硒葆的硒含量，因此，今后仍需對(duì)富硒土壤開展原位富硒葡萄有關(guān)作用機(jī)理、外源物質(zhì)優(yōu)化組配等方面的探究。

pH是影響硒生物有效性的因素之一，不同pH對(duì)其影響也各不相同，因?yàn)槲诓煌嵝詶l件下所表現(xiàn)的價(jià)態(tài)不同（酸性條件下硒主要以四價(jià)態(tài)形式存在，堿性條件下主要以六價(jià)態(tài)存在），植物對(duì)不同價(jià)態(tài)硒的吸收能力也不同。廖金鳳（1999）研究表明，pH為2～3時(shí)，粘粒對(duì)硒的吸附量最大，而后隨著pH的升高而減少，當(dāng)pH為6～8時(shí)，基本趨于穩(wěn)定。在本研究對(duì)石灰石粉單因素效應(yīng)的分析中，摩爾多瓦硒含量總體表現(xiàn)為低—高—低的凸型拋物線變化，該曲線各段對(duì)應(yīng)的pH分別為6.3～6.5、6.5～7.0和7.0～7.4。李嬌嬌等（2012）研究表明，pH為6時(shí)豆瓣菜的含硒量遠(yuǎn)高于pH為7時(shí)的含硒量，本研究結(jié)果與其不同，可能是由于物種類型差異所造成。

4 結(jié)論

摩爾多瓦葡萄分別開溝施用石灰石粉、AM菌劑和營(yíng)養(yǎng)調(diào)節(jié)劑2400.00、69.00和39.98 kg/ha，其果實(shí)硒含量達(dá)11.56 μg/kg，處于富硒水平。因此，施用石灰石粉、AM菌劑和營(yíng)養(yǎng)調(diào)節(jié)劑等外源物質(zhì)可作為生產(chǎn)富硒葡萄的新型技術(shù)手段。

參考文獻(xiàn)：

白厚義. 2003. 回歸設(shè)計(jì)及多元統(tǒng)計(jì)分析[M]. 南寧：廣西科學(xué)技術(shù)出版社. [Bai H Y. 2003. Regression Design and Multivariate Statistical Analysis[M]. Nanning：Guangxi Science and Technology Publishing House.]

白先進(jìn)，王舉兵，陳愛軍. 2010. 廣西葡萄產(chǎn)業(yè)發(fā)展的思考[J]. 廣西農(nóng)學(xué)報(bào)，25（1）：29-32. [Bai X J，Wang J B，Chen A J. 2010. Considerations on the grape industry development in Guangxi[J]. Journal of Guangxi Agriculture，25（1）：29-32.]

陳劍俠. 2013. 葡萄葉片和果實(shí)對(duì)硒的吸收和累積特性[J]. 江蘇農(nóng)業(yè)科學(xué)，41（7）：154-156. [Chen J X. 2013. Absorption and accumulation characteristics of selenium in grape leaves and fruits[J]. Jiangsu Agricultural Sciences，41（7）：154-156.]

陳劍俠，楊伯忠，范麗華，賴呈純，鄭京平. 2013. 富硒葡萄中硒限量標(biāo)準(zhǔn)的研究[J]. 福建分析測(cè)試，22（4）：31-33. [Chen J X，Yang B Z，F(xiàn)an L H，Lai C C，Zheng J P. 2013. A research on the selenium limit in the selenium-rich grape[J]. Fujian Analysis & Testing，22（4）：31-33.]

陳彥霖，范艷麗. 2008. 硒的生物學(xué)功能及富硒食品的研究進(jìn)展[J]. 河北農(nóng)業(yè)科學(xué)，12（9）：99-100. [Chen Y L，F(xiàn)an Y L. 2008. Research progress on the biological function of selenium and Se-enrich foodstuff[J]. Joumal of Hebei Agricultural Sciences，12（9）：99-100.]

何雪梅，孫健，李麗，盛金鳳，鄭鳳錦，李杰民，劉國(guó)明，李昌寶. 2015. 富硒金桔果糕生產(chǎn)工藝研究[J]. 西南農(nóng)業(yè)學(xué)報(bào)，28（6）：2736-2741. [He X M，Sun J，Li L，Sheng J F，Zheng F J，Li J M，Liu G M，Li C B. 2015. Study on production process of selenium-rich kumquat fruit cake[J]. Southwest China Journal of Agricultural Sciences，28（6）：2736-2741.]

何志山. 2013. 葡萄新品種摩爾多瓦在會(huì)澤縣的引種表現(xiàn)[J]. 現(xiàn)代農(nóng)業(yè)科技，（24）：119-123. [He Z S. 2013. Introduction of new grape variety in Moldova Huize county[J]. Modern Agricultural Science and Technology，（24）：119-123.]

李嬌嬌，原保忠，余展深，王明紅. 2012. 不同外界條件對(duì)豆瓣菜富硒能力的影響研究[J]. 長(zhǎng)江蔬菜，（16）：60-63. [Li J J，Yuan B Z，Yu Z S，Wang M H. 2012. Effects of different external conditions on selenium rich ability of Watercress （Nasturtium officinale R. Br.）[J]. Journal of Changjiang Vegetables，（16）：60-63.]

廖金鳳. 1999. 土壤粒級(jí)和酸度對(duì)硒吸附固定的影響[J]. 廣東微量元素科學(xué)，6（1）：47-49. [Liao J F. 1999. Effects of soil separate and soil acidity on selenium absorption and fixation[J]. Guangdond Trace Elements Science，6（1）：47-49.]

劉永賢，陽(yáng)繼輝，石玫莉，江澤普，邢穎，吳志林. 2015. 廣西富硒功能農(nóng)產(chǎn)品開發(fā)前景分析[J]. 農(nóng)業(yè)與技術(shù)，35（1）：176-178. [Liu Y X，Yang J H，Shi M L，Jiang Z P，Xing Y，Wu Z L. 2015. Prospects for selenium-enriched agriculture products development in Guangxi[J]. Agriculture & Technology，35（1）：176-178.]

彭羽，袁飛，周蕓蕓，王艷杰，馮金朝. 2010. 七個(gè)葡萄品種引種到銀川的植物學(xué)性狀表現(xiàn)初探[J]. 北方園藝，（9）： 44-46. [Peng Y，Yuan F，Zhou Y Y，Wang Y J，F(xiàn)eng J C. 2010. Primary study on plant characteristics of introduced 7 grapes varieties in Yinchuan[J]. Northern Horticulture，（9）：44-46.]

瞿建國(guó)，徐伯興，龔書椿. 1998. 氫化物發(fā)生—無(wú)色散原子熒光光度法測(cè)定土壤中有效態(tài)硒和總硒[J]. 土壤通報(bào)，29（1）：47. [Qu J G，Xu B X，Gong S C. 1998. Determination of available selenium and total selenium in soil by hydride generation atomic fluorescence spectrometry with non dispersive atomic fluorescence spectrometry[J]. Chinese Journal of Soil Science，29（1）：47.]

孫欣，李洪艷，李曉鵬，房經(jīng)貴，文仁德. 2015. 廣西葡萄一年兩收栽培模式及栽培技術(shù)[J]. 中外葡萄與葡萄酒，（3）：45-50. [Sun X，Li H Y，Li X P，F(xiàn)ang J G，Wen R D. 2015. Two harvest cultivation mode and cultivation techniques of Guangxi grape[J]. Sino-Overseas Grapevine & Wine，（3）：45-50.]

王海波，王孝娣，毋應(yīng)龍，王寶亮，謝計(jì)蒙，馮曉宇，劉鳳之. 2011. 設(shè)施葡萄對(duì)硒的吸收運(yùn)轉(zhuǎn)及積累特性[J]. 果樹學(xué)報(bào)，（6）：972-976. [Wang H B，Wang X D，Wu Y L，Wang B L，Xie J M，F(xiàn)eng X Y，Liu F Z. 2011. Selenium absorption，distribution and accumulation in the protected culti- vated grapevine[J]. Journal of Fruit Science，（6）：972-976.]

王忠. 2015. 施硒方式對(duì)葡萄吸收轉(zhuǎn)化硒的影響[J]. 農(nóng)業(yè)與科技，35（4）：22-23. [Wang Z. 2015. Effect of selenium on absorption and transformation of selenium in grape[J]. Agriculture and Technology，35（4）：22-23.]

張華華，康玉凡. 2013. 植物吸收和轉(zhuǎn)化硒的研究進(jìn)展[J]. 山地農(nóng)業(yè)生物學(xué)報(bào)，32（3）：270-275. [Zhang H H，Kang Y F. 2013. The advance in the uptake and transformation of selenium of plants[J]. Journal of Mountain Agriculture and Biology，32（3）：270-275.]

張潰珍，劉煥. 2015. 赤霉素對(duì)‘摩爾多瓦葡萄花序長(zhǎng)短及果粒大小的影響[J]. 中國(guó)園藝文摘，（11）：61. [Zhang K Z，Liu H. 2015. Effects of GA on the length and size of fruit grain of Moldova grape[J]. Chinese Horticulture Abstracts，（11）：61.]

周詠梅，黃羽，林玲，時(shí)曉芳，韋榮福，張瑛. 2015. 廣西葡萄育種研究進(jìn)展及對(duì)策[J]. 南方園藝，26（6）：49-51. [Zhou Y M，Huang Y，Lin L，Shi X F，Wei R F，Zhang Y. 2015. Research progress and countermeasures of grape breeding in Guangxi[J]. Southern Horticulture，26（6）：49-51.]

朱麗琴，魏欽平，許雪峰，韓振海，王小偉，李天忠. 2007. 葡萄對(duì)硒的吸收、分布和積累特性的初步研究[J]. 園藝學(xué)報(bào)，34（2）：325-328. [Zhu L Q，Wei Q P，Xu X F，Han Z H，Wang X W，Li T Z. 2007. Selenaxn absorbtion distribution and accumulation in grapeviie[J]. Acta Horticulturae Sinica，34（2）：325-328.]

（責(zé)任編輯 思利華）