Extended finite element-based cohesive zone method for modeling simultaneous hydraulic fracture height growth in layered reservoirs

In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.

The forecasting efficiency under different selected regions by Pattern Informatics Method and seismic potential estimation in the North-South Seismic Zone

In 2022,four earthquakes with M_(S)≥6.0 including the Menyuan M_(S)6.9 and Luding M_(S)6.8 earthquakes occurred in the North-South Seismic Zone(NSSZ),which demonstrated high and strong seismicity.Pattern Informatics(PI)method,as an effective long and medium term earthquake forecasting method,has been applied to the strong earthquake forecasting in Chinese mainland and results have shown the positive performance.The earthquake catalog with magnitude above M_(S)3.0 since 1970 provided by China Earthquake Networks Center was employed in this study and the Receiver Operating Characteristic(ROC)method was applied to test the forecasting efficiency of the PI method in each selected region related to the North-South Seismic Zone systematically.Based on this,we selected the area with the best ROC testing result and analyzed the evolution process of the PI hotspot map reflecting the small seismic activity pattern prior to the Menyuan M_(S)6.9 and Luding M_(S)6.8 earthquakes.A“forward”forecast for the area was carried out to assess seismic risk.The study shows the following.1)PI forecasting has higher forecasting efficiency in the selected study region where the difference of seismicity in any place of the region is smaller.2)In areas with smaller differences of seismicity,the activity pattern of small earthquakes prior to the Menyuan M_(S)6.9 and Luding M_(S)6.8 earthquakes can be obtained by analyzing the spatio-temporal evolution process of the PI hotspot map.3)The hotspot evolution in and around the southern Tazang fault in the study area is similar to that prior to the strong earthquakes,which suggests the possible seismic hazard in the future.This study could provide some ideas to the seismic hazard assessment in other regions with high seismicity,such as Japan,Californi,Turkey,and Indonesia.

Drilling-based measuring method for the c-φ parameter of rock and its field application

The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(RCZ) is built. The formula for an ultimate cutting force is established based on the limit equilibrium principle. The relationship between digital drilling parameters(DDP) and the c-φ parameter(DDP-cφ formula, where c refers to the cohesion and φ refers to the internal friction angle) is derived, and the response of drilling parameters and cutting ratio to the strength parameters is analyzed. The drillingbased measuring method for the c-φ parameter of rock is constructed. The laboratory verification test is then completed, and the difference in results between the drilling test and the compression test is less than 6%. On this basis, in-situ rock drilling tests in a traffic tunnel and a coal mine roadway are carried out, and the strength parameters of the surrounding rock are effectively tested. The average difference ratio of the results is less than 11%, which verifies the effectiveness of the proposed method for obtaining the strength parameters based on digital drilling. This study provides methodological support for field testing of rock strength parameters.

A method for predicting the water-flowing fractured zone height based on an improved key stratum theory

In the process of using the original key stratum theory to predict the height of a water-flowing fractured zone(WFZ),the influence of rock strata outside the calculation range on the rock strata within the calculation range as well as the fact that the shape of the overburden deformation area will change with the excavation length are ignored.In this paper,an improved key stratum theory(IKS theory)was proposed by fixing these two shortcomings.Then,a WFZ height prediction method based on IKS theory was established and applied.First,the range of overburden involved in the analysis was determined according to the tensile stress distribution range above the goaf.Second,the key stratum in the overburden involved in the analysis was identified through IKS theory.Finally,the tendency of the WFZ to develop upward was determined by judging whether or not the identified key stratum will break.The proposed method was applied and verified in a mining case study,and the reasons for the differences in the development patterns between the WFZs in coalfields in Northwest and East China were also fully explained by this method.

Numerical investigation of geostress influence on the grouting reinforcement effectiveness of tunnel surrounding rock mass in fault fracture zones

