Efficient reflection waveform inversion using a locally normalized zero-lag correlative objective function

Full-waveform inversion is characterized by cycle-skipping when the starting background model differs significantly from the true model and low-frequency data are unavailable. To mitigate this problem, reflection waveform inversion is applied to provide a background velocity model for full-waveform inversion. This technique attempts to extract background velocity updates along the reflection wavepath by matching the reflection waveforms. However, two issues arise during the implementation of reflection waveform inversion: amplitude and efficiency. The amplitude is always underestimated due to the complex subsurface parameter (i.e. the source signature, density, attenuation etc.). This makes it unreasonable to match the reflection amplitude involved in waveforms, especially in the filed data cases. In addition, generating the background velocity gradient requires the simulation of the reflection wavefield. However, simulating the reflection wavefield is time-consuming. To address the former, we introduced a locally normalized objective function, while for the latter, we used an efficient strategy by avoiding the explicit generation of the reflection wavefield. Results show that applying the proposed method to both synthetic and field data can provide a good background velocity model for full-waveform inversion with high efficiency.     

Theory of the electromagnetic transient response of the ionosphere for plane wave excitation

Summary The reflection and transmission of a transient plane wave by an ionosphere with a linear electron density profile and constant collision frequency are calculated. The exact frequency domain solution is inverted numerically, and the resultant reflected and transmitted waveforms are found to satisfy causality. A linearly varying, nondispersive dielectric half-space is also considered, but the nonphysical nature of the model causes the reflected waveform to be noncausal. Using convolution the waveforms for a step modulated sinusoid excitation are also obtained from the impulse response of the ionosphere. The results provide insight into the physical mechanism involved in the ionospheric propagation of HF electromagnetic pulsed signals.     

基于闪电双极性窄脉冲事件波形观测日出日落电离层D层变化

利用闪电作为辐射源来探测电离层D层是近年来国外学者研究的热点.本文基于江淮流域六站同步闪电观测站网,实现了一种利用闪电双极性窄脉冲事件（Narrow Bipolar Events,NBE）来探测电离层D层等效高度的方法.基于此方法,对两次分别发生在日出和日落阶段的雷暴分析结果显示,上述两个阶段D层反射高度变化特征呈现显著的不对称性：日出期间D层反射高度随时间线性降低速率为5.9 km/h;而日落期间D层反射高度随时间线性递增速率为8.6 km/h.综合日间太阳耀斑期间D层反射高度剧烈波动的观测事实,与日出、日落期间D层特征变化,结果表明太阳辐射电离中性大气分子的电子生成机制在日间D层电子密度变化中占主导地位.本文结果展现了利用NBE事件监测电离层D层变化特征的可行性,这一方法与基于地闪回击波形的D层探测方法结合在一起,有望把现有具有闪电事件定位和电磁波波形记录能力的闪电观测站网扩展为实时对电离层D层时空变化监测的网络.     

利用半参数核估计法预报全球电离层总电子含量

本文将半参数平差模型引入电离层球谐函数系数的预报中,建立了半参数球谐函数模型(Semiparametric-Spherical Harmonic,Semi-SH)来预测全球电离层总电子含量.首先,通过快速傅里叶变换获得球谐函数系数的周期和振幅,将振幅高的主周期归入趋势函数,振幅低的剩余周期归入随机信号,建立了半参数模型,同时利用核估计方法拟合趋势函数,解算随机信号,并在时间域上进行外推,得到了预报时间的球谐函数系数,代入15阶电离层球谐函数模型,最后得出电离层总电子含量(Total Electron Content,TEC)的预报值.本文基于欧洲定轨中心(CODE)发布的球谐函数系数进行电离层TEC长期预报和短期预报分析,其中长期预报采用四年预报两年的模式对球谐函数系数进行预报,短期预报设计了三个算例,采用前30天预报后一天的模式,分别预报1天、滑动预报7天和滑动预报30天.实验结果表明:长期预报能够较好地反映全球电离层TEC的变化趋势和波动情况,Semi-SH模型对全球电离层TEC平均值(Mean TEC global,MTECglobal)的拟合值和预报值与MTECglobal实际值的相关系数分别为0.8743和0.8010,呈现出高度相关性.短期预报中,在太阳活动高年和太阳活动低年,Semi-SH模型在中纬度地区预报精度较CODE发布的电离层TEC 1天预报产品(CODE′S 1-Day Predicted GIM,C1PG)有较大提升,在高纬度与低纬度地区两种模型预报精度相当;Semi-SH模型在太阳活动高年和太阳活动低年30天滑动预报精度的均值均高于C1PG模型.实验结果说明了Semi-SH模型预报电离层TEC值的有效性.     

