DreamBlooms/HansPub-15000-to-45000 huggingface dataset:
system: [SYSTEM]

user: 1. 引言 ESMD [ 1 ] [ 2 ] 是英文“Extreme-point Symmetric Mode Decomposition”的简称，其中文意思是“极点对称模态分解”。该方法是希尔伯特-黄变换方法 [ 3 ] [ 4 ] [ 5 ] 的新发展，由两部分组成：第一部分是模态分解，可以产生数个模态与一条最佳的自适应全局均线；第二部分是时-频分析 [ 6 ] 。本文主要针对ESMD方法的第一部分进行改进，原方法具体计算步骤如下： 该方法只考虑一维观测数据，即通常所谓的“时间序列”，记为 { Y ( t k ) } k = 1 N ，其中N是一个正整数。先是找出数据Y的所有极值点(极大值和极小值点)，对相邻极值点用线段连接并将线段中点依次记为 F i ( i = 1 ， 2 ， ⋯ ， n − 1 ) ，然后通过一定方式补充左、右边界中点 F 0 、 F n ，利用所获取的 n + 1 个中点构造p条插值线 L 1 ， ⋯ ， L p ( p ≥ 1 ) 并计算它们的均值曲线 L ∗ = ( L 1 + ⋯ + L p ) / p ，然后对 Y − L ∗ 重复上述步骤直到 | L ∗ | ≤ ε ，( ε 是预先设定的容许误差)或筛选次数达到了预先设定的最大值K。此时分解出第一个经验模 M 1 ，再对 Y − M 1 重复上述操作依次获得 M 2 ， M 3 ， ⋯ 直到最后余量R只剩一定数量的极值点。 最后的剩余模态R为“自适应全局均线(adaptive global mean curve)”，简称AGM曲线。事实上，只有当数据的拟合曲线是最佳的，去掉该曲线的剩余信号才可以被视为脉动量，进而可分解出一系列波动信号并认识其时-频变化特征。记原始信号为 Y = { y i } i = 1 N ，则其全局平均(数学期望)为 Y ¯ = 1 N ∑ i = 1 N y i (1) 相应的方差为 σ 0 2 = 1 N ∑ i = 1 N ( y i − Y ¯ ) 2 . (2) 这种常值形式的平均是最简单的全局均线形式。其实作为全局均线应当能够反应数据的总体变化趋势。这里的AGM曲线 R = { r i } i = 1 N 是由数据自动优选出来的，具有数据自适应特点。这里优选借用了“最小二乘”策略，此时的方差定义为 σ 2 = 1 N ∑ i = 1 N ( y i − r i ) 2 (3) 在实际应用中，我们通常选取 ε = 0.001 σ 0 ，采用标准差比率 v = σ σ 0 (4) 来反应AGM曲线R相对于全局平均 Y ¯ 的优化程度。 在数据分析中对称性与周期性关系密切，ESMD方法采用的是极点对称法则，这与我们以往所知道的周期的概念并不一样。一个极大值点与相邻的两个极小值点之间为一个周期，或者一个极小值点与相邻两个极大值点之间构成一个周期，也就是说两个相邻极值点之间是半个周期。极点对称反映的是连接极大值点和相邻极小值点的线段中点，所有对称中点都位于零值线上，这个中点称为局部中点。 ESMD方法适用于非线性、非平稳信号，如气候 [ 7 ] 、海气通量 [ 8 ] 等随机的时间序列数据，其在处理对称数据时取局部中点具有很好的代表性，分解之后的模态与自适应全局均线能准确的表达对称数据的信息与变化趋势。本文主要针对ESMD方法在处理非对称数据时所存在的缺陷对该方法进行改进。
system: [SYSTEM]

user: 2. ESMD方法在处理非对称数据时出现的问题及原因 基本上，符合时间序列的数据都可以用ESMD方法来处理，而且数据量越大效果越好。ESMD方法一般适应于在每半个周期中关于中点较对称(即极点对称)的数据。对于非对称数据(这时数据在形式上表现为峰谷不对称)来说，ESMD方法中选取的中点失去了代表性。峰谷不对称按照数学期望 [ 9 ] 来说，大部分数据点位于局部中点的一侧，只有一小部分位于另一侧，从而造成数据线峰谷不对称的现象。 本文针对非对称数据的一般示例进行讨论，即某一条河流近几十年的地表径流量数据。地表径流量的数据之所以不对称是因为春秋冬三季降雨量少，而夏季雨水较多。数据虽然具有较强的周期性，但是在每半个周期上并不关于局部中点对称，大部分都集中在局部中点的下方。这时如果用ESMD方法对数据进行处理，取中点并不能客观的代表原数据的趋势，反而会产生较大的偏差。所以最后求得的AGM曲线与原数据的标准差 σ 比全局平均与原数据的标准差 σ 0 要大，在筛选次数从1到40时，代入(4)式得到的标准差比率v都是大于100%的(如图1)。这种情况说明，数据的非对称性使得自适应全局均线反而比不上整体的平均值 [ 2 ] 。可以看出ESMD方法在处理类似于地表径流量这类非对称数据时存在极大的误差，从而严重影响最后结果的可信度。 图1. 筛选次数不同时的标准差比率(原方法)
system: [SYSTEM]

user: 3. ESMD改进算法的实验结果与分析 本文针对上述这种特殊情况，对ESMD方法进行改进，将这时不具有代表性的中点替换为能反应每半个周期数据情况的局部均值。每半个周期内的数据全部相加后除以相应的个数，即为局部均值。这时得到的自适应全局均线设为 R ′ = { r ′ i } i = 1 N ，设求得与原数据 Y = { y i } i = 1 N 的方差为 σ ′ 2 = 1 N ∑ i = 1 N ( y i − r ′ i ) 2 (5) 标准差比率为 v ′ = σ ′ / σ 0 . (6) 图1是原ESMD方法取局部中点时，求得筛选次数不同所对应不同的标准差比率。从图中可以看出最佳筛选次数为25次，对应的最小标准差比率为 v = σ / σ 0 = 104 .49% 。 图2是改进的ESMD方法取局部均值时，求得筛选次数不同所对应不同的标准差比率。从图中可以看出最佳筛选次数为23次，对应的最小标准差比率为 v ′ = σ ′ / σ 0 =98 .73% ，并且最大的方差比率也没有超过101.5%。 图2. 筛选次数不同时的标准差比率(改进的方法) 从1到40次筛选次数对应的方差比率可以看出，改进的ESMD方法明显比原方法的结果有了明显的优化，不仅方差比率降低了，筛选次数也相应减少了，从而节省了计算时间。 图3(a)表示原方法的模态分解，(b)表示原方法的数据与全局自适应均线的对比图，图4(a)表示改进方法的模态分解，(b)表示改进方法的数据与全局自适应均线的对比图。从图3(a)与图4(a)可以看出两种方法下的模态分解虽然有差距，但并不是很大。从图3(b)与图4(b)可以看出两种方法得到的全局自适应均线虽然形态近似，但是范围不同，原方法得到的全局自适应均线变化较剧烈，在极小值点几乎不变的情况下，受到极大值点的影响很明显，变化的范围在 4 × 10 3   m 3 到 8 × 10 3   m 3 之间；而改进的方法得到的自适应全局均线受极值点的影响较小，变化范围在 3 × 10 3   m 3 到 6 × 10 3   m 3 之间。 图3. 模态分解和自适应全局均线(原方法) 图4. 模态分解和自适应全局均线(改进的方法) 图5为原方法与改进方法得到的自适应全局均线的比较图放大后的部分图像，从图中可以看出改进方法后得到的AGM曲线明显比原方法的AGM曲线整体下移，且变化更加平缓。 图5. 原方法与改进方法自适应全局均线的比较 针对ESMD方法优化前与优化后，基于“最小二乘法”进行方差验证，比较结果如表1。全局平均与原数据的方差为13.32，局部中点插值后的数据与对应原数据的方差为12.50，局部均值插值后的数据与对应原数据的方差为8.06。 Table 1 System resulting data of standard experimen 全局平均 局部中点插值 局部均值插值 地表径流量 13.32 12.50 8.06 表1. 标准试验系统结果数据 从表1中可以看出，在非对称数据基础上，局部均值插值后的数据与原数据的方差小很多，从而证明局部均值比局部中点更具有代表性。 图6. 对原数据进行五点滑动平均处理 图6是对原数据进行五点滑动平均处理后的一部分图像，从图中可以看出五点滑动平均与原数据拟合较好，细节未丢失，但是并没有反映出数据的总体趋势，相比之下，ESMD方法既有数据不同模态下的变化趋势，并且最后得到的是自适应全局均线，反应数据的总体趋势，事实证明ESMD方法对数据分析方面非常全面。
system: [SYSTEM]

