Impacts of proppant distribution on development of tight oil reservoirs with threshold pressure gradient

Field evidence indicates that proppant distribution and threshold pressure gradient have great impacts on well productivity.Aiming at the development of unconventional oil reservoirs in Triassic Chang-7 Unit,Ordos Basin of China,we presented an integrated workflow to investigate how(1)proppant placement in induced fracture and(2)non-linear flow in reservoir matrix would affect well productivity and fluid flow in the reservoir.Compared with our research before(Yue et al.,2020),here we extended this study into the development of multi-stage fractured horizontal wells(MFHWs)with large-scale complicated fracture geometry.The integrated workflow is based on the finite element method and consists of simulation models for proppant-laden fluid flow,fracture flow,and non-linear seepage flow,respectively.Simulation results indicate that the distribution of proppant inside the induced cracks significantly affects the productivity of the MFHW.When we assign an idealized proppant distribution instead of the real distribution,there will be an overestimation of 44.98%in daily oil rate and 30.63%in cumulative oil production after continuous development of 1000 days.Besides,threshold pressure gradient(TPG)also significantly affects the well performance in tight oil reservoirs.If we simply apply linear Darcy’s law to the reservoir matrix,the overall cumulative oil production can be overrated by 77%after 1000 days of development.In general,this research provides new insights into the development of tight oil reservoirs with TPG and meanwhile reveals the significance of proppant distribution and non-linear fluid flow in the production scenario design.

Numerical investigation of a coupled moving boundary model of radial flow in low-permeable stress-sensitive reservoir with threshold pressure gradient

The threshold pressure gradient and formation stress-sensitive effect as the two prominent physical phenomena in the development of a low-permeable reservoir are both considered here for building a new coupled moving boundary model of radial flow in porous medium. Moreover, the wellbore storage and skin effect are both incorporated into the inner boundary conditions in the model. It is known that the new coupled moving boundary model has strong nonlinearity. A coordinate transformation based fully implicit finite difference method is adopted to obtain its numerical solutions. The involved coordinate transformation can equivalently transform the dynamic flow region for the moving boundary model into a fixed region as a unit circle, which is very convenient for the model computation by the finite difference method on fixed spatial grids. By comparing the numerical solution obtained from other different numerical method in the existing literature, its validity can be verified. Eventually, the effects of permeability modulus, threshold pressure gradient, wellbore storage coefficient, and skin factor on the transient wellbore pressure, the derivative, and the formation pressure distribution are analyzed respectively.

TRANSIENT PRESSURE OF PERCOLATION THROUGH ONE DIMENSION POROUS MEDIA WITH THRESHOLD PRESSURE GRADIENT

This paper studies the transient pressure of percolation during one production and one shutting in one dimension porous media with threshold pressure gradient, the differential equations are derived and solved with numerical computation. Basing on numerical solution, it is analyzed that: 1. the relation between the steady pressure at well bore (or at endpoint) and threshold pressure gradient, shut-in time, and the corresponding formulae are derived; 2, the regulation of transient pressure peak. The result is very useful and will help experiments and applications in the development of low permeability reservoirs with threshold pressure gradient.

Fractal model of spontaneous imbibition in low-permeability reservoirs coupled with heterogeneity of pore seepage channels and threshold pressure

