Effect of fracture fluid flowback on shale microfractures using CT scanning

The field data of shale fracturing demonstrate that the flowback performance of fracturing fluid is different from that of conventional reservoirs,where the flowback rate of shale fracturing fluid is lower than that of conventional reservoirs.At the early stage of flowback,there is no single-phase flow of the liquid phase in shale,but rather a gas-water two-phase flow,such that the single-phase flow model for tight oil and gas reservoirs is not applicable.In this study,pores and microfractures are extracted based on the experimental results of computed tomography(CT)scanning,and a spatial model of microfractures is established.Then,the influence of rough microfracture surfaces on the flow is corrected using the modified cubic law,which was modified by introducing the average deviation of the microfracture height as a roughness factor to consider the influence of microfracture surface roughness.The flow in the fracture network is simulated using the modified cubic law and the lattice Boltzmann method(LBM).The results obtained demonstrate that most of the fracturing fluid is retained in the shale microfractures,which explains the low fracturing fluid flowback rate in shale hydraulic fracturing.

Influences of clean fracturing fluid viscosity and horizontal in-situ stress difference on hydraulic fracture propagation and morphology in coal seam

The viscosity of fracturing fluid and in-situ stress difference are the two important factors that affect the hydraulic fracturing pressure and propagation morphology. In this study, raw coal was used to prepare coal samples for experiments, and clean fracturing fluid samples were prepared using CTAB surfactant. A series of hydraulic fracturing tests were conducted with an in-house developed triaxial hydraulic fracturing simulator and the fracturing process was monitored with an acoustic emission instrument to analyze the influences of fracturing fluid viscosity and horizontal in-situ stress difference on coal fracture propagation. The results show that the number of branched fractures decreased, the fracture pattern became simpler, the fractures width increased obviously, and the distribution of AE event points was concentrated with the increase of the fracturing fluid viscosity or the horizontal in-situ stress difference. The acoustic emission energy decreases with the increase of fracturing fluid viscosity and increases with the increase of horizontal in situ stress difference. The low viscosity clean fracturing fluid has strong elasticity and is easy to be compressed into the tip of fractures, resulting in complex fractures. The high viscosity clean fracturing fluids are the opposite. Our experimental results provide a reference and scientific basis for the design and optimization of field hydraulic fracturing parameters.

Fluid identification based on frequency-dependent AVO attribute inversion in multi-scale fracture media

A key problem in seismic inversion is the identification of the reservoir fluids. Elastic parameters, such as seismic wave velocity and formation density, do not have sufficient sensitivity, thus, the conventional amplitude-versus-offset（AVO） method is not applicable. The frequency-dependent AVO method considers the dependency of the seismic amplitude to frequency and uses this dependency to obtain information regarding the fluids in the reservoir fractures. We propose an improved Bayesian inversion method based on the parameterization of the Chapman model. The proposed method is based on 1） inelastic attribute inversion by the FDAVO method and 2） Bayesian statistics for fluid identification. First, we invert the inelastic fracture parameters by formulating an error function, which is used to match observations and model data. Second, we identify fluid types by using a Markov random field a priori model considering data from various sources, such as prestack inversion and well logs. We consider the inelastic parameters to take advantage of the viscosity differences among the different fluids possible. Finally, we use the maximum posteriori probability for obtaining the best lithology/fluid identification results.

Microscopic characteristics of tight sandstone reservoirs and their effects on the imbibition efficiency of fracturing fluids:A case study of the Linxing area,Ordos Basin

