Simulations of 3D DC borehole resistivity measurements with a goal-oriented <Emphasis Type="Italic">hp</Emphasis> finite-element method. Part II: through-casing resistivity instruments

We simulate direct current (DC) borehole resistivity measurements acquired in steel-cased deviated wells for the assessment of rock formation properties. The assumed data acquisition configuration considers one current (emitter) and three voltage (collector) electrodes that are utilized to measure the second difference of the electric potential along the well trajectory. We assume a homogeneous, 1.27-cm-thick steel casing with resistivity equal to 10^− 5 Ω· m. Simulations are performed with two different numerical methodologies. The first one is based on transferring two-dimensional (2D) axisymmetric optimal grids to a three-dimensional (3D) simulation software. The second one automatically produces optimal 3D grids yielded by a 3D self-adaptive goal-oriented algorithm. Both methodologies utilize high-order finite elements (FE) that are specially well-suited for problems with high-contrast coefficients and rapid spatial variations of the electric field, as it occurs in simulations that involve steel-cased wells. The method based on transferring 2D-optimal grids is efficient in terms of CPU time (few seconds per logging position). Unfortunately, it may produce inaccurate 3D simulations in deviated wells, even though the error remains below 1% for the axisymmetric (vertical) well. The method based on optimal 3D grids, although less efficient in terms of CPU time (few hours per logging position), produces more accurate results that are validated by a built-in a posteriori error estimator. This paper provides the first existing simulations of through-casing resistivity measurements in deviated wells. Simulated resistivity measurements indicate that, for a 30° deviated well, measurements in conductive layers 0.01 Ω· m) are similar to those obtained in vertical wells. However, in resistive layers (10,000 Ω· m), we observe 100% larger readings in the 30° deviated well. This difference becomes 3,000% for the case of a 60° deviated well. For this highly-deviated well, readings corresponding to the conductive formation layer are about 30% smaller in magnitude than those in a vertical well. Shoulder effects significantly vary in deviated wells.     

A coupling of mixed and discontinuous Galerkin finite-element methods for poroelasticity

In this paper, we formulate a finite-element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. We approximate the flow variables by a mixed finite-element space and the displacement by a family of discontinuous Galerkin methods. Theoretical convergence error estimates are derived and, in particular, are shown to be independent of the constrained specific storage coefficient, c _o . This suggests that our proposed algorithm is a potentially effective way to combat locking, or the nonphysical pressure oscillations, which sometimes arise in numerical algorithms for poroelasticity.     

The Glarus thrust: excursion guide and report of a field trip of the Swiss Tectonic Studies Group (Swiss Geological Society, 14.–16. 09. 2006)

This excursion guide results form a field trip to the Glarus nappe complex organized by the Swiss Tectonic Studies Group in 2006. The aim of the excursion was to discuss old and recent concepts related to the evolution of the Glarus thrust. The major aspects were (i) the interplay between deformation, fluid flow and geochemical alteration, (ii) episodic versus continuous deformation and fluid flow, and (iii) the link between large-scale structures, microstructures, and geochemical aspects. Despite 150 years of research in the Glarus nappe complex and the new results discussed during the excursion, there exist controversies that still are unsolved.     

西藏札达盆地控盆断裂有限元数值模拟

札达盆地是喜马拉雅构造带中的一个山间断陷盆地,其演化过程与盆地两侧的控盆构造密切相关。对控盆断裂的构造应力场进行模拟计算,将有助于进一步深化对本区构造控盆的认识。因此,在对盆地构造地质进行详细调查的基础上,结合本区的深部地质与地球物理资料,对札达盆地控盆断裂的构造应力场进行了模拟。计算结果表明,札达盆地的演化明显受盆地两侧边界断裂的控制,札达盆地是在整体南北向挤压应力的作用下,不同块体差异隆升作用的结果。其南侧的控盆断裂为北倾的正断层,北侧的控盆断裂为南倾的逆断层,二者共同形成了南降北升的翘板式断陷盆地运动过程,是喜马拉雅地块在陆内汇聚挤压构造环境中构造应力场调整的一种方式。     

