首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 156 毫秒
1.
采用易于模拟复杂地形起伏和倾斜界面的非结构三角单元剖分网格,并利用对偶加权后验误差估计指导网格自动细化过程,实现了2.5维直流电阻率法自适应有限元数值模拟。在实例模型分析中,分别计算了层状模型和垂直岩脉模型的直流电阻率响应,并与其解析解进行了比较。对比结果表明,该算法所得数值解精度很高,解的相对误差小于0.5%。最后,计算了起伏地形2.5维地电模型视电阻率异常,并利用比较法进行了地形改正。地形改正结果与水平地形时的结果对比表明,比较法可以较好地消除地形影响,突出局部地质体异常。  相似文献   

2.
起伏地形频率域可控源电磁二维快速正反演   总被引:1,自引:0,他引:1  
李静和 《地质与勘探》2018,54(2):325-331
起伏地形频率域可控源电磁正反演研究受到众多学者的关注和重视,目前普遍实用的为未考虑地形的二维正反演算法程序。本文基于积分方程法正演和多场源、多频率、对比源反演算法,开展起伏地形条件下频率域可控源电磁二维快速正反演研究。通过引入层状参考模型、将地形与目标体整体作为异常场剖分区域,实现起伏地形可控源电磁正反演计算。针对包含地形在内的大尺度剖分区域引起的大型计算代价问题,采用快速傅里叶算法提高正反演计算效率。通过与现有正反演算法进行模型算例对比,验证了本文采用方法的可行性与有效性。  相似文献   

3.
可地浸层间氧化带型砂岩铀矿自然电位异常计算   总被引:1,自引:0,他引:1  
理论和实践表明,在可地浸层间氧化带型砂岩型矿床上可观察到明显的自然电位异常。其起因不同于电子导电体形成的氧化还原电场,而在于发生在层间氧化带中一些非平衡现象之间相互影响所产生的混合电位。根据其自然极化原理,建立了等效数学模型,并对几种常见矿体形态的自然电位异常进行了模拟计算,并分析了曲线特征。最后,讨论了自然电位异常地形校正方法,并给出了几种常见地形条件下自然电位模拟曲线。  相似文献   

4.
基于计算机视觉的地形3D重建   总被引:1,自引:0,他引:1  
构造3D地形模型有很多种方法.本文提出一种用数码相机获得的未标定照片获取三维地形模型的简单方法.这种方法仅需要从两幅或多幅照片中选取对应点,根据这些对应点,求出基本矩阵,摄像机矩阵,进而进行三维地形表面的估计,实验结果及算法分析表明,本文提出的特征匹配算法精确度较高,实现地形重建方法切实可行,生成的模型可视效果较好。  相似文献   

5.
2.5维起伏地表条件下坐标变换法直流电场数值模拟   总被引:3,自引:1,他引:2  
研究起伏地表对视电阻率分布的影响是进行地形校正的基础。由于很难处理不规则边界,计算简洁且效率高的有限差分方法很少用于解决起伏地表问题。为了解决该问题,引入曲化平思想,实现了一种基于坐标变换法的起伏地表条件下的直流电场数值模拟方法。方法从传统的2.5维基本方程及边界条件出发,通过坐标变换将起伏地表问题转化为水平地表问题,并利用有限差分法求解水平地表问题,最后再通过坐标映射得到起伏地表条件下的电位和视电阻率分布。精度分析及计算实例表明:本方法实现简洁,计算误差主要集中在震源附近,整个计算区域内的平均相对误差为1.39%,计算结果满足起伏地表条件下的电位和视电阻率的分布规律。  相似文献   