user: 4. 结论 本文提出的用局部均值代替局部中点对非对称数据进行ESMD方法的模态分解，编入程序后得到的结果明显优于之前的算法。最后筛选出的自适应全局均线更能准确的反应数据的总体变化趋势与数值变化范围。由此可看出，对于非对称数据处理，局部均值比局部中点更具有代表性。
system: [SYSTEM]

user: 摘 要 本文通过聚类分析，对云南省规模以上工业企业财务数据进行研究，对各行业之间的关系，和整体的财务经济状况，以及产业类别之间财务状况是否有明显差异进行分析研究。通过分析了解到云南省规模以上工业企业中烟草制品业、化学纤维制造业两个工业企业的各种经济效益与财务状况都比较好。整体看来制造业的状况比采矿业、电力、热力生产和供应业的财务状况好。 关键词 :规模以上企业，财务数据，聚类分析 Copyright © 2018 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 企业是国民经济的重要组成部分，2015年云南省工业生产总值达到3848.26亿元，占全省生产总值的28.26%。云南省第二产业生产总值有71.05%是由工业企业贡献的。可见工业企业对经济增长起着举足轻重的作用。工业企业的经济财务指标不仅关系到工业企业的生存与发展，也是影响国民经济可持续发展的重要因素之一。企业管理经营中，企业的财务指标对企业未来发展至关重要，因此有必要对云南省规模以上工业企业的财务数据进行研究。对其财务指标进行分析，反映公司业务活动的变化情况，从而综合分析企业的财务状况和经营成果。利用上述指标，很多研究人员做出了丰硕的成果，金硕元等使用因子分析的方法对规模以上工业的财务情况进行分析 [ 1 ] 。余宏对国有企业的主要财务指标进行了分析 [ 2 ] 。
system: [SYSTEM]

user: 2. 数据来源及说明 本文采用数据来源于中华人民共和国统计局编制的《云南省统计年鉴(2016版)》。本组数据内容是2015年云南省规模以上工业企业主要财务指标，计量单位：%。其中原始变量为八个基础指标。 2.1. 其变量名及含义 Table 1 Variable descriptio 变量名 变量含义 X1 综合经济效益指数 X2 总资产贡献率 X3 资本保值增值率 X4 资产负债率 X5 流动资产周转率 X6 成本费用利润率 X7 全员劳动生产率 X8 产品销售率 表1. 变量说明 经济效益指数的选取必须可以正确体现出企业的经济效益，所以以利润为主的指标将被使用以实现上述目标。以利润和实物相关指标挂钩的企业，应降低实物指标比重，逐步提高以利润为主的指标的比重。总资产贡献率反映企业全部资产的获利能力，是企业经营业绩和管理水平的集中体现，是评价和考核企业盈利能力的核心指标。资本保值增值率是财政部制定的评价企业经济效益的十大指标之一，资本保值增值率反映了企业净资本的变动情况，是企业发展能力的集中体现。资产负债率又称举债经营比率，它反映债权人发放贷款的安全程度的指标，通过将企业的负债总额与资产总额相比较得出，反映在企业全部资产中属于负债比率。流动资产周转率指企业一定时期内流动资产的周转次数，反映企业经营的状况，也反映资金利用效果和再生产循环的速度。成本费用利润率是企业全部生产投入与实现利润之间的比率。全员劳动生产率指平均每个职工在单位时间内创造的工业生产的最终结果，反映企业的生产效率和劳动投入的经济效益。产品销售率：反映产品已实现销售的程度，是分析产销链接的情况，研究产品满足社会需求的重要指标。考虑到以上指标(表1)是反映企业运行的核心财务数据，本文着重对以上指标进行分析。 2.2. 各行业名称 行业细分见表2。 Table 2 Sub-sector 按行业类别分 制造业 采矿业(1) 造纸和纸制品业 黑色金属矿采行业 印刷和记录媒介复制业 有色金属矿采行业 文教、工美、体育和娱乐用品制造业 非金属矿采行业 石油加工、炼焦和核燃料加工业 开采辅助活动 化学原料和化学制品制造业 电力、热力、燃气及水生产和供应业(3) 医药制造业 电力、热力生产和供应业 化学纤维制造业 燃气生产和供应业 橡胶和塑料制品业 水的生产和供应业 非金属矿物制品业 制造业(2) 黑色金属冶炼和压延加工业 农副食品加工业 有色金属冶炼和压延加工业 食品制造业 金属制品业 酒、饮料和精制茶制造业 通用设备制造业 烟草制品业 专用设备制造业 纺织业 汽车制造业 纺织服装、服饰业 铁路、船舶、航空航天和其他运输设备制造业 皮革、毛皮、羽毛及其制品和制鞋业 电气机械和器材制造业 木材加工和木、竹、藤、棕、草制品业 计算机、通信和其他电子设备制造业 家具制造业 仪器仪表制造业 表2. 各细分行业
system: [SYSTEM]

