Effects of fracture evolution and non-Darcy flow on the thermal performance of enhanced geothermal system in 3D complex fractured rock

In fractured geothermal reservoirs,the fracture networks and internal fluid flow behaviors can significantly impact the thermal performance.In this study,we proposed a non-Darcy rough discrete fracture network(NR-DFN)model that can simultaneously consider the fracture evolution and non-Darcy flow dynamics in studying the thermo-hydro-mechanical(THM)coupling processes for heat extraction in geothermal reservoir.We further employed the model on the Habanero enhanced geothermal systems(EGS)project located in Australia.First,our findings illustrate a clear spatial-temporal variation in the thermal stress and pressure perturbations,as well as uneven spatial distribution of shear failure in 3D fracture networks.Activated shear failure is mainly concentrated in the first fracture cluster.Secondly,channeling flow have also been observed in DFNs during heat extraction and are further intensified by the expansion of fractures driven by thermal stresses.Moreover,the combined effect of non-Darcy flow and fracture evolution triggers a rapid decline in the resulting heat rate and temperature.The NR-DFN model framework and the Habanero EGS's results illustrate the importance of both fracture evolution and non-Darcy flow on the efficiency of EGS production and have the potential to promote the development of more sustainable and efficient EGS operations for stakeholders.

Non-Darcy Flow in Molding Sands

Darcy’s law is widely used to describe the flow in porous media in which there is a linear relationship between fluid velocity and pressure gradient. However, it has been found that for high numbers of Reynolds this law ceases to be valid. In this work, the Ergun equation is employed to consider the non-linearity of air velocity with the pressure gradient in casting sands. The contribution of non-linearity to the total flow in terms of a variable defined as a non-Darcy flow fraction is numerically quantified. In addition, the influence of the shape factor of the sand grains on the non-linear flow fraction is analyzed. It is found that for values of the Reynolds number less or equal than 1, the contribution of non-linearity for spherical particles is around 1.15%.

Gas Condensate Two Phase Flow Performance in Porous Media Considering Capillary Number and Non-Darcy Effects

Retrograde condensation frequently occurs during the development of gas condensate reservoirs. The loss of productivity is often observed due to the reduced relative permeability to gas as condensate accumulates near the well bore region. How to describe the condensate blockage effect exactly has been a continuous research topic. However, up to now, the present methods usually over-estimate or underestimate the productivity reduction due to an incorrect understanding of the mechanism of flow in porous medium, which inevitably results in an inaccurate prediction of production performance. It has been found in recent numerous theoretical and experimental studies that capillary number and non-Darcy flow have significant influence on relative permeability in regions near the well bore. The two effects impose opposite impacts on production performance, thus leading to gas condensate flow showing characteristics different from general understanding. It is significant for prediction of performance in gas condensate wells to understand the two effects exactly. The aim of the paper is to describe and analyze the flow dynamics in porous media accurately during the production of gas condensate reservoirs. Based on the description of three-zone flow mechanism, capillary number and non-Darcy effect are incorporated in the analysis of relative permeability, making it possible to describe the effect of condensate blockage. The effect of capillary number and inertial flow on gas and condensate relative permeability is analyzed in detail. Novel Inflow Performance Relation (IPR) models considering high velocity effects are formulated and the contrast analysis of different IPR models is conducted. The result shows that the proposed method can help predict the production performance and productivity more accurately than conventional methods.

Characteristics of non-Darcy flow in low-permeability reservoirs

Well testing is recognized as an effective means of accurately obtaining the formation parameters of low-permeability reservoirs and effectively analyzing the deliverability.Well test models must comply with the particular characteristics of flow in low-permeability reservoirs in order to obtain reasonable well test interpretation.At present,non-Darcy flow in low-permeability reservoirs is attracting much attention.In this study,displacement tests were conducted on typical cores taken from low-permeability reservoirs.Two dimensionless variables were introduced to analyze the collected experimental data.The results of the dimensionless analysis show whether non-Darcy flow happens or not depends on the properties of fluid and porous media and the pressure differential.The combination of the above three parameters was named as dimensionless criteria coefficient（DCC）.When the value of the DCC was lower than a critical Reynolds number（CRN）,the flow could not be well described by Darcy＇s law（so-called non-Darcy flow）,when the DCC was higher than CRN,the flow obeyed Darcy＇s law.Finally,this paper establishes a transient mathematical model considering Darcy flow and non-Darcy flow in low-permeability reservoirs,and proposes a methodology to solve the model.The solution technique,which is based on the Boltzmann transformation,is well suited for solving the flow model of low-permeability reservoirs.Based on the typical curves analysis,it was found that the pressure and its derivative curves were determined by such parameters as non-Darcy flow index and the flow characteristics.The results can be used for well test analysis of low-permeability reservoirs.

