AUV地形匹配导航快速收敛滤波

地形匹配定位（terrain aided position,TAP）的似然函数反映了AUV（autonomous underwater vehicle）的位置在空间中的分布概率,由于地形的强非线性、随机性以及测量误差的非高斯分布使得似然函数也表现出非高斯分布的特点。TAP的误差与局部地形特征和地形测量误差密切相关,由于现有的方法未考虑局部地形特征,仅考虑了测量误差的统计置信区间,使得TAP置信区间的估计结果明显偏小。为解决TAP置信区间的估计问题,建立了TAP定位点的跳变模型。设TAP定位点X^p可以向搜索区间内任一点跳变,且向某一点的跳变概率与该点的似然函数值正相关,X^p向某一点跳变的置信度小于α时,认为x^α不会向该点跳变,该点设为置信区间的边界点。另外,设地形匹配定位点的置信区间内匹配残差平方和函数为二次曲面,而X^p视为该曲面的待估计参数,则可以通过曲面参数的置信区间估计方法获得1-α置信度下的置信区间。新方法得到的置信区间范围大于现有的估计方法,试验结果表明,测量波束较少时,置信区间估计会出现异常,增加测量波束可以提高潮差和测量误差的估计精度,从而提高置信区间的估计精度,但测量误差非高斯分布条件下的补偿方法仍然需要进一步研究。     

顾及双差残差反演GPS信号方向的斜路径水汽含量

给出了顾及GPS双差残差反演斜路径水汽SWV的解算流程;然后详细给出了双差残差到非差残差的转化算法,并对算法进行了改进;最后利用并址的GPS和WVR实测数据对反演SWV算法进行了验证,结果证实改进的反演算法能以优于4 mm精度近实时估算SWV值,与目前国际研究精度在同一量级。     

利用残差统计检验评估区域GPS大地测量网解算方案

采用美国喷气推进实验室JPL发展的GIPSY软件解算区域GPS网，比较了固家精密星历和同时固定精密星历及卫星钟参数2种解算方案，残差统计检验表明前者的左中仍保留了部分未解出的信号而后者的残差接近于白噪声高斯分布，解算结果与ITRF96的比较和对重复率的统计表明，残差分布合理的解算较优，简单的线差统计检验的计算可提供重要的解算评估。     

Critical evaluation and geological correlation of teleseismic phase data from Indian seismological stations

We evaluated the quality of seismic phase data from Indian seismological stations through the analysis of teleseismic travel times reported during 1976–83 and infer that only WWSSN stations (NDI, SHL, POO, KOD) apart from GBA and HYB can be rated satisfactory while the majority of stations (more than 40) produce very poor quality data sets. Detailed analysis of teleseismic P-wave travel time residuals shows that while the average structure of the upper mantle beneath India has high velocity (negative residuals) there are marked lateral variations. In particular, three zones of anomalous positive residuals (low velocity) are observed: one beneath the north western part of the Deccan trap, the second covering the southernmost peninsula (granulite terrain) and a third rather localized one, to the north of Delhi coinciding with Delhi-Haridwar ridge. New Delhi exhibits strong negative residuals in the E-SE quadrant along with negative station anomaly, implying that it is underlain by an anomalous high velocity crust/upper mantle. The negative residuals observed over India, continue beneath the Himalaya till the south of Lhasa but change sign further northward, suggesting the northern limit of the Indian upper mantle structure.     

应用地震拓扑预测对华北地区进行地震危险趋势估计

本文应用地震拓扑预测方法,对华北地区进行了地震长期趋势估计。通过实例计算证明,其拓扑预测模型精度较高,可以预测该区未来强震的危险性。     

Leaching potential and redox transformations of arsenic and selenium in sediment microcosms with fly ash