The Linxing area within the Ordos Basin exhibits pronounced reservoir heterogeneity and intricate micro-pore structures,rendering it susceptible to water-blocking damage during imbibition extraction.This study delved into the traits of tight sandstone reservoirs in the 8th member of the Shihezi Formation(also referred to as the He 8 Member)in the study area,as well as their effects on fracturing fluid imbibition.Utilizing experimental techniques such as nuclear magnetic resonance(NMR),high-pressure mercury intrusion(HPMI),and gas adsorption,this study elucidated the reservoir characteristics and examined the factors affecting the imbibition through imbibition experiments.The findings reveal that:①The reservoir,with average porosity of 8.40%and average permeability of 0.642×10^(-3)μm^(2),consists principally of quartz,feldspar,and lithic fragments,with feldspathic litharenite serving as the primary rock type and illite as the chief clay mineral;②Nano-scale micro-pores and throats dominate the reservoir,with dissolution pores and intercrystalline pores serving as predominant pore types,exhibiting relatively high pore connectivity;③Imbibition efficiency is influenced by petrophysical properties,clay mineral content,and microscopic pore structure.Due to the heterogeneity of the tight sandstone reservoir,microscopic factors have a more significant impact on the imbibition efficiency of fracturing fluids;④A comparative analysis shows that average pore size correlates most strongly with imbibition efficiency,followed by petrophysical properties and clay mineral content.In contrast,the pore type has minimal impact.Micropores are vital in the imbibition process,while meso-pores and macro-pores offer primary spaces for imbibition.This study offers theoretical insights and guidance for enhancing the post-fracturing production of tight sandstone reservoirs by examining the effects of these factors on the imbibition efficiency of fracturing fluids in tight sandstones.

A laboratory acoustic emission experiment and numerical simulation of rock fracture driven by a high-pressure fluid source

In order to improve our understanding of rock fracture and fault instability driven by high-pressure fluid sources, the authors carried out rock fracture tests using granite under a confining pressure of 80 MPa with fluid injection in the laboratory. Furthermore, we tested a number of numerical models using the FLAC;modeling software to find the best model to represent the experimental results. The high-speed multichannel acoustic emission（AE） waveform recording system used in this study made it possible to examine the total fracture process through detailed monitoring of AE hypocenters and seismic velocity.The experimental results show that injecting high-pressure oil into the rock sample can induce AE activity at very low stress levels and can dramatically reduce the strength of the rock. The results of the numerical simulations show that major experimental results, including the strength, the temporal and spatial patterns of the AE events, and the role of the fluid can be represented fairly well by a model involving（1） randomly distributed defect elements to model pre-existing cracks,（2） random modification of rock properties to represent inhomogeneity introduced by different mineral grains, and（3）macroscopic inhomogeneity. Our study, which incorporates laboratory experiments and numerical simulations, indicates that such an approach is helpful in finding a better model not only for simulating experimental results but also for upscaling purposes.

LBM simulation of fluid flow in fractured porous media with permeable matrix

To analyze and depict complicated fluid behaviors in fractured porous media with variably permeable matrix,an integrated discrete computational algorithm is proposed based on lattice Boltzmann method（LBM）.This paper combines with the external force model and statistical material physics to effectively describe the feature changes while the fluid passes through the fractures within the permeable matrix.As an application example,a two dimensional rock sample is reconstructed using the digital image and characterized with different feature values at each LBM grid to distinguish pores,impermeable and permeable matrix by stating its local physical property.Compared with the conventional LBM,the results demonstrate the advantages of proposed algorithm in modeling fluid flow phenomenon in fractured porous media with variably permeable matrix.

Experimental study of reservoir damage of water-based fracturing fluids prepared by different polymers

