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.

Sensitivity of seismic attenuation and dispersion to dynamic elastic interactions of connected fractures: Quasi-static finite element modeling study

Prediction of seismic attenuation and dispersion that are inherently sensitive to hydraulic and elastic properties of the medium of interest in the presence of mesoscopic fractures and pores,is of great interest in the characterization of fractured formations.This has been very difficult,however,considering that stress interactions between fractures and pores,related to their spatial distributions,tend to play a crucial role on affecting overall dynamic elastic properties that are largely unexplored.We thus choose to quantitatively investigate frequency-dependent P-wave characteristics in fractured porous rocks at the scale of a representative sample using a numerical scale-up procedure via performing finite element modelling.Based on 2-D numerical quasi-static experiments,effects of fracture and fluid properties on energy dissipation in response to wave-induced fluid flow at the mesoscopic scale are quantified via solving Biot's equations of consolidation.We show that numerical results are sensitive to some key characteristics of probed synthetic rocks containing unconnected and connected fractures,demonstrating that connectivity,aperture and inclination of fractures as well as fracture infills exhibit strong impacts on the two manifestations of WIFF mechanisms in the connected scenario,and on resulting total wave attenuation and phase velocity.This,in turn,illustrates the importance of these two WIFF mechanisms in fractured rocks and thus,a deeper understanding of them may eventually allow for a better characterization of fracture systems using seismic methods.Moreover,this presented work combines rock physics predictions with seismic numerical simulations in frequency domain to illustrate the sensitivity of seismic signatures on the monitoring of an idealized geologic CO_(2) sequestration in fractured reservoirs.The simulation demonstrates that these two WIFF mechanisms can strongly modify seismic records and hence,indicating that incorporating the two energy dissipation mechanisms in the geophysical interpretation can potentially improving the monitoring and surveying of fluid variations in fractured formations.

Numerical simulation of two-phase flow in fractured porous media using streamline simulation and IMPES methods and comparing results with a commercial software

Streamline simulation is developed to simulate waterflooding in fractured reservoirs. Conventional reservoir simulation methods for fluid flow simulation in large and complex reservoirs are very costly and time consuming. In streamline method, transport equations are solved on one-dimensional streamlines to reduce the computation time with less memory for simulation. First, pressure equation is solved on an Eulerian grid and streamlines are traced. Defining the "time of flight", saturation equations are mapped and solved on streamlines. Finally, the results are mapped back on Eulerian grid and the process is repeated until the simulation end time. The waterflooding process is considered in a fractured reservoir using the dual porosity model. Afterwards, a computational code is developed to solve the same problem by the IMPES method and the results of streamline simulation are compared to those of the IMPES and a commercial software. Finally, the accuracy and efficiency of streamline simulator for simulation of two-phase flow in fractured reservoirs has been proved.

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.

Insight into the structural evolution of porous and fractured media by forced aeration during heap leaching

Despite lots of techniques in improving the heap leaching performance,many constraints on the industrial applications remain.We proposed a correspondingly effective and new idea of introducing forced aeration to improve the bad permeability and leaching effect of Yangla Copper Mine(YCM)during heap leaching.The dual-media theory was employed to study the impact mechanism of forced aeration on the variations of porous and fractured media during the column leaching experiments.An X-Ray Computed Tomography(CT)set was utilized to perform the pore imaging of the specimens and the fracture morphology of the particles within the columns was analyzed by Scanning Electron Microscope(SEM)as aeration rate(AR)changed.The results show that there exists copious fine particles within the heap of YCM,the particle size distribution of which is not reasonable.The forced aeration can not only promote the development of the porous and fractured structures but effectively break the blocked seepage paths.Then the leaching degree is improved and the seepage performance of the solute within the solution is enhanced.Therefore,the forced aeration is probable of making the leaching performance greatly improved.

P-wave attenuation anisotropy in TI media and its application in fracture parameters inversion

The existence of aligned fractures in fluid-saturated rocks leads to obvious attenuation anisotropy and velocity anisotropy. Attenuation anisotropy analysis can be applied to estimate fracture density and scale, which provide important information for reservoir identification. This paper derives P-wave attenuation anisotropy in the ATI media where the symmetry axis is in the arbitrary direction theoretically and modifies the spectral ratio method to measure attenuation anisotropy in the ATI media, thus avoiding a large measurement error when applied to wide azimuth or full azimuth data. Fracture dip and azimuth can be estimated through attenuation anisotropy analysis. For small-scale fractures, fracture scale and fracture density can be determined with enhanced convergence if velocity and attenuation information are both used. We also apply the modified spectralratio method to microseismic field data from an oilfield in East China and extract the fracture dip through attenuation anisotropy analysis. The result agrees with the microseismie monitoring.

