An adaptive finite element model for wellbore breakout simulation

A finite element formulation is proposed and implemented for analysing the stability of excavated wells using the DiMaggio-Sandler constitutive elastoplastic model with a typical carbonate reservoir configuration. The quality of the finite element approximation is ensured by applying smooth curved elements adapted to the wellbore geometry, and h − p adaptive finite element meshes in the plastic zone. General purpose procedures are defined to transfer the elastoplastic deformation history to newly created integration points. A breakout damage criterion is proposed based on the second invariant of the deviatoric plastic deformation tensor. This damage criterion is used to apply a mesh movement algorithm to represent material collapse. The automatic successive application of the breakout damage criterion results in elliptical realistically looking geometries obtained in experiments reported in the literature.     

Critical Top Tension for Static Equilibrium Configuration of A Steel Catenary Riser

This paper aims to present the critical top tension for static equilibrium configurations of a steel catenary riser（SCR） by using the finite element method. The critical top tension is the minimum top tension that can maintain the equilibrium of the SCR. If the top tension is smaller than the critical value, the equilibrium of the SCR does not exist. If the top tension is larger than the critical value, there are two possible equilibrium configurations. These two configurations exhibit the nonlinear large displacement. The configuration with the smaller displacement is stable, while the one with larger displacement is unstable. The numerical results show that the increases in the riser＇s vertical distances, horizontal offsets, riser＇s weights, internal flow velocities, and current velocities increase the critical top tensions of the SCR. In addition, the parametric studies are also performed in order to investigate the limit states for the analysis and design of the SCR.     

Improved algorithms for generalized finite element simulations of three-dimensional hydraulic fracture propagation

This paper reports improvements to algorithms for the simulation of 3-D hydraulic fracturing with the Generalized Finite Element Method (GFEM). Three optimizations are presented and analyzed. First, an improved initial guess based on solving a 3-D elastic problem with the pressure from the previous step is shown to decrease the number of Newton iterations and increase robustness. Second, an improved methodology to find the time step that leads to fracture propagation is proposed and shown to decrease significantly the number of iterations. Third, reduced computational cost is observed by properly recycling the linear part of the coupled stiffness matrix. Two representative examples are used to analyze these improvements. Additionally, a methodology to include the leak-off term is presented and verified against asymptotic analytical solutions. Conservation of mass is shown to be well satisfied in all examples.     

On the pore‐scale mechanisms leading to brittle and ductile deformation behavior of crystalline rocks

Deformation mechanisms at the pore scale are responsible for producing large strains in porous rocks. They include cataclastic flow, dislocation creep, dynamic recrystallization, diffusive mass transfer, and grain boundary sliding, among others. In this paper, we focus on two dominant pore‐scale mechanisms resulting from purely mechanical, isothermal loading: crystal plasticity and crofracturing. We examine the contributions of each mechanism to the overall behavior at a scale larger than the grains but smaller than the specimen, which is commonly referred to as the mesoscale. Crystal plasticity is assumed to occur as dislocations along the many crystallographic slip planes, whereas microfracturing entails slip and frictional sliding on microcracks. It is observed that under combined shear and tensile loading, microfracturing generates a softer response compared with crystal plasticity alone, which is attributed to slip weakening where the shear stress drops to a residual level determined by the frictional strength. For compressive loading, however, microfracturing produces a stiffer response than crystal plasticity because of the presence of frictional resistance on the slip surface. Behaviors under tensile, compressive, and shear loading invariably show that porosity plays a critical role in the initiation of the deformation mechanisms. Both crystal plasticity and microfracturing are observed to initiate at the peripheries of the pores, consistent with results of experimental studies. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical and analytical observations on long and infinite slopes

Interest in the mechanics of landslides has led to renewed evaluation of the infinite slope equations, and the need for a more general framework for estimating the factor of safety of long and infinite slopes involving non‐homogeneous soil profiles. The paper describes finite element methods that demonstrate the potential for predicting failure in long slope profiles where the critical mechanism is not necessarily at the base of the soil layer. The influence of slope angle is also examined in long slopes, leading to some counter‐intuitive conclusions about the impact of slope steepness on the factor of safety. Copyright © 2010 John Wiley & Sons, Ltd.     

Variations in elastic thickness and flexure of the Maracaibo block

We estimate the lateral variations of the elastic thickness of the Maracaibo block with a 3D numerical approach by using centered finite differences. The calculation is based on solving the fourth-order partial differential equation that governs the bending of a thin plate fixed on its boundaries (zero displacement) with variable thickness (or elastic thickness for this particular case). An initial plate-load model is built and is iteratively modified to fit the general basement configuration and gravity data. The final result is an elastic thickness map that covers the Maracaibo block and the surrounding sections of the South American plate. It shows that the elastic thickness ranges from 30 km to 18 km with a mean value of 23.6 km and a mode of 26 km. The largest elastic thickness values are associated with the location of the Santa Marta Mountains and the Barinas Apure Basin, while the smallest ones with the Mérida Andes-Maracaibo Basin flexural system. The current basement configuration within the Maracaibo basin, formed as a result of its geodynamic evolution, has affected the mechanical properties of the Maracaibo block near the current Mérida Andes position. The load of the Perijá Range is compensated by a complex stress tensor, and that of the Santa Marta Mountains does not have an isostatic root as it is held by a relatively strong lithosphere.     

A discontinuous finite element suspended sediment transport model for water quality assessments in river networks

Suspended sediment plays an important role in the distribution and transport of many pollutants (such as radionuclides) in rivers. Pollutants may adsorb on fine suspended particles (e.g. clay) and spread according to the suspended sediment movement. Hence, the simulation of the suspended sediment mechanism is indispensable for realistic transport modelling. This paper presents and tests a simple mathematical model for predicting the suspended sediment transport in river networks. The model is based on the van Rijn suspended load formula and the advection–diffusion equation with a source or sink term that represents the erosion or deposition fluxes. The transport equation is solved numerically with the discontinuous finite element method. The model evaluation was performed in two steps, first by comparing model simulations with the measured suspended sediment concentrations in the Grote Nete–Molse Nete River in Belgium, and second by a model intercomparison with the sediment transport model NST MIKE 11. The simulations reflect the measurements with a Nash‐Sutcliffe model efficiency of 0.6, while the efficiency between the proposed model and the NST MIKE 11 simulations is 0.96. Both evaluations indicate that the proposed sediment transport model, that is sufficiently simple to be practical, is providing realistic results.     

Analysis of shaft resistance of jacked piles in sands

In recent years, pile jacking has become a viable alternative installation method for displacement piles. Pile jacking produces minimal noise, vibration and air pollution during installation. In addition, it is possible, at the end of jacking, to have a good estimate of the ultimate static capacity of the pile. In this paper, the shaft resistance of piles jacked into sand is studied using one‐dimensional finite element analysis. The finite element simulations, using a two‐surface plasticity model, demonstrate the effects of relative density and confinement on the unit shaft resistance of piles jacked in sand. The impact of the number of jacking strokes on the unit shaft capacity is also assessed. Based on the numerical results, we developed equations for shaft resistance quantifying the effects of relative density, initial confinement and number of jacking strokes. Predictions using these equations are compared with data obtained from centrifuge tests and field tests. Copyright © 2010 John Wiley & Sons, Ltd.     

Steady‐state analytical solutions of flow and deformation coupling due to a point sink within a finite fluid‐saturated poroelastic layer

An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.