¹株洲市农业科学研究所，湖南 株洲

²株洲香之优农业科技发展有限责任公司，湖南 株洲

³湖南省水稻研究所，湖南 长沙

收稿日期：2022年6月8日；录用日期：2022年7月6日；发布日期：2022年7月19日

摘要

为筛选适合株洲市种植的优质、富硒、低镉的优质稻品种，引进15个优质稻品种进行随机区组试验，在灌浆初期施用生物纳米硒肥1次，硒肥用量(以Se计) 1.25 g/667m²，测定精米镉、总硒、有机硒含量。结果表明：1) 在重度镉污染土壤条件下(1.17 mg/kg)，15个优质稻品种对镉的吸收积累能力存在一定的差异，但未达到显著水平(P = 0.05)；其中，精米镉含量最低的是拾两优1817，说明该品种具有低镉积累特性。与不施硒肥相比，施用硒肥品种的精米镉含量平均降低3.05%，降镉效果不明显。2) 15个优质稻品种籽粒硒富集能力存在极显著差异(P = 0.01)，精米硒含量品种间变幅在18.3~121.6 μg/100g，平均为67.5 μg/100g，说明通过施用含量硒叶面肥，完全可以在低硒或缺硒区生产出高含量的富硒大米。其中：玉针香、泰丰优农39、泰优农39、巴斯马蒂等品种的硒富集能力较强，在富硒大米生产中具有较好的优势。3) 15个优质稻品种精米的有机硒占比(%)有极显著差异(P = 0.01)，品种间变幅在76%~100%，平均为92%，其中：香优52、悦两优美香新占、美香占2号、农香42的有机硒占比达到99%，说明这些品种的硒生物转化率高。4) 以精米总硒含量或富硒效率(权重系数0.4)和有机硒占比(权重系数0.6)为指标进行灰色关联度分析表明，15个优质稻品种中，富硒综合性状优良的常规稻品种是玉针香和巴斯马蒂、三系杂交稻品种是泰优农39、两系杂交稻品种是泰丰优农39和又香优浓39。玉针香和巴斯马蒂为1等优质米，香味浓郁，食味好；且垩白粒率和垩白度低，透明度1级，外观品质高；667 m²产量400 kg左右，又符合当地气候生产条件对品种的特征要求，在高档富硒香米开发中具有良好的推广应用前景。

关键词

优质稻，富硒，低镉，有机硒占比

Initial Study on Selenium and Cadmium Enrichment Characteristics of High Quality Rice Varieties

Min Kang¹, Jianhong Deng², Jianhua Gong^1*, Tao Huang¹, Shasha Mao¹, Zhenghong Zhao³, Yawen Deng¹

¹Zhuzhou Institute of Agricultural Sciences, Zhuzhou Hunan

²Zhuzhou Xiangzhiyou Agricultural Technology Development Co., Ltd., Zhuzhou Hunan

³Hunan Rice Research Institute, Changsha Hunan

Received: Jun. 8^th, 2022; accepted: Jul. 6^th, 2022; published: Jul. 19^th, 2022

ABSTRACT

