基于tcmgis的北苍术产地适宜性分析和定量评价

张玮玮1袁2袁张志鹏3袁郑司浩4袁蔡翠芳1袁2（1.山西药科职业学院，山西太原030031；2.山西中医药大学中药学院，山西晋中030600；3.山西国新晋药集团道地药材经营有限公司，山西太原030032；4.中国中药有限公司，北京100195）

摘要院利用中药材产地适宜性分析地理信息系统（TCMGIS）对311个北苍术野生标本地理分布信息进行研究，

北苍术适宜分布区主要位于我分析北苍术适宜产地空间分布状况，寻找适宜北苍术生长特性的区域。结果表明，生态相似系数100%区域分布在28个省区的1185个县市，国温暖湿润的东部季风区北部，适宜产区面积总和为陕西、辽宁、内蒙古、甘肃、黑龙江等省区。研究结果对北1055632.4km2，生态适宜区主要集中在河北、河南、山西、苍术的野生资源保护、人工引种栽培和资源可持续利用具有重要意义。生态区划关键词院北苍术；TCMGIS；产地适宜；相似系数；中图分类号院S567.21＋1

文献标识码院A

文章编号院1002-2481（2020）02-0208-05

QuantitativeEvaluationandRegionalSuitabilityAnalysisof

byTCMGIS

22，ZHANGZhipeng3，ZHENGSihao4，CAICuifang1，ZHANGWeiwei1，Taiyuan030031，China；2.CollegeofTraditionalChineseMedicine，Shanxi（1.ShanxiPharmaceuticalVocationalCollege，

UniversityofChineseMedicine，Jinzhong030600，China；3.ShanxiGuoxinLocalPharmaceuticalGroupGenuineMedicineCo.,

Ltd.，Taiyuan030032，China；4.ChinaNationalTraditionalChineseMedicineCo.，Ltd.，Beijing100195，China）Abstract：Inthisresearch,itdetermineddistributionareasandrangesfrom311samplesusingtheGeographicInformationSystemforTraditionalChineseMedicine（TCMGIS）,toprovidereferenceofthescientificplantingandproductionareaplanningandpredictnewgrowthareaforAtractylodeschinensis.TheresultsshowedthatAtractylodeschinensismainlydistributedinthenorthernofEastmonsoonareaofChina.Atotalof1185countiesin28provinces,coveringalandareaof1055632.4km2,sharedsimilarwithecologicalfactorswiththoseofnativehabitatsappropriateforAtractylodeschinensis.EcologicallysuitableareasweremainlyconcentratedinHebei,Henan,Shanxi,Shaanxi,Liaoning,Neimenggu,GansuandHeilongjiangprovinces.Itwillbeofbenefitfortheconservationofwildresources,artificialcultivationandsustainableutilization.

Keywords：Atractylodeschinensis;TCMGIS;regionalsuitability;similaritycoefficient;ecologicalregionalization

北苍术（Atractylodeschinensis（DC.）Koidz.）是菊科苍术属多年生草本植物，生于山坡草地、林下、灌丛及岩缝隙中[1-2]。其干燥的根茎是我国传统中药苍术的主要植物来源，具有燥湿健脾、祛风散寒的北苍术具有多种活性成分，功效[3]。药理研究发现，

抗炎和降血糖等作用[4-6]。有抗溃疡、保肝、抗肿瘤、随着大健康产业的发展，医药市场对药材苍术的需

求逐年增加，但由于以野生资源为主的药用苍术被连年过度采挖，以致资源储备量急剧下降[7]。为解决日益扩大的市场需求量，近年来加大了对其引种驯

东北部分地区开化和栽培繁育的工作，我国华北、始有了零星的人工栽培。

收稿日期院2019-08-26

（GZY-KJS-2018-004）基金项目院国家中医药管理局全国中药资源普查项目

中药苍术的另一基原植物为茅苍术Atracty原lodeslancea（Thunb.）DC.，近年来，有学者对其进行了一定的研究。陈士林等[8]对茅苍术的适宜种植区进行了系统研究；郭兰萍等[9]和ZHU等[10]对茅苍术的化学成分与产区生态特征相关性进行了研究。但关于北苍术种植适宜性方面的研究鲜有报道。

道地药材的形成要求合适的自然生境条件（温

[11]降水等）。本研究应用中药材产地适宜度、湿度、

（TCMGIS）性分析地理信息系统，依据气候和土壤

