福清市雷暴日分布趋势研究及其在雷电灾害风险评估中的应用

Ｊｏｕｒｎａｌ　ｏｆ　Ａｇｒｉｃｕｌｔｕｒａｌ　Ｃａｔａｓｔｒｏｐｈｏｌｏｇｙ　２０１３，Ｖｏｌ　３。Ｎｏ　１：３４—３６，４４　ｎｇ　Ｔｈｕｎ－　福清市雷暴日分布趋势研究及其在雷电灾害　Ｒｅｓｅａｒｃｈ　ｏｆ　Ｆｕｑｉｄｅｒｓｔｒｏｍ　Ｄａｙ’Ｓ　Ｄｉｓｔｒｉｂｕ－　风险评估中的应用　ｔｉｏｎ　ａｎｄ　Ｉｔｓ　Ａｐｐｌｉｃａｔｉｏｎ　ｉｎ　俞成标　，张烨方　，黄榕城。，冯真祯　１．福清市气象局，福建福清３５０００１；２．福建省防雷中心，福建福州３５０００１；３．福建　省气象局，福建福州３５０００１　Ｌｉｇｈｔｎｉｎｇ　Ｄｉｓａｓｔｅｒ　Ｒｉｓｋ　Ａｓｓｅｓｓｍｅｎｔ　ＹＵ　Ｃｈｅｎｇ－ｂｉａｏ　ｅｔ　ａｌ（Ｆｕｑｉｎｇ　Ｍｅｔｅｏ・　摘要以福清市１９６１—２００７年总雷暴日观测数据为样本，对福清市雷暴目的极端值和短　ｒｏｌｏｇｉｃａｌ　ａｄｍｉｎｉｓｔｒａｔｉｏｎ，Ｆｕｑｉｎｇ，Ｆｕｊｉａｎ　３５０００１）　期预测值进行研究与应用。采用频数分析的方法拟舍得到福清市年总雷暴日的概率密度、　分布函数，由此计算得到福清市年总雷暴日的极端数值，并对福清市年总雷暴日观测数据　进行趋势分析及多元回归预测，为雷电灾害风险计算中雷暴日参数的选取提供了参考。　关键词雷暴日；极端值；趋势分析　Ａｂｓｔｒａｃｔ　Ｂａｓｅｄ　ｏｎ　ｔｈｕｎｄｅｒｓｔｒｏｍ　ｄａｙ　（ｆｒｏｍ　１９６１　ｔｏ　２０１０）ｏｆ　Ｆｎｑｉｎｇ　Ｃｉｔｙ，ｔｈｅ　ｅｘｔｒｅｍｅ　ａｎｄ　ｓｈｏｒｔ—ｔｅｒｍ　ｐｒｅｄｉｃｔｉｏｎ　ｖａｌｕｅ　中图分类号Ｐ４２９　文献标识码Ａ　文章编号２０９５—３３０５（２０１３）０１—０３４—０３　ｏｆ　Ｆｕｑｉｎｇ　Ｃｉｔｙ　ｗｅｒｅ　ｄｉｓｃｕｓｓｅｄ．Ｔｈｅ　ｐｒｏｂ—　ａｂｉｌｉｔｙ　ｄｅｎｓｉｙ　ｔａｎｄ　ｄｉｓｔｒｉｂｕｔｉｏｎ　ｆｕｎｃｔｉｏｎ　ｏｆ　雷电灾害是联合国公布的十大自　然灾害之一，我国是一个发展中国家，　不可能对所有的建设项目都进行全面　市雷暴日的初、终日及绝对变率，年、　月分布及变化趋势，以及雷暴初、终日　的候平均气温阈值等方面进行了分　Ｆｕｑｉｎｇ’Ｓ　ｔｈｕｎｄｅｒｓｔｒｏｍ　ｄａｙ　ｗｅｒｅ　ｉｒｅｄ　ｆｂｙ　ｔｈｅ　ｍｅｔｈｏｄ　ｏｆ　ｒｅｑｕｅｎｃｙ　ａｎａｌｆｙｓｉｓ，ＳＯ　ｇｏｔ　ｉｔｓ　ｅｘｔｒｅｍｅ　ｖａｌｕｅ．Ｔｈｅ　ｔｒｅｎｄ　ａｎｄ　ｍｕｌｔｉｐｌｅ　的防雷工程设计，因此对重大、重点项　目进行雷电灾害风险评估是保障建设　项目雷电安全、社会和谐稳定的重要　析　”。笔者主要从雷电灾害风险评估　应用的角度，以福清市为例，在对区域　雷暴日常规数理统计分析的基础上，　对雷暴日的极端值、短期可能值等不　ｒｅｇｒｅｓｓｉｏｎ　ｐｒｅｄｉｃｔｉｏｎ　ｏｆ　Ｆｕｑｉｎｇ’Ｓ　ｔｈｕｎ—　ｄｅｒｓｔｒｏｍ　ｄａｙ　ｗａｓ　ａｎａｌｙｚｅｄ，ｆｏｒ　ｇｕｉｄｉｎｇ　ｌｉｇｈｍｉｎｇ　ｄｉｓａｓｔｅｒ　ｒｉｓｋ　ａｓｓｅｓｓｍｅｎｔ．　技术内容。在雷电灾害风险评估技术　服务中，评估项目所在地的年雷暴日　是反映地区雷电活动规律的重要参　数。目前在项目雷电灾害风险计算上　使用的雷暴日主要来源于一些标准规　范的附录给出的参考数据，或者由当　地气象部门提供的多年年雷暴日平均　值，这些数值具有一定的代表性，可以　定性指标进行研究，以指导具有精确　雷电灾害风险评估需求的项目。　１　福清市雷暴日分布规律研究与应用　１．１　分布规律探索与假设检验　以福清市１９６１—２０１０年的年总　Ｋｅｙ　ｗｏｒｄｓ　Ｔｈｕｎｄｅｒｓｔｒｏｍ　ｄａｙ；Ｅｘ—　ｔｒｅｍｅ　ｖａｌｕｅ；Ｔｒｅｎｄ　ａｎａｌｙｓｉｓ　雷暴日观测数据为样本，对福清市雷　暴日的概率分布规律进行研究。首先　对雷暴日数据进行正态性分布探索分　析与假设检验，如果福清市雷暴日服　从正态分布，则采用正态分布函数进　行规律分析；如果不服从正态分布，则　采用累积频数分布表进行累积概率函　应用于一般项目防雷工程的评估计　算，但对于一些重大、重点或对评估精　度要求较高的项目而言，就需要对区　域的雷暴日分布规律与趋势进行一定　的分析、研究，用更加科学的数据来指　导防雷工程、评估与检测。　目前全国范围内对雷暴日的研究　主要集中在对雷暴日的基本统计分　数拟合，并运用该函数进行规律分析。　首先绘制福清市１９６１—２Ｏ１０年　的年总雷暴日观测数据的正态ｐ－Ｐ　析、区域雷暴日的年季变化规律、雷暴　日的时空分布特征等方向０－６３，如易燕　明等对近５Ｏ年广东省雷暴、闪电时空　变化特征进行研究Ｅ７咽；郑燕等对莱芜　图、０一Ｑ图，对雷暴日的分布类型进行　判断（图１、２）。由图１、２可看出，通过　ｐ－ｐ、０—０检验，福清市雷暴日数据的　基金项目福建省气象局基层科技专项（２０１ｌｑｌ　ｏ　作者简介俞成标（１９５９一），男，福建福清人，工程师，从事雷电灾害防御与管理研究，Ｅ－ｍａｉｌ：２２８５３２１４８＠ｑｑ．ｃｏｍ。　收稿日期３４　２０１２—１２—２３　分位数　分啦数　图１福清市雷暴日正态Ｐ—Ｐ检验结果　图２福清市雷暴日正态Ｑ—Ｑ检验结果　ｐ－Ｐ、Ｑ－Ｑ与计算点基本处在同一条直　算过程如下。　线上，可认为福清市的雷暴日基本服　假设：　从正态分布。为了证明该假设，对福清　日０：Ｘ－Ｎ（＃，　。），　１，Ｘ２，…，‰为Ｘ样　市雷暴日是否服从正态分布进行假设　本，则：　检验，采用偏峰度的检验方法，具体计　偏度ｇ】　【１）　偏度方差Ｄ（舒　６ｎ　（ｎ　－＿１）　）（　（２）　峰度　而　ｎ（ｎ　＋１）　Ｅ（Ｘｉ　－￣）４　（３）　峰度方差Ｄ（　）＝　２　４ｎ　（ｎ　－１）　ｚ｜＋丽　（４　１　式中，五为样本均值；ｎ为样本　度　＝Ｏ．１　７９，Ｄ（ｇ２）　０．６６２　数目。　Ｊ偏度＝０．３８６　２ｆ≤０．６５９　７，通过　在假设风成立的前提下，近似有　预检　ｇｌＮＮ（Ｏ，Ｄ（ｇ１）），　一Ｎ（Ｏ，Ｄ（　））给定显著　Ｉ峰度：０．９０３　６ｆ≤１．２９７　３，通过　性水平　，依据Ｐ｛Ｕ＞ｕ　｝：　２，求得临　预检　界值“　。如果Ｉ　ｇ１　ｌ＜ｕ　×、／Ｄ　ｅ　ｇ　），　通过预检，接受总体服从正态分布　的假设。即在显著性水平ａ＝Ｏ．０５的前　且Ｉ　Ｉ＜　Ｘ、　，则接受　，反　提下，认为福清市年雷暴日的分布规律　之，则拒绝假设。　服从：　＝３７．０４，　＝８．３９３的正态分布。即　显著性水平Ｏｒ．＝０．０５，临界值　应有福清市年雷暴目的概率密度函　＝１．９５９　９６４，样本长度ｎ＝５０　数式：　偏度ｇｌ＝０．３８６　２，Ｄ（ｇ１）＝Ｏ．３３７，峰　ｐ（　）　１　面×ｅ　。一（　＝　Ｚ。：Ｑ　＝０．　０４７　５×ｅ（．　Ｔ￣”－３　７．鹋０４Ｚ）２＿　（５）　分布函数式为：　ｚ：唑　Ｐ（　）＝ｆ　ｐ（ｔ）ｄｔ＝ｌ　０．０４７　５ｘ　ｅ¨　ｄｔ　（６）　０　式中，　为雷暴日。　１．２年雷暴日概率密度、分布函数在　解析式来计算项目极端的雷暴日发生　雷电灾害风险评估中的应用　概率的计算。　在雷电灾害风险评估中，计算建　假设项目的建设性质为一般民用　设项目的预计年雷击次数是重要的评　建筑，设计使用年限为７０年。因此，按　估内容之一，应用计算得到的福清市　７０年内可能发生１次某极端雷暴目值　１９６１—２Ｏ１Ｏ年雷暴目正态分布函数　为累计概率值（即１／７０＝０．０１４　３），按　农业灾害研究２０１３。３（Ｏ１）：３４—３６，４４　～　３７．０４，８．３９３　）的正态分布有：　ｐ（Ｔ，Ｉ－３　７　．　０４一）＞０．９９２　８５　（７）　Ｈ查正态分布分位表得：　ＴｊＲ　ｑ　－　３７．０４　２：　‘．　４５０　（８）…　　求得Ｔ￣＝５７．６Ｏ　ｄ，即在现有统计数　据的基础上，可认为项目所在区域自项　目建成至项目设计使用年限内可能遇　到最大、最不利的雷暴日为５７．０６　ｄ，实　际设计过程中使用该极端计算值进行　建筑物防雷类别设计及其他防雷装置　参数的计算，相比于目前使用的平均值　３７　ｄ，通过统计分析得到的雷暴日数值　更能反映项目的雷电环境特点。　２福清市雷暴日趋势分析与应用　２．１基本趋势分析　笔者计算了福清市ｌ９７０—２０１０　年累计雷暴日平均值的数据并绘制成　折线图（图３），从福清市ｌ９７０—２０１Ｏ　年的总雷暴日数据折线可看出，福清市　年总雷暴日的分布规律并不明显，整体　上呈现出下降的波动趋势。福清市雷暴　曰平均值在持续下降后大概于１９８５　年以后逐渐趋于稳定，在３７￣３８　ｄ之　间，对福清市１９６１—２０１０年的年总雷　暴日数据进行线性拟合得到的直线方　程（图４中的直线）也证明福清市的雷　暴日呈现出下降的趋势，在福清市雷暴　日的波动趋势上，季节性的周期特征并　不明显。　２．２福清市雷暴日短期预测探索　通过对福清市雷暴日基本趋势的　分析，笔者认为福清市雷暴日在整体上　并不呈现出稳定的周期振动过程，即在　时间序列上的表现并不明显。通过采用　ＳＰＳＳＩ９．０程序对福清市雷暴日的　ＡＲＩＭＡ时间序列建模分析也可以得出　福清市的雷暴日分布规律不适合采用　时间序列的分析方法，需要根据样本数　据的特点自行建模。　该研究拟通过对福清市雷暴日的　分布规律进行更深入的分析，探求影响　雷暴日分布的可能组成因子，通过对各　个因子进行多元的线性拟合，采取多元　拟合的方式对福清市雷暴目的分布规　律进行研究。　首先，由于福清市雷暴曰在趋势上　３５　Ｊｏｕｒｎａｌ　ｏｆＡｇｒｉｃｕｌｔｕｒａｌ　Ｃａｔａｓｔｒｏｐｈｏｌｏｇｙ　２０１３，Ｖｏｌ　３。Ｎｏ　１：３４－３６，４４　口　岫　图３　１９７０—２０１０年福清市雷暴日平均值　叮　帕　年份　图４　１　９７Ｏ一２０１　０年福清市雷暴日真实值与拟合预测值　呈现出下降的趋势，因此根据最小二乘　据一阶差分因子的１　Ｎ４滞后序列加入　法的原理，图４中根据福清市已知雷　拟合方程中）。　暴日数据拟合得到的线性方程作为其　综上所述，拟建立如下拟合方程：　分布趋势的最小预测残差值可以作为　预测值Ｆ＝ｏｔ　×［线性拟合值］　福清市雷暴日下降趋势的一个代表因　＋ｄ２ｘ［平均值］＋ｄ　×［一阶差分０滞后　子。其次，从稳定性的角度分析，福清市　序列］＋　×［一阶差分１滞后序列］　雷暴日的累计平均值可以作为其基本　＋　５×［一阶差分２滞后序列］＋　×［一　数值波动中心的代表，也选择其作为一　阶差分３滞后序列］　（９）　个拟合因子。在对福清市雷暴日波动方　为了验证所设想的方程与实际是　面的表达上，笔者借鉴了一般时间序列　否相符，以福清市１９６１—２００６年的年　预测方法里的差分算法，利用对福清市　总雷暴日数据作为假设已知样本数据，　雷暴日的数据进行各个不同阶次的差　按照上述拟合方程进行最小二乘法回　分计算，用一阶差分数据的不同滞后序　归，得到预测方程后计算出２００７—　列作为另一类的拟合因子（通过多次拟　２０１０年的预测值，用其与实际值进行　合与比较，最后选择了福清市雷暴目数　对比来验证设想。采用Ｅｘｃｅｌ－２００３自　表１预测方程最小二乘回归计算结果　带的最小二乘回归程序计算回归系数，　具体计算与检验结果见表ｌ。　其判定系数Ｒ　＝０．９３７　６，综合得最　后的预测方程．　预测值产一２０．０７５＋３．１　５３　ｘ［线性　拟合值卜１．８３７×［平均值］一０．７８３×　［一阶差分０滞后序列卜Ｏ．５６８　ｘ［一阶　差分１滞后序列］－０．３７７×［一阶差分２　滞后序列卜Ｏ．１　７７×［一阶差分３滞后　序列］　（１０）　根据上述预测方程，笔者计算了　延后３年的福清市年总雷暴日数据并　与实际值进行了比较：　２００８年实际值３７，预测值　３６．３９７，预测偏差０．６０３；２００９年实　际值２７，预测值２９．１８２，预测偏差　一２．１８２；２０１０年实际值２６，预测值　３１．２ｌ　５，预测偏差一５－２１　５。　从预测数值与实际值比较结果来　看，笔者所采用根据样本趋势、规律总　结出与数据规律有关的统计量，再对　这些参数进行最小二乘回归的建模方　法，在实际中得到了很好的应用，特别　是对一些规律比较复杂、隐晦，数据关　系不明显，不适合使用常规预测、建模　方法的样本，采用该思想可能得到合　适的预测或拟合方程。　３结论　以福清市１９６１—２０ｌ０年年总雷　暴日数据为基础，对福清市年总雷暴　日的分布规律及短期趋势进行了讨论　与研究。结果认为，福清市年总雷暴日　服从置信度９５％的正态分布，其短期趋　势上的预测方程可通过对年总雷暴日　的线性拟合值、平均值、一阶差分０～３　滞后序列的回归计算拟合获得，并且　该拟合方程适用于福清市年总雷暴日　数据的短期预测。论文还提出对一些　分布规律不明显、数据关系复杂的样　本，可采用多元函数拟合的方式，通过　对样本的分析，获取可能影响样本数　据分布的统计量，采用最小二乘法进　行估计来获取样本的趋势方程。　Ｊｏｕｒｎａｌ　ｏｆ　Ａｇｒｉｃｕｌｔｕｒａｌ　Ｃａｔａｓｔｒｏｐｈｏｌｏｇｙ　２０１３，Ｖｏｌ　３。Ｎｏ　１：４１－４４　参考文献　…１虞统空气质量日报中的空气污染指数【Ｉ１．　城市管理与科技，２０００，２（１）：２３—２６．　【２】ＭＡＭＴＩＭＩＮ　Ｂ，ＭＥＩＸＮＥＲ　Ｆ　Ｘ．Ａｉｒ　ｐｏＨｕｆｉｏｎ　ａｎｄ　ｍｅｔｅｏｒｏｌｏｇｉｃａｌ　ｐｒｏｃｅｓｓｅｓ　ｉｎ　［８］李景林，郑玉萍，刘增强乌鲁木齐市低　空温度层结与采暖期大气污染的关系　干旱区地理，２００７，３０（４）：５ｔ９—５２５　空气质量的影响Ⅱ】环境科学学报，２０１　１，３１　（６）：１１６８—１１７４　［１８］唐燕秋，陈佳，熊强，等重庆市多年空　气污染指数分析及大气污染控制对策　四川环境，２００５，２４（６）：８０—８２　［９］张艳燕，孟凡，何友江，等鸟鲁木齐市冬　ｔｈｅ　ｇｒｏｗｉｎｇ　ｄｒｙｌａｎｄ　ｃｉｗ　ｏｆ　Ｕｒｕｍｑｉ（Ｘｉｎ—　季典型污染事件气象过程分析　环境　科学研究，２０１２，２５（１）：１０—１７．　ｊｉａｎｇ，Ｃｈｉｎａ）田．Ｓｃｉｅｎｃｅ　ｏｆ　ｔｈｅ　Ｔｏｔａｌ　Ｅｎｖｉ—　ｒｏｎｍｅｎｔ，２０１　１，４０９：１２７７－１２９０．　【１９］胡琳，林杨，何晓嫒基于ＡＰＩ方法的西　［１０１吴彦，王健，刘晖，等．乌鲁木齐大气污　染物的空间分布及地面风场效应Ｕ１ｌ中　国沙漠，２００８，２８（５）：９８６—９９１　安城市大气环境质量评价［『］．陕西气象，　２００７（１）：１８—２０．　［３］ＤＯＮＧ　Ｚ　Ｗ，ＬＩ　Ｚ　Ｑ，ＥＤＷＡＲＤＳ　Ｒ，ｅｔ　ａｌ　Ｔｅｍｐｏｒａｌ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｏｆ　ｍｉｎｅｒａｌ　ｄｕｓｔ　ｐａｒｔｉｃｌｅｓ　ｉｎ　ｐｒｅｃｉｐｉｔａｔｉｏｎ　ｏｆ　Ｕｒｕｍｑｉ　Ｒｉｖｅｒ　Ｖａｌｌｅｙ　ｉｎ　Ｔｉａｎ　Ｓｈａｎ，Ｃｈｉｎａ：Ａ　ｃｏｍｐａｒｉｓｏｎ　［２０】李文杰，张时煌，高庆先，等．京津石三　市空气污染指数（ＡＰＩ）的时空分布特　［１１１杨静，武疆艳，李霞，等鸟鲁木齐冬季　大气边界层结构特征及其对大气污染的　影响【Ｉ１．干旱区研究，２０１１，２８（４）：７１７　７２３．　征及其与气象要素的关系　资源科学，　２０１２，３４（８）：１３９２—１４００　［２１］Ｌ１　Ｊ，ＺＨＵＡＮＧ　Ｇ　Ｓ，ＨＵＡＮＧ　Ｋ，ｅｔ　ａｌ　Ｃｈａｒａｃｔｅｒｉｓｉｃｓｔ　ａｎｄ　ｓｏｕｒｃｅｓ　ｏｆ　ａｉｒ－ｂｏｒｎｅ　ｐａｒｔｉｃｕｌａｔｅ　ｉｎ　Ｕｒｕｍｑｉ，Ｃｈｉｎａ，ｔｈｅ　ｕｐ—　ｏｆａｌｐｉｎｅ　ｓｉｔｅ　ａｎｄ　ｍｒａｌ　ｓｉｔｅ　ＵＪ．Ａｔｍｏｓｐｈｅｒｉｃ　Ｒｅｓｅａｒｃｈ，２０１１，１０１：２９４—３０６　ｆ１２］刘增强，郑玉萍，李景林，等．乌鲁木齐　市低空大气逆温特征分析卟干旱区地　理，２００７，３０（３）：３５１—３５６　［４】ＣＨＡＮ　Ｃ　Ｋ，ＹＡＯ　Ｘ　Ｈ．Ａｉｒ　ｐｏｌｌｕｔｉｏｎ　ｉｎ　ｍｅｇａ　ｃｉｉｔｅｓ　ｉｎ　Ｃｈｉｎａ［Ｊ］．Ａｔｍｏｓｐｈｅｒｉｃ　Ｅｎｖｉ—　ｒｏｎｍｅｎｔ，２００８，４２：１－４２　ｆ１３］潘虹梅，李凤全，王俊荆，等基于ＡＰＩ　方法的城市大气污染评价Ⅱ］环境科学　与管理，２００８，３３（２）：１７８—１８４．　【１４］段玉森，魏海萍，伏晴艳，等．中国环保　重点城市ＡＰＩ指数的时空模态区域分异　ｓｔｒｅａｍ　ａｒｅａ　ｏｆ　Ａｓｉａ　ｄｕｓ　ｔ卟Ａｔｍｏｓｐｈｅｒｉｃ　Ｅｎｖｉｒｏｎｍｅｎｔ，２００８，４２：７７６—７８７．　【５】Ｊ．ｍ．ＢＡＬＤＡＳＡＮＯ，Ｅ　ＶＡＬＥＲＡ，Ｐ　ＪＩＭＥ　ＮＥＺＡｉｒ　ｑｕａｌｉｗ　ｄａｔａ　ｆｒｏｍ　ｌａｒｇｅ　ｃｉｉｔｅｓ　Ｕ１．　Ｔｈｅ　Ｓｃｉｅｎｃｅ　ｏｆ　ｔｈｅ　Ｔｏｔａｌ　Ｅｎｖｉｒｏｎｍｅｎｔ，　２００３，３０７：１４１—１６５　［２２】魏毅坞鲁木齐市空气污染成因及防治　对策研究ｏ］干旱区资源与环境，２０１０，２４　（９）：６８—７１．　【６］ｗＡＮＧ　Ｌ　Ｌ，ｗＡＮＧ　Ｙ　Ｓ，ＳＵＮ　Ｙ，ｅｔ　ａｌ　Ｕｓｉｎｇ　Ｓｙｎｏｐｄｃ　Ｃｌａｓｓｉｉｆｃａｔｉｏｎ　ａｎｄ　Ｔｒａｊｅｃ—　ｔｏｒｙ　Ａｎａｌｙｓｉｓ　ｔｏ　Ａｓｓｅｓｓ　Ａｉｒ　ＱｕａｌｉＷ　ｄｕｒｉｎｇ　ｔｈｅ　Ｗｉｎｔｅｒ　Ｈｅａｔｉｎｇ　Ｐｅｒｉｏｄ　ｉｎ　Ｕｒｕｍｑｉ，Ｃｈ　【Ｉ１．环境科学学报，２００８，２８（２）：３８４—３９１．　ｆ２３］魏毅，加旭辉，魏勇．乌鲁木齐市大气污　染总量控制和降低机制研究【『Ｉ．干旱区　资源与环境，２０１１，２５（２）：１　一１２３．　（１５１李向阳，丁晓妹，高宏，等．中国北方典　型城市ＡＰＩ特征分析【『１．干旱区资源与　环境，２０１１，２５（３）：９６　１０１．　［２４］刘新春，钟玉婷，何清，等．乌鲁木齐及　周边城市空气质量变化特征及影响因　ｉｎａⅡ］．Ａｄｖａｎｃｅｓ　Ｉｎ　Ａｔｍｏｓｐｈｅｒｉｃ　Ｓｃｉｅｎｃｅｓ，　［１６１廉丽姝，高军靖，束炯城市大气污染特　征及其与气象因子的关系——以济　２０１２，２９（２）：３０７—３１９．　素分析卟沙漠与绿洲气象，２０１０，４（４）：　１２—１７　【７］杨静，李霞，李秦，等乌鲁木齐近３０ａ大　气稳定度和混合层高度变化特征及与空　气污染的关系Ｄ】．干旱区地理，２０１１，３４（５）　：７４７－７５２　南、青岛市为例田环境污染与防治，　２０１１，３３（５）：２２—２６　［１７１马雁军，刘宁微，王扬锋，等沈阳及周　边城市大气细粒子的分布特征及其对　责任编辑：张杨林　（上接第３６页）　参考文献　【１】黎梓华，梁和生，张宇，等．桂林市雷电灾　害及风险区划研究口１．农业灾害研究，　２０１ｔ，２（１）：８３—８６．　［５】申元，王磊，马御棠，等．海拔高度对云南　的气候特征　安徽农业科学，２００７，３５　（２７）：８６０７—８６０８　某地雷电参数的影响　２０１２．３３（４）：３５—３７．　电力建设，　［１０】刘顺忠数理统计理论、方法、应用和软　［２】李丹，单振宇，韩卫东．丹东地区雷暴的　气候特征研究ｌＪ　ｌ＿宁夏农林科技，２０１１，　５２（３）：４５—４７　【６】戴小笠，段小凤，孙少华，等．宁夏红枣干　热风气象等级预报Ⅱ］．农业灾害研究，２０１１，　１（２）：７４—７６　件计算【Ｍ］．武汉：华中科技大学出版　社，２００５：１４９—１５５　［１１］李永福，司马文霞，陈林，等基于雷电　定位数据的雷电流参数随海拔变化规　律Ｕ］．高电压技术，２０１１，３７（７）：１６３４　【３】李政，肖稳安，李家启，等．区域海ｔ￣－ａ７度　对闪电特征影响的初步分析卟湖北大　学学｛艮自然科学版，２０１　１，３３（２）：１９７—２０１．　［４１张宇，张翼，廖慕科，等桂林雷暴时空分　布特征分析田农业灾害研究，２０１１，１＜２）：　８０－８２　【７］易燕明，杨兆礼，万齐林．广州市闪电密　度特征分析卟资源科学，２００６（１）：　１５１—１５６．　１６４１　［８】沈斐敏．安全系统工程理论与应用［Ｍ］　北京：煤炭工业出版社，２００１：２０１—２１５　责任编辑：张杨林　【９】郑燕，任学民，亓翠芸，等．莱芜市雷暴日