Non-Darcy flow in oil accumulation (oil displacing water) and relative permeability and oil saturation characteristics of low-permeability sandstones

Hydrocarbon resources in low-permeability sandstones are very abundant and are extensively distributed. Low-permeability reservoirs show several unique characteristics, including lack of a definite trap boundary or caprock, limited buoyancy effect, complex oil-gas-water distribution, without obvious oil-gas-water interfaces, and relatively low oil （gas） saturation. Based on the simulation experiments of oil accumulation in low-permeability sandstone （oil displacing water）, we study the migration and accumulation characteristics of non-Darcy oil flow, and discuss the values and influencing factors of relative permeability which is a key parameter characterizing oil migration and accumulation in low-permeability sandstone. The results indicate that： 1） Oil migration （oil displacing water） in low- permeability sandstone shows non-Darcy percolation characteristics, and there is a threshold pressure gradient during oil migration and accumulation, which has a good negative correlation with permeability and apparent fluidity; 2） With decrease of permeability and apparent fluidity and increase of fluid viscosity, the percolation curve is closer to the pressure gradient axis and the threshold pressure gradient increases. When the apparent fluidity is more than 1.0, the percolation curve shows modified Darcy flow characteristics, while when the apparent fluidity up＂ non-Darcy percolation curve; 3） Oil-water is less than 1.0, the percolation curve is a ＂concave- two-phase relative permeability is affected by core permeability, fluid viscosity, apparent fluidity, and injection drive force; 4） The oil saturation of low- permeability sandstone reservoirs is mostly within 35%-60%, and the oil saturation also has a good positive correlation with the permeability and apparent fluidity.

A fractal approach to low velocity non-Darcy flow in a low permeability porous medium

In this paper, the mechanism for fluid flow at low velocity in a porous medium is analyzed based on plastic flow of oil in a reservoir and the fractal approach. The analytical expressions for flow rate and velocity of non-Newtonian fluid flow in the low permeability porous medium are derived, and the threshold pressure gradient （TPG） is also obtained. It is notable that the TPG （J） and permeability （K） of the porous medium analytically exhibit the scaling behavior J ~ K-D＇r/（l＋Or）, where DT is the fractal dimension for tortuous capillaries. The fractal characteristics of tortuosity for capillaries should be considered in analysis of non-Darcy flow in a low permeability porous medium. The model predictions of TPG show good agreement with those obtained by the available expression and experimental data. The proposed model may be conducible to a better understanding of the mechanism for nonlinear flow in the low permeability porous medium.

Experimental Study of Seepage Properties of Non-Darcy Flow in Granular Coal Gangues

By using the steady-state seepage method, a patent seepage device together with the MTS815.02 Rock Mechanics Test System is used to test the seepage properties of non-Darcy flow in a granular gangue with five different grain sizes during the compaction. The experimental results show that the seepage properties are not only related to the stress or displacement level, but also to the grain size, the pore structure of the granular gangue, and the current porosity The permeability and the non-Darcy flow coefficient can be fitted respectively by the cubic polynomials and the power functions of the porosity, Formally, the flow in granular gangue satisfies the Forchheimer＇s binomial flow, but under the great axial and confining pressure and owing to the grain＇s crushing, the flow in granular gangues is different from that in rock-fills which are naturallv oiled un. As a result, the non-Darer flow coefficient may be negative.

Non-Darcy flows in layered porous media(LPMs)with contrasting pore space structures