Grouting is a widely used approach to reinforce broken surrounding rock mass during the construction of underground tunnels in fault fracture zones,and its reinforcement effectiveness is highly affected by geostress.In this study,a numerical manifold method(NMM)based simulator has been developed to examine the impact of geostress conditions on grouting reinforcement during tunnel excavation.To develop this simulator,a detection technique for identifying slurry migration channels and an improved fluid-solid coupling(FeS)framework,which considers the influence of fracture properties and geostress states,is developed and incorporated into a zero-thickness cohesive element(ZE)based NMM(Co-NMM)for simulating tunnel excavation.Additionally,to simulate coagulation of injected slurry,a bonding repair algorithm is further proposed based on the ZE model.To verify the accuracy of the proposed simulator,a series of simulations about slurry migration in single fractures and fracture networks are numerically reproduced,and the results align well with analytical and laboratory test results.Furthermore,these numerical results show that neglecting the influence of geostress condition can lead to a serious over-estimation of slurry migration range and reinforcement effectiveness.After validations,a series of simulations about tunnel grouting reinforcement and tunnel excavation in fault fracture zones with varying fracture densities under different geostress conditions are conducted.Based on these simula-tions,the influence of geostress conditions and the optimization of grouting schemes are discussed.

Selection of Water Transmission Method and Analysis of Pipe Network Zoning in Municipal Water Supply and Drainage Design

With the acceleration of urbanization,the demand for water supply and drainage pipe networks has increased significantly.In the planning of urban construction,it is necessary to optimize the design of the water supply and drainage system pipe network to effectively save energy while providing residents with more accessible water resources.Therefore,the municipal water supply and drainage system and the water transmission methods should be designed according to the geographical conditions of the city.In this paper,we mainly analyze the design of municipal water supply and drainage systems and the selection of water transmission methods.Besides,the optimization of the water supply and drainage network zoning process and pipe network maintenance is also discussed,so as to provide a reference for municipal water supply and drainage work.

A circular zone counting method of identifying a Duffing oscillator state transition and determining the critical value in weak signal detection

Identifying state transition and determining the critical value of the Duffing oscillator are crucial to indicating external signal existence and have a great influence on detection accuracy in weak signal detection. A circular zone counting （CZC） method is proposed in this paper, by combining the Duffing oscillator＇s phase trajectory feature and numerical calculation for quickly and accurately identifying state transition and determining the critical value, to realize a high- efficiency weak signal detection. Detailed model analysis and method construction of the CZC method are introduced. Numerical experiments into the reliability of the proposed CZC method compared with the maximum Lyapunov exponent （MLE） method are carried out. The CZC method is demonstrated to have better detecting ability than the MLE method, and furthermore it is simpler and clearer in calculation to extend to engineering application.

The efficiency and accuracy of probability diagram, spatial statistic and fractal methods in the identification of shear zone gold mineralization: a case study of the Saqqez gold ore district,NW Iran

In this study,geochemical anomaly separation was carried out with methods based on the distribution model,which includes probability diagram(MPD),fractal(concentration-area technique),and U-statistic methods.The main objective is to evaluate the efficiency and accuracy of the methods in separation of anomalies on the shear zone gold mineralization.For this purpose,samples were taken from the secondary lithogeochemical environment(stream sediment samples)on the gold mineralization in Saqqez,NW of Iran.Interpretation of the histograms and diagrams showed that the MPD is capable of identifying two phases of mineralization.The fractal method could separate only one phase of change based on the fractal dimension with high concentration areas of the Au element.The spatial analysis showed two mixed subpopulations after U=0 and another subpopulation with very high U values.The MPD analysis followed spatial analysis,which shows the detail of the variations.Six mineralized zones detected from local geochemical exploration results were used for validating the methods mentioned above.The MPD method was able to identify the anomalous areas higher than 90%,whereas the two other methods identified 60%(maximum)of the anomalous areas.The raw data without any estimation for the concentration was used by the MPD method using aminimum of calculations to determine the threshold values.Therefore,the MPD method is more robust than the other methods.The spatial analysis identified the detail soft hegeological and mineralization events that were affected in the study area.MPD is recommended as the best,and the spatial U-analysis is the next reliable method to be used.The fractal method could show more detail of the events and variations in the area with asymmetrical grid net and a higher density of sampling or at the detailed exploration stage.

