Using fracture-based continuum modeling of coupled geomechanical-hydrological processes for numerical simulation of hydraulic fracturing

This paper describes numerical simulation of hydraulic fracturing using fracture-based continuum modeling(FBCM)of coupled geomechanical-hydrological processes to evaluate a technique for high-density fracturing and fracture caging.The simulations are innovative because of modeling discrete fractures explicitly in continuum analysis.A key advantage of FBCM is that fracture initiation and propagation are modeled explicitly without changing the domain grid(i.e.no re-meshing).Further,multiple realizations of a preexisting fracture distribution can be analyzed using the same domain grid.The simulated hydraulic fracturing technique consists of pressurizing multiple wells simultaneously:initially without permeating fluids into the rock,to seed fractures uniformly and at high density in the wall rock of the wells;followed by fluid injection to propagate the seeded fracture density hydraulically.FBCM combines the ease of continuum modeling with the potential accuracy of modeling discrete fractures and fracturing explicitly.Fractures are modeled as piecewise planar based on intersections with domain elements;fracture geometry stored as continuum properties is used to calculate parameters needed to model individual fractures;and rock behavior is modeled through tensorial aggregation of the behavior of discrete fractures and unfractured rock.Simulations are presented for previously unfractured rock and for rock with preexisting fractures of horizontal,shallow-dipping,steeply dipping,or vertical orientation.Simulations of a single-well model are used to determine the pattern and spacing for a multiple-well design.The results illustrate high-density fracturing and fracture caging through simultaneous fluid injection in multiple wells:for previously unfractured rock or rock with preexisting shallow-dipping or horizontal fractures,and in situ vertical compressive stress greater than horizontal.If preexisting fractures are steeply dipping or vertical,and considering the same in situ stress condition,well pressurization without fluid permeation appears to be the only practical way to induce new fractures and contain fracturing within the target domain.

Mechanical and hydraulic properties of fault rocks under multi‑stage cyclic loading and unloading

The rock mass in fault zones is frequently subjected to cyclic loading and unloading during deep resource exploitation and tunnel excavation.Research on the mechanical and hydraulic characteristics of fault rock during the cyclic loading and unloading is of great signifcance for revealing the formation mechanism of water-conducting pathways in fault and preventing water inrush disasters.In this study,the mechanical and seepage tests of fault rock under the multi-stage cyclic loading and unloading of axial compression were carried out by using the fuid–solid coupling triaxial experimental device.The hysteresis loop of the stress–strain curve,peak strain rate,secant Young's modulus,and permeability of fault rock were obtained,and the evolution law of the dissipated energy of fault rock with the cyclic number of load and unloading was discussed.The experimental results show that with an increase in the cyclic number of loading and unloading,several changes occur.The hysteresis loop of the stress–strain curve of the fault rock shifts towards higher levels of strain.Additionally,both the peak strain rate and the secant Young's modulus of the fault rock increase,resulting in an increase in the secant Young's modulus of the fault rock mass.However,the growth rate of the secant Young's modulus gradually slows down with the increase of cyclic number of loading and unloading.The permeability evolution of fault rock under the multi-stage cyclic loading and unloading of axial compression can be divided into three stages:steady increase stage,cyclic decrease stage,and rapid increase stage.Besides,the calculation model of dissipated energy of fault rock considering the efective stress was established.The calculation results show that the relationship between the dissipated energy of fault rock and the cyclic number of loading and unloading conforms to an exponential function.

Comparative assessment of mechanical and chemical fluid diversion techniques during hydraulic fracturing in horizontal wells