Ionospheric response to the partial solar eclipse of March 29, 2006, according to the observations at Nizhni Novgorod and Murmansk

The results of observations of the solar eclipse ionospheric effects on March 29, 2006, are presented. The observations were conducted using the partial reflection method near Nizhni Novgorod and the vertical sounding method at the automatic ionospheric station near Murmansk. It has been obtained that the electron density at altitudes of 77 and 91 km decreases by a factor of more than 4; in this case the response of the ionosphere at an altitude of 91 km lags behind the eclipse maximum phase on the Earth by approximately 20 min. It has been established that the eclipse in the E and F1 regions of the polar ionosphere causes a change in the electron density by 15–20%. The delay time of this effect varies from 12 to 24 min depending on the altitude. It has been registered that the reflection virtual altitude at altitudes of the ionospheric F region increases in Murmansk and Nizhni Novgorod.     

Geomagnetic field effect on the ionospheric contribution to the error of navigation satellite systems

When using global navigation satellite systems (GNSSs) for high-precision measurements, one should consider high-order errors. The ionospheric second-order error caused by the geomagnetic field is approximately proportional to the total electron content. Therefore, this error can be taken into account by modifying the coefficients in an “ionosphere-free” combination of GNSS measurements at two frequencies. This study checks the approximations underlying this modification. We reveal that these approximations are valid and the results depend weakly on the accuracy of ionospheric parameters used a priori for calculating the coefficients of the modified two-frequency formula. In addition, we investigate how the choice of a model of the Earth’s magnetic field affects the second-order ionospheric error.     

基于电离层反射的袖珍云闪(CID)三维定位研究

袖珍云闪是一类区别于常规闪电放电过程的特殊放电现象,能够同时产生极强的高频和低频辐射信号,其低频辐射信号在电离层与地面之间反射后能够在电场变化波形上形成电离层反射脉冲对.电离层反射信号与原信号的时间差包含着放电三维位置和电离层高度的信息,而借助于多站闪电探测网络的同步观测就能够反演这些信息.基于这一规律,本文发展了一种对袖珍云闪实现三维定位的新方法.这种方法不仅能够对大范围内的袖珍云闪实现准确的三维定位,同时还能够反演电离层的高度,是一种潜在的研究电离层相关性质的有效手段.通过将定位结果与雷达回波比较,证明这种方法具有较高的精度.利用这种方法,计算了5489例正极性袖珍云闪和1400例负极性袖珍云闪的放电高度,发现正极性袖珍云闪主要集中在7~14 km,而负极性袖珍云闪达到了15~18 km.负极性袖珍云闪的放电高度总体上与对流层顶高度相当,其数量相比于正极性袖珍云闪明显偏少,因此很可能产生于较为罕见的极旺盛的雷暴过程中.     

Shear waves at a corrugated interface between anisotropic elastic and visco-elastic solid half-spaces

Rayleigh’s method of approximation is employed to find out the reflection and transmission coefficients due to an incident plane SH wave at a corrugated interface between a laterally and vertically inhomogeneous anisotropic elastic solid half-space and a laterally and vertically inhomogeneous isotropic visco-elastic solid half-space. The lateral and vertical inhomogeneities are described by the exponential variations of elastic parameters. The formulae of reflection and transmission coefficients are derived in closed form for the first-order approximation of the corrugation. The effects of the corrugation of the interface, the inhomogeneity, the anisotropy, the visco-elasticity and the frequency of the incident wave on these coefficients are studied analytically and numerically for a specific model containing a periodic interface. The results of earlier workers have been reduced as particular cases from the present formulation.     

Effects of anelasticity on reflection and transmission coefficients

It is known that the reflection and transmission coefficients used in the zeroth order approximation of asymptotic ray theory (ART) are identical to those obtained for the plane wave impinging on a plane interface separating two perfectly elastic half-spaces. We have used ART to compute reflection and transmission coefficients for two viscoelastic media separated by a plane interface. Our method is different from the plane-wave approach because the ART approach requires only a local application of the boundary conditions both for the eikonal and the ray amplitudes. Several types of viscoelastic media were studied. For a given model, the elastic case was emulated by setting all the quality factors Q equal to each other. Several anelastic cases were computed by keeping the same velocities and densities while changing the Qs. The quality factor is a relatively difficult parameter to measure exactly. Hence elastic coefficients are used in most synthetic seismogram computations, and the quality factors are chosen from experimental measurements or simply estimated. From these computations, amplitude and phase differences between elastic coefficients and coefficients for dissipative media are observed in some cases. These differences show the importance of knowing the exact values of Q. Incorrect Q values can lead to unrealistic moduli and to noticeable phase differences of these viscoelastic coefficients.     