因子相似原理，确定北苍术不同的生态相似度的适

引种驯宜产地[12-13]，以期为北苍术野生资源保护、进而化、筛选适宜种植区、扩大种植面积提供依据，

作者简介院张玮玮（1981-），女，山西应县人，讲师，硕士，主要从事药用植物资源及规范化栽培教学与研究工作。蔡翠芳为通信作者。

窑208窑

张玮玮等：基于TCMGIS的北苍术产地适宜性分析和定量评价

解决北苍术资源短缺及资源可持续发展的问题。

1

材料和方法

1.1

样点数据来源

北苍术分布点数据主要来自3个部分：

一是标本库，包括中国数字植物标本馆（http：//www.cvh.org.cn.）；二是文献，包括《中国植物志》、《中国高等植物图鉴》以及北苍术相关研究报道；三是部分实地调查数据，来自山西药科职业学院中药系师生对其的长期的实地考察。后续经过删除重复及信息不完整的数据，确定野生北苍术的采集地点共311个，主要集中分布在山西、河北、河南、内蒙古、辽宁、陕西等北部地区，湖北、安徽、江苏、云南等南部地区也有少量分布。1.2

数据库来源

TCMGIS是以地理信息系统（GIS）为平台，以气候因子、土壤、基础地理信息及第三次全国中药资源普查等数据库为支撑；其中，气候因子包括活动积温、年均日照时数、年均气温、年均湿度、年均降水量、1月平均气温、7月平均气温、最冷月和最热月的极端气温等指标，气候因子数据来源于全国气象站点30a地面气象数据，进行栅格数据计算1km2）

及衍生数据；土壤数据以全国土壤普查的数字化土壤类型数据为基础[8]。1.3

TCMGIS分析方法

根据311个野生北苍术样本点经纬度信息，通过TCMGIS系统提取标本点的生态因子[14]，分析确定影响植物分布的关键性生态地理因子，进行栅格数据分析并形成标准化的栅格数据，然后将栅格数据转化成矢量数据文件，进而叠加在全国行政区划图层上，获得各行政单位内不同等级生态相似度区域的相关信息[8]，最终运用Excel2010软件对数据进行分析。其中，栅格数据分析中对样本点生态数据采用绝对偏差均值标准化方法进行数据标准化处理。

Sf＝1n（x1f－mf＋x2f－mf＋…＋

xnf－mf）

（1）

其中，Sf是绝对偏差均值，mf＝1n（x1f＋x2f＋…＋

xnf），mf表示n个聚类对象第f个要素的均值；标准化值为Zif＝xif－Smff归一化到0～100；Zif是第i个

因子的标准化值[15]。

利用欧式距离公式dif＝

＝1aixik－xjk2标准化后生成的栅格数据进行蓘移

nk处理，

每个生态蓡12对

因子赋予等权，即ai＝1n（n为因子数目），

计算每个栅格与目标值的距离[8]。

栅格新值＝A×权重＋B×权重＋……＋E×权重

（2）

根据式（2）分析完成重分类的各个因子图层数

据，得到新的栅格值。

2

结果与分析

2.1

北苍术适宜因子范围

经过TCMGIS系统对311个野生北苍术标本

采集点生态因子的分析，得到适宜北苍术生长的生态因子范围为：年平均气温-5～22.2℃；7月最高气温32.7℃；1月最低气温-35.8℃；7月平均气温13.4～28.6℃；1月平均气温-30.9～-12.9℃；年均相对湿度48%～81.6%；年均降水量27.9～157.3mm；年日照时数1093～3027h；≥10℃活动积温为

1194～7634.8℃；适宜的土壤类型主要为壤土、褐土、黑垆土、粗骨土、潮土、沼泽土。2.2

北苍术适宜产地分布

利用TCMGIS系统分析得出北苍术生态适宜

分布区（图1），其主要位于我国温暖湿润的东部季风区北部，该区域受季风影响显著，

气候湿润，雨热同期，普遍高温多雨，冬季受北方冷气流影响，大部分地区寒冷干燥。北苍术主要位于华北暖温带野生家种中药区、内蒙古中温带野生中药区和西北中温带暖温带野生中药区[16]。本研究将生态因子相似系

数为100%的区域定义为最适生长区、98%～100%的区域为次适宜区，其他为不适宜区域。

窑209窑

（山西农业科学2020年第48卷第2期

