首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 703 毫秒
1.
A linear boundary element (BE) model is proposed for the uncoupied simulation of land subsidence due to gas, oil and hot water production over three-dimensional (3-D) arbitrarily shaped reservoirs. The pore pressure decline is assumed to be specified in advance, e.g. via a numerical model of flow. Use is made of the fundamental solution derived in 1885 by Boussinesq for a vertical load acting upon the traction-free surface of a semi-infinite medium. A straightforward application of Betti's (1872) reciprocal theorem allows for the development of a boundary integral whose numerical execution yields directly the downward settlement over the point of interest. The new procedure is applied to assess land sinking caused by an uniform pore pressure decline occurring within fields of elliptical shape and to explore the influence of the assumption of small reservoir thickness which underlies the ‘tension center’ or ‘strain nucleus’ approach previously developed by Geertsma in 1966. The results emphasize the numerical efficiency of the solution and the promising features of the BE method for the evaluation of ground subsidence in 3-D problems. The present model is based on the theory of the linear poroelasticity and is implemented for a mechanically homogeneous and isotropic half-space. It allows for any arbitrary geometry of the reservoir and for a non-uniform distribution of the pore pressure decline. It may easily be extended to other physical settings for which a vertical surface point load solution is available.  相似文献   

2.
A two-and-a-half-dimensional (2.5-D) coupled finite element–boundary element (FE–BE) model is presented to simulate the three-dimensional dynamic interaction between saturated soils and structures with longitudinally invariant geometries. A regularized 2.5-D boundary integral equation for saturated porous media is derived that avoids the evaluation of singular traction integrals. The 2.5-D coupled FE–BE model is established by using the continuity conditions on the soil–structure interface. The developed model is verified through comparison with an existing semi-analytical method. Two case studies of a tunnel embedded in a poroelastic half-space and the efficiency of a vibration isolating screen are presented.  相似文献   

3.
A century of land subsidence in Ravenna, Italy   总被引：4，自引：2，他引：4  
Over the last century the Ravenna area has been affected by a widespread land subsidence process of both natural and anthropogenic origin which has become a matter of increasingly great concern. Land settlement, initially of the order of a few mm/year, dramatically increased up to 110 mm/year after World War II primarily due to groundwater pumping and, subordinately, gas production from a number of deep on-shore and off-shore reservoirs. Spirit leveling was carried out in the area of the Ravenna Municipality by various authorities and agencies ever since the end of the nineteenth century using an increasingly refined network of benchmarks. The available information, spanning the time interval from 1897 to 2002, has been implemented and homogenized into a Geographic Information System (GIS). Measurements processed by kriging provide the pattern of the subsidence rates over the municipal territory during eight sub-periods. The reliability of data processing is discussed for both magnitude and areal distribution of the occurrence, showing a pronounced nonlinear behavior in time, as related to the major responsible factors. Cumulative land subsidence approaches almost 1 m over more than one third of the municipality area including the city and a significant fraction of the coastland, with peaks larger than 1.5 m over a 10 km2 zone located between the historical center and the coastline. Most recent records show that at present the mainland appears to be substantially stable, whereas some coastal areas are still subsiding at a rate of about 10 mm/year with a significant threat to the environment and the infrastructures as well, although to a lesser extent than in the past decades.  相似文献   

4.
The stress variation induced by gas/oil production may activate pre‐existing regional faults. This may enhance the expected land subsidence due to the generation of mechanically weak points close to the producing field. A class of elasto‐plastic interface elements (IE), specifically designed to address the mechanical behaviour of faults over a regional scale, is integrated into a finite element (FE) geomechanical model and used to investigate the role exerted by active faults in anthropogenic land subsidence. The importance of regional faults depends on a variety of factors including depth of the depleted reservoir, fault number, orientation and size, geomechanical properties of porous medium, pore pressure drawdown induced by fluid production, etc. With the aid of some representative examples, a useful indication is provided as to where and how fault activation may influence both magnitude and extent of the land subsidence bowl above producing gas/oil reservoirs, pointing to a generally limited impact on the ground surface. The simulation of a real faulted gas reservoir in a complex 3‐D setting shows that the proposed IE can be simply and efficiently incorporated into a FE geomechanical model, thus improving the quality of the stress and displacement prediction. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

