基于HadCM3模式的气候变化下北京地区冬小麦耗水研究

水

文章编号：0559-9350（2010）09-1101-07利学报基于HadCM3模式的气候变化下北京地区冬小麦耗水研究

丛振涛，辛儒，姚本智，雷志栋

（清华大学水利水电工程系，水沙科学与水利水电工程国家重点实验室，北京100084）

摘要：气候变化在影响作物生长及农业生产的同时，也对作物耗水过程发生影响，进而影响能量平衡与水分循环、影响作物耗水及灌溉需水量。本文采用IPCC第四次报告中给出的未来气候变化情景中的COMMIT、SRA1B、SRA2和SRB1，选用大气环流模式HadCM3，利用作物模型CERES，模拟不同的灌溉处理下、不同的预测期内北京地区冬小麦的生长和耗水过程，分析气候变化对冬小麦耗水的影响。结果表明，对于未来气候变化情景下的北京地区冬小麦，生长期内的日平均气温明显升高，辐射下降；生长期缩短，生长期长度与气温存在良好的负相关关系；充分灌溉下产量与腾发量均有不同程度的下降；设定灌溉下产量与腾发量规律不明显；无灌溉下产量与腾发量以下降为主；灌溉需水量有所减少，平均降幅8.4%；产量与耗水量大体上呈正相关关系。

关键词：气候变化；耗水；冬小麦；HadCM3；CERES

中图分类号：S162.5+文献标识码：A

1研究背景

全球气候变化日益引起世界各国政府与民众的广泛关注。。气温增加、CO2浓度变化、降水［1-3］

量及频率变化都会对作物生长产生直接影响。气候变化在影响作物生长及农业生产的同时，也对作物耗水过程发生影响，进而影响能量平衡与水分循环、影响作物耗水及灌溉需水量。气候变化对农业生产影响的研究，除了定性分析，主要通过气候模式模拟未来情景下的气象条件，利用作物生长模拟模型定量模拟未来气象条件下的生长过程与产量。

需水的影响［11-12］［4-10］。评估气候变化对作物耗水可能影响的研究主要有两类：一类以Penman-Monteith等潜在腾发量的计算公式为基础，分析气候变化对作物；一类通过气候模式模拟未来情景下的气候条件，然后利用作物生长模拟模型定量模拟未来气象条件下的作物耗水与灌溉需水量。气候变化情景，可以通过分析历史资料或进行气象条件的假定（如气温升高1℃）得到，大多数的研究采用IPCC提出SERS（SpecialReportonEmissionSce narios）系列情景。气候变化的模拟，以大气环流模式（GCMs，GeneralCirculationModels）为代表。GC Generator）得到逐日气象资料，应用广泛的天气发生器有WGEN、CLIGEN和LARS-WG等［1］。在气候变化对作物耗水及灌溉需水量的影响研究方面，HadCM模式得到了广泛应用

型以及美国农业部开发的CERES系列模型，国内有RCSODSWheat［18］［16］［5，13-15］Ms的空间尺度较大，一般为几百公里。GCMs的时间尺度为月，需要利用天气发生器（WG，Weather。作物生长与耗水的模拟，主要采用作物模型，具有代表性的作物模型是荷兰的WOFOST系列模、WheatGrow模型［17］和ThuSPAC－

SAT-CERES等。已经在分析气候变化对作物耗水影响中得到应用的作物生长与作物耗水模型包括DS

［19］、CropWat4［20］和SWAP［21］等。已有的研究，较多关注灌溉需水量的变化，这些研究一

［16，21-22］般认为未来气候变化情景下，灌溉需水量将有所增加，增加的幅度为10%～30%。

本文采用IPCC第四次报告中给出的未来气候变化的8种情景中的4种，选用大气环流模式Had 来稿日期：2008-10-03

基金项目：国家自然科学基金项目（50979039；50939004）

作者简介：丛振涛（1973-），男，辽宁瓦房店人，副教授，主要从事水文水资源研究。E-mail：congzht@。

—1101—

CM3，利用作物模型CERES，模拟不同的灌溉处理方式、不同的预测期内北京地区冬小麦的生长和耗水过程，分析气候变化对冬小麦耗水的影响。

2