2.3北苍术生态相似度100%的最适生长区

利用通过TCMGIS系统分析得出北苍术生态

（图1）相似系数100%最适生长区分布图，图中深绿色为野生北苍术生态相似系数100%区域。北苍术

生态相似数100%的最适生长区域较为广泛，总面（表1）共分布在1185个县区，积达1055632.4km2，

（137150.0km2）（102808.8km2）其中，河北、河南、

（956.6km2）（986.8km2）辽宁、山西分布较多，分别占全国生态相似系数100%区域总面积的

云南、13%、9.7%、9.3%、9.1%（图2），内蒙古、陕西、在福建、湖黑龙江、湖北、四川、甘肃等省偶有分布，南、宁夏、浙江也有零星分布。

表1北苍术不同生态相似系数区域分布情况

区）省（市、安徽

北京福建甘肃广东广西贵州河北河南黑龙江湖北湖南吉林江苏辽宁内蒙古宁夏青海山东山西陕西上海四川天津云南浙江重庆总计

生态相似度100%区域

县（市）56111214284761591206161332575011041026670121587291185

面积/km224755.816274.41696.534175.94217.913582.95262.9137150.0102808.857238.850686.0583.12328.8246.2986.8819.6186.491806.56.683914.1-39331.45742.611172.9-211.5

34.42845.31055632.4

占县（市）总面积比例/%

27.492.620.347.59.211.119.469.965.612.634.98.49.247.261.420.02.7

60.465.856.318.650.036.123.21.39.8

县（市）53513394963211199435219392742131916871317287322121128

生态相似度98%～100%面积/km2246.5887.6599.810815.914633.116550.83275.332628.230999.93618.335316.22204.017036.7133.937680.55959.7218.245073.423399.822768.9556.27976.445.53599.812361.615299.23253.334.6399716.5

占县（市）总面积比例/%

22.714.84.311.518.311.47.331.926.84.222.94.3

29.120.516.829.91.337.120.930.829.23.2.91.19.37.417.816.7

2.4北苍术生态相似度98%耀100%的次适宜区（图3）系数98%～100%和次适宜分布图，图中深绿色部分是生态相似度100%的区域，浅绿色是生态

通过TCMGIS系统分析得出北苍术生态相似

窑210窑

张玮玮等：基于TCMGIS的北苍术产地适宜性分析和定量评价

相似度98%～100%的区域。

由图3和表1可知，生态相似系数98%～100%的适宜生长区主要集中于我国28个省市的1128个县区，

面积达399716.5km2。其中，主要分布在山东、内蒙古、

湖北、河北等省，包括119个县（市），面积达32628.2km2，占所在县市总面积的31.9%；其次为河南、安徽、

山西等省。山东、内蒙古、湖北、河北、河南、安徽、山西次适宜区面积分别为45073.4、37680.5、35316.2、32628.2、30999.9、246.5、23399.8km2。

3结论与讨论

TCMGIS系统以地理信息系统平台（GIS）作为

分析平台，以影响药用植物生长的10个主要生态因子作为数据基础，通过计算机软件对该植物适宜产地分布状况进行基础评估[16]，目前至少已完成了

210种药材产地生态适宜性分析[8，17-19]，该系统作为

解决中药材适宜产地评价分析工具，操作简便、

适用性强，而且结果客观合理。本研究以北苍术为分

析对象，其在文献中记载的药材主产地是内蒙古、山西、河北、吉林、辽宁、黑龙江；

陕西、山东、甘肃等地亦产；原植物分布于东北、华北地区及陕西、宁夏、甘肃、山东、

河南等地[20]。本研究中北苍术生态相似度100%最适生长区的主要分布省份与文献记

载、资源普查等记载基本一致。由于样点数据中包含有自然分布于我国湖北、

江西、云南等南方地区的标本信息，研究结果表明，我国南方部分地区也适合北苍术的生长，尽管其与普遍认为的北苍术产于北方的认识不一致，但本研究中的北苍术是植物

学概念，而“北苍术产于北方”中的概念是基于商品学[1]。此外，研究中查阅到北苍术在湖北、

江西、云南等南方地区分布的标本信息及适宜在南方部分地

区生长的结果，也一定程度上佐证了《中国植物志》将北苍术Atractylodeschinensis（DC.）Koidz.并入Atractylodeslancea（Thunb.）DC.，并将其作异名处

