基于环形扫面测量的三维直流电阻率法任意各向异性模型响应特征     

殷长春  杨志龙  刘云鹤  张博  齐彦福  曹晓月  邱长凯  蔡晶 《吉林大学学报(地球科学版)》2018,48(3):872-880

高精度地下任意各向异性介质电性分布特征识别是当前的热门研究课题。本文以国内外前人工作为基础,基于非结构化网格,利用基于梯度恢复的后验误差估计指导网格自适应细化过程,实现了直流电阻率法三维有限元数值模拟。通过与一维模型半解析解的对比,验证了本文算法的有效性。考虑到电各向异性介质在观测中存在视电阻率反常现象,本文采用了环形扫面测量方法。通过对几种典型各向异性模型的模拟分析,得到了相应的各向异性影响规律和识别特征,各向异性主轴电阻率之间的比值大小决定了椭圆型视电阻率极性曲线长轴与短轴的比值大小,主轴电阻率的旋转方向会改变视电阻率极性曲线的形状。本文的算法研究及数值模拟技术可为直流电法数据精细处理和解释提供技术支撑。  相似文献   

Solving the linearized forward-speed radiation problem using a high-order finite difference method on overlapping grids     

《Applied Ocean Research》2017

The linearized potential flow approximation for the forward speed radiation problem is solved in the time domain using a high-order finite difference method. The finite-difference discretization is developed on overlapping, curvilinear body-fitted grids. To ensure numerical stability, the convective derivatives in the free-surface boundary conditions are treated using an upwind-biased stencil. Instead of solving for the radiation impulse response functions, a pseudo-impulsive Gaussian type displacement is employed in order to tailor the frequency-content to the discrete spatial resolution. Frequency-domain results are then obtained from a Fourier transform of the force and motion signals. In order to make a robust Fourier transform, and capture the response around the critical frequency, the tail of the force signal is asymptotically extrapolated assuming a linear decay rate. Fourth-order convergence of the calculations on simple geometries is demonstrated, along with a nearly linear scaling of the solution effort with increasing grid resolution. The code is validated by comparison with analytical and semi-analytical solutions using submerged and floating closed-form geometries. Calculations are also made for a modern bulk carrier, and good agreement is found with experimental measurements.  相似文献   

电导率各向异性的海洋电磁三维有限单元法正演   总被引：10，自引：8，他引：2       下载免费PDF全文

蔡红柱  熊彬  Michael Zhdanov 《地球物理学报》2015,58(8):2839-2850

本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究.  相似文献   

Multi‐method virtual electromagnetic experiments for developing suitable monitoring designs: A fictitious CO2 sequestration scenario in Northern Germany          下载免费PDF全文

Jana H. Börner  Feiyan Wang  Julia Weißflog  Matthias Bär  Ines Görz  Klaus Spitzer 《Geophysical Prospecting》2015,63(6):1430-1449

We present a numerical study for 3D time‐lapse electromagnetic monitoring of a fictitious CO₂ sequestration using the geometry of a real geological site and a suite of suitable electromagnetic methods with different source/receiver configurations and different sensitivity patterns. All available geological information is processed and directly implemented into the computational domain, which is discretized by unstructured tetrahedral grids. We thus demonstrate the performance capability of our numerical simulation techniques. The scenario considers a CO₂ injection in approximately 1100 m depth. The expected changes in conductivity were inferred from preceding laboratory measurements. A resistive anomaly is caused within the conductive brines of the undisturbed reservoir horizon. The resistive nature of the anomaly is enhanced by the CO₂ dissolution regime, which prevails in the high‐salinity environment. Due to the physicochemical properties of CO₂, the affected portion of the subsurface is laterally widespread but very thin. We combine controlled‐source electromagnetics, borehole transient electromagnetics, and the direct‐current resistivity method to perform a virtual experiment with the aim of scrutinizing a set of source/receiver configurations with respect to coverage, resolution, and detectability of the anomalous CO₂ plume prior to the field survey. Our simulation studies are carried out using the 3D codes developed in our working group. They are all based on linear and higher order Lagrange and Nédélec finite‐element formulations on unstructured grids, providing the necessary flexibility with respect to the complex real‐world geometry. We provide different strategies for addressing the accuracy of numerical simulations in the case of arbitrary structures. The presented computations demonstrate the expected great advantage of positioning transmitters or receivers close to the target. For direct‐current geoelectrics, 50% change in electric potential may be detected even at the Earth's surface. Monitoring with inductive methods is also promising. For a well‐positioned surface transmitter, more than 10% difference in the vertical electric field is predicted for a receiver located 200 m above the target. Our borehole transient electromagnetics results demonstrate that traditional transient electromagnetics with a vertical magnetic dipole source is not well suited for monitoring a thin horizontal resistive target. This is due to the mainly horizontal current system, which is induced by a vertical magnetic dipole.  相似文献   

油藏数值模拟中地质模型的建模流程与方法     