2.1研究方法IPCC情景与HadCM3模式在2007年IPCC第四次报告中，将未来气候变化情景归为8种，分别是：（1）1PTO2X：温室气体按1%速率增加至两倍后稳定；（2）1PTO4X：温室气体按1%速率增加至四倍后稳定；（3）20C3M：温室气体按20世纪增加速率增加；（4）COMMIT：温室气体维持在2000年水平；（5）PICTL：温室气体维持在工业化以前的水平；（6）SRA1B：经济高速增长，人口增长缓慢，新技术迅速得到应用；（7）SRA2：世界发展极其不均衡，自给自足和地方保护主义，人口持续增长，经济增长速度较慢；（8）SRB1：世界均衡发展，更加公平地在全球范围实现经济、社会和环境的可持续发展。这些情景是开展气候变化影响研究的基础。［1］

HadCM3在现状（1961—1990年）、2020s（2011—2030年）、2050s（2046—2065年）、2080s（2080—2099年）的逐月气象资料模拟结果，栅格大小为2.5°x3.75°（纬度×经度），共73×96个。HadCM3模型提供的原始数据为netcdf格式，可以用Matlab的MEXNC插件读取数据。

2.2试验站模型参数率定时采用的气象资料、土壤资料和作物生理资料均取自北京市水利科学研究所永乐店试验站结合国家自然科学重点基金项目“农田节水灌溉的增产与环境影响效应研究”进行的田间试验。试验站位于北京市通州区（地理坐标为116.8°E、39.7°N），土壤以粉壤土及砂壤土为主。试验区分为6行（编号A-F）、5列（编号1-5）共30个小区，每个大小为10m×5m，用混凝土板隔HadCM模型由英国气象局Hadley气候预测研究中心开发，目前版本为HadCM3。IPCC给出了开。小区南侧布置有小型气象站。试验开展时间为1998年10月—2001年9月，共播种3季冬小麦。试验测定的项目主要包括：气象资料、农艺资料、作物生理资料、土壤水热状况、土壤养分资料和作物含氮情况等。3年冬小麦品种均选用京冬8号（原代号为京农88-66）。3年的播种日期均为10月4日，收割日期分别为6月12日、6月12日和6月10日，播种量为300kg/hm。

根据HadCM3提供的数据，在栅格大小为2.5°×3.75°（纬度×经度）的栅格中选取北京永乐店所在位置的栅格：北纬（37.5°～40.0°）东经（116.25°～120.0°），见图1。以该栅格的数据作为未来气候变化情景下的气象数据。2

图1北京所在栅格位置永乐店位置

2.3CERES模型与参数率定CERES（CropEnvironmentREsourceSynthesis）由美国农业部开发，是最具代表性的作物模型之一。经不断改进，目前CERES-Wheat的最新版本为4.0。本研究采用美国农业部农业技术转移决策支持系统（DSSAT）v3.5中的CERES-Wheat模块。CERES-Wheat模型能够通过输入相对较少且容易获取的遗传参数和气象资料日值、土壤资料，模拟不同气候条件、土壤水分条件、氮肥条件、CO2浓度条件下冬小麦的生长过程，并可以计算出作物蒸发蒸腾量即水资源消耗量。—1102—

由计算春化作用、冬小麦冻害、分蘖和各生长阶段时期等模式组成；（2）生长子程序（GROSUB），包括计算叶片生长、叶片衰老、叶面积增加、光能量的截获、光合作用、植株不同器官生物量的分CERES-Wheat模型由主控程序和4个主要子程序组成，包括：（1）阶段发育子程序（PHENOL），配、茎杆、根和穗生长等模式；（3）土壤水分平衡子程序（WATBAL），有计算土壤径流、渗透、土壤蒸发、植物蒸腾、植物吸收的水分、水分胁迫和根分布等若干个模式；（4）氮素平衡子程序（NTRANS），包含计算土壤矿化速率、固氮、氮的硝化等模式。CERES-Wheat模型的输入文件包括试验管理文件FILEX、田间实测摘要文件FILEA、天气数据文件FILEW、土壤文件FILES、遗传参数文件FILEC。模型提供了一个GenCalc模块，可以进行遗传参数的率定。CERES-Wheat模型以积温作为冬小麦生育阶段的标志。CERES-Wheat模型的主要参数包括春化敏感性系数、光周期敏感性系数、灌浆持续期系数、每穗粒数系数、籽粒灌浆速率系数、单茎分蘖重系数、出叶积温。根据永乐店试验站作物生长过程实测资料，对CERES-Wheat模型中的遗传参数进行率定，各参数率定结果为2.0、4.9、-10.0、3.4、29.8、3.6和80.0℃，其中产量对籽粒灌浆速率系数G2最敏感。利用这些参数，模型对冬小麦产量、干物质积累和叶面积指数等均有较好的模拟效果

