更全的杂志信息网

Propagation and aperture of staged hydraulic fractures in unconventional resources in toughness-dominated regimes

更新时间:2016-07-05

1.Introduction

In unconventional completion,single-stage fracture treatments evolved to multistage stimulation treatments and fracturing of standalone wells was progressed to simultaneous fracturing of multilateral wells in order to increase reserves per well,enhance well productivity,and improve project economics(King et al.,2008).Some researchers mentioned that an optimization is required for such completion techniques.For example,it was reported that increase in spacing between the fractures induces less interference betweenpropagating fractures and hence requires lessbreakdown pressureto initiateafracture(Singh and Miskimins,2010).In addition,no deviation or collapse of hydraulic fractures occurs under this circumstance.It is the reason of considering stress shadow concept in designing of hydraulic fracturing pattern for this kind of reservoirs.Stress shadowconcept has been discussed by many researchers(e.g.Roussel and Sharma,2011;Morrill and Miskimins,2012;Taghichian et al.,2014).

The main purpose of hydraulic fracturing in unconventional reservoirs is to create parallel hydraulic fractures perpendicular to the horizontal wellbore and make a connected reservoir to the wellbore.In order to reach this goal,observing whether a hydraulic fracture propagates or is trapped is of paramount importance.In addition,the direction of propagation of hydraulic fractures plays a decisive role in having ideal straight hydraulic fractures with no deviation and collapse.Therefore,a fracturing pattern should be designed in such a way that simultaneously satisf i es both of the above-mentioned conditions,i.e.reasonable fracturing pressure(via stress intensity factor(SIF))and controlling the fracturing direction(via stress shadow effect(SSE)).

Propagation of hydraulic fractures depends on stress field around the tip.Hence,many researchers have tried to calculate the stress field around cracks which is influenced by crack geometry,fracturing scheme,and applied boundary conditions.In order to simplify the stress calculation close to the tip,a new term called SIF was defined.The generalized form of stress field equation for modes I,II and III of loading,valid in the vicinity of the crack tip,is given by

whereσijis the stress in the vicinity of the crack tip;f,g and h are the trigonometric functions;r andθare the cylindrical coordinate components;and KI,KIIand KIIIare the modes I,II and III SIFs,respectively.The concept of SIF was first proposed by Irwin(1957),by which the effects of geometry and boundary conditions are separated from spatial location of stress analysis.

1.不够主动,有部分学生对纠错的认识不够深刻,没能把它看作是一个自我检查、自我巩固与自我提高的过程,很多时候认为它可有可无,甚至是一个在教师强迫下的任务,在这样被动的心境下,纠错往往是应付了事;

基于德尔菲法的医院药品不良反应快速上报信息化系统评价指标体系的构建 …………………………… 温润龙等(13):1742

Assuming a planar crack in the x-z plane,one can have the SIFs as

One conspicuous point herein is the dependence of KIIand KIII on the shear stresses.When crack propagating,the angle of propagation in two dimensions depends on KIand KIIwhich can be obtained via the following relationship(Stone and Babuska,1998):

whereΘis the change angle of propagation direction.It is induced from Eq.(3)that no change in the propagation direction is observed in case of having no KII,which corresponds with the case of having no shear stress at the tip(see Eq.(2b)).Therefore,it can be induced that having straight hydraulic fractures with no deviation or collapse needs negligible shear stresses at the tip.It is noted that this is an ideal condition which is assumed for the current analyses in order to obtain the ideal distance between hydraulic fractures propagating merely under the circumstance of mode I.Therefore,under the above-mentioned condition,the only required term to determine stress field around the crack tip is KI,which changes with geometries of the crack,medium,boundary and loading conditions(Broek,1982).

In this way,having the SIF behavior of hydraulic fractures,one can judge how/when a hydraulic fracture propagates in the medium.Many two-dimensional(2D)crack geometries have been analytically solved and their SIFs have been reported in the literature(e.g.Tada et al.,2000).However,there have also been some problems in which the geometries of cracks/medium were challenging and stress field for these fractures was not easy to be analytically determined.The SIF of such problems was defined by utilizing boundary element(BE)and finite element(FE)methods(e.g.Sih,1973;Murakami,1987;Tada et al.,2000).Threedimensional(3D)cracks with different geometries,such as embedded cracks in an infinitely extended homogeneous,isotropic solid medium,opened up due to prescribed internal pressure,have also been analyticallysolved bya number of investigators(e.g.Keer,1964;Sneddon and Lowengrub,1969;Shah and Kobayashi,1971;Guidera and Lardner,1975).Mastrojannis et al.(1979)also developed a method for determination of SIF for a general-shaped crack with internal pressure in an infinite medium utilizing numerical integration.Furthermore,using the 2D Fourier transform method,Kassir(1981)succeeded in solving the SIFs around rectangular cracks.Nejati et al.(2015)also proposed a novel domain integral approach for SIF calculation of 3D cracks with tetrahedral elements and not requiring any structured mesh.It is worth mentioning that SIF of 3D cracks depends on spatial location around the crack edge.For instance,for an internally pressurized rectangular crack,SIF along the length is higher than that along the width.According to Kassir and Sih(1966),a basic characteristic of any 3D crack problem is the fact that the state of stress in a normal plane near a smooth crack front is essentially a plane-strain one.Therefore,for a rectangular crack internally pressurized,SIF is the highest along its length and it can be determined as