The application of fluid diversion during hydraulic fracturing is an evolving technology and has become popular amongst E&P operators over the past few years.The primary objective of the fluid diversion is to improve hydraulic fracturing treatment by increasing stimulated reservoir volume and improving hydrocarbon recovery.This is possible by achieving any of the following objectives:creating uniform distribution of treatment slurry within the target zone;treating unstimulated and under-stimulated zones;or by increasing fracture density by creating a complex fracture network.The fluid diversion application is also helpful in decreasing the number of stages(by increasing stage length)for multi-stage plug-n-perf(PnP)fracturing treatment.It is also applied to prevent fracture-driven interactions between adjacent wells,which is currently a major issue,especially in shale.In addition,for successful refracturing treatment,the diverter application is essential for isolating the existing fractures and redirecting the treatment slurry to the desired unstimulated zones.The diversion methods can be broadly categorized into the mechanical and chemical diversion.Several established mechanical diversion techniques are frac plugs,expandable casing patches,expandable liners,swellable packers,straddle packer assembly,sand plugs,frac sleeves,perforation ball sealers,and limited entry technique.The different chemical diversion techniques are particulates,fibers,gels,surfactants,perforation pods,and composite diverting.This paper describes the current status of established mechanical and chemical diverter technologies and examines their comparative advantages and challenges.Various techniques are suitable for diverter application,but the technique is selected based on the desired objective and conditions of the wellbore and reservoir.The general guidelines for selecting diversion techniques and operational considerations are also provided in the paper.The diagnosis of diversion treatment plays an essential role in diversion technique selection and optimization of selection parameters for the subsequent treatments.Therefore,the application of conventional surface pressure monitoring techniques and advanced diagnostic tools to evaluate diversion effectiveness are briefly described.Presently no standard laboratory testing method is established for the performance evaluation of diverting agents.Therefore,researchers have implemented various laboratory methods,which are briefly summarized in the paper.Significant insight into the diversion technology and guidelines for its selection and successful implementation is provided to help engineers to increase the effectiveness of hydraulic fracturing treatments.The limitations of individual diversion techniques are clarified,which provide the future scope of research for improvement in various diversion technologies.

Wellbore Cleaning Degree and Hydraulic Extension in Shale Oil Horizontal Wells

The efficient development and exploitation of shale oil depends on long-distance horizontal wells. As the degreeof cleaning of the wellbore plays a key role in these processes, in this study, this problem is investigated experimentallyby focusing on the dimensionless cuttings bed height. A method is proposed to calculate the horizontalwellhydraulic extension taking into account the influence of the wellbore cleaning degree on the wellborepressure distribution and assess the effect of a variety of factors such as the bottom hole pressure, the circulatingpressure drop, the drilling pump performance and the formation properties. The analysis shows that the hydraulicextension of horizontal wells decreases with an increase in the cuttings bed height, and the higher the displacementof drilling fluid, the faster the hydraulic extension declines. The annular pressure drop of the horizontalsection increases with the increase of the cuttings bed height, resulting in a higher bottom-hole pressure. Severalarguments are provided to guide the safe drilling of shale oil horizontal wells and overcome the limits of currenttechnological approaches.

Hydraulic properties and drought response of a tropical bamboo (Cephalostachyum pergracile)

Bamboo plants are an essential component of tropical ecosystems,yet their vulnerability to climate extremes,such as drought,is poorly understood due to limited knowledge of their hydraulic properties.Cephalostachyum pergracile,a commonly used tropical bamboo species,exhibited a substantially higher mortality rate than other co-occurring bamboos during a severe drought event in 2019,but the underlying mechanisms remain unclear.This study investigated the leaf and stem hydraulic traits related to drought responses,including leaf-stem embolism resistance(P50leaf;P50stem) estimated using optical and X-ray microtomography methods,leaf pressure-volume and water-releasing curves.Additionally,we investigated the seasonal water potentials,native embolism level(PLC) and xylem water source using stable isotope.We found that C.pergracile exhibited strong resistance to embolism,showing low P50leaf,P50stem,and turgor loss point,despite its rapid leaf water loss.Interestingly,its leaves displayed greater resistance to embolism than its stem,suggesting a lack of effective hydraulic vulnerability segmentation(HVS) to protect the stem from excessive xylem tension.During the dry season,approximately 49% of the water was absorbed from the upper 20-cm-deep soil layer.Consequently,significant diurnal variation in leaf water potentials and an increase in midday PLC from 5.87±2.33% in the wet season to 12.87±4.09%in the dry season were observed.In summary,this study demonstrated that the rapid leaf water loss,high reliance on surface water,and a lack of effective HVS in C.pergracile accelerated water depletion and increased xylem embolism even in the typical dry season,which may explain its high mortality rate during extreme drought events in 2019.

