利用高频天波返回散射反演电离层水平不均匀结构

高频天波返回散射探测作为重要的电离层探测手段,能够实现遥远区域电离层空间上的连续监测,探测获取的返回散射扫频电离图显示了探测频率-群路径-回波能量三者之间的关系.由于电离图包含了探测路径上的电离层状态信息,通过对其反演可以实时获取大面积范围的电离层参数.本文提出了一种基于解空间约束的返回散射前沿反演算法,能够重构电离层水平不均匀结构.针对反演非线性问题,采用Newton-Kontorovich方法进行求解,同时又引入了求解不适定问题的Tikhonov正则化方法,有益于解的稳定性和唯一性.利用模拟数据和实测数据分别对本文建立的算法进行了验证,并与Fridman和Fridman于1994年提出的反演方法进行了对比.结果表明,本文算法反演结果稳定,对返回散射前沿判读误差不敏感,与Fridman和Fridman 1994年方法相比,本文方法对电离层局部精细结构反演更加准确,具有较高的反演精度.本文提出的算法不但能够反演白天和夜间这种电离层较平稳时期的电离层状态,而且对于日出/日落时段等电子浓度分布变化较快情形下的电离层,也有很好的反演效果,表明了该算法在处理复杂多变的实际探测的返回散射电离图中的应用价值.     

1987年9月23日日环食的电离层波动现象

本文给出了利用高频天波返回散射雷达对1987年9月23日日环食电离层效应的观测结果.分析表明,日食期间在日环食带电离层存在明显的不均匀结构和小尺度的波动,而日食过后在见食电离层区域出现大尺度的波动现象.     

1987年9月23日日环食的电离层波动现象

本文给出了利用高频天波返回散射雷达对1987年9月23日日环食电离层效应的观测结果.分析表明,日食期间在日环食带电离层存在明显的不均匀结构和小尺度的波动,而日食过后在见食电离层区域出现大尺度的波动现象.     

1987年9月23日日环食的电离层波动现象

本文给出了利用高频天波返回散射雷达对1987年9月23日日环食电离层效应的观测结果.分析表明,日食期间在日环食带电离层存在明显的不均匀结构和小尺度的波动,而日食过后在见食电离层区域出现大尺度的波动现象.     

电离层赤道异常对返回散射电离图的影响

本文给出了太阳黑子高年冬季跨电离层赤道异常区传播试验的高频返回散射电离图.经过分类研究,指出了有一类双回波区组合型P ′（f）返回散射电离图是一种新的形态.通过射线追踪合成技术分析,表明了这类型实验返回散射电离图是沿探测路径存在双峰电子浓度分布的反映,时序电离图的变化特征显示了赤道异常区的电子浓度变化规律,因而证明了高频返回散射技术在研究赤道异常区电离层峰运动的可能性及其潜力.     

电离层赤道异常对返回散射电离图的影响

本文给出了太阳黑子高年冬季跨电离层赤道异常区传播试验的高频返回散射电离图.经过分类研究,指出了有一类双回波区组合型P ′（f）返回散射电离图是一种新的形态.通过射线追踪合成技术分析,表明了这类型实验返回散射电离图是沿探测路径存在双峰电子浓度分布的反映,时序电离图的变化特征显示了赤道异常区的电子浓度变化规律,因而证明了高频返回散射技术在研究赤道异常区电离层峰运动的可能性及其潜力.     

利用高频返回散射技术探测远地核爆炸电离层效应

本文叙述利用高频返回散射技术探测我国1976年的一次低空大当量核爆炸电离层效应的结果。返回散射探测设备设于离爆炸点2205km的新乡。本文得到了核爆炸产生的附加电离区和电离层扰动的参数。这些结果与爆炸区附近其他手段观测的结果相符。     

中国地区电离层TEC暴扰动研究

电离层总电子含量(TEC)是空间天气研究和监测预报的重要参量.本文引入了电离层TEC扰动指数DI, 对青岛等6个台站的DI数据进行分析,选取DI>0.35(DI≤-0.30)作为正(负)相电离层TEC扰动的强度标准,并以连续6 h及以上的DI满足该值来判定电离层TEC暴扰动事件.对电离层TEC暴扰动事件的统计分析表明,在地方时日落后至子夜前为发生高峰时段,正(负)相暴扰动事件平均持续时间约为10.9 h(10.5 h),正相暴发生率以冬季为多,夏季为少,而负相暴则以夏季略高.发现位于赤道异常驼峰区的广州站和位于高中纬度的海拉尔站比典型中纬地区的北京站电离层TEC暴扰动更易发生,且低纬地区以正相暴扰动为主.分析表明,约有70%的电离层TEC暴扰动伴随着有地磁扰动,但是电离层TEC暴扰动并不完全由地磁扰动所引起,强烈气象活动等局地环境因素也可能对电离层TEC暴扰动有着重要影响.     

中国地区电离层TEC暴扰动研究