where PHis the internal pressure;a is the half-length of the fracture;AR is the fracture aspect ratio;and the coefficients f1,f2and f3 are defined as 0.5415,1.8086 and 0.6943,respectively.

我国高校的师生数量都较大,上文指出的创新培养类的大数据信息都是由不同职能部门分散采集,如课程类数据多由教务部门采集、图书借阅类数据储存在图书馆、毕业以及第二课堂数据归学工和团委负责。这里的每个部门都有着自己的信息管理系统,基本都是供各个部门内部使用,信息系统之间没有实现信息数据共享和互联互通,也因此极大地制约了数据的整合与共享。

For almost all of the 3D crack problems,with respect to analytical solutions,SIFs were proposed for single fractures.Therefore,there seems to be lack of enough knowledge about the SIF values for the cases where multiple fractures exist in the medium.In hydraulic fracturing of unconventional reservoirs,due to the existence of many 3D fractures in the medium,the SIF change of multiple fractures placed between parallel lateral wells should be studied.Any fracturing scenario having a higher ratio of SIF with respect to the case of a single fracture in a standalone well can be considered as a fracturing technique with higher propagation potential in the target zone.Study of SIF behavior for different fracturing scenarios can be considered as a simultaneous tool required foroptimization ofhydraulic fracturestogetherwith SSE.Employing these two tools,the spacing and generally the pattern of the fracturing can be designed.

Table 1 Input variable range for the numerical simulation.

Note:c is the half length of the fracture.

Fracture aspect ratio Fracture spacing Fracture distance Fracture offset 0.2-1 (0.1-7)c (0.1-3)c (0.1-1)c

Regarding the SIF determination of 3D fractures,it is a challenging procedure to obtain acceptable values of SIF using FE method,since a very specific modeling procedure is required for this purpose.For the large number of numerical simulations required for different fracturing patterns(300 models;see Table 1),we tried to use an efficient way for SIFanalysis.Therefore,the ratios of stresses in elements around the crack tip with the same sizes and shapes at equal distances close to the tip are calculated for different scenarios of fracturing compared to that of a single hydraulic fracture.Using this method,any increase or decrease in SIF for any proposed fracturing scenario is reported compared to the case of a single fracture without dealing with the absolute values of SIFs.Considering the results of the current analyses,engineers can investigate whether or not the SIFs of the fractures in a fracturing scenario is changed with respect to a single fracture in an infinite domain qualitatively and quantitatively.As a result,any fracturing scenario can be checked using this methodology in addition to the SSE in order to understand whether there is any fracture trapped or deviated from its straight path.

2.Methodology and conventions

In this work,in order to study the effect of fracture geometry and different fracturing patterns on the SIF of the fractures,different aspect ratios,and multistage and simultaneous fracturing schemes are simulated using FE-based software,ABAQUS CAE 6.12.Modeling of 3D hydraulic fractures is done by assuming hydraulic fractures contained inside the target zone,perpendicular to the wellbore.Relative placement of the hydraulic fracture with the target zone and wellbore,and its geometry definition are shown in Fig. 1,in addition to the hydraulic fracture with the SIFs under study.SIF along the height of a hydraulic fracture is called SIFH which causes propagation in the horizontal direction,while SIF along the length of a hydraulic fracture is called SIFLwhich causes propagation in the vertical direction.

In addition to the SIF change,aperture of hydraulic fractures mayalso be influenced by fracturing pattern.According to Sneddon and Elliot(1946),aperture of a hydraulic fracture in toughnessdominated regime along the length can be represented by an elliptical function as

(8)All the modeling has been performed in three dimensions.

where νis the Poisson’s ratio,Enis the Young’s modulus of net play,wx)is the half aperture,and wmaxis the maximum half aperture(located at the center of the fracture(x=0)).From Eq.(5a),we observe that the displacement of crack tips(at x=c)is zero and the displacement of the edges increases from the tips to the center of the fracture where the maximum displacement occurs.When the maximum value of half aperture is known,by substituting it into Eq.(5a),the aperture values along the length and height can be determined.Therefore,having the maximum half aperture located in the center of the fracture,one can have the aperture distribution on the whole fracture surface.In this study,the term aperture change is the ratioof this maximumhalf aperturefordifferent cases with respect to the single fracture in a standalone well.Fig. 2 shows the scenarios considered for different hydraulic fracturing patterns.

《人民日报》2017年12月的一篇报道提到,东北某乡党委书记反映,有一天一早上要接待13个检查团,只好让乡里干部全部上阵。西部某市街道办张主任也表示,往往前脚刚把这个检查组送走,那个评比团眼看着就来了。每项检查都得准备完整资料,大量复印、装订的重复劳动耗费人力物力财力。去年仅迎接检查准备的纸质台账就装了100多盒;但上级检查往往走马观花,翻一下就完事了,部分材料还在仓库堆着。

It is evident in Fig. 2 that four different scenarios have been considered for investigation of propagation potential and aperture of hydraulic fractures in the target zone.Scenario 1 is the basic scenario inwhich we have a single stage fracture from a standalone well.Scenario 2 shows multistage hydraulic fracturing in a standalone well.Scenario 3 considers the effect of simultaneous single stage fracturing of the medium between two parallel wells.Finally,in Scenario 4,both of the fracturing strategies,i.e.simultaneous and multistage hydraulic fracturing,are considered.In order to study the change of propagation potential and aperture variation of hydraulic fractures as a result of fracturing pattern,SIFs and apertures of fractures in Scenarios 2-4 arecomparedwith those in Scenario 1.

