基于目标库容曲线的水库出流模拟模型.pdf

年 月水利学报 第 卷第 期文章编号：（）收稿日期：；网络首发日期：网络首发地址：?基金项目：国家自然科学基金项目（）；长江水科学研究联合基金项目（）作者简介：舒鹏（），博士生，主要从事水文模型研究。：基于目标库容曲线的水库出流模拟模型舒鹏，熊立华，陈杰，尹家波，程磊，周研来（武汉大学 水资源与水电工程科学国家重点实验室，湖北 武汉 ）摘要：随着水库的不断建设，水库对水循环的影响越发重要，水库出流过程已成为陆地水循环的重要一环。为了模拟水库对水循环的影响，需要构建水库出流模拟模型（，）模拟水库的出流过程，然而，目前提出的 仅适用于预设运行模式的水库，通用性不强。本文假设水库在年内不同季节的目标库容构成目标库容函数曲线（，），通过为不同水库设置不同的 提高 的通用性，并从两方面考虑 对水库出流的影响，提出了两种模型结构：（）水库在目标库容较高的季节倾向于蓄水，出流减少，在目标库容较低的季节则倾向于放水，出流增加；（）水库出流应使得实际蓄水量接近目标库容。根据两种模型结构，构建了 种基于 的 ，并将其与常用的 应用于长江、珠江、黄河流域的 座水库中。结果发现：（）考虑了两种模型结构的 的出流模拟、蓄水量模拟表现均为最优，说明所提出的两种模型结构合理有效；（）考虑了 的 模型在多年调节、年调节水库中的出流模拟精度、蓄水量模拟精度均远超其他 ，说明考虑 可以有效提高在这两类水库的模拟精度，提高了 的通用性。关键词：水库出流；模拟模型；水库蓄水量；水库目标库容曲线中图分类号：文献标识码：?研究背景在水库众多的大型流域，水库出流过程已经成为陆地水文循环的重要一环 。在流域的水文模拟过程中考虑水库的影响，不仅可以提高水文模拟精度，对研究气候变化、人类活动条件下水循环变化过程也有着重要的意义 。为了考虑水库调度过程的影响，有学者在集总式模型中增加相关参数以代表流域内部所有水库的调蓄作用 ，这种处理方式本质上是将水库的调蓄作用等效为流域下垫面性质的改变 ，在实际运用过程中具有较大限制，流域中有新的水库建成时则需要重新率定模型。同时，这种做法也难以探究单座水库对下游流量变化的影响。为了更好地模拟水库对水循环过程的影响，研究者往往需要概化水库运行模式，通过构建水库出流模拟模型（，），直接模拟流域内各水库的出流过程 。在构建 时，如何合理概化水库的运行模式，即如何确定水库在年内甚至年际间蓄水或放水的时机以及强度是最为关键的问题。归纳了常见的概化水库运行模式的方法，可以用到其中的一种或几种，这些方法包括多功能分区法、单一分区法、目标库容法以及条件曲线法。多功能分区法是将水库按功能分区，当水库水位位于不同分区时，采用不同的公式计算水库出流，如 、等所构建的 。单一分区法则是不进行分区，不考虑水库的不同功能，直接从年、甚至多年尺度上考虑水库库容的分配方式，如 所构建的 （，）假定水库在汛期蓄水而在非汛期均匀将汛期蓄水量放出，对于多年调节水库，还需要根据前两年流域的来水情况进一步调整放水强度。目标库容法则是基于目前水库的实际库容与目标库容，来计算水库的出流量。目标库容是水库调度计划中设定的，指水库在不同季节所需维持的特定库容，如 、均在模型中假设水库在汛期的目标库容较小而非汛期则回到较高水平，以考虑水库的防洪功能，考虑到我国水库普遍存在汛限水位 ，目标库容法对我国水库的出流模拟有着较好的适用性。条件曲线法类似于调度图法，如雷晓辉等 设计了一套归一化的水库通用防洪调度图以及智能优化技术，采用洪水调度图来模拟水库的出流过程。目前构建 遇到的难题是，传统 所概化的水库运行模式不具备通用性。出现这一问题的原因是大部分 在构建之初便对水库的运行模式进行了预设，使得 难以适用于具有非预设运行模式的水库。以 的 （）模型以及 为例，两个模型所概化的水库运行模式具有较大不同，适用于不同的水库。模型假设水库在汛期目标库容较低而非汛期较高，使得水库在汛期蓄水量较低，而在汛期结束后逐渐蓄水到最高水位。而 模型则假设水库在汛期蓄水，非汛期均匀将汛期蓄水放出，蓄水量在汛期逐渐增大，而在非汛期逐渐减少。为了提高 的通用性，本文基于目标库容法提供了一种新的构建 的思路。假设水库在年内不同季节的目标库容构成目标库容曲线 ，在某种程度上描述了水库的运行模式：在目标库容较大的季节，水库倾向于蓄水，出流减小；而在目标库容较小的季节，水库倾向于放水，出流增大。通过调整 的曲线形状，可以使 适用于具有不同运行模式的水库。基于此，本文的研究目标为构建基于 的 （，），通过为不同的水库设置相应 ，以提高模型通用性。将所提出的 应用于我国 座水库，以验证 是否可以有效应用于具有不同运行模式的水库的出流模拟。研究方法 基于目标库容曲线的水库出流模拟模型 目标库容曲线目标库容 是水库调度计划中设定的，水库在不同季节所需维持的特定库容。水库的兴利库容 和死库容 是水库的最为重要的特征库容之一，界定了目标库容的变化范围。如图 （）所示，在水库的正常使用过程中，水位一般尽可能保持在死水位和正常高水位之间，因此，水库的目标库容应满足 （）这一条件。图 目标库容曲线与 模型水库出流过程水库的目标库容曲线 描述了水库在年内的运行模式，如图 （）所示：在 较大的季节，水库倾向于蓄水；而在 较小的季节，水库倾向于放水。假设水库 的波动变化过程服从一个函数 （），称该函数为 函数。基于水库的特征库容、，根据 函数，可以得到水库在各个时刻的 ：，（），（）（）式中 为 时刻对应的日期在年内的日序数（，），不考虑闰年的情况，取值为 。水库作为调节径流时空分布不均的工程措施，总是存在蓄水期和放水期。考虑到 函数与 函数具有严格的周期性且仅有一个波峰和一个波谷，可以与水库每年的蓄水期、放水期相对应，因此，本文提出了一个以 函数为基础的通用 函数：（）（）（）式中：为 函数的最小值，范围为 ；为水库目标库容达到最大值的日序数，范围为 。如图 （）所示，函数的特点是当 时，（）达到最大值 ；若 ，则当 时，函数达到最小值 ，若 ，则当 时，函数达到最小值 。水库出流计算 模型的基本结构如图 （）所示，将水库出流量 分为得到发挥经济效益的有效出流量 以及未发挥经济效益的弃水 两部分。，（），主要受水库功能的影响，是水库为了满足灌溉、发电等功能而下泄的流量。参考 在 年提出的 （，），基于水库水的滞留时间（，）和一个根据水库实际蓄水量计算得到的水库库容校正系数进行 ，的计算。，（）式中：为水库水在 时刻的滞留时间；为水库在 时刻的实际蓄水量；为水库在 时刻的库容校正系数。假设 ，从两方面对 ，产生影响：（），较高时水库处于倾向于蓄水，较低，反之水库倾向于放水，较高；（）应尽可能向 ，靠近，当 ，时，减小 ，以蓄水，反之则增大 ，以放水。首先考虑第一点，表征了水库的放水速度，根据式（），在 与 保持不变的条件下，越大，越小，反之则 ，越大。因此，在 的基础上做出改进，假设 在年内并非保持不变，而是随着 ，的变化而变化：当水库的 ，较大时，水库在该季节倾向于蓄水，水库水体更新速度变慢，增加；而当水库的 ，较小时，水库在该季节倾向于放水，水库水体更新速度加快，下降。如式（）所示，函数 （）指示了 ，的波动变化过程，与之成正比。因此，基于 （）来计算 ：（）（）式中 、为模型参数，需要率定，参数范围为 。为了使 应尽可能向 ，靠近，在 的基础上对 进行重构，使其根据 与 ，的大小关系对 ，进行调整，当 ，时，减小 ，以蓄水，当 ，时则增大 ，以放水。所采用库容校正系数记为 ，水库的实际蓄水量越大，越大，的计算如式（）所示。根据 ，重构的库容校正系数记为 ，的计算如式（）所示。()（），()（）式中：为水库的最大库容；为模型参数，需要率定，参数范围为 。最后计算弃水 ，产生 ，的原因主要有两点，由于水库调度能力有限产生的弃水以及为了维持下游河道生态环境而下泄的生态流量。，在一定程度上与水库入流 ，存在正相关关系：在汛期，越大，因水库防洪库容与下泄能力有限而产生的弃水量越大；而在非汛期，虽然水库在此期间产生的弃水较少，但维持河流生态系统动态平衡所需要的生态流量也与 ，在一定程度上呈正比关系：在河道天然流量较大的季节，所需要的生态流量也相对较大。因此，假设 ，与 ，呈正比，即有固定比例的 ，化为弃水 ，：，（）式中 为比例系数，参数范围为 ，）。得到 ，后，基于水库水量平衡方程计算水库下一时刻的实际库容 。记时间步长为 ，忽略库区的降雨、蒸散发以及渗漏过程，则 为：（，）（）对比方案设置为了探究本文所提出的 的有效性，构建两类 以作对比：（）不同结构的 模型；（）常用的 模型。本文所提出的 模型引入了目标库容曲线 来描述水库年内 的变化过程，并从两方面考虑 对 的影响，构建了两种模型结构：（）在年内的相对大小描述了水库运行模式，水库水滞留时间 随 的变动而变化，具体见式（）；（）水库的出流过程应使得实际蓄水量接近目标库容，主要通过 来实现，具体见式（）。从模型机理上来看，模型结构（）体现了水库调度计划对水库出流的影响，影响程度仅和 有关，与外界条件无关。模型结构（）则体现了水库的实际蓄水量对水库出流的影响，其影响程度不仅与 有关，还与前期来水条件有关。基于这两种模型结构，构建了 组不同的 模型（表 ），以探究引入 后，是否能有效提高模型的模拟精度。模型未考虑 对水库出流的影响，、分别考虑了一种模型结构，则是结构最为全面的模型。另外，虽然没有考虑 对 的影响，与 模型相比，额外考虑了水库的弃水。表 对比模型方案 名称滞留时间 库容校正系数 （）（）除了不同内部结构的 间的对比以外，引入了常用的 以作模型间的对比（表 ）。其中，方案为不考虑水库调度过程的基准方案，直接采用水库的入流量模拟出流量，用以验证引入 后，是否能够提高水库出流的模拟精度 。表 其它水库出流模拟模型 参考文献模型结构模型参数 ()珚 珚()珚 珚 ()珚 ()，珚 ，为水库的最大库容、死库容；珚 为水库多年平均入流量，?；，为无量纲模型参数，需要率定，参数范围均为 ；为每年开始时水库的蓄水量；，为无量纲模型参数，需要率定，的范围为 ，的范围为 。数据与模型率定 水库数据将 应用于我国 座水库，以验证模型效果，其中有 座水库在长江流域，座水库在珠江流域，座水库在黄河流域，按蓄水系数（，）来划分水库的调节能力，的跨度为 ，包含季调节水库、年调节水库以及多年调节水库。并根据 将水库划分为三种类型：季调节水库（），年调节水库（），多年调节水库（）。长江流域、珠江流域水库的数据分别来自于长江水利委员会、珠江水利委员会，包括水库日出流数据以及根据出流、库容所反演的日入流数据。黄河流域水库的数据来自于由黄河水利委员会公布的水库日出流数据，以及水库上游作为水库入流的控制水文站点的日径流数据作为水库的入流数据（：?：?）。另外，座水库中，其中 座水库还有对应的实测蓄水量数据。各水库数据的时间跨度不一，少则 年，多则 年。本文 座水库具有好的代表性，其建设规模、功能、气候条件等均不一致。如表 所示，座水库中，最小的是总库容为 亿 的鲁布革水库，最大为 亿 的三峡水库。水库功能各异，包含发电、防洪、灌溉、供水、旅游等多种功能。水库位置纬度跨度接近 ，年平均降水量变幅为 。表 水库信息表水库名称 流域年降雨量?总库容?亿 蓄水系数 起始时间结束时间水库名称 流域年降雨量?总库容?亿 蓄水系数 起始时间结束时间三峡长江?岩滩珠江?水布垭长江?天一珠江?隔河岩长江?云鹏珠江?锦屏一级 长江?鲁布革珠江?丹江口长江?龙滩珠江?江口长江?长洲珠江?彭水长江?光照珠江?二滩长江?西津珠江?宝珠寺长江?红花珠江?碧口长江?龙羊峡黄河?草街长江?刘家峡黄河?亭子口长江?万家寨黄河?百色珠江?小浪底黄河?模型率定及评估将数据划分为两个时期，前三分之二为率定期，后三分之一为验证期。率定算法则采用 （）。优化算法 采用下山单纯形法进行最优解的搜索，为了避免算法受局部极值的影响，同时采用多个单纯形进行并行搜索，并实现了多个单纯形之间的信息共享，可以同时优化水文模型的多个参数，已在水文领域得到了广泛应用。以水库出流的纳什效率系数（，）为目标函数进行参数率定，以 模型为例，记 模型的 为 ，其计算公式为：，()，()（）式中：，、，分别为水库在 时刻的实测出流量与 模型的模拟出流量；，为水库实测出流量的均值。除了水库出流量以外，也能模拟水库的蓄水量。考虑到 的计算公式中涉及到误差的平方，使得峰值对 的影响较谷值更显著，而 所提出的 （）从均值、标准差、相关系数三个方面评价模拟精度。因此，采用蓄水量的 值 来评价各个 的蓄水量模拟精度，以 为例，的模拟蓄水量为 ，实测蓄水量为 ，基于 与 计算得到的 记为 ，的计算公式为：()()()槡（）式中：、与 、分别为 与 的均值、标准差；为 与 的线性相关系数。结果 水库出流模拟图 展示了各 验证期 的均值以及分布情况。从所有水库 的均值 ，来看，模型表现最差，其 ，仅为 ，与 持平。、则次之，在 左右。考虑了 的 的出流模拟精度显著高于其它模型，在 以上。另外从各 的箱体图来看，也是 的表现最优，分位线在 以上，远超其他模型，说明考虑 可以有效提高不同水库的出流模拟精度，模型适用性更高，在大部分水库中均能取得较好的出流模拟结果。图 不同 验证期出流 箱体图将各模型所模拟的各水库在验证期的 绘制成热力图，按 值对水库进行降序排列（图 ）。作为不考虑水库调度过程的基准模型，直接采用入流模拟出流，反映了水库对入流的改变程度，因此可以看到，多年调节水库的 显著低于季调节水库。对比不同 在各类水库的表现，可以发现，在季调节水库中，除了 相对较大的二滩、锦屏一级、云鹏三座水库以外，几乎所有 都难以提高其他季调节水库的出流模拟精度，主要是因为季调节水库的 较高，如草街、红花等水库的 已经达到 以上，提升空间有限。在年调节以及多年调节水库中，尤其是在多年调节水库中，的 显著高于其他模型，出流模拟精度更高。但是，与年调节水库相比，即使是表现最优的 模型，多年调节水库的 仍相对较低，造成这一现象的原因可能有以下三点。首先，多年调节水库对径流的改变程度大，水库行为复杂，且气候条件在年际间的变化规律相比于年内不确定性更强，模拟难度更大。其次，由于多年调节水库库容大，出流相对平缓，在其他指标相差无几时，相对平缓的流量序列计算得到的 值可能相对较低，实际上，若采用 评估出流模拟精度，多年调节水库和年调节水库的差异则会明显减小。最后，对 模型来说，由于多年调节水库调节周期大于 年，模型在每年均采用相同的 的做法较多年调节水库更适用于年调节水库。图 不同 验证期水库出流 热力图对比 这三个模型的出流模拟效果以探究本文提出的两种模型结构的合理性。在各水库中，若 为最优模型，则说明在该水库两种模型结构起协同作用模型结构合理，因此同时考虑了两项模型结构的 最优，反之则说明两种模型结构相互拮抗，模型结构有待进一步优化。在起协同作用的水库中，若 ，即 为次优模型，则认为模型结构（）在水库的出流模拟中起主导作用，反之则是模型结构（）起主导作用。考虑到部分水库三个模型的 相近，因此只对最优模型和次优模型的 相差大于 的水库。统计结果表明，满足三个模型 相差较大的水库共有 座，在其中 的水库中，为最优模型，说明在大部分情况下，两种结构起协同作用，均能提高水库的出流模拟效果。在起协同作用的座水库中，为次优模型的水库座数为 个，说明在水库出流模拟中，主要由模型结构（）占主导作用。以三峡水库为例，探究引入 如何提高水库的出流模拟精度。对比 年的模拟结果（图 （），可以发现，在每年的 月明显高于其它模型，而在 月则相对更低，统计 年 年间各模型的多年平均日出流量也是如此（图 （）。的这一变化规律与 是一致的，对比三峡水库 与观测日入流量 的多年平均值（图 （）黑色实线和虚线）可以发现，在 月均高于 ，说明水库不断泄水，并且在 月中旬至 月中旬泄水幅度最大，体现了三峡水库在枯水期增大出流发挥经济效益以及在汛前腾出防洪库容的运行模式。另外，在 月以及 月显著低于 ，体现了三峡水库削峰防洪以及汛后蓄水的运行模式。而 的 （图（）点划线）也较好的描述了这一运行模式：月 不断降低，随之减少，出流增加，水库不断泄水；月 处于较低水平，最低约为 亿，与三峡水库汛限水位对应的 亿 较为接近；月，不断上升，随之增加，出流减少，水库不断蓄水。图 不同 模拟的三峡水库出流过程与观测出流、入流过程图 不同 验证期蓄水量 热力图 水库蓄水量模拟将各个 在验证期模拟的各水库蓄水量的 绘制成热力图，并按 值对水库进行降序排列（图 ）。结果发现，在季调节水库中，几乎所有各 的 都小于 ，说明所有模型都无法有效模拟季调节水库的蓄水量过程。原因可能有两点：一是季调节水库库容小，库容变化波动大，模拟难度大；二是目前采用的 在年内仅有一次蓄放水过程，而季调节水库调 节 周 期 短，在 年 内 会 经 历 多 次 蓄 放 水 过 程。在 座 年 调 节 水 库 和 多 年 调 节 水 库 中，、模拟精度显著高于其他 ，除了丹江口水库和隔河岩水库以外，在其他水库中均能取得较好的模拟效果。、的蓄水量模拟表现显著优于其他模型，各有 、座水库的 大于 ，而 模型中 大于 的水库则有 座，仍强于 模型以及其它 。说明两种结构均能提高水库的蓄水量模拟表现，但模型结构（）占主导作用，对蓄水量的影响更大。这也与模型结构（）通过实际蓄水量与目标库容的差值对水库出流产生影响的模型机理相吻合。水库在年内的多年平均日蓄水量 ，过程显示了水库在不同季节的蓄放水情况，在一定程度上可以反映不同水库的运行模式，因此采用 ，来评估模型是否能够适用于不同运行模式的水库。挑选来自不同流域的具有不同调节能力和功能的 座典型水库，绘制了 、模拟的以及实际观测的水库在年内的 ，过程（图 ），来探究水库的功能、调节能力如何影响水库的运行模式，以及各个模型对不同运行模式的水库的适用性。首先看水库功能、调节能力对水库运行模式的影响。防洪、发电、航运是长江三峡水库的主要功能，而三峡水库实测多年平均日蓄水量 ，（图 （）黑色实线）在 月维持在较低水平，在 、月份有所上升，在 月再次回落，说明水库在汛期将水位维持在汛限水位以留出库容削峰。另外，三峡水库 ，在 月底达到最大值，且维持到了次年 月，说明水库在该时期尽可能维持高水位，以提高水电站发电效益以及上游河道航运效益。对小浪底水库来说，冲淤是水库承担的重要任务，一般在每年的 月 日至 月 日进行，因此，（图 （）黑色实线）在对应日期急剧下降。承担防凌任务的刘家峡水库，需要在凌汛期之前腾清库容，因此实测蓄水量在 月迅速下降，以腾出有效防凌库容。刘家峡水库的最低 ，低于死库容，可能是由于刘家峡水库自 年竣工以来已运行多年，长年淤积使得库容有所损失所导致的。水库的调节能力不同，多年平均日蓄水量的变化规律也有所不同，主要体现在 ，的变化幅度之上，三峡作为年调节水库每年需要 个月蓄满水库，水布垭水库则需要 个月才能蓄满，而 大于 的龙滩、小浪底水库则无法保证每年都能蓄水至正常蓄水位（图 （）（）黑色虚线）。图 典型水库多年平均日蓄水量过程模拟结果对比各 模拟和实测的 座水库的 ，过程，可以发现，在三峡、水布垭、小浪底、刘家峡 座水库中，仅 ，和 ，的变化保持一致。在龙滩、光照水库中，虽然在量级上更靠近，但是这种量级的接近也只是多年平均后的结果，实际上 与 相差较大，远低于 （图 ）。总的来看，在这 座水库中，与 ，最为接近，同起同落且起伏幅度接近，但总是出现 ，偏大的情况，这可能是因为在 的构建中，将 函数的最小值作为参数 来率定，但是却将 函数的最大值固定为 的原因。另外，在刘家峡水库中，由于刘家峡水库承担防凌任务，使得刘家峡水库的 ，在年内存在两个波谷，而本文所构建的 ，仅考虑了每年蓄水量仅存在一个波峰波谷的情况。但即便如此，在刘家峡水库的模拟表现仍最优，和 远高于其他 。结论本文采用目标库容曲线 来描述水库在年内的运行过程，在水库出流模拟模型 中，通过为不同的水库设置相应的 ，以提高 对不同水库的适用性。从两方面考虑 对水库出流的影响，提出了两种模型结构：（）水库在目标库容较高的季节倾向于蓄水，在目标库容较低的季节则倾向于放水；（）水库的出流过程应使得实际蓄水量接近目标库容。根据这两种模型结构构建了 种基于 的 模拟了分布在长江、珠江、黄河三大流域上的功能与调节能力各不相同的 座水库的出流过程与蓄水量过程，并与常用的 相对比，结果发现：在出流模拟方面，所有 中，考虑了 的 模型在大部分水库的表现优于其他 ，尤其是在年调节水库和多年调节水库中，出流量的 值远大于其它 。综合考虑两方面影响的 模拟表现最优，出流量 值大于 的水库占 、而大于 的占 ，适用于大部分水库。但是，包括 模型在内，对多年调节水库的出流模拟精度仍有不足之处，值均小于 ，仍需要进一步的深入研究。在蓄水量模拟方面，所有 均不能较好模拟季调节水库的蓄水量过程，但在 座多年调节、年调节水库中，和 模型分别有 、座水库蓄水量的 值超过 ，蓄水量模拟精度远超其他模型。另外，在 座典型水库中，模型所模拟的多年平均蓄水量过程与实际观测过程同起同落，说明 模型可以有效模拟出功能不同的水库库容的季节性变化规律。总的来看，考虑了两种模型结构的 模型无论是在出流模拟、蓄水量模拟两方面表现均为最优，说明所提出的目标库容对水库出流的两种模型 结 构合 理 有 效。考 虑 后 的 模型对多年调节、年调节水库的出流、蓄水量模拟精度显著高于不考虑 的模型，且可以有效还原出其在年内的运行模式，显著提高对不同水库的适用性。参考文献：张玮，刘攀，刘志武，等变化环境下水库适应性调度研究进展与展望 水利学报，（）：，（）：，（）：李致家，黄鹏年，张建中，等新安江 海河模型的构建与应用 河海大学学报（自然科学版），（）：，（）：赵胤懋，廖卫红，雷晓辉，等考虑人工活动影响的流域水文模拟研究 水电能源科学，（）：孙新国，彭勇，张小丽，等基于聚合水库蓄放水模拟的洪水预报研究 水利学报，（）：，张珂，张企诺，陈新宇，等栅格新安江 地表地下双人工调蓄分布式水文模型 水资源保护，（）：，（）：董宁澎，余钟波，王浩，等耦合水库群参数化方案的区域陆面水文模拟 水科学进展，（）：梁桂星，崔亚莉，黄奇波，等 模型在无水库出流量资料流域的研究与应用 人民珠江，（）：，：，（）：，：，（）：郭生练，熊丰，王俊，等三峡水库运行期设计洪水及汛控水位初探 水利学报，（）：，雷晓辉，王浩，廖卫红，等非一致条件下的水文模拟及预报关键技术研究 北京：中国水利水电出版社，王旭，雷晓辉，蒋云钟，等基于可行空间搜索遗传算法的水库调度图优化 水利学报，（）：?，：，：，：（），（）：，：，（）：，：，：（下转第 页）?，（，；，）：?，?（），（），：；?；（责任编辑：韩昆）（上接第 页），（，）：，（），（）：（），（），：（），；（），：；（责任编辑：于福亮）