Fracturing operations can effectively improve the production of low-permeable reservoirs. The performance of fracturing fluids directly affects the fracturing efficiency and back flow capacity. As polymerbased fracturing fluids(such as guar gum(GG), polyacrylamide(HPAM), etc.) are high-viscosity fluids formed by viscosifiers and crosslinking agents, the degree of gel breakage after the fracturing operation directly influences the damage degree to the reservoir matrix and the mobility of oil angd gas produced from the reservoir into the wellbore. This study compared the viscosity, molecular weight, and particle size of the fracturing fluid after gel breakage prepared by GG and HPAM as viscosifiers, as well as evaluate their damage to the core. Results show that the viscosities of the gel-breaking fluid increased with the concentration of the viscosifier for both the HPAM-based and GG-based fracturing fluids. For the breaking fluid with the same viscosity, the molecular weight in the HPAM-based gel-breaking fluid was much larger than that in the GG-based system. Moreover, for the gel-breaking fluid with the same viscosity, the molecular particle size of the residual polymers in the HPAM-based system was smaller than that in the GG-based system. The damage to the core with the permeability of 1 × 10^(-3)μm^(2) caused by both the HPAM-based and GG-based gel-breaking fluids decreased with the increase in the solution viscosity. For the gel-breaking fluid systems with the same viscosity(i.e., 2-4 mPa s), the damage of HPAM-based fracturing fluid to low-permeability cores was greater than the GG-based fracturing fluid(45.6%-80.2%) since it had a smaller molecular particle size, ranging from 66.2% to 77.0%. This paper proposed that the damage caused by hydraulic fracturing in rock cores was related to the partilce size of residual polymers in gel-breaking solution, rather than its molecular weight. It was helpful for screening and optimizing viscosifiers used in hydraulic fracturing process.

Probing the influence of secondary fracture connectivity on fracturing fluid flowback efficiency

A deep understanding of the geometric impacts of fracture on fracturing fluid flowback efficiency is essential for unconventional oil development. Using nuclear magnetic resonance and 2.5-dimensional matrix-fracture visualization microfluidic models, qualitative and quantitative descriptions of the influences of connectivity between primary fracture and secondary fracture on flowback were given from core scale to pore network scale. The flow patterns of oil-gel breaking fluid two-phase flow during flowback under different fracture connectivity were analyzed. We found some counterintuitive results that non-connected secondary fracture (NCSF, not connect with artificial primary fracture and embedded in the matrix) is detrimental to flowbackefficiency. The NCSF accelerates the formation of oil channeling during flowback, resulting in a large amount of fracturing fluid trapped in the matrix, which is not beneficial for flowback. Whereas the connected secondary fracture (CSF, connected with the artificial primary fracture) is conducive to flowback. The walls of CSF become part of primary fracture, which expands the drainage area with low resistance, and delays the formation of the oil flow channel. Thus, CSF increases the high-speed flowback stage duration, thereby enhancing the flowback efficiency. The fracturing fluid flowback efficiency investigated here follows the sequence of the connected secondary fracture model (72%) > the matrix model (66%) > the non-connected secondary fracture model (38%). Our results contribute to hydraulic fracturing design and the prediction of flowback efficiency.

A novel dandelion-based bionic proppant and its transportation mechanism in different types of fractures

Low-permeability reservoirs are generally characterized by low porosity and low permeability.Obtaining high production using the traditional method is technologically challenging because it yields a low reservoir recovery factor.In recent years,hydraulic fracturing technology is widely applied for efficiently exploiting and developing low-permeability reservoirs using a low-viscosity fluid as a fracturing fluid.However,the transportation of the proppant is inefficient in the low-viscosity fluid,and the proppant has a low piling-up height in fracture channels.These key challenges restrict the fluid(natural gas or oil)flow in fracture channels and their functional flow areas,reducing the profits of hydrocarbon exploitation.This study aimed to explore and develop a novel dandelion-bionic proppant by modifying the surface of the proppant and the fiber.Its structure was similar to that of dandelion seeds,and it had high transport and stacking efficiency in low-viscosity liquids compared with the traditional proppant.Moreover,the transportation efficiency of this newly developed proppant was investigated experimentally using six different types of fracture models(tortuous fracture model,rough fracture model,narrow fracture model,complex fracture model,large-scale single fracture model,and small-scale single fracture model).Experimental results indicated that,compared with the traditional proppant,the transportation efficiency and the packing area of the dandelion-based bionic proppant significantly improved in tap water or low-viscosity fluid.Compared with the traditional proppant,the dandelionbased bionic proppant had 0.1-4 times longer transportation length,0.3-5 times higher piling-up height,and 2-10 times larger placement area.The newly developed proppant also had some other extraordinary features.The tortuosity of the fracture did not influence the transportation of the novel proppant.This proppant could easily enter the branch fracture and narrow fracture with a high packing area in rough surface fractures.Based on the aforementioned characteristics,this novel proppant technique could improve the proppant transportation efficiency in the low-viscosity fracturing fluid and increase the ability of the proppant to enter the secondary fracture.This study might provide a new solution for effectively exploiting low-permeability hydrocarbon reservoirs.