Spontaneous imbibition(SI)is an important mechanism for enhancing oil recovery in low-permeability reservoirs.Due to the strong heterogeneity,and the non-Darcy flow,the construction of SI model for lowpermeability reservoirs is extremely challenging.Commonly,traditional SI models based on single or averaged capillary tortuosity ignore the influence of heterogeneity of pore seepage channels and the threshold pressure(TP)on imbibition.Therefore,in this work,based on capillary model and fractal theory,a mathematical model of characterizing SI considering heterogeneity of pore seepage channels is established.On this basis,the threshold pressure was introduced to determine the pore radius at which the wetted phase can displace oil.The proposed new SI model was verified by imbibition experimental data.The study shows that for weakly heterogeneous cores with permeability of 0-1 m D,the traditional SI model can characterize the imbibition process relatively accurately,and the new imbibition model can increase the coefficient of determination by 1.05 times.However,traditional model has serious deviations in predicting the imbibition recovery for cores with permeability of 10-50 m D.The new SI model coupling with heterogeneity of pore seepage channels and threshold pressure effectively solves this problem,and the determination coefficient is increased from 0.344 to 0.922,which is increased by2.68 times.For low-permeability reservoirs,the production of the oil in transitional pores(0.01-0.1μm)and mesopores(0.1-1μm)significantly affects the imbibition recovery,as the research shows that when the heterogeneity of pore seepage channels is ignored,the oil recovery in transitional pores and mesopores decreases by 7.54%and 4.26%,respectively.Sensitivity analysis shows that increasing interfacial tension,decreasing contact angle,oil-water viscosity ratio and threshold pressure will increase imbibition recovery.In addition,there are critical values for the influence of these factors on the imbibition recovery,which provides theoretical support for surfactant optimization.

NUMERICAL SIMULATION STUDY ON SURFACTANT FLOODING FOR LOW PERMEABILITY OILFIELD IN THE CONDITION OF THRESHOLD PRESSURE

Based on the non-Darcy flow characteristics of surfactant flooding in the low permeability oilfield, considering the changes of threshold pressure and influence of surfactant on convection, diffusion, adsorption and retention, a mathematical model is established for a three-dimensional, two-phase, three-component surfactant flooding. A new treatment for the changes of threshold pressure and a novel correction method for the relative permeability curve in the process of surfactant flooding are derived, which enhances the matching degree between the mathematical model and field practice. The mathematical model was used to perform the numerical simulation study for a pilot test of surfactant flooding in Chao 45 Block of Daqing Oilfield, a proper injection plan was optimized. After the optimized plan was carried out in oilfield, the desirable effects, like pressure-reducing, injection rate increase, and the increase of oil recovery, were achieved. The average oil increase for single well reaches 37%, the ratio of cost to revenue is above 1：4, so the economic effect of scale is promising.

WELL TESTING ANALYSIS FOR HORIZONTAL WELL WITH CONSIDERATION OF THRESHOLD PRESSURE GRADIENT IN TIGHT GAS RESERVOIRS

A fundamental solution for homogeneous reservoir in infinite space is derived by using the point source function with the consideration of the threshold pressure gradient. The fundamental solution of the continuous point source function is then derived based on the Green function. Various boundary conditions of the reservoirs are considered for this case and the corresponding solutions are obtained through the mirror image reflection and the principle of superimposition. The line source solution is obtained by integration. Subsequently, the horizontal-well bottom hole pressure response function for a non-linear gas flow in the homogeneous gas reservoir is obtained, and the response curve of the dimensionless bottom hole pressure and the derivative for a horizontal well in the homogeneous gas reservoir are obtained. In the end, the sensitivities of the relevant parameters are analyzed, The well test model presented in this paper can be used as the basis of the horizontal well test analysis for tight gas reservoirs.

Dual-porosity model of rate transient analysis for horizontal well in tight gas reservoirs with consideration of threshold pressure gradient

Most researches of the threshold pressure gradient in tight gas reservoirs are experimental and mainly focus on the transient pressure response, without paying much attention to the transient rate decline. This paper establishes a dual-porosity rate transient decline model for the horizontal well with consideration of the threshold pressure gradient, which represents the non-Darcy flow in a fracture system. The solution is obtained by employing the Laplace transform and the orthogonal transform. The bi-logarithmic type curves of the dimensionless production rate and derivative are plotted by the Stehfest numerical inversion method. Seven different flow regimes are identified and the effects of the influence factors such as the threshold pressure gradient, the elastic storativity ratio, and the cross flow coefficient are discussed. The presented research could interpret the production behavior more accurately and effectively for tight gas reservoirs.

Discussion of liquid threshold pressure gradient