In order to select high-quality rice varieties with high quality, low cadmium and rich selenium suitable for planting in Zhuzhou, 15 high-quality rice varieties were introduced to carry out a randomized block test. At the early stage of grain filling, biological nano selenium fertilizer was applied once, and the amount of selenium fertilizer (calculated as Se) was 1.25 g/667m². The contents of cadmium, total selenium and organic selenium in milled rice were determined. The results showed that: 1) under the condition of heavy cadmium pollution (1.17 mg/kg), there was a certain difference in the ability of cadmium absorption and accumulation among 15 high-quality rice varieties, but it did not reach a significant level (p = 0.05); Among them, the lowest cadmium content in milled rice is shiliangyou 1817, which indicates that this variety has the characteristics of low cadmium accumulation. Compared with no selenium fertilizer, the content of cadmium in milled rice with selenium fertilizer decreased by an average of 3.05%, and the effect of reducing cadmium was not obvious. 2) There was a very significant difference in grain selenium enrichment among 15 high-quality rice varieties (p = 0.01), and the selenium content in milled rice varied from 18.3 to 121.6 μG/100g, average 67.5 μG/100g, indicating that high selenium rich rice can be produced in low selenium or selenium deficient areas by applying selenium rich leaf fertilizer. Among them, Yu-zhen-xiang, Tai-feng-you-nong 39, Tai-you-nong 39, Basmati and other varieties have strong selenium enrichment ability, and have good advantages in the production of selenium rich rice. 3) The percentage of organic selenium in milled rice of 15 high-quality rice varieties was significantly different (p = 0.01), with a range of 76%~100% among varieties, with an average of 92%. Among them, the percentage of organic selenium in Xiang-you 52, Yue-liang-mei-mei-xiang-xin, Mei-xiang-zhan 2 and Nong-xiang 42 reached 99%, indicating that these varieties had a high selenium biological conversion rate. 4) Taking the total selenium content or selenium enrichment efficiency of milled rice (weight coefficient 0.4) and the proportion of organic selenium (weight coefficient 0.6) as indicators, the grey correlation analysis showed that among the 15 high-quality rice varieties, the conventional rice varieties with excellent selenium rich comprehensive characters were Yu-zhen-xiang and Basmati, the three-line hybrid rice variety was Tai-you-nong 39, and the two-line hybrid rice varieties were Tai-feng-you-nong 39 and You-xiang you-nong 39. Jade needle fragrance and basmati are first-class high-quality rice, with rich flavor and good taste; the chalky grain rate and chalkiness degree were low, the transparency was grade 1, and the appearance quality was high; the output of 667 m² is about 400 kg, which also meets the requirements of local climatic production conditions for varieties. It has a good prospect of popularization and application in the development of high-grade selenium rich fragrant rice.

Keywords:High Quality Rice, Selenium Enrichment, Low Cadmium, Proportion of Organic Selenium

Copyright © 2022 by author(s) and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

1. 引言

镉是联合国环境规划署确定的首位食品污染物，人类在食用镉超标大米后会造成呼吸道以及肝、肾脏功能的衰竭和病变，甚至引发癌症 [1] [2]。水稻是主要粮食作物，其吸收积累镉能力较强，但不同基因型之间存在显著差异 [3]。同时，稻米镉含量与土壤pH值呈显著的负相关 [4] [5] [6] [7]，与土壤有效态镉含量呈显著的正相关 [7]。为此，生产上采用了选育和推广镉低积累品种、降低土壤pH值和有效态镉含量等措施，以确保稻米镉含量达标。

硒是人和动物必需的微量营养元素，对人体具有抗癌作用 [8] [9]，缺硒可诱发的人体40余种疾病，如克山病、癌症等 [10] [11]。Se也是植物生长的必需微量元素，适宜浓度的硒可提高生物体的抗氧化酶活性、降低过氧化损伤，促进植物生长发育、提高抗逆境能力 [12] [13] [14]。硒还对镉、铅、汞等重金属元素具有拮抗作用，可缓解重金属对人体和植物的毒害 [15] [16]。我国是世界上缺硒最严重的国家之一，且缺硒多在以谷类为主食的居民中发生 [17] [18]，因而开发富硒大米成为硒营养健康的主要方向。研究表明，水稻生长期间施用硒肥，可降低大米中Cd、Pb、Cr等重金属的含量 [13]，显著提高大米中硒的含量 [19]，这为贫硒区开发富硒功能大米提供了理论基础。

株洲市属水稻主产区，是全国第一个过吨粮的地级市。株洲国民经济和社会发展第十四个五年规划和2035年远景目标纲要、以及近两年的株洲市1号文件要求：“实施优质粮食工程，推进种业和种植技术创新，扩大优质稻种植面积，发展功能性水稻。”但株洲土壤贫硒，发展富硒功能水稻需要通过外源补硒方式，株洲香之优农业科技发展有限责任公司与株洲市农业科学研究所合作研发的水稻富硒专利技术(ZL 2020 1 0289106.3)，不仅实现了中国科学院赵其国院士提出的从“富含”到“定向含有”的跨越，有效解决了外源补硒条件下富硒水稻生产中硒含量不稳定、有机硒占比不高的技术难题，而且可同时实现大米锌营养叠加，进一步提升了营养品质，为缺硒区高品质富硒大米生产提供了重要的科技支撑。但不同优质稻品种在外源补硒条件下的硒富集差异与硒品质差异鲜有报道。另一方面，株洲部分土壤存在轻中度镉污染，且优质稻生产中品种多、乱、杂现象突出。尽管不同水稻品种之间镉富集能力存在显著差异，但不同优质稻品种的差异很少有专项的研究报道。为此，开展优质稻品种硒镉富集特性研究，筛选出适宜株洲推广应用的优质、低镉、富集硒能力强且生物转化率高的优质稻品种，对引导企业、种粮大户发展高档优质稻，打造株洲优质稻米品牌，推进株洲粮食产业提质增效具有重要意义。