Optimization method of fracturing fluid volume intensity for SRV fracturing technique in shale oil reservoir based on forced imbibition:A case study of well X-1 in Biyang Sag of Nanxiang Basin,China

An optimization method of fracturing fluid volume strength was introduced taking well X-1 in Biyang Sag of Nanxiang Basin as an example.The characteristic curves of capillary pressure and relative permeability were obtained from history matching between forced imbibition experimental data and core-scale reservoir simulation results and taken into a large scale reservoir model to mimic the forced imbibition behavior during the well shut-in period after fracturing.The optimization of the stimulated reservoir volume(SRV)fracturing fluid volume strength should meet the requirements of estimated ultimate recovery(EUR),increased oil recovery by forced imbibition and enhancement of formation pressure and the fluid volume strength of fracturing fluid should be controlled around a critical value to avoid either insufficiency of imbibition displacement caused by insufficient fluid amount or increase of costs and potential formation damage caused by excessive fluid amount.Reservoir simulation results showed that SRV fracturing fluid volume strength positively correlated with single-well EUR and an optimal fluid volume strength existed,above which the single-well EUR increase rate kept decreasing.An optimized increase of SRV fracturing fluid volume and shut-in time would effectively increase the formation pressure and enhance well production.Field test results of well X-1 proved the practicality of established optimization method of SRV fracturing fluid volume strength on significant enhancement of shale oil well production.

Formation damage mechanism and control strategy of the compound function of drilling fluid and fracturing fluid in shale reservoirs

For the analysis of the formation damage caused by the compound function of drilling fluid and fracturing fluid,the prediction method for dynamic invasion depth of drilling fluid is developed considering the fracture extension due to shale minerals erosion by oil-based drilling fluid.With the evaluation for the damage of natural and hydraulic fractures caused by mechanical properties weakening of shale fracture surface,fracture closure and rock powder blocking,the formation damage pattern is proposed with consideration of the compound effect of drilling fluid and fracturing fluid.The formation damage mechanism during drilling and completion process in shale reservoir is revealed,and the protection measures are raised.The drilling fluid can deeply invade into the shale formation through natural and induced fractures,erode shale minerals and weaken the mechanical properties of shale during the drilling process.In the process of hydraulic fracturing,the compound effect of drilling fluid and fracturing fluid further weakens the mechanical properties of shale,results in fracture closure and rock powder shedding,and thus induces stress-sensitive damage and solid blocking damage of natural/hydraulic fractures.The damage can yield significant conductivity decrease of fractures,and restrict the high and stable production of shale oil and gas wells.The measures of anti-collapse and anti-blocking to accelerate the drilling of reservoir section,forming chemical membrane to prevent the weakening of the mechanical properties of shale fracture surface,strengthening the plugging of shale fracture and reducing the invasion range of drilling fluid,optimizing fracturing fluid system to protect fracture conductivity are put forward for reservoir protection.

Fluid identification and effective fracture prediction based on frequency-dependent AVOAz inversion for fractured reservoirs

