我国粮食安全评价指标体系的构建与实证研究 胡令 朱荣花摘要：从粮食产销链各环节出发，综合考虑粮食生产安全、粮食流通安全及粮食消费安全3个方面，构建了我国粮食安全综合评价指标体系，进而对我国粮食的安全运行现状进行了综合评估及详细分析在此基础上，利用投影寻踪模型定量分析了我国2008—2017年期间的粮食安全问题，发现我国目前的粮食安全仍然存在诸多风险，10年间粮食安全问题逐年凸显且安全问题有明显的地域特征，最后为我国粮食安全问题的解决提供意见和建议关键词：粮食安全;产销链;双链量子遗传算法;投影寻踪：F326.11 文献标志码： A：1002-1302（2019）20-0316-06粮食安全由联合国粮食及农业组织（FAO）在世界粮食首脑会议上第1次提出粮食安全概念的提出主要是为了确保广大群众在任何时间任何生活状态下都可以享有充足的食品供给，并得到健康安定的生活[1]对于中国这样的泱泱大国，有效地保障粮食供给和粮食安全是人民安定、社会稳定和经济发展的基础，同时也是中国作为FAO成员国的一项社会责任改革开放四十余年来，我国农业科技水平稳步提升，粮食产量基本可以实现供需平衡但随着工业污染、水污染、土壤污染和农药、化肥残留物的增加，导致市面上不断涌现出不安全粮食，致使粮食安全问题在社会上凸显出来。
我国依照《国际食品法典标准》出台的《食品中农药最大残留限量》就是为了降低农药的使用数量，从而在整体维度上来提升粮食的安全等级由此可见，我国粮食安全的问题尚在解决途中，还需要完善的相关法律法规以及政策的支持[2]针对粮食安全问题，国内外学术界已有大量研究，其中多集中于粮食安全评价方面，主要分为定量评估及定性评估2种方式关于定性评价的研究，Carter等通过统计各发育阶段的儿童生长情况，即营养摄取情况，间接反映本国粮食安全状况，概括了国家粮食安全的影响因素为粮食供需情况、国家人口粮食产量和某一地区人均粮食产量[3];Coates从粮食生产面积、产量和人均粮食拥有量入手，主要分析了影响粮食产量的内在因素[4];王禹通过梳理国内外相关文献得出，影响我国粮食安全的主要因素为市场政策、粮食贸易结构、粮食供求和粮食生产环节等4个方面[5];马述忠等通过研究全球各国粮食产量及安全情况，梳理了我国粮食安全所面临的挑战并做出综合分析，总结指出，粮食安全受贸易、流通和储备等维度的影响[6];居占杰根据联合国粮食及农业组织关于粮食安全的定义和国内外学者的研究成果，结合我国实际国情认为，资源稀缺、种粮相对收益低、科技贡献率不高、自然灾害、人口增长、国际贸易等是影响我国粮食安全的主要因素[1];张元红等构建了包括粮食供给、区域分配、粮食销售、使用效率、可靠性、保障、调控能力及持续性等多个层面的指标体系，对我国粮食安全保障的现状、趋势、问题进行了分析[7]。
关于定量评价的研究，Acevedo以粮食自给供给效率、产量的波动比率、储备率及粮食自给效率等构建了相应的安全评估体系，通过加权平均法计算出各指标对粮食安全影响大小，根据权重不同逐次分析各指标对粮食安全的影响情况[8]杨学利运用客观赋值法对影响粮食安全的生产、经济、社会、资源环境和技术进行评价[9]单哲等运用柯布——道格拉斯生产函数分析了粮食总、单产与农业机械投入量、有效农业灌溉面积、化肥施用量和粮食播种面积等指标对粮食安全综合评价的影响[10]李光泗等以政府视角运用主要成分分析方法，通过收据的量表数据，综合考察了粮食生产、贸易、储备、行业竞争以及消费趋势对我国粮食安全的影响，进一步为粮食安全的宏观调控提供参考依据[11]杨磊等运用熵权法、二阶模糊综合评价法，分析粮食生产安全、消费安全和流通安全对粮食安全的影响[12]通过梳理国内外相关文献发现，涉及粮食安全评价的定性分析沉淀出的指标较为主观化，指标赋值缺少客观科学依据;而定量研究的方法也存在较多不足，一方面影响因素指标划分不够细致全面，大多数停留在二级指标或深入至三级指标但不够细化;另一方面运用的赋值权重分析方法有待改进鉴于现有研究基础，本研究通过信息沉淀评价指标构建粮食安全综合评价体系，使用双链量子遗传算法对投影寻踪模型进行优化，定量分析我国粮食安全问题影响因素，旨在构建客观科学的粮食安全评价指标体系，为解决我国粮食安全问题提供可行依据。
1 粮食安全评价体系的影响因素Food security是在20世纪中期引入我国的，当时由于人们消费水平限制，主要粮食就只是豆类、麦类、粗粮类，这些粮食正是解决温饱问题的根本考虑到当时国情，Food security被翻译成为粮食安全如今大部分地区的温饱问题都已解决，从追求温饱的消费模式转变为追求营养均衡、健康安全的消费模式然而目前我国各个地区之间的差异仍然存在，许多偏远山区和贫困乡村的吃饭问题仍未得到根本解決，许多人仍然处于营养不良状态因此，粮食安全不仅仅须要关注粮食生产的安全，还应该对从粮食生产到最终流入市场再进入消费者手中这一全过程进行综合分析与评价基于此，本研究从粮食产销全过程的视角出发，将粮食生产安全、粮食流通安全和粮食消费安全作为影响我国粮食安全问题的3个关键性维度进行分析，并针对性的提出相关建议1.1 生产安全对于人们而言，粮食是维持生命的基本保障，占据着重要的地位保证粮食的产量，维持自给自足的粮食生产水平，是确保粮食生产安全的根本粮食的综合生产能力受到土地面积、土地产量、生产技术、投入资金、劳动力人口以及自然灾害等诸多客观因素的影响此外，政府财政支持力度也是影响粮食生产的重要因素，即通过财政拨款和最低限价政策，可以有效调节粮食产量。
本研究中的粮食生产安全指标主要包括粮食生产面积、生产资料价格指数率、人均粮食产量、灌溉面积比率、生产面积成灾率、财政拨付指数、单位面积劳动力数量、单位面积施肥量等1.2 流通安全粮食流通是联系粮食生产商和消费者的通道，流通环节的效率、成本直接影响着粮食安全高效、便捷的运输流通环节能够使粮食快速运达到消费者手中，这既保证了市场需要也降低了运输成本，实现了商品向货币的转换粮食流通的购、销、存、运等环节同样受到诸多因素影响，本研究中的粮食流通安全指标包括流通成本浮动率、粮食净进口量、粮食流通价格指数浮动比、粮食商品化比重、粮食产销比重、流通事业费比例等1.3 消费安全保证消费者营养均衡、健康安全的粮食消费是粮食保障的根本目的粮食消费主要包括2种：一种是直接消费，即消费者为满足自身生存需求而进行的粮食消费;另一种是间接消费，主要用于饲料、工业农业生产等随着生活水平的提高，消费者对粮食的直接消费量比例逐渐降低，间接消费量比例不断增长，工业、农业对粮食的需求量不断增加本研究中的粮食消费安全指标包括人均粮食消费量、人均粮食供给量、自给率、消费价格指数的浮动比率、粮食储备率、营养（蛋白质）摄取量等2 粮食安全影响因素指标体系与模型构建2.1 粮食安全影响因素指标体系构建粮食安全受多方面因素影响，本研究从粮食产销链的全过程出发，从粮食流通、生产及消费等多重视角，对粮食安全的影响因素进行分析。
同时根据前文文献综述归纳梳理，在现有科学量表的基础上，结合粮食产销链中粮食安全所涉及的各环节可能出现的风险，构建科学的粮食安全体系该体系包含的二级指标包括粮食流通、生产及消费安全;三级指标主要包括20个指标具体如表1所示2.2 粮食安全影响因素模型的构建2.2.1 投影寻踪模型1972年，美国科学家Kruscal首先使用投影寻踪的研究方法，将高维的数据处理成为低维数据，进而研究数据的结构特征、聚类程度等[13]1974年，Friedman等用数据的一维特征构造了投影指标，成功分析了莺尾花聚类问题等的一系列高维数据，并将此研究方法命名为投影寻踪[14]此后关于投影寻踪这一方法的研究和应用引起了广泛关注在研究高维数据方面，投影寻踪法通过低维投影数据的特征来研究高维数据的特征，在评价、预测等领域都有一定应用其建模过程主要包括以下3个步骤2.2.1.1 指标数值的归一化根据统计归纳得来的数据将评价指标分为低优指标和高优指标，分别进行处理，设包含m个指标、n个样本集的原始数据集为[A（i，j）]nm，通过下列公式，统一统计样本数据变化范围，使得数据量纲相同，为后续处理奠定基础对于高优指标（越大越优的指标）：对于低优指标（越小越优的指标）：式中：maxAj和minAj分别为第j个指标在所有n个样本集数据中的极大值和极小值;ai，j为经过公式处理后的量纲归一化指标数据。
2.2.1.2 投影函数的构造2.2.1.3 优化投影指标函数当各项统计指标样本集合一定时，投影函数Q（λ）只随着投影向量λ的变化而变化当投影指标函数取最大值时，则为该样本数据的最优投影方向，最能显示出高维数据某一方面特征，因此目标函数为2.2.2 双链量子遗传算法投影寻踪也有其局限性，在使用投影寻踪法解决高维数据最优化问题时，容易陷入局部最优而非全局最优对此本研究通过双链量子遗传算法对投影寻踪模型进行改进，进而对粮食安全各项指标进行评价对于量子遗传算法，它基于量子态及量子比特双重基础，在遗传编码中将量子态的向量注入，从而使得算法更加优化染色体双链上的2个同位基因分别代表1个量子位，每条染色体单链都代表1个解本研究引入的双链量子遗传算法，使投影寻踪优化过程中目标函数与投影寻踪模型一致需要优化的变量即各个样本数据投影值的具体求解过程如下2.2.2.1 量子编码2.2.2.2 变换解空间2.2.2.3 量子染色体的更新完成上述步骤后，基于量子旋转门的旋转矩阵来对量子比特相位进行更新，从而达到相应的目的这一步骤的目的是使得每个染色体都接近最优，从而实现种群的进化量子旋转门变换矩阵是一个可逆矩阵，通常采用如下矩阵：式中：θ为旋转角度。
由式（14）可知，经过量子门的旋转变换，染色体的长度并没有发生变化，而θ的取值关系到算法收敛的速度，通常θ取值范围为0.001π≤θ≤0.005π量子旋转角度示意如图1所示根据目标函数的变化率，采用的转角步长函数为2.2.2.4 量子染色体变异为了增加样本的多样性，降低过早收敛的概率，将双链量子遗传算法利用量子门对染色体进行变异为实现染色体双链同步变异，采用交换量子位对概率幅（η，ξ）T进行置换变异过程如下：通过双链量子遗传算法和投影寻踪模型，计算得出各个指标的最佳投影方向向量，并将其作为各个指标的权重而权重的大小，可在一定程度上表示指標的重要度基于权重对指标进行有序排列，依次对粮食安全影响因素进行分析，并提出相应的对策和建议3 粮食安全评价实证研究3.1 数据来源实证研究应遵循真实有效原则因此，为确保研究结果的可靠性，本研究构建的粮食安全指标体系中20个三级指标的数据均来源于权威机构出具的统计年鉴各三级指标面板数据来源如表2所示3.2 数据处理本研究采用双链量子遗传投影寻踪模型对粮食安全评价体系的20个三级指标的原始数据进行预处理由表1可知，粮食安全评价指标体系中的20个三级指标对粮食安全产生的作用包括正向作用和反向作用2种情况。
因此很难对不同作用趋势指标进行比较加之，各个指标直接量纲差距较大，进一步加剧了数据分析难度所以，在分析评价之前要对原始数据进行无量纲化处理，使各个指标可以被评价数据无量纲处理后，采取双链量子遗传投影寻踪模型消除指标差异性，得出粮食安全评价指标的投影值，以此为前提分析其时序变化特征设定样本维数（二级指标数）=3个;三级指标=20个;种群=100个;最大迭代次数=200次由此可得粮食安全时序变化最佳投影向量b，各三级指标赋值权重如下：与此同时，能够获得2008—2017年我国粮食安全各子系统最佳投影方向向量（表3）以及系统综合投影评价值（表4。