共查询到20条相似文献，搜索用时 531 毫秒
1.
为了解决三维复杂介质的射线追踪,本文改变了传统的三维层状地层的建模描述方式,提出了块状结构的建模描述方法,结合三角形面片来描述地质界面,可以构造非常复杂的三维地质模型．为了满足射线追踪的需要,本文对模型界面内的法向量进行光滑处理,光滑后的法向量在界面内是连续变化的．在块状模型的基础上,本文运用三角形的面积坐标,提出了几种试射角度的修正方法：细分三角形法、分割三角形法和子三角形法,计算表明子三角形法最好．文中给出了三维块状模型和射线追踪实例． 相似文献
2.
沥青混凝土是由骨料、沥青胶浆、空气按照一定的体积百分比混合而成的多相非匀质混合物,其骨料、沥青胶浆和空气的体积不等、形状各异、介电特性不同、空间位置随机分布,具有明显的多相、离散、随机介质特征.本文基于随机介质模型理论,(1)测量与统计了介电常数在典型沥青混凝土芯样空间上的随机分布统计特征;(2)估算了沥青混凝土介质的自相关函数及其特征参数(自相关长度、自相关角度等),确定其随机介质类型;(3)提出了量化约束下的多相离散随机介质建模算法,以混合型椭圆自相关函数为基础,构建了不同粗糙度因子的多相离散随机介质模型;(4)构建了不同空隙率的多相离散随机介质模型,正演模拟与对比分析了探地雷达波在均匀介质、连续型随机介质和多相离散随机介质中的传播特征.结果表明:多相离散随机介质模型不仅描述了沥青混凝土的多相、离散与空间随机分布统计特征,而且进一步描述了其各组成物质体积百分比,能更全面、准确地描述沥青混凝土的介质特征,同时也为描述其他类似材料或介质提供了新的方法和途径;在多相离散随机介质模型中,探地雷达波散射强烈,随机、无序传播的散射波相互叠加干涉,形成了明显的随机扰动和"噪声",致使异常体反射波扭曲变形、不连续,降低了探地雷达回波的信噪比和分辨率.研究探地雷达波的随机扰动特征与多相离散随机介质模型参数之间的关系,将为定量评价多相离散随机介质的属性参数提供参考和帮助. 相似文献
3.
4.
基于交错网格伪谱法和高阶精度有限差分方法,发展了模拟非均匀介质地震波传播的三维伪谱和有限差分混合算法.该方法在两个水平方向利用交错网格伪谱算子计算空间微分,保留了该方法高效、高精度的优势,在垂直方向采用交错网格高阶精度有限差分算子实现空间微分计算．利用有限差分方法的局部性特征,将三维计算区域在垂直方向上划分为一系列子区域,并分配给不同的处理器,实现了在并行计算机集群上的三维并行计算．通过模拟算例,与离散波数法比较,检验了该算法的精度．为了检验该方法的实用性,在64个处理器上,对三维沉积盆地模型进行了67108864个网格点的并行计算,模拟的波场主频率为1.25Hz,讨论了沉积盆地深度对三维沉积盆地地面运动的影响． 相似文献
5.
三维瞬变电磁正演算法作为研究处理与解释方法的重要基础,如何加速计算过程,减少计算机内存消耗尤为重要.为此,本文采用有限体积算法在八叉树(octree)网格上对时间域Maxwell方程组进行空间离散,相比于交错六面体网格,octree网格在局部区域网格细度相同的条件下可以显著减少细化区域外的网格数量,对复杂几何体边界的模拟更加灵活,而相比于非结构四面体网格,octree网格单元位置排列更加规律.通过octree网格离散三维正演模型使计算网格规模显著降低,减少了待求解方程的未知数,降低了物理内存消耗.空间离散后,瞬变电磁正演响应可以表示为关于初始磁场的矩阵指数函数,采用位移逆Krylov子空间模型降阶算法实现瞬变电磁场的求解,只需对系数矩阵进行一次矩阵分解和多次回代即可获得一系列时间序列的瞬变电磁场正演结果.本文算法在空间离散和方程求解两方面优化三维瞬变电磁正演过程,数值算例结果验证了本文算法的精度和高效性. 相似文献
6.
介绍一种基于离散粒子理论地震波传播数值模拟的网格剖分计算方法.根据离散粒子理论,将研究区域划分为由一系列相互作用的粒子组成的正六边形网格,这些粒子在它们的接触点处发生相互作用,并用Hooke定律和Newton定律描述.为解决六边形网格带来的网格交错而难以计算以及波场输出问题,将横向网格进行加密,加密处赋予假想的粒子,输出波场时选取偶数行偶数列点或奇数行奇数列点的波场值.均匀介质和层状介质模型的数值模拟结果表明,该网格剖分计算方法能够将离散粒子理论用于模拟弹性波在非均匀各向同性介质中地震波的传播. 相似文献
7.
一种新的三维速度结构反演模型参数化方法 总被引：3，自引：1，他引：3
本文针对三维速度反演问题中的一个重要步骤--模型的参数化作了探讨和研究,提出了一种新的参数化方法.同以往采用的矩形分块或划分网格点的方法不同,将待求扰动场按其空间频率展开,反演各阶频率系数.从三维傅氏级数理论出发,导出了一个简明且便于计算的三维空间频谱展开公式,并从理论上分析了在空间域内离散模型和空间频率域内离散模型的不同之处. 理论分析表明,在空间域内采用过参数化离散模型来反演,会在所得到的解中混入人为的虚假结构,从而使反演解发生畸变.采用在频率域内离散模型的频谱参数化方法反演,可以有效地避免这一点.频谱参数化方法也使得我们可以过参数化（over-parameterization）离散模型,便于最大程度地从资料中提取正确信息,对模型的分辨程度完全取决于观测系统的分辨能力.这种方法在得到稳定的反演解的同时,还可得到观测系统在空间三个方向上的分辨率.数值模拟实验结果证实了以上结论. 相似文献
8.
地震射线追踪是地震定位、层析成像、偏移等领域的重要正演环节.随着这些领域研究的深入,针对传统的网格结构和层状结构在描述复杂地质模型遇到的很大困难,我们采用大小不等、形状各异的地质块组成的集合体来描述三维复杂地质模型,并用三角形面片来描述地质块之间的物性间断面,理论上可以描述任意复杂的地质模型.为适应任意非均匀速度分布的地质模型,基于费马原理,本文发展了与之相适应的逐段迭代射线追踪方法.该方法属于弯曲法范畴,对路径点采用一阶显式增量修正,相对于传统的迭代法,高效省时.数值试验表明,联合逐段迭代法和伪弯曲法的射线追踪扰动修正方案在三维复杂非均匀块状模型中有适用性和高效性. 相似文献
9.
10.
研究了三维各向异性裂缝介质的正演模拟的三种不同有限差分法,即标准交错网格(SSG)、旋转交错网格(RSG)和Lebedev方法(LS),详细对比分析这三种交错网格方法在模拟复杂各向异性介质的优势与局限.提出一种新的改进方法,简化了LS有限差分法在任意各向异性介质中的正演模拟.为了模拟三维大规模复杂各向异性介质的地震响应,提出一个优化的正演模拟计算流程：将模型参数分为模型介质参数和模型构造参数.该计算流程适用于三种有限差分法中的任何一种.使用LS方法实现任意三维各向异性裂缝介质地震响应的三维全波场模拟,通过使用三种不同的有限差分模拟方法进行二维和三维模型数值模拟试验,验证了所提出方法有效. 相似文献
11.
随着数字化地震台网的普及和利用以及计算机分析处理技术的提高，近几年利用地震波研究地幔内部间断面有了一定的进展。本文概述了近年来国内外有关学者对地幔内部间断面的认识，采用的研究方法及研究动态。 相似文献
12.
Kerstin Ericson 《地球表面变化过程与地形》2004,29(4):495-509
Granite landforms in Sierra Nevada, California, and the southern part of Sweden were used to test whether a Schmidt hammer can be used to distinguish surfaces of different age and origin. Based on 3260 readings, statistically signi?cant differences were obtained from surfaces formed and/or affected by different geomorphological processes. Aeolian, ?uvial or glacially polished surfaces yield higher values than surfaces exposed to surface weathering, which are distinguishable from surfaces at a weathering front caused by deep weathering. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
13.
探讨了如何基于Midas/GTS来实现用FLAC3D建立复杂边坡模型。首先在Midas/GTS中建立几何模型和划分网格,然后根据两种软件的网格数据形式进行转换,把Midas/GTS的网格模型导入FLAC3D中。对于二维边坡问题,认为采用四边形单元比三角形单元的计算精度更高,在模型较为复杂的情况下建议在Midas/GTS中采用四边形+三角形的方式生成混合网格,然后导入FLAC3D中形成brick+wedge形式的混合网格进行分析,以保证计算精度。最后将上述建模方法应用于某边坡工程的稳定性分析,验证了该建模方法的可行性和实用性。 相似文献
14.
Abstract. A simple-to-use computer program (FLOWNS) has been developed for generating flow nets for any saturated rectangular domain with any combination of constant head or constant flux (including zero) boundary conditions. The program approximates with discrete values the continuous distributions of potential and stream function using finite-difference approximations of Laplace's equation. The hydraulic conductivity distribution may be anisotropic and/or heterogeneous. A contouring program is required to generate the final stream and equipotential lines. 相似文献
15.
采用有限元法精确模拟复杂介质条件下的地震波场,一般采用三角单元,但在节点数相同的情况下,三角单元的计算精度不如矩形单元高,采用三角单元模拟复杂界面时,编制确定结构刚度矩阵非零元素位置的程序也较为麻烦。采用矩形单元离散含有倾斜或起伏界面的地质模型时,无法避免绕射噪声,若加密网格又会增加计算量。为此,本文基于任意四边形单元模拟声波的传播,在倾斜或起伏界面条件下,可以有效避免因离散引起的“阶梯状”界面,在不增加计算量以及内存占用的前提下,有效地消除离散绕射噪声。采用对角的集中质量矩阵代替一致质量矩阵,避免矩阵的求逆运算,从而提高显式有限元法的计算效率;对结构刚度矩阵采用紧凑存储格式,每一行需要存储的元素最多为5个,同时零元素不参与运算,既减少内存的占用,又极大地提高计算效率。 相似文献
16.
Existing analytical procedures for nonsteady flow in a leaky confined aquifer assume that the aquifer system is areally infinite. A technique is presented that treats a leaky confined aquifer system of finite configuration. By means of a discrete space continuous time (DSCT) modeling approach, the partial differential equation governing the flow system is transformed into a set of ordinary differential equations that can be easily integrated numerically on a high speed digital computer using available scientific subroutines. The finite difference formulation is in effect an explicit scheme. A criterion is developed for which the scheme is computationally stable. A numerical example is presented. 相似文献
17.
Automated recognition of quasi‐planar ignimbrite sheets as paleosurfaces via robust segmentation of digital elevation models: an example from the Central Andes 
下载免费PDF全文
下载免费PDF全文 Balázs Székely Zsófia Koma Dávid Karátson Peter Dorninger Gerhard Wörner Melanie Brandmeier Clemens Nothegger 《地球表面变化过程与地形》2014,39(10):1386-1399
Quasi‐planar morphological surfaces may become dissected or degraded with time, but still retain original features related to their geologic‐geomorphic origin. To decipher the information hidden in the relief, recognition of such features is required, possibly in an automated manner. In our study, using Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM), an existing algorithm has been adapted to recognize quasi‐planar features fulfilling specified criteria. The method has been applied to a study area of the Central Andes with Miocene to Quaternary volcanic edifices, tilted ignimbrite surfaces, and basin‐filling sediments. The result is a surface segmentation, whereas non‐planar features (gullies, tectonic faults, etc.) are sorted out. The main types of geomorphic features that can be distinguished and interpreted are as follows. (1) The west‐dipping western margin of the Altiplano is differentiated into segments of the lower sedimentary cover that of increased erosion by tectonic steepening at intermediate levels, and an upper plane with limited erosion. (2) In the central part of the Western Cordillera, the Oxaya ignimbrite block shows a ‘striped’ bulging pattern that results from a smoothly changing surface dip. This pattern is due to continuous folding/warping of the ignimbrite block possibly related to gravitational movements. (3) To the west, large, uniform planes correspond to flat, smooth, tectonically undisturbed surfaces of young sedimentary cover of the Central Basin. (4) The evolution of Taapaca volcanoes with sector collapse events and cone‐building phases is shown by several segments with overlapping clastic aprons. (5) To the east, on the western margin of the Altiplano, young intermontane basins filled by Upper Miocene sediments show progressively increasing dip toward basin margins, reflected by a circular pattern of the segmentation planes. We show that the segmentation models provide meaningful images and additional information for geomorphometric analysis that can be interpreted in terms of geological and surface evolution models. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
18.
Elena Zhebel Sara Minisini Alexey Kononov Wim A. Mulder 《Geophysical Prospecting》2014,62(5):1111-1125
The finite‐difference method on rectangular meshes is widely used for time‐domain modelling of the wave equation. It is relatively easy to implement high‐order spatial discretization schemes and parallelization. Also, the method is computationally efficient. However, the use of finite elements on tetrahedral unstructured meshes is more accurate in complex geometries near sharp interfaces. We compared the standard eighth‐order finite‐difference method to fourth‐order continuous mass‐lumped finite elements in terms of accuracy and computational cost. The results show that, for simple models like a cube with constant density and velocity, the finite‐difference method outperforms the finite‐element method by at least an order of magnitude. Outside the application area of rectangular meshes, i.e., for a model with interior complexity and topography well described by tetrahedra, however, finite‐element methods are about two orders of magnitude faster than finite‐difference methods, for a given accuracy. 相似文献
19.
《Advances in water resources》2005,28(7):729-744
In this paper, we discuss the local discontinuous Galerkin (LDG) method applied to elliptic flow problems and give details on its implementation, focusing specifically on the case of piecewise linear approximating functions. The LDG method is one a family of discontinuous Galerkin (DG) methods proposed for diffusion models. These DG methods allow for very general hp finite element meshes, and produce locally conservative fluxes which can be used in coupling flow with transport. The drawback to DG methods, when compared to their continuous counterparts, is the number of degrees of freedom required to compute the solution. This motivates a coupled approach, discussed herein, where the solution is allowed to be continuous or discontinuous on a node-by-node basis. This coupled approximation is locally conservative in regions where the numerical solution is discontinuous. Numerical results for fully discontinuous, continuous and coupled discontinuous/continuous solutions are given, where we compare solution accuracy, matrix condition numbers and mass balance errors for the various approaches. 相似文献
20.
Hydrogeologists often are called upon to estimate surfaces from discrete, sparse data points. This estimation is often accomplished by manually drawing contours on maps using interpolation methods between points of known value while accounting for features known to influence the water table's surface. By contrast, geographic information systems (GIS) are good at creating smooth continuous surfaces from limited data points and allowing the user to represent the resulting surface resulting with contours, but these automated methods often fail to meet the expectations of many hydrogeologists because they do not include knowledge of other influences on the water table. In this study, we seek to fill this gap in the GIS‐based methodology for hydrogeologists through an interactive tool that shapes an interpolated surface based on additional knowledge of the water table inferred from gaining or losing streams. The modified surface is reflected in water table contours that, for example, “V” upstream for gaining streams, and can be interactively adjusted to fit the user's expectations. By modifying not only the contours but also the associated interpolated surface, additional contours will follow the same trend, and the modified surface can be used for other analyses like calculating average gradients and flow paths. The tool leverages Esri's ArcGIS Desktop software, building upon a robust suite of mapping tools. We see this as a prototype for other tools that could be developed for hydrogeologists to account for variations in the water table inferred from local topographic trends, pumping or injection wells, and other hydrogeologic features. 相似文献