Review of permeability evolution model for fractured porous media

The ability to capture permeability of fractured porous media plays a significant role in several engineering applications, including reservoir, mining, petroleum and geotechnical engineering. In order to solve fluid flow and coupled flow-deformation problems encountered in these engineering applications,both empirical and theoretical models had been proposed in the past few decades. Some of them are simple but still work in certain circumstances; others are complex but also need some modifications to be applicable. Thus, the understanding of state-of-the-art permeability evolution model would help researchers and engineers solve engineering problems through an appropriate approach. This paper summarizes permeability evolution models proposed by earlier and recent researchers with emphasis on their characteristics and limitations.

Experimental investigation on rainfall infiltration and solute transport in layered porous and fractured media

Layered structures with upper porous and lower fractured media are widely distributed in the world. An experimen- tal investigation on rainfall infiltration and solute transport in such layered structures can provide the necessary foundation for effectively preventing and forecasting water bursting in mines, controlling contamination of mine water, and accomplishing ecological restoration of mining areas. A typical physical model of the layered structures with porous and fractured media was created in this study. Then rainfall infiltration experiments were conducted after salt solution was sprayed on the surface of the layered structure. The volumetric water content and concentration of chlorine ions at different specified positions along the profile of the experiment system were measured in real-time. The experimental results showed that the lower fractured media, with a considerably higher permeability than that of the upper porous media, had significant effects on preventing water infil- tration. Moreover, although the porous media were homogeneous statistically in the whole domain, spatial variations in the features of effluent concentrations with regards to time, or so called breakthrough curves, at various sampling points located at the horizontal plane in the porous media near the porous-fractured interface were observed, indicating the diversity of solute transport at small scales. Furthermore, the breakthrough curves of the outflow at the bottom, located beneath the underlying fractured rock, were able to capture and integrate features of the breakthrough curves of both the upper porous and fractured media, which exhibited multiple peaks, while the peak values were reduced one by one with time.

Rock matrix-fractured media model for heterogeneous and fractured coal bed

Coal was considered rock matrix-fractured media composed of rock matrix and fractures, and the rock matrix-fractured media model for heterogeneous and fractured coal bed was presented. In this model the rock matrix is heterogeneous, and the mechanical parameters such as elastic modulus and strength follow Weibull distribution. Fractures in coal bed were generated with the discrete fracture network method, and the properties of fractures were simulated with Desai element. Then the virtual generating system (VGS) of natural heterogeneous and fractured coal bed was developed in Matlab 6.0. The coupled model of gas flow and deformation process based on the rock matrix-fractured media model method and VGS for heterogeneous and fractured coal bed was presented, and the numerical code was developed in Matlab 6.0. The gas flow process in the heterogeneous and fractured coal bed was simulated in a numerical case. The main conclusions are: 1) The natural heterogeneous and fractured coal bed could be simulated by the rock matrix-fractured media model and VGS; 2) The fractures connected with the well have much more effects on gas flow than those non-connected.

Hydrothermal fluid circulation in anisotropic permeable media associated with discrete fractures

A proper form of the Rayleigh number, containing the geometric mean of the vertical and horizontal permeabilities was obtained. The critical value for the onset of stable convection was found. The results proved analytically and numerically that anisotropy in permeability resists the initiation of hydrothermal convection. The equivalence between homogeneously anisotropic media and multiply fractured media was also investigated. It was confirmed that multiply fractured models are comparable to anisotropic models as long as they have the same averaged horizontal or vertical permeabilities and other physical parameters.

Fracture prediction approach for oil-bearing reservoirs based on AVAZ attributes in an orthorhombic medium

Fracture systems in nature are complicated. Normally vertical fractures develop in an isotropic background. However, the presence of horizontal fine layering or horizontal fractures in reservoirs makes the vertical fractures develop in a VTI（a transversely isotropic media with a vertical symmetry axis） background. In this case, reservoirs can be described better by using an orthorhombic medium instead of a traditional HTI（a transversely isotropic media with a horizontal symmetry axis） medium. In this paper, we focus on the fracture prediction study within an orthorhombic medium for oil-bearing reservoirs. Firstly, we simplify the reflection coefficient approximation in an orthorhombic medium. Secondly, the impact of horizontal fracturing on the reflection coefficient approximation is analyzed theoretically. Then based on that approximation, we compare and analyze the relative impact of vertical fracturing, horizontal fracturing and fluid indicative factor on traditional ellipse fitting results and the scaled B attributes. We find that scaled B attributes are more sensitive to vertical fractures, so scaled B attributes are proposed to predict vertical fractures. Finally, a test is developed to predict the fracture development intensity of an oil-bearing reservoir. The fracture development observed in cores is used to validate the study method. The findings of both theoretical analyses and practical application reveal that compared with traditional methods, this new approach has improved the prediction of fracture development intensity in oil-bearing reservoirs.