2.4日气象资料的生成［23］。

利用天气发生器，可以根据已有气候资料序列，计算序列的一些参数，生

成较为可靠的模拟序列来延长已有序列，天气发生器模拟的气候要素主要有降水量、最高气温、最低气温和太阳辐射等，本次研究利用北京永乐店站实测逐日气象资料，生成未来气候变化情景下的日气象资料，具体步骤为：（1）利用CERES模型中的WGEN模块，根据北京永乐店站1961—1990年的逐日气象资料，生成2001—2100年的日气象资料，包括降水量、最高气温、最低气温、太阳辐射；（2）以HadCM3给出的情景20C3M下1961—1990年的逐月气象资料为基础，分别计算各种气候变化情景在2020s、2050s、2080s的逐月气象资料与情景20C3M在1961—1990年的逐月气象资料的差值，将得到的差值按月叠加到步骤（1）得到的逐日气象资料，得到未来气候变化情景下的逐日气象资料，作为CERES模型的驱动。

2.52020s（2011—2030年）、2050s（2046—2065年）、2080s（2080—2099年）的逐日气象资料进行模拟，其中COMMIT情景无2080s资料。现状（1979—1998年）采用实测资料。

土层厚度设为90cm，土壤参数的选择见表1。

3模拟方案设定根据HadCM3的数据，选取COMMIT、SRA1B、SRA2、SRB1四种情景，利用

冬小麦生长的影响，假定养分充足即冬小麦生长不受养分胁迫。

表1

名称

土壤颜色

裸地反射率

根生长因子

排水速率参数值0.130.5033土壤水分的初始条件，根据2000年实测资料，取为0.343cm/cm。本研究重点考虑水分条件对3CERES模型中的土壤参数设定参数值0.25cm3/cm30.35cm3/cm3名称矿化因子干容重黏粒含量砾含量参数值1.40g/cm350.0%名称饱和含水率有机碳含量淤泥含量总氮参数值1.12%0.0%名称光合因子最适土壤含水率下限最适土壤含水率上限饱和含水率低于0.25cm/cm开始灌溉，高于0.35cm/cm停止灌溉；设定灌溉的灌溉制度为冬小麦生育期内灌水433在灌溉处理方式上，分别考虑充分灌溉，设定灌溉和不灌溉。充分灌溉指60cm平均土层含水率

32次，每次60mm（600m/hm），灌溉方式为喷灌，灌溉日期分别为11月23日（冬前）、3月29日（返青

期）、4月15日（拔节期）和5月9日（灌浆期）。这样，共有36种模拟方案，详见表2。

3

3.1结果与讨论模拟结果根据以上方案设定与研究方法，基于永乐店试验站的率定参数，利用CERES-Wheat

—1103—

模型，HadCM3模式下各种气候变化情景的北京地区冬小麦生长期长度、产量、腾发量与灌溉需水量模拟结果见表2。

表2

方案

2345678

10111213141516171819202122232425262728293031323334359气候变化情景20C3MCOMMITCOMMITCOMMITCOMMITCOMMITCOMMITSRA1BSRA1BSRA1BSRA1BSRA1BSRA1BSRA1BSRA1BSRA1BSRA2SRA2SRA2SRA2SRA2SRA2SRA2SRA2SRA2SRB1SRB1SRB1SRB1SRB1SRB1SRB1SRB120C3M北京地区冬小麦HadCM3模式下CERES模型模拟结果代表年现状现状现状2020s2050s2020s2050s2020s2050s2020s2050s2080s2020s2050s2080s2020s2050s2080s2020s2050s2080s2020s2050s2080s2020s2050s2080s2020s2050s2080s2020s2050s2080s2020s2050s灌溉处理充分灌溉设定灌溉无灌溉充分灌溉充分灌溉设定灌溉设定灌溉无灌溉无灌溉充分灌溉充分灌溉充分灌溉设定灌溉设定灌溉设定灌溉无灌溉无灌溉无灌溉充分灌溉充分灌溉充分灌溉设定灌溉设定灌溉设定灌溉无灌溉无灌溉无灌溉充分灌溉充分灌溉充分灌溉设定灌溉设定灌溉设定灌溉无灌溉无灌溉