6.
前陆盆地挠曲沉降和沉积过程3D模型研究   总被引:1,自引:0,他引:1  
胡明卿  刘少峰 《地质学报》2012,86(1):181-187
前陆盆地是在造山带负荷作用下岩石圈发生挠曲沉降而形成的,并且被主要从造山带搬运的沉积物所充填。为了更好地理解和认识前陆盆地的形成演化机制,特别是受控于周缘多个造山带活动所形成的前陆盆地的演化机制,本文通过建立前陆盆地挠曲沉降与沉积过程的3-D模型,模拟展示了造山带逆冲推覆作用、岩石圈挠曲沉降响应及在山盆体系中由于动力地形变化而导致的河流体系的发育变化及其产生的剥蚀和沉积过程。模型的建立和实验完整体现了逆冲推覆、弹性挠曲沉降和沉积物搬运这三者之间的耦合机制,为全面深入研究前陆盆地动力学提供了理论依据和方法。  相似文献   

7.
8.
This paper reports on the numerical modelling of flash flood propagation in urban areas after an excessive rainfall event or dam/dyke break wave. A two-dimensional (2-D) depth-averaged shallow-water model is used, with a refined grid of quadrilaterals and triangles for representing the urban area topography. The 2-D shallow-water equations are solved using the explicit second-order scheme that is adapted from MUSCL approach. Four applications are described to demonstrate the potential benefits and limits of 2-D modelling: (i) laboratory experimental dam-break wave in the presence of an isolated building; (ii) flash flood over a physical model of the urbanized Toce river valley in Italy; (iii) flash flood in October 1988 at the city of Nîmes (France) and (iv) dam-break flood in October 1982 at the town of Sumacárcel (Spain). Computed flow depths and velocities compare well with recorded data, although for the experimental study on dam-break wave some discrepancies are observed around buildings, where the flow is strongly 3-D in character. The numerical simulations show that the flow depths and flood wave celerity are significantly affected by the presence of buildings in comparison with the original floodplain. Further, this study confirms the importance of topography and roughness coefficient for flood propagation simulation.  相似文献   

9.
The accuracy (detail, resolution) of controlled-source electrical prospecting is a combination of several disparate elements. They include the model framework used in the interpretation. A real geologic medium has a complex structure. Sedimentary rocks have a layered structure with fractal properties: The layers split into smaller ones. Large-scale geoelectrical studies (for example, electrical prospecting) require a proper geoelectrical model. By necessity, the 1D geoelectrical model in electrical prospecting is horizontally layered, with thick (hundreds of meters) homogeneous layers, whereas a fine structure is neglected. To study some aspects of this problem, we performed a set of numerical experiments. They were aimed at studying the TEM response from a formation consisting of many thin layers with random geoelectric parameters, mainly resistivity.  相似文献   

10.
艾智勇  蔡建邦 《岩土力学》2015,36(Z2):685-688
将地基视为多层各向同性弹性体,对Euler-Bernoulli梁进行有限单元离散分析,对地基-梁接触面采用边界积分法求解,根据地基-梁接触面的竖向位移协调和光滑接触条件,应用有限元与边界元耦合的方法推导出各向同性成层弹性地基上的Euler-Bernoulli梁的半数值半解析解。基于文中理论,编制了相应的程序,通过与现有文献对比,验证了文中理论的正确性,对比分析了分层地基与等效均匀地基两种模型上的弹性地基梁。研究结果表明,分层地基与等效均匀地基两种模型上弹性地基梁性状差异较大,实际弹性地基梁计算中应采用分层地基模型。  相似文献   

11.
基于层面法向局部坐标系中的横观各向同性弹性本构模型,根据层面产状方位角,推导出整体坐标系中弹性应力的转换方法。针对层状岩体的各向异性特征,将岩体的破坏模式细分为5种,分别对每一种破坏模式建立各自的判别函数和屈服准则,根据塑性流动正交法则推导出相应的应力修正迭代计算方法。通过单轴、三轴压缩数值试验,对层状岩体的强度和变形性质进行研究,分析层面倾角和围压对层状岩体力学行为的影响。将所建模型与FLAC3D自带的横观各向同性弹性模型以及遍布节理模型进行对比,验证模型的正确性及精度。同时,将模型应用于三向等压理想圆形隧洞的开挖计算中,分析层状岩体破坏方式受层面方位角影响的变化规律,分析结果表明层面对围岩破坏模式和塑性区扩展方向起到了控制作用。  相似文献   