Some authors believe that a minimum pressure gradient(called threshold pressure gradient(TPG))is required before a liquid starts to flow in a porous medium.In a tight or shale oil formation,this TPG phenomenon becomes more important,as it is more difficult for a fluid to flow.In this paper,experimental data on TPG published in the literature are carefully reviewed.What we found is that a very low flow velocity corresponding to a very low pressure gradient cannot be measured in the experiments.Experiments can only be done above some measurable flow velocities.If these flow velocities and their corresponding pressure gradients are plotted in an XY plot and extrapolated to zero velocity,a non-zero pressure gradient corresponds to this zero velocity.This non-zero pressure gradient is called threshold pressure gradient in the literature.However,in the regime of very low velocity and very low pressure gradient,the data gradually approach to the origin of the plot,demonstrating a non-linear relationship between the pressure gradient and the velocity.But the data do not approach to a point of zero velocity and a threshold pressure gradient.Therefore,the concept of threshold pressure gradient is a result of data misinterpretation of available experimental data.The correct interpretation is that there are two flow regimes:nonlinear flow regime(non-Darcy flow regime)when the pressure gradients are low,and linear flow regime(Darcy flow regime)when the pressure gradient is intermediate or high.The nonlinear flow regime starts from the origin point.As the pressure gradient is increased,the curve becomes a straight line demonstrating the linear flow regime.We have verified our views by first analyzing the causes of non-Darcy flow,and then systematically analyzed typical experimental data and correlations in the literature.We conclude that TPG does not exist.We also use several counter examples to support our conclusion.

Exact analytical solutions for moving boundary problems of one-dimensional flow in semi-infinite porous media with consideration of threshold pressure gradient

By defining new dimensionless variables, nonlinear mathematical models for one-dimensional flow with unknown moving boundaries in semi-infinite porous media are modified to be solved analytically. The exact analytical solutions for both constant-rate and constant-pressure inner boundary constraint problems are obtained by applying the Green＇s function. Two transcendental equations for moving boundary problems are obtained and solved using the Newton-Raphson iteration. The exact analytical solutions are then compared with the approximate solutions. The Pascal＇s approximate formula in reference is fairly accurate for the moving boundary development under the constant-rate condition. But another Pascal＇s approximate formula given in reference is not very robust for constant-pressure condition problems during the early production period, and could lead to false results at the maximum moving boundary distance. Our results also show that, in presence of larger TPG, more pressure drop is required to maintain a constant-rate production. Under the constant-pressure producing condition, the flow rate may decline dramatically due to a large TPG. What＇s more, there exists a maximum distance for a given TPG, beyond which the porous media is not disturbed.

Lumbar interspinous pressure pain threshold values for healthy young men and women and the effect of prolonged fully flexed lumbar sitting posture:An observational study

BACKGROUND Low back pain(LBP)is a common condition with large burden worldwide.Exposure to prolonged sitting with a flexed lumbar posture has been suggested in the literature to be a potential risk factor for self-reported LBP.No study has previously investigated whether exposure to prolonged flexed sitting posture provokes discomfort/pain and decreased interspinous pressure pain thresholds for healthy young men and women without back pain,despite this being a suggested risk factor for LBP.AIM To investigate whether sitting in a prolonged flexed lumbar posture provokes discomfort and lowers interspinous pressure pain thresholds in the lumbar spine for healthy young men and women without previous LBP.METHODS This is a an observational before and after study of 26 participants(13 men,13 women)between 20-35 years old.Algometry was used to examine the pain threshold for pressure applied between spinous processes of the lumbar spine L1-L5.Pressure algometer measures were performed in prone before and after participants were instructed to sit in a fully flexed posture for a maximum of 15 min or until discomfort was experienced in the low back(Borg CR10=7/10).Wilcoxon signed-rank test was used for analyze values from the before and after test conditions.Mann-Whitney U test was used to investigate potential gender difference.RESULTS Fully flexed lumbar spine sitting posture up to 15 min provoked temporary discomfort but the proportion of participants experiencing discomfort 7/10 in the low back was 62%.For all pain pressure threshold locations tested,there was a significant difference for the study population with moderate-large decreased(r=-0.56)pressure pain threshold after exposure to prolonged flexed sitting posture(P<0.01).Comparisons between gender did not show any significant difference.CONCLUSION The result showed that exposure to fully flexed lumbar sitting posture for up to 15 min produced temporary discomfort in the low back in young healthy adults with no previous history of LBP and significantly reduced lumbar interspinous pressure pain thresholds.No gender-based differences were observed.

