甘肅道地藥材“岷歸”產量預測模型對比研究

摘要：本文基于《甘肅道地藥材志》等文獻中甘肅當歸和岷歸的產量數據，計算并對比了甘肅當歸和岷歸的畝產情況，對現有文獻中關于岷歸產量的預測模型進行了歸納和總結，并對預測模型的精度進行了述評。研究表明，岷縣作為甘肅省當歸的主要產區，始終承擔著當歸產量增長的重大作用。現有的岷歸產量預測模型不適用于預測當前的岷歸產量，特別是以年序為變量的趨勢產量預測模型具有較為嚴格的適用條件；當將年序拓展到其適用條件之外時，可能得到違背邏輯的計算結果。針對趨勢產量的后續研究中，需積極探索更為合理的變量，以達到較為準確的預測效果。

關鍵詞：甘肅道地藥材；岷歸；產量預測模型；趨勢產量；氣象產量

當歸為傘形科植物當歸Angelica sinensis（Oliv.）Diels的干燥根[1]，素有“婦科人參”或“婦科要藥”之稱，有補血活血、調經止痛和潤腸通便等功效[2]，常入當歸四逆湯[3]、七味抗體散[4]、補陽還五湯[5]等經典名方。當歸歷來以甘肅岷縣所產為最好[6]，故而岷縣當歸（以下簡稱“岷歸”）成為國家地理標志產品。

當前，栽培是岷歸當前最主要的種植方式，而岷歸的種植面積占全國當歸種植面積的80%以上。由于岷歸是當地農民主要收入來源，作為促進地方經濟發展和當地中醫藥產業的重要支柱[7]，是甘肅岷縣名副其實的“經濟引擎”，故而岷歸的產量研究是岷歸研究的重要內容之一。

早在1988年，路鵬南[8]即在基于岷縣當歸的地理條件分析的基礎上，提出了岷歸的產量預測模型。此后，嚴輝等[9]于2012年采用多項式預報法，也建立了岷歸的產量預測模型。鑒于岷歸產量的預測模型對于相關部門和機構掌握生產信息，提前謀劃生產、經營的管理較為重要，對堅定農戶種植信心也具有重要價值，而以上兩個預測模型已歷時已久，故本文擬以《甘肅道地藥材志》等文獻中甘肅當歸和岷歸的產量數據為基礎，通過對比分析，對甘肅當歸、岷歸的產量特征進行歸納總結，并對現有文獻中的當歸預測模型進行對比研究。

1 甘肅省當歸的產量及數據來源

在《甘肅道地藥材志》[10]中，給出了1978—2012年，

甘肅全省的當歸種植面積（萬畝）及總產量（萬t）的數據；而嚴輝等[9]則給出了岷縣中藥材產業發展局提供的1995—2009年岷縣當歸產量，其中實際單產為kg/hm2。

為對比方便，本文以kg/667 m2為單位，計算了《甘肅道地藥材志》[10]和文獻[9]中的667 m2產量。

計算發現，1978—2012年共35年間，甘肅省的當歸產量均值為133.16 kg/667 m2，而1995—2009年共

15年間，岷歸產量均值為157.16 kg/667 m2。此外，隨著年份的增長，盡管個別年份中甘肅當歸、岷歸的產量有所變化，但均呈現出隨年份而產量增加的趨勢。

1997年之前，甘肅省當歸的種植面積不斷波動，而1998年之后，則呈現出大幅增長的趨勢。與其對應的是，1997年之前甘肅省當歸的總產量則始終在1萬～2萬t，1998年之后則大幅攀升。

2 甘肅道地藥材“岷歸”的667 m2產量對比

產量是衡量當歸種植的重要參數之一，故得到了眾多學者的關注[11-12]。1954年前，岷縣當歸的產量100～125 kg/667 m2；將產量100 kg/667 m2與125 kg/667 m2進行平均可得，1954年之前的岷歸的平均產量約為112.5 kg/667 m2；而1961—1980年共20年間，岷縣當歸的平均產量117 kg/667 m2[8]。表明在1980年之前，岷歸的平均667 m2產量的變化幅度較小。