12.
13.
海岸带潮滩地形具有快速变化的特点;在盐沼植被覆盖的潮滩区域,地面激光扫描仪(terrestrial laser scanning,TLS)获取地形数据受到地表植被的影响。为了研究TLS在盐沼潮滩使用的地形测量精度与盐沼植被种类、盖度的关系,本文以芦苇群落、白茅群落、互花米草群落、海三棱藨草群落4种典型海岸带盐沼植被群落为研究对象,在移动窗口法的基础上辅以聚类分析的植被滤除算法,分别从点云原始数据中恢复地表地形地貌特征。研究结果表明:1)植被盖度越高,TLS反演地形精度越低,两者负相关。2)不同植被的激光穿透能力不同:盖度大于50%时,激光无法穿透原始盖度分别为70%、65%、65%的白茅群落、互花米草群落、海三棱藨草群落,均方根误差(root-mean-square error,RMSE)分别为22.0、22.0、8.6 cm;盖度等于50%时,白茅群落、海三棱藨草群落、芦苇群落和互花米草群落的RMSE分别为16.0、6.6、4.5、5.7 cm;盖度小于50%时,芦苇群落、互花米草群落、海三棱藨草群落地形反演精度小幅度提高,白茅群落地形反演精度提升较为明显。3)在盐沼潮滩地区使用TLS反演地形时,增加TLS架设高度、对同一区域多方位反复扫描可能有助于提高地形反演效果。  相似文献   

14.
长江口斜压诊断模式三维流场数值模拟   总被引:4,自引:0,他引:4       下载免费PDF全文
建立了σ坐标系下长江口斜压诊断模式三维流场数学模型,采用k-kl二方程紊流闭合模型求解垂向涡粘系数,计算域内恒定、非均匀盐度场反映了计算域密度斜压效应。验证结果表明,斜压诊断模式的模拟成果较切合实际。模型成功复演了长江口水域由往复流向旋转流过渡的流场特性;潮波从四条汊道传入后在分汊口附近相遇形成多个会潮点,其位置随着径流的大小、潮汐的强弱、汊道的形态等因素而变化。  相似文献   

15.

Satellite altimetry has been widely used to determine surface elevation changes in polar ice sheets. The original height measurements are irregularly distributed in space and time. Gridded surface elevation changes are commonly derived by repeat altimetry analysis (RAA) and subsequent spatial interpolation of height change estimates. This article assesses how methodological choices related to those two steps affect the accuracy of surface elevation changes, and how well this accuracy is represented by formal uncertainties. In a simulation environment resembling CryoSat-2 measurements acquired over a region in northeast Greenland between December 2010 and January 2014, different local topography modeling approaches and different cell sizes for RAA, and four interpolation approaches are tested. Among the simulated cases, the choice of either favorable or unfavorable RAA affects the accuracy of results by about a factor of 6, and the different accuracy levels are propagated into the results of interpolation. For RAA, correcting local topography by an external digital elevation model (DEM) is best, if a very precise DEM is available, which is not always the case. Yet the best DEM-independent local topography correction (nine-parameter model within a 3,000 m diameter cell) is comparable to the use of a perfect DEM, which exactly represents the ice sheet topography, on the same cell size. Interpolation by heterogeneous measurement-error-filtered kriging is significantly more accurate (on the order of 50% error reduction) than interpolation methods, which do not account for heterogeneous errors.

  相似文献   

16.
Ai  Zhi Yong  Ye  Zi  Song  Xiaoyu  Wang  Lu Jun 《Acta Geotechnica》2019,14(4):1143-1160

