共查询到20条相似文献,搜索用时 140 毫秒
1.
M. E. Everett 《Geophysical Journal International》1997,130(1):193-200
A new spherical mesh generator is described. It represents an efficient, deterministic packing of tetrahedra into a solid sphere, a spherical shell, or both. The mesh can be used for finite-element solutions to a wide variety of global numerical modelling problems in the geosciences. The nodes within the mesh are distributed uniformly, and long, thin tetrahedra are avoided. The method proposed here offers several advantages over 3-D Delaunay algorithms for finite-element mesh generation. For the related problem of trivariate scattered data interpolation, which is not considered here, the 3-D Delaunay algorithms are the method of choice. 相似文献
2.
Traction image method for irregular free surface boundaries in finite difference seismic wave simulation 总被引:16,自引:1,他引:16
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. 相似文献
3.
基于GIS的二维非结构化剖分网格优化 总被引:2,自引:1,他引:1
非结构化网格剖分是数值模拟的关键技术之一,网格质量直接影响到计算的收敛性和精确度。在GIS辅助建立非结构化网格空间拓扑关系的基础上,针对GIS和实际研究问题给出Spring-Laplace方法——一种新的单元尺度函数定义,在GIS空间算法下利用该方法优化节点位置,并基于推进阵面算法的思想,结合空间邻近拓扑关系实现了三角剖分节点和网格的重新编号算法,方便了开边界条件的赋值,提高了计算效率。实例表明,该方法大大提高了网格生成质量,能适应FVCOM数值模型对非结构化网格剖分的要求,其收敛速度快,具有较高的运算效率。 相似文献
4.
基于站点观测数据的气温空间化方法评述 总被引:1,自引:0,他引:1
基于统计学的插值方法是地理学、生态学领域研究气温空间化的主要方法之一,对获取精细化气温数据进行生态模拟具有重要意义。结合国内外气温空间插值的主要研究成果,对常用气温空间化方法进行了归纳、对比,探讨各种方法的适用性和不足之处,从而为涉及气温空间化的具体研究提供一定的参考,并探讨了各类方法优化的方向。不同方法的对比分析结果表明:各种气温空间化方法各有所长,在具体的应用中都取得过较好的效果,但并不存在普适性的方法,在实际应用时必须针对研究区域具体的地理特征进行方法适用性验证或对各类方法中的具体参数进行改进,才能实现区域气温的空间最优化模拟。根据气温场的物理分布特征,结合GIS技术,考虑地形等更多的相关因子以提高气温分布微观细节的模拟精度是未来重要的发展趋势。 相似文献
5.
Today's numerical methods like the Spectral Element Method (SEM) allow accurate simulation of the whole seismic field in complex 3-D geological media. However, the accuracy of such a method requires physical discontinuities to be matched by mesh interfaces. In many realistic earth models, the design of such a mesh is difficult and quite ineffective in terms of numerical cost. In this paper, we address a limited aspect of this problem: an earth model with a thin shallow layer below the free surface in which the elastic and density properties are different from the rest of the medium and in which rapid vertical variations are allowed. We only consider here smooth lateral variations of the thickness and elastic properties of the shallow layer. In the limit of a shallow layer thickness very small compared to the smallest wavelength of the wavefield, by resorting to a second order matching asymptotic approximation, the thin layer can be replaced by a vertically smooth effective medium without discontinuities together with a specific Dirichlet to Neumann (DtN) surface boundary condition. Such a formulation allows to accurately take into account complex thin shallow structures within the SEM without the classical mesh design and time step constraints. Corrections at receivers and source—when the source is located within the thin shallow layer—have been also derived. Accuracy and efficiency of this formulation are assessed on academic tests. The stability and limitations of this formulation are also discussed. 相似文献
6.
7.
Two-dimensional magnetotelluric inversion 总被引:2,自引:0,他引:2
Summary. When complex structure is encountered in magnetotelluric surveys, interpretation by locally fitted layered models is of questionable validity. However, when the processed data show two-dimensional structure, numerical inversion schemes for two-dimensional models may be constructed as an aid to regional data interpretation.
The two-dimensional magnetotellurics inversion problem is here formulated in a way that may be applied to many problems. A resulting computer program is analysed carefully in terms of its cost relative to that of simpler layered modelling.
As an example, the method is applied to some field data where the interpretive advantages of the program become evident. 相似文献
The two-dimensional magnetotellurics inversion problem is here formulated in a way that may be applied to many problems. A resulting computer program is analysed carefully in terms of its cost relative to that of simpler layered modelling.
As an example, the method is applied to some field data where the interpretive advantages of the program become evident. 相似文献
8.
针对基于全天空极光图像的极光事件自动分类问题,提出一种基于方向能量二元编码重组表征的自动分类方法。首先,通过对多个方向上能量分解来描述极光事件中的局部纹理和各个方向上的运动信息,并且结合分块策略获得极光事件的全局形态信息;然后,借鉴一种二元编码重组的方式对多个方向能量进行融合,从而使得极光事件的表征具有同时表征局部纹理、全局形态和运动信息的能力。该表征方法完全不依赖于极光事件的长度,可用于表征不同持续时间的极光事件,并且不需要复杂的训练过程。利用最近邻和支撑向量机分类器分别对从中国北极黄河站拍摄到的极光图像中挑选的特定极光事件进行自动分类,结果表明,与其他两种典型的动态纹理描述方法相比,本文所提出的表征方法结合最近邻分类器,得到了最好的分类效果,能有效用于极光事件的分析,为海量数据中的极光事件自动分类提供了一种新方法。 相似文献
9.
10.
11.
GIS技术支持下的洪水模型建模 总被引:11,自引:4,他引:11
在复杂区域建立洪水模型时,计算网格的手工生成方法容易出错甚至不可行,自动生成算法则可大大节省计算网格生成的工作量。洪水模型中的计算网格与GIS栅格数据及不规则三角网空间数据结构非常相似,因此,GIS中成熟的网格自动生成算法可用于生成洪水模型计算网格。文章详细讨论了GIS支持下的洪水模型自动建立步骤,并以黄河下游花园口~夹河滩河段为例,利用地形图、土地利用图、水利工程设施分布、水文站点图等资料,通过自动生成网格及其空间拓扑关系,建立了洪水过程数值模拟模型,并详细解释了计算网格数据格式。 相似文献
12.
Patrik Ottoson Hans Hauska 《International journal of geographical information science》2013,27(3):213-226
Systems for landscape visualization and geographical data handling require methods for efficient data access. Retrieval of data from large geographical databases, ten to thousands of Gbytes, is usually optimized with spatial indexing mechanisms. The simplest form of spatial indexing is achieved by dividing the database into congruent grid cells. The subsequent subdivision of the grid cells can be based on so-called quadtrees. Quadtrees for two-dimensional division and subdivision are appropriate for cartographical data. A geographical database, with objects stored in geocentric or geodetic (geographical) co-ordinates, requires indexing mechanisms that take into account the shape of the Earth. In this paper, we present a method for indexing of geographical data, named Ellipsoidal Quadtrees (EQT). In contrast to other global indexing methods, EQT is based on the Earth ellipsoid and not a spherical approximation. EQT division and subdivision make it possible to divide the Earth surface into a mesh of quadrangles with equal areas. We will demonstrate that EQT is flexible. It can be used for indexing databases of various sizes, including national and global databases. Tests on real data show that the performance of EQT is good. 相似文献
13.
区域尺度蒸散发遥感估算——反演与数据同化研究进展 总被引:3,自引:0,他引:3
遥感技术近年来在估算区域尺度蒸散发中应用广泛。不同方法在驱动数据、模型机理和适用范围往往存在很大差别。鉴于此,阐述了基于传统方法空间尺度扩展的遥感模型,经验统计公式,特征空间法,单源、双源垂向能量平衡余项法等几类的遥感蒸散发反演方法,简要介绍了三温模型、非参数化模型、半经验模型、集成模型等常用模型。同时,分析了遥感数据同化实现连续估算区域蒸散发的主要思路,综述了基于能量平衡和基于复杂过程模型的数据同化的原理、方法演进及常用同化算法等。最后,探讨了各类区域蒸散发遥感方法的优劣、展望了模型机理完善、不确定性研究、结果验证等与蒸散发直接反演和数据同化相关的研究方向。 相似文献
14.
流域水文模型计算域离散方法 总被引:7,自引:2,他引:7
常用的概念性水文模型 ,能够很好地模拟水文时间变化过程 ,但没有考虑水文变量和水文参数的空间变化与空间不均匀性。随着空间数据的获取手段的增多以及空间离散技术的发展 ,考虑水文参数和水文变量空间变化的分布式水文模型得到了极大的发展。本文详细介绍了分布式流域水文模型中用到的几种不同计算域离散方法 ,并讨论了河道汇流模型中常用到的有结构网格和无结构离散网格。地理信息系统技术对计算域离散有辅助作用 ,其有利于无结构离散网格的自动生成和交互修改 ,并可结合遥感技术 ,使水文模型能获取精确的空间分布的水文参数和水文变量。 相似文献
15.
We present a new numerical method to describe the internal dynamics of planetary mantles through the coupling of a dynamic model with the prediction of geoid and surface topography. Our tool is based on the simulation of thermal convection with variable viscosity in a spherical shell with a finite-volume formulation. The grid mesh is based on the 'cubed sphere' technique that divides the shell into six identical blocks. An investigation of various numerical advection schemes is proposed: we opted for a high-resolution, flux-limiter method. Benchmarks of thermal convection are then presented on steady-state tetrahedral and cubic solutions and time-dependent cases with a good agreement with the few recent programs developed to solve this problem.
A dimensionless framework is proposed for the calculation of geoid and topography introducing two dimensionless numbers: such a formulation provides a good basis for the systematic study of the geoid and surface dynamic topography associated to the convection calculations. The evaluation of geoid and surface dynamic topography from the gridded data is performed in the spectral domain. The flow solver is then tested extensively against a precise spectral program, producing response functions for geoid as well as bottom and surface topographies. For a grid mesh of a reasonable size (6 × 64 × 64 × 64) a very good agreement (to within ∼1 per cent) is found up to spherical harmonic degree 15. 相似文献
A dimensionless framework is proposed for the calculation of geoid and topography introducing two dimensionless numbers: such a formulation provides a good basis for the systematic study of the geoid and surface dynamic topography associated to the convection calculations. The evaluation of geoid and surface dynamic topography from the gridded data is performed in the spectral domain. The flow solver is then tested extensively against a precise spectral program, producing response functions for geoid as well as bottom and surface topographies. For a grid mesh of a reasonable size (6 × 64 × 64 × 64) a very good agreement (to within ∼1 per cent) is found up to spherical harmonic degree 15. 相似文献
16.
17.
Dengkui Li 《International journal of geographical information science》2018,32(9):1860-1883
Semi-parametric geographically weighted generalized linear models (S-GWGLMs) are a useful tool in modeling a regression relationship where the impact of certain explanatory variables on a non-Gaussian distributed response variable is global while that of others is spatially varying. In this article, we propose for S-GWGLMs a new estimation method, called two-stage geographically weighted maximum likelihood estimation, and further develop a likelihood ratio statistic-based bootstrap test to determine constant coefficients in the models. The performance of the estimation and test methods is then evaluated by simulations. The results show that the proposed estimation method performs as well as the existing method in estimating both constant and spatially varying coefficients but it is more efficient in terms of computation time; the bootstrap test is of accurate size under the null hypothesis and satisfactory power in identifying spatially varying coefficients. A real-world data set is finally analyzed to demonstrate the application of the proposed estimation and test methods. 相似文献
18.
Local Indicator of Colocation Quotient with a Statistical Significance Test: Examining Spatial Association of Crime and Facilities 总被引:6,自引:0,他引:6
Most existing point-based colocation methods are global measures (e.g., join count statistic, cross K function, and global colocation quotient). Most recently, a local indicator such as the local colocation quotient has been proposed to capture the variability of colocation across areas. Our research advances this line of work by developing a simulation-based statistical test for the local indicator of colocation quotient (LCLQ). The study applies the indicator to examine the association of land use facilities with crime patterns. Moreover, we use the street network distance in addition to the traditional Euclidean distance in defining neighbors because human activities (including facilities and crimes) usually occur along a street network. The method is applied to analyze the colocation of three types of crimes and three categories of facilities in a city in Jiangsu Province, China. The findings demonstrate the value of the proposed method in colocation analysis of crime and facilities and, in general, colocation analysis of point data. 相似文献
19.
水文循环模拟中下垫面参数化方法综述 总被引:2,自引:0,他引:2
针对水文循环模拟中地形、土地利用覆被等流域下垫面参数化方法众多,且模拟效果相差较大的现状。本文首先根据水文循环模拟中产汇流原理,对常用水文循环模拟中产汇流模拟方法进行汇总和分类;在此基础上,对产流模拟中的降水径流相关系数法、蓄满产流和超渗产流等及汇流模拟中的等流时线、单位线、圣维南方程、马斯京根法等主要模拟方法中地形、土地利用覆被和土壤类型参数化方法进行分析和讨论;根据其中流域地形、土地利用覆被和土壤类型参数化方法对机理过程的描述程度,将其分为无明确表示类、率定型参数类、确定型参数类、物理过程表达类;进而阐明不同参数化方法中流域地形、土地利用覆被和土壤类型对水文循环模拟结果的响应和贡献。最后回归模型本质,阐述水文循环模拟中流域下垫面参数化方法中存在经验关系对复杂机理简单表述的合理性和物理机理过程描述的欠缺性问题,并预估未来水文循环模拟中下垫面参数化方法朝着简洁实用化和复杂机理化两个方向发展。 相似文献
20.
黑河下游典型植被下垫面与大气间能量传输模拟研究 总被引:4,自引:3,他引:1
在分析过去地-气间相互作用的物理过程、研究进展基础上,探讨了土壤、植被、大气系统中水分与能量的传输过程在模型中不同的计算方法。同时对不同计算方法进行对比,确定了分析土壤、植被、大气系统中水分和能量传输方法的优缺点。利用黑河下游额济纳地区绿洲试验区2003年9月的大气资料作为陆面模式的强迫场,研究陆面过程模式(LSM)在极端干旱地区的模拟能力。模拟结果表明,在观测资料的强迫下,LSM能够较好地模拟出地表特征量的变化趋势。根据实际情况,定义额济纳地区绿洲的植被覆盖率为0.7, 叶冠高度为1.5 m,位移高度为0.64 h,Karman常数取0.4。通过对比试验发现,采用LSM模型模拟的空气温度和地面蒸发比实际蒸发少,但地表潜热通量、辐射和土壤热通量的模拟结果与实测值的比较吻合。如果地表和植被的参数选择较好,采用LSM模型模拟本地区的陆面过程将有较好的结果。 相似文献