In order to investigate the influence on shale gas well productivity caused by gas transport in nanometer- size pores, a mathematical model of multi-stage fractured horizontal wells in shale gas reservoirs is built, w...In order to investigate the influence on shale gas well productivity caused by gas transport in nanometer- size pores, a mathematical model of multi-stage fractured horizontal wells in shale gas reservoirs is built, which considers the influence of viscous flow, Knudsen diffusion, surface diffusion, and adsorption layer thickness. A dis- crete-fracture model is used to simplify the fracture mod- cling, and a finite element method is applied to solve the model. The numerical simulation results indicate that with a decrease in the intrinsic matrix permeability, Knudsen diffusion and surface diffusion contributions to production become large and cannot be ignored. The existence of an adsorption layer on the nanopore surfaces reduces the effective pore radius and the effective porosity, resulting in low production from fractured horizontal wells. With a decrease in the pore radius, considering the adsorption layer, the production reduction rate increases. When the pore radius is less than 10 nm, because of the combined impacts of Knudsen diffusion, surface diffusion, and adsorption layers, the production of multi-stage fractured horizontal wells increases with a decrease in the pore pressure. When the pore pressure is lower than 30 MPa, the rate of production increase becomes larger with a decrease in pore pressure.展开更多
The Gabes aquifer system,located in southeastern Tunisia,is a crucial resource for supporting local socio-economic activities.Due to its dual porosity structure,is particularly vulnerable to pollution.This study aims ...The Gabes aquifer system,located in southeastern Tunisia,is a crucial resource for supporting local socio-economic activities.Due to its dual porosity structure,is particularly vulnerable to pollution.This study aims to develop a hybrid model that combines the Fracture Aquifer Index(FAI)with the conventional GOD(Groundwater occurrence,Overall lithology,Depth to water table)method,to assess groundwater vulnerability in fractured aquifer.To develop the hybrid model,the classical GOD method was integrated with FAI to produce a single composite index.Each parameter within both GOD and FAI was scored,and a final index was calculated to delineate vulnerable areas.The results show that the study area can be classified into four vulnerability levels:Very low,low,moderate,and high,indicating that approximately 8%of the area exhibits very low vulnerability,29%has low vulnerability,25%falls into the moderate category,and 38%is considered highly vulnerable.The FAI-GOD model further incorporates fracture network characteristics.This refinement reduces the classification to three vulnerability classes:Low,medium,and high.The outcomes demonstrate that 46%of the area is highly vulnerable due to a dense concentration of fractures,while 17%represents an intermediate zone characterized by either shallow or deeper fractures.In contrast,37%corresponds to areas with lightly fractured rock,where the impact on vulnerability is minimal.Multivariate statistical analysis was employed using Principal Components Analysis(PCA)and Hierarchical Cluster Analysis(HCA)on 24 samples across six variables.The first three components account for over 76%of the total variance,reinforcing the significance of fracture dynamics in classifying vulnerability levels.The FAI-GOD model removes the very-low-vulnerability class and expands the spatial extent of low-and high-vulnerability zones,reflecting the dominant influence of fracture networks on aquifer sensitivity.While both indices use a five-class system,FAI-GOD redistributes vulnerability by eliminating very-low-vulnerability areas and amplifying low/high categories,highlighting the critical role of fractures.A strong correlation(R2=0.94)between the GOD and FAI-GOD indices,demonstrated through second-order polynomial regression,confirms the robustness of the FAI-GOD model in accurately predicting vulnerability to pollution.This model provides a useful framework for assessing the vulnerability of complex aquifers and serves as a decision-making tool for groundwater managers in similar areas.展开更多
The “well factory” mode's high-density well placement and multi-stage hydraulic fracturing technology enable efficient development of unconventional oil and gas resources.However,the deployment of platform wells...The “well factory” mode's high-density well placement and multi-stage hydraulic fracturing technology enable efficient development of unconventional oil and gas resources.However,the deployment of platform wells in the “well factory” model results in small wellbore spacing,and the stress disturbances caused by fracturing operations may affect neighboring wells,leading to inter-well interference phenomena that cause casing deformation.This study investigates the issue of inter-well interference causing casing deformation or even failure during multi-stage hydraulic fracturing in the “well factory”model,and predicts high-risk locations for casing failure.A flow-mechanics coupled geomechanical finite element model with retaining geological stratification characteristics was established.Based on the theory of hydraulic fracturing-induced rock fragmentation and fluid action leading to the degradation of rock mechanical properties,the model simulated the four-dimensional evolution of multi-well fracturing areas over time and space,calculating the disturbance in the regional stress field caused by fracturing operations.Subsequently,the stress distribution of multiple well casings at different time points was calculated to predict high-risk locations for casing failure.The research results show that the redistribution of the stress field in the fracturing area increases the stress on the casing.The overlapping fracturing zones between wells cause significant stress interference,greatly increasing the risk of deformation and failure.By analyzing the Mises stress distribution of multi-well casings,high-risk locations for casing failure can be identified.The conclusion is that the key to preventing casing failure in platform wells in the “well factory” model is to optimize the spatial distribution of fracturing zones between wells and reasonably arrange well spacing.The study provides new insights and methods for predicting casing failure in unconventional oil and gas reservoirs and offers references for optimizing drilling and fracturing designs.展开更多
To resolve the issue of design for multi-stage and multi-cluster fracturing in multi-zone reservoirs, a new efficient algorithm for the planar 3 D multi-fracture propagation model was proposed. The model considers flu...To resolve the issue of design for multi-stage and multi-cluster fracturing in multi-zone reservoirs, a new efficient algorithm for the planar 3 D multi-fracture propagation model was proposed. The model considers fluid flow in the wellbore-perforation-fracture system and fluid leak-off into the rock matrix, and uses a 3 D boundary integral equation to describe the solid deformation. The solid-fluid coupling equation is solved by an explicit integration algorithm, and the fracture front is determined by the uniform tip asymptotic solutions and shortest path algorithm. The accuracy of the algorithm is verified by the analytical solution of radial fracture, results of the implicit level set algorithm, and results of organic glass fracturing experiment. Compared with the implicit level set algorithm(ILSA), the new algorithm is much higher in computation speed. The numerical case study is conducted based on a horizontal well in shale gas formation of Zhejiang oilfield. The impact of stress heterogeneity among multiple clusters and perforation number distribution on multi-fracture growth and fluid distribution among multiple fractures are analyzed by numerical simulation. The results show that reducing perforation number in each cluster can counteract the effect of stress contrast among perforation clusters. Adjusting perforation number in each cluster can promote uniform flux among clusters, and the perforation number difference should better be 1-2 among clusters. Increasing perforation number in the cluster with high in situ stress is conducive to uniform fluid partitioning. However, uniform fluid partitioning is not equivalent to uniform fracture geometry. The fracture geometry is controlled by the stress interference and horizontal principal stress profile jointly.展开更多
A novel laboratory simulation method for modeling multi-staged fracturing in a horizontal well was established based on a true tri-axial hydraulic fracturing simulation system. Using this method, the influences of net...A novel laboratory simulation method for modeling multi-staged fracturing in a horizontal well was established based on a true tri-axial hydraulic fracturing simulation system. Using this method, the influences of net pressure in hydraulic fracture, stage spacing, perforation parameter, horizontal stress bias and well cementation quality on the propagation geometry of multiple fractures in a tight sandstone formation were studied in detail. The specimen splitting and analogy analysis of fracturing curve patterns reveals: Multiple fractures tend to merge under the condition of high horizontal stress bias and short stage spacing with pre-existing hydraulic fractures under critical closure situation, and the propagation of subsequent fractures is possibly suppressed because of high net pressure in pre-created fractures and asymmetric distribution of fracture width. And the subsequently created fractures are situated in the induced stress decreasing zone due to long stage spacing, leading to weak stress interference, and perforation with intense density and deep penetration facilitates the decrease of initiation fracture pressure. The deflection angle of subsequent fracture and horizontal stress variation tend to be amplified under low horizontal bias with constant net pressure in fractures. The longitudinal fracture is likely to be initiated at the interface of wellbore and concrete sample with poor cementation quality. The initiation fracture pressure of the different stages increases in turn, with the largest increase of 30%. Pressure quickly declines after initiation with low propagation pressure when the transverse hydraulic fracture is formed. The pressure reduces with fluctuation after the initiation of fracture when the fracture deflects, the extension pressure is high, and the fracture formed is tortuous and narrow. There is a violently fluctuant rise of pressure with multiple peak values when longitudinal fracture created, and it is hard to distinguish the features between the initiation stage and propagation stage.展开更多
Particle-fluid two-phase flows in rock fractures and fracture networks play a pivotal role in determining the efficiency and effectiveness of hydraulic fracturing operations,a vital component in unconventional oil and...Particle-fluid two-phase flows in rock fractures and fracture networks play a pivotal role in determining the efficiency and effectiveness of hydraulic fracturing operations,a vital component in unconventional oil and gas extraction.Central to this phenomenon is the transport of proppants,tiny solid particles injected into the fractures to prevent them from closing once the injection is stopped.However,effective transport and deposition of proppant is critical in keeping fracture pathways open,especially in lowpermeability reservoirs.This review explores,then quantifies,the important role of fluid inertia and turbulent flows in governing proppant transport.While traditional models predominantly assume and then characterise flow as laminar,this may not accurately capture the complexities inherent in realworld hydraulic fracturing and proppant emplacement.Recent investigations highlight the paramount importance of fluid inertia,especially at the high Reynolds numbers typically associated with fracturing operations.Fluid inertia,often overlooked,introduces crucial forces that influence particle settling velocities,particle-particle interactions,and the eventual deposition of proppants within fractures.With their inherent eddies and transient and chaotic nature,turbulent flows introduce additional complexities to proppant transport,crucially altering proppant settling velocities and dispersion patterns.The following comprehensive survey of experimental,numerical,and analytical studies elucidates controls on the intricate dynamics of proppant transport under fluid inertia and turbulence-towards providing a holistic understanding of the current state-of-the-art,guiding future research directions,and optimising hydraulic fracturing practices.展开更多
Multi-stage hydraulic fracturing of horizontal wells is the main stimulation method in recovering gas from tight shale gas reservoirs, and stage spacing deter- mination is one of the key issues in fracturing design. T...Multi-stage hydraulic fracturing of horizontal wells is the main stimulation method in recovering gas from tight shale gas reservoirs, and stage spacing deter- mination is one of the key issues in fracturing design. The initiation and propagation of hydraulic fractures will cause stress redistribution and may activate natural fractures in the reservoir. Due to the limitation of the analytical method in calculation of induced stresses, we propose a numerical method, which incorporates the interaction of hydraulic fractures and the wellbore, and analyzes the stress distri- bution in the reservoir under different stage spacing. Simulation results indicate the following： （1） The induced stress was overestimated from the analytical method because it did not take into account the interaction between hydraulic fractures and the horizontal wellbore. （2） The hydraulic fracture had a considerable effect on the redis- tribution of stresses in the direction of the horizontal wellbore in the reservoir. The stress in the direction per- pendicular to the horizontal wellbore after hydraulic frac- turing had a minor change compared with the original in situ stress. （3） Stress interferences among fractures were greatly connected with the stage spacing and the distance from the wellbore. When the fracture length was 200 m, and the stage spacing was 50 m, the stress redistribution due to stage fracturing may divert the original stress pat- tern, which might activate natural fractures so as to generate a complex fracture network.展开更多
The application of machine learning for pyrite discrimination establishes a robust foundation for constructing the ore-forming history of multi-stage deposits;however,published models face challenges related to limite...The application of machine learning for pyrite discrimination establishes a robust foundation for constructing the ore-forming history of multi-stage deposits;however,published models face challenges related to limited,imbalanced datasets and oversampling.In this study,the dataset was expanded to approximately 500 samples for each type,including 508 sedimentary,573 orogenic gold,548 sedimentary exhalative(SEDEX)deposits,and 364 volcanogenic massive sulfides(VMS)pyrites,utilizing random forest(RF)and support vector machine(SVM)methodologies to enhance the reliability of the classifier models.The RF classifier achieved an overall accuracy of 99.8%,and the SVM classifier attained an overall accuracy of 100%.The model was evaluated by a five-fold cross-validation approach with 93.8%accuracy for the RF and 94.9%for the SVM classifier.These results demonstrate the strong feasibility of pyrite classification,supported by a relatively large,balanced dataset and high accuracy rates.The classifier was employed to reveal the genesis of the controversial Keketale Pb-Zn deposit in NW China,which has been inconclusive among SEDEX,VMS,or a SEDEX-VMS transition.Petrographic investigations indicated that the deposit comprises early fine-grained layered pyrite(Py1)and late recrystallized pyrite(Py2).The majority voting classified Py1 as the VMS type,with an accuracy of RF and SVM being 72.2%and 75%,respectively,and confirmed Py2 as an orogenic type with 74.3% and 77.1%accuracy,respectively.The new findings indicated that the Keketale deposit originated from a submarine VMS mineralization system,followed by late orogenic-type overprinting of metamorphism and deformation,which is consistent with the geological and geochemical observations.This study further emphasizes the advantages of Machine learning(ML)methods in accurately and directly discriminating the deposit types and reconstructing the formation history of multi-stage deposits.展开更多
Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout a...Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout analyses were based on continuum mechanics.In addition to failure in intact rocks,wellbore breakouts can also be initiated along natural discontinuities,e.g.weak planes and fractures.Furthermore,the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation.This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin,China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method(DEM)and the discrete fracture network(DFN).The developed model was validated against caliper log measurement,and its stability study was carried out by stress and displacement analyses.A parametric study was performed to investigate the effects of the characteristics of fracture distribution(orientation and length)on borehole stability by sensitivity studies.Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability.Moreover,an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress(i.e.the direction of wellbore breakouts)towards alternative directions,ultimately leading to the whole wellbore failure.These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.展开更多
Predicting the productivity of multistage fractured horizontal wells plays an important role in exploiting unconventional resources.In recent years,machine learning(ML)models have emerged as a new approach for such st...Predicting the productivity of multistage fractured horizontal wells plays an important role in exploiting unconventional resources.In recent years,machine learning(ML)models have emerged as a new approach for such studies.However,the scarcity of sufficient real data for model training often leads to imprecise predictions,even though the models trained with real data better characterize geological and engineering features.To tackle this issue,we propose an ML model that can obtain reliable results even with a small amount of data samples.Our model integrates the synthetic minority oversampling technique(SMOTE)to expand the data volume,the support vector machine(SVM)for model training,and the particle swarm optimization(PSO)algorithm for optimizing hyperparameters.To enhance the model performance,we conduct feature fusion and dimensionality reduction.Additionally,we examine the influences of different sample sizes and ML models for training.The proposed model demonstrates higher prediction accuracy and generalization ability,achieving a predicted R^(2)value of up to 0.9 for the test set,compared to the traditional ML techniques with an R^(2)of 0.13.This model accurately predicts the production of fractured horizontal wells even with limited samples,supplying an efficient tool for optimizing the production of unconventional resources.Importantly,the model holds the potential applicability to address similar challenges in other fields constrained by scarce data samples.展开更多
Multi-stage SRV fracturing in horizontal wells is a new technology developed at home and abroad in recent years to effectively develop shale gas or low-permeability reservoirs,but on the other hand makes the mechanica...Multi-stage SRV fracturing in horizontal wells is a new technology developed at home and abroad in recent years to effectively develop shale gas or low-permeability reservoirs,but on the other hand makes the mechanical environment of fracturing strings more complicated at the same time.In view of this,based on the loading features of tubing strings during the multi-stage fracturing of a horizontal well,mechanical models were established for three working cases of multiple packer setting,open differential-pressure sliding sleeve,and open ball-injection sliding sleeve under a hold-down packer.Moreover,mathematical models were respectively built for the above three cases.According to the Lame formula and Von Mises stress calculation formula for the thick-walled cylinder in the theory of elastic mechanics,a mathematical model was also established to calculate the equivalent stress for tubing string safety evaluation when the fracturing string was under the combined action of inner pressure,external squeezing force and axial stress,and another mathematical model was built for the mechanical strength and safety evaluation of multi-stage fracturing strings.In addition,a practical software was developed for the mechanical safety evaluation of horizontal well multi-stage fracturing strings according to the mathematical model developed for the mechanical calculation of the multi-packer string in horizontal wells.The research results were applied and verified in a gas well of Tahe Oilfield in the Tarim Basin with excellent effects,providing a theoretical basis and a simple and reliable technical means for optimal design and safety evaluation of safe operational parameters of multistage fracturing strings in horizontal wells.展开更多
The high proportion of uncertain distributed power sources and the access to large-scale random electric vehicle(EV)charging resources further aggravate the voltage fluctuation of the distribution network,and the exis...The high proportion of uncertain distributed power sources and the access to large-scale random electric vehicle(EV)charging resources further aggravate the voltage fluctuation of the distribution network,and the existing research has not deeply explored the EV active-reactive synergistic regulating characteristics,and failed to realize themulti-timescale synergistic control with other regulatingmeans,For this reason,this paper proposes amultilevel linkage coordinated optimization strategy to reduce the voltage deviation of the distribution network.Firstly,a capacitor bank reactive power compensation voltage control model and a distributed photovoltaic(PV)activereactive power regulationmodel are established.Additionally,an external characteristicmodel of EVactive-reactive power regulation is developed considering the four-quadrant operational characteristics of the EVcharger.Amultiobjective optimization model of the distribution network is then constructed considering the time-series coupling constraints of multiple types of voltage regulators.A multi-timescale control strategy is proposed by considering the impact of voltage regulators on active-reactive EV energy consumption and PV energy consumption.Then,a four-stage voltage control optimization strategy is proposed for various types of voltage regulators with multiple time scales.Themulti-objective optimization is solved with the improvedDrosophila algorithmto realize the power fluctuation control of the distribution network and themulti-stage voltage control optimization.Simulation results validate that the proposed voltage control optimization strategy achieves the coordinated control of decentralized voltage control resources in the distribution network.It effectively reduces the voltage deviation of the distribution network while ensuring the energy demand of EV users and enhancing the stability and economic efficiency of the distribution network.展开更多
Due to complex geological structures and a narrow safe mud density window,offshore fractured formations frequently encounter severe lost circulation(LC)during drilling,significantly hindering oil and gas exploration a...Due to complex geological structures and a narrow safe mud density window,offshore fractured formations frequently encounter severe lost circulation(LC)during drilling,significantly hindering oil and gas exploration and development.Predicting LC risks enables the targeted implementation of mitigation strategies,thereby reducing the frequency of such incidents.To address the limitations of existing 3D geomechanical modeling in predicting LC,such as arbitrary factor selection,subjective weight assignment,and the inability to achieve pre-drilling prediction along the entire well section,an improved prediction method is proposed.This method integrates multi-source data and incorporates three LC-related sensitivity factors:fracture characteristics,rock brittleness,and in-situ stress conditions.A quantitative risk assessment model for LC is developed by combining the subjective analytic hierarchy process with the objective entropy weight method(EWM)to assign weights.Subsequently,3D geomechanical modeling is applied to identify regional risk zones,enabling digital visualization for pre-drilling risk prediction.The developed 3D LC risk prediction model was validated using actual LC incidents from drilled wells.Results were generally consistent with field-identified LC zones,with an average relative error of 19.08%,confirming its reliability.This method provides practical guidance for mitigating potential LC risks and optimizing drilling program designs in fractured formations.展开更多
Grouting has been the most effective approach to mitigate water inrush disasters in underground engineering due to its ability to plug groundwater and enhance rock strength.Nevertheless,there is a lack of potent numer...Grouting has been the most effective approach to mitigate water inrush disasters in underground engineering due to its ability to plug groundwater and enhance rock strength.Nevertheless,there is a lack of potent numerical tools for assessing the grouting effectiveness in water-rich fractured strata.In this study,the hydro-mechanical coupled discontinuous deformation analysis(HM-DDA)is inaugurally extended to simulate the grouting process in a water-rich discrete fracture network(DFN),including the slurry migration,fracture dilation,water plugging in a seepage field,and joint reinforcement after coagulation.To validate the capabilities of the developed method,several numerical examples are conducted incorporating the Newtonian fluid and Bingham slurry.The simulation results closely align with the analytical solutions.Additionally,a set of compression tests is conducted on the fresh and grouted rock specimens to verify the reinforcement method and calibrate the rational properties of reinforced joints.An engineering-scale model based on a real water inrush case of the Yonglian tunnel in a water-rich fractured zone has been established.The model demonstrates the effectiveness of grouting reinforcement in mitigating water inrush disaster.The results indicate that increased grouting pressure greatly affects the regulation of water outflow from the tunnel face and the prevention of rock detachment face after excavation.展开更多
Injecting CO_(2)when the gas reservoir of tight sandstone is depleted can achieve the dual purposes of greenhouse gas storage and enhanced gas recovery(CS-EGR).To evaluate the feasibility of CO_(2)injection to enhance...Injecting CO_(2)when the gas reservoir of tight sandstone is depleted can achieve the dual purposes of greenhouse gas storage and enhanced gas recovery(CS-EGR).To evaluate the feasibility of CO_(2)injection to enhance gas recovery and understand the production mechanism,a numerical simulation model of CS-EGR in multi-stage fracturing horizontal wells is established.The behavior of gas production and CO_(2)sequestration is then analyzed through numerical simulation,and the impact of fracture parameters on production performance is examined.Simulation results show that the production rate increases significantly and a large amount of CO_(2)is stored in the reservoir,proving the technical potential.However,hydraulic fractures accelerate CO_(2)breakthrough,resulting in lower gas recovery and lower CO_(2)storage than in gas reservoirs without fracturing.Increasing the length of hydraulic fractures can significantly increase CH4production,but CO_(2)breakthrough will advance.Staggered and spaced perforation of hydraulic fractures in injection wells and production wells changes the fluid flow path,which can delay CO_(2)breakthrough and benefit production efficiency.The fracture network of massive hydraulic fracturing has a positive effect on the CS-EGR.As a result,CH4production,gas recovery,and CO_(2)storage increase with the increase in stimulated reservoir volume.展开更多
BZ oilfield in Bohai Bay of China was a typical offshore low permeability oilfield, which was restricted by many factors such as environment and economy. In this paper, the development characteristics of BZ oilfield w...BZ oilfield in Bohai Bay of China was a typical offshore low permeability oilfield, which was restricted by many factors such as environment and economy. In this paper, the development characteristics of BZ oilfield were summarized in depth, and the new development mode of offshore low-permeability oilfield was explored from reservoir prediction, well spacing and fracturing technology. Taking BZ oilfield as an example, a set of technical system for the effective development of offshore low permeability oilfield had been formed through research, which mainly includes reservoir prediction and evaluation of offshore middle and deep low permeability oilfield, optimization of horizontal well pattern, multi-stage fracturing design of horizontal well and other technologies. The results show that improving the resolution of seismic data, strengthening the analysis of seismic reflection characteristics and carrying out the comprehensive study of seismic geology were the keys to solve the reservoir prediction of offshore low-permeability oil fields. Multi-stage fracturing horizontal well pattern is the main pattern of offshore low-permeability oilfield development. The parameters of multi-stage fracturing horizontal well together affect the development effect. Selecting the optimal fractured horizontal well pattern can greatly improve the development effect. The successful combination and application of new technology system was the foundation and core of conquering offshore low-permeability oil fields. On the basis of understanding the geological characteristics of oil reservoirs, it is an effective means of developing offshore low-permeability oil fields by selecting reasonable production methods, well types and well patterns. Using efficient perforation and fracturing technology to successfully control fracture parameters and form optimal injection and production well pattern was the key to improve low permeability offshore oil fields.展开更多
Coupled thermo-hydro-mechanical(THM)processes in fractured rock are playing a crucial role in geoscience and geoengineering applications.Diverse and conceptually distinct approaches have emerged over the past decades ...Coupled thermo-hydro-mechanical(THM)processes in fractured rock are playing a crucial role in geoscience and geoengineering applications.Diverse and conceptually distinct approaches have emerged over the past decades in both continuum and discontinuum perspectives leading to significant progress in their comprehending and modeling.This review paper offers an integrated perspective on existing modeling methodologies providing guidance for model selection based on the initial and boundary conditions.By comparing various models,one can better assess the uncertainties in predictions,particularly those related to the conceptual models.The review explores how these methodologies have significantlyenhanced the fundamental understanding of how fractures respond to fluid injection and production,and improved predictive capabilities pertaining to coupled processes within fractured systems.It emphasizes the importance of utilizing advanced computational technologies and thoroughly considering fundamental theories and principles established through past experimental evidence and practical experience.The selection and calibration of model parameters should be based on typical ranges and applied to the specificconditions of applications.The challenges arising from inherent heterogeneity and uncertainties,nonlinear THM coupled processes,scale dependence,and computational limitations in representing fieldscale fractures are discussed.Realizing potential advances on computational capacity calls for methodical conceptualization,mathematical modeling,selection of numerical solution strategies,implementation,and calibration to foster simulation outcomes that intricately reflectthe nuanced complexities of geological phenomena.Future research efforts should focus on innovative approaches to tackle the hurdles and advance the state-of-the-art in this critical fieldof study.展开更多
As the dominant seepage channel in rock masses,it is of great significance to study the influence of fracture roughness distribution on seepage and heat transfer in rock masses.In this paper,the fracture roughness dis...As the dominant seepage channel in rock masses,it is of great significance to study the influence of fracture roughness distribution on seepage and heat transfer in rock masses.In this paper,the fracture roughness distribution functions of the Bakhtiary dam site and Oskarshamn/Forsmark mountain were fitted using statistical methods.The COMSOL Multiphysics finite element software was utilized to analyze the effects of fracture roughness distribution types and empirical formulas for fracture hydraulic aperture on the seepage field and temperature field of rock masses.The results show that:(1)The fracture roughness at the Bakhtiary dam site and Oskarshamn/Forsmark mountain follows lognormal and normal distributions,respectively;(2)For rock masses with the same expected value and standard deviation of fracture roughness,the outflow from rock masses with lognormal distribution of fracture roughness is significantly larger than that of rock masses with normal distribution of fracture roughness;(3)The fracture hydraulic aperture,outflow,and cold front distance of the Li and Jiang model are significantly larger than those of the Barton model;(4)The outflow,hydraulic pressure distribution,and temperature distribution of the Barton model are more sensitive to the fracture roughness distribution type than those of the Li and Jiang model.展开更多
Deep underground excavation causes considerable unloading effects,leading to a pronounced bias pressure phenomenon.The deformation and seepage characteristics of rock masses under different gas and confining pressures...Deep underground excavation causes considerable unloading effects,leading to a pronounced bias pressure phenomenon.The deformation and seepage characteristics of rock masses under different gas and confining pressures were investigated via triaxial loading and unloading seepage tests.When the influential coefficient of effective confining pressure(β)is less than 0.065,the seepage force considerably weakens the strength of fractured rock masses.Conversely,whenβis greater than 0.065,the opposite is true.Moreover,the increase in the axial load leads to an increase in the precast fracture volumetric strain,which is the main reason for the increase in fracture permeability.This effect is particularly significant during the unloading stage.Based on the test results,a method for calculating the dynamic seepage evolution of rock masses,considering the effects of rock mass damage and fracture deformation,is introduced,and the effectiveness of the calculation is validated.The entire description of the seepage under loading and unloading was accomplished.The equivalent relationship between the lateral and normal stresses on fracture surfaces ranges from 0.001 to 0.1,showing an exponential variation between the lateral stress influence coefficient on normal deformation(χ)and seepage pressure.Before the failure of the rock mass,the seepage in the fractures was in a linear laminar flow state.However,after the failure,when the gas pressure reached 2 MPa,the flow state in the fractures transitioned to nonlinear laminar flow.The results are important for predicting hazardous gas leaks during deep underground engineering excavation.展开更多
Although previous researchers have attempted to decipher ore genesis and mineralization in the Erdaokan Ag-Pb-Zn deposit,some uncertainties regarding the mineralization process and evolution of both ore-forming fluids...Although previous researchers have attempted to decipher ore genesis and mineralization in the Erdaokan Ag-Pb-Zn deposit,some uncertainties regarding the mineralization process and evolution of both ore-forming fluids and magnetite types still need to be addressed.In this study,we obtained new EPMA,LA-ICP-MS,and in situ Fe isotope data from magnetite from the Erdaokan deposit,in order to better understand the mineralization mechanism and evolution of both magnetite and the ore-forming fluids.Our results identified seven types of magnetite at Erdaokan:disseminated magnetite(Mag1),coarse-grained magnetite(Mag2a),radial magnetite(Mag2b),fragmented fine-grained magnetite(Mag2c),vermicular gel magnetite(Mag3a1 and Mag3a2),colloidal magnetite(Mag3b)and dark gray magnetite(Mag4).All of the magnetite types were hydrothermal in origin and generally low in Ti(<400 ppm)and Ni(<800 ppm),while being enriched in light Fe isotopes(δ^(56)Fe ranging from−1.54‰to−0.06‰).However,they exhibit different geochemical signatures and are thus classified into high-manganese magnetite(Mag1,MnO>5 wt%),low-silicon magnetite(Mag2a-c,SiO_(2)<1 wt%),high-silicon magnetite(Mag3a-b,SiO_(2)from 1 to 7 wt%)and high-silicon-manganese magnetite(Mag4,SiO_(2)>1 wt%,MnO>0.2 wt%),each being formed within distinct hydrothermal environments.Based on mineralogy,elemental geochemistry,Fe isotopes,temperature trends,TMg-mag and(Ti+V)vs.(Al+Mn)diagrams,we propose that the Erdaokan Ag-Pb-Zn deposit underwent multi-stage mineralization,which can be broken down into four stages and nine sub-stages.Mag1,Mag2a-c,Mag3a-b and Mag4 were formed during the first sub-stage of each of the four stages,respectively.Additionally,fluid mixing,cooling and depressurization boiling were identified as the main mechanisms for mineral precipitation.The enrichment of Ag was significantly enhanced by the superposition of multi-stage ore-forming hydrothermal fluids in the Erdaokan Ag-Pb-Zn deposit.展开更多
基金supported by the National Natural Science Foundation of China (No. 51234007, No. 51490654, No. 51504276, and No. 51504277)Program for Changjiang Scholars and Innovative Research Team in University (IRT1294)+3 种基金the Natural Science Foundation of Shandong Province (ZR2014EL016, ZR2014EEP018)China Postdoctoral Science Foundation (No. 2014M551989 and No. 2015T80762)the Major Programs of Ministry of Education of China (No. 311009)Introducing Talents of Discipline to Universities (B08028)
文摘In order to investigate the influence on shale gas well productivity caused by gas transport in nanometer- size pores, a mathematical model of multi-stage fractured horizontal wells in shale gas reservoirs is built, which considers the influence of viscous flow, Knudsen diffusion, surface diffusion, and adsorption layer thickness. A dis- crete-fracture model is used to simplify the fracture mod- cling, and a finite element method is applied to solve the model. The numerical simulation results indicate that with a decrease in the intrinsic matrix permeability, Knudsen diffusion and surface diffusion contributions to production become large and cannot be ignored. The existence of an adsorption layer on the nanopore surfaces reduces the effective pore radius and the effective porosity, resulting in low production from fractured horizontal wells. With a decrease in the pore radius, considering the adsorption layer, the production reduction rate increases. When the pore radius is less than 10 nm, because of the combined impacts of Knudsen diffusion, surface diffusion, and adsorption layers, the production of multi-stage fractured horizontal wells increases with a decrease in the pore pressure. When the pore pressure is lower than 30 MPa, the rate of production increase becomes larger with a decrease in pore pressure.
文摘The Gabes aquifer system,located in southeastern Tunisia,is a crucial resource for supporting local socio-economic activities.Due to its dual porosity structure,is particularly vulnerable to pollution.This study aims to develop a hybrid model that combines the Fracture Aquifer Index(FAI)with the conventional GOD(Groundwater occurrence,Overall lithology,Depth to water table)method,to assess groundwater vulnerability in fractured aquifer.To develop the hybrid model,the classical GOD method was integrated with FAI to produce a single composite index.Each parameter within both GOD and FAI was scored,and a final index was calculated to delineate vulnerable areas.The results show that the study area can be classified into four vulnerability levels:Very low,low,moderate,and high,indicating that approximately 8%of the area exhibits very low vulnerability,29%has low vulnerability,25%falls into the moderate category,and 38%is considered highly vulnerable.The FAI-GOD model further incorporates fracture network characteristics.This refinement reduces the classification to three vulnerability classes:Low,medium,and high.The outcomes demonstrate that 46%of the area is highly vulnerable due to a dense concentration of fractures,while 17%represents an intermediate zone characterized by either shallow or deeper fractures.In contrast,37%corresponds to areas with lightly fractured rock,where the impact on vulnerability is minimal.Multivariate statistical analysis was employed using Principal Components Analysis(PCA)and Hierarchical Cluster Analysis(HCA)on 24 samples across six variables.The first three components account for over 76%of the total variance,reinforcing the significance of fracture dynamics in classifying vulnerability levels.The FAI-GOD model removes the very-low-vulnerability class and expands the spatial extent of low-and high-vulnerability zones,reflecting the dominant influence of fracture networks on aquifer sensitivity.While both indices use a five-class system,FAI-GOD redistributes vulnerability by eliminating very-low-vulnerability areas and amplifying low/high categories,highlighting the critical role of fractures.A strong correlation(R2=0.94)between the GOD and FAI-GOD indices,demonstrated through second-order polynomial regression,confirms the robustness of the FAI-GOD model in accurately predicting vulnerability to pollution.This model provides a useful framework for assessing the vulnerability of complex aquifers and serves as a decision-making tool for groundwater managers in similar areas.
基金supported by the National Natural Science Foundation of China (No.52104008&No.52274042)the Natural Science Foundation of Sichuan,China (No.2024NSFSC0963)。
文摘The “well factory” mode's high-density well placement and multi-stage hydraulic fracturing technology enable efficient development of unconventional oil and gas resources.However,the deployment of platform wells in the “well factory” model results in small wellbore spacing,and the stress disturbances caused by fracturing operations may affect neighboring wells,leading to inter-well interference phenomena that cause casing deformation.This study investigates the issue of inter-well interference causing casing deformation or even failure during multi-stage hydraulic fracturing in the “well factory”model,and predicts high-risk locations for casing failure.A flow-mechanics coupled geomechanical finite element model with retaining geological stratification characteristics was established.Based on the theory of hydraulic fracturing-induced rock fragmentation and fluid action leading to the degradation of rock mechanical properties,the model simulated the four-dimensional evolution of multi-well fracturing areas over time and space,calculating the disturbance in the regional stress field caused by fracturing operations.Subsequently,the stress distribution of multiple well casings at different time points was calculated to predict high-risk locations for casing failure.The research results show that the redistribution of the stress field in the fracturing area increases the stress on the casing.The overlapping fracturing zones between wells cause significant stress interference,greatly increasing the risk of deformation and failure.By analyzing the Mises stress distribution of multi-well casings,high-risk locations for casing failure can be identified.The conclusion is that the key to preventing casing failure in platform wells in the “well factory” model is to optimize the spatial distribution of fracturing zones between wells and reasonably arrange well spacing.The study provides new insights and methods for predicting casing failure in unconventional oil and gas reservoirs and offers references for optimizing drilling and fracturing designs.
基金China National Science and Technology Major Project(2016ZX05023).
文摘To resolve the issue of design for multi-stage and multi-cluster fracturing in multi-zone reservoirs, a new efficient algorithm for the planar 3 D multi-fracture propagation model was proposed. The model considers fluid flow in the wellbore-perforation-fracture system and fluid leak-off into the rock matrix, and uses a 3 D boundary integral equation to describe the solid deformation. The solid-fluid coupling equation is solved by an explicit integration algorithm, and the fracture front is determined by the uniform tip asymptotic solutions and shortest path algorithm. The accuracy of the algorithm is verified by the analytical solution of radial fracture, results of the implicit level set algorithm, and results of organic glass fracturing experiment. Compared with the implicit level set algorithm(ILSA), the new algorithm is much higher in computation speed. The numerical case study is conducted based on a horizontal well in shale gas formation of Zhejiang oilfield. The impact of stress heterogeneity among multiple clusters and perforation number distribution on multi-fracture growth and fluid distribution among multiple fractures are analyzed by numerical simulation. The results show that reducing perforation number in each cluster can counteract the effect of stress contrast among perforation clusters. Adjusting perforation number in each cluster can promote uniform flux among clusters, and the perforation number difference should better be 1-2 among clusters. Increasing perforation number in the cluster with high in situ stress is conducive to uniform fluid partitioning. However, uniform fluid partitioning is not equivalent to uniform fracture geometry. The fracture geometry is controlled by the stress interference and horizontal principal stress profile jointly.
文摘A novel laboratory simulation method for modeling multi-staged fracturing in a horizontal well was established based on a true tri-axial hydraulic fracturing simulation system. Using this method, the influences of net pressure in hydraulic fracture, stage spacing, perforation parameter, horizontal stress bias and well cementation quality on the propagation geometry of multiple fractures in a tight sandstone formation were studied in detail. The specimen splitting and analogy analysis of fracturing curve patterns reveals: Multiple fractures tend to merge under the condition of high horizontal stress bias and short stage spacing with pre-existing hydraulic fractures under critical closure situation, and the propagation of subsequent fractures is possibly suppressed because of high net pressure in pre-created fractures and asymmetric distribution of fracture width. And the subsequently created fractures are situated in the induced stress decreasing zone due to long stage spacing, leading to weak stress interference, and perforation with intense density and deep penetration facilitates the decrease of initiation fracture pressure. The deflection angle of subsequent fracture and horizontal stress variation tend to be amplified under low horizontal bias with constant net pressure in fractures. The longitudinal fracture is likely to be initiated at the interface of wellbore and concrete sample with poor cementation quality. The initiation fracture pressure of the different stages increases in turn, with the largest increase of 30%. Pressure quickly declines after initiation with low propagation pressure when the transverse hydraulic fracture is formed. The pressure reduces with fluctuation after the initiation of fracture when the fracture deflects, the extension pressure is high, and the fracture formed is tortuous and narrow. There is a violently fluctuant rise of pressure with multiple peak values when longitudinal fracture created, and it is hard to distinguish the features between the initiation stage and propagation stage.
基金the Australian Research Council Discovery Project(ARC DP 220100851)scheme and would acknowledge that.
文摘Particle-fluid two-phase flows in rock fractures and fracture networks play a pivotal role in determining the efficiency and effectiveness of hydraulic fracturing operations,a vital component in unconventional oil and gas extraction.Central to this phenomenon is the transport of proppants,tiny solid particles injected into the fractures to prevent them from closing once the injection is stopped.However,effective transport and deposition of proppant is critical in keeping fracture pathways open,especially in lowpermeability reservoirs.This review explores,then quantifies,the important role of fluid inertia and turbulent flows in governing proppant transport.While traditional models predominantly assume and then characterise flow as laminar,this may not accurately capture the complexities inherent in realworld hydraulic fracturing and proppant emplacement.Recent investigations highlight the paramount importance of fluid inertia,especially at the high Reynolds numbers typically associated with fracturing operations.Fluid inertia,often overlooked,introduces crucial forces that influence particle settling velocities,particle-particle interactions,and the eventual deposition of proppants within fractures.With their inherent eddies and transient and chaotic nature,turbulent flows introduce additional complexities to proppant transport,crucially altering proppant settling velocities and dispersion patterns.The following comprehensive survey of experimental,numerical,and analytical studies elucidates controls on the intricate dynamics of proppant transport under fluid inertia and turbulence-towards providing a holistic understanding of the current state-of-the-art,guiding future research directions,and optimising hydraulic fracturing practices.
基金supported by the Natural Science Foundation of China (Grant No. 51490653, Basic Theoretical Research of Shale Oil and Gas Effective Development)
文摘Multi-stage hydraulic fracturing of horizontal wells is the main stimulation method in recovering gas from tight shale gas reservoirs, and stage spacing deter- mination is one of the key issues in fracturing design. The initiation and propagation of hydraulic fractures will cause stress redistribution and may activate natural fractures in the reservoir. Due to the limitation of the analytical method in calculation of induced stresses, we propose a numerical method, which incorporates the interaction of hydraulic fractures and the wellbore, and analyzes the stress distri- bution in the reservoir under different stage spacing. Simulation results indicate the following： （1） The induced stress was overestimated from the analytical method because it did not take into account the interaction between hydraulic fractures and the horizontal wellbore. （2） The hydraulic fracture had a considerable effect on the redis- tribution of stresses in the direction of the horizontal wellbore in the reservoir. The stress in the direction per- pendicular to the horizontal wellbore after hydraulic frac- turing had a minor change compared with the original in situ stress. （3） Stress interferences among fractures were greatly connected with the stage spacing and the distance from the wellbore. When the fracture length was 200 m, and the stage spacing was 50 m, the stress redistribution due to stage fracturing may divert the original stress pat- tern, which might activate natural fractures so as to generate a complex fracture network.
基金the National Key Research and Development Program of China(2021YFC2900300)the Natural Science Foundation of Guangdong Province(2024A1515030216)+2 种基金MOST Special Fund from State Key Laboratory of Geological Processes and Mineral Resources,China University of Geosciences(GPMR202437)the Guangdong Province Introduced of Innovative R&D Team(2021ZT09H399)the Third Xinjiang Scientific Expedition Program(2022xjkk1301).
文摘The application of machine learning for pyrite discrimination establishes a robust foundation for constructing the ore-forming history of multi-stage deposits;however,published models face challenges related to limited,imbalanced datasets and oversampling.In this study,the dataset was expanded to approximately 500 samples for each type,including 508 sedimentary,573 orogenic gold,548 sedimentary exhalative(SEDEX)deposits,and 364 volcanogenic massive sulfides(VMS)pyrites,utilizing random forest(RF)and support vector machine(SVM)methodologies to enhance the reliability of the classifier models.The RF classifier achieved an overall accuracy of 99.8%,and the SVM classifier attained an overall accuracy of 100%.The model was evaluated by a five-fold cross-validation approach with 93.8%accuracy for the RF and 94.9%for the SVM classifier.These results demonstrate the strong feasibility of pyrite classification,supported by a relatively large,balanced dataset and high accuracy rates.The classifier was employed to reveal the genesis of the controversial Keketale Pb-Zn deposit in NW China,which has been inconclusive among SEDEX,VMS,or a SEDEX-VMS transition.Petrographic investigations indicated that the deposit comprises early fine-grained layered pyrite(Py1)and late recrystallized pyrite(Py2).The majority voting classified Py1 as the VMS type,with an accuracy of RF and SVM being 72.2%and 75%,respectively,and confirmed Py2 as an orogenic type with 74.3% and 77.1%accuracy,respectively.The new findings indicated that the Keketale deposit originated from a submarine VMS mineralization system,followed by late orogenic-type overprinting of metamorphism and deformation,which is consistent with the geological and geochemical observations.This study further emphasizes the advantages of Machine learning(ML)methods in accurately and directly discriminating the deposit types and reconstructing the formation history of multi-stage deposits.
基金supported by National Natural Science Foundation of China(Grant Nos.52074312 and 52211530097)CNPC Science and Technology Innovation Foundation(Grant No.2021DQ02-0505).
文摘Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout analyses were based on continuum mechanics.In addition to failure in intact rocks,wellbore breakouts can also be initiated along natural discontinuities,e.g.weak planes and fractures.Furthermore,the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation.This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin,China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method(DEM)and the discrete fracture network(DFN).The developed model was validated against caliper log measurement,and its stability study was carried out by stress and displacement analyses.A parametric study was performed to investigate the effects of the characteristics of fracture distribution(orientation and length)on borehole stability by sensitivity studies.Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability.Moreover,an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress(i.e.the direction of wellbore breakouts)towards alternative directions,ultimately leading to the whole wellbore failure.These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.
基金supported by the National Natural Science Foundation of China(52274055)the Shandong Provincial Natural Science Foundation(ZR2022YQ50)the Taishan Scholar Program of Shandong Province(tsqn202408088)。
文摘Predicting the productivity of multistage fractured horizontal wells plays an important role in exploiting unconventional resources.In recent years,machine learning(ML)models have emerged as a new approach for such studies.However,the scarcity of sufficient real data for model training often leads to imprecise predictions,even though the models trained with real data better characterize geological and engineering features.To tackle this issue,we propose an ML model that can obtain reliable results even with a small amount of data samples.Our model integrates the synthetic minority oversampling technique(SMOTE)to expand the data volume,the support vector machine(SVM)for model training,and the particle swarm optimization(PSO)algorithm for optimizing hyperparameters.To enhance the model performance,we conduct feature fusion and dimensionality reduction.Additionally,we examine the influences of different sample sizes and ML models for training.The proposed model demonstrates higher prediction accuracy and generalization ability,achieving a predicted R^(2)value of up to 0.9 for the test set,compared to the traditional ML techniques with an R^(2)of 0.13.This model accurately predicts the production of fractured horizontal wells even with limited samples,supplying an efficient tool for optimizing the production of unconventional resources.Importantly,the model holds the potential applicability to address similar challenges in other fields constrained by scarce data samples.
基金Foundation for Doctors from the Chinese Ministry of Education(No.20135121110005)National Major Science&Technology Project“Key technology for well bores in marine carbonate reservoirs”(No.2011ZX05005-006).
文摘Multi-stage SRV fracturing in horizontal wells is a new technology developed at home and abroad in recent years to effectively develop shale gas or low-permeability reservoirs,but on the other hand makes the mechanical environment of fracturing strings more complicated at the same time.In view of this,based on the loading features of tubing strings during the multi-stage fracturing of a horizontal well,mechanical models were established for three working cases of multiple packer setting,open differential-pressure sliding sleeve,and open ball-injection sliding sleeve under a hold-down packer.Moreover,mathematical models were respectively built for the above three cases.According to the Lame formula and Von Mises stress calculation formula for the thick-walled cylinder in the theory of elastic mechanics,a mathematical model was also established to calculate the equivalent stress for tubing string safety evaluation when the fracturing string was under the combined action of inner pressure,external squeezing force and axial stress,and another mathematical model was built for the mechanical strength and safety evaluation of multi-stage fracturing strings.In addition,a practical software was developed for the mechanical safety evaluation of horizontal well multi-stage fracturing strings according to the mathematical model developed for the mechanical calculation of the multi-packer string in horizontal wells.The research results were applied and verified in a gas well of Tahe Oilfield in the Tarim Basin with excellent effects,providing a theoretical basis and a simple and reliable technical means for optimal design and safety evaluation of safe operational parameters of multistage fracturing strings in horizontal wells.
基金funded by the State Grid Corporation Science and Technology Project(5108-202218280A-2-391-XG).
文摘The high proportion of uncertain distributed power sources and the access to large-scale random electric vehicle(EV)charging resources further aggravate the voltage fluctuation of the distribution network,and the existing research has not deeply explored the EV active-reactive synergistic regulating characteristics,and failed to realize themulti-timescale synergistic control with other regulatingmeans,For this reason,this paper proposes amultilevel linkage coordinated optimization strategy to reduce the voltage deviation of the distribution network.Firstly,a capacitor bank reactive power compensation voltage control model and a distributed photovoltaic(PV)activereactive power regulationmodel are established.Additionally,an external characteristicmodel of EVactive-reactive power regulation is developed considering the four-quadrant operational characteristics of the EVcharger.Amultiobjective optimization model of the distribution network is then constructed considering the time-series coupling constraints of multiple types of voltage regulators.A multi-timescale control strategy is proposed by considering the impact of voltage regulators on active-reactive EV energy consumption and PV energy consumption.Then,a four-stage voltage control optimization strategy is proposed for various types of voltage regulators with multiple time scales.Themulti-objective optimization is solved with the improvedDrosophila algorithmto realize the power fluctuation control of the distribution network and themulti-stage voltage control optimization.Simulation results validate that the proposed voltage control optimization strategy achieves the coordinated control of decentralized voltage control resources in the distribution network.It effectively reduces the voltage deviation of the distribution network while ensuring the energy demand of EV users and enhancing the stability and economic efficiency of the distribution network.
基金supported by the National Natural Science Foundation of China(No.52074312)the CNPC Science and Technology Innovation Foundation(No.2021DQ02-0505)+1 种基金the Open Fund Project of the National Key Laboratory for the Enrichment Mechanism and Efficient Development of Shale Oil and Gas(No.36650000-24-ZC0609-0006)the Major Science and Technology Project of Karamay City(No.20232023zdzx0003).
文摘Due to complex geological structures and a narrow safe mud density window,offshore fractured formations frequently encounter severe lost circulation(LC)during drilling,significantly hindering oil and gas exploration and development.Predicting LC risks enables the targeted implementation of mitigation strategies,thereby reducing the frequency of such incidents.To address the limitations of existing 3D geomechanical modeling in predicting LC,such as arbitrary factor selection,subjective weight assignment,and the inability to achieve pre-drilling prediction along the entire well section,an improved prediction method is proposed.This method integrates multi-source data and incorporates three LC-related sensitivity factors:fracture characteristics,rock brittleness,and in-situ stress conditions.A quantitative risk assessment model for LC is developed by combining the subjective analytic hierarchy process with the objective entropy weight method(EWM)to assign weights.Subsequently,3D geomechanical modeling is applied to identify regional risk zones,enabling digital visualization for pre-drilling risk prediction.The developed 3D LC risk prediction model was validated using actual LC incidents from drilled wells.Results were generally consistent with field-identified LC zones,with an average relative error of 19.08%,confirming its reliability.This method provides practical guidance for mitigating potential LC risks and optimizing drilling program designs in fractured formations.
基金supported by the China Scholarship Council(CSC,Grant No.202108050072)JSPS KAKENHI(Grant No.JP19KK0121)。
文摘Grouting has been the most effective approach to mitigate water inrush disasters in underground engineering due to its ability to plug groundwater and enhance rock strength.Nevertheless,there is a lack of potent numerical tools for assessing the grouting effectiveness in water-rich fractured strata.In this study,the hydro-mechanical coupled discontinuous deformation analysis(HM-DDA)is inaugurally extended to simulate the grouting process in a water-rich discrete fracture network(DFN),including the slurry migration,fracture dilation,water plugging in a seepage field,and joint reinforcement after coagulation.To validate the capabilities of the developed method,several numerical examples are conducted incorporating the Newtonian fluid and Bingham slurry.The simulation results closely align with the analytical solutions.Additionally,a set of compression tests is conducted on the fresh and grouted rock specimens to verify the reinforcement method and calibrate the rational properties of reinforced joints.An engineering-scale model based on a real water inrush case of the Yonglian tunnel in a water-rich fractured zone has been established.The model demonstrates the effectiveness of grouting reinforcement in mitigating water inrush disaster.The results indicate that increased grouting pressure greatly affects the regulation of water outflow from the tunnel face and the prevention of rock detachment face after excavation.
基金the financial support of the National Natural Science Foundation of China(Grant No.52304018)China Postdoctoral Science Foundation(Grant No.2023TQ0014,Grant No.2023M730088)。
文摘Injecting CO_(2)when the gas reservoir of tight sandstone is depleted can achieve the dual purposes of greenhouse gas storage and enhanced gas recovery(CS-EGR).To evaluate the feasibility of CO_(2)injection to enhance gas recovery and understand the production mechanism,a numerical simulation model of CS-EGR in multi-stage fracturing horizontal wells is established.The behavior of gas production and CO_(2)sequestration is then analyzed through numerical simulation,and the impact of fracture parameters on production performance is examined.Simulation results show that the production rate increases significantly and a large amount of CO_(2)is stored in the reservoir,proving the technical potential.However,hydraulic fractures accelerate CO_(2)breakthrough,resulting in lower gas recovery and lower CO_(2)storage than in gas reservoirs without fracturing.Increasing the length of hydraulic fractures can significantly increase CH4production,but CO_(2)breakthrough will advance.Staggered and spaced perforation of hydraulic fractures in injection wells and production wells changes the fluid flow path,which can delay CO_(2)breakthrough and benefit production efficiency.The fracture network of massive hydraulic fracturing has a positive effect on the CS-EGR.As a result,CH4production,gas recovery,and CO_(2)storage increase with the increase in stimulated reservoir volume.
文摘BZ oilfield in Bohai Bay of China was a typical offshore low permeability oilfield, which was restricted by many factors such as environment and economy. In this paper, the development characteristics of BZ oilfield were summarized in depth, and the new development mode of offshore low-permeability oilfield was explored from reservoir prediction, well spacing and fracturing technology. Taking BZ oilfield as an example, a set of technical system for the effective development of offshore low permeability oilfield had been formed through research, which mainly includes reservoir prediction and evaluation of offshore middle and deep low permeability oilfield, optimization of horizontal well pattern, multi-stage fracturing design of horizontal well and other technologies. The results show that improving the resolution of seismic data, strengthening the analysis of seismic reflection characteristics and carrying out the comprehensive study of seismic geology were the keys to solve the reservoir prediction of offshore low-permeability oil fields. Multi-stage fracturing horizontal well pattern is the main pattern of offshore low-permeability oilfield development. The parameters of multi-stage fracturing horizontal well together affect the development effect. Selecting the optimal fractured horizontal well pattern can greatly improve the development effect. The successful combination and application of new technology system was the foundation and core of conquering offshore low-permeability oil fields. On the basis of understanding the geological characteristics of oil reservoirs, it is an effective means of developing offshore low-permeability oil fields by selecting reasonable production methods, well types and well patterns. Using efficient perforation and fracturing technology to successfully control fracture parameters and form optimal injection and production well pattern was the key to improve low permeability offshore oil fields.
基金funding from the European Research Council(ERC)under the European Union’s Horizon 2020 Research and Innovation Program through the Starting Grant GEoREST(grant agreement No.801809)support by MICIU/AEI/10.13039/501100011033 and by"European Union Next Generation EU/PRTR"through the‘Ramón y Cajal’fellowship(reference RYC2021-032780-I)+9 种基金funding by MICIU/AEI/10.13039/501100011033 and by“ERDF,EU”through the‘HydroPoreII’project(reference PID2022-137652NBC44)support by the Institute for Korea Spent Nuclear Fuel(iKSNF)National Research Foundation of Korea(NRF)grant funded by the Korea government(Ministry of Science and ICT,MSIT)(2021M2E1A1085196)support by the Swedish Radiation Safety(SSM),Swedish Transport Administration(Trafikverket),Swedish Rock Engineering Foundation(BeFo),and Nordic Energy Research(Grant 187658)supported by the US Department of Energy(DOE),the Officeof Nuclear Energy,Spent Fuel and Waste Science and Technology Campaign,and by the US Department of Energy(DOE),the Office of Basic Energy Sciences,Chemical Sciences,Geosciences,and Biosciences Division both under Contract Number DE-AC02-05CH11231 with Lawrence Berkeley National Laboratorysupport from the US National Science Foundation(grant CMMI-2239630)funding from the European Research Council(ERC)under the European Union’s Horizon 2020 research and innovation programme(grant agreement No.101002507)the UK Natural Environment Research Council(NERC)for funding SeisGreen Project(Grant No.NE/W009293/1)which supported this workthe Royal Society UK for supporting this research through fellowship UF160443IMEDEA is an accredited"Maria de Maeztu Excellence Unit"(Grant CEX2021-001198,funded by MICIU/AEI/10.13039/501100011033).
文摘Coupled thermo-hydro-mechanical(THM)processes in fractured rock are playing a crucial role in geoscience and geoengineering applications.Diverse and conceptually distinct approaches have emerged over the past decades in both continuum and discontinuum perspectives leading to significant progress in their comprehending and modeling.This review paper offers an integrated perspective on existing modeling methodologies providing guidance for model selection based on the initial and boundary conditions.By comparing various models,one can better assess the uncertainties in predictions,particularly those related to the conceptual models.The review explores how these methodologies have significantlyenhanced the fundamental understanding of how fractures respond to fluid injection and production,and improved predictive capabilities pertaining to coupled processes within fractured systems.It emphasizes the importance of utilizing advanced computational technologies and thoroughly considering fundamental theories and principles established through past experimental evidence and practical experience.The selection and calibration of model parameters should be based on typical ranges and applied to the specificconditions of applications.The challenges arising from inherent heterogeneity and uncertainties,nonlinear THM coupled processes,scale dependence,and computational limitations in representing fieldscale fractures are discussed.Realizing potential advances on computational capacity calls for methodical conceptualization,mathematical modeling,selection of numerical solution strategies,implementation,and calibration to foster simulation outcomes that intricately reflectthe nuanced complexities of geological phenomena.Future research efforts should focus on innovative approaches to tackle the hurdles and advance the state-of-the-art in this critical fieldof study.
基金College Students Innovation and Entrepreneurship Project of Guangzhou Railway Polytechnic(2025CXCY015)。
文摘As the dominant seepage channel in rock masses,it is of great significance to study the influence of fracture roughness distribution on seepage and heat transfer in rock masses.In this paper,the fracture roughness distribution functions of the Bakhtiary dam site and Oskarshamn/Forsmark mountain were fitted using statistical methods.The COMSOL Multiphysics finite element software was utilized to analyze the effects of fracture roughness distribution types and empirical formulas for fracture hydraulic aperture on the seepage field and temperature field of rock masses.The results show that:(1)The fracture roughness at the Bakhtiary dam site and Oskarshamn/Forsmark mountain follows lognormal and normal distributions,respectively;(2)For rock masses with the same expected value and standard deviation of fracture roughness,the outflow from rock masses with lognormal distribution of fracture roughness is significantly larger than that of rock masses with normal distribution of fracture roughness;(3)The fracture hydraulic aperture,outflow,and cold front distance of the Li and Jiang model are significantly larger than those of the Barton model;(4)The outflow,hydraulic pressure distribution,and temperature distribution of the Barton model are more sensitive to the fracture roughness distribution type than those of the Li and Jiang model.
基金supported by the National Natural Science Foundation of China(Grant No.52374079)Chongqing Graduate Research Innovation Project(Grant No.CYB22032)Chongqing Talents and Outstanding Scientists Project(Grant No.cstc2024ycjhbgzxm0032).
文摘Deep underground excavation causes considerable unloading effects,leading to a pronounced bias pressure phenomenon.The deformation and seepage characteristics of rock masses under different gas and confining pressures were investigated via triaxial loading and unloading seepage tests.When the influential coefficient of effective confining pressure(β)is less than 0.065,the seepage force considerably weakens the strength of fractured rock masses.Conversely,whenβis greater than 0.065,the opposite is true.Moreover,the increase in the axial load leads to an increase in the precast fracture volumetric strain,which is the main reason for the increase in fracture permeability.This effect is particularly significant during the unloading stage.Based on the test results,a method for calculating the dynamic seepage evolution of rock masses,considering the effects of rock mass damage and fracture deformation,is introduced,and the effectiveness of the calculation is validated.The entire description of the seepage under loading and unloading was accomplished.The equivalent relationship between the lateral and normal stresses on fracture surfaces ranges from 0.001 to 0.1,showing an exponential variation between the lateral stress influence coefficient on normal deformation(χ)and seepage pressure.Before the failure of the rock mass,the seepage in the fractures was in a linear laminar flow state.However,after the failure,when the gas pressure reached 2 MPa,the flow state in the fractures transitioned to nonlinear laminar flow.The results are important for predicting hazardous gas leaks during deep underground engineering excavation.
基金financially supported by the Heilongjiang Provincial Key R&D Program Project(No.GA21A204)Heilongjiang Provincial Natural Science Foundation of China(No.LH2022D031)the Research Project of Heilongjiang Province Bureau of Geology and Mineral Resources(No.HKY202302).
文摘Although previous researchers have attempted to decipher ore genesis and mineralization in the Erdaokan Ag-Pb-Zn deposit,some uncertainties regarding the mineralization process and evolution of both ore-forming fluids and magnetite types still need to be addressed.In this study,we obtained new EPMA,LA-ICP-MS,and in situ Fe isotope data from magnetite from the Erdaokan deposit,in order to better understand the mineralization mechanism and evolution of both magnetite and the ore-forming fluids.Our results identified seven types of magnetite at Erdaokan:disseminated magnetite(Mag1),coarse-grained magnetite(Mag2a),radial magnetite(Mag2b),fragmented fine-grained magnetite(Mag2c),vermicular gel magnetite(Mag3a1 and Mag3a2),colloidal magnetite(Mag3b)and dark gray magnetite(Mag4).All of the magnetite types were hydrothermal in origin and generally low in Ti(<400 ppm)and Ni(<800 ppm),while being enriched in light Fe isotopes(δ^(56)Fe ranging from−1.54‰to−0.06‰).However,they exhibit different geochemical signatures and are thus classified into high-manganese magnetite(Mag1,MnO>5 wt%),low-silicon magnetite(Mag2a-c,SiO_(2)<1 wt%),high-silicon magnetite(Mag3a-b,SiO_(2)from 1 to 7 wt%)and high-silicon-manganese magnetite(Mag4,SiO_(2)>1 wt%,MnO>0.2 wt%),each being formed within distinct hydrothermal environments.Based on mineralogy,elemental geochemistry,Fe isotopes,temperature trends,TMg-mag and(Ti+V)vs.(Al+Mn)diagrams,we propose that the Erdaokan Ag-Pb-Zn deposit underwent multi-stage mineralization,which can be broken down into four stages and nine sub-stages.Mag1,Mag2a-c,Mag3a-b and Mag4 were formed during the first sub-stage of each of the four stages,respectively.Additionally,fluid mixing,cooling and depressurization boiling were identified as the main mechanisms for mineral precipitation.The enrichment of Ag was significantly enhanced by the superposition of multi-stage ore-forming hydrothermal fluids in the Erdaokan Ag-Pb-Zn deposit.