2. 材料与方法

2.1. 试验材料

水稻品种：参试品种共15个，其中：三系杂交优质稻品种2个，两系杂交优质稻品种8个，常规优质稻品种5个。本试验以1等优质米品种为主(详见表1)。

含硒叶面肥：“生物纳米硒肥”由武汉华硒生物科技有限公司提供，硒含量(以Se计) 5 g/L。

2.2. 试验设计

试验地点：株洲市农业科学研究所科研基地(芦淞区燎原村)。土壤镉含量1.17 mg/kg。

试验方法：每小区面积15.36 m²，株行距20 cm × 20 cm，小区间距50 cm。随机区组设计，三次重复。在镉积累比较中，设置独立对照(不施硒肥)，品种为改良型巴斯马蒂，独立对照不参与随机排列。于2021年6月18日播种，7月18日插秧。在灌浆初期施用生物纳米硒肥1次，硒肥用量(以Se计) 1.25 g/667m²，按每667 m²兑水30 kg人工喷雾。因参试品种的生育期不同，分别于9月9日、9月14日、9月18日喷施含硒叶面肥(详见表2)。肥水管理、病虫防治等同大田生产。

表1. 参试优质稻品种及米质情况表

表2. 含硒叶面肥施用时间表

2.3. 考察指标

品种镉富集能力：田间取稻谷样，检测精米镉含量。

品种硒富集能力：田间取稻谷样，检测精米总硒含量。

品种硒生物转化率：田间取稻谷样，检测精米有机硒含量并计算有机硒占比(%)。

2.4. 检测方法

土壤Cd含量按GB/T17141-1997方法检测，检测仪器为原子吸收分光光度计，设备型号为安捷伦240DUO。精米Cd含量按GB5009.15-2014方法检测；精米总Se含量按GB5009.93-2017第一法检测；精米有机Se含量按DBS42/002-2021方法检测；检测仪器均为原子荧光形态分析仪，设备型号为PF32。

2.5. 数据处理方法

试验数据采用Excel 2003进行数据整理和统计分析，采用Excel的TDIST函数进行随机区组试验的多重比较(LSD)。采用灰色关联度分析方法综合评价品种的富硒性状。

3. 结果与分析

3.1. 优质稻品种籽粒镉积累差异比较

由表3可知，参试品种中精米镉含量较高的品种有巴斯马蒂(P5)、荔香优丝苗(P15)、农香42 (P10)、悦两优美香新占(P8)、214A/制24 (P3)，较低的品种有拾两优1817 (P14)、明泰优丝苗(P12)，但不同品种之间精米镉含量的差异无统计学意义(P = 0.05)。精米镉含量最低的是拾两优1817 (P14)，说明该品种具有低镉积累特性，可能在轻中度污染区种植具有一定的优势。精米镉含量最高的是巴斯马蒂(P5)，说明该品种在生产中应采取严格的降镉措施。在土壤重度镉污染条件下(1.17 mg/kg)，与不施硒肥相比(独立对照)，施用硒肥品种的精米镉含量平均降低3.05%。

表3. 精米镉含量检测结果表(mg/kg)

3.2. 优质稻品种籽粒硒富集差异比较