Electrical Activities and Pressure Pain Threshold in Oral Contraceptives Users and Nonusers

The aim is to evaluate the influence of oral contraceptive intake and menstrual cycle on the electrical activity and pressure pain threshold from anterior temporal and masseter muscles. Twenty-eight women on reproductive age were selected, 13 OC users and 15 nonusers. They were weekly submitted to electromyography and algometry of the anterior temporal and masseter muscles during three consecutive menstrual cycles. Electrical activities at rest position and PPTs of temporal and masseter muscles were not affected by menstrual cycle or by OCs uses. Comparison between groups demonstrated that working side electrical activity was increased in OC users in both muscles, except during lutheal phase for the anterior temporal. However, comparison within weeks did not demonstrate statistical difference. It was suggested that, in healthy women, oral contraceptive use may influence electrical activity, but different phases of the cycle may not.

Tight gas production model considering TPG as a function of pore pressure,permeability and water saturation

Threshold pressure gradient has great importance in efficient tight gas field development as well as for research and laboratory experiments.This experimental study is carried out to investigate the threshold pressure gradient in detail.Experiments are carried out with and without back pressure so that the effect of pore pressure on threshold pressure gradient may be observed.The trend of increasing or decreasing the threshold pressure gradient is totally opposite in the cases of considering and not considering the pore pressure.The results demonstrate that the pore pressure of tight gas reservoirs has great influence on threshold pressure gradient.The effects of other parameters like permeability and water saturation,in the presence of pore pressure,on threshold pressure gradient are also examined which show that the threshold pressure gradient increases with either a decrease in permeability or an increase in water saturation.Two new correlations of threshold pressure gradient on the basis of pore pressure and permeability,and pore pressure and water saturation,are also introduced.Based on these equations,new models for tight gas production are proposed.The gas slip correction factor is also considered during derivation of this proposed tight gas production models.Inflow performance relationship curves based on these proposed models show that production rates and absolute open flow potential are always be overestimated while ignoring the threshold pressure gradients.

The Differing Effects of Nociception and Pain Memory on Pain Thresholds in Participants with and without a History of Injury: A Pretest-Posttest Quasi Experimental Study

Purpose: Memory has been identified as an important protective feature to prevent future injury, but its role has yet to be ascertained. The current study aimed to determine whether there was a difference in pressure pain threshold (PPT) responses between participants with a prior history of injury of lower extremity injury (PSI) and those without (NPSI) when exposed to 1) experimental mechanical pain, 2) short-term memory recall of a painful stimulus, or 3) long-term memory of the pain associated with a prior injury. Subjects and Methods: The study used a pretest-posttest quasi-experimental design. A convenience sample of 59 pain-free participants was recruited from an urban university. Twenty-nine PSI and 30 NPSI were stratified into two groups based on their injury history with PPT values measured at baseline and immediately following each of the three experimental conditions. A repeated measure ANCOVA analysis was conducted for each condition to determine whether there was a difference in PPT responses between the two groups. Results: There was a statistically significant difference in PPT values between the two groups when exposed to experimental pain, F(1,57) = 6.010, p = 0.017, partial η2 = 0.095 and with long-term pain memory, F(1,57) = 4.886, p = 0.031, partial η2 = 0.079. There was no statistically significant difference between groups with short-term pain memory, F(1,57) = 3.925, p = 0.052, partial η2 = 0.064. Conclusions: These findings suggest that pain processing may be altered by pain memory, highlighting the role of experience and memory in the rehabilitation process.

Study of the Seepage Mechanism in Thick Heterogeneous Gas Reservoirs