Compared to single layer porous media,fluid flow through layered porous media(LPMs)with contrasting pore space structures is more complex.This study constructed three-dimensional(3-D)pore-scale LPMs with different grain size ratios of 1.20,1.47,and 1.76.The flow behavior in the constructed LPMs and single layer porous media was numerically investigated.A total of 178 numerical experimental data were collected in LPMs and single layer porous media.In all cases,two different flow regimes(i.e.,Darcy and Non-Darcy)were observed.The influence of the interface of layers on Non-Darcy flow behavior in LPMs was analyzed based pore-scale flow data.It was found that the available correlations based on single layer porous media fail to predict the flow behavior in LPMs,especially for LPM with large grain size ratio.The effective permeability,which incorporated the influence of the interface is more accurate than the Kozeny-Carman equation for estimating the Darcy permeability of LPMs.The inertial pressure loss in LPMs,which determines the onset of the Non-Darcy flow,was underestimated when using a power law expression of mean grain size.The constant B,an empirical value in the classical Ergun equation,typically equals 1.75.The inertial pressure loss in LPMs can be significantly different from it in single lager porous media.For Non-Darcy flow in LPMs,it is necessary to consider a modified larger constant B to improve the accuracy of the Ergun empirical equation.

One-dimensional non-Darcy flow in a semi-infinite porous media:a multiphase implicit Stefan problem with phases divided by hydraulic gradients

One-dimensional non-Darcy flow in a semi-infinite porous media is investigated. We indicate that the non-Darcy relation which is usually determined from experimental results can always be described by a piecewise linear function, and the problem can be equivalently transformed to a multiphase implicit Stefan problem. The novel feature of this Stefan problem is that the phases of the porous media are divided by hydraulic gradients, not the excess pore water pressures. Using the similarity transformation technique, an exact solution for the situation that the external load increases in proportion to the square root of time is developed. The study on the existence and uniqueness of the solution leads to the requirement of a group of inequalities. A similar Stefan problem considering constant surface seepage velocity is also investigated, and the solution, which we indicate to be uniquely existent under all conditions, is established. Meanwhile, the relation between our Stefan problem and the traditional multiphase Stefan problem is demonstrated. In the end, computational examples of the solution are presented and discussed. The solution provides a useful benchmark for verifying the accuracy of general approximate algorithms of Stefan problems, and it is also attractive in the context of inverse problem analysis.

Research advances in non-Darcy flow in low permeability media

More and more experimental results show that Darcy’s law is not fully applicable in low permeability media,and non-Darcy flow has been identified.In this paper we reviewed the research of non-Darcy flow experiments in low-permeability media in recent decades,discuss the existence of non-Darcy flow,and summarize its constitutive equations.The reasons for the threshold gradient were also discussed and summarized for the criterion of the critical point of non-Darcy flow.On this basis,the future development of non-Darcy flow experiments in the rock and clay media were discussed,in order to provide a certain reference for subsequent research on seepage laws in low permeability media.

Efect of double dispersion on non-Darcy mixed convective flow over vertical surface embedded in porous medium

A numerical study of a non-Darcy mixed convective heat and mass transfer flow over a vertical surface embedded in a dispersion, melting, and thermal radiation is porous medium under the effects of double investigated. The set of governing boundary layer equations and the boundary conditions is transformed into a set of coupled nonlinear ordinary differential equations with the relevant boundary conditions. The transformed equations are solved numerically by using the Chebyshev pseudospectral method. Comparisons of the present results with the existing results in the literature are made, and good agreement is found. Numerical results for the velocity, temperature, concentration profiles, and local Nusselt and Sherwood numbers are discussed for various values of physical parameters.

Investigation of low water recovery based on gas-water two-phase low-velocity Non-Darcy flow model for hydraulically fractured horizontal wells in shale

Various mechanisms are employed to interpret the low water recovery during the shale-gas production period,such as extra-trapped water in the fracture network,water imbibition due to osmotic pressure and capillary pressure.These lead to the difficulty of water flow,which could be described by lowvelocity non-Darcy's law known as threshold pressure gradient(TPG).In this paper we firstly employ the low-velocity non-Darcy's law to describe the water flow and use Darcy flow accounting for slip flow and free molecular flow mechanisms to model gas flow in the shale formation.The sensitive study using numerical simulation shows that the proposed flow model could model the low fracturing liquid recovery and that large pseudo TPG leads to lower fracturing liquid recovery.Thus,the proposed model would give new insight to model the low water recovery in shale formations.

