Production decline analysis has been considered as an important method to obtain the flow parameters, reservoir properties and original gas in place. Although advanced Blasingame production decline analysis methods fo...Production decline analysis has been considered as an important method to obtain the flow parameters, reservoir properties and original gas in place. Although advanced Blasingame production decline analysis methods for vertical wells, fractured wells and horizontal wells are widely used, limited study has conducted on Blasingame production decline type curves for multi-fractured horizontal well(MFHW). Based on the perpendicular bisection(PEBI) grids, a numerical model was developed and the solution was obtained using control volume finite element method and the fully implicit method. Blasingame production decline-type curves of the infinitely conductive MFHW were plotted through computer programming. A field case was presented to analyse and verify the model developed. Five flow regimes, including early formation linear flow, early radial flow, compound linear flow, transient flow and pseudo-radial flow, are recognized. Fracture spacing is the main factor that affects early radial flow, compound linear flow and transient flow, the distance from the well to the circular boundary affects the pseudo-radial flow, and the type curves are also significantly affected by the formation permeability, fracture number and fracture half-length. The validation of field case suggests that the Blasingame production decline type curves proposed in this work can be applied to the production decline analysis for MFHW in tight gas reservoirs.展开更多
The Michelson Interferometer for Global High-resolution Thermospheric Imaging(MIGHTI)onboard the Ionospheric Connection Explorer(ICON)satellite offers the opportunity to investigate the altitude profile of thermospher...The Michelson Interferometer for Global High-resolution Thermospheric Imaging(MIGHTI)onboard the Ionospheric Connection Explorer(ICON)satellite offers the opportunity to investigate the altitude profile of thermospheric winds.In this study,we used the red-line measurements of MIGHTI to compare with the results estimated by Horizontal Wind Model 14(HWM14).The data selected included both the geomagnetic quiet period(December 2019 to August 2022)and the geomagnetic storm on August 26-28,2021.During the geomagnetic quiet period,the estimations of neutral winds from HWM14 showed relatively good agreement with the observations from ICON.According to the ICON observations,near the equator,zonal winds reverse from westward to eastward at around 06:00 local time(LT)at higher altitudes,and the stronger westward winds appear at later LTs at lower altitudes.At around 16:00 LT,eastward winds at 300 km reverse to westward,and vertical gradients of zonal winds similar to those at sunrise hours can be observed.In the middle latitudes,zonal winds reverse about 2-4 h earlier.Meridional winds vary more significantly than zonal winds with seasonal and latitudinal variations.According to the ICON observations,in the northern low latitudes,vertical reversals of meridional winds are found at 08:00-13:00 LT from 300 to 160 km and at around 18:00 LT from 300 to 200 km during the June solstice.Similar reversals of meridional winds are found at 04:00-07:00 LT from 300 to 160 km and at 22:00-02:00 LT from 270 to 200 km during the December solstice.In the southern low latitudes,meridional wind reversals occur at 08:00-11:00 LT from 200 to 160 km and at 21:00-02:00 LT from 300 to 200 km during the June solstice.During the December solstice,reversals of the meridional wind appear at 20:00-01:00 LT below 200 km and at 06:00-11:00 LT from 300 to 160 km.In the northern middle latitudes,the northward winds are dominant at 08:00-14:00 LT at 230 km during the June solstice.Northward winds persist until 16:00 LT at 160 and 300 km.During the December solstice,the northward winds are dominant from 06:00 to 21:00 LT.The vertical variations in neutral winds during the geomagnetic storm on August 26-28 were analyzed in detail.Both meridional and zonal winds during the active geomagnetic period observed by ICON show distinguishable vertical shear structures at different stages of the storm.On the dayside,during the main phase,the peak velocities of westward winds extend from a higher altitude to a lower altitude,whereas during the recovery phase,the peak velocities of the westward winds extend from lower altitudes to higher altitudes.The velocities of the southward winds are stronger at lower altitudes during the storm.These vertical structures of horizontal winds during the storm could not be reproduced by the HWM14 wind estimations,and the overall response to the storm of the horizontal winds in the low and middle latitudes is underestimated by HWM14.The ICON observations provide a good dataset for improving the HWM wind estimations in the middle and upper atmosphere,especially the vertical variations.展开更多
Understanding the mechanical properties of the lithologies is crucial to accurately determine the horizontal stress magnitude.To investigate the correlation between the rock mass properties and maximum horizontal stre...Understanding the mechanical properties of the lithologies is crucial to accurately determine the horizontal stress magnitude.To investigate the correlation between the rock mass properties and maximum horizontal stress,the three-dimensional(3D)stress tensors at 89 measuring points determined using an improved overcoring technique in nine mines in China were adopted,a newly defined characteristic parameter C_(ERP)was proposed as an indicator for evaluating the structural properties of rock masses,and a fuzzy relation matrix was established using the information distribution method.The results indicate that both the vertical stress and horizontal stress exhibit a good linear growth relationship with depth.There is no remarkable correlation between the elastic modulus,Poisson's ratio and depth,and the distribution of data points is scattered and messy.Moreover,there is no obvious relationship between the rock quality designation(RQD)and depth.The maximum horizontal stress σ_(H) is a function of rock properties,showing a certain linear relationship with the C_(ERP)at the same depth.In addition,the overall change trend of σ_(H) determined by the established fuzzy identification method is to increase with the increase of C_(ERP).The fuzzy identification method also demonstrates a relatively detailed local relationship betweenσ_H and C_(ERP),and the predicted curve rises in a fluctuating way,which is in accord well with the measured stress data.展开更多
Porous carbonate reservoirs,prevalent in the Middle East,are lithologically dominated by bioclastic limestones,exhibiting high porosity,low permeability,intricate pore structure,and strong heterogeneity.Waterflooding ...Porous carbonate reservoirs,prevalent in the Middle East,are lithologically dominated by bioclastic limestones,exhibiting high porosity,low permeability,intricate pore structure,and strong heterogeneity.Waterflooding through horizontal wells is commonly used for exploiting these reservoirs.However,challenges persist,such as significant uncertainty and complex operational procedures regarding adjustment effects during the exploitation.The USH reservoir of the Cretaceous D oilfield,Oman exemplifies typical porous carbonate reservoirs.It initially underwent depletion drive using vertical wells,followed by horizontal well waterflooding in the late stage.Currently,the oilfield is confronted with substantial developmental challenges,involving the understanding of residual oil distribution,effective water cut control,and sustaining oil production since it has entered the late development stage.Employing a microscopic visualization displacement system equipped with electrodes,this study elucidated the waterflooding mechanisms and residual oil distribution in the late-stage development of the USH reservoir.The results reveal that the reservoir's vertical stacking patterns act as the main factor affecting the horizontal well waterflooding efficacy.Distinct water flow channels emerge under varying reservoir stacking patterns,with post-waterflooding residual oil predominantly distributed at the reservoir's top and bottom.The oil recovery can be enhanced by adjusting the waterflooding's flow line and intensity.The findings of this study will provide theoretical insights of waterflooding mechanisms and injection-production adjustments for exploiting other porous carbonate reservoirs in the Middle East through horizontal wells.展开更多
Rock slope with horizontal-layered fractured structure(HLFS)has high stability in its natural state.However,a strong earthquake can induce rock fissure expansion,ultimately leading to slope failure.In this study,the d...Rock slope with horizontal-layered fractured structure(HLFS)has high stability in its natural state.However,a strong earthquake can induce rock fissure expansion,ultimately leading to slope failure.In this study,the dynamic response,failure mode,and spectral characteristics of rock slope with HLFS under strong earthquake conditions were investigated based on the large-scale shaking table model test.On this basis,multiple sets of numerical calculation models were further established by UDEC discrete element program.Five influencing factors were considered in the parametric study of numerical simulations,including slope height,slope angle,bedding-plane spacing and secondary joint spacing as well as bedrock dip angle.The results showed that the failure process of rock slope with HLFS under earthquake action is mainly divided into four phases,i.e.,the tensile crack of the slope shoulder joints and shear dislocation at the top bedding plane,the extension of vertical joint cracks and increase of shear displacement,the formation of step-through sliding surfaces and the instability,and finally collapse of fractured rock mass.The acceleration response of slopes exhibits elevation amplification effect and surface effect.Numerical simulations indicate that the seismic stability of slopes with HLFS exhibits a negative correlation with slope height and angle,but a positive correlation with bedding-plane spacing,joint spacing,and bedrock dip angle.The results of this study can provide a reference for seismic stability evaluation of weathered rock slopes.展开更多
Due to the high potential risk and many influencing factors of subsea horizontal X-tree installation,to guarantee the successful completion of sea trials of domestic subsea horizontal X-trees,this paper established a ...Due to the high potential risk and many influencing factors of subsea horizontal X-tree installation,to guarantee the successful completion of sea trials of domestic subsea horizontal X-trees,this paper established a modular risk evaluation model based on a fuzzy fault tree.First,through the analysis of the main process oftree down and combining the Offshore&Onshore Reliability Data(OREDA)failure statistics and the operation procedure and the data provided by the job,the fault tree model of risk analysis of the tree down installation was established.Then,by introducing the natural language of expert comprehensive evaluation and combining fuzzy principles,quantitative analysis was carried out,and the fuzzy number was used to calculate the failure probability of a basic event and the occurrence probability of a top event.Finally,through a sensitivity analysis of basic events,the basic events of top events significantly affected were determined,and risk control and prevention measures for the corresponding high-risk factors were proposed for subsea horizontal X-tree down installation.展开更多
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 forward model for optical fiber strain was established based on a planar 3D multi-fracture model. Then the forward method calculating distributed fiber strain induced by multi-fracture growth was proposed. Based on ...A forward model for optical fiber strain was established based on a planar 3D multi-fracture model. Then the forward method calculating distributed fiber strain induced by multi-fracture growth was proposed. Based on this method, fiber strain evolution during fracturing of the horizontal well was numerically simulated. Fiber strain evolution induced by fracture growth can be divided into three stages: strain increasing, shrinkage convergence, and straight-line convergence, whereas the evolution of fiber strain rate has four stages: strain rate increasing, shrinkage convergence, straight-line convergence, and strain rate reversal after pumping stops. Fiber strain does not flip after pumping stop, while the strain rate flips after pumping stop so that strain rate can reflect injection dynamics. The time when the fracture extends to the fiber and inter-well pressure channeling can be identified by the straight-line convergence band of distributed fiber strain or strain rate, and the non-uniform growth of multiple fractures can be evaluated by using the instants of fractures reaching the fiber monitoring well.When the horizontal section of the fiber monitoring well is within the height range of a hydraulic fracture, the instant of the fracture reaching the fiber can be identified;otherwise, the converging band is not apparent. In multi-stage fracturing, under the influence of stress shadow from previous fracturing stages, the tensile region of fiber strain may not appear, but the fiber strain rate can effectively show the fracture growth behavior in each stage. The evolution law of fiber strain rate in single-stage fracturing can be applied to multi-stage fracturing.展开更多
Aiming at the problems of large load of rotation drive system,low efficiency of torque transmission and high cost for operation and maintenance of liner steering drilling system for the horizontal well,a new method of...Aiming at the problems of large load of rotation drive system,low efficiency of torque transmission and high cost for operation and maintenance of liner steering drilling system for the horizontal well,a new method of liner differential rotary drilling with double tubular strings in the horizontal well is proposed.The technical principle of this method is revealed,supporting tools such as the differential rotation transducer,composite rotary steering system and the hanger are designed,and technological process is optimized.A tool face control technique of steering drilling assembly is proposed and the calculation model of extension limit of liner differential rotary drilling with double tubular strings in horizontal well is established.These results show that the liner differential rotary drilling with double tubular strings is equipped with measurement while drilling(MWD)and positive displacement motor(PDM),and directional drilling of horizontal well is realized by adjusting rotary speed of drill pipe to control the tool face of PDM.Based on the engineering case of deep coalbed methane horizontal well in the eastern margin of Ordos Basin,the extension limit of horizontal drilling with double tubular strings is calculated.Compared with the conventional liner drilling method,the liner differential rotary drilling with double tubular strings increases the extension limit value of horizontal well significantly.The research findings provide useful reference for the integrated design and control of liner completion and drilling of horizontal wells.展开更多
The boundary condition is a crucial factor affecting the permeability variation due to suffusion.An experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering the boundar...The boundary condition is a crucial factor affecting the permeability variation due to suffusion.An experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering the boundary effect is conducted,where the hydraulic head difference(DH)varies,and the boundary includes non-loss and soil-loss conditions.Soil samples are filled into seven soil storerooms connected in turn.After evaluation,the variation in content of fine sand(ΔR_(f))and the hydraulic conductivity of soils in each storeroom(C_(i))are analyzed.In the non-loss test,the soil sample filling area is divided into runoff,transited,and accumulated areas according to the negative or positive ΔR_(f) values.ΔR_(f) increases from negative to positive along the seepage path,and Ci decreases from runoff area to transited area and then rebounds in accumulated area.In the soil-loss test,all soil sample filling areas belong to the runoff area,where the gentle-loss,strengthened-loss,and alleviated-loss parts are further divided.ΔR_(f) decreases from the gentle-loss part to the strengthened-loss part and then rebounds in the alleviated-loss part,and C_(i) increases and then decreases along the seepage path.The relationship between ΔR_(f) and Ci is different with the boundary condition.Ci exponentially decreases with ΔR_(f) in the non-loss test and increases with ΔR_(f) generally in the soil-loss test.展开更多
In this paper,we consider pseudoharmonic heat flow with small initial horizontal energy and give the existence of pseudoharmonic maps from closed pseudo-Hermitian manifolds into closed Riemannian manifolds.
Multistage fracturing of horizontal wells is a critical technology for unconventional oil and gas reservoir stimulation. Ball-throwing temporary plugging fracturing is a new method for realizing uniform fracturing alo...Multistage fracturing of horizontal wells is a critical technology for unconventional oil and gas reservoir stimulation. Ball-throwing temporary plugging fracturing is a new method for realizing uniform fracturing along horizontal wells and plays an important role in increasing oil and gas production. However,the transportation and sealing law of temporary plugging balls(TPBs) in the perforation section of horizontal wells is still unclear. Using COMSOL computational fluid dynamics and a particle tracking module, we simulate the transportation process of TPBs in a horizontal wellbore and analyse the effects of the ball density, ball diameter, ball number, fracturing fluid injection rate, and viscosity on the plugging efficiency of TPB transportation. This study reveals that when the density of TPBs is close to that of the fracturing fluid and a moderate diameter of the TPB is used, the plugging efficiency can be substantially enhanced. The plugging efficiency is greater when the TPB number is close to twice the number of perforations and is lower when the number of TPBs is three times the number of perforations.Adjusting the fracturing fluid injection rate from low to high can control the position of the TPBs,improving plugging efficiency. As the viscosity of the fracturing fluid increases, the plugging efficiency of the perforations decreases near the borehole heel and increases near the borehole toe. In contrast, the plugging efficiency of the central perforation is almost unaffected by the fracturing fluid viscosity. This study can serve as a valuable reference for establishing the parameters for temporary plugging and fracturing.展开更多
The horizontal bearing capacity of the screw pile and monopile was analyzed by model tests.Results showed that the horizontal bearing capacity of the screw pile was significantly greater than that of the monopile unde...The horizontal bearing capacity of the screw pile and monopile was analyzed by model tests.Results showed that the horizontal bearing capacity of the screw pile was significantly greater than that of the monopile under the same loading conditions.With the increase in horizontal loading speed,the ultimate horizontal bearing capacity of the two piles also increases,and the difference decreases gradually.Moreover,the influence of vertical loading on the horizontal bearing capacity of screw pile and monopile is studied at the horizontal loading speed of 2 mm s-1.The findings indicate that vertical load evidently affects the horizontal bearing capacity of common piles,but slightly influences the horizontal bearing capacity of screw piles.展开更多
Many research findings have proven that the system of porous medium reservoirs exhibits different heterogeneous structures at various scales,demonstrating some form of self-similarity with fractal characteristics.In t...Many research findings have proven that the system of porous medium reservoirs exhibits different heterogeneous structures at various scales,demonstrating some form of self-similarity with fractal characteristics.In this paper,fractal theory is incorporated into the reservoir to investigate coupled flow between reservoir and horizontal well.By examining the pore structure of highly heterogeneous reservoirs,the fractal dimension can be determined.Analytical methods are utilized to solve the Green function of a point source in a reservoir with fractal characteristics.Employing Green's function and the principle of spatial superposition,a finite flow model for a horizontal well coupled with a fractal reservoir is developed to calculate the flow rate and flow profile of the horizontal well.The model also accounts for the impact of wellbore friction and is solved numerically.A specific example is used for calculation to analyze the influence of fractal parameters on the production and flow rate of the horizontal well.When considering the fractal characteristics of oil reservoirs,the flow rate of the horizontal well is lower than that in Euclidean space.As the fractal dimension increases,the connectivity of pores in the reservoir improves,making it easier to drive the fluid into the wellbore,and the flow distribution along the wellbore becomes more uniform.Conversely,as the anomalous diffusion index increases,the connectivity between pores deteriorates,thus the distribution of flow rate along the wellbore becomes more uneven.展开更多
This study proposes a novel open-type rectangular breakwater combined with horizontal perforated plates on both sides to enhance the sheltering effect of the rectangular box-type breakwaters against longer waves.The h...This study proposes a novel open-type rectangular breakwater combined with horizontal perforated plates on both sides to enhance the sheltering effect of the rectangular box-type breakwaters against longer waves.The hydrodynamic characteristics of this breakwater are analyzed through analytical potential solutions and experimental tests.The quadratic pressure drop conditions are exerted on the horizontal perforated plates to facilitate assessing the effect of wave height on the dissipated wave energy of breakwater through the analytical solution.The hydrodynamic quantities of the breakwater,including the reflection,transmission,and energyloss coefficients,together with vertical and horizontal wave forces,are calculated using the velocity potential decomposition method as well as an iterative algorithm.Furthermore,the reflection and transmission coefficients of the breakwater are measured by conducting experimental tests at various wave periods,wave heights,and both porosities and widths of the horizontal perforated plates.The analytical predicted results demonstrate good agreement with the iterative boundary element method solution and measured data.The influences of variable incident waves and structure parameters on the hydrodynamic characteristics of the breakwater are investigated through further calculations based on analytical solutions.Results indicate that horizontal perforated plates placed on the water surface for both sides of the rectangular breakwater can enhance the wave dissipation ability of the breakwater while effectively decreasing the transmission and reflection coefficients.展开更多
The vacuum-assisted prefabricated horizontal drain offers a promising method for strengthening soil slurry,allowing simultaneous filling and vacuum-dewatering via staged construction.However,there is limited research ...The vacuum-assisted prefabricated horizontal drain offers a promising method for strengthening soil slurry,allowing simultaneous filling and vacuum-dewatering via staged construction.However,there is limited research on the unique characteristics of staged filling.This study aims to investigate the vacuum consolidation process of staged-filled soil slurry through laboratory model tests and numerical simulations,also assessing the impact of anionic polyacrylamide.Comparative analyses are conducted between vacuum consolidation with and without anionic polyacrylamide,as well as self-weight consolidation without anionic polyacrylamide.Results reveal contour lines of excess pore pressure,water content,and soil strength forming an ellipse around the prefabricated horizontal drain board.During the consolidation process,a higher degree of consolidation,lower water content,and higher soil strength were observed closer to the prefabricated horizontal drain board.After treatment,the uppermost filling layer exhibits an average water content that was approximately 40%higher than the lower filling layer,and its average strength was about 60%lower.This discrepancy is primarily due to the absence of sealing on the top surface and the relatively short vacuum consolidation time caused by staged filling.The introduction of anionic polyacrylamide-induced flocculation significantly improves the initial consolidation rate but minimally affects the dewatering capacity of vacuum preloading.Using flocculant can enhance both the staged filling rate and soil strength(by 1e2 times).Additionally,employing a staggered arrangement between different prefabricated horizontal drain layers is advisable to prevent top-down penetration in areas with low soil strength.展开更多
A coupled numerical calculation method combining smooth particle hydrodynamics(SPH)and the finite element method(FEM)was implemented to investigate the seismic response of horizontal storage tanks.Anumericalmodel of a...A coupled numerical calculation method combining smooth particle hydrodynamics(SPH)and the finite element method(FEM)was implemented to investigate the seismic response of horizontal storage tanks.Anumericalmodel of a horizontal storage tank featuring a free liquid surface under seismic action was constructed using the SPH–FEM coupling method.The stored liquid was discretized using SPH particles,while the tank and supports were discretized using the FEM.The interaction between the stored liquid and the tank was simulated by using the meshless particle contact method.Then,the numerical simulation results were compared and analyzed against seismic simulation shaking table test data to validate the method.Subsequently,a series of numerical models,considering different liquid storage volumes and seismic effects,were constructed to obtain time history data of base shear and top center displacement,which revealed the seismic performance of horizontal storage tanks.Numerical simulation results and experimental data showed good agreement,with an error rate of less than 18.85%.And this conformity signifies the rationality of the SPH-FEM coupling method.The base shear and top center displacement values obtained by the coupled SPH-FEM method were only 53.3% to 69.1% of those calculated by the equivalent mass method employed in the current code.As the stored liquid volume increased,the seismic response of the horizontal storage tank exhibited a gradual upward trend,with the seismic response increasing from 73% to 388% for every 35% increase in stored liquid volume.The maximum von Mises stress of the tank and the supports remained below the steel yield strength during the earthquake.The coupled SPH-FEM method holds certain advantages in studying the seismic problems of tanks with complex structural forms,particularly due to the representation of the flow field distribution during earthquakes by involving reservoir fluid participation.展开更多
A genuine technological issue–the thermal convection of liquid in a rotating cavity–is investigated experimentally.The experiments are conducted within a horizontal annulus with isothermal boundaries. The inner boun...A genuine technological issue–the thermal convection of liquid in a rotating cavity–is investigated experimentally.The experiments are conducted within a horizontal annulus with isothermal boundaries. The inner boundaryof the annulus has a higher temperature, thus exerting a stabilising influence on the system. It is shown that whenthe layer rotation velocity diminishes, two-dimensional azimuthally periodic convective rolls, rotating togetherwith the cavity, emerge in a threshold manner. The development of convection is accompanied by a significantintensification of heat transfer through the layer. It is shown that the averaged thermal convection excitation inthe form of a system of two-dimensional rolls occurs against the background of oscillations of a non-isothermalfluid in the cavity reference frame caused by the gravity field. The excitation threshold and the structure ofconvective rolls are consistent with the results of the earlier theoretical studies by the authors performed usingthe equations of “vibrational” convection obtained by the averaging method. Furthermore, the experiments haverevealed a new type of averaged flow in the form of a spatially periodic system of toroidal vortices. It is shown thata steady streaming, excited by the inertial oscillations of the fluid, is responsible for the generation of the toroidalvortices. These flows develop in a non-threshold manner and are most clearly manifested in a case of resonantexcitation of one of the inertial modes.展开更多
A comprehensive method to evaluate the factors affecting the production capacity of horizontal wells in Carboniferous volcanic rocks after fracturing is investigated.A systematic approach combining gray correlation an...A comprehensive method to evaluate the factors affecting the production capacity of horizontal wells in Carboniferous volcanic rocks after fracturing is investigated.A systematic approach combining gray correlation analysis,hierarchical analysis and fuzzy evaluation is proposed.In particular,first the incidence of reservoir properties and fracturing parameters on production capacity is assessed.These parameters include reservoir base geological parameters(porosity,permeability,oil saturation,waterproof height)as well as engineering parameters(fracture halflength,fracture height,fracture conductivity,fracture distance).Afterwards,a two-by-two comparison judgment matrix of sensitive parameters is constructed by means of hierarchical analysis,and the weighting coefficients of the factors are determined,where oil saturation,fracture conductivity and fracture half-length are weighted higher.Finally,the horizontal wells in the target block are categorized in terms of production capacity based on the fuzzy evaluation method,and split accordingly into high-producing,relatively high-producing,medium-producing and low-producing wells.Such a categorization is intended to provide parametric guidance for reservoir fracturing and modification.展开更多
The segmented water control technology for bottom water reservoirs can effectively delay the entry of bottom water and adjust the production profile.To clarify the impact of different methods on horizontal well produc...The segmented water control technology for bottom water reservoirs can effectively delay the entry of bottom water and adjust the production profile.To clarify the impact of different methods on horizontal well production with different reservoir conditions and to provide theoretical support for the scientific selection of methods for bottom water reservoirs,a numerical simulation method is presented in this study,which is able to deal with wellbore reservoir coupling under screen tube,perforation,and ICD(Inflow Control Device)completion.Assuming the geological characteristics of the bottom-water conglomerate reservoir in the Triassic Formation of the Tahe Block 9 as a test case,the three aforementioned completion methods are tested to predict the transient production characteristics.The impact of completion parameters,reservoir permeability,bottom-water energy,and individual well control on the time to encounter water in horizontal wells(during a water-free production period)is discussed.A boundary chart for the selection of completion methods is introduced accordingly.The results show that the optimized ICD completion development effect for heterogeneous reservoirs is the best,followed by optimized perforation completion.Permeability is the main factor affecting the performances of completion methods,while bottom water energy and single well controlled reserves have a scarce impact.The average permeability of the reservoir is less than 500 mD,and ICD has the best water control effect.If the permeability is greater than 500 mD,the water control effect of perforation completion becomes a better option.展开更多
基金Project(2013CB228005)supported by the National Basic Research Program of China
文摘Production decline analysis has been considered as an important method to obtain the flow parameters, reservoir properties and original gas in place. Although advanced Blasingame production decline analysis methods for vertical wells, fractured wells and horizontal wells are widely used, limited study has conducted on Blasingame production decline type curves for multi-fractured horizontal well(MFHW). Based on the perpendicular bisection(PEBI) grids, a numerical model was developed and the solution was obtained using control volume finite element method and the fully implicit method. Blasingame production decline-type curves of the infinitely conductive MFHW were plotted through computer programming. A field case was presented to analyse and verify the model developed. Five flow regimes, including early formation linear flow, early radial flow, compound linear flow, transient flow and pseudo-radial flow, are recognized. Fracture spacing is the main factor that affects early radial flow, compound linear flow and transient flow, the distance from the well to the circular boundary affects the pseudo-radial flow, and the type curves are also significantly affected by the formation permeability, fracture number and fracture half-length. The validation of field case suggests that the Blasingame production decline type curves proposed in this work can be applied to the production decline analysis for MFHW in tight gas reservoirs.
基金supported by the National Key R&D Program of China (Grant No.2022YFF0503700)the special funds of Hubei Luojia Laboratory (Grant No.220100011)+1 种基金supported by the International Space Science Institute–Beijing(ISSI-BJ) project“The Electromagnetic Data Validation and Scientific Application Research based on CSES Satellite”and ISSI/ISSI-BJ project,“Multi-Scale Magnetosphere–Ionosphere–Thermosphere Interaction.”
文摘The Michelson Interferometer for Global High-resolution Thermospheric Imaging(MIGHTI)onboard the Ionospheric Connection Explorer(ICON)satellite offers the opportunity to investigate the altitude profile of thermospheric winds.In this study,we used the red-line measurements of MIGHTI to compare with the results estimated by Horizontal Wind Model 14(HWM14).The data selected included both the geomagnetic quiet period(December 2019 to August 2022)and the geomagnetic storm on August 26-28,2021.During the geomagnetic quiet period,the estimations of neutral winds from HWM14 showed relatively good agreement with the observations from ICON.According to the ICON observations,near the equator,zonal winds reverse from westward to eastward at around 06:00 local time(LT)at higher altitudes,and the stronger westward winds appear at later LTs at lower altitudes.At around 16:00 LT,eastward winds at 300 km reverse to westward,and vertical gradients of zonal winds similar to those at sunrise hours can be observed.In the middle latitudes,zonal winds reverse about 2-4 h earlier.Meridional winds vary more significantly than zonal winds with seasonal and latitudinal variations.According to the ICON observations,in the northern low latitudes,vertical reversals of meridional winds are found at 08:00-13:00 LT from 300 to 160 km and at around 18:00 LT from 300 to 200 km during the June solstice.Similar reversals of meridional winds are found at 04:00-07:00 LT from 300 to 160 km and at 22:00-02:00 LT from 270 to 200 km during the December solstice.In the southern low latitudes,meridional wind reversals occur at 08:00-11:00 LT from 200 to 160 km and at 21:00-02:00 LT from 300 to 200 km during the June solstice.During the December solstice,reversals of the meridional wind appear at 20:00-01:00 LT below 200 km and at 06:00-11:00 LT from 300 to 160 km.In the northern middle latitudes,the northward winds are dominant at 08:00-14:00 LT at 230 km during the June solstice.Northward winds persist until 16:00 LT at 160 and 300 km.During the December solstice,the northward winds are dominant from 06:00 to 21:00 LT.The vertical variations in neutral winds during the geomagnetic storm on August 26-28 were analyzed in detail.Both meridional and zonal winds during the active geomagnetic period observed by ICON show distinguishable vertical shear structures at different stages of the storm.On the dayside,during the main phase,the peak velocities of westward winds extend from a higher altitude to a lower altitude,whereas during the recovery phase,the peak velocities of the westward winds extend from lower altitudes to higher altitudes.The velocities of the southward winds are stronger at lower altitudes during the storm.These vertical structures of horizontal winds during the storm could not be reproduced by the HWM14 wind estimations,and the overall response to the storm of the horizontal winds in the low and middle latitudes is underestimated by HWM14.The ICON observations provide a good dataset for improving the HWM wind estimations in the middle and upper atmosphere,especially the vertical variations.
基金financially supported by the National Natural Science Foundation of China(No.52204084)the Open Research Fund of the State Key Laboratory of Coal Resources and safe Mining,CUMT,China(No.SKLCRSM 23KF004)+3 种基金the Interdisciplinary Research Project for Young Teachers of USTB(Fundamental Research Funds for the Central Universities),China(No.FRF-IDRY-GD22-002)the Fundamental Research Funds for the Central Universities and the Youth Teacher International Exchange and Growth Program,China(No.QNXM20220009)the National Key R&D Program of China(Nos.2022YFC2905600 and 2022 YFC3004601)the Science,Technology&Innovation Project of Xiongan New Area,China(No.2023XAGG0061)。
文摘Understanding the mechanical properties of the lithologies is crucial to accurately determine the horizontal stress magnitude.To investigate the correlation between the rock mass properties and maximum horizontal stress,the three-dimensional(3D)stress tensors at 89 measuring points determined using an improved overcoring technique in nine mines in China were adopted,a newly defined characteristic parameter C_(ERP)was proposed as an indicator for evaluating the structural properties of rock masses,and a fuzzy relation matrix was established using the information distribution method.The results indicate that both the vertical stress and horizontal stress exhibit a good linear growth relationship with depth.There is no remarkable correlation between the elastic modulus,Poisson's ratio and depth,and the distribution of data points is scattered and messy.Moreover,there is no obvious relationship between the rock quality designation(RQD)and depth.The maximum horizontal stress σ_(H) is a function of rock properties,showing a certain linear relationship with the C_(ERP)at the same depth.In addition,the overall change trend of σ_(H) determined by the established fuzzy identification method is to increase with the increase of C_(ERP).The fuzzy identification method also demonstrates a relatively detailed local relationship betweenσ_H and C_(ERP),and the predicted curve rises in a fluctuating way,which is in accord well with the measured stress data.
基金funded by a Major Science and Technology Project of China National Petroleum Corporation(CNPC)entitled Research on Key Technologies for Efficient Production of Overseas Large Carbonate Reservoir”(2023ZZ19-08).
文摘Porous carbonate reservoirs,prevalent in the Middle East,are lithologically dominated by bioclastic limestones,exhibiting high porosity,low permeability,intricate pore structure,and strong heterogeneity.Waterflooding through horizontal wells is commonly used for exploiting these reservoirs.However,challenges persist,such as significant uncertainty and complex operational procedures regarding adjustment effects during the exploitation.The USH reservoir of the Cretaceous D oilfield,Oman exemplifies typical porous carbonate reservoirs.It initially underwent depletion drive using vertical wells,followed by horizontal well waterflooding in the late stage.Currently,the oilfield is confronted with substantial developmental challenges,involving the understanding of residual oil distribution,effective water cut control,and sustaining oil production since it has entered the late development stage.Employing a microscopic visualization displacement system equipped with electrodes,this study elucidated the waterflooding mechanisms and residual oil distribution in the late-stage development of the USH reservoir.The results reveal that the reservoir's vertical stacking patterns act as the main factor affecting the horizontal well waterflooding efficacy.Distinct water flow channels emerge under varying reservoir stacking patterns,with post-waterflooding residual oil predominantly distributed at the reservoir's top and bottom.The oil recovery can be enhanced by adjusting the waterflooding's flow line and intensity.The findings of this study will provide theoretical insights of waterflooding mechanisms and injection-production adjustments for exploiting other porous carbonate reservoirs in the Middle East through horizontal wells.
基金supported by Central Guiding Local Science and Technology Development Special Fund Project(No.ZYYD2023B02)the National Natural Science Foundation of China(Nos.52078432 and 52168066)the Scientific Research Project of China Railway First Survey and Design Institute Group Co.(No.20-06).
文摘Rock slope with horizontal-layered fractured structure(HLFS)has high stability in its natural state.However,a strong earthquake can induce rock fissure expansion,ultimately leading to slope failure.In this study,the dynamic response,failure mode,and spectral characteristics of rock slope with HLFS under strong earthquake conditions were investigated based on the large-scale shaking table model test.On this basis,multiple sets of numerical calculation models were further established by UDEC discrete element program.Five influencing factors were considered in the parametric study of numerical simulations,including slope height,slope angle,bedding-plane spacing and secondary joint spacing as well as bedrock dip angle.The results showed that the failure process of rock slope with HLFS under earthquake action is mainly divided into four phases,i.e.,the tensile crack of the slope shoulder joints and shear dislocation at the top bedding plane,the extension of vertical joint cracks and increase of shear displacement,the formation of step-through sliding surfaces and the instability,and finally collapse of fractured rock mass.The acceleration response of slopes exhibits elevation amplification effect and surface effect.Numerical simulations indicate that the seismic stability of slopes with HLFS exhibits a negative correlation with slope height and angle,but a positive correlation with bedding-plane spacing,joint spacing,and bedrock dip angle.The results of this study can provide a reference for seismic stability evaluation of weathered rock slopes.
基金financially supported by the National Ministry of Industry and Information Technology Innovation Special Project-Engineering Demonstration Application of Subsea Production System,Topic 4:Research on Subsea X-Tree and Wellhead Offshore Testing Technology(Grant No.MC-201901-S01-04)the Key Research and Development Program of Shandong Province(Major Innovation Project)(Grant Nos.2022CXGC020405,2023CXGC010415)。
文摘Due to the high potential risk and many influencing factors of subsea horizontal X-tree installation,to guarantee the successful completion of sea trials of domestic subsea horizontal X-trees,this paper established a modular risk evaluation model based on a fuzzy fault tree.First,through the analysis of the main process oftree down and combining the Offshore&Onshore Reliability Data(OREDA)failure statistics and the operation procedure and the data provided by the job,the fault tree model of risk analysis of the tree down installation was established.Then,by introducing the natural language of expert comprehensive evaluation and combining fuzzy principles,quantitative analysis was carried out,and the fuzzy number was used to calculate the failure probability of a basic event and the occurrence probability of a top event.Finally,through a sensitivity analysis of basic events,the basic events of top events significantly affected were determined,and risk control and prevention measures for the corresponding high-risk factors were proposed for subsea horizontal X-tree down installation.
基金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.
基金Supported by the National Natural Science Foundation of China(52104060)Project of Natural Science Outstanding Youth Foundation of Shandong(ZR2019QEE005)。
文摘A forward model for optical fiber strain was established based on a planar 3D multi-fracture model. Then the forward method calculating distributed fiber strain induced by multi-fracture growth was proposed. Based on this method, fiber strain evolution during fracturing of the horizontal well was numerically simulated. Fiber strain evolution induced by fracture growth can be divided into three stages: strain increasing, shrinkage convergence, and straight-line convergence, whereas the evolution of fiber strain rate has four stages: strain rate increasing, shrinkage convergence, straight-line convergence, and strain rate reversal after pumping stops. Fiber strain does not flip after pumping stop, while the strain rate flips after pumping stop so that strain rate can reflect injection dynamics. The time when the fracture extends to the fiber and inter-well pressure channeling can be identified by the straight-line convergence band of distributed fiber strain or strain rate, and the non-uniform growth of multiple fractures can be evaluated by using the instants of fractures reaching the fiber monitoring well.When the horizontal section of the fiber monitoring well is within the height range of a hydraulic fracture, the instant of the fracture reaching the fiber can be identified;otherwise, the converging band is not apparent. In multi-stage fracturing, under the influence of stress shadow from previous fracturing stages, the tensile region of fiber strain may not appear, but the fiber strain rate can effectively show the fracture growth behavior in each stage. The evolution law of fiber strain rate in single-stage fracturing can be applied to multi-stage fracturing.
基金Supported by the Project of National Natural Science Foundation of China(52234002,42230814)。
文摘Aiming at the problems of large load of rotation drive system,low efficiency of torque transmission and high cost for operation and maintenance of liner steering drilling system for the horizontal well,a new method of liner differential rotary drilling with double tubular strings in the horizontal well is proposed.The technical principle of this method is revealed,supporting tools such as the differential rotation transducer,composite rotary steering system and the hanger are designed,and technological process is optimized.A tool face control technique of steering drilling assembly is proposed and the calculation model of extension limit of liner differential rotary drilling with double tubular strings in horizontal well is established.These results show that the liner differential rotary drilling with double tubular strings is equipped with measurement while drilling(MWD)and positive displacement motor(PDM),and directional drilling of horizontal well is realized by adjusting rotary speed of drill pipe to control the tool face of PDM.Based on the engineering case of deep coalbed methane horizontal well in the eastern margin of Ordos Basin,the extension limit of horizontal drilling with double tubular strings is calculated.Compared with the conventional liner drilling method,the liner differential rotary drilling with double tubular strings increases the extension limit value of horizontal well significantly.The research findings provide useful reference for the integrated design and control of liner completion and drilling of horizontal wells.
基金The research work described herein was funded by the National Nature Science Foundation of China(Grant No.41877213).This financial support is gratefully acknowledged.
文摘The boundary condition is a crucial factor affecting the permeability variation due to suffusion.An experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering the boundary effect is conducted,where the hydraulic head difference(DH)varies,and the boundary includes non-loss and soil-loss conditions.Soil samples are filled into seven soil storerooms connected in turn.After evaluation,the variation in content of fine sand(ΔR_(f))and the hydraulic conductivity of soils in each storeroom(C_(i))are analyzed.In the non-loss test,the soil sample filling area is divided into runoff,transited,and accumulated areas according to the negative or positive ΔR_(f) values.ΔR_(f) increases from negative to positive along the seepage path,and Ci decreases from runoff area to transited area and then rebounds in accumulated area.In the soil-loss test,all soil sample filling areas belong to the runoff area,where the gentle-loss,strengthened-loss,and alleviated-loss parts are further divided.ΔR_(f) decreases from the gentle-loss part to the strengthened-loss part and then rebounds in the alleviated-loss part,and C_(i) increases and then decreases along the seepage path.The relationship between ΔR_(f) and Ci is different with the boundary condition.Ci exponentially decreases with ΔR_(f) in the non-loss test and increases with ΔR_(f) generally in the soil-loss test.
文摘In this paper,we consider pseudoharmonic heat flow with small initial horizontal energy and give the existence of pseudoharmonic maps from closed pseudo-Hermitian manifolds into closed Riemannian manifolds.
基金supported by the National Natural Science Foundation of China (No. 52074250)。
文摘Multistage fracturing of horizontal wells is a critical technology for unconventional oil and gas reservoir stimulation. Ball-throwing temporary plugging fracturing is a new method for realizing uniform fracturing along horizontal wells and plays an important role in increasing oil and gas production. However,the transportation and sealing law of temporary plugging balls(TPBs) in the perforation section of horizontal wells is still unclear. Using COMSOL computational fluid dynamics and a particle tracking module, we simulate the transportation process of TPBs in a horizontal wellbore and analyse the effects of the ball density, ball diameter, ball number, fracturing fluid injection rate, and viscosity on the plugging efficiency of TPB transportation. This study reveals that when the density of TPBs is close to that of the fracturing fluid and a moderate diameter of the TPB is used, the plugging efficiency can be substantially enhanced. The plugging efficiency is greater when the TPB number is close to twice the number of perforations and is lower when the number of TPBs is three times the number of perforations.Adjusting the fracturing fluid injection rate from low to high can control the position of the TPBs,improving plugging efficiency. As the viscosity of the fracturing fluid increases, the plugging efficiency of the perforations decreases near the borehole heel and increases near the borehole toe. In contrast, the plugging efficiency of the central perforation is almost unaffected by the fracturing fluid viscosity. This study can serve as a valuable reference for establishing the parameters for temporary plugging and fracturing.
基金funded by the National Natural Science Foundation of China(No.51779171)。
文摘The horizontal bearing capacity of the screw pile and monopile was analyzed by model tests.Results showed that the horizontal bearing capacity of the screw pile was significantly greater than that of the monopile under the same loading conditions.With the increase in horizontal loading speed,the ultimate horizontal bearing capacity of the two piles also increases,and the difference decreases gradually.Moreover,the influence of vertical loading on the horizontal bearing capacity of screw pile and monopile is studied at the horizontal loading speed of 2 mm s-1.The findings indicate that vertical load evidently affects the horizontal bearing capacity of common piles,but slightly influences the horizontal bearing capacity of screw piles.
基金project is funded by the SINOPEC Science and Technology Research Project(No.P24200).
文摘Many research findings have proven that the system of porous medium reservoirs exhibits different heterogeneous structures at various scales,demonstrating some form of self-similarity with fractal characteristics.In this paper,fractal theory is incorporated into the reservoir to investigate coupled flow between reservoir and horizontal well.By examining the pore structure of highly heterogeneous reservoirs,the fractal dimension can be determined.Analytical methods are utilized to solve the Green function of a point source in a reservoir with fractal characteristics.Employing Green's function and the principle of spatial superposition,a finite flow model for a horizontal well coupled with a fractal reservoir is developed to calculate the flow rate and flow profile of the horizontal well.The model also accounts for the impact of wellbore friction and is solved numerically.A specific example is used for calculation to analyze the influence of fractal parameters on the production and flow rate of the horizontal well.When considering the fractal characteristics of oil reservoirs,the flow rate of the horizontal well is lower than that in Euclidean space.As the fractal dimension increases,the connectivity of pores in the reservoir improves,making it easier to drive the fluid into the wellbore,and the flow distribution along the wellbore becomes more uniform.Conversely,as the anomalous diffusion index increases,the connectivity between pores deteriorates,thus the distribution of flow rate along the wellbore becomes more uneven.
基金supported by the National Natural Sci-ence Foundation of China(Nos.52201345,and 52001293)the New Cornerstone Science Foundation through the XPLORER PRIZE.
文摘This study proposes a novel open-type rectangular breakwater combined with horizontal perforated plates on both sides to enhance the sheltering effect of the rectangular box-type breakwaters against longer waves.The hydrodynamic characteristics of this breakwater are analyzed through analytical potential solutions and experimental tests.The quadratic pressure drop conditions are exerted on the horizontal perforated plates to facilitate assessing the effect of wave height on the dissipated wave energy of breakwater through the analytical solution.The hydrodynamic quantities of the breakwater,including the reflection,transmission,and energyloss coefficients,together with vertical and horizontal wave forces,are calculated using the velocity potential decomposition method as well as an iterative algorithm.Furthermore,the reflection and transmission coefficients of the breakwater are measured by conducting experimental tests at various wave periods,wave heights,and both porosities and widths of the horizontal perforated plates.The analytical predicted results demonstrate good agreement with the iterative boundary element method solution and measured data.The influences of variable incident waves and structure parameters on the hydrodynamic characteristics of the breakwater are investigated through further calculations based on analytical solutions.Results indicate that horizontal perforated plates placed on the water surface for both sides of the rectangular breakwater can enhance the wave dissipation ability of the breakwater while effectively decreasing the transmission and reflection coefficients.
基金supported by the Research Grants Council of Hong Kong Special Administrative Region Government of China(Grant Nos.15210322 and R5037-18)the financial support(Grant No.86902-00000240)from Shenzhen University.
文摘The vacuum-assisted prefabricated horizontal drain offers a promising method for strengthening soil slurry,allowing simultaneous filling and vacuum-dewatering via staged construction.However,there is limited research on the unique characteristics of staged filling.This study aims to investigate the vacuum consolidation process of staged-filled soil slurry through laboratory model tests and numerical simulations,also assessing the impact of anionic polyacrylamide.Comparative analyses are conducted between vacuum consolidation with and without anionic polyacrylamide,as well as self-weight consolidation without anionic polyacrylamide.Results reveal contour lines of excess pore pressure,water content,and soil strength forming an ellipse around the prefabricated horizontal drain board.During the consolidation process,a higher degree of consolidation,lower water content,and higher soil strength were observed closer to the prefabricated horizontal drain board.After treatment,the uppermost filling layer exhibits an average water content that was approximately 40%higher than the lower filling layer,and its average strength was about 60%lower.This discrepancy is primarily due to the absence of sealing on the top surface and the relatively short vacuum consolidation time caused by staged filling.The introduction of anionic polyacrylamide-induced flocculation significantly improves the initial consolidation rate but minimally affects the dewatering capacity of vacuum preloading.Using flocculant can enhance both the staged filling rate and soil strength(by 1e2 times).Additionally,employing a staggered arrangement between different prefabricated horizontal drain layers is advisable to prevent top-down penetration in areas with low soil strength.
基金supported by Scientific Research Fund of Institute of Engineering Mechanics,China Earthquake Administration(Grant Nos.2021B06,2021C05)Heilongjiang Natural Science Foundation Joint Guidance Project(Grant No.LH2021E122).
文摘A coupled numerical calculation method combining smooth particle hydrodynamics(SPH)and the finite element method(FEM)was implemented to investigate the seismic response of horizontal storage tanks.Anumericalmodel of a horizontal storage tank featuring a free liquid surface under seismic action was constructed using the SPH–FEM coupling method.The stored liquid was discretized using SPH particles,while the tank and supports were discretized using the FEM.The interaction between the stored liquid and the tank was simulated by using the meshless particle contact method.Then,the numerical simulation results were compared and analyzed against seismic simulation shaking table test data to validate the method.Subsequently,a series of numerical models,considering different liquid storage volumes and seismic effects,were constructed to obtain time history data of base shear and top center displacement,which revealed the seismic performance of horizontal storage tanks.Numerical simulation results and experimental data showed good agreement,with an error rate of less than 18.85%.And this conformity signifies the rationality of the SPH-FEM coupling method.The base shear and top center displacement values obtained by the coupled SPH-FEM method were only 53.3% to 69.1% of those calculated by the equivalent mass method employed in the current code.As the stored liquid volume increased,the seismic response of the horizontal storage tank exhibited a gradual upward trend,with the seismic response increasing from 73% to 388% for every 35% increase in stored liquid volume.The maximum von Mises stress of the tank and the supports remained below the steel yield strength during the earthquake.The coupled SPH-FEM method holds certain advantages in studying the seismic problems of tanks with complex structural forms,particularly due to the representation of the flow field distribution during earthquakes by involving reservoir fluid participation.
基金supported by the Ministry of Education of the Russian Federation(Project KPZU-2023-0002).
文摘A genuine technological issue–the thermal convection of liquid in a rotating cavity–is investigated experimentally.The experiments are conducted within a horizontal annulus with isothermal boundaries. The inner boundaryof the annulus has a higher temperature, thus exerting a stabilising influence on the system. It is shown that whenthe layer rotation velocity diminishes, two-dimensional azimuthally periodic convective rolls, rotating togetherwith the cavity, emerge in a threshold manner. The development of convection is accompanied by a significantintensification of heat transfer through the layer. It is shown that the averaged thermal convection excitation inthe form of a system of two-dimensional rolls occurs against the background of oscillations of a non-isothermalfluid in the cavity reference frame caused by the gravity field. The excitation threshold and the structure ofconvective rolls are consistent with the results of the earlier theoretical studies by the authors performed usingthe equations of “vibrational” convection obtained by the averaging method. Furthermore, the experiments haverevealed a new type of averaged flow in the form of a spatially periodic system of toroidal vortices. It is shown thata steady streaming, excited by the inertial oscillations of the fluid, is responsible for the generation of the toroidalvortices. These flows develop in a non-threshold manner and are most clearly manifested in a case of resonantexcitation of one of the inertial modes.
基金Natural Science Foundation of China(NSFC),Grant No.52404016(Lijuan Huang)Natural Science Foundation of Hubei Province,Grant No.2024AFB322(Lijuan Huang),Open Fund of Hubei Key Laboratory of Oil and Gas Drilling and Production Engineering(Yangtze University),Grant No.YQZC202405(Lijuan Huang).
文摘A comprehensive method to evaluate the factors affecting the production capacity of horizontal wells in Carboniferous volcanic rocks after fracturing is investigated.A systematic approach combining gray correlation analysis,hierarchical analysis and fuzzy evaluation is proposed.In particular,first the incidence of reservoir properties and fracturing parameters on production capacity is assessed.These parameters include reservoir base geological parameters(porosity,permeability,oil saturation,waterproof height)as well as engineering parameters(fracture halflength,fracture height,fracture conductivity,fracture distance).Afterwards,a two-by-two comparison judgment matrix of sensitive parameters is constructed by means of hierarchical analysis,and the weighting coefficients of the factors are determined,where oil saturation,fracture conductivity and fracture half-length are weighted higher.Finally,the horizontal wells in the target block are categorized in terms of production capacity based on the fuzzy evaluation method,and split accordingly into high-producing,relatively high-producing,medium-producing and low-producing wells.Such a categorization is intended to provide parametric guidance for reservoir fracturing and modification.
文摘The segmented water control technology for bottom water reservoirs can effectively delay the entry of bottom water and adjust the production profile.To clarify the impact of different methods on horizontal well production with different reservoir conditions and to provide theoretical support for the scientific selection of methods for bottom water reservoirs,a numerical simulation method is presented in this study,which is able to deal with wellbore reservoir coupling under screen tube,perforation,and ICD(Inflow Control Device)completion.Assuming the geological characteristics of the bottom-water conglomerate reservoir in the Triassic Formation of the Tahe Block 9 as a test case,the three aforementioned completion methods are tested to predict the transient production characteristics.The impact of completion parameters,reservoir permeability,bottom-water energy,and individual well control on the time to encounter water in horizontal wells(during a water-free production period)is discussed.A boundary chart for the selection of completion methods is introduced accordingly.The results show that the optimized ICD completion development effect for heterogeneous reservoirs is the best,followed by optimized perforation completion.Permeability is the main factor affecting the performances of completion methods,while bottom water energy and single well controlled reserves have a scarce impact.The average permeability of the reservoir is less than 500 mD,and ICD has the best water control effect.If the permeability is greater than 500 mD,the water control effect of perforation completion becomes a better option.