user: 3. 理论与实证分析 3.1. 聚类分析 聚类就是把在未知分类的情况下，相似的样品聚集成一类。聚类中有一种方法叫系统聚类法(或层次聚类法)，它的基本思想是：第一步，将n个样品各自聚为一类；第二步，用规定的距离方法将最小的两个距离的样品聚为一个新类，将把新类与其他类的距离计算出来，进行合并；第三步，重复上述合并，直到合并完全为止 [ 3 ] [ 4 ] [ 5 ] [ 6 ] 。 聚类使用的距离的方法最常见的是一下几种，即最短距离法，最长距离法，类平均法，重心法，离差平方和法 [ 7 ] [ 8 ] 。 最短距离法也就是类与类之间的距离定义为两类最近的距离，即 D K L = min i ∈ G k ， j ∈ G l d i j (1) 其中： d i j 表示第i个样品与第j个样品之间的距离； G i 表示类别； D K L 表示 G K 与 G L 之间的距离。 最长距离法也就是类与类之间的距离定义为两类最长的距离，即 D K L = max i ∈ G k ， j ∈ G l d i j (2) 其中公式(2)中各字母含义与公式(1)相同。 类平均法也就是把类与类之间的距离定义为所有物之间的距离的平均，即 D K L = 1 n k n l ∑ i ∈ G k ， j ∈ G l d i j (3) 其中 n K ， n L 分别表示类 G K ，类 G L 的个数。 重心法也就是把类别的重心之间的欧式距离，即 D M L = ( x ¯ K − x ¯ L ) ′ ( x ¯ K − x ¯ L ) (4) 其中： x ¯ K 是的 G K 重心； x ¯ L 是的 G L 重心。 在论文中，我们将采用的是类平均法，重心法和ward法。它的聚类步骤如下： 1) 定义样品与样品之间的距离，计算样品的初始距离矩阵 D ( 0 ) ； 2) 把 D ( 0 ) 中最短的两类合并，即 G K 与 G L 合并为一个新类，记作 G M ，即 G M = G K ∪ G L ； 3) 计算新类 G M 与其他任意类 G j 的距离的递推公式如下： D M J = min { min i ∈ G K ， j ∈ G J d i j ， min i ∈ G L ， j ∈ G J d i j } = min { D K J ， D L J } (5) 4) 将上述合并的距离矩阵算出，记作 D ( 1 ) ，重复上述步骤，直到所有元素聚为一类为止。 其中： d i j 表示第i个样品与第j个样品的距离； G 1 ， G 2 等表示类； D K L 为 G K 与 G L 之间的距离。 3.2. 实证分析 由类平均法聚类结果如图1，我们可以把类别整理为五类，第一类为烟草制造业；第二类为化学纤维制造；第三类包含热力电力娱乐；第四类铁路航空；剩下为第五类。 图1. 类平均法聚类 图2. 重心法聚类 由重心法聚类结果如图2所示，与类平均法相同，可以把上面结果分为五类。所包含的行业有一些细小差异，但是总体来讲与原有的类别出入都比较大。 Ward法聚类结果如图3，用ward法聚类与类平均法与重心法聚类结果有很大差异，ward法聚类，可以明显的把聚类结果分为三类：第一类包含烟草、制品业、化学纤维制造业、电力热力行业、娱乐用品制造业、铁路航空及其他运输设备制造业；第二类包含黑色金属采矿业、水的生产和供应业、有色金属采矿行业，汽车制造业，非金属矿采行业，酒水饮料制造业；剩下的为第三类。与原来的分类也不太相符。 图3. Word法聚类 总结三种方法聚类，发现聚类的结果与原来的类别都不符合，原来的类别是按照行业来分，通过各类财务指标再来聚类，可以看到烟草制造在各种方法中都单独被分为一类。从经济学的角度看，云南省的烟草制造行业是云南省的较为核心的产业，自1987年起连续13年，云南烟草上缴的税利一直占云南财政收入的70%以上。2016年全国烟草行业工作会上指出，“十二五”期间，全国烟草行业累计实现工商税利47，680亿元，累计上缴国家财政41，323亿元；云南省“十二五”云南省烟草产业创造的增加值约为8880亿元，占全行业的18.6%。根据财务指标来聚类，把烟草分为一类是比较合理的。从2014年开始云南省电力热力供应行业形式越来越好，各项财务指标反应的企业的财务状况较好，被分为一类比较合理。
system: [SYSTEM]

user: 4. 总结 本文通过使用多元统计的三种不同的聚类分析方法对云南省2015年规模以上工业企业的财务数据进行分析。结果显示云南省规模以上工业企业烟草制品与化学纤维行业在经济效应状况以及资本情况较其他行业都好，可以把他们归为一类。对于铁路、航空航天和其他设备制造业在资产保值增值方面较其他行业要突出，市场的认可度较高。对于一些生产日常生活用品的工业企业，它们的流动资产周转情况较好，这与生产的产品性质有很大关系。综合各项财务指标烟草制品业、化学纤维制造业两个工业企业的各种经济效益与财务状况都比较好。整体看来制造业的状况比采矿业与电力、热力生产和供应业的财务状况好。
system: [SYSTEM]

user: 摘 要 本文全面分析研究了国内较早开展生态补偿的新安江流域、九洲江流域、汀江–韩江流域和东江流域的生态补偿实践。研究发现目前国内流域生态补偿仍以上一级财政“输血式补偿”为主，缺乏流域水资源产权制度安排和流域生态保护产业优化升级。据此提出了穆棱河流域要通过设立完整的流域考核指标体系、设置“双指标”考核制度、制定“双向补偿制度”和构建流域水权制度建立完整的生态补偿制度，更好地推进穆棱河流域生态文明建设和全面实现小康社会目标。 关键词 :生态文明，生态补偿，指标体系，穆棱河流域 Copyright © 2018 by authors and Wuhan University. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 根据2015年《中国环境状况公报》，全国主要江河湖泊IV类~劣V类水质断面仍占35.5%，流域水污染形势依然十分严峻。为确保实现生态文明建设和全面建成小康社会的总体目标，国内各流域已经开始流域生态补偿的实践，其中比较有代表性的是新安江流域、九洲江流域、汀江–韩江流域和东江流域。 穆棱河作为国际界河——乌苏里江和国际界湖——兴凯湖的源头，如何保证其流域水资源生态安全具有十分重要的意义。科学合理地借鉴国内其他流域生态补偿的实践经验，有助于协调穆棱河流域上中下游之间的利益关系，保护流域生态系统服务功能，实现流域的可持续健康发展，推进流域生态文明建设。
system: [SYSTEM]

user: 2. 穆棱河流域概况 穆棱河发源于长白山支脉老爷岭山脉东坡穆棱市窝集岭北麓，是黑龙江省十大江河之一。穆棱河全长834 km，流域面积17，490 km 2 ，年径流量41.2亿m 3 。 穆棱河流域流经穆棱市和鸡西市。穆棱市2016年实现地区生产总值188亿元，三次产业结构为16:50.8:33.2 1 。鸡西市2016年实现地区生产总值530.1亿元，三次产业结构为35:24.8:40.2，第一、二、三产业对GDP增长的贡献率分别为38.5%、23%和38.5%。全市人均地区生产总值实现29，800元 2 。
system: [SYSTEM]

