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.展开更多
As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accura...As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.展开更多
Borehole instability in naturally fractured rocks poses significant challenges to drilling.Drilling mud invades the surrounding formations through natural fractures under the difference between the wellbore pressure(P...Borehole instability in naturally fractured rocks poses significant challenges to drilling.Drilling mud invades the surrounding formations through natural fractures under the difference between the wellbore pressure(P w)and pore pressure(P p)during drilling,which may cause wellbore instability.However,the weakening of fracture strength due to mud intrusion is not considered in most existing borehole stability analyses,which may yield significant errors and misleading predictions.In addition,only limited factors were analyzed,and the fracture distribution was oversimplified.In this paper,the impacts of mud intrusion and associated fracture strength weakening on borehole stability in fractured rocks under both isotropic and anisotropic stress states are investigated using a coupled DEM(distinct element method)and DFN(discrete fracture network)method.It provides estimates of the effect of fracture strength weakening,wellbore pressure,in situ stresses,and sealing efficiency on borehole stability.The results show that mud intrusion and weakening of fracture strength can damage the borehole.This is demonstrated by the large displacement around the borehole,shear displacement on natural fractures,and the generation of fracture at shear limit.Mud intrusion reduces the shear strength of the fracture surface and leads to shear failure,which explains that the increase in mud weight may worsen borehole stability during overbalanced drilling in fractured formations.A higher in situ stress anisotropy exerts a significant influence on the mechanism of shear failure distribution around the wellbore.Moreover,the effect of sealing natural fractures on maintaining borehole stability is verified in this study,and the increase in sealing efficiency reduces the radial invasion distance of drilling mud.This study provides a directly quantitative prediction method of borehole instability in naturally fractured formations,which can consider the discrete fracture network,mud intrusion,and associated weakening of fracture strength.The information provided by the numerical approach(e.g.displacement around the borehole,shear displacement on fracture,and fracture at shear limit)is helpful for managing wellbore stability and designing wellbore-strengthening operations.展开更多
Complexities in mechanical behaviours of rock masses mainly stem from inherent discontinuities,which calls for advanced bolt-grouting techniques for stability enhancement.Understanding the mechanical properties of bol...Complexities in mechanical behaviours of rock masses mainly stem from inherent discontinuities,which calls for advanced bolt-grouting techniques for stability enhancement.Understanding the mechanical properties of bolt-grouted fractured rock mass(BGFR)and developing accurate prediction methods are crucial to optimize the BGFR support strategies.This paper establishes a new elastoplastic(E-P)model based on the orthotropic and the Mohr-Coulomb(M-C)plastic-yielding criteria.The elastic parameters of the model were derived through a meso-mechanical analysis of composite materials mechanics(CMM).Laboratory BGFR specimens were prepared and uniaxial compression test and variable-angle shear test considering different bolt arrangements were carried out to obtain the mechanical parameters of the specimens.Results showed that the anisotropy of BGFR mainly depends on the relative volume content of each component material in a certain direction.Moreover,the mechanical parameters deduced from the theory of composite materials which consider the short fibre effect are shown to be in good agreement with those determined by laboratory experiments,and the variation rules maintained good consistency.Last,a case study of a real tunnel project is provided to highlight the effectiveness,validity and robustness of the developed E-P model in prediction of stresses and deformations.展开更多
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.展开更多
The gas field in the Bohai Bay Basin is a fractured metamorphic buried-hill reservoir with dual-media characteristics. The retrograde vaporization mechanism observed in this type of gas condensate reservoir differs si...The gas field in the Bohai Bay Basin is a fractured metamorphic buried-hill reservoir with dual-media characteristics. The retrograde vaporization mechanism observed in this type of gas condensate reservoir differs significantly from that observed in sand gas condensate reservoirs. However, studies on improving the recovery of fractured gas condensate reservoirs are limited;thus, the impact of retrograde vaporization on condensate within fractured metamorphic buried-hill reservoirs remains unclear. To address this gap, a series of gas injection experiments are conducted in pressure-volume-temperature(PVT) cells and long-cores to investigate the retrograde vaporization effect of condensate using different gas injection media in fractured gas condensate reservoirs. We analyze the variation in condensate volume, gas-to-oil ratio, and condensate recovery during gas injection and examine the influence of various gas injection media(CO_(2), N_(2), and dry gas) under different reservoir properties and varying gas injection times. The results demonstrate that the exchange of components between injected gas and condensate significantly influences condensate retrograde vaporization in the formation. Compared with dry gas injection and N_(2) injection,CO_(2) injection exhibits a superior retrograde vaporization effect. At a CO_(2) injection volume of 1 PV, the percentage shrinkage volume of condensate is 13.82%. Additionally, at the maximum retrograde condensation pressure, CO_(2) injection can increase the recovery of condensate by 22.4%. However, the condensate recovery is notably lower in fractured gas condensate reservoirs than in homogeneous reservoirs, owing to the creation of dominant gas channeling by fractures, which leads to decreased condensate recovery. Regarding gas injection timing, the effect of gas injection at reservoir pressure on improving condensate recovery is superior to that of gas injection at the maximum retrograde condensation pressure. This research provides valuable guidance for designing gas injection development plans and dynamic tracking adjustments for fractured gas condensate 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 s...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.展开更多
For the case of a fractured reservoir surrounded by deformable rocks, the appropriateness and applicability of the two common methods of coupling of flow and deformation, explicit (coupled) and implicit (uncoupled) me...For the case of a fractured reservoir surrounded by deformable rocks, the appropriateness and applicability of the two common methods of coupling of flow and deformation, explicit (coupled) and implicit (uncoupled) methods are investigated. The explicit formulation is capable of modelling surrounding media;while the implicit coupling is unable to do so as deformation vector does not appear as a primary variable in the formulation. The governing differential equations and the finite element approximation of the governing equations for each of the methods are presented. Spatial discretization is achieved using the Galerkin method, and temporal discretisation using the finite difference technique. In the explicit model, coupling between flow and deformation is captured through volumetric strain compatibility amongst the phases within the system. In the implicit model, this is achieved by defining the pore space storativity as a function of the formation compressibility and the compressibility of the fluid phases within the pore space. The impact of rock deformability on early, intermediate and late time responses of fractured reservoir is investigated through several numerical examples. Salient features of each formulation are discussed and highlighted. It is shown that the implicit model is unable to capture the constraining effects of a non-yielding, surrounding rock, leading to incorrect projections of reservoir production irrespective of the history matching strategy adopted.展开更多
Natural fractures(NFs)are common in shale and tight reservoirs,where staged multi-cluster fracturing of horizontal wells is a prevalent technique for reservoir stimulation.While NFs and stress interference are recogni...Natural fractures(NFs)are common in shale and tight reservoirs,where staged multi-cluster fracturing of horizontal wells is a prevalent technique for reservoir stimulation.While NFs and stress interference are recognized as significant factors affecting hydraulic fracture(HF)propagation,the combined influence of these factors remains poorly understood.To address this knowledge gap,a novel coupled hydromechanical-damage(HMD)model based on the phase field method is developed to investigate the propagation of multi-cluster HFs in fractured reservoirs.The comprehensive energy functional and control functions are established,while incorporating dynamic fluid distribution between multiple perforation clusters and refined changes in rock mechanical parameters during hydraulic fracturing.The HMD coupled multi-cluster HF propagation model investigates various scenarios,including single HF and single NF,reservoir heterogeneity,single HF and NF clusters,and multi-cluster HFs with NF clusters.The results show that the HMD coupling model can accurately capture the impact of approach angle(θ),stress difference and cementation strength on the interaction of HF and NF.The criterion of the open and cross zones is not fixed.The NF angle(a)is not a decisive parameter to discriminate the interaction.According to the relationship between approach angle(θ)and NF angle(a),the contact relationship of HF can be divided into three categories(θ=a,θ<a,andθ>a).The connected NF can increase the complexity of HF by inducing it to form branch fracture,resulting in a fractal dimension of HF as high as2.1280 at angles of±45°.Inter-fracture interference from the heel to the toe of HF shows the phenomenon of no,strong and weak interference.Interestingly,under the influence of NFs,distant HFs from the injection can become dominant fractures.However,as a gradually increases,inter-fracture stress interference becomes the primary factor influencing HF propagation,gradually superseding the dominance of NF induced fractures.展开更多
The main objective is to optimize the development of shale gas-rich areas by predicting seismic sweet spot parameters in shale reservoirs. We systematically assessed the fracture development, fracture gas content, and...The main objective is to optimize the development of shale gas-rich areas by predicting seismic sweet spot parameters in shale reservoirs. We systematically assessed the fracture development, fracture gas content, and rock brittleness in fractured gas-bearing shale reservoirs. To better characterize gas-bearing shale reservoirs with tilted fractures, we optimized the petrophysical modeling based on the equivalent medium theory. Based on the advantages of shale petrophysical modeling, we not only considered the brittle mineral fraction but also the combined effect of shale porosity, gas saturation, and total organic carbon(TOC) when optimizing the brittleness index. Due to fractures generally functioning as essential channels for fluid storage and movement, fracture density and fracture fluid identification factors are critical geophysical parameters for fractured reservoir prediction. We defined a new fracture gas indication factor(GFI) to detect fracture-effective gas content. A new linear PP-wave reflection coefficient equation for a tilted transversely isotropic(TTI) medium was rederived, realizing the direct prediction of anisotropic fracture parameters and the isotropic elasticity parameters from offset vector tile(OVT)-domain seismic data. Synthetic seismic data experiments demonstrated that the inversion algorithm based on the L_P quasinorm sparsity constraint and the split-component inversion strategy exhibits high stability and noise resistance. Finally, we applied our new prediction method to evaluate fractured gas-bearing shale reservoirs in the Sichuan Basin of China, demonstrating its effectiveness.展开更多
A novel phase-field model for the propagation of mixed-mode hydraulic fractures,characterized by the formation of mixed-mode fractures due to the interactions between fluids and solids,is proposed.In this model,the dr...A novel phase-field model for the propagation of mixed-mode hydraulic fractures,characterized by the formation of mixed-mode fractures due to the interactions between fluids and solids,is proposed.In this model,the driving force for the phase field consists of both tensile and shear components,with the fluid contribution primarily manifesting in the tension driving force.The displacement and pressure are solved simultaneously by an implicit method.The numerical solution's iterative format is established by the finite element discretization and Newton-Raphson(NR)iterative methods.The correctness of the model is verified through the uniaxial compression physical experiments on fluid-pressurized rocks,and the limitations of the hydraulic fracture expansion phase-field model,which only considers mode I fractures,are revealed.In addition,the influence of matrix mode II fracture toughness value,natural fracture mode II toughness value,and fracturing fluid injection rate on the hydraulic fracture propagation in porous media with natural fractures is studied.展开更多
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.展开更多
Recent studies have indicated that the injection of carbon dioxide(CO_(2))can lead to increased oil recovery in fractured shale reservoirs following natural depletion.Despite advancements in understanding mass exchang...Recent studies have indicated that the injection of carbon dioxide(CO_(2))can lead to increased oil recovery in fractured shale reservoirs following natural depletion.Despite advancements in understanding mass exchange processes in subsurface formations,there remains a knowledge gap concerning the disparities in these processes between the matrix and fractures at the pore scale in formations with varying permeability.This study aims to experimentally investigate the CO_(2) diffusion behaviors and in situ oil recovery through a CO_(2) huff‘n’puff process in the Jimsar shale oil reservoir.To achieve this,we designed three matrix-fracture models with different permeabilities(0.074 mD,0.170 mD,and 0.466 mD)and experimented at 30 MPa and 91℃.The oil concentration in both the matrix and fracture was monitored using a low-field nuclear magnetic resonance(LF-NMR)technique to quantify in situ oil recovery and elucidate mass-exchange behaviors.The results showed that after three cycles of CO_(2) huff‘n’puff,the total recovery degree increased from 30.28%to 34.95%as the matrix permeability of the core samples increased from 0.074 to 0.466 mD,indicating a positive correlation between CO_(2) extraction efficiency and matrix permeability.Under similar fracture conditions,the increase in matrix permeability further promoted CO_(2) extraction efficiency during CO_(2) huff‘n’puff.Specifically,the increase in matrix permeability of the core had the greatest effect on the extraction of the first-cycle injection in large pores,which increased from 16.42%to 36.64%.The findings from our research provide valuable insights into the CO_(2) huff‘n’puff effects in different pore sizes following fracturing under varying permeability conditions,shedding light on the mechanisms of CO_(2)-enhanced oil recovery in fractured shale oil reservoirs.展开更多
Under the policy background and advocacy of carbon capture,utilization,and storage(CCUS),CO_(2)-EOR has become a promising direction in the shale oil reservoir industry.The multi-scale pore structure distribution and ...Under the policy background and advocacy of carbon capture,utilization,and storage(CCUS),CO_(2)-EOR has become a promising direction in the shale oil reservoir industry.The multi-scale pore structure distribution and fracture structure lead to complex multiphase flow,comprehensively considering multiple mechanisms is crucial for development and CO_(2) storage in fractured shale reservoirs.In this paper,a multi-mechanism coupled model is developed by MATLAB.Compared to the traditional Eclipse300 and MATLAB Reservoir Simulation Toolbox(MRST),this model considers the impact of pore structure on fluid phase behavior by the modified Peng—Robinson equation of state(PR-EOS),and the effect simultaneously radiate to Maxwell—Stefan(M—S)diffusion,stress sensitivity,the nano-confinement(NC)effect.Moreover,a modified embedded discrete fracture model(EDFM)is used to model the complex fractures,which optimizes connection types and half-transmissibility calculation approaches between non-neighboring connections(NNCs).The full implicit equation adopts the finite volume method(FVM)and Newton—Raphson iteration for discretization and solution.The model verification with the Eclipse300 and MRST is satisfactory.The results show that the interaction between the mechanisms significantly affects the production performance and storage characteristics.The effect of molecular diffusion may be overestimated in oil-dominated(liquid-dominated)shale reservoirs.The well spacing and injection gas rate are the most crucial factors affecting the production by sensitivity analysis.Moreover,the potential gas invasion risk is mentioned.This model provides a reliable theoretical basis for CO_(2)-EOR and sequestration in shale oil reservoirs.展开更多
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.展开更多
In this paper,we implement three scales of fracture integrated prediction study by classifying it to macro-( 1/4/λ),meso-( 1/100λ and 1/4λ) and micro-( 1/100λ) scales.Based on the multi-scales rock physics ...In this paper,we implement three scales of fracture integrated prediction study by classifying it to macro-( 1/4/λ),meso-( 1/100λ and 1/4λ) and micro-( 1/100λ) scales.Based on the multi-scales rock physics modelling technique,the seismic azimuthal anisotropy characteristic is analyzed for distinguishing the fractures of meso-scale.Furthermore,by integrating geological core fracture description,image well-logging fracture interpretation,seismic attributes macro-scale fracture prediction and core slice micro-scale fracture characterization,an comprehensive multi-scale fracture prediction methodology and technique workflow are proposed by using geology,well-logging and seismic multi-attributes.Firstly,utilizing the geology core slice observation(Fractures description) and image well-logging data interpretation results,the main governing factors of fracture development are obtained,and then the control factors of the development of regional macro-scale fractures are carried out via modelling of the tectonic stress field.For the meso-scale fracture description,the poststack geometric attributes are used to describe the macro-scale fracture as well,the prestack attenuation seismic attribute is used to predict the meso-scale fracture.Finally,by combining lithological statistic inversion with superposed results of faults,the relationship of the meso-scale fractures,lithology and faults can be reasonably interpreted and the cause of meso-scale fractures can be verified.The micro-scale fracture description is mainly implemented by using the electron microscope scanning of cores.Therefore,the development of fractures in reservoirs is assessed by valuating three classes of fracture prediction results.An integrated fracture prediction application to a real field in Sichuan basin,where limestone reservoir fractures developed,is implemented.The application results in the study area indicates that the proposed multi-scales integrated fracture prediction method and the technique procedureare able to deal with the strong heterogeneity and multi-scales problems in fracture prediction.Moreover,the multi-scale fracture prediction technique integrated with geology,well-logging and seismic multi-information can help improve the reservoir characterization and sweet-spots prediction for the fractured hydrocarbon reservoirs.展开更多
This study mainly investigates the mechanical mechanism of overlying strata breaking and the development of fractured zones during close-distance coal seam group mining in the Gaojialiang coal mine.First,a mechanical ...This study mainly investigates the mechanical mechanism of overlying strata breaking and the development of fractured zones during close-distance coal seam group mining in the Gaojialiang coal mine.First,a mechanical model for the second"activation"of broken overlying strata is established,and the related mechanical"activation"conditions are obtained.A recursive formula for calculating the separation distance of overlying strata is deduced.Second,a height determining method for predicting the height of fractured zones during close-distance coal seam group mining is proposed based on two values,namely,the separation distance and ultimate subsidence value of overlying strata.This method is applied to calculate the fractured zone heights in nos.20107 and 20307 mining faces.The calculated results are almost equal to the field observation results.Third,a modified formula for calculating the height of a waterflowing fractured zone is proposed.A comparison of the calculated and observed results shows that the errors are small.The height determining method and modified formula not only build a theoretical foundation for water conservation mining at the Gaojialiang coal mine,but also provide a reference for estimating the height of water-flowing fractured zones in other coal mines with similar conditions.展开更多
The damage process of fractured rock mass showed that the fracture in rocks induced roof collapse in Yangchangwan Coal Mine, China. The rock mass was particularly weak and fractured. There occurred 6 large-scale dynam...The damage process of fractured rock mass showed that the fracture in rocks induced roof collapse in Yangchangwan Coal Mine, China. The rock mass was particularly weak and fractured. There occurred 6 large-scale dynamical roof falls in the excavation disturbed zone (EDZ) with the collapsing volume of 216 m^3. First, the field detailed geological environment, regional seismic dynamics, and dynamic instability of roadways were generally investigated. Second, the field multiple-index monitoring measurements for detecting the deep delamination of the roof, convergence deformation, bolt-cable load, acoustic emission (AE) characteristic parameters, total AE events, AE energy-releasing rate, rock mass fracture, and damage were arranged. Finally, according to the time-space-strength relations, a quantitative assessment of the influence of rock-mass damage on the dynamic roof instability was accomplished.展开更多
The objective of this work is to study the gas desorption characteristics of the high-rank intact coal and fractured coal.The gas adsorption,mercury porosimetry and gas desorption experiments were carried out in this ...The objective of this work is to study the gas desorption characteristics of the high-rank intact coal and fractured coal.The gas adsorption,mercury porosimetry and gas desorption experiments were carried out in this study.Then,the theories of thermodynamics,diffusion mechanism and desorption kinetics were used to estimate the gas desorption characteristics.The results of gas adsorption experiments show that the initial isosteric adsorption heat of the intact coal is greater than that of the fractured coal,indicating that the gas molecules desorb more easily from fractured coal than intact coal.Using the mercury porosimetry,we find that the diffusion channels of fractured coal are more developed than those of intact coal.The difficult diffusion form dominates in the intact coal during the gas diffusing,while the easy diffusion form dominates in the fractured coal.The results of gas desorption experiments show that the initial gas desorption volume and velocity of the fractured coal are both greater than those of the intact coal.Using the Fick diffusion law,the study calculates the gas diffusion coefficients of the intact coal and fractured coal.The diffusion coefficients of the fractured coal are 2 times and 10 times greater than those of the intact coal at the time of 0-120 and 0-10 min,respectively.展开更多
Inter-crystalline pores, cavities, and fractures created from diagenetic shrinkage of dolomite are inter-connected each other, forming fine oil- and gas-bearing reservoirs. It is hard to predict these complex fracture...Inter-crystalline pores, cavities, and fractures created from diagenetic shrinkage of dolomite are inter-connected each other, forming fine oil- and gas-bearing reservoirs. It is hard to predict these complex fracture-cavity reservoirs because of their random distribution, different growth timing, and so on. Taking the lacustrine dolomite fracture-pore reservoir in the Lower Cretaceous Xiagou Formation in the Qingxi oilfield within the Jiuquan basin as an example, we put forward a comprehensive geophysical method to predict carbonate fractures.展开更多
基金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.
基金supported by the National Natural Science Foundation of China(Grant No.42277165)the Fundamental Research Funds for the Central Universities,China University of Geosciences(Wuhan)(Grant No.CUGCJ1821)the National Overseas Study Fund(Grant No.202106410040).
文摘As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.
基金financially supported by National Natural Science Foundation of China(Grant Nos.52074312 and 52211530097)CNPC Science and Technology Innovation Foundation(Grant No.2021DQ02-0505).
文摘Borehole instability in naturally fractured rocks poses significant challenges to drilling.Drilling mud invades the surrounding formations through natural fractures under the difference between the wellbore pressure(P w)and pore pressure(P p)during drilling,which may cause wellbore instability.However,the weakening of fracture strength due to mud intrusion is not considered in most existing borehole stability analyses,which may yield significant errors and misleading predictions.In addition,only limited factors were analyzed,and the fracture distribution was oversimplified.In this paper,the impacts of mud intrusion and associated fracture strength weakening on borehole stability in fractured rocks under both isotropic and anisotropic stress states are investigated using a coupled DEM(distinct element method)and DFN(discrete fracture network)method.It provides estimates of the effect of fracture strength weakening,wellbore pressure,in situ stresses,and sealing efficiency on borehole stability.The results show that mud intrusion and weakening of fracture strength can damage the borehole.This is demonstrated by the large displacement around the borehole,shear displacement on natural fractures,and the generation of fracture at shear limit.Mud intrusion reduces the shear strength of the fracture surface and leads to shear failure,which explains that the increase in mud weight may worsen borehole stability during overbalanced drilling in fractured formations.A higher in situ stress anisotropy exerts a significant influence on the mechanism of shear failure distribution around the wellbore.Moreover,the effect of sealing natural fractures on maintaining borehole stability is verified in this study,and the increase in sealing efficiency reduces the radial invasion distance of drilling mud.This study provides a directly quantitative prediction method of borehole instability in naturally fractured formations,which can consider the discrete fracture network,mud intrusion,and associated weakening of fracture strength.The information provided by the numerical approach(e.g.displacement around the borehole,shear displacement on fracture,and fracture at shear limit)is helpful for managing wellbore stability and designing wellbore-strengthening operations.
基金funded by the National Key Research and Development Plan(No.2022YFC3203200)Department of Science and Technology of Guangdong Province(No.2021ZT09G087)the National Natural Science Foundation Project of China(No.42167025).
文摘Complexities in mechanical behaviours of rock masses mainly stem from inherent discontinuities,which calls for advanced bolt-grouting techniques for stability enhancement.Understanding the mechanical properties of bolt-grouted fractured rock mass(BGFR)and developing accurate prediction methods are crucial to optimize the BGFR support strategies.This paper establishes a new elastoplastic(E-P)model based on the orthotropic and the Mohr-Coulomb(M-C)plastic-yielding criteria.The elastic parameters of the model were derived through a meso-mechanical analysis of composite materials mechanics(CMM).Laboratory BGFR specimens were prepared and uniaxial compression test and variable-angle shear test considering different bolt arrangements were carried out to obtain the mechanical parameters of the specimens.Results showed that the anisotropy of BGFR mainly depends on the relative volume content of each component material in a certain direction.Moreover,the mechanical parameters deduced from the theory of composite materials which consider the short fibre effect are shown to be in good agreement with those determined by laboratory experiments,and the variation rules maintained good consistency.Last,a case study of a real tunnel project is provided to highlight the effectiveness,validity and robustness of the developed E-P model in prediction of stresses and deformations.
基金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.
文摘The gas field in the Bohai Bay Basin is a fractured metamorphic buried-hill reservoir with dual-media characteristics. The retrograde vaporization mechanism observed in this type of gas condensate reservoir differs significantly from that observed in sand gas condensate reservoirs. However, studies on improving the recovery of fractured gas condensate reservoirs are limited;thus, the impact of retrograde vaporization on condensate within fractured metamorphic buried-hill reservoirs remains unclear. To address this gap, a series of gas injection experiments are conducted in pressure-volume-temperature(PVT) cells and long-cores to investigate the retrograde vaporization effect of condensate using different gas injection media in fractured gas condensate reservoirs. We analyze the variation in condensate volume, gas-to-oil ratio, and condensate recovery during gas injection and examine the influence of various gas injection media(CO_(2), N_(2), and dry gas) under different reservoir properties and varying gas injection times. The results demonstrate that the exchange of components between injected gas and condensate significantly influences condensate retrograde vaporization in the formation. Compared with dry gas injection and N_(2) injection,CO_(2) injection exhibits a superior retrograde vaporization effect. At a CO_(2) injection volume of 1 PV, the percentage shrinkage volume of condensate is 13.82%. Additionally, at the maximum retrograde condensation pressure, CO_(2) injection can increase the recovery of condensate by 22.4%. However, the condensate recovery is notably lower in fractured gas condensate reservoirs than in homogeneous reservoirs, owing to the creation of dominant gas channeling by fractures, which leads to decreased condensate recovery. Regarding gas injection timing, the effect of gas injection at reservoir pressure on improving condensate recovery is superior to that of gas injection at the maximum retrograde condensation pressure. This research provides valuable guidance for designing gas injection development plans and dynamic tracking adjustments for fractured gas condensate reservoirs.
基金the Major Science and Technology Project of Southwest Oil and Gas Field Company(2022ZD01-02).
文摘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.
文摘For the case of a fractured reservoir surrounded by deformable rocks, the appropriateness and applicability of the two common methods of coupling of flow and deformation, explicit (coupled) and implicit (uncoupled) methods are investigated. The explicit formulation is capable of modelling surrounding media;while the implicit coupling is unable to do so as deformation vector does not appear as a primary variable in the formulation. The governing differential equations and the finite element approximation of the governing equations for each of the methods are presented. Spatial discretization is achieved using the Galerkin method, and temporal discretisation using the finite difference technique. In the explicit model, coupling between flow and deformation is captured through volumetric strain compatibility amongst the phases within the system. In the implicit model, this is achieved by defining the pore space storativity as a function of the formation compressibility and the compressibility of the fluid phases within the pore space. The impact of rock deformability on early, intermediate and late time responses of fractured reservoir is investigated through several numerical examples. Salient features of each formulation are discussed and highlighted. It is shown that the implicit model is unable to capture the constraining effects of a non-yielding, surrounding rock, leading to incorrect projections of reservoir production irrespective of the history matching strategy adopted.
基金supported by the National Natural Science Foundation of China(No.52174045)。
文摘Natural fractures(NFs)are common in shale and tight reservoirs,where staged multi-cluster fracturing of horizontal wells is a prevalent technique for reservoir stimulation.While NFs and stress interference are recognized as significant factors affecting hydraulic fracture(HF)propagation,the combined influence of these factors remains poorly understood.To address this knowledge gap,a novel coupled hydromechanical-damage(HMD)model based on the phase field method is developed to investigate the propagation of multi-cluster HFs in fractured reservoirs.The comprehensive energy functional and control functions are established,while incorporating dynamic fluid distribution between multiple perforation clusters and refined changes in rock mechanical parameters during hydraulic fracturing.The HMD coupled multi-cluster HF propagation model investigates various scenarios,including single HF and single NF,reservoir heterogeneity,single HF and NF clusters,and multi-cluster HFs with NF clusters.The results show that the HMD coupling model can accurately capture the impact of approach angle(θ),stress difference and cementation strength on the interaction of HF and NF.The criterion of the open and cross zones is not fixed.The NF angle(a)is not a decisive parameter to discriminate the interaction.According to the relationship between approach angle(θ)and NF angle(a),the contact relationship of HF can be divided into three categories(θ=a,θ<a,andθ>a).The connected NF can increase the complexity of HF by inducing it to form branch fracture,resulting in a fractal dimension of HF as high as2.1280 at angles of±45°.Inter-fracture interference from the heel to the toe of HF shows the phenomenon of no,strong and weak interference.Interestingly,under the influence of NFs,distant HFs from the injection can become dominant fractures.However,as a gradually increases,inter-fracture stress interference becomes the primary factor influencing HF propagation,gradually superseding the dominance of NF induced fractures.
基金financially supported by the Sichuan Science and Technology Program (Grant No. 2023ZYD0158)the National Natural Science Foundation of China (Grant Nos. 42304147 and 42304076)。
文摘The main objective is to optimize the development of shale gas-rich areas by predicting seismic sweet spot parameters in shale reservoirs. We systematically assessed the fracture development, fracture gas content, and rock brittleness in fractured gas-bearing shale reservoirs. To better characterize gas-bearing shale reservoirs with tilted fractures, we optimized the petrophysical modeling based on the equivalent medium theory. Based on the advantages of shale petrophysical modeling, we not only considered the brittle mineral fraction but also the combined effect of shale porosity, gas saturation, and total organic carbon(TOC) when optimizing the brittleness index. Due to fractures generally functioning as essential channels for fluid storage and movement, fracture density and fracture fluid identification factors are critical geophysical parameters for fractured reservoir prediction. We defined a new fracture gas indication factor(GFI) to detect fracture-effective gas content. A new linear PP-wave reflection coefficient equation for a tilted transversely isotropic(TTI) medium was rederived, realizing the direct prediction of anisotropic fracture parameters and the isotropic elasticity parameters from offset vector tile(OVT)-domain seismic data. Synthetic seismic data experiments demonstrated that the inversion algorithm based on the L_P quasinorm sparsity constraint and the split-component inversion strategy exhibits high stability and noise resistance. Finally, we applied our new prediction method to evaluate fractured gas-bearing shale reservoirs in the Sichuan Basin of China, demonstrating its effectiveness.
基金Project supported by the National Natural Science Foundation of China(No.42202314)。
文摘A novel phase-field model for the propagation of mixed-mode hydraulic fractures,characterized by the formation of mixed-mode fractures due to the interactions between fluids and solids,is proposed.In this model,the driving force for the phase field consists of both tensile and shear components,with the fluid contribution primarily manifesting in the tension driving force.The displacement and pressure are solved simultaneously by an implicit method.The numerical solution's iterative format is established by the finite element discretization and Newton-Raphson(NR)iterative methods.The correctness of the model is verified through the uniaxial compression physical experiments on fluid-pressurized rocks,and the limitations of the hydraulic fracture expansion phase-field model,which only considers mode I fractures,are revealed.In addition,the influence of matrix mode II fracture toughness value,natural fracture mode II toughness value,and fracturing fluid injection rate on the hydraulic fracture propagation in porous media with natural fractures is studied.
基金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.
基金National Natural Science Foundation of China via grant number 52174035,52304048China Postdoctoral Science Foundation(2022M722637)Research and Innovation Fund for Graduate Students of Southwest Petroleum University(2022KYCX026).
文摘Recent studies have indicated that the injection of carbon dioxide(CO_(2))can lead to increased oil recovery in fractured shale reservoirs following natural depletion.Despite advancements in understanding mass exchange processes in subsurface formations,there remains a knowledge gap concerning the disparities in these processes between the matrix and fractures at the pore scale in formations with varying permeability.This study aims to experimentally investigate the CO_(2) diffusion behaviors and in situ oil recovery through a CO_(2) huff‘n’puff process in the Jimsar shale oil reservoir.To achieve this,we designed three matrix-fracture models with different permeabilities(0.074 mD,0.170 mD,and 0.466 mD)and experimented at 30 MPa and 91℃.The oil concentration in both the matrix and fracture was monitored using a low-field nuclear magnetic resonance(LF-NMR)technique to quantify in situ oil recovery and elucidate mass-exchange behaviors.The results showed that after three cycles of CO_(2) huff‘n’puff,the total recovery degree increased from 30.28%to 34.95%as the matrix permeability of the core samples increased from 0.074 to 0.466 mD,indicating a positive correlation between CO_(2) extraction efficiency and matrix permeability.Under similar fracture conditions,the increase in matrix permeability further promoted CO_(2) extraction efficiency during CO_(2) huff‘n’puff.Specifically,the increase in matrix permeability of the core had the greatest effect on the extraction of the first-cycle injection in large pores,which increased from 16.42%to 36.64%.The findings from our research provide valuable insights into the CO_(2) huff‘n’puff effects in different pore sizes following fracturing under varying permeability conditions,shedding light on the mechanisms of CO_(2)-enhanced oil recovery in fractured shale oil reservoirs.
基金supported by the National Natural Science Foundation of China(No.52174038 and No.52004307)China Petroleum Science and Technology Project-Major Project-Research on Tight Oil-Shale Oil Reservoir Engineering Methods and Key Technologies in Ordos Basin(No.ZLZX2020-02-04)Science Foundation of China University of Petroleum,Beijing(No.2462018YJRC015)。
文摘Under the policy background and advocacy of carbon capture,utilization,and storage(CCUS),CO_(2)-EOR has become a promising direction in the shale oil reservoir industry.The multi-scale pore structure distribution and fracture structure lead to complex multiphase flow,comprehensively considering multiple mechanisms is crucial for development and CO_(2) storage in fractured shale reservoirs.In this paper,a multi-mechanism coupled model is developed by MATLAB.Compared to the traditional Eclipse300 and MATLAB Reservoir Simulation Toolbox(MRST),this model considers the impact of pore structure on fluid phase behavior by the modified Peng—Robinson equation of state(PR-EOS),and the effect simultaneously radiate to Maxwell—Stefan(M—S)diffusion,stress sensitivity,the nano-confinement(NC)effect.Moreover,a modified embedded discrete fracture model(EDFM)is used to model the complex fractures,which optimizes connection types and half-transmissibility calculation approaches between non-neighboring connections(NNCs).The full implicit equation adopts the finite volume method(FVM)and Newton—Raphson iteration for discretization and solution.The model verification with the Eclipse300 and MRST is satisfactory.The results show that the interaction between the mechanisms significantly affects the production performance and storage characteristics.The effect of molecular diffusion may be overestimated in oil-dominated(liquid-dominated)shale reservoirs.The well spacing and injection gas rate are the most crucial factors affecting the production by sensitivity analysis.Moreover,the potential gas invasion risk is mentioned.This model provides a reliable theoretical basis for CO_(2)-EOR and sequestration in shale oil reservoirs.
文摘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.
基金supported by the national oil and gas major project(No.2011ZX05019-008)National Natural Science Foundation of China(No.41574108 and U1262208)presented at the Exploration Geophysics Symposium 2015 of the EAGE Local Chapter China
文摘In this paper,we implement three scales of fracture integrated prediction study by classifying it to macro-( 1/4/λ),meso-( 1/100λ and 1/4λ) and micro-( 1/100λ) scales.Based on the multi-scales rock physics modelling technique,the seismic azimuthal anisotropy characteristic is analyzed for distinguishing the fractures of meso-scale.Furthermore,by integrating geological core fracture description,image well-logging fracture interpretation,seismic attributes macro-scale fracture prediction and core slice micro-scale fracture characterization,an comprehensive multi-scale fracture prediction methodology and technique workflow are proposed by using geology,well-logging and seismic multi-attributes.Firstly,utilizing the geology core slice observation(Fractures description) and image well-logging data interpretation results,the main governing factors of fracture development are obtained,and then the control factors of the development of regional macro-scale fractures are carried out via modelling of the tectonic stress field.For the meso-scale fracture description,the poststack geometric attributes are used to describe the macro-scale fracture as well,the prestack attenuation seismic attribute is used to predict the meso-scale fracture.Finally,by combining lithological statistic inversion with superposed results of faults,the relationship of the meso-scale fractures,lithology and faults can be reasonably interpreted and the cause of meso-scale fractures can be verified.The micro-scale fracture description is mainly implemented by using the electron microscope scanning of cores.Therefore,the development of fractures in reservoirs is assessed by valuating three classes of fracture prediction results.An integrated fracture prediction application to a real field in Sichuan basin,where limestone reservoir fractures developed,is implemented.The application results in the study area indicates that the proposed multi-scales integrated fracture prediction method and the technique procedureare able to deal with the strong heterogeneity and multi-scales problems in fracture prediction.Moreover,the multi-scale fracture prediction technique integrated with geology,well-logging and seismic multi-information can help improve the reservoir characterization and sweet-spots prediction for the fractured hydrocarbon reservoirs.
基金supported by the National Natural Science Foundation of China(Nos.51474137,and 51574154)Shandong Province Natural Science Fund(No.ZR201709180101)+1 种基金Tai’shan Scholar Engineering Construction Fund of Shandong Province of ChinaPostgraduate Technology Innovation Project of Shandong University of Science and Technology(No.SDKDYC 180103).
文摘This study mainly investigates the mechanical mechanism of overlying strata breaking and the development of fractured zones during close-distance coal seam group mining in the Gaojialiang coal mine.First,a mechanical model for the second"activation"of broken overlying strata is established,and the related mechanical"activation"conditions are obtained.A recursive formula for calculating the separation distance of overlying strata is deduced.Second,a height determining method for predicting the height of fractured zones during close-distance coal seam group mining is proposed based on two values,namely,the separation distance and ultimate subsidence value of overlying strata.This method is applied to calculate the fractured zone heights in nos.20107 and 20307 mining faces.The calculated results are almost equal to the field observation results.Third,a modified formula for calculating the height of a waterflowing fractured zone is proposed.A comparison of the calculated and observed results shows that the errors are small.The height determining method and modified formula not only build a theoretical foundation for water conservation mining at the Gaojialiang coal mine,but also provide a reference for estimating the height of water-flowing fractured zones in other coal mines with similar conditions.
基金supported by the National Natural Science Foundation of China (No.10402033 and No.10772144)
文摘The damage process of fractured rock mass showed that the fracture in rocks induced roof collapse in Yangchangwan Coal Mine, China. The rock mass was particularly weak and fractured. There occurred 6 large-scale dynamical roof falls in the excavation disturbed zone (EDZ) with the collapsing volume of 216 m^3. First, the field detailed geological environment, regional seismic dynamics, and dynamic instability of roadways were generally investigated. Second, the field multiple-index monitoring measurements for detecting the deep delamination of the roof, convergence deformation, bolt-cable load, acoustic emission (AE) characteristic parameters, total AE events, AE energy-releasing rate, rock mass fracture, and damage were arranged. Finally, according to the time-space-strength relations, a quantitative assessment of the influence of rock-mass damage on the dynamic roof instability was accomplished.
基金provided by the National Basic Research Program of China(No.2011CB201204)the Natural Science Foundation for the Youth of China(Nos.41202118 and 51204173)
文摘The objective of this work is to study the gas desorption characteristics of the high-rank intact coal and fractured coal.The gas adsorption,mercury porosimetry and gas desorption experiments were carried out in this study.Then,the theories of thermodynamics,diffusion mechanism and desorption kinetics were used to estimate the gas desorption characteristics.The results of gas adsorption experiments show that the initial isosteric adsorption heat of the intact coal is greater than that of the fractured coal,indicating that the gas molecules desorb more easily from fractured coal than intact coal.Using the mercury porosimetry,we find that the diffusion channels of fractured coal are more developed than those of intact coal.The difficult diffusion form dominates in the intact coal during the gas diffusing,while the easy diffusion form dominates in the fractured coal.The results of gas desorption experiments show that the initial gas desorption volume and velocity of the fractured coal are both greater than those of the intact coal.Using the Fick diffusion law,the study calculates the gas diffusion coefficients of the intact coal and fractured coal.The diffusion coefficients of the fractured coal are 2 times and 10 times greater than those of the intact coal at the time of 0-120 and 0-10 min,respectively.
文摘Inter-crystalline pores, cavities, and fractures created from diagenetic shrinkage of dolomite are inter-connected each other, forming fine oil- and gas-bearing reservoirs. It is hard to predict these complex fracture-cavity reservoirs because of their random distribution, different growth timing, and so on. Taking the lacustrine dolomite fracture-pore reservoir in the Lower Cretaceous Xiagou Formation in the Qingxi oilfield within the Jiuquan basin as an example, we put forward a comprehensive geophysical method to predict carbonate fractures.