理[1，21]。因此，TCMGIS系统分析结果是可靠的，客观合理地对北苍术进行种植区划，为引种驯化、扩大人工种植区提供了理论依据。

生态因子是道地药材形成的重要因素[22]，盲目种植会对药材品质造成很大的影响[23-24]。通过采用TCMGIS系统对北苍术的产地适宜性进行分析，得出北苍术的适宜产区分布在28个省区的1185个县市，适宜产区面积总和为1055632.4km2，生态适宜区主要集中在河北、河南、山西、陕西、辽宁、内蒙古、甘肃、黑龙江等省区。虽然TCMGIS系统从理论上通过综合生态因子针对北苍术进行统计分析后得到适宜产区分布区域，为北苍术药材的生产提供了科学依据，然而植物个体的生长特性一般都有一定的生境或小气候要求，北苍术自然生长的生境多为山坡草地、林下、灌丛及岩缝隙中[2]，因此，要求在具体的工作中应以TCMGIS系统分析结果作为引种的基本依据，经过实地调研、田间预试验、药材品质检测评估等研究，并结合当地经济、文化、历史等因素对产区的可用性进行综合测评，

才能制定出更好的适宜北苍术的栽培种植产区区划及规范化种植计划。

中药苍术的市场供应均为野生资源，随着中药市场对原料药需求量的加大，

其野生资源受破坏程度较大，资源量逐年减少，因此对北苍术进行野生变家种、仿野生栽培是未来发展的必然趋势，合理布局发展北苍术生产、确定其区划以避免盲目引种、保证北苍术药材品质显得尤为重要。本研究通过TCMGIS系统分析得到了与野生北苍术标本分布点生态因子相似的区域，

该区域可作为最佳种植区选取和新产区预测的基本依据，进而解决北苍术资源可持续发展的问题。

参考文献院

1］中国科学研究院植物研究所.中国高等植物图鉴（第四册）[M].

北京：科学出版社，1994：602.

2］中国科学院中国植物志编辑委员会.中国植物志[M].北京：科

学出版社，1959：192-196.

3］中华人民共和国药典委员会.中华人民共和国药典：一部[M].

北京：化学工业出版，2015.

4］邓爱平，李颖，

吴志涛，等.苍术化学成分和药理的研究进展[J].中国中药杂志，2016，41（21）：3904-3913.

5］KOONRUNGSESOMBOONN，NA-BANGCHANGK，KARBWANG

窑211窑

［［［［［山西农业科学2020年第48卷第2期

J.TherapeuticpotentialandpharmacologicalactivitiesofAtracty原lodeslancea（Thunb.）DC[J].AsianPac.J.TropMed.，2014，7（6）：421.

［6］KIMURAY，SUMIYOSHIM.EffectsofanAtractylodeslancearhi-zomeextractandavolatilecomponent茁-eudesmolongastrointesti-nalmotilityinmice[J].J.Ethnopharmacol，2012，141（1）：530.郭利霄，等.北苍术价格波动及其影响因素分［7］温子帅，王建岭，

析[J].农村经济与科技，2017，28（17）：156-157.

［8］陈士林.中国药材产地生态适宜性区划[M].北京：科学出版社，

2011.

［9］郭兰萍，黄璐琦，阎洪，等.基于地理信息系统的苍术道地药材

气候生态特征研究[J].中国中药杂志，2005（8）：565-569.［10］ZHUSD，PENGHS，GUOLP，etal.RegionalizationofChinese

materialmedicalqualitybasedonmaximumentropymodel：AcasestudyofAtractylodeslancea[J].ScientificReports，2017（7）：42417.

王珍，等.濒危高原植物羌活化学成分与生态［11］黄林芳，李文涛，

因子的相关性[J].生态学报，2013，33（24）：7667-7678.曾凡琳，等.药用和观赏植物常山产地适宜性［12］张娜娜，谢彩香，

分析和定量评价[J].植物科学学报，2016，34（2）：230-237.［13］全健，索风梅，谢彩香，等.使君子药材产地生态适宜性研究

[J].世界科学技术—中医药现代化，2014，16（2）：339-345.索风梅，等.人参药材中人参皂苷与生态因子［14］贾光林，黄林芳，

的相关性及人参生态区划[J].植物生态学报，2012，36（4）：302-312.

任德权.道地药材的产地适宜性研究[J].中国［15］索风梅，陈士林，

