利用球体剖分瓦块构建真三维数字地球平台

全球离散网格是面向空间大数据的模型框架,常用于构建数字地球平台。基于球体的剖分瓦块不仅可以构建真三维的数字地球模型,而且可以实现天地一体化的空间数据集成、融合、表达和应用。本文详细论述了球体大圆弧QTM八叉树网格的剖分原理、网格几何特征分析和编解码方法等理论体系,并利用剖分瓦块实现了球体的任意分割以及地下、地表和空中实体的可视化建模。研究表明,球体QTM网格具有剖分规则简单、体系规整、几何特征明晰,适用性强等特点,尤其是可以推广到椭球。因而,该方案可用于天地一体化的空间数据的组织、管理与应用。     

球面离散网格与球体离散网格的比较研究

测地线在球体剖分中扮演着十分重要的角色,它保证了球体剖分网格单元的形体简单性和空间关系的一致性。本文基于对欧式空间和流空间的数学认识,探讨了球体剖分的评价标准,提出了基于大圆弧QTM的球体八叉树剖分方法,并指出作为一种大尺度的集成三维空间数据模型,全球三维离散网格是地理信息科学的一个新的发展方向。     

球面四元三角网的复杂拓扑关系计算

针对已有球面四元三角网拓扑关系计算方法不能处理带空洞复杂目标的情况,提出了基于公共三角格网来计算球面四元三角网复杂拓扑关系,该方法能准确判断出带空洞复杂目标的拓扑关系,为全球空间分析和应用提供了新的思路。     

地图数据缩编更新中的道路自动选取研究

总结道路选取的研究现状,分析在道路的缩编更新过程中,道路要素选取需要考虑的要素指标。介绍一种道路选取方法,综合考虑道路的等级、长度、网眼面积、区域道路密度、是否是悬挂链等几何、语义和拓扑信息,模拟人类在缩编更新时的思维过程,实现道路的自动选取,并进行试验,验证方法的有效性。     

An Adaptive Nonhydrostatic Atmospheric Dynamical Core Using a Multi-Moment Constrained Finite Volume Method

An adaptive 2 D nonhydrostatic dynamical core is proposed by using the multi-moment constrained finite-volume(MCV) scheme and the Berger-Oliger adaptive mesh refinement(AMR) algorithm. The MCV scheme takes several pointwise values within each computational cell as the predicted variables to build high-order schemes based on single-cell reconstruction. Two types of moments, such as the volume-integrated average(VIA) and point value(PV), are defined as constraint conditions to derive the updating formulations of the unknowns, and the constraint condition on VIA guarantees the rigorous conservation of the proposed model. In this study, the MCV scheme is implemented on a height-based, terrainfollowing grid with variable resolution to solve the nonhydrostatic governing equations of atmospheric dynamics. The AMR grid of Berger-Oliger consists of several groups of blocks with different resolutions, where the MCV model developed on a fixed structured mesh can be used directly. Numerical formulations are designed to implement the coarsefine interpolation and the flux correction for properly exchanging the solution information among different blocks. Widely used benchmark tests are carried out to evaluate the proposed model. The numerical experiments on uniform and AMR grids indicate that the adaptive model has promising potential for improving computational efficiency without losing accuracy.     

Stream tube integration method for analysing subsurface fluid flow

A stream tube integration method is introduced to solve transient subsurface fluid flow problems. The method combines a geometry-embedded form of Darcy's Law and the notion of location of average. Two types of problems, transient radial flow to a well of finite radius in an areally infinite aquifer and in a double porosity system, are solved by the stream tube integration method and the integral finite difference method. Results of the solutions show that the stream tube integration method, with fixed coarse mesh, are more accurate and better behaved than the integral finite difference method, with fine mesh. The fixed mesh stream tube integration method is readily extended to the moving mesh method. With much coarse mesh, the moving mesh technique can obtain the same accurate results as the fixed mesh stream tube integration method. It is suggested that the stream tube integration method is a viable way to state, solve, interpret and verify numerical solutions. The method provides efficient computation and improved accuracy for analysing subsurface fluid flow. © 1998 John Wiley & Sons, Ltd.     

Poynting and polarization vectors based wavefield decomposition and their application on elastic reverse time migration in 2D transversely isotropic media