通過對比1995—2009年甘肅省當歸的667 m2產量與“岷歸”的667 m2產量可以發現，除1995年、1997年、1998年、2005年和2006年此5年之外，其余的10年中，全省的當歸產量與岷歸的產量均有一定差距。特別是1999年、2003年、2004年、2007—2009年共6年間，全省內當歸的667 m2產量均高于岷歸的667 m2產量。

而對于岷歸的畝產量，1961—1980年平均為117 kg/667 m2，1995年、1996年的數值分別為104.06 kg/667 m2、113.39 kg/667 m2；但1997年之后則始終保持在130.19 kg/667 m2之上，而2005年則高達217.39 kg/667 m2。因此，1995—2009年岷歸平均產量則達到了157.66 kg/667 m2。

此外，在2000年全省當歸畝產量驟降的情況下，岷歸的產量則降幅較小，由此可知，岷縣作為甘肅省當歸的主要產區，始終承擔著產量增長的重大作用。這從1978年、1979年全省當歸產量低于80 kg/667 m2，而1961—1980年岷歸平均產量為117 kg/667 m2上也可看出。

3 岷歸產量預測模型及其檢驗

3.1 現有的岷歸產量預測模型

為探究影響當歸產量的因素，建立其產量預測模型，路鵬南[8]將岷縣1961—1980年共20年間的7個因素納入研究范圍，在從實際產量中剔除了時間趨勢后回歸分析得到了當歸產量的計算式：

式中，T為4月上旬～10月下旬期間的平均氣溫，H為

相對濕度，R為年降水量，T7為七月平均最高氣溫。

嚴輝等[9]根據1995～2009年岷縣當歸產量，將趨勢產量YN與氣象產量YTR進行分離，得到當歸產量模擬預報模型的計算式。其中，對于趨勢產量YN，由于嚴輝等[9]以1995～2009年岷縣當歸產量為研究依據，故令1995年的N=1，1996年的N=2，……，2008年的N=14，給出了YN的計算式為：

王慧珍等[13]指出，溫度和降雨量是影響產量的關鍵生態因子。而對于氣象產量，嚴輝等[9]進一步用SPSS進行分析，得到了氣象產量的模擬預報方程。為更為清晰地表達，本文基于嚴輝等[9]中的符號說明，采用表示n月m旬的平均氣溫，采用Rn，m表示n月m旬的降水量。其中，對于“m旬”，分別采用f、s、t表示上旬、中旬和下旬。則嚴輝等[9]的氣象產量的模擬預報方程可表示為：

式中，為7月下旬平均氣溫，為8月上旬平均氣溫，為10月下旬平均氣溫，為5月下旬降水量，為6月下旬降水量，為7月中旬降水量，為10月下旬降水量。

3.2 岷歸產量預測模型的檢驗

對于路鵬南[8]的預測模型，將4月上旬—10月下旬的平均氣溫T=4.9℃、相對濕度H=68％、七月平均最高氣溫T7=16℃以及年降水量R=596.5 mm代入式（1）后，得到的Y=－11.28 kg。而對于嚴輝等[9]的預測模型中的趨勢產量，將年序N=1995～2025代入式（2）并繪圖可以發現，盡管在2012年之前岷歸667 m2產量的變化情況仍然與趨勢產量相符，但2012年之后則驟降，直至2019年及其之后成為負值（圖1）。

顯然，由于路鵬南[8]的預測模型、嚴輝等[9]的預測模型得到的畝產量出現了負值，故而這兩個岷歸的預測模型在當前的岷歸產量預測中將不再適用。

需要指出的是，影響岷歸產量的因素較為復雜，氣候因素、自然災害等均可較大程度地影響岷歸的產量。如根據嚴輝等[9]的預測模型，在1995—2009年共15年間，氣象產量為正值的僅有

