共查询到20条相似文献，搜索用时 46 毫秒
1.
Using the extended finite element method, faults can be introduced into a three-dimensional reservoir-geomechanical model without meshing to assess the potential for fault reactivation associated with industrial activities such as disposal of CO2 or fluid extraction. The residuals to the governing equations include basis functions formed from the enrichment functions for strong and weak discontinuities. The traditional spatial integration scheme is based on an algorithm that partitions each enriched element into a collection of subtetrahedra, which is time and storage intensive. To avoid element partitioning, we adopt the homogeneous numerical integration scheme (referred as HNI hereafter) where the integration of homogeneous monomials over each polyhedron is converted into the integration of the same monomials over the one-dimensional edges of the polyhedron by using Stokes's theorem and Euler's homogeneous function theorem. The integrands for the strong and weak discontinuities are derived to implement the integration scheme in the three-dimensional reservoir-geomechanical model. An example of a jagged fault is presented to reveal the advantages of HNI scheme over the standard element partitioning approach. Several other examples that involve fluid flow, fault sliding, and fault sliding triggered by injection pressure are also presented to demonstrate that accurate and efficient computations are realized by the new integration scheme. 相似文献
2.
In this paper, a series of multimaterial benchmark problems in saturated and partially saturated two‐phase and three‐phase deforming porous media are addressed. To solve the process of fluid flow in partially saturated porous media, a fully coupled three‐phase formulation is developed on the basis of available experimental relations for updating saturation and permeabilities during the analysis. The well‐known element free Galerkin mesh‐free method is adopted. The partition of unity property of MLS shape functions allows for the field variables to be extrinsically enriched by appropriate functions that introduce existing discontinuities in the solution field. Enrichment of the main unknowns including solid displacement, water phase pressure, and gas phase pressure are accounted for, and a suitable enrichment strategy for different discontinuity types are discussed. In the case of weak discontinuity, the enrichment technique previously used by Krongauz and Belytschko [Int. J. Numer. Meth. Engng., 1998; 41:1215–1233] is selected. As these functions possess discontinuity in their first derivatives, they can be used for modeling material interfaces, generating only minor oscillations in derivative fields (strain and pressure gradients for multiphase porous media), as opposed to unenriched and constrained mesh‐free methods. Different problems of multimaterial poro‐elasticity including fully saturated, partially saturated one, and two‐phase flows under the assumption of fully coupled extended formulation of Biot are examined. As a further development, problems involved with both material interface and impermeable discontinuities, where no fluid exchange is permitted across the discontinuity, are considered and numerically discussed. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
3.
Reconstruction of geological structures has the potential to provide additional insight into the effect of the depositional history on the current-day geomechanical and hydro-geologic state. Accurate modeling of the reconstruction process is, however, complex, necessitating advanced procedures for the prediction of fault formation and evolution within fully coupled geomechanical, fluid flow and temperature fields. In this paper, a 3-D computational approach is presented that is able to forward model complex structural evolution with multiple intersecting faults that exhibit large relative movement within a coupled geomechanical/flow environment. The approach adopts the Lagrangian method, complemented by robust and efficient automated adaptive meshing techniques, an elasto-plastic constitutive model based on critical state concepts, and global energy dissipation regularized by inclusion of fracture energy in the equations governing state variable evolution. The proposed model is validated by comparison of 2-D plane strain and 3-D thin-slice predictions of a bench-scale experiment, and then applied to two conceptual coupled geomechanical/fluid flow field-scale benchmarks. 相似文献
4.
汶川8级地震地质灾害的类型及实例 总被引：11，自引：2，他引：9
龙门山地区活动断裂右旋斜冲运动与汶川8级强烈地震存在成因联系。汶川8级地震造成了惨重的人员伤亡和巨大的财产损失,地震地质灾害主要类型有地震灾害、地震触发地质灾害和地震引发地质灾害隐患。严重地震灾害包括房屋倒塌与部分坍塌、房屋平移、房屋倾斜变形、墙体破裂与结构破坏、桥梁垮塌等。地震触发严重地质灾害包括山体滑坡、山体滑塌、岩块崩塌,局部产生泥石流与沙土液化。地震引发地质灾害隐患包括潜在滑坡、不稳定边坡与滑塌隐患、危岩体与崩塌隐患及泥石流灾害隐患。地震地质灾害分布与活动断层存在密切关系,沿北川映秀断裂地震地质灾害最为严重,沿汉旺漩口断裂、茂县汶川断裂、青川断裂地震地质灾害也比较严重。在活动断裂之间相对稳定地块远离活动断裂超过3~5 km,地震灾害和地震触发地质灾害显著减轻。 相似文献
5.
This paper presents a fracture mapping (FM) approach combined with the extended finite element method (XFEM) to simulate coupled deformation and fluid flow in fractured porous media. Specifically, the method accurately represents the impact of discrete fractures on flow and deformation, although the individual fractures are not part of the finite element mesh. A key feature of FM‐XFEM is its ability to model discontinuities in the domain independently of the computational mesh. The proposed FM approach is a continuum‐based approach that is used to model the flow interaction between the porous matrix and existing fractures via a transfer function. Fracture geometry is defined using the level set method. Therefore, in contrast to the discrete fracture flow model, the fracture representation is not meshed along with the computational domain. Consequently, the method is able to determine the influence of fractures on fluid flow within a fractured domain without the complexity of meshing the fractures within the domain. The XFEM component of the scheme addresses the discontinuous displacement field within elements that are intersected by existing fractures. In XFEM, enrichment functions are added to the standard finite element approximation to adequately resolve discontinuous fields within the simulation domain. Numerical tests illustrate the ability of the method to adequately describe the displacement and fluid pressure fields within a fractured domain at significantly less computational expense than explicitly resolving the fracture within the finite element mesh. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
6.
Impact of single inclined faults on the fluid flow and heat transport: results from 3-D finite element simulations 总被引：1，自引：1，他引：0
Yvonne Cherubini Mauro Cacace Guido Blöcher Magdalena Scheck-Wenderoth 《Environmental Earth Sciences》2013,70(8):3603-3618
The impact of inclined faults on the hydrothermal field is assessed by adding simplified structural settings to synthetic models. This study is innovative in carrying out numerical simulations because it integrates the real 3-D nature of flow influenced by a fault in a porous medium, thereby providing a useful tool for complex geothermal modelling. The 3-D simulations for the coupled fluid flow and heat transport processes are based on the finite element method. In the model, one geological layer is dissected by a dipping fault. Sensitivity analyses are conducted to quantify the effects of the fault’s transmissivity on the fluid flow and thermal field. Different fault models are compared with a model where no fault is present to evaluate the effect of varying fault transmissivity. The results show that faults have a significant impact on the hydrothermal field. Varying either the fault zone width or the fault permeability will result in relevant differences in the pressure, velocity and temperature field. A linear relationship between fault zone width and fluid velocity is found, indicating that velocities increase with decreasing widths. The faults act as preferential pathways for advective heat transport in case of highly transmissive faults, whereas almost no fluid may be transported through poorly transmissive faults. 相似文献
7.
异源点集下三维断层建模方法 总被引：2，自引：0，他引：2
用于断层建模的异源点集,因其具有不同特征而影响到建模途径和方法。以三维断层模型的3D GIS应用为目标,分析了原始地质观测点集、钻孔点集、剖面点集和平面点集的信息构成特征,提出了异源点集下三维断层数据模型的构建方法和流程。原始地质观测点集和钻孔点集间接表达了断层的空间特征,通过空间平面、曲面和曲线构建了断层面、地层面和断层交线的数学模型,并借助于连续模型的离散化,实现了独立断层建模点集的重建;剖面点集和平面点集可直接用于断层建模。以双TIN面模型为例阐述了断层三维模型的构建方法。 相似文献
8.
断裂剪切带厘米级磨砾和纳米级磨粒的发现、命名及其油气地质意义 总被引：2，自引：0，他引：2
上世纪70年代初,在江苏江宁逆断层中发现了厘米（cm）级磨砾,这种磨粒通常在韧脆性断裂剪切带中都可能找到。80年代末,对美国加州深钻岩芯进行扫描电子显微镜观测时,在晚中生代花岗质糜棱岩滑移叶片中发现了纳米（nm）级超微球粒结构,球粒（2—5nm）紧密聚集、定向明显。因其形成同滚动、研磨和剪切作用过程有关,故将球形微粒称“磨粒”。嗣后,在陕北和黄海深钻岩芯石炭-二叠系泥岩滑移叶片中,也发现类同现象（微米-纳米级）。近年来又在郯庐断裂、日本中央构造线大型断裂带,以及南京幕府山寒武系至三叠系不同地层的各种剪切滑移构造中都发现了nm级磨粒构造。泥岩中含磨粒滑移层具有涂抹效应,起到油气封盖作用,这在陕北油气田已得到了验证。对其磨粒构造蕴藏的丰富地质信息,特别是在微观动力学方面（包括动热变质作用、流体渗散作用、颗粒自旋作用和地质催化作用等）应作进一步研究。 相似文献
9.
青藏高原东缘活动断裂地质灾害效应研究 总被引：6，自引：0，他引：6
活动断裂的地质灾害效应是工程地质与地质灾害领域研究的重要内容。本文以第四纪以来构造活动最强烈的青藏高原东缘为例,阐述了活动断裂地质灾害效应的主要表现形式,包括:(1)活动断裂对地形地貌和岩体结构的影响;(2)断裂剧烈活动(地震)诱发地质灾害;(3)断裂蠕滑作用对斜坡应力场和稳定性的影响;(4)断裂活动是地质灾害链的源头,为地质灾害提供物源。上述表现形式及灾害成因机理和分布规律是活动构造区地质灾害防治中需要关注的关键问题。根据青藏高原东缘典型地质灾害案例研究提出,内外动力耦合作用成灾机理是未来地质灾害研究方向,将为活动构造区地质灾害早期识别和防灾减灾提供理论依据。 相似文献
10.
Understanding the way fluids flow in fault zones is of prime importance to develop correct models of earthquake mechanics, especially in the case of the abnormally weak San Andreas Fault (SAF) system. Because fluid flow can leave detectable signatures in rocks, geochemistry is essential to bring light on this topic. The present detailed study combines, for the first time, C–O isotope analyses with a comprehensive trace element data set to examine the geometry of fluid flow within a significant fault system hosted by a carbonate sequence, from a single locality across the Little Pine Fault–SAF system. Such a fault zone contains veins, deformation zones, and their host rocks. Stable isotope geochemistry is used to establish a relative scale of integrated fluid–rock ratios. Carbonate δ18O varies between 28‰ and 15‰ and δ13C between 5‰ and −7‰. From highest to lowest delta values, thus from least to most infiltrated, are the host rocks, the vein fillings, and the deformation zone fillings, respectively. Infiltration increases toward fault core. The fluids are H2O–CO2 mixtures. Two fluid sources, one internal and the other external, are found. The external fluid is inferred to come essentially from metamorphism of the Franciscan formation underneath. The internal (local) fluid is provided by a 30% volume reduction of the host limestones resulting from pressure solution and pore size reduction. Most trace elements, including the lanthanides, show enrichment at the 100-m scale in host carbonate rocks as fluid–rock ratios increase toward the fault core. In contrast, the same trace element concentrations are low, relative to host rocks, in veins and deformation zone carbonate fillings, and this difference in concentration increases as fluid–rock ratio increases toward the fault core. We suggest that the fluid trace elements are scavenged by complexation with organic matter in the host rocks. Elemental complexation is especially illustrated by large fractionation of Y–Ho and Nb–Ta geochemical pairs. Complexation associated with external fluid flow has a significant effect on trace element enrichment (up to 700% relative enrichment) while concentration by pressure solution associated with volume decrease of host rocks has a more limited effect (up to 40% relative enrichment). Our observations from the millimeter to the kilometer scale call for the partitioning of fluid sources and pathways, and for a mixed focused–pervasive fluid flow mechanism. The fluid mainly flows within veins and deformation zones and, simultaneously, within at least 10 cm from these channels, part of the fluid flows pervasively in the host rock, which controls the fluid composition. Scavenging of the fluid rare earth elements (REE) by host rocks is responsible for the formation of REE-depleted vein and deformation zone carbonate fillings. Fluid flow is not only restricted to veins or deformation zones as commonly believed. An important part of fluid flow takes place in host rocks near fault zones. Hence, the nature of the lithologies hosting fault zones must be considered in order to take into account the role of fluids in the seismic cycle. 相似文献
11.
Dislocation modelling of an earthquake fault is of great importance due to the fact that ground surface response may be predicted by the model. However, geological features of a fault cannot be measured exactly, and therefore these features and data involve uncertainties. This paper presents a Monte Carlo based random model of faults with finite element method incorporating split node technique to impose the effects of discontinuities. Length and orientation of the fault are selected as random parameters in the domain model, and hence geometrical uncertainties are encountered. Mean and standard deviation values, as well as probability density function of ground surface responses due to the dislocation are computed. Based on analytical and numerical calculation of dislocation, two approaches of Monte Carlo simulations are proposed. Various comparisons are examined to illustrate the capability of both methods for random simulation of faults. 相似文献
12.
以北京地区地热流体中温度和矿物质含量较高的东坝凹陷南部地区作为实例,在完整的地热钻孔及相关地质资料基础上,通过对研究区断裂构造、热储层、热储盖层、地温场、流体水化学特性和富水性等地热地质条件进行分析,研究造成该区地热流体中温度、矿物质含量升高的主要原因。研究表明：受太阳宫断裂北段、良乡—前门断裂和楼梓庄断裂的切割作用,东坝凹陷南部地区热储层下落近千米,形成了一个较为封闭的“黑箱子”区域。这种独特的地质结构导致东坝凹陷南部地区接受东南城区地热田侧向补给较少,长期处于较为封闭的还原环境;太阳宫断裂北段存在深部的热源通道,地层深部的高温热流顺断裂上涌在东坝凹陷南部地区的热储层富集,因此该区地热流体表现出温度高、矿物质含量高、富水性差的特征,该研究对我国今后地热资源的开发具有积极的借鉴意义。 相似文献
13.
随着油气勘探程度的加深,岩性油气藏的重要性日益显现,而查明油气富集规律是岩性油气藏勘探的关键。渤海湾盆地东营凹陷牛庄洼陷沙三中亚段发育大量岩性油气藏,文章以此为研究目标,着重分析其平面上油气分布的主要控制因素。研究结果表明,牛庄洼陷岩性油气藏的分布受控于洼陷内流体性质、断层作用和沉积相展布3种地质因素。它们对洼陷区岩性油气藏富集规律的影响分别在于：洼陷内流体运聚单元控制着油气在平面上的富集区,牛庄洼陷南北两侧汇聚流发育区是油气聚集的有利地带,而洼陷北部平行流和东西两侧发散流发育区不利于油气成藏;同生断层泵吸作用使得油气在牛庄洼陷平行流和发散流发育区可形成局部汇聚流,从而聚集成藏,断层的侧向遮挡作用则使得牛庄洼陷生成的油气多富集在靠近边界断层的洼陷内侧;物性较好的浊积水道微相最有利于油气的富集,朵叶体和水道间次之,外缘最差。 相似文献
14.
针对长江中游洞庭湖防洪系统规模庞大、水流复杂、资料短缺和预报时限紧迫的实际条件,提出了具有层次和模块结构特点、一维与二维水流模拟、水力学与水文学方法、理论模型与补充信息相结合的建模途径.所建模型的湖泊部分采用无结构网格二维非恒定流高性能有限体积格式,以适应湖区复杂的边界形状和保持水量平衡;河网部分采用一维非恒定流守恒型显格式,避免隐格式矩阵算法的复杂性,同时有利于与二维模型的耦合及与各种复杂连通关系的显式连接.这种一、二维混合非恒定流模型可用于长江干流、洞庭湖河网及湖泊、堤垸区的洪水演进和防洪调度的水流仿真. 相似文献
15.
16.
从地震信息看阿尔金断裂带构造和塔里木盆地花彩弧断裂体系 总被引：1，自引：0，他引：1
阿尔金断裂带由多条断裂组成,主要有阿尔金断裂、且末断裂、三危山断裂。其中阿尔金断裂为主断裂,它呈左旋走滑兼具逆冲性质,中生代—古近纪为左旋走滑,新近纪由东南向西北逆冲推覆。且末断裂和三危山断裂均具左旋走滑性质。且末断裂受统一的阿尔金断裂带左旋应力场控制,但又叠加了塔里木台盆区向南挤压的应力场,从而具有双重属性。塔里木盆地的断裂总体上组成古生界塔北花彩弧断裂束和塔南花彩弧断裂束,展布成全盆地的菱形断裂系统,且末断裂构成其东南边界。在该菱形断裂系统的北弧顶和菱形内的中央轴部为背冲式的构造断裂带,显示挤压特征;在花彩弧两翼转弯处展布正花状构造样式,显示走滑特征。阿尔金断裂带及其两侧,主要在柴达木、塔里木两大盆地发现了大油气田,两者都是由断层控制油气的垂向运移与分布。柴达木盆地具有双重断—坳的特点,但油气田只分布在中—新生界构造层内;塔里木盆地,南北翘板式的构造运动是其形成复式油气区的最重要的地质构造条件。 相似文献
17.
煤炭地质勘查的研究对象、勘探方法、数据资料在空间上都是三维的,适合于采用三维地质建模技术进行管理。根据三维地质建模与煤炭地质勘查是否同时,可将煤炭地质勘探三维地质建模分为两类,一类为勘探工作中,主要以钻孔数据为主,用于制作三维地质模型的原型,生成煤层底板等高线、剖面图草图,供地质工程师进一步加工,构建精细的三维地质模型;另一类为勘探工作后,即以已经做好的煤层底板等高线图、剖面图为主要数据源,构建精细的三维地质模型。精细的三维地质模型可用于煤矿设计、数字矿山建设。三维地质建模技术的应用有利于煤炭地质构造的研究,主要表现在:第一,可对地面、钻探、物探、化探多种手段获得的探测资料进行集成与同化处理,提高信息空间的一致性;第二,地质构造要素表达具有直观性与三维立体性,如断层的表达,其在三维环境中不仅可表达为一个面,并且能清晰反映断层与地层、以及断层与断层之间的相互切割关系,弥补了平、剖面图上仅表现为一条断层线的不足;第三,可平、剖面三维联动编辑,很好地保证平面图与剖面图的空间一致性;第四,可以在剖面上以画草图的方式确定剖面上煤层对比与断层切割方案。 相似文献
18.
Modelling of coupled fluid‐mechanical problems in fractured geological media using enriched finite elements 
下载免费PDF全文
下载免费PDF全文 Jose Roberto Silvestre Euripedes do Amaral Vargas Jr. Luiz Eloy Vaz Antonio Claudio Soares 《国际地质力学数值与分析法杂志》2015,39(10):1104-1140
Geological environments, such as petroleum reservoirs, normally exhibit physical discontinuities, for example, fractures and faults. Because of the reduced thickness of these discontinuities, finite element formulations with strong discontinuity have been applied to the numerical modelling of geological environments. Until now, two relevant characteristics of petroleum reservoirs have not been addressed by these formulations. The first is the pore pressure jump in the direction normal to a discontinuity in a fluid‐mechanical coupling condition, which is present primarily in sealing faults owing to the contrast of permeability with the porous medium. The absence of this jump can affect the prediction of the deformability of a physical discontinuity. Furthermore, reservoir models frequently use coarse meshes. Thus, the method used to evaluate the pore pressure in the discontinuity may exhibit a strong dependence relative to the mesh refinement. Based on these characteristics, in this study, a formulation of an enriched finite element for application to coupled fluid‐mechanical problems with pre‐existing physical discontinuities saturated by a single fluid is presented. The formulation employs discontinuous interpolation functions and enables the reproduction of jumps of displacement and pore pressure associated with a discontinuity inside the element without the need to discretise it. An approximation to estimate the pore pressure in the discontinuity was developed, one which seeks to minimise the influence of refinement. The element's response is verified by comparison with a one‐dimensional analytical solution and simple examples that are simulated using commercial software. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
19.
20.
应用二维经验模分解（BEMD）法从滇东Pt-Pd含量数据中提取2类Pt-Pd异常,即由二维内蕴模函数3（BIMF3）表征的区域Pt-Pd异常和由二维内蕴模函数2（BIMF2）表征的局部Pt-Pd异常。根据异常产出的地质特征,区域Pt-Pd异常可进一步划分为2个亚类：其一是位于受SN向小江断裂和NE向寻甸-宣威深大断裂控制的二叠纪峨眉山玄武岩分布区的异常,其特征是具有明显的浓集中心,异常强度大;其二是分布于研究区中南部的黑色岩系中的异常,异常强度弱。其中前者是滇东Pt-Pd找矿的有利地段。局部Pt-Pd异常可进一步划分为3个亚类：①分布于玄武岩区,受SN-NE向深大断裂与其次级断裂交会域控制的异常;②分布于不同地质时期黑色岩系中的异常;③哀牢山断裂北侧与基性、超基性侵入岩有关的异常。其中,分布于玄武岩区的Pt-Pd局部异常强度高、规模大,是进一步找寻Pt-Pd资源的远景地段。 相似文献