Spherical Slepian functions and the polar gap in geodesy

The estimation of potential fields such as the gravitational or magnetic potential at the surface of a spherical planet from noisy observations taken at an altitude over an incomplete portion of the globe is a classic example of an ill-posed inverse problem. We show that this potential-field estimation problem has deep-seated connections to Slepian's spatiospectral localization problem which seeks bandlimited spherical functions whose energy is optimally concentrated in some closed portion of the unit sphere. This allows us to formulate an alternative solution to the traditional damped least-squares spherical harmonic approach in geodesy, whereby the source field is now expanded in a truncated Slepian function basis set. We discuss the relative performance of both methods with regard to standard statistical measures such as bias, variance and mean squared error, and pay special attention to the algorithmic efficiency of computing the Slepian functions on the region complementary to the axisymmetric polar gap characteristic of satellite surveys. The ease, speed, and accuracy of our method make the use of spherical Slepian functions in earth and planetary geodesy practical.     

VLBI资料处理软件系统OCCAM V6.0简介

OCCAM V6．0是在PC机上开发的一套VLBI大地测量数据处理软件。该软件用标准FORTRAN77(90)语言编写，可在Windows95以上版本和LINUX环境下的PC机上运行，其主要特点是操作简便、兼容性好、处理容量大、速度快和结果输出多样化等。随着该软件的不断改进和发展，利用OCCAM进行VLBI数据分析的学者逐渐增多，OCCAM软件的使用也日益广泛。本文简要介绍了OCCAM软件的安装、各个模块的主要功能和在运行中的注意事项，供今后OCCAM用户在使用过程中参考。     

少量地面控制点的遥感图像精纠正与定位技术探讨

本文主要讨论了利用空间投影一空间斜墨卡托（SOM）投影进行卫星遥感图像精纠正与定位，解决缺少控制点地区卫星遥感图像精纠正与定位问题，并根据正反算公式讨论了卫星遥感图像纠正的空间投影数学基础。     

Numerical evaluation of the damping‐solvent extraction method in the frequency domain

The damping‐solvent extraction method for the analysis of unbounded visco‐elastic media is evaluated numerically in the frequency domain in order to investigate the influence of the computational parameters—domain size, amount of artificial damping, and mesh density—on the accuracy of results. An analytical estimate of this influence is presented, and specific questions regarding the influence of the parameters on the results are answered using the analytical estimate and numerical results for two classical problems: the rigid strip and rigid disc footings on a visco‐elastic half‐space with constant hysteretic material damping. As the domain size is increased, the results become more accurate only at lower frequencies, but are essentially unaffected at higher frequencies. Choosing the domain size to ensure that the static stiffness is computed accurately leads to an unnecessarily large domain for analysis at higher frequencies. The results improve by increasing artificial damping but at a slower rate as the total (material plus artificial) damping ratio ζ_t gets closer to 0.866. However, the results do not deteriorate significantly for the larger amounts of artificial damping, suggesting that ζ_t≈0.6 is appropriate; a larger value is not likely to influence the accuracy of results. Presented results do not support the earlier suggestion that similar accuracy can be achieved by a large bounded domain with small damping or by a small domain with larger damping. Copyright © 2002 John Wiley & Sons, Ltd.     

渤海湾盆地埕岛东斜坡地区东三段油气成岩成藏模式

通过分析埕岛东斜坡地区东三段砂岩成岩作用特征 ,以及对该地区油气成藏系统的研究 ,总结出东三段砂岩体的成岩成藏模式。东三段上部砂体成藏于馆陶组沉积期 ,对应于早成岩 B期阶段 ,油气过早进入砂岩体 ,充填了孔隙 ,抑制了成岩作用的进行 ,使原生孔隙得以良好的保存。目前该砂体位于晚成岩 A期阶段 ,伴随形成了次生孔隙 ,并正处于第二次成藏期     

Combination of Dependent Realizations Within the Gradual Deformation Method

Gradual deformation is a parameterization method that reduces considerably the unknown parameter space of stochastic models. This method can be used in an iterative optimization procedure for constraining stochastic simulations to data that are complex, nonanalytical functions of the simulated variables. This method is based on the fact that linear combinations of multi-Gaussian random functions remain multi-Gaussian random functions. During the past few years, we developed the gradual deformation method by combining independent realizations. This paper investigates another alternative: the combination of dependent realizations. One of our motivations for combining dependent realizations was to improve the numerical stability of the gradual deformation method. Because of limitations both in the size of simulation grids and in the precision of simulation algorithms, numerical realizations of a stochastic model are never perfectly independent. It was shown that the accumulation of very small dependence between realizations might result in significant structural drift from the initial stochastic model. From the combination of random functions whose covariance and cross-covariance are proportional to each other, we derived a new formulation of the gradual deformation method that can explicitly take into account the numerical dependence between realizations. This new formulation allows us to reduce the structural deterioration during the iterative optimization. The problem of combining dependent realizations also arises when deforming conditional realizations of a stochastic model. As opposed to the combination of independent realizations, combining conditional realizations avoids the additional conditioning step during the optimization process. However, this procedure is limited to global deformations with fixed structural parameters.     

半无限空间中圆形孔洞周围SH波的散射

建立了求解在含有圆形孔洞的弹性半空间中SH波散射与圆形孔洞附近动应力集中问题的解析方法。利用SH波散射的对称性和多极坐标的方法，构造了一个可以预先满足半空间自由表面上应力自由边界条件的圆形孔洞对SH波散射的波函数。利用这一波函数，则可将该问题转化成对一个圆形孔洞散射的求解问题。该问题的解答最终又可归结为对一组无究代数方程组的求解问题，并可利用截断有限项的方法对其进行计算，最后给出了有关圆形也洞附近动应力集中问题的算例和数值结果，并讨论了波数与圆孔至自由边界距离变化对动应力集中的影响。     

The newest observational evidence on asymmetrical deformation of the Earth

IntroductionWhat is the shape of the Earth? Does it change continuously? It is a scientific question since the ancient times and is still being observed and explored at present. In 250 BC, Greek scholar Eratosthene supposed the shape of the Earth to be spherical according to the observations to the Sun and estimated the perimeter of the Earth to be 4 000 km (King-Hele, 1976) according to the camel-walking distance. Until the 16th century, the Earth was considered to be a very symmetrical …     

Relation of the Wave-Propagation Metric Tensor to the Curvatures of the Slowness and Ray-Velocity Surfaces

The contravariant components of the wave-propagation metric tensor equal half the second-order partial derivatives of the selected eigenvalue of the Christoffel matrix with respect to the slowness-vector components. The relations of the wave-propagation metric tensor to the curvature matrix and Gaussian curvature of the slowness surface and to the curvature matrix and Gaussian curvature of the ray-velocity surface are demonstrated with the help of ray-centred coordinates.     

城市群演化的空间过程及土地利用优化配置

城市群演化包含了诸多内容 ,如人口、产业、城市类型、城市职能等 ,而空间过程是最直接、最综合的表现。