Stress tensor determination by modified hydraulic tests on pre-existing fractures:Method and stress constraints

The hydraulic testing of pre-existing fractures(HTPF)is one of the most promising in situ stress measurement methods,particularly for three-dimensional stress tensor determination.However,the stress tensor determination based on the HTPF method requires at least six tests or a minimum of 14-15 tests(under different conditions)for reliable results.In this study,we modified the HTPF method by considering the shear stress on each pre-existing fracture,which increased the number of equations for the stress tensor determination and decreased the number of tests required.Different shear stresses were attributed to different fractures by random sampling;therefore,the stress tensors were obtained by searching for the optimal solution using the least squares criterion based on the Monte Carlo method.Thereafter,we constrained the stress tensor based on the tensile strength criterion,compressive strength criterion,and vertical stress constraints.The inverted stress tensors were presented and analyzed based on the tensorial nature of the stress using the Euclidean mean stress tensor.Two stress-measurement campaigns in Weifang(Shandong Province,China)and Mercantour road tunnel(France)were implemented to highlight the validity and efficiency of the modified HTPF(M-HTPF)method.The results showed that the M-HTPF method can be applied for stress tensor inversion using only three to four tests on pre-existing fractures,neglecting the stress gradient.The inversion results were confined to relatively small distribution dispersions and were significantly reliable and stable due to the shear stresses on the fractures and the stress constraints employed.The M-HTPF method is highly feasible and efficient for complete stress tensor determination in a single borehole.

Pumping-induced Well Hydraulics and Groundwater Budget in a Leaky Aquifer System with Vertical Heterogeneity in Aquitard Hydraulic Properties

In groundwater hydrology,aquitard heterogeneity is often less considered compared to aquifers,despite its significant impact on groundwater hydraulics and groundwater resources evaluation.A semi-analytical solution is derived for pumping-induced well hydraulics and groundwater budget with consideration of vertical heterogeneity in aquitard hydraulic conductivity(K)and specific storage(S_(s)).The proposed new solution is innovative in its partitioning of the aquitard into multiple homogeneous sub-layers to enable consideration of various forms of vertically heterogeneous K or S_(s).Two scenarios of analytical investigations are explored:one is the presence of aquitard interlayers with distinct K or S_(s) values,a common field-scale occurrence;another is an exponentially depth-decaying aquitard S_(s),a regional-scale phenomenon supported by statistical analysis.Analytical investigations reveal that a low-K interlayer can significantly increase aquifer drawdown and enhance aquifer/aquitard depletion;a high-S_(s) interlayer can noticeably reduce aquifer drawdown and increase aquitard depletion.Locations of low-K or high-S_(s) interlayers also significantly impact well hydraulics and groundwater budget.In the context of an exponentially depth-decaying aquitard S_(s),a larger decay exponent can enhance aquifer drawdown.When using current models with a vertically homogeneous aquitard,half the sum of the geometric and harmonic means of exponentially depth-decaying aquitard S_(s) should be used to calculate aquitard depletion and unconfined aquifer leakage.

Mechanical Behavior and Microstructure Evolution during Tensile Deformation of Twinning Induced Plasticity Steel Processed by Warm Forgings

