Comparison of methods to model the gravitational gradients from topographic data bases   总被引：7，自引：0，他引：7  

Christopher Jekeli  Lizhi Zhu 《Geophysical Journal International》2006,166(3):999-1014

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   总被引：16，自引：1，他引：16  

Wei Zhang  Xiaofei Chen 《Geophysical Journal International》2006,167(1):337-353

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     

Heather Purdie  Martin Brook 《New Zealand geographer》2006,62(2):97-104

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.  相似文献   

基于预测的边缘检测方法   总被引：3，自引：0，他引：3  

王晖  卢健 《武汉大学学报(信息科学版)》2005,30(1):23-26

提出了一种基于预测的、有一定自适应性的边缘检测方法。介绍了其原理和算法实现过程。通过计算两个相邻像素的灰度值的均值、均方差和梯度来预测下一个像素的灰度值，比较预测值和真实值来判断下一个像素是否是边界点。通过与现有同类算法的比较，证明了该算法的优越性。  相似文献   

玛纳斯河一级电站引水渠沿线提取地下水进行融冰的可行性研究     

王春得 《地下水》2005,27(5):347-349

通过对向斜洼地含水层水文地质特征的勘探及研究,对研究区的水文地质边界条件进行了概化,并用水均衡法估算了地下水资源量,对拟建一级电站引水渠沿线提取地下水进行融冰的可行性进行了研究,并提出了取水方案.  相似文献   

层序的测井、地震响应特征研究   总被引：3，自引：5，他引：3  

田景春  张翔  聂永生  夏青松  赵强  韦东晓  李强  李国忠  吴琦 《物探化探计算技术》2005,27(4):321-326

在层序地层学研究中,关键是层序划分和对比。而层序划分、对比的关键是层序识别。层序的识别包括层序界面(层序的底界面、初始海泛面和最大海泛面)的识别和构成层序的体系域识别。这里,在众多前人研究成果的基础上,详细研究了层序的测井、地震响应特征,即层序界面和体系域在测井曲线上和地震剖面上的特征。其中,层序底界面在测井上表现为突变的钟型、箱型或侧积式曲线的底界;在地震剖面上表现为剥蚀、顶超、上超、下超;而体系域在测井曲线上的响应为:低水位体系域的海底扇以漏斗形中、高幅的前积式,或钟型中、低幅的后积式模式为特征,陆坡扇成钟型、正向齿形,自下而上幅度由中高幅→低幅,即具后积式测井模式,低水位楔的测井曲线表现为旋回性进积模式特征,其特征表现为锯齿状箱型。海侵体系域的测井曲线呈现向上变细、变深序列,并表现为钟型、正向齿形或齿化状,幅度由高幅变化为低幅,包络线具后积式特征。高水位体系域相应的测井曲线呈现中幅箱形或桶形,不同体系域在地震剖面上的响应特征明显不同。  相似文献   

豫陕交界地区钾长石矿产地质特征及开发前景     

张斌  陈全树  何文平  杨福新 《西北地质》2005,38(2):82-85

本文研究证明区内有南、北两条钾长石成矿带，其产出形态、矿化特征、矿石类型均有差异，且导致了它们在工业应用上的不同。区内矿体以岩脉或岩株产出，有红色、浅红色、白色3种矿石。白色矿石富钾低铁，达到国标一级品，其余矿石为二级以下品级，均具有一定的销售市场。从全国的供需形势分析，今后钾长石矿产品市场前景广阔，而区内钾长石资源进一步找矿的潜力仍然很大。  相似文献   

地球表层系统浅论     

张猛刚  雷祥义 《西北地质》2005,38(2):99-101

地球表层系统是岩土圈、大气圈、水圈、人类圈和生物圈构成的地表自然社会综合体，笔者从系统论的观点出发，阐述了地球表层系统的边界、环境、要素、结构和功能；并分析了地球表层系统的开放性、有序性和耗散结构的特征；最后阐明了协调人地关系是地球表层系统研究的主要内容。  相似文献   

北京夏季城市热岛现状及楔形绿地规划对缓解城市热岛的作用   总被引：20，自引：2，他引：20       下载免费PDF全文

佟华  刘辉志李延明桑建国胡非 《应用气象学报》2005,16(3):357-366

通过遥感技术与地面测定相结合的方法，对北京城市热岛现状作观测研究，得到北京城市地面的温度分布特点。使用北京大学城市边界层模式从气象观点就“楔形绿地”规划对北京城市气候的影响进行研究和评价，模式通过对城市地表复杂性和多样性的特征进行细致描述，建立了一个细致模拟城市特点的城市边界层能量平衡模式，并用此能量平衡模式得到的地面温度作为下边界条件，中尺度气象模式MM5做初始场和侧边界条件，建立一个最小分辨率为500 m的城市边界层模式系统，来研究城市边界层在中尺度背景场作用下的精细结构。通过个例模拟，模式能够较准确地模拟城市边界层的风温场分布情况，可以用来对楔形绿地规划进行模拟试验。通过对规划后的气象场在特定的气象条件下进行模拟，结果显示，建造大型的楔形绿地后，绿地区域及绿地周围约1 km以内的地区温度有所降低，降低的程度由规划前后的地表类型改变的剧烈程度、风速大小及与绿地的距离决定，但是这种规划方案却会因城市的下风方向的风速减小而导致通风不畅。  相似文献   

大气边界层湍流温度序列的信息熵分析   总被引：2，自引：0，他引：2  

汪景烨  付遵涛  张霖  刘式达 《高原气象》2005,24(1):38-42

利用大气边界层内近地面的大气湍流温度时间序列，运用功率谱分析、信息熵分析等方法，分析了大气边界层内近地面的大气湍流特点，并对稳定层和不稳定层的大气湍流进行了对比。结果表明，信息熵和功率谱指数是区别稳定层结和不稳定层结大气边界层湍流特征的指标，对造成两者之间的差别做出了对应的解释。  相似文献