电离层总电子含量(TEC)是空间天气研究和监测预报的重要参量.本文引入了电离层TEC扰动指数DI,对青岛等6个台站的DI数据进行分析,选取DI＞0.35(DI≤-0.30)作为正(负)相电离层TEC扰动的强度标准,并以连续6h及以上的DI满足该值来判定电离层TEC暴扰动事件.对电离层TEC暴扰动事件的统计分析表明,在地方时日落后至子夜前为发生高峰时段,正(负)相暴扰动事件平均持续时间约为10.9 h(10.5 h),正相暴发生率以冬季为多,夏季为少,而负相暴则以夏季略高.发现位于赤道异常驼峰区的广州站和位于高中纬度的海拉尔站比典型中纬地区的北京站电离层TEC暴扰动更易发生,且低纬地区以正相暴扰动为主.分析表明,约有70％的电离层TEC暴扰动伴随着有地磁扰动,但是电离层TEC暴扰动并不完全由地磁扰动所引起,强烈气象活动等局地环境因素也可能对电离层TEC暴扰动有着重要影响.     

武汉地区中尺度电离层声重波扰动的变化特性

利用武汉电离层观象台高频多普勒台阵的覆盖太阳活动高、低年份,长达５年的连续观测数据,采用小波分析等方法估算电离层声重波扰动（ＴＩＤ）的传播参量,通过这些参量对武汉地区电离层扰动形态和变化规律进行了系统分析研究．结果表明,观测到的中尺度电离层声重波扰动（ＭＳＴＩＤ）存在二个显著季变化,在传播速度和周期上有明显差异的优势传播方向：一个指向东北方,传播的方位角主要分布在３０°─７０°之间（０°为正北,以顺时针方向表示传播方位角）,它在夏季出现率最大,冬季基本消失;另一个优势方向指向正南,方位角主要分布在１５０°─２２０°范围,主要出现在冬季．文中还给出了ＭＳＴＩＤ的年、日变化,并进一步探讨了其变化特性的可能形成机理．     

一种改进的线性走时插值射线追踪算法

线性走时插值法（LTI）在走时的计算中,由于射线方向考虑不全,计算得到的节点走时不一定最小,导致追踪的射线路径无法满足最小走时.针对这一问题,本文提出了一种改进的射线追踪算法,通过采用多方向的循环计算,得到所有计算节点的最小走时,使追踪到的射线路径能真正满足最小走时,以确保射线追踪的精度.模拟实验结果表明,在介质速度变化剧烈的结构中,该算法与传统的LTI算法相比,有效地提高了射线追踪的精度.     

电离层人工调制激发的VLF波在磁层的传播特性及应用研究

电离层人工调制可以激发甚低频（VLF）波,其中向上传播进入磁层的VLF波,不但能够用来研究磁层中的各种物理现象,且具有人工沉降高能粒子,消除辐射带等实际用途.本文使用射线追踪方法,模拟电离层调制激发的VLF波在磁层的传播路径,分析激发纬度和调制频率对传播路径和传播特性的影响;并基于低频波的色散方程和波粒共振条件,分析VLF波传播路径上与磁层高能粒子的最低共振能及其分布.研究表明,VLF波通过在磁层来回反射向更高的L-shell传播,最终稳定在某一L-shell附近.以较低的调制频率或者从较高的纬度激发的VLF波能够传播到更高的L-shell,但是,当激发纬度过高时,低频波也可能不发生磁层反射而直接进入电离层和大气层.低频波在磁层的传播过程中,在较高的纬度或者较低的L-shell能够与较高能量的电子发生共振相互作用,在较高的L-shell并且低纬地区,能够与较低能量的电子发生共振相互作用.共振谐数越高,能发生波粒共振的电子能量越高.     

The synthesis of travelling ionospheric disturbance (TID) signatures in HF radar observations using ray tracing

Characteristic signatures are often observed in HF radar range-time-intensity plots when travelling ionospheric disturbances (TIDs) are present. These signatures, in particular the variation of the F-region skip distance, have been synthesised using a ray tracing model. The magnitude of the skip variation is found to be a function of the peak electron density perturbation associated with the TID and radar frequency. Examination of experimental observations leads to an estimate of the peak electron density perturbation amplitude of around 25% for those TIDs observed by the CUTLASS radar system. The advantage of using the skip variation over the radar return amplitude as an indicator of density perturbation is also discussed. An example of a dual radar frequency experiment has been given. The investigation of the effect of radar frequency on the observations will aid the optimisation of future experiments.     

Nonlinear ray perturbation theory with its applications to ray tracing and inversion in anisotropic media