The mechanical behavior and microstructural evolution of an Fe-30Mn-3Al-3Si twinninginduced plasticity(TWIP)steel processed using warm forging was investigated.It is found that steel processed via warm forging improves comprehensive mechanical properties compared to the TWIP steel processed via cold rolling,with a high tensile strength(R_(m))of 793 MPa,a yield strength(R_(P))of 682 MPa,an extremely large R_(P)/R_(m)ratio as high as 0.86 as well as an excellent elongation rate of 46.8%.The microstructure observation demonstrates that steel processed by warm forging consists of large and elongated grains together with fine,equiaxed grains.Complicated micro-defect configurations were also observed within the steel,including dense dislocation networks and a few coarse deformation twins.As the plastic deformation proceeds,the densities of dislocations and deformation twins significantly increase.Moreover,a great number of slip lines could be observed in the elongated grains.These findings reveal that a much more dramatic interaction between microstructural defect and dislocations glide takes place in the forging sample,wherein the fine and equiaxed grains propagated dislocations more rapidly,together with initial defect configurations,are responsible for enhanced strength properties.Meanwhile,larger,elongated grains with more prevalently activated deformation twins result in high plasticity.

Microdynamic mechanical properties and fracture evolution mechanism of monzogabbro with a true triaxial multilevel disturbance method

The far-field microdynamic disturbance caused by the excavation of deep mineral resources and underground engineering can induce surrounding rock damage in high-stress conditions and even lead to disasters.However,the mechanical properties and damage/fracture evolution mechanisms of deep rock induced by microdynamic disturbance under three-dimensional stress states are unclear.Therefore,a true triaxial multilevel disturbance test method is proposed,which can completely simulate natural geostress,excavation stress redistribution(such as stress unloading,concentration and rotation),and subsequently the microdynamic disturbance triggering damaged rock failure.Based on a dynamic true triaxial test platform,true triaxial microdynamic disturbance tests under different frequency and amplitudes were carried out on monzogabbro.The results show that increasing amplitude or decreasing frequency diminishes the failure strength of monzogabbro.Deformation modulus gradually decreases during disturbance failure.As frequency and amplitude increase,the degradation rate of deformation modulus decreases slightly,disturbance dissipated energy increases significantly,and disturbance deformation anisotropy strengthens obviously.A damage model has been proposed to quantitatively characterize the disturbance-induced damage evolution at different frequency and amplitude under true triaxial stress.Before disturbance failure,the micro-tensile crack mechanism is dominant,and the micro-shear crack mechanism increases significantly at failure.With the increase of amplitude and frequency,the micro-shear crack mechanism increases.When approaching disturbance failure,the acoustic emission fractal dimension changes from a stable value to local large oscillation,and finally increases sharply to a high value at failure.Finally,the disturbance-induced failure mechanism of surrounding rock in deep engineering is clearly elucidated.

Experimental Investigation on Vertical Hydraulic Transport of Ores in Deepsea Mining

Deepsea mining has been proposed since the 1960s to alleviate the lack of resources on land.Vertical hydraulic transport of collected ores from the seabed to the sea surface is considered the most promising method for industrial applications.In the present study,an indoor model test of the vertical hydraulic transport of particles was conducted.A noncontact optical method has been proposed to measure the local characteristics of the particles inside a vertical pipe,including the local concentration and particle velocity.The hydraulic gradient of ore transport was evaluated with various particle size distributions,particle densities,feeding concentrations and mixture flow velocities.During transport,the local concentration is larger than the feeding concentration,whereas the particle velocity is less than the mixture velocity.The qualitative effects of the local concentration and local fluid velocity on the particle velocity and slip velocity were investigated.The local fluid velocity contributes significantly to particle velocity and slip velocity,whereas the effect of the local concentration is marginal.A higher feeding concentration and mixture flow velocity result in an increased hydraulic gradient.The effect of the particle size gradation is slight,whereas the particle density plays a crucial role in the transport.

Effects of seepage pressure on the mechanical behaviors and microstructure of sandstone