无灌溉生长期长度/d251.9251.9251.9252.6251.9252.6251.9252.6248.6237.4234.3248.6237.4234.3248.6237.4234.3246.2239.3233.4246.2239.3233.4246.2239.3233.4247.7241.4239.2247.7241.4239.2247.7241.4产量（/kg/hm2）447910676515688843144366483515腾发量/mm401.8196.5531.1537.4398.6387.7188.9185.7520.1538.1525.1403.0381.1409.6202.7184.9226.6538.3519.1505.0380.9386.4401.5166.1195.1223.7518.7533.7524.3411.0396.0385.5211.1213.3383.4393.3405.7240.0240.0240.0442.4398.4356.8240.0240.0240.0灌溉需水量/mm240.0411.5423.5240.0240.0657668336258490550005411822482391.3422.2371.8240.0240.0240.01385666467216085422049355301261516105067136767629956115083487510939403.2气候要素分析根据HadCM3提供的未来气候变化情景下的气象数据，得到不同情景下的北京地区冬小麦生育期平均气温、平均辐射与降水量。除了COMMIT情景，未来气候变化情景下，北京地区冬小麦生长期内气温不同程度升高（见图2），辐射不同程度下降（见图3），降水量变化无明显规律。SRB1情景下的气温增幅较小，辐射降幅也较小。

3.3生长期由于CERES模型对生长期的模拟主要考虑光（辐射）和热（气温）的影响，没有考虑土壤水分、土壤肥分、大气中CO2浓度等对生育时间的影响，因此，不同灌溉处理对作物生长期没有影响。与现状相比，除了COMMIT情景，各种气候变化情景下北京地区冬小麦生长期存在不同程度地—1104—

图2不同气候变化情景的生育期平均气温2生育期辐射/MJ（/m·d）3不同气候变化情景的生育期平均辐射

缩短，见表2。冬小麦生长期长度主要受气温的控制，两者相关性较好（见图4），气温每升高1℃，生长期缩短约6d。

3.4产量与耗水量（1）充分灌溉条件下的产量与腾发量。充分灌溉条件下，各种气候变化情景的北京地区冬小麦产量与腾发量均有不同程度的下降，见表1。产量平均降幅5.1%，最高降幅16.4%；腾发量平均降幅2.9%，最高降幅7.3%。气温的升高有利于增加逐日的光合作用，但同时引起生长期的缩短，辐射的下降将引起光合作用的下降，三方面共同作用的结果是产量下降。各种情景下，产量在2050s高于2020s和2080s，原因可能是气温增加显著而辐射下降不明显（见图2、图3）。在充分灌溉时，日腾发量主要受辐射影响，辐射下降将引起腾发量下降，同时生长期缩短也将导致腾发量下降；（2）设定灌溉条件下的产量与腾发量。设定灌溉条件下，水分供给不足，北京地区冬小麦产量与腾发量的模拟结果均小于充分供水条件下的产量与腾发量。尽管同样有气温升高与辐射下降的影响，设定灌溉条件下，部分气候变化情景的冬小麦产量有所增加（见表2），这可能是由于降雨过程与灌溉配合更好的原因，水分条件起着决定性作用。腾发量与产量的规律类似；（3）无灌溉条件下的产量与腾发量。无灌溉条件下的产量只有充分灌溉的1/8左右，各种气候变化情景的北京地区冬小麦产量与腾发量以下降为主，个别情景增加，见表2，主要受降水过程控制。

3.5需水分析根据充分灌溉的模拟结果，可以得到各情景下的灌溉需水量。北京地区冬小麦未来情景的灌溉需水量模拟结果见表2，灌溉需水量平均降幅8.4%，最高降幅18.1%。实际上，灌溉需水量还受到降雨过程的影响，本文利用天气发生器得到未来气候情景下的降雨过程可能无法反映降雨过程的变化，如果降雨更加不均匀，将引起灌溉需水量的增加。