A Coastal Zone Segmentation Variational Model and Its Accelerated ADMM Method

Effective and efficient SAR image segmentation has a significant role in coastal zone interpretation. In this paper, a coastal zone segmentation model is proposed based on Potts model. By introducing edge self-adaption parameter and modifying noisy data term, the proposed variational model provides a good solution for the coastal zone SAR image with common characteristics of inherent speckle noise and complicated geometrical details. However, the proposed model is difficult to solve due to to its nonlinear, non-convex and non-smooth characteristics. Followed by curve evolution theory and operator splitting method, the minimization problem is reformulated as a constrained minimization problem. A fast alternating minimization iterative scheme is designed to implement coastal zone segmentation. Finally, various two-stage and multiphase experimental results illustrate the advantage of the proposed segmentation model, and indicate the high computation efficiency of designed numerical approximation algorithm.

Vibration analysis of steel structures including the effect of panel zone flexibility based on the energy method

New approximate formulas are proposed to determine the natural frequencies of structures considering the effects of panel zone flexibility and soil-structure interaction. Several structures with various earthquake resisting systems are idealized as prismatic cantilever flexural-shear beams. Floor masses are considered as lumped masses at each story level and masses of columns are evenly distributed along the cantilever beam. Soil-structure interaction is considered as axial and rotational springs, whose potential energy are formulated and incorporated into overall potential energy of the structure. Subsequently, natural frequency equations are derived on the basis of energy conservation principle. The effect of axial forces on natural frequency is also considered in the proposed formulas. Using the method presented in this study, natural frequencies are computed using a simplified method with no complex numerical modeling. The proposed formulas are verified via experimental and numerical methods. Close agreement between the results from these three approaches are observed. Furthermore, the effects of panel zone flexibility, continuity plates and doubler plates on the natural frequencies of buildings are investigated.

Sintering zone prediction in direct metal laser sintering by finite element method

A three-dimensional finite element thermal model in direct metal laser sintering(DMLS) including the effect of powder-to-solid transition were established to predict sintering zone, which benefited the determination of suitable process parameters in DMLS. The nonlinear transient model of the metals thermal conductivity for powder-to-solid transition was developed. The model uses solid thermal properties of material in liquid-phase zone, transitional ones in sintering or sintered zone and powder ones in unsintered zones of powder bed to predict, respectively. Sintering zone boundary was estimated by maximum temperature history profile. Experiments were carried out using multi-component Cu-based metal powder. Compared experimental and predicted results, the mean error of sintering depth and width are 7.8% and 14.4%, respectively, which confirms the accuracy of the FEM prediction.

Study of the segmentation of active fault by the boundary element method──analysis of the Xianshuihe fault zone

On the basis of the recent geological surveys and the relevant studies of the Xianshuihe fault zone, this paper analyzes the seismogenic mechanism of some faults and characteristic morphology on the fault zone by the boundary element method and discusses the fault segmentation and the related distribution of the earthquake ruptures. The main conclusions are: For the first order segmentation, the Xianshuihe fault zone can be divided into three major segments (the northwestern Luhuo-Qianning segment, the middle linking fracture region and the southeastern Kangding segment). Among them, the differences are not only in geometry and structure, but also in mechanical property and dynamics. Some of the characteristic morphologies on the Xianshuihe fault zone reflect the effects in cumulative deformation due to long-term fault movement, and reveal the fault segmentation in different orders. Such morphologies control, to some extent, the developments and the distributions of the earthquake ruptures on the fault zone.