由表4可知，在含硒叶面肥用量1.25 g/亩(以Se计)条件下，参试品种精米硒含量平均为67.5 μg/100g，达到高富硒水平，说明缺硒区通过施用硒肥可以实现富硒水稻生产。方差分析表明，不同优质稻品种精米总硒含量存在极显著差异(P = 0.01)，说明不同品种对硒元素的吸收、转运与籽粒积累能力明显不同。其中：籽粒富硒效率高的品种有玉针香(P7)、泰丰优农39 (P2)、泰优农39 (P11)、巴斯马蒂(P5)，较低的品种是香优52 (P1)、悦两优美香新占(P8)、美香占2号(P9)、农香42 (P10)、荔香优丝苗(P15)。进一步的多重比较表明，玉针香(A7)极显著高于香优52 (A1)、悦两优美香新占(A8)、农香42 (A10)、美香占2号(A9)、荔香优丝苗(A15)、明泰优丝苗(A12)、拾两优1817 (A14)、耘两优玖48 (A6)、野香优莉丝(A13)；泰优农39 (A11)和泰丰优农39 (A2)均极显著高于香优52 (A1)、悦两优美香新占(A8)、农香42 (A10)、美香占2号(A9)、荔香优丝苗(A15)、明泰优丝苗(A12)、拾两优1817 (A14)、耘两优玖48 (A6)；巴斯马蒂(A5)、214A/制24 (A3)和又香优浓39 (A4)均极显著高于香优52 (A1)、悦两优美香新占(A8)、农香42 (A10)、美香占2号(A9)、荔香优丝苗(A15)、明泰优丝苗(A12)；耘两优玖48 (A6)和拾两优1817 (A14)均极显著高于香优52 (A1)、悦两优美香新占(A8)、农香42 (A10)、美香占2号(A9)、荔香优丝苗(A15)。显然，在富硒优质稻生产中，应优先选择硒富集能力强且综合性状较好的品种，如玉针香、巴斯马蒂等。泰丰优农39、214A/制24和又香优浓39为两系优质稻苗头品种，具有较好的富硒品种开发潜力。

表4. 精米总硒含量检测结果与分析表(μg/100g)

^*富硒效率：是指每667 m²水稻补硒(以Se计) 1 g时，每100 g精米的硒含量(μg)。

3.3. 优质稻品种籽粒硒生物转化率差异比较

研究表明，在含硒叶面肥用量1.25 g/亩(以Se计)条件下，不同优质稻品种对硒元素的生物转化率存在较大差异(生物转化率是指植物通过根茎叶吸收硒元素并生物转化为硒蛋白、硒肽等有机形态的比例，水稻硒生物转化率 = (籽粒有机硒 + 茎有机硒 + 叶有机硒 + 根有机硒)/环境硒量)。对水稻籽粒而言，通常以有机硒含量占总硒含量的百分比即有机硒占比(%)来表示。由表5可知，施用含硒叶面肥的15个水稻品种的平均有机硒占比为92%，达到富有机硒标准(DBS42002-2021)，其中：香优52 (P1)、悦两优美香新占(P8)、美香占2号(P9)、农香42 (P10)的有机硒占比达到99%，说明其生物转化率很高，生产的富硒大米硒品质优良；有机硒占比低于90%的品种6个，其中最低的品种是耘两优玖48 (P6)、拾两优1817 (P14)和玉针香(P7)，这类品种在富硒水稻生产中应研究和推广提高硒生物转化率的有效措施，如株洲香之优农业科技发展有限责任公司与株洲市农业科学研究所合作研发的硒高生物转化技术，多年多点生产的巴斯马蒂富硒香米的有机硒占比达到99%以上，较本试验巴斯马蒂提高了十个百分点。由于有机硒占比是硒品质评价的关键指标，因此籽粒硒生物转化率高的优质稻品种在高档富硒米开发中具有较大优势。

表5. 精米有机硒占比(%)结果表

方差分析表明，精米有机硒占比(%)品种间存在极显著差异(P = 0.01)。进一步的多重比较表明，香优52 (P1)、悦两优美香新占(P8)、美香占2号(P9)、农香42 (P10)均极显著高于耘两优玖48 (P6)、拾两优1817 (P14)、玉针香(P7)、泰优农39 (P11)、泰丰优农39 (P2)、巴斯马蒂(P5)、214A/制24 (P3)；荔香优丝苗(P15)极显著高于耘两优玖48 (P6)、拾两优1817 (P14)、玉针香(P7)、泰优农39 (P11)、泰丰优农39 (P2)、巴斯马蒂(P5)；又香优浓39 (P4)、野香优莉丝(P13)、明泰优丝苗(P12)均极显著高于耘两优玖48 (P6)、拾两优1817 (P14)、玉针香(P7)。

4. 结论与讨论

4.1. 结论