5.
In this paper, a new approach to applying confining stress to flexible boundaries in the smoothed particle hydrodynamics (SPH) method is developed to facilitate its applications in geomechanics. Unlike the conventional SPH methods that impose confining boundary conditions by creating extra boundary particles, the proposed approach makes use of kernel truncation properties of SPH approximations that occur naturally at free-surface boundaries. Therefore, it does not require extra boundary particles and, as a consequence, can be utilised to apply confining stresses onto any boundary with arbitrary geometry without the need for tracking the curvature change during the computation. This enables more complicated problems that involve moving confining boundaries, such as confining triaxial tests, to be simulated in SPH without difficulties. To further enhance SPH applications in elasto-plastic computations of geomaterials, a robust numerical procedure to implement Mohr-Coulomb plasticity model in SPH is presented for the first time to avoid difficulties associated with corner singularities in Mohr-Coulomb model. The proposed approach was first validated against two-dimensional finite element (FE) solutions for confining biaxial compression tests to demonstrate its predictive capability at small deformation range when FE solutions are still valid. It is then further extended to three-dimensional conditions and utilised to simulate triaxial compression experiments. Simulation results predicted by SPH show good agreement with experiments, FE solutions, and other numerical results available in the literature. This suggests that the proposed approach of imposing confining stress boundaries is promising and can handle complex problems that involve moving confining boundary conditions.  相似文献   

6.
A sinking of the land surface due to the pumping of groundwater has long been recognized as an environmental issue in the Shiroishi plain of Saga, Japan. Land subsidence can have several negative economic and social implications such as changes in groundwater and surface water flow patterns, restrictions on pumping in land subsidence prone areas, localized flooding, failure of well casings as well as shearing of structures. To minimize such an environmental effect, groundwater management should be considered in this area. In this study, a new integrated numerical model that integrates a three-dimensional numerical groundwater flow model coupled with a one-dimensional soil consolidation model and a groundwater optimization model was developed to simulate groundwater movement, to predict ground settlement and to search for optimal safe yield of groundwater without violating physical, environmental and socio-economic constraints. It is found that groundwater levels in the aquifers greatly vary from season to season in response to the varying climatic and pumping conditions. Consequently, land subsidence has occurred rapidly throughout the area with the Shiroishi plain being the most prone. The predicted optimal safe yield of the pumping amount is about 5 million m3. The study also suggests that pumping with this optimal amount will minimize the rate of land subsidence over the entire area. An erratum to this article can be found at  相似文献   

7.
The solution of the poroelastic equations for predicting land subsidence above productive gas/oil fields may be addressed by the principle of virtual works using either the effective intergranular stress, with the pore pressure gradient regarded as a distributed body force, or the total stress incorporating the pore pressure. In the finite element (FE) method both approaches prove equivalent at the global assembled level. However, at the element level apparently the equivalence does not hold, and the strength source related to the pore pressure seems to generate different local forces on the element nodes. The two formulations are briefly reviewed and discussed for triangular and tetrahedral finite elements. They are shown to yield different results at the global level as well in a three‐dimensional axisymmetric porous medium if the FE integration is performed using the average element‐wise radius. A modification to both formulations is suggested which allows to correctly solve the problem of a finite reservoir with an infinite pressure gradient, i.e. with a pore pressure discontinuity on its boundary. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

8.
Coal seam gas (CSG) is an increasingly important source of natural gas all over the world. Although the influence of conventional oil and gas extraction on surface subsidence has been widely recognized and studied, few studies are carried out on the surface subsidence in coal seam gas fields and its impact on surface infrastructure and the environment. This paper discusses modelling of the surface subsidence associated with coal seam gas production by applying both analytical and numerical methods. By comparison of results from the numerical model and two analytical models, i.e. the disc-shaped reservoir model and the uniaxial compaction model, the analytical solutions cannot describe the complex process of water and gas extraction and have the limitations to predict the surface subsidence, while the numerical model can be better used in prediction of subsidence. After applying the numerical model in numerical analysis, the deformation characteristics of coupled fluid flow, and the effects of permeability change of coal seam, associated overlying and underlying layers, and depressurization rates on surface subsidence are investigated. The results demonstrate that the proposed model can simulate the production of water and gas from coal seams and the associated surface subsidence.  相似文献   

