GB2624744A

GB2624744A - Method for optimizing fracturing cluster spacing near deep shale gas fault

Info

Publication number: GB2624744A
Application number: GB2313497.6A
Authority: GB
Inventors: ren Lan; Yang Ziqing; Zhao Jinzhou; Lin Ran; Yu Zhihao; Wu Jianfa
Original assignee: Southwest Petroleum University
Current assignee: Southwest Petroleum University
Priority date: 2022-11-23
Filing date: 2023-09-05
Publication date: 2024-05-29
Also published as: GB202313497D0; CN116306330A

Abstract

Data is acquired. Next a characterization model of an inhomogeneous in-situ stress field is established, calculate a stress field near a pre-fractured fault and produce stress distribution maps (fig.2-3). A hydraulic fracture network is established, hydraulic fracture extension, in-fracture pressure, filtration, and fracture location is calculated. A current reservoir pressure, the in-fracture pressure and the fracture position are solved using a finite difference method for the near fault pressure. A Warpinski two-dimensional criterion establishes a failure criterion for determining natural fractures and grid element positions for tensile and shear failures. Shear failure volume and tensile failure volume are summed using spatial numerical integration to obtain stimulated reservoir volume (SRV). A hydraulic fracture extension map and SRV of clusters is obtained at a set time step. Optimal cluster spacing is determined from dynamic expansion of the SRV and a graph of SRV with cluster spacing is produced (fig.8).

Description

METHOD FOR OPTIMIZING FRACTURING CLUSTER

SPACING NEAR DEEP SHALE GAS FAULT

TECHNICAL FIELD

[0001] The present invention relates to a method for optimizing fracturing cluster spacing near a deep shale gas fault, and falls within the technical field of shale gas development.

BACKGROUND

[0002] Deep shale gas is the main battlefield of building and scale production, and is an important direction of shale gas exploration and development at present. However, the shale experienced multi-stage tectonic movement, the uplift and denudation effect was obvious, the stratum strength anisotropy was strong, and the faults in the shale stratum were generally developed, which results in the inhomogeneous distribution of the in-situ stress field around the shale, the non-planar, asymmetric and irregular extension behavior of the hydraulic fracture, and then affects the SRV (stimulated reservoir volume). Therefore, under the inhomogeneous in-situ stress field, the appropriate cluster spacing design has important significance to enhance the distribution area of the deep shale pressure fracture network and improve the effect of the deep shale gas after reservoir pressure.

[0003] Some scholars at home and abroad have carried out a series of researches on the characterization of in-situ stress near faults and the design of horizontal well cluster spacing of the shale gas. In terms of in-situ stress characterization, Lin et al, found that the changes of in-situ stress and pore pressure would be affected by fault slip, and provided a pore elasticity calculation method, which could quantitatively predict the 3D in-situ stress and pore pressure of shale reservoirs with fault structure. NI, Wei et al, divided the faults into parallel faults and cross faults, and used the finite element method to simulate the in-situ stress under the two faults. In the cluster spacing design of shale gas horizontal well, Guo et al. found that when the cluster fracture extends, the middle fracture was squeezed and the opening was limited and the extension was insufficient, resulting in failure of sand adding to avoid such situation, the cluster spacing was qualitatively optimized. Liu et al. considered the change of the in-situ stress caused by multiple fractures, and formed the in-situ stress diverting area, which was regarded as the potential reservoir reconstruction area and then quantitatively optimized the cluster spacing.

[0004] In conclusion, the in-situ stress characterization, and the cluster spacing design of horizontal wells with shale gas near faults are more sufficient, but none of the above studies has further analyzed the influence of inhomogeneous in-situ stress field generated by faults on the propagation of pressure fracture network. None of the above cluster spacing designs has targeted reservoir SRV, a key parameter affecting shale fractured yield.

[0005] Therefore, it is urgently necessary to establish a method for optimizing fracturing cluster spacing near a deep shale gas fault, which will help to better deploy the horizontal well position, optimize the interval of perforation clusters, and improve the fracturing effect of shale reservoir near faults.

SUM MARY

[0006] In order to overcome the problems of the prior art, the present invention provides a method for optimizing fracturing cluster spacing near a deep shale gas fault.

[0007] The technical solution provided by the present invention for solving the above-mentioned technical problem is: a method for optimizing fracturing cluster spacing near a deep shale gas fault, which includes the steps of: [0008] Step 1: acquiring fault parameters, reservoir geological parameters, fracturing construction parameters, natural fracture parameters and model grid parameters; [0009] Step 2: establishing a characterization model of an inhomogeneous in-situ stress field near a deep shale gas fault, calculating an induced in-situ stress field generated by a hydraulic fracture by a displacement discontinuity method (DDM) based on an elastic mechanics theory model, calculating the induced in-situ stress by linear superposition of the original in-situ stress field and the fault induced in-situ stress field, then calculating the inhomogeneous in-situ stress field near the pre-fractured fault by a superposition principle, and drawing the in-situ stress deflection angle distribution map and the horizontal in-situ stress difference distribution map near the fault respectively; [0010] Step 3: establishing a hydraulic fracture network expansion model of horizontal well subsection multi-cluster fracturing near a deep shale gas fault, and calculating hydraulic fracture extension parameters, in-fracture pressure, filtration, and fracture location parameters of each cluster combining initial and boundary conditions; [0011] Step 4: taking the current reservoir pressure field distribution as an initial condition, taking the in-fracture pressure calculated in the hydraulic fracture extension part and the fracture position parameters as internal boundary conditions, and solving the reservoir pressure field near the fault using a finite difference method; [0012] Step 5: according to a Warpinski two-dimensional criterion, establishing a failure criterion for determining a natural fracture of any occurrence using a tensor calculation method, and respectively determining grid elements of natural fracture positions of tensile failure and shear failure; [0013] Step 6: calculating the sum of the reservoir shear failure reconstruction volume and the tensile failure reconstruction volume through spatial numerical integration according to the grid element of the natural fracture position of the tensile failure and the shear failure to obtain an SRV (stimulated reservoir volume); [0014] Step 7: repeating steps 4 to 6 until a set time step is reached, and obtaining a hydraulic fracture extension map and an SRV of each cluster at a perforation by numerical simulation; and [0015] Step 8: determining the optimal perforation cluster spacing by simulating the dynamic expansion of the SRV near the inhomogeneous in-situ stress field with an objective of maximizing the SRV, and drawing a graph of the variation of the SRV with the cluster spacing. [0016] In a further technical solution, the fault parameters include a fault type, a fault length, a fault height, a fault tilt angle, a fault strike, and a fault displacement; [0017] the reservoir geological parameters include a stratum maximum horizontal principal in-situ stress, a stratum minimum horizontal principal in-situ stress, a stratum horizontal in-situ stress difference, a stratum rock fracturing toughness, a stratum rock Young's modulus, and a stratum rock Poisson' s ratio; [0018] the fracturing construction parameters include a fracture displacement, a fracturing fluid amount, a fracturing time, a fracturing fluid viscosity, a fracturing leak off factor, a fracturing fluid density, an average concentration in fracturing injection propping agent, the number of perforation clusters, a single cluster perforation number, a perforation hole diameter and a cluster spacing; [0019] the natural fracture parameters include a natural fracture average tilt angle, a natural fracture average approximation angle, a natural fracture average length, and a natural fracture average height; [0020] the model grid parameters include model x-direction boundaries, y-direction boundaries, and z-direction boundaries.