Surrounding rocks of underground engineering are subjected to long-term seepage pressure,which can deteriorate the mechanical properties and cause serious disasters.In order to understand the impact of seepage pressure on the mechanical property of sandstone,uniaxial compression tests,P-wave velocity measurements,and nuclear magnetic resonance(NMR)tests were conducted on saturated sandstone samples with varied seepage pressures(i.e.0 MPa,3 MPa,4 MPa,5 MPa,6 MPa,7 MPa).The results demonstrate that the mechanical parameters(uniaxial compressive strength,peak strain,elastic modulus,and brittleness index),total energy,elastic strain energy,as well as elastic strain energy ratio,decrease with increasing seepage pressure,while the dissipation energy and dissipation energy ratio increase.Moreover,as seepage pressure increases,the micro-pores gradually transform into meso-pores and macro-pores.This increases the cumulative porosity of sandstone and decreases P-wave velocity.The numerical results indicate that as seepage pressure rises,the number of tensile cracks increases progressively,the angle range of microcracks is basically from 50-120to 80-100,and as a result,the failure mode transforms to the tensile-shear mixed failure mode.Finally,the effects of seepage pressure on mechanical properties were discussed.The results show that decrease in the effective stress and cohesion under the action of seepage pressure could lead to deterioration of strength behaviors of sandstone.

Intermittent disturbance mechanical behavior and fractional deterioration mechanical model of rock under complex true triaxial stress paths

Mechanical excavation,blasting,adjacent rockburst and fracture slip that occur during mining excavation impose dynamic loads on the rock mass,leading to further fracture of damaged surrounding rock in three-dimensional high-stress and even causing disasters.Therefore,a novel complex true triaxial static-dynamic combined loading method reflecting underground excavation damage and then frequent intermittent disturbance failure is proposed.True triaxial static compression and intermittent disturbance tests are carried out on monzogabbro.The effects of intermediate principal stress and amplitude on the strength characteristics,deformation characteristics,failure characteristics,and precursors of monzogabbro are analyzed,intermediate principal stress and amplitude increase monzogabbro strength and tensile fracture mechanism.Rapid increases in microseismic parameters during rock loading can be precursors for intermittent rock disturbance.Based on the experimental result,the new damage fractional elements and method with considering crack initiation stress and crack unstable stress as initiation and acceleration condition of intermittent disturbance irreversible deformation are proposed.A novel three-dimensional disturbance fractional deterioration model considering the intermediate principal stress effect and intermittent disturbance damage effect is established,and the model predicted results align well with the experimental results.The sensitivity of stress states and model parameters is further explored,and the intermittent disturbance behaviors at different f are predicted.This study provides valuable theoretical bases for the stability analysis of deep mining engineering under dynamic loads.

Jet Characteristics and Optimization of a Cavitation Nozzle for Hydraulic Fracturing Applications

Hydraulic jetting is a form of fracturing that involves using a high-pressure jet of water to create fractures in the reservoir rock with a nozzle serving as the central component of the hydraulic sandblasting perforation tool.In this study,the flow behavior of the nozzle is simulated numerically in the framework of a SST k-ωturbulence model.The results show that the nozzle structure can significantly influence the jet performance and related cavitation effect.Through orthogonal experiments,the nozzle geometric parameters are optimized,and the following configuration is found accordingly:contraction angle 20°,contraction segment length 6 mm,cylindrical segment diameter 6 mm,cylindrical segment length 12 mm,spread segment length 10 mm,and spread angle 55°.

Unraveling the hydraulic properties of loess for landslide prediction:A study on variations in loess landslides in Lanzhou,Dingxi,and Tianshui,China