Computation of Protection Zone of a Lightning Rod Using Method of Moments and Monte Carlo Integration Technique

An accurate determination of lightning protection zone is an important issue in the analysis and design of an appropri-ate lightning protection system. This paper presents a fast and accurate protection zone determination methodology for metallic lightning rod. The methodology is based on Quasi Monte Carlo Integration technique applied to Method of Moments (MoM) solution of Integral Equations. As an example, solution of the integral equation for unknown charge distribution on lightning rod is obtained. The electric field in the region surrounding the rod is then computed and the protection zone plotted accordingly. The effect of the thickness of the rod on the protection zone is also studied.

Radial integral boundary element method for simulating phase change problem with mushy zone

A radial integral boundary element method(BEM)is used to simulate the phase change problem with a mushy zone in this paper.Three phases,including the solid phase,the liquid phase,and the mushy zone,are considered in the phase change problem.First,according to the continuity conditions of temperature and its gradient on the liquid-mushy interface,the mushy zone and the liquid phase in the simulation can be considered as a whole part,namely,the non-solid phase,and the change of latent heat is approximated by heat source which is dependent on temperature.Then,the precise integration BEM is used to obtain the differential equations in the solid phase zone and the non-solid phase zone,respectively.Moreover,an iterative predictor-corrector precise integration method(PIM)is needed to solve the differential equations and obtain the temperature field and the heat flux on the boundary.According to an energy balance equation and the velocity of the interface between the solid phase and the mushy zone,the front-tracking method is used to track the move of the interface.The interface between the liquid phase and the mushy zone is obtained by interpolation of the temperature field.Finally,four numerical examples are provided to assess the performance of the proposed numerical method.

Single-borehole measuring method for broken rock zone in gently inclined thin layer weakness structure

According to the structural characteristics of gently inclined thin layer rock mass in which lots of weak interlayer existed,the concept of gently inclined thin layer weakness structure was proposed.If single-borehole measuring method of the acoustic along the conventional arrangement mode was used in measuring the broken rock zone in this structure,the change of the relationship curves (Vp-L) between acoustic p-wave velocity (Vp) and borehole depth (L) would present the irregular feature due to the mechanical characteristics of layered rock mass and harmful effects of weak interlayers,and the scope of broken rock zone couldn't be defined quickly.Based on the analysis of the me- chanical characteristics of layered rock mass,the propagation rule of acoustic and distri- butions characteristics of plastic zone and slip zone in layered rock mass,new arrange- ment mode of acoustic measuring boreholes for broken rock zone in gently inclined thin layer weakness structure was proposed.Namely,the measuring boreholes in two sides were parallel to the strata,the measuring boreholes in the roof and floor perpendicular to the strata.Besides the controlling depth of the measuring boreholes in the scope of the large plastic zones or the large slip zones should be increased.Engineering example showed that new acoustic measuring boreholes arrangement mode had the better appli- cability and could determine the scope of the broken rock zone in the gently inclined thin layer weakness structure quickly.

一种基于“四区”的低渗透油藏CO_(2)埋存量计算方法及应用

CO_(2)驱油和埋存能有效减少温室气体排放量达到碳中和目标,已有的CO_(2)埋存量计算方法主要针对CO_(2)的静态埋存量进行粗略计算,未考虑实际生产过程中CO_(2)埋存量的变化。针对上述问题,运用CO_(2)溶解、CO_(2)波及体积和驱油机理,将CO_(2)驱油与埋存过程分为气相区、两相或近混相区、扩散区和油相区,并基于“四区法”计算CO_(2)埋存量,得到了不同烃类的注入孔隙体积倍数、注入压力、注气速度下的CO_(2)动态埋存量的变化规律。将研究成果应用于W油田低渗储层,结果表明:注入烃类孔隙体积倍数、压力、注气速度与总埋存量呈正相关性,当压力由12 MPa升至30 MPa,CO_(2)埋存总量增加15.53×10^(4)t;当注气速度由5 t/d增加至30 t/d,峰值CO_(2)埋存总量由3.51×10^(4)t提高至12.62×10^(4)t。研究成果可为同类油藏开展CO_(2)驱油与埋存项目提供新的思路。