With the progress in computational power and seismic acquisition, elastic reverse time migration is becoming increasingly feasible and helpful in characterizing the physical properties of subsurface structures. To achieve high-resolution seismic imaging using elastic reverse time migration, it is necessary to separate the compressional (P-wave) and shear (S-wave) waves for both isotropic and anisotropic media. In elastic isotropic media, the conventional method for wave-mode separation is to use the divergence and curl operators. However, in anisotropic media, the polarization direction of P waves is not exactly parallel to the direction of wave propagation. Also, the polarization direction of S-waves is not totally perpendicular to the direction of wave propagation. For this reason, the conventional divergence and curl operators show poor performance in anisotropic media. Moreover, conventional methods only perform well in the space domain of regular grids, and they are not suitable for elastic numerical simulation algorithms based on non-regular grids. Besides, these methods distort the original wavefield by taking spatial derivatives. In this case, a new anisotropic wave-mode separation scheme is developed using Poynting vectors. This scheme can be performed in the angle domain by constructing the relationship between group and polarization angles of different wave modes. Also, it is performed pointwise, independent of adjacent space points, suitable for parallel computing. Moreover, there is no need to correct the changes in phase and amplitude caused by the derivative operators. By using this scheme, the anisotropic elastic reverse time migration is more efficiently performed on the unstructured mesh. The effectiveness of our scheme is verified by several numerical examples.     

Comparison of numerical solution and semi-empirical formulas to predict the effects of important design parameters on porpoising region of a planing vessel

Stability of the marine vessels in different conditions is one of the most important problems in the design of a planing vessel. In this research, the effects of some important design parameters (mass, longitudinal center of mass, deadrise angle, and length) of DTMB 62 model 4667-1 planing hull on the drag and also on the longitudinal dynamic stability (porpoising) are investigated in the velocity range of 2.12–8.486 m/s in calm water. In this paper, both numerical simulation of Reynolds Average Naiver Stokes (RANS) equations and semi-empirical formulas of Savitsky are used to analyze the motion of a 4667 planing vessel in calm water with two degrees of freedom (2DOF). For this purpose a finite volume, ANSYS-FLUENT, code is used to solve the Navier-Stokes equations for the simulation of the flow field around the vessel. In addition, an explicit VOF scheme and SST- Kω model is used with dynamic mesh scheme to capture the interface of a two-phase flow and to model the turbulence respectively, in 2DOF model (heave and pitch). Also, the results of both methods are compared with each other. According to the present results, changing the aspect ratio of the vessel and also the longitudinal center of gravity have the most effect on the porpoising region.     

电导率各向异性的海洋电磁三维有限单元法正演

本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究.     

Multi‐method virtual electromagnetic experiments for developing suitable monitoring designs: A fictitious CO2 sequestration scenario in Northern Germany

We present a numerical study for 3D time‐lapse electromagnetic monitoring of a fictitious CO₂ sequestration using the geometry of a real geological site and a suite of suitable electromagnetic methods with different source/receiver configurations and different sensitivity patterns. All available geological information is processed and directly implemented into the computational domain, which is discretized by unstructured tetrahedral grids. We thus demonstrate the performance capability of our numerical simulation techniques. The scenario considers a CO₂ injection in approximately 1100 m depth. The expected changes in conductivity were inferred from preceding laboratory measurements. A resistive anomaly is caused within the conductive brines of the undisturbed reservoir horizon. The resistive nature of the anomaly is enhanced by the CO₂ dissolution regime, which prevails in the high‐salinity environment. Due to the physicochemical properties of CO₂, the affected portion of the subsurface is laterally widespread but very thin. We combine controlled‐source electromagnetics, borehole transient electromagnetics, and the direct‐current resistivity method to perform a virtual experiment with the aim of scrutinizing a set of source/receiver configurations with respect to coverage, resolution, and detectability of the anomalous CO₂ plume prior to the field survey. Our simulation studies are carried out using the 3D codes developed in our working group. They are all based on linear and higher order Lagrange and Nédélec finite‐element formulations on unstructured grids, providing the necessary flexibility with respect to the complex real‐world geometry. We provide different strategies for addressing the accuracy of numerical simulations in the case of arbitrary structures. The presented computations demonstrate the expected great advantage of positioning transmitters or receivers close to the target. For direct‐current geoelectrics, 50% change in electric potential may be detected even at the Earth's surface. Monitoring with inductive methods is also promising. For a well‐positioned surface transmitter, more than 10% difference in the vertical electric field is predicted for a receiver located 200 m above the target. Our borehole transient electromagnetics results demonstrate that traditional transient electromagnetics with a vertical magnetic dipole source is not well suited for monitoring a thin horizontal resistive target. This is due to the mainly horizontal current system, which is induced by a vertical magnetic dipole.