Permeability Heterogeneity in a Fractured Sandstone-Mudstone Rock Mass in Xiaolangdi Dam Site,Central China

Heterogeneity of permeability in fractured media is a hot research topic in hydrogeology. Numerous approaches have been proposed to characterize heterogeneity in the last several decades. However, little attention has been paid to correlate permeability heterogeneity with geological information. In the present study, several causes of permeability heterogeneity, that is, lithology, tectonism, and depth, are identified. The unit absorption values （denoted as ω）, which are results obtained from the packer test, are employed to represent permeability. The variability of permeability in sandstone-mudstone is so significant that the value of unit absorptions span 3-4 orders of magnitude at any depth with several test sections. By declustering, it has been found that under a similar tectonic history, the means of permeability differ greatly at different formations as a result of different mudrock contents. It has also been found that in the same formation, permeability can be significantly increased as a result of faulting. The well-known phenomenon, the decrease in permeability with depth, is found to be caused by the fractures in the rock mass, and the relationship between permeability and depth can be established in the form of logoω-logd. After subtracting the trend of ω with absolute depth, the mean of the residual value at each relative depth can be well correlated with the distribution of mudstone. The methods proposed in this paper can be utilized to research in similar study areas.

Estimation of fracture density and orientation from azimuthal elastic impedance difference through singular value decomposition

Accurate estimation of fracture density and orientation is of great significance for seismic characterization of fractured reservoirs.Here,we propose a novel methodology to estimate fracture density and orientation from azimuthal elastic impedance(AEI)difference using singular value decomposition(SVD).Based on Hudson's model,we first derive the AEI equation containing fracture density in HTI media,and then obtain basis functions and singular values from the normalized AEI difference utilizing SVD.Analysis shows that the basis function changing with azimuth is related to fracture orientation,fracture density is the linearly weighted sum of singular values,and the first singular value contributes the most to fracture density.Thus,we develop an SVD-based fracture density and orientation inversion approach constrained by smooth prior elastic parameters.Synthetic example shows that fracture density and orientation can be stably estimated,and the correlation coefficient between the true value and the estimated fracture density is above 0.85 even when an S/N ratio of 2.Field data example shows that the estimated fracture orientation is consistent with the interpretation of image log data,and the estimated fracture density reliably indicates fractured gas-bearing reservoir,which could help to guide the exploration and development of fractured reservoirs.

DYNAMIC PRODUCTION PREDICTION AND PARAMETER IDENTIFICATION FOR GAS WELL WITH VERTICAL FRACTURE

In order to devoid the hard work and factitious error in selecting charts while analyzing and interpreting hydraulic fracturing fracture parameters, on the basis of the non-Darcy flow factor, this paper put out the non-Darcy flow mathematical model of real gas in the formation and fracture, established the production history automatic matching model to identify fracture parameters, and offered the numerical solutions of those models, which took the variation of fracture conductivity in production process. These results offered a precise and reliable method to understand formation, analyze and evaluate the fracturing treatment quality of gas well.

Modeling the efects of fracture infll on frequency‑dependent anisotropy and AVO response of a fractured porous layer

In a fractured porous hydrocarbon reservoir,wave velocities and refections depend on frequency and incident angle.A proper description of the frequency dependence of amplitude variations with ofset(AVO)signatures should allow efects of fracture inflls and attenuation and dispersion of fractured media.The novelty of this study lies in the introduction of an improved approach for the investigation of incident-angle and frequency variations-associated refection responses.The improved AVO modeling method,using a frequency-domain propagator matrix method,is feasible to accurately consider velocity dispersion predicted from frequency-dependent elasticities from a rock physics modeling.And hence,the method is suitable for use in the case of an anisotropic medium with aligned fractures.Additionally,the proposed modeling approach allows the combined contributions of layer thickness,interbedded structure,impedance contrast and interferences to frequency-dependent refection coefcients and,hence,yielding seismograms of a layered model with a dispersive and attenuative reservoir.Our numerical results show bulk modulus of fracture fuid signifcantly afects anisotropic attenuation,hence causing frequencydependent refection abnormalities.These implications indicate the study of amplitude versus angle and frequency(AVAF)variations provides insights for better interpretation of refection anomalies and hydrocarbon identifcation in a layered reservoir with vertical transverse isotropy(VTI)dispersive media.