Fluid and effective fracture identification in reservoirs is a crucial part of reservoir prediction.The frequency-dependent AVO inversion algorithms have proven to be effective for identifying fluid through its dispersion property.However,the conventional frequency-dependent AVO inversion algorithms based on Smith&Gidlow and Aki&Richards approximations do not consider the acquisition azimuth of seismic data and neglect the effect of seismic anisotropic dispersion in the actual medium.The aligned fractures in the subsurface medium induce anisotropy.The seismic anisotropy should be considered while accounting for the seismic dispersion properties through fluid-saturated fractured reservoirs.Anisotropy in such reservoirs is frequency-related due to wave-induced fluid-flow(WIFF)between interconnected fractures and pores.It can be used to identify fluid and effective fractures(fluid-saturated)by using azimuthal seismic data via anisotropic dispersion properties.In this paper,based on Rüger’s equation,we derived an analytical expression in the frequency domain for the frequencydependent AVOAz inversion in terms of fracture orientation,dispersion gradient of isotropic background rock,anisotropic dispersion gradient,and the dispersion at a normal incident angle.The frequency-dependent AVOAz equation utilizes azimuthal seismic data and considers the effect of both isotropic and anisotropic dispersion.Reassigned Gabor Transform(RGT)is used to achieve highresolution frequency division data.We then propose the frequency-dependent AVOAz inversion method to identify fluid and characterize effective fractures in fractured porous reservoirs.Through application to high-qualified seismic data of dolomite and carbonate reservoirs,the results show that the method is useful for identifying fluid and effective fractures in fluid-saturated fractured rocks.

Supramolecular polymer-based gel fracturing fluid with a double network applied in ultra-deep hydraulic fracturing

A gel based on polyacrylamide,exhibiting delayed crosslinking characteristics,emerges as the preferred solution for mitigating degradation under conditions of high temperature and extended shear in ultralong wellbores.High viscosity/viscoelasticity of the fracturing fluid was required to maintain excellent proppant suspension properties before gelling.Taking into account both the cost and the potential damage to reservoirs,polymers with lower concentrations and molecular weights are generally preferred.In this work,the supramolecular action was integrated into the polymer,resulting in significant increases in the viscosity and viscoelasticity of the synthesized supramolecular polymer system.The double network gel,which is formed by the combination of the supramolecular polymer system and a small quantity of Zr-crosslinker,effectively resists temperature while minimizing permeability damage to the reservoir.The results indicate that the supramolecular polymer system with a molecular weight of(268—380)×10^(4)g/mol can achieve the same viscosity and viscoelasticity at 0.4 wt%due to the supramolecular interaction between polymers,compared to the 0.6 wt%traditional polymer(hydrolyzed polyacrylamide,molecular weight of 1078×10^(4)g/mol).The supramolecular polymer system possessed excellent proppant suspension properties with a 0.55 cm/min sedimentation rate at 0.4 wt%,whereas the0.6 wt%traditional polymer had a rate of 0.57 cm/min.In comparison to the traditional gel with a Zrcrosslinker concentration of 0.6 wt%and an elastic modulus of 7.77 Pa,the double network gel with a higher elastic modulus(9.00 Pa)could be formed only at 0.1 wt%Zr-crosslinker,which greatly reduced the amount of residue of the fluid after gel-breaking.The viscosity of the double network gel was66 m Pa s after 2 h shearing,whereas the traditional gel only reached 27 m Pa s.

Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

A Novel Fracturing Fluid with High-Temperature Resistance for Ultra-Deep Reservoirs

Ultra-deep reservoirs play an important role at present in fossil energy exploitation.Due to the related high temperature,high pressure,and high formation fracture pressure,however,methods for oil well stimulation do not produce satisfactory results when conventional fracturing fluids with a low pumping rate are used.In response to the above problem,a fracturing fluid with a density of 1.2~1.4 g/cm^(3)was developed by using Potassium formatted,hydroxypropyl guanidine gum and zirconium crosslinking agents.The fracturing fluid was tested and its ability to maintain a viscosity of 100 mPa.s over more than 60 min was verified under a shear rate of 1701/s and at a temperature of 175℃.This fluid has good sand-carrying performances,a low viscosity after breaking the rubber,and the residue content is less than 200 mg/L.Compared with ordinary reconstruction fluid,it can increase the density by 30%~40%and reduce the wellhead pressure of 8000 m level reconstruction wells.Moreover,the new fracturing fluid can significantly mitigate safety risks.