Loess has distinctive characteristics,leading to frequent landslide disasters and posing serious threats to the lives and properties of local re sidents.The involvement of water repre sents a critical factor in inducing loess landslides.This study focuses on three neighboring cities sequentially situated on the Loess Plateau along the direction of aeolian deposition of loess,namely Lanzhou,Dingxi,and Tianshui,which are densely populated and prone to landslide disasters.The variations in hydraulic properties,including water retention capacity and permeability,are investigated through Soil Water Characteristic Curve(SWCC)test and hydraulic conductivity test.The experimental findings revealed that Tianshui loess exhibited the highest water retention capacity,followed by Dingxi loess,while Lanzhou loess demonstrated the lowest water retention capacity.Contrastingly,the results for the saturated permeability coefficient were found to be the opposite:Tianshui loess showed the lowest permeability,whereas Lanzhou loess displayed the highest permeability.These results are supported and analyzed by scanning electron microscopy(SEM)observation.In addition,the water retention capacity is mathematically expressed using the van Genuchten model and extended to predict unsaturated hydraulic properties of loess.The experimental results exhibit a strong accordance with one another and align with the regional distribution patterns of disasters.

Spatiotemporal variations of sand hydraulic conductivity by microbial application methods

The spatiotemporal distributions of microbes in soil by different methods could affect the efficacy of the microbes to reduce the soil hydraulic conductivity.In this study,the specimens of bio-mediated sands were prepared using three different methods,i.e.injecting,mixing,and pouring a given microbial so-lution onto compacted sand specimens.The hydraulic conductivity was measured by constant-head tests,while any soil microstructural changes due to addition of the microbes were observed by scan-ning electron microscope(SEM)and mercury intrusion porosimetry(MIP)tests.The amount of dextran concentration produced by microbes in each type of specimen was quantified by a refractometer.Results show that dextran production increased exponentially after 5-7 d of microbial settling with the supply of culture medium.The injection and mixing methods resulted in a similar amount and uniform dis-tribution of dextran in the specimens.The pouring method,however,produced a nonuniform distri-bution,with a higher concentration near the specimen surface.As the supply of culture medium discontinued,the dextran content near the surface produced by the pouring method decreased dramatically due to high competition for nutrients with foreign colonies.Average dextran concentration was negatively and correlated with hydraulic conductivity of bio-mediated soils exponentially,due to the clogging of large soil pores by dextran.The hydraulic conductivity of the injection and mixing cases did not change significantly when the supply of culture medium was absent.

Mechanical properties and energy evolution of Beishan shallow-layer granite under different unloading paths

Rock has mechanical characteristics and a fracture damage mechanism that are closely related to its loading history and loading path. The mechanical properties, fracture damage features, acoustic emission(AE) characteristics, and strain energy evolution of the Beishan shallow-layer granite used in triaxial unloading tests were investigated in this study. Three groups of triaxial tests, namely, conventional triaxial compression test(Group Ⅰ), maintaining deviatoric stress synchronously unloading confining pressure test(Group Ⅱ), and loading axial pressure synchronously unloading confining pressure test(Group Ⅲ), were carried out for the cylindrical granite specimens. AE monitoring device was utilized in these tests to determine the degree to which the AE waves and AE events reflected the degree of rock damage. In addition, the crack stress thresholds of the specimens were determined by volumetric strain method and AE parameter method, and strain energy evolution of the rock was explored in different damage stages. The results show that the shallow-layer granite experiences brittle failure during the triaxial loading test and unloading test, and the rock has a greater damage degree during the unloading test. The crack stress thresholds of these samples vary greatly between tests, but the threshold ratios of all samples are similar in the same crack damage stage. The Mogi-Coulomb strength criterion can better describe the unloading failure strength of the rock. The evolution of the AE parameter characteristics and strain energy differs between the specimens used in different stress path tests. The dissipative strain energy is the largest in Group Ⅱ and the smallest in Group Ⅰ.

Mechanical behaviours of bedded sandstone under hydromechanical coupling