Fig. 1.(a)Position of the hydraulic fracture with respect to the target zone and wellbore associated with its geometry definition,and(b)Propagation direction according to the SIFs along the length and height of the fracture.

Fig. 2.Different fracturing techniques/patterns in this study.Lpis the half-spacing between multistage fractures,and Lsis the distance between fracture tips.

Changing the fracture geometry by its aspect ratio,the spacing between fractures in multistage fracturing,and the distance between the meeting fractures in simultaneous fracturing scheme,one can study the changes of SIF and aperture of the fractures compared to Scenario 1 which is the basic scenario as a single fracture in an infinite medium.Using this method,the effect of these dimensional variables on the SIFand the aperture of fractures can be quantif i ed.Moreover,the qualitative study for the effect of offset between the meeting fractures on the SIF change is conducted as well.The term offset means the distance between two planes in which hydraulic fractures are located in a simultaneous fracturing scheme.The propagation potential and aperture of hydraulicfracturesare investigated considering the following assumptions:

据了解,今年以来江岸区食药监局为迎接军运会的召开,积极加强辖区餐饮食品安全的规范管理,切实完善餐饮食品安全长效监督机制,积极联合辖区街道在全区建成10条餐饮服务示范街,重点对小餐饮店的出店经营、内部布局、硬件设施等进行整体提档升级,通过示范引领、以点带面、推动餐饮店完善设施、健全制度、规范管理,促进餐饮店向规范化、标准化、特色化和品牌化方向发展,全面实现餐饮店持证亮证经营,落实食品安全管理制度,前厅后厨合理布局,消毒制度落实,“三防”设施配备,达到餐饮服务食品安全年度量化等级提升的目标,满足市民健康安全饮食需求。

(1)Constant pressure insidehydraulicfractures has been assumed in all the numerical models.

九点刚过,丁小强接了一个电话。是布雅兰的电话。丁小强告诉布雅兰晚上不回去了,因为喝多了。当时,他正和杜一朵打牌。他输得只剩下一条裤衩了。看着丁肚腩下的小布包,杜一朵说:“这就是你勤政廉洁的成果吗?”杜一朵还说:“就这么一点小成果,还好意思勾引美女?”

(7)The geometry of the hydraulic fracture has been assumed as a square(aspect ratio of unity)and as rectangles(aspect ratios less than unity).

综上所述,核磁共振诊断膝关节半月板损伤的效果显著,通过膝关节半月板形态、信号改变等判断损伤程度,为临床诊治提供依据。

(3)Modeling has been done in a completely elastic medium without considering any plasticity constitutive law.

(4)Hydraulic fracture has been assumed as a stationary crack without considering any propagation.

(5)The considered hydraulic fracturing regimes are assumed as toughness-dominated rather than viscosity-dominated(see Detournay,2004;Bao et al.,2017).

(6)The aspect ratio is defined as the ratio of height to length of the crack and is always equal to or less than unity in this study.

It can be seen from Fig. 6 that the SIF theoretically decreases to 0.21 of its original value when the spacing between the fractures is as low as 0.17c.This means that too closely-spaced fractures causethe SIF to decrease significantly.Therefore,this fracturing pattern may not be suitable and thus a much higher fracturing pressure is required for the fracture to propagate.In addition,when this fracture is located in the shadow region of the first fracture,it may deviate from its original path.It is also seen that the aspect ratio of the hydraulic fracture plays an important role in SIF change of the multistage fractures.In fact,going back to the original state(SIF(H,L)M/SIF(H,L)O→1)occurs in a shorter spacing for fractures with lower aspect ratios.Moreover,comparing the changes of SIFH and SIFLof hydraulic fractures,it is induced that SIFLis more influenced by multistage fracturing than SIFH(with AR convention;see Section 3).This causes the propagation potential to decrease more in vertical directions compared to that in horizontal directions,which results the fracture to stay in the net play and be contained rather than moving inside the bounding layers.

(2)In modeling of media containing more than one hydraulic fracture,similar hydraulicpressuresand aspectratios together with the same mesh type and size have been assumed for all the hydraulic fractures.

地西泮(Diazepam纯度99.4%)、甲硝唑(Metronidazole,MNZ纯度99.9%)、地美硝唑(Dimetridazole,DMZ纯度99.9%)、洛硝哒唑(Ronidazole纯度99%)、羟甲基甲硝咪唑(HMMNI纯度99%)和羟基甲硝唑(MNZOH纯度99.4%)。

在新媒体环境下,民族高校图书馆如何依托新媒体技术,通过信息服务理念的更新和服务模式的拓展,不断改善信息服务质量,依据民族高校读者特色进一步提升信息服务能力,成为亟需解决的问题。基于此,笔者细致深入探讨了新媒体时代信息服务的特点,针对民族高校图书馆信息服务工作开展的障碍以及具体开展的措施进行了一些探讨,以期为民族高校图书馆利用新媒体开展信息服务工提供参考依据。

(9)In the scenario of having simultaneous fractures,the distance between fracture tips is controlled by changing the distance between the wells but not the fracture geometry or its aspect ratio.

