汶川ms80地震前电离层电子密度变化.docxVIP

汶川ms80地震前电离层电子密度变化 地震反演与地震勘探 云表示系统（colt-key-meteoron）由六个相同的微卫星组成。每组卫星包括gps接收机、小型电离层光度计（tip）和三段星信息仪（tri-束缚）等科学勘探仪。这六个卫星分布在大约800公里的轨道高度上，倾角为72，相邻卫星轨道平面之间的闭合角度为30。通过接收gps卫星l段的信号，科奇卫星可以计算电离层通信过程中相位的振幅。利用反演技术，可以计算波前电离层的波前电离层的异常程度。在这项研究中，许多科学家提供了地震电离层耦合机和地震电离层异常的研究方法（fenoilo等人，1995；shen；duetal，2002；pulms（2004），本文从近地表过程、电离层d效应、e层效应和磁层效应四个方面探讨了地震电离层的物理模型。hayakawa等人（2004）提出了两种地震电离层电子显微镜下的d效应、e层效应和磁层效应等方面探讨了地震电离层耦合的物理模型。freud等人（2004）根据实验室岩石实验的基础提出了离子空穴运行假设模型。然而，与地震相关的各种电离作用效应的机制目前还没有统一的解释。在本研究中，基于上述地震透射现象，利用云杉星gps传播的电子层电子密度数据，研究了川地震前后电子密度的显著变化。基于12阶次球氟子层函数的最小二乘模拟值，规划了电离层模型的电离模型，为研究地震前后电子密度的变化提供了基础依据。同时，利用主成分分析法研究了中心附近最大电子密度的各个方面的能量份额随时间的推移而变化。 1 数据分析的方法 1.1 数据分离和相关分析 主成分分析法 (principal component analysis, 简写为PCA) (Xu, 2004;Zhao, Wan, 2005) 是将多个变量通过线性变换重新组合成一组新的线性无关的变量, 同时根据分析的需求选取前几个变量反映原来变量信息的一种多元统计方法.PCA方法可以有效地分离原始数据信号, 剔除噪音和冗余.因此本文使用PCA方法分离可能与太阳活动和地磁活动相关的信号, 进而考察可能与地震活动相关的信号. 1.2 模型估计误差 本文采用12阶次球谐函数拟合电离层模型.首先, 设计网格划分, 网格大小按2°×2°分割成一定数量的网格 (如果网格内有多个即值取均值) ;然后, 基于12阶次球谐函数 (袁运斌, 欧吉坤, 2002;章红平, 2006) : 式中, anm和bnm为待估系数, Pnm(cosθ) 为缔合勒让德函数, λ为经度, θ为余纬, n=12, 共有待估系数169个, 使用最小二乘拟合法计算球谐系数;最后, 计算设计区域模型值, 再以模型值作为背景依据, 计算δ, δ= (x0-xm) /xm(x0表示观测值, xm表示模型值) , 以估算模型误差的影响;以2011年1月1—15日为例, 这段时间地磁场平静, 太阳活动水平很低, 计算得到平均误差为0.17, 最大误差为0.42, 最小误差接近于0, 误差数据为83.3%, 分布在0—0.20区间内. 2 数据处理的结果 2.1 日性关系和日际形态模型 2008年5月12日14时28分 (地方时) 四川省汶川地区发生MS8.0地震, 震中位于103.4°E、31°N, 震源深度约14km.本文使用COSMIC掩星2级数据, 应用PCA方法及最小二乘拟合法拟合设计区域电离层模型, 尝试探测汶川地震震前电离层扰动. 应用PCA方法选取2008年3月1日—5月14日地方时最大电子密度数据.首先, 网格设计为经度98°—108°E和纬度26°—36°N, 网格大小5°×5°, 分割成4个网格, 即Y1, Y2, Y3和Y4, 如图1所示.然后, 每个网格内3月1日—5月14日数据以15天为间隔, 分割成5个时段, 再以3小时为间隔将每段数据分割成8个时段, 每个时段如果有多个值则取均值, 共计40个点, 每个数据点为3h间隔的时均值, 按时间序列表示为 从而得到矩阵 最后, 计算Y的主成分信号, 显示时间分割点的各主成分能量百分比变化. 模型建立选取2008年4月27日—5月11日12:00—14:00 (地方时) 、200—500km高度电子密度数据, 以50km分层拟合模型计算δ.δ局部显示60°—140°E和10°—60°N, 掩星事件分布如图2所示.同时选取2008年2月12日—5月11日12:00—14:00 (地方时) 、300—350km高度电子密度数据拟合模型计算δ. 通常情况影响电离层的因素包括很多方面:一是太阳活动性, 包括太阳电磁辐射、太阳黑子数、耀斑等;第二是地磁活动水平, 地磁活动引起的电离层扰动一般具有全球性;此外, 一些火山、地震活动也会引起电离层扰动 (Du et al, 2008;刘立波等, 2011;尼鲁帕