中药杂志，2005（19）：9-12.

［16］中国药材公司.中国中药区划[M].北京：科学出版社，1995.陈君，等.基于野生银柴胡的产地适宜性分析［17］马伟宝，谢彩香，

[J].中国现代中药，2017，19（5）：684-687.

关小羽，等.娑罗子原植物产地生态适宜性研［18］石召华，谢彩香，

究[J].中国现代中药，2013，15（4）：303-307.

［19］张海龙，陈乐.基于生态位模型的连翘产地适宜性定量分析[J].

山西农业科学，2017，45（8）：1321-1324.

编委会.中华本草：［20］国家中医药管理局（中华本草）第二十一卷

[M].上海：上海科学技术出版社，1999：709.

［21］石铸.关于苍术植物的学名问题[J].植物分类学报，1981，19

（3）：23-29.

张卫，等.谈［22］梁飞，李健，“道地药材”的形成原因[J].中国中药

杂志，2013，38（3）：466-468.

李晓娟，等.黄芩生态适宜性评价及生态因子［23］周国富，刘金欣，

对5种主要指标性成分的影响[J].中国实验方剂学杂志，2016，22（20）：28-32.

［24］卢燕林，张强，李彩萍，等.地道潞党参营养特性与产地土壤特

征研究[J].山西农业科学，2007，35（5）：42-45.

渊上接第197页冤

［5］李隆.间套作强化农田生态系统服务功能的研究进展与应用展

望[J].中国生态农业学报，2016，24（4）：403-415.

刘祖贵，等.间作种植模式对玉米和大豆干物质［6］高阳，段爱旺，

积累与产量组成的影响[J].中国农学通报，2009，25（2）：214-221.

白伟，等.玉米/花生间作系统作物产量及根［7］高砚亮，孙占祥，

系空间分布特征的影响[J].玉米科学，2016（6）：79-87.［8］王根全，王节之，郝晓芬，等.谷子轻简化栽培技术研究[J].山西

农业科学，2015，43（11）：1438-1441.

［9］方萍.玉/豆间作对菜用大豆干物质积累、产量及品质的影响

[D].雅安：四川农业大学，2016.

［10］雍太文，刘小明，宋春，等.种植方式对玉米/大豆套作体系中

作物产量、养分吸收和种间竞争的影响[J].中国生态农业学报，2015，23（6）：659-667.

［11］刘朝茂，李成云.玉米与大豆间作对玉米叶片衰老的影响[J].

江苏农业学报，2017，33（2）：322-326．

［12］AD魺NISMOREIRA，LARISSAALEXANDRACARDOSOMORA

ES，G魻TZSCHROTH.Copperfertilizationinsoybean-wheatin-tercroppingunderno-tillmanagement[J].Soil&TillageResearch，2019，193：71-84.

［13］MUHAMMADALIRAZA，LINGYANGFENG，WOPKEWERF，

etal.Narrowwiderowplantingpatternincreasestheradiationuseefficiencyandseedyieldofintercropspeciesinrelayintercroppingsystem[J].FoodandEnergySecurity，2019，8（3）：44-46.李进，等.不同氮水平下夏玉米夏大豆间作对［14］孔玮琳，薛燕慧，

其农艺性状及产量的影响[J].山东农业科学，2018，50（7）：116-120.

陈志权，等.林草间作系统土壤微生物及其活［15］白昌军，虞道耿，

性研究[J].热带作物学报，2015，36（5）：841-846.

杨萌珂，等.不同氮水平下夏玉米夏大豆间作［16］焦念元，宁堂原，

对其农艺性状及产量的影响[J].生态学报，2013，33（14）：4325-4330.

王锐，等.玉米－大豆带状间作对大豆生长、［17］范元芳，刘沁林，

光合荧光特性及产量的影响[J].核农学报，2017，31（5）：972-978.

吕娇.大豆/玉米间作下干物质累积动态变［18］李淑敏，乔海涛，

化与产量分析[J].天津农业科学，2011，17（4）：38-41.［19］张蕙琪，王宏富，王彦雯.氮磷钾配施对谷子农艺性状及光合

特性的影响[J].山西农业科学，2018，46（1）：42-49.

［20］郜庆炉，宋留轩，杨怀森，等.不同种植方式对玉米大豆间作群

体产量形成的影响[J].河南农业科学，1996（6）：6-9.

窑212窑