Comparison of methods to model the gravitational gradients from topographic data bases

A number of methods have been developed over the last few decades to model the gravitational gradients using digital elevation data. All methods are based on second-order derivatives of the Newtonian mass integral for the gravitational potential. Foremost are algorithms that divide the topographic masses into prisms or more general polyhedra and sum the corresponding gradient contributions. Other methods are designed for computational speed and make use of the fast Fourier transform (FFT), require a regular rectangular grid of data, and yield gradients on the entire grid, but only at constant altitude. We add to these the ordinary numerical integration (in horizontal coordinates) of the gradient integrals. In total we compare two prism, two FFT and two ordinary numerical integration methods using 1" elevation data in two topographic regimes (rough and moderate terrain). Prism methods depend on the type of finite elements that are generated with the elevation data; in particular, alternative triangulations can yield significant differences in the gradients (up to tens of Eötvös). The FFT methods depend on a series development of the topographic heights, requiring terms up to 14th order in rough terrain; and, one popular method has significant bias errors (e.g. 13 Eötvös in the vertical–vertical gradient) embedded in its practical realization. The straightforward numerical integrations, whether on a rectangular or triangulated grid, yield sub-Eötvös differences in the gradients when compared to the other methods (except near the edges of the integration area) and they are as efficient computationally as the finite element methods.     

Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.     

Drainage spacing regularity on a fault-block: A case study from the eastern Ruahine Range

Abstract:  Recent research has indicated river basin outlets draining linear sections of large, uplifting mountain belts often show a regularity of spacing, transverse to the main structural trend. A morphometric analysis of part of the Ruahine Range, on the North Island was undertaken to test whether drainage regularity may exist in smaller, younger mountain ranges. The ratio, R , of the half-width of the mountain belt, W , and the outlet spacing, S , was used to characterize drainage networks on the eastern side of the range. The spacing ratio for the range of 1.31 is lower than R results from studies of larger mountain belts ( R  = 1.91–2.23). We suggest the cause of this lower ratio is related to eastward migration of the Ruahine drainage divide.     

油气化探数据库系统

建立油气化探数据库是油气化探信息资源共享、科学管理、宏观决策的基础和依据,是油气化探数据处理及异常评价的基础。本系统采用动态自适应结构,收集、录入了几十年来我国油气化探数据,可以完成油气化探数据登录、管理、浏览、检索、报表、数据预处理。文章中论述了油气化探数据库的结构、特点及功能,反映了油气化探数据库研究领域的新进展和研究水平。     

青海湖地区生态环境动态变化遥感监测

由于青海湖地区的生态环境较为脆弱，且人类活动进一步加剧，人口、资源与环境的矛盾日渐突出，因此，近年来．青海湖及其周边地区的生态环境出现了明显变化，主要表现在水位下降及水域面积减小、草原退化、沙质荒漠化土地面积扩大等。文章采用1975年MSS卫星图像及1987、2000年TM卫星图像作为遥感信息源，并结合地理信息系统方法，旨在查明青海满地区耕地、沙质荒漠化土地和水域等生态环境要素的时空演化规律，为青海湖地区实现资源开发与环境协调发展提供科学依据。监测结果表明，25a来，青海湖地区的耕地及沙质荒漠化土地面积出现明显的扩大，而水域面积出现明显缩小，同时由于湖周各河流土壤侵蚀的加剧，在部分河流入湖处泥沙淤积较为严重，生态环境出现明显恶化。     

基于工作流技术的动态GIS服务链研究

提出了一种基于工作流技术的动态GIS服务链的实现方法,该方法使用抽象GIS服务组合服务链,由GIS服务动态选择机制将抽象服务映射到具体服务,并结合工作流技术,将工作流描述语言描述的GIS服务链流程送入工作流引擎,由工作流引擎管理和监控服务链的运行。     

地震作用下的灰坝液化特征及其动力稳定性分析--以安阳电厂为例

根据安阳电厂灰坝岩土工程性质及其动力学特性,建立了适合灰坝工程的动力分析模型并采用有限元法进行了动力分析.在此基础上,在不同工况条件下对灰坝进行了抗液化安全评价及抗震稳定性分析.研究表明,在不设碎石桩、无排渗体条件下,粉煤灰子坝的抗液化安全系数Ks＜1.25,将发生液化;在设碎石桩、有排渗体条件下,粉煤灰子坝的抗液化安全系数明显提高,Ks≥1.25,不会发生液化.抗震稳定性分析表明,在上述两种工况条件下灰坝是稳定的.     

高峰公司100号矿体矿床建模与矿山生产动态管理的研究

重点论述高峰锡矿 10 0号矿体的矿床数学经济模型的建立 ,分析模型在采矿、选矿等矿山生产动态管理中的作用及意义     

对永胜6级地震的预报及依据

2001年10月27日在云南永胜发生6.0级地震，震前四川省地震局有较准确的中期和短临预报，加强了川滇交界地区的地震监测预报工作，并且向四川省政府作了汇报，取得了明显的社会效益。经验性地震预报重视川滇菱形地块西界断裂带动态活动图像的分析，尤其对历史强震破裂空段附近出现的中小震活动带或新活跃区段的研究；同时，重视4级前震群、川滇交界地区地下水位、水温、形变观测异常的跟踪分析和震情预测。