Numerical Simulation of a Two-Phase Flow with Low Permeability anda Start-Up Pressure Gradient

A new numerical model for low-permeability reservoirs is developed.The model incorporates the nonlinear characteristics of oil-water two-phase flows while taking into account the initiation pressure gradient.Related numerical solutions are obtained using a finite difference method.The correctness of the method is demonstrated using a two-dimensional inhomogeneous low permeability example.Then,the differences in the cumulative oil and water production are investigated for different starting water saturations.It is shown that when the initial water saturation grows,the water content of the block continues to rise and the cumulative oil production gradually decreases.

Investigation of Barree-Conway non-Darcy flow effects on coalbed methane production

Coalbed gas non-Darcy flow has been observed in high permeable fracture systems,and some mathematical and numerical models have been proposed to study the effects of non-Darcy flow using Forchheimer non-Darcy model.However,experimental results show that the assumption of a constant Forchheimer factor may cause some limitations in using Forchheimer model to describe non-Darcy flow in porous media.In order to investigate the effects of non-Darcy flow on coalbed methane production,this work presents a more general coalbed gas non-Darcy flow model according to Barree-Conway equation,which could describe the entire range of relationships between flow velocity and pressure gradient from low to high flow velocity.An expanded mixed finite element method is introduced to solve the coalbed gas non-Darcy flow model,in which the gas pressure and velocity can be approximated simultaneously.Error estimate results indicate that pressure and velocity could achieve first-order convergence rate.Non-Darcy simulation results indicate that the non-Darcy effect is significant in the zone near the wellbore,and with the distance from the wellbore increasing,the non-Darcy effect becomes weak gradually.From simulation results,we have also found that the non-Darcy effect is more significant at a lower bottom-hole pressure,and the gas production from non-Darcy flow is lower than the production from Darcy flow under the same permeable condition.

Peristaltic Blood Flow with Gold Nanoparticles on a Carreau Nanofluid through a Non-Darcian Porous Medium

The article investigates the influences of a variable thermal conductivity and wall slip on a peristaltic motion of Carreau nanofluid. The model is concerned with heat and mass transfer inside asymmetric channel. The blood is considered as the base Carreau non-Newtonian fluid and gold (Au) as nanoparticles stressed upon. The Fronchiener effect of the non-Darcian medium is taken in consideration. The system is stressed upon a strong magnetic field and the Hall currents are completed. The problem is modulated mathematically by a system of non-linear partial differential equations which describe the fluid velocity, temperature and concentration. The system is reformulated under the approximation of long wavelength and low Reynolds number. It is solved on using multi-step differential transform method (Ms-DTM) as a semi-analytical method. A gold nanoparticle has increased the temperature distribution which is of great importance in destroying the cancer cells.

Two thin middle-crust low-velocity zones imaged in the Chuan-Dian region of southeastern Tibetan Plateau and their tectonic implications

Intracrustal low-velocity zones(LVZs)indicate a mechanically weak crust and are widely observed in the southeast margin of the Tibetan Plateau.However,their spatial distributions and formation mechanisms remain controversial.To investigate their distribution and detailed morphology of the LVZs in the southeastern Tibetan Plateau,here we used teleseismic events and continuous waveform data recorded by 40 broadband seismic stations newly deployed in the Sichuan-Yunnan region from December 2018 to October 2020.A total of 12,924 high-quality P-wave receiver functions and 5–40 s fundamental Rayleigh surface wave phase velocity dispersion curves from ambient noise cross-correlation functions were obtained.The Swave velocity model at a depth interval of 0–100 km in the study area was inverted by using the trans-dimensional Markov chain Monte Carlo strategy to jointly invert the complementary data of the receiver function waveform and Rayleigh surface wave phase velocity dispersion.Our results show that there are two separate LVZs(~3.5 km/s)surrounding the rigid Daliangshan subblock at crustal depths of approximately 30–40 km,providing new constraints on the geometry of the LVZs in our study region.The two LVZs obtained in this study may represent the middle crustal flow channels,through which the material in the center of the Tibetan Plateau extrudes to its southeast margin.Blocked by the rigid Sichuan Basin and the spindle-like Daliangshan subblock,the material continues to flow southward through the mechanically weak middle crustal channels surrounding the Daliangshan subblock.In addition,the existence of thin LVZs in the middle crust plays an important role in understanding the decoupling between the upper and lower crust in the study area.It also provides new constraint on the complex tectonic deformation process of the southeastern margin of the Tibetan Plateau caused by the collision and compression of the Indian and the Eurasian plates.

