共查询到20条相似文献，搜索用时 15 毫秒
1.
One of the major causes of instability in geotechnical structures such as dikes or earth dams is internal erosion, an insidious process that occurs over a long period of time. Research on this topic is still fairly new and much more needs to be understood in order to solve the problems posed by this phenomenon. This paper proposes a hydromechanical model based on porous continuous medium theory to assess how internal erosion impacts the safety of earthen structures. The saturated soil is considered as a mixture of four interacting constituents: soil skeleton, erodible fines, fluidized fine particles, and fluid. The detachment and transport of the fine particles are described by a mass exchange model between the solid and the fluid phases. An elastoplastic constitutive model for sand-silt mixtures has been developed to monitor the effect of the evolution of both porosity and fines content induced by internal erosion upon the behavior of the soil skeleton. The model has been numerically solved with the finite element method. It has then been applied to the specific case study of a dike foundation subjected to internal erosion induced by the presence of a karstic cavity beneath the alluvium layer. The numerical results show the onset of erosion, the time-space evolution of the eroded zone, and the hydromechanical response of the soil constituting the dike, all of which highlights the effects of the cavity location, the erosion rate, and the fines content. 相似文献
2.
Internal instability is a phenomenon of fine particle redistribution in granular materials under the seepage action and consequent change in the soil’s internal structure and hydraulic and mechanical properties. It is one of the primary causes of failures of sand-gravel foundations and embankment dams. The criteria establishment is considered the key to solving the erosion problems, so the existing internal stability criteria need a review and classification to study the recent development trends in soil seepage and erosion. Therefore, this paper aims at reviewing the internal stability factors of gap-graded soil with a focus on the internal erosion mechanism and internal stability evaluation based on geometric and hydraulic criteria. Firstly, the paper compared the effect of several commonly used geometric criteria for gap-graded soil evaluation, such as particle size, fine content, void ratio, and fractal dimension. Furthermore, it provided a hydraulic criteria overview and analyzed the effects of the hydraulic gradient, hydraulic shear stress, confining pressure, and pore velocity on internal erosion. The geometric–hydraulic coupling methods were introduced, with a detailed elaboration of the erosion resistance index method based on accumulated dissipated energy. The capabilities and limitations of these criteria were discussed throughout the paper. It was found that combined Kezdi’s criterion and Kenney and Lau’s criterion is more reliable to evaluate internal stability of soil. The gap-graded soil with fine particle content higher than 35% is not necessarily internally stable. Finally, the energy-based method (erosion resistance index method) can effectively reproduce the total amount of erosion mass and the final spatial distribution of fine particles and identifies erosion. The review's outcome can be used as a basis to evaluate the internal erosion risk for gap-graded soils. The evaluation methods discussed here can help identify the zones of relatively high erosion potential.
相似文献3.
含油气盆地普遍发育一个或多个剥蚀面,揭示剥蚀面的形成和油气藏的关系,不仅可以完善和发展已有的油气成藏理论,也必将对我国东部老油气区的勘探起到指导和促进作用。在对含油气盆地剥蚀次数和剥蚀量统计分析的基础上,建立了地层抬升剥蚀与油气成藏的空间和时间匹配关系;通过地质分析,从成藏机理的角度初步探讨了地层抬升剥蚀对油气成藏的贡献。统计分析和地质分析的结果显示:剥蚀层位本身或紧邻其下的层位就是储层;地层抬升剥蚀期多为油气聚集的时期;地层抬升过程中的流体温度降低及孔隙扩容可造成地层低压,从而为油气聚集提供动力,孔隙扩容还可为油气聚集提供更多的空间。 相似文献
4.
Currently, numerical studies at the real scale of an entire engineering structure considering internal erosion are still rare. This paper presents a three-dimensional (3D) numerical simulation of the effects of internal erosion within a linear dike located on a foundation. A two-dimensional (2D) finite element code has been extended to 3D in order to analyze the impact of internal erosion under more realistic hydromechanical conditions. The saturated soil has been considered as a mixture of four interacting constituents: soil skeleton, erodible fines, fluidized fine particles, and fluid. The detachment and transport of the fine particles have been modeled with a mass exchange model between the solid and the fluid phases. An elastoplastic constitutive model for sand-silt mixtures has been developed to monitor the effect of the evolution of both the porosity and the fines content induced by internal erosion upon the behavior of the soil skeleton. An unsaturated flow condition has been implemented into this coupled hydromechanical model to describe more accurately the seepage within the dike and the foundation. A stabilized finite element method was used to eliminate spurious numerical oscillations in solving the convection-dominated transport of fluidized particles. This numerical tool was then applied to a specific dike-on-foundation case subjected to internal erosion induced by a leakage located at the bottom of the foundation. Different failure modes were observed and analyzed for different boundary conditions, including the significant influence of the leakage cavity size and the elevation of the water level at the upstream and downstream sides of the dike. 相似文献
5.
This paper presents a computational method able to effectively model both the simultaneous processes typically observed in backward erosion piping, ie, the pipe tip propagation and the conduit cross-section enlargement. The numerical method is based on the novel formulation of a problem of localized erosion along a line propagating in a multidimensional porous medium. In this line, a conduit with evolving transverse size is embedded, which conveys a multiphase flow. The two systems, porous medium and pipe, are bridged by exchange terms of multiphase fluid mass and by a shared fluid pressure field. On the contrary, different fields are considered to describe flows, which are assumed as Darcian in the porous medium and turbulent in the conduit. These two flows drive pipe propagation and enlargement, respectively, as modeled by means of proper erosion kinetic laws. The corresponding numerical formulation is based on the combination between one- and multidimensional finite elements, to model the erosion conduit and the porous medium, respectively. Several simulations are proposed to demonstrate the ability of the proposed approach in reproducing available experimental data of real-scale tests on levees. Our results point out the crucial role played by the combined influence of pipe propagation and enlargement, as well as of three-dimensional (3D) effects. We also assess the mesh independence of the proposed numerical solution, particularly as concerns the calculated pipe propagation history. 相似文献
6.
Landslides can block mountainous streams and form landslide dams to threaten downstream residents. It is necessary for reliable methods to predict landslide dam dynamic for risk assessment. In this paper, we present a two-layer model of Savage–Hutter type to simulate the dynamic evolution of landslide dam which take account of the erosion of river bed. The two-layer shallow water system is derived by depth-averaging the incompressible Navier–Stokes equations with the hydrostatic assumption integrated of the erosion model of river bed. The effect of excess pore water pressure is considered in the erosion process. A high order accuracy scheme based on Roe-type solver is used to discretize the present model. Finally, several numerical tests are performed to verify the stability of the algorithm and reliability of the model. Numerical results indicate that the erosion effect enhances the huge destructiveness of landslide and increase the possibility of river blocked by landslides. The impact of excess pore water pressure on erosion process should be considered. 相似文献
7.
A general thermo-hydro-mechanical framework for the modelling of internal erosion is proposed based on the theory of mixtures applied to two-phase porous media. The erodible soil is partitioned in two phases: one solid phase and one fluid phase. The solid phase is composed of nonerodible grains and erodible particles. The fluid phase is composed of water and fluidized particles. Within the fluid phase, species diffuse. Across phases, species transfer. The modelling of internal erosion is contributed directly by mass transfer from the solid phase towards the fluid phase. The constitutive relations governing the thermomechanical behaviour, generalised diffusion, and transfer are structured by the dissipation inequality. The particular case of soil suffusion is investigated with a focus on constitutive laws. A new constitutive law for suffusion is constructed based on thermodynamic conditions and experimental investigations. This erosion law is linearly related to the power of seepage flow and to the erosion resistance index. Owing to its simplicity, this law tackles the overall trend of the suffusion process and permits the formulation of an analytical solution. This new model is then applied to simulate laboratory experiments, by both analytical and numerical methods. The comparison shows that the newly developed model, which is theoretically consistent, can reproduce correctly the overall trend of the cumulated eroded mass when the permeability evolution is small. In addition, the results are provided for four different materials, two different specimen sizes, and various hydraulic loading paths to demonstrate the applicability of the new proposed law. 相似文献
8.
Finite element simulations of two centrifuge tests on the same cantilever retaining wall model holding liquefiable backfill were conducted using the Biot formulation‐based program DIANA–SWANDYNE II. To demonstrate the effects due to different pore fluids in seismic centrifuge experiments, water was used as the pore fluid in one experiment whereas a substitute pore fluid was used in the second experiment. The cantilever wall model parameters were determined by comparing simulations with measurements from free‐vibration tests performed on the model wall without backfill. The initial stress conditions for dynamic analysis for the soil backfill were obtained by simulating static loads on the retaining wall from the soil backfill. Level‐ground centrifuge model results were used to select the parameters of the Pastor–Zienkiewicz mark III constitutive model used in the dynamic simulations of the soil. The effects due to different pore fluids were captured well by the simulations. The magnitudes of excess pore pressures in the soil, lateral thrust and its line of action on the wall, and wall bending strains, deflections, and accelerations were predicted well. Predictions of settlements and accelerations in the backfill were less satisfactory. Relatively high levels of Rayleigh damping were needed to be used in the retaining wall simulations in order to obtain numerically stable results, which is one of the shortcomings of the model. The procedure may be used for engineering purpose dealing with seismic analysis of flexible retaining walls where lateral pressures, bending strains and deflections in the wall are typically of importance. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
9.
A fully coupled transient two‐dimensional model was employed to study fundamentals of flood‐induced surface erosion in a particle bed. The interaction of the liquid and solid phases is the key mechanism related to surface erosion. The solid phase was idealized at a particle scale by using the discrete element method. The fluid phase was modeled at a mesoscale level and solved using the lattice Boltzmann method. The fluid forces applied on the particles were calculated on the basis of the momentum the fluid exchanges with the particle. The proposed approach was used to model both single particles and particle beds subjected to Couette flow conditions. The behavior of both the single particle and the particle bed depended on particle diameter and surface shear fluid velocity. The conducted simulations show that the fluid flow profile penetrates the bed for a small distance. This penetration initiates sheet‐flow and surface erosion as the fluid interacts with particles. The effect of suppressing particle rotation on the fluid‐induced forces on the particle was also examined. Suppressing particle spinning may lead to underestimated erosion rate. Results of fluid and particle velocities were compared against experimental results and appeared to agree with the observed trends.Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
10.
A novel two-dimensional mixed fracture–pore seepage model for fluid flow in fractured porous media is presented based on the computational framework of finite-discrete element method (FDEM). The model consists of a porous seepage model in triangular elements bonded by unbroken joint elements, as well as a fracture seepage model in broken joint elements. The principle for determining the fluid exchange coefficient of the unbroken joint element is provided to ensure numerical accuracy and efficiency. The mixed fracture–pore seepage model provides a simple but effective tool for solving fluid flow in fractured porous media. In this paper, examples of 1D and 2D seepage flow in porous media and porous media with a single fracture or multiple fractures are studied. The simulation results of the model match well with theoretical solutions or results obtained by commercial software, which verifies the correctness of the mixed fracture–pore seepage model. Furthermore, combining FDEM mechanical calculation and the mixed fracture–pore seepage model, a coupled hydromechanical model is built to simulate fluid-driven dynamic propagation of cracks in the porous media, as well as its influence on pore seepage and fracture seepage.
相似文献11.
This paper presents numerical simulations of Cone Penetration Test (CPT) in water-saturated soft soils taking into account pore pressure dissipation during installation. Besides modelling interaction between soil skeleton and pore fluid, the problem involves large soil deformations in the vicinity of the penetrometer, soil–structure interaction, and complex non-linear response of soil. This makes such simulations challenging. Depending on the soil’s permeability and compressibility, undrained, partially drained or drained conditions might occur. Partially drained conditions are commonly encountered in soils such as silts and sand–clay mixtures. However, this is often neglected in CPT interpretation, which may lead to inaccurate estimates of soil properties. This paper aims at improving the understanding of the penetration process in different drainage conditions through advanced numerical analyses. A two-phase Material Point Method is applied to simulate large soil deformations and generation and dissipation of excess pore pressures during penetration. The constitutive behaviour of soil is modelled with the Modified Cam Clay model. Numerical results are compared with experimental data showing good agreement. 相似文献
12.
烃源岩生排烃组分法模型研究与应用 总被引：4，自引：0，他引：4
本文提出了一种新的烃源岩生排烃组分法模拟模型。组分法生烃模型根据热模拟实验得到的组分产率曲线模拟烃源的生烃史。排烃模型包括开放体系模型和间歇开放体系模型，分别用于模拟以压实作用为主要动力的孔隙排烃和以异常高压为主要动力的微裂缝排烃。这一模型根据孔隙温度、压力和烃类组成，采用热力学方法模拟孔隙体系流体体积和相态的变化、不同组分在各相中的分布，并由此求出各组分的排出量和总排烃量。对陕甘宁盆地古生界源岩生排烃史的模拟结果表明，该模型具有实用价值。 相似文献
13.
S. Y. Wang S. W. Sloan S. G. Fityus D. V. Griffiths C. A. Tang 《Rock Mechanics and Rock Engineering》2013,46(5):1165-1182
Rock is a heterogeneous geological material. When rock is subjected to internal hydraulic pressure and external mechanical loading, the fluid flow properties will be altered by closing, opening, or other interaction of pre-existing weaknesses or by induced new fractures. Meanwhile, the pore pressure can influence the fracture behavior on both a local and global scale. A finite element model that can consider the coupled effects of seepage, damage and stress field in heterogeneous rock is described. First, two series of numerical tests in relatively homogeneous and heterogeneous rocks were performed to investigate the influence of pore pressure magnitude and gradient on initiation and propagation of tensile fractures. Second, to examine the initiation of hydraulic fractures and their subsequent propagation, a series of numerical simulations of the behavior of two injection holes inside a saturated rock mass are carried out. The rock is subjected to different initial in situ stress ratios and to an internal injection (pore) pressure at the two injection holes. Numerically, simulated results indicate that tensile fracture is strongly influenced by both pore pressure magnitude and pore pressure gradient. In addition, the heterogeneity of rock, the initial in situ stress ratio (K), the distance between two injection holes, and the difference of the pore pressure in the two injection holes all play important roles in the initiation and propagation of hydraulic fractures. At relatively close spacing and when the two principal stresses are of similar magnitude, the proximity of adjacent injection holes can cause fracturing to occur in a direction perpendicular to the maximum principal stress. 相似文献
14.
The purpose of this paper is to simulate the coupled dynamic deformation and water flow that occur in saturated soils when subjected to traffic loads, which is a problem with several practical applications. The wave propagation causes vibrations leading to discomfort for passengers and people in the surroundings and increase wear on both the vehicle and road structure. The water flow may cause internal erosion and material transport in the soil. Further, the increased pore water pressure could reduce the bearing capacity of embankments. The saturated soil is modelled as a water‐saturated porous medium. The traffic is modelled as a number of moving wheel contact loads. Dynamic effects are accounted for, which lead to a coupled problem with solid displacements, water velocity and pressure as primary unknowns. A finite element program has been developed to perform simulations. The simulations clearly demonstrate the induced wave propagation and water flow in the soil. The simulation technique is applicable to railway as well as road traffic. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
15.
孔隙流体压力与流体排驱的关系 总被引：2，自引：3，他引：2
本文讨论了在压实作用下孔隙流体压力的形成以及与流体排驱的关系，指出异常孔隙流体压力是因岩石渗透率变化引起毛细管力增加而产生的，且二者之间在流体压力孕育过程中一直维持一个动态平衡状态。对于连续沉积的盆地，只有当异常孔隙流体压力增加到超过岩石的抗剪强度时，因岩石发生剪切破裂导致毛细管力降低，流体才被排出；对于强烈构造变动的盆地，因地层大量剥蚀引起负荷压力降低，其降低幅度达到或超过岩石的抗张强度时，岩石 相似文献
16.
孔隙液体的可压缩性对非饱和孔隙介质弹塑性分析的影响 总被引：2，自引：0，他引：2
首先在多相孔隙介质的混合物理论框架内探讨了如何利用固-液两相孔隙介质模型描述介质的非饱和问题。其次,讨论了一个实用的描述孔隙液体可压缩性的方法,并建立了小变形条件下固-液两相孔隙介质的弹塑性本构方程。最后,采用俞茂宏统一强度理论,针对一个具体的结构例子,分析了孔隙液体的可压缩性对结构弹塑性分析的影响。 相似文献
17.
Visualizing complex pore structure and fluid flow in porous media using 3D printing technology and LBM simulation 
下载免费PDF全文
下载免费PDF全文 Pore structures of porous media and properties of fluid flow are key factors for the study of non-Darcy groundwater flow. However, it is difficult to directly observe pore structures and flow properties, resulting in a “black box” problem of porous media. This problem has hindered the in-depth study of the groundwater flow mechanism at the pore scale. In recent years, 3D rapid prototyping technology has seen tremendous development. 3D printing provides digital models and printing models of porous media with clear internal structure. Thus, Lattice Boltzmann Method can be used to simulate the flow processes at the pore scale based on real pore structures. In this study, 3D printing cores and Lattice Boltzmann Method were coupled to conduct both laboratory and numerical experiments in spherical porous media with different sphere diameters and periodic arrays. The LBM simulation results show a good agreement with laboratory experimental results. With the advantages of LBM and 3D printing, this approach provides a visualization of the complex pore structure and fluid flow in pores, which is a promising method for studies of non-Darcy groundwater flow at the pore scale. 相似文献
18.
Pirnia Pouyan Duhaime François Ethier Yannic Dubé Jean-Sébastien 《Acta Geotechnica》2020,15(10):2877-2889
This paper presents a coupled finite and discrete-element model (FEM and DEM) to simulate internal erosion. The model is based on ICY, an interface between COMSOL, an FEM engine, and YADE, a DEM code. With this model, smaller DEM subdomains are generated to simulate particle displacements at the grain scale. Particles in these small subdomains are subjected to buoyancy, gravity, drag and contact forces for short time steps (0.1 s). The DEM subdomains provide the macroscale (continuum) model with a particle flux distribution. Through a mass conservation equation, the flux distribution allows changes in porosity, hydraulic conductivity and hydraulic gradient to be evaluated for the same time steps at a larger, continuum scale. The updated hydraulic gradients from the continuum model provide the DEM subdomains with updated hydrodynamic forces based on a coarse-grid method. The number of particles in the DEM subdomains is also updated based on the new porosity distribution. The hierarchical multiscale model (HMM) was validated with the simulation of suffusion. Results for the proposed HMM algorithm are consistent with results based on a DEM model incorporating the full sample and simulation duration. The proposed HMM algorithm could enable the modelling of internal erosion for soil volumes that are too large to be modelled with a single DEM subdomain.
相似文献19.
A coupled discrete element model for the simulation of soil and water flow through an orifice 下载免费PDF全文
下载免费PDF全文 Soil erosion around defective underground pipes can cause ground collapses and sinkholes in urban areas. Most of these soil erosion events are caused by fluidization of the surrounding soil with subsequent washing into defective sewer pipes. In this study, this soil erosion process is simplified as the gradual washout of sand particles mixed with water through an orifice. The discrete element method is used to simulate the large deformation behavior of the sand particles, and the Darcy fluid model is coupled with this approach to simulate fluid flow through porous sand media. A coupled 3D discrete element model is developed and implemented based on this scheme. To simulate previous experiments using this coupled model considering the current computing capacity, we incorporated a ‘supply layer’ to study the continuous erosion process. The coupled model can predict the erosion flow rates of sand and water and the shape of erosion void. Thus, the model can be used as an effective and efficient tool to investigate the soil erosion process around defective pipes. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献
20.
Marine sedimentary basin is an important tectonic unit in the earth, and the evolution of marine sedimentary basin involves a series of the coupling and evolution of geodynamic mechanism such as the crust-mantle, the lithosphere-asthenosphere, the strata-fluid deposition. Therefore, the study of marine sedimentary basin dynamics includes deep structure state of earth, material composition and regional tectonic evolution, and also internal structure, tectonic characteristics and pore fluid characteristics strata of the basin. Wide angle Ocean Bottom Seismometer (OBS) seismic exploration is a marine geophysical survey method originated and developed since 1980’s and 1990’s, which has the advantages of strong penetration capability, high seismic imaging precision and reception of both P-wave and S-wave, and playing an increasing significant role in the research of marine sedimentary basin at the aspect of regional tectonic evolution, internal structure and pore fluid development characteristics of strata in recent years. In the study of passive continental margin, the crustal structure acquired from wide angle OBS seismic data provides the direct evidence that divides the passive continental margins into magma-poor and magma-dominated ones, and the degree of thinning and spatial variation characteristics of crust provide constraints for dynamics simulations of tectonic evolution in marine sedimentary basin. In the study of the structure features of basin, wide angle OBS seismic exploration fills in gaps at the aspect of investigation depth and complex geological structure in conventional multi-channel seismic survey, and acquires overpressure distribution status of basin according to the velocity structure characteristics of strata, and then infer the basin sedimentation velocity and pore fluid characteristics. In the study of internal fluid system in marine sedimentary basin, it reveals the velocity structure of natural gas hydrate reservoir through the analysis and processing of wide angle OBS seismic data, and calculates the thickness of natural gas hydrate reservoir and the content of hydrate and free gas in pore based on the velocity variation. Of course, the future wide angle OBS seismic exploration in the aspect of dynamic evolution and mechanism research in marine sedimentary basin will play a more important role with the development of marine seismic exploration technology, the improvement of data processing and instrument parameters. 相似文献