Investigation on microscopic invasion characteristics and retention mechanism of fracturing fluid in fractured porous media

Reservoir damage caused by guar gum fracturing fluid and slick water seriously affects the subsequent oil and gas production. However, the invasion characteristics and retention mechanisms of fracturing fluids in the fracture-matrix zone are still unclear. In this work, a microscopic model reflecting the characteristics of the fracture-matrix zone was designed. Based on the microfluidic experimental method, the process of fracturing fluid invasion, flowback and retention in the fracture-matrix zone was investigated visually and characterized quantitatively. The factors and mechanisms affecting fracturing fluid retention in the fracture-matrix zone were analyzed and clarified. The results indicated that in the invasion process, the frontal swept range of slick water was larger than that of the guar gum fracturing fluid, and the oil displacement efficiency and damage rate were lower than those of the guar gum fracturing fluid under the same invasion pressure. With the increase in invasion pressure, the damage rate of slick water increased from 61.09% to 82.77%, and that of the guar gum fracturing fluid decreased from 93.45% to83.36%. Before subsequent oil production, the invaded fracturing fluid was mainly concentrated in the medium-high permeability area of the fracture-matrix zone. The main resistance of slick water was capillary force, while that of the guar fracturing fluid was mainly viscous resistance. The fracturing fluid retention was most serious in the low permeability region and the region near the end of the fracture.The experimental and numerical simulation results showed that increasing the production pressure difference could improve the velocity field distribution of the fracture-matrix zone, increase the flowback swept range and finally reduce the retention rate of the fracture fluid. The retention mechanisms of slick water in the fracture-matrix zone include emulsion retention and flow field retention, while those of the guar gum fracturing fluid include viscous retention and flow field retention. Emulsion retention is caused by capillary force and flow interception effect. Viscous retention is caused by the viscous resistance of polymer, while flow-field retention is caused by uneven distribution of flowback velocity.

Numerical analysis of seismic wave propagation in fluid-saturated porous multifractured media

Elastic wave propagation and attenuation in porous rock layers with oriented sets of fractures, especially in carbonate reservoirs, are anisotropic owing to fracture sealing, fracture size, fracture density, filling fluid, and fracture strike orientation. To address this problem, we adopt the Chapman effective medium model and carry out numerical experiments to assess the variation in P-wave velocity and attenuation, and the shear-wave splitting anisotropy with the frequency and azimuth of the incident wave. The results suggest that velocity, attenuation, and anisotropy vary as function of azimuth and frequency. The azimuths of the minimum attenuation and maximum P-wave velocity are nearly coincident with the average strike of the two sets of open fractures. P-wave velocity is greater in sealed fractures than open fractures, whereas the attenuation of energy and anisotropy is stronger in open fractures than sealed fractures. For fractures of different sizes, the maximum velocity together with the minimum attenuation correspond to the average orientation of the fracture sets. Small fractures affect the wave propagation less. Azimuth-dependent anisotropy is low and varies more than the other attributes. Fracture density strongly affects the P-wave velocity, attenuation, and shear-wave anisotropy. The attenuation is more sensitive to the variation of fracture size than that of velocity and anisotropy. In the seismic frequency band, the effect of oil and gas saturation on attenuation is very different from that for brine saturation and varies weakly over azimuth. It is demonstrated that for two sets of fractures with the same density, the fast shear-wave polarization angle is almost linearly related with the orientation of one of the fracture sets.

Nonlinear fluid flow through three-dimensional rough fracture networks:Insights from 3D-printing,CT-scanning,and high-resolution numerical simulations

