The existence of aligned fractures in fluid-saturated rocks leads to obvious attenuation anisotropy and velocity anisotropy. Attenuation anisotropy analysis can be applied to estimate fracture density and scale, which...The existence of aligned fractures in fluid-saturated rocks leads to obvious attenuation anisotropy and velocity anisotropy. Attenuation anisotropy analysis can be applied to estimate fracture density and scale, which provide important information for reservoir identification. This paper derives P-wave attenuation anisotropy in the ATI media where the symmetry axis is in the arbitrary direction theoretically and modifies the spectral ratio method to measure attenuation anisotropy in the ATI media, thus avoiding a large measurement error when applied to wide azimuth or full azimuth data. Fracture dip and azimuth can be estimated through attenuation anisotropy analysis. For small-scale fractures, fracture scale and fracture density can be determined with enhanced convergence if velocity and attenuation information are both used. We also apply the modified spectralratio method to microseismic field data from an oilfield in East China and extract the fracture dip through attenuation anisotropy analysis. The result agrees with the microseismie monitoring.展开更多
Study on failure of soft stratum of roller compacted concrete (RCC) is an important aspect of stability of high RCC dam. Six kinds of specimens with different interfaces were investigated by wedge splitting method. Do...Study on failure of soft stratum of roller compacted concrete (RCC) is an important aspect of stability of high RCC dam. Six kinds of specimens with different interfaces were investigated by wedge splitting method. Double-K fracture parameters (initial fracture parameter and unsteady fracture parameter) were calculated by the concrete double-K fracture theory. It is indicated that the approach of construction joint or old joint after RCC final set is the most efficient among the six cases, and its fracture parameter is the largest among them. Its propagation path is sinuous. Its failure surface is scraggly. Bedding plane crack fails at the underside of the concrete surface, bond course or the surface between them for each approach. So disturbance on the underside of the concrete surface should be avoided or decreased at best during RCC construction.展开更多
Based on the actual data collected from the tight sandstone development zone, correlation analysis using theSpearman method was conducted to determine the main factors influencing the gas production rate of tightsands...Based on the actual data collected from the tight sandstone development zone, correlation analysis using theSpearman method was conducted to determine the main factors influencing the gas production rate of tightsandstone fracturing. An integrated model combining geological engineering and numerical simulation of fracturepropagation and production was completed. Based on data analysis, the hydraulic fracture parameters wereoptimized to develop a differentiated fracturing treatment adjustment plan. The results indicate that the influenceof geological and engineering factors in the X1 and X2 development zones in the study area differs significantly.Therefore, it is challenging to adopt a uniform development strategy to achieve rapid production increase. Thedata analysis reveals that the variation in gas production rate is primarily affected by the reservoir thickness andpermeability parameters as geological factors. On the other hand, the amount of treatment fluid and proppantaddition significantly impact the gas production rate as engineering factors. Among these factors, the influence ofgeological factors is more pronounced in block X1. Therefore, the main focus should be on further optimizing thefracturing interval and adjusting the geological development well location. Given the existing well location, thereis limited potential for further optimizing fracture parameters to increase production. For block X2, the fracturingparameters should be optimized. Data screening was conducted to identify outliers in the entire dataset, and adata-driven fracturing parameter optimization method was employed to determine the basic adjustment directionfor reservoir stimulation in the target block. This approach provides insights into the influence of geological,stimulation, and completion parameters on gas production rate. Consequently, the subsequent fracturing parameteroptimization design can significantly reduce the modeling and simulation workload and guide field operations toimprove and optimize hydraulic fracturing efficiency.展开更多
Multi-fractured horizontal wells are commonly employed to improve the productivity of low and ultra-low permeability gas reservoirs.However,conventional productivity models for open-hole multi-fractured horizontal wel...Multi-fractured horizontal wells are commonly employed to improve the productivity of low and ultra-low permeability gas reservoirs.However,conventional productivity models for open-hole multi-fractured horizontal wells do not consider the interferences between hydraulic fractures and the open-hole segments,resulting in significant errors in calculation results.In this article,a novel productivity prediction model for gas reservoirs with open-hole multi-fractured horizontal wells was proposed based on complex potential theories,potential superimposition,and numerical analysis.Herein,an open-hole segment between two adjacent fractures was regarded as an equivalent fracture,which was discretized as in cases of artificial fractures.The proposed model was then applied to investigate the effects of various parameters,such as the angle between the fracture and horizontal shaft,fracture quantity,fracture length,diversion capacity of fractures,horizontal well length,and inter-fracture distance,on the productivity of low permeability gas reservoirs with multi-fractured horizontal wells.Simulation results revealed that the quantity,length,and distribution of fractures had significant effects on the productivity of low permeability gas reservoirs while the effects of the diversion capacity of fractures and the angle between the fracture and horizontal shaft were negligible.Additionally,a U-shaped distribution of fracture lengths was preferential as the quantity of fractures at shaft ends was twice that in the middle area.展开更多
The accurate estimation of fracture geometry parameters and the characterization of rock mass structure are two important topics in the geological disposal system of high-level radioactive waste(HLW).The Beishan area,...The accurate estimation of fracture geometry parameters and the characterization of rock mass structure are two important topics in the geological disposal system of high-level radioactive waste(HLW).The Beishan area,as the current preselected area for China’s HLW disposal,has three subareas considered to be the key survey area at the stage of site selection.In this paper,a comprehensive survey method conducted on the outcrop is developed to estimate fracture geometry parameters.Results show that fracture occurrence obeys a Fisher distribution,fracture trace length obeys a normal distribution,and the distribution of spacing obeys a negative exponential distribution.An evaluation index,Rock Mass Structure Rating(RMSR),is proposed to characterize rock mass structure for the three subareas.The results show that the Xinchang area is more suitable to act as China’s HLW disposal repository site.At the same time,the index can also be applied to characterize surface rock mass structure and rock mass integrity at the site selection phase of HLW disposal.展开更多
Fracture parameters of welded joints with different strength matching and crack depth in weld metal are investigated by using the methods of elastoplastic finite element analysis and three point bend specimen test. Th...Fracture parameters of welded joints with different strength matching and crack depth in weld metal are investigated by using the methods of elastoplastic finite element analysis and three point bend specimen test. The results show that for shallow crack, the plastic zone turns large in loading process, and the fracture toughness turns high. The extent of the plastic zone of overmatched joint is larger than that of undermatched joint because it will extends to parent metal from the weld metal in loading process for the same CTOD value. The plastic zone of undermatched joint is restricted within the weld, and the size of that is small. Overmatched joint shows the fracture behaviour of shallow crack may more easily than the undermatched joint, while the two sorts of joint specimens have the same crack depth. Therefore, the fracture-resistant capability of overmatched weld is better than that of undermatched weld. when the toughness of weld metals is similar for both overmatched and undermatched joints.展开更多
An open-source MATLAB application(app)named Discontinuity Intensity Calculator and Estimator(DICE)was developed in order to quantitatively characterize the fractures,or in more general,discontinuities within a rocky o...An open-source MATLAB application(app)named Discontinuity Intensity Calculator and Estimator(DICE)was developed in order to quantitatively characterize the fractures,or in more general,discontinuities within a rocky outcrop in three-dimensional(3D)digital data,such as digital outcrop model(DOM).The workflow proposed for the parametrization of the discontinuities consists of the following steps:(1)Analysis and mapping of the fractures detected within the 3D DOMs;(2)Calculation of the orientation,position and dimensions of discontinuities that are represented by best-fit circular planes;(3)Determining the discontinuity parameters(dimension,distribution,spacing and intensity)by the DICE algorithm using different 3D oriented sampling techniques(3D oriented scanline,3D oriented circular scan window and spherical scan volume).Different sampling methods were bench tested with a synthetic,as well as a natural case study,and compared in order to understand the advantages and limitations of each technique.The 3D oriented circular scan window appears to be the most effective method for fracture intensity estimation with high accuracy(error 0.4%)and stability with variations in scan radius.展开更多
This paper establishes a 3D multi-well pad fracturing numerical model coupled with fracture propagation and proppant migration based on the displacement discontinuity method and Eulerian-Eulerian frameworks,and the fr...This paper establishes a 3D multi-well pad fracturing numerical model coupled with fracture propagation and proppant migration based on the displacement discontinuity method and Eulerian-Eulerian frameworks,and the fracture propagation and proppant distribution during multi-well fracturing are investigated by taking the actual multi-well pad parameters as an example.Fracture initiation and propagation during multi-well pad fracturing are jointly affected by a variety of stress interference mechanisms such as inter-cluster,inter-stage,and inter-well,and the fracture extension is unbalanced among clusters,asymmetric on both wings,and dipping at heels.Due to the significant influence of fracture morphology and width on the migration capacity of proppant in the fracture,proppant is mainly placed in the area near the wellbore with large fracture width,while a high-concentration sandwash may easily occur in the area with narrow fracture width as a result of quick bridging.On the whole,the proppant placement range is limited.Increasing the well-spacing can reduce the stress interference of adjacent wells and promote the uniform distribution of fractures and proppant on both wings.The maximum stimulated reservoir volume or multi-fracture uniform propagation can be achieved by optimizing the well spacing.Although reducing the perforation-cluster spacing also can improve the stimulated reservoir area,a too low cluster spacing is not conducive to effectively increasing the propped fracture area.Since increasing the stage time lag is beneficial to reduce inter-stage stress interference,zipper fracturing produces more uniform fracture propagation and proppant distribution.展开更多
Ultra-low permeability reservoirs are characterized by small pore throats and poor physical properties, which areat the root of well-known problems related to injection and production. In this study, a gas injection f...Ultra-low permeability reservoirs are characterized by small pore throats and poor physical properties, which areat the root of well-known problems related to injection and production. In this study, a gas injection floodingapproach is analyzed in the framework of numerical simulations. In particular, the sequence and timing of fracturechanneling and the related impact on production are considered for horizontal wells with different fracturemorphologies. Useful data and information are provided about the regulation of gas channeling and possible strategiesto delay gas channeling and optimize the gas injection volume and fracture parameters. It is shown that inorder to mitigate gas channeling and ensure high production, fracture length on the sides can be controlled andlonger fractures can be created in the middle by which full gas flooding is obtained at the fracture location in themiddle of the horizontal well. A Differential Evolution (DE) algorithm is provided by which the gas injectionvolume and the fracture parameters of gas injection flooding can be optimized. It is shown that an improvedoil recovery factor as high as 6% can be obtained.展开更多
By means of fracture testing on roller-compacted concrete (RCC) three-point bending beams with two different specimen sizes, the P-CMOD complete curve for RCC was gained. Furthermore, by applying double-K fracture t...By means of fracture testing on roller-compacted concrete (RCC) three-point bending beams with two different specimen sizes, the P-CMOD complete curve for RCC was gained. Furthermore, by applying double-K fracture theory, KiniⅠC,KunⅠC, as well as the critical effective crack length and the critical crack tip opening displacement, were evaluated. Based on the double-K fracture parameters above, the calculation model of equivalent strength for induced crack was established, thus the calculation method on its initiation, stable propagation and unstable fracture was ascertained. Moreover, the finite element simulation analysis of stress field in ShaPai arch dam and the on-site observational splaying points of induced crack at different altitudes validated the reliability of the model. Finally, crack inducer′s optimal setting in RCC arch dam was studied. It improves the design level of induced crack in RCC arch dam and satisfies the necessity of engineering practice.展开更多
CO_(2) pre-injection during hydraulic fracturing is an important method for the development of medium to deep heavy oil reservoirs.It reduces the interfacial tension and viscosity of crude oil,enhances its flowability...CO_(2) pre-injection during hydraulic fracturing is an important method for the development of medium to deep heavy oil reservoirs.It reduces the interfacial tension and viscosity of crude oil,enhances its flowability,maintains reservoir pressure,and increases reservoir drainage capacity.Taking the Badaowan Formation as an example,in this study a detailed three-dimensional geomechanical model based on static data from well logging interpretations is elaborated,which can take into account both vertical and horizontal geological variations and mechanical characteristics.A comprehensive analysis of the impact of key construction parameters on Pre-CO_(2) based fracturing(such as cluster spacing and injection volume),is therefore conducted.Thereafter,using optimized construction parameters,a non-structured grid for dynamic development prediction is introduced,and the capacity variations of different production scenarios are assessed.On the basis of the simulation results,reasonable fracturing parameters are finally determined,including cluster spacing,fracturing fluid volume,proppant concentration,and well spacing.展开更多
An accurate and straightforward symplectic method is presented for the fracture analysis of fractional two-dimensional(2D)viscoelastic media.The fractional Kelvin-Zener constitutive model is used to describe the time-...An accurate and straightforward symplectic method is presented for the fracture analysis of fractional two-dimensional(2D)viscoelastic media.The fractional Kelvin-Zener constitutive model is used to describe the time-dependent behavior of viscoelastic materials.Within the framework of symplectic elasticity,the governing equations in the Hamiltonian form for the frequency domain(s-domain)can be directly and rigorously calculated.In the s-domain,the analytical solutions of the displacement and stress fields are constructed by superposing the symplectic eigensolutions without any trial function,and the explicit expressions of the intensity factors and J-integral are derived simultaneously.Comparison studies are provided to validate the accuracy and effectiveness of the present solutions.A detailed analysis is made to reveal the effects of viscoelastic parameters and applied loads on the intensity factors and J-integral.展开更多
To deal with the stress interference caused by simultaneous propagation of multiple fractures and the wettability reversal and physical property changes of the reservoir caused by fracturing fluid getting in during la...To deal with the stress interference caused by simultaneous propagation of multiple fractures and the wettability reversal and physical property changes of the reservoir caused by fracturing fluid getting in during large-volume fracturing of tight oil reservoirs through a horizontal well, a non-planar 3D fracture growth model was built, wettability reversal characterizing parameters and change of relative permeability curve were introduced to correct the production prediction model of fractured horizontal well, a fracturing design optimization software(Fr Smart) by integrating geological and engineering data was developed, and a fracturing design optimization approach for tight oil reservoirs based on fracture control was worked out. The adaptability of the method was analyzed and the fracture parameters of horizontal wells in tight oil reservoirs were optimized. The simulation results show that fracturing technology based on fracture control is suitable for tight oil reservoirs, and by optimizing fracture parameters, this technology makes it possible to produce the maximum amount of reserves in the well-controlled unit of unconventional reservoirs. The key points of fracturing design optimization based on fracture control include increasing lateral length of and reducing the row spacing between horizontal wells, increasing perforation clusters in one stage to decrease the spacing of neighboring fractures, and also avoiding interference of old and new fracturing wells. Field tests show that this technology can increase single well production and ultimate recovery. Using this technology in developing unconventional resources such as tight oil reservoirs in China will enhance the economics significantly.展开更多
In part I and II of this series, experimental investigation in both EPFM and LEFM had been discussed. In this part, further theoretical analysis is given. The theoretical development of Two Parameter Fracture Mechanic...In part I and II of this series, experimental investigation in both EPFM and LEFM had been discussed. In this part, further theoretical analysis is given. The theoretical development of Two Parameter Fracture Mechanics by Hancock etc, has rationalized our experimental results. This method can be applied to engineering practice, and will allow the advantage of enhanced toughness for specimens with low levels of constraint to be taken into account for defect assessment.展开更多
Supercritical carbon dioxide(ScCO_(2))fracturing has great advantages and prospects in both shale gas exploitation and CO_(2)storage.This paper reviews current laboratory experimental methods and results for sedimenta...Supercritical carbon dioxide(ScCO_(2))fracturing has great advantages and prospects in both shale gas exploitation and CO_(2)storage.This paper reviews current laboratory experimental methods and results for sedimentary rocks fractured by ScCO_(2).The breakdown pressure,fracture parameters,mineral composition,bedding plane angle and permeability are discussed.We also compare the differences between sedimentary rock and granite fractured by ScCO_(2),ultimately noting problems and suggesting solutions and strategies for the future.The analysis found that the breakdown pressure of ScCO_(2)was reduced 6.52%–52.31%compared with that of using water.ScCO_(2)tends to produce a complex fracture morphology with significantly higher permeability.When compared with water,the fracture aperture of ScCO_(2)was decreased by 4.10%–72.33%,the tortuosity of ScCO_(2)was increased by 5.41%–70.98%and the fractal dimension of ScCO_(2)was increased by 4.55%–8.41%.The breakdown pressure of sandstone is more sensitive to the nature of the fracturing fluid,but fracture aperture is less sensitive to fracturing fluid than for shale and coal.Compared with granite,the tortuosity of sedimentary rock is more sensitive to the fracturing fluid and the fracture fractal dimension is less sensitive to the fracturing fluid.Existing research shows that ScCO_(2)has the advantages of low breakdown pressure,good fracture creation and environmental protection.It is recommended that research be conducted in terms of sample terms,experimental conditions,effectiveness evaluation and theoretical derivation in order to promote the application of ScCO_(2)reformed reservoirs in the future.展开更多
By systematically summarizing horizontal well fracturing technology abroad for shale oil and gas reservoirs since the “13th Five-Year Plan”, this article elaborates new horizontal well fracturing features in 3D deve...By systematically summarizing horizontal well fracturing technology abroad for shale oil and gas reservoirs since the “13th Five-Year Plan”, this article elaborates new horizontal well fracturing features in 3D development of stacked shale reservoirs, small well spacing and dense well pattern, horizontal well re-fracturing, fracturing parameters optimization and cost control. In light of requirements on horizontal well fracturing technology in China, we have summarized the technological progress in simulation of multi-fracture propagation, horizontal well frac-design, electric-drive fracturing equipment, soluble tools and low-cost downhole materials and factory-like operation. On this basis, combined with the demand analysis of horizontal well fracturing technology in the “14th Five-Year Plan” for unconventional shale oil and gas, we suggest strengthening the research and development in the following 7 aspects:(1) geology-engineering integration;(2) basic theory and design optimization of fracturing for shale oil and gas reservoirs;(3) development of high-power electric-drive fracturing equipment;(4) fracturing tool and supporting equipment for long horizontal section;(5) horizontal well flexible-sidetracking drilling technology for tapping remaining oil;(6) post-frac workover technology for long horizontal well;(7) intelligent fracturing technology.展开更多
A modified polarization saturation model is proposed and addressed math- ematically using a complex variable approach in two-dimensional (2D) semipermeable piezoelectric media. In this model, an existing polarizatio...A modified polarization saturation model is proposed and addressed math- ematically using a complex variable approach in two-dimensional (2D) semipermeable piezoelectric media. In this model, an existing polarization saturation (PS) model in 2D piezoelectric media is modified by considering a linearly varying saturated normal electric displacement load in place of a constant normal electric displacement load, applied on a saturated electric zone. A centre cracked infinite 2D piezoelectric domain subject to an arbitrary poling direction and in-plane electromechanical loadings is considered for the analytical and numerical studies. Here, the problem is mathematically modeled as a non-homogeneous Riemann-Hilbert problem in terms of unknown complex potential functions representing electric displacement and stress components. Having solved the Hilbert problem, the solutions to the saturated zone length, the crack opening displace- ment (COD), the crack opening potential (COP), and the local stress intensity factors (SIFs) are obtained in explicit forms. A numerical study is also presented for the proposed modified model, showing the effects of the saturation condition on the applied electrical loading, the saturation zone length, and the COP. The results of fracture parameters obtained from the proposed model are compared with the existing PS model subject to electrical loading, crack face conditions, and polarization angles.展开更多
It is difficult to determine the main controlling factors of tight oil production.In addition to the problem of uncontrollable prediction accuracy,the numerical prediction model established by the main controlling fac...It is difficult to determine the main controlling factors of tight oil production.In addition to the problem of uncontrollable prediction accuracy,the numerical prediction model established by the main controlling factors will also make the correctly predicted low production samples lose the value of development.Applying the optimization algorithm with fast convergence speed and global optimization to optimize the controllable parameters in the high-precision numerical prediction model can effectively improve the productivity of low production wells with timeliness,and bring greater economic value while saving development cost.Using PCA-GRA method,the sample weight and the weighted correlation ranking results of parameters affecting tight oil production were obtained.Thereupon then the main controlling factors of tight oil production were determined.Then we set up a BP neural network model with by taking the main controlling factors as input and tight oil production as output.The prediction effect of the network was good and can be put into use.The results of sensitivity analysis showed that the network was stable,and the total fracturing fluid volume had the greatest impact on the production of tight oil.Finally,by using genetic algorithm,we optimized the fracturing parameters of all low production well samples in the field data.Combined with the fracturing parameters of all high production well samples and the optimized fracturing parameters of low production wells,the optimal interval of fracturing parameters was given,which can provide guidance for the field fracturing operation.展开更多
基金supported by 973 Program of China(No.2013CB228602)National Science and Technology Major Project of China(No.2016ZX05004003-002)863 Program of China(No.2013AA064202)
文摘The existence of aligned fractures in fluid-saturated rocks leads to obvious attenuation anisotropy and velocity anisotropy. Attenuation anisotropy analysis can be applied to estimate fracture density and scale, which provide important information for reservoir identification. This paper derives P-wave attenuation anisotropy in the ATI media where the symmetry axis is in the arbitrary direction theoretically and modifies the spectral ratio method to measure attenuation anisotropy in the ATI media, thus avoiding a large measurement error when applied to wide azimuth or full azimuth data. Fracture dip and azimuth can be estimated through attenuation anisotropy analysis. For small-scale fractures, fracture scale and fracture density can be determined with enhanced convergence if velocity and attenuation information are both used. We also apply the modified spectralratio method to microseismic field data from an oilfield in East China and extract the fracture dip through attenuation anisotropy analysis. The result agrees with the microseismie monitoring.
基金the National Natural Science Foundation of China (Grant No.50179002)Liaoning Province Dr. Fund Opening Foundation of China.
文摘Study on failure of soft stratum of roller compacted concrete (RCC) is an important aspect of stability of high RCC dam. Six kinds of specimens with different interfaces were investigated by wedge splitting method. Double-K fracture parameters (initial fracture parameter and unsteady fracture parameter) were calculated by the concrete double-K fracture theory. It is indicated that the approach of construction joint or old joint after RCC final set is the most efficient among the six cases, and its fracture parameter is the largest among them. Its propagation path is sinuous. Its failure surface is scraggly. Bedding plane crack fails at the underside of the concrete surface, bond course or the surface between them for each approach. So disturbance on the underside of the concrete surface should be avoided or decreased at best during RCC construction.
基金Research and Application of Key Technologies for Tight Gas Production Improvement and Rehabilitation of Linxing Shenfu(YXKY-ZL-01-2021)。
文摘Based on the actual data collected from the tight sandstone development zone, correlation analysis using theSpearman method was conducted to determine the main factors influencing the gas production rate of tightsandstone fracturing. An integrated model combining geological engineering and numerical simulation of fracturepropagation and production was completed. Based on data analysis, the hydraulic fracture parameters wereoptimized to develop a differentiated fracturing treatment adjustment plan. The results indicate that the influenceof geological and engineering factors in the X1 and X2 development zones in the study area differs significantly.Therefore, it is challenging to adopt a uniform development strategy to achieve rapid production increase. Thedata analysis reveals that the variation in gas production rate is primarily affected by the reservoir thickness andpermeability parameters as geological factors. On the other hand, the amount of treatment fluid and proppantaddition significantly impact the gas production rate as engineering factors. Among these factors, the influence ofgeological factors is more pronounced in block X1. Therefore, the main focus should be on further optimizing thefracturing interval and adjusting the geological development well location. Given the existing well location, thereis limited potential for further optimizing fracture parameters to increase production. For block X2, the fracturingparameters should be optimized. Data screening was conducted to identify outliers in the entire dataset, and adata-driven fracturing parameter optimization method was employed to determine the basic adjustment directionfor reservoir stimulation in the target block. This approach provides insights into the influence of geological,stimulation, and completion parameters on gas production rate. Consequently, the subsequent fracturing parameteroptimization design can significantly reduce the modeling and simulation workload and guide field operations toimprove and optimize hydraulic fracturing efficiency.
基金This work was supported by grants from the National Natural Science Foundation of China(51574197)Educational Commission of Sichuan Province of China(16ZA0071).
文摘Multi-fractured horizontal wells are commonly employed to improve the productivity of low and ultra-low permeability gas reservoirs.However,conventional productivity models for open-hole multi-fractured horizontal wells do not consider the interferences between hydraulic fractures and the open-hole segments,resulting in significant errors in calculation results.In this article,a novel productivity prediction model for gas reservoirs with open-hole multi-fractured horizontal wells was proposed based on complex potential theories,potential superimposition,and numerical analysis.Herein,an open-hole segment between two adjacent fractures was regarded as an equivalent fracture,which was discretized as in cases of artificial fractures.The proposed model was then applied to investigate the effects of various parameters,such as the angle between the fracture and horizontal shaft,fracture quantity,fracture length,diversion capacity of fractures,horizontal well length,and inter-fracture distance,on the productivity of low permeability gas reservoirs with multi-fractured horizontal wells.Simulation results revealed that the quantity,length,and distribution of fractures had significant effects on the productivity of low permeability gas reservoirs while the effects of the diversion capacity of fractures and the angle between the fracture and horizontal shaft were negligible.Additionally,a U-shaped distribution of fracture lengths was preferential as the quantity of fractures at shaft ends was twice that in the middle area.
基金supported by the National Key Research and Development Program of China,under grant No.2018YFC1504903the Chongqing Natural Science Foundation,under grant No.cstc2020jcyj-msxm X0743 and cstc 2020jcyj-bsh0142+3 种基金the Open Research Fund of the State Key Laboratory of Geomechanics and Geotechnical Engineering,Institute of Rock and Soil Mechanics,Chinese Academy of Sciences,under grant No.Z019018China postdoctoral science foundation Grant No.2019M662918 and 2020M673152Regional Joint Fund for Basic and Applied Basic Research Fund of Guangdong Province,No.2019A1515110836the National Natural Science Foundation of China,under grant No.41688103。
文摘The accurate estimation of fracture geometry parameters and the characterization of rock mass structure are two important topics in the geological disposal system of high-level radioactive waste(HLW).The Beishan area,as the current preselected area for China’s HLW disposal,has three subareas considered to be the key survey area at the stage of site selection.In this paper,a comprehensive survey method conducted on the outcrop is developed to estimate fracture geometry parameters.Results show that fracture occurrence obeys a Fisher distribution,fracture trace length obeys a normal distribution,and the distribution of spacing obeys a negative exponential distribution.An evaluation index,Rock Mass Structure Rating(RMSR),is proposed to characterize rock mass structure for the three subareas.The results show that the Xinchang area is more suitable to act as China’s HLW disposal repository site.At the same time,the index can also be applied to characterize surface rock mass structure and rock mass integrity at the site selection phase of HLW disposal.
文摘Fracture parameters of welded joints with different strength matching and crack depth in weld metal are investigated by using the methods of elastoplastic finite element analysis and three point bend specimen test. The results show that for shallow crack, the plastic zone turns large in loading process, and the fracture toughness turns high. The extent of the plastic zone of overmatched joint is larger than that of undermatched joint because it will extends to parent metal from the weld metal in loading process for the same CTOD value. The plastic zone of undermatched joint is restricted within the weld, and the size of that is small. Overmatched joint shows the fracture behaviour of shallow crack may more easily than the undermatched joint, while the two sorts of joint specimens have the same crack depth. Therefore, the fracture-resistant capability of overmatched weld is better than that of undermatched weld. when the toughness of weld metals is similar for both overmatched and undermatched joints.
文摘An open-source MATLAB application(app)named Discontinuity Intensity Calculator and Estimator(DICE)was developed in order to quantitatively characterize the fractures,or in more general,discontinuities within a rocky outcrop in three-dimensional(3D)digital data,such as digital outcrop model(DOM).The workflow proposed for the parametrization of the discontinuities consists of the following steps:(1)Analysis and mapping of the fractures detected within the 3D DOMs;(2)Calculation of the orientation,position and dimensions of discontinuities that are represented by best-fit circular planes;(3)Determining the discontinuity parameters(dimension,distribution,spacing and intensity)by the DICE algorithm using different 3D oriented sampling techniques(3D oriented scanline,3D oriented circular scan window and spherical scan volume).Different sampling methods were bench tested with a synthetic,as well as a natural case study,and compared in order to understand the advantages and limitations of each technique.The 3D oriented circular scan window appears to be the most effective method for fracture intensity estimation with high accuracy(error 0.4%)and stability with variations in scan radius.
基金Supported by National Natural Science Foundation of China(51974332)Strategic Cooperation Project Between PetroChina and China University of Petroleum(Beijing)(ZLZX2020-07).
文摘This paper establishes a 3D multi-well pad fracturing numerical model coupled with fracture propagation and proppant migration based on the displacement discontinuity method and Eulerian-Eulerian frameworks,and the fracture propagation and proppant distribution during multi-well fracturing are investigated by taking the actual multi-well pad parameters as an example.Fracture initiation and propagation during multi-well pad fracturing are jointly affected by a variety of stress interference mechanisms such as inter-cluster,inter-stage,and inter-well,and the fracture extension is unbalanced among clusters,asymmetric on both wings,and dipping at heels.Due to the significant influence of fracture morphology and width on the migration capacity of proppant in the fracture,proppant is mainly placed in the area near the wellbore with large fracture width,while a high-concentration sandwash may easily occur in the area with narrow fracture width as a result of quick bridging.On the whole,the proppant placement range is limited.Increasing the well-spacing can reduce the stress interference of adjacent wells and promote the uniform distribution of fractures and proppant on both wings.The maximum stimulated reservoir volume or multi-fracture uniform propagation can be achieved by optimizing the well spacing.Although reducing the perforation-cluster spacing also can improve the stimulated reservoir area,a too low cluster spacing is not conducive to effectively increasing the propped fracture area.Since increasing the stage time lag is beneficial to reduce inter-stage stress interference,zipper fracturing produces more uniform fracture propagation and proppant distribution.
基金supported by the Forward Looking Basic Major Scientific and Technological Projects of CNPC (Grant No.2021DJ2202).
文摘Ultra-low permeability reservoirs are characterized by small pore throats and poor physical properties, which areat the root of well-known problems related to injection and production. In this study, a gas injection floodingapproach is analyzed in the framework of numerical simulations. In particular, the sequence and timing of fracturechanneling and the related impact on production are considered for horizontal wells with different fracturemorphologies. Useful data and information are provided about the regulation of gas channeling and possible strategiesto delay gas channeling and optimize the gas injection volume and fracture parameters. It is shown that inorder to mitigate gas channeling and ensure high production, fracture length on the sides can be controlled andlonger fractures can be created in the middle by which full gas flooding is obtained at the fracture location in themiddle of the horizontal well. A Differential Evolution (DE) algorithm is provided by which the gas injectionvolume and the fracture parameters of gas injection flooding can be optimized. It is shown that an improvedoil recovery factor as high as 6% can be obtained.
文摘By means of fracture testing on roller-compacted concrete (RCC) three-point bending beams with two different specimen sizes, the P-CMOD complete curve for RCC was gained. Furthermore, by applying double-K fracture theory, KiniⅠC,KunⅠC, as well as the critical effective crack length and the critical crack tip opening displacement, were evaluated. Based on the double-K fracture parameters above, the calculation model of equivalent strength for induced crack was established, thus the calculation method on its initiation, stable propagation and unstable fracture was ascertained. Moreover, the finite element simulation analysis of stress field in ShaPai arch dam and the on-site observational splaying points of induced crack at different altitudes validated the reliability of the model. Finally, crack inducer′s optimal setting in RCC arch dam was studied. It improves the design level of induced crack in RCC arch dam and satisfies the necessity of engineering practice.
基金supported by the Cutting-Edge Project Foundation of Petro-China(Cold-Based Method to Enhance Heavy Oil Recovery)(Grant No.2021DJ1406)Open Fund(PLN201802)of National Key Laboratory of Oil and Gas Reservoir Geology and Exploitation(Southwest Petroleum University).
文摘CO_(2) pre-injection during hydraulic fracturing is an important method for the development of medium to deep heavy oil reservoirs.It reduces the interfacial tension and viscosity of crude oil,enhances its flowability,maintains reservoir pressure,and increases reservoir drainage capacity.Taking the Badaowan Formation as an example,in this study a detailed three-dimensional geomechanical model based on static data from well logging interpretations is elaborated,which can take into account both vertical and horizontal geological variations and mechanical characteristics.A comprehensive analysis of the impact of key construction parameters on Pre-CO_(2) based fracturing(such as cluster spacing and injection volume),is therefore conducted.Thereafter,using optimized construction parameters,a non-structured grid for dynamic development prediction is introduced,and the capacity variations of different production scenarios are assessed.On the basis of the simulation results,reasonable fracturing parameters are finally determined,including cluster spacing,fracturing fluid volume,proppant concentration,and well spacing.
基金Project supported by the National Natural Science Foundation of China(Nos.11872303 and 11702221)the China Postdoctoral Science Foundation(No.2017M613198)the Fundamental Research Funds for the Central Universities of China(No.G2020KY05402)
文摘An accurate and straightforward symplectic method is presented for the fracture analysis of fractional two-dimensional(2D)viscoelastic media.The fractional Kelvin-Zener constitutive model is used to describe the time-dependent behavior of viscoelastic materials.Within the framework of symplectic elasticity,the governing equations in the Hamiltonian form for the frequency domain(s-domain)can be directly and rigorously calculated.In the s-domain,the analytical solutions of the displacement and stress fields are constructed by superposing the symplectic eigensolutions without any trial function,and the explicit expressions of the intensity factors and J-integral are derived simultaneously.Comparison studies are provided to validate the accuracy and effectiveness of the present solutions.A detailed analysis is made to reveal the effects of viscoelastic parameters and applied loads on the intensity factors and J-integral.
基金Supported by China National Science and Technology Major Project(2016ZX05023,2017ZX05013-005)
文摘To deal with the stress interference caused by simultaneous propagation of multiple fractures and the wettability reversal and physical property changes of the reservoir caused by fracturing fluid getting in during large-volume fracturing of tight oil reservoirs through a horizontal well, a non-planar 3D fracture growth model was built, wettability reversal characterizing parameters and change of relative permeability curve were introduced to correct the production prediction model of fractured horizontal well, a fracturing design optimization software(Fr Smart) by integrating geological and engineering data was developed, and a fracturing design optimization approach for tight oil reservoirs based on fracture control was worked out. The adaptability of the method was analyzed and the fracture parameters of horizontal wells in tight oil reservoirs were optimized. The simulation results show that fracturing technology based on fracture control is suitable for tight oil reservoirs, and by optimizing fracture parameters, this technology makes it possible to produce the maximum amount of reserves in the well-controlled unit of unconventional reservoirs. The key points of fracturing design optimization based on fracture control include increasing lateral length of and reducing the row spacing between horizontal wells, increasing perforation clusters in one stage to decrease the spacing of neighboring fractures, and also avoiding interference of old and new fracturing wells. Field tests show that this technology can increase single well production and ultimate recovery. Using this technology in developing unconventional resources such as tight oil reservoirs in China will enhance the economics significantly.
文摘In part I and II of this series, experimental investigation in both EPFM and LEFM had been discussed. In this part, further theoretical analysis is given. The theoretical development of Two Parameter Fracture Mechanics by Hancock etc, has rationalized our experimental results. This method can be applied to engineering practice, and will allow the advantage of enhanced toughness for specimens with low levels of constraint to be taken into account for defect assessment.
基金funded by the National Natural Science Foundation of China(Grant Nos.42141009,41825018,41888101 and 41902289)the Key Research Program of the Institute of Geology and Geophysics,CAS(Grant No.IGGCAS-202201)the Second Tibetan Plateau Scientific Expedition and Research Program(STEP)(Grant No.2019QZKK0904)。
文摘Supercritical carbon dioxide(ScCO_(2))fracturing has great advantages and prospects in both shale gas exploitation and CO_(2)storage.This paper reviews current laboratory experimental methods and results for sedimentary rocks fractured by ScCO_(2).The breakdown pressure,fracture parameters,mineral composition,bedding plane angle and permeability are discussed.We also compare the differences between sedimentary rock and granite fractured by ScCO_(2),ultimately noting problems and suggesting solutions and strategies for the future.The analysis found that the breakdown pressure of ScCO_(2)was reduced 6.52%–52.31%compared with that of using water.ScCO_(2)tends to produce a complex fracture morphology with significantly higher permeability.When compared with water,the fracture aperture of ScCO_(2)was decreased by 4.10%–72.33%,the tortuosity of ScCO_(2)was increased by 5.41%–70.98%and the fractal dimension of ScCO_(2)was increased by 4.55%–8.41%.The breakdown pressure of sandstone is more sensitive to the nature of the fracturing fluid,but fracture aperture is less sensitive to fracturing fluid than for shale and coal.Compared with granite,the tortuosity of sedimentary rock is more sensitive to the fracturing fluid and the fracture fractal dimension is less sensitive to the fracturing fluid.Existing research shows that ScCO_(2)has the advantages of low breakdown pressure,good fracture creation and environmental protection.It is recommended that research be conducted in terms of sample terms,experimental conditions,effectiveness evaluation and theoretical derivation in order to promote the application of ScCO_(2)reformed reservoirs in the future.
基金Supported by the National Science and Technology Major Project(2016ZX05023)。
文摘By systematically summarizing horizontal well fracturing technology abroad for shale oil and gas reservoirs since the “13th Five-Year Plan”, this article elaborates new horizontal well fracturing features in 3D development of stacked shale reservoirs, small well spacing and dense well pattern, horizontal well re-fracturing, fracturing parameters optimization and cost control. In light of requirements on horizontal well fracturing technology in China, we have summarized the technological progress in simulation of multi-fracture propagation, horizontal well frac-design, electric-drive fracturing equipment, soluble tools and low-cost downhole materials and factory-like operation. On this basis, combined with the demand analysis of horizontal well fracturing technology in the “14th Five-Year Plan” for unconventional shale oil and gas, we suggest strengthening the research and development in the following 7 aspects:(1) geology-engineering integration;(2) basic theory and design optimization of fracturing for shale oil and gas reservoirs;(3) development of high-power electric-drive fracturing equipment;(4) fracturing tool and supporting equipment for long horizontal section;(5) horizontal well flexible-sidetracking drilling technology for tapping remaining oil;(6) post-frac workover technology for long horizontal well;(7) intelligent fracturing technology.
文摘A modified polarization saturation model is proposed and addressed math- ematically using a complex variable approach in two-dimensional (2D) semipermeable piezoelectric media. In this model, an existing polarization saturation (PS) model in 2D piezoelectric media is modified by considering a linearly varying saturated normal electric displacement load in place of a constant normal electric displacement load, applied on a saturated electric zone. A centre cracked infinite 2D piezoelectric domain subject to an arbitrary poling direction and in-plane electromechanical loadings is considered for the analytical and numerical studies. Here, the problem is mathematically modeled as a non-homogeneous Riemann-Hilbert problem in terms of unknown complex potential functions representing electric displacement and stress components. Having solved the Hilbert problem, the solutions to the saturated zone length, the crack opening displace- ment (COD), the crack opening potential (COP), and the local stress intensity factors (SIFs) are obtained in explicit forms. A numerical study is also presented for the proposed modified model, showing the effects of the saturation condition on the applied electrical loading, the saturation zone length, and the COP. The results of fracture parameters obtained from the proposed model are compared with the existing PS model subject to electrical loading, crack face conditions, and polarization angles.
基金The authors gratefully acknowledge the financial support of the National Science and Technology Major Projects of China(2016ZX05065 and 2016ZX05042-003).
文摘It is difficult to determine the main controlling factors of tight oil production.In addition to the problem of uncontrollable prediction accuracy,the numerical prediction model established by the main controlling factors will also make the correctly predicted low production samples lose the value of development.Applying the optimization algorithm with fast convergence speed and global optimization to optimize the controllable parameters in the high-precision numerical prediction model can effectively improve the productivity of low production wells with timeliness,and bring greater economic value while saving development cost.Using PCA-GRA method,the sample weight and the weighted correlation ranking results of parameters affecting tight oil production were obtained.Thereupon then the main controlling factors of tight oil production were determined.Then we set up a BP neural network model with by taking the main controlling factors as input and tight oil production as output.The prediction effect of the network was good and can be put into use.The results of sensitivity analysis showed that the network was stable,and the total fracturing fluid volume had the greatest impact on the production of tight oil.Finally,by using genetic algorithm,we optimized the fracturing parameters of all low production well samples in the field data.Combined with the fracturing parameters of all high production well samples and the optimized fracturing parameters of low production wells,the optimal interval of fracturing parameters was given,which can provide guidance for the field fracturing operation.