PROSPECTING EXPLORATION

正20141888 Feng Lingang(Geological Survey Institute of Inner Mongolia,Hohhot 010020,China)The Conversion between Maps of Different Coordinate Systems Based on MapGIS(Geophysical and Geochemical Exploration,ISSN1000-8918,CN11-1906/P,37(5),2013,p.934-936,2 illus.,2 tables,7 refs.) Key words:geodetic coordinates,digital cartography     

变异函数模型参数的非线性加权总体最小二乘法

针对球状模型和指数模型两种变异函数模型线性化后的系数矩阵具有非线性形式,系数矩阵元素含有随机误差等问题,该文使用非线性加权总体最小二乘法估计变异函数模型参数。以高程异常数据为例,利用变异函数值的点对数定权法和系数矩阵中距离值的方差-协方差传播律定权法迭代解算参数,并与最小二乘法,加权最小二乘法和非线性总体最小二乘法进行对比分析。实验结果表明:非线性加权总体最小二乘法能够得到更高精度的变异函数模型参数。     

应用Slepian局部谱方法解算中国大陆重力场球谐模型

根据经典的球谐函数方法,为满足正交化要求,观测数据需要覆盖整个球面,而对于地表局部测量数据,则无法应用球谐方法解算重力场模型。针对此问题,采用Slepian局部谱分析方法解算中国大陆范围内的实测重力场变化数据,并以GOCE卫星球谐函数解作为已知模型,评估由于实际陆地重力测点的非均匀分布对球谐函数解的误差影响。通过计算多个阶次中国大陆局部范围的Slepian基函数分布;采用GOCE卫星获得重力场模型的前72阶球谐系数作为已知结果,评价实际测点非均匀分布的解算有效性,并针对中国大陆地区采用Slepian基函数进行解算,通过模型对比选择最优截段项数;针对2005—2008年中国大陆地区流动重力测量获得的重力场变化信号进行解算,获得了72阶重力场变化模型。     

扰动重力梯度的球冠谐分析建模

从球冠谐理论出发,详细推导了球冠坐标系下扰动重力梯度的无奇异性计算公式。基于Tikhonov正则化方法,利用GOCE卫星实际观测数据解算局部重力场球冠谐模型。数值计算表明,基于扰动重力梯度的球冠谐分析建模方法能够有效地恢复局部重力场中的短波信号,与GO_CONS_GCF_2_DIR_R5模型的差异在±0.3×10~(-5) m/s~2水平。     

A discrete global grid system for earth system modeling

To support Earth system modeling, we propose a discrete global grid system that expresses multi-resolution spatial data. Specifically, a unified coding model that expresses a grid of nodes, edges, and cells is constructed for a triangular discrete global grid system. To fulfill the requirements of practical applications, we design a code-based topological query method for this grid system and an algorithm to transform between grid codes and geographic coordinates. We evaluate the Global Finite Volume Community Ocean Model (Global-FVCOM) on the triangular discrete global grid system in the proposed uniform coding model. The ocean tidal waves simulated by the Global-FVCOM running on the coded grid are then compared with results obtained using a traditional irregular spherical grid system, and the results display comparable accuracy. The uniform coding model proposed in this paper provides a triangular discrete global grid system that can represent multi-resolution spatial data and can be used in Earth system models. This unified coding model can also be applied to the geographic coordinate system made up of latitudes and longitudes, as well as diamond and hexagonal grids.     

球谐分析及合成并行计算

本文介绍了球谐合成和分析问题的数值计算原理,研究分析了其内在的并行性,提出了区域划分和聚合通信的并行计算策略。基于傅立叶变换技术,应用MPI(消息传递接口)并行编程模式,实现了球谐分析及合成计算的并行化。实验结果表明,较之串行实现,并行实现大大缩短了计算时间,获得了近似线形加速比,取得很高的并行效率。     

根据三斜距确定点的三维坐标及精度

本文绘出一种根据三斜距解算点的三维坐标的直接方法,无需初值、迭代及多值判断,非常适合计算机编程;给出的实用精度估算式能方便地实现根据三斜里估算待定点三维点位精度之目的;对其应用条件的讨论,提出了应用中的具体倾角限值,对指导实际工作具有重要意义。     

GPS网与常规网三维联合平差的一种简便模型

提出了GPS控制网与常规控制网三维联合平差时分别固定平面坐标与高程的简便方法,从而简化了三维联合平差的模型,并给出算例验证了模型。     

C-RESPONSE OF GEOMAGNETIC DEPTH SOUNDING ON A 1D THIN SHELL MODEL

This paper tries to formulate the C-response of geomagnetic depth sounding(GDS)on an Earth model with finite electrical conductivity. The computation is performed in a spherical coordinate system. The Earth is divided into a series of thin spherical shells. The source is approximated by a single spherical harmonic P₁⁰ due to the spatial structure of electrical currents in the magnetosphere. The whole solution space is separated into inner and external parts by the Earth surface. Omitting displacement current, the magnetic field in the external space obeys Laplacian equation, while in the inner part, due to the finite conductivity, the electromagnetic fields obey Helmholtz equation. To connect the magnetic fields in the inner and external space, the continuity condition of magnetic fields is used on the Earth surface. The external magnetic fields are expressed by the inner and external source coefficients, from which a new parameter called C-response is computed from the inner coefficient divided by the external coefficient, thus normalizing the actual source strength. The inner magnetic fields in each layer can be recursively derived by the continuity boundary condition of both normal and tangential components of the magnetic field from the initial boundary condition at core-mantle-boundary. The consistency of our C-responses with that from a typical 1-D global model validates the accuracy of the proposed algorithm. Numerical results also show that the C-response estimated from the geomagnetic transfer function method will deviate exceeding 5%from the actual response at longer periods than about 10⁶s, which means that ignoring the curvature of the Earth at extreme long periods will make inversion result unreliable. Therefore, an accurate C-response should be computed in order to lay a solid foundation for reliable inversion.