Soft-bionic-fishtail structured triboelectric nanogenerator driven by flow-induced vibration for low-velocity water flow energy harvesting

To adapt to the low-velocity water flow closely related to human life,the natural energy can be efficiently harvested and used to power monitoring devices.Herein,a triboelectric soft fishtail(TE-SFT)driven by flow-induced vibration(FIV)effect is proposed based on the soft material synthesis technology.Specifically,inspired by the fishtail fin,a bluff body with the cross-section of fishtail-shaped is designed,and has a preferable vortex effect by fluid simulation.In power generation part,the triboelectric nanogenerator(TENG)is designed to act as an inertial pendulum structure by geometric method.Under the FIV effect,the TESFT driven by fishtail-shaped bluff body swings like a fish in the water and then brings the inertial pendulum to acquire the oscillation for harvesting energy from low-velocity water flow.The TE-SFT attains an open-circuit voltage(VOC)of 200 V to 313 V at the flow velocities of 0.24 to 0.89 m/s.Additionally,after 30 days of water immersion,the VOC of TE-SFT retains 96.81%.In demonstration,the TE-SFT is applied to power the temperature and humidity sensor through harvesting water flow energy.This work also provides a way for self-powered system based on the TENG and soft bionic fish in water environment.

Bifurcation and catastrophe of seepage flow system in broken rock

The study of dynamical behavior of water or gas flows in broken rock is a basic research topic among a series of key projects about stability control of the surrounding rocks in mines and the prevention of some disasters such as water inrush or gas outburst and the protection of the groundwater resource. It is of great theoretical and engineering importance in respect of promo- tion of security in mine production and sustainable development of the coal industry. According to the non-Darcy property of seepage flow in broken rock dynamic equations of non-Darcy and non-steady flows in broken rock are established. By dimensionless transformation, the solution diagram of steady-states satisfying the given boundary conditions is obtained. By numerical analysis of low relaxation iteration, the dynamic responses corresponding to the different flow parameters have been obtained. The stability analysis of the steady-states indicate that a saddle-node bifurcaton exists in the seepage flow system of broken rock. Consequently, using catastrophe theory, the fold catastrophe model of seepage flow instability has been obtained. As a result, the bifurcation curves of the seepage flow systems with different control parameters are presented and the standard potential function is also given with respect to the generalized state variable for the fold catastrophe of a dynamic system of seepage flow in broken rock.

Characteristics of unsteady flow in porous media while considering threshold pressure gradient with Green＇s function

The flow behavior in porous media with threshold pressure gradient(TPG) is more complex than Darcy flow and the equations of motion, and outer boundary and inner boundary with TPG are also different from Darcy flow for unsteady flow of a producing well in a reservoir. An analytic method to solve this kind of problem is in a need of reestablishment. The classical method of Green's function and Newman product principle in a new way are used to solve the unsteady state flow problems of various shapes of well and reservoir while considering the TPG. Four Green's functions of point, line, band and circle while considering the TPG are achieved. Then, two well models of vertical well and horizontal well are built and simultaneously the function to calculate the moving boundary of each well model is provided. The results show that when considering TPG the pressure field is much different, which has a sudden pressure change, with a moving boundary in it. And the moving boundary of each well model increases with time but slows down rapidly, especially when the TGP is large.

SENSITIVITY COEFFICIENTS OF SINGLE-PHASE FLOW IN LOW-PERMEABILITY HETEROGENEOUS RESERVOIRS

Theoretical equations for computing sensitivity coefficients of wellbore pressures to estimate the reservoir parameters in low-permeability reservoirs conditioning to non-Darcy flow data at low velocity were obtained. It is shown by a lot of numerical calculations that the wellbore pressures are much more sensitive to permeability very near the well than to permeability a few gridblocks away from the well. When an initial pressure gradient existent sensitivity coefficients in the region are closer to the active well than to the observation well. Sensitivity coefficients of observation well at the line between the active well and the observation well are influenced greatly by the initial pressure gradient.