捷联式航空重力扰动矢量解算及误差补偿

介绍捷联式航空重力矢量测量的基本原理,利用波数相关滤波（WCF）、比力线性校正的方法对重力扰动的水平和垂向分量进行误差补偿,评定内符合精度。对国产捷联式航空重力矢量仪的某次试验数据进行处理,以验证和评估航空重力矢量仪的性能。结果表明,在半波长分辨率为7.5 km时,6条重复测线的重力扰动水平分量经波数相关滤波处理后,东、北向分量的平均内符合精度分别从9.77 mGal、9.18 mGal提高到5.95 mGal、3.83 mGal;对比力的垂向分量线性校正后,将解算的重力扰动垂向分量再用WCF方法处理,其平均内符合精度从1.27 mGal提高到0.59 mGal。     

速度传感反馈地震计的动态范围上限

在开环参数一定的条件下,速度传感反馈地震计中的环路滤波是影响地震计主要技术指标的关键环节之一.本文探讨了一阶极点高通、一阶极点低通、一阶零点高通、双一阶和零阶环路滤波在扩展地震计动态范围上限方面的优缺点.研究指出,对于用闭环反馈直接生成主导二阶极点的速度传感反馈地震计,一阶零点高通环路滤波的优点较多;若允许在环外生成主导二阶极点,则零阶环路滤波也不失为一个较好的方案.     

快速序贯算法

本文讨论了序贯算法的基本原理，与ｋａｌｍａｎ滤波的关系，和序贯算法的渐近性态，得到无先验信息的序贯算法不受初值的影响。实例分析表明，此序贯算法可适于方程系数阵的求逆、求解和病态方程的快速解算。     

Exploratory spectral analysis of hydrological times series

Current methods of estimation of the univariate spectral density are reviewed and some improvements are made. It is suggested that spectral analysis may perhaps be best thought of as another exploratory data analysis (EDA) tool which complements, rather than competes with, the popular ARMA model building approach. A new diagnostic check for ARMA model adequacy based on the nonparametric spectral density is introduced. Additionally, two new algorithms for fast computation of the autoregressive spectral density function are presented. For improving interpretation of results, a new style of plotting the spectral density function is suggested. Exploratory spectral analyses of a number of hydrological time series are performed and some interesting periodicities are suggested for further investigation. The application of spectral analysis to determine the possible existence of long memory in natural time series is discussed with respect to long riverflow, treering and mud varve series. Moreover, a comparison of the estimated spectral densities suggests the ARMA models fitted previously to these datasets adequately describe the low frequency component. Finally, the software and data used in this paper are available by anonymous ftp from fisher.stats.uwo.ca.     

Transient deformation following the 30 January 1997 dike intrusion at Kīlauea volcano,Hawai’i

On 30 January 1997 an intrusion on Kīlauea volcano opened a new fissure within the East Rift Zone (ERZ) at Nāpau Crater, 3 km uprift from the ongoing eruptions at Pu’u ’Ō’ō. The fissure eruption lasted 22 h and opened a 5.1 km long, nearly vertical dike 1.9 m, extending from the surface to a depth of 2.4 km (Owen et al. 2000b). During the eruption, the lava pond at Pu’u ’Ō’ō drained, and eruptions ceased there. Pu’u ’Ō’ō eventually refilled in late February and eruptions resumed there on 28 March 1997. Continuous GPS data show a large transient following the 30 January 1997 dike intrusion. After lengthening 40 cm during the initial eruption, the baseline between two stations spanning the ERZ lengthened an additional 10 cm over the following 6 months. A coastal station KAEP also exhibited transient deformation, as it continued to move southward (5 cm) over the same 6-month period. The baseline between two stations spanning Kīlauea’s summit caldera contracted sharply during the eruption, but gradually recovered to slightly longer than its previous length 2 months after the intrusion. We use the extended network inversion filter (McGuire and Segall 2003) to invert continuous GPS data for volume change of a spherical pressure source under Kīlauea’s summit, opening distribution on a nearly vertical dike in the ERZ and potential slip on a decollement 9 km beneath the south flank. Following the 30 January intrusion, rift extension continued below the initial dike intrusion for the duration of the transient. Decollement slip, regardless of its assumed depth, is not required to fit the data. The modeled transient summit reinflation and rift opening patterns under Nāpau crater coincide with changes in observed behavior of Pu’u ’Ō’ō’s lava pond. Rift opening accelerated while Pu’u ’Ō’ō eruptions paused and began to decelerate after the lava pond reappeared nearly a month after the Nāpau eruption. The transient deformation is interpreted as resulting from shallow accommodation of the new dike volume.     

对重力势二阶垂直偏导数的一种新的取样方法

从目前的研究来看，二阶垂直导数（X和Y两个方向）的求解方法较简单，它采用的五个因子的数字滤波的最小二乘法来优化数据。但是在较小区域面积时（对工程和环境应用中的重力数据），重力异常很容易放大和被区域重力场掩盖，我们认为产生的缘由是表示重力场的x和y方向上的二次多项式引起的。因此，为了逼近重力水平方向的偏导可以设计一个好的滤波器来拟合在xy平面的重力场值，再得出想要的重力垂直方向偏导，文中最后对提出的方法在前人经典模型上和野外实践进行了验证。     

基于启发式图搜索的遥感影像道路半自动提取

启发式图搜索法用于线状目标识别的原理是：用图结构表示边缘点和边缘段，根据启发函数计算顶点权值，在图的路径上建立相应的代价函数，通过在图中搜索对应的最小代价的通道以找到最优路径。图搜索法是一种全局最优方法，它在受噪声影响较大时效果仍然较好。文中使用了启发式图搜索法（A*算法）实现了道路的半自动跟踪。它的基本思路是：首先利用自适应平滑滤波算子进行道路信息增强，然后对传统的道路数学模型进行了进一步的扩展，突出了对道路几何特性和辐射特性的描述，并依此构建图搜索的代价函数，实现了基于启发式图搜索法A*算法的道路半自动跟踪。经实验证明，该方法进行遥感影像的道路半自动提取效果较好。     

华北平原深层淡水在开采条件下接受上覆咸水越流补给——以天津平原为例

分析物探测井曲线和地下水位动态发现，天津平原深层孔隙地下淡水在开采条件下接受上覆咸水的越流补给，而且咸水底界未大量向下位移。本文讨论了咸质含水组底部钙化粘土微观粘粒与水溶液的相互作用，提出了“滤盐层”的概念。     

多级汇水盆地的自动提取

利用多级汇水盆地在频率城（或波数域）中的频率差异，通过滤波方法逐级分离，并提取它们，实验证明，该方法效果较好。