接收函数确定Moho面速度和密度跃变的方法研究及应用

通过对单层模型反射和透射系数的推导,提出了利用接收函数一次转换波和多次波确定Moho面速度和密度跃变的速度-密度跃变(δβ-δρ)扫描叠加方法.利用反射率法计算了不同模型的远震理论地震图,按照与处理实际观测波形一致的方法和流程计算了理论接收函数;根据不同模型数值试验结果,深入分析了界面速度和密度跃变对接收函数震相幅度的影响.利用(δβ-δρ)扫描叠加方法,对理论接收函数进行了数值试验,结果证明了该方法的可行性.最后将该方法应用于位于青藏高原东北缘的高台(GTA)台和兰州(LZH)台,确定了两个台站下方Moho面的速度跃变分别约为(19±1)%和(20±1)%,密度跃变最小值为(4±2)%和(6±2)%.     

Removing multiple reflections from the F2 layer to improve Autoscala performance

This work describes a simple filter able to eliminate the traces caused by second-order reflection from ionograms. This filter is applied to Autoscala in order to smooth out cases in which the autoscaling of the ionogram is mistaken because the second-order F2 layer reflection is identified as first order. A dataset of 32,626 ionograms recorded at the ionospheric observatory of Rome was used to test the filter.     

Numerical simulation of the structure of the high-latitude ionospheric <Emphasis Type="Italic">F</Emphasis> region during meridional HF propagation

A previously developed model of the high-latitude ionosphere is used to calculate the distribution of the ionospheric parameters in the polar region. A specific method for specifying input parameters of the mathematical model, using the experimental data obtained by the method of satellite radio tomography, is used in this case. The spatial distributions of the ionospheric parameters characterized by a complex inhomogeneous structure in the high-latitude region, calculated with the help of the mathematical model, are used to simulate the HF propagation along the meridionally oriented radio paths extending from middle to high latitudes. The method for improving the HF communication between a midlatitude transmitter and a polar-cap receiver is proposed.     

Modeling of vertical sounding ionograms for nonmonotonous ionospheric layers

The problem of radiosignal reflection from a plane layered isotropic ionosphere has been considered. The vertical sounding ionograms for ionospheric layers with a complicated structure have been modeled. It has been demonstrated that the structure of reflected signals depends on the degree of ionospheric irregularities. Diffuse reflection of radiosignals, caused by the vertical irregularity of the ionosphere, has been detected.     

利用GPS计算TEC的方法及其对电离层扰动的观测

在总结用ＧＰＳ研究电离层电子总量ＴＥＣ的数据处理方法基础上,分析了利用伪距观测量和载波相位观测量计算电离层ＴＥＣ的特点及误差来源．在处理过程中考虑了卫星的硬件延迟偏差,分析了应用ＩＲＩ模型进行接收机硬件延迟偏差修正的可能性,发现利用少量ＧＰＳ数据和ＩＲＩ模型修正接收机硬件延迟偏差有一定的困难．最后,利用一些ＧＰＳ观测数据有针对性地研究了电离层对若干次扰动事件的响应．包括一次大的太阳耀斑期间的电离层ＴＥＣ变化、一次较典型的电离层行扰以及日食期间的电离层ＴＥＣ的相对变化等电离层物理问题．结果表明,利用该方法计算ＴＥＣ的精度可满足电离层扰动现象的研究．     

Inversion of ocean-bottom seismometer (OBS) waveforms for oceanic crust structure: a synthetic study

The waveform inversion method is applied—using synthetic ocean-bottom seismometer(OBS) data—to study oceanic crust structure. A niching genetic algorithm(NGA) is used to implement the inversion for the thickness and P-wave velocity of each layer, and to update the model by minimizing the objective function, which consists of the misfit and cross-correlation of observed and synthetic waveforms. The influence of specific NGA method parameters is discussed, and suitable values are presented.The NGA method works well for various observation systems, such as those with irregular and sparse distribution of receivers as well as single receiver systems. A strategy is proposed to accelerate the convergence rate by a factor of five with no increase in computational complexity; this is achieved using a first inversion with several generations to impose a restriction on the preset range of each parameter and then conducting a second inversion with the new range. Despite the successes of this method,its usage is limited. A shallow water layer is not favored because the direct wave in water will suppress the useful reflection signals from the crust. A more precise calculation of the air-gun source signal should be considered in order to better simulate waveforms generated in realistic situations; further studies are required to investigate this issue.     