The combination of the dipping effect and hydromechanical(H-M)coupling effect can easily lead to water inrush disasters in water-rich roadways with different dip angles in coal mines.Therefore,H-M coupling tests of bedded sandstones under identical osmotic pressure and various confining pressures were conducted.Then,the evolution curves of stress-strain,permeability and damage,macro-and mesoscopic failure characteristics were obtained.Subsequently,the mechanical behaviour was characterized,and finally the failure mechanism was revealed.The results showed that:(1)The failure of the sandstone with the bedding angle of 45°or 60°was the structure-dominant type,while that with the bedding angle of 0°,30°or 90°was the force-dominant type.(2)When the bedding angle was in the range of(0°,30°)or(45°,90°),the confining pressure played a dominant role in influencing the peak strength.However,withinβ∈(30°,45°),the bedding effect played a dominant role in the peak strength.(3)With the increase in bedding angle,the cohesion increased first,then decreased and finally increased,while the internal friction angle was the opposite.(4)When the bedding angle was 0°or 30°,the“water wedging”effect and the“bedding buckling”effect would lead to the forking or converging shear failure.When the bedding angle was 45°or 60°,the sliding friction effect would lead to the shear slipping failure.When the bedding angle was 90°,the combination of the“bedding buckling”effect and shear effect would lead to the mixed tension-shear failure.The above conclusions obtained are helpful for the prevention of water inrush disasters in water-rich roadways with different dips in coal mines.

A phase-field model for simulating the propagation behavior of mixed-mode cracks during the hydraulic fracturing process in fractured reservoirs

A novel phase-field model for the propagation of mixed-mode hydraulic fractures,characterized by the formation of mixed-mode fractures due to the interactions between fluids and solids,is proposed.In this model,the driving force for the phase field consists of both tensile and shear components,with the fluid contribution primarily manifesting in the tension driving force.The displacement and pressure are solved simultaneously by an implicit method.The numerical solution's iterative format is established by the finite element discretization and Newton-Raphson(NR)iterative methods.The correctness of the model is verified through the uniaxial compression physical experiments on fluid-pressurized rocks,and the limitations of the hydraulic fracture expansion phase-field model,which only considers mode I fractures,are revealed.In addition,the influence of matrix mode II fracture toughness value,natural fracture mode II toughness value,and fracturing fluid injection rate on the hydraulic fracture propagation in porous media with natural fractures is studied.

Study on the low mechanical anisotropy of extruded Mg-Zn-Mn-Ce-Ca alloy tube in the compression process

In this study,the extruded Mg-Zn-Mn-Ce-Ca alloy tube with a low compression anisotropy along the ED,45ED and TD was prepared.The effect of the second phases,initial texture and deformation behavior on this low mechanical anisotropy was investigated.The results revealed that the alloy tube contains the high content(Mg1-xZnx)11Ce phase and the low content of Mg12Ce phase.These second phases are respectively incoherent and coherent with the Mg matrix,and their influence can be ignored.Additionally,the alloy tube exhibited a weak basal fiber texture,where the c-axis was aligned along the 0°∼30°tilt from TD to ED.Such a texture made the initial deformation(at 1.0%∼1.6%strain)of the three samples controlled by comparable basalslip.As deformation progressed(1.6∼9.0%strain),larger amounts of ETWs nucleated and gradually approached saturation in the three samples,re-orienting the c-axis to a 0°∼±30°deviation with respect to the loading directions.Meanwhile,the prismatic and pyramidal<c+a>slips replaced the dominant deformation progressively until fracture.Eventually,the similar deformation mechanisms determined by the weak initial texture in the three samples contribute to the comparable strain hardening rates,resulting in the low compressive anisotropy of the alloy tube.

Evaluating the stability and volumetric flowback rate of proppant packs in hydraulic fractures using the lattice Boltzmann-discrete element coupling method

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 towards an improved understanding of proppant flowback in hydraulic fractures,which also supports an optimised proppant size selection for hydraulic fracturing operations.