基于fdd的lemp在地下空间的分布特性研究.docxVIP

基于fdd的lemp在地下空间的分布特性研究 1 lemp在地下空间的分布特性 作为一种高频率的辐射源，lemp可能会干扰和破坏电子设备、电源线路和通信电缆。随着越来越多的易受攻击的电子设备和线路被使用，lemp越来越受到重视和研究。总的来说，lemp的研究方向主要包括lemp的观测和分析；lemp计算方法的研究和改进；lemp-fis-ton和设备的连接实验及其计算方法的研究；复杂的地表环境和土壤电离作用对lemp的影响。大多数设备、线路、天线和传感器配置在地面空间内，因此大部分研究集中在地面以上。随着人类经济的需求、科学技术的发展和世界极端天气的恶化，越来越多的管道、天然气和输油管、隧道、罐体、地下埋体和各种管道出现在地下空间。例如，lemp作为世界上最大、易受攻击的接收地点，对地下设施和人员的影响是显而易见的。作为地下设备和连接线的保护基础，研究lemp在地下空间的分布特性具有重要意义。 此前有关于LEMP在地下空间的研究不多,主要集中在LEMP在地下的计算方法和对于埋地线缆感应过电压的计算.就作者所知,只有文献对LEMP在大地中的分布做了专门的探讨.但它与其他已发表的研究类似的是,这些研究都是通过对比几个参考点上的LEMP变化,来推测LEMP在整个区域大略的变化趋势和分布特点.导致这种情况的一个重要原因就是过于巨大的计算和存储量.对于单个点位置上LEMP的计算,现有大部分计算方法(如sommerfle积分法及其快速等效算法、Cooray-Rubinstein(CR)公式法等)都有很好的表现.但是如果需要得出整个计算空间中每个细小单元上的LEMP,这些算法需要付出的计算时间将是不可承受的.另一方面,如果需要将空间中所有单元上的LEMP不失真的存储起来并进一步提取其特征,则对于所有已发表的算法(包括传统的FDTD方法),其数据存储量都是不可接受的.但若要得出LEMP在地下空间的分布特性,以上两方面的需求却恰恰都是现有计算方法无法避免的. 综上所述,为了满足这些需求,作者将传统的FDTD方法改造为一种实时中断处理的FDTD方法.通过这种改进的FDTD方法,整个空间每个单元网格上的LEMP和其特征参量(峰值、最大变化率、功率密度和能量密度)被计算得到并被表示在多组地下空间分布彩图中.通过分析各个参量在地下空间的分布图和传播过程,得到了LEMP在地下空间的衰减规律和分布特点.这些规律和结论将为各种地下管线及线缆的布置规划,地下设施的电磁防护设计提供重要的参考. 2 计算方法和验证 2.1 基于fdd的lemp空间分布特性 一般而言,计算LEMP的方法可以分为三类:1)基于sommerfled积分的数值等效方法;2)传统电磁场数值方法(如FDTD,MOM方法);3)CR公式及其简化算法.如上文所述,若要得出LEMP在地下空间的分布特性,整个地下空间中的每一个单元网格上的LEMP需要计算、保存后抽取其特征参量.而这样的计算量和存储量对于现有方法是不可承受的.举例来说,对于一个长3 km的闪电通道,如果需要得到其2 km范围的LEMP分布特性.当整个二维空间(轴对称空间的简化)被分为1 m2的单元网格时,共需要计算6百万个网格点上的LEMP.而就作者所知,在保证足够精确度的前提下,最快的计算单个观测点上LEMP的算法需要1.4 s左右.按此计算,完成整个空间的计算大约需要2333 h,即97 d.如果要进一步得到LEMP的空间分布特性,以文献的抽样率保存前20μs每个网格点上的LEMP波形,仅三个电磁场分量需要约550 GB的存储空间(每个单精度数占用4Byte存储空间). 传统的FDTD方法在计算整个空间LEMP方面有着别的算法无可比拟的优势.它通过求解Maxwell方程组计算整个空间中每一个网格上每一个时间步进的电磁场.但是由于闪电放电通道的跨度大,同时其产生的LEMP又具有很宽的频谱分量,这导致采用FDTD计算LEMP时的网格数量较多.并且当需要获得LEMP的空间分布特性时,为了得到精确的分布图,其网格也要求划分的较细这导致采用FDTD计算LEMP时的网格数量更多所以,如果为了提取每个网格上LEMP特征来得到空间分布特性,而将空间中每个网格上的LEMP波形信息都存储起来,存储量是十分巨大的.作者通过在FDTD程序的时间步进间加入了一种“处理窗”程序,包括中断和处理程序以及各种标志位.通过“处理窗”程序(窗的具体时间宽度视精度要求而定)得到LEMP的瞬时场值、时间导数等特征参量,然后经过对比、筛选、乘积和求和等运算得到在整个计算范围内每个网格上的特征参量.换言之就是将整个LEMP提取特征的过程分解后加入到了FDTD的不同时间步进之间,绕开了巨大的存储量,直接得到了LEMP在空间中的分布特性. 在本文中,作者在二维圆柱坐标系中采用