覆岩导水裂隙带发育高度动态演化规律研究

为了研究覆岩导水裂隙带发育高度动态演化规律,以小保当矿区2^(-2)煤层为研究对象,通过理论分析、相似模拟实验及实例验证的方法,对导水裂隙带动态发育高度进行研究;借助概率积分法预计其上部岩层移动变形的理论公式,给出了1种基于覆岩曲率变形的导水裂隙带动态发育高度预计方法。研究表明:导水裂隙带上部岩层下沉系数是1个与挠度及下部自由空间高度相关的分段函数,上部岩层曲率变形大小是决定导水裂隙带发育高度的关键因素;导水裂隙带发育高度与工作面推进长度相关;理论预计小保当矿区2^(-2)煤层导水裂隙带发育高度为160.8 m,相似模拟实验发育高度为155 m,现场实测发育高度为152.01~175.57 m。

掏土孔应力扰动法下的塑性区范围研究

针对普通浅层掏土法下倾斜建筑回倾速度慢、沉降不均匀的问题,提出掏土孔应力扰动法。在岩石力学的弹性圈和塑性圈应力解中引入修正系数,得到掏土孔应力扰动法下的塑性区半径计算公式。在进行现场掏土孔周应力扰动试验时,采用地锚反压的形式模拟上部建筑荷载,利用长钉在掏土孔周围进行深度为50 mm的螺旋扰动,通过测量应力扰动前后的孔周土体应力并结合有限元模拟结果,得到实际塑性区范围,验证了塑性区半径公式的正确性。结果表明:当掏土孔半径为100 mm,应力扰动范围为孔径的1/2时,塑性区范围扩大30%~40%,掏土孔应力扰动法对于扩大塑性区范围、加快建筑纠倾速度具有工程应用价值。

太原盆地交城断裂带北段工程避让范围研究

在系统分析和整理交城断裂带北段探测成果的基础上,采用地质探槽剖面法、经验统计法及数值模拟法,预测评估交城断裂带北段在发生最大潜在震级M_(S)7.2地震时的同震地表破裂带宽度,探讨不同方法结果的精度和可靠性,并提出采用综合加权方法确定交城断裂带北段地表破裂带宽度。在此基础上综合考虑探测方法精度、断层上下盘效应等影响因素,确定交城断裂带北段工程避让区为F_(1)^(1)迹线西侧外延15.0m至F_(1)^(2)迹线东侧外延86.9m范围。该结果对交城断裂带沿线的国土规划利用、建(构)筑抗震设计具有一定的指导意义,该方法可为同类研究提供借鉴。

大凉山地区不同马铃薯品种产量和营养品质的综合评价

为筛选出适宜四川大凉山地区种植的马铃薯品种,对14个马铃薯品种的12个产量与营养品质指标进行测定,并结合主成分分析、隶属函数法、系统聚类分析对马铃薯产量和品质表现进行综合评价。结果表明,单株结薯数、单薯鲜质量、单株块茎鲜质量、单产、还原糖含量的变异系数均超过30%;单株块茎鲜质量与单薯鲜质量、单产和蛋白质含量呈极显著正相关,干物质含量与淀粉含量、蛋白质含量呈极显著正相关,单薯鲜质量与VC含量、锌含量呈显著负相关,蛋白质含量与钾含量呈极显著负相关;主成分分析结果表明,12个指标可用4个主成分来表示,方差累积贡献率达到86.040%。进一步采用隶属函数法和系统聚类分析将14个品种分为3类,筛选出6个综合表现较优的品种,分别为川凉薯10号、青薯9号、川芋50、川凉芋13、云薯108、川芋22号。