利用DInSAR和GPS综合方法估计地表3维形变速率

在分析基于贝叶斯统计和马尔科夫随机场理论的解析优化法的基础上,提出以GPS提供的高精度的水平形变速率作为约束,利用直接分解法,将DInSAR的斜距向形变速率分解成垂直形变速率。试验结果表明,直接分解法有效地综合了GPS的高水平精度与DInSAR的高垂直精度优点,且计算简单,克服了解析优化法可能出现的数值不稳定问题。     

基于GPS技术监测高层建筑物变形的研究

以武汉大学信息学部教学实验大楼为实验对象,分析了大楼在N、E、U方向上的每h变化、d变化及变形规律,阐明GPS技术在高层建筑变形监测中的应用。     

钻孔灌注桩施工对邻近运营隧道变形影响研究

桩基施工将引起周边地层位移和邻近隧道结构变形,以南京龙津桥改建项目桩基工程为背景,通过现场实测,分析研究钻孔灌注桩施工全过程各工况条件对地层位移场和邻近既有隧道结构变形的影响规律。结果表明:由于钢套管的护壁作用,钻孔灌注桩施工过程中产生的最大地层位移和隧道结构变形较小,说明采用“钢套管边旋压边取土”、“群桩间跳施工”等工艺,对周围土体扰动影响程度较小,建议近隧桩基工程采用全套管灌注桩的施工方法,监测分析结果可为类似工程提供技术参考。     

Research on temporal and spatial distribution, evolutionary character and mechanism of crustal deformation field before and after the Tangshan earthquake

Ｒｅｓｅａｒｃｈｏｎｔｅｍｐｏｒａｌａｎｄｓｐａｔｉａｌｄｉｓｔｒｉｂｕ┐ｔｉｏｎ,ｅｖｏｌｕｔｉｏｎａｒｙｃｈａｒａｃｔｅｒａｎｄｍｅｃｈａｎｉｓｍｏｆｃｒｕｓｔａｌｄｅｆｏｒｍａｔｉｏｎｆｉｅｌｄｂｅｆｏｒｅａｎｄａｆｔｅｒｔｈｅＴａｎｇｓｈａｎ...     

1786年康定地震形变特征的初步研究

本文对鲜水河断裂带南东段,康定断裂地震的形变带进行了分析。认为,地震形变带主要由发育于地形斜坡上的线性坡中槽或垄岗组成。形变带具分段特点,单条长850—1500米,呈右阶“斜列式”展布,中段(极震区)一带为现状型。其中的破裂面具正断兼扭动特点。空间特征上,坡中槽一侧的交替上升变化是依次、轮换出现的,它是地震断层运动屈曲作用(Fault buckling)导致地表变形的反映。这种形变现象与该带北西段(炉霍段)走滑型地震的形变带相比有明显的差异,也表明鲜水河断裂带北西、南东两段的地震破裂方式是不尽相同的,它为同一走滑带不同地段运动特征的差异提供了证据。     

弹性波的三维有限元模拟

本文讨论了弹性波瞬态传播问题的三维有限元计算方法及当前存在的实际困难.针对要求计算机内存大和计算时间长的问题,采取了改进措施.由于采用了集中质量矩阵和修正的中心差分时间积分显格式相结合的方法,可以使计算机内存和计算时间大为减少;由于采用结点定位法,最适合用于目前发展的并行计算机系统,可使计算速度大大增快;还采用了有效激发震源法,有效激发区是随时间步进的增加而逐步增大,这不仅能节省计算时间,而且使波场的传播过程一目了然,本文计算了由两种介质组成的三维楔形问题,得到若干典型剖面的瞬时波场图及随时间变化的合成地震图.