Uncertainty Analysis of Seepage-Induced Consolidation in a Fractured Porous Medium

Numerical modeling of seepage-induced consolidation process usually encounters significant uncertainty in the properties of geotechnical materials.Assessing the effect of uncertain parameters on the performance variability of the seepage consolidation model is of critical importance to the simulation and tests of this process.To this end,the uncertainty and sensitivity analyses are performed on a seepage consolidation model in a fractured porous medium using the Bayesian sparse polynomial chaos expansion(SPCE)method.Five uncertain parameters including Young’s modulus,Poisson’s ratio,and the permeability of the porous matrix,the permeability within the fracture,and Biot’s constant are studied.Bayesian SPCE models for displacement,flow velocity magnitude,and fluid pressure at several reference points are constructed to represent the input-output relationship of the numerical model.Based on these SPCE models,the total and first-order Sobol’indices are computed to quantify the contribution of each uncertain input parameter to the uncertainty of model responses.The results show that at different locations of the porous domain,the uncertain parameters show different effects on the output quantities.At the beginning of the seepage consolidation process,the hydraulic parameters make major contributions to the uncertainty of the model responses.As the process progresses,the effect of hydraulic parameters decreases and is gradually surpassed by the mechanical parameters.This work demonstrates the feasibility to apply Bayesian SPCE approach to the uncertainty and sensitivity analyses of seepage-induced consolidation problems and provides guidelines to the numerical modelling and experimental testing of such problems.

Multi-fracture interactions during two-phase flow of oil and water in deformable tight sandstone oil reservoirs

Tight oil reservoirs are complex geological materials composed of solid matrix,pore structure,and mixed multiple phases of fluids,particularly for oil reservoirs suffering from high content of in situ pressurized water found in China.In this regard,a coupled model considering two-phase flow of oil and water,as well as deformation and damage evolution of porous media,is proposed and validated using associated results,including the oil depletion process,analytical solution of stress shadow effect,and physical experiments on multi-fracture interactions and fracture propagation in unsaturated seepage fields.Then,the proposed model is used to study the behavior of multi-fracture interactions in an unsaturated reservoir in presence of water and oil.The results show that conspicuous interactions exist among multiple induced fractures.Interaction behavior varies from extracted geological profiles of the reservoir due to in situ stress anisotropy.The differential pressures of water and that of oil in different regions of reservoir affect interactions and trajectories of multi-fractures to a considerable degree.The absolute value of reservoir average pressure is a dominant factor affecting fracture interactions and in favor of enhancing fracture network complexity.In addition,difference of reservoir average pressures in different regions of reservoir would promote the fracturing effectiveness.Factors affecting fracture interactions and reservoir treatment effectiveness are quantitatively estimated through stimulated reservoir area.This study confirms the significance of incorporating the two-phase flow process in analyses of multifracture interactions and fracture trajectory predictions during tight sandstone oil reservoir developments.

Fully coupled fluid-solid productivity numerical simulation of multistage fractured horizontal well in tight oil reservoirs

A mathematical model, fully coupling multiple porous media deformation and fluid flow, was established based on the elastic theory of porous media and fluid-solid coupling mechanism in tight oil reservoirs. The finite element method was used to determine the numerical solution and the accuracy of the model was verified. On this basis, the model was used to simulate productivity of multistage fractured horizontal wells in tight oil reservoirs. The results show that during the production of tight oil wells, the reservoir region close to artificial fractures deteriorated in physical properties significantly, e.g. the aperture and conductivity of artificial fractures dropped by 52.12% and 89.02% respectively. The simulations of 3000-day production of a horizontal well in tight oil reservoir showed that the predicted productivity by the uncoupled model had an error of 38.30% from that by the fully-coupled model. Apparently, ignoring the influence of fluid-solid interaction effect led to serious deviations of the productivity prediction results. The productivity of horizontal well in tight oil reservoir was most sensitive to the start-up pressure gradient, and second most sensitive to the opening of artificial fractures. Enhancing the initial conductivity of artificial fractures was helpful to improve the productivity of tight oil wells. The influence of conductivity, spacing, number and length of artificial fractures should be considered comprehensively in fracturing design. Increasing the number of artificial fractures unilaterally could not achieve the expected increase in production.

Geological Causes of Permeability Heterogeneity in a Fractured Sandstone-Mudstone Rock Mass in the Xiaolangdi Dam Site,Central China

Heterogeneity of permeability in fractured media is a hot research topic in hydrogeology and numerous approaches had been proposed to characterize heterogeneity in the past several decades.However。