9.
A 3D elasto-plastic rate-dependent model for rock mechanics is formulated and implemented into a Finite Element (FE) numerical code. The model is based on the approach proposed by Vermeer and Neher (A soft soil model that accounts for creep. In: Proceedings of the International Symposium “Beyond 2000 in Computational Geotechnics,” pages 249-261, 1999). An original strain-driven algorithm with an Inexact Newton iterative scheme is used to compute the state variables for a given strain increment.The model is validated against laboratory measurements, checked on a simplified test case, and used to simulate land subsidence due to groundwater and hydrocarbon production. The numerical results prove computationally effective and robust, thus allowing for the use of the model on real complex geological settings.  相似文献   

10.
The proposed injection of seawater in a brackish, sandy aquifer lying 600–800 m under the lagoon of Venice to lift Venice itself uniformly up to 30 cm in 10 years raises many problems. We discuss here whether a smooth uplift can be obtained as claimed by the proponents. We first underline the similarities between the Ravenna and Venice case. We then examine the subsidence around Ravenna because of fluid extraction and the measured surface displacements for the period when there is pressure recovery both in the depleted gas reservoir close to Ravenna and in the upper aquifer system. Uniform surface settlements are nowhere observed and strong influence of local heterogeneities cannot be ruled out. The whole picture suggests great caution and the need for much more investigation before undertaking such an injection.  相似文献   

11.
基于部分耦合原理,采用TOUGH2和FLAC3D建立抽水引起的三维地面沉降弹塑性模型,模型中综合考虑土体的弹塑性变形特征、渗流-应力的双向耦合作用以及参数的非线性,探讨了持续抽水和脉冲抽水两种抽水过程中地面沉降发展演化过程。研究结果表明:（1）集中抽水停止后地面沉降会发生回弹,抽水中心沉降量不断减小。由于水平方向存在水力梯度,地下水继续向地下水位漏斗中心渗流从而导致沉降漏斗的范围仍继续扩大;（2）脉冲抽水导致土体的孔隙水压力、渗透系数以及沉降量均呈周期性波动变化,地面沉降会局部回弹,但总体仍随着抽水的持续,沉降量不断增加;（3）在抽水量相同前提下,对比持续抽水与脉冲抽水两种方式引发的塑形沉降量可知,抽水速率小、脉冲式多次开采导致的塑性沉降量较小,持续抽水的抽水速率越小、脉冲抽水间隔越短越有利于控制地面沉降。研究成果为地面沉降数值模拟提供了一种新方法,其中算例研究能为抽水条件下地面沉降的控制提供参考。  相似文献   

12.
The multi-disciplinary work described in the paper was aimed at analyzing and predicting the cyclical ground surface movements induced by underground gas storage (UGS) activities in a depleted gas field located in the Po Plain (Italy). The field has been operated as a storage facility for nearly three decades. Currently, the possibility of delta-pressuring the reservoir (i.e. to increase the maximum operating pressure above the initial reservoir pressure) to enhance the storage performance is being considered. Significant information was collected over time: 2/3D seismic surveys, geological and sedimentological studies, 60+ logged wells, geotechnical lab tests and 50+ years of production history and monitoring were available for the development of a fully integrated static–dynamic–geomechanical analysis. The mechanical aspects of the study are the focus of this paper. The data coming from different sources at different scales were analysed and integrated to set up and characterize a 3D finite element method mechanical model to calculate the surface movements induced by UGS activity by adopting an elasto-plastic constitutive law. The model was then calibrated via a back analysis approach, i.e. the model parameters were fine-tuned so that the simulated subsidence/uplift would compare satisfactorily with the ground movements collected over nearly 10 years of monitoring via interferometric synthetic aperture radar analysis in the region under investigation. Eventually, the calibrated model was used as a forecasting tool for subsidence evaluation under different future storage strategies, including delta-pressuring conditions. Results proved that no significant subsidence is expected even if the maximum operating pressure reached 120 % of the initial formation pressure.  相似文献   