(10)Symmetry has been assumed for fractures for more efficient simulation.

In this study,a range of fracture spacing and distance has been assumed by which interaction between neighboring fractures is observed for all the simulations.Input variables and their ranges are listed in Table 1:five different aspect ratios with interval of 0.2,ten different fracture spacings in multistage fracturing,and six different fracture distances in simultaneous fracturing.

(11)Due to the low permeability of the reservoir,fluid diffusion time-scale is assumed much longer than the fracturing timescale and only undrained response is considered herein.

In this section,two horizontal wells are assumed to be placed parallel to each other and only one hydraulic fracture for each well is considered(Scenario 3;Fig. 2).It is also assumed that both of the hydraulic fracture faces align on a single plane without any offset.The distance between fracture tips is changed by assuming different distances between the two lateral wells.It was observed that the meeting edges are influenced by each other and other edges do notshowanychange.This means thatas theverticaledges(height)of simultaneous fractures are meeting each other,SIFHis much affected(increased)by the simultaneous fracturing.This is the reason that SIFLdoes not show significant change by this fracturing pattern.Therefore,SIFH ofsimultaneousfracturesis compared to that of a single stage fracture from a standalone well(Scenario 1;Fig. 2).The results of SIFHchange are shown in Fig. 7 for different aspect ratios.As it can be seen from Fig. 7,SIFHof simultaneous fractures is controlled by two key variables:distance between the tips and aspect ratio.

The methodology of estimating the SIF ratio of any hydraulic fracturing scheme over that of a single hydraulic fracture in an infinitedomain is explained herein.It is known that stressat a point in the vicinity of the crack tip in the plane where the fracture is located inside(r=R,θ=0?)is related to the SIF using the following equation:

where R is the distance from tip;KIis linearly related to the internal pressure(PH),which is directly related to a function of the geometry of the fracture,and a function incorporating the effect of boundaries on the fracture.In this way,having the ratio of stress for a fracturing scenario under study over that of the basic Scenario 1,one can have the ratio of the SIF for the fracturing technique over that of the basic Scenario 1.Therefore,we have

Fig. 3.(a)Numerical mesh and(b)model configuration for simulation of one-eighth of a hydraulic fracture.

Fig. 4.Stress variations along the fracture edges for AR=1 and 0.2.

Eq.(7)shows the direct relationship between the ratio of stresses along fracture edge and the ratio of SIF for the scenario under study over that of the basic Scenario 1.Therefore,stresses along the length and height of the fracture,showing an exponential behavior(see Fig. 4),can be indicative of SIF of the fracture edge.One important point in this regard is the singularity of stress on the crack edge which means that the comparison of stresses should be done along a structured mesh with the same type and size.Only under this circumstance,the stresses are comparable due to the precision of stress calculation and the same distance to the tip for both of the scenarios.The other important point is that SIF changes along the length and height of the fracture.This means that the SIF of a rectangular crack along every point on the half-length and halfheight of the fractures is different.Hence,in order to have a representing value of the SIF for the length or height of a fracture,a fi tting function is used to be fi tted on the stresses along the length of the fracture(σzz-x)and the other on stresses along the height of the fracture(σzz-y).The ratio of fi tting function coefficients for the designed scenario over a single fracture is used as a representative of SIF change(see solid lines in Fig. 4).The proposed functions,which give satisfactory fi t on the numerical stress values,are given by

测试采用的图1所示的带有Rogowski线圈(以下简称罗克CT)的限流电抗器,参数取优化设计以后的方案。测试得到的电流波形如图4所示。