Nonlinear flow behavior of fluids through three-dimensional(3D)discrete fracture networks(DFNs)considering effects of fracture number,surface roughness and fracture aperture was experimentally and numerically investigated.Three physical models of DFNs were 3D-printed and then computed tomography(CT)-scanned to obtain the specific geometry of fractures.The validity of numerically simulating the fluid flow through DFNs was verified via comparison with flow tests on the 3D-printed models.A parametric study was then implemented to establish quantitative relations between the coefficients/parameters in Forchheimer’s law and geometrical parameters.The results showed that the 3D-printing technique can well reproduce the geometry of single fractures with less precision when preparing complex fracture networks,numerical modeling precision of which can be improved via CT-scanning as evidenced by the well fitted results between fluid flow tests and numerical simulations using CT-scanned digital models.Streamlines in DFNs become increasingly tortuous as the fracture number and roughness increase,resulting in stronger inertial effects and greater curvatures of hydraulic pressure-low rate relations,which can be well characterized by the Forchheimer’s law.The critical hydraulic gradient for the onset of nonlinear flow decreases with the increasing aperture,fracture number and roughness,following a power function.The increases in fracture aperture and number provide more paths for fluid flow,increasing both the viscous and inertial permeabilities.The value of the inertial permeability is approximately four orders of magnitude greater than the viscous permeability,following a power function with an exponent a of 3,and a proportional coefficient b mathematically correlated with the geometrical parameters.

Research on the Performance of New Weighted Slippery Water Fracturing Fluid System

Deep and ultra-deep reservoirs have dense matrix and high fracture pressure, which leads to high pressure and difficulty in fracturing construction. Conventional aggravated fracturing fluids have the problems of low aggravation efficiency, high friction resistance, etc., and the reduction of construction pressure cannot reach the theoretical effect. In view of the above problems, this paper adopts the weighting agent HD160 and the drag reducing agent JHFR-2 to form a new type of weighted slippery water fracturing fluid system. And the weighting performance, drag reduction performance, corrosion performance, anti-expansion performance and reservoir damage of this system were studied. The results show that the density of the system is adjustable within 1.1 - 1.6 g·cm−3, and the drag reduction rate can be up to 68% at 1.5 g·cm−3, with low corrosion rate, surface tension less than 28 mN·m−1, anti-expansion rate as high as 94.5%, and the damage rate of the reservoir permeability is less than 10%, which is of good application prospect.

Application of Stochastic Fracture Network with Numerical Fluid Flow Simulations to Groundwater Flow Modeling in Fractured Rocks

The continuum approach in fluid flow modeling is generally applied to porous geological media, but has limited applicability to fractured rocks. With the presence of a discrete fracture network relatively sparsely distributed in the matrix, it may be difficult or erroneous to use a porous medium fluid flow model with continuum assumptions to describe the fluid flow in fractured rocks at small or even large field scales. A discrete fracture fluid flow approach incorporating a stochastic fracture network with numerical fluid flow simulations could have the capability of capturing fluid flow behaviors such as inhomogeneity and anisotropy while reflecting the changes of hydraulic features at different scales. Moreover, this approach can be implemented to estimate the size of the representative elementary volume (REV) in order to find out the scales at which a porous medium flow model could be applied, and then to determine the hydraulic conductivity tensor for fractured rocks. The following topics are focused on in this study: (a) conceptual discrete fracture fluid flow modeling incorporating a stochastic fracture network with numerical flow simulations; (b) estimation of REV and hydraulic conductivity tensor for fractured rocks utilizing a stochastic fracture network with numerical fluid flow simulations; (c) investigation of the effect of fracture orientation and density on the hydraulic conductivity and REV by implementing a stochastic fracture network with numerical fluid flow simulations, and (d) fluid flow conceptual models accounting for major and minor fractures in the 2 D or 3 D flow fields incorporating a stochastic fracture network with numerical fluid flow simulations.