Modelling erosion of a single rock block using a coupled CFD-DEM approach

Rock block removal is the prevalent physical mechanism for rock erosion and could affect the stability of dam foundations and spillways.Despite this,understanding of block removal is still inadequate because of the complex interactions among block characteristics,hydraulic forces,and erosive processes acting on the block.Herein,based on a previously conducted physical experiment of erosion of a single rock block,the removal processes of two different protruding blocks are represented by a coupled computational fluid dynamics-discrete element model(CFD-DEM)approach under varied flow conditions.Additionally,the blocks could be rotated with respect to the flow direction to consider the effect of the discontinuity orientation on the block removal process.Simulation results visualize the entire block removal process.The simulations reproduce the effects of the discontinuity orientation on the critical flow velocity inducing block incipient motion and the trajectory of the block motion observed in the physical experiments.The numerical results present a similar tendency of the critical velocities at different discontinuity orientations but have slightly lower values.The trajectory of the block in the simulations fits well with the experimental measurements.The relationship between the dimensionless critical shear stress and discontinuity orientation observed from the simulations shows that the effect of block protrusion becomes more dominant on the block incipient motion with the increase of relative protrusion height.To our knowledge,this present study is the first attempt to use the coupled finite volume method(FVM)-DEM approach for modelling the interaction behavior between the block and the flowing water so that the block removal process can be reproduced and analyzed.

Mechanisms of shale hydration and water block removal

Shale samples of Silurian Longmaxi Formation in the Changning area of the Sichuan Basin, SW China, were selected to carry out scanning electron microscopy, CT imaging, high-pressure mercury injection, low-temperature nitrogen adsorption and imbibition experiments to compare the hydration characteristics of montmorillonite and illite, analyze the main factors affecting the water block removal of shale, and reveal the mechanisms of pore structure evolution during shale hydration. The hydration characteristics of shale are closely related to the composition of clay minerals, the shale with high illite content is not susceptible to hydration and thus has limited room for pore structure improvement;the shale with high montmorillonite is susceptible to hydration expansion and thus has higher potential of pore structure improvement by stimulation;the shale with high illite content has stronger imbibition in the initial stage, but insufficient diffusion ability, and thus is likely to have water block;the shale with high montmorillonite content has weaker imbibition in the initial stage but better water diffusion, so water blocking in this kind of shale can be removed to some degree;the shale reservoir has an optimal hydration time, when it is best in physical properties, but hydration time too long would cause damage to the reservoir, and the shale with high illite content has a shorter optimal hydration time;inorganic cations can inhibit the hydration of clay minerals and have stronger inhibition to illite expansion, especially K^(+);for the reservoir with high content of montmorillonite, the cation content of fracturing fluid can be lowered to promote the shale hydration;fracturing fluid with high K^(+) content can be injected into reservoirs with high illite content to suppress hydration.