查看更多>>摘要：The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a preliminary proppant settling test,from which a solid volume fraction of 0.575 is calibrated for the proppant pack in the fracture.In the established workflow to investigate proppant flowback,a displacement is applied to the fracture surfaces to compact the generated proppant pack as well as further mimicking proppant embedment under closure stress.When a pressure gradient is applied to drive the fluid-particle flow,a critical aperture-to-diameter ratio of 4 is observed,above which the proppant pack would collapse.The results also show that the volumetric proppant flowback rate increases quadratically with the fracture aperture,while a linear variation between the particle flux and the pressure gradient is exhibited for a fixed fracture aperture.The research outcome contributes to-wards an improved understanding of proppant flowback in hydraulic fractures,which also supports an optimised proppant size selection for hydraulic fracturing operations.
查看更多>>摘要：Rockburst are often encountered in tunnel construction due to the complex geological conditions.To study the influence of unloading rate on rockburst,gneiss rockburst experiments were conducted under three groups of unloading rates.A high-speed photography system and acoustic emission(AE)system were used to monitor the entire process of rockburst process in real-time.The results show that the intensity of gneiss rockburst decreases with decrease of unloading rate,which is manifested as the reduction of AE energy and fragments ejection velocity.The mechanisms are proposed to explain this effect:(ⅰ)The reduction of unloading rate changes the crack propagation mechanism in the process of rockburst.This makes the rockbursts change from the tensile failure mechanism at high unloading rate to the tension-shear mixed failure mechanism at low unloading rate,and more energy released in the form of shear crack propagation.Then,less strain energy is converted into kinetic energy of fragments ejec-tion.(ⅱ )Less plate cracking degree of gneiss has taken shape due to decrease of unloading rate,resulting in the destruction of rockburst incubation process.The enlightenments of reducing the unloading rate for the project are also described quantitatively.The rockburst magnitude is reduced from the medium magnitude at the unloading rate of 0.1 MPa/s to the slight magnitude at the unloading rate of 0.025 MPa/s,which was judged by the ejection velocity.
查看更多>>摘要：This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accu-rately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.
查看更多>>摘要：In the context of repositories for nuclear waste,understanding the behavior of gas migration through clayey rocks with inherent anisotropy is crucial for assessing the safety of geological disposal facilities.The primary mechanism for gas breakthrough is the opening of micro-fractures due to high gas pressure.This occurs at gas pressures lower than the combined strength of the rock and its minimum principal stress under external loading conditions.To investigate the mechanism of microscale mode-Ⅰ ruptures,it is essential to incorporate a multiscale approach that includes subcritical microcracks in the modeling framework.In this contribution,we derive the model from microstructures that contain periodically distributed microcracks within a porous material.The damage evolution law is coupled with the macroscopic poroelastic system by employing the asymptotic homogenization method and considering the inherent hydro-mechanical(HM)anisotropy at the microscale.The resulting permeability change induced by fracture opening is implicitly integrated into the gas flow equation.Verification examples are presented to validate the developed model step by step.An analysis of local macroscopic response is undertaken to underscore the influence of factors such as strain rate,initial damage,and applied stress,on the gas migration process.Numerical examples of direct tension tests are used to demonstrate the model's efficacy in describing localized failure characteristics.Finally,the simulation results for prefer-ential gas flow reveal the robustness of the two-scale model in explicitly depicting gas-induced frac-turing in anisotropic clayey rocks.The model successfully captures the common behaviors observed in laboratory experiments,such as a sudden drop in gas injection pressure,rapid build-up of downstream gas pressure,and steady-state gas flow following gas breakthrough.
查看更多>>摘要：Geomechanical assessment using coupled reservoir-geomechanical simulation is becoming increasingly important for analyzing the potential geomechanical risks in subsurface geological developments.However,a robust and efficient geomechanical upscaling technique for heterogeneous geological res-ervoirs is lacking to advance the applications of three-dimensional(3D)reservoir-scale geomechanical simulation considering detailed geological heterogeneities.Here,we develop convolutional neural network(CNN)proxies that reproduce the anisotropic nonlinear geomechanical response caused by lithological heterogeneity,and compute upscaled geomechanical properties from CNN proxies.The CNN proxies are trained using a large dataset of randomly generated spatially correlated sand-shale re-alizations as inputs and simulation results of their macroscopic geomechanical response as outputs.The trained CNN models can provide the upscaled shear strength(R2＞0.949),stress-strain behavior(R2＞0.925),and volumetric strain changes(R2＞0.958)that highly agree with the numerical simulation results while saving over two orders of magnitude of computational time.This is a major advantage in computing the upscaled geomechanical properties directly from geological realizations without the need to perform local numerical simulations to obtain the geomechanical response.The proposed CNN proxy-based upscaling technique has the ability to(1)bridge the gap between the fine-scale geocellular models considering geological uncertainties and computationally efficient geomechanical models used to assess the geomechanical risks of large-scale subsurface development,and(2)improve the efficiency of nu-merical upscaling techniques that rely on local numerical simulations,leading to significantly increased computational time for uncertainty quantification using numerous geological realizations.
查看更多>>摘要：Porous materials present significant advantages for absorbing radioactive isotopes in nuclear waste streams.To improve absorption efficiency in nuclear waste treatment,a thorough understanding of the diffusion-advection process within porous structures is essential for material design.In this study,we present advancements in the volumetric lattice Boltzmann method(VLBM)for modeling and simulating pore-scale diffusion-advection of radioactive isotopes within geopolymer porous structures.These structures are created using the phase field method(PFM)to precisely control pore architectures.In our VLBM approach,we introduce a concentration field of an isotope seamlessly coupled with the velocity field and solve it by the time evolution of its particle population function.To address the computational intensity inherent in the coupled lattice Boltzmann equations for velocity and concentration fields,we implement graphics processing unit(GPU)parallelization.Validation of the developed model involves examining the flow and diffusion fields in porous structures.Remarkably,good agreement is observed for both the velocity field from VLBM and multiphysics object-oriented simulation environment(MOOSE),and the concentration field from VLBM and the finite difference method(FDM).Furthermore,we investigate the effects of background flow,species diffusivity,and porosity on the diffusion-advection behavior by varying the background flow velocity,diffusion coefficient,and pore volume fraction,respectively.Notably,all three parameters exert an influence on the diffusion-advection process.Increased background flow and diffusivity markedly accelerate the process due to increased advection intensity and enhanced diffusion capability,respectively.Conversely,increasing the porosity has a less significant effect,causing a slight slowdown of the diffusion-advection process due to the expanded pore volume.This comprehensive parametric study provides valuable insights into the kinetics of isotope uptake in porous structures,facilitating the development of porous materials for nuclear waste treatment applications.
查看更多>>摘要：A nuclear explosion in the rock mass medium can produce strong shock waves,seismic shocks,and other destructive effects,which can cause extreme damage to the underground protection infrastructures.With the increase in nuclear explosion power,underground protection engineering enabled by explosion-proof impact theory and technology ushered in a new challenge.This paper proposes to simulate nuclear explosion tests with on-site chemical explosion tests in the form of multi-hole ex-plosions.First,the mechanism of using multi-hole simultaneous blasting to simulate a nuclear explosion to generate approximate plane waves was analyzed.The plane pressure curve at the vault of the un-derground protective tunnel under the action of the multi-hole simultaneous blasting was then obtained using the impact test in the rock mass at the site.According to the peak pressure at the vault plane,it was divided into three regions:the stress superposition region,the superposition region after surface reflection,and the approximate plane stress wave zone.A numerical simulation approach was developed using PFC and FLAC to study the peak particle velocity in the surrounding rock of the underground protective cave under the action of multi-hole blasting.The time-history curves of pressure and peak pressure partition obtained by the on-site multi-hole simultaneous blasting test and numerical simu-lation were compared and analyzed,to verify the correctness and rationality of the formation of an approximate plane wave in the simulated nuclear explosion.This comparison and analysis also provided a theoretical foundation and some research ideas for the ensuing study on the impact of a nuclear explosion.
查看更多>>摘要：Natural slopes usually display complicated exposed rock surfaces that are characterized by complex and substantial terrain undulation and ubiquitous undesirable phenomena such as vegetation cover and rockfalls.This study presents a systematic outcrop research of fracture pattern variations in a compli-cated rock slope,and the qualitative and quantitative study of the complex phenomena impact on three-dimensional(3D)discrete fracture network(DFN)modeling.As the studies of the outcrop fracture pattern have been so far focused on local variations,thus,we put forward a statistical analysis of global variations.The entire outcrop is partitioned into several subzones,and the subzone-scale variability of fracture geometric properties is analyzed(including the orientation,the density,and the trace length).The results reveal significant variations in fracture characteristics(such as the concentrative degree,the average orientation,the density,and the trace length)among different subzones.Moreover,the density of fracture sets,which is approximately parallel to the slope surface,exhibits a notably higher value compared to other fracture sets across all subzones.To improve the accuracy of the DFN modeling,the effects of three common phenomena resulting from vegetation and rockfalls are qualitatively analyzed and the corresponding quantitative data processing solutions are proposed.Subsequently,the 3D frac-ture geometric parameters are determined for different areas of the high-steep rock slope in terms of the subzone dimensions.The results show significant variations in the same set of 3D fracture parameters across different regions with density differing by up to tenfold and mean trace length exhibiting dif-ferences of 3-4 times.The study results present precise geological structural information,improve modeling accuracy,and provide practical solutions for addressing complex outcrop issues.
查看更多>>摘要：In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying satu-rated soil dynamics.
查看更多>>摘要：An elastoplastic constitutive model based on the Modified Cam Clay(MCC)model is developed to describe the mechanical behaviour of soils cemented via microbially induced calcite precipitation(MICP).It considers the increase of the elastic stiffness,the change of the yield surface due to MICP cementation and the degradation of calcium carbonate bonds during shearing.Specifically,to capture the typical contraction-dilation transition in MICP soils,the original volumetric hardening rule in the MCC model is modified to a combined deviatoric and volumetric hardening rule.The model could reproduce a series of drained triaxial tests on MICP-treated soils with different calcium carbonate contents.Further,we carry out a parametric study and observe numerical instability in some cases.In combination with an analytical analysis,our numerical modelling has identified the benefits and limitations of using MCC-based models in the simulation of MICP-cemented soils,leading to suggestions for further model development.
NETL
NSTL
万方数据