We develop a new numerical model based on a precise integration method to investigate the coupled thermo-mechanical performance of layered transversely isotropic media around a cylindrical/tubular heat source. To obtain the relational matrices of the extended precise integration method, we first convert the governing equations of the problem into ordinary differential matrix equations through the Laplace–Hankel transform. Then, the cylindrical heat source is divided into a series of plane heat sources, and the plane temperature load term is added to the state vector between layer elements. By combining the layer elements, we build a layered transversely isotropic numerical model containing a cylindrical heat source in the transformed domain. Finally, we solve the model in the transformed domain and obtain the solution of the problem in the real domain through the Laplace–Hankel transform inversion. The accuracy of this method is verified by comparing the solutions with the results of the analytical method and the finite element method. Then, we study the influence of the anisotropy of thermal parameters, the embedded depth, the length/radius ratio, the type of heat source and the stratification of the medium on the thermo-mechanical coupled performance.

  相似文献   

17.
董岩  谭捍东  付兴 《现代地质》2023,37(1):74-83
为了研究复杂地电模型的航空瞬变电磁法全波形响应特征,需要开发考虑发射波形的三维数值模拟算法。本研究基于非结构四面体网格和位移逆Krylov子空间(Shift-and-Invert Krylov,简称SAI Krylov)方法,采用基于电偶极子离散的场源处理方法模拟场源,在时间域进行计算实现了全波形航空瞬变电磁法矢量有限元三维数值模拟。使用均匀半空间模型在阶跃波、半正弦波、三角波和梯形波激发下的全波形解析解、VTEM实际激发波形的后推欧拉算法计算结果,检验了本研究开发的数值模拟算法的正确性。设计地表起伏异常体模型,计算和分析了航空瞬变电磁响应特征。开发的基于位移逆Krylov子空间的全波形航空瞬变电磁法三维数值模拟算法适合模拟复杂地电模型的响应,具有较高的计算精度。  相似文献   

18.
Flood basalts represent large outpourings of lavas which often cover, and interact with, sedimentary basins. For this reason areas with significant flood basalt cover are often targets for hydrocarbon exploration. Problems exist, however, when trying to image sediments and structures in offshore regions covered by basalts. Here we present preliminary 3-D models of the Etendeka flood basalt province from NW Namibia, which can act as an aid in understanding the internal and external architecture of the flood basalt cover. Satellite images, digital elevation models, measured geological logs, sections and maps, are used to create the geological model. Models are presented in 2 parts; 1) models created using topography with images such as Land Sat and geological maps draped over them, and 2) a 3-D model of key lava and sediment surfaces in the basin as defined by measured geological sections. Initial results show a palaeo-volcanic feature early in the flood basalt history which is onlapped by later Iavas. The modelling also allows a simple correction for post emplacement subsidence by assuming an original sub horizontal position for the first basin wide silicic flow unit.  相似文献   

19.
In this paper we advocate the use of Newton's law of gravitational attraction to ensure perfect consistency between gravity and height data. Starting with the absolute gravity on the topography we decompose this signal into a number of quantities associated with physics of the system. To model gravitational attraction from topography we use DTM and Newton's law of gravitational attraction. A residual part of the gravity signal is interpreted as inconsistency between gravity and heights. In the paper we discuss a method by which such inconsistency (at least in principle) can be decomposed into a “gravity error” and a “terrain error”. In practice such separation is not possible because the two types of error are nearly 100% correlated. The inconsistency can be interpreted as a measure of ambiguity of the gravity-terrain models which are consistent with a set of measured/interpolated data. We discuss the influence of such ambiguity on the accuracy of the geoid for the investigated area of Jutland, Denmark.  相似文献   

20.
山区河流平面二维流场的数值模拟   总被引:11,自引:0,他引:11       下载免费PDF全文
针对山区河流边界地形复杂、水位及流速变化剧烈的特点,用"混合五对角法"建立平面二维数学模型并求解,进行水位和流速的数值模拟.利用嘉陵江重庆金刚碑-朝阳桥河段地形及5个横断面的水位和水深平均流速资料对此模型进行了验证,模型快速、准确地再现了该河段的水位和水深平均流速情况.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号