The seepage mechanism plays a crucial role in low-permeability gas reservoirs.Compared with conventional gas reservoirs,low-permeability sandstone gas reservoirs are characterized by low porosity,low permeability,strong heterogeneity,and high water saturation.Moreover,their percolation mechanisms are more complex.The present work describes a series of experiments conducted considering low-permeability sandstone cores under pressuredepletion conditions(from the Xihu Depression in the East China Sea Basin).It is shown that the threshold pressure gradient of a low-permeability gas reservoir in thick layers is positively correlated with water saturation and negatively correlated with permeability and porosity.The reservoir stress sensitivity is related to permeability and rock composition.Stress sensitivity is generally low when permeability is high or in the early stage of gas reservoir development.It is also shown that in sand conglomerates,especially the more sparsely filled parts,the interstitial materials among the conglomerates can be rapidly dislodged from the skeleton particles under stress.This material can therefore disperse,migrate,and block the pore throat producing serious,stress-sensitive damage.

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.

Analytical solution of a double moving boundary problem for nonlinear flows in one-dimensional semi-infinite long porous media with low permeability

Based on Huang＇s accurate tri-sectional nonlin- ear kinematic equation （1997）, a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction： The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation （PDE） sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient （TPG） has a great effect on the external moving boundary but has little effect on the internal moving boundary.

A static resistance model and the discontinuous pattern of hydrocarbon accumulation in tight oil reservoirs

In exploration for tight oil, the content and saturation of hydrocarbon in the tight reservoir is a key factor for evaluating the reserve. Therefore, it is necessary to study the geological history of hydrocarbon accumulation and the tight oil charging process. However, kinetic models used for petroleum development are not applicable for petroleum exploration. In this study, a static resistance model[ is proposed after analyzing resistances in ultra-slow flow in porous media. Using this model, the disco^atinuous pattern of oil charging is reproduced through incompressible Navier-Stokes equations, the phase field method and the finite element method. This study also explains macroscopic percolation behavior with microscopic flow mechanisms and discusses some issues in ultra-slow flow in a micro/nano pore-throat network. The resistance analysis reveals that capillary resistance and dissipation resistance are dominant factors in the mechanism of oil accumulation in tight reservoirs. Numerical simulations show that pressure thresholds exist and result in discontinuous oil charging. Generally, it is proven that the static model is more applicable than kinetic models in describing oil accumulation in tight reservoirs.

Evaluation of Well Spacing for Primary Development of Fractured Horizontal Wells in Tight Sandstone Gas Reservoirs

Methods for horizontal well spacing calculation in tight gas reservoirs are still adversely affected by the complexity of related control factors,such as strong reservoir heterogeneity and seepage mechanisms.In this study,the stress sensitivity and threshold pressure gradient of various types of reservoirs are quantitatively evaluated through reservoir seepage experiments.On the basis of these experiments,a numerical simulation model(based on the special seepage mechanism)and an inverse dynamic reserve algorithm(with different equivalent drainage areas)were developed.The well spacing ranges of Classes I,II,and III wells in the Q gas field are determined to be 802–1,000,600–662,and 285–400 m,respectively,with their average ranges as 901,631,and 342.5 m,respectively.By considering both the pairs of parallel well groups and series well groups as examples,the reliability of the calculation results is verified.It is shown that the combination of the two models can reduce errors and provide accurate results.

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.

EFFECT OF A MOVING BOUNDARY ON THE FLUID TRANSIENT FLOW IN LOW PERMEABILITY RESERVOIRS

In a low permeability reservoir, the existence of a moving boundary is considered in the study of the transient porous flow with threshold pressure gradient. The transmission of the moving boundary directly indicates the size of the drainage area as well as the apparent influences on the pressure behavior. The nonlinear transient flow mathematical model in which the threshold pressure gradient and the moving boundary are incorporated is solved by advanced mathematical methods. This paper presents some new analytical solutions describing the pressure distribution at a constant rate and the production decline in a constant pressure production with the boundary propagation. It is shown that the greater the threshold pressure gradient, the slower the transmission of the moving boundary, the larger the pressure loss will be, and there is no radial flow in the middle and later phases of the wellface pressure for a well at a constant rate. We have the the maximum moving boundary at a specific drawdown pressure for a low permeability reservoir The greater the threshold pressure gradient, the smaller the maximum moving boundary distance, the quicker the production decline for a well in a constant pressure production will be. The type curve charts for the modern well test analysis and the rate transient analysis with a moving boundary are obtained and the field test and the production data are interpreted as examples to illustrate how to use our new results.