user: 3. 穆棱河流域水生态危机 3.1. 流域中下游水质严重污染 根据穆棱河流域的水域功能和水体标准，该水域应该属于II类水质，即主要适用于集中式生活饮水地表水源地一级保护区、珍惜水生生物栖息地、鱼虾类产卵场、仔稚幼鱼的索饵场等。但是，多年来的研究表明，穆棱河水质总体上仍属于IV类水体，一些河段的水体为V类和劣V类水质；穆棱河90%的河段水质达不到功能区水质标准 [ 1 ] [ 2 ] 。 3.2. 流域水资源过度开发利用 穆棱河上游的三岔村以下至兴凯湖和乌苏里江口，仅干流就有拦河坝27道，流域内有大中型水库23座。水资源过度开发利用导致大、小兴凯湖和乌苏里江的鱼类无法上溯穆棱河及其支流各个河段。众多的拦河坝不仅截断鱼类产卵、洄游和索饵的通道，而且阻止了乌苏里江和兴凯湖生物多样性向穆棱河流域扩散，导致穆棱河流域生物多样无法自然修复 3 。 3.3. 流域水土流失状况堪忧 穆棱河流域多年来对土地、植物资源实行掠夺性经营，形成大量坡耕地，流域水土流失面积达到230万亩。植被破坏造成土地质量下降，极易形成山洪造成水土流失，河道泥沙量增大。流域河道无序采砂导致河水浑浊，水质下降。河道内泥沙淤积严重，造成大量泥沙下泄，致使下游小兴凯湖底泥沙量剧增，可利用水资源率下降。
system: [SYSTEM]

user: 4. 国内流域水生态补偿实践 4.1. 新安江流域 新安江干流长度373 km，流域面积11，674 km 2 ，年平均降雨量900~1700 mm，年径流量72.3亿m 3 。 新安江流域生态补偿机制是从2004年开始酝酿的，2006年十届全国人大四次会议上正式提出《关于在新安江流域建设国家级生态示范区和构建“和谐流域”试点的建议》。 2012年9月安徽、浙江两省正式签订《新安江流域水环境补偿协议》，期限为2012年至2014年。2012年中央财政生态补偿资金3亿元，浙江省和安徽省生态补偿资金各1亿元全部到位。跨省界面水质达到考核标准后，浙江省补偿安徽省1亿元，达不到考核标准时，安徽省赔偿浙江省1亿元。中国首个国家层面的跨省界流域水生态补偿项目正式启动。 为建立新安江跨省流域生态修复的长效机制，2016年12月8日，安徽省和浙江省开始第二轮试点，三年补偿资金21亿元，安徽和浙江两省在水质考核指标提高7%的前提下，把各自的生态补偿资金额度每年增加到2亿元，中央三年补偿资金9亿元，按4亿、3亿、2亿递减方式补助。 2018年4月12日，《新安江流域上下游横向生态补偿试点绩效评估报告(2012~2017)》通过专家评审。该报告显示，2012年新安江流域生态补偿试点实施以来，流域上游水质为优，连年达到补偿标准，整个流域水环境质量均得到提升。 4.2. 九洲江流域 九洲江发源于广西省玉林市，在广东省湛江市的廉江注入北部湾，为南海北部湾水系。九洲江流域干流长度162 km，流域面积3337 km 2 ，年平均径流量为28亿m 3 。 4 解放后九州江流域相继得到规划开发，流域内共建有50余座大中小型水库，十几座拦河坝和多处大型泵站 5 。九州江流域的开发治理导致干支流出现不同程度的水环境生态危机。 2016年3月21日，在财政部、环保部和20余省市代表的见证下，广西壮族自治区政府与广东省政府签署了《广西壮族自治区人民政府广东省人民政府关于九洲江流域上下游横向生态补偿的协议》。补偿协议主要内容为：协议有效期为2015年至2017年；设立九洲江流域生态补偿资金6亿元(两省各出资3亿元)；中央财政设立奖励资金(广西省完成水环境保护目标后，中央财政奖励9亿元)；水质监测指标合格率由中国环境监测总站认定。九洲江流域水质监测以石角作为考核监测断面，跨省界断面年均值达到Ⅲ类水质，采用《地表水环境质量标准》(GB3838-02)中五项指标考核，即pH值、高锰酸盐指数、氨氮、总磷、五日生化需氧量。2015年、2016年和2017年考核目标为水质达标率分别达到60%、80%、100%。通过年均值和达标率双指标的设置，体现水质稳定达标、持续改善的要求 6 。 九洲江生态补偿有效遏制了流域水环境恶化的趋势，流域内的违法采砂、侵占湿地、过度养殖、畜禽污染、城镇生活污染、工业点源污染和农业面源污染等损害水生态环境的行为都得到较好控制。 4.3. 汀江–韩江流域 汀江发源于福建省西部，在广东境内与大埔县的梅江汇合形成韩江。汀江干流长度285 km，流域面积9022 km 2 ，年平均径流量115亿m 3 。汀江流域建设水利工程28，156处，流域水土流失比较严重，径流含沙量大 [ 3 ] 。 韩江是广东省除珠江之外的第二大河流，干流总长度470 km，流域面积17，851 km 2 ，年平均径流量245亿m 3 。韩江流域建设大中小型水库10，041座，流域水生态问题主要是水土流失 7 。 鉴于上游汀江的水质关系到下游1000多万人的饮水安全，2016年3月21日，福建省与广东省签订了《福建省人民政府广东省人民政府关于汀江-韩江流域上下游横向生态补偿的协议》。协议主要内容为：汀江-韩江流域水环境补偿考核期限暂定为2016~2017年；设立汀江-韩江流域水环境补偿资金4亿元，两省各承担50%；中央财政按照考核目标完成情况拨付给流域上游奖励资金；水质监测指标合格率由中国环境监测总站认定。考核目标为，采用《地表水环境质量标准》(GB3838-02)中五项指标考核，即pH值、高锰酸盐指数、氨氮、总磷、五日生化需氧量，2016年和2017年跨省断面的梅潭河、石窟河、汀江的水质年均值达III类标准，水质达标率100%；象洞溪跨省界断面水质年均值达V类标准，达到V类水质的比例为50%~100% 8 。 汀江—韩江流域水环境补偿协议创新之处为“双指标考核”、“双向补偿”，即采用污染物浓度和水质达标率双指标考核；福建省水质监测指标达标后由广东省补偿资金，福建省水质监测指标超标后，福建省赔偿广东省，跨省断面水质部分达标时按达标河流来水量比例和不达标河流来水量比例分别计算补偿金额。 4.4. 东江流域 东江是珠江水系三大河流之一，河流长度562公里，流域面积35，340平方公里，平均年径流量257亿立方米。东江流域是香港特别行政区和珠江三角洲经济圈的深圳、广州、东莞等城市的重要饮用水水源地。 鉴于东江流域水资源保护的极端重要性，2016年10月19日，在国家环境保护部和财政部见证下，江西和广东两省签订了《江西省人民政府广东省人民政府东江流域上下游横向生态补偿协议》，共同推进东江流域的水生态环境保护。协议确定了东江流域江西省和广东省横向生态补偿期限暂定三年。广西广东省跨省断面水质年均值达到III类水质标准，水质达标率逐年改善。 该生态补偿协议明确以庙咀里两个跨省界断面为考核监测断面。采用《地表水环境质量标准》(GB3838-02)中五项指标考核，即pH值、高锰酸盐指数、氨氮、总磷、五日生化需氧量，如出现其他特征污染物，广西和广东协商后再纳入考核指标。中国环境监测总站负责组织对跨界断面水质开展联合监测。两省共同设立水生态环境保护补偿资金6亿元，各自承担50%，两省共同监管生态补偿资金使用。资金补偿与水质监测考核结果相关，中央财政依据考核目标完成情况拨付生态补偿资金给东江流域上游，用于东江流域上游水污染防治和生态环境保护与建设工作 9 。
system: [SYSTEM]