充分灌溉、设定灌溉和无灌溉条件下，产量与降水量的相关系数分别为0.59、0.72和0.90；3种条件下，腾发量与降水量的相关系数分别为0.73、0.86和0.90。总体而言，灌溉的减少，使作物生长更多地依赖于降水。另外，产量与耗水量大体上呈正相关关系（见图5），一般而言，同时实现节水与增产的困难较大。

图4生长期长度与生育期平均气温关系图5产量与耗水量关系

—1105—

4结论

基于不同气候变化情景（COMMIT、SRA1B、SRA2、SRB1）的大气环流模式（HadCM3）的输出，可以作为作物生长模拟模型（CERES）的驱动，从而模拟各气候变化情景下不同预见期的冬小麦产量与耗水，形成以下认识：（1）未来气候变化情景下（除了COMMIT），北京地区冬小麦生长期内的日平均气温明显升高，辐射下降，这是引起冬小麦产量与耗水量变化的主要气候变化要素；（2）未来气候变化情景下（除了COMMIT），各种模拟方案的北京地区冬小麦生长期长度缩短；冬小麦生长期长度与气温存在良好的负相关关系，气温每升高1℃，生长期缩短约6d；（3）充分灌溉条件下，受气温升高、辐射下降的综合作用，各种气候变化情景的北京地区冬小麦产量与腾发量均有不同程度的下降，产量平均降幅5.1%，腾发量平均降幅2.9%；设定灌溉条件下，受气温、辐射、降水的综合作用，气候变化情景的北京地区冬小麦产量与腾发量规律不明显；无灌溉条件下主要受降水过程控制，各种气候变化情景的北京地区冬小麦产量与腾发量以下降为主，个别情景增加；（4）不考虑未来情景下降水过程的可能变化，充分灌溉条件时未来气候变化情景下北京地区冬小麦的灌溉需水量有所减少，平均降幅8.4%，最高降幅18.1%。无灌溉时，冬小麦产量及耗水与降水相关性大。总体上，产量与耗水量呈正相关关系；（5）气候变化情景下的日气象资料生成，未考虑降水过程变化的影响，降水过程的愈加不均匀将引起灌溉需水量的增加；作物生长模拟中，未考虑CO2浓度、土壤养分未考虑品种变化带来的影响；这些方面需要开展进一步的研究工作。的影响，这些也是影响作物产量与耗水的重要因素；此外，未来情景下的冬小麦生长与耗水模拟也

参考文献：

［1］IPCC.ClimateChange：ThephysicalscienceBasis.contributionofworkinggroupitothefourthassessment［M］.［2］ThomsonAM，RosenbergNJ，IzaurraldeRC，etal.ClimatechangeimpactsfortheconterminousUSA：Aninte

69（1）：89-105.NewYork：CambridgeUniversityPress，2007.gratedassessment-Part5.Irrigatedagricultureandnationalgraincropproduction［J］.ClimaticChange，2005，［3］《气候变化国家评估报告》编写委员会.气候变化国家评估报告［M］.北京：科学出版社，2007.

［4］PapajorgjiP，JonesJW，PeartRM，ingcropmodelsandgeographicinformation-systemstostudythe

1994，53：82-86.impactofclimate-changeinthesoutheasternUSA［J］.SoilandCropScienceSocietyofFloridaProceedings，changeprojectionsonagriculturalproductivityandirrigationwatersupplyintheconterminousUnitedStates-I.turalandForestMeteorology，2003，117（1-2）：73-96.［5］RosenbergNJ，BrownRA，IzaurraldeRC，etal.IntegratedassessmentofHadleyCentre（HadCM2）climateClimatechangescenariosandimpactsonirrigationwatersupplysimulatedwiththeHUMUSmodel［J］.Agricul ［6］PeirisDR，CrawfordJW，GrashoffC，etal.Asimulationstudyofcropgrowthanddevelopmentunderclimate［7］金之庆，葛道阔，石春林，等.东北平原适应全球气候变化的若干粮食生产对策的模拟研究［J］.作物学［8］ErdaL，WeiX，HuiJ，etal.ClimatechangeimpactsoncropyieldandqualitywithCO2fertilizationinChina［J］.PhilosTransRSocLondBBiolSci，2005，360（1463）：2149-2154.

究，2007（5）：208-213.