1) 本研究表明，在土壤重度镉污染条件下(1.17 mg/kg)，不同优质稻品种的精米镉含量存在一定的差异，但未达到显著水平。其中，精米镉含量最低的是拾两优1817 (P14)，该品种米质1等，抗稻瘟病，在轻度镉污染区与稻瘟病区种植具有较大的优势。鉴于株洲农田土壤镉污染以轻中度污染为主，因此，不同优质稻品种在轻度或中度污染水平下的镉富集差异有待进一步研究。

2) 施用含硒叶面肥可降低精米镉含量。含硒叶面肥的降镉效应与施用的稻田土壤条件关系密切， 在低硒中低度镉污染的稻田中施用，能显著降低水稻精米中的镉含量；但在富硒高镉污染的稻田施用并无降镉效果 [20]。本试验因土壤镉含量高(1.17 mg/kg)，与不施硒肥的对照相比，施用硒肥品种的精米镉含量平均降低3.05%，降镉效果不明显，与黄太庆等 [20] 的研究结果一致。

3) 本研究表明，不同优质稻品种的硒富集能力和生物转化率均存在极显著差异，说明富硒优质稻品种选育具有可行性。但二者存在显著负相关(r = −0.8159, P = 0.05)，目前相关基因及位点尚不清楚，这意味着选育出硒富集效率高、有机硒占比也高的富硒优良品种存在较大的难度。根据富硒效率和有机硒占比两项指标，我们可以将参试品种分为三类，一类：硒富集效率高但有机硒占比较低。品种有玉针香(P7)泰丰优农39 (P2)、泰优农39 (P11)、巴斯马蒂(P5)、又香优浓39 (P4)和214A/制24 (P3)，这类品种生产富硒大米时的硒肥用量相对较少，但需采用提高硒生物转化率的技术措施，才能开发出高硒品质的富硒米产品。二类：硒富集效率低但有机硒占比高。品种有香优52 (P1)、悦两优美香新占(P8)、美香占2号(P9)、农香42 (P10)和荔香优丝苗(P15)，这类品种生产的富硒大米硒品质高，但较一类品种需提高硒肥用量。其中，香优52、农香42、美香占2号为常规优质稻品种，米质1级，是富硒品种选育的良好资源。三类：硒富集效率中等，且有机硒占比低或中等。品种有耘两优玖48 (P6)、拾两优1817 (P14)、明泰优丝苗(P12)、野香优莉丝(P13)，这类品种不宜做高品质富硒优质米开发。

4.2. 讨论

1) 品种富硒特性的综合评价方法。筛选富硒优良品种是富硒功能大米生产的重要基础，为综合评价品种的富硒特性，我们以精米总硒含量或富硒效率(权重系数0.4)、有机硒占比(权重系数0.6)为指标进行灰色关联度分析(表6)，结果表明：富硒性状最好的常规稻品种是玉针香和巴斯马蒂、三系杂交稻品种是泰优农39、两系杂交稻品种是泰丰优农39和又香优浓39。玉针香(排序第1)和巴斯马蒂(排序第4)为1等优质米，香味浓郁，食味好；且垩白粒率和垩白度低，透明度1级，外观品质高；亩产量400 kg左右，是高档富硒优质米开发的首选品种。泰优农39 (排序第6)为1等优质米，稻瘟病抗性较强(综合抗性指数3.8)，产量高(500 kg/667m²左右)，整精米率高(55.4%)，在高效富硒优质米开发中具有较好的前景。泰丰优农39 (排序第2)、又香优浓39 (排序第3)、214A/制24 (排序第5)为湖南省水稻研究所选育的两系杂交稻苗头品种，具有良好的硒富集特性。

表6. 优质稻品种富硒性状灰色关联度分析表

2) 富硒品种选育的目标性状。由于硒酸纳、亚硒酸纳等无机硒有较大的毒副作用，而有机硒无毒无害且吸收利用率高 [21]，因此，提高精米有机硒含量或籽粒硒生物转化率是今后富硒功能水稻品种选育的主要目标性状。另一方面，需要提高水稻品种对硒的吸收、转动能力，即富硒效率，这对富硒区提高大米硒含量而言十分重要。目前，富硒品种选育主要集中在提高富硒效率方面，而对硒生物转化率性状研究尚未得到足够的重视。