user: 5. 国内流域生态补偿实践的创新和不足 5.1. 流域生态补偿实践的创新 5.1.1. 生态补偿为流域生态文明建设奠定了坚实的基础 党的十八届三中全会提出：“建设生态文明，必须建立系统完整的生态文明制度体系，用制度保护生态环境。”新安江流域、九洲江流域、汀江–韩江流域和东江流域开启了国内流域水生态补偿的先河，为建设流域生态文明奠定了坚实的基础，这些流域生态补偿的实践丰富了生态文明建设的理论内容，为生态文明建设提供了重要基础和保障。 5.1.2. 设置双指标考核制度 九洲江流域跨省界断面水质标准年均值设定为III类水质，除监测pH值、高锰酸盐指数、氨氮、总磷、五日生化需氧量等五项考核指标外，还根据流域污染现状、上游地区经济发展水平和环境综合整治能力设置了年度水质达标率，即2015年、2016年和2017年考核目标为水质达标率分别达到60%、80%和100%。通过年均值和达标率双指标的设置，实现水质持续改善、稳定达标的目的。 5.1.3. 建立双向补偿制度 汀江–韩江流域水生态补偿创新之处为“双向补偿”，在采用污染物浓度和水质达标率双指标考核时，上游的福建省水质监测指标达标后，由下游的广东省补偿资金；上游的福建省水质监测指标超标后，上游的福建省赔偿下游的广东省。跨省断面水质部分达标时按达标河流来水量比例和不达标河流来水量比例分别计算上下游补偿和赔偿金额。 5.2. 流域生态补偿实践的不足 5.2.1. 流域生态补偿全部是“输血式生态补偿” 以新安江流域、九洲江流域和汀江–韩江流域为例，生态补偿行动全部是政府主导，所需资金全部由上级政府承担，其中中央财政承担补偿资金的60%，省级财政承担40%。生态补偿资金仅仅用于保证部分污染源头的治理，集中在工业点源污染治理、农村面源污染治理、城镇污水及垃圾处理和生态修复工程，是典型的“输血式生态补偿”。没有在流域内进行“造血式生态补偿”，难以保证流域生态补偿具有可持续性，流域生态环境的改变缺乏长期性和稳定性。 5.2.2. 流域生态补偿的产权制度安排空白 水资源是一种具有外部效应性、公共属性和使用博弈属性等多重特性的稀缺资源，流域水资源补偿只有建立在水权制度的基础上，才能使水资源配置真正达到帕累托最优。目前国内流域水资源生态补偿实践仅仅考虑水质方面，没有将水资源数量因素纳入生态补偿范围。中共中央、国务院印发的《生态文明体制改革总体方案》要求“构建归属清晰、权责明确、监管有效的自然资源资产产权制度，着力解决自然资源所有者不到位、所有权边界模糊等问题。”流域水资源生态补偿产权制度安排，要在水资源所有权归属国家所有的前提下，在水质和水量两方面明确流域上下游各利益方的水资源使用权和国家代表的流域生态用水使用权，并对水资源使用权确权登记。 5.2.3. 流域生态补偿缺少生物多样性监测指标 目前国内流域生态补偿实践中，考核依据仅仅局限于pH值、高锰酸盐指数、氨氮、总磷、五日生化需氧量等几项化学监测指标，这并不能客观反映水生态系统的健康程度。流域生态补偿的最终目标是保护和修复自然生态系统，恢复水生态系统的生物多样性。党的十八大明确提出“要保护生物多样性。”全国生态保护“十三五”规划纲要提出“把加强生物多样性保护作为工作主线。”“到2020年，生态空间得到保障，生态质量有所提升，生态功能有所增强，生物多样性下降速度得到遏制。” 10 为全面反映水生生态系统保护和修复程度，必须把水生生物多样性指标纳入流域生态补偿考核指标体系之中。 5.2.4. 流域生态补偿缺乏调动地方积极性机制 目前国内流域生态补偿实践主要是上级行政部门，即流域省级和中央层次倡导和推动，缺乏流域基层地方政府和公众参与机制。缺乏地方政府和公众积极参与机制，很容易导致污染治理不力或者敷衍了事。中央环保督察组在抽查交办污染案件时屡屡发现地方政府对污染事件虚假整改、敷衍整改，从而导致污染反弹的问题 11 。流域上下游的市(县)、乡(镇)、村(屯)是流域自然生态系统保护和修复的直接参与和实践者，整个流域内从事经济和社会生产活动的人口是流域自然生态系统的有机组成部分，只有建立调动流域内基层政府和公众积极主动参与流域生态系统保护的机制，才能确保流域生态系统修复和保护的长久性和可持续性。 5.2.5. 流域生态补偿的产业转型升级空缺 目前的流域生态补偿实践注重的是污染的工程治理，生态补偿资金仅用于生态修复工程、城镇污水及垃圾处理、工业点源污染治理和农村面源污染治理，忽视了“造血式生态补偿”，即应该从流域整体上调整优化产业结构。为实现流域生态系统保护的可持续发展，流域生态补偿的重点应该是推动流域内“战略性新兴产业和先进制造业健康发展，采用先进适用节能低碳环保技术改造提升传统产业，发展壮大服务业，合理布局建设基础设施和基础产业。积极化解产能严重过剩矛盾，加强预警调控。” 12 要发挥流域内资源优势，发展再生能源和清洁能源，化解产能过剩和淘汰落后产能企业。科学配置流域资源要素，支持鼓励流域内优势产业走向国内，走向世界，提高参与国际化产业分工水平。
system: [SYSTEM]