A comprehensive approach, based on the general nonlinear ray perturbation theory (Druzhinin, 1991), is proposed for both a fast and accurate uniform asymptotic solution of forward and inverse kinematic problems in anisotropic media. It has been developed to modify the standard ray linearization procedures when they become inconsistent, by providing a predictable truncation error of ray perturbation series. The theoretical background consists in a set of recurrent expressions for the perturbations of all orders for calculating approximately the body wave phase and group velocities, polarization, travel times, ray trajectories, paraxial rays and also the slowness vectors or reflected/transmitted waves in terms of elastic tensor perturbations. We assume that any elastic medium can be used as an unperturbed medium. A total 2-D numerical testing of these expressions has been established within the transverse isotropy to verify the accuracy and convergence of perturbation series when the elastic constants are perturbed. Seismological applications to determine crack-induced anisotropy parameters on VSP travel times for the different wave types in homogeneous and horizontally layered, transversally isotropic and orthorhombic structures are also presented. A number of numerical tests shows that this method is in general stable with respect to the choice of the reference model and the errors in the input data. A proof of uniqueness is provided by an interactive analysis of the sensitivity functions, which are also used for choosing optimum source/receiver locations. Finally, software has been developed for a desktop computer and applied to interpreting specific real VSP observations as well as explaining the results of physical modelling for a 3-D crack model with the estimation of crack parameters.     

Two-point ray tracing in 3-D

Algorithm for determination of all two-point rays of a given elementary wave by means of the shooting method is presented. The algorithm is designed for general 3-D models composed of inhomogeneous geological blocks separated by curved interfaces. It is independent of the initial conditions for rays and of the initial-value ray tracer. The algorithm described has been coded in Fortran 77, using subroutine packages MODEL and CRT for model specification and for initial-value ray tracing.     

Theory of anisotropic dynamic ray tracing in ray-centred coordinates

A 4×4-propagator matrix formalism is presented for anisotropic dynamic ray tracing, including the propagation across curved interfaces. The computations are organised in the same way as in ervený's well-known isotropic propagator matrix formalism. Attention is paid to cases where double eigenvalues of the Christoffel matrix result in unstable expressions in the dynamic ray tracing system, but where geometrical spreading is well-behaved.     

弯曲射线追踪中Dijkstra算法的改进与实现

文章针对图论中寻找最短路径的Dijkstra算法内存占用量大,效率低的缺点,对该算法进行了改进,修改后的算法计算效率是原来的四倍,内存使用量和图中节点数呈线性关系.在此基础上,用新算法求出了激发点和接收点的最短走时路径,并由激发接收点的旅行时结合联合迭代法对理论模型和实际场地进行了反演.结果表明：和直射线追踪相比,弯曲射线路径能更好地反演出地质体内部的速度场分布.理论模型和实际探测结果证实改进后的算法是有效的.     

Common-ray tracing and dynamic ray tracing for S waves in a smooth elastic anisotropic medium

Anisotropic common S-wave rays are traced using the averaged Hamiltonian of both S-wave polarizations. They represent very practical reference rays for calculating S waves by means of the coupling ray theory. They eliminate problems with anisotropic-ray-theory ray tracing through some S-wave slowness-surface singularities and also considerably simplify the numerical algorithm of the coupling ray theory for S waves. The equations required for anisotropic-common-ray tracing for S waves in a smooth elastic anisotropic medium, and for corresponding dynamic ray tracing in Cartesian or ray-centred coordinates, are presented. The equations, for the most part generally known, are summarized in a form which represents a complete algorithm suitable for coding and numerical applications.     

A note on dynamic ray tracing in ray-centered coordinates in anisotropic inhomogeneous media

Dynamic ray tracing plays an important role in paraxial ray methods. In this paper, dynamic ray tracing systems for inhomogeneous anisotropic media, consisting of four linear ordinary differential equations of the first order along the reference ray, are studied. The main attention is devoted to systems expressed in a particularly simple choice of ray-centered coordinates, here referred to as the standard ray-centered coordinates, and in wavefront orthonormal coordinates. These two systems, known from the literature, were derived independently and were given in different forms. In this paper it is proved that both systems are fully equivalent. Consequently, the dynamic ray tracing system, consisting of four equations in wavefront orthonormal coordinates, can also be used if we work in ray-centered coordinates, and vice versa. vcerveny@seis.karlov.mff.cuni.cz     

A fast and global two point low storage optimization technique for tracing rays in 2D and 3D isotropic media

We present the problem of tracing rays in 2D and 3D heterogeneous isotropic media as a set of optimization problems. Each optimization problem is obtained by applying Fermat's principle to an approximation of the travel time equation from a fixed source to a fixed receiver. We assume a piecewise linear ray path that simplifies the computations of the problem, in the same way Mao and Stuart suggested in a very recent paper. Here, instead, the reflector geometry and the velocity function are computed by using nonuniformly biharmonic splines. On the other hand, to solve the optimization problem we use the Global Spectral Gradient method. This recent developed optimization scheme is a low storage optimization technique that requires very few floating point operations. It only requires the gradient of the travel time function, and it is global because it converges independently of the initial guess, that is, it does not require a close initial ray path. These three properties of the optimization method and the assumption of piecewise linear rays make this ray tracing scheme a very fast, global and effective method when estimating velocities via tomography. Moreover, in a homogeneous stratified or dipped media, any solution of the optimization problem is the best solution, i.e., it is the global minimum, no matter what numerical approach is used. We present some numerical results that show the computational advantages and the performance of this ray tracing in homogeneous and heterogeneous media.