3) 富硒大米标准亟待规范。有机硒占比是富硒大米等农产品硒品质评价的关键指标。有机硒占比越高，说明硒品质越好。目前，只有湖北省制订了富有机硒农产品地方标准，有机硒占比≥80%即为富有机硒农产品。但从我们的实践来看，湖北省制订的富有机硒指标偏低，且没有出台国家标准。对于富硒大米的硒含量标准，地方标准与国家标准相差很大，管理比较混乱。在2008年《富硒稻谷》国家标准(GB/T22499-2008)中，规定“加工成符合GB1354规定的三级大米中的硒含量在0.04~0.30 mg/kg之间的稻谷为富硒稻谷”，但根据《食品安全国家标准 预包装食品营养标签通则》(GB28050-2011)，富硒稻米最低硒含量标准为≥0.15 mg/kg，且《食品安全国家标准 食品中污染物限量》(GB2762-2012)取消了对硒的限量规定，意味着硒已被认定为“非食品中的污染物”。但GB/T22499-2008至今并未进行修订，以至几乎所有的技术研究论文至今仍是以该标准为依据，客观上阻碍了水稻富硒技术的发展。湖南、陕西、湖北等省对富硒稻米的硒含量规定也差异很大，三省富硒大米的硒含量分别规定为0.150~1.000 mg/kg (T/HNFX001-2017)、>0.100 mg/kg (DB61/T556-2012)、0.20~0.500 mg/kg (DBS42/002-2021)。由此可见，正是由于富硒稻米标准制订的不统一，在很大程度上导致人们对富硒稻米和补硒的认识出现严重偏差，富硒稻米的技术研发、富硒产品开发者和消费者都无所适从。

4) 大力发展富有机硒优质大米的市场前景。中国是世界上缺硒最严重的国家之一，低硒地区占全国总面积的72%，部分地区属于严重缺硒区 [22]，约7亿人硒摄入量长期严重不足 [23]。因此，研究和应用外源补硒技术，提高农产品尤其是主粮的硒含量，对低硒或严重缺硒地区而言显得尤为重要。研究表明，植物可将环境中的无机硒(Se (IV)和Se (VI))转化易于被人体吸收的SeMet、SeCys等有机硒并参与植物蛋白的构成 [24]。有机硒是大米中硒形态的主要存在形式(约为78.4%~89.1%) [25] [26]，无论自然高硒土壤或通过外源施硒均可提高水稻籽粒中有机硒的含量及比例 [25] [27]，这为低硒或缺硒区富硒水稻开发奠定了理论基础。我国约三分之二的人口以大米为主食，因此大力发展高品质富有机硒优质大米具有广阔的市场需求。株洲香之优农业科技发展有限责任公司与株洲市农业科学研究所创新研发的稻米有机硒“定向含有”技术，可精准定制大米硒含量，以满足不同人群的补硒需要，如2019年开发的亚健康型富硒香米(品种为巴斯马蒂，下同)，硒含量定制为45 μg/100g，实际2个生产基地精米硒含量平均43 μg/100g (变幅为41 μg/100g~45 μg/100g)；2020年开发的健康型富硒香米，硒含量定制为30 μg/100g，实际3个生产基地精米硒含量平均31 μg/100g (变幅为30 μg/100g~32 μg/100g)；2021年开发的康复型富硒香米，硒含量定制为65 μg/100g，实际4个生产基地精米硒含量平均66.5 μg/100g (变幅为61 μg/100g~72 μg/100g)；且有机硒占比大幅提高到99.5%，大米硒品质居国内领先水平。该关键核心技术突破与应用，成功解决了富硒水稻生产中硒含量不稳定、有机硒占比不高的技术瓶颈，为高品质富有机硒优质米开发提供了重要的科技支撑，对实现全民科学补硒、定量补硒具有重要的现实意义。

基金项目

湖南省科技特派员服务乡村振兴专项(2021NK4137)，株洲市农业科学研究所科技攻关项目(2021-5)。

文章引用

康 敏,邓建红,龚建华,黄 涛,毛莎莎,赵正洪,邓雅文. 优质稻品种硒镉富集特性研究初报
Initial Study on Selenium and Cadmium Enrichment Characteristics of High Quality Rice Varieties[J]. 农业科学, 2022, 12(07): 551-561. https://doi.org/10.12677/HJAS.2022.127078

参考文献