Ionospheric support of HF radiocommunication using maps of total electron content

The possibility of determining the field of critical frequencies of the ionospheric F2 layer (foF2) using the maps of the total electron content, constructed based on the registration of signals from satellite radio navigation system of the GPS and GLONASS types, is considered. The calculation of foF2 is based on the SPIM (Standard Plasmasphere-Ionosphere Model) model specifying the ionospheric index of solar activity, which is determined at grid points of the map of the total electron content. The proposed method has been verified using the data of the hourly maps of the total electron content in the North American region during September 1–7, 2005. The variations in the critical frequencies for Boulder and Dyess sites, selected from the reconstructed foF2 maps, were compared with the data of the vertical sounding. The average error is ～10% during the entire period of measurements. The conclusion has been drawn that the proposed method can be used as an ionospheric support of HF radiocommunication in the cases when errors of tenths of MHz in foF2 values are permissible.     

The effect of the electron sound speed on wave propagation in the ionospheric plasma

In this study, the effect of the electron sound speed on the extraordinary wave propagation is calculated without an approximation for either collisional or collisionless cases in the ionospheric plasma by using the real geometry of the Earth’s magnetic field for the Northern Hemisphere. It is observed that there is no remarkable effect on the propagation of the extraordinary wave, especially at reflection altitudes. But it is also observed that the magnitudes of k ² (the square of the wave number) have changed every season, and the phase velocity of wave in warm ionospheric plasma has increased.     

电离层GPS掩星观测改正TEC反演方法

电离层掩星观测中,当低轨卫星（LEO）轨道高度较低时,轨道以上的电离层电子总含量（TEC）对掩星反演的影响不能忽略.目前,一般采用指数函数等外推方法来处理该问题,对反演结果可能引起较大误差.为提高电离层掩星反演精度,本文研究利用LEO处于非掩星一侧GPS观测数据的改正TEC新反演方法.用三维射线追踪程序计算出电离层掩星观测模拟数据,分别应用改正TEC方法和外推方法进行反演,将反演结果与所用模式值进行比较.结果表明：对于轨道高度约800km的GPS/MET掩星模拟数据,外推方法和改正TEC方法反演结果都与模式值基本一致;对于轨道高度约400km的CHAMP掩星模拟数据,外推方法误差较大,改正TEC方法反演结果与模式值相符得较好.将改正TEC方法应用于GPS/MET实测数据的反演,取得了合理的结果.这些说明,改正TEC算法是一种有效的电离层掩星反演方法,尤其是对于轨道较低的LEO的电离层掩星观测反演特别有用.     

STABILIZATION OF NORMAL-INCIDENCE SEISMOGRAM INVERSION REMOVING THE NOISE-INDUCED BIAS*

A main problem in computing reflection coefficients from seismograms is the instability of the inversion procedure due to noise. This problem is attacked for two well-known inversion schemes for normal-incidence reflection seismograms. The crustal model consists of a stack of elastic, laterally homogeneous layers between two elastic half-spaces. The first method, which directly computes the reflection coefficients from the seismogram is called “Dynamic Deconvolution”. The second method, here called “Inversion Filtering”, is a two-stage procedure. The first stage is the construction of a causal filter by factorization of the spectral function via Levinson-recursion. Filtering the seismogram is the second stage. The filtered seismogram is a good approximation for the reflection coefficients sequence (unless the coefficients are too large). In the non-linear terms of dynamic deconvolution and Levinson-recursion the noise could play havoc with the computation. In order to stabilize the algorithms, the bias of these terms is estimated and removed. Additionally incorporated is a statistical test for the reflection coefficients in dynamic deconvolution and the partial correlation coefficients in Levinson-recursion, which are set to zero if they are not significantly different from noise. The result of stabilization is demonstrated on synthetic seismograms. For unit spike source pulse and white noise, dynamic deconvolution outperforms inversion filtering due to its exact nature and lesser computational burden. On the other hand, especially in the more realistic bandlimited case, inversion filtering has the great advantage that the second stage acts linearly on the seismogram, which allows the calculation of the effect of the inversion procedure on the wavelet shape and the noise spectrum.