[0021] In a further technical solution, the characterization model of the inhomogeneous in-situ stress field near a deep shale gas fault includes an equation for calculating the induced in-situ stress component and an equation for calculating the inhomogeneous in-situ stress field near the pre-fractured fault; [0022] wherein the equation for calculating the induced in-situ stress component is: Au GD" [2eg1:3 + (e2 -g2)1,11+ ;(g1:, + e1:6)] n 27(1-v) ± CD' e2 F3 -2egE4 + (eF5 -gE6)] 27(1-v) -Au = GD" [2eg173 + (e2 -g2)114 -;(gF5 +e)] 27r (1-v) GT), [2g21-13 + 2eg1-14 + C(e1-15 -gE")] [0023] [0024] 27(1-v) [0025] Aa-= 27(1-v) GDn;(gFE, ep,)+ 27rGOD-, v)EF4 (gr, er)] [0026] A azz = v(Ac Ac) [0027] wherein Ao-," Au).,, Ao-i-and Atrx,), respectively represent the fault-induced in-situ stress component, INSPa, G represents the shear modulus of the reservoir rock, Pa; represents Poisson's ratio of reservoir rocks, dimensionless; Ai and A represent a normal fault distance and tangential fault displacement of fault respectively, m; 4-represents a value of y in a global coordinate system converted into a local coordinate system; e and g are the cosine values of an included angle between a E axis of the local coordinate system and an x-axis and y-axis of the global coordinate system, respectively; and FA represents a partial derivative equation of Papkovitch function, k E {3-6); [0028] the equation for calculating the in-situ stress field near the pre-fractured fault is (0) (0 (0) a-r; + Act, o-) a-, + A a-" (0)(0) (0) Cryv aA ± a-yy xv Cr vz (0) (0) rr(o) + Aa + A rrn o-, -zz. --zz [0029] axv Crxv Cr xz a ay)) Cr vz

C zz

