阿尔芬波在微电离大气中的二维MHD模拟.pdfVIP

2010 年第二十三届全国空间探测学术交流会论文 太阳风扫过磁层或高能粒子通过磁层运动，在波-粒子的相互作用下，可激发地磁微 脉动。地面上观测的地磁微脉动表现为地磁场的振荡，与磁层、辐射带和电离层状 态及磁活动有关，在部分在外观上为准正弦型，周期为十分之几秒到几分钟，频谱 范围一般为3MHz—5Hz，辐度与周期随观测位置与时间变化，在高纬度较常发生。 放在地面的磁强计能够探测到这些电磁信号，这也证明ULF 波在低空微电离大气中 的传播通道是存在的。在日常地磁观测中附近汽车、地铁的运动就会造成磁场的扰 动。在北京国家地震局白家疃地磁台就因为苹果园地铁站的存在而受到巨大扰动。 这些扰动并不完全是由于汽车、地铁的剩磁而造成的，从扰动的震荡衰减波形判断 是阿尔芬波动。我们将地球低层大气看做是微弱电离的磁化等离子体，其电子密度 和电导率远远小于电离层。可以设想，阿尔芬波的衰减将主要集中在这部分区域。 另外，在阿尔芬波波长范围内地球磁场和空间电离等离子体密度都有非常巨大的变 化，线性的等离子体波激发和传播理论不再适用。为了建立一个比较简单的数学模 型，研究阿尔芬波从地面向电离层区域的传播，我们目前采取层状均匀近似的方法， 即在一定厚度的大气或电离层中采用均匀磁化等离子体的模型，研究阿尔芬波穿透 这些薄层所产生的衰减。 研究地面产生阿尔芬波向近地空间区域的传播，主要是为空间磁场探测中消除 地面人为干扰的因素。随着现代社会工业化的发展，地面电网、大功率发射台也可 能产生大幅度的磁场扰动。这些扰动传播到近地空间而被在轨磁场探测器探测到， 就会干扰电离层动力学探测的研究结果。为了更好地实施近地空间物理探测，有必 要估算地面产生的磁扰动随地球磁力线向近地空间的传播及衰减性质。为今后更好 地分析近地空间磁场探测数据，奠定一些理论基础。 2 近地空间基本介电性质 2.1 大气的微电离 三十公里以下的地球大气基本上是低层大气，电离源主要由 3 部分组成。（1）土壤 中镭、铀和钍等放射性元素辐射的α、β和γ射线；（2 ）由土壤逸入大气的氡和钍等放 射性元素辐射的α、β和γ射线；（3 ）来自宇宙空间的宇宙线。在陆地近地层大气中， 电离源主要来源于土壤中以及由土壤逸入大气的放射性元素，因此大气电离率先随 高度递减，到 2 公里高度附近，电离率达极小值，此后的电离源主要是宇宙线，因 此电离率随高度递增。在海洋近地层大气中，电离源几乎全部来源于宇宙线，因此 电离率随高度单调递增。大约在2～3 公里以上，陆地和海洋上空电离率随高度的变 化趋于一致。到 10～16km 高度附近，电离率达极大值，其值约变动于 2～50cm3s-1 [10] 之间 。这是由于中性大气密度随高度增加呈指数减小，而宇宙线的电离本领则随 着高度递增而增加。然后电离率随高度递减。电子由于质量更小，因此它的迁移率 [11] 比离子大很多，但在中、低层大气中电子的寿命极短，可不考虑 。为方便计算， 本文采用全球平均离子密度。 2.2 大气的分层方法 由于在阿尔芬波波长范围内地球磁场和空间电离等离子体密度都有非常巨大的 变化，因此初步采用局部均匀分层的方法对各层的衰减情况进行计算。大气的分层 2 2010 年第二十三届全国空间探测学术交流会论文 方法如下图所示：将0～50 公里的微电离大气均匀分层，每层取0.1 公里。 对于一定频率的扰动，从参考资料[10] [11]可以得到电子密度和电导率的一些离散 值，然后采用插值的方法求出其它所有点的电导率和电子密度数密度值，将其作为 计算区域的初始值，每层采用中间位置的数值。对于磁场的初值，将地磁场看成一 个偶极子场，磁场大小与距地心距离的三次方成反比：B B r 3 (h +r )3 .其中B 为 0 0 地面的磁场强度，为简单计算取北京地区的平均值 40000nT ；r