第１７卷２０ｌ０年第３期５月安全与环境工程ＳａｆｅｔｙａｎｄＥｎｖｉｒｏｎｍｅｎｔａｌＥｎｇｉｎｅｅｒｉｎｇＶ０１．１７ＮＯ．３２０１０Ｍａｖ大风覆冰条件下输电杆塔可靠性模型的研究李致宇（武汉大学测绘学院，武汉４３００７９）摘要：大风及冻雨天气严重影响电力系统的安全运行，由此引发的线路覆冰现象会加重导线与输电杆塔的外力负荷，导致线路舞动甚至塔架垮塌。本文分析对比了若干有关输电塔架安全性与可靠度的研究成果．结合可靠性理论，提出了一种考虑大风与覆冰因素影响的塔架功能函数，并建立了输电杆塔的可靠性分析模型，选用基于一次二阶矩的映射变化法解算其可靠性指标，以为电力系统的安全分析和风险评估提供理论依据。关键词：输电杆塔；大风覆冰；功能函数；可靠性中图分类号：Ｘ９３；ＴＭ７５文献标识码：Ａ文章编号：１６７１—１５５６（２０１０）０３—００５９—０５ＳｔｕｄｙｏｎｔｈｅＲｅｌｉａｂｉｌｉｔｙＭｏｄｅｌｏｆＴｒａｎｓｍｉｓｓｉｏｎＴｏｗｅｒＩｍｐａｃｔｅｄｂｙＷｉｎｄａｎｄＩｃｅＬＩＺｈｉ—ｙｕ（ＳｃｈｏｏｌｏｆＧｅｏｄｅｓｙａｎｄＧｅｏｍａｔｉｃｓ，ＷｕｈａｎＵｎｉｖｅｒｓｉｔｙ，Ｗｕｈａｎ４３００７９，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｏｕｒｐｏｗｅｒｓｙｓｔｅｍｉｓｔｈｒｅａｔｅｎｅｄｂｙｔｈｅｅｘｔｒｅｍｅｌｙｂａｄｗｅａｔｈｅｒｌｉｋｅｓｔｒｏｎｇｗｉｎｄａｎｄｉｃｅｒａｉｎ．Ｔｈｅｉｃｅｗｒａｐｐｅｄａｒｏｕｎｄｔｈｅｆｒｏｚｅｎｅｌｅｃｔｒｉｃａｌｃａｂｌｅａｇｇｒａｖａｔｅｓｅｘｔｒａｌｏａｄｔｏｔｈｅｃｏｎｄｕｃｔｏｒｓ，ａｎｄｔｈｉｓｃｏｕｌｄｒｅｓｕｌｔｉｎｗｉｒｅｗａｖｉｎｇａｎｄｔｒａｎｓｍｉｓｓｉｏｎｔｏｗｅｒｃｏｌｌａｐｓｅ．Ｔｈｉｓｐａｐｅｒｃｏｍｐａｒｅｓｓｏｍｅｃｏｎｃｌｕｓｉｏｎｓｏｆｓｅｖｅｒａｌｒｅｌａｔｅｄｔｈｅｓｅｓ，ａｎｄｐｒｅｓｅｎｔｓｏｎａｐｅｒｆｏｒｍａｎｃｅｆｕｎｃｔｉｏｎｃｏｎｓｉｄｅｒｉｎｇｔｈｅｔｗｏｓｉｍｕｌｔａｎｅｏｕｓｆａｃｔｏｒｓ，ｗｉｎｄａｎｄｉｃｅ，ｂａｓｅｄａｎｔｈｅｒｅｌｉａｂｉｌｉｔｙｔｈｅｏｒｙ．Ａｆｔｅｒａｎａｌｙｓｉｓｏｆｔｈｅｐｅｒｆｏｒｍａｎｃｅｆｕｎｃｔｉｏｎ，ｔｈｅＦＯＲＭｐｒｏｃｅｄｕｒｅｗａｓｕｓｅｄｔｏｃａｌｃｕｌａｔｅｔｈｅｒｅｌｉａｂｉｌｉｔｙｉｎｄｅｘ，Ｋｅｙｗｏｒｄｓ：ｔｒａｎｓｍｉｓｓｉｏｎｔｏｗｅｒ；ｗｉｎｄａｎｄｉｃｅ；ｐｅｒｆｏｒｍａｎｃｅｆｕｎｃｔｉｏｎ；ｒｅｌｉａｂｉｌｉｔｙ电塔结构功能函数，并建立输电杆塔的可靠度分析Ｏ引言在对电力系统进行安全分析和风险评估时，恶模型，以为电力系统的可靠性分析和风险评估提供理论依据。劣天气的影响是不容忽视的，特别是大风及冻雨往往会使线路产生覆冰。覆冰会增加线路的迎风面积和外力负荷，轻则线路摆动，重则塔架垮塌，造成输电中断，甚至会造成人员伤亡的惨剧发生。如２００８年初我国南方突遭冰灾侵袭，ｌｏ多个省的电力供应局部中断，严重破坏了电力系统的正常运转，带来了巨大的经济损失。笔者基于可靠性原理，对比分析了国内外若干针对大风、覆冰条件下输电杆塔可靠性理论的研究成果，尝试提出了一种考虑大风与覆冰因素影响的１结构可靠性理论１．１结构的可靠性指标与可靠度结构的可靠性是针对结构功能而言的。结构的功能包括负荷能力、适用性能、耐久性能等指标【１］。一个结构的功能极限状态可用其功能函数等于零予以表达：Ｚ一厂（Ｘ）＝Ｏ（１）式中：ｘ＝［Ｘ，，Ｘ。，…，ｘ。］１’为随机向量。结构的可靠度就是结构处于正常运行状态的概收稿日期：２００９—１１—０３修回日期：２００９—１２—１０作者简介：李致宇（１９８５一），男。硕士研究生，主要研究方向为地学分析在安全应急领域的应用。Ｅ‘ｍａｉｌ：ｇｒｅａｔｌｚｙ＠ｓｉｎａ．ｃｏｒｎ万方数据６０安全与环境工程第１７卷率，即ＰＲ＝Ｐ（Ｚ＞０）（２）当满足如下条件时：①随机变量Ｘ・，Ｘｚ，…，墨均服从正态分布；②随机变量Ｘ，，Ｘｚ，…，Ｘ。相互；③Ｚ一，（Ｘ）为线性函数。结构的可靠性指标可表示为［１３（３）式中：∥。为随机变量Ｚ的期望值；屹为随机变量Ｚ的标准差。根据结构可靠性指标口，则结构的可靠度与结构失效概率分别为ＰＲ。Ｐ（Ｚ＞０）一西（ｐ）（４）Ｐｖ２Ｐ（Ｚ＜Ｏ）一西（一Ｄ（５）１．２结构可靠性指标的计算方法如上所述，对于满足上述３个条件的结构功能函数，可以按照公式（３）计算结构可靠性指标口。但实际情况中为线性函数的要求很难满足，所以要对功能函数ｚ在某点Ｘ＝Ｉｘ。，ｚ。，…，ｚ。］Ｔ进行泰勒展开并省略二次及更高次项，再计算结构可靠性指标Ｊ９。根据选取的展开点不同，一次二阶矩法可分为均值法和设计验算点法。均值法选取均值点ｚ＝匝。，，Ｐ屯，…，户‘］Ｔ作为展开点，也称中心点法，但由于均值点往往不在极限状态面上，且对于意义相同但表现形式不同的功能函数会得到不同的可靠性指标，所以误差较大［２３；设计验算点法选取位于极限状态面上的点ｚ。＝［ｚｆ，ｚ；，…，ｚ：］Ｔ作为展开点，很好地解决了均值法存在的问题，所以也称改进一次二阶矩法。随机变量为非正态分布的情况是经常遇到的，如结构的截面抗力一般服从对数正态分布，风速、风压和雪载等服从极限Ｉ型分布（也称Ｇｕｍｂｅｌ分布）［１１。对于由诸如此类非正态随机变量组成的结本文的研究对象为大风和覆冰条件下电塔的可设，２７。２［ｚ？，ｚ；，…，ｚ：］Ｔ满足功能极限状态万方数据万程Ｚ２／（ｚ’）。０，即Ｘ’征敬限状态回上。将么在Ｘ。泰勒展开并省略二次及更高次项，则有磊叫ｚ。）＋；耋『甏ｌ州（ｘ，叫）］（６）此时ＺＬ是由相互的正态随机变量Ｘ。，Ｘｚ，…，Ｘ。线性组合而成，所以Ｚ。的均值与方差分别为旷ｍ计喜［爰ｂ（旷刺］（７）噍２ｉ耋（爰ｌ㈣ｘ；＝Ｃａｘｉ）２利用构建状态极限方程Ｚ。＝０，并以标准化变量ｙｉ一！竺二垒２代换ｘ，并拆分，有（９）式两边同除Ｏ＇Ｚｋ，利用式（７）和式（８），并根据十攀塑一。㈣，。店ｃ甏ｋ？叫２定义Ｘｒ的灵敏度系数口ｘ．如下：％气。８铲藤＂Ｄ惹Ｚ蒂ａｘ；）ｚ一爰Ｌ．盯＿ａ×Ｉ。一．－。“‘Ｑ。ｘ；：。？为极限状态函数ｚ对其第ｉ个变量式（１１）则可写成（１２）ｆ＝１ｃｏｓＯｙ．Ｙｆ－－８＝０式（１２）可理解为标准正态随机变量ｙ空间的ｙ，２８ｃｏｓ＠（１３）在Ｘ空间中ｚ’＝Ｅｘ；，ｚ≠，…，ｚ：］Ｔ为Ｘｉ＊２ｐｘ：＋ｐｘｃｏｓＯｘ（１４）此时ｚ，与卢共行＋１个未知数，方程只有，ｚ驴ｍ计量矧州‰ｒ叫’］＋㈣；耋（爰ｋ掣ｔ）２。公式（３），计算功能函数磊的可靠性指标卢为构功能函数，在利用改进一次二阶矩法求解可靠性指标口时，要将非正态随机变量当量正态化或变换为正态随机变量。这些方法与改进一次二阶矩法结合，就形成了ＪＣ法、映射变换法和实用分析法等。对于相关随机变量不的情况，可用正交变换法或广义随机空间理论解决。靠性的计算方法，其结构功能函数多为含有非正态随机变量的非线性形式‘“。针对这一特点，下面将重点阐述改进一次二阶矩法和映射变换法。１．２．１改进一次二阶矩法式中：爰ｌＸｉ求偏导后并带人Ｘｉ。计算后所得的系数；％为随法线式平面方程，原点Ｏ到极限状态面上点Ｙ。是长度为卢的法线，方向余弦为ｃｏｓＯｖ二［２ｌ。可靠性指标卢是原点到极限状态面的最短距离，此时极限状态面上的点ｙ。（Ｘ空间中该点为ｚ’）称为设计验算点。验算点Ｙ。一Ｅｙ；，ｙ≠，…，ｙ：］Ｔ的坐标为个。所以采用迭代法求解，迭代计算步骤如下：①给ｙＯｓｏｃ１（第３期李致宇：大风覆冰条件下输电杆塔可靠性模型的研究定初始验算点（ｚ。）“’，可取（卫。）“’＝，ｕｘ（ｉ＝０）；②利用（１１）式计算Ｘ；的灵敏度系数ａｘ（即方向余弦ｃｏｓｏｘ．）；③利用（７）式和（８）式，结合（３）式计算可靠性指标口；④利用（１４）式计算新的ｚ。；⑤以新的ｚ。进行下一次迭代（ｚ’）“＋１’＝ｚ。，重复第②至第④步骤直至满足两次ＪＩｚ。ｌ｜之差小于允许的ｅ；⑥输出最终的验算点Ｘ’与可靠性指标口。得到可靠性指标卢后可利用公式（４）、（５）计算结构可靠度Ｒ和失效概率Ｐ，。１．２．２映射变换法映射变换法针对随机变量为非正态的情况，进一步完善了一次二阶矩法，其计算量与国际结构安全度联合委员会推荐的ＪＣ法相当，但数学理论上更加完整［１］。功能函数Ｚ＝厂（Ｘ）含有非正态随机变量Ｘ；，直接利用标准正态随机变量ｙ做映射变换，则有Ｅ（Ｘ）一西（Ｙ；）（１５）ｘ；一Ｆ７１嘧（ｙ；）］（１６）Ｙｆ一町１ＥＦ（Ｘｉ）］（１７）那么原功能函数Ｚ＝厂（Ｘ）转化为Ｚ—ｆ（Ｆ－１［①（ｙ）］）（１８）功能函数的自变量由非正态Ｘ映射为标准正态ｙ，映射后的功能函数记为Ｇ（ｙ）。注意标准正态随机变量ｙ；的前两阶矩分别为０和１，则相应地改造迭代公式（１１），则有一ａ＿ｑｑ｝．＝—：＝＝＝＝＝＝＝＝＝＝＝＝：＿＝兰型尘垒一－＂＂＂ｉＪ（）９）‘店ｃ簧Ｌ？训２其中求偏导因子为ＯＧ—ａｆａＸｉ（２０）ａＹｆａＸｉａｙｉ迭代过程的其他步骤和计算公式与１．２．１节所述相同。２电塔结构功能函数模型及可靠性指标计算电力网中对输电线路起到支撑与连接作用的设备大致分为输电塔架与电杆两种。其中塔架主要用于远程高压输电，大型塔架高度可达５０ＩＴｌ以上，一侧连接导线平均长度约为４００～５００１Ｔｌ，承载电压多为２２０ｋＶ、５００ｋＶ中高压，甚至８００ｋＶ、１０００ｋＶ特高压，布设环境远离市区；电杆主要用于近距离输电，杆高ｉ０ＩＴｌ左右，连接导线长度数十米，往往布万方数据设在城市中与最终用户的变压器相连。下面以输电塔架中使用最多的干字形塔架为例，研究大风及线路覆冰对于塔架结构功能的影响。２．１电塔功能函数模型电塔结构功能函数模型一般形式为［１］Ｇ（ｚ）一Ｃ（ｚ。，ｄ，）一Ｄ（ｚｄ，ｄｄ）（２１）即为承载能力Ｃ和实际荷载Ｄ的差值。其中ｚ。和ｚ。表示所有随机变量和不确定性因素，ｄ。和ｄ一表示所有确定变量。Ｇ（ｚ）＝０就是式（１）所述的结构的极限状态方程。对于干字形输电塔架，Ｒ．Ｏ．ＦｏｓｃｈｉＥ如曾提出以承载能力、风速、覆冰厚度为相互随机变量的结构功能函数，即Ｇ（Ｒ，口，ｉ）一Ｒ一［Ｋ∞＿ｒＴ（口，ｉ）Ｓ＋口ｖ（Ｖ（ｉ）＋Ｄ）Ｓ・ｐ＋ｄ＋ＫｒＣｓａＰＰ（ｕ）］（２２）式中：Ｒ为塔架承载能力随机变量（ｋＮ），服从对数正态分布；口为风速随机变量（ｍ／ｓ），服从极限Ｉ型分布；ｉ为覆冰厚度随机变量（ｒａｍ），服从极限Ｉ型分布；Ｋｃ为高度、地形、风偏角对导线载荷的修正系数；口ｒ为单位横向荷载系数；丁（７．／，ｉ）为单位长度导线因风速和覆冰产生的横向荷载（ｋＮ／ｍ）；Ｓ为水平档距（ｍ）；ａｙ为单位纵向载荷系数；Ｖ（ｉ）为单位长度导线因覆冰产生的竖向荷载（ｋＮ／ｍ）；Ｄ为单位长度导线产生的静荷载（ｋＮ／ｍ）；ｐ为重量跨度与风跨度的比值；ｄ为由构件产生的静荷载（ｋＮ）；Ｋｒ为高度、地形、风偏角对塔架载荷修正系数；Ｃｓ为塔架形状、防护效果修正系数；卿为单位风压在塔架竖直投影面上产生的荷载（ｋＮ／ｍ２）；Ｐ（ｕ）为风速可下的理论风压（ｋＮ／ｍ２）。在Ｒ．０．Ｆｏｓｃｈｉ的研究基础上，文献［４］给出了Ｔ（ｔ，，ｉ）的具体形式：＿ｔ２丁（ｕ，ｉ）＝七（ｄｃ＋２ｉ）ａ杀ｇ（２３）上Ｕ式中：ｕ为风速（ｍ／ｓ）；ｉ为覆冰厚度（ｒａｍ）；ｄｃ为导线直径（ｃｍ）；ｋ为风载体型系数；口为风速不均匀系数；ｇ为重力加速度（ｍ／ｓ２）。但文献［４］既未将覆冰厚度ｉ作为随机变量看待，也未将其作为普通变量看待，而是设覆冰厚度ｉ＝０，忽略了覆冰因素，相应的结构功能函数Ｇ（Ｒ，ｔ，）仅考虑了塔架承载能力与风速的统计性质，且认为结构承载能力Ｒ与风速ｕ服从正态分布，而不是Ｒ．０。Ｆｏｓｃｈｉ提出的对数正态分布与极限Ｉ型分布［３］。这样的简化处理可以利用最简单的均值一次二阶矩计算结构的可靠性指标口，但如上所述，均值一次二阶矩法在计算精度上有时不能满足要求，６２安全与环境工程第１７卷文献［２］中指出均值一次二阶矩法对结构可靠性指标口在１至２时的计算结果比较适用；再者，如将随机变量统一认为是服从正态分布，也与实际情况有差异，对计算结果的准确性无疑会产生一定影响。为此，将综合考虑Ｒ．０．Ｆｏｓｃｈｉ与文献［４］的研究成果，并结合大风与覆冰的简单模型对塔架的结构功能函数进行扩充。２．２考虑覆冰因素的塔架功能函数模型ｉ之间的简单关系模型［ｓ］，即ｉ≈鬲‰・［Ｐ２＋（ｏ．２４１２∥４ｓ）２］ｌ／２刘和云等［５］提出了瞬时风速ｕ与线路覆冰厚度（２４）式中：Ｐ为单位时间降水量（ｍｍ／ｈ）；口为瞬时风速（ｍ／ｓ）；Ｎ为冻雨时长（ｈ）；ｄ为冷却水滴直径（肛ｍ）。上式表示在冻雨过程中，由于风的作用使降水逐渐结为线路覆冰，最终覆冰厚度积累到ｉ。近似计算时利用研究地点的平均冻雨降水率和时长，以极大风速统计数据作为风速随机变量，简化式（２４）为覆冰厚度ｉ关于风速口的函数，利用ｉ与口的关系，进而将由覆冰ｉ对原塔架功能函数的影响因子用风速口的函数表示。即改写功能函数的形式如下：Ｇ（Ｒ，口）＝Ｒ一［Ｋ凹ｒＴ７（口）Ｓ＋口ｖ（ｙ７（口）＋Ｄ）・Ｓｐ＋ｄ＋ＫｒＣｓａＰＰ（ｕ）］（２５）式中：Ｒ为塔架承载能力随机变量（ｋＮ），服从对数正态分布；口为风速随机变量（ｍ／ｓ），服从极限Ｉ型分布；Ｔ７（可）与Ｖ７（ｕ）分别为利用式（２４）和当地冻雨气象参数改写Ｒ．Ｏ．Ｆｏｓｃｈｉ模型中Ｔ（ｖ，ｉ）和Ｖ（ｉ）的相应结果，单位均为ｋＮ／ｍ；其他参数意义同式（２２）。２．３电塔可靠性指标的计算上述构建的塔架功能函数涉及两个随机变量，服从对数正态分布的Ｒ和服从极限Ｉ型分布的口。根据第１．２．１节所述，利用在改进一次二阶矩法（设计验算点法）基础上扩展的映射变换法，既可以很好地处理功能函数非线性问题，也兼顾了随机变量非正态的情况。在利用映射变换法进行计算时，针对功能函数中非正态随机变量的具体分布分别予以如下处理：（１）Ｒ服从对数正态分布，并已知均值和方差分别为Ｅ（Ｒ）和Ｄ（Ｒ）。根据对数正态随机变量的性质，则ｌｎＲ服从正态分布。将Ｒ映射为标准正态随机变量ｙ的函数。引，即万方数据Ｒ＝ｅｘｐ（ｐ（１ｎＲ）＋Ｙ矗（１ｎＲ））（２６）式中：卢ｃ㈤为ＩｎＲ的期望值；仃（㈣为ＩｎＲ的标准差。卢ｃ雌，和鳓氓，的计算式分别为…乩ｃ杀舞，∽，盯（Ｉ根）＝￣／ｌｎ（１十Ｄ（Ｒ））（２８）对应式（２０）中求偏导，第二个因式为孤８Ｒ＝Ｒ刚哟（２９）（２）７．／服从极限Ｉ型分布，并已知均值和方差分别为Ｅ（口）和Ｄ（ｕ）。将口映射为标准正态随机变量ｚ的函数‘引，有口一Ｕ一！翌｛二！旦匦兰业口２一——￣＿———＝————‘（３０）Ｌ３ＵＪ其中：“一Ｅ（ｕ）－－０．４５棚灭万（３１）口一罢．—兰（３２）４６ＣＤ（ｕ）对应式（２０）中求偏导，第二个因式为８ｚ垫一—』二竺鱼）＿订ａｃＰ（ｚ）ｌｎ［ｑｂ（ｚ）］（３３）…～则原结构功能函数Ｇ（Ｒ，础）可简化为Ｇ７（ｙ，２）＝ｅｘｐ（ｐ（１氓）＋％（Ｉ删）一ｆ（ｕ一！翌｛二！翌［虫堕３１）可以根据第１．２．２节所述的利用映射变换法迭可靠度和失效概率。例利用文献［７３中关于我国某城市平均风速极大ｍ／ｓ的极限Ｉ型分布，干字形ｍ，结构承载能力Ｒ服从均值为７５０ｋＮ、ｋＮ的对数正态分布［３］，采用映射变ｍ、４００ｍ、５００ｍ表１年度输电塔架失效概率ＴａｂｌｅｌＹｅａｒｌｙｆａｉｌｕｒｅｐｒｏｂａｂｉｌｉｔｙｏｆｔｒａｎｓｍｉｓｓｉｏｎｔｏｗｅｒ式中：厂（ｕ）表示原功能函数关于风速口的函数部分；Ｇ７（ｙ，２）为服从标准正态分布且相互的随机变量ｙ和ｚ的函数；其他符号意义同前。代计算可靠性指标ｐ，进而利用公式（４）和（５）计算３算值的统计结果，如月极大风速口服从均值为２６．３１７ｍ／ｓ、标准差为２．９２６塔架高４０换法进行试算。表１为档距在３００标准差为０．１第３期李致宇：大风覆冰条件下输电杆塔可靠性模型的研究６３的，且各个影响因素及其相互之间统计性质的准确０．００２５０．００２０度直接依赖于观测资料收集的详实程度。当一个结构功能函数确定后，需要根据其特点选取适当的可靠性计算方法进行处理，才能保证分析结果满足所需精度要求。（２）通过分析比较几种输电杆塔功能函数的研１月２月３月４月５月６月７月８月９月１０月１１月１２月月份鍪ｏ．ｏｏｌ絮０．００１抵０．０００５０５Ｏ究成果，提出了一种考虑大风和覆冰因素的杆塔功能函数模型，并根据功能函数特点并结合可靠性指标计算理论选取适当方法进行了试算，其结果可为电力系统安全分析和风险评估提供依据。参考文献：图１各月塔架失效概率Ｆｉｇ．１Ｍｏｎｔｈｌｙｆａｉｌｕｒｅｐｒｏｂａｂｉｌｉｔｙｏｆｔｒａｎｓｍｉｓｓｉｏｎｔｏｗｅｒ时，全年发生因大风和覆冰导致塔架失效的概率，图１为４００ｍ档距下各月发生因大风和覆冰导致塔架发生失效的概率。表１反映出随着塔架档距的增加，即塔架所连接导线长度越长，失效概率越大；图１直观地反映出１２月至次年２月因大风和覆冰造成的塔架失效率偏高，这也正是大风寒潮与冻雨经常发生的月份。［１］赵国藩，金伟良，贡金鑫．结构可靠度理论［Ｍ］．北京：中国建筑工业出版社，２０００：ｚ２—５２．［２］张明．结构可靠度分析：方法与程序［Ｍ］．北京：科学出版社，２００９：３０—４５．［３］Ｆｏｓｅｈｉ，Ｒ．０．Ｒｅｌｉａｂｉｌｉｔｙｔｈｅｏｒｙａｎｄｏｆｐｏｗｅｒｃｏｍｐｏｎｅｎｔｓａｎｄａｐｐｌｉｃａｔｉｏｎｓｔｏｒｉｓｋａｎａｌｙｓｉｓｓｙｓｔｅｍｓ［Ｊ］．ＩｎｔｅｒｎａｔｉｏｎａｌＪｏｕｒｎａｌｏｆＥｌｅｃｔｒｉｃａｌＰｏＴｗｅｒ＆ＥｎｅｒｇｙＳｙｓｔｅｍｓ．２００４，２６（４）：２４９－－２５６．［４］黄超．典型城市生命线系统风险评估方法研究［Ｒ］．北京：清华大学，硕士学位论文，２００９．４结论［５］刘和云，周迪．导线雨淞覆冰预测简单模型研究［刀．中国电机工程学报，２００１，２１（４）ｔ１８—２２．［６］李云贯，赵国藩．广义随机空间内的一次可靠度分析方法［刀．大（１）功能函数可以准确地描述一个结构的功能运行状态，对于一个复杂结构安全可靠度的研究，若要考虑所有影响结构安全的因素，通常是比较困难连理工大学学报，１９９３，３３（增刊１）。１—５．［７］胡章文，陈国华，李其明．暴风雪下塔设备应力分布参数的确定［门．压力容器，２００９，２６（１）：３Ｚ～３４．（上接第５８页）［４３国家环境保护总局．地表水环境质量评价技术规范［ｓ］．［５］胡远志．关乎生命的水，到底该不该一票定优劣［ＥＢ／ＯＬ］．（２００９—０７－２７）．ｈｔｔｐ：／／ｉｎｆｏ．ｗａｔｅｒ．ｈｅ１３３２１４００４７．ｈｔｍｌ．３６０．ｃｏｒｎ／２００９０７２７／ＯＬ］．（２００９－０６—１９）．ｈｔｔｐ；／／ｗｗｗ．ｃｅ．ｅｎ／ｅｙｓｃ／ｓｐ／ｉｎｆｏ／２００９０６１９／１９２８９１６８．ｈｔｍｌ．［１１］周倩倩．农夫山泉开始攻造大势千岛湖寻源只是一场秀？［ＥＢ／ＯＬ］．（２００９—０７－２９）．ｈｔｔｐ：／／ｔｖ３１５．ｚｊ０１．ｃｏｒｎ．ｃｎ／ｐｌａｙ．ａｓｐｘ？ｐｉｄ＝２０７０ｂｄｄ２—５ｂａｅ．４４８２一ａｇｄｆ－ａ７ｃｄｂ２７ｅ８６ｅｅ＆ｔｙｐｅ一３．［６３李玉英，胡兰群，贾红霞。等．南水北调中线工程水源区粪大肠菌群和异养细菌丰度与分布［刀．南阳师范学院学报，２００６，５（３）：５３—５６．［１２］崔保山，杨志峰．湿地学［Ｍ］．北京；北京师范大学出版社，２００６．［１３］陈永根，李香华，胡志新，等．中国湖泊冬季水一气界面Ｃ０２通量［Ｊ］．生态环境。２００６，１５（４）：６６５—６６９．［１４］尹海龙，徐租信．我国单因子水质评价方法改进探讨［Ｊ］．净水技术，２００８，２７（２）；１—３．通讯作者：罗岳平（１９７１一），男，高级工程师，主要从事环境监测管理与环境质量综合评价工作。Ｅ－ｍａｉｌ：ｓｈａｎｅｕｎｙｅｌａｏ＠ｔｏｎＬｏＤｎｌ［７３邹志勇，黄颖媛．简析漓江地表水粪大肠菌群状况及综合整治［Ｊ］．农业环境与发展，２００７，２４（１）：７８—７９．［８］林美爱．闽江流域福州段粪大肠菌群分布的初步研究［Ｊ］．海峡科学，２００８，（９）：５９—６０．［９３中华人民共和国环境保护部．２００７年中国环境质量报告［Ｍ］．北京：中国环境科学出版社，２００８．［１０］商蕾．千岛湖水质惹争议千岛湖的水到底有没有问题？［ＥＢ／万方数据