[0030] wherein ai.i(o), 0"i) azz(o), (zyz(°), and ar" represent the original in-situ stress value component, Pa; (T, ari, (Tr, (7,,, mr, and (T, represent the current in-situ stress value component, Pa.

[0031] In a further technical solution, the fracture network extension model of the horizontal well subsection multi-cluster fracturing hydraulic fracture near the deep shale gas fault includes a material balance equation, an in-fracture flow equation a fracture opening equation, a fracture height equation, a fracture diverting equation, a fracture flow distribution equation, a fracture extension boundary condition, and an initial condition equation; [0032] wherein the material balance equation is: [0033] aqt-(s, t) 2qhf (S) &Of (S, hi (s) aS \it -r(s) at [0034] qf =1 i=1 [0035] wherein qf represents an in-fracture flow, m3/s; t represents time, s; Iv represents a height of the fracture, m; col-represents an opening of the fracture, m; CI represents a leak off factor of an injected fracturing fluid, m/s0.5; i(s) represents the time at which the fracturing fluid begins to drain at a position along a length direction of the fracture, s; /V represents the number of hydraulic fracture; qr represents a total flow rate of fracturing fluid injection, m3/s; and qi represents the flow rate obtained from the distribution of the iLI1 fracture, m3/s; [0036] wherein the in-fracture flow equation is: R-h f a apf)= qL(s,t)lif + [0037] (0) (s t) 64# as as at hf [0038] wherein pf represents a fluid pressure within the fracture, Pa; ft represents a viscosity of fluid in the fracture, Pa. s; and represents a leak off rate, m/s; [0039] wherein the fracture opening equation is: [0040] u.-( t). +E( [0041] A A nn)** (Un) * 11) = ±( ).(ut).+E( Y

Y J

7=1 J=1 [0042] W1 (s) n (in)7 se/ [0043] wherein A represents a total number of discrete elements of the fracture; i and represent fracture units, with a value of 1-A; (14),J represents a tangential in-situ stress component on an i-unit caused by the j-unit tangential displacement discontinuity, (Mm), represents a tangential in-situ stress component on the i-unit caused by the j-unit normal displacement discontinuity, (Min),, represents a normal in-situ stress component on the i-unit caused by the j-unit tangential displacement discontinuity, and (M"")" represents a normal in-situ stress component on the i-unit caused by the j-unit normal displacement discontinuity; (GO, and (04-respectively represent the shear in-situ stress and normal in-situ stress of the unit i in a local coordinate system, Pa; (LA), and (G), respectively represent tangential strain and normal strain in the local coordinate system, m; (G), i represents a normal offset of the unit i, m; and s represents a fracture length direction coordinate, m, [0044] wherein the fracture height equation is:

K

[0045] -2 ( )2 it Pf ciclose [0046] wherein crelose represents closing in-situ stress acting on a fracture wall, Pa Kic represents fracturing toughness of a reservoir rock, Pa/m°5; [0047] wherein the fracture diverting equation is: [0048] Ki sin + (3 cos RIF 0 = 127EG [0049]- , 4(1-v)V a n

VTTEG

[005o] K1 - [0051] wherein a represents a half-length of a discrete fracture unit, m; thir represents a hydraulic fracture extension diverting angle, 0; [0052] wherein the fracture flow distribution equation is: [0053] ph = ph,/ + ps,1 + AT)" Ap fp j 1=1 [0054] wherein / represents each cluster number; k represents a horizontal section number; ph represents a heel end pressure of a horizontal well, Pa; ph./ represents the fracturing pressure at the l-th cluster fracture hole, Pa; Psi represents the net pressure at a fracture of the first cluster of fractures, Pa; Apr,i represents a flow pressure drop of the klh horizontal section, Pa; Aprp.i represents a frictional pressure drop at the Ph cluster fracture hole; [0055] wherein the equation for the fracture extension boundary conditions and the initial conditions is: 4(1 -v),la 4 1 t _ 0 = 0 Wf 18-4 =0

_

Pf s=Lf _ -" dose qf It 0 = qt [0056] s=14 [0057] wherein Lc represents a hydraulic fracture length, m.

[0058] In a further technical solution, the solving equation in step four includes a continuity equation and an equation for the equivalent permeability after natural fracture failure; [0059] wherein the continuity equation is: a ( , ap) a, [0060] Kt y,t)-+-k (x,y,t)-)= ax ax ay ay at [0061] wherein yo represents porosity, dimensionless; p represents fluid pressure, IV1Pa; CL represents reservoir comprehensive compressibility, MPa1; and ki. represent permeability in the x and y directions, respectively, D; [0062] wherein the equation for the equivalent permeability after natural fracture failure is:

CO

[0063] k", = 12L [0064] kv(x, y, t)= sin 2 0 koF + kc [0065] icy (X, y, t)= COS2 0-kDF+ kc [0066] wherein Ltf represents a spacing of natural fractures, m; Off represents opening of the natural fracture, m; and Ice represents initial permeability, D. [0067] In a further technical solution, the failure criteria for determining any occurrence of a natural fracture includes a unit normal vector of the natural fracture, a normal in-situ stress on a surface of the natural fracture, a shear in-situ stress on a surface of the natural fracture, a determination equation for a tensile failure of the natural fracture and a determination equation for a shear failure of the natural fracture; [0068] wherein the unit normal vector of the natural fracture is: [0069] n =InJ e.1= (n n J x [0070] wherein the normal in-situ stress at a natural fracture face is: [0071] Crin=F*n=n.o-.n jk k [0072] wherein the shear in-situ stress at the natural fracture face is: [0073] = -crnn crnn nn nn jk- jk [0074] wherein the determination equation of the natural fracture tensile failure is: [0075] p > o-" + Ts [0076] wherein Ts represents the natural fracture tensile strength, Pa; and pf represents a fluid pressure within the fracture, IVIPa; [0077] wherein the determination equation of the natural fracture shear failure is: [0078] C > Fni + f. (a,111 -p [0079] wherein.fims represents a friction factor of natural fractures, dimensionless; represents cohesion of a natural fracture, MPa.

[0080] In a further technical solution, the equation for calculating the SRV in Step 6 is.

[0081] V = V ± V = AX(E) * Ay(s)* Az(a) twat s t SEE,1JEES, [0082] wherein represents a total reservoir reconstruction volume, m3, Vs represents a reservoir shear failure reconstruction volume, m3; Ij represents a reservoir tensile failure reconstruction volume, m3; es represents a grid element of the shear failure, ei represents the grid element of the tensile failure; Ax(e), Ay(s) an Az(e) represent a side length of the grid element in the x, y, and z directions, m.

[0083] In a further aspect, the hydraulic fracture extension parameters include a fracture extension length, a fracture extension height, a fracture opening, in-fracture pressure and an in-fracture flow.

[0084] In a further technical solution, the hydraulic fracture extension parameter obtained in step 7 converges with the hydraulic fracture extension parameter obtained in the previous step; wherein the fracture opening is a fracture width, and if the calculated fracture width does not converge, then the pressure in the fracture is changed to perform iterative calculation until the fracture width converges, if the calculated in-fracture flow rate does not converge, then the fracture length is changed to carry out iterative calculation until the in-fracture flow rate converges; if the calculated leak off does not converge, the in-fracture pressure is changed for iterative calculation until the leak off converges [0085] The present invention has the following advantageous effects: [0086] 1, the method is targeted at the mechanical characteristics of shale reservoir near faults, a characterization model of inhomogeneous in-situ stress field near a deep shale gas fault, a hydraulic fracture network expansion model of multi-cluster fracturing in horizontal wells near a deep shale gas fault, and a numerical calculation model of SRV near a deep shale gas fault are established, and thus a method for optimizing fracturing cluster spacing near a deep shale gas fault is provided; [0087] 2. due to the inhomogeneous distribution of in-situ stress near deep shale reservoir faults, the hydraulic fractures may have non-planar, asymmetric and irregular extension behavior, and the influence behavior of cluster spacing on the SRV of shale gas reservoir pressure in horizontal well sub-cluster fracturing is comprehensively considered. An ultimate goal of cluster spacing design is to maximize fractured yield Therefore, according to the inhomogeneous distribution of in-situ stress near deep shale reservoir faults, the hydraulic fracture network expansion model can be established, the SRV can be calculated numerically, the SRV is taken as the optimization objective, the cluster spacing can be optimized to solve the problem of optimizing the horizontal well fracturing cluster spacing near the lack of a deep shale gas fault.

BRIEF DESCRIPTION OF DRAWINGS

[0088] FIG. IA and FIG 1B is a computational flow diagram according to the present invention; [0089] FIG 2 is a graph of the in-situ stress deflection angle distribution near a deep shale gas faults; [0090] FIG 3 is a graph of horizontal in-situ stress difference distribution near a deep shale gas fault; [0091] FIG. 4 is a graph of the inhomogeneous in-situ stress field distribution near a deep shale gas fault; [0092] FIG. 5 is a graph of the fracture network spread at a perforation cluster spacing of 5 m; [0093] FIG 6 is a graph of the fracture network spread at a perforation cluster spacing of 7m; [0094] FIG. 7 is a graph of the fracture network spread at a perforation cluster spacing of 8m; [0095] FIG. 8 is a graph of SRV (stimulated reservoir volume) as a function of cluster spacing.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0096] The technical solutions of the present invention will be described clearly and completely with reference to the drawings, and it should be apparent that the described embodiments are some, but not all, embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by a person of ordinary skill in the art without inventive effort fall within the scope of the present invention. [0097] A method for optimizing fracturing cluster spacing near a deep shale gas fault according to the present invention includes the steps of: [0098] Step 1: acquiring fault parameters, reservoir geological parameters, fracturing construction parameters, natural fracture parameters and model grid parameters; [0099] wherein the fault parameters include a fault type, a fault length, a fault height, a fault tilt angle, a fault strike, and a fault displacement; [00100] the reservoir geological parameters include a stratum maximum horizontal principal in-situ stress, a stratum minimum horizontal principal in-situ stress, a stratum horizontal in-situ stress difference, a stratum rock fracturing toughness, a stratum rock Young's modulus and a stratum rock Poisson' s ratio; [00101] the fracturing construction parameters include a fracture displacement, a fracturing fluid amount, a fracturing time, a fracturing fluid viscosity, a fracturing leak off factor, a fracturing fluid density, an average concentration in fracturing injection propping agent, the number of perforation clusters, a single cluster perforation number, a perforation hole diameter and a cluster spacing; [00102] the natural fracture parameters include a natural fracture average tilt angle, a natural fracture average approximation angle, a natural fracture average length, and a natural fracture average height, etc.; [00103] the model grid parameters include model x-directi on boundaries, y-directi on boundaries, and z-direction boundaries; [00104] Step 2: establishing a characterization model of an inhomogeneousin-situ stress field near a deep shale gas fault, calculating an induced in-situ stress field generated by a hydraulic fracture by a displacement discontinuity method (DDM) based on an elastic mechanics theory model, calculating the induced in-situ stress by linear superposition of the original in-situ stress field and the fault induced in-situ stress field, then calculating the inhomogeneous in-situ stress field near the pre-fractured fault by a superposition principle, and drawing the in-situ stress deflection angle distribution map and the horizontal in-situ stress difference distribution map near the fault respectively; [00105] the characterization model of the inhomogeneous in-situ stress field near a deep shale gas fault includes an equation for calculating the induced in-situ stress component and an equation for calculating the inhomogeneous in-situ stress field near the pre-fractured fault; [00106] wherein the equation for calculating the induced in-situ stress component is:

GD

A o-= [2egF3+(e2. -g2)F4+ C(gFs+cF6)] 27(1-v) GI)/ [2e2E, -2egt:4 -gE6)] 27(1-v) - GDn [2egli;+(e2 -g2)1,:$-;(glis+ el-16)] 27(1-v) GI)/ [2g2E3 + 2egF +C(el-/5 -gF6)] 27(1-v) [00109] Ao-- GD " ch GD;)+ [F, + + eF6)] 27r(1-v) 27-/-(1-v) [00110] Au= = v(Acr" + Acr,}) (4) [00111] wherein Ao, Ao-n, Ao---and Ao-,ky respectively represent the fault-induced in-situ stress component, MPa; G represents the shear modulus of the reservoir rock, Pa; i' represents Poisson's ratio of reservoir rocks, dimensionless; Dr, and TA represent a normal fault distance and tangential fault displacement of fault respectively, m; Crepresents a value of y in a global coordinate system converted into a local coordinate system; e and g are the cosine values of an included angle between a E axis of the local coordinate system and an x-axis and y-axis of the global coordinate system, respectively, and Fic represents a partial derivative equation of Papkovitch function, k e 13-6 [00112] wherein the equation for calculating the pre-fractured inhomogeneous in-situ stress field is: [00107] [00108] (1) (2) (3) [00113] yy (0) (5) cr(0) + Ao- a o-(0) +Acr o_(0)(0) (0) o-+ Ao- 0-lz (0) (0) + Ao- (0) az, rz a-crzz -A 0-zz [00114] wherein cryyon, azzo axv(0) (0), and air(°) represent the original in-situ stress value component, Pa; an, an, azz, Oiv, 05z, and an represent the current in-situ stress value component, Pa; [00115] Step 3: establishing a hydraulic fracture network expansion model of horizontal well subsection multi-cluster fracturing near a deep shale gas fault, and calculating hydraulic fracture extension parameters, in-fracture pressure, filtration and fracture location parameters of each cluster combining initial and boundary conditions; [00116] the fracture network extension model of the horizontal well subsection multi-cluster fracturing hydraulic fracture near the deep shale gas fault includes a material balance equation, an in-fracture flow equation, a fracture opening equation, a fracture height equation, a fracture boundary diverting equation, a fracture flow distribution equation, a fracture extension condition arid an initial condition equation; [00117] wherein the material balance equation is: aqf (s, t) 2Chhf (s) a/Of (s, t) hi (s) [00118] as -r(s) at (6) [00119] qi (7) 1=1 [00120] wherein qt represents an in-fracture flow, m3/s; t represents time, s; /it represents a height of the fracture, m; cur represents an opening of the fracture, m; C1. represents a leak off factor of an injected fracturing fluid, m/s0.5; r(s) represents the time at which the fracturing fluid begins to drain at a position along a length direction of the fracture, s; N represents the number of hydraulic fracture; Tr represents a total flow rate of fracturing fluid injection, m3/s; and qi represents the flow rate obtained from the distribution of the ith fracture, m3/s; [00121] wherein the in-fracture flow equation is: 7-chf a (co3 -ap)= q, (s,t)hf +acor(s hf (8) [00122] 64,ti as f8s at [00123] whereinLk represents a fluid pressure within the fracture, Pa; p represents a viscosity of fluid in the fracture, Pa. s; and represents a leak off rate, m/s; [00124] wherein the fracture opening equation is: [00125] (Un).] (2) [00126] ) (10) [0W27] (CFA = ) (U3i + ( J=1 = ) (Of(s) = (un), j=1 (11) [00128] wherein A represents a total number of discrete elements of the fracture; i and represent fracture units, with a value of 1-A; (Mtt)o represents a tangential in-situ stress component on an i-unit caused by the j-unit tangential displacement discontinuity, (Mot),i represents a tangential in-situ stress component on the i-unit caused by the j-unit normal displacement discontinuity, (Mtn),, represents a normal in-situ stress component on the i-unit caused by the j-unit tangential displacement discontinuity, and (Moo)" represents a normal in-situ stress component on the i-unit caused by the j-unit normal displacement discontinuity; (ch), and (co)Trespectively represent the shear in-situ stress and normal in-situ stress of the unit i in a local coordinatc systcm, Pa; (U), and (Uo), respectively represent tangential strain and normal strain in the local coordinate system, m; (Uo), i represents a normal offset of the unit i, m; and s represents a fracture length direction coordinate, m; [00129] wherein the fracture height equation s: 2 KIC [00130] (12) It P -aclose [00131] wherein o-dose represents closing in-situ stress acting on a fracture wall, Pa; /Cte represents fracturing toughness of a reservoir rock, Paima5; [00132] wherein the fracture diverting equation is: [00E3] K1 sin 0,11, + K. (3 cost91 -1)=0 03)

ITEG

[00134] K1 c tin (14) 4(1-v),/a slTnG [00135] Kll = ,-Ut 05) -±,la [00136] wherein a represents a half-length of a discrete fracture unit, m, Ora represents a hydraulic fracture extension diverting angle, 0; [00137] wherein the fracture flow distribution equation is: [00138] Ph= Phi ± Psj ± APr k APfp,1 (16) )=I [00139] wherein / represents each cluster number; k represents a horizontal section number; ph represents a heel end pressure of a horizontal well, Pa, pb./ represents the fracturing pressure at the 1th cluster fracture hole, Pa; Ps., represents the net pressure at a fracture of the first cluster of fractures, Pa; Apr,k represents a flow pressure drop of the klb horizontal section, Pa; /Imp./ represents a frictional pressure drop at the Ph cluster fracture hole; [00140] wherein the equation for the fracture extension boundary conditions and the initial conditions is: 1-Lt. 11=0 =0 Wf 18=4 =0 Pf 15=4, qf 1 s=Lf C1 = 13 1 I t=0 = qi = close [00141] (17) [00142] wherein Li represents a hydraulic fracture length, m; [00143] Step 4: taking the current reservoir pressure field distribution as an initial condition, taking the in-fracture pressure calculated in the hydraulic fracture extension part and the fracture position parameters as internal boundary conditions, and solving the reservoir pressure field near the fault using a finite difference method; [00144] the solving equation in step four includes a continuity equation and an equivalent permeability equation after natural fracture failure; [00145] wherein the continuity equation is: a op\ a [00146] ,'k(x,y,t)-+-k \ Op (18), , t y \ci [00147] wherein co represents porosity, dimensionless; p represents fluid pressure, MPa; represents reservoir comprehensive compressibility, MiPa-1; Ict and ky, represent permeability in the x and y directions, respectively, D; [00148] wherein the equation for the equivalent permeability after natural fracture failure is: [00149] kir = ' (19) 12LE/, [00150] kx(x, y, = sin2 0 * kDF (20) [00151] Ic(x,y,t)=cos2 8k1. +k (21) [00152] wherein Ltf represents a spacing of natural fractures, m; (off represents opening of the natural fracture, m; and kC represents initial permeability, D [00153] Step 5: according to a Warpinski two-dimensional criterion, establishing a failure criterion for determining a natural fracture of any occurrence using a tensor calculation method, and respectively determining grid elements of natural fracture positions of tensile failure and shear failure; [00154] wherein the failure criteria for determining any occurrence of a natural fracture includes a unit normal vector of the natural fracture, a normal in-situ stress on a surface of the natural fracture, a shear in-situ stress on a surface of the natural fracture, a determination equation for a tensile failure of the natural fracture and a determination equation for a shear failure of the natural fracture; [00155] wherein the unit normal vector of the natural fracture s: [00156] n= fn e.1} (rc z) (22)

J

[00157] wherein the normal in-situ stress at a natural fracture face is: [00b8] Crnn = F * n=nio-Ank (23) [00159] wherein the shear in-situ stress at the natural fracture face s: [00160] Cr" = \IF F -o-"n * o-", = Vo-jknko-jknk -o-"n * o-n, (24) [00161] wherein the determination equation of the natural fracture tensile failure is: [00162] pi > Cr nn ± Ts (25) [00163] wherein I:, represents the natural fracture tensile strength, Pa; and pr represents a fluid pressure within the fracture, MPa, [00164] wherein the determination equation of the natural fracture shear failure is: [00165] Cyr > Fit fmc * (0 -pf) (26) [00166] wherein fine represents a friction factor of natural fractures, dimensionless; FAi represents cohesion of a natural fracture, MPa.

[00167] Step 6: calculating the sum of the reservoir shear failure reconstruction volume and the tensile failure reconstruction volume through spatial numerical integration according to the grid element of the natural fracture position of the tensile failure and the shear failure to obtain an SRV (stimulated reservoir volume); [00168] the equation for calculating the SRV in Step 6 is: [00169] Vora/ = Vs ± Ft = I Ax(e)*Ay(e)* Az(s) (27) CEC,USEC, [00170] wherein r4,,iai represents a total reservoir reconstruction volume, m3; r; represents a reservoir shear failure reconstruction volume, m3; (it represents a reservoir tensile failure reconstruction volume, m'; es represents a grid element of the shear failure, et represents the grid element of the tensile failure; Ax(8), 4(0 an Az(s) represent a side length of the grid element in the x, y, and z directions, m.

[00171] Step 7: repeating steps 4 to 6 until a set time step is reached, and obtaining a hydraulic fracture extension map and an SRV of each cluster at a perforation by numerical simulation; and [00172] the hydraulic fracture extension parameter obtained in step 7 converges with the hydraulic fracture extension parameter obtained in the previous step; wherein the fracture opening is a fracture width, and if the calculated fracture width does not converge, then the pressure in the fracture is changed to perform iterative calculation until the fracture width converges; if the calculated in-fracture flow rate does not converge, then the fracture length is changed to carry out iterative calculation until the in-fracture flow rate converges; if the calculated leak off does not converge, the in-fracture pressure is changed for iterative calculation until the leak off converges; [00173] wherein the difference between the hydraulic fracture extension parameter calculated in this cycle and the hydraulic fracture extension parameter calculated in the previous cycle is not more than 1%, i.e., convergence; [00174] Step 8: determining the optimal perforation cluster spacing by simulating the dynamic expansion of the SRV near the inhomogeneous in-situ stress field with an objective of maximizing the SRV, and drawing a graph of the variation of the SRV with the cluster spacing.

[00175] Examples

[00176] It is known that the actual parameters of a deep shale gas well near a typical fault are as shown in table 1, and an example calculation is carried out according to the flow chart in FIG lA and FIG 1B.

[00177] Table 1 List of actual parameters of deep shale gas well near a fault Parameter Parameter Numerical Unit type value Fault Tangential fault displacement 0.2 m parameter Length 2121.3 m Reservoir Minimum horizontal principal 94.35 MPa geological in-situ stress parameters Stratum maximum horizontal 117.16 MPa principal in-situ stress Stratum rock Young's modulus 47.2 GPa Stratum rock Poisson's ratio 0.24 Dimensionless Fracturing Fracture displacement 16.0 m3/min construction parameters Fracturing fluid amount 1800 m3 Fracturing time 110 min Fracturing fluid viscosity 3 mPa. s Fracturing leak off factor 0.000084 m/(s°.5) Fracturing fluid density 1000 kg/m3 Average concentration in 82.5 kg/m3 fracturing injection propping agent Number of perforation clusters 6 Cluster Single cluster perforation number 8 Perforation hole diameter 9.5 mm Natural Natural fracture tilt angle 60 0 fracture parameter Natural fracture approximation 15 0 angle Natural fracture length 0.9 m Natural fracture height 0.5 m [00178] Firstly, the inhomogeneous in-situ stress field characterization model near a deep shale gas fault is established, and the displacement discontinuity method is used to solve and calculate the stratum induced in-situ stress component caused by fault structure in any point of Example 1 by combining the following equation [00179] the specific steps are: with simultaneous equations (1)-(5), the induced in-situ stress component caused by fault structure at any point in the reservoir is calculated by a displacement discontinuity method; [00180] Secondly, the in-situ stress deflection angle distribution map near the deep shale gas fault (as shown in FIG. 2) and the horizontal in-situ stress difference distribution map near the deep shale gas fault (as shown in FIG. 3) are plotted; [00181] then, a model for the hydraulic fracture network expansion of horizontal well multi-cluster fracturing near the deep shale gas fault is established using a fluid-solid coupling theory, and the hydraulic fracture extension in Example 1 is calculated with the following equation: [00182] the specific steps are: with simultaneous equations (6)-(11), the extension length, height, opening of each hydraulic fracture and the pressure inside the fracture in the fracturing process near deep shale gas horizontal well fault are calculated by using a finite difference method.

[00183] Then, the deep shale gas reservoir pressure field with the fault structure is solved, and the pressure at any point in the reservoir in Example 1 is calculated in combination with the following equation: [00184] the specific steps axe: with simultaneous equations (12)-(21), the pressure field equation of deep shale gas reservoir under fault structure is solved by a finite difference method, and the pressure at any point in the reservoir is calculated.

[00185] Then, according to a Warpinski's two-dimensional criterion, using the tensor calculation method, in combination with the following equation, the failure criterion for determining the natural fracture of any occurrence in Example 1 is derived: [00186] the specific steps are: based on the simultaneous equations (22)-(26), the failure criteria of any occurrence natural fractures in deep shale reservoirs affected by faults are derived by using a tensor calculation method.

[00187] Then, combined with the following equation, the change of the SRV with the cluster spacing in the horizontal well of shale gas in Example 1 was calculated; [00188] The specific steps are: in conjunction with equation (27), the SRV of a horizontal shale gas well at different cluster spacings in Example 1 is calculated using numerical integration; [00189] then, an example calculation is carried out by using the numerical calculation flow chart of the method for the present invention as shown in FIG. IA and FIG. 1B, and according to the calculation result, the distribution map of the inhomogeneous in-situ stress field near a deep shale gas fault (as shown in FIG. 4) and the fracture network distribution under different perforation cluster spacings are respectively drawn (as shown in FIGS. 5-7: 5 m, 7 m, and 9 m); a plot of SRV as a function of cluster spacing (as shown in FIG. 8). When the cluster spacing of a shale gas horizontal well is 7.1m, the maximum SRI! is 53.1x104m3.

[00190] The above description is not intended to limit the present invention in any way. Although the present invention has been disclosed by the above embodiments, it is not intended to limit the present invention. A person skilled in the art will recognize that changes and modifications may be made to the disclosed embodiments without departing from the spirit and scope of the invention. However, any simple modification, equivalent change and modification to the above embodiments according to the technical spirit of the present invention are still within the scope of the technical solution of the present invention.

Claims

What is claimed is: 1. A method for optimizing fracturing cluster spacing near a deep shale gas fault, comprising the steps of: Step 1: acquiring fault parameters, reservoir geological parameters fracturing construction parameters, natural fracture parameters and model grid parameters; Step 2: establishing a characterization model of an inhomogeneous in-situ stress field near a deep shale gas fault, calculating an induced in-situ stress field generated by a hydraulic fracture by a displacement discontinuity method based on an elastic mechanics theory model, calculating the induced in-situ stress by linear superposition of the original in-situ stress field and the fault induced in-situ stress field, then calculating the inhomogeneous in-situ stress field near the pre-fractured fault by a superposition principle, and drawing the in-situ stress deflection angle distribution map and the horizontal in-situ stress difference distribution map near the fault respectively; Step 3: establishing a hydraulic fracture network expansion model of horizontal well subsection multi-cluster fracturing near a deep shale gas fault, and calculating hydraulic fracture extension parameters, in-fracture pressure, filtration and fracture location parameters of each cluster combining initial and boundary conditions; Step 4: taking the current reservoir pressure field distribution as an initial condition, taking the in-fracture pressure calculated in the hydraulic fracture extension part and the fracture position parameters as internal boundary conditions, and solving the reservoir pressure field near the fault using a finite difference method; Step 5: according to a Warpinski two-dimensional criterion, establishing a failure criterion for determining a natural fracture of any occurrence using a tensor calculation method, and respectively determining grid elements of natural fracture positions of tensile failure and shear failure; Step 6: calculating the sum of the reservoir shear failure reconstruction volume and the tensile failure reconstruction volume through spatial numerical integration according to the grid element of the natural fracture position of the tensile failure and the shear failure to obtain a stimulated reservoir volume; Step 7: repeating steps 4 to 6 until a set time step is reached, and obtaining a hydraulic fracture extension map and a stimulated reservoir volume of each cluster at a perforation by numerical simulation; and Step 8: determining the optimal perforation cluster spacing by simulating the dynamic expansion of the stimulated reservoir volume near the inhomogeneous in-situ stress field with an objective of maximizing the stimulated reservoir volume, and drawing a graph of the variation of the stimulated reservoir volume with the cluster spacing 2. The method for optimizing fracturing cluster spacing near a deep shale gas fault according to claim I, wherein the fault parameters comprise a fault type, a fault length, a fault height, a fault tilt angle, a fault strike, and a fault displacement, the reservoir geological parameters comprise a stratum maximum horizontal principal in-situ stress, a stratum minimum horizontal principal in-situ stress, a stratum horizontal in-situ stress difference, a stratum rock fracturing toughness, a stratum rock Young's modulus, and a stratum rock Poisson' s ratio; the fracturing construction parameters comprise a fracture displacement, a fracturing fluid amount, a fracturing time, a fracturing fluid viscosity, a fracturing leak off factor, a fracturing fluid density, an average concentration in fracturing injection propping agent, the number of perforation clusters, a single cluster perforation number, a perforation hole diameter and a cluster spacing, the natural fracture parameters comprise a natural fracture average tilt angle, a natural fracture average approximation angle, a natural fracture average length, and a natural fracture average height; the model grid parameters comprise model x-direction boundaries, y-direction boundaries, and z-directi on boundaries.3. The method for optimizing fracturing cluster spacing near a deep shale gas fault according to claim 1, wherein the characterization model of the inhomogeneous in-situ stress field near a deep shale gas fault comprises an equation for calculating the induced in-situ stress component and an equation for calculating the inhomogeneous in-situ stress field near the pre-fractured fault; wherein the equation for calculating the inducedin-situ stress component s: o-- GD" [2egr, + (e2 -g2)1/4 + + 2,71-(1-OR [2e2F, -2e gF, + "(eF, -gF,)] 2z(1-v) Aa =GI-1)n 2egF, + (e2 -g2)r, -"(gF, + eP6)] 3)/ 241-v) CD' [2g21-1, +2eghl, + (e1-1, -27-1-(1-v) D, [F4 + 4-(gF5+ eFo)] 2zG(1D C 2z(1-t') eF5)± 27-t-(1-v) Acrzz = v(Acrn + ) yy-wherein AOIXT, Aa, Ao-, and Ao)", respectively represent the fault-induced in-situ stress component, MiPa; G represents the shear modulus of the reservoir rock, Pa; I' represents Poisson's ratio of reservoir rocks, dimensionless; Dn and A represent a normal fault distance and tangential fault displacement of fault respectively, m; r epr esents a value of y in a global coordinate system converted into a local coordinate system; e and g are the cosine values of an included angle between a i axis of the local coordinate system and an x-axis and y-axis of the global coordinate system, respectively; arid 1-,k represents a partial derivative equation of Papkovitch function, k E 13-6); the equation for calculating the in-situ stress field near the pre-fractured fault is: + Cr) (0 (0) A Crxz -Cricz xx xx xy A C (0) Crify FI; (0 + a 0n) -;I) CL0 Ak) Aczz wherein o-no), 0-,y0), cry-0), and cf=0) represent the original in-situ stress value component, Pa; o-n, ayy, ass, ayy, ayz, and orr. represent the current in-situ stress value component, Pa.4. The method for optimizing fracturing cluster spacing near a deep shale gas fault according to claim 1, wherein the fracture network extension model of the horizontal well subsection multi-cluster fracturing hydraulic fracture near the deep shale gas fault comprises a material balance equation, an in-fracture flow equation, a fracture opening equation, a fracture height equation, a fracture diverting equation, a fracture flow distribution equation, a fracture extension boundary condition and an initial condition equation; wherein the material balance equation is: aq,(s,t) _ 2chf(s) ±acof(s,t) as Vt -r(s) i=1 wherein qf represents an in-fracture flow, m3/s, i represents time, s; hf represents a height of the fracture, m; (of represents an opening of the fracture, m; CL represents a leak off factor of an injected fracturing fluid, m/s°.5; t(s) represents the time at which the fracturing fluid begins to drain at a position along a length direction of the fracture, s; N represents the number of hydraulic fracture; eqr represents a total flow rate of fracturing fluid injection, m3/s; and qi represents the flow rate obtained from the distribution of the ill' fracture, m3/s; wherein the in-fracture flow equation is: aop ow (s t) (03,. ), qL(s,011 + f 64p as ' as at wherein pc represents a fluid pressure within the fracture, Pa p represents a viscosity of fluid in the fracture, Pa s; and qi, represents a leak off rate, m/s; wherein the fracture opening equation is:A(0-3, = I ( ) (u,), WLfl) ( ), 1=1A((In = ( A/Ifit (Ut) * + ( M) ) j 11 j=1 j=1 w1 (s) = (tin)/ Ise/ wherein A represents a total number of discrete elements of the fracture; i and/ represent fracture units, with a value of 1-A, (Mii),, represents a tangential in-situ stress component on an i-unit caused by the j-unit tangential displacement discontinuity, (MA), represents a tangential in-situ stress component on the i-unit caused by the j-unit normal displacement discontinuity, represents a normal in-situ stress component on the i-unit caused by the j-unit tangential displacement discontinuity, and (114i11i)i, represents a normal in-situ stress component on the i-unit caused by the j-unit normal displacement discontinuity; (cit.), and (c,i), respectively represent the shear in-situ stress and normal in-situ stress of the unit i in a local coordinate system, Pa; WO, and (G), respectively represent tangential strain and normal strain in the local coordinate system, m; (1/n), i represents a normal offset of the unit i, m; and s represents a fracture length direction coordinate, m; wherein the fracture height equation is: 2( )2 11 Pt-aclose wherein o-ciose represents closing in-situ stress acting on a fracture wall, Pa; Kw represents fracturing toughness of a reservoir rock, Pa/tuft': wherein the fracture diverting equation is: K/ sin OHF K (3 cos Opw -1) = 0K-NITIEG t II 4(1 -ONFI " "TEC; LI wherein a represents a half-length of a discrete fracture unit, m; OuF represents a hydraulic fracture extension diverting angle, 0; wherein the fracture flow distribution equation is: Ph = Ph,/ ± Ps,i +IAA, ± 1=1 wherein / represents each cluster number; k represents a horizontal section number; ph represents a heel end pressure of a horizontal well, Pa; ph,/ represents the fracturing pressure at the 1-th cluster fracture hole, Pa, Po represents the net pressure at a fracture of the first cluster of fractures, Pa, Apr.A represents a flow pressure drop of the V' horizontal section, Pa, Aprp.r represents a frictional pressure drop at the Ph cluster fracture hole; wherein the equation for the fracture extension boundary conditions and the initial conditions is.Pf1 IT 11=0 = ° 6.°1-s=Lf = 0 s=Lf = close 1 qf 1 s=Lf = ° qf 1 t=0 = qi wherein Lc represents a hydraulic fracture length, m.5. The method for optimizing fracturing cluster spacing near a deep shale gas fault according to claim 1, wherein the equation for solving the reservoir pressure field in step 4 comprises a continuity equation and an equation for the equivalent permeability after natural fracture failure; wherein the continuity equation is: a * iDp. a' 5pap kjx,y,t)-+-k(xY,t) = POC axt at ax," ay wherein co represents porosity, dimensionless; p represents fluid pressure, MPa; CL represents reservoir comprehensive compressibility, MPa-1; kx and Icy represent permeability in the x and y directions, respectively, D; wherein the equation for the equivalent permeability after natural fracture failure is: k =DT 12L tfkr(x, y, t)= sin' 0.k" + kr y, t)=cos2 0.k" +kc wherein La represents a spacing of natural fractures, m cow represents opening of the natural fracture, m; and kc represents initial permeability, D. 6. The method for optimizing fracturing cluster spacing near a deep shale gas fault according to claim 1, wherein the failure criteria for determining any occurrence of a natural fracture comprises a unit normal vector of the natural fracture, a normal in-situ stress on a surface of the natural fracture, a shear in-situ stress on a surface of the natural fracture, a determination equation for a tensile failure of the natural fracture and a determination equation for a shear failure of the natural fracture; wherein the unit normal vector of the natural fracture is: n=lniej}=(n. n nz) wherein the normal in-situ stress at a natural fracture face is: arm F * n = njo-jknk wherein the shear in-situ stress at the natural fracture face is: Cr rr =VF'*E a, * Cc? = Valk kcr iknk nn * a nn wherein the determination equation of the natural fracture tensile failure is: Pf Cr. + Ts wherein T, represents the natural fracture tensile strength, Pa; and pf represents a fluid pressure within the fracture, MPa; wherein the determination equation of the natural fracture shear failure is: a rz-Ft?' fmc ((inn P wherein j0 represents a friction factor of natural fractures, dimensionless; F, represents cohesion of a natural fracture, MiPa.7. The method for optimizing fracturing cluster spacing near a deep shale gas fault according to claim 6, wherein the equation for calculating the stimulated reservoir volume in Step 6 is: Viola/ = Vs + SeSslJsest wherein Vtotai represents a total reservoir reconstruction volume, m3; V, represents a reservoir shear failure reconstruction volume, m3; represents a reservoir tensile failure reconstruction volume, m3; s represents a grid element of the shear failure, si represents the grid element of the tensile failure, Auc(s), Ay(c) an Az(s) represent a side length of the grid element in the x, y, and z directions, m.8. The method for optimizing fracturing cluster spacing near a deep shale gas fault according to claim 1, wherein the hydraulic fracture extension parameters comprise a fracture extension length, a fracture extension height, a fracture opening, an in-fracture pressure, and an in-fracture flow.9. The method for optimizing fracturing cluster spacing near a deep shale gas fault according to claim 8, wherein the hydraulic fracture extension parameter obtained in step 7 converges with the hydraulic fracture extension parameter obtained in the previous step; wherein the fracture opening is a fracture width, and if the calculated fracture width does not converge, then the pressure in the fracture is changed to perform iterative calculation until the fracture width converges, if the calculated in-fracture flow rate does not converge, then the fracture length is changed to carry out iterative calculation until the in-fracture flow rate converges; if the calculated leak off does not converge, the in-fracture pressure is changed for iterative calculation until the leak off converges.