［10］XiaoGJ，ZhangQ，YaoYB，etal.Effectsoftemperatureincreaseonwateruseandcropyieldsinapea-spring［11］刘晓英，林而达.气候变化对华北地区主要作物需水量的影响［J］.水利学报，2004（2）：77-82.wheat-potatorotation［J］.AgriculturalWaterManagement，2007，91（1-3）：86-91.［9］张建平，赵艳霞，王春乙，等.未来气候变化情景下我国主要粮食作物产量变化模拟［J］.干旱地区农业研报，2002（1）：24-31.change［J］.AgriculturalandForestMeteorology，1996，79（4）：271-287.

—1106—

［12］佟玲，康绍忠，粟晓玲.石羊河流域气候变化对参考作物蒸发蒸腾量的影响［J］.农业工程学报，2004（2）：［13］IzaurraldeRC，RosenbergNJ，BrownRA，etal.IntegratedassessmentofHadleyCenter（HadCM2）cli mate-changeimpactsonagriculturalproductivityandirrigationwatersupplyintheconterminousUnitedStates-PartII.Regionalagriculturalproductionin2030and2095［J］.AgriculturalandForestMeteorology，2003，117

rigationwaterdemand［J］.ClimaticChange，2007，84（3-4）：441-461.

culturalEngineering，2004：13-29.15-18.（1-2）：97-122.［14］ElgaaliE，GarciaLA，OjimaDS，etal.Highresolutionmodelingoftheregionalimpactsofclimatechangeonir ［15］DeWD.Climatechangeandirrigationrequirementsinthemediterraneanenvironment［J］.ActualTasksonAgri ［16］高亮之，金之庆，黄耀，陈华.作物模拟与栽培优化原理的结合－RCSODS［J］.作物杂志，1994（3）.［17］曹卫星，罗卫红.作物系统模拟及智能管理［M］.北京：华文出版社，2000.

［18］丛振涛，倪广恒，雷志栋.用于田间作物—水分关系研究的ThuSPAC模型［J］.沈阳农业大学学报，2004，［19］PeartRM，CurryRB，RosenzweigC，etal.EnergyandirrigationinsoutheasternUSagricultureunderclimate［20］JorgeJ，FerreresE.Irrigationscenariovsclimatechangescenario［M］.DetectingandModellingRegionalCli ［21］YanoT，AydinM，HaraguchiT，etal.ImpactofclimatechangeonirrigationdemandandcropgrowthinaMedi ［22］DiazJA，WeatherheadEK，KnoxJW，etal.Climatechangeimpactsonirrigationwaterrequirementsinthe［23］丛振涛，王舒展，倪广恒.气候变化对冬小麦潜在产量影响的模型模拟分析［J］.清华大学学报（自然科学

版），2008，48（9）：46-50.GuadalquivirriverbasininSpain［J］.RegionalEnvironmentalChange，2007，7：149-159.terraneanenvironmentofTurkey［J］.Sensors，2007，7：2297-2315.mateChange，2001.change［J］.JournalofBiogeography，1995，22（4-5）：635-642.35（5-6）：459-461.

ImpactofclimatechangeonwateruseofwinterwheatwithHadCM3Model

CONGZhen-tao，XINRu，YAOBen-zhi，LEIZhi-dong（TsinghuaUniversity，Beijing100084，China）

Abstract：Climatechangeaffectscropgrowththatresultsinvariationincropwateruseandconsequentlyirrigationrequirement.ThispaperinvestigatedtheimpactofthepredictedclimatewithHadCM3modelus ingtheIPCCscenariosofCOMMIT，SRA1B，SRA2andSRB1onwinterwheatgrowthandwaterusenearBeijingunderdifferentirrigationtreatments.GrowthofwinterwheatissimulatedusingCERESmodel.Theconclusionincludes（1）theairtemperaturewillincreaseandtheradiationwilldecreaseduringthegrowthseason；（2）thegrowthdurationwillbeshortenedanditcorrelatesnegativelywithairtemperature；（3）un dersufficientornoirrigationcondition，yieldandtheevapotranspirationwilldecreaseunderdifferentscenar ios；（4）irrigationdemandwillreduceby8.4%undersufficientirrigation；（5）wheatyieldcorrelatespositive lywithevapotranspiration.Keywords：climatechange；wateruse；winterwheat；HadCM3；CERES

（责任编辑：吕斌秀）

—1107—

本文来源：https://www.bwwdw.com/article/ojx4.html