The unintended release of coal ash to the environment is a concern due to the enrichment of contaminants such as arsenic (As) and selenium (Se) in this solid waste material. Current risk assessments of coal ash disposal focus on pH as the primary driver of leaching from coal ash. However, redox speciation of As and Se is a major factor for their mobilization potential and has received much less attention for risk assessments, particularly in disposal scenarios where coal ash will likely be exposed to microbially-driven redox gradients. The aim of this study was to demonstrate the differences of aerobic and anaerobic conditions for the leaching of As and Se from coal ash. Batch sediment-ash slurry microcosms were performed to mimic an ash spill scenario and were monitored for changes in As and Se speciation and mobilization potential. The results showed that the dissolved As concentrations were up to 50 times greater in the anaerobic microcosms relative to the aerobic microcosms during the two week incubation. This trend was consistent with As redox speciation determined by X-ray absorption spectroscopy, which indicated that 55% of the As in the solid phase at the end of the experiment was present as As(III) (a more leachable form of arsenic relative to As(V)). In the aerobic microcosms, only 13% of the As was As(III) and the rest was As(V). More than half of the Se was present as Se(IV) in the original fly ash and in the aerobic microcosms, while in the anaerobic microcosms Se was gradually transformed to less soluble Se(0) species. Likewise, dissolved Se concentrations were up to 25 times greater in the aerobic microcosms relative to anaerobic conditions. While the overall observations of As and Se mobilization potential from coal ash were consistent with expectations for aqueous and solid phase speciation of these elements, the findings directly show the relevance of these processes for coal ash disposal. These results highlight the need to select appropriate environmental parameters to include in risk assessments as well as provide potential geochemical monitoring tools through the use of dissolved Se/As ratios to determine the redox conditions of ash storage and spill sites.     

利用Allan方差分析GPS非差随机模型特性

非差GPS定位中,通常采用以观测值服从高斯白噪声条件的估计准则进行参数求解。研究表明卫星端误差、传播路径误差、测站环境误差等会破坏观测值的白噪声特性,并且未模型化误差同样具有不利影响。这不仅破坏了估计准则的假设条件,而且部分非白噪声有可能被状态参数吸收,影响估计的准确性。本文将观测值白噪声、有色噪声和未模型化误差一同纳入GPS非差随机模型,以验后残差来表征GPS数据的随机特性,进行Allan方差分析,研究噪声成分及其参数。结果表明,GPS非差噪声组合主要为WN+GM,相位白噪声为2.392mm,GM过程噪声为4.450mm/s,相关时间为52.074s,伪距白噪声为0.936m,GM过程噪声为0.833m/s,相关时间为14.737s,相位的GM过程噪声与卫星相关性较大,而其余噪声则与测站相关性较大,大量分析结果表明GPS非差随机模型并不服从高斯白噪声假设,有待精化。     

抗差加权整体最小二乘模型的牛顿-高斯算法

基于加权整体最小二乘的牛顿-高斯迭代算法,提出了一种抗差加权整体最小二乘模型。利用标准化残差构造权因子函数,并采用中位数法获得具有抗差性的单位权中误差估值,能同时实现观测空间和结构空间抗差。为获得标准化残差,利用线性近似的协因数传播律推导了加权整体最小二乘残差协因数阵的表达式,并给出模型的迭代计算方法。试验结果表明:对于加权整体最小二乘的粗差处理问题,本文提出的方法具有良好的抗差性能,参数估值与不含粗差时加权整体最小二乘的结果没有显著的差异,性能优于直接由残差构造的稳健加权整体最小二乘模型。     

Investigation of the lithosphere beneath the Vogelsberg volcanic complex with P-wave travel time residuals

With the aim of investigating the P-wave velocity structure below the Tertiary volcano Vogelsberg, a network of 10 mobile short period seismograph stations was installed in May 1987 for a period of 20 months. P-Wave travel time residuals relative to the station Kleiner Feldberg/Taunus (TNS) were determined for 168 seismic events using the Jeffreys - Bullen travel time tables. At all stations the relative residuals showed a positive sign, indicating a low velocity zone beneath the Vogelsberg. Maxima were found in the northern part of the Vogelsberg (station VAD +0.5 s) and in the region of the Amöneburger Basin (station RAU +0.28 s).The travel time residuals were inverted using the tomographic inversion method of Aki et al. (1977). The slowness perturbations of the single blocks were calculated relative to a crustal and upper mantle model of the Rhenish Massif. The results show an intracrustal low velocity body (about –9%) striking in a Variscan direction and underlying the north-eastern part of the Vogelsberg, and another velocity minimum (about – 6%) in the region of the Am6neburger Basin. In the lower crust and the upper mantle the velocities are reduced by about 4% relative to the starting model.The Variscan alignment of the low velocity zone under the Vogelsberg correlates with results of other geological studies. It can be assumed that during the rifting phase of the Upper Rhinegraben Variscan lineations have been reactivated, favouring uprising of magma along these old structures. The position and extension of the low velocity zone correlate with the assumed sediment distributions in the area of investigation. This may account for about one-half of the observed anomaly. The reason for the velocity reduction of about 4% in the entire underground region of the Vogelsberg down to a depth of about 70 km can be explained by the intensive fracturing of the lithosphere, caused by thermal and pressure gradients during the magma eruption process.