张剑波  李谢清  石阳  朱建勃 《吉林大学学报(地球科学版)》2015,45(3):860-868

围绕油藏数值模拟过程中三维地质模型的建模技术进行了分析与探讨,详细介绍了基于角点网格模型的建模方法,给出了相应的实现步骤。其主要流程是:首先根据断层数据构造断层模型,在断层模型的基础上构建骨架模型;然后在骨架模型约束下采用地层恢复技术实现含断层的地层模型;最后基于结构模型插值物性参数完成属性模型。以塔河油田缝洞型油藏为例,对建模流程和技术的可行性进行验证,结果表明,其建模结果与专业地质建模软件Petrel相符。  相似文献   

顾及空间自相关的地理国情信息统计格网尺度选择——以植被覆盖信息统计为例   总被引：1，自引：0，他引：1  

连世忠  丁霖  陈江平 《测绘通报》2017,(3):46

统计格网尺度的不同会带来统计结果的差异,如何选择统计格网是地理国情信息统计的重要工作。本文提出了一种顾及空间自相关的地理国情信息统计格网尺度选择方法。采用地理国情普查数据,在50 m、60 m、70 m、80 m、90 m、100 m、250 m、500 m和1000 m几个尺度下,以植被覆盖信息统计为例,利用面积占优法和中心点归属法两种方法分别进行格网化,得到了不同尺度的植被格网数据;计算植被覆盖面积统计误差,分析不同尺度下植被覆盖信息的空间自相关的变化特征,并利用Moran’s I系数差值进行尺度选择,得到了植被覆盖信息统计格网的适宜尺度。以龙沙区和清涧县作为研究区域,结果表明,在地理国情植被覆盖信息统计时,不同地区的格网统计适宜尺度是不一样的,植被覆盖度中低的龙沙区的适宜尺度为100 m,而植被覆盖度高的清涧县的适宜尺度为250 m。  相似文献   

Higher-order compositional modeling with Fickian diffusion in unstructured and anisotropic media     

Joachim Moortgat  Abbas Firoozabadi 《Advances in water resources》2010

We present advances in compositional modeling of two-phase multi-component flow through highly complex porous media. Higher-order methods are used to approximate both mass transport and the velocity and pressure fields. We employ the Mixed Hybrid Finite Element (MHFE) method to simultaneously solve, to the same order, the pressure equation and Darcy's law for the velocity. The species balance equation is approximated by the discontinuous Galerkin (DG) approach, combined with a slope limiter. In this work we present an improved DG scheme where phase splitting is analyzed at all element vertices in the two-phase regions, rather than only as element averages. This approximation is higher-order than the commonly employed finite volume method and earlier DG approximations. The method reduces numerical dispersion, allowing for an accurate capture of shock fronts and lower dependence on mesh quality and orientation. Further new features are the extension to unstructured grids and support for arbitrary permeability tensors (allowing for both scalar heterogeneity, and shear anisotropy). The most important advancement in this work is the self-consistent modeling of two-phase multi-component Fickian diffusion. We present several numerical examples to illustrate the powerful features of our combined MHFE–dg method with respect to lower-order calculations, ranging from simple two component fluids to more challenging real problems regarding CO₂ injection into a vertical domain saturated with a multi-component petroleum fluid.  相似文献   

A three-dimensional model for suspended sediment transport based on the compact discontinuous Galerkin method     

《国际泥沙研究》2016,(1):36-43

  相似文献   

On generalised and specialised spatial information grids: are geo-services ready?     

《International Journal of Digital Earth》2013,6(4):315-325

Abstract

After the set-up of a spatial data infrastructure (SDI) and a national information infrastructure (NII) in many countries, the provision of geo-services became one of the most important and attractive tasks. With the integration of global positioning system (GPS), geographic information system (GIS) and remote sensing (RS), we can, in principle, answer any geo-spatial related question: when and where what object has which changes? An intelligent geo-service agent could provide end-users with the most necessary information in the shortest time and at the lowest cost. Unfortunately there is still a long way to go to achieve such goals. The central component in such geo-services is the integration of the spatial information system with a computing grid via wire- and wireless communication networks. This paper will mainly discuss the grid technology and its integration with spatial information technology, expounding potential problems and possible resolutions. A novel categorising of information grids in the context of geo-spatial information is proposed: generalised and specialised spatial information grids.  相似文献   

GIS科普资源发布系统实现的关键技术     

袁家勇  公茂玉  高诗亮 《测绘科学》2011,36(3):239-240,218

科普资源发布系统利用GIS强大的空间查询与分析功能,方便市民浏览、查询与统计科普资源情况,掌握更多科普知识.本文主要论述了系统设计与实现过程中应用的关键技术,详细给出了栅格影像自动选择的方法,矢量栅格匹配的原理,以及统计图的绘制技术.最后介绍了系统的开发模式与软件环境,给出了统计效果图.  相似文献