Stability analysis of underground constructions requires a model study of rock masses’ long-term performance. Creep tests under different stress conditions was conducted on intact granite and granite samples fracture...Stability analysis of underground constructions requires a model study of rock masses’ long-term performance. Creep tests under different stress conditions was conducted on intact granite and granite samples fractured at 30° and 45° angles. The experimental results indicate that the steady creep strain rates of intact and fractured rock present an exponential increase trend with the increase of stress level. A nonlinear creep model is developed based on the experimental results, in which the initial damage caused by fracture together with the damage caused by constant load have been taken into consideration. The fitting analysis results indicated that the model proposed is more accurate at identifying the full creep regions in fractured granite, especially the accelerated stage of creep deformation. The least-square fit error of the proposed creep model is significantly lower than that of Nishihara model by almost an order of magnitude. An analysis of the effects of elastic modulus, viscosity coefficient, and damage factors on fractured rock strain rate and creep strain is conducted. If no consideration is given to the effects of the damage, the proposed nonlinear creep model can degenerate into to the classical Nishihara model.展开更多
Considering the importance of fractured rock aquifers in the hydrogeologic process,this research aimed to analyze the flow regime,internal degree of karstification,and estimate storage volume in fractured rock aquifer...Considering the importance of fractured rock aquifers in the hydrogeologic process,this research aimed to analyze the flow regime,internal degree of karstification,and estimate storage volume in fractured rock aquifers of the Germi Chai Basin in northwest Iran,which is attributed to its active tectonics,erosion,and the lithological diversity.Given the geological setting,the hypothesis is that this basin is characterized by a high degree of karstification and diffuse or intermediate flow regime leading to variation in discharge flow rate.The hydrodynamic and hadrochemical analysis was conducted on 9 well distributed springs across the basin from 2019 to 2020.The maximum flow rate in most of the springs appeared in the early wet season despite their different levels of fluctuations on the monthly discharge time series.Analyzing the spring recession curve form revealed an aquifer containing multiple micro-regimes withαrecession coefficients and a degree of karstification ranging between 0.001 to 0.06 and 0.55 to 2.61,respectively.These findings indicated a dominant diffuse and intermediate flow system resulting from the development of a high density of fractures in this area.The electrical conductivity of the spring changes inversely proportional to the change in flow discharge,indicating the reasonable hydrological response of the aquifer to rainfall events.Hydrograph analysis revealed that the delay time of spring discharge after rainfall events mostly varies between 10 to 30 days.The total dynamic storage volume of the spring for a given period(2019-2020)was estimated to be approximately 1324 million cubic meters reflecting the long-term drainage potential and high perdurability of dynamic storage.Estimating the maximum and minimum ratio revealed that the springs recharging system in Germi Chai Basin comes under the slow aquifers category.This finding provides valuable insight into the hydrogeological properties of fractured rock aquifers contributing to effective water management strategy.展开更多
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 thinni...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.展开更多
BACKGROUND Although percutaneous vertebral augmentation(PVA)is a commonly used procedure for treating vertebral compression fracture(VCF),the risk of vertebral refracture should be considered.Chronic kidney disease-mi...BACKGROUND Although percutaneous vertebral augmentation(PVA)is a commonly used procedure for treating vertebral compression fracture(VCF),the risk of vertebral refracture should be considered.Chronic kidney disease-mineral and bone disorder(CKD-MBD)is a systemic disease of mineral and bone metabolism.It is associated with an increased risk of fracture.Few studies have reported the use of PVA in patients with CKD-MBD.We herein report a rare case wherein the cemented vertebra and the adjacent vertebra refractured simultaneously in a CKD-MBD patient after PVA.CASE SUMMARY A 74-year-old man suffered from low back pain after taking a fall about 3 wk ago.According to physical examination,imaging and laboratory findings,diagnoses of T12 VCF,CKD-MBD,and chronic kidney disease stage 5 were established.He then received percutaneous vertebroplasty at T12 vertebra.Fourteen weeks later,he presented with T12 and L1 vertebral refractures caused by lumbar sprain.Once again,he was given PVA which was optimized for the refractured vertebrae.Although the short-term postoperative effect was satisfactory,he reported chronic low back pain again at the 3-month follow-up.CONCLUSION It is necessary that patients with CKD-MBD who have received PVA are aware of the adverse effects of CKD-MBD.It may increase the risk of vertebral refracture.Furthermore,the PVA surgical technique needs to be optimized according to the condition of the patient.The medium-and long-term effects of PVA remain uncertain in patients with CKD-MBD.展开更多
Grouting is an effective method to improve the integrity and stability of fractured rocks that surround deep roadways.After years of research and practice,various theories and a complete set of grouting technologies f...Grouting is an effective method to improve the integrity and stability of fractured rocks that surround deep roadways.After years of research and practice,various theories and a complete set of grouting technologies for deep roadways with fractured rocks have been developed and are widely applied in Chinese coal mining production.This paper systematically summarizes and analyzes the research results concerning the theory,design,materials,processes,and equipment for the grouting and reinforcement of fractured rocks surrounding deep roadways.Specifically,in terms of grouting methods,pregrouting,groutingwhile-excavation,and postgrouting methods are explored;in terms of grouting theory,backfill grouting,compaction grouting,infiltration grouting,and fracture grouting theories are studied.In addition,this paper also studies grouting borehole arrangement,water-cement ratio,grouting pressure,grouting volume,grout diffusion radius,and other grouting parameters and their determination methods.On this basis,this paper explores the physical and mechanical properties of organic and organic-inorganic composite grouting materials,and assess grouting reinforcement quality testing methods and instruments.Taken as the field cases,the application of pregrouting in front of heading faces,groutingwhile-excavation,and postgrouting in the Kouzidong coal mine are then introduced,and the effects of the grouting reinforcements are evaluated.This paper proposes a development direction for grouting technology based on problems existing in the grouting reinforcement of fractured rocks surrounding deep roadways.展开更多
Naturally fractured reservoirs make important contributions to global oil and gas reserves and production.The modeling and simulation of naturally fractured reservoirs are different from conventional reservoirs as the...Naturally fractured reservoirs make important contributions to global oil and gas reserves and production.The modeling and simulation of naturally fractured reservoirs are different from conventional reservoirs as the existence of natural fractures.To address the development optimization problem of naturally fractured reservoirs,we propose an optimization workflow by coupling the optimization methods with the embedded discrete fracture model(EDFM).Firstly,the effective and superior performance of the workflow is verified based on the conceptual model.The stochastic simplex approximate gradient(StoSAG)algorithm,the ensemble optimization(EnOpt)algorithm,and the particle swarm optimization(PSO)algorithm are implemented for the production optimization of naturally fractured reservoirs based on the improved versions of the Egg model and the PUNQ-S3 model.The results of the two cases demonstrate the effectiveness of this optimization workflow by finding the optimal well controls which yield the maximum net present value(NPV).Compared to the initial well control guess,the final NPV obtained from the production optimization of fractured reservoirs based on all three optimization algorithms is significantly enhanced.Compared with the optimization results of the PSO algorithm,StoSAG and EnOpt have significant advantages in terms of final NPV and computational efficiency.The results also show that fractures have a significant impact on reservoir production.The economic efficiency of fractured reservoir development can be significantly improved by the optimization workflow.展开更多
The double-network prepared with an in-situ monomer gel and a fast-crosslinked Cr(III) gel is introduced to develop a thixotropic and high-strength gel (THSG), which is found to have many advantages over the tradition...The double-network prepared with an in-situ monomer gel and a fast-crosslinked Cr(III) gel is introduced to develop a thixotropic and high-strength gel (THSG), which is found to have many advantages over the traditional gels. The THSG gel demonstrates remarkable thermal stability, and no syneresis is observed after 12 months with high salinity brine (95,500 mg/L). Moreover, the SEM and XRD results indicate that the gel is intercalated into the lamellar structures of Na-MMT, where the gel can form a uniform and compact structure. In addition, the THSG gel has an excellent swelling behavior, even in the high salinity brine. In the slim tube experiments, the THSG gel exhibits high rupture pressure and improves blocking capacity after being ruptured. The core flooding results show that a layer of gel filter cake is formed on the face of the fracture, which may be promoted by a high matrix permeability, a small aperture fracture, and a high injection rate. After the gel treatment, the fracture can be completely blocked by the THSG gel. It is found that a high incremental oil recovery (65.3%) can be achieved when the fracture was completely blocked, compared to 40.2% if the gel is ruptured. Although the swelling of ruptured gel can improve oil recovery, part of the injected brine may be channeled through the gel-filled fractures, resulting in a decrease in the sweep efficiency. Therefore, the improved blocking ability by gel swelling (e.g., in fresh water) may be less efficient to contribute to an enhancement of oil recovery. It is also found that the pressure gradient and residual resistance factor to water (Frrw) are higher if the matrix is less permeable, indicating that the fractured reservoir with lower matrix permeability may require a higher gel strength for treatment. The findings of this study may provide novel insights on designing robust double network gels for water shutoff in fractured reservoirs.展开更多
Fractured hydrate-bearing reservoirs show significantly anisotropic geophysical properties. The joint application of seismic and electromagnetic explorations is expected to accurately assess hydrate resources in the f...Fractured hydrate-bearing reservoirs show significantly anisotropic geophysical properties. The joint application of seismic and electromagnetic explorations is expected to accurately assess hydrate resources in the fractured reservoirs. However, the anisotropic joint elastic-electrical properties in such reservoirs that are the key to the successful application of the joint explorations, remain poorly understood. To obtain such knowledge, we designed and implemented dedicated laboratory experiments to study the anisotropic joint elastic-electrical properties in fractured artificial silica sandstones (with fracture density of about 6.2%, porosity of approximately 25.7%, and mean grainsize of 0.089 mm) with evolving methane hydrate. The experimental results showed that the anisotropic compressional wave velocities respectively increased and decreased with the forming and dissociating hydrate, and the variation in the increasing trend and the decreasing extent of the velocity perpendicular to the fractures were more significant than that parallel to the fractures, respectively. The experimental results also showed that the overall decreasing trend of the electrical conductivity parallel to the fractures was steeper than that perpendicular to the fractures during hydrate formation, and the general variations of the two conductivities with complex trend were similar during hydrate dissociation. The variations in the elastic and electrical anisotropic parameters with forming and dissociating hydrate were also found to be distinct. Interpretation of the experimental results suggested that the hydrate binding to the grains evolved to bridge the surfaces of fractures when saturation exceeded 10% during hydrate formation, and the bridging hydrate gradually evolved to floating in fractures during dissociation. The experimental results further showed that the anisotropic velocities and electrical conductivities were correlated with approximately consistent trends of different slopes during hydrate formation, and the joint elastic-electrical anisotropic parameters exhibited a sharp peak at the hydrate saturation of about 10%. The results suggested that the anisotropic joint properties can be employed not only to accurately estimate hydrate saturation but also possibly to identify hydrate distribution in the fractures.展开更多
Reliable estimation of deformation and failure behaviors of fractured rock mass is important for practical engineering design.This study proposes a multi-domain equivalent method for fracture network to estimate the d...Reliable estimation of deformation and failure behaviors of fractured rock mass is important for practical engineering design.This study proposes a multi-domain equivalent method for fracture network to estimate the deformation properties of complex fractured rock mass.It comprehends both the advantages of the discrete fracture network model and the equivalent continuum model to capture the features of discontinuities explicitly while reducing computational intensity.The complex fracture network is stochastically split into a number of subfracture networks according to the domain,length or angle.An analytical solution is derived to infer theoretically the relationship between the elastic moduli of the original complex fractured rock mass and the split subfractured rock masses by introducing a correction term based on the deformation superposition principle.Numerical simulations are conducted to determine the elastic moduli of split subfractured rock masses using universal distinct element code(UDEC),while the elastic modulus of the original model is estimated based on the currently proposed analytical relationship.The results show that the estimation accuracy with the current domainbased splitting model is far superior compared to those with the other two splitting models.Thus,the estimation method of elastic modulus of complex fractured rock mass based on domain splitting mode of fracture network is identified as the multi-domain equivalent method proposed in this paper.The reliability of this method is evaluated,and its high computational efficiency is demonstrated through exemplification with regard to different geometric configurations for stochastically artificial discrete fracture network.The proposed multi-domain equivalent method constructs the theoretical framework except for the regression analysis hypothesis compared to the density-reduced model equivalent method.展开更多
In this paper,a viscoelasticity-plastic damage constitutive equation for naturally fractured shale is deduced,coupling nonlinear tensile-shear mixed fracture mode.Dynamic perforation-erosion on fluid re-distribution a...In this paper,a viscoelasticity-plastic damage constitutive equation for naturally fractured shale is deduced,coupling nonlinear tensile-shear mixed fracture mode.Dynamic perforation-erosion on fluid re-distribution among multi-clusters are considered as well.DFN-FEM(discrete fracture network combined with finite element method)was developed to simulate the multi-cluster complex fractures propagation within temporary plugging fracturing(TPF).Numerical results are matched with field injection and micro-seismic monitoring data.Based on geomechanical characteristics of Weiyuan deep shale gas reservoir in Sichuan Basin,SW China,a multi-cluster complex fractures propagation model is built for TPF.To study complex fractures propagation and the permeability-enhanced region evolution,intersecting and competition mechanisms between the fractures before and after TPF treatment are revealed.Simulation results show that:fracture from middle cluster is restricted by the fractures from side-clusters,and side-clusters plugging is benefit for multi fractures propagation in uniformity;optimized TPF timing should be delayed within a higher density or strike of natural fractures;Within a reservoir-featured natural fractures distribution,optimized TPF timing for most clustered method is 2/3 of total fluid injection time as the optimal plugging time under different clustering modes.展开更多
Groundwater flow through fractured rocks has been recognized as an important issue in many geotechnical engineering practices.Several key aspects of fundamental mechanisms,numerical modeling and engineering applicatio...Groundwater flow through fractured rocks has been recognized as an important issue in many geotechnical engineering practices.Several key aspects of fundamental mechanisms,numerical modeling and engineering applications of flow in fractured rocks are discussed.First,the microscopic mechanisms of fluid flow in fractured rocks,especially under the complex conditions of non-Darcian flow,multiphase flow,rock dissolution,and particle transport,have been revealed through a com-bined effort of visualized experiments and theoretical analysis.Then,laboratory and field methods of characterizing hydraulic properties(e.g.intrinsic permeability,inertial permeability,and unsaturated flow parameters)of fractured rocks in different flow regimes have been proposed.Subsequently,high-performance numerical simulation approaches for large-scale modeling of groundwater flow in frac-tured rocks and aquifers have been developed.Numerical procedures for optimization design of seepage control systems in various settings have also been proposed.Mechanisms of coupled hydro-mechanical processes and control of flow-induced deformation have been discussed.Finally,three case studies are presented to illustrate the applications of the improved theoretical understanding,characterization methods,modeling approaches,and seepage and deformation control strategies to geotechnical engi-neering projects.展开更多
In the process of using the original key stratum theory to predict the height of a water-flowing fractured zone(WFZ),the influence of rock strata outside the calculation range on the rock strata within the calculation...In the process of using the original key stratum theory to predict the height of a water-flowing fractured zone(WFZ),the influence of rock strata outside the calculation range on the rock strata within the calculation range as well as the fact that the shape of the overburden deformation area will change with the excavation length are ignored.In this paper,an improved key stratum theory(IKS theory)was proposed by fixing these two shortcomings.Then,a WFZ height prediction method based on IKS theory was established and applied.First,the range of overburden involved in the analysis was determined according to the tensile stress distribution range above the goaf.Second,the key stratum in the overburden involved in the analysis was identified through IKS theory.Finally,the tendency of the WFZ to develop upward was determined by judging whether or not the identified key stratum will break.The proposed method was applied and verified in a mining case study,and the reasons for the differences in the development patterns between the WFZs in coalfields in Northwest and East China were also fully explained by this method.展开更多
The production performances of a well with a shale gas reservoir displaying a complex fracture network are simulated.In particular,a micro-seismic cloud diagram is used to describe the fracture network,and accordingly...The production performances of a well with a shale gas reservoir displaying a complex fracture network are simulated.In particular,a micro-seismic cloud diagram is used to describe the fracture network,and accordingly,a production model is introduced based on a multi-scale flow mechanism.A finite volume method is then exploited for the integration of the model equations.The effects of apparent permeability,conductivity,Langmuir volume,and bottom hole pressure on gas well production are studied accordingly.The simulation results show that ignoring the micro-scale flow mechanism of the shale gas leads to underestimating the well gas production.It is shown that after ten years of production,the cumulative gas production difference between the two scenarios with and without considering the micro-scale flow mechanisms is 19.5%.The greater the fracture conductivity,the higher the initial gas production of the gas well and the cumulative gas production.The larger the Langmuir volume,the higher the gas production rate and the cumulative gas production.With the reduction of the bottom hole pressure,the cumulative gas production increases,but the growth rate gradually decreases.展开更多
The oil-water two-phase flow pressure-transient analysis model for polymer flooding fractured well is established by considering the comprehensive effects of polymer shear thinning,shear thickening,convection,diffusio...The oil-water two-phase flow pressure-transient analysis model for polymer flooding fractured well is established by considering the comprehensive effects of polymer shear thinning,shear thickening,convection,diffusion,adsorption retention,inaccessible pore volume and effective permeability reduction.The finite volume difference and Newton iteration methods are applied to solve the model,and the effects of fracture conductivity coefficient,injected polymer mass concentration,initial polymer mass concentration and water saturation on the well-test type curves of polymer flooding fractured wells are discussed.The results show that with the increase of fracture conductivity coefficient,the pressure conduction becomes faster and the pressure drop becomes smaller,so the pressure curve of transitional flow goes downward,the duration of bilinear flow becomes shorter,and the linear flow appears earlier and lasts longer.As the injected polymer mass concentration increases,the effective water phase viscosity increases,and the pressure loss increases,so the pressure and pressure derivative curves go upward,and the bilinear flow segment becomes shorter.As the initial polymer mass concentration increases,the effective water phase viscosity increases,so the pressure curve after the wellbore storage segment moves upward as a whole.As the water saturation increases,the relative permeability of water increases,the relative permeability of oil decreases,the total oil-water two-phase mobility becomes larger,and the pressure loss is reduced,so the pressure curve after the wellbore storage segment moves downward as a whole.The reliability and practicability of this new model are verified by the comparison of the results from simplified model and commercial well test software,and the actual well test data.展开更多
This study presents an avant-garde approach for predicting and optimizing production in tight reservoirs,employing a dual-medium unsteady seepage model specifically fashioned for volumetrically fractured horizontal we...This study presents an avant-garde approach for predicting and optimizing production in tight reservoirs,employing a dual-medium unsteady seepage model specifically fashioned for volumetrically fractured horizontal wells.Traditional models often fail to fully capture the complex dynamics associated with these unconventional reservoirs.In a significant departure from these models,our approach incorporates an initiation pressure gradient and a discrete fracture seepage network,providing a more realistic representation of the seepage process.The model also integrates an enhanced fluid-solid interaction,which allows for a more comprehensive understanding of the fluid-structure interactions in the reservoir.This is achieved through the incorporation of improved permeability and stress coupling,leading to more precise predictions of reservoir behavior.The numerical solutions derived from the model are obtained through the sophisticated finite element method,ensuring high accuracy and computational efficiency.To ensure the model’s reliability and accuracy,the outcomes were tested against a real-world case,with results demonstrating strong alignment.A key revelation from the study is the significant difference between uncoupled and fully coupled volumetrically fractured horizontal wells,challenging conventional wisdom in the field.Additionally,the study delves into the effects of stress,fracture length,and fracture number on reservoir production,contributing valuable insights for the design and optimization of tight reservoirs.The findings from this study have the potential to revolutionize the field of tight reservoir prediction and management,offering significant advancements in petroleum engineering.The proposed approach brings forth a more nuanced understanding of tight reservoir systems and opens up new avenues for optimizing reservoir management and production.展开更多
Tightening the well spacing for unconventional tight reservoirs is an efficient technique to enhance oil and gas recoveries.Infill well-caused fracture connection between wells is widely reported in the field with sma...Tightening the well spacing for unconventional tight reservoirs is an efficient technique to enhance oil and gas recoveries.Infill well-caused fracture connection between wells is widely reported in the field with small well spacing.This will make it difficult to make formation evaluation and fracture characterization between wells compared to single well cases.In this paper,a novel production data analysis(PDA)method is proposed for fracture characterization with the consideration of interwell fracture connections after the hydraulic fracturing of the infill.The PDA method is based on a semianalytical model,in which the small-scaled fractures are treated with the concept of stimulated reservoir volume(SRV).Thus,the fracture connections between wells are classified into three types,including SRV,fractures,and both SRV and fractures.The physical model is discretized into several linear flow regions,so the mathematical model can be solved semianalytically.An integrated workflow is proposed to analyze the production data for the wellpad,and three steps are mainly included in the workflow,including PDA for the parent well before infill,PDA for the parent well after infill,and PDA for the infill well.In each step,the production performance in the early linear and bilinear flow regimes are analyzed with approximate solutions in the square and fourth root-of-time plots.Because only the relationship between unknown model parameters can be obtained with the approximate solutions,history matching to the production data in log-log plots is further used to determine each unknown parameter.The PDA method is benchmarked with a synthetic case generated by the numerical simulator tNavigator and a field case from Southwestern China.The results show that both good matches and precise parameters can be obtained with the proposed PDA method.The connected fracture number will not be sensitive in PDA when the wells are connected with high-conductive dSRV.The innovation of this paper is that a practical method is provided for PDA analysis of well groups with fracture connection,and it will be a good technique for fracture characterization and well-interference analysis for tight formations.展开更多
Spontaneous imbibition is an important phenomenon in tight reservoirs.The existence of a large number of fractures and micro-nano pores is the key factor affecting the spontaneous imbibition of tight reservoirs.In thi...Spontaneous imbibition is an important phenomenon in tight reservoirs.The existence of a large number of fractures and micro-nano pores is the key factor affecting the spontaneous imbibition of tight reservoirs.In this study,based on high-pressure mercury injection and nuclear magnetic resonance experiments,the pore distribution of tight sandstone is described.The influence of fractures,core porosity and permeability,and surfactants on the spontaneous imbibition of tight sandstone are studied by physical fracturing,interfacial tension test,wettability test and imbibition experiments.The results show that:the pore radius of tight sandstone is concentrated in 0.01-1 mm.Fractures can effectively reduce the oil drop adsorption on the core surface,enhancing the imbibition recovery of the tight sandstone with an increase of about 10%.As the number of fractures increases,the number of oil droplets adsorbed on the core surface decrease and the imbibition rate increases.The imbibition recovery increases with the increase in pore connectivity,while the imbibition rate increases with the increases in core porosity and permeability.The surfactant can improve the core water wettability and reduce the oilwater interfacial tension,reducing the adsorption of oil droplets on the core surface,and improving the core imbibition recovery with an increase of about 15%.In a word,the existence of fractures and surfactants can enhance the pore connectivity of the reservoir,reduce the adsorption of oil droplets on the core surface,and improve the imbibition rate and recovery rate of the tight oil reservoir.展开更多
Hydraulic fracturing is widely used in geothermal resource exploitation, and many natural fractures exist in hot dry rock reservoirs due to in-situ stress and faults. However, the infuence of natural fractures on hydr...Hydraulic fracturing is widely used in geothermal resource exploitation, and many natural fractures exist in hot dry rock reservoirs due to in-situ stress and faults. However, the infuence of natural fractures on hydraulic fracture propagation is not considered in the current study. In this paper, based on the phase feld model, a thermo-hydro-mechanical coupled hydraulic fracture propagation model was established to reveal the infuence of injection time, fracturing method, injection fow rate, and natural fracture distribution on the fracture propagation mechanism. The results show that fracture complexity increases with an increase in injection time. The stress disturbance causes the fracture initiation pressure of the second cluster signifcantly higher than that of the frst and third clusters. The zipper-type fracturing method can reduce the degree of stress disturbance and increase fracture complexity by 7.2% compared to simultaneous hydraulic fracturing. Both low and high injection fow rate lead to a decrease in fracture propagation time, which is not conducive to an increase in fracture complexity. An increase in the natural fracture angle leads to hydraulic fracture crossing natural fracture, but has a lesser efect on fracture complexity. In this paper, we analyzed the infuence of diferent factors on initiation pressure and fracture complexity, providing valuable guidance for the exploitation of geothermal resources.展开更多
In coal mining,rock strata are fractured under cyclic loading and unloading to form fracture channels.Fracture channels are the main flow narrows for gas.Therefore,expounding the flow conductivity of fracture channels...In coal mining,rock strata are fractured under cyclic loading and unloading to form fracture channels.Fracture channels are the main flow narrows for gas.Therefore,expounding the flow conductivity of fracture channels in rocks on fluids is significant for gas flow in rock strata.In this regard,graded incremental cyclic loading and unloading experiments were conducted on sandstones with different initial stress levels.Then,the three-dimensional models for fracture channels in sandstones were established.Finally,the fracture channel percentages were used to reflect the flow conductivity of fracture channels.The study revealed how the particle size distribution of fractured sandstone affects the formation and expansion of fracture channels.It was found that a smaller proportion of large blocks and a higher proportion of small blocks after sandstone fails contribute more to the formation of fracture channels.The proportion of fracture channels in fractured rock can indicate the flow conductivity of those channels.When the proportion of fracture channels varies gently,fluids flow evenly through those channels.However,if the proportion of fracture channels varies significantly,it can greatly affect the flow rate of fluids.The research results contribute to revealing the morphological evolution and flow conductivity of fracture channels in sandstone and then provide a theoretical basis for clarifying the gas flow pattern in the rock strata of coal mines.展开更多
Based on imbibition replacement of shut-in well in tight oil reservoirs, this paper expounds the principle of saturation rebalancing during the shut-in process after fracturing, establishes an optimization method for ...Based on imbibition replacement of shut-in well in tight oil reservoirs, this paper expounds the principle of saturation rebalancing during the shut-in process after fracturing, establishes an optimization method for shut-in time after horizontal well volume fracturing with the goal of shortening oil breakthrough time and achieving rapid oil breakthrough, and analyzes the influences of permeability, porosity, fracture half-length and fracturing fluid volume on the shut-in time. The oil and water imbibition displacement in the matrix and fractures occurs during the shut-in process of wells after fracturing. If the shut-in time is too short, the oil-water displacement is not sufficient, and the oil breakthrough time is long after the well is put into production. If the shut-in time is too long, the oil and water displacement is sufficient, but the energy dissipation in the formation near the bottom of the well is severe, and the flowing period is short and the production is low after the well is put into production. A rational shut-in time can help shorten the oil breakthrough time, extend the flowing period and increase the production of the well. The rational shut-in time is influenced by factors such as permeability, porosity, fracture half-length and fracturing fluid volume. The shortest and longest shut-in times are negatively correlated with porosity, permeability, and fracture half-length, and positively correlated with fracturing fluid volume. The pilot test in tight oil horizontal wells in the Songliao Basin, NE China, has confirmed that the proposed optimization method can effectively improve the development effect of horizontal well volume fracturing.展开更多
基金supported by the National Natural Science Foundation of China(No.42307258)the technological research projects in Sichuan Province(No.2022YFSY0007)the China Atomic Energy Authority(CAEA)through the Geological Disposal Program.
文摘Stability analysis of underground constructions requires a model study of rock masses’ long-term performance. Creep tests under different stress conditions was conducted on intact granite and granite samples fractured at 30° and 45° angles. The experimental results indicate that the steady creep strain rates of intact and fractured rock present an exponential increase trend with the increase of stress level. A nonlinear creep model is developed based on the experimental results, in which the initial damage caused by fracture together with the damage caused by constant load have been taken into consideration. The fitting analysis results indicated that the model proposed is more accurate at identifying the full creep regions in fractured granite, especially the accelerated stage of creep deformation. The least-square fit error of the proposed creep model is significantly lower than that of Nishihara model by almost an order of magnitude. An analysis of the effects of elastic modulus, viscosity coefficient, and damage factors on fractured rock strain rate and creep strain is conducted. If no consideration is given to the effects of the damage, the proposed nonlinear creep model can degenerate into to the classical Nishihara model.
基金the Water and Wastewater Company of East Azarbaijan Province for providing the funding for this research
文摘Considering the importance of fractured rock aquifers in the hydrogeologic process,this research aimed to analyze the flow regime,internal degree of karstification,and estimate storage volume in fractured rock aquifers of the Germi Chai Basin in northwest Iran,which is attributed to its active tectonics,erosion,and the lithological diversity.Given the geological setting,the hypothesis is that this basin is characterized by a high degree of karstification and diffuse or intermediate flow regime leading to variation in discharge flow rate.The hydrodynamic and hadrochemical analysis was conducted on 9 well distributed springs across the basin from 2019 to 2020.The maximum flow rate in most of the springs appeared in the early wet season despite their different levels of fluctuations on the monthly discharge time series.Analyzing the spring recession curve form revealed an aquifer containing multiple micro-regimes withαrecession coefficients and a degree of karstification ranging between 0.001 to 0.06 and 0.55 to 2.61,respectively.These findings indicated a dominant diffuse and intermediate flow system resulting from the development of a high density of fractures in this area.The electrical conductivity of the spring changes inversely proportional to the change in flow discharge,indicating the reasonable hydrological response of the aquifer to rainfall events.Hydrograph analysis revealed that the delay time of spring discharge after rainfall events mostly varies between 10 to 30 days.The total dynamic storage volume of the spring for a given period(2019-2020)was estimated to be approximately 1324 million cubic meters reflecting the long-term drainage potential and high perdurability of dynamic storage.Estimating the maximum and minimum ratio revealed that the springs recharging system in Germi Chai Basin comes under the slow aquifers category.This finding provides valuable insight into the hydrogeological properties of fractured rock aquifers contributing to effective water management strategy.
基金This work is supported by the National Natural Science Foundation of China(No.52104049)the Young Elite Scientist Sponsorship Program by Beijing Association for Science and Technology(No.BYESS2023262)Science Foundation of China University of Petroleum,Beijing(No.2462022BJRC004).
文摘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.
文摘BACKGROUND Although percutaneous vertebral augmentation(PVA)is a commonly used procedure for treating vertebral compression fracture(VCF),the risk of vertebral refracture should be considered.Chronic kidney disease-mineral and bone disorder(CKD-MBD)is a systemic disease of mineral and bone metabolism.It is associated with an increased risk of fracture.Few studies have reported the use of PVA in patients with CKD-MBD.We herein report a rare case wherein the cemented vertebra and the adjacent vertebra refractured simultaneously in a CKD-MBD patient after PVA.CASE SUMMARY A 74-year-old man suffered from low back pain after taking a fall about 3 wk ago.According to physical examination,imaging and laboratory findings,diagnoses of T12 VCF,CKD-MBD,and chronic kidney disease stage 5 were established.He then received percutaneous vertebroplasty at T12 vertebra.Fourteen weeks later,he presented with T12 and L1 vertebral refractures caused by lumbar sprain.Once again,he was given PVA which was optimized for the refractured vertebrae.Although the short-term postoperative effect was satisfactory,he reported chronic low back pain again at the 3-month follow-up.CONCLUSION It is necessary that patients with CKD-MBD who have received PVA are aware of the adverse effects of CKD-MBD.It may increase the risk of vertebral refracture.Furthermore,the PVA surgical technique needs to be optimized according to the condition of the patient.The medium-and long-term effects of PVA remain uncertain in patients with CKD-MBD.
基金Innovation and Entrepreneurship Funds of Tiandi Science&Technology Co.Ltd.,Grant/Award Number:2022-2-TD-MS013。
文摘Grouting is an effective method to improve the integrity and stability of fractured rocks that surround deep roadways.After years of research and practice,various theories and a complete set of grouting technologies for deep roadways with fractured rocks have been developed and are widely applied in Chinese coal mining production.This paper systematically summarizes and analyzes the research results concerning the theory,design,materials,processes,and equipment for the grouting and reinforcement of fractured rocks surrounding deep roadways.Specifically,in terms of grouting methods,pregrouting,groutingwhile-excavation,and postgrouting methods are explored;in terms of grouting theory,backfill grouting,compaction grouting,infiltration grouting,and fracture grouting theories are studied.In addition,this paper also studies grouting borehole arrangement,water-cement ratio,grouting pressure,grouting volume,grout diffusion radius,and other grouting parameters and their determination methods.On this basis,this paper explores the physical and mechanical properties of organic and organic-inorganic composite grouting materials,and assess grouting reinforcement quality testing methods and instruments.Taken as the field cases,the application of pregrouting in front of heading faces,groutingwhile-excavation,and postgrouting in the Kouzidong coal mine are then introduced,and the effects of the grouting reinforcements are evaluated.This paper proposes a development direction for grouting technology based on problems existing in the grouting reinforcement of fractured rocks surrounding deep roadways.
基金This study was supported by the National Natural Science Foundation of China(51904323,52174052).
文摘Naturally fractured reservoirs make important contributions to global oil and gas reserves and production.The modeling and simulation of naturally fractured reservoirs are different from conventional reservoirs as the existence of natural fractures.To address the development optimization problem of naturally fractured reservoirs,we propose an optimization workflow by coupling the optimization methods with the embedded discrete fracture model(EDFM).Firstly,the effective and superior performance of the workflow is verified based on the conceptual model.The stochastic simplex approximate gradient(StoSAG)algorithm,the ensemble optimization(EnOpt)algorithm,and the particle swarm optimization(PSO)algorithm are implemented for the production optimization of naturally fractured reservoirs based on the improved versions of the Egg model and the PUNQ-S3 model.The results of the two cases demonstrate the effectiveness of this optimization workflow by finding the optimal well controls which yield the maximum net present value(NPV).Compared to the initial well control guess,the final NPV obtained from the production optimization of fractured reservoirs based on all three optimization algorithms is significantly enhanced.Compared with the optimization results of the PSO algorithm,StoSAG and EnOpt have significant advantages in terms of final NPV and computational efficiency.The results also show that fractures have a significant impact on reservoir production.The economic efficiency of fractured reservoir development can be significantly improved by the optimization workflow.
基金financial support from the Major Scientific and Technological Projects of CNPC under Grant(ZD2019-183-007)is gratefully acknowledge.
文摘The double-network prepared with an in-situ monomer gel and a fast-crosslinked Cr(III) gel is introduced to develop a thixotropic and high-strength gel (THSG), which is found to have many advantages over the traditional gels. The THSG gel demonstrates remarkable thermal stability, and no syneresis is observed after 12 months with high salinity brine (95,500 mg/L). Moreover, the SEM and XRD results indicate that the gel is intercalated into the lamellar structures of Na-MMT, where the gel can form a uniform and compact structure. In addition, the THSG gel has an excellent swelling behavior, even in the high salinity brine. In the slim tube experiments, the THSG gel exhibits high rupture pressure and improves blocking capacity after being ruptured. The core flooding results show that a layer of gel filter cake is formed on the face of the fracture, which may be promoted by a high matrix permeability, a small aperture fracture, and a high injection rate. After the gel treatment, the fracture can be completely blocked by the THSG gel. It is found that a high incremental oil recovery (65.3%) can be achieved when the fracture was completely blocked, compared to 40.2% if the gel is ruptured. Although the swelling of ruptured gel can improve oil recovery, part of the injected brine may be channeled through the gel-filled fractures, resulting in a decrease in the sweep efficiency. Therefore, the improved blocking ability by gel swelling (e.g., in fresh water) may be less efficient to contribute to an enhancement of oil recovery. It is also found that the pressure gradient and residual resistance factor to water (Frrw) are higher if the matrix is less permeable, indicating that the fractured reservoir with lower matrix permeability may require a higher gel strength for treatment. The findings of this study may provide novel insights on designing robust double network gels for water shutoff in fractured reservoirs.
基金financial supports received from the National Natural Science Foundation of China(42174136,41821002 and 41874151)the Shandong Provincial Natural Science Foundation,China(ZR2021JQ14).
文摘Fractured hydrate-bearing reservoirs show significantly anisotropic geophysical properties. The joint application of seismic and electromagnetic explorations is expected to accurately assess hydrate resources in the fractured reservoirs. However, the anisotropic joint elastic-electrical properties in such reservoirs that are the key to the successful application of the joint explorations, remain poorly understood. To obtain such knowledge, we designed and implemented dedicated laboratory experiments to study the anisotropic joint elastic-electrical properties in fractured artificial silica sandstones (with fracture density of about 6.2%, porosity of approximately 25.7%, and mean grainsize of 0.089 mm) with evolving methane hydrate. The experimental results showed that the anisotropic compressional wave velocities respectively increased and decreased with the forming and dissociating hydrate, and the variation in the increasing trend and the decreasing extent of the velocity perpendicular to the fractures were more significant than that parallel to the fractures, respectively. The experimental results also showed that the overall decreasing trend of the electrical conductivity parallel to the fractures was steeper than that perpendicular to the fractures during hydrate formation, and the general variations of the two conductivities with complex trend were similar during hydrate dissociation. The variations in the elastic and electrical anisotropic parameters with forming and dissociating hydrate were also found to be distinct. Interpretation of the experimental results suggested that the hydrate binding to the grains evolved to bridge the surfaces of fractures when saturation exceeded 10% during hydrate formation, and the bridging hydrate gradually evolved to floating in fractures during dissociation. The experimental results further showed that the anisotropic velocities and electrical conductivities were correlated with approximately consistent trends of different slopes during hydrate formation, and the joint elastic-electrical anisotropic parameters exhibited a sharp peak at the hydrate saturation of about 10%. The results suggested that the anisotropic joint properties can be employed not only to accurately estimate hydrate saturation but also possibly to identify hydrate distribution in the fractures.
基金financial support by the National Natural Science Foundation of China(Grant Nos.52008152,U1965204,52061160367,U2067203 and 52008153)Natural Science Foundation of Hebei Province of China(Grant No.E2021202087)Hebei Department of Human Resource(Grant No.E2020050015)。
文摘Reliable estimation of deformation and failure behaviors of fractured rock mass is important for practical engineering design.This study proposes a multi-domain equivalent method for fracture network to estimate the deformation properties of complex fractured rock mass.It comprehends both the advantages of the discrete fracture network model and the equivalent continuum model to capture the features of discontinuities explicitly while reducing computational intensity.The complex fracture network is stochastically split into a number of subfracture networks according to the domain,length or angle.An analytical solution is derived to infer theoretically the relationship between the elastic moduli of the original complex fractured rock mass and the split subfractured rock masses by introducing a correction term based on the deformation superposition principle.Numerical simulations are conducted to determine the elastic moduli of split subfractured rock masses using universal distinct element code(UDEC),while the elastic modulus of the original model is estimated based on the currently proposed analytical relationship.The results show that the estimation accuracy with the current domainbased splitting model is far superior compared to those with the other two splitting models.Thus,the estimation method of elastic modulus of complex fractured rock mass based on domain splitting mode of fracture network is identified as the multi-domain equivalent method proposed in this paper.The reliability of this method is evaluated,and its high computational efficiency is demonstrated through exemplification with regard to different geometric configurations for stochastically artificial discrete fracture network.The proposed multi-domain equivalent method constructs the theoretical framework except for the regression analysis hypothesis compared to the density-reduced model equivalent method.
基金Supported by the National Natural Science Foundation of China(52192622,52204005,U20A20265)Sichuan Outstanding Young Scientific and Technological Talents Project(2022JDJQ0007).
文摘In this paper,a viscoelasticity-plastic damage constitutive equation for naturally fractured shale is deduced,coupling nonlinear tensile-shear mixed fracture mode.Dynamic perforation-erosion on fluid re-distribution among multi-clusters are considered as well.DFN-FEM(discrete fracture network combined with finite element method)was developed to simulate the multi-cluster complex fractures propagation within temporary plugging fracturing(TPF).Numerical results are matched with field injection and micro-seismic monitoring data.Based on geomechanical characteristics of Weiyuan deep shale gas reservoir in Sichuan Basin,SW China,a multi-cluster complex fractures propagation model is built for TPF.To study complex fractures propagation and the permeability-enhanced region evolution,intersecting and competition mechanisms between the fractures before and after TPF treatment are revealed.Simulation results show that:fracture from middle cluster is restricted by the fractures from side-clusters,and side-clusters plugging is benefit for multi fractures propagation in uniformity;optimized TPF timing should be delayed within a higher density or strike of natural fractures;Within a reservoir-featured natural fractures distribution,optimized TPF timing for most clustered method is 2/3 of total fluid injection time as the optimal plugging time under different clustering modes.
基金The financial supports from the National Natural Science Foundation of China(Grant Nos.51988101,51925906 and 52122905)are gratefully acknowledged.
文摘Groundwater flow through fractured rocks has been recognized as an important issue in many geotechnical engineering practices.Several key aspects of fundamental mechanisms,numerical modeling and engineering applications of flow in fractured rocks are discussed.First,the microscopic mechanisms of fluid flow in fractured rocks,especially under the complex conditions of non-Darcian flow,multiphase flow,rock dissolution,and particle transport,have been revealed through a com-bined effort of visualized experiments and theoretical analysis.Then,laboratory and field methods of characterizing hydraulic properties(e.g.intrinsic permeability,inertial permeability,and unsaturated flow parameters)of fractured rocks in different flow regimes have been proposed.Subsequently,high-performance numerical simulation approaches for large-scale modeling of groundwater flow in frac-tured rocks and aquifers have been developed.Numerical procedures for optimization design of seepage control systems in various settings have also been proposed.Mechanisms of coupled hydro-mechanical processes and control of flow-induced deformation have been discussed.Finally,three case studies are presented to illustrate the applications of the improved theoretical understanding,characterization methods,modeling approaches,and seepage and deformation control strategies to geotechnical engi-neering projects.
基金supported by the Key Projects of Natural Science Foundation of China(No.41931284)the Scientific Research Start-Up Fund for High-Level Introduced Talents of Anhui University of Science and Technology(No.2022yjrc21).
文摘In the process of using the original key stratum theory to predict the height of a water-flowing fractured zone(WFZ),the influence of rock strata outside the calculation range on the rock strata within the calculation range as well as the fact that the shape of the overburden deformation area will change with the excavation length are ignored.In this paper,an improved key stratum theory(IKS theory)was proposed by fixing these two shortcomings.Then,a WFZ height prediction method based on IKS theory was established and applied.First,the range of overburden involved in the analysis was determined according to the tensile stress distribution range above the goaf.Second,the key stratum in the overburden involved in the analysis was identified through IKS theory.Finally,the tendency of the WFZ to develop upward was determined by judging whether or not the identified key stratum will break.The proposed method was applied and verified in a mining case study,and the reasons for the differences in the development patterns between the WFZs in coalfields in Northwest and East China were also fully explained by this method.
基金This work was supported by the National Natural Science Foundation of China(Grant No.52004237)Science and Technology Cooperation Project of the CNPC-SWPU Innovation Alliance(Grant No.2020CX020202)the Sichuan Science and Technology Program(No.2022JDJQ0009).
文摘The production performances of a well with a shale gas reservoir displaying a complex fracture network are simulated.In particular,a micro-seismic cloud diagram is used to describe the fracture network,and accordingly,a production model is introduced based on a multi-scale flow mechanism.A finite volume method is then exploited for the integration of the model equations.The effects of apparent permeability,conductivity,Langmuir volume,and bottom hole pressure on gas well production are studied accordingly.The simulation results show that ignoring the micro-scale flow mechanism of the shale gas leads to underestimating the well gas production.It is shown that after ten years of production,the cumulative gas production difference between the two scenarios with and without considering the micro-scale flow mechanisms is 19.5%.The greater the fracture conductivity,the higher the initial gas production of the gas well and the cumulative gas production.The larger the Langmuir volume,the higher the gas production rate and the cumulative gas production.With the reduction of the bottom hole pressure,the cumulative gas production increases,but the growth rate gradually decreases.
基金Supported by the National Natural Science Foundation of China(52104049)Science Foundation of China University of Petroleum,Beijing(2462022BJRC004)。
文摘The oil-water two-phase flow pressure-transient analysis model for polymer flooding fractured well is established by considering the comprehensive effects of polymer shear thinning,shear thickening,convection,diffusion,adsorption retention,inaccessible pore volume and effective permeability reduction.The finite volume difference and Newton iteration methods are applied to solve the model,and the effects of fracture conductivity coefficient,injected polymer mass concentration,initial polymer mass concentration and water saturation on the well-test type curves of polymer flooding fractured wells are discussed.The results show that with the increase of fracture conductivity coefficient,the pressure conduction becomes faster and the pressure drop becomes smaller,so the pressure curve of transitional flow goes downward,the duration of bilinear flow becomes shorter,and the linear flow appears earlier and lasts longer.As the injected polymer mass concentration increases,the effective water phase viscosity increases,and the pressure loss increases,so the pressure and pressure derivative curves go upward,and the bilinear flow segment becomes shorter.As the initial polymer mass concentration increases,the effective water phase viscosity increases,so the pressure curve after the wellbore storage segment moves upward as a whole.As the water saturation increases,the relative permeability of water increases,the relative permeability of oil decreases,the total oil-water two-phase mobility becomes larger,and the pressure loss is reduced,so the pressure curve after the wellbore storage segment moves downward as a whole.The reliability and practicability of this new model are verified by the comparison of the results from simplified model and commercial well test software,and the actual well test data.
文摘This study presents an avant-garde approach for predicting and optimizing production in tight reservoirs,employing a dual-medium unsteady seepage model specifically fashioned for volumetrically fractured horizontal wells.Traditional models often fail to fully capture the complex dynamics associated with these unconventional reservoirs.In a significant departure from these models,our approach incorporates an initiation pressure gradient and a discrete fracture seepage network,providing a more realistic representation of the seepage process.The model also integrates an enhanced fluid-solid interaction,which allows for a more comprehensive understanding of the fluid-structure interactions in the reservoir.This is achieved through the incorporation of improved permeability and stress coupling,leading to more precise predictions of reservoir behavior.The numerical solutions derived from the model are obtained through the sophisticated finite element method,ensuring high accuracy and computational efficiency.To ensure the model’s reliability and accuracy,the outcomes were tested against a real-world case,with results demonstrating strong alignment.A key revelation from the study is the significant difference between uncoupled and fully coupled volumetrically fractured horizontal wells,challenging conventional wisdom in the field.Additionally,the study delves into the effects of stress,fracture length,and fracture number on reservoir production,contributing valuable insights for the design and optimization of tight reservoirs.The findings from this study have the potential to revolutionize the field of tight reservoir prediction and management,offering significant advancements in petroleum engineering.The proposed approach brings forth a more nuanced understanding of tight reservoir systems and opens up new avenues for optimizing reservoir management and production.
基金supported by the Jiangsu Province Carbon Peak Carbon Neutral Technology Innovation Project in China(BE2022034-2)the foundation of the National Natural Science Foundation of China(51974328)+2 种基金the Fundamental Research Funds for the Central Universities(2021QN1005)the Natural Science Foundation of Jiangsu Province(BK20210520)China Postdoctoral Science Foundation(2022M713372).
文摘Tightening the well spacing for unconventional tight reservoirs is an efficient technique to enhance oil and gas recoveries.Infill well-caused fracture connection between wells is widely reported in the field with small well spacing.This will make it difficult to make formation evaluation and fracture characterization between wells compared to single well cases.In this paper,a novel production data analysis(PDA)method is proposed for fracture characterization with the consideration of interwell fracture connections after the hydraulic fracturing of the infill.The PDA method is based on a semianalytical model,in which the small-scaled fractures are treated with the concept of stimulated reservoir volume(SRV).Thus,the fracture connections between wells are classified into three types,including SRV,fractures,and both SRV and fractures.The physical model is discretized into several linear flow regions,so the mathematical model can be solved semianalytically.An integrated workflow is proposed to analyze the production data for the wellpad,and three steps are mainly included in the workflow,including PDA for the parent well before infill,PDA for the parent well after infill,and PDA for the infill well.In each step,the production performance in the early linear and bilinear flow regimes are analyzed with approximate solutions in the square and fourth root-of-time plots.Because only the relationship between unknown model parameters can be obtained with the approximate solutions,history matching to the production data in log-log plots is further used to determine each unknown parameter.The PDA method is benchmarked with a synthetic case generated by the numerical simulator tNavigator and a field case from Southwestern China.The results show that both good matches and precise parameters can be obtained with the proposed PDA method.The connected fracture number will not be sensitive in PDA when the wells are connected with high-conductive dSRV.The innovation of this paper is that a practical method is provided for PDA analysis of well groups with fracture connection,and it will be a good technique for fracture characterization and well-interference analysis for tight formations.
基金This work was supported by the National Natural Science Foundation of China(No.51874320).
文摘Spontaneous imbibition is an important phenomenon in tight reservoirs.The existence of a large number of fractures and micro-nano pores is the key factor affecting the spontaneous imbibition of tight reservoirs.In this study,based on high-pressure mercury injection and nuclear magnetic resonance experiments,the pore distribution of tight sandstone is described.The influence of fractures,core porosity and permeability,and surfactants on the spontaneous imbibition of tight sandstone are studied by physical fracturing,interfacial tension test,wettability test and imbibition experiments.The results show that:the pore radius of tight sandstone is concentrated in 0.01-1 mm.Fractures can effectively reduce the oil drop adsorption on the core surface,enhancing the imbibition recovery of the tight sandstone with an increase of about 10%.As the number of fractures increases,the number of oil droplets adsorbed on the core surface decrease and the imbibition rate increases.The imbibition recovery increases with the increase in pore connectivity,while the imbibition rate increases with the increases in core porosity and permeability.The surfactant can improve the core water wettability and reduce the oilwater interfacial tension,reducing the adsorption of oil droplets on the core surface,and improving the core imbibition recovery with an increase of about 15%.In a word,the existence of fractures and surfactants can enhance the pore connectivity of the reservoir,reduce the adsorption of oil droplets on the core surface,and improve the imbibition rate and recovery rate of the tight oil reservoir.
基金supported by the National Natural Science Foundation of China(52174024).
文摘Hydraulic fracturing is widely used in geothermal resource exploitation, and many natural fractures exist in hot dry rock reservoirs due to in-situ stress and faults. However, the infuence of natural fractures on hydraulic fracture propagation is not considered in the current study. In this paper, based on the phase feld model, a thermo-hydro-mechanical coupled hydraulic fracture propagation model was established to reveal the infuence of injection time, fracturing method, injection fow rate, and natural fracture distribution on the fracture propagation mechanism. The results show that fracture complexity increases with an increase in injection time. The stress disturbance causes the fracture initiation pressure of the second cluster signifcantly higher than that of the frst and third clusters. The zipper-type fracturing method can reduce the degree of stress disturbance and increase fracture complexity by 7.2% compared to simultaneous hydraulic fracturing. Both low and high injection fow rate lead to a decrease in fracture propagation time, which is not conducive to an increase in fracture complexity. An increase in the natural fracture angle leads to hydraulic fracture crossing natural fracture, but has a lesser efect on fracture complexity. In this paper, we analyzed the infuence of diferent factors on initiation pressure and fracture complexity, providing valuable guidance for the exploitation of geothermal resources.
基金This work was financially supported by the National Natural Science Foundation of China(No.52074041)the Chongqing Talent Program(No.cstc2022ycjh-bgzxm0077)the Postgraduate Research and Innovation Foundation of Chongqing,China(No.CYS23060).
文摘In coal mining,rock strata are fractured under cyclic loading and unloading to form fracture channels.Fracture channels are the main flow narrows for gas.Therefore,expounding the flow conductivity of fracture channels in rocks on fluids is significant for gas flow in rock strata.In this regard,graded incremental cyclic loading and unloading experiments were conducted on sandstones with different initial stress levels.Then,the three-dimensional models for fracture channels in sandstones were established.Finally,the fracture channel percentages were used to reflect the flow conductivity of fracture channels.The study revealed how the particle size distribution of fractured sandstone affects the formation and expansion of fracture channels.It was found that a smaller proportion of large blocks and a higher proportion of small blocks after sandstone fails contribute more to the formation of fracture channels.The proportion of fracture channels in fractured rock can indicate the flow conductivity of those channels.When the proportion of fracture channels varies gently,fluids flow evenly through those channels.However,if the proportion of fracture channels varies significantly,it can greatly affect the flow rate of fluids.The research results contribute to revealing the morphological evolution and flow conductivity of fracture channels in sandstone and then provide a theoretical basis for clarifying the gas flow pattern in the rock strata of coal mines.
基金Supported by China National Major Project of Science and Technology(2016ZX05046-004)PetroChina Major Project of Science and Technology(2017B-4905)PetroChina Jilin Oilfield Company Major Project of Science and Technology(JY21A2-12).
文摘Based on imbibition replacement of shut-in well in tight oil reservoirs, this paper expounds the principle of saturation rebalancing during the shut-in process after fracturing, establishes an optimization method for shut-in time after horizontal well volume fracturing with the goal of shortening oil breakthrough time and achieving rapid oil breakthrough, and analyzes the influences of permeability, porosity, fracture half-length and fracturing fluid volume on the shut-in time. The oil and water imbibition displacement in the matrix and fractures occurs during the shut-in process of wells after fracturing. If the shut-in time is too short, the oil-water displacement is not sufficient, and the oil breakthrough time is long after the well is put into production. If the shut-in time is too long, the oil and water displacement is sufficient, but the energy dissipation in the formation near the bottom of the well is severe, and the flowing period is short and the production is low after the well is put into production. A rational shut-in time can help shorten the oil breakthrough time, extend the flowing period and increase the production of the well. The rational shut-in time is influenced by factors such as permeability, porosity, fracture half-length and fracturing fluid volume. The shortest and longest shut-in times are negatively correlated with porosity, permeability, and fracture half-length, and positively correlated with fracturing fluid volume. The pilot test in tight oil horizontal wells in the Songliao Basin, NE China, has confirmed that the proposed optimization method can effectively improve the development effect of horizontal well volume fracturing.