user: 6. 基于生态补偿的穆棱河流域生态文明建设 为加强穆棱河流域污染治理，推进流域生态文明建设，2016年12月16日，黑龙江省财政厅、黑龙江省环境保护厅印发了《黑龙江省穆棱河和呼兰河流域跨行政区界水环境生态补偿办法》(黑财规审[ 2016 ]38号)，文件规定了穆棱河流域水环境生态补偿内涵，指省政府对造成水环境恶化的市县人民政府扣缴其下一年度财力，并将扣缴资金作为生态补偿资金，对主动改善水环境的市县人民政府给予横向补偿的活动。文件还规定水体污染物考核指标仅限于化学需氧量和氨氮两种指标 13 。相比国内其他流域生态补偿实践，穆棱河流域生态补偿办法还需要进一步补充完善，其他流域生态补偿实践经验为穆棱河流域治理和生态文明建设提供了很好的借鉴。 6.1. 完善流域生态补偿制度 目前的化学需氧量和氨氮两项考核指标过于单一和片面，无法全面客观反映流域水生态健康状况，考核指标体系应包含pH值、高锰酸盐指数、氨氮、总磷、五日生化需氧量、生物多样性等六项指标。 根据《地表水环境质量标准》(GB3838-02)规定，穆棱河流域应该达到II类水质，而目前水质现状令人堪忧，部分监测指标已经达到劣V类水质。为稳步推进流域生态系统改善，除监测pH值、高锰酸盐指数、氨氮、总磷、五日生化需氧量、生物多样性等六项考核指标外，还要根据流域污染现状、流域市(县)经济发展水平和环境综合整治能力设置穆棱河流域年度水质达标率，通过水质年均值和年度达标率双指标的设置，实现水质持续改善、稳定达到《地表水环境质量标准》设定的II类水质目标。 穆棱河流域涉及的市(县)为穆棱市、鸡西市区、鸡东县、密山市和虎林市等五个市(县)，目前生态补偿方式为扣缴水质超标市(县)的财力，补偿水质改善的市(县)。在采用污染物浓度和水质达标率双指标考核时，上游市(县)水质监测指标达标后，由下游的市(县)补偿资金；上游的市(县)水质监测指标超标后，上游的市(县)赔偿下游的市(县)。跨市(县)断面水质部分达标时，按达标河流来水量比例和不达标河流来水量比例分别计算上下游补偿和赔偿金额。 穆棱河流域水质污染、生态系统恶化，归根结底是由流域水资源的公共属性、外部效应性和使用博弈属性产生的“公地悲剧”导致的，从根本上解决流域水资源危机问题必须对流域水资源产权进行科学配置。应根据穆棱市、鸡西市区、鸡东县、密山市和虎林市等五个市(县)的人口数量、工农业生产情况和流域生态用水情况确定各个行政主体的水权配置，在水质和水量方面进行确权登记，并建立流域水市场，允许五个行政主体对各自的水权进行市场交易，通过市场主体实现改善流域生态环境的目标。 6.2. 实施流域环境整治工程 流域内城镇生活污水全部处理达标后才能排放进入穆棱河。加强工业循环水和再生水利用。农业面源污染治理中要禁止使用高浓度高残留农药，加强农用残膜清理回收、废弃农药包装物和废弃食用菌包装监管，积极推广测土配方和精准施肥技术，减少化肥施用数量和种类。全面进行农村环境整治，新建、改建和扩建的规模化畜禽养殖场要实行雨污分流，实现粪便污水资源化利用；统一规划、建设、管理农村污水处理设施。加强小流域治理，打击破坏河岸和沿河植被的违法行为；禁止在河道内洗砂、采金和洗煤行为。 6.3. 流域产业结构优化升级 优先发展优质高效益经济作物。充分发挥穆棱河流域山区及平原自然资源优势和对俄经济的区位优势，形成特色鲜明、优势突出、结构合理的穆棱河流域高效经济作物特色产业经济带。穆棱河上游是流域水生态涵养保护重点区域，优先发展中药材、向日葵、白瓜子、烤烟、林下经济和食用菌产业。穆棱河中下游是穆棱河生态环境保护和治理区域，优先发展瓜果、万寿菊、蔬菜等出口外向型的农业生产基地。 优先发展流域生态保护渔业。利用水生态系统的食物链和能量传递原理延伸渔业产业链，打造穆棱河流域保护水质、修复生态系统、优质高效的生态渔业。利用虎林和密山的口岸区位优势，发挥乌苏里江和兴凯湖的地缘优势，建设以兴凯湖大白鱼等淡水鱼为主题的特色外贸口岸，实现生态鱼养殖、加工和贸易一体化的特色产业。 推广综合节水利用农业。全流域稳步增加旱田节水灌溉面积，提高水资源重复利用率。 发挥区位优势，优化产业布局、加快产业升级，鼓励企业实施节水改造。穆棱河流域上游优先发展材料包装、木材加工、烟草和亚麻加工生产；穆棱河中游优先发展石墨和煤机制造产业，穆棱河下游优先发展水稻加工等绿色食品行业。 在冶金、煤炭、木材和石墨新材料等行业提高资源循环利用水平，实现废渣、废水、废气、能源和副产品的有效交换利用；发展清洁能源和新能源产业，利用流域丰富的水资源优势发展绿色食品产业，利用流域充分的太阳能资源发展光伏产业；推进流域循环经济园区和循环经济试点建设。 穆棱河流域具备丰富的旅游资源，可以依托著名的自然保护区(珍宝岛湿地国家级自然保护区、月牙湖国家级自然保护区、虎口湿地自然保护区和兴凯湖国家级自然保护区)资源、兴凯湖大白鱼、虎林和密山对俄口岸、北大荒开垦史、乌苏里江特产鱼类、虎头要塞遗址、北大荒书法长廊、循环经济产业园区等一大批特色旅游资源，着力打造穆棱河流域生态旅游品牌，把穆棱河流域整体打造成全国知名的生态旅游区、国际生态旅游养生度假区。
system: [SYSTEM]

user: 摘 要 在危化品零售企业，库存成本占企业物流成本的很大一部分。因此，如何有效的进行控制化学危险品库存已成了每个危险化学品企业不可避免的问题。本文在传统的(Q，R)模型中加入“风险成本”和“过期处理费用”建立符合危险品特性的库存控制模型，并对模型中各参数做出灵敏度分析。最后通过一个实例，应用matlab求解，验证了模型的实用性。 关键词 :危险化学品，风险成本，过期处理成本，库存控制，(Q，R)模型 Copyright © 2018 by author and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 危险化学品企业储存成本比普通货物高昂这已经是行业内的共识，其中乙类化学危险品仓库的租赁价格甚至是普通的几倍。这一部分的原因是化学危险品对物理储存条件要求非常高，特别是在安全、环保这两方面，这就使得仓库的使用条件更加严格和苛刻；另一部分的原因出在处理呆滞库存成本上面，由于是化学危险品的呆滞库存，这就必需专业的具备相应资质的处理公司接手，并且一定要在当地环保部门的全面监控下完成，这就使得成本进一步的增加了，有些处理费用甚至超过了采购价格的20%。由此可见，库存控制方法的探讨对化学危险品企业具有十分重要的借鉴及实践意义。 据统计，全球范围内每年生产的化学危险品约有40亿，危化品在生产、运输、储存、消耗过程中，存在着很大的安全隐患，一旦发生危险事故，将会对社会、经济、环境带来不可磨灭的影响。在存储方面，由于危险品的特性，其存储成本和要求比普通货物的要求高得多，尤其是在安全方面。因此制定出符合危险品特性的库存管理方法对危险品企业有着重要的实际意义。然而，国内目前对危险化学品库存优化问题的研究很少，大多都是定性的分析，很少有人针对危险化学品的库存控制提出定量的模型。本文设计了一种考虑包括风险成本和过期处理成本在内的总成本的化学危险品库存控制模型。
system: [SYSTEM]