where σzz0is the normal stress at the corners;(σzzLand (σzzHare the normal stresses along the length and height of the fracture,respectively;x and y axes originate from the corner toward the fracture length and height;b is the half height of fracture;and aHi and aLi(i=1 and 2)are the coefficients of the function.Observing the behavior of the proposed function,it is revealed that stress change is controlled by a1rather than a2.This is because the exponential part is directly multiplied by aH1and aL1,while aH2 and aL2are only the components of the exponential part with negative sign.Therefore,aH1and aL1were considered as sufficient variables for showing the SIF change along the edge of a rectangular hydraulic fracture.Fig. 4 shows a typical example of normal stress change along the crack edge.Considering the representative coefficient of the function(aL1,aH1)proposed in Eqs.(8a)and(8b),one can make a good comparison between the SIF of different fracturing scenarios and the case of a single fracture in an infinite domain.

Fig. 5.Stress validation in the vertical direction away from the fracture center.(a)Plane-strain crack,and(b)Penny-shaped crack.

Fig. 6.SIF changes along the(a)height and(b)length of the fracture as a result of having multistage hydraulic fractures from a standalone well.

It is evident from Fig. 4(solid lines are predictions and dotted lines are numerical values from the software)that the behavior of stress along crack length or height has satisfactorily been predicted using the proposed Eq.(8).It is also seen that normal stress is the highest in the middle point of the crack length and it is reduced to the minimum value at the crack corner.In addition,stress change along the height and length of fractures with aspect ratio of unity is similar as expected(aL1=aH1),but for cracks with aspect ratios less than unity,stress decreases in both the length and height of the fracture.The magnitude of decrease is more severe in the fracture height compared to that in the length.Considering Eq.(6),one can see that higher stress means higher SIF which leads to higher propagation potential.Therefore,as it is seen in Fig. 4,hydraulic fractures with aspect ratios less than unity have more propagation potential in up/down directions than that in left/right directions(SIFLSIFH).

3.verification of numerical simulation with analytical solutions

In order to perform a valid numerical simulation,it is first required to numerically model some specific crack problems using the proposed methodology and make a comparison between available analytical solutions and the obtained numerical simulations.Any crack can be specif i ed in space by three geometrical parameters:length,width(called aperture),and height.The ratio of any of these three parameters over the other can be called aspect ratio of the fracture.In this paper,the ratio of height to length is called aspect ratio(AR=height/length).In this section,two cracks with different geometries,i.e.a plane-strain crack with infinite height(AR→+∞),and a penny-shaped crack(AR=1),both with internal pressure,are numerically simulated and the induced stresses around the crack are compared with those of the analytical solutions.Sneddon and Elliot(1946)reported the analytical solution for stress around a crack in an infinite 2D medium with internal pressure.The formulation for penny-shaped cracks was also derived by Sneddon(1946).

The validation approach described here compares the horizontal and vertical stresses around the simulated crack with the analytical results presented by Sneddon and Elliot(1946).Fig. 5 shows the horizontal and vertical stresses along a line perpendicular to the face of the hydraulic fracture.Fig. 5 demonstrates good agreement between the analytical solutions and numerical results for plane-strain and penny-shaped cracks.Therefore,the numerical modeling strategy adopted here can be used for modeling of hydraulic fractures.

4.Effect of multistage fracturing on SIF and aperture in standalone wells

Due to the fact that multistage fracturing is one of the effective techniques for fracturing of unconventional reservoirs,in this section,we investigated the effect of adjacent hydraulic fractures from a standalone well(Scenario 2;Fig. 2)compared to the case of having a single fracture in an infinite medium(Scenario 1;Fig. 2)on the SIF of each scenario.To do so,parallel hydraulic fractures with varying spacing and aspect ratio are assumed,and the SIF change with all the available configurations is studied along the length and height of the fractures.Fig. 6 shows the behaviors of SIFHand SIFLof the hydraulic fractures of multistage fracturing technique with respect to a single stage fracture,in which SIF(H,L)Mand SIF(H,L)O stand for SIFs for the height and length of multistage fractures and the single stage fractures in an infinite medium,respectively.

党的十九大报告强调:“党的基层组织是确保党的路线方针政策和决策部署贯彻落实的基础。新时代医院党支部工作创新离不开“互联网+”的运用,医院党支部工作从完善制度体制、加强顶层设计、搭建“互联网+”党建平台,培养人才等方面来利用“互联网+”,进一步发挥基层党支部的战斗堡垒作用。

Table 2 Coefficients of the function predicting SIF reduction with distance between multistage fractures.

Usage Coefficients AR=1 AR=0.8 AR=0.6 AR=0.4 AR=0.2 RSIFHM m1H 0.16447 0.16582 0.16785 0.12873 -0.2065 m2H 0.42449 0.31392 0.20294 0.10429 0.0209 m3H 1.0168 1.01308 1.0082 1.00309 1.00179 m4H 1.92373 1.93651 1.94993 1.87369 1.85405 RSIFLM m1L 0.16447 0.17497 0.18725 0.19689 0.19505 m2L 0.42449 0.37417 0.29515 0.17889 0.05834 m3L 1.0168 1.01373 1.01024 1.00477 1.00275 m4L 1.92373 2.02108 2.12896 2.22962 2.28194 Rw M m1w 0.09109 0.09643 0.10355 0.10695 0.17979 m2w 0.46656 0.35414 0.23938 0.13539 0.03863 m3w 1.00792 1.00662 1.00542 1.00338 1.00603 m4w 2.41794 2.42928 2.40773 2.29029 2.35485

Fig. 7.SIFHratio of simultaneous fractures with respect to standalone fractures.

Table 3 Function coefficients for predicting SIFHratio of simultaneous over standalone fracturing.

Aspect ratio Function coefficient,n 1 0.08337 0.8 0.07161 0.6 0.05899 0.4 0.04337 0.2 0.03642

In order to have a quantif i ed SIF change in multistage fracturing technique,based on the behavior of the numerical results and the proposed technique for SIF ratio calculation,the following equation is proposed:

由于铸件结构较复杂,生产工艺参数若把握不精确,会造成铸件表面缺陷较多,主要出现浇注不足、裂纹和变形等缺陷.缺陷主要出现在铸件加强筋处,通过切取铸件缺陷易产生的同一部位本体试样,来研究工艺参数对其组织和性能的影响.

5.Effect of simultaneous fracturing on SIF in parallel wells

Table 1 shows that a number of numerical models(300 models for the quantitative simulation)are required to be built and analyzed to predict the SIF changes and aperture variation in different fracturing scenarios.Fig. 3 shows model geometry,assigned mesh,and symmetry that have been used.It is evident thatonly one-eighth of the hydraulic fracture has been simulated in three dimensions.

Fig. 8.Effects of distance(Ls)and offset(Lo)between fracture tips on the SIF of the fractures.(a)AR=1,and(b)AR=0.4.

Fig. 9.Effects of multistage simultaneous fracturing on(a,b)SIFHand(c,d)SIFLof hydraulic fractures.

In addition,it is observed that aspect ratio plays an important role in this SIFHincrease in such a way that higher aspect ratio results in higher SIFHincrease.The reason for this observation is that the region of inf l uencing stress field is larger for the case of higher aspect ratios.Based on the behavior of the observed numerical results,this increasing effect can be quantif i ed by fracture aspect ratio and distance between the tips.In order to have a quantifying equation for SIF ratio of simultaneous fractures(Scenario 3)versus single stage fractures(Scenario 1),the following relationship is proposed:

where is the ratio determining the increasing effect ofthe meeting tips,and n is the coefficient varying with aspect ratio of the fractures.coefficient n can be determined considering Table 3.

Of course,all these values are only valid for single stage simultaneous fractures without having any multistage ones.Considering the numerical results,it is observed that there is no aperture change in simultaneous fracturing scenario compared to the case ofa single fracture.This means that any change in the aperture is merely due to the multistage fracturing.

Table 4 Coefficients of the predicting function for SIFHchange as a result of simultaneous multistage fracturing.

Tip distance,Ls/c CoefficientsAR=1 AR=0.8AR=0.6AR=0.4AR=0.2 0.25 q1H 0.172010.18063 0.18685 0.20748 0.10236 q2H 0.900470.64194 0.39699 0.19147 0.05582 q3H 1.341131.27754 1.20242 1.12411 1.03763 q4H 1.982341.98857 1.96326 1.93938 1.70901 0.5 q1H 0.189980.19526 0.19393 0.19476 -0.24732 q2H 0.805020.57514 0.35572 0.16592 0.03463 q3H 1.16 1.12493 1.08375 1.04449 1.00091 q4H 1.845731.83626 1.79059 1.77249 1.5543 0.75 q1H 0.1797 0.18134 0.1736 0.16485 -0.19061 q2H 0.694290.4966 0.30539 0.13928 0.02931 q3H 1.102091.0776 1.05131 1.02491 0.99304 q4H 1.705691.70689 1.6986 1.73495 1.665 1 q1H 0.165540.16431 0.16126 0.15787 -0.74749 q2H 0.610760.43819 0.26749 0.1218 0.02059 q3H 1.074911.05268 1.03472 1.01397 0.99449 q4H 1.641141.65266 1.68927 1.77491 1.54847 2 q1H 0.152210.16056 0.17045 0.19478 -0.03603 q2H 0.462490.33927 0.2124 0.09934 0.02634 q3H 1.028571.0214 1.00862 1.00104 0.99125 q4H 1.695181.7541 1.83902 1.98255 1.83895 3 q1H 0.161550.17036 0.17857 0.21098 0.01983 q2H 0.430190.31744 0.20265 0.094 0.02545 q3H 1.020191.01411 1.00731 1.00177 0.9993 q4H 1.784981.83897 1.91251 2.05990 1.91793

The final point regarding simultaneous fracturing between parallel wells is that fractures may propagate in between the wells with some offset between their tips.Consequently,the effect of tip offset between the meeting fractures should also be studied.Fig. 8 depicts the results for aspect ratios of 1 and 0.4.

不知不觉,已是腊八了,中午戴主任请办公室几个人吃饭,易非脚疼,走在后面,戴主任看见了,拍了一下她的肩膀,说:“易非,怎么了?路滑摔了?”

Fig. 8 shows that the effect of offset is similar to the effect of tip distance.This means that SIFHraise is reduced by increasing the tip distance and offset.For the case of offset,however,the decrease of the SIFHis smoother by increasing the offset.This means that increasing the distance between the tips removes the effect of simultaneous fracturing in a sharper fashion than the increase of offset.Comparing the results shown in Fig. 8,it is evident that aspect ratio of the fractures also plays the same role,i.e.the lower the aspect ratio of the fracture,the lower the raise in SIFHof the meeting tips.

6.Effect of simultaneous fracturing on SIF of multistage fractures in parallel wells

In simultaneous multistage hydraulic fracturing,effect of meeting fracture tips is mixed with the effect of multistage fracturing.In this section,we intend to investigate the coupled effect of simultaneous multistage hydraulic fracturing by comparing the SIF of this scenario(SIF(H,L)SM)with thatof a singlestagefracture froma standalonewell(SIF(H,L)O).It is clear from the previous sections that multistage fracturing has negative influence on both SIFLand SIFH(i.e.decreasing effect),while simultaneous fracturing has positive influence on SIFH(i.e.increasing effect).Therefore,the two effects are coupled according to the fracturing configuration whether positive(increasing SIF)or negative(decreasing SIF)effect is observed for the height of the fracture.The coupled effect of simultaneous multistage fracturing on the SIFHand SIFLis shown in Fig. 9.It was observed that SIF shows the same behavior for all aspect ratios,typically for two aspect ratios of 1 and 0.4.The depicted surfaces shown in Fig. 9 are for SIFH(Fig. 9a and b)and for SIFL(Fig. 9c and d)of the fractures.It is observed that SIFHincreaserelates to the smaller distance between the tips and the larger spacing between multistage fractures.There are also Lsand Lp bounding values,beyond which there is no change of SIF during fracturing.

Table 5 Coefficients of the predicting function for SIFLchange as a result of simultaneous multistage fracturing.

Tip distance,Ls/c CoefficientsAR=1 AR=0.8 AR=0.6 AR=0.4 AR=0.2 0.25 q1L 0.1003 0.13374 0.15955 0.20065 0.19623 q2L 0.5016 0.42357 0.31912 0.18309 0.06155 q3L 1.05607 1.02335 1.00338 0.99307 0.98463 q4L 1.35512 1.52942 1.71935 1.98448 2.11508 0.5 q1L 0.11448 0.14275 0.16577 0.20523 0.2136 q2L 0.45853 0.39382 0.29901 0.17514 0.06038 q3L 1.04595 1.02318 1.01119 1.00808 0.99511 q4L 1.4979 1.65748 1.8443 2.09448 2.22139 0.75 q1L 0.12955 0.1548 0.17444 0.2103 0.22215 q2L 0.44295 0.38171 0.29288 0.172 0.06069 q3L 1.03794 1.01984 1.01349 1.00691 0.999 q4L 1.59708 1.74768 1.92115 2.15319 2.25454 1 q1L 0.142 0.16354 0.18298 0.21455 0.22451 q2L 0.43476 0.37644 0.28962 0.16983 0.06114 q3L 1.03249 1.01855 1.01175 1.0062 1.00006 q4L 1.67282 1.80997 1.9678 2.18542 2.25689 2 q1L 0.16136 0.17859 0.19345 0.22131 0.23672 q2L 0.42004 0.36649 0.28339 0.16834 0.05922 q3L 1.02062 1.01412 1.00835 1.00857 0.99849 q4L 1.796 1.91828 2.04427 2.23341 2.30633 3 q1L 0.16544 0.18171 0.19298 0.22098 0.23209 q2L 0.418 0.36501 0.28388 0.16886 0.06044 q3L 1.01861 1.0153 1.01114 1.00562 0.99843 q4L 1.8275 1.94703 2.06327 2.22831 2.29117

Another point which is seen from Fig. 9a and b is that when spacing between multistage fractures is too small,simultaneous fracturing has no influence to keep the SIFHin an acceptable range.This means that multistage fracturing effect dominates the SIF change for too close spacing and causes its substantial reduction,while simultaneous fracturing,even with really low tip distances,cannot remedy this reduction.

Comparing the results shown in Fig. 9c and d with those in Fig. 9a and b,it is revealed that simultaneous fracturing technique has no significant effect on the SIFLof the hydraulic fracture.This is benef i cial for hydraulic fracturing optimization,because by an efficient selection of spacing between the fractures,SIFHis raised as a result of simultaneous fracturing,while SIFLis reduced by multistage fracturing.This causes the fracture to remain contained in the net play and does not penetrate in the bounding layers.

Similar to the previous sections,it is also important to quantify the SIF change according to the distance between fracture tips and spacing between adjacent fractures in simultaneous multistage fracturing technique.Based on the behaviors of the numerical results,SIF changes can be predicted using the following equation:

whereis the simultaneous multistage SIF ratio which determines the SIF change with respect to the case of single-stage standalone fracturing.This ratio is estimated using four coefficients qifor i=1 to 4 that are reported in Tables 4 and 5,respectively,for SIFHand SIFL.One important point in using simultaneous multistage fracturing is that shorter fracture spacing can be selected,because simultaneous fracturing remedies the negative effect of spacing decrease.

7.Conclusions

A comprehensive simulation framework was proposed in order to investigate the interaction between fractures and the influence of closely spaced fractures on SIF.Different scenarios were considered,i.e.specifically fracture aspect ratio,spacing between multistage fractures,and distance between simultaneous fractures.Firstly,it was shown that multistage fracturing can reduce the SIF and aperture of the propagating fractures.The level of this decrease is more severe for higher aspect ratios and shorter spacing between fractures.Secondly,it was observed that the SIF of the meeting hydraulic fractures increases noticeably as a result of simultaneous fracturing.The magnitude of change is again more severe for higher aspect ratios and shorter distances between the tips.It was also observed that the aperture of fractures is not influenced by simultaneous fracturing scheme.Thirdly,it was observed that by having simultaneous multistage fracturing of parallel wells,the increasing effect of simultaneous fracturing is coupled with decreasing effect of multistage fracturing.Decreasing of the SIF and lower propagation potential are seen for propagation in vertical direction,while depending on the spacing(between multistage fractures)and distance(between simultaneous fractures),the resultant effect can be increasing or decreasing for propagation in horizontal direction.Fourthly,it was observed that offset reduces the intensifying effect of simultaneous hydraulic fracturing.All the above-mentioned decreasing and increasing effects on SIF and aperture in different fracturing scenarios were quantif i ed using appropriate fitting equations and the final equations together with their associated coefficients were reported.Using the proposed equations,engineers can predict the SIF change as a result of the designed fracturing pattern.It enables improved planning and placement of productive hydraulic fracture treatments;it offers the potential for considerable cost reductions in completion design and implementation;and it allows for an optimal simultaneous multistage hydraulic fracture treatment that drains larger volumes of the reservoir.

Conflict of interest

The authors wish to confirm that there are no known Conflicts of interest associated with this publication and there has been no significant financial support for this work that could have influenced its outcome.

Acknowledgement

The authors would like to thank Mr.Timothy Beard(Manager-ETG Operations,Chesapeake Energy Corporations),and Prof.Arul Britto(Emeritus faculty,University of Cambridge)for their productive advice.Final thanks are also given to Oklahoma Department of Transportation(ODOT)and Oklahoma Transportation Center for their financial support during the course of this study.

References

Andrews A,Folger P,Humphries M,Copeland C,Tiemann M,Meltz R,Brougher C.Unconventional gas shales:development,technology,and policy issues.Congressional Research Service Report for Congress,2009.Available at http://www.fas.org/sgp/crs/misc/R40894.pdf.

Bao J,Liu H,Zhang G,Jin J,Cheng W,Liu J.Fracture propagation laws in staged hydraulic fracturing and their effects on fracture conductivities.Petroleum Exploration and Development 2017;44(2):306-14.

Broek D.Elementary Engineering fracture mechanics.3rd ed.Boston:Martinus Nijhoff Publishers;1982.

Detournay E.Propagation regimes of fluid-driven fractures in impermeable rocks.International Journal of Geomechanics 2004;4(1):35-45.

Guidera JT,Lardner RW.Penny-shaped cracks.Journal of Elasticity 1975;5:59-73.

Irwin G.Analysis of stresses and strains near the end of a crack traversing a plate.Journal of Applied Mechanics 1957;24:361-4.

Kassir MK,Sih GC.Three-dimensional stress distribution around an elliptical crack under arbitrary loadings.Journal of Applied Mechanics 1966;33(3):601-11.

Kassir MK.Stress-intensity factor for a three-dimensional rectangular crack.Journal of Applied Mechanics 1981;48(2):309-12.

Keer LM.A class of non-symmetrical punch and crack problems.Quarterly Journal of Mechanics and Applied Mathematics 1964;17(4):423-36.

King GE,Haile L,Shuss J,Dobkins TA.Increasing fracture path complexity and controlling downward fracture growth in the Barnett shale.In:Proceedings of the society of petroleum engineers(SPE)shale gas production conference.SPE;2008.

Mastrojannis EN,Keer LM,Mura T.Stress intensity factor for a plane crack under normal pressure.International Journal of Fracture 1979;15(3).1979.

Morrill JC,Miskimins JL.Optimizing hydraulic fracture spacing in unconventional shales.In:SPE Hydraulic Fracturing Technology Conference.SPE;2012.

Murakami Y.Stress intensity factors handbook,vol.1.Pergamon;1987.

Nejati M,Paluszny A,Zimmerman RW.A disk-shaped domain integral method for the computation of stress intensity factors using tetrahedral meshes.International Journal of Solids and Structures 2015;69-70:230-51.

Roussel NP,Sharma MM.Strategies to minimize frac spacing and stimulate natural fractures in horizontal completions.In:SPE Annual Technical Conference and Exhibition.SPE;2011.

Shah RC,Kobayashi AS.Stress intensity factor for an elliptical crack under arbitrary normal loading.Engineering Fracture Mechanics 1971;3(1):71-96.

Sih GC.Handbook of stress-intensity factors:stress-intensity factor solutions and formulas for references.Institute of Fracture and Solid Mechanics,Lehigh University;1973.

Singh I,Miskimins JL.A numerical study of the effects of packer-induced stresses and stress shadowing on fracture initiation and stimulation of horizontal wells.In:Canadian unconventional resources and international petroleum conference;2010.

Sneddon IN.The Distribution of stress in the neighbourhood of a crack in an elastic solid.Proceedings of the Royal Society of London.Series A,Mathematical and Physical Sciences 1946;187(1009):229-60.

Sneddon IN,Elliot HA.The opening of a Griffith crack under internal pressure.Quarterly of Applied Mathematics 1946;4(3):262-7.

Sneddon IN,Lowengrub M.Crack problems in the classical theory of elasticity.New York:John Wiley and Sons;1969.

Soliman MY,Boonen P.Review of fractured horizontal wells technology.In:Proceedings of the Abu Dhabi international petroleum exhibition and conference.SPE;1997.

Stone TJ,Babuska I.A numerical method with a posteriori error estimation for determining the path taken by a propagating crack.Computer Methods in Applied Mechanics and Engineering 1998;160(3-4):245-71.

Tada H,Paris PC,Irwin GR.The stress analysis of cracks handbook.3rd ed.New York:ASME Press;2000.

Taghichian A,Zaman M,Devegowda D.Stress shadow size and aperture of hydraulic fractures in unconventional shales.Journal of Petroleum Science and Engineering 2014;124:209-21. http://www.sciencedirect.com/science/journal/09204105/124/supp/C.

US Energy Information Administration(EIA).Annual energy outlook.2011.http://www.eia.gov.

Ali Taghichian,Hamid Hashemalhoseini,Musharraf Zaman,Saied Beheshti Zavareh
《Journal of Rock Mechanics and Geotechnical Engineering》2018年第2期文献

服务严谨可靠 7×14小时在线支持 支持宝特邀商家 不满意退款

本站非杂志社官网,上千家国家级期刊、省级期刊、北大核心、南大核心、专业的职称论文发表网站。
职称论文发表、杂志论文发表、期刊征稿、期刊投稿,论文发表指导正规机构。是您首选最可靠,最快速的期刊论文发表网站。
免责声明:本网站部分资源、信息来源于网络,完全免费共享,仅供学习和研究使用,版权和著作权归原作者所有
如有不愿意被转载的情况,请通知我们删除已转载的信息 粤ICP备2023046998号