13.
This paper presents a combined spheropolyhedral discrete element (DE)–finite element (FE) computational approach to simulating vertical plate loading on cohesionless soils such as gravels. The gravel particles are modeled as discrete elements, and the plate is modeled as a deformable FE continuum. The simulations provide a meaningful step toward better understanding how deformable bodies transmit loads to granular materials. The DE–FE contact algorithm is verified through comparison with an analytical solution for impact between two symmetric bars. A parametric study is conducted to ensure boundary effects are not significantly influencing the simulations. Numerical simulations are compared to experimental test results of lightweight deflectometer loading on a gravel base course with satisfactory agreement. Future developments of the approach intend to simulate wheel loading of military aircraft on unsurfaced airfields.  相似文献   

14.
王刚  杨鑫祥  张孝强  薛娇  李文鑫 《岩土力学》2015,36(11):3322-3328
为了建立能够表征煤体真实孔隙结构的数字模型,进而进行数值模拟研究,通过μCT225kVFCB高精度CT系统扫描得到了大柳塔煤矿长焰煤CT数据,可观测到的最小孔隙直径为1.94 μm。使用基于Matlab语言的三维重建程序结合逆向工程技术,提出了一种将煤体CT三维数据转换为CAD数字模型的方法。以Ansys建立的煤体有限元模型为例,对煤体瓦斯渗流进行了模拟,分析了煤体孔隙内的速度及压力分布规律,并计算了沿X、Y、Z方向的渗透系数。计算结果表明：在微观尺度下（< 100 μm）煤体渗透系数呈现各向异性,其受煤体结构的影响较明显。提出的基于CT三维重建结合逆向工程技术构建的煤体CAD数字模型不仅可进行有限元分析,同时也可以被EDEM等离散元分析软件所使用,拓宽了煤体CT三维数据的应用范围,丰富了煤体在微观尺度上的研究方法。  相似文献   

15.
Shanghai, in China, has experienced two periods of rapid land subsidence mainly caused by groundwater exploitation related to economic and population growth. The first period occurred during 1956–1965 and was characterized by an average land subsidence rate of 83 mm/yr, and the second period occurred during 1990–1998 with an average subsidence rate of 16 mm/yr. Owing to the establishment of monitoring networks for groundwater levels and land subsidence, a valuable dataset has been collected since the 1960s and used to develop regional land subsidence models applied to manage groundwater resources and mitigate land subsidence. The previous geomechanical modeling approaches to simulate land subsidence were based on one-dimensional (1D) vertical stress and deformation. In this study, a numerical model of land subsidence is developed to simulate explicitly coupled three-dimensional (3D) groundwater flow and 3D aquifer-system displacements in downtown Shanghai from 30 December 1979 to 30 December 1995. The model is calibrated using piezometric, geodetic-leveling, and borehole extensometer measurements made during the 16-year simulation period. The 3D model satisfactorily reproduces the measured piezometric and deformation observations. For the first time, the capability exists to provide some preliminary estimations on the horizontal displacement field associated with the well-known land subsidence in Shanghai and for which no measurements are available. The simulated horizontal displacements peak at 11 mm, i.e. less than 10 % of the simulated maximum land subsidence, and seems too small to seriously damage infrastructure such as the subways (metro lines) in the center area of Shanghai.  相似文献   