user: 2. 文献综述 从经济订货批量模型被提出以来，随着应用的深入，人们考虑了更多周围环境变化所产生的影响，扩充了多种类的库存控制模型。Hadley和Whitin等人假设了供应商为刺激消费者需求，对产品采取价格折扣的策略，以此建立了具有价格折扣的EOQ模型。Crowther第一次从供应商和零售商双方角度出发，分析了一种批量购买时享受价格折扣的订货模型，从而使供应商的利益最大化的同时降低零售商的成本。Banerjee则考虑了两级供应链下的库存控制策略，以供需双方的库存费用最小为目标建立数学模型，并且具体分析了产品售价与订购量对库存成本的影响。陈宪章建立了季节性需求、变化生产率最优生产存储库存策略。研究发现，在不允许缺货的条件下，为满足季节性需求，单纯追求提前存储或增加机器数量，均不一定能够取得最佳效益，要根据存储费用和启动费用以及生产周期共同确定。在有限生产率的条件下，得出一种寻求最优生产–库存策略，为此类库存系统的管理决策提供了理论依据 [ 1 ]。施文武等建立了一种多周期随机需求生产/库存模型。该模型采用(s，Q)策略，设计了一种最优生产控制算法，根据该算法可以得出系统的最优生产准备点和最优生产量。并计算验证该方法可以有效地减小系统生产和库存的平均费用 [ 2 ]。王小斌讨论了含模糊参数的不确定环境下的经济订购批量问题与经济生产批量问题，运用可信性理论、不确定过程理论以及数学规划理论，分别建立了这些不确定库存问题的数学模型，并设计了相应的算法进行求解。张广霞通过联合库存在供应链库存控制中应用的研究，找出联合库存管理法在供应链库存控制中的应用对策，为实现供应链的同步化运作、零库存、准时化采购以及精细供应链管理提供条件和保证 [ 3 ]。张广霞通过联合库存在供应链库存控制中应用的研究，找出联合库存管理法在供应链库存控制中的应用对策，为实现供应链的同步化运作、零库存、准时化采购以及精细供应链管理提供条件和保证 [ 4 ]。蒋志明等将随机提前期作为变量，建立了一种基于随机提前期的(Q，R)模型，并用实例进行验 [ 5 ]。 目前针对危险化学品的研究集中在运输管理、安全风险评估、成本构成与管理等方面，而对危险化学品库存优化问题的研究很少，大多都是定性的分析。例如，王练以深圳的危险品库的库存管理为例，提出了一些危险化学品存储的建议 [ 6 ]。柴保身具体描述了我国危险化学品仓储业发展现状与展望 [ 7 ]。蔡翔根据危化品的特性，采用分库、分类、分区储存养护的办法对危险化学品进行库存管理，从而将风险降到最低 [ 8 ]。韩学飞将危害分析与关键控制点方法的预防性管理思想应用在危险品的仓储管理中，研究在危险品仓储过程中可能发生的危害，提出可行性的危害控制方法，对此方法在未来的发展方向进行了展望 [ 9 ]。刘立祥以唐山市危险品库存管理为例，给出了一些针对性的措施和建议 [ 10 ]。这些研究都是从危险化学品的特性出发定性的进行分析，得出库存管理策略。很少有将危化品的特性进行分析量化后加入库存模型中进行分析的，而准确的模型往往在现实生活更具有指导意义。本文将危险品的事故风险成本加入到库存总成本中，以库存总成本最小化为目标，建立库存控制模型，对库存控制进行了定量计算，求得最优订货量。
system: [SYSTEM]

user: 3. 数学建模 本文在(Q，R)模型中加入风险成本和过期处理成本，求解使得包括风险成本和过期处理成本在内的库存总成本最小的订货批量、再订购点和安全库存。我们考虑连续盘点策略，对企业的现有库存量进行连续盘点，一旦发现库存量到达或低于再订货点R，立刻实施补货策略，补货量为Q，以保证生产的连续。本模型以订货量Q和再订购点R作为决策变量，目标函数是最小化库存总成本。 3.1. 模型假设 1) 单位时间需求量服从正态分布。 2) 每次补完货的库存水平总是高于再订货点R，也就是一个提前期内只会订一次货，同时也保证了连续订货，并且不会发生交叉订货。 3.2. 模型符号意义 Q：订购量 v：单位危化品的价值 μ M D ：单位时间平均需求均值 σ M D ：单位时间需求标准差 L：提前期 ：提前期内需求均值， μ L D = L × μ M D σ L D ：提前期内需求标准差， σ L D = L × σ M D x：提前期内的需求量 f ( x ) ：提前期内需求量概率密度函数 F ( x ) ：提前期内需求量概率分布函数 P r ：单位时间内危化品发生事故的概率 C H ：每件危化品单位时间存储费 C O ：每次补货订购成本 C V ：单位时间单位危化品缺货成本 C r ：单位危化品风险成本 C m ：单位产品过期危化品的处理费用 3.3. 模型建立 1) 存储成本HC 由以上符号可得安全库存为： s s = R − μ L D (1) 则平均存储量为： I = Q 2 + s s = Q 2 + R − μ L D (2) 得单位时间存储成本为： H C = C H × I = C H × ( Q 2 + R − μ L D ) (3) 2) 订货成本OC 单位时间订货次数(周期数)为 ，可得单位时间订货成本为： (4) 3) 缺货成本VC 如果D1为提前期的需求量，则提前期内的缺货量为 max ( D 1 − R ， 0 ) ，由此可得一个周期内的期望缺货量为： B ( R ) = ∫ R ∞ ( x − R ) f ( x ) d x (5) 则单位时间的缺货成本为： V C = μ M D Q × C V × B ( R ) = C V μ M D Q × ∫ R ∞ ( x − R ) f ( x ) d x (6) 4) 风险成本RC R C = P r × C r × I = P r C r ( Q 2 + R − μ L D ) (7) 5) 过期危化品处理成本MC 危化品一旦过期只能作销毁处理，企业不但要承担由于过期所带来的损失，还要支付相应的处理费用。 (8) 6) 单位时间平均库存总成本 由以上计算步骤，可得年平均库存总成本为： T C ( Q ， R ) = H C + V C + O C + R C + M C = C H ( Q 2 + R − μ L D ) + C V μ M D Q × ∫ R ∞ ( x − R ) f ( x ) d x + μ M D Q C O     + P r C r ( Q 2 + R − μ L D ) + ( C m + v ) ( Q − μ M D μ M D Q ) (9) 7) 再订购点，订货批量和安全库存 对(8)中的Q和R分别求一阶偏导并令其等于0得： Q = 2 μ M D [ C o + C V ∫ R ∞ ( x − R ) f ( x ) d x ] C H + P r C r + 2 ( C m + v ) = 2 μ M D [ C o + C V B ( R ) ] C H + P r C r + 2 ( C m + v ) (10) R = F − 1 ( 1 − ( C H + P r C r ) Q C V μ M D ) (11) 对(9)和(10)交替迭代求解可得最优订货批量 Q * 、再订购点 R * ，步骤如下： ① 在式(9)中令 B ( R ) = 0 作为初始解，求解得 Q 1 ； ② 将 的值代入式(10)，求解得 R 1 ； ③ 将 R 1 的值代入式(5)计算 B ( R 1 ) ； ④ 再将 B ( R 1 ) 的值代入式(9)求得 Q 2 ； ⑤ 重复②~④步，直到 Q i 和 R i 的值基本上不再有较大的变化为止；此时的 Q i ， R i 即为最优解 Q * ， R * 。 再由 R * 可得到安全库存为： s s * = R * − μ L D (12)
system: [SYSTEM]

