基于ucar数据的nmf2局部地区nmf2模型构建.docxVIP

基于ucar数据的nmf2局部地区nmf2模型构建 1 实施opagation神经网络技术 根据地震电离层可视化分析，在发生大地震之前，婚震区空间电离层v2层的最大电子密度为nmf2（2拉离层），以及电子层电子密度（各层）的异常变迁[1、2、3、4、5、6、7、8、9、10]。但在这些研究成果中,电离层异常判断的背景值是基于震前多天观测值的均值,而这种异常判断方法无法保证背景值不会受到地震活动的影响,所以其准确性值得商榷。本文将采用BP(Back Propagation)神经网络技术,利用2006—2010年UCAR公布的NmF2数据,建立局部地区电离层模型,解算电离层NmF2,提供“干净”的电离层NmF2参考值,以避免异常判断方法的弊端。 掩星观测与电离层测高仪、卫星电离层顶部探测、地基TEC观测等电离层探测手段相比,掩星观测能够提供准实时、高精度、全天候的全球电离层三维电子密度数据,而且,相关研究结果表明掩星探测数据反演的电子密度剖面与非相关散射雷达和电离层测高仪观测结果一致。目前,美国、德国、阿根廷等多个国家均有掩星计划在研或运行,而中国台湾地区与美国合作的COSMIC(Constellation Observing System for Meteorology, Ionosphere and Climate)是迄今为止最成功的掩星计划,它由6颗在800km处运行的低轨小卫星组成观测星座,在轨期间能提供每天最多约3 000个掩星探测剖面,虽然现在卫星已到使用寿命,但其在轨6年累计的观测数据能为电离层模型的建立提供丰富的电离层电子密度数据,而且正在研制的COSMIC-2系统,预计其数据观测能力将十倍于COSMIC-1。 为了研究2008年5月12日汶川7.9级地震前孕震区电离层电子密度的变化情况,我们选取UCAR公布的2006-08—2010-12月电离层F2层最大电子密度NmF2,建立以震中为中心,经纬向跨度为80°×80°(10°S～70°N,60°E～140°E)的区域电离层NmF2模型,解算震前孕震区上空各掩星NmF2探测点的模型值,通过观测值与模型值的比较,研究汶川地震前电离层是否出现异常变化。 2 区域电离层模型的不确定性 选取年积日DOY、当地时LT、经度LON、纬度LAT和F10.7太阳活动指数FLUX为网络输入,掩星观测点NmF2为网络输出,利用MATLAB神经网络工具箱构建BP网络模型,以2006-08—2010-12月考察范围内的数据集训练网络。另外,由于在磁暴发生时期,COSMIC数据误差较大,为了使模型更加精确,我们剔除磁扰时段(kp≤3)的数据。而且,在日出日落时段,电离层电子密度变化剧烈,故也将此时段的数据剔除。目前的研究显示,汶川地震前的电离层异常主要出现在当地时的下午时段,同时,为了避免过多的数据导致MATLAB运算内存溢出错误,故只选取考察区域当地时12:00LT—17:00LT的掩星探测数据来训练网络。 另外,为了确保区域电离层模型的准确性,我们采用均方根(RMS,Root mean Square Error)描述网络的整体误差。其中,Oi为观测值,Mi为模型值,N为符合筛选条件的掩星NmF2观测点的总数。经过比较各种网络参数的组合,例如:隐层层数及神经元个数、激活函数、学习速率、期望误差等,得出最优的网络参数组合,根据 RMS=1N∑ni=1(Oi?Mi)2????????????????√RΜS=1Ν∑i=1n(Οi-Μi)2 得到对应的最小网络误差为:RMS=1.1×105el/cm3。 3 两组nmf2观测点的变化情况 针对汶川地震前出现的电离层异常,目前的研究成果均基于TEC数据,而电离层NmF2与TEC有明显的正相关性,于是可以将汶川震前TEC的异常变化情况与NmF2的变化情况相比对。相关研究表明,汶川地震前第3天(5月9日)当地时下午17:00LT震中附近2 495km半径范围内,震中南部及东南部上空TEC出现明显的增强现象,幅度增大约40%～50%。图1(a)黑色虚线圈内为NmF2明显增大的区域,半径约为2 000km,两条蓝色实线之间的NmF2观测点将用来代表震中所在经度带内NmF2在不同纬度的变化情况;图1(b)表示震中附近3个NmF2观测点对应掩星事件的电子密度剖面显示的2008年5月9日12:00LT—17:00LT考察区域内的电离层NmF2观测值与由BP神经网络模型解算相应的NmF2模型值之间的偏离程度,表示为:(O-M)/M,其中,O为NmF2观测值,M为对应NmF2观测点的模型值。从图1(a)可以看出,震中南侧大约475km(95km/°×5°=475km)处,电离层NmF2出现明显增强的现象(图1(a)中的3号NmF2点),幅度约为50%,而离3号Nm