16.
双列块法用于有限元分析计算中求解大型线代数方程组。目前在有限元计算中最通常使用的求解大型线代数方程组的方法是分块的三角分解法。由于这种方法对方程组系数矩阵的分块受带宽的限制‘以至目前一般计算机的内存贮量限制了大带宽问题的求解,使一些大型的有限元课题,特别是三维问题的解算难以实现。双列块法解决了这个西难。双列块法对方程组系数矩阵实行按计算机许可的容量划分列块,按列分解,可以完全不受带宽的限制。方法是成功的、有效的。用双列块法在PE-3220小型机上仅用686K字节的内存解算了最大半带宽为1746,系数矩阵存贮量为5000K以上字节的三维(124个20节点等参元,8个16节点等参节理元、831个节点)有限元课题,这种方法使小型和高档微机用于大型有限元的解算成为可能。  相似文献   

17.
For a potential geological barrier of high-level radioactive waste repositories in France, the long-term mechanical behavior of the Callovo-Oxfordian (COx) argillite is the most concern for engineers. In this paper, a micromechanical-based elasto-viscoplastic model is proposed, and its numerical realization is our main object. The COx argillite is considered as a three-phase composite consists of porous clay, quartz, and calcite. By assigning appropriate constitutive laws to those constituents, the macroscopic elasto-viscoplastic behavior of the COx argillite is determined with an extended Hill's incremental approach. The numerical aspects includes (a) a new formulation is proposed for the plastic multiplier when adopting the overstress (Perzyna) model to define the viscoplastic strain. Meanwhile, a new formulation is also proposed to solve it within the framework of an implicit returning mapping scheme. (b) The corresponding consistent tangent moduli are strictly derived by extending the method proposed for solving plastic problems; (c) the efficiency of the proposed integration algorithms for the local constitutive equations and the homogenization procedure are validated, receptively, by a built-in porous plasticity model of a commercial finite element (FE) program ABAQUS and by FE computations of a two-phase unit cell; and (d) the proposed micromechanical model is finally applied to simulate experiment data in short-term triaxial compression tests and long-term triaxial creep tests. And the numerical results show that it is able to reflect the variation of the mechanical behavior with respect to the varied mineralogical compositions.  相似文献   

18.
周万欢  殷建华 《岩土力学》2011,32(Z1):691-0696
灌浆土钉在边坡加固工程中有着广泛的应用。土钉在不同条件下的抗拔机制仍是进一步研究的课题。通过建立一个三维有限元模型,研究土钉在不同上覆压力和剪胀作用下的土钉抗拔行为。采用摩尔-库仑模型模拟花岗岩土,土与土钉界面采用“接触对”和库仑摩擦模型描述。模型可以模拟土钉发生较大抗拔位移下发生滑移破坏的过程。模型参数通过三轴试验结果进行校正。通过与大尺寸室内抗拔试验结果相比较,验证了模型的有效性。模型展现了土钉在钻孔后和抗拔过程中土钉周围的土压力的变化。最后通过参数分析,研究了上覆压力和剪胀对土钉抗拔力的影响。结果显示,当上覆压力较小时,剪胀角对抗拔应力的影响较小。随着上覆压力的增加,剪胀角对抗拔应力的影响变得更加明显。  相似文献   

19.
The double Column block method applies to solute of large sparse linear simultaneous equations in finite element analysis. At present the portioned triangular decomposition method is generally used for solving large systems of linear equations in FE. With the partition of the coefficient matrix dependent on its bandwidth, the solution of problems with a large bandwidth is restricted owing to the computer core storage and hence it is difficult for this method to solve some large-scal FE problems, especially the three-dimensional problems. This difficulty has been overcome by the double column block method, completely independent of the bandwidth, by means of our method the coefficient matrix is partitioned according to computer core storage capacity. It is proved that this approach is successful and effective. Based on this method, the solution has been completed of the three-dimensional FE problem with 124 twenty-node and 8 sixteen-node isoparametric elements, and 831 nodes, having a maximum half-bandwidth of 1746 and global stiffness matrix storage of more than 5000 K bytes. The computation was performed on PE-3220 minicomputer and only a core storage of 685 K bytes was used. The double column block method makes it possible for minicomputer and high-level microcomputer to be applied to calculations in large-scale FE problems.  相似文献   

20.
本文将给出垂直断层附近三维不均匀体电阻率法的数值解，所用方法是边界单元法，但基本解的选择与以往的边界单元法不同。  相似文献   

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

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