user: 4. 算例验证 本小节给出一个具体实例验证上述模型和求解算法的有效性。 已知某公司危险品每月需求量服从均值等于30，923，标准差等于6918的正态分布。此条件满足(Q，R)模型假设。已知订货提前期L平均为7.5天，即0.25个月，单位危化品的价值v为68元；单位危化品每月储存费为2元，单位危化品缺货费为10元，每次订货费为500元，单位危化品发生事故后造成的损失为3000元；单位时间内危化品发生事故的概率为0.0005；单位过期危化品的处理费用为5元。 因为该产品每月需求量的均值 μ M D = 30923 ，标准差 L = 0.25 ， L = 0.25 ，所以，提前期内需求量服从均值 μ L D = L × μ M D = 7731 ，标准差 σ L D = L × σ M D = 3459 的正态分布。 计算过程如表1所示。 Table 1 Process of iterative calculation of B ( R i ) ， Q i ， R B ( R i ) Q i R i 0 3016.667982 14，542.57619 46.64653 4194.069709 13，996.65561 68.11737 4636.597491 13，823.50327 76.51286 4798.548057 13，763.42684 表1. B ( R i ) ， Q i ， R i 迭代计算过程 Continued 由上表我们得到最优订货量 Q ∗ = 4890 ，最优再订购点 R ∗ = 13730 ，安全库存 s s ∗ = 5893 。
system: [SYSTEM]

user: 5. 灵敏度分析 通过Matlab数值计算对模型中各个参数( C H 、 C V 、 C O 、 C r 、 P r )进行灵敏度分析，得到了每个参数在选取不同值时库存总成本TC的变化情况，结果如下。 5.1. 存储费C H 的灵敏度分析 图1. 存储费C H 灵敏度曲线 图1中的曲线揭示了不同存储费C H 对应的库存总费用TC。从图中我们可以看出，单位产品存储费增加时，库存总成本随之增加。且存储费增加时，企业更愿意进行小批量订货，再订购点和安全库存水平也会设置为更低水平。由于危险化学品本身的特性，其存储费会比普通货物高得多，所以企业可以通过控制平均库存量在一个较低的水平来降低库存成本。 5.2. 缺货成本C V 的灵敏度分析 由图2可知，单位产品的缺货成本增加时，库存总成本随之增大，且缺货成本增加时，企业会设置较高的再订购点和安全库存水平以期降低因缺货给企业带来的损失。 5.3. 订购费C O 的灵敏度分析 由图3可知，当订购费C O 增加时，库存总费用随之增大。且当每次的订购费用较大时，企业不愿频繁订货，并会将再订购点和安全库存水平设置较低水平。 5.4. 风险损失C r 的灵敏度分 由图4可知，随着单位原材料风险损失C r 的增加，库存总费用相应增加。且风险损失增加时，企业更愿意订购较少的批量并设置较低的再订购点和安全库存水平。 图2. 缺货成本C V 灵敏度曲线 图3. 订购费C O 的灵敏度曲线 图4. 风险损失C r 的灵敏度曲线 5.5. 发生事故概率P r 的灵敏度分析 图5. 发生事故概率P r 的灵敏度曲线 由图5可知，当危险事故发生的概率P r 逐渐增大时，库存总费用随之增大。且发生事故的概率越大，企业为尽可能降低发生危险事故时的经济损失，会降低订购量并保持较低的安全库存水平。 5.6. 单位过期危化品处理费用C m 的敏感度分析 图6. 单位过期危化品处理费用C m 的灵敏度曲线 由图6可知，当单位过期危化品处理费用C m 逐渐增大时，库存总费用随之增大。 5.7. 比较得出最为敏感的参数 下面我们把不同参数的单位变化导致的库存总费用的平均变化率进行比较，得出相对敏感的参数。 Table 2 Average change rate of total inventory cost caused by unit changing of parameter 参数 C H C V C O C r P r C m 库存总费用平均变化率 7255.3 338.2 5.3 3.3 15，252，250.2 1118.4 表2. 不同参数单位变化导致库存总成本变化率 由上表2我们可以看出，随着参数P r 的变化，库存总成本TC的变化率最大，即发生事故的概率P r 是相比之下最为敏感的参数。因此P r 对企业的库存费用的控制至关重要，企业应关注产品的仓储管理以降低危险事故发生概率，从而降低库存总成本。
system: [SYSTEM]

user: 6. 结论与展望 相比于一般物资，化学危险品对储存环境、安全以及环保等多方面的诸多要求，导致化学危险品库存持有成本很高。另外，报废化学危险品的处理程序复杂、费用可观，高效的库存管理体系对化学危险品企业以及整个社会均有着重要的经济意义、社会意义和环境意义。危险化学品的事故风险成本和过期处理成本是与其他库存控制策略相区别的重要方面，是在进行库存优化过程中必须考虑的成本，本文设计的库存控制模型加入其相关内容，建立了更加符合化学品特性的库存控制模型，并带入数值进行验证，证明了模型的实用性。关于如何合理评估危险品库存的事故风险，将是今后研究的重点问题。
system: [SYSTEM]

user: 摘 要 本文通过使用美国国家环境预报中心能源部(National Centers for Environmental Prediction-Department of Energy，简称：NCEP-DOE)提供的1979~2014年的再分析资料，利用谐波分析法将北半球中纬度大气按波数分类，再借助国际上通用的T & M方法检索欧亚关键区阻高，系统研究中纬度大气1~3波(超长波)的时空变化特征及其与欧亚关键区阻高的关系，主要结论如下：北半球中纬度地区(30˚N~60˚N) 1~8波波谱比重在0.93~0.96变化，其中1~3波波谱比重超过0.5，表明超长波是大气运动中最主要的组成部分，且30˚N和60˚N两条纬圈上超长波的波谱比重有准2年的周期变化，冬季振荡更为显著；当1~3波波谱比重偏强时，乌山地区位势高度(geopotential height，简称：GPH)明显偏低，贝湖地区的偏高，鄂海地区GPH变化不明显，反之，鄂海地区GPH偏强，其他两个区域GPH变化不显著，尤其60˚N上1~3波波谱比重异常时，欧亚地区GPH异常表现更为突出；1波(3波)振幅与乌山阻高指数呈正(负)相关关系，与鄂海的呈负(正)相关，2波振幅的激增可能会导致乌山阻高持续时间增长。 关键词 :阻塞高压，谐波分析，超长波，波谱比重，振幅 Copyright © 2018 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/