5年，占全部年份的1/3；而氣象產量為負值的則有10年，占全部年份的2/3。其中，氣象產量的最大值是2005年的30.89 kg/667 m2，最小值則為2009年的—20.60 kg/667 m2。2009年減產的主要原因是當年甘肅岷縣發生了特大山洪[14]。

盡管路鵬南[8]的預測模型、嚴輝等[9]的預測模型在預測當前的岷歸產量上的適用范圍有限，但通過對比仍可得到如下啟示：

（1）隨著研究的深入，考慮的氣象因素越來越細致。如式（1）中僅包含4月上旬—10月下旬期間的平均氣溫和七月平均最高氣溫，而式（3）則包含了三個氣溫因素；又如式（1）中僅納入了年降水量，而式（3）則包含了三個降水量因素。

（2）以年序N為變量的趨勢產量預測模型具有較為嚴格的適用條件；當將年序拓展到其適用條件之外時，可能得到違背邏輯的計算結果。因此，針對趨勢產量的后續研究中，需積極探索更為合理的變量，以達到較為準確的預測效果。

4 結論

（1）1978—2012年共35年間，甘肅當歸的產量為133.16 kg/667 m2，1995—2009年共15年期間，岷歸的產量為157.66 kg/667 m2，均高于1961—1980年岷歸平均產量117 kg/667 m2。

（2）自然災害是造成當歸減產的主要因素之一。2009年岷歸產量和2012年甘肅當歸產量驟降均是由于產區發生了特大山洪或特大冰雹山洪泥石流災害。剔除自然災害的影響，氣象產量占實際產量的比值相對較小。

（3）以年序N為變量的趨勢產量預測模型具有較為嚴格的適用條件；當將年序拓展到其適用條件之外時，可能得到違背邏輯的計算結果。針對趨勢產量的后續研究中，需積極探索更為合理的變量，以達到較為準確的預測效果。

參考文獻

[1] 國家藥典委員會.中華人民共和國藥典：1部[M].北京：中國醫藥科技出版社，2015.

[2] 漆琚濤，許彩荷，紀瑛，等.當歸種子直播栽培對其產量和質量的影響研究[J].中藥材，2018，41（8）：1804-1808.

[3] 曾妮，吳承艷，王進.經典名方當歸四逆湯的藥物劑量考證[J].中草藥，2021，52（21）：6648-6655.

[4] 范榮，薛曉會.七味抗體散質量標準研究[J].中國藥房，2009，20（3）：215-216.

[5] 蔣健.補陽還五湯臨床運用發揮[J].時珍國醫國藥，2009，20（3）：755-756.

[6] 魏建和，屠鵬飛，李剛，等.我國中藥農業現狀分析與發展趨勢思考[J].中國現代中藥，2015，17（2）：94-98.

[7] 張玲.熱河黃芩栽培管理技術探討[J].新農民，2024（27）：67-69.

[8] 路鵬南.岷縣當歸的地理條件分析[J].地理學與國土研究，1988，4（4）：34-38.

[9] 嚴輝，段金廒，王瑞琪，等.岷縣當歸產量預報模型研究[J].中藥材，2012，35（4）：511-514.

[10] 宋平順，丁永輝，楊平榮.甘肅道地藥材志[M].蘭州：甘肅科學技術出版社，2016.

[11] YANG L，LI J，XIAO Y，et al.Early Bolting，Yield，and Quality of Angelica sinensis（Oliv.）Diels Responses

to Intercropping Patterns[J].Plants，2022，11（21）：2950.

[12] LIN H L U C，ZHUANG Q，FU H.Habitat niche-fitness

and radiXyield prediction models for Angelica sinensis

cultivated in the alpine area of the southeastern region of

Gansu Province，China[J].Plant Production Science，2008，11（1）：42-58.

[13] 王惠珍，張恩和，高素芳，等.海拔對岷縣當歸產量調控及關鍵影響因子分析[J].中草藥，2013，44（14）：1990-1994.

[14] 盧育霞，王偉鋒，石玉成，等.甘肅岷縣災害之地理解析[J].甘肅科學學報，2014，26（4）：70-76.